Technical note: LA–ICP-MS U–Pb dating of unetched and etched apatites
- 1CONACyT–Centro de Geociencias, Campus Juriquilla, UNAM, Querétaro, 76230, Mexico
- 2Centro de Geociencias, Campus Juriquilla, UNAM, Querétaro, 76230, Mexico
- 3LANGEM, Instituto de Geología, UNAM, Ciudad Universitaria, Mexico City, 04510, Mexico
Correspondence: Fanis Abdullin (firstname.lastname@example.org)
The same unetched and chemically etched apatite crystals from five rock samples were dated by the U–Pb method via laser ablation inductively coupled plasma mass spectrometry (LA–ICP-MS). The objective of this study is to test whether chemical etching required for apatite fission track analysis impacts the precision and accuracy of apatite U–Pb geochronology. The results of this experiment suggest that etching has insignificant effects on the accuracy of apatite U–Pb ages obtained by LA–ICP-MS. Therefore, LA–ICP-MS is reliable for U–Pb analysis as part of apatite fission track and U–Pb double dating.
Apatite, Ca5(PO4)3[F,Cl,OH], is the most common phosphate mineral in the Earth's crust and can be found in practically all igneous and metamorphic rocks, in many ancient and recent sediments as well as in certain mineral deposits (Piccoli and Candela, 2002; Morton and Yaxley, 2007; Webster and Piccoli, 2015). This accessory mineral is often used as a natural thermochronometer for fission track, helium, U–Th, and U–Pb dating (e.g., Zeitler et al., 1987; Wolf et al., 1996; Ehlers and Farley, 2003; Hasebe et al., 2004; Donelick et al., 2005; Chew and Donelick, 2012; Chew et al., 2014; Cochrane et al., 2014; Liu et al., 2014; Spikings et al., 2015; Glorie et al., 2017). Presently, apatite fission track (AFT) ages can be obtained rapidly by using laser ablation inductively coupled plasma mass spectrometry (LA–ICP-MS) for direct measurement of “parent nuclides”, i.e., 238U contents (Cox et al., 2000; Svojtka and Košler, 2002; Hasebe et al., 2004, 2009; Donelick et al., 2005; Abdullin et al., 2014, 2016, 2018; Vermeesch, 2017). The LA–ICP-MS technique may be used to measure 238U for AFT dating, together with Pb isotopes needed for U–Pb dating (e.g., Chew and Donelick, 2012; Liu et al., 2014; Glorie et al., 2017; Bonilla et al., 2020; Nieto-Samaniego et al., 2020).
Hasebe et al. (2009) previously performed an important experimental study, during which they demonstrated that chemical etching required for apatite and zircon fission track dating does not interfere with U analysis by LA–ICP-MS. The influence of etching needed for AFT dating on the precision and accuracy of dating the same crystals by U–Pb using LA–ICP-MS remains to be quantified. To investigate this issue, the same unetched and etched apatite grains extracted from five rock samples were analyzed via LA–ICP-MS for U–Pb dating. The chosen samples have either emplacement or metamorphic ages ranging from the Cretaceous to the Neoproterozoic (see Table 1 for further details).
2.1 OV-0421 (Tres Sabanas Pluton, Guatemala)
This sample is a two-mica-bearing deformed granite belonging to the Tres Sabanas Pluton, which is located northwest of Guatemala City, Guatemala. For sample OV-0421, an emplacement age of 115±4 (2σ) Ma was proposed based on zircon U–Pb data (Torres de León, 2016). A cooling age of 102±1 (2σ) Ma, obtained with K–Ar (on biotite), was also reported by the same author.
2.2 MCH-38 (Chiapas Massif Complex, Mexico)
MCH-38 is an orthogneiss from the Permian Chiapas Massif Complex. This rock was sampled to the west of Unión Agrarista, the State of Chiapas, southeastern Mexico. There is no reported age for this sample. Some zircon U–Pb dates obtained for the Chiapas Massif Complex (Weber et al., 2007, 2008; Ortega-Obregón et al., 2019) suggest that a Lopingian (260–252 Ma) crystallization or metamorphic age may be assumed for sample MCH-38.
2.3 TO-AM (Totoltepec Pluton, Mexico)
TO-AM is a granitic rock, sampled ca. 5 km west of Totoltepec de Guerrero, the State of Puebla, southern Mexico. There is no reported radiometric data for sample TO-AM. Previous geological studies indicate that the Pennsylvanian–Cisuralian Totoltepec Pluton was emplaced over a ca. 23 million year period (from ca. 308 to ca. 285 Ma; e.g., Kirsch et al., 2013).
2.4 CH-0403 (Altos Cuchumatanes, Guatemala)
CH-0403 was collected 5 km ESE of Barillas, in Altos Cuchumatanes, Guatemala. It consists of a gray to green granodiorite. Five zircon aliquots of sample CH-0403 were dated using isotope dilution thermal-ionization mass spectrometry, yielding a lower intercept date of 391±8 (2σ) Ma that is interpreted as its approximate crystallization age (Solari et al., 2009).
2.5 OC-1008 (Oaxacan Complex, Mexico)
This sample is a paragneiss from the Grenvillian Oaxacan Complex, southern Mexico. OC-1008 was collected in the federal road which connects Nochixtlán to Oaxaca. It was demonstrated that this sample underwent granulite facies metamorphism at 1000–980 Ma (Solari et al., 2014).
Accessory minerals were concentrated using conventional mineral separation techniques such as rock crushing, sieving, Wilfley table, Frantz magnetic separator, and bromoform. Approximately 300 apatite grains were extracted from each rock sample and mounted with their surfaces parallel to the crystallographic c axis in a 2.5 cm diameter epoxy mount. Mounted crystals were polished to expose their internal surfaces (i.e., up to 4π geometry). For this experiment, complete crystals lacking visible inclusions and other defects, such as cracks, were carefully selected for analysis. Sample preparation was performed at Taller de Molienda and Taller de Laminación, Centro de Geociencias (CGEO), Campus Juriquilla, Universidad Nacional Autónoma de Mexico (UNAM).
Single spot analyses were performed with a Resonetics RESOlution™ LPX Pro (193 nm, ArF excimer) laser ablation system, coupled to a Thermo Scientific iCAP™ Qc quadrupole ICP-MS at Laboratorio de Estudios Isotópicos (LEI), CGEO, UNAM. During this experimental work, LA–ICP-MS-based sampling was performed in central parts of the selected apatite grains before and after chemical etching (in 5.5 M HNO3 at 21 ∘C for 20 s to reveal spontaneous fission tracks), as shown schematically in Fig. 1. The LA–ICP-MS protocol used for apatite analyses, as given in Table 2, was established on the basis of numerous experiments carried out at LEI during the past 5 years and can be used for U–Pb and fission track double dating plus multielemental analysis (Abdullin et al., 2018; Ortega-Obregón et al., 2019). Corrected isotopic ratios and errors were calculated using Iolite 3.5 (Paton et al., 2011) and the VizualAge data reduction scheme (Petrus and Kamber, 2012). UcomPbine (Chew et al., 2014) was used to model 207Pb ∕ 206Pb initial values and thus force a 207Pb correction that considers the common Pb (non-radiogenic Pb) incorporated by apatite standards at the moment of their crystallization (see also Ortega-Obregón et al., 2019). The “First Mine Discovery” apatite from Madagascar, with a mean U–Pb age of ca. 480 Ma (Thomson et al., 2012; Chew et al., 2014), was used as a primary reference material. The results for measured isotopes using NIST-612 (Pearce et al., 1997) were normalized using 43Ca as an internal standard and taking an average CaO content of 55 %.
Tera–Wasserburg Concordia diagrams (T–W; Tera and Wasserburg, 1972) are used in apatite U–Pb dating, because the LA–ICP-MS-derived U–Pb results are generally discordant. The lower intercept in the T–W plot is considered a mean apatite U–Pb age that should have geological significance (crystallization or cooling age, the age of mineralization or metamorphic event). Apatite U–Pb ages were calculated with IsoplotR (Vermeesch, 2017, 2018) and described below. Detailed information on U–Pb experiments is given in Table S1 in the Supplement.
Note: MAD – “First Mine Discovery” U–Pb apatite standard from Madagascar; DUR – Durango apatite from Cerro de Mercado mine (Mexico); apt – unknown apatite crystals. * Laser pulse energy of 4 J cm−2, which was measured directly on target with a Coherent™ laser energy meter.
For rock sample OV-0421, 41 unetched apatites yielded a lower intercept age of 106±4 (2σ) Ma with a mean square weighted deviation (MSWD) of 1.07, passing the chi-squared test with the P(χ2) value of 0.35 (see in Fig. 2). Practically the same U–Pb date, 107±5 (2σ) Ma, was obtained after chemical etching of the same apatite grains, yielding a MSWD of 1.13 and a P(χ2) of 0.27. Both these apatite U–Pb ages lie between the zircon U–Pb date of 115±4 (2σ) Ma (i.e., crystallization age) and the biotite K–Ar age of 102±1 (2σ) Ma (i.e., cooling age), which were previously obtained for the same granite sample by Torres de León (2016).
For orthogneiss sample MCH-38, the lower intercept in T–W yielded a U–Pb age of 245±6 (2σ) Ma (obtained from 41 unetched apatites) with a MSWD of 0.28 and a P(χ2) of 1. Etched apatite grains from MCH-38 yielded an age of 240±4 (2σ) Ma with a MSWD of 0.36 and a P(χ2) of 1 (Fig. 2). Our U–Pb results are in close agreement with geochronological data reported from the Chiapas Massif Complex in previous studies (Damon et al., 1981; Torres et al., 1999; Schaaf et al., 2002; Ortega-Obregón et al., 2019). For instance, Torres et al. (1999) compiled biotite K–Ar ages, most of which lie within Early–Middle Triassic period. Triassic cooling ages in the Chiapas Massif Complex were also detected by Rb–Sr in mica–whole rock pairs that range from 244±12 (2σ) to 214±11 (2σ) Ma (Schaaf et al., 2002).
Unetched apatites (32 crystals; Fig. 2) from granite TO-AM yielded a lower intercept date of 303±5 (2σ) Ma with a MSWD of 0.6 and a P(χ2) of 0.96. After etching, a slightly younger age of 299±3 (2σ) Ma was obtained, with a MSWD of 0.89 and a P(χ2) of 0.65. These apatite U–Pb ages are in line with the zircon U–Pb ages of 306±2 (2σ) Ma to 287±2 (2σ) Ma reported for the Pennsylvanian–Cisuralian Totoltepec Pluton (e.g., see details in Kirsch et al., 2013).
A total of 36 unetched apatite grains from sample CH-0403 yielded a lower intercept U–Pb age of 345±10 (2σ) Ma with a MSWD of 0.7 and a P(χ2) of 0.9, whereas etched grains yielded an age of 334±8 (2σ) Ma with a MSWD of 1.37 and a P(χ2) of 0.08 (Fig. 2). These cooling dates are considerably younger if compared to the CH-0403 emplacement age of 391±8 (2σ) Ma (Solari et al., 2009).
A total of 41 unetched apatites belonging to sample OC-1008 yielded a U–Pb age of 839±12 (2σ) Ma with a MSWD of 0.98 and a P(χ2) of 0.50. After etching, the same apatite crystals yielded an age of 830±10 (2σ) Ma with a MSWD of 1.24 and a P(χ2) of 0.14 (Fig. 2). Both these apatite U–Pb ages are significantly younger than the age of granulite facies metamorphism in the Grenville-aged Oaxacan Complex (1 Ga to 980 Ma, Solari et al., 2014) and, thus, should be considered as cooling ages.
Most rock samples, except OV-0421, yielded slightly younger apatite U–Pb ages after chemical etching (up to 3.3 % in sample CH-0403). However, the lower intercept U–Pb ages obtained from unetched apatite grains are indistinguishable within error from the U–Pb ages obtained on the same etched grains (see diagram in Fig. 3). The results of this experiment demonstrate that chemical etching required for AFT analysis has negligible effects on the accuracy of apatite U–Pb ages determined via LA–ICP-MS. Thus, as a main conclusion of this study, LA–ICP-MS can be used for simultaneous AFT and U–Pb double dating, as it was already done in some previous studies (e.g., Chew and Donelick, 2012; Liu et al., 2014; Glorie et al., 2017; Bonilla et al., 2020; Nieto-Samaniego et al., 2020).
The authors declare that all the data supporting the findings of this study are available within the article (see Supplement).
The supplement related to this article is available online at: https://doi.org/10.5194/gchron-3-59-2021-supplement.
Conceptualization, investigation, and writing of the original draft were done by FA. LS and COO provided technical support. LS and JS acquired funding and resources, supervised the study, and reviewed the manuscript.
The authors declare that they have no conflict of interest.
The authors are grateful to Juan Tomás Vázquez Ramírez and Ofelia Pérez Arvizu for their help with sample preparation for this study. Stuart Thomson is acknowledged for sharing the Madagascar apatite. Michelangelo Martini kindly provided the sample TO-AM that was useful for our experimental study. Ziva Shulaker, Jakub Sliwinski, and Axel Schmitt are acknowledged for their constructive comments that improved our manuscript significantly.
This research has been supported by the PAPIIT DGAPA UNAM (grant no. IN101520).
This paper was edited by Axel Schmitt and reviewed by Jakub Sliwinski and Ziva Shulaker.
Abdullin, F., Solé, J., and Solari, L.: Datación mediante trazas de fisión y análisis multielemental con LA-ICP-MS del fluorapatito de Cerro de Mercado (Durango, México), Revista Mexicana de Ciencias Geológicas, 31, 395–406, 2014.
Abdullin, F., Solé, J., Meneses-Rocha, J. D. J., Solari, L., Shchepetilnikova, V., and Ortega-Obregón, C.: LA-ICP-MS-based apatite fission track dating of the Todos Santos Formation sandstones from the Sierra de Chiapas (SE Mexico) and its tectonic significance, Int. Geol. Rev., 58, 32–48, https://doi.org/10.1080/00206814.2015.1055596, 2016.
Abdullin, F., Solari, L., Ortega-Obregón, C., and Solé, J.: New fission-track results from the northern Chiapas Massif area, SE Mexico: trying to reconstruct its complex thermo-tectonic history, Revista Mexicana de Ciencias Geológicas, 35, 79–92, https://doi.org/10.22201/cgeo.20072902e.2018.1.523, 2018.
Bonilla, A., Franco, J. A., Cramer, T., Poujol, M., Cogné, N., Nachtergaele, S., and De Grave, J.: Apatite LA-ICP-MS U–Pb and fission-track geochronology of the Caño Viejita gabbro in E-Colombia: Evidence for Grenvillian intraplate rifting and Jurassic exhumation in the NW Amazonian Craton, J. S. Am. Earth Sci., 98, 102438, https://doi.org/10.1016/j.jsames.2019.102438, 2020.
Chew, D. M. and Donelick, R. A.: Combined apatite fission track and U-Pb dating by LA-ICP-MS and its application in apatite provenance analysis, Mineralogical Association of Canada, Short Course, 42, 219–247, 2012.
Chew, D. M., Petrus, J. A., and Kamber, B. S.: U–Pb LA–ICPMS dating using accessory mineral standards with variable common Pb, Chemical Geology, 363, 185–199, https://doi.org/10.1016/j.chemgeo.2013.11.006, 2014.
Cochrane, R., Spikings, R. A., Chew, D., Wotzlaw, J. F., Chiaradia, M., Tyrrell, S., Schaltegger, U., and Van der Lelij, R.: High temperature (> 350 ∘C) thermochronology and mechanisms of Pb loss in apatite, Geochim. Cosmochim. Acta, 127, 39–56, https://doi.org/10.1016/j.gca.2013.11.028, 2014.
Cox, R., Košler, J., Sylvester, P., and Hodych, P.: Apatite fission-track (FT) dating by LAM-ICP-MS analysis, Goldschmidt Conference, Oxford, UK, Journal of Conference Abstracts, 5, 322, 2000.
Damon, P. E., Shafiqullah, M., and Clark, K. F.: Age trends of igneous activity in relation to metallogenesis in the southern Cordillera, Tucson, Arizona, Arizona Geological Society Digest, 14, 137–153, 1981.
Donelick, R. A., O'Sullivan, P. B., and Ketcham, R. A.: Apatite fission-track analysis, Rev. Mineral. Geochem., 58, 49–94, https://doi.org/10.2138/rmg.2005.58.3, 2005.
Ehlers, T. A. and Farley, K. A.: Apatite (U–Th)/He thermochronometry: methods and applications to problems in tectonic and surface processes, Earth Planet. Sci. Lett., 206, 1–14, https://doi.org/10.1016/S0012-821X(02)01069-5, 2003.
Glorie, S., Alexandrov, I., Nixon, A., Jepson, G., Gillespie, J., and Jahn, B. M.: Thermal and exhumation history of Sakhalin Island (Russia) constrained by apatite U-Pb and fission track thermochronology, J. Asian Earth Sci., 143, 326–342, https://doi.org/10.1016/j.jseaes.2017.05.011, 2017.
Hasebe, N., Barbarand, J., Jarvis, K., Carter, A., and Hurford, A. J.: Apatite fission-track chronometry using laser ablation ICP-MS, Chem. Geol., 207, 135–145, https://doi.org/10.1016/j.chemgeo.2004.01.007, 2004.
Hasebe, N., Carter, A., Hurford, A. J., and Arai, S.: The effect of chemical etching on LA–ICP-MS analysis in determining uranium concentration for fission-track chronometry, Geol. Soc. Spec. Publ., 324, 37–46, https://doi.org/10.1144/SP324.3, 2009.
Kirsch, M., Keppie, J. D., Murphy, J. B., and Lee, J. K.: Arc plutonism in a transtensional regime: the late Palaeozoic Totoltepec pluton, Acatlán Complex, southern Mexico, Int. Geol. Rev., 55, 263–286, https://doi.org/10.1080/00206814.2012.693247, 2013.
Liu, W., Zhang, J., Sun, T., and Wang, J.: Application of apatite U–Pb and fission-track double dating to determine the preservation potential of magnetite–apatite deposits in the Luzong and Ningwu volcanic basins, eastern China, J. Geochem. Explor., 138, 22–32, https://doi.org/10.1016/j.gexplo.2013.12.006, 2014.
Morton, A. and Yaxley, G.: Detrital apatite geochemistry and its application in provenance studies, Special Papers, Geol. Soc. America, 420, 319–344, https://doi.org/10.1130/2006.2420(19), 2007.
Nieto-Samaniego, A. F., Olmos-Moya, M. D. J. P., Levresse, G., Alaniz-Alvarez, S. A., Abdullin, F., del Pilar-Martínez, A., and Xu, S.: Thermochronology and exhumation rates of granitic intrusions at Mesa Central, Mexico, Int. Geol. Rev., 62, 311–319, https://doi.org/10.1080/00206814.2019.1602789, 2020.
Ortega-Obregón, C., Abdullin, F., Solari, L., Schaaf, P., and Solís-Pichardo, G.: Apatite U-Pb dating at UNAM laboratories: analytical protocols and examples of its application, Revista Mexicana de Ciencias Geológicas, 36, 27–37, https://doi.org/10.22201/cgeo.20072902e.2019.1.749, 2019.
Paton, C., Hellstrom, J., Paul, B., Woodhead, J., and Hergt, J.: Iolite: Freeware for the visualisation and processing of mass spectrometric data, J. Anal. Atom. Spectrom., 26, 2508–2518, https://doi.org/10.1039/C1JA10172B, 2011.
Pearce, N. J., Perkins, W. T., Westgate, J. A., Gorton, M. P., Jackson, S. E., Neal, C. R., and Chenery, S. P.: A compilation of new and published major and trace element data for NIST SRM 610 and NIST SRM 612 glass reference materials, Geostandards Newslett., 21, 115–144, https://doi.org/10.1111/j.1751-908X.1997.tb00538.x, 1997.
Petrus, J. A. and Kamber, B. S.: VizualAge: A novel approach to laser ablation ICP-MS U-Pb geochronology data reduction, Geostand. Geoanal. Res., 36, 247–270, https://doi.org/10.1111/j.1751-908X.2012.00158.x, 2012.
Piccoli, P. M. and Candela, P. A.: Apatite in igneous systems, Rev. Mineral. Geochem., 48, 255–292, https://doi.org/10.2138/rmg.2002.48.6, 2002.
Schaaf, P., Weber, B., Weis, P., Groβ, A., Ortega-Gutiérrez, F., and Kohler, H.: The Chiapas Massif (Mexico) revised: New geologic and isotopic data and basement characteristics, Neues Jahrbuch fur Geologie und Paläontologie, Abhandlungen, 225, 1–23, 2002.
Solari, L. A., Ortega-Gutiérrez, F., Elías-Herrera, M., Schaaf, P., Norman, M., Ortega-Obregón, C., and Chiquín, M.: U-Pb zircon geochronology of Palaeozoic units in western and central Guatemala: Insights into the tectonic evolution of Middle America, Geol. Soc. Sp. Publ., 328, 295–313, https://doi.org/10.1144/SP328.12, 2009.
Solari, L. A., Ortega-Gutiérrez, F., Elías-Herrera, M., Ortega-Obregón, C., Macías-Romo, C., and Reyes-Salas, M.: Detrital provenance of the Grenvillian Oaxacan Complex, southern Mexico: a zircon perspective, Int. J. Earth Sci., 103, 1301–1315, https://doi.org/10.1007/s00531-013-0938-9, 2014.
Spikings, R., Cochrane, R., Villagomez, D., Van der Lelij, R., Vallejo, C., Winkler, W., and Beate, B.: The geological history of northwestern South America: From Pangaea to the early collision of the Caribbean large igneous province (290–75 Ma), Gondwana Res., 27, 95–139, https://doi.org/10.1016/j.gr.2014.06.004, 2015.
Svojtka, M. and Košler, J.: Fission-track dating of zircon by LA-ICP-MS, Goldschmidt Conference, Davos, Switzerland, Geochim. Cosmochim. Acta, 66, p. A756, 2002.
Tera, F. and Wasserburg, G. J.: U-Th-Pb systematics in three Apollo 14 basalts and the problem of initial Pb in lunar rocks, Earth Planet. Sc. Lett., 14, 281–304, https://doi.org/10.1016/0012-821X(72)90128-8, 1972.
Thomson, S. N., Gehrels, G. E., Ruiz, J., and Buchwaldt, R.: Routine low-damage apatite U-Pb dating using laser ablation–multicollector–ICPMS, Geochem. Geophy. Geosy., 13, Q0AA21, https://doi.org/10.1029/2011GC003928, 2012.
Torres, R., Ruiz, J., Patchett, P. J., Grajales, J. M., Bartolini, C., Wilson, J. L., and Lawton, T. F.: Permo-Triassic continental arc in eastern Mexico: Tectonic implications for reconstruction of southern North America, Geol. S. Am. S., 340, 191–196, https://doi.org/10.1130/0-8137-2340-X.191, 1999.
Torres de León, R.: Caracterización geológica y geocronológica de unidades metamórficas e intrusivas de la región centro–Oeste de la Cuenca del Rio Motagua, Sureste de Guatemala, Centroamerica: implicaciones en las conexiones Sur de México–Bloque Chortís, Universidad Nacional Autónoma de México, Posgrado en Ciencias de la Tierra, PhD Thesis, 221 pp., 2016.
Vermeesch, P.: Statistics for LA-ICP-MS based fission track dating, Chem. Geol., 456, 19–27, https://doi.org/10.1016/j.chemgeo.2017.03.002, 2017.
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, 2018.
Weber, B., Iriondo, A., Premo, W. R., Hecht, L., and Schaaf, P.: New insights into the history and origin of the southern Maya block, SE Mexico: U–Pb–SHRIMP zircon geochronology from metamorphic rocks of the Chiapas massif, Int. J. Earth Sci., 96, 253–269, https://doi.org/10.1007/s00531-006-0093-7, 2007.
Weber, B., Valencia, V. A., Schaaf, P., Pompa-Mera, V., and Ruiz, J.: Significance of provenance ages from the Chiapas Massif Complex (southeastern Mexico): redefining the Paleozoic basement of the Maya Block and its evolution in a peri-Gondwanan realm, J. Geol., 116, 619–639, https://doi.org/10.1086/591994, 2008.
Webster, J. D. and Piccoli, P. M.: Magmatic apatite: A powerful, yet deceptive, mineral, Elements, 11, 177–182, https://doi.org/10.2113/gselements.11.3.177, 2015.
Wolf, R. A., Farley, K. A., and Silver, L. T.: Helium diffusion and low-temperature thermochronometry of apatite, Geochim. Cosmochim. Acta, 60, 4231–4240, https://doi.org/10.1016/S0016-7037(96)00192-5, 1996.
Zeitler, P. K., Herczeg, A. L., McDougall, I., and Honda, M.: U-Th-He dating of apatite: A potential thermochronometer, Geochim. Cosmochim. Acta, 51, 2865–2868, https://doi.org/10.1016/0016-7037(87)90164-5, 1987.