ALBERT

Description: Dimethyl sulfide (DMS), dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) were measured at the Boknis Eck Time Series Station (BE, Eckernförde Bay, SW Baltic Sea) during the period February 2009–December 2018. Our results show considerable interannual and seasonal variabilities in the mixed-layer concentrations of DMS, total DMSP (DMSPt) and total DMSO (DMSOt). Positive correlations were found between particulate DMSP (DMSPp) and particulate DMSO (DMSOp) as well as DMSPt and DMSOt in the mixed layer, suggesting a similar source for both compounds. The decreasing long-term trends, observed for DMSPt and DMS in the mixed layer, were linked to the concurrent trend of the sum of 19′-hexanoyloxyfucoxanthin and 19′-butanoyloxy-fucoxanthin, which are the marker pigments of prymnesiophytes and chrysophytes, respectively. Major Baltic inflow (MBI) events influenced the distribution of sulfur compounds due to phytoplankton community changes, and sediment might be a potential source for DMS in the bottom layer during seasonal hypoxia/anoxia at BE. A modified algorithm based on the phytoplankton pigments reproduces the DMSPp : Chl a ratios well during this study and could be used to estimate future surface (5 m) DMSPp concentrations at BE.

Description: Coastal areas contribute significantly to the emissions of methane (CH4) from the ocean. In order to decipher its temporal variability in the whole water column, dissolved CH4 was measured on a monthly basis at the Boknis Eck Time-series Station (BE) located in the Eckernförde Bay (SW Baltic Sea) from 2006 to 2017. BE has a water depth of about 28 m and dissolved CH4 was measured at six water depths ranging from 0 to 25 m. In general CH4 concentrations increased with depth, indicating a sedimentary release of CH4. Pronounced enhancement of the CH4 concentrations in the bottom layer (15–25 m) was found during February, May–June and October. CH4 was not correlated with Chlorophyll a or O2 over the measurement period. Unusually high CH4 concentrations (of up to 696 nM) were sporadically observed in the upper layer (0–10 m) (e.g. in November 2013 and December 2014) and were coinciding with Major Baltic Inflow (MBI) events. Surface CH4 concentrations were always supersaturated throughout the monitoring period, indicating that the Eckernförde Bay is an intense but highly variable source of atmospheric CH4. We did not detect significant temporal trends in CH4 concentrations or emissions, despite of ongoing environmental changes such as warming and deoxygenation in the Eckernförde Bay. Overall, the CH4 variability at BE is driven by a complex interplay of various biological and physical processes.

Description: The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated, simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993-2012). The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences due to emissions and transport processes.Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA’s long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements - including high altitude observations from the NASA Global Hawk platform. The models generally capture the seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model-measurement correlation (r ≥ 0.7) and a low sensitivity to the choice of emission inventory, at most sites. In a given model, the absolute model-measurement agreement is highly sensitive to the choice of emissions and inter-model differences are also apparent, even when using the same inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve optimal agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2 Br2). In general, the models are able to reproduce well observations of CHBr3 and CH2 Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2 Br2) most elevated over the tropical West Pacific during boreal winter. The models also indicate the Asian Monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models. We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2 Br2 of 2.0 (1.2-2.5) ppt, ≫ 57% larger than the best estimate from the most recent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. However, transport-driven inter-annual variability in the annual mean bromine SGI is of the order of a ±5%, with SGI exhibiting a strong positive correlation with ENSO in the East Pacific.

Description: Nitrous oxide (N2O) is a potent greenhouse gas, and it is involved in stratospheric ozone depletion. Its oceanic production is mainly influenced by dissolved nutrient and oxygen (O2) concentrations in the water column. Here we examined the seasonal and annual variations in dissolved N2O at the Boknis Eck (BE) Time Series Station located in Eckernförde Bay (southwestern Baltic Sea). Monthly measurements of N2O started in July 2005. We found a pronounced seasonal pattern for N2O with high concentrations (supersaturations) in winter and early spring and low concentrations (undersaturations) in autumn when hypoxic or anoxic conditions prevail. Unusually low N2O concentrations were observed during October 2016–April 2017, which was presumably a result of prolonged anoxia and the subsequent nutrient deficiency. Unusually high N2O concentrations were found in November 2017 and this event was linked to the occurrence of upwelling which interrupted N2O consumption via denitrification and potentially promoted ammonium oxidation (nitrification) at the oxic–anoxic interface. Nutrient concentrations (such as nitrate, nitrite and phosphate) at BE have been decreasing since the 1980s, but oxygen concentrations in the water column are still decreasing. Our results indicate a close coupling of N2O anomalies to O2 concentration, nutrients, and stratification. Given the long-term trends of declining nutrient and oxygen concentrations at BE, a decrease in N2O concentration, and thus emissions, seems likely due to an increasing number of events with low N2O concentrations.

Description: Ice-nucleating particles (INPs) have a large impact on the climate-relevant properties of clouds over the oceans. Studies have shown that sea spray aerosols (SSAs), produced upon bursting of bubbles at the ocean surface, can be an important source of marine INPs, particularly during periods of enhanced biological productivity. Recent mesocosm experiments using natural seawater spiked with nutrients have revealed that marine INPs are derived from two separate classes of organic matter in SSAs. Despite this finding, existing parameterizations for marine INP abundance are based solely on single variables such as SSA organic carbon (OC) or SSA surface area, which may mask specific trends in the separate classes of INP. The goal of this paper is to improve the understanding of the connection between ocean biology and marine INP abundance by reporting results from a field study and proposing a new parameterization of marine INPs that accounts for the two associated classes of organic matter. The PEACETIME cruise took place from 10 May to 10 June 2017 in the Mediterranean Sea. Throughout the cruise, INP concentrations in the surface microlayer (INPSML) and in SSAs (INPSSA) produced using a plunging aquarium apparatus were continuously monitored while surface seawater (SSW) and SML biological properties were measured in parallel. The organic content of artificially generated SSAs was also evaluated. INPSML concentrations were found to be lower than those reported in the literature, presumably due to the oligotrophic nature of the Mediterranean Sea. A dust wet deposition event that occurred during the cruise increased the INP concentrations measured in the SML by an order of magnitude, in line with increases in iron in the SML and bacterial abundances. Increases in INPSSA were not observed until after a delay of 3 days compared to increases in the SML and are likely a result of a strong influence of bulk SSW INPs for the temperatures investigated (T=−18 ∘C for SSAs, T=−15 ∘C for SSW). Results confirmed that INPSSA are divided into two classes depending on their associated organic matter. Here we find that warm (T≥−22 ∘C) INPSSA concentrations are correlated with water-soluble organic matter (WSOC) in the SSAs, but also with SSW parameters (particulate organic carbon, POCSSW and INPSSW,−16C) while cold INPSSA (T〈−22 ∘C) are correlated with SSA water-insoluble organic carbon (WIOC) and SML dissolved organic carbon (DOC) concentrations. A relationship was also found between cold INPSSA and SSW nano- and microphytoplankton cell abundances, indicating that these species might be a source of water-insoluble organic matter with surfactant properties and specific IN activities. Guided by these results, we formulated and tested multiple parameterizations for the abundance of INPs in marine SSAs, including a single-component model based on POCSSW and a two-component model based on SSA WIOC and OC. We also altered a previous model based on OCSSA content to account for oligotrophy of the Mediterranean Sea. We then compared this formulation with the previous models. This new parameterization should improve attempts to incorporate marine INP emissions into numerical models.

Description: Nitric oxide (NO) is a short-lived intermediate of the oceanic nitrogen cycle, however, due to its high reactivity, measurements of dissolved NO in seawater are rare. Here we present an improved method to determine NO concentrations in discrete seawater samples. The set-up of our system consisted of a chemiluminescence NO analyser connected to a stripping unit. The limit of detection for our method was 5 pmol NO in aqueous solution which translates into 0.25 nmol L−1 when using a 20 mL seawater sample volume. Our method was applied to measure high resolution depth profiles of dissolved NO during a cruise to the eastern tropical South Pacific Ocean. Our method is fast and comparably easy to handle thus it opens the door for deciphering the distribution of NO in the ocean and it facilitates laboratory studies on NO pathways.

Description: A compilation of the published literature on dust content in terrestrial and marine sediment cores was synchronized with pollen data and speleothem growth phases on the Greenland Ice Core Chronology 2005 (GICC05) time axis. Aridity patterns for eight key areas of the global climate system have been reconstructed for the last 60 000 years. These records have different time resolutions and different dating methods, i.e. different types of stratigraphy. Nevertheless, all regions analysed in this study show humid conditions during early Marine Isotope Stage 3 (MIS3) and the early Holocene or deglaciation, but not always at the same time. Such discrepancies have been interpreted as regional effects, although stratigraphic uncertainties may affect some of the proposed interpretations. In comparison, most of the MIS2 interval becomes arid in all of the Northern Hemisphere records, but the peak arid conditions of the Last Glacial Maximum (LGM) and Heinrich event 1 differ in duration and intensity among regions. In addition, we also compare the aridity synthesis with modelling results using a global climate model (GCM). Indeed, geological archives and GCMs show agreement on the aridity pattern for the Holocene or deglaciation, for the LGM and for late MIS3.

Description: We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.

Description: To better understand possible reasons for the diverse modeling results and large discrepancies of the detected solar fingerprints, we took one step back and assessed the "initial" solar signals in the middle atmosphere based on large ensemble simulations with multiple climate models — FOCI, EMAC, and MPI-ESM-HR. Consistent with previous work, we find that the 11-year solar cycle signals in the short wave heating rate (SWHR) and ozone anomalies are robust and statistically significant in all three models. These "initial" solar cycle signals in SWHR, ozone, and temperature anomalies are sensitive to the strength of the solar forcing. Correlation coefficients of the solar cycle with the SWHR, ozone, and temperature anomalies linearly increase along with the enhancement of the solar cycle amplitude, and this reliance becomes more complex when the solar cycle amplitude exceeds a certain threshold. In addition, the cold bias in the tropical stratopause of EMAC dampens the subsequent results of the "initial" solar signal. The warm pole bias in MPI-ESM-HR leads to a weak polar night jet (PNJ), which may limit the top-down propagation of the initial solar signal. Although FOCI simulated a so-called top-down response as revealed in previous studies in a period with large solar cycle amplitudes, its warm bias in the tropical upper stratosphere results in a positive bias in PNJ and can lead to a "reversed" response in some extreme cases. We suggest a careful interpretation of the single model result and further re-examination of the solar signal based on more climate models.

Description: Reanalysis data sets are widely used to understand atmospheric processes and past variability, and are often used to stand in as "observations" for comparisons with climate model output. Because of the central role of water vapor (WV) and ozone (O3) in climate change, it is important to understand how accurately and consistently these species are represented in existing global reanalyses. In this paper, we present the results of WV and O3 intercomparisons that have been performed as part of the SPARC (Stratosphere–troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The comparisons cover a range of timescales and evaluate both inter-reanalysis and observation-reanalysis differences. We also provide a systematic documentation of the treatment of WV and O3 in current reanalyses to aid future research and guide the interpretation of differences amongst reanalysis fields. The assimilation of total column ozone (TCO) observations in newer reanalyses results in realistic representations of TCO in reanalyses except when data coverage is lacking, such as during polar night. The vertical distribution of ozone is also relatively well represented in the stratosphere in reanalyses, particularly given the relatively weak constraints on ozone vertical structure provided by most assimilated observations and the simplistic representations of ozone photochemical processes in most of the reanalysis forecast models. However, significant biases in the vertical distribution of ozone are found in the upper troposphere and lower stratosphere in all reanalyses. In contrast to O3, reanalysis estimates of stratospheric WV are not directly constrained by assimilated data. Observations of atmospheric humidity are typically used only in the troposphere, below a specified vertical level at or near the tropopause. The fidelity of reanalysis stratospheric WV products is therefore mainly dependent on the reanalyses' representation of the physical drivers that influence stratospheric WV, such as temperatures in the tropical tropopause layer, methane oxidation, and the stratospheric overturning circulation. The lack of assimilated observations and known deficiencies in the representation of stratospheric transport in reanalyses result in much poorer agreement amongst observational and reanalysis estimates of stratospheric WV. Hence, stratospheric WV products from the current generation of reanalyses should generally not be used in scientific studies.

Description: The Indian Ocean is coupled to atmospheric dynamics and chemical composition via several unique mechanisms, such as the seasonally varying monsoon circulation. During the winter monsoon season, high pollution levels are regularly observed over the entire northern Indian Ocean, while during the summer monsoon, clean air dominates the atmospheric composition, leading to distinct chemical regimes. The changing atmospheric composition over the Indian Ocean can interact with oceanic biogeochemical cycles and impact marine ecosystems, resulting in potential climate feedbacks. Here, we review current progress in detecting and understanding atmospheric gas-phase composition over the Indian Ocean and its local and global impacts. The review considers results from recent Indian Ocean ship campaigns, satellite measurements, station data, and information on continental and oceanic trace gas emissions. The distribution of all major pollutants and greenhouse gases shows pronounced differences between the landmass source regions and the Indian Ocean, with strong gradients over the coastal areas. Surface pollution and ozone are highest during the winter monsoon over the Bay of Bengal and the Arabian Sea coastal waters due to air mass advection from the Indo-Gangetic Plain and continental outflow from Southeast Asia. We observe, however, that unusual types of wind patterns can lead to pronounced deviations of the typical trace gas distributions. For example, the ozone distribution maxima shift to different regions under wind scenarios that differ from the regular seasonal transport patterns. The distribution of greenhouse gases over the Indian Ocean shows many similarities when compared to the pollution fields, but also some differences of the latitudinal and seasonal variations resulting from their long lifetimes and biogenic sources. Mixing ratios of greenhouse gases such as methane show positive trends over the Indian Ocean, but long-term changes in pollution and ozone due to changing emissions and transport patterns require further investigation. Although we know that changing atmospheric composition and perturbations within the Indian Ocean affect each other, the impacts of atmospheric pollution on oceanic biogeochemistry and trace gas cycling are severely understudied. We highlight potential mechanisms, future research topics, and observational requirements that need to be explored in order to fully understand such interactions and feedbacks in the Indian Ocean region.

Description: We present consistent annual mean atmospheric histories and growth rates for the mainly anthropogenic halogenated compounds HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14 and PFC-116, which are all potentially useful oceanic transient tracers (tracers of water transport within the ocean), for the Northern and Southern Hemisphere with the aim of providing input histories of these compounds for the equilibrium between the atmosphere and surface ocean. We use observations of these halogenated compounds made by the Advanced Global Atmospheric Gases Experiment (AGAGE), the Scripps Institution of Oceanography (SIO), the Commonwealth Scientific and Industrial Research Organization (CSIRO), the National Oceanic and Atmospheric Administration (NOAA) and the University of East Anglia (UEA). Prior to the direct observational record, we use archived air measurements, firn air measurements and published model calculations to estimate the atmospheric mole fraction histories. The results show that the atmospheric mole fractions for each species, except HCFC-141b and HCFC-142b, have been increasing since they were initially produced. Recently, the atmospheric growth rates have been decreasing for the HCFCs (HCFC-22, HCFC-141b and HCFC-142b), increasing for the HFCs (HFC-134a, HFC-125, HFC-23) and stable with little fluctuation for the PFCs (PFC-14 and PFC-116) investigated here. The atmospheric histories (source functions) and natural background mole fractions show that HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125 and HFC-23 have the potential to be oceanic transient tracers for the next few decades only because of the recently imposed bans on production and consumption. When the atmospheric histories of the compounds are not monotonically changing, the equilibrium atmospheric mole fraction (and ultimately the age associated with that mole fraction) calculated from their concentration in the ocean is not unique, reducing their potential as transient tracers. Moreover, HFCs have potential to be oceanic transient tracers for a longer period in the future than HCFCs as the growth rates of HFCs are increasing and those of HCFCs are decreasing in the background atmosphere. PFC-14 and PFC-116, however, have the potential to be tracers for longer periods into the future due to their extremely long lifetimes, steady atmospheric growth rates and no explicit ban on their emissions. In this work, we also derive solubility functions for HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14 and PFC-116 in water and seawater to facilitate their use as oceanic transient tracers. These functions are based on the Clark–Glew–Weiss (CGW) water solubility function fit and salting-out coefficients estimated by the poly-parameter linear free-energy relationships (pp-LFERs). Here we also provide three methods of seawater solubility estimation for more compounds. Even though our intention is for application in oceanic research, the work described in this paper is potentially useful for tracer studies in a wide range of natural waters, including freshwater and saline lakes, and, for the more stable compounds, groundwaters.

Description: Halocarbons are produced naturally in the oceans by biological and chemical processes. They are emitted from surface seawater into the atmosphere, where they take part in numerous chemical processes such as ozone destruction and the oxidation of mercury and dimethyl sulfide. Here we present oceanic and atmospheric halocarbon data for the Peruvian upwelling zone obtained during the M91 cruise onboard the research vessel METEOR in December 2012. Surface waters during the cruise were characterized by moderate concentrations of bromoform (CHBr3) and dibromomethane (CH2Br2) correlating with diatom biomass derived from marker pigment concentrations, which suggests this phytoplankton group is a likely source. Concentrations measured for the iodinated compounds methyl iodide (CH3I) of up to 35.4 pmol L−1, chloroiodomethane (CH2ClI) of up to 58.1 pmol L−1 and diiodomethane (CH2I2) of up to 32.4 pmol L−1 in water samples were much higher than previously reported for the tropical Atlantic upwelling systems. Iodocarbons also correlated with the diatom biomass and even more significantly with dissolved organic matter (DOM) components measured in the surface water. Our results suggest a biological source of these compounds as a significant driving factor for the observed large iodocarbon concentrations. Elevated atmospheric mixing ratios of CH3I (up to 3.2 ppt), CH2ClI (up to 2.5 ppt) and CH2I2 (3.3 ppt) above the upwelling were correlated with seawater concentrations and high sea-to-air fluxes. During the first part of the cruise, the enhanced iodocarbon production in the Peruvian upwelling contributed significantly to tropospheric iodine levels, while this contribution was considerably smaller during the second part.

Description: Very short-lived substances (VSLS) contribute as source gases significantly to the tropospheric and stratospheric bromine loading. At present, an estimated 25% of stratospheric bromine is of oceanic origin. In this study, we investigate how climate change may impact the ocean- atmosphere flux of brominated VSLS, their atmospheric transport, and chemical transformations and evaluate how these changes will affect stratospheric ozone over the 21st century. Under the assumption of fixed ocean water concentrations and RCP6.0 scenario, we find an increase of the ocean- atmosphere flux of brominated VSLS of about 8-10% by the end of the 21st century compared to present day. A decrease in the tropospheric mixing ratios of VSLS and an increase in the lower stratosphere are attributed to changes in atmospheric chemistry and transport. Our model simulations reveal that this increase is counteracted by a corresponding reduction of inorganic bromine. Therefore the total amount of bromine from VSLS in the stratosphere will not be changed by an increase in upwelling. Part of the increase of VSLS in the tropical lower stratosphere results from an increase in the corresponding tropopause height. As the depletion of stratospheric ozone due to bromine depends also on the availability of chlorine, we find the impact of bromine on stratospheric ozone at the end of the 21st century reduced compared to present day. Thus, these studies highlight the different factors influencing the role of brominated VSLS in a future climate

Description: We use isoprene and related field measurements from three different ocean data sets together with remotely sensed satellite data to model global marine isoprene emissions. We show that using monthly mean satellite-derived chl a concentrations to parameterize isoprene with a constant chl a normalized isoprene production rate underpredicts the measured oceanic isoprene concentration by a mean factor of 19 ± 12. Improving the model by using phytoplankton functional type dependent production values and by decreasing the bacterial degradation rate of isoprene in the water column results in only a slight underestimation (factor 1.7 ± 1.2). We calculate global isoprene emissions of 0.21 Tg C for 2014 using this improved model, which is twice the value calculated using the original model. Nonetheless, the sea-to-air fluxes have to be at least 1 order of magnitude higher to account for measured atmospheric isoprene mixing ratios. These findings suggest that there is at least one missing oceanic source of isoprene and, possibly, other unknown factors in the ocean or atmosphere influencing the atmospheric values. The discrepancy between calculated fluxes and atmospheric observations must be reconciled in order to fully understand the importance of marine-derived isoprene as a precursor to remote marine boundary layer particle formation.

Description: In this study, we investigate the maximum physical and biogeochemical potential of macroalgae open-ocean mariculture and sinking (MOS) as an ocean-based carbon dioxide removal (CDR) method. Embedding a macroalgae model into an Earth system model, we simulate macroalgae mariculture in the open-ocean surface layer followed by fast sinking of the carbon-rich macroalgal biomass to the deep seafloor (depth〉3000 m), which assumes no remineralization of the harvested biomass during the quick sinking. We also test the combination of MOS with artificial upwelling (AU), which fertilizes the macroalgae by pumping nutrient-rich deeper water to the surface. The simulations are done under RCP 4.5, a moderate-emissions pathway. When deployed globally between years 2020 and 2100, the carbon captured and exported by MOS is 270 PgC, which is further boosted by AU of 447 PgC. Because of feedbacks in the Earth system, the oceanic carbon inventory only increases by 171.8 PgC (283.9 PgC with AU) in the idealized simulations. More than half of this carbon remains in the ocean after cessation at year 2100 until year 3000. The major side effect of MOS on pelagic ecosystems is the reduction of phytoplankton net primary production (PNPP) due to the competition for nutrients with macroalgae and due to canopy shading. MOS shrinks the mid-layer oxygen-minimum zones (OMZs) by reducing the organic matter export to, and remineralization in, subsurface and intermediate waters, while it creates new OMZs on the seafloor by oxygen consumption from remineralization of sunken biomass. MOS also impacts the global carbon cycle by reducing the atmospheric and terrestrial carbon reservoirs when enhancing the ocean carbon reservoir. MOS also enriches dissolved inorganic carbon in the deep ocean. Effects are mostly reversible after cessation of MOS, though recovery is not complete by year 3000. In a sensitivity experiment without remineralization of sunken MOS biomass, the whole of the MOS-captured carbon is permanently stored in the ocean, but the lack of remineralized nutrients causes a long-term nutrient decline in the surface layers and thus reduces PNPP. Our results suggest that MOS has, theoretically, considerable CDR potential as an ocean-based CDR method. However, our simulations also suggest that such large-scale deployment of MOS would have substantial side effects on marine ecosystems and biogeochemistry, up to a reorganization of food webs over large parts of the ocean.

Description: Based on the Baltic Earth Assessment Reports of this thematic issue in Earth System Dynamics and recent peer-reviewed literature, current knowledge of the effects of global warming on past and future changes in climate of the Baltic Sea region is summarised and assessed. The study is an update of the Second Assessment of Climate Change (BACC II) published in 2015 and focuses on the atmosphere, land, cryosphere, ocean, sediments, and the terrestrial and marine biosphere. Based on the summaries of the recent knowledge gained in palaeo-, historical, and future regional climate research, we find that the main conclusions from earlier assessments still remain valid. However, new long-term, homogenous observational records, for example, for Scandinavian glacier inventories, sea-level-driven saltwater inflows, so-called Major Baltic Inflows, and phytoplankton species distribution, and new scenario simulations with improved models, for example, for glaciers, lake ice, and marine food web, have become available. In many cases, uncertainties can now be better estimated than before because more models were included in the ensembles, especially for the Baltic Sea. With the help of coupled models, feedbacks between several components of the Earth system have been studied, and multiple driver studies were performed, e.g. projections of the food web that include fisheries, eutrophication, and climate change. New datasets and projections have led to a revised understanding of changes in some variables such as salinity. Furthermore, it has become evident that natural variability, in particular for the ocean on multidecadal timescales, is greater than previously estimated, challenging our ability to detect observed and projected changes in climate. In this context, the first palaeoclimate simulations regionalised for the Baltic Sea region are instructive. Hence, estimated uncertainties for the projections of many variables increased. In addition to the well-known influence of the North Atlantic Oscillation, it was found that also other low-frequency modes of internal variability, such as the Atlantic Multidecadal Variability, have profound effects on the climate of the Baltic Sea region. Challenges were also identified, such as the systematic discrepancy between future cloudiness trends in global and regional models and the difficulty of confidently attributing large observed changes in marine ecosystems to climate change. Finally, we compare our results with other coastal sea assessments, such as the North Sea Region Climate Change Assessment (NOSCCA), and find that the effects of climate change on the Baltic Sea differ from those on the North Sea, since Baltic Sea oceanography and ecosystems are very different from other coastal seas such as the North Sea. While the North Sea dynamics are dominated by tides, the Baltic Sea is characterised by brackish water, a perennial vertical stratification in the southern subbasins, and a seasonal sea ice cover in the northern subbasins.

Description: This paper presents a modelling study on the fate of CHBr3 and its product gases in the troposphere within the context of tropical deep convection. A cloud-scale case study was conducted along the west coast of Borneo, where several deep convective systems were triggered on the afternoon and early evening of 19 November 2011. These systems were sampled by the Falcon aircraft during the field campaign of the SHIVA project and analysed using a simulation with the cloud-resolving meteorological model C-CATT-BRAMS at 2x2 km resolution that represents the emissions, transport by large-scale flow, convection, photochemistry, and washout of CHBr3 and its product gases (PGs). We find that simulated CHBr3 mixing ratios and the observed values in the boundary layer and the outflow of the convective systems agree. However, the model underestimates the background CHBr3 mixing ratios in the upper troposphere, which suggests a missing source at the regional scale. An analysis of the simulated chemical speciation of bromine within and around each simulated convective system during the mature convective stage reveals that 〉 85% of the bromine derived from CHBr3 and its PGs is transported vertically to the point of convective detrainment in the form of CHBr3 and that the remaining small fraction is in the form of organic PGs, principally insoluble brominated carbonyls produced from the photo-oxidation of CHBr3. The model simulates that within the boundary layer and free troposphere, the inorganic PGs are only present in soluble forms, i.e. HBr, HOBr, and BrONO2, and, consequently, within the convective clouds, the inorganic PGs are almost entirely removed by wet scavenging. We find that HBr is the most abundant PG in background lower-tropospheric air and that this prevalence of HBr is a result of the relatively low background tropospheric ozone levels at the regional scale. Contrary to a previous study in a different environment, for the conditions in the simulation, the insoluble Br-2 species is hardly formed within the convective systems and therefore plays no significant role in the vertical transport of bromine. This likely results from the relatively small quantities of simulated inorganic bromine involved, the presence of HBr in large excess compared to HOBr and BrO, and the relatively efficient removal of soluble compounds within the convective column.

Description: Particle sinking is a major form of transport for photosynthetically fixed carbon to below the euphotic zone via the biological carbon pump (BCP). Oxygen (O2) depletion may improve the efficiency of the BCP. However, the mechanisms by which O2 deficiency can enhance particulate organic matter (POM) vertical fluxes are not well understood. Here, we investigate the composition and vertical fluxes of POM in two deep basins of the Baltic Sea (GB: Gotland Basin and LD: Landsort Deep). The two basins showed different O2 regimes resulting from the intrusion of oxygen-rich water from the North Sea that ventilated the water column below 140 m in GB, but not in LD, during the time of sampling. In June 2015, we deployed surface-tethered drifting sediment traps in oxic surface waters (GB: 40 and 60 m; LD: 40 and 55 m), within the oxygen minimum zone (OMZ; GB: 110 m and LD: 110 and 180 m) and at recently oxygenated waters by the North Sea inflow in GB (180 m). The primary objective of this study was to test the hypothesis that the different O2 conditions in the water column of GB and LD affected the composition and vertical flux of sinking particles and caused differences in export efficiency between those two basins. The composition and vertical flux of sinking particles were different in GB and LD. In GB, particulate organic carbon (POC) flux was 18 % lower in the shallowest trap (40 m) than in the deepest sediment trap (at 180 m). Particulate nitrogen (PN) and Coomassie stainable particle (CSP) fluxes decreased with depth, while particulate organic phosphorus (POP), biogenic silicate (BSi), chlorophyll a (Chl a) and transparent exopolymeric particle (TEP) fluxes peaked within the core of the OMZ (110 m); this coincided with the presence of manganese oxide-like (MnOx-like) particles aggregated with organic matter. In LD, vertical fluxes of POC, PN and CSPs decreased by 28 %, 42 % and 56 %, respectively, from the surface to deep waters. POP, BSi and TEP fluxes did not decrease continuously with depth, but they were higher at 110 m. Although we observe a higher vertical flux of POP, BSi and TEPs coinciding with abundant MnOx-like particles at 110 m in both basins, the peak in the vertical flux of POM and MnOx-like particles was much higher in GB than in LD. Sinking particles were remarkably enriched in BSi, indicating that diatoms were preferentially included in sinking aggregates and/or there was an inclusion of lithogenic Si (scavenged into sinking particles) in our analysis. During this study, the POC transfer efficiency (POC flux at 180 m over 40 m) was higher in GB (115 %) than in LD (69 %), suggesting that under anoxic conditions a smaller portion of the POC exported below the euphotic zone was transferred to 180 m than under reoxygenated conditions present in GB. In addition, the vertical fluxes of MnOx-like particles were 2 orders of magnitude higher in GB than LD. Our results suggest that POM aggregates with MnOx-like particles formed after the inflow of oxygen-rich water into GB, and the formation of those MnOx–OM-rich particles may alter the composition and vertical flux of POM, potentially contributing to a higher transfer efficiency of POC in GB. This idea is consistent with observations of fresher and less degraded organic matter in deep waters of GB than LD.

Description: The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

Description: The Baltic Sea has a salinity gradient decreasing from fully marine (〉 25) in the west to below 7 in the central Baltic Proper. Habitat-forming and ecologically dominant mytilid mussels exhibit decreasing growth when salinity 〈 11; however, the mechanisms underlying reduced calcification rates in dilute seawater are not fully understood. Both [HCO−3] and [Ca2+] also decrease with salinity, challenging calcifying organisms through CaCO3 undersaturation (Ω≤1) and unfavourable ratios of calcification substrates ([Ca2+] and [HCO−3]) to the inhibitor (H+), expressed as the extended substrate–inhibitor ratio (ESIR). This study combined in situ monitoring of three southwest Baltic mussel reefs with two laboratory experiments to assess how various environmental conditions and isolated abiotic factors (salinity, [Ca2+], [HCO−3] and pH) impact calcification in mytilid mussels along the Baltic salinity gradient. Laboratory experiments rearing juvenile Baltic Mytilus at a range of salinities (6, 11 and 16), HCO−3 concentrations (300–2100 µmol kg−1) and Ca2+ concentrations (0.5–4 mmol kg−1) reveal that as individual factors, low [HCO−3], pH and salinity cannot explain low calcification rates in the Baltic Sea. Calcification rates are impeded when Ωaragonite ≤ 1 or ESIR ≤ 0.7 primarily due to [Ca2+] limitation which becomes relevant at a salinity of ca. 11 in the Baltic Sea. Field monitoring of carbonate chemistry and calcification rates suggest increased food availability may be able to mask the negative impacts of periodic sub-optimal carbonate chemistry, but not when seawater conditions are permanently adverse, as observed in two Baltic reefs at salinities 〈 11. Regional climate models predict a rapid desalination of the southwest and central Baltic over the next century and potentially a reduction in [Ca2+] which may shift the distribution of marine calcifiers westward. It is therefore vital to understand the mechanisms by which the ionic composition of seawater impacts bivalve calcification for better predicting the future of benthic Baltic ecosystems.

Description: Halogenated very short-lived substances (VSLSs) are naturally produced in the ocean and emitted to the atmosphere. When transported to the stratosphere, these compounds can have a significant influence on the ozone layer and climate. During a research cruise on RV Sonne in the subtropical and tropical west Indian Ocean in July and August 2014, we measured the VSLSs, methyl iodide (CH3I) and for the first time bromoform (CHBr3) and dibromomethane (CH2Br2), in surface seawater and the marine atmosphere to derive their emission strengths. Using the Lagrangian particle dispersion model FLEXPART with ERA-Interim meteorological fields, we calculated the direct contribution of observed VSLS emissions to the stratospheric halogen burden during the Asian summer monsoon. Furthermore, we compare the in situ calculations with the interannual variability of transport from a larger area of the west Indian Ocean surface to the stratosphere for July 2000–2015. We found that the west Indian Ocean is a strong source for CHBr3 (910 pmol m−2 h−1), very strong source for CH2Br2 (930 pmol m−2 h−1), and an average source for CH3I (460 pmol m−2 h−1). The atmospheric transport from the tropical west Indian Ocean surface to the stratosphere experiences two main pathways. On very short timescales, especially relevant for the shortest-lived compound CH3I (3.5 days lifetime), convection above the Indian Ocean lifts oceanic air masses and VSLSs towards the tropopause. On a longer timescale, the Asian summer monsoon circulation transports oceanic VSLSs towards India and the Bay of Bengal, where they are lifted with the monsoon convection and reach stratospheric levels in the southeastern part of the Asian monsoon anticyclone. This transport pathway is more important for the longer-lived brominated compounds (17 and 150 days lifetime for CHBr3 and CH2Br2). The entrainment of CHBr3 and CH3I from the west Indian Ocean to the stratosphere during the Asian summer monsoon is lower than from previous cruises in the tropical west Pacific Ocean during boreal autumn and early winter but higher than from the tropical Atlantic during boreal summer. In contrast, the projected CH2Br2 entrainment was very high because of the high emissions during the west Indian Ocean cruise. The 16-year July time series shows highest interannual variability for the shortest-lived CH3I and lowest for the longest-lived CH2Br2. During this time period, a small increase in VSLS entrainment from the west Indian Ocean through the Asian monsoon to the stratosphere is found. Overall, this study confirms that the subtropical and tropical west Indian Ocean is an important source region of halogenated VSLSs, especially CH2Br2, to the troposphere and stratosphere during the Asian summer monsoon.

Description: Marine particulate organic carbon-13 stable isotope ratios (δ13CPOC) provide insights in understanding carbon cycling through the atmosphere, ocean, and biosphere. They have been used to trace the input of anthropogenic carbon in the marine ecosystem due to the distinct isotopically light signature of anthropogenic emissions. However, δ13CPOC is also significantly altered during photosynthesis by phytoplankton, which complicates its interpretation. For such purposes, robust spatio-temporal coverage of δ13CP OC observations is essential. We collected all such available data sets, merged and homogenized them to provide the largest available marine δ13CPOC data set (Verwega et al., 2021). The data set consists of 4732 data points covering all major ocean basins beginning in the 1960s. We describe the compiled raw data, compare different observational methods, and provide key insights in the temporal and spatial distribution that is consistent with previously observed patterns. The main different sample collection methods (bottle, intake, net, trap) are generally consistent with each other when comparing within regions. An analysis of 1990s mean δ13CP OC values in an meridional section accross the Atlantic Ocean shows relatively high values (≥ −22 ‰) in the low latitudes (〈 30°) trending towards lower values in the Arctic Ocean (∼ −24 ‰) and Southern Ocean (≤ −28 ‰). The temporal trend since the 1960s shows a decrease of mean δ13CPOC by more than 3 ‰ in all basins except for the Southern Ocean which shows a weaker trend but contains relatively poor multi-decadal coverage.

Description: The eastern tropical South Pacific (ETSP) hosts the Peruvian upwelling system, which represents one of the most productive areas in the world ocean. High primary production followed by rapid heterotrophic utilization of organic matter supports the formation of one of the most intense oxygen minimum zones (OMZs) in the world ocean, where dissolved oxygen (O2) concentrations reach less than 1 µmol kg−1. The high productivity leads to an accumulation of dissolved organic matter (DOM) in the surface layers that may serve as a substrate for heterotrophic respiration. However, the importance of DOM utilization for O2 respiration in the Peruvian upwelling system in general and for shaping the upper oxycline in particular remains unclear so far. This study reports the first estimates of diapycnal fluxes and supply of O2, dissolved organic carbon (DOC), dissolved organic nitrogen, dissolved hydrolysable amino acids (DHAA) and dissolved combined carbohydrates (DCCHO) for the ETSP off Peru. Diapycnal flux and supply estimates were obtained by combining measured vertical diffusivities and solute concentration gradients. They were analysed together with the molecular composition of DCCHO and DHAA to infer the transport of labile DOM into the upper OMZ and the potential role of DOM utilization for the attenuation of the diapycnal O2 flux that ventilates the OMZ. The observed diapycnal O2 flux (50 mmol O2 m−2 d−1 at maximum) was limited to the upper 80 m of the water column; the O2 supply of ∼1 µmol kg−1 d−1 was comparable to previously published O2 consumption rates for the North and South Pacific OMZs. The diapycnal DOM flux (31 mmol C m−2 d−1 at maximum) was limited to ∼30 m water depth, suggesting that the labile DOM is extensively consumed within the upper part of the shallow oxycline off Peru. The analyses of DCCHO and DHAA composition support this finding, suggesting that DOM undergoes comprehensive remineralization within the upper part of the oxycline, as the DOM within the core of the OMZ was found to be largely altered. Estimated by a simple equation for carbon combustion, aerobic respiration of DCCHO and DHAA, supplied by diapycnal mixing (0.46 µmol kg−1 d−1 at maximum), could account for up to 38 % of the diapycnal O2 supply in the upper oxycline, which suggests that DOM utilization plays a significant role for shaping the upper oxycline in the ETSP.

Description: Our understanding of the biogeochemical cycling of the climate-relevant trace gas dimethyl sulfide (DMS) in the Peruvian upwelling system is still limited. Here we present oceanic and atmospheric DMS measurements which were made during two shipborne cruises in December 2012 (M91) and October 2015 (SO243) in the Peruvian upwelling region. Dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) were also measured during M91. DMS concentrations were 1.9 ± 0.9 and 2.5 ± 1.9 nmol L−1 in surface waters in October 2015 and December 2012, respectively. Nutrient availability appeared to be the main driver of the observed variability in the surface DMS distributions in the coastal areas. DMS, DMSP, and DMSO showed maxima in the surface layer, and no elevated concentrations associated with the oxygen minimum zone off Peru were measured. The possible role of DMS, DMSP, and DMSO as radical scavengers (stimulated by nitrogen limitation) is supported by their negative correlations with N:P (sum of nitrate and nitrite : dissolved phosphate) ratios. Large variations in atmospheric DMS mole fractions were measured during M91 (144.6 ± 95.0 ppt) and SO243 (91.4 ± 55.8 ppt); however, the atmospheric mole fractions were generally low, and the sea-to-air flux was primarily driven by seawater DMS. The Peruvian upwelling region was identified as a source of atmospheric DMS in December 2012 and October 2015. However, in comparison to the previous measurements in the adjacent regions, the Peru upwelling was a moderate source of DMS emissions at either time (M91: 5.9 ± 5.3 µmol m−2 d−1; SO243: 3.8 ± 2.7 µmol m−2 d−1).

Description: The climate active trace-gas carbonyl sulfide (OCS) is the most abundant sulfur gas in the atmosphere. A missing source in its atmospheric budget is currently suggested, resulting from an upward revision of the vegetation sink. Tropical oceanic emissions have been proposed to close the resulting gap in the atmospheric budget. We present a bottom-up approach including (i) new observations of OCS in surface waters of the tropical Atlantic, Pacific and Indian oceans and (ii) a further improved global box model to show that direct OCS emissions are unlikely to account for the missing source. The box model suggests an undersaturation of the surface water with respect to OCS integrated over the entire tropical ocean area and, further, global annual direct emissions of OCS well below that suggested by top-down estimates. In addition, we discuss the potential of indirect emission from CS2 and dimethylsulfide (DMS) to account for the gap in the atmospheric budget. This bottom-up estimate of oceanic emissions has implications for using OCS as a proxy for global terrestrial CO2 uptake, which is currently impeded by the inadequate quantification of atmospheric OCS sources and sinks.

Description: Circulation anomalies accompanying Sudden Stratospheric Warmings (SSWs) can have a significant impact on the troposphere. This surface response is observed for some but not all SSWs, and their downward coupling is not fully understood. We use an existing classification method to separate downward- and non-propagating SSWs (d/nSSWs) in ERA5 reanalysis data for the years 1979–2019. The differences in SSW downward propagation in composites of spatial patterns clearly show that dSSWs dominate the surface regional impacts following SSWs. During dSSWs, the upper-tropospheric jet stream is significantly displaced equatorward. Wave activity analysis shows remarkable differences between d/nSSWs for planetary and synoptic-scale waves. Enhanced stratospheric planetary eddy kinetic energy (EKE) and heat fluxes around the central date of dSSWs are followed by increased synoptic-scale wave activity and even surface coupling for synoptic-scale EKE. An observed significant reduction in upper-tropospheric synoptic-scale momentum fluxes following dSSWs confirms the important role of tropospheric eddy feedbacks for coupling to the surface. Our findings emphasize the role of the lower stratosphere and synoptic-scale waves in coupling the SSW signal to the surface and agree with mechanisms suggested in earlier modeling studies.

Description: Stratospheric sulfate aerosols from volcanic eruptions have a significant impact on the Earth's climate. To include the effects of volcanic eruptions in climate model simulations, the Easy Volcanic Aerosol (EVA) forcing generator provides stratospheric aerosol optical properties as a function of time, latitude, height, and wavelength for a given input list of volcanic eruption attributes. EVA is based on a parameterized three-box model of stratospheric transport and simple scaling relationships used to derive mid-visible (550 nm) aerosol optical depth and aerosol effective radius from stratospheric sulfate mass. Precalculated look-up tables computed from Mie theory are used to produce wavelength-dependent aerosol extinction, single scattering albedo, and scattering asymmetry factor values. The structural form of EVA and the tuning of its parameters are chosen to produce best agreement with the satellite-based reconstruction of stratospheric aerosol properties following the 1991 Pinatubo eruption, and with prior millennial-timescale forcing reconstructions, including the 1815 eruption of Tambora. EVA can be used to produce volcanic forcing for climate models which is based on recent observations and physical understanding but internally self-consistent over any timescale of choice. In addition, EVA is constructed so as to allow for easy modification of different aspects of aerosol properties, in order to be used in model experiments to help advance understanding of what aspects of the volcanic aerosol are important for the climate system.

Description: Coastal seas may account for more than 75 % of global oceanic methane emissions. There, methane is mainly produced microbially in anoxic sediments from where it can escape to the overlying water column. Aerobic methane oxidation (MOx) in the water column acts as a biological filter reducing the amount of methane that eventually evades to the atmosphere. The efficiency of the MOx filter is potentially controlled by the availability of dissolved methane and oxygen, as well as temperature, salinity, and hydrographic dynamics, and all of these factors undergo strong temporal fluctuations in coastal ecosystems. In order to elucidate the key environmental controls, specifically the effect of oxygen availability, on MOx in a seasonally stratified and hypoxic coastal marine setting, we conducted a 2-year time-series study with measurements of MOx and physico-chemical water column parameters in a coastal inlet in the southwestern Baltic Sea (Eckernförde Bay). We found that MOx rates always increased toward the seafloor, but were not directly linked to methane concentrations. MOx exhibited a strong seasonal variability, with maximum rates (up to 11.6 nmol l−1 d−1) during summer stratification when oxygen concentrations were lowest and bottom-water temperatures were highest. Under these conditions, 70–95 % of the sediment-released methane was oxidized, whereas only 40–60 % were consumed during the mixed and oxygenated periods. Laboratory experiments with manipulated oxygen concentrations in the range of 0.2–220 µmol l−1 revealed a sub-micromolar oxygen-optimum for MOx at the study site. In contrast, the fraction of methane-carbon incorporation into the bacterial biomass (compared to the total amount of oxidised methane) was up to 38-fold higher at saturated oxygen concentrations, suggesting a different partitioning of catabolic and anabolic processes under oxygen-replete and oxygen-starved conditions, respectively. Our results underscore the importance of MOx in mitigating methane emission from coastal waters and indicate an organism-level adaptation of the water column methanotrophs to hypoxic conditions.

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: Large amounts of methane are delivered by fluids through the erosive forearc of the convergent margin offshore Costa Rica and lead to the formation of cold seeps at the sediment surface. Besides mud extrusion, numerous cold seeps are created by landslides induced by seamount subduction or fluid migration along major faults. Most of the dissolved methane reaching the seafloor at cold seeps is oxidized within the benthic microbial methane filter by anaerobic oxidation of methane (AOM). Measurements of AOM and sulfate reduction as well as numerical modeling of porewater profiles revealed a highly active and efficient benthic methane filter at Quepos Slide site; a landslide on the continental slope between the Nicoya and Osa Peninsula. Integrated areal rates of AOM ranged from 12.9 ± 6.0 to 45.2 ± 11.5 mmol m-2 d-1, with only 1 to 2.5% of the upward methane flux being released into the water column. Additionally, two parallel sediment cores from Quepos Slide were used for in vitro experiments in a recently developed Sediment-F low-Through (SLOT) system to simulate an increased fluid and methane flux from the bottom of the sediment core. The benthic methane filter revealed a high adaptability whereby the methane oxidation efficiency responded to the increased fluid flow within 150–170 days. To our knowledge, this study provides the first estimation of the natural biogeochemical response of seep sediments to changes in fluid flow.

Description: The El Niño Southern Oscillation (ENSO) with its warm (El Niño) and cold (La Niña) phase has strong impacts on marine ecosystems off Peru. This influence extends from changes in nutrient availability to productivity and oxygen levels. While several studies have demonstrated the influence of ENSO events on biological productivity, less is known about their impact on oxygen concentrations. In situ observations along the Peruvian and Chilean coast have shown a strong water column oxygenation during the 1997/1998 strong El Niño event. These observations suggest a deepening of the oxygen minimum zone (OMZ) along the continental shelf. However, due to reduced spatial coverage of the existing in situ observations, no studies have yet demonstrated the OMZ response to El Niño events in the whole Eastern Tropical South Pacific (ETSP). Furthermore, most studies have focused on El Niño events. Much less attention was given to the oxygen dynamics under La Niña influence. Here, we provide a comprehensive analysis of the ENSO influence on OMZ dynamics. Interannual variability of the OMZ during the period 1990–2010 is derived from a regional coupled physical-biogeochemical model forced with realistic atmospheric and lateral boundary conditions. Our results show a reduction of the vertical extent and a deepening of suboxic waters (SW) during the El Niño phase. During the La Niña phase, there is a vertical expansion of SW. These fluctuations in OMZ extent are due to changes in oxygen supply into its core depth mainly from lateral margins. During the El Niño phase, the enhanced lateral oxygen supply from the subtropics is the main reason for the reduction of SW in both coastal and offshore regions. During the La Niña phase, the oxygenated subtropical waters are blocked by the poleward transport along the southern margin of the OMZ. Consequently, oxygen concentrations within the OMZ are reduced and suboxic conditions expand during La Niña. The detailed analysis of transport pathways presented here provides new insights into how ENSO variability affects the oxygen-sensitive marine biogeochemistry of the ETSP.

Description: The science guiding the EUREC4A campaign and its measurements is presented. EUREC4A comprised roughly 5 weeks of measurements in the downstream winter trades of the North Atlantic – eastward and southeastward of Barbados. Through its ability to characterize processes operating across a wide range of scales, EUREC4A marked a turning point in our ability to observationally study factors influencing clouds in the trades, how they will respond to warming, and their link to other components of the earth system, such as upper-ocean processes or the life cycle of particulate matter. This characterization was made possible by thousands (2500) of sondes distributed to measure circulations on meso- (200 km) and larger (500 km) scales, roughly 400 h of flight time by four heavily instrumented research aircraft; four global-class research vessels; an advanced ground-based cloud observatory; scores of autonomous observing platforms operating in the upper ocean (nearly 10 000 profiles), lower atmosphere (continuous profiling), and along the air–sea interface; a network of water stable isotopologue measurements; targeted tasking of satellite remote sensing; and modeling with a new generation of weather and climate models. In addition to providing an outline of the novel measurements and their composition into a unified and coordinated campaign, the six distinct scientific facets that EUREC4A explored – from North Brazil Current rings to turbulence-induced clustering of cloud droplets and its influence on warm-rain formation – are presented along with an overview of EUREC4A's outreach activities, environmental impact, and guidelines for scientific practice. Track data for all platforms are standardized and accessible at https://doi.org/10.25326/165 (Stevens, 2021), and a film documenting the campaign is provided as a video supplement.

Description: Chloromethane (CH3Cl) is the most important natural input of reactive chlorine to the stratosphere, contributing about 16 % to stratospheric ozone depletion. Due to the phase-out of anthropogenic emissions of chlorofluorocarbons, CH3Cl will largely control future levels of stratospheric chlorine. The tropical rainforest is commonly assumed to be the strongest single CH3Cl source, contributing over half of the global annual emissions of about 4000 to 5000 Gg (1 Gg = 109 g). This source shows a characteristic carbon isotope fingerprint, making isotopic investigations a promising tool for improving its atmospheric budget. Applying carbon isotopes to better constrain the atmospheric budget of CH3Cl requires sound information on the kinetic isotope effects for the main sink processes: the reaction with OH and Cl in the troposphere. We conducted photochemical CH3Cl degradation experiments in a 3500 dm3 smog chamber to determine the carbon isotope effect (ε=k13C/k12C−1 ) for the reaction of CH3Cl with OH and Cl. For the reaction of CH3Cl with OH, we determined an ε value of (−11.2±0.8) ‰ (n=3) and for the reaction with Cl we found an ε value of (−10.2±0.5 ) ‰ (n=1), which is 5 to 6 times smaller than previously reported. Our smaller isotope effects are strongly supported by the lack of any significant seasonal covariation in previously reported tropospheric δ13C(CH3Cl) values with the OH-driven seasonal cycle in tropospheric mixing ratios. Applying these new values for the carbon isotope effect to the global CH3Cl budget using a simple two hemispheric box model, we derive a tropical rainforest CH3Cl source of (670±200) Gg a−1, which is considerably smaller than previous estimates. A revision of previous bottom-up estimates, using above-ground biomass instead of rainforest area, strongly supports this lower estimate. Finally, our results suggest a large unknown CH3Cl source of (1530±200) Gg a−1.

Description: Depth profiles of nitrous oxide (N2O) were measured during six cruises to the upwelling area and oxygen minimum zone (OMZ) off Peru in 2009 and 2012/2013, covering both the coastal shelf region and the adjacent open ocean. N2O profiles displayed a strong sensitivity towards oxygen concentrations. Open ocean profiles with distances to the shelf break larger than the first baroclinic Rossby radius of deformation showed a transition from a broad maximum close to the Equator to a double-peak structure south of 5° S where the oxygen minimum was more pronounced. Maximum N2O concentrations in the open ocean were about 80 nM. A linear relationship between ΔN2O and apparent oxygen utilization (AOU) could be found for measurements within the upper oxycline, with a slope similar to studies in other oceanic regions. In contrast, N2O profiles close to the shelf revealed a much higher variability, and N2O concentrations higher than 100 nM were often observed. The highest N2O concentration measured at the shelf was  ∼  850 nM. Due to the extremely sharp oxygen gradients at the shelf, N2O maxima occurred in very shallow water depths of less than 50 m. In the coastal area, a linear relationship between ΔN2O and AOU could not be observed as extremely high ΔN2O values were scattered over the full range of oxygen concentrations. The data points that showed the strongest deviation from a linear ΔN2O ∕ AOU relationship also showed signals of intense nitrogen loss. These results indicate that the coastal upwelling at the Peruvian coast and the subsequent strong remineralization in the water column causes conditions that lead to extreme N2O accumulation, most likely due to the interplay of intense mixing and high rates of remineralization which lead to a rapid switching of the OMZ waters between anoxic and oxic conditions. This, in turn, could trigger incomplete denitrification or pulses of increased nitrification with extreme N2O production.

Description: The consideration of marine biogeochemistry is essential for simulating the carbon cycle in an Earth system model. Here we present the implementation and evaluation of a marine biogeochemical model, Model of Oceanic Pelagic Stoichiometry (MOPS) in the Flexible Ocean and Climate Infrastructure (FOCI) climate model. FOCI-MOPS enables the simulation of marine biological processes, the marine carbon, nitrogen and oxygen cycles, air-sea gas exchange of CO2 and O2, and simulations with prescribed atmospheric CO2 or CO2 emissions. A series of experiments covering the historical period (1850–2014) were performed following the DECK (Diagnostic, Evaluation and Characterization of Klima) and CMIP6 (Coupled Model Intercomparison Project 6) protocols. Overall, modelled biogeochemical tracer distributions and fluxes, as well as transient evolution in surface air temperature, air-sea CO2 fluxes, and changes of ocean carbon and heat, are in good agreement with observations. Modelled inorganic and organic tracer distributions are quantitatively evaluated by statistically-derived metrics. Results of the FOCI-MOPS model, also including sea surface temperature, surface pH, oxygen (100–600 m), nitrate (0–100 m), and primary production, are within the range of other CMIP6 model results. Overall, the evaluation of FOCI-MOPS indicates its suitability for Earth climate system simulations

Description: The coastal upwelling regime off Peru in December 2012 showed considerable vertical concentration gradients of dissolved nitrous oxide (N2O) across the top few meters of the ocean. The gradients were predominantly downward, i.e., concentrations decreased toward the surface. Ignoring these gradients causes a systematic error in regionally integrated gas exchange estimates, when using observed concentrations at several meters below the surface as input for bulk flux parameterizations – as is routinely practiced. Here we propose that multi-day near-surface stratification events are responsible for the observed near-surface N2O gradients, and that the gradients induce the strongest bias in gas exchange estimates at winds of about 3 to 6 m s−1. Glider hydrographic time series reveal that events of multi-day near-surface stratification are a common feature in the study region. In the same way as shorter events of near-surface stratification (e.g., the diurnal warm layer cycle), they preferentially exist under calm to moderate wind conditions, suppress turbulent mixing, and thus lead to isolation of the top layer from the waters below (surface trapping). Our observational data in combination with a simple gas-transfer model of the surface trapping mechanism show that multi-day near-surface stratification can produce near-surface N2O gradients comparable to observations. They further indicate that N2O gradients created by diurnal or shorter stratification cycles are weaker and do not substantially impact bulk emission estimates. Quantitatively, we estimate that the integrated bias for the entire Peruvian upwelling region in December 2012 represents an overestimation of the total N2O emission by about a third, if concentrations at 5 or 10 m depth are used as surrogate for bulk water N2O concentration. Locally, gradients exist which would lead to emission rates overestimated by a factor of two or more. As the Peruvian upwelling region is an N2O source of global importance, and other strong N2O source regions could tend to develop multi-day near-surface stratification as well, the bias resulting from multi-day near-surface stratification may also impact global oceanic N2O emission estimates.

Description: The ice–substrate interface is an important boundary condition for ice sheet modelling. The substrate affects the ice sheet by allowing sliding through sediment deformation and accommodating the storage and drainage of subglacial water. We present three datasets on a 1 : 5 000 000 scale with different geological parameters for the region that was covered by the ice sheets in North America, including Greenland and Iceland. The first dataset includes the distribution of surficial sediments, which is separated into continuous, discontinuous and predominantly rock categories. The second dataset includes sediment grain size properties, which is divided into three classes: clay, silt and sand, based on the dominant grain size of the fine fraction of the glacial sediments. The third dataset is the generalized bedrock geology. We demonstrate the utility of these datasets for governing ice sheet dynamics by using an ice sheet model with a simulation that extends through the last glacial cycle. In order to demonstrate the importance of the basal boundary conditions for ice sheet modelling, we changed the shear friction angle to account for a weaker substrate and found changes up to 40 % in ice thickness compared to a reference run. Although incorporation of the ice–bed boundary remains model dependent, our dataset provides an observational baseline for improving a critical weakness in current ice sheet modelling (https://doi.org/10.1594/PANGAEA.895889, Gowan et al., 2018b).

Description: Simulations of the glacial–interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a “cookbook” for the growing community of PISM users and paleo-ice sheet modelers in general. For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.

Description: The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation–age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model–data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016). This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4±4.1 m (or 6.5±2.0×106km3 ), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.

Description: Ocean alkalinity enhancement (OAE) is a marine carbon dioxide removal (CDR) approach. Publicly funded research projects have begun, and philanthropic funding and start-ups are collectively pushing the field forward. This rapid progress in research activities has created an urgent need to learn if and how OAE can work at scale. The Best Practices Guide to OAE research contains 7 topics broken down into 13 chapters that compare and synthesise previously published methods and offer guidance for future research.

Description: Halocarbons from oceanic sources contribute to halogens in the troposphere, and can be transported into the stratosphere where they take part in ozone depletion. This paper presents distribution and sources in the equatorial Atlantic from June and July 2011 of the four compounds bromoform (CHBr3), dibromomethane (CH2Br2), methyl iodide (CH3I) and diiodomethane (CH2I2). Enhanced biological production during the Atlantic Cold Tongue (ACT) season, indicated by phytoplankton pigment concentrations, led to elevated concentrations of CHBr3 of up to 44.7 and up to 9.2 pmol L−1 for CH2Br2 in surface water, which is comparable to other tropical upwelling systems. While both compounds correlated very well with each other in the surface water, CH2Br2 was often more elevated in greater depth than CHBr3, which showed maxima in the vicinity of the deep chlorophyll maximum. The deeper maximum of CH2Br2 indicates an additional source in comparison to CHBr3 or a slower degradation of CH2Br2. Concentrations of CH3I of up to 12.8 pmol L−1 in the surface water were measured. In contrary to expectations of a predominantly photochemical source in the tropical ocean, its distribution was mostly in agreement with biological parameters, indicating a biological source. CH2I2 was very low in the near surface water with maximum concentrations of only 3.7 pmol L−1. CH2I2 showed distinct maxima in deeper waters similar to CH2Br2. For the first time, diapycnal fluxes of the four halocarbons from the upper thermocline into and out of the mixed layer were determined. These fluxes were low in comparison to the halocarbon sea-to-air fluxes. This indicates that despite the observed maximum concentrations at depth, production in the surface mixed layer is the main oceanic source for all four compounds and one of the main driving factors of their emissions into the atmosphere in the ACT-region. The calculated production rates of the compounds in the mixed layer are 34 ± 65 pmol m−3 h−1 for CHBr3, 10 ± 12 pmol m−3 h−1 for CH2Br2, 21 ± 24 pmol m−3 h−1 for CH3I and 384 ± 318 pmol m−3 h−1 for CH2I2 determined from 13 depth profiles.

Description: We present a transient simulation of global vegetation and climate patterns of the mid- and late Holocene using the MPI-ESM (Max Planck Institute for Meteorology Earth System Model) at T63 resolution. The simulated vegetation trend is discussed in the context of the simulated Holocene climate change. Our model captures the main trends found in reconstructions. Most prominent are the southward retreat of the northern treeline that is combined with the strong decrease of forest in the high northern latitudes during the Holocene and the vast increase of the Saharan desert, embedded in a general decrease in precipitation and vegetation in the Northern Hemisphere monsoon margin regions. The Southern Hemisphere experiences weaker changes in total vegetation cover during the last 8000 years. However, the monsoon-related increase in precipitation and the insolation-induced cooling of the winter climate lead to shifts in the vegetation composition, mainly between the woody plant functional types (PFTs). The large-scale global patterns of vegetation almost linearly follow the subtle, approximately linear, orbital forcing. In some regions, however, non-linear, more rapid changes in vegetation are found in the simulation. The most striking region is the Sahel–Sahara domain with rapid vegetation transitions to a rather desertic state, despite a gradual insolation forcing. Rapid shifts in the simulated vegetation also occur in the high northern latitudes, in South Asia and in the monsoon margins of the Southern Hemisphere. These rapid changes are mainly triggered by changes in the winter temperatures, which go into, or move out of, the bioclimatic tolerance range of individual PFTs. The dynamics of the transitions are determined by dynamics of the net primary production (NPP) and the competition between PFTs. These changes mainly occur on timescales of centuries. More rapid changes in PFTs that occur within a few decades are mainly associated with the timescales of mortality and the bioclimatic thresholds implicit in the dynamic vegetation model, which have to be interpreted with caution. Most of the simulated Holocene vegetation changes outside the high northern latitudes are associated with modifications in the intensity of the global summer monsoon dynamics that also affect the circulation in the extra tropics via teleconnections. Based on our simulations, we thus identify the global monsoons as the key player in Holocene climate and vegetation change.

Description: It is well established that variable wintertime planetary wave forcing in the stratosphere controls the variability of Arctic stratospheric ozone through changes in the strength of the polar vortex and the residual circulation. While previous studies focused on the variations in upward wave flux entering the lower stratosphere, here the impact of downward planetary wave reflection on ozone is investigated for the first time. Utilizing the MERRA2 reanalysis and a fully coupled chemistry–climate simulation with the Community Earth System Model (CESM1(WACCM)) of the National Center for Atmospheric Research (NCAR), we find two downward wave reflection effects on ozone: (1) the direct effect in which the residual circulation is weakened during winter, reducing the typical increase of ozone due to upward planetary wave events and (2) the indirect effect in which the modification of polar temperature during winter affects the amount of ozone destruction in spring. Winter seasons dominated by downward wave reflection events (i.e., reflective winters) are characterized by lower Arctic ozone concentration, while seasons dominated by increased upward wave events (i.e., absorptive winters) are characterized by relatively higher ozone concentration. This behavior is consistent with the cumulative effects of downward and upward planetary wave events on polar stratospheric ozone via the residual circulation and the polar temperature in winter. The results establish a new perspective on dynamical processes controlling stratospheric ozone variability in the Arctic by highlighting the key role of wave reflection.

Description: It is an open question how localized elevated emissions of bromoform (CHBr3) and other very short-lived halocarbons (VSLHs), found in coastal and upwelling regions, and low background emissions, typically found over the open ocean, impact the atmospheric VSLH distribution. In this study, we use the Lagrangian dispersion model FLEXPART to simulate atmospheric CHBr3 resulting from assumed uniform background emissions, and from elevated emissions consistent with those derived during three tropical cruise campaigns. The simulations demonstrate that the atmospheric CHBr3 distributions in the uniform background emissions scenario are highly variable with high mixing ratios appearing in regions of convergence or low wind speed. This relation holds on regional and global scales. The impact of localized elevated emissions on the atmospheric CHBr3 distribution varies significantly from campaign to campaign. The estimated impact depends on the strength of the emissions and the meteorological conditions. In the open waters of the western Pacific and Indian oceans, localized elevated emissions only slightly increase the background concentrations of atmospheric CHBr3, even when 1∘ wide source regions along the cruise tracks are assumed. Near the coast, elevated emissions, including hot spots up to 100 times larger than the uniform background emissions, can be strong enough to be distinguished from the atmospheric background. However, it is not necessarily the highest hot spot emission that produces the largest enhancement, since the tug-of-war between fast advective transport and local accumulation at the time of emission is also important. Our results demonstrate that transport variations in the atmosphere itself are sufficient to produce highly variable VSLH distributions, and elevated VSLHs in the atmosphere do not always reflect a strong localized source. Localized elevated emissions can be obliterated by the highly variable atmospheric background, even if they are orders of magnitude larger than the average open ocean emissions.

Description: State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

Description: Oceanic very short-lived substances (VSLSs), such as bromoform (CHBr3), contribute to stratospheric halogen loading and, thus, to ozone depletion. However, the amount, timing, and region of bromine delivery to the stratosphere through one of the main entrance gates, the Indian summer monsoon circulation, are still uncertain. In this study, we created two bromoform emission inventories with monthly resolution for the tropical Indian Ocean and west Pacific based on new in situ bromoform measurements and novel ocean biogeochemistry modeling. The mass transport and atmospheric mixing ratios of bromoform were modeled for the year 2014 with the particle dispersion model FLEXPART driven by ERA-Interim reanalysis. We compare results between two emission scenarios: (1) monthly averaged and (2) annually averaged emissions. Both simulations reproduce the atmospheric distribution of bromoform from ship- and aircraft-based observations in the boundary layer and upper troposphere above the Indian Ocean reasonably well. Using monthly resolved emissions, the main oceanic source regions for the stratosphere include the Arabian Sea and Bay of Bengal in boreal summer and the tropical west Pacific Ocean in boreal winter. The main stratospheric injection in boreal summer occurs over the southern tip of India associated with the high local oceanic sources and strong convection of the summer monsoon. In boreal winter more bromoform is entrained over the west Pacific than over the Indian Ocean. The annually averaged stratospheric injection of bromoform is in the same range whether using monthly averaged or annually averaged emissions in our Lagrangian calculations. However, monthly averaged emissions result in the highest mixing ratios within the Asian monsoon anticyclone in boreal summer and above the central Indian Ocean in boreal winter, while annually averaged emissions display a maximum above the west Indian Ocean in boreal spring. In the Asian summer monsoon anticyclone bromoform atmospheric mixing ratios vary by up to 50% between using monthly averaged and annually averaged oceanic emissions. Our results underline that the seasonal and regional stratospheric bromine injection from the tropical Indian Ocean and west Pacific critically depend on the seasonality and spatial distribution of the VSLS emissions.

Description: Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Description: The climate research community uses atmospheric reanalysis data sets to understand a wide range of processes and variability in the atmosphere, yet different reanalyses may give very different results for the same diagnostics. The Stratosphere–troposphere Processes And their Role in Climate (SPARC) Reanalysis Intercomparison Project (S-RIP) is a coordinated activity to compare reanalysis data sets using a variety of key diagnostics. The objectives of this project are to identify differences among reanalyses and understand their underlying causes, to provide guidance on appropriate usage of various reanalysis products in scientific studies, particularly those of relevance to SPARC, and to contribute to future improvements in the reanalysis products by establishing collaborative links between reanalysis centres and data users. The project focuses predominantly on differences among reanalyses, although studies that include operational analyses and studies comparing reanalyses with observations are also included when appropriate. The emphasis is on diagnostics of the upper troposphere, stratosphere, and lower mesosphere. This paper summarizes the motivation and goals of the S-RIP activity and extensively reviews key technical aspects of the reanalysis data sets that are the focus of this activity. The special issue "The SPARC Reanalysis Intercomparison Project (S-RIP)" in this journal serves to collect research with relevance to the S-RIP in preparation for the publication of the planned two (interim and full) S-RIP reports.

Description: Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and inform on the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimated the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes and the forcings employed. This study presents results from 18 simulations from 15 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100, forced with different scenarios from the Coupled Model Intercomparison Project Phase 5 (CMIP5) representative of the spread in climate model results. The contribution of the Antarctic ice sheet in response to increased warming during this period varies between −7.8 and 30.0 cm of Sea Level Equivalent (SLE). The evolution of the West Antarctic Ice Sheet varies widely among models, with an overall mass loss up to 21.0 cm SLE in response to changes in oceanic conditions. East Antarctica mass change varies between −6.5 and 16.5 cm SLE, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional mass loss of 8 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 AOGCMs show an overall mass loss of 10 mm SLE compared to simulations done under present-day conditions, with limited mass gain in East Antarctica.

Description: Ocean eddies can both trigger mixing (during their formation and decay) and effectively shield water encompassed from being exchanged with ambient water (throughout their lifetimes). These antagonistic effects of eddies complicate the interpretation of synoptic snapshots typically obtained by ship-based oceanographic measurement campaigns. Here we use a coupled physical–biogeochemical model to explore biogeochemical dynamics within anticyclonic eddies in the eastern tropical South Pacific Ocean. The goal is to understand the diverse biogeochemical patterns that have been observed at the subsurface layers of the anticyclonic eddies in this region. Our model results suggest that the diverse subsurface nutrient patterns within eddies are associated with the presence of water masses of different origins at different depths.

Description: Carbonyl sulfide (OCS) and carbon disulfide (CS2) are volatile sulfur gases that are naturally formed in seawater and exchanged with the atmosphere. OCS is the most abundant sulfur gas in the atmosphere, and CS2 is its most important precursor. They have gained interest due to their direct (OCS) or indirect (CS2 via oxidation to OCS) contribution to the stratospheric sulfate aerosol layer. Furthermore, OCS serves as a proxy to constrain terrestrial CO2 uptake by vegetation. Oceanic emissions of both gases contribute a major part to their atmospheric concentration. Here we present a database of previously published and unpublished, mainly ship-borne measurements in seawater and the marine boundary layer for both gases, available at https://doi.pangaea.de/10.1594/PANGAEA.905430 (Lennartz et al., 2019). The database contains original measurements as well as data digitalized from figures in publications from 42 measurement campaigns, i.e. cruises or time series stations, ranging from 1982 to 2019. OCS data cover all ocean basins except for the Arctic Ocean, as well as all months of the year, while the CS2 dataset shows large gaps in spatial and temporal coverage. Concentrations are consistent across different sampling and analysis techniques for OCS. The database is intended to support the identification of global spatial and temporal patterns and to facilitate the evaluation of model simulations.

Description: We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture.

Description: Parameterizations of surface ocean isoprene concentrations are numerous, despite the lack of source/sink process understanding. Here we present isoprene and related field measurements in the mixed layer from the Indian Ocean and the eastern Pacific Ocean to investigate the production and consumption rates in two contrasting regions, namely oligotrophic open ocean and the coastal upwelling region. Our data show that the ability of different phytoplankton functional types (PFTs) to produce isoprene seems to be mainly influenced by light, ocean temperature, and salinity. Our field measurements also demonstrate that nutrient availability seems to have a direct influence on the isoprene production. With the help of pigment data, we calculate in-field isoprene production rates for different PFTs under varying biogeochemical and physical conditions. Using these new calculated production rates, we demonstrate that an additional significant and variable loss, besides a known chemical loss and a loss due to air–sea gas exchange, is needed to explain the measured isoprene concentration. We hypothesize that this loss, with a lifetime for isoprene between 10 and 100 days depending on the ocean region, is potentially due to degradation or consumption by bacteria.

Description: The project MarParCloud (Marine biological production, organic aerosol Particles and marine Clouds: a process chain) aims to improve our understanding of the genesis, modification and impact of marine organic matter (OM) from its biological production, to its export to marine aerosol particles and, finally, to its ability to act as ice-nucleating particles (INPs) and cloud condensation nuclei (CCN). A field campaign at the Cape Verde Atmospheric Observatory (CVAO) in the tropics in September–October 2017 formed the core of this project that was jointly performed with the project MARSU (MARine atmospheric Science Unravelled). A suite of chemical, physical, biological and meteorological techniques was applied, and comprehensive measurements of bulk water, the sea surface microlayer (SML), cloud water and ambient aerosol particles collected at a ground-based and a mountain station took place. Key variables comprised the chemical characterization of the atmospherically relevant OM components in the ocean and the atmosphere as well as measurements of INPs and CCN. Moreover, bacterial cell counts, mercury species and trace gases were analyzed. To interpret the results, the measurements were accompanied by various auxiliary parameters such as air mass back-trajectory analysis, vertical atmospheric profile analysis, cloud observations and pigment measurements in seawater. Additional modeling studies supported the experimental analysis. During the campaign, the CVAO exhibited marine air masses with low and partly moderate dust influences. The marine boundary layer was well mixed as indicated by an almost uniform particle number size distribution within the boundary layer. Lipid biomarkers were present in the aerosol particles in typical concentrations of marine background conditions. Accumulation- and coarse-mode particles served as CCN and were efficiently transferred to the cloud water. The ascent of ocean-derived compounds, such as sea salt and sugar-like compounds, to the cloud level, as derived from chemical analysis and atmospheric transfer modeling results, denotes an influence of marine emissions on cloud formation. Organic nitrogen compounds (free amino acids) were enriched by several orders of magnitude in submicron aerosol particles and in cloud water compared to seawater. However, INP measurements also indicated a significant contribution of other non-marine sources to the local INP concentration, as (biologically active) INPs were mainly present in supermicron aerosol particles that are not suggested to undergo strong enrichment during ocean–atmosphere transfer. In addition, the number of CCN at the supersaturation of 0.30 % was about 2.5 times higher during dust periods compared to marine periods. Lipids, sugar-like compounds, UV-absorbing (UV: ultraviolet) humic-like substances and low-molecular-weight neutral components were important organic compounds in the seawater, and highly surface-active lipids were enriched within the SML. The selective enrichment of specific organic compounds in the SML needs to be studied in further detail and implemented in an OM source function for emission modeling to better understand transfer patterns, the mechanisms of marine OM transformation in the atmosphere and the role of additional sources. In summary, when looking at particulate mass, we see oceanic compounds transferred to the atmospheric aerosol and to the cloud level, while from a perspective of particle number concentrations, sea spray aerosol (i.e., primary marine aerosol) contributions to both CCN and INPs are rather limited.

Description: Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

Description: During a cruise of R/V METEOR in December 2012 the oceanic sources and emissions of various halogenated trace gases and their mixing ratios in the marine atmospheric boundary layer (MABL) were investigated above the Peruvian upwelling. This study presents novel observations of the three very short lived substances (VSLSs) – bromoform, dibromomethane and methyl iodide – together with high-resolution meteorological measurements, Lagrangian transport and source–loss calculations. Oceanic emissions of bromoform and dibromomethane were relatively low compared to other upwelling regions, while those for methyl iodide were very high. Radiosonde launches during the cruise revealed a low, stable MABL and a distinct trade inversion above acting as strong barriers for convection and vertical transport of trace gases in this region. Observed atmospheric VSLS abundances, sea surface temperature, relative humidity and MABL height correlated well during the cruise. We used a simple source–loss estimate to quantify the contribution of oceanic emissions along the cruise track to the observed atmospheric concentrations. This analysis showed that averaged, instantaneous emissions could not support the observed atmospheric mixing ratios of VSLSs and that the marine background abundances below the trade inversion were significantly influenced by advection of regional sources. Adding to this background, the observed maximum emissions of halocarbons in the coastal upwelling could explain the high atmospheric VSLS concentrations in combination with their accumulation under the distinct MABL and trade inversions. Stronger emissions along the nearshore coastline likely added to the elevated abundances under the steady atmospheric conditions. This study underscores the importance of oceanic upwelling and trade wind systems on the atmospheric distribution of marine VSLS emissions.

Description: Benthic microbial methanogenesis is a known source of methane in marine systems. In most sediments, the majority of methanogenesis is located below the sulfate-reducing zone, as sulfate reducers outcompete methanogens for the major substrates hydrogen and acetate. The coexistence of methanogenesis and sulfate reduction has been shown before and is possible through the usage of noncompetitive substrates by methanogens such as methanol or methylated amines. However, knowledge about the magnitude, seasonality, and environmental controls of this noncompetitive methane production is sparse. In the present study, the presence of methanogenesis within the sulfate reduction zone (SRZ methanogenesis) was investigated in sediments (0–30 cm below seafloor, cm b.s.f.) of the seasonally hypoxic Eckernförde Bay in the southwestern Baltic Sea. Water column parameters such as oxygen, temperature, and salinity together with porewater geochemistry and benthic methanogenesis rates were determined in the sampling area "Boknis Eck" quarterly from March 2013 to September 2014 to investigate the effect of seasonal environmental changes on the rate and distribution of SRZ methanogenesis, to estimate its potential contribution to benthic methane emissions, and to identify the potential methanogenic groups responsible for SRZ methanogenesis. The metabolic pathway of methanogenesis in the presence or absence of sulfate reducers, which after the addition of a noncompetitive substrate was studied in four experimental setups: (1) unaltered sediment batch incubations (net methanogenesis), (2) 14C-bicarbonate labeling experiments (hydrogenotrophic methanogenesis), (3) manipulated experiments with the addition of either molybdate (sulfate reducer inhibitor), 2-bromoethanesulfonate (methanogen inhibitor), or methanol (noncompetitive substrate, potential methanogenesis), and (4) the addition of 13C-labeled methanol (potential methylotrophic methanogenesis). After incubation with methanol, molecular analyses were conducted to identify key functional methanogenic groups during methylotrophic methanogenesis. To also compare the magnitudes of SRZ methanogenesis with methanogenesis below the sulfate reduction zone (〉 30 cm b.s.f.), hydrogenotrophic methanogenesis was determined by 14C-bicarbonate radiotracer incubation in samples collected in September 2013. SRZ methanogenesis changed seasonally in the upper 30 cm b.s.f. with rates increasing from March (0.2 nmol cm−3 d−1) to November (1.3 nmol cm−3 d−1) 2013 and March (0.2 nmol cm−3 d−1) to September (0.4 nmol cm−3 d−1) 2014. Its magnitude and distribution appeared to be controlled by organic matter availability, C / N, temperature, and oxygen in the water column, revealing higher rates in the warm, stratified, hypoxic seasons (September–November) compared to the colder, oxygenated seasons (March–June) of each year. The majority of SRZ methanogenesis was likely driven by the usage of noncompetitive substrates (e.g., methanol and methylated compounds) to avoid competition with sulfate reducers, as was indicated by the 1000–3000-fold increase in potential methanogenesis activity observed after methanol addition. Accordingly, competitive hydrogenotrophic methanogenesis increased in the sediment only below the depth of sulfate penetration (〉 30 cm b.s.f.). Members of the family Methanosarcinaceae, which are known for methylotrophic methanogenesis, were detected by PCR using Methanosarcinaceae-specific primers and are likely to be responsible for the observed SRZ methanogenesis. The present study indicates that SRZ methanogenesis is an important component of the benthic methane budget and carbon cycling in Eckernförde Bay. Although its contributions to methane emissions from the sediment into the water column are probably minor, SRZ methanogenesis could directly feed into methane oxidation above the sulfate–methane transition zone.

Description: We studied the concurrence of methanogenesis and sulfate reduction in surface sediments (0–25 cm below sea floor, cmbsf) at six stations (70, 145, 253, 407, 770 and 1024 m) along the Peruvian margin (12° S). This oceanographic region is characterized by high carbon export to the seafloor, creating an extensive oxygen minimum zone (OMZ) on the shelf, both factors that could favor surface methanogenesis. Sediments sampled along the depth transect traversed areas of anoxic and oxic conditions in the bottom-near water. Net methane production (batch incubations) and sulfate reduction (35S-sulfate radiotracer incubation) were determined in the upper 0–25 cmbsf of multicorer cores from all stations, while deep hydrogenotrophic methanogenesis (〉 30 cmbsf, 14C-bicarbonate radiotracer incubation) was determined in two gravity cores at selected sites (78 and 407 m). Furthermore, stimulation (methanol addition) and inhibition (molybdate addition) experiments were carried out to investigate the relationship between sulfate reduction and methanogenesis. Highest rates of methanogenesis and sulfate reduction in the surface sediments, integrated over 0–25 cmbsf, were observed on the shelf (70–253 m, 0.06–0.1 and 0.5–4.7 mmol m−2 d−1, respectively), while lowest rates were discovered at the deepest site (1024 m, 0.03 and 0.2 mmol m−2 d−1, respectively). The addition of methanol resulted in significantly higher surface methanogenesis activity, suggesting that the process was mostly based on non-competitive substrates, i.e., substrates not used by sulfate reducers. In the deeper sediment horizons, where competition was probably relieved due to the decline of sulfate, the usage of competitive substrates was confirmed by the detection of hydrogenotrophic activity in the sulfate-depleted zone at the shallow shelf station (70 m). Surface methanogenesis appeared to be correlated to the availability of labile organic matter (C / N ratio) and organic carbon degradation (DIC production), both of which support the supply of methanogenic substrates. A negative correlation of methanogenesis rates with dissolved oxygen in the bottom-near water was not obvious, however, anoxic conditions within the OMZ might be advantageous for methanogenic organisms at the sediment–water interface. Our results revealed a high relevance of surface methanogenesis on the shelf, where the ratio between surface to deep (below sulfate penetration) methanogenic activity ranged between 0.13 and 105. In addition, methane concentration profiles indicate a partial release of surface methane into the water column as well as a partial consumption of methane by anaerobic methane oxidation (AOM) in the surface sediment. The present study suggests that surface methanogenesis might play a greater role in benthic methane budgeting than previously thought, especially for fueling AOM above the sulfate-methane transition zone.

Description: Marine-produced short-lived trace gases such as dibromomethane (CH2Br2), bromoform (CHBr3), methyliodide (CH3I) and dimethyl sulfide (DMS) significantly impact tropospheric and stratospheric chemistry. Describing their marine emissions in atmospheric chemistry models as accurately as possible is necessary to quantify their impact on ozone depletion and Earth's radiative budget. So far, marine emissions of trace gases have mainly been prescribed from emission climatologies, thus lacking the interaction between the actual state of the atmosphere and the ocean. Here we present simulations with the chemistry climate model EMAC (ECHAM5/MESSy Atmospheric Chemistry) with online calculation of emissions based on surface water concentrations, in contrast to directly prescribed emissions. Considering the actual state of the model atmosphere results in a concentration gradient consistent with model real-time conditions at the ocean surface and in the atmosphere, which determine the direction and magnitude of the computed flux. This method has a number of conceptual and practical benefits, as the modelled emission can respond consistently to changes in sea surface temperature, surface wind speed, sea ice cover and especially atmospheric mixing ratio. This online calculation could enhance, dampen or even invert the fluxes (i.e. deposition instead of emissions) of very short-lived substances (VSLS). We show that differences between prescribing emissions and prescribing concentrations (−28 % for CH2Br2 to +11 % for CHBr3) result mainly from consideration of the actual, time-varying state of the atmosphere. The absolute magnitude of the differences depends mainly on the surface ocean saturation of each particular gas. Comparison to observations from aircraft, ships and ground stations reveals that computing the air–sea flux interactively leads in most of the cases to more accurate atmospheric mixing ratios in the model compared to the computation from prescribed emissions. Calculating emissions online also enables effective testing of different air–sea transfer velocity (k) parameterizations, which was performed here for eight different parameterizations. The testing of these different k values is of special interest for DMS, as recently published parameterizations derived by direct flux measurements using eddy covariance measurements suggest decreasing k values at high wind speeds or a linear relationship with wind speed. Implementing these parameterizations reduces discrepancies in modelled DMS atmospheric mixing ratios and observations by a factor of 1.5 compared to parameterizations with a quadratic or cubic relationship to wind speed.

Description: Carbonyl sulfide (OCS) is the most abundant, long-lived sulphur gas in the atmosphere and a major supplier of sulfur to the stratospheric sulfate aerosol layer. The short-lived gas carbon disulfide (CS2) is oxidized to OCS and constitutes a major indirect source to the atmospheric OCS budget. The atmospheric budget of OCS is not well constrained due to a large missing source needed to compensate for substantial evidence that was provided for significantly higher sinks. Oceanic emissions are associated with major uncertainties. Here we provide a first, monthly resolved ocean emission inventory of both gases for the period 2000–2019 (available at https://doi.org/10.5281/zenodo.4297010) (Lennartz et al., 2020a). Emissions are calculated with a numerical box model (resolution 2.8° × 2.8° at equator, T42 grid) for the surface mixed layer. We find that interannual variability in OCS emissions is smaller than seasonal variability, and is mainly driven by variations in chromophoric dissolved organic matter (CDOM), which influences both photochemical and light-independent production. A comparison with a global database of more than 2500 measurements reveals overall good agreement. Emissions of CS2 constitute a larger sulfur source to the atmosphere than OCS, and equally show interannual variability connected to variability of CDOM. The emission estimate of CS2 is associated with higher uncertainties, as process understanding of the marine cycling of CS2 is incomplete. We encourage the use of the data provided here as input for atmospheric modelling studies to further assess the atmospheric OCS budget and the role of OCS in climate.

Description: Global models of the oceanic nitrogen cycle are subject to many uncertainties, among them type and form of biogeochemical processes involved in the fixed nitrogen cycle, and the spatial and temporal scales, on which the global nitrogen budget is regulated. We investigate these aspects using a global model of ocean biogeochemistry, that explicitly considers phosphorus and nitrogen, including pelagic denitrification and nitrogen fixation as sink and source terms of fixed nitrogen, respectively. The model explores different parameterizations of organic matter sinking speed, oxidant affinity of oxic and suboxic remineralization, and regulation of nitrogen fixation by temperature and different stoichiometric ratios. Examination of the initial transient behaviour of different model setups initialized from observed tracer distributions reveal changes in simulated nitrogen inventories and fluxes particularly during the first centuries. Millennial timescales have to be resolved in order to bring all biogeochemical and physical processes into a dynamically consistent steady state, for which global patterns of biogeochemical tracers and fluxes are reproduced quite well. Analysis of global properties suggests that particularly particle sinking speed, but also the parameterization of denitrification determines the extent of oxygen minimum zones, global nitrogen fluxes, and hence the oceanic nitrogen inventory. However, the ways and directions, in which different parameterizations of particle sinking, nitrogen fixation and denitrification affect the global diagnostics, are different, suggesting that these may, in principle, be constrained independently from each other. Analysis of the model misfit suggests a particle flux profile close to the one suggested by Martin et al. (1987). Simulated pelagic denitrification best agrees with the lower values between 59 and 84 Tg N yr−1 recently estimated by other authors.

Description: The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, al-though there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-corederived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Description: O2 minimum zones (OMZ) of the world's oceans are important locations for microbial dissimilatory NO3- reduction and subsequent loss of combined nitrogen (N) to biogenic N2 gas. This is particularly so when the OMZ is coupled to a region of high productivity leading to high rates of N-loss as found in the coastal upwelling region off Peru. Stable N isotope ratios (and O in the case of NO3- and NO2-) can be used as natural tracers of OMZ N-cycling because of distinct kinetic isotope effects associated with microbially-mediated N-cycle transformations. Here we present NO2- and NO3- stable isotope data from the nearshore upwelling region off Callao, Peru. Subsurface O2 was generally depleted below about 30 m depth with O2 less than 10 μM, while NO2- concentrations were high, ranging from 6 to 10 μM and NO3- was in places strongly depleted to near 0 μM. We observed for the first time, a positive linear relationship between NO2- δ15N and δ18O at our coastal stations, analogous to that of NO3- N and O isotopes during assimilatory and dissimilatory reduction. This relationship is likely the result of rapid NO2- turnover due to higher organic matter flux in these coastal upwelling waters. No such relationship was observed at offshore stations where slower turnover of NO2- facilitates dominance of isotope exchange with water. We also evaluate the overall isotope fractionation effect for N-loss in this system using several approaches that vary in their underlying assumptions. While there are differences in apparent fractionation factor (ε) for N-loss as calculated from the δ15N of [NO3-], DIN, or biogenic N2, values for ε are generally much lower than previously reported, reaching as low as 6.5‰. A possible explanation is the influence of sedimentary N-loss at our inshore stations which incurs highly suppressed isotope fractionation.

Description: Upwelling of nutrient-rich deep waters make eastern boundary upwelling systems (EBUSs), such as the Humboldt Current system, hot spots of marine productivity. Associated settling of organic matter to depth and consecutive aerobic decomposition results in large subsurface water volumes being oxygen depleted. Under these circumstances, organic matter remineralisation can continue via denitrification, which represents a major loss pathway for bioavailable nitrogen. Additionally, anaerobic ammonium oxidation can remove significant amounts of nitrogen in these areas. Here we assess the interplay of suboxic water upwelling and nitrogen cycling in a manipulative offshore mesocosm experiment. Measured denitrification rates in incubations with water from the oxygen-depleted bottom layer of the mesocosms (via 15N label incubations) mostly ranged between 5.5 and 20 nmol N2 L−1 h−1 (interquartile range), reaching up to 80 nmol N2 L−1 h−1. However, actual in situ rates in the mesocosms, estimated via Michaelis–Menten kinetic scaling, did most likely not exceed 0.2–4.2 nmol N2 L−1 h−1 (interquartile range) due to substrate limitation. In the surrounding Pacific, measured denitrification rates were similar, although indications of substrate limitation were detected only once. In contrast, anammox (anaerobic ammonium oxidation) made only a minor contribution to the overall nitrogen loss when encountered in both the mesocosms and the Pacific Ocean. This was potentially related to organic matter C / N stoichiometry and/or process-specific oxygen and hydrogen sulfide sensitivities. Over the first 38 d of the experiment, total nitrogen loss calculated from in situ rates of denitrification and anammox was comparable to estimates from a full nitrogen budget in the mesocosms and ranged between ∼ 1 and 5.5 µmol N L−1. This represents up to ∼  20 % of the initially bioavailable inorganic and organic nitrogen standing stocks. Interestingly, this loss is comparable to the total amount of particulate organic nitrogen that was exported into the sediment traps at the bottom of the mesocosms at about 20 m depth. Altogether, this suggests that a significant portion, if not the majority of nitrogen that could be exported to depth, is already lost, i.e. converted to N2 in a relatively shallow layer of the surface ocean, provided that there are oxygen-deficient conditions like those during coastal upwelling in our study. Published data for primary productivity and nitrogen loss in all EBUSs reinforce such conclusion.

Description: Estuaries are sources of nitrous oxide (N2O) and methane (CH4) to the atmosphere. However, our present knowledge of N2O and CH4 emissions from estuaries in the tropics is very limited because data are scarce. In this study, we present first measurements of dissolved N2O and CH4 from two estuaries in a peat-dominated region of northwestern Borneo. Two campaigns (during the dry season in June 2013 and during the wet season in March 2014) were conducted in the estuaries of the Lupar and Saribas rivers. Median N2O concentrations ranged between 7.2 and 12.3 nmol L−1 and were higher in the marine end-member (13.0 ± 7.0 nmol L−1). CH4 concentrations were low in the coastal ocean (3.6 ± 0.2 nmol L−1) and higher in the estuaries (medians between 10.6 and 64.0 nmol L−1). The respiration of abundant organic matter and presumably anthropogenic input caused slight eutrophication, which did not lead to hypoxia or enhanced N2O concentrations, however. Generally, N2O concentrations were not related to dissolved inorganic nitrogen concentrations. Thus, the use of an emission factor for the calculation of N2O emissions from the inorganic nitrogen load leads to an overestimation of the flux from the Lupar and Saribas estuaries. N2O was negatively correlated with salinity during the dry season, which suggests a riverine source. In contrast, N2O concentrations during the wet season were not correlated with salinity but locally enhanced within the estuaries, implying that there were additional estuarine sources during the wet (i.e., monsoon) season. Estuarine CH4 distributions were not driven by freshwater input but rather by tidal variations. Both N2O and CH4 concentrations were more variable during the wet season. We infer that the wet season dominates the variability of the N2O and CH4 concentrations and subsequent emissions from tropical estuaries. Thus, we speculate that any changes in the Southeast Asian monsoon system will lead to changes in the N2O and CH4 emissions from these systems. We also suggest that the ongoing cultivation of peat soil in Borneo is likely to increase N2O emissions from these estuaries, while the effect on CH4 remains uncertain.

Description: The northwestern Tropical Atlantic Ocean is a turbulent region, filled with mesoscale eddies and regional currents. In this intense dynamical context, several water masses with thermohaline characteristics of different origins are advected, mixed, and stirred at the surface and at depth. The EUREC4A-OA/ATOMIC experiment that took place in January and February 2020 was dedicated to assessing the processes at play in this region, especially the interaction between the ocean and the atmosphere. For that reason, four oceanographic vessels and different autonomous platforms measured properties near the air–sea interface and acquired thousands of upper-ocean (up to 400–2000 m depth) profiles. However, each device had its own observing capability, varying from deep measurements acquired during vessel stations to shipboard underway near-surface observations and measurements from autonomous and uncrewed systems (such as Saildrones). These observations were undertaken with a specific sampling strategy guided by near-real-time satellite maps and adapted every half day, based on the process that was investigated. These processes were characterized by different spatiotemporal scales, from mesoscale eddies, with diameters exceeding 100 km, to submesoscale filaments of 1 km width. This article describes the datasets gathered from the different devices and how the data were calibrated and validated. In order to ensure an overall consistency, the platforms' datasets are cross-validated using a hierarchy of instruments defined by their own specificity and calibration procedures. This has enabled the quantification of the uncertainty in the measured parameters when different datasets are used together, e.g., https://doi.org/10.17882/92071 (L'Hégaret et al., 2020a).

Description: At present, anthropogenic halogens and oceanic emissions of Very Short-Lived Substances (VSLS) are responsible for stratospheric ozone destruction. Emissions of the, mostly long-lived, anthropogenic halogens have been reduced, and as a consequence, their atmospheric abundance has started to decline since the beginning of the 21st century. Emissions of VSLS are, on the other hand, expected to increase in the future. VSLS are known to have large natural sources; however increasing evidence arises that their oceanic production and emission is enhanced by anthropogenic activities. Here, we introduce a new approach of assessing the overall impact of all oceanic halogen emissions on stratospheric ozone by calculating Ozone Depletion Potential (ODP)-weighted emissions of VSLS. Seasonally and spatially dependent, global distributions are derived exemplary for CHBr3 for the period 1999–2006. At present, ODP-weighted emissions of CHBr3 amount up to 50% of ODP-weighted anthropogenic emissions of CFC-11 and to 9% of all long-lived ozone depleting substances. The ODP-weighted emissions are large where strong oceanic emissions coincide with high-reaching convective activity and show pronounced peaks at the equator and the coasts with largest contributions from the Maritime Continent and West Pacific. Variations of tropical convective activity lead to seasonal shifts in the spatial distribution of the ODP with the updraught mass flux explaining 71% of the variance of the ODP distribution. Future climate projections based on RCP8.5 scenario suggest a 31% increase of the ODP-weighted CHBr3 emissions until 2100 compared to present values. This increase is related to larger convective activity and increasing emissions in a future climate; however, is reduced at the same time by less effective bromine-related ozone depletion. The comparison of the ODP-weighted emissions of short and long-lived halocarbons provides a new concept for assessing the overall impact of oceanic bromine emissions on stratospheric ozone depletion for current conditions and future projections.

Description: Magmatic sill intrusions into organic-rich sediments cause the release of thermogenic CH4 and CO2. Pore fluids from the Guaymas Basin (Gulf of California), a sedimentary basin with recent magmatic activity, were investigated to constrain the link between sill intrusions and fluid seepage as well as the timing of sill-induced hydrothermal activity. Sampling sites were close to a hydrothermal vent field at the northern rift axis and at cold seeps located up to 30km away from the rift. Pore fluids close to the active hydrothermal vent field showed a slight imprint by hydrothermal fluids and indicated a shallow circulation system transporting seawater to the hydrothermal catchment area. Geochemical data of pore fluids at cold seeps showed a mainly ambient diagenetic fluid composition without any imprint related to high temperature processes at greater depth. Seep communities at the seafloor were mainly sustained by microbial methane, which rose along pathways formed earlier by hydrothermal activity, driving the anaerobic oxidation of methane (AOM) and the formation of authigenic carbonates. Overall, our data from the cold seep sites suggest that at present, sill-induced hydrothermalism is not active away from the ridge axis, and the vigorous venting of hydrothermal fluids is restricted to the ridge axis. Using the sediment thickness above extinct conduits and carbonate dating, we calculated that deep fluid and thermogenic gas flow ceased 28 to 7kyr ago. These findings imply a short lifetime of hydrothermal systems, limiting the time of unhindered carbon release as suggested in previous modeling studies. Consequently, activation and deactivation mechanisms of these systems need to be better constrained for the use in climate modeling approaches.

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: In recent updates of the HITRAN water vapour H2O spectroscopic compilation covering the blue spectral region (here: 394–480 nm) significant changes for the absorption bands at 416 and 426 nm were reported. In order to investigate the consistency of the different cross-sections calculated from these compilations, H2O vapour column density ratios for different spectral intervals were retrieved from long-path and multi-axis differential optical absorption spectroscopy (DOAS) measurements. We observed a significant improvement of the DOAS evaluation when using the updated HITRAN water vapour absorption cross-sections for the calculation of the reference spectra. In particular the magnitudes of the residual spectra as well as the fit errors were reduced. However, we also found that the best match between measurement and model is reached when the absorption cross-section of groups of lines are scaled by factors ranging from 0.5 to 1.9, suggesting that the HITRAN water vapour absorption compilation still needs significant corrections. For this spectral region we present correction factors for HITRAN 2009, HITRAN 2012, HITEMP and BT2 derived from field measurements. Additionally, upper limits for water vapour absorption in the UV-A range from 330 to 390 nm are given.

Description: A suite of oxygenated volatile organic compounds (OVOCs – acetaldehyde, acetone, propanal, butanal and butanone) were measured concurrently in the surface water and atmosphere of the South China Sea and Sulu Sea in November 2011. A strong correlation was observed between all OVOC concentrations in the surface seawater along the entire cruise track, except for acetaldehyde, suggesting similar sources and sinks in the surface ocean. Additionally, several phytoplankton groups, such as haptophytes or pelagophytes, were also correlated to all OVOCs, indicating that phytoplankton may be an important source of marine OVOCs in the South China and Sulu seas. Humic- and protein-like fluorescent dissolved organic matter (FDOM) components seemed to be additional precursors for butanone and acetaldehyde. The measurement-inferred OVOC fluxes generally showed an uptake of atmospheric OVOCs by the ocean for all gases, except for butanal. A few important exceptions were found along the Borneo coast, where OVOC fluxes from the ocean to the atmosphere were inferred. The atmospheric OVOC mixing ratios over the northern coast of Borneo were relatively high compared with literature values, suggesting that this coastal region is a local hotspot for atmospheric OVOCs. The calculated amount of OVOCs entrained into the ocean seemed to be an important source of OVOCs to the surface ocean. When the fluxes were out of the ocean, marine OVOCs were found to be enough to control the locally measured OVOC distribution in the atmosphere. Based on our model calculations, at least 0.4 ppb of marine-derived acetone and butanone can reach the upper troposphere, where they may have an important influence on hydrogen oxide radical formation over the western Pacific Ocean.

Description: We evaluated NOVAC (Network for Observation of Volcanic and Atmospheric Change) gas emission data from the 2015 eruption of Cotopaxi volcano (Ecuador) for BrO/SO2 molar ratios. Statistical analysis of the data revealed a conspicuous periodic pattern with a periodicity of about two weeks in a three month time series. While the time series is too short to rule out a chance recurrence of transient geological or meteorological events as a possible origin for the periodic signal, we nevertheless took this observation as a motivation to examine the influence of natural forcings with periodicities of around two weeks on volcanic gas emissions. One strong aspirant with such a periodicity are the Earth tides, which are thus central in this study. We present the BrO/SO2 data, analyse the reliability of the periodic signal, discuss a possible meteorological or eruption-induced origin of this signal, and compare the signal with the theoretical ground surface displacement pattern caused by the Earth tides. Central result is the observation of a significant correlation between the BrO/SO2 molar ratios with the North-South and vertical components of the calculated tide-induced surface displacement with correlation coefficients of 47 % and 36 %, respectively. From all other investigated parameters, only the correlation between the BrO/SO2 molar ratios and the relative humidity in the local atmosphere resulted in a comparable correlation coefficient of about 33 %.

Description: Although coccolithophore physiological responses to CO2-induced changes in seawater carbonate chemistry have been widely studied in the past, there is limited knowledge on the variability of physiological responses between populations from different areas. In the present study, we investigated the specific responses of growth, particulate organic (POC) and inorganic carbon (PIC) production rates of three populations of the coccolithophore Emiliania huxleyi from three regions in the North Atlantic Ocean (Azores: six strains, Canary Islands: five strains, and Norwegian coast near Bergen: six strains) to a CO2 partial pressure (pCO2) range from 120 to 2630 µatm. Physiological rates of each population and individual strain increased with rising pCO2 levels, reached a maximum and declined thereafter. Optimal pCO2 for growth, POC production rates, and tolerance to low pH (i.e., high proton concentration) was significantly higher in an E. huxleyi population isolated from the Norwegian coast than in those isolated near the Azores and Canary Islands. This may be due to the large environmental variability including large pCO2 and pH fluctuations in coastal waters off Bergen compared to the rather stable oceanic conditions at the other two sites. Maximum growth and POC production rates of the Azores and Bergen populations were similar and significantly higher than that of the Canary Islands population. This pattern could be driven by temperature–CO2 interactions where the chosen incubation temperature (16 °C) was slightly below what strains isolated near the Canary Islands normally experience. Our results indicate adaptation of E. huxleyi to their local environmental conditions and the existence of distinct E. huxleyi populations. Within each population, different growth, POC, and PIC production rates at different pCO2 levels indicated strain-specific phenotypic plasticity. Accounting for this variability is important to understand how or whether E. huxleyi might adapt to rising CO2 levels.

Description: To mitigate the rumen enteric methane (CH4) produced by ruminant livestock, Asparagopsis taxiformis is proposed as an additive to ruminant feed. During the cultivation of Asparagopsis taxiformis in the sea or in terrestrially based systems, this macroalgae, like most seaweeds and phytoplankton, produces a large amount of bromoform (CHBr3), which contributes to ozone depletion once released into the atmosphere. In this study, we focus on the impact of CHBr3 on the stratospheric ozone layer resulting from potential emissions from proposed Asparagopsis cultivation in Australia. The impact is assessed by weighting the emissions of CHBr3 with its ozone depletion potential (ODP), which is traditionally defined for long-lived halocarbons but has also been applied to very short-lived substances (VSLSs). An annual yield of ∼3.5 × 104 Mg dry weight is required to meet the needs of 50 % of the beef feedlot and dairy cattle in Australia. Our study shows that the intensity and impact of CHBr3 emissions vary, depending on location and cultivation scenarios. Of the proposed locations, tropical farms near the Darwin region are associated with the largest CHBr3 ODP values. However, farming of Asparagopsis using either ocean or terrestrial cultivation systems at any of the proposed locations does not have the potential to significantly impact the ozone layer. Even if all Asparagopsis farming were performed in Darwin, the CHBr3 emitted into the atmosphere would amount to less than 0.02 % of the global ODP-weighted emissions. The impact of remaining farming scenarios is also relatively small even if the intended annual yield in Darwin is scaled by a factor of 30 to meet the global requirements, which will increase the global ODP-weighted emissions up to ∼0.5 %.

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: For the past decade, observations of carbonyl sulfide (OCS or COS) have been investigated as a proxy for carbon uptake by plants. OCS is destroyed by enzymes that interact with CO2 during photosynthesis, namely carbonic anhydrase (CA) and RuBisCO, where CA is the more important one. The majority of sources of OCS to the atmosphere are geographically separated from this large plant sink, whereas the sources and sinks of CO2 are co-located in ecosystems. The drawdown of OCS can therefore be related to the uptake of CO2 without the added complication of co-located emissions comparable in magnitude. Here we review the state of our understanding of the global OCS cycle and its applications to ecosystem carbon cycle science. OCS uptake is correlated well to plant carbon uptake, especially at the regional scale. OCS can be used in conjunction with other independent measures of ecosystem function, like solar-induced fluorescence and carbon and water isotope studies. More work needs to be done to generate global coverage for OCS observations and to link this powerful atmospheric tracer to systems where fundamental questions concerning the carbon and water cycle remain.

Description: In the Baltic Sea, salinity and its large variability, both horizontal and vertical, are key physical factors in determining the overall stratification conditions. In addition to that, salinity and its changes also have large effects on various ecosystem processes. Several factors determine the observed two-layer vertical structure of salinity. Due to the excess of river runoff to the sea, there is a continuous outflow of water masses in the surface layer with a compensating inflow to the Baltic in the lower layer. Also, the net precipitation plays a role in the water balance and consequently in the salinity dynamics. The salinity conditions in the sea are also coupled with the changes in the meteorological conditions. The ecosystem is adapted to the current salinity level: a change in the salinity balance would lead to ecological stress of flora and fauna, and related negative effects on possibilities to carry on sustainable development of the ecosystem. The Baltic Sea salinity regime has been studied for more than 100 years. In spite of that, there are still gaps in our knowledge of the changes of salinity in space and time. An important part of our understanding of salinity are its long-term changes. However, the available scenarios for the future development of salinity are still inaccurate. We still need more studies on various factors related to salinity dynamics. Among others more knowledge is needed, e.g. from meteorological patterns in various space and time scales and mesoscale variability in precipitation. Also, updated information on river runoff and inflows of saline water is needed to close the water budget. We still do not understand accurately enough the water mass exchange between North Sea and Baltic Sea and within its sub-basins. Scientific investigations of the complicated vertical mixing processes are additionally required. This paper is a continuation and update of the BACC II book which was published in 2015, including information from articles issued until 2012. After that, there have been many new publications on the salinity dynamics, not least because of the Major Baltic Inflow which took place in December 2014. Several key topics have been investigated, including the coupling of long-term variations of climate with the observed salinity changes. Here the focus is on observing and indicating the role of climate change for salinity dynamics. New results of MBI-dynamics and related water mass interchange between the Baltic Sea and the North Sea have been published. Those studies also included results from the MBI-related meteorological conditions, variability in salinity and exchange of water masses between various scales. All these processes are in turn coupled with changes in the Baltic Sea circulation dynamics.

Description: Methane (CH4) is a climate-relevant atmospheric trace gas which is emitted to the atmosphere from coastal areas such as the Baltic Sea. The oceanic CH4 emission estimates are still associated with a high degree of uncertainty partly because the temporal and spatial variability in the CH4 distribution in the ocean surface layer is usually not known. In order to determine the small-scale variability in dissolved CH4 we set up a purge and trap system with a significantly improved precision for the CH4 concentration measurements compared to static headspace equilibration measurements. We measured the distribution of dissolved CH4 in the water column of the western Kiel Bight and Eckernförde Bay in June and September 2018. The top 1 m was sampled in high resolution to determine potential small-scale CH4 concentration gradients within the mixed layer. CH4 concentrations throughout the water column of the western Kiel Bight and Eckernförde Bay were generally higher in September than in June. The increase in the CH4 concentrations in the bottom water was accompanied by a strong decrease in O2 concentrations which led to anoxic conditions favourable for microbial CH4 production in September. In summer 2018, northwestern Europe experienced a pronounced heatwave. However, we found no relationship between the anomalies of water temperature and excess CH4 in both the surface and the bottom layer at the site of the Boknis Eck Time Series Station (Eckernförde Bay). Therefore, the 2018 European heatwave most likely did not affect the observed increase in the CH4 concentrations in the western Kiel Bight from June to September 2018. The high-resolution measurements of the CH4 concentrations in the upper 1 m of the water column were highly variable and showed no uniform decreasing or increasing gradients with water depth. Overall, our results show that the CH4 distribution in the water column of the western Kiel Bight and Eckernförde Bay is strongly affected by both large-scale temporal (i.e. seasonal) and small-scale spatial variabilities which need to be considered when quantifying the exchange of CH4 across the ocean–atmosphere interface.

Description: The reaction between ozone and iodide at the sea surface is now known to be an important part of atmospheric ozone cycling, causing ozone deposition and the release of ozone-depleting reactive iodine to the atmosphere. The importance of this reaction is reflected by its inclusion in chemical transport models (CTMs). Such models depend on accurate sea surface iodide fields, but measurements are spatially and temporally limited. Hence, the ability to predict current and future sea surface iodide fields, i.e. sea surface iodide concentration on a narrow global grid, requires the development of process-based models. These models require a thorough understanding of the key processes that control sea surface iodide. The aim of this study was to explore if there are common features of iodate-to-iodide reduction amongst diverse marine phytoplankton in order to develop models that focus on sea surface iodine and iodine release to the troposphere. In order to achieve this, rates and patterns of changes in inorganic iodine speciation were determined in 10 phytoplankton cultures grown at ambient iodate concentrations. Where possible these data were analysed alongside results from previous studies. Iodate loss and some iodide production were observed in all cultures studied, confirming that this is a widespread feature amongst marine phytoplankton. We found no significant difference in log-phase, cell-normalised iodide production rates between key phytoplankton groups (diatoms, prymnesiophytes including coccolithophores and phaeocystales), suggesting that a phytoplankton functional type (PFT) approach would not be appropriate for building an ocean iodine cycling model. Iodate loss was greater than iodide formation in the majority of the cultures studied, indicating the presence of an as-yet-unidentified “missing iodine” fraction. Iodide yield at the end of the experiment was significantly greater in cultures that had reached a later senescence stage. This suggests that models should incorporate a lag between peak phytoplankton biomass and maximum iodide production and that cell mortality terms in biogeochemical models could be used to parameterise iodide production.

Description: Southern hemisphere lower stratospheric ozone depletion has been shown to lead to a poleward shift of the tropospheric jet stream during austral summer, influencing surface atmosphere and ocean conditions, such as surface temperatures and sea ice extent. The characteristics of stratospheric and tropospheric responses to ozone depletion, however, differ largely among climate models depending on the representation of ozone in the models. The most accurate way to represent ozone in a model is to calculate it interactively. However, due to computational costs, in particular for long-term coupled ocean-atmosphere model integrations, the more common way is to prescribe ozone from observations or calculated model fields. Here, we investigate the difference between an interactive and a specified chemistry version of the same atmospheric model in a fully-coupled setup using a 9-member chemistry-climate model ensemble. In the specified chemistry version of the model the ozone fields are prescribed using the output from the interactive chemistry model version. In contrast to earlier studies, we use daily-resolved ozone fields in the specified chemistry simulations to achieve a better comparability between the ozone forcing with and without interactive chemistry. We find that although the short-wave heating rate trend in response to ozone depletion is the same in the different chemistry settings, the interactive chemistry ensemble shows a stronger trend in polar cap stratospheric temperatures (by about 0.7 K per decade) and circumpolar stratospheric zonal mean zonal winds (by about 1.6 m/s per decade) as compared to the specified chemistry ensemble. This difference between interactive and specified chemistry in the stratospheric response to ozone depletion also affects the tropospheric response, namely the poleward shift of the tropospheric jet stream. We attribute part of these differences to the missing representation of feedbacks between chemistry and dynamics in the specified chemistry ensemble, which affect the dynamical heating rates, and part of it to the lack of spatial asymmetries in the prescribed ozone fields. This effect is investigated using a sensitivity ensemble that was forced by a three-dimensional instead of a two–dimensional ozone field. This study emphasizes the value of interactive chemistry for the representation of the southern hemisphere tropospheric jet response to ozone depletion and infers that for periods with strong ozone variability (trends) the details of the ozone forcing can be crucial for representing southern hemispheric climate variability.

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: 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: This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMIP5 predecessor by using time-slice experiments of two chemistry–climate models and a reference radiative transfer model. The differences in the long-term mean SSI in the CMIP6 dataset, compared to CMIP5, impact on climatological stratospheric conditions (lower shortwave heating rates of −0.35 K day−1 at the stratopause), cooler stratospheric temperatures (−1.5 K in the upper stratosphere), lower ozone abundances in the lower stratosphere (−3 %), and higher ozone abundances (+1.5 % in the upper stratosphere and lower mesosphere). Between the maximum and minimum phases of the 11-year solar cycle, there is an increase in shortwave heating rates (+0.2 K day−1 at the stratopause), temperatures ( ∼  1 K at the stratopause), and ozone (+2.5 % in the upper stratosphere) in the tropical upper stratosphere using the CMIP6 forcing dataset. This solar-cycle response is slightly larger, but not statistically significantly different from that for the CMIP5 forcing dataset. CMIP6 models with a well-resolved shortwave radiation scheme are encouraged to prescribe SSI changes and include solar-induced stratospheric ozone variations, in order to better represent solar climate variability compared to models that only prescribe TSI and/or exclude the solar-ozone response. We show that monthly-mean solar-induced ozone variations are implicitly included in the SPARC/CCMI CMIP6 Ozone Database for historical simulations, which is derived from transient chemistry–climate model simulations and has been developed for climate models that do not calculate ozone interactively. CMIP6 models without chemistry that perform a preindustrial control simulation with time-varying solar forcing will need to use a modified version of the SPARC/CCMI Ozone Database that includes solar variability. CMIP6 models with interactive chemistry are also encouraged to use the particle forcing datasets, which will allow the potential long-term effects of particles to be addressed for the first time. The consideration of particle forcing has been shown to significantly improve the representation of reactive nitrogen and ozone variability in the polar middle atmosphere, eventually resulting in further improvements in the representation of solar climate variability in global models.

Description: Bulk aerosol samples collected during cruise M91 of FS Meteor off the coast of Peru in December 2012 were analysed for their soluble trace metal (Fe, Al, Mn, Ti, Zn, V, Ni, Cu, Co, Cd, Pb, Th) and major ion (including NO3− and NH4+) content. These data are among the first recorded for trace metals in this relatively poorly studied region of the global marine atmosphere. To the north of ∼ 13° S, the concentrations of several elements (Fe, Ti, Zn, V, Ni, Pb) appear to be related to distance from the coast. At the south of the transect (∼ 15–16° S), elevated concentrations of Fe, Cu, Co and Ni were observed. These may be related to the activities of the large smelting facilities in the south of Peru or northern Chile. Calculated dry deposition fluxes (3370–17 800 and 16–107 nmol m−2 d−1 for inorganic nitrogen and soluble Fe respectively) indicated that atmospheric input to the waters of the Peru upwelling system contains an excess of Fe over N, with respect to phytoplankton requirements. This may be significant as primary production in these waters has been reported to be limited by Fe availability, but atmospheric deposition is unlikely to be the dominant source of Fe to the system

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: Volcanic plumes are common and far-reaching manifestations of volcanic activity during and between eruptions. Observations of the rate of emission and composition of volcanic plumes are essential to recognize and, in some cases, predict the state of volcanic activity. Measurements of the size and location of the plumes are important to assess the impact of the emission from sporadic or localized events to persistent or widespread processes of climatic and environmental importance. These observations provide information on volatile budgets on Earth, chemical evolution of magmas, and atmospheric circulation and dynamics. Space-based observations during the last decades have given us a global view of Earth's volcanic emission, particularly of sulfur dioxide (SO2). Although none of the satellite missions were intended to be used for measurement of volcanic gas emission, specially adapted algorithms have produced time-averaged global emission budgets. These have confirmed that tropospheric plumes, produced from persistent degassing of weak sources, dominate the total emission of volcanic SO2. Although space-based observations have provided this global insight into some aspects of Earth's volcanism, it still has important limitations. The magnitude and short-term variability of lower-atmosphere emissions, historically less accessible from space, remain largely uncertain. Operational monitoring of volcanic plumes, at scales relevant for adequate surveillance, has been facilitated through the use of ground-based scanning differential optical absorption spectrometer (ScanDOAS) instruments since the beginning of this century, largely due to the coordinated effort of the Network for Observation of Volcanic and Atmospheric Change (NOVAC). In this study, we present a compilation of results of homogenized post-analysis of measurements of SO2 flux and plume parameters obtained during the period March 2005 to January 2017 of 32 volcanoes in NOVAC. This inventory opens a window into the short-term emission patterns of a diverse set of volcanoes in terms of magma composition, geographical location, magnitude of emission, and style of eruptive activity. We find that passive volcanic degassing is by no means a stationary process in time and that large sub-daily variability is observed in the flux of volcanic gases, which has implications for emission budgets produced using short-term, sporadic observations. The use of a standard evaluation method allows for intercomparison between different volcanoes and between ground- and space-based measurements of the same volcanoes. The emission of several weakly degassing volcanoes, undetected by satellites, is presented for the first time. We also compare our results with those reported in the literature, providing ranges of variability in emission not accessible in the past. The open-access data repository introduced in this article will enable further exploitation of this unique dataset, with a focus on volcanological research, risk assessment, satellite-sensor validation, and improved quantification of the prevalent tropospheric component of global volcanic emission.

Description: The tropical tropopause layer (TTL) is the transition region between the well-mixed convective troposphere and the radiatively controlled stratosphere with air masses showing chemical and dynamical properties of both regions. The representation of the TTL in meteorological reanalysis data sets is important for studying the complex interactions of circulation, convection, trace gases, clouds, and radiation. In this paper, we present the evaluation of climatological and long-term TTL temperature and tropopause characteristics in the reanalysis data sets ERA-Interim, ERA5, JRA-25, JRA-55, MERRA, MERRA-2, NCEP-NCAR (R1), and CFSR. The evaluation has been performed as part of the SPARC (Stratosphere–troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The most recent atmospheric reanalysis data sets (ERA-Interim, ERA5, JRA-55, MERRA-2, and CFSR) all provide realistic representations of the major characteristics of the temperature structure within the TTL. There is good agreement between reanalysis estimates of tropical mean temperatures and radio occultation data, with relatively small cold biases for most data sets. Temperatures at the cold point and lapse rate tropopause levels, on the other hand, show warm biases in reanalyses when compared to observations. This tropopause-level warm bias is related to the vertical resolution of the reanalysis data, with the smallest bias found for data sets with the highest vertical resolution around the tropopause. Differences in the cold point temperature maximize over equatorial Africa, related to Kelvin wave activity and associated disturbances in TTL temperatures. Interannual variability in reanalysis temperatures is best constrained in the upper TTL, with larger differences at levels below the cold point. The reanalyses reproduce the temperature responses to major dynamical and radiative signals such as volcanic eruptions and the quasi-biennial oscillation (QBO). Long-term reanalysis trends in temperature in the upper TTL show good agreement with trends derived from adjusted radiosonde data sets indicating significant stratospheric cooling of around −0.5 to −1 K per decade. At 100 hPa and the cold point, most of the reanalyses suggest small but significant cooling trends of −0.3 to −0.6 K per decade that are statistically consistent with trends based on the adjusted radiosonde data sets. Advances of the reanalysis and observational systems over the last decades have led to a clear improvement in the TTL reanalysis products over time. Biases of the temperature profiles and differences in interannual variability clearly decreased in 2006, when densely sampled radio occultation data started being assimilated by the reanalyses. While there is an overall good agreement, different reanalyses offer different advantages in the TTL such as realistic profile and cold point temperature, continuous time series, or a realistic representation of signals of interannual variability. Their use in model simulations and in comparisons with climate model output should be tailored to their specific strengths and weaknesses.

Description: Oxygen minimum zones (OMZs) are characterized by enhanced carbon dioxide (CO2) levels and low pH and are being further acidified by uptake of anthropogenic atmospheric CO2. With ongoing intensification and expansion of OMZs due to global warming, carbonate chemistry conditions may become more variable and extreme, particularly in the eastern boundary upwelling systems. In austral summer (February–April) 2017, a large-scale mesocosm experiment was conducted in the coastal upwelling area off Callao (Peru) to investigate the impacts of ongoing ocean deoxygenation on biogeochemical processes, coinciding with a rare coastal El Niño event. Here we report on the temporal dynamics of carbonate chemistry in the mesocosms and surrounding Pacific waters over a continuous period of 50 d with high-temporal-resolution observations (every second day). The mesocosm experiment simulated an upwelling event in the mesocosms by addition of nitrogen (N)-deficient and CO2-enriched OMZ water. Surface water in the mesocosms was acidified by the OMZ water addition, with pHT lowered by 0.1–0.2 and pCO2 elevated to above 900 µatm. Thereafter, surface pCO2 quickly dropped to near or below the atmospheric level (405.22 µatm in 2017; Dlugokencky and Tans, 2021; NOAA/Global Monitoring Laboratory (GML)) mainly due to enhanced phytoplankton production with rapid CO2 consumption. Further observations revealed that the dominance of the dinoflagellate Akashiwo sanguinea and contamination of bird excrements played important roles in the dynamics of carbonate chemistry in the mesocosms. Compared to the simulated upwelling, natural upwelling events in the surrounding Pacific waters occurred more frequently with sea-to-air CO2 fluxes of 4.2–14.0 mmol C m−2 d−1. The positive CO2 fluxes indicated our site was a local CO2 source during our study, which may have been impacted by the coastal El Niño. However, our observations of dissolved inorganic carbon (DIC) drawdown in the mesocosms suggest that CO2 fluxes to the atmosphere can be largely dampened by biological processes. Overall, our study characterized carbonate chemistry in nearshore Pacific waters that are rarely sampled in such a temporal resolution and hence provided unique insights into the CO2 dynamics during a rare coastal El Niño event.

Description: Emissions of halogenated very short lived substances (VSLS) from the tropical oceans contribute to the atmospheric halogen budget and affect tropospheric and stratospheric ozone. Here we investigate the contribution of natural oceanic VSLS emissions to the Marine Atmospheric Boundary Layer (MABL) and their transport into the Free Troposphere (FT) over the tropical West Pacific. The study concentrates in particular on ship and aircraft measurements of the VSLS bromoform, dibromomethane and methyl iodide and meteorological parameters during the SHIVA (Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere) campaign in the South China and Sulu Seas in November 2011. Elevated oceanic concentrations of 19.9 (2.80–136.91) pmol L−1 for bromoform, 5.0 (2.43–21.82) pmol L−1 for dibromomethane and 3.8 (0.55–18.83) pmol L−1 for methyl iodide in particular close to Singapore and at the coast of Borneo with high corresponding oceanic emissions of 1486 ± 1718 pmol m−2 h−1 for bromoform, 405 ± 349 pmol m−2 h−1 for dibromomethane and 433 ± 482 pmol m−2 h−1 for methyl iodide characterize this tropical region as a strong source of these compounds. Unexpectedly atmospheric mixing ratios in the MABL were relatively low with 2.08 ± 2.08 ppt for bromoform, 1.17 ± 1.17 ppt for dibromomethane and 0.39 ± 0.09 ppt for methyl iodide. We use meteorological and chemical ship and aircraft observations, FLEXPART trajectory calculations and source-loss estimates to identify the oceanic VSLS contribution to the MABL and to the FT. Our results show that a convective, well-ventilated MABL and intense convection led to the low atmospheric mixing ratios in the MABL despite the high oceanic emissions in coastal areas of the South-China and Sulu Seas. While the accumulated bromoform in the FT above the region origins almost entirely from the local South China Sea area, dibromomethane is largely advected from distant source regions. The accumulated FT mixing ratio of methyl iodide is higher than can be explained with the local oceanic or MABL contributions. Possible reasons, uncertainties and consequences of our observations and model estimates are discussed.

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

Description: The southern African climate is strongly impacted by climate change. Precipitation is a key variable in this region, as it is linked to agriculture and water supply. Simulations with a regional atmospheric model over the past decades and the 21st century display a decrease in the past precipitation over some coastal areas of South Africa and an increase over the rest of southern Africa. However, precipitation is projected to decrease over the whole southern part of the domain in the future. This study shows that the Agulhas Current system, including the current and the leakage, which surrounds the continent in the east and south, impacts this precipitation trend. A reduction in the strength of the Agulhas Current is linked to a reduction in precipitation along the southeast coast. The Agulhas leakage, the part of the Agulhas Current that leaves the system and flows into the South Atlantic, impacts winter precipitation in the southwest of the continent. A more intense Agulhas leakage is linked to a reduction in precipitation in this region.

Description: In remote sensing applications, such as differential optical absorption spectroscopy (DOAS), atmospheric scattering processes need to be considered. After inelastic scattering on N2 and O2 molecules, the scattered photons occur as additional intensity at a different wavelength, effectively leading to filling-in of both solar Fraunhofer lines and absorptions of atmospheric constituents. Measured spectra in passive DOAS applications are typically corrected for rotational Raman scattering (RRS), also called Ring effect, which represents the main contribution to inelastic scattering. In contrast to that, vibrational Raman scattering (VRS) of N2 and O2 has often been thought to be negligible, but also contributes. Consequences of VRS are red-shifted Fraunhofer structures in scattered light spectra and filling-in of Fraunhofer lines, additional to RRS. We describe how to calculate VRS correction spectra in analogy to the Ring spectrum. We discuss further the impact of VRS cross-sections for O2 and N2 on passive DOAS measurements. The relevance of VRS is shown for the first time in spectral evaluations of Multi-Axis DOAS data. This measurement data yields in agreement with calculated scattering cross-sections, that the observed VRS cross-section amounts to 2.2 ± 0.4% of the cross-section of RRS under tropospheric conditions. It is concluded, that this phenomenon has to be included in the spectral evaluation of weak absorbers as it reduces the measurement error significantly and can cause apparent differential optical depth of up to 2.5 × 10−4. Its influence on the spectral retrieval of IO, Glyoxal, water vapour and NO2 in the blue wavelength range is evaluated. For measurements with a large Ring signal a significant and systematic bias of NO2 dSCDs up to (−3.8 ± 0.4) × 1014 molec cm−2 at low elevation angles is observed if this effect is not considered.