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.     

Flow patterns of lobate debris aprons and lineated valley fill north of Ismeniae Fossae, Mars: Evidence for extensive mid-latitude glaciation in the Late Amazonian

A variety of Late Amazonian landforms on Mars have been attributed to the dynamics of ice-related processes. Evidence for large-scale, mid-latitude glacial episodes existing within the last 100 million to 1 billion years on Mars has been presented from analyses of lobate debris aprons (LDA) and lineated valley fill (LVF) in the northern and southern mid-latitudes. We test the glacial hypothesis for LDA and LVF along the dichotomy boundary in the northern mid-latitudes by examining the morphological characteristics of LDA and LVF surrounding two large plateaus, proximal massifs, and the dichotomy boundary escarpment north of Ismeniae Fossae (centered at 45.3°N and 39.2°E). Lineations and flow directions within LDA and LVF were mapped using images from the Context (CTX) camera, the Thermal Emission Imaging Spectrometer (THEMIS), and the High Resolution Stereo Camera (HRSC). Flow directions were then compared to topographic contours derived from the Mars Orbiter Laser Altimeter (MOLA) to determine the down-gradient components of LDA and LVF flow. Observations indicate that flow patterns emerge from numerous alcoves within the plateau walls, are integrated over distances of up to tens of kilometers, and have down-gradient flow directions. Smaller lobes confined within alcoves and superposed on the main LDA and LVF represent a later, less extensive glacial phase. Crater size-frequency distributions of LDA and LVF suggest a minimum (youngest) age of 100 Ma. The presence of ring-mold crater morphologies is suggestive that LDA and LVF are formed of near-surface ice-rich bodies. From these observations, we interpret LDA and LVF within our study region to result from formerly active debris-covered glacial flow, consistent with similar observations in the northern mid-latitudes of Mars. Glacial flow was likely initiated from the accumulation and compaction of snow and ice on plateaus and in alcoves within the plateau walls as volatiles were mobilized to the mid-latitudes during higher obliquity excursions. Together with similar analyses elsewhere along the dichotomy boundary, these observations suggest that multiple glacial episodes occurred in the Late Amazonian and that LDA and LVF represent significant reservoirs of non-polar ice sequestered below a surface lag for hundreds of millions of years.     

A hybrid simulation of Mercury’s magnetosphere for the MESSENGER encounters in year 2008

The latest measurements from the two encounters of the MESSENGER spacecraft in year 2008 have discovered several interesting features of the magnetosphere of Mercury. We have performed high-resolution 3D hybrid model calculations to simulate the solar wind interaction with the Hermean magnetosphere during the first two Mercury encounters of the MESSENGER spacecraft in 2008. It is found that the global structure of the Hermean magnetosphere is significantly controlled by the direction of the interplanetary magnetic field. The bow shock size and shape and the magnetotail configuration have very large differences in these two encounters with northward-pointing and southward-pointing interplanetary magnetic field, respectively. Comparisons are also given with the observed magnetic field profiles and the computational results. In general, good agreement can be found including the interesting feature of the relatively thick magnetopause current layer at outbound measurements. Our work shows that 3D hybrid simulation is a promising method to study in detail the Hermean magnetosphere in parallel with the post-MOI observations of the MESSENGER spacecraft and the Bepi-Colombo mission in future.     

Evaporation effects on the formation of martian gullies

In order to investigate the formation of martian gullies and the stability of fluids on Mars, we examined about 120 gully images. Twelve HiRISE images contained a sufficient number of Transverse Aeolian Ridges (TARs) associated with the gullies to make the following measurements: overall gully length, length of the alcove, channel and apron, and we also measured the frequency of nearby TARs. Six of the 12 images examined showed a statistically significant negative correlation between overall gully length (alcove, channel and apron length) and TAR frequency. Previous experimental work from our group has shown that at temperatures below ∼200 K, evaporation rate increases by about an order of magnitude as wind speed increases from 0 to ∼15 m/s. Thus the negative correlations we observe between gully length and dune frequency can be explained by formation at temperatures below ∼200 K where wind speed/evaporation is a factor governing gully length. In these cases evaporation of the fluid carving the gully was a constraint on their dimensions. Cases where there is no correlation between gully length and TAR frequency, can be explained by formation at temperatures >200 K. The temperatures are consistent with Global Circulation Model and Thermal Emission Spectrometer (TES) data for these latitudes. The temperatures suggested by these trends are consistent with the fluid responsible for gully formation being a strong brine, such as Fe₂(SO₄)₃ which has a eutectic temperature of ∼200 K. We also find that formation timescales for gullies are 10⁵-10⁶ years.     

The analytic properties of the flapping current sheets in the earth magnetotail

The current sheet in Earth’s magnetotail often flaps, and the flapping waves could be induced propagating towards the dawn and dusk flanks, which could make the current sheet dynamic. To explore the dynamic characteristics of current sheet associated with the flapping motion holistically and provide reasonable physical interpretations, detailed direct calculation and analysis have been applied to one approximate analytic model of magnetic field in the flapping current sheet. The main results from the model demonstrate: (1) the magnetic fluctuation amplitude is attenuated from the center of current sheet to the lobe regions; The larger wave amplitude would induce the larger magnetic amplitude; (2) the curvature of magnetic field lines (MFLs), with maximum at the center of current sheet, is only dependent on the displacement Z along the south-north direction from the center of current sheet, regardless of the tilt of current sheet; (3) the half-thickness of neutral sheet, h, the minimum curvature radius of MFLs, R_cmin, and the tilt angle of current sheet, δ, satisfies h=R_cmin cos δ; (4) the gradient of magnetic strength forms a double-peak profile, and the peak value would be more intense if the local current sheet is more tilted; (5) current density j and its j_y, j_z components reach the extremum at the center of CS. j and j_z would be more intense if the local current sheet is more tilted, but it is not the case for j_y; and (6) the field-aligned component of current density mainly appears in the neutral sheet, and the sign of it would change alternatively as the flapping waves passing by. To check the validity of the model, one simulation on the virtual measurements has been made, and the results are in well consistence with actual observations of Cluster.     

The opposition and tilt effects of Saturn’s rings from HST observations

The two major factors contributing to the opposition brightening of Saturn’s rings are (i) the intrinsic brightening of particles due to coherent backscattering and/or shadow hiding on their surfaces, and (ii) the reduced interparticle shadowing when the solar phase angle α → 0°. We utilize the extensive set of Hubble Space Telescope observations (Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 199–223) for different elevation angles B and wavelengths λ to disentangle these contributions. We assume that the intrinsic contribution is independent of B, so that any B dependence of the phase curves is due to interparticle shadowing, which must also act similarly for all λ’s. Our study complements that of Poulet et al. (Poulet, F., Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 224), who used a subset of data for a single B ∼ 10°, and the French et al. (French, R.G., Verbiscer, A., Salo, H., McGhee, C.A., Dones, L. [2007b] PASP 119, 623–642) study for the B ∼ 23° data set that included exact opposition. We construct a grid of dynamical/photometric simulation models, with the method of Salo and Karjalainen (Salo and Karjalainen [2003]. Icarus 164, 428–460), and use these simulations to fit the elevation-dependent part of opposition brightening. Eliminating the modeled interparticle component yields the intrinsic contribution to the opposition effect: for the B and A rings it is almost entirely due to coherent backscattering; for the C ring, an intraparticle shadow hiding contribution may also be present.Based on our simulations, the width of the interparticle shadowing effect is roughly proportional to B. This follows from the observation that as B decreases, the scattering is primarily from the rarefied low filling factor upper ring layers, whereas at larger B’s the dense inner parts are visible. Vertical segregation of particle sizes further enhances this effect. The elevation angle dependence of interparticle shadowing also explains most of the B ring tilt effect (the increase of brightness with elevation). From comparison of the magnitude of the tilt effect at different filters, we show that multiple scattering can account for at most a 10% brightness increase as B → 26°, whereas the remaining 20% brightening is due to a variable degree of interparticle shadowing. The negative tilt effect of the middle A ring is well explained by the the same self-gravity wake models that account for the observed A ring azimuthal brightness asymmetry (Salo, H., Karjalainen, R., French, R.G. [2004]. Icarus 170, 70–90; French, R.G., Salo, H., McGhee, C.A., Dones, L. [2007]. Icarus 189, 493–522).     

Formation of an eroded lava channel within an Elysium Planitia impact crater: Distinguishing between a mechanical and thermal origin

A lava channel identified on the wall of an Elysium Planitia impact crater is investigated to identify the dominant erosion mechanism, mechanical vs. thermal, acting during channel formation. Observations of channel morphology are used to supplement analytical models of lava channel formation in order to calculate the duration of channel formation, the velocity of the lava flowing through the channel, and the erosion rate in each erosion regime considered. Results demonstrate that the channel observed in the Elysium Planitia impact crater formed primarily due to mechanical erosion. In a more general sense, results of this study suggest that lava channels can form primarily due to thermal erosion in the presence of more gradual slopes and more consolidated substrates whereas lava channels can form primarily due to mechanical erosion in the presence of more energetic flows on steeper slopes and more poorly consolidated substrates. Therefore, both erosion regimes must be considered when analyzing origins of eroded lava channels that cut through strata of different strengths.     

Catastrophe of Coronal Magnetic Flux Ropes Caused by Photospheric Motions

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, we carried out numerical simulations to investigate the equilibrium and evolution properties of a magnetic configuration that consists of a coronal magnetic flux rope and a partly open photospheric background field, which is equivalent to that produced by a two-patch magnetic source on the photospheric surface. The axial and annular magnetic fluxes of the flux rope are given and fixed. The global magnetic configuration evolves in response to three types of changes of the background field: decreasing of the distance between the two sources, shrinking of the size of each source, and increasing of the shear in the closed component of the background field. As a result, the geometrical parameters of the flux rope, i.e. the height of the rope axis, the half-width of the rope and the length of the vertical current sheet below the rope, change due to the variation of the background field. It is shown that for a given coronal magnetic flux rope in a partly open background field, the variation of the geometrical parameters of the flux rope displays a catastrophic behavior, namely, there exists a critical point for each case, at which an infinitesimal change of the background field leads to a loss of equilibrium, and thus a jump of the flux rope. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

Low Energy Transfer to the Moon

In 1991, the Japanese Hiten mission used a low energy transfer with a ballistic capture at the Moon which required less Vthan a standard Hohmann transfer. In this paper, we apply the dynamical systems techniques developed in our earlier work to reproduce systematically a Hiten-like mission. We approximate the Sun–Earth–Moon-spacecraft 4-body system as two 3-body systems. Using the invariant manifold structures of the Lagrange points of the 3-body systems, we are able to construct low energy transfer trajectories from the Earth which execute ballistic capture at the Moon. The techniques used in the design and construction of this trajectory may be applied in many situations.     

The Fossil Starburst in M82

We present high-resolution optical and near-infrared HST observations of two adjacent regions in the fossil starburst region inM82, M82 B1 and B2. The presence of both the active and the fossil starburst in M82 provides a unique physical environment to study the stellar and dynamical evolution of star cluster systems. The cluster population in B2 is more heavily affected by internal extinction than that in B1, amounting to an excess extinction in B2 of A_V,excess≃1.1±0.3 mag. Preliminary age estimates date the cluster population in the fossil starburst between ∼2× 10⁸ and ∼ 10⁹ years. The radial luminosity profiles of the brightest clusters are more closely approximated by power laws than by a Gaussian model, in particular in their wings, which favours a slow star formation scenario. This revised version was published online in July 2006 with corrections to the Cover Date.