Context. The origin of water in the stratospheres of Giant Planets has been an outstanding question ever since its first detection by ISO some 20 years ago. Water can originate from interplanetary dust particles, icy rings and satellites and large comet impacts. Analysis of Herschel Space Observatory observations have proven that the bulk of Jupiter's stratospheric water was delivered by the Shoemaker-Levy 9 impacts in 1994. In 2006, the Cassini mission detected water plumes at the South Pole of Enceladus, placing the moon as a serious candidate for Saturn's stratospheric water. Further evidence was found in 2011, when Herschel demonstrated the presence of a water torus at the orbital distance of Enceladus, fed by the moon's plumes. Finally, water falling from the rings onto Saturn's uppermost atmospheric layers at low latitudes was detected during the final orbits of Cassini's end-of-mission plunge into the atmosphere. Aims. In this paper, we use Herschel mapping observations of water in Saturn's stratosphere to identify its source. Methods. Several empirical models are tested against the Herschel-HIFI and -PACS observations, which were collected on December 30, 2010, and January 2nd, 2011 (respectively). Results. We demonstrate that Saturn's stratospheric water is not uniformly mixed as a function of latitude, but peaking at the equator and decreasing poleward with a Gaussian distribution. We obtain our best fit with an equatorial mole fraction 1.1 ppb and a half-width at half-maximum of 25{\deg}, when accounting for a temperature increase in the two warm stratospheric vortices produced by Saturn's Great Storm of 2010-2011. Conclusions. This work demonstrates that Enceladus is the main source of Saturn's stratospheric water.

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.

Neutrals sourced directly from Enceladus's plumes are initially confined to a dense neutral torus in Enceladus's orbit around Saturn. This neutral torus is redistributed by charge exchange, impact/photodissociation, and neutral-neutral collisions to produce Saturn's neutral clouds. Here we consider the former processes in greater detail than in previous studies. In the case of dissociation, models have assumed that OH is produced with a single speed of 1 km/s, whereas laboratory measurements suggest a range of speeds between 1 and 1.6 km/s. We show that the high-speed case increases dissociation's range of influence from 9 to 15 Rs. For charge exchange, we present a new modeling approach, where the ions are followed within a neutral background, whereas neutral cloud models are conventionally constructed from the neutrals' point of view. This approach allows us to comment on the significance of the ions' gyrophase at the moment charge exchange occurs. Accounting for gyrophase: (1) has no consequence on the H2O cloud; (2) doubles the local density of OH at the orbit of Enceladus; and (3) decreases the oxygen densities at Enceladus's orbit by less than 10%. Finally, we consider velocity-dependent, as well as species-dependent cross sections and find that the oxygen cloud produced from charge exchange is spread out more than H2O, whereas the OH cloud is the most confined.

In this paper, we look at space weathering processes on the icy surface of Jupiter's moon Europa. The heavy energetic ions of the Jovian plasma (H+, O+, S+, C+) can erode the surface of Europa via ion sputtering (IS), ejecting up to 1000 H2O molecules per ion. UV Photons impinging the Europa's surface can also result in neutral atom release via photon-stimulated desorption (PSD) and chemical change (photolysis). In this work, we study the efficiency of the IS and PSD processes for ejecting water molecules, simulating the resulting neutral H2O density. We also estimate the contribution to the total neutral atom release by the Ion Backscattering (IBS) process. Moreover, we estimate the possibility of detecting the sputtered high energy atoms, in order to distinguish the action of the IS process from other surface release mechanisms. Our main results are: 1) The most significant sputtered-particle flux and the largest contribution to the neutral H2O-density come from the incident S+ ions; 2) The H2O density produced via PSD is lower than that due to sputtering by ~1.5 orders of magnitude; 3) In the energy range below 1 keV, the IBS can be considered negligible for the production of neutrals, whereas in the higher energy range it becomes the dominant neutral emission mechanism; 4) the total sputtering rate for Europa is 2.0\cdot 1027 H2O s-1; 5) the fraction of escaping H2O via IS is 22% of the total sputtered population, while the escape fraction for H2O produced by PSD is 30% of the total PSD population. Since the PSD exosphere is lower than the IS one, the major agent for Europa's surface erosion is IS on both the non-illuminated and illuminated side. Lastly, the exospheric neutral density, estimated from the Galileo electron density measurements appears to be higher than that calculated for H2O alone; this favours the scenario of the presence of O2 produced by radiolysis and photolysis.

We study the orbits, tidal heating and mass loss from satellites around close-in gas giant exoplanets. The focus is on large satellites which are potentially observable by their transit signature. We argue that even Earth-size satellites around hot Jupiters may be immune to destruction by orbital decay; detection of such a massive satellite would strongly constrain theories of tidal dissipation in gas giants, in a manner complementary to orbital circularization. The star's gravity induces significant periodic eccentricity in the satellite's orbit. The resulting tidal heating rates, per unit mass, are far in excess of Io's and dominate radioactive heating out to planet orbital periods of months for reasonable satellite tidal $Q$. Inside planet orbital periods of about a week, tidal heating can completely melt the satellite. Lastly, we compute an upper limit to the satellite mass loss rate due to thermal evaporation from the surface, valid if the satellite's atmosphere is thin and vapor pressure is negligible. Using this upper limit, we find that although rocky satellites around hot Jupiters with orbital periods less than a few days can be significantly evaporated in their lifetimes, detectable satellites suffer negligible mass loss at longer orbital periods.