Cryovolcanism on the icy satellites

Evidence of past cryovolcanism is widespread and extremely varied on the icy satellites. Some cryovolcanic landscapes, notably on Triton, are similar to many silicate volcanic terrains, including what appear to be volcanic rifts, calderas and solidified lava lakes, flow fields, breached cinder cones or stratovolcanoes, viscous lava domes, and sinuous rilles. Most other satellites have terrains that are different in the important respect that no obvious volcanoes are present. The preserved record of cryovolcanism generally is believed to have formed by eruptions of aqueous solutions and slurries. Even Triton's volcanic crust, which is covered by nitrogen-rich frost, is probably dominated by water ice. Nonpolar and weakly polar molecular liquids (mainly N₂, CH₄, CO, CO₂, and Ar), may originate by decomposition of gas-clathrate hydrates and may have been erupted on some icy satellites, but without water these substances do not form rigid solids that are stable against sublimation or melting over geologic time. Triton's plumes, active at the time of Voyager 2's flyby, may consist of multicomponent nonpolar gas mixtures. The plumes may be volcanogenic fumaroles or geyserlike emissions powered by deep internal heating, and, thus, the plumes may be indicating an interior that is still cryomagmatically active; or Triton's plumes may be powered by solar heating of translucent ices very near the surface. The Uranian and Neptunian satellites Miranda, Ariel, and Triton have flow deposits that are hundreds to thousands of meters thick (implying highly viscous lavas); by contrast, the Jovian and Saturnian satellites generally have plains-forming deposits composed of relatively thin flows whose thicknesses have not been resolved in Voyager images (thus implying relatively low-viscosity lavas). One possible explanation for this inferred rheological distinction involves a difference in volatile composition of the Uranian and Neptunian satellites on one hand and of the Jovian and Saturnian satellites on the other hand. Perhaps the Jovian and Saturnian satellites tend to have relatively "clean" compositions with water ice as the main volatile (ammonia and water-soluble salts may also be present). The Uranian and Neptunian satellites may possess large amounts of a chemically unequilibrated comet-like volatile assemblage, including methanol, formaldehyde, and a host of other highly water- and ammonia-water-soluble constituents and gas clathrate hydrates. These two volatile mixtures would produce melts that differ enormously in viscosity The geomorphologic similarity in the products of volcanism on Earth and Triton may arise partly from a rheological similarity of the ammonia-water-methanol series of liquids and the silicate series ranging from basalt to dacite. An abundance of gas clathrate hydrates hypothesized to be contained by the satellites of Uranus and Neptune could contribute to evidence of explosive volcanism on those objects.     

Radar observations of the icy Galilean satellites

We report 12.6-cm-wavelength radar observations of Europa, Ganymede, and Callisto made at the Arecibo Observatory in November 1977 and February 1979. When combined with previous observations, our results establish firmly the distinguishing radar properties of these satellites: (i) high geometric albedos, α; (ii) circular polarization ratios, μ_C, which anomalously exceed unity; (iii) linear polarization ratios, μ_L, which are approximately 0.5; and (iv) diffuse scattering which varies as cosⁿθ, where θ is angle of incidence and 1 ? n ? 2. We tabulate weighted-mean values of α, μ_C, μ_L, and n derived from observations between 1975 and 1979. The values of μ_C for Ganymede and Europa are nearly identical and significantly larger than that for Callisto. The values of n for Ganymede and Callisto are nearly identical and significantly smaller than that for Europa. Although significant albedo and/or polarization features are common in the radar spectra, the fractional rms fluctuation in disk-integrated properties is only ～10%. No time variation in the radar properties has been evident during 1976–1979.     

Large impact features on middle-sized icy satellites

All of the large impact features of the middle-sized icy satellites of Saturn and Uranus that were clearly observed by the Voyager spacecraft are described. New image mosaics and stereo-and-photoclinometrically-derived digital elevation models are presented. Landforms related to large impact features, such as secondary craters and possible antipodal effects are examined and evaluated. Of the large impacts, Odysseus on Tethys appears to have had the most profound effect on its “target” satellite of any of the impact features we examined. Our modeling suggests that the Odysseus impact may have caused the prompt formation of Ithaca Chasma, a belt of tectonic troughs that roughly follow a great circle normal to the center of Odysseus, although other hypotheses remain viable. We identify probable secondary cratering from Tirawa on Rhea. We attribute a number of converging coalescing crater chains on Rhea to a putative, possibly relatively fresh, ∼350 km-diameter impact feature. We examine the antipodes of Odysseus, the putative ∼350 km-diameter Rhean impact feature, and Tirawa, and conclude that evidence from Voyager data for damage from seismic focusing is equivocal, although our modeling results indicate that such damage may have occurred. We propose a number of observations and tests for Cassini that offer the opportunity to differentiate among the various explanations and speculations reviewed and evaluated in this study.     

Thermal- and plasma-induced molecular redistribution on the icy satellites

The molecular transport of condensed gas species across the surfaces of the icy satellites of Jupiter and Saturn is examined with the view of describing, in part, certain gross visual features associated with these bodies. Molecular redistribution induced by thermal sublimation and magnetospheric plasma-ion impact on satellites with negligible atmospheres is calculated by assuming that the molecules follow ballistic trajectories and by statistically selecting initial molecular velocities and points of origin. Erosion/deposition profiles so calculated are compared for a variety of satellite sizes and environments in order to understand the relative importance of sublimation and cold corotating plasma-ion- and fast plasma-ion-induced transport. The results are scaled to make them useful as new data is available for the icy satellites and their plasma environment. The erosion/deposition profiles are then used to discuss the appearance of a polar frost on Ganymede. A balance of magnetospheric-ion implantation and ion-induced molecular redistribution is used to discuss the observation of embedded SO₂ and the darkening of the trailing side on Europa. Ion-induced molecular transport may also limit the deposition of SO₂ frost in the polar regions of Io and may be a source of heavy particles in the Jovian and Saturnian magnetospheres.     

Experimental impact shock chemistry on planetary icy satellites

We present new experimental results on impact shock chemistry into icy satellites of the outer planets. Icy mixtures of pure water ice with CO₂, Na₂CO₃, CH₃OH, and CH₃OH/(NH₄)₂SO₄ at 77 K were ablated with a powerful pulsed laser—a new technique used to simulate shock processes which can occur during impacts. New products were identified by GC-MS and FTIR analyses after laser ablation. Our results show that hydrogen peroxide is formed in irradiated H₂O/CO₂ ices with a final concentration of 0.23%. CO and CH₃OH were also detected as main products. The laser ablation of frozen H₂O/Na₂CO₃ generates only CO and CO₂ as destruction products from the salt. Pulsed irradiation of water ice containing methanol leads also to the formation of CO and CO₂, generates methane and more complex molecules containing carbonyl groups like acetaldehyde, acetone, methyl formate, and a diether, dimethyl formal. The last three compounds are also produced when adding ammonium sulfate to H₂O/CH₃OH ice, but acetone is more abundant. The formation of two hydrocarbons, CH₄ and C₂H₆ is observed as well as the production of three nitrogen compounds, nitrous oxide, hydrogen cyanide, and acetonitrile.     

Modeling the morphological diversity of impact craters on icy satellites

Impact craters on icy satellites display a wide range of morphologies, some of which have no counterpart on rocky bodies. Numerical simulation studies have struggled to reproduce the diversity of features, such as central pits and transitions in crater depth with increasing diameter, observed on the icy Galilean satellites. The transitions in crater depth (at diameters of about 26 and 150 km on Ganymede and Callisto) have been interpreted as reflecting subsurface structure. Using the CTH shock physics code, we model the formation of craters with diameters between 400 m and about 200 km on Ganymede using different subsurface temperature profiles. Our calculations include recent improvements in the model equation of state for H₂O and quasi-static strength parameters for ice. We find that the shock-induced formation of dense high-pressure polymorphs (ices VI and VII) creates a gap in the crater excavation flow, which we call discontinuous excavation. For craters larger than about 20 km, discontinuous excavation concentrates a hot plug of material (>270 K and mostly on the melting curve) in the center of the crater floor. The size and occurrence of the hot plug are in good agreement with the observed characteristics of central pit craters, and we propose that a genetic link exists between them. We also derive depth versus diameter curves for different internal temperature profiles. In a 120 K isothermal crust, calculated craters larger than about 30 km diameter are deeper than observed and do not reproduce the transition at about 26 km diameter. Calculated crater depths are shallower and in good agreement with observations between about 30 and 150 km diameter using a warm thermal gradient representing a convective interior. Hence, the depth-to-diameter transition at about 26 km reflects thermal weakening of ice. Finally, simulation results generally support the hypothesis that the anomalous interior morphologies for craters larger than 100 km are related to the presence of a subsurface ocean.     

Shapes of the saturnian icy satellites and their significance

The sizes and shapes of six icy saturnian satellites have been measured from Cassini Imaging Science Subsystem (ISS) data, employing limb coordinates and stereogrammetric control points. Mimas, Enceladus, Tethys, Dione and Rhea are well described by triaxial ellipsoids; Iapetus is best represented by an oblate spheroid. All satellites appear to have approached relaxed, equilibrium shapes at some point in their evolution, but all support at least 300 m of global-wavelength topography. The shape of Enceladus is most consistent with a homogeneous interior. If Enceladus is differentiated, its shape and apparent relaxation require either lateral inhomogeneities in an icy mantle and/or an irregularly shaped core. Iapetus supports a fossil bulge of over 30 km, and provides a benchmark for impact modification of shapes after global relaxation. Satellites such as Mimas that have smoother limbs than Iapetus, and are expected to have higher impact rates, must have relaxed after the shape of Iapetus was frozen.     

Oceans in the icy Galilean satellites of Jupiter?

Equilibrium models of heat transfer by heat conduction and thermal convection show that the three satellites of Jupiter—Europa, Ganymede, and Callisto—may have internal oceans underneath ice shells tens of kilometers to more than a hundred kilometers thick. A wide range of rheology and heat transfer parameter values and present-day heat production rates have been considered. The rheology was cast in terms of a reference viscosity ν₀ calculated at the melting temperature and the rate of change A of viscosity with inverse homologous temperature. The temperature dependence of the thermal conductivity k of ice I has been taken into account by calculating the average conductivity along the temperature profile. Heating rates are based on a chondritic radiogenic heating rate of 4.5 pW kg⁻¹ but have been varied around this value over a wide range. The phase diagrams of H₂O (ice I) and H₂O + 5 wt% NH₃ ice have been considered. The ice I models are worst-case scenarios for the existence of a subsurface liquid water ocean because ice I has the highest possible melting temperature and the highest thermal conductivity of candidate ices and the assumption of equilibrium ignores the contribution to ice shell heating from deep interior cooling. In the context of ice I models, we find that Europa is the satellite most likely to have a subsurface liquid ocean. Even with radiogenic heating alone the ocean is tens of kilometers thick in the nominal model. If tidal heating is invoked, the ocean will be much thicker and the ice shell will be a few tens of kilometers thick. Ganymede and Callisto have frozen their oceans in the nominal ice I models, but since these models represent the worst-case scenario, it is conceivable that these satellites also have oceans at the present time. The most important factor working against the existence of subsurface oceans is contamination of the outer ice shell by rock. Rock increases the density and the pressure gradient and shifts the triple point of ice I to shallower depths where the temperature is likely to be lower then the triple point temperature. According to present knowledge of ice phase diagrams, ammonia produces one of the largest reductions of the melting temperature. If we assume a bulk concentration of 5 wt% ammonia we find that all the satellites have substantial oceans. For a model of Europa heated only by radiogenic decay, the ice shell will be a few tens of kilometers thinner than in the ice I case. The underlying rock mantle will limit the depth of the ocean to 80-100 km. For Ganymede and Callisto, the ice I shell on top of the H₂O-NH₃ ocean will be around 60- to 80-km thick and the oceans may be 200- to 350-km deep. Previous models have suggested that efficient convection in the ice will freeze any existing ocean. The present conclusions are different mainly because they are based on a parameterization of convective heat transport in fluids with strongly temperature dependent viscosity rather than a parameterization derived from constant-viscosity convection models. The present parameterization introduces a conductive stagnant lid at the expense of the thickness of the convecting sublayer, if the latter exists at all. The stagnant lid causes the temperature in the sublayer to be warmer than in a comparable constant-viscosity convecting layer. We have further modified the parameterization to account for the strong increase in homologous temperature, and therefore decrease in viscosity, with depth along an adiabat. This modification causes even thicker stagnant lids and further elevated temperatures in the well-mixed sublayer. It is the stagnant lid and the comparatively large temperature in the sublayer that frustrates ocean freezing.     

Mass-radius relationships in icy satellites after voyager

Using improved data for the masses and radii of the satellites of Jupiter and Saturn, models accounting for self-compression effects are presented for the interiors of Europa, Ganymede, Callisto, Rhea, and Titan. For the differentiated models, two different possible scenarios for heat transport are treated, as well as two different compositions for the silicate component. Undifferentiated models are also treated. In each case, the models of Ganymede, Callisto, and Titan show noticeable similarities. It is found that estimates of the ice-rock ratio are dependent upon the assumptions made about the heat transport mechanisms, the rock composition, and on the distribution of rock and ice in the satellite's interior.     

High-resolution CASSINI-VIMS mosaics of Titan and the icy Saturnian satellites

The Visual Infrared Mapping Spectrometer (VIMS) onboard the CASSINI spacecraft obtained new spectral data of the icy satellites of Saturn after its arrival at Saturn in June 2004. VIMS operates in a spectral range from 0.35 to 5.2 μm, generating image cubes in which each pixel represents a spectrum consisting of 352 contiguous wavebands.As an imaging spectrometer VIMS combines the characteristics of both a spectrometer and an imaging instrument. This makes it possible to analyze the spectrum of each pixel separately and to map the spectral characteristics spatially, which is important to study the relationships between spectral information and geological and geomorphologic surface features.The spatial analysis of the spectral data requires the determination of the exact geographic position of each pixel on the specific surface and that all 352 spectral elements of each pixel show the same region of the target. We developed a method to reproject each pixel geometrically and to convert the spectral data into map projected image cubes. This method can also be applied to mosaic different VIMS observations. Based on these mosaics, maps of the spectral properties for each Saturnian satellite can be derived and attributed to geographic positions as well as to geological and geomorphologic surface features. These map-projected mosaics are the basis for all further investigations.     

Impact cratering records of the mid-sized, icy saturnian satellites

Resolution of Voyager 1 and 2 images of the mid-sized, icy saturnian satellites was generally not much better than 1 km per line pair, except for a few, isolated higher resolution images. Therefore, analyses of impact crater distributions were generally limited to diameters (D) of tens of kilometers. Even with the limitation, however, these analyses demonstrated that studying impact crater distributions could expand understanding of the geology of the saturnian satellites and impact cratering in the outer Solar System. Thus to gain further insight into Saturn’s mid-sized satellites and impact cratering in the outer Solar System, we have compiled cratering records of these satellites using higher resolution CassiniISS images. Images from Cassini of the satellites range in resolution from tens m/pixel to hundreds m/pixel. These high-resolution images provide a look at the impact cratering records of these satellites never seen before, expanding the observable craters down to diameters of hundreds of meters. The diameters and locations of all observable craters are recorded for regions of Mimas, Tethys, Dione, Rhea, Iapetus, and Phoebe. These impact crater data are then analyzed and compared using cumulative, differential and relative (R) size-frequency distributions. Results indicate that the heavily cratered terrains on Rhea and Iapetus have similar distributions implying one common impactor population bombarded these two satellites. The distributions for Mimas and Dione, however, are different from Rhea and Iapetus, but are similar to one another, possibly implying another impactor population common to those two satellites. The difference between these two populations is a relative increase of craters with diameters between 10 and 30 km and a relative deficiency of craters with diameters between 30 and 80 km for Mimas and Dione compared with Rhea and Iapetus. This may support the result from Voyager images of two distinct impactor populations. One population was suggested to have a greater number of large impactors, most likely heliocentric comets (Saturn Population I in the Voyager literature), and the other a relative deficiency of large impactors and a greater number of small impactors, most likely planetocentric debris (Saturn Population II). Meanwhile, Tethys’ impact crater size-frequency distribution, which has some similarity to the distributions of Mimas, Dione, Rhea, and Iapetus, may be transitional between the two populations. Furthermore, when the impact crater distributions from these older cratered terrains are compared to younger ones like Dione’s smooth plains, the distributions have some similarities and differences. Therefore, it is uncertain whether the size-frequency distribution of the impactor population(s) changed over time. Finally, we find that Phoebe has a unique impact crater distribution. Phoebe appears to be lacking craters in a narrow diameter range around 1 km. The explanation for this confined “dip” at D = 1 km is not yet clear, but may have something to do with the interaction of Saturn’s irregular satellites or the capture of Phoebe.     

Stagnant lid convection in the mid-sized icy satellites of Saturn

Thermal history models for the mid-sized saturnian satellites Mimas, Tethys, Dione, Iapetus, and Rhea have been calculated assuming stagnant lid convection in undifferentiated satellites and varying parameter values over broad ranges. Of all five satellites under consideration, only Dione, Rhea and Iapetus do show significant internal activities related to convective overturn for extended periods of time. The interiors of Mimas and Tethys do not convect or do so only for brief periods of time early in their thermal histories. Although we use lower densities than previous models, our calculations suggest higher interior temperatures but also thicker rigid shells above the convecting regions. Temperatures in the stagnant lid will allow melting of ammonia-dihydrate. Dione, Rhea and Iapetus may differentiate early and form early oceans, Iapetus only if ammonia is present. Mimas and Tethys with ammonia may differentiate if they accreted in an optically thick nebula with ambient temperatures around 250 K. Our models suggest that the outer shells of the satellites are largely primordial in composition even if the satellites differentiated. In these cases the deep interior may be layered with a pure ice shell underlain by an ammonia dihydrate layer and a rock core.     

Observational constraints on the identification and distribution of chaotic terrain on icy satellites

We outline the observational constraints required to identify chaos regions on Europa. Large incidence angle, rather than high resolution, appears to be the primary observational requirement for identifying chaos. At incidence angles >70°, chaos can be identified on Europa at image resolutions as low as 1.5 km/pixel. Similar images obtained at moderate or low incidence angles (<50°) require image resolutions upwards of ～250 m/pixel to identify chaos. If global images of Europa can be acquired at high incidence angles, the majority of its chaotic terrain can be identified, helping to constrain models of chaos formation and distribution. Furthermore, our results indicate that the areal coverage of chaos may be more uncertain than previously reported, representing as little as 10% or as much as 50% of the non-polar regions of Europa. These guidelines will aid in the development of optical instruments for future Europa missions, as well as other icy bodies, such as Triton.     

Time-dependent flexure of the lithospheres on the icy satellites of Jupiter and Saturn

Previous analyses into flexural deformation on the icy satellites of Jupiter and Saturn have assumed static, elastic lithospheres. Viscous creep within the lithosphere, however, can cause evolution over time. Here, we apply a finite-element model that employs a time-dependent elastic–viscous-plastic rheology in order to investigate flexure on icy satellites. Factors that affect this time-dependent response are those that control creep rates; surface temperature, heat flow, and grain size. Our results show that surface temperature is by far the dominant factor. At higher surface temperatures (100–130 K), the evolution of the deformation is such that the thickness of a modeled elastic lithosphere could vary by up to an order of magnitude, depending on the time scale over which the deformation occurred. Because the flexure observed on icy satellites generally indicates transient high heat flow events, our results indicate that the duration of the heat pulse is an important factor. For the icy worlds of Jupiter and Saturn, static models of lithospheric flexure should be used with caution.     

Application of a radiative transfer model to bright icy satellites

A radiative transfer model, derived largely from the work of B.W. Hapke (1981, J. Geophys. Res.86, 3039–3054) and J.D. Goguen (1981, Ph.D. thesis, Cornell University, Ithaca, N.Y.), is fit to Voyager imaging observations of Europa, Mimas, Enceladus, and Rhea. It is possible to place constraints on the single-scattering albedo, the porosity of the optically active upper regolith, the single-particle phase functions, and, in the cases of Europa and Mimas, the mean slope angle of macroscopic surface features. The texture of the surfaces of the Saturnian satellites appears to be similar to the Earth's moon. However, Europa is found to have a distinctly more compact regolith and a more forward-scattering single-particle phase function.     

Mapping of the icy Saturnian satellites: First results from Cassini-ISS

Images of the icy Saturnian satellites Mimas, Enceladus, Tethys, Dione, Rhea, Iapetus, and Phoebe, derived by the Voyager and Cassini cameras are used to produce new local high-resolution image mosaics as well as global mosaics [http://ciclops.org, http://photojournal.jpl.nasa.gov]. These global mosaics are valuable both for scientific interpretation and for the planning of future flybys later in the ongoing Cassini orbital tour. Furthermore, these global mosaics can be extended to standard cartographic products.     

The icy Galilean satellites: ultraviolet phase curve analysis

Ultraviolet disk-integrated solar phase curves of the icy galilean satellites Europa, Ganymede, and Callisto are presented, using combined data sets from the International Ultraviolet Explorer (IUE), Hubble Space Telescope (HST), and the Galileo Ultraviolet Spectrometer. Global, disk-integrated solar phase curves for all three satellites, in addition to disk-integrated solar phase curves for Europa's leading, trailing, jovian, and anti-jovian hemispheres, are modeled using Hapke's equations for 7 broadband UV wavelengths between 260 and 320 nm. The sparse coverage in solar phase angle, particularly for Ganymede and Callisto, and the noise in the data sets poorly constrain some of the photometric parameter values in the model. However, the results are sufficient for forming a preliminary relationship between the effects of particle bombardment on icy surfaces and photometric scattering properties at ultraviolet wavelengths. Callisto exhibits a large UV opposition surge and a surface comprised of relatively low-backward scattering particles. Europa's surface displays a dichotomy between the jovian and anti-jovian hemispheres (the anti-jovian hemisphere is more backward scattering), while a less pronounced hemispherical variation was detected between the leading and trailing hemispheres. Europa's surface, with the exception of the trailing hemisphere region, appears to have become less backscattering between the late-1970s-early-1980s and the mid-1990s. These results are commensurate with the bombardment history of these surfaces by magnetospheric charged particles.     

Model of explosive eruptions of water vapor and dust on icy satellites

The possibility of generating water vapor and other gaseous products during nonvolcanic explosive eruptions in lithospheres of icy satellites is discussed. Explosive eruptions of ice, with its fragmentation into micro-and nanofragments, can occur in the extensive deep layers of such icy satellites as Europa, Ganymede, Enceladus, etc., if giant cracks are episodically formed in the lithospheres of these satellites. Such cracks can be produced by tidal forces, synchronous resonances of satellites, or especially powerful impacts. The model is based on the recently obtained experimental evidence that explosive ice instability (Bridgman effect) is formed at a strong nonuniform compression in the regions of high pressures and low temperatures. Water films, the thicknesses of which reach several microns, can be formed during the process of the mutual friction of ice fragments during their quasi-liquid flow at the instant of an explosive eruption. About 1–10 dm³ of a water film can be produced in 1 m³ of erupted ice fragments. Water vapor can be formed from a water film when this water boils up after a rapid pressure drop as a result of an ice-water mixture eruption from cracks. A certain amount of gaseous products in the form of hydrogen, oxygen, and ammonia molecules and radicals on their basis can be generated during the sputtering induced by electrons and ions and the dissociation of nanofragments of ice during the process of ice explosive fragmentation as a result of fracto-, tribo-, and secondary emission. The estimates indicate that the volume of water vapor erupted on satellites can be larger than that of discharged ionized gases by a factor of not less than 10⁵–10⁷. Water vapor and microscopic ice fragments can be erupted from cracks in the lithospheres of small Enceladus-type satellites at velocities higher than the second cosmic velocity. Gaseous products generated in such episodic processes can, most probably, substantially contribute to the density of the atmosphere that exists on small icy satellites, but can only insignificantly contribute to this density on large satellites. The stick-slip motions of the most condensed plumes of water vapor and dust, normal to the satellite surface, along the mouths of gigantic cracks may indicate that the proposed model is realistic. Such wanderings of water vapor plumes can result in the synchronous motions of thermal patches on the satellite surface along crack mouths at velocities of about 10 km/h.     

Models,figures, and gravitational moments of the Galilean satellites of Jupiter and icy satellites of Saturn

A general theory for the figures of satellites, which are synchronously rotating in the gravitational field of a planet, is developed to the first approximation. Love numbers, figure parameters, and gravitational moments for two- and three-layer models of the Galilean satellites, Titan, and Saturn's icy satellites are calculated. With the assumed accuracy for flyby measurements of gravitational moments it should be possible to determine the degree of differentiation of Ganymede. The differences between equatorial a and polar c semiaxes, as derived from the observational data, appear to be exaggerated for Io and Mimas (although better agreement between calculated and observed values of (a?c) could be obtained if this satellite had a larger mass). For Enceladus the observed value of (a?c) is in satisfactory agreement with calculations, based on different types of trial models. However, in order to discriminate between different Enceladus trial models, it is necessary to determine the figure parameters more precisely.