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 Wei River Basin is the largest tributary of the Yellow River and a lifeline for tens of millions of people. It irrigates extensive cropland, sustains major urban centers, and has historically mediated the delicate balance between agricultural output and ecological health on the semiarid Loess Plateau. Decades of intensified human activity have already reshaped the basin’s hydrological processes, and questions about how much water will be available in the coming decades have become a central concern for planners and policymakers. Previous research has often examined climate change or land-use change in isolation, which leaves a critical gap: the two drivers interact, and their combined effects can differ substantially from what either would produce alone.

To close that gap, the research team—led by Yating Gao, Ganggang Zuo, Jiancang Xie, Ni Wang, Zheng Liu, and Tianfan Wang—built a framework that chains together three complementary modeling tools. The first is the Taylor diagram, a widely used statistical visualization developed by climate scientist Karl Taylor that summarizes how well a model reproduces observed patterns by comparing correlation, variance, and root-mean-square error in a single plot. In this study, the Taylor diagram served as a rigorous screening device for general circulation models, allowing the team to identify which global climate models best captured the basin’s historical climate behavior before trusting their future projections. This step addresses one of the persistent weaknesses in scenario studies: model uncertainty, which can propagate from coarse global simulations all the way into local water-resource estimates.

The second component is the Markov-PLUS model, a land-use simulation approach that merges a Markov chain’s ability to quantify transition probabilities between land categories with the PLUS model’s strength in generating spatially realistic land-change patterns. PLUS, short for patch-generating land use simulation, uses machine learning to understand the drivers behind historical land conversions and then produces future landscapes patch by patch, respecting both neighborhood effects and the underlying suitability of terrain. By coupling Markov-chain projections of how much land will change with PLUS’s determination of where that change will occur, the team generated land-use maps for the future under multiple development trajectories aligned with the shared socioeconomic pathways.

The third and final component is the Soil and Water Assessment Tool, or SWAT, a physically based, semi-distributed hydrological model that has become a global standard for watershed analysis. SWAT divides a basin into sub-basins and further into hydrological response units defined by soil type, land cover, and slope, then simulates the full water balance—including precipitation inputs, evapotranspiration, infiltration, surface runoff, and lateral and groundwater flows. Running SWAT with downscaled climate projections and the simulated future land-use maps allowed the researchers to quantify how water yield evolves across space and time under each scenario combination.

The scenarios examined follow the coupled SSP-RCP framework, which links socioeconomic storylines with representative concentration pathways describing different levels of future radiative forcing. The results on the climate side are unambiguous. Across all scenarios, the study finds increasing trends in precipitation, maximum temperature, and minimum temperature within the basin, with the largest temperature increases occurring under the high-emission SSP585 scenario. While rising precipitation might seem like good news for a water-stressed region, warmer temperatures drive up evapotranspiration—the return of water from soil and vegetation to the atmosphere—so more rainfall does not automatically translate into more available water. The interplay between these competing effects lies at the heart of the water-yield question.

On the land side, the Markov-PLUS simulations captured a consistent structural transformation across all development trajectories: continuous expansion of built-up land at the expense of cropland, with the most pronounced land-use changes again appearing under SSP585. Urbanization seals surfaces, alters infiltration, and changes the routing of water through the landscape, which is precisely why including realistic land dynamics matters for hydrological forecasting. The model’s ability to reproduce the basin’s historical land-use patterns gave the researchers confidence that its future simulations were grounded in credible transition dynamics rather than arbitrary assumptions.

Perhaps the most consequential finding comes from the attribution analysis. When the team separated the effects of climate change from those of land-use change, they found that variations in future water yield are primarily dominated by climatic effects, while land-use effects remain relatively limited in comparison. However, the interaction between the two drivers becomes increasingly significant under the SSP585 scenario, suggesting that in a high-emission world, the way land is managed will matter more as a modulator of hydrological outcomes than it does under milder pathways. This asymmetry carries a practical message: mitigation of greenhouse gas emissions remains the dominant lever for protecting the basin’s water resources, but land-use planning retains a meaningful, and growing, secondary role.

The spatial anatomy of the projections is equally revealing. Water yield in the Wei River Basin follows a clear decreasing gradient from south to north, reflecting the basin’s climatic transition from wetter mountainous headwaters in the south to the drier Loess Plateau in the north. Sub-basins in the central and lower reaches exhibit relatively higher water yield, whereas tributary and upstream areas show lower values. This geographic heterogeneity means that the impacts of declining yield will not be felt uniformly: communities and ecosystems in the northern and upstream portions of the basin, already operating closer to their hydrological margins, face the greatest relative vulnerability.

The temporal projections add urgency to the diagnosis. Annual hydrological water yield is projected to decline under all scenarios over the coming decades, with the greatest reduction occurring under SSP585 and the most pronounced monthly decreases concentrated between February and July. That seasonal window is far from arbitrary—it spans the late winter recession and the critical early growing season, when crops depend on soil moisture and streamflow and when reservoir operations must balance storage against downstream demands. A shrinking yield precisely when agricultural and ecological water needs ramp up compounds the challenge of adapting to climate change in one of China’s most historically water-constrained regions.

The authors frame their work as a contribution to climate-adaptation planning and watershed-scale water-resource assessment, and the integrated design of the framework is its central innovation. By screening climate models with Taylor diagrams, simulating land futures with Markov-PLUS, and translating both into hydrological outcomes with SWAT, the approach systematically captures coupled dynamics that single-driver studies miss. The findings offer scientific grounding for decisions about where to prioritize water conservation, how to schedule reservoir releases, and which sub-basins deserve the most attention in adaptation strategies. They also underscore a sobering reality for the Yellow River system and semiarid basins worldwide: even with somewhat increased precipitation, warming may overwhelm gains, leaving less water flowing through the landscape than the region has come to rely on. For the millions who depend on the Wei River, the study’s message is that the coming decades demand not just awareness of change, but deliberate, spatially informed preparation for it.

Subject of Research: Coupled effects of future climate and land-use change on hydrological water yield in the Wei River Basin, China, assessed under SSP-RCP scenarios

Article Title: Coupled Effects of Climate and Land-Use Changes on Hydrological Water Yield in the Wei River Basin of China under SSP-RCP Scenarios

Article References: Gao, Y., Zuo, G., Xie, J., Wang, N., Liu, Z., & Wang, T. (2026). Coupled Effects of Climate and Land-Use Changes on Hydrological Water Yield in the Wei River Basin of China under SSP-RCP Scenarios. Natural Resources Research. https://doi.org/10.1007/s11053-026-10775-z

Image Credits: AI Generated

DOI: 10.1007/s11053-026-10775-z

Keywords: Wei River Basin, water yield, climate change, land-use change, SSP-RCP scenarios, SWAT model, Markov-PLUS, Taylor diagram, hydrological modeling, Yellow River, climate adaptation, CMIP6