Search

The upper wind-driven circulation in the tropical Atlantic Ocean plays a key role in the basin-wide distribution of water mass properties and affects the transport of heat, freshwater, and biogeochemical tracers such as oxygen or nutrients. It is crucial to improve our understanding of its long-term behaviour, which largely relies on model simulations and applied forcing due to sparse observational data coverage, especially before the mid-2000s. Here, we apply two different forcing products, the Coordinated Ocean-ice Reference Experiments (CORE) v2 and the Japanese 55-year Reanalysis (JRA55-do) surface dataset, to a high-resolution ocean model. Where possible, we compare the simulated results to long-term observations. We find large discrepancies between the two simulations regarding the wind and current field. In the CORE simulation, strong, large-scale wind stress curl amplitudes above the upwelling regions of the eastern tropical North Atlantic seem to cause an overestimation of the mean and seasonal variability in the eastward subsurface current just north of the Equator. The wind stress curl of JRA55-do forcing shows much finer structures, and the JRA55-do simulation is in better agreement with the mean and intraseasonal fluctuations in the subsurface current found in observations. The northern branch of the South Equatorial Current flows westward at the surface just north of the Equator. On interannual to decadal timescales, it shows a high correlation of R=0.9 with the zonal wind stress in the CORE simulation but only a weak correlation of R=0.35 in the JRA55-do simulation. We also identify similarities between the two simulations. The strength of the eastward-flowing North Equatorial Counter Current located between 3 and 10° N covaries with the strength of the meridional wind stress just north of the Equator on interannual to decadal timescales in the two simulations. Both simulations present a comparable mean, seasonal cycle and trend of the eastward off-equatorial subsurface current south of the Equator but underestimate the current strength by half compared to observations. In both simulations, the eastward-flowing Equatorial Undercurrent weakened between 1990 and 2009. In the JRA simulation, which covers the modern period of observations, the Equatorial Undercurrent strengthened again between 2008 to 2018, which agrees with observations, although the simulation underestimates the strengthening by over a third. We propose that long-term observations, once they have reached a critical length, need to be used to test the quality of wind-driven simulations. This study presents one step in this direction.

Abstract Changes in the Atlantic Meridional Overturning Circulation (AMOC) represent a crucial component of Northern Hemisphere climate variability. In modelling studies decadal overturning variability has been attributed to the intensity of deep winter convection in the Labrador Sea. This linkage is challenged by transport observations at sections across the subpolar gyre. Here we report simulations with an eddy-rich ocean model which captures the observed concentration of downwelling in the northeastern Atlantic and the negligible impact of interannual variations in Labrador Sea convection during the last decade. However, the exceptionally cold winters in the Labrador Sea during the first half of the 1990s induced a positive AMOC anomaly of more than 20%, mainly by augmenting the downwelling in the northeastern North Atlantic. The remote effect of excessive Labrador Sea buoyancy forcing is related to rapid spreading of mid-depth density anomalies into the Irminger Sea and their entrainment into the deep boundary current off Greenland.

Increasing Greenland Ice Sheet melting is anticipated to impact water mass transformation in the subpolar North Atlantic and ultimately the meridional overturning circulation. Complex ocean and climate models are widely applied to estimate magnitude and timing of related impacts under global warming. We discuss the role of the ocean mean state, subpolar water mass transformation, 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 and 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 make the Atlantic meridional overturning circulation less susceptible to the freshwater perturbation applied, and (3) a more even spreading and deeper mixing 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 feedback 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.

Capacity building and engagement are crucial to support science-led sustainable ocean policies that are integrated, interdisciplinary and internationally collaborative as well as being effective at a range of temporal and spatial scales.

Agulhas leakage increased in the 1960s through the 1980s but, despite the continued strengthening of the Southern Hemisphere westerlies, has not substantially increased in the 1990s and 2000s, according to observation- and model-based estimates.

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.

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.

Abstract The Makassar Strait, the main passageway of the Indonesian Throughflow (ITF), is an important component of Indo‐Pacific climate through its inter‐basin redistribution of heat and freshwater. Observational studies suggest that wind‐driven freshwater advection from the marginal seas into the Makassar Strait modulates the strait's surface transport. However, direct observations are too short (

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.

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.

In highly fragmented and relatively stable cold-seep ecosystems, species are expected to exhibit high migration rates and long-distance dispersal of long-lived pelagic larvae to maintain genetic integrity over their range. Accordingly, several species inhabiting cold seeps are widely distributed across the whole Atlantic Ocean, with low genetic divergence between metapopulations on both sides of the Atlantic Equatorial Belt (AEB, i.e. Barbados and African/European margins). Two hypotheses may explain such patterns: (i) the occurrence of present-day gene flow or (ii) incomplete lineage sorting due to large population sizes and low mutation rates. Here, we evaluated the first hypothesis using the cold seep mussels Gigantidas childressi, G. mauritanicus, Bathymodiolus heckerae and B. boomerang. We combined COI barcoding of 763 individuals with VIKING20X larval dispersal modelling at a large spatial scale not previously investigated. Population genetics supported the parallel evolution of Gigantidas and Bathymodiolus genera in the Atlantic Ocean and the occurrence of a 1-3 Million-year-old vicariance effect that isolated populations across the Caribbean Sea. Both population genetics and larval dispersal modelling suggested that contemporary gene flow and larval exchanges are possible across the AEB and the Caribbean Sea, although probably rare. When occurring, larval flow was eastward (AEB - only for B. boomerang) or northward (Caribbean Sea - only for G. mauritanicus). Caution is nevertheless required since we focused on only one mitochondrial gene, which may underestimate gene flow if a genetic barrier exists. Non-negligible genetic differentiation occurred between Barbados and African populations, so we could not discount the incomplete lineage sorting hypothesis. Larval dispersal modelling simulations supported the genetic findings along the American coast with high amounts of larval flow between the Gulf of Mexico (GoM) and the US Atlantic Margin, although the Blake Ridge population of B. heckerae appeared genetically differentiated. Overall, our results suggest that additional studies using nuclear genetic markers and population genomics approaches are needed to clarify the evolutionary history of the Atlantic bathymodioline mussels and to distinguish between ongoing and past processes.

Juvenile sea turtles can disperse thousands of kilometers from nesting beaches to oceanic development habitats, aided by ocean currents. In the North Atlantic, turtles dispersing from American beaches risk being advected out of warm nursery grounds in the North Atlantic Gyre into lethally cold northern European waters (e.g. around the UK). We used an ocean model simulation to compare simulated numbers of turtles that were advected to cold waters around the UK with observed numbers of turtles reported in the same area over ~5 decades. Rates of virtual turtles predicted to encounter lethal (10°C) or detrimental (15°C) temperatures (mean 19% ± 2.7 SD) and reach the UK were consistently low (median 0.83%, lower quartile 0.67%, upper quartile 1.02%), whereas there was high inter-annual variability in the numbers of dead or critically ill turtles reported in the UK. Generalized additive models suggest inter-annual variability in the North Atlantic Oscillation (NAO) index to be a good indicator of annual numbers of turtle strandings reported in the UK. We demonstrate that NAO variability drives variability in the dispersion scenarios of juvenile turtles from key nesting regions into the North Atlantic. Coastal effects, such as the number of storms and mean sea surface temperatures in the UK, were significant but weak predictors, with a weak effect on turtle strandings. Further understanding how changing environmental conditions such as NAO variability and storms affect the fate of juvenile turtles is vital for understanding the distribution and population dynamics of sea turtles.

Changes in stratospheric ozone concentrations and increasing concentrations of greenhouse gases (GHGs) alter the temperature structure of the atmosphere and drive changes in the atmospheric and oceanic circulation. We systematically investigate the impacts of ozone recovery and increasing GHGs on the atmospheric and oceanic circulation in the Southern Hemisphere during the twenty-first century using a unique coupled ocean–atmosphere climate model with interactive ozone chemistry and enhanced oceanic resolution. We use the high-emission scenario SSP5-8.5 for GHGs under which the springtime Antarctic total column ozone returns to 1980s levels by 2048 in our model, warming the lower stratosphere and strengthening the stratospheric westerly winds. We perform a spatial analysis and show for the first time that the austral spring stratospheric response to GHGs exhibits a marked planetary wavenumber 1 (PW1) pattern, which reinforces the response to ozone recovery over the Western Hemisphere and weakens it over the Eastern Hemisphere. These changes, which imply an eastward phase shift in the PW1, largely cancel out in the zonal mean. The Southern Hemisphere residual circulation strengthens during most of the year due to the increase in GHGs and weakens in spring due to ozone recovery. However, we find that in November the GHGs also drive a weakening of the residual circulation, reinforcing the effect of ozone recovery, which represents another novel result. At the surface, the westerly winds weaken and shift equatorward due to ozone recovery, driving a weak decrease in the transport of the Antarctic Circumpolar Current and in the Agulhas leakage and a cooling of the upper ocean, which is most pronounced in the latitudinal band 35–45∘ S. The increasing GHGs drive changes in the opposite direction that overwhelm the ozone effect. The total changes at the surface and in the oceanic circulation are nevertheless weaker in the presence of ozone recovery than those induced by GHGs alone, highlighting the importance of the Montreal Protocol in mitigating some of the impacts of climate change. We additionally compare the combined effect of interactively calculated ozone recovery and increasing GHGs with their combined effect in an ensemble in which we prescribe the CMIP6 ozone field. This second ensemble simulates a weaker ozone effect in all the examined fields, consistent with its weaker increase in ozone. The magnitude of the difference between the simulated changes at the surface and in the oceanic circulation in the two ensembles is as large as the ozone effect itself. This shows the large uncertainty that is associated with the choice of the ozone field and how the ozone is treated.

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 of the Atlantic Meridional Overturning Circulation (AMOC), and in particular the long-term temporal evolution, 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. This points to a dominant role of the wind forcing. The ability of the model to represent regional observations in western boundary current (WBC) systems at 53∘ N, 26.5∘ N and 11∘ S is explored. The question is investigated of whether WBC systems are able to represent the AMOC, and in particular whether these WBC systems exhibit similar temporal evolution to that of the zonally integrated 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, once 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 of representing realistic forcing-related and ocean-intrinsic trends.

Leakage of warm, salty waters from the Indian Ocean into the Atlantic increases by up to 40 % in high-resolution numerical ocean model simulations, suggesting that low-resolution models underestimate this key part of the global meridional overturning circulation.

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.

Landscape maps based on multivariate cluster analyses provide an objective and comprehensive view on the (marine) environment. They can hence support decision making regarding sustainable ocean resource handling and protection schemes. Across a large number of scales, input parameters and classification methods, numerous studies categorize the ocean into seascapes, hydro-morphological provinces or clusters. Many of them are regional, however, while only a few are on a basin scale. This study presents an automated cluster analysis of the entire Atlantic seafloor environment, based on eight global datasets and their derivatives: Bathymetry, slope, terrain ruggedness index, topographic position index, sediment thickness, POC flux, salinity, dissolved oxygen, temperature, current velocity, and phytoplankton abundance in surface waters along with seasonal variabilities. As a result, we obtained nine seabed areas (SBAs) that portray the Atlantic seafloor. Some SBAs have a clear geological and geomorphological nature, while others are defined by a mixture of terrain and water body characteristics. The majority of the SBAs, especially those covering the deep ocean areas, are coherent and show little seasonal and hydrographic variation, whereas other, nearshore SBAs, are smaller sized and dominated by high seasonal changes. To demonstrate the potential use of the marine landscape map for marine spatial planning purposes, we mapped out local SBA diversity using the patch richness index developed in landscape ecology. It identifies areas of high landscape diversity, and is a practical way of defining potential areas of interest, e.g. for designation as protected areas, or for further research. Clustering probabilities are highest (100%) in the center of SBA patches and decrease towards the edges (< 98%). On the SBA point cloud which was reduced for probabilities

Marine Heatwaves (MHWs) are ocean extreme events, characterized by anomalously high temperatures, which can have significant ecological impacts. The Northeast U.S. continental shelf is of great economical importance as it is home to a highly productive ecosystem. Local warming rates exceed the global average and the region experienced multiple MHWs in the last decade with severe consequences for regional fisheries. Due to the lack of subsurface observations, the depth-extent of MHWs is not well-known, which hampers the assessment of impacts on pelagic and benthic ecosystems. This study utilizes a global ocean circulation model with a high-resolution (1/20°) nest in the Atlantic to investigate the depth structure of MHWs and associated drivers on the Northeast U.S. continental shelf. It is shown that MHWs exhibit varying spatial extents, with some only occurring at depth. The highest intensities are found around 100 m depth with temperatures exceeding the climatological mean by up to 7°C, while surface intensities are typically smaller (around 3°C). Distinct vertical structures are associated with different spatial MHW patterns and drivers. Investigation of the co-variability of temperature and salinity reveals that over 80% of MHWs at depth (>50 m) coincide with extreme salinity anomalies. Two case studies provide insight into opposing MHW patterns at the surface and at depth, being forced by anomalous air-sea heat fluxes and Gulf Stream warm core ring interaction, respectively. The results highlight the importance of local ocean dynamics and the need to realistically represent them in climate models.

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.

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.

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.

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.