We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

Monte Carlo sampling for error propagation in linear regression and applications in isochron geochronology.

Li, Y. and Zhang, S. and Hobbs, R. and Caiado, C. and Sproson, A.D. and Selby, D. and Rooney, A.D. (2019) 'Monte Carlo sampling for error propagation in linear regression and applications in isochron geochronology.', Science bulletin., 64 (3). pp. 189-197.


Geochronology is essential for understanding Earth’s history. The availability of precise and accurate isotopic data is increasing; hence it is crucial to develop transparent and accessible data reduction techniques and tools to transform raw mass spectrometry data into robust chronological data. Here we present a Monte Carlo sampling approach to fully propagate uncertainties from linear regressions for isochron dating. Our new approach makes no prior assumption about the causes of variability in the derived chronological results and propagates uncertainties from both experimental measurements (analytical uncertainties) and underlying assumptions (model uncertainties) into the final age determination. Using synthetic examples, we find that although the estimates of the slope and y-intercept (hence age and initial isotopic ratios) are comparable between the Monte Carlo method and the benchmark “Isoplot” algorithm, uncertainties from the later could be underestimated by up to 60%, which are likely due to an incomplete propagation of model uncertainties. An additional advantage of the new method is its ability to integrate with geological information to yield refined chronological constraints. The new method presented here is specifically designed to fully propagate errors in linear regressions especially in geochronological applications involves linear regressions such as Rb-Sr, Sm-Nd, Re-Os, Pt-Os, Lu-Hf, U-Pb (with discordant points), Pb-Pb and Ar-Ar.

Item Type:Article
Full text:(AM) Accepted Manuscript
Available under License - Creative Commons Attribution Non-commercial No Derivatives.
Download PDF
Publisher Web site:
Publisher statement:© 2018 This manuscript version is made available under the CC-BY-NC-ND 4.0 license
Date accepted:26 November 2018
Date deposited:04 December 2018
Date of first online publication:22 December 2018
Date first made open access:22 December 2019

Save or Share this output

Look up in GoogleScholar