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.     

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.     

Empirical constraints on the salinity of the europan ocean and implications for a thin ice shell

Induced electrical currents within Europa inferred from Galileo spacecraft magnetometer instrument data have been interpreted as due to a salty europan ocean. Published compositional models for Europa's ocean, based on aqueous leaching of carbonaceous chondrites, range over five orders of magnitude in predicted magnesium sulfate concentrations. We combine the Galileo spacecraft magnetometer-derived oceanic conductivities and radio Doppler data-derived interior models with laboratory conductivity vs concentration data for both magnesium sulfate solutions and terrestrial seawater to determine empirically the range of salt concentrations permitted for Europa's ocean. Solutions for both a three-layer spherical model, and a five-layer half-space model, that satisfy current preferred best fits to magnetometer data imply high, near-saturation salt concentrations and require a europan ice shell of less than 15 km thick, with a best fit at 4 km ice thickness. Adding a conductive core and mantle has a negligible effect on the amplitude when ocean conductivities are greater than a few Siemens per meter. Similarly, we find that including a realistic ionosphere has a negligible effect. We examine the implications of these results for the subsurface habitability of Europa.     

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-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.     

Cometary Delivery of Biogenic Elements to Europa

Jupiter's moon Europa may harbor an ocean beneath its ice cover, but the composition of that ocean and the overlying ice is nearly entirely unknown. Regardless of uncertainties in models for Europa's formation, we estimate lower limits for Europa's inventory of biogenic elements (such as C, N, O, and P) by investigating the contribution to the inventory of impact events over Europa's geologic history. A series of high-resolution hydrocode simulations were carried out over a range of comet densities (1.1, 0.8, and 0.6 g/cm³, corresponding to porosities between 0 and 45%) and impact velocities (16, 21.5, 26.5, and 30.5 km/s). We found that at typical impact velocities on Europa most impactor material reaches escape velocity, and it is assumed to be lost from Europa. For a nonporous comet, some fraction (20% or higher) of the projectile is retained by Europa even at the highest impact velocity modeled, 30.5 km/s. For porous comets, however, a significant fraction of the projectile (above 25%) is retained only for the lowest impact velocity modeled, 16 km/s. Integrated over solar system history, this suggests that 1 to 10 Gt of carbon could have been successfully delivered to Europa's surface by impacts of large comets (around 1 km in diameter). This is a few times more carbon than is contained in the procaryotic biomass of the upper 200 meters of the Earth's oceans, but about 2 orders of magnitude less if the whole depth of the oceans is considered. Therefore, regardless of its initial formation conditions, Europa should have a substantial inventory of “biogenic” elements, with implications for the chemistry of its oceans, ice cover, and the possibility of life.     

A model for the interior structure, evolution, and differentiation of Callisto

The recently measured dimensionless moment of inertia (MoI) factor for Callisto of 0.3549±0.0042 (Anderson et al., 2001, Icarus, 153, 157-161) poses a problem: its value cannot be explained by a model in which Callisto is completely differentiated into an ice shell above a rock shell and an iron core such as its neighboring satellite Ganymede nor can it be explained by a model of a homogeneous, undifferentiated ice-rock satellite. We show that Callisto may be incompletely differentiated into an outer ice-rock shell in which the volumetric rock concentration is close to the primordial one at the surface and decreases approximately linearly with depth, an ice mantle mostly depleted of rock, and an about 1800 km rock-ice core in which the rock concentration is close to the close-packing limit. The ice-rock shell thickness depends on uncertain rheology parameters and the heat flow and can be roughly 50 to 150 km thick. We show that if Callisto accreted from a mix of metal bearing rock and ice and if the average size of the rocks was of the order of meters to tens of meters, then Callisto may have experienced a gradual, but still incomplete unmixing of the two components. An ocean in Callisto at a depth of 100-200 km is difficult to obtain if the ice is pure H₂O and if the ice-rock lithosphere is 100 km or more thick; a water ocean is more plausible for ice contaminated by ammonia, methane or salts; or for pure H₂O at a depth of 400-600 km.     

Dim prospects for radar detection of Europa's ocean

The thickness of a porous regolith on Europa is estimated to be ∼1 km. Provided that this regolith contains as little as 1% by volume of cavities with sizes comparable to the radar wavelength, its presence poses an insurmountable obstacle in the search for Europa's ocean by means of a sounding radar.     

L- and LL-Chondritic Models of the Chemical Composition of Io

A geologically admissible range for the density of Io's mantle, its chemical and mineral compositions in the system, and the size of an Fe–FeS core are determined by mathematical simulation based on the fitting of the calculated mass and moment of inertia of Io to the experimental data. It is demonstrated that the bulk chemical composition of Io (crust + mantle + core) most closely resembles the composition of ordinary L and LL chondrites. The densities of mineral assemblages in the mantle and the iron concentration in the core, calculated on the basis of the L- and LL-chondritic models, meet the geophysical constraints on the mass and the moment of inertia of Io. The core radius is estimated to be 590–640 km for the Fe core (8–10% of Io's mass) and 820–890 km for the Fe–FeS eutectic core (13.5–17%).     

The temperature of Europa's subsurface water ocean

A 100 km deep liquid water ocean probably underlies the icy exterior of Jupiter's satellite Europa. The long-term persistence of a liquid ocean beneath an ice shell presents a thermal conundrum: Is the temperature of the ocean equal to the freezing point of water at the bottom of the ice shell, or is it equal to the somewhat warmer temperature at which water attains its maximum density? We argue that most of the ocean is at the temperature of maximum density and that the bulk of the vigorously convecting ocean is separated from the bottom of the ice shell by a thin “stratosphere” of stably stratified water which is at the freezing point, and therefore buoyant. If Europa's subsurface water ocean is warm, it could explain the widespread geologic evidence for apparent melt-through events observed on its surface and may constrain the overall age of its surface.     

Energy Distributions for Desorption of Sodium and Potassium from Ice: The Na/K Ratio at Europa

The sputtering and decomposition of the surface of Europa by fast ions and electrons lead to the production of an atomosphere containing sodium and potassium atoms. Here time-of-flight energy distributions are measured for Na and K sputtered from a vapor-deposited ice by 200-eV electrons. These data are then used in a Monte Carlo simulation for alkalis in Europa's atmosphere. Na/K ratios versus distance from Europa are calculated and compared to the recent observations in the range 6 to 18 Europan radii from the surface. Normalizing to the observations, the Na/K ratio for the loss rates is ∼27 and the ratio for the average surface source rates is ∼20. These ratios are very different from the Na/K ratio at Io and are larger than the Na/K ratio suggested for Europa's putative subsurface ocean, consistent with fractionation on freezing and upwelling of ocean material.     

A search for temporal variability in the surface chemistry of the icy Galilean satellites

The purpose of this study was to determine if any temporal variability in the broadband UV spectral properties of the icy Galilean satellites exists, and if so, to characterize its spatial distribution as a function of longitude in order to attempt to correlate any temporal changes with satellite surface interactions with the space environment. The temporal time period examined is between 1978-1984 (referred to as the 1980s data) and those from 1995-1996 (referred to as the 1990s data). The plausible temporal characteristics detected appear to vary from one satellite to the next. For Europa possible temporal variations are concentrated on the leading, anti-jovian quadrant. Example broadband UV spectra show Europa's spectral slope decreases (and darkens) with time on the leading and anti-jovian hemispheres, but remains essentially constant with time on the trailing hemisphere. The data quality does not support any definitive temporal changes for Ganymede. Possible temporal changes seen in the Callisto data set are concentrated on the jovian hemisphere. Example broadband UV spectra for Callisto show no definitive change in slope with time. The hypothesis is that these temporal differences in UV spectral properties are caused by variations in the surface ice chemistry due to temporal variability in the space environment. It is postulated that the UV spectral changes suggested for Europa may be linked to changes in H₂O₂ concentrations, whereas the changes on Callisto may be linked to variability in SO₂ concentration.     

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.     

Highly siderophile element and 187Re‐187Os isotopic systematics of ungrouped achondrite Northwest Africa 7325: Evidence for complex planetary processes

The abundances of highly siderophile elements (HSE; including Re, Os, Ir, Ru, Pt, and Pd) and ¹⁸⁷Re‐¹⁸⁷Os isotopic systematics were determined for two fragments from ungrouped achondrite NWA 7325. Rhenium‐Os systematics are consistent with closed‐system behavior since formation or soon after. The abundances of the HSE were therefore largely unaffected by late‐stage secondary processes such as shock or terrestrial weathering. As an olivine gabbro cumulate, this meteorite has a bulk composition consistent with derivation from a body that produced a core, mantle, and crust. Also consistent with derivation from a body that produced a core, both fragments of NWA 7325 have HSE abundances that are highly depleted compared to bulk chondrites. One fragment has ~0.002× CI chondrite Ir and relative HSE abundances similar to bulk chondrites. The other fragment has ~0.0002× CI chondrite Ir and relative HSE abundances that are fractionated compared to bulk chondrites. The chondritic relative HSE abundances of the fragment characterized by higher HSE abundances most likely reflect the addition of exogenous chondritic material during or after crystallization by surface impacts. The HSE in the other fragment is likely more representative of the parent body crust. One formation model that can broadly account for the HSE abundances in this fragment is multiple episodes of low‐pressure metal‐silicate equilibration, followed by limited late accretion and mantle homogenization. Given the different HSE compositions of the two adjoining fragments, this meteorite provides an example of the overprint of global processes (differentiation and late accretion) by localized impact contamination.     

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.     

Thermal and topographic tests of Europa chaos formation models from Galileo E15 observations

Stereo topography of an area near Tyre impact crater, Europa, reveals chaos regions characterised by marginal cliffs and domical topography, rising to 100-200 m above the background plains. The regions contain blocks which have both rotated and tilted. We tested two models of chaos formation: a hybrid diapir model, in which chaos topography is caused by thermal or compositional buoyancy, and block motion occurs due to the presence of near-surface (1-3 km) melt; and a melt-through model, in which chaos regions are caused by melting and refreezing of the ice shell. None of the hybrid diapir models tested generate any melt within 1-3 km of the surface, owing to the low surface temperature. A model of ocean refreezing following melt-through gives effective elastic thicknesses and ice shell thicknesses of 0.1-0.3 and 0.5-2 km, respectively. However, for such low shell thicknesses the refreezing model requires implausibly large lateral density contrasts (50-100 kg m⁻³) to explain the elevation of the centres of the chaos regions. Although a global equilibrium ice shell thickness of ≈2 km is possible if Europa's mantle resembles that of Io, it is unclear whether local melt-through events are energetically possible. Thus, neither of the models tested here gives a completely satisfactory explanation for the formation of chaos regions. We suggest that surface extrusion of warm ice may be an important component of chaos terrain formation, and demonstrate that such extrusion is possible for likely ice parameters.     

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.     

Secondary and sesquinary craters on Europa

We address impact cratering on Io and Europa, with the emphasis on the origin of small craters on Europa as secondary to the primary impacts of comets on Io, Europa, and Ganymede. In passing we also address the origin of secondary craters generated by Zunil, a well-studied impact crater on Mars that is a plausible analog to impact craters on Io. At nominal impact rates, and taking volcanic resurfacing into account, we find that there should be 1.3 impact craters on Io, equally likely to be of any diameter between 100 m and 20 km. The corresponding model age of Europa's surface is between 60 and 100 Ma. This range of ages does not include a factor three uncertainty stemming from the uncertain sizes and numbers of comets. The mass of basaltic impact ejecta from Io to reach Europa is found to meet or exceed the micrometeoroid flux as a source of rock-forming elements to Europa's ice crust. To describe impact ejecta in more detail we adapt models for impact-generated spalls and Grady-Kipp fragments originally developed by Melosh. Our model successfully reproduces the observed size-number distributions of small craters on both Mars and Europa. However, the model predicts that a significant fraction of the 200-500 m diameter craters on Europa are not traditional secondary craters but are instead sesquinary craters caused by impact ejecta from Io that had gone into orbit about Jupiter. This prediction is not supported by observation, which implies that high speed spalls usually break up into smaller fragments that make smaller sesquinary craters. Iogenic basalts are also interesting because they provide stratigraphic horizons on Europa that in principle could be used to track historic motions of the ice, and they provide materials suitable to radiometric dating of Europa's surface.     

Chondritic models of 4 Vesta: Implications for geochemical and geophysical properties

Simple mass‐balance and thermodynamic constraints are used to illustrate the potential geochemical and geophysical diversity of a fully differentiated Vesta‐sized parent body with a eucrite crust (e.g., core size and density, crustal thickness). The results of this analysis are then combined with data from the howardite–eucrite–diogenite (HED) meteorites and the Dawn mission to constrain Vesta's bulk composition. Twelve chondritic compositions are considered, comprising seven carbonaceous, three ordinary, and two enstatite chondrite groups. Our analysis excludes CI and LL compositions as plausible Vesta analogs, as these are predicted to have a negative metal fraction. Second, the MELTS thermodynamic calculator is used to show that the enstatite chondrites, the CV, CK and L‐groups cannot produce Juvinas‐like liquids, and that even for the other groups, depletion in sodium is necessary to produce liquids of appropriate silica content. This conclusion is consistent with the documented volatile‐poor nature of eucrites. Furthermore, carbonaceous chondrites are predicted to have a mantle too rich in olivine to produce typical howardites and to have Fe/Mn ratios generally well in excess of those of the HEDs. On the other hand, an Na‐depleted H‐chondrite bulk composition is capable of producing Juvinas‐like liquids, has a mantle rich enough in pyroxene to produce abundant howardite/diogenite, and has a Fe/Mn ratio compatible with eucrites. In addition, its predicted bulk‐silicate density is within 100 kg m^?3 of solutions constrained by data of the Dawn mission. However, oxidation state and oxygen isotopes are not perfectly reproduced and it is deduced that bulk Vesta may contain approximately 25% of a CM‐like component. Values for the bulk‐silicate composition of Vesta and a preliminary phase diagram are proposed.