Fact-Checked Science

Paper check · Subduction and the deep Earth

One rock went down warm, then cold: a refinement of subduction science, not its overthrow

A team led by the University of Göttingen found inclusions of coesite about 1 to 4 µm across, one imaged at about 0.8 µm [1 §2][1 Fig. 2b], in the garnets of one eclogite from Goodenough Island, Papua New Guinea, and dated its descent to between about 5.1 and 4.4 million years ago [1 §4]. The rock went down warm, about 470 °C per gigapascal to a depth of about 45 km, then cold and fast to beyond 90 km, and the authors propose that such warm shallow paths “might be common”, “challenging the use of peak metamorphic conditions as indicator for subduction zone thermal structures” [1 Abstract][1 §5]. S0 News said a “tiny freckle on a micro crystal” had beaten “a century of science” and “your entire concept of rock solid subduction, temperature by depth” [2 00:02:07–00:02:21]. The rock’s deep conditions agree with the models of mature subduction zones [1 §5], and the shallow mismatch it documents is an old one, put in global numbers in 2015 [4 §3.4][5].

Our analysis marks our own reasoning and calculations, and the verdicts are our assessment. Everything else comes from the numbered sources; for the paper [1], § gives one of its five main-text sections, Fig. a figure or panel, ED Fig. an extended-data figure and Methods its methods section.

Download as PDFSourcesHow this was checkedCorrections

A real finding about one rock, inflated into the fall of subduction science

Our assessment

~0.8 µm

the coesite inclusion imaged in the paper’s Fig. 2b, in the middle growth zone of a garnet; the press release calls it “a bit smaller than 1 µm” [1 Fig. 2b][3 garnet caption]

700 °C at ~45 km

the garnet core’s conditions, 1.5 ± 0.2 GPa and 700 ± 20 °C: warmer than thermal models of mature subduction zones predict at that depth [1 §3][1 §5]

4.4 Ma

million years ago, the age of the rock’s ultrahigh-pressure stage by 17 zircon dates, “among Earth’s youngest coesite eclogites” [1 §4]

The claims

Overstated

Speck beats a century of science

The presenter: “Could you imagine a century of science, the world’s best, billions and grant money, and your entire concept of rock solid subduction, temperature by depth, gets beat like they’re the Miami Dolphins by a tiny freckle on a micro crystal.” [2 00:02:07–00:02:21]

The rock’s deep conditions “agree with model predictions of the thermal structure of mature, modern-style subduction zones”; only its shallow stage is warmer, and the authors explain that by shear heating or a young subduction zone, “both being an inherent feature of all subduction systems” [1 §5]. What they challenge is reading a whole path from the peak point alone [1 Abstract]. Shear heating has been invoked for subduction zones hotter than expected since at least 1978, and a 2015 compilation found exhumed rocks 100–200 °C warmer on average than the models [4 §3.4][5]. The modern subduction concept dates from 1968 [11], and this is one sample [1 §2]. See sections 4 and 5.

Overstated

No idea what’s beneath our feet

The presenter: “They really don’t know what’s going on beneath our feet any more than they do in the vastness of space.” [2 00:02:24–00:02:30]

The finding was made with settled knowledge: at this rock’s temperatures coesite forms only above 2.7–3.0 GPa [1 §3], zircon can be dated by uranium–lead decay [1 §4], and the rock’s peak conditions landed where the models of mature subduction zones put them [1 §5]. The paper names its open questions, the size of shear heating and the thermal state of young subduction zones [1 §5], and cites a 2024 paper titled “We don’t know when plate tectonics began” for the debate on when plates began [1 §1][14]. Open questions at the edge are not ignorance of the whole. See section 4.

True

One-micrometre coesite speck

The presenter: “We are going down to one micrometer. That’s the size of many bacteria. That’s way smaller than a red blood cell.” [2 00:01:52–00:02:02]

The imaged coesite inclusion is “~0.8-µm sized”, and the single-mineral coesite inclusions range from about 4 to 1 µm [1 Fig. 2b][1 §2]. An E. coli cell is 1.26 µm wide and 1.6–3.9 µm long [15]; a red blood cell is about 7–8 µm across [16]. See section 2.

Context

Hard to be a scientist today

The presenter: “Oof. Hard to be a scientist these days.” [2 00:02:21–00:02:24]

An aside. The authors chose the region because it is “ideally suited to test such assumptions”, tested one and reported the result [1 §1]; the journal’s named reviewers include the first author of the 2015 compilation subtitled “Rocks are hotter than models” [1 Peer review][5]. Testing an assumption, and reporting where it fails, is the job. See section 6.

True

Warm, then cold: one rock’s path

The paper: “the prograde subduction path was initially warm with an apparent average gradient of approximately 470 °C GPa−1 until a depth of about 45 km, before switching to a cold thermal regime with ultrafast subduction rates of greater than 60 mm yr−1 (vertical component)” [1 Abstract]

Quartz-only garnet cores and coesite-bearing garnet mantles, melt inclusions, zirconium-in-rutile thermometry and 17 zircon dates support the path, and the arithmetic checks out [1 §3][1 §4]. The caveats, mostly the authors’ own: one lens [1 §2]; a core pressure estimated from the garnet’s composition, under an upper limit of 2.5 GPa from its melt inclusion, and only a lower bound for the mantle [1 §3]; lithostatic depths; and ages tied to garnet zones by trace elements [1 §4]. See section 3.

Partly true

Contradicts the thermal models

Jan Schönig, in the University of Göttingen’s press release: “Finding this change from hot to cold subduction was unexpected as this contradicts our current thermal models of subduction zones.” [3 ¶2]

True of the steady-state models of mature subduction zones: the core sits “clearly offset to higher temperatures” from them [1 §5][1 Fig. 5a], and models with shear heating stay, in the paper’s words, “still lower than those observed” [1 Fig. 5b]. But a time-dependent model of a young subduction zone, by co-author Holt and Condit, gives paths with “strong similarities” to the rock’s [1 Fig. 5c][8], and the paper words its challenge more narrowly, as one to the use of peak conditions [1 Abstract]. See section 6.

1 · The paper

One eclogite from Goodenough Island, read inclusion by inclusion

Eight authors from Göttingen, Cagliari, Potsdam, Graz, Perth, Syracuse and Miami studied a single partially retrogressed eclogite from the D’Entrecasteaux Islands of eastern Papua New Guinea, the youngest place on Earth where ultrahigh-pressure rocks are known to be coming back to the surface [1 §1][1 Author information]. The paper is a peer-reviewed article in Nature Geoscience, published on October 1, 2026, open access under CC BY 4.0 [1]. Its question is whether a rock that reached the depths of coesite, more than 90 km down, got there along a cold path, as the field assumes from such rocks’ peak conditions: “direct evidence for such a cold prograde evolution is lacking” [1 Abstract].

  1. Step 1

    Sampling

    Sample 03118b is “a small lens (~90 cm × 30 cm) of mafic bulk rock composition hosted by felsic gneiss” at eastern Taleba Bay on southern Goodenough Island (9° 28′ 55.84′′ S, 150° 14′ 46.97′′ E), collected by co-author Suzanne Baldwin [1 §2][1 Contributions]. It is a retrogressed eclogite: about 38% garnet and 22% quartz, with omphacitic clinopyroxene, plagioclase, amphibole, rutile and minor phases [1 §2].

  2. Step 2

    Looking inside garnet

    Raman spectroscopy mapped the inclusions in the garnets at steps of 0.20 µm with a 100× objective [1 Methods]. Quartz, rutile, zircon, apatite and crystallized melt droplets occur in all three growth zones of the garnets; coesite only in the middle zone, the mantle; omphacite and titanite only in the core [1 §2].

  3. Step 3

    Pressure and temperature

    Coesite in the mantle but only quartz in the core means the core grew below, and the mantle above, “the quartz-to-coesite transition located at 2.7–3.0 GPa” [1 §3]. A melt inclusion from the core re-homogenized to glass at 750 °C and 2.5 GPa, an upper limit for the core’s pressure; zirconium in rutile gave temperatures for both zones, and the garnet’s composition fixed the core’s pressure: 1.5 ± 0.2 GPa and 700 ± 20 °C for the core, above 2.9 GPa and 800 ± 25 °C for the mantle [1 §3][1 Fig. 3].

  4. Step 4

    Dating

    Zircon occurs only as micrometre-sized inclusions in garnet, freed by dissolving the host; 26 of the largest grains, about 10 µm across, were dated by ion microprobe and 17 kept, after excluding eight with too little radiogenic lead and one probably foreign grain of 362 million years [1 §4][1 Methods]. They span 5.1 ± 0.2 to 3.1 ± 0.1 million years ago, and their trace elements tie the 5.1–4.4 interval to the climb from ≤1.5 to >2.9 GPa, “corresponding to vertical rates of >60 mm yr−1 assuming lithostatic pressure” [1 §4].

Zircon grains dated

17 of 26

eight excluded for radiogenic lead below 90%, and one grain of 362 ± 6 million years, “probably xenocrystic”; uncertainties are two standard errors [1 Methods]

The descent

5.1 → 4.4 Ma

million years ago: from garnet-core growth at ≤1.5 GPa to the ultrahigh-pressure stage at >2.9 GPa, by the zircons’ heavy rare earths [1 §4]

Vertical rate

>60 mm/yr

“ultrafast” in the paper’s words: the vertical component of the descent, assuming lithostatic pressure [1 Abstract][1 §4]

In the authors’ words

UHP is ultrahigh pressure; HP high pressure; T/P the ratio of temperature to pressure, a geothermal gradient; GPa a gigapascal, the pressure under about 30 km of rock.

“This assumption stems from peak metamorphic conditions of ultrahigh-pressure rocks, which indicate low subduction gradients of less than 330 °C GPa−1, but direct evidence for such a cold prograde evolution is lacking.”

“Our most recent example of HP and UHP rock formation and exhumation occurs in eastern Papua New Guinea, making the region ideally suited to test such assumptions.”

“The age of about 4.4 Ma for the UHP stage is considerably younger than the about 7.1 Ma UHP stage for the coesite eclogite at Tumagabuna Island, making it among Earth’s youngest coesite eclogites.”

Schönig et al. (2026) [1 §4], on the earlier find at Tumagabuna Island [9]

What the paper says it challenges, in full

The abstract ends: “We propose that high geothermal gradients at shallow levels of deep subduction systems might be common, challenging the use of peak metamorphic conditions as indicator for subduction zone thermal structures.” [1 Abstract] The introduction puts it as a caution: the authors “caution that peak metamorphic conditions alone may be insufficient to deduce the prograde evolution and geodynamic setting of UHP rock formation” [1 §1]. And the conclusion: the identification of UHP metamorphism “is not unequivocal evidence that subduction was exclusively cold and evolved at a thermal steady state” [1 §5]. The target is an inference, reading a rock’s whole descent from its deepest point; not subduction, and not the models of its deep thermal structure, with which the rock agrees (section 4).

2 · The speck

A micrometre of coesite, and how it was found

The presenter opened with the scale bar: “Folks, can you see those scales there? We are going down to one micrometer. That’s the size of many bacteria. That’s way smaller than a red blood cell.” [2 00:01:49–00:02:02] That part is right. The paper’s Fig. 2b shows “an ~0.8-µm sized monomineralic coesite below the polished surface” of a garnet, and the press release’s caption calls the inclusion “a bit smaller than 1 µm – approximately one hundredth the size of the thickness of a human hair” [1 Fig. 2b][3 garnet caption].

The imaged coesite

~0.8 µm

single-mineral coesite inclusions are about 4 to 1 µm; two-mineral coesite-and-quartz inclusions about 5 µm [1 Fig. 2b][1 §2]

An E. coli cell

1.26 µm

wide, in every growth condition measured (± 0.16 µm), and 1.6 to 3.9 µm long depending on how fast the cells grow [15]

A red blood cell

7–8 µm

“quite small cells, with a mean diameter of only about 7–8 micrometers”, a standard textbook value [16]

How a speck that small is found

The presenter skipped “the obvious question of how these guys even spotted this tiny little speck” [2 00:02:00–00:02:07]. The Methods answer it. A Raman spectrometer with a 532-nm laser and a 100× objective scanned the garnets point by point, “a step size of 0.20 µm × 0.20 µm (Fig. 2b)” with one second per spectrum, to map inclusions sitting below the polished surface [1 Methods]. Coesite has its own Raman spectrum, shown beside the image [1 Fig. 2b]. The team knew where to look: co-authors had already found coesite inclusions in garnet from the beach placer beside the sample site, in 2021 [10 Abstract][1 §1].

Our analysis

Why a speck can carry a result

Coesite is a form of silica that, at this rock’s temperatures, forms only at pressures above the quartz-to-coesite transition, 2.7–3.0 GPa [1 §3], which is why the paper calls it “Unequivocal evidence for UHP metamorphism” [1 §1]. Sealed inside a garnet, a grain of it tends to survive the trip back up when it is small: the paper reads the larger size of the two-mineral coesite-and-quartz inclusions, about 5 µm, against 4 to 1 µm for the all-coesite ones, as agreeing “with the hypothesis that inclusion size is an important factor for coesite preservation” [1 §2]. So smallness is not a weakness of the evidence but a condition of its survival. What the paper adds to the speck is where it sits: in the garnet’s middle zone, with quartz-only cores inside it, which turns a single pressure marker into a sequence [1 §3].

3 · The path

Warm to about 45 km, then cold and fast to beyond 90 km

Put the two garnet zones on a pressure–temperature diagram and the shape of the descent appears: the core formed at about 700 °C and 1.5 GPa, warmer than thermal models of mature subduction zones allow at that depth; the mantle formed at about 800 °C above 2.9 GPa, where those models put it [1 §3][1 §5]. The chart is adapted from the paper’s Figs. 3 and 5a, from the numbers in its text.

The rock’s pressure–temperature path beside the models

  • Garnet core and garnet mantle: the rock’s own record, and the path between them
  • About 470 °C/GPa: the apparent gradient of the shallow path
  • About 280 °C/GPa: what the peak conditions alone would imply
  • Thermal models of mature subduction zones, as the paper summarizes them: gradients below 300 °C/GPa (darker), up to 400 for the warmest (lighter)
  • The quartz–coesite transition at this rock’s temperatures, 2.7–3.0 GPa; coesite is stable above it

Scroll sideways to see the whole chart →

thermal models of mature subduction zones (ref. 6 in the paper): gradients below 300 °C/GPa; the warmest, Mexico and Cascadia, below 400 coesite stable above the quartz → coesite transition, 2.7–3.0 GPa ≈280 °C/GPa: what the peak conditions alone imply ≈470 °C/GPa: the apparent gradient of the shallow path ≈45 km ≈90 km garnet core: 1.5 ± 0.2 GPa, 700 ± 20 °C from about 5.1 million years ago garnet mantle: >2.9 GPa, 800 ± 25 °C by about 4.4 million years ago 0 200 400 600 800 1,000 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 Pressure (GPa) Temperature (°C)

Temperature in degrees Celsius; pressure in gigapascals, with the paper’s depths for lithostatic pressure

~470 °C/GPa

the apparent average gradient of the shallow path, to about 1.5 GPa [1 Abstract][1 §5]

~280 °C/GPa

the gradient “implied by the peak conditions alone” [1 §5]

<300 °C/GPa

“typical” gradients in models of “mature, modern-style subduction zones” [1 §1]

<400 °C/GPa

“even the warmest subduction zones (Mexico and Cascadia)” in those models [1 §1]

Adapted from Fig. 3 and Fig. 5a of [1], CC BY 4.0, and drawn from the numbers in its text rather than traced: the garnet core at 1.5 ± 0.2 GPa and 700 ± 20 °C, the garnet mantle above 2.9 GPa at 800 ± 25 °C, the quartz-to-coesite transition at 2.7–3.0 GPa, the apparent shallow gradient of about 470 °C GPa−1 and the peak-only gradient of about 280 [1 §3][1 §5]. The transition is the paper’s range at this rock’s temperatures, drawn flat across the chart; in the paper’s Fig. 3 the boundary rises slightly with temperature [1 Fig. 3]. The gray wedges stand for the paper’s summary of thermal models of mature, modern-style subduction zones, “typical temperature−pressure (T/P) gradients <300 °C GPa−1” with the warmest, Mexico and Cascadia, “<400 °C GPa−1” [1 §1]; they are drawn as straight lines of constant gradient, where the paper’s Fig. 5a shows the models as curved paths that warm sharply near 2.5 GPa. The mantle’s shape is open at the top because its pressure is a minimum, and the dates come from the zircons [1 §4]. Depths are the paper’s [1 Abstract][1 §5].
Our analysis

The arithmetic, and what the numbers rest on

The paper’s gradients follow from its two points: 700 °C over 1.5 GPa is 467 °C per gigapascal, the “~470” of the text, and 800 °C over 2.9 GPa is 276, its “~280” [1 §3][1 §5]. For a rock density of 3,300 kg/m³ under lithostatic pressure, 1.5 GPa is about 46 km and 2.9 GPa about 90 km, the paper’s “about 45 km” and “depths >90 km” [1 Abstract][1 §5]. The climb from 1.5 to 2.9 GPa, about 43 km, took 0.7 million years (zircon ages of 5.1 and 4.4 Ma): about 62 mm per year, the paper’s “>60 mm yr−1” [1 §4]. Everything checks.

The limits are the authors’ own. The melt-inclusion experiment, run at 750 °C and 2.5 GPa, gives only an upper limit for the core’s pressure: “the pressure represents a maximum value”. The core’s best estimate of 1.5 ± 0.2 GPa comes from where the garnet’s composition isopleths cross the core’s temperatures [1 §3]. The mantle’s pressure is only a lower bound, because the isopleths “do not enable conditions of garnet mantle formation to be determined” above 1.8 GPa, and the solidus’s location “is uncertain at pressures >1.8 GPa” [1 §3]. The paper calls the gradient an “apparent average” [1 Abstract]. The depths and the rate assume lithostatic pressure [1 §4]. The zircons were dated as separates, so their ages are tied to the garnet zones by their trace elements, “with reasonable certainty” [1 §4]. None of this undoes the sequence, which rests on the simplest observation in the paper: quartz only in the cores, coesite in the mantles [1 §3].

4 · New and standing

What is new, and what stands

The paper joins a long-running debate, and cites it: subduction-related rocks “often record warmer conditions than models predict” [1 §5], a mismatch argued over for decades and put in global numbers in 2015 (section 5). What it brings is a rock that holds both a warm shallow stage and a cold deep stage, with dates.

New

A direct record of the shallow part of an ultrahigh-pressure rock’s descent, where “mineralogical evidence clearly ascribed to the prograde path is generally lacking” [1 §1]. A change of gradient inside one rock, warm then cold, rather than a single peak point; the paper reads the same shape in prograde paths reported from two older ultrahigh-pressure terranes, Sanbagawa and Dabie–Sulu [1 §5][1 Fig. 5a]. One of the youngest coesite eclogites known, at about 4.4 million years, with zircon ages that time its descent between the two stages [1 §4]. The press release’s own placing of the novelty: “Rocks that record higher temperature conditions than expected have been reported before, but this is the first time a journey from a shallow to a deep depth with an accompanied change to relatively colder gradients has been shown so clearly.” [3 ¶3]

Standing

Subduction itself: the paper’s setting is “the subduction of a rifted continental fragment of Australian provenance” [1 §1], and it notes that the Pacific and Australian plates converge obliquely today “at a rate of ~ 110 mm yr−1”, above the 100 mm per year beyond which, it says, temperatures “might have been even higher” with shear heating [1 §5]. The deep thermal structure: the mantle’s conditions “agree with model predictions of the thermal structure of mature, modern-style subduction zones” [1 §5]. Coesite as a depth marker and zircon as a clock, the tools of the result [1 §3][1 §4]. And the two explanations the authors offer are not new physics but mechanisms “both being an inherent feature of all subduction systems” [1 §5]; both were already weighed in 1992 [18 Abstract].

The debate the paper joins

Three statements of the same mismatch, from 2018 and 2023.

“Some commonly referenced thermal-mechanical models of current subduction zones imply temperatures that are 100-500 °C colder at 30-80-km depth than pressure-temperature conditions determined thermobarometrically from exhumed metamorphic rocks.”

Kohn et al. (2018), whose title says shear heating “reconciles” them [6 Abstract]

“Compilations of the maximum pressure-temperature conditions in exhumed rocks indicate significantly warmer conditions than those predicted by thermal models. This could be due to preferential exhumation of rocks from hotter conditions … Alternatively, the models might underestimate the forearc temperatures by neglecting certain heat sources.”

“The lawsonite eclogite data and a significant proportion of the blueschist data fall within the range of global models, but the near-steady-state models cannot explain the warmest exhumed rock data.”

van Keken and Wilson (2023), the review the paper takes its model gradients from [4 §3.4]

Explanation 1: friction on the plate boundary

Shear heating: rock sliding past rock at the plate interface dissipates energy as heat. The paper notes that “it has been argued that thermal models of mature subduction zones incorrectly neglect such important heat sources” [1 §5]; Kohn and colleagues found that forearc heat flow needs “An additional heat source--likely shear heating”, and that “Adding such shear stresses to thermal models quantitatively reproduces the pressure-temperature conditions recorded by exhumed metamorphic rocks” [6 Abstract]. Others find that “reasonable amounts of shear heating” warm the oceanic crust by “less than 50 °C”, an effect “rather skin deep”, and infer that typical exhumed blueschists and eclogites come from warm settings, young oceanic plates or “the warmer initial stages of subduction” [7 Abstract][4 §3.4]. In the paper’s Fig. 5b, the thermal structure with shear heating “more closely resembles the shape of the observed P–T path, but modelled temperatures are still lower than those observed”; they “might have been even higher” if convergence exceeded 100 mm per year, friction was unexpectedly high or the sinking plate was young [1 §5][1 Fig. 5b].

Explanation 2: a subduction zone still young

Thermal models of mature subduction zones describe a steady state. A time-dependent model by co-author Holt and Condit finds that pressure–temperature conditions on the slab “exhibit substantial variability throughout the various phases of subduction: initiation, free sinking, and mature subduction” [8 Abstract]. In the paper’s Fig. 5c the observed path has “strong similarities to the numerically predicted paths for rocks entering the subduction zone ~ 5–12 Myr after subduction initiation”, and the rate of over 60 mm per year fits rocks entering 8–13 million years after it [1 §5][1 ED Fig. 4]. The authors note small differences, “the observed path is slightly steeper and the modelled timescale somewhat longer”, which they expect from an “idealized, non-site-specific” model [1 §5]. The idea is not new either: in 1992 Peacock allowed for blueschist terrains that “represent material metamorphosed during the early stages of subduction (first 5–10 m.y.) while the thermal structure is still evolving” [18 Abstract].

Our analysis

Which explanation? The paper does not choose, and does not need to

The authors call the two “reasonable and non-mutually exclusive mechanisms” and leave “quantifying the relative contributions of the two mechanisms” outside the study [1 §5]; other work disputes how much heat friction can deliver to the slab’s crust [7 Abstract][4 §3.4]. Either way, the shallow warmth has an explanation in known physics, which is the opposite of the episode’s picture. What the paper does change is a habit of inference: a cold peak point no longer proves a cold path, so “considering T/P gradients of peak metamorphic conditions for reconstructing Earth’s geodynamic evolution might be misleading” [1 §5]. That matters for arguments about when modern, cold subduction began on Earth, which lean on such gradients [1 §1]; it does not question that plates subduct today, and at depth this rock agrees with the models.

5 · A century?

“A century of science”: a 58-year-old idea and a decades-old debate

The presenter counted “a century of science, the world’s best, billions and grant money” on the losing side [2 00:02:07–00:02:12]. In the dated record, from the sources we could resolve, the subduction concept and its thermal models are younger than a century, the observations behind them a little older, and the mismatch this paper adds to has been argued over for decades.

The years behind the claim

YearWhat happenedSource
1928–35Kiyoo Wadati shows earthquakes deeper than 300 km beneath central Honshu in 1928, and in 1935 draws “lines representing equal focal depths of the intermediate and deep earthquakes, which would later be called the Wadati-Benioff zone”[20]
1968“Seismology and the new global tectonics”, by Isacks, Oliver and Sykes, in the Journal of Geophysical Research: the plate-tectonics revolution in which the modern concept of subduction was framed[11]
1970–71Early thermal models of a descending slab: Minear and Toksöz’s “Thermal regime of a downgoing slab and new global tectonics” and Oxburgh and Turcotte’s “Thermal Structure of Island Arcs” in 1970, and “Temperature field and geophysical effects of a downgoing slab” in 1971[19][12], [13]
1978Shear heating is proposed as the cause of arc volcanoes, an early case of what the 2023 review calls “a long-time favorite ad hoc explanation” for subduction zones hotter than expected; arc volcanoes are now explained by the hot mantle wedge and fluids from the slab[4 §3.4]
1990Molnar and England, “Temperatures, heat flux, and frictional stress near major thrust faults”: the paper’s source for shear heating along the plate interface[17][1 §5]
1992Peacock sets the temperatures of blueschists against model paths with and without shear heating: without it, the paths are “very cold”; with low stresses, the rocks could come from “the early stages of subduction”[18 Abstract]
1996–2009Global databases of rocks exhumed from oceanic subduction zones are compiled, later set against the models’ predicted temperatures[4 §3.4]
2008A coesite eclogite at Tumagabuna Island, Papua New Guinea, later dated to an ultrahigh-pressure stage of about 7.1 million years ago[9][1 §4]
2015A global compilation of exhumed blueschists and eclogites, subtitled “Rocks are hotter than models”: on average, the rock record is 100–200 °C warmer below the forearc than the models predict[5][4 §3.4]
2018Forearc heat flow is found to need “An additional heat source--likely shear heating”; a second study finds that reasonable shear heating adds “less than 50 °C” and infers that exhumed rocks come preferentially from warm settings, young plates or the first stages of subduction[6 Abstract][7 Abstract]
2021Time-dependent models show slab conditions varying through “initiation, free sinking, and mature subduction”; coesite is found in garnet sand on the beach beside this paper’s sample site[8 Abstract][10 Abstract]
2023The review the paper draws its model gradients from finds that “the near-steady-state models cannot explain the warmest exhumed rock data”[4 §3.4]
2026This paper: one rock with a warm shallow stage and a cold deep stage, dated, and a proposal that such paths “might be common”[1 Abstract]

For the papers of 1968–71 we read titles and records only: [11], [12], [13] and [19]. For those of 1990 and 1992 we read the abstracts as deposited with Crossref [17][18]. The 2015 paper is cited by its title [5], and the 1978 work as the 2023 review describes it [4 §3.4].

Our analysis

The modern concept of subducting plates belongs to the late 1960s, 58 years before this paper; early thermal models of a downgoing slab to 1970–71, 55 to 56 years [11][19]. Some roots are older, and the presenter’s “century” can fairly reach them: deep earthquakes beneath Japan were shown in 1928, and by 1935 Kiyoo Wadati had drawn the inclined zone they form [20], the zone now used to trace subducting slabs [21 §2.2]. But the thermal models this paper tests descend from those of the 1970s, and at depth the rock agrees with them [1 Fig. 5a][1 §5].

Nor is the rock’s shallow warmth a sudden upset. Shear heating has been invoked for subduction zones hotter than expected since 1978 [4 §3.4], blueschists were set against model paths with and without it in 1992 [18 Abstract], and in 2015 the rock record was found 100–200 °C warmer than the models on average [5][4 §3.4]. The paper’s contribution is to show both stages in one rock, with the time between them [1 §4].

As for “billions and grant money”: the billions were the presenter’s figure for the whole field, which we did not check. The speck itself was found by that same grant-funded science: the study discloses one German Research Foundation grant and the Geosciences Center Göttingen’s support [1 Funding].

6 · Paper to episode

From “challenging the use of peak metamorphic conditions” to “gets beat”, in three steps

The episode linked the press release, not the paper [2][3]. Each step along the way kept the finding and dropped a qualifier.

  1. Step 1

    The paper, October 1, 2026

    “We propose that high geothermal gradients at shallow levels of deep subduction systems might be common, challenging the use of peak metamorphic conditions as indicator for subduction zone thermal structures.” [1 Abstract]

    “… and caution that peak metamorphic conditions alone may be insufficient to deduce the prograde evolution and geodynamic setting of UHP rock formation.” [1 §1]

    A proposal (“might be common”) and a caution, aimed at a method of inference.

  2. Step 2

    The press release, October 7, 2026

    “This change from very hot subduction at this depth to the more common colder subduction at its final depth challenges scientists’ view of the thermal structure of subduction zones.” [3 ¶1]

    Dr Jan Schönig: “Finding this change from hot to cold subduction was unexpected as this contradicts our current thermal models of subduction zones.” [3 ¶2]

    “The use of peak metamorphic conditions” has become “scientists’ view of the thermal structure”, and “challenging” has become “contradicts”. The release still records that “Rocks that record higher temperature conditions than expected have been reported before”, and it gives the two explanations [3 ¶3].

  3. Step 3

    The episode, October 23, 2026

    “Could you imagine a century of science, the world’s best, billions and grant money, and your entire concept of rock solid subduction, temperature by depth, gets beat like they’re the Miami Dolphins by a tiny freckle on a micro crystal. Oof. Hard to be a scientist these days. They really don’t know what’s going on beneath our feet any more than they do in the vastness of space.” [2 00:02:07–00:02:30]

    In the episode, “thermal models” has become “your entire concept of rock solid subduction”, “contradicts” has become “gets beat”, the earlier reports and the explanations are gone, and a debate decades old has become the defeat of “a century of science”.

The episode’s statements beside the paper

What was saidThe paperOur verdict
The presenter: “We are going down to one micrometer. That’s the size of many bacteria. That’s way smaller than a red blood cell.” [2 00:01:52–00:02:02]The imaged inclusion is “~0.8-µm sized” [1 Fig. 2b]; an E. coli cell is 1.26 µm wide [15] and a red blood cell about 7–8 µm across [16].True
The presenter: “a century of science … your entire concept of rock solid subduction, temperature by depth, gets beat … by a tiny freckle on a micro crystal” [2 00:02:07–00:02:21]The deep conditions agree with the models, and the shallow mismatch has known explanations [1 §5]; it has been argued over for decades [4 §3.4]; the challenge is to one inference, from one sample [1 Abstract][1 §2]. The subduction concept dates from 1968 [11].Overstated
The presenter: “Oof. Hard to be a scientist these days.” [2 00:02:21–00:02:24]The authors set out to test an assumption in a region “ideally suited” to it, and did [1 §1].Context
The presenter: “They really don’t know what’s going on beneath our feet any more than they do in the vastness of space.” [2 00:02:24–00:02:30]Built on settled knowledge, and in agreement with the models at depth [1 §3][1 §5]; the open questions it names are the size of shear heating and the thermal state of young subduction zones, and, citing a 2024 paper, when plate tectonics began [1 §1][14].Overstated
The press release, quoting Dr Schönig: “this contradicts our current thermal models of subduction zones” [3 ¶2]True of the steady-state models of mature subduction zones (Fig. 5a), and models with shear heating stay “still lower than those observed” (Fig. 5b) [1 §5]; but the time-dependent model of a young subduction zone gives paths with “strong similarities” to the rock’s (Fig. 5c), and the paper words its challenge more narrowly [8][1 Abstract].Partly true

Quotations from the episode are from its automatically generated transcript, with the time span in which each is spoken [2].

Nothing in the paper is about Earth’s magnetic field

The episode billed the item among “a micro particle doing crazy things and a reminder of why Earth’s magnetic changes we’re seeing are so bad” [2 00:00:17–00:00:27]. The paper has no bearing on that: the words “magnetic” (other than the mineral magnetite), “excursion”, “pole”, “disaster” and “biosphere” do not occur in its text or methods, and “life” only in the titles of works it cites [1]. The press release mentions earthquakes and volcanoes only in passing, as things subduction zones host and as effects to understand better [3 ¶2][3 ¶3]. The episode’s magnetic claims are checked on our pages for the episode of October 23, 2026 and for its other paper, on the rates of change in excursions and reversals.

7 · Read with care

A careful paper about one rock, and a proposal its authors label as such

The measurements are detailed and the writing is hedged. The limits below are our reading, most of which the paper states itself.

What would change these verdicts

  1. 01

    Another rock with a dated prograde record as clear as this one, from Papua New Guinea or elsewhere: the same warm-then-cold path would support “might be common”; a cold path throughout would limit it.

  2. 02

    Site-specific models for Papua New Guinea that weigh shear heating against the young-subduction effect, the quantification the authors leave “beyond the scope of this study” [1 §5].

  3. 03

    For the episode’s version: a result that put the deep, peak conditions outside the thermal models too, or that undercut coesite as a pressure marker or zircon as a clock. This paper does the opposite: it relies on both and agrees with the models at depth [1 §3][1 §5].

Our analysis

Bottom line: a speck of coesite told a careful story about one rock, and the episode told a different story about science.

The rock went down warm, then cold, and the paper shows it well, with its caveats in view [1 §3][1 §4]. That refines how a rock’s whole path should be read from its deepest point [1 Abstract], in a debate the field has been having for decades [4 §3.4], and in global numbers since 2015 [5]. It agrees with the models where the rock went deepest, relies on the very knowledge of the deep Earth that the presenter says scientists lack, and explains the surprise with mechanisms every subduction zone has [1 §5]. Nothing in it is about Earth’s magnetic field.

Sources

Bracketed numbers refer to this list. For the paper [1], Abstract is its abstract; § a section of its main text in the order printed (§1 the untitled introduction, §2 “Retrogressed eclogite from Goodenough Island”, §3 “Pressure–temperature conditions of garnet crystallization”, §4 “U–Pb dating and composition of zircon inclusions”, §5 “Geothermal gradient change during subduction”); Methods its methods section; Fig. a figure or panel; ED Fig. an extended-data figure; Funding, Author information, Contributions and Peer review its end-matter statements. For the episode [2], the time is a position in the video. For the press release [3], ¶ is a paragraph of its text and garnet caption the caption of its photograph of a garnet. For the other sources, § names a section and Abstract the abstract; where only an abstract or a title was read, the entry says so. Links checked October 24, 2026.

  1. 1

    Schönig, J., Ferrero, S., Gresky, K., Sorger, D., Schmitt, A. K., Wunder, B., Baldwin, S. L., and Holt, A. F. (2026). Geothermal gradient change during subduction recorded by ultrahigh-pressure eclogite. Nature Geoscience, published October 1, 2026 (received February 12, accepted August 25, 2026). Open access, CC BY 4.0. We read the full text in the publisher’s HTML and 13-page PDF, with Figs. 1–5 and Extended Data Figs. 1–4, and the figshare record of its Supplementary Table 1 (10.6084/m9.figshare.33306585, CC BY 4.0); its peer-review file is an image-only PDF we did not read. Its Fig. 1 is third-party material and is not reproduced here.

  2. 2

    SpaceWeatherNews (S0) (October 23, 2026). Rapid Pole Shifts Confirmed, Micro-Disruptor | S0 News Oct.23.2026. YouTube video, https://www.youtube.com/watch?v=CTJz20MnklM. Quoted from the automatically generated transcript; times are positions in the video. The label “Micro Rock Mystery” is the episode’s own, on its link to the press release [3] in the video’s description.

  3. 3

    Georg-August-Universität Göttingen (October 7, 2026). Two kinds of subduction recorded in one rock. Research team identify unusual journey deep into Earth for most recently formed ultrahigh-pressure eclogite. Press release No. 132, with a quotation from Dr Jan Schönig and four captioned photographs. The link the episode gave.

  4. 4

    van Keken, P. E., and Wilson, C. R. (2023). An introductory review of the thermal structure of subduction zones: III—Comparison between models and observations. Progress in Earth and Planetary Science 10, 57. CC BY 4.0. The paper’s ref. 6, the source of its model gradients. We read the full text; §3.4 is its comparison with exhumed rocks.

  5. 5

    Penniston-Dorland, S. C., Kohn, M. J., and Manning, C. E. (2015). The global range of subduction zone thermal structures from exhumed blueschists and eclogites: Rocks are hotter than models. Earth and Planetary Science Letters 428, 243–254. The paper’s ref. 39. Title and record only: its publisher’s page showed a bot check to our tools and its Crossref and OpenAlex records carry no abstract; its content is described by [1] and [4].

  6. 6

    Kohn, M. J., Castro, A. E., Kerswell, B. C., Ranero, C. R., and Spear, F. S. (2018). Shear heating reconciles thermal models with the metamorphic rock record of subduction. Proceedings of the National Academy of Sciences 115(46), 11706–11711. The paper’s ref. 40. We read the abstract, from its PubMed record (PMID 30373832) via Europe PMC; the full text is in PubMed Central (PMC6243260), and the journal’s page showed a bot check to our tools.

  7. 7

    van Keken, P. E., Wada, I., Abers, G. A., Hacker, B. R., and Wang, K. (2018). Mafic high-pressure rocks are preferentially exhumed from warm subduction settings. Geochemistry, Geophysics, Geosystems 19(9), 2934–2961. The paper’s ref. 46. Abstract only, as deposited with Crossref; the journal’s page showed a bot check to our tools.

  8. 8

    Holt, A. F., and Condit, C. B. (2021). Slab temperature evolution over the lifetime of a subduction zone. Geochemistry, Geophysics, Geosystems 22(6), e2020GC009476. The paper’s ref. 38, the time-dependent model behind its Fig. 5c; Holt is a co-author of [1]. Abstract only, as deposited with Crossref; the journal’s page showed a bot check to our tools.

  9. 9

    Baldwin, S. L., Webb, L. E., and Monteleone, B. D. (2008). Late Miocene coesite-eclogite exhumed in the Woodlark Rift. Geology 36(9), 735–738. The paper’s ref. 13, the earlier coesite eclogite at Tumagabuna Island. Title and record only.

  10. 10

    Baldwin, S. L., Schönig, J., Gonzalez, J. P., Davies, H., and von Eynatten, H. (2021). Garnet sand reveals rock recycling processes in the youngest exhumed high- and ultrahigh-pressure terrane on Earth. Proceedings of the National Academy of Sciences 118(3), e2017231118. The paper’s ref. 14. We read the abstract, from its PubMed record (PMID 33431682) via Europe PMC; the full text is in PubMed Central (PMC7826392), and the journal’s page showed a bot check to our tools.

  11. 11

    Isacks, B., Oliver, J., and Sykes, L. R. (1968). Seismology and the new global tectonics. Journal of Geophysical Research 73(18), 5855–5899. Title and record only, resolved through Crossref; cited for its date.

  12. 12

    Oxburgh, E. R., and Turcotte, D. L. (1970). Thermal Structure of Island Arcs. Geological Society of America Bulletin 81(6), 1665–1688. Title and record only; cited for its date.

  13. 13

    Toksöz, M. N., Minear, J. W., and Julian, B. R. (1971). Temperature field and geophysical effects of a downgoing slab. Journal of Geophysical Research 76(5), 1113–1138. Title and record only; cited for its date.

  14. 14

    Harrison, T. M. (2024). We don’t know when plate tectonics began. Journal of the Geological Society 181(4), jgs2023-212. The paper’s ref. 3. Abstract only, as deposited with Crossref; the journal’s page showed a bot check to our tools.

  15. 15

    Volkmer, B., and Heinemann, M. (2011). Condition-dependent cell volume and concentration of Escherichia coli to facilitate data conversion for systems biology modeling. PLoS ONE 6(7), e23126. CC BY. We read the full text (Europe PMC, PMC3146540); its Results give the cell dimensions.

  16. 16

    Betts, J. G., et al. (2022). Anatomy and Physiology 2e, §18.3 Erythrocytes. OpenStax, Rice University. A textbook; its page states a Creative Commons Attribution-NonCommercial-ShareAlike licence.

  17. 17

    Molnar, P., and England, P. (1990). Temperatures, heat flux, and frictional stress near major thrust faults. Journal of Geophysical Research: Solid Earth 95(B4), 4833–4856. The paper’s ref. 37, its source for shear heating along the plate interface. Abstract only, as deposited with Crossref; the journal’s page showed a bot check to our tools.

  18. 18

    Peacock, S. M. (1992). Blueschist-facies metamorphism, shear heating, and P-T-t paths in subduction shear zones. Journal of Geophysical Research: Solid Earth 97(B12), 17693–17707. Abstract only, as deposited with Crossref; the journal’s page showed a bot check to our tools.

  19. 19

    Minear, J. W., and Toksöz, M. N. (1970). Thermal regime of a downgoing slab and new global tectonics. Journal of Geophysical Research 75(8), 1397–1419. Title and record only; cited for its date.

  20. 20

    Suzuki, Y. (2001). Kiyoo Wadati and the path to the discovery of the intermediate-deep earthquake zone. Episodes 24(2), 118–123. CC BY-NC 4.0, as its Crossref record states. A history; we read the full text, the journal’s 6-page PDF.

  21. 21

    van Keken, P. E., and Wilson, C. R. (2023). An introductory review of the thermal structure of subduction zones: I—motivation and selected examples. Progress in Earth and Planetary Science 10, 42. CC BY 4.0. The first part of the review in [4]; we read its abstract and §2.2, on seismology.

Method

How this was checked

On October 24, 2026, we checked every statement, number, quotation and citation on this page against its source: the paper’s full text in the publisher’s HTML and PDF, with its five figures and four extended-data figures; its Crossref record and the figshare record of its supplementary table; the University of Göttingen’s press release; the episode’s automatically generated transcript [2]; parts III and I of van Keken and Wilson’s open-access review and Suzuki’s history of the deep earthquake zone, read in full or in the sections cited; an open-access study of E. coli cell size and an open textbook, for the size comparisons; and thirteen other works resolved through Crossref, read as abstracts or titles as the sources list says: about 120 items in all. Eight publisher pages showed bot checks to our tools (the AGU journal pages of Holt and Condit 2021, van Keken et al. 2018, Molnar and England 1990 and Peacock 1992, the Lyell Collection page of Harrison 2024, PNAS’s pages of Kohn et al. 2018 and Baldwin et al. 2021, and ScienceDirect’s page of Penniston-Dorland et al. 2015); we did not bypass them. The two PNAS abstracts were re-read from their PubMed records, the five others are read as deposited with Crossref, and Penniston-Dorland et al. is cited by title only. The paper’s peer-review file is an image-only PDF we did not read. A fresh reviewer then read the built page against the sources, re-resolving every DOI. That check corrected our account of the paper’s Fig. 5b, where the models with shear heating stay too cold and only the time-dependent model resembles the rock’s path; dated the debate the paper joins to decades before the 2015 compilation; moved the paper’s “maximum value” from the core’s 1.5 GPa, where we had put it, to the melt-inclusion experiment’s 2.5 GPa; restored the paper’s comparison with prograde paths from two other ultrahigh-pressure terranes; added the work that disputes how much heat friction delivers to the slab, and the older seismological roots of the subduction concept; and fixed two quotations, several locators and the chart’s outline of the garnet mantle.

The verdicts and the items tagged “Our analysis” are our own assessment, built on the cited sources. They include the arithmetic behind the paper’s gradients, depths and subduction rate; the pressure–temperature chart, drawn from numbers in the paper’s text; the timeline from 1928 to 2026; the cell-size comparison; and our reading of what each step from paper to press release to episode added or dropped.

No measurement was repeated. From the paper’s own numbers we recomputed its apparent gradients (700 °C over 1.5 GPa is 467 °C per gigapascal; 800 °C over 2.9 GPa is 276), its depths (1.5 and 2.9 GPa are about 46 and 90 km for a rock density of 3,300 kg/m³ under lithostatic pressure) and its vertical subduction rate (about 43 km in 0.7 million years is about 62 mm per year). We did not recompute the thermodynamic models, the zircon ages or the thermometry, which depend on the supplementary data and the authors’ software, and we did not digitise the figures: the chart’s positions come from the numbers in the text, and the thermal models of mature subduction zones appear on it as lines of constant gradient, not as the curved model paths of the paper’s Fig. 5a.

This page was researched and written with Claude, an AI assistant made by Anthropic, and then re-checked against its sources in a separate verification pass.

The verdicts on this page

Overstated
Has some basis but claims much more than the evidence supports.
Partly true
Part of the claim holds up and part doesn’t.
True
The evidence supports the claim.
Context
A question or aside rather than a claim. We add the facts it leaves out.

All ten verdicts, and how we check

Corrections

Corrections and changes

If we get something wrong, we fix it and record the change below with the date. If a correction changes a verdict, the log says what the old verdict was.

Change log

  1. Published on factchecked.science after a full check against the sources. Built from our verification of the paper and of the S0 News episode of October 23, 2026: six claims rated on the site’s scale, four of them the episode’s, one the paper’s own and one the press release’s; the pressure–temperature chart is adapted from the paper’s Figs. 3 and 5a (CC BY 4.0) using the numbers in its text, and the 1928–2026 timeline is our own, from the sources’ dates.

Suggest a correction

Spotted an error, or know a better source? Tell us here. Submissions aren’t published.