A refined method for calculating equivalent effective stratospheric chlorine

Chlorine and bromine atoms lead to catalytic depletion of ozone in the stratosphere. Therefore the use and production of ozone-depleting substances (ODSs) containing chlorine and bromine is regulated by the Montreal Protocol to protect the ozone layer. Equivalent effective stratospheric chlorine (EESC) has been adopted as an appropriate metric to describe the combined effects of chlorine and bromine released from halocarbons on stratospheric ozone. Here we revisit the concept of calculating EESC. We derive a refined formulation of EESC based on an advanced concept of ODS propagation into the stratosphere and reactive halogen release. A new transit time distribution is introduced in which the age spectrum for an inert tracer is weighted with the release function for inorganic halogen from the source gases. This distribution is termed the "release time distribution". We show that a much better agreement with inorganic halogen loading from the chemistry transport model TOMCAT is achieved compared with using the current formulation. The refined formulation shows EESC levels in the year 1980 for the mid-latitude lower stratosphere, which are significantly lower than previously calculated. The year 1980 is commonly used as a benchmark to which EESC must return in order to reach significant progress towards halogen and ozone recovery. Assuming that – under otherwise unchanged conditions – the EESC value must return to the same level in order for ozone to fully recover, we show that it will take more than 10 years longer than estimated in this region of the stratosphere with the current method for calculation of EESC. We also present a range of sensitivity studies to investigate the effect of changes and uncertainties in the fractional release factors and in the assumptions on the shape of the release time distributions. We further discuss the value of EESC as a proxy for future evolution of inorganic halogen loading under changing atmospheric dynamics using simulations from the EMAC model. We show that while the expected changes in stratospheric transport lead to significant differences between EESC and modelled inorganic halogen loading at constant mean age, EESC is a reasonable proxy for modelled inorganic halogen on a constant pressure level

Transport timescales and tracer properties in the extratropical UTLS

A comprehensive evaluation of seasonal backward trajectories initialized in the northern hemisphere lowermost stratosphere (LMS) has been performed to investigate the factors that determine the temporal and spatial structure of troposphere-to-stratosphere-transport (TST) and it's impact on the LMS. In particular we explain the fundamental role of the transit time since last TST (tTST) for the chemical composition of the LMS. According to our results the structure of the LMS can be characterized by a layer with tTST<40 days forming a narrow band around the local tropopause. This layer extends about 30 K above the local dynamical tropopause, corresponding to the extratropical tropopause transition layer (ExTL) as identified by CO. The LMS beyond this layer shows a relatively well defined separation as marked by an aprupt transition to longer tTST indicating less frequent mixing and a smaller fraction of tropospheric air. Thus the LMS constitutes a region of two well defined regimes of tropospheric influence. These can be characterized mainly by different transport times from the troposphere and different fractions of tropospheric air. Carbon monoxide (CO) mirrors this structure of tTST due to it's finite lifetime on the order of three months. Water vapour isopleths, on the other hand, do not uniquely indicate TST and are independent of tTST, but are determined by the Lagrangian Cold Point (LCP) of air parcels. Most of the backward trajectories from the LMS experienced their LCP in the tropics and sub-tropics, and TST often occurs 20 days after trajectories have encountered their LCP. Therefore, ExTL properties deduced from CO and H2O provide totally different informations on transport and particular TST for the LMS

Corrigendum to "Fractional release factors of long-lived halogenated organic compounds in the tropical stratosphere" published in Atmos. Chem. Phys., 10, 1093–1103, 2010

We have noted a computational error in the calculation of the averaged Fractional Release Factors (FRFs) relative to the averaged FRF of CFCl3 (CFC-11). ..

Observation-based assessment of stratospheric fractional release, lifetimes, and ozone depletion potentials of ten important source gases

Estimates of the recovery time of stratospheric ozone heavily rely on the exact knowledge of the processes that lead to the decomposition of the relevant halogenated source gases. Crucial parameters in this context are fractional release factors (FRFs) as well as stratospheric lifetimes and ozone depletion potentials (ODPs). We here present data from the analysis of air samples collected between 2009 and 2011 on board research aircraft flying in the mid- and high-latitude stratosphere and infer the above-mentioned parameters for ten major source gases: CFCl3 (CFC-11), CF2Cl2 (CFC-12), CF2ClCFCl2 (CFC-113), CCl4 (carbon tetrachloride), CH3CCl3 (methyl chloroform), CHF2Cl (HCFC-22), CH3CFCl2 (HCFC-141b), CH3CF2Cl (HCFC-142b), CF2ClBr (H-1211), and CF3Br (H-1301). The inferred correlations of their FRFs with mean ages of air reveal less decomposition as compared to previous studies for most compounds. When using the calculated set of FRFs to infer equivalent stratospheric chlorine, we find a reduction of more than 20% as compared to the values inferred in the most recent Scientific Assessment of Ozone Depletion by the World Meteorological Organisation (WMO, 2011). We also note that FRFs and their correlations with mean age are not generally time-independent as often assumed. The stratospheric lifetimes were calculated relative to that of CFC-11. Within our uncertainties the ratios between stratospheric lifetimes inferred here agree with the values in recent WMO reports except for CFC-11, CFC-12 and CH3CCl3. Finally, we calculate lower ODPs than recommended by WMO for six out of ten compounds, with changes most pronounced for the three HCFCs. Collectively these newly calculated values may have important implications for the severity and recovery time of stratospheric ozone loss

Short-lived brominated hydrocarbons – observations in the source regions and the tropical tropopause layer

We conducted measurements of the five important short-lived organic bromine species in the marine boundary layer (MBL). Measurements were made in the Northern Hemisphere mid-latitudes (Sylt Island, North Sea) in June 2009 and in the tropical Western Pacific during the TransBrom ship campaign in October 2009. For the one-week time series on Sylt Island, mean mixing ratios of CHBr3, CH2Br2, CHBr2Cl and CH2BrCl were 2.0, 1.1, 0.2, 0.1 ppt, respectively. We found maxima of 5.8 and 1.6 ppt for the two main components CHBr3 and CH2Br2. Along the cruise track in the Western Pacific (between 41° N and 13° S) we measured mean mixing ratios of 0.9, 0.9, 0.2, 0.1 and 0.1 ppt for CHBr3, CH2Br2, CHBrCl2, CHBr2Cl and CH2BrCl. Air samples with coastal influence showed considerably higher mixing ratios than the samples with open ocean origin. Correlation analyses of the two data sets yielded strong linear relationships between the mixing ratios of four of the five species (except for CH2BrCl). Using a combined data set from the two campaigns and a comparison with the results from two former studies, rough estimates of the molar emission ratios between the correlated substances were: 9/1/0.35/0.35 for CHBr3/CH2Br2/CHBrCl2/CHBr2Cl. Additional measurements were made in the tropical tropopause layer (TTL) above Teresina (Brazil, 5° S) in June 2008, using balloon-borne cryogenic whole air sampling technique. Near the level of zero clear-sky net radiative heating (LZRH) at 14.8 km about 2.25 ppt organic bromine was bound to the five short-lived species, making up 13% of total organic bromine (17.82 ppt). CH2Br2 (1.45 ppt) and CHBr3 (0.56 ppt) accounted for 90% of the budget of short-lived compounds in that region. Near the tropopause (at 17.5 km) organic bromine from these substances was reduced to 1.35 ppt, with 1.07 and 0.12 ppt attributed to CH2Br2 and CHBr3, respectively

Observation of mesospheric air inside the arctic stratospheric polar vortex in early 2003

During several balloon flights inside the Arctic polar vortex in early 2003, unusual trace gas distributions were observed, which indicate a strong influence of mesospheric air in the stratosphere. The tuneable diode laser (TDL) instrument SPIRALE (Spectroscopie InFrarouge par Absorption de Lasers Embarqués) measured unusually high CO values (up to 600 ppb) on 27 January at about 30 km altitude. The cryosampler BONBON sampled air masses with very high molecular Hydrogen, extremely low SF6 and enhanced CO values on 6 March at about 25 km altitude. Finally, the MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) Fourier Transform Infra-Red (FTIR) spectrometer showed NOy values which are significantly higher than NOy* (the NOy derived from a correlation between N2O and NOy under undisturbed conditions), on 21 and 22 March in a layer centred at 22 km altitude. Thus, the mesospheric air seems to have been present in a layer descending from about 30 km in late January to 25 km altitude in early March and about 22 km altitude on 20 March. We present corroborating evidence from a model study using the KASIMA (KArlsruhe Simulation model of the Middle Atmosphere) model that also shows a layer of mesospheric air, which descended into the stratosphere in November and early December 2002, before the minor warming which occurred in late December 2002 lead to a descent of upper stratospheric air, cutting of a layer in which mesospheric air is present. This layer then descended inside the vortex over the course of the winter. The same feature is found in trajectory calculations, based on a large number of trajectories started in the vicinity of the observations on 6 March. Based on the difference between the mean age derived from SF6 (which has an irreversible mesospheric loss) and from CO2 (whose mesospheric loss is much smaller and reversible) we estimate that the fraction of mesospheric air in the layer observed on 6 March, must have been somewhere between 35% and 100%

Sensitivity of age of air trends to the derivation method for non-linear increasing inert SF6

Mean age of air (AoA) is a diagnostic of transport along the stratospheric Brewer–Dobson circulation. While models consistently show negative trends, long-term time series (1975–2016) of AoA derived from observations show non-significant positive trends in mean AoA in the Northern Hemisphere. This discrepancy between observed and modelled mean AoA trends is still not resolved. There are uncertainties and assumptions required when deriving AoA from trace gas observations. At the same time, AoA from climate models is subject to uncertainties, too. In this paper, we focus on the uncertainties due to the parameter selection in the method that is used to derive mean AoA from SF$_{6}$ measurements in Engel et al. (2009, 2017). To correct for the non-linear increase in SF$_{6}$ concentrations, a quadratic fit to the time series at the reference location, i.e. the tropical surface, is used. For this derivation, the width of the AoA distribution (age spectrum) has to be assumed. In addition, to choose the number of years the quadratic fit is performed for, the fraction of the age spectrum to be considered has to be assumed. Even though the uncertainty range due to all different aspects has already been taken into account for the total errors in the AoA values, the systematic influence of the parameter selection on AoA trends is described for the first time in the present study. For this, we use the EMAC (ECHAM MESSy Atmospheric Chemistry) climate model as a test bed, where AoA derived from a linear tracer is available as a reference and modelled age spectra exist to diagnose the actual spatial age spectra widths. The comparison of mean AoA from the linear tracer with mean AoA from a SF$_{6}$ tracer shows systematic deviations specifically in the trends due to the selection of the parameters. However, for an appropriate parameter selection, good agreement for both mean AoA and its trend can be found, with deviations of about 1 % in mean AoA and 12 % in AoA trend. In addition, a method to derive mean AoA is evaluated that applies a convolution to the reference time series. The resulting mean AoA and its trend only depend on an assumption about the ratio of moments. Also in that case, it is found that the larger the ratio of moments, the more the AoA trend gravitates towards the negative. The linear tracer and SF$_{6}$ AoA are found to agree within 0.3 % in the mean and 6 % in the trend. The different methods and parameter selections were then applied to the balloon-borne SF$_{6}$ and CO$_{2}$ observations. We found the same systematic changes in mean AoA trend dependent on the specific selection. When applying a parameter choice that is suggested by the model results, the AoA trend is reduced from 0.15 to 0.07 years per decade. It illustrates that correctly constraining those parameters is crucial for correct mean AoA and trend estimates and still remains a challenge in the real atmosphere

A convolution of observational and model data to estimate age of air spectra in the northern hemispheric lower stratosphere

Derivation of mean age of air (AoA) and age spectra from atmospheric measurements remains a challenge and often requires output from atmospheric models. This study tries to minimize the direct influence of model output and presents an extension and application of a previously established inversion method to derive age spectra from mixing ratios of long- and short-lived trace gases. For a precise description of cross-tropopause transport processes, the inverse method is extended to incorporate air entrainment into the stratosphere across the tropical and extratropical tropopause. We first use simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) to provide a general proof of concept of the extended principle in a controllable and consistent environment, where the method is applied to an idealized set of 10 trace gases with predefined constant lifetimes and compared to reference model age spectra. In the second part of the study we apply the extended inverse method to atmospheric measurements of multiple long- and short-lived trace gases measured aboard the High Altitude and Long Range (HALO) research aircraft during the two research campaigns POLSTRACC–GW-LCYCLE–SALSA (PGS) and Wave-driven Isentropic Exchange (WISE). As some of the observed species undergo significant loss processes in the stratosphere, a Monte Carlo simulation is introduced to retrieve age spectra and chemical lifetimes in stepwise fashion and to account for the large uncertainties. Results show that in the idealized model scenario the inverse method retrieves age spectra robustly on annual and seasonal scales. The extension to multiple entry regions proves reasonable as our CLaMS simulations reveal that in the model between 50 % and 70 % of air in the lowermost stratosphere has entered through the extratropical tropopause (30–90∘ N and S) on annual average. When applied to observational data of PGS and WISE, the method derives age spectra and mean AoA with meaningful spatial distributions and quantitative range, yet large uncertainties. Results indicate that entrainment of fresh tropospheric air across both the extratropical and tropical tropopause peaked prior to both campaigns, but with lower mean AoA for WISE than PGS data. The ratio of moments for all retrieved age spectra for PGS and WISE is found to range between 0.52 and 2.81 years. We conclude that the method derives reasonable and consistent age spectra using observations of chemically active trace gases. Our findings might contribute to an improved assessment of transport with age spectra in future studies

Modelling marine emissions and atmospheric distributions of halocarbons and dimethyl sulfide: the influence of prescribed water concentration vs. prescribed emissions

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 spee

Age spectra and other transport diagnostics in the North American monsoon UTLS from SEAC4^{4}RS in situ trace gas measurements

The upper troposphere and lower stratosphere (UTLS) region during the summer monsoon season over North America (NAM) is influenced by the transport of air from a variety of source regions over a wide range of timescales (hours to years). Age spectra are useful for characterizing the transport into such a region, and in this study we use and build on recently developed techniques to infer age spectra from trace gas measurements with photochemical lifetimes from days to centuries. We show that the measurements taken by the whole-air sampler instrument during the SEAC$^{4}$RS campaign can be used to derive not only age spectra, but also path-integrated lifetimes of each of the trace gases and partitioning between North American and tropical surface source origins. The method used here can also clearly identify and adjust for measurement outliers that were influenced by polluted surface source regions. The results are generally consistent with expected transport features of the NAM but also provide a range of transport diagnostics (age spectra, trace gas lifetimes and surface source regions) that have not previously been computed solely from in situ measurements. These methods may be applied to many other existing in situ datasets, and the transport diagnostics can be compared with chemistry–climate model transport in the UTLS