AI Stacking Model Pinpoints Copper Deposits in Iran With Striking Accuracy

Artificial intelligence has just delivered one of its most impressive performances yet in the hunt for buried treasure beneath the Earth’s surface. In a study published in Natural Resources Research, researchers Elnaz Geravandi of Kharazmi University and Reza Ghezelbash of the University of Tehran unveiled a multi-level stacking ensemble architecture that maps the likelihood of hidden porphyry copper-gold deposits across the Ahar-Arasbaran metallogenic belt in northwest Iran. The model achieved an area under the curve of 0.99, an accuracy of 0.95, precision of 0.94, recall of 0.98, and an F1-score of 0.96, while flagging roughly 86.8 percent of known porphyry occurrences within just 11.36 percent of the highest-ranked prospectivity zones. Those numbers translate into a remarkably efficient targeting tool: exploration teams could concentrate their expensive drilling and field campaigns on a small fraction of the landscape and still capture the overwhelming majority of known mineralized sites.

Porphyry copper deposits are the world’s principal source of copper and a major source of gold and molybdenum, forming when metal-rich magmatic fluids rise from deep intrusions and precipitate ore minerals in large, diffuse zones near the surface. Finding new ones is notoriously difficult because the signatures they leave behind are subtle, overlapping, and nonlinear. Geochemical anomalies in stream sediments interact with geology, fault networks, and hydrothermal alteration patterns in ways that simple statistical methods struggle to untangle. Datasets also suffer from multicollinearity, where different evidence layers carry redundant information, and from spatial dependence, meaning that samples collected close together are not truly independent. These are precisely the conditions under which machine learning, and ensemble methods in particular, tend to outperform traditional approaches.

The heart of the new framework is a technique called stacking, an idea that dates back to David Wolpert’s 1992 work on stacked generalization. Instead of betting on a single algorithm, stacking trains several base learners on the same problem and then uses their predictions as inputs to a higher-level model that learns how best to combine them. Geravandi and Ghezelbash pushed this concept further with a hierarchical, multi-level architecture. Five base learners were deployed: random forest, Light Gradient Boosting Machine (LightGBM), Extreme Gradient Boosting (XGBoost), support vector regression, and a multilayer perceptron neural network. Each brings a different inductive bias to the table. Random forests average many decorrelated decision trees to suppress variance, gradient boosting machines sequentially correct the errors of weak learners to reduce bias, support vector regression finds flexible boundaries in high-dimensional feature space, and neural networks capture intricate nonlinear relationships among evidential layers.

The raw material feeding these algorithms was as important as the algorithms themselves. The researchers processed stream sediment geochemical data from 2,716 samples, a dataset capable of revealing spatially coherent multi-element anomalies that reflect the upstream footprints of porphyry mineralization. Stream sediments act as natural sampling nets: metals eroded from mineralized zones are transported downstream and concentrated in drainage sediments, so anomalous concentrations of copper, gold, and pathfinder elements can point prospectors back toward their sources. These geochemical layers were integrated with geological, structural, and hydrothermal alteration evidence layers within a unified geospatial machine learning environment, allowing the models to weigh lithology, fault density, and alteration minerals alongside chemistry.

A critical methodological innovation was the use of spatial block cross-validation rather than conventional random data splitting. Because neighboring locations share similar conditions, randomly splitting spatial data into training and test sets can leak information across the boundary and inflate performance estimates, a phenomenon known as spatial autocorrelation bias. By dividing the study area into spatial blocks and validating across them, the researchers ensured that the reported metrics reflect genuine generalization to unseen terrain. The data were split 70 percent for training and 30 percent for validation under this spatially constrained scheme, and distance-based spatial analysis and prediction-area (P-A) plots were used to evaluate how well each model balanced the proportion of deposits correctly predicted against the area of land flagged as prospective.

The P-A plot also guided feature engineering. Two feature configurations were constructed based on the quantitative importance of the evidential layers, and the comparison produced a nuanced finding. Refining features using P-A plot guidance did improve the performance of individual base models, trimming away layers that added noise rather than signal. Yet the full multi-level stacking framework demonstrated that comprehensive integration of all evidence improved predictive balance and spatial coherence more than aggressive feature reduction. In other words, when a well-designed ensemble learns how to weight diverse information, seemingly redundant or weak layers can still contribute to a more geologically plausible final map. This challenges a common instinct in applied machine learning, where pruning inputs is often assumed to be inherently beneficial.

The hierarchical stacking stage then fused the predictions of the five base learners into a single consensus prospectivity map. The result was not merely a statistical improvement but a spatially more coherent one: high-prospectivity zones aligned more cleanly with the known architecture of the Ahar-Arasbaran belt, a Cenozoic volcanic arc that hosts significant porphyry copper-molybdenum-gold systems, including the well-studied Sungun deposit. Importantly, the authors emphasize that the 86.8 percent capture rate within 11.36 percent of the map area reflects enhanced spatial targeting efficiency rather than predictive certainty, a careful framing that distinguishes exploration prioritization from guarantees of discovery.

Concerns about overfitting, the perennial bogeyman of high-performing machine learning models, were addressed directly. An AUC of 0.99 might raise eyebrows in fields where such scores often signal data leakage, but the combination of spatial block cross-validation, distance-based analysis, and P-A plot evaluation provides converging lines of evidence that the model’s performance is robust rather than artifactual. The authors also report that the framework demonstrates strong generalizability and can be transferred to other regions with different scales and mineralization types, suggesting the architecture is not tailored to the quirks of a single belt. Reinforcing that claim, the Python scripts and anonymized demonstration datasets needed to reproduce the entire workflow have been released publicly on GitHub, an unusually transparent step that allows other researchers to stress-test and adapt the method.

The broader implications extend well beyond northwest Iran. Global copper demand is projected to surge as electrification, renewable energy infrastructure, and grid expansion accelerate, yet discovery rates for new porphyry deposits have lagged for decades because the easy targets near the surface have largely been found. Machine learning prospectivity mapping offers a way to re-examine vast archives of legacy geochemical, geological, and remote sensing data through a fresh computational lens, prioritizing ground that previous generations of explorers may have undervalued. The study builds on a growing body of work applying random forests, gradient boosting, deep learning, and ensemble strategies to mineral exploration, but its multi-level stacking design, spatial validation rigor, and open release of code set a benchmark for how such studies should be conducted. If the framework transfers as well as its authors suggest, the dusty stream sediments of other mountain belts around the world may soon be whispering the locations of the next generation of copper mines, and they will be whispering it through the mathematics of stacked ensembles.

Subject of Research: Multi-level stacking ensemble machine learning for porphyry copper-gold mineral prospectivity mapping in northwest Iran

Article Title: A Multi-Level Stacking Ensemble Architecture: Advantages of Random Forest, Extreme Gradient Boosting, and Light Gradient Boosting Machine for Porphyry-Related Al-Based Mineral Prospectivity Mapping

Article References: A Multi-Level Stacking Ensemble Architecture: Advantages of Random Forest, Extreme Gradient Boosting, and Light Gradient Boosting Machine for Porphyry-Related Al-Based Mineral Prospectivity Mapping. (n.d.). https://doi.org/10.1007/s11053-026-10763-3

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10763-3

Keywords: mineral prospectivity mapping, machine learning, random forest, XGBoost, LightGBM, multi-level stacking, porphyry copper deposits, spatial cross-validation, stream sediment geochemistry, Ahar-Arasbaran belt, prediction-area plot, exploration targeting