Title: Using reduced-complexity volcanic aerosol and climate models to produce large ensemble simulations of Holocene temperature
Volcanic eruptions are one of the most important drivers of climate variability, but climate model simulations typically show stronger surface cooling than proxy-based reconstructions. Uncertainties associated with eruption source parameters, aerosol–climate modelling, and internal climate variability might explain those discrepancies, but their quantification using complex global climate models is computationally expensive. In this study, we combine a reduced-complexity volcanic aerosol model (EVA_H) and a climate model (FaIR) to simulate global-mean surface temperature from 6755 BCE to 1900 CE (8705 to 50 BP) accounting for volcanic forcing, solar irradiance, orbital, ice sheet, greenhouse gases, land-use forcing, and anthropogenic aerosols and ozone forcing for the historical period (1750–1900 CE). The negligible computational cost of the models enables us to use a Monte Carlo approach to propagate uncertainties associated with eruption source parameters, aerosol and climate modelling, and internal climate variability. Averaging over the last 9000 years, we obtain a global-mean volcanic forcing of −0.15 W m−2 and an associated surface cooling of 0.12 K. Averaged over the 14 largest eruptions (injecting more than 20 Tg of SO2) of 1250–1900 CE, the mean temperature response in tree-ring-based reconstructions is in good agreement with the our simulations, scaled to Northern Hemisphere summer temperature. For individual eruptions, discrepancies between the simulated and reconstructed surface temperature response are almost always within uncertainties. At multimillennial timescales, our simulations reproduce the Holocene global warming trend typically derived from simulations and data assimilation products but exhibit some discrepancies on centennial to millennial timescales. In particular, the Medieval Climate Anomaly to Little Ice Age transition is weaker in our simulations, and we also do not capture a relatively cool period between 3000 and 1000 BCE (5000 and 3000 BP), visible in climate reanalyses. We discuss how uncertainties in land-use forcing and model limitations might explain these differences. Our study demonstrates the value of reduced-complexity volcanic aerosol–climate models to simulate climate at annual to multimillennial timescales.  more » « less
Award ID(s):
1803946 2102993
PAR ID:
10661559
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
European Geophysical Union
Date Published:
Journal Name:
Climate of the Past
Volume:
21
Issue:
10
ISSN:
1814-9332
Page Range / eLocation ID:
1755 to 1778
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Climate variability in the last millennium (past 1000 years) is dominated by the effects of large-magnitude volcanic eruptions; however, a long-standing mismatch exists between model-simulated and tree-ring-derived surface cooling. Accounting for the self-limiting effects of large sulfur dioxide (SO2) injections and the limitations in tree-ring records, such as lagged responses due to biological memory, reconciles some of the discrepancy, but uncertainties remain, particularly for the largest tropical eruptions. The representation of volcanic forcing in the latest generation of climate models has improved significantly, but most models prescribe the aerosol optical properties rather than using SO2 emissions directly and including interactions between the aerosol, chemistry, and dynamics. Here, we use the UK Earth System Model (UKESM) to simulate the climate of the last millennium (1250–1850 CE) using volcanic SO2 emissions. Averaged across all large-magnitude eruptions, we find similar Northern Hemisphere (NH) summer cooling compared with other last-millennium climate simulations from the Paleoclimate Modelling Intercomparison Project Phase 4 (PMIP4), run with both SO2 emissions and prescribed forcing, and a continued overestimation of surface cooling compared with tree-ring reconstructions. However, for the largest-magnitude tropical eruptions in 1257 (Mt. Samalas) and 1815 (Mt. Tambora), some models, including UKESM1, suggest a smaller NH summer cooling that is in better agreement with tree-ring records. In UKESM1, we find that the simulated volcanic forcing differs considerably from the PMIP4 dataset used in models without interactive aerosol schemes, with marked differences in the hemispheric spread of the aerosol, resulting in lower forcing in the NH when SO2 emissions are used. Our results suggest that, for the largest tropical eruptions, the spatial distribution of aerosol can account for some of the discrepancies between model-simulated and tree-ring-derived cooling. Further work should therefore focus on better resolving the spatial distribution of aerosol forcing for past eruptions. 
    more » « less
  2. Abstract. The injection of sulfur into the stratosphere by volcanic eruptions is thedominant driver of natural climate variability oninterannual to multidecadal timescales. Based on a set of continuous sulfateand sulfur records from a suite of ice cores from Greenland and Antarctica,the HolVol v.1.0 database includes estimates of the magnitudes andapproximate source latitudes of major volcanic stratospheric sulfurinjection (VSSI) events for the Holocene (from 9500 BCE or 11 500 years BP to1900 CE), constituting an extension of the previous record by 7000 years.The database incorporates new-generation ice-core aerosol records with asub-annual temporal resolution and a demonstrated sub-decadal dating accuracyand precision. By tightly aligning and stacking the ice-core records on theWD2014 chronology from Antarctica, we resolve long-standing inconsistenciesin the dating of ancient volcanic eruptions that arise from biased (i.e.,dated too old) ice-core chronologies over the Holocene for Greenland. Wereconstruct 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 bipolarsulfate 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 locationsis in agreement with prior reconstructions for the past 2500 years. Intotal, these eruptions injected 7410 Tg S into thestratosphere: 70 % from tropical eruptions and 25 % from NHextratropical eruptions. A long-term latitudinally and monthly resolvedstratospheric aerosol optical depth (SAOD) time series is reconstructed fromthe HolVol VSSI estimates, representing the first Holocene-scalereconstruction constrained by Greenland and Antarctica ice cores. These newlong-term reconstructions of past VSSI and SAOD variability confirm evidencefrom regional volcanic eruption chronologies (e.g., from Iceland) in showingthat the Early Holocene (9500–7000 BCE) experienced a higher number ofvolcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared withthe past 2500 years. This increase coincides with the rapid retreat of icesheets during deglaciation, providing context for potential future increasesin volcanic activity in regions under projected glacier melting in the 21stcentury. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/PANGAEA.928646 (Sigl et al., 2021). 
    more » « less
  3. null (Ed.)
    Abstract. Simulations of climate of the last millennium (LM) show that external forcing had a major contribution to the evolution of temperatures; warmer and colder periods like the Medieval Climate Anomaly (MCA; ca. 950–1250 CE) and the Little Ice Age (LIA; ca. 1450–1850 CE) were critically influenced by changes in solar and volcanic activity. Even if this influence is mainly observed in terms of temperatures, evidence from simulations and reconstructions shows that other variables related to atmospheric dynamics and hydroclimate were also influenced by external forcing over some regions. In this work, simulations from the Coupled Model Intercomparison Project Phase 5 and Paleoclimate Modelling Intercomparison Project Phase 3 (CMIP5/PMIP3) are analyzed to explore the influence of external forcings on the dynamical and hydrological changes during the LM at different spatial and temporal scales. Principal component (PC) analysis is used to obtain the modes of variability governing the global evolution of climate and to assess their correlation with the total external forcing at multidecadal to multicentennial timescales. For shorter timescales, a composite analysis is used to address the response to specific events of external forcing like volcanic eruptions. The results show coordinated long-term changes in global circulation patterns, which suggest expansions and contractions of the Hadley cells and latitudinal displacements of westerlies in response to external forcing. For hydroclimate, spatial patterns of drier and wetter conditions in areas influenced by the North Atlantic Oscillation (NAO), Northern Annular Mode (NAM), and Southern Annular Mode (SAM) and alterations in the intensity and distribution of monsoons and convergence zones are consistently found. Similarly, a clear short-term response is found in the years following volcanic eruptions. Although external forcing has a greater influence on temperatures, the results suggest that dynamical and hydrological variations over the LM exhibit a direct response to external forcing at both long and short timescales that is highly dependent on the particular simulation and model. 
    more » « less
  4. Abstract. The 852/3 CE eruption of Mount Churchill, Alaska, was one of the largestfirst-millennium volcanic events, with a magnitude of 6.7 (VEI 6) and atephra volume of 39.4–61.9 km3 (95 % confidence). The spatial extent of the ash fallout from this event is considerable and the cryptotephra (White River Ash east; WRAe) extends as far as Finland and Poland. Proximal ecosystem and societal disturbances have been linked with this eruption; however, wider eruption impacts on climate and society are unknown. Greenland ice core records show that the eruption occurred in winter 852/3 ± 1 CE and that the eruption is associated with a relatively moderate sulfate aerosol loading but large abundances of volcanic ash and chlorine. Here we assess the potential broader impact of this eruption using palaeoenvironmental reconstructions, historical records and climate model simulations. We also use the fortuitous timing of the 852/3 CE Churchill eruption and its extensively widespread tephra deposition of the White River Ash (east) (WRAe) to examine the climatic expression of the warm Medieval Climate Anomaly period (MCA; ca. 950–1250 CE) from precisely linked peatlands in the North Atlantic region. The reconstructed climate forcing potential of the 852/3 CE Churchill eruptionis moderate compared with the eruption magnitude, but tree-ring-inferredtemperatures report a significant atmospheric cooling of 0.8 ∘Cin summer 853 CE. Modelled climate scenarios also show a cooling in 853 CE, although the average magnitude of cooling is smaller (0.3 ∘C). The simulated spatial patterns of cooling are generally similar to those generated using the tree-ring-inferred temperature reconstructions. Tree-ring-inferred cooling begins prior to the date of the eruption suggesting that natural internal climate variability may have increased the climate system's susceptibility to further cooling. The magnitude of the reconstructed cooling could also suggest that the climate forcing potential of this eruption may be underestimated, thereby highlighting the need for greater insight into, and consideration of, the role of halogens and volcanic ash when estimating eruption climate forcing potential. Precise comparisons of palaeoenvironmental records from peatlands acrossNorth America and Europe, facilitated by the presence of the WRAe isochron,reveal no consistent MCA signal. These findings contribute to the growingbody of evidence that characterises the MCA hydroclimate astime-transgressive and heterogeneous rather than a well-defined climaticperiod. The presence of the WRAe isochron also demonstrates that nolong-term (multidecadal) climatic or societal impacts from the 852/3 CEChurchill eruption were identified beyond areas proximal to the eruption.Historical evidence in Europe for subsistence crises demonstrate a degree of temporal correspondence on interannual timescales, but similar events were reported outside of the eruption period and were common in the 9thcentury. The 852/3 CE Churchill eruption exemplifies the difficulties ofidentifying and confirming volcanic impacts for a single eruption, even whenthe eruption has a small age uncertainty. 
    more » « less
  5. null (Ed.)
    Abstract. Changes in glacier length reflect the integrated response to local fluctuations in temperature and precipitation resulting from both external forcing (e.g., volcanic eruptions or anthropogenic CO2) and internal climate variability. In order to interpret the climate history reflected in the glacier moraine record, the influence of both sources of climate variability must therefore be considered. Here we study the last millennium of glacier-length variability across the globe using a simple dynamic glacier model, which we force with temperature and precipitation time series from a 13-member ensemble of simulations from a global climate model. The ensemble allows us to quantify the contributions to glacier-length variability from external forcing (given by the ensemble mean) and internal variability (given by the ensemble spread). Within this framework, we find that internal variability is the predominant source of length fluctuations for glaciers with a shorter response time (less than a few decades). However, for glaciers with longer response timescales (more than a few decades) external forcing has a greater influence than internal variability. We further find that external forcing also dominates when the response of glaciers from widely separated regions is averaged. Single-forcing simulations indicate that, for this climate model, most of the forced response over the last millennium, pre-anthropogenic warming, has been driven by global-scale temperature change associated with volcanic aerosols. 
    more » « less