The exploitation of underground mineral resources brings great material wealth but also damages the natural environment at the same time. For example, it can cause mine subsidence, which brings a series of severe consequences, such as decrease in cultivated fields, collapse of buildings, and damage to roads. Accordingly, it is urgent and necessary to effectively monitor mine subsidence.12.–3 Differential interferometric synthetic aperture radar (DInSAR) has proven very useful in monitoring mine subsidence. According to the different methods of removing topographical information, DInSAR can be divided into two-, three-, and four-pass. For two-pass DInSAR, it is necessary to simulate topographic phase by an external digital elevation model (DEM) and remove it from SAR interferometric phase to extract land subsidence phases and values. Therefore, external DEM error will inevitably lead to land subsidence error.4,5 In other words, due to the variation of error of different ground resolution in each external DEM, the land subsidence phases and values extracted from the same SAR interferometric pair will vary. Liu6 reported that the deformation error introduced by external DEM error was mainly systematic. Long et al.7,8 determined that the effect of orbit error on deformation phases was greater than the external DEM error. Hirt et al.9 pointed out that the accuracy of SRTM DEM () was higher than ASTER GDEM () in sparsely populated areas, and SRTM DEM was recommended for use in most applications. Taking typical flat and mountainous areas of Qinghai Tibet Plateau as research areas, Li et al.10 compared and analyzed three land subsidence results monitored by two-pass DInSAR with SRTM, ASTER GDEM, and external DEM, and then concluded that SRTM DEM was more suitable for the external reference DEM.
In this study, the effect of external DEM on monitored mine subsidence is systematically investigated from theory and in practical applications. The relationship of external DEM errors and monitored deformation errors is derived based on the principles of interferometric synthetic aperture radar (InSAR) and two-pass DInSAR. We then determine the effect of different DEM errors on monitored deformation results under the constant image parameters in theory. In this study, 25-, 30-, and 90-m ground resolution external DEMs covering Dongtan and Yangcun mine areas of Jining are resampled to the same SAR image coordinate system to generate three DEMs with the same coordinate system and resolution. The elevation values of three external DEMs before and after resampling are compared with each other and with the leveling-measured elevation values. Changes in DEM value and accuracy before and after resampling were found to be small, but larger differences existed between three external DEMs, although in the flat test areas of this study. The Advanced Land Observing Satellite Phased Array L-band SAR (ALOS PALSAR) interferometric pairs covering the Dongtan and Yangcun mine areas are processed using two-pass DInSAR with three external DEMs. Some key images in the processing are given, including original and enhanced differential interferograms, coherence images, unwrapped phase images, and subsidence images. All of these aim to compare and analyze the effect of three external DEMs in high- and low-coherence regions with different subsidence magnitude. In addition, the mine subsidence monitored by DInSAR using three external DEMs is validated further by leveling-measured mine subsidence at 23 leveling benchmarks in Dongtan mine and 29 leveling benchmarks in Yangcun mine at the nearest time interval. Research results show that the effect of external DEM on the monitored results of two-pass DInSAR is not only related to imaging parameters, such as radar look angle, perpendicular baseline, slant range, and external DEM error, but also related to the subsidence magnitude and coherence of test areas. Generally, in high-coherence regions with larger subsidence magnitude, the subsidence trends and values obtained by three external DEMs are more consistent, but compared with leveling data, the subsidence trends are consistent while subsidence values are quite different. In low-coherence regions with slight subsidence, the subsidence values obtained by three external DEMs vary greatly and are clearly different from the leveling-measured subsidence values.
InSAR acquires two images over the same test area from slightly different viewpoints and determines the relative height of the meeting point using the interferometric phase of each pixel in the two images (Fig. 1).11 In Fig. 1, and are different spatial positions of the satellite when the ground point is observed during repeat flights, is the spatial interferometric baseline vector, is the inclination angle of baseline, is the height difference between the platform and ground reference plane, and are the slant ranges of ground point to the platform, respectively, and are the respective look angles, and is the height of relative to the reference plane.
Assuming , the interferometric phase can be written as12
According to the geometric relations shown in Fig. 1, and considering that and are much smaller than , the approximate value of can be derived by the cosine theorem as follows:
So there are
The DInSAR technique is developed on the basis of InSAR. If the land subsidence occurs during the acquisition period of the SAR image pair, the interferometric phase of the target point is composed of two parts and can be rewritten as13
From Eqs. (2) and (6), it can be shown thatFig. 2.
Study Area and Data Source
Jining is one of China’s main coal-producing regions and is located in the north of the country. Overexploitation of this resource to meet industrial needs has led to severe land subsidence in many areas of Jining. The Dongtan and Yangcun mine areas are selected as the test areas and shown in Fig. 3. Mine subsidence has caused a series of problems, such as cultivated fields decrease, buildings collapse, roads crack, and other problems developing across the two test areas.
To analyze the effect of external DEM on mine subsidence monitored by two-pass DInSAR, two interferometric pairs of ALOS PALSAR single look complex (SLC) images (see Fig. 3) and three external DEMs with different resolutions are collected.
ALOS PALSAR SLC Images
The ALOS was Japan’s new-generation Earth Observation satellite, launched in January 2006 by the Japan Aerospace Exploration Agency. The Phased Array L-band SAR (PALSAR) carried by the ALOS satellite was an active microwave imaging radar that could image in all weather and at all times, and had been widely applied in military and civilian fields. The design life of ALOS was three years, but it was in service for five years and stopped on May 12, 2011, due to a power generation anomaly that led to the loss of communication functions. Considering the time span of the leveling data collected in the test areas, four ALOS PALSAR images were selected. Their characteristics are HH polarization modes, descending orbits, pixel spacing of 4.68 m in slant range, and 3.17 m in azimuth. Table 1 shows the basic parameters of four ALOS PALSAR images
Basic parameters of four ALOS PALSAR images.
|Pairs||Date||Orbit||Center longitude (°)||Center latitude (°)||Center slant range (m)||Center look angle (deg)||Perpendicular baseline (m)|
|1||Master||January 10, 2009||15791||116.492||35.656||886,390||38.735||195|
|Slave||February 25, 2009||16462||116.497||35.649||886,235||38.737|
|2||Master||February 11, 2010||21582||116.843||35.450||853,576||38.714||181|
|Slave||March 29, 2010||22253||116.844||35.450||853,473||38.307|
1∶5 million DEM
A million DEM covering test areas was obtained, and a mosaic and coordinate transformation between Xi’an 80 coordinate system and WGS84 coordinate system was implemented. The ground resolution of the DEM is 25 m with accuracy of 4 m on flat ground and 11 m in mountainous areas.
ASTER GDEM V2 DEM
An Advanced Spaceborne Thermal Emission and Reflection Radiometer Global DEM (ASTER GDEM) covering the test areas was downloaded from the website (Ref. 14), the geospatial data cloud of the Computer Network Information Center (CNIC). The DEM is publicly released by the Japanese Ministry of Economy, Trade, and Industry and the American National Aeronautics and Space Administration (NASA), and is the only elevation images data set with high resolution covering the global land surface. Its ground resolution is 30 m with average accuracy of 20 m.15,16
A Shuttle Radar Topography Mission (SRTM) DEM covering test areas was also downloaded from the geospatial data cloud of CNIC.14 It was jointly measured by American NASA and the National Imagery and Mapping Agency in February 2000. Its ground resolution is 90 m with an average accuracy of 16 m.17,18
Figure 4 shows the external DEMs covering test areas.
The above three external DEMs are resampled to the same SAR image coordinate system to give them a consistent coordinate system and resolution. Therefore, their difference and accuracy can be compared effectively to better assess the effect of using different external DEMs on mine subsidence monitored by two-pass DInSAR.
Figure 5 shows the three external DEMs after resampling. Table 2 shows the difference of resampled DEM elevation values of 12 even-selected pixels marked by red plus symbols in Fig. 5. Table 3 shows the difference of elevation values derived from leveling measurement and the three external DEMs before and after resampling. The locations of 12 leveling benchmarks are marked by black points in Fig. 5.
Elevation values of 12 pixels over three resampled external DEMs.
|Number||Elevation value (m)||Number||Elevation value (m)||Number||Elevation value (m)|
|25 m||30 m||90 m||25 m||30 m||90 m||25 m||30 m||90 m|
Elevation values derived from leveling measurements and three external DEMs before and after resampling.
|Number||Elevation value (m)||Number||Elevation value (m)|
|Leveling||25-m B/A||30-m B/A||90-m B/A||Leveling||25-m B/A||30-m B/A||90-m B/A|
In Tables 2 and 3, “25,” “30,” and “90 m” represent the million DEM with 25-m ground resolution, the ASTER GDEM V2 DEM with 30-m ground resolution, and the SRTM DEM with 90-m ground resolution, respectively. In Table 3, “B/A” is the abbreviation of “before resampling/after resampling” and represents the elevation values of three external DEMs before and after resampling.
As can be seen from Figs. 4 and 5, and Tables 2 and 3, the terrain of the test areas slopes gently, but three main external DEMs of million DEM, ASTER GDEM V2 DEM, and SRTM DEM show very large differences. The distributed intervals of elevation values of 25-, 30-, and 90-m ground resolution DEM after resampling are (42, 54 m), (, 72 m), and (20, 70 m), respectively. Comparatively speaking, the DEM elevation values of million DEM are closer to those of SRTM DEM with the maximum and average deviation of 12 and 3 m among 12 pixels, but for ASTER GDEM V2 DEM, the deviations are larger. Specifically, the maximum and average deviations are about 15 and 10 m compared with 1:5 million DEM, and about 17 and 7 m compared with SRTM DEM. The elevation values of the 25-m resolution external DEM are nearest to the 12 leveling-measured elevation values, followed by the 90-m resolution external DEM, while the values of the 30-m resolution external DEM are most different from the leveling-measured elevation values and therefore have the lowest accuracy. It is true that the value and accuracy (i.e., the DEM error) of a DEM will change before and after resampling. However, the degree of change is different in each case. The value and accuracy of the 90-m resolution external DEM before and after resampling change most, followed by the 30-m resolution external DEM, and there is no difference between the value and accuracy of the 25-m resolution external DEM before and after resampling.
Data Processing Method
For differential processing of two-pass DInSAR, the slave image and external DEM are first coregistered and resampled into the master image space to ensure that they share exactly the same SAR coordinates. The second step of the method is to deal with differential processing and remove the flat earth effect to obtain the differential interferogram. Then adaptive filtering must be used to acquire the enhanced differential interferogram and coherence image. Next, the minimum cost flow phase unwrapping method is used to unwrap the enhanced differential interferometric phase. This step is followed by transforming phase to height, geocoding, and then determining the amount of mine subsidence.19 The primary procedure for the first interferometric pair is shown in Fig. 6.
In addition to DEM error, the subsidence measurements derived by two-pass DInSAR are often degraded by baseline estimation, atmospheric delay, orbital parameters, noise, and other errors. The corresponding parameters and processing methods are consistent when two interferometric pairs are processed using two-pass DInSAR with three external DEMs, to ensure that each interferogram has the same baseline estimation, atmospheric delay, orbital parameters, noise, and other errors. This can guarantee that the final monitored subsidence difference depends only on the use of three external DEMs. Moreover, as we all know that the atmospheric effect is one of the most intractable problems in the processing of two-pass DInSAR to monitor mine subsidence, it is difficult to implement accurate estimation and removal of the atmospheric effect in the absence of priori data. Our two interferometric pairs of test areas in this study are processed using standard processing procedure of two-pass DInSAR without separately considering the atmospheric effect.
Test Results and Analysis
Land Subsidence in the Dongtan Mine Area
According to the processing procedure of test data in Fig. 6, the first interferometric pair is processed. Figure 7 shows the enhanced differential interferogram derived from two-pass DInSAR with 25-m ground resolution DEM. Its background is the SAR intensity image of the test area.
Considering the collected leveling-measured data at 23 leveling benchmarks, subsidence of region A in the Dongtan mine area is obtained and analyzed to study the effect of three external DEMs on monitored mine subsidence in the high-coherence region with larger subsidence magnitude. Figure 8 shows the local enlarged images of region A, from left to right sequentially corresponding to the original and enhanced differential interferograms, coherence images, unwrapped phase images, subsidence images, and from top to bottom, corresponding to the differential results of three external DEMs with 25-, 30-, and 90-m ground resolution.
As can be seen from Fig. 8, the coherences of region A are higher (average coherence coefficients are all ), and three complete and clear interferometric fringes are formed, which means that the largest subsidence magnitude, in the vertical direction, can be monitored by two-pass DInSAR. Moreover, the three subsidence results monitored are continuous without intermittent spots and thus are all highly reliable and accurate.
We also can see from Fig. 8, the uplift phenomenon occurred during this imaging period in our test area. Actually, the uplift phenomenon usually occurs at the edge of subsidence regions and it is normal. Mining of coal resources will lead to crustal movement, so that some parts of the ground uplift while some parts maybe subside, and the subsidence parts will be surrounded by those uplift parts. The amount of uplift is related to the strata structure and it is generally not large. Only when the quaternary loose layer is relatively thick, there will be no uplift phenomenon.
As can be seen from Fig. 9, three fan-shaped subsidence areas (broad in the west and narrow in the east) are monitored by two-pass DInSAR with three external DEMs, which are consistent in subsidence locations, distributions, and trends. These display differences mainly in three aspects: (1) values of monitored subsidence are different. For example, the largest monitored subsidence values are , , and , respectively. (2) Coverage areas of monitored subsidence are obviously different. Specifically, the areas from largest to smallest are the differential results of three external DEMs with 30-, 25-, and 90-m ground resolution. (3) Continuities of monitored subsidence are different. The continuity of subsidence monitored by 90-m resolution external DEM is worst, especially at the edge of the subsidence areas, and spots are more obvious, intermittent, and scattered. It follows that DEMs with high resolution are beneficial to maintain the continuity and integrity of subsidence values, especially in the edge regions.
We collected leveling-measured mine subsidence values at 23 leveling benchmarks between January 9, 2009, and February 24, 2009, in the Dongtan mine area, shown in Fig. 9 with black points and numbered in sequence of from left to right and from bottom to top. Table 4 and Fig. 10 show the comparisons between four types of subsidence values derived from leveling and from two-pass DInSAR with three external DEMs.
Comparisons between four types of subsidence values monitored by DInSAR and leveling of Dongtan mine.
|Number||Leveling (cm)||DInSAR (cm)||Number||Leveling (cm)||DInSAR (cm)|
|25 m||30 m||90 m||25 m||30 m||90 m|
As can be seen from Table 4 and Fig. 10, the 18 benchmarks numbered by 1 to 12 and 18 to 23, the subsidence trends, and values monitored by two-pass DInSAR with three external DEMs show good agreement. For example, the statistical calculation indicates that the largest and mean differences between two types of subsidence monitored by two-pass DInSAR with 25- and 30-m ground resolution DEM are and . Compared with leveling-measured subsidence, the subsidence trends gradually increasing from the edge toward the center of subsiding regions are monitored, but the values of subsidence given from the DInSAR data are much larger than those determined by leveling. The distributed intervals of subsidence values monitored by two-pass DInSAR with three external DEMs and leveling are , , , and , respectively. In addition, for the five benchmarks numbered 13 to 17 with larger subsidence values distributing in , the subsidence trends and values monitored by DInSAR with three external DEMs show obvious differences, especially in the 16th benchmark. Compared with leveling-measured subsidence, the subsidence trend monitored by DInSAR with 30-m resolution external DEM is basically consistent and the difference of values is relatively smaller, while the difference of subsidence trend and values of 25-m resolution DEM is larger, However, the subsidence trend of 90-m resolution DEM is completely inconsistent and the subsidence values are also quite different. Finally, leveling benchmarks are not arranged along the center line of subsidence, so it is not possible to determine the real central location, spatial distribution, maximum magnitude of mine subsidence, and other details using only the subsidence values of these leveling points. DInSAR data, however, address the above shortcomings and reconstruct intuitively the distribution, trend, and other subsidence information during the imaging period. However, two types of subsidence magnitude derived from two-pass DInSAR and from leveling cannot agree with each other very well mainly because, first, the period of leveling-measured subsidence cannot be completely consistent with that of DInSAR, and, second, leveling-measured subsidence can only reflect the subsidence of a “point,” whereas DInSAR-monitored subsidence reflects the subsidence of a “surface,” of in this case. Thus, the comparison between the two types of subsidence in this test is to compare a point with a surface.
Land Subsidence in the Yangcun Mine Area
The second interferometric pair was also processed using the processing procedure for test data shown in Fig. 6. Figure 11 shows the enlarged images of the Yangcun mine area, from left to right sequentially corresponding to the original and enhanced differential interferograms, coherence images, unwrapped phase images, subsidence images, and from top to bottom, respectively, corresponding to the differential results of three external DEMs with 25-, 30-, and 90-m ground resolution.
As can be seen from Fig. 11, the coherences of the subsidence area are lower (average coherence coefficients are all about 0.3), no interferometric fringes are formed, the monitored subsidence results are spotty and discontinuous, and all are, therefore, less reliable and accurate than those of the Dongtan area.
As can be seen from Fig. 12, the two subsidence areas on the left and right monitored by two-pass DInSAR with three external DEMs agree with the two locations of low-coherence regions of Fig. 11(c). It is therefore considered that the areas of subsidence can be detected by DInSAR in the low-coherence regions. In addition, the differences of subsidence values monitored in the low-coherence regions by two-pass DInSAR with three external DEMs are greater than those in high-coherence regions. For example, the largest monitored subsidence values are , , and , respectively.
We collected leveling-measured mine subsidence values at 29 leveling benchmarks between February 12, 2010, and March 22, 2010, in the area of the Yangcun mine, and the layout of these benchmarks is shown in Fig. 12 with black points and numbered in sequence of from top to bottom and from left to right. Table 5 and Fig. 13 show the comparisons between four types of subsidence derived from leveling and from two-pass DInSAR with three external DEMs.
Comparisons between the four types of subsidence of Yangcun mine.
|Number||Leveling (cm)||DInSAR (cm)||Number||Leveling (cm)||DInSAR/cm|
|25 m||30 m||90 m||25 m||30 m||90 m|
As can be seen from Table 5 and Fig. 13, the monitored values of 29 benchmarks show that smaller subsidence and uplift appear in this area. The statistical calculation indicates that subsidence appears in 17 benchmarks, with a maximum magnitude of , uplift appears in 10 benchmarks, with a maximum magnitude of 0.9 cm, and 2 benchmarks show no measured change in elevation. For three differential results of two-pass DInSAR with 25-, 30-, and 90-m ground resolution DEMs at locations corresponding to the 29 benchmarks, the result determined from the 25-m DEM indicates that subsidence occurs around 10 points, with a maximum of and uplift appears around 19 points, with a maximum of 6.4 cm. The result determined from the 30-m DEM indicates that subsidence occurs around 9 points, with a maximum of and uplift appears around 19 points, with a maximum of 11.5 cm and there is no measurable change around one point. The result determined from the 90-m DEM indicates that uplift occurs around all 29 points, with a maximum value of 5.5 cm. In low-coherence regions with small subsidence magnitude, the subsidence values monitored by two-pass DInSAR with three external DEMs show larger differences with each other, and all of them are relatively inconsistent with leveling-measured subsidence values. In other words, two-pass DInSAR with three external DEMs did not derive the relatively reliable subsidence values in these areas. Comparatively speaking, however, the subsidence trends and values monitored by two-pass DInSAR with 25-m ground resolution DEM show best agreement with those of the 29 benchmarks and the differences are relatively smaller in benchmarks 1, 2, 4, 9, and 16 to 29. The subsidence values measured using the 30-m ground resolution DEM show smaller differences in the benchmarks numbered 14 and 16 to 29. The subsidence values measured using the 90-m ground resolution DEM change gently and show large differences from values measured at the 29 benchmarks.
This study used different external DEMs in monitoring mine subsidence by two-pass DInSAR, leading to the following specific conclusions:
1. The values and errors of DEM are changed little before and after resampling. Three external DEMs show very large differences in flat test areas. Comparatively speaking, million DEM are closer to SRTM DEM and both show highest precision, ASTER GDEM V2 DEM before and after resampling is most different from million DEM, SRTM DEM, and leveling, with lowest accuracy.
2. In high-coherence regions, the magnitudes of subsidence monitored by two-pass DInSAR with three external DEMs are consistent in locations, distributions, and trends, but are different in values, coverage areas, and continuities. Generally, DEMs with high resolution are beneficial to maintain the continuity and reliability of subsidence values and provide richer subsidence information.
3. In high-coherence regions with smaller subsidence magnitudes, the monitored subsidence values of three external DEMs show good agreement. However, all of the subsidence values are much larger than those measured by leveling. In high-coherence regions with larger subsidence magnitudes, the subsidence values measured using the external DEMs are not consistent and all of them are basically smaller than those measured by leveling. Due to the different monitored periods and comparison between a point and a surface we consider that the differences between two types of subsidence values monitored by two-pass DInSAR and leveling are reasonable.
4. In the low-coherence regions with slight subsidence, the monitored subsidence values of three external DEMs show relatively large inconsistency and differ from the leveling-measured subsidence values. So we can think that they are therefore relatively unreliable. However, the monitored subsidence locations of three external DEMs are relatively consistent and all agree with locations of low-coherence regions. Then, we conclude that subsidence locations can be detected by DInSAR in these low-coherence regions. This conclusion may be helpful for the detection of illegal mining by means of two-pass DInSAR.
5. In theory, the deformation error is related only to radar look angle, perpendicular baseline, slant range, and external DEM error. However, the actual application tests show that it is also related to the ground resolution of the DEM, the magnitude of subsidence, and the coherence of test areas. Specifically, DEMs with high resolution are most useful and provide richer subsidence information. In high-coherence regions with smaller subsidence magnitudes, the effect of different external DEMs on monitored subsidence is relatively small and negligible. In high-coherence regions with larger subsidence magnitudes and low-coherence regions with small subsidence magnitudes, the effect is larger and researchers should use the DEM with higher ground resolution and precision.
This study was supported by the National Natural Science Foundation of China (Nos. 41274007 and 41404003), the Taishan Scholars Program of Shandong Province, China (No. TSXZ201509), the Scientific Research Foundation for the introduction of talent of Shandong University of Science and Technology (No. 2014RCJJ047), the Innovation Project for the Graduates of Shandong University of Science and Technology (Nos. SDKDYC170206 and SDKDYC170208).
Qiuxiang Tao has received her BEng degree in surveying engineering, her MEng degree in geodesy and surveying engineering, and PhD in photogrammetry and remote sensing from the College of Geomatics, Shandong University of Science and Technology, Qingdao, China, in 2001, 2004, and 2009, respectively. She is currently an associate professor in the College of Geomatics, Shandong University of Science and Technology. Her research interests include interferometric synthetic aperture radar (InSAR), differential interferometric synthetic aperture radar, short spatio-temporal PS InSAR technology, and their applications.
Tengfei Gao is currently pursuing her MEng degree in geodesy and surveying engineering at Shandong University of Science and Technology, Qingdao, China. Her research interests include InSAR technology and its application.
Guolin Liu has received his PhD in surveying engineering from China University of Mining and Technology, Xuzhou, China, in 1999. He is currently a professor in the College of Geomatics, Shandong University of Science and Technology. His research interests include the theories and applications of modern surveying data processing, InSAR technology, and its application.
Zhiwei Wang has received his BEng degree in surveying engineering and his MEng degree in geodesy and surveying engineering from Shandong University of Technology, Zibo, China, in 2011 and 2012, respectively. He is currently pursuing his PhD in geodesy and surveying engineering, Shandong University of Science and Technology, Qingdao, China. His current research interests include InSAR technology and its application.