Empirical Cosmic Shoreline: Atmospheric Threshold

• Empirical Cosmic Shoreline (ECS) is a data-driven boundary defining the transition between planets that can retain substantial atmospheres and those subject to atmospheric loss.
• It utilizes a power-law relation between planetary escape velocity and XUV irradiation, calibrated using Solar System benchmarks and exoplanet data.
• The ECS framework guides exoplanet target selection for missions like JWST and ELT while informing population-level studies of atmospheric evolution.

The Empirical Cosmic Shoreline (ECS) is a data-driven, power-law boundary in planet parameter space that empirically demarcates the transition between worlds expected to retain substantial atmospheres and those susceptible to atmospheric loss. Originally proposed in the context of the Solar System and subsequently generalized to exoplanet populations, the ECS encapsulates the cumulative effect of high-energy stellar irradiation (notably XUV flux) and planetary gravity (as parameterized by escape velocity) on the fate of planetary atmospheres. This construct is widely used for interpreting population-level trends in exoplanet atmospheric occurrence, prioritizing JWST and ELT targets, and evaluating theoretical models of hydrodynamic escape and volatile evolution.

1. Mathematical Formulation and Evolution of the ECS

The canonical ECS was introduced by Zahnle & Catling (2017) as a simple power law relating insolation $F$ or time-integrated XUV fluence $F_{\rm XUV}$ to planetary escape velocity $v_{\rm esc}$ through

$F \propto v_{\rm esc}^\alpha$

with an empirical exponent $\alpha \sim 4$ based on Solar System bodies. In normalized variables, this reads

$$\frac{F}{F_\oplus} = c\,\left(\frac{v_{\rm esc}}{v_{\rm esc,\oplus}}\right)^\alpha$$

where $F_\oplus = 1361$ W m$^{-2}$ and $v_{\rm esc,\oplus} = 11.2$ km s$^{-1}$ set the Earth reference point. More recent analyses, employing expanded exoplanet samples and atmospheric detections from the ExoAtmospheres database, yield a steeper slope. The 2024 ECS best-fit anchored to both Mars and 55 Cnc e is

$F_{\rm XUV}$0

or

$F_{\rm XUV}$1

(Meni-Gallardo et al., 18 Aug 2025). A 3D generalization includes explicit stellar luminosity dependence: $F_{\rm XUV}$2 with $F_{\rm XUV}$3 and $F_{\rm XUV}$4, calibrated on Solar System and exoplanet constraints (Berta-Thompson et al., 2 Jul 2025).

2. Physical Basis: Stellar Irradiation, Escape Velocity, and Atmospheric Loss

The ECS encapsulates the interplay between gravitational binding energy and loss drivers, mainly stellar XUV/EUV irradiation and impact erosion. High-energy photons drive hydrodynamic escape, with the critical threshold determined by the planet's $F_{\rm XUV}$5 and its volatile inventory. Under energy-limited escape,

$F_{\rm XUV}$6

(where $F_{\rm XUV}$7 is an efficiency factor), the critical XUV fluence for significant mass loss scales with $F_{\rm XUV}$8 ($F_{\rm XUV}$9, depending on atmospheric composition and cooling efficiency). Hydrodynamic models reveal non-linearities due to molecular cooling (especially in secondary N$v_{\rm esc}$0-CO$v_{\rm esc}$1 atmospheres), leading to a transition from energy-limited (low $v_{\rm esc}$2) to cooling-limited (high $v_{\rm esc}$3) regimes and a possible "knee" in the ECS (Chatterjee et al., 2024). Additional loss channels such as impact erosion (with $v_{\rm esc}$4 marking significant losses) further modulate atmospheric survival (Zahnle et al., 2017).

3. ECS Derivation and Empirical Calibration

Modern ECS studies leverage ensembles of atmospherically characterized exoplanets together with Solar System bodies. The current empirical ECS is typically defined by a Support Vector Machine (SVM) or logistic regression in $v_{\rm esc}$5 space, anchored on Solar System benchmarks (Mars) and high-irradiation exoplanets with confirmed atmospheres (e.g., 55 Cnc e) (Meni-Gallardo et al., 18 Aug 2025, Radica, 31 Mar 2026). When extrapolated to M-dwarf planets, the ECS is severely affected by uncertainties in historic XUV evolution. Incorporating corrections for prolonged activity and pre-main-sequence overluminosity (e.g., Pass et al. 2025), the ECS boundary for mid-to-late M-dwarfs shifts so that only the most massive rocky planets (e.g., LHS 1140 b, TOI-715 b) are expected to retain secondary atmospheres (Pass et al., 1 Apr 2025, Robertson et al., 15 May 2026).

Empirical ECS expressions for various contexts:

Context	ECS Expression	Paper
Solar System (original)	$v_{\rm esc}$6	(Zahnle et al., 2017)
Exoplanets (empirical ECS)	$v_{\rm esc}$7	(Meni-Gallardo et al., 18 Aug 2025)
M dwarfs (Galactic+exgal)	$v_{\rm esc}$8	(Radica, 31 Mar 2026)
3D ECS (with $v_{\rm esc}$9)	$F \propto v_{\rm esc}^\alpha$0	(Berta-Thompson et al., 2 Jul 2025)

Atmosphere retention is commonly quantified by an Atmosphere Retention Metric (ARM): $F \propto v_{\rm esc}^\alpha$1 or, for specific parameterizations and scalings (e.g., Pass et al., Radica et al.): $F \propto v_{\rm esc}^\alpha$2 with ARM = 0 corresponding to the shoreline; ARM > 0 implies atmosphere retention and ARM < 0 predicts loss (Pass et al., 1 Apr 2025, Robertson et al., 15 May 2026, Radica, 31 Mar 2026).

4. Modeling Uncertainties and Physical Regime Transitions

Key sources of ECS uncertainty and transition broadening include:

• Variation in initial volatile inventory ($F \propto v_{\rm esc}^\alpha$3–$F \propto v_{\rm esc}^\alpha$4).
• Time-dependent stellar XUV output: activity lifetimes, flare statistics, pre-main-sequence luminosity.
• Planetary age, composition (H$F \propto v_{\rm esc}^\alpha$5, CO$F \propto v_{\rm esc}^\alpha$6, CH$F \propto v_{\rm esc}^\alpha$7, N$F \propto v_{\rm esc}^\alpha$8), and thermal/cooling properties.
• Non-linear escape physics, such as line-cooling–limited mass loss at high irradiation (collisional-radiative regime), which can flatten or "knee" the $F \propto v_{\rm esc}^\alpha$9–$\alpha \sim 4$0 boundary for super-Earths (Chatterjee et al., 2024, Ji et al., 28 Apr 2025).
• Observational ambiguities; e.g., Venus-like high-altitude aerosol clouds can mimic the ECS signature in transmission spectra, necessitating emission/phase-curve diagnostics for robust atmosphere loss inferences (Lustig-Yaeger et al., 2019).

These factors broaden the ECS to an extended transition region, particularly relevant near the habitable zone thresholds and for planets "straddling" the shoreline, as demonstrated for GJ 3378 b ($\alpha \sim 4$1; atmospheric fate unresolved (Robertson et al., 15 May 2026)).

5. Application to Exoplanet Target Selection and Implications for Population Studies

The ECS constitutes an essential framework for prioritizing exoplanet atmosphere searches. Key applications include:

• Identification of JWST and ELT targets: ECS position and ARM/instellation distance ($\alpha \sim 4$2) metrics optimize searches for retained vs. airless rocky planets (Meni-Gallardo et al., 18 Aug 2025, Ji et al., 28 Apr 2025).
• Sample curation: For low-mass M-dwarfs ($\alpha \sim 4$3), only a handful of sub-$\alpha \sim 4$4 planets (e.g., TOI-1452 b, TOI-715 b) are predicted to reside securely in the retention zone, implying that most similar worlds will be airless (Meni-Gallardo et al., 18 Aug 2025).
• Population-level diagnostics: ECS predictions are statistically tested using density–instellation trends in large samples ($\alpha \sim 4$5 mass/radius measurements), with high-density planets lying above the shoreline and low-density, volatiles-bearing planets below (Ji et al., 28 Apr 2025).
• M-dwarf habitability prospects: ECS slopes and intercepts for M-dwarfs are as steep or steeper than for Sun-like stars, but the threshold for atmospheric loss occurs at much lower $\alpha \sim 4$6 due to persistent high-energy output, constraining the habitable planet inventory (Radica, 31 Mar 2026, Berta-Thompson et al., 2 Jul 2025).

6. Limitations, Degeneracies, and Future Directions

The ECS, while empirically successful, is subject to several caveats:

• The assumption of a universal power law is an oversimplification; regime transitions, "knees," or multi-dimensional thresholds are likely given atmospheric composition and mass–radius dependencies (Chatterjee et al., 2024, Berta-Thompson et al., 2 Jul 2025).
• Systematic uncertainties in stellar XUV histories, atmospheric initial conditions, and upper-atmosphere cooling physics limit predictive precision.
• Observational degeneracies between cloudy and thin (stripped) atmospheres—especially for Venus analogs—can produce statistical false positives that align with the ECS in parameter space (Lustig-Yaeger et al., 2019).
• The ECS is inherently a population-level construct; individual exceptions can arise due to stochastic impacts, magnetic shielding, or planetary evolutionary histories.

Forthcoming JWST programs (e.g., Rocky Worlds DDT) and future ELT surveys are anticipated to considerably sharpen the ECS via a dramatic expansion in the number of planets subjected to direct atmospheric investigations. This will enable improved calibration of ECS slopes, intercepts, and transition widths, test for physical breaks in scaling, and yield rigorous statistical validation for exoplanet habitability criteria (Berta-Thompson et al., 2 Jul 2025, Meni-Gallardo et al., 18 Aug 2025).

7. Empirical ECS as a Predictive and Interpretive Tool

The ECS serves as a quantitative, observer-anchored guide for atmospheric retention:

• Directs observing time toward high-probability retained worlds and interprets null detections as confirmation of XUV-controlled loss.
• Connects robustly with theoretical hydrodynamic escape and atmospheric chemistry models (energy-limited to cooling-limited regimes).
• Provides a unifying principle, buttressed by Solar System and exoplanet population evidence, for interpreting the volatile status and evolutionary fate of rocky exoplanets across stellar environments.

Ongoing advances in high-precision stellar characterization, time-resolved X-ray/UV monitoring, and mass–radius measurement accuracy are critical for progressively refining what constitutes the true empirical shoreline. As a result, the ECS remains a central organizing framework in comparative planetology and the search for habitable exoplanets.