skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium
Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.  more » « less
Award ID(s):
1743738 1702439 1923957
PAR ID:
10450224
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Paleoceanography and Paleoclimatology
Volume:
36
Issue:
4
ISSN:
2572-4517
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Large tropical volcanic eruptions can affect the climate of many regions on Earth, yet it is uncertain how the largest eruptions over the past millennium may have altered Earth’s hydroclimate. Here, we analyze the global hydroclimatic response to all the tropical volcanic eruptions over the past millennium that were larger than the Mount Pinatubo eruption of 1991. Using the Paleo Hydrodynamics Data Assimilation product (PHYDA), we find that these large volcanic eruptions tended to produce dry conditions over tropical Africa, Central Asia and the Middle East and wet conditions over much of Oceania and the South American monsoon region. These anomalies are statistically significant, and they persisted for more than a decade in some regions. The persistence of the anomalies is associated with southward shifts in the Intertropical Convergence Zone and sea surface temperature changes in the Pacific and Atlantic oceans. We compare the PHYDA results with the stand-alone model response of the Community Earth System Model (CESM)-Last Millennium Ensemble. We find that the proxy-constrained PHYDA estimates are larger and more persistent than the responses simulated by CESM. Understanding which of these estimates is more realistic is critical for accurately characterizing the hydroclimate risks of future volcanic eruptions. 
    more » « less
  2. Abstract The North American Southwest (NASW) and South American Southwest (SASW) are regions susceptible to prolonged and intense droughts that can span a decade or more (i.e., megadroughts). Although the drivers and impacts of megadroughts in each region and their co-occurrence have been examined in paleoclimate reconstructions, it is not known whether climate models simulate co-occurring megadroughts in these regions with characteristics and drivers that are similar to the real world. We compare the temporal characteristics of concurrent megadroughts and the Pacific Ocean conditions associated with these events in the Paleo Hydrodynamics Data Assimilation (PHYDA) product and the Community Earth System Model Last Millennium Ensemble (CESM-LME). We find that concurrent megadroughts in PHYDA and CESM-LME have similar temporal characteristics, but the relationship between hydroclimate conditions in the NASW and SASW is different between proxy-based estimates and the climate model. Further analyses reveal that changes in the tropical Pacific Ocean are weaker during concurrent megadroughts in the CESM-LME compared to those in PHYDA and that their teleconnection patterns and strengths are different. Reconstruction methodology is also found to be a factor in how the relationship between the tropical Pacific and each region is characterized. These results together indicate that while the CESM-LME simulates concurrent megadroughts with temporal characteristics similar to PHYDA, it does so for different reasons; this result leaves open the question of whether climate models used for future projections can accurately capture the risk of concurrent megadroughts in future projections. 
    more » « less
  3. Abstract. The Medieval Climate Anomaly (MCA; ca. 950–1250 CE) and the Little Ice Age (LIA; ca. 1450–1850 CE) were periods generally characterized by respectively higher and lower temperatures in many regions. However, they have also been associated with drier and wetter conditions in areas around the Intertropical Convergence Zone (ITCZ) and the Asian Monsoon region and in areas impacted by large-scale climatic modes like the Northern Annular Mode and Southern Annular Mode (NAM and SAM respectively). To analyze coordinated changes in large-scale hydroclimate patterns and whether similar changes also extend to other periods of the Last Millennium (LM) outside the MCA and the LIA, reconstruction-based products have been analyzed. This includes the collection of tree-ring-based drought atlases (DAs), the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR). These analyses have shown coherent changes in the hydroclimate of tropical and extratropical regions, such as northern and central South America, East Africa, western North America, western Europe, the Middle East, Southeast Asia, and the Indo-Pacific, during the MCA, the LIA and other periods of the LM. Comparisons with model simulations from the Community Earth System Model – Last Millennium Ensemble (CESM-LME) and phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6) show that both external forcing and internal variability contributed to these changes, with the contribution of internal variability being particularly important in the Indo-Pacific basin and that of external forcing in the Atlantic basin. These results may help to identify not only those areas showing coordinated changes, but also those regions more impacted by the internal variability, where forced model simulations would not be expected to successfully reproduce the evolution of past actual hydroclimate changes. 
    more » « less
  4. 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
  5. Abstract The tropical response to explosive volcanism remains underconstrained in the paleoclimate record. While the atmosphere cools due to aerosol forcing following volcanic eruptions, modeling evidence suggests that the tropical Pacific responds with compensatory warming. Given the rarity of large volcanic eruptions and the short instrumental record, these modeling results require independent verification. Here, we test for links between volcanism and tropical Pacific dynamics using the newly developed Paleo Hydrodynamics Data Assimilation product (PHYDA), which spans the past 2,000 years. Using Pacific sea surface temperature fields from PHYDA and coincident volcanic eruptions, we test the response of the El Niño–Southern Oscillation (ENSO) to large, tropical volcanic eruptions. We identify a weak El Niño‐like response of the tropical Pacific in the year following sufficiently large, tropical volcanic eruptions. While the response is not significant at the 95% confidence level using superposed epoch analysis (SEA) and self‐organizing maps, a significant result does emerge when employing probability density functions. Our results indicate that the widely used SEA approach, based on composite averaging, may not be sufficiently sensitive to capture an ENSO response in the presence of large internal variability. We additionally conclude that inconsistencies in both the spatial patterns and magnitudes between climate models and PHYDA results indicate that current models overestimate the regional tropical response to volcanic forcing. 
    more » « less