SPIE Membership Get updates from SPIE Newsroom
  • Newsroom Home
  • Astronomy
  • Biomedical Optics & Medical Imaging
  • Defense & Security
  • Electronic Imaging & Signal Processing
  • Illumination & Displays
  • Lasers & Sources
  • Micro/Nano Lithography
  • Nanotechnology
  • Optical Design & Engineering
  • Optoelectronics & Communications
  • Remote Sensing
  • Sensing & Measurement
  • Solar & Alternative Energy
  • Sign up for Newsroom E-Alerts
  • Information for:
    Advertisers
SPIE Photonics West 2018 | Call for Papers

SPIE Journals OPEN ACCESS

SPIE PRESS

SPIE PRESS

Print PageEmail PageView PDF

Remote Sensing

Remote sensing of ground deformation

An unconventional synthetic aperture radar interferometry technique can be used to measure deformation of vegetated regions.
18 March 2013, SPIE Newsroom. DOI: 10.1117/2.1201303.004772

Synthetic Aperture Radar Interferometry (InSAR) is a remote sensing technique that provides precise measurements of land surface deformation over large areas and at high spatial resolution.1–5 Land surface deformation is the changing shape of the ground (topographic changes) along fault lines caused by earthquakes. Deformation maps, which are important from a disaster mitigation perspective, of an area can be obtained by taking the difference of the phase measurements from InSAR images acquired at two different points in times. The phase difference is known as an interferometric phase or interferogram.

Different components of interferograms are produced by different sources. One such component is introduced by the topography of an area and can be removed using an external digital elevation model (DEM) to give a differential interferogram. The effect of the Earth's curvature also needs to be modeled and removed to produce a flattened interferogram. However, InSAR fails to measure the deformation in vegetated areas because these areas suffer from significant backscattering effect changes over time. This effect leads to a decrease in coherence between the different phase measurements and is known as temporal decorrelation.

The technique of persistent scatterer interferometry (PSI)6 addresses the decorrelation problem by identifying persistent scatterer (PS) points at which the deformation is estimated over long periods of time. PS points tend to be man-made structures, such as buildings and bridges, or natural features, such as large rocks or facets, within an area of vegetation. The returned signal amplitude from the PS points is high and constant over time. Here, we use a PSI approach whereby time series of interferograms are analyzed to investigate temporal variations in the backscattering behavior of the PS features. This is in contrast to conventional InSAR techniques in which only a few images are used. Most PSI methods involve an initial selection of PS points—known as candidate pixels—that are chosen based on amplitude information. Subsequent and detailed processing of phase information is used to identify the final PS points.

Our study area—in the La González basin of the central Venezuelan Andes (see Figure 1)—is covered by dense vegetation ideal for PSI-based measurements. The region is a typical pull-apart basin that lies on a releasing bend of the right-lateral strike-slip Boconó Fault System.7 We focused on measuring surface deformation caused by the interaction between the right-lateral movement of the Boconó fault and the active regional northwest to southeast shortening across the mountain range.


Figure 1. The study area located in the La González basin, Venezuela. The black rectangle (61×52 km) denotes the area from which radar data was processed.

We used PSI, specifically the Stanford Method for PS (StaMPS), to address the temporal decorrelation problem in our study area by identifying features for which deformation was estimated over a long period of time.8 This algorithm uses the spatial correlation between pixels and works well even in areas that lack man-made structures. The final PS pixels are identified by making an estimate of the phase stability of candidate pixels, through phase analysis.8 The ground deformation rate is then estimated using the results of the time series analysis and by applying the phase information from the PS pixels.

The available radar dataset (see Figure 2) for our study consisted of 10 Envisat Advanced Synthetic Aperture Radar (ASAR) images collected between 3 January 2003 and 5 August 2005. The subset of the data frame that covers our study area is indicated in Figure 1. In our PSI processing, an image that maximizes the stack coherence (a function of temporal and perpendicular baselines and of Doppler centroid frequency) was selected as the master.9, 10 Nine single master interferograms were then generated. We used the Shuttle Radar Topography Mission (SRTM) DEM to remove the main topographic contribution to the interferometric phase. We identified 53,864 PS pixels during processing and we constructed a velocity map (see Figure 3) from the time series analysis results. We estimated and removed all sources of atmospheric, topographic, and orbital errors from our results, and show that the area undergoes significant deformation that is likely caused by tectonic activity. We estimated the maximum deformation rate in the area to be 40mm per year (see Figure 3). However, deformation in this area is complicated, as illustrated in Figure 4. Points A and B (locations shown in Figure 3) are moving away from, and towards, the observing satellite at rates of 17 and 5mm per year, respectively.


Figure 2. Acquisition geometry of the available radar data in the area. Acquisition time versus spatial baseline. Solid lines represent the generated interferograms.

Figure 3. Estimated line-of-sight deformation rate for the area within the black rectangle shown in Figure 1. The deformation time series of points A and B are shown in Figure 4.

Figure 4. Deformation time series at points A (top) and B (bottom) whose locations are shown in Figure 3. InSAR: Synthetic Aperture Radar Interferometry. LOS: line-of-sight.

We have successfully used the StaMPS PSI method to detect and measure ground deformation caused by tectonic activity in the La González basin, Venezuela, which is covered by dense vegetation. We plan to decompose the line-of-sight deformation rate into horizontal and vertical components by using at least two different imaging geometries. This may allow us to separate different tectonic regimes that are active in this region.

We thank the European Space Agency (ESA) for providing Envisat ASAR data.


Maryam Dehghani
Shiraz University
Shiraz, Iran

Maryam Dehgahani received her BS degree in geodesy and geomatics engineering from Iran University of Science & Technology. She obtained her MSc and PhD degrees in remote sensing from K. N. Toosi University of Technology in 2004 and 2009, respectively. She is currently an assistant professor in the Department of Civil and Environmental Engineering.

Hamid Reza Javadi
Geological Survey of Iran
Tehran, Iran

Hamid Reza Javadi received his BS degree in geology from Gorgan University, and his MSc degree in structural geology and tectonics from the Research Institute for Earth Sciences in 2006. His main research interests are structural geology and active tectonics. He is currently a geological expert for the Geological Survey of Iran.


References:
1. D. Massonnet, M. Rossi, C. Carmona, F. Adragna, G. Peltzer, K. Feigl, T. Rabaute, The displacement field of the Landers earthquake mapped by radar interferometry, Nature 364, p. 138-142, 1993. doi:10.1038/364138a0
2. G. Peltzer, P. Rosen, F. Rogez, K. Hudnut, Poroelastic rebound along the Landers 1992 earthquake surface rupture, J. Geophys. Res. 103(B12), p. 30131-30145, 1998. doi:10.1029/98JB02302
3. B. Fruneau, F. Sarti, Detection of ground subsidence in the city of Paris using radar interferometry: isolation from atmospheric artifacts using correlation, Geophys. Res. Lett. 27(24), p. 3981-3984, 2000. doi:10.1029/2000GL008489
4. M. Tesauro, P. Beradino, R. Lanari, E. Sansoti, G. Fornaro, G. Franceschetti, Urban subsidence inside the City of Napoli (Italy) observed with synthetic aperture radar interferometry at Campi Flegrei caldera, Geophys. Res. Lett. 27(13), p. 1961-1964, 2000. doi:10.1029/2000GL008481
5. M. Crosetto, C. C. Tscherning, B. Crippa, M. Castillo, Subsidence monitoring using SAR interferometry: reduction of the atmospheric effects using stochastic filtering, Geophys. Res. Lett. 29(9), p. 26.1-26.4, 2002. doi:10.1029/2001GL013544
6. A. Ferretti, C. Prati, F. Rocca, Permanent scatterers in SAR interferometry, IEEE Trans. Geoscience Remote Sensing 39(1), p. 8-20, 2001. doi:10.1109/36.898661
7. C. Schubert, Late-Cenozoic pull-apart basins, Boconó fault zone, Venezuelan Andes, J. Struct. Geol. 2(4), p. 463-468, 1980. doi:10.1016/0191-8141(80)90008-5
8. A. Hooper, P. Segall, A. Zebker, Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcan Alcedo, Galapagos, J. Geophys. Res. 112(B7), p. B07407, 2007. doi:10.1029/2006JB004763
9. B. M. Kampes, Displacements Parameter Estimation using Permanent Scatterer Interferometry, PhD thesis, DLR, Bibliotheks- und Informationswesen, Stuttgart, 2005.
10. G. Ketelaar, Satellite Radar Interferometry Subsidence Monitoring Techniques, Delft University of Technology, the Netherlands, 2009.