Non-parametric calibration of multiple related radiocarbon determinations and their calendar age summarisation

Due to fluctuations in past radiocarbon ($^{14}$C) levels, calibration is required to convert $^{14}$C determinations $X_i$ into calendar ages $\theta_i$. In many studies, we wish to calibrate a set of related samples taken from the same site or context, which have calendar ages drawn from the same shared, but unknown, density $f(\theta)$. Calibration of $X_1, \ldots, X_n$ can be improved significantly by incorporating the knowledge that the samples are related. Furthermore, summary estimates of the underlying shared $f(\theta)$ can provide valuable information on changes in population size/activity over time. Most current approaches require a parametric specification for $f(\theta)$ which is often not appropriate. We develop a rigorous non-parametric Bayesian approach using a Dirichlet process mixture model, with slice sampling to address the multimodality typical within $^{14}$C calibration. Our approach simultaneously calibrates the set of $^{14}$C determinations and provides a predictive estimate for the underlying calendar age of a future sample. We show, in a simulation study, the improvement in calendar age estimation when jointly calibrating related samples using our approach, compared with calibration of each $^{14}$C determination independently. We also illustrate the use of the predictive calendar age estimate to provide insight on activity levels over time using three real-life case studies.Comment: Submitted to Journal Royal Statistical Society Series C (in revision). Version 2 is the revised version in response to referee reports. Code is available at https://github.com/TJHeaton/NonparametricCalibratio

Updated Cariaco Basin C-14 calibration dataset from 0-60 cal kyr BP

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Hughen, K. A., & Heaton, T. J. Updated Cariaco Basin C-14 calibration dataset from 0-60 cal kyr BP. Radiocarbon, 62(4), (2020): 1001-1043, doi:10.1017/RDC.2020.53.We present new updates to the calendar and radiocarbon (14C) chronologies for the Cariaco Basin, Venezuela. Calendar ages were generated by tuning abrupt climate shifts in Cariaco Basin sediments to those in speleothems from Hulu Cave. After the original Cariaco-Hulu calendar age model was published, Hulu Cave δ18O records have been augmented with increased temporal resolution and a greater number of U/Th dates. These updated Hulu Cave records provide increased accuracy as well as precision in the final Cariaco calendar age model. The depth scale for the Ocean Drilling Program Site 1002D sediment core, the primary source of samples for 14C dating, has been corrected to account for missing sediment from a core break, eliminating age-depth anomalies that afflicted the earlier calendar age models. Individual 14C dates for the Cariaco Basin remain unchanged from previous papers, although detailed comparisons of the Cariaco calibration dataset to those from Hulu Cave and Lake Suigetsu suggest that the Cariaco marine reservoir age may have shifted systematically during the past. We describe these recent changes to the Cariaco datasets and provide the data in a comprehensive format that will facilitate use by the community.K.A. Hughen was supported by funds from U.S. NSF grant #OCE-1657191, and by the Investment in Science Fund at WHOI. T.J. Heaton is supported by a Leverhulme Trust Fellowship RF-2019-140\9, “Improving the Measurement of Time Using Radiocarbon”

Isotopic fingerprint for phosphorus in drinking water supplies

Phosphate dosing of drinking water supplies, coupled with leakage from distribution networks, represents a significant input of phosphorus to the environment. The oxygen isotope composition of phosphate (δ18OPO4), a novel stable isotope tracer for phosphorus, offers new opportunities to understand the importance of phosphorus derived from sources such as drinking water. We report the first assessment of δ18OPO4 within drinking water supplies. A total of 40 samples from phosphate-dosed distribution networks were analyzed from across England and Wales. In addition, samples of the source orthophosphoric acid used for dosing were also analyzed. Two distinct isotopic signatures for drinking water were identified (average = +13.2 or +19.7‰), primarily determined by δ18OPO4 of the source acid (average = +12.4 or +19.7‰). Dependent upon the source acid used, drinking water δ18OPO4 appears isotopically distinct from a number of other phosphorus sources. Isotopic offsets from the source acid ranging from −0.9 to +2.8‰ were observed. There was little evidence that equilibrium isotope fractionation dominated within the networks, with offsets from temperature-dependent equilibrium ranging from −4.8 to +4.2‰. While partial equilibrium fractionation may have occurred, kinetic effects associated with microbial uptake of phosphorus or abiotic sorption and dissolution reactions may also contribute to δ18OPO4 within drinking water supplies

Extended dilation of the radiocarbon time scale between 40,000 and 48,000 y BP and the overlap between Neanderthals and Homo sapiens

The new radiocarbon calibration curve (IntCal20) allows us to calculate the gradient of the relationship between 14C age and calendar age over the past 55 millennia before the present (55 ka BP). The new gradient curve exhibits a prolonged and prominent maximum between 48 and 40 ka BP during which the radiocarbon clock runs almost twice as fast as it should. This radiocarbon time dilation is due to the increase in the atmospheric 14C/12C ratio caused by the 14C production rise linked to the transition into the Laschamp geomagnetic excursion centered around 41 ka BP. The major maximum in the gradient from 48 to 40 ka BP is a new feature of the IntCal20 calibration curve, with far-reaching impacts for scientific communities, such as prehistory and paleoclimatology, relying on accurate ages in this time range. To illustrate, we consider the duration of the overlap between Neanderthals and Homo sapiens in Eurasia

SHCal13 Southern Hemisphere calibration, 0–50,000 years cal BP

The Southern Hemisphere SHCal04 radiocarbon calibration curve has been updated with the addition of new data sets extending measurements to 2145 cal BP and including the ANSTO Younger Dryas Huon pine data set. Outside the range of measured data, the curve is based upon the Northern Hemisphere data sets as presented in IntCal13, with an interhemispheric offset averaging 43 ± 23 yr modeled by an autoregressive process to represent the short-term correlations in the offset

IntCal09 and Marine09 radiocarbon age calibration curves, 0-50,000yeats cal BP

The IntCal04 and Marine04 radiocarbon calibration curves have been updated from 12 cal kBP (cal kBP is here defined as thousands of calibrated years before AD 1950), and extended to 50 cal kBP, utilizing newly available data sets that meet the IntCal Working Group criteria for pristine corals and other carbonates and for quantification of uncertainty in both the 14C and calendar timescales as established in 2002. No change was made to the curves from 0–12 cal kBP. The curves were constructed using a Markov chain Monte Carlo (MCMC) implementation of the random walk model used for IntCal04 and Marine04. The new curves were ratified at the 20th International Radiocarbon Conference in June 2009 and are available in the Supplemental Material at www.radiocarbon.org

Marine20—the marine radiocarbon age calibration curve (0 – 55,000 cal BP)

T.J. Heaton is supported by a Leverhulme Trust Fellowship RF-2019-140\9, “Improving the Measurement of Time Using Radiocarbon”. M Butzin is supported by the German Federal Ministry of Education and Research (BMBF), as Research for Sustainability initiative (FONA); www.fona.de through the PalMod project (grant numbers: 01LP1505B, 01LP1919A). E. Bard is supported by EQUIPEX ASTER-CEREGE and ANR CARBOTRYDH. Meetings of the IntCal Marine Focus group have been supported by Collège de France.The concentration of radiocarbon (14C) differs between ocean and atmosphere. Radiocarbon determinations from samples which obtained their 14C in the marine environment therefore need a marine-specific calibration curve and cannot be calibrated directly against the atmospheric-based IntCal20 curve. This paper presents Marine20, an update to the internationally agreed marine radiocarbon age calibration curve that provides a non-polar global-average marine record of radiocarbon from 0–55 cal kBP and serves as a baseline for regional oceanic variation. Marine20 is intended for calibration of marine radiocarbon samples from non-polar regions; it is not suitable for calibration in polar regions where variability in sea ice extent, ocean upwelling and air-sea gas exchange may have caused larger changes to concentrations of marine radiocarbon. The Marine20 curve is based upon 500 simulations with an ocean/atmosphere/biosphere box-model of the global carbon cycle that has been forced by posterior realizations of our Northern Hemispheric atmospheric IntCal20 14C curve and reconstructed changes in CO2 obtained from ice core data. These forcings enable us to incorporate carbon cycle dynamics and temporal changes in the atmospheric 14C level. The box-model simulations of the global-average marine radiocarbon reservoir age are similar to those of a more complex three-dimensional ocean general circulation model. However, simplicity and speed of the box model allow us to use a Monte Carlo approach to rigorously propagate the uncertainty in both the historic concentration of atmospheric 14C and other key parameters of the carbon cycle through to our final Marine20 calibration curve. This robust propagation of uncertainty is fundamental to providing reliable precision for the radiocarbon age calibration of marine based samples. We make a first step towards deconvolving the contributions of different processes to the total uncertainty; discuss the main differences of Marine20 from the previous age calibration curve Marine13; and identify the limitations of our approach together with key areas for further work. The updated values for ΔR, the regional marine radiocarbon reservoir age corrections required to calibrate against Marine20, can be found at the data base http://calib.org/marine/.Publisher PDFPeer reviewe

Gaur genome reveals expansion of sperm odorant receptors in domesticated cattle

publishedVersio

Marine20—The Marine Radiocarbon Age Calibration Curve (0–55,000 cal BP)

The concentration of radiocarbon (14C) differs between ocean and atmosphere. Radiocarbon determinations from samples which obtained their 14C in the marine environment therefore need a marine-specific calibration curve and cannot be calibrated directly against the atmospheric-based IntCal20 curve. This paper presents Marine20, an update to the internationally agreed marine radiocarbon age calibration curve that provides a non-polar global-average marine record of radiocarbon from 0–55 cal kBP and serves as a baseline for regional oceanic variation. Marine20 is intended for calibration of marine radiocarbon samples from non-polar regions; it is not suitable for calibration in polar regions where variability in sea ice extent, ocean upwelling and air-sea gas exchange may have caused larger changes to concentrations of marine radiocarbon. The Marine20 curve is based upon 500 simulations with an ocean/atmosphere/biosphere box-model of the global carbon cycle that has been forced by posterior realizations of our Northern Hemispheric atmospheric IntCal20 14C curve and reconstructed changes in CO2 obtained from ice core data. These forcings enable us to incorporate carbon cycle dynamics and temporal changes in the atmospheric 14C level. The box-model simulations of the global-average marine radiocarbon reservoir age are similar to those of a more complex three-dimensional ocean general circulation model. However, simplicity and speed of the box model allow us to use a Monte Carlo approach to rigorously propagate the uncertainty in both the historic concentration of atmospheric 14C and other key parameters of the carbon cycle through to our final Marine20 calibration curve. This robust propagation of uncertainty is fundamental to providing reliable precision for the radiocarbon age calibration of marine based samples. We make a first step towards deconvolving the contributions of different processes to the total uncertainty; discuss the main differences of Marine20 from the previous age calibration curve Marine13; and identify the limitations of our approach together with key areas for further work. The updated values for ΔR, the regional marine radiocarbon reservoir age corrections required to calibrate against Marine20, can be found at the data bas