Description: Highlights: • Lagrangian ocean analysis is a powerful way to analyse the output of ocean circulation models • We present a review of the Kinematic framework, available tools, and applications of Lagrangian ocean analysis • While there are unresolved questions, the framework is robust enough to be used widely in ocean modelling Abstract: Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. Over several decades, a variety of tools and methods for this purpose have emerged. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolved physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. The overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.

Description: Highlights: • Comparison of global NEMO and FESOM configurations with emphasis on the Agulhas system. • Both models simulate a reasonable and comparable large-scale circulation. • Both models have individual strengths and weaknesses to match the observations of the WBC system. • The numerical cost of FESOM is twice the one of NEMO. Abstract: Many questions in ocean and climate modelling require the combined use of high resolution, global coverage and multi-decadal integration length. For this combination, even modern resources limit the use of traditional structured-mesh grids. Here we compare two approaches: A high-resolution grid nested into a global model at coarser resolution (NEMO with AGRIF) and an unstructured-mesh grid (FESOM) which allows to variably enhance resolution where desired. The Agulhas system around South Africa is used as a testcase, providing an energetic interplay of a strong western boundary current and mesoscale dynamics. Its open setting into the horizontal and global overturning circulations also requires global coverage. Both model configurations simulate a reasonable large-scale circulation. Distribution and temporal variability of the wind-driven circulation are quite comparable due to the same atmospheric forcing. However, the overturning circulation differs, owing each model's ability to represent formation and spreading of deep water masses. In terms of regional, high-resolution dynamics, all elements of the Agulhas system are well represented. Owing to the strong nonlinearity in the system, Agulhas Current transports of both configurations and in comparison with observations differ in strength and temporal variability. Similar decadal trends in Agulhas Current transport and Agulhas leakage are linked to the trends in wind forcing. Although the number of 3D wet grid points used in FESOM is similar to that in the nested NEMO, FESOM uses about two times the number of CPUs to obtain the same model throughput (in terms of simulated model years per day). This is feasible due to the high scalability of the FESOM code.

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: 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: Highlights: • OGCM simulations of the AMOC are highly sensitive to the subarctic freshwater forcing. • Trends in the simulated AMOC are linked to the salinity of the DSOW. • DSOW salinity trends can be traced back to the freshwater transport by the NAC. • The NAC freshwater budget is highly affected by the salinity restoring used in OGCMs. • Modifications in the subarctic precipitation can help to minimize the restoring flux. Global ocean sea-ice models with an atmospheric forcing based on bulk formulations of the air-sea fluxes exhibit spurious trends in key flow indices like the Atlantic Meridional Overturning Circulation (AMOC), constraining their use in investigations of multi-decadal ocean variability. To identify the critical model factors affecting the temporal evolution of the AMOC on time scales of up to 60 years, a series of experiments with both eddy-permitting (0.25°) and non-eddying (0.5°) ocean-ice models has been performed, focusing on the influence of artificial choices for the freshwater forcing, in particular the restoring of sea surface salinity towards climatological values. The atmospheric forcing builds on the proposal for Coordinated Ocean-ice Reference Experiments (CORE), utilizing the refined atmospheric reanalysis products for 1948–2006 compiled by Large and Yeager. Sensitivity experiments with small variations in precipitation (within the observational uncertainty) and sea surface salinity restoring in the subarctic Atlantic produce a wide range of AMOC transports, between upward drifts to more than 22 Sv and nearly-collapsed states with less than 7 Sv, reflecting the excessive role of the salinity feedback in such simulations. In all cases the AMOC is tightly related to the density of the Denmark Strait overflow; changes in that density are governed by the salinity in the Nordic Seas; and in turn, that salinity is strongly affected by the properties of the inflowing North Atlantic water.

Description: Highlights: • We combine high-resolution ocean models with population genetics • Variation in wind-driven ocean currents mediates the collapse of A. anguilla • Female eels are philopatric within the Sargasso Sea, while males maintain gene flow • We present first evidence of the role of ocean currents in shaping species’ evolution Summary: Worldwide, exploited marine fish stocks are under threat of collapse [1]. Although the drivers behind such collapses are diverse, it is becoming evident that failure to consider evolutionary processes in fisheries management can have drastic consequences on a species’ long-term viability [2]. The European eel (Anguilla anguilla; Linnaeus, 1758) is no exception: not only does the steep decline in recruitment observed in the 1980s [ 3 and 4] remain largely unexplained, the punctual detection of genetic structure also raises questions regarding the existence of a single panmictic population [ 5, 6 and 7]. With its extended Transatlantic dispersal, pinpointing the role of ocean dynamics is crucial to understand both the population structure and the widespread decline of this species. Hence, we combined dispersal simulations using a half century of high-resolution ocean model data with population genetics tools. We show that regional atmospherically driven ocean current variations in the Sargasso Sea were the major driver of the onset of the sharp decline in eel recruitment in the beginning of the 1980s. The simulations combined with genotyping of natural coastal eel populations furthermore suggest that unexpected evidence of coastal genetic differentiation is consistent with cryptic female philopatric behavior within the Sargasso Sea. Such results demonstrate the key constraint of the variable oceanic environment on the European eel population.

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