Description: Highlights: • We compare the simulated Arctic Ocean in 15 global ocean–sea ice models. • There is a large spread in temperature bias in the Arctic Ocean between the models. • Warm bias models have a strong temperature anomaly of inflow of Atlantic Water. • Dense outflows formed on Arctic shelves are not captured accurately in the models. In this paper we compare the simulated Arctic Ocean in 15 global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II). Most of these models are the ocean and sea-ice components of the coupled climate models used in the Coupled Model Intercomparison Project Phase 5 (CMIP5) experiments. We mainly focus on the hydrography of the Arctic interior, the state of Atlantic Water layer and heat and volume transports at the gateways of the Davis Strait, the Bering Strait, the Fram Strait and the Barents Sea Opening. We found that there is a large spread in temperature in the Arctic Ocean between the models, and generally large differences compared to the observed temperature at intermediate depths. Warm bias models have a strong temperature anomaly of inflow of the Atlantic Water entering the Arctic Ocean through the Fram Strait. Another process that is not represented accurately in the CORE-II models is the formation of cold and dense water, originating on the eastern shelves. In the cold bias models, excessive cold water forms in the Barents Sea and spreads into the Arctic Ocean through the St. Anna Through. There is a large spread in the simulated mean heat and volume transports through the Fram Strait and the Barents Sea Opening. The models agree more on the decadal variability, to a large degree dictated by the common atmospheric forcing. We conclude that the CORE-II model study helps us to understand the crucial biases in the Arctic Ocean. The current coarse resolution state-of-the-art ocean models need to be improved in accurate representation of the Atlantic Water inflow into the Arctic and density currents coming from the shelves.

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: This paper develops an analysis framework to identify how physical processes, as represented in ocean climate models, impact the evolution of global mean sea level. The formulation utilizes the coarse grained equations appropriate for an ocean model, and starts from the vertically integrated mass conservation equation in its Lagrangian form. Global integration of this kinematic equation results in an evolution equation for global mean sea level that depends on two physical processes: boundary fluxes of mass and the non-Boussinesq steric effect. The non-Boussinesq steric effect itself contains contributions from boundary fluxes of buoyancy; interior buoyancy changes associated with parameterized subgrid scale processes; and motion across pressure surfaces. The non-Boussinesq steric effect can be diagnosed in either volume conserving Boussinesq or mass conserving non-Boussinesq ocean circulation models, with differences found to be negligible. We find that surface heating is the dominant term affecting sea level arising from buoyancy fluxes, contributing to a net positive tendency to global mean sea level, largely due to low latitude heating and because the thermal expansion coefficient is much larger in the tropics than high latitudes. Subgrid scale effects from parameterized quasi-Stokes transport, vertical diffusion, cabbeling, and thermobaricity are also found to be significant, each resulting in a reduction of global mean sea level. Sea level rise through low latitude heating is largely compensated by a sea level drop from poleward eddy heat transport and ocean mixing. Spatial variations in the thermal expansion coefficient provide an essential modulation of how physical effects from mixing and eddy induced advective transport impact global mean sea level. Highlights: ► Theoretical framework for how physical processes impact global mean sea level in ocean models. ► Mathematical and physical specification of the non-Boussinesq steric effect. ► How boundary buoyancy fluxes and interior processes impact global mean sea level. ► Global model examples of the non-Boussinesq steric effect with associated budget for global mean sea level.

Description: We characterise the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988-2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the interannual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available. The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/warm and salty water column biases centred near 50 degrees S. Over the 1933-2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200-700 in ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice

Description: The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations

Description: The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlantic

Description: Highlights: • Inter-annual to decadal variability in AMOC from CORE-II simulations is presented. • AMOC variability shows three stages, with maximum transports in mid- to late-1990s. • North Atlantic temporal variability features are in good agreement among simulations. • Such agreements suggest variability is dictated by the atmospheric data sets. • Simulations differ in spatial structures of variability due to ocean dynamics. Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958–2007 period from twenty global ocean – sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958–2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their temporal representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.

Description: Highlights: • Mean circulation patterns are assessed and Kuroshio transport is underestimated. • Water mass distribution is compared and analyzed within COREII models. • Main biases of deep MLDs result from the inaccurate Kuroshio separation. • Reasonable modeled tropical dynamics but a discrepancy from the surface wind. Abstract: We evaluate the mean circulation patterns, water mass distributions, and tropical dynamics of the North and Equatorial Pacific Ocean based on a suite of global ocean-sea ice simulations driven by the CORE-II atmospheric forcing from 1963-2007. The first three moments (mean, standard deviation and skewness) of sea surface height and surface temperature variability are assessed against observations. Large discrepancies are found in the variance and skewness of sea surface height and in the skewness of sea surface temperature. Comparing with the observation, most models underestimate the Kuroshio transport in the Asian Marginal seas due to the missing influence of the unresolved western boundary current and meso-scale eddies. In terms of the Mixed Layer Depths (MLDs) in the North Pacific, the two observed maxima associated with Subtropical Mode Water and Central Mode Water formation coalesce into a large pool of deep MLDs in all participating models, but another local maximum associated with the formation of Eastern Subtropical Mode Water can be found in all models with different magnitudes. The main model bias of deep MLDs results from excessive Subtropical Mode Water formation due to inaccurate representation of the Kuroshio separation and of the associated excessively warm and salty Kuroshio water. Further water mass analysis shows that the North Pacific Intermediate Water can penetrate southward in most models, but its distribution greatly varies among models depending not only on grid resolution and vertical coordinate but also on the model dynamics. All simulations show overall similar large scale tropical current system, but with differences in the structures of the Equatorial Undercurrent. We also confirm the key role of the meridional gradient of the wind stress curl in driving the equatorial transport, leading to a generally weak North Equatorial Counter Current in all models due to inaccurate CORE-II equatorial wind fields. Most models show a larger interior transport of Pacific subtropical cells than the observation due to the overestimated transport in the Northern Hemisphere likely resulting from the deep pycnocline