The Bou Azzer–El Graara inlier occupies a special place in the geology of North Africa. It lies within the Pan-African suture zone, the scar left behind when an ancient ocean closed during the assembly of the West African Craton more than half a billion years ago. The inlier preserves one of the most representative Neoproterozoic ophiolite successions known, a fragment of oceanic crust and mantle rocks thrust onto continental margins during arc–continent collision. These same serpentinite bodies later became the host of the famous cobalt–nickel arsenide ores that have been mined at Bou Azzer for nearly a century, making the district a globally significant source of cobalt and a natural laboratory for studying how ancient plate boundaries control the distribution of metals.

Despite decades of geological mapping, structural analysis and geophysical surveying in the region, no previous study had combined quantitative machine learning classification of magnetic edges with three-dimensional magnetic inversion inside a single multiscale structural framework. That gap is what motivated Ayoub Soulaimani of Morocco’s Directorate of Geology, Mines and Hydrocarbons, Saâd Soulaimani of the Rabat School of Mines, and Abdelhalim Miftah of Hassan First University of Settat to design a four-stage workflow built around high-resolution aeromagnetic data acquired under the Moroccan national airborne survey program. The result is a methodology that turns raw magnetic measurements into a ranked, probability-based map of where hidden structures, and potentially hidden ores, are most likely to lie.

The first stage of the workflow involved computing a suite of potential-field derivatives from the total magnetic intensity data. These included the horizontal gradient, the analytic signal, three orthogonal first-order directional gradients, the tilt derivative, the horizontal gradient of the tilt angle, and a residual map reduced to the pole. Each of these filters emphasizes different aspects of the magnetic field. The horizontal gradient and analytic signal highlight the edges of magnetized bodies, the directional gradients accentuate lineaments oriented in specific azimuths, and the tilt derivative balances the responses of shallow and deep sources so that weak signals from deeply buried contacts remain visible. Together they form a multidimensional description of where the magnetic fabric of the region changes abruptly, which in turn corresponds to lithological boundaries and faults.

Those derivative maps then served as the input features for training a feed-forward artificial neural network, optimized through the Levenberg–Marquardt backpropagation algorithm, a classical and computationally efficient training method for small to medium networks. The researchers tested two weighting schemes for the input data: one in which all derivatives contributed equally, and another in which derivatives known to be especially sensitive to lithological contacts received selective weighting. Both configurations reached an accuracy of roughly 86 percent, and the very similar area-under-the-curve values obtained in the two cases suggested to the authors that the solution is robust rather than an artifact of any particular parameterization. In other words, the network was genuinely learning the structural signature of the region, not merely fitting noise.

The output of the trained network is a set of weighted probability surfaces that quantify, pixel by pixel, the likelihood that a structural discontinuity exists at that location. When compared against existing geological maps, these surfaces reproduced all of the major structural corridors already documented in the Bou Azzer region, providing an independent validation of the approach. More intriguingly, they also revealed zones of structural complexity that had never been mapped before, places where the magnetic fabric implies intersecting fault sets and rotated blocks that surface geology alone had not resolved. It is precisely such corridors of intersecting structures that tend to channel hydrothermal fluids and localize mineralization, which makes the newly identified zones prime candidates for future exploration.

In parallel with the machine learning analysis, the team performed a three-dimensional inversion of the magnetic data over the Foum Zguid dyke zone, one of the most striking magmatic features of the Anti-Atlas. The inversion reconstructed the volume distribution of magnetic susceptibility and allowed the researchers to recover the geometry of the dyke structure down to a depth of about two kilometers. This kind of quantitative depth constraint is rare in regional structural studies and demonstrates how modern inversion algorithms, originally developed for mineral deposit scale problems, can be scaled up to interrogate the deep architecture of entire tectonic corridors.

The final synthesis combined gradient analysis, Euler deconvolution for estimating source depths, the neural network prospectivity results, and the three-dimensional susceptibility model into a single high-resolution structural interpretation. The dominant fault orientations that emerged are north–south and WNW–ESE, a pattern that carries direct geodynamic significance. The north–south set is consistent with structures related to the Pan-African collision and its later reactivation, while the WNW–ESE set aligns with the major shear zones that dissect the Anti-Atlas belt. The study thus links the fine-scale magnetic fabric of the inlier to the broad tectonic evolution of the northern margin of the West African Craton, from Neoproterozoic ocean closure through Paleozoic deformation.

The economic implications are immediate. Bou Azzer remains one of the few operating cobalt arsenide districts in the world, and cobalt has become a strategically critical metal for battery technologies. The cobalt–nickel–arsenic ores at Bou Azzer are hosted in serpentinites along fault-controlled hydrothermal veins, so the identification of previously unrecognized structural corridors directly translates into new underexplored target zones for cobalt, copper, silver and gold. Rather than drilling blindly, exploration geologists can now prioritize the areas where the neural network’s probability surfaces intersect favorable host rocks and where Euler deconvolution indicates structures extending to mineable depths.

Beyond Morocco, the study offers a template that other exploration geologists and structural mappers can adapt. The workflow requires nothing more exotic than a standard aeromagnetic dataset, a set of well-established potential-field filters, and a modest feed-forward neural network, yet it delivers validated structural probability maps and quantitative three-dimensional models in a single coherent framework. As airborne magnetic surveys become cheaper and machine learning tools become more accessible, the approach demonstrated in the Central Anti-Atlas suggests that many of the world’s poorly mapped Precambrian terranes could be re-examined at high resolution, revealing both the deep tectonic history recorded in their magnetic fabrics and the mineral deposits that history has left behind.

Subject of Research: Machine learning analysis of aeromagnetic data for structural mapping and mineral prospectivity in the Moroccan Anti-Atlas

Article Title: Structural and tectonic analysis of the Moroccan Central anti-atlas using predictive targeting with machine learning and inversion of aeromagnetic data: geodynamic and mining implications

Article References: Soulaimani, A., Soulaimani, S., & Miftah, A. (2026). Structural and tectonic analysis of the Moroccan Central anti-atlas using predictive targeting with machine learning and inversion of aeromagnetic data: geodynamic and mining implications. Earth Science Informatics, 19(11), Article 203. https://doi.org/10.1007/s12145-026-02257-w

Image Credits: AI Generated

DOI: 10.1007/s12145-026-02257-w

Keywords: aeromagnetic data, machine learning, artificial neural network, Bou Azzer, Anti-Atlas, Morocco, cobalt mineralization, magnetic inversion, structural geology, Pan-African suture, mineral prospectivity, geophysics

The Mahalanobis distance, introduced by statistician P. C. Mahalanobis in 1936, measures how far a point sits from the center of a data cloud, accounting for the shape and orientation of that cloud through its covariance structure. When a multivariate dataset follows a single multivariate normal distribution, the squared Mahalanobis distance, MD², should obey a chi-square distribution with degrees of freedom equal to the number of variables. That theoretical link is the backbone of multivariate outlier detection: points whose MD² exceeds a chosen chi-square quantile are declared anomalous. In exploration geochemistry, this logic underlies the identification of anomalies—samples whose elemental signatures depart so far from background that they may betray concealed mineralization beneath the surface.

The problem, as the new study demonstrates through simulation experiments, is that real geochemical surveys almost never consist of a single homogeneous population. Stream sediment samples aggregate material from different rock types, weathering regimes and geological histories, so the dataset is a blend of several overlapping populations. The researchers showed that this multiple population mixing distorts the global hyperellipsoid that the Mahalanobis distance implicitly fits to the data in three characteristic ways: the centroid shifts away from any single population’s true center, the ellipsoid stretches anisotropically along directions driven by between-population differences rather than within-population variability, and a thin-shell effect emerges in which samples pile up at intermediate distances from the center. Each distortion biases MD² away from its theoretical chi-square distribution, so the thresholds that practitioners use to declare anomalies are calibrated against a distribution the data no longer follow.

The consequence is subtle but consequential. Because the global covariance matrix inflates to accommodate the spread between populations, genuinely anomalous samples can be swallowed into an artificially widened ellipsoid, while ordinary samples from an underrepresented subpopulation can be unfairly flagged. In geochemical terms, lithology-driven differences—such as the chemical fractionation of granitic magmas—can manufacture false anomalies that have nothing to do with mineralization, while true hydrothermal signatures may be suppressed. The team’s simulations made these failure modes visible: whenever mixed populations were injected into otherwise well-behaved synthetic data, the histogram of MD² values peeled away from the chi-square reference curve in predictable, geometry-dependent patterns.

To counter this, the authors propose a framework that attacks the problem at its root by first unmixing the populations before any distances are computed. The first stage employs a Gaussian mixture model estimated with the expectation–maximization algorithm, EM-GMM, to decompose the heterogeneous dataset into geochemically coherent subpopulations. The expectation–maximization algorithm, originally formalized by Dempster, Laird and Rubin in 1977, iteratively assigns samples probabilistic memberships across a set of Gaussian components and refines each component’s mean and covariance until convergence. In effect, the algorithm reverses the mixing process, recovering the constituent distributions whose overlap had corrupted the global statistics.

The second stage then applies the fast minimum covariance determinant method, Fast-MCD, within each identified subpopulation. Fast-MCD, developed by Rousseeuw and Van Driessen in 1999, searches for the subset of roughly half the data whose covariance matrix has the smallest determinant, yielding robust estimates of location and dispersion that resist contamination by outliers. Computing local, robust MD² values within each unmixed subpopulation substantially restores agreement with the chi-square distribution, because each subpopulation now satisfies the single-population assumption the classical statistic depends on. The framework thus pairs a probabilistic unmixing tool with a robust distance estimator, addressing both the cause and the symptom of the geometric distortion.

The researchers tested the framework on lithium exploration in northeastern Hunan Province, China, a region within the Jiangnan orogenic belt known for rare-metal pegmatites, including the giant Renli Nb-Ta deposit. Using 16 lithology-indicating elements—SiO₂, Al₂O₃, Fe₂O₃, MgO, CaO, K₂O, Na₂O, Ba, Ni, Sr, Th, Ti, U, V, Y and Zr—they clustered 2,447 stream sediment samples into six subpopulations. For each subpopulation, they computed local robust MD² values for the lithium metallogenic association of Li, Be, Sn, F and Bi. The resulting local MD² values aligned closely with the theoretical chi-square distribution with five degrees of freedom, confirming that the unmixing step had done its job and that anomaly thresholds could once again be set on solid statistical ground.

The performance gains were substantial. A prediction–area plot, a standard tool for evaluating prospectivity maps, showed that the proposed framework captured 92 percent of known deposits within just 7 percent of the study area. It outperformed global Fast-MCD applied without unmixing, classical Mahalanobis distance, adaptive MD², a one-class k-nearest-neighbor method known as AMSD-kNN, and a deep autoencoder model. When the team applied the conventional chi-square threshold at the 97.5th percentile for five degrees of freedom, the delineated anomalies shrank the target area to 12 percent of the study region while still identifying 100 percent of the known deposits—a combination of focus and completeness that is exceptionally valuable in mineral exploration, where drilling and fieldwork are expensive.

Beyond the raw numbers, the method demonstrated geological judgment. It suppressed false anomalies induced by lithology-driven processes, specifically the magmatic fractionation of granites that can mimic elemental enrichment patterns, and instead highlighted anomalies consistent with concealed hydrothermal mineralization controlled by fault structures and contact zones. That spatial pattern matters because lithium in the region is associated with rare-element pegmatites whose emplacement is structurally controlled, so anomalies tracing faults and granite contacts carry genuine exploration significance rather than lithological noise.

The study’s broader lesson reaches past geochemistry. Mahalanobis distance and chi-square thresholds are workhorses across the sciences, from structural health monitoring to chemometrics to machine-learning anomaly detection, and all of these applications inherit the same fragility when data are mixtures rather than single populations. By naming the distortions—centroid shift, anisotropic stretching and the thin-shell effect—and showing that a disciplined unmixing-then-robust-estimation pipeline restores theoretical behavior, the researchers have turned an abstract geometric insight into a practical workflow. For explorers hunting lithium and other critical metals in complex terrains, it means the maps that guide their next drill campaign can finally be drawn by a statistician’s tool that no longer lies about the shape of the data beneath it.

Subject of Research: Geometric distortion of the Mahalanobis distance under mixed populations and its mitigation for geochemical anomaly identification.

Article Title: Geometric Distortion Mechanisms of Mahalanobis Distance Under Multiple Population Mixing and Mitigation Strategies: A Case Study on Geochemical Anomaly Identification

Article References: Zhang, Q., Lu, J., Tian, M., Sun, L., Zhao, X., & Shi, Y. (2026). Geometric Distortion Mechanisms of Mahalanobis Distance Under Multiple Population Mixing and Mitigation Strategies: A Case Study on Geochemical Anomaly Identification. Natural Resources Research. https://doi.org/10.1007/s11053-026-10772-2

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10772-2

Keywords: Mahalanobis distance, geochemical anomaly, chi-square distribution, Gaussian mixture model, Fast-MCD, lithium exploration, mineral prospectivity, population mixing, outlier detection, Natural Resources Research, Geometric, Distortion

The core challenge the researchers set out to solve is one that plagues mineral exploration worldwide. Stream sediment geochemical surveys, which measure the concentrations of dozens of elements in sediment samples collected from drainage networks, produce vast datasets. Yet only a tiny fraction of the sampled locations can be confidently labeled as mineralized, because confirmed deposits are rare and expensive to verify. Standard supervised classifiers, which learn from labeled examples, struggle when the positive class is so sparsely represented. Meanwhile, the overwhelming majority of samples remain unlabeled, a reservoir of information that conventional methods simply ignore.

Zheng and Chen’s framework attacks this problem from three directions at once. First, LightGBM serves as the base classifier, chosen for its ability to capture the nonlinear relationships between element concentrations and the spatial distribution pattern of mineral deposits. LightGBM is a gradient boosting decision tree algorithm known for its speed and memory efficiency, achieved through techniques such as histogram-based splitting and leaf-wise tree growth. In mineral exploration, where the interplay between trace element signatures and ore-forming processes is anything but linear, this flexibility matters enormously.

Second, the framework employs a self-training algorithm, a semi-supervised technique in which a model is initially trained on the small labeled set and then iteratively predicts labels for the unlabeled data. The most confident predictions are added to the training pool, and the model is retrained, gradually bootstrapping its way toward a more complete understanding of the data. Self-training has a long pedigree in machine learning, dating back to the 1960s, but its application to geochemical anomaly recognition is relatively recent. By exploiting the vast unlabeled portion of stream sediment data, the method effectively converts a data-poor problem into a data-rich one.

Third, the researchers introduced the Synthetic Minority Oversampling Technique, or SMOTE, to address the severe class imbalance that would otherwise cripple the initial training stage. Because mineralized locations are so rare, a classifier trained naively would learn to simply predict that nothing is mineralized, achieving high accuracy while being useless in practice. SMOTE works by generating synthetic examples of the minority class, interpolating between existing positive samples in feature space rather than merely duplicating them. This enlarges the number of positive samples enough for LightGBM to establish a meaningful model during the first round of self-training, after which the iterative process takes over.

To test the framework, the team built four models on interpolated 1:50,000-scale stream sediment geochemical data from Moridawa: the self-training LightGBM, a self-training support vector classifier, a standalone LightGBM, and a standalone support vector classifier. This head-to-head comparison allowed the researchers to isolate the contributions of both the self-training strategy and the choice of base algorithm. The results were clear. Among the four models, the self-training LightGBM achieved the most favorable overall balance between classification performance, spatial prediction efficiency, and computational efficiency, outperforming its rivals across the evaluation metrics used in the study.

The spatial results are perhaps the most striking. The geochemical anomalies recognized by the self-training LightGBM model cover only 2.74 percent of the entire study area, yet they encompass all of the known molybdenum deposits. In exploration terms, this is exactly what one wants: a small, high-confidence footprint that directs drilling and follow-up fieldwork toward the most promising ground. Anomalies that blanket huge swaths of terrain may technically capture the deposits, but they offer little practical guidance and inflate exploration costs.

Equally important is how the predicted anomalies align with the region’s geology. The team found close spatial consistency between the recognized anomalies and major faults, Mesozoic intermediate to acidic intrusions, and the contact zones where those intrusions meet the surrounding country rock. This makes geological sense, since such intrusion-country-rock contact zones are classic loci for ore-forming fluids, and faults often serve as conduits for mineralizing fluids. The fact that the model’s predictions converge on these structurally and magmatically favorable settings suggests that the algorithm is genuinely learning the fingerprints of mineralization rather than fitting statistical noise.

The implications extend beyond a single case study in Inner Mongolia. Mineral exploration increasingly relies on machine learning to sift through ever-larger geochemical and geophysical datasets, but the field has been hampered by the same recurring obstacles: complex geological settings, extreme class imbalance, and sparse labels. By combining a fast, nonlinear base learner with semi-supervised self-training and targeted oversampling, the new framework offers a template that could be adapted to other deposit types, other elements, and other survey scales. The authors note that the method provides a useful application approach for recognizing mineralization-caused anomalies from sparse labeled and vast unlabeled geochemical data collected in complex geological settings.

There are, of course, caveats. The framework was validated in one area, and its performance elsewhere will depend on data quality, the representativeness of the labeled samples, and the specifics of local geology. The researchers also acknowledge that no datasets were generated or analyzed beyond those used in the case study, meaning broader benchmarking remains future work. Still, the study, funded by the National Natural Science Foundation of China under grant number 42472361, represents a meaningful step toward smarter exploration. As the global demand for critical metals such as molybdenum continues to climb, tools that can squeeze more predictive power out of existing survey data, without demanding expensive new field campaigns, are likely to find an eager audience across the mining industry and academic geochemistry alike.

Subject of Research: A self-training LightGBM machine learning framework for recognizing mineralization-caused geochemical anomalies in stream sediment data.

Article Title: A self-training framework based on LightGBM for recognizing mineralization-caused geochemical anomalies

Article References: Zheng, C., & Chen, Y. (2026). A self-training framework based on LightGBM for recognizing mineralization-caused geochemical anomalies. Earth Science Informatics, 19(10), Article 181. https://doi.org/10.1007/s12145-026-02239-y

Image Credits: AI Generated

DOI: 10.1007/s12145-026-02239-y

Keywords: LightGBM, self-training, geochemical anomalies, mineral exploration, SMOTE, semi-supervised learning, molybdenum deposits, Inner Mongolia, stream sediment data, machine learning, mineral prospectivity, Jilin University

The tool, called HGCTM-S—short for hierarchical geological context topic model with sparse deviations—was developed by Muhammad Atif Bilal of Jilin University’s College of Geoexploration Science and Technology and Kateryna Hlyniana of Jilin University’s School of Mathematics and the Institute of Mathematics of the National Academy of Sciences of Ukraine. Rather than trying to force every deposit into a rigid classification box, the model treats each copper deposit as a mixture of recurring “assemblage modes,” statistical themes that capture groups of minerals that tend to appear together. The approach borrows its core logic from topic modeling, a technique originally devised to discover latent themes in large collections of text, and adapts it to the language of rocks: instead of words in documents, the model reads mineral families in deposit records.

The scale of the analysis is considerable. The researchers drew on the global copper deposit dataset, an open-source compilation covering 1,335 deposits, and organized 1,205 distinct mineral species into 35 geologically defined families. Grouping species into families was a deliberate choice to counteract sparse and inconsistent documentation, since individual rare minerals appear in too few records to support robust statistics on their own. The hierarchical structure of the model allows geological context—information about where and in what kind of rocks a deposit sits—to inform how the assemblage modes are expressed, while a sparse-deviation component captures localized anomalies that depart from the broader patterns.

When the model was fitted to the full dataset, it recovered seven assemblage modes, a number derived from the data itself rather than imposed in advance. The most geologically meaningful of these were validated against independent deposit type labels. A copper–molybdenum mixed mode emerged as strongly enriched in porphyry deposits, the giant intrusion-related systems that supply much of the world’s copper, while a nickel–cobalt–arsenic mode aligned closely with magmatic sulfide deposits, which form when sulfide liquids segregate from cooling magmas. These correspondences matter because the model was never told which minerals should characterize which deposit types; the associations emerged from the raw mineralogical records alone, and the deposit type labels served only as an independent check.

Just as telling were the modes that did not correspond neatly to genetic classes. Several of the remaining themes represented shared sulfide backgrounds common to many deposit styles, secondary overprints imposed by later weathering and alteration, or residual components that likely reflect the idiosyncrasies of the dataset rather than genuine ore-forming processes. The authors are explicit on this point: HGCTM-S is a tool for comparing overlapping mineral assemblage components and their regional geological associations, not a universal deposit classifier or a regional predictor. That restraint is rare and refreshing in a field where machine learning results are sometimes oversold as oracle-like prediction engines.

One of the study’s central questions concerned the relative influence of regional lithology versus broad geological age on mineral assemblage composition. Across alternative priors and multiple lithology proxies, the model consistently found that lithology-associated effective deviations were larger than age-associated deviations, suggesting that the kinds of host and country rocks surrounding a deposit shape its mineralogy more powerfully than the era in which it formed. The magnitude of this effect varied, however, and the lithology proxies proved unreliable at reproducing finer details such as deposit scale or specific host rock types—a reminder that coarse global datasets can constrain broad patterns but stumble at deposit-level resolution.

The team also probed the stability of their results. Progressive initialization, a strategy in which model fits are seeded sequentially to encourage convergence toward consistent solutions, improved the aggregate stability of the recovered topics. Yet geographic performance remained heterogeneous: in some regions the model’s topic assignments added value beyond what geological context alone could provide, while in others they did not reliably outperform a baseline built purely from contextual information. This heterogeneity is itself informative, pointing to regions where mineralogical records are rich and internally consistent, and others where documentation gaps or sampling biases dominate the signal.

The significance of the work extends beyond copper. Mineral informatics, the emerging discipline that applies data science to mineralogical databases, has matured rapidly over the past decade, with network analyses and association-mining studies revealing deep structure in how minerals co-occur through Earth history. HGCTM-S adds a probabilistic, context-aware ingredient to that toolkit, one that explicitly models uncertainty and mixture rather than demanding clean categories from messy reality. For exploration geologists, the framework offers a way to compare deposits in terms of their full assemblage fingerprints, potentially highlighting overlooked analogs and guiding targeting in data-rich terranes.

The study is also a candid case study in the limits of big-data geoscience. Public mineral databases are treasures, but they are treasures assembled by many hands over many decades, with unequal documentation, variable data quality, and overlapping mineralogical signals baked in. By quantifying where the model succeeds and where it falls short, the authors provide a template for honest evaluation that the broader community can adopt. As the energy transition drives unprecedented demand for copper—and for the cobalt, nickel, and molybdenum that often accompany it—tools that can faithfully extract geological meaning from imperfect global datasets will only grow in value. HGCTM-S does not replace the trained eye of the field geologist, but it gives that eye a new way of seeing the planet’s copper endowment all at once, one statistical theme at a time.

Subject of Research: Statistical topic modeling of mineral assemblages in global copper deposit databases

Article Title: HGCTM-S: Modeling Regional Lithology-Associated Mineral Assemblage Modes in Global Copper Deposits

Article References: Bilal, M. A., & Hlyniana, K. (2026). HGCTM-S: Modeling Regional Lithology-Associated Mineral Assemblage Modes in Global Copper Deposits. Natural Resources Research. https://doi.org/10.1007/s11053-026-10767-z

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10767-z

Keywords: copper deposits, mineral assemblages, topic modeling, machine learning, mineral informatics, porphyry deposits, magmatic sulfide deposits, lithology, geochemistry, mineral prospectivity, global copper dataset, Natural Resources Research