ALBERT

Description: The Lagrangian analysis of sets of particles advected with the flow fields of ocean models is used to study connectivity, that is, exchange pathways, time scales, 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. Because of spatial and temporal discretization, turbulence is not fully resolved in modeled velocities, and the concept of eddy diffusivity is used to parameterize the impact of unresolved processes. However, the relations between observation- 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 it 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 noneddying ocean model (ORCA05), Lagrangian eddy diffusivities are greatly reduced, particularly when the Gent and McWilliams parameterization of mesoscale eddies is employed.

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–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 300 m 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 part originate from the Agulhas, leading to higher SSTs in the tropics.

Description: The Agulhas Current plays a crucial role in the thermohaline circulation through its leakage into the South Atlantic Ocean. 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 intensity of westerlies 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 one or two decades. Contrary to previous assertions, these results further show that an equatorward (poleward) shift in westerlies increases (decreases) leakage. This occurs because of 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: The relationship between the Agulhas Current and the Agulhas leakage is not well understood. Here, this is investigated using two basin-scale and two global ocean models of incrementally increasing resolution. The response of the Agulhas Current is evaluated under a series of sensitivity experiments, in which idealized anomalies, designed to geometrically modulate zonal trade wind stress, are applied across the Indian Ocean Basin. The imposed wind stress changes exceed plus or minus two standard deviations from the annual-mean trade winds and, in the case of intensification, are partially representative of recently observed trends. The Agulhas leakage is quantified using complimentary techniques based on Lagrangian virtual floats and Eulerian passive tracer flux. As resolution increases, model behavior converges and the sensitivity of the leakage to Agulhas Current transport anomalies is reduced. In the two eddy-resolving configurations tested, the leakage is insensitive to changes in Agulhas Current transport at 32°S, though substantial eddy kinetic energy anomalies are evident. Consistent with observations, the position of the retroflection remains stable. The decoupling of Agulhas Current variability from the Agulhas leakage suggests that while correlations between the two may exist, they may not have a clear dynamical basis. It is suggested that present and future Agulhas leakage proxies should be considered in the context of potentially transient forcing regimes.

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 ± 560 nmol L-1 (226 ± 142 μg L-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 1000 pmol m-2 h-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 495 kmol Br (∼42 t CHBr3) a-1. From Singapore ports the additional atmospheric Br input is calculated as 312 kmol Br (∼26 t CHBr3) a-1. We estimate the global anthropogenic Br input from ballast water into the atmosphere of up to 13 Mmol a-1. This is 0.1 % global Br input from background bromoform emissions and thus probably not relevant for stratospheric ozone depletion.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence onto the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding nests covering the South Atlantic and the western Indian oceans at 1/10° (INALT10) and 1/20° (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host, but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60° resolution (INALT60) has been developed to gain insights into sub-mesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both, the southward flowing Agulhas Current and the northward flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

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 DBPs from oxidative ballast water treatment is bromoform (CHBr3), for which we find an average concentration of 894±560 nmol L−1 (226±142 µg L−1) in the undiluted ballast water from measurements and the 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 quantify the oceanic bromoform concentration and emissions 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 DBPs and anthropogenic bromoform concentrations in the sea surface. Based on our simulations, we expect DBPs to spread into the open ocean, along the coast and through 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 measurements. In a moderate scenario in which 70 % of the ballast water is chemically treated, bromoform emissions to the atmosphere can locally exceed 1000 pmol m−2 h−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 495 kmol Br a−1 (∼42 t CHBr3). For Singapore ports the additional atmospheric Br input is calculated as 312 kmol Br a−1 (∼ 26 t CHBr3). We estimate a global anthropogenic Br input from ballast water into the atmosphere of up to 13 Mmol a−1. This is 0.1 % of global Br input from background bromoform emissions and thus not relevant for stratospheric ozone depletion.

Description: Bromoform is the major by-product from chlorination of cooling water in coastal power plants. The number of power plants in East and Southeast Asian economies has increased rapidly, exceeding mean global growth. Bottom-up estimates of bromoform emissions based on few measurements appear to under-represent the industrial sources of bromoform from East Asia. Using oceanic Lagrangian analyses, we assess the amount of bromoform produced from power plant cooling-water treatment in East and Southeast Asia. The spread of bromoform is simulated as passive particles that are advected using the three-dimensional velocity fields over the years 2005/2006 from the high-resolution NEMO-ORCA0083 ocean general circulation model. Simulations are run for three scenarios with varying initial bromoform concentrations based on the range of bromoform measurements in cooling-water discharge. Comparing the modelled anthropogenic bromoform to in situ observations in the surface ocean and atmosphere, the two lower scenarios show the best agreement, suggesting initial bromoform concentrations in cooling water to be around 20–60 µg L−1. Based on these two scenarios, the model produces elevated bromoform in coastal waters of East Asia with average concentrations of 23 and 68 pmol L−1 and maximum values in the Yellow Sea, Sea of Japan and East China Sea. The industrially produced bromoform is quickly emitted into the atmosphere with average air–sea flux of 3.1 and 9.1 nmolm-2h-1, respectively. Atmospheric abundances of anthropogenic bromoform are derived from simulations with the Lagrangian particle dispersion model FLEXPART based on ERA-Interim wind fields in 2016. In the marine boundary layer of East Asia, the FLEXPART simulations show mean anthropogenic bromoform mixing ratios of 0.4–1.3 ppt, which are 2–6 times larger compared to the climatological bromoform estimate. During boreal winter, the simulations show that some part of the anthropogenic bromoform is transported by the northeasterly winter monsoon towards the tropical regions, whereas during boreal summer anthropogenic bromoform is confined to the Northern Hemisphere subtropics. Convective events in the tropics entrain an additional 0.04–0.05 ppt of anthropogenic bromoform into the stratosphere, averaged over tropical Southeast Asia. In our simulations, only about 10 % of anthropogenic bromoform is outgassed from power plants located in the tropics south of 20∘ N, so that only a small fraction of the anthropogenic bromoform reaches the stratosphere. We conclude that bromoform from cooling-water treatment in East Asia is a significant source of atmospheric bromine and might be responsible for annual emissions of 100–300 Mmol of Br in this region. These anthropogenic bromoform sources from industrial water treatment might be a missing factor in global flux estimates of organic bromine. While the current emissions of industrial bromoform provide a significant contribution to regional tropospheric budgets, they provide only a minor contribution to the stratospheric bromine budget of 0.24–0.30 ppt of Br.

Description: Bottom pressure observations on both sides of the Atlantic basin, combined with satellite measurements of sea level anomalies and wind stress data, are utilized to estimate variations of the Atlantic Meridional Overturning Circulation (AMOC) at 11∘ S. Over the period 2013–2018, the AMOC and its components are dominated by seasonal variability, with peak-to-peak amplitudes of 12 Sv for the upper-ocean geostrophic transport, 7 Sv for the Ekman and 14 Sv for the AMOC transport. The characteristics of the observed seasonal cycles of the AMOC and its components are compared to results from an ocean general circulation model, which is known to reproduce the variability of the Western Boundary Current on longer timescales. The observed seasonal variability of zonally integrated geostrophic velocity in the upper 300 m is controlled by pressure variations at the eastern boundary, while at 500 m depth contributions from the western and eastern boundaries are similar. The model tends to underestimate the seasonal pressure variability at 300 and 500 m depth, especially at the western boundary, which translates into the estimate of the upper-ocean geostrophic transport. In the model, seasonal AMOC variability at 11∘ S is governed, besides the Ekman transport, by the geostrophic transport variability in the eastern basin. The geostrophic contribution of the western basin to the seasonal cycle of the AMOC is instead comparably weak, as transport variability in the western basin interior related to local wind curl forcing is mainly compensated by the Western Boundary Current. Our analyses indicate that while some of the uncertainties of our estimates result from the technical aspects of the observational strategy or processes not being properly represented in the model, uncertainties in the wind forcing are particularly relevant for the resulting uncertainties of AMOC estimates at 11∘ S.