The nucleus of Deep Impact target Comet 9P/Tempel 1

On UT 2000 August 21 we obtained simultaneous visible and mid-infrared observations of Comet 9P/Tempel 1, the target of the upcoming NASA Discovery Program mission Deep Impact. The comet was still quite active while 2.55 AU from the Sun (post-perihelion). Two independent analyses of our data, one parameterizing the coma morphology and the other modeling infrared spectrophotometry, show that the nucleus's cross section at the time the data were taken corresponds to an effective radius of 3.0±0.2 km. Based on visible-wavelength photometry of the comet taken during this observing run and others in the summer of 2000, all of which show the rotational modulation of the nucleus's brightness, we find that the infrared data were obtained near the maximum of the light curve. If we assume that the nucleus's light curve had a peak-to-valley range of 0.6±0.2 mag, then the mean effective radius is 2.6±0.2 km. Visible-wavelength photometry of the nucleus, including data published by other groups, lets us constrain the nucleus's R-band geometric albedo: 0.072±0.016. The nucleus's flux contributed about 85% of the light in the mid-infrared images.     

Non-LTE radiative transfer for sub-millimeter water lines in Comet 67P/Churyumov-Gerasimenko

The European Space Agency (ESA) Rosetta spacecraft (Schulz, R., Alexander, C., Boehnhardt, H., Glassmeier, K.H. (Eds.) [2009]. “ROSETTA - ESA”) will encounter Comet 67P/Churyumov-Gerasimenko in 2014 and spend the next 18 months in the vicinity of the comet, permitting very high spatial and spectral resolution observations of the coma and nucleus. During this time, the heliocentric distance of the comet will change from ∼3.5 AU to ∼1.3 AU, accompanied by an increasing temperature of the nucleus and the development of the coma. The Microwave Instrument for the Rosetta Orbiter (MIRO) will observe the ground-state rotational transition (1₁₀-1₀₁) of H₂¹⁶O at 556.936 GHz, the two isotopologues H₂¹⁷O and H₂¹⁸O and other molecular transitions in the coma during this time (Gulkis, S. et al., [2007]. MIRO: Microwave Instrument for Rosetta Orbiter. Space Sci. Rev. 128, 561-597).The aim of this study is to simulate the water line spectra that could be obtained with the MIRO instrument and to understand how the observed line spectra with various viewing geometries can be used to study the physical conditions of the coma and the water excitation processes throughout the coma. We applied an accelerated Monte Carlo method to compute the excitations of the seven lowest rotational levels (1₀₁, 1₁₀, 2₁₂, 2₂₁, 3₀₃, 3₁₂, and 3₂₁) of ortho-water using a comet model with spherically symmetric water outgassing, density, temperature and expansion velocity at three different heliocentric distances 1.3 AU, 2.5 AU, and 3.5 AU. Mechanisms for the water excitation include water-water collisions, water-electron collisions, and infrared pumping by solar radiation.Synthetic line spectra are calculated at various observational locations and directions using the MIRO instrument parameters. We show that observations at varying viewing distances from the nucleus and directions have the potential to give diagnostic information on the continuum temperature and water outgassing rates at the surface of the nucleus, and the gas density, expansion velocity, and temperature of the coma as a function of distance from the nucleus. The gas expansion velocity and temperature affect the spectral line width and frequency shift of the line from the rest frequency, while the gas density (which is directly related to the outgassing rate) and the line excitation temperature determine the antenna temperature of the absorption and emission signal in the line profile.     

Chandra observations of Comet 9P/Tempel 1 during the Deep Impact campaign

We present results from the Chandra X-ray Observatory's extensive campaign studying Comet 9P/Tempel 1 (T1) in support of NASA's Deep Impact (DI) mission. T1 was observed for ∼295 ks between 30th June and 24th July 2005, and continuously for ∼64 ks on July 4th during the impact event. X-ray emission qualitatively similar to that observed for the collisionally thin Comet 2P/Encke system [Lisse, C.M., Christian, D.J., Dennerl, K., Wolk, S.J., Bodewits, D., Hoekstra, R., Combi, M.R., Mäkinen, T., Dryer, M., Fry, C.D., Weaver, H., 2005b. Astrophys. J. 635 (2005) 1329-1347] was found, with emission morphology centered on the nucleus and emission lines due to C, N, O, and Ne solar wind minor ions. The comet was relatively faint on July 4th, and the total increase in X-ray flux due to the Deep Impact event was small, ∼20% of the immediate pre-impact value, consistent with estimates that the total coma neutral gas release due to the impact was 5×¹⁰⁶ kg (∼10 h of normal emission). No obvious prompt X-ray flash due to the impact was seen. Extension of the emission in the direction of outflow of the ejecta was observed, suggesting the presence of continued outgassing of this material. Variable spectral features due to changing solar wind flux densities and charge states were clearly seen. Two peaks, much stronger than the man-made increase due to Deep Impact, were found in the observed X-rays on June 30th and July 8th, 2005, and are coincident with increases in the solar wind flux arriving at the comet. Modeling of the Chandra data using observed gas production rates and ACE solar wind ion fluxes with a CXE mechanism for the emission is consistent, overall, with the temporal and spectral behavior expected for a slow, hot wind typical of low latitude emission from the solar corona interacting with the comet's neutral coma, with intermittent impulsive events due to solar flares and coronal mass ejections.     

SOHO/SWAN observations of comets with small perihelia: C/2002 V1 (NEAT), C/2002 X5 (Kudo–Fujikawa), 2006 P1 (McNaught) and 96P/Machholz 1

SWAN, the all-sky hydrogen Lyman-alpha camera on the SOHO spacecraft, designed primarily to image the interplanetary neutral hydrogen around the Sun, also observes comets continuously over large portions of their apparitions to the north and south of the ecliptic and at small solar elongation angles. Because of SOHO’s location at the L1 Lagrange point, analysis of SWAN images provides excellent temporal coverage of water production. We report here our results of observations of some interesting target comets selected from the extensive SWAN archive. These include three Oort Cloud Comets C/2002 V1 (NEAT), C/2002 X5 (Kudo–Fujikawa), C/2006 P1 (McNaught) and three apparitions of atypical short-period Comet 96P/Machholz 1. The common aspect of these four comets is their small perihelion distances, which are 0.19, 0.09, 0.17, and 0.12 AU, respectively. Their water production rates over their whole apparitions can be approximated by power laws in heliocentric distance (r in AU) as follows: 1.3 × 10²⁹ r^−2.1 s⁻¹ for C/2002 V1 (NEAT), 7.5 × 10²⁸ r^−2.0 s⁻¹ for C/2002 X5 (Kudo–Fujikawa), 5.4 × 10²⁹ r^−2.4 s⁻¹ for C/2006 (P1 McNaught) and 4.6 × 10²⁷ r^−2.1 s⁻¹ for 96P/Machholz 1. We also present daily-average water production rates for the long-period comets over long continuous time periods. We examine these results in light of our growing survey of comets that is yielding some interesting comparisons of water production rate variations with heliocentric distance and taxonomic classes.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

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.     

Comet classification with new methods for gas and dust spectroscopy

We present the results of a program of comet long-slit spectroscopy with the Kast Dual Spectrograph on the 3-m Shane Telescope at Lick Observatory. A total of 26 comets, from a variety of dynamical families, were observed on 39 different nights from 1996 to 2007. A new statistical method extracted the twilight sky from comet frames, because traditional sky subtraction techniques were inadequate. Because previously published Haser model parent and daughter scale lengths did not fit the data well, unbiased ranges of scale lengths were searched for the best-fitting pairs. Coma gas production rates for OH, CN, C₂, C₃, NH, NH₂, and OH confirmed the widely reported carbon-chain depletion for a sub-class of comets, most notably high-perihelion Jupiter-family comets observed at r_h > 1.5 AU, with different behaviors for C₂ and C₃. Our long-slit spectroscopy data was also adapted for the A(θ)fρ dust production parameter. The assumption that A(θ)fρ is constant throughout the nucleus was not upheld. High dust-to-gas ratios for comets with large perihelia were not a selection effect, and suggest that the dust was released earlier in the formation of the coma than the gas. The dust-to-gas ratio did not exhibit any evolutionary traces between different comet dynamical families. The comet survey illuminates the diversity among comets, including the unusually carbon poor Comet 96P/Machholz.     

Thermal model of water and CO activity of Comet C/1995 O1 (Hale-Bopp)

An investigation of the activity of Comet C/1995 O1 (Hale-Bopp) with a thermophysical nucleus model that does not rely on the existence of amorphous ice is presented. Our approach incorporates recent observations allowing to constrain important parameters that control cometary activity. The model accounts for heat conduction, heat advection, gas diffusion, sublimation, and condensation in a porous ice-dust matrix with moving boundaries. Erosion due to surface sublimation of water ice leads to a moving boundary. The movement of the boundary is modeled by applying a temperature remapping technique which allows us to account for the loss in the internal energy of the eroded surface material. These kind of problems are commonly referred to as Stefan problems. The model takes into account the diurnal rotation of the nucleus and seasonal effects due to the strong obliquity of Hale-Bopp as reported by Jorda et al. (Jorda, L., Rembor, K., Lecacheux, J., Colom, P., Colas, F., Frappa, E., Lara, L.M. [1997]. Earth Moon Planets 77, 167-180). Only bulk sublimation of water and CO ice are considered without further assumptions such as amorphous ices with certain amount of occluded CO gas. Confined and localized activity patterns are investigated following the reports of Lederer and Campins (Lederer, S.M., Campins, H. [2002]. Earth Moon Planets 90, 381-389) about the chemical heterogeneity of Hale-Bopp and of Bockelée-Morvan et al. (Bockelée-Morvan, D., Henry, F., Biver, N., Boissier, J., Colom, P., Crovisier, J., Despois, D., Moreno, R., Wink, J. [2009]. Astron. Astrophys. 505, 825-843) about a strong CO source at a latitude of 20°. The best fit to the observations of Biver et al. (Biver, N. et al. [2002]. Earth Moon Planets 90, 5-14) is obtained with a low thermal conductivity of 0.01 W m⁻¹ K⁻¹. This is in agreement with recent results of the Deep Impact mission to 9P/Tempel 1 (Groussin, O., A’Hearn, M.F., Li, J.-Y., Thomas, P.C., Sunshine, J.M., Lisse, C.M., Meech, K.J., Farnham, T.L., Feaga, L.M., Delamere, W.A. [2007]. Icarus 187, 16-25) and with previous thermal simulations (Kührt, E. [1999]. Space Sci. Rev. 90, 75-82). The water production curve matches the production rates well from −4 AU pre-perihelion to the outgoing leg while the model does not reproduce so well the water production beyond 4 AU pre-perihelion. The CO production curve is a good fit to the measurements of Biver et al. (2002) over the whole measured heliocentric range from −7 AU pre- to 15 AU post-perihelion.     

The outgassing and composition of Comet 19P/Borrelly from radio observations

Radio spectroscopic observations of Comet 19P/Borrelly were performed during the 1994 apparition and at, and near, the time of the Deep Space 1 flyby in 2001. HCN, CS, CH₃OH, and H₂CO were detected using the 30-m telescope of the Institut de Radioastronomie Millimétrique and the James Clerk Maxwell Telescope, and their production rates relative to water are estimated to be 0.06-0.11, 0.07, 1.7, and 0.4%, respectively. Only upper limits are derived for H₂S and CO. The upper limit for CO/H₂O (<15%) is not very constraining, while the upper limit for the H₂S/H₂O ratio of 0.45% is near the bottom of the range of values measured for other comets. Observations of the OH radical at the Nançay radio telescope provide water production rates a few weeks before the 1994 and 2001 perihelia. Observations of the 1₁₀-1₀₁ water line at 557 GHz with the Odin satellite yield a water production rate of (2.5±0.5)×10²⁸ s⁻¹ on September 22, 2001, at the time of the Deep Space 1 encounter, and (3.3±0.6)×10²⁸ s⁻¹ averaged over the September 22-24, 2001 period. The line shapes are asymmetric and blueshifted by V₀∼−0.18 km s⁻¹ for the best observed HCN lines recorded one week after perihelion. The HCN line shapes, and the similar OH and HCN velocity shifts over the September-November 1994 and August-September 2001 periods, favor anisotropic outgassing towards the Sun. Strong outgassing directed along the primary dust jet seen on visible images is not excluded by the HCN line shapes, but unrealistically high gas expansion velocities are required to explain the line shapes in that case.     

Small-scale dust structures in Halley's coma: II. Disintegration of large dust bodies

Small-scale dust structures, SDSs, altogether ∼35 events with extent ∼30-220 km, have been recognized owing to electric field records, mostly near the closest approach of Vega-2 to Halley's nucleus. Several (8-9) morphological forms of SDS have been identified, and all they make one family. Among the family members, the key form (with respect to which, all other forms can be regarded as degenerate) is a sequence of 3-5 dust clouds. The morphological forms represent various Vega-2 passes through SDSs at different stages of development. SDSs observable as the key form consisted of several fairly regularly spaced dust subpopulations, whose plane of symmetry was parallel to the comet orbit plane. That regularity together with specific features of morphological forms strongly constrain disintegration scenarios and dynamics of fragments, and allow to draw a number of conclusions, the main of which are: SDS parent bodies were ice-free dust aggregates lifted from the nucleus near the comet perihelion, whose masses were in the range ∼0.1-1 of the biggest emitted mass (mass of a body accelerated to the escape velocity, i.e., ∼300-1500 kg); the disintegration scenario comprised a few steps, and the first-step disintegration consisted mainly in consecutive detachments of biggest first-step fragments (BF-SFs) from the parent body; a SDS observable as the key form included the dust minitail of parent body and a few BF-SF minitails, the former one being longer than the latter ones; SDS parent bodies had a fractal-like internal structure, and the BF-SF mass was a few percent of the parent body mass; the thermal conductivity of SDS parent body was less than ∼0.4 W m⁻¹ K⁻¹ or so, while the latent heat of gluing organics was roughly 80 kJ mol⁻¹; the disintegration mechanism was a combination of sintering and sublimation of organics. The multistep disintegration of SDS parent bodies can be reconciled with the basically one-step disintegration of aggregates responsible for the dust boundary (Oberc, P., Icarus 1996, 124, 195-208). The fractal-like structure and the relation between BF-SF mass and parent body mass are in agreement with predictions from the Weidenschilling model of comet formation. Large ice-free dust bodies, in particular SDS parent bodies, can be identified with refractory boulders postulated by some comet nucleus models.     

Secular light curve of Comet 9P/Tempel 1

In support of the Deep Impact Mission, we have updated the secular light curve of 9P/Tempel 1 presented in Paper I [Ferrín, I., 2005. Icarus 178, 493-516], with new data sets. The secular light curves (SLC) of the comet are presented in the log and time plots (Figs. 1 and 2) and provide a clear profile of the overall shape of the envelope. We arrive at the following conclusions: (1) Improved values of 18 photometric parameters are derived including the turn on and turn off points, RON=−3.47±0.05 AU, ROFF=+4.20±0.05 AU, and TON=−410±25 d, TOFF=+555±25 d. (2) The improved SLC shows a most interesting and peculiar shape, with a linear power law of slope n=7.7±0.1 from RON=−3.47 AU to RBP=−2.08±0.05 AU, and then converts to a law with curvature. The break point of the power law at RBP=−2.08 AU, mV(1,R)=14.0±0.1 mag, is interpreted as a change in sublimating something more volatile than water ice (most probably CO₂), to water ice sublimation. In other words, the comet's sublimation is controlled by two different substances. (3) The photometric-age (defined in Paper I) and the time-age of the comet [Ferrín, I., 2006. Icarus. In press] are recomputed, and results in a value P-AGE=21±2 and T-AGE=11±2 comet years. Thus 9P is a young comet. (4) The comet is active almost up to aphelion since the turn off point has been determined at ROFF=+4.20±0.05 AU while aphelion takes place at Q=+4.74 AU. (5) The comet exhibits activity post-aphelion which is not understood. Two hypothesis are advanced to explain this behavior.     

Investigations of the pre-Deep Impact morphology of the dust coma of Comet Tempel-1

The pre-Deep Impact images of Comet Tempel-1 obtained at the Indian Astronomical Observatory are used to investigate the morphology of the dust coma of the comet. We show that the trajectory of a cometary grain under the influence of solar radiation pressure is a reliable diagnostic to estimate its initial velocity. Four main active regions at mean latitudes +45° ± 5°(D), 0° ± 5° (E),−30° ± 5°(A) and−60° ± 5°(F) are found to explain the morphology of the dust coma in the ground-based and published images obtained by the High Resolution Instrument(HRI) cameras aboard the Deep Impact flyby spacecraft. From a χ² fit of the intensity distribution in the observed and the simulated images, we derive the fraction of the productivity of the active vents to the total dust emission of the comet to be 27%. Of this the southern source alone accounts for 19.8%. The grains are found to be ejected with a velocity distribution with an upper limit of 70 ± 7 m s⁻¹. However, the broad region ‘A’ appears to eject slower grains with an upper limit of 24 ± 2.5 m s⁻¹. This source, that is active throughout the cycle is likely to be driven by CO₂ sublimation. We compute the dependence of the percentage contribution of the southern source on the heliocentric distance and show that this ratio varies over the apparition and reaches a maximum at around 260 days before perihelion. The published images of the nucleus of Comet Tempel-1 show significant departure from sphericity. Therefore, the torque exerted by the enhanced activity of the southern region may be significant enough to produce changes in the rotational state of the nucleus before each perihelion passage.     

Compositional and physical results for Rosetta's new target Comet 67P/Churyumov-Gerasimenko from narrowband photometry and imaging

We present compositional and physical results of Comet 67P/Churyumov-Gerasimenko, the new target of ESA's Rosetta mission. A total of 16 nights of narrowband photometry were obtained at Lowell Observatory during the 1982/83 and 1995/96 apparitions, along with one night of imaging near perihelion in 1996. These data encompass an interval of −61 to +118 days from perihelion, corresponding to a range of heliocentric distances before perihelion from 1.48 to 1.34 AU, and an outbound range from 1.30 to 1.86 AU. Production rates were determined for OH, NH, CN, C₃, and C₂, along with A(θ)fρ, a proxy of the dust production. Water production, based on OH, has a steep () power-law rH-dependence post-perihelion and the minor species are somewhat less steep ( to −4), while the dust is quite shallow (), possibly due to a lingering population of large, slow-moving grains. All species exhibit larger production rates after perihelion, with water having a ∼2×pre/post-perihelion asymmetry, while minor species and dust have larger asymmetries. These asymmetries imply a strong seasonal effect and probable high obliquity of the rotational axis, along with one or more isolated source regions coming into sunlight near perihelion. Peak water production (which occurred about 1 month after perihelion) was and, when combined with a standard water vaporization model, implies an effective active area on the surface of the nucleus of ∼1.5-2.2 km² or an active fraction of only about 3-4%. Abundances of carbon-chain molecules yield a classification of slightly “depleted” in the A'Hearn et al. [A'Hearn, M.F., Millis, R.L., Schleicher, D.G., Osip, D.J., Birch, P.V., 1995. Icarus 118, 223-270] database. The peak dust production (as measured by A(θ)fρ, and uncorrected for phase angle) was ∼450 cm, while the color of the dust is moderately reddened, and the mean radial profile has a power-law slope of −1.3. Large night-to-night variability is also present, presumably due to the source region(s) rotating in and out of sunlight along with effects due to the use of differently sized apertures. A strong sunward radial feature was detected in images obtained near perihelion, along with a significant asymmetry between the two perpendicular directions from the Sun/tail line. These features may be the result of a mid-latitude source region sweeping out a cone with each rotation, which we are viewing from the side and where the sunward radial feature is one edge of the cone seen in projection. When combined with other constraints on the pole orientation, a possible pole solution is found having an obliquity of about 134° at an RA of about 223° and a Dec of −65°, with a source region located near +50° and in overall agreement with the photometric results. In comparison to the original Rosetta target Comet 46P/Wirtanen, Comet Churyumov-Gerasimenko has essentially the same peak water production but a peak dust production about 3 times greater than does Wirtanen based on A(θ)fρ (i.e., if one assumes that the properties of the dust grains are similar) (cf. Farnham and Schleicher [1998. Astron. Astrophys. 335, L50-L55]).     

Physical properties of morphological units on Comet 9P/Tempel 1 derived from near-IR Deep Impact spectra

In this paper we analyze near-infrared thermal emission spectra of the spatially resolved nucleus of Comet 9P/Tempel 1 obtained by the NASA spacecraft Deep Impact. Maps of spectral reddening, the product X^′ between the beaming function and directional emissivity, as well as surface temperature are constructed. Thermophysical modeling is used to estimate the degree of small scale surface roughness and thermal inertia by detailed reproduction of the empirical temperature map. Mie and Hapke theories are used in combination with numerically calculated beaming functions to analyze the X^′ map and place constraints on composition and grain size of the surface material. We show that it is absolutely mandatory to include small scale surface roughness in thermophysical modeling of this object, since the resulting self heating is vital for reproducing the measured temperatures. A small scale self heating parameter in the range 0.6?ξ?0.75 is common, but smoother areas where 0.2?ξ?0.3 are also found. Contrary to models neglecting small scale surface roughness, we find that the thermal inertia of Comet 9P/Tempel 1 generally is high (1000-3000 J m⁻² K⁻¹ s^−1/2), although it may be substantially lower (40-380 J m⁻² K⁻¹ s^−1/2) in specific areas. We obtain a disk-averaged reddening of 3.5% kÅ⁻¹, with statistically significant local variations around that value on a ±1.0% kÅ−1 level. Vast regions appear covered by small (∼0.1 μm) highly absorbing grains such as carbon or iron-rich silicates. Other regions appear dominated by somewhat larger (∼0.5 μm) and/or less absorbing grains such as troilite or magnesium-rich silicates. Surface variations in reddening, roughness, thermal inertia, composition and/or grain size are moderately to strongly correlated to the locations of morphological units on the surface. The existence of morphological units with differing physical properties may be primordial, hence reflecting a diversity in the building block cometesimals, or resulting from evolutionary processes.     

Explosion of Comet 17P/Holmes as revealed by the Spitzer Space Telescope

An explosion on Comet 17P/Holmes occurred on 2007 October 23, projecting particulate debris of a wide range of sizes into the interplanetary medium. We observed the comet using the mid-Infrared Spectrograph (5-40 μm), on 2007 November 10 and 2008 February 27, and the imaging photometer (24 and 70 μm), on 2008 March 13, on board the Spitzer Space Telescope. The 2007 November 10 spectral mapping revealed spatially diffuse emission with detailed mineralogical features, primarily from small crystalline olivine grains. The 2008 February 27 spectra, and the central core of the 2007 November 10 spectral map, reveal nearly featureless spectra, due to much larger grains that were ejected from the nucleus more slowly. Optical images were obtained on multiple dates spanning 2007 October 27-2008 March 10 at the Holloway Comet Observatory and 1.5-m telescope at Palomar Observatory. The images and spectra can be segmented into three components: (1) a hemispherical shell fully 28′ on the sky in 2008 March, due to the fastest (262 m s⁻¹), smallest (2 μm) debris, with a mass ; (2) a ‘blob’ or ‘pseudonucleus’ offset from the true nucleus and subtending some 10′ on the sky, due to intermediate speed (93 m s⁻¹) and size (8 μm) particles, with a total mass ; and (3) a ‘core’ centered on the nucleus due to slower (9 m s⁻¹), larger (200 μm) ejecta, with a total mass . This decomposition of the mid-infrared observations can also explain the temporal evolution of the millimeter-wave flux. The orientation of the leading edge of the ejecta shell and the ejecta ‘blob,’ relative to the nucleus, do not change as the orientation of the Sun changes; instead, the configuration was imprinted by the orientation of the initial explosion. The distribution and speed of ejecta implies an explosion in a conical pattern directed approximately in the solar direction on the date of explosion. The kinetic energy of the ejecta >10²¹ erg is greater than the gravitational binding energy of the nucleus. We model the explosion as being due to crystallization and release of volatiles from interior amorphous ice within a subsurface cavity; once the pressure in the cavity exceeded the surface strength, the material above the cavity was propelled from the comet. The size of the cavity and the tensile strength of the upper layer of the nucleus are constrained by the observed properties of the ejecta; tensile strengths on >10 m scale must be greater than 10 kPa (or else the ejecta energy exceeds the binding energy of the nucleus) and they are plausibly 200 kPa. The appearance of the 2007 outburst is similar to that witnessed in 1892, but the 1892 explosion was less energetic by a factor of about 20.     

Submillimeter Wave Astronomy Satellite observations of Comet 9P/Tempel 1 and Deep Impact

On 4 July 2005 at 5:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the ₁₁₀-₁₀₁ ortho-water ground-state rotational transition in Comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS and we derive an upper limit of 1.8×¹⁰⁷ kg on the water ice evaporated by the impact. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before and after the impact. Episodes of increased activity with alternated with periods with low outgassing (). We estimate that 9P/Tempel 1 vaporized a total of N∼4.5×¹⁰³⁴ water molecules (∼1.3×¹⁰⁹ kg) during June-September 2005. Our observations indicate that only a small fraction of the nucleus of Tempel 1 appears to be covered with active areas. Water vapor is expected to emanate predominantly from topographic features periodically facing the Sun as the comet rotates. We calculate that appreciable asymmetries of these features could lead to a spin-down or spin-up of the nucleus at observable rates.     

Energy balance of the Deep Impact experiment

We present results on the energy balance of the Deep Impact experiment based on analysis of 180 infrared spectra of the ejecta obtained by the Deep Impact spacecraft. We derive an output energy of 16.5 (+9.1/−4.1) GJ. With an input energy of 19.7 GJ, the error bars are large enough so that there may or may not be a balance between the kinetic energy of the impact and that of outflowing materials. Although possible, no other source of energy other than the impactor or the Sun is needed to explain the observations. Most of the energy (85%) goes into the hot plume in the first few seconds, which only represents a very small fraction (<0.01%) of the total ejected mass. The hot plume contains 190 (+263/−71) kg of H₂O, 1.6 ± 0.5 kg of CO₂, 8.2 (+11.3/3.1) kg of CO (assuming a CO/H₂O ratio of 4.3%), 27.9 (+25.0/−8.9) kg of organic material and 255 ± 128 kg of dust, while the ejecta contains ∼10⁷ kg of materials. About 12% of the energy goes into the ejecta (mostly water) and 3% to destroy the impactor. Volatiles species other than H₂O (CO₂, CO or organic molecules) contribute to <7% of the energy balance. In terms of physical processes, 68% of the energy is used to accelerate grains (kinetic energy), 16% to heat them, 6% to sublimate or melt them and 10% (upper limit) to break and compress dust and/or water ice aggregates into small micron size particles. For the hot plume, we derive a dust/H₂O ratio of 1.3 (+1.9/−1.0), a CO₂/H₂O ratio of 0.008 (+0.009/−0.006), an organics/H₂O ratio of 0.15 (+0.29/−0.11) and an organics/dust ratio of 0.11 (+0.30/−0.07). This composition refers to the impact site and is different from that of the bulk nucleus, consistent with the idea of layers of different composition in the nucleus sub-surface. Our results emphasize the importance of laboratory impact experiments to understand the physical processes involved at such a large scale.     

Water production and release in Comet 153P/Ikeya-Zhang (C/2002 C1): accurate rotational temperature retrievals from hot-band lines near 2.9-μm

Multiple non-resonance fluorescence lines of water (H₂O) were detected in Comet 153/P Ikeya-Zhang (2002 C1) between UT 2002 March 21.9 (R_h=0.51 AU) and April 13.9 (R_h=0.78 AU), using the Cryogenic Echelle Spectrometer (CSHELL) at the NASA Infrared Telescope Facility. Analysis of 2.9-μm water lines enabled accurate determination of rotational temperatures on three dates. The derived H₂O rotational temperatures were 138⁺⁶₋₅, 141⁺¹⁰₋₉, and 94±3 K on UT 2002 March 22.0, March 23.0, and April 13.8, respectively. Water production rates were retrieved from spectral lines measured in nineteen separate grating settings over seven observing periods. The derived heliocentric dependence of the water production rate was Q=(9.2±1.1)×10²⁸[R_h^{(−3.21±0.26)}] molecules s⁻¹. The spatial distribution of H₂O in the coma was consistent with its release directly from the nucleus (as a native source) on all dates.     

CN, C2 radicals and dust in Comets Shoemaker-Levy 1991 T2 and P/deVico 1995 S1

The results of the multiaperture photometry of Comet Shoemaker-Levy 1991 T2 in the pre-perihelion and P/deVico in the post-perihelion period with the narrowband CN, C₂ and Blue Continuum (BC) IHW filters are presented. A Haser model of the molecular coma was used for the determination of the parent and daughter scale-lengths and production rates of the radicals. The comets showed some substantial differences between their parent scale-lengths. The CN parent scale-length (at 1.0 AU) was 16×¹⁰³ km for Comet Shoemaker-Levy and 39×¹⁰³ for P/deVico, the C₂ parent scale-lengths were respectively 29×¹⁰³ and 54×¹⁰³ km. Such divergences could be interpreted in the frame of different scenarios of emission of cometary parents, either from a nucleus or from a volume source. The daughter scale-lengths for these comets were quite similar, namely: 306×¹⁰³ and 318×¹⁰³ km for CN and 69×¹⁰³ and 66×¹⁰³ km for C₂. We determined the Afρ parameter for apertures of different radii. A Monte Carlo model of the dust coma was used to obtain the dust ejection velocity. It was of the order of 0.1 km s⁻¹ for both comets. The power index of the distribution of the β-parameter of dust particles (ratio of light pressure to the solar gravitation) was of the order of 3 for C/Shoemaker-Levy and close to 2 for P/deVico. The dependence on heliocentric distance (rh) of the radical and dust production rates for P/deVico in the range of 0.7-1.0 AU was described by the power law function with a power index equal to: 5.55±0.14 for CN, 5.70±0.24 for C₂ and 5.22±0.19 for dust. Relative abundances of the dynamically new Comet Shoemaker-Levy and short-period P/deVico were quite similar with an enhancement of C₂ comparing with standard values taken from A'Hearn et al. (1995).