# Time-Series Displacement of Land Subsidence in Fuzhou Downtown, Monitored by SBAS-InSAR Technique

Authors: Bo Hu; Bin Yang; Xingfu Zhang; Xiongle Chen; Yang Wu
Journal: Journal of Sensors
 (2019)

---

## Abstract

The land subsidence in response to the construction of underground engineering, e.g., subway, has caused many geological hazards and impeded the sustainable development of urbanization, in particular of China. Analysis of land subsidence with high temporal-spatial resolution is necessary and can help to assess the risk of geohazards. In this study, we apply the SBAS-InSAR technique to monitor the land subsidence in Fuzhou downtown after the program of metro construction. 24 scenes of X-band TerraSAR data from July 2013 to August 2015 and 32 scenes of C-band Sentinel-1 data from July 2015 to February 2018 were used in this experiment. Our results show a maximum subsidence rate of -12mm/yr, and eight subsidence funnels have been found during the observed period. After analyzing the subsidence of these regions in a long time span, it can be concluded that there are three regions which have a relatively stable disastrous subsidence effect, and there is a possibility of further intensified subsidence.


---

## Body

## 1. Introduction

Fuzhou, as the provincial capital of Fujian, China, has made an aggressive effort in expanding its urban construction engineering. As of 2018, the Fuzhou city has undergone major expansions:(1) Line 1 started construction on December 27, 2009. It has a total length of 28.8 kilometers and is connected to the central area of Jinan, Drum Tower, Taijiang and Cangshan.(2) Line 2 started construction on November 28, 2014, and is still under construction.The rapid development of underground engineering has also caused many negative impacts, such as land subsidence, soil erosion, building collapse, and underground pipeline rupture. Land subsidence is one of the main hazards of engineering construction. Urban land subsidence is a worldwide problem, and many cities in the world are endangered by land subsidence [1]. Land subsidence disasters will lead to changes in ground elevation, damage to municipal facilities and buildings, and threat to the safety of railway transportation, thus affecting the development and security of the city [2]. For the prevention and control of the impact and damage caused by land subsidence, it is urgent to establish an effective urban subway construction project area deformation monitoring and early warning system. In the field of subway construction engineering in densely populated commercial blocks, deformation monitoring technology is proposed to meet the requirements of high precision, high efficiency, and full time coverage.At present, the monitoring methods and techniques of traditional ground settlement mainly include high-precision traditional leveling, hierarchical dynamic monitoring (bedrock, stratification, groundwater dynamic monitoring), and GPS measurement technology. However, with the rapid development of urban construction, the scope of settlement has been continuously expanded and the network layout problem has been monitored, which has led to an increase in measurement difficulty: leveling point coverage gap, long measurement cycle, and waste of manpower and material resources; hierarchical dynamic monitoring of construction site is difficult to choose, one-time investment of a large amount of funds and limited outlets cloth; GPS measurement points are easily damaged and difficult to maintain.In recent years, Interferometric Synthetic Aperture Radar (InSAR) technology, having the advantages of being weather-independent and having high ground sensitivity, has received extensive attention in ground deformation measurement [3].In general, Differential InSAR (D-InSAR) technology firstly generates an interferogram containing both surface deformation information and terrain information from two radar images through traditional InSAR technology and then inverts the terrain phase based on the external DEM and removes it from the interference phase [4–7]. Finally, an interferogram containing only topographical deformation information is obtained, but, in fact, the interferogram also has an atmospheric phase effect, making it difficult to achieve high-precision measurement [8, 9].In order to solve the atmospheric phase effect, the scholars have conducted in-depth research and development on Permanent Scattering InSAR (PS-InSAR)[10, 11] and Small Baseline Subset InSAR (SBAS-InSAR)[12]. Among them, the SBAS-InSAR technology performs a suitable combination of the acquired data to obtain a series of small baseline subsets of differential interferograms, which can better overcome the atmospheric phase effects. Then the singular value decomposition (SVD) method is used to solve the deformation rate and solve the problem of rank loss in the solution process [13].Urban areas are susceptible to surface deformations caused by many factors, such as nature and human activities. In 2012, Xu Wenbin used SBAS-InSAR technology combined with GPS data and groundwater level data to study the surface time series deformation in Los Angeles, USA, and successfully revealed the surface deformation of the area from 2003 to 2009 [14]. In addition, in 2014, Estelle Chaussard used ALOS data for time series deformation studies, analyzed land subsidence problems throughout Mexico City through SBAS-InSAR technology, and successfully identified land subsidence in 21 regions (including 17 cities). Mexico City has a linear vertical speed of more than 30 cm/year, while other places detect a speed of 5-10 cm/year [15]. Similarly, Shaochun Dong used SBAS-InSAR to perform time series deformation analysis of Shanghai ground subsidence. The results showed that most of the study area sank evenly, but along the subway lines, elevated roads and highways, the skyscrapers quickly subsided around them [2].In summary, although relevant scholars have proved the scientific reliability of InSAR technology on urban ground deformation, the relevant research results are mostly short-term deformation monitoring based on single-satellite InSAR technology, and it is difficult to continuously monitor the construction process with high precision.In this research, two different radar satellite data, TerraSAR and Sentinel-1A, performed long-term surface deformation monitoring for Fuzhou downtown from 2013 to 2018, increasing the monitoring frequency and expanding the monitoring range. Combined with the construction of Fuzhou subway construction and related disaster information, the results of surface deformation during the construction of Fuzhou subway were studied to evaluate the feasibility and practicability of SBAS-InSAR technology in urban building deformation monitoring and to prove that this technology can be one of the convenient and effective monitoring methods in the construction and maintenance of urban buildings.


## 2. Study Area and Dataset

### 2.1. Study Area

The study area is located in the core area of Fuzhou downtown in Fujian province (Gulou District, Jinan District, Taijiang District, Cangshan District); these areas are with dense high buildings, large population flow, complex geological conditions, and traffic facilities. The main construction and operation of line 1 and line 2 were carried out during this period. The addition to focusing on the areas where the two metro lines pass through, the study also concerned on the urban areas affected by various large-scale urban construction engineering. The coverage of radar image data is shown in Figure1.Figure 1
The coverage of the radar image used in this research overlaid on Fuzhou administrative map.


### 2.2. Data

External DEM data, the “ALOS World 3D-30m” data with a spatial resolution of 30 meters provided by the Japan Aerospace Agency (JAXA), is used to simulate and eliminate terrain phase contributions. The DEM data is also used to reshape and geocode the resulting InSAR product from the running Ge-Doppler coordinates to the map geometry corresponding to the Universal Transverse Mercator (UTM) coordinate system.In this research, 24 scene TerraSAR-X images and 32 scene Sentinel-1A images were selected as research data sources. The former covers July 18, 2013, to August 27, 2015, and the latter covers the period from July 6, 2015, to February 8, 2018.Based on the selection principle of short baseline set, the optimal space baseline value and the optimal time baseline value are established, which not only ensures the continuity of time and space interference, but also ensures the correlation of interference results. The spatial baseline value of this experiment was 0-600 m, and the time baseline value was 0-200 days. The system automatically selects the super-master image in the master image in the small baseline subset according to the three-baseline minimum principle (spatial baseline, time baseline and Doppler centroid frequency baseline). Finally, 62 pairs of effective interference pairs were obtained between July 18, 2013, and August 27, 2015. 162 pairs of effective interference pairs were obtained between July 6, 2015, and February 8, 2018 (see Figure2).Figure 2
Time-baseline and time-position plots relevant to the radar images and interferograms used in this study. (a) refers to the TerraSAR-X (201307-201508) and (b) refers to the Sentinel-1(201507-201802).
(a) (b)


## 2.1. Study Area

The study area is located in the core area of Fuzhou downtown in Fujian province (Gulou District, Jinan District, Taijiang District, Cangshan District); these areas are with dense high buildings, large population flow, complex geological conditions, and traffic facilities. The main construction and operation of line 1 and line 2 were carried out during this period. The addition to focusing on the areas where the two metro lines pass through, the study also concerned on the urban areas affected by various large-scale urban construction engineering. The coverage of radar image data is shown in Figure1.Figure 1
The coverage of the radar image used in this research overlaid on Fuzhou administrative map.


## 2.2. Data

External DEM data, the “ALOS World 3D-30m” data with a spatial resolution of 30 meters provided by the Japan Aerospace Agency (JAXA), is used to simulate and eliminate terrain phase contributions. The DEM data is also used to reshape and geocode the resulting InSAR product from the running Ge-Doppler coordinates to the map geometry corresponding to the Universal Transverse Mercator (UTM) coordinate system.In this research, 24 scene TerraSAR-X images and 32 scene Sentinel-1A images were selected as research data sources. The former covers July 18, 2013, to August 27, 2015, and the latter covers the period from July 6, 2015, to February 8, 2018.Based on the selection principle of short baseline set, the optimal space baseline value and the optimal time baseline value are established, which not only ensures the continuity of time and space interference, but also ensures the correlation of interference results. The spatial baseline value of this experiment was 0-600 m, and the time baseline value was 0-200 days. The system automatically selects the super-master image in the master image in the small baseline subset according to the three-baseline minimum principle (spatial baseline, time baseline and Doppler centroid frequency baseline). Finally, 62 pairs of effective interference pairs were obtained between July 18, 2013, and August 27, 2015. 162 pairs of effective interference pairs were obtained between July 6, 2015, and February 8, 2018 (see Figure2).Figure 2
Time-baseline and time-position plots relevant to the radar images and interferograms used in this study. (a) refers to the TerraSAR-X (201307-201508) and (b) refers to the Sentinel-1(201507-201802).
(a) (b)


## 3. Method and Data Processing

SBAS-InSAR is a newly developed D-InSAR time series analysis method, first proposed by Berardino et al. [16]. It combines all the acquired SAR data to form interference pairs and then sets the spatial baseline and time baseline threshold and obtains several small baseline sets whose spatiotemporal baseline is smaller than the threshold. The surface deformation time series of each set can be solved by the least square method, and then several singular baseline sets are solved by singular value decomposition (SVD). This effectively solve the different small baseline between sets due to baseline caused by the long time discontinuity problem, improve the monitoring time resolution, effectively solve the general equation rank deficient problem, and finally get the sequence covering the whole observation time of ground subsidence.Assuming that the SAR image in the same area is N+1, the time of shooting ist0,t1,t2,t3,……tn; at least it is combined into M interference fringe images.For any two differential interference diagrams (ta<tb), the unwrapped phase can be expressed as [17]:(1)δϕx,i=ϕx,itb-ϕx,ita=ϕflat,x,i+ϕele,x,i+ϕdisp,x,i+ϕatm,x,i+ϕnoiseIn formula (1), δϕx,i is the i image and the interference phase of the x pixel. ϕx,itb and ϕx,ita is the phase of interference at two different times.Theϕflat,x,i is the phase difference of the flat phase, the ϕele,x,i is the phase difference of terrain phase, the ϕdisp,x,i is the phase difference of ground deformation, the ϕatm,x,i is the error of atmospheric phase, and ϕnoise is the error of noise.The flat phase, the phase of the terrain, the phase of the deformation, and the phase of the atmosphere in the formula (1) can be further expressed as(2)ϕflat,x,i=-4πλBnRtanθ(3)ϕele,x,i=-ΔqsinθBnR4πλ(4)ϕdisp,x,i=4πλdtb,x,i-dta,x,i(5)ϕatm,x,i=ϕatmtb,x,i-ϕatmta,x,iIn the formula,λ is the radar wavelength and R is the slope distance, Bn is the vertical baseline, Δq is the DEM elevation difference, and θ is the incidence angle.Suppose the deformation rate between two different interferograms isVi(i=1,……,n), then the cumulative deformation between ta and tb can be expressed as [18](6)ϕdisp,x,i=4πλ∑i=tatbVi,i+1ti+1-tiThe integral of the deformation rate is in the time interval between the master and slave images of the interferogram.Immediately the integration of speed is at each time interval between master image and slave image intervals. Written in matrix form:(7)BV=δϕBy using the SVD method, a generalized inverse matrix of the coefficient matrix B can be obtained, so that the minimum norm solution of the data vector can be obtained. Finally, the speed integral for each time period gives the amount of deformation for each time period.


## 4. Result

According to SBAS-InSAR technology, the annual average settlement map of Fuzhou city from July 2013 to August 2015 and from July 2015 to February 2018 was obtained (see Figure3), and the time series settlement displacement is shown in Figure 4. From the results, it can be found that there are a large number of discrete results and some abnormal uplift areas in the north of the city. This is because there are a large number of vegetation areas in the area. The growth and withering of plants tends to lead to low coherence or even complete decorrelation, which leads to unwrapping error and even unable to unwrapping. In addition, eight different large subsidence areas can be found (named I~VIII area in order, and the approximate settlement position and range are divided by two-year average settlement).Figure 3
Mean displacement velocity map over the study area derived by SBAS-InSAR technique. (a) reflects the LOS velocity of the Fuzhou downtown in 2013-2015, while (b) reflects the LOS velocity of the Fuzhou downtown in 2015-2018. The red rectangles indicate the area in which violent subsidence takes place (I-VIII). The local deformation map after enlargement can be seen in Figure6.
(a) (b)Figure 4
Interferograms of SBAS-InSAR in study area. (a) describes the time-series displacement in Fuzhou downtown in 2013-2015, and (b) shows the long-term subsidence in Fuzhou downtown in 2015-2018.After statistical analysis of the annual average settlement rate map at two different time periods, the average annual rate of Fuzhou downtown from 2013 to 2015 is +1.32mm/year, the average annual rate is between -123mm/year and +103mm/year, and the standard deviation among all pixel values is ±9.0mm/year. From 2015 to 2018, the average annual rate is -3.7mm/year, the average annual rate is -120mm/year to +80 mm/year, and the standard deviation was ±10 mm/year.It can be seen that the entire Fuzhou is in stable state from 2013 to 2018. The standard deviations of the results of the two data sources are very close, which proves that the monitoring accuracy of the two data sources is also very good. In addition, the subsidence area is basically located near the subway line, and it is speculated that the construction and operation of the subway has a huge impact on the land subsidence in Fuzhou.The interval value of the deformation point rate along the subway is calculated by the interval method. Taking the geographical location of two different subway stations as the center of the left and right boundaries of the interval, the range of influence of the selection interval is 100 meters. The rate of all the deformation points in the statistical interval is worth the average rate value of each interval of line 1 and line 2 in different periods. According to engineering geological experience, areas with deformation rate values between -5mm/year and +5mm/year in urban areas can be regarded as relatively stable areas, while deformation areas exceeding ±5mm/y years are classified as potential disaster areas. Based on the above statistical principle, the regional average rate value is obtained by counting the deformation point velocity values in the I to VIII region. The settlement rates for the first to sixth districts from 2013 to 2015 are -7 mm/year, -5.5 mm/year, -6.7 mm/year, -5.2 mm/year, -8.2 mm/year, and -12.2 mm/years, while the V, VII, and VIII zones reach -8. 7 mm/year, -5. 4 mm/year, and -7 and from 2015 to 2018, 4 mm per year (Zones I, II, III, VII, and VIII belong to Metro Line 2; Area IV belongs to Metro Line 1; Areas V and VI are road engineering areas and residential building areas).Since the two types of data vary greatly in monitoring time coverage, only one month of time overlap can be obtained. By selecting 8 relevant feature points (A01-A08) in the I-VIII area, the monitored shape variables of these points in the overlap time are analyzed for difference (see Figure5). It can be seen from Figure 5 that the two sets of data monitoring results maintained a high correlation between July 2015 and August 2015, and the difference was within 0-±2 mm. Therefore, it can be considered that the monitoring accuracy between the two sets of data in this experiment can reach the millimeter level.Figure 5
Analysis of TerraSAR data and Sentinel data shape variables from July 2015 to August 2015.Figure 6
Location of the selected deformation point on the deformation map (I-VIII).


## 5. Discussion

### 5.1. Deformation Characteristic Point

After determining the scope of the eight subsidence areas, the subsidence feature points (e.g., P1, P2, P3……P27; see in Figure6) in the area are searched to analyze the deformation characteristics of these feature points on the time series to research the subsidence trend and dynamic development.


### 5.2. Analysis

From 2013 to 2015, the area I (Yangli Station-Gushan) was in an irregular state of slow sinking, and the severe subsidence section exceeded 30 mm even in more than 24 months (see Figure7); from 2016 to 2018, the overall deformation rate of the area was -2.9 mm/year, which was relatively stable. Through data review, it was found that the construction of road works caused the destruction of geological soft soil layers, resulting in the occurrence of staged settlement effects. From Figure 7, the deformation evolution of the feature points in the time series could be found. The feature point had a stepped sedimentation effect before February 2014, and, after this time node, it was oscillating and sinking. It could be considered that before this time node, the road engineering was in the foundation pit pumping construction, which caused the stagnation of the soft soil layer on the surface. After that, the surface was in the recovery phase. However, due to irreversible damage to the surface soil structure, a certain degree of sinking effect will still occur.Figure 7
The figure status of the time-series displacement of P1, P2, and P3 within rectangle I.From 2013 to February 2014, the area II (Jinshan-Hongwan) experienced continuous subsidence with a settlement of up to 20 mm, and, after February 2014, it showed periodic uplift and settlement deformation rate of +1.6 mm/year (see Figure8). It is speculated that the reason for the settlement in the area was the pumping of groundwater, and the periodic uplift and settlement were caused by seasonal rainfall.Figure 8
The figure status of the time-series displacement of P4, P5, and P6 within rectangle II.Area III (Hongwan-Juyuanzhou) from 2013 to 2015 was between -10 mm and -20 mm, and the main settlement occurred before January 2015 (see Figure9). The subway construction plan showed that the Juyuanzhou site was in the fence construction state around January 2015. It could be reasonably considered that the excavation of the foundation pit destroyed the stable state of the nearby soil layer, resulting in uneven settlement.Figure 9
The figure status of the time-series displacement of P7, P8, and P9 within rectangle III.From 2013 to 2015, the area IV (Lulei-Triangle City) evolved between 0 mm and -30 mm, and the overall settlement trend was very consistent, which was slow settlement (see Figure10). The deformation rate of this interval is -1.5 mm/year after 2015 and was stable. During the monitoring period, the interval was in the construction phase, and the regional settlement characteristics were very close.Figure 10
The figure status of the time-series displacement of P10, P11, and P12 within rectangle IV.The sedimentation rate values of the area V between 2013 and 2015 and 2015 to 2018 exceeded 8 mm/year. The cumulative settlement in the region between 2013 and 2015 was between -15 mm and -65 mm, and between 2015 and 2018 it was between -45 mm and -80 mm. It was worth noting that the regional settlement characteristics were analyzed from the geographical distribution of the feature points. The settlement amplitude of the feature points in the northeast direction was generally larger than that of the southwest direction feature points. From 2013 to 2015, the settlement of P13 and P14 was larger than that of P15; between 2015 and 2018, the settlement of P20 and P21 was also larger than that of P19 (see Figure11).Figure 11
The figure status of the time-series displacement of P13, P14, P15, P19, P20, and P21 within rectangle V.In summary, it could be considered that the settlement trend of the area V was from the southwest to the northeast, and the subsidence area is further increased.After investigation, it was found that the area V belongs to the construction scope of the Yixu Highway Project. The excavation of the road works disturbed the structural stability of the soft soil layer and caused the influence of land subsidence. Moreover, with the segmental advancement of road engineering, the scope of settlement would be further shifted and expanded.The sedimentation rate in the area VI was the largest in all regions, reaching -12.2 mm/year between 2013 and 2015, and the cumulative deformation between -25 mm and -60 mm (see Figure12). The most serious settling funnel was formed at Longhua Apartment, where ground subsidence occurring during this period was confirmed. This was closely related to the construction of the project near the apartment. After 2015, the settlement trend would gradually stabilize and there would be no longer the impact of large-scale settlement disasters.Figure 12
The figure status of the time-series displacement of P16, P17, and P18 within rectangle VI.The accumulative subsidence of area VII (Gushan-Shangyang) from 2015 to 2018 is between -40mm and -60mm, and the main settlement occurred after April 2016 (see Figure13). This section started construction in February 2017, and the regional settlement node as in good agreement with the time. The Shangyang apartment was located on the side of the Shangyang station, and the settlement disaster of the residential building has been affected.Figure 13
The figure status of the time-series displacement of P22, P23, and P24 within rectangle VII.The deformation rate of area VIII reached -7.4mm/year from 2015 to 2018, and the subsidence point was mainly located at the intersection next to the metro internal. The settlement characteristic points showed that the cumulative settlement of this crossroads was between -30mm and -50mm, and it would continue to increase after August 2016 (see Figure14). The crossroads have collapsed badly in 2017.Figure 14
The figure status of the time-series displacement of P25, P26, and P27 within rectangle VIII.In summary, this leads to the following inference:(1) The areas VII and VIII belonged to the subway line of line 2, and the shield construction in the subway section caused the loss of the surface soil layer and the structural damage of the soil around the tunnel caused the sedimentation effect.(2) The area V belonged to the road construction area. The foundation pit excavation of the road construction destroyed the stable state of the groundwater layer, and the soil layer near the foundation pit also had the effect of sinking. The settlement effect of the Nanmendou subway station was mainly due to the excessive load.(3) The subway construction period was long, and the use of a large amount of groundwater in irregular conditions would cause surface subsidence. When the construction was completed, the sedimentation rate was obviously reduced, but the settlement was still present, which takes time from the consolidation of the rock and soil.


## 5.1. Deformation Characteristic Point

After determining the scope of the eight subsidence areas, the subsidence feature points (e.g., P1, P2, P3……P27; see in Figure6) in the area are searched to analyze the deformation characteristics of these feature points on the time series to research the subsidence trend and dynamic development.


## 5.2. Analysis

From 2013 to 2015, the area I (Yangli Station-Gushan) was in an irregular state of slow sinking, and the severe subsidence section exceeded 30 mm even in more than 24 months (see Figure7); from 2016 to 2018, the overall deformation rate of the area was -2.9 mm/year, which was relatively stable. Through data review, it was found that the construction of road works caused the destruction of geological soft soil layers, resulting in the occurrence of staged settlement effects. From Figure 7, the deformation evolution of the feature points in the time series could be found. The feature point had a stepped sedimentation effect before February 2014, and, after this time node, it was oscillating and sinking. It could be considered that before this time node, the road engineering was in the foundation pit pumping construction, which caused the stagnation of the soft soil layer on the surface. After that, the surface was in the recovery phase. However, due to irreversible damage to the surface soil structure, a certain degree of sinking effect will still occur.Figure 7
The figure status of the time-series displacement of P1, P2, and P3 within rectangle I.From 2013 to February 2014, the area II (Jinshan-Hongwan) experienced continuous subsidence with a settlement of up to 20 mm, and, after February 2014, it showed periodic uplift and settlement deformation rate of +1.6 mm/year (see Figure8). It is speculated that the reason for the settlement in the area was the pumping of groundwater, and the periodic uplift and settlement were caused by seasonal rainfall.Figure 8
The figure status of the time-series displacement of P4, P5, and P6 within rectangle II.Area III (Hongwan-Juyuanzhou) from 2013 to 2015 was between -10 mm and -20 mm, and the main settlement occurred before January 2015 (see Figure9). The subway construction plan showed that the Juyuanzhou site was in the fence construction state around January 2015. It could be reasonably considered that the excavation of the foundation pit destroyed the stable state of the nearby soil layer, resulting in uneven settlement.Figure 9
The figure status of the time-series displacement of P7, P8, and P9 within rectangle III.From 2013 to 2015, the area IV (Lulei-Triangle City) evolved between 0 mm and -30 mm, and the overall settlement trend was very consistent, which was slow settlement (see Figure10). The deformation rate of this interval is -1.5 mm/year after 2015 and was stable. During the monitoring period, the interval was in the construction phase, and the regional settlement characteristics were very close.Figure 10
The figure status of the time-series displacement of P10, P11, and P12 within rectangle IV.The sedimentation rate values of the area V between 2013 and 2015 and 2015 to 2018 exceeded 8 mm/year. The cumulative settlement in the region between 2013 and 2015 was between -15 mm and -65 mm, and between 2015 and 2018 it was between -45 mm and -80 mm. It was worth noting that the regional settlement characteristics were analyzed from the geographical distribution of the feature points. The settlement amplitude of the feature points in the northeast direction was generally larger than that of the southwest direction feature points. From 2013 to 2015, the settlement of P13 and P14 was larger than that of P15; between 2015 and 2018, the settlement of P20 and P21 was also larger than that of P19 (see Figure11).Figure 11
The figure status of the time-series displacement of P13, P14, P15, P19, P20, and P21 within rectangle V.In summary, it could be considered that the settlement trend of the area V was from the southwest to the northeast, and the subsidence area is further increased.After investigation, it was found that the area V belongs to the construction scope of the Yixu Highway Project. The excavation of the road works disturbed the structural stability of the soft soil layer and caused the influence of land subsidence. Moreover, with the segmental advancement of road engineering, the scope of settlement would be further shifted and expanded.The sedimentation rate in the area VI was the largest in all regions, reaching -12.2 mm/year between 2013 and 2015, and the cumulative deformation between -25 mm and -60 mm (see Figure12). The most serious settling funnel was formed at Longhua Apartment, where ground subsidence occurring during this period was confirmed. This was closely related to the construction of the project near the apartment. After 2015, the settlement trend would gradually stabilize and there would be no longer the impact of large-scale settlement disasters.Figure 12
The figure status of the time-series displacement of P16, P17, and P18 within rectangle VI.The accumulative subsidence of area VII (Gushan-Shangyang) from 2015 to 2018 is between -40mm and -60mm, and the main settlement occurred after April 2016 (see Figure13). This section started construction in February 2017, and the regional settlement node as in good agreement with the time. The Shangyang apartment was located on the side of the Shangyang station, and the settlement disaster of the residential building has been affected.Figure 13
The figure status of the time-series displacement of P22, P23, and P24 within rectangle VII.The deformation rate of area VIII reached -7.4mm/year from 2015 to 2018, and the subsidence point was mainly located at the intersection next to the metro internal. The settlement characteristic points showed that the cumulative settlement of this crossroads was between -30mm and -50mm, and it would continue to increase after August 2016 (see Figure14). The crossroads have collapsed badly in 2017.Figure 14
The figure status of the time-series displacement of P25, P26, and P27 within rectangle VIII.In summary, this leads to the following inference:(1) The areas VII and VIII belonged to the subway line of line 2, and the shield construction in the subway section caused the loss of the surface soil layer and the structural damage of the soil around the tunnel caused the sedimentation effect.(2) The area V belonged to the road construction area. The foundation pit excavation of the road construction destroyed the stable state of the groundwater layer, and the soil layer near the foundation pit also had the effect of sinking. The settlement effect of the Nanmendou subway station was mainly due to the excessive load.(3) The subway construction period was long, and the use of a large amount of groundwater in irregular conditions would cause surface subsidence. When the construction was completed, the sedimentation rate was obviously reduced, but the settlement was still present, which takes time from the consolidation of the rock and soil.


## 6. Conclusion

In this experiment, the land deformation of Fuzhou urban area from 2013 to 2018 was successfully performed by applying SBAS-InSAR technology combined with two different radar satellite data (TerraSAR and Sentine-l). By analyzing the difference between the two datasets of results in the overlap time, the data maintains a high degree of correlation, and the monitoring accuracy between them reaches the millimeter level.Through the method of time series analysis, the influences and development trends of various subsidence areas were discussed, some stable deformation areas were eliminated, and some potential subsidence disaster areas were also discovered, thereby preventing the occurrence and expansion of subsidence disasters. Among them, the settlement of the areas VII and VIII was caused by the shield construction, which has great reference significance for the future subway project. The settlement of the area V was caused by the excessive load. Besides, excessive use of groundwater in construction projects requires timely prevention and control measures and is prone to disasters. In summary, the feasibility and practicability of SBAS-InSAR technology in urban building deformation monitoring can be one of the convenient and effective monitoring methods for urban construction and maintenance.



---

*Source: 3162652-2019-05-19.xml*