ALBERT

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 (European Centre Hamburg general circulation model; ECHAM6.3) and an ocean model (Nucleus for European Modelling of the Ocean v3.6; NEMO3.6) as well as sea-ice (Louvain-la-Neuve sea Ice Model version 2; LIM2) and land surface model components (Jena Scheme for Biosphere Atmosphere Coupling in Hamburg; 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 (Model of Oceanic Pelagic Stoichiometry; 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 the Helmholtz Centre for Ocean Research Kiel (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, and 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 multidecadal eddy-rich configurations.

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: A set of experiments utilising different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-)tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Nino3 region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-)decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Ni\~no episodes. The behavior is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Description: The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

Description: A set of experiments utilizing different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-) tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Niño3-region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-) decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Niño episodes. The behaviour is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Description: The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

Description: We provide a time series of Agulhas leakage anomalies over the last twenty years from satellite altimetry. Until now, measuring the interannual variability of Indo-Atlantic exchange has been the major barrier in the investigation of the dynamics and large scale impact of Agulhas leakage. We compute the difference of transport between Agulhas Current and Agulhas Return Current which allows to deduce Agulhas leakage. The main difficulty is to separate the Agulhas Return Current from the southern limb of the subtropical "supergyre" south of Africa. For this purpose, an algorithm that uses absolute dynamic topography data is developed. The algorithm is applied to a state-of-the-art ocean model. The comparison with a Lagrangian method to measure the leakage allows to validate the new method. An important result is that it is possible to measure Agulhas leakage in this model using the velocity field along a section that crosses both the Agulhas Current and the Agulhas Return Current. In the model a good correlation is found between measuring leakage using the full depth velocities and using only the surface geostrophic velocities. This allows us to extend the method to along-track absolute dynamic topography from satellites. It is shown that the accuracy of the mean dynamic topography does not allow to determine the mean leakage but that leakage anomalies can be accurately computed.

Description: We provide a time series of Agulhas leakage anomalies over the last 20-years from satellite altimetry. Until now, measuring the interannual variability of Indo-Atlantic exchange has been the major barrier in the investigation of the dynamics and large scale impact of Agulhas leakage. We compute the difference of transport between the Agulhas Current and Agulhas Return Current, which allows us to deduce Agulhas leakage. The main difficulty is to separate the Agulhas Return Current from the southern limb of the subtropical "supergyre" south of Africa. For this purpose, an algorithm that uses absolute dynamic topography data is developed. The algorithm is applied to a state-of-the-art ocean model. The comparison with a Lagrangian method to measure the leakage allows us to validate the new method. An important result is that it is possible to measure Agulhas leakage in this model using the velocity field along a section that crosses both the Agulhas Current and the Agulhas Return Current. In the model a good correlation is found between measuring leakage using the full depth velocities and using only the surface geostrophic velocities. This allows us to extend the method to along-track absolute dynamic topography from satellites. It is shown that the accuracy of the mean dynamic topography does not allow us to determine the mean leakage but that leakage anomalies can be accurately computed.

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 analyzes 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. However, they also suggest a non-negligible DP contribution of at least 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 analyzes of decadal changes in the DP and AC contributions indicate that the ongoing increase in Agulhas leakage indeed may have evoked an increase in the AC contribution to the upper limb of the AMOC while the DP contribution decreased. In terms of thermohaline properties, our study highlights 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 those 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 ± 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.