Optimal fault resolution in geodetic inversion of coseismic data

SUMMARY With the continued growth in availability of differential interferometry from synthetic aperture radar and GPS data, space based geodesy has been widely applied to image the coseismic displacement field and to retrieve the static dislocation over the fault plane for almost all the significant earthquakes of the past two decades. This is performed by linear data inversion over a set of subfaults, generally characterized by a constant and predefined or manually adjusted dimensions. In this paper, we propose a new algorithm to automatically retrieve an optimized fault subdivision in the linear inversion of coseismic geodetic data. The code iteratively keeps the parameter resolution close to a predefined high value. We first discuss the rationale supporting ouralgorithmand,afteradetaileddescriptionofitsimplementation,weanalysetheadvantages of its introduction in the data inversion. The algorithm was tested against an exhaustive range of synthetic and real data sets and fault mechanisms. Among them, we present the results for theMw 6.2, 2009 L’Aquila (Central Italy) earthquake and compare the new and previously published slip distributions showing the disappearance of misleading slip pattern and the increased resolution for shallower zones.

[1]  Henriette Sudhaus,et al.  Improved source modelling through combined use of InSAR and GPS under consideration of correlated data errors: application to the June 2000 Kleifarvatn earthquake, Iceland , 2009 .

[2]  Christian Bignami,et al.  X-, C-, and L-Band DInSAR Investigation of the April 6, 2009, Abruzzi Earthquake , 2011, IEEE Geoscience and Remote Sensing Letters.

[3]  E. Trasatti,et al.  Numerical inversion of deformation caused by pressure sources: Application to Mount Etna (Italy) , 2008 .

[4]  L. Rivera,et al.  Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations , 2002 .

[5]  Rowena B. Lohman,et al.  Some thoughts on the use of InSAR data to constrain models of surface deformation: Noise structure and data downsampling , 2005 .

[6]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[7]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[8]  Y. Takane,et al.  Generalized Inverse Matrices , 2011 .

[9]  Alberto Michelini,et al.  The 2009 L'Aquila (central Italy) MW6.3 earthquake: Main shock and aftershocks , 2009 .

[10]  R. Goldstein,et al.  Mapping small elevation changes over large areas: Differential radar interferometry , 1989 .

[11]  Rongjiang Wang,et al.  Erratum to: "Computation of deformation induced by earthquakes in a multi-layered elastic crust - FORTRAN programs EDGRN/EDCMP": [Computers & Geosciences, 29(2) (2003) 195-207] , 2006, Comput. Geosci..

[12]  F. Radicioni,et al.  FAST TRACK PAPER: Coseismic and initial post-seismic slip of the 2009 Mw 6.3 L'Aquila earthquake, Italy, from GPS measurements , 2010 .

[13]  R. Hanssen.  Radar Interferometry: Data Interpretation and Error Analysis , 2001 .

[14]  H. Zebker,et al.  Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements , 2002 .

[15]  Morgan T. Page,et al.  Constraining earthquake source inversions with GPS data: 1. Resolution‐based removal of artifacts , 2009 .

[16]  Zhong Lu,et al.  Source model for the Mw 6.7, 23 October 2002, Nenana Mountain Earthquake (Alaska) from InSAR , 2003 .

[17]  K. Feigl,et al.  The displacement field of the Landers earthquake mapped by radar interferometry , 1993, Nature.

[18]  T. Sagiya,et al.  Coseismic slip resolution along a plate boundary megathrust: The Nankai Trough, southwest Japan , 1999 .

[19]  T. Wright,et al.  Fault identification for buried strike-slip earthquakes using InSAR: The 1994 and 2004 Al Hoceima, Morocco earthquakes , 2006 .

[20]  James Jackson,et al.  Surface displacements and source parameters of the 2003 Bam (Iran) earthquake from Envisat advanced synthetic aperture radar imagery , 2005 .

[21]  D. Sandwell,et al.  Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit , 2005, Nature.

[22]  Yukitoshi Fukahata,et al.  General expressions for internal deformation fields due to a dislocation source in a multilayered elastic half-space , 2005 .

[23]  Paul Segall,et al.  Fault geometry and slip distribution of the 1999 Chi‐Chi, Taiwan Earthquake imaged from inversion of GPS data , 2001 .

[24]  Paul A. Rosen,et al.  Co-seismic slip from the 1995 July 30 Mw= 8.1 Antofagasta, Chile, earthquake as constrained by InSAR and GPS observations , 2002 .

[25]  Grazia Pietrantonio,et al.  Coseismic deformation of the destructive April 6, 2009 L'Aquila earthquake (central Italy) from GPS data , 2009 .

[26]  Christian Bignami,et al.  Finite fault inversion of DInSAR coseismic displacement of the 2009 L'Aquila earthquake (central Italy) , 2009 .

[27]  WangRongjiang,et al.  Computation of deformation induced by earthquakes in a multi-layered elastic crust , 2003 .

[28]  K. Feigl,et al.  Coseismic and Postseismic Fault Slip for the 17 August 1999, M = 7.5, Izmit, Turkey Earthquake. , 2000, Science.

[29]  R. Penrose.  A Generalized inverse for matrices , 1955 .

[30]  M. Page,et al.  Constraining earthquake source inversions with GPS data: 2. A two‐step approach to combine seismic and geodetic data sets , 2009 .

[31]  Riccardo Lanari,et al.  Synthetic Aperture Radar Processing , 1999 .

[32]  Anthony Lomax,et al.  Rupture history of the 2009 L'Aquila (Italy) earthquake from non‐linear joint inversion of strong motion and GPS data , 2009 .

[33]  Xiaojing Li,et al.  Fault Geometry and Slip Distribution of the 2010 Yushu Earthquakes Inferred from InSAR Measurement , 2011 .