Inferring Shallow Subsurface Density Structure from Surface and Underground Gravity Measurements: Calibrating Models for Relatively Undeformed Volcanic Strata at the Jemez Volcanic Field, New Mexico, USA

Imaging shallow subsurface density structure is an important goal in a variety of applications, from hydrogeology to seismic and volcanic hazard assessment. We assess the effectiveness of surface and subsurface gravity measurements in estimating the density structure of a well-characterized rock volume: the mesa (a small, flat-topped plateau) upon which the town of Los Alamos, New Mexico, USA is located. Our gravity measurements were made on the mesa surface above a horizontal tunnel and underground, within the tunnel. We demonstrate that, in the absence of other geophysical data such as seismic data or muon attenuation, subsurface (tunnel) gravity measurements are critical to accurately recovering geologic structure. Without the tunnel data, our resolution is limited to roughly the surface gravity station spacing, but by including the tunnel data we can resolve structure to a depth of ~ 10 times the surface gravity station spacing. Densities were obtained using both forward modeling and a Bayesian inverse modeling approach, incorporating relevant constraints from geologic observations. We find that Bayesian inversion, with geologically relevant prior, is a superior approach to the forward models in terms of both robustness and efficiency and correctly predicts the orientation and elevation of important geologic features.