Rainfall is the single most important input to any urban flood model, yet cities typically measure it at only a handful of points. Rain gauges record precipitation where they stand, but a storm’s intensity can vary dramatically from one neighborhood to the next. To bridge those gaps, hydrologists use spatial interpolation, a family of statistical techniques that estimate rainfall values at unmeasured locations based on nearby observations. The problem, as the new research demonstrates, is that different interpolation methods — and even different parameter settings within the same method — can produce noticeably different rainfall fields from the same gauge data, and those differences cascade directly into flood predictions.

The research team systematically compared three widely used interpolation techniques: Ordinary Kriging, Inverse Distance Weighting, and Spline interpolation. Each method was tested under three search radius settings, meaning the number of surrounding gauge stations considered when estimating rainfall at any given point. That produced nine distinct configurations, labeled K3, K6, and K9 for Kriging; I3, I6, and I9 for Inverse Distance Weighting; and S3, S6, and S9 for Spline, where the numeral indicates the search radius. Each interpolated rainfall dataset was then used to drive a coupled hydrologic–hydrodynamic model of urban inundation, allowing the researchers to trace exactly how interpolation choices propagate through the entire flood simulation chain.

One of the study’s most striking results concerns a subtle but consequential artifact: negative rainfall values. Some interpolation methods, particularly Spline and certain Kriging configurations, can mathematically produce rainfall estimates below zero at locations between gauges — a physical impossibility that nonetheless appears in the interpolated output. The researchers found that K3, I3, I6, and I9 successfully avoided generating these spurious negative values, while other configurations did not. Negative rainfall, even in small amounts, can distort the water balance in a flood model and introduce errors that are difficult to trace back to their source. The fact that method selection alone determines whether this artifact appears underscores how sensitive the modeling pipeline is to choices that many practitioners make almost routinely.

When the team ran the urban inundation simulations with each of the nine rainfall fields, the differences were substantial. At the same locations within the study area, simulated inundation varied by more than 10 percent depending solely on which interpolation method had been used to prepare the rainfall input. In practical terms, this means that two flood forecasters working with identical gauge observations could produce flood maps that differ meaningfully in extent and depth — not because of model physics or terrain data, but because of how they chose to interpolate rainfall between measurement points. For emergency managers deciding which neighborhoods to evacuate or where to position response resources, a 10 percent discrepancy in predicted flood depth is far from trivial.

To move beyond simple comparison and establish which method actually performs best, the researchers employed spatial cross-validation, a rigorous statistical technique that systematically holds out rain gauge stations from the interpolation process and then evaluates how well each method predicts rainfall at the excluded locations. Combined with the urban flood hydrologic–hydrodynamic coupled model, this approach allowed the team to assess not just how well each method reproduced point rainfall measurements, but how those differences translated into flood simulation accuracy across the urban terrain.

The verdict, at least for the type of urban environment studied, was clear. For cities with flat terrain and minimal elevation variation — conditions that characterize many rapidly urbanizing areas — the Inverse Distance Weighting method with a search radius of six (I6) emerged as the optimal spatial interpolation technique for rainfall. This result is notable because Inverse Distance Weighting is often considered a simpler, less sophisticated alternative to geostatistical methods like Kriging, which explicitly model the spatial correlation structure of rainfall. Yet in the flat urban setting examined here, the simpler method with the right search radius outperformed its more mathematically elaborate competitors, suggesting that the assumption that more complex methods are inherently superior does not always hold in practice.

The study’s findings carry particular weight in the context of China’s ongoing efforts to manage urban flooding. Many Chinese cities have experienced catastrophic inundation events in recent years, prompting national investment in so-called sponge city infrastructure designed to absorb and slow stormwater. Accurate flood simulation is essential both for designing such interventions and for issuing timely warnings when extreme rainfall strikes. If the rainfall field driving a simulation is distorted by an inappropriate interpolation choice, the resulting flood maps — and the decisions based on them — may be systematically wrong in ways that are invisible to the analysts using them.

The research also highlights a broader challenge in hydrological modeling: the chain of assumptions and processing steps that stand between raw observations and model predictions. Much attention in flood forecasting focuses on improving model physics, increasing computational resolution, or incorporating new data sources such as weather radar and satellite precipitation. This study reminds the community that upstream data preparation — something as seemingly mundane as choosing an interpolation algorithm and its search radius — can exert influence comparable to those more glamorous improvements. A model fed with a poorly interpolated rainfall field will produce misleading results no matter how sophisticated its equations or how fine its computational grid.

Looking forward, the authors identify several promising directions for improving rainfall spatial interpolation in urban flood applications. Multi-source data fusion, which combines rain gauge measurements with radar estimates and other observational streams, could reduce reliance on any single sparse network. Hybrid interpolation techniques that blend the strengths of different methods may offer better performance across varied terrain and storm conditions. Rainfall merging algorithms, which optimally combine information from multiple sources while accounting for their respective uncertainties, represent another frontier. As cities worldwide face intensifying extreme rainfall under a changing climate, the seemingly small decision of how to interpolate between rain gauges may prove to be one of the most consequential choices in the entire flood forecasting pipeline.

Subject of Research: The influence of rainfall spatial interpolation methods on numerical simulation of urban inundation

Article Title: Research On the Influence of Different Rainfall Spatial Interpolation Methods On Numerical Simulation of Urban Inundation

Article References: Xiao, S., Hou, J., Wang, T., Li, D., Zhou, X., Pan, X., Chen, G., Xu, L., & Fei, K. (2026). Research On the Influence of Different Rainfall Spatial Interpolation Methods On Numerical Simulation of Urban Inundation. Water Resources Management, 40(13), Article 539. https://doi.org/10.1007/s11269-025-04402-4

Image Credits: AI Generated

DOI: 10.1007/s11269-025-04402-4

Keywords: rainfall interpolation, urban flooding, numerical simulation, Ordinary Kriging, Inverse Distance Weighting, Spline interpolation, spatial cross-validation, hydrologic-hydrodynamic model, flood forecasting, rain gauge data, Water Resources Management, Xi'an

Reference evapotranspiration, abbreviated ETo, is the backbone of irrigation planning worldwide. It represents the rate at which water would be lost from an idealized, uniformly green surface with ample moisture, and it is the quantity that farmers, water managers, and hydrologists multiply by crop coefficients to estimate real irrigation demand. The internationally accepted method for computing it is the FAO Penman–Monteith equation, a physically based formula that combines temperature, humidity, wind, and solar radiation into a single estimate. Because the equation draws on so many inputs simultaneously, a long-standing puzzle has been which of those inputs matters most, and whether the answer shifts from place to place and decade to decade. That question has become urgent as climate change alters the individual variables at different rates.

Zarei and Mahmoudi assembled a dataset spanning 1967 to 2024, one of the longest windows ever applied to this problem in the Middle East. From the records of 34 synoptic stations, they extracted nine candidate drivers: maximum, minimum, and mean temperature; maximum, minimum, and mean relative humidity; wind speed; sunshine hours; and precipitation. For thirteen representative years distributed across the nearly six decades, they computed ETo at every station using the FAO Penman–Monteith equation and then interpolated the results into continuous spatial maps using the Inverse Distance Weighting method, a technique that estimates values at unsampled locations by weighting nearby observations according to their distance.

To move beyond station-by-station analysis, the researchers overlaid 2,000 randomly distributed points across Iran’s territory and extracted both climatic data and interpolated ETo values at each one. This created a vast statistical sample capturing the country’s extraordinary climatic diversity, from the humid Caspian coast to the arid central deserts and the mountainous northwest. The random-point approach meant that the subsequent modeling reflected the spatial character of the whole country rather than the arbitrary placement of weather stations, which tend to cluster in populated areas.

The heart of the study lies in its dual statistical framework, and this is where the methodological story becomes genuinely interesting for anyone who follows climate science. The first tool, the Corrplot approach, examines simple linear correlations between each climatic variable and ETo, identifying which variable shows the strongest individual association in each year. The second, the Lindeman, Merenda, and Gold method, usually abbreviated LMG, belongs to a family of techniques called variance decomposition. Rather than looking at each variable in isolation, LMG averages the contribution of each predictor across all possible orderings in which it could enter a regression model, thereby partitioning the explained variance fairly among predictors that are correlated with one another.

Why does that distinction matter so much? Climatic variables are notoriously interdependent. Maximum temperature and sunshine hours rise and fall together; humidity tracks inversely with temperature; precipitation suppresses both sunshine and temperature. In statistical language, this is multicollinearity, and it is the bane of attribution studies. A naive correlation analysis can assign importance to a variable that is merely riding on the coattails of another. The LMG method, by decomposing R-squared across all orderings of predictors, corrects for this shared variance and reveals which variable carries genuine, independent explanatory power. The contrast between the two methods in this study is a textbook demonstration of why the correction matters.

According to the Corrplot results, maximum temperature was the most frequently dominant variable, leading in 46.16 percent of the years examined, followed by sunshine hours at 23.08 percent, with precipitation and wind speed each dominant in 15.38 percent of years. That picture is broadly consistent with intuition in a warming country: hotter days, more evaporative demand. But when the LMG variance decomposition was applied, the hierarchy flipped dramatically. Sunshine hours emerged as the most frequently dominant factor in 92.31 percent of the study period, with maximum temperature leading in only 15.38 percent of cases. In other words, much of the apparent importance of temperature in the simple correlation analysis appears to reflect its correlation with sunshine duration rather than an independent physical influence of that magnitude.

The authors validated both the spatial interpolations and the statistical models rigorously. The quality of the IDW-generated maps was assessed using the Nash–Sutcliffe Efficiency coefficient, a standard hydrological goodness-of-fit measure, and Willmott’s Index of Agreement, which evaluates how closely predicted values track observed ones. Both metrics confirmed the reliability of the spatial products at the 1 percent significance level, and the LMG model’s performance was likewise statistically verified at the 1 percent level for every year investigated. These are demanding thresholds, and their satisfaction across nearly six decades of data lends considerable weight to the central finding.

The practical implications for Iran are substantial. The country is already one of the most water-stressed in the world, with shrinking aquifers, recurrent drought, and intense competition between agriculture and urban demand. If sunshine hours, a proxy for the solar radiation that supplies the energy for evaporation, are the true dominant lever on ETo, then projections of future water demand need to track cloud cover, aerosol loading, and radiation trends as closely as they track temperature. It also suggests that the well-documented global phenomenon of dimming and brightening, in which the solar radiation reaching the surface has changed over decades due to aerosols and cloudiness, may have played an underappreciated role in Iran’s evolving water balance. Irrigation scheduling tools, drought indices, and crop water models calibrated around temperature sensitivity alone could be systematically misattributing the drivers of demand.

The study also offers a methodological lesson that extends well beyond Iran. Sensitivity analyses of the Penman–Monteith equation have been published for regions around the world, from the Yangtze basin to Australia’s climatic zones to tropical savannas, and their conclusions have often varied. Some of that variation is surely real, reflecting genuinely different climates. But this study demonstrates that part of the disagreement may be an artifact of method: analyses relying on simple correlations can rank drivers differently than analyses that properly account for multicollinearity. As climate change reshapes the individual variables at unequal speeds, knowing which lever actually moves evapotranspiration will determine whether water managers in arid lands anticipate the future accurately or are surprised by it. For Iran, at least, the answer now points upward, toward the sun itself.

Subject of Research: Spatiotemporal drivers of reference evapotranspiration in Iran over 58 years

Article Title: A multi-decadal spatiotemporal evaluation of reference evapotranspiration drivers in Iran: Leveraging an integrated LMG and Corrplot statistical framework

Article References: Zarei, A. R., & Mahmoudi, M. R. (2026). A multi-decadal spatiotemporal evaluation of reference evapotranspiration drivers in Iran: Leveraging an integrated LMG and Corrplot statistical framework. Theoretical and Applied Climatology, 157(11), Article 703. https://doi.org/10.1007/s00704-026-06635-8

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06635-8

Keywords: reference evapotranspiration, Iran, FAO Penman-Monteith, sunshine hours, LMG variance decomposition, multicollinearity, Corrplot, Inverse Distance Weighting, climate variables, water resources, irrigation planning, Theoretical and Applied Climatology

Construction and demolition waste, commonly abbreviated as CDW, is generated in staggering quantities wherever buildings rise and fall. It includes concrete, bricks, mortar, tiles, timber, metals, and a heterogeneous mixture of residual materials that end up stockpiled on vacant lots, dumped along roadsides, or spread informally across land at the edges of expanding neighborhoods. While much of this material appears inert, decades of research have shown that rainwater percolating through rubble can dissolve and carry away a suite of toxic elements. Lead, arsenic, cadmium, copper, zinc, and manganese are among the most frequently detected, and each behaves differently once it enters the soil environment. The new study led by Francis Agyapong of the Department of Materials Engineering at Kwame Nkrumah University of Science and Technology set out to quantify exactly which of these elements are escaping from waste piles in Kumasi, how concentrated they are, and how far they travel through the ground.

The research team collected samples of construction and demolition waste from an active construction site in Kumasi, along with samples of the soil lying directly beneath and around the waste deposits. Back in the laboratory, the samples were digested and analyzed to determine the concentrations of six target heavy metals: lead, arsenic, zinc, copper, manganese, and cadmium. These six were chosen deliberately. Lead and arsenic are notorious for their toxicity even at low exposures, cadmium accumulates in the kidneys and bones of humans and animals, and copper, zinc, and manganese, although essential micronutrients at trace levels, become damaging contaminants when concentrations climb beyond natural background ranges. By comparing metal levels in the waste itself with levels in the underlying soil, the researchers could infer which elements were actively leaching out of the debris and entering the soil column.

One of the study’s most consequential findings concerns the physical character of the waste itself. Analysis showed that the construction and demolition waste in Kumasi is dominated by sand-sized particles. This composition matters enormously for contaminant transport. Coarse, sandy material has large pore spaces between grains, allowing water to infiltrate quickly and percolate downward with minimal filtration. In contrast, fine-grained clays and silts tend to slow water movement and can adsorb metals onto particle surfaces, effectively immobilizing them. The sandy nature of the Kumasi waste therefore acts like a permeable conduit, helping leachate, the contaminated liquid that forms as water passes through waste, move easily through the deposit and into the underlying ground. In practical terms, the rubble piles scattered across the city are not merely unsightly; they function as diffuse, unlined chemical reactors that release their contents whenever it rains.

To move beyond simple concentration measurements, the team employed two complementary modeling approaches to map how contamination spreads through space. The first, inverse distance weighting, or IDW, is a geostatistical interpolation technique that estimates values at unsampled locations based on the measured values at surrounding sampling points, weighting nearby observations more heavily than distant ones. IDW is widely used in environmental monitoring because it is computationally straightforward and produces intuitive contamination surfaces. The second approach was more sophisticated: a radial basis function neural network, or RBF network, a machine learning architecture that uses radial basis functions as activation functions to approximate complex, nonlinear relationships in data. Radial basis function networks have a long history in engineering, having been applied to problems ranging from power system stabilization to function approximation, and they excel at fitting irregular multidimensional surfaces. By training the network on the measured metal concentrations, the researchers could predict both the vertical and horizontal spread of contaminants through the site, generating predictive maps of where pollution is heading rather than merely where it has already been detected.

The modeling results revealed a clear hierarchy of mobility among the six metals. Copper, zinc, and manganese emerged as the most mobile contaminants, spreading appreciably from their source and dispersing through the soil profile. Lead and arsenic, by contrast, tended to remain localized near the waste deposits, showing limited migration. This contrast reflects fundamental differences in geochemical behavior. Zinc and manganese form relatively soluble species under the mildly acidic conditions that often develop in decomposing waste, while copper, though it binds strongly to organic matter, can still be transported when complexing agents are present. Lead and arsenic, meanwhile, tend to adsorb strongly onto iron oxides and clay minerals or precipitate as insoluble compounds, anchoring them in place. The distinction is critical for risk assessment: mobile metals threaten groundwater and wider soil ecosystems, whereas immobile ones create persistent hotspots that can poison anything rooted directly in them and can be remobilized if soil chemistry changes.

The environmental consequences of this mobility are far-reaching. The study highlights soil degradation as a primary concern, since accumulating metals alter soil chemistry, disrupt nutrient cycling, and diminish fertility. Plants growing in contaminated ground take up metals through their roots, introducing toxins into the food chain and suppressing growth. Soil microbial communities, the invisible engines that decompose organic matter and sustain soil health, are also vulnerable, as elevated metal concentrations reduce microbial diversity and activity. In a city like Kumasi, where urban agriculture and informal gardening are common and where many residents rely on shallow wells, the migration of copper, zinc, and manganese through sandy soils raises the specter of contaminated drinking water and contaminated food produced on or near waste dumps. Because these impacts unfold gradually and invisibly, they often escape regulatory attention until damage is well advanced.

What makes the study particularly timely is the trajectory of urbanization in Ghana and across sub-Saharan Africa. Construction booms in rapidly growing cities generate enormous volumes of demolition debris, yet formal recycling infrastructure remains limited, and enforcement of disposal regulations is inconsistent. Much of the region’s CDW is managed informally, with little consideration of where waste is placed or what it contains. Previous international research, including life cycle assessments and leaching studies of recycled concrete aggregates in Europe, has documented similar contamination pathways, but localized data from African cities has been scarce. The Kumasi study helps fill that gap, demonstrating that the tools of modern environmental science, from laboratory digestion analysis to GIS-based interpolation and neural network modeling, can be deployed effectively in data-scarce settings to characterize risks that would otherwise go unmeasured.

The authors do not stop at diagnosis. Their recommendations form a three-pronged response to the problem. First, they call for regulatory reform, arguing that existing frameworks governing construction and demolition waste disposal are inadequate to the scale of the challenge and need stronger standards for where and how debris may be stockpiled. Second, they advocate sustainable construction practices, including better sorting of waste streams, reuse of clean materials, and reduction of hazardous components in building products, so that less toxic material enters the waste cycle in the first place. Third, they recommend chemical stabilization, a remediation technique in which binding agents such as lime, cement, or phosphate amendments are mixed into contaminated soil or waste to immobilize metals chemically, reducing their solubility and preventing further leaching. Applied together, these measures could transform CDW from an unmanaged hazard into a managed resource stream.

The broader lesson of the Kumasi research is that the environmental cost of a building does not end when the building comes down. Rubble piles are active chemical systems whose behavior depends on their composition, their particle size, the local climate, and the chemistry of the ground beneath them. As cities worldwide pursue circular economy goals and seek to recycle more construction waste into new roads, embankments, and aggregates, studies like this one serve as a caution: recycled does not automatically mean safe. Understanding what leaches from waste, how far it travels, and who is exposed downstream is a prerequisite for any genuinely sustainable construction future. For Kumasi, and for the countless cities following the same growth trajectory, the metals hidden in the rubble are now on the scientific record, and the clock on remediation is running.

Subject of Research: Leaching and spatial migration of heavy metals from construction and demolition waste in Kumasi, Ghana

Article Title: Assessment of leached heavy metals from construction and demolition waste

Article References: Agyapong, F., Koomson, B., Kwofie, S., Gidigasu, S. S. R., & Gyamfi, C. (2026). Assessment of leached heavy metals from construction and demolition waste. Environmental Monitoring and Assessment, 198(10), Article 1097. https://doi.org/10.1007/s10661-026-15868-z

Image Credits: AI Generated

DOI: 10.1007/s10661-026-15868-z

Keywords: heavy metals, construction and demolition waste, leaching, soil contamination, Ghana, Kumasi, groundwater, inverse distance weighting, radial basis function neural network, GIS mapping, waste management, chemical stabilization

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