A Bayesian algorithm for reconstructing climate anomalies in space and time. Part I

Development and applications to paleoclimate reconstruction problems

Martin Patrick Tingley, Peter Huybers

Research output: Contribution to journalArticle

98 Citations (Scopus)

Abstract

Reconstructing the spatial pattern of a climate field through time from a dataset of overlapping instrumental and climate proxy time series is a nontrivial statistical problem. The need to transform the proxy observations into estimates of the climate field, and the fact that the observed time series are not uniformly distributed in space, further complicate the analysis. Current leading approaches to this problem are based on estimating the full covariance matrix between the proxy time series and instrumental time series over a "calibration" interval and then using this covariance matrix in the context of a linear regression to predict the missing instrumental values from the proxy observations for years prior to instrumental coverage. A fundamentally different approach to this problem is formulated by specifying parametric forms for the spatial covariance and temporal evolution of the climate field, as well as "observation equations" describing the relationship between the data types and the corresponding true values of the climate field. A hierarchical Bayesian model is used to assimilate both proxy and instrumental datasets and to estimate the probability distribution of all model parameters and the climate field through time on a regular spatial grid. The output from this approach includes an estimate of the full covariance structure of the climate field and model parameters as well as diagnostics that estimate the utility of the different proxy time series. This methodology is demonstrated using an instrumental surface temperature dataset after corrupting a number of the time series to mimic proxy observations. The results are compared to those achieved using the regularized expectation-maximization algorithm, and in these experiments the Bayesian algorithm produces reconstructions with greater skill. The assumptions underlying these two methodologies and the results of applying each to simple surrogate datasets are explored in greater detail in Part II.

Original languageEnglish (US)
Pages (from-to)2759-2781
Number of pages23
JournalJournal of Climate
Volume23
Issue number10
DOIs
StatePublished - May 1 2010

Fingerprint

paleoclimate
anomaly
time series
climate
matrix
methodology
temporal evolution
surface temperature
transform
calibration
experiment

All Science Journal Classification (ASJC) codes

  • Atmospheric Science

Cite this

@article{aeec3857e0a048ee886b0459b4b11486,
title = "A Bayesian algorithm for reconstructing climate anomalies in space and time. Part I: Development and applications to paleoclimate reconstruction problems",
abstract = "Reconstructing the spatial pattern of a climate field through time from a dataset of overlapping instrumental and climate proxy time series is a nontrivial statistical problem. The need to transform the proxy observations into estimates of the climate field, and the fact that the observed time series are not uniformly distributed in space, further complicate the analysis. Current leading approaches to this problem are based on estimating the full covariance matrix between the proxy time series and instrumental time series over a {"}calibration{"} interval and then using this covariance matrix in the context of a linear regression to predict the missing instrumental values from the proxy observations for years prior to instrumental coverage. A fundamentally different approach to this problem is formulated by specifying parametric forms for the spatial covariance and temporal evolution of the climate field, as well as {"}observation equations{"} describing the relationship between the data types and the corresponding true values of the climate field. A hierarchical Bayesian model is used to assimilate both proxy and instrumental datasets and to estimate the probability distribution of all model parameters and the climate field through time on a regular spatial grid. The output from this approach includes an estimate of the full covariance structure of the climate field and model parameters as well as diagnostics that estimate the utility of the different proxy time series. This methodology is demonstrated using an instrumental surface temperature dataset after corrupting a number of the time series to mimic proxy observations. The results are compared to those achieved using the regularized expectation-maximization algorithm, and in these experiments the Bayesian algorithm produces reconstructions with greater skill. The assumptions underlying these two methodologies and the results of applying each to simple surrogate datasets are explored in greater detail in Part II.",
author = "Tingley, {Martin Patrick} and Peter Huybers",
year = "2010",
month = "5",
day = "1",
doi = "10.1175/2009JCLI3015.1",
language = "English (US)",
volume = "23",
pages = "2759--2781",
journal = "Journal of Climate",
issn = "0894-8755",
publisher = "American Meteorological Society",
number = "10",

}

A Bayesian algorithm for reconstructing climate anomalies in space and time. Part I : Development and applications to paleoclimate reconstruction problems. / Tingley, Martin Patrick; Huybers, Peter.

In: Journal of Climate, Vol. 23, No. 10, 01.05.2010, p. 2759-2781.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A Bayesian algorithm for reconstructing climate anomalies in space and time. Part I

T2 - Development and applications to paleoclimate reconstruction problems

AU - Tingley, Martin Patrick

AU - Huybers, Peter

PY - 2010/5/1

Y1 - 2010/5/1

N2 - Reconstructing the spatial pattern of a climate field through time from a dataset of overlapping instrumental and climate proxy time series is a nontrivial statistical problem. The need to transform the proxy observations into estimates of the climate field, and the fact that the observed time series are not uniformly distributed in space, further complicate the analysis. Current leading approaches to this problem are based on estimating the full covariance matrix between the proxy time series and instrumental time series over a "calibration" interval and then using this covariance matrix in the context of a linear regression to predict the missing instrumental values from the proxy observations for years prior to instrumental coverage. A fundamentally different approach to this problem is formulated by specifying parametric forms for the spatial covariance and temporal evolution of the climate field, as well as "observation equations" describing the relationship between the data types and the corresponding true values of the climate field. A hierarchical Bayesian model is used to assimilate both proxy and instrumental datasets and to estimate the probability distribution of all model parameters and the climate field through time on a regular spatial grid. The output from this approach includes an estimate of the full covariance structure of the climate field and model parameters as well as diagnostics that estimate the utility of the different proxy time series. This methodology is demonstrated using an instrumental surface temperature dataset after corrupting a number of the time series to mimic proxy observations. The results are compared to those achieved using the regularized expectation-maximization algorithm, and in these experiments the Bayesian algorithm produces reconstructions with greater skill. The assumptions underlying these two methodologies and the results of applying each to simple surrogate datasets are explored in greater detail in Part II.

AB - Reconstructing the spatial pattern of a climate field through time from a dataset of overlapping instrumental and climate proxy time series is a nontrivial statistical problem. The need to transform the proxy observations into estimates of the climate field, and the fact that the observed time series are not uniformly distributed in space, further complicate the analysis. Current leading approaches to this problem are based on estimating the full covariance matrix between the proxy time series and instrumental time series over a "calibration" interval and then using this covariance matrix in the context of a linear regression to predict the missing instrumental values from the proxy observations for years prior to instrumental coverage. A fundamentally different approach to this problem is formulated by specifying parametric forms for the spatial covariance and temporal evolution of the climate field, as well as "observation equations" describing the relationship between the data types and the corresponding true values of the climate field. A hierarchical Bayesian model is used to assimilate both proxy and instrumental datasets and to estimate the probability distribution of all model parameters and the climate field through time on a regular spatial grid. The output from this approach includes an estimate of the full covariance structure of the climate field and model parameters as well as diagnostics that estimate the utility of the different proxy time series. This methodology is demonstrated using an instrumental surface temperature dataset after corrupting a number of the time series to mimic proxy observations. The results are compared to those achieved using the regularized expectation-maximization algorithm, and in these experiments the Bayesian algorithm produces reconstructions with greater skill. The assumptions underlying these two methodologies and the results of applying each to simple surrogate datasets are explored in greater detail in Part II.

UR - http://www.scopus.com/inward/record.url?scp=77954378941&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77954378941&partnerID=8YFLogxK

U2 - 10.1175/2009JCLI3015.1

DO - 10.1175/2009JCLI3015.1

M3 - Article

VL - 23

SP - 2759

EP - 2781

JO - Journal of Climate

JF - Journal of Climate

SN - 0894-8755

IS - 10

ER -