The impact of sulfur hexafluoride (SF₆) 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 (SF$_{6}$), a fairly long-lived gas with a near-linear increase in emissions during recent decades. However, SF$_{6}$ has chemical sinks in the mesosphere that are not considered in most model studies. In this study, we explicitly compute the chemical SF$_{6}$ 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 SF$_{6}$ 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 SF$_{6}$ 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 SF$_{6}$ sinks. We illustrate that even a constant SF$_{6}$ destruction rate causes a positive trend in the derived AoA, as the amount of depleted SF$_{6}$ scales with increasing SF$_{6}$ 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 SF$_{6}$ 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 SF$_{6}$ loss and deduce that a linear correction might be applicable to values of AoA of up to 4 years