Negative ion chemistry in Titan's upper atmosphere

The Electron Spectrometer (ELS), one of the sensors making up the Cassini Plasma Spectrometer (CAPS) revealed the existence of numerous negative ions in Titan's upper atmosphere. The observations at closest approach (∼1000 km) show evidence for negatively charged ions up to ∼10,000 amu/q, as well as two distinct peaks at 22±4 and 44±8 amu/q, and maybe a third one at 82±14 amu/q. We present the first ionospheric model of Titan including negative ion chemistry. We find that dissociative electron attachment to neutral molecules (mostly HCN) initiates the formation of negative ions. The negative charge is then transferred to more acidic molecules such as HC₃N, HC₅N or C₄H₂. Loss occurs through associative detachment with radicals (H and CH₃). We attribute the three low mass peaks observed by ELS to CN⁻, C₃N⁻/C₄H⁻ and C₅N⁻. These species are the first intermediates in the formation of the even larger negative ions observed by ELS, which are most likely the precursors to the aerosols observed at lower altitudes.     

Cassini Langmuir probe measurements in the inner magnetosphere of Saturn

In the inner magnetosphere of Saturn, the plasma density and drift velocity are high enough, and the photoelectron current low enough, for a Langmuir probe to produce useful data on ion parameters. Plasma density and velocity are found by analyzing the current due to collected ions and emitted photoelectrons for a negative probe potential. In order to correctly analyze the data, the current caused by photoelectrons emitted from the probe must be known. For a spherical probe at negative bias this should be a constant current, but for Cassini's probe it varies with attitude. A likely cause of this is a leakage current from the stub to the probe. The plasma drift velocities derived from Langmuir probe measurements did not agree with those found by the Cassini plasma spectrometer in the inner magnetosphere, but did so elsewhere. A possible solution to this is a two-component plasma where the components have different drift velocities.     

What makes a planet habitable?

This work reviews factors which are important for the evolution of habitable Earth-like planets such as the effects of the host star dependent radiation and particle fluxes on the evolution of atmospheres and initial water inventories. We discuss the geodynamical and geophysical environments which are necessary for planets where plate tectonics remain active over geological time scales and for planets which evolve to one-plate planets. The discoveries of methane–ethane surface lakes on Saturn’s large moon Titan, subsurface water oceans or reservoirs inside the moons of Solar System gas giants such as Europa, Ganymede, Titan and Enceladus and more than 335 exoplanets, indicate that the classical definition of the habitable zone concept neglects more exotic habitats and may fail to be adequate for stars which are different from our Sun. A classification of four habitat types is proposed. Class I habitats represent bodies on which stellar and geophysical conditions allow Earth-analog planets to evolve so that complex multi-cellular life forms may originate. Class II habitats includes bodies on which life may evolve but due to stellar and geophysical conditions that are different from the class I habitats, the planets rather evolve toward Venus- or Mars-type worlds where complex life-forms may not develop. Class III habitats are planetary bodies where subsurface water oceans exist which interact directly with a silicate-rich core, while class IV habitats have liquid water layers between two ice layers, or liquids above ice. Furthermore, we discuss from the present viewpoint how life may have originated on early Earth, the possibilities that life may evolve on such Earth-like bodies and how future space missions may discover manifestations of extraterrestrial life.     

Characteristics of the dust-plasma interaction near Enceladus’ South Pole

We present RPWS Langmuir probe data from the third Enceladus flyby (E3) showing the presence of dusty plasma near Enceladus’ South Pole. There is a sharp rise in both the electron and ion number densities when the spacecraft traverses through Enceladus plume. The ion density near Enceladus is found to increase abruptly from about 10² cm⁻³ before the closest approach to 10⁵ cm⁻³ just 30 s after the closest approach, an amount two orders of magnitude higher than the electron density. Assuming that the inconsistency between the electron and ion number densities is due to the presence of dust particles that are collecting the missing electron charges, we present dusty plasma characteristics down to sub-micron particle sizes. By assuming a differential dust number density for a range in dust sizes and by making use of Langmuir probe data, the dust densities for certain lower limits in dust size distribution were estimated. In order to achieve the dust densities of micrometer and larger sized grains comparable to the ones reported in the literature, we show that the power law size distribution must hold down to at least 0.03 μm such that the total differential number density is dominated by the smallest sub-micron sized grains. The total dust number density in Enceladus’ plume is of the order of 10² cm⁻³ reducing to 1 cm⁻³ in the E-ring. The dust density for micrometer and larger sized grains is estimated to be about 10⁻⁴ cm⁻³ in the plume while it is about 10⁻⁶-10⁻⁷ cm⁻³ in the E-ring. Dust charge for micron sized grains is estimated to be about eight thousand electron charges reducing to below one hundred electron charges for 0.03 μm sized grains. The effective dusty plasma Debye length is estimated and compared with inter-grain distance as well as the electron Debye length. The maximum dust charging time of 1.4 h is found for 0.03 μm sized grains just 1 min before the closest approach. The charging time decreases substantially in the plume where it is only a fraction of a second for 1 μm sized grains, 1 s for 0.1 μm sized grains and about 10 s for 0.03 μm sized grains.     

Comparisons of Cassini flybys of the Titan magnetospheric interaction with an MHD model: Evidence for organized behavior at high altitudes

Recent papers suggest the significant variability of conditions in Saturn’s magnetosphere at the orbit of Titan. Because of this variability, it was expected that models would generally have a difficult time regularly comparing to data from the Titan flybys. However, we find that in contrast to this expectation, it appears that there is underlying organization of the interaction features roughly above ～1800 km (1.7 Rt) altitude by the average external field due to Saturn’s dipole moment. In this study, we analyze Cassini’s plasma and magnetic field data collected at 9 Titan encounters during which the external field is close to the ideal southward direction and compare these observations to the results from a 2-fluid (1 ion, 1 electron) 7-species MHD model simulations obtained under noon SLT conditions. Our comparative analysis shows that under noon SLT conditions the Titan plasma interaction can be viewed in two layers: an outer layer between 6400 and 1800 km where interaction features observed in the magnetic field are in basic agreement with a purely southward external field interaction and an inner layer below 1800 km where the magnetic field measurements show strong variations and deviate from the model predictions. Thus the basic features inferred from the Voyager 1 flyby seem to be generally present above ～1800 km in spite of the ongoing external variations from SLT excursions, time variability and magnetospheric current systems as long as a significant southward external field component is present. At around ～1800 km kinetic effects (such as mass loading and heavy ion pickup) and below 1800 km ionospheric effects (such as drag of ionospheric plasma due to coupling with neutral winds and/or magnetic memory of Titan’s ionosphere) complicate what is observed.     

The electromagnetic pickup of submicron-sized dust above Enceladus’s northern hemisphere

As the saturnian magnetoplasma sweeps past Enceladus, it experiences both a decrease in electron content and sharp slowdown in the northern hemisphere region within ～5 Enceladus Radii (R_e). This slowdown is observed by Cassini in regions not obviously associated with the southern directed plume-originating ions. We suggest herein that the decrease in northern hemisphere electron content and plasma slowdown could both be related to the presence of fine dust grains that are being accelerated by the Lorentz force created within the saturnian magnetic field system.     

Energetic electron observations of Rhea’s magnetospheric interaction

Saturn’s moon Rhea is thought to be a simple plasma absorber, however, energetic particle observations in its vicinity show a variety of unexpected and complex interaction features that do not conform with our current understanding about plasma absorbing interactions. Energetic electron data are especially interesting, as they contain a series of broad and narrow flux depletions on either side of the moon’s wake. The association of these dropouts with absorption by dust and boulders orbiting within Rhea’s Hill sphere was suggested but subsequently not confirmed, so in this study we review data from all four Cassini flybys of Rhea to date seeking evidence for alternative processes operating within the moon’s interaction region. We focus on energetic electron observations, which we put in context with magnetometer, cold plasma density and energetic ion data. All flybys have unique features, but here we only focus on several structures that are consistently observed. The most interesting common feature is that of narrow dropouts in energetic electron fluxes, visible near the wake flanks. These are typically seen together with narrow flux enhancements inside the wake. A phase-space-density analysis for these structures from the first Rhea flyby (R1) shows that Liouville’s theorem holds, suggesting that they may be forming due to rapid transport of energetic electrons from the magnetosphere to the wake, through narrow channels. A series of possibilities are considered to explain this transport process. We examined whether complex energetic electron drifts in the interaction region of a plasma absorbing moon (modeled through a hybrid simulation code) may allow such a transport. With the exception of several features (e.g. broadening of the central wake with increasing electron energy), most of the commonly observed interaction signatures in energetic electrons (including the narrow structures) were not reproduced. Additional dynamical processes, not simulated by the hybrid code, should be considered in order to explain the data. For the small scale features, the possibility that a flute (interchange) instability acts on the electrons is discussed. This instability is probably driven by strong gradients in the plasma pressure and the magnetic field magnitude: magnetometer observations show clearly signatures consistent with the (expected) plasma pressure loss due to ion absorption at Rhea. Another potential driver of the instability could have been gradients in the cold plasma density, which are, however, surprisingly absent from most crossings of Rhea’s plasma wake. The lack of a density depletion in Rhea’s wake suggests the presence of a local cold plasma source region. Hybrid plasma simulations show that this source cannot be the ionized component of Rhea’s weak exosphere. It is probably related to accelerated photoelectrons from the moon’s negatively charged surface, indicating that surface charging may play a very important role in shaping Rhea’s magnetospheric interaction region.     

Investigation of the force balance in the Titan ionosphere: Cassini T5 flyby model/data comparisons

Cassini’s Titan flyby on 16 April, 2005 (T5) is the only encounter when the two main ionizing sources of the moon’s atmosphere, solar radiation and corotating plasma, align almost anti-parallel. In this paper a single-fluid multi-species 3D MHD model of the magnetospheric plasma interaction for T5 conditions is analyzed. Model results are compared to observations to investigate the ionospheric dynamics at Titan as well as to understand the deviations from a typical solar wind interaction, such as Venus’ interaction with the solar wind. Model results suggest that for the T5 interaction configuration, corotating plasma is the dominant driver determining the global interaction features at high altitudes. In the lower ionosphere below ～1500 km altitude – where the control of the ionospheric composition transfers from dynamic to chemical processes – magnetic and thermal pressure gradients oppose each other locally, complicating the ionospheric dynamics. Model results also imply that the nightside ionosphere – produced only by the impact ionization in the model – does not provide enough thermal pressure to balance the incident plasma dynamic pressure. As a result, the induced magnetic barrier penetrates into the ionosphere by plasma convection down to ～1000 km altitude and by magnetic diffusion below this altitude. Moreover, strong horizontal drag forces due to ion-neutral collisions and comparable drag forces estimated from possible neutral winds in the lower ionosphere below ～1400 km altitude oppose over local regions, implying that the Titan interaction must be treated as a 3D problem. Ion and electron densities calculated from the model generally agree with the Cassini Ion Neutral Mass Spectrometer and Langmuir probe measurements; however, there are significant differences between the calculated and measured magnetic fields. We discuss possible explanations for the discrepancy in the magnetic field predictions.