Description: Highlights: • Phase II of the Coordinated Ocean-ice Reference Experiments (CORE-II) is introduced. • Solutions from CORE-II simulations from eighteen participating models are presented. • Mean states in the North Atlantic with a focus on AMOC are examined. • The North Atlantic solutions differ substantially among the models. • Many factors, including parameterization choices, contribute to these differences. Simulation characteristics from eighteen global ocean–sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort.

Description: Highlights: • Global mean sea level simulated in interannual CORE simulations. • Regional sea level patterns simulated in interannual CORE simulations. • Theoretical foundation for analysis of global mean sea level and regional patterns. Abstract: We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993–2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.

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 processlevel 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 earliest part of the Holocene, from 11.5k to 7k (k = 1000 years before present), is a critical transition period between the relatively cold last deglaciation and the warm middle Holocene. It is marked by more pronounced seasonality and reduced greenhouse gases (GHGs) than the present state, as well as by the presence of the Laurentide Ice Sheet (LIS) and glacial meltwater perturbation. This paper performs experiments under pre-industrial and different early-Holocene regimes with AWI-ESM (Alfred Wegener Institute–Earth System Model), a state-of-the-art climate model with unstructured mesh and varying resolutions, to examine the sensitivity of the simulated Atlantic meridional overturning circulation (AMOC) to early-Holocene insolation, GHGs, topography (including properties of the ice sheet), and glacial meltwater perturbation. In the experiments with early-Holocene Earth orbital parameters and GHGs applied, the AWI-ESM simulation shows a JJA (June–July–August) warming and DJF (December–January–February) cooling over the mid and high latitudes compared with pre-industrial conditions, with amplification over the continents. The presence of the LIS leads to an additional regional cooling over the North America. We also simulate the meltwater event around 8.2k. Big discrepancies are found in the oceanic responses to different locations and magnitudes of freshwater discharge. Our experiments, which compare the effects of freshwater release evenly across the Labrador Sea to a more precise injection along the western boundary of the North Atlantic (the coastal region of LIS), show significant differences in the ocean circulation response, as the former produces a major decline of the AMOC and the latter yields no obvious effect on the strength of the thermohaline circulation. Furthermore, proglacial drainage of Lakes Agassiz and Ojibway leads to a fast spin-down of the AMOC, followed, however, by a gradual recovery. Most hosing experiments lead to a warming over the Nordic Sea and Barents Sea of varying magnitudes, because of an enhanced inflow from lower latitudes and a northward displacement of the North Atlantic deep convection. These processes exist in both of our high- and low-resolution experiments, but with some local discrepancies such as (1) the hosing-induced subpolar warming is much less pronounced in the high-resolution simulations; (2) LIS coastal melting in the high-resolution model leads to a slight decrease in the AMOC; and (3) the convection formation site in the low- and high-resolution experiments differs, in the former mainly over northeastern North Atlantic Ocean, but in the latter over a very shallow subpolar region along the northern edge of the North Atlantic Ocean. In conclusion, we find that our simulations capture spatially heterogeneous responses of the early-Holocene climate.

Description: The Atlantic meridional overturning circulation (AMOC) represents the zonally integrated stream function of meridional volume transport in the Atlantic Basin. The AMOC plays an important role in transporting heat meridionally in the climate system. Observations suggest a heat transport by the AMOC of 1.3 PW at 26°N—a latitude which is close to where the Atlantic northward heat transport is thought to reach its maximum. This shapes the climate of the North Atlantic region as we know it today. In recent years there has been significant progress both in our ability to observe the AMOC in nature and to simulate it in numerical models. Most previous modeling investigations of the AMOC and its impact on climate have relied on models with horizontal resolution that does not resolve ocean mesoscale eddies and the dynamics of the Gulf Stream/North Atlantic Current system. As a result of recent increases in computing power, models are now being run that are able to represent mesoscale ocean dynamics and the circulation features that rely on them. The aim of this review is to describe new insights into the AMOC provided by high-resolution models. Furthermore, we will describe how high-resolution model simulations can help resolve outstanding challenges in our understanding of the AMOC.

Description: Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18–17 ka and 6–5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland’s climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately.

Description: 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.

Description: Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values for surface air temperature and precipitation are very similar to the day-length adjusted values, and therefore correcting the calendar based on the monthly model results can largely reduce the artificial bias. In addition, we examine the calendar effects in three transient simulations for 6–0 ka by AWI-ESM, MPI-ESM, and IPSL-CM. We find significant discrepancies between adjusted and unadjusted temperature values over continents for both hemispheres in boreal autumn, while for other seasons the deviations are relatively small. A drying bias can be found in the summer monsoon precipitation in Africa (in the classical calendar), whereby the magnitude of bias becomes smaller over time. Overall, our study underlines the importance of the application of calendar transformation in the analysis of climate simulations. Neglecting the calendar effects could lead to a profound artificial distortion of the calculated seasonal cycle of surface air temperature and precipitation.

Description: We estimate 3-dimensional ocean currents from hydrographic (ARGO) data by determining the associated circulation in an inverse model. While velocities are treated diagnostic as an instantaneous steady state response temperature and salinity are allowed to change slowly with their inter-annual variability. Annual mean solutions are presented for 1999 to 2008. Altimetry referenced to a geoid provides a mean dynamic topography that determines the large scale surface circulation. AGRO data extend this information further into the ocean. It appeared useful to regularize the solution by constraining deep velocities to be small or as in our case to be close to a prognostic model simulation. Altimetry alone already improves temperature and salinity fields while ARGO data are less useful in constraining the dynamic topography. Heat, volume and overturning transports are in general agreement with previous work. Their inter-annual variability appears to be large in comparison to possible trends. Transport variances are estimated by perturbing the input data in a Monte Carlo simulation. They are smaller than the changes between consecutive years. ARGO data coverage seems to be reliable for this work after the year 2002.

Description: A new climate model supporting multi-resolution meshes in the ocean component has been established at the Alfred Wegener Institute (AWI) in Bremerhaven. The atmospheric component is ECHAM6 with T63L47 setting, while the ocean is simulated by the AWI multi-resolution model FESOM, supporting triangular unstructured meshes. Two multi-century simulations with ECHAM6-FESOM, REF and TRO, document the beneficial role of an increased tropical ocean resolution for ENSO simulations. REF features a tropical ocean resolution of about 1°, TRO employs more than 0.25° in a narrow equatorial band, with resolution gradually decreasing to 1° as in REF. Outside the tropical belt (15°N to 15°S), both meshes are identical. REF and TRO simulate a mean climate comparable to some of the best CMIP5 models. In TRO, however, both the cold tongue SST bias and the western Pacific SST standard deviation bias appear to improve along with the Nino-3 index statistics. Also, advanced ENSO diagnostics including the Nino-3.4 seasonal variance, the annual cycle representation, and its interaction with ENSO tend to improve. The robustness of these improvements is analyzed and their physical explanations are explored.