We present analysis of {\it Hubble Space Telescope} (HST) observations of Centaur 29P/Schwassmann-Wachmann 1 (SW1) while in outburst to characterize the outburst coma and place constraints on the nucleus' spin state. The observations consist of Wide Field and Planetary Camera 2 (WFPC2) images from Cycle 5, GO-5829 (Feldman (1995)) acquired on UT 1996 March 11.3 and 12.1, which serendipitously imaged the Centaur shortly after a major outburst. A multi-component coma was detected consisting of: an expanding outburst dust coma with complex morphology possessing an east-west asymmetry and north-south symmetry contained within 5$''$ ($\sim$19,000 km) of the nucleus, the residual dust shell of an earlier UT 1996 February outburst, and a nearly circular underlying quiescent activity level coma detectable to $\sim$70$''$ ($\sim$267,000 km) away from the nucleus. Photometry of the calibrated WFPC2 images resulted in a measured 5$''$ radius aperture equivalent R-band magnitude of 12.86 $\pm$ 0.02 and an estimated (2.79 $\pm$ 0.05)$\times 10^8$ kg for the lower limit of dust material emitted during the outburst. No appreciable evolution of morphologic features, indicating signatures of nucleus rotation, were detected between the two imaging epochs. The observations were modeled using a 3-D Monte Carlo coma model (Samarasinha (2000)) to place constraints on the nucleus' rotation state. Modeling indicated the morphology is representative of a non-isotropic ejection of dust emitted during a single outburst event with a duration on the order of hours from a single source region corresponding to $\sim$1\% of the surface area. A spin period with lower limit on the order of days is suggested to reproduce the coma morphology seen in the observations.

Samarasinha & Mueller (2013) related changes of cometary rotation to other physical parameters for four Jupiter family comets defining a parameter $X$, which is approximately constant within a factor of two irrespective of the active fraction of a comet. Two additional comets are added to this sample in this paper and the claim of a nearly constant parameter $X$ for these six comets is confirmed, albeit with a larger scatter. Taking the geometric mean of $X$ for all the comets above excluding 2P/Encke (as $X$ for each comet was determined with respect to that of 2P/Encke), the expected changes in the rotation periods for a sample of 24 periodic comets are derived. We identify comets from this sample that are most likely to show observationally detectable changes in their rotation periods. Using this sample and including the six comets used to determine $X$, we find a correlation between the parameter $\zeta$ (i.e. the total water production per unit surface area per orbit approximated by that inside of 4 au) and the perihelion distance $q$; specifically we derive $\zeta$ $\propto$ $q^{-0.8}$ and provide a theoretical basis for this in Appendix A. This relationship between $\zeta$ and $q$ enables ready comparisons of activity due to insolation between comets. Additionally, a relationship between the nuclear radius $R$ and the rotation period $P$ is found. Specifically, we find that on average smaller nuclei have smaller rotation periods compared to the rotation periods of larger nuclei. This is consistent with expectations for rotational evolution and spin-up of comet nuclei, providing strong observational evidence for sublimation-driven rotational changes in comets.

Sublimative outgassing of comets produces torques that alter the rotation state of their nuclei. Recently, parameterized sublimative torque models have been developed to study rotation state changes of individual comet nuclei and populations of cometary bodies. However, these models simplify the interactions between the escaping gas and cometary surface into only a few parameters that hide the details of these complex interactions. Here we directly compare the X-parameter model (Samarasinha & Mueller, 2013) with the SYORP model (Steckloff & Jacobson, 2016) to tease out insights into the details of the gas-surface interactions driving sublimative torques. We find that, for both of these models to accurately model sublimative torques, the number of sublimating molecules that contribute to the net torque is largely independent of the detailed shape and activity of the nucleus, but rather depends primarily on the size of the nucleus and the effective heliocentric distance of the comet. We suggest that cometary activity must be largely restricted to regions of steep gravitational surface slopes (above the angle of repose), where mass wasting can refresh activity by shedding mantles of refractory materials and exposing fresh volatiles. We propose a new classification scheme for comets based on the frequency of this mass-wasting process (relative to the timescale of activity fading): quasi-equilibrium, episodic, quasi-dormant, and extinct.

We imaged Comet 252P/2000 G1 (LINEAR) (hereafter 252P) with the Hubble Space Telescope and both 252P and P/2016 BA$_{14}$ (PanSTARRS) (hereafter BA$_{14}$) with the Discovery Channel Telescope in March and April 2016, surrounding its close encounter to Earth. The r'-band $Af\rho$ of 252P in a 0.2"-radius aperture were $16.8\pm0.3$ and $57\pm1$ cm on March 14 and April 4, respectively, and its gas production rates were: $Q$(OH) = $(5.8\pm0.1)\times10^{27}$ s$^{-1}$, and $Q$(CN) = $(1.25\pm0.01)\times10^{25}$ s$^{-1}$ on April 17. The r'-band upper limit $Af\rho$ of BA1$_{14}$ was $0.19\pm0.01$ cm in a 19.2"-radius aperture, and $Q$(CN) = $(1.4\pm0.1)10^{22}$ s$^{-1}$ on April 17, 2017. 252P shows a bright and narrow jet of a few hundred kilometers long in the sunward direction, changing its projected position angle in the sky with a periodicity consistent with 7.24 hours. However, its photometric lightcurve is consistent with a periodicity of 5.41 hours. We suggest that the nucleus of 252P is likely in a non-principal axis rotation. The nucleus radius of 252P is estimated to be about $0.3\pm0.03$ km, indicating an active fraction of 40% to >100% in its 2016 apparition. Evidence implies a possible cloud of slow-moving grains surrounding the nucleus. The activity level of 252P in the 2016 apparition increased by two orders of magnitude from its previous apparitions, making this apparition unusual. On the other hand, the activity level of BA14 appears to be at least three orders of magnitude lower than that of 252P, despite its ten times or larger surface area.

We observed Comet C/Siding Spring using the Hubble Space Telescope (HST) during its close approach to Mars. The high spatial resolution images obtained through the F689M, F775W, and F845M filters reveal the characteristics of the dust coma. The dust production rate of C/Siding Spring, quantified by $Af\rho$, is 590$\pm$30, 640$\pm$30, and 670$\pm$30 cm in a 420 km-radius aperture at 38$^\circ$ solar phase angle through the three filters, respectively, consistent with other observations at similar time and geometry, and with model predictions based on earlier measurements. The dust expansion velocity is ~150-250 m s$^{-1}$ for micron-sized dust grains, similar to the speeds found for other comets. The coma has a color slope of (5.5$\pm$1.5)%/100 nm between 689 and 845 nm, similar to previous HST measurements at comparable aperture sizes, consistent with the lack of color dependence on heliocentric distance for almost all previously observed active comets. The rotational period of the nucleus of C/Siding Spring is determined from the periodic brightness variation in the coma to be 8.00$\pm$0.08 hours, with no excited rotational state detected. The dust coma shows a broad and diffuse fan-shaped feature in the sunward direction, with no temporal morphological variation observed. The projected orientation of the dust feature, combined with the previous analysis of the coma morphology and other characteristics, suggests secular activity evolution of the comet in its inner solar system passage as one previously observed active region turns off whereas new regions exposed to sunlight due to seasonal illumination change.

We present the results of a global coma morphology campaign for comet C/2012 S1 (ISON), which was organized to involve both professional and amateur observers. In response to the campaign, many hundreds of images, from nearly two dozen groups were collected. Images were taken primarily in the continuum, which help to characterize the behavior of dust in the coma of comet ISON. The campaign received images from January 12 through November 22, 2013 (an interval over which the heliocentric distance decreased from 5.1 AU to 0.35 AU), allowing monitoring of the long-term evolution of coma morphology during the pre-perihelion leg of comet ISON. Data were contributed by observers spread around the world, resulting in particularly good temporal coverage during November when comet ISON was brightest but its visibility was limited from any one location due to the small solar elongation. We analyze the northwestern sunward continuum coma feature observed in comet ISON during the first half of 2013, finding that it was likely present from at least February through May and did not show variations on diurnal time scales. From these images we constrain the grain velocities to ~10 m/s, and we find that the grains spent 2-4 weeks in the sunward side prior to merging with the dust tail. We present a rationale for the lack of continuum coma features from September until mid-November 2013, determining that if the feature from the first half of 2013 was present, it was likely too small to be clearly detected. We also analyze the continuum coma morphology observed subsequent to the November 12 outburst, and constrain the first appearance of new features in the continuum to later than November 13.99 UT.

We present a new analysis of ${\it Spitzer}$ observations of comet 29P/Schwassmann-Wachmann 1 taken on UT 2003 November 21, 23, and 24, similar to a previous investigation of the observations (Stansberry et al. 2004), but using the most recent ${\it Spitzer}$ data pipeline products and intensive image processing techniques. Analysis of images from the IRAC 5.8 & 8.0 $\mu$m bands and the MIPS 24.0 & 70.0 $\mu$m bands resulted in photometry measurements of the nucleus after a suite of coma modeling and removal processes were implemented. SW1 was not identified in the 5.8 $\mu$m image from the previous work so its incorporation into this analysis is entirely new. Using the Near Earth Asteroid Thermal Model (Harris 1998) resulted in a nucleus radius measurement of $R$ = 30.2 $^{+3.7}_{-2.9}$ km and an infrared beaming parameter value of $\eta = 0.99$ $^{+0.26}_{-0.19}$. We also measured an infrared geometric albedo, $p_{5.8}$ = 0.5 $\pm$ 0.5. Extrapolating a 0.04 V-band albedo and using a normalized reflectivity gradient $S' = 14.94 \pm 1.09$ [% (1000 \r{A})$^{-1}$] (Duffard et al. 2014) we recover an infrared albedo of $p_{5.8}$ = 0.31 in the near infrared consistent with the value recovered from thermal modeling. The dust composition extracted from IRS spectra are very comet-like, containing mainly amorphous ferromagnesian silicates (but with a minority of crystalline silicates as well), water ice, and metal sulfides.

We present an analysis of Kitt Peak National Observatory and Lowell Observatory observations of comet 103P/Hartley 2 obtained from August through December 2010. The results are then compared with contemporaneous observations made by the EPOXI spacecraft. Each ground-based dataset has previously been investigated individually; the combined dataset has complementary coverage that reduces the time between observing runs and allows us to determine additional apparent periods at intermediate times. We compare CN coma morphology between ground-based datasets, making nine new measurements of apparent periods. The first five are consistent with the roughly linearly increasing apparent period during the apparition found by previous authors. The final four suggest that the change in apparent period slowed or stopped by late November. We also measure an inner coma lightcurve in both CN and R-band ground-based images, finding a single-peaked lightcurve which repeats in phase with the coma morphology. The apparent period from the lightcurve had significantly larger uncertainties than from the coma morphology, but varied over the apparition in a similar manner. Our ground-based lightcurve aligns with the published EPOXI lightcurve, indicating that the lightcurve represents changing activity rather than viewing geometry of structures in the coma. The EPOXI lightcurve can best be phased by a triple-peaked period near 54-55 hr that increases from October to November. This phasing reveals that the spacing between maxima is not constant, and that the overall lightcurve shape evolves from one triple-peaked cycle to the next. These behaviors suggest that much of the scatter in apparent periods derived from ground-based datasets acquired at similar epochs are likely due to limited sampling of the data.

Lightcurve observations of asteroids and bare cometary nuclei are the most widely used observational tool to derive the rotational parameters. Therefore, an in-depth understanding of how component periods of dynamically excited non-principal axis (NPA) rotators manifest in lightcurves is a crucial step in this process. We investigated this with the help of numerically generated lightcurves of NPA rotators with component periods known a priori. The component periods of NPA rotation were defined in terms of two widely used yet complementary conventions. We derive the relationships correlating the component rotation periods in the two conventions. These relationships were then used to interpret the periodicity signatures present in the simulated lightcurves and rationalize them in either convention.

The close encounter of Comet C/2013 A1 (Siding Spring) with Mars on October 19, 2014 presented an extremely rare opportunity to obtain the first flyby quality data of the nucleus and inner coma of a dynamically new comet. However, the comet's dust tail potentially posed an impact hazard to those spacecraft. To characterize the comet at large heliocentric distances, study its long-term evolution, and provide critical inputs to hazard modeling, we imaged C/Siding Spring with the Hubble Space Telescope when the comet was at 4.58, 3.77, and 3.28 AU from the Sun. The dust production rate, parameterized by the quantity Af$\rho$, was 2500, 2100, and 1700 cm (5000-km radius aperture) for the three epochs, respectively. The color of the dust coma is 5.0$\pm$0.3$\%$/100 nm for the first two epochs, and 9.0$\pm$0.3$\%$/100 nm for the last epoch, and reddens with increasing cometocentric distance out to ~3000 km from the nucleus. The spatial distribution and the temporal evolution of the dust color are most consistent with the existence of icy grains in the coma. Two jet-like dust features appear in the north-northwest and southeast directions projected in the sky plane. Within each epoch of 1-2 hour duration, no temporal variations were observed for either feature, but the PA of the southeastern feature varied between the three epochs by ~30$^\circ$. The dust feature morphology suggests two possible orientations for the rotational pole of the nucleus, (RA, Dec) = (295$^\circ\pm$5$^\circ$, +43$^\circ\pm$2$^\circ$) and (190$^\circ\pm$10$^\circ$, 50$^\circ\pm$5$^\circ$), or their diametrically opposite orientations.

Many cometary coma features are only a few percent above the ambient coma (i.e., the background) and therefore coma enhancement techniques are needed to discern the morphological structures present in cometary comae. A range of image enhancement techniques widely used by cometary scientists is discussed by categorizing them and carrying out a comparative analysis. The enhancement techniques and the corresponding characteristics are described in detail and the respective mathematical representations are provided. As the comparative analyses presented in this paper make use of simulated images with known coma features, the feature identifications as well as the artifacts caused by enhancement provide an objective and definitive assessment of the various techniques. Examples are provided which highlight contrasting capabilities of different techniques to pick out qualitatively distinct features of widely different strengths and spatial scales. On account of this as well as serious image artifacts and spurious features associated with certain enhancement techniques, confirmation of the presence of coma features using qualitatively different techniques is strongly recommended.

Comet C/2013 A1 (Siding Spring) will have a close encounter with Mars on October 19, 2014. We model the dynamical evolution of dust grains from the time of their ejection from the comet nucleus to the Mars close encounter, and determine the flux at Mars. Constraints on the ejection velocity from Hubble Space Telescope observations indicate that the bulk of the grains will likely miss Mars, although it is possible that a few-percent of grains with higher velocities will reach Mars, peaking approximately 90--100 minutes after the close approach of the nucleus, and consisting mostly of millimeter-radius grains ejected from the comet nucleus at a heliocentric distance of approximately 9~AU or larger. At higher velocities, younger grains from sub-millimeter to several millimeter can reach Mars too, although an even smaller fraction of grains is expected have these velocities, with negligible effect on the peak timing. Using NEOWISE observations of the comet, we can estimate that the maximum fluence will be of the order of $10^{-7}$ grains/m$^2$. We include a detailed analysis of how the expected fluence depends on the grain density, ejection velocity, and size-frequency distribution, to account for current model uncertainties and in preparation of possible refined model values in the near future.

We introduce a parameter, X, to predict the changes in the rotational period of a comet in terms of the rotational period itself, the nuclear radius, and the orbital characteristics. We show that X should be a constant if the bulk densities and shapes of nuclei are nearly identical and the activity patterns are similar for all comets. For four nuclei for which rotational changes are well documented, despite the nearly factor 30 variation observed among the effective active fractions of these comets, X is constant to within a factor two. We present an analysis for the sungrazing comet C/2012 S1 (ISON) to explore what rotational changes it could undergo during the upcoming perihelion passage where its perihelion distance will be ~2.7 solar radii. When close to the sun, barring a catastrophic disruption of the nucleus, the activity of ISON will be sufficiently strong to put the nucleus into a non-principal-axis rotational state and observable changes to the rotational period should also occur. Additional causes for rotational state changes near perihelion for ISON are tidal torques caused by the sun and the significant mass loss due to a number of mechanisms resulting in alterations to the moments of inertia of the nucleus.

The CN coma structure of the EPOXI mission target, comet 103P/Hartley 2, was observed during twenty nights from September to December 2010. These CN images probe the rotational state of the comet's nucleus and provide a ground-based observational context to complement the EPOXI observations. A dynamically excited cometary nucleus with a changing rotational rate is observed, a characteristic not seen in any comet in the past. The lack of rotational damping during the four-month observing interval places constraints on the interior structure of the nucleus.