Dynamics, stratospheric ozone, and climate change

Dynamics affects the distribution and abundance of stratospheric ozone directly through transport of ozone itself and indirectly through its effect on ozone chemistry via temperature and transport of other chemical species. Dynamical processes must be considered in order to understand past ozone changes, especially in the northern hemisphere where there appears to be significant low-frequency variability which can look “trend-like” on decadal time scales. A major challenge is to quantify the predictable, or deterministic, component of past ozone changes. Over the coming century, changes in climate will affect the expected recovery of ozone. For policy reasons it is important to be able to distinguish and separately attribute the effects of ozone-depleting substances and greenhouse gases on both ozone and climate. While the radiative-chemical effects can be relatively easily identified, this is not so evident for dynamics — yet dynamical changes (e.g., changes in the Brewer-Dobson circulation) could have a first-order effect on ozone over particular regions. Understanding the predictability and robustness of such dynamical changes represents another major challenge. Chemistry-climate models have recently emerged as useful tools for addressing these questions, as they provide a self-consistent representation of dynamical aspects of climate and their coupling to ozone chemistry. We can expect such models to play an increasingly central role in the study of ozone and climate in the future, analogous to the central role of global climate models in the study of tropospheric climate change.

On the attribution of stratospheric ozone and temperature changes to changes in ozone-depleting substances and well-mixed greenhouse gases

The vertical profile of global-mean stratospheric temperature changes has traditionally represented an important diagnostic for the attribution of the cooling effects of stratospheric ozone depletion and CO2 increases. However, CO2-induced cooling alters ozone abundance by perturbing ozone chemistry, thereby coupling the stratospheric ozone and temperature responses to changes in CO2 and ozone-depleting substances (ODSs). Here we untangle the ozone-temperature coupling and show that the attribution of global-mean stratospheric temperature changes to CO2 and ODS changes (which are the true anthropogenic forcing agents) can be quite different from the traditional attribution to CO2 and ozone changes. The significance of these effects is quantified empirically using simulations from a three-dimensional chemistry-climate model. The results confirm the essential validity of the traditional approach in attributing changes during the past period of rapid ODS increases, although we find that about 10% of the upper stratospheric ozone decrease from ODS increases over the period 1975–1995 was offset by the increase in CO2, and the CO2-induced cooling in the upper stratosphere has been somewhat overestimated. When considering ozone recovery, however, the ozone-temperature coupling is a first-order effect; fully 2/5 of the upper stratospheric ozone increase projected to occur from 2010–2040 is attributable to CO2 increases. Thus, it has now become necessary to base attribution of global-mean stratospheric temperature changes on CO2 and ODS changes rather than on CO2 and ozone changes.

Impact of long-range correlations on trend detection in total ozone

Total ozone trends are typically studied using linear regression models that assume a first-order autoregression of the residuals [so-called AR(1) models]. We consider total ozone time series over 60°S–60°N from 1979 to 2005 and show that most latitude bands exhibit long-range correlated (LRC) behavior, meaning that ozone autocorrelation functions decay by a power law rather than exponentially as in AR(1). At such latitudes the uncertainties of total ozone trends are greater than those obtained from AR(1) models and the expected time required to detect ozone recovery correspondingly longer. We find no evidence of LRC behavior in southern middle-and high-subpolar latitudes (45°–60°S), where the long-term ozone decline attributable to anthropogenic chlorine is the greatest. We thus confirm an earlier prediction based on an AR(1) analysis that this region (especially the highest latitudes, and especially the South Atlantic) is the optimal location for the detection of ozone recovery, with a statistically significant ozone increase attributable to chlorine likely to be detectable by the end of the next decade. In northern middle and high latitudes, on the other hand, there is clear evidence of LRC behavior. This increases the uncertainties on the long-term trend attributable to anthropogenic chlorine by about a factor of 1.5 and lengthens the expected time to detect ozone recovery by a similar amount (from ∼2030 to ∼2045). If the long-term changes in ozone are instead fit by a piecewise-linear trend rather than by stratospheric chlorine loading, then the strong decrease of northern middle- and high-latitude ozone during the first half of the 1990s and its subsequent increase in the second half of the 1990s projects more strongly on the trend and makes a smaller contribution to the noise. This both increases the trend and weakens the LRC behavior at these latitudes, to the extent that ozone recovery (according to this model, and in the sense of a statistically significant ozone increase) is already on the verge of being detected. The implications of this rather controversial interpretation are discussed.

Multimodel projections of stratospheric ozone in the 21st century

Simulations from eleven coupled chemistry-climate models (CCMs) employing nearly identical forcings have been used to project the evolution of stratospheric ozone throughout the 21st century. The model-to-model agreement in projected temperature trends is good, and all CCMs predict continued, global mean cooling of the stratosphere over the next 5 decades, increasing from around 0.25 K/decade at 50 hPa to around 1 K/ decade at 1 hPa under the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A1B scenario. In general, the simulated ozone evolution is mainly determined by decreases in halogen concentrations and continued cooling of the global stratosphere due to increases in greenhouse gases (GHGs). Column ozone is projected to increase as stratospheric halogen concentrations return to 1980s levels. Because of ozone increases in the middle and upper stratosphere due to GHGinduced cooling, total ozone averaged over midlatitudes, outside the polar regions, and globally, is projected to increase to 1980 values between 2035 and 2050 and before lower stratospheric halogen amounts decrease to 1980 values. In the polar regions the CCMs simulate small temperature trends in the first and second half of the 21st century in midwinter. Differences in stratospheric inorganic chlorine (Cly) among the CCMs are key to diagnosing the intermodel differences in simulated ozone recovery, in particular in the Antarctic. It is found that there are substantial quantitative differences in the simulated Cly, with the October mean Antarctic Cly peak value varying from less than 2 ppb to over 3.5 ppb in the CCMs, and the date at which the Cly returns to 1980 values varying from before 2030 to after 2050. There is a similar variation in the timing of recovery of Antarctic springtime column ozone back to 1980 values. As most models underestimate peak Cly near 2000, ozone recovery in the Antarctic could occur even later, between 2060 and 2070. In the Arctic the column ozone increase in spring does not follow halogen decreases as closely as in the Antarctic, reaching 1980 values before Arctic halogen amounts decrease to 1980 values and before the Antarctic. None of the CCMs predict future large decreases in the Arctic column ozone. By 2100, total column ozone is projected to be substantially above 1980 values in all regions except in the tropics.

Using transport diagnostics to understand chemistry climate model ozone simulations

We use observations of N2O and mean age to identify realistic transport in models in order to explain their ozone predictions. The results are applied to 15 chemistry climate models (CCMs) participating in the 2010 World Meteorological Organization ozone assessment. Comparison of the observed and simulated N2O, mean age and their compact correlation identifies models with fast or slow circulations and reveals details of model ascent and tropical isolation. This process‐oriented diagnostic is more useful than mean age alone because it identifies models with compensating transport deficiencies that produce fortuitous agreement with mean age. The diagnosed model transport behavior is related to a model’s ability to produce realistic lower stratosphere (LS) O3 profiles. Models with the greatest tropical transport problems compare poorly with O3 observations. Models with the most realistic LS transport agree more closely with LS observations and each other. We incorporate the results of the chemistry evaluations in the Stratospheric Processes and their Role in Climate (SPARC) CCMVal Report to explain the range of CCM predictions for the return‐to‐1980 dates for global (60°S–60°N) and Antarctic column ozone. Antarctic O3 return dates are generally correlated with vortex Cly levels, and vortex Cly is generally correlated with the model’s circulation, although model Cl chemistry and conservation problems also have a significant effect on return date. In both regions, models with good LS transport and chemistry produce a smaller range of predictions for the return‐to‐1980 ozone values. This study suggests that the current range of predicted return dates is unnecessarily broad due to identifiable model deficiencies.

Multimodel assessment of the factors driving stratospheric ozone evolution over the 21st century

The evolution of stratospheric ozone from 1960 to 2100 is examined in simulations from 14 chemistry‐climate models, driven by prescribed levels of halogens and greenhouse gases. There is general agreement among the models that total column ozone reached a minimum around year 2000 at all latitudes, projected to be followed by an increase over the first half of the 21st century. In the second half of the 21st century, ozone is projected to continue increasing, level off, or even decrease depending on the latitude. Separation into partial columns above and below 20 hPa reveals that these latitudinal differences are almost completely caused by differences in the model projections of ozone in the lower stratosphere. At all latitudes, upper stratospheric ozone increases throughout the 21st century and is projected to return to 1960 levels well before the end of the century, although there is a spread among models in the dates that ozone returns to specific historical values. We find decreasing halogens and declining upper atmospheric temperatures, driven by increasing greenhouse gases, contribute almost equally to increases in upper stratospheric ozone. In the tropical lower stratosphere, an increase in upwelling causes a steady decrease in ozone through the 21st century, and total column ozone does not return to 1960 levels in most of the models. In contrast, lower stratospheric and total column ozone in middle and high latitudes increases during the 21st century, returning to 1960 levels well before the end of the century in most models.

Impact of sudden Arctic sea-ice loss on stratospheric polar ozone recovery

We investigate the sensitivity of Northern Hemisphere polar ozone recovery to a scenario in which there is rapid loss of Arctic summer sea ice in the first half of the 21st century. The issue is addressed by coupling a chemistry climate model to an ocean general circulation model and performing simulations of ozone recovery with, and without, an external perturbation designed to cause a rapid and complete loss of summertime Arctic sea ice. Under this extreme perturbation, the stratospheric response takes the form of a springtime polar cooling which is dynamical rather than radiative in origin, and is caused by reduced wave forcing from the troposphere. The response lags the onset of the sea-ice perturbation by about one decade and lasts for more than two decades, and is associated with an enhanced weakening of the North Atlantic meridional overturning circulation. The stratospheric dynamical response leads to a 10 DU reduction in polar column ozone, which is statistically robust. While this represents a modest loss, it has the potential to induce a delay of roughly one decade in Arctic ozone recovery estimates made in the 2006 Scientific Assessment of Ozone Depletion.

Summertime total ozone variations over middle and polar latitudes

The statistical relationship between springtime and summertime ozone over middle and polar latitudes is analyzed using zonally averaged total ozone data. Shortterm variations in springtime midlatitude ozone demonstrate only a modest correlation with springtime polar ozone variations. However by early summer, ozone variations throughout the extratropics are highly correlated. Analysis of correlation functions indicates that springtime midlatitude ozone, not polar ozone, is the best predictor for summertime polar ozone. Long-term total ozone trends at middle and high latitudes are also different for spring and nearly identical for summer. About 39% of the observed southern midlatitude ozone decline in December can be attributed to the polar ozone depletion up to November. In the Northern Hemisphere, the corresponding contribution is about 15%, but the error bars are too large to make an accurate estimate.

Ozone database in support of CMIP5 simulations: results and corresponding radiative forcing

A continuous tropospheric and stratospheric vertically resolved ozone time series, from 1850 to 2099, has been generated to be used as forcing in global climate models that do not include interactive chemistry. A multiple linear regression analysis of SAGE I+II satellite observations and polar ozonesonde measurements is used for the stratospheric zonal mean dataset during the well-observed period from 1979 to 2009. In addition to terms describing the mean annual cycle, the regression includes terms representing equivalent effective stratospheric chlorine (EESC) and the 11-yr solar cycle variability. The EESC regression fit coefficients, together with pre-1979 EESC values, are used to extrapolate the stratospheric ozone time series backward to 1850. While a similar procedure could be used to extrapolate into the future, coupled chemistry climate model (CCM) simulations indicate that future stratospheric ozone abundances are likely to be significantly affected by climate change, and capturing such effects through a regression model approach is not feasible. Therefore, the stratospheric ozone dataset is extended into the future (merged in 2009) with multimodel mean projections from 13 CCMs that performed a simulation until 2099 under the SRES (Special Report on Emission Scenarios) A1B greenhouse gas scenario and the A1 adjusted halogen scenario in the second round of the Chemistry-Climate Model Validation (CCMVal-2) Activity. The stratospheric zonal mean ozone time series is merged with a three-dimensional tropospheric data set extracted from simulations of the past by two CCMs (CAM3.5 and GISSPUCCINI)and of the future by one CCM (CAM3.5). The future tropospheric ozone time series continues the historical CAM3.5 simulation until 2099 following the four different Representative Concentration Pathways (RCPs). Generally good agreement is found between the historical segment of the ozone database and satellite observations, although it should be noted that total column ozone is overestimated in the southern polar latitudes during spring and tropospheric column ozone is slightly underestimated. Vertical profiles of tropospheric ozone are broadly consistent with ozonesondes and in-situ measurements, with some deviations in regions of biomass burning. The tropospheric ozone radiative forcing (RF) from the 1850s to the 2000s is 0.23Wm−2, lower than previous results. The lower value is mainly due to (i) a smaller increase in biomass burning emissions; (ii) a larger influence of stratospheric ozone depletion on upper tropospheric ozone at high southern latitudes; and possibly (iii) a larger influence of clouds (which act to reduce the net forcing) compared to previous radiative forcing calculations. Over the same period, decreases in stratospheric ozone, mainly at high latitudes, produce a RF of −0.08Wm−2, which is more negative than the central Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) value of −0.05Wm−2, but which is within the stated range of −0.15 to +0.05Wm−2. The more negative value is explained by the fact that the regression model simulates significant ozone depletion prior to 1979, in line with the increase in EESC and as confirmed by CCMs, while the AR4 assumed no change in stratospheric RF prior to 1979. A negative RF of similar magnitude persists into the future, although its location shifts from high latitudes to the tropics. This shift is due to increases in polar stratospheric ozone, but decreases in tropical lower stratospheric ozone, related to a strengthening of the Brewer-Dobson circulation, particularly through the latter half of the 21st century. Differences in trends in tropospheric ozone among the four RCPs are mainly driven by different methane concentrations, resulting in a range of tropospheric ozone RFs between 0.4 and 0.1Wm−2 by 2100. The ozone dataset described here has been released for the Coupled Model Intercomparison Project (CMIP5) model simulations in netCDF Climate and Forecast (CF) Metadata Convention at the PCMDI website (http://cmip-pcmdi.llnl.gov/).

Northern winter stratospheric temperature and ozone responses to ENSO inferred from an ensemble of Chemistry Climate Models

The connection between the El Ni˜no Southern Oscillation (ENSO) and the Northern polar stratosphere has been established from observations and atmospheric modeling. Here a systematic inter-comparison of the sensitivity of the modeled stratosphere to ENSO in Chemistry Climate Models (CCMs) is reported. This work uses results from a number of the CCMs included in the 2006 ozone assessment. In the lower stratosphere, the mean of all model simulations reports a warming of the polar vortex during strong ENSO events in February–March, consistent with but smaller than the estimate from satellite observations and ERA40 reanalysis. The anomalous warming is associated with an anomalous dynamical increase of column ozone north of 70� N that is accompanied by coherent column ozone decrease in the Tropics, in agreement with that deduced from the NIWA column ozone database, implying an increased residual circulation in the mean of all model simulations during ENSO. The spread in the model responses is partly due to the large internal stratospheric variability and it is shown that it crucially depends on the representation of the tropospheric ENSO teleconnection in the models.

Persistence and photochemical decay of springtime total ozone anomalies in the Canadian Middle Atmosphere Model

The persistence and decay of springtime total ozone anomalies over the entire extratropics (midlatitudes plus polar regions) is analysed using results from the Canadian Middle Atmosphere Model (CMAM), a comprehensive chemistry-climate model. As in the observations, interannual anomalies established through winter and spring persist with very high correlation coefficients (above 0.8) through summer until early autumn, while decaying in amplitude as a result of photochemical relaxation in the quiescent summertime stratosphere. The persistence and decay of the ozone anomalies in CMAM agrees extremely well with observations, even in the southern hemisphere when the model is run without heterogeneous chemistry (in which case there is no ozone hole and the seasonal cycle of ozone is quite different from observations). However in a version of CMAM with strong vertical diffusion, the northern hemisphere anomalies decay far too rapidly compared to observations. This shows that ozone anomaly persistence and decay does not depend on how the springtime anomalies are created or on their magnitude, but reflects the transport and photochemical decay in the model. The seasonality of the long-term trends over the entire extratropics is found to be explained by the persistence of the interannual anomalies, as in the observations, demonstrating that summertime ozone trends reflect winter/spring trends rather than any change in summertime ozone chemistry. However this mechanism fails in the northern hemisphere midlatitudes because of the relatively large impact, compared to observations, of the CMAM polar anomalies. As in the southern hemisphere, the influence of polar ozone loss in CMAM increases the midlatitude summertime loss, leading to a relatively weak seasonal dependence of ozone loss in the Northern Hemisphere compared to the observations.

Assessment of temperature, trace species, and ozone in chemistry-climate model simulations of the recent past

Simulations of the stratosphere from thirteen coupled chemistry-climate models (CCMs) are evaluated to provide guidance for the interpretation of ozone predictions made by the same CCMs. The focus of the evaluation is on how well the fields and processes that are important for determining the ozone distribution are represented in the simulations of the recent past. The core period of the evaluation is from 1980 to 1999 but long-term trends are compared for an extended period (1960–2004). Comparisons of polar high-latitude temperatures show that most CCMs have only small biases in the Northern Hemisphere in winter and spring, but still have cold biases in the Southern Hemisphere spring below 10 hPa. Most CCMs display the correct stratospheric response of polar temperatures to wave forcing in the Northern, but not in the Southern Hemisphere. Global long-term stratospheric temperature trends are in reasonable agreement with satellite and radiosonde observations. Comparisons of simulations of methane, mean age of air, and propagation of the annual cycle in water vapor show a wide spread in the results, indicating differences in transport. However, for around half the models there is reasonable agreement with observations. In these models the mean age of air and the water vapor tape recorder signal are generally better than reported in previous model intercomparisons. Comparisons of the water vapor and inorganic chlorine (Cly) fields also show a large intermodel spread. Differences in tropical water vapor mixing ratios in the lower stratosphere are primarily related to biases in the simulated tropical tropopause temperatures and not transport. The spread in Cly, which is largest in the polar lower stratosphere, appears to be primarily related to transport differences. In general the amplitude and phase of the annual cycle in total ozone is well simulated apart from the southern high latitudes. Most CCMs show reasonable agreement with observed total ozone trends and variability on a global scale, but a greater spread in the ozone trends in polar regions in spring, especially in the Arctic. In conclusion, despite the wide range of skills in representing different processes assessed here, there is sufficient agreement between the majority of the CCMs and the observations that some confidence can be placed in their predictions.

Ozone and climate: a review of interconnections

Stratospheric ozone has been depleted over the last 25 years following anthropogenic emissions of a number of chlorine- and bromine-containing compounds (ozone-depleting substances, ODSs), which are now regulated under the Montreal Protocol. The Protocol has been effective in controlling the net growth of these compounds in the atmosphere. As chlorine and bromine slowly decrease in the future, ozone levels are expected to increase in the coming decades, although the evolution will also depend on the changing climate system.

Seasonal persistence of midlatitude total ozone anomalies

Temporal autocorrelations of monthly mean total ozone anomalies over the 35–60°S and 35–60°N latitude bands reveal that anomalies established in the wintertime midlatitude ozone buildup persist (with photochemical decay) until the end of the following autumn, and then are rapidly erased once the next winter's buildup begins. The photochemical decay rate is found to be identical between the two hemispheres. High predictability of ozone through late summer exists based on the late-spring values. In the northern hemisphere, extending the 1979–2001 springtime ozone trend to other months through regression based on the seasonal persistence of anomalies captures the seasonality of the ozone trends remarkably well. In the southern hemisphere, the springtime trend only accounts for part of the summertime trends. There is a strong correlation between the ozone anomalies in northern hemisphere spring and those in the subsequent southern hemisphere spring, but not the converse.

Arctic ozone depletion in 2002–2003 measured by ASUR and comparison with POAM observations

We present ozone loss estimated from airborne measurements taken during January–February and March in the Arctic winter 2002/2003. The first half of the winter was characterized by unusually cold temperatures and the second half by a major stratospheric sudden warming around 15–18 January 2003. The potential vorticity maps show a vortex split in the lower stratosphere during the major warming (MW) in late January and during the minor warming in mid-February due to wave 1 amplification. However, the warming can be termed as a vortex displacement event as there was no vortex split during the MW period at 10 hPa. Very low temperatures, large areas of polar stratospheric clouds (PSCs), and high chlorine activation triggered significant ozone loss in the early winter, as the vortex moved to the midlatitude regions. The ozone depletion derived from the ASUR measurements sampled inside the vortex, in conjunction with the Mimosa-Chim model tracer, shows a maximum of 1.3 ± 0.2 ppmv at 450–500 K by late March. The partial column loss derived from the ASUR ozone profiles reaches up to 61 ± 4 DU in 400–550 K in the same period. The evolution of ozone and ozone loss assessed from the ASUR measurements is in very good agreement with POAM observations. The reduction in ozone estimated from the POAM measurements shows a similar maximum of 1.3 ± 0.2 ppmv at 400–500 K or 63 ± 4 DU in 400–550 K in late March. Our study reveals that the Arctic winter 2002/2003 was unique as it had three minor warmings and a MW, yet showed large loss in ozone. No such feature was observed in any other Arctic winter in the 1989–2010 period. In addition, an unusually large ozone loss in December, around 0.5 ± 0.2 ppmv at 450–500 K or 12 ± 1 DU in 400–550 K, was estimated for the first time in the Arctic. A careful and detailed diagnosis with all available published results for this winter exhibits an average ozone loss of 1.5 ± 0.3 ppmv at 450–500 K or 65 ± 5 DU in 400–550 K by the end of March, which exactly matches the ozone depletion derived from the ASUR, POAM and model data. The early ozone loss together with considerable loss afterwards put the warm Arctic winter 2002/2003 amongst the moderately cold winters in terms of the significance of the ozone loss.