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.     

Solar phase curves and phase integrals for the leading and trailing hemispheres of Iapetus from the Cassini Visual Infrared Mapping Spectrometer

We performed photometry of Cassini Visual Infrared Mapping Spectrometer observations of Iapetus to produce the first phase integrals calculated directly from solar phase curves of Iapetus for the leading hemisphere and to estimate the phase integrals for the trailing hemisphere. We also explored the phase integral dependence on wavelength and geometric albedo. The extreme dichotomy of the brightness of the leading and trailing sides of Iapetus is reflected in their phase integrals. Our phase integrals, which are lower than the results of Morrison et al. (Morrison, D., Jones, T.J., Cruikshank, D.P., Murphy, R.E. [1975]. Icarus 24, 157-171) and Squyres et al. (Squyres, S.W., Buratti, B.J., Veverka, J., Sagan, C. [1984]. Icarus 59, 426-435), have profound implications on the energy balance and volatile transport on this icy satellite.     

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.     

Cassini RADAR observations of Enceladus, Tethys, Dione, Rhea, Iapetus, Hyperion, and Phoebe

Cassini 2.2-cm radar and radiometric observations of seven of Saturn's icy satellites yield properties that apparently are dominated by subsurface volume scattering and are similar to those of the icy Galilean satellites. Average radar albedos decrease in the order Enceladus/Tethys, Hyperion, Rhea, Dione, Iapetus, and Phoebe. This sequence most likely corresponds to increasing contamination of near-surface water ice, which is intrinsically very transparent at radio wavelengths. Plausible candidates for contaminants include ammonia, silicates, metallic oxides, and polar organics (ranging from nitriles like HCN to complex tholins). There is correlation of our targets' radar and optical albedos, probably due to variations in the concentration of optically dark contaminants in near-surface water ice and the resulting variable attenuation of the high-order multiple scattering responsible for high radar albedos. Our highest radar albedos, for Enceladus and Tethys, probably require that at least the uppermost one to several decimeters of the surface be extremely clean water ice regolith that is structurally complex (i.e., mature) enough for there to be high-order multiple scattering within it. At the other extreme, Phoebe has an asteroidal radar reflectivity that may be due to a combination of single and volume scattering. Iapetus' 2.2-cm radar albedo is dramatically higher on the optically bright trailing side than the optically dark leading side, whereas 13-cm results reported by Black et al. [Black, G.J., Campbell, D.B., Carter, L.M., Ostro, S.J., 2004. Science 304, 553] show hardly any hemispheric asymmetry and give a mean radar reflectivity several times lower than the reflectivity measured at 2.2 cm. These Iapetus results are understandable if ammonia is much less abundant on both sides within the upper one to several decimeters than at greater depths, and if the leading side's optically dark contaminant is present to depths of at least one to several decimeters. As argued by Lanzerotti et al. [Lanzerotti, L.J., Brown, W.L., Marcantonio, K.J., Johnson, R.E., 1984. Nature 312, 139-140], a combination of ion erosion and micrometeoroid gardening may have depleted ammonia from the surfaces of Saturn's icy satellites. Given the hypersensitivity of water ice's absorption length to ammonia concentration, an increase in ammonia with depth could allow efficient 2.2-cm scattering from within the top one to several decimeters while attenuating 13-cm echoes, which would require a six-fold thicker scattering layer. If so, we would expect each of the icy satellites' average radar albedos to be higher at 2.2 cm than at 13 cm, as is the case so far with Rhea [Black, G., Campbell, D., 2004. Bull. Am. Astron. Soc. 36, 1123] as well as Iapetus.     

The variability of volcanic activity at Zamama, Culann, and Tupan Patera on Io as seen by the Galileo Near Infrared Mapping Spectrometer

Zamama, Culann, and Tupan Patera are three large, persistent volcanic centers on the jovian moon Io. As part of an ongoing project to quantify contributions from individual volcanic centers to Io’s thermal budget, we have quantified the radiant flux from all suitable observations made by the Galileo Near Infrared Mapping Spectrometer (NIMS) of these volcanoes, in some cases filling omissions in previous analyses. At Zamama, after a long period of cooling, we see a peak in thermal emission that corresponds with new plume activity. Subsequently, toward the end of the Galileo epoch, thermal emission from Zamama drops off in a manner consistent with a greatly reduced eruption rate and the cooling of emplaced flows. Culann exhibits possible episodic activity. We present the full Tupan Patera NIMS dataset and derive new estimates of thermal output and temporal behavior. Eruption rates at these three volcanoes are on the order of 30 m³ s⁻¹, consistent with a previous analysis of NIMS observations of Prometheus, and nearly an order of magnitude greater than Kilauea volcano, Hawai’i, Earth’s most active volcano. We propose that future missions to the jovian system could better constrain activity at these volcanoes and others where similar styles of activity are taking place by obtaining data on a time scale of, ideally, at least one observation per day. Observations at similar or even shorter timescales are desirable during initial waxing phases of eruption episodes. These eruptions are identifiable from their characteristic spectral signatures and temporal behavior.     

Tidal evolution of Mimas, Enceladus, and Dione

The tidal evolution through several resonances involving Mimas, Enceladus, and/or Dione is studied numerically with an averaged resonance model. We find that, in the Enceladus-Dione 2:1 e-Enceladus type resonance, Enceladus evolves chaotically in the future for some values of k₂/Q. Past evolution of the system is marked by temporary capture into the Enceladus-Dione 4:2 ee^′-mixed resonance. We find that the free libration of the Enceladus-Dione 2:1 e-Enceladus resonance angle of 1.5° can be explained by a recent passage of the system through a secondary resonance. In simulations with passage through the secondary resonance, the system enters the current Enceladus-Dione resonance close to tidal equilibrium and thus the equilibrium value of tidal heating of 1.1(18,000/QS) GW applies. We find that the current anomalously large eccentricity of Mimas can be explained by passage through several past resonances. In all cases, escape from the resonance occurs by unstable growth of the libration angle, sometimes with the help of a secondary resonance. Explanation of the current eccentricity of Mimas by evolution through these resonances implies that the Q of Saturn is below 100,000. Though the eccentricity of Enceladus can be excited to moderate values by capture in the Mimas-Enceladus 3:2 e-Enceladus resonance, the libration amplitude damps and the system does not escape. Thus past occupancy of this resonance and consequent tidal heating of Enceladus is excluded. The construction of a coherent history places constraints on the allowed values of k₂/Q for the satellites.     

Water ice crystallinity and grain sizes on Dione

Saturn’s satellite Dione is becoming an increasingly important object in the outer Solar System, as evidence for its current activity accumulates. Infrared observations of the surface can provide clues to the history of the body and currently active processes. Using data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS), we perform three sets of analyses that are sensitive to the ice state, temperature, thermal history, grain size and composition of surface ice. These are calculation of a “crystallinity factor”, spectral ratios and water ice band depths. In our analysis, we focus on the dichotomy between the wispy and dark terrain on Dione’s trailing hemisphere, to better understand the source of the different materials and their current properties. Our results suggest two different scenarios: (1) the ice from the wispy region has a higher crystallinity and water ice content than the dark region or (2) the wispy region contains larger grains. Both of these models imply recent geologic activity on Dione.     

Compositional mapping of Saturn's satellite Dione with Cassini VIMS and implications of dark material in the Saturn system

Cassini VIMS has obtained spatially resolved imaging spectroscopy data on numerous satellites of Saturn. A very close fly-by of Dione provided key information for solving the riddle of the origin of the dark material in the Saturn system. The Dione VIMS data show a pattern of bombardment of fine, sub-0.5-μm diameter particles impacting the satellite from the trailing side direction. Multiple lines of evidence point to an external origin for the dark material on Dione, including the global spatial pattern of dark material, local patterns including crater and cliff walls shielding implantation on slopes facing away from the trailing side, exposing clean ice, and slopes facing the trailing direction which show higher abundances of dark material. Multiple spectral features of the dark material match those seen on Phoebe, Iapetus, Hyperion, Epimetheus and the F-ring, implying the material has a common composition throughout the Saturn system. However, the exact composition of the dark material remains a mystery, except that bound water and, tentatively, ammonia are detected, and there is evidence both for and against cyanide compounds. Exact identification of composition requires additional laboratory work. A blue scattering peak with a strong UV-visible absorption is observed in spectra of all satellites which contain dark material, and the cause is Rayleigh scattering, again pointing to a common origin. The Rayleigh scattering effect is confirmed with laboratory experiments using ice and 0.2-μm diameter carbon grains when the carbon abundance is less than about 2% by weight. Rayleigh scattering in solids is also confirmed in naturally occurring terrestrial rocks, and in previously published reflectance studies. The spatial pattern, Rayleigh scattering effect, and spectral properties argue that the dark material is only a thin coating on Dione's surface, and by extension is only a thin coating on Phoebe, Hyperion, and Iapetus, although the dark material abundance appears higher on Iapetus, and may be locally thick. As previously concluded for Phoebe, the dark material appears to be external to the Saturn system and may be cometary in origin. We also report a possible detection of material around Dione which may indicate Dione is active and contributes material to the E-ring, but this observation must be confirmed.     

Thermal inertia and bolometric Bond albedo values for Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus as derived from Cassini/CIRS measurements

Spectra taken by Cassini’s Composite Infrared Spectrometer (CIRS) between 10 and 600 cm⁻¹ (17-1000 μm) of surface thermal emission of Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus have been used to derive the thermal inertia and bolometric Bond albedo values. Only an upper limit for the bolometric Bond albedo of Iapetus’ dark leading side could be determined due to the insensitivity of the thermal model to albedo when albedos are very low. The thermal inertia in this region however is better constrained. The CIRS coverage of Enceladus is extensive enough that the latitudinal variation in these values from 60°S to 70°N has been determined in 10° wide bins. The bolometric Bond albedos determined here are consistent with literature values which show the surface of the saturnian icy moons to be covered in ice contaminated to varying degrees. The thermal inertia of the moons is shown to be in the range 9-, approximately 2-6 times lower than that of the Galilean satellites, implying a less well consolidated and more porous surface. The thermal inertias of Iapetus and Phoebe are somewhat higher, suggesting that the very low thermal inertias of satellites from Rhea inwards may be related to their probable coating of E-ring material. Latitudinal variations on the surface of Enceladus show that the bolometric Bond albedo and thermal inertia increase towards the active plume source at the south pole.     

Titan's surface: Search for spectral diversity and composition using the Cassini VIMS investigation

The surface composition of Titan is of great importance for understanding both the internal evolution of Titan and its atmosphere. The Visual and Infrared Mapping Spectrometer (VIMS) investigation on Cassini is observing Titan from 0.35 to 5.11 μm with spatial resolution down to a few kilometers during each flyby of the spacecraft as it orbits Saturn. Our search for spectral diversity using seven methane transmission windows in the near infrared suggests that spectrally distinct units exist on the surface of Titan and that most of the surface can be modeled using only a few distinct spectral units: water frost, CO₂ frost, atmospheric scattering, and an unknown material bright at 2 μm. A dark, spectrally neutral material is also implied. Use of an atmospheric scattering component with spectral mixing analysis may provide a method for partially removing atmospheric effects. In some locations, atmospheric scattering accounts for the majority of the signal. There are also small regions with unusual spectra that may be due to low signal and high noise and/or may be exotic materials of interest. Further, we searched within the methane windows for spectral features associated with Titan's surface. Only the 5-μm and, to a lesser extent, the 2-μm window provide a reasonable opportunity for this, as the shorter-wavelength windows are too narrow and the 2.8-μm window is cluttered with an unknown atmospheric constituent. We find evidence for only one spectral feature: near 4.92 μm for the 5-μm bright Tui Regio region. CO₂ frost with grains smaller than about 10 μm is the best candidate we have found so far to explain this absorption as well as the feature's spectral contrast between the 2.7- and the 2.8-μm atmosphere subwindows. This suggested CO₂ identification is supported by the presence of an endmember in the spectral mixture analysis that is consistent with CO₂ frost with large grain sizes. We find no other absorption features that are statistically significant, including those reported earlier by others. These results are consistent with but greatly extend our early analysis that treated only the Ta data set [McCord, T.B., et al., 2006a. Planet. Space Sci. 54, 1524-1539]. In the spectral feature search process, we explored in detail the noise characteristics of the VIMS data within the 5-μm window, which has generally very low signal (4-20 DN), due to the measurement conditions and low illumination levels. We find noise of nearly Gaussian statistics except for some erratic darks and noise spikes, and the data set seems generally well behaved. We present examples of our attempt to improve on the standard VIMS pipeline data calibration.     

Spectroscopy, morphometry, and photoclinometry of Titan's dunefields from Cassini/VIMS

Fine-resolution (500 m/pixel) Cassini Visual and Infrared Mapping Spectrometer (VIMS) T20 observations of Titan resolve that moon's sand dunes. The spectral variability in some dune regions shows that there are sand-free interdune areas, wherein VIMS spectra reveal the exposed dune substrate. The interdunes from T20 are, variously, materials that correspond to the equatorial bright, 5-μm-bright, and dark blue spectral units. Our observations show that an enigmatic “dark red” spectral unit seen in T5 in fact represents a macroscopic mixture with 5-μm-bright material and dunes as its spectral endmembers. Looking more broadly, similar mixtures of varying amounts of dune and interdune units of varying composition can explain the spectral and albedo variability within the dark brown dune global spectral unit that is associated with dunes. The presence of interdunes indicates that Titan's dunefields are both mature and recently active. The spectrum of the dune endmember reveals the sand to be composed of less water ice than the rest of Titan; various organics are consistent with the dunes' measured reflectivity. We measure a mean dune spacing of 2.1 km, and find that the dunes are oriented on the average in an east-west direction, but angling up to 10° from parallel to the equator in specific cases. Where no interdunes are present, we determine the height of one set of dunes photoclinometrically to be between 30 and 70 m. These results pave the way for future exploration and interpretation of Titan's sand dunes.     

Sizes, shapes, and derived properties of the saturnian satellites after the Cassini nominal mission

The Cassini spacecraft has provided data on over twenty satellites orbiting Saturn. In this study we update the shape measurements derived from imaging data of nearly all the observed regular satellites (those with low eccentricities and inclinations) and briefly discuss some of the implications of their properties. In particular, the improved data show that Rhea’s shape is hydrostatic.     

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.     

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.     

Observations with the Visual and Infrared Mapping Spectrometer (VIMS) during Cassini's flyby of Jupiter

The Cassini Visual and Infrared Mapping Spectrometer (VIMS) is an imaging spectrometer covering the wavelength range 0.3-5.2 μm in 352 spectral channels, with a nominal instantaneous field of view of 0.5 mrad. The Cassini flyby of Jupiter represented a unique opportunity to accomplish two important goals: scientific observations of the jovian system and functional tests of the VIMS instrument under conditions similar to those expected to obtain during Cassini's 4-year tour of the saturnian system. Results acquired over a complete range of visual to near-infrared wavelengths from 0.3 to 5.2 μm are presented. First detections include methane fluorescence on Jupiter, a surprisingly high opposition surge on Europa, the first visual-near-IR spectra of Himalia and Jupiter's optically-thin ring system, and the first near-infrared observations of the rings over an extensive range of phase angles (0-120°). Similarities in the center-to-limb profiles of H⁺₃ and CH₄ emissions indicate that the H⁺₃ ionospheric density is solar-controlled outside of the auroral regions. The existence of jovian NH₃ absorption at 0.93 μm is confirmed. Himalia has a slightly reddish spectrum, an apparent absorption near 3 μm, and a geometric albedo of 0.06±0.01 at 2.2 μm (assuming an 85-km radius). If the 3-μm feature in Himalia's spectrum is eventually confirmed, it would be suggestive of the presence of water in some form, either free, bound, or incorporated in layer-lattice silicates. Finally, a mean ring-particle radius of 10 μm is found to be consistent with Mie-scattering models fit to VIMS near-infrared observations acquired over 0-120° phase angle.     

A hydrocode calculation coupled with reaction kinetics of carbon compounds within an impact vapor plume and its implications for cometary impacts on Galilean satellites

The synthesis of organic molecules via chemical reactions within impact vapor plumes has been proposed as a mechanism to supply organics on a planet. However, the kinetics of chemical reactions within a rapidly expanding vapor plume or quenching process of the reactions has not been studied extensively. In this study, we constructed a new numerical model that calculates kinetics of the entire chemical reactions within an impact vapor plume. Numerical results revealed that the semi-analytical models proposed so far, in which the final amount of a chemical species was given by the equilibrium abundance at the quenching temperature of the fastest reaction path involving the species, underestimates the yield of organic molecules, such as HCN, by up to a factor of 10. This is because the previously used assumption that a species can achieve equilibrium with the rest of the reaction system via the fastest reaction path involving the species is not necessarily valid. Our analysis of the high-temperature H/C/N/O reaction system suggests that the quenching of slow reactions divides the reaction network into smaller reaction sub-systems isolated from the rest of the reaction system. Then, the fastest reaction path cannot equilibrate an isolated reaction sub-system with the rest of the reaction system. Simulation of this actual disequilibrium mechanism requires a simultaneous numerical calculation of the entire reaction network, which is equivalent to conducting a full kinetic model calculation, such as our model. Our numerical code makes it possible to discuss quantitatively the impact chemistry for various situations, such as the Galilean satellites. In this study, our numerical model is applied to the delivery of organic molecules via cometary impact on the Galilean satellites. Our numerical results indicate that small-particle impacts would produce HCN efficiently. Resulting HCN may freeze out immediately and be deposited on satellite surfaces, where it may be eventually converted into complex organics via irradiation of charged particle. On the other hand, large-size impacts may form transient CH₄-N₂ atmospheres, in which complex organics (tholin) may be formed via energy deposition of UV and/or charged particle. Resulting complex organics may subsequently precipitate on the satellite surfaces without clear correlation with the locations of impact craters. Such distribution of complex organics created by chemical reactions within vapor plumes due to cometary impacts may explain an absorption (4.57 μm) on Galilean satellites nonassociated with observable (moderate- and large-size) impact craters.     

Experimental study on the rheological properties of polycrystalline solid nitrogen and methane: Implications for tectonic processes on Triton

To investigate the evolution of any processes on planetary surfaces in the outer Solar System, the rheological properties of non-water ices were studied by means of a sound velocity measurement system and a uniaxial deformation apparatus. A pulse transmission method was used to obtain longitudinal (V_p) and transverse (V_s) wave velocities through solid nitrogen and methane at temperatures ranging from 5 to 64 K and from 5 to 90 K, respectively. The measured velocities confirmed that the solid methane and solid nitrogen samples were non-porous polycrystalline samples without any cracks and bubbles inside. Compression tests at constant strain-rate were performed for solid nitrogen and methane at temperatures of 5-56 K and 5-77 K, respectively, at strain-rates of 10⁻⁴-10⁻² s⁻¹. Both brittle and ductile behavior was observed for solid nitrogen and methane under these conditions. The maximum strength of solid nitrogen was observed to be 9 MPa in the brittle failure mode, and that of solid methane was 10 MPa. These low strengths cannot support cantaloupe structures with the topographic undulation larger than several kilometers found on Triton’s surface, suggesting that other materials such as H₂O ice could underlay solid methane and nitrogen and support these structures.     

Morphology, temperature, and eruption dynamics at Pele

The Pele region of Io has been the site of vigorous volcanic activity from the time of the first Voyager I observations in 1979 up through the final Galileo ones in 2001. There is high-temperature thermal emission from a visibly dark area that is thought to be a rapidly overturning lava lake, and is also the source of a large sulfur-rich plume. We present a new analysis of Voyager I visible wavelength images, and Galileo Solid State Imager (SSI) and Near Infrared Mapping Spectrometer (NIMS) thermal emission observations which better define the morphology of the region and the intensity of the emission. The observations show remarkable correlations between the locations of the emission and the features seen in the Voyager images, which provide insight into eruption mechanisms and constrain the longevity of the activity. We also analyze an additional wavelength channel of NIMS data (1.87 μm) which paradoxically, because of reduced sensitivity, allows us to estimate temperatures at the peak locations of emission. Measurements of eruption temperatures on Io are crucial because they provide our best clues to the composition of the magma. High color temperatures indicative of ultramafic composition have been reported for the Pillan hot spot and possibly for Pele, although recent work has called into question the requirement for magma temperatures above those expected for ordinary basalts. Our new analysis of the Pele emission near the peak of the hot spot shows color temperatures near the upper end of the basalt range during the I27 and I32 encounters. In order to analyze the observed color temperatures we also present an analytical model for the thermal emission from fire-fountains, which should prove generally useful for analyzing similar data. This is a modification of the lava flow emission model presented in Howell (Howell, R.R. [1997]. Icarus 127, 394-407), adapted to the fire-fountain cooling curves first discussed in Keszthelyi et al. (Keszthelyi, L., Jaeger, W., Milazzo, M., Radebaugh, J., Davies, A.G., Mitchell, K.L. [2007]. Icarus 192, 491-502). When applied to the I32 observations we obtain a fire-fountain mass eruption rate of 5.1 × 10⁵ kg s⁻¹ for the main vent area and 1.4 × 10⁴ kg s⁻¹ for each of two smaller vent regions to the west. These fire-fountain rates suggest a solution to the puzzling lack of extensive lava flows in the Pele region. Much of the erupted lava may be ejected at high speed into the fire-fountains and plumes, creating dispersed pyroclastic deposits rather than flows. We compare gas and silicate mass eruption rates and discuss briefly the dynamics of this ejection model and the observational evidence.     

The Zamama-Thor region of Io: Insights from a synthesis of mapping, topography, and Galileo spacecraft data

We have studied data from the Galileo spacecraft's three remote sensing instruments (Solid-State Imager (SSI), Near-Infrared Mapping Spectrometer (NIMS), and Photopolarimeter-Radiometer (PPR)) covering the Zamama-Thor region of Io's antijovian hemisphere, and produced a geomorphological map of this region. This is the third of three regional maps we are producing from the Galileo spacecraft data. Our goal is to assess the variety of volcanic and tectonic materials and their interrelationships on Io using planetary mapping techniques, supplemented with all available Galileo remote sensing data. Based on the Galileo data analysis and our mapping, we have determined that the most recent geologic activity in the Zamama-Thor region has been dominated by two sites of large-scale volcanic surface changes. The Zamama Eruptive Center is a site of both explosive and effusive eruptions, which emanate from two relatively steep edifices (Zamama Tholi A and B) that appear to be built by both silicate and sulfur volcanism. A ∼100-km long flow field formed sometime after the 1979 Voyager flybys, which appears to be a site of promethean-style compound flows, flow-front SO₂ plumes, and adjacent sulfur flows. Larger, possibly stealthy, plumes have on at least one occasion during the Galileo mission tapped a source that probably includes S and/or Cl to produce a red pyroclastic deposit from the same vent from which silicate lavas were erupted. The Thor Eruptive Center, which may have been active prior to Voyager, became active again during the Galileo mission between May and August 2001. A pillanian-style eruption at Thor included the tallest plume observed to date on Io (at least 500 km high) and new dark lava flows. The plume produced a central dark pyroclastic deposit (probably silicate-rich) and an outlying white diffuse ring that is SO₂-rich. Mapping shows that several of the new dark lava flows around the plume vent have reoccupied sites of earlier flows. Unlike most of the other pillanian eruptions observed during the Galileo mission, the 2001 Thor eruption did not produce a large red ring deposit, indicating a relative lack of S and/or Cl gases interacting with the magma during that eruption. Between these two eruptive centers are two paterae, Thomagata and Reshef. Thomagata Patera is located on a large shield-like mesa and shows no signs of activity. In contrast, Reshef Patera is located on a large, irregular mesa that is apparently undergoing degradation through erosion (perhaps from SO₂-sapping or chemical decomposition of sulfur-rich material) from multiple secondary volcanic centers.     

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.