Next Article in Journal
Remote Sensing-Based Biomass Estimation and Its Spatio-Temporal Variations in Temperate Grassland, Northern China
Next Article in Special Issue
Ensemble Empirical Mode Decomposition Parameters Optimization for Spectral Distance Measurement in Hyperspectral Remote Sensing Data
Previous Article in Journal
A Nine-Year Climatology of Arctic Sea Ice Lead Orientation and Frequency from AMSR-E
Previous Article in Special Issue
Validation of the Two Standard MODIS Satellite Burned-Area Products and an Empirically-Derived Merged Product in South Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of InSAR and TomoSAR for Monitoring Deformations Caused by Mining in a Mountainous Area with High Resolution Satellite-Based SAR

1
Institute of Geodesy and Photogrammetry, TU Braunschweig, Braunschweig D-38106, Germany
2
Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China
3
School of Environment Science & Spatial Informatics, China University of Mining and Technology, Xuzhou 221008, China
4
Nottingham Geospatial Institute, University of Nottingham, Nottingham NG7 2TU, UK
*
Authors to whom correspondence should be addressed.
Remote Sens. 2014, 6(2), 1476-1495; https://doi.org/10.3390/rs6021476
Submission received: 20 December 2013 / Revised: 8 February 2014 / Accepted: 10 February 2014 / Published: 19 February 2014
(This article belongs to the Special Issue Analysis of Remote Sensing Image Data)

Abstract

:
Interferometric Synthetic Aperture Radar (InSAR) and Differential Interferometric Synthetic Aperture Radar (DInSAR) have shown numerous applications for subsidence monitoring. In the past 10 years, the Persistent Scatterer InSAR (PSI) and Small BAseline Subset (SBAS) approaches were developed to overcome the problem of decorrelation and atmospheric effects, which are common in interferograms. However, DInSAR or PSI applications in rural areas, especially in mountainous regions, can be extremely challenging. In this study we have employed a combined technique, i.e., SBAS-DInSAR, to a mountainous area that is severely affected by mining activities. In addition, L-band (ALOS) and C-band (ENVISAT) data sets, 21 TerraSAR-X images provided by German Aerospace Center (DLR) with a high resolution have been used. In order to evaluate the ability of TerraSAR-X for mining monitoring, we present a case study of TerraSAR-X SAR images for Subsidence Hazard Boundary (SHB) extraction. The resulting data analysis gives an initial evaluation of InSAR applications within a mountainous region where fast movements and big phase gradients are common. Moreover, the experiment of four-dimension (4-D) Tomography SAR (TomoSAR) for structure monitoring inside the mining area indicates a potential near all-wave monitoring, which is an extension of conventional InSAR.

1. Introduction

Synthetic aperture radar (SAR) is a coherent imaging method by active microwave with a long aperture [1]. In other words, SAR is a modification of a basic imaging radar system that creates an artificially long antenna synthetically from a small moving antenna to improve spatial resolution in the azimuth (along-track) direction. Interferometric SAR (InSAR) measures the phase differences between two SAR images, which are primarily related to the surface topographic components and deformation. Differential interferometric synthetic aperture radar (DInSAR) is a technique based on InSAR. It removes the topographic phase with an external digital elevation model (DEM) or another interferogram, which leaves the phase component related to the land motion between the two images. DInSAR is an appropriate method for fast and large deformations [2]. In addition, in order to detect slow deformation over a long time span, some other InSAR techniques have been developed. The first type of method, which we will refer to here as Persistent Scatterers Interferometry (PSI), performs an analysis of ‘persistent scatterer’ pixels on a succession of time-ordered images. Another approach, the SBAS (Small BAseline Subset) method, forms many interferograms with short baselines and derives movements with time [3,4].
The advantage of the use of InSAR techniques, such as DInSAR or PSI/SBAS methods, is that the deformation of many thousands of points over a scene can be determined very quickly and at very low cost when compared to the effort needed to perform a ground survey of points at the same density. Fruitful results with regard to wide applications have been reported by InSAR researchers in the past 15 years [57]. However, in this study, both DInSAR and PSI/SBAS methods revealed their limitations when applied to a mountainous mining area with fast earth surface movements [68]. Based on the requirements from mining engineering and the different features of various InSAR techniques, a combined technique SBAS-DInSAR is applied to this mountainous area that is severely affected by mining activities.

2. Application to Coal Mining

The extraction of underground coal deposits through mining always results in damages to geological structures. Mountainous areas are susceptible to not only surface subsidence, but, also, landslides and mud-rock flow. From a safety point-of-view, it is therefore necessary to monitor surface movements during and after the coal mining activities, in order to detect and warn of potential and actual geological hazards.
Normally, subsidence at coal mines occur at large scales, on the order of tens of square kilometers [8]. Therefore, InSAR techniques are appropriate for the monitoring of the surface effects of coal mining. Some excellent results from the application of InSAR for operating coal mines have been shown [9]. In [9], 20 Japanese Earth Resources Satellite (JERS) SAR scenes and 37 Environmental Satellite Advanced Synthetic Aperture Radar (ENVISAT ASAR) scenes covered Fengfeng area were analyzed, which is 560 square kilometers in size and has 14 coal mines. DInSAR results were shown to be comparable to leveling and historical coal mining data. On the other hand, one test of monitoring abandoned mines using PSI was carried out in Nottingham [10]. In that research, the Nottinghamshire Corner Reflector Array (NCRA) was established. Data processing showed the rates for both PSI and Global Positioning System (GPS) derived results are comparable in all cases. In addition, the results showed that a corner reflector network may provide some possibilities for calibration and improvements to a PSI approach to subsidence monitoring.

3. Xishan Coal Mine

3.1. Xishan Coal Mine Geography

Xishan coal mine is located in Shanxi Province, China (Figure 1). There are, in total, eight mines in this area, with a total coverage of 304.8 square kilometers. Five of them are of interests to us in this study due to their strong excavations. Three-quarters of this area feature large mountains with steep slopes and ravines, which are mainly composed of Carboniferous, Permian sandstone, shale, and Quaternary loess. Landslides, debris flows, and avalanches are quite common in this area [11].
As the relative altitudes among these mountains are 500 m to 800 m, SAR images of such an area will be greatly distorted by the oblique imaging geometry of side-looking radar, often causing layover and shadow effects [12]. Layover and shadow areas are anomalous areas in a SAR image containing no useful information and should not be considered in any precise analysis. Furthermore, the presence of extensive vegetation and very few man-made structures (which usually provide strong scatterers) make mountainous regions a very challenging environment to apply InSAR techniques. In this particular location, repeated mining in Xishan is being carrying out in the lower coal seams because the upper seams will be exhausted soon. It is predicted that such repeated mining activities will induce bigger and faster deformation than before. As a result, decorrelation may be a big problem in data processing.

3.2. Topographic Analysis

Rather than trees, grasses and bushes are common in this area, and mountains are of large relative altitudes (Figure 2). They may affect the signal of TerraSAR-X as its wavelength is relative short. Based on some analysis using the ENVISAT images, it is found that there were some layover effects, but the shadow affected almost a negligible number of pixels when such a small incidence angle (23°) was used in SAR observation.

4. SAR Data and Software Packages

The Japan Aerospace Exploration Agency (JAXA) launched the Advanced Land Observing Satellite (ALOS) satellite in 2006 which contained a Phased-Array L-band SAR (PALSAR). The characteristics of Fine Beam Single Polarisation (FBS) mode of PALSAR are summarized in Table 1, as FBS is the most common for interferometric applications. Four Level 1.1 Single Look Complex (SLC) ALOS images, provided in CEOS data format were collected from 2008 to 2010. Data specifications are shown in Table 1.
The baseline characteristics of the four images are shown in Table 2, the spatial baseline being the distance between the satellite and master positions and the temporal baseline being the time separation between the two.
Unlike other sensors, ALOS L1.1 products are not zero-Doppler. The consequence is that the geometric model used during the processing is slightly different and varies from product-to-product. Although the variation can be negligible [12], in general this is not the case and, thus, a slightly modified processing chain has to be applied. An improved algorithm in software, PUNNET, was developed to cope with this [12].
In addition, eight ENVISAT ASAR scenes were collected, as listed in Tables 3 and 4. The baselines are controlled so that only pairs with short baselines will be processed.
Germany launched the high resolution SAR satellite TerraSAR-X on 15 June 2007, and TanDEM-X on 21 June 2010. As for this study, 21 high-resolution TerraSAR-X scenes in strip map mode with the time span of 4 April 2012, to 21 November 2012, were programmed. Steep incidence angles were used so as to avoid strong shadow effects in this mountainous area (Table 5).
Figure 3 indicates the baseline information of 21 TerraSAR-X images used in this study. We select the image acquired on 18 May 2012, as the master image, and calculate baselines of other 20 scenes related to the master image. According to the estimation, most perpendicular baselines were controlled between −100 m to 150 m, which are obviously advantageous for deformation monitoring.
Apart from the modified method SBAS-DInSAR, we also applied conventional DInSAR and PSI/SBAS techniques by using various software packages for the data processing in this paper:
  • A shell based on Delft Object-oriented Radar Interferometric Software (DORIS), so-called Automated DORIS Environment (ADORE) was used to generate inteferograms [13,14];
  • The PUNNET software was used for some topographic analysis. In PUNNET, SBAS is based on the method used by and implemented by Andrew Sowter at the University of Nottingham and includes some additional modules for geometric analyses [3,4,12];
  • StaMPS package [15]; and
  • GAMMA processed results were provided by CUMT.

5. Analysis of Interferometric Results

5.1. Analysis of ENVISAT and ALOS Archive Data

In the following Figure 4, DInSAR analysis is applied to archive data, which includes ENVISAT and ALOS scenes. We processed all possible pairs and then abandoned inteferograms with low quality due to temporal and/or spatial decorrelations, or with large baseline errors. The analysis of individual interferograms illustrated the influence of topography on the interferometric phase. Large scale relief maps were collected and used to generate high resolution digital elevation model (DEM), which was then used during DInSAR processing to compensate topography phase, and aid the interpretation and unwrapping of the differential interferograms [16]. An adaptive filtering technique was applied to differential interferograms to clarify the fringes and to reduce noise [17]. However, filters must always be used with care and we should keep in mind that filtering can remove more subtle effects. For phase unwrapping, the Statistical-cost, Network-flow Algorithm for PHase Unwrapping (SNAPHU) approach was applied [18].
Considering similar time spans, subsidence information derived from ALOS data using pair 25 December 2009–9 February 2010 and from ENVISAT data using pair 4 October 2009–21 February 2010, could be compared. In this study, the ALOS data are ascending, while ENVISAT are descending. Therefore, in spite of detecting the same displacement, view angles are different from two images. However, L band detected decimeter-level but ENVISAT showed only centimeter-level movements. The reason might be the large phase gradients induced by strong mining activities at this area, when the maximum detectable displacement in radian is π between two pixels. One solution to such problem may be increasing the wavelength and/or reducing the revisit time. The other solution could be using higher resolution SAR data where smaller pixels are available.

5.2. Analysis of TerraSAR-X Data

In the archive data analysis, we have obtained a general impression of the rapid and large deformation in Xishan. Due to the large phase gradients, deformation signals could not be derived correctly by phase unwrapping. According to the PSI results in Figure 5 produced by StaMPS, persistent scatterers candidates are absent in many regions. It is easily explicable that many blank areas (without PS point) appeared in Figure 5, as such areas were suffering strong subsidence, and the ground surface changes severely so that no phase-stable point is selectable. In these bland areas, no information can be obtained though we know such areas may coincide with the mining activities.
In this study, we combined DInSAR and time series analysis using SBAS method, i.e., SBAS-DInSAR, to 21 scenes time-ordered images (Table 6). In such SBAS-DInSAR processing, we calculated time series based on a stack of DInSAR images, which were strictly selected according to their high qualities. Unlike normal SBAS method will do, we do not select ‘good points’ by amplitude and/or phase evaluation but only work with ‘surface’ here, as with conventional DInSAR processing. The main aim is to find out the precise locations of deformation areas, which is an essential parameter for mining engineering. The basis flow of data analysis here includes combination of interferometric pairs, interferogram generation and time series analysis. Time series analysis was carried out when 22 scenes of interferogram were ready. The first acquisition, 4 April 2012, was used as the reference (the start point of time series analysis). Therefore, the deformation phases on all other dates are related to that on 4 April 2012. A sequence of cumulative subsidence were estimated and then covered with the mining map (Section 5.3).
A control point which locates in the subsidence area was chosen for this comparative analysis. According to GPS data processing by Precise Point Positioning (PPP) approach [19], the control point moved −0.36 m in vertical direction from 4 April 2012 to 21 November 2012 (see Figure 6). The InSAR derived time series has the quite similar trend though the maximum subsidence is −0.16 m, which is smaller than GPS observed. The reason could be the maximum Detectable Deformation Gradient (DDG), while the maximum detectable displacement in radian is π between two pixels as discussed before. In other words, the phase cannot be unwrapped correctly if the actual relative displacement between two pixels is more than π. During the repeated mining progress, the ground surface was suffering from strong subsidence with large gradient. Therefore, the deformation gradient is far beyond the frontiers of InSAR’s maximum DDG.
Although the maximum subsidence is hard to derive here, the area with displacement could be detected by InSAR. That means some other important parameters in mining engineering are still able to be estimated using InSAR technique. An apparent advantage is InSAR’s object is ‘surface’ rather than discrete ‘point’ in GPS or other techniques, which is extremely necessary in this study (in Section 5.3).

5.3. Correlation Analysis of InSAR Results with Mining Data

Concurrent with underground mining activities, subsidence usually occurs synchronously due to cavities created by underground mining. The extraction of coal seam removes support from the overlying rock strata and, therefore, makes them collapse into the void space that is created by the excavation. Surface deformation varies with the type of mine, the geometry of the mine, and manner of excavations. After excavation, subsidence may continue for some time in a gradual manner. In some cases, subsidence would even stop for a period, to be followed by sudden or steady deformation at some later date.
One product, Subsidence Hazard Boundary (SHB) or subsidence hazard map, provided by geodesy engineer is requisite for mining production. Subsidence hazard boundary indicates areas where the ground surface suffered from subsidence related to underground mining activities. Two terms have to be discussed when talking about SHB: Advanced influence angle—The ground surface starts to move (10 mm as defined) before the underground working face arrives. In the sagittal plane, the angle between the horizontal line and the line connecting the working face and the ground starting moving point on the coal pillar side is, so-called, the advanced influence angle; and advanced influence distance—The horizontal distance between this starting moving point and the working face is a advanced influence distance.
As mentioned, the underground mining will result in a large subsidence area at the ground surface and then comes into being a trough. In Figure 7 (Adapted, [20]), a trough W1 will appear when the working face arrives at point 1. Along the working face carries on, the trough will expand to be W2 with a larger depth. The maximum value of subsidence will not increase any more when the working face reaches to point 3, as this is then critical extraction. Following, W4 is wider than W3 but its depth will keep the same as that of W3. In order to obtain SHB and relevant parameters, we carried out correlation analysis of mining and InSAR measurements at a working face 18207, where the mining schedule is available.
In Figure 8, there are several goafs, which are marked as 12212, 12208, 12206, and 12202 in one layer. In another layer, 18205 and 18203 are goafs. Working face 18a203 started to mining in the latter half of the excavation of working face 18207. Apparently, surrounding ground will be affected by repeated mining influence.
The extraction of coal in working face 18207 began in May 2011, and was completed in October 2012. The distance of cumulative tunneling is 870 m until 5 April 2012. During the period with SAR data acquisitions (from 5 April 2012 to 22 November 2012), the distance of cumulative tunneling is 370 m. Time ordered InSAR results are then overlaid by mining maps in ArcGIS [11] with the aim of parameter extraction (Figure 9).
In the displacement map 5 April 2012–2 July 2012, another trough appeared on the eastern side of the trough caused by the working face 18207. This is due to the started mining in the 18a203 working face. As time went on, these two subsidence troughs expanded wider and wider, and then merged with each other. If viewed from the vertical direction, 18207 and 18a203 working faces have different depths. The interconnecting area between these two troughs was suffering from combined subsidence induced by the falling and collapsing in two working faces underground. According to the mining record and the subsidence maps we derived already, advanced influence angles and distances in working face 18207 during the observation time were estimated and shown in Figure 10.
From 8 May 2012 to 28 September 2012, the advanced influence angle mainly varies from 43° to 46° when 41° on 4 August 2012 dispersed around 2° (see Figure 10). After 28 September 2012, no clear expansion of subsidence boundary was detected, although the maximum displacement increased. According to the displacement map based on pair 2012.04.05 and 2012.11.22, the final subsidence hazard boundary around this working face was shown in Figure 11.
The accuracy of SHB estimation is mainly subject on the accuracy of geocoding. In order to assess the accuracy of the geocoding, the displacement maps with the corner reflectors were transformed to the World Geodetic System 1984 (WGS84) in UTM. Moreover, GPS processed results of the corner reflectors were used to compare with and evaluate the coordinates of corner reflectors in the displacement maps. The comparison demonstrated that the average error vector was 4 m.
Comparing to subsidence maps generated by archive ENVISAT and ALOS images, TerraSAR data sets provided approximate 40 cm subsidence in half-a-year even though its wavelength is much shorter. This mainly benefits from TerraSAR-X’s short revisit cycle and the reduced pixel size by studying of DDG.

6. Monitoring by 4-D TomoSAR

Tomography SAR (TomoSAR) has been applied in urban area or construction monitoring successfully [2124]. In this mountainous mining area with both complex terrain and buildings, we wonder the performance of this new technique, and how TomoSAR can deal with, for instance, layover phenomenon.
Traditional SAR imaging process projects a three-dimensional scene into two-dimensional plane (azimuth-range plane). Due to the particular imaging geometry of SAR, the target with shorter slant range is imaged nearer the sensor than the target with longer slant range. Especially, two targets will be imaged at the same range cell when the slant ranges of the two targets are the same. This character causes the top and bottom inverted phenomenon and layover phenomenon. The phase discontinuity at the top and bottom inverted area causes the unwrapping phase process to fail, which is an important step in traditional InSAR technologies and decides the precision of the result. Furthermore, the traditional InSAR technologies are completely ineffective at the layover area, as all traditional InSAR technologies are based on the assumption that there is only one scatter in each cell. Two targets at the same cell maybe deform at different velocity, which causes the DInSAR technology out of work [21]. The top and bottom inverted phenomenon and layover phenomenon are extremely common in the urban area, as there are many high buildings [22]. The slant range of the target at the roof of the high building is shorter or equivalent to the range of the target at the floor around that building. However, the urban area is much more important in the environmental risk monitoring field.
TomoSAR is an effective method, which restrains the top and bottom inverted SAR phenomenon or the layover phenomenon without phase unwrapping process, by using the synthetic aperture along elevation direction, as shown in Figure 12 [2125]. Given a definite radar range r(0) for the master SAR image, one can find a virtual point Pref on the ellipsoidal reference surface for the real point P [25]. The slant range of the point Pref is the same as the slant range of the point P and |R0| = r(0). For point P, the direction from the sensor to the virtual point Pref is normally defined as the range direction, when the direction perpendicular to the range direction is set to be the elevation direction, shown as r and s in Figure 12, respectively. Bv and Bh are the baselines, perpendicular and parallel to the reference radar look vector R0, respectively.
The repeat-pass mode can obtain a series of SAR images of the same scene with a temporal baseline. 4-D and extended 4-D TomoSAR can estimate the subsidence and the thermal dilation. As the samplings along the elevation direction are very few and the targets in the echo cell are sparse, the Compressive Sensing (CS) method is more attractive than Backprojection, Singular value decomposition, and Wiener Regularization methods, as the CS method can get a higher resolution with a lower sample ratio. The traditional signals reconstructing approach follows the Shannon sampling theorem, which requires that the sampling rate be twice the highest frequency of the signals. The CS method provides an approach of data reconstruction which overcomes this theorem [25]. Comparing with the 4-D TomoSAR approach using CS method in [25], the proposed algorithm has one modification that the deformation velocity difference of two targets at the same cell is limited in dealing with large subsidence and serious layover phenomenon area. In this study, we chose two test sites in this mining area to see if there is any subsidence. There is no ground monitoring carried out, such as by leveling or GPS along with the mining activities. Figure 13 shows the 4-D TomoSAR results at AOI NO.1, a small town named Malan, where the layover phenomenon is not serious. There are many signal parameter scatters can be used to estimate the subsidence. Figure 13a,b show the area we analyzed. The analyzed area shown in the red rectangular is near a mine. Figure 13a is the layover map, which shows the layover effects in this area. Figure 13c,e are the height estimation results, while Figure 13d,f present the subsidence estimation results. The estimated results show that the subsidence is about 6.5 cm/year.
Figure 14 shows the 4-D TomoSAR results in AOI NO.2, the city of Gujiao, at where the layover phenomenon is serious. Both a playground and a business quarter are used to analyze the subsidence, as shown in Figure 14b. Figure 14c shows the subsidence of signal scatter and Figure 14d and Figure 14e show the subsidence of double scatters of the playground. Figure 14f plots the subsidence of signal scatter and Figure 14g and Figure 14h plot the subsidence of double scatters of the business quarter. The average deformation velocity is 8 cm/year.
Due to the serious layover phenomenon, the assumption that there is only one scatter in each cell is not met, which brings signification estimation error, the height estimation result, and deformation estimation result. However, the TomoSAR, including 3-D, 4-D, and extend 4-D TomoSAR, is not based on that assumption, and profiles the scattering distribution along the elevation direction and estimates the height and subsidence of the each scatter in one cell. Hence, layover phenomenon affects TomoSAR a little. With the subsidence analysis of Gujiao city, using 4-D TomoSAR technology, the estimated average subsidence was 7 cm/year.

7. Conclusions

Multiple platform SAR datasets were processed with DInSAR, PSI and SBAS-DInSAR techniques. The results showed the advantages of InSAR technique for large scale monitoring. However, applications to mountainous areas lead to more challenges than in other cases, which require us to consider not only spatial or/and temporal de-correlation, due to mountains terrain and baselines (spatial and temporal), but also topographic effects. In this study of Xishan coal mine monitoring, InSAR fulfils the aim of SHB derivation. Moreover, TerraSAR-X provides a new possibility for InSAR application in mountainous mining area, due to its high resolution and short revisit cycle. In addition, we also employed the extension of conventional InSAR, i.e., an extended 4-D TomoSAR, to measure the structure movements and indicated promising results.

Acknowledgments

This research work has been partly supported by NSFC (the National Natural Science Foundation of China, No.: 51174191). TerraSAR-X images were provided by DLR under the scientific proposal LAN1185. Thanks for some plotting supports provided by Yahya Ghassoun in IGP and Jianlin Zhao in Catholic University of Louvain. The authors would like to thank three anonymous reviewers for their meticulous works in enhancing the quality of this paper.

Author Contributions

Donglie Liu developed the main idea that led to this paper. Yunfeng Shao developed the TomoSAR methods and made discussions of relevant problems. Zhenguo Liu and Zhengfu Bian provided DInSAR and GPS data processing. Andrew Sowter gave support to topographic analysis. Björn Riedel and Wolfgang Niemeier contributed to the final discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bamler, R.; Hartl, P. Synthetic aperture radar interferometry. Inverse Probl 1998, 14. [Google Scholar] [CrossRef]
  2. Raucoules, D.; Colesanti, C.; Carnec, C. Use of SAR interferometry for detecting and assessing ground subsidence. Comptes Rendus Geosci 2007, 339, 289–302. [Google Scholar]
  3. Berardino, P.; Fornaro, G.; Lanari, R.; Sansosti, E. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms. IEEE Trans Geosci. Remote Sens 2002, 40, 2375–2383. [Google Scholar]
  4. Lanari, R.; Lundgren, P.; Manzo, M.; Casu, F. Satellite radar interferometry time series analysis of surface deformation for Los Angeles, California. Geophys. Res. Lett 2004, 31. [Google Scholar] [CrossRef]
  5. Lu, Z.; Rykhus, R.; Masterlark, T.; Dean, K.G. Mapping recent lava flows at Westdahl Volcano, Alaska, using radar and optical satellite imagery. Remote Sens. Environ 2004, 91, 345–353. [Google Scholar]
  6. Colesanti, C.; Mouelic, S.L.; Bennani, M.; Raucoules, D.; Carnec, C.; Ferretti, A. Detection of mining related ground instabilities using the Permanent Scatterers technique—A case study in the east of France. Int. J. Remote Sens 2005, 26, 201–207. [Google Scholar]
  7. Perski, Z.; Jura, D. ERS SAR interferometry for land subsidence detection in coal mining areas. Earth Obs. Q 1999, 63, 25–29. [Google Scholar]
  8. Raucoules, D.; Bourgine, B.; de Michele, M.; Le Cozanet, G.; Closset, L.; Bremmer, C.; Veldkamp, H.; Tragheim, D.; Bateson, L.; Crosetto, M. Persistent Scatterers Interferometry Independent Validation and Intercomparison of Results–Final Report; BRGM/RP-55649-FR; ESA: Frascati, Italy, 2007. [Google Scholar]
  9. Cao, L.; Zhang, Y. Coal mine land subsidence monitoring by using spaceborne insar data: A case study in fengfeng, hebei province, China. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci 2008, 37, 255–261. [Google Scholar]
  10. Leighton, J.; Sowter, A.; Tragheim, D.; Bingley, R.; Teferle, F. Land motion in the urban area of Nottingham observed by ENVISAT-1. Int. J. Remote Sens 2013, 34, 982–1003. [Google Scholar]
  11. Liu, Z.; Bian, Z.; Lv, F.; Dong, B. Monitoring on subsidence due to repeated excavation with DInSAR technology. Int. J. Min. Sci. Technol 2013, 23, 173–178. [Google Scholar]
  12. Sowter, A. Orthorectification and interpretation of differential InSAR data over mountainous areas: A case study of the May 2008 Wenchuan Earthquake. Int. J. Remote Sens 2010, 31, 3435–3448. [Google Scholar]
  13. Kampes, B.; Usai, S. Doris: The Delft Object-Oriented Radar Interferometric Software. Proceedings of the 2nd International Symposium on Operationalization of Remote Sensing, Enschede, The Netherlands, 16–20 August 1999.
  14. Osmanoğlu, B.; Dixon, T.H.; Wdowinski, S.; Cabral-Cano, E.; Jiang, Y. Mexico City subsidence observed with persistent scatterer InSAR. Int. J. Appl. Earth Obs. Geoinf 2011, 13, 1–12. [Google Scholar]
  15. Hooper, A.; Segall, P.; Zebker, H. Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos. J. Geophys. Res 2007, 112. [Google Scholar] [CrossRef]
  16. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L. The shuttle radar topography mission. Rev. Geophys 2007, 45. [Google Scholar] [CrossRef]
  17. Goldstein, R.M.; Werner, C.L. Radar interferogram filtering for geophysical applications. Geophys. Res. Lett 1998, 25, 4035–4038. [Google Scholar]
  18. Chen, C.W.; Zebker, H.A. Network approaches to two-dimensional phase unwrapping: Intractability and two new algorithms. JOSA A 2000, 17, 401–414. [Google Scholar]
  19. Geng, J.; Teferle, F.; Shi, C.; Meng, X.; Dodson, A.; Liu, J. Ambiguity resolution in precise point positioning with hourly data. GPS Solut 2009, 13, 263–270. [Google Scholar]
  20. He, G. Mining Subsidence; CUMT Publisher: Xuzhou, China, 1991. (In Chinese) [Google Scholar]
  21. Fornaro, G.; Serafino, F. Imaging of single and double scatterers in urban areas via SAR tomography. IEEE Trans Geosci. Remote Sens 2006, 44, 3497–3505. [Google Scholar]
  22. Fornaro, G.; Reale, D.; Verde, S. Bridge thermal dilation monitoring with millimeter sensitivity via multidimensional SAR imaging. IEEE Geosci. Romote Sens. Lett 2013, 10, 677–681. [Google Scholar]
  23. Zhu, X.X.; Bamler, R. Very high resolution spaceborne SAR tomography in urban environment. IEEE Trans. Geosci. Remote Sens 2010, 48, 4296–4308. [Google Scholar]
  24. Reale, D.; Fornaro, G.; Pauciullo, A.; Zhu, X.; Bamler, R. Tomographic imaging and monitoring of buildings with very high resolution SAR data. IEEE Geosci. Remote Sens. Lett 2011, 8, 661–665. [Google Scholar]
  25. Zhu, X.X.; Bamler, R. Tomographic SAR Inversion by L1 Norm Regularization—The compressive sensing approach. IEEE Trans. Geosci. Remote Sens 2010, 48, 3839–3846. [Google Scholar]
Figure 1. Five active coal mines are shown below in red polygon on Google Earth (UTM). Blue and yellow frame indicate coverage of ENVISAT and ALOS images, respectively. The SAR amplitude image indicates the coverage of TerraSAR-X acquisitions.
Figure 1. Five active coal mines are shown below in red polygon on Google Earth (UTM). Blue and yellow frame indicate coverage of ENVISAT and ALOS images, respectively. The SAR amplitude image indicates the coverage of TerraSAR-X acquisitions.
Remotesensing 06 01476f1
Figure 2. Geography of Xishan shown by 90-m SRTM DEM. The picture at the upper right corner gives a general view of Xishan’s geography. Tunlan mine and Malan town are the AOIs (Area Of Interest) of SBAS-DInSAR processing and TomoSAR processing that will be dealt with later, respectively. Gujiao city is also one object of TomoSAR test.
Figure 2. Geography of Xishan shown by 90-m SRTM DEM. The picture at the upper right corner gives a general view of Xishan’s geography. Tunlan mine and Malan town are the AOIs (Area Of Interest) of SBAS-DInSAR processing and TomoSAR processing that will be dealt with later, respectively. Gujiao city is also one object of TomoSAR test.
Remotesensing 06 01476f2
Figure 3. Baselines of TerraSAR-X data set.
Figure 3. Baselines of TerraSAR-X data set.
Remotesensing 06 01476f3
Figure 4. Displacement maps with clear deformation signals generated from ENVISAT and ALOS data. (a) ALOS. 2009.12.25–2010.02.09 B = 467 m; (b) ENVISAT. 2008.01.13–2008.12.28 B = −12 m; (c) ENVISAT. 2008.12.28–2009.05.17 B = 46 m; (d) ENVISAT. 2009.10.04–2010.02.21 B = −42 m. Red color indicates area with subsidence. Image centre coordinates: 37°51′27.90″N, 112°7′35.30″E.
Figure 4. Displacement maps with clear deformation signals generated from ENVISAT and ALOS data. (a) ALOS. 2009.12.25–2010.02.09 B = 467 m; (b) ENVISAT. 2008.01.13–2008.12.28 B = −12 m; (c) ENVISAT. 2008.12.28–2009.05.17 B = 46 m; (d) ENVISAT. 2009.10.04–2010.02.21 B = −42 m. Red color indicates area with subsidence. Image centre coordinates: 37°51′27.90″N, 112°7′35.30″E.
Remotesensing 06 01476f4
Figure 5. Full-scene PSI results of TerraSAR-X data. The abscissa is longitude when the ordinate is latitude.
Figure 5. Full-scene PSI results of TerraSAR-X data. The abscissa is longitude when the ordinate is latitude.
Remotesensing 06 01476f5
Figure 6. Time series of a control point displacement during 220 days derived by GPS and InSAR, respectively.
Figure 6. Time series of a control point displacement during 220 days derived by GPS and InSAR, respectively.
Remotesensing 06 01476f6
Figure 7. Surface subsidence trough (Adapted, [20]). ω is the advanced influence angle.
Figure 7. Surface subsidence trough (Adapted, [20]). ω is the advanced influence angle.
Remotesensing 06 01476f7
Figure 8. Working face 18207 and other vicinity working faces. Green points are corner reflectors. Blue lines indicate the working face while green lines are the progress line. Several abandoned mines (goaf) in other layers were shown as well. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Figure 8. Working face 18207 and other vicinity working faces. Green points are corner reflectors. Blue lines indicate the working face while green lines are the progress line. Several abandoned mines (goaf) in other layers were shown as well. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Remotesensing 06 01476f8
Figure 9. Time ordered subsidence maps around 18207 working face. Green points are corner reflectors and Green lines are progress lines. The unit of subsidence is mm. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Figure 9. Time ordered subsidence maps around 18207 working face. Green points are corner reflectors and Green lines are progress lines. The unit of subsidence is mm. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Remotesensing 06 01476f9
Figure 10. Expansion of subsidence trough caused by the 18207 working face.
Figure 10. Expansion of subsidence trough caused by the 18207 working face.
Remotesensing 06 01476f10
Figure 11. Final SHB caused by 18207 working face during the observation time. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Figure 11. Final SHB caused by 18207 working face during the observation time. Image center coordinates: 37°53′32.51″N, 112°7′16.85″E.
Remotesensing 06 01476f11
Figure 12. TomoSAR configuration. Pink circles are the positions of SAR, while black circles are points on the ground.
Figure 12. TomoSAR configuration. Pink circles are the positions of SAR, while black circles are points on the ground.
Remotesensing 06 01476f12
Figure 13. The 4-D TomoSAR result at AOI NO.1. (a) is the layover map with AOI NO.1 in the red rectangle; red on this map indicates layover. (b) shows two objects in AOI NO.1 by the red rectangular. (c) and (e) are the height estimation results, and (d) and (f) are the subsidence estimation results. Note that (a) was flipped and geocoded.
Figure 13. The 4-D TomoSAR result at AOI NO.1. (a) is the layover map with AOI NO.1 in the red rectangle; red on this map indicates layover. (b) shows two objects in AOI NO.1 by the red rectangular. (c) and (e) are the height estimation results, and (d) and (f) are the subsidence estimation results. Note that (a) was flipped and geocoded.
Remotesensing 06 01476f13aRemotesensing 06 01476f13b
Figure 14. The 4-D TomoSAR result at area with serious layover phenomenon. (a) is the layover map of AOI NO.2 defined by a red circle; red color indicates layover effect. (b) indicates two objects in AOI NO.2. (c) shows the subsidence of signal scatter, and (d) and (e) show the subsidence of double scatters of the playground. (f) shows the subsidence of signal scatter and (g) and (h) show the subsidence of double scatters of the business quarter. Note that (a) was flipped and geocoded.
Figure 14. The 4-D TomoSAR result at area with serious layover phenomenon. (a) is the layover map of AOI NO.2 defined by a red circle; red color indicates layover effect. (b) indicates two objects in AOI NO.2. (c) shows the subsidence of signal scatter, and (d) and (e) show the subsidence of double scatters of the playground. (f) shows the subsidence of signal scatter and (g) and (h) show the subsidence of double scatters of the business quarter. Note that (a) was flipped and geocoded.
Remotesensing 06 01476f14aRemotesensing 06 01476f14b
Table 1. Basic characteristics for ALOS PALSAR data in FBS mode used in this study.
Table 1. Basic characteristics for ALOS PALSAR data in FBS mode used in this study.
ALOS PALSARParameters
Frequency1,270 MHz (L-Band)
Wavelength23.6 cm
PolarisationHH
Swath Width60 km
Incidence Angle∼38°
Range Pixel Spacing4.68 m
Azimuth Pixel Spacing3.14 m
Orbit Repeat Cycle46 days
Precise Orbit Accuracy<1 m
Table 2. Specifications of ALOS data used.
Table 2. Specifications of ALOS data used.
Acquisition DateSpatial Baseline (m)Temp-Baseline (days)
2008/2/41,554.62−690
2008/12/22−2,363.30−368
2009/12/25 (master)00
2010/2/9600.2046
Table 3. Basic characteristics for ENVISAT ASAR data used in the study (Strip Map imaging).
Table 3. Basic characteristics for ENVISAT ASAR data used in the study (Strip Map imaging).
ENVISAT ASARParameters
Frequency5.3 GHz (C-Band)
Wavelength5.6 cm
PolarisationHH
Swath Width100 km
Incidence Angle∼23°
Range Pixel Spacing7.8 m
Azimuth Pixel Spacing4.05 m
Orbit Repeat Cycle35 days
Precise Orbit AccuracyCentimeters
Table 4. Specifications of ENVISAT data used.
Table 4. Specifications of ENVISAT data used.
No.Acquisition DateSpatial Baseline (m)Temp-Baseline (days)
12008/1/13167.8−630
22008/12/28159.2−280
32009/5/17204.3−140
42009/7/26127.4−70
52009/10/4 (master)00
62010/2/21−38.8140
72010/7/1162.2280
82010/8/15175.4315
Table 5. TerraSAR-X parameters used in this study.
Table 5. TerraSAR-X parameters used in this study.
TerraSAR-Xparameters
Frequency9.6 GHz (L-Band)
Wavelength3.1 cm
PolarisationHH
Swath Width50 km
Incidence Angle∼26°
Range Pixel Spacing0.9 m
Azimuth Pixel Spacing1.9 m
Orbit Repeat Cycle11 days
Precise Orbit Accuracy∼10 cm
Table 6. TerraSAR-X pairs formed in this study.
Table 6. TerraSAR-X pairs formed in this study.
PairsMasterSlaveBtemporal/dayBperp/m
12012/4/42012/4/151118
22012/4/152012/4/2611141
32012/4/152012/5/722−101
42012/4/152012/5/183310
52012/5/72012/5/1811113
62012/5/182012/5/291122
72012/5/292012/6/91112
82012/6/92012/6/2011−182
92012/6/202012/7/11178
102012/7/12012/7/121168
112012/7/122012/7/2311−155
122012/7/122012/8/43312
132012/8/32012/8/411152
142012/8/32012/8/252296
152012/8/142012/8/2511−58
162012/8/252012/9/511107
172012/8/252012/9/1622−46
182012/9/162012/9/27113
192012/9/272012/10/8117
202012/10/82012/10/191171
212012/10/192012/11/102230
222012/11/102012/11/2111−61

Share and Cite

MDPI and ACS Style

Liu, D.; Shao, Y.; Liu, Z.; Riedel, B.; Sowter, A.; Niemeier, W.; Bian, Z. Evaluation of InSAR and TomoSAR for Monitoring Deformations Caused by Mining in a Mountainous Area with High Resolution Satellite-Based SAR. Remote Sens. 2014, 6, 1476-1495. https://doi.org/10.3390/rs6021476

AMA Style

Liu D, Shao Y, Liu Z, Riedel B, Sowter A, Niemeier W, Bian Z. Evaluation of InSAR and TomoSAR for Monitoring Deformations Caused by Mining in a Mountainous Area with High Resolution Satellite-Based SAR. Remote Sensing. 2014; 6(2):1476-1495. https://doi.org/10.3390/rs6021476

Chicago/Turabian Style

Liu, Donglie, Yunfeng Shao, Zhenguo Liu, Björn Riedel, Andrew Sowter, Wolfgang Niemeier, and Zhengfu Bian. 2014. "Evaluation of InSAR and TomoSAR for Monitoring Deformations Caused by Mining in a Mountainous Area with High Resolution Satellite-Based SAR" Remote Sensing 6, no. 2: 1476-1495. https://doi.org/10.3390/rs6021476

APA Style

Liu, D., Shao, Y., Liu, Z., Riedel, B., Sowter, A., Niemeier, W., & Bian, Z. (2014). Evaluation of InSAR and TomoSAR for Monitoring Deformations Caused by Mining in a Mountainous Area with High Resolution Satellite-Based SAR. Remote Sensing, 6(2), 1476-1495. https://doi.org/10.3390/rs6021476

Article Metrics

Back to TopTop