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

The research team, led by Mukesh P. Chaudhari of Gujarat University’s Department of Chemistry together with Ruchi Nair, Pratik Chavda, Dharmik Patel and Divya R. Mishra, combined systematic field sampling with Geographic Information System (GIS) mapping and multivariate statistics. Their goal was straightforward but ambitious: to measure the spatial distribution of the major soil fertility indicators across all fourteen talukas of Banaskantha and to translate those measurements into practical guidance for sustainable land management in a region where salinity, nutrient depletion and drought are intensifying pressures.

To capture the district’s variability, the researchers collected 46 geo-referenced composite soil samples from agricultural fields using a systematic 15 by 15 kilometer grid, an approach aligned with internationally accepted spatial sampling guidelines from the Food and Agriculture Organization and the USDA-NRCS Soil Survey Manual. At each sampling point, five sub-samples were taken in a zig-zag pattern within a 10 to 15 meter radius and homogenized into a single representative composite. All samples came from the 0 to 15 centimeter plough layer, the zone most responsive to nutrient availability and management. Each sample was then analyzed in the laboratory for soil moisture, pH, electrical conductivity (EC), organic carbon (OC), and available nitrogen, phosphorus and potassium.

Turning those numbers into pictures required spatial interpolation. Using ArcGIS version 10.8.1, the team applied the Inverse Distance Weighting (IDW) technique, a deterministic method that estimates unknown values from the weighted average of nearby sampling points, assuming that closer observations carry more influence than distant ones. The resulting raster layers, classified into concentration ranges and overlaid on the district boundary, produced district-scale thematic maps of every fertility indicator. The authors are candid that these maps are intended for regional visualization rather than precise field-level prediction, and that formal cross-validation statistics were not performed; they recommend denser sampling networks and geostatistical validation in future work.

The findings reveal dramatic heterogeneity. Soil moisture ranged from less than 1 percent to more than 20 percent, with the highest values in Suigam (24.74 percent), Vav (13.92 percent) and Lakhani (12.87 percent), where finer clay- and silt-rich soils retain water through greater surface area and capillary action, and where shallow groundwater or capillary rise may contribute. In contrast, the coarse sandy soils of Dhanera, Deesa, Kankrej and parts of Tharad and Palanpur drain rapidly and lose moisture to intense evapotranspiration, leaving almost nothing in reserve. Soil pH spanned from slightly acidic (6.34) in the forest-influenced taluka of Danta, where higher rainfall leaches basic cations, to strongly alkaline (8.65) in Vav, Suigam and parts of Lakhani, where aridity and carbonate-rich parent materials drive salt accumulation.

Electrical conductivity told perhaps the starkest story. While some soils in Deesa, Palanpur and Vadgam showed EC values as low as 0.01 dS/m, samples from Suigam and Lakhani exceeded 200 dS/m, indicating severe salinity likely driven by saline groundwater, an arid climate and evaporitic concentration of salts at the surface. Such salinity lowers osmotic potential, hampers water uptake by plants and damages soil structure. The study also identified a persistent district-wide phosphorus deficit: available phosphorus ranged from just 0.25 to 7.97 mg/kg, with most talukas below 3 mg/kg, because alkaline conditions cause phosphorus to precipitate as insoluble calcium phosphates. Available potassium, by contrast, swung from 49 to 884 mg/kg, with very high values in Bhabhar, Deesa and Danta, likely reflecting mica-rich parent material or heavy fertilizer input. Available nitrogen ranged even more widely, from 78 to 2,100 mg/kg, with the highest levels in Bhabhar, Kankrej and Lakhani, probably tied to farmyard manure, nitrogen fertilizers and clay that retains ammonium.

Organic carbon emerged as a central thread running through the district’s fertility story. Levels ranged from a very low 0.10 percent to a high 3.50 percent, with the richest soils found in Danta’s forest-edge environments, where cooler microclimates, better moisture and greater biomass return slow decomposition and build humus. The intensively farmed alluvial plains of Deesa, Dhanera, Palanpur and Kankrej showed very low organic carbon, a consequence of sandy textures, rapid oxidation under high temperatures and continuous cropping with minimal residue return. Correlation analysis reinforced carbon’s pivotal role: organic carbon was positively associated with available phosphorus (r = 0.499) and available potassium (r = 0.444), while the strongest relationship in the entire matrix linked available phosphorus and potassium (r = 0.778), suggesting shared parent materials or similar fertilization histories.

Principal Component Analysis then distilled the district’s complexity into two dominant processes. Validated by a Kaiser-Meyer-Olkin statistic of 0.61 and a highly significant Bartlett’s test (p < 0.001), the PCA extracted two components with eigenvalues greater than 1 that together explained over 72 percent of the total variance. The first, accounting for 55.42 percent, loaded heavily on electrical conductivity, available potassium, available nitrogen and soil moisture, representing a salinity-nutrient enrichment factor typical of semi-arid regions where evapotranspiration exceeds precipitation and soluble ions accumulate at the surface. The second, explaining 17.32 percent, was dominated by pH, organic carbon and available phosphorus, capturing the organic matter-fertility relationship in which carbon-rich soils maintain more stable pH and better phosphorus availability.

The practical implications are as uneven as the soils themselves. The authors conclude that a single, uniform fertilizer recommendation is unlikely to optimize productivity across the district. Eastern and northeastern areas, including Danta, Bhabhar, Amirgadh and parts of Vav, show comparatively better soil health and should focus on conservation practices to maintain fertility. Nutrient-deficient zones in the central plains call for integrated nutrient management combining organic amendments, residue retention, biochar or farmyard manure with balanced fertilization. The severely saline talukas of Suigam and Lakhani require reclamation through gypsum application, improved drainage, optimized irrigation and salt-tolerant crops. In high-pH zones, phosphorus-solubilizing biofertilizers, split phosphorus doses and organic matter additions could unlock trapped nutrients, while high-nitrogen areas should be monitored to prevent nitrate leaching into groundwater.

The researchers emphasize that their nutrient management strategies remain preliminary decision-support recommendations. The agronomic effectiveness and economic feasibility of the proposals have not yet been tested through crop response experiments, fertilizer trials or cost-benefit analyses, and the study did not include indicators such as soil texture, cation exchange capacity, micronutrients or biological properties. Still, by fusing laboratory chemistry, GIS interpolation and multivariate statistics into a single framework, the study establishes an important baseline dataset for regional soil fertility assessment. In a state that contains some of India’s largest extents of saline and alkaline land, and a country where roughly 175 million hectares of agricultural area face soil-related constraints, showing precisely where a district’s soils are thriving and where they are failing may prove to be the first, indispensable step toward farming that fits the ground it stands on.

Subject of Research: Spatial variability of soil fertility indicators assessed with GIS and multivariate statistics for sustainable land management in a semi-arid Indian district

Article Title: Spatial variability of soil fertility indicators using GIS for sustainable land management in Banaskantha district western India

Article References: Chaudhari, M. P., Nair, R., Chavda, P., Patel, D., & Mishra, D. R. (2026). Spatial variability of soil fertility indicators using GIS for sustainable land management in Banaskantha district western India. Discover Soil, 3(1), Article 161. https://doi.org/10.1007/s44378-026-00320-y

Image Credits: AI Generated

DOI: 10.1007/s44378-026-00320-y

Keywords: soil fertility, spatial variability, GIS, Inverse Distance Weighting, soil salinity, organic carbon, soil pH, available nutrients, principal component analysis, precision agriculture, sustainable land management, semi-arid India

Zhejiang presents an ideal and demanding test case. This coastal province in southeastern China is home to roughly 66.7 million people and an economy exceeding 9 trillion CNY in 2024, yet nearly 75 percent of its land is covered by hills and mountains, with only about 20 percent plains and a sliver of rivers and lakes. Rapid urbanization compounds the challenge, since weather stations are sparse, unevenly distributed, and subject to relocations, instrumentation changes, and gaps in the historical record. Ground observations offer accuracy but limited spatial coverage; satellite land surface temperature products offer detail but generally begin only in the early 2000s and are vulnerable to cloud cover, terrain shading, and atmospheric interference; reanalysis products such as ERA5 provide continuity but at coarse spatial resolutions, typically around 0.25 degrees or coarser, far too blunt to resolve the fine thermal texture of mountainous landscapes.

The team’s solution is an elegant two-part construction they call a spatial-background residual framework. First, they built a high-resolution climatological baseline from an existing hourly, 1-kilometer gridded temperature dataset covering 2008 to 2018, which had itself been developed using the INCA data-fusion framework with reanalysis fields and dense automatic weather station observations. This baseline serves purely as a spatial background, encoding how temperature varies across the terrain on each calendar day of the year. Second, daily temperature residuals—the departures of each station observation from that climatological expectation—were calculated for every meteorological station across the full 1961 to 2020 period. Because these residuals are computed directly from observed temperatures, they inherently preserve the long-term warming trend and interannual variability, while the baseline contributes the terrain-driven spatial detail. Summing the interpolated residual field with the baseline yields the finished reconstruction.

A critical methodological question was how best to interpolate those daily residuals across space. The researchers evaluated three schemes representing different levels of topographic correction and complexity: plain inverse distance weighting, or IDW, which relies only on spatial proximity; a lapse-rate-adjusted version of IDW, in which station temperatures are first corrected to grid-cell elevation using a fixed adiabatic lapse rate of 6.0 degrees Celsius per kilometer; and a multiple linear regression incorporating longitude, latitude, and elevation as predictors. Using leave-one-out cross-validation, in which each station is successively withheld and predicted from the others, the team assessed performance with mean absolute error, root-mean-square error, and the coefficient of determination across decades, seasons, and elevation zones.

The verdict was striking: the simplest method won. IDW without any lapse-rate correction consistently delivered the lowest errors and highest skill, achieving the best performance at roughly 68 percent of stations for daily minimum temperature, 70.7 percent for mean temperature, and 74.7 percent for maximum temperature. In a representative example from 1971 to 1980, IDW reconstructed minimum temperatures with a mean absolute error of just 0.73 degrees Celsius and an R-squared of 0.987, comfortably beating both rivals. The reason lies in the behavior of the lapse rate itself. Analysis of 60 years of observations revealed that near-surface temperature lapse rates in Zhejiang are strongly non-stationary: they peak in summer, with minimum-temperature lapse rates exceeding 7.0 degrees Celsius per kilometer in mountainous areas during July and August, yet collapse toward zero or even turn negative in winter lowlands, where temperature inversions prevail. Applying a fixed correction therefore risks systematic, elevation-related biases—a caution with implications well beyond Zhejiang.

The errors that do remain follow clear and intelligible patterns. Reconstruction accuracy improved steadily from the 1960s onward as station density grew, and summer months outperformed winter months because spatial temperature gradients are weaker in warm weather. Low-elevation areas below 400 meters consistently yielded smaller errors than high-elevation zones, where complex terrain and sparse instrumentation conspire against interpolation. Spatially, larger uncertainties cluster in the mountainous southwest, including parts of Lishui and western Wenzhou, while the plains around Hangzhou, Shaoxing, and Jinhua show excellent agreement, with most stations achieving R-squared values above 0.90 and many above 0.95. Across the entire 60-year span, the annual mean error for all three temperature variables stayed within plus or minus 0.1 degrees Celsius, with no systematic drift across decades—a testament to the temporal stability of the method.

Perhaps the most compelling validation came from real disasters. The team tested ZJ-DAT against two extreme events from 2007, using more than a thousand automatic weather stations as independent ground truth while deliberately excluding any stations that had contributed to the reconstruction. During the cold wave of 4 to 9 March 2007, ZJ-DAT tracked the south-to-north advance of the cold air, accurately reproducing the observed cold centers around Lishui, with R-squared values of 0.65 to 0.79 and root-mean-square errors of 1.01 to 1.67 degrees Celsius. By comparison, the CDAT national dataset managed only moderate agreement, while ERA5 performed poorly, with near-zero or negative correlations and errors approaching 3 degrees Celsius. The heat wave of 30 June to 10 July 2007 told the same story: ZJ-DAT best resolved the core hot zones above 37 degrees Celsius over Jinhua, Shaoxing, and Ningbo and the inland-coastal thermal contrast, while CDAT smoothed away local extremes and ERA5 drifted with warm biases and excessive homogenization. Case studies of cold and heat events in January and July 2020 at four environmentally distinct stations—an island, a mountain site, and two plain stations—further confirmed the reconstruction’s fidelity, with discrepancies generally under 2 degrees Celsius.

Beyond validation, the dataset enabled a first-of-its-kind hazard assessment for the province. Using Gumbel distribution analysis of return periods, the researchers mapped the intensity of extreme cold and heat expected at 5-, 20-, and 50-year recurrence intervals. The results expose stark geographic contrasts in climate risk. Extreme low-temperature hazards concentrate in the northwestern inland regions, where 50-year minimum temperatures plunge below minus 15 degrees Celsius, while the southeastern coast stays comparatively mild. Extreme heat hazards show the opposite pattern, dominated by low-altitude basins in central and northern Zhejiang, where 50-year maximum temperatures climb above 43 to 44 degrees Celsius—figures that carry sobering weight given projections of accelerating heatwave duration under global warming. These maps, grounded in kilometer-scale temperature data rather than coarse reanalysis, offer planners a far sharper picture of where adaptation investments are most needed.

The authors are candid about limitations. Anchoring the reconstruction to a climatology drawn from 2008 to 2018 means the reference field does not reflect the climate state of earlier decades, though because it functions only as a spatial scaffold while temporal signals come from station residuals, warming trends and variability remain intact. The team suggests that future refinements could employ temporally adaptive reference fields. The broader significance, however, is clear: ZJ-DAT demonstrates that a simple, computationally efficient interpolation of station anomalies, layered onto a modern high-resolution climatology, can outperform more elaborate schemes in complex terrain—provided the scheme is chosen with local lapse-rate physics in mind. The framework, and the publicly available dataset released through Zenodo, is designed to be transferable to other topographically complex, observation-limited regions, offering a practical foundation for climate monitoring, extreme-event risk assessment, and adaptation planning as the planet continues to warm.

Subject of Research: High-resolution daily near-surface air temperature reconstruction for Zhejiang Province, China, from 1961 to 2020 using statistical residual interpolation

Article Title: A high-resolution daily temperature reconstruction for Zhejiang Province during 1961–2020 using statistical residual interpolation

Article References: Li, Y., Guo, H., Dong, M., Wu, J., Deng, F., Chen, Y., & Chen, F. (2026). A high-resolution daily temperature reconstruction for Zhejiang Province during 1961–2020 using statistical residual interpolation. Theoretical and Applied Climatology, 157(10), Article 632. https://doi.org/10.1007/s00704-026-06521-3

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06521-3

Keywords: Zhejiang Province, air temperature reconstruction, inverse distance weighting, spatial interpolation, temperature lapse rate, extreme temperature events, heat wave, cold wave, climate monitoring, hazard assessment, ERA5, high-resolution gridded dataset