Permafrost is conventionally defined as soil or rock that has remained at or below 0 degrees Celsius for at least two consecutive years, and it blankets more than 20 percent of the northern hemisphere’s land area. Yet the definition conceals a subtlety that has fascinated cryospheric scientists for decades: permafrost is not necessarily completely frozen. Salts dissolved in pore water depress the freezing point, and capillary forces plus adsorption onto fine-grained mineral surfaces allow thin films of liquid water to persist at temperatures well below zero. Those unfrozen films can become highways for ions and fluids, permitting redistribution and migration of water within the cryotic pore space. Near-surface saline permafrost, in which soluble salts are abundant enough to significantly lower the freezing point, may cover as much as 35 percent of the continuous permafrost zone of the circumpolar north.

The research team drilled three cores at the ADE site in Adventdalen, a U-shaped fjord valley carved by glaciers and later flooded by seawater during the Early Holocene. Two boreholes spaced half a meter apart were drilled in March 2017 and treated as a single core, ADE-17, while two more, ADE-1 and ADE-2, were drilled in March 2022 within roughly 30 meters of the first. The valley fill at the site consists of fine aeolian deposits over fluviatile gravel and roughly twelve meters of deltaic sediments laid down when the sea reached this far up the valley, all resting on Lower Cretaceous shales. After the sea retreated following rapid glacial rebound, the exposed marine sediments froze from the top down, while later fluviatile and aeolian layers froze in place as they accumulated, creating a shallow syngenetic cap above the deeper epigenetic permafrost.

What the cores contained stunned the researchers. Ground ice chemistry differed dramatically between boreholes separated by no more than a few tens of meters. In ADE-17 and ADE-1, the shallow syngenetic ice was essentially fresh water, while the deeper epigenetic sections were brackish to saline, with chloride concentrations between 450 and 3230 milligrams per liter. ADE-2 was another story entirely: salinity reached 7960 milligrams of chloride per liter, with one sample exceeding 23,000 milligrams per liter and total dissolved solids above 40,000 milligrams per liter, saltier than seawater. Even the syngenetic section of ADE-2, deposited only two to four thousand years ago, carried total dissolved solids above 6000 milligrams per liter, a clue that something had pushed saline fluid upward into sediments that had never known the sea.

The ionic ratios told a story of two different waters. In ADE-17 and ADE-1, sodium-to-chloride and sulfate-to-chloride ratios far exceeded those of seawater, indicating that porewater had undergone substantial water-rock interaction, dissolving minerals or exchanging ions during long residence in the valley-fill sediments. In ADE-2, by contrast, the ratios closely matched seawater, with sodium-to-chloride values between 0.84 and 1.02, suggesting a marine-derived saline component that had experienced little chemical modification on its journey. Because freezing and salt expulsion would have produced sodium-to-chloride ratios lower than seawater rather than higher, the team ruled out brine formation by freezing, pointing instead to an intrusion of saline fluid that mixed with the original, fresher ground ice.

The decisive evidence came from the decay chains of uranium and thorium. Radium has four naturally occurring isotopes with half-lives spanning an extraordinary range: radium-223 and radium-224 decay within days, radium-228 within 5.75 years, and radium-226 within 1601 years. All are produced by the alpha decay of thorium isotopes, which bind tightly to mineral grains, while radium is comparatively mobile and tends to enter the dissolved phase. Alpha recoil during decay ejects radium from the mineral lattice into pore space. Because short-lived isotopes reach steady-state activities within weeks, the ratio of a long-lived isotope such as radium-226 to a short-lived one such as radium-223 depends almost entirely on how long the ice or fluid has sat in place, growing toward equilibrium with its parent over five to six half-lives, or roughly eight to ten thousand years for radium-226.

In ADE-17 and ADE-1, the long-to-short-lived radium ratios were high, mostly between 20 and 70, close to or above the secular equilibrium value of 21.7 and approaching the parent ratios measured in the bulk sediment. This implies that the ground ice in these cores is ancient, potentially as old as the Early Holocene permafrost that formed roughly 9000 to 9500 years ago after the sediments emerged from the sea. The team attributes the elevated long-lived radium activities to slow diffusion out of the mineral grains themselves, likely through partially liquidized nano-pores, where thin films of unfrozen water along grain boundaries allow ions to escape over years to decades. Simulations using an effective diffusion coefficient of 10 to the minus 13 square meters per second showed that radium could escape silt-sized grains within years by this pathway, enriching the pore space in radium-226 and radium-228 without requiring any contribution from the short-lived isotopes.

ADE-2 defied this pattern. Its long-to-short-lived radium ratios were strikingly low, between 4 and 14, and matched the parent isotope ratios found on grain surfaces in the exchangeable fraction of the sediment rather than in the bulk grains. In other words, the radium clock in this core had been reset recently, before the longer-lived isotopes had any chance to diffuse out from inside the grains or accumulate significantly from recoil. Modeling the buildup of the radium-226 to radium-223 ratio showed that with no adsorption, reaching equilibrium would take several thousand years, but with adsorption coefficients typical of seawater salinity, the timescale collapses to centuries or decades. The team concludes that the saline ground ice in ADE-2 is younger than 100 years, and possibly only years to decades old, the product of a Late Holocene, and likely very recent, intrusion of saline fluids that mixed with the original ground ice in a fingering style.

Where did those fluids come from? A direct connection to Adventfjorden roughly ten kilometers away seems implausible given the low hydraulic gradient and shale-dominated geology, which would produce residence times of thousands of years. The researchers instead point to a local sandstone aquifer beneath the permafrost, most likely the Helvetiafjellet Formation more than 70 meters below the surface, or possibly the deeper fractured Geerdalen Formation. Fractured sandstone typically carries a high water-to-rock ratio and low thorium content, which would explain both the seawater-like chemistry of the intruding fluid and its capacity to reset the radium clock through enhanced desorption. Artesian discharge through pingo structures nearby demonstrates that water does ascend from the sub-permafrost zone in this region, driven by hydraulic heads from surrounding mountains or by basal permafrost aggradation, and the same pressure may have forced saline fluids upward into the shallow permafrost at the ADE site.

The findings carry weight well beyond one Arctic valley. They demonstrate that saline permafrost is heterogeneous, dynamic, and permeable to fluid migration on human timescales, not the inert deep-frozen archive it is often assumed to be. As Arctic air temperatures continue to climb, thawing of saline permafrost may mobilize subsurface brines, threatening infrastructure built on frozen ground and potentially accelerating the release of stored carbon to the atmosphere. The same radioactive clocks that revealed a decades-old intrusion in Adventdalen now offer scientists a tool to detect such hidden flow elsewhere in the circumpolar north, turning the quiet chemistry of decaying atoms into an early-warning system for a warming Arctic.

Subject of Research: Ground ice residence times and intra-permafrost fluid flow in Svalbard permafrost traced with radium and thorium isotopes

Article Title: The history of ground ice formation and intra-permafrost fluid flow as documented by Ra and Th isotopes

Article References: Rotem, D., Weinstein, Y., Harlavan, Y., Torfstein, A., & Christiansen, H. H. (2026). The history of ground ice formation and intra-permafrost fluid flow as documented by Ra and Th isotopes. The Cryosphere, 20(9), 5303-5326. https://doi.org/10.5194/tc-20-5303-2026

Image Credits: AI Generated

DOI: 10.5194/tc-20-5303-2026

Keywords: permafrost, ground ice, radium isotopes, thorium isotopes, Svalbard, Adventdalen, saline permafrost, fluid migration, geochronology, cryosphere, Arctic warming, hydrogeology

The technique at the heart of the study is post-infrared infrared stimulated luminescence, measured at 200 degrees Celsius, known as post-IR IRSL200. When feldspar grains are buried in sediment, background radiation from their surroundings and from within the crystals themselves steadily accumulates trapped electrical charge. Exposing the grains to light releases that charge as a faint glow, effectively resetting the clock. For luminescence dating to yield accurate ages, this resetting — called bleaching — must be complete before burial. In rivers, however, grains can travel in turbid water, in landslides, or in shadowed pulses, and sunlight may only partially wipe the signal. The leftover dose at the moment of burial, called the remnant dose, can inflate calculated ages, sometimes dramatically for young deposits.

Previous work on this problem relied mostly on multi-grain measurements, in which thousands of crystals are stimulated together and their signals averaged. That approach risks masking the true diversity of the grain population: well-bleached, partially bleached, and saturated grains all blur into one number. The new research, led by Arindam Biswas of the University of Cologne and colleagues at Toulouse, took the systematic single-grain route that had never before been applied to modern fluvial analogues at this scale. The team measured up to 2,900 individual grains per sample for equivalent dose determination and up to 2,200 grains for residual dose analysis, mounting crystals in discs holding a ten-by-ten grid of 300-micrometre holes and stimulating each with an 830-nanometre infrared laser on an automated Risø reader.

The samples came from eleven catchments draining perpendicularly off the Andean range, with outlets at the foot of the main Cordillera to minimise the confounding effects of downstream storage and reworking. This positioning lets floodplain sediments record the millennial-scale erosion processes shaping the landscape. The region is a natural laboratory: active tectonics, volcanism, glaciers, steep slopes, and a pronounced climatic gradient from desert to Mediterranean conditions all vary across the study area, offering a wide range of lithologies and environments in which to test what controls bleaching.

In the laboratory, the team ran controlled bleaching experiments on four samples chosen for their diverse catchment characteristics and higher yields of luminescence-sensitive grains. Grains were first reset under a Hönle Sol2 solar simulator, then given a fixed regenerated dose of 30 grays before being exposed to simulated sunlight for durations ranging from one minute to 30,000 minutes. The results showed a strikingly uniform pattern at the sample-average level. After just one minute of exposure, the signal dropped by about 24 percent. By 100 minutes, roughly 79 percent of the initial signal had been depleted, and after 1,000 minutes only 10 percent remained. At 2,880 minutes — two full days — the signal stabilised at approximately 4 percent of its starting value, equivalent to about 1.2 grays, marking a bleaching plateau.

Individual grains, however, told a more complicated story. When the researchers classified grains into fast-, medium-, and slow-bleaching groups based on their response after one minute of light, all three populations behaved distinctly at first. Yet, remarkably, all converged toward the same asymptotic level between two days and 30,000 minutes of exposure. This convergence suggests a fundamental limit to signal reduction rather than insufficient bleaching time. Even more intriguing, prolonged light exposure homogenised the grain population: the interquartile range of normalised signal values shrank four- to six-fold, from about 0.16 to 0.18 at one minute down to 0.03 to 0.05 at 30,000 minutes. Extended sunlight, in other words, drives diverse grains toward increasingly uniform signal levels — a dose-homogenising effect with practical implications for how long sediments must be exposed to reset their clocks reliably.

What explains the variability in residual doses, which ranged from near zero to about 23 grays among individual grains and from 2.44 to 9.49 grays on average across samples? The team tested three candidate culprits. First, grain chemistry: using an electron microprobe, they measured major oxide concentrations in 121 individual grains, finding potassium oxide levels spanning from 0.20 to 16.60 weight percent. Despite this wide compositional range, no statistically significant correlation emerged between any major oxide — potassium, sodium, calcium, iron, barium, silicon, or aluminium — and either residual dose magnitude or bleaching rate. Mineral composition, including the potassium content that affects internal radiation dose, simply does not control how efficiently a grain bleaches. Second, catchment lithology: samples from catchments with similar dominant rock types showed markedly different doses, while samples from geologically different basins sometimes matched, ruling out provenance as a direct control.

The third candidate proved decisive. When the team compared the modal equivalent dose — the peak of the single-grain dose distribution, representing the best-bleached grain population — with the modal residual dose measured after two days of laboratory bleaching, a strong positive linear correlation emerged, with a coefficient of determination of 0.89. Grains carrying larger accumulated natural doses are progressively harder to bleach. Crucially, the intercept of the linear fit was indistinguishable from zero within uncertainty, implying a negligible unbleachable component in these samples — a contrast with earlier studies that reported small irreducible doses. The youngest samples, with low natural remnant doses, could plausibly reach zero residual dose, while older samples carrying more than 10 grays of remnant dose retained significant residuals even after two days of intense simulated sunlight.

These findings carry direct consequences for how scientists correct palaeodoses when dating sedimentary archives. The team found that remnant doses in their modern analogues, ranging from 0.88 to 25.50 grays, systematically exceeded both laboratory residual doses and the unbleachable component — meaning that subtracting a laboratory residual from a palaeodose would undercorrect the true inherited dose, potentially by 0.5 to 4 grays in the best cases and more than 20 grays in the worst. Subtracting the remnant dose of a single randomly chosen modern analogue is equally unreliable given the scatter. Instead, the authors evaluate three context-sensitive strategies: making no correction, using a distribution-based remnant dose estimate from multiple modern analogues, or applying the unbleachable component. Their recommendations depend on depositional setting — unbleachable-component corrections may suit well-bleached aeolian and coastal deposits, while fluvial and glacio-fluvial settings demand modern-analogue approaches — and on measurement protocol, aliquot size, and the statistical age model applied.

Beyond dating, the results open a path for luminescence-based sediment tracing, a technique that uses incomplete bleaching as a fingerprint of how sediment moved through a landscape. Identifying well-bleached grains is essential for tracing sediment pathways, and previous studies have used fixed thresholds, such as a 10-gray cutoff, to classify them. The new work suggests a more refined, sample-specific approach: classify grains as well-bleached if their natural dose falls within two standard deviations of the mean of that sample’s single-grain residual dose distribution, provided a clear bleaching plateau has been established. Though time-intensive and still requiring validation in field applications, this method promises more realistic thresholds tailored to each geomorphic setting. From the arid slopes of the high Andes to the rain-soaked basins of the south, the humble feldspar grain is proving to be not just a clock, but a recorder of its own journey through sunlight and stone.

Subject of Research: Bleaching behaviour of the post-IR IRSL200 luminescence signal in single-grain K-feldspars from modern fluvial deposits in the Southern Central Andes, Chile

Article Title: Novel insights into the post-IR IRSL200 signal bleachability of single-grain K-feldspars in fluvial modern analogues from the Southern Central Andes, Chile

Article References: Biswas, A., Riedesel, S., Karman-Besson, L., Hellers, M., Guyez, A., Bonnet, S., & Reimann, T. (2026). Novel insights into the post-IR IRSL 200 signal bleachability of single-grain K-feldspars in fluvial modern analogues from the Southern Central Andes, Chile. Geochronology, 8(2), 297-312. https://doi.org/10.5194/gchron-8-297-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-297-2026

Keywords: luminescence dating, K-feldspar, post-IR IRSL, single-grain analysis, bleaching, residual dose, fluvial sediment, Southern Central Andes, Chile, sediment tracing, geochronology, modern analogues

For decades, the Christiana Archipelago has been the neglected sibling of the volcanic field that contains Santorini, one of the most intensively studied volcanoes on the planet. Only two scientific papers on the islands had appeared in the past half century. Because the chemistry of Christiana’s lavas resembles the oldest exposed rocks on Santorini’s Akrotiri peninsula, dated to about 600,000 years ago, researchers simply assumed the two were contemporaneous. Other studies took a different route, correlating submarine volcanic deposits with Pliocene sediment layers in seismic reflection profiles, which suggested ages of roughly 1.6 to 1.7 million years. Without isotopic ages from the rocks themselves, the true age of the archipelago remained an open question.

To settle the matter, Pieter Z. Vroon of VU University Amsterdam and colleagues collected twenty-two rock samples weighing one to three kilograms each from sixteen outcrops along the coasts and higher ground of Christiani island. Weathered rinds were sliced away with diamond saws, and after careful petrographic screening twelve samples had to be discarded because their groundmass had been too heavily altered by time and seawater. The remaining ten samples were crushed, sieved and purified using heavy-liquid centrifuge separation and magnetic sorting until only fresh, inclusion-free groundmass grains remained, ready for one of the most precise dating techniques in geology.

The method of choice was argon-argon dating, a refinement of the classic potassium-argon clock. Volcanic rocks contain potassium-40, which decays slowly into argon-40, a gas that becomes trapped in the rock’s crystal lattice the moment the lava cools. In the argon-argon variant, samples are bombarded with neutrons in a reactor, converting potassium-39 into argon-39 so that both parent and daughter isotopes can be measured from the same grain. The team irradiated their samples at Oregon State University’s TRIGA reactor alongside the Fish Canyon Tuff sanidine standard, then step-heated them with a carbon dioxide laser through ten to eighteen incremental temperature steps, analysing the released gas on a high-sensitivity multi-collector mass spectrometer.

The ages that emerged were remarkably consistent. Five reliable samples from the island’s Upper Lava formation clustered tightly between 2.57 and 2.69 million years old, with analytical uncertainties of as little as 20,000 to 30,000 years, and a weighted mean of 2.63 plus or minus 0.05 million years. Inverse isochron calculations, which correct for any trapped atmospheric or excess argon, agreed with the plateau ages within error. One sample gave a slightly older and less precise age of about 2.23 million years, possibly reflecting minor argon loss, while three samples were rejected outright on statistical grounds following established quality-control protocols for argon-argon data.

Then came the surprise. A single obsidian sample collected from the island’s east coast yielded an age of just 133,000 plus or minus 14,000 years, a million times younger than its neighbours. The researchers concluded that this glassy clast, a lithic fragment embedded in a pyroclastic deposit, almost certainly did not erupt on Christiana at all. The volcano has likely been extinct for some 1.5 million years. Instead, the obsidian was probably ripped from older rock, completely reset by the searing heat of a passing pyroclastic flow, and carried to Christiani by one of Santorini’s great Plinian eruptions. Its age overlaps with the Middle Pumice eruption of approximately 145,000 years ago, one of at least eleven large explosive events that Santorini produced between 350,000 years ago and the famous Late Bronze Age catastrophe.

The new dates transform our picture of how the Christiana-Santorini-Kolumbo volcanic field evolved. The Christiana volcano, a conical edifice rising some 500 metres from the seafloor with an estimated minimum volume of 25 cubic kilometres, was built during a period when the regional stress field was stretching the crust in a north-northeast to south-southwest direction. When that stress field rotated toward a northwest-southeast orientation between roughly 2.5 and 1.5 million years ago, Christiana died. What followed was an extraordinarily long volcanic silence of more than a million years, until the submarine Poseidon volcano and the Kolumbo seamount chain, along with the ancestors of modern Santorini, finally reawakened the rift around 1.2 million years ago.

Perhaps the most consequential implication reaches far beyond one forgotten island. Comparing the new ages with published data from the other four volcanic fields of the South Aegean Volcanic Arc, from Sousaki and Aegina-Poros-Methana in the west to Milos and Kos-Nisyros-Yali in the east, the authors show that all five fields were simultaneously active around 3 million years ago. This synchronised ignition coincides with the moment that subducted oceanic crust of the African plate first reached a depth of about 100 kilometres beneath the arc, the depth at which dehydration of the slab triggers melting in the overlying mantle wedge. Before that, the subduction zone had been grinding down continental material for tens of millions of years without producing the modern arc volcanoes.

The study also reveals a striking arc-wide pattern of volcanic famine. Between roughly 2 and 1 million years ago, every volcanic field along the arc experienced quiescence or sharply reduced output, even as the orientation of the tectonic stress field was changing in each of them. The authors argue that this synchronised lull reflects the difficulty parental magmas faced in reaching the surface while the stress regime was reorganising, rather than any fundamental change in the subduction system itself. Where volcanoes erupt, and when they fall silent, appears to be dictated by the local architecture of faults and extension, a reminder that even the largest volcanic systems are prisoners of the crust above them.

For a volcanic field that includes one of the most notorious caldera-formers in human history, pinning down its deep past is more than an academic exercise. The work demonstrates that the Christiana-Santorini-Kolumbo rift has hosted at least two distinct volcanic eras separated by a million-year gap, and that the tectonic switch between them set the stage for everything from the giant 89-cubic-kilometre Archaeos Tuff of 520,000 years ago to the Minoan eruption that may have helped end Minoan civilisation. The three lonely islands that started it all, meanwhile, now hold a secure place in the story: not young cousins of Santorini, but the ancient founding fathers of one of the Mediterranean’s most dangerous volcanic neighbourhoods.

Subject of Research: Argon-argon geochronology of the Christiana volcano in the South Aegean Volcanic Arc, Greece

Article Title: 40Ar ∕ 39Ar constraints on the eruption history of the Christiana Volcano of the Christiana-Santorini-Kolumbo volcanic field, Greece

Article References: Vroon, P. Z., Beemster, T., Zhou, X., Nomikou, P., Klaver, M., Wijbrans, J. R., & Kuiper, K. F. (2026). 40 Ar ∕ 39 Ar constraints on the eruption history of the Christiana Volcano of the Christiana-Santorini-Kolumbo volcanic field, Greece. Geochronology, 8(2), 279-295. https://doi.org/10.5194/gchron-8-279-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-279-2026

Keywords: Christiana volcano, Santorini, argon-argon dating, geochronology, South Aegean Volcanic Arc, volcanology, Kolumbo, pyroclastic deposits, subduction, Aegean Sea, obsidian, tectonics

Biotite, the dark, flaky mica that gives many metamorphic rocks their shimmer, is prized by geochronologists because it readily incorporates rubidium into its crystal lattice. Rubidium-87 decays slowly to strontium-87, so measuring the ratio of parent to daughter isotopes in biotite has, for decades, been a standard way to date the cooling of mountain belts. The traditional framework, built on the concept of closure temperature introduced by Martin Dodson in 1973, holds that each mineral system effectively stops exchanging isotopes with its surroundings below a characteristic temperature. For biotite’s rubidium–strontium system, community consensus has hovered somewhere between roughly 300 and 500 degrees Celsius, a range derived largely from studies that crushed rocks, separated out the mica, and analyzed bulk powders.

The problem with those conventional approaches is that they destroy the very information that matters most: the textural context of each grain. McFarlane’s study exploits a technological leap, laser ablation combined with tandem inductively coupled plasma mass spectrometry, or LA ICP-MS/MS, which allows rubidium and strontium isotopes to be measured in situ, on spots as small as tens of micrometers, while the grain remains embedded in its thin section. Using sulfur hexafluoride as a reaction gas in a triple-quadrupole mass spectrometer, the instrument chemically separates strontium from rubidium inside the machine itself, converting strontium ions to strontium fluoride species that can be measured free of the isobaric interference from rubidium-87 that plagued earlier microanalytical attempts.

The setting could hardly be better suited to the task. The Makhavinekh Lake Pluton, part of the Nain Plutonic Suite, is a 30-kilometer-wide composite intrusion assembled around 1322 million years ago from troctolite, anorthosite, and ferrodiorite, ringed by granite. Its host rocks, the Tasiuyak Gneiss, are migmatitic metasediments that had already experienced intense regional metamorphism around 1850 million years ago, at temperatures near 850 degrees Celsius and pressures of 6 to 9 kilobars. When the pluton intruded, it reheated the gneiss in a graded thermal halo, or aureole, reaching more than 800 degrees Celsius within meters of the contact but only about 550 degrees Celsius some 6 kilometers away. Numerical thermal models show that rocks in the outer aureole took roughly 15 million years to cool back below 500 degrees Celsius, while rocks nearest the contact needed more than 25 million years. Crucially, the timing of the contact event is pinned independently at 1322 million years by uranium–lead dating of zircon and of new monazite that grew during the heating pulse.

What McFarlane found in the outer aureole, more than 4 kilometers from the contact, was striking. Biotite grains that had been completely enclosed, or armoured, inside large garnet crystals preserved rubidium–strontium isochron ages spanning from about 1850 million years down toward the 1322-million-year mark of the intrusion. This spread defines what geochronologists call a zone of partial retention, a statistical signature visible as an S-shaped distribution on a linearized probability plot. The key to the pattern lies in the garnet itself: microscopic fractures cutting through the host garnet acted as short-circuit pathways, allowing radiogenic strontium to leak out of connected biotite inclusions and into newly formed plagioclase in the contact metamorphic assemblages. Inclusions that happened to be isolated from these microfractures, sealed off in intact garnet, retained their full 1850-million-year regional metamorphic age, with a filtered subset of analyses yielding an anchored isochron of 1824 plus or minus 32 million years, overlapping independently known uranium–lead ages for the regional event.

Moving inward, into the central aureole between roughly 3.7 and 1.1 kilometers from the contact, the story intensifies. Here the regional garnet is progressively eaten away by diffusion-controlled coronas of orthopyroxene and cordierite, minerals that grew as the garnet reacted with the surrounding matrix during the thermal pulse. Biotite grains that were once armoured but later intersected by these coronal assemblages show ragged, reacted boundaries decorated with tiny ilmenite grains. Only fully armoured inclusions still record regional ages, in this case an anchored isochron of 1870 plus or minus 25 million years. Biotite that was partly or fully liberated into the coronas yields ages clustered near 1330 million years, and these reset domains show even higher rubidium-to-strontium ratios than their outer-aureole counterparts, because strontium loss was more extreme at the higher temperatures involved. The youngest fully liberated grains give an anchored age of 1317 plus or minus 19 million years, consistent with thermal models predicting cooling below 500 degrees Celsius after more than 20 million years.

In the inner aureole, within about 100 meters of the contact, garnet and sillimanite were completely consumed, and a new generation of biotite, called neoblasts, grew in textural equilibrium with the orthopyroxene and cordierite assemblages. These neoblasts tell a subtler tale. Their pooled isochron gives 1298 plus or minus 15 million years, but individual grains display dramatic internal variability in rubidium-to-strontium ratios, sometimes varying by nearly a factor of two across a single plate less than 200 micrometers across. That heterogeneity means single grains can yield their own internal isochrons, and six such grains produced a weighted mean age of 1294 plus or minus 14 million years. Yet McFarlane cautions against taking these numbers at face value. The initial strontium isotope composition of each neoblast likely depended on the precursor mineral in the microscopic domain where it nucleated: old regional biotite had accumulated highly radiogenic strontium over 530 million years, while low-rubidium phases like garnet, apatite, and plagioclase had barely evolved at all. Some neoblasts may have inherited radiogenic strontium, and slow intergranular diffusion in the dry, granulite-facies environment may have armoured others from further exchange, making a simple averaged isochron age potentially meaningless.

Perhaps the most evocative textures come from the immediate contact zone, where the Tasiuyak Gneiss partially melted. Biotite and potassium feldspar form cuspate, interstitial films with low dihedral angles wrapping rounded quartz grains, the classic signature of crystallized melt, so-called melt pseudomorphs. These biotite films show a much narrower spread of rubidium-to-strontium ratios than the neoblasts, which McFarlane attributes to the three-dimensional connectivity of melt pockets and the fast diffusion of rubidium and strontium through silicate liquid, homogenizing the system at the centimeter rather than micrometer scale. Their anchored isochron of 1329 plus or minus 17 million years sits within uncertainty of the pluton’s crystallization age, though the data hint at open-system behavior and an initial strontium ratio above 0.80, possibly reflecting preferential melting of radiogenic biotite and potassium feldspar, a phenomenon also reported from Himalayan leucogranites.

The upshot is a fundamental challenge to the closure-temperature paradigm as applied to biotite. Strontium diffusion in biotite was demonstrably efficient at temperatures of 550 degrees Celsius and above, on timescales of just a few million years, placing the effective closure temperature at the upper end of, or above, the classic 300-to-500-degree estimates. But there is no single number to quote, because armoured grains resisted resetting even where their liberated neighbors were fully re-equilibrated. Strontium diffuses roughly eight orders of magnitude more slowly in garnet than in biotite at 600 degrees Celsius, so a garnet host can act as an essentially impermeable shield. The practical consequence is that biotite rubidium–strontium dating, far from being unreliable in high-grade rocks, may actually recover prograde metamorphic ages in upper-amphibolite and granulite-facies terranes, provided analysts systematically hunt for armoured grains and recognize the diagnostic zones of partial retention that require large datasets, often more than a hundred spots, to reveal.

McFarlane points toward an even more tantalizing frontier: ultrahigh-temperature rocks, where biotite inclusions armoured in peritectic minerals such as cordierite, sapphirine, orthopyroxene, and garnet record dehydration reactions at conditions exceeding 900 degrees Celsius. Because armoured biotite ages can be tied directly to prograde reactions in the major mineral assemblage, unlike uranium–lead ages from accessory zircon or monazite, they could offer an entirely new window into the deepest thermal episodes in Earth’s crust. For now, the message from the Labrador aureole is clear: before trusting any mineral date, look at the rock, grain by grain, because the architecture of a crystal’s prison determines whether its memory of deep time survives.

Subject of Research: Rb–Sr isotope systematics of biotite during high-temperature contact metamorphism

Article Title: Response of the Rb–Sr system in biotite during contact metamorphism in the aureole of the Makhavinekh Lake Pluton, Labrador

Article References: McFarlane, C. R. M. (2026). Response of the Rb–Sr system in biotite during contact metamorphism in the aureole of the Makhavinekh Lake Pluton, Labrador. Geochronology, 8(2), 313-327. https://doi.org/10.5194/gchron-8-313-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-313-2026

Keywords: biotite, Rb–Sr geochronology, contact metamorphism, closure temperature, laser ablation ICP-MS/MS, Makhavinekh Lake Pluton, Labrador, Tasiuyak Gneiss, garnet, strontium diffusion, metamorphic rocks, isotope dating

The physics behind surface-exposure dating is elegant. High-energy cosmic rays striking the upper atmosphere shower the Earth’s surface with secondary particles that split atomic nuclei within minerals such as quartz. One product of this spallation is beryllium-10, a rare isotope that accumulates in the top few centimetres of an exposed rock at a known rate. Measure the concentration of beryllium-10 in a glacially polished boulder, divide by the production rate, and you obtain the time since the ice retreated. The catch is that the production rate is not constant across the planet. It varies with altitude, latitude, atmospheric pressure and the strength of Earth’s geomagnetic field, which shields the surface from cosmic radiation. Calibration experiments therefore measure beryllium-10 in surfaces whose true exposure age is known independently, typically from radiocarbon dating, and then use numerical scaling models to extrapolate the site-specific rate to anywhere else on Earth. Most modern calibrations converge on sea-level, high-latitude values of roughly 3.8 to 4.1 atoms per gram of quartz per year, but the residual spread among them is large enough to matter for short-lived events.

Bromley and his colleagues chose the Redpoint Peninsula as their natural laboratory, a low-relief coastal landscape between Loch Torridon and Loch Gairloch underlain by Neoproterozoic Torridonian sandstones. The peninsula carries a conspicuous suite of moraine ridges, first described in 1979, that record a stabilisation or readvance of the decaying last British ice sheet during Termination 1, the final deglaciation that ended the last glacial period. This episode, known as the Wester Ross Readvance, has been dated repeatedly over the years with conflicting results, with published beryllium-10 ages ranging from about 13.5 to 16.3 thousand years depending on which production rate was applied. The key to resolving the ambiguity lay in a small lake called Loch Bad na h-Achlaise, sitting just inside the mapped readvance limits. Sediment cores from the lake basin record the transition from marine to freshwater conditions as the ice retreated and the land rebounded, and the basal organic sediments yielded radiocarbon ages showing the peninsula was ice free by at least 15.9 thousand years ago.

That radiocarbon benchmark is what gives the new study its power. Because the boulders sampled for beryllium-10 analysis all lie above the local marine limit and close to the lake, the onset of organic sedimentation provides a minimum-limiting age for their exposure: the rocks must be at least as old as the lake sediments, and probably somewhat older, since vegetation needed time to colonise the freshly deglaciated basin. The team collected samples from the top few centimetres of eleven glacial boulders at three sites, including perched boulders on ice-moulded bedrock highs and boulders embedded in the crest of the outermost Redpoint moraine. Back in the laboratory, quartz was painstakingly isolated through acid leaching and froth flotation, beryllium was extracted by ion chromatography, and isotope ratios were measured by accelerator mass spectrometry at Lawrence Livermore National Laboratory. Procedural blanks contributed only about one percent of typical sample signals, and two statistical outliers were identified and set aside, leaving nine tightly clustered measurements.

The internal consistency of the nine remaining ages is striking. Treated as a single population, they form a normal distribution whose scatter can be explained by analytical uncertainty alone, with a low reduced chi-squared value. That statistical coherence means the landforms sampled, from terrain outboard of the moraines to the moraine crest itself, were all exposed within a very narrow window of time. It also means the dataset is an ideal yardstick for testing production rates: whatever rate is applied, the relative pattern stays the same, so any systematic offset from the radiocarbon benchmark points directly at the calibration rather than at the geology. Using version 3 of the University of Washington online calculator and three standard scaling schemes, the team calculated apparent exposure ages with eight production-rate datasets drawn from Scotland, Switzerland, Sweden, Norway, north-eastern North America and a global compilation.

The results were decisive. Only two of the eight calibrations, Rannoch Moor and the Swiss Chironico landslide, produced mean ages overlapping the radiocarbon control within one standard deviation, and only Rannoch Moor yielded a mean age, 16.1 thousand years, that actually predates the onset of lacustrine sedimentation. The remaining six rates produced ages younger than the radiocarbon benchmark, under-predicting the true deglaciation age by up to seven to eight percent. The worst performers shared a common weakness: their calibration surfaces lack direct, independent age control. The Isle of Skye and Highlands dataset, for example, relies on presumed ages for the deglaciation of cirques correlated with distant climate records, while the Glen Roy calibration ties a wave-cut shoreline to a varve sequence twenty-five kilometres away without unequivocal dating. The Mount Billingen calibration assumes bedrock erosion features correspond to the indirectly dated drainage of the Baltic Ice Lake. Even the widely used Borchers global compilation, the default in popular online calculators, technically passed a statistical test only because of its large internal scatter, and the authors warn it too would yield ages that are stratigraphically too young.

The study also delivers a hard number. Using the radiocarbon control as an independent minimum-limiting age, the team derived a maximum-limiting sea-level, high-latitude beryllium-10 production rate of 3.925 plus or minus 0.070 atoms per gram per year under time-independent scaling. Any production rate higher than this ceiling will produce exposure ages that are unrealistically young relative to the radiocarbon chronology. The finding carries weight beyond Scotland because it independently vindicates the Rannoch Moor calibration, which had come under challenge. A 2019 critique questioned the radiocarbon ages of moss fragments underpinning the calibration, suggesting inbuilt age errors from aquatic carbon or carbonate contamination. Bromley and colleagues systematically dismantle that critique, noting that the re-sampling campaign examined only one of seven original sites, that the new basal ages from that site are virtually indistinguishable from the originals, and that the carbon-isotope values of the dated mosses fall squarely within the range seen in other well-validated British radiocarbon datasets. No carbonate sources exist up-ice of the original coring locations.

Beyond the metrology, the new chronology rewrites the story of the Wester Ross Readvance itself. Because the nine exposure ages from terrain on both sides of the Redpoint moraines are statistically indistinguishable, the ice-margin stabilisation they record must have been brief, lasting no longer than roughly two hundred years. The readvance, previously envisaged by some workers as a spatially extensive resurgence of ice, now appears as a short-lived hiccup in the wholesale collapse of the British ice sheet during Heinrich Stadial 1, a cold interval between about 17.8 and 14.6 thousand years ago. The timing also places the event at least five hundred years earlier than previous estimates. Brief but prominent moraine-building episodes at around this time appear in records from elsewhere in Britain and Ireland and as far afield as New Zealand and South America, hinting at a widespread, if subtle, glacial response to North Atlantic climate variability during the deglaciation.

The broader lesson is one that resonates across the geosciences: a dating method is only as good as its calibration. Surface-exposure dating underpins reconstructions of glacier and ice-sheet behaviour that in turn inform projections of future sea-level rise, so minimising production-rate uncertainty is not an academic nicety but a practical necessity. By anchoring beryllium-10 chronology to an independent radiocarbon benchmark in the same landscape, the Redpoint study shows that calibrations built on directly dated surfaces consistently outperform those resting on inferred ages, and it hands researchers working at mid-latitude Northern Hemisphere sites a validated tool. The boulders of Redpoint, battered by Atlantic gales for sixteen millennia, have finally told the time with precision, and in doing so they have sharpened the clock that scientists around the world use to read the history of a melting planet.

Subject of Research: Calibration and testing of cosmogenic beryllium-10 production rates for surface-exposure dating of the last British ice sheet deglaciation

Article Title: Testing current estimates of the in situ cosmogenic 10Be production rate in the north-western British Isles, with implications for ice sheet behaviour during Termination 1

Article References: Bromley, G. R. M., Hall, B. L., Putnam, A. E., & Lowell, T. V. (2026). Testing current estimates of the in situ cosmogenic 10 Be production rate in the north-western British Isles, with implications for ice sheet behaviour during Termination 1. Geochronology, 8(2), 329-349. https://doi.org/10.5194/gchron-8-329-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-329-2026

Keywords: cosmogenic nuclide dating, beryllium-10, surface-exposure dating, production rate calibration, British ice sheet, Wester Ross Readvance, Termination 1, Heinrich Stadial 1, radiocarbon dating, Scotland, deglaciation, geochronology

The traditional method is elegantly simple. Researchers count the density of tracks crossing a polished mineral surface, measure the uranium concentration, and combine these with the decay constant in a standard age equation. For minerals that cooled quickly, such as ash or volcanic extrusions, this produces a precise crystallisation age. But the equation is also routinely applied to apatite that cooled slowly through the so-called partial annealing window, the temperature range between roughly 120 and 60 degrees Celsius in which tracks form but are gradually shortened and erased. In that regime the resulting single apparent age is at best a rough average of a complicated thermal story, and Jensen argues it deserves far more scrutiny.

The key insight behind the new work is that track length carries information that the conventional equation ignores. Fission tracks are born at a standard initial length of about 16 micrometres and then shorten over time as heat allows the crystal to heal itself, a process called annealing. Crucially, the oldest tracks in a grain have generally experienced the most annealing and are therefore the shortest, while the youngest tracks remain the longest. In principle, then, a track of a given length can be dated simply by counting how many shorter tracks exist below it, adding one, and dividing by the rate at which tracks are generated per unit volume of crystal. That deceptively clean relationship was established theoretically decades ago, but turning it into a practical dating tool has proved difficult.

The difficulty lies in blurring. Tracks of the same true age do not all have the same length. The fission fragments carry a spread of kinetic energies, the crystal lattice shortens tracks at different rates depending on their orientation relative to the crystallographic c-axis, and the etching and measurement procedures used to make tracks visible under the microscope introduce further biases. The result is that each column of a measured track-length histogram mixes tracks from several different time intervals, and the tidy one-to-one link between length and age is smeared out. Without correcting for this smearing, any attempt to assign ages to individual length classes would be meaningless.

Jensen’s solution borrows a mathematical technique from signal processing: deconvolution, the same family of methods used to strip noise from telecommunications signals. The idea is that the measured histogram can be treated as a weighted sum of reference histograms, each representing tracks generated in a short time interval and then annealed to a particular degree. Those reference distributions, called kernels, come from laboratory experiments in which apatite samples are irradiated with neutrons in a reactor and then heated at controlled temperatures and durations, producing known track-length spreads. If the measured histogram can be decomposed into the right weighted combination of these kernels, the weights reveal how many tracks were produced in each time interval, restoring the link between histogram columns and time.

The new study advances this approach in two important ways. First, it applies the deconvolution to c-axis projected tracks, a standardised measurement protocol in which track lengths are corrected along the crystal’s principal axis. Because tracks of the same age tend to cluster along an elliptical curve in the length-versus-angle diagram, projecting them onto the c-axis already gathers contemporaneous tracks together, acting as a partial deconvolution in itself. Second, the inversion is formulated in a fully probabilistic framework based on the work of geophysicist Albert Tarantola, using an extended least-squares method that incorporates both the uncertainties of the measurements and prior information about the expected answer. Because the total number of counted tracks is fixed, the histogram columns are strongly covariant, and the method handles this explicitly through variance-covariance matrices rather than treating each bin as independent.

One stubborn mathematical problem is that least-squares inversion can return negative numbers of tracks, which is physically impossible. Jensen tames this by choosing positive prior values with variances large enough that the priors guide the solution without dominating the data, ensuring that the deconvolved histogram columns remain mainly positive. The variances are also chosen generously enough that, in principle, all tracks could end up in a single bin, so the prior does not impose an artificially rigid shape on the result. The full calculation, implemented in the open-source Octave platform and released with the paper, converts the deblurred length histogram into a histogram of time intervals and then cumulates those intervals from the youngest tracks backwards to assign an age to the oldest track in each column.

To test the method, Jensen applied it to twenty apatite samples from Ellesmere Island and Northwest Greenland, part of a large dataset assembled to reconstruct the Phanerozoic tectonic history of the Arctic. The results are striking. In a Late Devonian sandstone sample, the deconvolution extracted three statistically distinct age nodes at roughly 134, 274 and 443 million years, with non-overlapping uncertainties. The oldest of these predates the deposition of the sandstone itself, revealing that the apatite grains carried inherited tracks from earlier mountain-building episodes before they were eroded and redeposited. By comparing track ages with the known depositional age, the method could even separate post-depositional tracks from inherited ones, showing in this case that about 70 percent of the tracks formed after burial.

The deconvolution also sharpens features that are nearly invisible in raw data. For several samples, histograms that appeared broad and unimodal in the original measurements revealed clear bimodality after deblurring, a signature of distinct thermal episodes. One sample told a three-act story: a long stretch between about 193 and 90 million years ago when maximum temperatures reached the upper part of the partial annealing window, a shorter cooling phase through the lower part of the window, and then a prolonged period of moderate temperatures below 60 degrees Celsius lasting to the present. The method even proved useful as a quality check, flagging samples whose histograms showed unrealistically sharp peaks caused by too few measured tracks, since the forward calculation from the deconvolved histogram could not reproduce such features.

The broader payoff is a new layer of information for reconstructing how mountains rose, basins filled and plates shifted. Across the sample suite, the oldest track ages typically came out about 20 percent older than the conventional central ages, and they broadly align with the timing of the last thermal peaks above 100 degrees Celsius inferred from independent thermal history modelling, albeit with large uncertainties. Most samples yielded between two and five reliable age nodes, averaging three, from a single histogram. Jensen suggests these multiple age estimates could guide the choice of age-temperature constraints in thermal history simulation models, potentially tightening reconstructions of past temperatures and tectonic events. For a technique that has spent half a century squeezing a single number out of each mineral grain, the prospect of reading an entire chronological sequence from one crystal marks a quiet but potentially transformative shift in how geologists read Earth’s thermal diary.

Subject of Research: Extraction of multiple thermal history ages from c-axis projected apatite fission tracks using probabilistic deconvolution

Article Title: Extraction of multiple ages from c-axis projected fission tracks

Article References: Jensen, P. K. (2026). Extraction of multiple ages from c -axis projected fission tracks. Geochronology, 8(2), 373-386. https://doi.org/10.5194/gchron-8-373-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-373-2026

Keywords: fission track dating, apatite, thermochronology, deconvolution, uranium-238 fission, partial annealing zone, c-axis projection, geochronology, tectonic history, inverse theory, track length distribution, Ellesmere Island

The (U-Th)/He dating method rests on a straightforward nuclear clock. Uranium-238, uranium-235, thorium-232, and samarium-147 decay by emitting alpha particles, and each alpha particle is, in essence, a helium-4 nucleus. Over geological time, these nuclei accumulate in minerals such as apatite and zircon. Determining a date requires measuring both the accumulated radiogenic helium and the remaining parent isotopes. The parent isotopes are routinely quantified by inductively coupled plasma mass spectrometry, but that technique cannot touch noble gases. Helium, therefore, must be measured separately on a gas-source mass spectrometer attached to a high-vacuum extraction line, and the quality of that helium measurement directly controls the quality of the resulting date.

Commercial helium analysis lines exist and serve many laboratories well, but the Colorado team faced a familiar research dilemma. Their existing commercial instrument, an Alphachron, had reliably produced whole-grain helium data for more than a decade, yet its proprietary control software made it poorly suited for developing brand-new protocols, and its busy schedule left no machine time for the experimental iterations that laser-ablation work demands. Closed systems also complicate mundane realities: replacing failed hardware, updating software and operating systems to satisfy institutional security requirements, and adapting analytical routines as methods evolve. Rather than fight those constraints, the researchers decided in 2018 to build their own line from scratch, and to document the process in unusual detail because, as they note, published guidance on constructing a helium line is remarkably scarce, and institutional knowledge of this kind tends to vanish when experienced scientists retire.

The central design philosophy of the Jimbochron is ruthless volume minimization. When a laser ablates a minute pit in a crystal, it liberates roughly a femtomole of helium, and that gas expands into whatever volume is available. A smaller expansion volume means a higher helium concentration, a larger signal at the mass spectrometer, and consequently more stable and more precise measurements. Smaller volumes also evacuate faster, driving down the background helium blanks that would otherwise swamp tiny samples. To squeeze the internal volume down, the team arranged pneumatic valves around two hexagonal vacuum manifolds, a configuration that packs many valves closely together with minimal stainless-steel tubing while leaving spare ports for future additions such as a diode laser, extra pumps, or additional getters. The frame itself is built from T-slotted aluminum bars, chosen because it is inexpensive, easy to cut and tap in-house, light enough to mount on wheels, and directly compatible with the mounting holes on the Swagelok pneumatic valve bodies.

Measuring femtomole helium also demands an uncompromising vacuum. Although helium makes up only about five parts per million of dry air, the helium in a mere hundredth of a microliter of air rivals the radiogenic helium in a typical apatite crystal, so atmospheric contamination must be virtually eliminated. The Jimbochron is conceptually split into a process side, where samples are loaded and which is regularly exposed to atmosphere, and an analysis side, which houses the mass spectrometer and almost never sees unconditioned gas. A dry scroll pump backs a turbo pump on the process side, while an ion pump maintains the analysis side, where pressures typically sit below five times ten to the minus ten torr. Baking the entire line at 200 degrees Celsius for roughly 48 hours drives adsorbed gases from the metal surfaces, and all components are rated for bakeout temperatures of at least 225 degrees. Vibrations from the scroll and turbo pumps, which travel readily through rigid aluminum framing, are isolated by mounting the pumps on a separate wheeled cart and connecting them only through flexible bellows.

Gas handling follows a three-stage logic of liberation, conditioning, and analysis. A short-wavelength excimer laser operating at 193 nanometers ablates small regions of a grain, shattering the crystal lattice and releasing trapped helium in situ; a long-wavelength diode laser will later be added for conventional whole-grain heating and continuous ramped-heating experiments. The liberated gas is then cleaned by getters, devices packed with reactive zirconium-based compounds that sequester every non-noble gas, and simultaneously spiked with a known aliquot of helium-3. Finally, the mixed helium-3 and helium-4 charge is analyzed in static mode by a high-sensitivity Hiden 3F quadrupole mass spectrometer, chosen because manufacturer data indicate it is roughly an order of magnitude more sensitive than the quadrupoles common in (U-Th)/He laboratories, and because its triple-quadrupole filters resolve the helium peaks more effectively. A second, less sensitive residual gas analyzer with a wider mass range is reserved for future ramped-heating work that would expose a spectrometer to large gas loads.

Automation was treated as a first-class design goal rather than an afterthought. Custom LabView software, written as modular subroutines and shared through an online repository, controls valves, pumps, gauges, and the laser through a network of USB and serial hubs. Before every measurement the software checks the health of the system, confirming that the mass spectrometer filament is alive, that pressures are below thresholds, and that helium-3 and helium-4 signals sit at proper blank levels, adjusting in real time to changing conditions. Every action, manual or automated, is logged and time-stamped for quality control, and analytical sequences are loaded as text files, which permits runs of more than a hundred analyses, including automated helium mapping of individual zircon crystals, without the human error of manual re-entry. The design even fails safely: valves default closed, and an isolation valve on the scroll pump inlet snaps shut during power loss, preventing a rush of air into the turbo pump, while an inexpensive uninterruptible power supply and a backup compressed-air reservoir keep the line running through most outages.

Perhaps the most scientifically interesting part of the project is its calibration from first principles. Mass spectrometers measure ion-beam sizes, not moles, so converting signals into absolute helium amounts requires knowing exactly how much gas each spike aliquot delivers. The team built three ten-liter reservoir tanks feeding roughly 0.2-cubic-centimeter pipettes, and determined their volumes through a chain of metrological reasoning. First, the volumes of ultra-pure aluminum rods and silicon cubes were computed from repeated weighings on a calibrated balance and literature densities known to parts per million. Those standard objects were then used in gas-expansion experiments, monitored by an absolute capacitance manometer and corrected for thermal transpiration, to calibrate a known reference volume. That known volume, in turn, calibrated the pipettes and tanks via Boyle’s law, with uncertainties propagated at every step. Because each helium-3 shot slightly depletes its reservoir, the team also tracks tank depletion directly with a dedicated correction tank, preventing small errors from compounding over the roughly five thousand measurements a typical year brings.

The payoff is already evident. The Jimbochron now routinely measures helium amounts below 0.1 femtomoles accurately and precisely, and analyses of secondary standards such as Durango fluorapatite and Juina zircon confirm its performance across the 0.1-to-1000-femtomole range typical of laser-ablation work. Full laser-ablation (U-Th)/He dates agree with expected values, and automated sequences have produced helium abundance maps within single zircon crystals. That sensitivity translates into smaller ablation pits, finer spatial resolution, smaller grains, and higher-precision dates on younger samples. Remarkably, the total component cost came in under 160,000 US dollars, a fraction of what closed commercial systems can cost, and the open documentation, from parts suppliers to archived code, offers a template other laboratories can follow. As emerging techniques like laser-ablation helium chronology mature, instruments built to be understood, modified, and shared may prove just as important as the discoveries they enable.

Subject of Research: Design, automation, and first-principles calibration of a low-volume helium measurement line for (U-Th)/He thermochronology

Article Title: Technical note: Design, construction, automation, and first-principles calibration of a low volume He measurement line

Article References: Metcalf, J. R., & Flowers, R. M. (2026). Technical note: Design, construction, automation, and first-principles calibration of a low volume He measurement line. Geochronology, 8(3), 447-462. https://doi.org/10.5194/gchron-8-447-2026

Image Credits: AI Generated

DOI: 10.5194/gchron-8-447-2026

Keywords: helium dating, U-Th/He chronology, thermochronology, mass spectrometry, laser ablation, vacuum technology, geochronology, isotope dilution, instrumentation, metrology, apatite, zircon

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 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

The project grew out of the Hominin Sites and Paleolakes Drilling Project, or HSPDP, an initiative launched in November 2008 with the goal of assembling continuous, high-resolution paleoenvironmental records from lake sediments across eastern Africa, close to the fossil and archaeological sites that document human evolutionary history. During the first ICDP-funded workshop, a pivotal moment came when Asfawossen Asrat and the late Mohammed Umer proposed the Chew Bahir basin as a new drilling location, even though it was not on the original list of candidate sites. Their case rested on intimate knowledge of the regional geology, which comprises thick fluvio-lacustrine deposits within a tectonic basin, and on unpublished industry gravity survey data from the Ethiopian Ministry of Mines indicating that the basin contained two to three kilometres of sediment. The site also offered a decisive scientific advantage: proximity to key archaeological localities such as Omo-Kibish, where some of the oldest known Homo sapiens fossils have been found, along with important sites in the Ethiopian Highlands.

Before any deep drilling could begin, the team needed proof that the archive was worth the considerable expense. In late 2009, researchers from the universities of Cologne, Addis Ababa and Aberystwyth recovered a pilot core, CB01, from the western margin of the basin using a rotary single-tube drill. The core extended nearly 19 metres into the lake deposits with a recovery rate of 81 percent, and an organic-rich layer sent for radiocarbon dating returned an age of 13,600 years before present at about 7.35 metres depth. That single date, included in the first full HSPDP proposal submitted in January 2010, demonstrated the site’s potential. A second campaign in November 2010 added five more short cores, CB02 through CB06, collected with a vibro-corer along a roughly 20-kilometre transect running from the western margin toward the geographic centre of the basin. Together, these six cores revealed fine-grained clays, silts and carbonate-rich horizons, along with diatoms, molluscs and fish bones that pointed to shifting lake stages through time.

The short cores also delivered a crucial quantitative insight: sedimentation rates in the basin vary dramatically with climate. During wet episodes, sediment accumulated at roughly 0.5 to 0.7 millimetres per year, while drier intervals such as the Last Glacial Maximum saw rates drop to about 0.1 millimetres per year, a six- to seven-fold difference. Rates also declined from the basin margin toward its centre. This spatial gradient in accumulation would later prove decisive for choosing where to drill. It also shaped an early scientific surprise: the team discovered that potassium in the sediments, which they had expected to indicate humid conditions and enhanced feldspar weathering, actually tracked aridity through a smectite-to-illite transformation driven by increasing pore-water alkalinity. Defending this reinterpretation of the potassium proxy against entrenched prior knowledge became one of the project’s formative intellectual battles.

The deep-drilling campaign itself was originally scheduled for 2013, but it was cancelled on short notice when the contracted drilling rig failed to meet contractual specifications. Forecasts of an impending ENSO and Indian Ocean Dipole extreme event, which threatened severe flooding in eastern Africa, prompted a further one-year postponement. To keep the team’s momentum and satisfy funding requirements, the researchers improvised a shorter campaign in March 2014, recovering a roughly 41-metre core, CHB14-1, from the basin centre. That operation came within days of disaster: intense overnight rainfall flooded the saline mudflat, rendering the fine lacustrine clays slippery and unstable, and the drilling crew had to camp at the basin margin for two days until the water receded. The episode taught the team a lesson they would never forget about the volatility of working on a playa surface, even in the supposedly dry season.

When the deep drilling finally commenced in November and December 2014, the team applied everything the earlier campaigns had taught them. Two parallel boreholes, HSPDP-CHB14-2A and 2B, were drilled about 20 metres apart, roughly halfway along the original transect and close to the western margin of the basin. The location was chosen using seismic lines provided by the oil company Tullow Oil, which helped identify undisturbed deep sediment layers free of tectonic displacement. Drilling proceeded around the clock in day and night shifts, using a wireline diamond coring rig fitted with specialised soft-sediment tools adapted from ocean drilling: a hydraulic piston corer for the softest sediments, an extended nose corer with a non-rotating shoe for firmer material, and an alien corer functioning like a standard rock-coring tool for the most lithified intervals. Each tool could collect three-metre cores of 6.4-centimetre diameter in plastic liners, and all three shared the same bottom hole assembly, allowing rapid switching as lithology changed with depth.

The duplicate coring strategy, which the team had to defend during planning because of its higher cost, paid off handsomely. Incomplete recovery is inevitable in unconsolidated sediments, with gaps especially common between coring runs, and the gaps in hole 2A, mostly associated with sandy layers, could largely be filled from the overlapping sections of hole 2B. Drilling in the second hole was continuously optimised based on what the first hole revealed, with the hydraulic piston corer deployed at depths where sandy deposits had caused core loss, recovering them in short, slow drives rather than rotation-based drilling. The combined recovery rate reached approximately 90 percent, and the two holes yielded cores of about 279 and 266 metres that were merged into a composite profile of roughly 293 metres. Logistics were formidable: drilling fluid demand, initially estimated at under 20,000 litres per day, more than doubled due to lost circulation into the formation, requiring two 20,000-litre water trucks to shuttle continuously from the Weyto River some 35 kilometres away, with the water stored in six hand-dug sumps at the site.

Merging the twin cores into a single continuous record proved to be the project’s most technically demanding challenge. The cores were shipped to the Continental Scientific Drilling Facility at the University of Minnesota, where they were logged, split, photographed and described section by section, with 306 sections in hole 2A and 276 in hole 2B. Correlation between the holes relied on identifying matching stratigraphic features, such as colour changes and magnetic susceptibility peaks, to define tie points where the record could switch from one core to the other. The splice went through three major versions as new datasets became available, from low-resolution whole-core logging through high-resolution scanning to X-ray fluorescence elemental data. Each revision shifted the composite depths, causing roughly ten percent of the approximately 14,000 discrete samples to fall off-splice and requiring manual re-attribution of 1,362 samples to their equivalent depths in the parallel core. In the final version, about five percent of tie points carried depth uncertainties of 10 to 20 centimetres, tracked through a colour-coded quality flag system. An automated alternative using a dynamic time warping algorithm was tested but abandoned because it could not exclude known data artifacts that the human analysts could recognise and disregard.

Dating the composite record demanded an equally multi-pronged approach. Radiocarbon dating of ostracods and bulk organic material constrained the upper sections, optically stimulated luminescence dating of fine-silt quartz grains extended the chronology deeper, and single-crystal argon-40/argon-39 dating of feldspars from volcanic ash layers anchored the oldest intervals. A tephra layer in the core was geochemically matched to the Silver Tuff outcropping at Konso, and an independent check came from identifying ash from the Corbetti Ignimbrite, dated to 177,000 years, at a depth corresponding to 178,000 years in the core, a striking confirmation of the age model’s accuracy. Confidence intervals grow with depth, from under 10,000 years in the uppermost 50 metres to 10,000 to 35,000 years across the oldest 100,000 to 620,000 years of the record. The resulting archive, documented across more than 25 publications, has revealed how Pleistocene climate variability in eastern Africa influenced hominin evolution, including hydroclimate shifts that may have shaped early human dispersal.

The lessons the team draws are both practical and philosophical. Drilling close to the basin margin minimised logistical risks while maximising sedimentation rates and time resolution, and the pre-study transect of short cores proved essential for testing hypotheses about proxies and chronometers before committing to expensive deep drilling. Duplicate coring, now regarded as best practice for scientific drilling in modern sedimentary basins, delivered higher recovery, better tool selection and additional material for dating, but the compositing process introduced its own hidden uncertainties with every tie point chosen and every splice revision. The team recommends keeping tie points to an absolute minimum, clustering splice revisions closely in time, and exploring hybrid approaches that combine automated pattern recognition with manual oversight. They also caution that small drill rigs and lightweight supply chains impose real limits, as the boreholes fell short of the original 400-metre target due to borehole instability, strong groundwater inflow and logistical constraints, factors that future projects in the Chew Bahir and Turkana basins will need to address with larger rigs and more detailed engineering plans. The paper is dedicated to Andrew Cohen and Mohammed Umer, two of the project’s founding visionaries, whose insistence on drilling twin cores rather than one ultimately transformed how paleoclimate archives of human evolution are built.

Subject of Research: Scientific deep drilling of overlapping sediment cores in the Chew Bahir basin, southern Ethiopia, to reconstruct paleoclimate during human evolution

Article Title: Scientific deep drilling in the Chew Bahir basin: advantages and pitfalls of two overlapping sediment cores

Article References: Foerster, V., Asrat, A., Bronk Ramsey, C., Brown, E. T., Deino, A., Erbello, A., Fischer, M. L., Gebregiorgis, D., Junginger, A., Kaboth-Bahr, S., Lane, C. S., Opitz, S., Noren, A., Roberts, H. M., Tiedemann, R., Vidal, C.-M., Viehberg, F., Vogelsang, R., Zachow, C., … Trauth, M. H. (2026). Scientific deep drilling in the Chew Bahir basin: advantages and pitfalls of two overlapping sediment cores. Scientific Drilling, 35(1), 61-81. https://doi.org/10.5194/sd-35-61-2026

Image Credits: AI Generated

DOI: 10.5194/sd-35-61-2026

Keywords: Chew Bahir, scientific drilling, paleoclimate, Homo sapiens, Ethiopian Rift, sediment cores, HSPDP, ICDP, duplicate coring, geochronology, tephrochronology, human evolution