Morphology driven density distribution estimation for small bodies

Abstract We explore methods to detect and characterize the internal mass distribution of small bodies using the gravity field and shape of the body as data, both of which are determined from orbit determination process. The discrepancies in the spherical harmonic coefficients are compared between the measured gravity field and the gravity field generated by homogeneous density assumption. The discrepancies are shown for six different heterogeneous density distribution models and two small bodies, namely 1999 KW4 and Castalia. Using these differences, a constraint is enforced on the internal density distribution of an asteroid, creating an archive of characteristics associated with the same-degree spherical harmonic coefficients. Following the initial characterization of the heterogeneous density distribution models, a generalized density estimation method to recover the hypothetical (i.e., nominal) density distribution of the body is considered. We propose this method as the block density estimation, which dissects the entire body into small slivers and blocks, each homogeneous within itself, to estimate their density values. Significant similarities are observed between the block model and mass concentrations. However, the block model does not suffer errors from shape mismodeling, and the number of blocks can be controlled with ease to yield a unique solution to the density distribution. The results show that the block density estimation approximates the given gravity field well, yielding higher accuracy as the resolution of the density map is increased. The estimated density distribution also computes the surface potential and acceleration within 10 % for the particular cases tested in the simulations, the accuracy that is not achievable with the conventional spherical harmonic gravity field. The block density estimation can be a useful tool for recovering the internal density distribution of small bodies for scientific reasons and for mapping out the gravity field environment in close proximity to small body’s surface for accurate trajectory/safe navigation purposes to be used for future missions.

[1]  Robert A. Werner Spherical harmonic coefficients for the potential of a constant-density polyhedron , 1997 .

[2]  Daniel J. Scheeres,et al.  Radar Imaging of Binary Near-Earth Asteroid (66391) 1999 KW4 , 2006, Science.

[3]  P. Kuchynka,et al.  Constraining Ceres’ interior from its rotational motion , 2011, 1107.2051.

[4]  Zuber,et al.  The shape of 433 eros from the NEAR-shoemaker laser rangefinder , 2000, Science.

[5]  Daniel J. Scheeres,et al.  Surface Gravity Fields for Asteroids and Comets , 2013 .

[6]  T. McCord,et al.  Ceres’ evolution and present state constrained by shape data , 2010 .

[7]  S. Ostro,et al.  Dynamical Configuration of Binary Near-Earth Asteroid (66391) 1999 KW4 , 2006, Science.

[8]  Nitin Arora,et al.  Global Point Mascon Models for Simple, Accurate, and Parallel Geopotential Computation , 2012 .

[9]  Yu Takahashi,et al.  Gravity Field Characterization around Small Bodies , 2013 .

[10]  R. S. Hudson,et al.  Shape of Asteroid 4769 Castalia (1989 PB) from Inversion of Radar Images , 1994, Science.

[11]  W. M. Kaula Theory of satellite geodesy , 1966 .

[12]  Estimating asteroid density distributions from shape and gravity information , 2000 .

[13]  B. Tapley,et al.  Statistical Orbit Determination , 2004 .

[14]  B. Jones Efficient models for the evaluation and estimation of the gravity field , 2010 .

[15]  T N Titus,et al.  Dawn at Vesta: Testing the Protoplanetary Paradigm , 2012, Science.

[16]  R. Park,et al.  Estimating Small-Body Gravity Field from Shape Model and Navigation Data , 2008 .

[17]  D. Scheeres,et al.  Exterior gravitation of a polyhedron derived and compared with harmonic and mascon gravitation representations of asteroid 4769 Castalia , 1996 .

[18]  Radar Image of Asteroid 1989 PB , 1990, Science.