Titan trace gaseous composition from CIRS at the end of the Cassini-Huygens prime mission

This paper reports on the results from an extensive study of all nadir-looking spectra acquired by Cassini/CIRS during the 44 flybys performed in the course of the nominal mission (2004-2008). With respect to the previous study (Coustenis, A., and 24 colleagues [2007]. Icarus 189, 35-62, on flybys TB-T10) we present here a significantly richer dataset with, in particular, more data at high northern and southern latitudes so that the abundances inferred here at these regions are more reliable. Our enhanced high-resolution dataset allows us to infer more precisely the chemical composition of Titan all over the disk. We also include improved spectroscopic data for some molecules and updated temperature profiles. The latitudinal distributions of all of the gaseous species are inferred. We furthermore test vertical distributions essentially for acetylene (C₂H₂) from CIRS limb-inferred data and from current General Circulation Models for Titan and compare our results on all the gaseous abundances with predictions from 1-D photochemical-radiative models to check the reliability of the chemical reactions and pathways.     

The photochemical products of benzene in Titan’s upper atmosphere

The Cassini spacecraft detected benzene high in Titan’s atmosphere as well as the presence of large mass positive and negative ions. Previous work has suggested that these large mass ions could be composed of fused-ring polycyclic aromatic hydrocarbon compounds. These fused-ring PAHs, such as naphthalene and anthracene, are usually the result of high temperature processes that may not occur in Titan’s thin, cold, upper thermosphere. Here we suggest that a different class of aromatic compounds, polyphenyls, may be a better explanation of the data. Polyphenyls can grow to be large polymeric structures and could condense to form the aerosols seen in Titan’s cloud and hazes. They have similar properties to fused-ring PAHs (for example, electron affinity, ionization potential) and could be the negative ion species seen in the CAPS instrument data from the Cassini spacecraft.     

Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

Coupled convection and tidal dissipation in Europa’s ice shell

We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 10¹³ Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other features.     

Theory of the rotation of Janus and Epimetheus

The saturnian coorbital satellites Janus and Epimetheus present a unique dynamical configuration in the Solar System, because of high-amplitude horseshoe orbits, due to a mass ratio of order unity. As a consequence, they swap their orbits every 4 years, while their orbital periods is about 0.695 days. Recently, Tiscareno et al. (Tiscareno, M.S., Thomas, P.C., Burns, J.A. [2009]. Icarus 204, 254-261) got observational informations on the shapes and the rotational states of these satellites. In particular, they detected an offset in the expected equilibrium position of Janus, and a large libration of Epimetheus.We here propose to give a three-dimensional theory of the rotation of these satellites in using these observed data, and to compare it to the observed rotations. We consider the two satellites as triaxial rigid bodies, and we perform numerical integrations of the system in assuming the free librations as damped.The periods of the three free librations we get, associated with the three dimensions, are respectively 1.267, 2.179 and 2.098 days for Janus, and 0.747, 1.804 and 5.542 days for Epimetheus. The proximity of 0.747 days to the orbital period causes a high sensitivity of the librations of Epimetheus to the moments of inertia. Our theory explains the amplitude of the librations of Janus and the error bars of the librations of Epimetheus, but not an observed offset in the orientation of Janus.     

Hydrate sulfuric acid after sulfur implantation in water ice

For many years an ongoing research program performed at our laboratory has had the aim to investigate the implantation of reactive ions in ices relevant to planetology by using IR spectroscopy. We present new results obtained by implanting 200 keV sulfur ions into water ice at 80 K. We have looked at the formation of sulfur-bearing molecules such as sulfuric acid, sulfur dioxide and hydrogen sulfide. We find that hydrated sulfuric acid is formed with high yield (0.65±0.1 molecules/ion). An upper limit to the production yield of SO₂ (Y0.025 molecules/ion) has been estimated; no hydrogen sulfide has been detected. The formation of hydrogen peroxide is confirmed. Ozone is not detected. The results are discussed relevant to the inquiry on the radiolytic sulfur cycle considered responsible for the formation of sulfur-bearing molecules on the surfaces of the Galilean satellites. We demonstrate that sulfur implantation efficiently forms hydrated sulfuric acid whose observed abundance is explained as caused by an exogenic process. It is more difficult to say if the observed sulfur dioxide is quantitatively supported by only sulfur implantation; additional experimental studies are necessary along with direct observations, especially at UV wavelengths such as those that could be performed by instruments on board Hubble Space Telescope or by the forthcoming World Space Observatory (WSO/UV).     

Ice-silicate fractionation among icy bodies due to the difference of impact strength between ice and ice-silicate mixture

Laboratory experiments on the impact disruption of ice-silicate mixtures were conducted to clarify the accretion process of small icy bodies. Since the icy bodies are composed of ice and silicates with various porosities, we investigated the effect of porosity on the impact disruption of mixtures. We tested the mixture target with the mass ratio of ice to silicate, 0.5 and with 5 different porosities (0, 12.5, 25, 32, 37%) at the impact velocities of 150 to 670 m/s. The silicate mass ratio was changed from 0 to 0.5 in steps of 0.1 at a porosity of 12.5% and a constant impact velocity of about 300 m/s. The impact strength of the mixture was found to decrease with increasing porosity and the silicate mass ratio between 0.1 and 0.5 could enhance the strength of the icy target. The observed dependence of the impact strength on the porosity is opposite to that observed for pure ice. This difference could play an important role in ice-silicate fractionation during the accretion process. Because, ice rich bodies are easily broken as the porosity decreases in their evolution, the collisional growth could be prohibited. On the other hand, among the silicate rich bodies the collisional growth could be enhanced.     

Tidal pattern instabilities on multi-moon planets

The equilibrium tide-generating forces in the lunar orbital plane of a planet of radius R are calculated for the case of N moons of mass Mi orbiting the planet at instantaneous polar coordinates (Di, αi). For the case of a single moon, there are only two high tides. For the case of two moons, it is found that there can exist a critical lunar orbital distance at which two high tides become unstable with respect to formation of three high tides. Bifurcation diagrams are presented which depict how the angular positions of the high and low tides on the planet vary with the lunar distances and lunar separation angle. Tidal stability diagrams, which illustrate the stability regions for various tidal patterns as a function of lunar distances and lunar separation angle, are presented for various values of D₂/D₁ and M₂/M₁. Generally speaking, the aforementioned tidal instability, and hence the propensity for formation of three high tides on a two-moon planet, exists over a significant range of lunar distances and separation angles provided that M₂/M₁∼³(D₂/D₁). For the case of N>2 moons, the tidal stability diagram becomes more complex, revealing a diversity of potential tidal patterns.     

The formation of Ganymede's grooved terrain: Numerical modeling of extensional necking instabilities

Ganymede's grooved terrain likely formed during an epoch of global expansion, when unstable extension of the lithosphere resulted in the development of periodic necking instabilities. Linear, infinitesimal-strain models of extensional necking support this model of groove formation, finding that the fastest growing modes of an instability have wavelengths and growth rates consistent with Ganymede's grooves. However, several questions remain unanswered, including how nonlinearities affect instability growth at large strains, and what role instabilities play in tectonically resurfacing preexisting terrain. To address these questions we numerically model the extension of an icy lithosphere to examine the growth of periodic necking instabilities over a broad range of strain rates and temperature gradients. We explored thermal gradients up to 45 K km⁻¹ and found that, at infinitesimal strain, maximum growth rates occur at high temperature gradients (45 K km⁻¹) and moderate strain rates (10⁻¹³ s⁻¹). Dominant wavelengths range from 1.8 to 16.4 km (post extension). Our infinitesimal growth rates are qualitatively consistent with, but an order of magnitude lower than, previous linearized calculations. When strain exceeds ∼10% growth rates decrease, limiting the total amount of amplification that can result from unstable extension. This fall-off in growth occurs at lower groove amplitudes for high-temperature-gradient, thin-lithosphere simulations than for low-temperature-gradient, thick-lithosphere simulations. At large strains, this shifts the ideal conditions for producing large amplitude grooves from high temperature gradients to more moderate temperature gradients (15 K km⁻¹). We find that the formation of periodic necking instabilities can modify preexisting terrain, replacing semi-random topography up to 100 m in amplitude with periodic ridges and troughs, assisting the tectonic resurfacing process. Despite this success, the small topographic amplification produced by our model presents a formidable challenge to the necking instability mechanism for groove formation. Success of the necking instability mechanism may require rheological weakening or strain localization by faulting, effects not included in our analysis.