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.

Peatlands are a major terrestrial carbon store and a persistent natural carbon sink during the Holocene, but there is considerable uncertainty over the fate of peatland carbon in a changing climate. It is generally assumed that higher temperatures will increase peat decay, causing a positive feedback to climate warming and contributing to the global positive carbon cycle feedback. Here we use a new extensive database of peat profiles across northern high latitudes to examine spatial and temporal patterns of carbon accumulation over the past millennium. Opposite to expectations, our results indicate a small negative carbon cycle feedback from past changes in the long-term accumulation rates of northern peatlands. Total carbon accumulated over the last 1000 yr is linearly related to contemporary growing season length and photosynthetically active radiation, suggesting that variability in net primary productivity is more important than decomposition in determining long-term carbon accumulation. Furthermore, northern peatland carbon sequestration rate declined over the climate transition from the Medieval Climate Anomaly (MCA) to the Little Ice Age (LIA), probably because of lower LIA temperatures combined with increased cloudiness suppressing net primary productivity. Other factors including changing moisture status, peatland distribution, fire, nitrogen deposition, permafrost thaw and methane emissions will also influence future peatland carbon cycle feedbacks, but our data suggest that the carbon sequestration rate could increase over many areas of northern peatlands in a warmer future.