Coastal Subsidence Monitoring Associated with Land Reclamation Using the Point Target Based SBAS-InSAR Method : A Case Study of Shenzhen , China

Shenzhen, the first special economic zone of China, has witnessed earth-shaking changes since the late 1980s. In the past 35 years, about 80 km2 of land has been reclaimed from the sea in Shenzhen. In order to investigate coastal vertical land motions associated with land reclamation, we proposed an elaborated Point Target (PT) based Small Baseline Subset InSAR (SBAS-InSAR) strategy to process an ENVISAT ASAR ascending and descending orbits dataset both acquired from 2007 to 2010. This new strategy can not only select high density PTs but also generate a reliable InSAR measurement with full spatial resolution. The inter-comparison between InSAR-derived deformation velocities from different orbits shows a good self-consistency of the results extracted by the elaborated PT-based SBAS-InSAR strategy. The InSAR measurements show that the reclaimed land is undergoing remarkable coastal subsidence (up to 25 mm/year), especially at the Shenzhen Airport, Bao’an Center, Qianhai Bay and Shenzhen Bay. By analyzing the results together with land reclamation evolution, we conclude that the ground deformation is expected to continue in the near future, which will amplify the regional sea level rise.


Introduction
The Coastal Zone [1] plays a very important role in transportation, the circulation of resources and funds, etc.To promote urbanization and economic development, land reclamation [2] from the sea is a common practice for many countries in the world, such as USA [3], Japan [4], Italy [5], the Netherlands [6,7], and China [8][9][10].
Over the past decades, Synthetic Aperture Radar Interferometry (InSAR) [11,12] has been proven to be a powerful geodetic technique, and been extensively exploited to investigate ground deformation induced, for example, by earthquake [13][14][15] or volcano activities [16,17].Furthermore, recently remarkable improvements in Multi-Temporal InSAR (MT-InSAR), including the Small Baseline Subset InSAR (SBAS-InSAR) [18,19] and the Persistent Scatterer InSAR (PS-InSAR) [20][21][22][23], enable us to monitor long period surface deformation caused by permafrost [24], ground water extraction [25][26][27][28], mining [29], and land reclamation [30].The MT-InSAR exploits a set of high phase-quality pixels [20], other than the original full two-dimensional image grid, to map the deformation history and the corresponding average deformation velocity.It can overcome most intrinsic temporal and geometric decorrelations [25] of conventional InSAR, and simultaneously estimate the atmospheric effect [26][27][28][29] to generate high accuracy deformation results.In recent years, researches that investigated the ground deformation related to land reclamation by MT-InSAR have been carried out in reclaimed lands around the world, such as the Chek Lap Kok Airport in Hong Kong [30], Sibari in Southern Italy [31], and the Venice coastland [32].Besides the MT-InSAR techniques, geodetic techniques [33], such as GPS and leveling, have been widely used in ground deformation monitoring.However, the MT-InSAR could provide a better spatial resolution to capture the deformation details compared with GPS or leveling.
Shenzhen is a city built on the weak alluvial clay [34] of sludge and reclaimed land.Due to high compressibility, low permeability, and weak strength of the alluvial clay [35], the soil becomes compacted under the pressure from constructions and clay self-consolidation.The compaction process causes ground deformation [36], threatening the stability of ground constructions and underground man-made facilities [37].Therefore, ground deformation monitoring is very important for the reclaimed land utilization in Shenzhen.Dense Point Targets (PTs) for deformation monitoring can be selected from built-up area [20] that maintains long-time high coherence, while for the non-built-up area, e.g., reclaimed land, for the area surrounding buildings built on the reclaimed land, few PTs can be selected [38].Although previous studies [36,39] have tried to capture the deformation in urban areas with MT-InSAR, the following two issues remain unsolved: (1) When traditional MT-InSAR methods are used, the ground deformation is not well determined due to the low PTs density in the reclaimed land; (2) The deformation associated with land reclamation has not yet been analyzed.So how to select the appropriate number of usable PTs for MT-InSAR deformation monitoring and the analysis of the deformation results needs to be further studied for the reclaimed land of Shenzhen.
In this study, an elaborated PT-based SBAS-InSAR strategy is proposed to select usable PTs and generate an InSAR measurement with full spatial resolution.Based on the selection of PTs, two independent ENVISAT ASAR datasets along the ascending and descending orbits were used to map the coastal land deformation in Shenzhen between 2007 and 2010.Inter-comparison of the InSAR measurements from different orbits was carried out to check the self-consistency of the measurements.Finally, we discussed the implications and potential subsidence hazards caused by land reclamation in this area.

Geological Environment
Shenzhen is located at the southeast of the Pearl River Delta (PRD) region in China.The city has complicated topography and is surrounded by the sea in the south and the west.Its elevation ranges between 1 and 10 m in the southwest of Shenzhen, reaches 30 m in the heart of the city, and goes up to 80 m at the northeast (see Figure 1).The marine sediments of the Quaternary period widely spread across the low-lying coastal land areas.This sediment layer is of high water content and high compressibility, big void ratio, and weak shearing strength [35].The thickness of this layer is generally 3-10 m, but is up to 20 m in some particular areas [40].

Historical Land Reclamation
Due to its geographical location, Shenzhen is the direct link between the China mainland and Hong Kong.With its economic development and boom [43], the demand for land use [44] is increasing.The comparison of Landsat satellite images (see Figure A1) shows the dramatic landscape changes in Shenzhen from 1979 to 2015.Small hills were leveled, and thousands of high buildings have replaced previous agricultural and vegetated areas.Meanwhile, 80 km 2 of land has been reclaimed from the sea along the coastal zone in the past 35 years [45].The temporal and spatial evolution of land reclamation [43] is shown in Figure 2. Land reclamation has brought considerable economic benefits [46] to the government and led to further economic development in Shenzhen.However, long-time compaction of the reclaimed land causes ground deformation [30], which affects the stability of the surface and underground constructions, thereby threatens the safety of people's lives and assets.
According to an investigation conducted by the Geological Bureau of Shenzhen, Shenzhen is in a stable tectonic environment with low tectonic hazard risk [41].However, local ground deformation [1,42] related to anthropogenic activity has become the major geohazard in Shenzhen, especially in the coastal reclaimed areas.

Historical Land Reclamation
Due to its geographical location, Shenzhen is the direct link between the China mainland and Hong Kong.With its economic development and boom [43], the demand for land use [44] is increasing.The comparison of Landsat satellite images (see Figure A1) shows the dramatic landscape changes in Shenzhen from 1979 to 2015.Small hills were leveled, and thousands of high buildings have replaced previous agricultural and vegetated areas.Meanwhile, 80 km 2 of land has been reclaimed from the sea along the coastal zone in the past 35 years [45].The temporal and spatial evolution of land reclamation [43] is shown in Figure 2. Land reclamation has brought considerable economic benefits [46] to the government and led to further economic development in Shenzhen.However, long-time compaction of the reclaimed land causes ground deformation [30], which affects the stability of the surface and underground constructions, thereby threatens the safety of people's lives and assets.
According to an investigation conducted by the Geological Bureau of Shenzhen, Shenzhen is in a stable tectonic environment with low tectonic hazard risk [41].However, local ground deformation [1,42] related to anthropogenic activity has become the major geohazard in Shenzhen, especially in the coastal reclaimed areas.

Dataset
Two sets of ENVISAT ASAR data acquired between 2007 and 2010, along the ascending (31 scenes, see Table A1 for detail) and descending orbits (28 scenes, see Table A2 for detail), were employed to analyze the coastal deformation.Two dense networks of SBAS interferometric pairs (see Figure 3) with both short perpendicular and temporal baselines were used in the data processing.The available multi-orbit datasets allow us to validate the self-consistency of the monitoring results by inter-comparison.

Dataset
Two sets of ENVISAT ASAR data acquired between 2007 and 2010, along the ascending (31 scenes, see Table A1 for detail) and descending orbits (28 scenes, see Table A2 for detail), were employed to analyze the coastal deformation.Two dense networks of SBAS interferometric pairs (see Figure 3) with both short perpendicular and temporal baselines were used in the data processing.The available multi-orbit datasets allow us to validate the self-consistency of the monitoring results by inter-comparison.

Dataset
Two sets of ENVISAT ASAR data acquired between 2007 and 2010, along the ascending (31 scenes, see Table A1 for detail) and descending orbits (28 scenes, see Table A2 for detail), were employed to analyze the coastal deformation.Two dense networks of SBAS interferometric pairs (see Figure 3) with both short perpendicular and temporal baselines were used in the data processing.The available multi-orbit datasets allow us to validate the self-consistency of the monitoring results by inter-comparison.

PT-Based SBAS-InSAR Strategy
In our study area, due to the vegetation covering and humid subtropical climate, the application of SAR interferometry is mainly limited by temporal decorrelation [25,47] and inhomogeneity in the tropospheric path delay [27,28].Besides, the land reclamation related deformation has limited spatial range, thus we do not multi-look the interferograms to retain deformation details.We processed the datasets with an elaborated PT-based SBAS-InSAR, and below (Figure 4) is a flow chart showing the strategy, which consists of four major processing steps.

PT-Based SBAS-InSAR Strategy
In our study area, due to the vegetation covering and humid subtropical climate, the application of SAR interferometry is mainly limited by temporal decorrelation [25,47] and inhomogeneity in the tropospheric path delay [27,28].Besides, the land reclamation related deformation has limited spatial range, thus we do not multi-look the interferograms to retain deformation details.We processed the datasets with an elaborated PT-based SBAS-InSAR, and below (Figure 4) is a flow chart showing the strategy, which consists of four major processing steps.

Pre-Processing
Data processing starts with a stack of Single Look Complex (SLC) images [48].To minimize the effect of mis-coregistration caused by geometric, temporal, and Doppler decorrelations [25], the images acquired on 28 January 2009 and 17 August 2008 for track 025 and track 175, respectively, were chosen as geometry references for SLC coregistration using the Formula (4.1) of reference [49].Interferometric pairs with both short perpendicular and temporal baselines were selected to construct networks of the SBAS interferometric pairs.Also, two dense networks of SBAS interferometric pairs (for details see Table 1 and Figure 3) were established.The topographic phase in each interferogram was removed with the 3-arc SRTM DEM [50].To preserve the detail deformation pattern, the multilooking process was not applied to any interferogram pairs.A reference point, black star in Figures 5 and 6, was chosen from an inland low-rise building in a stable area.

Pre-Processing
Data processing starts with a stack of Single Look Complex (SLC) images [48].To minimize the effect of mis-coregistration caused by geometric, temporal, and Doppler decorrelations [25], the images acquired on 28 January 2009 and 17 August 2008 for track 025 and track 175, respectively, were chosen as geometry references for SLC coregistration using the Formula (4.1) of reference [49].Interferometric pairs with both short perpendicular and temporal baselines were selected to construct networks of the SBAS interferometric pairs.Also, two dense networks of SBAS interferometric pairs (for details see Table 1 and Figure 3) were established.The topographic phase in each interferogram was removed with the 3-arc SRTM DEM [50].To preserve the detail deformation pattern, the multi-looking process was not applied to any interferogram pairs.A reference point, black star in Figures 5 and 6, was chosen from an inland low-rise building in a stable area.

Improved Point Target Selection
Recently, many researches demonstrated that exploiting point-like scatterers [20,21] instead of a two-dimensional grid for deformation monitoring can overcome the effect of interferometric decorrelation [25,51].The selection of more usable PTs in a low-coherence area [21] for deformation retrieval is a key issue in MT-InSAR data processing.In our study, the coastal reclaimed land may undergo significant deformation and the backscattering coefficients may vary significantly in the time domain.As a consequence, these two factors limit the application of the dispersion of amplitude (DA) method [20] or the distributed scatterers (DS) method [38] to select sufficient usable PTs.Besides, as a phase quality indicator, the coherence maps associated with the interferometric phase could be exploited to identify usable targets.Point target with coherence greater than a given threshold could be selected as a usable PT.However, the coherence is underestimated, as the phase values are non-stationary.Thus non-stationary signals, including DEM, atmospheric effect [26] and possibly the target motion should be estimated and thus be removed.The space averaging of the data inside a certain window (e.g., 50 m × 50 m) is used as the estimation of these non-stationary signals [21] for the center pixel of the window.However, this is not a common case because the average value may not be the genuine estimate [51].
To select adequate usable PTs, an elaborated strategy based on Spectral Diversity of SLC images [48] and the SBAS network [18] is exploited.The concept of spectral diversity originates from Interferometric Point Target Analysis (IPTA) [52], whose principle is that the energy of usable PTs remains almost the same when processing different looks with fractional azimuth and different range bandwidths.We calculated the spectral diversity map γ for each co-registered SLC image, and selected a point target as a usable PT, which has the weighted average of the time series spectral diversity values greater than a given threshold (e.g., γ > 0.38).Given M SLC images, and N interferometric pairs in the SBAS network, then γ can be expressed as where γ m (r, a) is the spectral diversity value at pixel location (r, a) (range, azimuth) for the m-th SLC image, n m denotes the frequency of the m-th SLC image used in the SBAS network, and Our PT selection strategy is consistent with the one based on the mean interferometric coherence of interferogram stack for SBAS-InSAR.In the coherence-based method, the more frequently a SLC image is used to combine an interferometric pair, the larger the weight factor would be assigned to this image in calculating the mean interferometric coherence for the SBAS network interferogram stack.Similarly, using the frequency of a SLC image as the weight factor in the calculation of mean spectral diversity map γ, would also adjust the influence of that image's spectral diversity map accordingly.The spectral diversity map could be treated as the Signal-to-Noise Ratio (SNR) indicator of the SLC pixel's phase.Two pixels with high spectral diversity values would generate high quality interferometric phases, and vice versa.As a consequence, the weighted spectral diversity map γ is the equivalence of the SNR of the interferometric phases for the given SBAS network, which allows a more accurate PT selection in SBAS-InSAR.
After the PT selection, the data information, including the interferometric phases, spatial coherence etc., is extracted from the two-dimensional grid and stored in vector format for subsequent data processing.

Phase Filtering, Orbital Phase Error Removal and Phase Unwrapping
In order to reduce the phase noise, the differential interferometric phase for PTs, was spatially filtered by an improved SUSAN filter [53,54] with a search window size of 11 × 11 (range × azimuth) pixels.This method can avoid noise contamination from unusable PTs.The spatial coherence for each PT was estimated with the same window.Subsequently, the smoothed phase of each interferometric pair from the SBAS network was unwrapped by the method based on network programming [55] with the spatial coherence as the weight factor.
In InSAR processing, the flat-earth phase component is removed with precise orbits.For ENVISAT, it has an accuracy of 5-7 cm radially and 10-15 cm cross-track [56], and its residual orbital phase error [57] with long wavelength character may remain in the interferogram.Such error is easy to discriminate from the local deformation in our study area.Hence, we estimated the orbital error using quadratic polynomial models [57] and removed it the error from the unwrapped interferogram.It needs to be noted that the quadratic models would also absorb long-wavelength components of atmospheric effect [27,28].

Modeling of Deformation and Unwrapping Error Checking
The Singular Value Decomposition (SVD) method was employed to solve the rank-deficient problem [18] in SBAS-InSAR deformation modeling.The raw time series deformations that are contaminated by atmospheric delay and phase noise and DEM error correction can be obtained based on the unwrapped phase obtained in Section 3.2.2.The reliability of phase unwrapping was evaluated via the phase difference between the unwrapped phase and the one simulated with raw time series deformation, especially for Qianhai Bay and Shenzhen Bay where sparsely distributed PTs are identified.If the interferometric phase was unwrapped with error, then the phase difference would be significantly large.We then computed the histograms of phase difference for statistics (see Tables A3 and A4).We can see that over 82% of PTs were unwrapped with an error less than 0.5 rad for Qianhai Bay, and 83.3% and 86.4% of the PTs were unwrapped for track 025 and track 175.For these two areas, over 95% of the PTs were unwrapped with error less than 1.0 rad.Additionally, the result showed that over 95% of PTs were unwrapped with error less than 1.0 rad for the whole study area.As a result, we believe that the interferometric phase was unwrapped without significant error.
In order to remove the error of atmospheric delay and phase noise from the raw time series deformations, the linear regression with the raw time series deformations was firstly carried out to inverse the mean deformation velocity map [18].Then the inverted linear deformation was subsequently subtracted from the raw time series deformations calculated previously.The remaining components are the non-linear ground deformation, atmospheric delay, and random phase noise.Non-linear ground deformations might have some degree of correlation (low-pass) in the temporal domain, while atmospheric delay and phase noise are both high-pass, thus non-linear deformation can be separated by the temporal triangle weighted low-pass filter [20] with a certain length of time window (e.g., 180 days).Finally, we obtained the total deformation by adding up the linear deformation and non-linear deformation.

Results
Due to the side-looking geometry, the InSAR measurement is the sum of the projections of ground three-dimensional (3-D) ground deformations in the (Line-of-Sight) LOS direction.In Shenzhen coastal zone, the subsidence of the reclaimed land and sediment layer dominates the ground deformation.Hence, the deformation in our study area was assumed mainly in the vertical direction and the LOS displacement was directly back-projected into the vertical direction using the local incidence angle [36].

Deformation Rate Maps
The deformation rate maps generated from the ascending and the descending data are between -25 mm/year and 15 mm/year, as shown in Figures 5 and 6, respectively.The deformation patterns are very similar, suggesting that vertical subsidence dominates the coastal deformation.The deformation rate maps demonstrate that most of the inland area within the coastline in 1979 is stable, which means the selected reference point did not undergo obvious deformation.However, the reclaimed land area along the coastline in 1979 (the yellow line in Figures 5 and 6) experienced significant subsidence, especially at Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay, marked by black rectangles in Figures 5 and 6.Obvious deformation was not observed in the Chiwan and Shekou harbors, because the foundation of the wharf apron space was enhanced by using special rammer of small base [35].Besides, we also observed non-uniform deformation in the mean deformation rate maps.The phenomenon is particularly evident in the developed urban region B, Bao'an Center, where many roads and high rise buildings cover the reclaimed land.In addition to the ground subsidence, there are several patches with remarkable uplift located around the right bottom corner in region A. The positive deformation signal may be related to the rise of groundwater level [46].This area is an old industrial district, where a great amount of groundwater was extracted for production.Now the factories are moving out, decreasing the groundwater extraction [46].As a result, the groundwater level is recovering.
To validate the non-uniform subsidence, field work was carried out at Bao'an Center (around point TP4) in 2011.Figure 7a,b show remarkable effects of the differential subsidence rates observed at a road and a building.This is because they were built on different foundations.The black arrows in Figure 7 point to the side with solid foundations, which may have slight subsidence.While, red arrows point to the side with unreinforced foundations, where remarkable subsidence occurred.
Remote Sens. 2016, 8, 652 9 of 22 significant subsidence, especially at Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay, marked by black rectangles in Figures 5 and 6.Obvious deformation was not observed in the Chiwan and Shekou harbors, because the foundation of the wharf apron space was enhanced by using special rammer of small base [35].Besides, we also observed non-uniform deformation in the mean deformation rate maps.The phenomenon is particularly evident in the developed urban region B, Bao'an Center, where many roads and high rise buildings cover the reclaimed land.In addition to the ground subsidence, there are several patches with remarkable uplift located around the right bottom corner in region A. The positive deformation signal may be related to the rise of groundwater level [46].This area is an old industrial district, where a great amount of groundwater was extracted for production.Now the factories are moving out, decreasing the groundwater extraction [46].As a result, the groundwater level is recovering.
To validate the non-uniform subsidence, field work was carried out at Bao'an Center (around point TP4) in 2011.Figure 7a,b show remarkable effects of the differential subsidence rates observed at a road and a building.This is because they were built on different foundations.The black arrows in Figure 7 point to the side with solid foundations, which may have slight subsidence.While, red arrows point to the side with unreinforced foundations, where remarkable subsidence occurred.

Time Series of Selected Points
Several points were selected to show the time series deformation (see Figure 8).The image dataset start time is usually taken as the reference time of the time series in the SBAS-InSAR technique.As seen in Table 1, the two datasets have different start times.For comparison, we set one reference time for these two datasets as 28 February 2007.As a consequence, the time series deformations mapped in the descending orbit was shifted with a constant, which was calculated according to the deformation rate and time difference (116 days).Clearly, in the reclaimed land, TP1, TP2, TP4, PT5, and TP6, are subsiding linearly (see Figure 8).TP3, located inland according to the coastline in 1979, has a different temporal pattern, showing that the TP3 was in rapid deformation (−31.7 mm/year) before January 2009, then its deformation was slowed down (−18.6 mm/year) between January 2009 and January 2010, and finally the area became stable.This deformation is very consistent with the pattern from predictions of the consolidation theory of unsaturated soil [58].Although the subsidence rate at each selected point is slightly different, the time series obtained from the ascending and descending orbits matched well.

Time Series of Selected Points
Several points were selected to show the time series deformation (see Figure 8).The image dataset start time is usually taken as the reference time of the time series in the SBAS-InSAR technique.As seen in Table 1, the two datasets have different start times.For comparison, we set one reference time for these two datasets as 28 February 2007.As a consequence, the time series deformations mapped in the descending orbit was shifted with a constant, which was calculated according to the deformation rate and time difference (116 days).Clearly, in the reclaimed land, TP1, TP2, TP4, PT5, and TP6, are subsiding linearly (see Figure 8).TP3, located inland according to the coastline in 1979, has a different temporal pattern, showing that the TP3 was in rapid deformation (−31.7 mm/year) before January 2009, then its deformation was slowed down (−18.6 mm/year) between January 2009 and January 2010, and finally the area became stable.This deformation is very consistent with the pattern from predictions of the consolidation theory of unsaturated soil [58].Although the subsidence rate at each selected point is slightly different, the time series obtained from the ascending and descending orbits matched well.

Self-Consistency Checking by Inter-Comparison between the Ascending and Descending InSAR Measurements
In order to quantitatively check the self-consistency of InSAR measurements, inter-comparison of the deformation rates between ascending and descending results was carried out.Due to the uncertainty in SAR geo-location [59], the PTs from different orbits were unlikely to be at the exact same locations.The resolution of the results is down-sampled to reduce the effect of the geo-location uncertainty [60].Firstly, a data grid with a resolution of 50 m × 50 m was created with the tradeoff between the deformation detail and resolutions.Then InSAR measurements from the ascending and descending orbits were resampled respectively to the defined data grid by Pixel Aggregate [61].Finally, the InSAR measurements from different orbits were compared.We assumed that the InSAR measurements on the defined grid had a common deformation signal.The inter-comparison results between the ascending and descending InSAR-derived deformation rates showed that correlation between the two InSAR measurements is 0.92, revealing that the InSAR measurements mapped by our PT-based SBAS-InSAR strategy are in good agreement with each other.The Root-Mean-Square Error (RMSE) of the difference between the ascending and descending orbits is 1.8 mm/year, implying that the InSAR measurements are highly self-consistent.However, there are points lying

Self-Consistency Checking by Inter-Comparison between the Ascending and Descending InSAR Measurements
In order to quantitatively check the self-consistency of InSAR measurements, inter-comparison of the deformation rates between ascending and descending results was carried out.Due to the uncertainty in SAR geo-location [59], the PTs from different orbits were unlikely to be at the exact same locations.The resolution of the results is down-sampled to reduce the effect of the geo-location uncertainty [60].Firstly, a data grid with a resolution of 50 m × 50 m was created with the tradeoff between the deformation detail and resolutions.Then InSAR measurements from the ascending and descending orbits were resampled respectively to the defined data grid by Pixel Aggregate [61].Finally, the InSAR measurements from different orbits were compared.We assumed that the InSAR measurements on the defined grid had a common deformation signal.The inter-comparison results between the ascending and descending InSAR-derived deformation rates showed that correlation between the two InSAR measurements is 0.92, revealing that the InSAR measurements mapped by our PT-based SBAS-InSAR strategy are in good agreement with each other.The Root-Mean-Square Error (RMSE) of the difference between the ascending and descending orbits is 1.8 mm/year, implying that the InSAR measurements are highly self-consistent.However, there are points lying away from the fitted model (the purple line in Figure 9), which may be caused by the errors in InSAR measurements [57], errors in SAR geo-location [59] (mismatch of the measurements from different viewing geometries), horizontal displacement and possibly the non-linear deformation process [20].away from the fitted model (the purple line in Figure 9), which may be caused by the errors in InSAR measurements [57], errors in SAR geo-location [59] (mismatch of the measurements from different viewing geometries), horizontal displacement and possibly the non-linear deformation process [20].

Deformation Associated with Land Reclamation
As the first special economic zone of China, Shenzhen has witnessed great changes since the 1980s, accompanied by extensive land reclamation over the past 30 years.The InSAR measurements of regions A, B, C, and D (Figures 5 and 6) are zoomed in Figures 10-13.According to the InSAR measurements, we can observe remarkable deformation in most of the reclaimed land, and the maximum rate is up to 25 mm/year.Generally, ground deformation is induced by primary compaction and secondary compression of the alluvial clay beneath the reclamation [30].Previous investigations [30,58] showed that 90% primary compaction of alluvial clay in Hong Kong Chek Lap Kok Airport, a site that has very similar land reclamation activity to our study area, took about 10 to 20 years, and longer for a thickness sediment layer.As shown in Figure 2, about 34.3 km 2 of land, accounting for 42.6% of the total, was reclaimed from the sea during the past 15 years (period of [2000][2001][2002][2003][2004][2005][2006][2007][2008][2009][2010][2011][2012][2013][2014][2015], on which large infrastructure has been built.Thus, ground settlement is expected in most of the reclaimed area. Figure 10 shows different subsidence patterns for the new airport built on reclaimed land and the old airport built inland.The new airport was built between years 2000 and 2010 (See Figure A1), thus few PTs were selected from that region.Numerous PTs were identified along the dam to the west part of the new airport.Although the new airport was reclaimed gradually (see the coastline change in Figure 10, the yellow line, yellow dash line, magenta dash line, and the white dash line), remarkable deformation was observed during different phases.For region B (Figure 11), Bao'an center, where many high-rise buildings stand, remarkable subsidence signals in the reclaimed land and along the Metro Lines were detected in both ascending and descending datasets (Figure 11).The Bao'an Center was mainly reclaimed in the period 1979-1990 (the area along the coastlines between 1979 and 1990, see Figure 11).However, it still undergoes remarkable ground deformation, which may be induced by engineering construction.Slight subsidence was observed on the newly reclaimed land (the area outside the coastline of 1990), as little urban construction is in progress.If a new construction were to be built, the temporary equilibrium could be broken, which may cause

Deformation Associated with Land Reclamation
As the first special economic zone of China, Shenzhen has witnessed great changes since the 1980s, accompanied by extensive land reclamation over the past 30 years.The InSAR measurements of regions A, B, C, and D (Figures 5 and 6) are zoomed in Figures 10-13.According to the InSAR measurements, we can observe remarkable deformation in most of the reclaimed land, and the maximum rate is up to 25 mm/year.Generally, ground deformation is induced by primary compaction and secondary compression of the alluvial clay beneath the reclamation [30].Previous investigations [30,58] showed that 90% primary compaction of alluvial clay in Hong Kong Chek Lap Kok Airport, a site that has very similar land reclamation activity to our study area, took about 10 to 20 years, and longer for a thickness sediment layer.As shown in Figure 2, about 34.3 km 2 of land, accounting for 42.6% of the total, was reclaimed from the sea during the past 15 years (period of 2000-2015), on which large infrastructure has been built.Thus, ground settlement is expected in most of the reclaimed area.
Figure 10 shows different subsidence patterns for the new airport built on reclaimed land and the old airport built inland.The new airport was built between years 2000 and 2010 (See Figure 11), thus few PTs were selected from that region.Numerous PTs were identified along the dam to the west part of the new airport.Although the new airport was reclaimed gradually (see the coastline change in Figure 10, the yellow line, yellow dash line, magenta dash line, and the white dash line), remarkable deformation was observed during different phases.For region B (Figure 11), Bao'an center, where many high-rise buildings stand, remarkable subsidence signals in the reclaimed land and along the Metro Lines were detected in both ascending and descending datasets (Figure 11).The Bao'an Center was mainly reclaimed in the period 1979-1990 (the area along the coastlines between 1979 and 1990, see Figure 11).However, it still undergoes remarkable ground deformation, which may be induced by engineering construction.Slight subsidence was observed on the newly reclaimed land (the area outside the coastline of 1990), as little urban construction is in progress.If a new construction were to be built, the temporary equilibrium could be broken, which may cause additional ground deformation.The Qianhai Bay, an unexploited reclaimed land area, has a few sites with high deformation rates (see Figure 12) caused by compaction of the reclaimed foundation and the alluvial soil beneath.The detail subsidence rate maps of Shenzhen Bay are shown in Figure 13.The results revealed that the body of Shenzhen Bay Check Point is in a weak deformation state ([−6, −3] mm/year), while evident settlements were observed in some sites of the reclaimed land, especially the north part of the Check Point.High deformation rates were also detected along the Shenzhen Metro Line (black dash line in Figure 13), which was centered on the reclaimed land.Based on these, we can conclude that remarkable subsidence has occurred on the lands reclaimed 15 years before, such as the Bao'an Center.
Remote Sens. 2016, 8, 652 12 of 22 additional ground deformation.The Qianhai Bay, an unexploited reclaimed land area, has a few sites with high deformation rates (see Figure 12) caused by compaction of the reclaimed foundation and the alluvial soil beneath.The detail subsidence rate maps of Shenzhen Bay are shown in Figure 13.
The results revealed that the body of Shenzhen Bay Check Point is in a weak deformation state ([−6, −3] mm/year), while evident settlements were observed in some sites of the reclaimed land, especially the north part of the Check Point.High deformation rates were also detected along the Shenzhen Metro Line (black dash line in Figure 13), which was centered on the reclaimed land.Based on these, we can conclude that remarkable subsidence has occurred on the lands reclaimed 15 years before, such as the Bao'an Center.
In addition, the engineering constructions, e.g., the construction of Shenzhen Airport and urbanization of Bao'an, would also contribute to subsidence.Considering the deformation maps and time series at selected points, we believe that the four major land reclamation regions are undergoing remarkable subsidence, which can be expected to continue in the near future.Work to further monitor the ground subsidence shall focus on the coastal land reclamation region, especially Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay.additional ground deformation.The Qianhai Bay, an unexploited reclaimed land area, has a few sites with high deformation rates (see Figure 12) caused by compaction of the reclaimed foundation and the alluvial soil beneath.The detail subsidence rate maps of Shenzhen Bay are shown in Figure 13.
The results revealed that the body of Shenzhen Bay Check Point is in a weak deformation state ([−6, −3] mm/year), while evident settlements were observed in some sites of the reclaimed land, especially the north part of the Check Point.High deformation rates were also detected along the Shenzhen Metro Line (black dash line in Figure 13), which was centered on the reclaimed land.Based on these, we can conclude that remarkable subsidence has occurred on the lands reclaimed 15 years before, such as the Bao'an Center.
In addition, the engineering constructions, e.g., the construction of Shenzhen Airport and urbanization of Bao'an, would also contribute to subsidence.Considering the deformation maps and time series at selected points, we believe that the four major land reclamation regions are undergoing remarkable subsidence, which can be expected to continue in the near future.Work to further monitor the ground subsidence shall focus on the coastal land reclamation region, especially Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay.In addition, the engineering constructions, e.g., the construction of Shenzhen Airport and urbanization of Bao'an, would also contribute to subsidence.Considering the deformation maps and time series at selected points, we believe that the four major land reclamation regions are undergoing remarkable subsidence, which can be expected to continue in the near future.Work to further monitor the ground subsidence shall focus on the coastal land reclamation region, especially Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay.Other remakes are the same as in Figure 10.Other remakes are the same as in Figure 10.

Effects of Coastal Subsidence Coupled with Sea Level Changes
The subsidence of the reclaimed land threatens infrastructure, and contributes to the relative sea level changes [62], which may cause saltwater intrusion [63,64] and destroy coastal biodiversity [37].Assuming a mean sea level rise of 2.5 mm/year [33,65], the selected regions A, B, C, and D (shown in Figure 10-13) would fall to sea level in 114, 121, 225, and 313 years (see Table 2).Compared with Other remakes are the same as in Figure 10.Other remakes are the same as in Figure 10.

Effects of Coastal Subsidence Coupled with Sea Level Changes
The subsidence of the reclaimed land threatens infrastructure, and contributes to the relative sea level changes [62], which may cause saltwater intrusion [63,64] and destroy coastal biodiversity [37].Assuming a mean sea level rise of 2.5 mm/year [33,65], the selected regions A, B, C, and D (shown in

Effects of Coastal Subsidence Coupled with Sea Level Changes
The subsidence of the reclaimed land threatens infrastructure, and contributes to the relative sea level changes [62], which may cause saltwater intrusion [63,64] and destroy coastal biodiversity [37].Assuming a mean sea level rise of 2.5 mm/year [33,65], the selected regions A, B, C, and D (shown in Figures 10-13) would fall to sea level in 114, 121, 225, and 313 years (see Table 2).Compared with previous studies [38], our deformation rate is slightly larger due to the fact that more PTs have been selected from significant deformation areas in our method.As a consequence, the large deformation rate areas have no measurement of their results whereas appropriate PTs have been identified by our method for deformation mapping.So the difference between their results and ours is expected.The risk assessment based only on sea level rise had been underestimated in the Shenzhen coastal zone because of the amplification effect from coastal subsidence.To decelerate the relative sea level rise, effective measures, such as reinforcing the ground foundation, need to be taken to prevent subsidence in the coastal zone [66].

Potential and Limits of InSAR in Mapping Coastal Subsidence
The InSAR measurements from both ascending (track 025) and descending (track 175) orbits showed a very similar pattern in the deformation rate maps (Figures 5, 6 and 10-13) derived by PT-based SBAS-InSAR, and the displacement time-series of the selected points (Figure 8).In order to properly map the coastal deformation associated with land reclamation, the effect of low-coherence will be considered in the implementation of data processing.To illustrate this problem, the dataset from track 175 was also processed with the STAMPS PS method.The comparison (see Appendix B for detail) between STAMPS result and ours showed that few PTs (see Figures B1-B3) have been selected by STAMPS due to long time decorrelation.Although the two results (see Figures B2 and B3) have a similar deformation trend in slight deformation areas, our method has much denser measurements in both significant deformation and stable urban areas, because the PT-based SBAS-InSAR method utilizes multi-reference to overcome the decorrelation effect of single reference interferometric phases, whereas the STAMPS could only select coherent PTs over a long time.As a consequence, such difference is expected.The results show that our method has great potential in mapping coastal deformation.However, our PT-based SBAS-InSAR technique also has its limitations.First, the threshold for PT selection needs to be carefully determined.Second, for a reliable InSAR measurement, a network of dense SBAS interferometric pairs needs to be constructed.Besides, the phase unwrapping error needs to be checked carefully, because such an error would become significant on the data margin for phase unwrapping on the network of irregular PTs.
It is worth remarking that not all coastal zones present subsidence signals in the ENVISAT InSAR measurements.It can be explained by temporal decorrelation [25,47], causing few usable PTs to be available.For example, there are several parts in our study area which tend to decorrelate due to a construction process.Compared with the long revisit cycle (35 days) of ENVISAT, the new generation of SAR platforms have shorter revisits, such as the ALOS-2 (14 days) and Sentinel-1A (12 days).SAR images from these new platforms generally maintain much higher coherence.In addition, only vertical displacement is retrieved, based on the assumption of non-horizontal displacement in this study.However the horizontal displacement cannot be neglected, especially for areas with a large subsidence rate.So combined multiple-pass InSAR datasets with different viewing angles can be considered to reveal ground 3-D deformations further.

Conclusions
To meet the demand of urbanization and economic development, about 80 km 2 of land has been reclaimed from the sea in Shenzhen over the past three decades.The reclaimed land undergoes local heterogeneous deformation that affects the stability of both ground constructions and underground man-made facilities.In order to monitor the reclaimed land deformation, we proposed a PTs-based SBAS-InSAR method, which has the ability to detect large amplitude deformation and choose high PTs density in a low coherence area.The spatial and temporal variability of the deformation was discussed over the lands reclaimed over the past decades in Shenzhen.Inter-comparison of the InSAR measurements from different orbits showed that the correlation between these measurements is 0.92, and the RMSE of the difference between these measurements is 1.8 mm/year, implying a good self-consistency of the InSAR measurements.The InSAR measurements showed that the reclaimed land underwent remarkable deformation (up to 25 mm/year in the LOS direction), especially at the Shenzhen Airport, Bao'an Center, Qianhai Bay, and Shenzhen Bay, and the deformation is expected to continue in the near future.Therefore, continuous deformation monitoring is necessary in these reclaimed area.

Appendix B
To evaluate the effectiveness of the proposed strategy, the dataset from track 175 is processed with the STAMPS PSI method [21,49] (see Figure B1).Numerous PTs in built-up areas have been detected by the STAMPS method and the proposed method.However, in low-coherence areas, the STAMPS method can select few PTs due to the decorrelation effect, while our strategy can select an appropriate number of PTs.The density of PTs selected by STAMPS is 301 points/km 2 , however, it is 779 points/km 2 by our method.Such difference is particularly evident in reclaimed land areas, such as Bao'an Center, the north to the Shenzhen Bay Checkpoint, and the dam of Shenzhen Airport.
For more detailed comparison, the regions A and B marked with a black rectangle (see Figure B1) were further analyzed with zoom-in view shown in Figures B2 and B3.The deformation rate maps obtained by STAMPS and our method have a very similar deformation trend (see Figures B2 and B3).However, our deformation rate can detect much larger deformation than that of the STAMPS results.As shown in Figures B2 and B3, few PTs have been selected from areas with large deformations by the STAMPS method.As a consequence, the area that undergoes rapid deformation has few measurements on STAMPS results, while appropriate PTs have been identified by our method for deformation mapping.appropriate number of PTs.The density of PTs selected by STAMPS is 301 points/km 2 , however, it is 779 points/km 2 by our method.Such difference is particularly evident in reclaimed land areas, such as Bao'an Center, the north to the Shenzhen Bay Checkpoint, and the dam of Shenzhen Airport.
For more detailed comparison, the regions A and B marked with a black rectangle (see Figure B1) were further analyzed with zoom-in view shown in Figures B2 and B3.The deformation rate maps obtained by STAMPS and our method have a very similar deformation trend (see Figures B2  and B3).However, our deformation rate can detect much larger deformation than that of the STAMPS results.As shown in Figures B2 and B3, few PTs have been selected from areas with large deformations by the STAMPS method.As a consequence, the area that undergoes rapid deformation has few measurements on STAMPS results, while appropriate PTs have been identified by our method for deformation mapping.However, the density of PTs selected by our method is larger than that of STAMPS, especially in areas with large deformation.However, in (b) the significant deformation could not be well mapped, due to the sparse PTs density of the STAMPS method.However, the density of PTs selected by our method is larger than that of STAMPS, especially in areas with large deformation.However, in (b) the significant deformation could not be well mapped, due to the sparse PTs density of the STAMPS method.

Figure 1 .
Figure 1.Study area and SAR data coverage.The study area is outlined by the black rectangle.Blue rectangles represent the coverage of the SAR dataset used.Inset, the approximate location of Shenzhen in China.The right sub-figure is the distribution of the sediment layer in Shenzhen.The colors indicate the thickness.

Figure 1 .
Figure 1.Study area and SAR data coverage.The study area is outlined by the black rectangle.Blue rectangles represent the coverage of the SAR dataset used.Red lines denote the fault lines surrounding [41,42].Inset, the approximate location of Shenzhen in China.The right sub-figure is the distribution of the sediment layer in Shenzhen.The colors indicate the thickness.

Figure 2 .
Figure 2. Land reclamation and its evolution in Shenzhen.The reclaimed land for every period is spatially overlapped on the Landsat-8 background.Inset shows the temporal evolution.The coastlines are extracted from a set of co-registered Landsat images shown in Figure A1.Firstly, coastlines are automatically extracted by image segmentation with an appropriate threshold.Then the initial coastlines are overlaid to the Landsat image in ArcGIS software-we edited the inaccurate coastline manually.

Figure 3 .
Figure 3. SBAS networks of the datasets used.(a) Track 025; (b) track 175.Each black dot denotes one image.The line between each two dots represents one interferometric pair, whose color indicates the average coherence of the corresponding interferogram.

Figure 2 .
Figure 2. Land reclamation and its evolution in Shenzhen.The reclaimed land for every period is spatially overlapped on the Landsat-8 background.Inset shows the temporal evolution.The coastlines are extracted from a set of co-registered Landsat images shown in Figure A1.Firstly, coastlines are automatically extracted by image segmentation with an appropriate threshold.Then the initial coastlines are overlaid to the Landsat image in ArcGIS software-we edited the inaccurate coastline manually.

Figure 2 .
Figure 2. Land reclamation and its evolution in Shenzhen.The reclaimed land for every period is spatially overlapped on the Landsat-8 background.Inset shows the temporal evolution.The coastlines are extracted from a set of co-registered Landsat images shown in Figure A1.Firstly, coastlines are automatically extracted by image segmentation with an appropriate threshold.Then the initial coastlines are overlaid to the Landsat image in ArcGIS software-we edited the inaccurate coastline manually.

Figure 3 .
Figure 3. SBAS networks of the datasets used.(a) Track 025; (b) track 175.Each black dot denotes one image.The line between each two dots represents one interferometric pair, whose color indicates the average coherence of the corresponding interferogram.

Figure 3 .
Figure 3. SBAS networks of the datasets used.(a) Track 025; (b) track 175.Each black dot denotes one image.The line between each two dots represents one interferometric pair, whose color indicates the average coherence of the corresponding interferogram.

Figure 5 .
Figure 5. Deformation rates generated with track 025 dataset for the whole study area.The black star locates the reference point, and the inset shows a more detail location of it.Black dashed lines represent the Shenzhen Metro Lines.The yellow line denotes the coastline in 1979.The background is a Landsat-8 image acquired in 2015.Regions A, B, C, and D marked with black rectangles will be further analyzed together with land reclamation in the discussion section.Points marked with white cross (×) are selected to show time series deformation.

Figure 6 .
Figure 6.Deformation rates generated with track 175 dataset for the whole study area.Other remakes are the same as those in Figure 5.

Figure 5 . 22 Figure 5 .
Figure 5. Deformation rates generated with track 025 dataset for the whole study area.The black star locates the reference point, and the inset shows a more detail location of it.Black dashed lines represent the Shenzhen Metro Lines.The yellow line denotes the coastline in 1979.The background is a Landsat-8 image acquired in 2015.Regions A, B, C, and D marked with black rectangles will be further analyzed together with land reclamation in the discussion section.Points marked with white cross (×) are selected to show time series deformation.

Figure 6 .
Figure 6.Deformation rates generated with track 175 dataset for the whole study area.Other remakes are the same as those in Figure 5.Figure 6. Deformation rates generated with track 175 dataset for the whole study area.Other remakes are the same as those in Figure 5.

Figure 6 .
Figure 6.Deformation rates generated with track 175 dataset for the whole study area.Other remakes are the same as those in Figure 5.Figure 6. Deformation rates generated with track 175 dataset for the whole study area.Other remakes are the same as those in Figure 5.

Figure 7 .
Figure 7. Field work pictures showing non-uniform subsidence.(a) Site of a road; (b) site of a building.

Figure 7 .
Figure 7. Field work pictures showing non-uniform subsidence.(a) Site of a road; (b) site of a building.

Figure 8 .
Figure 8.Time series of the deformation at selected points.(a) TP1; (b) TP2; (c) TP3; (d) TP4; (e) TP5; (f) TP6.Purple dots and blue triangles represent the ascending (025) and descending (175) tracks.Purple lines and blue lines represent the linear deformations.The corresponding slopes are given in the lower left corner

Figure 8 .
Figure 8.Time series of the deformation at selected points.(a) TP1; (b) TP2; (c) TP3; (d) TP4; (e) TP5; (f) TP6.Purple dots and blue triangles represent the ascending (025) and descending (175) tracks.Purple lines and blue lines represent the linear deformations.The corresponding slopes are given in the lower left corner.

Figure 9 .
Figure 9. Inter-comparison of the deformation rates between track 025 and track 175.Blue dots denote the measured points on track 025 and track 175.The purple line is the fitted linear model, i.e., f(x) = 0.9245x − 0.07292.

Figure 9 .
Figure 9. Inter-comparison of the deformation rates between track 025 and track 175.Blue dots denote the measured points on track 025 and track 175.The purple line is the fitted linear model, i.e., f(x) = 0.9245x − 0.07292.

Figure 10 .
Figure 10.Subsidence rate of Region A of Figures 5 and 6, Shenzhen Airport.(a) Track 025; (b) track 175.Black dash line with black triangles represents the Shenzhen Metro Line, and the black triangle denotes the Metro station.The yellow line denotes the coastline in 1979.Background is a Landsat-8 image acquired in 2015.

Figure 11 .
Figure 11.Subsidence rate of Region B of Figures 5 and 6, Bao'an center.(a) Track 025; (b) track 175.Other remakes are the same as in Figure10.

Figure 10 .
Figure 10.Subsidence rate of Region A of Figures 5 and 6, Shenzhen Airport.(a) Track 025; (b) track 175.Black dash line with black triangles represents the Shenzhen Metro Line, and the black triangle denotes the Metro station.The yellow line denotes the coastline in 1979.Background is a Landsat-8 image acquired in 2015.

Figure 10 .
Figure 10.Subsidence rate of Region A of Figures 5 and 6, Shenzhen Airport.(a) Track 025; (b) track 175.Black dash line with black triangles represents the Shenzhen Metro Line, and the black triangle denotes the Metro station.The yellow line denotes the coastline in 1979.Background is a Landsat-8 image acquired in 2015.

Figure 11 .
Figure 11.Subsidence rate of Region B of Figures 5 and 6, Bao'an center.(a) Track 025; (b) track 175.Other remakes are the same as in Figure 10.Figure 11.Subsidence rate of Region B of Figures 5 and 6, Bao'an center.(a) Track 025; (b) track 175.Other remakes are the same as in Figure 10.

Figure 11 .
Figure 11.Subsidence rate of Region B of Figures 5 and 6, Bao'an center.(a) Track 025; (b) track 175.Other remakes are the same as in Figure 10.Figure 11.Subsidence rate of Region B of Figures 5 and 6, Bao'an center.(a) Track 025; (b) track 175.Other remakes are the same as in Figure 10.

Figure 12 .
Figure 12.Subsidence rate of Region C of Figures 5 and 6, Qianhai Bay.(a) Track 025; (b) track 175.Other remakes are the same as in Figure10.

Figure 13 .
Figure 13.Subsidence rate of Region D of Figures 5 and 6, Shenzhen Bay.(a) Track 025; (b) track 175.Other remakes are the same as in Figure10.

Figure 13 .
Figure 13.Subsidence rate of Region D of Figures 5 and 6, Shenzhen Bay.(a) Track 025; (b) track 175.Other remakes are the same as in Figure10.

Figure 13 .
Figure 13.Subsidence rate of Region D of Figures 5 and 6, Shenzhen Bay.(a) Track 025; (b) track 175.Other remakes are the same as in Figure10.

Figure B1 .
Figure B1.Deformation rate maps extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.We can see that more PTs have been selected by our method (PT-based SBAS-InSAR).Regions A and B are further analyzed.

Figure B1 .
Figure B1.Deformation rate maps extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.We can see that more PTs have been selected by our method (PT-based SBAS-InSAR).Regions A and B are further analyzed.

Figure B1 .
Figure B1.Deformation rate maps extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.We can see that more PTs have been selected by our method (PT-based SBAS-InSAR).Regions A and B are further analyzed.

Figure B2 .
Figure B2.Deformation rate maps for region A, Bao'an Center extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.Numerous PTs have been selected in the urban area, upper-right part of the Figure.However, the density of PTs selected by our method is larger than that of STAMPS, especially in areas with large deformation.However, in (b) the significant deformation could not be well mapped, due to the sparse PTs density of the STAMPS method.

Figure B2 .
Figure B2.Deformation rate maps for region A, Bao'an Center extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.Numerous PTs have been selected in the urban area, upper-right part of the Figure.However, the density of PTs selected by our method is larger than that of STAMPS, especially in areas with large deformation.However, in (b) the significant deformation could not be well mapped, due to the sparse PTs density of the STAMPS method.

Figure B3 .
Figure B3.Deformation rate maps for region B, Shenzhen Bay extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.Numerous PTs have been selected in the urban area.However, our method can select more PTs than the STAMPS method in areas with a large deformation rate, especially the northern part of the Shenzhen Bay Check Point, which enables us to map the significant deformation.

Figure B3 .
Figure B3.Deformation rate maps for region B, Shenzhen Bay extracted by (a) PT-based SBAS-InSAR and (b) STAMPS.Numerous PTs have been selected in the urban area.However, our method can select more PTs than the STAMPS method in areas with a large deformation rate, especially the northern part of the Shenzhen Bay Check Point, which enables us to map the significant deformation.

Table 1 .
Summary information of the networks of SBAS interferometric pairs.

Table 1 .
Summary information of the networks of SBAS interferometric pairs.

Table 2 .
Estimation of the elapsed time of land surface falling sea level regarding coastal subsidence in selected regions.

Sea Level Rise v 1 (mm/year) Subsidence Rate 2 v 2 (mm/year) Elapsed Time 3 Estimated T (year)
1 average elevation of region; 2 average subsidence rate of region;3

Table A4 .
Statistics of the phase unwrapping error in track 175.