skip to main content


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
NSF-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. 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
  3. null (Ed.)
    Abstract We use theNorthern Hemisphere Tree-RingNetwork Development (NTREND) tree-ring database to examine the effects of using a small, highly-sensitive proxy network for paleotemperature data assimilation over the last millennium. We first evaluate our methods using pseudo-proxy experiments. These indicate that spatial assimilations using this network are skillful in the extratropical Northern Hemisphere and improve on previous NTREND reconstructions based on Point-by-Point regression. We also find our method is sensitive to climate model biases when the number of sites becomes small. Based on these experiments, we then assimilate the real NTREND network. To quantify model prior uncertainty, we produce 10 separate reconstructions, each assimilating a different climate model. These reconstructions are most dissimilar prior to 1100 CE, when the network becomes sparse, but show greater consistency as the network grows. Temporal variability is also underestimated before 1100 CE. Our assimilation method produces spatial uncertainty estimates and these identify treeline North America and eastern Siberia as regions that would most benefit from development of new millennial-length temperature-sensitive tree-ring records. We compare our multi-model mean reconstruction to five existing paleo-temperature products to examine the range of reconstructed responses to radiative forcing. We find substantial differences in the spatial patterns and magnitudes of reconstructed responses to volcanic eruptions and in the transition between the Medieval epoch and Little Ice Age. These extant uncertainties call for the development of a paleoclimate reconstruction intercomparison framework for systematically examining the consequences of proxy network composition and reconstruction methodology and for continued expansion of tree-ring proxy networks. 
    more » « less
  4. Abstract

    The effects of volcanic eruptions on hurricane statistics are examined using two long simulations from the Community Earth System Model (CESM) Last Millennium Ensemble (LME). The first is an unforced control simulation, wherein all boundary conditions were held constant at their 850 CE values (LMEcontrol). The second is a “fully forced” simulation with time evolving radiative changes from volcanic, solar, and land use changes from 850 CE through present (LMEforced). Large tropical volcanic eruptions produce the greatest change in radiative forcing during this time period, which comprise the focus of this study. The Weather Research and Forecasting (WRF) model is used to dynamically downscale 150 control years of LMEcontroland an additional 84 years of LMEforcedfor all mid-latitude volcanic eruptions between 1100 and 1850 CE. This time period was selected based on computational considerations. For each eruption, 2 years are dynamically downscaled. 23 of these volcanic eruptions are in the Northern Hemisphere and 19 are in the Southern Hemisphere. The effectiveness of the downscaling methodology is examined by applying the same downscaling approach to historical ERA-I reanalysis data and comparing the downscaled storm tracks and intensities to the International Best Track Archive for Climate Stewardship (IBTrACS) database. Hurricane statistics are then computed from both the downscaled control and downscaled forced LME simulations. Results suggest moderate effects on hurricanes from the average of all northern hemisphere eruptions, with the largest effects being from the volcanoes with the most aerosol forcing. More specifically, reductions in hurricane frequency, intensity, and lifetime following northern hemisphere eruptions are apparent. Strong evidence is also shown for correlation between eruption strength and changes in these diagnostics. The aggregate effect from both northern and southern hemisphere eruptions is minor. While reductions in frequency, intensity, and lifetime from northern hemisphere eruptions occur, the opposite effect is observed from southern hemisphere eruptions.

     
    more » « less
  5. The effects of volcanic eruptions on hurricane statistics are examined using two long simulations from the Community Earth System Model (CESM) Last Millennium Ensemble (LME). The first is an unforced control simulation, wherein all boundary conditions were held constant at their 850 CE values (LMEcontrol). The second is a “fully forced” simulation with time evolving radiative changes from volcanic, solar, and land use changes from 850 CE through present (LMEforced). Large tropical volcanic eruptions produce the greatest change in radiative forcing during this time period, which comprise the focus of this study. The Weather Research and Forecasting (WRF) model is used to dynamically downscale 150 control years of LMEcontrol and an additional 84 years of LMEforced for all mid-latitude volcanic eruptions between 1100 and 1850 CE. This time period was selected based on computational considerations. For each eruption, 2 years are dynamically downscaled. 23 of these volcanic eruptions are in the Northern Hemisphere and 19 are in the Southern Hemisphere. The effectiveness of the downscaling methodology is examined by applying the same downscaling approach to historical ERA-I reanalysis data and comparing the downscaled storm tracks and intensities to the International Best Track Archive for Climate Stewardship (IBTrACS) database. Hurricane statistics are then computed from both the downscaled control and downscaled forced LME simulations. Results suggest moderate effects on hurricanes from the average of all northern hemisphere eruptions, with the largest effects being from the volcanoes with the most aerosol forcing. More specifically, reductions in hurricane frequency, intensity, and lifetime following northern hemisphere eruptions are apparent. Strong evidence is also shown for correlation between eruption strength and changes in these diagnostics. The aggregate effect from both northern and southern hemisphere eruptions is minor. While reductions in frequency, intensity, and lifetime from northern hemisphere eruptions occur, the opposite effect is observed from southern hemisphere eruptions. 
    more » « less