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.

Sudden-stratospheric-warming (SSW) events are often followed by significant weather and climate impacts at the surface. By affecting the North Atlantic Oscillation (NAO), SSWs can lead to periods of extreme cold in parts of Europe and North America. Previous studies have used observations and free-running climate models to try to identify features of the atmosphere prior to an SSW that can determine the subsequent impact at the surface. However, the limited observational record makes it difficult to accurately quantify these relationships. Here, we instead use a large ensemble of seasonal hindcasts. We first test whether the hindcasts reproduce the observed characteristics of SSWs and their surface signature. We find that the simulations are statistically indistinguishable from the observations, in terms of the overall risk of an SSW per winter (56 %), the frequency of SSWs with negative NAO responses (65 %), the magnitude of the NAO responses, and the frequency of wavenumber-2-dominated SSWs (26 %). We also assess the relationships between prior conditions and the NAO response in the 30 d following an SSW. We find that there is little information in the precursor state to guide differences in the subsequent NAO behaviour between one SSW and another, reflecting the substantial natural variability between SSW events. The strongest relationships with the NAO response are from pre-SSW sea level pressure anomalies over the polar cap and from zonal-wind anomalies in the lower stratosphere, both exhibiting correlations of around 0.3. The pre-SSW NAO has little bearing on its post-SSW state. The strength of the pre-SSW zonal-wind anomalies at 10 hPa is also not significantly correlated with the NAO response. Finally, we find that the mean NAO response in the first 10 d following wave-2-dominated SSWs is much more strongly negative than in wave-1 cases. However, the subsequent response in days 11–30 is very similar regardless of the dominant wavenumber. In all cases, the composite mean responses are the result of very broad distributions from individual SSW events, necessitating a probabilistic analysis using large ensembles.

The degree to which Arctic sea ice decline influences the mid-latitude atmospheric circulation is widely debated. Here, the authors use a coordinated multi-model experiment to show that Arctic sea ice loss causes a weakening of the mid-latitude westerly winds, but the effect is overall small.

Over recent years there have been concomitant advances in the development of stratosphere-resolving numerical models, our understanding of stratosphere–troposphere interaction, and the extension of long-range forecasts to explicitly include the stratosphere. These advances are now allowing for new and improved capability in long-range prediction. We present an overview of this development and show how the inclusion of the stratosphere in forecast systems aids monthly, seasonal, and annual-to-decadal climate predictions and multidecadal projections. We end with an outlook towards the future and identify areas of improvement that could further benefit these rapidly evolving predictions.

Abstract It is well established that the positive phase of El Niño Southern Oscillation (ENSO) tends to weaken the Northern Hemisphere stratospheric polar vortex (SPV), promoting a negative North Atlantic Oscillation (NAO). Pacific Decadal Variability (PDV) is characterized by a pattern of sea surface temperatures similar to ENSO, but its impacts are more uncertain: some studies suggest similar impacts of ENSO and PDV on the SPV and NAO, while others find the opposite. We use climate model experiments and reanalysis to find further evidence supporting opposite interannual and decadal impacts of Pacific variability on the extratropics. We propose that the decadal strengthening of the SPV in response to positive PDV is caused by a build‐up of stratospheric water vapor leading to enhanced cooling at the poles, an increased meridional temperature gradient and a strengthened extratropical jet. Our results are important for understanding decadal variability, seasonal to decadal forecasts and climate projections.

The Brewer–Dobson circulation (BDC) is a key feature of the stratosphere that models need to accurately represent in order to simulate surface climate variability and change adequately. For the first time, the Climate Model Intercomparison Project includes in its phase 6 (CMIP6) a set of diagnostics that allow for careful evaluation of the BDC. Here, the BDC is evaluated against observations and reanalyses using historical simulations. CMIP6 results confirm the well-known inconsistency in the sign of BDC trends between observations and models in the middle and upper stratosphere. Nevertheless, the large uncertainty in the observational trend estimates opens the door to compatibility. In particular, when accounting for the limited sampling of the observations, model and observational trend error bars overlap in 40 % of the simulations with available output. The increasing CO2 simulations feature an acceleration of the BDC but reveal a large spread in the middle-to-upper stratospheric trends, possibly related to the parameterized gravity wave forcing. The very close connection between the shallow branch of the residual circulation and surface temperature is highlighted, which is absent in the deep branch. The trends in mean age of air are shown to be more robust throughout the stratosphere than those in the residual circulation.

Major mid-winter stratospheric sudden warmings (SSWs) are the largest instance of wintertime variability in the Arctic stratosphere. Because SSWs are able to cause significant surface weather anomalies on intra-seasonal timescales, several previous studies have focused on their potential future change, as might be induced by anthropogenic forcings. However, a wide range of results have been reported, from a future increase in the frequency of SSWs to an actual decrease. Several factors might explain these contradictory results, notably the use of different metrics for the identification of SSWs and the impact of large climatological biases in single-model studies. To bring some clarity, we here revisit the question of future SSW changes, using an identical set of metrics applied consistently across 12 different models participating in the Chemistry–Climate Model Initiative. Our analysis reveals that no statistically significant change in the frequency of SSWs will occur over the 21st century, irrespective of the metric used for the identification of the event. Changes in other SSW characteristics – such as their duration, deceleration of the polar night jet, and the tropospheric forcing – are also assessed: again, we find no evidence of future changes over the 21st century.

>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.

Free-running and nudged versions of a Met Office chemistry–climate model are evaluated and used to investigate the impact of dynamics versus transport and chemistry within the model on the simulated evolution of stratospheric ozone. Metrics of the dynamical processes relevant for simulating stratospheric ozone are calculated, and the free-running model is found to outperform the previous model version in 10 of the 14 metrics. In particular, large biases in stratospheric transport and tropical tropopause temperature, which existed in the previous model version, are substantially reduced, making the current model more suitable for the simulation of stratospheric ozone. The spatial structure of the ozone hole, the area of polar stratospheric clouds, and the increased ozone concentrations in the Northern Hemisphere winter stratosphere following sudden stratospheric warmings, were all found to be sensitive to the accuracy of the dynamics and were better simulated in the nudged model than in the free-running model. Whilst nudging can, in general, provide a useful tool for removing the influence of dynamical biases from the evolution of chemical fields, this study shows that issues can remain in the climatology of nudged models. Significant biases in stratospheric vertical velocities, age of air, water vapour, and total column ozone still exist in the Met Office nudged model. Further, these can lead to biases in the downward flux of ozone into the troposphere.

We present an overview of state-of-the-art chemistry–climate and chemistry transport models that are used within phase 1 of the Chemistry–Climate Model Initiative (CCMI-1). The CCMI aims to conduct a detailed evaluation of participating models using process-oriented diagnostics derived from observations in order to gain confidence in the models' projections of the stratospheric ozone layer, tropospheric composition, air quality, where applicable global climate change, and the interactions between them. Interpretation of these diagnostics requires detailed knowledge of the radiative, chemical, dynamical, and physical processes incorporated in the models. Also an understanding of the degree to which CCMI-1 recommendations for simulations have been followed is necessary to understand model responses to anthropogenic and natural forcing and also to explain inter-model differences. This becomes even more important given the ongoing development and the ever-growing complexity of these models. This paper also provides an overview of the available CCMI-1 simulations with the aim of informing CCMI data users.

We describe Global Atmosphere 4.0 (GA4.0) and Global Land 4.0 (GL4.0): configurations of the Met Office Unified Model and JULES (Joint UK Land Environment Simulator) community land surface model developed for use in global and regional climate research and weather prediction activities. GA4.0 and GL4.0 are based on the previous GA3.0 and GL3.0 configurations, with the inclusion of developments made by the Met Office and its collaborators during its annual development cycle. This paper provides a comprehensive technical and scientific description of GA4.0 and GL4.0 as well as details of how these differ from their predecessors. We also present the results of some initial evaluations of their performance. Overall, performance is comparable with that of GA3.0/GL3.0; the updated configurations include improvements to the science of several parametrisation schemes, however, and will form a baseline for further ongoing development.

We describe the HadGEM2 family of climate configurations of the Met Office Unified Model, MetUM. The concept of a model "family" comprises a range of specific model configurations incorporating different levels of complexity but with a common physical framework. The HadGEM2 family of configurations includes atmosphere and ocean components, with and without a vertical extension to include a well-resolved stratosphere, and an Earth-System (ES) component which includes dynamic vegetation, ocean biology and atmospheric chemistry. The HadGEM2 physical model includes improvements designed to address specific systematic errors encountered in the previous climate configuration, HadGEM1, namely Northern Hemisphere continental temperature biases and tropical sea surface temperature biases and poor variability. Targeting these biases was crucial in order that the ES configuration could represent important biogeochemical climate feedbacks. Detailed descriptions and evaluations of particular HadGEM2 family members are included in a number of other publications, and the discussion here is limited to a summary of the overall performance using a set of model metrics which compare the way in which the various configurations simulate present-day climate and its variability.

The scientific understanding of the Earth's climate system, including the central question of how the climate system is likely to respond to human-induced perturbations, is comprehensively captured in GCMs and Earth System Models (ESM). Diagnosing the simulated climate response, and comparing responses across different models, is crucially dependent on transparent assumptions of how the GCM/ESM has been driven – especially because the implementation can involve subjective decisions and may differ between modelling groups performing the same experiment. This paper outlines the climate forcings and setup of the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document the prescribed greenhouse gas concentrations, aerosol precursors, stratospheric and tropospheric ozone assumptions, as well as implementation of land-use change and natural forcings for the HadGEM2-ES historical and future experiments following the Representative Concentration Pathways. In addition, we provide details of how HadGEM2-ES ensemble members were initialised from the control run and how the palaeoclimate and AMIP experiments, as well as the "emission-driven" RCP experiments were performed.

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.

Sudden stratospheric warmings (SSWs) are amongst the most dramatic events in the Earth’s atmosphere and they drive extreme surface weather conditions. They have been recently linked to the hot and dry weather conditions that favour wildfires over Australia. However, the chance of a southern hemisphere event is unknown because it has only been observed once. Legitimate estimation of the frequency of SSW events requires a large sample of realistic model simulations. Here we show that the chance of an event is close to 4% per year, implying that an event will occur, on average, every 25 yr, using a state-of-the-art model that simulates SSWs accurately. It is thus not surprising that there was a near miss in the September prior to the Australian wildfire of 2019, given the 40 yr of comprehensive satellite records and just one observed Antarctic event. According to this new estimate, it would also not be surprising to see a second SSW event in the coming years in the southern hemisphere. Such a stratospheric warming event might bring further extreme surface weather conditions and natural hazards, as it may raise the risk of increased rainfall in the latitudinal band of 35–50°S. Meanwhile, the associated hot and dry weather conditions over austral subtropical continents might increase the risk of wildfires over these regions.