Coma expansion and photometry of comet Bowell (1980b)

Optical and infrared observations of comet Bowell are presented. The optical observations indicate that the solid grain coma is expanding at only 0.9 ± 0.2 m sec^?1. This is two orders of magnitude slower than the local gas sound speed and may suggest that gas drag is not responsible for stripping the grains from the nucleus. The hypothesis of “electrostatic snap-off” is tentatively advanced to account for the ejection of the grains. Alternatively, the grains may have an unusual size distribution. The extrapolated motion of the grains suggests that the bulk of the coma was formed when the comet was at a heliocentric distance $\text{R ? 10}\text{AU}$. Any water ice in the nucleus would be too cold to give rise to the observed grain coma by equilibrium sublimation at this R. Further evidence against the production of the grain coma by equilibrium sublimation of the nucleus is provided by broadband (J) photometric observations. Almost all of the observed photometric variations of comet Bowell can be ascribed to geometric effects. Simple models indicate that the total grain cross section has been nearly constant since the time of the earliest observations. The present observations, which suggest that water ice sublimation does not control either the optical morphology or the near infrared photometric behavior of comet Bowell, are contrasted with reported high OH production rates. It is concluded that the grain coma may be largely a relic of activity occurring on the nucleus at $\text{R ? 10}\text{AU}$ while the OH may indicate sublimation from the nucleus near perihelion and from coma grains near $\text{R ? 4.6}\text{AU}$.     

A spectroscopic study of the surfaces of Saturn's large satellites: H2O ice, tholins, and minor constituents

We present spectra of Saturn's icy satellites Mimas, Enceladus, Tethys, Dione, Rhea, and Hyperion, 1.0-2.5 μm, with data extending to shorter (Mimas and Enceladus) and longer (Rhea and Dione) wavelengths for certain objects. The spectral resolution (R=λ/Δλ) of the data shown here is in the range 800-1000, depending on the specific instrument and configuration used; this is higher than the resolution (R=225 at 3 μm) afforded by the Visual-Infrared Mapping Spectrometer on the Cassini spacecraft. All of the spectra are dominated by water ice absorption bands and no other features are clearly identified. Spectra of all of these satellites show the characteristic signature of hexagonal H₂O ice at 1.65 μm. We model the leading hemisphere of Rhea in the wavelength range 0.3-3.6 μm with the Hapke and the Shkuratov radiative transfer codes and discuss the relative merits of the two approaches to fitting the spectrum. In calculations with both codes, the only components used are H₂O ice, which is the dominant constituent, and a small amount of tholin (Ice Tholin II). Tholin in small quantities (few percent, depending on the mixing mechanism) appears to be an essential component to give the basic red color of the satellite in the region 0.3-1.0 μm. The quantity and mode of mixing of tholin that can produce the intense coloration of Rhea and other icy satellites has bearing on its likely presence in many other icy bodies of the outer Solar System, both of high and low geometric albedos. Using the modeling codes, we also establish detection limits for the ices of CO₂ (a few weight percent, depending on particle size and mixing), CH₄ (same), and NH₄OH (0.5 weight percent) in our globally averaged spectra of Rhea's leading hemisphere. New laboratory spectral data for NH₄OH are presented for the purpose of detection on icy bodies. These limits for CO₂, CH₄, and NH₄OH on Rhea are also applicable to the other icy satellites for which spectra are presented here. The reflectance spectrum of Hyperion shows evidence for a broad, unidentified absorption band centered at 1.75 μm.     

On the dectectability of icy grains in the comae of comets

Evaporation of icy grains over the distance scale of the visible cometary coma sets very specific limits on their temperature. Unless the grains are very pure water ice, the maximum size of an icy grain halo will be limited to a few hundred kilometers at heliocentric distances ?2.5 AU. It is unlikely that the 1.5- or 2-μm ice band could be detected in the scattering by icy grains. Detection of the 3?μm ice band might be possible in comets which display a coma at large heliocentric distances.     

How to make metal-poor stars,redden ob associations and grow mantles on grains

Three consequences of the existence of grains with metal-rich ice mantles are considered: (i) The production of metal-poor stars by expulsion of protostellar grains by radiation pressure during star formation. (ii) The effects of these expelled grains in reddening massive stars in an OB association. (iii) The production of the icy mantles on grains in OB associations.     

Long- and Short-Term Photometric Behavior of Comet Hyakutake (1996 B2)

Narrowband filter photometry observations of Comet Hyakutake (1996 B2) were used to investigate this comet's short-term variability as well as its behavior for the apparition as a whole. Utilizing measurements obtained on a total of 13 nights between February 9, 1996, and April 14, 1996, we find that the heliocentric distance (r_H) dependence of the production rates of OH and NH were much shallower than those for either the carbon-bearing species or the visible dust. Based on the OH measurements, the derived water r_H-dependence was also significantly less steep than expected from a basic water vaporization model and required an effective active surface area of about 29 km² at r_H=1.8 AU, 16 km² at r_H=1 AU, and only 13 km² at r_H=0.6 AU. This decrease in active area may be due to seasonally induced variations of a heterogeneous surface, or due to a decreasing contribution of gas from icy grains in the innermost coma. The relative abundances of the minor gas species place Hyakutake into the “typical” category of comets in the A'Hearn et al. (1995, Icarus118, 223-270) taxonomic classification system. The spectrum is generally redder at shorter wavelengths throughout the apparition; however, the dust color progressively changes from being significantly reddened (37%/1000 Å) at large r_H to near-solar at small r_H. This change of color with distance implies a significant change in grain sizes or a changing proportion between two or more grain populations.A major outburst was initiated near March 19.9, just prior to the comet's close approach to Earth. The characteristic recovery from the outburst differed among the observed species, with OH recovering most rapidly, essentially returning to its baseline values by March 25. The spatial radial fall-off of OH throughout this interval was consistent with the expected nominal spatial distribution, while CN and C₂ displayed fall-offs consistent with a distributed source, and the dust fall-off was significantly less steep than 1/ρ, possibly due to fragmenting grains. Rotational lightcurve amplitudes were largest for the OH, NH, and dust, again consistent with the carbon-bearing species primarily originating from a distributed source. Significant variations were observed in the lightcurve amplitude and phase shifts as functions of aperture size. Finally, a refined value for the rotation period of 0.2614±0.0003 day was determined.     

Icy Galilean satellite reflectance spectra: Less ice on Ganymede and Callisto?

Recent interpretations of the reflectance spectra of the icy Galilean satellites (Europa, Ganymede, and Callisto) have implied very ice-rich surfaces, as high as 90 wt% ice even on the dark surface of Callisto. A reevaluation of the spectra, taking into account the depth of the 3-μm fundamental water ice absorption feature as well as the shorter wavelength bands, suggests that the spectra of at least Ganymede and Callisto may also be consistent with much lower ice abundances if the ice is segregated from the nonicy material. Reasonable fits to all band depths (including the shallow 1.04- and 1.25-μm bands) are obtained with around 50% areal coverage of ice on Ganymede and 10% on Callisto, the rest of the surface being occupied by carbonaceous chondrite-like material which has a strong 3-μm absorption due to bound water. Europa's spectrum probably indicates a homogeneous icy surface. The darkness beyond 3 μm, and lack of a 3.6-μm peak, for all three objects may be consistent with the presence of small quantities of sulfuric acid on the satellite surfaces.     

Dust ring formation due to ice sublimation of radially drifting dust particles under the Poynting-Robertson effect in debris disks

In a disk with a low optical depth, dust particles drift radially inward by the Poynting-Robertson (P-R) drag rather than are blown out by stellar radiation pressure following destructive collisions. We investigate the radial distribution of icy dust composed of pure ice and refractory materials in dust-debris disks taking into account the P-R drag and ice sublimation. We find that icy dust particles form a dust ring by their pile-ups at the edge of their sublimation zone, where they sublime substantially at the temperature 100-110 K. The distance of the dust ring is 20-35 AU from the central star with its luminosity L_??30L_⊙ and 65(L_?/100L_⊙)^1/2 AU for L_??30L_⊙, where L_⊙ is the solar luminosity. The effective optical depth is enhanced by a factor of 2 for L_??100L_⊙ and more than 10 for L_??100L_⊙. The optical depth of the outer icy dust disk exceeds that of the inner disk filled with refractory particles, namely, the residue of ice sublimation, which are further subjected to the P-R effect. As a result, an inner hole is formed inside the sublimation zone together with a dust ring along the outer edge of the hole.     

A VLA search for 2-cm continuum radiation from Comet Halley

We used the VLA² to search for continuum emission from icy grains in a halo around Comet Halley in mid November 1985. We found the 3σ upper limit to the 2-cm flux density from the comet to be 1 × 10⁻⁴ Jy, which is consistent with the detections at 1.3 and 3.5 mm by W.J. Altenhoff et al. (1986, Astron. Astrophys., in press) only if the emission comes from particles which do not radiate efficiently at centimeter wavelengths. These particles could be slightly dirty submillimeter-sized icy grains or small refractory grains.     

Ganymede,Europa, Callisto,and Saturn's rings: Compositional analysis from reflectance spectroscopy

The reflectance spectra of Ganymede, Europa, Callisto, and Saturn's rings are analyzed using recent laboratory reflectance studies of water frost, water ice, and water and mineral mixtures. It is found that the spectra of the icy Galilean satellites are characteristic of water ice (e.g., ice blocks or possibly very large ice crystals ? 1 cm) or frost on ice rather than pure water frost, and that the decrease in reflectance at visible wavelengths is caused by other mineral grains in the surface. The spectra of Saturn's rings are more characteristic of water frost with some other mineral grains mixed in the frost but not on the surface. The impurities on all these objects are not in spectrally isolated patches but appear to be intimately mixed with the water. The impurity grains appear to have reflectance spectra typical of minerals containing Fe³⁺. Some carbonaceous chondrite meteorite spectra show the necessary spectral shape. Ganymede is found to have more water ice on the surface than previously thought (～90 wt%), as is Callisto (30–90 wt%). The surface of Europa has a vast frozen water surface with only a few percent impurities. Saturn's rings also have only a few percent impurities. The amount of bound water or bound OH for these objects is 5 ± 5 wt% averaged over the entire surface. Thus with the small amount of nonicy material present on these objects, no hydrated minerals can be ruled out. A new absorption feature is identified in Ganymede, Callisto, and probably Europa at 1.5 μm which is also seen in the spectra of Io but not in Saturn's rings. This feature has not been seen in laboratory studies and its cause is unknown.     

Radiogenic heating of comets by 26Al and implications for their time of formation

The effect of radiogenic heating on the thermal evolution of spherical icy bodies with radii 1 km < R < 100 km was investigated. The radioisotopes considered were 26Al, 40K, 232Th, 235U, and 238U. Except for the 26Al abundance, which was varied, the other initial abundances were kept fixed, at values derived from those of chondritic meteorites and corresponding to a gas-to-dust ratio of 1. The initial models were homogeneous and isothermal (To = 10 K) amorphous ice spheres, in a circular orbit at 10(4) AU from the Sun. The main object of this study was to examine the conditions under which the transition temperature from amorphous into cubic ice (Ta = 137 K) would be reached. It was shown that the influence of the short-lived radionuclide 26Al dominates the effect of other radioactive species for bodies of radii up to approximately 50 km. Consequently, if we require comets to retain their ice in amorphous form, as suggested by observations, an upper limit of approximately 4 x 10(-9) is obtained for the initial 26Al abundance in comets, a factor of 100 lower than that of the inclusions in the Allende meteorite. A lower limit for the formation time of comets may thus be derived. The possibility of a coexistence of molten cometary cores and extended amorphous ice mantles is ruled out. Larger icy spheres (R > 100 km) reached Ta even in the absence of 26Al, due to the decay of the other radionuclides. As a result, a crystalline core formed whose relative size depended on the composition assumed. Thus the outermost icy satellites in the solar system, which might have been formed of ice in the amorphous state, have probably undergone crystallization and may have exhibited eruptive activity when the gas trapped in the amorphous ice was released (e.g., Miranda).     

A process to make massive ice in the martian regolith using long-term diffusion and thermal cracking

This paper describes a “simple standard” model of water transport through regolith that includes diffusive migration and phase changes driven by damped seasonal temperature waves. A hitherto unused first-order process is then added, that can produce ice densities much greater than those allowed by the initial dry porosity. Voids are produced in cooling icy regolith when tensile stresses exceed the cracking threshold . These stresses build up through an interaction of thermal contraction and elastic-plastic response. When the cracks open up after tensile failure there is purely thermal void enhancement and subsequent reduction as the regolith warms again. When the cracks are open the porosity is increased and they partially fill with ice crystals. Thus the void reduction on warming cannot go back to the original zero point and the bulk density of ice is increased with each temperature cycle. The cracking and thermal adjustment happen at scales of meters to millimeters. The large cracks can occur in pure ice and/or homogeneous icy material and the smaller cracks are produced by rock cobbles, pebbles, and grains having a much smaller coefficient of thermal expansion than ice. Thus a hierarchy of cracks and voids forms each temperature cycle and augments the ice content. The process can take the upper few meters of a pore-saturated icy soil from 28% by mass ice content to 70% in 10 Ma. This mechanism and the seasonal temperature cycle can plausibly produce massive ice deposits in the upper few meters of Mars' high-latitude regolith by diffusion and also keep the massive-ice regolith effectively porous to water vapor transport. The obliquity cycle can produce tensile stresses nearing 2 MPa down to depth so even deeper cracking could be happening.     

Radiation chemistry of H2O + O2 ices

The chemistry and spectroscopy of proton-irradiated H₂O + O₂ ices have been investigated in relation to the production of oxidants in icy satellite surfaces. Hydrogen peroxide (H₂O₂), ozone (O₃), and the hydroperoxy (HO₂) and hydrogen trioxide (HO₃) radicals have all been observed, and their temperature and dose dependent production trends have been measured. We find that O₂ aggregates form during the growth of H₂O + O₂ ice films, and the presence of these aggregates greatly affects the HO₂ and H₂O₂ yields. In addition, we have found that the position of the spectral maximum of the ν₃ vibration of O₃ shifts with ice composition, giving an indication of the degree of dispersion of O₃ molecules within the ice. We discuss the relevance of these measurements to icy satellite surfaces.     

Ice in the Southern Coalsack

To better understand the conditions under which ice mantles form on grains in molecular clouds, three globules in the Southern Coalsack have been searched for the presence of H₂O ice. Given the total lack of star formation in the Coalsack, it is an ideal site for studying unprocessed icy molecular mantles. In our sample of eight field stars lying behind the Coalsack we detect strong H₂O ice absorption in the lines of sight to two stars and possible weak absorption in four others. We estimate H₂O ice column densities or upper limits for these lines of sight. Compared to dark clouds such as Taurus, the Coalsack H₂O ice column densities are lower than expected given the quiescent nature of the Coalsack region. It is possible that the chemical evolution of the Coalsack may simply be at too early a stage for significant ice mantles to appear on the grains, except perhaps in the densest parts of some of the globules. Alternatively, the presence or absence of ice absorption may be related to the distribution of dust along each line of sight, specifically, the relative contributions of dense globules and a more extended diffuse component. For example, our observations are consistent with an ice threshold extinction similar to that observed in the Taurus dark cloud if extinction amounting to   A _V∼5  towards Globules 2 and 3 arises in the extended component. Globule 1 appears to have no extended component.     

Circular polarimetry and the line of sight to the Becklin–Neugebauer object

The 3.1-μm absorption feature of water-ice has been observed spectroscopically in many molecular clouds and, when it has been observed spectropolarimetrically, usually a corresponding polarization feature is seen. Typically, on these occasions, and particularly for the Becklin–Neugebauer (BN) object, a distinct position angle shift between the feature and continuum is seen, which indicates both a fractionation of the icy material and a changing alignment direction along the line of sight.
Here, the dependence of circular polarimetry on fractionation along the line of sight is investigated and it is shown that the form of its spectrum, together with the sign of the position angle shift, indicates where along the line of sight the icy material lies. More specifically, a coincidence between the sign of the position angle displacement in the ice feature, measured north through east, and that of the circular polarization ice feature means that the icy grains are overlaid by bare grains. Some preliminary circular polarimetry of BN has this characteristic, and a similar situation is found in the only two other cases for which relevant observations so far exist.     

Thermal convection in ice-I shells of Titan and Enceladus

Cassini-Huygens observations have shown that Titan and Enceladus are geologically active icy satellites. Mitri and Showman [Mitri, G., Showman, A.P., 2005. Icarus 177, 447-460] and McKinnon [McKinnon, W.B., 2006. Icarus 183, 435-450] investigated the dynamics of an ice shell overlying a pure liquid-water ocean and showed that transitions from a conductive state to a convective state have major implications for the surface tectonics. We extend this analysis to the case of ice shells overlying ammonia-water oceans. We explore the thermal state of Titan and Enceladus ice-I shells, and also we investigate the consequences of the ice-I shell conductive-convective switch for the geology. We show that thermal convection can occur, under a range of conditions, in the ice-I shells of Titan and Enceladus. Because the Rayleigh number Ra scales with δ³/ηb, where δ is the thickness of the ice shell and ηb is the viscosity at the base of the ice-I shell, and because ammonia in the liquid layer (if any) strongly depresses the melting temperature of the water ice, Ra equals its critical value for two ice-I shell thicknesses: for relatively thin ice shell with warm, low-viscosity base (Onset I) and for thick ice shell with cold, high-viscosity base (Onset II). At Onset I, for a range of heat fluxes, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. Switches between these states can cause large, rapid changes in the ice-shell thickness. For Enceladus, we demonstrate that an Onset I transition can produce tectonic stress of ∼500 bars and fractures of several tens of km depth. At Onset II, in contrast, we demonstrate that zero equilibrium states exist for a range of heat fluxes. For a mean heat flux within this range, the satellite experiences oscillations in surface heat flux and satellite volume with periods of ∼50-800 Myr even when the interior heat production is constant or monotonically declining in time; these oscillations in the thermal state of the ice-I shell would cause repeated episodes of extensional and compressional tectonism.     

Do explosive ice ejections occur on Jupiter’s and Saturn’s satellites?

The possibility of an explosive mechanical instability of ice (the Bridgman effect) in the thick icy shells of Jupiter’s and Saturn’s satellites is discussed in principle. The Bridgman effect is an explosive instability of dielectric solid bodies, which disintegrate into microscopic fragments under a quasistatic uniaxial loading in open compression systems at high pressures. The explosive instabilities of ice recently discovered in laboratory experiments with the Bridgman effect are also expected to occur in the extensive deep layers of the shells of icy planetary satellites (for example, in the case of episodical formation of major cracks in their lithospheres due to tidal forces, nonsynchronous rotation of the satellites, or extremely powerful impacts). The depths of occurrence of mechanically unstable ice in the thick crusts of Ganymede, Europa, and Titan, taken as examples, are crudely estimated using a pure-ice model without a possible ammonia admixture. The estimated thickness of the explosive-instability zone in the icy crust of Ganymede (under the assumption that the crust is ～75 km thick) ranges from ～7 to ～27 km at depths from ～40 to ～67 km, depending on the scaling parameter E = 0.2–1. This parameter relates the experimentally determined thicknesses of the ice samples in which the Bridgman effect occurs under laboratory conditions to the expected thicknesses of the explosively unstable layers in the envelopes of the icy satellites. Explosive effects are possible not throughout the entire thickness of the unstableice layer but only within some part of it, several centimeters to several tens of meters in thickness. According to the estimated location of the unstable layer in the crust of Europa (for an assumed crust thickness of ～30 km), such a layer can exist only at scaling factors E < 0.6 at depths ranging from ～21 to ～28 km. For Titan, if its crust is ～100 km thick, the thickness of the unstable layer is similarly estimated to range from ～15 to ～55 km at depths from ～37 to ～92 km for a scaling parameter E lying within the range 0.2–1. At E 0.2, which is quite possible, explosive instabilities of ice could also be expected on the Earth, in the icy shells of Antarctica and Greenland at depths from ～1 to ～1.5 km.     

Far ultraviolet spectral properties of Saturn’s rings from Cassini UVIS

Spectra taken by the Cassini Ultraviolet Imaging Spectrograph (UVIS) of Saturn’s C ring, B ring, Cassini Division, and A ring have been analyzed in order to characterize ring particle surface properties and water ice abundance in the rings. UVIS spectra sense the outer few microns of the ring particles. Spectra of the normalized reflectance (I/F) in all four regions show a characteristic water ice absorption feature near 165 nm. Our analysis shows that the fractional abundance of surface water ice is largest in the outer B ring and decreases by over a factor of 2 across the inner C ring. We calculate the mean path length of UV photons through icy ring particle regolith and the scattering asymmetry parameter using a Hapke reflectance model and a Shkuratov reflectance model to match the location of the water ice absorption edge in the data. Both models give similar retrieved values of the photon mean length, however the retrieved asymmetry (g) values are different. The photon mean path lengths are nearly uniform across the B and A rings. Shortward of 165 nm the rings exhibit a slope that turns up towards shorter wavelengths, while the UV slope of 180/150 nm (reflectance outside the water absorption ratioed to that inside the absorption band) tracks I/F with maxima in the outer B ring and in the central A ring. Retrieved values of the scattering asymmetry parameter show the regolith grains to be highly backscattering in the FUV spectral regime.     

The Water-Ammonia Phase Diagram up to 300 MPa: Application to Icy Satellites

We performed high-pressure experiments on the crystallization of water ice I and III in the ammonia-water (NH₃)_x(H₂O)_(1−x) system, and apply the results to the interiors of icy bodies in the Solar System. Phase equilibrium lines between an entirely liquid solution and a liquid solution in which water ice forms (liquidus lines) were determined for ammonia concentration by mass X equal to 0.034, 0.0472, 0.111, 0.176, and 0.229. Growth-melting of ice I as well as ice III crystals were observed. Application of the results to icy satellites that are potential bearers of ammonia shows that ammonia admixture decreases the depth of the liquidus surface. A shift of the liquidus temperature within a satellite depends on three parameters: the ammonia concentration, X; the temperature gradient, α; and the product of density and gravity, ρg.     

First experimental studies of large samples of gas-laden amorphous “cometary” ices

In a unique machine, the first of its kind, large (200 cm² × 10 cm) samples of gas-laden amorphous ice were prepared at 80 K and 10⁻⁵ Torr. The sample consisted of a fluffy agglomerate of 200-μm ice grains, similar to what is presumed to be the structure of comet nuclei. The sample was heated from above by IR radiation. The properties studied were gas content in the ice and its emanation from the ice upon warming and bearing on the gas/water vapor ratio observed in cometary comae vs this ratio in cometary nuclei and the effect of internal trapped gas on the thermal conductivity of the ice and the density and mechanical properties of pure ice vs gas-laden ice. These findings might have significance for the interpretation of comet observations, the forthcoming ESA’s Rosetta space mission to Comet 46P/Wirtanen in 2012, and to other comet missions.     

PFS-MEX observation of ices in the residual south polar cap of Mars

The Mars Express spacecraft has a highly inclined orbit around Mars and so has been able to observe the south pole of Mars in illuminated conditions at the end of the southern summer (L_s=330). Spectra from the planetary Fourier spectrometer (PFS) short wavelength (SW) channel were recorded over the permanent ice cap to study its composition in terms of CO₂ ice and H₂O ice. Models are fitted to the observed data, which include a spatial mixture of soil (not covered by ice) and CO₂ frost (with a specific grain size and a small amount of included dust and H₂O ice). Two different kinds of spectra were observed: those over the permanent polar cap with almost pure CO₂ ice, negligible water ice, no soil fraction required, and bright; and those over mixed terrain (at the edge of the cap or near troughs) containing a significant soil spatial fraction, more water ice and smaller CO₂ grain size. The amount of water ice given by fits to scaled albedo models is less than 10 ppm by weight. When using multi-stream reflectance models with the appropriate lighting geometry, the water amount must be 2-5 times greater than the albedo fit (less than 50 ppm). At the periphery of the residual polar cap, we found a region almost completely covered by water frost, modeled as a mixture of micron-sized and sub-mm sized grains. Our result using a granular mixture of micron-sized grains of water ice and dust with the CO₂ grains is different from the modeling of OMEGA polar cap observations using molecular mixtures.