Geostatistics · Final-Year Research

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.

Role
Lead Researcher
Context
Final-Year Project, KNUST
Location
Abuakwa South, Ghana
Status
Ongoing
Contour map visualization of spatial data

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.

Why kriging over naive methods

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.

QGIS Python Ordinary Kriging Spline Interpolation Variogram Modeling Environmental Geochemistry

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