The ultraviolet reflectance of Enceladus: Implications for surface composition

The reflectance of Saturn’s moon Enceladus has been measured at far ultraviolet (FUV) wavelengths (115-190 nm) by Cassini’s Ultraviolet Imaging Spectrograph (UVIS). At visible and near infrared (VNIR) wavelengths Enceladus’ reflectance spectrum is very bright, consistent with a surface composed primarily of H₂O ice. At FUV wavelengths, however, Enceladus is surprisingly dark - darker than would be expected for pure water ice. Previous analyses have focused on the VNIR spectrum, comparing it to pure water ice (Cruikshank, D.P., Owen, T.C., Dalle Ore, C., Geballe, T.R., Roush, T.L., de Bergh, C., Sandford, S.A., Poulet, F., Benedix, G.K., Emery, J.P. [2005] Icarus, 175, 268-283) or pure water ice plus a small amount of NH₃ (Emery, J.P., Burr, D.M., Cruikshank, D.P., Brown, R.H., Dalton, J.B. [2005] Astron. Astrophys., 435, 353-362) or NH₃ hydrate (Verbiscer, A.J., Peterson, D.E., Skrutskie, M.F., Cushing, M., Helfenstein, P., Nelson, M.J., Smith, J.D., Wilson, J.C. [2006] Icarus, 182, 211-223). We compare Enceladus’ FUV spectrum to existing laboratory measurements of the reflectance spectra of candidate species, and to spectral models. We find that the low FUV reflectance of Enceladus can be explained by the presence of a small amount of NH₃ and a small amount of a tholin in addition to H₂O ice on the surface. The presence of these three species (H₂O, NH₃, and a tholin) appears to satisfy not only the low FUV reflectance and spectral shape, but also the middle-ultraviolet to visible wavelength brightness and spectral shape. We expect that ammonia in the Enceladus plume is transported across the surface to provide a global coating.     

An assessment and test of Enceladus as an important source of Saturn’s ring atmosphere and ionosphere

The existence of an oxygen exosphere and ionosphere in Saturn’s main ring region has been confirmed by the Saturn Orbital Insertion (SOI) observations of the Cassini spacecraft. Through the ion-molecule collisions, the ring atmosphere could serve as a source of ions throughout Saturn’s magnetosphere. If photolysis of ice in the main rings is the dominant source of O₂, then the complex structure of the ring atmosphere/ionosphere and the injection rate of neutral O₂ will be subject to modulation by the seasonal variation of Saturn along its orbit (Tseng, Wei-Ling, Ip, W.-H., Johnson, R.E., Cassidy, T.A., Erlod, M.K. [2010]. Icarus 206, 382-389). In addition, the radio and plasma wave science (RPWS) instrument onboard Cassini found that a large amount of the Enceladus-originated water-group plasma would be deposited on the outer edge of the A ring (Farrell, W.M., Kaiser, M.L., Gurnett, D.A., Kurth, W.S., Persoon, A.M., Wahlund, J.E., Canu, P. [2008]. Geophys. Res. Lett. 35, L02203). A large amount of Enceladus’ plume neutrals (water-group neutrals) would collide with the main rings through collisional interaction with the ambient neutrals and plasma ions (Jurac, S., Richardson, J.D. [2007]. Geophys. Res. Lett. 34, L08102; Cassidy, T.A., Johnson, R.E. [2010]. Icarus, in press). These absorbed ions and neutrals could be recycled to neutral oxygen molecules via grain-surface chemistry to contribute the ring oxygen atmosphere. In this work, we have examined the mass budget of the ring oxygen atmosphere of Saturn taking into account such an “exogenic” source. The maximum O₂ source rate from recycling of Enceladus-originated plasma and neutrals is probably comparable or higher to the one from photolytic decomposition of ices. In the above case, the neutral O₂ source rate would be independent of the solar insolation angle. Therefore, even at Saturn’s Equinox, the extended oxygen atmosphere still could be an important supplier of oxygen ions in the saturnian magnetosphere. We have performed several studies for different recycling source rates from Enceladus. These predictions need further the Cassini Plasma Spectrometer (CAPS) and the Magnetospheric Imaging Instrument (MIMI) observations to be verified in future.     

Obliquity tides do not significantly heat Enceladus

Recently, Tyler [Tyler, R.H., 2009. Geophys. Res. Lett. 36, L15205; Tyler, R., 2011. Icarus, 211, 770-779] proposed that the tide due to an obliquity of greater than 0.1° might drive resonant flow in a liquid ocean at Enceladus, and that dissipation of the ocean’s kinetic energy may be an alternate source for the observed global heat flux. While there is currently no measurement of Enceladus’ obliquity, dissipation is expected to drive the spin pole to a Cassini state. Under this assumption, we find that Enceladus should occupy Cassini state 1 and that the obliquity of Enceladus should be less than 0.0015° for values of the degree-2 gravity coefficient C_2,2 between 1.0 × 10⁻³ and 2.5 × 10⁻³. Unless there is a significant free obliquity or the gravity coefficient C_2,2 has been significantly overestimated, it is unlikely that obliquity-driven flow in a subsurface ocean is the source of the extreme heat on Enceladus.     

Geysers of Enceladus: Quantitative analysis of qualitative models

Aspects of two qualitative models of Enceladus’ dust plume—the so-called “Cold Faithful” [Porco, C.C., et al., 2006. Cassini observes the active south pole of Enceladus. Science 311, 1393-1401; Ingersoll, A.P., et al., 2006. Models of the Enceladus plumes. In: Bulletin of the American Astronomical Society, vol. 38, p. 508] and “Frigid Faithful” [Kieffer, S.W., et al., 2006. A clathrate reservoir hypothesis for Enceladus’ south polar plume. Science 314, 1764; Gioia, G., et al., 2007. Unified model of tectonics and heat transport in a Frigid Enceladus. Proc. Natl. Acad. Sci. 104, 13578-13591] models—are analyzed quantitatively. The former model assumes an explosive boiling of subsurface liquid water, when pressure exerted by the ice crust is suddenly released due to an opening crack. In the latter model the existence of a deep shell of clathrates below Enceladus’ south pole is conjectured; clathrates can decompose explosively when exposed to vacuum through a fracture in the outer icy shell. For the Cold Faithful model we estimate the maximal velocity of ice grains, originating from water splashing in explosive boiling. We find that for water near the triple point this velocity is far too small to explain the observed plume properties. For the Frigid Faithful model we consider the problem of momentum transfer from gas to ice particles. It arises since any change in the direction of the gas flow in the cracks of the shell requires re-acceleration of the entrained grains. While this effect may explain the observed speed difference of gas and grains if the gas evaporates from triple point temperature (273.15 K) [Schmidt, J., et al., 2008. Formation of Enceladus dust plume. Nature 451, 685], the low temperatures of the Frigid Faithful model imply a too dilute vapor to support the observed high particle fluxes in Enceladus’ plume.     

How the Enceladus dust plume feeds Saturn’s E ring

Pre-Cassini models of Saturn’s E ring [Horányi, M., Burns, J., Hamilton, D., 1992. Icarus 97, 248-259; Juhász, A., Horányi, M., 2002. J. Geophys. Res. 107, 1-10] failed to reproduce its peculiar vertical structure inferred from Earth-bound observations [de Pater, I., Martin, S.C., Showalter, M.R., 2004. Icarus 172, 446-454]. After the discovery of an active ice-volcanism of Saturn’s icy moon Enceladus the relevance of the directed injection of particles for the vertical ring structure of the E ring was swiftly recognised [Juhász, A., Horányi, M., Morfill, G.E., 2007. Geophys. Res. Lett. 34, L09104; Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., Economou, T., Schmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. However, simple models for the delivery of particles from the plume to the ring predict a too small vertical ring thickness and overestimate the amount of the injected dust.Here we report on numerical simulations of grains leaving the plume and populating the dust torus of Enceladus. We run a large number of dynamical simulations including gravity and Lorentz force to investigate the earliest phase of the ring particle life span. The evolution of the electrostatic charge carried by the initially uncharged grains is treated selfconsistently. Freshly ejected plume particles are moving in almost circular orbits because the Enceladus orbital speed exceeds the particles’ ejection speeds by far. Only a small fraction of grains that leave the Hill sphere of Enceladus survive the next encounter with the moon. Thus, the flux and size distribution of the surviving grains, replenishing the ring particle reservoir, differs significantly from the flux and size distribution of the particles freshly ejected from the plume. Our numerical simulations reproduce the vertical ring profile measured by the Cassini Cosmic Dust Analyzer (CDA) [Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., EconoDmou, T., Smchmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. From our simulations we calculate the deposition rates of plume particles hitting Enceladus’ surface. We find that at a distance of 100 m from a jet a 10 m sized ice boulder should be covered by plume particles in 10⁵-10⁶ years.     

On the thermal history of Saturn’s satellites Titan and Enceladus

The thermal histories of two geologically active satellites of Saturn—Titan and Enceladus—are discussed. During the Cassini mission, it was found that there are both nitrogen-containing compounds—NH₃ and N₂-and CO₂ and CH₄ in the water plumes of Enceladus; at that, ammonia is the prevailing form. This may testify that during evolution, the material of the satellite was warmed up to T ∼ 500–600 K, when NH₃ (the form of nitrogen capable of being accreted) could only be partly converted into N₂. Contrary to Enceladus, the temperature inside Titan probably reached values higher than 800 K or even higher than 1000 K, since the process of the chemical dissociation of ammonia was completely finished on this satellite and its atmosphere contains only molecular nitrogen. While the internal heating of Titan up to high temperatures can be explained by its large mass, the heating source for Enceladus’ interior is far from evident. Such traditional heating sources as the energy of gravitational differentiation and the radiogenic heating due to shortliving ²⁶Al and ⁶⁰Fe could not be effective. The first one is because of the small size of Enceladus (RE ≈ 250 km), and the inefficiency of the second one is caused by the fact that the satellite was formed not earlier than 8–10 Myr after the formation of calcium and aluminum-enriched inclusions in carbonaceous chondrites (CAI), i.e., after ²⁶Al had completely decayed. In the present paper, we propose other heating mechanisms-the heat of long-living radioactive elements and tidal heat, which could provide the observed chemical composition of the water plumes of Enceladus rather than only the differentiation of its protomatter into the ironstone core and the ice mantle.     

Recent orbital evolution and the internal structures of Enceladus and Dione

We study the orbital behavior of Saturn’s satellites Enceladus and Dione during their passage through the 2:1 mean-motion resonances to constrain their interior structures, parameterized by the quantity k₂/Q (assumed constant). Enceladus’ evolution after escape from the second-order e-Enceladus e-Dione resonance requires that (k₂/Q)_Enceladus<8×10^-4, for that Q_Saturn>18,000. This result is in agreement with [Meyer, J., Wisdom, J., 2008b. Icarus 193, 213-223]. The present-day libration amplitude of Enceladus requires that (k₂/Q)_Enceladus>1.2×10^-4, assuming that Q_Saturn<10⁵. Dione’s present-day eccentricity indicates that (k₂/Q)_Dione?3×10^-4 for Q_Saturn>18,000. Assuming Maxwellian viscoelastic behavior, we find that for Enceladus a convective ice shell overlying an ocean is too dissipative to match the orbital constraints. We conclude that a conductive shell overlying an ocean is more likely, and discuss the implications of this result. Dione’s ice shell is also likely to be conductive, but our results are less constraining.     

Enceladus: Present internal structure and differentiation by early and long-term radiogenic heating

Pre-Cassini images of Saturn's small icy moon Enceladus provided the first indication that this satellite has undergone extensive resurfacing and tectonism. Data returned by the Cassini spacecraft have proven Enceladus to be one of the most geologically dynamic bodies in the Solar System. Given that the diameter of Enceladus is only about 500 km, this is a surprising discovery and has made Enceladus an object of much interest. Determining Enceladus' interior structure is key to understanding its current activity. Here we use the mean density of Enceladus (as determined by the Cassini mission to Saturn), Cassini observations of endogenic activity on Enceladus, and numerical simulations of Enceladus' thermal evolution to infer that this satellite is most likely a differentiated body with a large rock-metal core of radius about 150 to 170 km surrounded by a liquid water-ice shell. With a silicate mass fraction of 50% or more, long-term radiogenic heating alone might melt most of the ice in a homogeneous Enceladus after about 500 Myr assuming an initial accretion temperature of about 200 K, no subsolidus convection of the ice, and either a surface temperature higher than at present or a porous, insulating surface. Short-lived radioactivity, e.g., the decay of ²⁶Al, would melt all of the ice and differentiate Enceladus within a few million years of accretion assuming formation of Enceladus at a propitious time prior to the decay of ²⁶Al. Long-lived radioactivity facilitates tidal heating as a source of energy for differentiation by warming the ice in Enceladus so that tidal deformation can become effective. This could explain the difference between Enceladus and Mimas. Mimas, with only a small rock fraction, has experienced relatively little long-term radiogenic heating; it has remained cold and stiff and less susceptible to tidal heating despite its proximity to Saturn and larger eccentricity than Enceladus. It is shown that the shape of Enceladus is not that of a body in hydrostatic equilibrium at its present orbital location and rotation rate. The present shape could be an equilibrium shape corresponding to a time when Enceladus was closer to Saturn and spinning more rapidly, or more likely, to a time when Enceladus was spinning more rapidly at its present orbital location. A liquid water layer on Enceladus is a possible source for the plume in the south polar region assuming the survivability of such a layer to the present. These results could place Enceladus in a category similar to the large satellites of Jupiter, with the core having a rock-metal composition similar to Io, and with a deep overlying ice shell similar to Europa and Ganymede. Indeed, the moment of inertia factor of a differentiated Enceladus, C/MR², could be as small as that of Ganymede, about 0.31.     

Ring torque on janus and the melting of Enceladus

Voyager 2 images show parts of Enceladus' surface to be very smooth, lacking craters down to the resolution limit of 4 km. This absence of craters indicates geologically recent resurfacing, probably due to internal melting. However, calculations of current heating mechanisms, including radioactive decay and tidal heating due to Enceladus' resonance with Dione, yield heating rates too small to cause melting. The orbital mean motion of Janus (1980S1) is slightly less than twice that of Enceladus and, according to theoretical calculations, is currently decreasing as Janus' orbit evolves outward due to resonant torques from Saturn's rings. If Janus were ever locked into a stable 2:1 orbital commensurability with Enceladus, the resulting angular momentum transfer could have sufficiently enhanced the eccentricity of Enceladus' orbit for the ensuing tidal heating to have melted Enceladus' interior. The existence of a Laplace-like three-body resonance including Dione, although unlikely, would have increased heating. If Janus were indeed held in resonance with Enceladus until recently (10⁷–10⁸ years B.P.) when the lock was disrupted by an unspecified event (possibly a catastrophic collision which simultaneously created the coorbital pair, or by the influence of Dione) both the recent internal activity of Enceladus and the proximity of Janus to Saturn's rings may be explained. However, the predicted rapid time scale for ring evolution due to resonant torques from Saturn's inner moons remains a major problem.     

Saturn's E ring: I. CCD observations of March 1980

The tenuous E ring of Saturn is found to commence abruptly at 3 Saturn radii, to peak sharply in the vicinity of the orbit of the satellite Enceladus (about 4 radii), and to spread out thinly to more than 8 radii. This distribution strongly suggests it to be associated with Enceladus and perhaps to be material ejected from Enceladus. The spread of E-ring material above and below the ring plane is greater in its tenuous outskirts than in its denser inner region, suggesting that the E ring may be at an early stage in its evolution. Thus far, our analysis reveals only a marginal variation of the ring with time or Enceladus azimuth. In this paper we describe the special instrumentation used for photometric observations of the E ring, and we present some of the data obtained in March 1980. In Paper II we shall derive the three-dimensional distribution of material in the E ring and discuss its cosmogonic implications.     

Accretion of Saturn’s mid-sized moons during the viscous spreading of young massive rings: Solving the paradox of silicate-poor rings versus silicate-rich moons

The origin of Saturn’s inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn’s rings is debated. Charnoz et al. [Charnoz, S., Salmon J., Crida A., 2010. Nature 465, 752–754] introduced the idea that the smallest inner moons could form from the spreading of the rings’ edge while Salmon et al. [Salmon, J., Charnoz, S., Crida, A., Brahic, A., 2010. Icarus 209, 771–785] showed that the rings could have been initially massive, and so was the ring’s progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested by Canup [Canup, R., 2010. Nature 468, 943–946]. However, the process driving mid-sized moon accretion from the icy debris disks has not been investigated in details. In particular, Canup’s (2010) model does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea’s mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Q_p < 2000) is required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites’ cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process). These moons may be either mostly icy, or, if they contain a significant amount of rock, already differentiated from the ice without the need for radiogenic heating. The resulting inner mid-sized moons may be significantly younger than the Solar System and a ∼1 Gyr formation delay is possible between Mimas and Rhea. The rings resulting from this process would evolve to a state compatible with current mass estimates of Saturn’s rings, and nearly devoid of silicates, apart from isolated silicate chunks coated with ice, interpreted as today Saturn’s rings’ propellers and ring-moons (like Pan or Daphnis).     

Europa's atmosphere, gas tori, and magnetospheric implications

A two-dimensional kinetic model calculation for the water group species (H₂O, H₂, O₂, OH, O, H) in Europa's atmosphere is undertaken to determine its basic compositional structure, gas escape rates, and velocity distribution information to initialize neutral cloud model calculations for the most important gas tori. The dominant atmospheric species is O₂ at low altitudes and H₂ at higher altitudes with average day-night column densities of 4.5×¹⁰¹⁴ and 7.7×¹⁰¹³ cm⁻², respectively. H₂ forms the most important gas torus with an escape rate of ∼2×¹⁰²⁷ s⁻¹ followed by O with an escape rate of ∼5×¹⁰²⁶ s⁻¹, created primarily as exothermic O products from O₂ dissociation by magnetospheric electrons. The circumplanetary distributions of H₂ and O are highly peaked about the satellite location and asymmetrically distributed near Europa's orbit about Jupiter, have substantial forward clouds extending radially inward to Io's orbit, and have spatially integrated cloud populations of 4.2×¹⁰³³ molecules for H₂ and 4.0×¹⁰³² atoms for O that are larger than their corresponding populations in Europa's local atmosphere by a factor of ∼200 and ∼1000, respectively. The cloud population for H₂ is a factor of ∼3 times larger than that for the combined cloud population of Io's O and S neutral clouds and provides the dominant neutral population beyond the so-called ramp region at 7.4-7.8 RJ in the plasma torus. The calculated brightness of Europa's O cloud on the sky plane is very dim at the sub-Rayleigh level. The H₂ and O tori provide a new source of europagenic molecular and atomic pickup ions for the thermal plasma and introduce a neutral barrier in which new plasma sinks are created for the cooler iogenic plasma as it is transported radially outward and in which new sinks are created to alter the population and pitch angle distribution of the energetic plasma as it is transported radially inward. The europagenic instantaneous pickup ion rates are peaked at Europa's orbit, dominate the iogenic pickup ion rates beyond the ramp region, and introduce new secondary plasma source peaks in the solution of the plasma transport problem. The H₂ torus is identified as the unknown Europa gas torus that creates both the observed loss of energetic H⁺ ions at Europa's orbit and the corresponding measured ENA production rate for H.     

Cassini spectra and photometry 0.25-5.1 μm of the small inner satellites of Saturn

The nominal tour of the Cassini mission enabled the first spectra and solar phase curves of the small inner satellites of Saturn. We present spectra from the Visual Infrared Mapping Spectrometer (VIMS) and the Imaging Science Subsystem (ISS) that span the 0.25-5.1 μm spectral range. The composition of Atlas, Pandora, Janus, Epimetheus, Calypso, and Telesto is primarily water ice, with a small amount (∼5%) of contaminant, which most likely consists of hydrocarbons. The optical properties of the “shepherd” satellites and the coorbitals are tied to the A-ring, while those of the Tethys Lagrangians are tied to the E-ring of Saturn. The color of the satellites becomes progressively bluer with distance from Saturn, presumably from the increased influence of the E-ring; Telesto is as blue as Enceladus. Janus and Epimetheus have very similar spectra, although the latter appears to have a thicker coating of ring material. For at least four of the satellites, we find evidence for the spectral line at 0.68 μm that Vilas et al. [Vilas, F., Larsen, S.M., Stockstill, K.R., Gaffley, M.J., 1996. Icarus 124, 262-267] attributed to hydrated iron minerals on Iapetus and Hyperion. However, it is difficult to produce a spectral mixing model that includes this component. We find no evidence for CO₂ on any of the small satellites. There was a sufficient excursion in solar phase angle to create solar phase curves for Janus and Telesto. They bear a close similarity to the solar phase curves of the medium-sized inner icy satellites. Preliminary spectral modeling suggests that the contaminant on these bodies is not the same as the exogenously placed low-albedo material on Iapetus, but is rather a native material. The lack of CO₂ on the small inner satellites also suggests that their low-albedo material is distinct from that on Iapetus, Phoebe, and Hyperion.     

Aegaeon (Saturn LIII), a G-ring object

Aegaeon (Saturn LIII, S/2008 S1) is a small satellite of Saturn that orbits within a bright arc of material near the inner edge of Saturn’s G-ring. This object was observed in 21 images with Cassini’s Narrow-Angle Camera between June 15 (DOY 166), 2007 and February 20 (DOY 051), 2009. If Aegaeon has similar surface scattering properties as other nearby small saturnian satellites (Pallene, Methone and Anthe), then its diameter is approximately 500 m. Orbit models based on numerical integrations of the full equations of motion show that Aegaeon’s orbital motion is strongly influenced by multiple resonances with Mimas. In particular, like the G-ring arc it inhabits, Aegaeon is trapped in the 7:6 corotation eccentricity resonance with Mimas. Aegaeon, Anthe and Methone therefore form a distinctive class of objects in the Saturn system: small moons in corotation eccentricity resonances with Mimas associated with arcs of debris. Comparisons among these different ring-arc systems reveal that Aegaeon’s orbit is closer to the exact resonance than Anthe’s and Methone’s orbits are. This could indicate that Aegaeon has undergone significant orbital evolution via its interactions with the other objects in its arc, which would be consistent with the evidence that Aegaeon’s mass is much smaller relative to the total mass in its arc than Anthe’s and Methone’s masses are.     

Ion and neutral sources and sinks within Saturn's inner magnetosphere: Cassini results

Using ion-electron fluid parameters derived from Cassini Plasma Spectrometer (CAPS) observations within Saturn's inner magnetosphere as presented in Sittler et al. [2006a. Cassini observations of Saturn's inner plasmasphere: Saturn orbit insertion results. Planet. Space Sci., 54, 1197-1210], one can estimate the ion total flux tube content, N_IONL², for protons, H⁺, and water group ions, W⁺, as a function of radial distance or dipole L shell. In Sittler et al. [2005. Preliminary results on Saturn's inner plasmasphere as observed by Cassini: comparison with Voyager. Geophys. Res. Lett. 32(14), L14S04), it was shown that protons and water group ions dominated the plasmasphere composition. Using the ion-electron fluid parameters as boundary condition for each L shell traversed by the Cassini spacecraft, we self-consistently solve for the ambipolar electric field and the ion distribution along each of those field lines. Temperature anisotropies from Voyager plasma observations are used with (T_⊥/T_∥)_W⁺∼5 and (T_⊥/T_∥)_H⁺∼2. The radio and plasma wave science (RPWS) electron density observations from previous publications are used to indirectly confirm usage of the above temperature anisotropies for water group ions and protons. In the case of electrons we assume they are isotropic due to their short scattering time scales. When the above is done, our calculation show N_IONL² for H⁺ and W⁺ peaking near Dione's L shell with values similar to that found from Voyager plasma observations. We are able to show that water molecules are the dominant source of ions within Saturn's inner magnetosphere. We estimate the ion production rate S_ION∼10²⁷ ions/s as function of dipole L using N_H⁺, N_W⁺ and the time scale for ion loss due to radial transport τ_D and ion-electron recombination τ_REC. The ion production shows localized peaks near the L shells of Tethys, Dione and Rhea, but not Enceladus. We then estimate the neutral production rate, S_W, from our ion production rate, S_ION, and the time scale for loss of neutrals by ionization, τ_ION, and charge exchange, τ_CH. The estimated source rate for water molecules shows a pronounced peak near Enceladus’ L shell L∼4, with a value S_W∼2×10²⁸ mol/s.     

Revised vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ Northern Spring Equinox from 2006 to 2008: Application of new methane absorption data and comparison with Neptune

Long-slit spectroscopy observations of Uranus by the United Kingdom InfraRed Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s Northern Spring Equinox in December 2007.These spectra were analysed and presented by Irwin et al. (Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2009]. Icarus 203, 287-302), but since publication, a new set of methane absorption data has become available (Karkoschka, E., Tomasko, M. [2010]. Methane absorption coefficients for the jovian planets from laboratory, Huygens, and HST data. Icarus 205, 674-694.), which appears to be more reliable at the cold temperatures and high pressures of Uranus’ deep atmosphere. We have fitted k-coefficients to these new methane absorption data and we find that although the latitudinal variation and inter-annual changes reported by Irwin et al. (2009) stand, the new k-data place the main cloud deck at lower pressures (2-3 bars) than derived previously in the H-band of ∼3-4 bars and ∼3 bars compared with ∼6 bars in the J-band. Indeed, we find that using the new k-data it is possible to reproduce satisfactorily the entire observed centre-of-disc Uranus spectrum from 1 to 1.75 μm with a single cloud at 2-3 bars provided that we make the particles more back-scattering at wavelengths less than 1.2 μm by, for example, increasing the assumed single-scattering albedo from 0.75 (assumed in the J and H-bands) to near 1.0. In addition, we find that using a deep methane mole fraction of 4% in combination with the associated warm ‘F’ temperature profile of Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92, 14987-15001), the retrieved cloud deck using the new (Karkoschka and Tomasko, 2010) methane absorption data moves to between 1 and 2 bars.The same methane absorption data and retrieval algorithm were applied to observations of Neptune made during the same programme and we find that we can again fit the entire 1-1.75 μm centre-of-disc spectrum with a single cloud model, providing that we make the stratospheric haze particles (of much greater opacity than for Uranus) conservatively scattering (i.e. ω = 1) and we also make the deeper cloud particles, again at around the 2 bar level more reflective for wavelengths less than 1.2 μm. Hence, apart from the increased opacity of stratospheric hazes in Neptune’s atmosphere, the deeper cloud structure and cloud composition of Uranus and Neptune would appear to be very similar.     

Uranus’ cloud structure and seasonal variability from Gemini-North and UKIRT observations

Observations of Uranus were made in September 2009 with the Gemini-North telescope in Hawaii, using both the NIFS and NIRI instruments. Observations were acquired in Adaptive Optics mode and have a spatial resolution of approximately 0.1″.NIRI images were recorded with three spectral filters to constrain the overall appearance of the planet: J, H-continuum and CH₄(long), and long slit spectroscopy measurements were also made (1.49-1.79 μm) with the entrance slit aligned on Uranus’ central meridian. To acquire spectra from other points on the planet, the NIFS instrument was used and its 3″ × 3″ field of view stepped across Uranus’ disc. These observations were combined to yield complete images of Uranus at 2040 wavelengths between 1.476 and 1.803 μm.The observed spectra along Uranus central meridian were analysed with the NEMESIS retrieval tool and used to infer the vertical/latitudinal variation in cloud optical depth. We find that the 2009 Gemini data perfectly complement our observations/conclusions from UKIRT/UIST observations made in 2006-2008 and show that the north polar zone at 45°N has continued to steadily brighten while that at 45°S has continued to fade. The improved spatial resolution of the Gemini observations compared with the non-AO UKIRT/UIST data removes some of the earlier ambiguities with our previous analyses and shows that the opacity of clouds deeper than the 2-bar level does indeed diminish towards the poles and also reveals a darkening of the deeper cloud deck near the equator, perhaps coinciding with a region of subduction. We find that the clouds at 45°N,S lie at slightly lower pressures than the clouds at more equatorial latitudes, which suggests that they might possibly be composed of a different condensate, presumably CH₄ ice, rather than H₂S or NH₃ ice, which is assumed for the deeper cloud. In addition, analysis of the centre-to-limb curves of both the Gemini/NIFS and earlier UKIRT/UIST IFU observations shows that the main cloud deck has a well-defined top, and also allows us to better constrain the particle scattering properties.Overall, Uranus appeared to be less convectively active in 2009 than in the previous 3 years, which suggests that now the northern spring equinox (which occurred in 2007) is passed the atmosphere is settling back into the quiescent state seen by Voyager 2 in 1986. However, a number of discrete clouds were still observed, with one at 15°N found to lie near the 500 mb level, while another at 30°N, was seen to be much higher at near the 200 mb level. Such high clouds are assumed to be composed of CH₄ ice.     

The absence of endogenic methane on Titan and its implications for the origin of atmospheric nitrogen

We calculate the D/H ratio of CH₄ from serpentinization on Titan to determine whether Titan’s atmospheric CH₄ was originally produced inside the giant satellite. This is done by performing equilibrium isotopic fractionation calculations in the CH₄-H₂O-H₂ system, with the assumption that the bulk D/H ratio of the system is equivalent to that of the H₂O in the plume of Enceladus. These calculations show that the D/H ratio of hydrothermally produced CH₄ would be markedly higher than that of atmospheric CH₄ on Titan. The implication is that Titan’s CH₄ is a primordial chemical species that was accreted by the moon during its formation. There are two evolutionary scenarios that are consistent with the apparent absence of endogenic CH₄ in Titan’s atmosphere. The first is that hydrothermal systems capable of making CH₄ never existed on Titan because Titan’s interior has always been too cold. The second is that hydrothermal systems on Titan were sufficiently oxidized so that C existed in them predominately in the form of CO₂. The latter scenario naturally predicts the formation of endogenic N₂, providing a new hypothesis for the origin of Titan’s atmospheric N₂: the hydrothermal oxidation of ¹⁵N-enriched NH₃. A primordial origin for CH₄ and an endogenic origin for N₂ are self-consistent, but both hypotheses need to be tested further by acquiring isotopic data, especially the D/H ratio of CH₄ in comets, and the ¹⁵N/¹⁴N ratio of NH₃ in comets and that of N₂ in one of Enceladus’ plumes.     

A comprehensive numerical simulation of Io’s sublimation-driven atmosphere

Io’s sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo (DSMC) method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma heating, planetary rotation, inhomogeneous surface frost, molecular residence time of SO₂ on the exposed (non-volatile) rocky surface, and surface temperature distribution are investigated. Circumplanetary flow is predicted to develop from the warm dayside toward the cooler nightside. Io’s rotation leads to a highly asymmetric frost surface temperature distribution (due to the frost’s high thermal inertia) which results in circumplanetary flow that is not axi-symmetric about the subsolar point. The non-equilibrium thermal structure of the atmosphere, specifically vibrational and rotational temperatures, is also examined. Plasma heating is found to significantly inflate the atmosphere on both the dayside and nightside. The plasma energy flux causes high temperatures at high altitudes but plasma energy depletion through the dense gas column above the warmest frost permits gas temperatures cooler than the surface at low altitudes. A frost map (Douté, S., Schmitt, B., Lopes-Gautier, R., Carlson, R., Soderblom, L., Shirley, J., and the Galileo NIMS Team [2001]. Icarus 149, 107-132) is used to control the sublimated flux of SO₂ which can result in inhomogeneous column densities that vary by nearly a factor of four for the same surface temperature. A short residence time for SO₂ molecules on the “rock” component is found to smooth lateral atmospheric inhomogeneities caused by variations in the surface frost distribution, creating an atmosphere that looks nearly identical to one with uniform frost coverage. A longer residence time is found to agree better with mid-infrared observations (Spencer, J.R., Lellouch, E., Richter, M.J., López-Valverde, M.A., Jessup, K.L, Greathouse, T.K., Flaud, J. [2005]. Icarus 176, 283-304) and reproduce the observed anti-jovian/sub-jovian column density asymmetry. The computed peak dayside column density for Io assuming a surface frost temperature of 115 K agrees with those suggested by Lyman-α observations (Feaga, L.M., McGrath, M., Feldman, P.D. [2009]. Icarus 201, 570-584). On the other hand, the peak dayside column density at 120 K is a factor of five larger and is higher than the upper range of observations (Jessup, K.L., Spencer, J.R., Ballester, G.E., Howell, R.R., Roesler, F., Vigel, M., Yelle, R. [2004]. Icarus 169, 197-215; Spencer et al., 2005).