Bayesian analysis of individual and systematic multiplicative errors for estimating ages with stratigraphic constraints in optically stimulated luminescence dating

Many dating techniques include significant error terms which are not independent between samples to date. This is typically the case in Optically Stimulated Luminescence (OSL) dating where the conversion from characteristic equivalent doses to the corresponding ages using the annual dosimetry data includes error terms that are common to all produced datings. Dealing with these errors is essential to estimate ages from a set of datings whose chronological ordering is known. In this work, we propose and we study a Bayesian model to address this problem. For this purpose, we first consider a multivariate model with multiplicative Gaussian errors in a Bayesian framework. This model relates a set of characteristic equivalent doses to the corresponding ages while taking into account for the systematic and non-systematic errors associated to the dosimetry. It thus offers the opportunity to deal properly with stratigraphic constraints within OSL datings, but also with other datings possessing errors which are independent from systematic errors of OSL (e.g. radiocarbon). Then, we use this model to extend an existing Bayesian model for the assessment of characteristic equivalent doses from Single Aliquot and Regenerative (SAR) dose measurements. The overall Bayesian model leads to the joint estimation of all the variables (which include all the dose-response functions and characteristic equivalent doses) of a sequence of, possibly heterogeneous, datings. We also consider a more generic solution consisting in using directly the age model from a set of characteristic equivalent dose estimates and their associated standard errors. We finally give an example of application on a set of five OSL datings with stratigraphic constraints and observe a good adequacy between the two approaches.