Hidden Geometry of Mahalanobis Distance Warps Geochemical Anomaly Maps

A statistical tool that geologists have trusted for nearly a century to flag hidden mineral deposits is quietly misleading them, and a team of Chinese researchers has now shown exactly why. In a study published in Natural Resources Research, Qinghao Zhang, Jilong Lu, Xinyun Zhao and colleagues at Jilin University, together with Mi Tian of the Chinese Academy of Geological Sciences and Liji Sun of the Hunan Center of Natural Resources Affairs, dissected the geometric failure modes of the Mahalanobis distance when a dataset secretly contains multiple mixed populations. Their diagnosis, and the two-stage remedy they built from it, offers a striking example of how thinking in terms of high-dimensional geometry can rescue a classical method from its own assumptions—and, in a case study from northeastern Hunan Province, helped capture 92 percent of known lithium deposits within just 7 percent of the surveyed area.

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