Modeling the interior of Haumea

The Kuiper Belt Object Haumea is one of the most fascinating objects in the solar system. Spectral reflectance observations reveal a surface of almost pure water ice, yet it has a mass of 4.006 × 1021 kg, measured from orbits of its moons, along with an inferred mean radius of 715 km, and these imply a mean density of around 2600 kg m−3. Thus the surface ice must be a veneer over a rocky core. This model is supported by observations of Haumea's light curve, which shows large photometric variations over an anomalously rapid 3.9154-hour rotational period. Haumea's surface composition is uniform, therefore the light curve must be due to a varying area presented to the observer, implying that Haumea has an oblong, ellipsoidal shape. If Haumea's rotation axis is normal to our line of sight, and Haumea reflects with a lunar-like scattering function, then its axis ratios are p = b/a = 0.80 (in the equatorial cross section) and q = c/a = 0.52 (in the polar cross section). In this work, I assume that Haumea is in hydrostatic equilibrium, and I model it as a two-phase ellipsoid with an ice mantle and a rocky core. I model the core assuming it has a given density in the range between 2700–3300 kg m−3 with axis ratios that are free to vary. The metric which my code uses calculates the angle between the gravity vector and the surface normal, then averages this over both the outer surface and the core-mantle boundary. When this fit angle is minimized, it allows an interpretation of the size and shape of the core, as well as the thickness of the ice mantle. Results of my calculations show that Haumea's most likely core density is 2700–2800 kg m−3, with ice thicknesses anywhere from 12–32 km over the poles and as thin as 4–18 km over the equator.