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: 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
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 “Online” data assimilation (DA) is used to generate a seasonal-resolution reanalysis dataset over the last millennium by combining forecasts from an ocean–atmosphere–sea-ice coupled linear inverse model with climate proxy records. Instrumental verification reveals that this reconstruction achieves the highest correlation skill, while using fewer proxies, in surface temperature reconstructions compared to other paleo-DA products, particularly during boreal winter when proxy data are scarce. Reconstructed ocean and sea-ice variables also have high correlation with instrumental and satellite datasets. Verification against independent proxy records shows that reconstruction skill is robust throughout the last millennium. Analysis of the results reveals that the method effectively captures the seasonal evolution and amplitude of El Niño events, seasonal temperature trends that are consistent with orbital forcing over the last millennium, and polar-amplified cooling in the transition from the Medieval Climate Anomaly to the Little Ice Age. 
    more » « less
  3. Abstract. Paleoclimate data assimilation (DA) is a tool for reconstructing past climates that directly integrates proxy records with climate model output. Despite the potential for DA to expand the scope of quantitative paleoclimatology, these methods remain difficult to implement in practice due to the multi-faceted requirements and data handling necessary for DA reconstructions, the diversity of DA methods, and the need for computationally efficient algorithms. Here, we present DASH, a MATLAB toolbox designed to facilitate paleoclimate DA analyses. DASH provides command line and scripting tools that implement common tasks in DA workflows. The toolbox is highly modular and is not built around any specific analysis, and thus DASH supports paleoclimate DA for a wide variety of time periods, spatial regions, proxy networks, and algorithms. DASH includes tools for integrating and cataloguing data stored in disparate formats, building state vector ensembles, and running proxy (system) forward models. The toolbox also provides optimized algorithms for implementing ensemble Kalman filters, particle filters, and optimal sensor analyses with variable and modular parameters. This paper reviews the key components of the DASH toolbox and presents examples illustrating DASH's use for paleoclimate DA applications. 
    more » « less
  4. Global hydroclimate significantly differed from modern climate during the mid-Holocene (6 ka) and Last Glacial Maximum (21 ka). Consequently, both periods have been described as either a partial or reverse analogue for current climate change. To reconstruct past hydroclimate, an offline paleoclimate data assimilation methodology is applied to a dataset of 216 lake status records which provide relative estimates of water level change. The proxy observations are integrated with the climate dynamics of two transient simulations (TraCE-21ka and HadCM3) using a multivariate proxy system model (PSM) which estimates relative lake status from available climate simulation variables. The resulting DAMP-21ka (Data Assimilation of Moisture Patterns 21 000–0 BP) reanalysis reconstructs annual lake status and precipitation values at 500-year resolution and represents the first application of the methodology to global hydroclimate on timescales spanning the Holocene and longer. Validation using Pearson's correlation coefficients indicates that the reconstruction (0.24) is more skillful, on average, than model simulations (0.09), particularly in portions of North America and east Africa, where data density is high and proxy–model disagreement is prominent during the Holocene. Results of the PSM and assimilation are used to evaluate climatic controls on lake status, spatiotemporal patterns of moisture variability, and proxy–model disagreement. During the mid-Holocene, wetter conditions are reconstructed for northern and eastern Africa, Asia, and southern Australia, but in contrast to the model prior, negative anomalies are observed in North America, resulting in drier-than-modern conditions throughout the Northern Hemisphere midlatitudes. Proxy–model disagreement in western North America may reflect a bias in model simulations to stronger sea level pressure gradients in the North Pacific during the mid-Holocene. The data assimilation framework is able to reconcile these differences by integrating the constraints of proxy observations with the dynamics of the model prior to produce a more robust estimation of hydroclimate variability during the past 21 000 years. 
    more » « less
  5. Abstract. Climate field reconstruction (CFR) refers to the estimation of spatiotemporal climate fields (such as surface temperature) from a collection of pointwise paleoclimate proxy datasets. Such reconstructions can provide rich information on climate dynamics and provide an out-of-sample validation of climate models. However, most CFR workflows are complex and time-consuming, as they involve (i) preprocessing of the proxy records, climate model simulations, and instrumental observations; (ii) application of one or more statistical methods; and (iii) analysis and visualization of the reconstruction results. Historically, this process has lacked transparency and accessibility, limiting reproducibility and experimentation by non-specialists. This article presents an open-source and object-oriented Python package called cfr that aims to make CFR workflows easy to understand and conduct, saving climatologists from technical details and facilitating efficient and reproducible research. cfr provides user-friendly utilities for common CFR tasks such as proxy and climate data analysis and visualization, proxy system modeling, and modularized workflows for multiple reconstruction methods, enabling methodological intercomparisons within the same framework. The package is supported with extensive documentation of the application programming interface (API) and a growing number of tutorial notebooks illustrating its usage. As an example, we present two cfr-driven reconstruction experiments using the PAGES 2k temperature database applying the last millennium reanalysis (LMR) paleoclimate data assimilation (PDA) framework and the graphical expectation–maximization (GraphEM) algorithm, respectively. 
    more » « less