A Theory of transmission spectroscopy of planetary winds: Spectral-line saturation and limits on mass-loss inference

TL;DR

This work develops an analytic theory for transmission spectroscopy of hydrodynamically escaping exoplanet atmospheres by coupling standard transit geometry to a steady, isothermal Parker wind. It yields a closed-form expression for the effective transit radius using a Lambert-$W$ inversion, revealing two regimes: an opacity-limited regime where $R_{\rm eff}$ constrains the mass-loss rate, and a saturation-limited regime where line cores saturate and the inversion breaks down, leaving the observable signal governed by geometric extent. The key contributions include the explicit saturation boundary $\sigma(\lambda)\dot{M} \le C_{\rm sat}$, the branch choice $W_{-1}$ for physical solutions, and a physically transparent link between transmission spectra and escape theory via the sonic radius $r_s$ and Jeans parameter. These results explain why strong line cores often fail to constrain $\dot{M}$ and why weaker lines and wings remain diagnostic, providing practical guidance for observations and retrievals and a bridge between analytic insight and numerical radiative transfer in planetary winds.

Abstract

Transmission spectroscopy is a key technique in the characterization of exoplanet atmospheres and has been widely applied to planets undergoing hydrodynamic escape. While a robust analytic theory exists for transmission spectra of hydrostatic atmospheres, the corresponding interpretation for escaping atmospheres has so far relied on numerical modeling. In this work, we develop a theory of transmission spectroscopy in hydrodynamically escaping atmospheres by coupling the standard transmission geometry to a steady-state, spherically symmetric, isothermal outflow. This approach yields closed-form expressions and allows the optical depth inversion problem to be examined. The analytic solution reveals that transmission spectroscopy of planetary winds naturally separates into two regimes. In an opacity-limited regime, transmission depths retain sensitivity to the atmospheric mass-loss rate. Beyond a critical threshold, however, spectral-line cores become saturated and no longer provide a unique constraint on the mass flux. This transition is marked by a sharp analytic boundary of the form $σ(λ)\times \dot M \le C_{sat}$, where $C_{sat}$ is a constant set by the thermodynamic and geometric properties of the wind. This condition specifies when the inversion between transmission depth and mass-loss rate admits a real solution. Once it is violated, the effective transit radius is no longer controlled by opacity or mass loss, but by the geometric extent of the absorbing wind. These results demonstrate that spectral-line saturation in transmission spectroscopy corresponds to a fundamental loss of invertibility between absorption and atmospheric mass loss, rather than a gradual weakening of sensitivity. The theory provides a physically transparent explanation for why strong transmission line cores often fail to constrain mass-loss rates, while weaker lines and line wings remain diagnostic.

Figures (6)

• Figure 1: Geometry of transmission spectroscopy in a planetary wind. A planet of radius $R_p$ is surrounded by a bound hydrostatic atmosphere and an outer hydrodynamically escaping region. Stellar rays intersect the atmosphere along chords of impact parameter $b$, accumulating a wavelength-dependent slant optical depth $\tau(b)$ along the line of sight toward the observer. The effective transit radius $R_{\mathrm{eff}}(\lambda)$ is defined by the impact parameter for which $\tau(b)$ reaches the reference value $\tau_\ast$. The dashed blue circle marks the sonic radius $r_s$ and the approximate location of the XUV photosphere, $R_{\mathrm{XUV}}$, is also indicated.
• Figure 2: Radial velocity and density profiles of the Parker wind for the hot Jupiter HD 209458 b. With black solid lines we show the profiles as calculated by the Parker wind formulation while in dashed green we show the subsonic approximation (see Section \ref{['sec:derivation']}). For comparison, we overplot the radii where the optical depth of the NUV, FUV, EUV are reaching unity, $\tau\approx 1$.
• Figure 3: Optical depth as a function of planet radius for EUV, FUV and NUV wavelength bands as well as approximations for the same bands overplotted. With shaded stripes we show the range $\tau=0.56$--1 expressed in terms of the corresponding planet radii for the ultraviolet bands shown.
• Figure 4: Up: Transmission spectrum of HD209458b under the hydrostatic/stable thermosphere hypothetical scenario. With shaded regions we indicate the dominant absorption species at each wavelength. Bottom: Comparison of HD209458b hydrodynamic atmosphere's transmission spectrum from numerical integration including line broadening effects from bulk flow and thermal motion, compared to the analytic model derived in Eq. \ref{['eq:Reff_W_minus1']}.
• Figure 5: We define the analytic coverage as the fraction of wavelength points for which the Lambert-$W$ argument satisfies $-1/e \le z(\lambda) < 0$, ensuring a real-valued effective transit radius.