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.

Large-scale dynamics of the mesosphere and lower thermosphere: an analysis using the extended Canadian Middle Atmosphere Model

The extended Canadian Middle Atmosphere Model is used to investigate the large-scale dynamics of the mesosphere and lower thermosphere (MLT). It is shown that the 4-day wave is substantially amplified in southern polar winter in the presence of instabilities arising from strong vertical shears in the MLT zonal mean zonal winds brought about by parameterized nonorographic gravity wave drag. A weaker 4-day wave in northern polar winter is attributed to the weaker wind shears that result from weaker parameterized wave drag. The 2-day wave also exhibits a strong dependence on zonal wind shears, in agreement with previous modeling studies. In the equatorial upper mesosphere, the migrating diurnal tide provides most of the resolved westward wave forcing, which varies semiannually in conjunction with the tide itself; resolved forcing by eastward traveling disturbances is dominated by smaller scales. Nonmigrating tides and other planetary-scale waves play only a minor role in the zonal mean zonal momentum budget in the tropics at these heights. Resolved waves are shown to play a significant role in the zonal mean meridional momentum budget in the MLT, impacting significantly on gradient wind balance. Balance fails at low latitudes as a result of a strong Reynolds stress associated with the migrating diurnal tide, an effect which is most pronounced at equinox when the tide is strongest. Resolved and parameterized waves account for most of the imbalance at higher latitudes in summer. This results in the gradient wind underestimating the actual eastward wind reversal by up to 40%.

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.

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.

The tropopause inversion layer in models and analyses

Recent high-resolution radiosonde climatologies have revealed a tropopause inversion layer (TIL) in the extratropics: temperature strongly increases just above a sharp local cold point tropopause. Here, it is asked to what extent a TIL exists in current general circulation models (GCMs) and meteorological analyses. Only a weak hint of a TIL exists in NCEP/NCAR reanalysis data. In contrast, the Canadian Middle Atmosphere Model (CMAM), a comprehensive GCM, exhibits a TIL of realistic strength. However, in data assimilation mode CMAM exhibits a much weaker TIL, especially in the Southern Hemisphere where only coarse satellite data are available. The discrepancy between the analyses and the GCM is thus hypothesized to be mainly due to data assimilation acting to smooth the observed strong curvature in temperature around the tropopause. This is confirmed in the reanalysis where the stratification around the tropopause exhibits a strong discontinuity at the start of the satellite era.

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.

Correlations of long-lived chemical species in a middle atmosphere general circulation model

Correlations between various chemical species simulated by the Canadian Middle Atmosphere Model, a general circulation model with fully interactive chemistry, are considered in order to investigate the general conditions under which compact correlations can be expected to form. At the same time, the analysis serves to validate the model. The results are compared to previous work on this subject, both from theoretical studies and from atmospheric measurements made from space and from aircraft. The results highlight the importance of having a data set with good spatial coverage when working with correlations and provide a background against which the compactness of correlations obtained from atmospheric measurements can be confirmed. It is shown that for long-lived species, distinct correlations are found in the model in the tropics, the extratropics, and the Antarctic winter vortex. Under these conditions, sparse sampling such as arises from occultation instruments is nevertheless suitable to define a chemical correlation within each region even from a single day of measurements, provided a sufficient range of mixing ratio values is sampled. In practice, this means a large vertical extent, though the requirements are less stringent at more poleward latitudes.

Flow-dependent predictability of the North Atlantic Jet

The North Atlantic eddy-driven jet is a major component of the large-scale flow in the northern hemisphere. Here we present evidence from reanalysis and ensemble forecast data for systematic flow-dependent predictability of the jet during northern hemisphere winter (DJF). It is found that when the jet is weakened or split it is both less persistent and less predictable. The lack of predictability manifests itself as the onset of an anomalously large instantaneous rate of spread of ensemble forecast members as the jet becomes weakened. This suggests that as the jet weakens or splits it enters into a state more sensitive to small differences between ensemble forecast members, rather like the sensitive region between the wings of the Lorenz attractor.

Sensitivity of biogenic isoprene emissions to past, present and future environmental conditions and implications for atmospheric chemistry.

[1] We have implemented a process-based isoprene emission model in the HadGEM2 Earth-system model with coupled atmospheric chemistry in order to examine the feedback between isoprene emission and climate. Isoprene emissions and their impact on atmospheric chemistry and climate are estimated for preindustrial (1860–1869), present-day (2000–2009), and future (2100–2109) climate conditions. The estimate of 460 TgC/yr for present-day global total isoprene emission is consistent with previous estimates. Preindustrial isoprene emissions are estimated to be 26% higher than present-day. Future isoprene emissions using the RCP8.5 scenario are similar to present-day because increased emissions resulting from climate warming are countered by CO2 inhibition of isoprene emissions. The impact of biogenic isoprene emissions on the global O3 burden and CH4 lifetime is small but locally significant, and the impact of changes in isoprene emissions on atmospheric chemistry depends strongly on the state of climate and chemistry.

Extended Canadian middle atmosphere model: zonal-mean climatology and physical parameterizations

This paper describes the energetics and zonal-mean state of the upward extension of the Canadian Middle Atmosphere Model, which extends from the ground to ~210 km. The model includes realistic parameterizations of the major physical processes from the ground up to the lower thermosphere and exhibits a broad spectrum of geophysical variability. The rationale for the extended model is to examine the nature of the physical and dynamical processes in the mesosphere/lower thermosphere (MLT) region without the artificial effects of an imposed sponge layer which can modify the circulation in an unrealistic manner. The zonal-mean distributions of temperature and zonal wind are found to be in reasonable agreement with observations in most parts of the model domain below ~150 km. Analysis of the global-average energy and momentum budgets reveals a balance between solar extreme ultraviolet heating and molecular diffusion and a thermally direct viscous meridional circulation above 130 km, with the viscosity coming from molecular diffusion and ion drag. Below 70 km, radiative equilibrium prevails in the global mean. In the MLT region between ~70 and 120 km, many processes contribute to the global energy budget. At solstice, there is a thermally indirect meridional circulation driven mainly by parameterized nonorographic gravity-wave drag. This circulation provides a net global cooling of up to 25 K d^-1.

Consistent large-scale temperature responses in warm and cold climates

Climate-model simulations of the large-scale temperature responses to increased radiative forcing include enhanced land-sea contrast, stronger response at higher latitudes than in the tropics, and differential responsesin warm and cool season climates to uniform forcing. Here we show that these patterns are also characteristic of model simulations of past climates. The differences in the responses over land as opposed to over the ocean, between high and low latitudes, and between summer and winter are remarkably consistent (proportional and nearly linear) across simulations of both cold and warm climates. Similar patterns also appear in historical observations and paleoclimatic reconstructions, implying that such responses are characteristic features of the climate system, and not simple model artifacts, thereby increasing our confidence in the ability of climate models to correctly simulate different climatic states.

A high-resolution near-infrared extraterrestrial solar spectrum derived from ground-based Fourier transform spectrometer measurements

A detailed spectrally-resolved extraterrestrial solar spectrum (ESS) is important for line-by-line radiative transfer modeling in the near-infrared (near-IR). Very few observationally-based high-resolution ESS are available in this spectral region. Consequently the theoretically-calculated ESS by Kurucz has been widely adopted. We present the CAVIAR (Continuum Absorption at Visible and Infrared Wavelengths and its Atmospheric Relevance) ESS which is derived using the Langley technique applied to calibrated observations using a ground-based high-resolution Fourier transform spectrometer (FTS) in atmospheric windows from 2000–10000 cm-1 (1–5 μm). There is good agreement between the strengths and positions of solar lines between the CAVIAR and the satellite-based ACE-FTS (Atmospheric Chemistry Experiment-FTS) ESS, in the spectral region where they overlap, and good agreement with other ground-based FTS measurements in two near-IR windows. However there are significant differences in the structure between the CAVIAR ESS and spectra from semi-empirical models. In addition, we found a difference of up to 8 % in the absolute (and hence the wavelength-integrated) irradiance between the CAVIAR ESS and that of Thuillier et al., which was based on measurements from the Atmospheric Laboratory for Applications and Science satellite and other sources. In many spectral regions, this difference is significant, as the coverage factor k = 2 (or 95 % confidence limit) uncertainties in the two sets of observations do not overlap. Since the total solar irradiance is relatively well constrained, if the CAVIAR ESS is correct, then this would indicate an integrated “loss” of solar irradiance of about 30 W m-2 in the near-IR that would have to be compensated by an increase at other wavelengths.