Home / Science / Earth sciences / Geophysics and seismology / Gravity survey correctionsTool #880

Gravity corrections: drift, latitude (normal gravity), free-air, Bouguer and terrain corrections, with free-air and Bouguer anomalies

Paste a day of gravity meter readings with base station repeats. Get drift-corrected observed gravity, normal gravity from GRS80 or WGS84, free-air and Bouguer slab corrections at your density, ring terrain corrections, and free-air, simple and complete Bouguer anomalies as a table, a profile plot and a CSV.

Keep it one tap away for everyday use.

Bookmark Gravity survey corrections?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Readings
  2. Drift
  3. Corrections
  4. Anomalies
Nothing is uploadedWorks offline once open

Survey readings

Read the base at the start, during and at the end of the day; drift is taken as a straight line between base readings. The example survey is made up.

Anomalies

Download

Terrain correction by rings

One ring per line: inner radius, outer radius (m), then the mean height difference between the station and the ground in each compartment of the ring (m, above or below; the ring is split into as many equal compartments as heights you give). Uses the density above.

Each compartment adds (2πGρ/n)[r2 − r1 + √(r1² + h²) − √(r2² + h²)]. Hills and valleys both lower gravity, so terrain corrections are always added. Put the total in the station's terrain column above. The ring radii and heights here are made up.

What it does

A gravity meter measures tiny differences in the pull of the Earth, a few parts in a hundred million. To see the pull of what is under the ground (an ore body, a buried valley, a salt dome), you must first remove everything else: the meter’s drift, the Earth’s shape and spin, the height of each station and the hills around it. This page does those reductions:

  • Drift. Gravity meters creep during the day. You read a base station several times; the page draws straight lines between those base readings and takes the drift off every station read in between.
  • Observed gravity. Each drift-corrected reading is tied to the base station’s known gravity value. A meter factor turns reading units into mGal if your meter needs one.
  • Latitude correction. Gravity is stronger at the poles. The page computes normal gravity on the reference ellipsoid with Somigliana’s exact formula and the GRS80 or WGS84 constants.
  • Free-air correction. Gravity falls by about 0.3086 mGal for every metre you climb. You can also use the second-order formula, which allows for latitude and height squared.
  • Bouguer correction. The rock between the station and the datum pulls too. It is taken as a flat slab: 2πGρh, about 0.1119 mGal per metre at 2670 kg/m³.
  • Terrain correction. Hills above and valleys below the station both reduce gravity, so a terrain correction is always added. The ring calculator works it out from the mean height of each compartment of rings around the station.
  • Anomalies. Free-air anomaly = observed − normal + free-air correction. Simple Bouguer anomaly = free-air anomaly − Bouguer correction. Complete Bouguer anomaly = simple Bouguer anomaly + terrain correction.

How to use it

  1. Enter the base station’s name as it appears in your table, and its gravity value in mGal.
  2. Enter the Bouguer density. 2670 kg/m³ is the usual convention for crustal rock; use your own if you know the local rock.
  3. Paste or open the readings in time order: station, time (hh:mm or decimal hours), reading, latitude, elevation in metres, and optionally the distance along the line and a terrain correction in mGal.
  4. Read the anomalies in the table and the profile plot, and download the CSV.
  5. For terrain, type one ring per line in the ring calculator, read off the total and put it in the station’s terrain column.

Worked example

The example is a made-up line of 11 stations, 200 m apart, over a buried dense body, read between 08:00 and 11:30 with the base read three times. The base readings rise from 2512.300 to 2512.440, so the drift is +0.040 mGal/h. At station G04 (800 m along the line, 152.8 m high), the free-air anomaly is 4.519 mGal, the simple Bouguer anomaly −12.590 mGal and, with its 0.21 mGal terrain correction, the complete Bouguer anomaly −12.380 mGal. Along the line the complete Bouguer anomaly runs from −13.967 to −10.945 mGal, a regional slope with a bump of about 0.8 mGal over the body at 1100 m. The example ring set gives a terrain correction of 0.949 mGal.

Good to know

  • Elevations matter most. An error of 1 m in height is about 0.2 mGal in the Bouguer anomaly, so survey heights carefully.
  • The density you pick changes the Bouguer anomaly wherever the ground is high. If the anomaly mirrors the topography, the density is probably wrong (the Nettleton method).
  • Earth tides move gravity by up to about 0.3 mGal over a day. Frequent base readings take most of the tide out with the drift.
  • For long lines or mountains, surveys also add a curvature (Bullard B) correction to the slab; it is not included here.

FAQ

What is a mGal? A milligal is 10⁻⁵ m/s², about a millionth of the Earth’s gravity.

Should elevations be above sea level or the ellipsoid? Modern practice uses heights above the ellipsoid (from GNSS) with normal gravity on the ellipsoid. Older surveys used sea level. Be consistent within a survey.

Is my data sent anywhere? No. All the sums are done on this page and nothing leaves your device.

Where the numbers come from

  • Normal gravity: Somigliana's closed formula with the GRS80 constants (Moritz 1980, Geodetic Reference System 1980) and the WGS84 constants (NIMA TR8350.2); general geodetic constants.
  • Free-air gradient 0.3086 mGal/m, the second-order free-air formula (as given by Hinze et al. 2005, Geophysics 70 J25), the Bouguer slab 2πGρh and the ring-compartment formula of Hammer (1939) are standard results, written by fibub from the formulas; no zone charts or tables are reproduced.
  • Checked against an independent NumPy reduction of the same survey and SciPy numerical integration of each ring sector.
  • The example survey, base gravity value, ring radii and heights are made up by fibub.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-09.

Related tools

Traverse calculator #132 Paste or open a survey traverse as bearings and distances or as angles and distances. Get the angular and linear misclosure, the precision ratio, Bowditch or transit adjusted coordinates, the area inside a loop, and the bearing and distance between points and where two lines or two distances meet. Level book reducer #133 Paste or open a levelling field book (back, intermediate and fore sights). Get reduced levels by both rise-and-fall and height of collimation, every arithmetic check, the closing error against a benchmark or loop, the allowed misclosure from your own k, adjusted levels and a profile. Publication plot maker #27 Paste columns or open a CSV or Excel file and get a clean scientific figure with error bars, sized for a journal column. Download vector SVG or PDF, or PNG at 300 or 600 dpi. Your data stays in your browser. Earthquake locator and magnitude calculator #878 Paste P and S arrival times from three or more seismic stations. Get each station's S−P distance, the epicentre and origin time by least squares (Geiger's method, on a flat map or the globe), optional depth, error estimates, a map with distance circles, local magnitude ML from amplitudes, and Mw from seismic moment or fault size. Magnetic anomaly modeller #881 Model the total-field magnetic anomaly of a buried sphere, horizontal cylinder, thin dipping dyke or vertical-sided prism in any field inclination and declination, with optional remanent magnetisation. See the profile and a map, and compare with pasted survey readings by RMS misfit. Seismic refraction interpreter #882 Paste first-arrival picks from a forward shot and, if you have one, a reverse shot. The page splits them into 2 or 3 layers, fits each line, and gives speeds, depths by the intercept-time and crossover methods, a dipping-layer solution from the two shots, reciprocal time check and a forward model to compare. Resistivity sounding (VES) inversion #883 Paste Schlumberger or Wenner sounding readings (AB/2, MN/2, ΔV and I, or apparent resistivity). Get apparent resistivities, a layered-earth model of 2 to 6 layers by damped least-squares inversion, the log-log sounding curve, a resistivity-depth plot, error factors and a CSV.