Description: Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.

Description: The Eddy Kinetic Energy (EKE) associated with the Subtropical Countercurrent (STCC) in the western subtropical South Pacific is known to exhibit substantial seasonal and decadal variability. Using an eddy-permitting ocean general circulation model, which is able to reproduce the observed, salient features of the seasonal cycles of shear, stratification, baroclinic production and the associated EKE, we investigate the decadal changes of EKE. We show that the STCC region exhibits, uniquely among the subtropical gyres of the world’s oceans, significant, atmospherically forced, decadal EKE variability. The decadal variations are driven by changing vertical shear between the STCC in the upper 300 m and the South Equatorial Current below, predominantly caused by variations in STCC strength associated with a changing meridional density gradient. In the 1970s, an increased meridional density gradient results in EKE twice as large as in later decades in the model. Utilizing sensitivity experiments, decadal variations in the wind field are shown to be the essential driver. Local wind stress curl anomalies associated with the Interdecadal Pacific Oscillation (IPO) lead to up- and downwelling of the thermocline, inducing strengthening or weakening of the STCC and the associated EKE. Additionally, remote wind stress curl anomalies in the eastern subtropical South Pacific, which are not related to the IPO, generate density anomalies that propagate westward as Rossby waves and can account for up to 30–40 % of the density anomalies in the investigated region.

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: • 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: Highlights: • A joint analysis of deep current meter records in the western North Atlantic. • Intra-seasonal variability dominates the deep boundary current. • Topographic waves near 10d periods trapped over steep topography. • Basin centers are showing longer periods (50d) caused by the eddy field. • Observed variability characteristics compared to high resolution model simulation. Abstract The Deep Western Boundary Current (DWBC) along the western margin of the subpolar North Atlantic is an important component of the deep limb of the Meridional Overturning near its northern origins. A network of moored arrays from Denmark Strait to the tail of the Grand Banks has been installed for almost two decades to observe the boundary currents and transports of North Atlantic Deep Water as part of an internationally coordinated observatory for the Atlantic Meridional Overturning Circulation. The dominant variability in all of the moored velocity time series is in the week-to-month period range. While the temporal characteristics of this variability change only gradually between Denmark Strait and Flemish Cap, a broad band of longer term variability is present farther along the path of the DWBC at the Grand Banks and in the interior basins (Labrador and Irminger Seas). The vigorous intra-seasonal variability may well mask possible interannual to decadal variability that is typically an order of magnitude smaller than the high-frequency fluctuations. Here, the intra-seasonal variability is quantified at key positions along the DWBC path using both, observations and high resolution model data. The results are used to evaluate the model circulation, and in turn the model is used to relate the discrete measurements to the overall pattern of the subpolar circulation. Topographic waves are found to be trapped by the steep topography all around the western basins, the Labrador and Irminger Seas. In the Labrador Sea, the high intra-seasonal variability of the boundary current regime is separated by a region of extremely low variability in narrow recirculation cells from the basin interior. There, the variability is also on intra-seasonal timescales, but at much longer periods around 50 days.

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