Over-exploitation of groundwater has caused severe land subsidence in Beijing during the past two decades. Since the middle route of South-to-North Water Diversion Project (SNWDP), the biggest water diversion project in China, started to deliver water to Beijing in December 2014, the groundwater shortage has been greatly alleviated. This study aims to analyze the impact of SNWDP on the spatiotemporal evolution of land subsidence in Beijing. Change in surface displacement in Beijing after SNWDP was retrieved and the spatiotemporal patterns of the change were analyzed based on long time-series Envisat Advanced Synthetic Aperture Radar (ASAR) (2004–2010), Radarsat-2 (2011–2014), and Sentinel-1 (2015–2017) satellite datasets using Permanent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) techniques. Land subsidence unevenness index (LSUI) was proposed to represent the spatial unevenness of surface displacement. PS-Time approach was then adapted to examine the time series evolution of LSUI. The results showed that the InSAR measurements agree well with leveling measurements with R² over 0.96. Although the maximum annual displacement rate reached −159.7 mm/year by 2017, over 57% of the area within 25 mm/year contour line showed decreasing or unchanged displacement rate after the south-north water delivered to Beijing. The settlement rate in Chaoyang-Dongbalizhuang (CD) subsidence center has decreased for 26 mm/year from 2011–2014 to 2015–2017. Only around 15% of the area experienced continued accelerating settlement rate through the three time periods, which was mainly located in the area with the compressible layer thickness over 190 m, while the magnitude of velocity increment considerably decreased after SNWDP. Land subsidence unevenness, represented by LSUI, developed more slowly after SNWDP than that during 2011–2014. However, LSUI at the edge of settlement funnel has kept developing and reached 1.7‰ in 2017. Decreasing groundwater level decline after SNWDP and the positive relationship (R² > 0.74) between land subsidence and groundwater level clearly showed impacts of SNWDP on the alleviating land subsidence. Other reasons include geological background, increasing precipitation, and strict water management policies implemented during these years.

In urban areas, deformation of transportation infrastructures may lead to serious safety accidents. Timely and accurate monitoring of the structural deformation is critical for prevention of transportation accidents and assurance of construction quality, particularly in areas with regional land subsidence, such as the city of Beijing. In this study, we proposed a method for the detection of seasonal deformation of highway overpasses using the integration of persistent scatterers Interferometric Synthetic Aperture Radar (PS-InSAR) techniques and seasonal indices, i.e., deformation concentration degree (DCD) and deformation concentration period (DCP) indices. Taking eastern Beijing urban area as a case study area, we first used the PS-InSAR technique to derive time series surface deformation based on 55 TerraSAR-X images during 2010–2016. Then, we proposed DCD and DCP indices to characterize seasonal deformation of 25 highway overpasses in the study area, with DCD representing to what degree the annual deformation is distributed in a year, and DCP representing the period on which deformation concentrates in the year. Our results showed that the maximum annual deformation rate reached −141.3 mm/year in Beijing urban area, and the PS-InSAR measurements agreed well with levelling measurements (R² > 0.97). For PS pixels with DCD ≥ 0.3, the monthly deformation showed obvious seasonal patterns with deformation values during some months greater than those during the other months. DCP revealed that the settlement during autumn and winter was more serious than that in spring and summer. The seasonal patterns seemed to be related to the location, structure, and construction age of the overpasses. The upper-level overpasses, the newly constructed overpasses, and those located in the subsidence area (rate < −40 mm/year) tended to show a greater seasonal pattern. The seasonal deformation variations were also affected by groundwater-level fluctuation, temperature, and compressible layer.

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.

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.

The construction of subways exacerbates the non-uniformity of surface deformation, which in turn poses a potential threat to the safe construction and stable operation of urban rail transit systems. Beijing, the city with the most extensive subway network in China, has long been affected by land subsidence. Utilizing data from Envisat ASAR, Radarsat-2, and Sentinel-1 satellites, this study employs PS-InSAR technology to monitor and analyze land subsidence within a 2 km buffer zone along Beijing Metro Line 22 over a span of 20 years (from January 2004 to November 2024). The results indicate that land subsidence at Guanzhuang Station and Yanjiao Station along Metro Line 22 is particularly pronounced, forming two distinct subsidence zones. After 2016, the overall rate of subsidence along the subway line began to stabilize, with noticeable ground rebound emerging around 2020. This study further reveals a strong correlation between land subsidence and confined groundwater levels, while geological structures and building construction also exert a significant influence on subsidence development. These findings provide a crucial scientific foundation for the formulation of effective prevention and mitigation strategies for land subsidence along urban rail transit lines.