Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

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: 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.