The tool in question is zircon (U-Th)/He thermochronology, often abbreviated ZHe. The method exploits a simple but powerful principle. Zircon crystals incorporate uranium and thorium when they grow, and these radioactive elements decay over time, producing helium atoms as daughter products. At high temperatures, helium diffuses out of the crystal almost as fast as it is produced, so the clock stays open. As a rock cools below a characteristic temperature range, diffusion slows dramatically and helium begins to accumulate. The measured age therefore records how long the crystal has resided at relatively low temperatures, typically below about 200 to 250 degrees Celsius for zircon. Because erosion brings rocks closer to the cool surface, these ages can be read as a thermal archive of mountain building, burial, and denudation over billions of years.

What makes the zircon helium system both remarkably powerful and notoriously complicated is radiation damage. Every decay event that produces a helium atom also sends a heavy alpha particle recoiling through the crystal lattice, leaving a trail of atomic-scale disruption. Over hundreds of millions of years, this damage accumulates at a rate governed by the crystal’s uranium and thorium content. Crucially, damage is not permanent: at elevated temperatures it can anneal, healing the lattice in a process described by fission-track annealing kinetics. The result is that two zircon grains from the same hand sample, sharing an identical thermal history, can behave entirely differently as thermochronometers simply because they contain different amounts of uranium and thorium and therefore different degrees of radiation damage. At low damage levels, increasing damage raises the closure temperature; at high damage levels, the trend reverses and the closure temperature falls.

To capture this complex behaviour, researchers rely on the Zircon Radiation Damage and Annealing Model, or ZRDAAM, developed by William Guenthner and colleagues in 2013. The model treats a crystal as a mixture of undamaged and damaged domains whose combined diffusion kinetics shift as damage accumulates and anneals through time. Both camps in the Great Unconformity debate have used ZRDAAM to interpret their zircon helium data. One influential analysis argued for a broadly synchronous, glacially driven origin tied to Neoproterozoic snowball Earth events, while another concluded that the unconformity developed diachronously through regional tectonic processes long before the glaciations. In the Eastern Grand Canyon, one team even reconstructed a thermal history spanning 1700 million years, reporting temperature constraints to within less than 10 degrees Celsius for intervals between 700 and 250 million years ago. The new study, led by Matthew Fox of University College London together with Adam Smith, Pieter Vermeesch, Kerry Gallagher, and Andrew Carter, asks a deceptively simple question: how well do we actually know the parameters that make such precise reconstructions possible?

The answer hinges on the Arrhenius relationship, the fundamental equation describing how diffusivity depends on temperature. Two parameters define it: the activation energy, which sets the slope of the relationship, and the frequency factor, which sets its intercept. These are extracted from step-heating experiments in which zircon crystals are progressively degassed under vacuum, and the two parameters are strongly correlated with one another. The original ZRDAAM calibration combined these parameters into single diffusivity or closure-temperature values, discarding information about their correlation. Worse, the kinetics of the theoretical minimally damaged crystal were obtained by extrapolating measured frequency factors down two orders of magnitude using a power-law relationship, a procedure that introduces uncertainties with no obvious way to quantify them. The kinetics of the extremely damaged end member rested on a single crystal, sample N17.

Fox and colleagues re-analysed the very same diffusion dataset using a fundamentally different statistical approach. Rather than fitting straight lines through individual Arrhenius plots and extrapolating, they modelled the measured diffusivities directly through the radiation damage framework, fitting four parameters representing the two hypothetical end-member crystals. They then sampled the full range of parameter values consistent with the data using a Bayesian Markov Chain Monte Carlo algorithm, tuning the proposals so that roughly 20 percent of candidate models were accepted and running the chain until one million parameter sets had been collected. Because the degassing experiments contained different numbers of heating steps, the team weighted the misfit of each experiment so that longer experiments would not dominate the calibration. The result is not just a single best-fit model but a full posterior probability distribution, complete with a covariance matrix that preserves the correlations between parameters.

The correlations turned out to be extreme. Within each end-member crystal, the activation energy and frequency factor are almost perfectly correlated, with correlation coefficients of 0.999 for the damaged end member and 0.991 for the pristine one. In contrast, the parameters of the two end members are essentially independent of each other. This structure matters enormously, because the strong correlations mean that many combinations of parameters produce nearly identical diffusion behaviour, while the residual scatter around the best fit reveals genuine irreducible uncertainty. Neither the original calibration nor the new one fits every data point perfectly, and the new maximum-likelihood parameters for the minimally damaged crystal differ noticeably from the original values, even though the original values still fall along the same correlation trend.

When the team propagated these uncertainties forward through a representative billion-year thermal history, mimicking conditions recorded by ancient rocks in Minnesota, the consequences were striking. Twenty simulated zircon crystals spanning effective uranium concentrations from 31 to 2828 parts per million were aged a thousand times over, each time with different diffusion parameters drawn from the posterior distribution. For a specific amount of radiation damage, the predicted ages dispersed by hundreds of millions of years. At an effective uranium concentration of about 1600 parts per million, predicted ages ranged from roughly 50 to 550 million years for the same thermal history. Adding realistic variation in grain size, parent-isotope zonation, inclusions, or broken grains would widen that spread further. In other words, the age-eU patterns that researchers routinely interpret as fingerprints of specific cooling events carry an intrinsic scatter that has largely gone unaccounted for.

The team also tested the effect inside a full thermal-history inversion by modifying the widely used QTQt software to resample the four diffusion parameters from their joint posterior distribution during the search. Applying the method to the Grand Canyon dataset of McDannell and colleagues, they found that the broad shape of the inferred thermal history survived, but the credible intervals widened, and the posterior distribution split into two distinct families of solutions, one including and one excluding a second cooling episode around 400 to 350 million years ago. The practical implication is sobering: thermal histories reconstructed with zircon helium data may be less well resolved than previously appreciated, and the tight temperature constraints sometimes reported for deep time may partly reflect overconfidence in the adopted kinetic parameters rather than genuine resolving power in the data.

What can be done? The authors outline a path forward. More diffusion experiments are needed to characterise helium kinetics across a wider range of radiation damage values, including replicate measurements to quantify dispersion. Natural laboratories, regions whose thermal histories are independently well constrained, could be used to calibrate diffusion parameters against observed ages and complementary thermochronometers. In the meantime, forward modelling offers a pragmatic workaround: predicting the expected age spread for a plausible thermal history and folding that spread into the uncertainties used in inversions. Screening crystals with ramped heating experiments or zircon helium-4/helium-3 thermochronology could identify anomalous diffusion behaviour before it contaminates a dataset. And because capturing the true age spread at a given uranium concentration may require 40 or more grains per sample, far more than the 5 to 30 typically analysed, the field may need to rethink its sampling strategies. The Great Unconformity’s missing billion years remain as enigmatic as ever, but this study makes clear that resolving them will demand a more honest accounting of the uncertainties hidden inside the crystals themselves.

Subject of Research: Uncertainty in helium diffusion kinetics in zircon and its impact on thermochronologic reconstruction of the Great Unconformity

Article Title: Uncertainty in helium diffusion in zircon limits thermochronologic resolution: application to the Great Unconformity

Article References: Fox, M., Smith, A. G. G., Vermeesch, P., Gallagher, K., & Carter, A. (2026). Uncertainty in helium diffusion in zircon limits thermochronologic resolution: application to the Great Unconformity. Geochronology, 8(3), 463-474. https://doi.org/10.5194/gchron-8-463-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-463-2026

Keywords: thermochronology, zircon, helium diffusion, Great Unconformity, radiation damage, ZRDAAM, Geochronology, thermal history modeling, Grand Canyon, Precambrian, Markov Chain Monte Carlo, closure temperature

The setting alone makes the work compelling. Bangladesh sits atop the Bengal Basin, one of the world’s great sedimentary repositories, where strata eroded from the rising Himalaya and carried by the Ganges, Brahmaputra and Meghna rivers have buried the ancient crust under hundreds to thousands of metres of Late Permian to Quaternary sediment. In the northwest, however, the Precambrian basement rises unusually close to the surface. The shallowest known occurrence lies at just 128 metres depth at Maddhapara, where diorite-tonalite-granodiorite rocks have been dated at 1.64 to 1.73 billion years using U-Pb SHRIMP analysis of zircon. That age places the Dinajpur Block firmly within the Paleoproterozoic, an era when the planet’s landmasses were assembling into the supercontinent Columbia, also known as Nuna.

Until now, drill hole GDH-62 had never been examined petrographically. The team collected nine core samples from the Geological Survey of Bangladesh repository and prepared twenty-seven thin sections, two or three per sample to capture mineralogical variation along and across the core axis. Under a standard polarizing microscope, they identified minerals and counted them along twenty-five to thirty-five grid-based traverses per section, estimating modal compositions visually with an uncertainty of about plus or minus five percent for major minerals. The approach is deliberately classical, but it is exactly the kind of careful, quantitative petrography that turns a stack of grey core boxes into a readable geological history.

The results split cleanly into two rock families. The felsic rocks, which occupy about 20.2 metres of the core or roughly 73 percent of the studied interval, are granodiorite: massive, grey to greenish grey, leucocratic to mesocratic, and fine to medium grained. Their modal composition ranges from 33.4 to 49.2 percent plagioclase, 11.8 to 27.0 percent quartz, 6.6 to 27.2 percent biotite and 6.7 to 11.7 percent K-feldspar, with secondary chlorite and epidote and accessory muscovite, zircon, hematite and opaque minerals. Plotted on the standard Quartz-Alkali Feldspar-Plagioclase discrimination diagram, every felsic sample falls within or adjacent to the granodiorite field, some near the granodiorite-tonalite boundary. The textures are hypidiomorphic and inequigranular, with fine- to medium-grained interlocking fabrics that speak of slow crystallization in a deep plutonic environment.

The microscopic detail is where the rocks begin to tell their deformation story. Quartz grains show wavy, undulose extinction, a classic signature of strain accumulated in the crystal lattice, and locally the quartz has recrystallized into finer equant aggregates. Biotite occurs as blade-like, euhedral to subhedral crystals with light brown to yellowish brown pleochroism, and many grains are twisted, bent and partially foliated. Aligned biotite and chlorite along recrystallized quartz grain boundaries define a weak foliation whose intensity, notably, decreases with depth in the section, hinting at a gradient in deformation intensity that the authors caution cannot be quantified from petrology alone. K-feldspar displays textbook Carlsbad twinning and hosts inclusions of acicular muscovite, biotite and rare epidote, the latter possibly pointing to hydrothermal alteration. Plagioclase shows weak to moderate sericitization, and much of the biotite is altered to chlorite, evidence that fluids, possibly of hydrothermal origin, percolated through these rocks long after they crystallized.

Cutting through the granodiorite near the base of the hole is the study’s most eye-catching discovery: an appinite dyke, massive, greenish grey to very dark green, and dominated by hornblende at 67 to 73 percent of its mode, accompanied by biotite, plagioclase, chlorite, hematite and accessory opaque minerals, with calcite forming micro-veins and veinlets. Appinites are an unusual and petrogenetically significant class of rock, typically hornblende-rich mafic intrusions associated with hydrous, mantle-derived magmas in convergent settings. In thin section the appinite shows hypidiomorphic, inequigranular and locally poikilitic textures, with hornblende enclosing plagioclase crystals. Hornblende exhibits the characteristic amphibole cleavage angles of 56 and 124 degrees with extinction angles between 12 and 34 degrees, while reddish brown biotite suggests a titanium-rich composition that the authors flag as potentially significant for petrogenesis. A tentative identification of orthopyroxene, based on pleochroism and parallel extinction, remains unconfirmed without optical or chemical constraint on its 2V angle.

The contact between the dyke and its host is sharp and non-gradational, with no chilled margin, and the authors interpret it as evidence that the appinite was emplaced after the felsic host had consolidated, most likely as a hypabyssal dyke within the main acid to intermediate association. From granodiorite to appinite, quartz and plagioclase decrease while biotite increases, tracking the shift from felsic to mafic chemistry. What petrography cannot constrain, the researchers emphasize, is the emplacement depth, a limitation they state plainly rather than paper over. That honesty matters in a field where over-interpretation of thin-section evidence is a chronic hazard.

Perhaps the most technically rich part of the study is its analysis of vein microstructures, which function as a kind of strain gauge for events long after the magmas cooled. In the granodiorite, lensoidal quartz veins with smooth vein-wall interfaces and lensoid wall-rock fragments point to a fracture-controlled vein system formed by syntaxial growth, in which crystals grow from both walls toward the vein centre. The small angle between the growth phase and the crack wall suggests a low growth rate, and the absence of fluid inclusions indicates that the vein crystals were deformed under non-hydrostatic stress during or after formation. In the appinite, by contrast, monomineralic calcite veins show slightly curved growth and selvages at their margins, hallmarks of antitaxial growth in which crystals expand outward from the vein centre, recording progressive, stress-related deformation over time. Minor calcite veins oriented obliquely and sub-parallel to the major vein hint at changing stress directions and possible fracture reactivation across the life of the intrusive complex.

The comparison the authors draw with the West African Craton gives the work its global resonance. Paleoproterozoic granitoids of the Bas Draa inlier in Morocco’s Anti-Atlas, dated around 2.1 billion years, share calc-alkaline, subduction-related characteristics with the Dinajpur Block, and both regions feature sodic TTG-style suites derived from partial melting of hydrated mafic sources. Mafic and appinitic dykes in GDH-62 resemble the lamprophyric dyke swarms of the West African Craton, though the Bangladeshi examples point more strongly to metasomatic and mantle-wedge influences than to plume-related processes. The similarities support the interpretation of GDH-62 as a fragment of Paleoproterozoic arc crust, a piece of a convergent-to-collisional margin that helped stitch together Columbia. For a country whose geology is usually described in terms of deltas, floods and gas fields, the message is striking: beneath the youngest sediments lies a two-billion-year-old arc, complete with dykes, veins and deformation histories, waiting in core boxes for anyone patient enough to grind the thin sections.

Subject of Research: Petrographic characterization of Precambrian crystalline basement rocks in drill hole GDH-62, Dighipara, northwest Bangladesh

Article Title: Petrographic characterization of the crystalline basement in Dighipara, Northwest Bangladesh

Article References: Zahan, H., Ameen, S. M. M., Hossain, M. S., & Hossain, M. S. (2026). Petrographic characterization of the crystalline basement in Dighipara, Northwest Bangladesh. Discover Geoscience, 4(1), Article 332. https://doi.org/10.1007/s44288-026-00703-2

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00703-2

Keywords: petrography, crystalline basement, granodiorite, appinite, Bangladesh, Dinajpur Block, Bengal Basin, Precambrian, vein microstructures, Columbia supercontinent, drill core, hornblende