A simultaneous inversion for deformation rates and topographic errors of DInSAR data utilizing linear least square inversion technique

We demonstrate here a computer code for calculation of time series and also mean and linear deformation rates from a set of coregistered unwrapped differential interferograms using a linear least-squares inversion technique based on the small baseline subset (SBAS) algorithm. The computer code is written in C and uses a singular value decomposition (SVD) routine from the LAPACK library and the fast Fourier transform for spatial filtering from the FFTW library. Various offset estimation and topographic correction algorithms are implemented, including simultaneous inversion for deformation rates and residual topographic error. This approach is particularly useful when applied to ALOS PALSAR interferograms that are coherent even at large perpendicular baselines and acquired with orbital parameters correlated with the time of acquisition. This methodology is applied to produce time series of ground deformation at Tauhara and Wairakei geothermal fields (Taupo Volcanic Zone, North Island, New Zealand) from 12 ALOS PALSAR images acquired between July 2007 and December 2009. We also present here a high-resolution deformation map of the ground subsidence caused by the extraction of geothermal groundwater for power generation, with maximum rates of subsidence of about 7cm/y.

[1]  Sergey V. Samsonov,et al.  Application of DInSAR-GPS optimization for derivation of three-dimensional surface motion of the southern California region along the San Andreas fault , 2008, Comput. Geosci..

[2]  Riccardo Lanari,et al.  Satellite radar interferometry time series analysis of surface deformation for Los Angeles, California , 2004 .

[3]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[4]  H. Zebker,et al.  Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos , 2007 .

[5]  K. Feigl,et al.  Radar interferometry and its application to changes in the Earth's surface , 1998 .

[6]  D. Sandwell,et al.  Phase gradient approach to stacking interferograms , 1998 .

[7]  Gene H. Golub,et al.  Matrix computations , 1983 .

[8]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[9]  Manabu Watanabe,et al.  ALOS PALSAR: A Pathfinder Mission for Global-Scale Monitoring of the Environment , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[10]  Yuri Fialko,et al.  Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system , 2006, Nature.

[11]  D. Schmidt.  Time-dependent land uplift and subsidence in the Santa Clara Valley , 2003 .

[12]  Andrew Hooper,et al.  A multi‐temporal InSAR method incorporating both persistent scatterer and small baseline approaches , 2008 .

[13]  Rafael C. González,et al.  Local Determination of a Moving Contrast Edge , 1985, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[14]  R. Bürgmann,et al.  Earthquake potential along the northern hayward fault, california , 2000, Science.

[15]  Kristy F. Tiampo,et al.  Detection of displacements on Tenerife Island, Canaries, using radar interferometry , 2004 .

[16]  Gianfranco Fornaro,et al.  A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms , 2002, IEEE Trans. Geosci. Remote. Sens..

[17]  Chris J. Bromley,et al.  Subsidence in the geothermal fields of the Taupo Volcanic Zone, New Zealand from 1996 to 2005 measured by InSAR , 2007 .

[18]  Fuk K. Li,et al.  Synthetic aperture radar interferometry , 2000, Proceedings of the IEEE.

[19]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[20]  T. Farr,et al.  Shuttle radar topography mission produces a wealth of data , 2000 .

[21]  Sergey V. Samsonov,et al.  Analytical optimization of a DInSAR and GPS dataset for derivation of three-dimensional surface motion , 2006, IEEE Geoscience and Remote Sensing Letters.

[22]  Fabio Rocca,et al.  Permanent scatterers in SAR interferometry , 1999, Remote Sensing.

[23]  H. Zebker,et al.  Persistent Scatterer InSAR for Crustal Deformation Analysis , with Application to Volcán Alcedo , Galápagos , 2022 .

[24]  K. Berryman,et al.  A late Quaternary extension rate in the Taupo Volcanic Zone, New Zealand, derived from fault slip data , 2001 .

[25]  Laura M. Wallace,et al.  Subduction zone coupling and tectonic block rotations in the North Island, New Zealand , 2004 .

[26]  Zhong Lu,et al.  Quiescent deformation of the Aniakchak Caldera, Alaska, mapped by InSAR , 2006 .

[27]  Sergey V. Samsonov,et al.  Application of DInSAR-GPS Optimization for Derivation of Fine-Scale Surface Motion Maps of Southern California , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[28]  Riccardo Lanari,et al.  Mining‐related ground deformation in Crescent Valley, Nevada: Implications for sparse GPS networks , 2007 .

[29]  Sergey V. Samsonov,et al.  Topographic Correction for ALOS PALSAR Interferometry , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[30]  Michele Manunta,et al.  A small-baseline approach for investigating deformations on full-resolution differential SAR interferograms , 2004, IEEE Transactions on Geoscience and Remote Sensing.