Ammonium sulfate on Titan: Possible origin and role in cryovolcanism

We model the chemical evolution of Titan, wherein primordial NH₃ reacts with sulfate-rich brines leached from the silicate core during its hydration. The resulting differentiated body consists of a serpentinite core overlain by a high-pressure ice VI mantle, a liquid layer of aqueous ammonium sulfate, and a heterogeneous shell of methane clathrate, low-pressure ice Ih and solid ammonium sulfate. Cooling of the subsurface ocean results in underplating of the outer shell with ice Ih; this gravitationally unstable system can produce compositional plumes as ice Ih ascends buoyantly. Ice plumes may aid in advection of melt pockets through the shell and, in combination with surface topography, provide the necessary hydraulic pressure gradients to drive such melts to the surface. Moreover, contact between the magma and wall rock (methane clathrate) will allow some methane to dissolve in the magma, as well as eroding fragments of wall rock that can be transported as xenoliths. Upon rising to the clathrate decomposition depth (∼2 MPa, or 1700 m), the entrained xenoliths will break down to ice + methane gas, powering highly explosive eruptions with lava fountains up to several kilometers high. Hence we predict that Titan is being resurfaced by cryoclastic ash consisting of ice and ammonium sulfate (or its tetrahydrate), providing an abundance of sedimentary grains, a potential source of bedload for fluvial transport and erosion, and of sand-sized material for aeolian transport and dune-building. The infrared reflectance spectrum of ammonium sulfate makes it a plausible candidate for the 5 μm-bright material on Titan's surface.     

Thermal convection in the porous methane-soaked regolith of Titan: Investigation of stability

Titan is the only body, other than the Earth where liquid is present on the surface. In the present work we consider behavior of methane in the pores of Titan's regolith. Using numerical model we investigate quantitative conditions necessary for the onset of convection. We have found that the methane convection in Titan's regolith is possible. It can be expected in regions where the regolith has sufficiently high porosity, independently of the geothermal heat flux.     

Resurfacing of Titan by ammonia-water cryomagma

The Cassini Titan Radar Mapper observed on Titan several large features interpreted as cryovolcanic during the October 26, 2004 pass at high northern latitudes [Lopes, R.M.C., and 43 colleagues, 2007. Icarus 186, 395-412]. To date, models of ammonia-water resurfacing have not been tied to specific events or evolutionary stages of Titan. We propose a model of cryovolcanism that involves cracking at the base of the ice shell and formation of ammonia-water pockets in the ice. As these ammonia-water pockets undergo partial freezing in the cold ice shell, the ammonia concentration in the pockets increases, decreasing the negative buoyancy of the ammonia-water mixture. If the ice shell is contaminated by silicates delivered in impacts, the liquid-solid density difference would be even less. While the liquid cannot easily become buoyant relative to the surrounding ice, these concentrated ammonia-water pockets are sufficiently close to the neutral buoyancy point that large-scale tectonic stress patterns (tides, non-synchronous rotation, satellite volume changes, solid state convection, or subsurface pressure gradients associated with topography) would enable the ammonia to erupt effusively onto the surface. Rather than suggesting steady-state volcanism over the history of the Solar System, we favor a scenario where the cryovolcanic features could have been associated with episodic (potentially late) geological activity.     

Modelling of possible mud volcanism on Titan

Possible sedimentary basins on Titan are potential sites for the formation of mud volcanoes. In order to constrain the appearance of such features in remotely sensed imagery being acquired by the Cassini spacecraft, we have modelled the formation of mud volcanoes on Titan for a series of plausible mud compositions, climatic conditions and geological settings, as well as addressing the full range of eruption variables; mud viscosity, conduit diameter and eruption duration. We find that for an acetylene mud source containing 20 wt% liquid methane in pore spaces, overlain by a sheet of water ice 500-m thick, a mud volcano can grow as high as 140 m. Assuming reasonable eruption parameters, such an edifice may develop into a pancake-like dome several kilometres in diameter. If observed and properly characterised, mud volcanoes would provide an important window on the subsurface distribution and dynamics of solids and liquids in sedimentary basins on Titan.     

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.     

Coupling of thermal evolution and despinning of early Iapetus

The Cassini mission revealed two spectacular characteristics of Iapetus: (1) a geologically old and high equatorial ridge, which is unique in the Solar System and (2) a large flattening of 35 km consistent with the equilibrium figure for a hydrostatic body rotating with a period of 16 h, whereas the current spin period is 79.33 days. This study describes three-dimensional simulations of solid-state convection within an undifferentiated Iapetus. It investigates the implications for the evolution of the interior thermal structure and its spin rate and global shape using radially layered viscoelastic models. The role of the concentration in the short-lived radiogenic element [²⁶Al], just after accretion is completed, is specifically addressed. The first result is to show that whatever the [²⁶Al] value, convection occurs. As suggested by Castillo-Rogez et al. [Castillo-Rogez, J., Matson, D., Sotin, C., Johnson, T., Lunine, J., Thomas, P. [2007] Icarus, 190, 179-202], convection reduces the warming of the interior compared to the conductive evolution and therefore limits the conditions for despinning. In our calculations, two conceptual linear viscoelastic models are used. When considering a Maxwell rheology, the interior temperature (viscosity) never reaches a value high (low) enough to induce despinning. In order to promote dissipation at low temperature, a Burgers rheology, which includes an additional dissipation peak, is introduced. For favorable parameter values, this latter rheology leads to despinning. However, only models associated with large amounts of short-lived radiogenic elements lead to the observed flattening. This suggests that the accretion process needs to be completed shortly after the formation of CAIs (Calcium-Aluminum-rich Inclusions) (?4 Myr). For [²⁶Al] varying between 72 and 46 ppb, the observed flattening is obtained only for a limited range of initial spin period, between 9.5 and 10.2 h. For [²⁶Al] ranging between 30 and 15 ppb, initial spin rates smaller than 8.5 h are required. For smaller values of [²⁶Al], the body is too cold and viscous to acquire a significant flattening even if a rotation period close to the body disruption limit is considered. Even with a thin lithosphere during the early stage, our simulations show that Iapetus never reaches the equilibrium figure for a hydrostatic body due to the non-zero rigidity of the lithosphere. The 35 km value of the flattening is the result of the partial relaxation of an ancient larger flattening ranging between 45 and 80 km, depending on the evolution of the lithosphere thickness mainly controlled by the radiogenic content. A thin lithosphere is consistent with an early building of the equatorial ridge. The lithosphere thickening due to interior cooling can explain the preservation of the ridge throughout the remaining evolution of Iapetus.     

The effect of gravitational and pressure torques on Titan's length-of-day variations

Cassini radar observations show that Titan's spin is slightly faster than synchronous spin. Angular momentum exchange between Titan's surface and the atmosphere over seasonal time scales corresponding to Saturn's orbital period of 29.5 year is the most likely cause of the observed non-synchronous rotation. We study the effect of Saturn's gravitational torque and torques between internal layers on the length-of-day (LOD) variations driven by the atmosphere. Because static tides deform Titan into an ellipsoid with the long axis approximately in the direction to Saturn, non-zero gravitational and pressure torques exist that can change the rotation rate of Titan. For the torque calculation, we estimate the flattening of Titan and its interior layers under the assumption of hydrostatic equilibrium. The gravitational forcing by Saturn, due to misalignment of the long axis of Titan with the line joining the mass centers of Titan and Saturn, reduces the LOD variations with respect to those for a spherical Titan by an order of magnitude. Internal gravitational and pressure coupling between the ice shell and the interior beneath a putative ocean tends to reduce any differential rotation between shell and interior and reduces further the LOD variations by a few times. For the current estimate of the atmospheric torque, we obtain LOD variations of a hydrostatic Titan that are more than 100 times smaller than the observations indicate when Titan has no ocean as well as when a subsurface ocean exists. Moreover, Saturn's torque causes the rotation to be slower than synchronous in contrast to the Cassini observations. The calculated LOD variations could be increased if the atmospheric torque is larger than predicted and or if fast viscous relaxation of the ice shell could reduce the gravitational coupling, but it remains to be studied if a two order of magnitude increase is possible and if these effects can explain the phase difference of the predicted rotation variations. Alternatively, the large differences with the observations may suggest that non-hydrostatic effects in Titan are important. In particular, we show that the amplitude and phase of the calculated rotation variations are similar to the observed values if non-hydrostatic effects could strongly reduce the equatorial flattening of the ice shell above an internal ocean.     

Thermal interactions of the Huygens probe with the Titan environment: Constraint on near-surface wind

The Huygens probe lost heat to its cold environment during its descent through Titan's atmosphere and after landing. Here I report measurements of the probe's thermal behavior and comparison with ground tests (1) to provide a context for other scientific investigations, such as the release of volatiles from the landing site, and (2) to place constraints on Titan environmental parameters directly, such as the thermal conductivity of the surface material and the strength of winds at the surface. Near-surface winds are constrained to be less than 0.2 m s⁻¹, and probably much less.     

Heat flow and depth to a possible internal ocean on Triton

The Raz Fossae, a pair of ≈15-km wide trough en echelon interpreted as grabens, can be used to propose an estimation of the depth to the brittle-ductile transition on Triton. This estimation may in turn give an idea of the thermal state of Triton's icy lithosphere when these features formed. Given the young age of its surface, the conclusions obtained could be roughly applicable to the present state of this satellite of Neptune. Considering water or ammonia dihydrate as possible components of the lithosphere and a feasible range of strain rates, it was estimated that surface heat flow is greater than that inferred from radiogenic heating, especially for a lithosphere dominated by water. Also, an internal ocean could lie at a depth of only ∼20 km beneath the surface. The presence over the surface of an insulating layer of ice of low thermal conductivity (e.g., nitrogen) or of regolith would only substantially alter these estimates if the effective surface temperature were considerably higher than the observed value of 38 K.     

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.     

Convective-conductive transitions and sensitivity of a convecting ice shell to perturbations in heat flux and tidal-heating rate: Implications for Europa

We investigate the response of conductive and convective ice shells on Europa to variations of heat flux and interior tidal-heating rate. We present numerical simulations of convection in Europa's ice shell with Newtonian, temperature-dependent viscosity and tidal heating. Modest variations in the heat flux supplied to the base of a convective ice shell, ΔF, can cause large variations of the ice-shell thickness Δδ. In contrast, for a conductive ice shell, large ΔF involves relatively small Δδ. We demonstrate that, for a fluid with temperature-dependent viscosity, the heat flux undergoes a finite-amplitude jump at the critical Rayleigh number Racr. This jump implies that, for a range of heat fluxes relevant to Europa, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. We show that, as a result, modest variations in heat flux near the critical Rayleigh number can force the ice shell to switch between the thin, conductive and thick, convective configurations over a ∼10⁷-year interval, with thickness changes of up to ∼10-30 km. Depending on the orbital and thermal history, such switches might occur repeatedly. However, existing evolution models based on parameterized-convection schemes have to date not allowed these transitions to occur. Rapid thickening of the ice shell would cause radial expansion of Europa, which could produce extensional tectonic features such as fractures or bands. Furthermore, based on interpretations for how features such as chaos and ridges are formed, several authors have suggested that Europa's ice shell has recently undergone changes in thickness. Our model provides a mechanism for such changes to occur.     

The internal activity and thermal evolution of Earth-like planets

We model the internal thermal evolution of planets with Earth-like composition and masses ranging from 0.1 to 10 Earth masses over a period of 10 billion years. We also characterize the internal activity of the planets by the velocity of putative tectonic plates, the rate at which mantle material is processed through melting zones, and the time taken to process one mantle mass. The more massive the planet the larger its processing rate (?), which scales approximately as ?∝M0.8-1.0. The processing times for all the planets increase with time as they cool and become less active. As would be expected, the surface heat flow scales with planet mass. All planets have similar declines in mantle temperature except for the largest, in which pressure effects cause a larger decline. The larger planets have higher mantle temperatures over all times. The less massive the planet, the larger the decrease in core temperature with time. The core heat flow is also found to decrease more rapidly for smaller planet masses. Finally, rough predictions are made for the time required to generate an atmosphere from estimates of the time to degas water and carbon dioxide in mantle melting zones. The degassing times depend strongly on the initial temperature of the planet, but for the temperatures used in our model all the planets degas within ∼32 Ma after their formation.     

The effect of a dense atmosphere on the tidally induced potential of Titan

The effect of the dense atmosphere of Titan on the tidal variations of the external gravitational potential of degree two is quantified. The atmospheric tides perturb the external gravitational potential of Titan in two ways. First, the atmosphere itself contributes directly to the external gravitational potential with a period of 15.945 days. Second, the variable loading of the atmosphere induces mass redistribution within Titan, which also changes the external gravitational potential. It is shown that the relative atmospheric contributions to the tides are most likely less than 2% and vanish almost completely for the most plausible models with a subsurface ocean. This suggest that atmospheric tidal perturbations will contribute only negligibly to Cassini measurements of Titan's gravitational field so that the tidal Love numbers derived from these observations can be directly interpreted in terms of the satellite's interior.     

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.     

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.     

The orbital-thermal evolution and global expansion of Ganymede

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

Despinning plus global contraction and the orientation of lobate scarps on Mercury: Predictions for MESSENGER

Observations by the Mariner 10 spacecraft suggest that the lobate scarps on Mercury, which have been interpreted to record at most 1-2 km of radial contraction of the planet after the end of the Late Heavy Bombardment, possess a global, preferred N-S orientation but lack a strong latitudinal dependence on their surface expression. Here, we reexamine the idea that a decrease in the planetary rotation rate (despinning) coupled with global contraction of at least 3-5.5 km prior to the end of Late Heavy Bombardment resulted in global N-S oriented thrust faults. The surface expression of these faults is assumed to have been erased by the end of the Late Heavy Bombardment, and the faults were subsequently reactivated by later global contraction, producing generally N-S oriented thrust faults from an isotropic stress field. We use the estimate of >3-5.5 km contraction prior to ∼4 Ga as an additional constraint to thermomechanical simulations of the evolution of Mercury, finding that a wide range of models are consistent with this observation. The fact that a wide range of states are consistent with the contraction of Mercury prior to the end of Late Heavy Bombardment but only a restricted set of states are consistent with the at most 1-2 km of subsequent contraction bolsters the idea that there may be hidden strain on Mercury, features unseen by Mariner 10 but likely visible to the MESSENGER spacecraft.     

The rheology of cryovolcanic slurries: Motivation and phenomenology of methanol-water slurries with implications for Titan

The Cassini spacecraft has revealed landforms on the surface of Titan suggested to be viscous cryovolcanic flows and possibly eruptive domes. In order to relate those surface features to the processes and chemistries that produced them, it is necessary to construct flow models, which rely on characterization of the rheological properties of the eruptants. This paper describes our initial exploratory attempts to understand the rheological characteristics of cryogenic slurries, using a 40% methanol-water mixture, as a precursor to more detailed experiments. We have devised a new automated cryogenic rotational viscometer system to more fully characterize cryovolcanic slurry rheologies. A series of measurements were performed, varying first temperature, and then strain rate, which revealed development of yield stress-like behaviors, shear-rate dependence, and thixotropic behavior, even at relatively low crystal fractions, not previously reported. At fixed shear rate our data are fit well by the Andrade equation, with the activation energy modified by a solid volume fraction. At fixed temperature, depending on shearing history, a Cross model could describe our data over a wide shear rate range. A Bingham plastic model appears to be a good constitutive model for the data measured at high shear rates when the shear was global. The yield stress like behavior implies that levee formation on cryolava flows is more likely than would be inferred from the previous studies, and may provide a partial explanation for features interpreted as steep-sided volcanic constructs on Titan.     

Testing the survival of microfossils in artificial martian sedimentary meteorites during entry into Earth’s atmosphere: The STONE 6 experiment

If life ever appeared on Mars, could we find traces of primitive life embedded in sedimentary meteorites? To answer this question, a 3.5-byr-old volcanic sediment containing microfossils was embedded in the heat shield of a space capsule in order to test survival of the rock and the microfossils during entry into the Earth’s atmosphere (the STONE 6 experiment). The silicified volcanic sediment from the Kitty’s Gap Chert (Pilbara, Australia) is considered to be an excellent analogue for Noachian-age volcanic sediments. The microfossils in the chert are also analogues for potential martian life. An additional goal was to investigate the survival of living microorganisms (Chroococcidiopsis) protected by a 2-cm thick layer of rock in order to test whether living endolithic organisms could survive atmospheric entry when protected by a rocky coating.Mineralogical alteration of the sediment due to shock heating was manifested by the formation of a fusion crust, cracks in the chert due to prograde and retrograde changes of α quartz to β quartz, increase in the size of the fluid inclusions, and dewatering of the hydromuscovite-replaced volcanic protoliths. The carbonaceous microfossils embedded in the chert matrix survived in the rock away from the fusion crust but there was an increase in the maturity index of the kerogen towards the crust. We conclude that this kind of sediment can survive atmospheric entry and, if it contains microfossils, they could also survive. The living microorganisms were, however, completely carbonised by flame leakage to the back of the sample and therefore non-viable. However, using an analytical model to estimate the temperature reached within the sample thickness, we conclude that, even without flame leakage, the living organisms probably need to be protected by at least 5 cm of rock in order to be shielded from the intense heat of entry.