Search

We develop a novel technique for characterising the latitude (ϕ‾), tilt (α) and intensity (Umean) of the North Atlantic eddy-driven jet using a feature identification method and two-dimensional moment analysis. Applying this technique to the ERA5 reanalysis, the distribution of the daily winter ϕ‾ is unimodal, which is in contrast to the trimodal distribution of the daily jet latitude index (JLI) described by Woollings et al. (2010). We show that our method exhibits a higher persistence than the JLI, casting doubt on the previous interpretations of the trimodal distribution as evidence for regime behaviour of the North Atlantic jet. It also explicitly and straightforwardly handles days where the jet is split. Although climatologically α is positive, indicating a tilt from south to east, around a fifth of winter days show negative α. When plotted as a function of the North Atlantic Oscillation and East Atlantic pattern indices, there is a higher fraction of explained variance in the daily ϕ‾ within each quadrant of the phase space than is found for JLI, supporting the conclusion that ϕ‾ has smoother variations than JLI and has a closer relationship with these indices. Our method is simple, requiring only the daily 850 hPa zonal wind as input data, and diagnoses the jet in a more informative and robust way than other methods using low-level wind fields.

The underestimation in multidecadal variability in the wintertime North Atlantic Oscillation (NAO) by global climate models remains poorly understood. Understanding the origins of this weak NAO variability is important for making model projections more reliable. Past studies have linked the weak multidecadal NAO variability in models to an underestimated atmospheric response to the Atlantic Multidecadal Variability (AMV). We investigate historical simulations from Coupled Model Intercomparison Project Phase 6 (CMIP6) large-ensemble models and find that most of the models do not reproduce observed multidecadal NAO variability, as found in previous generations of climate models. We explore statistical relationships with physical drivers that may contribute to inter-model spread in NAO variability. There is a significant anticorrelation across models between the AMV–NAO coupling parameter and multidecadal NAO variability over the full historical period (r=-0.55, p<0.05). However, this relationship is relatively weak and becomes obscured when using a common period (1900–2010) and de-trending the data in a consistent way, with observations to enable a model–data comparison. This suggests that the representation of NAO–AMV coupling contributes to a modest proportion of inter-model spread in multidecadal NAO variability, although the importance of this process for model spread could be underestimated, given evidence of a systematically poor representation of the coupling in the models. We find a significant inter-model correlation between multidecadal NAO variability and multidecadal stratospheric polar vortex variability and a stratosphere–troposphere coupling parameter, which quantifies the relationship between stratospheric winds and the NAO. The models with the lowest NAO variance are associated with weaker polar vortex variability and a weaker stratosphere–troposphere coupling parameter. The two stratospheric indices are uncorrelated across models and together give a pooled R2 with an NAO variability of 0.7, which is larger than the fraction of inter-model spread related to the AMV (R2=0.3). The identification of this relationship suggests that modelled spread in multidecadal NAO variability has the potential to be reduced by improved knowledge of observed multidecadal stratospheric variability; however, observational records are currently too short to provide a robust constraint on these indices.

It has been proposed that externally forced trends in the Aleutian Low can induce a basin-wide Pacific sea surface temperature (SST) response that projects onto the pattern of the Pacific Decadal Oscillation (PDO). To investigate this hypothesis, we apply local atmospheric nudging in an intermediate-complexity climate model to isolate the effects of an intensified winter Aleutian Low sustained over several decades. An intensification of the Aleutian Low produces a basin-wide SST response with a similar pattern to the model's internally generated PDO. The amplitude of the SST response in the North Pacific is comparable to the PDO, but in the tropics and southern subtropics the anomalies induced by the imposed Aleutian Low anomaly are a factor of 3 weaker than for the internally generated PDO. The tropical Pacific warming peaks in boreal spring, though anomalies persist year-round. A heat budget analysis shows the northern subtropical Pacific SST response is predominantly driven by anomalous surface turbulent heat fluxes in boreal winter, while in the equatorial Pacific the response is mainly due to meridional heat advection in boreal spring. The propagation of anomalies from the extratropics to the tropics can be explained by the seasonal footprinting mechanism, involving the wind–evaporation–SST feedback. The results show that low-frequency variability and trends in the Aleutian Low could contribute to basin-wide anomalous Pacific SST, but the magnitude of the effect in the tropical Pacific, even for the extreme Aleutian Low forcing applied here, is small. Therefore, external forcing of the Aleutian Low is unlikely to account for observed decadal SST trends in the tropical Pacific in the late 20th and early 21st centuries.

The regional climate impacts of hypothetical future emissions scenarios can be estimated by combining Earth system model simulations with a linear pattern scaling model such as MESMER (Modular Earth System Model Emulator with spatially Resolved output), which uses estimated patterns of the local response per degree of global temperature change. Here we use the mean trend component of MESMER to emulate the regional pattern of the surface temperature response based on historical single-forcer and future Shared Socioeconomic Pathway (SSP) CMIP6 (Coupled Model Intercomparison Project Phase 6) simulations. Errors in the emulations for selected target scenarios (SSP1–1.9, SSP1–2.6, SSP2–4.5, SSP3–7.0, and SSP5–8.5) are decomposed into two components, namely (1) the differences in scaling patterns between scenarios as a consequence of varying combinations of external forcings and (2) the intrinsic time series differences between the local and global responses in the target scenario. The time series error is relatively small for high-emissions scenarios, contributing around 20 % of the total error, but is similar in magnitude to the pattern error for lower-emissions scenarios. This irreducible time series error limits the efficacy of linear pattern scaling for emulating strong mitigation pathways and reduces the dependence on the predictor pattern used. The results help guide the choice of predictor scenarios for simple climate models and where to target for the introduction of other dependent variables beyond global surface temperature into pattern scaling models.

Concentrations of water vapour entering the tropical lower stratosphere are primarily determined by conditions that air parcels encounter as they are transported through the tropical tropopause layer (TTL). Here we quantify the relative roles of variations in TTL temperatures and transport in determining seasonal and interannual variations of stratospheric water vapour. Following previous studies, we use trajectory calculations with the water vapour concentration set by the Lagrangian dry point (LDP) along trajectories. To assess the separate roles of transport and temperatures, the LDP calculations are modified by replacing either the winds or the temperatures with those from different years to investigate the wind or temperature sensitivity of water vapour to interannual variations and, correspondingly, with those from different months to investigate the wind or temperature sensitivity to seasonal variations. Both ERA-Interim reanalysis data for the 1999–2009 period and data generated by a chemistry–climate model (UM-UKCA) are investigated. Variations in temperatures, rather than transport, dominate interannual variability, typically explaining more than 70 % of variability, including individual events such as the 2000 stratospheric water vapour drop. Similarly seasonal variation of temperatures, rather than transport, is shown to be the dominant driver of the annual cycle in lower stratospheric water vapour concentrations in both the model and reanalysis, but it is also shown that seasonal variation of transport plays an important role in reducing the seasonal cycle maximum (reducing the annual range by about 30 %). The quantitative role in dehydration of sub-seasonal and sub-monthly Eulerian temperature variability is also examined by using time-filtered temperature fields in the trajectory calculations. Sub-monthly temperature variability reduces annual mean water vapour concentrations by 40 % in the reanalysis calculation and 30 % in the model calculation. As with other aspects of dehydration, simple Eulerian measures of variability are not sufficient to quantify the implications for dehydration, and the Lagrangian sampling of the variability must be taken into account. These results indicate that, whilst capturing seasonal and interannual variation of temperature is a major factor in modelling realistic stratospheric water vapour concentrations, simulation of seasonal variation of transport and of sub-seasonal and sub-monthly temperature variability are also important and cannot be ignored.

We perform time-slice experiments using HadGEM3-A to decompose the long-term (years 101–150) response of the Brewer–Dobson circulation (BDC) to an abrupt quadrupling in CO2 (4×CO2) into (1) a rapid atmospheric adjustment, (2) a contribution from the global-average sea surface temperature (SST) change (+3.4 K), and (3) an SST pattern effect. The SST fields are derived from the CMIP5 multi-model ensemble. Two further experiments explore the impact on the BDC of the spread in global-average SST response to 4×CO2 across the CMIP5 models (range 2.1–4.9 K). At 70 hPa (10 hPa) the annual-mean tropical upward mass flux increases by 45 % (35 %) due to the 4×CO2 perturbation. At 70 hPa, around 70 % of the increase is from the global-uniform SST warming, with the remainder coming in similar contributions from the rapid adjustment and SST pattern effect. In contrast, at 10 hPa the increase in mass flux comes mainly from the rapid adjustment (∼40 %) and the uniform SST warming (∼45 %), with a small contribution from the SST pattern. At 10 hPa, the effect of the multi-model spread in global-mean SST is comparable in magnitude to the rapid adjustment. Conversely, at 70 hPa the effect of spread in global-mean SST is substantially larger than both the rapid adjustment and the SST pattern effect. We derive an approximately linear sensitivity of the tropical upward mass flux to global surface air temperature change of 0.62×109 kg s−1 K−1 (9 % K−1) at 70 hPa and 0.10×109 kg s−1 K−1 (6 % K−1) at 10 hPa. The results confirm the most important factor for the acceleration of the BDC in the lower stratosphere under increased CO2 is global SST changes. We also quantify for the first time that the rapid adjustment to CO2 is of similar importance to SSTs for the increased BDC in the upper stratosphere. This demonstrates a potential for a fast and slow timescale of the response of the BDC to greenhouse gas forcing, with the relative prominence of those timescales being height dependent.

We perform the first multi-model intercomparison of the impact of nudged meteorology on the stratospheric residual circulation using hindcast simulations from the Chemistry–Climate Model Initiative (CCMI). We examine simulations over the period 1980–2009 from seven models in which the meteorological fields are nudged towards a reanalysis dataset and compare these with their equivalent free-running simulations and the reanalyses themselves. We show that for the current implementations, nudging meteorology does not constrain the mean strength of the stratospheric residual circulation and that the inter-model spread is similar, or even larger, than in the free-running simulations. The nudged models generally show slightly stronger upwelling in the tropical lower stratosphere compared to the free-running versions and exhibit marked differences compared to the directly estimated residual circulation from the reanalysis dataset they are nudged towards. Downward control calculations applied to the nudged simulations reveal substantial differences between the climatological lower-stratospheric tropical upward mass flux (TUMF) computed from the modelled wave forcing and that calculated directly from the residual circulation. This explicitly shows that nudging decouples the wave forcing and the residual circulation so that the divergence of the angular momentum flux due to the mean motion is not balanced by eddy motions, as would typically be expected in the time mean. Overall, nudging meteorological fields leads to increased inter-model spread for most of the measures of the mean climatological stratospheric residual circulation assessed in this study. In contrast, the nudged simulations show a high degree of consistency in the inter-annual variability in the TUMF in the lower stratosphere, which is primarily related to the contribution to variability from the resolved wave forcing. The more consistent inter-annual variability in TUMF in the nudged models also compares more closely with the variability found in the reanalyses, particularly in boreal winter. We apply a multiple linear regression (MLR) model to separate the drivers of inter-annual and long-term variations in the simulated TUMF; this explains up to ∼75 % of the variance in TUMF in the nudged simulations. The MLR model reveals a statistically significant positive trend in TUMF for most models over the period 1980–2009. The TUMF trend magnitude is generally larger in the nudged models compared to their free-running counterparts, but the intermodel range of trends doubles from around a factor of 2 to a factor of 4 due to nudging. Furthermore, the nudged models generally do not match the TUMF trends in the reanalysis they are nudged towards for trends over different periods in the interval 1980–2009. Hence, we conclude that nudging does not strongly constrain long-term trends simulated by the chemistry–climate model (CCM) in the residual circulation. Our findings show that while nudged simulations may, by construction, produce accurate temperatures and realistic representations of fast horizontal transport, this is not typically the case for the slower zonal mean vertical transport in the stratosphere. Consequently, caution is required when using nudged simulations to interpret the behaviour of stratospheric tracers that are affected by the residual circulation.

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.

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 impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry–climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry–Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1–5 hPa) decreases by more than a factor of 2, from around 5 to 2 %, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55 K, or around 80 % of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper.

Over the 21st century changes in both tropospheric and stratospheric ozone are likely to have important consequences for the Earth's radiative balance. In this study, we investigate the radiative forcing from future ozone changes using the Community Earth System Model (CESM1), with the Whole Atmosphere Community Climate Model (WACCM), and including fully coupled radiation and chemistry schemes. Using year 2100 conditions from the Representative Concentration Pathway 8.5 (RCP8.5) scenario, we quantify the individual contributions to ozone radiative forcing of (1) climate change, (2) reduced concentrations of ozone depleting substances (ODSs), and (3) methane increases. We calculate future ozone radiative forcings and their standard error (SE; associated with inter-annual variability of ozone) relative to year 2000 of (1) 33 ± 104 m Wm−2, (2) 163 ± 109 m Wm−2, and (3) 238 ± 113 m Wm−2 due to climate change, ODSs, and methane, respectively. Our best estimate of net ozone forcing in this set of simulations is 430 ± 130 m Wm−2 relative to year 2000 and 760 ± 230 m Wm−2 relative to year 1750, with the 95 % confidence interval given by ±30 %. We find that the overall long-term tropospheric ozone forcing from methane chemistry–climate feedbacks related to OH and methane lifetime is relatively small (46 m Wm−2). Ozone radiative forcing associated with climate change and stratospheric ozone recovery are robust with regard to background climate conditions, even though the ozone response is sensitive to both changes in atmospheric composition and climate. Changes in stratospheric-produced ozone account for ∼ 50 % of the overall radiative forcing for the 2000–2100 period in this set of simulations, highlighting the key role of the stratosphere in determining future ozone radiative forcing.

The ozone radiative forcings (RFs) resulting from projected changes in climate, ozone-depleting substances (ODSs), non-methane ozone precursor emissions and methane between the years 2000 and 2100 are calculated using simulations from the UM-UKCA chemistry–climate model (UK Met Office's Unified Model containing the United Kingdom Chemistry and Aerosols sub-model). Projected measures to improve air-quality through reductions in non-methane tropospheric ozone precursor emissions present a co-benefit for climate, with a net global mean ozone RF of −0.09 W m−2. This is opposed by a positive ozone RF of 0.05 W m−2 due to future decreases in ODSs, which is driven by an increase in tropospheric ozone through stratosphere-to-troposphere transport of air containing higher ozone amounts. An increase in methane abundance by more than a factor of 2 (as projected by the RCP8.5 scenario) is found to drive an ozone RF of 0.18 W m−2, which would greatly outweigh the climate benefits of non-methane tropospheric ozone precursor reductions. A small fraction (∼ 15 %) of the ozone RF due to the projected increase in methane results from increases in stratospheric ozone. The sign of the ozone RF due to future changes in climate (including the radiative effects of greenhouse gases, sea surface temperatures and sea ice changes) is shown to be dependent on the greenhouse gas emissions pathway, with a positive RF (0.05 W m−2) for RCP4.5 and a negative RF (−0.07 W m−2) for the RCP8.5 scenario. This dependence arises mainly from differences in the contribution to RF from stratospheric ozone changes. Considering the increases in tropopause height under climate change causes only small differences (≤ |0.02| W m−2) for the stratospheric, tropospheric and whole-atmosphere RFs.

The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models.For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %).We compare the climatic effects of the CMIP6 solar forcing dataset to its CMIP5 predecessor by using time-slice experiments of two chemistry–climate models and a reference radiative transfer model. The differences in the long-term mean SSI in the CMIP6 dataset, compared to CMIP5, impact on climatological stratospheric conditions (lower shortwave heating rates of −0.35 K day−1 at the stratopause), cooler stratospheric temperatures (−1.5 K in the upper stratosphere), lower ozone abundances in the lower stratosphere (−3 %), and higher ozone abundances (+1.5 % in the upper stratosphere and lower mesosphere). Between the maximum and minimum phases of the 11-year solar cycle, there is an increase in shortwave heating rates (+0.2 K day−1 at the stratopause), temperatures ( ∼ 1 K at the stratopause), and ozone (+2.5 % in the upper stratosphere) in the tropical upper stratosphere using the CMIP6 forcing dataset. This solar-cycle response is slightly larger, but not statistically significantly different from that for the CMIP5 forcing dataset.CMIP6 models with a well-resolved shortwave radiation scheme are encouraged to prescribe SSI changes and include solar-induced stratospheric ozone variations, in order to better represent solar climate variability compared to models that only prescribe TSI and/or exclude the solar-ozone response. We show that monthly-mean solar-induced ozone variations are implicitly included in the SPARC/CCMI CMIP6 Ozone Database for historical simulations, which is derived from transient chemistry–climate model simulations and has been developed for climate models that do not calculate ozone interactively. CMIP6 models without chemistry that perform a preindustrial control simulation with time-varying solar forcing will need to use a modified version of the SPARC/CCMI Ozone Database that includes solar variability. CMIP6 models with interactive chemistry are also encouraged to use the particle forcing datasets, which will allow the potential long-term effects of particles to be addressed for the first time. The consideration of particle forcing has been shown to significantly improve the representation of reactive nitrogen and ozone variability in the polar middle atmosphere, eventually resulting in further improvements in the representation of solar climate variability in global models.

The structure and amplitude of the radiative contributions of the annual cycles in ozone and water vapour to the prominent annual cycle in temperatures in the tropical tropopause layer (TTL) are considered. This is done initially through a seasonally evolving fixed dynamical heating (SEFDH) calculation. The annual cycle in ozone is found to drive significant temperature changes predominantly locally (in the vertical) and roughly in phase with the observed TTL annual cycle. In contrast, temperature changes driven by the annual cycle in water vapour are out of phase with the latter. The effects are weaker than those of ozone but still quantitatively significant, particularly near the cold point (100 to 90 hPa) where there are substantial non-local effects from variations in water vapour in lower layers of the TTL. The combined radiative heating effect of the annual cycles in ozone and water vapour maximizes above the cold point and is one factor contributing to the vertical structure of the amplitude of the annual cycle in lower-stratospheric temperatures, which has a relatively localized maximum around 70 hPa. Other important factors are identified here: radiative damping timescales, which are shown to maximize over a deep layer centred on the cold point; the vertical structure of the dynamical heating; and non-radiative processes in the upper troposphere that are inferred to impose a strong constraint on tropical temperature perturbations below 130 hPa. The latitudinal structure of the radiative contributions to the annual cycle in temperatures is found to be substantially modified when the SEFDH assumption is relaxed and the dynamical response, as represented by a zonally symmetric calculation, is taken into account. The effect of the dynamical response is to reduce the strong latitudinal gradients and inter-hemispheric asymmetry seen in the purely radiative SEFDH temperature response, while leaving the 20° N–20° S average response relatively unchanged. The net contribution of the annual ozone and water vapour cycles to the peak-to-peak amplitude in the annual cycle of TTL temperatures is found to be around 35 % of the observed 8 K at 70 hPa, 40 % of 6 K at 90 hPa, and 45 % of 3 K at 100 hPa. The primary sensitivity of the calculated magnitude of the temperature response is identified as the assumed annual mean ozone mixing ratio in the TTL.

The Aerosol Chemistry Model Intercomparison Project (AerChemMIP) is endorsed by the Coupled-Model Intercomparison Project 6 (CMIP6) and is designed to quantify the climate and air quality impacts of aerosols and chemically reactive gases. These are specifically near-term climate forcers (NTCFs: methane, tropospheric ozone and aerosols, and their precursors), nitrous oxide and ozone-depleting halocarbons. The aim of AerChemMIP is to answer four scientific questions. 1. How have anthropogenic emissions contributed to global radiative forcing and affected regional climate over the historical period? 2. How might future policies (on climate, air quality and land use) affect the abundances of NTCFs and their climate impacts? 3.How do uncertainties in historical NTCF emissions affect radiative forcing estimates? 4. How important are climate feedbacks to natural NTCF emissions, atmospheric composition, and radiative effects? These questions will be addressed through targeted simulations with CMIP6 climate models that include an interactive representation of tropospheric aerosols and atmospheric chemistry. These simulations build on the CMIP6 Diagnostic, Evaluation and Characterization of Klima (DECK) experiments, the CMIP6 historical simulations, and future projections performed elsewhere in CMIP6, allowing the contributions from aerosols and/or chemistry to be quantified. Specific diagnostics are requested as part of the CMIP6 data request to highlight the chemical composition of the atmosphere, to evaluate the performance of the models, and to understand differences in behaviour between them.

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.

Changes in incoming solar ultraviolet radiation over the 11-year solar cycle affect stratospheric ozone abundances. It is important to quantify the magnitude, structure, and seasonality of the associated solar-ozone response (SOR) to understand the impact of the 11-year solar cycle on climate. Part 1 of this two-part study uses multiple linear regression analysis to extract the SOR in a number of recently updated satellite ozone datasets covering different periods within the epoch 1970 to 2013. The annual mean SOR in the updated version 7.0 (v7.0) Stratospheric Aerosol and Gas Experiment (SAGE) II number density dataset (1984–2004) is very consistent with that found in the previous v6.2. In contrast, we find a substantial decrease in the magnitude of the SOR in the tropical upper stratosphere in the SAGE II v7.0 mixing ratio dataset (∼ 1 %) compared to the v6.2 (∼ 4 %). This difference is shown to be largely attributable to the change in the independent stratospheric temperature dataset used to convert SAGE II ozone number densities to mixing ratios. Since these temperature records contain substantial uncertainties, we suggest that datasets based on SAGE II number densities are currently most reliable for evaluating the SOR. We further analyse three extended ozone datasets that combine SAGE II v7.0 number densities with more recent GOMOS (Global Ozone Monitoring by Occultation of Stars) or OSIRIS (Optical Spectrograph and Infrared Imager System) measurements. The extended SAGE–OSIRIS dataset (1984–2013) shows a smaller and less statistically significant SOR across much of the tropical upper stratosphere compared to the SAGE II data alone. In contrast, the two SAGE–GOMOS datasets (1984–2011) show SORs that are in closer agreement with the original SAGE II data and therefore appear to provide a more reliable estimate of the SOR. We also analyse the SOR in the recent Solar Backscatter Ultraviolet Instrument (SBUV) Merged Ozone Dataset (SBUVMOD) version 8.6 (VN8.6) (1970–2012) and SBUV Merged Cohesive VN8.6 (1978–2012) datasets and compare them to the previous SBUVMOD VN8.0 (1970–2009). Over their full lengths, the three records generally agree in terms of the broad magnitude and structure of the annual mean SOR. The main difference is that SBUVMOD VN8.6 shows a smaller and less significant SOR in the tropical upper stratosphere and therefore more closely resembles the SAGE II v7.0 mixing ratio data than does the SBUV Merged Cohesive VN8.6, which has a more continuous SOR of ∼ 2 % in this region. The sparse spatial and temporal sampling of limb satellite instruments prohibits the extraction of sub-annual variations in the SOR from SAGE-based datasets. However, the SBUVMOD VN8.6 dataset suggests substantial month-to-month variations in the SOR, particularly in the winter extratropics, which may be important for the proposed high-latitude dynamical response to the solar cycle. Overall, the results highlight substantial uncertainties in the magnitude and structure of the observed SOR from different satellite records. The implications of these uncertainties for understanding and modelling the effects of solar variability on climate should be explored.

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.

Trends in the vertical distribution of ozone are reported and compared for a number of new and recently revised data sets. The amount of ozone-depleting compounds in the stratosphere (as measured by equivalent effective stratospheric chlorine – EESC) was maximised in the second half of the 1990s. We examine the periods before and after the peak to see if any change in trend is discernible in the ozone record that might be attributable to a change in the EESC trend, though no attribution is attempted. Prior to 1998, trends in the upper stratosphere (~ 45 km, 4 hPa) are found to be −5 to −10 % per decade at mid-latitudes and closer to −5 % per decade in the tropics. No trends are found in the mid-stratosphere (28 km, 30 hPa). Negative trends are seen in the lower stratosphere at mid-latitudes in both hemispheres and in the deep tropics. However, it is hard to be categorical about the trends in the lower stratosphere for three reasons: (i) there are fewer measurements, (ii) the data quality is poorer, and (iii) the measurements in the 1990s are perturbed by aerosols from the Mt Pinatubo eruption in 1991. These findings are similar to those reported previously even though the measurements for the main satellite and ground-based records have been revised. There is no sign of a continued negative trend in the upper stratosphere since 1998: instead there is a hint of an average positive trend of ~ 2 % per decade in mid-latitudes and ~ 3 % per decade in the tropics. The significance of these upward trends is investigated using different assumptions of the independence of the trend estimates found from different data sets. The averaged upward trends are significant if the trends derived from various data sets are assumed to be independent (as in Pawson et al., 2014) but are generally not significant if the trends are not independent. This occurs because many of the underlying measurement records are used in more than one merged data set. At this point it is not possible to say which assumption is best. Including an estimate of the drift of the overall ozone observing system decreases the significance of the trends. The significance will become clearer as (i) more years are added to the observational record, (ii) further improvements are made to the historic ozone record (e.g. through algorithm development), and (iii) the data merging techniques are refined, particularly through a more rigorous treatment of uncertainties.

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.

It has been shown previously that one member of the Met Office Hadley Centre single-parameter perturbed physics ensemble – the so-called "low entrainment parameter" member – has a much higher climate sensitivity than other individual parameter perturbations. Here we show that the concentration of stratospheric water vapour in this member is over three times higher than observations, and, more importantly for climate sensitivity, increases significantly when climate warms. The large surface temperature response of this ensemble member is more consistent with stratospheric humidity change, rather than upper tropospheric clouds as has been previously suggested. The direct relationship between the bias in the control state (elevated stratospheric humidity) and the cause of the high climate sensitivity (a further increase in stratospheric humidity) lends further doubt as to the realism of this particular integration. This, together with other evidence, lowers the likelihood that the climate system's physical sensitivity is significantly higher than the likely upper range quoted in the Intergovernmental Panel on Climate Change's Fourth Assessment Report.