ALBERT

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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.

Description: Recent observational and modeling studies suggest that stratospheric ozone depletion not only influences the surface climate in the Southern Hemisphere (SH), but also impacts Northern Hemisphere (NH) spring, which implies a strong interaction between dynamics and chemistry. Here, we systematically analyze the importance of interactive chemistry with respect to the representation of stratosphere–troposphere coupling and in particular the effects on NH surface climate during the recent past. We use the interactive and specified chemistry version of NCAR's Whole Atmosphere Community Climate Model coupled to an ocean model to investigate differences in the mean state of the NH stratosphere as well as in stratospheric extreme events, namely sudden stratospheric warmings (SSWs), and their surface impacts. To be able to focus on differences that arise from two-way interactions between chemistry and dynamics in the model, the specified chemistry model version uses a time-evolving, model-consistent ozone field generated by the interactive chemistry model version. We also test the effects of zonally symmetric versus asymmetric prescribed ozone, evaluating the importance of ozone waves in the representation of stratospheric mean state and variability. The interactive chemistry simulation is characterized by a significantly stronger and colder polar night jet (PNJ) during spring when ozone depletion becomes important. We identify a negative feedback between lower stratospheric ozone and atmospheric dynamics during the breakdown of the stratospheric polar vortex in the NH, which contributes to the different characteristics of the PNJ between the simulations. Not only the mean state, but also stratospheric variability is better represented in the interactive chemistry simulation, which shows a more realistic distribution of SSWs as well as a more persistent surface impact afterwards compared with the simulation where the feedback between chemistry and dynamics is switched off. We hypothesize that this is also related to the feedback between ozone and dynamics via the intrusion of ozone-rich air into polar latitudes during SSWs. The results from the zonally asymmetric ozone simulation are closer to the interactive chemistry simulations, implying that under a model-consistent ozone forcing, a three-dimensional (3-D) representation of the prescribed ozone field is desirable. This suggests that a 3-D ozone forcing, as recommended for the upcoming CMIP6 simulations, has the potential to improve the representation of stratospheric dynamics and chemistry. Our findings underline the importance of the representation of interactive chemistry and its feedback on the stratospheric mean state and variability not only in the SH but also in the NH during the recent past.

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

Description: Oceanic emissions of the climate-relevant trace gases carbonyl sulfide (OCS) and carbon disulfide (CS2) are a major source to their atmospheric budget. Their current and future emission estimates are still uncertain due to incomplete process understanding and therefore inexact quantification across different biogeochemical regimes. Here we present the first concurrent measurements of both gases together with related fractions of the dissolved organic matter (DOM) pool, i.e., solid-phase extractable dissolved organic sulfur (DOSSPE, n=24, 0.16±0.04 µmol L−1), chromophoric (CDOM, n=76, 0.152±0.03), and fluorescent dissolved organic matter (FDOM, n=35), from the Peruvian upwelling region (Guayaquil, Ecuador to Antofagasta, Chile, October 2015). OCS was measured continuously with an equilibrator connected to an off-axis integrated cavity output spectrometer at the surface (29.8±19.8 pmol L−1) and at four profiles ranging down to 136 m. CS2 was measured at the surface (n=143, 17.8±9.0 pmol L−1) and below, ranging down to 1000 m (24 profiles). These observations were used to estimate in situ production rates and identify their drivers. We find different limiting factors of marine photoproduction: while OCS production is limited by the humic-like DOM fraction that can act as a photosensitizer, high CS2 production coincides with high DOSSPE concentration. Quantifying OCS photoproduction using a specific humic-like FDOM component as proxy, together with an updated parameterization for dark production, improves agreement with observations in a 1-D biogeochemical model. Our results will help to better predict oceanic concentrations and emissions of both gases on regional and, potentially, global scales.

Description: Eddy covariance measurements show gas transfer velocity suppression at medium to high wind speed. A wind-wave interaction described by the transformed Reynolds number is used to characterize environmental conditions favoring this suppression. We take the transformed Reynolds number parameterization to review the two most cited wind speed gas transfer velocity parameterizations: Nightingale et al. (2000) and Wanninkhof (1992, 2014). We propose an algorithm to adjust k values for the effect of gas transfer suppression and validate it with two directly measured dimethyl sulfide (DMS) gas transfer velocity data sets that experienced gas transfer suppression. We also show that the data set used in the Nightingale 2000 parameterization experienced gas transfer suppression. A compensation of the suppression effect leads to an average increase of 22% in the k vs. u relationship. Performing the same correction for Wanninkhof 2014 leads to an increase of 9.85 %. Additionally, we applied our gas transfer suppression algorithm to global air-sea flux climatologies of CO2 and DMS. The global application of gas transfer suppression leads to a decrease of 11% in DMS outgassing. We expect the magnitude of Reynolds suppression on any global air-sea gas exchange to be about 10%

Description: The Stratospheric Sulfur and its Role in Climate (SSiRC) interactive stratospheric aerosol model intercomparison project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulphur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present 4 co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The "Background" (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The "Transient Aerosol Record" (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulphur emissions and transport processes over the period 1998–2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulphate aerosol evolution after major volcanic eruptions. The "Historical Eruptions SO2 Emission Assessment" (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the "Pinatubo Emulation in Multiple models" (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt. Pinatubo eruption.

Description: The enhancement of the stratospheric aerosol layer by volcanic eruptions induces a complex set of responses causing global and regional climate effects on a broad range of timescales. Uncertainties exist regarding the climatic response to strong volcanic forcing identified in coupled climate simulations that contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). In order to better understand the sources of these model diversities, the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) has defined a coordinated set of idealized volcanic perturbation experiments to be carried out in alignment with the CMIP6 protocol. VolMIP provides a common stratospheric aerosol data set for each experiment to minimize differences in the applied volcanic forcing. It defines a set of initial conditions to assess how internal climate variability contributes to determining the response. VolMIP will assess to what extent volcanically forced responses of the coupled ocean–atmosphere system are robustly simulated by state-of-the-art coupled climate models and identify the causes that limit robust simulated behavior, especially differences in the treatment of physical processes. This paper illustrates the design of the idealized volcanic perturbation experiments in the VolMIP protocol and describes the common aerosol forcing input data sets to be used.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: Changes in incoming solar ultraviolet radiation over the 11-year solar cycle affect stratospheric ozone abundances. It is important to quantify the magnitude, structure, and seasonality of the associated solar-ozone response (SOR) to understand the impact of the 11-year solar cycle on climate. Part 1 of this two-part study uses multiple linear regression analysis to extract the SOR in a number of recently updated satellite ozone datasets covering different periods within the epoch 1970 to 2013. The annual mean SOR in the updated version 7.0 (v7.0) Stratospheric Aerosol and Gas Experiment (SAGE) II number density dataset (1984–2004) is very consistent with that found in the previous v6.2. In contrast, we find a substantial decrease in the magnitude of the SOR in the tropical upper stratosphere in the SAGE II v7.0 mixing ratio dataset (∼ 1 %) compared to the v6.2 (∼ 4 %). This difference is shown to be largely attributable to the change in the independent stratospheric temperature dataset used to convert SAGE II ozone number densities to mixing ratios. Since these temperature records contain substantial uncertainties, we suggest that datasets based on SAGE II number densities are currently most reliable for evaluating the SOR. We further analyse three extended ozone datasets that combine SAGE II v7.0 number densities with more recent GOMOS (Global Ozone Monitoring by Occultation of Stars) or OSIRIS (Optical Spectrograph and Infrared Imager System) measurements. The extended SAGE–OSIRIS dataset (1984–2013) shows a smaller and less statistically significant SOR across much of the tropical upper stratosphere compared to the SAGE II data alone. In contrast, the two SAGE–GOMOS datasets (1984–2011) show SORs that are in closer agreement with the original SAGE II data and therefore appear to provide a more reliable estimate of the SOR. We also analyse the SOR in the recent Solar Backscatter Ultraviolet Instrument (SBUV) Merged Ozone Dataset (SBUVMOD) version 8.6 (VN8.6) (1970–2012) and SBUV Merged Cohesive VN8.6 (1978–2012) datasets and compare them to the previous SBUVMOD VN8.0 (1970–2009). Over their full lengths, the three records generally agree in terms of the broad magnitude and structure of the annual mean SOR. The main difference is that SBUVMOD VN8.6 shows a smaller and less significant SOR in the tropical upper stratosphere and therefore more closely resembles the SAGE II v7.0 mixing ratio data than does the SBUV Merged Cohesive VN8.6, which has a more continuous SOR of ∼ 2 % in this region. The sparse spatial and temporal sampling of limb satellite instruments prohibits the extraction of sub-annual variations in the SOR from SAGE-based datasets. However, the SBUVMOD VN8.6 dataset suggests substantial month-to-month variations in the SOR, particularly in the winter extratropics, which may be important for the proposed high-latitude dynamical response to the solar cycle. Overall, the results highlight substantial uncertainties in the magnitude and structure of the observed SOR from different satellite records. The implications of these uncertainties for understanding and modelling the effects of solar variability on climate should be explored.

Description: A quality assessment of the CFC-11 (CCl3F), CFC-12 (CCl2F2), HF, and SF6 products from limb-viewing satellite instruments is provided by means of a detailed intercomparison. The climatologies in the form of monthly zonal mean time series are obtained from HALOE, MIPAS, ACE-FTS, and HIRDLS within the time period 1991–2010. The intercomparisons focus on the mean biases of the monthly and annual zonal mean fields and aim to identify their vertical, latitudinal and temporal structure. The CFC evaluations (based on MIPAS, ACE-FTS and HIRDLS) reveal that the uncertainty in our knowledge of the atmospheric CFC-11 and CFC-12 mean state, as given by satellite data sets, is smallest in the tropics and mid-latitudes at altitudes below 50 and 20 hPa, respectively, with a 1σ multi-instrument spread of up to ±5 %. For HF, the situation is reversed. The two available data sets (HALOE and ACE-FTS) agree well above 100 hPa, with a spread in this region of ±5 to ±10 %, while at altitudes below 100 hPa the HF annual mean state is less well known, with a spread ±30 % and larger. The atmospheric SF6 annual mean states derived from two satellite data sets (MIPAS and ACE-FTS) show only very small differences with a spread of less than ±5 % and often below ±2.5 %. While the overall agreement among the climatological data sets is very good for large parts of the upper troposphere and lower stratosphere (CFCs, SF6) or middle stratosphere (HF), individual discrepancies have been identified. Pronounced deviations between the instrument climatologies exist for particular atmospheric regions which differ from gas to gas. Notable features are differently shaped isopleths in the subtropics, deviations in the vertical gradients in the lower stratosphere and in the meridional gradients in the upper troposphere, and inconsistencies in the seasonal cycle. Additionally, long-term drifts between the instruments have been identified for the CFC-11 and CFC-12 time series. The evaluations as a whole provide guidance on what data sets are the most reliable for applications such as studies of atmospheric transport and variability, model–measurement comparisons and detection of long-term trends. The data sets will be publicly available from the SPARC Data Centre and through PANGAEA (doi:10.1594/PANGAEA.849223).

Description: A new numerical model for the on-demand computation of optimal ship routes based on sea-state forecasts has been developed. The model, named VISIR (discoVerIng Safe and effIcient Routes) is designed to support decision-makers when planning a marine voyage. The first version of the system, VISIR-I, considers medium and small motor vessels with lengths of up to a few tens of metres and a displacement hull. The model is comprised of three components: a route optimization algorithm, a mechanical model of the ship, and a processor of the environmental fields. The optimization algorithm is based on a graph-search method with time-dependent edge weights. The algorithm is also able to compute a voluntary ship speed reduction. The ship model accounts for calm water and added wave resistance by making use of just the principal particulars of the vessel as input parameters. It also checks the optimal route for parametric roll, pure loss of stability, and surfriding/broaching-to hazard conditions. The processor of the environmental fields employs significant wave height, wave spectrum peak period, and wave direction forecast fields as input. The topological issues of coastal navigation (islands, peninsulas, narrow passages) are addressed. Examples of VISIR-I routes in the Mediterranean Sea are provided. The optimal route may be longer in terms of miles sailed and yet it is faster and safer than the geodetic route between the same departure and arrival locations. Time savings up to 2.7 % and route lengthening up to 3.2 % are found for the case studies analysed. However, there is no upper bound for the magnitude of the changes of such route metrics, which especially in case of extreme sea states can be much greater. Route diversions result from the safety constraints and the fact that the algorithm takes into account the full temporal evolution and spatial variability of the environmental fields.

Description: The injection of sulphur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulphate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulphur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The VSSI estimates incorporate improvements to the ice core records in terms of synchronization and dating, refinements to the methods used to estimate VSSI from ice core records, and includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815) and Laki (1783) are within 10% of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682 and 1108 CE. The long term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr−1, ∼ 50 % greater than a prior reconstruction, due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term, latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33%) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates and including a constant "background" injection of stratospheric sulphur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2〉.

Description: Local coupling between nitrogen fixation and denitrification in current biogeochemical models could result in runaway feedback in open-ocean oxygen minimum zones (OMZs), eventually stripping OMZ waters of all fixed nitrogen. This feedback does not seem to operate at full strength in the ocean, as nitrate does not generally become depleted in open-ocean OMZs. To explore in detail the possible mechanisms that prevent nitrogen depletion in the OMZ of the eastern tropical South Pacific (ETSP), we develop a box model with fully prognostic cycles of carbon, nutrients and oxygen in the upwelling region and its adjacent open ocean. Ocean circulation is calibrated with Δ14C data of the ETSP. The sensitivity of the simulated nitrogen cycle to nutrient and oxygen exchange and ventilation from outside the model domain and to remineralization scales inside an OMZ is analysed. For the entire range of model configurations explored, we find that the fixed-N inventory can be stabilized at non-zero levels in the ETSP OMZ only if the remineralization rate via denitrification is slower than that via aerobic respiration. In our optimum model configuration, lateral oxygen supply into the model domain is required at rates sufficient to oxidize at least about one fifth of the export production in the model domain to prevent anoxia in the deep ocean. Under these conditions, our model is in line with the view of phosphate as the ultimate limiting nutrient for phytoplankton, and implies that for the current notion of nitrogen fixation being favoured in N-deficit waters, the water column of the ETSP could even be a small net source of nitrate.