The challenge that motivated the study is deceptively simple to state but remarkably hard to solve. Groundwater-induced land subsidence is one of the most widespread geohazards associated with excessive groundwater exploitation, and it produces broad, uneven patterns of sinking that carry serious economic and social consequences. Traditional monitoring tools such as leveling benchmarks and GPS receivers provide accurate measurements, but only at discrete points scattered across the landscape. That sparse coverage makes it nearly impossible to characterize the full regional footprint of subsidence, let alone to work out which buried layers are compressing and driving the deformation. Remote sensing techniques such as Interferometric Synthetic Aperture Radar, or InSAR, can map surface deformation over huge areas in all weather, but by themselves they reveal little about what is happening hundreds of meters underground, where the real mechanical action unfolds within layered aquifer and aquitard systems.

The new framework, developed by Hongwei Sang, Ke Fang, Qimeng Liu, Liang Yuan, and Bin Shi, closes that gap in three sequential steps. First, the researchers applied Persistent Scatterer InSAR, known as PS-InSAR, to derive detailed regional maps of surface deformation. Second, they quantified the temporal relationship between the satellite-derived deformation and regional terrestrial water storage, or TWS, using Pearson correlation analysis, allowing them to flag areas where sinking was strongly tied to hydrological variation. Third, they turned to distributed fiber optic sensing, DFOS, installed in a borehole, to identify the specific compressible strata responsible for the deformation at depth. Each step covers a different scale, and together they build a continuous chain of evidence from regional surface patterns all the way down to individual geological layers.

For the satellite component, the team processed 171 images acquired by the Sentinel-1A satellite between January 2017 and December 2022. The C-band radar instrument revisits the area every twelve days, and the researchers used the SNAP, ISCE, and StaMPS processing chain to extract line-of-sight deformation velocities at roughly half a million monitoring points across the study area. Careful corrections were applied throughout, including precise orbit correction, removal of the topographic phase using the 30-meter Shuttle Radar Topography Mission digital elevation model, three-dimensional phase unwrapping, and atmospheric correction based on the Generic Atmospheric Correction Online Service, GACOS, which uses high-resolution weather-model products to strip out delays introduced by the atmosphere. Because the study area sits on an exceptionally flat alluvial plain at an average elevation of less than six meters, the line-of-sight measurements serve as a reasonable first-order proxy for vertical settlement.

The satellite data painted a telling picture. Line-of-sight velocities ranged from minus 22 millimeters per year, indicating movement away from the sensor, to plus 10 millimeters per year. The most pronounced deformation anomalies appeared in the central sector of the study area, coinciding with zones of land subsidence documented in earlier investigations. But surface maps alone cannot answer the crucial question of cause. To address that, the researchers turned to terrestrial water storage data from the GLDAS V2.2 CLSM product, obtained through the Google Earth Engine platform at a spatial resolution of about 25 kilometers and a daily temporal resolution. TWS aggregates groundwater, soil moisture, surface water, snow and ice, canopy interception, and wet biomass, and thus serves as a regional-scale indicator of overall water storage conditions, including groundwater variability.

The correlation analysis revealed a striking story of changing hydro-mechanical behavior. During 2017, areas with Pearson correlation coefficients between 0.4 and 1.0, the threshold the team adopted to flag moderate to strong correspondence, clustered in the central part of the study area, precisely where deformation rates were highest. Two representative InSAR points, P1 and P2, showed correlation coefficients of 0.958 and 0.811 respectively in 2017, indicating an almost synchronized dance between falling water storage and accumulating subsidence during a period when groundwater exploitation still dominated the system. Similar positive correlations persisted in 2018 and 2019, though generally weaker.

Then, after 2020, something changed. The spatial distribution of correlation coefficients shifted toward lower and even negative values, and the percentage of InSAR pixels exceeding the 0.4 threshold fluctuated downward, from 16.9 percent in 2017 to 4.3, 9.1, and 7.1 percent in 2020 through 2022. Regional terrestrial water storage began to rise, reflecting the gradual recovery of water storage following the enforcement of groundwater extraction restrictions, yet surface deformation continued to accumulate regardless. The explanation lies in the physics of consolidation. According to Terzaghi’s consolidation theory, when groundwater recovers, pore-water pressure increases within the permeable confined aquifer, but the surrounding low-permeability aquitards dissipate their excess pore pressure far more slowly. Residual consolidation of these compressible strata therefore continues long after the hydrological turnaround, producing a time lag between hydrological recovery and surface deformation. The subsurface, in effect, keeps replaying the memory of past over-extraction.

To pin down exactly where that lingering deformation resides, the team exploited the borehole fiber optic monitoring system at borehole SZ1, established in 2013 in Shengze Town, a silk-textile hub in southeastern Suzhou whose township industries once drove severe aquifer over-extraction. The system uses a Brillouin Optical Time Domain Reflectometer, BOTDR, which measures the Brillouin frequency shift generated by the interaction between incident light and acoustic phonons within the optical fiber, converting that shift into axial strain. Configured with a 20-kilometer measurement range, one-meter spatial resolution, and a strain accuracy of plus or minus 50 microstrain, the fiber laid continuously along the borehole records strain profiles across the full vertical extent of the layered aquifer system, something no array of point sensors could achieve.

The fiber optic observations delivered a clear verdict on depth. Negative strain, representing compressive deformation, concentrated within the Af2 confined aquifer and the adjacent Ad2 and Ad3 aquitards. Intriguingly, although groundwater extraction occurs primarily within the Af2 aquifer itself, the largest cumulative deformation occurred in the neighboring aquitards. This counterintuitive pattern is exactly what classical hydro-mechanics predicts. Groundwater withdrawal first reduces pore-water pressure in the high-conductivity aquifer, transmitting pressure changes rapidly and producing only limited elastic deformation. The low-permeability aquitards, by contrast, drain slowly, sustain delayed consolidation, and accumulate disproportionately large inelastic compression. The Ad2 and Ad3 layers therefore contribute far more to the observed subsidence than the aquifer that is actually being pumped, a conclusion that matches Terzaghi’s theory and previous investigations of multilayer aquifer systems.

The broader significance of the study extends well beyond one town in the Yangtze River Delta. Suzhou itself has a long history with subsidence; by the end of 2005, deep groundwater extraction had been completely prohibited in the Suzhou-Wuxi-Changzhou region, and subsidence rates slowed as water levels rose, yet localized zones of concentrated sinking persist. The authors acknowledge limitations, including the absence of descending-orbit radar acquisitions, the fact that TWS is not a direct substitute for groundwater-level measurements, and the reliance on a single representative borehole. Even so, the framework is designed to be general rather than site-specific. Similar combinations of regional surface deformation and subsurface compaction afflict the North China Plain, California’s Central Valley, Jakarta, Mexico City, and the Konya basin in Türkiye, and the researchers argue that pairing regional InSAR with representative borehole fiber optic monitoring offers a practical, scalable strategy for identifying groundwater-induced subsidence and understanding its stratigraphic mechanisms. For water managers confronting the slow-motion disaster of sinking cities, the message is clear: the ground remembers what we pump, and now we can read that memory from orbit and from deep within the earth itself.

Subject of Research: Identification of groundwater-induced land subsidence through integration of satellite InSAR, terrestrial water storage, and borehole distributed fiber optic sensing

Article Title: Regional identification of groundwater-induced land subsidence using integrated surface and subsurface observations

Article References: Sang, H., Fang, K., Liu, Q., Yuan, L., & Shi, B. (2026). Regional identification of groundwater-induced land subsidence using integrated surface and subsurface observations. Environmental Earth Sciences, 85(15), Article 390. https://doi.org/10.1007/s12665-026-13113-x

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13113-x

Keywords: land subsidence, groundwater extraction, InSAR, PS-InSAR, terrestrial water storage, distributed fiber optic sensing, aquitard consolidation, Sentinel-1, BOTDR, Yangtze River Delta, Suzhou, geohazard monitoring

The core problem the researchers identify is structural. For decades, flood risk identification and resilience assessment have been treated as separate exercises, typically published in parallel literatures, using different indicator sets, and feeding into different branches of policy. Risk mapping tells planners where hazard, exposure, and vulnerability converge; resilience assessment tells them how quickly a community or infrastructure network can bounce back after an event. But because the two are rarely analyzed together, management strategies can develop internal contradictions: a city may invest heavily in defenses for high-risk zones while neglecting the recovery capacity of those same zones, or bolster resilience in areas where the hazard itself is comparatively modest. Climate change and rapid urbanization have intensified both the frequency and the severity of flood disasters, and the study argues that this segmented approach now generates structural inconsistencies that undermine regional resilience.

To close that gap, the team developed a full-cycle framework in which flood resilience is embedded directly into flood risk management rather than appended to it. The quantitative engine of the framework is the extension catastrophe progression method, or ECPM, a multi-criteria evaluation technique derived from catastrophe theory. Catastrophe progression methods are well suited to problems where several indicator systems must be combined without arbitrary weighting, because they use the mathematical structure of catastrophe models to aggregate indicators in a standardized way. The extension component widens the set of relationships the method can handle, allowing the researchers to score flood risk and flood resilience across the study region on comparable scales. Validation drew on receiver operating characteristic analysis, with the area under the curve used to test how well the modeled risk surfaces discriminated between locations with and without recorded flood problems.

The indicator architecture is deliberately comprehensive. Flood risk was evaluated across hazard, exposure, and vulnerability dimensions, incorporating variables such as the concentration index of daily precipitation, which captures how violently rainfall is packed into short episodes, along with the concentration index of monthly precipitation, typhoon frequency, distance to rivers, digital elevation model data, land-use type, population density, and gross domestic product density. Resilience was organized around a pressure-state-response logic and drew on concepts from the sustainable livelihoods framework, adding indicators such as road density, railway density, and distance to hospitals to represent the infrastructure and service capacity that determines how quickly a flooded area can be served, evacuated, and rebuilt. All layers were assembled in a geographic information system so that every indicator could be mapped, overlaid, and compared cell by cell across the delta.

The study region is the Yangtze River Delta, an urban agglomeration in eastern China where megacities such as Shanghai, Nanjing, Hangzhou, and Suzhou sit on low-lying ground threaded by rivers and canals and exposed to typhoons arriving from the western Pacific. The region’s eastward slope toward the coast, its extraordinary concentration of population and economic assets, and its history of compound flooding driven by both rainfall and storm surge make it an ideal laboratory for a risk-resilience coupling analysis. It is also a region where urbanization has reshaped the hydrological cycle itself, sealing surfaces, channelizing rivers, and amplifying the intensity of extreme precipitation events, trends documented extensively in prior research on Chinese deltas.

The headline findings are striking. Flood risk in the delta rises from west to east, with high-risk and highest-risk zones together covering 58.0 percent of the region. Resilience displays a broadly similar eastward gradient, with high and highest resilience areas accounting for 73.9 percent of the territory. At first glance that symmetry might look reassuring, but the joint analysis reveals that the match is far from uniform. When the two surfaces are crossed, the largest zoning category is high risk paired with high resilience, covering 35.60 percent of the delta, meaning that more than a third of the region faces serious flood danger but possesses substantial capacity to cope and recover. The most alarming category is the inverse: high-risk, low-resilience areas account for 12.35 percent of the region and are concentrated mainly in highly urbanized districts, where dense built environments, high exposure, and constrained drainage combine to produce danger without adequate defensive depth.

These zoning categories are not merely cartographic curiosities; they translate directly into differentiated prescriptions. In high-risk, high-resilience zones, the priority is to protect and maintain existing coping capacity while monitoring whether intensifying hazards gradually erode it. In high-risk, low-resilience zones, the framework calls for simultaneous risk reduction and resilience enhancement, combining engineered defenses with investments in emergency services, transport redundancy, and social preparedness. Areas of low risk but high resilience can absorb redirected resources with lower urgency, while low-risk, low-resilience zones represent latent vulnerabilities where relatively modest, early investments could prevent future mismatches from forming. The authors frame this as precision flood management, an analogy to precision medicine in which treatment is tailored to the specific profile of each zone rather than applied uniformly across the region.

The study’s methodological contribution lies in showing that the coupling itself carries information that neither risk maps nor resilience maps provide alone. A resilience score of 73.9 percent for the delta sounds impressive until it is laid over a risk surface showing that 58.0 percent of the same territory faces high or highest danger; the residual mismatch, concentrated in exactly the fast-growing urban cores where people and assets pile up, is where the next generation of flood losses is most likely to accumulate. By quantifying the overlap and the divergence, the framework gives governments a diagnostic tool for identifying which districts need protection, which need recovery capacity, and which need both at once, moving flood governance away from one-size-fits-all defense spending toward targeted, evidence-based regulation.

The implications extend well beyond the Yangtze River Delta. Rapidly urbanizing deltas across Asia, Africa, and the Americas face the same collision of intensifying hydro-climatic hazards and constrained adaptive capacity, and many lack any systematic way to see where their risk and resilience profiles have slipped out of alignment. The integrated framework, validated with ROC analysis and grounded in openly available indicator data, offers a replicable template for regional and city-level assessments elsewhere. As extreme precipitation becomes more concentrated and tropical cyclones reach further inland under a warming climate, the study suggests that the most dangerous places will not necessarily be those with the highest flood risk on paper, but those where high risk and low resilience coincide, hidden in plain sight until a coupled analysis makes the divide visible.

Subject of Research: An integrated flood risk and resilience coupling framework applied to the Yangtze River Delta to identify spatial mismatches for precision flood management

Article Title: Bridging the divide: an integrated risk-resilience coupling framework to decode spatial mismatches for precision flood management

Article References: Yu, W., Jiang, M., Wang, X., Yin, L., & Zhang, B. (2026). Bridging the divide: an integrated risk-resilience coupling framework to decode spatial mismatches for precision flood management. Natural Hazards, 122(19), Article 639. https://doi.org/10.1007/s11069-026-08401-5

Image Credits: AI Generated

DOI: 10.1007/s11069-026-08401-5

Keywords: flood risk, flood resilience, Yangtze River Delta, extension catastrophe progression method, spatial mismatch, precision flood management, urbanization, climate change, natural hazards, GIS, disaster risk reduction, resilience assessment