Subsidence Patterns Research Articles

Vertical Displacements and Sea‐Level Changes in Eastern North America Driven by Glacial Isostatic Adjustment: An Ensemble Modeling Approach

AbstractGlacial isostatic adjustment (GIA) describes the response of the solid Earth, oceans, and gravitational field to the spatio‐temporal evolution of ice sheets during a glacial cycle. Present‐day vertical displacements and sea‐level changes vary throughout eastern North America in response to the melting of the Laurentide Ice Sheet following the Last Glacial Maximum. We use the open‐source software SELEN4.0 (a SealEveL EquatioN solver) to investigate the influence of GIA on vertical land motions and sea‐level changes in eastern North America. Further, we evaluate the uncertainties associated with the lithospheric thickness and viscosity structure using an ensemble modeling approach (129,956 total simulations). We identify the best‐fitting rheological profiles by comparing modeled vertical displacements to vertical velocity rates derived from Global Positioning System (GPS). We find a general pattern of subsidence (causing accelerated relative sea‐level rise) in the eastern United States region and uplift (causing relative sea‐level fall) in the eastern Canada region consistent with previous studies for two tested ice sheet models (ICE‐6G(VM5a) and ICE‐7G(VM7)). Overall, we find lower rates of modeled vertical displacement using ICE‐6G(VM5a) compared with ICE‐7G(VM7), which produces lower residuals when compared with the GPS‐derived vertical velocity rates. Our ensemble analysis identifies adjustments to the nominal VM5a and VM7 viscosity models that improve fits to the GPS‐imaged vertical velocity rates throughout eastern North America and on the North American Atlantic Coast. The differences in our best‐fitting models for inland versus coastal regions highlight the importance of exploring lateral viscosity variations for GIA modeling throughout North America and elsewhere.

Controls on the Stratigraphic Architecture of the US Atlantic Margin: Processes Forming the Accommodation Space

AbstractAccommodation space governs the spatial and temporal distributions of sediments in continental margins. Mapping the sedimentation patterns, therefore, offers insights into the solid‐Earth processes that shape accommodation space. We assembled an unprecedented amount of seismic and borehole data along the Eastern North American Margin and used it to divide the margin's sedimentary package into eight chronostratigraphic intervals, identifying temporal shifts in depocenters under the continental shelf, slope, and rise. The Jurassic depocenters follow the syn‐rift structure and its thermal subsidence loci. The Long Island Platform is the only margin segment where the early post‐rift sediment thickness matches subsidence predictions from uniform‐stretching models, whereas in Georges Bank Basin (GBB) and Baltimore Canyon Trough (BCT), sediment thickness is 1.5–3 times higher than predicted, pointing to other factors at play. A margin‐wide Jurassic transient shoulder uplift is inferred from the occurrence of stratigraphic onlaps above thinned crust. Unlike the Jurassic, the Cretaceous and Cenozoic depocenters disregard the inherited subsidence pattern. The accommodation space over the shelf and coastal plain during the Cretaceous was affected by regional isostatic compensation of the sedimentary loads accumulated on the shelf and rise. Accommodation space development in the GBB was interrupted during the Cretaceous after the margin crossed the Great Meteor Hotspot track, resulting in a widespread permanent uplift, erosion, and sediment redistribution. The distribution of anomalous Neogene subsidence in the BCT challenges previous suggestions of mantle dynamic control on the accommodation space and favors flexural downwarping of the shelf by sediment accumulation on the rise.

Assessment of Land Subsidence Pattern in Pokhara Valley Using Sentinel-1 InSAR Processing

This research examines the use of Sentinel-1's enhanced remote method and Interferometric Synthetic Aperture Radar (InSAR) processing to monitor ground subsidence in Pokhara Valley City, Nepal. It aims to comprehend the reciprocal effects of land sinking on environmental sustainability, infrastructure stability, and urban growth. The methodology encompasses the acquisition and processing of Sentinel-1 SAR data, along with rigorous interferometric analysis to detect and quantify ground deformation. Spanning from 2019 to 2022 AD and, -analyzing annual displacement trends, initially, a subsiding trend was observed in 2019-20 AD with an average rate of 0.037 m/year followed by 0.072 m/year in 2020-21 AD, and then escalating to 0.165 m/year in 2021-22 AD. These fluctuations in rates illustrate the intricate interplay of urbanization and geological dynamics affecting land stability in the region. The implications of this study extend to urban planning and infrastructure management, emphasizing the necessity of integrating land subsidence monitoring into the developmental policies and strategies of Pokhara Valley. Additionally, the potential of Sentinel-1 InSAR as a valuable tool for geotechnical investigations offers a non-invasive, cost-effective means to assess and monitor ground deformation risks. By leveraging the capabilities of Sentinel-1 InSAR, the study contributes to the body of knowledge in geotechnical and civil engineering fields, particularly within the context of sustainable urban planning and disaster risk reduction.

Stratigraphic overview of the Neogene Transylvanian Basin, Romania

Although the Transylvanian Basin is often referred as part of the Miocene Pannonian Basin System, we point out important tectonic and sedimentary differences reflected in the stratigraphy. The Neogene of the Transylvanian Basin can be split into two tectonic megasequences: Lower Miocene and Middle to Upper Miocene. The Early Miocene onset of the back arc extension in the Pannonian Basin was not observed in the Transylvanian Basin, where shallow to deep marine sediments were deposited in a flexural basin. Diverse types of facies developed in the continental (e.g., Sânmihai) and marine (e.g., Coruş, Chechiş, Hida) Lower Miocene formations. Although these are poorly constrained to the standard biostratigraphy, recent 87 Sr- 86 Sr data support their limited correlation with the Burdigalian. The Middle Miocene extension in the Transylvanian Basin is very weak compared to that the Pannonian Basin. As the high-resolution Middle Miocene regional biostratigraphy of the Badenian and Sarmatian was calibrated by magnetostratigraphy and radiometric dating of the volcanic tuffs up to the Serravallian/Tortonian boundary, the influence of the regional tectonics on the basin's evolution and stratigraphy can be better understood. The regional tectonics produced a progressive basin isolation starting with the late part of the Middle Miocene. The uplift of the Carpathians isolated the basin at the beginning of the Pannonian and initiated the fan delta, lacustrine, and fluvial sedimentation. The basin was rapidly uplifted and eroded, while post-extensional subsidence was still ongoing in the Pannonian Basin. Endemic fossil assemblages allow only regional-scale correlations for this interval. The uppermost part of the stratigraphic record is represented by the volcano-sedimentary formations generated by the late evolution of the Eastern Carpathians. The sedimentation and subsidence pattern of the Transylvanian Basin reflect the effect of retreating subduction zone in the Eastern Carpathians.

Depositional, stratigraphic and structural controls on the geometry of Montney Formation reservoirs: Lower Triassic, northeast British Columbia, northwest Alberta

Abstract The Lower Triassic Montney Formation is a world-class unconventional hydrocarbon reservoir that spans the border of northeast British Columbia and northwest Alberta. Integrated sedimentological and ichnological analyses suggest that the Montney Formation was deposited in a predominantly prodeltaic setting that was commonly influenced by both river flood and storm processes. Riverine processes include the local occurrence of normally- and inversely-graded siltstone and sandstone interbeds, interpreted as the record of sediment-laden hyperpycnal flows. Episodic storms are recorded by the local occurrence of hummocky cross-stratified interbeds, as well as normally graded, parallel and wave-ripple laminated beds, which are commonly draped by structureless carbonaceous mudstones. These mudstone drapes represent river-flood-derived fluid muds that collapsed from hypopycnal plumes following the passage of the storms and suggest the intimate linkage between storms and heightened fluvial discharge. The basin in which the Montney Formation was deposited began as three sub-basins, northern, central and southern, each of which experienced different patterns of subsidence and uplift that varied both spatially and temporally during sedimentation. Subsidence of the northern sub-basin was compartmentalized by deep seated faults that provided accommodation for the basal sequence of the Montney Formation, Sequence 1. Subsequent fault movement resulted in reduced accommodation in the northern sub-basin and accelerated subsidence in the central sub-basin, providing accommodation for deposition of Sequence 2, the basal sequence in the central sub-basin. The northern sub-basin remained paleotopographically elevated throughout deposition of Sequence 2 and early Sequence 3 successions, and these younger units did not extend beyond the northern limit of the central sub-basin. During subsequent deposition of the late Sequences 3 and Sequence 4 successions, compartmentalization between the north and central sub-basin was greatly reduced and these intervals were deposited relatively more uniformly along-strike. The unconventional and conventional reservoirs in both the northern and central sub-basins occur at the scale of systems tracts within the four depositional sequences of the Montney Formation. Reconstructing the depositional history of the component systems tracts of Sequence 1 through early Sequence 3 reveals the dynamic interplay between syndepositional tectonism, fluctuating relative sea-level and deltaic depositional processes in controlling the character, geometry and distribution of prolific Montney Formation reservoirs. The recognition of these key geological controls enhances predictability in this valuable resource and may be applicable to other fine-grained reservoirs elsewhere.

Tectonic evolution of the eastern margin of the Northern South Yellow Sea Basin in the Yellow Sea since the Late Cretaceous

The South Yellow Sea Basin is the largest sedimentary basin in the Yellow Sea. The Gunsan Basin in the central eastern part of the Northern South Yellow Sea Basin comprises the Western, Central, and Eastern Subbasins. The Eastern Subbasin marks the eastern margin of the South Yellow Sea Basin. Interpretation of multi-channel seismic profiles and balanced cross-section restoration of depth-converted seismic profiles reveal the structural characteristics and evolution of the Eastern Subbasin. The subbasin comprises three groups of faults: NW-SE, NE-SW, and NNE-SSW trending faults. The subsidence pattern of the subbasin, derived from the cross-section restoration analysis, indicates five distinguished evolution phases: (i) rapid subsidence in the Late Cretaceous, (ii) slow subsidence from the Paleocene to the Eocene, (iii) accelerated subsidence in the Oligocene, (iv) alternation of the uplift and subsidence in the Early Miocene, and (v) gradual subsidence since the Middle Miocene. The main and moderate subsidence can be explained by the combination of extension in the SE and NE-SW directions that formed double duplex structures. We suggest that the NW oblique subduction of the Pacific Plate under the NE Asia margin induced both extension toward the trench and dextral strike-slip parallel to the margin. The extension toward the trench caused SE transtension in a local pull-apart setting, whereas the dextral strike-slip parallel to the margin caused NE-SW extension, inducing more significant subsidence. The combined processes resulted in the progressive development of nested duplex structures. The evolution of the Eastern Subbasin is not compatible with previously suggested models for the western part of the South Yellow Sea Basin including foreland basin formation and transtension induced by branch faults of the Tan-Lu Fault.

Understanding land subsidence in the Pearl River Delta region of China based on InSAR observations

Land subsidence occurs in river delta areas. A comprehensive understanding of the factors contributing to land subsidence is crucial for effective geohazard mitigation and land planning. This study investigates the distribution characteristics and spatiotemporal evolution pattern of the land subsidence in the Pearl River Delta based on the InSAR-derived ground deformation data spanning 2016 to 2021; and, the effects of the influencing factors, including the Quaternary sediment thickness, land use classification, land reclamation, and subway construction and operation, on the subsidence are analyzed. The results unveil the presence of four major subsidence zones and several minor coastal subsidence bowls in the research area, with a total area of about 1100 km, all of which exhibit prolonged and intense subsidence. In the remaining regions, land subsidence predominantly exhibits weak and short-term occurrences. The subsidence in the study area follows a distinct linear evolutionary pattern with a slight rebound, potentially linked to the rise of groundwater levels. Regions with Quaternary sediment thickness surpassing 20 m, engaging in aquaculture or agricultural activities, and featuring reclaimed land are predisposed to prolonged and intense ground subsidence. In addition, the construction and operation of subway lines may induce ground subsidence, especially in geologically vulnerable areas. This study provides valuable insights into land subsidence research in delta regions.

3D deformation velocity field analysis and TEM method to detect the tectonic influence on the land subsidence in Zamora, Mexico

The city of Zamora is situated in a tectonic basin and, according to previous interferometric studies, has been experiencing sinking up to 13 cm/yr (2007-2011). Although the reported subsidence pattern (WNW-ESE) is similar to the orientation of the regional fault system, there is a lack of detailed studies to establish the connection between this phenomenon, the thickness of sediments, and the basement geometry. Therefore, this paper presents a 3D deformation velocity field analysis for the period of 2014–2021, using 130 images from the Sentinel-2 satellite to calculate the vertical and east-west components of the subsidence. Likewise, a campaign of 28 Transient Electromagnetic soundings was carried out around the city to determine the thickness of sediments and bedrock geometry. The results of the interferometric analysis reveal a WNW-ESE sinking pattern with the maximum vertical velocities observed near the Zamora Museum, amounting to 108 mm/yr (2014–2017) and 102 mm/yr (2018–2021). Meanwhile, the horizontal component has a heterogeneous behavior in both periods with rates reaching up to 1 cm/year. The Transient Electromagnetic survey results show the presence of two grabens, a horst, and a series of half-grabens. The thickness of sediments varies between 94 and 314 m, with the latter corresponding to the main graben located beneath the city. The axis of this structure and the largest sediment deposits are linked to the area of maximum sinking and the spatial pattern of subsidence.

Spatiotemporal characterization of the subsidence and change detection in Tehran plain (Iran) using InSAR observations and Landsat 8 satellite imagery

Urban areas worldwide are increasingly facing challenges related to land subsidence, a phenomenon exacerbated by uncontrolled groundwater extraction and urban expansion. This research focuses on the Tehran plain, Iran's capital city, where significant subsidence has been observed due to uncontrolled migrations influenced by various economic and political factors. This expansion has increased demand for energy, notably water, leading to irregular water withdrawals from underground sources and, consequently, land subsidence. Monitoring this subsidence, particularly its effects on urban infrastructure, has become a critical challenge. This research first reviewed the existing body of knowledge related to subsidence measurement in the Tehran plain with an emphasis on their findings and limitations and then used radar images to study the subsidence patterns in the Tehran plain from 2016 to the end of 2020. Finally, the results collaborated by optical imagery analysis to find the relationship between surface change detection and spatiotemporal distribution of subsidence. As a result, through processing Sentinel-1A SAR images, consistent vertical displacements (subsidence) were observed, especially in areas heavily reliant on groundwater from wells, with some areas experiencing a rate of more than −20 mm/year. Horizontal displacement, however, was approximately about ±8 mm/year. Also, our results show that the subsidence rate in this plain has decreased in recent years. Therefore, the study integrated multispectral satellite data to clarify this issue and compensate for missing groundwater level data, specifically the Normalized-Difference Vegetation Index (NDVI) and Normalized-Difference Moisture Index (NDMI). These datasets were used to monitor changes in vegetation cover distribution and moisture in response to the variations of groundwater depth over time. The results of this research can be beneficial in adequately managing groundwater resource utilization to reduce the potential damage to infrastructure and the environment.

Asymmetric analogous hyperbola model of overburden movement and its verification

The extraction of underground coal resources induces the fracture and movement of overlying strata, leading to geological hazards such as surface deformation, cracks, and even subsidence. Utilizing the analogous hyperbola model of overlying strata movement, we conducted a mechanical analysis to examine the asymmetric fracture mechanism resulting from coal seam mining in thick loose strata. An asymmetric analogous hyperbola model was established by introducing distinct virtual half-axis lengths (b). The thickness impact of thick loose layers (H) and bedrock layer (h) on the asymmetric movement of overlying rock during mining was also discussed. Similarity model tests were conducted to research the migration characteristics and surface subsidence patterns of overburdened rock and thick loose layers at different mining stages and validate the hypothesis of asymmetric overburdened rock migration. Additionally, the discrete element numerical model for thick and loose layers mining was established by using UDEC and discussed the asymmetric analogous hyperbola behaviour of overburden movement and surface subsidence. The comparison results show that the established asymmetric hyperbolic model can effectively predict the movement law of overlying strata and surface subsidence characteristics. Therefore, the proposed model can provide valuable theoretical support for predicting the movement patterns of overburden under thick loose layers and mitigating surface subsidence disasters.

Editorial: Towards a better understanding of the correlation between the subsidence pattern and land use type

Editorial: Towards a better understanding of the correlation between the subsidence pattern and land use type

Detectability of low-viscosity zone along lithosphere–asthenosphere boundary beneath the Nankai Trough, Japan, based on high-fidelity viscoelastic simulation

After the 2011 Tohoku-oki earthquake, a seafloor observation system observed rapid landward deformation. These seafloor observation data are potentially explicable via modeling of a low-viscosity zone under the subducting plate, called the lithosphere–asthenosphere boundary (LAB). However, the effect of a low-viscosity LAB has not been empirically examined in past earthquakes because of the absence of seafloor observation data, which are expected to have high sensitivity in the detection of a low-viscosity LAB. The Nankai Trough, southwest Japan, is one location where seafloor geodetic stations are in operation that could detect a low-viscosity LAB after a future earthquake. Therefore, we quantitatively model how the low-viscosity of the LAB is expected to affect postseismic crustal deformation following an anticipated future large earthquake under the Nankai Trough. Calculating viscoelastic deformation using a model that takes into account the realistic crustal structure in the southwestern Japanese Archipelago, we examine the effects of LAB viscosity on surface deformation patterns. The results show for very low LAB viscosities (2.5×1017\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2.5 \ imes 10^{17}$$\\end{document} Pa s), rapid landward displacement as large as 38 cm/yr occurs in the offshore area, with uplift and subsidence patterns highly dependent on the viscosity structure. Especially in the first year after the earthquake, the difference in deformation between the cases with and without a low-viscosity LAB was much larger than the observational error level of the current seafloor observation system. For such low LAB viscosities, the probability of detection is therefore high. On the other hand, for moderately low LAB viscosities (2.5×1018\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2.5 \ imes 10^{18}$$\\end{document} Pa s), the deformation patterns in the cases with and without a low-viscosity LAB are similar. Indeed, the difference in displacement lies within the observational error level of the current observational network, so detection is expected to be difficult for a LAB of only moderately low-viscosity. However, such a LAB may be detected by improving observational accuracy in the Global Navigation Satellite System–acoustic observation technique by performing more frequent measurements. Additional detection potential lies in expected future accuracy improvements in vertical displacement estimates at the DONET stations.Graphical

Investigating the Mechanism of Land Subsidence Due to Water Network Integration at the Guangzhou Longgui Salt Mine and Its Impact on Adjacent Subway

Water-soluble mining was invariably associated with surface subsidence, which in some cases escalated to the movement, deformation, and even collapse of the overlying rock layers, triggering grave subsidence calamities. The caprock of the salt-bearing strata in the Longgui salt rock mining area was closely adjacent to the third aquifer, which mainly consisted of fractured, porous, high-permeability materials such as mudstone conglomerates, rendering the geological conditions highly complex. Years of water-soluble mining had led to significant surface subsidence in the mining area, with a trend toward accelerated subsidence. In this study, the geological conditions of the Longgui salt rock mining area were analyzed, and through simulated experiments of pillar dissolution mining, the mechanisms of surface subsidence in the area were examined. Over time, the dissolution gradually perforated the pillars and caprock, with the pillars ceasing to support the caprock, ultimately transforming small cavities into a large single cavity. Utilizing subsidence data, this research employed numerical simulation to inverse and predict subsidence patterns from 2019 to 2025, revealing that the maximum subsidence reached 1367.6 mm in mining area I and 1879.5 mm in mining area II, with subsidence rates of 12.05 mm/y and 44.78 mm/y, respectively. Moreover, the impact of ground subsidence on the construction of adjacent subways was assessed by establishing monitoring points and evaluating subsidence along subway cross-sections and longitudinal directions. The findings provided valuable insights for guiding the prevention and control of surface subsidence calamities in the Longgui salt rock mine and similar mining areas in Guangzhou, China.

Exploring groundwater depletion and land subsidence dynamics in Taiwan’s Choushui river alluvial fan: insights from integrated GNSS and hydrogeological data analysis

The Choushui River Alluvial Fan (CRAF) is a major agricultural area in Taiwan with heavy groundwater usage. The extraction of groundwater here has caused land subsidence, which is now a significant global environmental issue. This study analyzes land subsidence in the CRAF by integrating hydrogeological data from 233 groundwater monitoring stations across four aquifers (CRAF Groundwater_NET) and data from 50 continuous GNSS stations (CRAF GNSS_NET). We developed an automated processing flow for GNSS static surveying within CRAF GNSS_NET, and further employed a time-series fitting method to examine the long-term trends and annual changes for both GNSS and groundwater level data. Our analysis of the time-series data from the past decade identifies areas of significant groundwater level depletion and subsidence hotspots. We explore the relationship between groundwater level variations and surface displacements within CRAF, utilizing GNSS data to analyze horizontal and vertical displacement trends, as well as annual changes. We integrate these findings with hydrogeological data to understand regional subsidence patterns. Our results indicate that CRAF is characterized by distinct hydrogeological features. The study finds that the amplitudes of annual changes in both groundwater level and vertical displacement generally increase from northeast to southwest in the analyzed region. One particular area shows significant groundwater level decline, with the most severe rate recorded at 0.54 m/year. Similarly, GNSS analysis indicates pronounced subsidence trends in the same area, with rates ranging from 4.2 to 5.2 cm/year. These findings highlight the critical need for the development of effective groundwater management strategies to ensure sustainable use of groundwater resources and to implement mitigation measures against land subsidence in similar multiple-aquifer settings.

Testing Mantle Convection Simulations With Paleobiology and Other Stratigraphic Observations: Examples From Western North America

AbstractMantle convection plays a fundamental role in driving evolution of oceanic and continental lithosphere. In turn it impacts a broad suite of processes operating at or close to Earth's surface including landscape evolution, glacio‐eustasy, magmatism, and climate. A variety of theoretical approaches now exist to simulate mantle convection. Outputs from such simulations are being used to parameterize models of landscape evolution and basin formation. However, the substantial body of existing simulations has generated a variety of conflicting views on the history of dynamic topography, its evolution and key parameters for modeling mantle flow. The focus of this study is on developing strategies to use large‐scale quantitative stratigraphic observations to assess model predictions and identify simulation parameters that generate realistic predictions of Earth surface evolution. Spot measurements of uplift or subsidence provide useful target observations for models of dynamic topography, but finding areas where tectonics have not also influenced vertical motions is challenging. To address this issue, we use large inventories of stratigraphic data from across North America with contextual geophysical and geodetic data to constrain the regional uplift and subsidence history. We demonstrate that a suite of typical geodynamic simulations struggle to match the amplitude, polarity and timing of observed vertical motions. Building on recent seismological advances, we then explore strategies for understanding patterns of continental uplift and subsidence that incorporate (and test) predicted evolution of the lithosphere, asthenosphere and deep mantle. Our results demonstrate the importance of contributions from the uppermost mantle in driving vertical motions of continental interiors.

Quantifying land subsidence and its nexus with groundwater depletion in isfahan-borkhar plain: An integrated approach using radar interferometry and spatial bivariate relationships

Land subsidence has emerged as a significant global concern, particularly in arid regions like Iran. Among these vulnerable areas, the Isfahan-Borkhar plain stands out due to its strategic importance and critical high-population settlements. In this study, we focus specifically on land subsidence within this plain. This study specifically examines land subsidence in this plain. Utilizing the precise radar interferometry technique, considered the most accurate method, this study investigates the vertical changes of the Earth's surface in the Isfahan-Borkhar plain between 2015 and 2022 by analyzing Sentinel-1 satellite data. To establish causality, we correlate these subsidence patterns with spatial water extraction data, employing a spatial bivariate relationship model. The results obtained from this seven-year monitoring period reveal that average annual subsidence in different areas of the Isfahan-Borkhar plain ranges from 25 to 45 cm, while the average annual decline in the water table in this aquifer is 230 cm. The pronounced subsidence observed aligns with areas experiencing the highest water level decline, emphasizing the critical role of water table depletion. In summary, our study underscores the need for sustainable water management practices to mitigate land subsidence risks in this strategically vital region.

Triggering of Land Subsidence in and Surrounding the Hangjiahu Plain Based on Interferometric Synthetic Aperture Radar Monitoring

In the early stages, uncontrolled groundwater extraction led to the Hangjiahu (HJH) Plain becoming one of the areas with the most severe land subsidence in China. Since the beginning of this century, comprehensive measures have been taken to control the continuous aggravation of large land subsidence patterns in some areas; however, urban land subsidence issues, influenced by various factors, still persist and exhibit complex geographical distribution characteristics. In this study, we utilized Sentinel-1A images and the SBAS-InSAR technique to capture surface deformation over the HJH Plain in Zhejiang from 16 March 2017 to 20 January 2023. Through a comparative analysis with geological conditions, changes in surface mass loading, rainfall and groundwater, and land use types, we discussed the contributions of natural and anthropogenic factors to land subsidence. Augmented with optical remote sensing images and field investigations, we conducted a correlation analysis of the land subsidence status. The preliminary findings suggest that changes in surface mass loading and short-term heavy rainfall under extreme weather conditions can lead to periodic land subsidence changes in the region. Additionally, extensive infrastructure construction triggered by urbanization has resulted in significant and sustained land subsidence deformation. The research findings play an important guiding role in formulating scientifically effective strategies for land subsidence prevention and control, as well as urban planning and construction.

Spatio-temporal evolution of land subsidence and susceptibility zonation of associated ground fissures in the urban area of loess plateau: a case in Xianyang city, China

For over 50 years, the Xianyang city on the Chinese loess plateau, has been deeply affected by land subsidence and its associated ground fissures. In this study, 67 Sentinel-1A images from 2015 to 2022 were analysed to obtain the spatio-temporal evolution of land subsidence. Subsidence centers identified in Yunyang town and Luqiao town exhibit annual subsidence rates of −31 mm/a and −26.7 mm/a, respectively, culminating in total subsidence of −258 mm in Yunyang town and −139 mm in Luqiao town. Moreover, time series analysis revealed that both towns exhibit pronounced seasonal subsidence patterns. It is believed that human cultivation and groundwater overuse, combined with the effects of active faults, have led to uneven land subsidence. Then, the gradient of land subsidence and its potential relation to ground fissure hazards were evaluated. The results highlight a strong spatial correlation between high subsidence gradient areas and dense ground fissure distribution. Therefore, we proposed for the first time to introduce the land subsidence gradient factors for the susceptibility mapping of ground fissures using the artificial neural network algorithm. The susceptibility zonation results show that integrating the gradient factor can improve the rationality of the susceptibility evaluation. The above understanding provides valuable support for disaster prevention and mitigation in urban areas affected by land subsidence and ground fissures in loess regions.

Unravelling the influencing hydrogeological factors contributing to land subsidence in the Tianjin Plain of China using a multi-scale geographically weighted regression model and monitoring data

The establishment of a quantitative relationship between land subsidence and its influencing factors is a crucial task in developing prevention strategies of land subsidence in specific areas. In this study we examined the dynamic patterns of land subsidence before and after establishing the South-to-North Water Division Project (SNWDP) in the Tianjin Plain (TJP) that aims to reduce groundwater extractions in certain areas. Statistical analyses were used to determine the key attributing factors contributing to land subsidence that were subsequently used to derive the so-called ‘critical water level of land subsidence’(CWLS), which were used to manage groundwater extractions with the objective of minimizing further land subsidence. The main results obtained were as follows: under the influence of the SNWDP, the groundwater level in the TJP stopped declining and the water level rose in most areas. The area of ‘strong’ subsidence, defined as the area with a subsidence rate greater than 50 mm a −1 , decreased from 16.0% before the SNWDP to 2.5% after. The results of the multi-scale geographically weighted regression (MGWR) model revealed that groundwater drawdown in the second and third confined aquifer was the main contributor to the areas classified as ‘strong’ and ‘medium’ land subsidence areas. Following the reduction in groundwater extractions, the behaviour of surface deformation can be grouped into two categories: first, a rebound in surface elevation; and, second, continuous compression with a lower subsidence rate. The first mode occurs mainly in the area of ‘medium’ and ‘weak’ subsidence, and the second mode occurs mainly in the ‘strong’ subsidence area. Supplementary material : The result of the MGWR model is available at https://doi.org/10.6084/m9.figshare.c.7103611

Multi-Scale Analysis of Surface Building Density and Land Subsidence Using a Combination of Wavelet Transform and Spatial Autocorrelation in the Plains of Beijing

Land subsidence is a major issue in the Beijing Plain in China, caused by the construction of new buildings and infrastructure combined with groundwater extraction. This study employs a multi-level two-dimensional wavelet decomposition to decompose land subsidence into high- and low-frequency components, and Moran’s I index to analyze the spatial distribution of land subsidence and its main influencing factors. By comparing the spatial distributions of the high- and low-frequency components, we estimate the correlation between land subsidence and influencing factors at different scales. Utilizing a combination of wavelet decomposition and Moran’s I analysis, our study establishes a clear spatial correlation between continuously varying factors such as groundwater and clay layer thickness, and the low-frequency components of land subsidence, allowing for a focused analysis of the relationship between surface building density and the high-frequency components of land subsidence. Quantitatively, the study identifies a significant correlation at specific granularities, particularly at 480 m and 960 m, underscoring the nuanced interaction between urban development and land subsidence patterns. These insights into the spatial distribution of land subsidence and its contributing factors can inform the development of effective strategies to address this issue.