---
title: 'M Dwarf Cosmic Shoreline: Atmospheric Retention Boundaries'
url: https://www.emergentmind.com/topics/m-dwarf-cosmic-shoreline
type: topic
---

# M Dwarf Cosmic Shoreline: Atmospheric Retention Boundaries

The M dwarf cosmic shoreline is the extension of the cosmic shoreline concept to rocky planets orbiting M stars, especially mid-to-late, fully convective M dwarfs. It is usually formulated in the plane of cumulative stellar X-ray plus extreme-ultraviolet fluence, \(I_{\rm XUV}\), versus planetary escape velocity, \(v_{\rm esc}\), where planets below the boundary are more likely to retain substantial atmospheres and planets above it are more likely to be stripped. Around M dwarfs, the boundary is inferred to shift because these stars sustain elevated high-energy emission for long intervals, so atmospheric survival depends more stringently on a planet’s gravity and irradiation history than in Solar-System-based calibrations [2504.01182, 1702.03386].

## 1. Classical shoreline and physical interpretation

In its classical form, the shoreline is an empirical relation between irradiation and escape speed. Zahnle and Catling found that the Solar System dividing line follows
\[
I \propto v_{\rm esc}^4,
\]
or, in logarithmic form,
\[
\log_{10} I = 4\,\log_{10} v_{\rm esc} + \mathrm{const}.
\]
Bodies with atmospheres lie to the lower right of this relation and air-free worlds to the upper left. In the M-dwarf context, the same diagram is commonly recast using cumulative XUV fluence rather than bolometric instellation, because X-ray and EUV photons are the proximate drivers of upper-atmosphere escape [1702.03386].

The physical meaning of the shoreline is not unique. In the simple energy-limited picture, atmospheric escape scales with the incident high-energy flux relative to gravitational binding, yielding a characteristic scaling
\[
I_{\rm xuv}\propto v_{\rm esc}^3\sqrt{\rho},
\]
which differs from the empirical \(v_{\rm esc}^4\) law. Zahnle and Catling also emphasized diffusion-limited escape and impact erosion as additional channels: diffusion-limited escape bounds H\(_2\) loss to \(\sim0.5\%\) of planetary mass over 5 Gyr for \(T\sim1000\) K, while impact erosion empirically corresponds to \(v_{\rm imp}\simeq4\mbox{–}5\,v_{\rm esc}\) [1702.03386].

Recent analytic work on secondary atmospheres has added another layer. For N\(_2\)- and CO\(_2\)-dominated atmospheres, hydrodynamic loss is described as an onset problem: a critical XUV flux must first be exceeded before transonic outflow begins. In this framework, escape rates remain approximately linear in \(F_{\rm XUV}\) in a weakly ionized, energy-limited regime, but flatten to \(\dot M\propto F_{\rm XUV}^{1/2}\) once forbidden-line cooling produces a collisional-radiative thermostat. This places the M-dwarf shoreline in a regime where high molecular weight, ionization state, and line cooling are intrinsic to the boundary rather than secondary corrections [2412.05188].

## 2. Why the shoreline shifts around mid-to-late M dwarfs

Earlier M-dwarf shoreline estimates often extrapolated from more massive stars using a canonical XUV scaling,
\[
L_{\rm XUV}\propto L_{\rm bol}^{0.4},
\]
together with the historical XUV instellation relation
\[
\frac{I_{\rm XUV}}{I_{\rm XUV,\oplus}}
=
\frac{a_\oplus^2}{a^2}
\left(\frac{L_{\rm bol}}{L_\odot}\right)^{0.4}.
\]
Pass et al. argued that this prescription is inadequate for mid-to-late M dwarfs because fully convective stars exhibit rotation and activity histories that differ from Sun-like stars and early M dwarfs [2504.01182].

For \(0.1\le M_\star/M_\odot\le0.3\), the revised treatment adopts a saturated fractional X-ray luminosity
\[
\left(\frac{L_{\rm X}}{L_{\rm bol}}\right)_{\rm sat}=10^{-3.05},
\]
and a stellar-mass-dependent saturation lifetime
\[
t_{\rm sat}(M_\star)=5.9-15.4\,\frac{M_\star}{M_\odot}\quad{\rm [Gyr]}.
\]
Neglecting pre-main-sequence and flaring corrections, the historical XUV fluence is written as
\[
I_{\rm XUV}
=
\frac{10^{-3.05}\,L_{\rm bol}\,t_{\rm sat}}{4\pi a^2}.
\]
Including the pre-main-sequence luminosity correction \(f_L\approx1.13\mbox{–}1.22\) and the energetic-flare factor \(f_{\rm flare}=1.25\) gives
\[
I_{\rm XUV}
=
\frac{\bigl(f_{\rm flare}\,10^{-3.05}\bigr)\,\bigl(f_L\,L_{\rm bol}\bigr)\,t_{\rm sat}}{4\pi a^2}.
\]
These prescriptions formalize the statement that the M-dwarf shoreline “recedes”: the same orbital location corresponds to a much larger time-integrated high-energy dose than canonical scaling would imply [2504.01182].

## 3. Revised fluence histories and the receding shoreline

Using the revised activity-lifetime synthesis, the base mid-to-late M-dwarf model yields historical \(I_{\rm XUV}\) values that are on average \(2.1\times\) the Zahnle et al. prescription. When pre-main-sequence evolution and energetic flares are included, the enhancement rises to \(\sim3.1\times\) on average. The largest discrepancy, \(2.3\mbox{–}3.3\times\), occurs near \(M_\star\approx0.2\,M_\odot\), where a canonical \(t_{\rm sat}\approx0.6\) Gyr is a poor match to the revised \(\sim2.8\) Gyr saturation lifetime [2504.01182].

Within the cosmic shoreline paradigm, this shift is consequential for known rocky planets. Pass et al. found that almost all small planets with \(v_{\rm esc}\lesssim20\ {\rm km\,s^{-1}}\) around inactive mid-to-late M dwarfs fall above the updated shoreline, implying likely atmospheric loss. The best atmosphere-retention candidates are the largest terrestrial worlds, with \(v_{\rm esc}\gtrsim20\ {\rm km\,s^{-1}}\). In the same interpretation, small Earth-mass planets around mid-to-late M dwarfs are generally predicted to be airless, while super-Earths with \(R_p\sim1.3\mbox{–}1.8\,R_\oplus\) and \(v_{\rm esc}\approx20\mbox{–}25\ {\rm km\,s^{-1}}\) are the principal candidate class for retained atmospheres [2504.01182].

This conclusion is consistent with analytic secondary-atmosphere escape models. In those models, planets with \(v_{\rm esc}\lesssim v_\oplus\) are expected to exceed the critical XUV threshold during the long saturated phases of active M dwarfs, whereas more massive super-Earths can remain below that threshold. The examples explicitly discussed are TRAPPIST-1 b, placed well into the rapid-loss regime; LHS 1140 c, placed near the boundary; and LHS 1140 b, placed safely in the retention regime [2412.05188].

## 4. Quantitative descriptors and probabilistic generalizations

To rank individual planets, Pass et al. introduced the Atmosphere Retention Metric (ARM), defined as the log-space distance from the shoreline:
\[
\mathrm{ARM}
\equiv
4\,\log_{10}(v_{\rm esc})
-
\log_{10}(I_{\rm XUV})
-
3.16,
\]
with \(v_{\rm esc}\) in \({\rm km\,s^{-1}}\) and \(I_{\rm XUV}\) in Earth units. By construction, \(\mathrm{ARM}<0\) places a planet above the shoreline and implies likely complete atmospheric loss, while \(\mathrm{ARM}>0\) places it below the shoreline and makes retention more probable [2504.01182].

A separate strand of the literature replaces reconstructed XUV histories with a three-dimensional boundary in escape velocity, bolometric flux, and host luminosity:
\[
S_{\rm crit}(v_{\rm esc},L_\star)
=
S_0
\left(\frac{v_{\rm esc}}{v_0}\right)^p
\left(\frac{L_\star}{L_0}\right)^q.
\]
Under the restrictive assumption that one planar boundary applies across a wide parameter space, the inferred indices are
\[
p=6.08^{+0.69}_{-0.48},
\qquad
q=1.25^{+0.31}_{-0.22}.
\]
The condition \(S_{\rm crit}(v=1\,v_\oplus,L_\star)=1\,f_\oplus\) intersects at
\[
\log_{10}(L_\star/L_\odot)\approx-2.22\pm0.21,
\]
roughly spectral type M4.5V, below which Earth-sized planets at nominal habitable-zone instellation are inferred to be more likely to have lost atmospheres [2507.02136].

This three-dimensional framework is explicitly probabilistic rather than deterministic. The atmosphere probability is written as
\[
P(\mathrm{atmosphere})
=
\left[
1+\exp\left(\frac{\log_{10}f-\log_{10}f_{\rm crit}}{w}\right)
\right]^{-1},
\]
where \(w\) is the intrinsic width or “fuzziness” of the boundary. The fit was performed with NUTS in numpyro, using uninformative priors on the slopes and latent true parameters for \(v_{\rm esc}\), \(f\), and \(L_\star\). In this formulation, luminosity serves as an ensemble-level proxy for integrated XUV history, partially avoiding star-by-star historical reconstruction [2507.02136].

## 5. Known planets, boundary cases, and survey design

Applying the revised mid-to-late M-dwarf fluence history to nearby terrestrial planets leaves only a small set of strong retention candidates. Pass et al. reported that the planets with \(\mathrm{ARM}>0\) are the largest known terrestrial worlds around mid-to-late M dwarfs, and that all other known \(R_p<1.8\,R_\oplus\) mid-to-late M-dwarf planets within 50 pc have \(\mathrm{ARM}<0\) [2504.01182].

| Planet | \(v_{\rm esc}\) (km s\(^{-1}\)) | ARM |
|---|---:|---:|
| TOI-715 b | 20.2 | +0.28 |
| TOI-6002 b | 22.2 | +0.12 |
| TOI-1452 b | 22.7 | +0.21 |
| LHS 1140 b | 23.9 | +0.72 |

For these same planets, the corrected fluences tabulated by Pass et al. are \(61\), \(130\), \(110\), and \(43\) in Earth XUV units, respectively [2504.01182].

A temperate example near the boundary is GJ 3378 b. The revised orbital solution gives \(P=21.45\pm0.01\) d, \(m\sin i=2.3\pm0.4\,M_\oplus\), \(a=0.09673\pm0.0008\) AU, \(S=0.91\pm0.09\,S_\oplus\), an estimated \(R_p\simeq1.29\,R_\oplus\), and \(v_{\rm esc}\simeq15\ {\rm km\,s^{-1}}\). Its integrated XUV exposure is \(\approx48\times\) Earth’s, and its best-estimate metric, \(\mathrm{ARM}\approx-0.13\), places it essentially on the dividing line, making it a literal case of a planet “straddling the cosmic shoreline” [2605.16499].

At survey scale, the shoreline has also been framed as a population-level hypothesis test for the JWST Rocky Worlds DDT program. Ih et al. compared Pessimist, Random, Bolometric Shoreline, and XUV Shoreline models using secondary-eclipse depths, endmember bare-rock and shallow-CO\(_2\)-atmosphere templates, and Bayesian Information Criterion model comparison. In that framework, a \(\sim500\) h MIRI/F1500W campaign can either place a 95% upper limit of \(\sim12\%\) on the atmosphere occurrence rate if all planets are bare, or recover strong evidence \((\Delta{\rm BIC}>5)\) for a shoreline trend using a “wide and shallow” strategy that observes \(\sim50\) targets, with 1–2 eclipses on \(\sim40\) dry candidates and 4–18 eclipses on \(\sim8\) wet candidates. The same study identifies F1500W as the optimal single-band choice for distinguishing thick CO\(_2\)-bearing atmospheres from bare rocks [2508.08253].

## 6. Caveats, degeneracies, and competing formulations

A central caveat is that an apparent shoreline can emerge from observational bias rather than escape. Lustig-Yaeger et al. showed that sulfuric-acid clouds in Venus-like atmospheres can generate a “mirage of the cosmic shoreline”: if cloud tops are mistaken for solid surfaces, transmission spectra can imply a monotonic thinning of atmospheres with increasing flux even when the underlying planets retain substantial atmospheres. In their TRAPPIST-1 calculations, transmission spectroscopy probes only \(P_{\rm trans}\sim10^2\mbox{–}10^3\) Pa, whereas JWST/MIRI LRS secondary-eclipse emission in the \(6\,\mu{\rm m}\) window reaches \(P_{\rm emiss}\sim10^5\mbox{–}10^6\) Pa, probing \(3\mbox{–}4\) orders of magnitude deeper. Proposed discriminants include secondary-eclipse emission, pressure-broadened CO\(_2\) line wings, the sign of the \(P_{\rm trans}(\max)\)–semi-major-axis slope, and direct aerosol signatures [1911.09132].

The numerical form of the M-dwarf shoreline is also unsettled. The classical Solar-System relation has slope \(4\) [1702.03386]. The three-dimensional luminosity-aware shoreline gives \(p=6.08^{+0.69}_{-0.48}\) [2507.02136]. The Empirical Cosmic Shoreline anchored to Mars and 55 Cnc e gives
\[
\log_{10}(I_{\rm XUV})\approx5.89\,\log_{10}(v_{\rm esc})-4.49,
\]
and, when revised XUV histories are applied for hosts with \(M<0.35\,M_\odot\), leaves only seven known small planets \((R\lesssim1.7\,R_\oplus)\) securely in the retention zone [2508.12865]. These alternatives differ not only in slope and intercept, but also in whether the independent variable is cumulative XUV fluence, bolometric flux, or a luminosity-augmented proxy.

An additional 2026 study proposed an M-dwarf shoreline
\[
\log_{10} I_{\rm XUV}
=
(6.04\pm0.30)\,\log_{10} v_{\rm esc}
-
(5.35\pm0.15),
\]
together with a mixed-population line
\[
\log_{10} I_{\rm XUV}
=
(4.02\pm0.25)\,\log_{10} v_{\rm esc}
-
(3.21\pm0.20),
\]
derived with a support-vector-machine classifier. That construction explicitly included eight “extragalactic” planets sourced from a “Galactic Senate” survey via Wookiepedia, with Kamino providing the sole clear M-dwarf atmosphere detection in the sample [2603.29743]. This suggests that M-dwarf shoreline inferences remain highly sensitive to sample definition, atmosphere labels, and the treatment of stellar activity histories.

Across these formulations, the recurring empirical result is narrower than the diversity of models might suggest: prolonged M-dwarf activity shifts the atmospheric-retention boundary against small rocky planets, and the most favorable atmosphere-bearing candidates are the highest-\(v_{\rm esc}\) terrestrial worlds, especially those that remain on the retention side after revised fluence histories are applied.

Source: https://www.emergentmind.com/topics/m-dwarf-cosmic-shoreline