On convection in ice I shells of outer Solar System bodies, with detailed application to Callisto

It has been argued that the dominant non-Newtonian creep mechanisms of water ice make the ice shell above Callisto's ocean, and by inference all radiogenically heated ice I shells in the outer Solar System, stable against solid-state convective overturn. Conductive heat transport and internal melting (oceans) are therefore predicted to be, or have been, widespread among midsize and larger icy satellites and Kuiper Belt objects. Alternatively, at low stresses (where non-Newtonian viscosities can be arbitrarily large), convective instabilities may arise in the diffusional creep regime for arbitrarily small temperature perturbations. For Callisto, ice viscosities are low enough that convection is expected over most of geologic time above the internal liquid layer for plausible ice grain sizes (?a few mm); the alternative for early Callisto, a conducting shell over a very deep ocean (>450 km), is not compatible with Callisto's present partially differentiated state. Moreover, if convection is occurring today, the stagnant lid would be quite thick (∼100 km) and compatible with the lack of active geology. Nevertheless, Callisto's steady-state heat flow may have fallen below the convective minimum for its ice I shell late in Solar System history. In this case convection ends, the ice shell melts back at its base, and the internal ocean widens considerably. The presence of such an ocean, of order 200 km thick, is compatible with Callisto's moment-of-inertia, but its formation would have caused an ∼0.25% radial expansion. The tectonic effects of such a late, slow expansion are not observed, so convection likely persists in Callisto, possibly subcritically. Ganymede, due to its greater size, rock fraction and full differentiation, has a substantially higher heat flow than Callisto and has not reached this tectonic end state. Titan, if differentiated, and Triton should be more similar to Ganymede in this regard. Pluto, like Callisto, may be near the tipping point for convective shutdown, but uncertainties in its size and rock fraction prevent a more definitive assessment.     

Models of the Internal Structure of Callisto

Models of the internal structure of Callisto were constructed and the extent of its differentiation was determined based on geophysical information from the Galileo spacecraft (the mass, the radius, the mean density, and the moment of inertia), geochemical data (the chemical composition of meteorites), and the equations of state of water, ices, and meteoritic material. The thickness and the phase state of the water-ice shell were defined as well as the ice concentrations in the rock-ice mantle and the bulk concentration of H₂O. The constraints on the density distribution in the mantle and the size of the rock-iron core were derived. We considered models of the internal structure of Callisto in which the presence of a continuous ice shell was assumed (models without ocean) and models with an internal ocean. We demonstrated that it is possible to apply three-layer models with an icy shell up to 320 km in thickness and a rock-iron core in different combinations with a rock-ice mantle. These models do not reject a two-layer structure of Callisto (an ice lithosphere plus a rock-ice mantle or a rock-ice mantle plus a rock-iron core) and a one-layer model of the satellite composed only of a rock-ice mantle with an ice concentration that is variable in depth. Taking into account the chemically bound water, the bulk content of H₂O in the satellite is found to be 49–55 wt %. For the model with an internal ocean, the geophysically allowed thickness of the water-ice shell of Callisto was estimated to be 270–315 km with thicknesses of the icy crust and the underlying water layer of 135–150 and 120–180 km, respectively. The results of reconstruction of the composition and structure of the regular satellites of Jupiter allow us to conclude that they were possibly formed from material whose composition was close to ordinary L/LL chondrites at relatively low temperatures, lower than the temperature of evaporation of iron and Fe-Mg silicates.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 321–341.Original Russian Text Copyright © 2005 by Kuskov, Kronrod.     

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.     

Internal structure of Europa and Callisto

Models of the internal structure of completely differentiated Europa and partially differentiated Callisto have been constructed on the basis of Galileo gravity measurements, geochemical constraints on composition of ordinary and carbonaceous chondrites, and thermodynamic data on the equations of state of water, high-pressure ices, and meteoritic material. We assume thermal and mechanical equilibrium for the interiors of the satellites. A geophysically and geochemically permissible thickness of Europa's outer water-ice shell lies between 105 and 160 km (6.2-9.2% of total mass). Our results show that the bulk composition of the rock-iron core of Europa may be described by material approaching the L/LL-type chondrites in composition, but cannot be correlated either with the material of CI chondrites or H chondrites. For Europa's L/LL-chondritic models, core radii are estimated to be 470-640 km (5.3-12.5% of total mass). The allowed thickness of Europa's H₂O layer ranges from 115±10 km for a differentiated L/LL-type chondritic mantle with a crust to 135±10 km for an undifferentiated mantle. We show that Callisto must only be partially differentiated into an outer ice-I layer, a water ocean, a rock-ice mantle, and a rock-iron core (mixture of anhydrous silicates and/or hydrous silicates + FeFeS alloy). We accept that the composition of the rock-iron material of Callisto is similar to the bulk composition of L/LL-type chondritic material containing up to 10-15% of iron and iron sulfide. Assuming conductive heat transfer through the ice-I crust [Ruiz, 2001. The stability against freezing of an internal liquid-water ocean on Gallisto. Nature, 412, 409-411], heat flows were estimated and the possibility of the existence of a water ocean in Callisto was evaluated. The liquid phase is stable (not freezing) beneath the ice crust, if the heat flow is between 3.3 and 3.7 mW m⁻², which corresponds to the heat flow from radiogenic sources. The thickness of the ice-I crust is 135-150 km, and that of the underlying water layer, 120-180 km. The results of modeling support the hypothesis that Callisto may have an internal liquid-water ocean. The allowed total (maximum) thickness of the outer water-ice shell is up to 270-315 km. Rock-iron core radii, depending on the presence or absence of hydrous silicates, do not exceed 500-700 km, the thickness of an intermediate ice-rock mantle is not less than 1400 km, and its density is in the range of 1960-2500 kg m⁻³. The surface temperature of Callisto is expected to be 100-112 K. The total amount of H₂O in Callisto is found to be 49-55 wt%. The correspondence between the density and moment of inertia values for bulk ice-free Io, rock-iron core of ice-poor Europa, and rock-iron cores of Ganymede and Callisto shows that their bulk compositions may be, in general, similar and may be described by the composition close to a material of the L/LL-type chondrites with the (Fe_tot/Si) weight ratios ranging from 0.9 to 1.3. Planetesimals composed of these types of ordinary chondrites could be considered as analogues of building material for the rock-iron cores of the Galilean satellites. Similarity of bulk composition of the rock-iron cores of the inner and outer satellites implies the absence of iron-silicon fractionation in the protojovian nebula.     

Internal structures of the Galilean satellites

New models for the interiors of Io, Ganymede, and Callisto are proposed. The model of Io consists of a thin, high-rigidity outer layer separated from a solid interior by a thin, molten or partially molten shell. The modulus of rigidity of the outer layer must be at least 100 times larger than that of the underlying partially molten shell. These layers have thicknesses of order 100 km or less. The near-surface partially molten layer was most likely produced early in Io's history as a consequence of accretional heating; enhanced tidal heating in the outer rigid layer has kept the underlying region partially molten to the present day. The model of Ganymede consists of an ice outer layer, a shell of undifferentiated, primordial ice-silicate mixture, and a rock core. Accretional heating is responsible for melting the ice in the outer layers of Ganymede's initially homogeneous ice-silicate interior. Most of the rock in this outer layer accumulates in a shell on top of Ganymede's early cold and rigid central region; the water in the outer layer quickly refreezes. Heating of the undifferentiated region by the decay of radioactive elements in the silicate fraction would gradually warm it and reduce its viscosity. The rock layer would become gravitationally unstable and sink through the undifferentiated materials to form a rock core. Callisto's heavily cratered surface strongly suggests that relatively little, if any, ice-rock differentiation has occured in its interior.     

Subsurface oceans and deep interiors of medium-sized outer planet satellites and large trans-neptunian objects

The detection of induced magnetic fields in the vicinity of the jovian satellites Europa, Ganymede, and Callisto is one of the most surprising findings of the Galileo mission to Jupiter. The observed magnetic signature cannot be generated in solid ice or in silicate rock. It rather suggests the existence of electrically conducting reservoirs of liquid water beneath the satellites' outermost icy shells that may contain even more water than all terrestrial oceans combined. The maintenance of liquid water layers is closely related to the internal structure, composition, and thermal state of the corresponding satellite interior. In this study we investigate the possibility of subsurface oceans in the medium-sized icy satellites and the largest trans-neptunian objects (TNO's). Controlling parameters for subsurface ocean formation are the radiogenic heating rate of the silicate component and the effectiveness of the heat transfer to the surface. Furthermore, the melting temperature of ice will be significantly reduced by small amounts of salts and/or incorporated volatiles such as methane and ammonia that are highly abundant in the outer Solar System. Based on the assumption that the satellites are differentiated and using an equilibrium condition between the heat production rate in the rocky cores and the heat loss through the ice shell, we find that subsurface oceans are possible on Rhea, Titania, Oberon, Triton, and Pluto and on the largest TNO's 2003 UB₃₁₃, Sedna, and 2004 DW. Subsurface oceans can even exist if only small amounts of ammonia are available. The liquid subsurface reservoirs are located deeply underneath an ice-I shell of more than 100 km thickness. However, they may be indirectly detectable by their interaction with the surrounding magnetic fields and charged particles and by the magnitude of a satellite's response to tides exerted by the primary. The latter is strongly dependent on the occurrence of a subsurface ocean which provides greater flexibility to a satellite's rigid outer ice shell.     

Implications from Galileo Observations on the Interior Structure and Chemistry of the Galilean Satellites

Data from the recent gravity measurements by the Galileo mission are used to construct wide ranges of interior structure and composition models for the Galilean satellites of Jupiter. These models show that mantle densities of Io and Europa are consistent with an olivine-dominated mineralogy with the ratios of Mg to Fe components depending on mantle temperature for Io and on ice shell thickness for Europa. The mantle density and composition depend relatively little on core composition. The size of the core is largely determined by the core's composition with core radius increasing with the concentration of a light component such as sulfur. For Io, the range of possible core sizes is between 38 and 53% of the satellite's radius. For Europa, there is also a substantial effect of the thickness of the ice layer which is varied between 120 and 170 km on the core size. Core sizes are between 10 and 45% of Europa's radius. The core size of Ganymede ranges between one-quarter and one-third of the surface radius depending on its sulfur content and the thickness of the ice shell. A subset of the Ganymede models is consistent with an olivine-dominated mantle mineralogy. The thickness of the silicate mantle above the core varies between 900 and 1100 km. The outermost ice shell is about 900 km in thickness and is further subdivided by pressure-induced phase transitions into ice I, ice III, ice V, and ice VI layers. Callisto should be differentiated, albeit incompletely. It is proposed that this satellite was never molten at a large scale but differentiated through the convective gradual unmixing of the ice and the metal/rock component. Bulk iron-to-silicon ratios Fe/Si calculated for the inner pair of satellites, Io and Europa, are less than the CI carbonaceous chondrite value of 1.7±0.1, whereas ratios for the outer pair, Ganymede and Callisto, cover a broad range above the chondritic value. Although the ratios are uncertain, in particular for Ganymede and Callisto, the values are sufficiently distinct to suggest a difference in composition between these two pairs of satellites. This may indicate a difference in iron-silicon fractionation during the formation of both classes of satellites in the protojovian nebula.     

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.     

Origin,Bulk Chemical Composition And Physical Structure Of The Galilean Satellites Of Jupiter: A Post-Galileo Analysis

The origin of Jupiter and the Galilean satellite system is examinedin the light of the new data that has been obtained by the NASA Galileo Project. In particular, special attention is given to a theory of satellite origin which was put forward at the start of the Galileo Mission and on the basis of which several predictions have now been proven successful (Prentice, 1996a–c). These predictions concern the chemical composition of Jupiter's atmosphere and the physical structure of the satellites. According to the proposed theory of satellite origin, each of the Galilean satellites formed by chemical condensation and gravitational accumulation of solid grains within a concentricfamily of orbiting gas rings. These rings were cast off equatorially by the rotating proto-Jovian cloud (PJC) which contracted gravitationally to form Jupiter some 4 billion years ago. The PJC formed from the gas and grains left over from the gas ring that had been shed at Jupiter's orbit by the contracting proto-solar cloud (PSC). Supersonic turbulentconvection provides the means for shedding discrete gas rings.The temperatures T_n of the system of gas rings shed by the PSCand PJC vary with their respective mean orbital radii R_n (n = 0, 1, 2, Ϊ ) according as T_n ∝ R_n ^-0.9. If the planet Mercury condenses at 1640 K, so accounting for the high density ofthat planet via a process of chemical fractionation between iron and silicates, then T_n at Jupiter's orbit is 158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock and ice, together with an enhancement in the abundance of solids relative to gas which takes place through gravitational sedimentation of solids onto the mean orbit of the gas ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. Io and Europa condense above the H₂O ice point and consist solely of hydrated rock (h-rock). The Ganymedan condensate consists of h-rock and H₂O ice. For Callisto, NH₃ ice makes up ∼5% of the condensate mass next to h-rock (∼50%) and H₂O ice (∼45%). Detailed thermal and structural models for each of Europa, Ganymedeand Callisto are constructed on the basis of the above initial bulk chemicalcompositions. For Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912M_E and a 150 km thick frozen mantle of salty H₂O yields a moment-of-inertiacoefficient which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone model possessing an inner core of solid FeS and mass ∼0.116M_G, and an outer H₂O ice mantle of mass ∼0.502M_G is needed to explain the gravity data.Ganymede's native magnetic field was formed by thermoremanent magnetization of Fe₃O₄. A new Callisto (C) model is proposed consisting of a core of mass 0.826M_C containing a uniform mixture of h-rock (60% by mass) and H₂O and NH₃ ices, and capped by a mantle of pure ice. This model may have the capacity to yield a thin layer of liquid NH₃ċ2H₂O at the core boundary, in line with Galileo's discovery of an induced magnetic field This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Fracture penetration in planetary ice shells

The proposed past eruption of liquid water on Europa and ongoing eruption of water vapor and ice on Enceladus have led to discussion about the feasibility of cracking a planetary ice shell. We use a boundary element method to model crack penetration in an ice shell subjected to tension and hydrostatic compression. We consider the presence of a region at the base of the ice shell in which the far-field extensional stresses vanish due to viscoelastic relaxation, impeding the penetration of fractures towards a subsurface ocean. The maximum extent of fracture penetration can be limited by hydrostatic pressure or by the presence of the unstressed basal layer, depending on its thickness. Our results indicate that Europa's ice shell is likely to be cracked under 1-3 MPa tension only if it is ?2.5 km thick. Enceladus' ice shell may be completely cracked if it is capable of supporting ∼1-3 MPa tension and is less than 25 km thick.     

The tidal response of Ganymede and Callisto with and without liquid water oceans

Calculations of the tidal responses of Ganymede and Callisto reveal that tidal amplitudes on these bodies may be as large as a few meters if a liquid ocean exists to decouple the surface ice from the interior. Tides on Ganymede's surface can exceed 7 m peak-to-peak variation, while on Callisto the tidal amplitude can exceed 5 m in the presence of a liquid ocean. Without an ocean, tidal amplitudes are less than 0.5 m on Ganymede and less than 0.3 m on Callisto. An orbiting spacecraft using an altimeter for crossover analysis and Doppler tracking from Earth should be able to achieve sufficient accuracy to identify the tidal amplitude to within about a meter over the course of a few months (observing tens of tidal cycles).     

Heat flow, lenticulae spacing, and possibility of convection in the ice shell of europa

Two opposing models to explain the geological features observed on Europa’s surface have been proposed. The thin-shell model states that the ice shell is only a few kilometers thick, transfers heat by conduction only, and can become locally thinner until it exposes an underlying ocean on the satellite’s surface. According to the thick-shell model, the ice shell may be several tens of kilometers thick and have a lower convective layer, above which there is a cold stagnant lid that dissipates heat by conduction. Whichever the case, from magnetic data there is strong support for the presence of a layer of salty liquid water under the ice. The present study was performed to examine whether the possibility of convection is theoretically consistent with surface heat flows of ∼100-200 mW m⁻², deduced from a thin brittle lithosphere, and with the typical spacing of 15-23 km proposed for the features usually known as lenticulae. It was obtained that under Europa’s ice shell conditions convection could occur and also account for high heat flows due to tidal heating of the convective (nearly isothermal) interior, but only if the dominant water ice rheology is superplastic flow (with activation energy of 49 kJ mol⁻¹; this is the rheology thought dominant in the warm interior of the ice shell). In this case the ice shell would be ∼15-50 km thick. Furthermore, in this scenario explaining the origin of the lenticulae related to convective processes requires ice grain size close to 1 mm and ice thickness around 15-20 km.     

Librations of the Galilean satellites: The influence of global internal liquid layers

The four Galilean satellites are thought to harbor one or even two global internal liquid layers beneath their surface layer. The iron core of Io and Ganymede is most likely (partially) liquid and also the core of Europa may be liquid. Furthermore, there are strong indications for the existence of a subsurface ocean in Europa, Ganymede, and Callisto. Here, we investigate whether libration observations can be used to prove the existence of these liquid layers and to constrain the thickness of the overlying solid layers. For Io, the presence of a small liquid core increases the libration of the mantle by a few percent with respect to an entirely solid Io and mantle libration observations could be used to determine the mantle thickness with a precision of several tens of kilometers given that the libration amplitude can be measured with a precision of 1 m. For Europa, Ganymede, and Callisto, the presence of a water ocean close to the surface increases by at least an order of magnitude the ice shell libration amplitude with respect to an entirely solid satellite. The shell libration depends essentially on the shell thickness and to a minor extent on the density difference between the ocean and the ice shell. The possible presence of a liquid core inside Europa and Ganymede has no noticeable influence on their shell libration. For a precision of several meters on the libration measurements, in agreement with the expected accuracy with the NASA/ESA EJSM orbiter mission to Europa and Ganymede, an error on the shell thickness of a few tens kilometers is expected. Therefore, libration measurements can be used to detect liquid layers such as Io’s core or water subsurface oceans in Europa, Ganymede, and Callisto and to constrain the thickness of the overlying solid surface layers.     

Thermal convection in ice-I shells of Titan and Enceladus

Cassini-Huygens observations have shown that Titan and Enceladus are geologically active icy satellites. Mitri and Showman [Mitri, G., Showman, A.P., 2005. Icarus 177, 447-460] and McKinnon [McKinnon, W.B., 2006. Icarus 183, 435-450] investigated the dynamics of an ice shell overlying a pure liquid-water ocean and showed that transitions from a conductive state to a convective state have major implications for the surface tectonics. We extend this analysis to the case of ice shells overlying ammonia-water oceans. We explore the thermal state of Titan and Enceladus ice-I shells, and also we investigate the consequences of the ice-I shell conductive-convective switch for the geology. We show that thermal convection can occur, under a range of conditions, in the ice-I shells of Titan and Enceladus. Because the Rayleigh number Ra scales with δ³/ηb, where δ is the thickness of the ice shell and ηb is the viscosity at the base of the ice-I shell, and because ammonia in the liquid layer (if any) strongly depresses the melting temperature of the water ice, Ra equals its critical value for two ice-I shell thicknesses: for relatively thin ice shell with warm, low-viscosity base (Onset I) and for thick ice shell with cold, high-viscosity base (Onset II). At Onset I, for a range of heat fluxes, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. Switches between these states can cause large, rapid changes in the ice-shell thickness. For Enceladus, we demonstrate that an Onset I transition can produce tectonic stress of ∼500 bars and fractures of several tens of km depth. At Onset II, in contrast, we demonstrate that zero equilibrium states exist for a range of heat fluxes. For a mean heat flux within this range, the satellite experiences oscillations in surface heat flux and satellite volume with periods of ∼50-800 Myr even when the interior heat production is constant or monotonically declining in time; these oscillations in the thermal state of the ice-I shell would cause repeated episodes of extensional and compressional tectonism.     

Formation of the galilean satellites in a gaseous nebula

A model for Galilean satellite formation was analyzed in which the satellites accrete in the presence of a dense, gaseous disk-shaped nebula and rapidly form optically thick, gravitationally bound primordial atmospheres. Upper-bound temperatures expected during accretion lead to partially differentiated structures for both Ganymede and Callisto, although with Ganymede much more differentiated than Callisto. When allowance is made for the aerodynamic breaking of infalling planetesimal fragments, lower surface temperatures result, and the amount of partial differentiation of Callisto is small, possibly approaching zero for a narrow size distribution of infalling planetesimals. The model is chosen to be consistent with the observed densities of the Galilean satellites and our current understanding of Jupiter formation. The retention of ices more volatile than H₂O is considered but not modeled in detail. A nominal nebula of ～0.1 Jupiter masses is constructed by consideration of likely surface density profiles and existing Jupiter collapse calculations. This nebula is optically thick (even if grain opacity is ignored) in both radial and vertical directions and has a temperature profile T ～ 3600 (R_J/R), where R_J is Jupiter's radius and R is the radial distance in the disk midplane. Satellites accrete very rapidly (dynamical time scales being 10²–10⁴ years) and their optically thick gaseous envelopes are unable to eliminate the heat of accretion by radiation. Water-saturated, convective, adiabatic envelopes form, through which planetesimals fall, break up, and partially disseminate their mass. The resulting satellite surface temperatures during accretion are calculated. Possible implications of these models for the subsequent evolution of Ganymede and Callisto are explored and it is suggested that the extensive differentiation undergone by Ganymede may provide the right environment for subsequent resurfacing, whereas the relative lack of extensive differentiation for Callisto may explain the inferred absence of endogenic tectonism.     

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.     

Origin and bulk chemical composition of the Galilean satellites and the primitive atmosphere of Jupiter: A pre-Galileo analysis

A theory for the origin and bulk chemical composition of the Galilean satellites is presented — to coincide with the start of the 2-year orbital tour of this satellite system by the Galileo Orbiter. The theory is based on the author's modern Laplacian theory of solar system origin (Prentice 1978a). The nub of the work reported here is that the Jupiter system is indeed a miniature planetary system that formed by much the same physical and chemical processes that were responsible for the condensation of the sun's own family of planets. In particular, a phenomenon of supersonic turbulent convection which I claim caused the proto-solar cloud to rid excess spin angular momentum, by shedding a concentric family of orbiting gas rings at the present planetary orbits, may also have operated with similar effect within the proto-Jovian cloud.Several predictions are made for the bulk chemical composition and physical structure of the icy Galilean satellites which, it is hoped, can be tested by the Galileo Orbiter. The mean density of Callisto is consistent with that of a chemically homogeneous body consisting of about 50% rock, 45% water ice, and 5% ammonia ice, incorporated as the hydrate NH₃·H₂O. Such a higher-than-solar mass abundance ratio of rock to ice arises naturally within the proto-Jovian cloud since (i) only 34% of the available H₂O vapor within the gas ring shed by the proto-solar cloud at Jupiter's orbit was condensed in solid form, and (ii) gravitational sedimentation of solids onto the mean orbit of the proto-solar gas ring leads to an enhancement in the heavy element fraction of the captured primitive Jovian atmosphere. All in all, I predict Jupiter's primitive atmosphere to be enhanced by a factor _en 2 in its rock mass fraction (including S) and by a factor 1.3 in its water content, relative to solar abundances. NH₃ and CH₄4 are present in almost solar proportions.Initially, Ganymede consisted of a chemically uniform mixture of rock and water ice in the proportions 0.524 : 0.476. The observed mean density of this satellite, however, lies midway between the mean densities expected for homogeneous and fully differentiated rock/ice bodies. The calculations presented here suggest that this body is about half-differentiated. I predict that the Galileo Orbiter will find the mean axial moment-of-inertia factor of Ganymede to be 0.35 ± 0.01.The circum-Jovian gas ring from which Europa condensed had a temperature of 302 K and a mean orbit gas pressure of 2.8 bar. Initially, this satellite consisted of a uniform mix of hydrated rocks, of which brucite Mg(OH)₂ was the principal constituent. The observed mean density of Europa coincides with that expected for this mix, provided that its 9.4% native H₂O content is now fractionated from the rock and resides at the satellite surface, forming a frozen mantle some 155 km thick. Regretfully, the mean density of Io cannot be matched by the solid composition reported here. Perhaps this satellite has a molten interior.     

Comment on “Mechanics of tidally driven fractures in Europa's ice shell” by S. Lee, R.T. Pappalardo, and N.C. Makris [2005. Icarus 177, 367-379]

We show Lee, Pappalardo, and Makris' [2005. Icarus 177, 367-379] argument that surface cracks in Europa's icy shell penetrate 3-10 times deeper in the presence of subsurface ocean is not correct. We use numerical calculations to demonstrate that there is at most 50% increase in penetration depth for a crack opening in a shell of finite thickness compared to a half-space. We also propose a simple equation based on force balances to estimate the maximum thickness of an ice shell that can be opened under tensile stress. Our calculations show that a crack can only penetrate 330-m-thick ice shell under 200 kPa far-field tensile stress and half of that if the stress is 100 kPa. But the presence of water would allow crack penetrate ∼4.0 km into the ice shell with zero porosity.     

Formation of the regular satellites of giant planets in an extended gaseous nebula I: subnebula model and accretion of satellites

We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect each giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet’s centrifugal radius (where the specific angular momentum of the collapsing giant planet gaseous envelope achieves centrifugal balance, located at rCJ ∼ 15R_J for Jupiter and rCS ∼ 22R_S for Saturn) and an optically thin, extended outer disk out to a fraction of the planet’s Roche-lobe (R_H), which we choose to be ∼R_H/5 (located at ∼150 R_J near the inner irregular satellites for Jupiter, and ∼200R_S near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk may result from the nebula gas flowing into the protoplanet during the time of giant planet gap-opening (or cessation of gas accretion). For the sake of specificity, we use a solar composition “minimum mass” model to constrain the gas densities of the inner and outer disks of Jupiter and Saturn (and also Uranus). Our model has Ganymede at a subnebula temperature of ∼250 K and Titan at ∼100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 and 90 K, respectively.Our model has Callisto forming in a time scale ∼10⁶ years, Iapetus in 10⁶-10⁷ years, Ganymede in 10³-10⁴ years, and Titan in 10⁴-10⁵ years. Callisto takes much longer to form than Ganymede because it draws materials from the extended, low density portion of the disk; its accretion time scale is set by the inward drift times of satellitesimals with sizes 300-500 km from distances ∼100R_J. This accretion history may be consistent with a partially differentiated Callisto with a ∼300-km clean ice outer shell overlying a mixed ice and rock-metal interior as suggested by Anderson et al. (2001), which may explain the Ganymede-Callisto dichotomy without resorting to fine-tuning poorly known model parameters. It is also possible that particulate matter coupled to the high specific angular momentum gas flowing through the gap after giant planet gap-opening, capture of heliocentric planetesimals by the extended gas disk, or ablation of planetesimals passing through the disk contributes to the solid content of the disk and lengthens the time scale for Callisto’s formation. Furthermore, this model has Hyperion forming just outside Saturn’s centrifugal radius, captured into resonance by proto-Titan in the presence of a strong gas density gradient as proposed by Lee and Peale (2000). While Titan may have taken significantly longer to form than Ganymede, it still formed fast enough that we would expect it to be fully differentiated. In this sense, it is more like Ganymede than like Callisto (Saturn’s analog of Callisto, we expect, is Iapetus). An alternative starved disk model whose satellite accretion time scale for all the regular satellites is set by the feeding of planetesimals or gas from the planet’s Roche-lobe after gap-opening is likely to imply a long accretion time scale for Titan with small quantities of NH₃ present, leading to a partially differentiated (Callisto-like) Titan. The Cassini mission may resolve this issue conclusively. We briefly discuss the retention of elements more volatile than H₂O as well as other issues that may help to test our model.