Search

Solar coronal mass ejections can accelerate charged particles, mostly protons, to high energies, causing solar proton events (SPEs). Such energetic particles can precipitate upon the Earth's atmosphere, mostly in polar regions because of geomagnetic shielding. Here, SPE-induced chlorine activation due to ion chemistry can occur, and the activated chlorine depletes ozone in the polar middle atmosphere. We use the state-of-the-art 1D stacked-box Exoplanetary Terrestrial Ion Chemistry (ExoTIC) model of atmospheric ion and neutral composition to investigate such events in the Northern Hemisphere (NH). The Halloween SPE that occurred in late October 2003 is used as a test field for our study. This event has been extensively studied before using different 3D models and satellite observations. Our main purpose is to use such a large event that has been recorded by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on the Environmental Satellite (ENVISAT) to evaluate the performance of the ion chemistry model. Sensitivity tests were carried out for different model settings with a focus on the chlorine species of HOCl and ClONO2 as well as O3 and reactive nitrogen, NOy. The model simulations were performed in the Northern Hemisphere at a high latitude of 67.5∘ N, inside the polar cap. Comparison of the simulated effects against MIPAS observations for the Halloween SPE revealed rather good temporal agreement, also in terms of altitude range for HOCl, O3 and NOy. For ClONO2, good agreement was found in terms of altitude range. The model showed ClONO2 enhancements after the peak of the event. The best model setting was the one with full ion chemistry where O(1D) was set to photo-chemical equilibrium. HOCl and ozone changes are very well reproduced by the model, especially for nighttime. HOCl was found to be the main active chlorine species under nighttime conditions, resulting in an increase of more than 0.2 ppbv. Further, ClONO2 enhancements of 0.2–0.3 ppbv have been observed during both daytime and nighttime. Model settings that compared best with MIPAS observations were applied to an extreme solar event that occurred in AD 775, presumably once in a 1000-year event. With the model applied to this scenario, an assessment can be made about what is to be expected at worst for the effects of a SPE on the middle atmosphere, concentrating on the effects of ion chemistry compared to crude parameterizations. Here, a systematic analysis comparing the impact of the Halloween SPE and the extreme event on the Earth's middle atmosphere is presented. As seen from the model simulations, both events were able to perturb the polar stratosphere and mesosphere with a high production of NOy and HOx. Longer-lasting and stronger stratospheric ozone loss was seen for the extreme event. A qualitative difference between the two events and a long-lasting impact on HOCl and HCl for the extreme event were found. Chlorine ion chemistry contributed to stratospheric ozone losses of 2.4 % for daytime and 10 % for nighttime during the Halloween SPE, as seen with time-dependent ionization rates applied to the model. Furthermore, while comparing the Halloween SPE and the extreme scenario, with ionization rate profiles applied just for the event day, the inclusion of chlorine ion chemistry added ozone losses of 10 % and 20 % respectively.

The number of large urban agglomerations is steadily increasing worldwide. At a local scale, their emissions lead to air pollution, directly affecting people's health. On a global scale, their emissions lead to an increase of greenhouse gases, affecting climate. In this context, in 2017 and 2018, the airborne campaign EMeRGe (Effect of Megacities on the transport and transformation of pollutants on the Regional to Global scales) investigated emissions of European and Asian major population centres (MPCs) to improve the understanding and predictability of pollution outflows. Here, we present two methods to identify and characterise pollution outflows probed during EMeRGe. First, we use a set of volatile organic compounds (VOCs) as chemical tracers to characterise air masses by specific source signals, i.e. benzene from anthropogenic pollution of targeted regions, acetonitrile from biomass burning (BB, primarily during EMeRGe-Asia), and isoprene from fresh biogenic signals (primarily during EMeRGe-Europe. Second, we attribute probed air masses to source regions and estimate their individual contribution by constructing and applying a simple emission uptake scheme for the boundary layer which combines FLEXTRA back trajectories and EDGAR carbon monoxide (CO) emission rates (acronyms are provided in the Appendix). During EMeRGe-Europe, we identified anthropogenic pollution outflows from northern Italy, southern Great Britain, the Belgium–Netherlands–Ruhr (BNR) area and the Iberian Peninsula. Additionally, our uptake scheme indicates significant long-range transport of pollution from the USA and Canada. During EMeRGe-Asia, the pollution outflow is dominated by sources in China and Taiwan, but BB signals from Southeast Asia and India contribute as well. Outflows of pre-selected MPC targets are identified in less than 20 % of the sampling time, due to restrictions in flight planning and constraints of the measurement platform itself. Still, EMeRGe combines in a unique way near- and far-field measurements, which show signatures of local and distant sources, transport and conversion fingerprints, and complex air mass compositions. Our approach provides a valuable classification and characterisation of the EMeRGe dataset, e.g. for BB and anthropogenic influence of potential source regions and paves the way for a more comprehensive analysis and various model studies.

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

Impacts of chlorinated very short-lived substances (Cl-VSLS) on stratospheric chlorine budget over the first two decades of the 21st century are assessed using the Met Office’s Unified Model coupled to the United Kingdom Chemistry and Aerosol (UM-UKCA) chemistry–climate model; this constitutes the most up-to-date assessment and the first study to simulate Cl-VSLS impacts using a whole atmosphere chemistry–climate model. We examine the Cl-VSLS responses using a small ensemble of free-running simulations and two pairs of integrations where the meteorology was “nudged” to either ERA5 or ERA-Interim reanalysis. The stratospheric chlorine source gas injection due to Cl-VSLS estimated from the free-running integrations doubled from ∼40 ppt Cl injected into the stratosphere in 2000 to ∼80 ppt Cl injected in 2019. Combined with chlorine product gas injection, the integrations showed a total of ∼130 ppt Cl injected into the stratosphere in 2019 due to Cl-VSLS. The use of the nudged model significantly increased the abundance of Cl-VSLS simulated in the lower stratosphere relative to the free-running model. Averaged over 2010–2018, simulations nudged to ERAI-Interim and ERA5 showed 20 ppt (i.e. a factor of 2) and 10 ppt (i.e. ∼50 %) more Cl, respectively, in the tropical lower stratosphere at 20 km in the form of Cl-VSLS source gases compared to the free-running case. These differences can be explained by the corresponding differences in the speed of the large-scale circulation. The results illustrate the strong dependence of the simulated stratospheric Cl-VSLS levels on the model dynamical fields. In UM-UKCA, this corresponds to the choice between free-running versus nudged set-up, and to the reanalysis dataset used for nudging. Temporal changes in Cl-VSLS are found to have significantly impacted recent HCl and COCl2 trends in the model. In the tropical lower stratosphere, the inclusion of Cl-VSLS reduced the magnitude of the negative HCl and COCl2 trends (e.g. from ∼-8%(HCl)/decade and ∼-4 ppt(COCl2)/decade at ∼20 km to ∼-6%(HCl)/decade and ∼ −2 ppt(COCl2)/decade in the free running simulations) and gave rise to positive tropospheric trends in both tracers. In the tropics, both the free-running and nudged integrations with Cl-VSLS included compared much better to the observed trends from the ACE-FTS satellite record than the analogous simulations without Cl-VSLS. Since observed HCl trends provide information on the evolution of total stratospheric chlorine and, thus, the effectiveness of the Montreal Protocol, our results demonstrate that Cl-VSLS are a confounding factor in the interpretation of such data and should be factored into future analysis. Unlike the nudged model runs, the ensemble mean free-running integrations did not reproduce the hemispheric asymmetry in the observed mid-latitude HCl and COCl2 trends related to short-term dynamical variability. The individual ensemble members also showed a considerable spread of the diagnosed tracer trends, illustrating the role of natural interannual variability in modulating the diagnosed responses and the need for caution when interpreting both model and observed tracer trends derived over a relatively short time period.

Water vapour and ozone are important for the thermal and radiative balance of the upper troposphere (UT) and lowermost stratosphere (LMS). Both species are modulated by transport processes. Chemical and microphysical processes affect them differently. Thus, representing the different processes and their interactions is a challenging task for dynamical cores, chemical modules and microphysical parameterisations of state-of-the-art atmospheric model components. To test and improve the models, high-resolution measurements of the UT–LMS are required. Here, we use measurements taken in a flight of the GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) instrument on HALO (High Altitude and LOng Range Research Aircraft). The German research aircraft HALO performed a research flight on 26 February 2016 that covered deeply subsided air masses of the aged 2015/16 Arctic vortex, high-latitude LMS air masses, a highly textured region affected by troposphere-to-stratosphere exchange and high-altitude cirrus clouds. Therefore, it provides a challenging multifaceted case study for comparing GLORIA observations with state-of-the-art atmospheric model simulations in a complex UT–LMS region at a late stage of the Arctic winter 2015/16. Using GLORIA observations in this manifold scenario, we test the ability of the numerical weather prediction (NWP) model ICON (ICOsahedral Nonhydrostatic) with the extension ART (Aerosols and Reactive Trace gases) and the chemistry–climate model (CCM) EMAC (ECHAM5/MESSy Atmospheric Chemistry – fifth-generation European Centre Hamburg general circulation model/Modular Earth Submodel System) to model the UT–LMS composition of water vapour (H2O), ozone (O3), nitric acid (HNO3) and clouds. Within the scales resolved by the respective model, we find good overall agreement of both models with GLORIA. The applied high-resolution ICON-ART set-up involving an R2B7 nest (local grid refinement with a horizontal resolution of about 20 km), covering the HALO flight region, reproduces mesoscale dynamical structures well. Narrow moist filaments in the LMS observed by GLORIA at tropopause gradients in the context of a Rossby wave breaking event and in the vicinity of an occluded Icelandic low are clearly reproduced by the model. Using ICON-ART, we show that a larger filament in the west was transported horizontally into the Arctic LMS in connection with a jet stream split associated with poleward breaking of a cyclonically sheared Rossby wave. Further weaker filaments are associated with an older tropopause fold in the east. Given the lower resolution (T106) of the nudged simulation of the EMAC model, we find that this model also reproduces these features well. Overall, trace gas mixing ratios simulated by both models are in a realistic range, and major cloud systems observed by GLORIA are mostly reproduced. However, we find both models to be affected by a well-known systematic moist bias in the LMS. Further biases are diagnosed in the ICON-ART O3, EMAC H2O and EMAC HNO3 distributions. Finally, we use sensitivity simulations to investigate (i) short-term cirrus cloud impacts on the H2O distribution (ICON-ART), (ii) the overall impact of polar winter chemistry and microphysical processing on O3 and HNO3 (ICON-ART and EMAC), (iii) the impact of the model resolution on simulated parameters (EMAC), and (iv) consequences of scavenging processes by cloud particles (EMAC). We find that changing the horizontal model resolution results in notable systematic changes for all species in the LMS, while scavenging processes play a role only in the case of HNO3. We discuss the model biases and deficits found in this case study that potentially affect forecasts and projections (adversely) and provide suggestions for further model improvements.

We present a global and multi-annual space-borne dataset of tropospheric {H2O, δD} pairs that is based on radiance measurements from the nadir thermal infrared sensor IASI (Infrared Atmospheric Sounding Interferometer) on board the Metop satellites of EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites). This dataset is an a posteriori processed extension of the MUSICA (MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water) IASI full product dataset as presented in Schneider et al. (2021b). From the independently retrieved H2O and δD proxy states, their a priori settings and constraints, and their error covariances provided by the IASI full product dataset, we generate an optimal estimation product for pairs of H2O and δD. Here, this standard MUSICA method for deriving {H2O, δD} pairs is extended using an a posteriori reduction of the constraints for improving the retrieval sensitivity at dry conditions. By applying this improved water isotopologue post-processing for all cloud-free MUSICA IASI retrievals, this yields a {H2O, δD} pair dataset for the whole period from October 2014 to December 2020 with global coverage twice per day (local morning and evening overpass times). In total, the dataset covers more than 1500 million individually processed observations. The retrievals are most sensitive to variations in {H2O, δD} pairs within the free troposphere, with up to 30 % of all retrievals containing vertical profile information in the {H2O, δD} pair product. After applying appropriate quality filters, the largest number of reliable pair data arises for tropical and subtropical summer regions, but higher latitudes also show a considerable amount of reliable data. Exemplary time series over the tropical Atlantic and West Africa are chosen to illustrate the potential of the MUSICA IASI {H2O, δD} pair data for atmospheric moisture pathway studies. Furthermore, in order to facilitate the application of this rather comprehensive MUSICA IASI {H2O, δD} pair dataset (referred to as Level-2), we additionally provide the data in a re-gridded and simplified format (Level-3) with focus on the quality-filtered {H2O, δD} pairs in the free troposphere. A technical documentation for guiding the use of both datasets is attached as the Supplement. Finally, the Level-2 dataset is referenced with the DOI https://doi.org/10.35097/415 (Diekmann et al., 2021a) and the Level-3 dataset with DOI https://doi.org/10.35097/495 (Diekmann et al., 2021b).

Aerosols, ozone, surface reflectivity, and clouds are, among other factors, important for the modulation of UV radiation levels at the Earth’s surface. In this study, these variables were extracted from climate model integrations that contributed to CMIP6 and from MACC global reanalysis provided by the CAMS. From the 1950s until the end of the 21st century and for various shared socioeconomic pathways considered by CMIP6, we conclude that total ozone will increase globally in 2100 by up to 11% under SSP5–8.5 relative to the 1950s, while under SSP1–2.6, ozone is not projected to recover to pre-ozone depletion levels. AOD at 550 nm shows reductions in 2090–2100 over the Northern Hemisphere of up to −0.38. Compared to CAMS, CMIP6 models show a general overestimation for total ozone of up to 2.5% in extra-polar regions for models with interactive chemistry. This implies systematically lower UV radiation from models based on these projections.

Polar stratospheric clouds (PSCs) are a driver for ozone depletion in the lower polar stratosphere. They provide surface for heterogeneous reactions activating chlorine and bromine reservoir species during the polar night. The large-scale effects of PSCs are represented by means of parameterisations in current global chemistry–climate models, but one process is still a challenge: the representation of PSCs formed locally in conjunction with unresolved mountain waves. In this study, we investigate direct simulations of PSCs formed by mountain waves with the ICOsahedral Nonhydrostatic modelling framework (ICON) with its extension for Aerosols and Reactive Trace gases (ART) including local grid refinements (nesting) with two-way interaction. Here, the nesting is set up around the Antarctic Peninsula, which is a well-known hot spot for the generation of mountain waves in the Southern Hemisphere. We compare our model results with satellite measurements of PSCs from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and gravity wave observations of the Atmospheric Infrared Sounder (AIRS). For a mountain wave event from 19 to 29 July 2008 we find similar structures of PSCs as well as a fairly realistic development of the mountain wave between the satellite data and the ICON-ART simulations in the Antarctic Peninsula nest. We compare a global simulation without nesting with the nested configuration to show the benefits of adding the nesting. Although the mountain waves cannot be resolved explicitly at the global resolution used (about 160 km), their effect from the nested regions (about 80 and 40 km) on the global domain is represented. Thus, we show in this study that the ICON-ART model has the potential to bridge the gap between directly resolved mountain-wave-induced PSCs and their representation and effect on chemistry at coarse global resolutions.

An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH3CCl3) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom-up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long-lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH-SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long-term trend and emissions derived from the measured hemispheric gradient, the combination of HFC-32 (CH2F2), HFC-134a (CH2FCF3, HFC-152a (CH3CHF2), and HCFC-22 (CHClF2), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF.

An important source of polar stratospheric clouds (PSCs), which play a crucial role in controlling polar stratospheric ozone depletion, is the temperature fluctuations induced by mountain waves. These enable stratospheric temperatures to fall below the threshold value for PSC formation in regions of negative temperature perturbations or cooling phases induced by the waves even if the synoptic-scale temperatures are too high. However, this formation mechanism is usually missing in global chemistry–climate models because these temperature fluctuations are neither resolved nor parameterised. Here, we investigate in detail the episodic and localised wintertime stratospheric cooling events produced over the Antarctic Peninsula by a parameterisation of mountain-wave-induced temperature fluctuations inserted into a 30-year run of the global chemistry–climate configuration of the UM-UKCA (Unified Model – United Kingdom Chemistry and Aerosol) model. Comparison of the probability distribution of the parameterised cooling phases with those derived from climatologies of satellite-derived AIRS brightness temperature measurements and high-resolution radiosonde temperature soundings from Rothera Research Station on the Antarctic Peninsula shows that they broadly agree with the AIRS observations and agree well with the radiosonde observations, particularly in both cases for the “cold tails” of the distributions. It is further shown that adding the parameterised cooling phase to the resolved and synoptic-scale temperatures in the UM-UKCA model results in a considerable increase in the number of instances when minimum temperatures fall below the formation temperature for PSCs made from ice water during late austral autumn and early austral winter and early austral spring, and without the additional cooling phase the temperature rarely falls below the ice frost point temperature above the Antarctic Peninsula in the model. Similarly, it was found that the formation potential for PSCs made from ice water was many times larger if the additional cooling is included. For PSCs made from nitric acid trihydrate (NAT) particles it was only during October that the additional cooling is required for temperatures to fall below the NAT formation temperature threshold (despite more NAT PSCs occurring during other months). The additional cooling phases also resulted in an increase in the surface area density of NAT particles throughout the winter and early spring, which is important for chlorine activation. The parameterisation scheme was finally shown to make substantial differences to the distribution of total column ozone during October, resulting from a shift in the position of the polar vortex.

The annual minimum Antarctic sea ice extent (SIE) in February 2022 hits a record low in the satellite era, with less than 2 million square kilometres observed on 25 February 2022, contrasting with the slightly positive trend in the Antarctic SIE prior to 2014. However, the preceding Amundsen Sea Low (ASL) in austral spring 2021 was the deepest since 1950. According to a linear regression model, the very low ASL contributed about 60% to the record low SIE in 2022. This study further investigates the underlying mechanism. The investigation of the lagged impact of the ASL on Antarctic SIE is based on observational data and state-of-the-art simulations. We found that (a) the deepened ASL associated with strengthened southerly winds accelerates the sea ice export away from the western Antarctic continent in spring, leading to the expansion of coastal polynyas (open water areas); (b) through the positive ice-ocean albedo feedback, the lack of the sea ice off the coastline enhances solar heating in the upper ocean and further sea ice melting in summer can occur. Specifically, in spring 2021, the deepest ASL is accompanied by a large sea-ice area flux of about 17.6 × 10 ^3 km ^2 across 70° S over the Ross Sea in October and November, contributing to a significant increase in net surface radiation of 20–30 W m ^−2 and upper ocean warming of about 0.5 K in summer. Therefore, the deepened ASL in spring 2021 plays a crucial role for the record low Antarctic SIE in February 2022. In addition, it is found that both the La Niña conditions and the strong stratospheric polar vortex contributed significantly to the very strong ASL in 2021. Currently, nearly 2/3 of Earth system models in the Coupled Model Intercomparison Project Phase 6 have difficulties capturing the relationship between the ASL and the Antarctic SIE.

The Arctic winter 2015–2016 was characterized by exceptionally low stratospheric temperatures, favouring the formation of polar stratospheric clouds (PSCs) from mid-December until the end of February down to low stratospheric altitudes. Observations by GLORIA (Gimballed Limb Observer for Radiance Imaging of the Atmosphere) on HALO (High Altitude and LOng range research aircraft) during the PGS (POLSTRACC–GW-LCYCLE II–SALSA) campaign from December 2015 to March 2016 allow the investigation of the influence of denitrification on the lowermost stratosphere (LMS) with a high spatial resolution. Two-dimensional vertical cross sections of nitric acid (HNO3) along the flight track and tracer–tracer correlations derived from the GLORIA observations document detailed pictures of wide-spread nitrification of the Arctic LMS during the course of an entire winter. GLORIA observations show large-scale structures and local fine structures with enhanced absolute HNO3 volume mixing ratios reaching up to 11 ppbv at altitudes of 13 km in January and nitrified filaments persisting until the middle of March. Narrow coherent structures tilted with altitude of enhanced HNO3, observed in mid-January, are interpreted as regions recently nitrified by sublimating HNO3-containing particles. Overall, extensive nitrification of the LMS between 5.0 and 7.0 ppbv at potential temperature levels between 350 and 380 K is estimated. The GLORIA observations are compared with CLaMS (Chemical Lagrangian Model of the Stratosphere) simulations. The fundamental structures observed by GLORIA are well reproduced, but differences in the fine structures are diagnosed. Further, CLaMS predominantly underestimates the spatial extent of HNO3 maxima derived from the GLORIA observations as well as the overall nitrification of the LMS. Sensitivity simulations with CLaMS including (i) enhanced sedimentation rates in case of ice supersaturation (to resemble ice nucleation on nitric acid trihydrate (NAT)), (ii) a global temperature offset, (iii) modified growth rates (to resemble aspherical particles with larger surfaces) and (iv) temperature fluctuations (to resemble the impact of small-scale mountain waves) slightly improved the agreement with the GLORIA observations of individual flights. However, no parameter could be isolated which resulted in a general improvement for all flights. Still, the sensitivity simulations suggest that details of particle microphysics play a significant role for simulated LMS nitrification in January, while air subsidence, transport and mixing become increasingly important for the simulated HNO3 distributions towards the end of the winter.

The atmospheric response to the 11-year solar cycle is separated into the contributions from changes in direct radiative heating and photolysis rates using specially designed sensitivity simulations with the UM-UKCA (Unified Model coupled to the United Kingdom Chemistry and Aerosol model) chemistry–climate model. We perform a number of idealised time-slice experiments under perpetual solar maximum (SMAX) and minimum conditions (SMIN), and we find that contributions from changes in direct heating and photolysis rates are both important for determining the stratospheric shortwave heating, temperature and ozone responses to the amplitude of the 11-year solar cycle. The combined effects of the processes are found to be largely additive in the tropics but nonadditive in the Southern Hemisphere (SH) high latitudes during the dynamically active season. Our results indicate that, in contrast to the original mechanism proposed in the literature, the solar-induced changes in the horizontal shortwave heating rate gradients not only in autumn/early winter but throughout the dynamically active season are important for modulating the dynamical response to changes in solar forcing. In spring, these gradients are strongly influenced by the shortwave heating anomalies at higher southern latitudes, which are closely linked to the concurrent changes in ozone. In addition, our simulations indicate differences in the winter SH dynamical responses between the experiments. We suggest a couple of potential drivers of the simulated differences, i.e. the role of enhanced zonally asymmetric ozone heating brought about by the increased solar-induced ozone levels under SMAX and/or sensitivity of the polar dynamical response to the altitude of the anomalous radiative tendencies. All in all, our results suggest that solar-induced changes in ozone, both in the tropics/mid-latitudes and the polar regions, are important for modulating the SH dynamical response to the 11-year solar cycle. In addition, the markedly nonadditive character of the SH polar vortex response simulated in austral spring highlights the need for consistent model implementation of the solar cycle forcing in both the radiative heating and photolysis schemes.

Inference of the lapse rate tropopause or the cold point from temperature profiles of finite vertical resolution entails an uncertainty of the tropopause altitude. For tropical radiosonde profiles the tropopause altitude inferred from coarse-grid profiles was found to be lower than that inferred from the original profiles. The mean (median) displacements of the lapse rate tropopause altitude when inferred from a temperature profile of 3 km vertical resolution and a Gaussian kernel are −130, −400, −730, and −590 m (−70, −230, −390, and −280 m) for Nairobi, Hilo, Munich, and Greifswald, respectively. In the case of a Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) averaging kernel, the displacement of the lapse rate tropopause altitude is −640 m. The mean (median) displacement of the cold point tropopause inferred from a temperature profile of 3 km vertical resolution (Gaussian kernels) was found to be −510, −610, −530, and −390 m (−460, −510, −370, and −280 m) for the stations mentioned above. Unsurprisingly, the tropopause altitude displacement is larger for coarser resolutions. The effect of the tropopause displacement on the water vapor saturation mixing ratio is roughly proportional to the vertical resolution. In tropical latitudes the resulting error is about 1 to 2 ppmv per vertical resolution in kilometers. The spread of the tropopause displacements within each sample of profiles seems too large as to recommend a correction scheme for tropical temperature profiles, while for midlatitudinal temperature profiles of vertical resolutions of 1 to 5 km a lapse rate of −1.3 K km−1 reproduces tropopause altitudes determined from high-resolution temperature profiles with the nominal lapse rate criterion of −2 K km−1 fairly well.

The 11-year solar cycle forcing is recognised as an important atmospheric forcing; however, there remain uncertainties in characterising the effects of solar variability on the atmosphere from observations and models. Here we present the first detailed assessment of the atmospheric response to the 11-year solar cycle in the UM-UKCA (Unified Model coupled to the United Kingdom Chemistry and Aerosol model) chemistry–climate model (CCM) using a three-member ensemble over the recent past (1966–2010). Comparison of the model simulations is made with satellite observations and reanalysis datasets. The UM-UKCA model produces a statistically significant response to the 11-year solar cycle in stratospheric temperatures, ozone and zonal winds. However, there are also differences in magnitude, spatial structure and timing of the signals compared to observational and reanalysis estimates. This could be due to deficiencies in the model performance, and so we include a critical discussion of the model limitations, and/or uncertainties in the current observational estimates of the solar cycle signals. Importantly, in contrast to many previous studies of the solar cycle impacts, we pay particular attention to the role of the chosen analysis method in UM-UKCA by comparing the model composite and a multiple linear regression (MLR) results. We show that the stratospheric solar responses diagnosed using both techniques largely agree with each other within the associated uncertainties; however, the results show that apparently different signals can be identified by the methods in the troposphere and in the tropical lower stratosphere. Lastly, we examine how internal atmospheric variability affects the detection of the 11-year solar responses in the model by comparing the results diagnosed from the three individual ensemble members (as opposed to those diagnosed from the full ensemble). We show overall agreement between the responses diagnosed for the ensemble members in the tropical and mid-latitude mid-stratosphere to lower mesosphere but larger apparent differences at Northern Hemisphere (NH) high latitudes during the dynamically active season. Our UM-UKCA results suggest the need for long data sets for confident detection of solar cycle impacts in the atmosphere, as well as for more research on possible interdependence of the solar cycle forcing with other atmospheric forcings and processes (e.g. Quasi-Biennial Oscillation, QBO; El Niño–Southern Oscillation, ENSO).

The European Climate Research Alliance (ECRA) is an association of leading European research institutions in the field of climate research (http://www.ecra-climate.eu/, last access: 6 December 2018). ECRA is a bottom-up initiative and helps to facilitate the development of climate change research, combining the capacities of national research institutions, and inducing closer ties between existing national research initiatives, projects and infrastructures. ECRA works as an open platform to bring together climate researchers, providing excellent scientific expertise for policy makers and of societal relevance. The ECRA Board consists of representatives of ECRA partners and decides on governance, scientific priorities, and organisational matters. Currently organized into four Collaborative Programmes, climate scientists share their knowledge, experience and expertise to identify the most important research requirements for the future, thus developing a foresight approach. The CPs cover the topics: (1) Arctic variability and change, (2) Sea level changes and coastal impacts, (3) Changes in the hydrological cycle and (4) High impact events. The CP activities are planned in workshops and participation is open to all interested scientists from the relevant research fields. In particular, young researchers are actively encouraged to join the network. Each CP develops its joint research priorities for shaping European research into the future. Because scientific themes are interconnected, the four Collaborative Programmes interact with each other, e.g. through the organization of common workshops or joint applications. In addition, the Collaborative Programme leads attend the Board meetings. The different formats of ECRA meetings range from scientific workshops to briefing events and side events at conferences to involve different groups of interests. This facilitates the interaction of scientists, various stakeholder groups and politicians. A biennial open ECRA General Assembly that is organised in Brussels represents an umbrella event and acts as a platform for discussion and contact with stakeholders. This event is an excellent opportunity to jointly discuss research priorities of high societal relevance.

We present the new isotope-enabled model ICON-ART-Iso. The physics package of the global ICOsahedral Nonhydrostatic (ICON) modeling framework has been extended to simulate passive moisture tracers and the stable isotopologues HDO and H218O. The extension builds on the infrastructure provided by ICON-ART, which allows for high flexibility with respect to the number of related water tracers that are simulated. The physics of isotopologue fractionation follow the model COSMOiso. We first present a detailed description of the physics of fractionation that have been implemented in the model. The model is then evaluated on a range of temporal scales by comparing with measurements of precipitation and vapor. A multi-annual simulation is compared to observations of the isotopologues in precipitation taken from the station network GNIP (Global Network for Isotopes in Precipitation). ICON-ART-Iso is able to simulate the main features of the seasonal cycles in δD and δ18O as observed at the GNIP stations. In a comparison with IASI satellite retrievals, the seasonal and daily cycles in the isotopologue content of vapor are examined for different regions in the free troposphere. On a small spatial and temporal scale, ICON-ART-Iso is used to simulate the period of two flights of the IAGOS-CARIBIC aircraft in September 2010, which sampled air in the tropopause region influenced by Hurricane Igor. The general features of this sample as well as those of all tropical data available from IAGOS-CARIBIC are captured by the model. The study demonstrates that ICON-ART-Iso is a flexible tool to analyze the water cycle of ICON. It is capable of simulating tagged water as well as the isotopologues HDO and H218O.

Atmospheric composition studies on weather and climate timescales require flexible, scalable models. The ICOsahedral Nonhydrostatic model with Aerosols and Reactive Trace gases (ICON-ART) provides such an environment. Here, we introduce the most up-to-date version of the flexible tracer framework for ICON-ART and explain its application in one numerical weather forecast and one climate related case study. We demonstrate the implementation of idealised tracers and chemistry tendencies of different complexity using the ART infrastructure. Using different ICON physics configurations for weather and climate with ART, we perform integrations on different timescales, illustrating the model's performance. First, we present a hindcast experiment for the 2002 ozone hole split with two different ozone chemistry schemes using the numerical weather prediction physics configuration. We compare the hindcast with observations and discuss the confinement of the vortex split using an idealised tracer diagnostic. Secondly, we study AMIP-type integrations using a simplified chemistry scheme in conjunction with the climate physics configuration. We use two different simulations: the interactive simulation, where modelled ozone is coupled back to the radiation scheme, and the non-interactive simulation that uses a default background climatology of ozone. Additionally, we introduce changes of water vapour by methane oxidation for the interactive simulation. We discuss the impact of stratospheric ozone and water vapour variations in the interactive and non-interactive integrations on the water vapour tape recorder, as a measure of tropical upwelling changes. Additionally we explain the seasonal evolution and latitudinal distribution of the age of air. The age of air is a measure of the strength of the meridional overturning circulation with young air in the tropical upwelling region and older air in polar winter downwelling regions. We conclude that our flexible tracer framework allows for tailor-made configurations of ICON-ART in weather and climate applications that are easy to configure and run well.

We present model simulations with the atmospheric chemistry–climate model ECHAM5/MESSy Atmospheric Chemistry (EMAC) nudged toward European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalyses for the Arctic winters 2009/2010 and 2010/2011. This study is the first to perform an extensive assessment of the performance of the EMAC model for Arctic winters as previous studies have only made limited evaluations of EMAC simulations which also were mainly focused on the Antarctic winter stratosphere. We have chosen the two extreme Arctic winters 2009/2010 and 2010/2011 to evaluate the formation of polar stratospheric clouds (PSCs) and the representation of the chemistry and dynamics of the polar winter stratosphere in EMAC. The EMAC simulations are compared to observations by the Michelson Interferometer for Passive Atmospheric Soundings (Envisat/MIPAS) and the observations from the Aura Microwave Limb Sounder (Aura/MLS). The Arctic winter 2010/2011 was one of the coldest stratospheric winters on record, leading to the strongest depletion of ozone measured in the Arctic. The Arctic winter 2009/2010 was, from the climatological perspective, one of the warmest stratospheric winters on record. However, it was distinguished by an exceptionally cold stratosphere (colder than the climatological mean) from mid-December 2009 to mid-January 2010, leading to prolonged PSC formation and existence. Significant denitrification, the removal of HNO3 from the stratosphere by sedimentation of HNO3-containing polar stratospheric cloud particles, occurred in that winter. In our comparison, we focus on PSC formation and denitrification. The comparisons between EMAC simulations and satellite observations show that model and measurements compare well for these two Arctic winters (differences for HNO3 generally within ±20 %) and thus that EMAC nudged toward ECMWF ERA-Interim reanalyses is capable of giving a realistic representation of the evolution of PSCs and associated sequestration of gas-phase HNO3 in the polar winter stratosphere. However, simulated PSC volume densities are smaller than the ones derived from Envisat/MIPAS observations by a factor of 3–7. Further, PSCs in EMAC are not simulated as high up (in altitude) as they are observed. This underestimation of PSC volume density and vertical extension of the PSCs results in an underestimation of the vertical redistribution of HNO3 due to denitrification/re-nitrification. The differences found here between model simulations and observations stipulate further improvements in the EMAC set-up for simulating PSCs.

>We analyse simulations performed for the Chemistry-Climate Model Initiative (CCMI) to estimate the return dates of the stratospheric ozone layer from depletion caused by anthropogenic stratospheric chlorine and bromine. We consider a total of 155 simulations from 20 models, including a range of sensitivity studies which examine the impact of climate change on ozone recovery. For the control simulations (unconstrained by nudging towards analysed meteorology) there is a large spread (±20 DU in the global average) in the predictions of the absolute ozone column. Therefore, the model results need to be adjusted for biases against historical data. Also, the interannual variability in the model results need to be smoothed in order to provide a reasonably narrow estimate of the range of ozone return dates. Consistent with previous studies, but here for a Representative Concentration Pathway (RCP) of 6.0, these new CCMI simulations project that global total column ozone will return to 1980 values in 2049 (with a 1σ uncertainty of 2043–2055). At Southern Hemisphere mid-latitudes column ozone is projected to return to 1980 values in 2045 (2039–2050), and at Northern Hemisphere mid-latitudes in 2032 (2020–2044). In the polar regions, the return dates are 2060 (2055–2066) in the Antarctic in October and 2034 (2025–2043) in the Arctic in March. The earlier return dates in the Northern Hemisphere reflect the larger sensitivity to dynamical changes. Our estimates of return dates are later than those presented in the 2014 Ozone Assessment by approximately 5–17 years, depending on the region, with the previous best estimates often falling outside of our uncertainty range. In the tropics only around half the models predict a return of ozone to 1980 values, around 2040, while the other half do not reach the 1980 value. All models show a negative trend in tropical total column ozone towards the end of the 21st century. The CCMI models generally agree in their simulation of the time evolution of stratospheric chlorine and bromine, which are the main drivers of ozone loss and recovery. However, there are a few outliers which show that the multi-model mean results for ozone recovery are not as tightly constrained as possible. Throughout the stratosphere the spread of ozone return dates to 1980 values between models tends to correlate with the spread of the return of inorganic chlorine to 1980 values. In the upper stratosphere, greenhouse gas-induced cooling speeds up the return by about 10–20 years. In the lower stratosphere, and for the column, there is a more direct link in the timing of the return dates of ozone and chlorine, especially for the large Antarctic depletion. Comparisons of total column ozone between the models is affected by different predictions of the evolution of tropospheric ozone within the same scenario, presumably due to differing treatment of tropospheric chemistry. Therefore, for many scenarios, clear conclusions can only be drawn for stratospheric ozone columns rather than the total column. As noted by previous studies, the timing of ozone recovery is affected by the evolution of N2O and CH4. However, quantifying the effect in the simulations analysed here is limited by the few realisations available for these experiments compared to internal model variability. The large increase in N2O given in RCP 6.0 extends the ozone return globally by ∼ 15 years relative to N2O fixed at 1960 abundances, mainly because it allows tropical column ozone to be depleted. The effect in extratropical latitudes is much smaller. The large increase in CH4 given in the RCP 8.5 scenario compared to RCP 6.0 also lengthens ozone return by ∼ 15 years, again mainly through its impact in the tropics. Overall, our estimates of ozone return dates are uncertain due to both uncertainties in future scenarios, in particular those of greenhouse gases, and uncertainties in models. The scenario uncertainty is small in the short term but increases with time, and becomes large by the end of the century. There are still some model–model differences related to well-known processes which affect ozone recovery. Efforts need to continue to ensure that models used for assessment purposes accurately represent stratospheric chemistry and the prescribed scenarios of ozone-depleting substances, and only those models are used to calculate return dates. For future assessments of single forcing or combined effects of CO2, CH4, and N2O on the stratospheric column ozone return dates, this work suggests that it is more important to have multi-member (at least three) ensembles for each scenario from every established participating model, rather than a large number of individual models.

The Stratosphere–troposphere Processes And their Role in Climate (SPARC) Quasi-Biennial Oscillation initiative (QBOi) aims to improve the fidelity of tropical stratospheric variability in general circulation and Earth system models by conducting coordinated numerical experiments and analysis. In the equatorial stratosphere, the QBO is the most conspicuous mode of variability. Five coordinated experiments have therefore been designed to (i) evaluate and compare the verisimilitude of modelled QBOs under present-day conditions, (ii) identify robustness (or alternatively the spread and uncertainty) in the simulated QBO response to commonly imposed changes in model climate forcings (e.g. a doubling of CO2 amounts), and (iii) examine model dependence of QBO predictability. This paper documents these experiments and the recommended output diagnostics. The rationale behind the experimental design and choice of diagnostics is presented. To facilitate scientific interpretation of the results in other planned QBOi studies, consistent descriptions of the models performing each experiment set are given, with those aspects particularly relevant for simulating the QBO tabulated for easy comparison.

The 2015/2016 Arctic winter was one of the coldest stratospheric winters in recent years. A stable vortex formed by early December and the early winter was exceptionally cold. Cold pool temperatures dropped below the nitric acid trihydrate (NAT) existence temperature of about 195 K, thus allowing polar stratospheric clouds (PSCs) to form. The low temperatures in the polar stratosphere persisted until early March, allowing chlorine activation and catalytic ozone destruction. Satellite observations indicate that sedimentation of PSC particles led to denitrification as well as dehydration of stratospheric layers. Model simulations of the 2015/2016 Arctic winter nudged toward European Centre for Medium-Range Weather Forecasts (ECMWF) analysis data were performed with the atmospheric chemistry–climate model ECHAM5/MESSy Atmospheric Chemistry (EMAC) for the Polar Stratosphere in a Changing Climate (POLSTRACC) campaign. POLSTRACC is a High Altitude and Long Range Research Aircraft (HALO) mission aimed at the investigation of the structure, composition and evolution of the Arctic upper troposphere and lower stratosphere (UTLS). The chemical and physical processes involved in Arctic stratospheric ozone depletion, transport and mixing processes in the UTLS at high latitudes, PSCs and cirrus clouds are investigated. In this study, an overview of the chemistry and dynamics of the 2015/2016 Arctic winter as simulated with EMAC is given. Further, chemical–dynamical processes such as denitrification, dehydration and ozone loss during the 2015/2016 Arctic winter are investigated. Comparisons to satellite observations by the Aura Microwave Limb Sounder (Aura/MLS) as well as to airborne measurements with the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) performed aboard HALO during the POLSTRACC campaign show that the EMAC simulations nudged toward ECMWF analysis generally agree well with observations. We derive a maximum polar stratospheric O3 loss of ∼ 2 ppmv or 117 DU in terms of column ozone in mid-March. The stratosphere was denitrified by about 4–8 ppbv HNO3 and dehydrated by about 0.6–1 ppmv H2O from the middle to the end of February. While ozone loss was quite strong, but not as strong as in 2010/2011, denitrification and dehydration were so far the strongest observed in the Arctic stratosphere in at least the past 10 years.

We present a recently developed emission module for the ICON (ICOsahedral Non-hydrostatic)-ART (Aerosols and Reactive Trace gases) modelling framework. The emission module processes external flux data sets and increments the tracer volume mixing ratios in the boundary layer accordingly. The performance of the emission module is illustrated with simulations of acetone, using a simplified chemical depletion mechanism based on a reaction with OH and photolysis only. In our model setup, we calculate a tropospheric acetone lifetime of 33 days, which is in good agreement with the literature. We compare our results with ground-based as well as with airborne IAGOS-CARIBIC measurements in the upper troposphere and lowermost stratosphere (UTLS) in terms of phase and amplitude of the annual cycle. In all our ICON-ART simulations the general seasonal variability is well represented but uncertainties remain concerning the magnitude of the acetone mixing ratio in the UTLS region. In addition, the module for online calculations of biogenic emissions (MEGAN2.1) is implemented in ICON-ART and can replace the offline biogenic emission data sets. In a sensitivity study we show how different parametrisations of the leaf area index (LAI) change the emission fluxes calculated by MEGAN2.1 and demonstrate the importance of an adequate treatment of the LAI within MEGAN2.1. We conclude that the emission module performs well with offline and online emission fluxes and allows the simulation of the annual cycles of emissions-dominated substances.

Measurement data from the long-term passenger aircraft project IAGOS-CARIBIC are often used to derive climatologies of trace gases in the upper troposphere and lower stratosphere (UTLS). We investigate to what extent such climatologies are representative of the true state of the atmosphere. Climatologies are considered relative to the tropopause in mid-latitudes (35 to 75° N) for trace gases with different atmospheric lifetimes. Using the chemistry–climate model EMAC, we sample the modeled trace gases along CARIBIC flight tracks. Representativeness is then assessed by comparing the CARIBIC sampled model data to the full climatological model state. Three statistical methods are applied for the investigation of representativeness: the Kolmogorov–Smirnov test and two scores based on the variability and relative differences. Two requirements for any score describing representativeness are essential: representativeness is expected to increase (i) with the number of samples and (ii) with decreasing variability of the species considered. Based on these two requirements, we investigate the suitability of the different statistical measures for investigating representativeness. The Kolmogorov–Smirnov test is very strict and does not identify any trace-gas climatology as representative – not even of long-lived trace gases. In contrast, the two scores based on either variability or relative differences show the expected behavior and thus appear applicable for investigating representativeness. For the final analysis of climatological representativeness, we use the relative difference score and calculate a representativeness uncertainty for each trace gas in percent. In order to justify the transfer of conclusions about representativeness of individual trace gases from the model to measurements, we compare the trace gas variability between model and measurements. We find that the model reaches 50–100 % of the measurement variability. The tendency of the model to underestimate the variability is caused by the relatively coarse spatial and temporal model resolution. In conclusion, we provide representativeness uncertainties for several species for tropopause-referenced climatologies. Long-lived species like CO2 have low uncertainties ( ≤ 0.4 %), while shorter-lived species like O3 have larger uncertainties (10–15 %). Finally, we translate the representativeness score into a number of flights that are necessary to achieve a certain degree of representativeness. For example, increasing the number of flights from 334 to 1000 would reduce the uncertainty in CO to a mere 1 %, while the uncertainty for shorter-lived species like NO would drop from 80 to 10 %.

The year 1980 has often been used as a benchmark for the return of Antarctic ozone to conditions assumed to be unaffected by emissions of ozone-depleting substances (ODSs), implying that anthropogenic ozone depletion in Antarctica started around 1980. Here, the extent of anthropogenically driven Antarctic ozone depletion prior to 1980 is examined using output from transient chemistry–climate model (CCM) simulations from 1960 to 2000 with prescribed changes of ozone-depleting substance concentrations in conjunction with observations. A regression model is used to attribute CCM modelled and observed changes in Antarctic total column ozone to halogen-driven chemistry prior to 1980. Wintertime Antarctic ozone is strongly affected by dynamical processes that vary in amplitude from year to year and from model to model. However, when the dynamical and chemical impacts on ozone are separated, all models consistently show a long-term, halogen-induced negative trend in Antarctic ozone from 1960 to 1980. The anthropogenically driven ozone loss from 1960 to 1980 ranges between 26.4 ± 3.4 and 49.8 ± 6.2 % of the total anthropogenic ozone depletion from 1960 to 2000. An even stronger ozone decline of 56.4 ± 6.8 % was estimated from ozone observations. This analysis of the observations and simulations from 17 CCMs clarifies that while the return of Antarctic ozone to 1980 values remains a valid milestone, achieving that milestone is not indicative of full recovery of the Antarctic ozone layer from the effects of ODSs.

Deliberate injection of aerosol particles into the stratosphere is a potential climate engineering scheme. Particles injected into the stratosphere would scatter solar radiation back to space, thereby reducing the temperature at the Earth's surface and hence the impacts of global warming. Minerals such as TiO2 or SiO2 are among the potentially suitable aerosol materials for stratospheric particle injection due to their greater light-scattering ability than stratospheric sulfuric acid particles. However, the heterogeneous reactivity of mineral particles towards trace gases important for stratospheric chemistry largely remains unknown, precluding reliable assessment of their impacts on stratospheric ozone, which is of key environmental significance. In this work we have investigated for the first time the heterogeneous hydrolysis of ClONO2 on TiO2 and SiO2 aerosol particles at room temperature and at different relative humidities (RHs), using an aerosol flow tube. The uptake coefficient, γ(ClONO2), on TiO2 was ∼ 1.2 × 10−3 at 7 % RH and remained unchanged at 33 % RH, and increased for SiO2 from ∼ 2 × 10−4 at 7 % RH to ∼ 5 × 10−4 at 35 % RH, reaching a value of ∼ 6 × 10−4 at 59 % RH. We have also examined the impacts of a hypothetical TiO2 injection on stratospheric chemistry using the UKCA (United Kingdom Chemistry and Aerosol) chemistry–climate model, in which heterogeneous hydrolysis of N2O5 and ClONO2 on TiO2 particles is considered. A TiO2 injection scenario with a solar-radiation scattering effect very similar to the eruption of Mt Pinatubo was constructed. It is found that, compared to the eruption of Mt Pinatubo, TiO2 injection causes less ClOx activation and less ozone destruction in the lowermost stratosphere, while reduced depletion of N2O5 and NOx in the middle stratosphere results in decreased ozone levels. Overall, no significant difference in the vertically integrated ozone abundances is found between TiO2 injection and the eruption of Mt Pinatubo. Future work required to further assess the impacts of TiO2 injection on stratospheric chemistry is also discussed.

Future trends in Arctic springtime total column ozone, and its chemical and dynamical drivers, are assessed using a seven-member ensemble from the Met Office Unified Model with United Kingdom Chemistry and Aerosols (UM-UKCA) simulating the period 1960–2100. The Arctic mean March total column ozone increases throughout the 21st century at a rate of ∼ 11.5 DU decade−1, and is projected to return to the 1980 level in the late 2030s. However, the integrations show that even past 2060 springtime Arctic ozone can episodically drop by ∼ 50–100 DU below the corresponding long-term ensemble mean for that period, reaching values characteristic of the near-present-day average level. Consistent with the global decline in inorganic chlorine (Cly) over the century, the estimated mean halogen-induced chemical ozone loss in the Arctic lower atmosphere in spring decreases by around a factor of 2 between the periods 2001–2020 and 2061–2080. However, in the presence of a cold and strong polar vortex, elevated halogen-induced ozone losses well above the corresponding long-term mean continue to occur in the simulations into the second part of the century. The ensemble shows a significant cooling trend in the Arctic winter mid- and upper stratosphere, but there is less confidence in the projected temperature trends in the lower stratosphere (100–50 hPa). This is partly due to an increase in downwelling over the Arctic polar cap in winter, which increases transport of ozone into the polar region as well as drives adiabatic warming that partly offsets the radiatively driven stratospheric cooling. However, individual winters characterised by significantly suppressed downwelling, reduced transport and anomalously low temperatures continue to occur in the future. We conclude that, despite the projected long-term recovery of Arctic ozone, the large interannual dynamical variability is expected to continue in the future, thereby facilitating episodic reductions in springtime ozone columns. Whilst our results suggest that the relative role of dynamical processes for determining Arctic springtime ozone will increase in the future, halogen chemistry will remain a smaller but non-negligible contributor for many decades to come.

During boreal summer the upper troposphere/lower stratosphere (UTLS) in the Northern Hemisphere shows a distinct maximum in water vapour (H2O) mixing ratios and a coincident minimum in ozone (O3) mixing ratios, both confined within the Asian monsoon anticyclone (AMA). This well-known feature has been related to transport processes emerging above the convective systems during the Asian summer monsoon (ASM), further modified by the dynamics of the AMA. We compare the ability of chemistry–climate models (CCMs) to reproduce the climatological characteristics and variability of H2O, O3, and temperature in the UTLS during the boreal summer with MIPAS satellite observations and ERA-Interim reanalyses. By using a multiple linear regression model the main driving factors, the strength of the ASM, the quasi-biennial oscillation (QBO), and the El Niño–Southern Oscillation (ENSO), are separated. The regression patterns related to ENSO show a coherent, zonally asymmetric signal for temperatures and H2O mixing ratios for ERA-Interim and the CCMs and suggest a weakening of the ASM during ENSO warm events. The QBO modulation of the lower-stratospheric temperature near the Equator is well represented as a zonally symmetric pattern in the CCMs. Changes in H2O and O3 mixing ratios are consistent with the QBO-induced temperature and circulation anomalies but less zonally symmetric than the temperature pattern. Regarding the ASM, the results of the regression analysis show for ERA-Interim and the CCMs enhanced H2O and reduced O3 mixing ratios within the AMA for stronger ASM seasons. The CCM results can further confirm earlier studies which emphasize the importance of the Tibetan Plateau/southern slope of the Himalayas as the main source region for H2O in the AMA. The results suggest that H2O is transported towards higher latitudes at the north-eastern edge of the AMA rather than towards low equatorial latitudes to be fed into the tropical pipe.

We introduce a diagnostic tool to assess a climatological framework of the optimal propagation conditions for stationary planetary waves. Analyzing 50 winters using NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric Research) reanalysis data we derive probability density functions (PDFs) of positive vertical wave number as a function of zonal and meridional wave numbers. We contrast this quantity with classical climatological means of the vertical wave number. Introducing a membership value function (MVF) based on fuzzy logic, we objectively generate a modified set of PDFs (mPDFs) and demonstrate their superior performance compared to the climatological mean of vertical wave number and the original PDFs. We argue that mPDFs allow an even better understanding of how background conditions impact wave propagation in a climatological sense. As expected, probabilities are decreasing with increasing zonal wave numbers. In addition we discuss the meridional wave number dependency of the PDFs which is usually neglected, highlighting the contribution of meridional wave numbers 2 and 3 in the stratosphere. We also describe how mPDFs change in response to strong vortex regime (SVR) and weak vortex regime (WVR) conditions, with increased probabilities of the wave propagation during WVR than SVR in the stratosphere. We conclude that the mPDFs are a convenient way to summarize climatological information about planetary wave propagation in reanalysis and climate model data.

A stratosphere-resolving configuration of the Met Office's Unified Model (UM) with the United Kingdom Chemistry and Aerosols (UKCA) scheme is used to investigate the atmospheric response to changes in (a) greenhouse gases and climate, (b) ozone-depleting substances (ODSs) and (c) non-methane ozone precursor emissions. A suite of time-slice experiments show the separate, as well as pairwise, impacts of these perturbations between the years 2000 and 2100. Sensitivity to uncertainties in future greenhouse gases and aerosols is explored through the use of the Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios. The results highlight an important role for the stratosphere in determining the annual mean tropospheric ozone response, primarily through stratosphere–troposphere exchange (STE) of ozone. Under both climate change and reductions in ODSs, increases in STE offset decreases in net chemical production and act to increase the tropospheric ozone burden. This opposes the effects of projected decreases in ozone precursors through measures to improve air quality, which act to reduce the ozone burden. The global tropospheric lifetime of ozone (τO3) does not change significantly under climate change at RCP4.5, but it decreases at RCP8.5. This opposes the increases in τO3 simulated under reductions in ODSs and ozone precursor emissions. The additivity of the changes in ozone is examined by comparing the sum of the responses in the single-forcing experiments to those from equivalent combined-forcing experiments. Whilst the ozone responses to most forcing combinations are found to be approximately additive, non-additive changes are found in both the stratosphere and troposphere when a large climate forcing (RCP8.5) is combined with the effects of ODSs.

Various forms of geoengineering have been proposed to counter anthropogenic climate change. Methods which aim to modify the Earth's energy balance by reducing insolation are often subsumed under the term solar radiation management (SRM). Here, we present results of a standard SRM modelling experiment in which the incoming solar irradiance is reduced to offset the global mean warming induced by a quadrupling of atmospheric carbon dioxide. For the first time in an atmosphere–ocean coupled climate model, we include atmospheric composition feedbacks for this experiment. While the SRM scheme considered here could offset greenhouse gas induced global mean surface warming, it leads to important changes in atmospheric composition. We find large stratospheric ozone increases that induce significant reductions in surface UV-B irradiance, which would have implications for vitamin D production. In addition, the higher stratospheric ozone levels lead to decreased ozone photolysis in the troposphere. In combination with lower atmospheric specific humidity under SRM, this results in overall surface ozone concentration increases in the idealized G1 experiment. Both UV-B and surface ozone changes are important for human health. We therefore highlight that both stratospheric and tropospheric ozone changes must be considered in the assessment of any SRM scheme, due to their important roles in regulating UV exposure and air quality.

More than a decade ago it was suggested that a cooling of stratospheric temperatures by 1 K or an increase of 1 ppmv of stratospheric water vapour could promote denitrification, the permanent removal of nitrogen species from the stratosphere by solid polar stratospheric cloud (PSC) particles. In fact, during the two Arctic winters 2009/10 and 2010/11 the strongest denitrification in the recent decade was observed. Sensitivity studies along air parcel trajectories are performed to test how a future stratospheric water vapour (H2O) increase of 1 ppmv or a temperature decrease of 1 K would affect PSC formation. We perform our study based on measurements made during the Arctic winter 2010/11. Air parcel trajectories were calculated 6 days backward in time based on PSCs detected by CALIPSO (Cloud Aerosol Lidar and Infrared Pathfinder satellite observations). The sensitivity study was performed on single trajectories as well as on a trajectory ensemble. The sensitivity study shows a clear prolongation of the potential for PSC formation and PSC existence when the temperature in the stratosphere is decreased by 1 K and water vapour is increased by 1 ppmv. Based on 15 years of satellite measurements (2000–2014) from UARS/HALOE, Envisat/MIPAS, Odin/SMR, Aura/MLS, Envisat/SCIAMACHY and SCISAT/ACE-FTS it is further investigated if there is a decrease in temperature and/or increase of water vapour (H2O) observed in the polar regions similar to that observed at midlatitudes and in the tropics. Performing linear regression analyses we derive from the Envisat/MIPAS (2002–2012) and Aura/MLS (2004–2014) observations predominantly positive changes in the potential temperature range 350 to 1000 K. The linear changes in water vapour derived from Envisat/MIPAS observations are largely insignificant, while those from Aura/MLS are mostly significant. For the temperature neither of the two instruments indicate any significant changes. Given the strong inter-annual variation observed in water vapour and particular temperature the severe denitrification observed in 2010/11 cannot be directly related to any changes in water vapour and temperature since the millennium. However, the observations indicate a clear correlation between cold winters and enhanced water vapour mixing ratios. This indicates a connection between dynamical and radiative processes that govern water vapour and temperature in the Arctic lower stratosphere.

An important source of polar stratospheric clouds (PSCs), which play a crucial role in controlling polar stratospheric ozone depletion, is from the temperature fluctuations induced by mountain waves. However, this formation mechanism is usually missing in chemistry–climate models because these temperature fluctuations are neither resolved nor parameterised. Here, we investigate the representation of stratospheric mountain-wave-induced temperature fluctuations by the UK Met Office Unified Model (UM) at climate scale and mesoscale against Atmospheric Infrared Sounder satellite observations for three case studies over the Antarctic Peninsula. At a high horizontal resolution (4 km) the regional mesoscale configuration of the UM correctly simulates the magnitude, timing, and location of the measured temperature fluctuations. By comparison, at a low horizontal resolution (2.5° × 3.75°) the global climate configuration fails to resolve such disturbances. However, it is demonstrated that the temperature fluctuations computed by a mountain wave parameterisation scheme inserted into the climate configuration (which computes the temperature fluctuations due to unresolved mountain waves) are in relatively good agreement with the mesoscale configuration responses for two of the three case studies. The parameterisation was used to include the simulation of mountain-wave-induced PSCs in the global chemistry–climate configuration of the UM. A subsequent sensitivity study demonstrated that regional PSCs increased by up to 50% during July over the Antarctic Peninsula following the inclusion of the local mountain-wave-induced cooling phase.

The impact of polar stratospheric ozone loss resulting from chlorine activation on polar stratospheric clouds is examined using a pair of model integrations run with the fully coupled chemistry climate model UM-UKCA. Suppressing chlorine activation through heterogeneous reactions is found to produce modelled ozone differences consistent with observed ozone differences between the present and pre-ozone hole period. Statistically significant high-latitude Southern Hemisphere (SH) ozone loss begins in August and peaks in October–November, with > 75% of ozone destroyed at 50 hPa. Associated with this ozone destruction is a > 12 K decrease of the lower polar stratospheric temperatures and an increase of > 6 K in the upper stratosphere. The heating components of this temperature change are diagnosed and it is found that the temperature dipole is the result of decreased short-wave heating in the lower stratosphere and increased dynamical heating in the upper stratosphere. The cooling of the polar lower stratosphere leads, through thermal wind balance, to an acceleration of the polar vortex and delays its breakdown by ~ 2 weeks. A link between lower stratospheric zonal wind speed, the vertical component of the Eliassen–Palm (EP) flux, Fz and the residual mean vertical circulation, w*, is identified. In November and December, increased westerly winds and a delay in the breakup of the polar vortex lead to increases in Fz, indicating increased wave activity entering the stratosphere and propagating to higher altitudes. The resulting increase in wave breaking, diagnosed by decreases to the EP flux divergence, drives enhanced downwelling over the polar cap. Many of the stratospheric signals modelled in this study propagate down to the troposphere, and lead to significant surface changes in December.

Naturally produced very short-lived substances (VSLS) account for almost a quarter of the current stratospheric inorganic bromine, Bry. Following VSLS oxidation, bromine radicals (Br and BrO) can catalytically destroy ozone. The extent to which possible increases in surface emissions or transport of these VSLS bromocarbons to the stratosphere could counteract the effect of halogen reductions under the Montreal Protocol is an important policy question. Here, by using a chemistry–climate model, UM-UKCA, we investigate the impact of a hypothetical doubling (an increase of 5 ppt Bry) of VSLS bromocarbons on ozone and how the resulting ozone changes depend on the background concentrations of chlorine and bromine. Our model experiments indicate that for the 5 ppt increase in Bry from VSLS, the ozone decrease in the lowermost stratosphere of the Southern Hemisphere (SH) may reach up to 10% in the annual mean; the ozone decrease in the Northern Hemisphere (NH) is smaller (4–6%). The largest impact on the ozone column is found in the Antarctic spring. There is a significantly larger ozone decrease following the doubling of the VSLS burden under a high stratospheric chlorine background than under a low chlorine background, indicating the importance of the inter-halogen reactions. For example, the decline in the high-latitude, lower-stratospheric ozone concentration as a function of Bry is higher by about 30–40% when stratospheric Cly is ~ 3 ppb (present day), compared with Cly of ~ 0.8 ppb (a pre-industrial or projected future situation). Bromine will play an important role in the future ozone layer. However, even if bromine levels from natural VSLS were to increase significantly later this century, changes in the concentration of ozone will likely be dominated by the decrease in anthropogenic chlorine. Our calculation suggests that for a 5 ppt increase in Bry from VSLS, the Antarctic ozone hole recovery date could be delayed by approximately 6–8 years, depending on Cly levels.

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.

Lightning is one of the major natural sources of NOx in the atmosphere. A suite of time slice experiments using a stratosphere-resolving configuration of the Unified Model (UM), containing the United Kingdom Chemistry and Aerosols sub-model (UKCA), has been performed to investigate the impact of climate change on emissions of NOx from lightning (LNOx) and to highlight its critical impacts on photochemical ozone production and the oxidising capacity of the troposphere. Two Representative Concentration Pathway (RCP) scenarios (RCP4.5 and RCP8.5) are explored. LNOx is simulated to increase in a year-2100 climate by 33% (RCP4.5) and 78% (RCP8.5), primarily as a result of increases in the depth of convection. The total tropospheric chemical odd oxygen production (P(Ox)) increases linearly with increases in total LNOx and consequently, tropospheric ozone burdens of 29 ± 4 Tg(O3) (RCP4.5) and 46 ± 4 Tg(O3) (RCP8.5) are calculated here. By prescribing a uniform surface boundary concentration for methane in these simulations, methane-driven feedbacks are essentially neglected. A simple estimate of the contribution of the feedback reduces the increase in ozone burden to 24 and 33 Tg(O3), respectively. We thus show that, through changes in LNOx, the effects of climate change counteract the simulated mitigation of the ozone burden, which results from reductions in ozone precursor emissions as part of air quality controls projected in the RCP scenarios. Without the driver of increased LNOx, our simulations suggest that the net effect of climate change would be to lower free tropospheric ozone. In addition, we identify large climate-change-induced enhancements in the concentration of the hydroxyl radical (OH) in the tropical upper troposphere (UT), particularly over the Maritime Continent, primarily as a consequence of greater LNOx. The OH enhancement in the tropics increases oxidation of both methane (with feedbacks onto chemistry and climate) and very short-lived substances (VSLS) (with implications for stratospheric ozone depletion). We emphasise that it is important to improve our understanding of LNOx in order to gain confidence in model projections of composition change under future climate.

Injection of aerosol particles (or their precursors) into the stratosphere to scatter solar radiation back into space has been suggested as a solar-radiation management scheme for the mitigation of global warming. TiO2 has recently been highlighted as a possible candidate particle because of its high refractive index, but its impact on stratospheric chemistry via heterogeneous reactions is as yet unknown. In this work the heterogeneous reaction of airborne sub-micrometre TiO2 particles with N2O5 has been investigated for the first time, at room temperature and different relative humidities (RH), using an atmospheric pressure aerosol flow tube. The uptake coefficient of N2O5 onto TiO2, γ(N2O5), was determined to be ~1.0 × 10−3 at low RH, increasing to ~3 × 10−3 at 60% RH. The uptake of N2O5 onto TiO2 is then included in the UKCA chemistry–climate model to assess the impact of this reaction on stratospheric chemistry. While the impact of TiO2 on the scattering of solar radiation is chosen to be similar to the aerosol from the Mt Pinatubo eruption, the impact of TiO2 injection on stratospheric N2O5 is much smaller.

In this paper, we present a description of the tropospheric chemistry component of the UK Chemistry and Aerosols (UKCA) model which has been coupled to the Met Office Hadley Centre's HadGEM family of climate models. We assess the model's transport and scavenging processes, in particular focussing on convective transport, boundary layer mixing, wet scavenging and inter-hemispheric exchange. Simulations with UKCA of the short-lived radon tracer suggest that modelled distributions are comparable to those of other models and the comparison with observations indicate that apart from a few locations, boundary layer mixing and convective transport are effective in the model as a means of vertically redistributing surface emissions of radon. Comparisons of modelled lead tracer concentrations with observations suggest that UKCA captures surface concentrations in both hemispheres very well, although there is a tendency to underestimate the observed geographical and interannual variability in the Northern Hemisphere. In particular, UKCA replicates the shape and absolute concentrations of observed lead profiles, a key test in the evaluation of a model's wet scavenging scheme. The timescale for inter-hemispheric transport, calculated in the model using a simple krypton tracer experiment, does appear to be long relative to other models and could indicate deficiencies in tropical deep convection and/or insufficient boundary layer mixing. We also describe the main components of the tropospheric chemistry and evaluate it against observations and other tropospheric chemistry models. In particular, from a climate forcing perspective, present-day observed surface methane concentrations and tropospheric ozone concentrations are reproduced very well by the model, thereby making it suitable for long centennial integrations as well as studies of biogeochemical feedbacks. Results from both historical and future simulations with UKCA tropospheric chemistry are presented. Future projections of tropospheric ozone vary with the Representative Concentration Pathway (RCP). In RCP2.6, for example, tropospheric ozone increases up to 2010 and then declines by 13% of its year-2000 global mean by the end of the century. In RCP8.5, tropospheric ozone continues to rise steadily throughout the 21st century, with methane being the main driving factor. Finally, we highlight aspects of the UKCA model which are undergoing and/or have undergone recent developments and are suitable for inclusion in a next-generation Earth System Model.

We report results from two pairs of chemistry-climate model simulations using the same climate model but different chemical perturbations. In each pair of experiments an ozone change was triggered by a simple change in the chemistry. One pair of model experiments looked at the impact of polar stratospheric clouds (PSCs) and the other pair at the impact of short-lived halogenated species on composition and circulation. The model response is complex with both positive and negative changes in ozone concentration, depending on location. These changes result from coupling between composition, temperature and circulation. Even though the causes of the modelled ozone changes are different, the high latitude Southern Hemisphere response in the lower stratosphere is similar. In both pairs of experiments the high-latitude circulation changes, as evidenced by N2O differences, are suggesting a slightly longer-lasting/stronger stratospheric descent in runs with higher ozone destruction (a manifestation of a seasonal shift in the circulation). We contrast the idealised model behaviour with interannual variability in ozone and N2O as observed by the MIPAS instrument on ENVISAT, highlighting similarities of the modelled climate equilibrium changes to the year 2006–2007 in observations. We conclude that the climate system can respond quite sensitively in its seasonal evolution to small chemical perturbations, that circulation adjustments seen in the model can occur in reality, and that coupled chemistry-climate models allow a better assessment of future ozone and climate change than recent CMIP-type models with prescribed ozone fields.

Atmospheric chemistry is driven by photolytic reactions, making their modelling a crucial component of atmospheric models. We describe the implementation and validation of Fast-JX, a state of the art model of interactive photolysis, into the MetUM chemistry-climate model. This allows for interactive photolysis rates to be calculated in the troposphere and augments the calculation of the rates in the stratosphere by accounting for clouds and aerosols in addition to ozone. In order to demonstrate the effectiveness of this new photolysis scheme we employ new methods of validating the model, including techniques for sampling the model to compare to flight track and satellite data.

We introduce a methodology to visualise rapid vertical and zonal tropical transport pathways. Using prescribed sea-surface temperatures in four monthly model integrations for 2005, we characterise preferred transport routes from the troposphere to the stratosphere in a high resolution climate model. Most efficient transport is modelled over the Maritime Continent (MC) in November and February, i.e., boreal winter. In these months, the ascending branch of the Walker Circulation over the MC is formed in conjunction with strong deep convection, allowing fast transport into the stratosphere. In the model the upper tropospheric zonal winds associated with the Walker Circulation are also greatest in these months in agreement with ERA-Interim reanalysis data. We conclude that the Walker circulation plays an important role in the seasonality of fast tropical transport from the lower and middle troposphere to the upper troposphere and so impacts at the same time the potential supply of surface emissions to the tropical tropopause layer (TTL) and subsequently to the stratosphere.

Monthly averaged surface erythemal solar irradiance (UV-Ery) for local noon from 1960 to 2100 has been derived using radiative transfer calculations and projections of ozone, temperature and cloud change from 14 chemistry climate models (CCM), as part of the CCMVal-2 activity of SPARC. Our calculations show the influence of ozone depletion and recovery on erythemal irradiance. In addition, we investigate UV-Ery changes caused by climate change due to increasing greenhouse gas concentrations. The latter include effects of both stratospheric ozone and cloud changes. The derived estimates provide a global picture of the likely changes in erythemal irradiance during the 21st century. Uncertainties arise from the assumed scenarios, different parameterizations – particularly of cloud effects on UV-Ery – and the spread in the CCM projections. The calculations suggest that relative to 1980, annually mean UV-Ery in the 2090s will be on average ~12 % lower at high latitudes in both hemispheres, ~3 % lower at mid latitudes, and marginally higher (~1 %) in the tropics. The largest reduction (~16 %) is projected for Antarctica in October. Cloud effects are responsible for 2–3 % of the reduction in UV-Ery at high latitudes, but they slightly moderate it at mid-latitudes (~1 %). The year of return of erythemal irradiance to values of certain milestones (1965 and 1980) depends largely on the return of column ozone to the corresponding levels and is associated with large uncertainties mainly due to the spread of the model projections. The inclusion of cloud effects in the calculations has only a small effect of the return years. At mid and high latitudes, changes in clouds and stratospheric ozone transport by global circulation changes due to greenhouse gases will sustain the erythemal irradiance at levels below those in 1965, despite the removal of ozone depleting substances. At northern high latitudes (60°–90°), the projected decreases in cloud transmittance towards the end of the 21st century will reduce the yearly average surface erythemal irradiance by ~5 % with respect to the 1960s.

The correlation between measured tropospheric ozone (O3) and carbon monoxide (CO) has been used extensively in tropospheric chemistry studies to explore the photochemical characteristics of different regions and to evaluate the ability of models to capture these characteristics. Here, we present the first study that uses multi-year, global, vertically resolved, simultaneous and collocated O3 and CO satellite (Tropospheric Emission Spectrometer) measurements, to determine this correlation in the middle/lower free troposphere for two different seasons, and to evaluate two chemistry-climate models. We find results that are fairly robust across different years, altitudes and timescales considered, which indicates that the correlation maps presented here could be used in future model evaluations. The highest positive correlations (around 0.8) are found in the northern Pacific during summer, which is a common feature in the observations and the G-PUCCINI model. We make quantitative comparisons between the models using a single-figure metric (C), which we define as the correlation coefficient between the modeled and the observed O3-CO correlations for different regions of the globe. On a global scale, the G-PUCCINI model shows a good performance in the summer (C=0.71) and a satisfactory performance in the winter (C=0.52). It captures midlatitude features very well, especially in the summer, whereas the performance in regions like South America or Central Africa is weaker. The UKCA model (C=0.46/0.15 for July–August/December–January on a global scale) performs better in certain regions, such as the tropics in winter, and it captures some of the broad characteristics of summer extratropical correlations, but it systematically underestimates the O3-CO correlations over much of the globe. It is noteworthy that the correlations look very different in the two models, even though the ozone distributions are similar. This demonstrates that this technique provides a powerful global constraint for understanding modeled tropospheric chemical processes. We investigated the sources of the correlations by performing a series of sensitivity experiments. In these, the sign of the correlation is, in most cases, insensitive to removing different individual emissions, but its magnitude changes downwind of emission regions when applying such perturbations. Interestingly, we find that the O3-CO correlation does not solely reflect the strength of O3 photochemical production, as often assumed by earlier studies, but is more complicated and may reflect a mixture of different processes such as transport.

Three recently-completed sets of simulations of multiple chemistry-climate models with greenhouse gases only, with all anthropogenic forcings, and with anthropogenic and natural forcings, allow the causes of observed stratospheric changes to be quantitatively assessed using detection and attribution techniques. The total column ozone response to halogenated ozone depleting substances and to natural forcings is detectable in observations, but the total column ozone response to greenhouse gas changes is not separately detectable. In the middle and upper stratosphere, simulated and observed SBUV/SAGE ozone changes are broadly consistent, and separate anthropogenic and natural responses are detectable in observations. The influence of ozone depleting substances and natural forcings can also be detected separately in observed lower stratospheric temperature, and the magnitudes of the simulated and observed responses to these forcings and to greenhouse gas changes are found to be consistent. In the mid and upper stratosphere the simulated natural and combined anthropogenic responses are detectable and consistent with observations, but the influences of greenhouse gases and ozone-depleting substances could not be separately detected in our analysis.

The UK Met Office's Unified Model is used at a climate resolution (N216, ~0.83°×~0.56°, ~60 km) to assess the impact of deep tropical convection on the structure of the tropical tropopause layer (TTL). We focus on the potential for rapid transport of short-lived ozone depleting species to the stratosphere by rapid convective uplift. The modelled horizontal structure of organised convection is shown to match closely with signatures found in the OLR satellite data. In the model, deep convective elevators rapidly lift air from 4–5 km up to 12–14 km. The influx of tropospheric air entering the TTL (11–12 km) is similar for all tropical regions with most convection stopping below ~14 km. The tropical tropopause is coldest and driest between November and February, coinciding with the greatest upwelling over the tropical warm pool. As this deep convection is co-located with bromine-rich biogenic coastal emissions, this period and location could potentially be the preferential gateway for stratospheric bromine.

Projections of stratospheric ozone from a suite of chemistry-climate models (CCMs) have been analyzed. In addition to a reference simulation where anthropogenic halogenated ozone depleting substances (ODSs) and greenhouse gases (GHGs) vary with time, sensitivity simulations with either ODS or GHG concentrations fixed at 1960 levels were performed to disaggregate the drivers of projected ozone changes. These simulations were also used to assess the two distinct milestones of ozone returning to historical values (ozone return dates) and ozone no longer being influenced by ODSs (full ozone recovery). The date of ozone returning to historical values does not indicate complete recovery from ODSs in most cases, because GHG-induced changes accelerate or decelerate ozone changes in many regions. In the upper stratosphere where CO2-induced stratospheric cooling increases ozone, full ozone recovery is projected to not likely have occurred by 2100 even though ozone returns to its 1980 or even 1960 levels well before (~2025 and 2040, respectively). In contrast, in the tropical lower stratosphere ozone decreases continuously from 1960 to 2100 due to projected increases in tropical upwelling, while by around 2040 it is already very likely that full recovery from the effects of ODSs has occurred, although ODS concentrations are still elevated by this date. In the midlatitude lower stratosphere the evolution differs from that in the tropics, and rather than a steady decrease in ozone, first a decrease in ozone is simulated from 1960 to 2000, which is then followed by a steady increase through the 21st century. Ozone in the midlatitude lower stratosphere returns to 1980 levels by ~2045 in the Northern Hemisphere (NH) and by ~2055 in the Southern Hemisphere (SH), and full ozone recovery is likely reached by 2100 in both hemispheres. Overall, in all regions except the tropical lower stratosphere, full ozone recovery from ODSs occurs significantly later than the return of total column ozone to its 1980 level. The latest return of total column ozone is projected to occur over Antarctica (~2045–2060) whereas it is not likely that full ozone recovery is reached by the end of the 21st century in this region. Arctic total column ozone is projected to return to 1980 levels well before polar stratospheric halogen loading does so (~2025–2030 for total column ozone, cf. 2050–2070 for Cly+60×Bry) and it is likely that full recovery of total column ozone from the effects of ODSs has occurred by ~2035. In contrast to the Antarctic, by 2100 Arctic total column ozone is projected to be above 1960 levels, but not in the fixed GHG simulation, indicating that climate change plays a significant role.

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an "ensemble of opportunity" of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21st century, up-to and after the time when ozone concentrations return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

In the 1990s the rates of increase of greenhouse gas concentrations, most notably of methane, were observed to change, for reasons that have yet to be fully determined. This period included the eruption of Mt. Pinatubo and an El Niño warm event, both of which affect biogeochemical processes, by changes in temperature, precipitation and radiation. We examine the impact of these changes in climate on global isoprene emissions and the effect these climate dependent emissions have on the hydroxy radical, OH, the dominant sink for methane. We model a reduction of isoprene emissions in the early 1990s, with a maximum decrease of 40 Tg(C)/yr in late 1992 and early 1993, a change of 9%. This reduction is caused by the cooler, drier conditions following the eruption of Mt. Pinatubo. Isoprene emissions are reduced both directly, by changes in temperature and a soil moisture dependent suppression factor, and indirectly, through reductions in the total biomass. The reduction in isoprene emissions causes increases of tropospheric OH which lead to an increased sink for methane of up to 5 Tg(CH4)/year, comparable to estimated source changes over the time period studied. There remain many uncertainties in the emission and oxidation of isoprene which may affect the exact size of this effect, but its magnitude is large enough that it should remain important.

We have run a chemistry transport model (CTM) to systematically examine the drivers of interannual variability of tropospheric composition during 1996–2000. This period was characterised by anomalous meteorological conditions associated with the strong El Niño of 1997–1998 and intense wildfires, which produced a large amount of pollution. On a global scale, changing meteorology (winds, temperatures, humidity and clouds) is found to be the most important factor driving interannual variability of NO2 and ozone on the timescales considered. Changes in stratosphere-troposphere exchange, which are largely driven by meteorological variability, are found to play a particularly important role in driving ozone changes. The strong influence of emissions on NO2 and ozone interannual variability is largely confined to areas where intense biomass burning events occur. For CO, interannual variability is almost solely driven by emission changes, while for OH meteorology dominates, with the radiative influence of clouds being a very strong contributor. Through a simple attribution analysis for 1996–2000 we conclude that changing cloudiness drives 25% of the interannual variability of OH over Europe by affecting shortwave radiation. Over Indonesia this figure is as high as 71%. Changes in cloudiness contribute a small but non-negligible amount (up to 6%) to the interannual variability of ozone over Europe and Indonesia. This suggests that future assessments of trends in tropospheric oxidizing capacity should account for interannual variability in cloudiness, a factor neglected in many previous studies.