Question 6 of 10: Magnetic Anomaly Dependence on Magnetic Latitude and Reduction Methods

Question 6: Magnetic Anomaly Dependence on Magnetic Latitude and Reduction Methods (16.66% of paper)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

A magnetic anomaly's shape (not just its amplitude) depends on the inclination and declination of the ambient inducing field, which in turn depends on magnetic latitude — this is a fundamentally different behaviour from gravity, whose anomaly shape over a given body is the same everywhere on Earth.

Why latitude matters. An induced magnetization is aligned parallel to the Earth's field, whose inclination varies from 0° at the magnetic equator to ±90° at the magnetic poles. A body's dipolar magnetic field, observed at the surface, is projected onto the direction of the ambient field to produce the measured total-field anomaly; the resulting anomaly shape (symmetric vs. skewed, single-peaked vs. dipolar with a low flanking the high) depends directly on this inclination, as well as on declination and the body's own strike relative to magnetic north. At the magnetic equator (inclination ≈0°), a compact body on a north–south profile produces a symmetric anomaly with a central low directly over it flanked by weaker highs to the north and south; at high magnetic latitude (inclination near ±90°), the same body produces a single, roughly symmetric positive total-field anomaly directly over it (positive at both the north and south magnetic poles, because the anomaly is projected onto the field direction); at intermediate latitudes the anomaly is skewed, with the peak displaced toward the magnetic equator and a low on the poleward side.

How this complicates interpretation. Because the same source body produces a different-shaped anomaly at different latitudes (and different strike orientations), a raw total-field anomaly map cannot be interpreted by shape alone using a single, latitude-independent rule of thumb — the interpreter must either know and account for the local inclination/declination explicitly, or transform the data into a latitude-independent form before applying standard interpretation templates (depth rules, Euler deconvolution, forward modelling).

Reduction methods.

Reduction to the pole (RTP). A filter applied in the frequency domain that mathematically recomputes the anomaly as if both the inducing field and the body's induced magnetization were vertical (as at the pole), removing the skew and centring the anomaly directly over the causative body. Strengths: produces anomalies that are simple, symmetric, and directly comparable across surveys at different latitudes; standard, widely available processing step. Weaknesses: becomes numerically unstable near the magnetic equator (the RTP filter is singular as inclination →0°), and assumes purely induced magnetization — it distorts anomalies that carry a significant remanent-magnetization component whose direction differs from the present-day field.

Reduction to the equator (RTE) / pseudo-equator transforms. An alternative low-latitude filter designed specifically to be numerically stable near the magnetic equator (where RTP fails), producing a comparably simplified anomaly shape. Strengths: stable at low latitude, where RTP cannot be used. Weaknesses: less universally implemented than RTP, still assumes purely induced magnetization.

Analytic signal. Rather than removing the latitude dependence directly, this transform (the amplitude of the 3-D gradient of the field) produces maxima over the edges/centre of the causative body that are largely independent of the direction of magnetization (induced or remanent) and of latitude — strictly so for 2-D bodies, approximately for 3-D ones. (The simpler total horizontal gradient also maps edges but, on raw total-field data, still inherits the inclination skew, so it is normally computed on RTP or pseudogravity data.) Strengths: robust to unknown remanence, works at any latitude, good for locating body edges. Weaknesses: gives only body location/edges, not the true anomaly shape or amplitude needed for a quantitative susceptibility/mass estimate.