Experimental simulation of the formation of non-circular active depressions on Comet Wild-2 and of ice grain ejection from cometary surfaces

Following the tracing of jets emanating from Comet Wild-2 to depressions in the ice by Brownlee et al. [2004. The Stardust—A successful encounter with the remarkable Comet Wild 2. Lunar Planet. Sci. 35. Abstract 1981], we demonstrated experimentally the formation of depressions and chaotic terrain on comet analogs when gas is released from underlying ice pockets. We also demonstrated experimentally the ejection of ice grains into the experimental cometary “coma.”     

Comparison between the findings of Deep Impact and our experimental results on large samples of gas-laden amorphous ice

The findings of Deep Impact on the structure and composition of Tempel-1 are compared with our experimental results on large (20 cm diameter and up to 10 cm high) samples of gas-laden amorphous ice. The mechanical ∼tensile strength inferred for Tempel-1: ∼65 Pa is 30 to 60 times smaller than our experimental findings of 2-4 kPa. This means that Tempel-1 is even fluffier than our very fluffy, talcum like, ice sample. The thermal inertia: is very close to our value of 80. The density of , is close to our value of 250-300 kg m⁻³, taking into account an ice/silicate ratio of 1 in the comet, while we study pure ice. Surface morphological features, such as non-circular depressions, chaotic terrain and smooth surfaces, were observed in our experiments. The only small increase in the gas/water vapor ratio pre- and post-impact, suggest that in the area excavated by the impactor, the 135 K front did not penetrate deeper than a few meters. Altogether, the agreement between the findings of Deep Impact and our experimental results point to a loose agglomerate of ice grains (with a silicate-organic core), which was formed by a very gentle aggregation of the ice grains, without compaction.     

Activity of comets at large heliocentric distances pre-perihelion

We present observational data for two long-period and three dynamically new comets observed at heliocentric distances between 5.8 to 14.0 AU. All of the comets exhibited activity beyond the distance at which water ice sublimation can be significant. We have conducted experiments on gas-laden amorphous ice samples and show that considerable gas emission occurs when the ice is heated below the temperature of the amorphous-crystalline ice phase transition (T∼137 K). We propose that annealing of amorphous water ice is the driver of activity in comets as they first enter the inner Solar System. Experimental data show that large grains can be ejected at low velocity during annealing and that the rate of brightening of the comet should decrease as the heliocentric distance decreases. These results are consistent with both historical observations of distant comet activity and with the data presented here. If observations of the onset of activity in a dynamically new comet are ever made, the distance at which this occurs would be a sensitive indicator of the temperature at which the comet had formed or represents the maximum temperature that it has experienced. New surveys such as Pan STARRS, may be able to detect these comets while they are still inactive.     

Gas trapping in water ice at very low deposition rates and implications for comets

The effect of water ice formation temperature and rate of ice deposition on a cold plate on the amount of trapped argon (equivalent to CO), and the ratios of Ar/Kr/Xe trapped in the water ice were studied at 50, 27 and 22 K and at ice formation rates ranging over four orders of magnitude, from 10⁻¹ to 10⁻⁵ μm min⁻¹. Contrary to our previous conclusions that cometary ices were formed at 50-60 K, we now conclude that these ices were formed at about 25 K. At 25 K the enrichment ratios for Ar, Kr, and Xe remained the same as those at 50 K, reinforcing our suggestion of cometary contribution of these noble gases to the atmospheres of Earth and Mars.     

The composition and size distribution of the dust in the coma of Comet Hale-Bopp

We discuss the composition and size distribution of the dust in the coma of Comet Hale-Bopp. We do this using a model fit for the infrared emission measured by the Infrared Space Observatory (ISO) and the measured degree of linear polarization of scattered light at various phase angles and wavelengths. The effects of particle shape on the modeled optical properties of the dust grains are taken into account. Both the short wavelength (7-44 μm) and the long wavelength (44-120 μm) infrared spectrum are fitted using the same dust parameters, as well as the degree of linear polarization at twelve different wavelengths in the optical to near-infrared domains. We constrain our fit by forcing the abundances of the major rock forming chemical elements to be equal to those observed in meteorites. The infrared spectrum at long wavelengths reveals that large grains are needed in order to fit the spectral slope. The size and shape distribution we employ allows us to estimate the sizes of the crystalline silicates. The ratios of the strength of various forsterite features show that the crystalline silicate grains in Hale-Bopp must be submicrometer-sized. On the basis of our analysis the presence of large crystalline silicate grains in the coma can be excluded. Because of this lack of large crystalline grains combined with the fact that we do need large amorphous grains to fit the emission spectrum at long wavelengths, we need only approximately 4% of crystalline silicates by mass (forsterite and enstatite) to reproduce the observed spectral features. After correcting for possible hidden crystalline material included in large amorphous grains, our best estimate of the total mass fraction of crystalline material is ∼7.5%, which is significantly lower than deduced in previous studies in which the typical derived crystallinity is ∼20-30%. The implications of this low abundance of crystalline material on the possible origin and evolution of the comet are discussed. We conclude that the crystallinity we observe in Hale-Bopp is consistent with the production of crystalline silicates in the inner Solar System by thermal annealing and subsequent radial mixing to the comet forming region (∼30 AU).     

An experimental study of the formation of an ice crust and migration of water vapor in a comet's upper layers

From recent close encounters with Comets Wild-2 and Tempel 1 we learned that their surfaces are very rugged and no simple uniform layers model can be applied to them. Rather, a glaciological approach should be applied for describing their surface features and behavior. Such intrinsically rugged surface is formed in our large scale experiments, where an agglomerate of ∼200 μm gas-laden amorphous ice particles is accumulated to form a 20 cm diameter and few cm high ice sample. The density, tensile strength and thermal inertia of our ice sample were found to be very close to those found by Deep Impact for Comet Tempel 1: density 250-300 kg m⁻³ vs DI 350-400 kg m⁻³; tensile strength 2-4 kPa vs DI 1-10 kPa; thermal inertia 80 W K⁻¹ m⁻² s^1/2 vs <100 W K⁻¹ m⁻² s^1/2 and <50 W K⁻¹ m⁻² s^1/2. From the close agreement between the thermal inertias measured in our ice sample, which had no dust coverage and that of Comet Tempel 1, we deduce that the low thermal inertia is an intrinsic property of the fluffy structure of the ice as a result of its low density, with an addition by the broken terrain and not due to the formation of a dust layer. Upon warming up of the ice, water vapor migrates both outward into the coma and inward. Reaching cooler layers, the water vapor condenses, forming a denser ice crust, as we show experimentally. We also demonstrate the inward and outward flow of water vapor in the outer ice layers through the exchange between layers of D₂O ice and H₂O ice, to form HDO.     

Crystallization of ice in Comet 17P/Holmes: Probably not responsible for the explosive 2007 megaburst

Our work was inspired by the recent brightening of Comet 17P/Holmes. The recently observed increase in brightness of this comet was correlated with emission of dust, probably larger in mass than the dust mantle of the nucleus. We analyzed the hypothesis that the comet can eject a large mass of dust due to non-uniform crystallization of amorphous water ice. For this purpose, we simulated the evolution of a model nucleus on the orbit of Comet 17P/Holmes. The nucleus is composed of water ice and dust and has the shape of an elongated ellipsoid. The simulations include crystallization of amorphous ice in the nucleus, changes in the dust mantle thickness, and changes in the nucleus orientation in space. Our computations indicate that: (i) ejection of the dust cover triggers crystallization of ice independently on the material properties of the nucleus; (ii) moderate changes in the nucleus orientation (∼50°) may result in an acceleration of the crystallization of ice in the northern hemisphere, while a rather large change in the orientation (∼120°) is needed to cause a significant jump of the crystallization front in the southern hemisphere, where the emission of dust during the recent brightening was strongest. We investigated the possible reason for an explosion and we have found that the crystallization of the water ice itself is probably not sufficient.     

The contribution of grains to the activity of comets: II. The brightness of the coma

We use the model of grain behavior in the coma developed by Beer et al. [Beer, E.H., Podolak, M., Prialnik, P., 2006. Icarus 180, 473-486] to compute the contribution of ice grains to the brightness of the coma. The motion of an ice grain along the comet-Sun axis is computed, taking into account gas drag, the gravity of the nucleus, and radiation pressure of sunlight. The sublimation of the grains is also included. We assume that the maximum distance that a grain travels along this axis is indicative of the size of the coma, and we compute the resultant brightness as a function of heliocentric distance. The results are then compared to observations.     

The United Kingdom Infrared Telescope Deep Impact observations: Light curve, ejecta expansion rates and water spectral features

We present results from the United Kingdom Infrared Telescope observations of the impact of Deep Impact with Comet 9P/Tempel 1, on July 4, 2005 UT. These observations were carried out in conjunction with the worldwide observing campaign co-ordinated by K.J. Meech [Meech, K.J., and 208 colleagues, 2005. Science 310, 265-269]. The UKIRT team was the first to observe and announce the successful impact. At 05:50:52 (±2.5 s) UT the visible camera that is used to guide the telescope on the comet showed the start of a rapid rise in intensity, such that the visible brightness of Tempel 1 approximately doubled in 70 s. After that time there was a steady increase in the visible flux from the comet until it reached a maximum around 35 min post-impact, at which point it was more than ten times its original intensity. From an average of the time to maximum brightness and the time to noticeable intensity decline, we deduce that the material ejected by the impact expanded with a range of velocities between ∼125 and ∼390 m/s. We also observed water emission lines in the spectral region from 2.8945 to 2.8985 μm. We noted several water lines, which are known to be pumped by sunlight. But there was a lower intensity spectral component, which we propose may result from solar heating of icy grains freshly exposed by the impact.     

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.     

The contribution of icy grains to the activity of comets: I. Grain lifetime and distribution

We have developed a computer code (GEM—grain evolution model) to simulate the behavior of ice grains in a comet coma. The grains are assumed to be composed of water-ice with an admixture of dark material (“dirt”). An initial size distribution of grains is assumed to be ejected from the nucleus. The ejected mass is taken to be proportional to the rate of gas production by the nucleus. The efficiency for absorption and re-radiation of sunlight is computed from Mie scattering theory. The grain temperature and sublimation rate at a given heliocentric distance is then derived from energy balance considerations. The evolution of the grain size distribution is followed as a function of distance from the nucleus.     

Detection of large grains in the coma of Comet C/2001 A2 (LINEAR) from Arecibo radar observations

Arecibo S-band () radar observations of Comet C/2001 A2 (LINEAR) on 2001 July 7-9 showed a strong echo from large coma grains. This echo was significantly depolarized. This is the first firm detection of depolarization in a grain-coma radar echo and indicates that the largest grains are at least λ/2π or 2 cm in radius. The grains are moving at tens of m s⁻¹ with respect to the nucleus. The nondetection of the nucleus places an upper limit of 3 km on its diameter. The broad, asymmetric echo power spectrum suggests a fan of grains that have a steep (differential number ∼a⁻⁴) size distribution at cm-scales, though the observed fragmentation of this comet complicates that picture.     

Formation of smooth terrains on Comet Tempel 1

We suggest that the regions of smooth terrain which were observed on Comet 9P/Tempel 1 by the Deep Impact spacecraft were formed by blowing ice grains in an outburst of gas from the comet interior. When gas is released from 10 to 20 m deep layers which were heated to 135 K, it is released quiescently onto the surface by individual conduits. If large amounts of gas are released, the drainage system cannot release them fast enough and wider interconnected channels are formed, leading to sudden outburst of gas. Instability triggering a sudden shift of flow is well known in subglacial drainage of water. The ballistic trajectory of the ice particles reach a distance of 3 km in the atmosphereless comet, whose gravity is 0.034 cm s⁻¹, if ejected at an angle of 45° at a speed of 95 cm s⁻¹. This speed is close to the speeds measured in laboratory experiments: 167, 140×sini and 167 cm s⁻¹, for particles of 0.3, 1000 and 14-650 μm, respectively. Blowing of ice grains can overcome the 1650 m long horizontal section of smooth terrain i1 (Fig. 1), whereas simple flow of material downhill would stop close to the foot of the hill. The ice particles at the end of their trajectory have a horizontal velocity component and this low velocity ballistic sedimentation would lead to formation of lineaments on the smooth terrain, like in solid-particulate volcanic eruptions.     

Ice grain size and the rheology of the martian polar deposits

The history and dynamics of the martian polar deposits (MPD), the largest known water reservoirs on Mars, are of great interest, but estimates of ice grain size are required before detailed modeling can be performed. We clarify the microphysical processes that may control grain size in the MPD. If the MPD are ∼2% dust by mass, the maximum ice grain size is ∼1 mm due to grain boundary pinning by silicate microparticles. Relatively dusty layers in the MPD will have smaller grain sizes. If MPD ice has a very low impurity content and has experienced a significant amount of strain, grains may reach a steady state size of ∼1.5 to 3 mm due to dynamic recrystallization, wherein a steady state grain size is maintained due to the balance of grain growth and destruction during flow. If the near-bed ice in the MPD is warmed close to its melting point and has been extensively sheared, grain sizes at its base may be between 10 and 40 mm, by analogy with warm, dirty, near-bed ice in terrestrial ice sheets.     

Comet Hale-Bopp in outburst: Imaging the dynamics of icy particles with HST/NICMOS

Comet Hale-Bopp was imaged at wavelengths from 1.87 to 2.22 μm by HST/NICMOS in post-perihelion observations starting on UT 1997 August 27.95. Diffraction-limited (∼02) images were obtained at high signal-to-noise (∼1500) to probe the composition and dynamics of the inner coma and also the size and activity of the nucleus. The velocities of several unusual morphological features over a 1.7 h period, indicate that a significant outburst occurred 7.4 h prior to these images while the comet was at a heliocentric distance of 2.49 AU. Similar features are also apparent after re-analysis of pre-perihelion ground-based images. The inner coma (radius ?2500 km) is dominated by an “arc” feature, which expanded and became more diffuse with time. This feature can be modeled as the bright central portion of a “jet of outburst” from a near-equatorial region of the nucleus. Less prominent, time-variable linear and circular morphologies are also apparent. The expansion rates of both the arc feature and the circular morphologies imply a common origin and also suggest a grain size distribution with two broad maxima. In addition, several static linear features extend to the edge of the field of view (21,100 km). Radial brightness profiles are highly asymmetric and only approach a ρ⁻¹ decline at distances ?15,000 km. Images in a narrow-band filter at 2.04 μm exhibit a ∼4% absorption feature relative to nearly simultaneous images at wavelengths of 2.22, 1.90, and 1.87 μm. This absorption is attributed to H₂O ice in the coma grains. The spatial distribution and expansion velocity of the absorption at 2.04 μm indicate that these grains are associated with the outburst. The constancy of the absorption feature indicates no appreciable sublimation over 1.7 h. The unresolved nucleus has a flux density consistent with a 40±10 km diameter assuming a 4% geometric albedo.     

Radar observations of Comet P/2005 JQ5 (Catalina)

Arecibo radar observations of Comet P/2005 JQ5 (Catalina) have produced the first delay-Doppler images of a comet nucleus and the first radar detection of large-grain ejection from a Jupiter-family comet. The nucleus is small (1.4 km diameter), rough, and rapidly rotating. The large (>cm) grains have low velocities (∼1 m/s) and a low production rate.     

Secular light curves of comets, II: 133P/Elst-Pizarro, an asteroidal belt comet

We present the secular light curve (SLC) of 133P/Elst-Pizarro, and show ample and sufficient evidence to conclude that it is evolving into a dormant phase. The SLC provides a great deal of information to characterize the object, the most important being that it exhibits outburst-like activity without a corresponding detectable coma. 133P will return to perihelion in July of 2007 when some of our findings may be corroborated. The most significant findings of this investigation are: (1) We have compiled from 127 literature references, extensive databases of visual colors (37 comets), rotational periods and peak-to-valley amplitudes (64 comets). 2-Dimensional plots are created from these databases, which show that comets do not lie on a linear trend but in well defined areas of these phase spaces. When 133P is plotted in the above diagrams, its location is entirely compatible with those of comets. (2) A positive correlation is found between cometary rotational periods and diameters. One possible interpretation suggest the existence of rotational evolution predicted by several theoretical models. (3) A plot of the historical evolution of cometary nuclei density estimates shows no trend with time, suggesting that perhaps a consensus is being reached. We also find a mean bulk density for comets of 〈ρ〉=0.52±0.06 g/cm³. This value includes the recently determined spacecraft density of Comet 9P/Tempel 1, derived by the Deep Impact team. (4) We have derived values for over 18 physical parameters, listed in the SLC plots, Figs. 6-9. (5) The secular light curve of 133P/Elst-Pizarro exhibits a single outburst starting at +42±4 d (after perihelion), peaking at LAG=+155±10 d, duration 191±11 d, and amplitude 2.3±0.2 mag. These properties are compatible with those of other low activity comets. (6) To explain the large time delay in maximum brightness, LAG, two hypothesis are advanced: (a) the existence of a deep ice layer that the thermal wave has to reach before sublimation is possible, or (b) the existence of a sharp polar active region pointing to the Sun at time = LAG, that may take the form of a polar ice cap, a polar fissure or even a polar crater. The diameter of this zone is calculated at ∼1.8 km. (7) A new time-age is defined and it its found that T-AGE = 80 cy for 133P, a moderately old comet. (8) We propose that the object has its origin in the main belt of asteroids, thus being an asteroid-comet hybrid transition object, an asteroidal belt comet (ABC), proven by its large density. (9) Concerning the final evolutionary state of this object, to be a truly extinct comet the radius must be less than the thermal wave depth, which at 1 AU is ∼250 m (at the perihelion distance of 133P the thermal wave penetrates only ∼130 m). Comets with radius larger than this value cannot become extinct but dormant. Thus we conclude that 133P cannot evolve into a truly extinct comet because it has too large a diameter. Instead it is shown to be entering a dormant phase. (10) We predict the existence of truly extinct comets in the main belt of asteroids (MBA) beginning at absolute magnitude ∼21.5 (diameter smaller than ∼190 m). (11) The object demonstrates that a comet may have an outburst of ∼2.3 mag, and not show any detectable coma. (12) Departure from a photometric R⁺² law is a more sensitive method (by a factor of 10) to detect activity than star profile fitting or spectroscopy. (13) Sufficient evidence is presented to conclude that 133P is the first member of a new class of objects, an old asteroidal belt comet, ABC, entering a dormant phase.     

Crystalline silicates in comets: How did they form?

Two processes have been proposed to explain observations of crystalline silicate minerals in comets and in protostellar sources, both of which rely on the thermal annealing of amorphous grains. First, high temperatures generated by nebular shock processes can rapidly produce crystalline magnesium silicate grains and will simultaneously produce a population of crystalline iron silicates whose average grain size is ∼10-15% that of the magnesium silicate minerals. Second, exposure of amorphous silicate grains to hot nebular environments can produce crystalline magnesium silicates that might then be transported outward to regions of comet formation. At the higher temperatures required for annealing amorphous iron silicates to crystallinity the evaporative lifetime of the grains is much shorter than a single orbital period where such temperatures are found in the nebula. Thermal annealing is therefore unable to produce crystalline iron silicate grains for inclusion into comets unless such grains are very quickly transported away from the hot inner nebula. It follows that observation of pure crystalline magnesium silicate minerals in comets or protostars is a direct measure of the importance of simple thermal annealing of grains in the innermost regions of protostellar nebulae followed by dust and gas transport to the outer nebula. The presence of crystalline iron silicates would signal the action of transient processes such as shock heating that can produce crystalline iron, magnesium and mixed iron-magnesium silicate minerals. These different scenarios result in very different predictions for the organic content of protostellar systems.     

Trapping of N2, CO and Ar in amorphous ice—Application to comets

Recent attempts using high resolution spectra to detect N⁺₂ in several comets were unsuccessful [Cochran, A.L., Cochran, W.D., Baker, E.S., 2000. Icarus 146, 583-593; Cochran, A.L., 2002. Astrophys. J. 576, L165-L168]. The upper limits on N⁺₂ in comparison with the positively detected CO⁺ for Comets C/1995 O1 Hale-Bopp, 122P/1995 S1 de Vico and 153P/2002 C1 Ikeya-Zhang range between . Ar was not detected in three recent comets [Weaver, H.A., Feldman, P.D., Combi, M.R., Krasnopolsky, V., Lisse, C.M., Shemansky, D.E., 2002. Astrophys. J. 576, L95-L98], with upper limits of Ar/CO<(3.4-7.8)×¹⁰⁻² for Comets C/1999 T1 McNaught-Hartley, C/2001 A2 LINEAR and C/2000 WM1 LINEAR. The Ar detected by Stern et al. [Stern, S.A., Slater, D.C., Festou, M.C., Parker, J.Wm., Gladstone, G.R., A'Hearn, M.F., Wilkinson, E., 2000. Astrophys. J. 544, L169-L172] for Comet C/1995 O1 Hale-Bopp, gives a ratio Ar/CO=7.25×¹⁰⁻², which was not confirmed by Cosmovici et al. [Cosmovici, C.B., Bratina, V., Schwarz, G., Tozzi, G., Mumma, M.J., Stalio, R., 2006. Astrophys. Space Sci. 301, 135-143]. Trying to solve the two problems, we studied experimentally the trapping of N₂+CO+Ar in amorphous water ice, at 24-30 K. CO was found to be trapped in the ice 20-70 times more efficiently than N₂ and with the same efficiency as Ar. The resulting Ar/CO ratio of 1.2×¹⁰⁻² is consistent with Weaver et al.'s [Weaver, H.A., Feldman, P.D., Combi, M.R., Krasnopolsky, V., Lisse, C.M., Shemansky, D.E., 2002. Astrophys. J. 576, L95-L98] non-detection of Ar. However, with an extreme starting value for N₂/CO = 0.22 in the region where the ice grains which agglomerated to produce comet nuclei were formed, the expected N₂/CO ratio in the cometary ice should be 6.6×¹⁰⁻³, much higher than its non-detection limit.     

CO clathrate hydrate: Near to mid-IR spectroscopic signatures

Carbon monoxide is the second most abundant molecule after H₂ in the molecular universe, and as such an abundant constituent of interstellar and Solar System ices. To trace the possibility of this molecule to be found in a clathrate hydrate inclusion compound, its pure phase FTIR spectrum is investigated. We confirm the formation of a type I clathrate structure whereas simple guest size estimates would favour a type II clathrate hydrate, revealing interactions of this molecule with its water network during clathrate formation. The observed cage vibrational downshift with respect to pure CO ice is within 5 cm⁻¹. The temperature dependent wavenumber separation between the two enclathrated CO vibrational transitions in the two distinct type I clathrate cages is less than a wavenumber below 140 K, implying that the spectral simplification for detailed spectroscopic analysis of the individual profiles is a difficult task. The dynamics of the CO molecules in its cage change considerably from 5 K to 140 K. At temperatures above 30 K, the molecule is extremely mobile in the cages, as revealed by the infrared profile, significantly different from CO entrapped in water ice and different from observed profiles in astrophysical objects.