Cassini observations of Saturn's inner plasmasphere: Saturn orbit insertion results

We present new and definitive results of Cassini plasma spectrometer (CAPS) data acquired during passage through Saturn's inner plasmasphere by the Cassini spacecraft during the approach phase of the Saturn orbit insertion period. This analysis extends the original analysis of Sittler et al. [2005. Preliminary results on Saturn's inner plasmasphere as observed by Cassini: comparison with Voyager. Geophys. Res. Lett. 32, L14S07, doi:10.1029/2005GL022653] to L∼10 along with also providing a more comprehensive study of the interrelationship of the various fluid parameters. Coincidence data are sub-divided into protons and water group ions. Our revised analysis uses an improved convergence algorithm which provides a more definitive and independent estimate of the spacecraft potential Φ_SC for which we enforce the protons and water group ions to co-move with each other. This has allowed us to include spacecraft charging corrections to our fluid parameter estimations and allow accurate estimations of fluctuations in the fluid parameters for future correlative studies. In the appendix we describe the ion moments algorithm, and minor corrections introduced by not weighting the moments with sinθ term in Sittler et al. [2005] (Correction offset by revisions to instruments geometric factor). Estimates of the spacecraft potential and revised proton densities are presented. Our total ion densities are in close agreement with the electron densities reported by Moncuquet et al. [2005. Quasi-thermal noise spectroscopy in the inner magnetosphere of Saturn with Cassini/RPWS: electron temperatures and density. Geophys. Res. Lett. 32, L20S02, doi:10.1029/2005GL022508] who used upper hybrid resonance (UHR) emission lines observed by the radio and plasma wave science (RPWS) instrument. We show a positive correlation between proton temperature and water group ion temperature. The proton and thermal electron temperatures track each with both having a positive radial gradient. These results are consistent with pickup ion energization via Saturn's rotational electric field. We see evidence for an anti-correlation between radial flow velocity V_R and azimuthal velocity V_φ, which is consistent with the magnetosphere tending to conserve angular momentum. Evidence for MHD waves is also present. We show clear evidence for outward transport of the plasma via flux tube interchange motions with the radial velocity of the flow showing positive radial gradient with functional dependence for 4<L<10 (i.e., if we assume to be diffusive transport then D_LL∼D₀L¹¹ for fixed stochastic time step δt). Previous models with centrifugal transport have used D_LL∼D₀L³ dependence. The radial transport seems to begin at Enceladus’ L shell, L∼4, where we also see a minimum in the W⁺ ion temperature . For the first time, we are measuring the actual flux tube interchange motions in the magnetosphere and how it varies with radial distance. These observations can be used as a constraint with regard to future transport models for Saturn's magnetosphere. Finally, we evaluate the thermodynamic properties of the plasma, which are all consistent with the pickup process being the dominant energy source for the plasma.     

Sputtering of ice grains and icy satellites in Saturn's inner magnetosphere

Icy grains and satellites orbiting in Saturn's magnetosphere are immersed in a plasma that sputters their surfaces. This limits the lifetime of the E-ring grains and ejects neutrals that orbit Saturn until they are ionized and populate its magnetosphere. Here we re-evaluate the sputtering rate of ice in Saturn's inner magnetosphere using the recent Cassini data on the plasma ion density, temperature and composition [Sittler Jr., E.C., et al., 2007a. Ion and neutral sources and sinks within Saturn's inner magnetosphere: Cassini results. Planet. Space Sci. 56, 3-18.] and a recent summary of the relevant sputtering data for ice [Famá, M., Shi, J., Baragiola, R.A., 2008. Sputtering of ice by low-energy ions. Surf. Sci. 602, 156-161.]. Although the energetic (>10 keV) ion component at Saturn is much smaller than was assumed to be the case after Voyager [Jurac, S., Johnson, R.E., Richardson, J.D., Paranicas, C., 2001a. Satellite sputtering in Saturn's magnetosphere. Planet. Space Sci. 49, 319-326; Jurac, S., Johnson, R.E., Richardson, J.D., 2001b. Saturn's E ring and production of the neutral torus. Icarus 149, 384-396.], we show that the sputtering rates are sensitive to the temperature of the thermal plasma and are still robust, so that sputtering likely determines the lifetime of the grains in Saturn's tenuous E-ring.     

Saturn's magnetospheric interaction with Titan as defined by Cassini encounters T9 and T18: New results

We present new results of Cassini's T9 flyby with complementary observations from T18. Based on Cassini plasma spectrometer (CAPS) and Cassini magnetometer (MAG), compositional evidence shows the upstream flow for both T9 and T18 appears composed of light ions (H⁺ and H₂⁺), with external pressures ∼30 times lower than that for the earlier TA flyby where heavy ions dominated the magnetospheric plasma. When describing the plasma heating and sputtering of Titan's atmosphere, T9 and T18 can be considered interactions of low magnetospheric energy input. On the other hand, T5, when heavy ion fluxes are observed to be higher than typical (i.e., TA), represents the limiting case of high magnetospheric energy input to Titan's upper atmosphere. Anisotropy estimates of the upstream flow are 1<T_⊥/T_‖<3 and the flow is perpendicular to B, indicative of local picked up ions from Titan's H and H₂ coronae extending to Titan's Hill sphere radius. Beyond this distance the corona forms a neutral torus that surrounds Saturn. The T9 flyby unexpectedly resulted in observation of two “wake” crossings referred to as Events 1 and 2. Event 2 was evidently caused by draped magnetosphere field lines, which are scavenging pickup ions from Titan's induced magnetopause boundary with outward flux ∼2×10⁶ ions/cm²/s. The composition of this out flow is dominated by H₂⁺ and H⁺ ions. Ionospheric flow away from Titan with ion flux ∼7×10⁶ ion/cm²/s is observed for Event 1. In between Events 1 and 2 are high energy field aligned flows of magnetosphere protons that may have been accelerated by the convective electric field across Titan's topside ionosphere. T18 observations are much closer to Titan than T9, allowing one to probe this type of interaction down to altitudes ∼950 km. Comparisons with previously reported hybrid simulations are made.     

Production, ionization and redistribution of O2 in Saturn's ring atmosphere

Molecular oxygen produced by the decomposition of icy surfaces is ubiquitous in Saturn's magnetosphere. A model is described for the toroidal O₂ atmosphere indicated by the detection of and O⁺ over the main rings. The O₂ ring atmosphere is produced primarily by UV photon-induced decomposition of ice on the sunlit side of the ring. Because O₂ has a long lifetime and interacts frequently with the ring particles, equivalent columns of O₂ exist above and below the ring plane with the scale height determined by the local ring temperature. Energetic particles also decompose ice, but estimates of their contribution over the main rings appear to be very low. In steady state, the O₂ column density over the rings also depends on the relative efficiency of hydrogen to oxygen loss from the ring/atmosphere system with oxygen being recycled on the grain surfaces. Unlike the neutral density, the ion densities can differ on the sunlit and shaded sides due to differences in the ionization rate, the quenching of ions by the interaction with the ring particles, and the northward shift of the magnetic equator relative to the ring plane. Although O⁺ is produced with a significant excess energy, is not. Therefore, should mirror well below those altitudes at which ions were detected. However, scattering by ion-molecule collisions results in much larger mirror altitudes, in ion temperatures that go through a minimum over the B-ring, and in the redistribution of both molecular hydrogen and oxygen throughout the magnetosphere. The proposed model is used to describe the measured oxygen ion densities in Saturn's toroidal ring atmosphere and its hydrogen content. The oxygen ion densities over the B-ring appear to require either significant levels of UV light scattering or ion transmission through the ring plane.     

Characteristics of Saturn's polar atmosphere and auroral electrons derived from HST/STIS, FUSE and Cassini/UVIS spectra

Ultraviolet (UV) spectra of Saturn's aurora obtained with the Hubble Space Telescope Imaging Spectrograph (STIS), the Cassini Ultraviolet Imaging Spectrograph (UVIS) and the Far Ultraviolet Spectroscopic Explorer (FUSE) have been analyzed. Comparisons between the observed spectra and synthetic models of electron-excited H₂ have been used to determine various auroral characteristics. Far ultraviolet (FUV: 1200-1700 Å) STIS and UVIS spectra exhibit, below 1400 Å, weak absorption due to methane, with a vertical column ranging between 1.4×¹⁰¹⁵ and . Using the low-latitude Moses et al. [Moses, J.I., Bézard, B., Lellouch, E., Feuchtgruber, H., Gladstone, G.R., Allen, M., 2000. Icarus, 143, 244-298] atmospheric model of Saturn and an electron energy-H₂ column relationship, these methane columns are converted into the mean energy of the primary precipitating electrons, estimated to lie in the range 10-18 keV. This result is confirmed by the study of self-absorption with UVIS and FUSE extreme ultraviolet (EUV: 900-1200 Å) spectra. Below 1200 Å, it is seen that transitions connecting to the v^″<2 vibrational levels of the H₂ electronic ground state are partially self-absorbed by H₂ molecules overlying the auroral emission. Because of its low spectral resolution (∼5.5 Å), the UVIS EUV spectrum we analyzed does not allow us to unequivocally determine reasonable ranges of temperatures and H₂ columns. On the other hand, the high spectral resolution (∼0.2 Å) of the FUSE LiF1a and LiF2a EUV spectra we examined resolve the H₂ rotational lines and makes it possible to determine the H₂ temperature. The modeled spectrum best fitting the FUSE LiF1a observation reveals a temperature of 500 K and self-absorption by a H₂ vertical column of . When converted to energy of precipitating electrons, this H₂ column corresponds to primary electrons of ∼10 keV. The model that best fits the LiF2a spectrum is characterized by a temperature of 400 K and is not self-absorbed, making this segment ideal to determine the H₂ temperature at the altitude of the auroral emission. The latter value is in agreement with temperatures obtained from infrared polar spectra. Self-absorption is detectable in the LiF2a segment for H₂ columns exceeding , which sets the maximum mean energy determined from the FUSE observations to ∼15 keV. The total electron energy range of 10-18 keV deduced from FUV and EUV observations places the auroral emission peak between the 0.1 and 0.3 μbar pressure levels. These values should be seen as an upper limit, since most of the Voyager UVS spectra of Saturn's aurora examined by Sandel et al. [Sandel, B.R., Shemansky, D.E., Broadfoot, A.L., Holberg, J.B., Smith, G.R., 1982. Science 215, 548] do not exhibit methane absorption. The auroral H₂ emission is thus likely located above but close to the methane homopause. The H₂ auroral brightness in the 800-1700 Å bandwidth varies from 2.9 kR to 139 kR, comparable to values derived from FUV Faint Object Camera (FOC) and STIS images.     

Infrared observations of Saturn's rings by Cassini CIRS : Phase angle and local time dependence

Since the Saturn orbit insertion (SOI) of the Cassini spacecraft, in July 2004, the Cassini Composite Infrared Spectrometer (CIRS) has obtained a large number of thermal infrared spectra of Saturn's rings. Over the two and a half years of observations to date, ring temperatures were retrieved for a large range of unique geometries, inaccessible from Earth. Understanding their dependencies with phase angle and local time is a clue to understanding the thermal properties and dynamics of Saturn's ring particles.Azimuthal scans of rings, which have been obtained by CIRS at constant radial distance from the planet, have been planned to measure ring temperature variations with local hour angle. Over 47 azimuthal scans for Saturn's main rings (A, B, C and Cassini Division) have been retrieved to date, on both lit and unlit sides, at different phase angles and spacecraft elevations. The first measurements of the transient thermal episode of eclipse cooling in the planetary shadow have also been obtained for all three rings.In this paper, we present an overview of all azimuthal scans obtained by the Cassini/CIRS instrument so far and the dependencies of the temperature and the filling factor with the phase angle and the local hour angle. The ring temperature varies with longitude as the input heating flux coming from Saturn and the Sun changes. The decrease in temperature with the increasing phase angle on both the lit and the unlit sides and for most of the local time also suggests the presence of slowly rotating particles. The crossing of the planet's shadow generates drastic azimuthal variations in temperature, up to 20 K in the C ring. The strong anisotropy of emission observed outside the shadow between low and high phase angles decreases when ring particles cross the shadow, suggesting that particles are almost isothermal in the shadow. This suggests a thermal inertia associated with a rotating rate of particles low enough to have a thermal contrast on their surface.The temperature in the B ring is less sensitive to the phase angle effect on the lit side, suggesting that particles are close enough to form a flat layer at a scale larger than the particle's radius. On the unlit side, particles in the B ring are less sensitive to the lack of solar input than in the C ring or in the A ring. Azimuthal variations of the filling factor in the A ring are also detected with changing ring local time. This effect might be created by the presence of gravitational instabilities (wakes).     

Electric fields within the martian magnetosphere and ion extraction: ASPERA-3 observations

Observations made by the ASPERA-3 experiment onboard the Mars Express spacecraft found within the martian magnetosphere beams of planetary ions. In the energy (E/q)-time spectrograms these beams are often displayed as dispersive-like, ascending or descending (whether the spacecraft moves away or approach the planet) structures. A linear dependence between energy gained by the beam ions and the altitude from the planet suggests their acceleration in the electric field. The values of the electric field evaluated from ion energization occur close to the typical values of the interplanetary motional electric field. This suggests an effective penetration of the solar wind electric field deep into the martian magnetosphere or generation of large fields within the magnetosphere. Two different classes of events are found. At the nominal solar wind conditions, a ‘penetration’ occurs near the terminator. At the extreme solar wind conditions, the boundary of the induced magnetosphere moves to a more dense upper atmosphere that leads to a strong scavenging of planetary ions from the dayside regions.     

Models of magnetic field generation in partly stable planetary cores: Applications to Mercury and Saturn

A substantial part of Mercury's iron core may be stably stratified because the temperature gradient is subadiabatic. A dynamo would operate only in a deep sublayer. We show that such a situation arises for a wide range of values for the heat flow and the sulfur content in the core. In Saturn the upper part of the metallic hydrogen core could be stably stratified because of helium depletion. The magnetic field is unusually weak in the case of Mercury and unusually axisymmetric at Saturn. We study numerical dynamo models in rotating spherical shells with a stable outer region. The control parameters are chosen such that the magnetic Reynolds number is in the range of expected Mercury values. Because of its slow rotation, Mercury may be in a regime where the dipole contribution to the internal magnetic field is weak. Most of our models are in this regime, where the dynamo field consists mainly of rapidly varying higher multipole components. They can hardly pass the stable conducting layer because of the skin effect. The weak low-degree components vary more slowly and control the structure of the field outside the core, whose strength matches the observed field strength at Mercury. In some models the axial dipole dominates at the planet's surface and in others the axial quadrupole is dominant. Differential rotation in the stable layer, representing a thermal wind, is important for attenuating non-axisymmetric components in the exterior field. In some models that we relate to Saturn the axial dipole is intrinsically strong inside the dynamo. The surface field strength is much larger than in the other cases, but the stable layer eliminates non-axisymmetric modes. The Messenger and Bepi Colombo space missions can test our predictions that Mercury's field is large-scaled, fairly axisymmetric, and shows no secular variations on the decadal time scale.     

Saturn's latitudinal C2H2 and C2H6 abundance profiles from Cassini/CIRS and ground-based observations

Hydrocarbons in the upper atmosphere of Saturn are known, from Voyager, ground-based, and early Cassini results, to vary in emission intensity with latitude. Of particular interest is the marked increase in hydrocarbon line intensity near the south pole during southern summer, as the increased line intensity cannot be simply explained by the increased temperatures observed in that region since the variations between C₂H₂ and C₂H₆ emission in the south pole region are different. In order to measure the latitudinal variations of hydrocarbons in Saturn's southern hemisphere we have used 3 cm⁻¹ resolution Cassini CIRS data from 2006 and combined this with measurements from the ground in October 2006 at NASA's IRTF using Celeste, an infrared high-resolution cryogenic grating spectrometer. These two data sets have been used to infer the molecular abundances of C₂H₂ and C₂H₆ across the southern hemisphere in the 1-10 mbar altitude region. We find that the latitudinal acetylene profile follows the yearly average mean daily insolation except at the southern pole where it peaks in abundance. Near the equator (5° S) the C₂H₂ abundance at the 1.2 mbar level is (1.6±0.19)×¹⁰−7 and it decreases by a factor of 2.7 from the equator toward the pole. However, at the pole (∼87° S) the C₂H₂ abundance jumps to (1.8±0.3)×¹⁰−7, approximately the equatorial value. The C₂H₆ abundance near the equator at the 2 mbar level is (0.7±0.1)×¹⁰−5 and stays approximately constant until mid-latitudes where it increases gradually toward the pole, attaining a value of (1.4±0.4)×¹⁰−5 there. The increase in ethane toward the pole with the corresponding decrease in acetylene is consistent with southern hemisphere meridional winds [Greathouse, T.K., Lacy, J.H., Bézard, B., Moses, J.I., Griffith, C.A., Richter, M.J., 2005. Icarus 177, 18-31]. The localized increase in acetylene at the pole provides evidence that there is dynamical transport of hydrocarbons from the equator to the southern pole.     

Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r∼118,800 km, 4.5 m at r∼120,700 km) and the inner A Ring (10-15 m for r<127,000 km).     

Jovian plasma torus interaction with Europa: 3D hybrid kinetic simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa ( [Cassidy et al., 2007] and [Shematovich et al., 2005]). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma: (a) with O⁺⁺ ions; (b) with O⁺⁺ and S⁺⁺ ions. The majority of O₂ atmosphere is thermal with an extended cold population (Cassidy et al., 2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.     

Dynamics of Saturn's magnetodisk near Titan's orbit: Comparison of Cassini magnetometer observations from real and virtual Titan flybys

We analyze the variability of the ambient magnetospheric field along Titan's orbit at 20.3 Saturn radii. However, while our preceding study (Simon et al., 2010) focused on Cassini magnetometer observations from the 62 Titan flybys (TA-T62) between October 2004 and October 2009, the present work discusses magnetic field data that were collected near Titan's orbit when the moon was far away. In analogy to the observations during TA-T62, the magnetospheric fields detected during these 79 “virtual” Titan flybys are strongly affected by the presence of Saturn's bowl-shaped and highly dynamic magnetodisk current sheet. We therefore provide a systematic classification of the magnetic field observations as magnetodisk current sheet or lobe-type scenarios. Among the 141 (62 real+79 virtual) crossings of Titan's orbit between July 2004 and December 2009, only 17 encounters (9 real+8 virtual) took place within quiet, magnetodisk lobe-type fields. During another 50 encounters (21 real+29 virtual), rapid transitions between current sheet and lobe fields were observed around the moon's orbital plane. Most of the encounters (54=22 real+32 virtual) occurred when Titan's orbit was embedded in highly distorted current sheet fields, thereby invalidating the frequently applied idealized picture of Titan interacting with a homogeneous and stationary magnetospheric background field. The locations of real and virtual Titan flybys are correlated to each other. Each of the 62 real Titan flybys possesses at least one virtual counterpart that occurred shortly before or after the real encounter and at nearly the same orbital position. A systematic comparison between Cassini magnetometer observations from the real Titan flybys and their virtual companions suggests that there is no clear evidence of Titan exerting a significant level of control on the vertical oscillatory motion of the magnetodisk near its orbit.     

Carbon dioxide on the satellites of Saturn: Results from the Cassini VIMS investigation and revisions to the VIMS wavelength scale

Several of the icy satellites of Saturn show the spectroscopic signature of the asymmetric stretching mode of C-O in carbon dioxide (CO₂) at or near the nominal solid-phase laboratory wavelength of 4.2675 μm (2343.3 cm⁻¹), discovered with the Visible-Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft. We report here on an analysis of the variation in wavelength and width of the CO₂ absorption band in the spectra of Phoebe, Iapetus, Hyperion, and Dione. Comparisons are made to laboratory spectra of pure CO₂, CO₂ clathrates, ternary mixtures of CO₂ with other volatiles, implanted and adsorbed CO₂ in non-volatile materials, and ab initio theoretical calculations of CO₂ * nH₂O. At the wavelength resolution of VIMS, the CO₂ on Phoebe is indistinguishable from pure CO₂ ice (each molecule’s nearby neighbors are also CO₂) or type II clathrate of CO₂ in H₂O. In contrast, the CO₂ band on Iapetus, Hyperion, and Dione is shifted to shorter wavelengths (typically ∼4.255 μm (∼2350.2 cm⁻¹)) and broadened. These wavelengths are characteristic of complexes of CO₂ with different near-neighbor molecules that are encountered in other volatile mixtures such as with H₂O and CH₃OH, and non-volatile host materials like silicates, some clays, and zeolites. We suggest that Phoebe’s CO₂ is native to the body as part of the initial inventory of condensates and now exposed on the surface, while CO₂ on the other three satellites results at least in part from particle or UV irradiation of native H₂O plus a source of C, implantation or accretion from external sources, or redistribution of native CO₂ from the interior.The analysis presented here depends on an accurate VIMS wavelength scale. In preparation for this work, the baseline wavelength calibration for the Cassini VIMS was found to be distorted around 4.3 μm, apparently as a consequence of telluric CO₂ gas absorption in the pre-launch calibration. The effect can be reproduced by convolving a sequence of model detector response profiles with a deep atmospheric CO₂ absorption profile, producing distorted detector profile shapes and shifted central positions. In a laboratory blackbody spectrum used for radiance calibration, close examination of the CO₂ absorption profile shows a similar deviation from that expected from a model. These modeled effects appear to be sufficient to explain the distortion in the existing wavelength calibration now in use. A modification to the wavelength calibration for 13 adjacent bands is provided. The affected channels span about 0.2 μm centered on 4.28 μm. The maximum wavelength change is about 10 nm toward longer wavelength. This adjustment has implications for interpretation of some of the spectral features observed in the affected wavelength interval, such as from CO₂, as discussed in this paper.     

Saturn's icy satellites investigated by Cassini-VIMS: I. Full-disk properties: 350-5100 nm reflectance spectra and phase curves

Saturn's icy satellites are among the main scientific objectives of the Cassini-VIMS (Visual and Infrared Mapping Spectrometer) experiment. This paper contains a first systematic and comparative analysis of the full-disk spectral properties of Dione, Enceladus, Epimetheus, Hyperion, Iapetus, Mimas, Phoebe, Rhea and Tethys as observed by VIMS from July 2004 to June 2005. The disk integrated properties (350-5100 nm reflectance spectra and phase curves at 550-2232 nm) and images of satellites are reported and discussed in detail together with the observed geometry. In general, the spectra in the visible spectral range are almost featureless and can be classified according to the spectral slopes: from the bluish Enceladus and Phoebe to the redder Iapetus, Hyperion and Epimetheus. In the 1000-1300 nm range the spectra of Enceladus, Tethys, Mimas and Rhea are characterized by a negative slope, consistent with a surface largely dominated by water ice, while the spectra of Iapetus, Hyperion and Phoebe show a considerable reddening pointing out the relevant role played by darkening materials present on the surface. In between these two classes are Dione and Epimetheus, which have a flat spectrum in this range. The main absorption bands identified in the infrared are the 1520, 2020, 3000 nm H₂O/OH bands (for all satellites), although Iapetus dark terrains show mostly a deep 3000 nm band while the 1520 and 2020 nm bands are very faint. In this spectral range, the Iapetus spectrum is characterized by a strong reddening. The CO₂ band at 4260 nm and the Fresnel ice peak around 3100 nm are evident only on Hyperion, Phoebe and Iapetus. The phase curves at 550 and at 2232 nm are reported for all the available observations in the 0°-144° range; Rhea shows an opposition surge at visible wavelengths in the 0.5°-1.17° interval. The improvement on the retrieval of the full-disk reflectance spectra can be appreciated by a direct comparison with ground-based telescopic data available from literature. Finally, data processing strategies and recent upgrades introduced in the VIMS-V calibration pipeline (flat-field and destriping-despiking algorithm) are discussed in appendices.     

Prometheus and Pandora: masses and orbital positions during the Cassini tour

Hubble Space Telescope (HST) images of Prometheus and Pandora show longitude discrepancies of about 20° with respect to the Voyager ephemerides, with an abrupt change in mean motion at the end of 2000 (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07). These discrepancies are anti-correlated and arise from chaotic interactions between the two moons, occurring at interval of 6.2 yr, when their apses are anti-aligned (Goldreich and Rappaport, 2003a, Icarus 162, 391-399). This behavior is attributed to the overlap of four 121:118 apse-type mean motion resonances (Goldreich and Rappaport, 2003b, Icarus 166, 320-327). We study the Prometheus-Pandora system using a Radau-type integrator taking into account Saturn's oblateness up to and including terms in J₆, plus the effects of the major satellites. We first confirm the chaotic behavior of Prometheus and Pandora. By fitting the numerical integrations to the HST data (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07), we derive the satellite masses. The resulting GM values (with their standard 3-σ errors) for Prometheus and Pandora are respectively and . Using the nominal shape of the two moons (Thomas, 1989, Icarus 77, 248-274), we derive Prometheus and Pandora's densities, 0.40^+0.03_−0.07 and 0.49^+0.05_−0.09 g cm⁻³, respectively. Our numerical fits also enable us to constrain the time of the latest apse anti-alignment in 2000. Finally, using our fit, we predict the orbital positions of the two satellites during the Cassini tour, and provide a lower limit of the uncertainties due to chaos. These uncertainties amount to about 0.2° in mean longitude at the arrival of the Cassini spacecraft in July 2004, and to about 3° in 2008, at the end of the nominal tour.     

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.     

Photometric modeling of Saturn's rings: I. Monte Carlo method and the effect of nonzero volume filling factor

The scattering properties of particulate rings with volume filling factors in the interval D=0.001-0.3 are studied, with photometric Monte Carlo ray tracing simulations combining the advantages of direct (photons followed from the source) and indirect methods (brightness as seen from the observing direction). Besides vertically homogeneous models, ranging from monolayers to classical many-particle thick rings, particle distributions obtained from dynamical simulations are studied, possessing a nonuniform vertical profile and a power law distribution of particle sizes. Self-gravity is not included to assure homogeneity in planar directions. Our main goal is to check whether the moderately flattened ring models predicted by dynamical simulations (with central plane D>0.1) are consistent with the basic photometric properties of Saturn's rings seen in ground-based observations, including the brightening near zero phase angle (opposition effect), and the brightening of the B-ring with increasing elevation angle (tilt effect). Our photometric simulations indicate that dense rings are typically brighter in reflected light than those with D→0, due to enhanced single scattering. For a vertically illuminated layer of identical particles this enhancement amounts at intermediate viewing elevations to roughly 1+2D. Increased single scattering is also obtained for low elevation illumination, further augmented at low phase angles α by the opposition brightening when D increases: the simulated opposition effect agrees very well with the Lumme and Bowell (1981, Astron. J. 86, 1694-1704) theoretical formula. For large α the total intensity may also decrease, due to reduced amount of multiple scattering. For the low (α=13°) and high (α=155°) phase angle geometries analyzed in Dones et al. (1993, Icarus 105, 184-215) the brightness change for D=0.1 amounts to 20% and −17%, respectively. In the case of an extended size distribution, dynamical simulations indicate that the smallest particles typically occupy a layer several times thicker than the largest particles. Even if the large particles form a dynamically dense system, a narrow opposition peak can arise due to mutual shadowing among the small particles: for example, a size distribution extending about two decades can account for the observed about 1° wide opposition peak, solely in terms of mutual shadowing. The reduced width of the opposition peak for extended size distribution is in accordance with Hapke's (1986, Icarus 67, 264-280) treatment for semi-infinite layers. Due to vertical profile and particle size distribution, the photometric behavior is sensitive to the viewing elevation: this can account for the tilt-effect of the B-ring, as dense and thus bright central parts of the ring become better visible for larger elevation, whereas in the case of smaller elevation, mainly low volume density upper layers are visible. Since multiple scattering is not involved, the explanation works also for albedo well below unity. Inclusion of nonzero volume density helps also to model some of the Voyager observations. For example, the discrepancy between predicted and observed brightness at large phase angles for much of the A-ring (Dones et al., 1993, Icarus 105, 184-215) is removed when the enhanced low α single scattering and reduced large α multiple scattering is allowed for. Also, a model with vertical thickness increasing with saturnocentric distance offers at least a qualitative explanation for the observed contrast reversal between the inner and outer A-ring in low and high phase Voyager images. Differences in local size distribution and thus on the effective D may also account for the contrast reversal in resonance sites.     

Laboratory kinetics of C2H radical reactions with ethane, propane, and n-butane at T = 96-296 K: implications for Titan

The kinetics of the reactions of C₂H radical with ethane (k₁), propane (k₂), and n-butane (k₃) are studied over the temperature range of T = 96-296 K with a pulsed Laval nozzle apparatus that utilizes a pulsed laser photolysis-chemiluminescence technique. The C₂H decay profiles in the presence of both the alkane reactant and O₂ are monitored by the CH(A²Δ) chemiluminescence tracer method. The results, together with available literature data, yield the following Arrhenius expressions: k₁(T) = (0.51 ± 0.06) × 10⁻¹⁰ exp[(−76 ± 30)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-800 K), k₂(T) = (0.98 ± 0.32) × 10⁻¹⁰exp[(−71 ± 60)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-361 K), and k₃(T) = (1.23 ± 0.26) × 10⁻¹⁰ cm³ molecule⁻¹ s⁻¹ (T = 96-297 K). At T = 296 K, k₁ is measured as a function of total pressure and has little or no pressure dependence. The results from this work support a direct hydrogen abstraction mechanism for the title reactions. Implications to the atmospheric chemistry of Titan are discussed.     

Planetary,solar and astrophysical relativistic electrons: Common energization mechanisms?

Direct observations or deduced analysis indicate clearly that formation of intense fluxes of relativistic electrons is an important ingredient in the evolution of numerous active magnetized plasma systems. Examples of relativistic electron energization include the recovery phase of a planetary magnetic storm, post solar flare coronal activity and the afterglow of gamma ray bursts. It is suggested that there exists a universal mechanism, which may explain electron energization at the vastly different magnetized plasma environments. The favorite configuration consists of an inhomogeneous magnetic field anchored at a given magnetic structure and excitation of whistler waves due to external injection of low-energy non-isotropic electrons. The energization proceeds as a bootstrap process due to interaction with the propagating whistler waves along the inhomogeneous magnetic field.