Constraints on Ceres internal strcuture from the Dawn gravity and shape data

Abstract:

Ceres is the largest body in the asteroid belt with a radius of approximately 470 km. It is large enough to attain a shape much closer to hydrostatic equilibrium than major asteroids. Pre-Dawn shape models of Ceres (e.g. Thomas et al., 2005; Carry et al., 2008) revealed that its shape is consistent with a hydrostatic ellipsoid. After the arrival of the Dawn spacecraft in Ceres orbit in March 2015, Framing Camera images were used to construct shape models of Ceres. Meanwhile, radio-tracking data are being used to develop gravity models. We use the Dawn-derived shape and gravity models to constrain Ceres’ internal structure. These data for the first time allow estimation of the degree to which Ceres is hydrostatic. Observed non-hydrostatic effects include a 2.1 km triaxiality (difference between the two equatorial axes) as well as an 660-m center-of-mass – center-of-figure offset. The Dawn gravity data from the Survey orbit shows that Ceres has a central density concentration. Second-degree sectorial gravity coefficients are negatively correlated with topography indicating a peculiar interior structure. We compute the relative crustal thickness based on the observed Bouguer anomaly. Hydrostatic models show that Ceres appears more differentiated based on its gravity than on its shape.

We expand the Ceres shape in spherical harmonics, observing that the power spectrum of topography deviates from the power law at low degrees (Fig. 1). We interpret the decrease of power at low degrees to be due to viscous relaxation. We suggest that relaxation happens on Ceres but, unlike modeled in Bland (2013), it is important only at the lowest degrees that correspond to scales of several hundreds of km. There are only a few features on Ceres of that size and at least one of them (an impact basin provisionally named Kerwan) appears relaxed. The simplest explanation is that Ceres’s outer shell is not pure ice or pure rock but an ice-rock mixture that allows some relaxation at the longest wavelengths. We use the deal.ii finite-element library (Bangerth 2007) to compute relaxed topography spectra. In out future work, we plan to model viscous relaxation to constrain the viscosity profile and thermal evolution.