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: Uncertainty in helium diffusion in zircon limits thermochronologic resolution: application to the Great Unconformity. (n.d.). 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 research centers on Milanković cycles, the periodic variations in Earth’s orbital eccentricity, axial tilt, and precession that redistribute the sunlight reaching the planet and thereby pace long-term climate change. In the modern Solar System, these cycles operate on well-known timescales, and they have been used to calibrate the geological time scale for the Cenozoic era with remarkable accuracy. Extending that approach into the Precambrian, however, has been hampered by a fundamental problem: numerical models of the Solar System’s orbital motion become chaotic and lose predictive power over tens of millions of years, and no astronomical solution can currently be trusted beyond roughly the last 100 million years. For intervals more than a billion years in the past, scientists have had to rely on theory alone to estimate how orbital frequencies differed from today’s values.

Sedimentary rocks offer a way around this limitation. When climate cycles driven by orbital variations imprint regular patterns on accumulating sediment—alternations between more resistant and more recessive beds, for example, or rhythmic changes in grain size and composition—the resulting cyclostratigraphy can be read as a recording of the astronomical forcing that produced it. The Hakatai Shale, deposited in shallow-water settings roughly 1.4 to 1.1 billion years ago during the Mesoproterozoic era, preserves such rhythms in striking detail. The team logged and analyzed stratigraphic sections at Red Canyon and Tapeats Creek within Grand Canyon National Park, conducting fieldwork under permit from the National Park Service, and measured the thickness and character of successive sedimentary cycles with centimeter-scale resolution.

The key to interpreting these rhythms lies in the physics of the Earth–Moon system. Tidal friction, the braking effect of lunar tides on Earth’s rotation, has steadily slowed the planet’s spin over geological time while pushing the Moon gradually farther away. As the day lengthens, the frequency of the climatic precession cycle—the wobble in Earth’s axis that changes how seasons align with the planet’s position around the Sun—changes in a predictable way. Because the precession signal modulates the amplitude of the eccentricity cycle, sedimentary records that capture both can be used to solve for the precession constant and, from it, the ancient Earth–Moon distance and length of day. This approach, known as TimeOpt and its Bayesian extension TimeOptBMCMC, was applied to the Hakatai Shale using the Astrochron software package, with 100,000 Monte Carlo samples used to constrain the statistical uncertainty of the reconstruction.

The analysis of the Tapeats Creek composite record, corrected for variations in sediment thickness, revealed a coherent suite of astronomical signals. The team identified cycles corresponding to climatic precession, orbital eccentricity, and obliquity, and used the ratios among them to test which cyclostratigraphic interpretation best fit the data. Among three competing interpretations of the dominant spectral peaks, the preferred option yielded sedimentation rates of a few centimeters per thousand years—values consistent with the quiet, low-energy depositional environments inferred independently from the rock’s lithology, which includes reworked microbial mat fabrics, wind-blown quartz grains, and pseudomorphs after evaporite minerals such as gypsum and anhydrite.

Beyond confirming that Milanković forcing operated in the Mesoproterozoic, the record delivered a surprise. The relative amplitudes of the astronomical forcing frequencies, particularly obliquity—the tilt of Earth’s spin axis—showed anomalous patterns compared with what present-day dynamics would predict. In the spectra of the Hakatai Shale, the strength of individual obliquity-related peaks shifted between different stratigraphic intervals in ways that mirror the behavior of state-of-the-art deep-time astronomical models, specifically the ZB23 solutions developed by Zeebe and colleagues, which extend orbital calculations back 3.5 billion years. In those models, the dominance of particular obliquity cycles changes through time as secular resonances among the planets drift in and out of critical configurations.

One such configuration involves the resonance angle associated with the motions of Mars and the inner planets, which can transiently disrupt the dominant obliquity cycle. Another involves a secular resonance that interferes with the main eccentricity cycle linked to the orbital frequencies of Earth and Jupiter. The Hakatai spectra show amplitude trends—weak expression of one eccentricity peak, enhanced power in a particular obliquity band—that are consistent with the models’ predictions for conditions around 1.2 billion years ago, including the possible influence of a resonance in which combinations of planetary orbital frequencies and Earth’s axial precession frequencies nearly coincide. The authors note that these patterns could also reflect a nonlinear climate response, in which interactions between multiple forcing frequencies generate combination tones that appear in the sedimentary record at sums and differences of the original periods.

Either interpretation carries weighty implications. If the amplitude anomalies record shifts in secular Solar System resonances, then the Grand Canyon strata provide the first empirical evidence from the rock record for how the gravitational architecture of the planetary system has changed over more than a billion years, complementing purely numerical approaches that are limited by chaos. If, instead, the signals arise from nonlinear climate dynamics, they illuminate how the Precambrian climate system responded to astronomical forcing in an atmosphere and ocean very different from today’s, before the rise of complex life and with substantially different greenhouse gas inventories. Distinguishing between these possibilities is a central goal of ongoing work, and the Bayesian inverse modeling framework applied here is designed to weigh such alternatives quantitatively.

The study builds on a growing effort to use geology as a probe of Solar System dynamics, sometimes described as mapping Solar System chaos with the geological record. Previous work by members of the team demonstrated that Milankovitch cycles preserved in 2.46-billion-year-old banded iron formations constrain the Earth–Moon system in the Paleoproterozoic, and theoretical studies have traced how tidal evolution reshaped the lunar orbit through resonant episodes. The Hakatai Shale now extends this empirical reach into the Mesoproterozoic with a record whose internal consistency—matching precession, eccentricity, and obliquity signals across two geographically separated sections—strengthens confidence that the rhythms are genuinely astronomical in origin rather than products of local tectonic or depositional noise.

The practical implications extend well beyond deep-time astronomy. Accurate knowledge of past astronomical frequencies underpins astrochronology, the dating method that uses orbital cycles to refine the geological time scale, and the new results constrain how those frequencies differed in the Precambrian, when shorter days and a closer Moon altered the pacing of climate cycles. The team’s cyclostratigraphic data and analyses have been made openly available through Zenodo, and the ZB23 astronomical solutions are publicly accessible, allowing other researchers to test and extend the reconstruction. As more ancient rhythmically deposited successions are examined with these tools, sedimentary rocks may continue to serve as long-term observatories of the heavens—recording, in ordinary mud, the slow gravitational conversation between Earth, the Moon, and the wandering planets.

Subject of Research: Reconstruction of Mesoproterozoic Earth–Moon dynamics and Solar System orbital evolution from Milanković cycles in Grand Canyon sedimentary strata

Article Title: Earth–Moon and Solar System history recorded in Mesoproterozoic Grand Canyon strata

Article References: Lantink, M. L., Eyster, A., Kocken, I. J., Meyers, S. R., & Zeebe, R. E. (2026). Earth–Moon and Solar System history recorded in Mesoproterozoic Grand Canyon strata. Nature Geoscience. https://doi.org/10.1038/s41561-026-02100-3

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02100-3

Keywords: Milanković cycles, Hakatai Shale, Grand Canyon, Earth–Moon system, Mesoproterozoic, cyclostratigraphy, solar system dynamics, obliquity, lunar distance, astrochronology, secular resonance, tidal evolution