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: A set of experiments utilizing different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-) tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Niño3-region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-) decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Niño episodes. The behaviour is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Description: Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.

Description: Model drift in the Labrador Sea in eddy permitting model simulations is examined using a series of configurations based on the NEMO numerical framework. There are two phases of the drift that we can identify, beginning with an initial rapid 3-year period, associated with the adjustment of the model from its initial conditions followed by an extended model drift/adjustment that continued for at least another decade. The drift controlled the model salinity in the Labrador Sea, over-riding the variability. Thus, during this initial period, similar behavior was observed between the inter-annually forced experiments as with perpetual year forcing. The results also did not depend on whether the configuration was global, or regional North Atlantic Ocean. The inclusion of an explicit sea-ice component did not seem to have a significant impact on the interior drift. Clear cut evidence for the drift having an advective nature was shown, based on two separate currents/flow regimes. We find, as expected, the representation of freshwater in the sub-polar gyre’s boundary currents important. But this study also points out another, equally important process and pathway: the input of high salinity mode water from the subtropical North Atlantic. The advective regime is dependent on the details of the model, such as the representation of the freshwater transport in the model’s East Greenland Current being very sensitive to the strength of the local sea surface salinity restoring (and the underlying field that the model is being restored to).

Description: The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

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.

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: • Mid-Atlantic vent mussel populations are contemporarily isolated • Population connectivity can only be maintained in a stepwise manner • Four mussel lineages exist on the Mid-Atlantic Ridge • Recolonization of perturbed vent localities is uncertain Summary: Deep-sea hydrothermal vents are patchily distributed ecosystems inhabited by specialized animal populations that are textbook meta-populations. Many vent-associated species have free-swimming, dispersive larvae that can establish connections between remote populations. However, connectivity patterns among hydrothermal vents are still poorly understood because the deep sea is undersampled, the molecular tools used to date are of limited resolution, and larval dispersal is difficult to measure directly. A better knowledge of connectivity is urgently needed to develop sound environmental management plans for deep-sea mining. Here, we investigated larval dispersal and contemporary connectivity of ecologically important vent mussels (Bathymodiolus spp.) from the Mid-Atlantic Ridge by using high-resolution ocean modeling and population genetic methods. Even when assuming a long pelagic larval duration, our physical model of larval drift suggested that arrival at localities more than 150 km from the source site is unlikely and that dispersal between populations requires intermediate habitats (“phantom” stepping stones). Dispersal patterns showed strong spatiotemporal variability, making predictions of population connectivity challenging. The assumption that mussel populations are only connected via additional stepping stones was supported by contemporary migration rates based on neutral genetic markers. Analyses of population structure confirmed the presence of two southern and two hybridizing northern mussel lineages that exhibited a substantial, though incomplete, genetic differentiation. Our study provides insights into how vent animals can disperse between widely separated vent habitats and shows that recolonization of perturbed vent sites will be subject to chance events, unless connectivity is explicitly considered in the selection of conservation areas.