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: • 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: • Assessment of the Indian Ocean simulation from global forced sea- ice models. • SST biases are 2 times smaller in forced simulations than the coupled simulations. • Coupled model shows large inter-model spread over the eastern equatorial Indian Ocean. • Refinement in model horizontal resolution does not significantly improve simulations. • Uncover a secondary pathway of northward cross-equatorial transport along 75 °E. • Models are unable to capture the observed thick barrier layer in the north Bay of Bengal. Abstract: We present an analysis of annual and seasonal mean characteristics of the Indian Ocean circulation and water masses from 16 global ocean–sea-ice model simulations that follow the Coordinated Ocean-ice Reference Experiments (CORE) interannual protocol (CORE-II). All simulations show a similar large-scale tropical current system, but with differences in the Equatorial Undercurrent. Most CORE-II models simulate the structure of the Cross Equatorial Cell (CEC) in the Indian Ocean. We uncover a previously unidentified secondary pathway of northward cross-equatorial transport along 75 °E, thus complementing the pathway near the Somali Coast. This secondary pathway is most prominent in the models which represent topography realistically, thus suggesting a need for realistic bathymetry in climate models. When probing the water mass structure in the upper ocean, we find that the salinity profiles are closer to observations in geopotential (level) models than in isopycnal models. More generally, we find that biases are model dependent, thus suggesting a grouping into model lineage, formulation of the surface boundary, vertical coordinate and surface salinity restoring. Refinement in model horizontal resolution (one degree versus degree) does not significantly improve simulations, though there are some marginal improvements in the salinity and barrier layer results. The results in turn suggest that a focus on improving physical parameterizations (e.g. boundary layer processes) may offer more near-term advances in Indian Ocean simulations than refined grid resolution.

Description: Many models of the large-scale thermohaline circulation in the ocean exhibit strong zonally integrated upwelling in the midlatitude North Atlantic that significantly decreases the amount of deep water that is carried from the formation regions in the subpolar North Atlantic toward low latitudes and across the equator. In an analysis of results from the Community Modeling Effort using a suite of models with different horizontal resolution, wind and thermohaline forcing, and mixing parameters, it is shown that the upwelling is always concentrated in the western boundary layer between roughly 30° and 40°N. The vertical transport across 1000 m appears to be controlled by local dynamics and strongly depends on the horizontal resolution and mixing parameters of the model. It is suggested that in models with a realistic deep-water formation rate in the subpolar North Atlantic, the excessive upwelling can be considered as the prime reason for the typically too low meridional overturning rates and northward heat transports in the subtropical North Atlantic. A new isopycnal advection and mixing parameterization of tracer transports by mesoscale eddies yield substantial improvements in these integral measures of the circulation.

Description: Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations. The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.

Description: A new field of study, "decadal prediction," is emerging in climate science. Decadal prediction lies between seasonal/interannual forecasting and longer-term climate change projections, and focuses on time-evolving regional climate conditions over the next 10-30 yr. Numerous assessments of climate information user needs have identified this time scale as being important to infrastructure planners, water resource managers, and many others. It is central to the information portfolio required to adapt effectively to and through climatic changes. At least three factors influence time-evolving regional climate at the decadal time scale: 1) climate change commitment (further warming as the coupled climate system comes into adjustment with increases of greenhouse gases that have already occurred), 2) external forcing, particularly from future increases of greenhouse gases and recovery of the ozone hole, and 3) internally generated variability. Some decadal prediction skill has been demonstrated to arise from the first two of these factors, and there is evidence that initialized coupled climate models can capture mechanisms of internally generated decadal climate variations, thus increasing predictive skill globally and particularly regionally. Several methods have been proposed for initializing global coupled climate models for decadal predictions, all of which involve global time-evolving threedimensional ocean data, including temperature and salinity. An experimental framework to address decadal predictability/prediction is described in this paper and has been incorporated into the coordinated Coupled Model Intercomparison Model, phase 5 (CMIP5) experiments, some of which will be assessed for the IPCC Fifth Assessment Report (AR5). These experiments will likely guide work in this emerging field over the next 5 yr.

Description: The effects of spatial variations of the thickness diffusivity (K) appropriate to the parameterisation of [Gent, P.R. and McWilliams, J.C., 1990. Isopycnal mixing in ocean circulation models. J. Phys. Oceanogr., 20, 150–155.] are assessed in a coarse resolution global ocean general circulation model. Simulations using three closures yielding different lateral and/or vertical variations in K are compared with a simulation using a constant value. Although the effects of changing K are in general small and all simulations remain biased compared to observations, we find systematic local sensitivities of the simulated circulation on K. In particular, increasing K near the surface in the tropical ocean lifts the depth of the equatorial thermocline, the strength of the Antarctic Circumpolar Current decreases while the subpolar and subtropical gyre transports in the North Atlantic increase by increasing K locally. We also find that the lateral and vertical structure of K given by a recently proposed closure reduces the negative temperature biases in the western North Atlantic by adjusting the pathways of the Gulf Stream and the North Atlantic Current to a more realistic position.

Description: We present an analysis of annual and seasonal mean characteristics of the Indian Ocean circulation and water masses from 16 global ocean–sea-ice model simulations that follow the Coordinated Ocean-ice Reference Experiments (CORE) interannual protocol (CORE-II). All simulations show a similar large-scale tropical current system, but with differences in the Equatorial Undercurrent. Most CORE-II models simulate the structure of the Cross Equatorial Cell (CEC) in the Indian Ocean. We uncover a previously unidentified secondary pathway of northward cross-equatorial transport along 75 °E, thus complementing the pathway near the Somali Coast. This secondary pathway is most prominent in the models which represent topography realistically, thus suggesting a need for realistic bathymetry in climate models. When probing the water mass structure in the upper ocean, we find that the salinity profiles are closer to observations in geopotential (level) models than in isopycnal models. More generally, we find that biases are model dependent, thus suggesting a grouping into model lineage, formulation of the surface boundary, vertical coordinate and surface salinity restoring. Refinement in model horizontal resolution (one degree versus degree) does not significantly improve simulations, though there are some marginal improvements in the salinity and barrier layer results. The results in turn suggest that a focus on improving physical parameterizations (e.g. boundary layer processes) may offer more near-term advances in Indian Ocean simulations than refined grid resolution.