Orbital history of the Martian satellites with inferences on their origin

Recent Viking results indicate the Martian satellites are composed of carbonaceous chondritic material, suggesting that Phobos and Deimos were once asteroids captured by Mars. On the other hand, the low eccentricities and inclinations of their orbits on the equator of Mars argue against that hypothesis. This paper presents detailed calculations of the tidal evolution of Phobos and Deimos, considering dissipation in both Mars and its satellites simultaneously and using a new method applicable for any value of the eccentricity. In particular, including precession of the satellites' orbits indicates that they have always remained close to their Laplacian plane, so that the orbital planes of Phobos and Deimos switched from near the Martian orbital plane to the Martian equator once the perturbations due to the planetary oblateness dominated the solar perturbations, as they do presently. The results show that Deimos has been little affected by tides, but several billion (10⁹) years ago, Phobos was in a highly eccentric orbit lying near the common plane of the solar system. This outcome is obtained for very reasonable values of dissipation inside Mars and inside Phobos. Implications for the origin of the Martian satellites are discussed.     

The orbital-thermal evolution and global expansion of Ganymede

The tectonically and cryovolcanically resurfaced terrains of Ganymede attest to the satellite's turbulent geologic history. Yet, the ultimate cause of its geologic violence remains unknown. One plausible scenario suggests that the Galilean satellites passed through one or more Laplace-like resonances before evolving into the current Laplace resonance. Passage through such a resonance can excite Ganymede's eccentricity, leading to tidal dissipation within the ice shell. To evaluate the effects of resonance passage on Ganymede's thermal history we model the coupled orbital-thermal evolution of Ganymede both with and without passage through a Laplace-like resonance. In the absence of tidal dissipation, radiogenic heating alone is capable of creating large internal oceans within Ganymede if the ice grain size is 1 mm or greater. For larger grain sizes, oceans will exist into the present epoch. The inclusion of tidal dissipation significantly alters Ganymede's thermal history, and for some parameters (e.g. ice grain size, tidal Q of Jupiter) a thin ice shell (5 to 20 km) can be maintained throughout the period of resonance passage. The pulse of tidal heating that accompanies Laplace-like resonance capture can cause up to 2.5% volumetric expansion of the satellite and contemporaneous formation of near surface partial melt. The presence of a thin ice shell and high satellite orbital eccentricity would generate moderate diurnal tidal stresses in Ganymede's ice shell. Larger stresses result if the ice shell rotates non-synchronously. The combined effects of satellite expansion, its associated tensile stress, rapid formation of near surface partial melt, and tidal stress due to an eccentric orbit may be responsible for creating Ganymede's unique surface features.     

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

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

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Did tidal deformation power the core dynamo of Mars?

We first show that 7 out of the 20 giant impact basins of Mars recently reported by Frey [Frey, H., 2008. Geophys. Res. Lett. 35. L13203] trace a great circle on Mars. The other five basins trace another great circle and still the other three basins trace yet another great circle. The latter great circle is in good agreement with the pre-Tharsis equator of Mars that is estimated from modeling crustal magnetic anomalies [Arkani-Hamed, J., 2001. Geophys. Res. Lett. 28, 3409-3412] and diagonalizing the moment of inertia of Mars after removing the loading effects of Tharsis bulge [Sprenke, K.F., Baker, L.L., Williams, A.F., 2005. Icarus 174, 486-489]. It is shown in this paper that the three great circles were likely the equatorial plane of Mars at certain times and Mars experienced appreciable polar wander. The great circles also indicate that the asteroids that created the basins were satellites of Mars whose orbits decayed in time through spin-orbit coupling with tidally deforming Mars, and eventually impacted on the planet creating the giant basins at around 4 Ga. The orbital dynamics of four largest asteroids show that they could have orbited Mars for several hundred million years if they were retrograde satellites. Continual elliptical straining of otherwise circular fluid streamlines of the liquid core of Mars by tidal deformation could have exerted a strong strain that was large enough to overcome dissipation and excite the elliptical instability inside the core. We investigate the physical properties of the martian core that are required to allow the tidal deformation to power the core dynamo, i.e., the growth time of the elliptical instability to become shorter than the dissipation time. The tidal energy dissipation rate inside Mars caused by even only one of the 4 largest asteroids is found to be over two orders of magnitude greater than the magnetic energy dissipation rate in the core, indicating that if only one of the 4 largest asteroids were orbiting in retrograde sense, it would have likely powered the core dynamo of Mars for several hundred million years.     

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.     

伽利略卫星的平衡形状参数和潮汐耗散因子

木星的四颗大卫星都是同步轨旋卫星，常被称为伽利略卫星，美国发射的伽利略飞船自1995年12月抵达木星系统后，的几年来对木星的四颗伽利略卫星进行了一系列的探测，利用飞船探测的最新资料作为约束条件，建立了伽利略卫星的一组内部结构模型，然后按照同步轨旋卫星的形态理论公式计算了它们的平衡形太参数及潮汐耗散因子等。     

Thermal Equilibrium States of Europa's Ice Shell: Implications for Internal Ocean Thickness and Surface Heat Flow

Ice-shell thickness and ocean depth are calculated for steady state models of tidal dissipation in Europa's ice shell using the present-day values of the orbital elements. The tidal dissipation rate is obtained using a viscoelastic Maxwell rheology for the ice, the viscosity of which has been varied over a wide range, and is found to strongly increase if an (inviscid) internal ocean is present. To determine steady state values, the tidal dissipation rate is equated to the heat-transfer rate through the ice shell calculated from a parameterized model of convective heat transfer or from a thermal conduction model, if the ice layer is found to be stable against convection. Although high dissipation rates and heat fluxes of up to 300 mWm⁻² are, in principle, possible for Europa, these values are unrealistic because the states for which they are obtained are thermodynamically unstable. Equilibrium models have surface heat flows around 20 mWm⁻² and ice-layer thicknesses around 30 km, which is significantly less than the total thickness of the H₂O-layer. These results support models of Europa with ice shells a few tens of kilometers thick and around 100-km-thick subsurface oceans.     

Ice-rich terrain in Gusev Crater, Mars?

The morphology of materials on the floor of Gusev Crater (14° S, 175° W), Mars, imply a history of volcanism and subsequent removal of an ice-rich deposit. Fluid lava flows observed in the western portion of Gusev Crater paradoxically terminate in a steep, thick (<60 m) flow front adjacent to hummocky terrain. The hummocky terrain is morphologically similar to deglaciated terrain on Earth, generated when glacial debris are left behind after the glacier has retreated. We propose the following scenario for the floor of Gusev Crater. First, ice-rich material was deposited adjacent to Thira Crater. Second, fluid lavas were emplaced and ponded against the ice-rich deposits. At some later time, the ice within the deposit sublimated, leaving hummocky terrain. Current age estimates for the Gusev flows are Hesperian, suggesting that the ice removal occurred in the upper Hesperian or more recently. If this hypothesis is correct, quench features (glassy rinds, columnar jointing) should be observed at the lava flow margin; the hummocky deposit should be poorly sorted, angular debris.     

Accurate Mars Express orbits to improve the determination of the mass and ephemeris of the Martian moons

The determination of the ephemeris of the Martian moons has benefited from observations of their plane-of-sky positions derived from images taken by cameras onboard spacecraft orbiting Mars. Images obtained by the Super Resolution Camera (SRC) onboard Mars Express (MEX) have been used to derive moon positions relative to Mars on the basis of a fit of a complete dynamical model of their motion around Mars. Since, these positions are computed from the relative position of the spacecraft when the images are taken, those positions need to be known as accurately as possible. An accurate MEX orbit is obtained by fitting two years of tracking data of the Mars Express Radio Science (MaRS) experiment onboard MEX. The average accuracy of the orbits has been estimated to be around 20–25 m. From these orbits, we have re-derived the positions of Phobos and Deimos at the epoch of the SRC observations and compared them with the positions derived by using the MEX orbits provided by the ESOC navigation team. After fit of the orbital model of Phobos and Deimos, the gain in precision in the Phobos position is roughly 30 m, corresponding to the estimated gain of accuracy of the MEX orbits. A new solution of the GM of the Martian moons has also been obtained from the accurate MEX orbits, which is consistent with previous solutions and, for Phobos, is more precise than the solution from the Mars Global Surveyor (MGS) and Mars Odyssey (ODY) tracking data. It will be further improved with data from MEX-Phobos closer encounters (at a distance less than 300 km). This study also demonstrates the advantage of combining observations of the moon positions from a spacecraft and from the Earth to assess the real accuracy of the spacecraft orbit. In turn, the natural satellite ephemerides can be improved and participate to a better knowledge of the origin and evolution of the Martian moons.     

Morphologic and topographic analyses of debris aprons in the eastern Hellas region, Mars

This paper presents new, detailed analyses of small-scale morphologic and topographic characteristics of martian debris aprons that support Viking-based interpretations of debris aprons as ice-rich flow features derived from local uplands. Fifty-four debris apron complexes in the eastern Hellas region of Mars were examined using Mars Global Surveyor data sets, including Mars Orbiter Camera images and Mars Orbiter Laser Altimeter topographic profiles. Consistent patterns in a suite of small-scale surface textures and geomorphic features observed throughout the population reflect a history of viscous flow and surface degradation through wind ablation and loss of contained ice. A wide variety of shapes seen in topographic profile reveal variations in distribution of contained ice and different stages of apron development and degradation. The results of this study provide new evidence consistent with multiple models of apron formation, including rock glacier, debris-covered glacier, and ice-rich landslide models. Typical eastern Hellas debris aprons formed from a series of large-scale events, emplacing debris that was enriched initially or later by ground ice, complemented by small-scale mass wasting of multiple styles and postemplacement flow of apron masses.     

A morphological view on potential niches for exobiology on Mars

The discovery of microbiota in the Dry Valleys of Antarctica has encouraged the construction of new models of Martian ecosystems in order to determine if life could have once existed on Mars. The Antarctic cyanobacteria reside just below the surface of sandstone rocks where they are protected from the extreme cold and dry environment. Analogy with the Antarctic Dry Valleys supports speculation that hypothetical micro-organisms existed on Mars in the early history of the planet and could have migrated into suitable rocks as the availability of liquid water decreased. Although evidence for sandstone layers on Mars has not been substantiated, the palaeohydrology of Martian fluvial valleys (MFVs) reveals the evidence of lake bed sediment depositions which have formed consolidated sediments. As the MFVs formation may result from underground drainage processes, the sediment material would be expected to contain debris such as pumice washload, and pumilith of volcanic and meteoritic origin. These materials may have formed consolidated porous terrains similar to the Antarctic sandstone. Therefore, the endolithic model is consistent with the Martian liquid water habitat model of perenially ice-covered lakes.     

Gas-drag-assisted capture of Himalia's family

To elucidate the capture of Jupiter's outer moons, we reverse-evolve satellites from their present orbits to their original heliocentric paths in the presence of Jupiter's primordial circumplanetary disk (Lubow et al., 1999, Astrophys. J. 526, 1001-1012; Canup and Ward, 2003, Astron. J. 124, 3404-3423). Our orbital histories use a symplectic integrator that allows dissipation. We assume that the present satellites Himalia, Elara, Lysithea, Leda, and S/2000 J11 are collisional fragments of a single parent. Our simulations show that this “prograde-cluster progenitor” (PCP) could be derived from objects with heliocentric orbits like those of the Hilda asteroid group. We show analytically that this capture is energetically possible. We also compare the spectroscopic characteristics of the prograde cluster members (Grav et al., 2003, Icarus, submitted for publication) with those of the Hildas, and conclude that the surface color of the prograde-cluster progenitor is consistent with an origin within the Hilda group. Accordingly, gas drag in the primordial jovian nebula is found to offer a plausible explanation for the origin of the prograde cluster. A similar capture mechanism is proposed for Saturn's Phoebe.     

Constraints to Uranus' Great Collision III: The Origin of the Outer Satellites

The theory of satellite loss resulting from a giant impact on Uranus (Parisi and Brunini 1997, Planet. Space Sci.45, 181-187) is revisited, in the light of the discovery of its five outer moons (Gladman et al. 1998, Nature392, 897-899; Gladman et al. 2000, Icarus147, 320-324; erratum 148, 320). Physical conditions and dynamical constraints in the great collision scenario and restrictions in the possible mechanisms for the origin of the outer uranian satellites are obtained from the knowledge of their actual orbital properties. We conclude that the existence of these moons implies that their origin must be connected to a breakup process. Other scenarios for their origin cast serious doubts on the occurrence of a giant collision at the end of Uranus' formation process to account for its large spin axis obliquity.     

Are Phobos and Deimos the result of a giant impact?

Despite many efforts an adequate theory describing the origin of Phobos and Deimos has not been realized. In recent years a number of separate observations suggest the possibility that the martian satellites may have been the result of giant impact. Similar to the Earth–Moon system, Mars has too much angular momentum. A planetesimal with 0.02 Mars masses must have collided with that planet early in its history in order for Mars to spin at its current rate (Dones, L., Tremaine, S. [1993]. Science 259, 350–354). Although subject to considerable error, current crater-scaling laws and an analysis of the largest known impact basins on the martian surface suggest that this planetesimal could have formed either the proposed 10,600 by 8500-km-diameter Borealis basin, the 4970-km-diameter Elysium basin, the 4500-km-diameter Daedalia basin or, alternatively, some other basin that is no longer identifiable. It is also probable that this object impacted Mars at a velocity great enough to vaporize rock (>7 km/s), which is necessary to place large amounts of material into orbit. If material vaporized from the collision with the Mars-spinning planetesimal were placed into orbit, an accretion disk would have resulted. It is possible that as material condensed and dissipated beyond the Roche limit forming small, low-mass satellites due to gravity instabilities within the disk. Once the accretion disk dissipated, tidal forces and libration would have pulled these satellites back down toward the martian surface. In this scenario, Phobos and Deimos would have been among the first two satellites to form, and Deimos the only satellite formed—and preserved—beyond synchronous rotation. The low mass of Phobos and Deimos is explained by the possibility that they are composed of loosely aggregated material from the accretion disk, which also implies that they do not contain any volatile elements. Their orbital eccentricity and inclination, which are the most difficult parameters to explain easily with the various capture scenarios, are the natural result of accretion from a circum-planetary disk.     

Mars geodesy,rotation and gravity

This review provides explanations of how geodesy, rotation and gravity can be addressed using radioscience data of an orbiter around a planet or of the lander on its surface.The planet Mars is the center of the discussion.The information one can get from orbitography and radioscience in general concerns the global static gravitational field, the time variation of the gravitational field induced by mass exchange between the atmosphere and the ice caps, the time variation of the gravitational field induced by the tides, the secular changes in the spacecraft's orbit induced by the little moons of Mars named Phobos and Deimos, the gravity induced by particular targets, the Martian ephemerides, and Mars' rotation and orientation.The paper addresses as well the determination of the geophysical parameters of Mars and, in particular, the state of Mars' core and its size, which is important for understanding the planet's evolution.Indeed, the state and dimension of the core determined from the moment of inertia and nutation depend in turn on the percentage of light elements in the core as well as on the core temperature, which is related to heat transport in the mantle.For example, the radius of the core has implications for possible mantle convection scenarios and, in particular, for the presence of a perovskite phase transition at the bottom of the mantle.This is also important for our understanding of the large volcanic province Tharsis on the surface of Mars.     

Coupled convection and tidal dissipation in Europa’s ice shell

We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 10¹³ Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other features.     

Tidal dissipation within large icy satellites: Applications to Europa and Titan

This paper describes a new approach based on variational principles to calculate the radial distribution of tidal energy dissipation in any satellite. The advantage of the model with respect to classical solutions, is that it relates in a straightforward way the radial distribution of the time-averaged dissipation rate to its sensitivity to the corresponding distribution of viscoelastic parameters. This method is applied to Io-, Europa-, and Titan-like interiors, and it is tested against the results obtained by two classical methods by determining global dissipation as well as radial and lateral distributions within satellite interiors. By exploring systematically the different parameters defining the interior models, we demonstrate that the presence of a deep ocean below an outer ice layer strongly influences the tidal dissipation distribution in both the outer ice layer and in the innermost part of the satellite. On the one hand, the ocean by imposing a large radial displacement at the base of the outer ice I layer, controls the distribution of tidal strain rate within the outer layer, making the tidal strain rate field very weakly sensitive to the viscosity variations. Conversely, in the high-pressure ice layer below the ocean, both tidal strain rate and dissipation are very sensitive to any variation of the ice viscosity. On the other hand, for identical structures of the mantle and of the core, the presence of a subsurface ocean reduces the strength of dissipation in the silicate mantle. The existence of a liquid layer within Europa makes models of the silicate mantle less dissipative than the predictions for Io.     

Collisional breakup of particles in a planetary ring

Jeffreys (1947) estimated the size of fragments resulting from breakup of a satellite inside the Roche limit, obtaining a result of ～100 km. This result does not allow for the further breakup of the fragments due to collisions among themselves, which should reduce the maximum size to ?3 km for rock, or ?1 km for ice. This result affects not only Jeffrey's speculations as to the origin of Saturn's rings, but also recent speculations on the origin of the moon by capture and the possible tidal destruction of satellites of Mercury or Venus.     

Thermal-orbital evolution of Io and Europa

Coupled thermal-orbital evolution models of Europa and Io are presented. It is assumed that Io, Europa, and Ganymede evolve in the Laplace resonance and that tidal dissipation of orbital energy is an internal heat source for both Io and Europa. While dissipation in Io occurs in the mantle as in the mantle dissipation model of Segatz et al. (1988, Icarus 75, 187), two models for Europa are considered. In the first model dissipation occurs in the silicate mantle while in the second model dissipation occurs in the ice shell. In the latter model, ice shell melting and variations of the shell thickness above an ocean are explicitly included. The rheology of both the ice and the rock is cast in terms of a viscoelastic Maxwell rheology with viscosity and shear modulus depending on the average temperature of the dissipating layer. Heat transfer by convection is calculated using a parameterization for strongly temperature-dependent viscosity convection. Both models are consistent with the present orbital elements of Io, Europa, and Ganymede. It is shown that there may be phases of quasi-steady evolution with large or small dissipation rates (in comparison with radiogenic heating), phases with runaway heating or cooling and oscillatory phases during which the eccentricity and the tidal heating rate will oscillate. Europa's ice thickness varies between roughly 3 and 70 km (dissipation in the silicate layer) or 10 and 60 km (dissipation in the ice layer), suggesting that Europa's ocean existed for geological timescales. The variation in ice thickness, including both convective and purely conductive phases, may be reflected in the formation of different geological surface features on Europa. Both models suggest that at present Europa's ice thickness is several tens of km thick and is increasing, while the eccentricity decreases, implying that the satellites evolve out of resonance. Including lithospheric growth in the models makes it impossible to match the high heat flux constraint for Io. Other heat transfer processes than conduction through the lithosphere must be important for the present Io.