Satellites and Game Theory Reveal What Really Controls Groundwater Quality in Iran

Deep beneath the semi-arid plains of western Iran, the water that sustains more than half of the country’s consumption has long been monitored the hard way: technicians driving to piezometer wells, hauling samples back to laboratories, and waiting days for chemical analyses. Now a research team has shown that a European satellite orbiting 700 kilometres overhead, combined with hybrid machine learning and a mathematical trick borrowed from cooperative game theory, can rank the chemical parameters that matter most in groundwater — and do it cheaply, quickly, and across entire landscapes. The study, published in Environmental Advances, focuses on the Khorramabad Plain in Lorestan Province and could reshape how water-stressed regions around the world keep watch over their hidden aquifers.

The stakes are enormous. More than two billion people worldwide depend on groundwater as their primary drinking source, and it supplies roughly 43 percent of global agricultural water demand. In Iran, where abstraction already exceeds natural recharge, groundwater accounts for over 55 percent of total water consumption. Population growth, intensified agriculture, urbanisation, and climate change are pushing these aquifers toward declining water tables, land subsidence, saltwater intrusion, and contamination that is often irreversible. Yet traditional monitoring — accurate as it is — remains slow, expensive, and too sparse to capture the spatial patterns that determine where problems actually emerge.

The core challenge the researchers faced is fundamental: groundwater is invisible from space. Unlike lakes and rivers, whose surfaces reflect light in ways satellites can decode directly, aquifers lie buried under soil and rock. The team, led by Mozhgan Yarahmadi with Ali Haghizadeh and Leila Ghasemi, circumvented this with an indirect, surrogate-based strategy. Their reasoning was that groundwater quality is tightly coupled to what happens at the surface — land use, vegetation, soil moisture, temperature, and rainfall. If satellites can see the surface, and machine learning can learn the hidden relationships between surface signals and underground chemistry, then the orbiting sensors become proxies for what lies below.

To build those proxies, the researchers turned to Sentinel-2, the European Space Agency’s twin-satellite Copernicus mission offering thirteen spectral bands at 10 to 60 metre resolution and a five-day revisit cycle. They extracted eight bands — from blue through near-infrared — for dates matching field sampling within a ten-day window, and computed five spectral indices designed to capture chlorophyll content, turbidity, and surface moisture. These were fused with climatic variables from the ERA5-Land reanalysis dataset, including precipitation and land surface temperature. Around each of five piezometer stations, a five-kilometre buffer zone averaged the pixel values, smoothing the mismatch between point-based well measurements and satellite footprints. The result was a dataset of 56 valid samples collected seasonally from 2015 to 2024, each paired with sixteen environmental predictor variables.

With the dataset assembled, the team deployed nine competing models: five base algorithms — Random Forest, XGBoost, Support Vector Regression, Multi-Layer Perceptron, and Linear Regression — and four hybrid architectures that combined them, including a Voting Ensemble averaging four base models and a Stacking Ensemble feeding base-model outputs into a Random Forest meta-model. The standouts were striking. The Voting Ensemble achieved a coefficient of determination of 0.986 for sodium at the Sabour station and 0.974 for electrical conductivity at Gosheh, while the SVR+MLP hybrid reached 0.900 for total dissolved solids at Robat. Hybrid models claimed five of the six top-performing station-model-parameter combinations, confirming that blending complementary algorithms consistently outperforms any single learner.

The pattern behind those numbers is physically meaningful. Parameters the researchers classify as optically active — sodium, electrical conductivity, total dissolved solids, and chloride — achieved the highest accuracies, because dissolved salts alter the refractive index and light absorption of water in ways Sentinel-2’s spectrometer can detect. Semi-optical parameters such as calcium, magnesium, and total hardness landed in a moderate range, while non-optical parameters like pH, sodium percentage, and sulfate proved hardest to predict, lacking any direct spectral signature. Cross-validation told an honest story: while point-specific performance soared above 0.97, stable regional performance settled between 0.52 and 0.55 for the main parameters — respectable given the modest sample size and the inherent spatial variability of aquifer chemistry.

The study’s most novel contribution lies in how it weighted the importance of each quality parameter. Rather than relying on a single importance metric, the team combined SHAP — SHapley Additive exPlanations, a technique rooted in cooperative game theory that fairly distributes each prediction among contributing features — with classical Shapley values derived from linear regression. The final weight for each parameter was the normalised average of the two. The verdict was unambiguous: electrical conductivity dominated with 50.47 percent of the total weight, followed by total dissolved solids at 32.40 percent. Total hardness and sodium percentage ranked third and fourth at 7.87 and 7.56 percent, while the remaining six parameters each fell below half a percent. In other words, just two indicators — both optically active and strongly correlated with satellite reflectance — carry more than 82 percent of the explanatory power.

The researchers then stress-tested their rankings. Seasonal analysis revealed that while conductivity and dissolved solids held steady across spring, summer, and autumn, total hardness surged in summer — climbing from 17.12 percent in spring to 21.15 percent — likely driven by higher temperatures intensifying carbonate dissolution in the alluvial aquifer. Sodium percentage also crept upward in summer, a plausible fingerprint of agricultural fertiliser use during the growing season. Bootstrap resampling with 1,000 iterations quantified the uncertainty: conductivity and dissolved solids showed the lowest coefficients of variation, around 9.7 to 9.9 percent, confirming their weights are stable and reliable, while chloride and pH carried the highest uncertainty, consistent with their weak spectral linkage.

The authors are candid about limitations. Fifty-six samples is a small dataset, and reanalysis-based climatic data may diverge from ground stations; atmospheric correction remains tricky under cloudy skies, and models trained on Khorramabad’s geology cannot simply be transplanted elsewhere without local retraining. Yet the practical implications are compelling. Because Sentinel-2 imagery and ERA5-Land data are freely available and continuously refreshed, the trained models can be extrapolated to unmonitored areas, generating spatial maps of conductivity, dissolved solids, sodium, and chloride where no wells exist. For water managers in semi-arid regions, the framework offers an interpretable, low-cost decision-support tool: prioritise conductivity and dissolved solids as headline indicators, intensify hardness monitoring during summer, and treat satellite-based assessment as a powerful complement — not a replacement — for laboratory sampling. As aquifers worldwide face mounting pressure, watching them from orbit may soon become standard practice.

Subject of Research: Prediction and weighting of groundwater quality parameters using Sentinel-2 remote sensing, hybrid machine learning models, and SHAP-game theory analysis in Khorramabad, Iran

Article Title: Groundwater Quality Parameter Weighting Using Hybrid Machine Learning and SHAP-Game Theory Based on Sentinel-2 Data (Case Study: Khorramabad, Iran)

Article References: Yarahmadi, M., Haghizadeh, A., & Ghasemi, L. (2026). Groundwater Quality Parameter Weighting Using Hybrid Machine Learning and SHAP-Game Theory Based on Sentinel-2 Data (Case Study: Khorramabad, Iran). Environmental Advances, Article 100770. https://doi.org/10.1016/j.envadv.2026.100770

Image Credits: AI Generated

DOI: 10.1016/j.envadv.2026.100770

Keywords: groundwater quality, Sentinel-2, machine learning, SHAP, game theory, remote sensing, hybrid models, Iran, water monitoring, electrical conductivity, total dissolved solids, semi-arid regions