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.     

Radar observations of Asteroids 64 Angelina and 69 Hesperia

We report new radar observations of E-class Asteroid 64 Angelina and M-class Asteroid 69 Hesperia obtained with the Arecibo Observatory S-band radar (2480 MHz, 12.6 cm). Our measurements of Angelina’s radar bandwidth are consistent with reported diameters and poles. We find Angelina’s circular polarization ratio to be 0.8 ± 0.1, tied with 434 Hungaria for the highest value observed for main-belt asteroids and consistent with the high values observed for all E-class asteroids (Benner, L.A.M., Ostro, S.J., Magri, C., Nolan, M.C., Howell, E.S., Giorgini, J.D., Jurgens, R.F., Margot, J.L., Taylor, P.A., Busch, M.W., Shepard, M.K. [2008]. Icarus 198, 294-304; Shepard, M.K., Kressler, K.M., Clark, B.E., Ockert-Bell, M.E., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J. [2008b]. Icarus 195, 220-225). Our radar observations of 69 Hesperia, combined with lightcurve-based shape models, lead to a diameter estimate, D_eff = 110 ± 15 km, approximately 20% smaller than the reported IRAS value. We estimate Hesperia to have a radar albedo of , consistent with a high-metal content. We therefore add 69 Hesperia to the Mm-class (high metal M) (Shepard, M.K., Clark, B.E., Ockert-Bell, M., Nolan, M.C., Howell, E.S., Magri, C., Giorgini, J.D., Benner, L.A.M., Ostro, S.J., Harris, A.W., Warner, B.D., Stephens, R.D., Mueller, M. [2010]. Icarus 208, 221-237), bringing the total number of Mm-class objects to eight; this is 40% of all M-class asteroids observed by radar to date.     

Looking a gift horse in the mouth: Evaluation of wide-field asteroid photometric surveys

It has recently become possible to do a photometric survey of many asteroids at once, rather than observing single asteroids one (or occasionally a couple) at a time. We evaluate two such surveys. Dermawan et al. (Dermawan et al. [2011]. Publ. Astron. Soc. Jpn. 63, S555–S576) observed one night on the Subaru 8.2 m telescope, and Masiero et al. (Masiero, J., Jedicke, R., Durech, J., Gwen, S., Denneau, L., Larsen, J. [2009]. Icarus 204, 145–171) observed six nights over 2 weeks with the 3.6 m CFHT. Dermawan claimed 83 rotation periods from 127 detected asteroids; Masiero et al. claimed 218 rotation periods from 828 detections. Both teams claim a number of super-fast rotators (P < 2.2 h) among main belt asteroids larger than 250 m diameter, some up to several km in diameter. This would imply that the spin rate distribution of main belt asteroids differs from like-sized NEAs, that there are larger super-fast rotators (monolithic asteroids) in the main belt than among NEAs. Here we evaluate these survey results, applying the same criteria for reliability of results that we apply to all results listed in our Lightcurve Database (Warner, B.D., Harris, A.W., Pravec, P. [2009a]. Icarus 202, 134–146). In doing so, we assigned reliability estimates judged sufficient for inclusion in statistical studies for only 27 out of 83 (33%) periods claimed by Dermawan, and only 87 out of 218 (40%) periods reported by Masiero et al.; none of the super-fast rotators larger than about 250 m diameter claimed by either survey received a reliability rating judged sufficient for analysis. We find no reliable basis for the claim of different rotation properties between main belt and near-Earth asteroids. Our analysis presents a cautionary message for future surveys.     

2007 Mutual events within the binary system of (22) Kalliope

In 2007, the asteroid Kalliope will reach one of its annual equinoxes. As a consequence, its small satellite Linus orbiting in the equatorial plane will undergo a season of mutual eclipses and occultations very similar to the one that the Galilean satellites undergo every 6 years. This paper is aimed at preparing a campaign of observations of these mutual events occurring from February to May 2007. This opportunity occurs only under favorable geometric conditions when the Sun and/or the Earth are close to the orbital plane of the system. This is the first international campaign devoted to the observation of photometric events within an asynchronous asteroidal binary system. We took advantage of a reliable orbit solution of Linus to predict a series of 24 mutual eclipses and 12 mutual occultations observable in the spring of 2007. Thanks to the brightness of Kalliope (m_v?11), these observations are easy to perform even with a small telescope. Anomalous attenuation events could be observed lasting for about 1-3 h with amplitude up to 0.09 mag. The attenuations are of two distinct types that can clearly be identified as primary and secondary eclipses similar to those that have been previously observed in other minor planet binary systems [Pravec, P., Scheirich, P., Kusnirák, P., Sarounová, L., Mottola, S., Hahn, G., Brown, P., Esquerdo, G., Kaiser, N., Krzeminski, Z., Pray, D.P., Warner, B.D., Harris, A.W., Nolan, M.C., Howell, E.S., Benner, L.A.M., Margot, J.-L., Galád, A., Holliday, W., Hicks, M.D., Krugly, Yu.N., Tholen, D., Whiteley, R., Marchis, F., Degraff, D.R., Grauer, A., Larson, S., Velichko, F.P., Cooney, W.R., Stephens, R., Zhu, J., Kirsch, K., Dyvig, R., Snyder, L., Reddy, V., Moore, S., Gajdos, S., Világi, J., Masi, G., Higgins, D., Funkhouser, G., Knight, B., Slivan, S., Behrend, R., Grenon, M., Burki, G., Roy, R., Demeautis, C., Matter, D., Waelchli, N., Revaz, Y., Klotz, A., Rieugné, M., Thierry, P., Cotrez, V., Brunetto, L., Kober, G., 2006. Photometric survey of binary near-Earth asteroids. Icarus 181, 63-93]. With these favorable circumstances, such photometric observations will provide us tight constraints regarding physical properties of Linus such as the size, shape and synodic spin period.     

Physical investigation of the potentially hazardous Asteroid (144898) 2004 VD17

In this paper we present the observational campaign carried out at ESO NTT and VLT in April and May 2006 to investigate the nature and the structure of the near-Earth object (144898) 2004 VD17. In spite of a great quantity of dynamical information, according to which it will have a close approach with the Earth in the next century, the physical properties of this asteroid are largely unknown. We performed visible and near-infrared photometry and spectroscopy, as well as polarimetric observations. Polarimetric and spectroscopic data allowed us to classify 2004 VD17 as an E-type asteroid. A good agreement was also found with the spectrum of the aubrite meteorite Mayo Belwa. On the basis of the polarimetric albedo (pv=0.45) and of photometric data, we estimated a diameter of about 320 m and a rotational period of about 2 h. The analysis of the results obtained by our complete survey have shown that (144898) 2004 VD17 is a peculiar NEO, since it is close to the breakup limits for fast rotator asteroids, as defined by Pravec and Harris [Pravec, P., Harris, A.W., 2000. Icarus 148, 12-20]. These results suggest that a more robust structure must be expected, as a fractured monolith or a rubble pile in a “strength regime” [Holsapple, K.A., 2002. Speed limits of rubble pile asteroids: Even fast rotators can be rubble piles. In: Workshop on Scientific Requirements for Mitigation of Hazardous Comets and Asteroids, Washington, September, 2002].     

Thermal infrared and optical observations of four near-Earth asteroids

We present thermal infrared photometry and spectrophotometry of four near-Earth asteroids (NEAs), namely (433) Eros, (66063) 1998 RO₁, (137032) 1998 UO₁, and (138258) 2000 GD₂, using the United Kingdom Infrared Telescope (UKIRT) in 2002. For two objects, i.e. (433) Eros and (137032) 1998 UO₁, quasi-simultaneous optical observations were also obtained, using the Jacobus Kapteyn Telescope (JKT). For (127032) 1998 UO₁, we obtain a rotation period P=3.0±0.1 h and an absolute visual magnitude HV=16.7±0.4. The Standard Thermal Model (STM), Fast Rotating Model (FRM) and near-Earth asteroid Thermal Model (NEATM) have been fitted to the IR fluxes to determine effective diameters Deff, geometric albedos pv, and beaming parameters η. The derived values are (433) Eros: Deff=23.3±3.5 km (at lightcurve maximum), pv=0.24±0.07, η=0.95±0.19; (66063) 1998 RO₁: , ; (137032) 1998 UO₁: Deff<1.13 km, pv>0.29; (138258) 2000 GD₂: Deff=0.27±0.04 km, , η=0.74±0.15. (66063) 1998 RO₁ is a binary asteroid from lightcurve characteristics [Pravec, P., and 56 colleagues, 2006. Icarus 181, 63-93] and we estimate the effective diameter of the primary (Dp) and secondary (Ds) components: and . The diameter and albedo of (138258) 2000 GD₂ are consistent with the trend of decreasing diameter for S- and Q-type asteroids found by Delbó et al. [Delbó, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. A possible trend of increasing beaming parameter with diameter for small (less than about 3 km) S- and Q-type asteroids is found.     

Limits on the size and orbit distribution of main belt comets

The first of a new class of objects now known as main belt comets (MBCs) or “activated asteroids” was identified in 1996. The seven known members of this class have orbital characteristics of main belt asteroids yet exhibit dust ejection like comets. In order to constrain their physical and orbital properties we searched the Thousand Asteroid Light Curve Survey (TALCS; Masiero, J.R., Jedicke, R., Durech, J., Gwyn, S., Denneau, L., Larsen, J. [2009]. Icarus 204, 145-171) for additional candidates using two diagnostics: tail and coma detection. This was the most sensitive MBC survey effort to date, extending the search from MBCs with H ∼ 18 (D ∼ 1 km) to MBCs as small as H ∼ 21 (D ∼ 150 m).We fit each of the 924 objects detected by TALCS to a PSF model incorporating both a coma and nuclear component to measure the fractional contribution of the coma to the total surface brightness. We determined the significance of the coma detection using the same algorithm on a sample of null detections of comparable magnitude and rate of motion. We did not identify any MBC candidates with this technique to a sensitivity limit on the order of cometary mass loss rate of about 0.1 kg/s.Our tail detection algorithm relied on identifying statistically significant flux in a segmented annulus around the candidate object. We show that the technique can detect tail activity throughout the asteroid belt to the level of the currently known MBCs. Although we did not identify any MBC candidates with this technique, we find a statistically significant detection of faint activity in the entire ensemble of TALCS asteroids. This suggests that many main belt asteroids are active at very low levels.Our null detection of MBCs allows us to set 90% upper confidence limits on the number distribution of MBCs as a function of absolute magnitude, semi-major axis, eccentricity, and inclination. There are ?400,000 MBCs in the main belt brighter than H_V = 21 (∼150-m in diameter) and the MBC:MBA ratio is ?1:400.We further comment on the ability of observations to meaningfully constrain the snow line’s location. Under some reasonable and simple assumptions we claim 85% confidence that the contemporary snow line lies beyond 2.5 AU.     

Initial sizes of planetesimals and accretion of the asteroids

The present size frequency distribution (SFD) of bodies in the asteroid belt appears to have preserved some record of the primordial population, with an excess of bodies of diameter D ∼ 100 km relative to a simple power law. The survival of Vesta’s basaltic crust also implies that the early SFD had a shallow slope in the range ∼10-100 km. (Morbidelli, A., Bottke, W.F., Nesvorny, D., Levison, H.F. [2009]. Icarus 204, 558-573) were unable to produce these features by accretion from an initial population of km-sized planetesimals. They concluded that bodies with sizes in the range ∼100-1000 km and a SFD similar to the current population were produced directly from solid particles of sub-meter scale, without experiencing accretion through intermediate sizes. We present results of new accretion simulations in the primordial asteroid region. The requisite SFD can be produced from an initial population of planetesimals of sizes ?0.1 km, smaller than the usual assumption of km-sized bodies. The bump at D ∼ 100 km is produced by a transition from dispersion-dominated runaway growth to a regime dominated by Keplerian shear, before the formation of large protoplanetary embryos. Thus, accretion of the asteroids from an initial population of small (sub-km) planetesimals cannot be ruled out.     

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in three primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations (Hubickyj, O., Bodenheimer, P., Lissauer, J.J. [2005]. Icarus 179, 415-431; Lissauer, J.J., Hubickyj, O., D’Angelo, G., Bodenheimer, P. [2009]. Icarus 199, 338-350). These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 M_⊕, for solid surface densities, σ, of 10, 6, and 4 g cm⁻², respectively. For σ = 10 and σ = 6 g cm⁻², respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.     

Redetermination of the space weathering rate using spectra of Iannini asteroid family members

We have obtained moderate S/N (∼85) spectra at a realized resolution of R∼100 for 11 members of the Iannini family, until recently the youngest known family at under 5 million years of age [Nesvorný, D., Bottke, W.F., Levison, H.F., Dones, L., 2003. Astrophys. J. 591, 486-497, 720-771]. The spectra were acquired using the Echellette Spectrograph and Imager in its low-resolution prism mode on the Keck II telescope. The family members belong to the S-complex of asteroids with perhaps some K class members. The Iannini family members' average spectral slope, defined as the slope of the best-fit line constrained to pivot about 1 at 550 nm, is (0.30±0.04)/μm, matching the (0.26±0.03)/μm reported by Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277] using SDSS [Ivezi?, ?., Juri?, M., Lupton, R.H., Tabachnik, S., Quinn, T., 2002. In: Tyson, J.A., Wolff, S. (Eds.), Survey and Other Telescope Technologies and Discoveries. In: Proc. SPIE, vol. 4836. SPIE, Bellingham, pp. 98-103] color photometry. Using our spectra for this family as well as new observations of Karin family members [Vernazza, P., Birlan, M., Rossi, A., Dotto, E., Nesvorný, D., Brunetto, R., Fornasier, S., Fulchignoni, M., Renner, S., 2006. Astron. Astrophys. 460, 945-951] and new classifications of some older families we have revised the space weathering rate of S-complex asteroids originally determined by Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277]. Following Jedicke et al. [Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M., 2004. Nature 429, 275-277] we parameterize the space weathering rate of the principal component color of the spectrum (PC₁), which is correlated with the spectral slope, as PC₁(t)=PC₁(0)+ΔPC₁[1−exp−α(t/τ)]. Our revised rate suggests that the characteristic time scale for space weathering is τ=570±220 Myr and that new S-complex clusters will have an initial color of PC₁(0)=0.31±0.04. The revised time scale is in better agreement with lab measurements and our measurements support the use of space weathering as a dating method. Under the assumption that all the spectra should be identical, since the members all derive from the same parent body and are presumably covered with similar regolith, we combined them to obtain a high-S/N composite spectrum for the family. The combined spectrum is within the S-complex.     

Space weathering of small Koronis family asteroids in the SDSS Moving Object Catalog

Following the work of Rivkin et al. (Rivkin, A.S., Thomas, C.A., Trilling, D.E., Enga, M., Grier, J.A. [2011]. Icarus 211, 1294–1297) and Thomas et al. (Thomas, C.A., Rivkin, A.S, Trilling, D.E., Enga, M., Grier, J.A. [2011a]. Icarus 212, 158–166), we investigate space weathering trends in the Koronis family using the larger sample size of the Sloan Digital Sky Survey Moving Object Catalog. We confirm the trend in spectral slope seen in our earlier work and extend our results by investigating the trend in band depth (i ? z color index) to show that Koronis family asteroids smaller than 4 km show the transition from fresh Q-type to weathered S-type surfaces.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.     

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 formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.     

Size-frequency distributions of fragments from SPH/N-body simulations of asteroid impacts: Comparison with observed asteroid families

We investigate the morphology of size-frequency distributions (SFDs) resulting from impacts into 100-km-diameter parent asteroids, represented by a suite of 161 SPH/N-body simulations conducted to study asteroid satellite formation [Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M., 2004. Icarus 170, 243-257]. The spherical basalt projectiles range in diameter from 10 to 46 km (in equally spaced mass increments in logarithmic space, covering six discrete sizes), impact speeds range from 2.5 to 7 km/s (generally in 1 km/s increments), and impact angles range from 15° to 75° (nearly head-on to very oblique) in 15° increments. These modeled SFD morphologies match very well the observed SFDs of many known asteroid families. We use these modeled SFDs to scale to targets both larger and smaller than 100 km in order to gain insights into the circumstances of the impacts that formed these families. Some discrepancies occur for families with parent bodies smaller than a few tens of kilometers in diameter (e.g., 832 Karin), however, so due caution should be used in applying our results to such small families. We find that ∼20 observed main-belt asteroid families are produced by the catastrophic disruption of D>100 km parent bodies. Using these data as constraints, collisional modeling work [Bottke Jr., W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H.F., 2005b. Icarus 179, 63-94] suggests that the threshold specific energy, , needed to eject 50% of the target body's mass is very close to that predicted by Benz and Asphaug [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20].     

Determination of the complex refractive indices of Titan haze analogs using photothermal deflection spectroscopy

The spectrometers of the Cassini mission to the Saturn system have detected haze layers reaching up to 800 km in Titan’s atmosphere. Knowledge of the complex refractive index (k) of the haze is important for modeling the surface and atmosphere of Titan and retrieving some information about the functional groups present in the aerosols. Plasma discharges or ultraviolet radiation are commonly used to drive the formation of solid organics assumed to be good analogs of the Titan aerosols. [Tran, B.N., Ferris, J.P., Chera, J.J., 2003a. The photochemical formation of a Titan haze analog. Structural analysis by X-ray photoelectron and infrared spectroscopy. Icarus 162, 114-124; Tran, B.N., Force, M., Briggs, R., Ferris J.P., Persans, P., Chera, J.J., 2008. Photochemical processes on Titan: Irradiation of mixtures of gases that simulate Titan’s atmosphere. Icarus 177, 106-115] reported the index of refraction of analogs synthesized by far ultraviolet irradiation of various gas mixtures. k was determined in the 200-800 nm wavelength range from transmission and reflection spectroscopy. However, this technique is limited by (i) uncertainties in the absorption values because of the small amounts of organics available, (ii) light scattering by the surface roughness and particulates in the sample. These limitations prompted us to perform new measurements using photothermal deflection spectroscopy (PDS), a technique based on the conversion of absorbed light into heat in the material of interest. By combining traditional spectroscopy (λ < 500 nm) and PDS (λ > 500 nm), we determined values of k over the 375-1550 nm range. k values as low as 10⁻⁴ above 1000 nm were determined. This is one order of magnitude lower than the measurements generally used as a reference for Titan’s aerosols analogs [Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callicott, T.A., Williams, M.W., 1984. Optical-constants of organic Tholins produced in a simulated Titanian atmosphere—from soft-X-ray to microwave-frequencies. Icarus 60(1), 127-137]. We recommend that these results were used in models to describe the optical properties of the aerosols produced in Titan’s stratosphere.     

Linking the collisional history of the main asteroid belt to its dynamical excitation and depletion

The main belt is believed to have originally contained an Earth mass or more of material, enough to allow the asteroids to accrete on relatively short timescales. The present-day main belt, however, only contains ∼5×¹⁰⁻⁴ Earth masses. Numerical simulations suggest that this mass loss can be explained by the dynamical depletion of main belt material via gravitational perturbations from planetary embryos and a newly-formed Jupiter. To explore this scenario, we combined dynamical results from Petit et al. [Petit, J. Morbidelli, A., Chambers, J., 2001. The primordial excitation and clearing of the asteroid belt. Icarus 153, 338-347] with a collisional evolution code capable of tracking how the main belt undergoes comminution and dynamical depletion over 4.6 Gyr [Bottke, W.F., Durda, D., Nesvorny, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H., 2005. The fossilized size distribution of the main asteroid belt. Icarus 175, 111-140]. Our results were constrained by the main belt's size-frequency distribution, the number of asteroid families produced by disruption events from diameter D>100 km parent bodies over the last 3-4 Gyr, the presence of a single large impact crater on Vesta's intact basaltic crust, and the relatively constant lunar and terrestrial impactor flux over the last 3 Gyr. We used our model to set limits on the initial size of the main belt as well as Jupiter's formation time. We find the most likely formation time for Jupiter was 3.3±2.6 Myr after the onset of fragmentation in the main belt. These results are consistent with the estimated mean disk lifetime of 3 Myr predicted by Haisch et al. [Haisch, K.E., Lada, E.A., Lada, C.J., 2001. Disk frequencies and lifetimes in young clusters. Astrophys. J. 553, L153-L156]. The post-accretion main belt population, in the form of diameter D?1000 km planetesimals, was likely to have been 160±40 times the current main belt's mass. This corresponds to 0.06-0.1 Earth masses, only a small fraction of the total mass thought to have existed in the main belt zone during planet formation. The remaining mass was most likely taken up by planetary embryos formed in the same region. Our results suggest that numerous D>200 km planetesimals disrupted early in Solar System history, but only a small fraction of their fragments survived the dynamical depletion event described above. We believe this may explain the limited presence of iron-rich M-type, olivine-rich A-type, and non-Vesta V-type asteroids in the main belt today. The collisional lifetimes determined for main belt asteroids agree with the cosmic ray exposure ages of stony meteorites and are consistent with the limited collisional evolution detected among large Koronis family members. Using the same model, we investigated the near-Earth object (NEO) population. We show the shape of the NEO size distribution is a reflection of the main belt population, with main belt asteroids driven to resonances by Yarkovsky thermal forces. We used our model of the NEO population over the last 3 Gyr, which is consistent with the current population determined by telescopic and satellite data, to explore whether the majority of small craters (D<0.1-1 km) formed on Mercury, the Moon, and Mars were produced by primary impacts or by secondary impacts generated by ejecta from large craters. Our results suggest that most small craters formed on these worlds were a by-product of secondary rather than primary impacts.     

How the Enceladus dust plume feeds Saturn’s E ring

Pre-Cassini models of Saturn’s E ring [Horányi, M., Burns, J., Hamilton, D., 1992. Icarus 97, 248-259; Juhász, A., Horányi, M., 2002. J. Geophys. Res. 107, 1-10] failed to reproduce its peculiar vertical structure inferred from Earth-bound observations [de Pater, I., Martin, S.C., Showalter, M.R., 2004. Icarus 172, 446-454]. After the discovery of an active ice-volcanism of Saturn’s icy moon Enceladus the relevance of the directed injection of particles for the vertical ring structure of the E ring was swiftly recognised [Juhász, A., Horányi, M., Morfill, G.E., 2007. Geophys. Res. Lett. 34, L09104; Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., Economou, T., Schmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. However, simple models for the delivery of particles from the plume to the ring predict a too small vertical ring thickness and overestimate the amount of the injected dust.Here we report on numerical simulations of grains leaving the plume and populating the dust torus of Enceladus. We run a large number of dynamical simulations including gravity and Lorentz force to investigate the earliest phase of the ring particle life span. The evolution of the electrostatic charge carried by the initially uncharged grains is treated selfconsistently. Freshly ejected plume particles are moving in almost circular orbits because the Enceladus orbital speed exceeds the particles’ ejection speeds by far. Only a small fraction of grains that leave the Hill sphere of Enceladus survive the next encounter with the moon. Thus, the flux and size distribution of the surviving grains, replenishing the ring particle reservoir, differs significantly from the flux and size distribution of the particles freshly ejected from the plume. Our numerical simulations reproduce the vertical ring profile measured by the Cassini Cosmic Dust Analyzer (CDA) [Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., EconoDmou, T., Smchmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. From our simulations we calculate the deposition rates of plume particles hitting Enceladus’ surface. We find that at a distance of 100 m from a jet a 10 m sized ice boulder should be covered by plume particles in 10⁵-10⁶ years.     

In search of water vapor on Jupiter: Laboratory measurements of the microwave properties of water vapor under simulated jovian conditions

Detection and measurement of atmospheric water vapor in the deep jovian atmosphere using microwave radiometry has been discussed extensively by Janssen et al. (Janssen, M.A., Hofstadter, M.D., Gulkis, S., Ingersoll, A.P., Allison, M., Bolton, S.J., Levin, S.M., Kamp, L.W. [2005]. Icarus 173 (2), 447-453.) and de Pater et al. (de Pater, I., Deboer, D., Marley, M., Freedman, R., Young, R. [2005]. Icarus 173 (2), 425-447). The NASA Juno mission will include a six-channel microwave radiometer system (MWR) operating in the 1.3-50 cm wavelength range in order to retrieve water vapor abundances from the microwave signature of Jupiter (see, e.g., Matousek, S. [2005]. The Juno new frontiers mission. Tech. Rep. IAC-05-A3.2.A.04, California Institute of Technology). In order to accurately interpret data from such observations, nearly 2000 laboratory measurements of the microwave opacity of H₂O vapor in a H₂/He atmosphere have been conducted in the 5-21 cm wavelength range (1.4-6 GHz) at pressures from 30 mbars to 101 bars and at temperatures from 330 to 525 K. The mole fraction of H₂O (at maximum pressure) ranged from 0.19% to 3.6% with some additional measurements of pure H₂O. These results have enabled development of the first model for the opacity of gaseous H₂O in a H₂/He atmosphere under jovian conditions developed from actual laboratory data. The new model is based on a terrestrial model of Rosenkranz et al. (Rosenkranz, P.W. [1998]. Radio Science 33, 919-928), with substantial modifications to reflect the effects of jovian conditions. The new model for water vapor opacity dramatically outperforms previous models and will provide reliable results for temperatures from 300 to 525 K, at pressures up to 100 bars and at frequencies up to 6 GHz. These results will significantly reduce the uncertainties in the retrieval of jovian atmospheric water vapor abundances from the microwave radiometric measurements from the upcoming NASA Juno mission, as well as provide a clearer understanding of the role deep atmospheric water vapor may play in the decimeter-wavelength spectrum of Saturn.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.