Description: The Lagrangian analysis of sets of particles advected with the flow fields of ocean models are used to study connectivity, i.e. exchange pathways, timescales and volume transports, between distinct oceanic regions. One important factor influencing the dispersion of fluid particles and hence connectivity is the Lagrangian eddy diffusivity, which quantifies the influence of turbulent processes on the rate of particle dispersal. Due to spatial and temporal discretization, turbulence is not fully resolved in modelled velocities, and the concept of eddy diffusivity is used to parametrize the impact of unresolved processes. However, the relations between observational- and model-based Lagrangian eddy diffusivity estimates as well as eddy parameterizations are not clear. This study presents an analysis of the spatially variable near-surface lateral eddy diffusivity estimates obtained from Lagrangian trajectories simulated with 5-day mean velocities from an eddy-resolving ocean model (INALT01) for the Agulhas system. INALT01 features diffusive regimes for dynamically different regions, some of which exhibit strong suppression of eddy mixing by mean flow, and is consistent with the pattern and magnitude of drifter-based eddy diffusivity estimates. Using monthly-mean velocities decreases the estimated diffusivities less than eddy kinetic energy, supporting the idea that large and persistent eddy features dominate eddy diffusivities. For a non-eddying ocean model (ORCA05), Lagrangian eddy diffusivities are greatly reduced, in particular when the Gent and McWilliams parameterization of mesoscale eddies is employed.

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: 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: The Indian Ocean has sustained robust surface warming in recent decades, but the role of multi-decadal variability remains unclear. Using ocean model hindcasts, characteristics of low-frequency Indian Ocean temperature variations are explored. Simulated upper-ocean temperature changes across the Indian Ocean in the hindcast are consistent with those recorded in observational products and ocean reanalyses. Indian Ocean temperatures exhibit strong warming trends since the 1950s limited to the surface and south of 30°S, while extensive subsurface cooling occurs over much of the tropical Indian Ocean. Previous work focused on diagnosing causes of these long-term trends in the Indian Ocean over the second half of the 20th Century. Instead, the temporal evolution of Indian Ocean subsurface heat content is shown here to reveal distinct multi-decadal variations associated with the Pacific Decadal Oscillation and the long-term trends are thus interpreted to result from aliasing of the low-frequency variability. Transmission of the multi-decadal signal occurs via an oceanic pathway through the Indonesian Throughflow and is manifest across the Indian Ocean centered along 12°S as westward propagating Rossby waves modulating thermocline and subsurface heat content variations. Resulting low-frequency changes in the eastern Indian Ocean thermocline depth are associated with decadal variations in the frequency of Indian Ocean Dipole (IOD) events, with positive IOD events unusually common in the 1960s and 1990s with a relatively shallow thermocline. In contrast, the deeper thermocline depth in the 1970s and 1980s is associated with frequent negative IOD and rare positive IOD events. Changes in Pacific wind forcing in recent decades and associated rapid increases in Indian Ocean subsurface heat content can thus affect the basin’s leading mode of variability, with implications for regional climate and vulnerable societies in surrounding countries.

Description: The upper tropical Atlantic Ocean has markedly warmed since the 1960s. It has been shown that this warming was not due to local heat fluxes, and that the trade winds that drive the coastal and equatorial upwelling have intensified rather than weakened. Remote forcing might thus have played an important role. Here model experiments are used to investigate the contribution from an increased inflow of warm Indian Ocean water through Agulhas leakage. A high-resolution hindcast experiment with interannually varying forcing for the time period 1948 to 2007, in which Agulhas leakage increases by about 45% from the 1960s to the early 2000s, reproduces the observed warming trend. To tease out the role of Agulhas leakage, a sensitivity experiment designed to only increase Agulhas leakage is used. Compared to a control simulation it shows a pronounced warming in the upper tropical Atlantic Ocean. A Lagrangian trajectory analysis confirms that a significant portion of Agulhas leakage water reaches the upper 300m of the tropical Atlantic Ocean within two decades, and that the tropical Atlantic warming in the sensitivity experiment is mainly due to water of Agulhas origin. Therefore, it is suggested that the increased trade winds since the 1960s favor upwelling of warmer subsurface waters, which in parts originate from the Agulhas, leading to higher SSTs in the tropics

Description: To model tracer spreading in the ocean, Lagrangian simulations in an offline framework are a practical and efficient alternative to solving the advective–diffusive tracer equations online. Differences in both approaches raise the question of whether both methods are comparable. Lagrangian simulations usually use model output averaged in time, and trajectories are not subject to parameterized subgrid diffusion, which is included in the advection–diffusion equations of ocean models. Previous studies focused on diffusivity estimates in idealized models but could show that both methods yield similar results as long as the deformations-scale dynamics are resolved and a sufficient amount of Lagrangian particles is used. This study compares spreading of an Eulerian tracer simulated online and a cloud of Lagrangian particles simulated offline with velocities from the same ocean model. We use a global, eddy-resolving ocean model featuring 1/20° horizontal resolution in the Agulhas region around South Africa. Tracer and particles were released at one time step in the Cape Basin and below the mixed layer and integrated for 3 years. Large-scale diagnostics, like mean pathways of floats and tracer, are almost identical and 1D horizontal distributions show no significant differences. Differences in vertical distributions, seen in a reduced vertical spreading and downward displacement of particles, are due to the combined effect of unresolved subdaily variability of the vertical velocities and the spatial variation of vertical diffusivity. This, in turn, has a small impact on the horizontal spreading behavior. The estimates of eddy diffusivity from particles and tracer yield comparable results of about 4000 m2 s−1 in the Cape Basin.

Description: Satellite observations and output from a high-resolution ocean model are used to investigate how the Loop Current in the Gulf of Mexico affects the Gulf Stream transport through the Florida Straits. We find that the expansion (contraction) of the Loop Current leads to lower (higher) transports through the Straits of Florida. The associated surface velocity anomalies are coherent from the southwestern tip of Florida to Cape Hatteras. A simple continuity-based argument can be used to explain the link between the Loop Current and the downstream Gulf Stream transport: as the Loop Current lengthens (shortens) its path in the Gulf of Mexico, the flow out of the Gulf decreases (increases). Anomalies in the surface velocity field are first seen to the southwest of Florida and within 4 weeks propagate through the Florida Straits up to Cape Hatteras and into the Gulf Stream Extension. In both the observations and the model this propagation can be seen as pulses in the surface velocities. We estimate that the Loop Current variability can be linked to a variability of several Sverdrups (1Sv = 10(6) m(3) s(-1)) through the Florida Straits. The exact timing of the Loop Current variability is largely unpredictable beyond a few weeks and its variability is therefore likely a major contributor to the chaotic/intrinsic variability of the Gulf Stream. However, the time lag between the Loop Current and the flow downstream of the Gulf of Mexico means that if a lengthening/shortening of the Loop Current is observed this introduces some predictability in the downstream flow for a few weeks.