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: Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semi-annual variability is also pronounced, despite weak forcing at that period. Here we use multi-year, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semi-annual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the 4th mode and the semi-annual cycle by the 2nd mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semi-annual cycles of the respective dominant baroclinic modes are associated with characteristic basin-wide structures. Using an idealized linear reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e. square basin versus realistic coastlines) or forcing (i.e. spatially uniform versus spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

Description: The response of eddy kinetic energy (EKE) to the strengthening of Southern Hemisphere winds occurring since the 1950s is investigated with a global ocean model having a resolution of 1/12° in the Antarctic Circumpolar Current domain. The simulations expose regional differences in the relative importance of stochastic and wind-related contributions to inter-annual EKE changes. In the Pacific and Indian sectors the model captures the EKE variability observed since 1993 and confirms previous hypotheses of a lagged response to regional wind stress anomalies. Here, the multi-decadal trend in wind stress is reflected in an increase in EKE typically exceeding 5 cm2 sec-2 decade-1. In the western Atlantic EKE variability is mostly stochastic, is weakly correlated with wind fluctuations, and its multi-decadal trends are close to zero. The non-uniform distribution of wind-related changes in the eddy activity could affect the regional patterns of ocean circulation and biogeochemical responses to future climate change.

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: • Arctic sea ice extent and solid freshwater in 14 CORE-II models are inter-compared. • The models better represent the variability than the mean state. • The September ice extent trend is reasonably represented by the model ensemble mean. • The descending trend of ice thickness is underestimated compared to observations. • The models underestimate the reduction in solid freshwater content in recent years. Abstract: The Arctic Ocean simulated in fourteen global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE II) is analyzed. The focus is on the Arctic sea ice extent, the solid freshwater (FW) sources and solid freshwater content (FWC). Available observations are used for model evaluation. The variability of sea ice extent and solid FW budget is more consistently reproduced than their mean state in the models. The descending trend of September sea ice extent is well simulated in terms of the model ensemble mean. Models overestimating sea ice thickness tend to underestimate the descending trend of September sea ice extent. The models underestimate the observed sea ice thinning trend by a factor of two. When averaged on decadal time scales, the variation of Arctic solid FWC is contributed by those of both sea ice production and sea ice transport, which are out of phase in time. The solid FWC decreased in the recent decades, caused mainly by the reduction in sea ice thickness. The models did not simulate the acceleration of sea ice thickness decline, leading to an underestimation of solid FWC trend after 2000. The common model behavior, including the tendency to underestimate the trend of sea ice thickness and March sea ice extent, remains to be improved.

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: • Arctic liquid freshwater budget simulated in 14 CORE-II models is studied. • The models better represent the temporal variability than the mean state. • Multi-model mean (MMM) FW fluxes compare well with observations. • MMM FWC shows an upward trend in the recent years, with an underestimated rate. • FW flux interannual variability is more consistent where volume flux determines it. Abstract: The Arctic Ocean simulated in 14 global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II) is analyzed in this study. The focus is on the Arctic liquid freshwater (FW) sources and freshwater content (FWC). The models agree on the interannual variability of liquid FW transport at the gateways where the ocean volume transport determines the FW transport variability. The variation of liquid FWC is induced by both the surface FW flux (associated with sea ice production) and lateral liquid FW transport, which are in phase when averaged on decadal time scales. The liquid FWC shows an increase starting from the mid-1990s, caused by the reduction of both sea ice formation and liquid FW export, with the former being more significant in most of the models. The mean state of the FW budget is less consistently simulated than the temporal variability. The model ensemble means of liquid FW transport through the Arctic gateways compare well with observations. On average, the models have too high mean FWC, weaker upward trends of FWC in the recent decade than the observation, and low consistency in the temporal variation of FWC spatial distribution, which needs to be further explored for the purpose of model development.

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

Description: The spatial and temporal variability of Antarctic Intermediate Water (AAIW) formation rates during 1948-2007 is investigated in a set of high-resolution experiments performed with the NEMO-LIM ocean sea-ice model. Two sets of experiments are analyzed: one having a global 1/4° horizontal resolution (ORCA025) and another having a resolution of 1/12° in the Southern Ocean, achieved using a two-way nesting technique. The hindcast experiments are forced with the CORE-II atmospheric reanalysis from 1948-2007, with companion climatological experiments used to correct for model spurious trends. The spatial distribution of AAIW formation rates is explored via available AAIW volume and via kinematic subduction rates. Maximum AAIW formation rates are found in the regions where the mixed layer depth (MLD) is deepest in winter, i.e. in the southeast Pacific and along the Antarctic Circumpolar Current fronts, whereas kinematic subduction rates are highest in regions dominated by lateral induction. A circumpolar multi-decadal decline of AAIW formation rates and kinematic subduction rates is found, with largest decline in the southeast Pacific. Here, a strong negative MLD trend is accompanied by freshening of the mixed layer, possibly linked to a multi-decadal increase in dominance of the freshwater versus the thermal buoyancy forcing as a response to large scale warming.