Search

Rapid simulation of land subsidence can provide an effective means of facilitating its management and control. This paper innovatively introduced a back-propagation (BP) neural network subsidence simulation model which was optimized by the sparrow search algorithm (SSA), hereinafter referred to as the SSA-BP model, to simulate land subsidence in Tongzhou District, Beijing. First, based on the time series interferometric synthetic aperture radar (InSAR) monitoring, different technologies such as spatial analysis, Google Earth Engine (GEE), and machine learning were utilized together with a variety of multi-source spatial data, including groundwater level, compressible layer thickness, Visible Infrared Imager Radiometer Suite (VIIRS) nighttime light images, and the OpenStreetMap (OSM) road distribution. Furthermore, we determined the optimal grid scale for land subsidence research and worked out a multifactor-driven SSA-BP land subsidence simulation model for which sensitivity analysis was performed with great care. Main conclusions: (1) From November 2010 to January 2020, the average annual surface displacement rate in Beijing's subcentre, Tongzhou, ranged from −133.9 to +3.9 mm/year. (2) Our SSA-BP land subsidence simulation model allows for a relatively high overall accuracy. The average root mean square error (RMSE) is 5.00 mm/year, the average mean absolute error (MAE) is 3.80 mm/year, and the average coefficient of determination (R2) is 0.83. (3) Sensitivity analysis shows that the SSA-BP model driven by multiple factors has a far better simulation performance than the model driven by any single weighting factor, and displays the highest sensitivity to the groundwater level factor among all the weighting factors. In terms of subdividing the study area, our SSA-BP land subsidence model runs with multifunctional zones displayed a higher accuracy than that without. This paper made use of a machine learning model to simulate land subsidence in Beijing's Tongzhou District and explored the sensitivity of different factors to land subsidence, which is helpful for its scientific prevention. [ABSTRACT FROM AUTHOR]

Xiong'an New Area is a national event and a project planned for a millennium of China. Its high-quality construction is of great significance to easing the noncapital functions of Beijing and the coordinated development of the Beijing-Tianjin-Hebei region. As an emerging city, the development and construction of Xiong'an New Area is bound to be restricted by geological and resource conditions. Therefore, geo-environment suitability analysis is the necessary basis of urban development and construction. Geo-environment suitability analysis of urban construction is a complex process that requires various geological indicator information, and relevant expertise to analyze their relevance. This paper focuses on the analytic hierarchy process (AHP) for the assessment of geo-environment suitability for urban construction in Rongcheng district, which is a Start Construction Region in Xiong'an New Area. Multiple factors, including the characteristic value of bearing capacity of foundation soil, land subsidence rate, geological faults, ground fissures, potential liquefied sands, quality of groundwater chemistry, quality of soil chemistry, chemical corrosion of concrete by groundwater, chemical corrosion of steel by groundwater, and enrichment of deep groundwater and geothermal resource, were used for the suitability assessments. From the evaluation achievements, the high and very high suitable lands for urban construction, with an acreage percentage of 89.2%, were located in most parts of the study area. Meanwhile, for another 9.1% of the land, the impacts of geological faults, land subsidence, and potential liquefied sands needed to be noted preferentially for urban construction. [ABSTRACT FROM AUTHOR]

Advanced techniques of multi-temporal InSAR (MT-InSAR) represent a valuable tool in ground subsidence studies allowing remote investigation of the behavior of mass movements in long time intervals by using large datasets of SAR images covering the same area and acquired at different epochs. The settlement in Beijing is very serious, producing adverse environmental impacts. Excessive groundwater withdrawal is the most direct cause of land subsidence.In this paper, a spatial-temporal analysis of the land subsidence in Beijing was performed using one of the MT-InSAR techniques, referred to as Small Baseline Subset (SBAS). This technique allows monitoring the temporal evolution of a deformation phenomenon, via the generation of mean deformation velocity maps and displacement time series from a data set of acquired SAR images. 52 C-band ENVISAT ASAR images acquired from June 2003 to August 2010 were used to produce a linear deformation rate map and to derive time series of ground deformation. The results show that there are three large subsidence funnels within this study area, which separately located in Balizhuang-Dajiaoting in Chaoyang district, Wangjing-Laiguangying Chaoyang district, Gaoliying Shunyi district. The maximum settlement center is Wangsiying-Tongzhou along the Beijing express; the subsidence velocity exceeds 110 mm/y in the LOS direction. In particular, we compared the achieved results with leveling measurements that are assumed as reference. The estimated long-term subsidence results obtained by SBAS approach agree well with the development of the over-exploitation of ground water, indicating that SBAS techniques is adequate for the retrieval of land subsidence in Beijing from multi-temporal SAR data. [ABSTRACT FROM AUTHOR]

Beijing, China's capital city, has experienced decades of severe land subsidence due to the long-term overexploitation of groundwater. The implementation of the South-to-North Water Diversion Project (SNWDP) and artificial ecological restoration have significantly changed Beijing's hydro-ecological and geological environment in recent years, leading to a widespread rise in groundwater levels. However, whether the related land subsidence has slowed down or reversed under these measures has not yet been effectively monitored and quantitatively analyzed in terms of time and space. Accordingly, in this study, we proposed using an improved time-series deformation method, which combines persistent scatterers and distributed scatterers, to process Sentinel-1 images from 2015 to 2022 in the Beijing Plain region. We performed a geospatial analysis to gain a better understanding of how the new hydrological conditions changed the pattern of deformation on the Beijing Plain. The results indicated that our combined PS and DS method provided more measurements both in total quantity and spatial density than conventional PSI methods. The land subsidence in the Beijing Plain area has been effectively alleviated from a subsidence region of approximately 1377 km2 in 2015 to only approximately 78 km2 in 2022. Ecological restoration areas in the northeastern part of the Plain have even rebounded over this period, at a maximum of approximately 40 mm in 2022. The overall pattern of ground deformation (subsidence and uplift) is negatively correlated with changes in the groundwater table (decline and rise). Local deformation is controlled by the thickness of the compressible layer and an active fault. The year 2015, when anthropogenic water transfers were eliminated and ecological measures to recharge groundwater were implemented, was the crucial turning point of the change in the deformation trend in the subsidence history of Beijing. Our findings carry significance, not only for China, but also for other areas where large-scale groundwater extractions are causing severe ground subsidence or rebound. [ABSTRACT FROM AUTHOR]

The coordinated development of the Beijing-Tianjin-Hebei has become a national strategy with Beijing and Tianjin as twin engines driving the regional development. However, the Beijing-Tianjin-Hebei region has suffered dramatic ground subsidence during last two to three decades, mainly due to long-term groundwater withdrawal. Although, annual spirit leveling has been conducted routinely in some parts of Beijing and Tianjin, and InSAR technique has also been used to monitor ground subsidence in some local areas of the region, there is a lack of a complete survey of ground subsidence over the whole region. In this paper, we report a research on mapping ground subsidence in the Beijing-Tianjin-Hebei region over a long time span from 1992 to 2014. Three SAR datasets from four satellites are used: ERS-1/2 SAR images from 1992 to 2000, ENVISAT ASAR images from 2003 to 2010, and RADARSAT-2 images from 2012 to 2014. An improved multi-temporal InSAR method, namely "Multiple-master Coherent Target Small-Baseline InSAR" (MCTSB-InSAR), has been developed to process the datasets. A unique feature of MCTSB-InSAR is the adjustment process useful for wide area monitoring which provides an integrated solution for both calibration of InSAR-derived deformation and the harmonization of the deformation estimates from overlapping SAR frames. Three maps of the subsidence rate corresponding to the three periods over the wide Beijing-Tianjin-Hebei region are generated, with respective accuracy of 8.7 mm/year (1992-2000), 4.7 mm/year (2003-2010), and 5.4 mm/year (2012-2014) validated by more than 120 leveling measurements. The spatial-temporal characteristics of the development of ground subsidence in Beijing and Tianjin are analyzed. This research represents a first-ever effort on mapping ground subsidence over very large area and over long time span in China. The result is of significance to serve the decision-making on ground subsidence mitigation in the Beijing-Tianjin-Hebei region. [ABSTRACT FROM AUTHOR]

Excessive exploitation of the groundwater has resulted in obvious three-dimensional (3D) deformation features on the surface of the Beijing Plain. This paper, by combining Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) technologies, has obtained time-series information of the 3D surface deformation in the Beijing Plain, analyzing its spatial distribution characteristics. On this basis, the relationship between different controlling factors with the 3D deformation of the surface has been analyzed as well. The following results are obtained: (1) From 2013 to 2018, the land subsidence, which generally showed the trend of slowing down, was mainly concentrated in the eastern, northern, and southern regions of Beijing Plain, with multiple subsidence centers. (2) Under the International Terrestrial Reference Frame 2005 (ITRF2005), the horizontal direction of all GPS points in the plain is basically the same, with the dominant movement direction being NE112.5°~NE113.8°. Under the Eurasian reference frame, the horizontal movement rate of GPS points significantly decreases. The movement rate and direction of each point are not characteristic of overall trend activity. (3) The distribution and extent of the 3D surface deformation in the Beijing Plain are controlled by the basement structure. Part of the subsided area corresponds to a Quaternary depression formed at the junction of active faults disrupting the area. Similarly, the distribution of horizontal deformation in the E-W and N-S directions of the plain is controlled by the regional basement structure comprising major faults bounding horizontal deformation. (4) Groundwater exploitation is the main cause of the 3D surface deformation in the Beijing Plain. The groundwater funnels of the second and third confined aquifer are in suitable agreement with the land subsidence. The horizontal movement in the Beijing Plain is either directed toward the center of the groundwater or the land subsidence funnel, and the deformation is directed from areas with higher to areas with lower groundwater levels. [ABSTRACT FROM AUTHOR]

Beijing-Tianjin-Hebei (BTH) has been suffering from severe groundwater storage (GWS) consumption and land subsidence (LS) for a long period. The overexploitation of groundwater brings about severe land subsidence, which affects the safety and development of BTH. In this paper, we utilized multi-frame synthetic aperture radar datasets obtained by the Rardarsat-2 satellite to monitor land subsidence's temporal and spatial distribution in the BTH from 2012 to 2016 based on multi-temporal interferometric synthetic aperture radar (MT-InSAR). In addition, we also employed the Gravity Recovery and Climate Experiment (GRACE) mascon datasets acquired by the Center for Space Research (CSR) and Jet Propulsion Laboratory (JPL) to obtain the GWS anomalies (GWSA) of BTH from 2003 to 2016. Then we evaluate the accuracy of the results obtained. Furthermore, we explored the relationship between the regional GWSA and the average cumulative subsidence in the BTH. The total volume change of subsidence is 59.46% of the total volume change of groundwater storage. Moreover, the long-term decreasing trend of the GWSA (14.221 mm/year) and average cumulative subsidence (17.382 mm/year) show a relatively high consistency. Finally, we analyze the heterogeneity of GWS change (GWSC) and LS change (LSC) in the four typical areas by the Lorenz curve model. The implementation of the South-to-North Water Diversion Project (MSWDP) affects the heterogeneity of GWSC and LSC. It can be seen that the largest heterogeneity of LSC lags behind the GWSC in the Tianjin-Langfang-Hengshui-Baoding area. The largest uneven subsidence in Beijing and Tianjin occurred in 2015, and the largest uneven subsidence in Hengshui-Baoding occurred in 2014. After that, the heterogeneity of subsidence gradually tends to stable. [ABSTRACT FROM AUTHOR]

Surface subsidence is a serious threat to human life, buildings and traffic in Beijing. Surface subsidence is closely related to human activities, and human activities in Beijing area showed a decreasing trend during the Corona Virus Disease 2019 (COVID-19). To study surface subsidence in Beijing before and after the COVID-19 outbreak and its causes, a total of 51 Sentinel-1A SAR images covering Beijing from January 2018 to April 2022 were selected to derive subsidence information by Time Series Interferometry Synthetic Aperture Radar (TS-InSAR). The results of surface subsidence in Beijing demonstrate that Changping, Chaoyang, Tongzhou and Daxing Districts exhibited the most serious subsidence phenomenon before the COVID-19 outbreak. The four main subsidence areas form an anti-Beijing Bay that surrounds other important urban areas. The maximum subsidence rate reached −57.0 mm/year. After the COVID-19 outbreak, the main subsidence area was separated into three giant subsidence funnels and several small subsidence funnels. During this period, the maximum subsidence rate was reduced to −43.0 mm/year. Human activity decrease with the COVID-19 outbreak. This study effectively analysed the influence of natural factors on surface subsidence after excluding most of the human factors. The following conclusions are obtained from the analysis: (1) Groundwater level changes, Beijing's geological structure and infrastructure construction are the main reasons for surface subsidence in Beijing. (2) Seasonal changes in rainfall and temperature indirectly affect groundwater level changes, thereby affecting surface subsidence in the area. (3) The COVID-19 outbreak in early 2020 reduced the payload of Beijing's transportation facilities. It also slowed down the progress of various infrastructure construction projects in Beijing. These scenarios affected the pressure on the soft land base in Beijing and reduced the surface subsidence trend to some extent. [ABSTRACT FROM AUTHOR]

Land subsidence threatens the stable operation of urban rail transit, including subways. Obtaining deformation information during the entire life-cycle of a subway becomes a necessary means to guarantee urban safety. Restricted by sensor life and cost, the single-sensor Multi-temporal Interferometric Synthetic Aperture Radar (MTI) technology has been unable to meet the needs of long-term sequence, high-resolution deformation monitoring, especially of linear objects. The multi-sensor MTI time-series fusion (MMTI-TSF) techniques has been proposed to solve this problem, but rarely mentioned. In this paper, an improved MMTI-TSF is systematically explained and its limitations are discussed. Taking the Beijing Subway Network (BSN) as a case study, through cross-validation and timing verification, we find that the improved MMTI-TSF results have higher accuracy (R2 of 98% and, Root Mean Squared Error (RMSE) of 4mm), and compared with 38 leveling points, the fitting effect of the time series is good. We analyzed the characteristics of deformation along the BSN over a 15-year periods. The results suggest that there is a higher risk of instability in the eastern section of Beijing Subway Line 6 (L6). The land subsidence characteristics along the subway lines are related to its position from the subsidence center, and the edge of the subsidence center presents a segmented feature. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The eastern Beijing plain has been suffering severe subsidence for the last decades, mainly associated with the long-term excessive extraction of groundwater resource. Since the end of 2014, the annual water supply in Beijing plain has reached several hundred million cubic meters because of the South-to-North Water Diversion (SNWD) Project, which has reduced the groundwater exploitation and changed the status of land subsidence. In this work, we first obtain the current spatiotemporal variations of land subsidence in the eastern Beijing plain by using progressive small baseline subsets (SBAS) InSAR time series analysis method with Sentinel-1 SAR data acquired from July 2015 to December 2021. Then, we analyze the correlations between InSAR-derived subsidence and groundwater level change by applying the cross wavelet method. The results show that two major subsidence zones are successfully detected with the maximum deformation rate of −150 mm/yr and maximum cumulative deformation of −950 mm. Besides, the ground deformation at different stages from 2016 to 2021 reveal that the area and magnitude of major deformation significantly slow down, even in the regions with severe subsidence, especially in the year of 2017, which is about two years later than the start time of SNWD Project in Beijing. Further, we find the InSAR-derived subsidence lags groundwater level change with about 1–2-month lagging time, indicating that the dynamic variation of groundwater level fluctuation may be the main factor affecting the uneven subsidence in the severe subsiding zones. Last, differential subsidence rates are identified at both sides of geological faults, such as Nankou-Sunhe fault and Nanyuan-Tongxian fault, from the observed deformation map, which could be explained that the groundwater flow is blocked when a fault is encountered. These findings can provide significant information to reveal the deformation mechanisms of land subsidence, establish the hydrogeological models and assist decision-making, early warning and hazard relief in Beijing, China. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 × 107 m3, of which the depletion volume of UGWS was −6.168 × 107 m3, 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. [ABSTRACT FROM AUTHOR]

Flooding is one of the most destructive natural events that severely damage the ground and inundate underground infrastructure. Subway systems in metropolitan areas are susceptible to flooding, which may be exacerbated when land subsidence occurs. However, previous studies have focused on flood risk evaluation on regional/watershed-scales and land subsidence monitoring in plains, instead of on subway flood risk evaluation and how land subsidence aggravates the flood risk in subway systems. Using the proposed risk indicators and field survey data, we present a method assessing the flood risk of metropolitan subway systems under a subsidence condition based on the fuzzy analytic hierarchy process (FAHP) combined with a geographic information system (GIS). We use the regional risk level within the 500 m buffer zone of the subway line to depict the flood risk of the subway system. The proposed method was used to evaluate the flood risk of the Beijing subway system. The results show that the flood risks of the Beijing subway show a ring-like distribution pattern—risk levels decreasing from the central urban area to the suburbs. Very high and high risks are mainly located within third and fourth ring roads, accounting for 63.58% (29.40 km2) and 63.83% (81.19 km2) of the total area. Land subsidence exacerbated the Beijing subway system's flood risk level—the moderate to very high risk increased by 46.88 km2 (16.33%), indicating that land subsidence is an essential factor affecting the flood risk level of subway systems. In addition to enhancing flood warnings, future subway flooding could be reduced by elevating the height of the stations' exit (entrance) and installing water stop plates and watertight doors. This study is of great significance for flood warning and prevention in the Beijing subway system; it provides a theoretical basis for flood risk evaluation in other metropolitan areas. [ABSTRACT FROM AUTHOR]

Beijing is severely affected by land subsidence, and rapid urbanisation and building construction might accelerate the land subsidence process. Based on 39 Envisat Advanced Synthetic Aperture Radar (ASAR) images acquired between 2003–2010, 55 TerraSAR-X images acquired between 2010–2016, and urban building information, we analysed the relationship between land subsidence and buildings at the regional, block, and building scales. The results show that the surface displacement rate in the Beijing urban area ranged from −109 mm/year to +13 mm/year between 2003–2010, and from −151 mm/year to +19 mm/year between 2010–2016; two subsidence bowls were mainly distributed in the eastern part of the Chaoyang District. The displacement rate agreed well with the levelling measurements, with an average bias of less than six mm/year. At the regional scale, the spatial pattern of land subsidence was mainly controlled by groundwater extraction, compressible layer thickness, and geological faults. Subsidence centres were located in the area around ground water funnels with a compressible layer depth of 50–70 m. The block-scale analysis demonstrated a clear correlation between the block construction age and the spatial unevenness of subsidence. The blocks constructed between 1998–2005 and after 2005 showed considerably more subsidence unevenness and temporal instability than the blocks constructed before 1998 during both time periods. The examination of the new blocks showed that the spatial unevenness increased with building volume variability. For the 16 blocks with a high building volume, variability, and subsidence unevenness, the building-scale analysis showed a positive relationship between building volume and settlement in most blocks, although the R2 was lower than 0.5. The results indicate that intense building construction in urban areas could cause differential settlement at the block scale in Beijing, while the settlement of single buildings could be influenced by the integrated effects of building volume, foundation structures, and the hydrogeological background. [ABSTRACT FROM AUTHOR]

Land subsidence is one of the most important geological hazards in Beijing, China, and its scope and magnitude have been growing rapidly over the past few decades, mainly due to long-term groundwater withdrawal. Interferometric Synthetic Aperture Radar (InSAR) has been used to monitor the deformation in Beijing, but there is a lack of analysis of the long-term spatiotemporal evolution of land subsidence. This study focused on detecting and characterizing spatiotemporal changes in subsidence in the Beijing Plain by using Persistent Scatterer Interferometry (PSI) and geographic spatial analysis. Land subsidence during 2003–2015 was monitored by using ENVISAT ASAR (2003–2010), RADARSAT-2 (2011–2015) and TerraSAR-X (2010–2015) images, with results that are consistent with independent leveling measurements. The radar-based deformation velocity ranged from −136.9 to +15.2 mm/year during 2003–2010, and −149.4 to +8.9 mm/year during 2011–2015 relative to the reference point. The main subsidence areas include Chaoyang, Tongzhou, Shunyi and Changping districts, where seven subsidence bowls were observed between 2003 and 2015. Equal Fan Analysis Method (EFAM) shows that the maximum extensive direction was eastward, with a growing speed of 11.30 km2/year. Areas of differential subsidence were mostly located at the boundaries of the seven subsidence bowls, as indicated by the subsidence rate slope. Notably, the area of greatest subsidence was generally consistent with the patterns of groundwater decline in the Beijing Plain. [ABSTRACT FROM AUTHOR]