skip to main content


Title: Evaluating Proxy Influence in Assimilated Paleoclimate Reconstructions—Testing the Exchangeability of Two Ensembles of Spatial Processes
Climate field reconstructions (CFRs) attempt to estimate spatiotemporal fields of climate variables in the past using climate proxies such as tree rings, ice cores, and corals. Data assimilation (DA) methods are a recent and promising new means of deriving CFRs that optimally fuse climate proxies with climate model output. Despite the growing application of DA-based CFRs, little is understood about how much the assimilated proxies change the statistical properties of the climate model data. To address this question, we propose a robust and computationally efficient method, based on functional data depth, to evaluate differences in the distributions of two spatiotemporal processes. We apply our test to study global and regional proxy influence in DA-based CFRs by comparing the background and analysis states, which are treated as two samples of spatiotemporal fields.We find that the analysis states are significantly altered from the climate-model-based background states due to the assimilation of proxies. Moreover, the difference between the analysis and background states increases with the number of proxies, even in regions far beyond proxy collection sites. Our approach allows us to characterize the added value of proxies, indicating where and when the analysis states are distinct from the background states. Supplementary materials for this article are available online.  more » « less
Award ID(s):
1922758 1830312 1811768
NSF-PAR ID:
10291059
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of the American Statistical Association
Volume:
00
Issue:
0
ISSN:
0162-1459
Page Range / eLocation ID:
1-14
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Climate field reconstructions (CFRs) combine modern observational data with paleoclimatic proxies to estimate climate variables over spatiotemporal grids during time periods when widespread observations of climatic conditions do not exist. The Common Era (CE) has been a period over which many seasonally‐ and annually‐resolved CFRs have been produced on regional to global scales. CFRs over the CE were first produced in the 1970s using dendroclimatic records and linear regression‐based approaches. Since that time, many new CFRs have been produced using a wide range of proxy data sets and reconstruction techniques. We assess the early history of research on CFRs for the CE, which provides context for our review of advances in CFR research over the last two decades. We review efforts to derive gridded hydroclimatic CFRs over continental regions using networks of tree‐ring proxies. We subsequently explore work to produce hemispheric‐ and global‐scale CFRs of surface temperature using multi‐proxy data sets, before specifically reviewing recently‐developed data assimilation techniques and how they have been used to produce simultaneous reconstructions of multiple climatic fields globally. We then review efforts to develop standardized and digitized databases of proxy networks for use in CFR research, before concluding with some thoughts on important next steps for CFR development.

     
    more » « less
  2. 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
  3. The Azores High (AH), a subtropical ridge in the atmosphere over the North Atlantic comprising one node of the North Atlantic Oscillation (NAO) system, has a dominant influence on the weather and climate of the Iberian Peninsula and northwest Africa. The behavior of the entire NAO system over the last millennium has been the subject of much debate in both proxy- and model-based studies. Many studies have focused on the behavior of the entire NAO system, but we focus solely on the behavior of the AH due to its proximity to this region. Other proxies from this region, mainly from Spain and Morocco, have provided details about atmospheric dynamics yet spatiotemporal gaps remain. In this study, we present a continuous, sub-decadally-resolved composite stalagmite carbon isotopic record from three partially overlapping stalagmites from Buraca Gloriosa (BG) cave, western Portugal, situated within the center of the AH, that preserves evidence of regional hydroclimate variability from approximately 800 CE to the present. This composite record, developed from U-Th dating and laminae counting paired with carbon isotopes, primarily reflects effective moisture in western Portugal. Given the close pairing of AH behavior (intensity, size, and location) and moisture transport in this region, the BG composite record allows for a thorough analysis of AH behavior over time. Multidecadal to centennial scale variability in the BG record and state-of-the-art last millennium climate model simulations show considerable coherence with precipitation-sensitive records from Spain and Morocco that, like BG, are strongly influenced by the intensity, size, and location of the AH. Synthesis of model output and proxy data suggests that western Portugal was persistently dry during much of the Medieval Climate Anomaly (MCA; ~850-1250 CE) and Modern era (1850 CE-present) and experienced wetter conditions during Little Ice Age (LIA; ~1400-1850 CE). Even considering age uncertainties from the Iberian Peninsula and northwest Africa proxy records, the apparent timing in the transition from a relatively dry MCA to a wetter LIA is spatially variable across this region, likely due to the non-stationary behavior of the AH system. 
    more » « less
  4. The Azores High (AH), a subtropical ridge in the atmosphere over the North Atlantic comprising one node of the North Atlantic Oscillation (NAO) system, has a dominant influence on the weather and climate of the Iberian Peninsula and northwest Africa. The behavior of the entire NAO system over the last millennium has been the subject of much debate in both proxy- and model-based studies. Many studies have focused on the behavior of the entire NAO system, but we focus solely on the behavior of the AH due to its proximity to this region. Other proxies from this region, mainly from Spain and Morocco, have provided details about atmospheric dynamics yet spatiotemporal gaps remain. In this study, we present a continuous, sub-decadally-resolved composite stalagmite carbon isotopic record from three partially overlapping stalagmites from Buraca Gloriosa (BG) cave, western Portugal, situated within the center of the AH, that preserves evidence of regional hydroclimate variability from approximately 800 CE to the present. This composite record, developed from U-Th dating and laminae counting paired with carbon isotopes, primarily reflects effective moisture in western Portugal. Given the close pairing of AH behavior (intensity, size, and location) and moisture transport in this region, the BG composite record allows for a thorough analysis of AH behavior over time. Multidecadal to centennial scale variability in the BG record and state-of-the-art last millennium climate model simulations show considerable coherence with precipitation-sensitive records from Spain and Morocco that, like BG, are strongly influenced by the intensity, size, and location of the AH. Synthesis of model output and proxy data suggests that western Portugal was persistently dry during much of the Medieval Climate Anomaly (MCA; ~850-1250 CE) and Modern era (1850 CE-present) and experienced wetter conditions during Little Ice Age (LIA; ~1400-1850 CE). Even considering age uncertainties from the Iberian Peninsula and northwest Africa proxy records, the apparent timing in the transition from a relatively dry MCA to a wetter LIA is spatially variable across this region, likely due to the non-stationary behavior of the AH system. 
    more » « less
  5. Abstract. Paleoclimatic records provide valuable information about Holocene climate, revealing aspects of climate variability for a multitude of sites around the world. However, such data also possess limitations. Proxy networks are spatially uneven, seasonally biased, uncertain in time, and present a variety of challenges when used in concert to illustrate the complex variations of past climate. Paleoclimatic data assimilation provides one approach to reconstructing past climate that can account for the diversenature of proxy records while maintaining the physics-based covariancestructures simulated by climate models. Here, we use paleoclimate dataassimilation to create a spatially complete reconstruction of temperatureover the past 12 000 years using proxy data from the Temperature 12k database and output from transient climate model simulations. Following the last glacial period, the reconstruction shows Holocene temperatures warming to a peak near 6400 years ago followed by a slow cooling toward the present day, supporting a mid-Holocene which is at least as warm as the preindustrial. Sensitivity tests show that if proxies have an overlooked summer bias, some apparent mid-Holocene warmth could actually represent summer trends rather than annual mean trends. Regardless, the potential effects of proxy seasonal biases are insufficient to align the reconstructed global mean temperature with the warming trends seen in transient model simulations. 
    more » « less