The Bouguer Anomaly: Gravity Reduction Workflow, Density Corrections, and Salt Body Detection in Regional Gravity Surveys
The Bouguer anomaly is the residual gravitational acceleration at a measurement point after removing the contributions of: the reference theoretical gravity at that latitude (the normal gravity accounting for the Earth's oblate spheroid shape and rotation); the free-air correction (the decrease in gravity with elevation above the reference ellipsoid, 3.086 µGal/cm of elevation); and the Bouguer slab correction (the gravitational attraction of the rock between the measurement point and the geoid, treated as an infinite horizontal slab of uniform density), with the resulting value representing the gravitational effect of lateral density variations in the subsurface below the datum — the signal of geological interest for petroleum exploration, salt body delineation, basement mapping, and density contrast characterization of sedimentary sequences. Named after French mathematician Pierre Bouguer who conducted gravity measurements in the Andes in the 1730s, the Bouguer anomaly is computed by the reduction sequence BA = g_obs - g_theoretical + FAC - BC, where g_obs is the raw measured gravity; g_theoretical is the International Gravity Formula for the observation latitude (Somigliana formula on the GRS80 reference ellipsoid: g(φ) = 9.7803267715 × (1 + 0.0052790414 sin²φ + ...) m/s²); FAC is the free-air correction (+3.086 µGal/cm of elevation above the geoid); and BC is the Bouguer slab correction (the gravity of an infinite horizontal slab of density ρ_Bouguer and thickness h = elevation above datum: BC = 2π G ρ h = 0.04193 × ρ_Bouguer × h µGal, with ρ in g/cc and h in metres). The standard Bouguer density in WCSB regional gravity surveys is 2.67 g/cc (representing the average density of upper-crustal granitic basement and mixed sedimentary sequence), chosen because it minimizes the correlation between topographic elevation and the Bouguer anomaly for typical continental crustal rocks — a correlation that would otherwise produce spurious anomalies where high-density near-surface rocks expose at topographic highs. In areas of significant topographic relief (WCSB foothills and Rocky Mountain front range, northeastern BC terrain), an additional terrain correction is applied to remove the gravity effect of topographic masses that cannot be approximated as an infinite slab, bringing the fully corrected value to the "complete Bouguer anomaly" that is the standard product for WCSB regional gravity interpretation. The Bouguer anomaly is the same mathematical reduction applied to borehole gravity data, where the Bouguer slab correction between successive measurement stations converts the measured gravity difference to the formation bulk density between those stations — the shared mathematical foundation that makes Bouguer gravity the conceptual link between regional surface surveys and downhole formation density measurement.
Key Takeaways
- Bouguer slab correction: density choice and its effect on the anomaly character: The Bouguer slab correction removes the gravitational effect of the rock between the measurement elevation and the datum plane, but this correction depends on the assumed density of that rock (ρ_Bouguer). If the assumed density is too low, high-elevation measurement points will have their correction under-applied and will show a positive Bouguer anomaly purely from topographic effect (not geological density contrast). If ρ_Bouguer is too high, the anomaly is over-corrected. Nettleton's method for selecting the optimal ρ_Bouguer minimizes the correlation between the Bouguer anomaly and the surface topography over the survey area — the density that makes the Bouguer anomaly profile flattest across topographic relief is the best estimate of the average near-surface density. In WCSB Foothills surveys where carbonate thrust sheets at surface (density 2.65-2.70 g/cc) alternate with valley fills of clastic sediment (density 2.20-2.40 g/cc), a single ρ_Bouguer of 2.67 g/cc may be inadequate and a variable Bouguer density terrain model is required to achieve a geologically meaningful anomaly map.
- Regional versus residual Bouguer anomaly: wavelength separation and WCSB petroleum applications: The Bouguer anomaly map contains superimposed signals from density contrasts at all depths — the deep crustal and lithospheric density variations produce long-wavelength anomalies (hundreds to thousands of km), the basement topography produces medium-wavelength anomalies (tens to hundreds of km), and the sedimentary basin density contrasts (salt bodies, carbonate reefs, fault-bounded density blocks) produce short-wavelength anomalies (1-50 km) that are the targets of petroleum-focused gravity interpretation. Separating the regional Bouguer anomaly (the long-wavelength field from deep crustal structure) from the residual anomaly (the short-wavelength signal from near-surface and shallow crustal structures) is done by polynomial fitting, upward continuation, or wavelength filtering of the gridded Bouguer anomaly data. In WCSB petroleum exploration, the residual Bouguer anomaly after removing a second-degree polynomial regional trend is used to map: Devonian salt dissolution anomalies (low-density collapse features left after halite dissolution); carbonate reef-mound buildups (slightly elevated density relative to surrounding shale); and basement fault geometry (basement relief that affects the thickness of overlying oil-productive formations).
- Evaporite and salt body detection from Bouguer anomaly lows in WCSB Devonian stratigraphy: WCSB Devonian and Mississippian evaporite sequences (Lotsberg Salt at 1,200-1,800 m depth in central Alberta, Cold Lake Evaporite in northeastern Alberta, Beaverhill Lake anhydrite, Muskeg Formation halite) produce Bouguer anomaly lows because their densities (halite: 2.16 g/cc; anhydrite: 2.98 g/cc; mixed evaporite: approximately 2.3-2.5 g/cc average) contrast with the surrounding Devonian carbonate and shale (typically 2.50-2.70 g/cc average). A pure halite body at 1,500 m depth with 200 m thickness and 5 km radius produces a Bouguer anomaly low of approximately 3-6 µGal above its center — small by geological survey standards but measurable with modern gravimeters (precision ±1-2 µGal at 500 m station spacing). Mapping Bouguer anomaly lows in the Alberta Basin has guided the delineation of Lotsberg and Muskeg salt subcrop edges where the salt is thin enough that drilling hazards from salt creep are minimal, and the interpretation of apparent salt dissolution features (collapse breccias and solution-collapse faults) that create structural traps for Devonian oil accumulations in the Rainbow, Zama, and Slave Lake area fields.
- Bouguer anomaly in basement structure mapping for WCSB deep basin exploration: The Precambrian basement beneath the WCSB varies in density from approximately 2.65 g/cc (granite, granodiorite) to 2.85 g/cc (mafic intrusions, high-grade metamorphic rocks) and its depth varies from less than 1 km in the east (Williston Basin shelf) to more than 6 km in the Rocky Mountain Foothills. Basement relief creates Bouguer anomaly gradients because higher-density basement rocks replace lower-density sedimentary section as basement rises. Bouguer anomaly mapping at 2 km station spacing (the standard for WCSB regional gravity surveys) can resolve basement relief features of 500 m amplitude at 5 km lateral scale — adequate for identifying major basement faults and accommodation zones that control the thickness and structural configuration of overlying Devonian and Carboniferous reservoir formations. In the WCSB Deep Basin of Alberta and northeastern BC, basement-controlled structural highs (the Peace River Arch and Athabasca Arch paleotopographic features) are identified in the Bouguer anomaly map as gravity highs corresponding to shallower, denser crystalline basement, and these features guided the early exploration of the Montney and Doig gas formations whose structural configuration follows the basement paleotopographic pattern.
- Terrain correction and the complete Bouguer anomaly in WCSB foothills gravity surveys: The simple Bouguer correction assumes the rock between the measurement elevation and the datum is an infinite horizontal slab — a poor approximation in areas with significant topographic relief where nearby terrain masses (valleys, ridges, peaks) also attract the gravimeter and distort the reading. The terrain correction (TC) removes the gravitational effect of these masses by modeling the actual topography around each station using digital elevation models (DEMs) in concentric zones from the station outward to approximately 200 km. In WCSB Rocky Mountain Foothills surveys where topographic relief can exceed 2,000 m within 10 km of a measurement station, the terrain correction can be 5-50 mGal — large enough that ignoring it would completely distort the interpreted Bouguer anomaly map. The complete Bouguer anomaly (CBA) adds the terrain correction to the simple Bouguer anomaly: CBA = g_obs - g_theoretical + FAC - BC + TC, and is the standard product used for petroleum-relevant geological interpretation in WCSB Foothills and Front Range gravity surveys, where the terrain correction is computed using the ArcGIS or Oasis Montaj software packages with a 10-metre resolution DEM from the Alberta government Spatial Discovery platform.
Bouguer Anomaly Mapping for Lotsberg Salt Delineation in Central Alberta
A regional gravity survey at 1 km station spacing covers a 2,500 km² block of central Alberta over the Lotsberg Salt subcrop edge (halite density 2.16 g/cc, background Cambrian to Devonian carbonate and shale density 2.55 g/cc). After simple Bouguer anomaly reduction at ρ = 2.67 g/cc and polynomial regional trend removal, the residual Bouguer anomaly map shows a 12 mGal low centered on the interpreted Lotsberg Salt depocenter (maximum salt thickness approximately 200 m at depth of 1,400 m), with a gradient of 0.5 mGal/km at the interpreted salt edge. The 12 mGal anomaly amplitude corresponds to a density contrast of approximately 0.41 g/cc (2.57 background minus 2.16 salt) over a 200 m thick body at 1,400 m depth — consistent with the forward model computed using Hammer chart gravity modeling of a disc-shaped halite body with the interpreted dimensions. The anomaly edge is mapped at 1 km precision, locating the 0 mGal residual contour (the interpreted halite pinch-out line) 4.5 km northeast of the position interpreted from the nearest well control (Devonian stratigraphic test that penetrated the base of the Lotsberg section). The well-to-gravity discrepancy of 4.5 km is within the expected structural uncertainty given 10 km well spacing, and the gravity-defined salt edge is used to constrain the regional geological model used for Devonian reef exploration licensing in the adjacent block.
Fast Facts
Pierre Bouguer's 1749 publication La Figure de la Terre, documenting his gravity measurements in the Peruvian Andes in 1735-1744 as part of the French Geodetic Mission, established the free-air correction and the principle of removing the topographic mass contribution from raw gravity readings — the direct conceptual ancestor of the Bouguer correction still applied in every modern gravity survey. The Bouguer anomaly was not named for him during his lifetime; the retrospective naming came from 19th-century geodesists who recognized that Bouguer's systematic gravity reduction workflow was the intellectual foundation of the entire field of applied geodetic and exploration gravity measurement.
Related Terms
The Bouguer correction applied to surface gravity measurements at each station is mathematically identical to the Bouguer slab correction used to convert measured gravity differences between borehole stations into formation bulk density, as described under borehole gravity, which covers the density formula ρbulk = (3.084 - Δg/Δz) / (4.192 × 10⁻³) g/cc and the measurement precision requirements of the downhole gravimeter needed to resolve formation density contrasts of practical interest in WCSB reservoir characterization. The terrain correction that converts the simple Bouguer anomaly to the complete Bouguer anomaly in areas of significant topographic relief is described in the context of the Bouguer correction derivation under Bouguer correction, where the mathematical derivation of the slab gravity formula, the effect of density choice on anomaly character, and the application to both surface and downhole gravity reduction are covered alongside the historical development of gravity reduction standards used in WCSB exploration geophysics.