Subsidence detection by TerraSAR-X interferometry on a network of natural persistent scatterers and artificial corner reflectors
Introduction
As a geodetic technique, spaceborne interferometric synthetic aperture radar (InSAR) (Rosen et al., 2000), especially persistent scatterer interferometry (PSI) (Ferretti et al., 2000, Ferretti et al., 2001, Liu et al., 2009, Hooper et al., 2012), has exhibited great potential in mapping subsidence with advantage in spatial resolution and coverage (Buckley, 2000, Strozzi et al., 2003, Helmut, 2009, Ge et al., 2010, Cigna et al., 2012). The recent radar sensor evolution has resulted in a remarkable improvement of spatial resolution, e.g., the typical X-band radar sensor onboard the German satellite TerraSAR-X (TSX) with capability of providing SAR images at a resolution of about 1–3 m (Pitz and Miller, 2010). This extends the data availability for subsidence detection by PSI (Wegmüller et al., 2010, Liu et al., 2011).
Although the PSI technique has been widely used to monitor land subsidence in various areas, it is a challenging task to detect subsidence in vegetated areas and cultivated lands where the random surface changes occur frequently (Xia et al., 2004, Marinkovic et al., 2004, Ferretti et al., 2007, Crosetto et al., 2010) and few “natural PSs (NPSs)” are available. The NPSs here mean the existing hard objects (e.g., rocks, pylons, buildings and bridges) with temporal stability of radar backscattering (Ferretti et al., 2001). Previous investigations indicate that the man-made corner reflectors (CRs) can be deployed as the artificial PSs in the difficult areas to increase the density of PSs (Xia et al., 2004, Marinkovic et al., 2004). With availability of high-resolution SAR images, it is desired to further investigate the combined use of NPSs and CRs for subsidence detection over the difficult areas. This is vital for tracking subsidence troughs by PSI in a huge area with vegetation coverage and cultivated lands. For example, severe land subsidence should be monitored around the North China Plain due to excessive exploitation of groundwater (Liu et al., 2001). It is reported that most of the subsidence troughs in the North China Plain are being expanded in the cultivated lands (Hu et al., 2004).
To extend the application of PSI in detecting subsidence in areas with frequent surface changes, this paper presents a method of TSX PSI with a network of NPSs and CRs deployed on site. The primary procedures include detection of NPSs, identification of CRs, NPS-CR networking, phase modeling and subsidence estimation. For validation purpose, we select a suburban area of southwest Tianjin (China) in the North China Plain as the testing site, and utilize 13 TSX images collected over this area between April of 2009 and July of 2010 to extract subsidence by the method proposed. 16 CRs and 10 leveling points (LPs) are deployed in the study area. For comparison analysis, we set two types of CRs with different size around the fishponds and crop parcels. 6 CRs are the conventional ones, i.e., fixed CRs (FCRs), while 10 CRs are the ones newly designed by our research group, i.e., so-called portable CRs (PCRs) with capability of repeatable installation. The performance of the CRs is evaluated by statistical analysis of radar signatures, and the quality of subsidence measurements derived by the TSX PSI method is assessed by comparing with the leveling data.
Section snippets
Study area and data source
Fig. 1 shows the study area and the location map around Tianjin, China. A suburban area of southwest Tianjin is selected as the testing site (Xiqing district of Tianjin, marked by a smaller rectangle). The entire coverage of the TSX scene used for subsidence detection is marked by a larger rectangle in Fig. 1. Tianjin is located at the eastern coastal region of the North China Plain, bordering Beijing, the Yanshan Mountains and the Bohai Bay to the northwest, the north and the east,
Identification of NPSs and CRs
Basically, we detect NPSs by following the method proposed by Ferretti et al. (2001). The NPSs are identified on a pixel-by-pixel basis. A pixel corresponds to a NPS if its amplitude dispersion index (ADI) is smaller than 0.25 (Ferretti et al., 2001). Suppose that all the CRs are temporally coherent during the period of 13 TSX acquisitions, they should be screened out as the artificial PSs together with the NPSs. However, the location of the CRs should be precisely determined to support further
Distribution of PSs
We determined the NPSs and CRs using the procedures described in Section 3.1. Fig. 7 shows the distribution of 696779 PSs (including 16 CRs) identified from the entire study area. 16 CRs, 10 LPs and 6 NPSs selected are annotated in Fig. 7 for subsequent analysis. The averaged density of PSs in the study area is about 6968/km2. Very dense PSs appear in the built-up areas due to the availability of many man-made and natural hard objects. Seven optical-image patches corresponding to the industrial
Conclusion
We present a method of TSX PSI on a network of NPSs and CRs to extend the application of PSI in detecting subsidence in areas with frequent surface changes. The suburban area of southwest Tianjin (China) is selected as the testing site in which 16 CRs and 10 LPs are deployed. The 13 TSX images collected over this area between 2009 and 2010 are used to extract subsidence by the method proposed. Both the instrumentation and precise identification of FCRs and PCRs are presented and discussed. The
Acknowledgments
This work was jointly supported by the National Natural Science Foundation of China under Grant 41074005, the R&D Program of Railway Ministry under Grant 2008G031-5, and the 2013 Doctoral Innovation Funds of Southwest Jiaotong University. We thank Infoterra GmbH and USGS for providing TerraSAR-X SAR images and SRTM DEM, respectively.
References (28)
- et al.
Monitoring land subsidence and its induced geological hazard with Synthetic Aperture Radar Interferometry: a case study in Morelia, Mexico
Remote Sensing of Environment
(2012) - et al.
Recent advances in SAR interferometry time series analysis for measuring crustal deformation
Tectonophysics
(2012) - et al.
Review on current status and challenging issues of land subsidence in China
Engineering Geology
(2004) - et al.
A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms
IEEE Transactions on Geoscience and Remote Sensing
(2002) Radar interferometry measurement of land subsidence
(2000)- et al.
Persistent scatterer interferometry: potential, limits and initial C- and X-band comparison
Photogrammetric Engineering and Remote Sensing
(2010) Location of pixels in spaceborne SAR imagery
IEEE Transactions on Geoscience and Remote Sensing
(1982)- et al.
Linking reference target properties to its perceived RCS in SAR images
Proceedings of 9th European Conference on Synthetic Aperture Radar, Nuremberg, Germany
(2012) - Enviro-Library - Challenges and prospects of sustainable water management in Tianjin, 2008. Sustainable Groundwater...
- et al.
Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry
IEEE Transactions on Geoscience and Remote Sensing
(2000)
Permanent scatterers in SAR interferometry
IEEE Transactions on Geoscience and Remote Sensing
Submillimeter accuracy of InSAR time series: experimental validation
IEEE Transactions on Geoscience and Remote Sensing
Impact of ground subsidence on the Beijing–Tianjin high-speed railway as mapped by radar interferometry
Annals of GIS
Advances in interferometric synthetic aperture radar (InSAR) in earth system science
Progress in Physical Geography
Cited by (24)
Review of satellite radar interferometry for subsidence analysis
2022, Earth-Science ReviewsCitation Excerpt :Zhang et al. (2016b) reported a complete tracking of ground subsidence for the Beijing-Tianjin-Hebei region over a long-time span from 1992 to 2014 using the existing C-band archives (ERS, Envisat and RADARSAT-2). Having suffered dramatic ground subsidence during the last two to three decades, the spatial-temporal characteristics of the development of ground subsidence in Beijing and Tianjin have been analysed by many researchers using multi-temporal InSAR technology (Du et al., 2018a; Guo et al., 2019), also using high-resolution X-band data to increase the detail of surface evolution of strategic infrastructures (Gao et al., 2019) and to extend the application of PSI in detecting subsidence in areas with frequent surface changes (Yu et al., 2013) with the deployment of artificial corner reflectors. Supported by accessible satellite datasets and by the flourishing of software packages, subsidence monitoring is feasible worldwide (Cian et al., 2019) covering a wide range of scenarios: deltaic environment (Parcharidis et al., 2013), tunnel excavation advancement (Roccheggiani et al., 2019; Strozzi et al., 2017), land reclamation (Sun et al., 2018), peat settlement (Marshall et al., 2018), natural gas (Gee et al., 2016) and oil extraction (Moghaddam et al., 2013) and civil infrastructures analysis (Cigna et al., 2017; Corsetti et al., 2018; Luo et al., 2017).
Coal fire identification and state assessment by integrating multitemporal thermal infrared and InSAR remote sensing data: A case study of Midong District, Urumqi, China
2022, ISPRS Journal of Photogrammetry and Remote SensingCitation Excerpt :Subsidence in coal fields where mining has ceased and coal fires are ongoing is mainly caused by coal fires. As a new type of active remote sensing technology, MT-InSAR can effectively detect surface subsidence and provide new perspectives for identifying coal fire areas (Ferretti et al. 2000, 2001; Hooper 2008; Hooper et al. 2012; Yu et al. 2013). To date, scholars have applied the MT-InSAR method in conjunction with GPS and field survey data and comprehensively verified that MT-InSAR can play a valuable role in identifying coal fires (Jiang et al. 2011a; Kumar et al. 2020; Liu et al. 2019).
Landslide monitoring for risk mitigation by using corner reflector and satellite SAR interferometry: The large landslide of Carlantino (Italy)
2017, CatenaCitation Excerpt :The SI values are listed in the legend of Fig. 10. For stable CRs, the SI values are comparable with those obtained in other studies using similar CR size and SAR data (e.g. Yu et al., 2013). CR-6 shows clearly a very poor stability index, due to the mentioned vandalism exposure, which induced strong tilting of the CR structure, thus destroying the CR alignment with the sensor LOS.
Persistent Scatterer Interferometry: A review
2016, ISPRS Journal of Photogrammetry and Remote SensingCitation Excerpt :CRs have been used in different studies to assess the performance of DInSAR and PSI (Ferretti et al., 2007; Marinkovic et al., 2007; Quin and Loreaux, 2013). They have also been employed in several PSI applications, e.g. see van Leijen and Hanssen (2007), Wegmüller et al. (2010), Daqing et al. (2011), Teatini et al. (2012) and Yu et al. (2013). The deformation time series represent the most advanced PSI product.
Research on deformation monitoring of reservoir slope in high-vegetation area based on CR-InSAR technology
2024, Proceedings of SPIE - The International Society for Optical EngineeringDisentangling Shallow Subsidence Sources by Data Assimilation in a Reclaimed Urbanized Coastal Plain, South Flevoland Polder, the Netherlands
2023, Journal of Geophysical Research: Earth Surface