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
Last-millennium volcanic forcing and climate response using SO 2 emissions
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
- PAR ID:
- 10618809
- Publisher / Repository:
- Copernicus Publications on behalf of the European Geosciences Union.
- Date Published:
- Journal Name:
- Climate of the Past
- Volume:
- 21
- Issue:
- 1
- ISSN:
- 1814-9332
- Page Range / eLocation ID:
- 161 to 184
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
The 540s, 1450s, and 1600s represent three of the five coldest decades in the Common Era (CE). In each of these cases, the cause of these cold pulses has been attributed to large volcanic eruptions. However, the provenance of the eruption and magnitude of the volcanic forcing remains uncertain. Here, we use high-resolution sulfur isotopes in Greenland and Antarctic ice cores measured across these events to provide a means of improving sulfur loading estimates for these eruptions. In each case, the largest reconstructed tree-ring cooling is associated with an extratropical eruption, and the high-altitude stratospheric sulfate loading of these events is substantially smaller than previous estimates (by up to a factor of two). These results suggest an increased sensitivity of the reconstructed Northern Hemisphere summer temperature response to extratropical eruptions. This highlights the importance of climate feedbacks and processes that amplify and prolong the cooling signal from high latitudes, such as changes in sea ice extent and ocean heat content.more » « less
-
Abstract We critically reexamine the question of whether volcanic eruptions cause surface warming over Eurasia in winter, in the light of recent modeling studies that have suggested internal variability may overwhelm any forced volcanic response, even for the very largest eruptions during the Common Era. Focusing on the last millennium, we combine model output, instrumental observations, tree-ring records, and ice cores to build a new temperature reconstruction that specifically targets the boreal winter season. We focus on 20 eruptions over the last millennium with volcanic stratospheric sulfur injections (VSSIs) larger than the 1991 Pinatubo eruption. We find that only 7 of these 20 large events are followed by warm surface temperature anomalies over Eurasia in the first posteruption winter. Examining the 13 events that show cold posteruption anomalies, we find no correlation between the amplitude of winter cooling and VSSI mass. We also find no evidence that the North Atlantic Oscillation is correlated with VSSI in winter, a key element of the proposed mechanism through which large, low-latitude eruptions might cause winter warming over Eurasia. Furthermore, by inspecting individual eruptions rather than combining events into a superposed epoch analysis, we are able to reconcile our findings with those of previous studies. Analysis of two additional paleoclimatic datasets corroborates the lack of posteruption Eurasian winter warming. Our findings, covering the entire last millennium, confirm the findings of most recent modeling studies and offer important new evidence that large, low-latitude eruptions are not, in general, followed by significant surface wintertime warming over Eurasia.more » « less
-
null (Ed.)Abstract Scientific understanding of low-frequency tropical Pacific variability, especially responses to perturbations in radiative forcing, suffers from short observational records, sparse proxy networks, and bias in model simulations. Here, we combine the strengths of proxies and models through coral-based paleoclimate data assimilation. We combine coral archives ( δ 18 O, Sr/Ca) with the dynamics, spatial teleconnections, and intervariable relationships of the CMIP5/PMIP3 Past1000 experiments using the Last Millennium Reanalysis data assimilation framework. This analysis creates skillful reconstructions of tropical Pacific temperatures over the observational era. However, during the period of intense volcanism in the early nineteenth century, southwestern Pacific corals produce El Niño–Southern Oscillation (ENSO) reconstructions that are of opposite sign from those from eastern Pacific corals and tree ring records. We systematically evaluate the source of this discrepancy using 1) single-proxy experiments, 2) varied proxy system models (PSMs), and 3) diverse covariance patterns from the Past1000 simulations. We find that individual proxy records and coral PSMs do not significantly contribute to the discrepancy. However, following major eruptions, the southwestern Pacific corals locally record more persistent cold anomalies than found in the Past1000 experiments and canonical ENSO teleconnections to the southwest Pacific strongly control the reconstruction response. Furthermore, using covariance patterns independent of ENSO yields reconstructions consistent with coral archives across the Pacific. These results show that model bias can strongly affect how proxy information is processed in paleoclimate data assimilation. As we illustrate here, model bias influences the magnitude and persistence of the response of the tropical Pacific to volcanic eruptions.more » « less
An official website of the United States government

