Bayesian Analysis

Flexible paleoclimate age-depth models using an autoregressive gamma process

Maarten Blaauw and J. Andrés Christen

Full-text: Open access


Radiocarbon dating is routinely used in paleoecology to build chronologies of lake and peat sediments, aiming at inferring a model that would relate the sediment depth with its age. We present a new approach for chronology building (called "Bacon") that has received enthusiastic attention by paleoecologists. Our methodology is based on controlling core accumulation rates using a gamma autoregressive semiparametric model with an arbitrary number of subdivisions along the sediment. Using prior knowledge about accumulation rates is crucial and informative priors are routinely used. Since many sediment cores are currently analyzed, using different data sets and prior distributions, a robust (adaptive) MCMC is very useful. We use the t-walk (Christen and Fox, 2010), a self adjusting, robust MCMC sampling algorithm, that works acceptably well in many situations. Outliers are also addressed using a recent approach that considers a Student-t model for radiocarbon data. Two examples are presented here, that of a peat core and a core from a lake, and our results are compared with other approaches.

Article information

Bayesian Anal., Volume 6, Number 3 (2011), 457-474.

First available in Project Euclid: 13 June 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Primary: 60K40: Other physical applications of random processes
Secondary: 62F15: Bayesian inference 62M10: Time series, auto-correlation, regression, etc. [See also 91B84] 62P12: Applications to environmental and related topics 86A60: Geological problems

AR Gamma process Radiocarbon Paleoecology Age-Depth models


Blaauw, Maarten; Christen, J. Andrés. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6 (2011), no. 3, 457--474. doi:10.1214/11-BA618.

Export citation


  • Andrieu, C. and Thoms, J. (2008). "A tutorial on adaptive MCMC." Statistics and Computing, 18(4): 343–373.
  • Barndorff-Nielsen, O. E. and Shephard, N. (2001). "Non-Gaussian Ornstein-Uhlenbeck-based models and some of their uses in financial economics (with discussion)." Journal of the Royal Statistical Society, Series B, 63: 167–241.
  • –- (2003). "Integrated OU processes and non-Gaussian OU-based stochastic volatility models." Scandinavian Journal of Statistics, 30: 277–295.
  • Belyea, L. W. and Clymo, R. S. (2001). "Feedback control of the rate of peat formation." Proceedings of the Royal Society of London. Series B, Biological sciences, 268: 1315–1321.
  • Blaauw, M. (2010). "Methods and code for 'classical' age-modelling of radiocarbon sequences." Quaternary Geochronology, 5: 512–518.
  • Blaauw, M., Bakker, R., Christen, J. A., Hall, V. A., and van der Plicht, J. (2007). "A Bayesian framework for age modeling of radiocarbon-dated peat deposits: case studies from the Netherlands." Radiocarbon, 49(2): 357–367.
  • Blaauw, M. and Christen, J. A. (2005). "Radiocarbon peat chronologies and environmental change." Applied Statistics, 54(5): 805–816.
  • Blaauw, M., Christen, J. A., Ampel, L., Veres, D., Hughen, K., and Preusser, F. (2010). "Were last glacial climate events simultaneous between Greenland and France? A quantitative comparison using non-tuned chronologies." Journal of Quaternary Science, 25: 387–394.
  • Blaauw, M., Christen, J. A., Mauquoy, D., van der Plicht, J., and Bennett, K. D. (2007). "Testing the timing of radiocarbon-dated events between proxy archives." The Holocene, 17: 283–288.
  • Blaauw, M., Heuvelink, G. B. M., Mauquoy, D., van der Plicht, J., and van Geel, B. (2003). "A numerical approach to 14C wiggle-match dating of organic deposits: best fits and confidence intervals." Quaternary Science Reviews, 22(14): 1485–1500.
  • Blockley, S. P. E., Blaauw, M., Bronk Ramsey, C., and van der Plicht, J. (2007). "Assessing uncertainties in age modelling sedimentary records in the Lateglacial and early Holocene." Quaternary Science Reviews, 26: 1915–1926.
  • Bronk Ramsey, C. (2008). "Deposition models for chronological records." Quaternary Science Reviews, 27(1-2): 42–60.
  • –- (2009). "Bayesian analysis of radiocarbon dates." Radiocarbon, 51(1): 337–360.
  • –- (2009). "Dealing with outliers and offsets in radiocarbon dating." Radiocarbon, 51(3): 1023–1045.
  • Cai, B. and Dunsun, D. B. (2007). "Bayesian Multivariate Isotonic Regression Splines: Applications to Carcinogenicity Studies." Journal of the American Statistical Association, 102(480): 1158–1171.
  • Chambers, F. M., Mauquoy, D., Brain, S. A., Blaauw, M., and Daniell, J. R. G. (2007). "Globally synchronous climate change 2800 years ago: proxy data from peat in South America." Earth and Planetary Science Letters, 253: 439–444.
  • Charman, D., Barber, K., Blaauw, M., Langdon, P., Mauquoy, D., Daley, T., Hughes, P., and Karofeld, E. (2009). "Climate drivers for peatland palaeoclimate records." Quaternary Science Reviews, 28: 1811–1819.
  • Charman, D., Beilman, D. W., Blaauw, M., Booth, B., Brewer, S., Chambers, F. M., Christen, J. A., Gallego-Sala, A., Harrison, S. P., Hughes, P. D. M., Jackson, S. T., Korhola, A., Mauquoy, D., Mitchell, F. J. G., Prentice, I. C., van der Linden, M., de Vleeschouwer, F., Yu, Z. C., Bauer, I. E., Garneau, M., Karofeld, E., Le Roux, G., Moschen, R., Nichols, J. E., Phadtare, N. R., Rausch, N., Swindles, G. T., Väliranta, M., van Bellen, S., and van Geel, B. (2011). "Carbon–cycle implications of climate–driven changes in peat accumulation during the last millennium." Nature, (Revision).
  • Christen, J., Clymo, D., and Litton, C. (1995). "A Bayesian approach to the use of C14 dates in the estimation of the age of peat." Radiocarbon, 37: 431–442.
  • Christen, J. A. (1994). "Summarising a set of radiocarbon determinations: a robust approach." Applied Statistics, 43(3): 489–503.
  • Christen, J. A. and Fox, C. (2010). "A general purpose sampling algorithm for continuous distributions (the t-walk)." Bayesian Analisis, 4(2): 263–282.
  • Christen, J. A. and Pérez, S. (2009). "A new robust statistical model for radiocarbon data." Radiocarbon, 51(3): 1047–1059.
  • de Vleeschouwer, F., Pazdur, A., Luthers, C., Streel, M., Mauquoy, D., Wastiaux, C., Le Roux, G., Moschen, R., Blaauw, M., Pawlyta, J., Sikorski, J., and Piotrowska, N. (2010). "A millennial record of environmental change recorded in peat deposits from the Misten Bog (east Belgium)." Quaternary International, (Submitted).
  • FiriGroup (2003). "The Third International Radiocarbon Intercomparison (Tiri) and the Fourth International Radiocarbon (Fir) – 1999-2002 – Results, analysis and conclusions." Radiocarbon, 45(2): 135–328.
  • Gamerman, D. and Lopes, H. (2006). Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference. London: Chapman and Hall.
  • Geyer, C. J. (1992). "Practical Markov Chain Monte Carlo." Statistical Science, 43: 179–189.
  • Harkness, D., Miller, B. F., and Tipping, R. M. (1997). "NERC radiocarbon measurements 1977-1988." Quaternary Science Reviews, 16: 925–927.
  • Haslett, J. and Parnell, A. (2008). "A simple monotone process with application to radiocarbon-dated depth chronologies." Applied Statistics, 57(4): 399–418.
  • Holmes, C. C. and Heard, N. A. (2003). "Generalized monotonic regression using random change points." Statistics in Medicine, 22: 623–838.
  • Jones, V. J., Stevenson, A. C., and Battarbee, R. W. (1989). "Acidification of lakes in Galloway, southwest Scotland: a diatom and pollen study of the post-glacial history of the Round Loch of Glenhead." Journal of Ecology, 77: 1–23.
  • Lavine, M. and Mokus, A. (1995). "A nonparametric Bayes method for isotonic regression." Journal of Statistical Planning and Inference, 46: 235–248.
  • Liu, J. S. (2001). Monte Carlo Strategies in Scientific Computing. New York: Springer.
  • Lowe, J. J. and Walker, M. J. C. (1997). Reconstructing Quaternary Environments (2nd edition). Prentice Hall.
  • Mena, R. H. and Walker, S. G. (2004). "A density function connected with a non–negative self-decomposable random variable." Journal of Statistical Computation and Simulation, 74(10): 765–775.
  • –- (2005). "Stationary autoregressive models via a Bayesian nonparametric approach." Journal of Time Series Analysis, 26: 789–805.
  • Parnell, A., Haslett, J., Allen, J., Buck, C., and Huntley, B. (2008). "A flexible approach to assessing synchroneity of past events using Bayesian reconstructions of sedimentation history." Quaternary Science Reviews, 27(19-20): 1872–1885.
  • Plunkett, G., Whitehouse, N. J., Charman, D. J., Blaauw, M., Kelly, E., and Mulhall, I. (2009). "A multi-proxy palaeoenvironmental investigation of the findspot of an Iron Age bog body from Oldcroghan, Co. Offaly, Ireland." Journal of Archaeological Science, 36: 265–277.
  • Reimer, P. J., Baillie, M. G. L., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Bronk Ramsey, C., Buck, C. E., Burr, G. S., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Hajdas, I., Heaton, T. J., Hogg, A. G., Hughen, K. A., Kaiser, K. F., Kromer, B., McCormac, F. G., Manning, S. W., Reimer, R. W., Richards, D. A., Southon, J. R., Talamo, S., Turney, C. S. M., van der Plicht, J., and Weyhenmeyer, C. E. (2009). "IntCal09 and Marine09 radiocarbon age calibration curves, 0–50,000 years cal BP." Radiocarbon, 51: 1111–1150.
  • Sillasoo, Ü., Mauquoy, D., Blundell, A., Charman, D., Blaauw, M., Daniell, J. G. R., Toms, P., Newberry, J., Chambers, F. M., and Karofeld, E. (2007). "Peat multi-proxy from Männikjärve bog as indicators of Late Holocene climate changes in Estonia." Boreas, 36: 20–37.
  • Stevenson, A. C., Jones, V. J., and Battarbee, R. W. (1990). "The cause of peat erosion: a palaeolimnological approach." New Phytologist, 114: 727–735.
  • Wohlfarth, B., Blaauw, M., Davies, S. M., Andersson, M., Wastegård, S., Hormes, A., and Possnert, G. (2006). "Constraining the age of Lateglacial and early Holocene pollen zones and tephra horizons in southern Sweden with Bayesian probability methods." Journal of Quaternary Science, 21: 321–334.
  • Wohlfarth, B., Veres, D., Ampel, L., Lacourse, T., Blaauw, M., Preusser, F., Andrieu-Ponel, V., Kéravis, D., Lallier-Vergès, E., Björck, S., Davies, S. M., de Beaulieu, J. L., Risberg, J., Hormes, A., Kasper, H. U., Possnert, G., Reille, M., Thouveny, N., and Zander, A. (2008). "Rapid ecosystem response to abrupt climate changes during the last glacial period in Western Europe, 40–16 kyr BP." Geology, 36: 407–410.
  • Yeloff, D., Bennett, K. D., Blaauw, M., Mauquoy, D., Sillasoo, Ü., van der Plicht, J., and van Geel, B. (2006). "High precision C-14 dating of Holocene peat deposits: a comparison of Bayesian calibration and wiggle-matching approaches." Quaternary Geochronology, 1: 222–235.