Determination of Celestial Body Principal Axes via Gravity Field Estimation

The determination of physical parameters of celestial bodies provides critical navigational and scientific information. Determining the mass, rotation state, and density distribution is an important task of the navigation team on an operational mission. One representation of the internal density distribution, spherical harmonics, may be leveraged to compute the orientation of the principal axis frame of the celestial body, which in turn informs the dynamics of the body’s motion. Using data from NASA’s Dawn mission to the massive asteroids (1) Ceres and (4) Vesta, a method is proposed that refines the knowledge of the principal axis frame by utilizing estimated spherical harmonic coefficients. Applying an iterative rotation scheme to the solved-for body-fixed frame leads to convergence on a frame that is very close to the actual dynamical principal axis frame.