Description: Models from phase 5 of the Coupled Model Inter-comparison Project (CMIP5) show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here, we analyze the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), the climate model used at the Al-fred Wegener Institute, Helmholtz Centre for Polar and Ma-rine Research, which employs a sea ice–ocean model com-ponent formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of ∼ 1◦ (CMIP5 type) up to locally eddy resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to reductions in deep ocean biases, although these im-provements are not strictly monotonic for the five different ocean grids. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two key re-gions at the surface are responsible for the development of the deep bias in the Atlantic Ocean: the northeastern North Atlantic and the region adjacent to the Strait of Gibraltar. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcrop-ping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used in the final runs for phase 6 of CMIP (CMIP6) and for the related flagship simulations in the High Resolution Model Intercomparison Project (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice–ocean model enables us to per-form the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day.

Description: The Finite Element Sea Ice-Ocean Model (FESOM) is the first global ocean general circulation model based on unstructured-mesh methods that has been developed for the purpose of climate research. The advantage of unstructured-mesh models is their flexible multi-resolution modelling functionality. In this study, an overview of the main features of FESOM will be given; based on sensitivity experiments a number of specific parameter choices will be explained; and directions of future developments will be outlined. It is argued that FESOM is sufficiently mature to explore the benefits of multi-resolution climate modelling and that its applications will provide information useful for the advancement of climate modelling on unstructured meshes.

Description: 〈jats:p〉Abstract. Computation of barotropic and meridional overturning streamfunctions for models formulated on unstructured meshes is commonly preceded by interpolation to a regular mesh. This operation destroys the original conservation, which can be then artificially imposed to make the computation possible. An elementary method is proposed that avoids interpolation and preserves conservation in a strict model sense. The method is described as applied to the discretization of the Finite volumE Sea ice – Ocean Model (FESOM2) on triangular meshes. It, however, is generalizable to colocated vertex-based discretization on triangular meshes and to both triangular and hexagonal C-grid discretizations. 〈/jats:p〉

Description: 〈jats:p〉Abstract. We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 (Finite-volumE Sea ice–Ocean Model) has the multi-resolution functionality typical of unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of the latest developments in the numerical-weather-prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable-resolution (25–125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other ongoing research activities with AWI-CM3. This includes the exploration of high and variable resolution and the development of a full Earth system model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above-CMIP6-average skills (where CMIP6 denotes Coupled Model Intercomparison Project phase 6) in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model. 〈/jats:p〉

Description: Arctic Ocean simulations in 19 global ocean-sea-ice models participating in the Ocean Model Intercomparison Project (OMIP) of the Coupled Model Intercomparison Project Phase 6 (CMIP6) are evaluated in this paper. Our findings show no significant improvements in Arctic Ocean simulations from the previous Coordinated Ocean-ice Reference Experiments phase II (CORE-II) to the current OMIP. Large model biases and inter-model spread exist in the simulated mean state of the halocline and Atlantic Water layer in the OMIP models. Most of the OMIP models suffer from a too thick and deep Atlantic Water layer, a too deep halocline base, and large fresh biases in the halocline. The OMIP models qualitatively agree on the variability and change of the Arctic Ocean freshwater content; sea surface height; stratification; and volume, heat, and freshwater transports through the Arctic Ocean gateways. They can reproduce the changes in the gateway transports observed in the early 21st century, with the exception of the Bering Strait. We also found that the OMIP models employing the NEMO ocean model simulate relatively larger volume and heat transports through the Barents Sea Opening. Overall, the performance of the Arctic Ocean simulations is similar between the CORE2-forced OMIP-1 and JRA55-do-forced OMIP-2 experiments.

Description: 〈jats:p〉Abstract. The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice and plays a key role in climate variability. It is thus critical that numerical models used in climate studies are capable of a good representation of the mixed layer, especially its depth. Here we evaluate the mixed-layer depth (MLD) in six pairs of non-eddying (1∘ grid spacing) and eddy-rich (up to 1/16∘) models from the Ocean Model Intercomparison Project (OMIP), forced by a common atmospheric state. For model evaluation, we use an updated MLD dataset computed from observations using the OMIP protocol (a constant density threshold). In winter, low-resolution models exhibit large biases in the deep-water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode-water formation regions of the Northern Hemisphere. Results in the Southern Ocean are more contrasted, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability in MLD in winter. Contrary to a hypothesis that the deepening of the mixed layer in anticyclones would make the MLD larger globally, eddy-rich models tend to have a shallower mixed layer at most latitudes than coarser models do. In addition, our study highlights the sensitivity of the MLD computation to the choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects. 〈/jats:p〉

Description: Numerical simulations employing prognostic sta- ble water isotopes can not only facilitate our understanding of hydrological processes and climate change but also al- low for a direct comparison between isotope signals obtained from models and various archives. In the current work, we describe the performance and explore the potential of a new version of the Earth system model AWI-ESM (Alfred We- gener Institute Earth System Model), labeled AWI-ESM-2.1- wiso, in which we incorporated three isotope tracers into all relevant components of the water cycle. We present here the results of pre-industrial (PI) and mid-Holocene (MH) simula- tions. The model reproduces the observed PI isotope compo- sitions in both precipitation and seawater well and captures their major differences from the MH conditions. The sim- ulated relationship between the isotope composition in precipitation (d18Op) and surface air temperature is very similar between the PI and MH conditions, and it is largely consis- tent with modern observations despite some regional model biases. The ratio of the MH–PI difference in δ18Op to the MH–PI difference in surface air temperature is comparable to proxy records over Greenland and Antarctica only when summertime air temperature is considered. An amount effect is evident over the North African monsoon domain, where a negative correlation between δ18Op and the amount of pre- cipitation is simulated. As an example of model applications, we studied the onset and withdrawal date of the MH West African summer monsoon (WASM) using daily variables. We find that defining the WASM onset based on precipitation alone may yield erroneous results due to the substantial daily variations in precipitation, which may obscure the dis- tinction between pre-monsoon and monsoon seasons. Com- bining precipitation and isotope indicators, we suggest in this work a novel method for identifying the commencement of the WASM. Moreover, we do not find an obvious difference between the MH and PI periods in terms of the mean onset of the WASM. However, an advancement in the WASM with- drawal is found in the MH compared to the PI period due to an earlier decline in insolation over the northern location of Intertropical Convergence Zone (ITCZ).

Description: 〈jats:p〉Abstract. This study evaluates the impact of increasing resolution on Arctic Ocean simulations using five pairs of matched low- and high-resolution models within the OMIP-2 (Ocean Model Intercomparison Project phase 2) framework. The primary objective is to assess whether a higher resolution can mitigate typical biases in low-resolution models and improve the representation of key climate-relevant variables. We reveal that increasing the horizontal resolution contributes to a reduction in biases in mean temperature and salinity and improves the simulation of the Atlantic water layer and its decadal warming events. A higher resolution also leads to better agreement with observed surface mixed-layer depth, cold halocline base depth and Arctic gateway transports in the Fram and Davis straits. However, the simulation of the mean state and temporal changes in Arctic freshwater content does not show improvement with increased resolution. Not all models achieve improvements for all analyzed ocean variables when spatial resolution is increased so it is crucial to recognize that model numerics and parameterizations also play an important role in faithful simulations. Overall, a higher resolution shows promise in improving the simulation of key Arctic Ocean features and processes, but efforts in model development are required to achieve more accurate representations across all climate-relevant variables. 〈/jats:p〉