Preprints
https://doi.org/10.5194/egusphere-2026-5845
https://doi.org/10.5194/egusphere-2026-5845
23 Oct 2026
 | 23 Oct 2026
Status: this preprint is open for discussion and under review for Geochronology (GChron).

A practical quantile approach to nonlinear uncertainty propagation in radiometric dating

Stefaan Pommé, Zsuzsanna Macsik, and Oana-Alexandra Dumitru

Abstract. Radiometric ages are commonly reported with a symmetric uncertainty (±u) obtained by first-order propagation. This is adequate when the age uncertainty is small relative to the characteristic timescale of the chronometer, but becomes unreliable when a measured isotope or activity ratio is transformed through a strongly nonlinear age equation. Equal positive and negative ratio deviations then produce unequal age deviations and asymmetric coverage intervals. Where the measured ratio approaches a physical bound, a finite upper age limit may cease to exist. Monte Carlo propagation can account for these effects, but is unnecessary for monotonic dating equations when the uncertainty budget can be adequately represented by a single equivalent perturbation of the measured ratio with standard deviation σeq. The degree of nonlinearity is then described by a single dimensionless parameter, β = ulin/τ, which compares the linearly propagated age uncertainty with the chronometer timescale (τ = 109–158 ka for 230Th/238U, 8.2 ka for 14C). Symmetric propagation is generally adequate for |β|<0.1. At larger |β|, asymmetric coverage limits follow from solving the dating equation at the measured ratio ±kσeq. Above |β|≈1/k, the ratio quantile near the phsycial bound reaches the limit of the age equation, so no finite upper coverage limit exists and the result must be reported as a lower age limit. The approach is developed for 230Th/238U dating and validated against Monte Carlo simulation. It is then extended to three families of monotonic chronometers: decaying parents, daughter/parent activity ratios approaching secular equilibrium, and accumulating daughter/parent atom ratios. Only the first two families have a dating limit, whereas classical geochronometers such as K–Ar, Rb–Sr, and U–Pb generally do not. Requiring only standard uncertainty budget parameters, the method can be implemented directly in a spreadsheet or data-reduction software.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Share
Stefaan Pommé, Zsuzsanna Macsik, and Oana-Alexandra Dumitru

Status: open (until 04 Dec 2026)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Stefaan Pommé, Zsuzsanna Macsik, and Oana-Alexandra Dumitru

Data sets

Spreadsheet workbook (Excel/VBA) for nonlinear uncertainty propagation in radiometric dating Stefaan Pommé https://doi.org/10.5281/zenodo.23021594

Model code and software

Spreadsheet workbook (Excel/VBA) for nonlinear uncertainty propagation in radiometric dating Stefaan Pommé https://doi.org/10.5281/zenodo.23021594

Stefaan Pommé, Zsuzsanna Macsik, and Oana-Alexandra Dumitru
Metrics will be available soon.
Latest update: 23 Oct 2026
Download
Short summary
Radiometric dating reveals the age of rocks, cave deposits, and other natural archives, but every age comes with an uncertainty. Because age does not grow in simple proportion to the quantity measured, symmetric error bars can misstate the result, especially for old or imprecise samples. We present a simple recipe, verified against computer simulations, that yields realistic asymmetric age ranges with ordinary spreadsheet tools and signals when a sample is too old to date reliably.
Share