Abstract:Conventional (U-Th(-Sm))/He thermochronology requires accurate determination of the grain shape and dimension for calculation of 4He loss by alpha particle ejection and diffusion. In routine practice, calculations of 4He loss and diffusion are calculated assuming that measured grains approximate symmetric grain shapes with perfectly smooth crystal surfaces. This may introduce uncertainties that are virtually impossible to quantify if grains that do not pass the latter criteria are used. Alternatively, this leads to a selection bias for the analyses of detrital minerals if only grains with perfect crystal surfaces are measured. To address this limitation, we introduce a flexible, open-source numerical approach in Matlab™ to accurately determine the three-dimensional grain shape and numerically model He loss by alpha ejection and diffusion on individual grains. The three-dimensional surface of a mineral grain is interpolated from photomicrographs taken orthogonal to the principal grain axes. He ejection and diffusion is simulated using a coupled Monte Carlo and Brownian motion algorithm. Comparing our approach to high-resolution X-ray microtomography, we demonstrate that our “3D-He” approach improves the accuracy of He loss corrections when mineral grains do not approximate perfect crystal shapes by avoiding constraints inherent to standard manual-analytical approaches. Importantly, this method mitigates any bias introduced by user-dependent grain selection, measurements and assigning grains to a mineral shape for which analytical solutions of He production have been derived. In routine practice, we find that our 3D-He approach improves the accuracy of geometric measures (volume and surface area) and associated sphere-equivalent radii (SER) significantly, reducing the average deviation from true values from ~7–8% (manual-analytical approach) to ~2.5% (3D-He approach). The difference in SER significantly impacts the resulting apatite (U-Th-Sm)/He ages, which may deviate up to 40% when long held at temperatures where He production roughly matches diffusive loss of He. Using the 3D-He approach, the average deviation and overall spread of He loss corrections due to alpha ejection (FT factors) could be reduced from 1.6 ± 4.6% to 0.6 ± 0.8%. In summary, our 3D-He approach is a significant refinement to (U-Th(-Sm))/He age calculation considering relatively low analytical uncertainty in U, Th, Sm and He measurement (~2%). Importantly, our 3D-He approach should mitigate the selection bias introduced through manual grain selection. Toward this end, we discuss additional advantages of our numerical approach to handling broken grains or grains with known parent isotope zonation.

A Statistical Method for Treating the “Negative-Age” and “Negative-Atom-number” Conundrums in Radioisotope Dating

Statistical approach for radioactivity detection: A brief review

Insufficient Thermalization Effects on Determining Fission-Track Ages

Nuclear data evaluation for decay heat analysis of spent nuclear fuel over 1–100 k year timescale

A Bayesian hierarchical age model for single-grain optical dating of feldspars

A novel method of identification of in-ring decay and its application in the half-life estimates of 94mRu44

A Modified Method to Calculate Dual-Isotope Slopes for the Natural Attenuation of Organic Pollutants in the Environment

Enhancing artificial radionuclides monitoring: a Bayesian statistical approach combined with the multi time-interval analysis method

Analysis and packaging of radiochemical solar neutrino data. 1. Bayesian approach

Estimation of Ru-97 Half-Life Using the Most Frequent Value Method and Bootstrapping Techniques

On the biases affecting water ages inferred from isotopic data

Neutron Counting Statistics calculations Using Deterministic Transport

Radioactive Methods for Determining Geological Age

Evaluation of K-Ar and (40)air/(39) Ar Dating and Q Value

Bayesian analysis of De distributions in optical dating: Towards a robust method for dealing with outliers

Decay/Ingrowth Uncertainty Correction Of Po-210/Pb-210 In Seawater

Increasing the accuracy of (U-Th(-Sm))/He dating with 3D grain modelling

Optimization of analytical conditions for precise and accurate isotope analyses of Li, Mg, Fe, Cu, and Zn by MC-ICPMS

Systematic Calculations on Cluster Radioactivity Half-Lives

The Issues of Application for the Fission Track Dating and Its Geological Significance

Quantifying the age peaks, age ranges and weights of detrital ages based on the EM algorithm