Co-seismic slip from the 1995 July 30 Mw= 8.1 Antofagasta, Chile, earthquake as constrained by InSAR and GPS observations

SUMMARY We analyse radar interferometric and GPS observations of the displacement field from the 1995 July 30 Mw = 8.1 Antofagasta, Chile, earthquake and invert for the distribution of slip along the co-seismic fault plane. Using a fixed fault geometry, we compare the use of singular-value decomposition and constrained linear inversion to invert for the slip distribution and find that the latter approach is better resolved and more physically reasonable. Separate inversions using only GPS data, only InSAR data from descending orbits, and InSAR data from both ascending and descending orbits without the GPS data illustrate the complimentary nature of GPS and the presently available InSAR data. The GPS data resolve slip near GPS benchmarks well, while the InSAR provides greater spatial sampling. The combination of ascending and descending InSAR data contributes greatly to the ability of InSAR to resolve the slip model, thereby emphasizing the need to acquire this data for future earthquakes. The rake, distribution of slip and seismic moment of our preferred model are generally consistent with previous seismic and geodetic inversions, although significant differences do exist. GPS data projected in the radar line-of-sight (LOS) and corresponding InSAR pixels have a root mean square (rms) difference of about 3 cm. Comparison of our predictions of vertical displacement and observed uplift from corraline algae have an rms of 10 cm. Our inversion and previous results reveal that the location of slip might be influenced by the 1987 Mw = 7.5 event. Our analysis further reveals that the 1995 slip distribution was affected by a 1988 Mw = 7.2 event, and might have influenced a 1998 Mw = 7.0 earthquake that occurred downdip of the 1995 rupture. Our slip inversion reveals a potential change in mechanism in the southern portion of the rupture, consistent with seismic results. Predictions of the satellite LOS displacement from a seismic inversion and a joint seismic/GPS inversion do not compare favourably with the InSAR observations.

[1]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[2]  F. Webb,et al.  Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR C radar interferometry , 1996 .

[3]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[4]  C. Reigber,et al.  The Antofagasta 1995 Earthquake : Crustal deformation pattern as observed by GPS and D-INSAR , 1997 .

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

[6]  L. Rivera,et al.  The Andean subduction zone between 22 and 25°S (northern Chile): precise geometry and state of stress , 1996 .

[7]  E. Kissling,et al.  Accurate hypocentre determination in the seismogenic zone of the subducting Nazca Plate in northern Chile using a combined on-/offshore network , 1999 .

[8]  Thomas F. Coleman,et al.  A Reflective Newton Method for Minimizing a Quadratic Function Subject to Bounds on Some of the Variables , 1992, SIAM J. Optim..

[9]  H. Philip,et al.  Recent crustal deformation in the Antofagasta region (northern Chile) and the subduction process , 1998 .

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

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

[12]  J. Mechie,et al.  Two‐dimensional velocity models of the Nazca Plate subduction zone between 19.5°S and 25°S from wide‐angle seismic measurements during the CINCA95 project , 1999 .

[13]  H. Kanamori,et al.  The asperity model and the nature of large subduction zone earthquakes. , 1982 .

[14]  J. Montagner,et al.  Self-consistent retrieval of source parameters using mantle waves , 1998, Bulletin of the Seismological Society of America.

[15]  S. Barrientos,et al.  Coseismic coastal uplift and coralline algae record in Northern Chile: The 1995 Antofagasta earthquake case , 1996 .

[16]  T. Lay,et al.  Rupture Process of the 1995 Antofagasta Subduction Earthquake (Mw = 8.1) , 1999 .

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

[18]  J. Ruegg,et al.  Source tomography by simulated annealing using broad‐band surface waves and geodetic data: application to the Mw= 8.1 Chile 1995 event , 1997 .

[19]  P. Segall,et al.  Comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 Borah Peak earthquake , 1992, Bulletin of the Seismological Society of America.

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

[21]  Didier Massonnet,et al.  A comparison between short term (Co‐Seismic) and long term (one year) slip for the Landers Earthquake: Measurements from strong motion and SAR interferometry , 1997 .

[22]  L. Ruff,et al.  Seismic coupling along the Chilean Subduction Zone , 1991 .

[23]  K. Satake.  Depth distribution of coseismic slip along the Nankai Trough, Japan, from joint inversion of geodetic and tsunami data , 1993 .

[24]  P. Segall,et al.  Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas Fault , 1987 .

[25]  Michel Campillo,et al.  Contribution of radar interferometry to a two-step inversion of the kinematic process of the 1992 Landers earthquake , 1999 .

[26]  Kenneth W. Hudnut,et al.  Stress loading from viscous flow in the lower crust and triggering of aftershocks following the 1994 Northridge, California, Earthquake , 1999 .

[27]  J. C. Savage.  A dislocation model of strain accumulation and release at a subduction zone , 1983 .

[28]  M. Sobiesiak.  Fault plane structure of the Antofagasta, Chile Earthquake of 1995 , 2000 .

[29]  Richard G. Gordon,et al.  Effect of recent revisions to the geomagnetic reversal time scale on estimates of current plate motions , 1994 .

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

[31]  P. Heinrich,et al.  Numerical modeling of the 1995 Chilean Tsunami. Impact on French Polynesia , 1997 .

[32]  L. Rivera,et al.  Determination of seismogenic interplate contact zone and crustal seismicity around Antofagasta, northern Chile using local data , 1994 .

[33]  P. Rosen,et al.  On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake , 1994, Proceedings of IGARSS '94 - 1994 IEEE International Geoscience and Remote Sensing Symposium.

[34]  L. Rivera,et al.  The Mw = 8.0 Antofagasta (northern Chile) earthquake of 30 July 1995: A precursor to the end of the large 1877 gap , 1997, Bulletin of the Seismological Society of America.

[35]  Mario Pardo,et al.  Reappraisal of great historical earthquakes in the northern Chile and southern Peru seismic gaps , 1991 .

[36]  Paul A. Rosen,et al.  Crustal deformation measurements using repeat‐pass JERS 1 synthetic aperture radar interferometry near the Izu Peninsula, Japan , 1998 .

[37]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

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

[39]  Philip E. Gill,et al.  Practical optimization , 1981 .

[40]  Wayne Thatcher,et al.  Resolution of fault slip along the 470‐km‐long rupture of the great 1906 San Francisco earthquake and its implications , 1997 .

[41]  P. Visser,et al.  Precise orbit determination and gravity field improvement for the ERS satellites , 1998 .

[42]  W. Menke.  Geophysical data analysis , 1984 .

[43]  Ricardo Thiele,et al.  Active faulting in northern Chile: ramp stacking and lateral decoupling along a subduction plate boundary? , 1990 .

[44]  K. Koper,et al.  The waveguide effect of metastable olivine in slabs , 2000 .

[45]  Detlef Angermann,et al.  GPS-derived Deformation of the Central Andes Including the 1995 Antofagasta Mw = 8.0 Earthquake , 1999 .

[46]  R. Armijo,et al.  The MW=8.1 Antofagasta (North Chile) Earthquake of July 30, 1995: First results from teleseismic and geodetic data , 1996 .

[47]  P. Segall,et al.  Earthquake deformation cycle on the San Andreas Fault near Parkfield, California , 1987 .

[48]  P. Segall,et al.  The 1989 Loma Prieta earthquake imaged from inversion of geodetic data , 1994 .