### 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 language | English (US) |
---|---|

Pages (from-to) | 2759-2781 |

Number of pages | 23 |

Journal | Journal of Climate |

Volume | 23 |

Issue number | 10 |

DOIs | |

State | Published - May 1 2010 |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Atmospheric Science

### Cite this

*Journal of Climate*,

*23*(10), 2759-2781. https://doi.org/10.1175/2009JCLI3015.1

}

*Journal of Climate*, vol. 23, no. 10, pp. 2759-2781. https://doi.org/10.1175/2009JCLI3015.1

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

Research output: Contribution to journal › Article

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 P.

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

AN - SCOPUS:77954378941

VL - 23

SP - 2759

EP - 2781

JO - Journal of Climate

JF - Journal of Climate

SN - 0894-8755

IS - 10

ER -