3D geologic modeling and mapping often relies on gravity modeling to identify key geologic structures, such as basin depth, fault offset, or fault dip. Such gravity models generally assume either homogeneous or spatially uncorrelated densities within modeled rock bodies and overlying sediments, with average densities typically derived from surface and drill-hole sampling. The noise contributed to the gravity anomaly by these density assumptions is zero in the homogeneous case and typically <200 μGal in the uncorrelated case. Rock bodies and sediments, however, show both a range of densities and spatial correlation of these densities, in both surface and drill-hole samples, and this correlation causes an increase in power in the low frequency content of the resulting gravity anomaly. Spatial correlation of densities can be modeled as a Gaussian random field (GRF), with the random field parameters derived from drill-hole and geologic map data. Data from alluvial fan sediments in southern Nevada indicate correlation lengths of up to 300 m in the vertical dimension and kilometers in the horizontal dimension. The resulting GRF density models show that the noise contributed to the measured gravity anomaly is of low frequency and can be several mGal in amplitude, contradicting the common attribution of lower frequencies to deeper sources. This low-frequency noise increases in power with an increase in sediment thickness. Its presence increases the ambiguity of interpretations of subsurface geologic body shape, such as basin analyses that attempt to quantify concealed basement fault depths, offsets, and dip angles. In the southwestern United States, where basin analyses are important for natural resource applications, such ambiguity increases the uncertainty of subsequent process modeling.