Archive for Virtual Thermo2020

Session 1: Tuesday November 10, 2020, 10:00 AM (CST); 5:00 PM (CET)

Name Institution Title Video/Chat
Remi Rateau Trinity College, Dublin A geo-thermochronology study (LAFT-U/Pb and AHe) offshore West of Ireland – Age and exhumation of the North Porcupine HighVideo
Ian Hillenbrand University of Massachusetts, Amherst The rise and fall of the Acadian altiplano: Evidence for a Paleozoic orogenic plateau in New England, USAVideo
Audrey Margirier GFZ Potsdam New constraints on the timing of Carnegie Ridge subduction from low-temperature thermochronology (Ecuadorian Andes)Video

Session 2: Tuesday November 10, 2020, 5:00 PM (CST); Wednesday November 11, 10:00 AM (AEDT)

Name Institution Preliminary TitleVideo/Chat
Honcheng Guo Lehigh University Using continuous ramped heating to assess dispersion in (U-Th)/He ages from Transantarctic Mountain apatitesVideo
Murat Tamer University of Texas at Austin On step-etching experiments on fission tracks in apatite and effective track etch times Video
Rich Ketcham University of Texas at Austin TOPIC: Quantifying fission-track etching, and how it affects efficiency and reproducibility.

Session 3: Tuesday December 15, 2020, 10:00 AM (CST); 5:00 PM (CET)

Name Institution Title Video/Chat
Birk Härtel TU Bergakademie Freiberg Zircon Raman dating: Age calculation and closure temperature(s)Video
Chris Mark University College DublinIsotopic forward modelling: A new tool to interpret U-Pb thermochronometry dataVideo
Catherine Ross University of Texas at Austin Zircon (U-Th)/He impact crater thermochronometry and the effects of shock microstructures on helium diffusion kineticsVideo

Session 4: Tuesday December 15, 2020, 5:00 PM (CST); Wednesday December 16, 10:00 AM (AEDT)

Name Institution TitleVideo/Chat
Li YangChina University of Geosciences, Beijing Application of fission track thermochronology to porphyry belts in Eastern Kunlun Mountains, Qinghai-Tibet PlateauVideo
Samuel Boone University of Melbourne A global platform solution for big data in low-temperature thermochronology: AusGeochem and LithoSurferVideo
Ling Chung University of Melbourne The FastTracks digital fission-track analysis training moduleVideo
Extra discussion
Murat TamerUniversity of Texas at AustinProposal for a new inter-laboratory experimentVideo

Session 5: Tuesday February 23, 2021, 10:00 AM (CST); 5:00 PM (CET)

Name Institution Title Video/Chat
Daan Beelen Colorado School of MinesEvidence for a Late Cretaceous impact structure in Black Mesa, ArizonaVideo
Rachel Hoar Texas A&M UniversityHigh-resolution fission track analysis of detrital apatite grains grouped using UPb ages and REE concentrationsVideo
Andrew Donelick Gone Fission Mining LLC A mobile fission track microscope system for STEM education in remote communitiesVideo

November 10, 2020, Session 1

Talk 1

A geo-thermochronology study (LAFT-U/Pb and AHe) offshore West of Ireland – Age and exhumation of the North Porcupine High

Rémi Rateau1*, David Chew1

1 Department of Geology, Trinity College Dublin – Irish Centre for Research in Applied Geosciences


The offshore basins and basement highs west of Ireland are part of the northeast Atlantic Margin and are the result of the protracted break-up of Pangea and opening of the North Atlantic during the Mesozoic to Cenozoic. We acquired new geo-thermochronological data (dual-dating LAFT-U/Pb, AHe and zircon U/Pb) from borehole cuttings and seabed dredges and modelled both new and legacy data (QTQt) to better understand the thermal history and evolution of this part of the margin.

One area of interest is the Northern Porcupine High (NPH) which is a source of sediments for both the hyper-extended Porcupine Basin and Rockall Basin. The new results show that protracted exhumation of the high was associated with pulsed Mesozoic rifting (1 to 3 km of total erosion in total) and that the SE margin of the NPH started to be buried during the Early Cretaceous while the summit was still being exhumed and intruded by rift-related magmatism with an associated geothermal gradient increase. However, the new data fail to constrain a base Eocene exhumation event (visible on seismic) because of conflicting track length and AHe ages and the presence of coeval magmatism that may have affected some samples.

Zircon and apatite LA-ICP-MS U/Pb dating of basement samples on the SE margin of the NPH revealed Neoproterozoic (Lower Dalradian Supergroup) metasediments affected by a Late Caledonian thermal event. This new basement geochronology data, combined with previous published ages for the basement samples at the summit of the NPH, confirms that the low magnetic anomaly associated with the NPH represents Dalradian basement and that the sharp and straight magnetic boundary to the south probably represents the offshore extension of the Fair Head-Clew Bay lineament (the Irish extension of the Scottish Highland Boundary Fault, a major Caledonian suture).

Talk 2

The rise and fall of the Acadian altiplano: Evidence for a Paleozoic orogenic plateau in New England, USA

Ian Hillebrand1*, Michael Williams1, Cong Li1, Haiying Gao1

1 Department of Geosciences, University of Massachusetts at Amherst


High elevation orogenic plateaus are formed by a complex interplay of deep and surficial processes and influence a variety of Earth systems. However, limited exposures of plateau mid-crust are presently recognized, hindering understanding of these deeper processes. We present evidence for the existence of an orogenic plateau during and after the Devonian Acadian orogeny whose mid-crustal roots are exposed in the New England Appalachians. The four-dimensional crustal evolution of this paleo-plateau is constrained by the integration of new petrologic and geochronologic databases with geophysical imaging. Doubly thickened crust, widespread amphibolite to granulite-facies metamorphic conditions, a paleo-isobaric surface, and protracted mid-crustal anatexis indicate the development of a high elevation, low relief plateau by 380 Ma.40Ar/39Ar thermochronology shows a distinct thermochronologic signature with very slow cooling rates of 2-4°C/m.y. following peak metamorphic conditions. Thermochronologic data, trace element and Nd isotope geochemistry, and monazite petrochronology suggest a 50 m.y. lifespan of the plateau. Orogen parallel ductile flow and extrusion of gneiss domes resulted in plateau collapse, crustal thinning, and block-like exhumation at ca. 330-310 Ma. Thinning of the plateau crust may have led to the sharp 12-15 km step in Moho depth in western New England, possibly by reactivating the suture between Laurentia and accreted Gondwanan-derived terranes. 

Talk 3

New constraints on the timing of Carnegie Ridge subduction from low-temperature thermochronology (Ecuadorian Andes)

Audrey Margirier1,2*, Manfred Strecker1, Peter Reiners3, Ismael Casado1, Stuart Thompson3, Sarah W.M. George3, Alexandra Alvarado4

1 Institut für Erd-und Umweltwissenschaften, Universität Potsdam, Potsdam, Germany
2 Helmholtz Centre Potsdam, GFZ German Research Center for Geosciences, Potsdam, Germany
3 Department of Geosciences, University of Arizona, Tucson, United States of America
4 Instituto Geofísico, Escuela Politécnica Nacional, Quito, Ecuador


The Cenozoic growth of the Andes has been strongly influenced by the subduction dynamics and the superposed effects of climate. Recent studies proposed that the Carnegie ridge subduction controlled the late Cenozoic tectonic activity in Ecuador with the formation of a crustal sliver escaping northward; however, the timing of the arrival of the ridge and its effect on topographic growth remain unclear. We provided new thermochronological data from the Western Cordillera to evaluate the possible role of ridge subduction in prompting the growth of the Ecuadorian Andes and to pin point the timing of the ridge subduction. Time-temperature inverse modeling of this thermochronological dataset evidenced two cooling phases. A first one following intrusions emplacement in the Western Cordillera and a second one initiating at ~6 Ma. We interpret the first cooling phase to be associated with magmatic cooling. The second cooling phase, synchronous with the last cooling phase evidenced in the Eastern Cordillera, is likely to be associated with exhumation of the Western Cordillera. Based on our thermal modeling and geological cross-sections we propose that recent crustal shortening and rock uplift triggered exhumation of Ecuadorian Andes at ~6 Ma. We suggest that the onset of Carnegie Ridge subduction at ~6 Ma increased the coupling at the subduction interface and promoted shortening and rock uplift.

November 10, 2020, Session 2

Talk 1

Using continuous ramped heating to assess dispersion in (U-Th)/He ages from Transantarctic Mountain apatites

Hongcheng Guo1*, Peter K. Zeitler1, Bruce D. Idleman1, Annia K. Fayon2, Paul G. Fitzgerald3

1 Department of Earth and Environmental Sciences, Lehigh University, Bethlehem, PA 18015
2 Department of Earth Sciences, University of Minnesota, Minneapolis, MN 55455
3 Department of Earth Sciences, Syracuse University, Syracuse, NY 13210


We performed continuous ramped heating (CRH) on apatite suites from the Transantarctic Mountains (TAM) that are known to have significant intra-sample dispersion in He ages. Examining > 100 apatite grains from a total of six samples, we confirmed prior results showing that measured AHe ages have at least three-fold intra-sample dispersion with no obvious relationships between ages and effective uranium (eU) concentration or grain size.

Apatites showing simple volume diffusion behavior, characterized by unimodal incremental 4He gas-release curves, have younger He ages. In contrast, older grains generally show anomalous gas release in the form of sharp spikes and / or extended gas-release at high temperatures (i.e., >= 800 °C). Considerable age dispersion still remains in well-behaved apatites and exceeds what grain size, radiation damage, and analytical uncertainty can explain, but this dispersion appears to be related to variations in 4He diffusion kinetics. The screened AHe ages from well-behaved apatites together with kinetic information obtained from these grains suggest that the sampled region experienced slow exhumation prior to rapid exhumation beginning ca.  ~35 Ma, an interpretation consistent with studies suggesting that increasing incision rates at this time was related to initiation of glaciation and the Eocene-Oligocene climate transition. We attempted to correct the anomalous apatite ages by simply removing anomalous gas-release components (spikes and high-temperature release), but this approach yields some ages that are too young for the samples’ geologic setting, suggesting that the factors that lead to anomalous laboratory release behavior can impact both the expected radiogenic component as well as those that are anomalous.

From our observations we conclude that: (1) CRH screening can serve as a routine tool for AHe dating and offers opportunities to reveal first-order kinetic variations; (2) model-dependent age correction might still be possible but would require some means of estimating the broad proportions of 4He components incorporated into grains before and after cooling, and (3) interpretation of highly dispersed AHe ages will require assessment of individual-grain diffusion kinetics beyond that predicted by radiation-damage models. We further infer that many apatite grains will contain imperfections of varying kinds that contribute significantly to kinetic variability beyond that associated with radiation damage.

Talk 2

Step-etch experiments on fission tracks in apatite: How much has been done and where are we going?

Murat T. Tamer1*, Richard A. Ketcham1

1 Jackson School of Geosciences, University of Texas at Austin


In step-etch experiments, fission tracks are etched and confined track lengths measured in a series of steps. For apatite, step-etch experiments were used to create the standard single-step etch protocols around three decades ago based the concept of maximum etchable length, by monitoring the increase of mean track lengths to reach an optimum etch time. In these early step-etch studies, different sets of tracks were measured in each etch step, due to practical limitations. Recent developments and innovations in hardware and software for fission track analysis have enabled a new generation of step-etch studies, where individual tracks are monitored throughout an experiment, instead of random track selection in each etching step. These have provided new insights affecting the fundamentals of fission-track dating and modelling, such as the need to revise the linear track etch model, a difference in etching behavior between fossil and induced tracks, and the failure to reach the maximum etchable length over the prescribed etch time in the current standard single-step etch protocols.

We conducted a new series of step-etch experiments on unannealed induced, lab-annealed induced, and fossil tracks in Durango apatite to decipher along-track etch rates (VT(x)). Additional etch-anneal-etch (EAE) experiments aimed to determine the bulk etch rate (VB). We estimate the average VB as ~0.02 µm/s in Durango apatite, and 98% of EAE experiment tracks had etching velocities under 0.06 µm/s.  Generally, etch rates along unannealed tracks did not fall to VB until the 25-30s etching step, even when considering only tracks that were visible within the first 10 s of etching. Using 0.06 µm/s as a threshold for under-etching, ~78% of induced and ~53% of fossil tracks in our data were under etched in Durango apatite using the 5.5M 20s 21°C etching protocol, once again considering only tracks that were visible after 10 s. This procedure eliminates analyst bias in track selection, as it only considers a group of tracks that have etched for well over 10 seconds, resulting in an overall increase in mean length.  Analyst selection (bias) by normal measurement after 20 s lowers this number to ~60% for induced tracks. We expect that different analysts with different track selection criteria would get different values, and that the combination of under-etching and variable criteria underlies the poor level of community reproducibility of these measurements.  A quantitative understanding of along-track etching, based on data such as this study, can provide an avenue to addressing this important problem, and optimizing etching procedures.


TOPIC: Quantifying fission-track etching, and how it affects efficiency and reproducibility.

Led by: Rich Ketcham, Jackson School of Geosciences, University of Texas at Austin

Why do we etch apatite fission tracks the way we do?    Why do different lab groups etch with different recipes?  Does it even matter, if labs using the same etching get different results?  What are the advantages we seek, and what are the penalties we suffer?

A contribution, now in its public review period at, proposes that differences in how analysts select tracks to measure is actually more important for reproducibility than precisely how we etch.  However, it also proposes we don’t etch as much as we should, and that etching can be optimized to reveal more and better-defined confined tracks, while still retaining a quantitative link to experimental annealing data sets used for modeling.  Is this credible?

There will be a short presentation to kick things off, and then the floor will be open for discussion.  The more different perspectives we get, the better!

December 15, 2020, Session 3

Talk 1

Zircon Raman dating: Age calculation and closure temperature(s)

Birk Härtel1*, Raymond Jonckheere1, Bastian Wauschkuhn1, Lothar Ratschbacher1

1 Institute of Geology, TU Bergakademie Freiberg


Zircon Raman dating based on radiation damage accumulation and annealing is a debated concept but not yet an established thermochronological method. The α-disintegration of U, Th and their decay products creates lattice disorder in zircon (ZrSiO4) that can be quantified using the broadening of the Raman bands. Measuring the radiation damage and the U and Th content allows to calculate a zircon Raman date. Radiation damage is annealed at high enough temperatures.

We present a new empirical calibration of Raman bandwidth on radiation damage dose based on fission-track age standards. We derived an age equation to calculate zircon Raman dates of single grains. For multi-grain dates, standard procedures such as weighted mean, pooled, and isochron dates are applicable. We carried out isothermal annealing experiments, and fitted a Johnson-Mehl-Avrami-Kolmogorov and a distributed activation energy model to the annealed fraction of the bandwidth. From the kinetic parameters of the models, we estimated closure temperatures of zircon Raman thermochronometry at ~330-370 °C for internal ν2(SiO4) and ν3(SiO4) bands and ~250-300 °C for the external rotation band with a Raman shift of ~356 cm-1. This difference in closure temperatures offers the prospect of multi-temperature zircon Raman dating as a part of multi-method zircon dating.

Talk 2

Isotopic forward modelling: A new tool to interpret U-Pb thermochronometry data

Chris Mark1*, David Chew2, Nathan Cogné3, Andrew Smye4, Pieter Vermeesch5, and J. Stephen Daly1

1 School of Earth Sciences, University College Dublin, Dublin, Ireland;
2 Department of Geology, Trinity College Dublin, Dublin, Ireland;
3 Géosciences Rennes, Université Rennes 1, Rennes, France;
4 Department of Geosciences, Pennsylvania State University, University Park, USA;
5 Department of Earth Sciences, University College London, London, UK.


The availability of high-temperature thermochronometers suitable for generation of continuous thermal histories at mid- to lower-crustal temperatures (i.e., ≥ 400 °C) is limited. Available thermochronometers include the recently developed apatite and rutile U-Pb thermochronometers ( ≤ 550 and 640 °C; Kooijman et al., 2010; Cochrane et al., 2014) and arguably also the K-Ar system in white mica (sensitive to temperatures ≤ 500 °C.

Recent work has focussed on U-Pb analysis of apatite and rutile by sector-field and multi-collector LA-ICPMS to generate single-crystal U-Pb age profiles. Such profiles can be inverted to yield continuous thermal histories for high-temperature processes (e.g., Smye et al., 2018). Both rutile and especially apatite routinely incorporate non-trivial amounts of common-Pb during crystallisation (as opposed to radiogenic Pb generated by in-situ radionuclide decay), rendering them discordant in U-Pb isotope space. This common-Pb must be corrected for during age calculation. Such a correction necessitates critical assumptions regarding the isotopic composition of the common-Pb, and also the diffusive response of common- and radiogenic-Pb to heating. Additionally, both apatite and rutile can exhibit crystal growth and dissolution-reprecipitation reactions in the same temperature ranges at which measurable Pb diffusion occurs: neither behaves as a pure thermochronometer in all circumstances (e.g., Chambers and Kohn, 2012; Harlov et al., 2005).

Here, we exploit the observation that common-Pb is isotopically distinct from radiogenic Pb to predict the U-Pb isotopic composition of a given crystal arising from a proposed thermal history by forward modelling. We show that common-Pb can be exploited to validate the assumption of Pb-loss by volume diffusion; the thermal history predicted by age profile inversion; and age profiles arising from volume diffusion (as opposed to (re)crystallisation). We present a case study from the Eastern Alps.

Chambers, J.A., & Kohn, M.J., Am. Mineral., 97, 543–555 (2012); Cochrane, R., et al., Geochim. Cosmochim. Acta, 127, 39–56, (2014); Harlov, D.E., et al., Contrib. Mineral. Petrol, 150, 268–286 (2005); Kooijman, E., et al., Earth Planet. Sci. Lett, 293, 321–330, (2010); Smye, A.J., et al., Chem. Geol., 494, 1–18 (2018).

Talk 3

Zircon (U-Th)/He impact crater thermochronometry and the effects of shock microstructures on helium diffusion kinetics

Catherine H Ross1,2*, Daniel F. Stockli1, Desmond B. Patterson1, Sean P. S. Gulick1,2, Timmons Erickson3

1 Department of Geological Sciences, The University of Texas at Austin
2 Institute for Geophysics, The University of Texas at Austin
3 Johnson Space Center, NASA


Accurate age determinations of impact and cratering events remains difficult and often controversial and less than half of all known impact craters are regarded as accurately and precisely dated. Besides 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has become novel technique to date large- to medium-sized impact structures. ZHe ages can be fully reset in minutes at 1000°C, T reached commonly reached in central part of impact targets, whereas complete resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-4 years. However, there is a critical need to test the reliability of (U-Th)/He impact dating, in light of shock-induced microstructures and impact metamorphism, and to quantify helium diffusion kinetics in well-characterized variably shocked zircon. For this purpose, we investigated samples from two impact structures, the Chicxulub multi-ring basin and the Ries crater, to compare zircon diffusion kinetics from impact structures with different size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to determination of diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incrementally 10°C steps from 300°C to 600°C. We find that zircon with low-level shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features and granular textures classified by SEM are characterized by a dramatic decrease in helium retentivity, similar to radiation damage, due to the reduction in the effective domain size and the introduction of fast diffusion pathways.  This likely renders shocked grains more susceptible to impact-induced hydrothermal resetting. Hence, characterization of impact microstructure is critical for determining accurate impact ages, but also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

December 15, 2020, Session 4

Talk 1

Application of fission track thermochronology to porphyry belts in Eastern Kunlun Mountains, Qinghai-Tibet Plateau

Li Yang1,*, Wanming Yuan1, Shujiong Hong1, Zirui Feng1, Yunlei Feng1, Pengfei Tian1, Congmiao Cao1

1 Institute of Earth Sciences, China University of Geosciences, Beijing 100083, China


Studies on the preservation of ore deposits increasingly raise the concern from researchers in recent years. It is an important component of ore deposit geology. However, the extent to which ore bodies are preserved remains a poorly understood area of ore deposit research. The preservation of ore deposits is usually linked with the processes of uplift and exhumation. Thus, it is a reasonable and feasible method to evaluate the preservation of ore deposits based on the analysis and reconstruction of the exhumation and cooling history of mining area. In this dissertation, the hydrothermal ore deposits in the East Kunlun Mountains are chosen to be worked on. The exhumation and cooling history of selected mining area are reconstructed by the combination of fission track themochronology, apatite (U-Th)/He thermochronology, zircon U-Pb geochronology and fluid inclusions analysis. The relationship between buried depth of ore bodies and the amount of exhumation is clearly revealed. A real model is presented to outline the preservation and exhumation history affected by faulting activities. The conclusions of this dissertation will provide reasonable clues for further exploration.

The East Kunlun Mountains, located in the north of the Tibetan Plateau, is one of the significant metallogenic belts in China. In the region, we work on four hydrothermal ore deposits, including Harizha copper deposit, Halongxiuma molybdenum (copper) deposit, Qingshuihe copper and molybdenum deposit, Dahongshan iron deposit and Aikengdelesite molybdenum and copper deposit, to evaluate the preservation of the ore deposits. The zircon U-Pb dating data suggests that all of these five hydrothermal deposits formed from 253 Ma to 215 Ma in Triassic. The fluid inclusions from these deposits are characterized of high-medium forming temperature, high salinity and high-medium density. The mean ore forming depths ranging from 4.20 km to 7.42 km. The forming and post-forming processes of these four porphyry deposits are generally consistent.

All deposits experienced rapid cooling after magma emplacement, and the temperature dropped to the temperature of the surrounding rock, then suffered different modes of denudation cooling process. Four deposits underwent similar rapid-stable-rapid cooling pattern in this stage. But the Aikengdelesite deposit in this stage experienced a unique rapid-stable pattern that different from other four deposits. As evaluated in this dissertation, the preservation of these four deposits are different. The ore bodies of Harizha and Halongxiuma deposits are well preserved, while the preservation of Qingshiuhe deposit is not so well but still positive. However, the Aikengdelesite deposit is not well preserved.

Key words: thermochronology, preservation of ore deposits, East Kunlun, cooling and exhumation

Talk 2

A global platform solution for big data in low-temperature thermochronology: AusGeochem and LithoSurfer

Samuel C Boone1,2*, Fabian Kohlmann3, Moritz Theile3, Wayne Noble3, Barry Kohn1,2, Andrew Gleadow1,2, Stijn Glorie4, Martin Danišík5 and Renjie Zhou6

1University of Melbourne, School of Earth Sciences, Melbourne, Australia
2AuScope Geochemistry Network, Curtin University, John de Laeter Centre, Perth, Australia
3Lithodat Pty. Ltd., Melbourne, Australia
4University of Adelaide, Centre for Tectonics, Resources and Exploration, Department of Earth Sciences, School of Physical Sciences, Adelaide, Australia
5Curtin University, John de Laeter Centre, Perth, Australia
6University of Queensland, School of Earth and Environmental Sciences, Brisbane, Australia


The AuScope Geochemistry Network (AGN; have partnered with Lithodat Pty Ltd ( to develop the dual AusGeochem-LithoSurfer geochemistry data platforms. The open and cloud-based dual-platform solution will allow laboratories to upload, archive, disseminate and publicise their datasets, as well as perform statistical analyses and data synthesis within the context of large volumes of global geochemical data aggregated by the AGN, Lithodat and other platform users. While AusGeochem is designed specifically to serve the Australian geochemistry community (, the coupled LithoSurfer platform is designed as a global data repository for laboratories across the globe. Geosample information and geochemistry data can be bulk uploaded into the AusGeochem-LithoSurfer relational database through a specially developed ETL (Extract, Transform and Load), allowing legacy data to be uploaded via an easy drag-and-drop of csv files. The AusGeochem-LithoSurfer platforms also utilise an Open REST API which can be used to automatically upload data from analytical software and in-lab hard drives, enabling a flexible and automated data migration process that accommodates diverse data types while minimising manual data handling. In addition to facilitating data management, the Open REST API is designed such that any developer can build clients to interact with the platform to, for example, automatically retrieve data from its database, add-in enhanced data visualisation and create direct links to analytical equipment.

As part of this endeavour, representatives from four Australian low-temperature thermochronology laboratories (University of Melbourne, University of Adelaide, Curtin University and University of Queensland) are advising the AGN and Lithodat on the development of low-temperature thermochronology (LTT)-specific data models for the relational AusGeochem-LithoSurfer database. Adopting established international data reporting best practices, the LTT advisory group has designed database schemas for the fission track and (U-Th-Sm)/He techniques. A separate data model enables the archiving of digitised thermal history modelling results and metadata, making time-temperature reconstructions geospatially comparable and quantifiable. In addition to recording the parameters required for LTT analyses, the schemas include fields for reference material results and error reporting, allowing platform users to independently perform QA/QC on data archived in the database.

The advent of a LTT relational database heralds the beginning of a new era of structured Big Data in the field of low-temperature thermochronology. By methodically archiving detailed LTT (meta-)data in structured schemas, intractably large datasets comprising 1000s of analyses produced by numerous laboratories can be readily interrogated in new and powerful ways. These include rapid derivation of inter-data relationships, facilitating on-the-fly age computation, statistical analysis and data visualisation. With the detailed LTT data stored in relational schemas, measurements can then be re-calculated and re-modelled using user-defined constants and kinetic algorithms. This enables LTT analyses determined using different parameters to be equated and compared across regional- to global scales. Once standardised and normalised, these data can then be compared directly with other datasets, such as thermal conductivity, heat flow and other geochemistry data.

Talk 3

The FastTracks digital fission-track analysis training module

Ling Chung1*, Samuel C Boone1, Andrew Gleadow1, Malcolm McMillan1 and Barry Kohn1

1 School of Earth Sciences, University of Melbourne


Training for fission tack analysis is normally a very slow and labour-intensive process.  The development of the digital fission-track training module aims to deliver a guided training routine to the analyst’s computer, equipping researchers with equivalent level of confidence and skill to produce reliable and reproducible External Detector Method (EDM) and LA-ICP-MS/Digital Fission Track (LAFT) analyses using an image-based protocol. Originally, there are five stages of hand-on practice involved in comprehensive training conducted by Melbourne Thermochronology Research Group (MTRG). Chiefly utilising the MTRG-developed Fission Track Studio (FTS), a cross-platform dual software suite that is specialized for microscope control and image acquisition (TrackWorks) and image analysis (FastTracks), trainees are led to complete a series of carefully designed exercises to acquire the various skills in fission track dating. (1) Preparing grain mount for Durango (DUR) apatite provided. (2) Practicing grain selection and track-counting (track identification and region of interest, ROI, determination) by analysing seven sets of pre-imaged mica, DUR and Fish Canyon Tuff (FCT) apatites using FastTracks. (3) Challenging ‘unknowns’ by using the full capacity of FTS, i.e. capturing images of >25 selected grains each from the DUR mount and 20 well-characterized apatite reference mounts from the MTRG collection using TrackWorks and performing semi-automated analysis on grain and length images using FastTracks. (4) Acquiring single grain/spot uranium concentration using LA-ICP-MS. (5) Practicing age calculation. Reproducibility of data obtained from exercise (2) to (5), including single-grain track density, 3D confined track length measurements and average Dpar values, are constantly reviewed throughout the course.

Experience learned from this laboratory-based routine is a powerful tool for efficiently and comfortably examining data quality and evaluating source(s) of analytical discrepancy between analysts. The digital training module is a streamlined version of the laboratory-based practice. It allows trainees remote access to exercises for steps (2) to (5) above through analysing a series of selected grain and length image sets using FastTracks. Steps  (3) and (4) require specialized laboratory equipment but can be practised by providing pre-captured images of 6 ‘unknown’ samples with a mixture of many suitable and unsuitable grains from which an appropriate selection of grains for analysis can be made.  Similarly, provision of LA-ICP-MS data for the grains enables calculation of fission track ages from the analysed grains. In addition, a sub-module allows EDM users calculating their own user-specific zeta-calibration through analysis of images of apatite-mica sample pairs and co-irradiated external detectors from standard glasses. Together with proposed analytical solutions and a list of recommended reading list, trainees are able to evaluate their progress by comparing their data with solution files on a grain-by-grain and track-by-track basis. Images of two samples used in previous inter-laboratory studies (Ketcham et al., 2015 & Ketcham et al., 2018) are part of the exercise for data quality evaluation. As the training package is cloud-stored, researchers, who do not have access to suitable facility to conduct full suite of analysis or cannot physically travel for training due to various restrictions, can become trained without face-to-face tutelage or specialised equipment. In collaboration with two international laboratories, the module is being tested on both experienced conventional fission track analysts and untrained students and augmented for improved content and usability.

Development of the digital training module enables a new coordinated digital fission track analysis stream, whereby researchers can outsource sample preparation and image capture to laboratories equipped with suitable equipment. Captured image stacks and parent isotope concentrations, in the case of the LA-ICP-MS technique, would then be returned electronically to the newly trained researcher for digital fission track analysis and interpretation. This advance will enhance the accessibility and affordability of this powerful technique and make digital fission track analysis achievable for geoscientists globally.

Ketcham, R. A., Carter, A., & Hurford, A. J. (2015). Inter‐laboratory comparison of fission track confined length and etch figure measurements in apatite. American Mineralogist, 100 (7), 1452– 1468.‐2015‐5167
Ketcham, R. A., van der Beek, P., Barbarand, J., Bernet, M., & Gautheron, C. (2018). Reproducibility of thermal history reconstruction from apatite fission‐track and (U‐Th)/He data. Geochemistry, Geophysics, Geosystems, 19, 2411– 2436.

February 16, 2021, Session 5

Talk 1

Evidence for a Late Cretaceous impact structure in Black Mesa, Arizona

Daan Beelen1*, Ray Donelick2

1 SAnD Consortium, Colorado School of Mines, 1500 Illinois St, Golden, Colorado, USA
2 Partners LLC, 1075 Matson Road, Viola, Idaho, USA


This study presents evidence for a hitherto undiscovered impact structure in Black Mesa, NE Arizona. The structure’s macrostructure is mapped on satellite data as a series of concentric catchments and landscape geometries approximately 7 kilometers in diameter, with upward tilting strata dipping from the center of the structure. Subsequent field research has yielded observations of meso-scale impact structures in the form of melt breccia’s and titled blocks with abundant pseudotachylite melt. Partially molten and shocked zircon grains inside the molten breccias were identified using a SEM, and X-ray diffraction indicates elevated amounts of ultra-refractory elements within the same samples. The structure lies within strata belonging to the Upper Cretaceous Toreva Member and Wepo Members, Mesa Verde Group, Upper Cretaceous.

Apatite fission track (AFT) data were measured for two samples showing evidence of partial melting. Apatite grains – and their associated fission track data – were sorted into sub-groups using rare earth element concentrations and uranium-lead ages. The pooled AFT ages for all sub-groups studied are shown in the figure below.


Seven apatite grain sub-groups are common to both samples studied. For sample P10014_004, apatite grain sub-groups 2A and 4A yield pooled AFT ages of 70.6 -12.4/+15.1 Ma (95% CI; 14 grains) and 63.2 -8.6/+9.9 (95% CI; 30 grains), respectively, and a weighted mean age of 65.5 ± 7.7 Ma (2σ). The evidence presented here insinuates to the presence of the largest impact structure in the state of Arizona and may have implications to the Cretaceous–Paleogene mass extinction, although more work is necessary to prove these theories.

Field work on the Navajo Nation was conducted under a permit from the Navajo Nation Minerals Department. Any person(s) wishing to conduct geologic investigations on the Navajo Nation must first apply for and receive a permit from the Navajo Nation Minerals Department, P.O. Box 1910, Window Rock, Arizona 86515 and Telephone No. +1 (928) 871-6587

Talk 2

High-resolution fission track analysis of detrital apatite grains grouped using UPb ages and REE concentrations

Ray Donelick1, Cleber Soares2, Liu Zhaoqian3, Rachel Hoar4,*, Daan Beelen5, Andrew Donelick6, and Richard Carlton7

1 Partners LLC, 1075 Matson Road, Viola, Idaho, USA
2 Chronuscamp Research, Rua João Rodrigues dos Santos 60, Itapira, São Paulo, Brazil
3 Key Laboratory of Tectonics and Petroleum Resources, China University of Geosciences, Ministry of Education, Wuhan 430074, China
4 Geology & Geophysics, Texas A&M University, College Station, Texas, USA
5 SAnD Consortium, Colorado School of Mines, 1500 Illinois St, Golden, Colorado, USA

6 Gone Fission Mining LLC, 238 E 5th Street, Walsenburg, Colorado, USA
7 Navajo Nation Minerals Department, PO Box 1910, Window Rock, Arizona, USA


O’Sullivan et al. (2018, G3, v. 19, no. 9; 2020, Earth-Science Reviews, v. 201, no. 103044) demonstrate how UPb age and rare earth element (REE) concentrations can provide useful genetic and provenance information for detrital apatite grains. Here, we use UPb age, REE concentrations, and parameter Rmr0 to bin apatite grains – and their associated fission track data – into kinetic populations for thermal history modeling. Each UPb age is based on a single LA-QICP-MS spot analysis with an assumed common-Pb component (after Chew and Donelick, 2012, MAC Short Course, v. 42). The REEs are C1-chondrite normalized using – as a matrix-matched standard – the same Durango crystal studied by Chew et al. (2016, Chemical Geology, v. 435). Kinetic parameter Rmr0 is determined using LA-QICP-MS data (Na, Mg, Al, Si, P, S, Cl, Ca, Mn, Fe, As, Br Y, 14 REEs, Hg, Pb, Th, U), assuming only chlorine (measured after Chew et al. 2014, Geostandards and Geoanalytical Research, v. 38; US Patent 8,901,485) and fluorine (calculated from chlorine) in the halogen site.

Conventional AFT analysis uses measured values such as Dpar (US Patent 5,267,274) or elemental compositions (e.g., chlorine concentration from EPMA) to distribute apatite fission track data into kinetic populations. The protocol used here offers a higher degree of resolution of the annealing kinetics, because the apatite fission track data are first distributed into provenance groups defined by UPb age and REE concentrations. Within these provenance groups, apatite fission track data can be further distributed into kinetic populations using Rmr0, much like the conventional analysis.

This work was funded in part by The National Science and Technology Major Project of China (No.2016ZX05035-001-003 and No.2016ZX05002-006-007) and National Natural Science Foundation of China (No.41302112).

Talk 3

A mobile fission track microscope system for STEM education in remote communities

Ray Donelick1, Andrew Donelick2,*, Cleber Soares3, Tais Fontes Pinto4, Zachary Dodds5, Daan Beelen6, Murat Tamer7

1 Partners LLC, 1075 Matson Road, Viola, Idaho, USA
2 Gone Fission Mining LLC, 238 E 5th Street, Walsenburg, Colorado, USA
3 Chronuscamp Research, Rua João Rodrigues dos Santos 60, Itapira, São Paulo, Brazil
4 Universidade Federal do Rio Grande do Sul, Porto Alegre, Rio Grande do Sul, Brazil
5 Computer Science Department, Harvey Mudd College, Claremont, California, USA
6 SAnD Consortium, Colorado School of Mines, 1500 Illinois St, Golden, Colorado, USA
7 Jackson School of Geosciences, University of Texas at Austin, Austin, Texas, USA


We tested and demonstrated components of a mobile fission track microscope system, at American Associate of Petroleum Geologists, Annual Conference and Exhibition, in San Antonio, Texas during May 2019.  The microscope system and the technology behind it can be hauled anywhere – such as from Idaho to our AAPG exhibitor’s booth in Texas – using a personal vehicle and requires simple climate-controlled space with electricity.

Fission track thermochronology – its calibration and its application – will benefit from datasets comprised of 1000s of datum points per sample. We envision employing young students (in USA, Grades 6-12 for example) to collect the bulk of these data for select interesting samples using our mobile fission track microscope system. We believe this system can be used to excite and educate young students in a range of science, technology, engineering, and mathematics (STEM) disciplines. The young students can do the heavy lifting of finding and initially characterizing mineral grains for fission track analysis using various AI tools. Results from the young students can be vetted by digitally connected and collaborating professional scientists. This mobile fission track laboratory enables distance learning with a personal touch while applying state-of-the-art science in collaboration with an international network of passionate and skilled scientists.

We envision the following value propositions:

  • For young students
    • introduction to STEM using state-of-the-art science over a range of subjects: natural sciences from isotope geology to tectonics to meteorites, chemistry and physics, computer science
    • hands-on experience in a true ‘peer laboratory’ on par with any advanced laboratory anywhere