Eclipse reappearances of Io: Time-resolved spectroscopy (1.9-4.2 μm)

We obtained time-resolved, near-infrared spectra of Io during the 60-90 min following its reappearance from eclipse by Jupiter on five occasions in 2004. The purpose was to search for spectral changes, particularly in the well-known SO₂ frost absorption bands, that would indicate surface-atmosphere exchange of gaseous SO₂ induced by temperature changes during eclipse. These observations were a follow-on to eclipse spectroscopy observations in which Bellucci et al. [Bellucci et al., 2004. Icarus 172, 141-148] reported significant changes in the strengths of two strong SO₂ bands in data acquired with the VIMS instrument aboard the Cassini spacecraft. One of the bands (4.07 μm [ν₁ + ν₃]) observed by Bellucci et al. is visible from ground-based observatories and is included in our data. We detected no changes in Io’s spectrum at any of the five observed events during the approximately 60-90 min during which spectra were obtained following Io’s emergence from Jupiter’s shadow. The areas of the three strongest SO₂ bands in the region 3.5-4.15 μm were measured for each spectrum; the variation of the band areas with time does not exceed that which can be explained by the Io’s few degrees of axial rotation during the intervals of observation, and in no case does the change in band strength approach that seen in the Cassini VIMS data. Our data are of sufficient quality and resolution to show the weak 2.198 μm (4549.6 cm⁻¹) 4ν₁ band of SO₂ frost on Io for what we believe is the first time. At one of the events (June 22, 2004), we began the acquisition of spectra ∼6 min before Io reappeared from Jupiter’s shadow, during which time it was detected through its own thermal emission. No SO₂ bands were superimposed on the purely thermal spectrum on this occasion, suggesting that the upper limit to condensed SO₂ in the vertical column above Io’s surface was ∼4 × 10⁻⁵ g cm⁻².     

Cassini spectra and photometry 0.25-5.1 μm of the small inner satellites of Saturn

The nominal tour of the Cassini mission enabled the first spectra and solar phase curves of the small inner satellites of Saturn. We present spectra from the Visual Infrared Mapping Spectrometer (VIMS) and the Imaging Science Subsystem (ISS) that span the 0.25-5.1 μm spectral range. The composition of Atlas, Pandora, Janus, Epimetheus, Calypso, and Telesto is primarily water ice, with a small amount (∼5%) of contaminant, which most likely consists of hydrocarbons. The optical properties of the “shepherd” satellites and the coorbitals are tied to the A-ring, while those of the Tethys Lagrangians are tied to the E-ring of Saturn. The color of the satellites becomes progressively bluer with distance from Saturn, presumably from the increased influence of the E-ring; Telesto is as blue as Enceladus. Janus and Epimetheus have very similar spectra, although the latter appears to have a thicker coating of ring material. For at least four of the satellites, we find evidence for the spectral line at 0.68 μm that Vilas et al. [Vilas, F., Larsen, S.M., Stockstill, K.R., Gaffley, M.J., 1996. Icarus 124, 262-267] attributed to hydrated iron minerals on Iapetus and Hyperion. However, it is difficult to produce a spectral mixing model that includes this component. We find no evidence for CO₂ on any of the small satellites. There was a sufficient excursion in solar phase angle to create solar phase curves for Janus and Telesto. They bear a close similarity to the solar phase curves of the medium-sized inner icy satellites. Preliminary spectral modeling suggests that the contaminant on these bodies is not the same as the exogenously placed low-albedo material on Iapetus, but is rather a native material. The lack of CO₂ on the small inner satellites also suggests that their low-albedo material is distinct from that on Iapetus, Phoebe, and Hyperion.     

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.     

Production and detection of carbon dioxide on Iapetus

Cassini VIMS detected carbon dioxide on the surface of Iapetus during its insertion orbit. We evaluated the CO₂ distribution on Iapetus and determined that it is concentrated almost exclusively on Iapetus’ dark material. VIMS spectra show a 4.27-μm feature with an absorption depth of 24%, which, if it were in the form of free ice, requires a layer 31 nm thick. Extrapolating for all dark material on Iapetus, the total observable CO₂ would be 2.3 × 10⁸ kg.Previous studies note that free CO₂ is unstable at 10 AU over geologic timescales. Carbon dioxide could, however, be stable if trapped or complexed, such as in inclusions or clathrates. While complexed CO₂ has a lower thermal volatility, loss due to photodissociation by UV radiation and gravitational escape would occur at a rate of 2.6 × 10⁷ kg year⁻¹. Thus, Iapetus’ entire inventory of surface CO₂ could be lost within a few decades.The high loss/destruction rate of CO₂ requires an active source. We conducted experiments that generated CO₂ by UV radiation of simulated icy regolith under Iapetus-like conditions. The simulated regolith was created by flash-freezing degassed water, crushing it into sub-millimeter sized particles, and then mixing it with isotopically labeled amorphous carbon (¹³C) dust. These samples were placed in a vacuum chamber and cooled to temperatures between 50 K and 160 K. The samples were irradiated with UV light, and the products were measured using a mass spectrometer, from which we measured ¹³CO₂ production at a rate of 2.0 × 10¹² mol s⁻¹. Extrapolating to Iapetus and adjusting for the solar UV intensity and Iapetus’ surface area, we calculated that CO₂ production for the entire surface would be 1.1 × 10⁷ kg year⁻¹, which is only a factor of two less than the loss rate. As such, UV photochemical generation of CO₂ is a plausible source of the detected CO₂.     

Revisiting the thermal inertia of Iapetus: Clues to the thickness of the dark material

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m⁻² K⁻¹ s^−1/2 while the light material varies between 15 and 25 J m⁻² K⁻¹ s^−1/2. Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.     

Saturn’s icy satellites investigated by Cassini-VIMS: II. Results at the end of nominal mission

We report the detailed analysis of the spectrophotometric properties of Saturn’s icy satellites as derived by full-disk observations obtained by visual and infrared mapping spectrometer (VIMS) experiment aboard Cassini. In this paper, we have extended the coverage until the end of the Cassini’s nominal mission (June 1st 2008), while a previous paper (Filacchione, G., and 28 colleagues [2007]. Icarus 186, 259-290, hereby referred to as Paper I) reported the preliminary results of this study.During the four years of nominal mission, VIMS has observed the entire population of Saturn’s icy satellites allowing us to make a comparative analysis of the VIS-NIR spectral properties of the major satellites (Mimas, Enceladus, Tethys, Dione, Rhea, Hyperion, Iapetus) and irregular moons (Atlas, Prometheus, Pandora, Janus, Epimetheus, Telesto, Calypso, Phoebe). The results we discuss here are derived from the entire dataset available at June 2008 which consists of 1417 full-disk observations acquired from a variety of distances and inclinations from the equatorial plane, with different phase angles and hemispheric coverage. The most important spectrophotometric indicators (as defined in Paper I: I/F continua at 0.55 μm, 1.822 μm and 3.547 μm, visible spectral slopes, water and carbon dioxide bands depths and positions) are calculated for each observation in order to investigate the disk-integrated composition of the satellites, the distribution of water ice respect to “contaminants” abundances and typical regolith grain properties. These quantities vary from the almost pure water ice surfaces of Enceladus and Calypso to the organic and carbon dioxide rich Hyperion, Iapetus and Phoebe. Janus visible colors are intermediate between these two classes having a slightly positive spectral slope. These results could help to decipher the origins and evolutionary history of the minor moons of the Saturn’s system. We introduce a polar representation of the spectrophotometric parameters as function of the solar phase angle (along radial distance) and of the effective longitude interval illuminated by the Sun and covered by VIMS during the observation (in azimuth) to better investigate the spatial distribution of the spectrophotometric quantities across the regular satellites hemispheres. Finally, we report the observed spectral positions of the 4.26 μm band of the carbon dioxide present in the surface material of three outermost moons Hyperion, Iapetus and Phoebe.     

The near-IR spectrum of Titan modeled with an improved methane line list

We have obtained spatially resolved spectra of Titan in the near-infrared J, H and K bands at a resolving power of ∼5000 using the near-infrared integral field spectrometer (NIFS) on the Gemini North 8 m telescope. Using recent data from the Cassini/Huygens mission on the atmospheric composition and surface and aerosol properties, we develop a multiple-scattering radiative transfer model for the Titan atmosphere. The Titan spectrum at these wavelengths is dominated by absorption due to methane with a series of strong absorption band systems separated by window regions where the surface of Titan can be seen. We use a line-by-line approach to derive the methane absorption coefficients. The methane spectrum is only accurately represented in standard line lists down to ∼2.1 μm. However, by making use of recent laboratory data and modeling of the methane spectrum we are able to construct a new line list that can be used down to 1.3 μm. The new line list allows us to generate spectra that are a good match to the observations at all wavelengths longer than 1.3 μm and allow us to model regions, such as the 1.55 μm window that could not be studied usefully with previous line lists such as HITRAN 2008. We point out the importance of the far-wing line shape of strong methane lines in determining the shape of the methane windows. Line shapes with Lorentzian, and sub-Lorentzian regions are needed to match the shape of the windows, but different shape parameters are needed for the 1.55 μm and 2 μm windows. After the methane lines are modeled our observations are sensitive to additional absorptions, and we use the data in the 1.55 μm region to determine a D/H ratio of 1.77 ± 0.20 × 10⁻⁴, and a CO mixing ratio of 50 ± 11 ppmv. In the 2 μm window we detect absorption features that can be identified with the ν₅ + 3ν₆ and 2ν₃ + 2ν₆ bands of CH₃D.     

Hydrocarbons on Saturn's satellites Iapetus and Phoebe

Material of low geometric albedo (pV?0.1) is found on many objects in the outer Solar System, but its distribution in the saturnian satellite system is of special interest because of its juxtaposition with high-albedo ice. In the absence of clear, diagnostic spectral features, the composition of this low-albedo (or “dark”) material is generally inferred to be carbon-rich, but the form(s) of the carbon is unknown. Near-infrared spectra of the low-albedo hemisphere of Saturn's satellite Iapetus were obtained with the Visible-Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft at the fly-by of that satellite of 31 December 2004, yielding a maximum spatial resolution on the satellite's surface of ∼65 km. The spectral region 3-3.6 μm reveals a broad absorption band, centered at 3.29 μm, and concentrated in a region comprising about 15% of the low-albedo surface area. This is identified as the CH stretching mode vibration in polycyclic aromatic hydrocarbon (PAH) molecules. Two weaker bands attributed to CH₂ stretching modes in aliphatic hydrocarbons are found in association with the aromatic band. The bands most likely arise from aromatic and aliphatic units in complex macromolecular carbonaceous material with a kerogen- or coal-like structure, similar to that in carbonaceous meteorites. VIMS spectra of Phoebe, encountered by Cassini on 11 June 2004, also show the aromatic hydrocarbon band, although somewhat weaker than on Iapetus. The origin of the PAH molecular material on these two satellites is unknown, but PAHs are found in carbonaceous meteorites, cometary dust particles, circumstellar dust, and interstellar dust.     

Dione’s spectral and geological properties

We present a detailed analysis of the variations in spectral properties across the surface of Saturn’s satellite Dione using Cassini/VIMS data and their relationships to geological and/or morphological characteristics as seen in the Cassini/ISS images. This analysis focuses on a local region on Dione’s anti-saturnian hemisphere that was observed by VIMS with high spatial resolution during orbit 16 in October 2005. The results are incorporated into a global context provided by VIMS data acquired within Cassini’s first 50 orbits. Our results show that Dione’s surface is dominated by at least one global process. Bombardment by magnetospheric particles is consistent with the concentration of dark material and enhanced CO₂ absorption on the trailing hemisphere of Dione independent of the geology. Local regions within this terrain indicate a special kind of resurfacing that probably is related to large-scale impact process. In contrast, the enhanced ice signature on the leading side is associated with the extended ejecta of the fresh impact crater Creusa (∼49°N/76°W). Although no geologically active regions could be identified, Dione’s tectonized regions observed with high spatial resolution partly show some clean H₂O ice implying that tectonic processes could have continued into more recent times.     

Sublimation’s impact on temporal change of albedo dichotomy on Iapetus

Iapetus, one of the saturnian moons, has an extreme albedo contrast between the leading and trailing hemispheres. The origin of this albedo dichotomy has led to several hypotheses, however it remains controversial. To clarify the origin of the dichotomy, the key approach is to investigate the detailed distribution of the dark material. Recent studies of impact craters and surface temperature from Cassini spacecraft data implied that sublimation of H₂O ice can occur on Iapetus’ surface. This ice sublimation can change the albedo distribution on the moon with time.In this study, we evaluate the effect of ice sublimation and simulate the temporal change of surface albedo. We assume the dark material and the bright ice on the surface to be uniformly mixed with a certain volume fraction, and the initial albedo distribution to incorporate the dark material deposits on the surface. That is, the albedo at the apex is lowest and concentrically increases in a sinusoidal pattern. This situation simulates that dark materials existed around the Iapetus’ orbit billions of years ago, and the synchronously rotating Iapetus swept the material and then deposited it on its surface. The evolution of the surface albedo during 4.0 Gyr is simulated by estimating the surface temperature from the insolation energy on Iapetus including the effect of Saturn’s eccentricity and Iapetus’ obliquity precession, and evaluating the sublimation rate of H₂O ice from the Iapetus’ surface.As a result, we found that the distribution of the surface albedo changed dramatically after 4.0 Gyr of evolution. The sublimation has three important effects on the resultant surface albedo. First, the albedo in the leading hemisphere has significantly decreased to approach the minimum value. Second, the albedo distribution has been elongated along the equator. Third, the edge of the low albedo region has become clear. Considering the effect of ice sublimation, the current albedo distribution can be reconstructed from the sinusoidal albedo distribution, suggesting the apex-antapex cratering asymmetry as a candidate for the origin of the albedo dichotomy. From the model analysis, we obtained an important aspect that the depth of the turn-over layer where the darkening process proceeded for 4 Gyr should be an order of 10 cm, which is consistent with evaluation from the Cassini radar observations.     

Latitudinal variations in Titan’s methane and haze from Cassini VIMS observations

We analyze observations taken with Cassini’s Visual and Infrared Mapping Spectrometer (VIMS), to determine the current methane and haze latitudinal distribution between 60°S and 40°N. The methane variation was measured primarily from its absorption band at 0.61 μm, which is optically thin enough to be sensitive to the methane abundance at 20-50 km altitude. Haze characteristics were determined from Titan’s 0.4-1.6 μm spectra, which sample Titan’s atmosphere from the surface to 200 km altitude. Radiative transfer models based on the haze properties and methane absorption profiles at the Huygens site reproduced the observed VIMS spectra and allowed us to retrieve latitude variations in the methane abundance and haze. We find the haze variations can be reproduced by varying only the density and single scattering albedo above 80 km altitude. There is an ambiguity between methane abundance and haze optical depth, because higher haze optical depth causes shallower methane bands; thus a family of solutions is allowed by the data. We find that haze variations alone, with a constant methane abundance, can reproduce the spatial variation in the methane bands if the haze density increases by 60% between 20°S and 10°S (roughly the sub-solar latitude) and single scattering absorption increases by 20% between 60°S and 40°N. On the other hand, a higher abundance of methane between 20 and 50 km in the summer hemisphere, as much as two times that of the winter hemisphere, is also possible, if the haze variations are minimized. The range of possible methane variations between 27°S and 19°N is consistent with condensation as a result of temperature variations of 0-1.5 K at 20-30 km. Our analysis indicates that the latitudinal variations in Titan’s visible to near-IR albedo, the north/south asymmetry (NSA), result primarily from variations in the thickness of the darker haze layer, detected by Huygens DISR, above 80 km altitude. If we assume little to no latitudinal methane variations we can reproduce the NSA wavelength signatures with the derived haze characteristics. We calculate the solar heating rate as a function of latitude and derive variations of ∼10-15% near the sub-solar latitude resulting from the NSA. Most of the latitudinal variations in the heating rate stem from changes in solar zenith angle rather than compositional variations.     

Disk-integrated bolometric Bond albedos and rotational light curves of saturnian satellites from Cassini Visual and Infrared Mapping Spectrometer

We present values from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) of four fundamental disk-integrated spectrophotometric properties (bolometric Bond albedo, solar phase curve, phase integral, and geometric albedo at 7-15 different wavelengths in the λ = 0.35-5.1 μm range) for five mid-sized saturnian icy satellites: Rhea, Dione, Tethys, Mimas, and Enceladus. These values, which include data from the period 2004-2008 and add to past VIMS phase curves, include opposition surge effects at down to fractions of a degree in solar phase angle for several moons and extend to over double the solar phase angle coverage of the Voyager mission. We also present new rotational light curves for Rhea and Dione at 7 near-infrared bands not previously available in ground-based or spacecraft studies. The bolometric Bond albedos we derive are as follows: 0.48 ± 0.09 (Rhea), 0.52 ± 0.08 (Dione), 0.61 ± 0.09 (Tethys), 0.67 ± 0.10 (Mimas), and 0.85 ± 0.11 (Enceladus). We also provide breakdowns of the major photometric quantities in both leading and trailing hemispheres. These refined parameters can be used to construct future bolometric Bond albedo maps that will contribute to surface composition identification studies, as well as models of volatile transport and sublimation. Through such applications, these data will help to determine the physical properties of surface particles, how the E-ring affects the inner saturnian moons, what is responsible for the dark albedo patterns seen on Tethys, and if these moons (e.g., Dione) are geologically active.     

Titan impacts and escape

We report on hydrodynamic calculations of impacts of large (multi-kilometer) objects on Saturn’s moon Titan. We assess escape from Titan, and evaluate the hypothesis that escaping ejecta blackened the leading hemisphere of Iapetus and peppered the surface of Hyperion.We carried out two- and three-dimensional simulations of impactors ranging in size from 4 to 100 km diameter, impact velocities between 7 and 15 km s⁻¹, and impact angles from 0° to 75° from the vertical. We used the ZEUSMP2 hydrocode for the calculations. Simulations were made using three different geometries: three-dimensional Cartesian, two-dimensional axisymmetric spherical polar, and two-dimensional plane polar. Three-dimensional Cartesian geometry calculations were carried out over a limited domain (e.g. 240 km on a side for an impactor of size d_i = 10 km), and the results compared to ones with the same parameters done by Artemieva and Lunine (2005); in general the comparison was good. Being computationally less demanding, two-dimensional calculations were possible for much larger domains, covering global regions of the satellite (from 800 km below Titan’s surface to the exobase altitude 1700 km above the surface). Axisymmetric spherical polar calculations were carried out for vertical impacts. Two-dimensional plane-polar geometry calculations were made for both vertical and oblique impacts. In general, calculations among all three geometries gave consistent results.Our basic result is that the amount of escaping material is less than or approximately equal to the impactor mass even for the most favorable cases. Amounts of escaping material scaled most strongly as a function of velocity, with high-velocity impacts generating the largest amount, as expected. Dependence of the relative amount of escaping mass f_esc = m_esc/M_i on impactor diameter d_i was weak. Oblique impacts (impact angle θ_i > 45°) were more effective than vertical or near-vertical impacts; ratios of m_esc/M_i ∼ 1-2 were found in the simulations.     

Europa’s ridged plains and smooth low albedo plains: Distinctive compositions and compositional gradients at the leading side-trailing side boundary

This investigation uses linear mixture modeling employing cryogenic laboratory reference spectra to estimate surface compositions and water ice grain sizes of Europa’s ridged plains and smooth low albedo plains. Near-infrared spectra for 23 exposures of ridged plains materials are analyzed along with 11 spectra representing low albedo plains. Modeling indicates that these geologic units differ both in the relative abundance of non-ice hydrated species and in the abundance and grain sizes of water ice. The background ridged plains in our study area appear to consist predominantly of water ice (∼46%) with approximately equal amounts (on average) of hydrated sulfuric acid (∼27%) and hydrated salts (∼27%). The solutions for the smooth low albedo plains are dominated by hydrated salts (∼62%), with a relatively low mean abundance of water ice (∼10%), and an abundance of hydrated sulfuric acid similar to that found in ridged plains (∼27%). The model yields larger water ice grain sizes (100 μm versus 50-75 μm) in the ridged plains. The 1.5-μm water ice absorption band minimum is found at shorter wavelengths in the low albedo plains deposits than in the ridged plains (1.498 ± .003 μm versus 1.504 ± .001 μm). The 2.0-μm band minimum in the low albedo plains exhibits a somewhat larger blueshift (1.964 ± .006 μm versus 1.983 ± .006 μm for the ridged plains).The study area spans longitudes from 168° to 185°W, which includes Europa’s leading side-trailing side boundary. A well-defined spatial gradient of sulfuric acid hydrate abundance is found for both geologic units, with concentrations increasing in the direction of the trailing side apex. We associate this distribution with the exogenic effects of magnetospheric charged particle bombardment and associated chemical processing of surface materials (the radiolytic sulfur cycle). However, one family of low albedo plains exposures exhibits sulfuric acid hydrate abundances up to 33% lower than found for adjacent exposures, suggesting that these materials have undergone less processing, thus implying that these deposits may have been emplaced more recently.Modeling identifies high abundances (to 30%) of magnesium sulfate brines in the low albedo plains exposures. Our investigation marks the first spectroscopic identification of MgSO₄ brine on Europa. We also find significantly higher abundances of sodium-bearing species (bloedite and mirabilite) in the low albedo plains. The results illuminate the role of radiolytic processes in modifying the surface composition of Europa, and may provide new constraints for models of the composition of Europa’s putative subsurface ocean.     

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.     

New Cassini RADAR results for Saturn’s icy satellites

Cassini radar tracks on Saturn’s icy satellites through the end of the Prime Mission in 2008 have increased the number of radar albedo estimates from 10 (Ostro et al., 2006) to 73. The measurements sample diverse subradar locations (and for Dione, Rhea, and Iapetus almost always use beamwidths less than half the target angular diameters), thereby constraining the satellites’ global radar albedo distributions. The echoes result predominantly from volume scattering, and their strength is thus strongly sensitive to ice purity and regolith maturity. The combination of the Cassini data set and Arecibo 13-cm observations of Enceladus, Tethys, Dione, Rhea (Black et al., 2007), and Iapetus (Black et al., 2004) discloses an unexpectedly complex pattern of 13-to-2-cm wavelength dependence. The 13-cm albedos are generally smaller than 2-cm albedos and lack the correlation seen between 2-cm and optical geometric albedos. Enceladus and Iapetus are the most interesting cases. We infer from hemispheric albedo variations that the E-ring has a prominent effect on the 13-cm radar “lightcurve”. The uppermost trailing-side regolith is too fresh for meteoroid bombardment to have developed larger-scale heterogeneities that would be necessary to elevate the 13-cm radar albedo, whereas all of Enceladus is clean and mature enough for the 2-cm albedo to be uniformly high. For, Iapetus, the 2-cm albedo is strongly correlated with optical albedo: low for the optically dark, leading-side material and high for the optically bright, trailing-side material. However, Iapetus’ 13-cm albedo values show no significant albedo dichotomy and are several times lower than 2-cm values, being indistinguishable from the weighted mean of 13-cm albedos for main-belt asteroids, 0.15 ± 0.10. The leading side’s optically dark contaminant must be present to depths of at least one to several decimeters, so 2-cm albedos can mimic the optical dichotomy; however, it does not have to extend any deeper than that. The fact that both hemispheres of Iapetus look Asteroid-like at 13 cm means that coherent backscattering itself is not nearly as effective as it is at 2 cm. Since Iapetus’ entire surface is mature regolith, the wavelength dependence must involve composition, not structure. Either the composition is a function of depth everywhere (with electrical loss much greater at depths greater than a decimeter or two), or the intrinsic electrical loss of some pervasive constituent is much higher at 13 cm than at 2 cm. Ammonia is a candidate for such a contaminant. If ammonia’s electrical properties do not depend on frequency, and if ammonia is globally much less abundant within the upper one or two decimeters than at greater depths, then coherent backscattering would effectively be shut down at 13 cm, explaining the Asteroid-like 13-cm albedo.     

Geology of the Selk crater region on Titan from Cassini VIMS observations

Observations of Titan obtained by the Cassini Visual and Infrared Mapping Spectrometer (VIMS) have revealed Selk crater, a geologically young, bright-rimmed, impact crater located ∼800 km north-northwest of the Huygens landing site. The crater rim-crest diameter is ∼90 km; its floor diameter is ∼60 km. A central pit/peak, 20-30 km in diameter, is seen; the ratio of the size of this feature to the crater diameter is consistent with similarly sized craters on Ganymede and Callisto, all of which are dome craters. The VIMS data, unfortunately, are not of sufficient resolution to detect such a dome. The inner rim of Selk crater is fluted, probably by eolian erosion, while the outer flank and presumed ejecta blanket appear dissected by drainages (particularly to the east), likely the result of fluvial erosion. Terracing is observed on the northern and western walls of Selk crater within a 10-15 km wide terrace zone identified in VIMS data; the terrace zone is bright in SAR data, consistent with it being a rough surface. The terrace zone is slightly wider than those observed on Ganymede and Callisto and may reflect differences in thermal structure and/or composition of the lithosphere. The polygonal appearance of the crater likely results from two preexisting planes of weakness (oriented at azimuths of 21° and 122° east of north). A unit of generally bright terrain that exhibits similar infrared-color variation and contrast to Selk crater extends east-southeast from the crater several hundred kilometers. We informally refer to this terrain as the Selk “bench.” Both Selk and the bench are surrounded by the infrared-dark Belet dune field. Hypotheses for the genesis of the optically bright terrain of the bench include: wind shadowing in the lee of Selk crater preventing the encroachment of dunes, impact-induced cryovolcanism, flow of a fluidized-ejecta blanket (similar to the bright crater outflows observed on Venus), and erosion of a streamlined upland formed in the lee of Selk crater by fluid flow. Vestigial circular outlines in this feature just east of Selk’s ejecta blanket suggest that this might be a remnant of an ancient, cratered crust. Evidently the southern margin of the feature has sufficient relief to prevent the encroachment of dunes from the Belet dune field. We conclude that this feature either represents a relatively high-viscosity, fluidized-ejecta flow (a class intermediate to ejecta blankets and long venusian-style ejecta flows) or a streamlined upland remnant that formed downstream from the crater by erosive fluid flow from the west-northwest.     

Impact-induced N2 production from ammonium sulfate: Implications for the origin and evolution of N2 in Titan’s atmosphere

Chemical reactions and volatile supply through hypervelocity impacts may have played a key role for the origin and evolution of both planetary and satellite atmospheres. In this study, we evaluate the role of impact-induced N₂ production from reduced nitrogen-bearing solids proposed to be contained in Titan’s crust, ammonium sulfate ((NH₄)₂SO₄), for the replenishment of N₂ to the atmosphere in Titan’s history. To investigate the conversion of (NH₄)₂SO₄ into N₂ by hypervelocity impacts, we measured gases released from (NH₄)₂SO₄ that was exposed to hypervelocity impacts created by a laser gun. The sensitivity and accuracy of the measurements were enhanced by using an isotope labeling technique for the target. We obtained the efficiency of N₂ production from (NH₄)₂SO₄ as a function of peak shock pressure ranging from ∼8 to ∼45 GPa. Our results indicate that the initial and complete shock pressures for N₂ degassing from (NH₄)₂SO₄ are ∼10 and ∼25 GPa, respectively. These results suggest that cometary impacts on Titan (i.e., impact velocity v_i > ∼8 km/s) produce N₂ efficiently; whereas satellitesimal impacts during the accretion (i.e., v_i < 4 km/s) produce N₂ only inefficiently. Even when using the proposed small amount of (NH₄)₂SO₄ content in the crust (∼4 wt.%) (Fortes, A.D. et al., 2007. Icarus 188, 139-153), the total amount of N₂ provided through cometary impacts over 4.5 Ga reaches ∼2-6 times the present atmospheric N₂ (i.e., ∼7 × 10²⁰-2 × 10²¹ [mol]) based on the measured production efficiency and results of a hydrodynamic simulation of cometary impacts onto Titan. This implies that cometary impacts onto Titan’s crust have the potential to account for a large part of the present N₂ through the atmospheric replenishment after the accretion.     

Channel morphometry, sediment transport, and implications for tectonic activity and surficial ages of Titan basins

Fluvial features on Titan and drainage basins on Earth are remarkably similar despite differences in gravity and surface composition. We determined network bifurcation (R_b) ratios for five Titan and three terrestrial analog basins. Tectonically-modified Earth basins have R_b values greater than the expected range (3.0-5.0) for dendritic networks; comparisons with R_b values determined for Titan basins, in conjunction with similarities in network patterns, suggest that portions of Titan’s north polar region are modified by tectonic forces. Sufficient elevation data existed to calculate bed slope and potential fluvial sediment transport rates in at least one Titan basin, indicating that 75 mm water ice grains (observed at the Huygens landing site) should be readily entrained given sufficient flow depths of liquid hydrocarbons. Volumetric sediment transport estimates suggest that ∼6700-10,000 Titan years (∼2.0-3.0 × 10⁵ Earth years) are required to erode this basin to its minimum relief (assuming constant 1 m and 1.5 m flows); these lowering rates increase to ∼27,000-41,000 Titan years (∼8.0-12.0 × 10⁵ Earth years) when flows in the north polar region are restricted to summer months.     

Finding the trigger to Iapetus’ odd global albedo pattern: Dynamics of dust from Saturn’s irregular satellites

The leading face of Saturn’s moon Iapetus, Cassini Regio, has an albedo only one tenth that on its trailing side. The origin of this enigmatic dichotomy has been debated for over 40 years, but with new data, a clearer picture is emerging. Motivated by Cassini radar and imaging observations, we investigate Soter’s model of dark exogenous dust striking an originally brighter Iapetus by modeling the dynamics of the dark dust from the ring of the exterior retrograde satellite Phoebe under the relevant perturbations. In particular, we study the particles’ probabilities of striking Iapetus, as well as their expected spatial distribution on the Iapetian surface. We find that, of the long-lived particles (?5 μm), most particle sizes (?10 μm) are virtually certain to strike Iapetus, and their calculated distribution on the surface matches up well with Cassini Regio’s extent in its longitudinal span. The satellite’s polar regions are observed to be bright, presumably because ice is deposited there. Thus, in the latitudinal direction we estimate polar dust deposition rates to help constrain models of thermal migration invoked to explain the bright poles (Spencer, J.R., Denk, T. [2010]. Science 327, 432-435). We also analyze dust originating from other irregular outer moons, determining that a significant fraction of that material will eventually coat Iapetus—perhaps explaining why the spectrum of Iapetus’ dark material differs somewhat from that of Phoebe. Finally we track the dust particles that do not strike Iapetus, and find that most land on Titan, with a smaller fraction hitting Hyperion. As has been previously conjectured, such exogenous dust, coupled with Hyperion’s chaotic rotation, could produce Hyperion’s roughly isotropic, moderate-albedo surface.