Time-dependent lunar density models

Simple models of the geochemical and geophysical evolution of the moon were constructed with the particular aim of investigating the mass density distribution within the moon as a function of time. The strongly inverted (densest material on top) density distribution resulting from fractional crystallization of the lunar magma ocean was found to rearrange itself into a highly stable stratification. Density stratification due to compositional variation dominated gravitational stability; thus, thermally driven convection was found to be suppressed in regions with relaxed (rearranged so as to minimize the gravitational potential energy) compositional density variation. The suppression of convection as a heat transfer process results in considerable modification of thermal evolution. The models produced sufficient density segregation from magma oceans of a few hundred km thickness, such that relaxation to the most stable rearrangement yielded a polar moment of inertia approximating that measured for the present-day moon. Thus, lunar models matching moment-of-inertia constraints but lacking a metallic core are motivated.