Cosmic dust fertilization of glacial prebiotic chemistry on early Earth

The study addresses how concentrated feedstocks of bioessential elements (P, S, N, C) could have been generated on early Earth to enable high-yield prebiotic chemistry. Common terrestrial rocks are relatively poor in reactive and soluble forms of these elements, and pre-enzymatic systems lacked biological mechanisms to access them efficiently. The authors hypothesize that fine-grained cosmic dust, continuously accreted and subsequently concentrated by sedimentary processes, could have supplied and focused these limiting elements at Earth's surface. They propose that arid environments, particularly glacial settings forming cryoconite sediments, provided efficient traps that could concentrate exogenous dust and its bioessential inventory to prebiotically relevant levels. To test this, they combine astrophysical simulations of early Solar System dust fluxes and sources with geological models of terrestrial sedimentation and concentration.

Prior work has highlighted the importance of moderate-to-high concentrations of simple species (for example, HCN, phosphate, bisulfite) for prebiotic synthesis of nucleic acids, lipids, and peptides, yet the geological means to generate such concentrations remain debated. Terrestrial sources of reactive P and S are scarce, and Earth's surface reservoirs of bioavailable P, S, N, and C are limited. Large impactors can deliver organics and minerals but do so episodically and destructively, with substantial thermal loss of volatiles. In contrast, cosmic dust provides a continuous, finer-grained flux, with some particles traversing the atmosphere gently enough to retain primitive CHNS. Historically, the prebiotic role of cosmic dust was discounted due to dilute global delivery, but field observations show that aeolian, fluvial, and glaciogenic processes can locally concentrate dust by orders of magnitude, including modern cryoconite deposits enriched in micrometeorites. However, quantitative models linking early Earth dust fluxes, sedimentary concentration, and resulting chemical inventories were lacking, motivating the present study.

The authors simulate dust production, transport, and accretion to Earth during the first 500 Myr after the Moon-forming impact (time zero at 50 Myr after CAIs). They model two compositional parent body groups: Jupiter-family comets (JFCs) and asteroids, the latter split into (i) main belt asteroids and (ii) dynamically unstable, rapidly depleted planetesimals left over from terrestrial planet formation. A kinetic collisional-evolution model tracks particle populations in size, pericentre, and eccentricity bins, including destructive collisions, Poynting–Robertson drag, and radiation pressure. For comets, an updated model of comet fragmentation is driven by elevated early scattering rates inferred from N-body simulations. For asteroids, initial masses and depletion histories are set to reflect a more massive early belt and a rapidly decaying leftover planetesimal population, constrained to deliver a Late Veneer equivalent mass of highly siderophile elements. The overlap of dust grain orbits with Earth's orbit is used to compute time-resolved accretion fluxes from each source. To evaluate the fraction of dust that survives atmospheric entry, melts, or vaporizes, they anchor to empirical constraints (e.g., Love & Brownlee) and observed modern size-frequency distributions, applying survival fractions as functions of grain size and entry conditions. They assume an early atmosphere of similar mass to today and consider endmember oxidizing vs reducing conditions for volatile survival, especially nitrogen. Sedimentary concentration is modeled with an all-else-equal scaling: the deposition flux of cosmic dust in a given environment scales linearly with the top-of-atmosphere accretion flux relative to modern. They consider three end-member environments—glacier surfaces (cryoconite in ablation zones), hot deserts, and deep-sea sediments—using empirically derived modern scaling factors for local dust concentration relative to endogenous sedimentation. The cosmic dust mass fraction in sediments (f_dust) is computed as a function of early-to-modern flux ratios and modern flux baselines (evaluated for 10 and 100 t d−1). Compositionally, resulting sediments are treated as mixtures of cometary and asteroidal dust (with entry heating losses applied) and terrestrial endmember sediment (modern upper continental crust for reference), using literature-reported average concentrations of C, N, P, and S for each dust class and terrestrial sediment. They assume concentration mechanisms do not saturate at the elevated early fluxes. Chemical outputs for sediments (C, N, P, S) are then calculated by mass balance across the contributing components with survival/melting/vaporization fractions applied. Sensitivity considerations include uncertainties in modern flux, dust source proportions, entry heating, and atmospheric redox for N retention.

• Early Earth accretion fluxes of cosmic dust were predicted to be ~100–10,000 times higher than modern, due to a more massive early main belt, giant planet migration-driven perturbations and comet injections, and collisional erosion of leftover planetesimals.
• Fluxes were dominated by dynamically unstable asteroidal fragments and comets, with episodic million-year-scale spikes in cometary material; at times >65% of accreted dust was cometary and volatile rich.
• Size-frequency distributions are similar to modern, bimodal with peaks near ~1–10 µm and ~100–500 µm; most particles <0.1 mm survive entry with limited heating, preserving organics and volatiles except for N under oxidizing conditions.
• Predicted cosmic dust proportions in sediments depend strongly on environment: deep-sea sediments remain dust poor even at peak fluxes, whereas deserts and especially glacial ablation zones (cryoconite) can exceed 50% cosmic dust by mass; cryoconite can surpass 80% under high-flux conditions. Post-breakup transient fluxes (~10× background, up to ~10,000× modern) would yield extreme dust-rich compositions.
• Dust-rich cryoconite sediments become highly enriched in bioessential elements relative to average upper continental crust, up to ~100× for P, S, N, and C, with S, N, and C peaking during comet-dominated intervals.
• Mineralogical inventory delivered includes fine-grained reactive P and S phases (e.g., schreibersite, apatite/merrillite, troilite, pyrrhotite, pentlandite, pyrite, chalcopyrite) facilitating rapid dissolution and availability for prebiotic reactions; UV plus sulfide chemistry can oxidize phosphide to phosphate in situ.
• Cosmic dust can stockpile largely inert, refractory C- and N-bearing organics (kerogen-like materials, N-heterocycles, amino acids, purines, etc.). Under reducing anoxic entry conditions, cryoconite blankets could accumulate ~0.2 wt% N and ~5 wt% C, with potential for much higher enrichment via transport fractionation of organics.
• Glacial systems (cryoconite holes and endorheic proglacial lakes) provide semi-closed, low water-rock ratio, freeze–thaw cycling, and UV/UV-shielding niches conducive to prebiotic chemistry, with nutrient recycling documented in Antarctic analogs.

The findings indicate that continuous, largely non-destructive accretion of volatile-rich cosmic dust, when focused by sedimentary processes, could have created localized early Earth environments with abundant, bioavailable P, S, and reservoirs of C and N, overcoming the scarcity in typical crustal materials. Glacial cryoconite systems emerge as the most effective terrestrial concentrators, producing sediments that are both dust-rich and chemically enriched. Such settings would have offered advantageous physical and chemical conditions: cyclical freeze–thaw and wet–dry processes, low dilution, potential UV-driven chemistry with selective shielding, and connectivity among multiple dust-rich microreactors. The study acknowledges that many modern ablation-zone cryoconite holes are transient and hydraulically connected, potentially diluting leached species. However, Antarctic dry valley analogs show ice-lidded, relatively closed cryoconite holes with strong nutrient enrichment and limited lateral transport, suggesting early Earth counterparts could maintain concentrated solutions. Even if cryoconite is transient, sediments and solutes can be delivered into endorheic proglacial lakes that act as longer-lived chemical repositories, enabling sustained prebiotic processes. Chemically, the availability of fine-grained reactive P and S from dust (including phosphide and sulfides) provides solutions to long-standing source problems for prebiotic phosphate and sulfidic anions. Concurrently, dust acts as a stockpiler of refractory C- and N-bearing organics which, upon thermal processing (e.g., impacts or magmatic remobilization), could release HCN and related species central to cyanosulfidic pathways. These mechanisms are synergistic: glacial meltwaters enriched in P and S could mix with HCN outgassed from thermally processed dust-rich sediments. The scenario’s relevance is conditioned by early Earth climate: extensive glaciation or Antarctic-like ice sheets would favor cryoconite formation. Models suggesting a cold, possibly alkaline Hadean surface environment, and geological evidence for early glaciation, are compatible but not definitive. The framework also extends to exoplanets, where observable dust accretion rates may predict the potential for dust-fertilized prebiotic niches.

This work integrates astrophysical dust-flux simulations with geological concentration models to show that early Earth could have hosted localized, cosmic-dust-rich sediments—especially cryoconite on glacier surfaces—that were highly enriched in bioessential elements (P, S, N, C). Early dust fluxes were orders of magnitude above modern, volatile rich, and punctuated by cometary spikes, enabling sediments to become >50–80% dust by mass and up to ~100× enriched relative to upper continental crust in key elements. These environments present compelling settings for prebiotic chemistry, offering both reactive mineral feedstocks and stockpiled organics, and favorable physical dynamics (freeze–thaw, semi-closed aqueous microreactors, UV regimes). Future research directions include: constraining early Earth’s glacial extent and climate to test environmental availability; refining modern flux baselines and early dust source proportions; experimentally quantifying leaching/alteration of dust minerals in cryoconite-like conditions; assessing nitrogen survival under plausible early atmospheric redox states; and searching geological archives for signatures of ancient dust-rich deposits. Extension of these models to exoplanetary contexts could link observed debris/dust environments to the potential for prebiotic fertilization.

• Environmental assumptions: an all-else-equal approach sets early Earth environmental parameters to modern values except dust flux; vegetation absence is assumed not to alter key concentration mechanisms in the modeled settings.
• Flux normalization: results depend on the chosen modern top-of-atmosphere dust flux (evaluated at 10 and 100 t d−1), introducing order-of-magnitude uncertainty in predicted sediment dust fractions.
• Atmospheric entry and composition: survival, melting, and vaporization fractions, and volatile loss (especially N), are based on modern empirical constraints; early atmospheric redox and entry velocities/mineralogy could alter these outcomes.
• Sediment model simplifications: concentration mechanisms are assumed unsaturated even at elevated fluxes; leaching during aeolian transport is neglected; terrestrial endmember composition is taken as modern upper continental crust (early crust could differ, generally making enrichments larger, but still uncertain).
• Temporal and stochastic variability: cometary dust spikes and parent body breakup events are stochastic and transient (Myr scale); cryoconite holes are often short-lived and hydrologically connected, risking dilution, although ice-lidded analogs mitigate this.
• Ground truth: direct constraints on extreme post-breakup fluxes and corresponding sediment compositions in deep time are limited; scaling from modern analog deposits introduces methodological uncertainties (~10% for measured sediment dust fractions, but larger from flux modeling assumptions).