Explore the vast range of titles available.

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

Rising emissions of anthropogenic greenhouse gases (GHG) have led to tropospheric warming and stratospheric cooling over recent decades. As a thermodynamic consequence, the troposphere has expanded and the rise of the tropopause, the boundary between the troposphere and stratosphere, has been suggested as one of the most robust fingerprints of anthropogenic climate change. Conversely, at altitudes above ∼55 km (in the mesosphere and thermosphere) observational and modeling evidence indicates a downward shift of the height of pressure levels or decreasing density at fixed altitudes. The layer in between, the stratosphere, has not been studied extensively with respect to changes of its global structure. Here we show that this atmospheric layer has contracted substantially over the last decades, and that the main driver for this are increasing concentrations of GHG. Using data from coupled chemistry-climate models we show that this trend will continue and the mean climatological thickness of the stratosphere will decrease by 1.3 km following representative concentration pathway 6.0 by 2080. We also demonstrate that the stratospheric contraction is not only a response to cooling, as changes in both tropopause and stratopause pressure contribute. Moreover, its short emergence time (less than 15 years) makes it a novel and independent indicator of GHG induced climate change.

Energy and momentum deposition from planetary-scale Rossby waves as well as from small-scale gravity waves (GWs) largely control stratospheric dynamics. Interactions between these different wave types, however, complicate the quantification of their individual contribution to the overall dynamical state of the middle atmosphere. In state-of-the-art general circulation models (GCMs), the majority of the GW spectrum cannot be resolved and therefore has to be parameterised. This is commonly implemented in two discrete schemes, one for GWs that originate from flow over orographic obstacles and one for all other kinds of GWs (non-orographic GWs). In this study, we attempt to gain a deeper understanding of the interactions of resolved with parameterised wave driving and of their influence on the stratospheric zonal winds and on the Brewer–Dobson circulation (BDC). For this, we set up a GCM time slice experiment with two sensitivity simulations: one without orographic GWs and one without non-orographic GWs. Our findings include an acceleration of the polar vortices, which has historically been one of the main reasons for including explicit GW parameterisations in GCMs. Further, we find inter-hemispheric differences in BDC changes when omitting GWs that can be explained by wave compensation and amplification effects. These are partly evoked through local changes in the refractive properties of the atmosphere caused by the omitted GW drag and a thereby increased planetary wave propagation. However, non-local effects on the flow can act to suppress vertical wave fluxes into the stratosphere for a very strong polar vortex. Moreover, we study mean age of stratospheric air to investigate the impact of missing GWs on tracer transport. On the basis of this analysis, we suggest that the larger ratio of planetary waves to GWs leads to enhanced horizontal mixing, which can have a large impact on stratospheric tracer distributions.

The columnar approach of gravity wave (GW) parameterisations in weather and climate models has been identified as a potential reason for dynamical biases in middle-atmospheric dynamics. For example, GW momentum flux (GWMF) discrepancies between models and observations at 60∘ S arising through the lack of horizontal orographic GW propagation are suspected to cause deficiencies in representing the Antarctic polar vortex. However, due to the decomposition of the model domains onto different computing tasks for parallelisation, communication between horizontal grid boxes is computationally extremely expensive, making horizontal propagation of GWs unfeasible for global chemistry–climate simulations.To overcome this issue, we present a simplified solution to approximate horizontal GW propagation through redistribution of the GWMF at one single altitude by means of tailor-made redistribution maps. To generate the global redistribution maps averaged for each grid box, we use a parameterisation describing orography as a set of mountain ridges with specified location, orientation and height combined with a ray-tracing model describing lateral propagation of so-generated mountain waves. In the global chemistry–climate model (CCM) EMAC (ECHAM MESSy Atmospheric Chemistry), these maps then allow us to redistribute the GW momentum flux horizontally at one level, obtaining an affordable overhead of computing resources. The results of our simulations show GWMF and drag patterns that are horizontally more spread out than with the purely columnar approach; GWs are now also present above the ocean and regions without mountains. In this paper, we provide a detailed description of how the redistribution maps are computed and how the GWMF redistribution is implemented in the CCM. Moreover, an analysis shows why 15 km is the ideal altitude for the redistribution. First results with the redistributed orographic GWMF provide clear evidence that the redistributed GW drag in the Southern Hemisphere has the potential to modify and improve Antarctic polar vortex dynamics, thereby paving the way for enhanced credibility of CCM simulations and projections of polar stratospheric ozone.

Climate models consistently predict an acceleration of the Brewer–Dobson circulation (BDC) due to climate change in the 21st century. However, the strength of this acceleration varies considerably among individual models, which constitutes a notable source of uncertainty for future climate projections. To shed more light upon the magnitude of this uncertainty and on its causes, we analyse the stratospheric mean age of air (AoA) of 10 climate projection simulations from the Chemistry-Climate Model Initiative phase 1 (CCMI-I), covering the period between 1960 and 2100. In agreement with previous multi-model studies, we find a large model spread in the magnitude of the AoA trend over the simulation period. Differences between future and past AoA are found to be predominantly due to differences in mixing (reduced aging by mixing and recirculation) rather than differences in residual mean transport. We furthermore analyse the mixing efficiency, a measure of the relative strength of mixing for given residual mean transport, which was previously hypothesised to be a model constant. Here, the mixing efficiency is found to vary not only across models, but also over time in all models. Changes in mixing efficiency are shown to be closely related to changes in AoA and quantified to roughly contribute 10 % to the long-term AoA decrease over the 21st century. Additionally, mixing efficiency variations are shown to considerably enhance model spread in AoA changes. To understand these mixing efficiency variations, we also present a consistent dynamical framework based on diffusive closure, which highlights the role of basic state potential vorticity gradients in controlling mixing efficiency and therefore aging by mixing.

Aviation emissions are responsible for climate impacts through both carbon dioxide emissions and other emissions, in particular, of nitrogen oxides, water vapour, particulates, and contrail formation. In December 2022, the European Commission, Parliament and Council agreed to revise the European Union Emission Trading System for aviation. As such, from January 1, 2025, aircraft operators must monitor non-carbon dioxide climate effects, but suitable metrics for climate impact, handling of uncertainties and practical implementation are still under discussion or at least heavily debated. In this perspective, we propose a procedure for how to include non-carbon dioxide aviation effects into political frameworks. The main goal must be to create incentives for climate change mitigation for the aviation industry. Uncertainties in atmospheric processes need to be appropriately incorporated to minimise risk, and pilot projects are required to test implementation capabilities. Analysing risk, employing consistent monitoring, and determining economic effects will provide scientific grounds for including non-carbon dioxide effects in the European Union Emission Trading System. For the inclusion of aviation non-carbon dioxide climate effects in the European Union Emission Trading System, uncertainties in atmospheric processes can be addressed through careful calculation and risk assessment.

Climate model simulations show an acceleration of the Brewer–Dobson circulation (BDC) in response to climate change. While the general mechanisms for the BDC strengthening are widely understood, there are still open questions concerning the influence of the details of the wave driving. Mean age of stratospheric air (AoA) is a useful transport diagnostic for assessing changes in the BDC. Analyzing AoA from a subset of Chemistry–Climate Model Initiative part 1 climate projection simulations, we find a remarkable agreement between most of the models in simulating the largest negative AoA trends in the extratropical lower to middle stratosphere of both hemispheres (approximately between 20 and 25 geopotential kilometers (gpkm) and 20–50∘ N and S). We show that the occurrence of AoA trend minima in those regions is directly related to the climatological AoA distribution, which is sensitive to an upward shift of the circulation in response to climate change. Also other factors like a reduction of aging by mixing (AbM) and residual circulation transit times (RCTTs) contribute to the AoA distribution changes by widening the AoA isolines. Furthermore, we analyze the time evolution of AbM and RCTT trends in the extratropics and examine the connection to possible drivers focusing on local residual circulation strength, net tropical upwelling and wave driving. However, after the correction for a vertical shift of pressure levels, we find only seasonally significant trends of residual circulation strength and zonal mean wave forcing (resolved and unresolved) without a clear relation between the trends of the analyzed quantities. This indicates that additional causative factors may influence the AoA, RCTT and AbM trends. In this study, we postulate that the shrinkage of the stratosphere has the potential to influence the RCTT and AbM trends and thereby cause additional AoA changes over time.

Recent observations show a significant decrease in lower-stratospheric (LS) ozone concentrations in tropical and mid-latitude regions since 1998. By analysing 31 chemistry climate model (CCM) simulations performed for the Chemistry Climate Model Initiative (CCMI; Morgenstern et al., 2017), we find a large spread in the 1998–2018 trend patterns between different CCMs and between different realizations performed with the same CCM. The latter in particular indicates that natural variability strongly influences LS ozone trends. However none of the model simulations reproduce the observed ozone trend structure of coherent negative trends in the LS. In contrast to the observations, most models show an LS trend pattern with negative trends in the tropics (20∘ S–20∘ N) and positive trends in the northern mid-latitudes (30–50∘ N) or vice versa. To investigate the influence of natural variability on recent LS ozone trends, we analyse the sensitivity of observational trends and the models' trend probability distributions for varying periods with start dates from 1995 to 2001 and end dates from 2013 to 2019. Generally, modelled and observed LS trends remain robust for these different periods; however observational data show a change towards weaker mid-latitude trends for certain periods, likely forced by natural variability. Moreover we show that in the tropics the observed trends agree well with the models' trend distribution, whereas in the mid-latitudes the observational trend is typically an extreme value of the models' distribution. We further investigate the LS ozone trends for extended periods reaching into the future and find that all models develop a positive ozone trend at mid-latitudes, and the trends converge to constant values by the period that spans 1998–2060. Inter-model correlations between ozone trends and transport-circulation trends confirm the dominant role of greenhouse gas (GHG)-driven tropical upwelling enhancement on the tropical LS ozone decrease. Mid-latitude ozone, on the other hand, appears to be influenced by multiple competing factors: an enhancement in the shallow branch decreases ozone, while an enhancement in the deep branch increases ozone, and, furthermore, mixing plays a role here too. Sensitivity simulations with fixed forcing of GHGs or ozone-depleting substances (ODSs) reveal that the GHG-driven increase in circulation strength does not lead to a net trend in LS mid-latitude column ozone. Rather, the positive ozone trends simulated consistently in the models in this region emerge from the decline in ODSs, i.e. the ozone recovery. Therefore, we hypothesize that next to the influence of natural variability, the disagreement of modelled and observed LS mid-latitude ozone trends could indicate a mismatch in the relative role of the response of ozone to ODS versus GHG forcing in the models.

Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis.

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis.Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.