Home / Science / Earth sciences / Geophysics and seismology / Magnetic anomaly modellerTool #881

Magnetic anomaly modelling: total-field anomaly of a sphere, horizontal cylinder, thin dyke or prism, with remanence, profile, map and RMS misfit

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.

Keep it one tap away for everyday use.

Bookmark Magnetic anomaly modeller?

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

  1. Field
  2. Body
  3. Profile
  4. Misfit
Nothing is uploadedWorks offline once open

Earth's field and profile

The field values are an example; look up the main field for your site and date from a geomagnetic model. Distances run along the profile; the map's other axis runs 90° clockwise from it.

Body

Total-field anomaly

Map of ΔT around the body, red above zero and blue below. The dashed line is the profile.

Download

Compare with observed data

Change the body until the curves agree. The second RMS figure ignores a constant difference, such as a regional level. The example readings are made up.

What it does

Magnetised rocks, such as magnetite-rich dykes, ore bodies, basalt flows or buried steel, add a small field to the Earth’s own. A magnetometer measures the total-field anomaly ΔT: how much stronger or weaker the field is than it would be without the body. Because the Earth’s field is tilted almost everywhere, the anomaly of even a simple body has a high on one side and a low on the other. This page computes that shape for four classic bodies so you can understand a survey line or check an interpretation:

  • Sphere: a compact body such as a boulder of ore, a buried drum or a pipe end, behaving as a single dipole.
  • Horizontal cylinder: a long body crossing the profile at right angles, such as a pipeline, a buried channel or a fold hinge.
  • Thin dyke or sheet: a tabular body that strikes across the profile, with any dip and any down-dip extent.
  • Vertical prism: a box with vertical sides, finite in all three directions, for plugs and ore lenses.

The body is magnetised by induction (susceptibility χ times the Earth’s field) and, if you tick the box, by remanence, a magnetisation frozen into the rock with its own direction. The page shows the Koenigsberger ratio Q, remanent over induced.

How to use it

  1. Enter the Earth’s field: strength F in nT, inclination I (positive downward, as in the northern hemisphere) and declination D. Enter the direction the profile runs, the sensor height and the profile’s start and end.
  2. Choose the body, its susceptibility in SI units, its position along the profile and its size and depth.
  3. Read the high, the low and the peak-to-trough size, look at the profile and the map, and download the profile as CSV.
  4. To interpret a survey line, paste or open the readings (distance and ΔT). Adjust the body until the RMS misfit is as small as the noise. The second RMS figure ignores a constant offset, because surveys often sit on a regional level.

Worked example

The page opens with a field of 50,000 nT at inclination 65° and declination 5°, a north-running profile from −300 to 300 m and a sphere of radius 30 m centred 80 m deep at 0 m, with χ = 0.05 SI. The induced magnetisation is 1.99 A/m. The anomaly has a high of +73.2 nT at −18 m, south of the body because the field dips north, and a low of −9.7 nT at 80 m, with a peak to trough of 82.9 nT.

The example readings were made from a sphere 15 m further north, 85 m deep and of radius 32 m, then shifted by 2 nT and given some noise. Against the starting model the RMS misfit is 7.45 nT (6.94 nT after a level shift). Move the sphere to 15 m, 85 m deep and radius 32 m, and the misfit drops to 2.09 nT, or 0.69 nT after the level shift: about the size of the noise.

Good to know

  • The same anomaly can come from different bodies. A deeper, stronger body can look like a shallower, weaker one; use geology and other data to choose.
  • The width of the anomaly says most about depth. As a rough rule, the half-width of a sphere’s peak is about its depth.
  • 2D bodies (cylinder and dyke) are infinitely long across the profile. Use the prism when the body is short.
  • Strong remanence can flip an anomaly upside down. Reversed remanence, as in many old basalts, often shows as a low.
  • ΔT here is the anomalous field projected on the main field direction, which is what total-field magnetometers record when the anomaly is small compared with F.

FAQ

What susceptibility should I use? It varies a lot: sedimentary rocks are often below 0.001 SI, basalts around 0.01 to 0.1 and magnetite ores much higher. Measure samples if you can.

Why is the high not over the body? Except at the magnetic poles, both the magnetisation and the measured component are tilted, so the high moves toward the equator side and a low appears on the pole side.

Is my data sent anywhere? No. Everything is worked out on this page and nothing leaves your device.

Where the numbers come from

  • Dipole and 2D line-dipole fields, the thin-sheet integral and the face-by-face (surface pole) solution for a uniformly magnetised prism are standard potential-field results (see Blakely 1995, Potential Theory in Gravity and Magnetic Applications), written by fibub from the physics.
  • Checked against independent NumPy and SciPy models: a dipole potential differentiated numerically, and point dipoles integrated along strike, over a sheet and through a box.
  • The field values, body and example readings are made up by fibub. Look up the real main field for a site from a geomagnetic reference model.

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

Stereonet plotter #33 Plot planes and lines on an equal-area (Schmidt) or equal-angle (Wulff) stereonet. Get great circles, poles, the mean orientation with its 95 % cone and a rose diagram, and download SVG or PNG. 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. GPR radargram viewer #513 Open ground-penetrating radar profiles in SEG-Y (rev 0 or 1, IBM or IEEE floats, integers, either byte order) or GSSI DZT format in your browser. Set time zero, dewow, remove the background, add linear, exponential or AGC gain, show depth from a wave speed, fit a hyperbola for the speed, view single traces and save the radargram as PNG. 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. Gravity survey corrections #880 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. Seismogram viewer #879 Open SAC and miniSEED seismogram files (Steim-1, Steim-2, integers and floats) in your browser. See several traces at once with zoom and pan, remove the trend, taper, apply a Butterworth bandpass (one pass or zero phase), see the amplitude spectrum and a spectrogram, tap to pick P and S, and download CSV.