the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Correction of stratospheric age of air (AoA) derived from sulfur hexafluoride (SF_{6}) for the effect of chemical sinks
Hella Garny
Roland Eichinger
Johannes C. Laube
Eric A. Ray
Gabriele P. Stiller
Harald Bönisch
Laura Saunders
Marianna Linz
Observational monitoring of the stratospheric transport circulation, the Brewer–Dobson circulation (BDC), is crucial to estimate any decadal to longterm changes therein, a prerequisite to interpret trends in stratospheric composition and to constrain the consequential impacts on climate. The transport time along the BDC (i.e. the mean stratospheric age of air, AoA) can best be deduced from trace gas measurements of tracers which increase linearly with time and are chemically passive. The gas sulfur hexafluoride (SF_{6}) is often used to deduce AoA because it has been increasing monotonically since the ∼1950s, and previously its chemical sinks in the mesosphere have been assumed to be negligible for AoA estimates. However, recent studies have shown that the chemical sinks of SF_{6} are stronger than assumed and become increasingly relevant with rising SF_{6} concentrations.
To adjust biases in AoA that result from the chemical SF_{6} sinks, we here propose a simple correction scheme for SF_{6}based AoA estimates accounting for the timedependent effects of chemical sinks. The correction scheme is based on theoretical considerations with idealized assumptions, resulting in a relation between ideal AoA and apparent AoA which is a function of the tropospheric reference time series of SF_{6} and of the AoAdependent effective lifetime of SF_{6}. The correction method is thoroughly tested within a selfconsistent data set from a climate model that includes explicit calculation of chemical SF_{6} sinks. It is shown within the model that the correction successfully reduces biases in SF_{6}based AoA to less than 5 % for mean ages below 5 years. Tests using only subsampled data for deriving the fit coefficients show that applying the correction scheme even with imperfect knowledge of the sink is far superior to not applying a sink correction.
Furthermore, we show that based on currently available measurements, we are not able to constrain the fit parameters of the correction scheme based on observational data alone. However, the modelbased correction curve lies within the observational uncertainty, and we thus recommend using the modelderived fit coefficients until more highquality measurements are able to further constrain the correction scheme. The application of the correction scheme to AoA from satellites and in situ data suggests that it is highly beneficial to reconcile different observational estimates of mean AoA.
 Article
(6760 KB)  Fulltext XML

Supplement
(352 KB)  BibTeX
 EndNote
The distributions of trace gases in the stratosphere are strongly influenced by the slow Equatortopole overturning circulation in the stratosphere, the Brewer–Dobson circulation (BDC). Most prominently, total column abundances of stratospheric ozone are maximized in mid to high latitudes despite the strongest chemical production being in the tropics. The drive of this overturning circulation by wave dissipation, acting to drive the slow adiabatic rising and sinking as well as inducing wave stirring that subsequently leads to mixing, has been understood for many decades (Haynes et al., 1991; Butchart, 2014). However, dynamical measures of the overturning circulation like the residual circulation velocities are not observable. The distributions of trace gases rather serve to deduce measures of the circulation strength – fittingly given the initial definition of the BDC as tracer transport circulation (Brewer, 1949; Dobson, 1956).
A wellestablished measure to quantify BDC strength is the mean transport time from a reference surface (e.g. the tropical tropopause or the Earth's surface) to a point in the stratosphere, referred to as the mean stratospheric age of air (AoA) (Waugh and Hall, 2002). One advantage of this AoA measure is that it can be directly deduced from trace gas abundances; given an inert tracer whose surface abundance increases linearly with time, AoA is simply the lag time between the mixing ratio at the surface and the mixing ratio at a given point in the stratosphere. Strictly, this measures the transport time both within the troposphere and in the stratosphere, but as transport times from the surface to the tropical tropopause can be neglected compared to stratospheric transport times, AoA relative to the surface is still a good measure of the stratospheric transport circulation. However, such an ideal AoA tracer does not exist in reality. One gas that has often been used to derive AoA from observations is sulfur hexafluoride (SF_{6}). SF_{6} is a wellsuited AoA tracer since its surface abundances have increased monotonically since the ∼1950s, and its chemical sinks are located in the mesosphere and above, which have previously been assumed to have little influence on SF_{6} abundances outside the winter polar vortices. A number of highquality measurements of SF_{6} from different platforms (in situ and satellite) have existed for many decades up until today (e.g. Andrews et al., 2001; Engel et al., 2009; Stiller et al., 2012; Haenel et al., 2015). While SF_{6} surface abundances do not increase linearly, it has been shown that the nonlinear increase in SF_{6} can be corrected for in the derivation of AoA (Stiller et al., 2012; Fritsch et al., 2019). SF_{6} has the advantage compared to the other common age tracer CO_{2} that it has almost no seasonal cycle in emissions. The strong seasonal cycle in CO_{2} complicates the AoA calculation in particular in the lower stratosphere (e.g. Andrews et al., 2001). Thus, SF_{6} is, in principle, a useful and oftenmeasured AoA tracer. However, it has become increasingly clear in recent years that the chemical sinks of SF_{6} introduce substantial biases in AoA (e.g. Stiller et al., 2012; Ray et al., 2017; Kovács et al., 2017; Leedham Elvidge et al., 2018; Loeffel et al., 2022). It is strongly desirable to be able to correct for the effects of chemical sinks in SF_{6} in order to use the existing and possible future observations of SF_{6} to monitor the stratospheric transport circulation.
Loeffel et al. (2022) investigated the effects of chemical sinks on AoA derived from SF_{6} in a global model and showed that the bias of AoA induced by the chemical sink grows with increasing mixing ratios. Thus, biases have been particularly strong since the 2000s. Furthermore, this effect can induce apparent positive trends in AoA, hampering efforts to deduce longterm circulation trends from SF_{6}based AoA. In Loeffel et al. (2022) it was discussed that there is a compact relationship between ideal AoA and SF_{6}based apparent AoA for AoA below a certain threshold. On this basis, it was suggested that it might be possible to correct observed SF_{6}based AoA for the effects of chemical sinks. Here we follow up on this idea and develop a correction scheme for SF_{6}based AoA. Firstly, theoretical considerations with simplifications are used to deduce a possible formulation of the correction scheme (Sect. 2). The strategy of the paper is to develop and test the correction scheme thoroughly within the selfconsistent model world (Sect. 4, with the model being described in Sect. 3.1). This includes subsampling the model data to test the robustness of the correction against incomplete knowledge of the ideal AoA to apparent AoA relationship (Sect. 4.3). The next step is a comparison of the relation of SF_{6}based AoA to AoA deduced from other tracers from observational data to the model relation. Ideally, the fit parameters for the correction scheme would be based on observational data alone; however, we show that current observations are not able to constrain the fit parameters (Sect. 5). Instead, we use fit parameters from the model and apply them to independent observations to test the performance of the correction scheme (Sect. 5). In Sect. 6 we summarize our findings, discuss the shortcomings and give recommendations on the application of the correction scheme.
In this section, we will lay out the theoretical basis for how chemical depletion of a tracer affects the derivation of mean AoA from this tracer. We will make various simplifying assumptions to derive analytical expressions for the effects of chemical sinks on tracer abundances and on the AoA derived from them. The climate model data introduced in Sect. 3.1 will serve to test how welljustified those assumptions are.
Generally, for any location, we can express the stratospheric tracer mixing ratio χ_{s} at time t for a tracer that is transported from a reference surface as
where t^{′} denotes the transit time, χ_{0}(t) is the time series of the tracer mixing ratio at a reference surface, and G(t^{′}) represents Green's function and is equivalent to the distribution of transit times (often called the age spectrum). Note that Green's function is generally dependent on location x and on time t (hereafter those indices will be dropped for brevity without loss of generality). The exponential function in the integral represents the effects of chemical sinks on the tracer with transittimedependent and pathintegrated lifetime τ. For simplicity, we will first consider a tracer with a linear increase at the reference surface before moving on to a tracer with a monotonic but nonlinear increase.
2.1 Linearly increasing tracer
We assume a linearly increasing tracer with ${\mathit{\chi}}_{\mathrm{0}}={\mathit{\delta}}_{{\mathit{\chi}}_{\mathrm{0}}}t$. Given the definition of mean AoA as the first moment of the age spectrum ($\mathrm{\Gamma}={\int}_{{t}^{\prime}=\mathrm{0}}^{\mathrm{\infty}}{t}^{\prime}G\left({t}^{\prime}\right)\mathrm{d}{t}^{\prime}$), AoA from a linearly increasing tracer can be calculated as the lag time between the tracer mixing ratio at the location of interest and the time series at the reference surface ($\mathrm{\Gamma}=t{\mathit{\chi}}_{\mathrm{s}}\left(t\right)/{\mathit{\delta}}_{{\mathit{\chi}}_{\mathrm{0}}}$). We will refer to such estimated AoA without taking the effects of chemical sinks into account as the apparent AoA $\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}$. In this case, the apparent AoA is given by
We will now make two assumptions. Firstly, we assume that the mean AoA Γ is small compared to lifetime τ, which is generally the case for SF_{6} within the stratosphere (in which case lifetime is on the order of 100 to 1000 years, and the mean AoA is less than 10 years). Thus, we can approximate $\mathrm{exp}\left(\frac{{t}^{\prime}}{\mathit{\tau}}\right)\approx \mathrm{1}\frac{{t}^{\prime}}{\mathit{\tau}}$. Secondly, we will make the assumption that the transittimedependent lifetime τ(t^{′}) can be replaced by the pathaverage lifetime τ_{eff}. This lifetime τ_{eff} is the pathaverage lifetime across all possible paths to the point of interest of the pathintegrated lifetime τ(t^{′}), and we refer to this as the “effective” lifetime. Note that similar formulations have been used previously by Holzer and Waugh (2015), where they refer to this quantity as the “pathdependent inverse loss frequency”, and by Schoeberl et al. (2000) in their “average path assumption”.
With the two assumptions stated above, we can simplify Eq. (2) to
Using the expression for the width of the age spectrum ${\mathrm{\Delta}}^{\mathrm{2}}=\frac{\mathrm{1}}{\mathrm{2}}{\int}_{{t}^{\prime}=\mathrm{0}}^{\mathrm{\infty}}({t}^{\prime}\mathrm{\Gamma}{)}^{\mathrm{2}}G({t}^{\prime})\mathrm{d}{t}^{\prime}$, one can rewrite the latter term as $\frac{\mathrm{1}}{{\mathit{\tau}}_{\mathrm{eff}}}\left(\mathrm{2}{\mathrm{\Delta}}^{\mathrm{2}}+{\mathrm{\Gamma}}^{\mathrm{2}}\right)$ and after some rearrangement obtain
The latter approximation is justified given that the time t since start of emission of the tracer (for SF_{6} about 1950s) is around 1 order of magnitude larger (on the order of 10 years) than the mean AoA Γ and the ratio of moments $\frac{{\mathrm{\Delta}}^{\mathrm{2}}}{\mathrm{\Gamma}}$ (which are both on the order of 1; see e.g. Fritsch et al., 2019). This approximation holds less well in earlier years (1960s to 1970s), but given that the correction at those times is small, this poses no major issue. Thus, this approximated expression states that the apparent AoA differs from the true AoA Γ by a factor of $(\mathrm{1}+\frac{t}{{\mathit{\tau}}_{\mathrm{eff}}})$, which increases as time progresses (due to the increasing tracer mixing ratios; see discussion in Loeffel et al., 2022). The other necessary parameter for estimating the effects of chemical sinks on AoA is the effective lifetime. Since chemical sinks of SF_{6} are located in the upper stratosphere and the mesosphere, the effective lifetime is locationdependent. For the correction scheme proposed here, we will parameterize the effective lifetime as a function of AoA. The relation of the effective lifetime to AoA is expected to bear latitudinal dependence, as is discussed in Sect. 4.1. The effective lifetime can be calculated from simultaneous data of SF_{6}derived apparent AoA and AoA based on another tracer that is unaffected by the chemical sinks via the relation derived above:
We deduce this relationship based on model data and then parameterize the effective lifetime as a function of ideal AoA (Γ) via different fit functions (exponential and polynomial), and the estimated parameters then serve as input for the correction scheme.
2.2 Nonlinearly increasing tracer
In the next step, we extend the formulations given above to correct for a nonlinearly increasing tracer, as is the case for SF_{6}. To account for the nonlinear increase in the reference time series, we linearize the tropospheric time series ${\mathit{\chi}}_{\mathrm{0}}(t{t}^{\prime})$ around the average transit time ${t}^{\prime}={\mathrm{\Gamma}}_{\mathrm{0}}$ so that
Inserting this expression into Eq. (1) for a passive tracer named ${\mathit{\chi}}_{\mathrm{s}}^{\mathrm{p}}$ (i.e. $\mathrm{exp}({t}^{\prime}/\mathit{\tau}({t}^{\prime}\left)\right)=\mathrm{1}$) and solving for Γ yields
Note that this is not the common way to account for nonlinearity when deriving AoA from a nonlinearly increasing tracer. Usually, either a convolution with the age spectrum or a secondorder polynomial fit to the reference time series is used (Engel et al., 2009; Stiller et al., 2012; Fritsch et al., 2019). The simple linearization of the tropospheric reference time series only approximates the local slope. For the mean age calculation, the timeevolving slope over all transit times that contribute to the age spectrum is important. However, for the sink correction we find that the approximation of linearization around a certain mean age proves to work well (as shown in the following in Sect. 4.2).
To obtain the relationship between apparent age and the true ideal age, we insert the chemically depleted tracer mixing ratio χ_{s}(t) instead of the passive tracer in Eq. (7), which after some rearrangement gives
We define a timedependent function F_{t}(t) as the ratio of the reference mixing ratio to its temporal gradient:
We make the same assumptions as for Eq. (4), yielding the following approximate relation of apparent AoA to ideal AoA:
The latter approximation is even slightly betterjustified than for the linear tracer (see Eq. 4) because of the near cancellation of Γ and Γ_{0}. Thus, we obtain a similar expression as for the linearly increasing tracer (Eq. 4) but replace the linear time dependence t with the timedependent function F_{t}(t). If the reference tracer mixing ratio increases linearly, Eq. (10) collapses to Eq. (4).
2.3 Correction scheme
The correction scheme we propose here is based on Eq. (10) and contains the following steps:

For any pair of given data of apparent and ideal AoA (${\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}}_{i}$, Γ_{i}) observed at a certain time t_{i}, the effective lifetime is calculated based on ${\mathit{\tau}}_{\mathrm{eff}}={F}_{t}\left({t}_{i}\right)/({\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}}_{i}/{\mathrm{\Gamma}}_{i}\mathrm{1})$.

A function is fitted to parameterize τ_{eff} as a function of ideal AoA (Γ). In Sect. 4.1, we test an exponential fit function and polynomial fit functions of different orders. The fits are performed over a limited range of ideal AoA values (in this case 1 to 5 years; see Sect. 4.1).

Given the parameters from the fit and the timedependent function based on the reference time series F_{t}(t), the relation between $\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}$ and Γ is given for any time and can be used to correct SF_{6}based apparent AoA. Since the relation between $\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}$ and Γ cannot be analytically inverted (due to the nonlinear dependence of τ_{eff} on Γ), we calculate $\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}$ for a range of ideal mean ages in practice, and Γ is estimated from $\stackrel{\mathrm{\u0303}}{\mathrm{\Gamma}}$ via linear interpolation. AoA values beyond the range for which the fit was performed in the previous step (in this case for ideal AoA > 5 years) are treated in the following way – in the case of the exponential fit, the exponential fit function is evaluated also beyond AoA > 5 years. In the case of the polynomial fit, this would lead to problems as the polynomial is not well constrained outside the fit range and does not increase monotonically. Therefore, the polynomial is not evaluated for AoA > 5 years, but the values are extrapolated linearly.
2.4 Timedependent function for nonlinear increase
The timedependent function F_{t}(t) takes the nonlinear increase in the reference mixing ratio time series into account and is defined as the ratio of the mixing ratio and its temporal slope at a given time as given by Eq. (9).
The reference SF_{6} time series at the surface is shown in Fig. 1a for both the best estimate based on surface observations (orange) and the time series taken from the global model data (blue). As can be seen, the two time series are close to identical. Nevertheless, in the remainder of the paper the model time series is used for model SF_{6} data, and the observational time series is used for correcting observed SF_{6}based AoA (in Sect. 5). Figure 1b shows the temporal gradient in the reference time series. Function F_{t}(t) (Fig. 1c) is calculated based on the annual mean time series of SF_{6}, where the gradient in the denominator of F_{t}(t) is calculated as the centred difference over 5 years. We find that in practice the offset mean AoA (Γ_{0}) of 4 years gives the best results, but generally results were found to be robust with respect to small variations in the choice of the offset mean AoA and also of the period for the gradient calculation. The values of F_{t} from the model versus observational time series only show differences in the earlier years (the 1950s and less so in the 1970s), which have negligible effects on the results presented in the remainder of the paper.
3.1 The ECHAM/MESSy Atmospheric Chemistry (EMAC) time slice simulation
We use the ECHAM/MESSy Atmospheric Chemistry (EMAC) model, which is a numerical chemistry and climate simulation system that includes submodels describing tropospheric and middleatmospheric processes and their interaction with oceans, land and human influences (Jöckel et al., 2010). It contains the fifthgeneration European Centre Hamburg general circulation model (ECHAM5, Roeckner et al., 2006) as the dynamical core atmospheric model and uses the second version of the Modular Earth Submodel System (MESSy2) to link multiinstitutional computer codes. The physics subroutines of the original ECHAM code have been modularized and reimplemented as MESSy submodels and have been continuously further developed.
For the present study we applied EMAC (MESSy version 2.54.0; Jöckel et al., 2010, 2016) in the T42L90MAresolution, i.e. with a spectral truncation of T42, which corresponds to a quadratic Gaussian grid of ∼ 2.8° × ∼2.8 ° in latitude and longitude and with 90 hybrid pressure levels in the vertical that reach up to 0.01 hPa. The applied model setup comprises the basic submodels for dynamics, radiation, clouds and diagnostics (AEROPT, CLOUD, CLOUDOPT, CVTRANS, E5VDIFF, GWAVE, ORBIT, OROGW, PTRAC, QBO, RAD, SURFACE, TNUDGE, TROPOP and VAXTRA; see Jöckel et al., 2005, 2010, for details) and the submodel SF_{6}. The submodel SF_{6} is used to explicitly calculate the chemical sinks of SF_{6} in the middle atmosphere. The calculations are based on the reaction scheme by Reddmann et al. (2001); details of the submodel are provided in Loeffel et al. (2022), but here we summarize the most important points. For electronattachmentbased SF_{6} degradation, the electron field by Brasseur and Solomon (1986) is used and UV photolysis of SF_{6} is not included as its loss rate is several orders of magnitude weaker than that of electron attachment up to about 100 km. Further reactions considered are photodetachment of SF${}_{\mathrm{6}}^{}$ (Datskos et al., 1995); destruction of SF${}_{\mathrm{6}}^{}$ by atomic hydrogen, hydrogen chloride and ozone (Huey et al., 1995); and stabilization of SF${}_{\mathrm{6}}^{}$ excited by collisions, as well as autodetachment of SF${}_{\mathrm{6}}^{}$. The autodetachment rate was set to 10^{6} s^{−1}, and to compute the photodetachment rate of SF${}_{\mathrm{6}}^{}$, we use the extraterrestrial solar photon flux with no attenuation of the UVphoton flux (see Reddmann et al., 2001), as provided by WMO (1986). The overall global lifetime of SF_{6} in the stratosphere and mesosphere is around 1900 years in the model (Loeffel et al., 2022). In general, the atmospheric SF_{6} lifetime is subject to large uncertainties, with recent estimates ranging from about 800 to 3000 years (Ray et al., 2017; Kovács et al., 2017; Kouznetsov et al., 2020), encompassing the SF_{6} lifetime in the EMAC model.
The simulation setup is designed as a time slice simulation with climate conditions of the year 2000. This means that climatologies of the 1995–2004 period of the greenhouse gases (GHGs) CO_{2}, CH_{4}, N_{2}O and O_{3}; the sea surface temperatures (SSTs); and the sea ice concentrations (SICs) are prescribed as monthly means. The Earth System Chemistry Integrated Modelling (ESCiMo) RC1base07 simulation (see Jöckel et al., 2016) data set was taken for the GHGs, and the Hadley Centre Sea Ice and Sea Surface Temperature (HadISST) data set (which was also used for the RC1base07 simulation) was taken for SSTs and SICs. The CCMI1 volcanic aerosol data (for their effect on infrared radiative heating; see Arfeuille et al., 2013; Morgenstern et al., 2017) are prescribed, and the quasibiennial oscillation (QBO) is nudged (see Jöckel et al., 2016). Moreover, the 1995–2004 average of the reactant species for the SF_{6} sinks, namely HCl, H, N_{2}, O_{2}, O(^{3}P) and O_{3} from the RC1base07 simulation, are prescribed. Other than the SF_{6} submodel, no interactive chemistry is activated in the simulations for this study. Note that in contrast to the climate conditions, the SF_{6} lower boundary conditions are prescribed transiently over the period of the simulation with values starting in year 1950. This simulation has also been used in Loeffel et al. (2022) and therein was referred to as the TS2000 simulation. The setup is designed to investigate the effects of the SF_{6} sinks under constant climate conditions and was chosen for this study to better isolate the effects of the SF_{6} sinks without the influence of modeldependent secular trends in circulation and composition. It was shown by Loeffel et al. (2022) that the transient evolution of circulation and composition only has minor influences on the difference in ideal AoA compared to apparent AoA so that a transient simulation would give similar results to the time slice simulation in terms of the sink correction scheme.
We use the same AoA tracers as described in Loeffel et al. (2022) and in particular the linearly increasing ideal AoA tracer and two SF_{6}tracers, one without chemical depletion and one with chemical depletion. AoA is calculated from the tracer mixing ratios as described in Loeffel et al. (2022), following Fritsch et al. (2019). In brief, we correct for the nonlinear increase in the reference time series in the calculation of AoA from the SF_{6} tracers via a polynomial fit method (with the parameter of the moment ratio of 1.0 and the fraction of input of 95 %; see Fritsch et al., 2019, for details). With this method, the passive SF_{6} tracer results in almost identical AoA to the ideal, linearly increasing AoA tracer. In the following sections, we use the two SF_{6} tracers with and without sinks and refer to the AoA from SF_{6} without sink as “ideal” AoA and to AoA from SF_{6} with sink as “apparent” or SF_{6}based AoA.
3.2 Observational data
3.2.1 Aircraft and balloon data
Concurrent measurements of SF_{6} and other, alternative age tracers are necessary to constrain the relationship between ideal and apparent AoA based on observations. Here, we use a number of such measurements available from aircraft and balloon campaigns spanning a period of 21 years (1997 to 2017) taken during different seasons and at different latitudes (see Fig. 2 for an overview).
Leedham Elvidge et al. (2018) described the stratospheric observations of a set of inert trace gases and demonstrated their usefulness as AoA tracers, including a comparison to SF_{6}based AoA. The five species are CF_{4}, C_{2}F_{6}, C_{3}F_{8}, CHF_{3} (HFC23) and C_{2}HF_{5} (HFC125), all of which result in slightly lower AoAs as compared to SF_{6}, with the discrepancy increasing with increasing AoA (see Fig. 2). This data set is based on air samples collected during several highaltitude balloon (labelled B34 and B44 in Fig. 2) and aircraft (labelled OB09, K2010, K2011AL and SC2016 in Fig. 2) campaigns in the tropics, midlatitudes and high latitudes. It was subsequently augmented by AoA data from a further aircraft campaign (SC2017; Adcock et al., 2021).
The balloon data from 1997, 1998 and 2000 come from in situ measurements of SF_{6} and CO_{2} taken during the Observations of the Middle Stratosphere (OMS) and the SAGE III Ozone Loss and Validation Experiment (SOLVE) campaigns (Volk et al., 1997; Andrews et al., 2001; Moore et al., 2003; Ray et al., 2017). The mean ages were calculated differently for each trace gas, although they both use a variation on the lag technique. The sampling rates were also different between the instruments measuring each trace gas so the points shown in Fig. 2 for these flights were chosen based on a time coincidence of ±30 s.
3.2.2 Satellite data
The Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) was an FTIR spectrometer on the Envisat satellite. It scanned the atmosphere from July 2002 to April 2012 in the altitude range between 6 and 70 km. During this period, about 2 million SF_{6} profiles covering the entire globe and measured during the day and night were obtained. The retrieval of SF_{6} profiles from the radiance spectra has been described previously (Stiller et al., 2008, 2012; Haenel et al., 2015). The data version used here has been retrieved following the approach described by Haenel et al. (2015), however, using newer absorption cross sections of SF_{6} (Harrison, 2020) and accounting for a trichlorofluoromethane (CFC11) band in the vicinity of the SF_{6} signature (Harrison, 2018). This SF_{6} data record is publicly available at https://doi.org/10.5445/IR/1000139453 (Stiller et al., 2021).
Moreover, we use AoA deduced from N_{2}O from the GOZCARDS (Global OZone Chemistry And Related trace gas Data records for the Stratosphere; Froidevaux et al., 2015) satellite product, as produced by Linz et al. (2017). The N_{2}O data are available for the years 2004–2012 and are based on observations by the Atmospheric Chemistry Experiment Fouriertransform spectrometer (ACEFTS) for 2004–2010 and by the microwave limb sounder aboard the Aura satellite (Aura MLS) for 2004–2012. The AoA product relies on an empirical relation between CO_{2}based AoA and N_{2}O as found by Andrews et al. (2001), accounting for the tropospheric trend in N_{2}O (Linz et al., 2017).
In this section, we use selfconsistent model data to develop a sink correction scheme based on the concept laid out in Sect. 2. We start with an overview of the simulated relation between ideal and SF_{6}based apparent AoA. Figure 3 shows the difference between ideal AoA and apparent AoA derived from the chemically depleted SF_{6} tracer as a function of ideal AoA. As found by Loeffel et al. (2022), AoA differences due to chemical sinks increase approximately linearly for AoA lower than about 4 years and increase very strongly for older AoA values. The AoA difference increases from the 1960s to the 2020s due to the increasing mixing ratios of SF_{6} (see Fig. 3a). The increase in differences is smaller in the first decades, when the SF_{6} growth rates are small (see Fig. 1), and these differences increase more strongly with the stronger growth rates after ∼1990. The seasonal and latitudinal variations in the AoA differences are highlighted in Fig. 3b for one specific year (2015). The relation of ideal to apparent AoA is generally found to be quite compact for AoA below about 4 years. Variations with latitude and season are found in the ideal AoA values for which the difference starts to strongly increase; the strong increase in the AoA difference is found at lower ideal AoA values for tropical air compared to high latitudes. The contrast of tropical to highlatitude air is discussed in more detail in the next section. Generally, the mostly compact relationship between ideal and apparent AoA is promising for our aim to deduce a correction method that is valid globally.
Our strategy in developing the correction scheme is as follows. The first step is to parameterize the effective lifetime as a function of AoA using different fit functions (exponential or polynomial; see Sect. 4.1). We put particular focus on testing whether the fit coefficients change over time and how strongly they vary with latitude. Note that we have tested all shown diagnostics first with a pair of linearly increasing tracers (one with and one without chemical sink), but since we found very similar results for the nonlinearly (SF_{6}like) increasing tracer, we only show results for the latter, more realistic case. We then apply the correction schemes with fit coefficients from the model to the model data themselves (Sect. 4.2) to quantify to what degree the biases in SF_{6}derived apparent AoA can be eliminated by the sink correction. This serves to decide which fit function performs best and how much the results improve through latitudedependent fit coefficients compared to global uniform fit coefficients. Finally, in Sect. 4.3 it is tested how the limited availability of data in space and time impacts the fit coefficients and the subsequent correction, as will be relevant for the application to observational data.
4.1 Fit coefficients for different fitting methods
We first present two exemplary profiles of apparent AoA and ideal AoA, one in the tropics (10° S–10° N) and one over the south pole (80–90° S), arbitrarily chosen for the year 2000. In Fig. 4, the difference between apparent and ideal AoA is shown as a function of ideal AoA (Fig. 4a; similar to Fig. 3 but restricted to lower AoA values), showing the strong, nonlinear increase in the AoA difference with increasing ideal AoA. For ideal AoA of up to about 4 years, the difference between apparent and ideal AoA is lower for tropical than for polar latitudes. This can be understood from observing a lower fraction of old air in the tropics compared to the high latitudes for a given AoA, i.e. less weight on the tail of the AoA spectrum (see e.g. Ploeger and Birner, 2016). The older air is more likely to be depleted in SF_{6}. For AoA above about 4 years, on the other hand, the difference between apparent and ideal AoA is larger in the tropics compared to at high latitudes. Ideal AoA of 4 years is located at a pressure level of around 7 hPa in the tropics but at a pressure level as low as around 60 hPa in high latitudes. Thus, the tropical air with ideal AoA above about 4 years is located much higher in the atmosphere and is therefore closer to the region where chemical sinks are active. This likely explains the strong increase in the apparent AoA already for lower ideal AoA compared to at high latitudes. Next to the annual mean profiles, shaded regions display the variation over the year. The seasonal cycle in the ideal–apparent AoA relation is very small in the tropics (see lightblue shading in Fig. 4a), whereas it is more pronounced at high latitudes for AoA values above about 4 years. Specifically, the AoA difference for a given ideal AoA value rises strongly in early spring (around November; not shown). This sudden increase in apparent AoA for ideal AoA of around 5 years is due to a sudden rise in the AoA isosurface – while the 5year contour is located at around 30 to 40 hPa in winter, it rises in altitude to above 10 hPa in spring (not shown); thus, we are looking at a different air mass. If we analyse the difference between apparent and ideal AoA on a fixed pressure level (e.g. 20 hPa), this difference is highest in late winter just before the polar vortex breaks down (October) and the strong downwelling resulted in the accumulation of SF_{6}depleted air over winter in the stratosphere. However, apart from in spring at high latitudes, the relationship between ideal and apparent AoA does not have a strong seasonal cycle (see also Fig. 3), and we will consider only annual mean profiles for performing the fits hereafter.
The effective lifetimes for the tropics and high latitudes calculated from the ideal and apparent AoA profiles via Eq. (5) are shown in the right panel of Fig. 4 and reflect the behaviour observed in the AoA differences. For both profiles, the lifetime of SF_{6} in young air is very long, and it drops rapidly for old AoA. For AoA between about 1 and 4 years, the effective lifetime of tropical air is longer as compared to high latitudes, but it drops to shorter effective lifetimes for ages above about 4 years. The values of the effective lifetimes lie around 200 years. As noted in Sect. 2, the effective lifetime is the pathintegrated lifetime averaged across all possible paths and should not be confused with the total middleatmospheric (stratospheric and mesospheric) lifetime of SF_{6}. The latter was estimated to be about 1900 years in the simulation used here (Loeffel et al., 2022; see also Sect. 3.1).
For the sink correction scheme, a fit to parameterize the effective lifetime as a function of AoA is sought. The fit aims at ages between about 1 and 5 years for ideal AoA. Below 1 year, the effects of chemical sinks are small (strong increase in effective lifetime; see Fig. 4b), and for air with an AoA of above 5 years, the relationship between ideal and apparent AoA gets very steep (see also Fig. 3). Therefore, we perform the fits of τ_{eff} for the range of AoA between 1 and 5 years. Tests with variations in this range showed that performing the fits for this interval gives generally the best results in terms of the bias correction of apparent AoA. For the AoA range between 1 and 5 years, τ_{eff} roughly decreases exponentially with AoA. Thus, we test an exponential function as the fit for τ_{eff}(Γ). However, in particular in the tropics, the turning point in the τ_{eff}–AoA relationship around 3 years cannot be captured by an exponential function, and we also test a thirdorder polynomial fit to capture these turning points better. As shown by the fit lines, the exponential fit roughly captures the decrease in lifetimes with AoA, but the thirdorder polynomial is bettersuited to capture the curvature of τ_{eff}. We further tested a fifthorder polynomial fit, which is only slightly superior to the thirdorder fit (not shown in Fig. 4) at the cost of more fitting parameters.
The resulting apparent–ideal AoA differences using the fit functions for the effective lifetime are shown as solid and dashed lines in the left panel of Fig. 4. Note that while the fit was performed only for 1 year < AoA < 5 years, the functions are still evaluated for lower AoA values. For AoA > 5 years, the exponential function is evaluated as well. For the polynomial fit, the effective lifetime approaches zero for high AoA values so that the apparent AoA would go to infinity and switch sign after the zero crossing of the effective lifetime. This is only a problem for ideal AoA above 5 years, i.e. the values above the range the fits were performed for. To avoid this nonphysical behaviour, the apparent AoA for ideal AoA above 5 years is extrapolated linearly instead of using the polynomial fit function (see Sect. 2). The effect of this can be seen in Fig. 3b and leads to a rather more conservative correction for high AoA values.
Consistent with the fits to the effective lifetime, the polynomial fits capture the AoA differences better, in particular for older AoA. The strong increase in the apparent–ideal AoA difference for older air is not fully captured by the exponential fits; however, they do still capture the general increase in the difference between apparent and ideal AoA reasonably well. When using global fit coefficients obtained by fitting τ_{eff} values from all latitudes simultaneously, the general shape of the difference between apparent and ideal AoA can be reproduced well (as shown for the exponential fit in Fig. 4 as black line).
Overall, the deviations between the fit lines for different regions or different fit functions are small (< 0.1 years) for AoA below 5 years, indicating that the specific details of the fit might be of little importance for the AoA correction. A more quantitative assessment of the correction with the fitted τ_{eff} functions for all latitudes is given in the next section.
As discussed in Sect. 2, by including the timedependent function F_{t}(t) the lifetimes ought to be timeindependent, provided that the simplified assumptions made in the derivations hold. Figure 5 shows the time series of the global fit coefficients for the exponential and the thirdorder polynomial fit. Coefficients vary more strongly from the year 1965 to about 1990 for the polynomial and to about the year 2000 for the exponential fit coefficients. After the year 2000, the fit coefficients show lower variability (with variations mostly smaller than 20%). We interpret the general decline in the fit coefficients as spinup effect, following the simulation start in 1950. While we estimate that ideal AoA is spun up sufficiently after about 15 years (i.e. in 1965), the effective lifetime is more sensitive to the long tail of the spectrum (contributing depleted air to the spectrum) and therefore takes longer to equilibrate. This is confirmed by a similar analysis with a linearly increasing SF_{6} tracer, for which the effective lifetime decreases strongly over the first ∼ 20 years and converges to within 10 % of its equilibrium value by the year 2000 (not shown). Note that this spinup effect of the SF_{6} effective lifetimes is not just an internal property of the model but that a similar effect is likely present in the atmosphere given the start of SF_{6} emissions in the ∼1950s. In the following section, we use mean values of the fit coefficients over the years 2000–2019 (period indicated by the dashed line in Fig. 5) for which the values have converged rather well. The deviations in effective lifetime towards the beginning of the simulation might induce errors in the correction, but it is shown in the following section that those are generally small compared to the overall beneficial effect of the sink correction.
4.2 Application of correction schemes to model data
As the next step, we apply the correction scheme with different fits to model data to quantify how well the correction scheme is able to remove biases due to the chemical sinks in SF_{6}based AoA. In the selfconsistent model context, where we have the full information of ideal and apparent AoA, we can test whether the assumptions and fits made in the formulation of the correction scheme hold. In particular, we test how well different fit functions perform and how much the correction improves using latitudinally dependent fit coefficients instead of globally averaged ones.
We first present an example time series at midlatitudes in the middle stratosphere (at 30 hPa and 46° N; see Fig. 6). Consistent with Loeffel et al. (2022), SF_{6}derived apparent AoA increases over time, while the ideal AoA remains constant, as expected in the time slice simulation with constant climate. The difference between apparent AoA and ideal AoA is around 0.5 years in the year 2000 and up to 1 year towards the end of the simulation in 2030 (see Fig. 6). Corrections with exponential and thirdorder polynomial fits are performed, both with either latitudedependent or global fit coefficients. All four methods are able to strongly reduce the difference between apparent and ideal AoA over the entire simulation period to less than about 0.15 years. All correction schemes overcorrect slightly before ∼ the year 2000, likely because of the change in effective lifetime over time in this period due to spinup effects (see Fig. 5). Furthermore, the remaining difference between the corrected AoA and ideal AoA has a seasonal cycle, indicative of seasonal variations in the effective lifetime which are not captured by the annual mean fit coefficients. The correction using an exponential fit performs slightly worse than the correction using the thirdorder polynomial, but differences between the methods are small. There is almost no difference between the corrected AoA using global versus latitudedependent fit coefficients at this location in the midlatitude middle stratosphere.
The gain of applying the correction at each latitude and pressure level is demonstrated in Fig. 7 for corrections with the exponential and the thirdorder polynomial, each of them for global and latitudedependent fit coefficients. For reference, the top panel in Fig. 7 shows the relative difference between apparent and ideal AoA for the uncorrected SF_{6}derived AoA, showing average differences of 2 %–5 % for young air aged 1–2 years, differences between 5 % and 10 % for AoA between 2 and 4 years of age and strongly rising differences for older air. All correction methods are capable of reducing the relative difference below 2 % for AoA younger than 4 years in almost all regions of the stratosphere. The polynomial fit even reduces the relative differences to less than 1 % in almost all regions for ideal AoA of up to 5 years, in particular when using latitudedependent fit coefficients. Above ideal AoA of 5 years, the bias in apparent AoA increases strongly. As can be seen in Fig. 3b, both the corrections based on the exponential and the linearly extrapolated polynomial fit (see discussion in Sect. 4.1) are for most seasons and latitude bands too conservative (i.e. undercorrect). Since the fits are performed only for ideal AoA below 5 years, the correction scheme is not designed to remove this bias for very old air completely.
The reduction in biases in apparent AoA through the sink correction is further quantified in Fig. 8 for different AoA bins. The probability distributions of relative differences are based on monthly zonal mean data for the decade 2000–2009. The relative differences between apparent corrected and ideal AoA are strongly reduced compared to the uncorrected SF_{6}derived AoA. This is true for all correction schemes. The improvement is smallest for AoA below 1 year, where the effects of sinks are small (median deviation of less than 5 %) and for which the correction scheme is not designed. Nevertheless, applying the correction scheme does not degrade AoA for those young air masses.
All correction methods reduce the relative differences to below 1 % in the median for AoA between 1 and 3 years and to below 2 % for AoA up to 5 years. This corresponds to a factor of 10 or more compared to the median of relative differences for the uncorrected AoA. The polynomial fit methods perform slightly better than the exponential fit method in particular for AoA between 3 and 5 years, as is expected from the better ability to fit the effective lifetimes (see Sect. 4.1). There is little gain when moving from a thirdorder to a fifthorder polynomial fit. For AoA outside the range of 1–5 years, the fifthorder polynomial performed even more poorly, likely because outside the fit region higherorder polynomials are poorly constrained. Thus, based on those results we conclude that a thirdorder polynomial fit performs overall the best.
Using latitudedependent fit coefficients versus global fit coefficients leads to slightly larger reductions in the relative AoA differences in particular for AoA of 2–3 years. The better performance of latitudinal fit coefficients stems mostly from tropical latitudes (see Fig. A1), where the effective lifetimes differ from the remaining latitudes in particular for AoA of 2–4 years (see Fig. 4).
For AoA beyond 5 years, relative differences between uncorrected apparent and ideal AoA become very large, exceeding 50 % in the median and reaching up to almost 250 % for extreme cases (note the scaled values for this AoA bin in Fig. 8). While the correction is not designed for this AoA range, it is still capable of reducing the relative difference from ideal AoA substantially. In particular, the polynomial fit reduces the median of relative differences to below 10 %. The latitudedependent fit coefficients perform worse than the global fits for those old air masses, stemming from an undercorrection at high latitudes (see Figs. A1 and 7).
The fit coefficients show substantial changes for earlier time periods (see Fig. 5), but we use the average fit coefficients for 2000–2019 for all times. Despite this, the relative differences between SF_{6}derived and ideal AoA can also be reduced in the 1980s data to below 2 % in the median for all ideal AoA ranges between 1 and 5 years and with all correction methods (not shown). However, the gain compared to uncorrected AoA is generally smaller since the effects of sinks on AoA are weaker in the 1980s data (ranging from 2 % to 9 % for AoA between 1 and 5 years).
The increasing effect of SF_{6} sinks over time induces apparent positive trends in AoA (see Loeffel et al., 2022, and Fig. 9b). Since we use a time slice simulation, trends in ideal AoA are not significant (except a few locations, highlighting the role of natural variability). For the corrected SF_{6}based AoA, trends are strongly reduced compared to the uncorrected apparent AoA trends and are mostly insignificant for AoA values below about 4 years. However, the correction scheme does not completely remove apparent trends; a weak significant trend remains in midlatitudes in the lower stratosphere in particular. Also for AoA values above about 5 years, a significant positive trend remains, in particular for the exponential scheme. This is consistent with the AoA correction not being able to reduce the very strong biases for those very old air masses (see Fig. 7).
Furthermore, we investigate the effects of the sink correction scheme on the difference between tropical and midlatitude AoA profiles. This is a useful diagnostic as it can be linked to the strength of the overturning circulation (e.g. Linz et al., 2017). The latitudinal age difference for the years 2000–2009 is displayed in Fig. 10, and ideal age shows the typical profile with the largest latitudinal differences in ideal AoA in the lower stratosphere around 60 hPa (e.g. Dietmüller et al., 2018). Apparent AoA strongly overestimates the AoA difference (by about 30 % at 60 hPa, increasing to more than 100 % above 10 hPa). The mass flux through a given surface is inversely related to the age difference between upwelling and downwelling regions on that surface (Linz et al., 2016), and so apparent age leads to a significant underestimate of the circulation strength throughout the stratosphere. The corrected AoA is able to remove this overestimation of the AoA difference. In particular, with the P3 correction, the AoA difference is reproduced remarkably well up to about 3 hPa. Interestingly, the latitudedependent fit coefficients are not beneficial for the AoA difference diagnostic but perform as well as the correction schemes with global mean coefficients. This is consistent with small differences between corrected mean AoA with global versus latitudedependent coefficients (of less than 5 %; see Fig. 7).
4.3 Subsampled data
In the previous sections, we examined the skill of the sink correction scheme when all model data are available for estimating the fit coefficients. However, the goal is to apply the sink correction to observational data and ideally constrain the coefficients for the correction based on observations. However, from observations we have much less information on the relation of “ideal” to apparent AoA. Therefore, we test here by how much the fits and the correction degrade if only limited and uncertain data pairs of ideal and apparent AoA are available to perform the fits on.
To do so, we select data from the model simulation representative of a limited number of observational data sets and choose data from 1 selected month and the latitude and height region of observations for each of the samples. We choose here 10 different data samples, representative of the observational data we have available for comparison (see Sect. 3.2 for a description of the data sets). The subsampled data are almost all located in the Northern Hemisphere, mostly in midlatitudes with a few tropical and highlatitude profiles. Furthermore, only four profiles reach or are just above the height of 30 km; i.e. older air is underrepresented. The differences between apparent and ideal AoA for the subsampled data are shown in Fig. 11a. Due to the limited data, only one global fit is performed, since the data are not sufficient for latitudedependent fits. However, as shown in the last section, using global instead of latitudedependent fit coefficients does not degrade the correction significantly.
We perform a bootstrap procedure to estimate the uncertainty in the fit coefficients. Specifically, we choose 141 samples from the 141 available data points with replacements and add a random error drawn from a normal distribution with standard deviation of 0.1 years to each apparent and ideal AoA. This error is roughly representative of the uncertainty in observational data (see Sect. 3.2). We repeat this procedure 1000 times and calculate the average and standard deviations of the fit coefficients from this suite.
As seen in Fig. 11, the resulting fit line for the exponential fit to the subsampled data (magenta) lies below the fit based on the full model data (black). The fit varies moderately within 2 standard deviations, resulting in a correction at 5year AoA spanning between about 0.65 and 0.95 years. The thirdorder polynomial fit, on the other hand, is very unconstrained for the subsampled data (dotted blue lines). The 2standarddeviation range spans essentially the whole domain; therefore, only a few example fits are shown to demonstrate the unconstrained fit.
The effects of the uncertainty in the fit coefficients on the correction of AoA are summarized in the error statistics in Fig. 12. The correction was performed for fit coefficients plus and minus 2 standard deviations from its best estimate for the fits to the subsampled data. The probability distributions show relative differences including the corrections with both the fit coefficients plus and minus 2 standard deviations. The relative differences between ideal AoA and corrected AoA with the thirdorder polynomial fit are very large and span a wide range, which is consistent with the unconstrained fit lines found above. Thus, we conclude that correction with a thirdorder polynomial fit based on a limited amount of data and/or uncertain data is not advisable, as the coefficients are not well constrained. The exponential fit, on the other hand, is much more robust. The relative difference between corrected and ideal AoA is larger when using the exponential fit coefficients obtained from the subsampled data compared to the full model data for all AoA bins. Nevertheless, the correction still reduces the error compared to uncorrected data substantially (to below 5 % for AoA between 1–4 years and below 10 % for AoA between 4–5 years). In conclusion, despite the fact that limited information of the relation of SF_{6}derived and ideal AoA leads to degraded fits, the correction scheme based on the global exponential fit is still very beneficial compared to using uncorrected SF_{6}based AoA.
In Sect. 4, the performance of the correction scheme was carefully evaluated based on model data, and it was found that the schemes are highly beneficial to remove biases in SF_{6}based AoA estimates. However, this was entirely based on model data, and one can question whether the chemical sink of SF_{6} is constrained well enough in the model to trust the derived coefficients. Ideally, we would perform similar fits between AoA derived from observed SF_{6} and from simultaneous observations of other age tracers that do not suffer from the sinks. However, as we will show in the following section, available observational data are currently not able to constrain the fit parameters well enough. Instead, we will show that the modelderived correction fit lines lie within the bounds given by observations, and we will test the performance of the correction with modelbased parameters on independent AoA observations.
We use the observational data sets described in Sect. 3.2 from a number of aircraft campaigns and balloon flights. Figure 13a shows differences between SF_{6}based AoA and AoA based on the alternative tracers as a function of the latter, which is similar to what is shown in Fig. 4a but for observational data. It is immediately clear from this figure that observations show a large spread in the difference between SF_{6} and alternativetracerbased AoA. In the left panel, the difference is compared to model data sampled in the same month and years as when the observational campaigns took place. The model data are still monthly and zonal mean profiles, explaining the muchreduced variability compared to observational data. Furthermore, note that the meteorological conditions in the freerunning model do not of course match the observed state of the atmosphere, so the sampling is rather meant to capture variations due to the latitude, season and the underlying trend in SF_{6}. Generally, the model data are found to lie within the large scatter of observational data for mean AoA below about 5 years. One highlatitude balloon flight sampled air within the polar vortex (Ray et al., 2017) with CO_{2}based mean AoA between 5 and 6 years. For this profile, the difference in CO_{2} and SF_{6}based AoA increases strongly to up to 4 years. The AoA differences in the model from a similar latitude band and month as this balloon flight show a similarly strong increase but occurring already at younger ideal mean AoA values.
The measurement campaigns took place during different years within a 21year period (between 1997 and 2017), so the effect of the sinks increased over this time frame. To eliminate this time dependence, the age differences are normalized by the timecorrection function F_{t} (see Sect. 2), and the resulting data points are shown in Fig. 13b. The scatter is somewhat reduced compared to the uncorrected differences on the left but is still substantial. Attempts to fit the parameters for the correction scheme result in a very large uncertainty. In particular for young AoA values below 2 years, some measurements indicate already a large difference from SF_{6}based AoA ranging between 0 and 1 year (see Fig. 13a). This might partly be due to local influence close to sources of the alternative AoA tracers, potentially leading to biases in the estimate of ideal AoA from alternative tracers and/or SF_{6}based AoA (Adcock et al., 2021). In Sect. 4.3, we find that the subsampled model data on the availability of observational data are not able to constrain the thirdorder polynomial fit. Therefore, we only use the exponential fit method in the following section, both for attempts to fit coefficients based on observations and for any corrections of the observational data. When discarding all mean AoA values below 2 years and above 5 years, the exponential fit for the correction scheme as displayed in grey in Fig. 13b is obtained (based on a bootstrapping algorithm to estimate the uncertainty; see the grey shading). However, the fit is very sensitive to including or excluding individual data points. When including data for values between 1 and 2 years, the resulting fit would indicate an increase in the effective lifetime with mean AoA, which can be ruled out as nonphysical. Therefore, we have to conclude that the currently available observational data that allow for an estimation of SF_{6}based AoA and AoA from an alternative tracer are not sufficient in quantity and quality to constrain the fit parameters for the sink correction.
Instead, we have to rely on the modelbased parameters for the correction scheme. The modelbased global mean correction curve for the exponential fit is included in Fig. 13b (black curve). While the large scatter of observations does not allow for a detailed evaluation, we find that the modelbased correction scheme lies within the observational uncertainty. Therefore, we suggest at this stage that the best practice is to apply the sink correction with the modelbased parameters. While fit functions which are of higher order than the exponential fit and the latitudedependent fit coefficient did lead to a slightly better performance of the correction within the model (see Sect. 4.2), we rather use the more conservative exponential fit and the global mean fit parameters in the following to correct observational data. This is, on the one hand, because of the larger sensitivity in the polynomial fits (see Sect. 4.3) and, on the other hand, because we want to reduce modeldependent information. In particular, the latitude dependence of fit coefficients likely is a function of the representation of circulation in the model (e.g. location and strength of the polar vortex), while we assume the global fit to be less sensitive to the representation of circulation in a particular model. Comparisons to other independent observational data and intercomparisons with other models with the implementation of chemical sinks of SF_{6} should be performed in the future to better constrain the fit coefficients of the sink correction.
The correction scheme with the modelderived parameters is applied to in situ and satellite data presented in Fig. 14. The latitudinal profiles of in situ mean AoA data compiled by Andrews et al. (2001) at around 20 km consist of AoA both based on CO_{2} and based on SF_{6}. The SF_{6}based mean AoA data consistently lie above AoA based on CO_{2}, in particular for mid to high latitudes. This difference is not necessarily caused only due to the chemical SF_{6} sink, but also by the different temporal slopes and seasonality of the tracers time series, affecting the derivation of mean AoA (see e.g. Andrews et al., 2001; Bönisch et al., 2009).
When applying the sink correction to the SF_{6}based mean AoA values, this consistent high bias of SF_{6}based mean AoA is strongly reduced. In particular at northern mid to high latitudes, CO_{2} and SF_{6}based AoA agree extremely well after the sink correction is applied, with the mean difference between CO_{2} and SF_{6} AoA for latitudes of > 40° being reduced from 0.43 to 0.03 years. In the Southern Hemisphere, the scatter is larger due to generally poorer data quality. Nevertheless, the bias between CO_{2} and SF_{6}based AoA is also reduced from 0.33 to 0.13 years for latitudes south of 40° S.
Furthermore, the latitudinal climatological profile of AoA around 20 km is shown for two satellite products, SF_{6}based mean AoA from MIPAS (Stiller et al., 2021, and see Sect. 3.2) and AoA derived from N_{2}O measured by Aura MLS and ACEFTS (the GOZCARDS data product; see Sect. 3.2). The SF_{6}based mean AoA from MIPAS is older compared to AoA based on N_{2}O and in situ CO_{2} observations, in particular at high latitudes, but agrees well overall with the in situ SF_{6}based mean AoA. After applying the sink correction to the MIPAS AoA data, agreement between MIPAS AoA and the nonSF_{6}based AoA data sets is much improved, in particular at high latitudes. At midlatitudes, corrected MIPAS mean AoA is slightly lower than the other estimates, while in the tropics the high bias is only slightly reduced.
Overall, mean AoA estimates from different trace gases and from different observational sources agree much better after applying the SF_{6} sink correction. This provides further evidence that the sink correction is a promising way forward to remove biases in mean AoA estimated based on SF_{6}.
In this paper, we have developed a correction scheme for SF_{6}based AoA to remove the biases generated by chemical SF_{6} sinks. The correction scheme is based on a simplified formulation derived from the theory that provides a relation of ideal to apparent AoA dependent on (1) a timedependent function of the reference SF_{6} surface mixing ratio (which would be a linear function for a linear reference time series) and (2) the pathintegrated average effective lifetime based on AoA. For the former, we suggest a simplified approximation based on the linearized reference time series. For the latter, we parameterize the effective lifetime as a function of ideal AoA. Based on global model data, we find that either an exponential or a thirdorder polynomial fit is suited for this parameterization. While the thirdorder polynomial fit captures the details of the AoA dependency of the effective lifetime better, it comes at the cost of high parameter sensitivity to uncertainty in the underlying data. Furthermore, we find that using latitudedependent coefficients for the parameterization of the effective lifetime as a function of ideal AoA is beneficial, but the overall ability to reduce biases in SF_{6}based AoA is only marginally improved. In the model, the simplified sink correction with global mean fit coefficients is able to reduce the biases in AoA by a factor of 5 to 10 to, on average, less than 2 % relative to ideal AoA for AoA values below 5 years. For AoA above 5 years, the relation between ideal and apparent AoA becomes weaker, and, therefore, the correction is not designed for those high ages above 5 years. Nevertheless, correcting AoA with the extrapolated fit lines for this AoA range also substantially reduced the bias. In addition to mean AoA, we showed that the correction scheme is able to remove apparent AoA trends in the lower stratosphere. Furthermore, the sink correction is able to correct the age difference between tropics and midlatitudes, a diagnostic of the residual circulation strength, even when based on global mean fit coefficients. This is a promising result to enable more confident estimates of the circulation from SF_{6} satellite data, which previously was limited only to the lower stratosphere (Linz et al., 2017).
In conclusion, the sink correction scheme performs very well within the selfconsistent “model world” with all information available. However, how useful is this scheme to correcting observational AoA estimates given imperfect knowledge of the ideal to apparent AoA relationship? Firstly, we approach this question by subsampling the model data leading to erroneous fit coefficients. When using the less sensitive exponential fit function for the relation of the effective lifetime to AoA, we could show that despite incomplete knowledge and thus uncertain fit coefficients, the sink correction is able to reduce the AoA biases in all cases. In other words, even applying the sink correction with uncertain parameters is better than not applying any corrections. Of course this statement hinges on staying within a certain range of fit uncertainty; if the sink in the model were completely different than the real world, this would naturally not hold. The total stratospheric and mesospheric lifetimes of SF_{6} estimated in the EMAC model simulations spanning 1900 years fall within the range of values provided by previous model and observational studies (e.g. Ray et al., 2017; Kovács et al., 2017; Kouznetsov et al., 2020, all providing ranges between around 800 and 3000 years), indicating the validity of EMAC's SF_{6} depletion mechanisms despite existing uncertainties in lifetime determination (see also discussion in Loeffel et al., 2022). Furthermore, comparisons of the model fit to available observational data of SF_{6}based AoA and AoA deduced from CO_{2} or the new alternative AoA tracers (Leedham Elvidge et al., 2018) indicate that the model sinks do lie within the observational spread. The spread is high, however, likely because of two reasons. Firstly, AoA deduced from realistic tracers is generally errorprone due to uncertainty in the reference time series and the nonlinearity therein (including the seasonal cycle, in particular for CO_{2}) and, for example, the influence of local sources of the tracers. Those factors induce the large error bars in observational AoA (seen in Fig. 13). Secondly, the observational data used here are point measurements, while in the model we used monthly mean and zonal mean data. While the averaged data fall on the compact relationship between ideal AoA and apparent AoA, point data might differ from this line due to local mixing between different air masses (as is wellknown for tracer–tracer relationships; see e.g. Plumb, 2002). It is not clear how strong the contribution of those two factors is to the large spread in the observed relationship between SF_{6}based AoA and AoA from other tracers. If the contribution from internal variability by mixing is substantial, this would imply that the compactness of the relationship between SF_{6}AoA and ideal AoA is limited for point measurements, and thus a correction is not easily possible. Until this issue is clarified, we recommend using the sink correction scheme on averaged data (over time periods of a month or so and/or large regions). Application of the sink correction to satellite and compiled in situ data confirms that the correction is able to reduce the discrepancy between different AoA estimates. In situ SF_{6} and CO_{2}based AoA values agree very well after the correction is applied, and the SF_{6}based AoA deduced from MIPAS satellite observations agrees well with the N_{2}Obased AoA from the merged GOZCARDS satellite product when the correction is applied. Note, however, that AoA derived from N_{2}O also bears uncertainties; AoA is calculated from the N_{2}O–AoA relationship based on specific concurrent observations of N_{2}O and CO_{2} in the northern midlatitudes (see Sect. 3.2), and improved quantification of the latitudinal and temporal variations in the N_{2}O–AoA relation will be necessary to more reliably derive mean AoA by this method.
Overall, we conclude that observational data are currently not able to constrain the parameters for the sink correction scheme proposed here. Until this is possible, we recommend using the modelbased global mean exponential fit parameters for the sink correction. While the model might be subject to biases, we could demonstrate that the relation of ideal to apparent AoA from the model lies within observational uncertainty and that the application of the modelbased sink correction is successfully reducing differences in observational AoA estimates. However, future work is needed to test and better constrain the parameters for the sink correction scheme and to generally enhance our knowledge on the chemical sinks of SF_{6} in the atmosphere.
The Modular Earth Submodel System (MESSy) is continuously further developed and applied by a consortium of institutions. The usage of MESSy and access to the source code are licensed to all affiliates of institutions which are members of the MESSy Consortium. Institutions can become members of the MESSy Consortium by signing the MESSy Memorandum of Understanding. More information can be found on the MESSy Consortium website (http://www.messyinterface.org, last access: 25 March 2024). The simulation presented here has been carried out using MESSy version 2.54.0. The Python code to perform the sink correction is provided as a Supplement to this paper. The ECHAM/MESSy Atmospheric Chemistry (EMAC) simulation data are archived at the German Climate Computation Centre (DKRZ) and are available upon request to the authors; likewise, observational data are available upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/acp2441932024supplement.
HG designed the study, performed the analysis and wrote the paper. RE performed the model simulations. JL, ER, GS, HB, LS and ML provided observational data and contributed to the design of the sink correction. All authors contributed to improving the paper.
At least one of the (co)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peerreview process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
This paper arises as the outcome of the International Space Science Institute (ISSI) project “Stratospheric AgeofAir: Reconciling Observations and Models”. We thank ISSI for the support in the form of two team meetings held in Bern. The authors thank all team members for stimulating discussion and comments on the work. We thank Matthias Nuetzel for comments on the manuscript and the two anonymous reviewers for constructive comments that helped to improve the manuscript.
Roland Eichinger received support from the Czech Science Foundation (grant no. 2103295S).
The article processing charges for this openaccess publication were covered by the German Aerospace Center (DLR).
This paper was edited by John Plane and reviewed by two anonymous referees.
Adcock, K., Fraser, P., Hall, B., Langenfelds, R., Lee, G., Montzka, S., Oram, D., Röckmann, T., Stroh, F., Sturges, W., Vogel, B., and Laube, J.: AircraftBased Observations of OzoneDepleting Substances in the Upper Troposphere and Lower Stratosphere in and Above the Asian Summer Monsoon, J. Geophys. Res.Atmos., 126, e2020JD033137, https://doi.org/10.1029/2020JD033137, 2021. a, b
Andrews, A. E., Boering, K. A., Daube, B. C., Wofsy, S. C., Loewenstein, M., Jost, H., Podolske, J. R., Webster, C. R., Herman, R. L., Scott, D. C., Flesch, G. J., Moyer, E. J., Elkins, J. W., Dutton, G. S., Hurst, D. F., Moore, F. L., Ray, E. A., Romashkin, P. A., and Strahan, S. E.: Mean ages of stratospheric air derived from in situ observations of CO_{2}, CH_{4}, and N_{2}O, J. Geophys. Res.Atmos., 106, 32295–32314, https://doi.org/10.1029/2001JD000465, 2001. a, b, c, d, e, f, g
Arfeuille, F., Luo, B. P., Heckendorn, P., Weisenstein, D., Sheng, J. X., Rozanov, E., Schraner, M., Brönnimann, S., Thomason, L. W., and Peter, T.: Modeling the stratospheric warming following the Mt. Pinatubo eruption: uncertainties in aerosol extinctions, Atmos. Chem. Phys., 13, 11221–11234, https://doi.org/10.5194/acp13112212013, 2013. a
Bönisch, H., Engel, A., Curtius, J., Birner, Th., and Hoor, P.: Quantifying transport into the lowermost stratosphere using simultaneous insitu measurements of SF_{6} and CO_{2}, Atmos. Chem. Phys., 9, 5905–5919, https://doi.org/10.5194/acp959052009, 2009. a
Brasseur, G. and Solomon, S.: Aeronomy of the middle atmosphere, Reidel, Dordrecht, the Netherlands, 2 edn., ISBN 9027723443, 1986. a
Brewer, A. W.: Evidence for a world circulation provided by the measurements of helium and water vapor distribution in the stratosphere, Q. J. Roy. Meteor. Soc., 75, 351–363, 1949. a
Butchart, N.: The BrewerDobson circulation, Rev. Geophys., 52, 157–184 https://doi.org/10.1002/2013RG000448, 2014. a
Datskos, P., Carter, J., and Christophorou, L.: Photodetachment of SF_{6}, Chem. Phys. Lett., 239, 37–43, https://doi.org/10.1016/00092614(95)004173, 1995. a
Dietmüller, S., Eichinger, R., Garny, H., Birner, T., Boenisch, H., Pitari, G., Mancini, E., Visioni, D., Stenke, A., Revell, L., Rozanov, E., Plummer, D. A., Scinocca, J., Jöckel, P., Oman, L., Deushi, M., Kiyotaka, S., Kinnison, D. E., Garcia, R., Morgenstern, O., Zeng, G., Stone, K. A., and Schofield, R.: Quantifying the effect of mixing on the mean age of air in CCMVal2 and CCMI1 models, Atmos. Chem. Phys., 18, 6699–6720, https://doi.org/10.5194/acp1866992018, 2018. a
Dobson, G. M. B.: Origin and Distribution of the Polyatomic Molecules in the Atmosphere, P. Roy. Soc. Lond. A Mat., 236, 187–193, 1956. a
Engel, A., Mobius, T., Bonisch, H., Schmidt, U., Heinz, R., Levin, I., Atlas, E., Aoki, S., Nakazawa, T., Sugawara, S., Moore, F., Hurst, D., Elkins, J., Schauffler, S., Andrews, A., and Boering, K.: Age of stratospheric air unchanged within uncertainties over the past 30 years, Nat. Geosci., 2, 28–31, https://doi.org/10.1038/ngeo388, 2009. a, b, c
Fritsch, F., Garny, H., Engel, A., Bönisch, H., and Eichinger, R.: Sensitivity of age of air trends to the derivation method for nonlinear increasing inert SF_{6}, Atmos. Chem. Phys., 20, 8709–8725, https://doi.org/10.5194/acp2087092020, 2020. a, b, c, d, e
Froidevaux, L., Anderson, J., Wang, H.J., Fuller, R. A., Schwartz, M. J., Santee, M. L., Livesey, N. J., Pumphrey, H. C., Bernath, P. F., Russell III, J. M., and McCormick, M. P.: Global OZone Chemistry And Related trace gas Data records for the Stratosphere (GOZCARDS): methodology and sample results with a focus on HCl, H_{2}O, and O_{3}, Atmos. Chem. Phys., 15, 10471–10507, https://doi.org/10.5194/acp15104712015, 2015. a
Haenel, F. J., Stiller, G. P., von Clarmann, T., Funke, B., Eckert, E., Glatthor, N., Grabowski, U., Kellmann, S., Kiefer, M., Linden, A., and Reddmann, T.: Reassessment of MIPAS age of air trends and variability, Atmos. Chem. Phys., 15, 13161–13176, https://doi.org/10.5194/acp15131612015, 2015. a, b, c
Harrison, J. J.: New and improved infrared absorption cross sections for trichlorofluoromethane (CFC11), Atmos. Meas. Tech., 11, 5827–5836, https://doi.org/10.5194/amt1158272018, 2018. a
Harrison, J. J.: New infrared absorption cross sections for the infrared limb sounding of sulfur hexafluoride (SF_{6}), J. Quant. Spectrosc. Ra., 260, 107432, https://doi.org/10.1016/j.jqsrt.2020.107432, 2020. a
Haynes, P., Marks, C., McIntyre, M., Shepherd, T., and Shine, K.: On the “downward control” of extratropical diabatic circulations by eddyinduced mean zonal forces, J. Atmos. Sci., 48, 651–678, 1991. a
Holzer, M. and Waugh, D. W.: Interhemispheric transit time distributions and pathdependent lifetimes constrained by measurements of SF_{6}, CFCs, and CFC replacements, Geophys. Res. Lett., 42, 4581–4589, https://doi.org/10.1002/2015GL064172, 2015. a
Huey, L. G., Hanson, D. R., and Howard, C. J.: Reactions of SF${}_{\mathrm{6}}^{}$ and I− with Atmospheric Trace Gases, J. Phys. Chem., 99, 5001–5008, https://doi.org/10.1021/j100014a021, 1995. a
Jöckel, P., Sander, R., Kerkweg, A., Tost, H., and Lelieveld, J.: Technical Note: The Modular Earth Submodel System (MESSy) – a new approach towards Earth System Modeling, Atmos. Chem. Phys., 5, 433–444, https://doi.org/10.5194/acp54332005, 2005. a
Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H., Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717–752, https://doi.org/10.5194/gmd37172010, 2010. a, b, c
Jöckel, P., Tost, H., Pozzer, A., Kunze, M., Kirner, O., Brenninkmeijer, C. A. M., Brinkop, S., Cai, D. S., Dyroff, C., Eckstein, J., Frank, F., Garny, H., Gottschaldt, K.D., Graf, P., Grewe, V., Kerkweg, A., Kern, B., Matthes, S., Mertens, M., Meul, S., Neumaier, M., Nützel, M., OberländerHayn, S., Ruhnke, R., Runde, T., Sander, R., Scharffe, D., and Zahn, A.: Earth System Chemistry integrated Modelling (ESCiMo) with the Modular Earth Submodel System (MESSy) version 2.51, Geosci. Model Dev., 9, 1153–1200, https://doi.org/10.5194/gmd911532016, 2016. a, b, c, d
Kouznetsov, R., Sofiev, M., Vira, J., and Stiller, G.: Simulating age of air and the distribution of SF_{6} in the stratosphere with the SILAM model, Atmos. Chem. Phys., 20, 5837–5859, https://doi.org/10.5194/acp2058372020, 2020. a, b
Kovács, T., Feng, W., Totterdill, A., Plane, J. M. C., Dhomse, S., GómezMartín, J. C., Stiller, G. P., Haenel, F. J., Smith, C., Forster, P. M., García, R. R., Marsh, D. R., and Chipperfield, M. P.: Determination of the atmospheric lifetime and global warming potential of sulfur hexafluoride using a threedimensional model, Atmos. Chem. Phys., 17, 883–898, https://doi.org/10.5194/acp178832017, 2017. a, b, c
Leedham Elvidge, E. C., Bönisch, H., Brenninkmeijer, C. A. M., Engel, A., Fraser, P. J., Gallacher, E., Langenfelds, R., Mühle, J., Oram, D. E., Ray, E. A., Ridley, A. R., Röckmann, T., Sturges, W. T., Weiss, R. F., and Laube, J. C.: Evaluation of stratospheric age of air from CF_{4}, C_{2}F_{6}, C_{3}F_{8}, CHF_{3}, HFC125, HFC227ea and SF_{6}; implications for the calculations of halocarbon lifetimes, fractional release factors and ozone depletion potentials, Atmos. Chem. Phys., 18, 3369–3385, https://doi.org/10.5194/acp1833692018, 2018. a, b, c
Linz, M., Plumb, R. A., Gerber, E. P., and Sheshadri, A.: The Relationship between Age of Air and the Diabatic Circulation of the Stratosphere, J. Atmos. Sci., 73, 4507–4518, https://doi.org/10.1175/JASD160125.1, 2016. a
Linz, M., Plumb, R., Gerber, E., Haenel, F., Stiller, S., Kinnison, D. E., Ming, A., and Neu, J. L.: The strength of the meridional overturning circulation of the stratosphere, Nat. Geosci., 10, 663–667, https://doi.org/10.1038/ngeo3013, 2017. a, b, c, d, e
Loeffel, S., Eichinger, R., Garny, H., Reddmann, T., Fritsch, F., Versick, S., Stiller, G., and Haenel, F.: The impact of sulfur hexafluoride (SF_{6}) sinks on age of air climatologies and trends, Atmos. Chem. Phys., 22, 1175–1193, https://doi.org/10.5194/acp2211752022, 2022. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o
Moore, F. L., Elkins, J. W., Ray, E. A., Dutton, G. S., Dunn, R. E., Fahey, D. W., McLaughlin, R. J., Thompson, T. L., Romashkin, P. A., Hurst, D. F., and Wamsley, P. R.: Balloonborne in situ gas chromatograph for measurements in the troposphere and stratosphere, J. Geophys. Res., 108, 8330, https://doi.org/10.1029/2001JD000891, 2003. a
Morgenstern, O., Hegglin, M. I., Rozanov, E., O'Connor, F. M., Abraham, N. L., Akiyoshi, H., Archibald, A. T., Bekki, S., Butchart, N., Chipperfield, M. P., Deushi, M., Dhomse, S. S., Garcia, R. R., Hardiman, S. C., Horowitz, L. W., Jöckel, P., Josse, B., Kinnison, D., Lin, M., Mancini, E., Manyin, M. E., Marchand, M., Marécal, V., Michou, M., Oman, L. D., Pitari, G., Plummer, D. A., Revell, L. E., SaintMartin, D., Schofield, R., Stenke, A., Stone, K., Sudo, K., Tanaka, T. Y., Tilmes, S., Yamashita, Y., Yoshida, K., and Zeng, G.: Review of the global models used within phase 1 of the Chemistry–Climate Model Initiative (CCMI), Geosci. Model Dev., 10, 639–671, https://doi.org/10.5194/gmd106392017, 2017. a
Ploeger, F. and Birner, T.: Seasonal and interannual variability of lower stratospheric age of air spectra, Atmos. Chem. Phys., 16, 10195–10213, https://doi.org/10.5194/acp16101952016, 2016. a
Plumb, R.: Stratospheric transport, J. Meteorol. Soc. Jpn., 80, 793–809, 2002. a
Ray, E. A., Moore, F. L., Elkins, J. W., Rosenlof, K. H., Laube, J. C., Röckmann, T., Marsh, D. R., and Andrews, A. E.: Quantification of the SF_{6} lifetime based on mesospheric loss measured in the stratospheric polar vortex, J. Geophys. Res.Atmos., 122, 4626–4638, https://doi.org/10.1002/2016JD026198, 2017. a, b, c, d, e
Reddmann, T., Ruhnke, R., and Kouker, W.: Threedimensional model simulations of SF_{6} with mesospheric chemistry, J. Geophys. Res., 106, 14525–14537, https://doi.org/10.1029/2000JD900700, 2001. a, b
Roeckner, E., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kornblueh, L., Manzini, E., Schlese, U., and Schulzweida, U.: Sensitivity of simulated climate to horizontal and vertical resolution in the ECHAM5 atmosphere model, J. Climate, 19, 3771–3791, https://doi.org/10.1175/JCLI3824.1, 2006. a
Schoeberl, M. R., Sparling, L. C., Jackman, C. H., and Fleming, E. L.: A Lagrangian view of stratospheric trace gas distributions, J. Geophys. Res.Atmos., 105, 1537–1552, https://doi.org/10.1029/1999JD900787, 2000. a
Stiller, G., Harrison, J. J., Haenel, F., Glatthor, N., and Kellmann, S.: MIPAS monthly zonal mean Age of Air data, based on MIPAS SF_{6} retrievals, KITIMK [data set], https://doi.org/10.5445/IR/1000139453, 2021. a, b, c
Stiller, G. P., von Clarmann, T., Höpfner, M., Glatthor, N., Grabowski, U., Kellmann, S., Kleinert, A., Linden, A., Milz, M., Reddmann, T., Steck, T., Fischer, H., Funke, B., LópezPuertas, M., and Engel, A.: Global distribution of mean age of stratospheric air from MIPAS SF_{6} measurements, Atmos. Chem. Phys., 8, 677–695, https://doi.org/10.5194/acp86772008, 2008. a
Stiller, G. P., von Clarmann, T., Haenel, F., Funke, B., Glatthor, N., Grabowski, U., Kellmann, S., Kiefer, M., Linden, A., Lossow, S., and LópezPuertas, M.: Observed temporal evolution of global mean age of stratospheric air for the 2002 to 2010 period, Atmos. Chem. Phys., 12, 3311–3331, https://doi.org/10.5194/acp1233112012, 2012. a, b, c, d, e
Volk, C., Elkins, J., Fahey, D., Dutton, G., Gilligan, J., Loewenstein, M., Podolske, J., Chan, K., and Gunson, M.: Evaluation of source gas lifetimes from stratospheric observations, J. Geophys. Res., 102, 25543–25564, 1997. a
Waugh, D. and Hall, T. M.: Age of stratospheric air: theory, observations, and models, Rev. Geophys., 40, 1010, https://doi.org/10.1029/2000RG000101, 2002. a
WMO: Atmospheric ozone 1985 – volume I–III: assessment of our understanding of the processes controlling its present distribution and change, World Meteorological Organization, https://library.wmo.int/idurl/4/60466 (last access: 25 March 2024), 1986. a
 Abstract
 Introduction
 Theoretical considerations for sink correction methods
 Data and model description
 Testing the correction scheme in a global model
 Comparison and application to observational data
 Summary and conclusion
 Appendix A: Error statistics for tropics and high latitudes
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Supplement
 Abstract
 Introduction
 Theoretical considerations for sink correction methods
 Data and model description
 Testing the correction scheme in a global model
 Comparison and application to observational data
 Summary and conclusion
 Appendix A: Error statistics for tropics and high latitudes
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Supplement