The rock in question is an eclogite, a strikingly beautiful red-and-green rock that forms when basaltic oceanic crust is buried to great depths and recrystallizes under immense pressure. The red comes from garnet, the green from the sodium-rich pyroxene omphacite, and together they signal that the rock has experienced pressures far beyond those at which ordinary crustal minerals are stable. Eclogites are therefore treated by geologists as messengers from the deep, carrying information about conditions in the uppermost mantle back to the surface when subducted material is exhumed. The specimen analyzed in the new study was collected from an outcrop in Taleba Bay on Goodenough Island, in the Milne Bay province of Papua New Guinea, a region where the collision of the Australian and Pacific plates has created one of the most tectonically energetic landscapes on Earth.

To interrogate the rock, the researchers deployed a combination of cutting-edge analytical techniques capable of resolving pressure and temperature conditions at very fine scales. The decisive evidence came from tiny inclusions of coesite, a high-pressure polymorph of silica, trapped inside a host mineral during the rock’s formation. Coesite is a diagnostic mineral for ultrahigh-pressure metamorphism: it only forms at pressures equivalent to burial depths of roughly 90 kilometers or more, and its preservation as minute inclusions requires that the host crystal acted as a sealed pressure vessel, protecting the coesite from transforming back to quartz as the rock rose toward the surface. Finding coesite therefore allowed the team to conclude with confidence that their sample had been subducted to at least 90 kilometers below the surface, where it experienced temperatures of around 800 degrees Celsius.

The truly unexpected result emerged when the researchers compared conditions at different stages of the rock’s burial history. By reconstructing the pressure and temperature pairs recorded at shallower and deeper levels of the subduction path, they found that temperatures at around 45 kilometers depth were only about 100 degrees Celsius cooler than those at the maximum depth of more than 90 kilometers. In other words, the rock passed through a zone of anomalously high temperature at intermediate depth before entering a regime with a much colder geothermal gradient. Because temperature normally increases steadily with depth in a subducting slab, a gradient that steepens only modestly over nearly 50 kilometers of additional burial implies a dramatic thermal transition, one that the team describes as a change from very hot subduction at shallow levels to the more common colder subduction at depth.

Dr Jan Schönig of the Göttingen University Geoscience Center, who led the study, explained the significance of the observation: finding this change from hot to cold subduction was unexpected, he noted, because it contradicts current thermal models of subduction zones. Standard numerical models generally depict subduction zones as cold features overall, because the downgoing slab is much cooler than the surrounding mantle and acts as a thermal sink, chilling the interface region along which rocks descend. Hot conditions are expected to develop gradually with depth as the slab warms and as mantle material flows around it, not to appear early in the burial history and then give way to colder conditions. The Papua New Guinea eclogite records precisely the opposite pattern, forcing geodynamicists to reconsider how heat is distributed within the upper reaches of subduction channels.

Rocks recording higher-than-expected temperatures have been documented before, and anomalous thermal gradients have been inferred from other ultrahigh-pressure localities around the world. What distinguishes the new study is the clarity with which it captures a complete journey from shallow to deep levels, accompanied by an unambiguous shift to a relatively colder gradient. Previous reports could often be explained by later heating events or by uncertainties in the pressure-temperature reconstructions. Here, the coesite inclusions and the associated mineral assemblages pin down the thermal evolution along a continuous path, making this the first time such a hot-to-cold transition during a single subduction episode has been demonstrated so convincingly in a single rock.

The researchers consider two principal explanations for the anomaly, which may also operate in combination. The first is that rocks at shallow subduction depths are heated far more intensely than previously assumed, potentially because of shearing stress along the plate contact. As the downgoing plate grinds against the overriding plate, frictional and mechanical work could generate localized heat in the subduction channel, elevating temperatures at depths of a few tens of kilometers well above what conductive models predict. The second possibility is temporal: the subduction zone in question may simply not yet have had time to cool and evolve toward the lower temperatures expected of a mature system. Young subduction zones, in which cold oceanic material has only recently begun to descend, may transiently record hotter conditions before steady-state thermal structure is established. A combination of shear heating and youthful subduction could plausibly account for the full pattern preserved in the eclogite.

Understanding which of these mechanisms dominates matters far beyond the academic satisfaction of refining a thermal model. Subduction zones are the sites of the planet’s largest earthquakes and much of its explosive volcanic activity, and both hazards are intimately controlled by temperature. Earthquake rupture depends on the frictional properties of the plate interface, which vary with thermal structure, while the melting processes that feed arc volcanoes are governed by the flux of water and heat in the subducting slab and the mantle wedge above it. A subduction zone that runs hotter at shallow depths than models assume could shift the boundaries of seismogenic behavior, altering where large earthquakes nucleate and how deep they extend. Accurate thermal models are therefore a prerequisite for meaningful seismic hazard assessment in convergent margins.

The implications extend to the chemistry of the Earth system as a whole. Subduction zones act as the planet’s long-term recycling valve, returning carbon dioxide and other volatile compounds locked in oceanic crust and sediments to the mantle, a process that has helped stabilize Earth’s climate over millions of years. The temperatures rocks experience during descent control how much carbon and other elements are released from the slab into the overlying mantle and atmosphere through arc volcanism, and how much is carried deeper into the mantle for long-term storage. If shallow levels of young or rapidly shearing subduction zones are systematically hotter than assumed, the fluxes of carbon and other geochemically important species through these systems could differ substantially from current estimates, with consequences for reconstructions of the long-term carbon cycle across geological time.

The study also speaks to one of the deepest questions in geology: how subduction has operated over millions to billions of years and how it continues to shape the planet today. Plate tectonics is the defining feature of Earth among the rocky planets of the solar system, and the thermal regime of subduction zones is central to how the mantle cools, how continents grow and recycle, and how the surface and interior exchange material. Every well-constrained pressure-temperature path from a natural laboratory such as Goodenough Island provides a ground-truth check on the numerical models that extrapolate these processes through deep time. The research, published in Nature Geoscience as a study led by Schönig and colleagues documenting a geothermal gradient change recorded by ultrahigh-pressure eclogite, was supported by funding from the German Research Foundation. As more young ultrahigh-pressure rocks are identified and analyzed with increasingly sensitive techniques, geologists may find that the hot-to-cold transition recorded in this four-million-year-old traveler is not an oddity, but a common and previously overlooked chapter in the life of every young subduction zone.

Subject of Research: Thermal evolution and geothermal gradient changes in subduction zones recorded by ultrahigh-pressure eclogite from Papua New Guinea

Article Title: Two kinds of subduction recorded in one rock

Article References: Two kinds of subduction recorded in one rock. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: subduction zones, eclogite, ultrahigh-pressure metamorphism, coesite, geothermal gradient, Papua New Guinea, plate tectonics, shear heating, carbon cycle, earthquakes, volcanism, Nature Geoscience

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

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

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

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

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

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

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

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

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

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

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

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

Image Credits: AI Generated

DOI: 10.5194/gchron-8-423-2026

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

Ultrahigh-pressure metamorphism is the process by which crustal rocks, the buoyant material of continents, are dragged down to depths greater than 90 kilometres, deep enough for the mineral coesite, a dense high-pressure form of quartz, to become stable. Because reaching those depths typically requires rapid, efficient transport into the mantle, geologists have long treated the presence of coesite as evidence for cold subduction, where the descending slab stays cooler than its surroundings. Peak metamorphic conditions recorded by such rocks usually imply low temperature-to-pressure gradients of less than 330 degrees Celsius per gigapascal, consistent with thermal models of mature subduction zones. Yet, as the new study points out, direct mineralogical evidence for a cold prograde journey, the actual path a rock takes on its way down, has been conspicuously lacking. Overprinting during burial and exhumation tends to erase the early record, so the cold prograde gradient has always been an assumption inferred from peak conditions rather than an observation.

The rock at the centre of the story, sample 03118b, was collected from eastern Taleba Bay on southern Goodenough Island, close to a beach placer where coesite-bearing detrital garnets had previously been discovered. It forms a small lens of mafic composition, roughly 90 by 30 centimetres, hosted within felsic gneiss. The eclogite is partially retrogressed, meaning it was partially altered as it returned toward the surface, but its garnet crystals remain pristine archives of the burial history. Garnet is a remarkable mineral for this purpose because it grows in concentric zones, each layer locking in whatever minerals, melts and fluids were present at the moment of crystallization. In this sample, the garnet displays three distinct compositional zones: a calcium-rich core, a mantle zone with lower calcium, and a rim with rising calcium again, each zone corresponding to a different stage of the subduction journey.

The inclusion inventory within the garnet zones provides the crucial mineralogical evidence. The cores contain quartz, rutile, zircon, apatite, omphacite and titanite, but no coesite. The mantle zones, by contrast, host monomineralic coesite inclusions measuring roughly one to four micrometres, along with slightly larger bimineralic coesite-quartz inclusions of about five micrometres. Because coesite is only stable above the quartz-to-coesite transition at approximately 2.7 to 3.0 gigapascals, the boundary between core and mantle marks the moment the rock crossed into the ultrahigh-pressure stability field. The size difference between the two inclusion types also supports the idea that larger inclusions preserve coesite better, since the surrounding garnet exerts protective overpressure on its trapped cargo.

Even more revealing are the polyphase inclusions found in both core and mantle. These tiny droplets, termed nanogranitoids, consist of silica polymorphs such as quartz or cristobalite, feldspar polymorphs including kumdykolite, kokchetavite and dmisteinbergite, and phyllosilicates such as muscovite, phlogopite and chlorite, sometimes accompanied by magnetite and carbonaceous matter. Their peculiar metastable mineralogy is the signature of a silicate melt that was trapped as a liquid and then crystallized in place within the confined space of the inclusion. In other words, the garnet grew in the presence of melt, meaning the rock was above its solidus, the temperature at which partial melting begins, throughout its burial. That single observation carries enormous thermal implications, because melting of basic rocks requires temperatures exceeding roughly 670 degrees Celsius.

To pin down the exact conditions, the team deployed an impressive analytical arsenal. Raman spectroscopy mapped the inclusions in two and three dimensions at sub-micrometre resolution. Zirconium-in-rutile thermometry on rutile inclusions yielded temperatures from the Zr content of the tiny crystals. Most dramatically, the researchers performed re-homogenization experiments using a piston-cylinder apparatus at GFZ Potsdam: garnet chips containing crack-free nanogranitoids were compressed to 2.5 gigapascals and heated to 750 degrees Celsius for 24 hours, successfully melting the crystallized inclusions back to glass without any reaction with the host garnet. This experiment demonstrated that the garnet core grew at approximately 700 to 750 degrees Celsius. Combined with thermodynamic modelling of the garnet composition, the team arrived at a best estimate for core formation of 1.5 plus or minus 0.2 gigapascals and 700 plus or minus 20 degrees Celsius, while the mantle zone, constrained by coesite and rutile thermometry, formed at pressures greater than 2.9 gigapascals and about 800 plus or minus 25 degrees Celsius.

Dating the journey required extracting micrometre-sized zircon inclusions from the garnet, a feat accomplished by dissolving the host mineral and handpicking the roughly ten-micrometre zircons under a stereo microscope. Using a large-geometry secondary ion mass spectrometer at Curtin University, the team analysed 26 grains, of which 17 yielded reliable uranium-lead ages spanning 5.1 to 3.1 million years. The trace-element chemistry of the zircons allowed the dates to be tied to specific garnet growth zones. All grains show low thorium-to-uranium ratios below 0.05, no negative europium anomaly and flat heavy rare-earth-element patterns, indicating growth in equilibrium with garnet, which preferentially partitions heavy rare earths. From 5.1 to 4.4 million years ago, heavy rare-earth concentrations declined as garnet grew in abundance, marking the interval of subduction from pressures of at most 1.5 gigapascals to greater than 2.9 gigapascals.

That timing yields a staggering result. Covering the pressure interval from 1.5 to 2.9 gigapascals, equivalent to roughly 45 to more than 90 kilometres of depth, in only about 700,000 years implies vertical subduction rates exceeding 60 millimetres per year, assuming lithostatic pressure. The ultrahigh-pressure stage itself, dated at about 4.4 million years, is considerably younger than the roughly 7.1-million-year coesite eclogite previously reported from Tumagabuna Island nearby, confirming that the Woodlark Rift region hosts Earth’s youngest known coesite eclogites. The rapid burial also explains why the earlier prograde record survived at all: the rock was whisked down and back up before diffusion and overprinting could erase the evidence.

When the pressure-temperature data are assembled into a path, the thermal transition becomes unmistakable. The shallow prograde path, from the surface to about 1.5 gigapascals, follows an apparent gradient of roughly 470 degrees Celsius per gigapascal, far warmer than any mature subduction zone model predicts, even for the warmest modern examples such as Mexico and Cascadia. Then, beyond about 45 kilometres of depth, the path pivots sharply into a cold regime consistent with gradients below 330 degrees Celsius per gigapascal, matching the thermal structure of mature subduction zones. The team considered whether the core and mantle might record two separate metamorphic events, but the tight zircon age range of 5.1 to 3.1 million years supports a single continuous burial episode. Intriguingly, published pressure-temperature paths from other ultrahigh-pressure terranes such as the Sanbagawa belt of Japan and the Dabie-Sulu orogen of China show similar warm-then-cold shapes, suggesting this may be a common feature rather than a local curiosity.

Two mechanisms, not mutually exclusive, can explain the warm start. The first is shear heating along the plate interface: as the descending plate grinds against the overriding plate, frictional energy dissipates as heat, warming the interface rocks. Thermal models that include shear heating with an effective coefficient of friction of about 0.07 produce interface temperatures closer to the observed path, though still somewhat cooler, and the present-day convergence between the Pacific and Australian plates at roughly 110 millimetres per year could push temperatures higher still. The second mechanism is the inherent time-dependence of young subduction zones. Numerical models show that in the first several million years after subduction initiates, before a quasi-steady-state thermal structure establishes, rocks entering the trench follow markedly warmer and more variable paths. The observed path closely resembles model predictions for rocks entering a subduction zone about 5 to 12 million years after initiation, and the inferred subduction rate of more than 60 millimetres per year matches models for rocks entering about 8 to 13 million years after initiation. The authors caution that peak metamorphic conditions alone, which in this rock imply a gradient of only about 280 degrees Celsius per gigapascal, can be profoundly misleading: some high-temperature rocks may have formed in zones that later turned cold, and some low-gradient rocks may never have experienced a cold, steady-state subduction system at all. For geologists using metamorphic rocks to reconstruct Earth’s tectonic past, the garnets of Goodenough Island deliver a clear message: the journey matters as much as the destination.

Subject of Research: Geothermal gradient evolution during early subduction recorded by a young ultrahigh-pressure eclogite from Papua New Guinea

Article Title: Geothermal gradient change during subduction recorded by ultrahigh-pressure eclogite

Article References: Schönig, J., Ferrero, S., Gresky, K., Sorger, D., Schmitt, A. K., Wunder, B., Baldwin, S. L., & Holt, A. F. (2026). Geothermal gradient change during subduction recorded by ultrahigh-pressure eclogite. Nature Geoscience. https://doi.org/10.1038/s41561-026-02110-1

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02110-1

Keywords: subduction zones, ultrahigh-pressure metamorphism, eclogite, coesite, garnet, geothermal gradient, zircon U-Pb dating, melt inclusions, shear heating, Papua New Guinea, thermobarometry, plate tectonics

The research, led by Samuel O. Sedara of the CSIR-National Geophysical Research Institute in Hyderabad and Adekunle Ajasin University in Nigeria, together with colleagues from Nigerian institutions, integrated an unusually broad set of data. The team compiled 175 heat-flow and geothermal gradient records alongside 1,300 radiometric measurements of the three heat-producing isotopes: uranium-238, potassium-40 and thorium-232. These measurements were organized across six geomorphic regions of the country, covering the northwest, northeast, north-central, southeast, south-south and southwest zones, allowing the authors to compare thermal behavior between Nigeria’s ancient crystalline basement complexes and its younger sedimentary basins for the first time in a single analytical framework.

What makes the study methodologically distinctive is its reliance on an independent multi-criteria decision analysis, or MCDA, framework designed to strip subjectivity out of geothermal favorability mapping. Rather than assigning weights to different geophysical parameters by expert judgment alone, the researchers derived weights objectively using Shannon entropy, a statistical measure of the information content of each variable, and cross-checked them with the analytic hierarchy process, which achieved a consistency ratio of just 0.003, far below the 0.10 threshold generally accepted as reliable. From this foundation the team constructed three composite indices: a thermal index, a heat-flow-weighted index, and a three-parameter geothermal favorability threshold index, each designed to explicitly separate conductive heat flow, the heat arriving from deeper sources, from radiogenic heat flow generated within the crust itself.

That separation is the analytical heart of the paper. Continental heat flow has two principal contributors: the slow decay of uranium, thorium and potassium in crustal rocks, and heat conducted upward from the mantle and lower lithosphere. By computing the ratio of conductive to radiogenic heat flow for each region, the authors could ask a deceptively simple question: is a region hot because its rocks are radioactive, or because something beneath the crust is supplying additional heat? In the northwest, the answer was unambiguous. The region recorded the dataset’s highest conductive heat flow at 151.4 milliwatts per square meter and a geothermal gradient of 56.0 degrees Celsius per kilometer, and its conductive-to-radiogenic ratio of roughly 86 means that 98.8 percent of the surface heat flow cannot be attributed to crustal radioactivity.

The implications of that number are profound. A purely radiogenic explanation for the northwest’s thermal anomaly is inconsistent with the data, the authors conclude, and a deeper heat source is required. While surface measurements and Curie point depth estimates, which track the depth at which magnetic minerals lose their magnetization and provide an independent proxy for crustal temperature, cannot distinguish among mantle upwelling, lithospheric thinning, or deep tectonic reactivation as the specific mechanism, the decoupling between conductive and radiogenic heat is real and quantified. This is, according to the study, the primary geoscientific result: Nigeria’s northwest is not merely warm, it is anomalously hot for reasons rooted in the deep Earth, placing it in a geodynamic category that few expected for a region far from active plate boundaries.

Crucially, the finding does not rest on a single line of evidence. The researchers subjected their rankings to a battery of independent validations. Bootstrap resampling with 1,000 iterations tested the stability of the favorability class boundaries at plus or minus 0.50 standard deviations. A random forest regression model, a machine learning algorithm that builds predictions from many decision trees, achieved a coefficient of determination of 0.702, a moderate but respectable performance given the small regional sample. K-means unsupervised clustering, which groups regions by similarity without reference to the MCDA labels, produced a silhouette score of 0.397 and an adjusted Rand index of 0.074, modest values the authors attribute candidly to the limited number of data points, since with only six regions even a single misclassification sharply depresses the agreement metric. The team is explicit that these machine learning results should be read as supporting evidence within the overall multi-method framework rather than as standalone proof.

The convergence, however, is what gives the northwest result its weight. Curie point depth-derived heat flow, random forest regression, bootstrap validation and prior independent studies all point to the same region. Across every composite index, the northwest ranks first, with a thermal index of plus 1.715, a heat-flow-weighted index of plus 1.864 and a geothermal favorability threshold index of plus 1.343. The Sokoto Basin, which spans Sokoto, Kebbi and Zamfara states, together with nearby northwest sedimentary terrains, is continuously identified as the most attractive exploration target in the country. The authors recommend urgent geophysical investigation of the basin and slim-hole drilling, along with a dedicated thermal model to establish actual temperature-at-depth profiles before any full exploration campaign begins.

The study also redraws the map of geothermal opportunity in a second, very different way. The southeast region, with a thermal index of plus 0.490, represents not a secondary version of the northwest prospect but an entirely distinct type of geothermal target. Its thermal regime is radiogenically sourced and dominated by potassium-40, meaning its heat comes from the decay of radioactive elements within crustal rocks rather than from below. That character makes the southeast better suited to direct-use applications, such as district heating, agricultural drying and greenhouse warming, than to power generation, and the authors caution that it cannot be directly compared with the northwest. The north-central and northeast regions emerge as moderately favorable, while the south-south and southwest rank as unfavorable in this screening.

The authors are careful to frame the scope of their conclusions. This was a regional-scale screening exercise, and it did not address the factors that ultimately determine whether a geothermal prospect can be developed: permeability, fluid content, reservoir depth, sustainable extraction rates and drilling feasibility. Those questions remain prerequisites for any actual development decision, and no amount of surface and magnetic data can substitute for them. Yet the study transforms the starting point for Nigerian geothermal exploration. By quantifying, for the first time at national scale, where the country’s heat comes from and how strongly each region scores across independently weighted and validated indices, it converts a diffuse sense of possibility into a ranked, evidence-based agenda. If subsequent geophysics and drilling confirm what the numbers suggest, the Sokoto Basin could place Nigeria on the geothermal map, adding a firm, weather-independent renewable resource to an energy economy in urgent need of diversification.

Subject of Research: National-scale geothermal resource assessment of Nigeria using multivariate and machine learning analysis of heat flow and radiometric data

Article Title: A New Multivariate Analysis of Geothermal Assessment of Nigeria’s Crystalline Basement and Sedimentary Basins: Thermal and Geodynamic Implications

Article References: Sedara, S. O., Alabi, O. O., Adekanle, O. J., Adekanye, O. O., & Christopher, O. A. (2026). A New Multivariate Analysis of Geothermal Assessment of Nigeria’s Crystalline Basement and Sedimentary Basins: Thermal and Geodynamic Implications. Natural Resources Research. https://doi.org/10.1007/s11053-026-10773-1

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10773-1

Keywords: geothermal energy, Nigeria, heat flow, radiogenic heat production, Sokoto Basin, multi-criteria decision analysis, random forest, Curie point depth, geothermal gradient, crystalline basement, sedimentary basins, geodynamics