Crack azimuths on Europa: time sequence in the southern leading face

The formation sequence of prominent ridges and other tectonic lineaments on the southern portion of the leading hemisphere of Europa is determined from cross-cutting relationships. These selected features formed fairly recently relative to most of the surface; older lineaments no longer retain clear evidence of cross-cutting. If we assume that this sequence represents the order of formation of cracks that underlie each lineament, and that the orientation of each crack was determined by tidal stress whose azimuth varies monotonically counter-clockwise with time, the azimuth must have rotated more than 740°, which would correspond to the change in tidal stress over two periods of nonsynchronous rotation (relative to the direction of Jupiter). However, that interpretation is not necessarily compelling, because the observed orientations of cross-cutting lineaments are not densely spaced over these cycles; in fact, the sequence would fit nearly as well into an arbitrary model with rotation in the opposite sense from that predicted by theory. This tectonic record may have formed over many more rotational cycles, such that typically only a few cracks form per cycle, which would be consistent with evidence from considerations of cycloidal crack patterns. Sets of cracks that cluster near certain azimuth orientations appear to be parts of globe-encircling lineament systems and may result from other effects, perhaps polar wander that occurred rapidly relative to nonsynchronous rotation.     

Modeling stresses on satellites due to nonsynchronous rotation and orbital eccentricity using gravitational potential theory

The tidal stress at the surface of a satellite is derived from the gravitational potential of the satellite's parent planet, assuming that the satellite is fully differentiated into a silicate core, a global subsurface ocean, and a decoupled, viscoelastic lithospheric shell. We consider two types of time variability for the tidal force acting on the shell: one caused by the satellite's eccentric orbit within the planet's gravitational field (diurnal tides), and one due to nonsynchronous rotation (NSR) of the shell relative to the satellite's core, which is presumed to be tidally locked. In calculating surface stresses, this method allows the Love numbers h and ?, describing the satellite's tidal response, to be specified independently; it allows the use of frequency-dependent viscoelastic rheologies (e.g. a Maxwell solid); and its mathematical form is amenable to the inclusion of stresses due to individual tides. The lithosphere can respond to NSR forcing either viscously or elastically depending on the value of the parameter , where μ and η are the shear modulus and viscosity of the shell respectively, and ω is the NSR forcing frequency. Δ is proportional to the ratio of the forcing period to the viscous relaxation time. When Δ?1 the response is nearly fluid; when Δ?1 it is nearly elastic. In the elastic case, tensile stresses due to NSR on Europa can be as large as ∼3.3 MPa, which dominate the ∼50 kPa stresses predicted to result from Europa's diurnal tides. The faster the viscous relaxation the smaller the NSR stresses, such that diurnal stresses dominate when Δ?100. Given the uncertainty in current estimates of the NSR period and of the viscosity of Europa's ice shell, it is unclear which tide should be dominant. For Europa, tidal stresses are relatively insensitive both to the rheological structure beneath the ice layer and to the thickness of the icy shell. The phase shift between the tidal potential and the resulting stresses increases with Δ. This shift can displace the NSR stresses longitudinally by as much as 45° in the direction opposite of the satellite's rotation.     

Elastic ice shells of synchronous moons: Implications for cracks on Europa and non-synchronous rotation of Titan

A number of synchronous moons are thought to harbor water oceans beneath their outer ice shells. A subsurface ocean frictionally decouples the shell from the interior. This has led to proposals that a weak tidal or atmospheric torque might cause the shell to rotate differentially with respect to the synchronously rotating interior. Applications along these lines have been made to Europa and Titan. However, the shell is coupled to the ocean by an elastic torque. As a result of centrifugal and tidal forces, the ocean would assume an ellipsoidal shape with its long axis aligned toward the parent planet. Any displacement of the shell away from its equilibrium position would induce strains thereby increasing its elastic energy and giving rise to an elastic restoring torque. In the investigation reported on here, the elastic torque is compared with the tidal torque acting on Europa and the atmospheric torque acting on Titan.Regarding Europa, it is shown that the tidal torque is far too weak to produce stresses that could fracture the ice shell, thus refuting an idea that has been widely advocated. Instead, it is suggested that the cracks arise from time-dependent stresses due to non-hydrostatic gravity anomalies from tidally driven, episodic convection in the satellite’s interior.Two years of Cassini RADAR observations of Titan’s surface have been interpreted as implying an angular displacement of ∼0.24° relative to synchronous rotation. Compatibility of the amplitude and phase of the observed non-synchronous rotation with estimates of the atmospheric torque requires that Titan’s shell be decoupled from its interior. We find that the elastic torque balances the seasonal atmospheric torque at an angular displacement ?0.05°, effectively coupling the shell to the interior. Moreover, if Titan’s surface were spinning faster than synchronous, the tidal torque tending to restore synchronous rotation would almost certainly be larger than the atmospheric torque. There must either be a problem with the interpretation of the radar observations, or with our basic understanding of Titan’s atmosphere and/or interior.     

Patterns of fracture and tidal stresses due to nonsynchronous rotation: Implications for fracturing on Europa

This study considers the global patterns of fracture that would result from nonsynchronous rotation of a tidally distorted planetary body. The incremental horizontal stresses in a thin elastic or viscous shell due to a small displacement of the axis of maximum tidal elongation are derived, and the resulting stress distributions are applied to interpret the observed pattern of fracture lineaments on Europa. The observed pattern of lineaments can be explained by nonsynchronous rotation if these features formed by tension fracturing and dike emplacement. Tension fracturing can occur for a small displacement of the tidal axis, so that the resulting lineaments may be consistent with other evidence suggesting a young age for the surface.     

Is There Evidence for Polar Wander on Europa?

The orientations of the albedo lineaments, bands, and lineations on Europa's surface have been compared in previous studies with the global stress fields set up by orbital eccentricity, orbital recession, and nonsynchronous rotation. Of these orbital and rotational effects, nonsynchronous rotation, combined with an offsetting of the tidal bulge, comes closest to providing agreement between the stress field generated and the lineation orientations, if the lineations trace tension or extension fractures (McEwen 1986.Nature321, 49–51). However, inferred minimum principal stress directions for a broad region of wedge-shaped bands near the anti-Jove point cannot satisfactorily be accounted for by any of the stress fields above, but are consistent with the stresses resulting from a rotation of Europa's ice shell about an axis through the sub- and anti-Jove points, clockwise as seen from the anti-Jove hemisphere (P. M. Schenk and W. B. McKinnon 1989.Icarus79, 75–100). Calculations by Ojakangas and Stevenson (1989.Icarus81, 220–241) of the thermal state of Europa's ice shell indicate that spatial variations in the thickness of the shell may cause it to undergo such a reorientation. We have investigated whether any reorientation of the shell about an axis through the sub- and anti-Jove points produces a stress field consistent with the full, global set of prominent lineations on Europa's surface. We find that no such reorientation provides a good fit between the lineations and plausible fracture orientations derived from the principal stress trajectories. Topographic ridges, identified in a limited zone south of the anti-Jove point, are roughly consistent with compression due to clockwise polar wander, but the orientations of these ridges may be strongly biased by illumination direction. Within the limitations of the presently available imagery, nonsynchronous rotation is still the most likely cause of the prominant fractures on Europa's surface, and the best specific, albeit regionally limited, tectonic evidence consistent with recent polar wander remains the wedge-shaped bands.     

Cycloid crack sequences on Europa: Relationship to stress history and constraints on growth mechanics based on cusp angles

The original model developed to explain cycloidal cracks on Europa interprets cycloids as tensile fractures that grow in a curved path in response to the constantly rotating diurnal tidal stress field. Cusps form when a new cycloid crack segment propagates at an angle to the first in response to a rotation of the principal tidal stress orientation during a period of no crack growth. A recent revised model states that a cycloid cusp forms through the creation of a secondary fracture called a tailcrack at the tip of an existing cycloid segment during shearing motion induced by the rotating tidal stress field. As the tailcrack propagates away from the cusp, it becomes the next cycloid segment in the chain. The qualitative tailcrack model uniquely accounts for the normal and shear stresses that mechanically must resolve onto the tip of an existing cycloid segment at the instant of cusp formation. In this work, we provide a quantitative framework and test of the hitherto purely conceptual tailcrack model. We first present a relative age sequence inferred from geologic mapping of multiply cross-cutting cycloids in Europa's trailing hemisphere and place this into the context of the global stress history. The age sequence requires a cumulative minimum of 630° of shell reorientation due to nonsynchronous rotation to account for the observed range of orientations of cycloids of different ages. We determined the back-rotated longitudes of formation of two cycloid chain examples and used mathematical modeling of europan tidal stresses to show that the tailcrack model for cusp formation is not only viable, but places constraints on the overall development of a cycloid chain by controlling the timing of cusp development within Europa's orbit. For all cusps analyzed, the exact ratio of resolved shear to normal stress required to form the cusp angles by a process of tailcracking, as governed by the principles of linear elastic fracture mechanics, is produced at the tip of a shearing cycloid segment during Europa's orbit. Cusp formation occurs after the point in the orbit at which the maximum tensile principal tidal stress occurs, implying that tensile tidal stresses are not directly responsible for cusp development. Instead, cusps develop when a tailcrack forms at the tip of a cycloid segment in response to the highly perturbed stress field induced during concomitant opening and shearing at the tip of the cycloid segment.     

How fast do Galilean satellites spin?

Each of the Galilean satellites, as well as most other satellites whose initial rotations have been substantially altered by tidal dissipation, has been widely assumed to rotate synchronously with its orbital mean motion. Such rotation would require a small permanent asymmetry in the mass distribution in order to overcome the small mean tidal torque. Since Io and Europa may be substantially fluid, they may not have the strenght to support the required permanent asymmetry. Thus, each may rotate at the unknown but slightly nonsynchronous rate that corresponds to zero mean tidal torque. This behaviour may be observable by Galileo spacecraft imaging. It may help explain the longitudinal variation of volcanism on Io and the cracking of Europa's crust.     

Habitats and taphonomy of Europa

Jupiter's moon Europa possesses an icy shell kilometers thick that may overlie a briny ocean. The inferred presence of water, tidal and volcanic energy, and nutrients suggests that Europa is potentially inhabited by some kind of life; indeed Europa is a primary target in the search for life in the Solar System although no evidence yet exists for any kind of life. The thickness of the icy crust would impose limits on life, but at least 15 broad kinds of habitats seem possible for Europa. They include several on the sea floor, at least 3 in the water column, and many in the ice itself. All of these habitats are in, or could be transported to, the icy shell where they could be exposed by geologic activity or impacts so they might be explored from the surface or orbit by future planetary missions. Taphonomic processes that transport, preserve, and expose habitats include buoyant ice removing bottom habitats and sediment to the underside of the ice, water currents depositing components of water column habitats on the ice bottom, cryovolcanoes depositing water on the surface, tidal pumping bringing water column and ice habitats to the near-surface ice, and subice freezing and diapiric action incorporating water column and bottom ice habitats into the lower parts of the icy shell. The preserved habitats could be exposed at or near the surface of Europa chiefly in newly-formed ice, tilted or rotated ice blocks, ridge debris, surface deposits, fault scarps, the sides of domes and pits, and impact craters and ejecta. Future exploration of Europa for life must consider careful targeting of sites where habitats are most likely preserved or exist close to the surface.     

Tidal Stress Patterns on Europa's Crust

The ice crust of Europa probably floats over a deep liquid-water ocean, and has been continually resurfaced by tectonic and thermal processes driven by tides. Tidal working causes rotational torque, surface stress, internal heating, and orbital evolution. The stress patterns expected on such a crust due to reorientation of the tidal bulge by non-synchronous rotation and due to orbital eccentricity, which introduces periodic ('diurnal') variations in the tide, are shown as global maps. By taking into account the finite rate of crack propagation, global maps are generated of cycloidal features and other distinctive patterns, including the crack shapes characteristic of the wedges region and its antipode on the sub-Jovian hemisphere. Theoretical maps of tidal stress and cracking can be compared with observed tectonics, with the possibility of reconstructing the rotational history of the satellite.     

Patterns of fracture and tidal stresses on Europa

the hypothesis that lineaments on Europa are fractures produced by tidal distortion and planetary volume change is examined by comparing the orientations of dark bands, triple bands, and cuspate ridges to fracture patterns predicted for tidal distortion due to orbital recession and orbital eccentricity. If short, reticulate dark band nnear the anti-Jove point are tension cracks which formed in response to tidal distortion, they could only have been produced by orbital eccentricity. Long, arcuate dark band and triple bands peripheral to the anti-Jove point orientations which suggest that they are strike-slip faults which formed in response to orbital recession. If cuspate ridges are compressional features, their orientations and distribution suggest that they formed in response to combined orbital recession and a decrease in planetary volume. Stresses due to orbital eccentricit could have produced tension cracks near the anti-Jove point only if tensile failure occurred either prior to the accumulation of orbital recession stresses or after they had relaxed. Surface fracturing, if a consequence of tidal deformation, places important constraints on the orbital evolution of Europa.     

Do explosive ice ejections occur on Jupiter’s and Saturn’s satellites?

The possibility of an explosive mechanical instability of ice (the Bridgman effect) in the thick icy shells of Jupiter’s and Saturn’s satellites is discussed in principle. The Bridgman effect is an explosive instability of dielectric solid bodies, which disintegrate into microscopic fragments under a quasistatic uniaxial loading in open compression systems at high pressures. The explosive instabilities of ice recently discovered in laboratory experiments with the Bridgman effect are also expected to occur in the extensive deep layers of the shells of icy planetary satellites (for example, in the case of episodical formation of major cracks in their lithospheres due to tidal forces, nonsynchronous rotation of the satellites, or extremely powerful impacts). The depths of occurrence of mechanically unstable ice in the thick crusts of Ganymede, Europa, and Titan, taken as examples, are crudely estimated using a pure-ice model without a possible ammonia admixture. The estimated thickness of the explosive-instability zone in the icy crust of Ganymede (under the assumption that the crust is ～75 km thick) ranges from ～7 to ～27 km at depths from ～40 to ～67 km, depending on the scaling parameter E = 0.2–1. This parameter relates the experimentally determined thicknesses of the ice samples in which the Bridgman effect occurs under laboratory conditions to the expected thicknesses of the explosively unstable layers in the envelopes of the icy satellites. Explosive effects are possible not throughout the entire thickness of the unstableice layer but only within some part of it, several centimeters to several tens of meters in thickness. According to the estimated location of the unstable layer in the crust of Europa (for an assumed crust thickness of ～30 km), such a layer can exist only at scaling factors E < 0.6 at depths ranging from ～21 to ～28 km. For Titan, if its crust is ～100 km thick, the thickness of the unstable layer is similarly estimated to range from ～15 to ～55 km at depths from ～37 to ～92 km for a scaling parameter E lying within the range 0.2–1. At E 0.2, which is quite possible, explosive instabilities of ice could also be expected on the Earth, in the icy shells of Antarctica and Greenland at depths from ～1 to ～1.5 km.     

Tectonic patterns on reoriented and despun planetary bodies

Several processes may produce global tectonic patterns on the surface of a planetary body. The stresses associated with distortions of biaxial figures due to despinning or reorientation were first calculated by Vening Meinesz [Vening Meinesz, F.A., 1947. Trans. Am. Geophys. Union 28 (1), 1-23]. We adopt a mathematically equivalent, but physically more meaningful treatment for distortions associated with rotation. The new approach allows us to find analytic solutions for the general case of stresses associated with distortions of biaxial or triaxial planetary figures. Distortions of biaxial figures may be driven by variations in rotation rate, rotation axis orientation, or the combination of both. Distortions of triaxial figures may be driven by the same mechanisms and/or variations in tidal axis orientation for tidally deformed satellites. While the magnitude of the resulting stresses depends on the adopted elastic and physical parameters, the expected tectonic pattern is independent of these parameters for these mechanisms. Reorientation of the rotation/tidal axis alone is expected to produce normal/thrust faulting provinces enclosing the initial rotation/tidal poles, and thrust/normal faulting provinces enclosing the final rotation/tidal poles. Reorientation of both the rotation and tidal axis results in a wide variety of tectonic patterns for different reorientation geometries. On Europa, the tidal axis reorientation which generally accompanies rotation axis reorientations may provide an alternative explanation for tectonic features that have been interpreted as evidence for non-synchronous rotation. The observed tectonic pattern on Enceladus is more easily explained by a large reorientation (∼90°) of the rotation axis, than by rotation rate variations.     

The librations, shape, and icy shell of Europa

Europa, the smallest of the Galilean satellites, has a young icy surface and most likely contains an internal ocean. The primary objective of possible future missions to Europa is the unambiguous detection and characterization of a subsurface ocean. The thickness of the overlying icy shell provides important information on the thermal evolution of the satellite and on the interaction between the ocean and the surface, the latter being fundamental for astrobiology. However, the thickness is not well known, and estimates range from several hundred of meters to some ten of kilometers. Here, we investigate the use of libration (rotation variation) observations to study the interior structure of Europa and in particular its icy shell. A dynamical libration model is developed, which includes gravitational coupling between the icy shell and the heavy solid interior. The amplitude of the main libration signal at 3.55 days (the orbital period) is shown to depend on Europa's shape and structure. Models of the interior structure of Europa are constructed and the equatorial flattening of the internal layers, which are key parameters for the libration, are calculated by assuming that Europa is in hydrostatic equilibrium. Europa's flattened shape is assumed to be due to rotation and permanent tides, and we extend the classical Radau equation for rotationally flattened bodies to include also tidal deformation. We show that the presence of an ocean increases the amplitude of libration by about 10%, depending mainly on the thickness of the icy shell. Therefore, libration observations offer possibility of detection of a subsurface ocean in Europa and estimation of the thickness of its overlying icy shell.     

Effects of low-viscous layers and a non-zero obliquity on surface stresses induced by diurnal tides and non-synchronous rotation: The case of Europa

In this study we present a semi-analytical Maxwell-viscoelastic model of the variable tidal stress field acting on Europa’s surface. In our analysis, we take into account surface stresses induced by the small eccentricity of Europa’s orbit, the non-zero obliquity of Europa’s spin axis - both acting on a diurnal 3.55-days timescale - and the reorientation of the ice shell as a result of non-synchronous rotation (NSR). We assume that Europa’s putative ocean is covered by an ice shell, which we subdivide in a low-viscous and warm lower ice layer (asthenosphere, viscosity 10¹²-10¹⁷ Pa s), and a high-viscous and cold upper ice layer (lithosphere, viscosity 10²¹ Pa s).Viscoelastic relaxation influences surface stresses in two ways: (1) through viscoelastic relaxation in the lithosphere and (2) through the viscoelastic tidal response of Europa’s interior. The amount of relaxation in the lithosphere is proportional to the ratio between the period of the forcing mechanism and the Maxwell time of the high-viscous lithosphere. As a result, this effect is only relevant to surface stresses caused by the slow NSR mechanism. On the other hand, the importance of the viscoelastic response on surface stresses is proportional to the ratio between the relaxation time (τ_j) of a given viscoelastic mode j and the period of the forcing function. On a diurnal timescale the fast relaxation of transient modes related to the low viscosity of the asthenosphere can alter the magnitude and phase shift of the diurnal stress field at Europa’s surface. The effects are largest, up to 20% in magnitude and 7° in phase for ice rigidities lower than 3.487 GPa, when the relaxation time of the aforementioned transient modes approaches the inverse of the average angular rate of Europa’s orbit. On timescales relevant for NSR (>10⁴ years) the magnitude and phase shift of NSR surface stresses can be affected by viscoelastic relaxation of the ocean-ice boundary. This effect, however, becomes only important when the behavior of the lithosphere w.r.t. NSR approaches the fluid limit, i.e. for strong relaxation in the lithosphere. The combination of NSR and diurnal stresses for different amounts of viscoelastic relaxation of NSR stresses in the lithosphere leads to a large variety of global stress fields that can explain the formation of the large diversity of lineament morphologies observed on Europa’s surface. Variation of the amount of relaxation in the lithosphere is likely due to changes in the spin rate of Europa and/or the rheological properties of the surface.In addition, we show that a small obliquity(<1°) can have a considerable effect on Europa’s diurnal stress field. A non-zero obliquity breaks the symmetric distribution of stress patterns with respect to the equator, thereby affecting the magnitude and orientation of the principal stresses at the surface. As expected, increasing the value of Europa’s obliquity leads to larger diurnal stresses at the surface, especially when Europa is located 90° away from the nodes formed by the intersection of its orbital and equatorial planes.     

Crack azimuths on Europa: the G1 lineament sequence revisited

The tectonic sequence in the anti-jovian area covered by regional mapping images from Galileo's orbit E15 is determined from a study of cross-cutting relationships among lineament features. The sequence is used to test earlier results from orbit G1, based on lower resolution images, which appeared to display a progressive change in azimuthal orientation over about 90° in a clockwise sense. Such a progression is consistent with expected stress variations that would accompany plausible non-synchronous rotation. The more recent data provide a more complete record than the G1 data did. We find that to fit the sequence into a continual clockwise change of orientation would require at least 1000° (>5 cycles) of azimuthal rotation. If due to non-synchronous rotation of Europa, this result implies that we are seeing back further into the tectonic record than the G1 results had suggested. The three sets of orientations found by Geissler et al. now appear to have been spaced out over several cycles, not during a fraction of one cycle. While our more complete sequence of lineament formation is consistent with non-synchronous rotation, a statistical test shows that it cannot be construed as independent evidence. Other lines of evidence do support non-synchronous rotation, but azimuths of crack sequences do not show it, probably because only a couple of cracks form in a given region in any given non-synchronous rotation period.     

Mechanics of tidally driven fractures in Europa's ice shell

A fracture mechanics model is developed for the initiation and propagation of a crack through a porous ice layer of finite thickness under gravitational overburden. It is found that surface cracks generated in response to a tidally induced stress field may penetrate through the entire outer brittle layer if a subsurface ocean is present on Europa. Such penetration is found to be very unlikely in the absence of an ocean. A cycloidal crack would then form as a sequence of near instantaneous discrete failures, each extending roughly the brittle layer thickness in range, linked with a much lower apparent propagation speed set by the moving tidal stress field. The implications of this porous ice fracture model for ice-penetrating radar scattering loss and seismic activity are quantified.     

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.     

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.     

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.