04-Geol-B10 · December 2017
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
EGBC National Exam — Geological Engineering, 04-Geol-B10-1 Gravity and Magnetic Fields, 2017-Dec. Closed book; no calculator permitted. All ten questions require an answer in essay format, with diagrams used wherever appropriate. The exam instructs "choose six (6) of the following ten (10) questions, the first six as they appear in the answer book will be marked, each of equal value".
Reference texts: Telford, Geldart & Sheriff, Applied Geophysics, 2nd ed. (physical properties ch.2 & 5; gravimeters, gravity reduction and terrain correction ch.2; magnetometers and magnetic surveying ch.4–5; anomaly interpretation throughout); Kearey, Brooks & Hill, An Introduction to Geophysical Exploration, 3rd ed. (survey design, diurnal correction, case-history applications ch.6 & 7); Blakely, Potential Theory in Gravity and Magnetic Applications (potential-field theory, Fourier-domain filters, reduction-to-pole, non-uniqueness ch.2, 5, 9 & 12).
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.
At the magnetic pole, the Earth's field is purely vertical (inclination I = 90°), and a simple, symmetric, isolated magnetic source there produces a symmetric, single-peak (monopole-like) anomaly directly over the body. At any other latitude, the field is oblique (inclined at some angle I < 90° and with a declination D), so the SAME source produces an anomaly that is asymmetric, skewed toward the pole-facing side of the true source location, and often dipolar (a positive lobe paired with a negative lobe), which makes anomalies difficult to correlate directly with mapped geology or with each other by eye. RTP is a mathematical filter applied to remove this obliquity, transforming the anomaly into what it would look like if it had been measured at the pole, directly over its source, so that anomaly patterns can be interpreted geologically and compared consistently across a survey.
RTP is performed in the wavenumber (Fourier) domain. The gridded magnetic anomaly is Fourier transformed, multiplied by an RTP operator that depends on the field's inclination I, declination D, and the wavenumber direction θ, of the general form
$$\Phi_{RTP}(\theta) = \frac{1}{\big[\sin I + i\cos I \cos(D-\theta)\big]^{2}}$$
and the result is inverse Fourier transformed back to the spatial domain, giving the reduced-to-pole grid. Multiplying by this operator effectively re-orients both the inducing field AND the induced magnetization to vertical, removing the skew and asymmetry that oblique magnetization otherwise imposes.
The RTP operator's denominator contains sin I, which approaches zero as the survey's magnetic latitude approaches the magnetic equator (I→0); the operator's amplitude then blows up, so RTP becomes numerically unstable and amplifies noise catastrophically at low magnetic latitudes — a well-known practical limit on RTP. RTP also assumes the source's magnetization is purely INDUCED, aligned with the present inducing field; if the source carries significant remanent magnetization (a high Koenigsberger ratio Q, common in mafic/volcanic rocks), the true magnetization direction differs from the inducing field direction and standard RTP, based on the field direction alone, will not fully symmetrize the anomaly, potentially producing an artificially distorted result.
At low magnetic latitude, reduction to the equator (RTE) or pseudo-reduction schemes use an amplitude-based (rather than direct phase-inversion) formulation that avoids the sin I singularity, trading some resolution for numerical stability. Where remanence is suspected or unknown, the analytic signal (the 3D total-gradient amplitude, combining the vertical and both horizontal derivatives) is a transform that is mathematically INDEPENDENT of the magnetization direction altogether — it locates source edges correctly whether the source is purely induced, purely remanent, or a mixture, avoiding the RTP assumption entirely, at the cost of a somewhat broader, less sharply peaked response than a correctly-executed RTP.
