Search

International audience; We examine differences among reanalysis high-cloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifth-generation European Centre for Medium-range Weather Forecasts-ECMWF-reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanaly-sis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Cli-mate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanaly-ses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of long-wave, shortwave , and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and shortwave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radia-tive heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERA-Interim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and Published by Copernicus Publications on behalf of the European Geosciences Union. 8990 J. S. Wright et al.: Tropical high clouds in reanalyses feed back to convective behaviour through their relationships with thermodynamic stability.

The Los Chocoyos (14.6°N, 91.2°W) supereruption happened ?75,000 years ago in Guatemala and was one of the largest eruptions of the past 100,000 years. It emitted enormous amounts of sulfur, chlorine, and bromine, with multi-decadal consequences for the global climate and environment. Here, we simulate the impact of a Los Chocoyos-like eruption on the quasi-biennial oscillation (QBO), an oscillation of zonal winds in the tropical stratosphere, with a comprehensive aerosol chemistry Earth System Model. We find a ?10-year disruption of the QBO starting 4 months post eruption, with anomalous easterly winds lasting ?5 years, followed by westerlies, before returning to QBO conditions with a slightly prolonged periodicity. Volcanic aerosol heating and ozone depletion cooling leads to the QBO disruption and anomalous wind regimes through radiative changes and wave-mean flow interactions. Different model ensembles, volcanic forcing scenarios and results of a second model back up the robustness of our results.

The supereruption of Los Chocoyos (14.6∘ N, 91.2∘ W) in Guatemala ∼84 kyr ago was one of the largest volcanic events of the past 100 000 years. Recent petrologic data show that the eruption released very large amounts of climate-relevant sulfur and ozone-destroying chlorine and bromine gases (523±94 Mt sulfur, 1200±156 Mt chlorine, and 2±0.46 Mt bromine). Using the Earth system model (ESM) of the Community Earth System Model version 2 (CESM2) coupled with the Whole Atmosphere Community Climate Model version 6 (WACCM6), we simulated the impacts of the sulfur- and halogen-rich Los Chocoyos eruption on the preindustrial Earth system. Our simulations show that elevated sulfate burden and aerosol optical depth (AOD) persists for 5 years in the model, while the volcanic halogens stay elevated for nearly 15 years. As a consequence, the eruption leads to a collapse of the ozone layer with global mean column ozone values dropping to 50 DU (80 % decrease) and leading to a 550 % increase in surface UV over the first 5 years, with potential impacts on the biosphere. The volcanic eruption shows an asymmetric-hemispheric response with enhanced aerosol, ozone, UV, and climate signals over the Northern Hemisphere. Surface climate is impacted globally due to peak AOD of >6, which leads to a maximum surface cooling of >6 K, precipitation and terrestrial net primary production decrease of >25 %, and sea ice area increases of 40 % in the first 3 years. Locally, a wetting (>100 %) and strong increase in net primary production (NPP) (>700 %) over northern Africa is simulated in the first 5 years and related to a southward shift of the Intertropical Convergence Zone (ITCZ) to the southern tropics. The ocean responds with pronounced El Niño conditions in the first 3 years that shift to the southern tropics and are coherent with the ITCZ change. Recovery to pre-eruption ozone levels and climate takes 15 years and 30 years, respectively. The long-lasting surface cooling is sustained by an immediate increase in the Arctic sea ice area, followed by a decrease in poleward ocean heat transport at 60∘ N which lasts up to 20 years. In contrast, when simulating Los Chocoyos conventionally by including sulfur and neglecting halogens, we simulate a larger sulfate burden and AOD, more pronounced surface climate changes, and an increase in column ozone. By comparing our aerosol chemistry ESM results to other supereruption simulations with aerosol climate models, we find a higher surface climate impact per injected sulfur amount than previous studies for our different sets of model experiments, since the CESM2(WACCM6) creates smaller aerosols with a longer lifetime, partly due to the interactive aerosol chemistry. As the model uncertainties for the climate response to supereruptions are very large, observational evidence from paleo archives and a coordinated model intercomparison would help to improve our understanding of the climate and environment response.

Uncertainty in the physicochemical and optical properties of volcanic ash particles creates errors in the detection and modeling of volcanic ash clouds and in quantification of their potential impacts. In this study, we provide a dataset that describes the physicochemical and optical properties of a representative selection of volcanic ash samples from nine different volcanic eruptions covering a wide range of silica contents (50 - 80 wt. % SiO2). We measured and calculated parameters describing the physical (size distribution, complex shape, and dense-rock equivalent mass density), chemical (bulk and surface composition) and optical (complex refractive index from ultra-violet to near-infrared wavelengths) properties of the volcanic ash and classified the samples according to their SiO2 and total alkali contents into the common igneous rock types Basalt to Rhyolite. We found that the mass density ranges between ρ = 2.49 and 2.98 g/cm3 for rhyolitic to basaltic ash types, and that the particle shape varies with changing particle size (d < 100 μm). The complex refractive indices in the wavelength range between λ = 300 nm and 1500 nm depends systematically on the composition of the samples. The real part values vary from n = 1.38 to 1.66 depending on ash type and wavelength and the imaginary part values from k = 0.00027 to 0.00268. We place our results into the context of existing data, and thus provide a comprehensive dataset that can be used for future and historic eruptions, when only basic information about the magma type producing the ash is known.

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.

Large explosive tropical volcanic eruptions inject high amounts of gases into the stratosphere, where they disperse globally through the large-scale meridional circulation. There is now increasing observational evidence that volcanic halogens can reach the upper troposphere and lower stratosphere. Here, we present the first study that combines measurement-based data of sulfur, chlorine and bromine releases from tropical volcanic eruptions with complex coupled chemistry climate model simulations taking radiative-dynamical-chemical feedbacks into account. Halogen model input parameters represent a size-time-region-wide average for the Central American eruptions over the last 200 ka ensuring a comprehensive perspective. The simulations reveal global, long-lasting impact on the ozone layer affecting atmospheric composition and circulation for a decade. Column ozone drops below 220 DU (ozone hole conditions) in the tropics, Arctic and Antarctica, increasing biologically active UV by 80 to 400%. Our model results could potentially be validated using high-resolution proxies from ice cores and pollen records.

Halogenated Very Short-lived Substances (VSLS), such as bromoform, dibromomethane and methyl iodide, are naturally produced in the oceans and are involved in ozone depletion in the troposphere and the stratosphere. The effect of climate change on the oceanic emissions of these compounds is not well quantified. Based on present-day observed global oceanic and atmospheric concentrations, and historic and future data from three CMIP5 models, past and future sea-to-air fluxes of these VSLS are calculated. The simulations are used to infer possible effects of projected changes of physical forcing on emissions in different oceanic regimes. CMIP5 model output for 1979–2100 from the historical scenario and the RCP scenarios 2.6 and 8.5 are used as input data for the emission calculations. Of the parameters that have the main influence on the sea-to-air fluxes, the global sea surface temperatures show a steady increase during the twenty-first century, while the projected changes of sea surface wind speed is very small. The calculated emissions based on the historical CMIP5 model runs (1979–2005) increased over the 26 year period and agree well with the emissions based on ERA-Interim data. The future sea-to-air fluxes of VSLS generally increase during the twenty-first century under the assumption of constant concentration fields in the ocean and atmosphere. The multi-model mean global emissions of bromoform increase by 29.4% (9.0%) between 1986 and 2005 and 2081–2100 under RCP 8.5 (2.6) and dibromomethane and methyl iodide emissions increase by 23.3% (6.4%) and 5.5% (1.5%), respectively. Uncertainties of the future emission estimates, driven by ongoing environmental changes such as changing oceanic productivity (not considered in this study) are discussed.

The current generation of Earth system models that participate in the Coupled Model Intercomparison Project phase 5 (CMIP5) does not, on average, produce a strengthened Northern Hemisphere (NH) polar vortex after large tropical volcanic eruptions as suggested by observational records. Here we investigate the impact of volcanic eruptions on the NH winter stratosphere with an ensemble of 20 model simulations of the Max Planck Institute Earth system model. We compare the dynamical impact in simulations of the very large 1815 Tambora eruption with the averaged dynamical response to the two largest eruptions of the CMIP5 historical simulations (the 1883 Krakatau and the 1991 Pinatubo eruptions). We find that for both the Tambora and the averaged Krakatau-Pinatubo eruptions the radiative perturbation only weakly affects the polar vortex directly. The position of the maximum temperature anomaly gradient is located at approximately 30°N, where we obtain significant westerly zonal wind anomalies between 10hPa and 30hPa. Under the very strong forcing of the Tambora eruption, the NH polar vortex is significantly strengthened because the subtropical westerly wind anomalies are sufficiently strong to robustly alter the propagation of planetary waves. The average response to the eruptions of Krakatau and Pinatubo reveals a slight strengthening of the polar vortex, but individual ensemble members differ substantially, indicating that internal variability plays a dominant role. For the Tambora eruption the ensemble variability of the zonal mean temperature and zonal wind anomalies during midwinter and late winter is significantly reduced compared to the volcanically unperturbed period.

Halogen- and sulfur-containing compounds are supersaturated in the surface ocean, which results in their emission to the atmosphere. These compounds can be transported to the stratosphere, where they impact ozone, the background aerosol layer, and climate. In this study we calculate the seasonal and interannual variability of transport from the West Indian Ocean (WIO) surface to the stratosphere for 2000–2016 with the Lagrangian transport model FLEXPART using ERA-Interim meteorological fields. We investigate the transport relevant for very short lived substances (VSLS) with tropospheric lifetimes corresponding to dimethylsulfide (1 day), methyl iodide (CH3I, 3.5 days), bromoform (CHBr3, 17 days), and dibromomethane (CH2Br2, 150 days). The stratospheric source gas injection of VSLS tracers from the WIO shows a distinct annual cycle associated with the Asian monsoon. Over the 16-year time series, a slight increase in source gas injection from the WIO to the stratosphere is found for all VSLS tracers and during all seasons. The interannual variability shows a relationship with sea surface temperatures in the WIO as well as the El Niño–Southern Oscillation. During boreal spring of El Niño, enhanced stratospheric injection of VSLS from the tropical WIO is caused by positive sea surface temperature anomalies and enhanced vertical uplift above the WIO. During boreal fall of La Niña, strong injection is related to enhanced atmospheric upward motion over the East Indian Ocean and a prolonged Indian summer monsoon season. Related physical mechanisms and uncertainties are discussed in this study. ©2018. American Geophysical Union. All Rights Reserved.

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.

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.

Assessments of climate sensitivity to projected greenhouse gas concentrations underpin environmental policy decisions, with such assessments often based on model simulations of climate during recent centuries and millennia1, 2, 3. These simulations depend critically on accurate records of past aerosol forcing from global-scale volcanic eruptions, reconstructed from measurements of sulphate deposition in ice cores4, 5, 6. Non-uniform transport and deposition of volcanic fallout mean that multiple records from a wide array of ice cores must be combined to create accurate reconstructions. Here we re-evaluated the record of volcanic sulphate deposition using a much more extensive array of Antarctic ice cores. In our new reconstruction, many additional records have been added and dating of previously published records corrected through precise synchronization to the annually dated West Antarctic Ice Sheet Divide ice core7, improving and extending the record throughout the Common Era. Whereas agreement with existing reconstructions is excellent after 1500, we found a substantially different history of volcanic aerosol deposition before 1500; for example, global aerosol forcing values from some of the largest eruptions (for example, 1257 and 1458) previously were overestimated by 20–30% and others underestimated by 20–50%.

We investigate the contribution of oceanic methyl iodide (CH3I) to the stratospheric iodine budget. Based on CH3I measurements from three tropical ship campaigns and the Lagrangian transport model FLEXPART, we provide a detailed analysis of CH3I transport from the ocean surface to the cold point in the upper tropical tropopause layer (TTL). While average oceanic emissions differ by less than 50% from campaign to campaign, the measurements show much stronger variations within each campaign. A positive correlation between the oceanic CH3I emissions and the efficiency of CH3I troposphere–stratosphere transport has been identified for some cruise sections. The mechanism of strong horizontal surface winds triggering large emissions on the one hand and being associated with tropical convective systems, such as developing typhoons, on the other hand, could explain the identified correlations. As a result of the simultaneous occurrence of large CH3I emissions and strong vertical uplift, localized maximum mixing ratios of 0.6 ppt CH3I at the cold point have been determined for observed peak emissions during the SHIVA (Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere)-Sonne research vessel campaign in the coastal western Pacific. The other two campaigns give considerably smaller maxima of 0.1 ppt CH3I in the open western Pacific and 0.03 ppt in the coastal eastern Atlantic. In order to assess the representativeness of the large local mixing ratios, we use climatological emission scenarios to derive global upper air estimates of CH3I abundances. The model results are compared with available upper air measurements, including data from the recent ATTREX and HIPPO2 aircraft campaigns. In the eastern Pacific region, the location of the available measurement campaigns in the upper TTL, the comparisons give a good agreement, indicating that around 0.01 to 0.02 ppt of CH3I enter the stratosphere. However, other tropical regions that are subject to stronger convective activity show larger CH3I entrainment, e.g., 0.08 ppt in the western Pacific. Overall our model results give a tropical contribution of 0.04 ppt CH3I to the stratospheric iodine budget. The strong variations in the geographical distribution of CH3I entrainment suggest that currently available upper air measurements are not representative of global estimates and further campaigns will be necessary in order to better understand the CH3I contribution to stratospheric iodine.

Sea surface and atmospheric measurements of dimethylsulphide (DMS) were performed during the TransBrom cruise in the western Pacific Ocean between Japan and Australia in October 2009. Air–sea DMS fluxes were computed between 0 and 30 μmol m−2 d−1, which are in agreement with those computed by the current climatology, and peak emissions of marine DMS into the atmosphere were found during the occurrence of tropical storm systems. Atmospheric variability in DMS, however, did not follow that of the computed fluxes and was more related to atmospheric transport processes. The computed emissions were used as input fields for the Lagrangian dispersion model FLEXPART, which was set up with actual meteorological fields from ERA-Interim data and different chemical lifetimes of DMS. A comparison with aircraft in situ data from the adjacent HIPPO2 campaign revealed an overall good agreement between modelled versus observed DMS profiles over the tropical western Pacific Ocean. Based on observed DMS emissions and meteorological fields along the cruise track, the model projected that up to 30 g S per month in the form of DMS, emitted from an area of 6 × 104 m2, can be transported above 17 km. This surprisingly large DMS entrainment into the stratosphere is disproportionate to the regional extent of the area of emissions and mainly due to the high convective activity in this region as simulated by the transport model. Thus, if DMS can cross the tropical tropopause layer (TTL), we suggest that the considerably larger area of the tropical western Pacific Ocean can be a source of sulphur to the stratosphere, which has not been considered as yet.

Reconstructions of the atmospheric sulfate aerosol burdens resulting from past volcanic eruptions are based on ice core-derived estimates of volcanic sulfate deposition and the assumption that the two quantities are directly proportional. We test this assumption within simulations of tropical volcanic stratospheric sulfur injections with the MAECHAM5-HAM aerosol-climate model. An ensemble of 70 simulations is analyzed, with SO2 injections ranging from 8.5 to 700 Tg, with eruptions in January and July. Modeled sulfate deposition flux to Antarctica shows excellent spatial correlation with ice core-derived estimates for Pinatubo and Tambora, although the comparison suggests the modeled flux to the ice sheets is 4–5 times too large. We find that Greenland and Antarctic deposition efficiencies (the ratio of sulfate flux to each ice sheet to the maximum hemispheric stratospheric sulfate aerosol burden) vary as a function of the magnitude and season of stratospheric sulfur injection. Changes in simulated sulfate deposition for large SO2 injections are connected to increases in aerosol particle size, which impact aerosol sedimentation velocity and radiative properties, the latter leading to strong dynamical changes including strengthening of the winter polar vortices, which inhibits the transport of stratospheric aerosols to high latitudes. The resulting relationship between Antarctic and Greenland volcanic sulfate deposition is nonlinear for very large eruptions, with significantly less sulfate deposition to Antarctica than to Greenland. These model results suggest that variability of deposition efficiency may be an important consideration in the interpretation of ice core sulfate signals for eruptions of Tambora-magnitude and larger.

We describe the main differences in simulations of stratospheric climate and variability by models within the fifth Coupled Model Intercomparison Project (CMIP5) that have a model top above the stratopause and relatively fine stratospheric vertical resolution (high-top), and those that have a model top below the stratopause (low-top). Although the simulation of mean stratospheric climate by the two model ensembles is similar, the low-top model ensemble has very weak stratospheric variability on daily and interannual time scales. The frequency of major sudden stratospheric warming events is strongly underestimated by the low-top models with less than half the frequency of events observed in the reanalysis data and high-top models. The lack of stratospheric variability in the low-top models affects their stratosphere-troposphere coupling, resulting in short-lived anomalies in the Northern Annular Mode, which do not produce long-lasting tropospheric impacts, as seen in observations. The lack of stratospheric variability, however, does not appear to have any impact on the ability of the low-top models to reproduce past stratospheric temperature trends. We find little improvement in the simulation of decadal variability for the high-top models compared to the low-top, which is likely related to the fact that neither ensemble produces a realistic dynamical response to volcanic eruptions.

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.

The stratospheric climate and variability from simulations of sixteen chemistry-climate models is evaluated. On average the polar night jet is well reproduced though its variability is less well reproduced with a large spread between models. Polar temperature biases are less than 5 K except in the Southern Hemisphere (SH) lower stratosphere in spring. The accumulated area of low temperatures responsible for polar stratospheric cloud formation is accurately reproduced for the Antarctic but underestimated for the Arctic. The shape and position of the polar vortex is well simulated, as is the tropical upwelling in the lower stratosphere. There is a wide model spread in the frequency of major sudden stratospheric warnings (SSWs), late biases in the breakup of the SH vortex, and a weak annual cycle in the zonal wind in the tropical upper stratosphere. Quantitatively, "metrics" indicate a wide spread in model performance for most diagnostics with systematic biases in many, and poorer performance in the SH than in the Northern Hemisphere (NH). Correlations were found in the SH between errors in the final warming, polar temperatures, the leading mode of variability, and jet strength, and in the NH between errors in polar temperatures, frequency of major SSWs, and jet strength. Models with a stronger QBO have stronger tropical upwelling and a colder NH vortex. Both the qualitative and quantitative analysis indicate a number of common and long-standing model problems, particularly related to the simulation of the SH and stratospheric variability. Copyright 2011 by the American Geophysical Union.

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.

Volcanic activity in and around the year 536 CE led to severe cold and famine, and has been speculatively linked to large-scale societal crises around the globe. Using a coupled aerosol-climate model, with eruption parameters constrained by recently re-dated ice core records and historical observations of the aerosol cloud, we reconstruct the radiative forcing resulting from a sequence of two major volcanic eruptions in 536 and 540 CE. We estimate that the decadal-scale Northern Hemisphere (NH) extra-tropical radiative forcing from this volcanic “double event” was larger than that of any period in existing reconstructions of the last 1200 years. Earth system model simulations including the volcanic forcing show peak NH mean temperature anomalies reaching more than −2 °C, and show agreement with the limited number of available maximum latewood density temperature reconstructions. The simulations also produce decadal-scale anomalies of Arctic sea ice. The simulated cooling is interpreted in terms of probable impacts on agricultural production in Europe, and implies a high likelihood of multiple years of significant decreases in crop production across Scandinavia, supporting the theory of a connection between the 536 and 540 eruptions and evidence of societal crisis dated to the mid-6th century. © 2016 The Author(s)

Simulations of tropical volcanic eruptions using a general circulation model with coupled aerosol microphysics are used to assess the influence of season of eruption on the aerosol evolution and radiative impacts at the Earth's surface. This analysis is presented for eruptions with SO2 injection magnitudes of 17 and 700 Tg, the former consistent with estimates of the 1991 Mt. Pinatubo eruption, the later a near-"super eruption". For each eruption magnitude, simulations are performed with eruptions at 15° N, at four equally spaced times of year. Sensitivity to eruption season of aerosol optical depth (AOD), clear-sky and all-sky shortwave (SW) radiative flux is quantified by first integrating each field for four years after the eruption, then calculating for each cumulative field the absolute or percent difference between the maximum and minimum response from the four eruption seasons. Eruption season has a significant influence on AOD and clear-sky SW radiative flux anomalies for both eruption magnitudes. The sensitivity to eruption season for both fields is generally weak in the tropics, but increases in the mid- and high latitudes, reaching maximum values of ~75 %. Global mean AOD and clear-sky SW anomalies show sensitivity to eruption season on the order of 15–20 %, which results from differences in aerosol effective radius for the different eruption seasons. Smallest aerosol size and largest cumulative impact result from a January eruption for Pinatubo-magnitude eruption, and from a July eruption for the near-super eruption. In contrast to AOD and clear-sky SW anomalies, all-sky SW anomalies are found to be insensitive to season of eruption for the Pinatubo-magnitude eruption experiment, due to the reflection of solar radiation by clouds in the mid- to high latitudes. However, differences in all-sky SW anomalies between eruptions in different seasons are significant for the larger eruption magnitude, and the ~15 % sensitivity to eruption season of the global mean all-sky SW anomalies is comparable to the sensitivity of global mean AOD and clear-sky SW anomalies. Our estimates of sensitivity to eruption season are larger than previously reported estimates: implications regarding volcanic AOD timeseries reconstructions and their use in climate models are discussed.

For the first time the long-term interannual and spatial variability of residence time (τ) is presented for the TTL between 360 K and 400 K potential temperature (~14 to 18 km altitude). The analysis is based on a Lagrangian approach using offline calculated diabatic heating rates as vertical velocities, covering Northern Hemisphere (NH) winters from 1962–2004. The residence time τLCP–400 K, being the duration time of air parcels in the layer between the Lagrangian Cold Point (LCP) and 400 K, varies spatially and is longer (>50 days) over the maritime continent as the LCP is lowest there (

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.

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.

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.

At present, anthropogenic halogens and oceanic emissions of very short-lived substances (VSLSs) both contribute to the observed stratospheric ozone depletion. Emissions of the long-lived anthropogenic halogens have been reduced and are currently declining, whereas emissions of the biogenic VSLSs are expected to increase in future climate due to anthropogenic activities affecting oceanic production and emissions. Here, we introduce a new approach for assessing the impact of oceanic halocarbons on stratospheric ozone by calculating their ozone depletion potential (ODP)-weighted emissions. Seasonally and spatially dependent, global distributions are derived within a case-study framework 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 halogens. 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 western Pacific Ocean. Variations of tropical convective activity lead to seasonal shifts in the spatial distribution of the trajectory-derived ODP with the updraught mass flux, used as a proxy for trajectory-derived ODP, explaining 71 % of the variance of the ODP distribution. Future climate projections based on the RCP 8.5 scenario suggest a 31 % increase of the ODP-weighted CHBr3 emissions by 2100 compared to present values. This increase is related to a larger convective updraught mass flux in the upper troposphere and increasing emissions in a future climate. However, at the same time, it is reduced by less effective bromine-related ozone depletion due to declining stratospheric chlorine concentrations. The comparison of the ODP-weighted emissions of short- and long-lived halocarbons provides a new concept for assessing the overall impact of oceanic halocarbon emissions on stratospheric ozone depletion for current conditions and future projections.

In late May 2005 unusual high levels of solar ultraviolet radiation were observed over central Europe. In Northern Germany the measured irradiance of erythemally effective radiation exceeded the climatological mean by more than about 20%. An extreme low ozone event for the season coincided with high solar elevation angles and high pressure induced clear sky conditions leading to the highest value of erythemal UV-radiation ever observed over this location in May since 1994. This hereafter called "ozone mini-hole" was caused by an elevation of tropopause height accompanied with a poleward advection of ozone-poor air from the tropics. The resultant increase in UV-radiation is of particular significance for human health. Dynamically induced low ozone episodes that happen in late spring can considerably enhance the solar UV-radiation in mid latitudes and therefore contribute to the UV-burden of people living in these regions.

The purpose of this study is to investigate horizontal transport processes in the winter stratosphere using data with a resolution relevant for chemistry and climate modeling. For this reason the Freie Universität Berlin Climate Middle Atmosphere Model (FUB-CMAM) with its model top at 83km altitude, increased horizontal resolution T42 and the semi-Lagrangian transport scheme for advecting passive tracers is used. A new approach of this paper is the classification of specific transport phenomena within the stratosphere into tropical-subtropical streamers (e.g. Offermann et al., 1999) and polar vortex extrusions hereafter called polar vortex streamers. To investigate the role played by these large-scale structures on the inter-annual and seasonal variability of transport processes in northern mid-latitudes, the global occurrence of such streamers was calculated based on a 10-year model climatology, concentrating on the existence of the Arctic polar vortex. For the identification and counting of streamers, the new method of zonal anomaly was chosen. The analysis of the months October-May yielded a maximum occurrence of tropical-subtropical streamers during Arctic winter and spring in the middle and upper stratosphere. Synoptic maps revealed highest intensities in the subtropics over East Asia with a secondary maximum over the Atlantic in the northern hemisphere. Furthermore, tropical-subtropical streamers exhibited a higher occurrence than polar vortex streamers, indicating that the subtropical barrier is more permeable than the polar vortex barrier (edge) in the model, which is in good correspondence with observations (e.g. Plumb, 2002; Neu et al., 2003). Interesting for the total ozone decrease in mid-latitudes is the consideration of the lower stratosphere for tropical-subtropical streamers and the stratosphere above ~20km altitude for polar vortex streamers, where strongest ozone depletion is observed at polar latitudes (WMO, 2003). In the lower stratosphere the FUB-CMAM simulated a climatological maximum of 10% occurrence of tropical-subtropical streamers over East-Asia/West Pacific and the Atlantic during early- and mid-winter. The results of this paper demonstrate that stratospheric streamers e.g. large-scale, tongue-like structures transporting tropical-subtropical and polar vortex air masses into mid-latitudes occur frequently during Arctic winter. They can therefore play a significant role on the strength and variability of the observed total ozone decrease at mid-latitudes and should not be neglected in future climate change studies.

Observations and simple theoretical arguments suggest that the Northern Hemisphere (NH) stratospheric polar vortex is stronger in winters following major volcanic eruptions. However, recent studies show that climate models forced by prescribed volcanic aerosol fields fail to reproduce this effect. We investigate the impact of volcanic aerosol forcing on stratospheric dynamics, including the strength of the NH polar vortex, in ensemble simulations with the Max Planck Institute Earth System Model. The model is forced by four different prescribed forcing sets representing the radiative properties of stratospheric aerosol following the 1991 eruption of Mt. Pinatubo: two forcing sets are based on observations, and are commonly used in climate model simulations, and two forcing sets are constructed based on coupled aerosol–climate model simulations. For all forcings, we find that temperature and zonal wind anomalies in the NH high latitudes are not directly impacted by anomalous volcanic aerosol heating. Instead, high latitude effects result from robust enhancements in stratospheric residual circulation, which in turn result, at least in part, from enhanced stratospheric wave activity. High latitude effects are therefore much less robust than would be expected if they were the direct result of aerosol heating. While there is significant ensemble variability in the high latitude response to each aerosol forcing set, the mean response is sensitive to the forcing set used. Significant differences, for example, are found in the NH polar stratosphere temperature and zonal wind response to two different forcing data sets constructed from different versions of SAGE II aerosol observations. Significant strengthening of the polar vortex, in rough agreement with the expected response, is achieved only using aerosol forcing extracted from prior coupled aerosol–climate model simulations. Differences in the dynamical response to the different forcing sets used imply that reproducing the polar vortex responses to past eruptions, or predicting the response to future eruptions, depends on accurate representation of the space-time structure of the volcanic aerosol forcing.

Most of the short-lived biogenic and anthropogenic chemical species that are emitted into the atmosphere break down efficiently by reaction with OH and do not reach the stratosphere. Here we show the existence of a pronounced minimum in the tropospheric column of ozone over the West Pacific, the main source region for stratospheric air, and suggest a corresponding minimum of the tropospheric column of OH. This has the potential to amplify the impact of surface emissions on the stratospheric composition compared to the impact when assuming globally uniform OH conditions. Specifically, the role of emissions of biogenic halogenated species for the stratospheric halogen budget and the role of increasing emissions of SO2 in Southeast Asia or from minor volcanic eruptions for the increasing stratospheric aerosol loading need to be reassessed in light of these findings. This is also important since climate change will further modify OH abundances and emissions of halogenated species. Our study is based on ozone sonde measurements carried out during the TransBrom cruise with the RV Sonne roughly along 140–150° E in October 2009 and corroborating ozone and OH measurements from satellites, aircraft campaigns and FTIR instruments. Model calculations with the GEOS-Chem Chemistry and Transport Model (CTM) and the ATLAS CTM are used to simulate the tropospheric OH distribution over the West Pacific and the transport pathways to the stratosphere. The potential effect of the OH minimum on species transported into the stratosphere is shown via modeling the transport and chemistry of CH2Br2 and SO2.

20 pags, figs13, 8 tabs, Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidising capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF (National Science Foundation) HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) campaign, in the tropical western Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ∼4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (∼4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics. © 2013 Author(s).f 0., This work was supported by the UK Natural Environment Research Council (NERC) and the EU Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) project (SHIVA-226224-FP7-ENV-2008-1). TOMCAT work is suported by the National Centre for Atmospheric Science (NCAS). NOAA measurements were supported in part by NOAA’s Atmospheric Chemistry, Carbon Cycle and Climate Program of its Climate Program Office. Technical assistance, standardisation, and programmatic support for NOAA flask measurements was provided by C. Siso, B. Hall, and J. W. Elkins. The HIPPO project was supported by the National Science Foundation (NSF).

Monthly zonal mean climatologies of atmospheric measurements from satellite instruments can have biases due to the non-uniform sampling of the atmosphere by the instruments. We characterize potential sampling biases in stratospheric trace gas climatologies of the Stratospheric Processes and their Role in Climate (SPARC) Data Initiative using chemical fields from a chemistry climate model simulation and sampling patterns from 16 satellite-borne instruments. The exercise is performed for the long-lived stratospheric trace gases O3 and H2O. Monthly sample biases for O3 exceed 10% for many instruments in the high latitude stratosphere and in the upper troposphere/lower stratosphere, while annual mean sampling biases reach values of up to 20% in the same regions for some instruments. Sampling biases for H2O are generally smaller than for O3, although still notable in the upper troposphere/lower stratosphere and Southern Hemisphere high latitudes. The most important mechanism leading to monthly sampling bias is the non-uniform temporal sampling of many instruments, i.e., the fact that for many instruments, monthly means are produced from measurements which span less than the full month in question. Similarly, annual mean sampling biases are well explained by non-uniformity in the month-to-month sampling by different instruments. Non-uniform sampling in latitude and longitude are shown to also lead to non-negligible sampling biases, which are most relevant for climatologies which are otherwise free of sampling biases due to non-uniform temporal sampling.

Methyl iodide (CH3I), bromoform (CHBr3) and dibromomethane (CH2Br2), which are produced naturally in the oceans, take part in ozone chemistry both in the troposphere and the stratosphere. The significance of oceanic upwelling regions for emissions of these trace gases in the global context is still uncertain although they have been identified as important source regions. To better quantify the role of upwelling areas in current and future climate, this paper analyzes major factors that influenced halocarbon emissions from the tropical North East Atlantic including the Mauritanian upwelling during the DRIVE expedition. Diel and regional variability of oceanic and atmospheric CH3I, CHBr3 and CH2Br2 was determined along with biological and physical parameters at six 24 h-stations. Low oceanic concentrations of CH3I from 0.1–5.4 pmol L−1 were equally distributed throughout the investigation area. CHBr3 and CH2Br2 from 1.0 to 42.4 pmol L−1 and to 9.4 pmol L−1, respectively were measured with maximum concentrations close to the Mauritanian coast. Atmospheric CH3I, CHBr3, and CH2Br2 of up to 3.3, 8.9, and 3.1 ppt, respectively were detected above the upwelling, as well as up to 1.8, 12.8, and 2.2 ppt at the Cape Verdean coast. While diel variability in CH3I emissions could be mainly ascribed to oceanic non-biological production, no main driver was identified for its emissions over the entire study region. In contrast, biological parameters showed the greatest influence on the regional distribution of sea-to-air fluxes of bromocarbons. The diel impact of wind speed on bromocarbon emissions increased with decreasing distance to the coast. The height of the marine atmospheric boundary layer (MABL) influenced halocarbon emissions via its influence on atmospheric mixing ratios. Oceanic and atmospheric halocarbons correlated well in the study region, and in combination with high oceanic CH3I, CHBr3 and CH2Br2 concentrations, local hot spots of atmospheric halocarbons could solely be explained by marine sources. This conclusion is in contrast to previous studies that hypothesized elevated atmospheric halocarbons above the eastern tropical Atlantic to be mainly originated from the West-African continent.

During the DRIVE~(Diurnal and Regional Variability of Halogen Emissions) ship campaign we investigated the variability of the halogenated very short-lived substances (VSLS) bromoform (CHBr3), dibromomethane (CH2Br2) and methyl iodide (CH3I) in the marine atmospheric boundary layer in the eastern tropical and subtropical North Atlantic Ocean during May/June 2010. The highest VSLS mixing ratios were found near the Mauritanian coast and close to Lisbon (Portugal). With backward trajectories we identified predominantly air masses from the open North Atlantic with some coastal influence in the Mauritanian upwelling area, due to the prevailing NW winds. The maximum VSLS mixing ratios above the Mauritanian upwelling were 8.92 ppt for bromoform, 3.14 ppt for dibromomethane and 3.29 ppt for methyl iodide, with an observed maximum range of the daily mean up to 50% for bromoform, 26% for dibromomethane and 56% for methyl iodide. The influence of various meteorological parameters – such as wind, surface air pressure, surface air and surface water temperature, humidity and marine atmospheric boundary layer (MABL) height – on VSLS concentrations and fluxes was investigated. The strongest relationship was found between the MABL height and bromoform, dibromomethane and methyl iodide abundances. Lowest MABL heights above the Mauritanian upwelling area coincide with highest VSLS mixing ratios and vice versa above the open ocean. Significant high anti-correlations confirm this relationship for the whole cruise. We conclude that especially above oceanic upwelling systems, in addition to sea–air fluxes, MABL height variations can influence atmospheric VSLS mixing ratios, occasionally leading to elevated atmospheric abundances. This may add to the postulated missing VSLS sources in the Mauritanian upwelling region (Quack et al., 2007).

We carry out a case study of the transport and chemistry of bromoform and its product gases (PGs) in a sea breeze driven convective episode on 19 November 2011 along the North West coast of Borneo during the "Stratospheric ozone: Halogen Impacts in a Varying Atmosphere" (SHIVA) campaign. We use ground based, ship, aircraft and balloon sonde observations made during the campaign, and a 3-D regional online transport and chemistry model capable of resolving clouds and convection explicitly that includes detailed bromine chemistry. The model simulates the temperature, wind speed, wind direction fairly well for the most part, and adequately captures the convection location, timing, and intensity. The simulated transport of bromoform from the boundary layer up to 12 km compares well to aircraft observations to support our conclusions. The model makes several predictions regarding bromine transport from the boundary layer to the level of convective detrainment (11 to 12 km). First, the majority of bromine undergoes this transport as bromoform. Second, insoluble organic bromine carbonyl species are transported to between 11 and 12 km, but only form a small proportion of the transported bromine. Third, soluble bromine species, which include bromine organic peroxides, hydrobromic acid (HBr), and hypobromous acid (HOBr), are washed out efficiently within the core of the convective column. Fourth, insoluble inorganic bromine species (principally Br2) are not washed out of the convective column, but are also not transported to the altitude of detrainment in large quantities. We expect that Br2 will make a larger relative contribution to the total vertical transport of bromine atoms in scenarios with higher CHBr3 mixing ratios in the boundary layer, which have been observed in other regions. Finally, given the highly detailed description of the chemistry, transport and washout of bromine compounds within our simulations, we make a series of recommendations about the physical and chemical processes that should be represented in 3-D chemical transport models (CTMs) and chemistry climate models (CCMs), which are the primary theoretical means of estimating the contribution made by CHBr3 and other very short-lived substances (VSLS) to the stratospheric bromine budget.

Volatile halogenated organic compounds containing bromine and iodine, which are naturally produced in the ocean, are involved in ozone depletion in both the troposphere and stratosphere. Three prominent compounds transporting large amounts of marine halogens into the atmosphere are bromoform (CHBr3), dibromomethane (CH2Br2) and methyl iodide (CH3I). The input of marine halogens to the stratosphere has been estimated from observations and modelling studies using low-resolution oceanic emission scenarios derived from top-down approaches. In order to improve emission inventory estimates, we calculate data-based high resolution global sea-to-air flux estimates of these compounds from surface observations within the HalOcAt (Halocarbons in the Ocean and Atmosphere) database (https://halocat.geomar.de/). Global maps of marine and atmospheric surface concentrations are derived from the data which are divided into coastal, shelf and open ocean regions. Considering physical and biogeochemical characteristics of ocean and atmosphere, the open ocean water and atmosphere data are classified into 21 regions. The available data are interpolated onto a 1°×1° grid while missing grid values are interpolated with latitudinal and longitudinal dependent regression techniques reflecting the compounds' distributions. With the generated surface concentration climatologies for the ocean and atmosphere, global sea-to-air concentration gradients and sea-to-air fluxes are calculated. Based on these calculations we estimate a total global flux of 1.5/2.5 Gmol Br yr−1 for CHBr3, 0.78/0.98 Gmol Br yr−1 for CH2Br2 and 1.24/1.45 Gmol Br yr−1 for CH3I (robust fit/ordinary least squares regression techniques). Contrary to recent studies, negative fluxes occur in each sea-to-air flux climatology, mainly in the Arctic and Antarctic regions. "Hot spots" for global polybromomethane emissions are located in the equatorial region, whereas methyl iodide emissions are enhanced in the subtropical gyre regions. Inter-annual and seasonal variation is contained within our flux calculations for all three compounds. Compared to earlier studies, our global fluxes are at the lower end of estimates, especially for bromoform. An under-representation of coastal emissions and of extreme events in our estimate might explain the mismatch between our bottom-up emission estimate and top-down approaches.

A latitudinal cross-section and vertical profiles of iodine monoxide (IO) are reported from the marine boundary layer of the Western Pacific. The measurements were taken using Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) during the TransBrom cruise of the German research vessel Sonne, which led from Tomakomai, Japan (42° N, 141° E) through the Western Pacific to Townsville, Australia (19° S, 146° E) in October 2009. In the marine boundary layer within the tropics (between 20° N and 5° S), IO mixing ratios ranged between 1 and 2.2 ppt, whereas in the subtropics and at mid-latitudes typical IO mixing ratios were around 1 ppt in the daytime. The profile retrieval reveals that the bulk of the IO was located in the lower part of the marine boundary layer. Photochemical simulations indicate that the organic iodine precursors observed during the cruise (CH3I, CH2I2, CH2ClI, CH2BrI) are not sufficient to explain the measured IO mixing ratios. Reasonable agreement between measured and modelled IO can only be achieved if an additional sea-air flux of inorganic iodine (e.g., I2) is assumed in the model. Our observations add further evidence to previous studies that reactive iodine is an important oxidant in the marine boundary layer.

This special section of Atmospheric Chemistry and Physics gives an overview of scientific results, collected during a West Pacific ship expedition in October 2009 with the Research Vessel (R/V) Sonne. The cruise focussed on chemical interactions between the ocean surface and the atmosphere above the tropical West Pacific and was planned within the national research project TransBrom (www.geomar.de/~transbrom). TransBrom aimed to particularly investigate very short lived bromine compounds in the ocean and their transport to and relevance for the stratosphere. For this purpose, chemical and biological parameters were analysed in the ocean and in the atmosphere, accompanied by a high frequency of meteorological measurements, to derive new insights into the multidisciplinary research field. This introduction paper presents the scientific goals and the meteorological and oceanographic background. The main research findings of the TransBrom Sonne expedition are highlighted.

Oceanic emissions of halogenated very short-lived substances (VSLS) are expected to contribute significantly to the stratospheric halogen loading and therefore to ozone depletion. The amount of VSLS transported into the stratosphere is estimated based on in-situ observations around the tropical tropopause layer (TTL) and on modeling studies which mostly use prescribed global emission scenarios to reproduce observed atmospheric concentrations. In addition to upper-air VSLS measurements, direct observations of oceanic VSLS emissions are available along ship cruise tracks. Here we use such in-situ observations of VSLS emissions from the West Pacific and tropical Atlantic together with an atmospheric Lagrangian transport model to estimate the direct contribution of bromoform (CHBr3), and dibromomethane (CH2Br2) to the stratospheric bromine loading as well as their ozone depletion potential. Our emission-based estimates of VSLS profiles are compared to upper-air observations and thus link observed oceanic emissions and in situ TTL measurements. This comparison determines how VSLS emissions and transport in the cruise track regions contribute to global upper-air VSLS estimates. The West Pacific emission-based profiles and the global upper-air observations of CHBr3 show a relatively good agreement indicating that emissions from the West Pacific provide an average contribution to the global CHBr3 budget. The tropical Atlantic, although also being a CHBr3 source region, is of less importance for global upper-air CHBr3 estimates as revealed by the small emission-based abundances in the TTL. Western Pacific CH2Br2 emission-based estimates are considerably smaller than upper-air observations as a result of the relatively low sea-to-air flux found in the West Pacific. Together, CHBr3 and CH2Br2 emissions from the West Pacific are projected to contribute to the stratospheric bromine budget with 0.4 pptv Br on average and 2.3 pptv Br for cases of maximum emissions through product and source gas injection. These relatively low estimates reveal that the tropical West Pacific, although characterized by strong convective transport, might overall contribute less VSLS to the stratospheric bromine budget than other regions as a result of only low CH2Br2 and moderate CHBr3 oceanic emissions.

A major stratospheric sudden warming (SSW) in January 2009 was the strongest and most prolonged on record. Aura Microwave Limb Sounder (MLS) observations are used to provide an overview of dynamics and transport during the 2009 SSW, and to compare with the intense, long-lasting SSW in January 2006. The Arctic polar vortex split during the 2009 SSW, whereas the 2006 SSW was a vortex displacement event. Winds reversed to easterly more rapidly and reverted to westerly more slowly in 2009 than in 2006. More mixing of trace gases out of the vortex during the decay of the vortex fragments, and less before the fulfillment of major SSW criteria, was seen in 2009 than in 2006; persistent well-defined fragments of vortex and anticyclone air were more prevalent in 2009. The 2009 SSW had a more profound impact on the lower stratosphere than any previously observed SSW, with no significant recovery of the vortex in that region. The stratopause breakdown and subsequent reformation at very high altitude, accompanied by enhanced descent into a rapidly strengthening upper stratospheric vortex, were similar in 2009 and 2006. Many differences between 2006 and 2009 appear to be related to the different character of the SSWs in the two years.

[1] We present the first complete budget of the interannual variability in Arctic springtime ozone taking into account anthropogenic chemical and natural dynamical processes. For the winters 1991/1992 to 2003/2004 the Arctic chemical ozone loss is available from observations. This work investigates the dynamical supply of ozone to the Arctic polar vortex due to mean transport processes for the same winters. The ozone supply is quantified in a vortex-averaged framework using estimates of diabatic descent over winter. We find that the interannual variability of both dynamical ozone supply and chemical ozone loss contribute, in equal shares, to the variability of the total ozone change. Moreover, together they explain nearly all of the interannual variability of Arctic springtime column ozone. Variability in planetary wave activity, characterized by the Eliassen-Palm flux at 100 hPa, contributes significantly to the variability of ozone supply, chemical ozone loss and total springtime ozone.

Microwave Limb Sounder and Sounding of the Atmosphere with Broadband Emission Radiometry data provide the first opportunity to characterize the four-dimensional stratopause evolution throughout the life-cycle of a major stratospheric sudden warming (SSW). The polar stratopause, usually higher than that at midlatitudes, dropped by ∼30 km and warmed during development of a major “wave 1” SSW in January 2006, with accompanying mesospheric cooling. When the polar vortex broke down, the stratopause cooled and became ill-defined, with a nearly isothermal stratosphere. After the polar vortex started to recover in the upper stratosphere/lower mesosphere (USLM), a cool stratopause reformed above 75 km, then dropped and warmed; both the mesosphere above and the stratosphere below cooled at this time. The polar stratopause remained separated from that at midlatitudes across the core of the polar night jet. In the early stages of the SSW, the strongly tilted (westward with increasing altitude) polar vortex extended into the mesosphere, and enclosed a secondary temperature maximum extending westward and slightly equatorward from the highest altitude part of the polar stratopause over the cool stratopause near the vortex edge. The temperature evolution in the USLM resulted in strongly enhanced radiative cooling in the mesosphere during the recovery from the SSW, but significantly reduced radiative cooling in the upper stratosphere. Assimilated meteorological analyses from the European Centre for Medium-Range weather Forecasts (ECMWF) and Goddard Earth Observing System Version 5.0.1 (GEOS-5), which are not constrained by data at polar stratopause altitudes and have model tops near 80 km, could not capture the secondary temperature maximum or the high stratopause after the SSW; they also misrepresent polar temperature structure during and after the stratopause breakdown, leading to large biases in their radiative heating rates. ECMWF analyses represent the stratospheric temperature structure more accurately, suggesting a better representation of vertical motion; GEOS-5 analyses more faithfully describe stratopause level wind and wave amplitudes. The high-quality satellite temperature data used here provide the first daily, global, multiannual data sets suitable for assessing and, eventually, improving representation of the USLM in models and assimilation systems.

This paper describes the retrievals algorithm used to determine temperature and height from radiance measurements by the Microwave Limb Sounder on EOS Aura. MLS is a "limbscanning" instrument, meaning that it views the atmosphere along paths that do not intersect the surface - it actually looks forwards from the Aura satellite. This means that the temperature retrievals are for a "profile" of the atmosphere somewhat ahead of the satellite. Because of the need to view a finite sample of the atmosphere, the sample spans a box about 1.5km deep and several tens of kilometers in width; the optical characteristics of the atmosphere mean that the sample is representative of a tube about 200-300km long in the direction of view. The retrievals use temperature analyses from NASA's Goddard Earth Observing System, Version 5 (GEOS-5) data assimilation system as a priori states. The temperature retrievals are somewhat deperrde~zt on these a priori states, especially in the lower stratosphere. An important part of the validation of any new dataset involves comparison with other, independent datasets. A large part of this study is concerned with such comparisons, using a number of independent space-based measurements obtained using different techniques, and with meteorological analyses. The MLS temperature data are shown to have biases that vary with height, but also depend on the validation dataset. MLS data are apparently biased slightly cold relative to correlative data in the upper troposphere and slightly warm in the middle stratosphere. A warm MLS bias in the upper stratosphere may be due to a cold bias in GEOS-5 temperatures.

A number of field-campaigns in the tropics have been conducted in the recent years with two different LIDAR systems at Paramaribo in Suriname (5.8° N, 55.2° W). The lidars detect particles in the atmosphere with high vertical and temporal resolution and are capable of detecting extremely thin cloud layers which frequently occur in the tropical tropopause layer (TTL). Radiosonde as well as operational ECMWF analysis show that temperature anomalies caused by equatorial Kelvin waves propagate downward, well below the cold point tropopause (CPT). We find a clear correlation between the temperature anomalies introduced by these waves and the occurrence of thin cirrus in the TTL. In particular we found that extremely thin ice clouds form regularly where cold anomalies shift the tropopause to high altitudes. This finding suggests an influence of Kelvin wave activity on the dehydration in the TTL and thus on the global stratospheric water vapour concentration.

The first three Arctic winters of the ACE mission represented two extremes of winter variability: Stratospheric sudden warmings (SSWs) in 2004 and 2006 were among the strongest, most prolonged on record; 2005 was a record cold winter. Canadian Arctic Atmospheric Chemistry Experiment (ACE) Validation Campaigns were conducted at Eureka (80° N, 86° W) during each of these winters. New satellite measurements from ACE-Fourier Transform Spectrometer (ACE-FTS), Sounding of the Atmosphere using Broadband Emission Radiometry (SABER), and Aura Microwave Limb Sounder (MLS), along with meteorological analyses and Eureka lidar temperatures, are used to detail the meteorology in these winters, to demonstrate its influence on transport, and to provide a context for interpretation of ACE-FTS and validation campaign observations. During the 2004 and 2006 SSWs, the vortex broke down throughout the stratosphere, reformed quickly in the upper stratosphere, and remained weak in the middle and lower stratosphere. The stratopause reformed at very high altitude, near 75 km. ACE measurements covered both vortex and extra-vortex conditions in each winter, except in late-February through mid-March 2004 and 2006, when the strong, pole-centered vortex that reformed after the SSWs resulted in ACE sampling only inside the vortex in the middle through upper stratosphere. The 2004 and 2006 Eureka campaigns were during the recovery from the SSWs, with the redeveloping vortex over Eureka. 2005 was the coldest winter on record in the lower stratosphere, but with an early final warming in mid-March. The vortex was over Eureka at the start of the 2005 campaign, but moved away as it broke up. Disparate temperature profile structure and vortex evolution resulted in much lower (higher) temperatures in the upper (lower) stratosphere in 2004 and 2006 than in 2005. Satellite temperatures agree well with lidar data up to 50–60 km, and ACE-FTS, MLS and SABER show good agreement in high-latitude temperatures throughout the winters. Consistent with a strong, cold upper stratospheric vortex and enhanced radiative cooling after the SSWs, MLS and ACE-FTS trace gas measurements show strongly enhanced descent in the upper stratospheric vortex in late January through March 2006 compared to that in 2005.

A multiyear time series of the vortex-averaged diabatic descent for 47 Arctic winters from 1957/1958 until 2003/2004 is presented. The climatology of diabatic descent is based on trajectory calculations coupled with diabatic heating rate calculations carried out in the polar lower stratosphere of the Northern Hemisphere winters. We demonstrate the improved performance of the approach based on diabatic heating rates compared to the approach based on vertical winds from meteorological analysis. The time series of the vortex-averaged diabatic descent gives a detailed picture of intensity and altitude dependence of the stratospheric vertical transport processes during the Arctic winter. In addition to the overall vortex-averaged diabatic descent, the spatial structure of the descent is analyzed for two different Arctic winters. We demonstrate for this case study that not only the intensity but also the zonal structure of the diabatic descent depends on the meteorological conditions in the polar vortex. The climatology is characterized by very pronounced interannual variability which is linked to the variability of temperature anomalies and to the variability of Eliassen-Palm (EP)-flux anomalies, wherein strong planetary wave activity leads to strong diabatic descent and vice versa. The correlation between EP-flux and descent shows that tropospheric dynamics have a strong influence on the strength of the polar branch of the residual circulation by means of the atmospheric wave activity.

In the framework of the European Project STAR the mobile Aerosol Raman Lidar (MARL) of the Alfred Wegener Institute (AWI) was operated in Paramaribo/Suriname (5.8°, 55.2°W) and carried out extensive observations of tropical cirrus clouds during the fall dry season from September 28 to November 15 2004. With 81% the coverage with ice clouds was very high in the upper troposphere (above 12 km). The frequency of occurrence of subvisible clouds was found to be clearly enhanced compared to the midlatitudes. The extinction-to-backscatter ratio of thin tropical cirrus is with 26 ± 7 sr significantly higher than that of mid-latitude cirrus (16 ± 9 sr).Subvisible cirrus occur mostly in the tropical tropopause layer (TTL) abovean upper tropospheric inversion. The formation conditions of these cloudsand the dehydration of the TTL was investigated by means of a newly developed trajectory model. Clouds were assumed to occur when the relativehumidity above ice (RHI) exceeds 100% due to adiabatic cooling and subsequently the humidity is reduced to the saturation vapor pressure. This simple dehydration scheme yields an amazingly good agreement with cloud andhumidity observations at Paramaribo. However, it challenges the way cirrusformation is generally conceived, since high supersaturation were thoughtto be required. Our results suggest, that the inclination of the TTL to contain solid particles is rather high, even under modestly moist conditions. Wealso detected extremely thin layers of solid particles slightly above the temperature minimum in sub-saturated air which are possibly stabilized by HNO3.

During the European heat wave summer 2003 with predominant high pressure conditions we performed a detailed study of upper tropospheric humidity and ice particles which yielded striking results concerning the occurrence of ice supersaturated regions (ISSR), cirrus, and contrails. Our study is based on lidar observations and meteorological data obtained at Lindenberg/Germany (52.2° N, 14.1° E) as well as the analysis of the European centre for medium range weather forecast (ECMWF). Cirrus clouds were detected in 55% of the lidar profiles and a large fraction of them were subvisible (optical depth

A long-term climatology of air mass transport through the tropical tropopause layer (TTL) is presented, covering the period from 1962–2005. The transport through the TTL is calculated with a Lagrangian approach using radiative heating rates as vertical velocities in an isentropic trajectory model. We demonstrate the improved performance of such an approach compared to previous studies using vertical winds from meteorological analyses. Within the upper part of the TTL, the averaged diabatic ascent is 0.5 K/day during Northern Hemisphere (NH) winters 1992–2001. Climatological maps show a cooling and strengthening of this part of the residual circulation during the 1990s and early 2000s compared to the long-term mean. Lagrangian cold point (LCP) fields show systematic differences for varying time periods and natural forcing components. The interannual variability of LCP temperature and density fields is found to be influenced by volcanic eruptions, El Niño Southern Oscillation (ENSO), Quasi-Biennial Oscillation (QBO) and the solar cycle. The coldest and driest TTL is reached during QBO easterly phase and La Niña over the western Pacific, whereas during volcanic eruptions, El Niño and QBO westerly phase it is warmer and less dry.

The 2003-2004 Arctic winter was remarkable in the 40-year record of meteorological analyses. A major warming beginning in early January 2004 led to nearly two months of vortex disruption with high-latitude easterlies in the middle to lower stratosphere. The upper stratospheric vortex broke up in late December, but began to recover by early January, and in February and March was the strongest since regular observations began in 1979. The lower stratospheric vortex broke up in late January. Comparison with two previous years, 1984-1985 and 1986-1987, with prolonged mid-winter warming periods shows unique characteristics of the 2003-2004 warming period: The length of the vortex disruption, the strong and rapid recovery in the upper stratosphere, and the slow progression of the warming from upper to lower stratosphere. January 2004 zonal mean winds in the middle and lower stratosphere were over two standard deviations below average. Examination of past variability shows that the recent frequency of major stratospheric warmings (seven in the past six years) is unprecedented. Lower stratospheric temperatures were unusually high during six of the past seven years, with five having much lower than usual potential for PSC formation and ozone loss (nearly none in 1998-1999, 2001-2002 and 2003-2004, and very little in 1997-1998 and 2000-2001). Middle and upper stratospheric temperatures, however, were unusually low during and after February. The pattern of five of the last seven years with very low PSC potential would be expected to occur randomly once every approximately 850 years. This cluster of warm winters, immediately following a period of unusually cold winters, may have important implications for possible changes in interannual variability and for determination and attribution of trends in stratospheric temperatures and ozone.