Global Characteristics of Porosity and Density Stratification Within the Lunar Crust from GRAIL Gravity and Lunar Orbiter Laser Altimeter Topography Data

The Gravity Recovery and Interior Laboratory (GRAIL) mission is providing unprecedentedly high-resolution gravity data. The gravity signal in relation to topography decreases from 100 km to 30 km wavelength, equivalent to a uniform crustal density of 2450 kg/cu m that is 100 kg/cu m smaller than the density required at 100 km. To explain such frequency-dependent behavior, we introduce rock compaction models under lithostatic pressure that yield radially stratified porosity (and thus density) and examine the depth extent of porosity. Our modeling and analysis support the assertion that the crustal density must vary from surface to deep crust by up to 500 kg/cu m. We found that the surface density of mega regolith is around 2400 kg/cu m with an initial porosity of 10-20%, and this porosity is eliminated at 10-20 km depth due to lithostatic overburden pressure. Our stratified density models provide improved fits to both GRAIL primary and extended mission data.