Electrical, electromagnetic and ground-penetrating methods

Fictional unit investigation: A hypothetical valley aquifer is imaged by resistivity electrodes, frequency-domain induction and radar lines. Clay, salinity and moisture have overlapping responses; students must avoid declaring contamination from conductivity alone.

Electrical resistivity is a bulk property influenced by interconnected water, ions, temperature, pore geometry, clay surface conduction and mineral phases. Four-electrode apparent resistivity is not the resistivity at a single sampled depth. Electromagnetic surveys sense conductivity through inductive coupling, and radar travel time depends on effective electromagnetic velocity and attenuation. Different array configurations sample different sensitivity regions, so apparent contrasts require geometry-aware interpretation. Saltwater and clay can both yield high conductivity; dry rock, air voids and clean sand may give resistive signatures. Instrument coupling, buried utilities, cultural noise and inversion regularization add further uncertainty.

For each lesson the forward physics and observation type are explicit. The lesson is complete only when the student can name a competing explanation, check units in the worked calculation and explain a measurement that could reject the preferred interpretation. Every stated numerical value is invented for this textbook; none is a field station or an operational parameter.

The physical mechanism and its measurable consequence

The measured potential difference generated by injected current is converted using an electrode geometry factor to apparent resistivity. The result is an effective measure of a sensitivity volume, not a local rock sample. Contact resistance and near-surface heterogeneity affect data quality.

Acquisition protocol and observable variables

Record electrode coordinates, current, voltage, polarity reversal, stacking error, contact resistance and reciprocal measurements.

The inverse ambiguity: alternative explanations

Clay conduction, saline water or metallic mineralization can each lower apparent resistivity; high values may indicate dry soils or lithology changes.

Worked calculation: synthetic classroom data

Question — Wenner apparent resistivity. Wenner spacing 5 m, V=0.12 V, I=0.01 A; compute 2πa V/I.

Explicit computation: 2math.pi5*0.12/0.01 = 376.9911 ohm m. Express the result in the units requested; perform any conversion before substituting into the equation rather than concealing a factor of one thousand in a final sentence. A rounded answer should not imply accuracy better than the synthetic inputs. This arithmetic operation provides a physical scale, not a certified subsurface interpretation.

Why the calculation matters. Apparent resistivity uses a geometric factor. Boundary of validity: Earth assumed homogeneous for direct interpretation. For an independent sensitivity check, change one stated input by ten percent, leaving others fixed. Explain how the result responds, and whether a plausible alternative hypothesis could produce the same scalar.

Practical: an observation that could disprove the preferred model

Choose an array spacing sequence to distinguish a shallow clay lens from a deeper conductor, then explain the inevitable trade-off in resolution.

Data stewardship. Document synthetic status, provenance, coordinates, units, processing version, license and intended use. Keep measurements in a separate file from interpolated or inversion-derived quantities. Draw a table containing at least two testable predictions, the needed instrument reading, its physical units, and one reason the observation might be biased.

Analysis note: how to decide whether the interpretation survives

Inspect the observation footprint. Ask which wavelengths and parts of the Earth contribute to the instrument. Potential fields integrate over substantial volumes, seismic picks follow ray paths and electrode arrays have diffuse sensitivity; a pixel-sized interpretation usually exceeds the resolving claim. In the context of resistivity principles and apparent resistivity, explain what a reviewer should expect to see in raw data if the mechanism is credible, and what unexpected pattern would require a revised geological model.

Synthesis and interpretation checkpoint

Write a short technical conclusion that separates three items: (1) what the instrument directly records in the electrical, electromagnetic and ground-penetrating methods problem; (2) the material property or process inferred from it; and (3) a geological claim that would still be unsupported. The exercise is not passed merely by obtaining 376.9911 ohm m. Explain why the uncertainty in the assumptions can exceed calculator precision and identify a discriminating field observation.

The physical mechanism and its measurable consequence

ERT estimates a resistivity distribution from many electrode combinations using numerical forward models and regularization. Depth of investigation is gradual and sensitivity decreases, not a sharp maximum depth. Smooth and blocky inverse priors can produce different geometries from the same data.

Acquisition protocol and observable variables

Compare observed and predicted apparent resistivity, inversion residual maps, reciprocal error and model smoothness. Mask poorly resolved cells.

The inverse ambiguity: alternative explanations

An appealing low-resistivity plume may be an inversion artefact, clay layer or unmodelled topography; it is not proof of contamination.

Worked calculation: synthetic classroom data

Question — Voltage reciprocal check. Forward apparent-resistivity datum is 110, reciprocal is 100 ohm m; report symmetric percent difference.

Explicit computation: 100*abs(110-100)/((110+100)/2) = 9.52381 %. Express the result in the units requested; perform any conversion before substituting into the equation rather than concealing a factor of one thousand in a final sentence. A rounded answer should not imply accuracy better than the synthetic inputs. This arithmetic operation provides a physical scale, not a certified subsurface interpretation.

Why the calculation matters. Reciprocal difference is an empirical quality metric. Boundary of validity: It measures consistency, not all inversion bias. For an independent sensitivity check, change one stated input by ten percent, leaving others fixed. Explain how the result responds, and whether a plausible alternative hypothesis could produce the same scalar.

Practical: an observation that could disprove the preferred model

Run both smooth and robust inversions and design a simple borehole sample that discriminates between them.

Acquisition design. Define the physical source, sensor property, site geometry and independent calibration needed before collecting data. State which raw channels must survive processing and what a second researcher could reproduce. Draw a table containing at least two testable predictions, the needed instrument reading, its physical units, and one reason the observation might be biased.

Analysis note: how to decide whether the interpretation survives

State the forward equation and its domain. Write the physical equation together with geometry, units and assumed constants. List a realistic departure that changes the prediction and distinguish a useful approximation from an unjustified physical claim. In the context of electrical resistivity tomography and inversion, explain what a reviewer should expect to see in raw data if the mechanism is credible, and what unexpected pattern would require a revised geological model.

Synthesis and interpretation checkpoint

Write a short technical conclusion that separates three items: (1) what the instrument directly records in the electrical, electromagnetic and ground-penetrating methods problem; (2) the material property or process inferred from it; and (3) a geological claim that would still be unsupported. The exercise is not passed merely by obtaining 9.52381 %. Explain why the uncertainty in the assumptions can exceed calculator precision and identify a discriminating field observation.

The physical mechanism and its measurable consequence

Alternating fields induce electric currents whose secondary response depends on conductivity, magnetic permeability, frequency, geometry and coupling. Simple skin-depth approximations assume a uniform half-space and become misleading in strongly layered or anisotropic media.

Acquisition protocol and observable variables

Record frequency, transmitter–receiver separation, coil orientation, phase, height and calibration standard. Check nearby metallic objects.

The inverse ambiguity: alternative explanations

A conducting utility can overshadow a regional aquifer signal. A shallow conductive layer may shield deeper information.

Worked calculation: synthetic classroom data

Question — Resistivity to conductivity. For resistivity 25 ohm m, compute sigma=1/rho in S/m.

Explicit computation: 1/25 = 0.04 S/m. Express the result in the units requested; perform any conversion before substituting into the equation rather than concealing a factor of one thousand in a final sentence. A rounded answer should not imply accuracy better than the synthetic inputs. This arithmetic operation provides a physical scale, not a certified subsurface interpretation.

Why the calculation matters. Reciprocal relation is for scalar bulk quantities. Boundary of validity: Anisotropic materials need tensor descriptions. For an independent sensitivity check, change one stated input by ten percent, leaving others fixed. Explain how the result responds, and whether a plausible alternative hypothesis could produce the same scalar.

Practical: an observation that could disprove the preferred model

Compare ideal skin depth at two frequencies and decide which could, in principle, interrogate a deeper target under the assumed half-space conditions.

Competing models. Choose two different geologic or instrumental explanations for the observation and write one numerical or spatial prediction on which they disagree. Reserve a new observation as a prospective test. Draw a table containing at least two testable predictions, the needed instrument reading, its physical units, and one reason the observation might be biased.

Analysis note: how to decide whether the interpretation survives

Decide a safe intended use. Clarify whether the result is a numerical exercise, exploratory screening, academic interpretation or professional decision product. Educational calculations do not confer authority to make engineering, emergency or environmental determinations. In the context of electromagnetic induction and skin depth, explain what a reviewer should expect to see in raw data if the mechanism is credible, and what unexpected pattern would require a revised geological model.

Synthesis and interpretation checkpoint

Write a short technical conclusion that separates three items: (1) what the instrument directly records in the electrical, electromagnetic and ground-penetrating methods problem; (2) the material property or process inferred from it; and (3) a geological claim that would still be unsupported. The exercise is not passed merely by obtaining 0.04 S/m. Explain why the uncertainty in the assumptions can exceed calculator precision and identify a discriminating field observation.

The physical mechanism and its measurable consequence

Radar reflections arise mainly from contrasts in dielectric properties, and subsurface travel time must be converted to depth using velocity. Penetration varies dramatically with conductivity, moisture, frequency and antenna coupling; conductive clay strongly limits it.

Acquisition protocol and observable variables

Collect common-offset profiles with time-zero calibration, trace spacing and antenna frequency; estimate velocity from hyperbolae or a measured reflector depth.

The inverse ambiguity: alternative explanations

A hyperbola may come from a pipe, stone or diffraction rather than a continuous stratigraphic horizon. Depth can be severely wrong when velocity varies.

Worked calculation: synthetic classroom data

Question — EM skin depth. Approximate skin depth 503 sqrt(rho/f) m, rho=100 ohm m and f=1000 Hz.

Explicit computation: 503*math.sqrt(100/1000) = 159.0626 m. Express the result in the units requested; perform any conversion before substituting into the equation rather than concealing a factor of one thousand in a final sentence. A rounded answer should not imply accuracy better than the synthetic inputs. This arithmetic operation provides a physical scale, not a certified subsurface interpretation.

Why the calculation matters. Frequency changes ideal electromagnetic penetration. Boundary of validity: Formula assumes nonmagnetic homogeneous half-space. For an independent sensitivity check, change one stated input by ten percent, leaving others fixed. Explain how the result responds, and whether a plausible alternative hypothesis could produce the same scalar.

Practical: an observation that could disprove the preferred model

Repeat an interpreted buried target depth at two plausible velocities and report a range rather than a single exact number.

Error budget. Separate repeatability, calibration bias, position error and model uncertainty in a table with units. Decide which uncertainty dominates the interpretation and what extra information would shrink it. Draw a table containing at least two testable predictions, the needed instrument reading, its physical units, and one reason the observation might be biased.

Analysis note: how to decide whether the interpretation survives

Track the reference frame. Coordinates, vertical datums, event origin times, magnetic epochs and polarity conventions are parts of the measurement. Confusing reference frames can shift all stations coherently without producing conspicuous random residuals. In the context of ground-penetrating radar and dielectric velocity, explain what a reviewer should expect to see in raw data if the mechanism is credible, and what unexpected pattern would require a revised geological model.

Synthesis and interpretation checkpoint

Write a short technical conclusion that separates three items: (1) what the instrument directly records in the electrical, electromagnetic and ground-penetrating methods problem; (2) the material property or process inferred from it; and (3) a geological claim that would still be unsupported. The exercise is not passed merely by obtaining 159.0626 m. Explain why the uncertainty in the assumptions can exceed calculator precision and identify a discriminating field observation.

The physical mechanism and its measurable consequence

Induced polarization records delayed electrical response that can be influenced by pore geometry, metallic particles, clays and fluid chemistry. Self-potential and streaming potentials can accompany subsurface water movement but have multiple sources. Multi-physics interpretation can be helpful without being decisive.

Acquisition protocol and observable variables

Measure time windows or phase spectra, background drift, salinity, temperature and field control data. Compare against direct borehole water samples where appropriate.

The inverse ambiguity: alternative explanations

A strong chargeability anomaly is not automatically an ore body, and low resistivity is not an unambiguous salinity estimate.

Worked calculation: synthetic classroom data

Question — Radar depth. Radar speed 0.10 m/ns and two-way echo time 120 ns; compute depth.

Explicit computation: 0.1*120/2 = 6 m. Express the result in the units requested; perform any conversion before substituting into the equation rather than concealing a factor of one thousand in a final sentence. A rounded answer should not imply accuracy better than the synthetic inputs. This arithmetic operation provides a physical scale, not a certified subsurface interpretation.

Why the calculation matters. Two-way radar conversion resembles seismic travel-time conversion. Boundary of validity: Velocity changes with moisture and dielectric structure. For an independent sensitivity check, change one stated input by ten percent, leaving others fixed. Explain how the result responds, and whether a plausible alternative hypothesis could produce the same scalar.

Practical: an observation that could disprove the preferred model

Design a noninvasive follow-up that chooses a second observable with different physical sensitivity from bulk DC resistivity.

Sensitivity study. Alter the assumed physical-property contrast, source depth or acquisition spacing. Determine whether the sign of the conclusion persists or whether the interpretation is contingent on a narrow prior. Draw a table containing at least two testable predictions, the needed instrument reading, its physical units, and one reason the observation might be biased.

Analysis note: how to decide whether the interpretation survives

Use a withheld test. Fit the parameters with one subset and predict genuinely unobserved measurements. Spatial clustering, repeated use of validation stations or tuning after inspecting the holdout exaggerates predictive evidence. In the context of polarization, fluids and discrimination, explain what a reviewer should expect to see in raw data if the mechanism is credible, and what unexpected pattern would require a revised geological model.

Synthesis and interpretation checkpoint

Write a short technical conclusion that separates three items: (1) what the instrument directly records in the electrical, electromagnetic and ground-penetrating methods problem; (2) the material property or process inferred from it; and (3) a geological claim that would still be unsupported. The exercise is not passed merely by obtaining 6 m. Explain why the uncertainty in the assumptions can exceed calculator precision and identify a discriminating field observation.

Imagine a valley where one electrode line finds a conductive zone. One hypothesis is saline groundwater, another is clay-rich sediment, and a third is a cultural conductor such as buried metal. All can lower apparent resistivity, and an inverse model might display a blue lens regardless of which interpretation is correct. A resistivity map is not a chemical analysis of water and cannot establish contamination by itself.

The measured voltage/current ratio depends on electrode geometry; the familiar Wenner formula converts it to apparent resistivity for an ideal array. Because the sensitivity kernel occupies a broad volume, changing electrode spacing changes depth sensitivity but not in the manner of a perfectly focused probe. Inversion imposes smoothness or other constraints to produce a model whose predicted apparent resistivities resemble data. A different regularization can yield a different-looking structure with similar residuals. Check reciprocal measurements, electrode contact and topographic treatment before trusting the inverse image.

Electromagnetic data provide another sensitivity, but field conductivity still depends on clay, water chemistry, connected porosity, frequency and geometry. A simple skin-depth equation is a useful scaling relation under a homogeneous nonmagnetic half-space. It is not a guaranteed exploration depth in a layered conductive valley. High-frequency ground-penetrating radar may be attenuated so rapidly in wet clay that a blank section is evidence of measurement limitation, not absence of an interface.

Design an integrated survey to distinguish the three hypotheses with minimal disturbance. One shallow borehole may reveal clay type, while fluid sampling directly tests dissolved ions. Independently map buried infrastructure and repeat measurements in a second season if moisture may change. The final interpretation should map where conductivity is inferred, where fluid salinity is actually measured, and where model resolution is inadequate. Conflating those layers creates an actionable-sounding but unsupported environmental claim.

Research aim. The invented scenario is: A hypothetical valley aquifer is imaged by resistivity electrodes, frequency-domain induction and radar lines. Clay, salinity and moisture have overlapping responses; students must avoid declaring contamination from conductivity alone. Treat datasets/unit-06.csv as forty fictional rows deliberately designed to exercise quality control and model criticism. Read the file header before drawing any map. Plot the two relevant measured fields from unit 6; display the quality_flag column rather than deleting all flagged stations silently. The synthetic table contains no geographic coordinate reference system or real site identifiers and cannot be used to identify a place.

Part A — raw measurements and dimensions. Read each column as a quantity, not as a decorative label. For electrical, electromagnetic and ground-penetrating methods, check whether the stated units are compatible with its equation and record which columns represent observed signals, instrument corrections, independent controls or model-derived quantities. Identify missing information such as date, site geology, true coordinate system, instrument orientation or calibration. If those omissions prevent a physical inversion, state the limitation and limit the task to transparent exploratory analysis. This is a stronger result than assigning fictitious precision.

Part B — repeat and quality flags. Count the rows marked review. Compare summary statistics, plotted gradients and model parameters with all forty records against the subset that passes the flag. Document every change; flags indicate inspection priority, not a universal mandate to delete. Check whether questionable values cluster along position or signal amplitude, because exclusions can introduce bias. If the noise appears correlated, do not assume a root-sum-square independent-error formula without a justification.

Part C — model competition. Describe two geophysical models that explain the synthetic observations in electrical, electromagnetic and ground-penetrating methods; they must differ in physical property, geometry or forcing rather than just parameter precision. For each hypothesis make an observable prediction about station response, uncertainty pattern, external control or independent method. Select one new measurement that would best discriminate them, and explain its predicted signal under both hypotheses. A new dense sample from the same short profile may add less independent information than a crossing line or sample.

Part D — numerical sensitivity. Recompute one of the unit worked cases with the principal input increased by 10% and decreased by 10%. Check if the response is linear, inverse, logarithmic, quadratic or subject to another nonlinearity. Show a table with input, computed output, unit and physical assumption. Do not claim a sensitivity of the geological interpretation itself if only the toy forward equation was changed; the inverse solution may have a much broader family of possibilities.

Part E — out-of-sample prediction and reporting. Split the synthetic observations by physical position, retaining a contiguous group as a blind holdout. Fit a deliberately simple forward trend to remaining rows and predict the held-out quantity without retuning. Compare predictive error with a baseline that uses a constant mean. State why geographic separation can be a stronger test than random point withholding. Report all transformations, quality flags, fitted parameters, residual plots and failure modes in a compact methods appendix.

Unit-specific evidence synthesis. Electrical resistivity is a bulk property influenced by interconnected water, ions, temperature, pore geometry, clay surface conduction and mineral phases. Four-electrode apparent resistivity is not the resistivity at a single sampled depth. Electromagnetic surveys sense conductivity through inductive coupling, and radar travel time depends on effective electromagnetic velocity and attenuation. Different array configurations sample different sensitivity regions, so apparent contrasts require geometry-aware interpretation. Saltwater and clay can both yield high conductivity; dry rock, air voids and clean sand may give resistive signatures. Instrument coupling, buried utilities, cultural noise and inversion regularization add further uncertainty. Map the chain of evidence from sensor to interpretation; highlight an uncertainty not reducible by extra decimal places. In your conclusion, indicate which measurement you would acquire next if the synthetic experiment were replaced with licensed real-world data. Keep any real data sourced from the references separate from these forty invented rows.

Suggested report. Submit a transparent CSV-reading script or spreadsheet calculation, one raw-data figure, a sensitivity plot, two candidate geological explanations, a withheld-test result and a one-page uncertainty ledger. Grade on conceptual correctness (25%), numerical units and reproducibility (25%), alternate hypotheses (20%), validation logic (20%) and scientific communication (10%). Never reward a map for appearing more certain than its underlying observations.

1. Identify which part of electrical, electromagnetic and ground-penetrating methods would remain ambiguous even if measurement noise were reduced to zero. Explain how the survey geometry or parameter relationship creates that ambiguity.

1. Find a case where a large change in the observed signal would not imply a proportionally large change in target size. Write the forward relation, name fixed parameters, and state why inversion differs from prediction.

1. Propose a new survey that would risk rejecting the preferred geological interpretation. Its value should be based on discriminating predictions rather than its technological sophistication.

1. Draw a data lineage chart from raw instrument data through correction, inversion and geological claim. Underline any operation that assumes a material property not measured in the original dataset.

1. Explain why the synthetic data and interactive sliders are educational examples and cannot support emergency warnings, ground-investigation certification, resource appraisal or environmental enforcement.

Use REFERENCES.md to find a genuine program or open repository relevant to electrical, electromagnetic and ground-penetrating methods. Before downloading any data, check the provider documentation, citation requirements, usage permissions and the time period represented. Record the instrument type, station metadata and acquisition geometry, then reproduce one simple physical calculation from scratch. The objective is not to match a published colored inversion image exactly; it is to understand the observable and give a defensible account of differences.

Open this chapter in the HTML reader to change one controlled physical variable at a time. Each model is a simplified forward calculation, not an inversion or a real-world estimator.

Wenner apparent resistivity

Input quantity: Electrode separation measured in m; demonstration range 1–20, initial setting 5. Output: ohm m. Computed from 2Math.PIx*0.12/0.01.

Unmodelled influences and interpretation boundary: Fixed current/voltage and ideal geometry; result apparent, not a layer property. Students should articulate a case where this formula yields a numerically precise but physically misleading answer and state an independently measurable quantity needed for a more complete model.

EM half-space skin depth

Input quantity: Frequency measured in Hz; demonstration range 10–10000, initial setting 1000. Output: m. Computed from 503*Math.sqrt(100/x).

Unmodelled influences and interpretation boundary: Nonmagnetic uniform half-space, resistivity 100 ohm m. Students should articulate a case where this formula yields a numerically precise but physically misleading answer and state an independently measurable quantity needed for a more complete model.