The impact of sulfur hexafluoride (SF6) sinks on age of air climatologies and trends
Mean age of air (AoA) is a common diagnostic for the strength of the stratospheric overturning circulation in both climate models and observations. AoA climatologies and AoA trends over the recent decades of model simulations and proxies derived from observations of long-lived tracers do not agree. Satellite observations show much older air than climate models, and while most models compute a clear decrease in AoA over the last decades, a 30-year time series from measurements shows a statistically nonsignificant positive trend in the Northern Hemisphere extratropical middle stratosphere. Measurement-based AoA derivations are often founded on observations of the trace gas sulfur hexafluoride (SF6), a fairly long-lived gas with a near-linear increase in emissions during recent decades. However, SF6 has chemical sinks in the mesosphere that are not considered in most model studies. In this study, we explicitly compute the chemical SF6 sinks based on chemical processes in the global chemistry climate model EMAC (ECHAM/MESSy Atmospheric Chemistry). We show that good agreement between stratospheric AoA in EMAC and MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) is reached through the inclusion of chemical SF6 sinks, as these sinks lead to a strong increase in the stratospheric AoA and, therefore, to a better agreement with MIPAS satellite observations. Remaining larger differences at high latitudes are addressed, and possible reasons for these differences are discussed. Subsequently, we demonstrate that the AoA trends are also strongly influenced by the chemical SF6 sinks. Under consideration of the SF6 sinks, the AoA trends over the recent decades reverse sign from negative to positive. We conduct sensitivity simulations which reveal that this sign reversal does not result from trends in the stratospheric circulation strength nor from changes in the strength of the SF6 sinks. We illustrate that even a constant SF6 destruction rate causes a positive trend in the derived AoA, as the amount of depleted SF6 scales with increasing SF6 abundance itself. In our simulations, this effect overcompensates for the impact of the accelerating stratospheric circulation which naturally decreases AoA. Although various sources of uncertainties cannot be quantified in detail in this study, our results suggest that the inclusion of SF6 depletion in models has the potential to reconcile the AoA trends of models and observations. We conclude the study with a first approach towards a correction to account for SF6 loss and deduce that a linear correction might be applicable to values of AoA of up to 4 years.
The Brewer–Dobson circulation (BDC) describes the stratospheric transport circulation, consisting of the mean overturning circulation of air ascending in the tropical pipe, moving poleward and descending in the extratropics (Brewer, 1949; Dobson and Massey, 1956), as well as isentropic mixing. A good measure to diagnose this transport circulation is the age of stratospheric air (AoA), which is defined as the mean transport time of an air parcel from its entry into the stratosphere (or from the surface) to any point therein (Hall and Plumb, 1994; Waugh and Hall, 2002). In general circulation models (GCMs), AoA is commonly represented by an inert tracer with a strictly linear temporally increasing surface mixing ratio and is calculated as the corresponding time lag between the local mixing ratio and the mixing ratio at a reference point (Hall and Plumb, 1994). The same method can be applied to real long-lived tracers with a linear trend in tropospheric concentration, and AoA has been derived, for example, from balloon-borne in situ measurements of sulfur hexafluoride (SF6) (Andrews et al., 2001; Engel et al., 2009, 2017). This trace gas is particularly suitable for these studies, as it is stable in the troposphere and stratosphere and its tropospheric concentrations have increased nearly linearly over recent decades. Along with observations of other trace gases, these measurements form a long-term dataset of observationally based AoA restricted to the Northern Hemisphere midlatitudes that covers more than 40 years. A near-global dataset of AoA covering 10 years from 2002 to 2012 was derived in Stiller et al. (2012), Haenel et al. (2015), and Stiller et al. (2017), who retrieved SF6 distributions from MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) satellite observations, but these cover a much shorter time period.
Observations and model simulations of AoA often disagree. AoA derived from observations is mainly older than simulated AoA (see, e.g., SPARC, 2010, and Dietmüller et al., 2018), and the AoA trend over recent decades even differs in sign between observations and models. While most climate models show a clear decrease in AoA over time (see, e.g., Butchart and Scaife, 2001, Garcia et al., 2011, and Eichinger et al., 2019), consistent with the simulated acceleration of the BDC in the course of climate change (see, e.g., Garcia and Randel, 2008), the time series of the observations presented in the studies by Engel et al. (2009), Ray et al. (2014), and Engel et al. (2017) show a (statistically nonsignificant) positive trend (note that Ray et al., 2014, also shows negative balloon AoA trends in the lower extratropical stratosphere). This discrepancy has been addressed in numerous studies. For example, Garcia et al. (2011) showed that, due to the concave growth rate of tropospheric SF6 concentrations, the AoA trends derived from an SF6 tracer are smaller than the trends derived from a synthetic, linearly growing AoA tracer (also after accounting for the nonlinear growth rates of SF6). They noted that the sparse sampling of in situ observations can be the reason for the abovementioned trend discrepancies. Birner and Bönisch (2011) as well as Bönisch et al. (2011) argued that differences in the changes between the deep and the shallow BDC branch can possibly explain these. Ploeger et al. (2015) showed that the residual circulation transit time cannot explain the AoA trends and that the integrated effect of mixing (which is coupled to residual circulation changes; see Garny et al., 2014) is crucial. Moreover, Stiller et al. (2017) could explain a hemispheric asymmetry by a shift of subtropical transport barriers. In a study based on a chemistry transport model, Kouznetsov et al. (2020) showed that changes in SF6-derived apparent AoA over 1 decade are highly influenced by the SF6 sink and can even turn positive. However, a comprehensive understanding of the magnitude of the individual effects on the AoA trend depending on altitude and latitude is still missing.
SF6 sinks lead to an older apparent AoA (see, e.g., Waugh and Hall, 2002, and Kouznetsov et al., 2020) as well as shorter lifetimes. Leedham Elvidge et al. (2018) evaluated AoA from several tracers including SF6 and found clear differences between them, indicating a shorter SF6 lifetime than previously assumed. The strongest chemical SF6 removal reactions take place in the mesosphere; the most important removal processes are electron attachment and UV photolysis, but these processes have not yet been precisely quantified. Ravishankara et al. (1993) estimated an SF6 lifetime of 3200 years, and Reddmann et al. (2001) found a lifetime of between 400 and 10000 years, depending on the assumed loss reactions and electron density. A more recent model study by Kovács et al. (2017), who used the Whole Atmosphere Community Climate Model (WACCM) to determine the atmospheric lifetime of SF6, reported a mean SF6 lifetime of 1278 years, and Ray et al. (2017) provided a range between 580 and 1400 years based on in situ measurements in the stratospheric polar vortex. The most recent model study of Kouznetsov et al. (2020), who performed simulations of tracer transport with a chemical transport model, shows an SF6 lifetime ranging between 600 and 2900 years. Due to these uncertainties and the complex computation of the chemical reactions, most model studies do not consider any SF6 sinks for the calculation of AoA from SF6 mixing ratios. This can explain why most climate models generally show younger stratospheric air than observations, in particular within the polar vortices (e.g., Haenel et al., 2015; Ray et al., 2017).
In the present study, we apply the chemistry climate model EMAC (ECHAM MESSy Atmospheric Chemistry; Jöckel et al., 2010, 2016) with the aim of understanding the effects of SF6 sinks on tracer-derived AoA and its long-term trends. Specifically, for the first time, we calculate the effect of the sinks on the long-term trend in SF6-derived AoA and quantify how this effect is modulated by circulation changes (recent climate change), specified model dynamics, or by changes in the abundance of relevant species for SF6 chemistry. Furthermore, we analyze the contribution of the SF6 sinks themselves on the long-term trend in SF6-based AoA. As an outlook, we thereupon provide first thoughts on how to apply an AoA correction to observations taking SF6 sinks into account. The chemistry climate model uses the second version of the Modular Earth Submodel System (MESSy2) to link multi-institutional computer codes. In our simulations, we employed the MESSy submodel “SF6” which explicitly calculates SF6 sinks based on physical processes (based on Reddmann et al., 2001), rather than on crude parameterizations. We apply a correction for the nonlinear growth of SF6 in the calculation of AoA, based on Fritsch et al. (2020), which allows for the quantification of the effect of SF6 sinks on SF6-based AoA in isolation. In Sect. 2, we describe the EMAC model and the SF6 submodel as well as the observational data that we use for comparison. Section 3 contains a comparison of the EMAC climatologies with MIPAS data, a comparison of the EMAC trends with MIPAS and balloon-borne measurements, and an analysis of the results of two sensitivity simulations. The model results are discussed in the following using theoretical considerations of the effects of sinks on AoA trends (Sect. 4), including first thoughts on possible correction methods for the sinks, that are highly desirable for the use of observational data. In Sect. 5, we discuss the results and provide some concluding remarks.
2.1 EMAC model
For this study, we use the EMAC (ECHAM MESSy Atmospheric Chemistry, v2.54.0; Jöckel et al., 2010, 2016) model, a numerical chemistry and climate model (CCM) system. It contains the general circulation model (GCM) ECHAM5 (ECMWF Hamburg; Roeckner et al., 2003), with its spectral dynamical core, as well as the MESSy (Modular Earth Submodel System; Jöckel et al., 2005, 2010) submodel coupling interface. The latter is a modular interface structure for the standardized control of process-based modules (submodels) and their interconnections. We apply the model in a T42 horizontal (∼ 2.8∘ × 2.8∘) resolution with 90 layers in the vertical and explicitly resolved middle-atmosphere dynamics (T42L90MA). In this setup, the uppermost model layer is centered at 0.01 hPa, and the vertical resolution in the upper-troposphere–lower-stratosphere (UTLS) region is 500–600 m. In the standard reference setup, we use the basic EMAC modules for dynamics, radiation, clouds, and diagnostics (AEROPT, CLOUD, CLOUDOPT, CVTRANS, E5VDIFF, GWAVE, ORBIT, OROGW, PTRAC, QBO, RAD, SURFACE, TNUDGE, TROPOP, VAXTRA; the reader is referred to Jöckel et al., 2005, 2010, for details on these submodels). Additionally we included the new submodel SF6.
2.2 Submodel SF6
The submodel SF6 is used to calculate the lifetime of SF6 by explicitly accounting for the sinks of SF6 in the mesosphere. The submodel is operationally available for all users in MESSy from version 2.54.0 onward. The calculation method for this is based on the reaction scheme of Reddmann et al. (2001). The most important reaction involved in the chemical degradation of SF6, namely electron attachment, is included in the SF6 submodel. The configuration of the submodel allows for a simple exponential profile for the electron field and a more complex field based on Brasseur and Solomon (1986); in the present study, we use the latter option. It depends on altitude, latitude, solar zenith angle, air density, and day of year. In contrast to Reddmann et al. (2001), UV photolysis of SF6 is not included in the submodel. The loss rate by photolysis is several orders of magnitude weaker than that of electron attachment up to altitudes of about 100 km (see, e.g., Fig. 9 in Totterdill et al., 2015) and is, therefore, not relevant for the focus of our study. Further reactions considered are the photodetachment of (Datskos et al., 1995); the destruction of by atomic hydrogen, hydrogen chloride, and ozone (Huey et al., 1995); the stabilization of excited by collisions; and the autodetachment of . An overview is provided in Table 1. Reddmann et al. (2001) used climatological profiles for the aforementioned gases, whereas channel objects (see Jöckel et al., 2016) are used in our submodel. Such channel objects can be calculated in other submodules (e.g., interactive chemistry), prescribed as external time series (in this study), or just be simple climatologies. The autodetachment rate can be chosen in the namelist and was set to 106 s−1 (see Reddmann et al., 2001). For a general overview of the various reactions, see Fig. S1 in the Supplement.(Jöckel et al., 2016)(Dee et al., 2011)(Jöckel et al., 2016)
2.3 Simulation setup
The simulations performed in this study include a more comprehensive approach for the calculation of the SF6 sinks. We use a climate chemistry model (as opposed to studies based on chemistry transport models as, e.g., in Kouznetsov et al., 2020) and use a more comprehensive SF6 submodel than previous chemistry climate model studies (see, e.g., Marsh et al., 2013, for WACCM). Other than the SF6 submodel, no interactive chemistry is activated in the simulations for this study. The reactant species involved in the SF6 chemistry (HCl, H, N2, O2, O(3P), and O3) and the radiatively active gases (CO2, CH4, N2O, and O3) are transiently prescribed from the ESCiMo RC1-base-07 simulation (see Jöckel et al., 2016) as monthly and zonal means. Moreover, we prescribe the Hadley Centre Sea Ice and Sea Surface Temperature (HadISST) dataset, the CCMI-1 volcanic aerosol dataset (for its effect on infrared radiative heating, see Arfeuille et al., 2013, and Morgenstern et al., 2017), and quasi-biennial oscillation (QBO) nudging (see Jöckel et al., 2016). To compute the photodetachment rate of , we follow Reddmann et al. (2001) and use the extraterrestrial solar photon flux with no attenuation of the UV-photon flux, as provided by WMO (1986). Our simulations range from 1950 to 2011; however, at least the first 10 years have to be considered as a spin-up period. The projection simulation runs from 1950 to 2100 with the SF6 reactant species and greenhouse gases (GHGs) prescribed from the ESCiMo RC2-base-04 simulation (see Jöckel et al., 2016) as monthly and zonal means. In addition to the reference (REF) simulation, we performed two sensitivity simulations and one specified dynamics simulation. The two sensitivity simulations are as follows: the CSS (constant reaction partners for SF6 sinks) sensitivity simulation differs from the REF simulation only with respect to the constant mixing ratios of the reactant species (see above) that influence the SF6 sinks. For that purpose, we kept the mixing ratios at the level of the start of the simulation (year 1950 on repeat). With this simulation, we aim to address the influence of the reactant species involved in the SF6 sink reactions. The second sensitivity simulation, referred to as TS2000, is not a transient simulation, as is the REF simulation, but is instead a time slice simulation with climate conditions (GHGs; sea surface temperatures, SSTs; and SICs) of the year 2000 (climatology of the 1995–2004 period). Furthermore, the reactant species for the SF6 sinks have been averaged over the 1995–2004 period for the TS2000 simulation. This will allow us to investigate the effects of the SF6 sinks under a constant climate. In the specified dynamics (SD) simulation, we apply Newtonian relaxation (“nudging”) towards ERA-Interim (Dee et al., 2011) reanalysis data of potential vorticity, divergence, temperature, and the logarithmic surface pressure up to 1 hPa. This assures that the meteorological situation largely resembles the ERA-Interim data. The flexible structure of MESSy allows us to use the same executable for all simulations, with the differences between them realized through changes in aforementioned namelist settings (see Jöckel et al., 2005). A summary of the simulations used in this study can be found in Table 2.
2.4 Satellite and in situ data
The MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) instrument on Envisat (Environmental Satellite) allowed for the retrieval of SF6 by measuring the thermal emission in the mid-infrared, while orbiting the Earth sun-synchronously 14 times a day. This high-resolution Fourier transform spectrometer measured at the atmospheric limb and provided data for SF6 retrievals in full spectral resolution from 2002 to 2004 and in reduced resolution from 2005 to 2012 between 6 and 40 km of altitude (Stiller et al., 2012; Haenel et al., 2015). In this study, a newer version of the MIPAS dataset, existing as of 2019, will be shown (Stiller, 2021) in which new SF6 absorption cross sections have been used for the SF6 retrieval (Stiller et al., 2020; Harrison, 2020). Except for the newer absorption cross sections for SF6 and the consideration of a trichlorofluoromethane (CFC-11) band in the vicinity of the SF6 signature, the SF6 retrieval and conversion into AoA was done according to the description by Haenel et al. (2015). In particular, the Level-1b data version is still V5.
Engel et al. (2009) collected available air samples of SF6 and CO2 from a balloon-borne cryogenic whole-air sampler flown during 27 balloon flights, with data up to 43 km, and reanalyzed these samples in a self-consistent manner. The derived SF6 data cover the years from 1975 to 2005 (with a gap between 1985 and 1994) and the midlatitudes between 32 and 51∘ N. As AoA profiles from midlatitudes above approximately 25 km or 30 hPa are constant over altitude, a mean midlatitude middle-stratospheric AoA value from each profile was determined by averaging the vertical profile between 30 hPa and the top balloon flight height. With this procedure, a time series of midlatitude middle-stratospheric AoA values could be determined back to 1975. It is important to note that part of the Engel et al. (2009) AoA time series is derived from CO2 measurements. In this paper, only the AoA data points derived from SF6 are used. Engel et al. (2017) extended the initial dataset from Engel et al. (2009) to 2016, but the AoA here is derived from CO2 measurements and is, thus, also excluded.
2.5 Analysis method
The basic concepts for the calculation of mean AoA are introduced in Hall and Plumb (1994). In the case of a tracer with a linear increasing lower boundary condition, AoA can be determined by the time lag between the mixing ratio at a given point in the atmosphere and the same mixing ratio of the reference time series. As for any realistic tracer, SF6 does not exhibit perfectly linear growth, for which adjustments in the AoA calculation are needed. This study follows the calculation method employed by studies such as Engel et al. (2009), which was introduced in Volk et al. (1997). The calculation uses a polynomial fit to the reference time series to approximate mean AoA. However, in our study, we modified the parameters compared with those used in Engel et al. (2009) to ensure that (passive) SF6-based AoA agrees with the ideal AoA derived from the linear tracer, following Fritsch et al. (2020). Specifically, in our calculations, we used a ratio of moments of 1.0 years and a fraction of input of 95 %. Further, we used the SF6 mixing ratio averaged over 20∘ S–20∘ N at the ground as the reference time series. Due to the availability of data, this reference region is also used in Engel et al. (2009) with balloon-borne observations.
For the derivation of AoA from MIPAS SF6 observations, Stiller et al. (2008, 2012) and Haenel et al. (2015) used a slightly smoothed version of the global mean of SF6 surface measurements as the reference time series instead of SF6 at the stratospheric entry point, which is not available from observations (see, e.g., Dlugokencky, 2005). The nonlinearity of the reference curve was considered by its convolution with an idealized age spectrum parameterized as a function of the mean age within an iterative approach. For more details, see Stiller et al. (2012) and Haenel et al. (2015).
In our simulations, the AoA calculations are applied to a total of four tracers, which can be organized into two groups. The first assumes a strict linear growth of SF6, producing a linear reference curve, whereas the second considers a realistic growth of SF6 based on observed emissions, creating a nonlinear SF6 reference curve. Technically, in our simulations these “emissions” are realized via lower boundary conditions, which are based on surface observations, as in Jöckel et al. (2016). As previously mentioned, SF6 undergoes chemical degradation predominantly in the mesosphere. Consequently, the absence or presence of mesospheric sinks is additionally considered, resulting in a total of four tracers: tr(WS, SF6), tr(NS, SF6), tr(WS, lin), and tr(NS, lin). The labeling of these depends on the chemistry involved (with sinks: “WS”; without (no) sinks: “NS”) and the growth assumed (linear: “lin”; nonlinear: “SF6”) and follows the pattern tr(Chemistry, Growth). When referring to simulations with a specific tracer, the labeling will follow the notation Simulation(Chemistry, Growth), and similarly we use the following notation for AoA inferred from the tracer used in a simulation: AoA(Chemistry, Growth)SIM. For example, AoA from the reference simulation based on the tracer with mesospheric sinks and nonlinear (SF6-emission-based) increase is referred to as “AoA(WS, SF6)REF”, and AoA inferred from the linear tracer without mesospheric sinks in the reference simulation is denoted as “AoA(NS, lin)REF”.
3.1 SF6 mixing ratios
In order to evaluate the SF6 mixing ratios simulated by the EMAC model, we first analyze the four tracers in the SD simulation in comparison to observational data. This study does not perform a detailed comparison of SF6 profiles, as the major aim is not an in-depth evaluation of the SF6 submodel but rather a quantification of the potential effects of the SF6 sinks on AoA and its long-term trends. However, to ensure that SF6 values in the EMAC model are within the range of observational estimates, we perform selected comparisons to data from Ray et al. (2017) and MIPAS SF6 (Stiller et al., 2020, paper in preparation). Figure 1a depicts the modeled SF6 vertical profile climatologies in comparison with MIPAS SF6. Figure 1b shows the modeled SF6 vertical profiles and balloon-borne measurements of SF6 (Ray et al., 2017) on a particular day. The former comparison is for zonal mean SF6 averaged over 30–50∘ N and 2007–2010. These years have been chosen as the dataset is complete in this period. The error bars represent the standard deviation of the zonal mean ensemble, which consists of the measurement noise error of MIPAS, further random error sources from the retrieval (e.g., temperature uncertainties), the natural variability over the longitudes of the latitude band, and the 4 years of averaging (2007–2010).
Figure 1b shows modeled SF6 mixing ratios of the day of the balloon flight. The balloon was launched on 5 March 2000 at 67∘ N in Kiruna, Sweden. To ensure that the SF6 profile is based on air masses from within the vortex, modeled SF6 values are averaged over 65–80∘ N and 0–100∘ E, which corresponds to the area of the vortex core for the given day. The standard deviation of SF6 for 0–100∘ E averaged over the respective latitude range is shown as error bars.
The tracers with nonlinear growth in the SD simulation show smaller tropospheric SF6 mixing ratios than the linear tracers (Fig. 1a, b). This can be explained by the two different growth scenarios of SF6 and the prescribed lower boundary conditions (see Fig. S2 in the Supplement). The sinks do not have a considerable effect in the troposphere; hence, the effect of the SF6 sinks only becomes noticeable higher up. Furthermore, the effect of the SF6 depletion becomes increasingly evident with altitude. This is portrayed in the growing differences with altitude between tr(WS, lin) and tr(NS, lin) and the nonlinear equivalent. The differences particularly increase for the tracers with linear emissions, as these exhibit higher SF6 mixing ratios and, hence, experience greater SF6 depletion than those with nonlinear boundary conditions. Due to the small turnaround times for air in the middle atmosphere, the tracers without sinks exhibit a very low decrease in the SF6 mixing ratios with altitude.
Figure 1a shows that the EMAC-simulated nonlinear SF6 is within the observed range of MIPAS SF6. Below 30 km, MIPAS SF6 mixing ratios are smaller, with a near-constant offset of approximately 0.5 pmol mol−1 up to 20 km. Above 30 km, MIPAS SF6 shows larger mixing ratios than EMAC. This means that EMAC SF6 (SD(WS, SF6)) shows a larger decrease with altitude than MIPAS SF6, suggesting that the sinks in EMAC are too strong. Another explanation could be overly strong vertical mixing in EMAC. However, the EMAC SF6 lies within the MIPAS uncertainty range throughout the atmosphere. The standard deviation increase with height in MIPAS SF6 can be attributed to the decrease in the SF6 signal with height, which leads to an increase in the noise error of SF6. Additionally, the natural variability in SF6 itself, as well as the evolution of SF6 over time, contribute to the increasing standard deviation in the MIPAS SF6 profile. The increase in the standard deviation with height can also be seen in the EMAC SF6 profiles, particularly in the tracers tr(WS, SF6) and tr(WS, lin). However, it is (by far) not as large as in the MIPAS data because the simulations have no measurement error and possibly show a smaller natural variability than the observations.
The balloon flight SF6 profile (Ray et al., 2014) in Fig. 1b largely resembles the profile of the realistically modeled tracer tr(WS, SF6). Below 25 km, the modeled SF6 profile shows a constant high bias of around 0.3 pmol mol−1, presumably due to the lower boundary conditions used. Larger discrepancies can be seen between 25 and 35 km altitude, with higher mixing ratios of the modeled SF6. As the data presented in Fig. 1b are only for a specific day and region, the particular meteorological situation and location of the balloon can be crucial for the comparison.
3.2 SF6 lifetimes
It is based on the reaction rates (ki) of the chemical reactions (Ri) marked in Table 1, the branching fraction ϵ (taken as 0.999; see Reddmann et al., 2001), and the efficiency of the SF6-recovery reactions (η), where η is calculated as
Only the realistic tracer tr(WS, SF6) is considered. The reference simulation yields an average lifetime of 2101 years. Reddmann et al. (2001) carried out sensitivity experiments with this scheme, using various options for the chemical mechanisms. In this way, the lifetime could be varied between 400 and 10 000 years. Our value is below the value of 3200 years calculated by Ravishankara et al. (1993) but above the values of 1278 and 850 years of the more recent studies by Kovács et al. (2017) and Ray et al. (2017), respectively. In another new modeling study, Kouznetsov et al. (2020) presented a range of 600–2900 years. Therefore, our value of around 2100 years is within, although rather at the upper range, of the uncertainties. In contrast to the comparison of the model results with SF6 observations shown in the previous section, our lifetime value points towards rather weak SF6 sinks in our scheme. To assess the variability in the atmospheric SF6 lifetime, we show the time series of the SF6 lifetimes for the four simulations that were described in Sect. 2 (Fig. 2). The lifetime of the REF simulation lies at about 2500 years in 1965 and decreases by approximately 25 % to 1900 years in 2011. The lifetime of the projection simulation behaves similarly and increases to about 2400 years by the year 2100. This shape of the lifetime resembles that of projected O3 (see, e.g., Eyring et al., 2007), which is reflected in the ESCiMo simulations (Jöckel et al., 2016) from which the O3 and the other SF6 reactant species are prescribed here. However, our SF6 degradation scheme includes a number of simplifications, which can modify the lifetimes. For example, a constant profile is used to prescribe the sinks through NO. This can simplify the long-term variability in the SF6 lifetimes by making it overly dependent on the species that are transiently prescribed in the simulations. Apart from seasonal variations, the lifetimes of the CSS and the TS2000 simulations are fairly constant, their lifetime values lie around 2600 and 1900 years, respectively. The fact that the lifetimes of the CSS simulations are fairly constant implies that the long-term trends of the SF6 lifetimes can mostly be attributed to the abundance of the species involved in the SF6 degradation. Variations in stratospheric temperatures or the circulation strength seem to only play a minor role.
3.3 Age of air climatologies
AoA climatologies averaged over the period from 2007 to 2010 are shown in Fig. 3: from the reference simulation for REF(NS, lin), REF(NS, SF6), REF(WS, lin), and REF(WS, SF6) in Fig. 3a–d, respectively; from the specified dynamics simulation (SD(WS, SF6)) in Fig. 3e; and from MIPAS SF6 observations in Fig. 3f (Stiller, 2021). The years 2007 to 2010 were chosen as these are the only years with complete MIPAS data.
In all cases, AoA increases with increasing altitude and latitude. The cases considering sinks show an older apparent AoA than those without, especially with increasing altitude and latitude. This apparent aging can be explained by the fact that the inclusion of mesospheric sinks results in smaller SF6 mixing ratios. The reduced mixing ratios lead to a seemingly older AoA as the corresponding reference value lies further in the past. Downwelling within the polar vortices transports old air from the mesosphere to the stratosphere. With the breakdown of the polar vortex at the end of the winter season, the old air is then mixed into lower latitudes. The relative effect of the sinks on AoA derived from the SF6 tracers with nonlinear growth can be seen in Fig. 3g and h. Mean AoA values derived from SF6 in the early period (1970–1980) are moderately affected by the sinks, with a difference of around 20 %–25 % in the polar middle stratosphere and above (i.e., for mean AoA values above 5.5 years). Differences are small (less than 10 %) for mean AoA values below 4 years. However, as will be discussed (see Sect. 3.5), the effect of the SF6 sinks increases over time, and for the later period (2000–2010), mean AoA derived from SF6 is considerably affected by the sinks. Differences greater than 20 % can be seen in Fig. 3h) for AoA above about 3 years, and in the extratropical lower stratosphere, differences are larger than 10 % even for mean AoA values of 2 years and above.
The patterns of modeled AoA from the linear tracers are similar to those of the SF6-emission-based tracers (compare Fig. 3a with Fig. 3b and Fig. 3c with Fig. 3d). However, the similarities are weaker in the polar regions, especially when the SF6 sinks are considered (Fig. 3c, d). In these regions, REF(WS, lin) reaches AoA values spanning 10–15 years or more, and REF(WS, SF6) reaches AoA values of 9–14 years. This difference can be attributed to the greater initial growth of the tracer with linear emissions than that with nonlinear emissions (Fig. S2). This leads to enhanced SF6 mixing ratios in the REF(WS, lin) case and, in turn, strengthens the influence of SF6 sinks on AoA. This is particularly relevant in the winter months when SF6-depleted mesospheric air is transported downward into the polar stratosphere.
AoA derived from SF6 emissions including chemical SF6 sinks (REF(WS, SF6); Fig. 3d) agrees best with MIPAS AoA (Fig. 3f). Overall good agreement between EMAC and MIPAS AoA is found in the tropics, but there is a large high bias in the high latitudes in EMAC: between 40 and 50∘ N, modeled AoA is up to 2 years older than MIPAS AoA in the stratosphere and up to 3 years older in the polar upper stratosphere (see Fig. S5 in the Supplement). In comparison to the MIPAS observations, the SF6 sinks therefore seem to be too strong in the model, as already mentioned above. However, in comparison with the previously published MIPAS data (Stiller et al., 2012; Haenel et al., 2015), the EMAC AoA was actually too young (i.e., the MIPAS AoA was much older in the previous versions in the polar regions). The spectroscopic data used for the SF6 retrieval in MIPAS cause a rather large bias (that has now been corrected by improved spectroscopy). The new spectroscopic data lead to a considerably younger AoA in the middle to upper stratosphere. There are, however, good reasons to believe that the most recent MIPAS data are improved compared with the previous ones: the spectroscopic data used are far better characterized than the previous ones (Harrison, 2020), and the new AoA data from MIPAS agree significantly better with independent measurements than the previous version, in particular at higher altitudes (Stiller et al., 2020; Stiller, 2021). On the other hand, free-running EMAC simulations generally have an overly weak Antarctic polar vortex (see Jöckel et al., 2016), which is, however, stronger than that of the reference simulation. Therefore, the more stable vortex in the SD simulation leads to enhanced isolation and aging of polar stratospheric air, especially in the Southern Hemisphere during austral spring (see Fig. S4 in the Supplement for details). This could somewhat resolve the issue for the comparison with the previous MIPAS AoA version (see Stiller et al., 2012, and Haenel et al., 2015, for further details); for the present MIPAS version, however, the discrepancies in the high altitudes and latitudes with the EMAC SD simulation are even larger than with the REF simulation. The comparison of Fig. 3e with Fig. 3f illustrates that the model cannot reproduce the tropical pipe and exhibits too much horizontal mixing, or overly slow upwelling. This could explain why the model does not reproduce the constant AoA with height in the midlatitudes. A detailed assessment of both the satellite data and the model simulations is necessary to resolve these discrepancies. In the model, dynamical effects like the strength of the polar vortex or the gravity wave parameterization can play important roles in the downwelling strength. Moreover, various processes of chemical SF6 removal can be revised and/or parameterized differently. Here, we showed that SF6 sinks have the potential to resolve the differences between simulated and observed climatologies of AoA and that EMAC AoA lies within the uncertainties of MIPAS AoA throughout the atmosphere. Therefore, we consider our simulations suitable for studying the temporal evolution of AoA.
3.4 Apparent age of air trends
In this section, we analyze the EMAC AoA trends and compare them with observation-based AoA trends. Figure 4 shows the AoA time series and linear regressions from the REF and the SD simulations as well as from MIPAS observations (Stiller, 2021) and from the SF6 measurements by Engel et al. (2009). As the latter were collected from balloon flights in the Northern Hemisphere midlatitudes at around 30 km altitude, the EMAC and MIPAS data are also taken from that height and averaged over 30–50∘ N for consistency.Engel et al. (2009)
a Trend calculated following methods of Haenel et al. (2015). b Trend calculated over 1980–2011 at 30 km altitude. c MIPAS trend calculated over 2002–2012. d Trend calculated over 1975–2005 between 32 and 51∘ N and between 24 and 35 km.
For better quantification of the trends, Table 3 provides the AoA trend values of the EMAC simulations for two periods. The trend in the entire simulation period from 1965 to 2011 is taken for long-term trend assessment and comparison with the measurements by Engel et al. (2009). For the comparison with MIPAS data, the EMAC AoA trends are shown for the 2002–2011 period between 30 and 50∘ N, for the realistic tracer tr(WS, SF6). The trend calculation follows that of Haenel et al. (2015).
The tracers without SF6 sinks lead to negative AoA trends, which are consistent with the simulated acceleration of the BDC in the course of climate change (e.g., Garcia and Randel, 2008). Positive AoA trends are obtained for all tracers that take SF6 chemistry into account. The trend of 0.19 ± 0.01 yr per decade in the REF simulation (WS, SF6) is within the limits of the uncertainties of the trend obtained by Engel et al. (2009), who calculated an AoA trend of 0.24 ± 0.22 yr per decade. This means that, in our simulations, the sinks help to reconcile the modeled and the measured AoA trends over the recent decades. Note that Engel et al. (2009, 2017) also obtained a positive trend for AoA derived from CO2 measurements.
Haenel et al. (2015) calculated MIPAS AoA trends for the period from 2002 to 2012 of 0.25 ± 0.11 yr per decade for 30–40∘ N at 30 km and of 0.24 ± 0.11 yr per decade for 40–50∘ N at 30 km. For the new MIPAS dataset, the AoA trend is 0.34 ± 0.13 yr per decade. Note that these trends are calculated by applying a bias correction for the discontinuity between the two different observational periods of MIPAS (see Fig. 4); a description of the method can be found in von Clarmann et al. (2010).
Following the trend calculation used in Haenel et al. (2015), the REF(WS, SF6) and SD(WS, SF6) time series show that EMAC AoA bears a good resemblance to that of the new MIPAS retrieval, with an AoA trend of 0.22 ± 0.12 and 0.50 ± 0.13 yr per decade, respectively. Consistent with the trend calculation of the MIPAS data, the variability due to the QBO is considered by a respective term in the multivariate linear analysis. However, this measure induces only small differences in the EMAC trend calculations (see Table S1 in the Supplement for further details). Note also that the rather short period of the MIPAS observations implies rather large uncertainties in the trends because interannual variability can have a large effect on the trend calculations. This is also apparent from the highly variable trend signals in the different simulations; for example, the TS2000(NS, lin) exhibits a strongly negative trend over this period despite no forced long-term trends. This is confirmed by a trend of −0.00 ± 0.01 yr per decade over the 1965–2011 period.
As shown in Fig. 5, the strong deviation of SF6-derived AoA trends holds almost throughout the stratosphere. Only below about 20 km altitude for the period from 1965 to 2011 are the effects of the SF6 sinks smaller than the trend uncertainty. For the trend in the shorter time period of 2000–2011, effects from the SF6 sinks become significant at about 22 km altitude and higher, which is mainly due to the larger uncertainty stemming from a smaller time period. This means that the SF6-based AoA trends with and without sinks are not distinguishable from each other up to 20 and 22 km altitude in the midlatitudes, depending on the period. Furthermore, the uncertainty in the trend calculated from SF6-based AoA with sinks increases with increasing altitude – this is to be expected, as the effect of the SF6 sinks increases with increasing altitude. The trend in AoA(NS,SF6)REF over the 2000–2011 period is positive at 25 km and negative over the longer 1965–2011 period. However, it is important to note that the time period of about 1 decade implies that the trend is strongly influenced by interannual variability (see, e.g., Dietmüller et al., 2021). The strong influence of interannual variability also explains the difference in trend values in Fig. 5, for which trends are calculated with a simple linear fit, versus the value for a similar period in Table 3. The latter trends were calculated using a regression model taking other variability modes into account to enable a comparison with MIPAS.
3.5 Explanations of apparent age of air trends
In this section, we will analyze the EMAC apparent AoA trends, in particular the sign change of the trend when SF6 sinks are switched on. For this, Fig. 6 shows the AoA (WS,SF6) time series of the sensitivity simulations CSS and TS2000 averaged between 30 and 50∘ N. For comparison, AoA values from the reference simulation (REF(WS, SF6) and REF(NS, lin)) are also included.
In the CSS simulation, the mixing ratios of the reactant species involved in the sink reactions of SF6 are held constant. The CSS simulation with the realistic tracer (CSS(WS, SF6)) shows an AoA trend of 0.11 ± 0.01 yr per decade over the 1965–2011 period. This is lower than the AoA trend in the reference simulation REF(WS, SF6) and means that while the changes in the SF6-depletive substances influence the magnitude of the positive trend, they cannot explain the positive sign of the AoA trend.
The TS2000 simulation is a time slice simulation with climate conditions of the year 2000. AoA of that simulation derived from the realistic tracer (TS2000(WS, SF6)) shows a positive trend of 0.23 ± 0.02 yr per decade. This is even stronger than the trend in the REF simulation (0.19 ± 0.01 yr per decade). By definition, the TS2000 simulation does not feature any changes in its climatic state nor in the composition of the atmosphere. This is reflected by the fact that the idealized tracer tr(NS, lin) does not show a trend in this simulation (see Table 3; TS2000(NS, lin): −0.00 ± 0.01). The temporal increase in apparent AoA rise in the TS2000(WS, SF6) simulation despite no climate changes therefore points to the fact that the SF6 sinks themselves lead to the positive trend. The difference between the TS2000 and the REF AoA trends for both tr(WS, SF6) and tr(NS, lin) reflect the negative contribution of the accelerating BDC to the trend, which can be seen in the idealized AoA trend from the REF simulation (NS, lin). Overall, neither changes in the SF6 sinks nor changes in the stratospheric circulation due to climate change are responsible for the positive trend found in AoA with sinks. Instead, the results indicate that the sinks themselves can generate a positive trend.
For a complete discussion of the features that can be seen in Fig. 6, we now also describe the sudden decrease in AoA(WS, SF6) shortly after 1982 seen in the REF simulation. It is also visible in the two sensitivity simulations (CSS and TS2000). This means that changes in the SF6-depleting substances as well as climate change and volcanic activity (in particular the volcanic eruption of El Chichón in 1982) can be ruled out as possible causes of the drop. The solar cycle is not taken into account in the time slice simulation. The time series of the realistic SF6 tracer emissions show a stronger increase in the 1980s (not shown). Fritsch et al. (2020) showed that, due to this increase, the calculation of AoA based on SF6-like tracers is more sensitive to the chosen parameters in the AoA derivation for this time. The drop in AoA that we see is caused by this limitation of the derivation from SF6-like tracers.
The increasing variability and trend in AoA in the last 2 decades, seen in the simulations with mesospheric SF6 chemistry, can be attributed at first order to the SF6 depletion reactions (see Fig. 6). In particular, the effect of mesospheric SF6 sinks is stronger with higher SF6 mixing ratios. Subsequent downward transport of SF6-depleted air into the vortex and in-mixing thereof into lower latitudes after the vortex breakup results in apparent older AoA. This can explain that the annual variability increases over time due to the increase in SF6 mixing ratios.
The following section focuses on the theoretical examination of the link between SF6 sinks and positive AoA trends. First, we will show, from a theoretical standpoint, that SF6 sinks with constant destruction rates lead to a positive trend in AoA. Based on the theoretical considerations and on the model data, we will then discuss the possibilities of a correction of AoA derived from SF6 data for the effects of the sinks. Many observational mean AoA estimates are based on measurements of SF6, and given that the relevance of the sinks of the mean AoA estimates increases over time, a correction method is required to obtain unbiased information on stratospheric transport strengths.
The aforementioned link between the positive AoA trend and mesospheric SF6-depletive chemistry, based on Hall and Waugh (1998), is illustrated below and follows the mathematical formulations put forward by Hall and Plumb (1994) and Schoeberl et al. (2000). To allow for analytical expressions, we will only consider the case of a linearly increasing tracer here, and we will further assume that the lifetime of SF6 is constant in time. We consider a tracer χ(t) experiencing relative loss with time-constant lifetime τ. Note that τ is equivalent to the inverse of the loss rate. We denote the reference mixing ratio as χo(t) and assume a constant growth rate of Δχoo:
For any location, we can then express the tracer mixing ratio as
The expression under the first integral corresponds to the arrival time distribution . This represents the transit time distribution of a chemically depleted tracer with lifetime τ (see, e.g., Plumb et al., 1999; however, note that G∗ is not normalized in our definition, following Engel et al., 2017). Note that τ is a measure of the lifetime an air parcel experiences on its path of length t′ and is referred to the path-integrated lifetime (not to be confused with the local lifetime). The second integral in Eq. (5) is the first moment of the arrival time distribution (i.e., the mean arrival time and is denoted as Γ∗). Using these terms, Eq. (5) can be expressed as
For the tracer without sinks (here referred to as the passive tracer χp(t)), the integral over equals 1, and Γ represents the mean AoA. Equation 6 can then be rearranged to
giving the common expression to derive mean AoA from a linear increasing tracer.
The change of apparent mean AoA with time can then be expressed as
In the case of a passive tracer (i.e., a tracer without sinks), is equal to the age spectrum G(t′); thus, its integral equals 1. Consequently, the AoA trend is zero in the absence of a trend in transport (constant G). However, for the depleted tracer (with sinks), a positive linear trend in the apparent mean AoA () is induced by the sinks, even if the circulation strength (i.e., G(t′)) and the lifetimes τ(t′) are constant. This is consistent with the results shown in the previous sections. In particular, we showed a positive trend in the linearly increasing SF6-like tracer in the time slice experiment (TS2000, see, e.g., Table 2), which satisfies all of the conditions made above. In the reference simulation, the BDC is accelerating over time, so that mean AoA from the passive tracer is decreasing. However, in our simulations, the positive trend induced by the sinks overcompensates for the impact of the BDC acceleration. It is important to note that the above is valid for a linearly increasing tracer and that variable growth rates will modify the influence of SF6 on apparent AoA.
A correction of the mean AoA estimate for the effects of the sinks would require knowledge of the (integrated) arrival time distribution G∗ (or, equivalently, the transit-time-dependent path-integrated lifetimes τ(t′)), a quantity that is not readily available and possibly nonlinearly dependent on Γ. As a thought experiment, with the aim of deriving an analytical concept for the correction of mean AoA for the sinks, we make the hypothetical assumption that the age spectrum is represented by a single, average path (i.e., we assume the age spectrum is a delta function ). With this assumption, the tracer mixing ratio of the depleted tracer (Eq. 5) simplifies to
where τeff is the path-integrated lifetime along this single path. To avoid confusion with either the local lifetime or the averaged transit-time-dependent path-integrated lifetime, we will refer to τeff as the “effective” lifetime.
Given that the mean age Γ is about 2 orders of magnitude lower than the lifetime τeff, we can approximate ; thus, . Again for small Γ compared with τeff, the second term can be neglected so that
In general, τeff will be dependent on Γ. However, if the effective lifetime were constant, the apparent mean age would be linearly related to the actual mean AoA. The slope between the latter two then is merely a function of the effective lifetime and increases linearly over time t. Thus, for constant circulation (Γ) and effective lifetime τeff, a linear trend in strength
Under these assumptions, a relatively simple, linear relation between the true mean AoA (Γ) and the apparent AoA () is obtained. While the assumptions will clearly be violated in the model, we investigate, based on data from our model simulations, whether the violations might effectively be small enough so that the linear relation still holds. In particular, the linearly increasing tracers from the time slice simulation (TS2000) with constant circulation strength satisfy the initial assumptions (linear increase and constant circulation), and we can test whether the approximations of a single average path hold.
Figure 7 shows a quasi-linear relation between the ideal AoA and AoA derived from the linearly increasing SF6-like tracer for mean AoA values below approximately 4 years. The slope increases over time (as apparent from the transition from yellow to red colors in Fig. 7), consistent with the simplified theoretical considerations for a constant effective lifetime τeff (see Eq. 12). However, for mean AoA values above about 5 years, this quasi-linear regime does not hold anymore. We instead find exponential growth of SF6-based apparent mean AoA. Those older mean AoA values are closer to the region of depletion in the mesosphere, and a constant effective lifetime (as assumed for the linear regime) does not hold anymore. Rather, we can expect τeff to be strongly dependent on Γ.
Based on Eq. (11) (or its linear approximation Eq. 12), we can estimate the effective lifetime from the model data of Γ (from the passive tracer) and (from the linearly increasing SF6-like tracer). The resulting lifetime is shown in Fig. 8 as a function of the ideal mean AoA. Consistent with the previous results, τeff varies comparatively little with Γ for a range between 1 and 4 years. Figure 7 also includes the theoretical values of according to Eq. (11) for t = 40 years and when using either a constant mean value for τeff of 140 years (obtained for Γ = 3.3 years) or the Γ-dependent values. The approximation of a constant effective lifetime matches the model data for ages below about 4 years reasonably well, but for older air, the assumption of a constant τeff obviously fails. When using the Γ-dependent value for τeff, on the other hand, the exponential increase in is generally well captured. This justifies the applicability of the simplified assumptions of one average path (Eqs. 11, 12).
Overall, the results derived in this section indicate that a correction of observational SF6-derived mean AoA for the effects of chemical sinks is likely possible by applying a time-dependent linear correction function. This linear relation between AoA from the ideal and from the chemically depleted SF6 tracer holds for mean ages below about 4 years. Here, we show this relation for the linearly increasing tracer in northern midlatitudes. However, further analysis indicates that this linear relation also holds for the realistic SF6 tracer with nonlinear growth over time (not shown). Furthermore, the linear relation seems to be nearly identical for different latitude bands (not shown), which is a very promising property for future applications of a correction method.
We emphasize that the strength of the SF6 sink in our model simulations is not well enough constrained to properly establish such a correction function. Deriving suitable values for the linear relation between Γ and (and, thus, the effective lifetime τeff) should be obtained by means of observational data. This could be achieved by using simultaneous measurements of SF6 and other age tracers, as previously shown by Leedham Elvidge et al. (2018) and Adcock et al. (2021). Furthermore, the concept needs to be evaluated vigorously on model data to assess its errors and limitations.
For over a decade, disagreements with regards to stratospheric AoA and, particularly, its trends between model simulations and observations have raised many questions from scientists. AoA from observations is mostly older than AoA from model simulations, and models simulate a decrease in AoA over recent decades, whereas trend estimates from observational data report a nonsignificant positive trend. In agreement with our results, previous studies (see, e.g., Kouznetsov et al., 2020) have shown that the chemical sinks strongly influence SF6-derived AoA in terms of absolute values and decadal changes. We investigate, for the first time, how longer-term trends are affected in a consistent manner, and we explore the different contributions from circulation changes, changes in abundance of reaction partners, and trends induced by constant destruction rates. Thus, to make this step towards understanding the reasons for the discrepancies, we study the impact of mesospheric SF6 sinks on AoA climatologies and trends using the chemistry climate model EMAC (Jöckel et al., 2010, 2016) with the SF6 submodel (Reddmann et al., 2001). This submodel allows for explicit calculation of SF6 sinks, and we applied a correction for the nonlinear growth of SF6 in the calculation of AoA (Fritsch et al., 2020).
The EMAC SF6 mixing ratio profiles show good agreement with balloon-borne measurements as well as with satellite observations. Some of the differences between the model and observations are within the uncertainty range of the observations. However, reasons for the quantitative differences in the high latitudes and altitudes can also be found in deficiencies in the representation of the SF6 sinks in the model or in the dynamics simulated by the model. The EMAC reference simulation yields a global stratospheric SF6 lifetime of 2100 years, varying between 2500 and 1900 years for the simulation period. This value lies within the range of 600–2900 years provided by the model study of Kouznetsov et al. (2020) and below the value 3200 years calculated by Ravishankara et al. (1993). Kovács et al. (2017) and Ray et al. (2017) recently found somewhat lower SF6 lifetimes of 1278 and 850 years, respectively. Although this shows that large uncertainties still exist in determining the SF6 lifetime, these results also confirm that the EMAC SF6 depletion mechanisms are reasonable. In our transient simulations (REF and PRO), the SF6 lifetimes vary by about 25 % following the abundances of the reactant species, basically resembling the pattern of the stratospheric ozone concentration. This behavior, however, may be a result of the fact that several effects that potentially influence SF6 lifetime variability are not implemented in full detail.
The inclusion of SF6 sinks translates into apparent older stratospheric air; thus, when SF6 sinks are enabled, EMAC AoA compares better with MIPAS satellite observations. In the tropics, good agreement can generally be found, but the results also indicate that EMAC does simulate an overly broad, less-isolated tropical pipe. In polar regions, however, EMAC AoA is higher than MIPAS AoA. In comparison with the previously published MIPAS data (Stiller et al., 2012; Haenel et al., 2015), the EMAC AoA in the polar regions was actually too low. More research on both models and observations is necessary to resolve these remaining discrepancies.
In this study, we show the effect of SF6 sinks on the tracer-derived AoA trend for both longer time series as well as for the short time series corresponding to the MIPAS time frame, although issues in variability are associated with the latter. Without SF6 sinks, EMAC shows a negative AoA trend over 1965–2011. This is consistent with the simulated acceleration of the Brewer–Dobson circulation resulting from climate change (see, e.g., Garcia and Randel, 2008, Butchart et al., 2011, and Eichinger et al., 2019). The inclusion of chemical SF6 sinks leads to positive AoA trends throughout the stratosphere (except the tropical lower stratosphere below 50 hPa) in our simulations, which, in turn, is consistent with the positive AoA trend derived from MIPAS in the Northern Hemisphere (Haenel et al., 2015). Moreover, the SF6 sinks help to improve the agreement of our model results with the AoA derived from the balloon-borne in situ measurements by Engel et al. (2009), from which a (nonsignificant) positive AoA trend was obtained. However, this only accounts for SF6-derived AoA, and our results cannot help explain the positive trend in the CO2-derived AoA in Engel et al. (2009, 2017). Furthermore, it has recently been shown that the AoA trend derived in Engel et al. (2009, 2017) was likely overestimated due to nonideal parameter choices in the calculation of AoA (Fritsch et al., 2020).
Our sensitivity studies quantified that the positive AoA trends are neither a result of climate change nor of changes in the substances involved in SF6 depletion. The SF6 sinks themselves are the reason for the increase in apparent AoA. The reason for that is the temporally increasing influence of the chemical SF6 sinks on AoA. In our simulations, this effect overcompensates for the effect of the simulated acceleration of the stratospheric circulation, leading to a net increase in AoA. Due to various sources of uncertainties, this result bears quantitative leeway and has to be assessed in finer detail. Nevertheless, for now, we can conclude that SF6 sinks have the potential to explain the long-lasting AoA trend discrepancies between models and observations. Furthermore, we put forward a first approach towards a method for SF6 loss correction, which shall be further developed and applied on observational data in the future. From our first analyses, we can conclude that a linear correction (that is dependent on both time and the effective lifetime of SF6) can likely be applied to AoA values up to 4 years. However, further studies with more comprehensive approaches are required for a precise quantification of these values.
The Modular Earth Submodel System (MESSy) is continuously developed and applied by a consortium of institutions. Use of MESSy and access to the source code are licensed to all affiliates of institutions that are members of the MESSy Consortium. Institutions can become a member of the MESSy Consortium by signing the MESSy Memorandum of Understanding. More information can be found on the MESSy Consortium website (http://www.messy-interface.org, last access: 27 March 2020). The exact code version used to produce the simulation results is archived at the German Climate Computing Center (DKRZ) and can be made available to members of the MESSy community upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/acp-22-1175-2022-supplement.
SL, RE, and HG designed the study. SL analyzed the data with support from RE, HG, and FF. RE performed the simulations. TR and SV implemented the SF6 submodel in MESSy. GS and FH provided the MIPAS data analyses. SL and RE drafted the manuscript, and all authors helped with discussions and with finalizing the paper.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “The Modular Earth Submodel System (MESSy) (ACP/GMD inter-journal SI)”. It is not associated with a conference.
This work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under project ID bd1022. Moreover, we thank Eric Ray and Andreas Engel for providing observational data, Axel Lauer for his helpful suggestions, and the reviewers (see the Review statement) for their thorough reviews of the paper.
This research has been supported by the Helmholtz-Gemeinschaft (grant no. VH-NG-1014) and the Grantová Agentura Ceské Republiky (grant nos. 21-20293J and 21-03295S).
The article processing charges for this open-access publication were covered by the German Aerospace Center (DLR).
This paper was edited by Jerome Brioude and reviewed by Eric Ray and Rostislav Kouznetsov.
Adcock, K. E., Fraser, P. J., Hall, B. D., Langenfelds, R. L., Lee, G., Montzka, S. A., Oram, D. E., Röckmann, T., Stroh, F., Sturges, W. T., Vogel, B., and Laube, J. C.: Aircraft-Based Observations of Ozone-Depleting 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
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 CO2, CH4, and N2O, J. Geophys. Res.-Atmos., 106, 32295–32314, https://doi.org/10.1029/2001JD000465, 2001. a
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/acp-13-11221-2013, 2013. a
Birner, T. and Bönisch, H.: Residual circulation trajectories and transit times into the extratropical lowermost stratosphere, Atmos. Chem. Phys., 11, 817–827, https://doi.org/10.5194/acp-11-817-2011, 2011. a
Bönisch, H., Engel, A., Birner, Th., Hoor, P., Tarasick, D. W., and Ray, E. A.: On the structural changes in the Brewer-Dobson circulation after 2000, Atmos. Chem. Phys., 11, 3937–3948, https://doi.org/10.5194/acp-11-3937-2011, 2011. a
Brasseur, G. and Solomon, S.: Aeronomy of the middle atmosphere, 2 edn., Reidel, Dordrecht, The Netherlands, 1986. a
Brewer, A. W.: Evidence for a world circulation provided by the measurements of helium and water vapour distribution in the stratosphere, Q. J. Roy. Meteor. Soc., 75, 351–363, https://doi.org/10.1002/qj.49707532603, 1949. a
Butchart, N. and Scaife, A. A.: Removal of chlorofluorocarbons by increased mass exchange between the stratosphere and troposphere in a changing climate, Nature, 410, 799–802, https://doi.org/10.1038/35071047, 2001. a
Butchart, N., Charlton-Perez, A. J., Cionni, I., Hardiman, S. C., Haynes, P. H., Krüger, K., Kushner, P. J., Newman, P. A., Osprey, S. M., Perlwitz, J., Sigmond, M., Wang, L., Akiyoshi, H., Austin, J., Bekki, S., Baumgaertner, A., Braesicke, P., Brühl, C., Chipperfield, M., Dameris, M., Dhomse, S., Eyring, V., Garcia, R., Garny, H., Jöckel, P., Lamarque, J.-F., Marchand, M., Michou, M., Morgenstern, O., Nakamura, T., Pawson, S., Plummer, D., Pyle, J., Rozanov, E., Scinocca, J., Shepherd, T. G., Shibata, K., Smale, D., Teyssèdre, H., Tian, W., Waugh, D., and Yamashita, Y.: Multimodel climate and variability of the stratosphere, J. Geophys. Res.-Atmos., 116, D5, https://doi.org/10.1029/2010JD014995, 2011. a
Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a, b
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 CCMVal-2 and CCMI-1 models, Atmos. Chem. Phys., 18, 6699–6720, https://doi.org/10.5194/acp-18-6699-2018, 2018. a
Dietmüller, S., Garny, H., Eichinger, R., and Ball, W. T.: Analysis of recent lower-stratospheric ozone trends in chemistry climate models, Atmos. Chem. Phys., 21, 6811–6837, https://doi.org/10.5194/acp-21-6811-2021, 2021. a
Dlugokencky, E.: Global Monitoring Laboratory – Carbon Cycle Greenhouse Gases – Trends in Atmospheric Sulpher Hexaflouride, available at: https://www.esrl.noaa.gov/gmd/ccgg/trends_sf6/ (last access: 3 September 2020), 2005. a
Eichinger, R., Dietmüller, S., Garny, H., Šácha, P., Birner, T., Bönisch, H., Pitari, G., Visioni, D., Stenke, A., Rozanov, E., Revell, L., Plummer, D. A., Jöckel, P., Oman, L., Deushi, M., Kinnison, D. E., Garcia, R., Morgenstern, O., Zeng, G., Stone, K. A., and Schofield, R.: The influence of mixing on the stratospheric age of air changes in the 21st century, Atmos. Chem. Phys., 19, 921–940, https://doi.org/10.5194/acp-19-921-2019, 2019. a, b
Engel, A., Möbius, T., Bönisch, 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, d, e, f, g, h, i, j, k, l, m, n, o, p, q
Engel, A., Bönisch, H., Ullrich, M., Sitals, R., Membrive, O., Danis, F., and Crevoisier, C.: Mean age of stratospheric air derived from AirCore observations, Atmos. Chem. Phys., 17, 6825–6838, https://doi.org/10.5194/acp-17-6825-2017, 2017. a, b, c, d, e, f, g
Eyring, V., Waugh, D. W., Bodeker, G. E., Cordero, E., Akiyoshi, H., Austin, J., Beagley, S. R., Boville, B. A., Braesicke, P., Brühl, C., Butchart, N., Chipperfield, M. P., Dameris, M., Deckert, R., Deushi, M., Frith, S. M., Garcia, R. R., Gettelman, A., Giorgetta, M. A., Kinnison, D. E., Mancini, E., Manzini, E., Marsh, D. R., Matthes, S., Nagashima, T., Newman, P. A., Nielsen, J. E., Pawson, S., Pitari, G., Plummer, D. A., Rozanov, E., Schraner, M., Scinocca, J. F., Semeniuk, K., Shepherd, T. G., Shibata, K., Steil, B., Stolarski, R. S., Tian, W., and Yoshiki, M.: Multimodel projections of stratospheric ozone in the 21st century, J. Geophys. Res.-Atmos., 112, D16, https://doi.org/10.1029/2006JD008332, 2007. a
Fritsch, F., Garny, H., Engel, A., Bönisch, H., and Eichinger, R.: Sensitivity of age of air trends to the derivation method for non-linear increasing inert SF6, Atmos. Chem. Phys., 20, 8709–8725, https://doi.org/10.5194/acp-20-8709-2020, 2020. a, b, c, d, e
Garcia, R. R., Randel, W. J., and Kinnison, D. E.: On the Determination of Age of Air Trends from Atmospheric Trace Species, J. Atmos. Sci., 68, 139–154, https://doi.org/10.1175/2010JAS3527.1, 2011. a, b
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/acp-15-13161-2015, 2015. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o
Hall, T. M. and Waugh, D. W.: The influence of nonlocal chemistry on tracer distributions: Inferring the mean age of air from SF6, J. Geophys. Res., 103, 13327–13336, https://doi.org/10.1029/98JD00170, 1998. a
Harrison, J. J.: New infrared absorption cross sections for the infrared limb sounding of sulfur hexafluoride (SF6), J. Quant. Spectrosc. Ra., 254, 107202, https://doi.org/10.1016/j.jqsrt.2020.107202, 2020. a, b
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/acp-5-433-2005, 2005. a, b, c
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/gmd-3-717-2010, 2010. a, b, c, d, e
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änder-Hayn, 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/gmd-9-1153-2016, 2016. a, b, c, d, e, f, g, h, i, j, k, l
Kouznetsov, R., Sofiev, M., Vira, J., and Stiller, G.: Simulating age of air and the distribution of SF6 in the stratosphere with the SILAM model, Atmos. Chem. Phys., 20, 5837–5859, https://doi.org/10.5194/acp-20-5837-2020, 2020. a, b, c, d, e, f, g
Kovács, T., Feng, W., Totterdill, A., Plane, J. M. C., Dhomse, S., Gómez-Martí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 three-dimensional model, Atmos. Chem. Phys., 17, 883–898, https://doi.org/10.5194/acp-17-883-2017, 2017. a, b, c
Leedham Elvidge, E., 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 CF4, C2F6, C3F8, CHF3, HFC-125, HFC-227ea and SF6; 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/acp-18-3369-2018, 2018. a, b
Marsh, D. R., Mills, M. J., Kinnison, D. E., Lamarque, J.-F., Calvo, N., and Polvani, L. M.: Climate Change from 1850 to 2005 Simulated in CESM1(WACCM), J. Climate, 26, 7372–7391, https://doi.org/10.1175/JCLI-D-12-00558.1, 2013. 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., Saint-Martin, 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/gmd-10-639-2017, 2017. a
Ploeger, F., Abalos, M., Birner, T., Konopka, P., Legras, B., Müller, R., and Riese, M.: Quantifying the effects of mixing and residual circulation on trends of stratospheric mean age of air, Geophys. Res. Lett., 42, 2047–2054, https://doi.org/10.1002/2014GL062927, 2014GL062927, 2015. a
Plumb, I. C., Vohralik, P. F., and Ryan, K. R.: Normalization of correlations for atmospheric species with chemical loss, J. Geophys. Res.-Atmos., 104, 11723–11732, https://doi.org/10.1029/1999JD900014, 1999. a
Ravishankara, A. R., Solomon, S., Turnipseed, A. A., and Warren, R. F.: Atmospheric lifetimes of long-lived halogenated species, Science, 17, 194–199, https://doi.org/10.1126/science.259.5092.194, 1993. a, b, c
Ray, E. A., Moore, F. L., Rosenlof, K. H., Davis, S. M., Sweeney, C., Tans, P., Wang, T., Elkins, J. W., Bönisch, H., Engel, A., Sugawara, S., Nakazawa, T., and Aoki, S.: Improving stratospheric transport trend analysis based on SF6 and CO2 measurements, J. Geophys. Res.-Atmos., 14, 110–128, https://doi.org/10.1002/2014JD021802, 2014. a, b, c
Ray, E. A., Moore, F., Elkins, J. W., Rosenlof, K., Laube, J., Röckmann, T., Marsh, D., and Andrews, A.: Quantification of the SF6 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, f, g
Reddmann, T., Ruhnke, R., and Kouker, W.: Three-dimensional model simulations of SF6 with mesospheric chemistry, J. Geophys. Res.-Atmos., 106, 14525–14537, https://doi.org/10.1029/2000JD900700, 2001. a, b, c, d, e, f, g, h, i, j, k, l
Roeckner, E., Bäuml, G., Bonaventura, L., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kirchner, I., Kornblueh, L., Manzini, E., Rhodin, A., Schlese, U., Schulzweida, U., and Tompkins, A.: The atmospheric general circulation model ECHAM 5. PART I: model description, Max Planck Institute for Meteorology Report, Max Planck Institute for Meteorology, 349, 2003. 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
SPARC: SPARC CCMVal Report on the Evaluation of Chemistry-Climate Models, Tech. rep., SPARC, available at: http://www.sparc-climate.org/publications/sparc-reports/ (last access: 15 November 2021), 2010. a
Stiller, G.: MIPAS monthly zonal mean Age of Air data, based on MIPAS SF6 retrievals, processed at KIT-IMK, Karlsruhe, 12.11.25; LK 01, KIT [data set], https://doi.org/10.5445/IR/1000139453, 2021. a, b, c, d, e, f, g
Stiller, G., Harrison, J., Haenel, F., Glatthor, N., Kellmann, S., and von Clarmann, T.: Improved global distributions of SF6 and mean age of stratospheric air by use of new spectroscopic data, in: EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-2660, https://doi.org/10.5194/egusphere-egu2020-2660, 2020. a, b, c, d, e
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ópez-Puertas, M., and Engel, A.: Global distribution of mean age of stratospheric air from MIPAS SF6 measurements, Atmos. Chem. Phys., 8, 677–695, https://doi.org/10.5194/acp-8-677-2008, 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ópez-Puertas, 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/acp-12-3311-2012, 2012. a, b, c, d, e, f, g
Stiller, G. P., Fierli, F., Ploeger, F., Cagnazzo, C., Funke, B., Haenel, F. J., Reddmann, T., Riese, M., and von Clarmann, T.: Shift of subtropical transport barriers explains observed hemispheric asymmetry of decadal trends of age of air, Atmos. Chem. Phys., 17, 11177–11192, https://doi.org/10.5194/acp-17-11177-2017, 2017. a, b
Stiller, G. P., Harrison, J. J., Haenel, F. J., Glatthor, N., and Kellmann, S.: Technical note: Improved global distributions of SF6 and mean age of stratospheric air by use of new spectroscopic data, Atmos. Chem. Phys. Discuss., in preparation, 2022. a
Totterdill, A., Kovács, T., Gómez Martín, J. C., Feng, W., and Plane, J. M. C.: Mesospheric Removal of Very Long-Lived Greenhouse Gases SF6 and CFC-115 by Metal Reactions, Lyman-α Photolysis, and Electron Attachment, J. Phys. Chem. A, 119, 2016–2025, https://doi.org/10.1021/jp5123344, pMID: 25647411, 2015. a
Volk, C. M., Elkins, J. W., Fahey, D. W., Dutton, G. S., Gilligan, J. M., Loewenstein, M., Podolske, J. R., Chan, K. R., and Gunson, M. R.: Evaluation of source gas lifetimes from stratospheric observations, J. Geophys. Res.-Atmos., 102, 25543–25564, https://doi.org/10.1029/97JD02215, 1997. a
von Clarmann, T., Stiller, G., Grabowski, U., Eckert, E., and Orphal, J.: Technical Note: Trend estimation from irregularly sampled, correlated data, Atmos. Chem. Phys., 10, 6737–6747, https://doi.org/10.5194/acp-10-6737-2010, 2010. a
WMO: Atmospheric ozone 1985 – volume I–III: assessment of our understanding of the processes controlling its present distribution and change, World Meteorological Organization, 1986. a