The injection of sulfur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulfate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The database incorporates improvements to the ice core records (in terms of synchronisation and dating) and refinements to the methods used to estimate VSSI from ice core records, and it includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815), and Laki (1783), are within 10 % of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682, and 1108 CE. The long-term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr⁻¹, ∼ 50 % greater than a prior reconstruction due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33 %) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–1900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates including a constant background injection of stratospheric sulfur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2.

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

The enhancement of the stratospheric aerosol layer by volcanic eruptions induces a complex set of responses causing global and regional climate effects on a broad range of timescales. Uncertainties exist regarding the climatic response to strong volcanic forcing identified in coupled climate simulations that contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). In order to better understand the sources of these model diversities, the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) has defined a coordinated set of idealized volcanic perturbation experiments to be carried out in alignment with the CMIP6 protocol. VolMIP provides a common stratospheric aerosol data set for each experiment to minimize differences in the applied volcanic forcing. It defines a set of initial conditions to assess how internal climate variability contributes to determining the response. VolMIP will assess to what extent volcanically forced responses of the coupled ocean–atmosphere system are robustly simulated by state-of-the-art coupled climate models and identify the causes that limit robust simulated behavior, especially differences in the treatment of physical processes. This paper illustrates the design of the idealized volcanic perturbation experiments in the VolMIP protocol and describes the common aerosol forcing input data sets to be used.

The injection of sulfur into the stratosphere by volcanic eruptions is the dominant driver of natural climate variability on interannual to multidecadal timescales. Based on a set of continuous sulfate and sulfur records from a suite of ice cores from Greenland and Antarctica, the HolVol v.1.0 database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events for the Holocene (from 9500 BCE or 11 500 years BP to 1900 CE), constituting an extension of the previous record by 7000 years. The database incorporates new-generation ice-core aerosol records with a sub-annual temporal resolution and a demonstrated sub-decadal dating accuracy and precision. By tightly aligning and stacking the ice-core records on the WD2014 chronology from Antarctica, we resolve long-standing inconsistencies in the dating of ancient volcanic eruptions that arise from biased (i.e., dated too old) ice-core chronologies over the Holocene for Greenland. We reconstruct a total of 850 volcanic eruptions with injections in excess of 1 teragram of sulfur (Tg S); of these eruptions, 329 (39 %) are located in the low latitudes with bipolar sulfate deposition, 426 (50 %) are located in the Northern Hemisphere extratropics (NHET) and 88 (10 %) are located in the Southern Hemisphere extratropics (SHET). The spatial distribution of the reconstructed eruption locations is in agreement with prior reconstructions for the past 2500 years. In total, these eruptions injected 7410 Tg S into the stratosphere: 70 % from tropical eruptions and 25 % from NH extratropical eruptions. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the HolVol VSSI estimates, representing the first Holocene-scale reconstruction constrained by Greenland and Antarctica ice cores. These new long-term reconstructions of past VSSI and SAOD variability confirm evidence from regional volcanic eruption chronologies (e.g., from Iceland) in showing that the Early Holocene (9500–7000 BCE) experienced a higher number of volcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared with the past 2500 years. This increase coincides with the rapid retreat of ice sheets during deglaciation, providing context for potential future increases in volcanic activity in regions under projected glacier melting in the 21st century. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/PANGAEA.928646 (Sigl et al., 2021).

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.

Volcanic eruptions have long been studied for their wide range of climatic effects. Although global-scale climatic impacts following the formation of stratospheric sulfate aerosol are well understood, many aspects of the evolution of the early volcanic aerosol cloud and regional impacts are uncertain. In the last twenty years, several advances have been made, mainly due to improved satellite measurements and observations enabling the effects of small-magnitude eruptions to be quantified, new proxy reconstructions used to investigate the impact of past eruptions, and state-of-the-art aerosol-climate modelling that has led to new insights on how volcanic eruptions affect the climate. Looking to the future, knowledge gaps include the role of co-emissions in volcanic plumes, the impact of eruptions on tropical hydroclimate and Northern Hemisphere winter climate, and the role of eruptions in long-term climate change. Future model development, dedicated model intercomparison projects, interdisciplinary collaborations, and the application of advanced statistical techniques will facilitate more complex and detailed studies. Ensuring that the next large-magnitude explosive eruption is well observed will be critical in providing invaluable observations that will bridge remaining gaps in our understanding.