Tiny Crystals Reveal When Magma Is Ready to Erupt

University of Queensland researchers have identified a chemical signature in tiny plagioclase crystals that signals when magma is ready to erupt.

Magma crystallization is the process by which molten rock (magma) cools and solidifies into solid minerals (crystals), eventually forming igneous rock.

Researchers analysed thin sections of 28 representative rocks from seven volcanoes plus two ancient super- eruptions using microscopy and laser-based chemical mapping. Similar microcryst signatures appear more widely around the Ring of Fire, pointing to a common triggering process.

Infographic showing plagioclase crystals over an erupting volcano, with text explaining how tiny crystals reveal magma pressure, temperature, and melt composition before eruption.
AI generated by Grok

University of Queensland researchers have identified a chemical signature in tiny plagioclase crystals that signals when magma is ready to erupt.

PhD candidate Daniela Parra- Encalada (School of the Environment) and senior author Associate Professor Teresa Ubide led the work, published in Nature Communications.

They studied plagioclase, a common white mineral that grows in magma a few kilometres below the surface and responds to changes in temperature, water content, and composition, from Andean volcanoes in Chile and Argentina (across a ~230 km stretch) spanning different ages, eruptive styles, and activity levels.

⁠Key findings

  • The smallest, newest crystals (microcrysts, roughly the width of a human hair) that form in the final stages before eruption share a consistent chemical signature across the studied volcanoes and other Pacific Ring of Fire datasets.
  • Modelling links this signature to hot, deeper magma injecting into a shallower reservoir. This destabilises the volcanic plumbing system and can trigger eruption.
  • The team named the geochemical pattern the “plagioclase eruptibility field”.
  • Earlier-formed plagioclase with lower calcium content often records non-eruptible magma that is cooling and accumulating in the crust (relevant for metal enrichment, such as copper and gold, and only rarely involved in super-eruptions).⁠

Researchers analysed thin sections of 28 representative rocks from seven volcanoes plus two ancient super- eruptions using microscopy and laser-based chemical mapping. Similar microcryst signatures appear more widely around the Ring of Fire, pointing to a common triggering process.⁠

This crystal “fingerprint” could help identify volcanoes (including seemingly quiet ones) that merit closer monitoring because their magma systems show signs of eruptibility. Future work will test whether the pattern holds in oceanic volcanoes and explore links to critical-mineral accumulation.⁠

Plagioclase acts like a recorder of a volcano’s recent history because its composition and zoning (growth layers) change with magmatic conditions.

Focusing on the tiniest crystals isolates the final pre- eruption stage more clearly than bulk- rock chemistry.

Black terrazzo surface with colorful recycled glass and stone fragments
A thin section of volcanic rock under the microscope, showing the white plagioclase crystals of varying sizes. This rock is from Cerro Tuzgle, a dormant volcano in northwestern Argentina. Credit: The University of Queensland

Magma crystallization is the process by which molten rock (magma) cools and solidifies into solid minerals (crystals), eventually forming igneous rock.

1. What is magma?

Magma is a hot, mobile mixture of:

  • Liquid silicate melt (the dominant component)
  • Suspended crystals
  • Dissolved gases (mainly water vapor, CO₂, sulfur gases)
  • Sometimes solid rock fragments

Temperatures typically range from ~700 °C (rhyolitic, silicic magma) to >1200 °C (basaltic and mafic magma).

2. How crystallization begins

Crystallization starts when the temperature drops below the liquidus (the temperature at which the first crystals begin to form). As cooling continues, the magma moves toward the solidus (the temperature at which it becomes completely solid).

The process involves two main stages:

  • Nucleation: Tiny clusters of atoms arrange into stable crystal “seeds.”
  • Crystal growth: Atoms attach to these seeds, causing crystals to enlarge.

3. Sequence of crystallization (Bowen’s Reaction Series)

Minerals do not all form at once. They crystallize in a predictable order based on temperature and composition:

Discontinuous series (mafic minerals):
Olivine → Pyroxene → Amphibole → Biotite mica

Continuous series (plagioclase feldspar):
Calcium- rich plagioclase (high anorthite, An-rich) → progressively more sodium -rich plagioclase (lower An)

Later- forming minerals (lower temperature):
Potassium feldspar, muscovite mica, and quartz.

This is why volcanic rocks often contain a mix of early- formed (larger, sometimes resorbed) crystals and late- formed microcrysts.

4. Key factors that control crystallization

  • Cooling rate: Slow cooling (deep underground) produces large crystals (phaneritic texture, e.g., granite). Fast cooling (eruption at the surface) produces small crystals or glass (aphanitic or glassy texture, e.g., basalt or obsidian).
  • Composition: Higher silica content raises viscosity and lowers crystallization temperature.
  • Pressure: Higher pressure generally raises melting points; decompression can trigger crystallization.
  • Water and volatiles: Dissolved water lowers the melting temperature and can suppress or promote certain minerals (e.g., plagioclase stability expands as magma degasses near the surface).
  • Magma mixing and recharge: Injection of hotter, more mafic magma can dissolve existing crystals or cause new ones to grow with different compositions (producing zoning).

5. Crystal zoning and records of history

Crystals often grow in layers (like tree rings). Changes in temperature, pressure, water content, or melt composition produce compositional zoning.
In plagioclase, for example, the anorthite (Ca) content decreases as the melt evolves or increases when hotter mafic magma recharges the system. These zones preserve a detailed timeline of magma conditions, exactly the “crystal records” used in the recent research on eruptibility.

6. Final outcomes

  • If magma fully crystallizes underground → intrusive (plutonic) rock.
  • If it erupts before full crystallization → volcanic (extrusive) rock containing a mixture of phenocrysts (early crystals), microcrysts (late crystals formed during ascent), and groundmass (fine crystals or glass).

In short, magma crystallization is a progressive, temperature- and composition- driven process that turns liquid melt into a solid crystalline rock, while recording the physical and chemical history of the magma system in the crystals themselves.

Crystal records of erupted vs. trapped magma in volcanic arcs

Arc andesites form by mixing of mafic and silicic magmas. The silicic end-member is hard to constrain because it is underrepresented in erupted products.

Using textural, elemental, and isotopic records in plagioclase, the authors distinguish trapped crystal mush from erupted magma. They define the plagioclase eruptibility field, a narrow, recurrent anorthite range (An ≈ 61), that marks the melt that actually erupts across Central Andes volcanoes after mafic recharge.

Recycled plagioclase cargo is dominated by lower- An compositions (An ≈ 43), recording prolonged crystallisation and assimilation in silicic mushes.

The same pattern appears in continental and composite arcs globally, indicating strong eruptive filtering: most felsic and mafic mushes remain dormant and contribute to crustal growth, though felsic mushes can occasionally remobilise to produce super- eruptions.

Key findings

  • Samples: 28 fresh lava and tephra samples from a ~230 km transect in the Central Andes (main-arc and back-arc settings). Includes polygenetic stratovolcanoes (Lascar, Socompa, Tuzgle), monogenetic centres (El Negrillar, Cerro Overo, Back-Arc Monogenetic Volcanism fields), and older flare-up ignimbrites.
  • Plagioclase populations:
    • Microcrysts (small crystals formed from the final erupted melt during ascent/degassing) and phenocryst rims cluster tightly around An 61 ± 7 (labradorite). This is the “plagioclase eruptibility field.”
    • These crystals also show elevated Mg and Fe, consistent with rapid growth from relatively mafic, recharged melt.
    • The broader crystal cargo (cores, mantles, mush fragments, lithics) is dominated by lower-An plagioclase (average ~An 43- 45), recording silicic (rhyolitic) mushes that largely remain trapped.
  • Interpretation: Mafic recharge into shallow reservoirs homogenises and mobilises intermediate (andesitic) melt that can erupt. The eruptibility field is remarkably consistent across volcano ages, styles (effusive/explosive), and tectonic settings in the study area and appears in other arc datasets.
  • Trapped vs. erupted: Most of the crystal cargo records non-eruptible (or rarely eruptible) felsic mush that contributes to crustal growth and can host metal enrichment. Only under special conditions (e.g., flare-ups) do these mushes remobilise for large eruptions.
  • Supporting data include EPMA, LA -Q- ICP- MS elemental mapping, Sr- isotope analyses, and thermodynamic modelling of equilibrium melts.

Broader implications

The plagioclase eruptibility field offers a petrographic and geochemical indicator of magma that is primed for eruption. It helps explain the prevalence of andesites in arcs (selective filtering of intermediate compositions) and highlights the large volume of unerupted mush.

Future work may test the pattern in oceanic arcs and explore links to critical-metal accumulation in trapped mushes.

The full paper (open access) includes detailed figures of crystal textures, compositional distributions, elemental maps, and global comparisons. A PDF download is available on the Nature page.

Journal: Nature Communications, volume 17, Article number: 11203 (2026)

DOI: 10.1038/s41467-026-78932-9

Provided: University of Queensland

Published: 09 November 2026

Authors: Daniela Parra-Encalada,
Teresa Ubide,
Gideon Rosenbaum,
Yujin Jegal,
Patricia Larrea,
Emilce Bustos &
Catharine Wang

Abstract

Arc andesites are hybrids generated by mixing of mafic and silicic magmas. The silicic endmember can contribute to crustal growth or trigger explosive eruptions, yet remains difficult to constrain due to limited representation in the erupted record. Here, we show that textural, elemental and isotopic records preserved in plagioclase can distinguish between trapped and erupted crystal mush, thus constraining conditions of magma storage vs. eruption in arc systems. We define the plagioclase eruptibility field, a narrow and recurrent anorthite range (An~61) that represents erupted melt across Central Andes volcanoes upon mafic recharge. The recycled plagioclase cargo is dominantly An~43, recording prolonged crystallisation and assimilation within silicic mushes. This pattern repeats in continental and composite arcs globally, indicating strong eruptive filtering of arc magmas and widespread occurrence of dormant felsic and mafic mushes that, in most cases, contribute to crustal growth; yet exceptionally, felsic mushes can rapidly remobilise to generate super-eruptions.

Earth’s Mantle Was “Breathing” Oxygen Billion Years Before the Atmosphere

Scientists found evidence that Earth’s mantle was interacting with (and being oxidized by) surface oxygen as early as ~3.4 billion years ago, roughly a billion years before the Great Oxidation Event permanently raised atmospheric oxygen levels around 2.3- 2.2 billion years ago.

Sanukitoid petrogenesis refers to the origin and formation processes of a distinctive suite of intermediate to felsic igneous rocks (mainly diorites, monzodiorites, granodiorites, and related rocks) that are characteristic of the late Archean to early Paleoproterozoic (~3.0- 2.5 Ga, with some earlier examples from ~3.4- 2.9 Ga).

Early subduction involving inorganic oxidation and intermittent “whiffs” of oxygen between 3.4 and 2.7 Ga may have driven an early rise in sub-arc mantle fugacity before the Great Oxidation Event (GOE, ~2.3-2.2 Ga), reaching conditions comparable to the modern Earth by ~2.1 Ga.

A surreal landscape featuring a glowing pool of lava surrounded by dark volcanic rock, with shimmering blue streams of energy and crystalline formations emerging from the ground.
AI generated by Grok

Scientists found evidence that Earth’s mantle was interacting with (and being oxidized by) surface oxygen as early as ~3.4 billion years ago, roughly a billion years before the Great Oxidation Event permanently raised atmospheric oxygen levels around 2.3- 2.2 billion years ago.

This comes from a University of Portsmouth- led international study published in Nature Geoscience (DOI: 10.1038/s41561-026-02119-6; lead author Dr. Hugo Moreira). It builds on earlier 2023 work by the same group linking mantle chemistry changes to the Great Oxidation Event itself.

Key methods and findings

Researchers analyzed tiny “time capsules”: apatite mineral inclusions sealed inside durable zircon crystals from sanukitoid rocks (a distinctive type of magma linked to early subduction- like processes). Samples spanned more than 3 billion years and came from sites including the Pilbara Craton (Australia; ~3.42 and 2.9 Ga), Karelian and Kaapvaal Cratons (Finland and South Africa; ~2.7 Ga), São Francisco Craton (Brazil; ~2.13 Ga), and Scotland (~425 Ma).⁠

  • Zircons were dated precisely (some to exactly 3.4 billion years old) using laser techniques at Portsmouth.
  • Sulfur oxidation states (S⁶⁺, S⁴⁺, S²⁻) in the protected apatite inclusions were measured via sulfur K-edge micro- X- ray absorption near- edge structure (XANES) spectroscopy at the European Synchrotron Radiation Facility (ESRF) in France. These states act as a sensitive proxy for oxygen fugacity (how oxidizing the magma and mantle environment was).⁠

Results show a progressive increase in oxygen fugacity in the sub-arc mantle starting between ~3.4 and 3 billion years ago. The signature remained present across the studied interval and approached modern levels by ~2.1 billion years ago. This indicates that subduction was already recycling oxidized surface material (from early weathering and intermittent “whiffs” of oxygen) deep into the mantle long before permanent atmospheric oxygenation.⁠

The team ruled out alternative explanations (temperature and pressure variations, fractional crystallization, degassing) and noted that even small amounts of oxidized sulfur would have had a strong effect on the more reduced early mantle.

They propose feedback: early oxidative recycling raised sub-arc mantle oxidation → more oxidizing volcanic gases → more felsic crust and intensified surface weathering → further contribution to long- term atmospheric oxygen buildup.

This provides the most continuous timeline yet of how surface oxygen and deep- Earth processes co-evolved. It supports the idea that plate tectonics (or early analogs) connected Earth’s exterior and interior reservoirs far earlier than often assumed, helping create conditions for a habitable planet.

The work is curiosity- driven basic research into Earth’s origins and the emergence of life, rather than an immediate practical application.⁠

In short, the mantle was effectively “breathing” (exchanging and responding to) oxygen via recycling processes a billion years before free oxygen dominated the atmosphere.


Earth’s First ‘Breath’

Sanukitoid petrogenesis refers to the origin and formation processes of a distinctive suite of intermediate to felsic igneous rocks (mainly diorites, monzodiorites, granodiorites, and related rocks) that are characteristic of the late Archean to early Paleoproterozoic (~3.0- 2.5 Ga, with some earlier examples from ~3.4- 2.9 Ga).

They mark a transitional phase in Earth’s crustal evolution, bridging the dominant Archean tonalite- trondhjemite- granodiorite (TTG) magmatism and more modern- style calc- alkaline arc magmatism.

Geochemical Characteristics

Sanukitoids are defined by a dual “mantle and crust” signature:

  • Mantle- compatible elements: High MgO (typically 1.5–9 wt%), Mg# (molecular Mg/ (Mg and Fe) ≈ 45–65 or higher in primitive members), Ni (>100 ppm in many definitions), Cr (>100- 200 ppm).
  • Incompatible elements: Strong enrichment in large- ion lithophile elements (LILE), especially K₂O (1.5–5 wt%), Ba, and Sr (Ba and Sr often >1,000–1,400 ppm); high light rare- earth elements (LREE); relatively high K₂O/Na₂O compared with TTGs; fractionated REE patterns with low HREE and no or minor Eu anomalies.⁠

Classic strict definition (Stern et al., 1989) focused on SiO₂ 55–- 60 wt% with the above traits, but the broader “sanukitoid suite and series” includes more evolved members produced by fractional crystallization. High- Ti and low- Ti varieties are sometimes distinguished.

They resemble high- Mg andesites (sanukites) from the Setouchi volcanic belt in Japan, hence the name.

Petrogenetic Models

The consensus is that sanukitoids form by interaction between mantle peridotite and a recycled crustal component (typically TTG- like or slab- derived felsic melts, possibly with fluids and/or sediments).

Two main (not mutually exclusive) pathways are debated:

  1. Two- stage (metasomatism and later melting) process (most widely invoked)
    • Stage 1: Subduction- related fluids and or silicate melts (from partial melting of subducted oceanic crust and metabasite, often eclogite- facies, ± sediments) metasomatize the mantle wedge or subcontinental lithospheric mantle. This creates enriched domains (e.g., amphibole- or phlogopite- bearing veins).
    • Stage 2: Later thermal perturbation (slab break-off, mantle plume, or continued subduction heating) causes partial melting of this metasomatized peridotite, generating the sanukitoid parental magma.⁠
  2. One- stage (assimilation and reaction) process
    • TTG- like or hydrous eclogite melts directly react with and assimilate mantle peridotite during ascent, producing hybrid melts with the characteristic dual geochemistry.⁠

Recent Ti- isotope studies (δ⁴⁹Ti) of primitive sanukitoids show elevated values requiring a hydrous metabasite (eclogite) melt component formed in equilibrium with residual rutile. This favors fluid- fluxed melting of eclogite in a subduction- like setting and indicates parental magmas were water- rich and relatively oxidized, comparable to modern arc magmas.⁠

Melting conditions are typically high- pressure (>1- 2.5 GPa, corresponding to ~35- 80 and km depth), hydrous, and oxidized. Differentiation often involves amphibole ± plagioclase ± orthopyroxene fractionation, sometimes at lower pressures where a distinct cotectic may control the series.

Geodynamic Significance

Sanukitoids appear mainly near the Archean- Proterozoic transition and are linked to the progressive cooling of Earth and the evolution of subduction:

  • Early Archean: High heat production favored extensive slab melting → pure TTG magmas with limited mantle interaction.
  • Late Archean: Reduced slab melting → more interaction with the mantle wedge → sanukitoids (hybrids).
  • Post ~2.5 Ga: Predominantly fluid- fluxed melting of the mantle wedge (modern- style arcs), with slab melts rarer.⁠

They are therefore key evidence for the onset or maturation of plate- tectonic (or proto- subduction) processes, deep crust- mantle recycling, and the changing nature of continental crust growth.

In the context of the recent Nature Geoscience study on mantle oxidation, sanukitoids serve as excellent monitors of sub- arc mantle redox state because their sources record metasomatism by recycled surface materials.

Modern Analogues and Distribution

Phanerozoic high- Mg diorites, vaugnerites, appinites, and some post-collisional I- type granitoids share similarities and may represent analogues. Archean examples occur worldwide (e.g., Superior Province, Pilbara, Kaapvaal, Karelian, North China, Amazonian/Carajás, Yilgarn, and São Francisco cratons).

In summary, sanukitoid petrogenesis records hybridization between depleted mantle and recycled crustal melts and fluids in a subduction- related setting under hotter Archean conditions.

This dual signature makes them invaluable for tracking early Earth geodynamics, mantle metasomatism, and the co- evolution of crust and mantle.

Close-up of an orange mineral rock on a textured surface with small pebbles and dust.
Ancient Zircons Reveal Oxygen Was Quietly Oxidizing Earth’s Mantle

Progressive oxidative weathering raised sub-arc mantle fugacity on the Archaean Earth

Abstract (summary)

In the modern Earth, sub- arc mantle oxygen fugacity is linked to addition of sulfur species by slab fluids. It has been unclear to what extent volatile recycling controlled redox changes in the early Earth.

Using sulfur K -edge micro- X- ray absorption near- edge structure (μ-XANES) spectroscopy, the authors measured the relative abundances of S⁶⁺, S⁴⁺, and S²⁻ in apatite inclusions hosted in geochemically characterized and precisely dated zircons from sub- arc sanukitoid magmas spanning 3.4 billion years of Earth history.

Key results:

  • The sub- arc mantle shows a marked increase in oxygen fugacity (fO₂) between 3.4 and 3.0 Ga.
  • It is substantially more oxidized than the unaltered ambient upper mantle.
  • This cannot be explained by fractional crystallization, degassing, or intrinsic mantle processes.

The authors propose that progressive oxidative weathering increased the proportion of oxidized sulfur recycled into mantle- wedge sources.

This established a feedback between subduction of surface materials and progressively more oxidized arc magmatism.

Early subduction involving inorganic oxidation and intermittent “whiffs” of oxygen between 3.4 and 2.7 Ga may have driven an early rise in sub-arc mantle fugacity before the Great Oxidation Event (GOE, ~2.3-2.2 Ga), reaching conditions comparable to the modern Earth by ~2.1 Ga.

Methods highlights

  • Samples: Geochemically coherent sanukitoids (metasomatized- mantle- derived magmas) from the Pilbara Craton (Australia; 3.42 and 2.9 Ga), Karelian Craton (Finland; 2.7 Ga), Kaapvaal Craton (South Africa; 2.7 Ga), São Francisco Craton (Brazil; 2.13 Ga), and Northern Highland Terrane (Scotland; 0.425 Ga).
  • Protected apatite inclusions inside zircon crystals acted as unaltered “time capsules”.
  • 337 analyses on 39 apatite inclusions + 213 analyses on 65 matrix apatites, performed at the European Synchrotron Radiation Facility (ESRF, ID21 beamline).
  • S⁶⁺/S_total ratios converted to ΔFMQ (deviation from the fayalite- magnetite- quartz buffer) using an experimental oxybarometer.
  • Supporting U- Pb dating, trace- element, and oxygen- isotope data on the host zircons.

Broader implication

This provides the first continuous record showing that surface oxidative weathering products were already being recycled into the deep Earth via early subduction- like processes more than a billion years before permanent atmospheric oxygenation. The mantle was not a passive spectator; it participated in a feedback loop that helped shape the conditions leading to a habitable, oxygenated Earth.

Journal: Nature Geoscience

DOI: 10.1038/s41561-026-02119-6

Provided: University of Portsmouth

Published: 14 October 2026 (Open Access)

Authors: Hugo Moreira,
Emilie Bruand,
Craig Storey,
James Darling,
Mike Fowler,
Marine Cotte,
Fleurice Parat,
Lorraine Tual &
Bruno Dhuime

Abstract

In the modern Earth, sub-arc mantle oxygen fugacity is intrinsically associated with addition of sulfur species by slab fluids. However, it is unclear to what extent recycling of volatiles controlled redox changes of the sub-arc mantle in the early Earth. Here we use sulfur K-edge micro-X-ray absorption near-edge structure spectroscopy to measure the relative abundances of S6+, S4+ and S2− in apatite inclusions hosted by geochemically characterized and dated zircons from sub-arc sanukitoid magmas that span 3.4 billion years of Earth’s geological history. The sub-arc mantle shows marked increase in oxygen fugacity between 3.4 and 3.0 billion years ago (Ga). It is substantially more oxidized than unaltered ambient upper mantle, a characteristic that cannot be explained by fractionation, degassing or intrinsic mantle processes. We propose that progressive oxidative weathering enhanced the proportion of oxidized sulfur recycled to mantle wedge sources, establishing feedback between subduction of surface materials and incrementally oxidized arc magmatism. Early subduction involving inorganic oxidation and ‘whiffs’ of oxygen between 3.4 and 2.7 Ga may have regulated an early increase in fugacity of the sub-arc mantle before the Great Oxidation Event, matching similar conditions to the modern Earth by 2.1 Ga.

Farmers’ Lime Is Quietly Pulling Carbon From the Sky

A Yale-led study finds that agricultural liming (adding crushed limestone) in the Mississippi River Basin has acted as a net carbon sink, removing more CO₂ than it emits.

Silicate weathering is the natural chemical breakdown of silicate minerals (e.g., those in basalt, olivine, or other igneous rocks) by carbonic acid formed when atmospheric CO₂ dissolves in water.

A tractor spreading fertilizer over a field at sunset, with lush green crops in the background and dramatic clouds in the sky.
AI generated by Grok

A Yale- led study finds that agricultural liming (adding crushed limestone) in the Mississippi River Basin has acted as a net carbon sink, removing more CO₂ than it emits.

The practice, long used by farmers to raise soil pH, reduce acidity, and boost crop yields, covers a vast region: the Mississippi River Basin spans about 41% of the contiguous U.S. and includes roughly 65% of U.S. cropland.

Researchers, led by geochemist Tim Jesper Suhrhoff (Yale Center for Natural Carbon Capture), analyzed more than a century of data (roughly 1900- 2015) on river chemistry, fertilizer use, atmospheric pollution, and liming rates, combined with geochemical modeling.

Conventional accounting (including IPCC guidelines) treats liming as a CO₂ source because the carbon in limestone can be released as CO₂ when it reacts with strong acids in soil.

The new work argues this is incomplete: soil acidity from nitrogen fertilizers and fossil- fuel pollution is the main driver of those emissions. Without liming, the acids would still react with natural soil alkalinity and release CO₂.

In the presence of liming, much of the reaction instead forms bicarbonate that is exported by rivers (including the Mississippi) and can ultimately sequester carbon in the ocean for long periods. Compared with a no- liming scenario, liming in the basin removed an estimated 300- 400 million tonnes of CO₂.

Efficacy has improved since the 1930s as liming rates rose; current removal is estimated at about 75% of the theoretical maximum.

The study (published in Nature) concludes that liming aligns climate benefits with agronomic ones (higher yields, potentially lower nitrous- oxide emissions). It suggests policy and emissions inventories should attribute emissions to acid inputs rather than to liming itself, and that strategies such as adding silicate rock first to neutralize acidity could further enhance carbon removal.

The linked phys.org article appears to cover this same research (the URL slug matches the topic and timing). Related earlier work from the Yale group has also highlighted liming’s dual role in soil management and carbon cycling.

Silicate weathering is the natural chemical breakdown of silicate minerals (e.g., those in basalt, olivine, or other igneous rocks) by carbonic acid formed when atmospheric CO₂ dissolves in water.

The simplified reaction consumes CO₂ and produces bicarbonate ions (HCO₃⁻) plus dissolved cations such as Ca²⁺, Mg²⁺, K⁺, and Si:

CaSiO3+2CO2+H2O→Ca2++2HCO3−+SiO2\mathrm{CaSiO_3 + 2CO_2 + H_2O \to Ca^{2+} + 2HCO_3^- + SiO_2}

(Similar reactions occur for magnesium and other silicates.) The bicarbonate is ultimately transported to rivers and oceans, where the carbon can remain sequestered for thousands to hundreds of thousands of years. Natural silicate weathering already removes roughly 0.5 Gt CO₂ per year globally; enhanced rock weathering (ERW) accelerates it by crushing silicate rocks into fine powder and spreading them on soils (especially agricultural land).

Primary climate benefit: Carbon dioxide removal (CDR)

ERW converts atmospheric CO₂ into stable bicarbonate that is stored in groundwater and the ocean on intergenerational to geological timescales.

Field trials (e.g., U.S. Corn Belt maize-soybean rotations) have measured cumulative CDR potentials of ~10.5 t CO₂ ha⁻¹ over four years at application rates of ~50 t ha⁻¹ yr⁻¹ of basalt.

Modeling studies project that widespread deployment on U.S. farmland could deliver 0.16–0.30 Gt CO₂ yr⁻¹ by 2050 (rising further by 2070), and global estimates often fall in the 0.5–2 Gt CO₂ yr⁻¹ range under ambitious scenarios (after accounting for mining, grinding, and transport emissions).

Unlike carbonate liming (which can release some of its own carbon under acidic conditions), pure silicate minerals contain no carbon of their own, so the alkalinity they generate is more fully atmospheric in origin. Silicates can therefore act as a clearer net sink and are sometimes proposed as a complement or partial replacement for agricultural lime.

Agronomic and soil- health co -benefits

Soil pH regulation: Silicate dissolution neutralizes acidity, raising pH in a manner similar to traditional liming. This improves nutrient availability, reduces aluminum toxicity, and supports healthier microbial communities.

Nutrient release: Weathering supplies plant- available calcium, magnesium, potassium, phosphorus, silicon, and micronutrients. Field studies show elevated concentrations of these elements in soil and crops, partially offsetting the need for synthetic fertilizers.

Crop yield increases: Trials report 12- 16% higher maize and soybean yields, with similar gains for other crops (e.g., sugarcane, miscanthus). Silicon uptake can also strengthen plant tissues, improving resistance to pests, pathogens, and drought.

Reduced nitrous oxide (N₂O) emissions: Better soil conditions and higher pH can lower emissions of this potent greenhouse gas, improving the overall greenhouse- gas balance of ERW by tens of percent in some simulations.

Additional ecosystem effects: Higher silicon in runoff can favor diatoms over harmful algal blooms (reducing eutrophication risk). Downstream alkalinity can help buffer ocean acidification. Some studies also note modest reductions in surface ozone in agricultural regions, which further supports yields.

Relation to the Mississippi River Basin and liming context

In the recent Nature study on agricultural liming, the authors note that silicates are already used for pH management and have attracted strong interest as enhanced- weathering feedstocks. Because silicates do not carry their own carbonate carbon, they avoid some of the accounting complexities of limestone while still providing the same agronomic benefits (pH correction and nutrient release). Modeling and field work in the U.S. Midwest (including the Corn Belt portion of the Mississippi basin) indicate that combining or sequencing silicate and carbonate amendments could optimize both yield and long- term carbon removal.

Practical considerations

Benefits are greatest in warm, wet, acidic soils with active biological communities (roots and microbes accelerate dissolution). Application rates of tens of tonnes per hectare are typical. Life- cycle emissions from grinding and transport must be minimized (renewable energy helps), and potential release of trace metals is monitored; well- chosen basalts generally show low risk in trials. Measurement, reporting, and verification of the actual CDR remain active research areas because of soil heterogeneity and time lags.

In short, silicate weathering (especially when enhanced on farmland) offers durable carbon removal while delivering direct benefits to soil fertility, crop productivity, and other greenhouse- gas reductions, aligning climate action with existing agricultural practices.

Agricultural liming is a carbon sink in the Mississippi River Basin

Application of carbonate minerals (limestone or dolomite) to farmland, agricultural liming, is a long-established practice that counters soil acidification and boosts crop yields. Conventional greenhouse- gas accounting treats it as a net source of agricultural CO₂ emissions.

Using century-scale records of liming rates and anthropogenic acidity inputs across the Mississippi River Basin (MRB), the authors show that liming has instead acted as a net carbon sink over the past century. River alkalinity flux records indicate that roughly 90% of the ideal CO₂- removal potential of the lime applied since 1900 (≈ 0.44 Gt CO₂) has been realized at the catchment scale. A multi- decadal time lag occurs because of soil cation- exchange processes and solute transport. Reactive- transport modelling corroborates this: an initial emissions pulse (from neutralizing existing soil acidity pools) is followed by net CO₂ removal.

Current accounting frameworks use an incomplete counterfactual, they attribute emissions to the lime itself rather than to the anthropogenic strong acids (from fertilizers and fossil- fuel pollution) that drive CO₂ release. When evaluated against the proper counterfactual (acidity inputs without liming), liming is a net sink on decade- to- century timescales. The results imply that optimized soil- pH management can simultaneously cut agricultural greenhouse- gas emissions, raise yields, and improve soil health.

Key scientific points from the paper

Chemical basis of the conventional view
At low soil pH, carbonate dissolution releases CO₂:

CaCO3+2H+→Ca2++CO2+H2O\mathrm{CaCO_3 + 2H^+ \to Ca^{2+} + CO_2 + H_2O}

IPCC default methods essentially treat the carbon in applied lime as emitted; U.S. inventories use a lower emission factor (~50%).

Revised framework
The dominant driver of CO₂ release is the input of strong acids (nitric and sulfuric) from nitrogen fertilizers, manure, soil- organic- nitrogen mineralization, biological nitrogen fixation, and atmospheric deposition of SO₂/NOₓ oxidation products. These acids titrate existing alkalinity (bicarbonate) and force CO₂ out of the system regardless of whether lime is added. Adding lime supplies alkalinity that can form bicarbonate, which is exported downstream and ultimately sequesters atmospheric carbon.

Mississippi River Basin as the study system
The MRB covers ~41% of the contiguous U.S. land area and ~65% of its cropland, receives most U.S. agricultural lime, and has well-documented river chemistry. The authors compiled spatially explicit, century-long time series of acidity inputs (fertilizer/manure oxidation, atmospheric deposition, etc.) and lime additions. Cumulative acidity has exceeded lime- derived alkalinity by ~20 Tmol by 2015, but the excess is concentrated outside intensively limed Corn- Belt croplands.

Evidence for net sink

  • Empirical river alkalinity fluxes.
  • Gridded reactive- transport modelling (SCEPTER) comparing “acidity- only” vs. “acidity and lime” scenarios.
    Both approaches show that, after an early emissions pulse, liming produces a cumulative carbon sink on the order of hundreds of millions of tonnes of CO₂.

Implications
Emissions inventories and climate policy should attribute CO₂ release to acidity inputs rather than to liming. Better -aligned pH management (including possible pre- treatment with silicate rocks) can deliver agronomic benefits while contributing to carbon- dioxide removal.

The paper is open- access; a PDF download link is available on the Nature page. Authors include researchers affiliated with Yale (among others).

Journal information: Nature (2026) 

DOI: 10.1038/s41586-026-11040-2. www.nature.com/articles/s41586-026-11040-2

Provided: Yale University

Authors: Tim Jesper Suhrhoff,
Christopher T. Reinhard,
Yoshiki Kanzaki,
Samuel Shou-En Tsao,
Beck Woollen,
Tom Reershemius,
Samuel Shaheen,
James Saiers,
Shuang Zhang,
Peter A. Raymond &
Noah J. Planavsky

Abstract

Application of carbonate minerals to arable lands, known as agricultural liming, is a long-standing practice for counteracting soil acidification1,2,3. Although liming boosts crop yields4, it is also considered a source of agricultural carbon dioxide (CO2) emissions5. Here we show, using century-scale records of agricultural liming and anthropogenic acidity inputs for the Mississippi River Basin, that agricultural liming has acted as a net carbon sink over the past century. Records of river alkalinity fluxes suggest that approximately 90% of the ideal CO2 removal potential of agricultural lime added since 1900 (approximately 0.44 GtCO2) has been realized at the catchment scale, with a decadal-scale time lag owing to soil cation exchange and solute transport. These results are consistent with reactive transport modelling of soil cation throughput, which indicates that net CO2 removal emerges after an initial emissions pulse associated with neutralization of soil acidity pools. Current accounting frameworks implicitly apply an incomplete counterfactual, attributing CO2 emissions to lime addition rather than to the anthropogenic acidity inputs that drive CO2 release. Evaluated against the counterfactual of anthropogenic acidity inputs, agricultural liming represents a net carbon sink in the Mississippi River Basin on decade-to-century timescales. These results suggest that optimized soil pH management can reduce agricultural greenhouse gas emissions while simultaneously improving crop yields and soil health.