Computing the Bouguer Correction: Slab Formula Derivation, Hammer Terrain Zones, and Free-Air Gradient in WCSB Gravity Data Reduction
The Bouguer correction is the specific computational step in a gravity data reduction workflow that removes the gravitational attraction of the rock mass between the measurement point and the chosen reference datum — calculated as the attraction of an infinite horizontal slab of uniform density ρ_B and thickness h equal to the station elevation above datum — to isolate the gravity anomaly caused by lateral density variations in the subsurface below the datum plane from the larger and dominant effect of topographic elevation differences between measurement stations. The Bouguer slab formula derives from the gravitational potential of an infinite horizontal slab: BC = 2π G ρ_B h, where G is the gravitational constant (6.674 × 10⁻¹¹ m³/kg/s²), ρ_B is the Bouguer density in kg/m³, and h is the station elevation in metres above the datum. Substituting numerical values for WCSB standard practice (ρ_B = 2,670 kg/m³, G = 6.674 × 10⁻¹¹): BC = 2π × 6.674 × 10⁻¹¹ × 2,670 × h = 0.04193 × h × 10⁻⁵ m/s² = 0.04193 × h µGal per metre of elevation — meaning a 100 m elevated station must have 4.19 mGal subtracted from its Bouguer-corrected gravity to remove the slab effect at the standard WCSB Bouguer density of 2.67 g/cc. The Bouguer correction is applied in conjunction with the free-air correction (FAC), which accounts for the decrease in gravitational acceleration with height above the reference ellipsoid at the rate of 3.086 µGal/cm (or 308.6 µGal/m) without considering the intervening rock mass — the opposite sign and a much larger magnitude than the Bouguer correction for typical WCSB survey elevations. The combined reduction sequence produces the Bouguer anomaly (described under the corresponding entry) by the formula: BA = g_obs + FAC - BC - g_normal(φ), where g_normal is the theoretical normal gravity at the station latitude. In areas of significant topographic relief exceeding 50-100 m within 10 km of a measurement station — typical of the WCSB Rocky Mountain Foothills and northeastern BC terrain — the infinite slab approximation of the Bouguer correction is inadequate, and a terrain correction must supplement it to account for the gravitational attraction of actual topographic masses above and below the slab datum plane, converting the simple Bouguer correction to a complete Bouguer correction and the simple Bouguer anomaly to the complete Bouguer anomaly required for geologically meaningful interpretation in the WCSB Foothills and Front Range gravity surveys.
Key Takeaways
- Free-air correction versus Bouguer correction: opposite signs, different physical content: The free-air correction (FAC = +3.086 µGal/cm of elevation above the geoid) adds gravity back that was lost due to the station's height above the reference level — accounting for the geometric decrease in gravitational acceleration with distance from the Earth's center, without any assumption about what fills the space between the datum and the station. The Bouguer correction (BC = 0.04193 × ρ × h µGal, with ρ in g/cc and h in metres) subtracts the gravitational attraction of the rock that actually fills that space. At 2.67 g/cc Bouguer density: FAC correction rate = 308.6 µGal/m; Bouguer correction rate = 41.9 µGal/m; net correction = FAC - BC = 308.6 - 41.9 = +266.7 µGal/m (the Bouguer plate gradient used in downhole gravity density calculations, approximately matching the theoretical value of 3.084 - 0.419 × 2.67 = 1.965 × 10⁻⁶ µGal/m × 1,000 = ... the combined effect is the standard 308.6 - 41.9 = 266.7 µGal/m positive free-air-minus-Bouguer gradient).
- Hammer terrain zones and the terrain correction computation for WCSB Foothills surveys: The terrain correction removes the gravitational effect of topographic masses that differ from the infinite horizontal slab assumed in the simple Bouguer correction — including the "attraction deficit" of valleys (where missing rock relative to the slab creates a positive terrain correction, because the slab correction over-subtracts for the rock that is not actually present) and the "attraction excess" of hills (where additional rock above the slab creates a negative terrain correction, since the slab correction under-subtracts for this extra mass). The Hammer zone method divides the area around each station into concentric annular compartments (zones A through M, with inner radii from 2 m to 895 m and outer radii from 50 m to 166 km) and computes the terrain correction for each compartment using the average elevation difference between the zone's terrain and the station elevation, weighted by the zone's geometric factor. Modern terrain corrections use digital elevation models (DEM at 10-25 m resolution) and software integration rather than Hammer zone tables, but the physical principle is identical: a valley 50 m deep at 500 m from the station adds +0.5-2 mGal to the terrain correction, and a ridge 200 m high at 2 km distance adds +1-3 mGal, in WCSB Foothills surveys with 500-1,500 m local relief.
- Optimal Bouguer density selection: Nettleton's correlation method: The correct Bouguer density for a survey area minimizes the correlation between the computed Bouguer anomaly and the surface topography — because if the density is wrong, elevated terrain will appear to have either a positive or negative anomaly that is purely an artifact of the elevation, not a real geological contrast. Nettleton's method plots the Bouguer anomaly profile across a topographic feature (a hill or ridge) for several trial Bouguer densities (typically 2.3, 2.4, 2.5, 2.6, 2.67, 2.7, and 2.8 g/cc) and selects the density that makes the anomaly profile flattest (zero correlation with topography) over the feature. In WCSB surveys across the Foothills thrust belt where Devonian carbonate thrust sheets (density 2.68-2.72 g/cc) outcrop, the Nettleton optimal density is 2.68-2.70 g/cc — slightly above the standard 2.67 g/cc used in basin-wide surveys, and the difference (0.01-0.03 g/cc) produces a terrain-correlated artifact of approximately 1-2 mGal across 200 m relief that is significant for detailed structural interpretation but acceptable for regional basin-scale anomaly mapping.
- Bouguer correction in downhole gravity density calculation: the derivation from surface to borehole: The same mathematical framework that defines the surface Bouguer correction is applied inverted to derive formation bulk density from borehole gravity measurements: the difference in gravitational acceleration between two downhole stations (Δg) is related to the bulk density of the formation between them (ρ_bulk) by the same slab formula. Rearranging: ρ_bulk = (FAC_gradient - Δg/Δz) / (2π G) = (3.084 - Δg/Δz) / (4.192 × 10⁻³) g/cc, where Δg is in µGal and Δz is in metres. This is the borehole gravity density equation used by the LaCoste-Romberg downhole gravimeter for WCSB formation density measurements with a 100-500 m radius of investigation, derived directly from the same Bouguer slab formula used to correct surface gravity measurements — demonstrating that borehole gravity is conceptually the inverse of surface gravity reduction, with the formation density as the unknown solved from the measured gravity gradient rather than the density assumed to perform the correction.
- WCSB gravity survey infrastructure: the Canadian Geodetic Survey and AER gravity database: The Bouguer correction requires accurate station elevations above a consistent datum. In Canada, gravity surveys reference their elevations to the Canadian Geodetic Vertical Datum 2013 (CGVD2013), which replaced the older CGVD28 in 2015. The elevation accuracy required for Bouguer-corrected gravity interpretation at 1 µGal precision (the measurement precision of a modern Scintrex CG-6 gravimeter) is approximately 0.3 m — achievable with differential GPS or total station levelling but not with standard handheld GPS (±3-5 m accuracy, introducing ±0.9-1.5 mGal elevation error). The AER maintains a provincial gravity database compiled from historical industry surveys, GSC surveys, and university research surveys, accessible through AER Finder as a public data service. Operators conducting WCSB gravity surveys for salt body delineation or structure mapping typically tie their stations to base stations from the GSC National Gravity Network (spacing approximately 100 km), established with tidal-drift-corrected absolute gravimeter measurements calibrated to the International Absolute Gravity Basestation Network (IAGBN), to ensure consistent datum reference across the survey.
Bouguer Correction Computation: A Single Station Example From the WCSB Foothills
A gravity station in the WCSB Foothills is at elevation 1,285 m above sea level (CGVD2013 datum), latitude 51.42°N. Observed gravity (g_obs): 979,342.8 µGal. Normal gravity at latitude 51.42°N (Somigliana formula): 980,899.2 µGal. Free-air correction: +3.086 µGal/cm × 128,500 cm = +396,554 µGal. Bouguer slab correction (ρ = 2.67 g/cc): -0.04193 × 1,285 = -53.9 mGal = -53,880 µGal. Terrain correction (computed from DEM, zone A through J): +14,230 µGal (large positive value due to 600 m valley 1.5 km to the east and 900 m ridge 3 km to the west). Complete Bouguer anomaly: 979,342.8 + 396,554 - 980,899.2 - 53,880 + 14,230 = -44,652.4 µGal = -44.65 mGal. This strongly negative Bouguer anomaly (typical of mountain ranges where isostatic compensation — crustal root below the mountain — creates a mass deficit) is the regional signal; after removing a second-degree polynomial regional fit, the residual Bouguer anomaly at this station is -2.8 mGal, consistent with a 150 m thick Devonian carbonate structural high beneath the Foothills thrust sheet at approximately 2,200 m depth, as confirmed by two offset wells in the same fold.
Fast Facts
Sigmund Hammer's 1939 paper "Terrain Corrections for Gravimeter Stations" (Geophysics, vol. 4, no. 3) introduced the zone chart method for terrain correction computation that bore his name — the "Hammer zones" — and remained the standard manual method for terrain correction calculation for more than 50 years until digital elevation models and GIS software made automated terrain correction computation practical in the early 1990s. Hammer's zone tables are still published in all standard exploration geophysics textbooks and are used by WCSB gravity survey operators as a check on software-computed terrain corrections at stations with complex near-field topography where DEM resolution may be inadequate to capture the precise local relief.
Related Terms
The Bouguer anomaly that results from applying the Bouguer correction (along with the free-air and terrain corrections) to the observed gravity measurements — and the geological interpretation of the resulting anomaly for salt body detection, basement structure mapping, and carbonate reef identification in WCSB exploration programs — is described under Bouguer anomaly, where the full reduction sequence, density contrast interpretation, and WCSB regional and residual anomaly separation are covered alongside the Devonian salt delineation and basement arch mapping applications that motivate gravity surveys in Alberta and British Columbia. The borehole gravity measurement that applies the identical Bouguer slab mathematics to compute formation bulk density from downhole gravity differences is described under borehole gravity, and the instrument used to make those measurements under borehole gravity meter.