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

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

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.     

Three-dimensional numerical modeling of refraction and reflection scattering from icy galilean satellites

Radar scattering from the icy galilean satellites is marked by unusually high backscatter cross sections and polarization ratios at wavelengths λ₀=3.5-70 cm. The persistence of exotic scattering behavior over this large a wavelength range suggests that the responsible mechanisms remain at least partially effective as the wavelength approaches or exceeds the size of individual scatterers. We examine two models previously analyzed in the geometrical optics limit—radar glory from buried craters (Eshleman, 1986, Science 234, 587-590) and refraction scattering from subsurface lenses (Hagfors et al., 1985, Nature 315, 637-640)—at wavelength scales using three-dimensional finite-difference time-domain (FDTD) numerical simulations. We include craters with rough walls and lenses with random inclusions of heterogeneous material. For hemispherical craters spanning up to 3λ₀ in diameter, we observe none of the exotic backscatter behavior attributed to the geometrical optics models. Nonspherical refraction scatterers can produce circular polarization ratios μ_C>1 and linear polarization ratios μ_L=0.5-0.8 at diameters as small as ∼λ₀, but the density of such inclusions must be high if refraction scattering alone is to account for the measured cross sections.     

Synthesis of hydrogen peroxide in water ice by ion irradiation

We present infrared absorption studies on the effects of 50-100 keV Ar⁺ and 100 keV H⁺ ion irradiation of water ice films at 20-120 K. The results support the view that energetic ions can produce hydrogen peroxide on the surface of icy satellites and rings in the outer Solar System, and on ice mantles on interstellar grains. The ion energies are characteristic of magnetospheric ions at Jupiter, and therefore the results support the idea that radiolysis by ion impact is the source of the H₂O₂ detected on Europa by the Galileo infrared spectrometer. We found that Ar⁺ ions, used to mimic S⁺ impacts, are roughly as efficient as H⁺ ions in producing H₂O₂, and that 100 keV H⁺ ions can produce hydrogen peroxide at 120 K. The synthesized hydrogen peroxide remained stable while warming the ice film after irradiation; the column density of the formed H₂O₂ is constant until the ice film begins to desorb, but the concentration of H₂O₂ increases with time during desorption because the water sublimes at a faster rate. Comparing the shape of the 3.5-μm absorption feature of H₂O₂ to the one measured on Europa shows excellent agreement in both shape and position, further indicating that the H₂O₂ detected on Europa is likely caused by radiolysis of water ice.     

Cassini ISS astrometric observations of the inner jovian satellites, Amalthea and Thebe

We present a total of 289 new astrometric observations of the inner jovian satellites, Amalthea and Thebe, obtained using the Cassini ISS narrow angle camera. Observations were made using image sequences from 2000 December 11-12 (inbound) and 2001 January 15-16 (outbound), at phase angles of approximately 2° and 122°, respectively. Target distances were of order 284 R_J, giving a maximum resolution of approximately 100 km/pixel. Centroided line and sample values for 239 observations of Amalthea and 50 of Thebe are provided, together with estimated camera pointing information for each image. Orbit fitting using a uniformly precessing Keplerian ellipse model, taking into account the oblateness of Jupiter up to terms in J₆, gave RMS fit residuals of 0.364 and 0.443 pixel for Amalthea and Thebe, respectively (equivalent to 0.450 and 0.547 arcsec). RMS residuals relative to the JPL JUP230 ephemeris were 0.306 and 0.604 pixel (equivalent to 0.378 and 0.746 arcsec), for Amalthea and Thebe. The fitted orbital parameters confirm the relatively high inclinations of these satellites (0.374°±0.002° and 1.076°±0.003°, respectively), equivalent to maximum vertical displacements above Jupiter's equatorial plane of 1188±6 and 4240±12 km, respectively, consistent with current estimates of the half-thicknesses of the Amalthea and Thebe gossamer rings [Ockert-Bell, M.E., Burns, J.A., Dauber, I.J., Thomas, P.C., Veverka, J., Belton, M.J.S., Klaasen, K.P., 1999. Icarus 138, 188-213].     

Constraints from deuterium on the formation of icy bodies in the Jovian system and beyond

We consider the role of deuterium as a potential marker of location and ambient conditions during the formation of small bodies in our Solar system. We concentrate in particular on the formation of the regular icy satellites of Jupiter and the other giant planets, but include a discussion of the implications for the Trojan asteroids and the irregular satellites. We examine in detail the formation of regular planetary satellites within the paradigm of a circum-Jovian subnebula. Particular attention is paid to the two extreme potential subnebulae—“hot” and “cold”. In particular, we show that, for the case of the “hot” subnebula model, the D:H ratio in water ice measured from the regular satellites would be expected to be near-Solar. In contrast, satellites which formed in a “cold” subnebula would be expected to display a D:H ratio that is distinctly over-Solar. We then compare the results obtained with the enrichment regimes which could be expected for other families of icy small bodies in the outer Solar system—the Trojan asteroids and the irregular satellites. In doing so, we demonstrate how measurements by Laplace, the James Webb Space Telescope, HERSCHEL and ALMA will play an important role in determining the true formation locations and mechanisms of these objects.     

A hydrocode calculation coupled with reaction kinetics of carbon compounds within an impact vapor plume and its implications for cometary impacts on Galilean satellites

The synthesis of organic molecules via chemical reactions within impact vapor plumes has been proposed as a mechanism to supply organics on a planet. However, the kinetics of chemical reactions within a rapidly expanding vapor plume or quenching process of the reactions has not been studied extensively. In this study, we constructed a new numerical model that calculates kinetics of the entire chemical reactions within an impact vapor plume. Numerical results revealed that the semi-analytical models proposed so far, in which the final amount of a chemical species was given by the equilibrium abundance at the quenching temperature of the fastest reaction path involving the species, underestimates the yield of organic molecules, such as HCN, by up to a factor of 10. This is because the previously used assumption that a species can achieve equilibrium with the rest of the reaction system via the fastest reaction path involving the species is not necessarily valid. Our analysis of the high-temperature H/C/N/O reaction system suggests that the quenching of slow reactions divides the reaction network into smaller reaction sub-systems isolated from the rest of the reaction system. Then, the fastest reaction path cannot equilibrate an isolated reaction sub-system with the rest of the reaction system. Simulation of this actual disequilibrium mechanism requires a simultaneous numerical calculation of the entire reaction network, which is equivalent to conducting a full kinetic model calculation, such as our model. Our numerical code makes it possible to discuss quantitatively the impact chemistry for various situations, such as the Galilean satellites. In this study, our numerical model is applied to the delivery of organic molecules via cometary impact on the Galilean satellites. Our numerical results indicate that small-particle impacts would produce HCN efficiently. Resulting HCN may freeze out immediately and be deposited on satellite surfaces, where it may be eventually converted into complex organics via irradiation of charged particle. On the other hand, large-size impacts may form transient CH₄-N₂ atmospheres, in which complex organics (tholin) may be formed via energy deposition of UV and/or charged particle. Resulting complex organics may subsequently precipitate on the satellite surfaces without clear correlation with the locations of impact craters. Such distribution of complex organics created by chemical reactions within vapor plumes due to cometary impacts may explain an absorption (4.57 μm) on Galilean satellites nonassociated with observable (moderate- and large-size) impact craters.     

Hydrogen peroxide dimers and the production of O2 in icy satellite surfaces

Oxygen is seen in the reflectance spectra of the icy surfaces of Ganymede and Europa via absorption bands at 627 and 573 nm. Here we show that the trapped O₂ associated with these spectral features can be produced via radiation-induced decomposition of hydrogen peroxide dimers contained in hydrogen peroxide inclusions in these icy satellite surfaces.     

The radiolysis of SO2 and H2S in water ice: Implications for the icy jovian satellites

Spectra of Europa, Ganymede, and Callisto reveal surfaces dominated by frozen water, hydrated materials, and minor amounts of SO₂, CO₂, and H₂O₂. These icy moons undergo significant bombardment by jovian magnetospheric radiation (protons, electrons, and sulfur and oxygen ions) which alters their surface compositions. In order to understand radiation-induced changes on icy moons, we have measured the mid-infrared spectra of 0.8 MeV proton-irradiated SO₂, H₂S, and H₂O-ice mixtures containing either SO₂ or H₂S. Samples with H₂O/SO₂ or H₂O/H₂S ratios in the 3-30 range have been irradiated at 86, 110, and 132 K, and the radiation half-lives of SO₂ and H₂S have been determined. New radiation products include the H₂S₂ molecule and HSO⁻₃, HSO⁻₄, and SO²⁻₄ ions, all with spectral features that make them candidates for future laboratory work and, perhaps, astronomical observations. Spectra of both unirradiated and irradiated ices have been recorded as a function of temperature, to examine thermal stability and phase changes. The formation of hydrated sulfuric acid in irradiated ice mixtures has been observed, along with the thermal evolution of hydrates to form pure sulfuric acid. These laboratory studies provide fundamental information on likely processes affecting the outer icy shells of Europa, Ganymede, and Callisto.     

Polarization opposition effect for the Galilean satellites of Jupiter

We present results of polarimetric observations of the Galilean satellites Io, Europa, Ganymede, and Callisto at phase angles ranging from 0.19° to 2.22°. The observations in the UBVR filters were performed using a one-channel photoelectric polarimeter attached to 70-cm telescope of the Chuguev Observation Station (Ukraine) on November 19-December 7, 2000. We have observed the polarization opposition effect for Io, Europa, and Ganymede to be a sharp secondary spike of negative polarization with an amplitude of about −0.4% centered at phase angles of 0.2°-0.7° and superimposed on the regular negative polarization branch. Although these minima for Io, Europa, and Ganymede show many similarities, they also exhibit a number of distinctions. The polarization opposition effect appears to be wavelength-dependent, at least for Europa and Ganymede. No polarization opposition effect was found for Callisto. The results obtained are discussed within the framework of different mechanisms of light scattering.     

H2O2 production by high-energy electrons on icy satellites as a function of surface temperature and electron flux

Chemistry on the icy surface of Europa is heavily influenced by the incident energetic particle flux from the jovian magnetosphere. The majority (>75%) of this energy is in the form of high energy electrons (extending to >10 MeV). We have simulated the electron irradiation environment of Europa with a vacuum system containing a high-energy electron gun for irradiation of ice samples formed on a gold mirror cooled with a cryostat. Pure water films of ∼2.6 μm thickness were grown at 100 K and then either cooled (to 80 K), warmed (to 120 K) or left at 100 K and subsequently irradiated with 10 keV electrons. The production of hydrogen peroxide (H₂O₂) was monitored by observation of the 2850 cm⁻¹ (3.5 μm) band. Equilibrium concentrations of H₂O₂, in units of percent by number H₂O₂ relative to water, were found to be 0.043% (80 K), 0.029% (100 K), and 0.0063% (120 K). These values are 33%, 22%, and 5%, respectively, that of the reported surface concentration on the leading hemisphere of Europa (Carlson, R.W., Anderson, M.S., Johnson, R.E., Smythe, W.D., Hendrix, A.R., Barth, C.A., et al. [1999]. Science 283(5410), 2062-2064) and less than the equilibrium concentrations formed by ion irradiation. In addition to the ice film temperature, the current of electrons was varied between different experiments to determine the production and destruction of H₂O₂ as a function of both electron flux and ice temperature. Variation in current was found to have little effect on the results other than accelerating arrival at radiolytic equilibrium.     

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.     

A gas-poor planetesimal capture model for the formation of giant planet satellite systems

Assuming that an unknown mechanism (e.g., gas turbulence) removes most of the subnebula gas disk in a timescale shorter than that for satellite formation, we develop a model for the formation of regular (and possibly at least some of the irregular) satellites around giant planets in a gas-poor environment. In this model, which follows along the lines of the work of Safronov et al. [1986. Satellites. Univ. of Arizona Press, Tucson, pp. 89-116], heliocentric planetesimals collide within the planet's Hill sphere and generate a circumplanetary disk of prograde and retrograde satellitesimals extending as far out as ∼RH/2. At first, the net angular momentum of this proto-satellite swarm is small, and collisions among satellitesimals leads to loss of mass from the outer disk, and delivers mass to the inner disk (where regular satellites form) in a timescale ?¹⁰⁵ years. This mass loss may be offset by continued collisional capture of sufficiently small <1 km interlopers resulting from the disruption of planetesimals in the feeding zone of the giant planet. As the planet's feeding zone is cleared in a timescale ?¹⁰⁵ years, enough angular momentum may be delivered to the proto-satellite swarm to account for the angular momentum of the regular satellites of Jupiter and Saturn. This feeding timescale is also roughly consistent with the independent constraint that the Galilean satellites formed in a timescale of ¹⁰⁵-¹⁰⁶ years, which may be long enough to accommodate Callisto's partially differentiated state [Anderson et al., 1998. Science 280, 1573; Anderson et al., 2001. Icarus 153, 157-161]. In turn, this formation timescale can be used to provide plausible constraints on the surface density of solids in the satellitesimal disk (excluding satellite embryos for satellitesimals of size ∼1 km), which yields a total disk mass smaller than the mass of the regular satellites, and means that the satellites must form in several ∼10 collisional cycles. However, much more work will need to be conducted concerning the collisional evolution both of the circumplanetary satellitesimals and of the heliocentric planetesimals following giant planet formation before one can assess the significance of this agreement. Furthermore, for enough mass to be delivered to form the regular satellites in the required timescale one may need to rely on (unproven) mechanisms to replenish the feeding zone of the giant planet. We compare this model to the solids-enhanced minimum mass (SEMM) model of Mosqueira and Estrada [2003a. Icarus 163, 198-231; 2003b. Icarus 163, 232-255], and discuss its main consequences for Cassini observations of the saturnian satellite system.     

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.     

Rate-dependent strength of porous ice–silica mixtures and its implications for the shape of small to middle-sized icy satellites

Porosity is one of the most important physical properties in the rheology of small icy satellites composed of ice–silicate mixtures. Deformation experiments involving ice and 1 μm silica bead mixtures were conducted to clarify the effect of porosity on the flow law of ice–silica mixtures. Mixtures with silica mass contents of 0, 30, and 50 wt.% were used for the experiments, and the porosity was changed from 0% to 25% in each mixture. The temperature ranged from −10 to −20 °C, and the strain rate was changed from 1.2 × 10⁻⁶ to 4.2 × 10⁻⁴ s⁻¹. As a result, it was found that the ice–silica mixtures deformed plastically, and that the relationship between the maximum stress, σ_max, on the stress–strain curve and the applied strain rate, , could be described by the following flow law: . The mixture became softer as the porosity or silica mass content increased, and the stress exponent n and activation energy Q were independent of porosity, depending only on the silica mass content. Furthermore, the parameter A₀ could be written as A₀ = B(1 − ?)^−α, where ? is the porosity. The constants B and α also depended only on the silica mass content, and they increased with the increase in this content. The Maxwell relaxation time was calculated in order to estimate the conditions for topographic relaxation of icy satellites, and it was found that topographic relaxation occurred at temperatures higher than 160 K in the case of icy satellites with mean radii of 200 km.     

The variability of volcanic activity at Zamama, Culann, and Tupan Patera on Io as seen by the Galileo Near Infrared Mapping Spectrometer

Zamama, Culann, and Tupan Patera are three large, persistent volcanic centers on the jovian moon Io. As part of an ongoing project to quantify contributions from individual volcanic centers to Io’s thermal budget, we have quantified the radiant flux from all suitable observations made by the Galileo Near Infrared Mapping Spectrometer (NIMS) of these volcanoes, in some cases filling omissions in previous analyses. At Zamama, after a long period of cooling, we see a peak in thermal emission that corresponds with new plume activity. Subsequently, toward the end of the Galileo epoch, thermal emission from Zamama drops off in a manner consistent with a greatly reduced eruption rate and the cooling of emplaced flows. Culann exhibits possible episodic activity. We present the full Tupan Patera NIMS dataset and derive new estimates of thermal output and temporal behavior. Eruption rates at these three volcanoes are on the order of 30 m³ s⁻¹, consistent with a previous analysis of NIMS observations of Prometheus, and nearly an order of magnitude greater than Kilauea volcano, Hawai’i, Earth’s most active volcano. We propose that future missions to the jovian system could better constrain activity at these volcanoes and others where similar styles of activity are taking place by obtaining data on a time scale of, ideally, at least one observation per day. Observations at similar or even shorter timescales are desirable during initial waxing phases of eruption episodes. These eruptions are identifiable from their characteristic spectral signatures and temporal behavior.     

CO2 production by ion irradiation of H2O ice on top of carbonaceous materials and its relevance to the Galilean satellites

In this work we report on new experiments of ion irradiation of water ice deposited on top of solid carbonaceous materials to study the production of CO₂ at the interface ice/refractory material and discuss the possibility that this mechanism accounts for the quantity of CO₂ ice detected on the surfaces of the Galilean satellites. The used experimental technique has been in situ infrared spectroscopy. We have irradiated thin films of H₂O frost on carbonaceous layers with 200 keV of He⁺ and Ar⁺, and 30 keV of He⁺ at 16 and 80 K. The used carbonaceous layers have been asphaltite, a natural bitumen, and solid organic residues obtained by irradiation of frozen benzene. In both cases the results show that CO₂ is produced very efficiently after irradiation obtaining a maximum quantity of the order of . These results are, also quantitatively similar, to those recently obtained for water ice deposited on amorphous carbon films [Mennella, V., Palumbo, M.E., Baratta, G.A., 2004. Formation of CO and CO₂ molecules by ion irradiation of water ice covered hydrogenated carbon grains. Astrophys. J. 615, 1073-1080]. Thus we suggest that, whatever is the carbonaceous residue, CO₂ will be produced efficiently by the studied process. These results have interest in the context of the surfaces of the icy Galilean satellites in which CO₂ has been detected mainly trapped in the non-ice material, not in the pure water ice. We suggest that radiolysis of mixtures of water ice and refractory carbonaceous materials is the primary formation mechanism responsible for the CO₂ formation on the surfaces of the Galilean satellites.     

Impact of electron chemistry on the structure and composition of Io's atmosphere

Two-dimensional model calculations (altitude and solar zenith angle) are performed to investigate the impact of electron chemistry on the composition and structure of Io's atmosphere. The calculations are based upon the model of Wong and Smyth (2000, Icarus 146, 60-74) for Io's SO₂ sublimation atmosphere with the addition of new electron chemistry, where the interactions of the electrons and neutrals are treated in a simple fashion. The model calculations are presented for Io's atmosphere at western elongation (dusk ansa) for both a low-density case (subsolar temperature of 113 K) and a high-density case (subsolar temperature of 120 K). The impact of electron-neutral chemistry on the composition and structure of Io's atmosphere is confined primarily to an interaction layer. The penetration depth of the interaction layer is limited to high altitudes in the thicker dayside atmosphere but reaches the surface in the thinner dayside and/or nightside atmosphere at larger solar zenith angles. Within most of the thicker dayside atmosphere, the column density of SO₂ is not significantly altered by electrons, but in the interaction layer all number densities are significantly altered: SO₂ is reduced, O, SO, S, and O₂ are greatly enhanced, and O, SO, and S become comparable to SO₂ at high altitudes. For the thinner nightside atmosphere, the species number densities are dramatically altered: SO₂ is drastically reduced to the least abundant species of the SO₂ family, SO and O₂ are significantly reduced at all altitudes, and O and S are dramatically enhanced and become the dominant species at all altitudes except near the surface. The interaction layer also defines the location of the emission layer for neutrals excited by electron impact and hence determines the fraction of the total neutral column density that is visible in remote observation. Electron chemistry may also impact the ratio of the equatorial to polar SO₂ column density deduced from Lyman-α images and the north-south alternating and System III longitude-dependent asymmetry observed in polar O and S emissions.     

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.     

A spectroscopic study of the surfaces of Saturn's large satellites: H2O ice, tholins, and minor constituents

We present spectra of Saturn's icy satellites Mimas, Enceladus, Tethys, Dione, Rhea, and Hyperion, 1.0-2.5 μm, with data extending to shorter (Mimas and Enceladus) and longer (Rhea and Dione) wavelengths for certain objects. The spectral resolution (R=λ/Δλ) of the data shown here is in the range 800-1000, depending on the specific instrument and configuration used; this is higher than the resolution (R=225 at 3 μm) afforded by the Visual-Infrared Mapping Spectrometer on the Cassini spacecraft. All of the spectra are dominated by water ice absorption bands and no other features are clearly identified. Spectra of all of these satellites show the characteristic signature of hexagonal H₂O ice at 1.65 μm. We model the leading hemisphere of Rhea in the wavelength range 0.3-3.6 μm with the Hapke and the Shkuratov radiative transfer codes and discuss the relative merits of the two approaches to fitting the spectrum. In calculations with both codes, the only components used are H₂O ice, which is the dominant constituent, and a small amount of tholin (Ice Tholin II). Tholin in small quantities (few percent, depending on the mixing mechanism) appears to be an essential component to give the basic red color of the satellite in the region 0.3-1.0 μm. The quantity and mode of mixing of tholin that can produce the intense coloration of Rhea and other icy satellites has bearing on its likely presence in many other icy bodies of the outer Solar System, both of high and low geometric albedos. Using the modeling codes, we also establish detection limits for the ices of CO₂ (a few weight percent, depending on particle size and mixing), CH₄ (same), and NH₄OH (0.5 weight percent) in our globally averaged spectra of Rhea's leading hemisphere. New laboratory spectral data for NH₄OH are presented for the purpose of detection on icy bodies. These limits for CO₂, CH₄, and NH₄OH on Rhea are also applicable to the other icy satellites for which spectra are presented here. The reflectance spectrum of Hyperion shows evidence for a broad, unidentified absorption band centered at 1.75 μm.