ALBERT

Description: The impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry–climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry–Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1–5hPa) decreases by more than a factor of 2, from around 5 to 2%, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55K, or around 80% of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper.

Description: The radiative effects induced by the zonally asymmetric part of the ozone field have been shown to significantly change the temperature of the NH winter polar cap, and correspondingly the strength of the polar vortex. In this paper, we aim to understand the physical processes behind these effects using the National Center for Atmospheric Research (NCAR)'s Whole Atmosphere Community Climate Model, run with 1960s ozone-depleting substances and greenhouse gases. We find a mid-winter polar vortex influence only when considering the quasi-biennial oscillation (QBO) phases separately, since ozone waves affect the vortex in an opposite manner. Specifically, the emergence of a midlatitude QBO signal is delayed by 1–2 months when radiative ozone-wave effects are removed. The influence of ozone waves on the winter polar vortex, via their modulation of shortwave heating, is not obvious, given that shortwave heating is largest during fall, when planetary stratospheric waves are weakest. Using a novel diagnostic of wave 1 temperature amplitude tendencies and a synoptic analysis of upward planetary wave pulses, we are able to show the chain of events that lead from a direct radiative effect on weak early fall upward-propagating planetary waves to a winter polar vortex modulation. We show that an important stage of this amplification is the modulation of individual wave life cycles, which accumulate during fall and early winter, before being amplified by wave–mean flow feedbacks. We find that the evolution of these early winter upward planetary wave pulses and their induced stratospheric zonal mean flow deceleration is qualitatively different between QBO phases, providing a new mechanistic view of the extratropical QBO signal. We further show how these differences result in opposite radiative ozone-wave effects between east and west QBOs.

Description: To understand solar cycle signals on the Earth's surface and identify the physical mechanisms responsible, surface temperature variations from observations as well as climate model data are analyzed to characterize their spatial structure. The solar signal in the annual mean surface temperature is characterized by (i) mid-latitude warming and (ii) no warming in the tropics. The mid-latitude warming during solar maxima in both hemispheres is associated with a downward penetration of zonal mean zonal wind anomalies from the upper stratosphere during late winter. During Northern Hemisphere winter this is manifested in a modulation of the polar-night jet whereas in the Southern Hemisphere the subtropical jet plays the major role. Warming signals are particularly apparent over the Eurasian continent and ocean frontal zones, including a previously reported lagged response over the North Atlantic. In the tropics, local warming occurs over the Indian and central Pacific oceans during high solar activity. However, this warming is counter balanced by cooling over the cold tongue sectors in the southeastern Pacific and the South Atlantic, and results in a very weak zonally averaged tropical mean signal. The cooling in the ocean basins is associated with stronger cross-equatorial winds resulting from a northward shift of the ascending branch of the Hadley circulation during solar maxima. To understand the complex processes involved in the solar signal transfer, results of an idealized middle atmosphere–ocean coupled model experiment on the impact of stratospheric zonal wind changes are compared with solar signals in observations. The model results suggest that both tropical and extra-tropical solar surface signals can result from circulation changes in the upper stratosphere through (i) a downward migration of wave–zonal mean flow interactions and (ii) changes in the stratospheric mean meridional circulation. These experiments support earlier evidence of an indirect solar influence from the stratosphere.

Description: Changes in incoming solar ultraviolet radiation over the 11-year solar cycle affect stratospheric ozone abundances. It is important to quantify the magnitude, structure, and seasonality of the associated solar-ozone response (SOR) to understand the impact of the 11-year solar cycle on climate. Part 1 of this two-part study uses multiple linear regression analysis to extract the SOR in a number of recently updated satellite ozone datasets covering different periods within the epoch 1970 to 2013. The annual mean SOR in the updated version 7.0 (v7.0) Stratospheric Aerosol and Gas Experiment (SAGE) II number density dataset (1984–2004) is very consistent with that found in the previous v6.2. In contrast, we find a substantial decrease in the magnitude of the SOR in the tropical upper stratosphere in the SAGE II v7.0 mixing ratio dataset (∼ 1 %) compared to the v6.2 (∼ 4 %). This difference is shown to be largely attributable to the change in the independent stratospheric temperature dataset used to convert SAGE II ozone number densities to mixing ratios. Since these temperature records contain substantial uncertainties, we suggest that datasets based on SAGE II number densities are currently most reliable for evaluating the SOR. We further analyse three extended ozone datasets that combine SAGE II v7.0 number densities with more recent GOMOS (Global Ozone Monitoring by Occultation of Stars) or OSIRIS (Optical Spectrograph and Infrared Imager System) measurements. The extended SAGE–OSIRIS dataset (1984–2013) shows a smaller and less statistically significant SOR across much of the tropical upper stratosphere compared to the SAGE II data alone. In contrast, the two SAGE–GOMOS datasets (1984–2011) show SORs that are in closer agreement with the original SAGE II data and therefore appear to provide a more reliable estimate of the SOR. We also analyse the SOR in the recent Solar Backscatter Ultraviolet Instrument (SBUV) Merged Ozone Dataset (SBUVMOD) version 8.6 (VN8.6) (1970–2012) and SBUV Merged Cohesive VN8.6 (1978–2012) datasets and compare them to the previous SBUVMOD VN8.0 (1970–2009). Over their full lengths, the three records generally agree in terms of the broad magnitude and structure of the annual mean SOR. The main difference is that SBUVMOD VN8.6 shows a smaller and less significant SOR in the tropical upper stratosphere and therefore more closely resembles the SAGE II v7.0 mixing ratio data than does the SBUV Merged Cohesive VN8.6, which has a more continuous SOR of ∼ 2 % in this region. The sparse spatial and temporal sampling of limb satellite instruments prohibits the extraction of sub-annual variations in the SOR from SAGE-based datasets. However, the SBUVMOD VN8.6 dataset suggests substantial month-to-month variations in the SOR, particularly in the winter extratropics, which may be important for the proposed high-latitude dynamical response to the solar cycle. Overall, the results highlight substantial uncertainties in the magnitude and structure of the observed SOR from different satellite records. The implications of these uncertainties for understanding and modelling the effects of solar variability on climate should be explored.

Description: This study aims to quantify how much of the extratropical Tropopause Inversion Layer (TIL) comes from the modulation by planetary and synoptic-scale waves. By analyzing high-resolution observations, it also puts other TIL enhancing mechanisms into context. Using gridded COSMIC GPS-RO temperature profiles from 2007–2013 we are able to extract the extratropical wave signal by a simplified wavenumber-frequency domain filtering method, and to quantify the resulting TIL enhancement. By subtracting the extratropical wave signal, we show how much of the TIL is associated with other processes, at mid and high latitudes, for both Hemispheres and all seasons. The instantaneous modulation by planetary and synoptic-scale waves is almost entirely responsible for the TIL in mid-latitudes. This means that wave-mean flow interactions, inertia-gravity waves or the residual circulation are of minor importance in mid-latitudes. At polar regions, the extratropical wave modulation is dominant for the TIL strength as well, but there is also a clear fingerprint from sudden stratospheric warmings (SSWs) and final warmings in both hemispheres. Therefore, polar vortex breakups are partially responsible for the observed polar TIL strength in winter (if SSWs occur) and spring. Also, part of the polar summer TIL strength cannot be explained by extratropical wave modulation. After many modelling studies that proposed different TIL enhancing mechanisms in the last decade, our study finally identifies which processes dominate the extratropical TIL strength and their relative contribution, by analyzing observations only. It remains to be determined, however, which roles the different planetary and synoptic-scale wave types play within the total extratropical wave modulation of the TIL; and what causes the observed amplification of extratropical waves near the tropopause.

Description: 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 (1850–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 Maunder-minimum-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 one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). 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 time-slice 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 with a well-resolved shortwave radiation scheme are encouraged to prescribe SSI changes and include solar-induced stratospheric ozone variations, in order to better represent solar climate variability compared to models that only prescribe TSI and/or exclude the solar-ozone response. We show that monthly-mean solar-induced ozone variations are implicitly included in the SPARC/CCMI CMIP6 Ozone Database for historical simulations, which is derived from transient chemistry–climate model simulations and has been developed for climate models that do not calculate ozone interactively. CMIP6 models without chemistry that perform a preindustrial control simulation with time-varying solar forcing will need to use a modified version of the SPARC/CCMI Ozone Database that includes solar variability. CMIP6 models with interactive chemistry are also encouraged to use the particle forcing datasets, which will allow the potential long-term effects of particles to be addressed for the first time. The consideration of particle forcing has been shown to significantly improve the representation of reactive nitrogen and ozone variability in the polar middle atmosphere, eventually resulting in further improvements in the representation of solar climate variability in global models.

Description: It is well established that variable wintertime planetary wave forcing in the stratosphere controls the variability of Arctic stratospheric ozone through changes in the strength of the polar vortex and the residual circulation. While previous studies focused on the variations in upward wave flux entering the lower stratosphere, here the impact of downward planetary wave reflection on ozone is investigated for the first time. Utilizing the MERRA2 reanalysis and a fully coupled chemistry–climate simulation with the Community Earth System Model (CESM1(WACCM)) of the National Center for Atmospheric Research (NCAR), we find two downward wave reflection effects on ozone: (1) the direct effect in which the residual circulation is weakened during winter, reducing the typical increase of ozone due to upward planetary wave events and (2) the indirect effect in which the modification of polar temperature during winter affects the amount of ozone destruction in spring. Winter seasons dominated by downward wave reflection events (i.e., reflective winters) are characterized by lower Arctic ozone concentration, while seasons dominated by increased upward wave events (i.e., absorptive winters) are characterized by relatively higher ozone concentration. This behavior is consistent with the cumulative effects of downward and upward planetary wave events on polar stratospheric ozone via the residual circulation and the polar temperature in winter. The results establish a new perspective on dynamical processes controlling stratospheric ozone variability in the Arctic by highlighting the key role of wave reflection.

Description: The recently observed variability in the tropical tropopause layer (TTL), which features a warming of 0.9 K over the past decade (2001–2011), is investigated with a number of sensitivity experiments from simulations with NCAR's CESM-WACCM chemistry–climate model. The experiments have been designed to specifically quantify the contributions from natural as well as anthropogenic factors, such as solar variability (Solar), sea surface temperatures (SSTs), the quasi-biennial oscillation (QBO), stratospheric aerosols (Aerosol), greenhouse gases (GHGs) and the dependence on the vertical resolution in the model. The results show that, in the TTL from 2001 through 2011, a cooling in tropical SSTs leads to a weakening of tropical upwelling around the tropical tropopause and hence relative downwelling and adiabatic warming of 0.3 K decade-1; stronger QBO westerlies result in a 0.2 K decade-1 warming; increasing aerosols in the lower stratosphere lead to a 0.2 K decade-1 warming; a prolonged solar minimum contributes about 0.2 K decade-1 to a cooling; and increased GHGs have no significant influence. Considering all the factors mentioned above, we compute a net 0.5 K decade-1 warming, which is less than the observed 0.9 K decade-1 warming over the past decade in the TTL. Two simulations with different vertical resolution show that, with higher vertical resolution, an extra 0.8 K decade-1 warming can be simulated through the last decade compared with results from the "standard" low vertical resolution simulation. Model results indicate that the recent warming in the TTL is partly caused by stratospheric aerosols and mainly due to internal variability, i.e. the QBO and tropical SSTs. The vertical resolution can also strongly influence the TTL temperature response in addition to variability in the QBO and SSTs.

Description: The tropical tropopause layer (TTL) acts as a transition layer between the troposphere and the stratosphere over several kilometers, where air has both tropospheric and stratospheric properties. Within this region, a fine-scale feature is located: the tropopause inversion layer (TIL), which consists of a sharp temperature inversion at the tropopause and the corresponding high static stability values right above, which theoretically affect the dispersion relations of atmospheric waves like Rossby or inertia–gravity waves and hamper stratosphere–troposphere exchange (STE). Therefore, the TIL receives increasing attention from the scientific community, mainly in the extratropics so far. Our goal is to give a detailed picture of the properties, variability and forcings of the tropical TIL, with special emphasis on small-scale equatorial waves and the quasi-biennial oscillation (QBO). We use high-resolution temperature profiles from the COSMIC satellite mission, i.e., ∼ 2000 measurements per day globally, between 2007 and 2013, to derive TIL properties and to study the fine-scale structures of static stability in the tropics. The situation at near tropopause level is described by the 100 hPa horizontal wind divergence fields, and the vertical structure of the QBO is provided by the equatorial winds at all levels, both from the ERA-Interim reanalysis. We describe a new feature of the equatorial static stability profile: a secondary stability maximum below the zero wind line within the easterly QBO wind regime at about 20–25 km altitude, which is forced by the descending westerly QBO phase and gives a double-TIL-like structure. In the lowermost stratosphere, the TIL is stronger with westerly winds. We provide the first evidence of a relationship between the tropical TIL strength and near-tropopause divergence, with stronger (weaker) TIL with near-tropopause divergent (convergent) flow, a relationship analogous to that of TIL strength with relative vorticity in the extratropics. To elucidate possible enhancing mechanisms of the tropical TIL, we quantify the signature of the different equatorial waves on the vertical structure of static stability in the tropics. All waves show, on average, maximum cold anomalies at the thermal tropopause, warm anomalies above and a net TIL enhancement close to the tropopause. The main drivers are Kelvin, inertia–gravity and Rossby waves. We suggest that a similar wave modulation will exist at mid- and polar latitudes from the extratropical wave modes

Description: Detection and attribution (D&A) simulations were important components of CMIP5 and underpinned the climate change detection and attribution assessments of the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. The primary goals of the Detection and Attribution Model Intercomparison Project (DAMIP) are to facilitate improved estimation of the contributions of anthropogenic and natural forcing changes to observed global warming as well as to observed global and regional changes in other climate variables; to contribute to the estimation of how historical emissions have altered and are altering contemporary climate risk; and to facilitate improved observationally constrained projections of future climate change. D&A studies typically require unforced control simulations and historical simulations including all major anthropogenic and natural forcings. Such simulations will be carried out as part of the DECK and the CMIP6 historical simulation. In addition D&A studies require simulations covering the historical period driven by individual forcings or subsets of forcings only: such simulations are proposed here. Key novel features of the experimental design presented here include firstly new historical simulations with aerosols-only, stratospheric-ozone-only, CO2-only, solar-only, and volcanic-only forcing, facilitating an improved estimation of the climate response to individual forcing, secondly future single forcing experiments, allowing observationally constrained projections of future climate change, and thirdly an experimental design which allows models with and without coupled atmospheric chemistry to be compared on an equal footing.