Description: The North Brazil Current (NBC) constitutes a bottleneck for the mean northward return flow of the Atlantic Meridional Overturning Circulation (AMOC) in the tropical South Atlantic. Previous studies suggested a link between interannual to multidecadal NBC and AMOC transport variability and proposed to use NBC observations as an index for the AMOC. Here we use a set of hindcast, sensitivity, and perturbation experiments performed within a hierarchy of ocean general circulation models to show that decadal to multidecadal buoyancy-forced changes in the basin-scale AMOC transport indeed manifest themselves in the NBC. The relation is, however, masked by a strong interannual to decadal wind-driven gyre variability of the NBC. While questioning the NBC transport as a direct index for the AMOC, the results support its potential merit for an AMOC monitoring system, provided that the wind-driven circulation variability is properly accounted for.

Description: A new Earth system model, the Flexible Ocean and Climate Infrastructure (FOCI), is introduced. A first version of FOCI consists of a global high-top atmosphere (ECHAM6.3) and an ocean model (NEMO3.6) as well as sea ice (LIM2) and land surface model components (JSBACH), which are coupled through the OASIS3-MCT software package. FOCI includes a number of optional modules which can be activated depending on the scientific question of interest. In the atmosphere, interactive stratospheric chemistry can be used (ECHAM6-HAMMOZ) to study, for example, the effects of the ozone hole on the climate system. In the ocean, a biogeochemistry model (MOPS) is available to study the global carbon cycle. A unique feature of FOCI is the ability to explicitly resolve mesoscale ocean eddies in specific regions. This is realized in the ocean through nesting; first examples for the Agulhas Current and the Gulf Stream systems are described here. FOCI therefore bridges the gap between coarse-resolution climate models and global high-resolution weather prediction and ocean-only models. It allows to study the evolution of the climate system on regional and seasonal to (multi-) decadal scales. The development of FOCI resulted from a combination of the long-standing expertise in ocean and climate modeling in several research units and divisions at GEOMAR. FOCI will thus be used to complement and interpret long-term observations in the Atlantic, enhance the process understanding of the role of mesoscale oceanic eddies for large-scale oceanic and atmospheric circulation patterns, study feedback mechanisms with stratospheric processes, estimate future ocean acidification, improve the simulation of the Atlantic Meridional Overturning Circulation changes and their influence on climate, ocean chemistry and biology. In this paper we present both the scientific vision for the development of FOCI as well as some technical details. This includes a first validation of the different model components using several configurations of FOCI. Results show that the model in its basic configuration runs stably under pre-industrial control as well as under historical forcing, and produces a mean climate and variability which compares well with observations, reanalysis products and other climate models. The nested configurations reduce some long-standing biases in climate models and are an important step forward to include the atmospheric response in multi-decadal eddy-rich configurations.

Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.

Description: The westerlies and trade winds over the South Atlantic and Indian Ocean are important drivers of the regional oceanography around southern Africa, including features such as the Agulhas Current, the Agulhas leakage, and the Benguela upwelling. Agulhas leakage constitutes a fraction of warm and saline water transport from the Indian Ocean into the South Atlantic. The leakage is stronger during intensified westerlies. Here, we analyze the wind stress of different observational and modeled atmospheric data sets (covering the last 2 millennia, the recent decades, and the 21st century) with regard to the intensity and position of the southeasterly trades and the westerlies. The analysis reveals that variations of both wind systems go hand in hand and that a poleward shift of the westerlies and trades and an intensification of westerlies took place during the recent decades. Furthermore, upwelling in South Benguela is slightly intensified when trades are shifted poleward. Projections for strength and position of the westerlies in the 21st century depend on assumed CO2 emissions and on their effect relative to the ozone forcing. In the strongest emission scenario (RCP8.5) the simulations show a further southward displacement, whereas in the weakest emission scenario (RCP2.6) a northward shift is modeled, possibly due to the effect of ozone recovery dominating the effect of anthropogenic greenhouse forcing. We conclude that the Agulhas leakage has intensified during the last decades and is projected to increase if greenhouse gas emissions are not reduced. This will have a small impact on Benguela upwelling strength and may also have consequences for water mass characteristics in the upwelling region. An increased contribution of Agulhas water to the upwelling water masses will import more preformed nutrients and oxygen into the upwelling region.

Description: Mesoscale dynamics of the Agulhas Current system determine the exchange between the Indian and Atlantic oceans, thereby influencing the global overturning circulation. Using a series of ocean model experiments compared to observations, we show that the representation of mesoscale eddies in the Agulhas ring path improves with increasing resolution of submesoscale flows. Simulated submesoscale dynamics are validated with time‐mean horizontal‐wavenumber spectra from satellite sea surface temperature measurements and mesoscale dynamics with spectra from sea surface height. While the Agulhas ring path in a nonsubmesoscale‐resolving (1/20)° configuration is associated with too less power spectral densities on all scales and too steep spectral slopes, the representation of the mesoscale dynamics improves when the diffusion and the dissipation of the model are reduced and some small‐scale features are resolved. Realistic power spectral densities over all scales are achieved when additionally the horizontal resolution is increased to (1/60)° and a larger portion of the submesoscale spectrum is resolved. Results of an eddy detection algorithm applied to the model outputs as well as to a gridded sea surface height satellite product show that in particular strong cyclones are much better represented when submesoscale flows are resolved by the model. The validation of the submesoscale dynamics with sea surface temperature spectra provides guidance for the choice of advection schemes and explicit diffusion and dissipation as well as for further subgrid‐scale parameterizations. For the Agulhas ring path, the use of upstream biased advection schemes without explicit diffusion and dissipation is found to be associated with realistically simulated submesoscales.

Description: Key Points: • A regional ocean model is used to examine multidecadal shelf temperature changes on the Agulhas Bank • There are distinct shelf temperature regime changes in 1966 and 1996 • These regime shifts are caused by changes in coastal upwelling linked to large-scale wind variability The Agulhas Bank is an important area for the spawning of small pelagic fish and other species. Here, within a NEMO ocean model, we investigate changes in temperature over the Bank on multidecadal time scales. In agreement with previous observational studies, a shift to colder temperatures is found in 1997. The model also simulates an earlier shift from colder to warmer temperatures in 1966. These shifts are coastally confined and shown, using a climatologically forced model run as a control, to be driven by a north‐south migration in the large‐scale wind belts, rather than by changes in downward heat fluxes or changes in the Agulhas Current itself. The zonal wind changes on the Agulhas Bank show a significant relationship with the Southern Annular Mode, showing some promise for future predictability of cold and warm regimes on the Agulhas Bank. Thus, while the Agulhas Current has been shown in previous work to have a large impact on intra‐annual and interannual temperature variability, this work shows that multidecadal variability in temperature on the shelf is likely to be wind forced.

Description: Aim: Invasive species are of increasing global concern. Nevertheless, the mechanisms driving further distribution after the initial establishment of non‐native species remain largely unresolved, especially in marine systems. Ocean currents can be a major driver governing range occupancy, but this has not been accounted for in most invasion ecology studies so far. We investigate how well initial establishment areas are interconnected to later occupancy regions to test for the potential role of ocean currents driving secondary spread dynamics in order to infer invasion corridors and the source–sink dynamics of a non‐native holoplanktonic biological probe species on a continental scale. Location: Western Eurasia. Time period: 1980s–2016. Major taxa studied: ‘Comb jelly’ Mnemiopsis leidyi. Methods: Based on 12,400 geo‐referenced occurrence data, we reconstruct the invasion history of M. leidyi in western Eurasia. We model ocean currents and calculate their stability to match the temporal and spatial spread dynamics with large‐scale connectivity patterns via ocean currents. Additionally, genetic markers are used to test the predicted connectivity between subpopulations. Results: Ocean currents can explain secondary spread dynamics, matching observed range expansions and the timing of first occurrence of our holoplanktonic non‐native biological probe species, leading to invasion corridors in western Eurasia. In northern Europe, regional extinctions after cold winters were followed by rapid recolonizations at a speed of up to 2,000 km per season. Source areas hosting year‐round populations in highly interconnected regions can re‐seed genotypes over large distances after local extinctions. Main conclusions: Although the release of ballast water from container ships may contribute to the dispersal of non‐native species, our results highlight the importance of ocean currents driving secondary spread dynamics. Highly interconnected areas hosting invasive species are crucial for secondary spread dynamics on a continental scale. Invasion risk assessments should consider large‐scale connectivity patterns and the potential source regions of non‐native marine species.

Description: The northward flow of the upper limb of the Atlantic Meridional Overturning Circulation (AMOC) is fed by waters entering the South Atlantic from the Indian Ocean mainly via the Agulhas Current (AC) system and by waters entering from the Pacific through Drake Passage (DP), commonly referred to as the “warm” and “cold” water routes, respectively. However, there is no final consensus on the relative importance of these two routes for the upper limb's volume transport and thermohaline properties. In this study we revisited the AC and DP contributions by performing Lagrangian analyses between the two source regions and the North Brazil Current (NBC) at 6∘ S in a realistically forced high-resolution (1∕20∘) ocean model. Our results agree with the prevailing conception that the AC contribution is the major source for the upper limb transport of the AMOC in the tropical South Atlantic. However, they also suggest a non-negligible DP contribution of around 40 %, which is substantially higher than estimates from previous Lagrangian studies with coarser-resolution models but now better matches estimates from Lagrangian observations. Moreover, idealized analyses of decadal changes in the DP and AC contributions indicate that the ongoing increase in Agulhas leakage indeed may have induced an increase in the AC contribution to the upper limb of the AMOC in the tropics, while the DP contribution decreased. In terms of thermohaline properties, our study highlights the fact that the AC and DP contributions cannot be unambiguously distinguished by their temperature, as the commonly adopted terminology may imply, but rather by their salinity when entering the South Atlantic. During their transit towards the NBC the bulk of DP waters experiences a net density loss through a net warming, whereas the bulk of AC waters experiences a slight net density gain through a net increase in salinity. Notably, these density changes are nearly completely captured by Lagrangian particle trajectories that reach the surface mixed layer at least once during their transit, which amount to 66 % and 49 % for DP and AC waters, respectively. This implies that more than half of the water masses supplying the upper limb of the AMOC are actually formed within the South Atlantic and do not get their characteristic properties in the Pacific and Indian Oceans.

Description: Ballast water treatment is required for vessels to prevent the introduction of potentially invasive neobiota. Some treatment methods use chemical disinfectants which produce a variety of halogenated compounds as disinfection by-products (DBPs). One of the most abundant DBP from oxidative ballast water treatment is bromoform (CHBr3) where we find an average concentration of 894±560nmolL-1 (226±142μgL-1) in the undiluted ballast water from measurements and literature. Bromoform is a relevant gas for atmospheric chemistry and ozone depletion, especially in the tropics where entrainment into the stratosphere is possible. The spread of DBPs in the tropics over months to years is assessed here for the first time. With Lagrangian trajectories based on the NEMO-ORCA12 model velocity field, we simulate DBP spread in the sea surface and try to quantify the oceanic bromoform concentration and emission to the atmosphere from ballast water discharge at major harbours in the tropical region of Southeast Asia. The exemplary simulations of two important regions, Singapore and the Pearl River Delta, reveal major transport pathways of the DBPs and the anthropogenic bromoform concentrations in the sea surface. Based on our simulations, we expect DBPs to spread into the open ocean, along the coast and also an advection with monsoon-driven currents into the North Pacific and Indian Ocean. Furthermore, anthropogenic bromoform concentrations and emissions are predicted to increase locally around large harbours. In the sea surface around Singapore we estimate an increase in bromoform concentration by 9% compared to recent measurement. In a moderate scenario where 70% of the ballast water is chemically treated bromoform emissions to the atmosphere can locally exceed 1000pmolm-2h-1 and double climatological emissions. In the Pearl River Delta all bromoform is directly outgassed which leads to an additional bromine (Br) input into the atmosphere of 495kmolBr (∼42tCHBr3) a-1. From Singapore ports the additional atmospheric Br input is calculated as 312kmolBr (∼26tCHBr3) a-1. We estimate the global anthropogenic Br input from ballast water into the atmosphere of up to 13Mmola-1. This is 0.1% global Br input from background bromoform emissions and thus probably not relevant for stratospheric ozone depletion.