Cassini RPWS observations of dust in Saturn's E Ring

The Cassini radio and plasma wave science (RPWS) instrument is sensitive to few-micron dust grains impacting on the spacecraft at relative speeds of order 10 km/s. Through the first year or so of operations in orbit at Saturn, the RPWS has made a number of both inclined and equatorial crossings of the E ring, particularly near the orbit of Enceladus. Assuming water ice grains, the typical size particle detected by the RPWS has a radius of a few microns. Peak impact rates of about 50 s⁻¹ are found near the orbit of Enceladus corresponding to densities of order 5×10⁻⁴ m⁻³. The variation of dust fluxes as a function of height above or below the equator is well described by a Gaussian distribution with a scale height of about 2800 km although there is usually some non-Gaussian variation near the peak fluxes suggesting some structure in the core of the ring. Offsets of the peak number densities are typically of the order of a few hundred km from the geometric equator. A near-equatorial radial profile through the orbit of Enceladus shows a sharply peaked distribution at the orbit of the moon. A size distribution averaged over several passes through the orbit of Enceladus is determined which varies as m^−2.80. The peak in dust number density at the orbit of Enceladus is consistent with previous optical measurements and strongly supports the suggestion that Enceladus is a primary source for E ring particles.     

Deconvoluting measurement uncertainty from the meteor speed distribution

Debiasing the velocity distribution of meteors observed by the Canadian Meteor Orbit Radar (CMOR) yields a distribution with large numbers of slow meteors. The distribution also contains significant numbers of hyperbolic meteors, in conflict with the expectation that interstellar meteors should be rare. In Moorhead et al. (2017a), we noted that measurement uncertainties were possibly smoothing the speed distribution and redistributing meteors to the extreme ends of the speed distribution. In this report, we use techniques analogous to image sharpening to remove the blurring caused by measurement uncertainties. The deconvolved speed distribution appears to have no meteors slower than 14 km s⁻¹ and none faster than 74 km s⁻¹. The result is to substantially raise the characteristic velocity of incoming meteoroids from 12.9 to 20.0 km s⁻¹.     

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.     

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.     

Decreased values of cosmic dust number density estimates in the Solar System

Experiments to investigate the effect of impacts on side-walls of dust detectors such as the present NASA/ESA Galileo/Ulysses instrument are reported. Side walls constitute 27% of the internal area of these instruments, and increase field of view from 140° to 180°. Impact of cosmic dust particles onto Galileo/Ulysses Al side walls was simulated by firing Fe particles, 0.5-5 μm diameter, 2-50 km s⁻¹, onto an Al plate, simulating the targets of Galileo and Ulysses dust instruments. Since side wall impacts affect the rise time of the target ionization signal, the degree to which particle fluxes are overestimated varies with velocity. Side-wall impacts at particle velocities of 2-20 km s⁻¹ yield rise times 10-30% longer than for direct impacts, so that derived impact velocity is reduced by a factor of ∼2. Impacts on side wall at 20-50 km s⁻¹ reduced rise times by a factor of ∼10 relative to direct impact data. This would result in serious overestimates of flux of particles intersecting the dust instrument at velocities of 20-50 km s⁻¹. Taking into account differences in laboratory calibration geometry we obtain the following percentages for previous overestimates of incident particle number density values from the Galileo instrument [Grün et al., 1992. The Galileo dust detector. Space Sci. Rev. 60, 317-340]: 55% for 2 km s⁻¹ impacts, 27% at 10 km s⁻¹ and 400% at 70 km s⁻¹. We predict that individual particle masses are overestimated by ∼10-90% when side-wall impacts occur at 2-20 km s⁻¹, and underestimated by ∼10-¹⁰² at 20-50 km s⁻¹. We predict that wall impacts at 20-50 km s⁻¹ can be identified in Galileo instrument data on account of their unusually short target rise times. The side-wall calibration is used to obtain new revised values [Krüger et al., 2000. A dust cloud of Ganymede maintained by hypervelocity impacts of interplanetary micrometeoroids. Planet. Space Sci. 48, 1457-1471; 2003. Impact-generated dust clouds surrounding the Galilean moons. Icarus 164, 170-187] of the Galilean satellite dust number densities of 9.4×¹⁰⁻⁵, 9.9×¹⁰⁻⁵, 4.1×¹⁰⁻⁵, and 6.8×¹⁰⁻⁵ m⁻³ at 1 satellite radius from Io, Europa, Ganymede, and Callisto, respectively. Additionally, interplanetary particle number densities detected by the Galileo mission are found to be 1.6×¹⁰⁻⁴, 7.9×¹⁰⁻⁴, 3.2×¹⁰⁻⁵, 3.2×¹⁰⁻⁵, and 7.9×¹⁰⁻⁴ m⁻³ at heliocentric distances of 0.7, 1, 2, 3, and 5 AU, respectively. Work by Burchell et al. [1999b. Acceleration of conducting polymer-coated latex particles as projectiles in hypervelocity impact experiments. J. Phys. D: Appl. Phys. 32, 1719-1728] suggests that low-density “fluffy” particles encountered by Ulysses will not significantly affect our results—further calibration would be useful to confirm this.     

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.     

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.     

Effect of target properties and impact velocity on ejection dynamics and ejecta deposition

Impact craters are formed by the displacement and ejection of target material. Ejection angles and speeds during the excavation process depend on specific target properties. In order to quantify the influence of the constitutive properties of the target and impact velocity on ejection trajectories, we present the results of a systematic numerical parameter study. We have carried out a suite of numerical simulations of impact scenarios with different coefficients of friction (0.0–1.0), porosities (0–42%), and cohesions (0–150 MPa). Furthermore, simulations with varying pairs of impact velocity (1–20 km s⁻¹) and projectile mass yielding craters of approximately equal volume are examined. We record ejection speed, ejection angle, and the mass of ejected material to determine parameters in scaling relationships, and to calculate the thickness of deposited ejecta by assuming analytical parabolic trajectories under Earth gravity. For the resulting deposits, we parameterize the thickness as a function of radial distance by a power law. We find that strength—that is, the coefficient of friction and target cohesion—has the strongest effect on the distribution of ejecta. In contrast, ejecta thickness as a function of distance is very similar for different target porosities and for varying impact velocities larger than ~6 km s⁻¹. We compare the derived ejecta deposits with observations from natural craters and experiments.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

Distribution of icy particles across Enceladus' surface as derived from Cassini-VIMS measurements

The surface of Enceladus consists almost completely of water ice. As the band depths of water ice absorptions are sensitive to the size of particles, absorptions can be used to map variations of icy particles across the surface. The Visual and Infrared Mapping Spectrometer (VIMS) observed Enceladus with a high spatial resolution during three Cassini flybys in 2005 (orbits EN 003, EN 004 and EN 011). Based on these data we measured the band depths of water ice absorptions at 1.04, 1.25, 1.5, and 2 μm. These band depths were compared to water ice models that represent theoretically calculated reflectance spectra for a range of particle diameters between 2 μm and 1 mm. The agreement between the experimental (VIMS) and model values supports the assumption that pure water ice characterizes the surface of Enceladus and therefore that variations in band depth correspond to variations in water ice particle diameters. Our measurements show that the particle diameter of water ice increases toward younger tectonically altered surface units with the largest particles exposed in relatively “fresh” surface material. The smallest particles were generally found in old densely cratered terrains. The largest particles (∼0.2 mm) are concentrated in the so called “tiger stripes” at the south pole. In general, the particle diameters are strongly correlated with geologic features and surface ages, indicating a stratigraphic evolution of the surface that is caused by cryovolcanic resurfacing and impact gardening.     

Unconfined shock experiments: A pilot study into the shock‐induced melting and devolatilization of calcite

We shocked calcite in an unconfined environment by launching small marble cylinders at 0.8–5.5 km s^?1 into aluminum or copper plates, producing shock stresses between 5 and 79 GPa. The resulting 5–20 mm craters contained intimately mixed clastic and molten projectile residues over the entire pressure range, with melting commencing already at 5 GPa. Stoichiometrically pure calcite melts were not observed as all melts contained target metal. Some of these residues were distinctly depleted in CO₂ and some contained even tiny CaO crystals, thus illustrating partial to complete loss of CO₂. We interpret a thin seam of finely crystalline calcite to be the product of back reactions between CaO and CO₂. The amount of carbonate residue in these craters, especially those at low velocities (<2 km s^?1), is dramatically less than that of silicate impactors in similar cratering experiments, and we suggest that this is due to substantial outgassing of CO₂. Similarly, the volume of carbonate melts relative to the volume of limestone or dolomite in many terrestrial crater structures seems insignificant as well, as is the volume of carbonate melt compared to the volume of impact melts derived from silicates. These volume considerations suggest that volatilization of CO₂ is the dominant process in carbonate‐containing targets. Because we have difficulties in explaining naturally occurring calcite melts by shock processes in dolomite‐dominated targets, we speculate—essentially via process of elimination—that such carbonate melt blebs might be condensation products from an impact‐produced vapor cloud.     

Extent of thermal ablation suffered by model organic microparticles during aerogel capture at hypervelocities

Abstract— New model organic microparticles are used to assess the thermal ablation that occurs during aerogel capture at speeds from 1 to 6 km s^?1. Commercial polystyrene particles (20 μm diameter) were coated with an ultrathin 20 nm overlayer of an organic conducting polymer, polypyrrole. This overlayer comprises only 0.8% by mass of the projectile but has a very strong Raman signature, hence its survival or destruction is a sensitive measure of the extent of chemical degradation suffered. After aerogel capture, microparticles were located via optical microscopy and their composition was analyzed in situ using Raman microscopy. The ultrathin polypyrrole overlayer survived essentially intact for impacts at ～1 km s^?1, but significant surface carbonization was found at 2 km s^?1, and major particle mass loss at ≥3 km s^?1. Particles impacting at ～6.1 km s^?1 (the speed at which cometary dust was collected in the NASA Stardust mission) were reduced to approximately half their original diameter during aerogel capture (i.e., a mass loss of 84%). Thus significant thermal ablation occurs at speeds above a few km s^?1. This suggests that during the Stardust mission the thermal history of the terminal dust grains during capture in aerogel may be sufficient to cause significant processing or loss of organic materials. Further, while Raman D and G bands of carbon can be obtained from captured grains, they may well reflect the thermal processing during capture rather than the pre‐impact particle's thermal history.     

Allan Sandage and the cosmic expansion

This is an account of Allan Sandage’s work on (1) The character of the expansion field. For many years he has been the strongest defender of an expanding Universe. He later explained the CMB dipole by a local velocity of 220±50 km s⁻¹ toward the Virgo cluster and by a bulk motion of the Local supercluster (extending out to ∼3500 km s⁻¹) of 450–500 km s⁻¹ toward an apex at l=275, b=12. Allowing for these streaming velocities he found linear expansion to hold down to local scales (∼300 km s⁻¹). (2) The calibration of the Hubble constant. Probing different methods he finally adopted—from Cepheid-calibrated SNe Ia and from independent RR Lyr-calibrated TRGBs—H ₀=62.3±1.3±5.0 km s⁻¹ Mpc⁻¹.     

The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

The E-ring in the vicinity of Enceladus: II. Probing the moon's interior—The composition of E-ring particles

The population of Saturn's outermost tenuous E-ring is dominated by tiny water ice particles. Active volcanism on the moon Enceladus, embedded in the E-ring, has since late 2005 been known to be a major source of particles replenishing the ring. Therefore particles in the vicinity of Enceladus may provide crucial information about the dynamical and chemical processes occurring below the moon's icy surface. Here we present a statistical evaluation of more than 2000 impact ionisation mass spectra of Saturn's E-ring particles, with sizes predominantly below 1 μm, detected by the Cosmic Dust Analyser onboard the Cassini spacecraft. We focus on the identification of non-water features in spectra otherwise dominated by water ice signatures. Here we specify the categorisation of two different spectrum types, which probably represent two particle populations. Type I spectra imply pure water ice particles, whereas in Type II spectra organic compounds and/or silicate minerals are identified as impurities within the icy particles. This finding supports the hypothesis of a dynamic interaction of Enceladus' rocky core with liquid water.     

Experimental shock metamorphism of terrestrial basalts: Agglutinate-like particle formation,petrology, and magnetism

Hypervelocity impacts occur on bodies throughout our solar system, and play an important role in altering the mineralogy, texture, and magnetic properties in target rocks at nanometer to planetary scales. Here we present the results of hypervelocity impact experiments conducted using a two-stage light-gas gun with 5 mm spherical copper projectiles accelerated toward basalt targets with ~6 km s⁻¹ impact velocities. Four different types of magnetite- and titanomagnetite-bearing basalts were used as targets for seven independent experiments. These laboratory impacts resulted in the formation of agglutinate-like particles similar in texture to lunar agglutinates, which are an important fraction of lunar soil. Materials recovered from the impacts were examined using a suite of complementary techniques, including optical and scanning electron microscopy, micro-Raman spectroscopy, and high- and low-temperature magnetometry, to investigate the texture, chemistry, and magnetic properties of newly formed agglutinate-like particles and were compared to unshocked basaltic parent materials. The use of Cu-projectiles, rather than Fe- and Ni-projectiles, avoids magnetic contamination in the final shock products and enables a clearer view of the magnetic properties of impact-generated agglutinates. Agglutinate-like particles show shock features, such as melting and planar deformation features, and demonstrate shock-induced magnetic hardening (two- to seven-fold increases in the coercivity of remanence B_cr compared to the initial target materials) and decreases in low-field magnetic susceptibility and saturation magnetization.     

The role of sticky interstellar organic material in the formation of asteroids

Abstract— Collision experiments and measurements of viscoelastic properties were performed involving an interstellar organic material analogue to investigate the growth of organic grains in the protosolar nebula. The organic material was found to be stickiest at a radius of between 2.3 and 3.0 AU, with a maximum sticking velocity of 5 m s^?1 for millimeter‐size organic grains. This stickiness is considered to have resulted in the very rapid coagulation of organic grain aggregates and subsequent formation of planetesimals in the early stage of the turbulent accretion disk. The planetesimals formed in this region appear to be represent achondrite parent bodies. In contrast, the formation of planetesimals at <2.1 and >3.0 AU begins with the establishment of a passive disk because silicate and ice grains are not as sticky as organic grains.     

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.     

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.