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.     

Release of N2, CH4, CO2, and H2O from surface ices on Enceladus

We vapor deposit at 20 K a mixture of gases with the specific Enceladus plume composition measured in situ by the Cassini INMS [Waite, J.H., Combi, M.R., Ip, W.H., Cravens, T.E., McNutt, R.L., Kasprzak, W., Yelle, R., Luhmann, J., Niemann, H., Gell, D., Magee, B., Fletcher, G., Lunine, J., Tseng, W.L., 2006. Science 311, 1419-1422] to form a mixed molecular ice. As the sample is slowly warmed, we monitor the escaping gas quantity and composition with a mass spectrometer. Pioneering studies [Schmitt, B., Klinger, J., 1987. Different trapping mechanisms of gases by water ice and their relevance for comet nuclei. In: Rolfe, E.J., Battrick, B. (Eds.), Diversity and Similarity of Comets. SP-278. ESA, Noordwijk, The Netherlands, pp. 613-619; Bar-Nun, A., Kleinfeld, I., Kochavi, E., 1988. Phys. Rev. B 38, 7749-7754; Bar-Nun, A., Kleinfeld, I., 1989. Icarus 80, 243-253] have shown that significant quantities of volatile gases can be trapped in a water ice matrix well above the temperature at which the pure volatile ice would sublime. For our Enceladus ice mixture, a composition of escaping gases similar to that detected by Cassini in the Enceladus plume can be generated by the sublimation of the H₂O:CO₂:CH₄:N₂ mixture at temperatures between 135 and 155 K, comparable to the high temperatures inferred from the CIRS measurements [Spencer, J.R., Pearl, J.C., Segura, M., Flasar, F.M., Mamoutkine, A., Romani, P., Buratti, B.J., Hendrix, A.R., Spilker, L.J., Lopes, R.M.C., 2006. Science 311, 1401-1405] of the Enceladus “tiger stripes.” This suggests that the gas escape phenomena that we measure in our experiments are an important process contributing to the gases emitted from Enceladus. A similar experiment for ice deposited at 70 K shows that both the processes of volatile trapping and release are temperature dependent over the temperature range relevant to Enceladus.     

A model for the temperature-dependence of tidal dissipation in convective plumes on icy satellites: Implications for Europa and Enceladus

To explain the formation of surface features on Europa, Enceladus, and other satellites, many authors have postulated the spatial localization of tidal heating within convective plumes. However, the concept that enhanced tidal heating can occur within a convective plume has not been rigorously tested. Most models of this phenomenon adopt a tidal heating with a temperature-dependence derived for an incompressible, homogeneous (zero-dimensional) Maxwell material, but it is unclear whether this formulation is relevant to the heterogeneous situation of a warm plume surrounded by cold ice. To determine whether concentrated dissipation can occur in convective plumes, we develop a two-dimensional model to compute the volumetric dissipation rate for an idealized, vertically oriented, isolated convective plume obeying a Maxwellian viscoelastic compressible rheology. We apply the model to the Europa and Enceladus ice shells, and we investigate the consequences for partial melting and resurfacing processes on these bodies. We find that the tidal heating is strongly temperature dependent in a convective ice plume and could produce elevated temperatures and local partial melting in the ice shells of Europa and Enceladus. Our calculation provides the first quantitative verification of the hypothesis by Sotin et al. [Sotin, C., Head, J.W., Tobie, G., 2002. Geophys. Res. Lett. 29. 74-1] and others that the tidal dissipation rate is a strong function of temperature inside a convective plume. On Europa, such localized heating could help allow the formation of domes and chaos terrains by convection. On Enceladus, localized tidal heating in a thermal plume could explain the concentrated activity at the south pole and its associated heat transport of 2-7 GW.     

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.     

Search for and limits on plume activity on Mimas, Tethys, and Dione with the Cassini Visual Infrared Mapping Spectrometer (VIMS)

Cassini Visual Infrared Mapping Spectrometer (VIMS) observations of Mimas, Tethys, and Dione obtained during the nominal and extended missions at large solar phase angles were analyzed to search for plume activity. No forward scattered peaks in the solar phase curves of these satellites were detected. The upper limit on water vapor production for Mimas and Tethys is one order of magnitude less than the production for Enceladus. For Dione, the upper limit is two orders of magnitude less, suggesting this world is as inert as Rhea (Pitman, K.M., Buratti, B.J., Mosher, J.A., Bauer, J.M., Momary, T., Brown, R.H., Nicholson, P.D., Hedman, M.M. [2008]. Astrophys. J. Lett. 680, L65-L68). Although the plumes are best seen at ∼2.0 μm, Imaging Science Subsystem (ISS) Narrow Angle Camera images obtained at the same time as the VIMS data were also inspected for these features. None of the Cassini ISS images shows evidence for plumes. The absence of evidence for any Enceladus-like plumes on the medium-sized saturnian satellites cannot absolutely rule out current geologic activity. The activity may below our threshold of detection, or it may be occurring but not captured on the handful of observations at large solar phase angles obtained for each moon. Many VIMS and ISS images of Enceladus at large solar phase angles, for example, do not contain plumes, as the active “tiger stripes” in the south pole region are pointed away from the spacecraft at these times. The 7-year Cassini Solstice Mission is scheduled to gather additional measurements at large solar phase angles that are capable of revealing activity on the saturnian moons.     

Appearance of Saturn’s F ring azimuthal channels for the anti-alignment configuration between the ring and Prometheus

In this article we explore the aspect of the F ring with respect to the anti-alignment configuration between the ring and Prometheus. We focus our attention on the shape of the F ring’s azimuthal channels which were first reported by Porco et al. (Porco, C.C., Baker, E., Barbara, J., Beurle, K., Brahic, A., Burns, J.A., Charnoz, S., Cooper, N., Dawson, D.D., Del Genio, A.D., Denk, T., Dones, L., Dyudina, U., Evans, M.W., Giese, B., Grazier, K., Helfenstein, P., Ingersoll, A.P., Jacobson, R.A., Johnson, T.V., McEwen, A., Murray, C.D., Neukum, G., Owen, W.M., Perry, J., Roatsch, T., Spitale, J., Squyres, S., Thomas, P., Tiscareno, M., Turtle, E., Vasavada, A.R., Veverka, J., Wagner, R., West, R. [2005] Science, 307, 1226-1236) and numerically explored by Murray et al. (Murray, C.D., Chavez, C., Beurle, K., Cooper, N., Evans, M.W., Burns, J.A., Porco, C.C. [2005] Nature 437, 1326-1329) who found excellent agreement between Cassini’s ISS reprojected images and their numerical model via a direct comparison. We find that for anti-alignment the channels are wider and go deeper inside the ring material. From our numerical model we find a new feature, an island in the middle of the channel. This island is made up of the particles that have been perturbed the most by Prometheus and only appears when this satellite is close to apoapsis. In addition, plots of the anti-alignment configuration for different orbital stages of Prometheus are obtained and discussed here.     

Enceladus: The likely dominant nitrogen source in Saturn's magnetosphere

The spatial distribution of N⁺ in Saturn's magnetosphere obtained from Cassini Plasma Spectrometer (CAPS) data can be used to determine the spatial distribution and relative importance of the nitrogen sources for Saturn's magnetosphere. We first summarize CAPS data from 15 orbits showing the spatial and energy distribution of the nitrogen component of the plasma. This analysis re-enforces our earlier discovery [Smith, H.T., Shappirio, M., Sittler, E.C., Reisenfeld, D., Johnson, R.E., Baragiola, R.A., Crary, F.J., McComas, D.J., Young, D.T., 2005. Geophys. Res. Lett. 32 (14). L14S03] that Enceladus is likely the dominant nitrogen source for Saturn's inner magnetosphere. We also find a sharp enhancement in the nitrogen ion to water ion ratio near the orbit of Enceladus which, we show, is consistent with the presence of a narrow Enceladus torus as described in [Johnson, R.E., Liu, M., Sittler Jr., E.C., 2005. Geophys. Res. Lett. 32. L24201]. The CAPS data and the model described below indicate that N⁺ ions are a significant fraction of the plasma in this narrow torus. We then simulated the combined Enceladus and Titan nitrogen sources using the CAPS data as a constraint. This simulation is an extension of the model we employed earlier to describe the neutral tori produced by the loss of nitrogen from Titan [Smith, H.T., Johnson, R.E., Shematovich, V.I., 2004. Geophys. Res. Lett. 31 (16). L16804]. We show that Enceladus is the principal nitrogen source in the inner magnetosphere but Titan might account for a fraction of the observed nitrogen ions at the largest distances discussed. We also show that the CAPS data is consistent with Enceladus being a molecular nitrogen source with a nitrogen to water ratio roughly consistent with INMS [Waite, J.H., and 13 colleagues, 2006. Science 311 (5766), 1419-1422], but out-gassing of other nitrogen-containing species, such as ammonia, cannot be ruled out.     

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.     

Solid tidal friction above a liquid water reservoir as the origin of the south pole hotspot on Enceladus

Earth, Jupiter's moon Io and Saturn's tiny moon Enceladus are the only solid objects in the Solar System to be sufficiently geologically active for their internal heat to be detected by remote sensing. Interestingly, the endogenic activity on Enceladus is only located on a specific region at the south pole, from which jets of water vapor and ice particles have been observed [Spencer, J.R., and 9 colleagues, 2006. Science 311, 1401-1405; Porco, C.C., and 24 colleagues, 2006. Science 311, 1393-1401]. The current polar location of the thermal anomaly can possibly be explained by diapir-induced reorientation of the satellite [Nimmo, F., Pappalardo, R.T., 2006. Nature 441, 614-616], but the thermal anomaly triggering and the heat power required to sustain it over geological timescales remain problematic. Using a three-dimensional viscoelastic numerical model simulating the response of Enceladus to tidal forcing, we explore the effect of a low-viscosity anomaly in the ice shell, localized to the south polar region, on the tidal dissipation patterns. We demonstrate that only interior models with a liquid water layer at depth can explain the observed magnitude of dissipation rate and its particular location at the south pole. Moreover, we show that tidally-induced heat must be generated over a relatively broad region in the southern hemisphere, and it is then transferred toward the south pole where it is episodically released during relatively short resurfacing events. As large tidal dissipation and internal melting cannot be induced in the south polar region in the absence of a pre-existing liquid decoupling layer, we propose that liquid water must have been present in the interior for a very long period of time, and possibly since the satellite formation. Owing to the orbital equilibrium requirement [Meyer, J., Wisdom, J., 2007. Icarus 188, 535-539], sustaining some liquid water at depth is impossible if heat is continuously emitted at a rate of 4-8 GW at the south pole. Based on that requirement, we propose that the current thermal emission is not in equilibrium with the heat production, and that the thermal emission rate is abnormally high at present time. Alternatively, continuous dissipation at a rate of 1-2 GW in the ice shell at the south pole should be sufficient to induce internal melting and it could sustain a layer of liquid water at depth over geologic timescales.     

Negative ions in the Enceladus plume

During Cassini’s Enceladus encounter on 12th March 2008, the Cassini Electron Spectrometer, part of the CAPS instrument, detected fluxes of negative ions in the plumes from Enceladus. It is thought that these ions include negatively charged water group cluster ions associated with the plume and forming part of the ‘plume ionosphere’. In this paper we present our observations, argue that these are negative ions, and present preliminary mass identifications. We also suggest mechanisms for production and loss of the ions as constrained by the observations. Due to their short lifetime, we suggest that the ions are produced in or near the water vapour plume, or from the extended source of ice grains in the plume. We suggest that Enceladus now joins the Earth, Comet Halley and Titan as locations in the Solar System where negative ions have been directly observed although the ions observed in each case have distinctly different characteristics.     

Experimental and theoretical simulations of ice sublimation with implications for the chemical,isotopic, and physical evolution of icy objects

Recent ground-based and spacecraft observations of comets in the inner solar system reveal two distinct features that provide important insight into their origin and evolution. The first of these is the observation that the D/H ratio of water vapor in cometary comae is significantly higher than that in Vienna Standard Mean Ocean Water (VSMOW). The second observation is that cometary jets are bursty (i.e, roughly steady state emissions that are punctuated with short-lived outbursts of water vapor and other materials) (Hughes, 1990, Soderblom et al., 2002, Soderblom et al., 2004a, Soderblom et al., 2004b). We present an experimental and theoretical study of ice sublimation in a vacuum that reveals several heretofore unknown and fundamental characteristics about the kinetics and mechanisms of ice sublimation that may explain both of these observed phenomena. In particular, we observe quasi-periodic sublimation cascades on time scales of hours to days, the D/H ratio in the vapor issuing from the sample is in general different from that of the sample, and in many cases, quasi-periodic changes in the D/H ratio of the vapor accompany the sublimation cascades. Changes are also observed in the infrared spectrum of the sample before, during and after a cascade that are consistent with our hypothesis that vacuum sublimation of water ice is a diffusive process that works to leave behind the most strongly bound molecules. Finally, we speculate as to whether the effects observed in the lab can be extrapolated to cometary-nucleus-scale phenomena.     

Sublimation coefficient of water ice under simulated cometary-like conditions

In papers dealing with evolution of cometary nuclei it is commonly assumed that the coefficients of sublimation _s and condensation _c of vapour are both equal to one. The experimental investigation of ice samples under simulated cometary-like conditions (Kossacki, K.J., Kömle, N.I., Leliwa-Kopysty ski, J., Kargl, G., 1997. Thermal and structural evolution of cometary subsurface layer: selfconsistent model and experimental verification. Icarus 128, 127–144) suggests, however, that the sublimation flux calculated with the Hertz–Knudsen formula and the above assumption is nearly an order of magnitude too high. This may imply that actual values of _s for the ice/dust sample used in these experiments are of the order of 0.1. A similar conclusion can be drawn for _c from the results of various experiments concerning growth of ice crystals from the vapour phase and their sublimation (Lamb, D., Scott, W.D., 1972. Linear growth rates of ice crystals grown from the vapor phase. Journal of Crystal Growth 12, 21–31; Beckmann, W., Lacmann, R., 1982. Interface kinetics of growth and evaporation of ice II. Journal of Crystal Growth 58, 433–442; Sei, T., Gonda, T., 1989. The growth mechanism and the habit change of ice crystals growing from the vapour phase. Journal of Crystal Growth 94, 697–707). The exact values of both of these coefficients depend on various parameters such as temperature, concentration of surface impurities and deviation of the vapour pressure from that of the phase equilibrium. In this work the temperature dependence of the sublimation and condensation coefficients is discussed and an appropriate formula is proposed to fit the experimental results. This new formulation is then used to analyse the implications for the thermal conductivity of a porous cometary-like ice and the rate of vapour flux from a cometary nucleus.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

Characterizing organic particle impacts on inert metal surfaces: Foundations for capturing organic molecules during hypervelocity transits of Enceladus plumes

The presence and accessibility of a sub‐ice‐surface saline ocean at Enceladus, together with geothermal activity and a rocky core, make it a compelling location to conduct further, in‐depth, astrobiological investigations to probe for organic molecules indicative of extraterrestrial life. Cryovolcanic plumes in the south polar region of Enceladus enable the use of remote in situ sampling and analysis techniques. However, efficient plume sampling and the transportation of captured organic materials to an organic analyzer present unique challenges for an Enceladus mission. A systematic study, accelerating organic ice‐particle simulants into soft inert metal targets at velocities ranging 0.5–3.0 km s⁻¹, was carried out using a light gas gun to explore the efficacy of a plume capture instrument. Capture efficiency varied for different metal targets as a function of impact velocity and particle size. Importantly, organic chemical compounds remained chemically intact in particles captured at speeds up to ~2 km s⁻¹. Calibration plots relating the velocity, crater, and particle diameter were established to facilitate future ice‐particle impact experiments where the size of individual ice particles is unknown.     

High-resolution Enceladus atlas derived from Cassini-ISS images

The Cassini Imaging Science Subsystem (ISS) acquired 377 high-resolution images (<1 km/pixel) during three close flybys of Enceladus in 2005 [Porco, C.C., et al., 2006. Cassini observes the active south pole of Enceladus. Science 311, 1393-1401.]. We combined these images with lower resolution Cassini images and four others taken by Voyager cameras to produce a high-resolution global controlled mosaic of Enceladus. This global mosaic is the baseline for a high-resolution Enceladus atlas that consists of 15 tiles mapped at a scale of 1:500,000. The nomenclature used in this atlas was proposed by the Cassini imaging team and was approved by the International Astronomical Union (IAU). The whole atlas is available to the public through the Imaging Team's website (http://ciclops.org/maps).     

Efficiency of neutrino-induced radio measurements to inspect local areas of Enceladus

The Cassini probe observed a young and smooth surface around the south pole of Enceladus, while around the north pole the surface was found to be relatively old and inactive (Porco, C.C. et al. [2006]. Science 311, 1393–1401). This heterogeneous surface implies that the ice thickness of Enceladus is not uniform between the north and south polar regions. Determining the thickness of the icy layer is important to confirm the existence of an internal ocean as well as to reveal the heating mechanism of Enceladus. We show that the measurement of radio waves induced by cosmic neutrinos can be an effective method to constrain the ice thickness of a localized area where conventional gravity or electromagnetic field measurements cannot be used. This method could be used to constrain the thickness of the icy layer on Enceladus even if the ice is a few tens of kilometers thick, measuring over a period of several years, which greatly exceeds the ability of radar sounding, and hence could be used in future orbiter missions.     

Total particulate mass in Enceladus plumes and mass of Saturn’s E ring inferred from Cassini ISS images

The eclipse mosaic (PIA08329) of the Saturn system, taken on September 15, 2006 when Cassini was in Saturn’s shadow, contains numerous color images of the Enceladus plume and the E ring at phase angles ranging from 173° to 179°. These forward-scattering observations sample the diffraction peak for particle radii in the 1–5 μm range. The phase angle dependence and total brightness are sensitive indicators of the total mass of solid material in the plume. We fit the data with a variety of particle shapes and size distributions, and find that the median radius of the equivalent-volume sphere is 3.1 μm, with an uncertainty of ±0.5 μm. The total mass of particles in the plume is (1.45 ± 0.5) × 10⁵ kg. We have not considered variations with altitude in the particle size and shape distribution, and we leave that for another paper. We find that the brightness of the E ring varies with position in the orbit, not only because of the viewing geometry, e.g., variations in phase angle, but also because of some unknown intrinsic variability. The total mass of solid material in the E ring is (12 ± 5.5) × 10⁸ kg. For the plume, the production rate of particles – the mass per unit time leaving the vents is 51 ± 18 kg s⁻¹. We estimate that 9% of these particles are escaping from Enceladus, implying lifetimes of ∼8 years for the E ring particles. Based on three comparisons with vapor amounts from ultraviolet spectroscopy, the ice/vapor ratio is in the range 0.35–0.70. This high ratio poses a problem for theories in which particles form by condensation from the gas phase, and could indicate that particles are formed as spray from a liquid reservoir.     

On the origin of south polar folds on Enceladus

Recent high-resolution Cassini images of the south polar terrain of Enceladus reveal regions of short-wavelength deformation, inferred to be compressional folds between the Baghdad and Damascus tiger stripes (Spencer, J.R., Barr, A.C., Esposito, L.W., Helfenstein, P., Ingersoll, A.P., Jaumann, R., McKay, C.P., Nimmo, F., Waite, J.H. [2009a]. Enceladus: An active cryovolcanic satellite. In: Saturn after Cassini-Huygens. Springer, New York, pp. 683-722). Here, we use Fourier analysis of the bright/dark variations to show that the folds have a dominant wavelength of 1.1 ± 0.4 km. We use the simple model of lava flow folding from Fink (Fink, J. [1980]. Geology 8, 250-254) to show that the folds could form in an ice shell with an upper high-viscosity boundary layer of thickness <400 m, with a driving stress of 40-80 kPa, and strain rate between 10⁻¹⁴ s⁻¹ and 10⁻¹² s⁻¹. Such deformation rates imply resurfacing of the SPT in 0.05-5 Myr, consistent with its estimated surface age. Measurements of fold topography and more sophisticated numerical modeling can narrow down the conditions of fold formation and provide valuable constraints on the thermal structure of the ice shell on Enceladus.     

Monte Carlo simulations of the water vapor plumes on Enceladus

Monte Carlo simulations are used to model the July 14, 2005 UVIS stellar occultation observations of the water vapor plumes on Enceladus. These simulations indicate that the observations can be best fit if the water molecules ejected along the Tiger Stripes in the South Polar region of Enceladus have a vertical surface velocity of 300-500 m/s at the surface. The high surface velocity suggests that the plumes on Enceladus originate from some depth beneath the surface. The total escape rate of water molecules is 4-6×¹⁰²⁷ s⁻¹, or 120-180 kg/s, consistent with previous works, and more than 100 times the estimated mass escape rate for ice particles. The average deposition rate in the South Polar region is on the order of 10¹¹ cm⁻² s⁻¹, yielding a resurfacing rate as high as 3×¹⁰⁻⁴ cm/yr. The globally averaged deposition rate of water molecules is about one order of magnitude lower.     

Photochemistry of methane-water ices

We report a study on the broadband ultraviolet photolysis of methane-water ice mixtures, at low methane concentrations and temperatures relevant to the icy satellites of the outer Solar System. The photochemistry of these mixtures is dominated by the action of hydroxyl radicals on methane and the resulting products. This implies that, given sufficient exposure time, the methane will eventually be completely oxidized to carbon dioxide. The presence of methane inhibits the formation of hydrogen peroxide by serving as a trap for hydroxyl radicals. The distribution of photochemical products is broadly similar to that previously conducted using ion and electron sources, with some differences possibly attributable to the difference in radiation source. The results are applicable to a variety of icy bodies in the Solar System. On Enceladus, where methane mixed with water is measured in the plumes, methane in the surface ices is subject to oxidation and will eventually be converted to CO₂. The CH stretch feature detected in the VIMS spectra of the Enceladus surface ice suggests that methane is currently being supplied to the surface ice, likely from re-condensation of the plume gas.