The model at the heart of the study, called GLOBGM v1.1, is a two-layer, transient groundwater flow model built on a parallelized prototype of the widely used MODFLOW6 code. To make a global simulation at 30-arcsecond resolution computationally feasible, the researchers split the planet into four independent but coupled continental-scale models covering Afro-Eurasia, the Americas, Australia, and the remaining islands. Each of these is further partitioned into non-overlapping sub-models that are solved together within the MODFLOW linear solver. The two-layer design allows the model to represent unconfined, confined, and semi-confined aquifer systems: where a confining layer exists, the upper layer represents that barrier and the lower layer the pressurized aquifer beneath it. Forcing data, including runoff, recharge, and groundwater abstraction, come from the global hydrological model PCR-GLOBWB2, resampled from its native 5-arcminute grid down to the finer GLOBGM grid.

Two technical innovations distinguish this version from its predecessor. The first is a dynamic drainage elevation routine, in which the elevation at which water drains from the groundwater system adjusts according to soil saturation. As soils become wetter, the drainage level rises toward the surface; as they dry, it drops. This coupling, based on the improved ARNO scheme for identifying saturated areas, keeps simulated water levels closer to the surface in wetland regions and suppresses unrealistic head fluctuations and intermittent discharge events. The second is a statistical downscaling of groundwater recharge. Because PCR-GLOBWB2 is known to underestimate recharge, particularly in arid and semi-arid regions where recharge is low and episodic, the team used a Generalized Additive Model trained on observed recharge data to produce a kilometer-scale recharge field, then applied a multiplicative correction factor to the model’s monthly recharge inputs, capped at local precipitation to preserve numerical stability.

Calibration was a major undertaking. The team tested 162 permutations of prefactors applied to hydraulic conductivity, anisotropy, and entrance resistance, varying the conductivity adjustments by lithological class drawn from the global GLiM dataset. Candidate settings were evaluated against long-term average water table depths from 34,800 observation wells in IGRAC’s Global Groundwater Monitoring Network, deliberately excluding known groundwater abstraction hotspots so that depletion would not be confused with model error. The best-performing configuration reduced the mean bias in water table depth from minus 4.8 meters in the uncalibrated predecessor to 3.6 meters, with the depth-weighted bias falling from 34.2 to 32.5 meters. The improvements were most pronounced for shallow water tables, which matter most for human use and for sustaining ecosystems through baseflow.

Validation against independent observations showed that the model produces skillful predictions, measured by a non-parametric Kling-Gupta Efficiency skill score, in approximately 83 percent of shallow to intermediate wells with water tables between zero and 60 meters deep when compared at the monthly scale, and in about 60 percent of deep wells exceeding 60 meters. At the annual scale, 71 percent of shallow wells and 66 percent of deep wells were skillful. The depth dependence is physically sensible: deeper aquifers respond slowly to recharge through thick unsaturated zones that the model does not explicitly resolve. To check performance at scales where point observations are absent, the team compared simulated total water storage anomalies with satellite gravimetry from the GRACE and GRACE-FO missions, finding agreement in trend direction over 62.3 percent of the global land surface and agreement in both direction and magnitude over 45.7 percent, consistent with other state-of-the-art global hydrological models.

The researchers also added a machine learning bias correction. A two-stage LightGBM ensemble, one classifier predicting whether the model overestimates or underestimates groundwater heads and one regressor predicting the magnitude of the error, was trained on nearly 97,000 observation wells using model parameters, Google’s AlphaEarth geospatial embeddings, and absolute latitude as a proxy for aridity. Training and testing were separated spatially using K-means clustering to prevent autocorrelation from inflating performance. The correction correctly predicted the direction of bias in 80 percent of holdout cases and shifted median bias toward zero across all depth classes, with the largest gains in deep wells. An asymmetric ramp-factor scheme prevents the correction from pushing simulated water tables unrealistically close to the land surface, preserving the model’s physics where it matters most.

The historical simulation from 1960 to 2019, driven by the ISIMIP3a observational forcing, reproduced the world’s best-documented groundwater depletion zones with striking fidelity: the southern and central High Plains of the United States, Iran, the Arabian Peninsula, and the Indo-Gangetic Plain all show the deepening water tables that independent observations have recorded. The model also identified rising water tables in regions known to be gaining storage, including the Guarani aquifer in South America, the northern Great Plains of North America, and the Karoo Basin in South Africa. Intriguingly, it flagged rising water tables across northern latitudes and the Arctic, potentially linked to climate-driven increases in recharge, though the authors caution that these areas overlap with permafrost zones where the model is known to be less reliable.

The future projections, run with five global climate models under three combined socioeconomic and climate scenarios from ISIMIP3b, paint a nuanced picture. Globally, more regions are projected to see rising than falling water tables through 2100, consistent with an intensifying hydrological cycle, and the pace of change increases from the low-emission SSP1-RCP2.6 scenario to the high-emission SSP5-RCP8.5. Water tables are projected to rise across much of the Sahel, the Horn of Africa, North America’s Pacific Coast, the Tibetan Plateau, and the Arabian Peninsula, while declines are projected for the Mediterranean, southeastern Australia, and the African and South American tropics. Europe stands out as the sole continent where water tables are expected to deepen markedly, particularly under high emissions. Known depletion hotspots are expected to persist, and the model projects continued declines over southern Europe, the Amazon, California’s Central Valley, and the High Plains, diverging from earlier studies that omitted groundwater abstraction, a difference the authors attribute to the inclusion of human water demand.

The team is transparent about the model’s limits. Simulations are flagged as less reliable in karst aquifers, where Darcian flow assumptions break down; in permafrost regions, where frozen ground complicates recharge; and in steep mountainous terrain, where the schematization is too coarse to resolve hillslope hydrology. Quality assurance maps accompany the dataset so users can interpret results appropriately. All simulations, along with the model code, are openly accessible, and because the outputs follow the ISIMIP framework, they can be combined with other Earth system datasets to study how changing groundwater affects ecosystems, food production, and society. In a world where the United Nations has warned of an era of global water bankruptcy, a kilometer-scale, globally consistent view of where water tables are falling and rising may prove one of the most consequential datasets yet produced for planning a sustainable water future.

Subject of Research: Global hyper-resolution modeling of historical and future groundwater dynamics at approximately 1 km resolution

Article Title: Global hyper-resolution modeling of historical and future groundwater dynamics

Article References: van Jaarsveld, B., Wanders, N., Otoo, N. G., Sutanudjaja, E. H., Verkaik, J., Zamrsky, D., & Bierkens, M. F. P. (2026). Global hyper-resolution modeling of historical and future groundwater dynamics. Earth System Dynamics, 17(5), 1201-1236. https://doi.org/10.5194/esd-17-1201-2026

Image Credits: AI Generated

DOI: 10.5194/esd-17-1201-2026

Keywords: groundwater, hyper-resolution modeling, GLOBGM, water table depth, climate change, groundwater depletion, MODFLOW, ISIMIP, machine learning bias correction, GRACE satellite, SSP-RCP scenarios, water resources management

The study, published in Environmental and Sustainability Indicators by Qing Lu, Baifei Ren, Liyan Xu, and Hang Yin, introduces what the authors call the National Ecological Security Pattern for Groundwater Preservation, or NESP-GP. The framework builds on the ecological security pattern theory first articulated in 1996, which has long guided landscape planning in China by identifying the spatially critical areas that must be protected for ecosystems to survive. But while that approach has been applied to lakes, typhoons, and broad ecological infrastructure, groundwater—the invisible backbone of water security—has largely been left out. The new work aims to close that gap by treating groundwater not as a marginal afterthought in regional planning, but as a central organizing principle.

Constructing the map was a formidable data challenge. Groundwater information in China is famously fragmented, scattered across ministries, provinces, and agencies, with open-source datasets extremely scarce. The researchers pulled together an extraordinary array of sources: groundwater level yearbooks from monitoring wells concentrated in key regions such as the North China Plain, the Sanjiang Plain, the Weihe Plain, and the Hexi Corridor; the Atlas of Groundwater Resources and Environment of China; provincial water resources bulletins; and layers from the Resource and Environment Data Cloud Platform. Everything was harmonized to a uniform one-kilometer grid, a resolution chosen as a pragmatic compromise between national coverage and computational feasibility. The authors are careful to stress that the result is a strategic spatial assessment, not a high-resolution portrait of local aquifer conditions.

The heart of the framework lies in its three-tiered indicator system. Nineteen indicators were organized into three categories: recharge, utilization, and socio-environmental impacts. The recharge tier captures the natural capacity of the land to replenish aquifers, incorporating average annual precipitation, topographic slope, lineament density, soil texture, depth to groundwater, aridity index, drainage density, geomorphology, and total groundwater resources. The utilization tier measures human demand, from population and GDP density to the proportion of groundwater in the total water supply and per capita water resources. The final tier captures environmental and policy pressures, including land use change, groundwater salinity, the rate of groundwater table decline, and the cone of groundwater depression—the telltale funnel formed beneath areas where pumping has dramatically outpaced replenishment.

To weight these indicators, the team assembled a panel of thirty experts from Peking University and the Chinese Research Academy of Environmental Sciences, spanning hydrogeology, ecology, landscape architecture, urban planning, geography, and environmental management. Each expert independently scored pairwise comparisons using the Analytic Hierarchy Process, a structured decision-making method developed by Thomas Saaty, with consistency ratios below 0.10 deemed acceptable. The resulting weights assigned just under half of the total influence to recharge conditions, about 31 percent to utilization pressures, and roughly 20 percent to socio-environmental impacts. A sensitivity analysis that perturbed the three top-level weights by plus or minus 20 percent showed the national pattern held firm, with Spearman correlations against the baseline ranging from 0.956 to 0.992 and priority-area overlap reaching as high as 96 percent. The map, in other words, is not an artifact of expert idiosyncrasy.

When the index was classified into five security levels using the Natural Breaks method, the geography of China’s groundwater vulnerability emerged with striking clarity. Low security areas covered the largest share of the country, 31.2 percent, followed by very low security areas at 25.3 percent and medium security areas at 21.7 percent. High security zones accounted for 15 percent, while the most alarming category—extremely low security—encompassed 6.8 percent of the study area. Taken together, nearly a third of the nation falls into the two worst tiers, where groundwater protection should command the highest management priority. The spatial logic differs sharply by region. In Northwest China, the problem is a brutal combination of meager natural recharge, driven by low precipitation and high aridity, and heavy dependence on groundwater for livelihoods. In the North China Plain, by contrast, recharge conditions are comparatively favorable, but intensive agriculture, dense populations, and a legacy of relentless extraction have imposed crushing utilization pressures.

The researchers validated their map in two independent ways. First, they compared GRACE satellite observations—which track changes in Earth’s gravity field to estimate shifts in terrestrial water storage—with groundwater monitoring well records over the period 2006 to 2016. The two datasets showed statistically significant but strongly scale-dependent agreement: Spearman’s correlation rose from a modest 0.127 at the native GRACE grid scale to 0.396 when observations were aggregated across 16-by-16 blocks, with directional agreement climbing to 68.1 percent. The lesson is that agreement between satellite and well data improves as their spatial footprints align, offering evidence of broad regional consistency rather than aquifer-level validation. Second, the team compared their security index against China’s official 2024 groundwater overexploitation delineation across eighteen provinces. Provinces with larger officially designated overexploitation areas showed significantly lower NESP-GP values, with a Spearman correlation of minus 0.505, lending independent external credibility to the national pattern.

The most provocative application of the new map concerns China’s sponge cities. The sponge city program, launched through a 2014 technical guide from the Ministry of Housing and Urban-Rural Development, aims to make urban landscapes absorb, store, purify, and reuse rainwater through low-impact development. Its core planning indicator, the volume capture ratio of annual rainfall, is currently zoned almost entirely by regional precipitation patterns—subsurface conditions, groundwater dynamics, and soil characteristics receive comparatively little consideration. By overlaying the NESP-GP onto the existing zoning, the researchers identified striking mismatches in parts of southern Hubei, northern Shandong, and the Jianghan Plain, where groundwater security conditions diverge from the rain-centric planning framework. In water-scarce northern China, they argue, sponge city construction should incorporate reasonable groundwater recharge indicators informed by local precipitation and groundwater level data, much as a case study in Jinan showed that a 75 percent volume capture ratio could help preserve the city’s famous springs by boosting recharge.

The authors are careful about the limits of their claims. The overlay analysis is explicitly a screening-level planning reference, not a recalibration of engineering targets, and any adjustments to capture ratios should be grounded in site-specific hydrological and hydrogeological assessment. The datasets underlying the map largely reflect conditions up to 2016, which means the notably low scores for the North China Plain do not capture more recent groundwater recovery attributed to the South-to-North Water Diversion Project and tightened extraction controls—developments that recent research suggests have significantly reversed the long-term decline in that region’s water table. The framework is meant to complement, not replace, process-based hydrogeological modeling, and the researchers recommend that future versions incorporate water quality data, pollutant monitoring, updated observations, and comparisons between expert-derived weights and data-driven alternatives.

Even with those caveats, the study represents a conceptual shift with implications far beyond China’s borders. Groundwater supplies a substantial share of the world’s irrigation and drinking water, and aquifer depletion from the North China Plain to California’s Central Valley has repeatedly shown that recovery, when possible at all, takes decades. By embedding groundwater into the ecological security pattern framework, the researchers offer governments a template for shifting water management away from narrow, single-source supply thinking and toward integrated spatial planning that treats aquifers as living ecological infrastructure. The invisible water beneath our feet, the study suggests, can no longer afford to remain invisible in the places where cities are planned—and as climate change tightens its grip on global water cycles, the maps that reveal it may become some of the most important documents of the century.

Subject of Research: A national ecological security pattern framework for assessing groundwater sustainability in China and integrating groundwater considerations into sponge city planning

Article Title: A National Ecological Security Pattern for Groundwater sustainability in China: Implications for sponge city planning

Article References: Lu, Q., Ren, B., Xu, L., & Yin, H. (2026). A National Ecological Security Pattern for Groundwater sustainability in China: Implications for sponge city planning. Environmental and Sustainability Indicators, 32, Article 101536. https://doi.org/10.1016/j.indic.2026.101536

Image Credits: AI Generated

DOI: 10.1016/j.indic.2026.101536

Keywords: groundwater, ecological security pattern, China, sponge city, GRACE satellite, aquifer depletion, North China Plain, analytic hierarchy process, water resources management, spatial planning, groundwater recharge, low impact development