Mineral evolution facilitated Earth's oxidation

The study addresses why and how Earth transitioned from an ancient O2-deficient atmosphere to the modern O2-rich world, focusing on whether mineral evolution contributed causally to oxygenation. Geological evidence documents two rapid, irreversible rises in atmospheric oxygen: the Great Oxidation Event (~2.4–2.3 Ga) and the Neoproterozoic Oxidation Event (0.8–0.54 Ga). While oxygenation clearly reshaped mineral diversity and abundance, the reverse influence—mineral evolution affecting atmospheric O2—has been understudied. Oxygen produced by oxygenic photosynthesis is consumed by aerobic respiration and oxidation of reduced compounds; long-term O2 accumulation requires a leakage in the carbon cycle via burial of organic matter. The persistence of organic matter arises from intrinsic recalcitrance and extrinsic environmental protections, notably physical protection by mineral association. This work investigates how mineral-mediated physical protection could have influenced organic burial and thereby atmospheric O2, proposing that mineral evolution in low-O2 ancient environments increased adsorption capacity, potentially destabilizing the oxygen cycle and facilitating oxygenation.

The paper synthesizes prior work on: (1) Earth’s oxygenation history and geochemical proxies for redox evolution; (2) mechanisms of organic matter preservation, including intrinsic recalcitrance (e.g., humified material, graphitized petrogenic carbon) and extrinsic factors such as dilution in the water column and mineral-associated physical protection via ligand exchange, cation bridging, van der Waals forces, and hydrogen bonding; (3) first-order kinetics of organic matter degradation and mineral adsorption/desorption; and (4) mineral evolution over geologic time, including the roles of iron and clay minerals in organic matter preservation. It highlights that while oxygenation impacts mineralogy is well documented, causal pathways from mineral evolution to atmospheric O2 change are less explored.

The study develops a minimal conceptual model partitioning organic matter into unprotected (y1) and mineral-protected (y2) pools. Dynamics are governed by first-order kinetics: dy1/dt = −k1 y1 − kp y1 + ka y2 and dy2/dt = kp y1 − ka y2. The protected pool decays via desorption followed by rapid oxidation; thus k2 ≈ ka. Dimensionless variables are introduced: κp = kp/k1, κa = ka/k1, time τ = k1 t; oxygen-exposure time τo scales with atmospheric O2 (τo ∝ pO2). The initial protected fraction f = y2(0)/y0 represents burial efficiency at τo = 0. A probabilistic framework is constructed by expressing adsorption/desorption rate constants as exponentials of an overall interaction factor α, itself a sum of independent contributions βi from physical, chemical, and biological variables (temperature, pH, polymer properties, mineral surface area/charge, enzymes). By the Central Limit Theorem, log κp and log κa approach normal distributions; transforming back yields a broad distribution P(κ) ≈ 1/κ away from distribution tails. The expected values of rate constants are computed over bounds set by αmin and αmax, yielding ⟨κp⟩ ≈ [exp(αmax)−exp(αmin)]/(αmax−αmin) and ⟨κa⟩ ≈ [exp(−αmin)−exp(−αmax)]/(αmax−αmin). For heterogeneous systems (αmax ≫ αmin), the ratio ⟨κp⟩/⟨κa⟩ ≈ exp(αmax+αmin) ≈ κp,max, indicating adsorption dominates desorption on average. Stability analysis links burial efficiency y2(τo)/y0 versus τo (and thus pO2) to two parameters: initial protection f and the maximum adsorption rate κp,max. When κp,max < f/(1−f), burial efficiency decreases monotonically with O2 (stabilizing negative feedback). When κp,max > f/(1−f), burial efficiency initially increases with O2 at low pO2 (feedback absent; unstable regime) until surface capacity saturates, after which burial declines with further O2 rise. To evaluate deep-time implications without direct paleo-parameter constraints, a Fermi-estimation approach benchmarks parameters using modern observations: compiled field data of burial efficiency vs oxygen-exposure time yield f ≈ 0.7 and κp ≈ f/(1−f) ≈ 2; characteristic times constrain κ2 = k2/k1 ∈ (10−9, 1) with ⟨κ2⟩ ≈ 0.05. Modern stability requires κp,max < χ ≈ 2; combined constraints give 0.05 < κp,max,Modern < 2. Deep-time scenarios are evaluated relative to these benchmarks, considering how mineral evolution alters κp,max and f.

• Adsorption dominates desorption: A probabilistic analysis with heterogeneous interaction factors yields P(κ) ≈ 1/κ for adsorption and desorption rate constants and implies ⟨κp⟩/⟨κa⟩ ≈ κp,max, indicating adsorption processes dominate average dynamics.
• Stability criterion: The oxygen cycle exhibits two regimes separated by κp,max = f/(1−f). If κp,max < f/(1−f), burial efficiency declines monotonically with increasing pO2 (stabilizing negative feedback). If κp,max > f/(1−f), burial efficiency initially increases with pO2 at low O2 (negative feedback absent; unstable regime), then declines after adsorption capacity saturates.
• Modern benchmarks: From compiled field data, fModern ≈ 0.7, κp,Modern ≈ 2; constraints on protected-pool degradation give ⟨κ2⟩ ≈ 0.05 and 0.05 < κp,max,Modern < 2. Modern conditions reside in the stable regime.
• Deep-time tendencies: Ancient environments likely had lower initial protection (fArchean, fProterozoic < 0.7) due to fewer terrestrial/riverine minerals and shorter residence times of organic matter on land, but higher κp,max due to lower initial surface coverage (κp,max ∝ exp(−f)). This combination makes it easier to enter the unstable regime at low O2.
• Threshold to destabilize: Using modern benchmarks, a ≳40-fold increase in κp,max would suffice to cross into the unstable regime and facilitate O2 accumulation.
• Iron mineral evolution and the GOE: Transforming abundant dissolved Fe(II) (low organic affinity) into ferrihydrite/Fe(III) oxides (high organic affinity; OM:Fe(III) ≳ 4) increases adsorption capacity by ~4×. Because κp,max scales exponentially with capacity, a fourfold capacity increase leads to an approximate e4 (~55-fold) rise in κp,max, exceeding the ~40-fold threshold and plausibly facilitating the GOE.
• Clay mineral evolution and the NOE: Enhanced Neoproterozoic production of weathering-derived clays (e.g., smectite/kaolinite) increased adsorption capacity by ~3–10× (nominally ~4×), implying an e4 (~55-fold) increase in κp,max, sufficient to destabilize the oxygen cycle and promote the NOE.

The findings link microscale organic matter–mineral interactions to macroscale oxygen-cycle dynamics. In low-O2 ancient settings with low initial mineral association (small f) but ample available mineral surface (high κp,max), adsorption can continue despite rising O2, removing the stabilizing negative feedback and allowing O2 to accumulate until adsorption capacity saturates. Mineral evolution—specifically the appearance and proliferation of Fe(III) minerals during the Neoarchean and increased production of high-surface-area clays during the Neoproterozoic—could have increased adsorption capacity sufficiently (by factors that exponentiate κp,max) to push the system into this unstable regime, facilitating the GOE and NOE. The analysis reconciles modern stable behavior (high f, limited κp,max) with ancient potential instability and provides a quantitative threshold (≈40-fold κp,max increase) for tipping behavior. It also addresses apparent contradictions such as the paucity of organic matter in iron formations by invoking diagenetic desorption and alternative fates of released organic matter, and highlights conditions under which ferrihydrite aging may enhance, not reduce, organic stabilization. The framework suggests testable predictions for mineral-associated organic carbon trends through time and invites integration with biological processes.

This work proposes and analyzes a minimal theoretical framework showing how mineral evolution could destabilize Earth’s oxygen cycle and facilitate atmospheric oxygenation. A probabilistic treatment of adsorption/desorption predicts adsorption dominance and a broad (1/κ) distribution of rates; stability analysis identifies a threshold κp,max = f/(1−f) separating stable from unstable regimes. Using modern environmental benchmarks (f ≈ 0.7; 0.05 < κp,max < 2), the study infers that deep-time increases in mineral adsorption capacity—via Neoarchean Fe(III) mineral proliferation and Neoproterozoic clay production—could yield e4-scale increases in κp,max, surpassing a ~40-fold tipping threshold and contributing to the GOE and NOE. Future work should parameterize mineral evolutionary variables explicitly, integrate biological influences, perform laboratory experiments under ancient analog conditions to constrain parameters, and analyze sedimentary records (e.g., via ramped pyrolysis/oxidation) for mineral-associated organic carbon across oxygenation intervals.

The model uses simplified, idealized assumptions (first-order kinetics; independent, identically distributed factors; broad but bounded α) and does not explicitly include variables directly characterizing mineral evolution or biological processes. The predicted P(κ) ≈ 1/κ distribution and stability boundary derive from these simplifications. The estimated ~40-fold κp,max threshold for destabilization is approximate; high-dimensional models could refine it. Benchmarks rely on modern observations extrapolated to deep time; paleo-parameter values remain uncertain. The treatment of diagenetic transformations (e.g., ferrihydrite aging) and their net effect on organic preservation is simplified and context dependent, requiring targeted experiments and field studies.