2015 Annual Science Report

NASA Goddard Space Flight Center Reporting  |  JAN 2015 – DEC 2015

Serial Measurements of the Volatile Composition of Comet D/2012 S1 (ISON) between 1.2 and 0.34 AU from the Sun

Project Summary

The composition of ices and rocky material in cometary nuclei is central to understanding their origins, and to assessing their possible roles in delivering water and prebiotic organic compounds to the young Earth. For most comets, measurements of primary volatiles (ices contained in the cometary nucleus) exist for only a single date or for very few dates, questioning whether such ‘momentary’ measurements represent the bulk content of the nucleus. The early discovery of the dynamically new, sun-grazing comet C/2012 S1 (ISON) was extremely rare in that it permitted measurements of the abundances of sublimed ices over a large range of heliocentric distance (Rh). As part of a world-wide observing campaign, world-class astronomical observatories provided large amounts of observing time dedicated to studying Comet ISON. Using high-resolution infrared spectroscopy at Keck-2 and the NASA-IRTF, the GCA Team measured production rates for H2O and eight trace gases (CO, C2H6, CH4, CH3OH, NH3, H2CO, HCN, C2H2) on ten pre-perihelion dates that spanned heliocentric distances ranging from 1.21 to 0.34 AU. This project addressed the evolution in molecular production and composition as the comet approached the Sun. GCA members also investigated the spatial distribution of H2O in the near-nucleus coma to identify modes of water release and of heat injection by release from icy grains, and they conducted a sensitive search for HDO to test the potential delivery of Earth’s oceans by such comets.

4 Institutions
3 Teams
2 Publications
2 Field Sites
Field Sites

Project Progress

We measured production rates for H2O and eight trace molecules (CO, C2H6, CH4, CH3OH, NH3, H2CO, HCN, and C2H2) using high-resolution spectroscopy at Keck-2 and the NASA-IRTF on ten pre-perihelion dates encompassing heliocentric distances (Rh) = 1.21 – 0.34 AU. During this period Comet ISON exhibited two major outbursts, the first near 0.6 – 0.7 AU, and the second near 0.4 AU. Our data span this evolution in activity and reveal significant changes in the apparent composition of certain primary volatiles, prior to the comet’s complete disruption as it neared perihelion (approaching and within 0.1 AU of the Sun). Our rich dataset contains as yet unexplored information on the detailed processes involved.

For each measurement, spectra were extracted by summing signal spatially along the slit in an aperture centered on the peak emission intensity. From these spectra, telluric and cometary contributions were established by means of an iterative approach we developed and that is detailed in the paper (Figure 1).

Figure 1: Our iterative approach to quantifying terrestrial atmospheric absorption features and cometary molecular emission. (a) We first compare the observed cometary spectrum and a modeled spectrum of atmospheric extinction (i.e., telluric transmittance), omitting channels containing (known) strong cometary emission lines and convolving the fully resolved transmittance spectrum to the spectral resolving power of the instrument (RP = λ/δλ ~ 2.5 x 104). We fit this to establish quantitative column burdens for absorbing species in the terrestrial atmosphere. The trace in panel 'a' labeled “cometary emission” represents the observed net emission in excess of the cometary continuum, which we then compare with similarly convolved molecular fluorescence models (panel 'a', bottom traces). For each spectral setting, the goal is to iteratively fit both telluric and cometary features accurately, such that the resulting “residual” spectrum (bottom trace in panel 'b') falls within the ±1σ stochastic noise envelope.

Measured water production rates spanned two orders of magnitude, and appeared to track periods of both “outburst” and “baseline” activity as reported by other observers. A fit to our baseline measurements was consistent with a heliocentric power law Q(H2O) ~ Rh-n, with n = -3.1 ± 0.1 (Figure 2).

Figure 2: Evolution of water production in Comet ISON prior to its complete disruption. The various NIRSPEC results shown for a given date were obtained simultaneously; the points on November 7 are separated horizontally only to visualize them more clearly. Of the 32 measurements, seven (in dashed boxes) were associated with periods of enhanced water production (i.e., with “outbursts”). A heliocentric fit excluding these seven points (dot-dashed line) results in a power law [Q(H2O) = (1.89 ± 0.11) x 1028 Rh-n, with n = 3.10 ± 0.1) molecules s-1] that represents the evolution of water production in Comet ISON within the terrestrial planets region.

Abundance ratios for CO, C2H6, and CH4 with respect to H2O were relatively constant with Rh, and were smaller than their corresponding mean ratios measured among a dominant sample of Oort cloud comets. CH3OH was also depleted for Rh > 0.5 AU, but was closer to its mean value for Rh ≤ 0.5 AU. The remaining four gases exhibited higher abundance ratios within 0.5 AU: for Rh > 0.8 AU, NH3 and C2H2 were consistent with their mean values while H2CO and HCN were depleted. For Rh < 0.5 AU all four were enriched, with NH3, H2CO, and HCN increasing most (Fig. 3).

Figure 3: Evolution of molecular abundance ratios relative to simultaneously or contemporaneously measured H2O in Comet ISON. Points shown in red correspond to measurements during periods of enhanced gas production (i.e., during “outburst”). No discernable dependence of composition on “baseline” versus “outburst” activity was seen.

Spatial profiles of gas emission in ISON consistently peaked sunward of the dust continuum, which peaked anti-sunward and displayed a single peak throughout our observing period. Within 0.5 AU, NH3 showed a broad spatial distribution, possibly indicating its release in the coma (if optical depth effects were unimportant). The column abundance ratio NH2/H2O at 0.83 AU was close to the “typical” value for NH/OH from optical wavelengths, but was higher within 0.5 AU. Establishing its production rate and testing its parent (e.g., NH3) requires modeling of coma outflow, but this is deferred to a future reporting period.

Figure 4: Spatial profiles of simultaneously measured emissions from individual gases and dust (continuum) in several CSHELL settings. The distributions of species in the coma are used to obtain total production rates and also to assess the nature of volatile release. The schematic compass shows the east-west slit orientation and the sunward direction projected onto the sky plane. The solar phase angle of ISON remained near 90° on the dates shown, thus the true sunward direction was nearly in the sky plane.