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

Human papillomavirus is the most common sexually transmitted infection worldwide and the necessary cause of nearly all cervical cancer cases. Yet the virus is far from uniform: its genotype distribution varies markedly between countries, regions, and even cities, and those differences matter enormously for vaccination policy. Prophylactic HPV vaccines target specific subsets of genotypes, so knowing which types actually prevail in a given population determines how much protection a vaccination program can realistically deliver. Because large-scale data on genotype prevalence and viral load characteristics in the Taihu region had been lacking, the research team set out to fill that gap with a three-year survey spanning January 2023 to December 2025.

The study design was straightforward but ambitious. The investigators analyzed cervical exfoliated cell samples collected from 23,196 women who underwent HPV testing at their institution during the study window. Each sample was screened for 21 HPV genotypes using real-time fluorescence quantitative polymerase chain reaction, a technique that both identifies the viral type present and quantifies how much viral DNA is in the specimen. This dual capability allowed the team to go beyond simple positivity rates and construct viral load profiles across age groups, a dimension that many prevalence surveys omit.

The headline finding was an overall HPV infection rate of 13.20 percent, with 3,061 of the 23,196 participants testing positive. Notably, this rate remained statistically stable across the three years of the study, suggesting an entrenched and steady level of viral circulation in the region rather than a rising or falling epidemic curve. That stability provides a useful baseline against which the future impact of vaccination and screening programs can be measured.

Perhaps the most striking pattern emerged when the researchers stratified infections by age. Rather than declining monotonically with age, as some classical models of HPV epidemiology predict, prevalence followed a bimodal distribution. The highest infection rate appeared in women aged 60 and older, while the lowest rate occurred among women aged 40 to 49. This shape echoes findings from other Chinese regions and reflects a complex interplay of factors: cohort effects in lifetime exposure, changes in immune competence with age, hormonal shifts after menopause, and possibly reduced participation in screening among older women, allowing persistent infections to accumulate undetected.

The genotype ranking carried direct implications for vaccine policy. Among high-risk types, HPV52 led with a prevalence of 2.70 percent, followed by HPV58 at 1.72 percent, HPV16 at 1.35 percent, HPV39 at 1.11 percent, and HPV56 at 0.96 percent. Among low-risk genotypes, HPV53 was the most common, detected at 1.78 percent. HPV16, the single most carcinogenic type globally, ranked only third locally, while HPV52 and HPV58, types that are notably overrepresented in Asian populations relative to Western ones, dominated. This regional signature matters because the bivalent and quadrivalent vaccines, which target HPV16 and HPV18 (with the quadrivalent also covering HPV6 and HPV11), would leave a substantial fraction of local infections uncovered. The nonavalent vaccine, which adds HPV31, HPV33, HPV45, HPV52, and HPV58 to the mix, aligns far better with the genotype spectrum observed in Taihu, and the authors concluded that their data show greater theoretical concordance with the nine-valent formulation.

Age also shaped the composition of infections in ways that compound the vulnerability of older women. Women aged 60 and above exhibited significantly higher proportions of high-risk mixed infections and multiple-genotype infections compared with all other age groups. Multiple infections, in which two or more HPV types coexist in the same cervical sample, are of particular clinical interest because co-infection can complicate risk stratification and may reflect broader exposure or diminished viral clearance. The elevated burden of mixed high-risk infections in older women underscores that HPV is not simply a concern of young adulthood and that screening programs must not neglect postmenopausal populations.

The viral load analysis added a further layer of nuance. Viral load, typically measured as the quantity of viral DNA per cellular sample, is thought to correlate with infection persistence and lesion risk, although its clinical interpretation remains debated. In this cohort, the distribution of viral loads varied systematically with age. High viral loads were most frequent among women under 30, accounting for 24.29 percent of infections in that group, whereas low viral loads predominated among women aged 60 and older, making up 46.82 percent of their infections. One plausible interpretation is that younger women carry more recently acquired, actively replicating infections that are often cleared by the immune system, while older women more frequently harbor long-standing, low-level persistent infections whose cumulative cellular damage may be clinically significant despite modest viral quantities.

Taken together, the results sketch a coherent picture of HPV epidemiology in the Taihu region: a stable overall prevalence near 13 percent, a genotype hierarchy led by HPV52, HPV58, and HPV16, a bimodal age distribution peaking in the oldest women, and an age-dependent viral load architecture. The authors argue that these findings argue for region-specific prevention and control strategies rather than one-size-fits-all policies. In practical terms, that means prioritizing vaccines whose genotype coverage matches local circulating types, sustaining and possibly extending cervical screening into older age brackets, and considering viral load and multiple-infection status when triaging positive results in clinical practice.

The study also demonstrates the value of quantitative PCR-based surveillance as a public health instrument. By combining genotyping with viral load quantification across a very large sample, the Taihu dataset offers a template for other regions in China and beyond that lack local HPV epidemiology. As vaccination programs scale up globally, longitudinal studies of this kind will be essential for tracking shifts in genotype prevalence, detecting vaccine-type declines, and identifying any replacement by non-vaccine types. For now, the message from Suzhou is clear: the local viral landscape favors the broadest available vaccine coverage, and older women deserve renewed attention in the campaign against cervical cancer.

Subject of Research: Epidemiology, genotype distribution, and viral load profiles of human papillomavirus infection among women in the Taihu region of Suzhou, China

Article Title: Epidemiological characteristics, genotype distribution, and viral load profiles of human papillomavirus infection among women in the Taihu region of Suzhou: a retrospective study

Article References: Epidemiological characteristics, genotype distribution, and viral load profiles of human papillomavirus infection among women in the Taihu region of Suzhou: a retrospective study. (n.d.). https://doi.org/10.1186/s12985-026-03301-3

Image Credits: AI Generated

DOI: 10.1186/s12985-026-03301-3

Keywords: human papillomavirus, HPV prevalence, HPV52, HPV58, HPV16, genotype distribution, viral load, cervical cancer, nonavalent vaccine, Suzhou, retrospective study, qPCR