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: 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: 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: Benthic storms are important for both the energy budget of the ocean and for sediment resuspension and transport. Using 30 years of output from a high-resolution model of the North Atlantic, it is found that most of the benthic storms in the model occur near the western boundary in association with the Gulf Stream and the North Atlantic Current, in regions that are generally co-located with the peak near-bottom eddy kinetic energy. A common feature are meander troughs in the near-surface jets that are accompanied by deep low pressure anomalies spinning up deep cyclones with near-bottom velocities of up to more than 0.5 m/s. A case study of one of these events shows the importance of both baroclinic and barotropic instability of the jet, with energy being extracted from the jet in the upstream part of the meander trough and partly returned to the jet in the downstream part of the meander trough. This motivates examining the 30-year time mean of the energy transfer from the (annual mean) background flow into the eddy kinetic energy. This quantity is shown to be co-located well with the region in which benthic storms and large increases in deep cyclonic relative vorticity occur most frequently, suggesting an important role for mixed barotropic-baroclinic instability driven cyclogenesis in generating benthic storms throughout the model simulation. Regions of largest energy transfer and most frequent benthic storms are found to be the Gulf Stream west of the New England Seamounts and the North Atlantic Current near Flemish Cap.

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: The Agulhas Current plays a crucial role in the thermohaline circulation through its leakage into the South Atlantic. Under both past and present climates, the trade winds and westerlies could have the ability to modulate the amount of Indian-Atlantic inflow. Compelling arguments have been put forward suggesting that trade winds alone have little impact on the magnitude of Agulhas leakage. Here, employing three ocean models for robust analysis – a global coarse resolution, a regional eddy-permitting and a nested high-resolution eddy-resolving configuration – and systematically altering the position and intensity of the westerly wind belt in a series of sensitivity experiments, it is shown that the westerlies, in particular their intensity, control the leakage. Leakage responds proportionally to the westerlies intensity up to a certain point. Beyond this, through the adjustment of the large-scale circulation, energetic interactions occur between the Agulhas Return Current and the Antarctic Circumpolar Current that result in a state where leakage no longer increases. This adjustment takes place within 1 to 2 decades. Contrary to previous assertions, our results further show that an equatorward (poleward) shift in westerlies increases (decreases) leakage. This occurs due to the redistribution of momentum input by the winds. It is concluded that the reported present-day leakage increase could therefore reflect an unadjusted oceanic response mainly to the strengthening westerlies over the last few decades.

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.