Land subsidence is a global environmental geological hazard caused by natural or human activities. The high spatial resolution and continuous time coverage of interferometric synthetic aperture radar (InSAR) time series analysis techniques provide data and a basis for the development of methods for the investigation and evolution mechanism study of regional land subsidence. Beijing, the capital city of China, has suffered from land subsidence for decades since it was first recorded in the 1950s. It was reported that uneven ground subsidence and fractures have seriously affected the normal operation of the Beijing Capital International Airport (BCIA) in recent years before the overhaul of the middle runway in April 2017. In this study, InSAR time series analysis was successfully used to detect the uneven local subsidence and ground fissure activity that has been gradually increasing in BCIA since 2010. A multi-temporal InSAR (MT-InSAR) technique was performed on 63 TerraSAR-X/TanDem-X (TSX/TDX) images acquired between 2010 and 2017, then deformation rate maps and time series for the airport area were generated. Comparisons of deformation rate and displacement time series from InSAR and ground-leveling were carried out in order to evaluate the accuracy of the InSAR-derived measurements. After an integrated analysis of the distribution characteristics of land subsidence, previous research results, and geological data was carried out, we found and located an active ground fissure. Then main cause of the ground fissures was preliminarily discussed. Finally, it can be conducted that InSAR technology can be used to identify and monitor geological processes, such as land subsidence and ground fissure activities, and can provide a scientific approach to better explore and study the cause and formation mechanism of regional subsidence and ground fissures.

The long-term overexploitation of groundwater leads to serious land subsidence and threatens the safety of Beijing-Tianjin-Hebei (BTH). In this paper, an interferometric point target analysis (IPTA) with small baseline subset InSAR (SBAS-InSAR) technique was used to derive the land subsidence in a typical BTH area from 2012 to 2018 with 126 Radarsat-2 and 184 Sentinel-1 images. The analysis reveals that the average subsidence rate reached 118 mm/year from 2012 to 2018. Eleven subsidence features were identified: Shangzhuang, Beijing Airport, Jinzhan and Heizhuanghu in Beijing, Guangyang and Shengfang in Langfang, Wangqingtuo in Tianjin, Dongguang in Cangzhou, Jingxian and Zaoqiang in Hengshui and Julu in Xingtai. Comparing the different types of land use in subsidence feature areas, the results show that when the land-use type is relatively more complex and superimposed with residential, industrial and agricultural land, the land subsidence is relatively more significant. Moreover, the land subsidence development patterns are different in the BTH areas because of the different methods adopted for their water resource development and utilization, with an imbalance in their economic development levels. Finally, we found that the subsidence changes are consistent with groundwater level changes and there is a lag period between land subsidence and groundwater level changes of approximately two months in Beijing Airport, Jinzhan, Jingxian and Zaoqiang, of three months in Shangzhuang, Heizhuanghu, Guangyang, Wangqingtuo and Dongguang and of four months in Shengfang.

In the Beijing Plain, land subsidence is one of the most prominent geological problems, which is affected by multiple factors. Groundwater exploitation, thickness of the Quaternary deposit and urban development and construction are important factors affecting the formation and development of land subsidence. Here we choose groundwater level change, thickness of the Quaternary deposit and index-based built-up index (IBI) as influencing factors, and we use the influence factors to predict the subsidence amount in the Beijing Plain. The Sentinel-1 radar images and the persistent scatters interferometry (PSI) were adopted to obtain the information of land subsidence. By using Google Earth Engine platform and Landsat8 optical images, IBI was extracted. Groundwater level change and thickness of the Quaternary deposit were obtained from hydrogeological data. Machine learning algorithms Linear Regression and Principal Component Analysis (PCA) were used to investigate the relationship between land subsidence and influencing factors. Based on the results obtained by Linear Regression and PCA, a suitable machine learning algorithm was selected to predict the subsidence amount in the Beijing Plain in 2018 through influencing factors. In this study, we found that the maximum subsidence rate in the Beijing Plain had reached 115.96 mm/y from 2016 to 2018. The land subsidence was serious in eastern Chaoyang and northwestern Tongzhou. In addition, the area where thickness of the Quaternary deposit reached 150–200 m was prone to more serious land subsidence in the Beijing Plain. In groundwater exploitation, the second confined aquifer had the greatest impact on land subsidence. Through Linear Regression and PCA, we found that the relationship between land subsidence and influencing factors was nonlinear. XGBoost was feasible to predict subsidence amount. The prediction accuracy of XGBoost on the subsidence amount reached 0.9431, and the mean square error was controlled at 15.97. By using XGBoost to predict the subsidence amount, our research provides a new idea for land subsidence prediction.

Geological disasters, including ground deformation, fractures and collapse, are serious problems in coal mining regions, which have threatened the sustainable development for local industry. The Ordos Basin is most known for its abundant coal resources. Over-mining the underground coal resources had induced land deformation. Detecting the evolution of the land deformation features and identifying the potential risk are important for decision-makers to prevent geological disasters. We analyzed land subsidence induced by coal mining in a 200 ${\mathrm{km}}^2$ area in the Ordos Basin for the time period 2006–2015. ALOS-1 PALSAR images from December 2006 to January 2011 and ALOS-2 PALSAR-2 images from December 2014 to July 2015, optical remotely sensed images and coal mining information were collected. The small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) method and differential interferometric synthetic aperture radar (D-InSAR) method, GIS and statistical analysis were adopted. Results show that the maximum subsidence rate and cumulative subsidence along the line of sight (LOS) were −65 mm/year and −246 mm, respectively, from December 2006 to January 2011. The maximum cumulative subsidence was −226 mm from December 2014 to July 2015. The new boundary of the mining goafs from 2014 to 2015 and the most dangerous risk region were mapped. Moreover, the effect of large-scale mining coal, with the production volume exceeds 1.2 million tons per year, with the operation time more than 20 years on land subsidence was found greater than small and medium-scale coal mines and reached −59 mm/year. The recently established small-sized and medium-sized coal mines show high land subsidence. This study will contribute to better understand the land subsidence process in mining region and provide scientific support for government to prevent land subsidence.

Groundwater resources have been exploited and utilized on a large scale in the North China Plain (NCP) since the 1970s. As a result of extensive groundwater depletion, the NCP has experienced significant land subsidence, which threatens geological stability and infrastructure health and exacerbates the risks of other geohazards. In this study, we employed multi-track Synthetic Aperture Radar (SAR) datasets acquired by the Sentinel-1A (S1A) satellite to detect spatial and temporal distributions of surface deformation in the NCP from 2016 to 2018 based on multi-temporal interferometric synthetic aperture radar (MT-InSAR). The results show that the overall ground displacement ranged from −165.4 mm/yr (subsidence) to 9.9 mm/yr (uplift) with a standard variance of 28.8 mm/yr. During the InSAR monitoring period, the temporal pattern of land subsidence was dominated by a decreasing tendency and the spatial pattern of land subsidence in the coastal plain exhibited an expansion trend. Validation results show that the S1A datasets agree well with levelling data, indicating the reliability of the InSAR results. With groundwater level data, we found that the distribution of subsidence in the NCP is spatially consistent with that of deep groundwater depression cones. A comparison with land use data shows that the agricultural usage of groundwater is the dominant mechanism responsible for land subsidence in the whole study area. Through an integrated analysis of land subsidence distribution characteristics, geological data, and previous research results, we found that other triggering factors, such as active faults, precipitation recharge, urbanization, and oil/gas extraction, have also impacted land subsidence in the NCP to different degrees.

The long-term over-exploitation of groundwater has not only caused the compaction of aquifer thickness and surface deformation but has also further aggravated the loss of groundwater storage (GWS) in Beijing plain. The South-to-North Water Diversion Project (SNWDP) furnishes a new source of water for Beijing. By reviewing related studies, it was found that there are few studies on the realization of GWS estimation based on InSAR technology considering the lag effect. Therefore, in this study, firstly, the long-time series deformation characteristics of Beijing plain were obtained from 46 RadarSat-2 images using small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technology. Secondly, the seasonal components of surface deformation and hydraulic head change were extracted by means of multichannel singular spectrum analysis (MSSA), verifying the separation accuracy by means of Monto Carlo-SSA (MC-SSA). Finally, for the hydrodynamic delay (aquifer water supply/drainage) of the complex aquifer system, we introduced the time lag cross-correlation (TLCC) approach to correct the hysteresis response of seasonal deformation relative to the variation of the aquifer system head, so as to realize the estimation of aquifer storage properties and GWS loss, even unrecoverable GWS (UGWS). The results showed that the average annual variation of total GWS (TGWS) in Beijing plain was −6.702 × 10⁷ m³, of which the depletion volume of UGWS was −6.168 × 10⁷ m³, accounting for 92.03% of the TGWS. On a temporal scale, the depletion of UGWS lagged behind the total head change, with about one year of lag time. On a spatial scale, in contrast to the north of Beijing plain, the depletion of UGWS in the south only recovered briefly after 2015 and then continued to decline. This further indicated that the process of the decline of middle-deep confined head and long-term GWS loss caused by over-exploitation of groundwater was irreversible. These findings are of great significance to optimize the allocation of groundwater resources, reduce the harm of land subsidence and protect groundwater resources.