Phosphate availability and implications for life on ocean worlds

The search for life in the solar system focuses on the presence of liquid water, availability of biogenic elements, and chemical energy for metabolism. Subsurface oceans are proposed for multiple icy moons (e.g., Europa, Enceladus) and dwarf planets, with spacecraft data indicating potentially habitable conditions on Europa (to be explored by Europa Clipper) and Enceladus. On Enceladus, Cassini plume measurements suggest ocean conditions within Earth-like temperature and pressure ranges, alkaline pH (>11), salinity 0.05–0.2 mol/kg, hydrothermal activity (>50 °C), and redox disequilibria (CO₂ and H₂) that could support life. For Europa, magnetometer data imply a saline ocean reacted with silicate crust, and surface radiolytic oxidants could supply energy when mixed with reductants if delivered to the ocean. A central question is whether phosphate, essential to all known life and often limiting in Earth’s oceans, could be limiting on ocean worlds. Given that bulk phosphorus abundances in carbonaceous chondrites (~0.01 wt%) and basalts are substantial, the key issue is not total P inventory but the processes governing its partitioning into oceans during water-rock interaction and subsequent sinks and sources over time.

Extensive work on Earth’s phosphorus cycle shows continental weathering as a major source and hydrothermal systems as important sinks over geologic timescales. Extending to ocean worlds, prior modeling suggested that small, ice-covered satellites lack granitic crust and rainfall-driven weathering, and that hydrothermal circulation may be active on Enceladus and Europa, potentially removing phosphate. Lingam and Loeb argued that alkaline oceans (e.g., Enceladus) may have low phosphate release by primary mineral dissolution, while more acidic conditions (postulated for Europa) could enhance phosphate via higher dissolution rates. Measurements from terrestrial hydrothermal vents (50–350 °C) show dissolved phosphate of 2.3–3.3 × 10⁻⁶ mol/kg, indicating that hydrothermal fluids can contain biologically useful phosphate and may be sources relative to phosphate-poor oceans. Hao et al. inferred abundant phosphate for Enceladus under equilibrium with bulk ocean conditions. Bulk P in chondrites is consistently near 0.01 wt%, and basaltic compositions are only slightly lower, reinforcing that availability is governed by partitioning processes (release and sequestration) rather than bulk inventory.

Equilibrium geochemical modeling was performed using EQ3/6 to simulate water-rock reactions across a wide range of compositions and conditions representative of ocean worlds. The initial aqueous phase was a “minimal fluid” (neutral, anoxic water) seeded with low elemental inventories for numerical stability; equilibrium speciation was computed at 50 MPa over 0–300 °C in 1 °C steps. Solid reactants were parameterized as a 1.00 kg “special reactant” with elemental abundances based on different carbonaceous chondrite classes (Al, Ca, Cr/Fe, K, Li, Mg, Mn, Na, O, P, S, Si; per Supplemental Table 1), and redox sensitivity analyses adjusted oxygen fractions to represent different bulk oxidation states while holding other elemental mass fractions constant. Reactions were simulated over a broad range of water/rock (W:R) ratios via the EQ6 progress variable, yielding thousands of W:R conditions at each temperature up to 50 MPa, representative of hydrothermal aquifers on Enceladus and Europa. Orthophosphate species (H2PO4−, HPO4²−, PO4³−) and relevant mineral phases were included; dissolved phosphate was defined as the sum of orthophosphate species, which accounted for 99.9% of dissolved P in all simulations. Multiple chondrite classes were evaluated, with CI chondrites emphasized for detailed analysis due to their close match to solar stoichiometry and higher modeled phosphate yields. Preliminary calculations were also conducted for MORB-like compositions and for magmatic processing scenarios including CO₂ volatiles. Biological implications were assessed by estimating upper bounds on cell abundance from modeled dissolved phosphate divided by mean per-cell P content from 336 aquatic cells (mean ~0.69 ± 0.38 fg P per cell), and by computing hypothetical P-limited doubling times using phosphate uptake kinetics parameterized for the picoplankton Prochlorococcus (following Lomas et al.), representing conservative uptake rates for phosphate-limited microbial populations.

• Modeled equilibrium dissolved phosphate from water-rock reaction spans a wide range across compositions and conditions: reported as 10⁻¹¹ to 10⁻¹⁴ mol/kg for the full parameter space considered, but generally exceeds 10⁻⁵ mol/kg for most plausible scenarios.
• Orthophosphate species constitute 99.9% of dissolved P in all simulations.
• Terrestrial hydrothermal fluids (50–350 °C) contain 2.3–3.3 × 10⁻⁶ mol/kg dissolved phosphate; despite being sinks relative to deep-ocean P-enriched waters, such fluids can be sources in phosphorus-poor oceans.
• For CI chondrite compositions and reaction conditions consistent with Enceladus and Europa constraints, modeled dissolved phosphate is elevated: for Enceladus, NaCl observations imply W:R ~0.1–6, within which modeled equilibrium phosphate is high (e.g., often >10⁻² mol/kg for nominal CI at W:R ~1); for Europa, salinity constraints implying W:R < 40 correspond to phosphate concentrations >10⁻¹ mol/kg for CI compositions (noting parameter-space variability and that highest temperatures yield lower phosphate).
• High-temperature (>200 °C) reactions tend to yield lower phosphate, but available constraints suggest ocean chemistry is unlikely dominated by such high-temperature inputs; inferred hydrothermal temperatures for Enceladus and Europa are typically ≤150 °C.
• Across temperatures <150 °C and a wide W:R range, total dissolved phosphate commonly exceeds 10⁻² mol/kg for nominal CI compositions, with concentrations consistently >10⁻¹ mol/kg for W:R < 100.
• Biological implications: using mean cellular P quotas, modeled phosphate supports potential cell densities of ~10⁶–10⁹ cells/m³, comparable to or greater than Earth’s deep ocean (10⁴–10⁶ cells/m³) under many modeled conditions. P-limited doubling times derived from conservative uptake kinetics are typically days to weeks across broad temperature and W:R ranges.
• Overall, phosphate availability is unlikely to limit the establishment of detectable populations of Earth-like cells on ocean worlds under most plausible conditions.

The study addresses whether phosphate scarcity could limit life on ocean worlds by modeling equilibrium phosphate generated during water-rock interaction across realistic compositions (multiple carbonaceous chondrites, MORB) and conditions (0–300 °C, broad W:R, ≤50 MPa, varying redox states). Results indicate that equilibrium with alteration mineral assemblages yields substantial dissolved phosphate over a large swath of parameter space, particularly at temperatures ≤150 °C and moderate-to-low W:R ratios inferred for Enceladus and Europa. While some processes (e.g., high-temperature reactions, unusual chondrite compositions) can depress phosphate to low values, ocean chemistry is unlikely to be dominated by such extremes. Observational constraints (e.g., NaCl in E-ring grains; hydrothermal nanoparticle inferences; Europa salinity bounds) are consistent with W:R and temperature regimes that produce abundant phosphate. Comparison with terrestrial microbial requirements and kinetics suggests that model concentrations would permit rapid growth and large biomass, implying phosphate abundance is not a primary barrier to habitability. The work reframes “sinks” (e.g., hydrothermal uptake) as equilibrium partitions that set nonzero dissolved concentrations, rather than driving phosphate to extinction, and highlights that dynamic processes mixing diverse reaction environments likely further elevate bulk ocean phosphate above minimum equilibrium values.

Equilibrium geochemical models of water-rock reactions indicate that most plausible ocean world scenarios yield dissolved phosphate concentrations sufficient to support Earth-like microbial growth. Phosphate is therefore unlikely to be the limiting nutrient controlling habitability on small icy ocean worlds such as Enceladus, and potentially Europa, given their expected compositions and hydrothermal conditions. The findings emphasize equilibrium partitioning with alteration assemblages, the importance of reaction conditions (temperature, W:R, redox), and the likelihood that dynamic mixing of diverse sources elevates bulk ocean phosphate. Future work should refine predictions for larger, differentiated bodies where magmatic processing and high-pressure mineralogy alter phosphate equilibria, incorporate additional volatile species and mineral precipitation pathways (e.g., gypsum affecting Ca–apatite equilibria), and explore feedbacks between biology and geochemistry that could raise or lower steady-state phosphate relative to abiotic baselines.

• Applicability is most direct to small bodies with chondritic, undifferentiated rocky interiors; extrapolation to larger differentiated worlds (e.g., Europa’s deeper crust, Ganymede, Titan) is limited, especially at high pressures and with different mineralogies.
• Models assume equilibrium with alteration mineral assemblages using EQ3/6; kinetic effects, transient phases, and transport limitations are not explicitly resolved.
• Pressure was limited to ≤50 MPa; higher-pressure regimes could change mineral equilibria relevant to phosphate partitioning.
• Initial aqueous composition used a simplified “minimal fluid”; real oceans may have more complex chemistries affecting phosphate speciation and mineral saturation.
• Biological implications are illustrative only: cell quotas and uptake kinetics are based on Earth organisms (e.g., Prochlorococcus) and do not predict alien life; energy and other nutrient limitations are not co-modeled.
• Some compositional scenarios (e.g., unusual chondrite classes, high-temperature reactions) yield lower phosphate; the global ocean’s integration of diverse sources is not explicitly time-resolved.
• Origin-of-life chemistry requirements are not addressed.