Description: We present a new surface-atmospheric dataset for driving ocean–sea-ice models based on Japanese 55-year atmospheric reanalysis (JRA-55), referred to here as JRA55-do. The JRA55-do dataset aims to replace the CORE interannual forcing version 2 (hereafter called the CORE dataset), which is currently used in the framework of the Coordinated Ocean-ice Reference Experiments (COREs) and the Ocean Model Intercomparison Project (OMIP). A major improvement in JRA55-do is the refined horizontal grid spacing (∼ 55 km) and temporal interval (3 hr). The data production method for JRA55-do essentially follows that of the CORE dataset, whereby the surface fields from an atmospheric reanalysis are adjusted relative to reference datasets. To improve the adjustment method, we use high-quality products derived from satellites and from several other atmospheric reanalysis projects, as well as feedback on the CORE dataset from the ocean modelling community. Notably, the surface air temperature and specific humidity are adjusted using multi-reanalysis ensemble means. In JRA55-do, the downwelling radiative fluxes and precipitation, which are affected by an ambiguous cloud parameterisation employed in the atmospheric model used for the reanalysis, are based on the reanalysis products. This approach represents a notable change from the CORE dataset, which imported independent observational products. Consequently, the JRA55-do dataset is more self-contained than the CORE dataset, and thus can be continually updated in near real-time. The JRA55-do dataset extends from 1958 to the present, with updates expected at least annually. This paper details the adjustments to the original JRA-55 fields, the scientific rationale for these adjustments, and the evaluation of JRA55-do. The adjustments successfully corrected the biases in the original JRA-55 fields. The globally averaged features are similar between the JRA55-do and CORE datasets, implying that JRA55-do can suitably replace the CORE dataset for use in driving global ocean–sea-ice models.

Description: Highlights: • We focus on ACC and Southern Ocean MOC during 1958–2007 in 17 CORE-II forced models. • Most CORE-II simulations are close to eddy saturation. • Most CORE-II simulations are far from showing signs of eddy compensation. • Constant in time or space k results in poor representation of mesoscale eddy effects. • MOC has larger sensitivity than ACC transport even in eddy saturated state. Abstract: In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean meridional overturning circulation (MOC) in a suite of seventeen global ocean–sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 1958–2007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variabilities are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient κ show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying κ show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space κ is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.

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 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: 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: 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: We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.

Description: The Earth system is accumulating energy due to human-induced activities. More than 90% of this energy has been stored in the ocean as heat since 1970, with similar to 60% of that in the upper 700 m. Differences in upper-ocean heat content anomaly (OHCA) estimates, however, exist. Here, we use a dataset protocol for 1970-2008-with six instrumental bias adjustments applied to expendable bathythermograph (XBT) data, and mapped by six research groups-to evaluate the spatiotemporal spread in upper OHCA estimates arising from two choices: 1) those arising from instrumental bias adjustments and 2) those arising from mathematical (i.e., mapping) techniques to interpolate and extrapolate data in space and time. We also examined the effect of a common ocean mask, which reveals that exclusion of shallow seas can reduce global OHCA estimates up to 13%. Spread due to mapping method is largest in the Indian Ocean and in the eddy-rich and frontal regions of all basins. Spread due to XBT bias adjustment is largest in the Pacific Ocean within 30 degrees N-30 degrees S. In both mapping and XBT cases, spread is higher for 1990-2004. Statistically different trends among mapping methods are found not only in the poorly observed Southern Ocean but also in the well-observed northwest Atlantic. Our results cannot determine the best mapping or bias adjustment schemes, but they identify where important sensitivities exist, and thus where further understanding will help to refine OHCA estimates. These results highlight the need for further coordinated OHCA studies to evaluate the performance of existing mapping methods along with comprehensive assessment of uncertainty estimates.