The study was conducted during a three-day dry-season campaign from 10 to 12 December 2025, a period chosen deliberately because dry road surfaces maximize the availability of loose dust that vehicles and wind can throw into the air. The team, led by Yemisi Aramide Tijani of Redeemer’s University with colleagues from Olabisi Onabanjo University, Redeemer’s University, and the Federal University of Technology, Ikot Abasi, selected seven paved and seven unpaved roads spanning major traffic arteries and quieter residential streets. At each site, measurements were taken at roughly two meters above the ground—the height of the human breathing zone—during both rush-hour and off-rush-hour periods, with seven replicate readings per period. In total, the campaign produced 296 valid observations, recorded using a calibrated Aeroqual Series 300 monitor operating on optical light-scattering principles, alongside a handheld WindMate weather station tracking temperature, relative humidity, and wind speed.

The numbers that emerged are striking. On paved roads, mean PM2.5 concentrations ranged from about 26 to 215 micrograms per cubic meter, while on unpaved roads they climbed from roughly 66 to nearly 248 micrograms per cubic meter. Coarse particles told an even more dramatic story: PM10 on paved roads spanned 79 to about 1,405 micrograms per cubic meter, but on unpaved roads reached as high as 1,871 micrograms per cubic meter, with a single extreme short-term reading of 2,742 micrograms per cubic meter recorded at one site. For context, the World Health Organization’s 24-hour guideline values sit at 15 micrograms per cubic meter for PM2.5 and 45 for PM10—orders of magnitude below what was measured. The highest concentrations clustered at busy junctions such as Adatan/Carwash and Elega/Bode-Olude, while the lowest values appeared on Okeero, a minor residential street with little commercial activity.

Statistical analysis confirmed that this spatial variation was not random. One-way analysis of variance revealed highly significant differences in both PM2.5 (F = 13.248, p < 0.001) and PM10 (F = 75.473, p < 0.001) across the fourteen sites, with the coarse fraction showing the strongest heterogeneity. Geospatial interpolation of the measurements produced concentration maps revealing distinct high-pollution zones in the central and southern portions of the metropolis, aligned with high-traffic corridors and mechanically disturbed surfaces, while the northeastern sector registered comparatively lower burdens. The authors caution that these interpolated hotspots are visual estimates based on only fourteen monitoring points rather than statistically confirmed pollution zones, but the broad pattern is clear: particulate pollution in Abeokuta is spatially clustered around key activity nodes.

One of the study’s most intriguing findings concerns the relationship between traffic timing and particle behavior. On paved roads, mean PM2.5 fell from 97.73 micrograms per cubic meter during rush hours to 61.04 during off-rush periods—the intuitive expectation that fewer vehicles mean cleaner air. But on unpaved roads the pattern reversed, with PM2.5 rising from 97.14 to 132.90 micrograms per cubic meter between the same periods. This opposing behavior produced a traffic-period by road-surface interaction that approached statistical significance with a substantial effect size (p = 0.069, partial eta squared = 0.249), meaning roughly a quarter of the explainable variance was tied to this interaction. The result suggests that on dusty roads, vehicle movement itself—through wheel-induced turbulence and mechanical resuspension of loose surface material—can dominate the particulate signal in ways that do not simply track the number of vehicles present.

To probe the origin of the particles, the researchers turned to the PM2.5/PM10 ratio, a diagnostic tool that distinguishes coarse-mode dominance from fine-particle contributions. Ratios below 0.5 indicate that coarse particles—typically generated by mechanical processes such as road-dust resuspension—predominate, while higher ratios point toward combustion-related fine particles from exhaust, biomass burning, or secondary aerosol formation. Several of the highest-concentration zones in Abeokuta corresponded to low-to-moderate ratios, supporting the inference that mechanical dust generation along unpaved and heavily trafficked corridors is a major driver of the city’s particulate burden. The authors are careful to note, however, that the ratio is not source-specific: without chemical tracers such as black carbon or elemental composition data, specific emission sources cannot be conclusively identified.

Meteorology emerged as a secondary but meaningful influence on the city’s aerosol field. Wind speed varied significantly across sites (F = 4.517, p < 0.001), ranging from 0.27 to 2.10 meters per second, and showed positive—though not statistically significant—associations with both particle fractions on unpaved roads, consistent with wind-assisted resuspension of loose material. Relative humidity also differed significantly between locations, with values at many unpaved sites exceeding 70 percent, a level at which hygroscopic particle growth and enhanced deposition can alter atmospheric residence times. Temperature, by contrast, remained relatively uniform across the study area and showed no significant spatial variation, ruling it out as a primary driver of the observed particulate differences. On unpaved roads, PM10 was significantly and negatively correlated with temperature (r = −0.788, p = 0.035), while fine and coarse particles were strongly co-varying on both road types, with correlation coefficients of 0.967 on paved and 0.927 on unpaved roads.

When the measured concentrations were converted into Air Quality Index values using the United States Environmental Protection Agency’s standardized method, the results painted a sobering picture of public health risk. PM2.5 AQI values ranged from 81 to 265 on paved roads and 156 to 298 on unpaved roads, spanning categories from moderate to very unhealthy for all population groups. PM10 fared worse still, with AQI values reaching beyond 500—the top of the scale—at multiple locations, corresponding to hazardous conditions. Most of the sampled roads, both paved and unpaved, fell within the very unhealthy to hazardous range, implying significant risks for children, the elderly, people with respiratory or cardiovascular conditions, and roadside vendors who spend extended hours exposed to traffic corridors. The authors emphasize that these are short-term dry-season roadside measurements rather than formal 24-hour guideline exceedances, but the magnitude of the values leaves little room for comfort.

The findings carry direct implications for how rapidly urbanizing African cities manage their air. The study’s practical recommendations center on progressively paving or chemically stabilizing heavily trafficked unpaved roads, deploying dust suppression during prolonged dry periods, enforcing lower vehicle speeds on dust-prone routes, and improving roadside cleaning practices in ways that avoid dry sweeping, which itself resuspends dust. Strengthened vehicle-emission inspection along major corridors is also urged. More broadly, the work demonstrates that road-surface condition, traffic activity, and meteorology must be considered together rather than in isolation when assessing urban particulate pollution—a lesson with resonance well beyond Nigeria, given that unpaved-road dust emissions are a recognized problem across the developing world. The authors acknowledge the limitations of their short campaign: no collocation with reference-grade instruments, no direct measurement of road-surface silt loading or moisture, and no continuous multi-day monitoring. They call for future work incorporating wet- and dry-season sampling, 24-hour measurements, and vehicle classification to build a fuller picture. For now, the message from Abeokuta is unambiguous: what a road is made of may matter as much as what drives on it.

Subject of Research: Spatial variability of PM2.5 and PM10 particulate pollution across paved and unpaved urban roads and its meteorological controls in Abeokuta, Nigeria

Article Title: Spatial variability and meteorological associations of pm₂.₅ and pm₁₀ across paved and unpaved roads in Abeokuta, Nigeria

Article References: Tijani, Y. A., Olukayode, O. O., Ojo, O. T., & Agbasi, O. E. (2026). Spatial variability and meteorological associations of pm₂.₅ and pm₁₀ across paved and unpaved roads in Abeokuta, Nigeria. Discover Cities, 3(1), Article 191. https://doi.org/10.1007/s44327-026-00371-4

Image Credits: AI Generated

DOI: 10.1007/s44327-026-00371-4

Keywords: particulate matter, PM2.5, PM10, unpaved roads, road dust resuspension, air quality index, Abeokuta, Nigeria, urban air pollution, meteorology, spatial interpolation, dry season

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

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

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

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

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

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

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

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

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

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

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

Image Credits: AI Generated

DOI: 10.1007/s12145-026-02243-2

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

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