Search

Here we present a study of the 11 yr sunspot cycle's imprint on the Northern Hemisphere atmospheric circulation, using three recently developed gridded upper-air data sets that extend back to the early twentieth century. We find a robust response of the tropospheric late-wintertime circulation to the sunspot cycle, independent from the data set. This response is particularly significant over Europe, although results show that it is not directly related to a North Atlantic Oscillation (NAO) modulation; instead, it reveals a significant connection to the more meridional Eurasian pattern (EU). The magnitude of mean seasonal temperature changes over the European land areas locally exceeds 1 K in the lower troposphere over a sunspot cycle. We also analyse surface data to address the question whether the solar signal over Europe is temporally stable for a longer 250 yr period. The results increase our confidence in the existence of an influence of the 11 yr cycle on the European climate, but the signal is much weaker in the first half of the period compared to the second half. The last solar minimum (2005 to 2010), which was not included in our analysis, shows anomalies that are consistent with our statistical results for earlier solar minima.

The chemistry-climate model (CCM) SOCOL has been used to evaluate the contributions of the main anthropogenic factors to the simulated changes of ozone and stratospheric dynamics during the 21st century. As the main anthropogenic factors we consider the atmospheric concentration of the greenhouse gases (GHG), ozone depleting substances (ODS) and sea surface temperature and sea ice (SST/SI). The latter is considered here as an independent factor because the majority of the CCMs prescribe its evolution. We have performed three sets of "time slice" numerical experiments for the years 2000, 2050, and 2100 taking into account all factors separately and all together. The total column ozone increase during the first half of the 21st century is caused by the ODS, especially in the middle and high latitudes of both hemispheres. In the tropics and the extra tropical region of the Northern Hemisphere (NH) the SST/SI forcing plays a very important role in the evolution of atmospheric ozone during the second half of the 21st century. The GHG affect the temperature and ozone mainly in the upper stratosphere and in the lower stratosphere of the high latitudes of the Southern Hemisphere (SH). In the lower tropical stratosphere of the NH, the long-term changes of the temperature, zonal wind and the meridional circulation are controlled mainly by the SST/SI. The strong contribution of the SST/SI to the ozone and circulation changes in the future implies that some differences between the results by different CCMs could be caused by the applied SST/SI rather than by the CCM's deficiencies. We suggest taking this issue into account for the planning of the future model evaluation campaigns.

Ozone depletion is caused by the anthropogenic increase of halogen-containing species in the atmosphere, which results in the enhancement of the concentration of reactive chlorine and bromine in the stratosphere. To reduce the influence of anthropogenic ozone-depleting substances (ODS), the Montreal Protocol was agreed by Governments in 1987, with several Amendments and Adjustments adopted later. In order to assess the benefits of the Montreal Protocol and its Amendments and Adjustments (MPA) on ozone and UV radiation, two different runs of the chemistry-climate model (CCM) SOCOL have been carried out. The first run was driven by the emission of ozone depleting substances (ODS) prescribed according to the restrictions of the MPA. For the second run we allow the ODS to grow by 3% annually. We find that the MPA would have saved up to 80% of the global annual total ozone by the end of the 21st century. Our calculations also show substantial changes of the stratospheric circulation pattern as well as in surface temperature and precipitations that could occur in the world without MPA implementations. To illustrate the changes in UV radiation at the surface and to emphasise certain features, which can only be seen for some particular regions if the influence of the cloud cover changes is accounted for, we calculate geographical distribution of the erythemally weighted irradiance (Eery). For the no Montreal Protocol simulation Eery increases by factor of 4 to 16 between the 1970s and 2100. For the scenario including the Montreal Protocol it is found that UV radiation starts to decrease in 2000, with continuous decline of 5% to 10% at middle latitudes in the both Northern and Southern Hemispheres.

The lack of long and reliable time series of solar spectral irradiance (SSI) measurements makes an accurate quantification of solar contributions to recent climate change difficult. Whereas earlier SSI observations and models provided a qualitatively consistent picture of the SSI variability, recent measurements by the SORCE (SOlar Radiation and Climate Experiment) satellite suggest a significantly stronger variability in the ultraviolet (UV) spectral range and changes in the visible and near-infrared (NIR) bands in anti-phase with the solar cycle. A number of recent chemistry-climate model (CCM) simulations have shown that this might have significant implications on the Earth's atmosphere. Motivated by these results, we summarize here our current knowledge of SSI variability and its impact on Earth's climate. We present a detailed overview of existing SSI measurements and provide thorough comparison of models available to date. SSI changes influence the Earth's atmosphere, both directly, through changes in shortwave (SW) heating and therefore, temperature and ozone distributions in the stratosphere, and indirectly, through dynamical feedbacks. We investigate these direct and indirect effects using several state-of-the art CCM simulations forced with measured and modelled SSI changes. A unique asset of this study is the use of a common comprehensive approach for an issue that is usually addressed separately by different communities. We show that the SORCE measurements are difficult to reconcile with earlier observations and with SSI models. Of the five SSI models discussed here, specifically NRLSSI (Naval Research Laboratory Solar Spectral Irradiance), SATIRE-S (Spectral And Total Irradiance REconstructions for the Satellite era), COSI (COde for Solar Irradiance), SRPM (Solar Radiation Physical Modelling), and OAR (Osservatorio Astronomico di Roma), only one shows a behaviour of the UV and visible irradiance qualitatively resembling that of the recent SORCE measurements. However, the integral of the SSI computed with this model over the entire spectral range does not reproduce the measured cyclical changes of the total solar irradiance, which is an essential requisite for realistic evaluations of solar effects on the Earth's climate in CCMs. We show that within the range provided by the recent SSI observations and semi-empirical models discussed here, the NRLSSI model and SORCE observations represent the lower and upper limits in the magnitude of the SSI solar cycle variation. The results of the CCM simulations, forced with the SSI solar cycle variations estimated from the NRLSSI model and from SORCE measurements, show that the direct solar response in the stratosphere is larger for the SORCE than for the NRLSSI data. Correspondingly, larger UV forcing also leads to a larger surface response. Finally, we discuss the reliability of the available data and we propose additional coordinated work, first to build composite SSI data sets out of scattered observations and to refine current SSI models, and second, to run coordinated CCM experiments.

Through the 21st century, anthropogenic emissions of the greenhouse gases N2O and CH4 are projected to increase, thus increasing their atmospheric concentrations. Consequently, reactive nitrogen species produced from N2O and reactive hydrogen species produced from CH4 are expected to play an increasingly important role in determining stratospheric ozone concentrations. Eight chemistry-climate model simulations were performed to assess the sensitivity of stratospheric ozone to different emissions scenarios for N2O and CH4. Global-mean total column ozone increases through the 21st century in all eight simulations as a result of CO2-induced stratospheric cooling and decreasing stratospheric halogen concentrations. Larger N2O concentrations were associated with smaller ozone increases, due to reactive nitrogen-mediated ozone destruction. In the simulation with the largest N2O increase, global-mean total column ozone increased by 4.3 DU through the 21st century, compared with 10.0 DU in the simulation with the smallest N2O increase. In contrast, larger CH4 concentrations were associated with larger ozone increases; global-mean total column ozone increased by 16.7 DU through the 21st century in the simulation with the largest CH4 concentrations and by 4.4 DU in the simulation with the lowest CH4 concentrations. CH4 leads to ozone loss in the upper and lower stratosphere by increasing the rate of reactive hydrogen-mediated ozone loss cycles, however in the lower stratosphere and troposphere, CH4 leads to ozone increases due to photochemical smog-type chemistry. In addition to this mechanism, total column ozone increases due to H2O-induced cooling of the stratosphere, and slowing of the chlorine-catalyzed ozone loss cycles due to an increased rate of the CH4 + Cl reaction. Stratospheric column ozone through the 21st century exhibits a near-linear response to changes in N2O and CH4 surface concentrations, which provides a simple parameterization for the ozone response to changes in these gases.

We have modeled the atmospheric impact of a major solar energetic particle event similar in intensity to what is thought of the Carrington Event of 1–2 September 1859. Ionization rates for the August 1972 solar proton event, which had an energy spectrum comparable to the Carrington Event, were scaled up in proportion to the fluence estimated for both events. We have assumed such an event to take place in the year 2020 in order to investigate the impact on the modern, near future atmosphere. Effects on atmospheric chemistry, temperature and dynamics were investigated using the 3-D Chemistry Climate Model SOCOL v2.0. We find significant responses of NOx, HOx, ozone, temperature and zonal wind. Ozone and NOx have in common an unusually strong and long-lived response to this solar proton event. The model suggests a 3-fold increase of NOx generated in the upper stratosphere lasting until the end of November, and an up to 10-fold increase in upper mesospheric HOx. Due to the NOx and HOx enhancements, ozone reduces by up to 60–80% in the mesosphere during the days after the event, and by up to 20–40% in the middle stratosphere lasting for several months after the event. Total ozone is reduced by up to 20 DU in the Northern Hemisphere and up to 10 DU in the Southern Hemisphere. Free tropospheric and surface air temperatures show a significant cooling of more than 3 K and zonal winds change significantly by 3–5 m s−1 in the UTLS region. In conclusion, a solar proton event, if it took place in the near future with an intensity similar to that ascribed to of the Carrington Event of 1859, must be expected to have a major impact on atmospheric composition throughout the middle atmosphere, resulting in significant and persistent decrease in total ozone.

The mesospheric hydroxyl radical (OH) is mainly produced by the water vapor (H2O) photolysis and could be considered as a proxy for the influence of the solar irradiance variability on the mesosphere. We analyze the tropical mean response of the mesospheric OH and H2O data as observed by the Aura Microwave Limb Sounder (MLS) to 27-day solar variability. The analysis is performed for two time periods corresponding to the different phases of the 11-yr cycle: from December 2004 to December 2005 (the period of "high activity" with a pronounced 27-day solar cycle) and from August 2008 to August 2009 ("solar minimum" period with a vague 27-day solar cycle). We demonstrate, for the first time, that in the mesosphere the daily time series of OH concentrations correlate well with the solar irradiance (correlation coefficients up to 0.79) at zero time-lag. At the same time H2O anticorrelates (correlation coefficients up to −0.74) with the solar irradiance at non-zero time-lag. We found that the response of OH and H2O to the 27-day variability of the solar irradiance is strong for the period of the high solar activity and negligible for the solar minimum conditions. It allows us to suggest that the 27-day cycle in the solar irradiance and in OH and H2O are physically connected.

Time series of total column abundances of hydrogen chloride (HCl), chlorine nitrate (ClONO2), and hydrogen fluoride (HF) were determined from ground-based Fourier transform infrared (FTIR) spectra recorded at 17 sites belonging to the Network for the Detection of Atmospheric Composition Change (NDACC) and located between 80.05° N and 77.82° S. By providing such a near-global overview on ground-based measurements of the two major stratospheric chlorine reservoir species, HCl and ClONO2, the present study is able to confirm the decrease of the atmospheric inorganic chlorine abundance during the last few years. This decrease is expected following the 1987 Montreal Protocol and its amendments and adjustments, where restrictions and a subsequent phase-out of the prominent anthropogenic chlorine source gases (solvents, chlorofluorocarbons) were agreed upon to enable a stabilisation and recovery of the stratospheric ozone layer. The atmospheric fluorine content is expected to be influenced by the Montreal Protocol, too, because most of the banned anthropogenic gases also represent important fluorine sources. But many of the substitutes to the banned gases also contain fluorine so that the HF total column abundance is expected to have continued to increase during the last few years. The measurements are compared with calculations from five different models: the two-dimensional Bremen model, the two chemistry-transport models KASIMA and SLIMCAT, and the two chemistry-climate models EMAC and SOCOL. Thereby, the ability of the models to reproduce the absolute total column amounts, the seasonal cycles, and the temporal evolution found in the FTIR measurements is investigated and inter-compared. This is especially interesting because the models have different architectures. The overall agreement between the measurements and models for the total column abundances and the seasonal cycles is good. Linear trends of HCl, ClONO2, and HF are calculated from both measurement and model time series data, with a focus on the time range 2000–2009. This period is chosen because from most of the measurement sites taking part in this study, data are available during these years. The precision of the trends is estimated with the bootstrap resampling method. The sensitivity of the trend results with respect to the fitting function, the time of year chosen and time series length is investigated, as well as a bias due to the irregular sampling of the measurements. The measurements and model results investigated here agree qualitatively on a decrease of the chlorine species by around 1% yr−1. The models simulate an increase of HF of around 1% yr−1. This also agrees well with most of the measurements, but some of the FTIR series in the Northern Hemisphere show a stabilisation or even a decrease in the last few years. In general, for all three gases, the measured trends vary more strongly with latitude and hemisphere than the modelled trends. Relative to the FTIR measurements, the models tend to underestimate the decreasing chlorine trends and to overestimate the fluorine increase in the Northern Hemisphere. At most sites, the models simulate a stronger decrease of ClONO2 than of HCl. In the FTIR measurements, this difference between the trends of HCl and ClONO2 depends strongly on latitude, especially in the Northern Hemisphere.

We have compared composition changes of NO, NO2, H2O2, O3, N2O, HNO3, N2O5, HNO4, ClO, HOCl, and ClONO2 as observed by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat in the aftermath of the "Halloween" solar proton event (SPE) in late October 2003 at 25–0.01 hPa in the Northern Hemisphere (40–90° N) and simulations performed by the following atmospheric models: the Bremen 2-D model (B2dM) and Bremen 3-D Chemical Transport Model (B3dCTM), the Central Aerological Observatory (CAO) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, the modeling tool for SOlar Climate Ozone Links studies (SOCOL and SOCOLi), and the Whole Atmosphere Community Climate Model (WACCM4). The large number of participating models allowed for an evaluation of the overall ability of atmospheric models to reproduce observed atmospheric perturbations generated by SPEs, particularly with respect to NOy and ozone changes. We have further assessed the meteorological conditions and their implications for the chemical response to the SPE in both the models and observations by comparing temperature and tracer (CH4 and CO) fields. Simulated SPE-induced ozone losses agree on average within 5 % with the observations. Simulated NOy enhancements around 1 hPa, however, are typically 30 % higher than indicated by the observations which are likely to be related to deficiencies in the used ionization rates, though other error sources related to the models' atmospheric background state and/or transport schemes cannot be excluded. The analysis of the observed and modeled NOy partitioning in the aftermath of the SPE has demonstrated the need to implement additional ion chemistry (HNO3 formation via ion-ion recombination and water cluster ions) into the chemical schemes. An overestimation of observed H2O2 enhancements by all models hints at an underestimation of the OH/HO2 ratio in the upper polar stratosphere during the SPE. The analysis of chlorine species perturbations has shown that the encountered differences between models and observations, particularly the underestimation of observed ClONO2 enhancements, are related to a smaller availability of ClO in the polar night region already before the SPE. In general, the intercomparison has demonstrated that differences in the meteorology and/or initial state of the atmosphere in the simulations cause a relevant variability of the model results, even on a short timescale of only a few days.

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

This study investigates the influence of the Galactic Cosmic Rays (GCRs) on the atmospheric composition, temperature and dynamics by means of the 3-D Chemistry Climate Model (CCM) SOCOL v2.0. Ionization rates were parameterized according to CRAC:CRII (Cosmic Ray induced Cascade: Application for Cosmic Ray Induced Ionization), a detailed state-of-the-art model describing the effects of GCRs in the entire altitude range of the CCM from 0–80 km. We find statistically significant effects of GCRs on tropospheric and stratospheric NOx, HOx, ozone, temperature and zonal wind, whereas NOx, HOx and ozone are annually averaged and the temperature and the zonal wind are monthly averaged. In the Southern Hemisphere, the model suggests the GCR-induced NOx increase to exceed 10 % in the tropopause region (peaking with 20 % at the pole), whereas HOx is showing a decrease of about 3 % caused by enhanced conversion into HNO3. As a consequence, ozone is increasing by up to 3 % in the relatively unpolluted southern troposphere, where its production is sensitive to additional NOx from GCRs. Conversely, in the northern polar lower stratosphere, GCRs are found to decrease O3 by up to 3 %, caused by the additional heterogeneous chlorine activation via ClONO2 + HCl following GCR-induced production of ClONO2. There is an apparent GCR-induced acceleration of the zonal wind of up to 5 m s−1 in the Northern Hemisphere below 40 km in February, and a deceleration at higher altitudes with peak values of 3 m s−1 around 70 km altitude. The model also indentifies GCR-induced changes in the surface air, with warming in the eastern part of Europe and in Russia (up to 2.25 K for March values) and cooling in Siberia and Greenland (by almost 2 K). We show that these surface temperature changes develop even when the GCR-induced ionization is taken into account only above 18 km, suggesting that the stratospherically driven strengthening of the polar night jet extends all the way down to the Earth's surface.

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

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

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

Stratospheric ozone recovery in the Southern Hemisphere is expected to drive pronounced trends in atmospheric temperature and circulation from the stratosphere to the troposphere in the 21st century; therefore ozone changes need to be accounted for in future climate simulations. Many climate models do not have interactive ozone chemistry and rely on prescribed ozone fields, which may be obtained from coupled chemistry-climate model (CCM) simulations. However CCMs vary widely in their predictions of ozone evolution, complicating the selection of ozone boundary conditions for future climate simulations. In order to assess which models might be expected to better simulate future ozone evolution, and thus provide more realistic ozone boundary conditions, we assess the ability of twelve CCMs to simulate observed ozone climatology and trends and rank the models according to their errors averaged across the individual diagnostics chosen. According to our analysis no one model performs better than the others in all the diagnostics; however, combining errors in individual diagnostics into one metric of model performance allows us to objectively rank the models. The multi-model average shows better overall agreement with the observations than any individual model. Based on this analysis we conclude that the multi-model average ozone projection presents the best estimate of future ozone evolution and recommend it for use as a boundary condition in future climate simulations. Our results also demonstrate a sensitivity of the analysis to the choice of reference data set for vertical ozone distribution over the Antarctic, highlighting the constraints that large observational uncertainty imposes on such model verification.

The connection between the El Niño 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.

While chemistry-climate models are able to reproduce many characteristics of the global total column ozone field and its long-term evolution, they have fared less well in simulating the commonly used diagnostic of the area of the Antarctic ozone hole i.e. the area within the 220 Dobson Unit (DU) contour. Two possible reasons for this are: (1) the underlying Global Climate Model (GCM) does not correctly simulate the size of the polar vortex, and (2) the stratospheric chemistry scheme incorporated into the GCM, and/or the model dynamics, results in systematic biases in the total column ozone fields such that the 220 DU contour is no longer appropriate for delineating the edge of the ozone hole. Both causes are examined here with a view to developing ozone hole area diagnostics that better suit measurement-model inter-comparisons. The interplay between the shape of the meridional mixing barrier at the edge of the vortex and the meridional gradients in total column ozone across the vortex edge is investigated in measurements and in 5 chemistry-climate models (CCMs). Analysis of the simulation of the polar vortex in the CCMs shows that the first of the two possible causes does play a role in some models. This in turn affects the ability of the models to simulate the large observed meridional gradients in total column ozone. The second of the two causes also strongly affects the ability of the CCMs to track the observed size of the ozone hole. It is shown that by applying a common algorithm to the CCMs for selecting a delineating threshold unique to each model, a more appropriate diagnostic of ozone hole area can be generated that shows better agreement with that derived from observations.

The representation of the Tropical Tropopause Layer (TTL) in 13 different Chemistry Climate Models (CCMs) designed to represent the stratosphere is analyzed. Simulations for 1960–2005 and 1980–2100 are analyzed. Simulations for 1960–2005 are compared to reanalysis model output. CCMs are able to reproduce the basic structure of the TTL. There is a large (10 K) spread in annual mean tropical cold point tropopause temperatures. CCMs are able to reproduce historical trends in tropopause pressure obtained from reanalysis products. Simulated historical trends in cold point tropopause temperatures are not consistent across models or reanalyses. The pressure of both the tropical tropopause and the level of main convective outflow appear to have decreased (increased altitude) in historical runs as well as in reanalyses. Decreasing pressure trends in the tropical tropopause and level of main convective outflow are also seen in the future. Models consistently predict decreasing tropopause and convective outflow pressure, by several hPa/decade. Tropical cold point temperatures are projected to increase by 0.09 K/decade. Tropopause anomalies are highly correlated with tropical surface temperature anomalies and with tropopause level ozone anomalies, less so with stratospheric temperature anomalies. Simulated stratospheric water vapor at 90 hPa increases by up to 0.5–1 ppmv by 2100. The result is consistent with the simulated increase in temperature, highlighting the correlation of tropopause temperatures with stratospheric water vapor.

We have estimated changes in surface solar ultraviolet (UV) radiation under cloud free conditions in the 21st century based on simulations of 11 coupled Chemistry-Climate Models (CCMs). The total ozone columns and vertical profiles of ozone and temperature projected from CCMs were used as input to a radiative transfer model in order to calculate the corresponding erythemal irradiance levels. Time series of monthly erythemal irradiance received at the surface during local noon are presented for the period 1960 to 2100. Starting from the first decade of the 21st century, the surface erythemal irradiance decreases globally as a result of the projected stratospheric ozone recovery at rates that are larger in the first half of the 21st century and smaller towards its end. This decreasing tendency varies with latitude, being more pronounced over areas where stratospheric ozone has been depleted the most after 1980. Between 2000 and 2100 surface erythemal irradiance is projected to decrease over midlatitudes by 5 to 15%, while at the southern high latitudes the decrease is twice as much. In this study we have not included effects from changes in cloudiness, surface reflectivity and tropospheric aerosol loading, which will likely be affected in the future due to climate change. Consequently, over some areas the actual changes in future UV radiation may be different depending on the evolution of these parameters.

Interannual-to-decadal variability in stratospheric ozone and climate have a number of common sources, such as variations in solar irradiance, stratospheric aerosol loading due to volcanic eruptions, El Niño Southern Oscillation variability and the quasi-biennial oscillation (QBO). Currently available data records as well as model simulations addressing stratospheric chemical climate variability mostly cover only the past few decades, which is often insufficient to address natural interannual-to-decadal variability. Here we make use of recently reconstructed and re-evaluated data products to force and validate transient ensemble model simulations (nine members) across the twentieth century computed by means of the chemistry-climate model SOCOL (SOlar Climate Ozone Links). The forcings include sea surface temperatures, sea ice, solar irradiance, stratospheric aerosols, QBO, changes in land properties, greenhouse gases, ozone depleting substances, and emissions of carbon monoxides, and nitrogen oxides. The transient simulations are in good agreement with observations, reconstructions and reanalyses and allow quantification of interannual-to-decadal variability during the 20th century. All ensemble members are able to capture the low-frequency variability in tropical and mid-latitude total ozone as well as in the strength of the subtropical jet, suggesting a realistic response to external forcings in this area. The region of the northern polar vortex exhibits a large internal variability that is found in the frequency, seasonality, and strength of major warmings as well as in the strength of the modeled polar vortex. Results from process-oriented analysis, such as correlation between the vertical Eliassen Palm flux (EP flux) component and polar variables as well as stratospheric ozone trends, are of comparable magnitude to those observed and are consistent in all analysed ensemble members. Yet, trend estimates of the vertical EP flux component vary greatly among ensemble members precluding any robust conclusions. This suggests that internal variability in models must be accounted for in order to quantify the atmospheric model response in wave energy upon external forcings.

We describe version 2.0 of the chemistry-climate model (CCM) SOCOL. The new version includes fundamental changes of the transport scheme such as transporting all chemical species of the model individually and applying a family-based correction scheme for mass conservation for species of the nitrogen, chlorine and bromine groups, a revised transport scheme for ozone, furthermore more detailed halogen reaction and deposition schemes, and a new cirrus parameterisation in the tropical tropopause region. By means of these changes the model manages to overcome or considerably reduce deficiencies recently identified in SOCOL version 1.1 within the CCM Validation activity of SPARC (CCMVal). In particular, as a consequence of these changes, regional mass loss or accumulation artificially caused by the semi-Lagrangian transport scheme can be significantly reduced, leading to much more realistic distributions of the modelled chemical species, most notably of the halogens and ozone.

The LYRA instrument onboard ESA PROBA2 satellite will provide 6-hourly solar irradiance at the Lyman-alpha (121.6 nm) and the Herzberg continuum (~200–220 nm wavelength range). Because the nowcasting of the neutral and ionic state of the middle atmosphere requires the solar irradiance for the wide spectral range (120–680 nm) we have developed the statistical tool for the reconstruction of the full spectrum from the LYRA measurements. The accuracy of the reconstructed irradiance has been evaluated with 1-D transient radiative-convective model with neutral and ion chemistry using the daily solar spectral irradiance measured with SUSIM and SOLSTICE instruments onboard UARS satellite. We compared the results of transient 1-year long model simulations for 2000 driven by the observed and reconstructed solar irradiance and showed that the reconstruction of the full spectrum using linear regression equation based on the solar irradiance in two LYRA channels can be successfully used for nowcasting of the middle atmosphere. We have also identified conditions when the proposed approach does not yield spectral reconstruction with sufficient accuracy.

A pronounced ENSO cycle occurred from 1986 to 1989, accompanied by distinct dynamical and chemical anomalies in the global troposphere and stratosphere. Reproducing these effects with current climate models not only provides a model test but also contributes to our still limited understanding of ENSO's effect on stratosphere-troposphere coupling. We performed several sets of ensemble simulations with a chemical climate model (SOCOL) forced with global sea surface temperatures. Results were compared with observations and with large-ensemble simulations performed with an atmospheric general circulation model (MRF9). We focus our analysis on the extratropical stratosphere and its coupling with the troposphere. In this context, the circulation over the North Atlantic sector is particularly important. Relative to the La Niña winter 1989, observations for the El Niño winter 1987 show a negative North Atlantic Oscillation index with corresponding changes in temperature and precipitation patterns, a weak polar vortex, a warm Arctic middle stratosphere, negative and positive total ozone anomalies in the tropics and at middle to high latitudes, respectively, as well as anomalous upward and poleward Eliassen-Palm (EP) flux in the midlatitude lower stratosphere. Most of the tropospheric features are well reproduced in the ensemble means in both models, though the amplitudes are underestimated. In the stratosphere, the SOCOL simulations compare well with observations with respect to zonal wind, temperature, EP flux, meridional mass streamfunction, and ozone, but magnitudes are underestimated in the middle stratosphere. With respect to the mechanisms relating ENSO to stratospheric circulation, the results suggest that both, upward and poleward components of anomalous EP flux are important for obtaining the stratospheric signal and that an increase in strength of the Brewer-Dobson circulation is part of that signal.

In this paper we document 'SOCOL', a new chemistry-climate model, which has been ported for regular PCs and shows good wall-clock performance. An extensive validation of the model results against present-day climate data obtained from observations and assimilation data sets shows that the model describes the climatological state of the atmosphere for the late 1990s with reasonable accuracy. The model has a significant temperature bias only in the upper stratosphere and near the tropopause at high latitudes. The latter is the result of the rather low vertical resolution of the model near the tropopause. The former can be attributed to a crude representation of radiation heating in the middle atmosphere. A comparison of the simulated and observed link between the tropical stratospheric structure and the strength of the polar vortex shows that in general, both observations and simulations reveal a higher temperature and ozone mixing ratio in the lower tropical stratosphere for the case with stronger Polar night jet (PNJ) and slower Brewer-Dobson circulation as predicted by theoretical studies.

In recent years a number of chemistry-climate models have been developed with an emphasis on the stratosphere. Such models cover a wide range of time scales of integration and vary considerably in complexity. The results of specific diagnostics are here analysed to examine the differences amongst individual models and observations, to assess the consistency of model predictions, with a particular focus on polar ozone. For example, many models indicate a significant cold bias in high latitudes, the "cold pole problem', particularly in the southern hemisphere during winter and spring. This is related to wave propagation from the troposphere which can be improved by improving model horizontal resolution and with the use of non-orographic gravity wave drag. As a result of the widely differing modelled polar temperatures, different amounts of polar stratospheric clouds are simulated which in turn result in varying ozone values in the models. The results are also compared to determine the possible future behaviour of ozone, with an emphasis on the polar regions and mid-latitudes. All models predict eventual ozone recovery, but give a range of results concerning its timing and extent. Differences in the simulation of gravity waves and planetary waves as well as model resolution are likely major sources of uncertainty for this issue. In the Antarctic, the ozone hole has probably reached almost its deepest although the vertical and horizontal extent of depletion may increase slightly further over the next few years. According to the model results, Antarctic ozone recovery could begin any year within the range 2001 to 2008. The limited number of models which have been integrated sufficiently far indicate that full recovery of ozone to 1980 levels may not occur in the Antarctic until about the year 2050. For the Arctic, most models indicate that small ozone losses may continue for a few more years and that recovery could begin any year within the range 2004 to 2019. The start of ozone recovery in the Arctic is therefore expected to appear later than in the Antarctic. Further, interannual variability will tend to mask the signal for longer than in the Antarctic, delaying still further the date at which ozone recovery may be said to have started. Because of this inherent variability of the system, the decadal evolution of Arctic ozone will not necessarily be a direct response to external forcing.

We have applied chemistry-climate model (CCM) SOCOL to simulate the distribution of the temperature and gas species in the upper stratosphere and mesosphere. As an input for the simulation, we employ daily spectral solar UV irradiance measured by SUSIM instrument onboard UARS satellite in January 1992. We have carried out an ensemble of nine 1-month long simulations using slightly different initial states of the atmosphere. We have compared the obtained time evolution of the simulated species and temperature with available satellite measurements. The obtained results allowed us to define the areas where the nowcast and short-term forecast of the atmospheric species with CCM SOCOL could be successful.