REVIEW 2 major objections 2 minor 4 cited by
Oxygen left behind: Atmospheric Enrichment due to Fractionation in Sub-Neptunes using BOREAS
T0 review · 2 major / 2 minor · reviewed 2026-05-16 · grok-4.3
Pith's one-line read Sub-Neptunes can become water-rich planets when hydrogen escapes efficiently while oxygen is mostly retained.
desk verdict BOREAS gives a concrete 200 Myr enrichment timeline for sub-Neptunes but layers fractionation on a single-fluid Parker wind, so the partial oxygen entrainment needs multi-fluid validation. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
BOREAS, the 1D Parker wind model coupled to a mass-dependent fractionation scheme that calculates separate escape rates for hydrogen and oxygen in H2+H2O mixtures.
What would settle it
Detection of a sub-Neptune near the radius valley that retains a thick hydrogen-dominated atmosphere with no measurable water enrichment, or direct measurement of oxygen escape rates far higher than the model predicts for the observed XUV flux.
Extended reading notes
Core claim
BOREAS couples a 1D Parker wind solution to a mass-dependent fractionation scheme and follows hydrogen and oxygen loss rates across planet mass, radius, temperature, and XUV flux. Oxygen remains largely bound except at the highest XUV levels or on the lowest-gravity worlds; elsewhere efficient hydrogen escape with only partial oxygen entrainment produces substantial atmospheric enrichment in roughly 200 Myr, turning sub-Neptunes near the radius valley into water-rich planets.
Load-bearing premise
A one-dimensional Parker wind plus simple mass-dependent fractionation accurately describes the coupled loss of hydrogen and oxygen from mixed atmospheres across the full range of planet masses, radii, temperatures, and XUV fluxes explored.
Editorial extensions
If this is right
- Sub-Neptunes near the radius valley can evolve into water-rich planets over 200 Myr.
- Present-day water-rich atmospheres can form from initially water-poor envelopes.
- Chemical fractionation must be included in long-term atmospheric evolution calculations.
- Oxygen loss is efficient only at high XUV fluxes or on low-gravity planets; elsewhere hydrogen loss dominates.
Reading between the lines
- Atmospheric processing may erase the original envelope composition for many close-in planets, so current water content need not reflect formation location.
- Similar fractionation could operate on other H2-dominated worlds and produce observable differences in C/O or D/H ratios.
- Radius-valley planets may record a common late-stage enrichment phase rather than two distinct formation channels.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces BOREAS, a model coupling a 1D Parker wind formulation with a mass-dependent fractionation scheme, to study hydrodynamic escape and chemical fractionation in mixed H2 + H2O atmospheres of sub-Neptunes. Across ranges of planet mass, radius, equilibrium temperature, and XUV flux, the model predicts that oxygen is largely retained while hydrogen escapes efficiently under intermediate conditions, producing substantial atmospheric enrichment over ~200 Myr and allowing sub-Neptunes near the radius valley to evolve into water-rich planets, consistent with observations of GJ 9827 d.
Significance. If the fractionation implementation is robust, the work supplies a concrete mechanism linking escape physics to the observed radius-valley demographics and the possible origin of water-rich envelopes from initially H2-dominated ones. Public release of the BOREAS code is a clear strength that supports reproducibility and community testing of the enrichment timeline.
major comments (2)
- [Model formulation (BOREAS description)] The central enrichment timeline (~200 Myr with partial O entrainment) rests on the BOREAS fractionation scheme being applied atop the single-fluid, isothermal Parker wind solution. This construction does not automatically satisfy multi-species momentum exchange or energy balance when H2 and H2O possess different scale heights and dissociation states; the transition from O retention to partial loss at intermediate XUV fluxes could therefore be sensitive to the specific drag or efficiency parameterization chosen rather than emerging from the coupled fluid equations.
- [Results and parameter exploration] No direct comparison is shown to multi-fluid hydrodynamic escape calculations or to laboratory/analytic benchmarks for H2-H2O fractionation under XUV-driven flows. Without such validation, the reported oxygen retention efficiency across the explored (M_p, R_p, T_eq, F_XUV) grid remains an untested modeling choice that directly controls the predicted water-rich outcome for radius-valley planets.
minor comments (2)
- [Abstract] The abstract states that 'oxygen is efficiently retained over most of the parameter space' but does not quantify the fraction of the grid or the precise XUV threshold separating retention from loss; adding a brief numerical summary would improve clarity.
- [Figures] Figure captions and axis labels should explicitly state the assumed H2O/H2 mixing ratio and the functional form of the fractionation factor (e.g., mass-ratio scaling or drag coefficient) so readers can reproduce the plotted escape rates.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed review. The comments raise valid points about model assumptions and validation that we address point by point below. We have revised the manuscript to strengthen the discussion of limitations and add supporting comparisons where feasible.
read point-by-point responses
-
Referee: The central enrichment timeline (~200 Myr with partial O entrainment) rests on the BOREAS fractionation scheme being applied atop the single-fluid, isothermal Parker wind solution. This construction does not automatically satisfy multi-species momentum exchange or energy balance when H2 and H2O possess different scale heights and dissociation states; the transition from O retention to partial loss at intermediate XUV fluxes could therefore be sensitive to the specific drag or efficiency parameterization chosen rather than emerging from the coupled fluid equations.
Authors: We acknowledge that BOREAS employs a single-fluid isothermal Parker wind base with a superimposed mass-dependent fractionation scheme rather than solving the full multi-fluid momentum and energy equations. This approximation was selected to permit broad parameter exploration while incorporating the dominant effect of differential escape velocities arising from the distinct molecular masses and scale heights of H2 and H2O. The fractionation factors are computed directly from the Parker solution at each altitude, ensuring consistency with the underlying wind structure. We agree that this approach can introduce sensitivity to the adopted drag and efficiency terms. In the revised manuscript we have expanded the methods and discussion sections to explicitly state these assumptions, quantify the range of drag coefficients explored, and note that full multi-fluid treatments may modify the precise transition fluxes. The core trends—efficient H escape with partial O retention at intermediate XUV—remain robust across the grid and align with the observed properties of GJ 9827 d. revision: partial
-
Referee: No direct comparison is shown to multi-fluid hydrodynamic escape calculations or to laboratory/analytic benchmarks for H2-H2O fractionation under XUV-driven flows. Without such validation, the reported oxygen retention efficiency across the explored (M_p, R_p, T_eq, F_XUV) grid remains an untested modeling choice that directly controls the predicted water-rich outcome for radius-valley planets.
Authors: We agree that explicit benchmarks would strengthen confidence in the quantitative retention efficiencies. The original manuscript relied on the well-established Parker wind framework and fractionation principles validated in prior single-fluid studies. For the revision we have added direct comparisons to analytic fractionation expressions (e.g., those derived from Jeans and hydrodynamic escape theory) and referenced published multi-fluid results for H/He and H2O mixtures that show comparable oxygen retention under similar irradiation levels. A full new suite of multi-fluid simulations lies outside the scope of the present work, but the public BOREAS code release enables such tests by the community. We have also clarified in the text that the reported enrichment timeline is model-dependent within the stated approximations and should be viewed as a prediction to be tested against more complete treatments. revision: partial
Circularity Check
No significant circularity; model outputs emerge from parameter sweeps on standard physics
full rationale
The paper constructs BOREAS by coupling the standard 1D Parker wind solution to a mass-dependent fractionation scheme and then integrates the resulting escape rates over a grid of planet mass, radius, T_eq and XUV flux. The reported ~200 Myr enrichment timeline and the implication that sub-Neptunes near the radius valley can become water-rich are direct numerical outcomes of those integrations, not quantities that were fitted to the same data or defined by the model equations themselves. No load-bearing self-citation, uniqueness theorem, or ansatz is invoked to force the central result; the derivation therefore remains self-contained against external hydrodynamic and fractionation benchmarks.
Assumptions & free parameters
free parameters (2)
- XUV flux scaling
- planet mass and radius
assumptions (2)
- domain assumption The 1D Parker wind formulation remains valid for the mixed H2+H2O composition and the explored XUV heating regime.
- domain assumption Fractionation occurs via mass-dependent drag in the outflow without additional chemical or radiative effects.
Cite this review
Pith. "Pith review of Oxygen left behind: Atmospheric Enrichment due to Fractionation in Sub-Neptunes using BOREAS." pith.science (2026). https://pith.science/paper/2602.12201
@misc{pith2026260212201,
author = {Pith},
title = {Pith review of: Oxygen left behind: Atmospheric Enrichment due to Fractionation in Sub-Neptunes using BOREAS},
year = {2026},
howpublished = {\url{https://pith.science/paper/2602.12201}},
note = {Machine review of arXiv:2602.12201}
}
read the original abstract
The evolution of exoplanetary atmospheres is strongly influenced by atmospheric escape, particularly for close-in planets. Fractionation during atmospheric loss can preferentially remove lighter elements such as hydrogen, while retaining heavier species like oxygen. In this study, we investigate how and under what conditions hydrodynamic escape and chemical fractionation jointly shape the mass and composition of exoplanet atmospheres, especially for mixed H2 + H2O atmospheres. We develop BOREAS, a self-consistent mass loss model coupling a 1D Parker wind formulation with a mass-dependent fractionation scheme, which we apply across a range of planet masses, radii, equilibrium temperatures, and incident XUV fluxes, allowing us to track hydrogen and oxygen escape rates at different snapshots in time. We find that oxygen is efficiently retained over most of the parameter space. Significant oxygen loss occurs under high incident XUV fluxes, while at intermediate fluxes oxygen loss is largely confined to low-gravity planets. Where oxygen is retained, irradiation is too weak to drive significant escape of hydrogen and thus limiting atmospheric enrichment. By contrast, our model predicts that sub-Neptunes undergo substantial atmospheric enrichment over approx. 200 Myr when hydrogen escape is efficient and accompanied by partial oxygen entrainment. Notably, our results imply that sub-Neptunes near the radius valley can evolve into water-rich planets, in agreement with GJ 9827 d. Present-day water-rich atmospheres may have originated from water-poor envelopes under some conditions, highlighting the need to include chemical fractionation in evolution models. BOREAS is publicly available.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We develop BOREAS, a self-consistent mass loss model coupling a 1D Parker wind formulation with a mass-dependent fractionation scheme... diffusion–drag framework of D. M. Hunten et al. (1987) and K. J. Zahnle & J. F. Kasting (1986)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The hydrodynamic solution is obtained by coupling an isothermal Parker wind above R_XUV to a hydrostatic, isothermal layer below R_XUV... energy-limited mass-loss rate MEL = η F_XUV ...
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 4 Pith papers
-
Reflation: redox-driven atmospheric inflation as tracer of super-Earth geochemistry
Reduced-mantle super-Earths can transiently re-inflate late in life, with bulk density drops up to ~60%, as escape-driven pressure loss turns dissolved water into a light H2 atmosphere.
-
KRONOS I: The $1{-}2.8\mu$m JWST Transmission Spectrum of the 23 Myr V1298 Tau c
JWST NIRISS/SOSS data yield a transmission spectrum of the 23 Myr V1298 Tau c showing H2O absorption and an atmospheric metallicity of 14.8^{+56}_{-12} times solar, lower than mature planets of similar mass.
-
Coupled Thermal-Chemical Evolution Models of Sub-Neptunes Reveal Atmospheric Signatures of Their Formation Location
Coupled thermal-chemical models indicate that sub-Neptunes formed outside the water-ice line exhibit high atmospheric CH4, H2O, and C/O ratios while those formed inside show suppressed CH4 and low C/O.
-
The HST/WFC3 Transmission Spectrum of AU Mic b Part I: An Atmosphere Obscured by Contamination and Systematics
The transmission spectrum of AU Mic b is dominated by the transit light source effect from stellar spots, yielding only weak atmospheric constraints with a preferred scale height below 185 km.
Reviewed May 16, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.