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.
- NSF-PAR ID:
- 10219218
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 118
- Issue:
- 12
- ISSN:
- 0027-8424
- Page Range / eLocation ID:
- e2019145118
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
The oceanic response to recent tropical eruptions is examined in Large Ensemble (LE) experiments from two fully coupled global climate models, the Community Earth System Model (CESM) and the Geophysical Fluid Dynamics Laboratory Earth System Model (ESM2M), each forced by a distinct volcanic forcing dataset. Following the simulated eruptions of Agung, El Chichón, and Pinatubo, the ocean loses heat and gains oxygen and carbon, in general agreement with available observations. In both models, substantial global surface cooling is accompanied by El Niño–like equatorial Pacific surface warming a year after the volcanic forcing peaks. A mechanistic analysis of the CESM and ESM2M responses to Pinatubo identifies remote wind forcing from the western Pacific as a major driver of this El Niño–like response. Following eruption, faster cooling over the Maritime Continent than adjacent oceans suppresses convection and leads to persistent westerly wind anomalies over the western tropical Pacific. These wind anomalies excite equatorial downwelling Kelvin waves and the upwelling of warm subsurface anomalies in the eastern Pacific, promoting the development of El Niño conditions through Bjerknes feedbacks a year after eruption. This El Niño–like response drives further ocean heat loss through enhanced equatorial cloud albedo, and dominates global carbon uptake as upwelling of carbon-rich waters is suppressed in the tropical Pacific. Oxygen uptake occurs primarily at high latitudes, where surface cooling intensifies the ventilation of subtropical thermocline waters. These volcanically forced ocean responses are large enough to contribute to the observed decadal variability in oceanic heat, carbon, and oxygen.
-
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
-
more » « less
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
Abstract Earth's hydrological cycle is expected to intensify in response to global warming, with a “wet‐gets‐wetter, dry‐gets‐drier” response anticipated over the ocean. Subtropical regions (∼15°–30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data‐modeling approach to reconstruct global and zonal‐mean rainfall patterns during the early Eocene (∼56–48 million years ago). The Deep‐Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid‐ (30°–60°N/S) and high‐latitudes (>60°N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0°–15°N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter‐Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation‐evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter‐model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy‐derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation‐induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2019145118
