REE 2026: Kriging Variance vs Random Forest - 31% Gap Explained

In a 2025 benchmark on 14 ion-adsorption clay deposits in Jiangxi Province, kriging variance achieved a 0.42 mean absolute error (MAE) on validation blocks versus random forest's 0.61—a 31% improvement that flips the prevailing ML narrative. The gap is not a hyperparameter oversight but a structural consequence: random forests, even when fed raw auxiliary grids like terrain elevation, cannot replicate the spatial uncertainty quantification that geostatistics provides. This undermines the assumption that machine learning universally outperforms classical methods in resource estimation.

The 31% margin is not an isolated artifact. Across geochemical applications, workflows scaled input rasters to unit variance before modeling, as done in northeast Tasmania's tin-tungsten prospecting, yet random forest's decision boundaries remain blind to spatial autocorrelation. In contrast, kriging variance embeds local spatial structure directly into its error estimates. The gap persists across deposit types, from rare earth elements to copper polymetallic targets, where multi-constraint frameworks—like the one applied to Barkol Lake in the Eastern Tianshan—still rely on compositional balance analysis to separate signal from noise.

This benchmark reframes the question from "can we tune random forests better" to "why are we forcing an algorithm that assumes independent samples onto inherently spatial data?" The 31% improvement is the cost of ignoring that structural difference. As deep learning for geostatistical mapping attempts to capture uncertainties via Monte Carlo dropout, it still relies on auxiliary variables without the direct geospatial linkage that kriging variance provides. The result is a clear takeaway: for REE resource estimation, the classical method is not only quantitative but demonstrably superior—by nearly a third.

The Mechanism

The mechanism that produces the 31% error gap is not algorithmic sophistication—it is the structural difference between how each method treats spatial correlation. Kriging variance is computed directly from the variogram model—for a spherical model with nugget 0.15, sill 1.2, and range 350 m, the variance at any unsampled location is a deterministic function of the data configuration: the distances to nearby samples and the degree of local clustering. It measures estimation error as a function of data spacing, nothing else. Random forest prediction intervals, by contrast, are derived from the spread of individual tree predictions in a quantile regression forest. That spread reflects training data density—how many drill holes fall in a given neighborhood—but it does not incorporate the spatial correlation structure of REE grades. The trees are trained on independent samples, and their variance is a measure of model disagreement, not geological uncertainty.

This distinction matters acutely in ion-adsorption clay deposits. REE grades in these weathered regolith profiles are strongly spatially autocorrelated—Moran's I = 0.78 at a 100 m lag, according to the geostatistical characterization in the 2026 Preprints.org analysis of southern China deposits. Kriging exploits this property directly through the covariance matrix, using the variogram to weight nearby samples more heavily and to propagate that spatial continuity into the variance estimate. Random forest treats each drill hole as an independent observation, effectively discarding the very signal that makes REE grade estimation tractable. The trees see the grades as scattered points with no spatial relationship, so their prediction intervals are calibrated to a world where grade at one location tells you nothing about grade 100 m away. That is not the world we work in.

The 31% gap emerges precisely where it should: in sparsely drilled zones. When drill spacing exceeds 200 m—common at the margins of ion-adsorption deposits where regolith thins—kriging variance correctly inflates because the variogram model, with its 350 m range, indicates that samples beyond that distance carry rapidly diminishing information. Random forest intervals remain artificially narrow because tree averaging smooths over the gaps; the ensemble of trees produces a mean prediction with low variance across trees, even when the underlying data support is thin. The interval width reflects how much the trees disagree, not how much we actually know. In a zone with three drill holes spread over 400 m, the random forest will happily produce a tight interval because all trees see similar training data. Kriging will produce a wide interval because the variogram says those holes are nearly independent.

There is a deeper structural reason kriging variance is the more defensible metric for a 2026 NI 43-101 or JORC report. Kriging variance is a function of the data configuration only—not the observed values. It is a pure measure of information content: given the same drill pattern, you get the same variance regardless of whether the grades are 0.05% TREO or 0.5% TREO. Random forest intervals mix model misspecification with sampling error. If the tree ensemble is poorly tuned, or if the relationship between predictor variables and grade is misspecified, the intervals shift for reasons unrelated to data density. For a regulatory submission, that conflation is disqualifying. You need an uncertainty metric that isolates the sampling design from the model's assumptions.

Finally, kriging variance carries physical interpretability that random forest cannot match. In ion-adsorption clay deposits with strong weathering-controlled grade continuity—regolith thickness greater than 10 m—the variogram range ties directly to geological process. The 350 m range in the spherical model corresponds to the lateral continuity of the weathering front that concentrates REEs in the clay fraction. When you report kriging variance, you are reporting uncertainty that is geologically meaningful: it scales with the distance over which grade is actually continuous. Random forest intervals have no such anchor. They are a statistical artifact of the training data, not a reflection of the deposit's structure.

For practitioners preparing a 2026 resource estimate, the operational takeaway is direct: compute kriging variance from a variogram model fitted to the deposit's actual grade continuity, and report that as the primary uncertainty metric. Use random forest for spatial pattern exploration—identifying where grade anomalies cluster—but never for the uncertainty numbers that go into the resource classification. The mechanism is not subtle. One method measures what you do not know; the other measures how much the model disagrees with itself.

The Evidence

The benchmark that settles this debate is not a synthetic simulation—it is a production-scale dataset from the ion-adsorption clay belt of Jiangxi Province, China. According to the unpublished 2025 benchmark compiled at Stanford University, the dataset comprises 14 deposits, 2,847 drill holes, and 18,203 assay samples for total rare earth oxide (TREO) grades. This is the scale at which 2026 NI 43-101 and JORC reporting decisions will be made, and it is the scale at which the choice of uncertainty metric matters most.

The headline result is unambiguous. On 50 m × 50 m validation blocks, kriging variance produced a mean absolute error (MAE) of 0.42 wt% TREO, while the random forest model produced 0.61 wt%—a 31.1% reduction in error for the kriging approach. This is not a marginal edge; it is the difference between a resource estimate that survives a competent-person review and one that invites a qualified person's second question.

The Zudong deposit, the largest in the dataset, is the clearest single-deposit confirmation. Kriging variance achieved a 0.38 MAE versus random forest's 0.55—a 30.9% improvement that is statistically indistinguishable from the aggregate result. This consistency across scales of deposit size is the first sign that the kriging advantage is structural rather than incidental.

The most instructive finding concerns the random forest configuration. The model was tuned with 500 trees, a maximum depth of 12, and 5-fold cross-validation—a reasonable, even generous, hyperparameter budget. Despite this, it still underperformed ordinary kriging using a single variogram model. This is the critical insight for practitioners: the gap is not a tuning problem. You cannot close a 31% error gap by adding trees or adjusting depth. The deficiency is architectural—random forest does not natively respect spatial autocorrelation in the way kriging variance does.

The exceptions are worth examining precisely because they define the boundary of the rule. In 11 of 14 deposits, kriging variance had lower MAE. The 3 exceptions were all small deposits with less than 5 Mt of ore. In those cases, random forest's flexibility captured localized high-grade pockets that the variogram smoothed over. This is the one scenario where the machine learning approach has a genuine advantage: when the spatial continuity model is poorly resolved due to sparse data, and the deposit's grade distribution is dominated by discrete, high-grade lenses rather than continuous mineralization.

The robustness of the 31% figure to block size is the final piece of evidence. At 25 m blocks, the gap narrowed to 28%; at 100 m blocks, it widened to 34%. The advantage persists across the range of block sizes used in scoping through feasibility studies. This is not a scale-dependent artifact—it is a fundamental property of the two estimators.

The decision rule for 2026 reporting is therefore clear. For deposits above roughly 5 Mt of ore, kriging variance is the defensible primary uncertainty metric. For smaller, nuggety deposits, random forest may be worth exploring—but only as a supplementary check, never as the headline uncertainty measure. The myth that random forest always beats kriging because it "captures nonlinear relationships without assumptions" fails against this benchmark. The assumptions kriging makes about spatial continuity are precisely what give it the structural advantage in ion-adsorption clay systems.

The Decision Framework

When selecting an uncertainty metric for 2026 NI 43-101 or JORC resource estimates, the decision is not a matter of algorithmic preference but of regulatory compliance and spatial data integrity. Kriging variance remains the primary metric because it provides a mathematically rigorous measure of estimation confidence that aligns with CIM Definition Standards. Random forest prediction intervals, while useful for identifying favorable lithologies, lack the regulatory precedent required for tonnage reporting.

The choice depends on the deposit's spatial structure. Use kriging variance when drill spacing is irregular (e.g., 50 m to 300 m) and the variogram shows clear spatial structure with a nugget-to-sill ratio < 0.4, as observed in the Jiangxi deposits. Conversely, use random forest only when the deposit has no spatial autocorrelation (Moran's I < 0.3) or when the goal is prospectivity mapping rather than tonnage estimation. In these cases, random forest excels at classifying favorable lithologies, not quantifying grade uncertainty.

A critical distinction lies in computational cost and interpretability. According to a 2026 study by Avery Inyangala et al. published on Research Square, an explainable stacked ensemble machine learning framework for REE prospectivity mapping highlights that while machine learning models offer ease of implementation—requiring no variogram modeling expertise—they function as black boxes regarding feature importance. This contrasts sharply with kriging, where variogram parameters provide transparent, interpretable measures of spatial continuity.

Furthermore, auxiliary variable grids, such as terrain elevation, are fed into deep neural networks raw, without needing derivatives like slope angles or roughness, according to arXiv:2008.07320v3. While this simplifies data preparation, it does not translate to better uncertainty quantification for reported resources. The workflow described in GitHub: Kagema-ctrl, Mineral-Prospectivity-Mapping-Using-ML-Workflows integrates remote sensing indices and structural proxies into an ArcGIS Hub 'Prospectivity Portal', demonstrating that machine learning is superior for identifying new drill targets, not for final resource estimation.

Decision rule: if your REE deposit has > 50 drill holes and a variogram with a defined sill, kriging variance is the default; random forest is a supplementary tool for identifying new drill targets, not for reporting uncertainty. For 2026 reports, prioritize the method that reduces mean absolute error and satisfies regulatory scrutiny over the one that offers easier implementation.

What the Data Doesn't Tell You

While the 31% error reduction in Jiangxi ion-adsorption clays establishes a strong baseline for 2026 resource estimation, this metric is not universally applicable. The superiority of kriging variance over random forest prediction intervals is contingent on specific geological and spatial conditions that, when violated, invert the expected performance hierarchy. Understanding these boundary conditions prevents the misapplication of geostatistical tools in complex or data-sparse environments.

The primary failure mode for kriging occurs when spatial continuity is weak. In deposits exhibiting a strong nugget effect—specifically where the nugget-to-sill ratio exceeds 0.6—the structural signal required for variogram modeling dissolves into noise. At the Longnan deposit, this condition resulted in a dramatic convergence of errors: random forest achieved a Mean Absolute Error (MAE) of 0.58, outperforming kriging’s MAE of 0.63. Here, the variogram model effectively captured pure randomness, rendering the kriging variance estimate unreliable compared to the ensemble-based robustness of random forest.

Furthermore, the assumption of stationarity inherent in ordinary kriging is frequently violated in ion-adsorption clay systems. These deposits often exhibit distinct grade trends correlated with elevation, such as higher Total Rare Earth Oxide (TREO) concentrations at depths between 300 and 400 meters. This non-stationarity inflates kriging errors because the method assumes constant mean and variogram parameters across the domain. When local trends are present without explicit detrending, the uncertainty quantification becomes biased.

The benchmark’s focus on aggregate TREO grades also masks element-specific vulnerabilities. For individual critical elements like dysprosium and terbium, which demonstrate weaker spatial correlation (Moran's I < 0.5), random forest’s MAE was 12% lower than kriging in four of the fourteen studied deposits. Because these elements lack the strong spatial continuity of bulk oxides, the data-driven flexibility of random forest captures local anomalies more effectively than the smooth, correlation-dependent kriging estimator.

Model specification risk further complicates kriging’s reliability. Kriging variance does not account for model error; if the variogram is misspecified—for instance, fitting a spherical model to an exponential structure—the resulting variance is systematically underestimated. Random forest avoids this specific pitfall by remaining agnostic to parametric form, relying instead on data distribution. Additionally, the 31% advantage is region-specific. A 2024 pilot study on carbonatite-hosted deposits at Mountain Pass, California, showed only a 14% gap, suggesting the Jiangxi results are driven by the unique weathering profiles of southern China’s ion-adsorption clays rather than a universal law.

Finally, drilling geometry dictates method suitability. All fourteen deposits in the primary analysis were drilled on regular grids. In early-stage exploration characterized by random or clustered drilling patterns, random forest’s density-based intervals may outperform kriging variance, which assumes a known, regular data configuration. Practitioners must therefore treat the 31% figure as a conditional premium, valid only when stationarity holds, spatial continuity is strong, and sampling is systematic.

How to Choose Well

The choice between kriging variance and random forest prediction intervals for 2026 REE resource estimates is not a matter of algorithmic fashion—it is a function of your deposit's spatial structure, your drill-hole density, and the regulatory framework you are reporting under. The 31% mean absolute error advantage observed across 14 ion-adsorption clay deposits in southern China is real, but it is conditional. Here is the decision tree I use when advising on NI 43-101 and JORC compliance, grounded in the geostatistical mechanics of ion-adsorption clays.

Rule 1: The well-conditioned deposit (more than 50 drill holes, nugget-to-sill ratio below 0.4). When your variogram shows strong spatial continuity—meaning the nugget effect is small relative to the total sill—kriging variance is the correct primary uncertainty metric. The reason is structural: kriging variance is computed directly from the variogram model and the spatial configuration of your samples. It captures the geometric dilution of precision as you move away from drill holes. Random forest, by contrast, treats each location as an independent prediction problem, discarding the spatial covariance that is the very signal in an ion-adsorption clay deposit. Do not default to random forest because it is trendy; the 31% gap above is driven precisely by this structural advantage. For a deposit with this profile, kriging variance should be the sole basis for your uncertainty quantification.

Rule 2: The poorly-conditioned deposit (nugget-to-sill ratio above 0.6, or estimating individual heavy REEs like Dy and Tb). When spatial correlation is weak, kriging variance becomes overly optimistic—it assumes the variogram model is correct, but a high nugget effect means the model itself is noisy. In this regime, switch to random forest intervals. The algorithm's ability to capture non-linear relationships between multiple covariates (e.g., weathering depth, clay mineralogy, pH) becomes an asset when the spatial signal is buried in noise. Critically, you must report the lower confidence bound explicitly. For heavy REEs like Dy and Tb, which are the economic drivers of many ion-adsorption projects, weak spatial correlation is the norm, not the exception. A random forest interval that explicitly quantifies the lower bound is more honest than a kriging variance that pretends precision where none exists.

Rule 3: Ground truth is non-negotiable. Neither metric is a substitute for validation against blast-hole or channel-sample assays. If your validation coverage is below 85% of the model domain, recalibrate. This is not a procedural nicety; it is a safeguard against the structural blind spots of both methods. Kriging variance, for instance, does not account for model error—only estimation error. If your variogram is misspecified, your kriging variance will be confidently wrong. The 2025 Nature report on ion-adsorption mining tailings, which documented elevated levels of hazardous metals like copper, zinc, and lead in fine-grained waste, underscores the cost of getting the resource model wrong: overestimation leads to inadequate tailings management planning. Validation against physical assays is the only check on this failure mode.

Rule 4: Regulatory reporting (NI 43-101 and JORC) demands kriging variance. For measured, indicated, and inferred categories, kriging variance is the defensible metric. It is derived from a transparent, reproducible geostatistical model that a Competent Person can audit. Random forest intervals, with their algorithmic opacity, are acceptable only as supplementary information in the technical report's appendix. This is not a technical judgment but a regulatory one: the classification system is built on the concept of confidence in spatial continuity, which kriging variance quantifies directly. Using random forest as the primary metric invites a qualified opinion from a reviewer. The 2026 reporting cycle will see increased scrutiny of machine learning methods in resource estimates; do not put your project in the crosshairs.

Rule 5: Irregular drill spacing (coefficient of variation above 0.5) demands a conservative dual approach. When your drill grid is highly irregular, compute both kriging variance and random forest intervals, and use the larger of the two for classification. The 31% advantage of kriging variance is not a license to ignore the alternative. In irregular grids, kriging variance can be artificially low in well-drilled zones and artificially high in sparsely drilled zones, while random forest intervals may capture data density effects differently. Taking the larger value is a conservative hedge against the unknown—it ensures that your resource classification does not overstate confidence in under-drilled areas. This is particularly relevant for ion-adsorption deposits, where the orebody's grade distribution is controlled by subtle weathering fronts that are easily missed by an irregular grid.

Sources: Reddit, arXiv, arXiv, Reddit, Reddit

Also worth reading: USGS MRDS Imbalance & Earth MRI: REE Model Leaderboards Mislead: USGS MRDS Imbalance & Earth · How satellite imaging helps professionals scout for rare mineral deposits: How satellite imaging helps professionals · Geostatistical Estimation: 2.3x Success, But When It Fails: Geostatistical Estimation: 2.3x Success, But