Evolution of depressions on Comet 67P/Churyumov-Gerasimenko: Role of ice metamorphism

This work is intended to investigate the influence of temperature-dependent metamorphism of ice on the shape of small depressions in the surface of cometary nuclei. We are mainly interested in the role of initial cohesivity of a nucleus. For this purpose we simulate sublimation of ice from the facets of initially cylindrical depressions in ice of different initial structure. The simulations account for the diurnal and orbital changes of insolation and its dependence on the current shape of the depressions. Our model includes heat transport in the cometary material and metamorphism of ice. We present the results obtained for the nucleus of the Comet 67P/Churyumov-Gerasimenko, target of the ESA cornerstone mission Rosetta.     

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 effect of internal inhomogeneity on the activity of comet nuclei - Application to Comet 67P/Churyumov-Gerasimenko

We present thermal evolution calculations of inhomogeneous asymmetric initial configurations of a spherical model of Comet 67P/Churyumov-Gerasimenko, using a fully 3-dimensional numerical code. The initial composition is amorphous H₂O ice and dust, in a “layered-pile” configuration, where layers differing in ice/dust ratio and thermal properties extend over a fraction of the surface area and about 10 m in depth and may overlap. We analyze the effect of one such layer, as well as the combined effect of many layers, randomly distributed. We find that internal inhomogeneities affect both the surface temperature and the activity pattern of the comet. In particular, they may lead to outbursts at large heliocentric distances and also to activity on the night-side of the nucleus. The rates of ablation and depths of dust mantle and crystalline ice outer layer as functions of longitude and latitude are shown to be affected as well.     

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.     

A ∼25 K temperature of formation for the submicron ice grains which formed comets

Following the observations of ice grains in cometary comae and their size distributions, we reexamined experimentally our previous conclusion that the ice grains which agglomerated to form comet nuclei were formed at ∼25 K. The suggestion of a ∼25 K formation temperature was confirmed experimentally. Moreover, we suggest that these ice grains had to be of submicron size.     

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.     

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.     

Comparison of the composition of the Tempel 1 ejecta to the dust in Comet C/Hale-Bopp 1995 O1 and YSO HD 100546

Spitzer Infrared Spectrograph observations of the Deep Impact experiment in July 2005 have created a new paradigm for understanding the infrared spectroscopy of primitive solar nebular (PSN) material—the ejecta spectrum is the most detailed ever observed in cometary material. Here we take the composition model for the material excavated from Comet 9P/Tempel 1's interior and successfully apply it to Infrared Space Observatory spectra of material emitted from Comet C/1995 O1 (Hale-Bopp) and the circumstellar material found around the young stellar object HD 100546. Comparison of our results with analyses of the cometary material returned by the Stardust spacecraft from Comet 81P/Wild 2, the in situ Halley flyby measurements, and the Deep Impact data return provides a fundamental cross-check for the spectral decomposition models presented here. We find similar emission signatures due to silicates, carbonates, phyllosilicates, water ice, amorphous carbon, and sulfides in the two ISO-observed systems but there are significant differences as well. Compared to Tempel 1, no Fe-rich olivines and few crystalline pyroxenes are found in Hale-Bopp and HD 100546. The YSO also lacks amorphous olivine, while being super-rich in amorphous pyroxene. All three systems show substantial emission due to polycyclic aromatic hydrocarbons. The silicate and PAH material in Hale-Bopp is clearly less processed than in Tempel 1, indicating an earlier age of formation for Hale-Bopp. The observed material around HD 100546 is located ∼13 AU from the central source, and demonstrates an unusual composition due to either a very different, non-solar starting mix of silicates or due to disk material processing during formation of the interior disk cavity and planet(s) in the system.     

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.     

The shape, topography, and geology of Tempel 1 from Deep Impact observations

Deep Impact images of the nucleus of Comet Tempel 1 reveal pervasive layering, possible impact craters, flows with smooth upper surfaces, and erosional stripping of material. There are at least 3 layers 50-200 m thick that appear to extend deep into the nucleus, and several layers 1-20 m thick that parallel the surface and are being eroded laterally. Circular depressions show geographical variation in their forms and suggest differences in erosion rates or style over scales >1 km. The stratigraphic arrangement of these features suggests that the comet experienced substantial periods of little erosion. Smooth surfaces trending downslope suggest some form of eruption of materials from this highly porous object. The Deep Impact images show that the nucleus of Tempel 1 cannot be modeled simply as either an onion-layer or rubble pile structure.     

Near-infrared spectroscopy of 133P/Elst-Pizarro

Comet 133P/Elst-Pizarro, known to be a main-belt asteroid with a cometary activity, was observed with the near-infrared integral field spectrograph SINFONI at the 8.2-m Very Large Telescope in J, H, and K bands during its 2007 perihelion passage. The goal of these observations was to attempt detection of water ice absorption bands. We present here the details of the data processing, the results of these observations, and our compositional modeling of the final spectrum. No water ice absorptions were detected within the noise of the spectrum but we show that this spectrum is compatible with a reasonable near-infrared albedo value of 7-10% and a mixture of water ice, black carbon, Tholin and silicates. This interpretation, nevertheless, is not unique.     

Circular polarization in comets: Observations of Comet C/1999 S4 (LINEAR) and tentative interpretation

Comet C/1999 S4 (LINEAR) was exceptional in many respects. Its nucleus underwent multiple fragmentations culminating in the complete disruption around July 20, 2000. We present circular polarization measurements along the cuts through the coma and nucleus of the comet during three separate observing runs, in June 28-July 2, July 8-9, and July 21-22, 2000. The circular polarization was detected at a rather high level, up to 0.8%. The left-handed as well as right-handed polarization was observed over the coma with the left circularly polarized light systematically observed in the sunward part of the coma. During our observations the phase angle of the comet varied from 61 up to 122°, which allowed us to reveal variations of circular polarization with the phase angle. Correlation between the degree of circular polarization, visual magnitude, water production rate, and linear polarization of Comet C/1999 S4 (LINEAR) during its final fragmentation in July 2000 was found. The mechanisms that may produce circular polarization in comets and specifically in Comet C/1999 S4 (LINEAR) are discussed and some tentative interpretation is presented.     

Optical polarimetry of Comet NEAT C/2001 Q4

Comet NEAT C/2001 Q4 was observed for linear polarization using the optical polarimeter mounted at the 1.2 m telescope at Mt. Abu Observatory, during the months of May and June 2004. Observations were conducted through the International Halley Watch narrow band (continuum) and BVR broad band filters. During the observing run the phase angle ranged from 85.6° in May to 55° in June. As expected, polarization increases with wavelength in this phase angle range. Polarization colour in the narrow bands changes at different epochs, perhaps related to cometary activity or molecular emission contamination. The polarization was also measured in the cometary coma at different locations along a line, in the direction of the tail. As expected, we notice minor decrease in the polarization as photocenter (nucleus) is traversed while brightness decreases sharply away from it. Based on these polarization observations we infer that the Comet NEAT C/2001 Q4 has high polarization and a typical grain composition—mixture of silicates and organics.     

A multi-wavelength study of parent volatile abundances in Comet C/2006 M4 (SWAN)

Volatile organic emissions were detected post-perihelion in the long-period Comet C/2006 M4 (SWAN) in October and November 2006. Our study combines target-of-opportunity infrared observations using the Cryogenic Echelle Spectrometer (CSHELL) at the NASA-IRTF 3-m telescope, and millimeter wavelength observations using the Arizona Radio Observatory (ARO) 12-m telescope. Five parent volatiles were measured with CSHELL (H₂O, CO, CH₃OH, CH₄, and C₂H₆), and two additional species (HCN and CS) were measured with the ARO 12-m. These revealed highly depleted CO and somewhat enriched CH₃OH compared with abundances observed in the dominant group of long-period (Oort cloud) comets in our sample and similar to those observed recently in Comet 8P/Tuttle. This may indicate highly efficient H-atom addition to CO at very low temperature (∼10-20 K) on the surfaces of interstellar (pre-cometary) grains. Comet C/2006 M4 had nearly “normal” C₂H₆ and CH₄, suggesting a processing history similar to that experienced by the dominant group. When compared with estimated water production at the time of the millimeter observations, HCN was slightly depleted compared with the normal abundance in comets based on IR observations but was consistent with the majority of values from the millimeter. The ratio CS/HCN in C/2006 M4 was within the range measured in ten comets at millimeter wavelengths. The higher apparent H-atom conversion efficiency compared with most comets may indicate that the icy grains incorporated into C/2006 M4 were exposed to higher H-atom densities, or alternatively to similar densities but for a longer period of time.     

Fully 3-dimensional calculations of dust mantle formation for a model of Comet 67P/Churyumov-Gerasimenko

A fully 3-dimensional implicit numerical model for comet nucleus evolution is presented, emphasizing dust mantle formation. A spherical configuration is considered with an initial composition of amorphous H₂O ice and dust, taking into account a discrete dust-grain size distribution. The model is applied to Comet 67P/Churyumov-Gerasimenko, adopting its orbital elements, rotation period and rotation axis inclination. We find that the dust mantle thickness varies over the surface from 1 cm to about 10 cm (thus lower and higher than the diurnal skin-depth, respectively). The size distribution of ejected grains varies along the orbit and is steeper than the initial one adopted for the nucleus. The crystallization front advances inward in spurts, and its depth varies between 1 and several meters. We test the effect of the thermal conductivity on the surface temperature distribution and depths of the dust mantle and crystallization front.     

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.     

Radar observations of 8P/Tuttle: A contact-binary comet

Arecibo radar imagery of Comet 8P/Tuttle reveals a 10-km-long nucleus with a highly bifurcated shape consistent with a contact binary. A separate echo component was also detected from large (>cm-size), slow-moving grains of the type expected to contribute to the Ursid meteor stream.     

The size and thermal properties of the nucleus of Comet 22P/Kopff

We detected the nucleus of Comet 22P/Kopff at 4.87 AU from the Sun with the two IRS peak-up cameras of the Spitzer Space Telescope on April 19, 2007. Using the thermal model of [Groussin, O., and 15 colleagues, 2007. Icarus 187, 16-25], we derive a nucleus size of 1.89±0.16 km, in agreement with [Lamy, P., Toth, I., Jorda, L., Groussin, O., A'Hearn, M.F., Weaver, H.A., 2002. Icarus 156, 442-455], and a thermal inertia .     

Surface temperature of the nucleus of Comet 9P/Tempel 1

The Deep Impact (DI) spacecraft encountered Comet 9P/Tempel 1 on July 4th, 2005 and observed it with several instruments. In particular, we obtained infrared spectra of the nucleus with the HRI-IR spectrometer in the wavelength range of 1.0-4.9 μm. The data were taken before impact, with a maximum resolution of ∼120 m per pixel at the time of observation. From these spectra, we derived the first directly observed temperature map of a comet nucleus. The surface temperature varied from 272±7 to 336±7 K on the sunlit hemisphere, matching the surface topography and incidence angle. The derived thermal inertia is low, most probably <50 W K⁻¹ m⁻² s^1/2. Combined with other arguments, it is consistent with the idea that most of rapidly varying thermal physical processes, in particular the sublimation of volatiles around perihelion, should occur close to the surface. Thermal inertia is sufficient to explain the temperature map of the nucleus of Comet Tempel 1 to first order, but other physical processes like roughness and self-radiation are required to explain the details of the temperature map. Finally, we evaluated that the Standard Thermal Model is a good approximation to derive the effective radius of a cometary nucleus with an uncertainty lower than ∼10% if combined with a thermal infrared light curve.