The NCAR Community Earth System Model (CESM) now includes an atmospheric component that extends in altitude to the lower thermosphere. This atmospheric model, known as the Whole Atmosphere Community Climate Model (WACCM), includes fully interactive chemistry, allowing, for example, a self-consistent representation of the development and recovery of the stratospheric ozone hole and its effect on the troposphere. This paper focuses on analysis of an ensemble of transient simulations using CESM1(WACCM), covering the period from the preindustrial era to present day, conducted as part of phase 5 of the Coupled Model Intercomparison Project. Variability in the stratosphere, such as that associated with stratospheric sudden warmings and the development of the ozone hole, is in good agreement with observations. The signals of these phenomena propagate into the troposphere, influencing near-surface winds, precipitation rates, and the extent of sea ice. In comparison of tropospheric climate change predictions with those from a version of CESM that does not fully resolve the stratosphere, the global-mean temperature trends are indistinguishable. However, systematic differences do exist in other climate variables, particularly in the extratropics. The magnitude of the difference can be as large as the climate change response itself. This indicates that the representation of stratospheretroposphere coupling could be a major source of uncertainty in climate change projections in CESM.

[1] We have used the Whole Atmosphere Community Climate Model to produce a small (three-member) ensemble of simulations of the period 1950-2003. Comparison of model results against available observations shows that for the most part, the model is able to reproduce well the observed trends in zonal mean temperature and ozone, both as regards their magnitude and their distribution in latitude and altitude. Calculated trends in water vapor, on the other hand, are not at all consistent with observations from either the HALOE satellite instrument or the Boulder, Colorado, hygrosonde data set. We show that such lack of agreement is actually to be expected because water vapor has various sources of low-frequency variability (heating due to volcanic eruptions, the quasi-biennial oscillation and El Niño-Southern Oscillation) that can confound the determination of secular trends. The simulations also reveal the presence of other interesting behavior, such as the lack of any significant temperature trend near the mesopause, a decrease in the stratospheric age of air, and the rare occurrence of an extremely disturbed Southern Hemisphere winter.

[1] 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 . 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 (Cl y ) 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 Cl y , 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 D223081 of 29 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. Citation: Eyring, V., et al. (2006), Assessment of temperature, trace species, and ozone in chemistry-climate model simulations of the recent past,

[1] The Model for Ozone and Related Chemical Tracers, version 3 (MOZART-3), which represents the chemical and physical processes from the troposphere through the lower mesosphere, was used to evaluate the representation of long-lived tracers and ozone using three different meteorological fields. The meteorological fields are based on (1) the Whole Atmosphere Community Climate Model, version 1b (WACCM1b), (2) the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis, and (3) a new reanalysis for year 2000 from ECMWF called EXP471. Model-derived tracers (methane, water vapor, and total inorganic nitrogen) and ozone are compared to data climatologies from satellites. Model mean age of air was also derived and compared to in situ CO 2 and SF 6 data. A detailed analysis of the chemical fields simulated by shows that even though the general features characterizing the three dynamical sets are rather similar, slight differences in winds and temperature can produce substantial differences in the calculated distributions of chemical tracers. The MOZART-3 simulations that use meteorological fields from WACCM1b and ECMWF EXP471 represented best the distribution of long-lived tracers and mean age of air in the stratosphere. There was a significant improvement using the ECMWF EXP471 reanalysis data product over the ECMWF operational data product. The effect of the quasi-biennial oscillation circulation on long-lived tracers and ozone is examined.

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 GHG‐induced 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.

Abstract. This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850CMIP6 preindustrial control, historical ( -2014, and future (2015-2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunderminimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models.For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical onePublished by Copernicus Publications on behalf of the European Geosciences Union. A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m −2 . The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m −2 . In the 200-400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %).We compare the climatic effects of the CMIP6 solar forcing dataset to its CMIP5 predecessor by using timeslice experiments of two chemistry-climate models and a reference radiative transfer model. The differences in the long-term mean SSI in the CMIP6 dataset, compared to CMIP5, impact on climatological stratospheric conditions (lower shortwave heating rates of −0.35 K day −1 at the stratopause), cooler stratospheric temperatures (−1.5 K in the upper stratosphere), lower ozone abundances in the lower stratosphere (−3 %), and higher ozone abundances (+1.5 % in the upper stratosphere and lower mesosphere). Between the maximum and minimum phases of the 11-year solar cycle, there is an increase in shortwave heating rates (+0.2 K day −1 at the stratopause), temperatures (∼ 1 K at the stratopause), and ozone (+2.5 % in the upper stratosphere) in the tropical upper stratosphere using the CMIP6 forcing dataset. This solar-cycle response is slightly larger, but not statistically significantly different from that for the CMIP5 forcing dataset.CMIP6 models wi...

[1] Ozone changes and associated climate impacts in the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations are analyzed over the historical and future (2006-2100) period under four Representative Concentration Pathways (RCP). In contrast to CMIP3, where half of the models prescribed constant stratospheric ozone, CMIP5 models all consider past ozone depletion and future ozone recovery. Multimodel mean climatologies and long-term changes in total and tropospheric column ozone calculated from CMIP5 models with either interactive or prescribed ozone are in reasonable agreement with observations. However, some large deviations from observations exist for individual models with interactive chemistry, and these models are excluded in the projections. Stratospheric ozone projections forced with a single halogen, but four greenhouse gas (GHG) scenarios show largest differences in the northern midlatitudes and in the Arctic in spring (~20 and 40 Dobson units (DU) by 2100, respectively). By 2050, these differences are much smaller and negligible over Antarctica in austral spring. Differences in future tropospheric column ozone are mainly caused by differences in methane concentrations and stratospheric input, leading to~10 DU increases compared to 2000 in RCP 8.5. Large variations in stratospheric ozone particularly in CMIP5 models with interactive chemistry drive correspondingly large variations in lower stratospheric temperature trends. The results also illustrate that future Southern Hemisphere summertime circulation changes are controlled by both the ozone recovery rate and the rate of GHG increases, emphasizing the importance of simulating and taking into account ozone forcings when examining future climate projections.

Accurate representation of global stratospheric aerosols from volcanic and nonvolcanic sulfur emissions is key to understanding the cooling effects and ozone losses that may be linked to volcanic activity. Attribution of climate variability to volcanic activity is of particular interest in relation to the post-2000 slowing in the rate of global average temperature increases. We have compiled a database of volcanic SO 2 emissions and plume altitudes for eruptions from 1990 to 2014 and developed a new prognostic capability for simulating stratospheric sulfate aerosols in the Community Earth System Model. We used these combined with other nonvolcanic emissions of sulfur sources to reconstruct global aerosol properties from 1990 to 2014. Our calculations show remarkable agreement with ground-based lidar observations of stratospheric aerosol optical depth (SAOD) and with in situ measurements of stratospheric aerosol surface area density (SAD). These properties are key parameters in calculating the radiative and chemical effects of stratospheric aerosols. Our SAOD calculations represent a clear improvement over available satellite-based analyses, which generally ignore aerosol extinction below 15 km, a region that can contain the vast majority of stratospheric aerosol extinction at middle and high latitudes. Our SAD calculations greatly improve on that provided for the Chemistry-Climate Model Initiative, which misses about 60% of the SAD measured in situ on average during both volcanically active and volcanically quiescent periods.