We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

We develop an evolutionary model of the long-period comet (LPC) population, starting from their birthplace in a massive trans-Neptunian disk that was dispersed by migrating giant planets. Most comets that remain bound to the Solar system are stored in the Oort cloud. Galactic tides and passing stars make some of these bodies evolve into observable comets in the inner Solar system. Our approach models each step in a full-fledged numerical framework. Subsequent analysis consists of applying plausible fading models and computing the original orbits to compare with observations. Our results match the observed semimajor axis distribution of LPCs when Whipple's power-law fading scheme with an exponent $\kappa= 0.6^{+0.1}_{-0.2}$ is adopted. The cumulative perihelion ($q$) distribution is fit well by a linear increase plus a weak quadratic term. Beyond $q = 15$~au, however, the population increases steeply and the isotropy of LPC orbital planes breaks. We find tentative evidence from the perihelion distribution of LPCs that the returning comets are depleted in supervolatiles and become active due to water ice sublimation for $q\leq 3$ au. Using an independent calibration of the population of the initial disk, our predicted LPC flux is smaller than observations suggest by a factor of $\simeq 2$. Current data only characterize comets from the outer Oort cloud (semimajor axes $\gtrsim 10^4$~au). A true boost in understanding the Oort cloud's structure should result from future surveys when they detect LPCs with perihelia beyond $15$~au. Our results provide observational predictions of what can be expected from these new data.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

Asteroid families are formed as the result of collisions. Large fragments are ejected with speeds of the order of the escape velocity from the parent body. After the family formation, the fragments' orbits evolve in the space of proper elements because of gravitational and non-gravitational perturbations, such as the Yarkovsky effect. Disentangling the contribution to the current orbital position of family members caused by the initial ejection velocity field and the subsequent orbital evolution is usually a difficult task. Among the more than 100 asteroid families currently known, some interact with secular resonances. Linear secular resonances occur when there is a commensurability between the precession frequency of the longitude of the pericenter (g) or of the longitude of node (s) of an asteroid and a planet, or a massive asteroid. The linear secular resonance most effective in increasing an asteroid eccentricity is the ${\nu}_6$, that corresponds to a commensurability between the precession frequency g of an asteroid and Saturn's $g_6$. Non-linear secular resonances involve commensurabilities of higher order, and can often be expressed as combinations of linear secular resonances. This is the case, for instance, of the $z_k=k(g-g_6)+(s-s_6)$ resonances. Asteroid families that are crossed by secular resonances are of particular interest in dynamical astronomy. First, they often provide a clear evidence of asteroid orbit evolution due to the Yarkovsky effect. Second, conserved quantities of secular dynamics can be used to set valuable constraints on the magnitude of the original ejection velocity field. Finally, by changing the value of inclination of family members nodal secular resonances with massive asteroids or dwarf planets can cause the i distribution to become more and more leptokurtic (i.e., more peaked and with larger tails than that of a Gaussian distribution).

The Nele group, formerly known as the Iannini family, is one of the youngest asteroid families in the main belt. Previously, it has been noted that the pericenter longitudes $\varpi$ and nodal longitudes $\Omega$ of its largest member asteroids are clustered at the present time, therefore suggesting that the collisional breakup of parent body must have happened recently. Here we verify this conclusion by detailed orbit-propagation of a synthetic Nele family and show that the current level of clustering of secular angles of the largest Nele family members requires an approximate age limit of $4.5$ Myr. Additionally, we make use of an updated and largely extended Nele membership to obtain, for the first time, an age estimate of this family using the Backward Integration Method (BIM). Convergence of the secular angles in a purely gravitational model and in a model including the non-gravitational forces caused by the Yarkovsky effect are both compatible with an age younger than $7$ Myr. More accurate determination of the Nele family age would require additional data about the spin state of its members.

Context: Multiple stellar systems are unique laboratories for astrophysics. Analysis of their orbital dynamics may reveal invaluable information about the physical properties of the participating stars. Unfortunately, there are only a few known and well described multiple systems, this is even more so for systems located outside the Milky Way galaxy. A particularly interesting situation occurs when the inner binary in a compact triple system is eclipsing. This is because the stellar interaction, typically resulting in precession of orbital planes, may be observable as a variation of depth of the eclipses on a long timescale. Aims: We aim to present a novel method to determine compact triples using publicly available photometric data from large surveys. Here we apply it to eclipsing binaries (EBs) in Magellanic Clouds from OGLE III database. Methods: We analyzed light curves (LCs) of 26121 LMC and 6138 SMC EBs with the goal to identify those for which the orbital inclination varies in time. Archival LCs of the selected systems, when complemented by our own observations with Danish 1.54m telescope, were thoroughly analyzed using the PHOEBE program. Time dependence of the EB's inclination was described using the theory of orbital-plane precession. By observing the parameter-dependence of the precession rate, we were able to constrain the third companion mass and its orbital period around EB. Results: We identified 58 candidates of new compact triples in Magellanic Clouds. This is the largest published sample of such systems so far. Eight of them were analyzed thoroughly and physical parameters of inner binary were determined together with an estimation of basic characteristics of the third star. These data may provide important clues about stellar formation mechanisms for objects with different metalicity than found in our galactic neighborhood.

The age of a young asteroid family can be determined by tracking the orbits of family members backward in time and showing that they converge at some time in the past. Here we consider the Veritas family. We find that the membership of the Veritas family increased enormously since the last detailed analysis of the family. Using backward integration, we confirm the convergence of nodal longitudes $\Omega$, and, for the first time, also obtain a simultaneous convergence of pericenter longitudes $\varpi$. The Veritas family is found to be $8.23^{+0.37}_{-0.31}$~Myr old. To obtain a tight convergence of $\Omega$ and $\varpi$, as expected from low ejection speeds of fragments, the Yarkovsky effect needs to be included in the modeling of the past orbital histories of Veritas family members. Using this method, we compute the Yarkovsky semi-major axis drift rates, ${\rm d}a/{\rm d}t$, for 274 member asteroids. The distribution of ${\rm d}a/{\rm d}t$ values is consistent with a population of C-type objects with low densities and low thermal conductivities. The accuracy of individual ${\rm d}a/{\rm d}t$ measurements is limited by the effect of close encounters of member asteroids to (1) Ceres and other massive asteroids, which cannot be evaluated with confidence.

Context: Large all-sky surveys provide us with a lot of photometric data that are sparse in time (typically few measurements per night) and can be potentially used for the determination of shapes and rotational states of asteroids. The method which is generally used to derive these parameters is the lightcurve inversion. However, for most asteroids their sparse data are not accurate enough to derive a unique model and the lightcurve inversion method is thus not very efficient. Aims: To fully utilize photometry sparse in time, we developed a new simplified model and applied it on the data from the Lowell photometric database. Our aim was to derive spin axis orientations and shape elongations of asteroids and to reconstruct distributions of these parameters for selected subpopulations to find if there are some differences. Methods: We model asteroids as geometrically scattering triaxial ellipsoids. Observed values of mean brightness and the dispersion of brightness are compared with computed values obtained from the parameters of the model -- ecliptical longitude $\lambda$ and latitude $\beta$ of the pole and the ratios $a/b$, $b/c$ of axes of the ellipsoid. These parameters are optimized to get the best agreement with the observation. Results: We found that the distribution of $\lambda$ for main-belt asteroids is not uniform (in agreement with findings of Bowell et al., 2014, M&PS, 49, 95) and is dependent on the inclination of orbit. Surprisingly, the non-uniformity of $\lambda$ distribution is larger for asteroids residing on low-inclination orbits. We also studied distributions of $a/b$ for several groups of asteroids and found that small asteroids ($D<25\,$km) are on average more elongated than large ones.

The Karin cluster is a young asteroid family thought to have formed only $\simeq 5.75$~My ago. The young age can be demonstrated by numerically integrating the orbits of Karin cluster members backward in time and showing the convergence of the perihelion and nodal longitudes (as well as other orbital elements). Previous work has pointed out that the convergence is not ideal if the backward integration only accounts for the gravitational perturbations from the Solar System planets. It improves when the thermal radiation force known as the Yarkovsky effect it is accounted for. This argument can be used to estimate the spin obliquities of the Karin cluster members. Here we take advantage of the fast growing membership of the Karin cluster and show that the obliquity distribution of diameter $D\simeq 1-2$ km Karin asteroids is bimodal, as expected if the YORP effect acted to move obliquities toward the extreme values ($0^\circ$ or $180^\circ$). The measured magnitude of the effect is consistent with the standard YORP model. The surface thermal conductivity is inferred to be $0.07$-0.2 W m$^{-1}$ K$^{-1}$ (thermal inertia $\simeq 300-500$ J m$^{-2}$ K$^{-1}$s$^{-1/2}$). We find that the strength of the YORP effect is roughly $\simeq 0.7$ of the nominal strength obtained for a collection of random Gaussian spheroids. These results are consistent with a surface composed of rough, rocky regolith. The obliquity values predicted here for 480 members of the Karin cluster can be validated by the light-curve inversion method.

We explored the statistical and compositional link between Chelyabinsk meteoroid and potentially hazardous asteroid (86039) 1999 NC43 to investigate their proposed relation proposed by Borovi\v{c}ka et al. (2013). Using detailed computation we confirm that the orbit of the Chelyabinsk impactor is anomalously close to 1999 NC43. We find about (1-3) x 10-4 likelihood of that to happen by chance. Taking the standpoint that the Chelyabinsk impactor indeed separated from 1999 NC43 by a cratering or rotational fission event, we run a forward probability calculation, which is an independent statistical test. However, we find this scenario is unlikely at the about (10-3 -10-2) level. We also verified compositional link between Chelyabinska and 1999NC43. Mineralogical analysis of Chelyabinsk (LL chondrite) and (8) Flora (the largest member of the presumed LL chondrite parent family) shows that their olivine and pyroxene chemistries are similar to LL chondrites. Similar analysis of 1999 NC43 shows that its olivine and pyroxene chemistries are more similar to L chondrites than LL chondrites (like Chelyabinsk). We also took photometric observations of 1999 NC43 over 54 nights during two apparitions (2000, 2014). The lightcurve of 1999 NC43 resembles simulated lightcurves of tumblers in Short-Axis Mode with the mean wobbling angle 20-30 deg. While, a mechanism of the non-principal axis rotation excitation is unclear, we can rule out the formation of asteroid in disruption of its parent body as a plausible cause, as it is unlikely that the rotation of an asteroid fragment from catastrophic disruption would be nearly completely halted. Considering all these facts, we find the proposed link between the Chelyabinsk meteoroid and the asteroid 1999 NC43 to be unlikely.

The Design Reference Asteroid (DRA) is a compilation of all that is known about the OSIRIS-REx mission target, asteroid (101955) Bennu. It contains our best knowledge of the properties of Bennu based on an extensive observational campaign that began shortly after its discovery, and has been used to inform mission plan development and flight system design. The DRA will also be compared with post-encounter science results to determine the accuracy of our Earth-based characterization efforts. The extensive observations of Bennu in 1999 has made it one of the best-characterized near-Earth asteroids. Many physical parameters are well determined, and span a number of categories: Orbital, Bulk, Rotational, Radar, Photometric, Spectroscopic, Thermal, Surface Analog, and Environment Properties. Some results described in the DRA have been published in peer-reviewed journals while others have been reviewed by OSIRIS-REx Science Team members and/or external reviewers. Some data, such as Surface Analog Properties, are based on our best knowledge of asteroid surfaces, in particular those of asteroids Eros and Itokawa. This public release of the OSIRIS-REx Design Reference Asteroid is a annotated version of the internal OSIRIS-REx document OREX-DOCS-04.00-00002, Rev 9 (accepted by the OSIRIS-REx project on 2014-April-14). The supplemental data products that accompany the official OSIRIS-REx version of the DRA are not included in this release. We are making this document available as a service to future mission planners in the hope that it will inform their efforts.

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m.

The rotational-fission of a rubble-pile asteroid can result in an "asteroid pair", two un-bound asteroids sharing similar orbits. This mechanism might exposes material that previously had never have been exposed to the weathering conditions of space. Therefore, the surfaces of asteroid pairs offer the opportunity to observe non-weathered fresh spectra. We report near-IR spectroscopic observations of 31 asteroids in pairs. We analyze their spectral slopes, 1 {\mu}m absorption band, taxonomy, and estimate the time elapsed since their separation. Analyzing the 19 S-complex objects in our sample, we find two fresh Q-type asteroids that are the first of their kind to be observed in the main-belt over the full visible and near-IR range. This solidly demonstrates that Q-type objects are not limited to the NEA population. The pairs in our sample present a range of fresh and weathered surfaces with no clear evidence for a correlation with the ages of the pairs. However, our sample includes old pairs (1 to 2 My) that present low spectral slopes. This illustrates a timescale of at least ~2 My before an object develops high spectral slope that is typical for S-type asteroids. We discuss mechanisms that explain the existence of weathered pairs with young dynamical ages and find that the "secondary fission" model (Jacobson & Scheeres 2011) is the most robust with our observations since: 1) the secondary members in our sample present fresh parameters that tend to be fresher than their weathered primaries; 2) most of the fresh pairs in our sample have low size ratios between the secondary and the primary; 3) 33% of the primaries in our sample are fresh, similar to the prediction set by this model; 4) known satellites orbit two of the pairs in our sample with low size ratio and fresh surface; 5) there is no correlation between the weathering state and the primary shape as predicted by other models.

The so-called Nysa-Polana complex of asteroids is a diverse and widespread group. It appears to be two overlapping families of different asteroid taxonomies: (44) Nysa is an E-type asteroid with the lowest number in the midst of a predominantly S-type cluster and (142) Polana is a B-type asteroid near the low-albedo B- and C-type cluster. Using the data from the Wide-field Infrared Survey Explorer (WISE) mission we have re-analyzed the region around the Nysa-Polana complex in the inner Main Belt, focusing on the low-albedo population. (142) Polana does not appear to be a member of the family of low-albedo asteroids in the Nysa-Polana complex. Rather, the largest is asteroid (495) Eulalia. This asteroid has never before been linked to this complex for an important dynamical reason: it currently has a proper eccentricity slightly below the range of most of the family members. However, its orbit is very close to the 3:1 mean motion resonance with Jupiter and is in a weak secular resonance. We show that its osculating eccentricity varies widely on short timescales and the averaged value diffuses over long timescales. The diffusive orbit, low-albedo, taxonomic similarity and semimajor axis strongly suggests that despite its current proper eccentricity, (495) Eulalia could have recently been at an orbit very central to the family. Hierarchical Clustering Method tests confirm that at an eccentricity of e=0.15, (495) Eulalia could be the parent of the family. The ``Eulalia family'' was formed between 900--1500 Myr ago, and likely resulted from the breakup of a 100--160 km parent body. There is also compelling evidence for an older and more widespread primitive family in the same region of the asteroid belt parented by asteroid (142) Polana.

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.

We have identified and quantified semi-major axis drifts in Near-Earth Asteroids (NEAs) by performing orbital fits to optical and radar astrometry of all numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the most reliable and strongest drift rates. Our selection criteria include a Yarkovsky sensitivity metric that quantifies the detectability of semi-major axis drift in any given data set, a signal-to-noise metric, and orbital coverage requirements. In 42 cases, the observed drifts (~10^-3 AU/Myr) agree well with numerical estimates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is the dominant non-gravitational process affecting these orbits, and allows us to derive constraints on asteroid physical properties. In 12 cases, the drifts exceed nominal Yarkovsky predictions, which could be due to inaccuracies in our knowledge of physical properties, faulty astrometry, or modeling errors. If these high rates cannot be ruled out by further observations or improvements in modeling, they would be indicative of the presence of an additional non-gravitational force, such as that resulting from a loss of mass of order a kilogram per second. We define the Yarkovsky efficiency f_Y as the ratio of the change in orbital energy to incident solar radiation energy, and we find that typical Yarkovsky efficiencies are ~10^-5.

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

Asteroids residing in the first-order mean motion resonances with Jupiter hold important information about the processes that set the final architecture of giant planets. Here we revise current populations of objects in the J2/1 (Hecuba-gap group), J3/2 (Hilda group) and J4/3 (Thule group) resonances. The number of multi-opposition asteroids found is 274 for J2/1, 1197 for J3/2 and 3 for J4/3. Using both hierarchical clustering technique and colour identification we characterise a collisionally-born asteroid family around the object (1911) Schubart in the J3/2 resonance. There is also a looser cluster around the largest asteroid (153) Hilda. Using N-body numerical simulations we prove that the Yarkovsky effect (infrared thermal emission from the surface of asteroids) causes a systematic drift in eccentricity for resonant asteroids, while their semimajor axis is almost fixed due to the strong coupling with Jupiter. This is a different mechanism from main belt families, where the Yarkovsky drift affects basically the semimajor axis. We use the eccentricity evolution to determine the following ages: (1.7+-0.7) Gyr for the Schubart family and >~4 Gyr for the Hilda family. We also find that collisionally-born clusters in the J2/1 resonance would efficiently dynamically disperse. The steep size distribution of the stable population inside this resonance could thus make sense if most of these bodies are fragments from an event older than ~1 Gyr. Finally, we test stability of resonant populations during Jupiter's and Saturn's crossing of their mutual mean motion resonances. In particular we find primordial objects in the J3/2 resonance were efficiently removed from their orbits when Jupiter and Saturn crossed their 1:2 mean motion resonance.

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10^6 years and to newly identified asteroid pairs with ages <5x10^5 years. We have identified three S-complex asteroid clusters with ages in the range 10^{5-6} years. The average color of the objects in these clusters agree with the prediction of Willman et al., 2008. SDSS photometry of the members of very young asteroid pairs with ages <10^5 years was used to determine their taxonomy. The average color of the S-complex pairs is PC_1=0.49+/-0.03, over 5-sigma redder than predicted by Willman et al., 2008. Therefore, the most likely pair formation mechanism is gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. In this case our color measurement allows us to set an upper limit of ~64% on the disturbed surface portion. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic timescales of tau_w=960+/-160My and tau_g=2000+/-290My respectively. The first principal component color for fresh S-complex material is 0.37+/-0.01 while the maximum amount of local reddening is 0.33+/-0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of tau_g~270My based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a `honeycomb' mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.