The core innovation is called vertically augmented three-dimensional Empirical Bayesian Kriging, or EBK3D. Classical kriging, a statistical interpolation technique borrowed from geostatistics, estimates values at unsampled locations by weighting nearby observations according to their spatial correlation. The Bayesian variant automates the difficult process of fitting the underlying semi-variogram, the mathematical function describing how similarity decays with distance, by simulating many plausible models and averaging the results. What the research team recognized, however, is that simply feeding three-dimensional coordinates into a standard kriging engine fails spectacularly for geotechnical data, because the statistics of soil change far more rapidly with depth than they do laterally.

To correct this anisotropy, the researchers introduced two empirically calibrated vertical adjustments. The first, an Elevation Inflation Factor, or EIF, rescales vertical separation distances relative to horizontal distances, so that a one-meter change in depth can be made statistically equivalent to tens or hundreds of meters of lateral separation. The second, a Vertical Trend Order, or VTO, captures systematic depth-related drift in soil properties, such as the steady increase in stiffness as overburden pressure compacts sediments. These two parameters are not chosen arbitrarily; the team calibrated them against a quality-controlled regional database compiled from field and laboratory investigation records, ensuring that the vertical geometry of the model reflects the actual behavior of the deposit rather than a convenient mathematical assumption.

Conventional two-dimensional and pseudo-three-dimensional approaches, the authors note, suffer from a cluster of related failures. They commonly assume planar stationarity, meaning they treat the statistical character of the soil as uniform across the mapped area, which is rarely true in layered alluvial environments. They discretize the subsurface into fixed layers, discarding information from observations that straddle layer boundaries. And they underrepresent directional covariance, ignoring the fact that correlations in soil properties often stretch further along a valley axis than across it. The consequences are over-smoothed predictions that wash out sharp geological contacts, dilution of lateral edges where data are sparse, and distorted volumetric estimates that can mislead foundation designers and zoning authorities alike.

The EBK3D framework attacks each of these weaknesses directly. Beyond the vertical adjustments, it couples directionally dependent semi-variogram parameterization, allowing the correlation structure to stretch and compress according to compass direction, with depth-aware spatial weighting that gives nearby observations in the vertical profile proportionally greater influence. It also employs covariance-optimized sector-based neighborhoods, dividing the space around each prediction point into sectors so that sample support is distributed evenly across directions rather than clustering where boreholes happen to be dense. Together, these components enable adaptive three-dimensional prediction and what the team calls digital stratigraphic reconstruction, essentially a probabilistic portrait of the stacked soil architecture at any point in the volume.

To demonstrate the power of the approach, the researchers generated computer-aided geotechnical soil maps, or cGSMs, for five critical engineering parameters: the Standard Penetration Test blow count, known as SPT-N, which indexes soil density and stiffness; lithology, the identity of the soil or rock type; shear strength, which governs bearing capacity and slope stability; the plasticity index, or PI, which characterizes how a fine-grained soil behaves when wet; and linear shrinkage, or LS, which reflects a clay’s tendency to crack and swell as moisture changes. These maps were then benchmarked against a suite of baseline two-dimensional and pseudo-three-dimensional models using identical data, providing a rigorous head-to-head comparison.

The results were dramatic. For the investigated dataset, a K-Bessel semi-variogram integrated with EBK3D proved the best performer, cutting predictive error by as much as 76 percent and raising the correlation between predicted and observed values by up to 40 percent. Focal standard deviation analysis, which measures how well a method preserves genuine localized variability rather than smearing it into regional averages, showed that EBK3D retained far more of the true spatial texture of the subsurface. Independent validation against field observations yielded correlation coefficients ranging from 0.87 to 0.98, an exceptionally tight band for spatial interpolation of geotechnical properties. Taken as a whole, the framework achieved improvements of up to 36 percent over the evaluated 2D and pseudo-3D competitors.

Just as important as accuracy is honesty about uncertainty, and here the study introduces a three-tier validation scheme spanning diagnostic, deterministic, and probabilistic metrics. Diagnostic checks confirm that the model’s internal error estimates are statistically well calibrated, comparing average standard errors against actual residuals. Deterministic metrics, including root mean square error, mean absolute error, Nash-Sutcliffe efficiency, Pearson correlation, and Kling-Gupta efficiency, quantify point-prediction skill. Probabilistic measures, such as the continuous ranked probability score, evaluate whether the full predictive distributions are trustworthy, which matters enormously when engineers must decide how much margin to build into a foundation design. This layered verification supports the practical applicability of the resulting three-dimensional cGSMs for preliminary geotechnical zoning and for targeting site investigations where they will yield the most information.

The implications reach well beyond one study region. Reliable regional subsurface characterization is a bottleneck for infrastructure planning worldwide, particularly in rapidly urbanizing areas where funding for dense drilling campaigns is scarce. By extracting maximal value from existing borehole archives and field investigation records, and by explicitly quantifying where predictions are confident and where they are not, uncertainty-aware 3D mapping can direct scarce investigation budgets to the locations of greatest risk and ignorance. It can also sharpen seismic site-response zonation, excavation design, and pile foundation assessment, all of which are acutely sensitive to vertical soil variability. The study’s authors, drawn from Tongji University, Quanzhou University of Information Engineering, Northern Border University, the University of Derby, and Khwaja Fareed University of Engineering and Information Technology, acknowledge support from the National Natural Science Foundation of China, the Shanghai Post-doctoral Excellence Program, the Natural Science Foundation of Fujian Province, and the Deanship of Scientific Research at Northern Border University. As cities dig deeper and climate change stresses the ground beneath our infrastructure, tools that reveal the third dimension of the subsurface, together with the confidence attached to every prediction, may become as indispensable to civil engineers as the maps they already hang on their walls. What began as a statistical refinement of kriging may end up changing how the ground itself is charted.

Subject of Research: Uncertainty-aware three-dimensional digital stratigraphic reconstruction of geotechnical soil variability using vertically augmented empirical Bayesian Kriging

Article Title: Uncertainty-aware digital stratigraphic reconstruction of geotechnical variability using vertically augmented empirical Bayesian Kriging

Article References: Ijaz, Z., Zhao, C., Ijaz, N., Hussain, W., Akbar, M., Rehman, Z. U., & Khalid, U. (2026). Uncertainty-aware digital stratigraphic reconstruction of geotechnical variability using vertically augmented empirical Bayesian Kriging. Earth Science Informatics, 19(11), Article 192. https://doi.org/10.1007/s12145-026-02243-2

Image Credits: AI Generated

DOI: 10.1007/s12145-026-02243-2

Keywords: geotechnical engineering, Empirical Bayesian Kriging, 3D subsurface mapping, spatial interpolation, semi-variogram, soil variability, digital stratigraphy, uncertainty quantification, infrastructure planning, geospatial modeling, SPT-N, shear strength