Abstract

Abstract. Reconstructions of the late-Holocene climate rely heavily upon proxies that are assumed to be accurately dated by layer counting, such as measurements of tree rings, ice cores, and varved lake sediments. Considerable advances could be achieved if time-uncertain proxies were able to be included within these multiproxy reconstructions, and if time uncertainties were recognized and correctly modeled for proxies commonly treated as free of age model errors. Current approaches for accounting for time uncertainty are generally limited to repeating the reconstruction using each one of an ensemble of age models, thereby inflating the final estimated uncertainty – in effect, each possible age model is given equal weighting. Uncertainties can be reduced by exploiting the inferred space–time covariance structure of the climate to re-weight the possible age models. Here, we demonstrate how Bayesian hierarchical climate reconstruction models can be augmented to account for time-uncertain proxies. Critically, although a priori all age models are given equal probability of being correct, the probabilities associated with the age models are formally updated within the Bayesian framework, thereby reducing uncertainties. Numerical experiments show that updating the age model probabilities decreases uncertainty in the resulting reconstructions, as compared with the current de facto standard of sampling over all age models, provided there is sufficient information from other data sources in the spatial region of the time-uncertain proxy. This approach can readily be generalized to non-layer-counted proxies, such as those derived from marine sediments.

Highlights

  • Large-scale climate reconstructions over the last two millennia rely heavily on climatic proxies that are annually resolved, assumed to be precisely dated by layer counting, and overlap with instrumental climate data – including tree rings, varved sediments, and annually layered ice cores (e.g., NRC, 2006; Jones et al, 2009)

  • The development below is specific for simplicity, the core idea of updating the probabilities associated with time-uncertain proxy data by including an age–depth models (ADMs) level within a Bayesian hierarchical model is general, and discussion provided in Sect. 5 focuses on how the core ideas can be extended, including to radiometrically derived ADMs

  • This results in a Bayesian age model selection (BARCAST+AMS), whereby at each iteration of the Markov chain Monte Carlo (MCMC), an ADM is selected conditional on the data and the current estimate of the climate state

Read more

Summary

Introduction

Large-scale climate reconstructions over the last two millennia rely heavily on climatic proxies that are annually resolved, assumed to be precisely dated by layer counting, and overlap with instrumental climate data – including tree rings, varved sediments, and annually layered ice cores (e.g., NRC, 2006; Jones et al, 2009). The possible age models for a set of proxy records can be constrained by assuming that well-documented global-scale events, such as magnetic reversals, glacial terminations, or tephra layers (Haflidason et al, 2000), are simultaneous within a specified tolerance across the different time-uncertain records Such assumptions reflect the intuition that large-scale climate features are recorded at different locations at “close” to the same time, where the extent to which events are permitted to be asynchronous is given by estimates of dating uncertainties, and reflected in tolerances or penalty functions. 4. the development below is specific for simplicity, the core idea of updating the probabilities associated with time-uncertain proxy data by including an ADM level within a Bayesian hierarchical model is general, and discussion provided in Sect. The development below is specific for simplicity, the core idea of updating the probabilities associated with time-uncertain proxy data by including an ADM level within a Bayesian hierarchical model is general, and discussion provided in Sect. 5 focuses on how the core ideas can be extended, including to radiometrically derived ADMs

Bayesian hierarchical models for climate field reconstructions
Extending BARCAST to include time-uncertain proxies
Addressing miscounted layers
Parallel tempering for ADM selection
Simulation experiments
One time-certain and one time-uncertain proxy
No instrumental data and all proxies time-uncertain
Discussions and extensions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call