The upper atmospheres of Mars and Titan, as well as those on many other planetary bodies, exhibit significant density variations vs. altitude that are interpreted as gravity waves. Such data is then used to extract vertical temperature profiles, even when such perturbations propagate through the transition region from a collision dominated regime and into a planet's exosphere. Since the temperature profile is critical for describing the upper atmospheric heating and evolution, we use molecular kinetic simulations to describe transient perturbations in a Mars-like upper atmosphere. We show that the standard methods for extracting the temperature profile can fail dramatically so that molecular kinetic simulations, calibrated to observed density profiles, are needed in this region of a planet's atmosphere.

Escape of an early atmosphere from Titan, during which time NH3 could be converted by photolysis into the present N2 dominated atmosphere, is an important problem in planetary science. Recently Gilliam and Lerman (2014) estimated escape driven by the surface temperature and pressure, which we show gave loss rates that are orders of magnitude too large. Their model, related to Jeans escape from an isothermal atmosphere, was used to show that escape driven only by surface heating would deplete the atmospheric inventory of N for a suggested Titan accretion temperature of ~355 K. Therefore, they concluded that the accretion temperature must be lower in order to retain the present Titan atmosphere. Here we show that the near surface atmospheric temperature is essentially irrelevant for determining the atmospheric loss rate from Titan and that escape is predominantly driven by solar heating of the upper atmosphere. We also give a rough estimate of the escape rate in the early solar system (~10^4 kg/s) consistent with an inventory of nitrogen being available over the time period suggested by Atreya et al. (1978) for conversion of NH3 into N2.

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.