Subsidence Bowl Research Articles

Unveiling subsidence patterns: Time series analysis for land deformation investigation in the west-Qurna oil field, Iraq

Land subsidence is a worldwide geological and environmental risk caused by natural occurrences and human actions. Its effects include a range of socio-economic, environmental, and hydrogeological consequences, such as damage to infrastructure like buildings, roads, bridges, and pipelines, as well as increased flooding and reduced groundwater storage capacity. Due to these diverse impacts, it is crucial to monitor the spatial and temporal scope of land subsidence. This study presents an investigation into land subsidence within the West-Qurna oil field, a large oil reservoir situated in Iraq's Basrah governorate. The study employs Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) analysis from the European Space Agency Sentinel 1A over six years, from June 2017 to May 2023. The Stanford Method for Persistent Scatterers (StaMPS) has been utilised to assess the scale and magnitude of land deformations in this region. Results revealed a notable subsidence within the central urban area of the oil field, forming an ellipsoidal subsidence bowl spanning 86 km2. The peak subsidence rate is identified at −13.2 ± 0.4 mm/yr within this bowl, with a cumulative vertical displacement of 75 mm throughout the six-year observation period. Furthermore, uplifting phenomena are also detected at the study area's peripheries, reaching a maximum rate of 12 ± 0.4 mm/yr and a cumulative shift of 54 mm. Temporal analysis showcases a significant alteration in subsidence rates, with rates of −18 mm/yr observed between 2017 and 2020, followed by −5 mm/yr post-2020. This change is attributed to COVID-19-related oil production reductions enacted by the government to boost prices. Our analysis points toward oil extraction as a probable primary driver of subsidence in the studied area, although a deeper probe into the impact of groundwater extraction for reservoir injection remains essential.

Quantitative Evaluations of Pumping-Induced Land Subsidence and Mitigation Strategies by Integrated Remote Sensing and Site-Specific Hydrogeological Observations

Land subsidence is an environmental hazard occurring gradually over time, potentially posing significant threats to the structural stability of civilian buildings and essential infrastructures. This study presented a workflow using the SBAS-PSInSAR approach to analyze surface deformation in the Choushui River Fluvial Plain (CRFP) based on Sentinel-1 SAR images and validated against precise leveling. Integrating the InSAR results with hydrogeological data, such as groundwater levels (GWLS), multilayer compactions, and borehole loggings, a straightforward model was proposed to estimate appropriate groundwater level drops to minimize further subsidence. The results showed a huge subsidence bowl centered in Yunlin, with maximal sinking at an average 60 mm/year rate. High-resolution subsidence maps enable the quantitative analyses of safety issues for Taiwan High-Speed Rail (THSR) across the areas with considerable subsidence. In addition, the analysis of hydrogeological data revealed that half of the major compaction in the study area occurred at shallow depths that mainly included the first and second aquifers. Based on a maximal subsidence control rate of 40 mm/year specified in the CRFP, the model results indicated that the groundwater level drops from wet to dry seasons needed to be maintained from 3 to 5 m for the shallowest aquifer and 4–6 m for Aquifers 3 and 4. The workflow demonstrated the compatibility of InSAR with traditional geodetic methods and the effectiveness of integrating multiple data sources to assess the complex nature of land subsidence in the CRFP.

InSAR data for detection and modelling of overexploitation-induced subsidence: application in the industrial area of Prato (Italy)

Spaceborne-based monitoring for environmental purposes has become a well-established practice. The recent progress of synthetic aperture radar (SAR) sensors, including through the European Space Agency’s (ESA) Sentinel-1 constellation, has enabled the scientific community to identify and monitor several geohazards, including subsidence ground deformations. A case study in the Tuscany Region, Italy, highlights the effectiveness of interferometric synthetic aperture radar (InSAR) in detecting abrupt increases in ground deformation rates in an industrial area of Montemurlo municipality. In this case, InSAR data enabled prompt identification of the phenomenon, supporting the authorities in charge of environmental management to thoroughly investigate the situation. First, an on-site validation was performed via field surveys confirming the presence of cracks and fissures on some edifices. Further analysis, including water pumping rates, settlement gauge and topographic levelling, corroborated the InSAR data's findings regarding vertical deformation. Integration of collected data allowed for spatial identification and assessment of the subsidence bowl and its source depth recognized by the remote sensing data. The Montemurlo case offers a procedural guideline for managing abrupt accelerations, identified by InSAR data in subsidence-prone areas due to fluid overexploitation. In fact, these data proved useful in helping local authorities responsible for hydrogeomorphological risk management. With the exacerbation of deformation issues in subsidence-prone regions due to climate change, early detection and monitoring of such phenomena are increasingly crucial, with InSAR data playing a central role in achieving this goal.

Three-dimensional surface deformation from multi-track InSAR and oil reservoir characterization: A case study in the Liaohe Oilfield, northeast China

Oil productions can result in pore pressure drop in the reservoir, generating an increase in effective stress and leading to reservoir compaction. The compaction in the subsurface reservoir translates to the earth's surface, which is manifested as a loss of elevation (land subsidence), causing damages to oil production facilities and surface infrastructures. The Liaohe oilfield, located in Liaohe River Delta (LRD), northeast of China, is one of the most significant subsidence areas in China as a direct consequence of oil extraction from the reservoir. Previous studies carried out in this area assumed the oil production-induced displacement retrieved from Interferometric Synthetic Aperture Radar (InSAR) corresponds only to vertical deformation. In this work, for the first time, we proposed a method to retrieve the full three-dimensional (3D) displacement field over the oilfield. We retrieved the vertical and east-west displacement components by combining the multiple-geometry InSAR line-of-sight (LOS) observations and retrieved the north-south component based on the assumption of a physical relationship between the horizontal and vertical displacement. Two ascending and two descending datasets from Sentinel-1 satellite covering the area were processed by an InSAR time series analysis over the 2017 to 2021 period, providing consistent displacement rate maps and displacement time series in the LOS direction. Spatial local-scale land subsidence was found in several producing fields over the deltaic region, including Shuguang, Huanxiling, and Jinzhou oilfields. The 3D displacement decomposition was then conducted in Shuguang oilfield. The derived 3D displacement field exhibit a circular subsidence bowl with a maximum subsiding rate reaching 212 mm/year, accompanied by a centripetal pattern of horizontal displacements with maximum rates up to 50–60 mm/year moving towards the subsidence center. The retrieved 3D displacements are in good agreement with predictions from the geomechanical modeling by assuming a disk-shaped reservoir subject to a uniform reduction in pore fluid pressure. Finally, we show the importance of knowing both the vertical and horizontal displacement in characterizing the lateral boundary of the subsurface reservoir.

InSAR time series analysis of natural and anthropogenic coastal plain subsidence: A case of Hangjiahu plain

The early excessive groundwater exploitation caused extensive land subsidence in the Hangjiahu Plain (HJHP), threatening geological stability and infrastructure. Government restrictions on groundwater exploitation led to considerable changes in surface subsidence. To explore the distribution and causes of surface deformation after restrictions on groundwater exploitation, we used 46 Sentinel-1A TOPS images to detect the spatiotemporal distribution of surface deformation in the HJHP area from January 2019 to October 2022 by SBAS-InSAR technique and verified it with leveling results. Results indicated the following: (1) The uneven surface subsidence in the HJHP was relatively apparent, with the localized subsidence concentrated in a few areas. During InSAR monitoring, the maximum subsidence velocity reached −21.0 mm/yr. The most severe subsidence area was in the Shengze town industrial park, with remarkable subsidence bowls and a maximum subsidence of −85.8 mm. (2) Results were verified by comparing the InSAR results with leveling, and the results were in good agreement, indicating the reliability of the InSAR results. (3) Comprehensive analysis of the spatial distribution characteristics of surface subsidence, geological data and rainfall changes revealed that fault activity and rainfall changes had varying degrees of effect on surface subsidence in the HJHP area. (4) The comparison with land use data indicated that the built area was the dominant mechanism of surface subsidence in the study area. In addition, the effect of industrial activities on the surface cannot be ignored. The study results provide scientific data support for emergency disaster prevention work in the HJHP.

Kalman filter-based integration of GNSS and InSAR observations for local nonlinear strong deformations

The continuous monitoring of ground deformations can be provided by various methods, such as leveling, photogrammetry, laser scanning, satellite navigation systems, Synthetic Aperture Radar (SAR), and many others. However, ensuring sufficient spatiotemporal resolution of high-accuracy measurements can be challenging using only one of the mentioned methods. The main goal of this research is to develop an integration methodology, sensitive to the capabilities and limitations of Differential Interferometry SAR (DInSAR) and Global Navigation Satellite Systems (GNSS) monitoring techniques. The fusion procedure is optimized for local nonlinear strong deformations using the forward Kalman filter algorithm. Due to the impact of unexpected observations discontinuity, a backward Kalman filter was also introduced to refine estimates of the previous system’s states. The current work conducted experiments in the Upper Silesian coal mining region (southern Poland), with strong vertical deformations of up to 1 m over 2 years and relatively small and horizontally moving subsidence bowls (200 m). The overall root-mean-square (RMS) errors reached 13, 17, and 35 mm for Kalman forward and 13, 17, and 34 mm for Kalman backward in North, East, and Up directions, respectively, in combination with an external data source - GNSS campaign measurements. The Kalman filter integration outperformed standard approaches of 3-D GNSS estimation and 2-D InSAR decomposition.

Revealing the Land Subsidence Deceleration in Beijing (China) by Gaofen-3 Time Series Interferometry

Revealing the land subsidence in Beijing, China, induced by the massive groundwater extraction in the past three decades, is important to mitigate the hazards and protect the residences and infrastructure. Many SAR (Synthetic Aperture Radar) datasets have been successfully applied to reveal the land subsidence over Beijing in previous research, while few works were achieved on land subsidence revealed by time-series InSAR (Interferometric Synthetic Aperture Radar) with Gaofen-3 SAR images. In this study, we successfully perform the time-series InSAR analysis with Gaofen-3 SAR images to extract the land subsidence in Beijing from 2020 to 2021. The Sentinel-1 SAR images were used to assess the accuracy of Gaofen-3 images. The subsidence scale and extent are consistent in detected major subsidence bowls between the two datasets. The spatial–temporal evolution and the deceleration of Beijing land subsidence were revealed by comparing with the Sentinel-1 results from 2017 to 2020. Moreover, we evaluated the interferometric performance of Gaofen-3 satellite SAR imagery and analyzed the main factors that mostly influence the coherence and quality of interferograms. Our results proved that the long perpendicular baselines decrease the coherence seriously over the study area, and the artifacts induced by inaccurate orbit information reduce the quality of the Gaofen-3 interferograms. Refining and removing the two main artifacts could improve the quality of interferograms formed by Gaofen-3 SAR images.

Land subsidence simulation based on Extremely Randomized Trees combined with Monte Carlo algorithm

Land subsidence is a threat to the living safety of urban residents. Land subsidence simulation helps make prevention and conservation efforts more effective. This work introduces an innovative model for land subsidence simulation which combines Extremely Randomized Trees (ET) with the Monte Carlo algorithm to conduct quantitative analysis on the importance of different factors in land subsidence. Long-time series land surface deformation results in the eastern part of Beijing Plain were obtained by time series Interferometric Synthetic Aperture Radar (InSAR) monitoring technology based on 91 scenes of RADARSAT-2 data collected from Sept. 2009 to Jan. 2020. In addition, multi-source data, such as layered compressible layer and layered groundwater level were introduced to analysis the evolutionary mechanism of land subsidence. The results showed that the average annual rates of land subsidence at the following three typical subsidence bowls in the Beijing Plain, Chaoyang-Jinzhan, Dong Balizhuang-Da Jiaoting and Tongzhou-Song Gezhuang were −66.4 mm/a, −64.2 mm/a and −53.2 mm/a from 2010 to 2019, indicating a rather alarming subsidence issue. The most important feature affecting land subsidence is the groundwater level which contributed 67.6% to 81.8% to land surface deformation. The Monte Carlo simulation results indicated that land subsidence occurs at Dong Balizhuang-Da Jiaoting Subsidence Bowl is estimated to be the worst, with a 95% confidence interval differing by 95.7 mm/a between the upper and lower limits. Such work contributes to better understand the complex relationship between the land subsidence and influencing factors. And the simulation results provide different scientific means to facilitate efficient management and control of three land subsidence bowls.

Automatic detection and classification of land subsidence in deltaic metropolitan areas using distributed scatterer InSAR and Oriented R-CNN

Multi-temporal interferometric synthetic aperture radar (InSAR) is an effective tool for measuring large-scale land subsidence. However, the measurement points generated by InSAR are too many to be manually analyzed, and automatic subsidence detection and classification methods are still lacking. In this study, we developed an oriented R-CNN deep learning network to automatically detect and classify subsidence bowls using InSAR measurements and multi-source ancillary data. We used 541 Sentinel-1 images acquired during 2015–2021 to map land subsidence of the Guangdong-Hong Kong-Macao Greater Bay Area by resolving persistent and distributed scatterers. Multi-source data related to land subsidence, including geological and lithological, land cover, topographic, and climatic data, were incorporated into deep learning, allowing the local subsidence to be classified into seven categories. The results showed that the oriented R-CNN achieved an average precision (AP) of 0.847 for subsidence detection and a mean AP (mAP) of 0.798 for subsidence classification, which outperformed the other three state-of-the-art methods (Rotated RetinaNet, R3Det, and ReDet). An independent effect analysis showed that incorporating all datasets improved the AP by 11.2% for detection and the mAP by 73.9% for classification, respectively, compared with using InSAR measurements only. Combining InSAR measurements with globally available land cover and digital elevation model data improved the AP for subsidence detection to 0.822, suggesting that our methods can be potentially transferred to other regions, which was further validated this using a new dataset in Shanghai. These results improve the understanding of deltaic subsidence and facilitate geohazard assessment and management for sustainable environments.

Land subsidence modeling and assessment in the West Pearl River Delta from combined InSAR time series, land use and geological data

West Pearl River Delta (WPRD) is sinking as a result of the jointed effect of natural and anthropogenic factors. Land subsidence has increasingly become a concern because of tremendous population growth and rapid urbanization over this region in the last few decades. In this study, sixty-seven Sentinel-1 images, acquired between 2016 and 2021, were analyzed with the persistent scatterer interferometry technique (PSI), to monitor and reveal the ground subsidence characteristics in the WPRD. It is found that the overall vertical deformation velocities observed in the WPRD ranged between −70 mm/year and 10 mm/year. Three subsidence bowls were found in the study area (Gaolan island of Zhuhai, the junction area of Zhuhai and Zhongshan, and the junction area of Zhongshan and Jiangmen). The spatial–temporal subsidence characteristics have been analyzed. It is discovered that the ground subsidence is mostly dispersed in Quaternary deposits and is highly relevant to the thickness of sediments, indicating that soft soil consolidation is one of the primary causes contributing to land subsidence. Furthermore, land use maps for 2016 and 2021 were generated using Landsat-8 images for the investigation on the relationship between land subsidence and land use. The results obtained from analysis demonstrated that the rapid subsiding areas mainly occurred in the land-use classes as follows: aquaculture, urban land, and agricultural land. The land use conversion pattern with more significant anthropogenic influence usually causes a higher subsidence rate. In addition, based on soft soil thickness, groundwater exploitation, land use, elevation, and strata lithology, a Random Forest Regression (RFR) model was used to predict subsidence rates (R2 = 0.631, RMSE = 2.7 mm/year). The importance of these influencing factors of land subsidence was calculated based on the RFR algorithm. The results indicated that soft soil thickness, elevation, groundwater exploitation, strata lithology, and landcover type are the most significant factors affecting subsidence. The applicability of geological data and land-use history for land subsidence prediction has been demonstrated with the use of the RFR algorithm.

Correlating severe land subsidence and confined brine aquifer compaction in the Yellow River Delta, China, with Sentinel-1A/1B satellite images

In recent years, noticeable subsidence depressions have occurred along the coastal zone of the Yellow River Delta. In some coastal areas, the average annual subsidence varies from tens of millimeters to hundreds of millimeters. Although some studies have discovered a significant land subsidence funnel in the coastal zone of the Yellow River Delta, it has rarely been reported in recent years. Using Sentinel-1A/1B images from the last five years and permanent scatterers interferometric synthetic aperture radar technology, we found a typical subsidence bowl in the northeastern part of the delta, with a cumulative settlement of nearly 1 m over five years. In addition, we used the approach of soil mechanics to simulate the settlement of confined brine layers with different thicknesses under drainage conditions. We found that the 15 m thick confined aquifer can produce 1 m of settlement for every 15% decrease in water content. The simulation results explain the large settlement of the brine industrial area. This study can provide guidance for brine mining in the future. If we continue to overexploit underground brine, there will be more severe land subsidence in the delta in the future.

Urban growth and land subsidence: Multi-decadal investigation using human settlement data and satellite InSAR in Morelia, Mexico

Limited attention is typically paid to the cause-effect relationship between land subsidence due to aquifers overexploitation in expanding metropolises and urban growth models and patterns. This paper implements an integrated urban and satellite Interferometric Synthetic Aperture Radar (InSAR) approach to investigate subsidence, multi-decadal urban growth and peopling trends in the Metropolitan Area of Morelia (ZMM) in the Mexican state of Michoacán. Stacking of JRC's Global Human Settlement Layer, DLR's World Settlement Footprint and INEGI's National Geostatistical Framework datasets reveals a predominant edge-expansion growth model, with urban densification in 1975–2020 and some sprawling in 1990–2000. Population of the ZMM doubled in the last 30 years, reaching over 1 million inhabitants. The ENVISAT and Sentinel-1 InSAR analysis confirms that subsidence is structurally-controlled by the main normal faults within the Cuitzeo half-graben. Differential sinking and ground discontinuities are aligned with buried tectonic faults and contrasting compressible sediment thickness. Non-linearly deforming subsidence bowls develop at extraction wells in both old and newly urbanized sectors of the ZMM. Maximum vertical displacement velocities increased from −2.5 cm/year in 2003−2010 to −9.0 cm/year in 2014–2021, with subsidence migrating towards recently urbanized zones. More than 250 new groundwater wells were added to the public registry since 2000, many of which within new urban sectors. Time-lapse InSAR reveals a 4 km2 rapidly subsiding bowl that formed at the largest social housing neighbourhood of Villas del Pedregal, as building lots were progressively completed and sold, and new wells registered. With angular distortions due to the differential subsidence reaching 0.12% in 2014–2021, new buildings and roads are exposed to fracturing and surface faulting risk of comparable level as the city historic building blocks located along the main faults. By providing useful insights into the relationship between urban growth and land subsidence in the ZMM, the approach proves valuable for application to other metropolises worldwide.

Formation mechanism of pumping-induced earth fissures associated with a pre-existing normal fault, Beijing, China

The Shunyi (SY) earth fissures originated due to groundwater pumping and have caused severe damage to the Beijing Capital International Airport (BCIA) and other infrastructure in Beijing, China. They are distributed along the pre-existing NE-SW trending normal Shunyi-Liangxiang (SL) fault zone, with 1–30 cm vertical offsets. Previous studies of this long-term ground displacement have concluded that a subsidence bowl has formed on the footwall (northwestern side) of the SL fault in the BCIA. However, based on our observations, the relative subsidence occurred exclusively on the southeastern side, behaving like the activation of the SL fault. This complex deformation pattern has led to the formation mechanism of the SY earth fissures remaining unclear. In this study, two groups of numerical simulations based on the geological and hydrogeological settings of the SY earth fissures were performed. In Scenario I, the pumping load was on the footwall, whereas in Scenario II, the pumping load was on the hanging wall. Our results show that there are always two areas in the faulted model that have undergone pumping-induced subsidence. The first subsidence area indicates that the actual regional subsidence bowl occurred extensively around the pumping well. This subsidence was blocked by the pre-existing normal fault zone. The second subsidence area is localized and occurred abruptly in the hanging wall adjacent to the fault zone regardless of the effect of the position of the pumping load. This localized differential subsidence, representing the secondary surface faulting, was accompanied by earth fissure formation due to concentration of the surface shear strain. These simulation results agree with the actual land deformation across the SY earth fissures. Furthermore, our results suggest that both the formation location and deformation pattern of the SY earth fissures were dependent upon the characteristics of the pre-existing SL fault zone. This is mainly attributed to the fact that the fault zone has a lower resistance than the surrounding materials. Our findings reveal the formation mechanism of the SY earth fissures and provide a new understanding of pumping-induced earth fissures in other faulted areas.

Mechanism of Land Subsidence Mutation in Beijing Plain under the Background of Urban Expansion

Under the background of over-exploitation of groundwater and urban expansion, the land subsidence in the Beijing Plain has dramatically increased recently, and has demonstrated obvious mutation characteristics. Firstly, this paper used the land-use transfer matrix (LUTM) to quantify the urban expansion of Beijing, from 1990 to 2015. Secondly, the gravity center migration model (GCM) and standard deviation ellipse (SDE) methods were employed in order to quantitatively reveal the response relationship between urban expansion and land subsidence in the study area. Finally, the research innovatively combines multi-disciplinary (remote sensing, geophysical prospecting, spatial analysis, and hydrogeology), to analyze the mechanism of land subsidence mutation in the Beijing Plain, at multiple scales. The results showed the following: 1. The development direction of the urban expansion and the subsidence bowl (subsidence rate > 50 mm/year) were highly consistent, with values of 116.8° and 113.3°, respectively. 2. At the regional scale, the overall spatial distribution of subsidence mutations is controlled by the geological conditions, and the subsidence mutation time was mainly in 2005 and 2015. The area where mutation occurred in 2005 was basically located in the subsidence bowls, and the correlation between the confined water level and the subsidence rate was relatively high (r > 0.62). The area where the settlement mutation occurred in 2015, was mainly located outside the subsidence bowls, and the correlation between the confined water level and the subsidence rate was relatively low (r < 0.71). 3. In the typical subsidence area, the subsidence mutation occurred mostly in the places where the stratigraphic density is reduced, due to human activities (such as groundwater exploitation). Human activities caused the reduction in stratigraphic density, at 20 m and 90 m vertical depth in urban and rural areas, respectively. 4. At the local scale, clusters of subsidence mutation were located in the fault buffer zone, with a lateral influence range of nearly 1 km in Tongzhou. The scattered settlement mutation is distributed as a spot pattern, and the affected area is relatively small, which basically includes high-rise buildings.

Evaluating the subsidence above gas reservoirs with an elasto-viscoplastic constitutive law. Laboratory evidences and case histories

Subsidence of hydrocarbon fields is induced by the depletion of reservoirs and connected aquifers, which causes the compaction of geologic layers. A fundamental requirement of subsidence models is then capturing the mechanical behaviour of such layers with proper constitutive laws.This paper refers to the case of two gas reservoirs from the Adriatic basin, where a time delay of subsidence with respect to production was observed. This phenomenon has been reported for other fields worldwide, also occurring when subsidence rates do not immediately vanish after the end of production. In those cases, numerical simulations based on elasto-plastic constitutive models did not provide a satisfactory match of subsidence history, even though very accurate reconstructions of pore pressure evolution were used. This motivated a further insight into the time-dependent behaviour of the sandy reservoir materials. Thus, the results of an exhaustive laboratory investigation were interpreted with an elasto-viscoplastic model from the literature, implemented by the authors in a commercial Finite Element code. The satisfactory reproduction of the experimental results proved the adequacy of the model for the materials of concern and allowed determining the parameters that were then used in the new one-way coupled field simulations. A very good agreement was finally obtained between simulation results and field monitoring data in terms of reservoir compaction, maximum settlements, settlement history and extent of the subsidence bowl.

Development of spatially varying groundwater-drawdown functions for land subsidence estimation

Study regionChoshui River alluvial fan, Taiwan. Study focusLand subsidence caused by groundwater overexploitation is a critical global problem. The spatial distribution of land subsidence is crucial for effective environmental management and land planning in subsidence prone areas. Because of the nonlinear relationship between subsidence and drawdown due to groundwater exploitation in heterogeneous aquifers, a spatial regression (SR) model is developed to effectively estimate nonlinear and spatially varying land subsidence. Considering various data inputs in the Choshui River alluvial fan, the SR model offers a robust method for accurately estimating the spatial patterns of subsidence using only drawdown as input data. New hydrological insights for the regionWithout requiring extensive calibration or an elaborate numerical groundwater flow and subsidence model, the model provides annual subsidence patterns using a spatially varying relationship between drawdown and resulting land subsidence. Results show that the largest water-level cone of depression occurs in the distal fan area. Nonetheless, the calculated subsidence bowl closely approximates the observed one located much farther inland. The root-mean-square-errors (RMSEs) of annual subsidence is less or equal to 0.76 cm for the SR. Results indicate that the SR model reasonably estimates the spatial distribution of the skeletal storage coefficient in the aquifer system. The large coefficient that represents high potential of inelastic compaction occurs in the southern inland area, whereas the small coefficient that represents elastic compaction occurs in the northern area and proximal fan. Furthermore, this method can be used efficiently for subsidence management/ regulation and might be widely used for subsidence estimation solely based on drawdown.

Adjacent-Track InSAR Processing for Large-Scale Land Subsidence Monitoring in the Hebei Plain

Large-scale land subsidence has threatened the safety of the Hebei Plain in China. For tens of thousands of square kilometers of the Hebei Plain, large-scale subsidence monitoring is still one of the most difficult problems to be solved. In this paper, we employed the small baseline subset (SBAS) and NSBAS technique to monitor the land subsidence in the Hebei Plain (45,000 km2). The 166 Sentinel-1A data of adjacent-track 40 and 142 collected from May 2017 to May 2019 were used to generate the average deformation velocity and deformation time-series. A novel data fusion flow for the generation of land subsidence velocity of adjacent-track is presented and tested, named as the fusion of time-series interferometric synthetic aperture radar (TS-InSAR) results of adjacent-track using synthetic aperture radar amplitude images (FTASA). A cross-comparison analysis between the two tracks results and two TS-InSAR results was carried out. In addition, the deformation results were validated by leveling measurements and benchmarks on bedrock results, reaching a precision 9 mm/year. Twenty-six typical subsidence bowls were identified in Handan, Xingtai, Shijiazhuang, Hengshui, Cangzhou, and Baoding. An average annual subsidence velocity over −79 mm/year was observed in Gaoyang County of Baoding City. Through the cause analysis of the typical subsidence bowls, the results showed that the shallow and deep groundwater funnels, three different land use types over the building construction, industrial area, and dense residential area, and faults had high spatial correlation related to land subsidence bowls.

Monitoring of land surface subsidence using persistent scatterer interferometry techniques and ground truth data in arid and semi-arid regions, the case of Remah, UAE

The United Arab Emirates (UAE) is located in an arid desert climate with very limited water resources and scarce rainfall. Along with the fast development of the country, the water demand for agriculture, industrial, and domestic purposes increased and led to diminishing groundwater resources. In this study, we explore the land surface deformations due to groundwater overexploitation in the agricultural area of Remah by analyzing Sentinel-1 data between 2015 and 2019 with the novel Parallelized-Persistent Scatterer Interferometry (P-PSI) technique. The detected land surface deformations have been correlated to the recorded groundwater levels at nearby water wells. This study detected land surface deformations in a form of an extensive subsidence bowl (with 28.5 km in diameter) with a maximum subsidence rate of 40 mm/year and a standard deviation within the bowl of less than 2 mm/year. The detected subsidence was associated with a 12 m drop in the water table level within the study area. The Persistent Scatterers with the highest deformations rate were spatially correlated with the depression cone of the groundwater level. These findings provide useful insights in understanding the groundwater regime of the area and have an important role in assessing regional hazards and driving mitigation measures towards managing uncontrolled groundwater overexploitation for sustainable management of groundwater resources.

Assessments of land subsidence along the Rizhao–Lankao high-speed railway at Heze, China, between 2015 and 2019 with Sentinel-1 data

Abstract. The Heze section of Rizhao–Lankao high-speed railway (RLHR-HZ) has been under construction since 2018 and will be in operation by the end of 2021. However, there is a concern that land subsidence in the Heze region may affect the regular operation of RLHR-HZ. In this study, we investigate the contemporary ground deformation in the region between 2015 and 2019 by using more than 350 C-band interferograms constructed from two tracks of Sentinel-1 data over the region. The small baseline subset (SBAS) technique is adopted to compile the time-series displacement. We find that the RLHR-HZ runs through two main subsidence areas: one is located east of the Heze region with rates ranging from −4 to −1 cm yr−1, and another one is located in the coalfield with rates ranging from −8 to −2 cm yr−1. A total length of 35 km of RLHR-HZ is affected by the two subsidence basins. Considering the previous investigation and the monthly precipitation, we infer that the subsidence bowl east of the Heze region is due to massive extraction of deep groundwater. Close inspections of the relative locations between the second subsidence area and the underground mining reveals that the subsidence there is probably caused by the groundwater outflow and fault instability due to mining, rather than being directly caused by mining. The InSAR-derived ground subsidence implies that it is necessary to continue monitoring the ground deformation along RLHR-HZ.

Multi-Temporal Satellite Interferometry for Fast-Motion Detection: An Application to Salt Solution Mining

Underground mining is one of the human activities with the highest impact in terms of induced ground motion. The excavation of the mining levels creates pillars, rooms and cavities that can evolve in chimney collapses and sinkholes. This is a major threat where the mining activity is carried out in an urban context. Thus, there is a clear need for tools and instruments able to precisely quantify mining-induced deformation. Topographic measurements certainly offer very high spatial accuracy and temporal repeatability, but they lack in spatial distribution of measurement points. In the past decades, Multi-Temporal Satellite Interferometry (MTInSAR) has become one of the most reliable techniques for monitoring ground motion, including mining-induced deformation. Although with well-known limitations when high deformation rates and frequently changing land surfaces are involved, MTInSAR has been exploited to evaluate the surface motion in several mining area worldwide. In this paper, a detailed scale MTInSAR approach was designed to characterize ground deformation in the salt solution mining area of Saline di Volterra (Tuscany Region, central Italy). This mining activity has a relevant environmental impact, depleting the water resource and inducing ground motion; sinkholes are a common consequence. The MTInSAR processing approach is based on the direct integration of interferograms derived from Sentinel-1 images and on the phase splitting between low (LF) and high (HF) frequency components. Phase unwrapping is performed for the LF and HF components on a set of points selected through a “triplets closure” method. The final deformation map is derived by combining again the components to avoid error accumulation and by applying a classical atmospheric phase filtering to remove the remaining low frequency signal. The results obtained reveal the presence of several subsidence bowls, sometimes corresponding to sinkholes formed in the recent past. Very high deformation rates, up to −250 mm/yr, and time series with clear trend changes are registered. In addition, the spatial and temporal distribution of velocities and time series is analyzed, with a focus on the correlation with sinkhole occurrence.