A probabilistic model of chronological errors in layer-counted climate proxies: Applications to annually banded coral archives

60Citations
Citations of this article
49Readers
Mendeley users who have this article in their library.

Abstract

The ability to precisely date climate proxies is central to the reconstruction of past climate variations. To a degree, all climate proxies are affected by age uncertainties, which are seldom quantified. This article proposes a probabilistic age model for proxies based on layer-counted chronologies, and explores its use for annually banded coral archives. The model considers both missing and doubly counted growth increments (represented as independent processes), accommodates various assumptions about error rates, and allows one to quantify the impact of chronological uncertainties on different diagnostics of variability. In the case of a single coral record, we find that time uncertainties primarily affect high-frequency signals but also significantly bias the estimate of decadal signals. We further explore tuning to an independent, tree-ring-based chronology as a way to identify an optimal age model. A synthetic pseudocoral network is used as testing ground to quantify uncertainties in the estimation of spatiotemporal patterns of variability. Even for small error rates, the amplitude of multidecadal variability is systematically overestimated at the expense of interannual variability (El Niño-Southern Oscillation, or ENSO, in this case), artificially flattening its spectrum at periods longer than 10 years. An optimization approach to correct chronological errors in coherent multivariate records is presented and validated in idealized cases, though it is found difficult to apply in practice due to the large number of solutions. We close with a discussion of possible extensions of this model and connections to existing strategies for modeling age uncertainties. ©Author(s) 2014. CC Attribution 3.0 License.

Figures

  • Fig. 1. Coral core X-radiographs: (a) top of a core from Mange Reef, Tanzania (Cole, unpublished). Each layer shows the changing density of the coral growth over a year. Observation xi could be the δ18O measured from the ith layer (most likely corresponding to 1i = 1 year). The collection date t1 and the oldest one tn are also shown to illustrate our notation; (b) top of a core from Onotoa, Republic of Kiribati (Thompson, unpublished) showing weak seasonality and annual density banding; (c) Fanning Island fossil coral V33, aged ∼ 6500 years old (Cobb et al., 2013), with a varying growth axis, which typically requires multiple sampling transects.
  • Fig. 2. Illustration of the phase distortion of coherent signals due to miscounted annual bands with θ1 = θ2 = 0.05 and a collection date t1 = 100. (a) Five hypothetical and identical harmonic signals shifted vertically for clarity; (b) the same signals after introducing age perturbations and making the (erroneous) assumption that 1 layer amounts to 1 year (note the loss of coherency reflected in the ensemble mean in gray, which is performed as part of composites in most climate reconstructions); (c) the same age-perturbed signals, after correcting for miscounted bands and interpolating missing data with SSAM.
  • Fig. 3. Ensemble realization (top) and power spectra (bottom) of age-perturbed coral data for (a, b) Havannah Island (Isdale et al., 1998) and (c, d) Palmyra Island δ18O record from 1636 to 1703 (Cobb et al., 2003), where the top layer is U/Th dated. The original data and corresponding spectral estimates are shown in red, and the 95 % confidence intervals (CIs) and the 25–75 % interquartile range (IQR) of the perturbed data ensemble and corresponding spectra are shown by the shaded areas. Power spectra are represented in variance-preserving coordinates (that is, the area under the curve is an exact measure of the variance contained in each frequency band).
  • Fig. 4. Palmyra relict coral δ18O records (top): published data SB3 in red, SB8 in orange and SB13 in brown (Cobb et al., 2003); the shaded area is the 95 % confidence interval of the ensemble of realizations with final time shifted tf ± 10 years and θ1 = θ2 = 0.05, and the blue chronology corresponds to the realizations that maximize the correlation with the first principal component of the NADA record (bottom) over each continuous section. For SB3, there were 2 added years and 2 removed years, and tf was shifted by −2 years. For SB8, there were 2 added years and 2 removed years, and tf was shifted by 5 years. For SB13, there were 1 added year and 2 removed years, and tf was shifted by −6 years. The correlation coefficients with NADA PC1 are also shown, wherein ρ is computed using the aligned Palmyra chronologies, and ρ0 using the original chronologies.
  • Table 1. List of coral locations used for the EOF analysis of Fig. 5. The records were obtained from the database of Emile-Geay and Eshleman (2013) after selecting data sets that record δ18O and cover the period 1860–1980.
  • Fig. 5. Spatiotemporal uncertainty quantification on a pseudocoral network. (a) EOF loadings (circles) corresponding to the ENSO mode of an ensemble of age-perturbed pseudocoral records with miscounting rate θ = 0.05. EOF loadings for error-free data are shown in light colors circled in white, while the median and 95 % quantile are shown by dark disks and black-circled disks, respectively. Contours depict the SST field associated with the mode’s principal component (PC) (b), whose power spectrum is shown in (c). Results for the time-uncertain ensemble are shown in blue: median (solid line), 95 % confidence interval (light-filled area) and interquartile range [25–75 %] (dark-filled area). Results for the original (error-free) data set are depicted by solid red lines. Dashed red lines denote χ2 error estimates for the MTM spectrum of the error-free data set.
  • Fig. 6. Spectral density of the ENSO mode PC for θ = 0.001. Compare to Fig. 5, bottom right.
  • Table 2. Complexity as a function of θ and the matrix dimension np. Infinity symbols denote numbers that exceed machine capacity.

References Powered by Scopus

A Pliocene-Pleistocene stack of 57 globally distributed benthic δ <sup>18</sup>O records

6746Citations
N/AReaders
Get full text

Error Detecting and Error Correcting Codes

3995Citations
N/AReaders
Get full text

Spectrum Estimation and Harmonic Analysis

3476Citations
N/AReaders
Get full text

Cited by Powered by Scopus

Holocene global mean surface temperature, a multi-method reconstruction approach

230Citations
N/AReaders
Get full text

Tropical sea surface temperatures for the past four centuries reconstructed from coral archives

210Citations
N/AReaders
Get full text

Towards a rigorous understanding of societal responses to climate change

195Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Comboul, M., Emile-Geay, J., Evans, M. N., Mirnateghi, N., Cobb, K. M., & Thompson, D. M. (2014). A probabilistic model of chronological errors in layer-counted climate proxies: Applications to annually banded coral archives. Climate of the Past, 10(2), 825–841. https://doi.org/10.5194/cp-10-825-2014

Readers over time

‘13‘14‘15‘16‘17‘18‘19‘20‘21‘22‘23‘240481216

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 21

57%

Researcher 14

38%

Professor / Associate Prof. 2

5%

Readers' Discipline

Tooltip

Earth and Planetary Sciences 24

67%

Environmental Science 7

19%

Agricultural and Biological Sciences 3

8%

Physics and Astronomy 2

6%

Save time finding and organizing research with Mendeley

Sign up for free
0