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.

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 North Atlantic Deep Water (NADW) is a crucial component of the Atlantic meridional overturning circulation and is therefore an important factor of the climate system. In order to estimate the mean relative contributions, sources, and pathways of the NADW at the southern exit of the Labrador Sea, a Lagrangian particle experiment is performed. The particles were seeded according to the strength of the velocity field along the 53∘ N section and traced 40 years backward in time in the three-dimensional velocity and hydrography field. The resulting transport pathways, their sources and corresponding transit timescales were inferred. Our experiment shows that, of the 30.1 Sv of NADW passing 53∘ N on average, the majority of this water is associated with a diapycnal mass flux without contact to the atmosphere, accounting for 14.3 Sv (48 %), where 6.2 Sv originate from the Labrador Sea, compared to 4.7 Sv from the Irminger Sea. The second-largest contribution originates from the mixed layer with 7.2 Sv (24 %), where the Labrador Sea contribution (5.9 Sv) dominates over the Irminger Sea contribution (1.0 Sv). Another 5.7 Sv (19 %) of NADW crosses the Greenland–Scotland Ridge within the NADW density class, where about two-thirds pass the Denmark Strait, while one-third crosses the Iceland–Scotland Ridge. The NADW exported at 53∘ N is hence dominated by entrainment through the diapycnal mass flux and mixed-layer origin in the Labrador Sea.

Description: With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.

Description: A hierarchy of global 1/4° (ORCA025) and Atlantic Ocean 1/20° nested (VIKING20X) ocean/sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation, and in particular the long-term temporal evolution, of the Atlantic Meridional Overturning Circulation (AMOC) strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5° N pointing at a dominant role of the forcing. Regional observations in western boundary current systems at 53° N, 26.5° N and 11° S are explored in respect to their ability to represent the AMOC and to monitor the temporal evolution of the AMOC. Apart from the basin-scale measurements at 26.5° N, it is shown that in particular the outflow of North Atlantic Deepwater at 53° N is a good indicator of the subpolar AMOC trend during the recent decades, if the latter is provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable in representing realistic forcing-related and ocean-intrinsic trends.

Description: Observations of the eastern subpolar North Atlantic in the 2010s show exceptional freshening and cooling of the upper ocean, peaking in 2016 with the lowest salinities recorded for 120 years. Published theories for the mechanisms driving the freshening include: reduced transport of saltier, warmer surface waters northwards from the subtropics associated with reduced meridional overturning; shifts in the pathways of fresher, cooler surface water from the Labrador Sea driven by changing patterns of wind stress; and the eastward expansion of the subpolar gyre. Using output from a high-resolution hindcast model simulation, we propose that the primary cause of the exceptional freshening and cooling is reduced surface heat loss in the Labrador Sea. Tracking virtual fluid particles in the model backwards from the eastern subpolar North Atlantic between 1990 and 2020 shows the major cause of the freshening and cooling to be an increased outflow of relatively fresh and cold surface waters from the Labrador Sea; with a minor contribution from reduced transport of warmer, saltier surface water northward from the subtropics. The cooling, but not the freshening, produced by these changing proportions of waters of subpolar and subtropical origin is mitigated by reduced along-track heat loss to the atmosphere in the North Atlantic Current. We analyse modelled boundary exchanges and water mass transformation in the Labrador Sea to show that since 2000, while inflows of lighter surface waters remain steady, the increasing output of these waters is due to reduced surface heat loss in the Labrador Sea beginning in the early 2000s. Tracking particles further upstream reveals that the primary source of the increased volume of lighter water transported out of the Labrador Sea is increased recirculation of water, and therefore longer residence times, in the upper 500–1000 m of the subpolar gyre.

Description: The inflow of relatively warm and salty water from the Indian Ocean into the South Atlantic via Agulhas leakage is important for the global overturning circulation and the global climate. In this study, we analyse the robustness of Agulhas leakage estimates as well as the thermohaline property modifications of Agulhas leakage south of Africa. Lagrangian experiments with both the newly developed tool Parcels and the well established tool Ariane were performed to simulate Agulhas leakage in the eddy-rich ocean–sea-ice model INALT20 (1/20∘ horizontal resolution) forced by the JRA55-do atmospheric boundary conditions. The average transport, its variability, trend and the transit time from the Agulhas Current to the Cape Basin of Agulhas leakage is simulated comparably with both Lagrangian tools, emphasizing the robustness of our method. Different designs of the Lagrangian experiment alter in particular the total transport of Agulhas leakage by up to 2 Sv, but the variability and trend of the transport are similar across these estimates. During the transit from the Agulhas Current at 32∘ S to the Cape Basin, a cooling and freshening of Agulhas leakage waters occurs especially at the location of the Agulhas Retroflection, resulting in a density increase as the thermal effect dominates. Beyond the strong air–sea exchange around South Africa, Agulhas leakage warms and salinifies the water masses below the thermocline in the South Atlantic.

Description: The Antarctic ozone hole has led to substantial changes in the Southern Hemisphere atmospheric circulation, such as the strengthening and poleward shift of the midlatitude westerly jet. Ozone recovery during the twenty-first century is expected to continue to affect the jet's strength and position, leading to changes in the opposite direction compared to the twentieth century and competing with the effect of increasing greenhouse gases. Simulations of the Earth's past and future climate, such as those performed for the Coupled Model Intercomparison Project Phase 6 (CMIP6), require an accurate representation of these ozone effects. Climate models that use prescribed ozone fields lack the important feedbacks between ozone chemistry, radiative heating, dynamics, and transport. In addition, when the prescribed ozone field was not generated by the same model to which it is prescribed, the imposed ozone hole is inconsistent with the simulated dynamics. These limitations ultimately affect the climate response to ozone depletion. This study investigates the impact of prescribing the ozone field recommended for CMIP6 on the simulated effects of ozone depletion in the Southern Hemisphere. We employ a new state-of-the-art coupled climate model, Flexible Ocean Climate Infrastructure (FOCI), to compare simulations in which the CMIP6 ozone is prescribed with simulations in which the ozone chemistry is calculated interactively. At the same time, we compare the roles played by ozone depletion and by increasing concentrations of greenhouse gases in driving changes in the Southern Hemisphere atmospheric circulation using a series of historical sensitivity simulations. FOCI captures the known effects of ozone depletion, simulating an austral spring and summer intensification of the midlatitude westerly winds and of the Brewer–Dobson circulation in the Southern Hemisphere. Ozone depletion is the primary driver of these historical circulation changes in FOCI. The austral spring cooling of the polar cap in the lower stratosphere in response to ozone depletion is weaker in the simulations that prescribe the CMIP6 ozone field. We attribute this weaker response to a prescribed ozone hole that is different to the model dynamics and is not collocated with the simulated polar vortex, altering the strength and position of the planetary wavenumber one. As a result, the dynamical contribution to the ozone-induced austral spring lower-stratospheric cooling is suppressed, leading to a weaker cooling trend. Consequently, the intensification of the polar night jet is also weaker in the simulations with prescribed CMIP6 ozone. In contrast, the differences in the tropospheric westerly jet response to ozone depletion fall within the internal variability present in the model. The persistence of the Southern Annular Mode is shorter in the prescribed ozone chemistry simulations. The results obtained with the FOCI model suggest that climate models that prescribe the CMIP6 ozone field still simulate a weaker Southern Hemisphere stratospheric response to ozone depletion compared to models that calculate the ozone chemistry interactively.

Description: Increasing Greenland Ice Sheet–melting is anticipated to impact watermass transformation in the subpolar North Atlantic and ultimately the meridional overturning circulation. Complex ocean and climate models are widely applied to predict magnitude and timing of related impacts under projected future climate. We discuss the role of the ocean mean state, subpolar gyre circulation, mesoscale eddies and atmospheric coupling in shaping the response of the subpolar North Atlantic Ocean to enhanced Greenland runoff. In a suite of eight dedicated 60 to 100-year long model experiments with and without atmospheric coupling, with eddy processes parameterized and explicitly simulated, with regular and significantly enlarged Greenland runoff, we find (1) a major impact by the interactive atmosphere in enabling a compensating temperature feedback, (2) a non-negligible influence by the ocean mean state biased towards greater stability in the coupled simulations, both of which making the Atlantic Merdional Overturning Circulation less susceptible to the freshwater perturbation applied, and (3) a more even spreading of the runoff tracer in the subpolar North Atlantic and enhanced inter-gyre exchange with the subtropics in the strongly eddying simulations. Overall, our experiments demonstrate the important role of mesoscale ocean dynamics and atmosphere feedbacks in projections of the climate system response to enhanced Greenland Ice Sheet–melting and hence underline the necessity to advance scale-aware eddy parameterizations for next-generation climate models.

Description: The Agulhas Current (AC) off the southern tip of Africa is one of the strongest western boundary currents and a crucial choke point of inter-ocean heat and salt exchange between the Indian Ocean and the southern Atlantic Ocean. However, large uncertainties remain concerning the sea surface temperature (SST) and salinity (SSS) variability in the AC region and their driving mechanisms over longer timescales, due to only short observational datasets being available and the highly dynamic nature of the region. Here, we present an annual coral skeletal Sr/Ca composite record paired with an established composite oxygen isotope record from Ifaty and Tulear reefs in southwestern Madagascar to obtain a 334-year (1661-1995) reconstruction of δ18Oseawater changes related to surface salinity variability in the wider Agulhas Current region. Our new annual δ18Oseawater composite record from Ifaty traces surface salinity of the southern Mozambique Channel and AC core region from the SODA reanalysis between 1958 and 1995. δ18Oseawater appears to be mainly driven by large-scale wind forcing in the southern Indian Ocean on interannual to decadal timescales. The δ18Oseawater and SST at Ifaty show characteristic interannual variability of between 2 and 4 years and interdecadal variability of 8 to 16 years, coherent with El Niño-Southern Oscillation (ENSO) records. Lagged correlations with the multivariate ENSO index reveals a 1-2-year lag of δ18Oseawater and salinity at Ifaty and the AC region, suggesting that propagation of anomalies by ocean Rossby waves may contribute to salinity changes in the wider southwestern Indian Ocean. The δ18Oseawater and SST reconstructions at Ifaty reveal the highest interannual variability during the Little Ice Age, especially around 1700 CE, which is in agreement with other Indo-Pacific coral studies. Our study demonstrates the huge potential to unlock past interannual and decadal changes in surface ocean hydrology and ocean transport dynamics from coral δ18Oseawater beyond the short instrumental record.