Near-nucleus photometry of comets using archived NEAT data

Though optimized to discover and track fast moving Near-Earth Objects (NEOs), the Near-Earth Asteroid Tracking (NEAT) survey dataset can be mined to obtain information on the comet population observed serendipitously during the asteroid survey. We have completed analysis of over 400 CCD images of comets obtained during the autonomous operations of two 1.2-m telescopes: the first on the summit of Haleakala on the Hawaiian island of Maui and the second on Palomar Mountain in southern California. Photometric calibrations of each frame were derived using background catalog stars and the near-nucleus comet photometry measured. We measured dust production and normalized magnitudes for the coma and nucleus in order to explore cometary activity and comet size-frequency distributions. Our data over an approximately two-year time frame (2001 August-2003 February) include 52 comets: 12 periodic, 19 numbered, and 21 non-periodic, obtained over a wide range of viewing geometries and helio/geocentric distances. Nuclear magnitudes were estimated for a subset of comets observed. We found that for low-activity comets (Afρ<100 cm) our model gave reasonable estimates for nuclear size and magnitude. The slope of the cumulative luminosity function of our sample of low-activity comets was 0.33 ± 0.04, consistent with the slope we measured for the Jupiter-family cometary nuclei collected by Fernández et al. [Fernández, J.A., Tancredi, G., Rickman, H., Licandro, J., 1999. Astron. Astrophys. 392, 327-340] of 0.38 ± 0.02. Our slopes of the cumulative size distribution α=1.50±0.08 agree well with the slopes measured by Whitman et al. [Whitman, K., Morbidelli, A., Jedicke, R., 2006. Icarus 183, 101-114], Meech et al. [Meech, K.J., Hainaut, O.R., Marsden, B.G., 2004. Icarus 170, 463-491], Lowry et al. [Lowry, S.C., Fitzsimmons, A., Collander-Brown, S., 2003. Astron. Astrophys. 397, 329-343], and Weissman and Lowry [Weissman, P.R., Lowry, S.C., 2003. Lunar Planet. Sci. 34. Abstract 34].     

The population of faint Jupiter family comets near the Earth

We study the population of faint Jupiter family comets (JFCs) that approach the Earth (perihelion distances q<1.3 AU) by applying a debiasing technique to the observed sample. We found for the debiased cumulative luminosity function (CLF) of absolute total magnitudes H₁₀ a bimodal distribution in which brighter comets (H₁₀?9) follow a linear relation with a steep slope α=0.65±0.14, while fainter comets follow a much shallower slope α=0.25±0.06 down to H₁₀∼18. The slope can be pushed up to α=0.35±0.09 if a second break in the H₁₀ distribution to a much shallower slope is introduced at H₁₀∼16. We estimate a population of about 10³ faint JFCs with q<1.3 AU and 10<H₁₀<15 (radii ∼0.1-0.5 km). The shallowness of the CLF for faint near-Earth JFCs may be explained either as: (i) the source population (the scattered disk) has an equally very shallow distribution in the considered size range, or (ii) the distribution is flattened by the disintegration of small objects before that they have a chance of being observed. The fact that the slope of the magnitude distribution of the faint active JFCs is very similar to that found for a sample of dormant JFCs candidates suggests that for a surviving (i.e., not disintegrated) object, the probability of becoming dormant versus keeping some activity is roughly size independent.     

Are There Many Inactive Jupiter-Family Comets among the Near-Earth Asteroid Population?

We analyze the dynamical evolution of Jupiter-family (JF) comets and near-Earth asteroids (NEAs) with aphelion distances Q>3.5 AU, paying special attention to the problem of mixing of both populations, such that inactive comets may be disguised as NEAs. From numerical integrations for 2×10⁶ years we find that the half lifetime (where the lifetime is defined against hyperbolic ejection or collision with the Sun or the planets) of near-Earth JF comets (perihelion distances q<1.3 AU) is about 1.5×10⁵ years but that they spend only a small fraction of this time (∼ a few 10³ years) with q<1.3 AU. From numerical integrations for 5×10⁶ years we find that the half lifetime of NEAs in “cometary” orbits (defined as those with aphelion distances Q>4.5 AU, i.e., that approach or cross Jupiter's orbit) is 4.2×10⁵ years, i.e., about three times longer than that for near-Earth JF comets. We also analyze the problem of decoupling JF comets from Jupiter to produce Encke-type comets. To this end we simulate the dynamical evolution of the sample of observed JF comets with the inclusion of nongravitational forces. While decoupling occurs very seldom when a purely gravitational motion is considered, the action of nongravitational forces (as strong as or greater than those acting on Encke) can produce a few Enckes. Furthermore, a few JF comets are transferred to low-eccentricity orbits entirely within the main asteroid belt (Q<4 AU and q>2 AU). The population of NEAs in cometary orbits is found to be adequately replenished with NEAs of smaller Q's diffusing outward, from which we can set an upper limit of ∼20% for the putative component of deactivated JF comets needed to maintain such a population in steady state. From this analysis, the upper limit for the average time that a JF comet in near-Earth orbit can spend as a dormant, asteroid-looking body can be estimated to be about 40% of the time spent as an active comet. More likely, JF comets in near-Earth orbits will disintegrate once (or shortly after) they end their active phases.     

The size-frequency distribution of dormant Jupiter family comets

We estimate the total number and the slope of the size-frequency distribution (SFD) of dormant Jupiter family comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 near-Earth objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J., Levison, H.F., Michel, P., Metcalfe, T.S., 2002a. Icarus 156, 399-433]. Then, from the numerical simulations of Levison and Duncan [1997. Icarus 127, 13-32], we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit-independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J., Bottke, W.F., 2003. Icarus 161, 17-33]. By fitting the biased model to the data, we estimate that there are ∼75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is −1.5±0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry [2003. Lunar Planet. Sci. 34. Abstract 2003]. Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size-independent. Our results imply that the size distribution of the JFC progenitors—the scattered disk trans-neptunian population—either (i) has a similar and shallow SFD or (i^′) is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size-dependent impact strength [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20], suggests that scattered disk bodies reached collisional equilibrium inside the protoplanetary disk prior to their removal from the planetary region.     

Comets in the near-Earth object population

Because the lifespan of near-Earth objects (NEOs) is shorter than the age of the Solar System, these objects originate elsewhere. Their most likely sources are the main asteroid belt and comets. Through physical observations we seek to identify potential dormant or extinct comets among “asteroids” catalogued as NEOs and thereby determine the fraction of “comet candidates” within the total NEO population. Both discovery statistics and dynamical models indicate that candidate cometary objects in near-Earth space are predominantly found among those having a jovian Tisserand parameter Tj<3. Therefore, we seek to identify comet candidates among asteroid-like NEOs using three criteria: Tj<3, spectral parameters (C, D, T, or P taxonomic types), and/or low (<0.075) albedos. We present new observations for 20 NEOs having Tj<3, consisting of visible spectra, near-infrared spectra, and/or albedo measurements obtained using the NASA Infrared Telescope Facility, the Kitt Peak National Observatory 4 m, and the Magellan Observatory 6.5-m. Four of our “asteroid” targets have been subsequently confirmed as low activity comets. Thus our sample includes spectra of the nuclei of Comets 2002 EX12 = 169P (NEAT), 2001 WF2 = 182P (LONEOS), 2003 WY25 = D/1891 W1 (Blanplain), and Halley Family Comet 2006 HR30 = P/2006 HR30 (Siding Spring). From the available literature, we tabulate physical properties for 55 NEOs having Tj<3, and after accounting for possible bias effects, we estimate that 54±10% of NEOs in Tj<3 orbits have “comet-like” spectra or albedos. Bias corrected discovery statistics [Stuart, J.S., Binzel, R.P., 2004. Icarus 170, 295-311] estimate 30±5% of the entire NEO population resides in orbits having Tj<3. Combining these two factors suggests that 16±5% of the total discovered “asteroid-like” NEO population has “comet-like” dynamical and physical properties. Outer main-belt asteroids typically have similar taxonomic and albedo properties as our “comet candidates.” Using the model of Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.M., Levison, H., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] to evaluate source region probabilities, we conclude that 8±5% of the total asteroid-like NEO population have the requisite orbital properties, physical properties, and dynamical likelihood to have originated as comets from the outer Solar System.     

The size distribution of trans-Neptunian bodies

The magnitude distribution of the trans-Neptunian bodies composed of the Kuiper Belt Objects (KBOs) and Scattered Disk Objects (SDOs) is determined for absolute magnitudes H?7, using maximum likelihood estimation methods. This is translated into a corresponding size distribution. This gave a differential size index of q=3.966±0.15 for KBOs and q=3.016±0.32 for SDOs. It was found that these two distributions were statistically different. The KBOs were further split into classical KBOs and Plutinos which had indices of q=4.074±0.18 and q=3.301±0.37, respectively. There was no statistical evidence that these are different populations. The classical KBOs were further split and examined for four different semi-major axis ranges and it was found that there was moderate evidence that the entire sample was not well represented by one index. The distribution indices of the SDOs were compared with the distributions of short period comets and found to be similar. It is likely that the scattered disk population is the source of the short period comets.     

Absolute magnitudes of asteroids and a revision of asteroid albedo estimates from WISE thermal observations

We obtained estimates of the Johnson V absolute magnitudes (H) and slope parameters (G) for 583 main-belt and near-Earth asteroids observed at Ond?ejov and Table Mountain Observatory from 1978 to 2011. Uncertainties of the absolute magnitudes in our sample are <0.21 mag, with a median value of 0.10 mag. We compared the H data with absolute magnitude values given in the MPCORB, Pisa AstDyS and JPL Horizons orbit catalogs. We found that while the catalog absolute magnitudes for large asteroids are relatively good on average, showing only little biases smaller than 0.1 mag, there is a systematic offset of the catalog values for smaller asteroids that becomes prominent in a range of H greater than ～10 and is particularly big above H ～ 12. The mean (H_catalog ? H) value is negative, i.e., the catalog H values are systematically too bright. This systematic negative offset of the catalog values reaches a maximum around H = 14 where the mean (H_catalog ? H) is ?0.4 to ?0.5. We found also smaller correlations of the offset of the catalog H values with taxonomic types and with lightcurve amplitude, up to ～0.1 mag or less. We discuss a few possible observational causes for the observed correlations, but the reason for the large bias of the catalog absolute magnitudes peaking around H = 14 is unknown; we suspect that the problem lies in the magnitude estimates reported by asteroid surveys. With our photometric H and G data, we revised the preliminary WISE albedo estimates made by Masiero et al. (Masired, J.R. et al. [2011]. Astrophys. J. 741, 68–89) and Mainzer et al. (Mainzer, A. et al. [2011b]. Astrophys. J. 743, 156–172) for asteroids in our sample. We found that the mean geometric albedo of Tholen/Bus/DeMeo C/G/B/F/P/D types with sizes of 25–300 km is p_V = 0.057 with the standard deviation (dispersion) of the sample of 0.013 and the mean albedo of S/A/L types with sizes 0.6–200 km is 0.197 with the standard deviation of the sample of 0.051. The standard errors of the mean albedos are 0.002 and 0.006, respectively; systematic observational or modeling errors can predominate over the quoted formal errors. There is apparent only a small, marginally significant difference of 0.031 ± 0.011 between the mean albedos of sub-samples of large and small (divided at diameter 25 km) S/A/L asteroids, with the smaller ones having a higher albedo. The difference will have to be confirmed and explained; we speculate that it may be either a real size dependence of surface properties of S type asteroids or a small size-dependent bias in the data (e.g., a bias towards higher albedos in the optically-selected sample of asteroids). A trend of the mean of the preliminary WISE albedo estimates increasing with asteroid size decreasing from D ～ 30 down to ～5 km (for S types) showed in Mainzer et al. (Mainzer, A. et al. [2011a]. Astrophys. J. 741, 90–114) appears to be mainly due to the systematic bias in the MPCORB absolute magnitudes that progressively increases with H in the corresponding range H = 10–14.     

Cometary nuclear magnitudes from sky survey observations

The determination of the nuclear magnitudes of comets, and with it nuclear size frequency distributions, is strongly complicated by cometary activity. By now, only nuclear size frequency distributions for Jupiter Family comets are available, and they are still subject of uncertainties. For comets of other dynamical classes, nuclear magnitudes are known for only a few comets. The size frequency distributions are thus not well constrained.In this work we study whether nuclear magnitudes of comets can be constrained from sky survey observations as published by the Minor Planet Center. Observations from sky survey programs in which the comet was classified as a point-like source are analyzed in this respect.From the available published observations from 1998 to 2008, we derive nuclear magnitudes, as well as nuclear radii, for 84 comets. Among these are comets of the Jupiter Family, dynamically old and new isotropic comets, Halley-type comets and Centaurs. For Jupiter Family comets and for isotropic comets, the size frequency distributions are presented.Uncertainties of derived nuclear magnitudes arise from photometry and from potentially undetected activity. However, a comparison with objects with well known nuclear parameters shows that, despite substantial observational uncertainties, nuclear magnitudes are constrained to ±0.6 mag, thereby providing first indications for nuclear sizes. This is particularly relevant for isotropic comets with so far ill-constrained size distributions. Exponents of the differential size frequency distributions of for Jupiter Family comets and for isotropic comets are presented. The values derived here form a basis for future, dedicated observational studies which provide higher measurement accuracy.     

Comet nucleus size distributions from HST and Keck telescopes

The Wide Field Camera (WFC) on the Hubble Space Telescope and the Low Resolution Imaging Spectrograph (LRIS) on the Keck II telescope have been used to image 21 distant dynamically new, long-period (LP) and short-period (SP) Jupiter-family (JF) comet nuclei (near aphelion), as part of a long-term program to search for physical differences between short-period comets and Oort cloud comets. WFC data were obtained on Comets C/1987 H1 (Shoemaker) and C/1984 K1 (Shoemaker) during Cycle 5 (1995 December) and on C/1988 B1 (Shoemaker), C/1987 F1 (Torres), and C/1983 O1 (?ernis) during Cycle 6 (1997 April, May, and June). The HST comets were at heliocentric distances 20.4 < r[AU] < 29.5. Each comet observation was allocated 7 orbits, for ≈3.6 hrs of integration. The most difficult part of the image reduction was the removal of cosmic rays. We present our scheme for cosmic ray removal. None of the HST comet nuclei was detected to the 3-σ level at m_R∼27. The inferred upper limits to the nucleus radii are . The SP comets range in radius between , with a median value of R_N∼1.61 km. The LP comets ranged in size between <4.0-56 km. Over a range of radii between 1-10 km, the nuclei can be fit with a cumulative distribution N(>R_N)∝R_N^−α with α=1.45±0.05, and for nuclei in the range 2-5 km, α=1.91±0.06. Statistical analysis and modeling shows that the slopes of the observed TNO and JF comet distributions are not compatible, suggesting that the intrinsic distribution of JF comet nuclei is a differential a^−3.5 power law truncated at small nucleus radii between 0.3 and 2.0 km.     

A Model for the Common Origin of Jupiter Family and Halley Type Comets

A numerical simulation of the Oort cloud is used to explain the observed orbital distributions and numbers of Jupiter-family (JF) and Halley-type (HT) short-period (SP) comets. Comets are given initial orbits with perihelion distances between 5 and 36 au, and evolve under planetary, stellar and Galactic perturbations for 4.5 Gyr. This process leads to the formation of an Oort cloud (which we define as the region of semimajor axes a > 1,000 au), and to a flux of cometary bodies from the Oort cloud returning to the planetary region at the present epoch. The results are consistent with the dynamical characteristics of SP comets and other observed cometary populations: the near-parabolic flux, Centaurs, and high-eccentricity trans-Neptunian objects. To achieve this consistency with observations, the model requires that the number of comets versus initial perihelion distance is concentrated towards the outer planetary region. Moreover, the mean physical lifetime of observable comets in the inner planetary region (q < 2.5 au) at the present epoch should be an increasing function of the comets’ initial perihelion distances. Virtually all observed HT comets and nearly half of observed JF comets come from the Oort cloud, and initially (4.5 Gyr ago) from orbits concentrated near the outer planetary region. Comets that have been in the Oort cloud also return to the Centaur (5 < q < 28 au, a < 1,000 au) and near-Neptune high-eccentricity regions. Such objects with perihelia near Neptune are hard to discover, but Centaurs with characteristics predicted by the model (e.g. large semimajor axes, above 60 au, or high inclinations, above 40°) are increasingly being found by observers. The model provides a unified picture for the origin of JF and HT comets. It predicts that the mean physical lifetime of all comets in the region q < 1.5 au is less than ～200 revolutions.     

Modeling of cometary evolution by kinetic theory: Method and first results

Physical evolution of Jupiter family (JF) comets is considered as a simultaneous process of erosion and fading. Dynamical effects are limited to discrete changes of the perihelion distance, that result in changes of the evaporation rate. Assuming that the JF comet population is in a steady state, a distribution function of this population in the two dimensional phase space consisting of radius and active fraction of the nucleus surface is found as the solution of a set of kinetic equations, each one of them for a different perihelion distance. With use of the distribution function some statistical properties of the comet population, like the total number of comets in the considered region of the phase space, the number of objects that evaporate or get dormant per unit time, etc., are obtained. The cumulative distribution function with respect to the absolute brightness is calculated and compared with the observed one as a check on the considered models.     

Modeling of cometary evolution by kinetic theory: Method and first results

Physical evolution of Jupiter family (JF) comets is considered as a simultaneous process of erosion and fading. Dynamical effects are limited to discrete changes of the perihelion distance, that result in changes of the evaporation rate. Assuming that the JF comet population is in a steady state, a distribution function of this population in the two dimensional phase space consisting of radius and active fraction of the nucleus surface is found as the solution of a set of kinetic equations, each one of them for a different perihelion distance. With use of the distribution function some statistical properties of the comet population, like the total number of comets in the considered region of the phase space, the number of objects that evaporate or get dormant per unit time, etc., are obtained. The cumulative distribution function with respect to the absolute brightness is calculated and compared with the observed one as a check on the considered models.     

Nuclear spin temperature of ammonia in Comet 9P/Tempel 1 before and after the Deep Impact event

The Deep Impact mission succeeded in excavating inner materials from the nucleus of Comet 9P/Tempel 1 on 2005 July 04 (at 05:52 UT). Comet 9P/Tempel 1 is one of Jupiter family short period comets, which might originate in the Kuiper belt region in the solar nebula. In order to characterize the comet and to support the mission from the ground-based observatory, optical high-dispersion spectroscopic observations were carried out with the echelle spectrograph (UVES) mounted on the 8-m telescope VLT (UT2) before and after the Deep Impact event. Ortho-to-para abundance ratios (OPRs) of cometary ammonia were determined from the NH₂ emission spectra. The OPRs of ammonia on July 3.996 UT and 4.997 UT were derived to be 1.28±0.07 (nuclear spin temperature: Tspin=24±2 K) and 1.26±0.08 (Tspin=25±2 K), respectively. There is no significant change between before and after the impact. Actually, most materials ejected from the impact site could have moved away from the nucleus on July 4.997 UT, about 17 h after the impact. However, a small fraction of the ejected materials might remain in the slit of UVES instrument at that time because an excess of about 20% in the NH₂ emission flux is observed above the normal activity level was found [Manfroid, J., Hutsemékers, D., Jehin, E., Cochran, A.L., Arpigny, C., Jackson, W.M., Meech, K.J., Schulz, R., Zucconi, J.-M., 2007. Icarus. This issue]. If the excess of NH₂ on July 04.997 UT was produced from icy materials excavated by the Deep Impact, then an upper-limit of the ammonia OPR would be 1.75 (Tspin>17 K) for those materials. On the other hand, the OPR of ammonia produced from the quiescent sources was similar to that of the Oort cloud comets observed so far. This fact may imply that physical conditions where cometary ices formed were similar between Comet 9P/Tempel 1 and the Oort cloud comets.     

Keck NIRC2 photometry of Uranus, uranian satellites, and Triton in August 2004

On August 11, 2004, we made adaptive optics observations of the Uranus and Neptune systems with the Keck II Near Infrared Camera. Uranus and Triton were observed in three broadband filters (J, H, and K-prime) and four narrowband filters [Hcont, FeII, He1_B, and H2(v=1-0)]. Miranda, Ariel, Umbriel, and Oberon were observed in the four narrowband filters only. To achieve the highest possible photometric accuracy, and thus the tightest possible constraints on atmospheric aerosol models and gas mixing ratios, we used aperture photometry that accounted for the dependence of point-spread functions and growth curves on the adaptive optics reference object, and accounted for recently determined phase curves of each object. The satellite albedos we determined were compared with published relative spectra to verify the relative consistency of our observations, which were subsequently used to convert published relative spectra to absolute spectra. We used these absolute spectra and synthetic photometry methods to compare published values in dissimilar photometric systems to each other and to our observations. We find our satellite albedos in best agreement with photometry from Karkoschka [Karkoschka, E., 2001. Icarus 151, 51-68], which is the best characterized and most contemporaneous data set. Our estimated absolute accuracy of 5% to 8% is consistent with these intercomparisons. We obtained the following ring-subtracted and discrete feature-free albedos of Uranus: J: (1.66±0.07)×¹⁰⁻², H: (1.09±0.05)×¹⁰⁻², K^′: (2.08±0.15)×¹⁰⁻⁴, Hcont: (3.71±0.23)×¹⁰⁻², FeII: (1.14±0.07)×¹⁰⁻³, He1_B: (2.06±0.16)×¹⁰⁻⁴, and H2: (1.27±0.10)×¹⁰⁻⁴.     

Deep Impact photometry of Comet 9P/Tempel 1

The photometric properties of the nucleus of Comet 9P/Tempel 1 are studied from the disk-resolved color images obtained by Deep Impact (DI). Comet Tempel 1 has typical photometric properties for comets and dark asteroids. The disk-integrated spectrum of the nucleus of Tempel 1 between 309 and 950 nm is linear without any features at the spectral resolution of the filtered images. At V-band, the red slope of the nucleus is 12.5±1% per 100 nm at 63° phase angle, translating to B-V=0.84±0.01, V-R=0.50±0.01, and R-I=0.49±0.02. No phase reddening is confirmed. The phase function of the nucleus of Tempel 1 is constructed from DI images and earlier ground-based observations found from the literature. The phase coefficient is determined to be β=0.046±0.007 mag/deg between 4° and 117° phase angle. Hapke's theoretical scattering model was used to model the photometric properties of this comet. Assuming a single Henyey-Greenstein function for the single-particle phase function, the asymmetry factor of Tempel 1 was fitted to be g=−0.49±0.02, and the corresponding single-scattering albedo (SSA) was modeled to be 0.039±0.005 at 550 nm wavelength. The SSA spectrum shows a similar linear slope to that of the disk-integrated spectrum. The roughness parameter is found to be 16°±8°, and independent of wavelength. The Minnaert k parameter is modeled to be 0.680±0.014. The photometric variations on Tempel 1 are relatively small compared to other comets and asteroids, with a ∼20% full width at half maximum of albedo variation histogram, and ∼3% for color. Roughness variations are evident in one small area, with a roughness parameter about twice the average and appearing to correlate with the complex morphological texture seen in high-resolution images.     

Observational Constraints On Surface Characteristics Of Comet Nuclei

Direct observations of the nuclear surfaces of comets have been difficult; however a growing number of studies are overcoming observational challenges and yielding new information on cometary surfaces. In this review, we focus on recent determinations of the albedos, reflectances, and thermal inertias of comet nuclei. There is not much diversity in the geometric albedo of the comet nuclei observed so far (a range of 0.025 to 0.06). There is a greater diversity of albedos among the Centaurs, and the sample of properly observed TNOs(2) is still too small. Based on their albedos and Tisser and invariants, Fernández et al. (2001) estimate that about 5% of the near-Earth asteroids have a cometary origin, and place an upper limit of 10%. The agreement between this estimate and two other independent methods provide the strongest constraint to date on the fraction of objects that comets contribute to the population of near-Earth asteroids. There is a diversity of visible colors among comets, extinct comet candidates, Centaurs and TNOs. Comet nuclei are clearly not as red as the reddest Centaurs and TNOs. What Jewitt (2002) calls ultra-red matter seems to be absent from the surfaces of comet nuclei. Rotationally resolved observations of both colors and albedos are needed to disentangle the effects of rotational variability from other intrinsic qualities. New constraints on thermal inertia of comets are consistent with previous independent estimates. The thermal inertia estimates for Centaurs 2060 Chiron and 8405 Asbolus are significantly lower than predicted by thermal models, and also lower than the few upper limits or constraints known for active, ordinary nuclei.     

Albedo and size determination of potentially hazardous asteroids: (99942) Apophis

The near-Earth object (99942) Apophis will make an extremely close approach to the Earth in 2029, and currently has approximately a one-in-45,000 chance of impacting our planet in 2036 (JPL Sentry, November 2006). Computation of the orbital evolution of this object is limited by insufficient knowledge of physical properties required to determine the role played by non-gravitational effects. Using polarimetric observations, we have obtained the first reliable determination of the albedo of Apophis, obtaining 0.33±0.08. We also derive an updated estimate of the asteroid's absolute magnitude: H=19.7±0.4. Using this albedo and H, we find that Apophis has a diameter of 270±60 m, slightly smaller than preliminary estimates based upon an assumed albedo. Our observations demonstrate the feasibility of polarimetric observations aimed at obtaining albedos and sizes of small, potentially hazardous asteroids.     

Disk-integrated bolometric Bond albedos and rotational light curves of saturnian satellites from Cassini Visual and Infrared Mapping Spectrometer

We present values from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) of four fundamental disk-integrated spectrophotometric properties (bolometric Bond albedo, solar phase curve, phase integral, and geometric albedo at 7-15 different wavelengths in the λ = 0.35-5.1 μm range) for five mid-sized saturnian icy satellites: Rhea, Dione, Tethys, Mimas, and Enceladus. These values, which include data from the period 2004-2008 and add to past VIMS phase curves, include opposition surge effects at down to fractions of a degree in solar phase angle for several moons and extend to over double the solar phase angle coverage of the Voyager mission. We also present new rotational light curves for Rhea and Dione at 7 near-infrared bands not previously available in ground-based or spacecraft studies. The bolometric Bond albedos we derive are as follows: 0.48 ± 0.09 (Rhea), 0.52 ± 0.08 (Dione), 0.61 ± 0.09 (Tethys), 0.67 ± 0.10 (Mimas), and 0.85 ± 0.11 (Enceladus). We also provide breakdowns of the major photometric quantities in both leading and trailing hemispheres. These refined parameters can be used to construct future bolometric Bond albedo maps that will contribute to surface composition identification studies, as well as models of volatile transport and sublimation. Through such applications, these data will help to determine the physical properties of surface particles, how the E-ring affects the inner saturnian moons, what is responsible for the dark albedo patterns seen on Tethys, and if these moons (e.g., Dione) are geologically active.     

Asteroid age distributions determined by space weathering and collisional evolution models

We provide evidence of consistency between the dynamical evolution of main belt asteroids and their color evolution due to space weathering. The dynamical age of an asteroid’s surface (Bottke, W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H. [2005]. Icarus 175 (1), 111-140; Nesvorný, D., Jedicke, R., Whiteley, R.J., Ivezi?, ?. [2005]. Icarus 173, 132-152) is the time since its last catastrophic disruption event which is a function of the object’s diameter. The age of an S-complex asteroid’s surface may also be determined from its color using a space weathering model (e.g. Willman, M., Jedicke, R., Moskovitz, N., Nesvorný, D., Vokrouhlický, D., Mothé-Diniz, T. [2010]. Icarus 208, 758-772; Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M. [2004]. Nature 429, 275-277; Willman, M., Jedicke, R., Nesvorny, D., Moskovitz, N., Ivezi?, ?., Fevig, R. [2008]. Icarus 195, 663-673. We used a sample of 95 S-complex asteroids from SMASS and obtained their absolute magnitudes and u, g, r, i, z filter magnitudes from SDSS. The absolute magnitudes yield a size-derived age distribution. The u, g, r, i, z filter magnitudes lead to the principal component color which yields a color-derived age distribution by inverting our color-age relationship, an enhanced version of the ‘dual τ’ space weathering model of Willman et al. (2010).We fit the size-age distribution to the enhanced dual τ model and found characteristic weathering and gardening times of τ_w = 2050 ± 80 Myr and respectively. The fit also suggests an initial principal component color of −0.05 ± 0.01 for fresh asteroid surface with a maximum possible change of the probable color due to weathering of ΔPC = 1.34 ± 0.04. Our predicted color of fresh asteroid surface matches the color of fresh ordinary chondritic surface of PC₁ = 0.17 ± 0.39.