We report Arecibo S-band (2380 MHz; 12.6 cm) radar observations of near-Earth asteroid (3200) Phaethon during the December 2017 apparition when Phaethon passed within 0.07 au of Earth. Radar images with a resolution of 75 m per pixel reveal a roughly spheroidal shape more than 6 km in diameter at the equator with several discernible surface features hundreds of meters in extent. These include a possible crater more than 1 km across located below 30$^{\circ}$ latitude and a roughly 600-m radar-dark region near one of the poles. Overall, the radar images of Phaethon are reminiscent of those of (101955) Bennu, target of the OSIRIS-REx mission. As such, the shape of Phaethon is suspected to have an equatorial ridge similar to the top-shaped models of several other radar-observed near-Earth asteroids as well as the optical images of (162173) Ryugu returned by the Hayabusa2 spacecraft. Preliminary analysis of the radar data finds no satellites and gives no indication of a dusty coma at the time of these observations.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion.

Near-Earth asteroid (NEA) 1566 Icarus ($a=1.08$ au, $e=0.83$, $i=22.8^{\circ}$) made a close approach to Earth in June 2015 at 22 lunar distances (LD). Its detection during the 1968 approach (16 LD) was the first in the history of asteroid radar astronomy. A subsequent approach in 1996 (40 LD) did not yield radar images. We describe analyses of our 2015 radar observations of Icarus obtained at the Arecibo Observatory and the DSS-14 antenna at Goldstone. These data show that the asteroid is a moderately flattened spheroid with an equivalent diameter of 1.44 km with 18% uncertainties, resolving long-standing questions about the asteroid size. We also solve for Icarus' spin axis orientation ($\lambda=270^{\circ}\pm10^{\circ}, \beta=-81^{\circ}\pm10^{\circ}$), which is not consistent with the estimates based on the 1968 lightcurve observations. Icarus has a strongly specular scattering behavior, among the highest ever measured in asteroid radar observations, and a radar albedo of $\sim$2\%, among the lowest ever measured in asteroid radar observations. The low cross-section suggests a high-porosity surface, presumably related to Icarus' cratering, spin, and thermal histories. Finally, we present the first use of our orbit determination software for the generation of observational ephemerides, and we demonstrate its ability to determine subtle perturbations on NEA orbits by measuring Icarus' orbit-averaged drift in semi-major axis ($(-4.62\pm0.48) \times 10^{-4}$ au/My, or $\sim$60 m per revolution). Our Yarkovsky rate measurement resolves a discrepancy between two published rates that did not include the 2015 radar astrometry.

We evaluated the planetary radar capabilities at Arecibo, the Goldstone 70-m DSS-14 and 34-m DSS-13 antennas, the 70-m DSS-43 antenna at Canberra, the Green Bank Telescope, and the Parkes Radio Telescope in terms of their relative sensitivities and the number of known near-Earth asteroids (NEAs) detectable per year in monostatic and bistatic configurations. In the 2015 calendar year, monostatic observations with Arecibo and DSS-14 were capable of detecting 253 and 131 NEAs respectively, with signal-to-noise ratios (SNRs) greater than 30/track. Combined, the two observatories were capable of detecting 276 NEAs. Of these, Arecibo detected 77 and Goldstone detected 32, or 30% and 24% the numbers that were possible. The two observatories detected an additional 18 and 7 NEAs respectively, with SNRs of less than 30/track. This indicates that a substantial number of potential targets are not being observed. The bistatic configuration with DSS-14 transmitting and the Green Bank Telescope receiving was capable of detecting about 195 NEAs, or ~50% more than with monostatic observations at DSS-14. Most of the detectable asteroids were targets of opportunity that were discovered less than 15 days before the end of their observing windows. About 50% of the detectable asteroids have absolute magnitudes > 25, which corresponds diameters < ~30 m.

An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes.

Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two components and to improve the estimates of the mutual orbit, component masses, and spin periods. The primary has a sidereal spin period of 2.7745 +/- 0.0007 h and is roughly spheroidal with an equivalent diameter of 863 m +/- 5 %. It has a mass of 4.656 +/- 0.43 x 10^11 kg and a density of 1381 +/- 244 kg m^{-3}. It exhibits an equatorial ridge similar to the (66391) 1999 KW4 primary, however the equatorial ridge in this case is not as regular and has a ~300 m diameter concavity on one side. The secondary has a sidereal spin period of 1.77 +/- 0.02 days commensurate with the orbital period. The secondary is slightly elongated and has overall dimensions of 377 x 314 x 268 m (6 % uncertainties). Its mass is 0.178 +/- 0.021 x 10^{11} kg and its density is 1047 +/- 230 kg m^{-3}. The mutual orbit has a semi-major axis of 2.659 +/- 0.08 km, an eccentricity of 0.019 +/- 0.01, and a period of 1.7556 +/- 0.0015 days. The normalized total angular momentum of this system exceeds the amount required for the expected spin-up formation mechanism. An increase of angular momentum from non-gravitational forces after binary formation is a possible explanation. The two components have similar radar reflectivity, suggesting a similar composition consistent with formation by spin-up. The secondary appears to exhibit a larger circular polarization ratio than the primary, suggesting a rougher surface or subsurface at radar wavelength scales.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar systems over a period of 12 days during its close approach to the Earth in February 2012. We obtained continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion of the radar images yields a detailed shape model with an effective spatial resolution of 100 m. The asteroid has overall dimensions of 2.6 km X 2.2 km X 2.1 km (5% uncertainties) and a surface rich with kilometer-scale ridges and concavities. This size, combined with absolute magnitude measurements, implies an extremely low albedo (~2%). It is a principal axis rotator and spins in a retrograde manner with a sidereal spin period of 8.96 +/- 0.01 hours. In terms of gravitational slopes evaluated at scales of 100 m, the surface seems mostly relaxed with over 99% of the surface having slopes less than 30 degrees, but there are some outcrops at the north pole that may have steeper slopes. Our precise measurements of the range and velocity of the asteroid, combined with optical astrometry, enables reliable trajectory predictions for this potentially hazardous asteroid in the interval 460-2813.

Three-body model fits to Arecibo and Goldstone radar data reveal the nature of two near-Earth asteroid triples. Triple-asteroid system 2001 SN263 is characterized by a primary of ~10^13 kg, an inner satellite ~1% as massive orbiting at ~3 primary radii in ~0.7 days, and an outer satellite ~2.5% as massive orbiting at ~13 primary radii in ~6.2 days. 1994 CC is a smaller system with a primary of mass ~2.6 \times 10^11 kg and two satellites ~2% and ~1% as massive orbiting at distances of ~5.5 and ~19.5 primary radii. Their orbital periods are ~1.2 and ~8.4 days. Examination of resonant arguments shows that the satellites are not currently in a mean-motion resonance. Precession of the apses and nodes are detected in both systems (2001 SN263 inner body: d{\varpi}/dt ~1.1 deg/day, 1994 CC inner body: d{\varpi}/dt ~ -0.2 deg/day), which is in agreement with analytical predictions of the secular evolution due to mutually interacting orbits and primary oblateness. Nonzero mutual inclinations between the orbital planes of the satellites provide the best fits to the data in both systems (2001 SN263: ~14 degrees, 1994 CC: ~16 degrees). Our best-fit orbits are consistent with nearly circular motion, except for 1994 CC's outer satellite which has an eccentric orbit of e ~ 0.19. We examine several processes that can generate the observed eccentricity and inclinations, including the Kozai and evection resonances, past mean-motion resonance crossings, and close encounters with terrestrial planets. In particular, we find that close planetary encounters can easily excite the eccentricities and mutual inclinations of the satellites' orbits to the currently observed values.

We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notably smooth on decameter scales. EV5's radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.