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.     

Angular momentum of binary asteroids: Implications for their possible origin

We describe in this work a thorough study of the physical and orbital characteristics of extensively observed main-belt and trojan binaries, mainly taken from the LAOSA (Large Adaptive Optics Survey of Asteroids [Marchis, F., Baek, M., Berthier, J., Descamps, P., Hestroffer, D., Kaasalainen, M., Vachier, F., 2006c. In: Workshop on Spacecraft Reconnaissance of Asteroid and Comet Interiors. Abstract #3042]) database, along with a selection of bifurcated objects. Dimensionless quantities, such as the specific angular momentum and the scaled primary spin rate, are computed and discussed for each system. They suggest that these asteroidal systems might be the outcome of rotational fission or mass shedding of a parent body presumably subjected to an external torque. One of the most striking features of separated binaries composed of a large primary (Rp>100 km) with a much smaller secondary (Rs<20 km) is that they all have total angular momentum of ∼0.27. This value is quite close to the Maclaurin-Jacobi bifurcation (0.308) of a spinning fluid body. Alternatively, contact binaries and tidally locked double asteroids, made of components of similar size, have an angular momentum larger than 0.48. They compare successfully with the fission equilibrium sequence of a rotating fluid mass. In conclusion, we find that total angular momentum is a useful proxy to assess the internal structure of such systems.     

Distribution of icy particles across Enceladus' surface as derived from Cassini-VIMS measurements

The surface of Enceladus consists almost completely of water ice. As the band depths of water ice absorptions are sensitive to the size of particles, absorptions can be used to map variations of icy particles across the surface. The Visual and Infrared Mapping Spectrometer (VIMS) observed Enceladus with a high spatial resolution during three Cassini flybys in 2005 (orbits EN 003, EN 004 and EN 011). Based on these data we measured the band depths of water ice absorptions at 1.04, 1.25, 1.5, and 2 μm. These band depths were compared to water ice models that represent theoretically calculated reflectance spectra for a range of particle diameters between 2 μm and 1 mm. The agreement between the experimental (VIMS) and model values supports the assumption that pure water ice characterizes the surface of Enceladus and therefore that variations in band depth correspond to variations in water ice particle diameters. Our measurements show that the particle diameter of water ice increases toward younger tectonically altered surface units with the largest particles exposed in relatively “fresh” surface material. The smallest particles were generally found in old densely cratered terrains. The largest particles (∼0.2 mm) are concentrated in the so called “tiger stripes” at the south pole. In general, the particle diameters are strongly correlated with geologic features and surface ages, indicating a stratigraphic evolution of the surface that is caused by cryovolcanic resurfacing and impact gardening.     

HST spectral imaging of Titan's haze and methane profile between 0.6 and 1 μm during the 2000 opposition

We report on the analysis of high spatial resolution visible to near-infrared spectral images of Titan at Ls=240° in November 2000, obtained with the Space Telescope Imaging Spectrograph instrument on board the Hubble Space Telescope as part of program GO-8580. We employ a radiative transfer fractal particle aerosol model with a Bayesian parameter estimation routine that computes Titan's absolute reflectivity per pixel for 122 wavelengths by modeling the vertical distribution of the lower atmosphere haze and tropospheric methane. Analysis of these data suggests that Titan's haze concentration in the lower atmosphere varies in strength with latitude. We find Titan's tropospheric methane profile to be fairly consistent with latitude and longitude, and we find evidence for local areas of a CH₄-N₂ binary saturation in Titan's troposphere. Our results suggest that a methane and haze profile at one location on Titan would not be representative of global conditions.     

IRCS/Subaru observations of water in the inner coma of Comet 73P-B/Schwassmann-Wachmann 3: Spatially resolved rotational temperatures and ortho-para ratios

Comet 73P-B/Schwassmann-Wachmann 3 was observed with IRCS/Subaru at geocentric distance of 0.074 AU on UT 10 May 2006. Multiple H₂O emission lines were detected in non-resonant fluorescence near 2.9 μm. No significant variation in total H₂O production rate was found during the (3 h) duration of our observations. H₂O rotational temperatures and ortho-to-para abundance ratios were measured for several positions in the coma. The temperatures extracted from two different time intervals show very similar spatial distributions. For both, the rotational temperature decreased from ∼110 to ∼90 K as the projected distance from the nucleus increased from ∼5 to ∼30 km. We see no evidence for OPR change in the coma. The H₂O ortho-para ratio is consistent with the statistical equilibrium value (3.0) for all spatially resolved measurements. This implies a nuclear spin temperature higher than ∼45 K.     

The long-term stability of a possible aqueous ammonium sulfate ocean inside Titan

We model the thermal evolution of a subsurface ocean of aqueous ammonium sulfate inside Titan using a parameterized convection scheme. The cooling and crystallization of such an ocean depends on its heat flux balance, and is governed by the pressure-dependent melting temperatures at the top and bottom of the ocean. Using recent observations and previous experimental data, we present a nominal model which predicts the thickness of the ocean throughout the evolution of Titan; after 4.5 Ga we expect an aqueous ammonium sulfate ocean 56 km thick, overlain by a thick (176 km) heterogeneous crust of methane clathrate, ice I and ammonium sulfate. Underplating of the crust by ice I will give rise to compositional diapirs that are capable of rising through the crust and providing a mechanism for cryovolcanism at the surface. We have conducted a parameter space survey to account for possible variations in the nominal model, and find that for a wide range of plausible conditions, an ocean of aqueous ammonium sulfate can survive to the present day, which is consistent with the recent observations of Titan's spin state from Cassini radar data [Lorenz, R.D., Stiles, B.W., Kirk, R.L., Allison, M.D., del Marmo, P.P., Iess, L., Lunine, J.I., Ostro, S.J., Hensley, S., 2008. Science 319, 1649-1651].     

Diagnostics of Titan's stratospheric dynamics using Cassini/CIRS data and the 2-dimensional IPSL circulation model

The dynamics of Titan's stratosphere is discussed in this study, based on a comparison between observations by the CIRS instrument on board the Cassini spacecraft, and results of the 2-dimensional circulation model developed at the Institute Pierre-Simon Laplace, available at http://www.lmd.jussieu.fr/titanDbase [Rannou, P., Lebonnois, S., Hourdin, F., Luz, D., 2005. Adv. Space Res. 36, 2194-2198]. The comparison aims at both evaluating the model's capabilities and interpreting the observations concerning: (1) dynamical and thermal structure using temperature retrievals from Cassini/CIRS and the vertical profile of zonal wind at the Huygens landing site obtained by Huygens/DWE; and (2) vertical and latitudinal profiles of stratospheric gases deduced from Cassini/CIRS data. The modeled thermal structure is similar to that inferred from observations (Cassini/CIRS and Earth-based observations). However, the upper stratosphere (above 0.05 mbar) is systematically too hot in the 2D-CM, and therefore the stratopause region is not well represented. This bias may be related to the haze structure and to misrepresented radiative effects in this region, such as the cooling effect of hydrogen cyanide (HCN). The 2D-CM produces a strong atmospheric superrotation, with zonal winds reaching 200 m s⁻¹ at high winter latitudes between 200 and 300 km altitude (0.1-1 mbar). The modeled zonal winds are in good agreement with retrieved wind fields from occultation observations, Cassini/CIRS and Huygens/DWE. Changes to the thermal structure are coupled to changes in the meridional circulation and polar vortex extension, and therefore affect chemical distributions, especially in winter polar regions. When a higher altitude haze production source is used, the resulting modeled meridional circulation is weaker and the vertical and horizontal mixing due to the polar vortex is less extended in latitude. There is an overall good agreement between modeled chemical distributions and observations in equatorial regions. The difference in observed vertical gradients of C₂H₂ and HCN may be an indicator of the relative strength of circulation and chemical loss of HCN. The negative vertical gradient of ethylene in the low stratosphere at 15° S, cannot be modeled with simple 1-dimensional models, where a strong photochemical sink in the middle stratosphere would be necessary. It is explained here by dynamical advection from the winter pole towards the equator in the low stratosphere and by the fact that ethylene does not condense. Near the winter pole (80° N), some compounds (C₄H₂, C₃H₄) exhibit an (interior) minimum in the observed abundance vertical profiles, whereas 2D-CM profiles are well mixed all along the atmospheric column. This minimum can be a diagnostic of the strength of the meridional circulation, and of the spatial extension of the winter polar vortex where strong descending motions are present. In the summer hemisphere, observed stratospheric abundances are uniform in latitude, whereas the model maintains a residual enrichment over the summer pole from the spring cell due to a secondary meridional overturning between 1 and 50 mbar, at latitudes south of 40-50° S. The strength, as well as spatial and temporal extensions of this structure are a difficulty, that may be linked to possible misrepresentation of horizontally mixing processes, due to the restricted 2-dimensional nature of the model. This restriction should also be kept in mind as a possible source of other discrepancies.     

Spin rate distribution of small asteroids

The spin rate distribution of main belt/Mars crossing (MB/MC) asteroids with diameters 3-15 km is uniform in the range from f=1 to 9.5 d⁻¹, and there is an excess of slow rotators with f<1 d⁻¹. The observed distribution appears to be controlled by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. The magnitude of the excess of slow rotators is related to the residence time of slowed down asteroids in the excess and the rate of spin rate change outside the excess. We estimated a median YORP spin rate change of ≈0.022 d⁻¹/Myr for asteroids in our sample (i.e., a median time in which the spin rate changes by 1 d⁻¹ is ≈45 Myr), thus the residence time of slowed down asteroids in the excess is ≈110 Myr. The spin rate distribution of near-Earth asteroids (NEAs) with sizes in the range 0.2-3 km (∼5 times smaller in median diameter than the MB/MC asteroids sample) shows a similar excess of slow rotators, but there is also a concentration of NEAs at fast spin rates with f=9-10 d⁻¹. The concentration at fast spin rates is correlated with a narrower distribution of spin rates of primaries of binary systems among NEAs; the difference may be due to the apparently more evolved population of binaries among MB/MC asteroids.     

Coupling photochemistry with haze formation in Titan's atmosphere, Part II: Results and validation with Cassini/Huygens data

The new one-dimensional radiative-convective/photochemical/microphysical model described in Part I is applied to the study of Titan's atmospheric processes that lead to haze formation. Our model generates the haze structure from the gaseous species photochemistry. Model results are presented for the species vertical concentration profiles, haze formation and its radiative properties, vertical temperature/density profiles and geometric albedo. These are validated against Cassini/Huygens observations and other ground-based and space-borne measurements. The model reproduces well most of the latest measurements from the Cassini/Huygens instruments for the chemical composition of Titan's atmosphere and the vertical profiles of the observed species. For the haze production we have included pathways that are based on pure hydrocarbons, pure nitriles and hydrocarbon/nitrile copolymers. From these, the nitrile and copolymer pathways provide the stronger contribution, in agreement with the results from the ACP instrument, which support the incorporation of nitrogen in the pyrolized haze structures. Our haze model reveals a new second major peak in the vertical profile of haze production rate between 500 and 900 km. This peak is produced by the copolymer family used and has important ramifications for the vertical atmospheric temperature profile and geometric albedo. In particular, the existence of this second peak determines the vertical profile of haze extinction. Our model results have been compared with the DISR retrieved haze extinction profiles and are found to be in very good agreement. We have also incorporated in our model heterogeneous chemistry on the haze particles that converts atomic hydrogen to molecular hydrogen. The resultant H₂ profile is closer to the INMS measurements, while the vertical profile of the diacetylene formed is found to be closer to that of the CIRS profile when this heterogenous chemistry is included.     

The Venusian induced magnetosphere: A case study of plasma and magnetic field measurements on the Venus Express mission

Plasma and magnetic field measurements made onboard the Venus Express on June 1, 2006, are analyzed and compared with predictions of a global model. It is shown that in the orbit studied, the plasma and magnetic field observations obtained near the North Pole under solar minimum conditions were qualitatively and, in many cases also, quantitatively in agreement with the general picture obtained using a global numerical quasi-neutral hybrid model of the solar wind interaction (HYB-Venus). In instances where the orbit of Venus Express crossed a boundary referred to as the magnetic pileup boundary (MPB), field line tracing supports the suggestion that the MPB separates the region that is magnetically connected to the fluctuating magnetosheath field from a region that is magnetically connected to the induced magnetotail lobes.