Terrain Correction

Key Takeaways

• The terrain correction is computed in zones extending outward from the gravity station, with inner zones requiring higher-precision topographic data than outer zones because the gravitational effect of terrain decreases with distance — directly adjacent terrain (within a few meters of the station) has the greatest per-unit-volume gravitational effect but also the most poorly characterized volume, because the local micro-topography around each measurement point is rarely captured in regional DEM data; the traditional Hammer chart method divided the terrain into concentric zones (labeled A through M in the original Hammer 1939 scheme) and annular sectors, with the field geophysicist estimating the average terrain elevation in each sector from topographic maps and applying tabulated correction values; the inner zones (within about 53 meters) required field estimation by visual inspection because no map captured this level of detail; modern digital terrain correction replaces the manual Hammer chart approach with automated computation from LiDAR-derived DEMs with 1-5 meter resolution, dramatically improving the accuracy of inner zone corrections and making the entire process reproducible; terrain correction error in inner zones remains the dominant accuracy limitation for gravity surveys in steep mountain terrain even with modern DEMs, because the corrections are highly sensitive to small errors in estimated terrain height near the station.
• In offshore gravity surveys, the equivalent of terrain correction is the bathymetric correction that accounts for the varying water depth and seafloor topography — water is less dense than rock (approximately 1.03 g/cc for seawater versus 2.2-2.9 g/cc for common sedimentary rocks), so seafloor topography creates density contrasts that produce gravity anomalies analogous to subaerial terrain effects; the free-air anomaly measured from a ship or satellite includes the combined gravitational effects of seafloor topography, the water-sediment interface density contrast, and the subsurface geological structures being sought; the marine equivalent of Bouguer correction replaces the water column with rock of assumed density to remove the seafloor topography effect, with an iterative terrain correction applied for areas of significant seafloor relief (mid-ocean ridges, continental margins, and seamounts where the seafloor relief can reach thousands of meters); satellite-derived gravity data from altimeters measuring sea surface height (which mirrors seafloor topography through isostatic effects on the ocean surface) provides regional terrain correction capability across the world's ocean basins at spatial resolution of approximately 1-3 km.
• The residual terrain correction for regional exploration gravity interpretation uses a reference surface (typically the mean elevation of the survey area) rather than the exact measurement station elevation, allowing the geophysicist to separate the regional terrain effect from the local geological signal — in a mountainous survey area, the terrain correction referenced to the actual station elevation removes all topographic effects relative to the infinite Bouguer slab; however, the resulting complete Bouguer anomaly still contains the regional isostatic gravity effect of mountain roots and crustal thickness variations, which can be the dominant signal even after terrain correction; a residual terrain correction approach removes only the high-frequency, local terrain roughness effect, preserving the regional isostatic signal for further analysis; the choice between complete Bouguer and residual terrain correction approaches depends on the geological question being asked — basin exploration prefers complete terrain correction to see basin-edge density contrasts, while crustal structure studies may prefer to preserve the regional isostatic signal.
• Terrain correction accuracy directly limits the detection of exploration-scale gravity anomalies in mountainous terrain — a salt dome or deep carbonate reef with a density contrast of 0.2 g/cc against surrounding shale over 1 km thickness at 3 km depth produces a Bouguer gravity anomaly of approximately 5-15 milligals; terrain corrections in the Alps, Andes, or Zagros can reach 50-200 milligals for stations in steep valley or ridge positions, meaning the terrain correction is 5-20 times larger in magnitude than the target signal; a 1% error in terrain correction at a mountain station therefore introduces a 0.5-2 milligal error in the corrected Bouguer anomaly, which is comparable to the geological signal of interest and can create false anomalies or completely mask real ones; this terrain correction sensitivity was a major limitation on the effectiveness of gravity exploration in mountainous regions before digital terrain models and modern computing made precise terrain corrections feasible; in the most extreme mountain terrain, even modern terrain corrections may have residual errors of 1-5 milligals that limit gravity exploration to identifying only the largest and shallowest structural anomalies.
• Borehole gravity surveys require an inner zone terrain correction that accounts for the topography immediately around the wellbore at each measurement depth — borehole gravimeters (small, highly sensitive gravimeters lowered into a wellbore on wireline) measure the vertical gravity gradient between measurement stations in the well, from which density profiles of the formation can be calculated; because the borehole gravimeter samples a large volume of formation (typically 5-15 meters radius versus the few centimeters radius of the conventional density log), it provides a bulk density measurement that is insensitive to borehole washout and is useful for measuring heavy minerals and low-porosity formations where the conventional density log performs poorly; the borehole gravity measurement must be corrected for the terrain at the surface above each measurement point (the overlying topography affects the measured gravity through the vertical derivative of the terrain correction), as well as for the mass of the wellbore fluid and the casing, which represent an artificial density contrast in the formation being measured; the combination of these corrections makes borehole gravity interpretation computationally intensive but provides density data that is unavailable from any surface or standard wireline measurement.