Reading the Earth's Density — Without Drilling a Single Foot
Microgravity tomography is one of geophysics’ most underappreciated tools. At GSR, it’s revealing a striking density contrast between deep competent carbonates and the formations above them.
Every rock formation has a density. And because density generates a gravitational field, variations in the density of subsurface rock produce measurable, if subtle, variations in the gravitational acceleration measured at the surface. Gravity anomalies (the spatial pattern of those variations) are one of the oldest tools in applied geophysics, used for more than a century to infer the lithological composition of rock bodies that no one has yet drilled. The technique is non-invasive, covers large areas quickly, and responds to physical properties that are directly relevant to reservoir quality: porosity, fluid content, lithological type, and the presence of structural voids.
Microgravity tomography, the modern evolution of gravity surveying, has substantially refined both the precision and the interpretive power of the method. Where traditional gravity surveys measured regional-scale density contrasts in units of milligals, microgravity instruments operating at the microgal level can resolve density variations across intervals of tens of meters rather than hundreds — approaching the resolution needed for reservoir-scale interpretation. Combined with inversion algorithms that convert surface gravity measurements into three-dimensional subsurface density models, microgravity tomography has become a legitimate pre-drill characterization tool for geothermal exploration, mining, and groundwater assessment. At the GSR Geothermal Energy Project site, the technique is revealing a density architecture that has direct implications for how the reservoir target should be understood and how the drilling program should be designed.
2.71
Approximate bulk density (g/cm³) of the Mississippian Leadville-Ouray Carbonate — significantly higher than the overlying clastic and evaporite formations it underlies
~1 µGal
Measurement precision achievable with modern CG-6 and similar superconducting gravimeters — sufficient to resolve meter-scale density contrasts at reservoir depths
3D
Dimensional output of constrained gravity inversion — a volumetric density model that guides well placement before the first drillbit turns
The physics of gravity surveying
Why the Earth’s pull varies and what that variation encodes
Gravity at the Earth’s surface is not uniform. It varies with latitude, elevation, topographic relief, tidal forcing from the Moon and Sun, and, most relevant to exploration geophysics, the density distribution of the rock beneath the measurement point. Newton’s law of gravitation is the governing principle: a mass concentration at depth exerts a slightly greater gravitational pull on the surface instrument above it than the surrounding lower-density material does. That excess pull is the gravity anomaly — the signal that geophysicists process, correct for confounding factors, and invert to infer subsurface density structure.
The corrections applied to raw gravity data before interpretation are extensive, and their rigor largely determines the quality of the final model. Latitude correction accounts for the Earth’s ellipsoidal shape and rotation. The free-air correction removes the effect of measurement elevation above a reference datum. The Bouguer correction subtracts the gravitational attraction of the rock mass between the measurement point and sea level. Terrain corrections account for the irregular topography surrounding each measurement station. This is a particularly important step in rugged terrain like the western Colorado Plateau, where topographic relief can produce gravity signatures that compete in magnitude with the subsurface lithological anomalies being sought. After these corrections, the residual gravity field (the Bouguer anomaly) reflects subsurface density variations with the regional geological trend mathematically separated from local structural and lithological features of interest.
The density inversion problem
From surface measurements to subsurface models
Gravity inversion — the mathematical process of converting a map of surface gravity anomalies into a three-dimensional model of subsurface density — is mathematically non-unique. An infinite number of subsurface density distributions can, in principle, produce the same surface gravity field. The practical resolution of this ambiguity requires incorporating geological constraints: known formation depths from any existing well control, seismic reflection or refraction data that defines structural geometry, and lithological density ranges from laboratory measurements on outcrop or core samples of the formations in question. When these constraints are applied to bound the inversion, the resulting density model is not a pure mathematical artifact but a geologically informed interpretation that reflects both the gravity data and the structural framework derived from independent observations.
At the GSR site, the inversion benefits from an unusually favorable set of external constraints. Regional well data from the western Colorado Plateau provides formation tops and thickness estimates for the major stratigraphic units in the section. The structural framework derived from surface mapping of the Ridgway and Johnson Spring fault systems defines the geometric boundaries within which the inversion is performed. Laboratory density measurements on outcrop samples of the Mississippian Leadville-Ouray Carbonate — the primary reservoir target — establish the density endmember for the inversion’s deepest resolvable horizon. Together, these constraints transform the inversion from an underdetermined mathematical problem into a geologically testable hypothesis about where dense, competent carbonate rock is present at depth and where it is disrupted by porosity, fluid, or structural modification.
“Constrained gravity inversion, when applied to a well-characterized stratigraphic section, can resolve lateral density variations of 0.05–0.10 g/cm³ at reservoir depths — sufficient to distinguish between tight, low-porosity carbonate and fractured or vuggy carbonate with fluid-filled secondary porosity.”
— Nabighian, M.N. et al., Geophysics, Society of Exploration Geophysicists, 2005
What the GSR gravity signature shows
A density contrast that points directly to the reservoir target
The near-surface gravity signature at the GSR project area exhibits a pattern that, when processed and inverted against the known stratigraphic framework, points to a striking lithological contrast at depth. The overlying formations — Mesozoic and Paleozoic clastic sequences, evaporite-bearing intervals, and relatively low-density sedimentary units that have been structurally cased off in the drilling program design — produce a comparatively muted gravity response. The deeper Mississippian Leadville-Ouray Carbonate sequence, by contrast, generates an anomalous positive gravity signature consistent with its substantially higher bulk density relative to the overlying section.
That density contrast is the microgravity tomography equivalent of a directional signal: the high-density carbonate at depth is physically distinct from the formations above it in ways that the gravity field is actively encoding. Where the anomaly is strongest, the carbonate sequence is most intact, most laterally continuous, and most likely to be competent and well-developed as a reservoir unit. Where the anomaly is interrupted or attenuated, the carbonate may be thinned by structural erosion, displaced by faulting, or modified by fluid-related diagenesis — all factors that the drilling program should anticipate and design around.
Critically, the gravity anomaly is not spatially uniform across the site. It exhibits lateral variations — local maxima and minima within the broader positive anomaly associated with the Leadville-Ouray sequence — that correspond, in the constrained inversion model, to variations in carbonate thickness, depth, and potentially porosity structure. Those lateral variations carry exploration significance: they point toward drilling locations that are geometrically optimal for intersecting the thickest, most continuous, and structurally least-disrupted portion of the reservoir target.
Microgravity and reservoir porosity
How density encodes pore structure — and why that matters for flow
The relationship between density and porosity is one of the most direct in petrophysics. Porosity — the fraction of rock volume occupied by void space, whether primary depositional pores or secondary fractures and vugs — reduces bulk density. A carbonate with 10% porosity is measurably less dense than the same rock at 2% porosity, and that difference propagates upward as a measurable gravity signal when the porous interval is thick enough and the density contrast is sharp enough. Microgravity surveys calibrated against laboratory density measurements on formation samples can, in favorable circumstances, distinguish between tight matrix carbonate and fractured or vuggy carbonate with fluid-filled secondary pore space — the very distinction that separates a reservoir unit from a non-reservoir unit in a geothermal context.
This porosity sensitivity is particularly valuable at the GSR site because the Leadville-Ouray Carbonate is expected to exhibit heterogeneous porosity — a mix of tight matrix intervals, fracture-enhanced zones, and potentially karstified horizons that vary spatially in ways that surface outcrop mapping cannot resolve at drilling scale. The microgravity inversion, where anomalies can be tied to the structural framework from Article 1 of this series, provides a pre-drill approximation of that porosity heterogeneity — not a high-resolution core description, but a three-dimensional density model that informs where within the relay accommodation zone the reservoir is likely to be most productive.
Integration with the structural framework
Gravity and structural geology read the same story from different directions
The microgravity dataset at GSR does not stand in isolation. Its power as an exploration tool derives from its integration with the structural interpretation described in the previous article in this series: the extensional relay dynamics of the Ridgway and Johnson Spring fault systems that define the accommodation zone as the primary permeability target. Where the gravity inversion identifies the thickest and most continuous expression of the high-density Leadville-Ouray Carbonate, and where that interval coincides spatially with the structurally predicted zone of maximum fracture density within the relay ramp, the convergence of two independent geophysical lines of evidence defines a drilling target with a higher confidence basis than either dataset alone would support.
This multi-method convergence is standard practice in mature hydrocarbon exploration and is increasingly applied in geothermal resource characterization as the industry adopts the subsurface analysis toolkit developed by the petroleum sector over decades. The principle is straightforward: independent geophysical observations that point to the same subsurface feature provide mutual corroboration that reduces the uncertainty in the pre-drill resource model. They do not eliminate subsurface risk — no surface geophysical method can substitute for the direct measurement that drilling provides — but they compress the range of likely outcomes and improve the probability that the first well intersects the target reservoir in a productive configuration.
Our perspective
Geophysics as a pre-drill risk management tool
At Mycelian & Co., the microgravity tomography work at the GSR site is part of a pre-drill characterization program designed to maximize the geological information available before capital is committed to a well program. The density inversion, integrated with the structural framework from fault mapping and the formation depth control from regional well data, provides a three-dimensional subsurface model that informs both well location selection and the expected formation sequence the drillbit will encounter. It is not a substitute for drilling — it is the analytical foundation on which a rigorous drilling program is built.
The anomalous near-surface gravity signatures at GSR, driven by the density contrast between the deep Leadville-Ouray Carbonate and the overlying cased-off formations, are telling a story about what lies beneath. Reading that story correctly — before the first well is permitted, before the first dollar of drilling capital is committed — is what separates exploration conducted on a geological basis from exploration conducted on hope alone.
By Manny Melendez · Founder & Principal, Mycelian & Co. · May 21, 2026
Interested in knowing how you can be a part of the GSR Geothermal Energy Project?
Our Team is selectively available for conversations with institutional partners, family offices, and project sponsors aligned with this space. Reach out here or connect with us at Mycelian & Co.