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: This paper presents global comparisons of fundamental global climate variables from a suite of four pairs of matched low- and high-resolution ocean and sea ice simulations that are obtained following the OMIP-2 protocol (Griffies et al., 2016) and integrated for one cycle (1958–2018) of the JRA55-do atmospheric state and runoff dataset (Tsujino et al., 2018). Our goal is to assess the robustness of climate-relevant improvements in ocean simulations (mean and variability) associated with moving from coarse (∼ 1∘) to eddy-resolving (∼ 0.1∘) horizontal resolutions. The models are diverse in their numerics and parameterizations, but each low-resolution and high-resolution pair of models is matched so as to isolate, to the extent possible, the effects of horizontal resolution. A variety of observational datasets are used to assess the fidelity of simulated temperature and salinity, sea surface height, kinetic energy, heat and volume transports, and sea ice distribution. This paper provides a crucial benchmark for future studies comparing and improving different schemes in any of the models used in this study or similar ones. The biases in the low-resolution simulations are familiar, and their gross features – position, strength, and variability of western boundary currents, equatorial currents, and the Antarctic Circumpolar Current – are significantly improved in the high-resolution models. However, despite the fact that the high-resolution models “resolve” most of these features, the improvements in temperature and salinity are inconsistent among the different model families, and some regions show increased bias over their low-resolution counterparts. Greatly enhanced horizontal resolution does not deliver unambiguous bias improvement in all regions for all models.

Description: The evaluation and model element description of the second version of the unstructured-mesh Finite-volumE Sea ice-Ocean Model (FESOM2.0) are presented. The new version of the model takes advantage of the finite-volume approach, whereas its predecessor version, FESOM1.4 was based on the finite-element approach. The model sensitivity to arbitrary Lagrangian–Eulerian (ALE) linear and nonlinear free-surface formulation, Gent–McWilliams eddy parameterization, isoneutral Redi diffusion and different vertical mixing schemes is documented. The hydrographic biases, large-scale circulation, numerical performance and scalability of FESOM2.0 are compared with its predecessor, FESOM1.4. FESOM2.0 shows biases with a magnitude comparable to FESOM1.4 and simulates a more realistic Atlantic meridional overturning circulation (AMOC). Compared to its predecessor, FESOM2.0 provides clearly defined fluxes and a 3 times higher throughput in terms of simulated years per day (SYPD). It is thus the first mature global unstructured-mesh ocean model with computational efficiency comparable to state-of-the-art structured-mesh ocean models. Other key elements of the model and new development will be described in follow-up papers.

Description: We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

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: Disintegration of ice shelves in the Amundsen Sea, in front of the West Antarctic Ice Sheet, has the potential to cause sea level rise by inducing an acceleration of ice discharge from upstream grounded ice. Moore et al. (2018) proposed that using a submarine wall to block the penetration of warm water into the subsurface cavities of these ice shelves could reduce this risk. We use a global sea ice–ocean model to show that a wall shielding the Amundsen Sea below 350 m depth successfully suppresses the inflow of warm water and reduces ice shelf melting. However, these warm water masses get redirected towards neighboring ice shelves, which reduces the net effectiveness of the wall. The ice loss is reduced by 10 %, integrated over the entire Antarctic continent.

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).