Articles | Volume 5, issue 1
https://doi.org/10.5194/gchron-5-263-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gchron-5-263-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Short communication: The Wasserstein distance as a dissimilarity metric for comparing detrital age spectra and other geological distributions
Merton College, University of Oxford, Oxford, UK
Pieter Vermeesch
Department of Earth Sciences, University College London, London, UK
Related authors
No articles found.
Pieter Vermeesch, Tim Band, Jiangping He, Rex Galbraith, and Andrew Carter
EGUsphere, https://doi.org/10.5194/egusphere-2025-4948, https://doi.org/10.5194/egusphere-2025-4948, 2025
This preprint is open for discussion and under review for Geochronology (GChron).
Short summary
Short summary
geochron@home is an open-source platform that makes fission track dating more transparent and reliable. It combines a virtual microscope with an online database to share images and data openly, following FAIR principles. Researchers can analyse tracks privately, archive data for peer review, teach students, or involve citizen scientists. By improving data access and reproducibility, geochron@home helps build trust and supports future advances in Earth science.
Pieter Vermeesch, Noah McLean, Anton Vaks, Tzahi Golan, Sebastian F. M. Breitenbach, and Randall Parrish
Geochronology, 7, 459–473, https://doi.org/10.5194/gchron-7-459-2025, https://doi.org/10.5194/gchron-7-459-2025, 2025
Short summary
Short summary
U–Pb dating of cave sediments has provided important new time constraints on the evolution of cave-dwelling organisms (including early humans) and of Earth's climate during the past 5 Myr. This paper shows that the most common type of U–Pb dating, which uses 238U and 206Pb, can be inaccurate beyond ca. 2 Myr ago. It proposes an alternative type of U–Pb dating, using 235U and 207Pb, as a more accurate alternative.
Pieter Vermeesch
Geochronology, 6, 397–407, https://doi.org/10.5194/gchron-6-397-2024, https://doi.org/10.5194/gchron-6-397-2024, 2024
Short summary
Short summary
The age of some geological materials can be estimated from the ratio of certain radiogenic "daughter" isotopes to their radioactive "parent". However, in many cases, the age estimation process is complicated by the presence of an inherited component of non-radiogenic daughter isotopes. This paper presents an improved algorithm to estimate the radiogenic and non-radiogenic components, either separately or jointly.
Pieter Vermeesch, Yuntao Tian, Jae Schwanethal, and Yannick Buret
Geochronology, 5, 323–332, https://doi.org/10.5194/gchron-5-323-2023, https://doi.org/10.5194/gchron-5-323-2023, 2023
Short summary
Short summary
The U–Th–He method is a technique to determine the cooling history of minerals. Traditional approaches to U–Th–He dating are time-consuming and require handling strong acids and radioactive solutions. This paper presents an alternative approach in which samples are irradiated with protons and subsequently analysed by laser ablation mass spectrometry. Unlike previous in situ U–Th–He dating attempts, the new method does not require any absolute concentration measurements of U, Th, or He.
Matthew Fox, Adam G. G. Smith, Pieter Vermeesch, Kerry Gallagher, and Andrew Carter
Geochronology Discuss., https://doi.org/10.5194/gchron-2022-23, https://doi.org/10.5194/gchron-2022-23, 2022
Publication in GChron not foreseen
Short summary
Short summary
The Great Unconformity represents an enormous amount of time lost from the sedimentary record. Its origin is debated, in part, due to different approaches used to interpret zircon (U–Th)/He ages. This thermochronometric system is ideal for this problem because the temperature sensitivity varies according to radiation damage. Here we explore the uncertainty associated with the radiation damage model and show how this limits our ability to resolve the origin of the Great Unconformity.
Pieter Vermeesch
Geochronology, 4, 561–576, https://doi.org/10.5194/gchron-4-561-2022, https://doi.org/10.5194/gchron-4-561-2022, 2022
Short summary
Short summary
Secondary ion mass spectrometry (SIMS) is the oldest and most sensitive analytical technique for in situ U–Pb geochronology. This paper introduces a new algorithm for SIMS data reduction that treats data as
compositional data, which means that the relative abundances of 204Pb, 206Pb, 207Pb, and 238Pb are processed within a tetrahedral data space or
simplex. The new method is implemented in an eponymous computer programme that is compatible with the two dominant types of SIMS instruments.
Yang Li and Pieter Vermeesch
Geochronology, 3, 415–420, https://doi.org/10.5194/gchron-3-415-2021, https://doi.org/10.5194/gchron-3-415-2021, 2021
Short summary
Short summary
A conventional isochron is a straight-line fit to two sets of isotopic ratios, D/d and P/d, where P is the radioactive parent, D is the radiogenic daughter, and d is a second isotope of the daughter element. The slope of this line is proportional to the age of the system. An inverse isochron is a linear fit through d/D and P/D. The horizontal intercept of this line is inversely proportional to the age. The latter approach is preferred when d<D, which is the case in Re–Os and K–Ca geochronology.
Pieter Vermeesch
Geochronology, 3, 247–257, https://doi.org/10.5194/gchron-3-247-2021, https://doi.org/10.5194/gchron-3-247-2021, 2021
Short summary
Short summary
This paper shows that the current practice of filtering discordant U–Pb data based on the relative difference between the 206Pb/238U and 207Pb/206Pb ages is just one of several possible approaches to the problem and demonstrably not the best one. An alternative approach is to define discordance in terms of isotopic composition, as a log ratio distance between the measurement and the concordia line. Application to real data indicates that this reduces the positive bias of filtered age spectra.
Cited articles
Amidon, W. H., Burbank, D. W., and Gehrels, G. E.: Construction of detrital
mineral populations: insights from mixing of U–Pb zircon ages in
Himalayan rivers, Basin Res., 17, 463–485,
https://doi.org/10.1111/j.1365-2117.2005.00279.x, 2005. a
Benamou, J.-D., Carlier, G., Cuturi, M., Nenna, L., and Peyré, G.: Iterative
Bregman Projections for Regularized Transportation Problems, SIAM
J. Sci. Comput., 2, A1111–A1138, https://doi.org/10.1137/141000439,
2015. a
Berry, R. F., Jenner, G. A., Meffre, S., and Tubrett, M. N.: A North
American provenance for Neoproterozoic to Cambrian sandstones in
Tasmania?, Earth Planet. Sc. Lett., 192, 207–222,
https://doi.org/10.1016/S0012-821X(01)00436-8, 2001. a
Cawood, P., Hawkesworth, C., and Dhuime, B.: Detrital zircon record and
tectonic setting, Geology, 40, 875–878, https://doi.org/10.1130/G32945.1, 2012. a
Condie, K. C., Belousova, E., Griffin, W. L., and Sircombe, K. N.: Granitoid
events in space and time: Constraints from igneous and detrital zircon age
spectra, Gondwana Res., 15, 228–242, https://doi.org/10.1016/j.gr.2008.06.001,
2009. a
De Doncker, F., Herman, F., and Fox, M.: Inversion of provenance data and
sediment load into spatially varying erosion rates, Earth Surf. Proc.
Land., 45, 3879–3901, https://doi.org/10.1002/esp.5008, 2020. a
DeGraaff-Surpless, K., Graham, S. A., Wooden, J. L., and McWilliams, M. O.:
Detrital zircon provenance analysis of the Great Valley Group,
California: Evolution of an arc-forearc system, GSA Bulletin, 114,
1564–1580, https://doi.org/10.1130/0016-7606(2002)114<1564:DZPAOT>2.0.CO;2, 2002. a, b, c
Dietze, E., and Dietze, M.: Grain-size distribution unmixing using the R package EMMAgeo, E&G Quaternary Sci. J., 68, 29–46, https://doi.org/10.5194/egqsj-68-29-2019, 2019. a
Engquist, B. and Froese, B. D.: Application of the Wasserstein metric to
seismic signals, Commun. Math. Sci., 12, 979–988,
https://doi.org/10.4310/CMS.2014.v12.n5.a7, 2014. a
Flamary, R., Courty, N., Gramfort, A., Alaya, M. Z., Boisbunon, A., Chambon,
S., Chapel, L., Corenflos, A., Fatras, K., Fournier, N., Gautheron, L.,
Gayraud, N. T. H., Janati, H., Rakotomamonjy, A., Redko, I., Rolet, A.,
Schutz, A., Seguy, V., Sutherland, D. J., Tavenard, R., Tong, A., and Vayer,
T.: POT: Python Optimal Transport, J. Mach. Learn.
Res., 22, 1–8, 2021. a
Irpino, A. and Romano, E.: Optimal histogram representation of large data sets:
Fisher vs piecewise linear approximation, in: Actes des cinquièmes
journées Extraction et Gestion des Connaissances, edited by: Noirhomme-Fraiture, M. and Venturini, G., Vol. E-9,
99–110, Namur, Belgium, https://editions-rnti.fr/?inprocid=1001314 (last access: 15 May 2023), 2007. a
Košler, J., Sláma, J., Belousova, E., Corfu, F., Gehrels, G. E., Gerdes, A.,
Horstwood, M. S. A., Sircombe, K. N., Sylvester, P. J., Tiepolo, M.,
Whitehouse, M. J., and Woodhead, J. D.: U-Pb Detrital Zircon Analysis
– Results of an Inter-laboratory Comparison, Geostand.
Geoanal. Res., 37, 243–259, https://doi.org/10.1111/j.1751-908X.2013.00245.x,
2013. a, b
Lipp, A.: AlexLipp/detrital-wasserstein: Acceptance (1.0), Zenodo [code, data set], https://doi.org/10.5281/zenodo.7937484, 2023. a
Magyar, J. C. and Sambridge, M.: Hydrological objective functions and ensemble averaging with the Wasserstein distance, Hydrol. Earth Syst. Sci., 27, 991–1010, https://doi.org/10.5194/hess-27-991-2023, 2023. a
Métivier, L., Brossier, R., Mérigot, Q., Oudet, E., and Virieux, J.: An
optimal transport approach for seismic tomography: application to 3D full
waveform inversion, Inverse Probl., 32, 115008,
https://doi.org/10.1088/0266-5611/32/11/115008, 2016. a
Morton, A., Fanning, M., and Milner, P.: Provenance characteristics of
Scandinavian basement terrains: Constraints from detrital zircon ages in
modern river sediments, Sediment. Geol., 210, 61–85,
https://doi.org/10.1016/j.sedgeo.2008.07.001, 2008. a
Peyré, G. and Cuturi, M.: Computational Optimal Transport, Foundations and
Trends in Machine Learning, 11, 355–607, 2019. a
Reimink, J. R., Davies, J. H. F. L., and Ielpi, A.: Global zircon analysis
records a gradual rise of continental crust throughout the Neoarchean,
Earth Planet. Sc. Lett., 554, 116654,
https://doi.org/10.1016/j.epsl.2020.116654, 2021. a
Sambridge, M., Jackson, A., and Valentine, A. P.: Geophysical inversion and
optimal transport, Geophys. J. Int., 231, 172–198,
https://doi.org/10.1093/gji/ggac151, 2022. a, b
Satkoski, A. M., Wilkinson, B. H., Hietpas, J., and Samson, S. D.: Likeness
among detrital zircon populations – An approach to the comparison of age
frequency data in time and space, GSA Bulletin, 125, 1783–1799,
https://doi.org/10.1130/B30888.1, 2013. a
Saylor, J., Stockli, D., Horton, B., Nie, J., and Mora, A.: Discriminating
rapid exhumation from syndepositional volcanism using detrital zircon double
dating: Implications for the tectonic history of the Eastern
Cordillera, Colombia, B. Geol. Soc. Am., 124,
762–779, https://doi.org/10.1130/B30534.1, 2012. a
Saylor, J. E. and Sundell, K. E.: Quantifying comparison of large detrital
geochronology data sets, Geosphere, 12, 203–220, https://doi.org/10.1130/GES01237.1,
2016. a
Schuhmacher, D., Bähre, B., Gottschlich, C., Hartmann, V., Heinemann, F., and
Schmitzer, B.: transport: Computation of Optimal Transport Plans and
Wasserstein Distances,
https://cran.r-project.org/package=transport (last access: 15 May 2023), 2022. a
Sharman, G. R. and Johnstone, S. A.: Sediment unmixing using detrital
geochronology, Earth Planet. Sc. Lett., 477, 183–194,
https://doi.org/10.1016/j.epsl.2017.07.044, 2017. a
Sharman, G. R., Sharman, J. P., and Sylvester, Z.: detritalPy: A
Python-based toolset for visualizing and analysing detrital
geo-thermochronologic data, The Depositional Record, 4, 202–215,
https://doi.org/10.1002/dep2.45, 2018. a
Sundell, K. E. and Saylor, J. E.: Two-Dimensional Quantitative Comparison
of Density Distributions in Detrital Geochronology and
Geochemistry, Geochem. Geophy. Geosy., 22, e2020GC009559,
https://doi.org/10.1029/2020GC009559, 2021. a
Vermeesch, P.: Multi-sample comparison of detrital age distributions, Chem.
Geol., 341, 140–146, https://doi.org/10.1016/j.chemgeo.2013.01.010, 2013. a
Vermeesch, P.: Dissimilarity measures in detrital geochronology, Earth-Sci.
Rev., 178, 310–321, https://doi.org/10.1016/j.earscirev.2017.11.027,
2018a. a, b, c
Vermeesch, P.: IsoplotR: A free and open toolbox for geochronology,
Geosci. Front., 9, 1479–1493, https://doi.org/10.1016/j.gsf.2018.04.001,
2018b. a
Vermeesch, P., Lipp, A. G., Hatzenbühler, D., Caracciolo, L., and Chew, D.:
Multidimensional scaling of varietal data in sedimentary provenance analysis,
J. Geophys. Res.-Earth, 128, e2022JF006992,
https://doi.org/10.1029/2022JF006992, 2023. a
Weltje, G. J.: End-member modeling of compositional data:
Numerical-statistical algorithms for solving the explicit mixing problem,
Math. Geol., 29, 503–549, https://doi.org/10.1007/BF02775085, 1997. a
Wobus, C. W., Hodges, K. V., and Whipple, K. X.: Has focused denudation
sustained active thrusting at the Himalayan topographic front?, Geology,
31, 861–864, https://doi.org/10.1130/G19730.1, 2003. a, b, c
Short summary
We propose using the Wasserstein-2 distance (W2) as an alternative to the widely used Kolmogorov–Smirnov (KS) statistic for analysing distributional data in geochronology. W2 measures the horizontal distance between observations, while KS measures vertical differences in cumulative distributions. Using case studies, we find that W2 is preferable in scenarios where the absolute age differences in observations provide important geological information. W2 has been added to the R package IsoplotR.
We propose using the Wasserstein-2 distance (W2) as an alternative to the widely used...