Crystalline silicates in comets: How did they form?

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

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

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

Spectrocam-10 Thermal Infrared Observations Of The Dust In Comet C/1995 O1 Hale-Bopp

We describe 5 to 18 μm broadband images and R ∼ 100 spectra of comet Hale-Bopp taken with SpectroCam-10 on the 5 m Hale telescope during six runs between 1996 June and 1997 April. Our data show the gradual warming of dust grains as the comet approached perihelion. In June, the 10 μm silicate emission feature was already stronger than observed in any other comet, and it increased to about 3 times the continuum level near perihelion. Spectral structure attributed to a crystalline olivine grain component remained relatively constant, but other features associated with pyroxenes appeared to vary with heliocentric distance. This revised version was published online in July 2006 with corrections to the Cover Date.     

The composition of dust in Jupiter-family comets inferred from infrared spectroscopy

We review the composition of Jupiter-family comet (JFC) dust as inferred from infrared spectroscopy. We find that JFCs have silicate emission features with fluxes roughly 20-25% over the dust continuum (emission strength 1.20-1.25), similar to the weakest silicate features in Oort Cloud (OC) comets. We discuss the grain properties that alter the silicate emission feature (composition, size, and structure/shape), and emphasize that thermal emission from the comet nucleus can have significant influence on the derived silicate emission strength. Recent evidence suggests that grain porosity is the is different between JFCs and OC comets, but more observations and models of silicates in JFCs are needed to determine if a consistent set of grain parameters can explain their weak silicate emission features. Models of 8 m telescope and Spitzer Space Telescope observations have shown that JFCs have crystalline silicates with abundances similar to or less than those found in OC comets, although the crystalline silicate mineralogy of comets 9P/Tempel and C/1995 O1 (Hale-Bopp) differ from each other in Mg and Fe content. The heterogeneity of comet nuclei can also be assessed with mid-infrared spectroscopy, and we review the evidence for heterogeneous dust properties in the nucleus of comet 9P/Tempel. Models of dust formation, mixing in the solar nebula, and comet formation must be able to explain the observed range of Mg and Fe content and the heterogeneity of comet 9P/Tempel, although more work is needed in order to understand to what extent do comets 9P/Tempel and Hale-Bopp represent comets as a whole.     

Aperture Polarimetry And Photometry Of Comet Hale-Bopp

We present results of polarimetric and photometric observations of bright comet C/1995 O1 (Hale-Bopp) obtained at the 0.7 m telescope of Kharkov University Observatory from June 18, 1996 to April 24, 1997. The IHW and HB comet filters were used. The C₂ and C₃ production rates for Hale-Bopp are more than one order of magnitude larger and the dust production rates are more than two orders of magnitude larger than the Halley ones at comparable distances. Hence, Hale-Bopp was one of the most dusty comets. The average UC-BC and BC-RC colours of the dust were −0.02 and 0.13 mag, respectively. The polarization of comet Hale-Bopp at small phase angles of 4.8–13.0° was in good agreement with the date for comet P1/Halley at the same phase angles in spite of the fact that the heliocentric distances of comments differed nearly twice. However, at intermediate phase angles of 34–49° the polarization of comet Hale-Bopp was significantly larger than the polarization of the other dusty comets. It is the first case of such a large difference found in the continuum polarization of comets. The wavelength dependence of polarization for Hale-Bopp was steeper than for other dusty comets. The observed degree of polarization for the anti-sunward side of the coma was permanently higher than that for the sunward shell side. The polarization phase dependence of Hale-Bopp is discussed and compared with the polarization curves for other dusty comets. The peculiar polarimetric properties of comet Hale-Bopp are most likely caused by an over-abundance of small or/and absorbing dust particles in the coma. This revised version was published online in July 2006 with corrections to the Cover Date.     

Thermal Emission From The Dust Coma Of Comet Hale-Bopp And The Composition Of The Silicate Grains

The dust coma of comet Hale-Bopp was observed in the thermal infrared over a wide range in solar heating (R = 4.9–0.9 AU) and over the full wavelength range from 3 μm to 160 μm. Unusual early activity produced an extensive coma containing small warm refractory grains; already at 4.9 AU, the 10 μm silicate emission feature was strong and the color temperature was 30% above the equilibrium blackbody temperature. Near perihelion the high color temperature, strong silicate feature, and high albedo indicated a smaller mean grain size than in other comets. The 8–13 μm spectra revealed a silicate emission feature similar in shape to that seen in P/Halley and several new and long period comets. Detailed spectral structure in the feature was consistent over time and with different instruments; the main peaks occur at 9.3, 10.0 and 11.2 μm. These peaks can be identified with olivine and pyroxene minerals, linking the comet dust to the anhydrous chondritic aggregate interplanetary dust particles. Spectra at 16–40 μm taken with the ISO SWS displayed pronounced emission peaks due to Mg-rich crystalline olivine, consistent with the 11.2 μm peak. This revised version was published online in July 2006 with corrections to the Cover Date.     

ISO Observations of Isolated Herbig AE/BE Stars

We report on SWS and LWS observations of the circumstellar disks of young stars of a few solar masses. The ISO spectra of these objects present a diversity of emission features of carbon-rich and oxygen-rich grains. The similarity of the forsterite spectra observed for Comet Hale-Bopp and the Haebe star HD100546 is particularly striking and provides a new argument that huge comet swarms are formed in the disks surrounding young stars. While the data suggest that the formation of crystalline silicates in the dust disks essentially occurs when a Haebe star has already reached the main sequence, no clear correlation with stellar age only is apparent. This revised version was published online in August 2006 with corrections to the Cover Date.     

Comet Hale-Bopp,Circumstellar Dust,and the Interstellar Medium

The infrared spectrum of comet Hale-Bopp, as observed with the Infrared Space Observatory (ISO), presents astonishing similarities with those of the dusty disks surrounding young stars. The prominent spectral features of crystalline silicates that occur in comets, have been observed in the circumstellar media of young and old stars, yet not in the interstellar medium. We suggest from the analogy with the young-star spectra that the crystalline nature of the silicates in comet Hale-Bopp dates from the early evolution of the solar system, and also that an Oort Cloud or Kuiper Belt is presently forming around young stars observed by ISO. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Infrared Spectral Energy Distribution And Polarization Of Comet C/1995 O1 (Hale–Bopp) During 1997

Comets, such as C/1995 O1 (Hale-Bopp), are important to studies of the origins of the solar system because they are believed to be frozen reservoirs of the most primitive pre-solar dust grains and ices. Here, we report 1.2–18.5 μm infrared (IR) spectrophotometric and polarimetric observations of comet Hale-Bopp. Our measurements of the spectral energy distribution (SED) and IR polarization near perhelion passage suggest that emission from the coma was dominated by scattering and thermal emission from sub-micron sized dust grains. Hale-Bopp's surprising brightness may have been largely a result of the properties of its coma grains rather than the size of its nucleus. The thermal emission continuum from the grains had a superheat of S = T_color/T_BB ≥ 1.7, the peak of the 10 μm silicate emission feature was 1.7 mags above the carbon grain continuum, and the albedo (reflectivity) of the grains was ≥ 0.4 at a scattering angles, θ ≥ 135° This revised version was published online in July 2006 with corrections to the Cover Date.     

X-Ray Emission From Comet Hale–Bopp

The discovery of X-ray emission from comets has created a number of questions about the physical mechanism producing the radiation. There are now a variety of explanations for the emission, from thermal bremsstrahlung of electrons off neutrals or dust, to charge exchange induced emission from solar wind ions, to scattering of solar X-rays from attogram dust, to reconnection of solar magnetic field lines. In an effort to understand this new phenomenon, we observed but failed to detect in the X-ray the very dusty and active comet C/Hale-Bopp 1995 O1 over a two year period, September 1996 to December 1997, using the ROSAT HRI imaging photometer at 0.1–2.0 keV and the ASCA SIS imaging spectrometer at 0.5–10.0 keV. The results of our Hale-Bopp non-detections, when combined with spectroscopic imaging 0.08–1.0 keV observations of the comet by EUVE and BeppoSAX, show that the emission has the same spectral shape and strong variability seen in other comets. Comparison of the ROSAT photometry of the comet to our ROSAT database of 8 comets strongly suggests that the overall X-ray faintness of the comet was due to an emission mechanism coupled to gas, and not dust, in the comet’s coma. This revised version was published online in July 2006 with corrections to the Cover Date.     

Chemical and physical properties of gas jets in comets: I. Monte Carlo model of an inner cometary coma

We describe a 3-dimensional, time-dependent Monte Carlo model developed to analyze the chemical and physical nature of a cometary gas coma. Our model includes the necessary physics and chemistry to recreate the conditions applicable to Comet Hale-Bopp when the comet was near 1 AU from the Sun. Two base models were designed and are described here. The first is an isotropic model that emits particles (parents of the observed gases) from the entire nucleus; the second is a jet model that ejects parent particles solely from discrete active areas on the surface of the comet nucleus, resulting in coma jets. The two models are combined to produce the final model, which is compared with observations. The physical processes incorporated in both base models include: (1) isotropic ejection of daughter molecules (the observed gases) in the parent's frame of reference, (2) solar radiation pressure, (3) solar insolation effects, (4) collisions of daughter products with other molecules in the coma, and (5) acceleration of the gas in the coma. The observed daughter molecules are produced when a parent decays, which is represented by either an exponential decay distribution (photodissociation of the parent gas) or a triangular distribution (production from a grain extended source). Application of this model to the analysis the OH, C₂ and CN gas jets observed in the coma of Comet Hale-Bopp is the focus of the accompanying paper [Lederer, S.M., Campins, H., Osip, D.J., 2008. Icarus, in press (this issue)].     

Dust tail of the active distant Comet C/2003 WT42 (LINEAR) studied with photometric and spectroscopic observations

We present the study of dust environment of dynamically new Comet C/2003 WT42 (LINEAR) based on spectroscopic and photometric observations. The comet was observed before and after the perihelion passage at heliocentric distances from 5.2 to 9.5 AU. Although the comet moved beyond the zone where water ice sublimation could be significant, its bright coma and extended dust tail evidenced the high level of physical activity. Afρ values exceeded 3000 cm likely reaching its maximum before the perihelion passage. At the same time, the spectrum of the comet did not reveal molecular emission features above the reflected continuum. Reddening of the continuum derived from the cometary spectrum is nonlinear along the dispersion with the steeper slop in the blue region. The pair of the blue and red continuum images was analyzed to estimate a color of the comet. The mean normalized reflectivity gradient derived from the innermost part of the cometary coma equals to 8% per 1000 Å that is typical for Oort cloud objects. However, the color map shows that the reddening of the cometary dust varies over the coma increasing to 15% per 1000 Å along the tail axis. The photometric images were fitted with a Monte Carlo model to construct the theoretical brightness distribution of the cometary coma and tail and to investigate the development of the cometary activity along the orbit. As the dust particles of distant comets are expected to be icy, we propose here the model, which describes the tail formation taking into account sublimation of grains along their orbits. The chemical composition and structure of these particles are assumed to correspond with Greenberg’s interstellar dust model of comet dust. All images were fitted with the close values of the model parameters. According to the results of the modeling, the physical activity of the comet is mainly determined by two active areas with outflows into the wide cones. The obliquity of the rotation axis of the nucleus equals to 20° relative to the comet’s orbital plane. The grains occupying the coma and tail are rather large amounting to 1 mm in size, with the exponential size distribution of a^−4.5. The outflow velocities of the dust particles vary from a few centimeters to tens of meters per second depending on their sizes. Our observations and the model findings evidence that the activity of the nucleus decreased sharply to a low-level phase at the end of April–beginning of May 2007. About 190 days later, in the first half of November 2007 the nucleus stopped any activity, however, the remnant tail did not disappear for more than 1.5 years at least.     

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.     

Crystalline comet dust: Laboratory experiments on a simple silicate system

Abstract— Spectra for certain comets show the presence of crystalline silicate dust grains believed to have been incorporated during comet formation. While grain crystallization is widely assumed to result from the thermal annealing of precursor amorphous grains, the physical processes behind the silicate amorphous‐to‐crystalline transition are poorly understood. This makes it difficult to place constraints on the evolutionary histories of both grains and comets, and consequently, on the nebular conditions in which they formed. It has, therefore, become necessary to study this process in the laboratory using simulated grain materials. In this paper, we discuss recent results from laboratory investigations into a basic amorphous MgSiO₃ silicate annealed in the region of 1000 K. Our object is not to model the behavior of dust grains per se, but to study the underlying process of crystallization and separate the physics of the material from the astrophysics of dust grains. In our experiments, we bring together spectroscopic measurements made in the infrared with the high resolution structural probing capabilities of synchrotron X‐ray powder diffraction. The combined use of these complementary techniques provides insights into the crystallization process that would not be easily obtained if each was used in isolation. In particular, we focus on the extent to which the identification of certain spectral features attributed to crystalline phases extends to the physical structure of the grain material itself. Specifically, we have identified several key features in the way amorphous MgSiO₃ behaves when annealed. Rather than crystallize directly to enstatite (MgSiO₃) structures, in crystallographic terms, amorphous MgSiO₃ can enter a mixed phase of crystalline forsterite (Mg₂SiO₄) and SiO₂‐rich amorphous silicate where structural evolution appears to stall. Spectroscopically, the evolution of the 10 μm band does not appear to correlate directly with structural evolution, and therefore, may be a poor indicator of the degree of crystallinity. Indeed, certain features in this band may not be indicators of crystal type. However, the 20 μm band is found to be a good indicator of crystal structure. We suggest that forsterite forms from the ordering of pre‐existing regions rich in SiO₄ and that this phase separation is aided by a dehydrogenation processes that results in the evolutionary stall. The implications of this work regarding future observations of comets are discussed.     

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.     

Cometary dust properties retrieved from polarization observations: Application to C/1995 O1 Hale-Bopp and 1P/Halley

The analysis of the polarized light scattered by cometary dust particles provides information on the physical properties of the solid component of cometary comae for C/1995 O1 Hale-Bopp and 1P/Halley. A model of light scattering by a size distribution of aggregates of up to 256 submicron-sized grains (spherical or spheroidal) mixed with single spheroidal particles has been developed, with its parameters adjusted to fit the phase angle and wavelength dependence of the polarization observations. The particles are built of two materials: a non-absorbing silicates-type material and a more absorbing organic-type material. The model reproduces accurately the inversion angle and the positive branch of the polarization phase curves from the visible to the near-infrared spectral domains. A negative branch of the polarization phase curves appears in our model, although the negative branch is not deep enough to reproduce accurately the observations. Significant differences are shown between the two comets, with dominance of small grains in the coma of Comet C/1995 O1 Hale-Bopp, well fitted by a distribution of the volume-equivalent diameter, a, following a^−3.0 with a lower cutoff around 0.20 μm and an upper cutoff of at least 40 μm. For 1P/Halley, the size distribution follows a^−2.8 with a lower cutoff around 0.26 μm and an upper cutoff of about 38 μm. The relative amount of organic-type particles is larger for 1P/Halley while the amount of aggregates, significant for both comets, is larger for C/1995 O1 Hale-Bopp.     

Mid-Infrared Spectro-Imaging Observations Of Comet Hale-Bopp

Comet C/1995 O1 (Hale-Bopp) was observed with camiras, the Saclay mid-infrared camera, mounted on the Nordic Optical Telescope (Roque de los Muchachos, La Palma), from April 14th to April 25th, 1997. This observing run was part of the European campaign devoted to Hale-Bopp observations from the Canary Islands (PI R. West). camiras spectro-imaging capabilities, achieved with a Circular Variable Filter (of spectral resolution R ≈ 50), were used to obtain spectra of the comet in the N atmospheric window (8–13 μm) over a large field of view (52″ × 78″), at a spatial resolution of ∼1 arcsec. Data were also collected at shorter wavelengths (in the L and M bands). The silicate dust feature around 10 μm is clearly apparent in the data. The shape of the feature varies continuously according to the position in the field. A dust model developed at Saclay in the framework of β-Pictoris dust disk studies was used to interpret the observations. The presence of both olivine and pyroxene with a high degree of crystallinity is needed to account for the observations. A change in the size distribution of the grains allows to reproduce the spatial change of the silicate feature according to the position in the field. This revised version was published online in July 2006 with corrections to the Cover Date.     

The dust trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of Comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2 m telescope at 4.7 AU before aphelion, and at with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter β<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust environment and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of −4.1, and emission speeds for a β=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of .     

No compelling evidence of distributed production of CO in Comet C/1995 O1 (Hale-Bopp) from millimeter interferometric data and a re-analysis of near-IR lines

Based on long-slit infrared spectroscopic observations, it has been suggested that half of the carbon monoxide present in the atmosphere of Comet C/1995 O1 (Hale-Bopp) close to perihelion was released by a distributed source in the coma, whose nature (dust or gas) remains unidentified. We re-assess the origin of CO in Hale-Bopp’s coma from millimeter interferometric data and a re-analysis of the IR lines.Simultaneous observations of the CO J(1–0) (115 GHz) and J(2–1) (230 GHz) lines were undertaken with the IRAM Plateau de Bure interferometer in single-dish and interferometric modes. The diversity of angular resolutions (from 1700 to 42,000 km diameter at the comet) is suitable to study the radial distribution of CO and detect the extended source observed in the infrared. We used excitation and radiative transfer models to simulate the single-dish and interferometric data. Various CO density distributions were considered, including 3D time-dependent hydrodynamical simulations which reproduce temporal variations caused by the presence of a CO rotating jet. The CO J(1–0) and J(2–1) observations can be consistently explained by a nuclear production of CO. Composite 50:50 nuclear/extended productions with characteristic scale lengths of CO parent L_p > 1500 km are rejected.Based on similar radiation transfer calculations, we show that the CO v = 1–0 ro-vibrational lines observed in Comet Hale-Bopp at heliocentric distances less than 1.5 AU are severely optically thick. The broad extent of the CO brightness distribution in the infrared is mainly due to optical depth effects entering in the emitted radiation. Additional factors can be found in the complex structure of the CO coma, and non-ideal slit positioning caused by the anisotropy of dust IR emission.We conclude that both CO millimeter and infrared lines do not provide compelling evidence for a distributed source of CO in Hale-Bopp’s atmosphere.