The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox.

Observations of water ice clouds and aerosols on Mars can provide important insights into the complexity of the water cycle. Recent observations have indicated an important link between dust activity and the water cycle, as intense dust activity can significantly raise the hygropause, and subsequently increase the escape of water after dissociation in the upper atmosphere. Here present observations from NOMAD/TGO that investigate the variation of water ice clouds in the perihelion season of Mars Year 34 (April 2018-19), their diurnal and seasonal behavior, and the vertical structure and microphysical properties of water ice and dust. These observations reveal the recurrent presence of a layer of mesospheric water ice clouds subsequent to the 2018 Global Dust Storm. We show that this layer rose from 45 to 80 km in altitude on a timescale of days from heating in the lower atmosphere due to the storm. In addition, we demonstrate that there is a strong dawn dusk asymmetry in water ice abundance, related to nighttime nucleation and subsequent daytime sublimation. Water ice particle sizes are retrieved consistently and exhibit sharp vertical gradients (from 0.1 to 4.0 um), as well as mesospheric differences between the Global Dust Storm (<0.5 um) and the 2019 regional dust storm (1.0 um), which suggests differing water ice nucleation efficiencies. These results form the basis to advance our understanding of mesospheric water ice clouds on Mars, and further constrain the interactions between water ice and dust in the middle atmosphere.

We adapted a multi-species escape model, developed for close-in extrasolar planets, to calculate the escape rates of CH4 and N2 from Pluto. In the absence of escape, CH4 should overtake N2 as the dominant species below the exobase. The CH4 profile depends strongly on the escape rate, however, and the typical escape rates predicted for Pluto lead to a nearly constant mixing ratio of less than 1 % below the exobase. In this case the CH4 escape rate is only 5-10 % of the N2 escape rate. Observations of the CH4 profile by the New Horizons/ALICE spectrograph can constrain the CH4 escape rate and provide a unique test for escape models.

The equations of gas dynamics are extensively used to describe atmospheric loss from solar system bodies and exoplanets even though the boundary conditions at infinity are not uniquely defined. Using molecular-kinetic simulations that correctly treat the transition from the continuum to the rarefied region, we confirm that the energy-limited escape approximation is valid when adiabatic expansion is the dominant cooling process. However, this does not imply that the outflow goes sonic. In fact in the sonic regime, the energy limited approximation can significantly under estimate the escape rate. Rather large escape rates and concomitant adiabatic cooling can produce atmospheres with subsonic flow that are highly extended. Since this affects the heating rate of the upper atmosphere and the interaction with external fields and plasmas, we give a criterion for estimating when the outflow goes transonic in the continuum region. This is applied to early terrestrial atmospheres, exoplanet atmospheres, and the atmosphere of the ex-planet, Pluto, all of which have large escape rates. The paper and its erratum, combined here, are published: ApJL 768, L4 (2013); ApJ, 779, L30 (2013).

Predicting the rate of escape and thermal structure of Pluto's upper atmosphere in preparation for the New Horizons Spacecraft encounter in 2015 is important for planning and interpreting the expected measurements. Having a moderate Jeans parameter Pluto's atmosphere does not fit the classic definition of Jeans escape for light species escaping from the terrestrial planets, nor does it fit the hydrodynamic outflow from comets and certain exoplanets. It has been proposed for some time that Pluto lies in the region of slow hydrodynamic escape. Using a hybrid fluid/molecular-kinetic model, we previously demonstrated the typical implementation of this model fails to correctly describe the appropriate temperature structure for the upper atmosphere for solar minimum conditions. Here we use a time-dependent solver to allow us to extend those simulations to higher heating rates and we examine fluid models in which Jeans-like escape expressions are used for the upper boundary conditions. We compare these to hybrid simulations of the atmosphere under heating conditions roughly representative of solar minimum and mean conditions as these bracket conditions expected during the New Horizon encounter. Although we find escape rates comparable to those previously estimated by the slow hydrodynamic escape model, and roughly consistent with energy limited escape, our model produces a much more extended atmosphere with higher temperatures roughly consistent with recent observations of CO. Such an extended atmosphere will be affected by Charon and will affect Pluto's interaction with the solar wind at the New Horizon encounter. Since we have previously shown that such models can be scaled, these results have implications for modeling exoplanet atmospheres for which the energy limited escape approximation is often used.

A combined fluid/kinetic model is developed to calculate thermally driven escape of N2 from Pluto's atmosphere for two solar heating conditions: no heating above 1450 km and solar minimum heating conditions. In the combined model, one-dimensional fluid equations are applied for the dense part of the atmosphere, while the exobase region is described by a kinetic model and calculated by the direct simulation Monte Carlo method. Fluid and kinetic parts of the model are iteratively solved in order to maintain constant total mass and energy fluxes through the simulation region. Although the atmosphere was found to be highly extended, with an exobase altitude at ~6000 km at solar minimum, the outflow remained subsonic and the escape rate was within a factor of two of the Jeans rate for the exobase temperatures determined. This picture is drastically different from recent predictions obtained solely using a fluid model which, in itself, requires assumptions about atmospheric density, flow velocity and energy flux carried away by escaping molecules at infinity. Gas temperature, density, velocity and heat flux versus radial distance are consistent between the hydrodynamic and kinetic model up to the exobase, only when the energy flux across the lower boundary and escape rate used to solve the hydrodynamic equations is obtained from the kinetic model. This limits the applicability of fluid models to atmospheric escape problems. Finally, the recent discovery of CO at high altitudes, the effect of Charon and the conditions at the New Horizon encounter are briefly considered.

Thermally-driven atmospheric escape evolves from an organized outflow (hydrodynamic escape) to escape on a molecule by molecules basis (Jeans escape) with increasing Jeans parameter, the ratio of the gravitational to thermal energy of molecules in a planet's atmosphere. This transition is described here using the direct simulation Monte Carlo method for a single component spherically symmetric atmosphere. When the heating is predominantly below the lower boundary of the simulation region, R0, and well below the exobase, this transition is shown to occur over a surprisingly narrow range of Jeans parameters evaluated at R0: {\lambda}0 ~ 2-3. The Jeans parameter {\lambda}0 ~ 2.1 roughly corresponds to the upper limit for isentropic, supersonic outflow and for {\lambda}0 >3 escape occurs on a molecule by molecule basis. For {\lambda}0 > ~6, it is shown that the escape rate does not deviate significantly from the familiar Jeans rate evaluated at the nominal exobase, contrary to what has been suggested. Scaling by the Jeans parameter and the Knudsen number, escape calculations for Pluto and an early Earth's atmosphere are evaluated, and the results presented here can be applied to thermally-induced escape from a number of solar and extrasolar planetary bodies.