The scale of the decline in some areas is striking. The team divided the district into four watershed groups using a 30-meter digital elevation model and the HydroBASINS Level-12 sub-basin topology, then summarized groundwater observations from wells measured before and after the monsoon in 2005, 2010, 2015, and 2018. In the first watershed group, the mean increase in groundwater depth over the 13-year period was 19.32 meters, and all seven monitored wells in that group showed the water table falling, with individual changes ranging from 5.25 to 31.26 meters. The 95 percent confidence interval around that group mean, from 9.11 to 29.54 meters, was the only one of the four groups that excluded zero, making Group 1 the clearest and most statistically robust hotspot of decline. Group 4 followed with a mean increase of 15.69 meters, but its wells told wildly different stories, ranging from a 29.50-meter shallowing to a 52.42-meter deepening, the widest spread of any group.

What makes these numbers genuinely puzzling is that the most urbanized watershed was not the one sinking fastest. Group 2, which had the highest urban fraction in every mapped year, rising from 20.54 percent of its area in 2010 to 27.32 percent in 2020, recorded a mean groundwater-depth increase of only 7.21 meters. Group 1, by contrast, saw its urban fraction climb from just 3.64 percent to 10.56 percent, yet it suffered the steepest decline. Group 3, with the lowest mean change at 5.76 meters, contained wells that both deepened by more than 40 meters and shallowed by more than 24 meters. The researchers deliberately refrained from running correlations or regressions on these four group-level values, arguing that with so few units, any statistical test would be misleading. Instead, they describe the patterns as non-monotonic: urbanization and groundwater decline co-occur in some places and diverge in others.

To probe whether changing land surfaces might still be shifting the hydrology, the team turned to the Curve Number method, a standard engineering framework that combines land-cover classes with Hydrologic Soil Groups to estimate how much rainfall becomes direct runoff in a storm. Soil data from SoilGrids at 250-meter resolution showed that Hydrologic Soil Group D, the category with the highest runoff potential and lowest infiltration capacity, dominated all four watershed groups. That is why mean Curve Number values were uniformly high, between 89.15 and 89.94 in 2010, creeping up by only 0.24 to 0.64 points by 2020. A scenario analysis using hypothetical rainfall events of 10, 25, 50, and 100 millimeters showed that these small Curve Number shifts would translate into only modest increases in modeled runoff, at most 1.56 millimeters for a 100-millimeter storm. In other words, the land surface became slightly more runoff-prone as cities grew, but not enough to constitute a major hydrological upheaval.

Climate, the other obvious suspect, also failed to provide a clean explanation. Drawing on NASA POWER gridded reanalysis data from 2003 to 2018, quality-controlled down to 204,540 complete daily records across 35 meteorological grid cells, the researchers found that mean annual precipitation varied only modestly among the watershed groups, from 166.4 to 195.2 millimeters, with a spatial coefficient of variation of just 7.7 percent. Group 1 was actually the wettest group and had the lowest estimated reference evapotranspiration, yet it experienced the deepest groundwater decline. Mann-Kendall trend tests found no statistically significant monotonic trends in annual precipitation, monsoon rainfall, mean temperature, or Hargreaves reference evapotranspiration over the 16-year record. The climatic water balance, precipitation minus reference evapotranspiration, was deeply negative in every year and every group, averaging minus 1,744.2 millimeters per year, a stark reminder of how far atmospheric demand exceeds rainfall in this hot desert climate.

The drought record added further nuance. The 12-month Standardized Precipitation Index, computed relative to the study period rather than a long-term baseline, flagged 2003 as moderately dry and 2018, the final monitoring year, as severely dry, with December values of minus 1.46 and minus 1.55 respectively, while 2015 was extremely wet at plus 2.76. Yet groundwater depth generally increased across the monitoring period despite these alternating wet and dry phases. Exploratory correlations between district-mean groundwater depth and antecedent 90-, 180-, and 365-day precipitation, as well as the previous month’s SPI-12, were all negative, hinting that wetter conditions coincided with shallower water tables, but none reached statistical significance. With only eight intermittent groundwater snapshots, the authors caution, such relationships cannot be pushed toward causal interpretation.

This is precisely where the study’s intellectual honesty becomes its most valuable feature. The authors are explicit that their analysis does not simulate groundwater flow, calibrate a water balance, or attribute the observed declines to any single mechanism. Pumping rates, canal seepage, irrigation return flows, aquifer transmissivity, well-screen depths, hydraulic gradients, and river-aquifer exchange were not consistently available across the district, and without them, the spatial associations remain descriptive. The Indus Basin Irrigation System is a labyrinth of conjunctive water use; the larger Lower Bari Doab system alone carries roughly 201 kilometers of main canals, 95 kilometers of link canals, 2,264 kilometers of distributaries, and more than 60,000 tubewells in freshwater areas. Canal leakage and recirculated irrigation water can contribute substantially to aquifer replenishment, meaning that recharge in this landscape is engineered as much as natural, and it can vary dramatically over short distances.

The satellite era has only sharpened the urgency of understanding these local dynamics. Gravity-recovery missions have identified significant terrestrial water-storage losses in the Indus Basin since the mid-2010s, and global assessments have documented rapid groundwater declines, with occasional recovery, in dry regions dominated by irrigated agriculture, from the North China Plain to the US High Plains and the Central Valley. A recent national-scale analysis in Pakistan found regionally heterogeneous groundwater responses to land-use transitions, reinforcing the message that similar surface transformations do not produce uniform subsurface consequences. The Multan study adds a district-level, well-by-well texture to that picture, showing that even within a single administrative area, the mismatch between land-surface indicators and aquifer response can be profound, and that Global Moran’s I tests found no statistically significant district-wide spatial clustering of well-level declines under a four-nearest-neighbor weighting scheme.

The practical payoff is a screening framework rather than a verdict. Group 1 emerges as the unambiguous management priority, with every monitored well declining and a confidence interval that excludes zero, while Group 4 demands attention for its high mean decline and extreme well-to-well variability, and Groups 2 and 3 require finer-scale investigation before uniform measures are imposed. The authors recommend denser monitoring networks with continuous water-level loggers, measurements referenced to a common elevation datum, pumping inventories, canal-delivery and seepage records, aquifer tests, and eventually calibrated groundwater-flow or integrated surface-water-groundwater models. Their broader warning deserves to travel well beyond Punjab: urban fraction and Curve Number alone should never be used to infer groundwater decline. In data-scarce semi-arid plains around the world, organizing whatever observations exist within common hydrological units is a sensible first step, but causal stories about dying aquifers must wait for the pumping, recharge, and subsurface measurements that reveal what is really happening underground.

Subject of Research: Spatial associations between urbanization, runoff potential, and groundwater-depth change in the semi-arid alluvial plain of Multan, Pakistan

Article Title: Spatial associations between urbanization, watershed runoff potential, and groundwater-depth change in a Semi-Arid Alluvial Plain, Multan, Pakistan

Article References: Ullah, A., Wang, Y., Wang, H., liu, J., Sajjad, M. M., Hussain, I., Abbas, H., & Nazim, H. (2026). Spatial associations between urbanization, watershed runoff potential, and groundwater-depth change in a Semi-Arid Alluvial Plain, Multan, Pakistan. Environmental Earth Sciences, 85(16), Article 419. https://doi.org/10.1007/s12665-026-13164-0

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13164-0

Keywords: groundwater, urbanization, Multan, Indus Basin, Curve Number, Hydrologic Soil Groups, NASA POWER, watershed delineation, semi-arid climate, land-use change, water-table decline, hydrogeology

The research, led by Hai Sun of the Ocean University of China together with colleagues at Qingdao University of Technology, Clemson University, and Western Sydney University, tackles one of the most persistent blind spots in flood management: the pathways linking land-use patterns to hydrological responses. Urbanization reshapes basin hydrology by replacing soils and vegetation with impervious surfaces such as roads, rooftops, and parking lots. These surfaces prevent infiltration, accelerate runoff generation, and amplify flood peaks, sending more water into rivers faster and overwhelming channels and drainage systems. Green infrastructure — permeable pavements, rain gardens, wetlands, and vegetated storage areas — counteracts this by improving infiltration, storage, and surface roughness. But until now, planners have lacked a rigorous, basin-scale framework for measuring how these benefits vary across space and under different storm conditions.

To close that gap, the team constructed three urban expansion scenarios for the year 2044 in the Poyang Lake Basin: a no-green-infrastructure baseline, a storage-oriented green infrastructure scenario, and a permeability-oriented scenario. The Poyang Lake Basin, China’s largest freshwater lake system and a critical node in the Yangtze River’s hydrology, is a natural laboratory for this question. Its low-lying floodplains, dense tributary network, and rapidly expanding urban platforms make it acutely sensitive to changes in land cover.

Land-use changes under each scenario were simulated using a cellular automaton–Markov (CA–Markov) model, a technique that combines transition probabilities derived from historical land-change data with spatial neighborhood constraints to project how urban footprints evolve. This allowed the researchers to generate realistic maps of where impervious surfaces would spread by 2044 and where green infrastructure would be deployed under each strategy. The hydrological consequences were then evaluated through a coupled one-dimensional and two-dimensional modeling framework, chaining HEC-HMS, a rainfall-runoff model, to HEC-RAS, a river hydraulics and flood inundation model. The coupling is significant: it enables a continuous simulation chain from land-use evolution to runoff generation to flood dynamics, rather than treating each link in isolation.

The baseline results are sobering. Under unconstrained urban expansion, impervious surface coverage in the basin rises from 3.93 percent to 7.37 percent — nearly a doubling of sealed ground. Correspondingly, the basin-averaged curve number, a standard parameter in the Soil Conservation Service runoff method that encapsulates how readily a landscape converts rainfall into runoff, increases from 68 to 72. That seemingly modest shift carries heavy consequences: the researchers calculate an approximately 18 percent decrease in potential maximum retention, the landscape’s capacity to absorb and store rainfall before it becomes floodwater. In plain terms, by mid-century the basin could surrender nearly a fifth of its natural buffering capacity to concrete and asphalt.

Green infrastructure partially blunts this trajectory, but the two strategies do so very differently. Permeability-oriented green infrastructure — designed to restore infiltration across the urban surface — achieves the strongest reduction in impervious coverage, limiting the rise to 6.68 percent, and effectively reverses the degradation of infiltration and storage capacity captured by the curve number. Storage-oriented green infrastructure, which concentrates water retention in discrete facilities, shows only limited improvement in these landscape-scale parameters. The reason is structural: storage works locally, behind berms and inside basins, while permeability works everywhere, beneath every street and rooftop it touches.

Those parameter-level differences propagate directly into flood behavior. Under small to moderate rainfall events, permeability-oriented green infrastructure reduces peak discharge by 6.2 percent and total runoff volume by 3.54 percent, and — critically — it maintains its effectiveness even under extreme conditions, because infiltration capacity does not fill up the way a storage basin does. Storage-oriented measures, by contrast, remain constrained by finite capacity: once a retention facility fills, additional rainfall passes through unattenuated. Yet the picture is not one-sided. The analysis reveals clear spatial heterogeneity: in areas with sufficient storage capacity, storage-based strategies actually outperform infiltration-based measures during large rainfall events, when rainfall intensity outpaces the soil’s ability to absorb water and detention volume becomes the deciding factor. The catch is that this advantage is limited by spatial and capacity constraints at the basin scale — there is simply not enough suitable land and storage volume to deploy it everywhere.

The two-dimensional hydraulic component of the framework sharpens this spatial story further. Permeability-oriented green infrastructure reduces the extent of high-depth and high-velocity flood hotspots, with the strongest benefits concentrated along river corridors and urban platforms — precisely the locations where people and assets cluster. This leads the authors to a practical siting logic built around the curve number itself: boost permeability and roughness on the slopes, maintain them through the channels, and add storage capacity at confluences where flows converge and backwater effects amplify. In other words, read the landscape’s hydrological fingerprints and match the intervention to the terrain.

The study’s most consequential recommendation is that no single strategy suffices. Because green infrastructure performance varies with location, terrain, and storm magnitude, the authors argue for differentiated strategies aimed at residual high-risk areas: infiltration-dominated measures for moderate rainfall, enhanced storage and conveyance capacity for extreme events, and spatially targeted deployment that combines infiltration and storage synergies in low-lying zones, storage-control combinations along flow corridors, and strict land-use regulation in high-risk areas. This is a departure from the one-size-fits-all deployments that have characterized much green infrastructure planning, including China’s high-profile sponge city program, and it provides theoretical support and decision-making guidance for coordinated planning and refined flood risk management at the basin scale.

The timing could hardly be more urgent. Recent global research has documented rapid urban growth inside flood zones since 1985 and projected substantial increases in future fluvial flood risk across China’s major urban agglomerations, with the Global South bearing disproportionately higher exposure. As climate change intensifies extreme rainfall and cities continue to seal their surfaces, the Poyang Lake findings offer a template that travels: simulate the land, couple it to the water, and let the spatial heterogeneity of performance — not generic best practice — dictate where every permeable meter and every storage basin goes. The difference, this research suggests, may be measured not just in percentage points of peak discharge, but in neighborhoods that stay dry.

Subject of Research: Spatial variability in the flood mitigation performance of green infrastructure under future urban expansion in the Poyang Lake Basin

Article Title: Spatial heterogeneity of green infrastructure performance in flood mitigation under urban expansion

Article References: Sun, H., Wang, H., Chu, Y., Yao, W., Fan, C., Chu, Z., & Liang, B. (2026). Spatial heterogeneity of green infrastructure performance in flood mitigation under urban expansion. Natural Hazards, 122(19), Article 637. https://doi.org/10.1007/s11069-026-08312-5

Image Credits: AI Generated

DOI: 10.1007/s11069-026-08312-5

Keywords: green infrastructure, flood mitigation, urban expansion, Poyang Lake Basin, CA-Markov model, HEC-HMS, HEC-RAS, curve number, impervious surfaces, runoff reduction, spatial heterogeneity, flood risk management