Evolution of H{}_{2}O Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission

Evolution of HO Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission

Adam J. McKay Anita L. Cochran Michael A. DiSanti Neil Dello Russo Harold Weaver Ronald J. Vervack Jr. Walter M. Harris Hideyo Kawakita NASA GSFC/USRA, 8800 Greenbelt Rd, Greenbelt, MD 20771 (U.S.A.); adam.mckay@nasa.gov University of Texas Austin/McDonald Observatory, 2512 Speedway, Stop C1402 , Austin, TX 78712, (U.S.A); anita@astro.as.utexas.edu NASA Goddard Center for Astrobiology, NASA GSFC, Mail Stop 690, Greenbelt, MD 20771 (U.S.A.); Michael.A.Disanti@nasa.gov Solar System Exploration Division, Mail Stop 690, Greenbelt, MD 20771 (U.S.A) Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Rd., Laurel, MD, 20723 (U.S.A.); neil.dello.russo@jhuapl.edu, Hal.Weaver@jhuapl.edu, Ron.Vervack@jhuapl.edu Lunar and Planetary Laboratory, University of Arizona, 1629 E University Blvd., Tucson, AZ, 85721 (U.S.A);wharris@lpl.arizona.edu Koyama Astronomical Observatory, Kyoto Sangyo University, Motoyama, Kamigamo, Kita-ku, Kyoto 603-8555, Japan; kawakthd@cc.kyoto-su.ac.jp

We present HO production rates for comet C/2012 S1 (ISON) derived from observations of [OI] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8-0.44 AU. Our production rates are in agreement with previous measurements using a variety of instruments and techniques and with data from the various observatories greatly differing in their projected fields of view. The consistent results across all data suggest the absence of an extended source of HO production, for example sublimation of icy grains in the coma, or a source with spatial extent confined to the dimensions of the smallest projected field of view (in this case 1,000 km). We find that ISON had an active area of around 10 km for heliocentric distances R 1.2 AU, which then decreased to about half this value from R=1.2-0.9 AU. This was followed by a rapid increase in active area at about R=0.6 AU, corresponding to the first of three major outbursts ISON experienced inside of 1 AU. The combination of a detected outburst in the light curve and rapid increase in active area likely indicates a major nucleus fragmentation event. The 5-10 km active area observed outside of R=0.6 AU is consistent with a 50-100% active fraction for the nucleus, larger than typically observed for cometary nuclei. Although the absolute value of the active area is somewhat dependent on the thermal model employed, the changes in active area observed are consistent among models. The conclusion of a 50-100+% active fraction is robust for realistic thermal models of the nucleus. However the possibility of a contribution of a spatially unresolved distribution of icy grains cannot be discounted. As our [OI]-derived HO production rates are consistent with values derived using other methods, we conclude that the contribution of O photodissociation to the observed [OI] emission is at most 5-10% that of the contribution of HO for ISON. This is consistent with the expected contribution of O photodissociation if O/HO 4%, meaning [OI] emission can still be utilized as a reliable proxy for HO production in comets as long as O/HO 4%, similar to the abundance measured by the ROSINA instrument on Rosetta at comet 67P/Churyumov-Gerasimenko.

Comets; Comets, Coma; Comets, Composition
journal: Icarusthanks: This paper includes data taken at The McDonald Observatory of The University of Texas at Austin, as well as data collected at the W.M. Keck Observatory, Maunakea, HI, USA, operated as a scientific partnership among Caltech, UCLA, and NASA, and made possible by the generous financial support of the W.M. Keck Foundation.

, , , , , , ,

Copyright © 2018 Adam J. McKay, Anita L. Cochran, Michael A. DiSanti, Neil Dello Russo, Harold Weaver, Ronald J. Vervack, Walter M. Harris, Hideyo Kawakita

Number of pages: Evolution of HO Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission

Number of tables: 5

Number of figures: 5

Proposed Running Head:
HO Production in Comet C/2012 S1 (ISON)

Please send Editorial Correspondence to:

Adam J. McKay
8800 Greenbelt Rd
Greenbelt, MD 20771, USA.

Email: adam.mckay@nasa.gov
Phone: 301-614-5701

1 Introduction

Comet C/2012 S1 (ISON) was a dynamically new comet from the Oort Cloud, meaning it was making its first passage through the inner Solar System (Novski et al., 2012). It passed within 3 solar radii of the Sun in late November 2013, making ISON a sungrazing comet. Many sungrazing comets have been discovered, most belonging to a dynamical family called the Kreutz group. However, ISON was unique in two respects: 1) it was dynamically new, whereas the Kreutz group comets are not, and 2) more importantly, ISON was discovered over a year before its perihelion passage, meaning observations could be planned well in advance to follow its plunge toward the Sun. This meant that ISON could be observed over a very large range of heliocentric distances, which was not possible for previous sun-grazing comets, potentially providing new insights into cometary composition and activity. Because of this unique opportunity, a world-wide observing campaign was organized, providing an unprecedented data set.

Comet ISON experienced very irregular activity levels as it approached the Sun, experiencing several outbursts that may or may not have been correlated with fragmentation or disruption events (e.g. Sekanina and Kracht, 2014). Measuring the HO production rate throughout its inbound trajectory provides a window into the sublimation processes driving the outbursts.

While direct observations of HO are possible at IR wavelengths (e.g. Dello Russo et al., 2011), HO production can also be measured indirectly through its dissociation products H, O, and OH. The main emission mechanism for the [OI]6300 Å line is prompt emission after photodissociation of an HO molecule (Festou and Feldman, 1981). Therefore observations of this line have been used as a reliable proxy for HO production in the past (e.g. Morgenthaler et al., 2001, 2007; Fink, 2009; McKay et al., 2015). Also, OH in the coma results from photodissociation of HO, and OH is often used as a proxy for HO at NUV (e.g. A’Hearn et al., 1995; Opitom et al., 2015), IR (e.g. Bonev et al., 2006; Dello Russo et al., 2011), and radio wavelengths (e.g. Biver et al., 2002). Lyman- emission has also been used to derive HO production rates, particularly with the SOHO satellite (e.g. Combi et al., 2013).

We present analysis of observations of [OI] and OH emission in comet ISON as a proxy for HO production. We obtained multiple observations spanning the time period from September-November 2013, covering a heliocentric distance range of 1.8-0.44 AU. Our observations of [OI] in particular, which have excellent sensitivity due to the abilities of optical detectors and the prompt nature of the [OI] emission mechanism, cover larger heliocentric distances than other studies, providing a characterization of ISON’s activity earlier in the apparition. This paper is organized as follows: in section 2 we describe our observations and reduction and analysis procedures. Section 3 presents our results. In section 4 we discuss our results in the context of the wider observing campaign, including comparison to other measurements of the HO production rate obtained with different methods. Section 5 presents a summary of our conclusions.

2 Observations and Data Analysis

2.1 Observations

We obtained spectra of C/2012 S1 (ISON) using three facilities. The ARCES instrument is mounted on the Astrophysical Research Consortium 3.5-m telescope at Apache Point Observatory (APO) in Sunspot, New Mexico. ARCES provides a spectral resolution of R = 31,500 and a spectral range of 3500-10,000 Å with no interorder gaps. More specifics for this instrument are discussed elsewhere (Wang et al., 2003). In mid-October we obtained spectra using the Tull Coude Spectrograph mounted on the 2.7-m Harlan J. Smith telescope at McDonald Observatory. The spectral range is the same as ARCES, but the Tull Coude provides twice the spectral resolution, with R 60,000. However, unlike ARCES, the Tull Coude spectrograph has interorder gaps redward of 5700 Å. To account for this, we used a spectrograph setup such that the interorder gaps did not fall on emission features of interest, which for the work presented here constitutes the [OI]6300 Å line. We obtained observations with the HIRES instrument on Keck I in late October. We utilized the HIRESb mode, which covers wavelengths from 3000-5700 Å at R 48,000 with small gaps in coverage at the edges of the CCD detectors. Although this mode does not cover the redder wavelengths sampled by the ARCES and Tull Coude data (and therefore not the [OI]6300 Å line), the additional extension blueward into the UV provides observations of the OH A-X (0-0) band at 3080 Å that is not sampled in the ARCES or Tull Coude spectra.

The observation dates and geometries, as well as standard stars used and observing conditions are described in Table 1. For all observations we centered the slit on the optocenter of the comet. All the instruments employed have relatively narrow and short slits, specifically 1.63.2(ARCES), 0.867(HIRES), and 1.28.2(Tull Coude). We used an ephemeris generated from JPL Horizons for non-sidereal tracking of the optocenter. For short time-scale tracking, the guiding software uses a boresight technique, which utilizes optocenter flux that falls outside the slit to keep the slit on the optocenter. For the ARCES and HIRES observations, we observed a G2V star in order to remove the underlying solar continuum and Fraunhofer absorption lines. For the McDonald observations, we obtained observations through a solar port which feeds sunlight directly into the spectrograph. We obtained spectra of a fast rotating (vsin(i) 150 km s), O, B, or A star to account for telluric features and spectra of a flux standard to establish absolute intensities of cometary emission lines. The calibration stars used for each observation date are given in Table 1. We obtained spectra of a quartz lamp for flat fielding and acquired spectra of a ThAr lamp for wavelength calibration.

2.2 Data Reduction

Spectra were extracted and calibrated using IRAF scripts that perform bias subtraction, cosmic ray removal, flat fielding, and wavelength calibration. We removed telluric absorption features, removed the reflected solar continuum from the dust coma, and flux calibrated the spectra employing our standard star observations. We assumed an exponential extinction law and extinction coefficients for the observatory site when flux calibrating the cometary spectra (e.g. Hogg et al., 2001). More details of this procedure can be found in McKay et al. (2012) and Cochran and Cochran (2002). We determined slit losses for the flux standard star observations by performing aperture photometry on the slit viewer images as described in McKay et al. (2014). Slit losses introduce a systematic error in the flux calibration of 10%. Slit viewer images are not available for the McDonald and Keck data sets, so for these data sets we adopted a slit loss value based on the measured seeing and assumption of a Gaussian PSF. Accounting for possible variability in seeing and the ideal nature of this assumption, we estimate the uncertainties in slit losses are 20% in these cases. For UT September 23 and UT November 20, flux standards taken from UT October 3 and UT November 15, respectively, were used to flux calibrate the spectra since a flux standard star was not obtained those nights. This may introduce additional uncertainty that is not accounted for in our quoted measurements.

With the high airmass ( 2-3) of some of our observations and the slits employed (widths of 0.9-1.6, lengths of 3.2-8.2), the effect of differential refraction on fluxes obtained at different wavelengths could be significant. This could be particularly important for OH, whose primary emission band is in the near UV (3100 Å) and therefore suffers from more significant differential refraction than redder spectral features. To evaluate the effect of differential refraction, we employ the computed differential refraction as a function of wavelength and airmass from Filippenko (1982), and interpolate to the relevant wavelengths and airmasses for our observations. We calculate the amount of differential refraction and use a Haser model profile convolved with a Gaussian point spread function to estimate the slit loss due to the differential refraction. For all our observations we estimate the effect of differential refraction is less than 10% and is insignificant compared to the uncertainties in flux calibration noted above (10-20%). The uncertainties quoted in this paper are dominated by the uncertainty in flux calibration and not the signal-to-noise ratio of the data nor differential refraction effects.

2.3 Analysis of [OI] Emission

The [OI]6300 Å line is among the emission features present in the ARCES and Tull Coude bandpasses and can be employed to derive the HO production rate. The [OI]6300 Å line is also present as a telluric emission feature, and so a combination of high spectral resolution and adequate geocentric velocity (and therefore Doppler shift) is needed to resolve the cometary line from the telluric feature. For all observation dates except November 20 the cometary and telluric lines are fully separated. On November 20 the comet was so bright and active that the cometary line is expected to dwarf the telluric feature, therefore any influence of the telluric feature on the measured flux is negligible. An example spectrum from November 15 is shown in the top panel of Fig. 1, with a more detailed view of the [OI]6300 Å and [OI]6364 Å lines shown in Fig. 2. To measure the line flux we fit a Gaussian function to the line profile. More details for our line fitting procedure can be found in McKay et al. (2012).

For all of our observations, to check our calibration we measured the flux ratio of the [OI]6300 Å line to the [OI]6364 Å line (hereafter the red line ratio). As these lines both originate from the D state, the red line ratio is determined directly by the branching ratio into the different transitions and is not dependent on coma chemistry or physics. This value is well established at a value of 3.0 (Cochran and Cochran, 2001; Sharpee and Slanger, 2006). However, our observations in September and October yield significantly lower values for the red line ratio, in the range 1.5-2.5. Examination of the spectra revealed that at the very large Doppler shift that ISON had at the time of these observations (-50 km s), the [OI]6300 Å line falls precisely on an O telluric absorption feature, as shown in Fig. 3. Although we remove telluric signatures during our reduction process, the removal may not be perfect. On nights for which we had multiple observations, [OI]6300 Å line intensities have a large amount of scatter, evidence that our [OI]6300 Å line intensities were affected by imperfect removal of this telluric absorption feature from spectrum to spectrum. There are no telluric or solar absorption features near the [OI]6364 Å line, suggesting this line should not be affected by the presence of these features (see Fig. 3). The [OI]6364 Å line intensities exhibit significantly less scatter, consistent with the stochastic uncertainties and supporting this conclusion. Therefore for the September and October data we use the [OI]6364 Å line flux and multiply it by the branching ratio of 3.0 to obtain the line flux that is used to calculate HO production rates. This was not required for our November observations, as the Doppler shift had decreased (in absolute value) so that the [OI]6300 Å line was no longer located on an O telluric absorption (see Fig. 3), and spectra obtained on these dates have red line ratios consistent with the known branching ratio of 3.0.

2.4 Spectral Fitting Model for OH

Unlike [OI] emission, the OH emission exhibits ro-vibrational structure superimposed on the electronic transition, which results in an intricate band structure for the emission features (see Fig. 1). To fit the OH emission we employ the same modeling technique used by McKay et al. (2016), which we briefly describe here. The program references a line list for the species of interest (in this case the OH A-X (0-0) band). Once this line list is compiled, we fit each spectral feature in the list to a Gaussian profile. We add all the line profile fits together to create an empirical fit to the spectrum and then integrate over this model fit to obtain a flux for the OH band. While in general the program handles multiple species simultaneously, the OH A-X (0-0) band does not suffer from contamination of spectral features due to other species, therefore their inclusion is not necessary. More details on the general program can be found in McKay et al. (2014), with the specific application to OH detailed in McKay et al. (2016).

2.5 Conversion of Observed Flux to Production Rate

We derived HO production rates based on two different spectral features: [OI]6300, 6364 Å emission and the OH A-X (0-0) band at 3080 Å. For [OI], we employ the same Haser model for the [OI] emission used in McKay et al. (2012) to infer HO production rates from our [OI] observations. For OH, we employ the same Haser model as in McKay et al. (2016), which has modifications that emulate the vectorial model (see McKay et al. (2014) for more details), and adopt fluorescence efficiencies from Schleicher and A’Hearn (1988) that account for the Swings effect (the dependence of the fluorescence efficiency on the heliocentric velocity of the comet). The OH production rates are converted to HO production rates using the relation from Cochran and Schleicher (1993):


where R is heliocentric distance in AU. The relation depends on heliocentric distance because the expansion velocity of HO in the coma is dependent on heliocentric distance while the ejection velocity of OH into the coma is not. Parameters for the Haser models are given in Table 2.

3 Results

The resulting average HO production rates for each observation date are given in Table 3 and plotted in Fig. 4, along with other measured HO production rates for ISON from the literature. In Table 3 we also present the surface area needed to provide the observed HO production (i.e. the active area). We employ the sublimation model of Cowan and A’Hearn (1979) with a visual albedo of 0.05 and thermal emissivity of unity and the slow-rotator approximation. The slow-rotator approximation is appropriate for objects with slow rotation rates or low thermal inertia, and is typically the model used to describe cometary nuclei (e.g. Bodewits et al., 2014). Table 3 also displays the active fraction of the nucleus based on our active area calculations and assuming a spherical nucleus with a radius of 0.75 0.15 km as derived by Lamy et al. (2014) from HST observations (although they quote a more precise value of 0.68 0.02 km, Lamy et al. (2014) also note that this precise number is model dependent, so for this work we will employ their more conservative value of 0.75 0.15 km). Size estimates for Oort Cloud comets like ISON are very rare, so the size estimate for ISON provides a unique opportunity to constrain the active fraction of an Oort Cloud comet. As both active area and active fraction are model dependent, Table 4 shows active areas and active fractions derived for isothermal, subsolar, and fast-rotator models in addition to the slow-rotator model to examine the effect of model assumptions on our results. Table 5 is the same as Table 4, except for a visual albedo of 0.5 instead of 0.05. The active area results are also plotted in Fig. 5. The derived active areas are as high as 100% and in some cases even higher. Implications for this will be discussed in the next section.

4 Discussion

4.1 Evolution of HO Production

Our HO production rates show a very slow increase from 1.8 to 0.9 AU, with a dramatic rise in HO production at around 0.6 AU, in agreement with previous measurements as shown in Fig. 4. Our observation at a heliocentric distance of 0.6 AU was only 12 hours after the first of several outbursts that were observed at heliocentric distances less than 0.6 AU (see Sekanina and Kracht, 2014, for an analysis of comet ISON’s lightcurve and outbursts). The difference in HO production before and after the outburst is about a factor of 15, which matches well the observed increase in brightness during the outburst (Sekanina and Kracht, 2014), suggesting that the outburst was indeed correlated to a massive increase in HO production. We will discuss what led to this strong increase in gas production in the following sections. The HO production rate seems to have leveled out between 0.6 and 0.44 AU, however we note that clouds may have influenced our derived fluxes on this date. Although our derived production rate is consistent with other measurements around this time, the uncertainty introduced by clouds means the uncertainty in this measurement could be a factor of two or maybe even higher, making it difficult to derive any strong conclusions on the behavior of the HO production between 0.6 and 0.44 AU.

4.2 Extended Sources of HO

Past observations of comets such as C/2009 P1 (Garradd) have shown trends in derived HO production rate with the projected field of view at the comet, with larger fields of view showing larger production rates. This has been cited as evidence for an icy grain halo with an extent of 10-10 km that serves as an extended source of HO production (Combi et al., 2013; Bodewits et al., 2014; McKay et al., 2015). The Lyman- observations of Combi et al. (2014a) have very large fields of view ( 10-10 km), while slit-based spectroscopy like the observations presented here and IR observations such as Dello Russo et al. (2016) and DiSanti et al. (2016) have much smaller projected fields (a few 10 km), and narrowband OH imaging (Knight and Schleicher, 2015; Opitom et al., 2014) has intermediate fields of view on the order of tens of thousands of km. However, throughout the apparition for ISON there appears to be no correlation between projected field of view and derived HO production, with the possible exception of our [OI] observations at R=1.8 AU and 1.6 AU and nearly concurrent measurements of OH from Knight and Schleicher (2015), which employ a field of view on the order of 510 km (see Fig. 4). Knight and Schleicher (2015) argue from analysis of the dust spatial profiles that before late October (i.e. at the time of our observations at R=1.8 AU and 1.6 AU) the coma did contain icy grains, so it is possible that there was an extended source of HO production at larger heliocentric distance that dissipated with the disappearance of the icy grains in late October.  Combi et al. (2014a) reached a similar conclusion about the presence of an extended source of HO production in comet ISON from late October onward. It is also possible that the extended source did not disappear, but rather became much smaller in spatial extent and was no longer spatially resolved by observations. With derived HO production rates being consistent across platforms, any extended source present at R 1.3 AU must have had a spatial extent smaller than the projected field of view of our observations ( 1000 km).

4.3 Analysis of the Active Area: Evidence for Fragmentation

4.3.1 Evolution of the Active Area

Our derived active areas (Table 3) show a relatively constant value of  10 km outside of R=1.2 AU, which then decreased to about half this value over the R=0.9-1.2 AU range, in agreement with calculations by Combi et al. (2014a), although Combi et al. (2014a) have a slightly later onset (1 AU) for their tentative detection of a decrease in active area. For the nucleus radius of 0.75 0.15 km determined by Lamy et al. (2014), the active areas derived for R 0.9 AU are consistent with a very high active fraction for the nucleus of 50% or even a completely active surface, though the error bars are quite large due to uncertainty in the actual size of the nucleus coupled with uncertainties in the HO production rate. In addition, our calculations of the active fraction for ISON assume a spherical nucleus, which is not necessarily a valid assumption. Although a precise measurement of the active fraction for ISON is difficult, both our results and those of Combi et al. (2014a) suggest that it is 50-100% before the first outburst at R 0.6 AU. However, at R=0.6 AU, both we and Combi et al. (2014a) find an active area of 30 km, much larger than the total surface area of a spherical nucleus with radius 0.75 km. This implies that at this time there was a large fragmentation event, or at least shedding of a large amount of ice-rich material, that resulted in more surface area being exposed directly to sunlight. There is also the possibility that the active area decreased to half this value by 0.44 AU, but as stated in Section 4.1 the water production rate derived on this date may be affected by clouds, making it difficult to arrive at any definitive conclusions about the active area at this time. Comparison to other measurements of the water production rates in Fig. 2 and the active areas they imply suggest that it seems unlikely that the active area was higher at 0.44 AU than at 0.6 AU.

4.3.2 Alternate Models

To test the robustness of this conclusion, we examine not only active areas from the slow-rotator model, but also isothermal, subsolar, and fast-rotator models. The slow-rotator model is valid for slow rotation periods and/or low thermal inertia of the surface material. These conditions mean that every point on the cometary surface is in thermal equilibrium with the solar radiation incident upon it. Comet rotation periods vary anywhere from hours to days, with a possible value of 10.4 hours determined for ISON (Lamy et al., 2014). This does not fulfill the slow-rotator approximation. However, the thermal inertias for comets where such measurements are available are quite low (e.g. Groussin et al., 2013; Choukroun et al., 2015), hence the tendency in the literature is to assume cometary nuclei follow the slow-rotator approximation. Therefore the slow-rotator is our preferred model. The isothermal model assumes the whole nucleus surface is the same temperature, while the subsolar model assumes the whole nucleus surface has the temperature of the subsolar point. Although these models are idealizations that likely do not occur in reality, they provide bounds on the sublimation rate coming off the surface. The fast-rotator model assumes that lines of latitude on the cometary nucleus are isotherms, which occurs for fast rotation periods and/or high thermal inertias of the surface material.

The active areas and active fractions derived for the different models with a visual albedo of 0.05 are shown in Table 4 and Fig. 5. At the heliocentric distances of our observations ( 2 AU, well within the HO sublimation line, such that HO can be assumed to be fully activated), the fast-rotator and isothermal approximations give very similar active areas, which are slightly higher than those derived from the slow-rotator model. This discrepancy is maximum at the largest heliocentric distances, and decreases as the comet moves towards the Sun. The subsolar point model produces significantly smaller values for the active area than the other models, and is the only model with active fractions less than 50%.

Because all of these models suggest ISON’s active fraction was much higher than the typical 5% seen in other comets (e.g. A’Hearn et al., 1995), this either suggests that a significant fraction of the comet’s surface was active, or there was an extended source of HO in the coma such as sublimation from icy grains. We discussed the possibility of an extended source in the previous section. For the following we will explore the possibility of a surface that was 20-100% active (encompassing values from all the thermal models). This high an active fraction could suggest that a large fraction of the cometary surface was exposed HO ice, meaning assuming a visual albedo of 0.05 may not be a reasonable assumption. Therefore we explored thermal models for the nucleus where both the visual albedo and thermal albedo of the nucleus have a much higher value of 0.50. The results of these models are presented in Table 5 and Fig. 5. Assuming a higher albedo shifts the derived active areas upward by about a factor of two compared to the lower albedo models. This only strengthens the conclusion that ISON’s active fraction was quite large and perhaps an icy grain source of HO is needed to account for the comet being hyperactive, though this source would have to have a small spatial extent (see section 4.2).

We will now examine which of our thermal models is most likely to describe ISON’s activity. The subsolar model is the only model with active fractions less than 50%, meaning this is the only model that does not require a nearly uniformly active surface or an extended source. The subsolar model is only valid if all the activity is dominated by the subsolar point, which is generally not observed in comets. A special case that would fit this description is where a polar source region is the dominant source of activity and the rotation pole is pointed directly at the Sun (i.e. the rotation pole is at the subsolar point). Early in the apparition (Spring 2013, R=4.0 AU) there was evidence for a polar source region and a rotation pole pointed at the Sun based on HST (Li et al., 2013) and ground-based observations (Knight and Schleicher, 2015) of a dust jet pointed in the solar direction, the orientation of which did not change over time. However, this feature dissipated as ISON moved towards the Sun and the CN jet morphology during October-November did not match well with a polar jet (Knight and Schleicher, 2015), making it likely that during the time of our observations the activity was not driven by a polar source region located at the subsolar point. Therefore we do not find the subsolar model a likely description of ISON’s activity. We cannot rule out the plausibility of the other models with current results.

4.3.3 Implications for ISON’s Activity

The arguments of the previous section mean that ISON’s surface was at least 50% active or there was an extended source of HO production in the coma, most likely due to the presence of icy grains. The high albedo (A=0.5) cases all favor an active area of more than 100%, implying an extended source, in which case there is a degeneracy between the active area of the extended source and the active area of the nucleus. The low albedo cases (A=0.05) are consistent with both a highly active nucleus or an extended source. If an extended source is present it must have a small spatial extent (less than 1000 km, the projected field of view of our slit spectroscopy) as argued in section 4.2. This is plausible, considering that 103P/Hartley 2 was shown to have a hyperactive nucleus and that this extra active area came from an extended source of icy grains in the coma (A’Hearn et al., 2011; Kelley et al., 2013). However, HO production rates derived from remote sensing observations agree over drastically different fields of view, from 500 km for slit spectroscopy (Dello Russo et al., 2011; Mumma et al., 2011; McKay et al., 2014) to 10-10 km for imaging of OH (Knight and Schleicher, 2013) and 10 km for Lyman- (Combi et al., 2011), suggesting that the spatial extent of the icy grain source was not much larger than the smallest fields of view used (in this case a few hundered km). This is consistent with the expected lifetimes of micron-sized icy grains at the heliocentric distance of 103P/Hartley 2 at the time of observation (R 1.1 AU) of only 100-1000 seconds (Beer et al., 2006), which is comparable to the crossing time for icy grains moving at 1 km s. Therefore for ISON to exhibit a similar phenomenon at similar and smaller heliocentric distances to 103P/Hartley 2 is certainly plausible, although uncertainties in the active area and nucleus size preclude this from being demonstrated conclusively. Using observations of the spatial distribution of HO and rotational temperature in the inner coma (nucleocentric distance 1000 km), Bonev et al. (2014) argued that ISON was releasing icy grains at heliocentric distances less than 0.6 AU, but as they obtained no observations at larger heliocentric distances (in particular before the first outburst in mid-November), their observations do not place a constraint on an icy grain source before the outburst.

All models for ISON on November 15 display active fractions greater than unity (Fig. 5, Tables 4 and 5), implying a large increase in the active area of the comet/extended coma source. This could be due to either a catastrophic fragmentation of the nucleus or a large outburst ejecting icy grains into the coma. As discussed above, Bonev et al. (2014) suggested that after the first outburst at 0.6 AU sublimation from icy grains was a significant source of water production. However,  Steckloff et al. (2015) argued that the coma morphology during and after the outburst in mid-November is inconsistent with a simple explosive event such as amorphous-crystalline water ice transition or trapped hypervolatiles triggering the outbursts and ejecting icy grains into the coma. They argued instead that ISON split into a swarm of fragments with a characteristic size of 100 m. As they are all the same size, drifting between fragments would not show up in Earth-based observations of the coma morphology until days after the event, consistent with observations. However, it is likely that both the fragmentation event and release of small (micron-sized) icy grains contributed to the increased water production.

Although the absolute values of the active area and active fraction are model dependent, the evolution of the active area over the course of the apparition is much less sensitive to the specific thermal model employed, as shown in Fig. 5. All models support a relatively constant active area during October (R 1.2 AU), a decrease in late October and early November, and a massive increase in active area in mid-November. This is very similar to the results presented in Combi et al. (2014a). The massive increase in active area in mid-November was likely associated with the first of several fragmentation events (Steckloff et al., 2015), which increased the sublimating surface area in the form of smaller fragments and may also have released micron-sized icy grains into the coma, further increasing the effective active area.

It is also possible in mid-November that both thermal emission from sub-micron grains in the coma and reflected solar radiation from the dust coma could provide significant heating to the surface, meaning our simple thermodynamical models for the nucleus temperature and sublimation would not be accurate. If this is the case, this would mean less active area is required to explain the observed HO production rates. However, quantifying this effect requires detailed calculations that are beyond the scope of this work.

4.4 Influence of O

Recently, the Rosetta mission discovered the presence of O in the coma of 67P/Churyumov-Gerasimenko at an abundance of 4% relative to H(Bieler et al., 2015). This discovery was unexpected and means that O photodissociation could be an unaccounted source of OI in cometary comae. Specifically of relevance to this work is that O has a high branching ratio for releasing OI in the D state upon photodissociation (0.85, Huebner et al., 1992), which will then radiatively decay and release [OI]6300 Å line emission. This could potentially influence attempts to use [OI]6300 Å line emission as a proxy for HO production. To test the potential influence of O on [OI]6300 Å line intensities we modified our Haser model for [OI] to have O be the parent of OI instead of HO. For this purpose we employed a branching ratio of 0.854 for photodissociation of O O(P) + O(D) and a photodestruction timescale of 2.0810 s at 1 AU (Huebner et al., 1992). We find that at the abundance of 4% relative to HO observed by Rosetta at 67P by the ROSINA instrument (Bieler et al., 2015), O can only account for 5-10% of the observed [OI]6300 line flux. As our uncertainties are 10-20%, O at an abundance of 4% relative to HO does not significantly influence our derived HO production rates. Additionally, the consistency between our production rates derived from [OI] and OH, as well as agreement with other values already in the literature from methods other than [OI], support the conclusion that any O present in the coma of ISON did not influence our derived HO production rates. Although this suggests [OI] studies cannot be used to infer O abundances from the ground (which is very challenging to do directly due to severe telluric absorption), it does show that [OI]6300 Å line emission can still be used as a proxy for HO in comets for O/HO ratios similar to that measured by ROSINA. This conclusion is valid for the small FOV employed for our observations, with the potential contribution of O becoming more important for larger FOV. Using Rosetta Alice UV stellar occultation data Keeney et al. (2017) found higher O/HO ratios than those measured by ROSINA (some as high as 50%). For O/HO of 50% we expect that O photodissociation would be a significant contributor to the [OI] emission, being responsible for approximately half of the emission. It is possible that such a high O abundance could artificially raise our derived HO production rates so that they agree with observations such as Combi et al. (2014b) with larger FOV’s, and could hide the presence of an extended source. However, the consistency between our OH and [OI] derived HO production rates (which have similar FOV) argues against this scenario and therefore against such a high O abundance for comet ISON.

5 Conclusions

In this work we presented HO production rates for comet C/2012 S1 (ISON) derived from observations of [OI] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8-0.44 AU. Our production rates are in agreement with previous measurements using a variety of different techniques. As our [OI]-derived HO production rates are consistent with other methods, we conclude that the contribution of O photodissociation to the [OI] emission is negligible compared to the contribution of HO (for abundances observed by the ROSINA instrument aboard Rosetta), meaning [OI] emission can still be utilized as a reliable proxy for HO production in comets (though for the much higher abundances observed by Rosetta Alice the presence of O could introduce systematic error). The lack of any apparent dependence of derived production rate with instrument field of view suggests that any extended source of HO production must have had a spatial extent smaller than the smallest projected field of view observation (in this case 1,000 km). We find that ISON had an evolving active area throughout the apparition, with a rapid increase at about R=0.6 AU corresponding to the first of three major outbursts ISON experienced inside of 1 AU, consistent with other studies. This rapid increase in active area likely indicates a major fragmentation event of the nucleus. The 5-10 km active area observed outside of R=0.6 AU is consistent with a 50-100% active fraction for the nucleus, much higher than typically observed for cometary nuclei. Although the absolute value of the active area is somewhat dependent on the thermal model employed, the changes in observed active area are consistent among models, and the conclusion of a 50-100+% active fraction is robust for realistic thermal models of the nucleus. However, the possibility of a spatially unresolved icy grain source cannot be discounted.

We are grateful to Lori Feaga and an anonymous reviewer for helpful comments that improved the quality of this manuscript. We thank John Barentine, Jurek Krzesinski, Chris Churchill, Pey Lian Lim, Paul Strycker, and Doug Hoffman for developing and optimizing the ARCES IRAF reduction script used to reduce these data. We thank the APO and Keck observing staffs for their invaluable help in conducting the observations. We are extremely grateful to Dr. Zlatan Tsvetanov for yielding some of his observing time so we could obtain the UT November 20 observations presented in this paper. We thank Dr. Cyrielle Opitom for sharing some of her unpublished water production rates for inclusion in Fig. 4. We would also like to acknowledge the JPL Horizons System, which was used to generate ephemerides for nonsidereal tracking of the comets during the observations, and the SIMBAD database, which was used for selection of reference stars. This work was supported by the NASA Postdoctoral Program, administered by USRA, as well as the NASA Planetary Atmospheres Program, the NASA Planetary Astronomy Program, and the NASA Emerging Worlds Program. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.


  • A’Hearn et al. (1995) M. F. A’Hearn, R. L. Millis, D. G. Schleicher, D. J. Osip, and P. V. Birch. The ensemble properties of comets: Results from narrowband photometry of 85 comets, 1976-1992. Icarus, 118:223–270, December 1995. doi: 10.1006/icar.1995.1190.
  • A’Hearn et al. (2011) M. F. A’Hearn, M. J. S. Belton, W. A. Delamere, L. M. Feaga, D. Hampton, J. Kissel, K. P. Klaasen, L. A. McFadden, K. J. Meech, H. J. Melosh, P. H. Schultz, J. M. Sunshine, P. C. Thomas, J. Veverka, D. D. Wellnitz, D. K. Yeomans, S. Besse, D. Bodewits, T. J. Bowling, B. T. Carcich, S. M. Collins, T. L. Farnham, O. Groussin, B. Hermalyn, M. S. Kelley, M. S. Kelley, J.-Y. Li, D. J. Lindler, C. M. Lisse, S. A. McLaughlin, F. Merlin, S. Protopapa, J. E. Richardson, and J. L. Williams. EPOXI at Comet Hartley 2. Science, 332:1396–, June 2011. doi: 10.1126/science.1204054.
  • Beer et al. (2006) E. H. Beer, M. Podolak, and D. Prialnik. The contribution of icy grains to the activity of comets. I. Grain lifetime and distribution. Icarus, 180:473–486, February 2006. doi: 10.1016/j.icarus.2005.10.018.
  • Bhardwaj and Raghuram (2012) A. Bhardwaj and S. Raghuram. A Coupled Chemistry-emission Model for Atomic Oxygen Green and Red-doublet Emissions in the Comet C/1996 B2 Hyakutake. Astrophysical Journal, 748:13, March 2012. doi: 10.1088/0004-637X/748/1/13.
  • Bieler et al. (2015) A. Bieler, K. Altwegg, H. Balsiger, A. Bar-Nun, J.-J. Berthelier, P. Bochsler, C. Briois, U. Calmonte, M. Combi, J. de Keyser, E. F. van Dishoeck, B. Fiethe, S. A. Fuselier, S. Gasc, T. I. Gombosi, K. C. Hansen, M. Hässig, A. Jäckel, E. Kopp, A. Korth, L. Le Roy, U. Mall, R. Maggiolo, B. Marty, O. Mousis, T. Owen, H. Rème, M. Rubin, T. Sémon, C.-Y. Tzou, J. H. Waite, C. Walsh, and P. Wurz. Abundant molecular oxygen in the coma of comet 67P/Churyumov-Gerasimenko. Nature, 526:678–681, October 2015. doi: 10.1038/nature15707.
  • Biver et al. (2002) N. Biver, D. Bockelée-Morvan, P. Colom, J. Crovisier, F. Henry, E. Lellouch, A. Winnberg, L.E.B. Johansson, M. Gunnarsson, H. Rickman, F. Rantakyrö, Davies. J.K., W.R.F. Dent, G. Paubert, R. Moreno, J. Wink, D. Despois, D. Benford, M. Gardner, D.C. Lis, D. Mehringer, T.G. Phillips, and H. Rauer. The 1995-2002 long-term monitoring of comet C/1995 O1 (Hale-Bopp) at radio wavelength. Earth, Moon, and Planets, 90(1):5–14, 2002.
  • Bodewits et al. (2014) D. Bodewits, T. L. Farnham, M. F. A’Hearn, L. M. Feaga, A. McKay, D. G. Schleicher, and J. M. Sunshine. The Evolving Activity of the Dynamically Young Comet C/2009 P1 (Garradd). Astrophysical Journal, 786:48, May 2014. doi: 10.1088/0004-637X/786/1/48.
  • Bonev et al. (2006) B. P. Bonev, M. J. Mumma, M. A. DiSanti, N. Dello Russo, K. Magee-Sauer, R. S. Ellis, and D. P. Stark. A Comprehensive Study of Infrared OH Prompt Emission in Two Comets. I. Observations and Effective g-Factors. Astrophysical Journal, 653:774–787, December 2006. doi: 10.1086/508452.
  • Bonev et al. (2014) B. P. Bonev, M. A. DiSanti, G. L. Villanueva, E. L. Gibb, L. Paganini, and M. J. Mumma. The Inner Coma of Comet C/2012 S1 (ISON) at 0.53 AU and 0.35 AU from the Sun. Astrophysical Journal Letters, 796:L6, November 2014. doi: 10.1088/2041-8205/796/1/L6.
  • Choukroun et al. (2015) M. Choukroun, S. Keihm, F. P. Schloerb, S. Gulkis, E. Lellouch, C. Leyrat, P. von Allmen, N. Biver, D. Bockelée-Morvan, J. Crovisier, P. Encrenaz, P. Hartogh, M. Hofstadter, W.-H. Ip, C. Jarchow, M. Janssen, S. Lee, L. Rezac, G. Beaudin, B. Gaskell, L. Jorda, H. U. Keller, and H. Sierks. Dark side of comet 67P/Churyumov-Gerasimenko in Aug.-Oct. 2014. MIRO/Rosetta continuum observations of polar night in the southern regions. Astronomy and Astrophysics, 583:A28, November 2015. doi: 10.1051/0004-6361/201526181.
  • Cochran and Schleicher (1993) A. L. Cochran and D. G. Schleicher. Observational Constraints on the Lifetime of Cometary H O. Icarus, 105:235–253, September 1993. doi: 10.1006/icar.1993.1121.
  • Cochran and Cochran (2001) A.L. Cochran and W.D. Cochran. Observations of O (S) and O (D) in spectra of C/1999 S4 (LINEAR). Icarus, 154(2):381–390, 2001.
  • Cochran and Cochran (2002) A.L. Cochran and W.D. Cochran. A high spectral resolution atlas of comet 122P/de Vico. Icarus, 157(2):297–308, 2002.
  • Combi et al. (2011) M. R. Combi, J.-L. Bertaux, E. Quémerais, S. Ferron, and J. T. T. Mäkinen. Water Production by Comet 103P/Hartley 2 Observed with the SWAN Instrument on the SOHO Spacecraft. Astrophysical Journal Letters, 734:L6, June 2011. doi: 10.1088/2041-8205/734/1/L6.
  • Combi et al. (2013) M. R. Combi, J. T. T. Mäkinen, J.-L. Bertaux, E. Quémerais, S. Ferron, and N. Fougere. Water production rate of Comet C/2009 P1 (Garradd) throughout the 2011-2012 apparition: Evidence for an icy grain halo. Icarus, 225:740–748, July 2013. doi: 10.1016/j.icarus.2013.04.030.
  • Combi et al. (2014a) M. R. Combi, N. Fougere, J. T. T. Mäkinen, J.-L. Bertaux, E. Quémerais, and S. Ferron. Unusual Water Production Activity of Comet C/2012 S1 (ISON): Outbursts and Continuous Fragmentation. Astrophysical Journal Letters, 788:L7, June 2014a. doi: 10.1088/2041-8205/788/1/L7.
  • Combi et al. (2014b) M. R. Combi, J. T. Mäkinen, J. L. Bertaux, E. Quémerais, and S. Ferron. The Water Production Rate of Recent Comets (2013-2014) by SOHO/SWAN: 2P/Encke (2013), C/2013 R1 (Lovejoy), and C/2013 A1 (Siding Spring). In AAS/Division for Planetary Sciences Meeting Abstracts, volume 46 of AAS/Division for Planetary Sciences Meeting Abstracts, page 110.09, November 2014b.
  • Cowan and A’Hearn (1979) J. J. Cowan and M. F. A’Hearn. Vaporization of comet nuclei - Light curves and life times. Moon and Planets, 21:155–171, October 1979. doi: 10.1007/BF00897085.
  • Dello Russo et al. (2011) N. Dello Russo, R. J. Vervack, Jr., C. M. Lisse, H. A. Weaver, H. Kawakita, H. Kobayashi, A. L. Cochran, W. M. Harris, A. J. McKay, N. Biver, D. Bockelée-Morvan, and J. Crovisier. The Volatile Composition and Activity of Comet 103P/Hartley 2 During the EPOXI Closest Approach. Astrophysical Journal Letters, 734:L8+, June 2011. doi: 10.1088/2041-8205/734/1/L8.
  • Dello Russo et al. (2016) N. Dello Russo, R. J. Vervack, H. Kawakita, A. Cochran, A. J. McKay, W. M. Harris, H. A. Weaver, C. M. Lisse, M. A. DiSanti, H. Kobayashi, N. Biver, D. Bockelée-Morvan, J. Crovisier, C. Opitom, and E. Jehin. The compositional evolution of C/2012 S1 (ISON) from ground-based high-resolution infrared spectroscopy as part of a worldwide observing campaign. Icarus, 266:152–172, March 2016. doi: 10.1016/j.icarus.2015.11.030.
  • DiSanti et al. (2016) M. A. DiSanti, B. P. Bonev, E. L. Gibb, L. Paganini, G. L. Villanueva, M. J. Mumma, J. V. Keane, G. A. Blake, N. Dello Russo, K. J. Meech, R. J. Vervack, Jr., and A. J. McKay. En Route to Destruction: The Evolution in Composition of Ices in Comet D/2012 S1 (ISON) between 1.2 and 0.34 AU from the Sun as Revealed at Infrared Wavelengths. Astrophysical Journal, 820:34, March 2016. doi: 10.3847/0004-637X/820/1/34.
  • Festou and Feldman (1981) M.C. Festou and P.D. Feldman. The forbidden oxygen lines in comets. Astronomy and Astrophysics, 103(1):154–159, 1981.
  • Filippenko (1982) A. V. Filippenko. The importance of atmospheric differential refraction in spectrophotometry. Publications of the Astronomical Society of the Pacific, 94:715–721, August 1982. doi: 10.1086/131052.
  • Fink (2009) U. Fink. A taxonomic survey of comet composition 1985-2004 using CCD spectroscopy. Icarus, 201:311–334, May 2009. doi: 10.1016/j.icarus.2008.12.044.
  • Groussin et al. (2013) O. Groussin, J. M. Sunshine, L. M. Feaga, L. Jorda, P. C. Thomas, J.-Y. Li, M. F. A’Hearn, M. J. S. Belton, S. Besse, B. Carcich, T. L. Farnham, D. Hampton, K. Klaasen, C. Lisse, F. Merlin, and S. Protopapa. The temperature, thermal inertia, roughness and color of the nuclei of Comets 103P/Hartley 2 and 9P/Tempel 1. Icarus, 222:580–594, February 2013. doi: 10.1016/j.icarus.2012.10.003.
  • Hogg et al. (2001) D. W. Hogg, D. P. Finkbeiner, D. J. Schlegel, and J. E. Gunn. A Photometricity and Extinction Monitor at the Apache Point Observatory. Astronomical Journal, 122:2129–2138, October 2001. doi: 10.1086/323103.
  • Huebner et al. (1992) W.F. Huebner, J.J. Keady, and S.P. Lyon. Solar photo rates for planetary atmospheres and atmospheric pollutants. Astrophysics and Space Science, 195(1):1–289, 291–294, 1992.
  • Keeney et al. (2017) B. A. Keeney, S. A. Stern, M. F. A’Hearn, J.-L. Bertaux, L. M. Feaga, P. D. Feldman, R. A. Medina, J. W. Parker, J. P. Pineau, E. Schindhelm, A. J. Steffl, M. Versteeg, and H. A. Weaver. HO and O absorption in the coma of comet 67P/Churyumov-Gerasimenko measured by the Alice far-ultraviolet spectrograph on Rosetta. Monthly Notices of the Royal Astronomical Society, 469:S158–S177, July 2017. doi: 10.1093/mnras/stx1426.
  • Kelley et al. (2013) M. S. Kelley, D. J. Lindler, D. Bodewits, M. F. A’Hearn, C. M. Lisse, L. Kolokolova, J. Kissel, and B. Hermalyn. A distribution of large particles in the coma of Comet 103P/Hartley 2. Icarus, 222:634–652, February 2013. doi: 10.1016/j.icarus.2012.09.037.
  • Knight and Schleicher (2013) M. M. Knight and D. G. Schleicher. The highly unusual outgassing of Comet 103P/Hartley 2 from narrowband photometry and imaging of the coma. Icarus, 222:691–706, February 2013. doi: 10.1016/j.icarus.2012.06.004.
  • Knight and Schleicher (2015) M. M. Knight and D. G. Schleicher. Observations of Comet ISON (C/2012 S1) from Lowell Observatory. Astronomical Journal, 149:19, January 2015. doi: 10.1088/0004-6256/149/1/19.
  • Lamy et al. (2014) P. L. Lamy, I. Toth, and H. A. Weaver. Hubble Space Telescope Observations of the Nucleus of Comet C/2012 S1 (ISON). Astrophysical Journal Letters, 794:L9, October 2014. doi: 10.1088/2041-8205/794/1/L9.
  • Li et al. (2013) J.-Y. Li, M. S. P. Kelley, M. M. Knight, T. L. Farnham, H. A. Weaver, M. F. A’Hearn, M. J. Mutchler, L. Kolokolova, P. Lamy, and I. Toth. Characterizing the Dust Coma of Comet C/2012 S1 (ISON) at 4.15 AU from the Sun. Astrophysical Journal Letters, 779:L3, December 2013. doi: 10.1088/2041-8205/779/1/L3.
  • McKay et al. (2012) A. J. McKay, N. J. Chanover, J. P. Morgenthaler, A. L. Cochran, W. M. Harris, and N. D. Russo. Forbidden oxygen lines in Comets C/2006 W3 Christensen and C/2007 Q3 Siding Spring at large heliocentric distance: Implications for the sublimation of volatile ices. Icarus, 220:277–285, July 2012. doi: 10.1016/j.icarus.2012.04.030.
  • McKay et al. (2014) A. J. McKay, N. J. Chanover, M. A. DiSanti, J. P. Morgenthaler, A. L. Cochran, W. M. Harris, and N. D. Russo. Rotational variation of daughter species production rates in Comet 103P/Hartley: Implications for the progeny of daughter species and the degree of chemical heterogeneity. Icarus, 231:193–205, March 2014. doi: 10.1016/j.icarus.2013.11.029.
  • McKay et al. (2015) A. J. McKay, A. L. Cochran, M. A. DiSanti, G. Villanueva, N. D. Russo, R. J. Vervack, J. P. Morgenthaler, W. M. Harris, and N. J. Chanover. Evolution of HO, CO, and CO production in Comet C/2009 P1 Garradd during the 2011-2012 apparition. Icarus, 250:504–515, April 2015. doi: 10.1016/j.icarus.2014.12.023.
  • McKay et al. (2016) A. J. McKay, M. S. P. Kelley, A. L. Cochran, D. Bodewits, M. A. DiSanti, N. D. Russo, and C. M. Lisse. The CO abundance in Comets C/2012 K1 (PanSTARRS), C/2012 K5 (LINEAR), and 290P/Jäger as measured with Spitzer. Icarus, 266:249–260, March 2016. doi: 10.1016/j.icarus.2015.11.004.
  • Morgenthaler et al. (2001) J.P. Morgenthaler, W.M. Harris, F. Scherb, C. Anderson, R.J. Oliversen, N.E. Doane, M.R. Combi, M.L. Marconi, and W.H. Smyth. Large aperture [O I] 6300 photometry of comet Hale-Bopp: Implications for the photochemistry of OH. Astrophysical Journal, 563(1):451–461, 2001.
  • Morgenthaler et al. (2007) J.P. Morgenthaler, W.M. Harris, and M.R. Combi. Large Aperture O I 6300 Å Observations of Comet Hyakutake: Implications for the Photochemistry of OH and O I Production in Comet Hale-Bopp. Astrophysical Journal, 657:1162–1171, March 2007. doi: 10.1086/511062.
  • Mumma et al. (2011) M. J. Mumma, B. P. Bonev, G. L. Villanueva, L. Paganini, M. A. DiSanti, E. L. Gibb, J. V. Keane, K. J. Meech, G. A. Blake, R. S. Ellis, M. Lippi, H. Boehnhardt, and K. Magee-Sauer. Temporal and Spatial Aspects of Gas Release During the 2010 Apparition of Comet 103P/Hartley 2. Astrophysical Journal Letters, 734:L7, June 2011. doi: 10.1088/2041-8205/734/1/L7.
  • Novski et al. (2012) V. Novski, A. Novichonok, O. Burhonov, W. H. Ryan, E. V. Ryan, H. Sato, E. Guido, G. Sostero, N. Howes, and G. V. Williams. Comet C/2012 S1 (Ison). Central Bureau Electronic Telegrams, 3238:1, September 2012.
  • Opitom et al. (2014) C. Opitom, E. Jehin, J. Manfroid, D. Hutsemékers, and M. Gillon. TRAPPIST monitoring of comets C/2012 S1 (ISON) and C/2013 R1 (Lovejoy). In K. Muinonen, A. Penttilä, M. Granvik, A. Virkki, G. Fedorets, O. Wilkman, and T. Kohout, editors, Asteroids, Comets, Meteors 2014, July 2014.
  • Opitom et al. (2015) C. Opitom, E. Jehin, J. Manfroid, D. Hutsemékers, M. Gillon, and P. Magain. TRAPPIST photometry and imaging monitoring of comet C/2013 R1 (Lovejoy): Implications for the origin of daughter species. Astronomy and Astrophysics, 584:A121, December 2015. doi: 10.1051/0004-6361/201526427.
  • Schleicher and A’Hearn (1988) D. G. Schleicher and M. F. A’Hearn. The fluorescence of cometary OH. Astrophysical Journal, 331:1058–1077, August 1988. doi: 10.1086/166622.
  • Sekanina and Kracht (2014) Z. Sekanina and R. Kracht. Disintegration of Comet C/2012 S1 (ISON) Shortly Before Perihelion: Evidence from Independent Data Sets. ArXiv e-prints, August 2014.
  • Sharpee and Slanger (2006) BD Sharpee and TG Slanger. O (1d2-3p2, 1, 0) 630.0, 636.4, and 639.2 nm forbidden emission line intensity ratios measured in the terrestrial nightglow. The Journal of Physical Chemistry A, 110(21):6707–6710, 2006.
  • Steckloff et al. (2015) J. K. Steckloff, B. C. Johnson, T. Bowling, H. Jay Melosh, D. Minton, C. M. Lisse, and K. Battams. Dynamic sublimation pressure and the catastrophic breakup of Comet ISON. Icarus, 258:430–437, September 2015. doi: 10.1016/j.icarus.2015.06.032.
  • Wang et al. (2003) S.-i. Wang, R. H. Hildebrand, L. M. Hobbs, S. J. Heimsath, G. Kelderhouse, R. F. Loewenstein, S. Lucero, C. M. Rockosi, D. Sandford, J. L. Sundwall, J. A. Thorburn, and D. G. York. ARCES: an echelle spectrograph for the Astrophysical Research Consortium (ARC) 3.5m telescope. In M. Iye and A. F. M. Moorwood, editors, Instrument Design and Performance for Optical/Infrared Ground-based Telescopes, volume 4841 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, pages 1145–1156, March 2003. doi: 10.1117/12.461447.

Observation Log

Date (UT) Instrument (AU) (km s) (AU) (km s) FOV (km) G2V A0V Flux Cal Conditions
9/23/2013 ARCES 1.80 -31.3 2.38 -48.9 2760 5540 HD30854 HR3224 HR 3454 Clear
10/3/2013 ARCES 1.61 -33.0 2.09 -50.9 2430 4860 HR 79078 HR 3711 HR 3454 Clear
10/18/2013 Tull Coude 1.31 -36.6 1.64 -51.9 1430 9780 Solar Port Gem Gem Clear
10/20/2013 Tull Coude 1.27 -37.2 1.58 -51.8 1380 9420 Solar Port Gem Gem Clear
10/21/2013 Tull Coude 1.25 -37.5 1.55 -51.7 1350 9240 Solar Port Gem Gem Clear
10/25/2013 HIRESb 1.15 -39.1 1.42 -50.8 890 7230 Hyades 64 HR 2207 Hilt 600 Clear
10/28/2013 HIRESb 1.08 -40.3 1.33 -49.8 830 6770 Hyades 64 HR 3134 Hilt 600 Clear
11/6/2013 ARCES 0.87 -44.8 1.10 -43.5 1280 2560 HD 95868 HR 4464 HR 4468 Clear
11/15/2013 ARCES 0.60 -53.5 0.90 -25.7 1050 2100 HD 110747 HR 4722 HR 4963 Clear
11/20/2013 ARCES 0.44 -62.8 0.86 -4.1 1000 2000 HD 110747 HR 4722 HR 4963 Passing Clouds

Flux Standard taken from UT October 3 observation
Flux Standard taken from UT November 15 observation

Table 1:

Parameter Values Used in the Haser Models

Molecule (s) (s) (km s) g-factor (ergs s molecule)
OH 1.3 10 2.1 10 1.02 3.5 10
OI 8.3 10 - - -
OI 1.3 10 - - -

Given for =1 AU. The listed g-factor also accounts for the Swings effect.
For [OI] from dissociation of HO into H and O; branching ratio employed is 0.07 [Bhardwaj and Raghuram, 2012]
For [OI] from dissociation of OH; branching ratio for HO to OH + H employed is 0.855 [Huebner et al., 1992] and the branching ratio for OH to O + H is 0.094 [Bhardwaj and Raghuram, 2012].

Table 2:

Production Rates

UT Date R (AU) (10 mol s) Active Area (km) Active Fraction (%)
9/23/2013 1.80 7.65 1.03 9.0 1.2 127 54
10/3/2013 1.61 9.10 1.15 7.9 1.0 112 47
10/18/2013 1.31 14.5 3.5 7.5 1.8 106 49
10/20/2013 1.27 21.8 4.9 10.5 2.4 149 68
10/21/2013 1.25 17.3 3.9 8.0 1.8 113 52
10/25/2013 1.15 12.7 2.4 4.8 0.9 68 30
10/28/2013 1.08 14.1 2.4 4.7 0.8 66 29
11/6/2013 0.87 21.3 4.3 4.3 1.0 61 28
11/15/2013 0.60 347 29.1 31.8 2.7 450 180
11/20/2013 0.44 400 29.0 19.2 1.4 270 110

Assuming an albedo of 0.05, thermal emissivity of unity, and a slow-rotator thermal model [Cowan and A’Hearn, 1979]. See Tables 4 and 5 for different parameter values. Assuming a spherical nucleus with radius of 0.75 0.15 km [Lamy et al., 2014].

Table 3:

Active Area (km) and Active Fraction (%) for 0.05 Albedo

Model Isothermal Fast-Rotator Slow-Rotator Subsolar
UT Date Active Area Active Fraction Active Area Active Fraction Active Area Active Fraction Active Area Active Fraction
9/22/2013 16.2 2.2 230 97 15.3 2.1 217 92 9.0 1.2 127 54 1.9 0.3 27 11
10/3/2013 12.2 1.5 173 73 11.8 1.5 168 70 7.9 1.0 112 47 1.7 0.2 24 10
10/18/2013 9.8 2.4 139 65 9.7 2.3 137 64 7.5 1.8 106 50 1.7 0.4 24 11
10/20/2013 13.5 3.0 190 87 13.3 3.0 188 86 10.5 2.4 148 68 2.4 0.5 34 16
10/21/2013 10.2 2.3 144 66 10.1 2.3 143 65 8.0 1.8 113 52 1.8 0.4 26 12
10/25/2013 5.9 1.1 84 37 5.9 1.1 83 37 4.8 0.9 69 30 1.1 0.2 16 7
10/28/2013 5.5 0.9 78 34 5.5 0.9 78 34 4.7 0.8 66 29 1.1 0.2 15 7
11/6/2013 4.8 1.0 68 31 4.8 1.0 68 31 4.3 0.9 61 27 1.0 0.2 15 7
11/15/2013 33.7 2.8 480 200 33.7 2.8 480 200 31.8 2.7 450 180 7.7 60 110 40
11/20/2013 19.8 1.4 280 110 19.9 1.4 280 110 19.2 1.4 270 110 4.7 30 67 27
Table 4:

Active Area (km) and Active Fraction (%) for 0.50 Albedo

Model Isothermal Fast-Rotator Slow-Rotator Subsolar
UT Date Active Area Active Fraction Active Area Active Fraction Active Area Active Fraction Active Area Active Fraction
9/22/2013 26.3 3.5 370 160 25.2 3.4 360 150 16.2 2.2 230 100 3.5 0.5 49 21
10/3/2013 20.7 2.6 290 120 20.3 2.6 290 120 14.4 1.8 200 90 3.2 0.4 45 19
10/18/2013 17.4 4.2 250 120 17.2 4.2 240 110 13.9 3.4 200 90 3.2 0.8 45 21
10/20/2013 23.9 5.4 340 160 23.7 5.3 340 150 19.5 4.4 280 130 4.5 1.0 63 29
10/21/2013 18.2 4.1 260 120 18.0 4.1 260 120 14.9 3.4 210 100 3.4 0.8 49 22
10/25/2013 10.7 2.0 150 70 10.6 2.0 150 70 9.0 1.7 130 60 2.1 0.4 30 13
10/28/2013 10.1 1.7 140 60 10.0 1.7 140 60 8.7 1.5 120 50 2.0 0.3 29 13
11/6/2013 8.9 1.8 130 60 8.9 1.8 130 60 8.1 1.6 120 50 2.0 0.4 28 12
11/15/2013 63.1 5.3 890 370 63.1 5.3 890 370 60.0 5.0 850 350 14.7 1.2 210 90
11/20/2013 37.4 2.7 530 210 37.4 2.7 530 210 36.4 2.6 510 210 9.0 0.7 130 50
Table 5:

Figure Captions

Fig 1: Top: Spectrum showing the [OI] line region on November 15. The [OI]6300 Å and [OI]6364 Å lines are labeled. Most other emission features present are due to NH. Bottom: Spectrum showing the OH A-X band on October 25. Error bars have been omitted from both plots for clarity.

Fig 2: Top: Spectrum showing the [OI]6300 Å line on November 15. The telluric line is the weaker feature redward of the cometary line. Bottom: Spectrum showing the [OI]6364 Å line on November 15. The telluric line is the weaker feature redward of the cometary line.

Fig 3: Raw extracted spectra (i.e. no telluric removal, solar subtraction, or flux calibration) showing the [OI]6300 Å line (top row) and the [OI]6364 Å line (bottom row) on UT October 3 (left column) and UT November 15 (right column), along with raw extracted spectra of the telluric standard on each date. The telluric standard has been shifted down and for October 3 the counts were scaled by a factor of 10 to facilitate comparison with the comet spectrum. On October 3 the [OI]6300 Å line is directly coincident with an O telluric absorption feature, while by November 15 the cometary line is no longer coincident with any O telluric absorption. The placement of the cometary line relative to the O absorption is similar for all our September and October observations. There are no strong telluric features underneath the [OI]6364 Å line on either date. Therefore for our September and October observations we calculate the HO production rate using the the [OI]6364 Å line rather than the [OI]6300 Å line.

Fig 4: Plot of HO production in comet ISON as a function of heliocentric distance. Production rates from this work are plotted as filled symbols, with circles denoting production rates based on [OI] observations and squares denoting production rates based on OH. Other results from the literature are plotted for comparison as empty symbols. Our values are in agreement with other measurements, and extend coverage of ISON’s HO production rate to larger heliocentric distances than many previously published results.

Fig 5: Derived active area as a function of heliocentric distance for different thermal models and albedos. The top panel shows results for an albedo of 0.05, while the bottom panel shows results for a much larger albedo of 0.50. For both albedos the isothermal and fast-rotator models give nearly identical values, so for clarity the isothermal model has not been plotted. The dashed line is the total surface area of a spherical nucleus with a radius of 0.75 km [Lamy et al., 2014] and the shaded area depicts the uncertainty in this measurement. Therefore active areas that fall in the shaded region are consistent with an active fraction of 100%, while those above the shaded region indicate ISON is hyperactive, and additional surface area besides that of the nucleus is required to explain the observed HO production rates. Changing the albedo and employing different thermal models changes the absolute value of the derived active area, but trends over the apparition are independent of the model employed. The fast-rotator and slow-rotator models suggest an active fraction of at least 50%, which could be indicative of a largely active surface or the presence of icy grains in the coma. The subsolar model gives much lower active areas, but we believe this model does not provide a realistic description of ISON’s activity (see Section 4.3.2).

Figure 1:
Figure 2:
Figure 3:
Figure 4:
Figure 5:
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description