The research, published in Discover Geoscience by Dionisio Uendro Carlos, César Augusto Moreira, and Guilherme Libardi Maestrelo, introduces a moving-window ordinary kriging (MW-OK) algorithm designed specifically for airborne and ground radiometric data. Conventional interpolators such as inverse distance weighting and minimum curvature are fast and widely used, but they share a fundamental weakness: they do not explicitly account for the local spatial structure of the data. In practice, that means they can oversmooth short-wavelength variability, blurring exactly the sharp radiometric gradients that mark boundaries between rock units or zones of hydrothermal alteration.

Kriging, by contrast, is a geostatistical method built on a statistical model of how a property varies in space. Because radioactive decay is inherently a stochastic process, the authors argue, an interpolator that respects this randomness is desirable. Yet standard kriging typically relies on a single global variogram, the mathematical function describing spatial correlation, assumed to hold across the entire survey area. Radiometric variables rarely obey that assumption at regional scale. Lithological boundaries, weathering domains, and survey-related effects all conspire to make the spatial structure change from place to place, so a global model can over-smooth local features or produce biased estimates where the structure shifts.

MW-OK sidesteps this problem by abandoning the global model altogether. For every prediction node on the output grid, the algorithm defines a circular window of investigation, in this case with a radius of 750 meters, roughly 1.5 times the 500-meter flight-line spacing, so that each window intersects two adjacent flight lines. Within that window, the method assumes second-order stationarity holds only approximately and locally. An empirical semivariogram is computed from all unique pairs of samples in the neighborhood, binned into eight lag classes, with any bin containing fewer than three valid pairs rejected. An exponential model is then fitted to each local semivariogram by bounded non-linear least squares, yielding a nugget effect, partial sill, and range parameter tailored to that specific patch of ground.

The computational scale is striking. The study area, a 40 by 40 kilometer subset of the Sul-Riograndense Shield Aerogeophysical Project flown by the Geological Survey of Brazil in 2010, contains 50,374 measurement points. Because the algorithm estimates a local variogram for each measurement point, it effectively generates 50,374 variograms. Neighbor retrieval is accelerated with a k-dimensional tree using ball queries, and safeguards keep the local kriging systems well conditioned: a minimum of eight neighbors, a maximum of thirty, removal of duplicate coordinates to prevent singular matrices, and a tiny ridge term added to the kriging matrix diagonal for numerical robustness. If a non-linear variogram fit fails, the workflow falls back on physically sensible defaults derived from the local data variance and median nearest-neighbor spacing.

The fitted local variograms themselves carry geological information. For potassium, the nugget effect was essentially zero at the scale of the data support, meaning microscale variability and measurement noise are minor compared with the spatially structured component, and correlation decayed over a practical range of about 1.3 kilometers. Equivalent uranium showed an even longer continuity, around 1.7 kilometers, consistent with smoother spatial behavior, while equivalent thorium displayed a larger partial sill, indicating higher structured variance, with a moderate range near 1.5 kilometers. These parameters confirm that, within each window, the radiometric variability is predominantly spatially structured rather than random noise.

When the researchers applied MW-OK to the survey data, the resulting potassium, uranium, and thorium maps reproduced the broad patterns that geologists expect. Potassium enrichment above 4 percent traced the São Sepé Granite in the north, elevated values of roughly 4 percent marked the Ramada Granite in the southwest, and low concentrations coincided with the metatuff-dominated Bossoroca Complex and the sandstones of the Maricá Group. Thorium highs of about 21 parts per million again outlined the São Sepé Granite, with additional anomalies over the Cerro da Cria Granite and nearby rhyolites. But the kriged surfaces went further than the conventional minimum-curvature grids: they better represented local gradient variations near lithological contacts and revealed internal heterogeneity within the felsic bodies that the smoother deterministic algorithm muted. Notably, three east-west linear artifacts aligned with the survey tie lines, visible in the minimum-curvature uranium map, were suppressed in the kriged surface.

A key advantage of kriging is that it also delivers an estimate of its own uncertainty. The kriging variance maps for all three radioelements bore the clear imprint of the acquisition geometry: uncertainty was minimized along and near the north-south flight lines and increased in the gaps between them, telling interpreters exactly where the estimates are most reliable. To quantify how closely MW-OK matched the industry-standard approach, the team performed pixel-wise comparisons with the minimum-curvature grids. The agreement was remarkably strong, with coefficients of determination of 0.9938 for potassium, 0.9727 for equivalent uranium, and 0.9946 for equivalent thorium, and residual histograms tightly centered near zero. Ordinary least-squares fits revealed a subtle but systematic amplitude scaling, with slopes slightly below unity, suggesting MW-OK modestly attenuates the highest values while remaining close to the identity overall. The authors are careful to frame this as evidence of methodological consistency rather than broad superiority: the two methods agree at the domain scale, and their differences emerge primarily at local scales.

Beyond the three radioelement grids, the workflow generates the interpretive products that exploration geologists actually use. It computes standard radiometric ratios (K/eTh, K/eU, eU/eTh) and the F-parameter, a widely used proxy for hydrothermal alteration that relates potassium abundance to the uranium-to-thorium ratio. It also produces ternary radioelement maps in three complementary forms: the classical K-eU-eTh display, normalized ternaries, and compositional-data transformations using additive, centered, and isometric log-ratios that respect the closure constraint and reduce spurious correlation when treating the three elements as parts of a composition. In the study area, a historically significant gold district in southern Brazil where mineralization is hydrothermal and structurally controlled in quartz veins hosted by schists and granites, these products delineated the main lithological contacts and showed spatial correspondence with known gold-mineralized zones and areas of gold enrichment identified by stream-sediment prospecting.

Perhaps most importantly for the wider community, the entire workflow is open source and exports standard deliverables, including GXF grids, GeoTIFF rasters, KMZ overlays, and PNG figures, with proper spatial metadata, making it directly interoperable with both open-source tools such as Python, QGIS, and Google Earth and commercial platforms like Oasis Montaj, Intrepid Geophysics, and ArcGIS. The authors acknowledge the method’s limitations: effective spatial resolution remains constrained by flight-line spacing and sampling interval, computational cost grows with finer grids and larger neighborhoods, and sparsely sampled sectors can yield under-populated windows with unstable variograms, which the workflow mitigates through minimum-neighbor thresholds and k-nearest-neighbor fallbacks. Looking ahead, they plan to incorporate co-regionalized cokriging constrained by lithology and structure, and deep-kriging approaches, within the same moving circular-window framework. For now, MW-OK stands as a robust, locally adaptive alternative that offers incremental rather than transformative gains over deterministic gridding, but gains that could make the difference between a blurred contact and a clearly resolved exploration target.

Subject of Research: A locally adaptive geostatistical interpolation method for processing airborne and ground gamma-ray spectrometric data

Article Title: A moving-window ordinary kriging approach for airborne and ground gamma-ray spectrometric data processing

Article References: Carlos, D. U., Moreira, C. A., & Maestrelo, G. L. (2026). A moving-window ordinary kriging approach for airborne and ground gamma-ray spectrometric data processing. Discover Geoscience, 4(1), Article 289. https://doi.org/10.1007/s44288-026-00674-4

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00674-4

Keywords: gamma-ray spectrometry, kriging, geostatistics, radiometric mapping, mineral exploration, semivariogram, interpolation, gold prospecting, geological mapping, Brazil, open-source software, compositional data

The research, published in Discover Geoscience by Khirsagar Patel, Prasoon Soni, and Pushpraj Singh of Guru Ghasidas Central University in Bilaspur, drew on groundwater level records from 69 observation wells maintained by the Central Ground Water Board, spanning the years 2001 to 2025. Rather than treating the water table as a single static surface, the team divided the record into three hydrological seasons: the rainy monsoon months from June to September, the winter period from October to January, and the summer months from February to May. For each well and each season, they computed average depths to water in meters below ground level, then screened the dataset for completeness, duplicates, and spatial validity before any mapping began. All well locations were projected into a common coordinate system, and the interpolated surfaces were generated at a uniform 30-meter resolution, providing a consistent framework across every test.

The core of the study was a comparison between two fundamentally different families of interpolation. Deterministic methods, including inverse distance weighting and regularized spline, estimate values at unsampled locations using purely mathematical rules. Inverse distance weighting assumes that nearby wells exert more influence than distant ones, with the influence decaying according to a distance power parameter, while spline fits a smooth curved surface that passes exactly through every measured point and minimizes overall curvature. These methods are simple and fast, but they carry no statistical model of how groundwater actually varies across the landscape, and they offer no way to quantify how confident one should be in a prediction between wells.

Geostatistical methods take a different route. Ordinary kriging, the technique tested here, explicitly models the spatial autocorrelation of the data using a semivariogram, a function that describes how the variance between paired observations grows with their separation distance. From the experimental semivariogram, computed from the well data, the researchers fitted four theoretical models: Circular, Spherical, Exponential, and Gaussian. Each model captures a different character of spatial continuity. The Spherical model, a staple of hydrogeological work, shows spatial dependence rising with distance and then flattening at a defined range. The Exponential model rises quickly at short distances and approaches its sill asymptotically, suiting datasets with short-range variability. The Gaussian model produces very smooth surfaces, while the Circular model describes moderate continuity with an abrupt change in behavior at the range. Kriging then uses the fitted model to assign optimal weights to surrounding wells, and crucially, it also yields an estimate of prediction uncertainty at every location.

The semivariogram analysis itself revealed a great deal about how the Arpa Basin’s aquifer behaves. Across the rainy season, semivariance increased steadily with distance, indicating strong spatial dependence, with a low nugget effect suggesting minimal measurement error. Groundwater levels remained spatially correlated over distances of roughly 20 to 35 kilometers, meaning the water table varies in a structured, basin-scale way rather than randomly from well to well. In winter, the fitted models showed ranges of about 21.9 to 33.5 kilometers, with the Spherical and Exponential models extending furthest. Summer, when reduced recharge and heavy pumping stress the aquifer, showed weaker but still statistically significant spatial dependence, with model ranges between roughly 18.4 and 31.9 kilometers. These seasonal shifts matter because they demonstrate that the spatial structure of the water table is not fixed: it changes as recharge, abstraction, and hydrogeological conditions change through the year.

To judge which method and model performed best, the team used cross-validation, a demanding test in which each observation is removed in turn and predicted from the remaining wells, allowing predicted and observed values to be compared across the entire dataset. Five statistics were examined together: Mean Error, which measures bias; Root Mean Square Error, which measures overall prediction accuracy; Average Standard Error and Mean Standardized Error, which assess whether the model’s own uncertainty estimates are honest; and Root Mean Square Standardized Error, which should ideally sit close to one. The ideal model would show a mean error near zero, the lowest possible RMSE, an average standard error close to the RMSE, and a standardized error near one. No single candidate achieved all of these targets simultaneously, which is itself an instructive finding.

The Exponential semivariogram consistently delivered the lowest prediction errors. In the rainy season it produced an RMSE of 3.787 meters and a mean error of just 0.0282; in winter, an RMSE of 3.9939 meters; and in summer, a mean error of 0.0486 with an RMSE of 4.3308 meters. Yet when the researchers turned to the uncertainty measures, the picture shifted with the seasons. In the rainy season, the Spherical model’s root mean square standardized error of 0.9953 came closest to the ideal value of one, with good agreement between its average standard error and RMSE. In winter, it was the Circular model that best calibrated its uncertainty, with a standardized error of 0.9968 and closely matched error statistics. In summer, the Gaussian model took that role. In other words, the model that predicts most accurately is not always the model that most honestly reports how sure it is, and the balance between the two shifts from season to season.

The mapped surfaces told a consistent visual story that reinforced the statistics. Ordinary kriging with the Spherical semivariogram produced the most realistic and coherent groundwater surfaces in all three seasons, with smooth transitions and clearly defined zones of shallow and deep water tables. The deterministic methods fared less well: inverse distance weighting generated localized exaggerations and bull’s-eye patterns around observation wells, while spline produced surfaces that were overly smooth and prone to overshooting. The seasonal maps also exposed a pattern of concern for water managers, with the deepest water tables concentrated in the southern and central portions of the basin, particularly during summer, signaling significant groundwater stress in the dry months in those areas.

The broader lesson of the study is that there is no universally superior interpolation model, and that model selection should rest on the combined interpretation of prediction accuracy and uncertainty measures rather than on any single statistic. For the Arpa River Basin, the Spherical semivariogram within ordinary kriging emerged as the most balanced overall choice for seasonal groundwater mapping, while the Exponential model offers an alternative where raw prediction accuracy is the priority. The geostatistical approach outperformed the deterministic alternatives in reliability and spatial consistency across every season tested. The resulting seasonal groundwater level maps provide a scientific foundation for monitoring network design, watershed planning, and sustainable abstraction policies in the basin. The authors suggest that future gains could come from incorporating additional hydrological variables, densifying the monitoring network, and applying advanced geostatistical or machine learning techniques, steps that would be especially valuable as climate variability and growing demand place ever greater pressure on Central India’s hidden water reserves.

Subject of Research: Comparative evaluation of geostatistical and deterministic interpolation methods for seasonal groundwater level mapping in the Arpa River Basin, Central India

Article Title: Comparative assessment of geostatistical and deterministic interpolation methods for seasonal groundwater level mapping in the Arpa River Basin, Central India

Article References: Patel, K., Soni, P., & Singh, P. (2026). Comparative assessment of geostatistical and deterministic interpolation methods for seasonal groundwater level mapping in the Arpa River Basin, Central India. Discover Geoscience, 4(1), Article 377. https://doi.org/10.1007/s44288-026-00736-7

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00736-7

Keywords: groundwater, spatial interpolation, ordinary kriging, semivariogram, inverse distance weighting, spline, ArcGIS, Arpa River Basin, Central India, hydrogeology, seasonal variability, cross-validation