Magnetic Anomaly Calculator
Our geology & geophysics calculator computes magnetic anomaly accurately. Enter measurements for results with formulas and error analysis.
Calculator
Adjust values & calculateFormula
Where Delta T is the total field anomaly, B_observed is the measured field, B_reference is the IGRF reference, M is magnetization, k is susceptibility, H is field strength, r is body radius, d is depth, and G(i) is the geometric factor dependent on inclination i.
Last reviewed: December 2025
Worked Examples
Example 1: Buried Magnetite Ore Body
Example 2: Sedimentary Basin Anomaly
Background & Theory
The Magnetic Anomaly Calculator applies the following established principles and formulas. Earth science calculators draw on a wide range of measurement scales and physical principles that quantify natural phenomena across geological, atmospheric, and hydrological systems. Earthquake magnitude is most precisely described by the Moment Magnitude Scale (Mw), which replaced the original Richter scale for larger events. Mw is calculated as Mw = (2/3) log10(M0) − 10.7, where M0 is the seismic moment in dyne-centimeters. The Richter scale, while still referenced colloquially, is a local magnitude (ML) measurement derived from peak seismograph amplitude at a standard 100 km distance. Wind intensity is classified using the Beaufort Scale, a 13-point empirical scale (0–12) relating wind speed in knots to observable sea and land effects, with Beaufort 12 corresponding to hurricane-force winds above 64 knots. Tropical cyclone intensity is further categorized by the Saffir-Simpson Hurricane Wind Scale, which assigns Categories 1 through 5 based on sustained wind speed, correlating with expected structural damage. Mineral hardness is quantified on the Mohs scale (1–10), comparing scratch resistance relative to reference minerals from talc (1) to diamond (10). Soil composition analysis measures the proportions of sand, silt, and clay by particle size, alongside organic matter content, bulk density, and porosity, which together determine engineering and agricultural suitability. Seismic wave velocity in rock varies by material: P-waves travel at approximately 5–7 km/s in granite and 1.5 km/s in water, while S-waves travel at roughly 60% of P-wave speeds. Atmospheric pressure decreases with altitude according to the barometric formula: P = P0 × exp(−Mgh / RT), where M is molar mass of air, g is gravitational acceleration, h is altitude, R is the universal gas constant, and T is temperature in Kelvin. Standard sea-level pressure is 101,325 Pa. Tidal calculations use harmonic analysis of gravitational forcing by the Moon and Sun, with the principal lunar semidiurnal tidal constituent (M2) having a period of approximately 12.42 hours.
History
The history behind the Magnetic Anomaly Calculator traces back through the following developments. The systematic study of Earth's structure and processes spans millennia, but the scientific foundations were laid in the seventeenth century. In 1669, Danish naturalist Nicolas Steno published his principles of stratigraphy, establishing the laws of superposition, original horizontality, and lateral continuity — foundational rules for reading rock layers that remain in use today. Scottish geologist James Hutton introduced the concept of uniformitarianism in 1788, proposing that geological processes observable in the present have operated throughout Earth's history at broadly consistent rates. This idea of deep time challenged prevailing biblical chronologies and set the stage for modern geology. Charles Lyell systematized these ideas in his landmark three-volume work Principles of Geology, published beginning in 1830, which directly influenced Charles Darwin's thinking on biological evolution during the voyage of the Beagle. The nineteenth century saw growing curiosity about continental shapes, but a coherent theory awaited Alfred Wegener, a German meteorologist who proposed continental drift in 1912, arguing that the continents had once formed a supercontinent he called Pangaea. His evidence included matching fossil records and geological formations across the Atlantic, but his mechanism was disputed for decades. The theory gained acceptance in the 1960s when seafloor spreading was confirmed through paleomagnetic studies, and plate tectonics emerged as the unifying framework of modern geoscience. The United States Geological Survey was established by Congress in 1879 to classify public lands and examine the geological structure, mineral resources, and products of the national domain. The twentieth century brought instrumental advances, including the global seismograph network deployed after World War II, initially to monitor nuclear tests, which dramatically improved earthquake detection and characterization. Satellite Earth observation began in earnest with the Landsat program launched in 1972, enabling continuous global monitoring of land use, glacier retreat, and vegetation patterns. Today, GPS networks, LIDAR scanning, and ocean-floor mapping provide centimeter-scale precision for tracking tectonic motion, sea level rise, and volcanic deformation in near real time.
Frequently Asked Questions
Formula
Delta T = B_observed - B_reference; M = k * H; Anomaly = (2/3) * mu0 * M * (r/d)^3 * G(i)
Where Delta T is the total field anomaly, B_observed is the measured field, B_reference is the IGRF reference, M is magnetization, k is susceptibility, H is field strength, r is body radius, d is depth, and G(i) is the geometric factor dependent on inclination i.
Worked Examples
Example 1: Buried Magnetite Ore Body
Problem: An observed magnetic field of 49,200 nT is measured over a location where the IGRF reference is 48,000 nT. A spherical ore body with radius 150 m is buried at 600 m depth with susceptibility 0.05 SI. Calculate the anomaly and theoretical response.
Solution: Total anomaly = 49,200 - 48,000 = 1,200 nT\nH = 48,000e-9 / (4pi x 1e-7) = 38,197 A/m\nM = 0.05 x 38,197 = 1,909.9 A/m\nVolume ratio = (150/600)^3 = 0.01563\nAt inclination 60 deg: geometric factor = 2sin^2(60) - cos^2(60) = 1.25
Result: Observed Anomaly: 1,200 nT | Magnetization: 1,909.9 A/m | Half-width: 459 m
Example 2: Sedimentary Basin Anomaly
Problem: A weak negative anomaly of -50 nT is observed. The reference field is 52,000 nT. Estimate parameters for a low-susceptibility body (k=0.001) at 200 m depth with radius 80 m at inclination 45 degrees.
Solution: H = 52,000e-9 / (4pi x 1e-7) = 41,380 A/m\nM = 0.001 x 41,380 = 41.38 A/m\nVolume ratio = (80/200)^3 = 0.064\nGeometric factor at 45 deg = 2sin^2(45) - cos^2(45) = 0.5\nTheoretical anomaly = (2/3)(4pi x 1e-7)(41.38)(0.064)(0.5) x 1e9
Result: Observed Anomaly: -50 nT | Magnetization: 41.38 A/m | Half-width: 153 m
Frequently Asked Questions
What is a magnetic anomaly in geophysics?
A magnetic anomaly is the difference between the observed magnetic field at a location and the expected theoretical or reference field value for that same location. The reference field is typically the International Geomagnetic Reference Field (IGRF), which models the main dipolar field generated by Earth's outer core. Anomalies arise from variations in the magnetization of crustal rocks, which can be caused by differences in mineral composition, particularly the concentration of ferromagnetic minerals like magnetite. Positive anomalies indicate regions where the local field is stronger than expected, often due to highly magnetic rock bodies, while negative anomalies indicate weaker fields. Magnetic anomaly analysis is fundamental to mineral exploration, geological mapping, and understanding tectonic plate history.
How does the Beer-Lambert law relate to magnetic anomaly calculations?
The Beer-Lambert law does not directly relate to magnetic anomaly calculations as it is a principle from optics and spectroscopy. However, both involve attenuation or decay relationships with distance. In magnetic anomaly calculations, the field strength from a buried magnetic body decays with the cube of the distance (for a dipole source), analogous to how light intensity decreases exponentially with path length in Beer-Lambert. The key equation for magnetic anomalies uses the inverse cube relationship where the anomaly amplitude is proportional to the ratio of the body radius cubed to depth cubed. Understanding these decay relationships is essential for estimating the depth and size of subsurface magnetic sources from surface measurements.
What is magnetic susceptibility and why does it matter?
Magnetic susceptibility is a dimensionless quantity that describes how easily a material can be magnetized when placed in an external magnetic field. It is defined as the ratio of induced magnetization to the applied magnetic field strength. In geophysics, magnetic susceptibility varies enormously between rock types: sedimentary rocks typically have very low susceptibility values between 0.0001 and 0.001 SI units, while mafic igneous rocks like basalt can have values of 0.01 to 0.1 SI units, and iron ore deposits can exceed 1.0 SI units. This property is critical for magnetic anomaly interpretation because it directly determines the strength of induced magnetization in crustal rocks and therefore the magnitude of the observed anomaly.
How do you determine the depth of a magnetic source from anomaly data?
Determining the depth of a magnetic source from surface anomaly data involves several analytical and computational techniques. The half-width method uses the horizontal distance at which the anomaly amplitude drops to half its peak value, which for a spherical body equals the depth times a geometric constant approximately equal to 0.766. Euler deconvolution is a more sophisticated approach that uses the spatial derivatives of the anomaly field to solve simultaneously for source position and structural index. The Werner deconvolution method fits line segments to the anomaly profile assuming simple geometric source shapes. Peters half-slope method uses the maximum gradient of the anomaly profile to estimate depth. Each method has strengths and limitations depending on the source geometry and data quality.
What instruments are used to measure magnetic anomalies?
Magnetic anomalies are measured using magnetometers, which come in several types suited to different applications. Proton precession magnetometers measure the total field strength by detecting the precession frequency of hydrogen protons in a fluid, offering accuracy of about 1 nanoTesla. Optically pumped cesium or potassium vapor magnetometers achieve higher sensitivity of 0.01 nanoTesla and faster sampling rates, making them ideal for airborne surveys. Fluxgate magnetometers measure individual components of the magnetic field vector and are commonly used in borehole logging and satellite missions. Superconducting quantum interference devices or SQUIDs provide the highest sensitivity for laboratory measurements. Surveys may be conducted on the ground, from aircraft, aboard ships, or using satellites depending on the spatial scale of the target anomalies.
What causes negative magnetic anomalies?
Negative magnetic anomalies occur when the observed magnetic field is weaker than the regional reference field at a given location. They can result from rock bodies with lower magnetic susceptibility than their surroundings, such as sedimentary basins embedded in more magnetic basement rock. Remanent magnetization oriented opposite to the present-day field can also produce negative anomalies, commonly seen in rocks formed during geomagnetic reversals. Additionally, the geometry of a magnetic body relative to the measurement point and the inclination of the ambient field can create negative lobes flanking a central positive peak. Understanding negative anomalies is important for interpreting geological structure and identifying demagnetized or reversely magnetized formations.
References
Reviewed by Daniel Agrici, Founder & Lead Developer · Editorial policy