New dating strategy untangles ancient burial heat from mountain cooling in old sand grains

Buried sandstones carry a hidden double record. Every zircon grain inside them remembers two things at once: the cooling of the mountains that produced it, long before it was eroded and swept into a basin, and the heat it endured after the sediment was buried. For decades, geologists have struggled to pull those two stories apart in so-called partially reset samples, where burial heating was strong enough to erase some of the original memory but not all of it. A new study by Birk P. Härtel, Eva Enkelmann, and Akeek Maitra of the University of Calgary, published in the journal Geochronology, presents a fresh algorithmic strategy that finally makes this disentangling possible without guessing what happened in the source region first.

The problem hinges on a temperature window that thermochronologists call the partial retention zone. In the zircon (U-Th)/He system, this is the range of temperatures, roughly between about 130 and 200 degrees Celsius depending on grain properties, within which radiogenic helium leaks out slowly rather than being fully retained or fully lost. A sedimentary sample that never entered this window after deposition is unreset: its single-grain helium ages all predate deposition and speak purely of source-region cooling. A sample heated well above the window is fully reset: its ages all postdate deposition and record cooling in the basin. But a partially reset sample, heated into the window, produces a spread of ages straddling the depositional age, a mixture of each grain’s individual pre-depositional history and the shared post-depositional one.

That mixture is precisely why partially reset samples have long been treated as a nuisance rather than a resource. In theory they are the most informative samples of all, because merely knowing that a sample spent time in the partial retention zone pins down the timing and temperature of maximum burial heating far more tightly than the loose upper or lower limits offered by unreset or fully reset samples. Correcting the measured ages for post-depositional helium loss would also recover when each grain’s source rock cooled, information that fully reset samples destroy entirely. The obstacle has always been that standard inverse thermal-history modeling, started at the depositional age, produces systematic biases: it underestimates burial temperature when many ages predate deposition, and overestimates it toward a full-reset scenario when most ages postdate deposition.

Earlier approaches tried to break the deadlock by combining two dating methods on the same grains. In the popular zircon U/Pb–(U-Th)/He double-dating scheme, the U/Pb age marks each grain’s crystallization, the high-temperature starting point of its thermal history, while the helium age records how much of the accumulated helium survived. Reiners and colleagues in 2005 devised a first-order correction based on the grain with the smallest total helium loss, and Fosdick and colleagues in 2015 introduced forward modeling of combined pre- and post-depositional paths evaluated with a Kolmogorov-Smirnov test. Both, however, rest on assumptions about the pre-depositional history, such as rapid post-crystallization cooling near the surface, that are often unknowable, especially for ancient sediments whose source rocks have long since eroded away.

The Calgary team’s innovation is elegant in its simplicity: stop assuming anything about the pre-depositional past, and instead exploit the few grains that record only the shared post-depositional story. These are syn-depositional grains, volcanic or shallow magmatic zircons that crystallized at or very near the time the sediment was laid down. Because they had almost no time to accumulate and lose helium before burial, their helium ages reflect essentially nothing but the basin’s thermal history. The researchers use inverse temperature-time modeling of these grains, performed in the widely used software HeFTy, to generate a large set of candidate post-depositional temperature-time paths.

Each candidate path is then put on trial. The team forward-models the post-depositional history to determine how much helium accumulated after deposition and what fraction of the pre-depositional helium survived. From the measured helium concentration they back-calculate a pre-depositional model age for every grain under that path, using a generic age equation that accounts for the decay of uranium and thorium through geologic time. The crucial test is a likelihood criterion: a plausible path must yield pre-depositional model ages that fall between each grain’s U/Pb crystallization age and the depositional age, because cooling through the partial retention zone in the source region can only have happened after the grain formed and before it was deposited. A Monte Carlo procedure incorporating the uncertainties on all three ages assigns a probability to every path, and a knee-detection algorithm selects the subset of highest-likelihood histories whose fits remain distinguishable from one another.

On synthetic data with known thermal histories, the method performs impressively. Where the classic correction of Reiners and colleagues pushed six of ten grains more than ten million years away from their true pre-depositional ages, the new strategy recovered maximum burial temperatures of 160 to 190 degrees Celsius at 25 to 20 million years ago, close to the true value of 168 degrees, and matched nearly all true pre-depositional ages within ten million years. Unlike the Fosdick approach, which evaluates the bulk sample and can overlook minor sources with distinct cooling histories, the new method delivers pre-depositional ages grain by grain, allowing the same kind of provenance interpretation applied to unreset samples, from clusters of grains that crystallized and cooled together to recycled grains from fully reset metamorphic sources.

To demonstrate the approach on real rocks, the team applied it to a Devonian sandstone from the Imperial Formation in the Mackenzie Plain of the Northern Canadian Cordillera, part of the vast Ellesmerian clastic wedge shed from terranes that collided with Laurentia in the Late Devonian to Mississippian. The zircon grains yielded U/Pb ages ranging from 364 million years to 2.9 billion years, while their helium ages spanned 57 million to 2.03 billion years, straddling the depositional age of about 364 million years, with roughly 73 percent of grains postdating deposition. That pattern confirmed a high degree of partial reset. Three young, near-syn-depositional grains with concordant U/Pb ages below 400 million years anchored the inverse modeling, and the knee-detection algorithm distilled six best-fitting thermal histories from 250 candidate paths.

The results were validated against an independent, higher-temperature chronometer measured on the very same grains: zircon Raman dating, which tracks radiation damage in the crystal lattice rather than helium. The pre-depositional Raman model ages preserved a bimodal distribution, separating Paleo- to Mesoproterozoic cooling of the Canadian Shield from Neoproterozoic to Paleozoic cooling in younger source terranes, and for one Raman band every single model age fell within the expected envelope. The modeled burial histories also agreed broadly with existing thermal models of the Imperial Formation built from apatite fission-track and helium data, pointing to Late Cretaceous to Paleogene heating at peak temperatures of roughly 100 to 160 degrees Celsius between 90 and 70 million years ago, and matching detrital muscovite argon-argon ages from the same formation.

The authors are candid about the method’s Achilles heel: it needs syn-depositional grains, ideally three to ten of them, which are plentiful in strata fed by active magmatic arcs but scarce in sediments derived from ancient cratons. They suggest screening strategies, laser-ablation double dating that can process dozens to hundreds of grains without bias, and fallbacks ranging from clast-based analysis to fully reset low-temperature thermochronometers. Near-fully reset samples also weaken the likelihood criterion, and simplifications such as uniform uranium-thorium zoning or an alpha-ejection-shaped initial helium profile can be refined in future work. Yet the payoff is a fundamental shift in perspective: rather than avoiding partially reset samples, geologists can now target them, extracting both the depth and timing of basin burial and the low-temperature cooling of vanished mountain belts from a single handful of sand. To spread the technique, the team released PaRACAS, an open-source Jupyter notebook that calculates, plots, and evaluates pre-depositional model ages from HeFTy output, and the same logic extends naturally to apatite U/Pb with helium or fission-track dating, promising a broader toolkit for reading Earth’s thermal memory.

Subject of Research: A new thermochronological modeling strategy for disentangling pre- and post-depositional thermal histories in partially reset detrital zircon samples

Article Title: Technical note: Disentangling pre- and post-depositional thermal histories in partially reset samples

Article References: Härtel, B. P., Enkelmann, E., & Maitra, A. (2026). Technical note: Disentangling pre- and post-depositional thermal histories in partially reset samples. Geochronology, 8(3), 475-494. https://doi.org/10.5194/gchron-8-475-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-475-2026

Keywords: thermochronology, zircon, U-Pb dating, helium dating, partial reset, sedimentary basins, thermal history modeling, detrital grains, Canadian Cordillera, geochronology, burial heating, provenance