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.
Despite the importance of interdecadal climate variability, we have a limited understanding of which geographic regions are associated with global temperature variability at these timescales. The instrumental record tends to be too short to develop sample statistics to study interdecadal climate variability, and Coupled Model Intercomparison Project, Phase 5 (CMIP5) climate models tend to disagree about which locations most strongly influence global mean interdecadal temperature variability. Here we use a new paleoclimate data assimilation product, the Last Millennium Reanalysis (LMR), to examine where local variability is associated with global mean temperature variability at interdecadal timescales. The LMR framework uses an ensemble Kalman filter data assimilation approach to combine the latest paleoclimate data and state‐of‐the‐art model data to generate annually resolved field reconstructions of surface temperature, which allow us to explore the timing and dynamics of preinstrumental climate variability in new ways. The LMR consistently shows that the middle‐ to high‐latitude north Pacific and the high‐latitude North Atlantic tend to lead global temperature variability on interdecadal timescales. These findings have important implications for understanding the dynamics of low‐frequency climate variability in the preindustrial era.
more » « less- Award ID(s):
- 1702423
- NSF-PAR ID:
- 10445730
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Atmospheres
- Volume:
- 124
- Issue:
- 17-18
- ISSN:
- 2169-897X
- Page Range / eLocation ID:
- p. 9905-9917
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We use online data assimilation to combine information from a linear inverse model of coupled atmosphere‐ocean dynamics with proxy records to create a new annual‐resolution reconstruction of atmosphere and ocean fields over the last millennium. Instrumental validation of reconstructed sea‐surface temperature and 0–700 m ocean heat content shows broad regions of positive spatial correlations, and high correlations (∼0.6–0.9) for global averages and indices of large‐scale modes of atmospheric variability. Compared to previous reconstructions, the online reconstructions show global and hemispheric averages with little‐to‐no millennial‐scale trend and global‐mean temperatures ∼0.25–0.5 K cooler during early periods (1000–1400 C.E.). The spatial anomaly differences of average temperature between an early (1000–1250 C.E.) and later (1400–1700 C.E.) period show warm anomalies over high‐latitude Europe and cool tropical conditions in partial agreement with previous assessments. The addition of online data assimilation, which provides dynamical memory to climate proxy information, is shown to be crucial for adequately characterizing decadal‐to‐centennial‐scale variability of 0–700 m ocean heat content. Furthermore, the climate forecasts provide model‐based physical constraints for atmosphere–ocean interaction, which become increasingly important during early periods when less proxy information is available for assimilation.
-
null (Ed.)Global hydroclimatic changes from 1950 to 2018 are analyzed using updated data of land precipitation, streamfow, and an improved form of the Palmer Drought Severity Index. The historical changes are then compared with climate model-simulated response to external forcing to determine how much of the recent change is forced response. It is found that precipitation has increased from 1950 to 2018 over mid-high latitude Eurasia, most North America, Southeast South America, and Northwest Australia, while it has decreased over most Africa, eastern Australia, the Mediterranean region, the Middle East, and parts of East Asia, central South America, and the Pacifc coasts of Canada. Streamfow records largely confrm these precipitation changes. The wetting trend over Northwest Australia and Southeast South America is most pronounced in austral summer while the drying over Africa and wetting trend over mid-high latitude Eurasia are seen in all seasons. Coupled with the drying caused by rising surface temperatures, these precipitation changes have greatly increased the risk of drought over Africa, southern Europe, East Asia, eastern Australia, Northwest Canada, and southern Brazil. Global land precipitation and continental freshwater discharge show large interannual and inter-decadal variations, with negative anomalies during El Niño and following major volcanic eruptions in 1963, 1982, and 1991; whereas their decadal variations are correlated with the Interdecadal Pacifc Oscillation (IPO) with IPO’s warm phase associated with low land precipitation and continental discharge. The IPO and Atlantic multidecadal variability also dominate multidecadal variations in land aridity, accounting for 90% of the multidecadal variance. CMIP5 multi-model ensemble mean shows decreased precipitation and runoff and increased risk of drought during 1950–2018 over Southwest North America, Central America, northern and central South America (including the Amazon), southern and West Africa, the Mediterranean region, and Southeast Asia; while the northern mid-high latitudes, Southeast South America, and Northwest Australia see increased precipitation and runoff. The consistent spatial patterns between the observed changes and the model-simulated response suggest that many of the observed drying and wetting trends since 1950 may have resulted at least partly from historical external forcing. However, the drying over Southeast Asia and wetting over Northwest Australia are absent in the 21st century projections.more » « less
-
more » « less
Abstract Large uncertainties exist in climate model projections of the Asian summer monsoon (ASM). The El Niño‐Southern Oscillation (ENSO) is an important modulator of the ASM, but the ENSO‐ASM teleconnection is not stationary. Furthermore, teleconnections between ENSO and the East Asian versus South Asian subcomponents of the ASM exhibit distinct characteristics. Therefore, understanding the variability of the ENSO‐ASM teleconnection is critical for anticipating future variations in ASM intensity. To this end, we here use paleoclimate records to extend temporal coverage beyond the instrumental era by millennia. Recently, data assimilation techniques have been applied for the last millennium, which facilitates physically consistent, globally gridded climate reconstructions informed by paleoclimate observations. We use these novel data assimilation products to investigate variations in the ENSO‐ASM relationship over the last 1,000 years. We find that correlations between ENSO and ASM indices are mostly negative in the last millennium, suggesting that strong ASM years are often associated with La Niña events. During periods of weak correlations between ENSO and the East Asian summer monsoon, we observe an El Niño‐like sea surface temperature (SST) pattern in the Pacific. Additionally, SST patterns associated with periods of weak correlations between ENSO and South Asian summer monsoon rainfall are not consistent among data assimilation products. This underscores the importance of developing more precipitation‐sensitive paleoclimate proxies in the Indian subcontinental realm over the last millennium. Our study serves as a baseline for future appraisals of paleoclimate assimilation products and an example of informing our understanding of decadal‐scale ENSO‐ASM teleconnection variability using paleoclimate data sets.
-
more » « less
Abstract Global Climate Models (GCMs) exhibit substantial biases in their simulation of tropical climate. One particularly problematic bias exists in GCMs' simulation of the tropical rainband known as the Intertropical Convergence Zone (ITCZ). Much of the precipitation on Earth falls within the ITCZ, which plays a key role in setting Earth's temperature by affecting global energy transports, and partially dictates dynamics of the largest interannual mode of climate variability: The El Niño‐Southern Oscillation (ENSO). Most GCMs fail to simulate the mean state of the ITCZ correctly, often exhibiting a “double ITCZ bias,” with rainbands both north and south rather than just north of the equator. These tropical mean state biases limit confidence in climate models' simulation of projected future and paleoclimate states, and reduce the utility of these models for understanding present climate dynamics. Adjusting GCM parameterizations of cloud processes and atmospheric convection can reduce tropical biases, as can artificially correcting sea surface temperatures through modifications to air‐sea fluxes (i.e., “flux adjustment”). Here, we argue that a significant portion of these rainfall and circulation biases are rooted in orographic height being biased low due to assumptions made in fitting observed orography onto GCM grids. We demonstrate that making different, and physically defensible, assumptions that raise the orographic height significantly improves model simulation of climatological features such as the ITCZ and North American rainfall as well as the simulation of ENSO. These findings suggest a simple, physically based, and computationally inexpensive method that can improve climate models and projections of future climate.
