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.

Ganymede's atmosphere is produced by radiative interactions with its surface, sourced by the Sun and the Jovian plasma. The sputtered and thermally desorbed molecules are tracked in our Exospheric Global Model (EGM), a 3-D parallelized collisional model. This program was developed to reconstruct the formation of the upper atmosphere/exosphere of planetary bodies interacting with solar photon flux and magnetospheric and/or the solar wind plasmas. Here, we describe the spatial distribution of the H$_2$O and O$_2$ components of Ganymede's atmosphere, and their temporal variability along Ganymede's rotation around Jupiter. In particular, we show that Ganymede's O$_2$ atmosphere is characterized by time scales of the order of Ganymede's rotational period with Jupiter's gravity being a significant driver of the spatial distribution of the heaviest exospheric components. Both the sourcing and the Jovian gravity are needed to explain some of the characteristics of the observed aurora emissions. As an example, the O$_2$ exosphere should peak at the equator with systematic maximum at the dusk equator terminator. The sputtering rate of the H$_2$O exosphere should be maximum on the leading hemisphere because of the shape of the open/close field lines boundary and displays some significant variability with longitude.

The evolution of Europa's water-product exosphere over its 85-hour day, based on current models, has not been shown to exhibit any diurnal asymmetries. Here we simulate Europa's exosphere using a 3-D Monte Carlo routine including, for the first time, the role of Europa's rotation on the evolution of exospheric molecules throughout the orbit. We focus on O$_2$, sputtered by a trailing hemisphere source with a temperature-dependence under isotropic plasma conditions. We find that the O$_2$ component, while global, is not homogenous in Europa local time. Rather, the O$_2$ accumulates at the dusk hemisphere. When rotation is explicitly excluded, no diurnal asymmetries exist. We find that the assumed thermal-dependence on the O$_2$ source is critical for a diurnal asymmetry: the diurnal surface temperature profile is imprinted on to the near-surface O$_2$ atmosphere, due to the small hop times of the non-adsorbing O$_2$ effectively rotating with Europa. Altogether, our simulations conclude that the dusk-over-dawn asymmetry is driven by Europa's day-night O$_2$ cycle synchronized with Europa's orbital period based on our model assumptions on O$_2$ production and loss. This conclusion is in agreement with the recent understanding that a non-adsorbing, rotating O$_2$ source peaking at noon will naturally accumulate from dawn-to-dusk, should the O$_2$ lifetime be sufficiently long compared to the orbital period. Lastly we compare hemispherically-averaged dusk-over-dawn ratios to the recently observed oxygen emission data by the Hubble Space Telescope. We find that while the simulations are globally consistent with brighter oxygen emission at dusk than at dawn, the orbital evolution of the asymmetries in our simulations can be improved by ameliorating the O$_2$ source & loss rates, and possibly adsorption onto the regolith.