Determining the size and orbital distribution of the population of near-Earth asteroids (NEAs) is the focus of intense research, with the most recent models converging to a population of approximately $1000$ NEAs larger than 1 km and up to approximately $10^9$ NEAs with absolute magnitude $H<30$. We present an analysis of the combined observations of nine of the leading asteroid surveys over the past two decades, and show that for an absolute magnitude $H<17.75$, which is often taken as proxy for an average diameter larger than 1 km, the population of NEAs is $920\pm10$, lower than other recent estimates. The population of small NEAs is estimated at $(4 \pm 1)\times 10^8$ for $H<30$, and the number of decameter NEAs is lower than other recent estimates. This population tracks accurately the orbital distribution of recently discovered large NEAs, and produces an estimated Earth impact rate for small NEAs in good agreement with the bolide data.

A comprehensive characterization of the detection efficiency of nine of the major asteroid surveys that have been active over the past two decades is presented. The detection efficiency is estimated on a nightly basis by comparing the detected asteroids with the complete catalog of known asteroids propagated to the same observing epoch. Results include a nightly estimate of the detection efficiency curves as a function of apparent magnitude and apparent velocity of the asteroids, as well as a cumulative analysis to estimate the overall performance of each survey. The limiting magnitude distribution is estimated for each survey, and it is then modeled as a function of telescope aperture to obtain an estimate over a wide range of apertures.

Although the nucleus of comet C/2013 A1 (Siding Spring) will safely pass Mars in October 2014, the dust in the coma and tail will more closely approach the planet. Using a dynamical model of comet dust, we estimate the impact fluence. Based on our nominal model no impacts are expected at Mars. Relaxing our nominal model's parameters, the fluence is no greater than ~10^-7 grains/m^2 for grain radii larger than 10 {\mu}m. Mars orbiting spacecraft are unlikely to be impacted by large dust grains, but Mars may receive as many as ~10^7 grains, or ~100 kg of total dust. We also estimate the flux of impacting gas molecules commonly observed in comet comae.

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

Comet C/2013 A1 (siding Spring) will experience a high velocity encounter with Mars on October 19, 2014 at a distance of 135,000 km +- 5000 km from the planet center. We present a comprehensive analysis of the trajectory of both the comet nucleus and the dust tail. The nucleus of C/2013 A1 cannot impact on Mars even in the case of unexpectedly large nongravitational perturbations. Furthermore, we compute the required ejection velocities for the dust grains of the tail to reach Mars as a function of particle radius and density and heliocentric distance of the ejection. A comparison between our results and the most current modeling of the ejection velocities suggests that impacts are possible only for millimeter to centimeter size particles released more than 13 au from the Sun. However, this level of cometary activity that far from the Sun is considered extremely unlikely. The arrival time of these particles spans a 20-minute time interval centered at October 19, 2014 at 20:09 TDB, i.e., around the time that Mars crosses the orbital plane of C/2013 A1. Ejection velocities larger than currently estimated by a factor >2 would allow impacts for smaller particles ejected as close as 3 au from the Sun. These particles would reach Mars from 43 to 130 min after the nominal close approach epoch of the purely gravitational trajectory of the nucleus.

The hydrostatic equilibrium of multi-layer bodies lacks a satisfactory theoretical treatment despite its wide range of applicability. Here we show that by using the exact analytical potential of homogeneous ellipsoids we can obtain recursive analytical solutions and an exact numerical method for the hydrostatic equilibrium shape problem of multi-layer planets and synchronous moons. The recursive solutions rely on the series expansion of the potential in terms of the polar and equatorial shape eccentricities, while the numerical method uses the exact potential expression. These solutions can be used to infer the interior structure of planets and synchronous moons from the observed shape, rotation, and gravity. When applied to dwarf planet Ceres, we show that it is most likely a differentiated body with an icy crust of equatorial thickness 30-90 km and a rocky core of density 2.4-3.1 g/cm$^3$. For synchronous moons, we show that the $J_2/C_{22} \simeq 10/3$ and the $(b-c)/(a-c) \simeq 1/4$ ratios have significant corrections of order $\Omega^2/(\pi G \rho)$, with important implications on how their gravitational coefficients are determined from flyby radio science data and on how we assess their hydrostatic equilibrium state.

Gravity inversion allows us to constrain the interior mass distribution of a planetary body using the observed shape, rotation, and gravity. Traditionally, techniques developed for gravity inversion can be divided into Monte Carlo methods, matrix inversion methods, and spectral methods. Here we employ both matrix inversion and Monte Carlo in order to explore the space of exact solutions, in a method which is particularly suited for arbitrary shape bodies. We expand the mass density function using orthogonal polynomials, and map the contribution of each term to the global gravitational field generated. This map is linear in the density terms, and can be pseudo-inverted in the under-determined regime using QR decomposition, to obtain a basis of the affine space of exact interior structure solutions. As the interior structure solutions are degenerate, assumptions have to be made in order to control their properties, and these assumptions can be transformed into scalar functions and used to explore the solutions space using Monte Carlo techniques. Sample applications show that the range of solutions tend to converge towards the nominal one as long as the generic assumptions made are correct, even in the presence of moderate noise. We present the underlying mathematical formalism and an analysis of how to impose specific features on the global solution, including uniform solutions, gradients, and layered models. Analytical formulas for the computation of the relevant quantities when the shape is represented using several common methods are included in the Appendix.

Dawn is the first NASA mission to operate in the vicinity of the two most massive asteroids in the main belt, Ceres and Vesta. This double-rendezvous mission is enabled by the use of low-thrust solar electric propulsion. Dawn will arrive at Vesta in 2011 and will operate in its vicinity for approximately one year. Vesta's mass and non-spherical shape, coupled with its rotational period, presents very interesting challenges to a spacecraft that depends principally upon low-thrust propulsion for trajectory-changing maneuvers. The details of Vesta's high-order gravitational terms will not be determined until after Dawn's arrival at Vesta, but it is clear that their effect on Dawn operations creates the most complex operational environment for a NASA mission to date. Gravitational perturbations give rise to oscillations in Dawn's orbital radius, and it is found that trapping of the spacecraft is possible near the 1:1 resonance between Dawn's orbital period and Vesta's rotational period, located approximately between 520 and 580 km orbital radius.This resonant trapping can be escaped by thrusting at the appropriate orbital phase. Having passed through the 1:1 resonance, gravitational perturbations ultimately limit the minimum radius for low-altitude operations to about 400 km,in order to safely prevent surface impact. The lowest practical orbit is desirable in order to maximize signal-to-noise and spatial resolution of the Gamma-Ray and Neutron Detector and to provide the highest spatial resolution observations by Dawn's Framing Camera and Visible InfraRed mapping spectrometer. Dawn dynamical behavior is modeled in the context of a wide range of Vesta gravity models. Many of these models are distinguishable during Dawn's High Altitude Mapping Orbit and the remainder are resolved during Dawn's Low Altitude Mapping Orbit, providing insight into Vesta's interior structure.

We derive an expression for the mutual gravitational force and torque of two bodies having arbitrary shapes and mass distributions, as an expansion in power series of their products of inertia and of the relative coordinates of their centres of mass. The absolute convergence of all the power series developed is rigorously demonstrated. The absence of transcendental functions makes this formalism suitable for fast numerical applications. The products of inertia used here are directly related to the spherical harmonics coefficients, and we provide a detailed analysis of this relationship.

Asteroids can be eclipsed by other bodies in the Solar System, but no direct observation of an asteroid eclipse has been reported to date. We describe a statistical method to predict an eclipse for an asteroid based on the analysis of the orbital elements covariance matrix. By propagating a set of Virtual Asteroids to an epoch correspondent to a close approach with a Solar System planet or natural satellite, it is possible to estimate the probability of a partial or total eclipse. The direct observation of an eclipse can provide data useful to improve the asteroid orbit, especially for dim asteroids typically observed only for a few days. We propose two different methods: the first, based on the inclusion of the apparent magnitude residuals into the orbit's least squares minimization process, capable of improving the asteroid's nominal orbit and the related covariance matrix; the second, based on weighting different Virtual Asteroids in relation to their apparent magnitude during the eclipse, useful for recovery purposes. As an application, we have numerically investigated the possibility of a Near Earth Asteroid eclipsed by the Moon or the Earth in the 1990-2050 period. A total of 74 distinct eclipses have been found, involving 59 asteroids. In particular, the asteroid (99942) Apophis has a probability of about 74% to enter the Moon's penumbra cone and a probability of about 6% to enter the umbra cone on April 14, 2029, less than six hours after a very close approach to Earth.