Simple 1D Model Reveals When Rock Cooling Ages Mislead Geologists

Every mountain range tells a story, and geologists have long relied on a set of clever dating techniques to read it. Low-temperature thermochronology measures how radioactive decay products are retained in minerals such as apatite and zircon, allowing scientists to reconstruct when rocks cooled as they rose through the upper crust. But a new modeling study warns that the cooling histories these methods produce may not mean what researchers often assume. Geological processes such as erosion, sediment burial, faulting, and even the wholesale removal of the mantle root beneath a continent can perturb the thermal state of the crust in ways that decouple a rock’s temperature history from its actual path to the surface.

The study, published in the journal Geochronology by Dawn A. Kellett of the Geological Survey of Canada-Atlantic and David M. Whipp of the University of Helsinki, explores these effects using a deliberately simple one-dimensional modeling code called Tc1D. Rather than attempting to simulate every detail of complex geodynamic events, the software solves the transient heat transfer equation through a full 125-kilometer-thick lithospheric column, tracking the temperatures experienced by individual rock particles as they travel toward the surface. Those recorded thermal histories are then fed into established annealing and diffusion models to predict ages for four widely used thermochronometers: apatite (U–Th)/He, apatite fission-track, zircon (U–Th)/He, and zircon fission-track.

The authors ran a large suite of scenarios spanning roughly 50 to 55 million years of model time, covering erosional exhumation at constant, stepped, and exponentially decaying rates; sedimentary burial followed by erosion; thrust and extensional faulting; and instantaneous delamination of the lithospheric mantle. Each scenario was run in two variants. In the fixed Moho case, erosional exhumation is balanced by crustal flux at depth, mimicking a syn-collisional setting where the crust maintains a steady thickness. In the moving Moho case, the entire lithosphere moves toward the surface as material is removed, thinning the crust in a post-tectonic fashion. The difference between these two assumptions, the researchers found, is far from trivial.

For identical depth histories, moving Moho models almost always produced younger cooling ages than their fixed Moho counterparts. The reason lies in the thermal structure: a thicker starting crust contains more radiogenic heat-producing elements, raising the geothermal gradient and the initial temperature of the tracked particle. In one comparison, a subtle difference in the upper-crustal gradient of just a few degrees per kilometer translated into zircon helium and fission-track ages that were 15 to 19 percent younger. Across all models, the average age difference between the two Moho variants was about 10 percent, but many cases exceeded 20 percent, with differences as large as 100 percent for zircon fission-track ages in scenarios involving large amounts of exhumation. Because typical analytical uncertainties for these methods run 10 to 20 percent, the authors argue these differences are genuinely detectable in measured data, yet Moho evolution is rarely considered in thermal history studies.

The models also quantify how heat moves through the crust using two dimensionless numbers. The Péclet number compares heat transfer by advection, the physical transport of warm rock toward the surface, against conductive heat transfer through stationary rock. Values well above one indicate that erosion is efficiently advecting heat, bending the geotherm and producing curved cooling paths that deviate from the rock’s depth history. The dimensionless heat production parameter plays a similar role for radiogenic heating, which can dominate in crust rich in uranium, thorium, and potassium. Together, these parameters explain why the same exhumation history can yield strikingly different thermal histories depending on crustal heat production and starting crustal thickness.

The sedimentary burial experiments revealed another temporal subtlety. When 9 kilometers of sediment were deposited over a basement rock, peak temperatures lagged peak burial depths by several million years, because conduction takes time to warm the buried rocks. With extremely rapid sedimentation of 9 kilometers in just 1 million years, the lag stretched to about 8 million years, meaning a rock’s thermal history can record a heating event that appears asynchronous with the burial that caused it. Radiogenic heat production in the sediments and underlying crust also proved critical: doubling the heat production raised peak temperatures by tens of degrees and shifted predicted zircon helium ages by roughly 10 million years.

Faulting scenarios showed that only some thermochronometer ages actually date fault slip. In thrust models where the hanging wall was rapidly uplifted during faulting, apatite fission-track, zircon helium, and zircon fission-track ages clustered around the faulting interval, nicely documenting the event. But in scenarios where the tracked particles sat at deeper structural levels, most ages were instead controlled by post-faulting erosional exhumation, and footwall rocks heated by juxtaposition against the hot hanging wall delayed significant cooling by around 10 million years after fault motion ceased. Frictional heating along the fault surface was ignored, so the models likely underestimate the thermal perturbation of real faults.

Perhaps the most striking result concerns lithospheric delamination, the density-driven foundering of the continental mantle root that juxtaposes hot asthenosphere against the base of the crust. Delamination produced younger thermochronometer ages than an identical exhumation history without it, with fixed Moho variants showing particles heated by roughly 100 degrees Celsius before cooling and ages up to 36 percent younger than equivalent erosion-only models. Yet crucially, those young ages did not record the timing of delamination. Conductive heating from below must dissipate before low-temperature thermochronometers can close, pushing the recorded dates millions of years away from the event itself. For researchers who have linked inflections in exhumation rates to delamination events, this is a cautionary finding: the timing of such deep-seated processes may be difficult or impossible to extract from thermochronometer data alone.

The study also identified a diagnostic signature worth celebrating. Rapidly decaying exhumation scenarios, in which aggressive erosion strips the crust quickly, produced tightly clustered to overlapping ages across all four thermochronometers and were largely insensitive to Moho position or bottom-up heating. The authors infer that closely spaced ages from multiple systems are most likely evidence of an erosion-dominated system rather than a thermal event like delamination or fault-related heat transfer. In contrast, even a perfectly known thermal history can be derived from multiple non-unique geological histories, as demonstrated by four model variants sharing an identical exhumation path yet producing cooling ages differing by up to 30 percent.

The authors are candid about the limitations of their approach. Collapsing three-dimensional geology into one dimension means topographic effects on near-surface isotherms, horizontal heat advection during faulting, sediment compaction, and magmatic processes are not simulated, and the code is best applied to regions with limited relief and moderate to slow exhumation rates. More sophisticated two- and three-dimensional tools such as Pecube remain necessary for detailed analyses in complex settings. Still, the message is clear and broadly applicable: the depth at which the Moho sits, the radiogenic character of the crust, and the timescales of conductive heat transfer all shape the thermal histories recorded in minerals. Simple models like Tc1D, freely available to the community, offer an accessible bridge between the thermal histories produced by inversion software and the geological processes geologists ultimately want to reconstruct, and they suggest that interpreting cooling ages demands as much attention to the crust’s thermal architecture as to its kinematic history.

Subject of Research: One-dimensional modeling of how lithospheric-scale geological processes affect crustal thermal histories and low-temperature thermochronometer ages

Article Title: Simplified modeling of the impact of lithospheric-scale geological processes on thermal histories and low-temperature thermochronometers

Article References: Kellett, D. A., & Whipp, D. M. (2026). Simplified modeling of the impact of lithospheric-scale geological processes on thermal histories and low-temperature thermochronometers. Geochronology, 8(3), 423-445. https://doi.org/10.5194/gchron-8-423-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-423-2026

Keywords: thermochronology, thermal history modeling, lithosphere, erosion, sedimentary burial, faulting, delamination, Moho, apatite (U-Th)/He, zircon fission-track, geothermal gradient, Tc1D