The award citation, written by Thomas Pettke of the University of Bern, traces a career that began in 2012, when Eberhard started her undergraduate studies in Bern. Her bachelor thesis, completed in 2015, dealt with the crystal-chemical characterisation of the barite-celestine solid-solution series, a system in which barium and strontium sulfate minerals substitute for one another in a continuous crystal lattice, a classic problem in mineral chemistry that requires precise control of composition and structure. For her master’s thesis she turned to a question that would shape much of her later work: the petrochemical fate of ophicarbonate, a rock mixture of serpentinite and carbonate, as it is progressively heated in the contact aureole of the Bergell Pluton in Valmalenco, Italy. Contact aureoles are the metamorphic halos that surround igneous intrusions, and they provide natural laboratories in which the thermal transformation of hydrated and carbonated rocks can be studied directly in the field. Her master’s work earned the University of Bern Faculty of Science Award for the best MSc thesis in the geosciences in 2017.

From Bern, Eberhard moved to the Bavarian Research Institute of Experimental Geochemistry and Geophysics in Bayreuth, where she carried out her doctoral research under the supervision of Dan Frost. Pettke’s citation notes that she arrived without prior experimental experience, yet went on to devise a series of multi-anvil experiments, the high-pressure apparatus in which a sample is squeezed between opposing tungsten carbide anvils to reproduce conditions deep within the Earth. Her experiments were designed to derive thermodynamic parameters for the ferric iron and aluminium endmembers of the serpentine minerals antigorite and lizardite, thereby quantifying the relationship between the oxidation state of iron in serpentinite and the oxygen fugacity of the system. Oxygen fugacity, the effective partial pressure of oxygen, is a master variable in geochemistry because it controls which elements change valence and therefore how much redox budget a rock can carry or release.

A second line of doctoral work quantified the effect of carbon-bearing phases on fluid release in subducted serpentinite. Serpentinites form when seawater alters mantle rocks at the ocean floor, binding enormous quantities of water and, where carbonation occurs, carbon dioxide into their crystal structures. As these rocks are dragged down subduction zones, heating drives dehydration and decarbonation reactions that release fluids into the overlying mantle wedge, triggering the melting that feeds arc volcanoes. A third, novel experimental approach allowed Eberhard to constrain how permeability evolves in natural serpentinite as it dehydrates at high pressures and temperatures, addressing one of the central uncertainties in subduction geophysics: how quickly and along which pathways fluids can escape from the descending slab. Taken together, Pettke writes, her findings offer prominent and quantitative insights into redox budget cycling from the subducting slab to the sources of arc magmas.

After completing her doctorate in 2020, Eberhard pursued postdoctoral research, first in Bayreuth and then in Utrecht in the group of Oliver Plümper. There she diversified her research into the complex interplay between rock deformation, prograde mineral reactions and fluid migration, working both as an independent lead scientist and as an active collaborator. Prograde reactions are those that proceed as temperature and pressure rise during burial, and their rates and pathways are strongly modulated by deformation, which creates new surface area, opens fluid pathways and can drive reaction overprint. The results of this quantitative work, according to the citation, already prominently advance the understanding of how deformation affects petrology and, conversely, how petrochemical processes influence rock mechanics, a coupling that remains one of the hardest problems in solid-Earth science because it spans scales from atomic diffusion to plate tectonics.

The Foundation’s citation also highlights Eberhard’s service to the scientific community. She was instrumental in revitalising scientific exchange in the serpentinite research community after the COVID-19 pandemic, launching and running the monthly webinar series SerpentineDays together with colleagues, and co-organising the SerpentineDays international conference in Granada, Spain, in 2024. Pettke describes her as a dedicated teacher and mentor who offers inspiration, sound knowledge and hands-on support to the next student generation, and credits her catching enthusiasm, scientific curiosity, rigour and independent, highly creative mindset as forces already shaping innovative research directions.

In her response, Eberhard expressed deep gratitude to the Foundation and to Pettke, who nominated her and supervised her master’s work a decade ago. She recalled how he introduced her to metamorphic petrology and the geochemistry of carbonated serpentinites in Val Malenco during two field campaigns, and how he encouraged her ideas with the Swiss German phrase Mach mou, meaning go for it, whenever she shyly asked whether a question was worth pursuing. She noted that he supported her long after she left Switzerland, and joked that in German the terms Doktorvater and Doktormutter for doctoral advisers should be extended in her case to include Mastervater as well.

She thanked Dan Frost for introducing her to experimental petrology and for helping her continue the serpentine-carbonate phase relations work that grew out of her master’s thesis, admitting that experimental work proved far harder than she had naively imagined and crediting his persistent questions with teaching her to evaluate her own results critically. At Bayreuth she also collaborated with Catherine McCammon, a Mössbauer spectroscopy specialist whose methods are central to determining iron oxidation states in minerals, and with geodynamicist Gregor Golabek, collaborations she says shaped both her research and her career. Her subsequent move to Oliver Plümper’s group in Utrecht gave her the scientific freedom to develop new research questions and test them experimentally, and she has now returned to work with him at the University of Bremen.

Eberhard’s international training included six months at Tohoku University in Sendai, Japan, during her doctorate as part of a research training group funded by the German Research Foundation, and research stays at the Institute of Earth Sciences in Orléans, France, where she learned from Holger Stünitz how to develop, run and interpret deformation experiments. Science, she reflected, is certainly not a one-woman show, and without the help of so many people she could not have achieved so much. Among the lessons of eight years of experimental research, she singled out the fact that experiments often do not work out as planned, and that the knowledge gained from failed experiments is invaluable, provided one keeps the curiosity and creativity to try again with a different approach. She closed by looking forward to exploring new ideas and to becoming a good mentor herself, bringing the fascination of science to the next generation.

The Paul Niggli Medal, named for the pioneering Swiss mineralogist and petrologist after whom the foundation is also named, continues to serve its founding purpose of honouring and supporting young ambassadors of Swiss geoscience as they build international careers. Eberhard’s trajectory, from a bachelor thesis on solid-solution crystal chemistry in Bern to high-pressure experiments that constrain the redox and fluid budget of the world’s subduction zones, exemplifies the kind of quantitative, experimentally grounded petrology the medal was created to recognise, and her combination of laboratory rigour, field grounding and community leadership marks her as one of the most versatile early-career Earth scientists the Swiss system has produced in recent years.

Subject of Research: Award of the 2025 Paul Niggli Medal to experimental petrologist Lisa Eberhard for research on fluid and redox budget transfer in subduction zones

Article Title: Lisa Eberhard receives the 2025 Paul Niggli Medal

Article References: Schönbächler, M., Pettke, T., & Eberhard, L. (2026). Lisa Eberhard receives the 2025 Paul Niggli Medal. Swiss Journal of Geosciences, 119, 61-63. https://doi.org/10.3897/sjg.119.183296

Image Credits: AI Generated

DOI: 10.3897/sjg.119.183296

Keywords: Paul Niggli Medal, Lisa Eberhard, subduction zones, serpentinite, experimental petrology, geochemistry, redox budget, fluid release, multi-anvil experiments, Swiss Journal of Geosciences, University of Bremen, metamorphic petrology

The key player is bridgmanite, the magnesium-iron silicate that makes up roughly 80 percent of the lower mantle in the standard pyrolite model. Like many minerals, bridgmanite crystals are intrinsically anisotropic: seismic waves travel faster along some crystallographic directions than others. When rocks deform slowly over millions of years, dislocation creep causes grains to rotate and align into a crystallographic preferred orientation, or CPO, and the aggregate as a whole becomes anisotropic. In the average lower mantle, diffusion creep is often cited as the dominant deformation mechanism, which would scramble any fabric. But numerical models and seismic observations agree that in high-stress regions such as subduction zones, dislocation creep remains active, allowing fabrics to persist and grow to great depths.

What makes the new work possible is a recent set of laser-heated diamond anvil cell experiments that deformed pyrolitic material at true lower-mantle pressures and temperatures. Those experiments revealed something unexpected: bridgmanite does not deform the same way at all depths. Below about 24 gigapascals, the easiest slip system is (100)[001], but as pressure climbs toward 50 gigapascals, corresponding to roughly 1,250 kilometers depth, the critically resolved shear stresses of the (100)[001] and (010)[100] systems cross at about 40 gigapascals, near 1,050 kilometers. Past that crossover, (010)[100] takes over as the dominant slip system. The transition is gradual, spread across the depth range from about 670 to 1,250 kilometers, yet its consequences for seismic imaging turn out to be sharp.

To translate this mineral physics into a seismic prediction, the team coupled a two-dimensional thermo-mechanical model of intraoceanic subduction, built with the open-source geodynamics code ASPECT, with viscoplastic self-consistent fabric calculations that track the orientations of 3,000-grain aggregates advected through the simulated mantle flow. They modeled a slab that flattens near 1,000 kilometers depth, the so-called deeply flattened configuration seen in regions such as the Western Java subduction zone, and compared it with a shallowly stagnating slab trapped near the 660-kilometer boundary. Bridgmanite was treated as the only phase that develops meaningful texture, since davemaoite and ferropericlase showed no coherent fabrics in the experiments and were assigned random orientations.

The result is a deformation-induced anisotropic discontinuity. As mantle material is entrained around the slab and accumulates strain, the pressure-dependent slip transition produces a layer where the anisotropy index, a scalar measure of the total anisotropy in the elastic tensor, shows a persistent local minimum. In the mantle wedge and subslab regions, this discontinuity emerges near 1,000 kilometers and progressively deepens toward the slab, reaching about 1,400 kilometers directly beneath it. Remarkably, even though the underlying change in slip systems is smooth, the resulting seismic boundary appears sharp: at locations where natural strain exceeds 2, the half-width of the discontinuity is only about 40 to 54 kilometers, compared with roughly 80 kilometers where strain is weaker. Where the slab stagnates shallowly and strains stay low, the discontinuity persists near 1,000 kilometers but becomes markedly fainter.

The crucial question is whether such a subtle feature could actually be seen by seismometers. The team computed reflection coefficients for underside reflections, the SS precursor phases that seismologists use to probe mantle discontinuities, at three locations along the modeled boundary. Because the medium is anisotropic, the reflection coefficients vary with both azimuth and epicentral distance, and their polarities flip depending on the local flow pattern: negative for all azimuths where poloidal flow dominates, predominantly positive where slab-driven flattening stretches entrained aggregates. At the most strongly strained location, reflection amplitudes reach about 0.02, comparable to coefficients computed for textural interfaces invoked to explain deep-mantle reflections and to those generated by the well-known transition-zone discontinuities. P-wave reflections, by contrast, were negligible, indicating that the boundary is essentially a shear-wave feature.

Full waveform simulations sealed the detectability argument. Using the spectral-element code AxiSEM3D, the researchers built a three-dimensional global model of the anisotropic mantle and recorded synthetic seismograms at a virtual array of 3,000 stations. Array processing with slowness vespagrams revealed a clear SS precursor generated by the anisotropic discontinuity alone, with no isotropic velocity contrast required. The precursor-to-SS amplitude ratio of 0.0084 exceeds the global detectability threshold of 0.0065 established in earlier surveys of mid-mantle reflectors. An isotropic model containing only the slab’s compositional velocity anomaly produced a much shallower and closer-in reflection from the slab-mantle interface, suggesting that where both a shallow and a deep reflector are observed, as in some global SS precursor studies, a flattened slab and its deformation halo could simultaneously supply both.

The pressure-dependent fabric also reconciles long-standing puzzles in mantle anisotropy. Previous forward models that imposed a single bridgmanite slip system could match some features of radial anisotropy, expressed as the ratio of horizontally to vertically polarized shear-wave speeds, but always failed somewhere: one fabric produced too little subslab anisotropy, another generated strong anisotropy everywhere including the mantle wedge where observations show the opposite sign. The new models, with slip systems shifting smoothly with depth, reproduce the fast horizontal shear beneath the slab and the faster vertical shear in the mantle wedge, matching tomographic images of the Western Java subduction zone. Modeled shear wave splitting also matches source-side observations beneath Western Java, with trench-parallel fast directions for downdip rays and trench-perpendicular fast directions for updip rays, although computed delay times remain somewhat low, plausibly because anisotropic phases such as akimotoite and phase D, shape-preferred orientation, and receiver-side complexity all add to the path-integrated signal.

The implications reach beyond subduction zones. The enigmatic 1,000-kilometer discontinuity has been variously attributed to a viscosity jump caused by grain-size variations in ancient bridgmanite-enriched domains, to stiff periclase-rich material, or to sluggish pyroxene-to-garnet transformations in cold slabs. The new mechanism does not replace these thermo-chemical hypotheses, the authors stress, but adds a purely deformation-driven alternative that requires no compositional contrast at all, only accumulated strain and the pressure sensitivity of bridgmanite’s crystal plasticity. Candidate testing grounds include the Cocos, Northern Peru, and Kamchatka slabs, which all show mid-mantle flattening, while vertically penetrating slabs such as the Marianas should instead produce narrow, slightly dipping anisotropic boundaries hugging the slab tip, an equally plausible explanation for the dipping reflectors imaged there. Looking ahead, the team proposes that fabric evolution in bridgmanite may even feed back on slab dynamics, softening the mantle around descending plates much as olivine fabric eases plate motion in the asthenosphere. If so, the deepest architecture of our planet may be written not only in chemistry and temperature, but in the strained, aligned grains of its most abundant mineral.

Subject of Research: Deformation-induced seismic anisotropy and mid-mantle discontinuities beneath subduction zones caused by depth-dependent bridgmanite deformation

Article Title: Mid-mantle discontinuities beneath subduction zones linked to depth-dependent deformation of bridgmanite

Article References: Magali, J. K., Merkel, S., Pisconti, A., Gay, J., Wookey, J., & Thomas, C. (2026). Mid-mantle discontinuities beneath subduction zones linked to depth-dependent deformation of bridgmanite. Nature Communications, 17(1), Article 10494. https://doi.org/10.1038/s41467-026-78212-6

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78212-6

Keywords: bridgmanite, subduction zones, mid-mantle discontinuity, seismic anisotropy, crystallographic preferred orientation, mantle dynamics, SS precursors, shear wave splitting, mineral physics, lower mantle, geodynamics, seismology

The debate at the heart of the study concerns the efficiency of deep carbon recycling in hot subduction systems. Hot slabs, which are young, warm and buoyant, were long expected to lose their carbonate cargo early. Contrasting models make opposite predictions: one holds that extensive decarbonation occurs at shallow depths, releasing carbon dioxide into the forearc and arc regions before the slab can carry carbon any deeper; the other holds that substantial carbonate survives within the subducted slab and is transported into the deep mantle. Which of these scenarios dominates matters enormously, because subduction zones are the principal route by which carbon returns from the surface reservoir to the mantle, balancing the carbon vented by volcanoes and mid-ocean ridges over geological time.

To test these competing pictures, the team examined olivine-hosted melt inclusions in two distinctive suites of volcanic rocks. The first are adakites from western Panama, erupted above a modern hot subduction system where young oceanic crust descends beneath Central America. The second are high-magnesium andesites from southwest Japan, which the researchers treat as a window into Archean-like hot subduction, the kind of tectonic regime thought to have operated on the early Earth more than two and a half billion years ago, when the mantle was hotter and slabs warmer. Melt inclusions are precious to geochemists because they seal a sample of the parental magma inside a host crystal as it grows, protecting it from later degassing and alteration that would otherwise erase the volatile record.

The key measurement is the ratio of boron’s two stable isotopes, expressed as the δ11B value. Boron has an unusually light isotope, boron-10, and a heavier one, boron-11, and the two are fractionated strongly at low temperatures near the seafloor, where seawater alteration loads the oceanic crust with boron enriched in the heavy isotope. When a slab is cold, this heavy boron is flushed out early by fluids, so arc magmas in cool subduction zones typically carry moderate δ11B signatures. In the melt inclusions from Panama and Japan, however, the team found strikingly elevated δ11B values reaching up to +23.2 per mil, among the heaviest recorded in arc magmas, indicating that the boron source had been heated to the point of melting rather than simply dewatering.

Crucially, the heavy boron isotope signal does not stand alone. The same melt inclusions show high ratios of cerium to lead, lanthanum to ytterbium, europium to titanium, and thorium to zircon. Each of these trace-element ratios is a diagnostic tracer of contributions from subducted sediment or carbonate-bearing lithologies that have been melted, rather than merely fluxed by aqueous fluids. Together, the elevated δ11B values and the enriched trace-element patterns indicate that the arc magmas incorporated melts derived from carbonate-bearing parts of the subducted slab. In other words, the slab was not shedding its carbon as fluid at shallow depth; it was carrying it down intact until temperatures grew high enough to melt the carbonate-bearing material directly.

The team then integrated these melt-inclusion observations with regional estimates of carbon dioxide flux from western Panama. The comparison yields a striking quantitative conclusion: more than roughly 60 percent of the carbonate in the subducted slab escapes shallow decarbonation. That figure challenges the intuition that hot slabs should decarbonate efficiently, and the authors propose a mechanism to explain the paradox. Extensive dehydration in the forearc and beneath the arc, they argue, causes such substantial fluid loss early in the subduction journey that the slab is left essentially dry by the time it reaches the deeper forearc and the arc front. Because fluid-mediated decarbonation depends on water-rich fluids infiltrating and reacting with carbonate, a dehydrated slab is a protected slab: with little fluid left to drive the reactions, the carbonate survives.

This mechanism reframes how scientists think about the fate of carbon in warm subduction zones. Rather than heat alone determining whether carbon is released, the balance between dehydration and decarbonation becomes decisive. A hot slab loses its water quickly and thoroughly, and in doing so it loses the very agent needed to strip away its carbonate. The result is that substantial carbon can bypass shallow decarbonation and continue its descent, in both modern hot subduction systems like Panama and Archean-like systems like the one recorded by the southwest Japan andesites. The finding extends the conclusion back into deep time, suggesting that early Earth’s hotter tectonics did not necessarily vent more carbon to the atmosphere; it may instead have buried more of it.

The implications ripple outward across the Earth sciences. Global carbon cycle models, which attempt to balance volcanic outgassing against subduction inputs over millions of years, depend on assumptions about how much carbon each slab delivers to the deep mantle. If more than half of slab carbonate survives shallow processing in hot subduction zones, the deep carbon flux may be considerably larger than many current estimates allow. That, in turn, bears on reconstructions of atmospheric carbon dioxide and climate across geological eras, including the Archean, when the interplay of a hotter mantle, faster plate recycling and a different biosphere shaped the planet’s early climate system. It also informs the study of carbonated mantle peridotites and other proposed hidden reservoirs of subducted carbon, which recent work has highlighted as potentially significant sinks.

Methodologically, the study demonstrates the power of melt inclusions as archives of slab inputs. Because the inclusions trap the magma before degassing, they preserve volatile and trace-element information that whole-rock samples, altered by eruption and weathering, often lose. Pairing boron isotopes with a suite of trace-element ratios gives the approach redundancy: the isotope signal identifies a melted, carbonate-bearing source, while the element ratios independently corroborate it. Applying the same toolkit to both a modern arc and an Archean analogue gives the conclusion a temporal reach that single-site studies cannot achieve, allowing the authors to argue that the behavior is a general property of hot subduction rather than a local peculiarity of Panama.

Questions remain, as they always do in deep-Earth science. The carbon budget of any single subduction zone is difficult to constrain, and regional flux estimates carry uncertainties that propagate into the 60 percent figure. Extending the melt-inclusion approach to other hot subduction systems around the world will be needed to test whether the Panama and Japan results are representative. But the core message of the study is clear and consequential: the hottest, most aggressive subduction zones may be among the most efficient at smuggling carbon into the deep mantle, not the least. In the planet’s long-term carbon ledger, the descending slab appears to be a far more faithful courier than many had assumed, carrying its carbonate cargo past the shallow gauntlet and into the depths where it may reside for hundreds of millions of years before returning, if ever, to the surface.

Subject of Research: Deep carbon recycling and carbonate retention in hot subduction zones traced by boron isotopes in melt inclusions

Article Title: Heavy boron isotopes reveal deep carbon recycling in hot subduction zones

Article References: Liu, H.-Q., Tian, F., He, M.-H., Cui, Z.-X., Zhang, L., He, P.-L., Hong, L.-B., Huang, X.-L., Bindeman, I., Hoernle, K., & Xu, Y.-G. (2026). Heavy boron isotopes reveal deep carbon recycling in hot subduction zones. Nature Communications. https://doi.org/10.1038/s41467-026-78096-6

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78096-6

Keywords: subduction zones, boron isotopes, carbon cycle, melt inclusions, adakites, decarbonation, deep carbon recycling, Archean, geochemistry, volcanic arcs, carbonate, Panama

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 research team, led by Christian Reinberg and Boriana Mihailova of the University of Hamburg, together with Stylianos Aspiotis of Hamburg University of Technology, Thomas Malcherek, and Stefan T. M. Peters, examined a natural fluorophlogopite crystal from the Cardiff Uranium Mine in Ontario, Canada. The specimen, a dark-brownish tabular crystal roughly 26 by 24 by 6 millimeters, came from the collection of the Museum of Nature Hamburg. It contained about 0.15 atoms per formula unit of divalent iron and approximately 0.40 atoms per formula unit of hydroxyl groups, making it an ideal candidate for probing how iron and water-bearing defects influence the thermal behavior of trioctahedral mica.

The experimental setup was elegant in its precision. The researchers cut thin rectangular flakes parallel to the crystal’s cleavage plane and heated them in air at 20 kelvin per minute on a calibrated heating stage, collecting Raman spectra with a green argon-ion laser focused to a spot about two micrometers across. Crucially, they extended their measurements down to 15 wavenumbers, capturing the lowest-energy vibrational mode of the crystal, a mode involving out-of-plane translations of the silicate layers that had never before been tracked at high temperature in this mineral. This low-frequency phonon turned out to be the key witness to a hidden structural drama.

Between roughly 550 and 650 kelvin, the team observed classic signatures of a phonon-driven structural instability. The wavenumber of the layer-translation mode near 92 wavenumbers showed a plateau-like minimum, while its linewidth peaked near 600 to 650 kelvin, a combination that typically heralds a phase transition driven by lattice vibrations. At the same time, the phonon expandability of this mode, a measure of how its frequency shifts with temperature, flipped from positive to negative, and the hard mode near 739 wavenumbers, known to be sensitive to the occupancy of the interlayer potassium site, began to rise anomalously with heating. Together, these signals pointed to a rearrangement of the interlayer region that mobilizes potassium cations.

This mobilization has profound geophysical implications. Because potassium ions carry positive charge, their delocalization within the interlayer space between 650 and roughly 1050 kelvin means they can act as charge carriers, contributing to the unusually high electrical conductivity that phlogopite exhibits above 900 kelvin. The process proved fully reversible up to about 1100 kelvin: when the sample was cooled, the spectral anomalies vanished and the crystal returned to its original state. This reversibility suggests that in the Earth’s crust and shallow mantle, where temperatures fall within this range, potassium migration through phlogopite could be an ongoing, dynamic process rather than a one-way chemical alteration.

Around 1150 kelvin, a second, more dramatic change occurred. The hydroxyl-stretching peaks in the Raman spectra, which arise from O-H bonds in the crystal, first intensified sharply, a sign that hydrogen cations were vibrating with greater amplitude just before breaking free from their oxygen anchors. Above that temperature, the peaks disappeared entirely, indicating that all hydrogen cations in the structure had delocalized. Remarkably, when the sample was cooled, the hydroxyl signals returned, showing that most of the hydrogen had not actually left the crystal but had merely been mobilized. Only after heating to 1450 kelvin did a permanent loss emerge, with about 65 percent of the hydroxyl groups recovering at room temperature, implying a 35 percent reduction attributable to dehydroxylation.

The third act of this thermal drama involved iron. Starting near 1300 kelvin, the divalent iron sitting in the octahedral sites of the crystal began to oxidize irreversibly to trivalent iron. The evidence came from a subtle but telling feature: a weak extra Raman signal at 602 wavenumbers that appeared only after the sample had cooled back to room temperature from 1450 kelvin. This signal is a resonance Raman scattering feature, produced when a polar optical phonon couples with an excited electron, a combination that breaks the normal symmetry selection rules of the centrosymmetric crystal. Its persistence after cooling confirmed that the iron oxidation was permanent. A minimum in the temperature trend of the tetrahedral ring-breathing mode at 1300 kelvin provided independent corroboration, reflecting the local puckering of silicate rings as they adjusted to the smaller iron-three-plus octahedra.

Above 1300 kelvin, the mineral began to fall apart, though only partially. Raman mapping of the heated crystal revealed tiny amounts of forsterite, the magnesium end-member of olivine, accounting for roughly one percent of the sample volume. Electron microprobe analysis and single-crystal X-ray diffraction confirmed the chemical picture: a slight shrinkage of the unit cell consistent with iron oxidation and hydroxyl loss, no detectable change in fluorine content, and a potassium loss of only about two percent. The confocal Raman maps delivered a surprise: the forsterite was concentrated in an ultrathin layer on the crystal surface parallel to the cleavage plane, decaying exponentially to negligible levels within about half a micrometer of the surface. The nanometer-scale grains were too small for X-ray diffraction to detect, explaining why earlier studies relying on diffraction alone had missed the spatial details of the decomposition.

The cumulative picture that emerges is one of a mineral that does not fail catastrophically at a single threshold but instead passes through a choreographed sequence of transformations, each with its own temperature window and its own consequences for charge transport. Between 500 and 650 kelvin, the bonds linking tetrahedral and octahedral sheets weaken, priming the interlayer for potassium mobilization. From 650 to 1050 kelvin, mobile potassium cations enhance conductivity along the cleavage plane. Near 1150 kelvin, hydrogen joins the ranks of mobile charge carriers. At 1300 kelvin, iron oxidation begins, and by 1400 to 1450 kelvin, dehydroxylation triggers partial structural collapse and the nucleation of forsterite. The authors note that the specific temperatures of each stage likely depend on the initial iron and hydroxyl contents of the sample and on the surrounding atmosphere, meaning that natural phlogopites with different chemistries may transform at somewhat different points.

For Earth scientists, the implications extend well beyond the laboratory. Phlogopite is a principal vehicle for transporting water and halogens into the upper mantle through subduction zones, reaching depths of roughly 180 to 200 kilometers before breaking down. The new results suggest that the mobility of charge carriers within this mineral begins at far lower temperatures than the established pressure-temperature stability fields would imply, meaning that mantle metasomatism, the chemical alteration of the mantle by fluids and melts, may operate over a wider temperature range and a larger depth profile than previously appreciated. Because fluorine substitution for hydroxyl groups extends the thermal stability of phlogopite, fluorine-rich varieties may carry these crystallochemical transformations to even greater depths, modulating the deep volatile cycle in ways that geophysical models are only beginning to capture. What happens to a single crystal on a heating stage in Hamburg, it turns out, echoes through the electrical structure of the deep Earth.

Subject of Research: Thermally activated structural and crystallochemical transformations in iron-bearing fluorophlogopite mica studied by in situ high-temperature Raman spectroscopy

Article Title: Thermally activated multistep alteration of Fe2+-bearing fluorophlogopite revealed by in situ Raman spectroscopy

Article References: Reinberg, C., Aspiotis, S., Malcherek, T., Peters, S. T. M., & Mihailova, B. (2026). Thermally activated multistep alteration of Fe 2+ -bearing fluorophlogopite revealed by in situ Raman spectroscopy. European Journal of Mineralogy, 38(4), 461-475. https://doi.org/10.5194/ejm-38-461-2026

Image Credits: AI Generated

DOI: 10.5194/ejm-38-461-2026

Keywords: phlogopite, fluorophlogopite, Raman spectroscopy, high temperature, potassium mobility, charge carriers, iron oxidation, dehydroxylation, forsterite, electrical conductivity, mantle metasomatism, subduction zones

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 physics underlying the approach centers on subduction zones, the convergent plate boundaries where one tectonic plate slides beneath another. These settings produce the planet’s largest earthquakes, including events exceeding magnitude 8.5, and frequently unleash devastating tsunamis when the seafloor lurches upward or sideways. But not every part of a subduction interface moves the same way. Some patches slip freely and harmlessly as the plates grind past one another, while others remain firmly stuck, a state scientists call locking. These locked patches, known as asperities, behave like patches of high friction that resist motion even as the plates around them continue to converge. Year after year, the surrounding crust deforms to accommodate the accumulating mismatch, storing elastic strain energy like a slowly winding spring. When the frictional resistance finally fails, that stored energy is released in seconds as a major earthquake. The Riverside team’s insight was that this years-long accumulation phase leaves a measurable fingerprint at the surface, and that fingerprint can be detected and mapped before rupture.

The fingerprint in question comes from GPS. Networks of continuously operating ground stations scattered across subduction-zone coastlines record millimeter-scale movements of the Earth’s surface as the crust flexes in response to a locked fault below. Where the interface is locked, coastal points are typically dragged landward or held back relative to the motion expected from steady plate convergence; where it creeps, the surface moves more freely. By inverting these subtle deformation patterns, the researchers developed a new algorithm that distinguishes locked, energy-storing portions of the fault from freely slipping ones. The technique does not require dense or perfect data. As Periollat noted, the team had identified the Kamchatka asperity based on a relatively limited data set, and the model still pinpointed the region that later broke. That robustness matters, because offshore faults are among the most poorly instrumented places on the planet, and any forecasting method must work with sparse observations to be practically useful.

The Kamchatka test case is the study’s most dramatic result. Before the rupture, the researchers’ analysis of GPS measurements had highlighted a specific locked region beneath Russia’s Kamchatka Peninsula as a zone of concentrated strain accumulation. The timing of any eventual earthquake was beyond the scope of their method, and they make no claim to have predicted when the event would occur. But when a giant earthquake did strike in 2025, the rupture occurred exactly where the model indicated strain had built up. “We had an idea where the strain was accumulating based on a relatively limited data set,” Periollat said. “Seeing it work so well confirmed that this approach has real potential.” For a field in which retrospective claims of predictive success are often contested, a documented forecast of a rupture location made before the event stands out as an unusually clean validation.

Equally instructive is what the comparison with history revealed about how differently earthquakes can behave at the same boundary. Kamchatka experienced giant earthquakes in both 1952 and 2025, yet the more recent event generated a much smaller tsunami. The researchers interpret this contrast as evidence that the shallowest portion of the fault slipped less during the 2025 rupture than it did during the earlier one. Because shallow slip on a subduction interface displaces the seafloor most efficiently and is a primary driver of tsunami generation, this difference in rupture behavior carries direct hazard implications. The finding underscores a key limitation the team is careful to state: the method identifies where strain is stored, not how a future rupture will distribute that slip, and it cannot predict tsunami size or earthquake timing. Even so, knowing which segments of a megathrust are primed to break narrows the search space for hazard assessment considerably.

The practical value of location-specific forecasting lies in long-term planning rather than short-term warning. Building codes, evacuation routes, tsunami inundation maps, and infrastructure investment all depend on knowing which stretches of coastline face the greatest risk. A method that can rank fault segments by their stored strain gives planners a physically grounded basis for prioritizing those decisions. The Riverside team is already extending the approach to other major subduction zones in Japan, Mexico, New Zealand, and the Pacific Northwest, each of which presents additional complications. Some of these boundaries host slow-slip events that release accumulated energy gradually rather than in sudden earthquakes, blurring the simple distinction between locked and creeping patches and demanding more sophisticated models of how strain builds and dissipates over time.

The researchers are also testing whether the same framework can be applied to continental transform faults closer to home. “In the Bay Area, the Hayward Fault has both creeping and locked sections, much like subduction zones,” Funning said. “We’re investigating whether we can identify the parts most likely to generate future earthquakes.” The parallel is scientifically appealing: if the algorithm can resolve locking along a strike-slip fault threaded through a densely populated urban corridor, it could sharpen hazard estimates for millions of residents. The Hayward Fault creeps visibly along much of its length, yet historical and paleoseismic evidence shows it is capable of large ruptures, making the question of which sections remain truly locked a matter of considerable public interest. Translating a method validated at a subduction zone to a different faulting regime will require careful adaptation, but the underlying principle of mapping surface deformation to infer subsurface locking applies broadly.

Scaling the technique globally, however, runs into an observational bottleneck. GPS stations on land provide valuable information, but many of the world’s most dangerous faults lie offshore, beneath the ocean, where conventional ground-based geodesy cannot reach. Researchers in Japan have begun deploying acoustic instruments on the seafloor that combine underwater sound travel times with satellite positioning to measure slow deformation over many years, and similar efforts are being proposed for Chile and the Pacific Northwest. A recently launched satellite could eventually provide additional geodetic measurements in regions that currently lack GPS coverage, extending the method’s reach. The work also highlights major gaps in global earthquake monitoring, particularly in parts of the Pacific where sparse data make it difficult to assess tsunami hazards. Closing those gaps, the authors argue, is a prerequisite for applying strain-based forecasting to the plate boundaries that need it most.

Throughout, the team is emphatic that improved forecasting of earthquake locations should complement, never replace, public preparedness. “Your peace of mind shouldn’t come from believing we can forecast the exact earthquake,” Funning said. “Especially where we live in Southern California, it’s not a matter of if, but when. There is no substitute for preparation.” That message frames the study’s contribution in sober terms: the algorithm cannot tell a coastal community the day or hour of the next megathrust rupture, but it can tell engineers and emergency managers where the spring is most tightly wound. “We can identify where large earthquakes are likely to occur, even if we can’t predict exactly when,” Periollat said. “With better observations and continued monitoring, we can learn much more about Earth’s most dangerous faults.” If the method performs as well at other plate boundaries as it did at Kamchatka, the map of future earthquake risk may become considerably sharper, one locked asperity at a time.

Subject of Research: A GPS-based method for identifying locked fault segments where giant subduction-zone earthquakes are likely to rupture

Article Title: New method predicts where massive earthquakes will strike

Article References: New method predicts where massive earthquakes will strike. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: earthquakes, subduction zones, GPS geodesy, fault locking, asperities, Kamchatka, tsunami hazard, seismic forecasting, tectonic strain, Hayward Fault, Geophysical Research Letters, disaster preparedness

The appeal of the roughness hypothesis has always been its visual plausibility. On bathymetric maps, the regions that hosted the largest recorded events, such as the 2004 Sumatra-Andaman earthquake, the 2011 Tohoku-oki earthquake and the 1960 Chile earthquake, seemed to correspond with segments of the trench where the incoming plate appeared comparatively smooth, while rougher margins like those off parts of Central America or Japan’s older, sediment-starved trenches appeared to host only smaller ruptures. The hypothesis gained traction because it offered a forecasting shortcut: map the bumps on the incoming seafloor before it subducts, and you have a proxy for how big the next earthquake might be. Given how difficult it is to observe the actual fault surface at depth, a measurable feature of the incoming plate seemed like a gift to hazard assessors.

Yet the new analysis finds that the correlation does not survive careful, systematic testing. Rather than selecting a few celebrated case studies, the researchers assembled a globally consistent dataset of subduction interface roughness measurements and paired them with the maximum magnitudes that each margin has produced, constrained by both the instrumental record and, where available, longer historical and paleoseismic evidence. When roughness is quantified uniformly, using the same topographic measures of relief and spectral character of the incoming plate across all margins, the expected relationship between rough interfaces and smaller maximum earthquakes largely dissolves. Margins with prominently rough seafloor have still hosted very large ruptures, and some comparatively smooth margins have not delivered the giant events the hypothesis predicts.

The authors point to several reasons why the roughness argument fails as a control on maximum magnitude. First, roughness measured on the seafloor before subduction is a poor guide to roughness on the fault at seismogenic depths. As the plate descends, sediments smear into the topographic lows, hydrothermal alteration and mineralization smooth the interface, and the accretionary wedge redistributes material. A seamount that looks like a formidable asperity at the trench may be largely underplated or subducted within a weaker sedimentary layer by the time it reaches the depth range where most seismic moment is released. The fault that ruptures in a great earthquake is therefore not the same rough surface that geodesists and marine geologists can map from shipboard sonar.

Second, and more fundamentally, the maximum magnitude of an earthquake is a geometric property of how large a patch can rupture in a single event, and that is governed by the lateral and downdip extent of coherent locking, the segmentation of the plate boundary, and the accumulated slip deficit, not necessarily by centimeter- to kilometer-scale frictional heterogeneity. A rupture can jump or creep past small roughness elements if the surrounding fault is sufficiently stressed and strongly coupled. Conversely, a smooth interface segmented by major structural boundaries such as fracture zones or tear faults may still be unable to host a rupture longer than a few hundred kilometers. In other words, the features that actually arrest or release ruptures may be structural discontinuities with tens of kilometers of offset, not the relief on the incoming plate.

The study also re-examines the physical reasoning behind the roughness hypothesis itself. Laboratory and theoretical work shows that roughness influences frictional behavior most strongly at small scales, affecting the onset of slip and the distribution of aftershocks, but its effect on the total possible rupture area at the scale of a magnitude 9 earthquake is weak. Scale considerations matter: an earthquake of magnitude 9 ruptures a fault patch on the order of 1,000 kilometers long. Bumps a few kilometers across are simply too small to act as persistent barriers to a rupture front carrying enormous elastic strain energy. The comparison is often made to tearing a sheet of paper: small wrinkles in the paper do not determine where the tear stops if the sheet is being pulled hard enough.

The implications for seismic hazard assessment are significant and uncomfortable. If roughness cannot be used as a proxy for maximum magnitude, then several regional hazard models that factor seafloor roughness into estimates of maximum credible earthquakes may need revision. The study suggests that hazard assessors should weight other evidence more heavily, including geodetic measurements of plate coupling, the distribution of past ruptures inferred from historical accounts, tsunami deposits and coral microatolls, and the structural segmentation of the margin. It also argues against using roughness to declare any margin incapable of producing a giant earthquake, a conclusion with direct consequences for coastal communities and infrastructure planning along subduction margins worldwide.

The researchers emphasize that their findings do not make roughness irrelevant. Interface roughness still shapes where within a rupture the largest slip occurs, how strong ground shaking is distributed, and possibly the frequency of smaller-to-moderate events. What the study removes is the assumption that roughness imposes a hard ceiling on rupture size. Distinguishing between influences on the distribution of slip and controls on maximum magnitude is, the authors argue, an essential refinement that the field has too often glossed over. Roughness may sculpt the earthquake, but it does not appear to cap it.

The work also speaks to a broader lesson in earthquake science: the danger of building predictive frameworks on a handful of spectacular case studies. The great earthquakes of the past century are few, and any global comparison involving fewer than a dozen giant events is statistically fragile. By compiling a comprehensive, uniformly processed dataset, the new study reduces the risk of overfitting a narrative to memorable examples. The result is a more sober picture of subduction zones: nearly any of them, given enough time to accumulate strain, may be capable of producing ruptures larger than their instrumental records suggest, and the size of the next great earthquake may be far less predictable from the shape of the seafloor than scientists once hoped.

Subject of Research: The relationship between subduction interface roughness and maximum earthquake magnitude

Article Title: Maximum earthquake magnitude unlikely to be controlled by subduction interface roughness

Article References: Yang, X., Bell, R. E., Whittaker, A. C., Xu, H., Han, X., Knowlson, A. R., & Locher, V. A. (2026). Maximum earthquake magnitude unlikely to be controlled by subduction interface roughness. Nature Geoscience. https://doi.org/10.1038/s41561-026-02093-z

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02093-z

Keywords: subduction zones, earthquake magnitude, fault roughness, seismic hazard, plate boundaries, rupture segmentation, great earthquakes, seismology, geodesy, marine geophysics, plate coupling, Nature Geoscience