Comet 17P/Holmes: Possibility of a CO driven explosion

This work is a continuation of our previous paper about brightening of Comet 17P/Holmes (Kossacki, K.J., Szutowicz, S. [2010]. Icarus 207, 320–340). In that paper we presented results of simulations indicating that the nonuniform crystallization of amorphous water ice itself is probably not sufficient for an explosion. In the present work we investigate the possibility that the explosion is caused by a rapid sublimation of the CO ice leading to the rise of gas pressure above the tensile strength of the nucleus. We simulated evolution of a model nucleus in the orbit of Comet 17P/Holmes. The nucleus is composed of water ice, carbon monoxide ice and dust and has the shape of an elongated ellipsoid. The simulations include crystallization of amorphous ice in the nucleus, changes of the dust mantle thickness, and sublimation of the CO ice. In our model CO is mantling grains composed of dust and amorphous water ice. Orientation of the nuclear spin axis in space is the same as derived in Moreno et al. (Moreno, F., Ortiz, J.L., Santos-Sanz, P., Morales, N., Vidal-Nunez, M.J., Lara, L.M., Gutierrez, P.J. [2008]. Astrophys. J. 677, L63–L66) for Comet Holmes during recent brightening event. Hence, the angle between the orbital and the equatorial planes of the comet is I = 95°, and the cometocentric solar longitude at perihelion is Φ = 210°. The calculations are performed for the south pole being the sub-solar point close to time of the outburst. Our computations indicate, that the CO pressure within the comet nucleus can rise to high values. When the layer between the dust mantle and the crystallization front of the amorphous water ice is very fine grained, few microns in radius, the CO pressure within the nucleus can exceed 10 kPa. This value is the lowest estimate for the tensile strength of the nucleus of Comet Holmes (Reach, W.T., Vaubaillon, J., Lisse, C.M., Holloway, M., Rho, J. [2010]. Icarus 208, 276–292). Hence, when the gas pressure reaches this value the nucleus may explode.     

Follow‐up observations of Comet 17P/Holmes after its extreme outburst in brightness end of October 2007

We present follow‐up observations of comet 17/P Holmes after its extreme outburst in brightness, which occurred end of October 2007. We obtained 58 V‐band images of the comet between October 2007 and February 2008, using the Cassegrain‐Teleskop‐Kamera (CTK) at the University Observatory Jena. We present precise astrometry of the comet, which yields its most recent Keplerian orbital elements. Furthermore, we show that the comet's coma expands quite linearly with a velocity of about 1650 km/s between October and December 2007. The photometric monitoring of comet 17/P Holmes shows that its photometric activity level decreased by about 5.9 mag within 105 days after its outburst (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

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.     

Comet P/Holmes, 1892III: A case of duplicity?

The observations of comet P/Holmes 1892III, exhibiting two 8- to 10-mag bursts, have been carefully analyzed. The phenomena are consistent with the grazing encounter of a small satellite with the nucleus on November 4.6, 1892, and the final encounter on January 16.3, 1893. The grazing encounter produced, besides the first great burst, an active area on the nucleus, which was rotating retrograde with a period of 16.3 hr and inclination of nearly 180°. After the final encounter, the spin period was essentially unchanged, but two areas became active, separated some 164° in longitude on the nucleus. After the first burst the total magnitude fell less than 2 mag from November 7 to 30 (barely naked eye) while the nuclear region remained diffuse or complex, rarely of ever showing a stellar appearance. The fading was much more rapid after the second burst (barely naked eye at maximum) while the nucleus frequently appeared stellar after the first day. It seems reasonable to conclude that the grazing encounter distributed a volume of large chunks in the neighborhood of the nucleus, maintaining activity for weeks. The final encounter activated a new area on the nucleus, the shock and fall back disturbing the area already exposed by the grazing encounter. Several details of this scenario are fitted rather well.     

A 3-D global MHD model for the effect of neutral jets during the Deep Space 1 Comet 19P/Borrelly flyby

The Deep Space 1 (DS1) spacecraft passed the sunward side of Comet 19P/Borrelly in 2001. Along its relatively north-south orbit, a set of plasma density and velocity measurement revealed a northward shift of the plasma boundaries and the mass loading peak. Both onboard and ground based telescopes found evidence for asymmetric distribution of the dust and neutrals. In this paper, five mass-loading patterns are studied to present the first study of the effect of non-spherical neutral distribution profiles on the solar wind-cometary plasma interaction environment. Using magnetohydrodynamic simulations, it is found that a combination of Gaussian and cosine neutral jet distribution, with cosine being the major part, can fit the DS1 general plasma measurement well, with a total gas production rate of around 5×¹⁰²⁸ s⁻¹. These model-data comparisons indicate that the general plasma distribution around Comet Borrelly can be explained with its aspherical neutral jet distribution. However, such neutral jets by themselves are insufficient to produce the density offset in the central peak. Kinetic effects, such as finite gyroradius may be required to create the offset plasma peak.     

SOHO/LASCO observation of an outburst of Comet 2P/Encke at its 2000 perihelion passage

Comet 2P/Encke was observed with the SOHO/LASCO C2 and C3 coronagraphs over a time interval of 11 days, starting 4 days before its September 2000 perihelion passage and through several broadband visible filters. The lightcurve reveals an outburst which started 4.9 days after perihelion, with the brightness of the coma increasing by 1.5 mag in just a few hours and progressively decreasing thereafter, probably going back to its original state in about 9 days. The color information indicates that an approximately solar color continuum was detected, implying that the observed signals were dominated by solar light scattered off submillimetric dust grains. We propose that the rapid migration of the subsolar point over the southern hemisphere during the perihelion passage activates one or several new active regions enriched in submillimetric grains, with the observed outburst corresponding to the initial blow-off of their mantle. This scenario is consistent with other observations and implies that the south polar region of the nucleus of 2P/Encke has very distinct properties.     

The 1991 outburst of Comet P/Metcalf-Brewington: a study of the light curve

The light curve of the outburst of Comet P/Metcalf-Brewington (1991 I) is investigated. The parameters of the light curve are quite uncertain as it may be fitted both in terms of variation with heliocentric distance and exponential decay with time however, after maximum a empirical fit ofm ₁ = 4.1 + 5 log Delta + 14.8 logr gives a good fit to the data. The absolute magnitude brightened at least 11 magnitudes in outburst. The outburst is most plausibly caused by the release of micron sized grains as, were the outbust due to dust emission, it would imply the expulsion of a large fraction of the total mass of the nucleus.     

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.     

CCD surface photometry of an outburst of P/Schwassmann-Wachmann 1 comet

We have observed an outburst of P/Schwassmann-Wachmann 1 (SW1) on 25–26 July 1987 using a CCD imager with R-band filter. The total brightness increased from 15.5 mg (25d) to 14.4 mag (26d) during tabout 24 hrs. The southward elongated coma of 25 was detected. The radial surface brightness (B) profiles are plotted against apparent distance p from the nucleus. The logarithmic derivative k = d ln B/d ln p for the inner coma is found to have steepened from k = –1.40 (25d) to k = –1.69 (26d), whereas that for the outer coma showed no appreciable change (k = –1.19 ~ –1.22). The ellipticity of the isophotal contour of the inner coma increased about 15% fro 25d to 6d. It is concluded that the scale of this outburst was smaller than the typical ones whose magnitude change is 5–8 mag. From recent findings on the outburst natur eof SW1 including ours, a working model of the nuclues is proposed.     

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.     

2006 Fragmentation of Comet 73P/Schwassmann-Wachmann 3B observed with Subaru/Suprime-Cam

The fragmentation of the split Comet 73P/Schwassmann-Wachmann 3 B was observed with the prime-focus camera Suprime-Cam attached to the Subaru 8.2-m telescope. The fragmentation revealed dozens of miniature comets [Fuse, T., Yamamoto, N., Kinoshita, D., Furusawa, H., Watanabe, J., 2007. Publ. Astron. Soc. Jpn. 59 (2), 381-386]. We analyzed the Subaru/Suprime-Cam images, detecting no fewer than 154 mini-comets, mostly extending to the southwest. Three were close to the projected orbit of fragment B. We applied synchrone-syndyne analysis, modified for rocket effect analysis, to the mini-fragment spatial distribution. We found that most of these mini-comets were ejected from fragment B by an outburst occurring around 1 April 2006, and three fragments on the leading side of nucleus B could have been released sunward on the previous return. Several fragments might have been released by successive outbursts around 24 April and 2 May 2006. The ratio of the rocket force to solar gravity was 7-23 times larger than that exerted on fragment B. No significant color variation was found. The mean color index, V-R = 0.50 ± 0.07, was slightly redder than that of the Sun and similar to that of the largest fragment, C, which suggests that these mini-fragments were detected mainly through sunlight reflected by dust particles and materials on the nuclei. We examined the surface brightness profiles of all detected fragments and estimated the sizes of 154 fragments. We found that the radius of these mini-fragments was in the 5- to 108-m range (equivalent size of Tunguska impactor). The power-law index of the differential size distribution was q = −3.34 ± 0.05. Based on this size distribution, we found that about 1-10% of the mass of fragment B was lost in the April 2006 outbursts. Modeling the cometary fragment dynamics [Desvoivres, E., Klinger, J., Levasseur-Regourd, A.C., Lecacheux, J., Jorda, L., Enzian, A., Colas, F., Frappa, E., Laques, P., 1999. Mon. Not. Roy. Astron. Soc. 303 (4), 826-834; Desvoivres, E., Klinger, J., Levasseur-Regourd, A.C., Jones, G.H., 2000. Icarus 144, 172-181] revealed that it is likely that mini-fragments smaller than ∼10-20 m could be depleted in water ice and become inactive, implying that decameter-sized comet fragments could survive against melting and remain as near-Earth objects. We attempted to detect the dust trail, which was clearly found in infrared wavelengths by Spitzer. No brightness enhancement brighter than 30.0 mag arcsec⁻² (3σ) was detected in the orbit of fragment B.     

Estimating the nucleus density of Comet 19P/Borrelly

The nucleus bulk density of Comet 19P/Borrelly has been estimated by modeling the sublimation-induced non-gravitational force acting upon the orbital motion, thereby reproducing the empirical perihelion advance (i.e., the shortening of the orbital period). The nucleus has been modeled as a prolate ellipsoid, covered by various surface activity maps which reproduce the observed water production rate. The theoretical water production rate of active areas has been obtained by applying a sophisticated thermophysical model. This model takes into account net sublimation of ice and thermal reradiation from the surface, solid state conductivity, sub-surface sublimation and recondensation, mass and heat transport by diffusing gas, layer absorption of solar energy, a full treatment of local time-dependent illumination conditions, and a detailed consideration of nucleus/coma interaction mechanisms. The outgassing properties of the modeled nucleus are physically consistent with the gas kinetic structure of the innermost coma since the molecular backflux and surface gas density required in the thermophysical model (as functions of the nucleus surface temperature and the sub-surface temperature profile) have been obtained from Direct Simulation Monte Carlo modeling of inelastic intermolecular collisions in the cometary Knudsen layer. The calculation of local normal forces acting on the nucleus due to outgassing has been made within the same framework—recoil and/or impact momentum transfer to the nucleus caused by sublimating molecules and by recondensing and/or scattered coma molecules is therefore evaluated in accordance with local nucleus/coma conditions. According to this model, the density is found to be 100-300 kg, depending on the applied spin axis orientation and surface activity map. This range can be narrowed down to 180-300 kg by also requiring that the empirical changes (per orbital revolution) of the argument of perihelion and the longitude of the ascending node are reproduced.     

A simulation experiment on the GM estimation for Comet 133P/Elst-Pizarro

In China’s asteroid mission to be launched around 2025,(7968)133 P/Elst-Pizarro(hereafter 133 P)will be the second target,after a visit to asteroid(469219)Kamo’oalewa.This paper describes a simulation of precise orbit determination for the spacecraft around comet 133 P,as well as estimation of its gravitational parameter(GM)value and the solar radiation pressure coefficient Cr for the spacecraft.Different cometocentric distances of 200,150 and 100 km orbits are considered,as well as two tracking modes:exclusive two-way range-rate mode(Earth station to spacecraft)and combinations of two-way range-rate and local spacecraft onboard ranging to the comet.Compared to exclusive two-way range-rate,the introduction of local ranging observables improves the final GM uncertainties by up to one order of magnitude.An ephemeris error in the orbit of 133 P is also considered,and we show that,to obtain a reliable estimate of the GM for 133 P,this error cannot exceed a one km range.     

Deep Impact, Stardust-NExT and the behavior of Comet 9P/Tempel 1 from 1997 to 2010

We present observational data for Comet 9P/Tempel 1 taken from 1997 through 2010 in an international collaboration in support of the Deep Impact and Stardust-NExT missions. The data were obtained to characterize the nucleus prior to the Deep Impact 2005 encounter, and to enable us to understand the rotation state in order to make a time of arrival adjustment in February 2010 that would allow us to image at least 25% of the nucleus seen by the Deep Impact spacecraft to better than 80 m/pixel, and to image the crater made during the encounter, if possible. In total, ∼500 whole or partial nights were allocated to this project at 14 observatories worldwide, utilizing 25 telescopes. Seventy percent of these nights yielded useful data. The data were used to determine the linear phase coefficient for the comet in the R-band to be 0.045 ± 0.001 mag deg⁻¹ from 1° to 16°. Cometary activity was observed to begin inbound near r ∼ 4.0 AU and the activity ended near r ∼ 4.6 AU as seen from the heliocentric secular light curves, water-sublimation models and from dust dynamical modeling. The light curve exhibits a significant pre- and post-perihelion brightness and activity asymmetry. There was a secular decrease in activity between the 2000 and 2005 perihelion passages of ∼20%. The post-perihelion light curve cannot be easily explained by a simple decrease in solar insolation or observing geometry. CN emission was detected in the comet at 2.43 AU pre-perihelion, and by r = 2.24 AU emission from C₂ and C₃ were evident. In December 2004 the production rate of CN increased from 1.8 × 10²³ mol s⁻¹ to Q_CN = 2.75 × 10²³ mol s⁻¹ in early January 2005 and 9.3 × 10²⁴ mol s⁻¹ on June 6, 2005 at r = 1.53 AU.     

The dust coma of Comet P/Giacobini-Zinner in the infrared

We present 1-20 micrometers photometry of P/Giacobini-Zinner obtained at the NASA Infrared Telescope Facility, during 1985 June-September (r = 1.57-1.03 AU). A broad, weak 10 micrometers silicate emission feature was detected on August 26.6; a similar weak emission feature could have been hidden in the broadband photometry on other dates. The total scattering and emitting cross section of dust in the inner coma was similar to that in other short-period comets, but a factor of 10 (r = 1.56 AU) to 100 (r = 1.03 AU) lower than the amount of dust in Comet Halley. The thermal emission continuum can be fit with models weighted toward either small or large absorbing grains. The dust production rate near perihelion was approximately 10(5) g/s (small-grain model) to approximately 10(6) g/s (large-grain model). The corresponding dust/gas mass ratio on August 26 was approximately 0.1-1. A silicate-rich heterogeneous grain model with an excess of large particles is compatible with the observed spectrum of Giacobini-Zinner on August 26. Thus, weak or absent silicate emission does not necessarily imply an absence of silicates in the dust, although the abundance of silicate particles < or = 1 micrometer radius must have been lower than in Comet Halley.     

A Multi-Wavelength Simultaneous Study of the Composition of the Halley Family Comet 8P/Tuttle

We report on simultaneous optical and infrared observations of the Halley Family comet 8P/Tuttle performed with the ESO Very Large Telescope. Such multi-wavelength and coordinated observations are a good example of what can be done to support space missions. From high resolution optical spectroscopy of the CN (0,0) 388 nm and NH₂ (0,9,0) 610 nm bands using UVES at UT2 we determined ¹²C/¹³C = 90 ± 10 and ¹⁴N/¹⁵N = 150 ± 20 in CN and we derived a nuclear spin temperature of NH₃ of 29 ± 1 K. These values are similar to those found in Oort-Cloud and Jupiter Family comets. From low resolution long slit spectroscopy with FORS1 at UT2 we determined the CN, C₃ and C₂ production rates and the parent and daughter scale lengths up to 5.2 10⁵ km tailward. From high resolution IR spectroscopy with CRIRES at UT1 we measured simultaneously the production rates and mixing ratios of H₂O, HCN, C₂H₂, CH₄, C₂H₆, and CH₃OH.     

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.     

On the origin of the unusual orbit of Comet 2P/Encke

The orbit of Comet 2P/Encke is difficult to understand because it is decoupled from Jupiter—its aphelion distance is only 4.1 AU. We present a series of orbital integrations designed to determine whether the orbit of Comet 2P/Encke can simply be the result of gravitational interactions between Jupiter-family comets and the terrestrial planets. To accomplish this, we integrated the orbits of a large number of objects from the trans-neptunian region, through the realm of the giant planets, and into the inner Solar System. We find that at any one time, our model predicts that there should be roughly 12 objects in Encke-like orbits. However, it takes roughly 200 times longer to evolve onto an orbit like this than the typical cometary physical lifetime. Thus, we suggest that (i) 2P/Encke became dormant soon after it was kicked inward by Jupiter, (ii) it spent a significant amount of time inactive while rattling around the inner Solar System, and (iii) it only became active again as the ν₆ secular resonance drove down its perihelion distance.