Spatial Distribution of Trace Metals in Soil
A geostatistical investigation mapping arsenic, lead, cadmium, and nickel across the Abuakwa South Municipal area, using Ordinary Kriging and spline interpolation to turn scattered soil samples into continuous contamination surfaces.
The problem
Soil contamination from trace metals is a real concern in areas touched by small-scale mining and shifting land use. But contamination is invisible between sample points: you can measure a handful of locations in a lab, yet a decision-maker needs to understand the whole area, including the ground nobody sampled.
This project asks a focused question: given a limited set of soil samples across Abuakwa South, how do you responsibly estimate contamination everywhere else, and how confident can you be in those estimates?
Approach
The work targets four trace metals of environmental concern: arsenic (As), lead (Pb), cadmium (Cd), and nickel (Ni). The core of the method is geostatistical interpolation, comparing two techniques rather than trusting one blindly.
- Ordinary Kriging, a geostatistical method that predicts values at unsampled locations by weighting nearby samples according to a modeled spatial correlation structure (the variogram). Crucially, it also returns a prediction variance, so you know where you are confident and where you are guessing.
- Spline interpolation, a smoother, curve-fitting approach used as a comparison, to see where each method agrees and where they diverge.
The comparison matters. Naive methods can smear a sharp contamination plume into a soft blur. Kriging, by honoring spatial autocorrelation and the sample values at their true locations, preserves structure that smoothing methods lose.
It accounts for spatial autocorrelation, returns a prediction variance so uncertainty is explicit, and honors the measured values at their exact locations, revealing patterns that inverse-distance and spline methods tend to smooth away.
Tools & workflow
The pipeline runs end to end: sample data preparation, coordinate reference alignment for a Ghanaian study area, variogram modeling, interpolation, and cartographic output. It combines the desktop GIS environment with scripted analysis.
What the maps show
The interpolated surfaces turn point measurements into readable contamination maps, one per metal, alongside their uncertainty. The value is interpretive: hot spots are not random. They track the underlying geology, the land use, and the small-scale mining footprint of the area.
Framed this way, the maps become a screening tool. They highlight where contamination is likely elevated, which supports prioritizing further sampling and remediation attention rather than treating the whole area uniformly.
Why it matters
Environmental data in Ghana is often underused. A methodology like this is repeatable, transparent, and grounded in open tools, which means the same workflow can be applied to other districts and other contaminants. The project is as much about demonstrating a defensible spatial method as it is about this one study area.
This research is ongoing. Full methodology, maps, and interpretation will be published in due course.
Want the technical detail or the maps as they finalize?
Get in touch