Ammonia-water ice laboratory studies relevant to outer Solar System surfaces

Although water- and ammonia-ices have been observed or postulated as important components of the icy surfaces of planetary satellites in the outer Solar System, significant gaps exist in our knowledge of the spectra and behavior of such mixtures under astrophysical conditions. To that end, we have completed low-temperature spectroscopic studies (1-20 μm) of water-ammonia mixtures, with an emphasis on features in the near-IR, a region which is accessible to ground-based observations. The influences of composition, formation temperature, thermal- and radiation-processing, and phase (crystalline or amorphous) of the components were examined. Spectra of both pure NH₃ and H₂O-NH₃ icy mixtures with ratios from 0.7 to 57 were measured at temperatures from 10 to 120 K. Conditions for the formation and thermal stability of the ammonia hemihydrate (2NH₃⋅H₂O) and the ammonia monohydrate (NH₃⋅H₂O) have been examined. Band positions of NH₃ in different H₂O-ices and major band positions of the hydrates were measured. We report spectral shifts that depend on concentration and temperature. The radiation-induced amorphization of the hemihydrate was observed and the radiation destruction of NH₃ in H₂O-ices was measured. Implications of these results for the formation, stability, and detection of ammonia on outer satellite surfaces are discussed.     

Understanding the behavior of Prometheus and Pandora

We revisit the dynamics of Prometheus and Pandora, two small moons flanking Saturn's F ring. Departures of their orbits from freely precessing ellipses result from mutual interactions via their 121:118 mean motion resonance. Motions are chaotic because the resonance is split into four overlapping components. Orbital longitudes were observed to drift away from predictions based on Voyager ephemerides. A sudden jump in mean motions took place close to the time at which the orbits' apses were antialigned in 2000. Numerical integrations reproduce both the longitude drifts and the jumps. The latter have been attributed to the greater strength of interactions near apse antialignment (every 6.2 yr), and it has been assumed that this drift-jump behavior will continue indefinitely. We re-examine the dynamics of the Prometheus-Pandora system by analogy with that of a nearly adiabatic, parametric pendulum. In terms of this analogy, the current value of the action of the satellite system is close to its maximum in the chaotic zone. Consequently, at present, the two separatrix crossings per precessional cycle occur close to apse antialignment. In this state libration only occurs when the potential's amplitude is nearly maximal, and the “jumps” in mean motion arise during the short intervals of libration that separate long stretches of circulation. Because chaotic systems explore the entire region of phase space available to them, we expect that at other times the Prometheus-Pandora system would be found in states of medium or low action. In a low action state it would spend most of the time in libration, and separatrix crossings would occur near apse alignment. We predict that transitions between these different states can happen in as little as a decade. Therefore, it is incorrect to assume that sudden changes in the orbits only happen near apse antialignment.     

H2O2 production by high-energy electrons on icy satellites as a function of surface temperature and electron flux

Chemistry on the icy surface of Europa is heavily influenced by the incident energetic particle flux from the jovian magnetosphere. The majority (>75%) of this energy is in the form of high energy electrons (extending to >10 MeV). We have simulated the electron irradiation environment of Europa with a vacuum system containing a high-energy electron gun for irradiation of ice samples formed on a gold mirror cooled with a cryostat. Pure water films of ∼2.6 μm thickness were grown at 100 K and then either cooled (to 80 K), warmed (to 120 K) or left at 100 K and subsequently irradiated with 10 keV electrons. The production of hydrogen peroxide (H₂O₂) was monitored by observation of the 2850 cm⁻¹ (3.5 μm) band. Equilibrium concentrations of H₂O₂, in units of percent by number H₂O₂ relative to water, were found to be 0.043% (80 K), 0.029% (100 K), and 0.0063% (120 K). These values are 33%, 22%, and 5%, respectively, that of the reported surface concentration on the leading hemisphere of Europa (Carlson, R.W., Anderson, M.S., Johnson, R.E., Smythe, W.D., Hendrix, A.R., Barth, C.A., et al. [1999]. Science 283(5410), 2062-2064) and less than the equilibrium concentrations formed by ion irradiation. In addition to the ice film temperature, the current of electrons was varied between different experiments to determine the production and destruction of H₂O₂ as a function of both electron flux and ice temperature. Variation in current was found to have little effect on the results other than accelerating arrival at radiolytic equilibrium.     

Constraints from deuterium on the formation of icy bodies in the Jovian system and beyond

We consider the role of deuterium as a potential marker of location and ambient conditions during the formation of small bodies in our Solar system. We concentrate in particular on the formation of the regular icy satellites of Jupiter and the other giant planets, but include a discussion of the implications for the Trojan asteroids and the irregular satellites. We examine in detail the formation of regular planetary satellites within the paradigm of a circum-Jovian subnebula. Particular attention is paid to the two extreme potential subnebulae—“hot” and “cold”. In particular, we show that, for the case of the “hot” subnebula model, the D:H ratio in water ice measured from the regular satellites would be expected to be near-Solar. In contrast, satellites which formed in a “cold” subnebula would be expected to display a D:H ratio that is distinctly over-Solar. We then compare the results obtained with the enrichment regimes which could be expected for other families of icy small bodies in the outer Solar system—the Trojan asteroids and the irregular satellites. In doing so, we demonstrate how measurements by Laplace, the James Webb Space Telescope, HERSCHEL and ALMA will play an important role in determining the true formation locations and mechanisms of these objects.     

Diagnostics of Titan's stratospheric dynamics using Cassini/CIRS data and the 2-dimensional IPSL circulation model

The dynamics of Titan's stratosphere is discussed in this study, based on a comparison between observations by the CIRS instrument on board the Cassini spacecraft, and results of the 2-dimensional circulation model developed at the Institute Pierre-Simon Laplace, available at http://www.lmd.jussieu.fr/titanDbase [Rannou, P., Lebonnois, S., Hourdin, F., Luz, D., 2005. Adv. Space Res. 36, 2194-2198]. The comparison aims at both evaluating the model's capabilities and interpreting the observations concerning: (1) dynamical and thermal structure using temperature retrievals from Cassini/CIRS and the vertical profile of zonal wind at the Huygens landing site obtained by Huygens/DWE; and (2) vertical and latitudinal profiles of stratospheric gases deduced from Cassini/CIRS data. The modeled thermal structure is similar to that inferred from observations (Cassini/CIRS and Earth-based observations). However, the upper stratosphere (above 0.05 mbar) is systematically too hot in the 2D-CM, and therefore the stratopause region is not well represented. This bias may be related to the haze structure and to misrepresented radiative effects in this region, such as the cooling effect of hydrogen cyanide (HCN). The 2D-CM produces a strong atmospheric superrotation, with zonal winds reaching 200 m s⁻¹ at high winter latitudes between 200 and 300 km altitude (0.1-1 mbar). The modeled zonal winds are in good agreement with retrieved wind fields from occultation observations, Cassini/CIRS and Huygens/DWE. Changes to the thermal structure are coupled to changes in the meridional circulation and polar vortex extension, and therefore affect chemical distributions, especially in winter polar regions. When a higher altitude haze production source is used, the resulting modeled meridional circulation is weaker and the vertical and horizontal mixing due to the polar vortex is less extended in latitude. There is an overall good agreement between modeled chemical distributions and observations in equatorial regions. The difference in observed vertical gradients of C₂H₂ and HCN may be an indicator of the relative strength of circulation and chemical loss of HCN. The negative vertical gradient of ethylene in the low stratosphere at 15° S, cannot be modeled with simple 1-dimensional models, where a strong photochemical sink in the middle stratosphere would be necessary. It is explained here by dynamical advection from the winter pole towards the equator in the low stratosphere and by the fact that ethylene does not condense. Near the winter pole (80° N), some compounds (C₄H₂, C₃H₄) exhibit an (interior) minimum in the observed abundance vertical profiles, whereas 2D-CM profiles are well mixed all along the atmospheric column. This minimum can be a diagnostic of the strength of the meridional circulation, and of the spatial extension of the winter polar vortex where strong descending motions are present. In the summer hemisphere, observed stratospheric abundances are uniform in latitude, whereas the model maintains a residual enrichment over the summer pole from the spring cell due to a secondary meridional overturning between 1 and 50 mbar, at latitudes south of 40-50° S. The strength, as well as spatial and temporal extensions of this structure are a difficulty, that may be linked to possible misrepresentation of horizontally mixing processes, due to the restricted 2-dimensional nature of the model. This restriction should also be kept in mind as a possible source of other discrepancies.     

Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

Titan trace gaseous composition from CIRS at the end of the Cassini-Huygens prime mission

This paper reports on the results from an extensive study of all nadir-looking spectra acquired by Cassini/CIRS during the 44 flybys performed in the course of the nominal mission (2004-2008). With respect to the previous study (Coustenis, A., and 24 colleagues [2007]. Icarus 189, 35-62, on flybys TB-T10) we present here a significantly richer dataset with, in particular, more data at high northern and southern latitudes so that the abundances inferred here at these regions are more reliable. Our enhanced high-resolution dataset allows us to infer more precisely the chemical composition of Titan all over the disk. We also include improved spectroscopic data for some molecules and updated temperature profiles. The latitudinal distributions of all of the gaseous species are inferred. We furthermore test vertical distributions essentially for acetylene (C₂H₂) from CIRS limb-inferred data and from current General Circulation Models for Titan and compare our results on all the gaseous abundances with predictions from 1-D photochemical-radiative models to check the reliability of the chemical reactions and pathways.     

The initial responses of hot liquid water released under low atmospheric pressures: Experimental insights

Experiments have been performed to simulate the shallow ascent and surface release of water and brines under low atmospheric pressure. Atmospheric pressure was treated as an independent variable and water temperature and vapor pressure were examined as a function of total pressure variation down to low pressures. The physical and thermal responses of water to reducing pressure were monitored with pressure transducers, temperature sensors and visible imaging. Data were obtained for pure water and for solutions with dissolved NaCl or CO₂. The experiments showed the pressure conditions under which the water remained liquid, underwent a rapid phase change to the gas state by boiling, and then solidified because of removal of latent heat. Liquid water is removed from phase equilibrium by decompression. Solid, liquid and gaseous water are present simultaneously, and not at the 611 Pa triple point, because dynamic interactions between the phases maintain unstable temperature gradients. After phase changes stop, the system reverts to equilibrium with its surroundings. Surface and shallow subsurface pressure conditions were simulated for Mars and the icy satellites of the outer Solar System. Freezing by evaporation in the absence of wind on Mars is shown to be unlikely for pure water at pressures greater than c. 670 Pa, and for saline solutions at pressures greater than c. 610 Pa. The physical nature of ice that forms depends on the salt content. Ice formed from saline water at pressures less than c. 610 Pa could be similar to terrestrial sea ice. Ice formed from pure water at pressures less than c. 100 Pa develops a low thermal conductivity and a ‘honeycomb’ structure created by sublimation. This ice could have a density as low as c. 450 kg m⁻³ and a thermal conductivity as low as 1.6 W m⁻¹ K⁻¹, and is highly reflective, more akin to snow than the clear ice from which it grew. The physical properties of ice formed from either pure or saline water at low pressures will act to reduce the surface temperature, and hence rate of sublimation, thereby prolonging the lifespan of any liquid water beneath.     

The rotation of Janus and Epimetheus

Epimetheus, a small moon of Saturn, has a rotational libration (an oscillation about synchronous rotation) of 5.9°±1.2°, placing Epimetheus in the company of Earth’s Moon and Mars’ Phobos as the only natural satellites for which forced rotational libration has been detected. The forced libration is caused by the satellite’s slightly eccentric orbit and non-spherical shape.Detection of a moon’s forced libration allows us to probe its interior by comparing the measured amplitude to that predicted by a shape model assuming constant density. A discrepancy between the two would indicate internal density asymmetries. For Epimetheus, the uncertainties in the shape model are large enough to account for the measured libration amplitude. For Janus, on the other hand, although we cannot rule out synchronous rotation, a permanent offset of several degrees between Janus’ minimum moment of inertia (long axis) and the equilibrium sub-Saturn point may indicate that Janus does have modest internal density asymmetries.The rotation states of Janus and Epimetheus experience a perturbation every 4 years, as the two moons “swap” orbits. The sudden change in the orbital periods produces a free libration about synchronous rotation that is subsequently damped by internal friction. We calculate that this free libration is small in amplitude (<0.1°) and decays quickly (a few weeks, at most), and is thus below the current limits for detection using Cassini images.     

The atmospheric signature of Io's Prometheus plume and anti-jovian hemisphere: evidence for a sublimation atmosphere

Using the Hubble Space Telescope's Space Telescope Imaging Spectrograph we have obtained for the first time spatially resolved 2000-3000 Å spectra of Io's Prometheus plume and adjoining regions on Io's anti-jovian hemisphere in the latitude range 60° N-60° S, using a 0.1″ slit centered on Prometheus and tilted roughly 45° to the spin axis. The SO₂ column density peaked at 1.25×10¹⁷ cm⁻² near the equator, with an additional 5×10¹⁶ cm⁻² enhancement over Prometheus corresponding to a model volcanic SO₂ output of 10⁵ kg s⁻¹. Apart from the Prometheus peak, the SO₂ column density dropped fairly smoothly away from the subsolar point, even over regions that included potential volcanic sources. At latitudes less than ±30°, the dropoff rate was consistent with control by vapor pressure equilibrium with surface frost with subsolar temperature 117.3±0.6 K, though SO₂ abundance was higher than predicted by vapor pressure control at mid-latitudes, especially in the northern hemisphere. We conclude that, at least at low latitudes on the anti-jovian hemisphere where there are extensive deposits of optically-thick SO₂ frost, the atmosphere is probably primarily supported by sublimation of surface frost. Although the 45° tilt of our slit prevents us from separating the dependence of atmospheric density on solar zenith angle from its dependence on latitude, the pattern is consistent with a sublimation atmosphere regardless of which parameter is the dominant control. The observed drop in gas abundance towards higher latitudes is consistent with the interpretation of previous Lyman alpha images of Io as indicating an atmosphere concentrated at low latitudes. Comparison with previous disk-resolved UV spectroscopy, Lyman-alpha images, and mid-infrared spectroscopy suggests that Io's atmosphere is denser and more widespread on the anti-jovian hemisphere than at other longitudes. SO₂ gas temperatures were in the range of 150-250 K over the majority of the anti-jovian hemisphere, consistent with previous observations. SO was not definitively detected in our spectra, with upper limits to the SO/SO₂ ratio in the range 1-10%, roughly consistent with previous observations. S₂ gas was not seen anywhere, with an upper limit of 7.5×10¹⁴ cm⁻² for the Prometheus plume, confirming that this plume is significantly poorer in S₂ than the Pele plume (S₂ /SO₂<0.005, compared to 0.08-0.3 at Pele). In addition to the gas absorption signatures, we have observed continuum emission in the near ultraviolet (near 2800 Å) for the first time. The brightness of the observed emission was directly correlated with the SO₂ abundance, strongly peaking in the equatorial region over Prometheus. Emission brightness was modestly anti-correlated with the jovian magnetic latitude, decreasing when Io intersected the torus centrifugal equator.