Long-term Evolution of the Outgassing of Comet Hale-Bopp From Radio Observations

C/1995 O1 (Hale-Bopp) has been observed on a regular basis since August 1995 at millimetre and submillimetre wavelengths using IRAM, JCMT, CSO and SEST radio telescopes. The production rates of eight molecular species (CO, HCN, CH₃OH, H₂CO,H₂S, CS, CH₃CN,HNC) have been monitored as a function of heliocentric distance(r_h from 7 AU pre-perihelion to 4 AU post-perihelion. As comet Hale-Bopp approached and receded from the Sun, these species displayed different behaviours. Far from the Sun, the most volatile species were found in general relatively more abundant in the coma. In comparison to other species, HNC, H₂CO and CS showed a much steeper increase of the production rate with decreasing r_h. Less than 1.5 AU from the Sun, the relative abundances were fairly stable and approached those found in other comets near 1 AU. The kinetic temperature of the coma, estimated from the relative intensities of the CH₃OH and CO lines, increased with decreasing r_h, from about10 K at 7 AU to 110 K around perihelion. The expansion velocity of the gaseous species, derived from the line shapes, also increased with a law close tor_h ³. This revised version was published online in July 2006 with corrections to the Cover Date.     

The 1995–2002 Long-Term Monitoring of Comet C/1995 O1 (HALE–BOPP) at Radio Wavelength

The bright comet Hale–Bopp provided the first opportunity to follow the outgassing rates of a number of molecular species over a large range of heliocentric distances. We present the results of our observing campaign at radio wavelengths which began in August 1995 and ended in January 2002. The observations were carried out with the telescopes of Nançay, IRAM, JCMT, CSO and, since September 1997, SEST. The lines of nine molecules (OH, CO, HCN, CH₃OH, H₂CO, H₂S, CS, CH₃CN and HNC) were monitored. CS, H₂S, H₂CO, CH₃CN were detected up to r_h= 3–4 AU from the Sun, while HCN and CH₃OH were detected up to 6 AU. CO, which is the main driver of cometary activity at heliocentric distances larger than 3–4 AU, was last detected in August 2001, at r_h= 14 AU. The gas production rates obtained from this programme contain important information on the nature of cometary ices, their thermal properties and sublimation mechanisms.Line shapes allow to measure gas expansion velocities, which, at large heliocentric distances, might be directly connected to the temperature of the nucleus surface. Inferred expansion velocity of the gas varied as r_h ^-0.4 within 7 AU from the Sun, but remained close to 0.4 km s^-1 further away. The CO spectra obtained at large r_hare strongly blueshifted and indicative of an important day-to-night asymmetry in outgassing and expansion velocity. The kinetic temperature of the coma, estimated from the relative intensities of the CH₃OH and CO lines, increased with decreasing r_h, from about 10 K at 7 AU to 110 K around perihelion.     

Thermal modeling of cometary nuclei

A new model of the sublimation of volatile ices from a cometary nucleus has been developed which includes the effects of diurnal heating and cooling, rotation period and pole orientation, and thermal properties of the ice and subsurface layers. The model also includes the contribution from coma opacity, scattering, and thermal emission, where the properties of the coma are derived from the integrated rate of volatile production by the nucleus. The model is applied to the specific case of the 1986 apparition of Halley's comet. It is found that the generation of a cometary dust coma actually increases the total energy reaching the Halley nucleus. This results because of the significantly greater geometrical cross section of the coma as compared with the bare nucleus, and because the coma provides an essentially isotropic source of multiply scattered sunlight and thermal emission over the entire nucleus surface. For Halley, the calculated coma opacity is approximately 0.2 at 1 AU from the Sun, and 1.2 at perihelion (0.587 AU). At 1 AU this has little effect on dayside temperatures (maximum ≈200°K) but raises nightside temperatures (minimum ≈150°K) by about 40°K. At perihelion the higher opacity results in a nearly isothermal nucleus with only small diurnal and latitudinal temperature variations. The general surface temperature is 205°K with a maximum of 209°K at local noon on the equator. Some possible consequences of the results with respect to the generation of nongravitational forces, observed volatile production rates for comets, and cometary lifetimes against sublimation are discussed.     

Observations of the inner coma of C/1995O1 (Hale-Bopp)—gas and dust production

Hale-Bopp (C/1995 O1) was the most productive recent comet observed in terms of gas and dust output. Since its discovery in 1995 at a distance of 7.14 AU from the Sun, the comet has been well observed, revealing the dynamics of a rare and large comet. Hale-Bopp showed strong emissions of the principle cometary gases CN, C₃, and C₂, as well as an abundance of dust. The production rates of these gases were found to be 1.45×10²⁸, 1.71×10²⁸, and , respectively, with dust production, in terms of Afρ, , as measured in the green continuum (5260 Å). The observations for this paper are presented in two groups spanning 10 days each, one group centered near 32 days prior to and the other 21 days after perihelion. The averages of dust and gas production rates show a slightly higher value for each prior to perihelion than after perihelion, consistent with a possible peak in production a few weeks prior to perihelion passage.     

Lessons of Comet Hale-Bopp for Coma Chemistry: Observations and Theory

Many new cometary molecules — both parents and daughters — were detected in the exceptionally productive comet C/1995 O1 (Hale-Bopp).The space distribution of several of these species could be investigated from radio interferometry or from long-slit spectroscopy in the infrared. The distinction between parent species — directly sublimated from nucleus ices — and secondary species — resulting from chemical processing in the coma or produced by a secondary source — is not always clear. It is important to assess whether or not observed minor species (HCOOCH₃, HCOOH...) could be synthesized by chemical reactions favoured by the high density of the coma of comet Hale-Bopp. Chemical modelling by Rodgers and Charnley suggests that this is notthe case. CO and H₂CO are abundant cometary species which partly come from distributed sources. The nature of these sources is still a mystery. A special case, now well documented, is that of HNC, for which the abundance evolution with heliocentric distance could be observed in comet Hale-Bopp and which was observed in several much less productive comets.     

Radio Line Observations Of Molecular And Isotopic Species In Comet C/1995 O1 (Hale-Bopp)

We present here a review of the radio observations of the remarkable comet Hale-Bopp C/1995 O1 in which most major radio astronomical facilities have been involved. These observations started in August 1995, soon after the discovery of the comet (it was then at ∼7 AU from the sun), and well before its perihelion on April 1st, 1997; they are still going on, hopefully up to end of 1998. Extended cartographies have been obtained using multibeam receivers and on-the-fly techniques. High spatial resolution (a few ″) has been achieved with interferometers. Submillimetric observations are playing an increasing role, and high resolution (R ∼ 10⁶−10⁷) spectroscopy of cometary lines is now performed from decimetric to submillimetric wavelengths. The number of species observed at radio wavelengths now reaches ∼28,when it was ∼14 for comet C/1996 B2 Hyakutake. Most of these species are parent molecules. However, ions have been observed for the first time at radio wavelengths, and their velocities measured. Several isotopic species (involving D,¹³C,³⁴S,¹⁵N) have been sought, allowing isotopic enrichment determinations. The abundances of cometary molecules present many similarities and some differences with the abundances of interstellar molecules in regions where grain mantles are believed to be evaporated to the gas phase (hot cores, bipolar flows). They will be discussed for their implications on the origin of cometary ices and of comets themselves. This revised version was published online in July 2006 with corrections to the Cover Date.     

MHD Simulation of Comets: The Plasma Environment of Comet Hale-Bopp

MHD simulation results of the interaction of the expanding atmosphere of comet Hale-Bopp with the magnetized solar wind are presented. At the upstream boundary a supersonic and superalfvénic solar wind enters into the simulation box 25 million km upstream of the nucleus. The solar wind is continuously mass loaded with cometary ions originating from the nucleus. The effects of photoionization, recombination and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our MUSCL-type upwind numerical technique, MAUS-MHD (Multiscale Adaptive Upwind Scheme for MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. Detailed simulation results for the plasma environment of comet Hale-Bopp for slow and fast solar wind conditions are presented. We also calculate synthetic H₂O⁺, CO⁺ and soft x-ray images for observing conditions on April 11, 1997. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Infrared Spectrum of Comet Hele–Bopp as Seen by the Infrared Space Observatory

Spectra of comet C/1995 O1 (Hale-Bopp) were obtained with the Infrared Space Observatory (ISO) at medium resolution with the grating spectrometer in the photometer (PHT-S) and/or at high resolution with the short wavelength spectrometer (SWS) and long wavelength spectrometer (LWS) in April 1996 (Crovisier et al., 1996), September–October 1996 (Crovisier et al., 1997a, b) and December 1997, at distances from the Sun of 4.6, 2.9 and 3.9 AU, respectively. For the first time, high-resolution spectra of a comet covering the entire 2.4 to 200 μm spectral range were obtained. The vibrational bands of H₂O, CO₂ and CO are detected in emission with PHT-S. Relative production rates of 100:22:70 are derived for H₂O:CO₂:CO at 3 AU pre-perihelion. H₂O is observed at high spectral resolution in the ν₃ group of bands around 2.7 μm and the ν₂ group around 6 μm with SWS, and in several rotational lines in the 100–180 μm region with LWS. The high signal-to-noise ratio of the ν₃ band observed on September–October 1996 allows accurate determinations of the water rotational temperature (28 K) and of its ortho-to-para ratio(2.45 ± 0.10, which significantly differs from the high temperature limit and corresponds to a spin temperature of 25 K). Longward of 6 μm the spectrum is dominated by dust thermal continuum emission, upon which broad emission features are superimposed. The wavelengths of the emission peaks correspond to those of Mg-rich crystalline olivine (forsterite). In the September–October 1996 spectra, emission features at 45 and 65 μm and possible absorption at 2.9–3.2 μm suggest that grains of water ice were present at 3 AU from the Sun. The observations made post-perihelion in late December 1997 led to the detections of H₂O, CO₂ and CO at 3.9 AU from the Sun (Figures 1 and 2). The production rates were ≈3.0 × 10²⁸,3.5 × 10²⁸ and ≈1.5 × 10²⁹ s^-1, respectively. This corresponds to H₂O:CO₂:CO = 100:110:500 and confirms that at such distances from the Sun, cometary activity is dominated by sublimation of CO and CO₂ rather than by H₂O. This revised version was published online in July 2006 with corrections to the Cover Date.     

Narrow-Band Photometry of Comet C/1995 O1 (Hale-Bopp)

Photometric observations of comet C/1995 O1 (Hale-Bopp) carried out at the Stará Lesná Observatory since February to April 1997 are analyzed and discussed. Emission band fluxes and continuum fluxes are presented, from which the total numbers of molecules in the columns of the coma encircled by diaphragms are calculated. The production rates are estimated from the conventional Haser model. We found that the photometric exponent of dust contribution two months prior perihelion was n = 5.2. The photometric exponent n of the cometary magnitude solely to the C₂ emission alone equals 3.3 and that of CN equals 2.5. These values can be explained by a fact that the maximums of production rates of the gases were reached between March 2and 12 and not at the perihelion as it is valid for dust. These results are compared with the values of 1P/Halley (1986 III) under the similar conditions, obtained with the same method and instrument. C/Hale-Bopp exhibited 4.1 times more molecules radiating the CN-emission than 1P/Halley in the same column of the coma. The continuum flux of C/Hale-Bopp was also very strong. The ratios (to 1P/Halley) are 94:1 (Cont. 484.5) and 74:1 (Cont. 365.0). The cometary colour was the same as that of the Sun. This revised version was published online in July 2006 with corrections to the Cover Date.     

A CCD Image Archive Of Comet C/1995 O1 (Hale-Bopp): Dust Expansion Velocities

An optical archive of comet Hale-Bopp was obtained over 28 nights from June 1996 to May 1997 covering a range of heliocentric distances of r = 4.12 − 0.91. The CCD images were obtained with the 0.8-m University of Washington Manastash Ridge Observatory and the 0.25-m Pennsylvania State University at Eriemobile observatory. The reduced images clearly document the evolution of features of the coma, which were dominated by radial jetsfrom r = 4.1 − 2.3 AU, curved fan-shaped jets from r = 1.2 − 1.0 AU, and arcs or shells for 0.9 < r < 1.0 AU. We present the images and derive an average dust expansion velocity of V_dust = 0.67 ± 0.07 km s⁻¹ from measurements of the arc features. The images in thearchive will be made available to other researchers for further study. This revised version was published online in July 2006 with corrections to the Cover Date.     

Infrared Observations Of Volatile Molecules In Comet Hale-Bopp

Infrared observations of comets C/1996 B2 (Hyakutake) and C/1995 O1 (Hale-Bopp) benefited from the high spectral resolution and sensitivity of échelle spectrometers now equipping ground-based telescopes and from the availability of the Infrared Space Observatory (ISO). From the ground, several hydrocarbons were unambiguously detected for the first time: CH₄, C₂H₂, C₂H₆. Water was observed through several of its hot vibrational bands, escaping telluric absorption. CO, HCN, NH₃ and OCS were also observed, as well as several radicals. This permitted the evaluation of molecular production rates, of rotational temperature, and — taking advantage of the 1-D imaging of long-slit spectroscopy — of the space distribution of these species. With ISO, carbon dioxide was directly observed for the second time in a comet (after its detection from the Vega probes in P/Halley). The spectrum of water was investigated in detail (several bands of vibration and far-infrared rotational lines), permitting the evaluation of the rotational temperature of water, and of it spin temperature from the ortho-to-para ratio. Water ice was identified in the grains of Hale-Bopp as far as 7 AU from the ground and possibly at 3 AU with ISO. The composition of cometary volatiles appears to be strikingly similar to that of interstellar ices. This revised version was published online in July 2006 with corrections to the Cover Date.     

Narrowband Electrophotometry of Comet Hale-Bopp (C/1995 O1) Near Perihelion

We present the results of narrowband photoelectric observations of comet Hale-Bopp near perihelion, which were made with the AZT-14telescope at the station Lisniki of Kyiv University Astronomical Observatory. The standard set of IHW cometary filters was used. The total number of cometary observations was more than 500 during the interval March 13–April 29, 1997. The comet's nuclear gas production rates Q of C₂ and C₃ were calculated using the stellar calibration for these filters. The calculations used the Haser model for a neutral cometary atmosphere. The value Afρ, which is characteristic of the dust production rate, was determined too. Mean values of the production rates near perihelion are log Q(C₂) = 28.4; log Q(C₃) = 27.2; (Afρ)_BC = 6.0. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Early Investigation Of The Striated Tail Of Comet Hale-Bopp (C/1995 O1)

The Sekanina-Farrell particle fragmentation model for the striated tails of dust comets is successfully applied to two images of comet Hale-Bopp to study the motions of 12 striae in a time span of March 12–15, 1997. There is evidence for recurring outbursts with a periodicity of 11^h21^m, consistent with results based on analysis of dust jets. The ejecta in all the striae appear to have been released from one source on the nucleus between the end of January and the second half of February 1997, some 60 to 40 days before perihelion. The parent particles were subjected to a radiation pressure acceleration of β_p ≃ 0.55 and their fragmentation lifetimes in 11 of the 12 striae were practically constant and equal to 13–15 days, when normalized to 1 AU from the Sun. Brief analysis of Watanabe et al.'s measurements of striae on their images from March 5–9, 1997 shows even shorter fragmentation lifetimes for the parent particles, mostly about 7–11 days at1 AU. This revised version was published online in July 2006 with corrections to the Cover Date.     

On liquid phases in cometary nuclei

In this paper we review the relevant literature and investigate conditions likely to lead to melting of H₂O ice, methanol (CH₃OH) ice, ethane (C₂H₆) ice and other volatile ices in cometary nuclei. On the basis of a heat balance model which takes account of volatiles loss, we predict the formation of occasional aqueous and hydrocarbon liquid phases in subsurface regions at heliocentric distances, r_h of 1–3 AU, and 5–12 AU, respectively. Low triple-point temperatures and low vapour pressures of C₂H₆, C₃H₈, and some higher-order alkanes and alkenes, favour liquid phase formation in cometary bodies at high r_h. Microporosity and the formation of a stabilization crust occluding the escape of volatiles facilitate liquid-phase formation. Characteristics of the near-surface which favour subsurface melting include; low effective surface emissivity (at low r_h), high amorphous carbon content, average pore sizes of ～10 μm or less, presence of solutes (e.g. CH₃OH), mixtures of C₂–C₆ hydrocarbons (for melting at high r_h), diurnal thermal cycling, and slow rotation rate. Applying the principles of soil mechanics, capillary forces are shown to initiate pre-melting phenomena and subsequent melting, which is expected to impart considerable strength of ～10⁴ Pa in partially saturated layers, reducing porosity and permeability, enhancing thermal conductivity and heat transfer. Diurnal thermal cycling is expected to have a marked effect on the composition and distribution of H₂O ice in the near-surface leading to frost heave-type phenomena even where little if any true melting occurs. Where melting does take place, capillary suction in the wetted zone has the potential to enhance heat transfer via capillary wetting in a low-gravity environment, and to modify surface topography creating relatively smooth flat-bottomed features, which have a tendency to be located within small depressions. An important aspect of the “wetted layer” model is the prediction that diurnal melt–freeze cycles alter the mixing ratio vs. depth of solutes present, or of other miscible components, largely through a process of fractional crystallization, but also potentially involving frost heave. Wetted layers are potentially durable and can involve significant mass transport of volatile materials in the near-surface, increasing in extent over many rotations of the nucleus prior to and just after perihelion passage, and causing stratification and trapping of the lowest-melting mixtures at depths of several metres. A possible mechanism for cometary outbursts is proposed involving a heat pulse reaching the liquid phase in the deepest wetted zone, leading to supersaturation and triggering the sudden release under pressure of dissolved gases, in particular CO₂, CO, CH₄ or N₂, contained beneath a consolidated near-surface layer. This study indicates that liquid water can persist for long periods of time in the near-surface of some intermediate-sized bodies (10²–10³ km radius) within protoplanetary discs.     

Discovery Of Mg-Rich Pyroxenes In Comet C/1995 O1 (Hale-Bopp): Pristine Grains Revealed At Perihelion

The NASA Ames HIFOGS spectrometer observed comet C/1995 O1 (Hale-Bopp) at epochs including 96 Oct 7–14 UT (2.8 AU), 97 Feb 14–15 UT (1.2 AU), 97 Apr 11 UT (0.93 AU), and 97 Jun 22, 25 UT (1.7 AU). The HIFOGS 7.5–13.5 μm spectrophotometry (R = 360 - 180) of the silicate feature at 2.8 AU is identical in shape to the ISO SWS spectra of comet Hale-Bopp (Crovisier et al., 1997); the strong 11.2 μm peak in the structured silicate feature is identified as olivine. Upon close passage to the sun, the HIFOGS spectra at 1.2 AU and 0.93 AU reveals strong peaks at 9.3 μm and 10.0 μm. The post-perihelion 10 μm silicate feature at 1.7 AU is weaker but has nearly the same shape as the pre-perihelion spectra at 1.2 AU, reverting to its pre-perihelion shape: there is no change in the dust chemistry by close passage to the sun. The appearance of the strong peaks at 9.3 μm and 10.0 μm at r_h ≲ 1.7 AU is attributed to the rise in the contribution of pryoxenes (clino-pyroxene and orthopyroxene crystals) to the shape of the feature, and leads to the hypothesis that the pyroxenes are significantly cooler than the olivines. The pyroxenes are radiating on the Wien side of the blackbody at 2.8 AU and transition to the Rayleigh-Jeans tail of the blackbody upon closer approach to the Sun. Composite fits to the observed 10 μm silicate features using IDPs and laboratory minerals shows that a good empirical fit to the spectra is obtained when the pryoxenes are about 150 K cooler than the olivines. The pyroxenes, because they are cooler and contribute signficantly at perihelion, are more abundant than the olivines. The perihelion temperature of the pyroxenes implies that the pyroxenes are more Mg-rich than the other minerals including the olivines, amorphous olivines, and amorphous pyroxenes. The PUMA-1 flyby measurements of comet P/Halley also indicated an overabundance of Mg-rich pryoxenes compared to olivines. Comet Hale-Bopp's pyroxenes are similar to pyroxere IDPs from the ’Spray‘ class, known for their D-richness and their unaltered morphologies: Hale-Bopp's Mg-rich pyroxenes may be pristine relic ISM grains. This revised version was published online in July 2006 with corrections to the Cover Date.     

The loss and depth of CO2 ice in comet nuclei

An analytical model has been developed to simulate the chemical differentiation of a homogeneous, initially unmantled cometary nucleus composed of water ice, putative unclathrated CO₂ ice, and silicate dust in specified proportions. Selective sublimation of any free CO₂ ice present in a new comet should produce a surface layer of water ice and dust overlying the undifferentiated core. This surface layer modifies the temperature of buried CO₂ ice and restricts the outflow of gaseous CO₂. On each orbit, water sublimation closer to perihelion temporarily reduces the thickness of the water ice and dust layer and liberates dust. Most of the dust is blown off the nucleus, but a small amount of residual dust remains on the surface (cf. H. L. F. Houpis, W. H. Ip, and D. A. Mendis, 1986, Astrophys. J., in press). Our model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the interior of the nucleus, restriction of CO₂ gas outflow by the water ice and dust layer, and the use of thermal conductivities for both amorphous and crystalline water ice as appropriate, featuresthat were not included in the Houpis et al. model. The model also accounts for the erosion of the water ice surface, which Houpis et al. appear to have accounted for and which is an important effect. Specifically, we investigate the effects of varying the permeability of the surface water ice layer, the mass fraction of CO₂, the orbit and the latitude, using the orbital parameters of Comets Halley and Tempel 2. It is found that CO₂ gas production should exceed H₂O gas production beyond ∼3 AU, and at 1 AU CO₂ gas production should be between 20 to 25% of H₂O gas production. The depth of CO₂ ice and the variation in the depth of CO₂ ice throughout an orbit are affected significantly by the perihelion of the orbit. The effects due to water ice permeability are significant but much less than expected on the basis of flow area. Latitude and CO₂ concentration produce relatively small effects. Under all conditions considered here, CO₂ ice should always be found within ∼1 m from the surface of comet nuclei if it is present as a free species to begin with. This result is probably generally valid for unmantled portions of most comets and qualitatively simulates the behavior of an abundant, highly volatile component in an H₂O/silicate matrix. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects. The method is applicable to other nonwater ices.     

Light curve behaviour in C/1995 O1 (Hale-Bopp) — II

Comet Hale-Bopp is, without doubt, one of the most important objects for cometary science which has ever been observed. The light curve is well-observed from a heliocentric distance 7.2 AU and some information is available at distances as great as 17 AU. This allows the photometric evolution of the comet to be studied as different volatiles dominate the activity. Three different phases of activity are seen as the heliocentric distance reduces during 1996: the initial phase of very fast brightening, characterized by a r ^–5 law; a standstill in the light curve when the brightening law reduced to r ^–1, which coincides with the initiation of water sublimation at r ~ 4 AU; and a further phase of more rapid brightening with an r ^–3.5 law, similar to the mean for comets classed as fairly new, which initiated at r ~ 4 AU.     

Evolution of cometary perihelion distances in oort cloud: Another statistical approach

The evolution of the perihelion distance distribution in the Oort cloud was studied over the age of the solar system, under the gravitational perturbations of random passing stars, using a statistical approach. These perturbations are accounted for through an empirical relation relating the change in cometary perihelion distance to the closest-approach comet-star distance; this relation is deduced from a previous study [H. Scholl, A. Cazenave, and A. Brahic, Astron. Astrophys.112, 157–166 (1982)]. Two kinds of initial perihelion distances are considered: (a) perihelion distances <2500 AU, associated with an origin of comets as icy planetesimals in the region of the giant planets, and (b) larger perihelion distances (up to 5 × 10⁴ AU), possibly representative of comet formation as satellite fragments in the accretion disk of the primitive solar nebula. Distant star-comet encounters, as well as rare close encounters, are considered. Several quantities are estimated: (i) number of “new” comets entering into the planetary region, (ii) number of comets escaping the Sun sphere of influence or lost by hyperbolic ejection and (iii) percentage of total comet loss over the age of the solar system. From these quantities, the current and original cloud populations are deduced, as well as the corresponding cloud mass, for the two types of formation scenarios.     

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.     

Modelling Gas And Dust Release From Comet Hale–Bopp

Numerical simulations of the evolving activity of comet Hale-Bopp are presented, assuming a porous, spherical nucleus, 20 km in radius, made of dust and gas-laden amorphous ice. The main effects included are: crystallization of amorphous ice and release of occluded gas, condensation, sublimation and flow of gases through the pores, changing pore sizes, and flow of dust grains. The model parameters, such as initial pore size and porosity, emissivity, dust grain size, are varied in order to match the observed activity. In all cases, a sharp rise in the activity of the nucleus occurs at a large heliocentric distance pre-perihelion, marked by a few orders of magnitude increase in the CO and the CO₂ fluxes and in the rate of dust emission. This is due to the onset of crystallization, advancing down to a few meters below the surface, accompanied by release of the trapped gases. A period of sustained, but variable, activity ensues. The emission of water molecules is found to surpass that of CO at a heliocentric distance of 3 AU. Thereafter the activity is largely determined by the behaviour of the dust. If a dust mantle is allowed to build up, the water production rate does not increase dramatically towards perihelion; if most of the dust is ejected, the surface activity increases rapidly, producing a very bright comet. This revised version was published online in July 2006 with corrections to the Cover Date.