Fault Parameters of the Nisqually Earthquake Determined from Moment Tensor Solutions and the Surface Deformation from GPS and InSAR

The magnitude 6.8 Nisqually earthquake occurred on 28 February 2001 at a depth of 50-60 km within the subducting Juan de Fuca plate. The fault parameters are estimated from moment tensor solutions and by comparing the surface displacements from Global Positioning System (GPS) data and from satellite interferometry, with predictions from elastic deformation models. The simple deformation model calculates the coseismic surface displacements caused by an earthquake on an extensional rectangular fault in a uniform half-space. Continuous GPS stations within 200 km of the epicenter resolved horizontal displacements as great as 9 mm (210°) and vertical displacements as great as 13 mm of subsidence near the epicenter. The near-vertical displacements were also determined from a differential interferogram created from synthetic aperture radar (InSAR) data from RADARSAT , the only satellite data available for the event. A maximum vertical subsidence of approximately 20 mm is observed 30 km east of the epicenter. The GPS, InSAR, and moment tensor solutions provide a consistent solution for the rupture parameters of the Nisqually earthquake. The fault has a strike of 180° and a dip of 20° and is centered at 47.10° N and 122.67° W (4 km east and 6 km south of the epicenter) or has a strike of 360° and a dip of 70° and is centered at 47.10° N and 122.69° W (2 km east and 6 km south). The trade-off between fault area and rupture displacement is not resolved by our data, but a good fit is found with the main rupture having an area of 230 km2 and an along-strike length of 23 km and downdip width of 10 km. The rupture area is centered at a depth of 51 km with a scalar moment of 2.0 × 1019 N m and, with the given area, an average rupture displacement of 1.4 m. To refine the interpretation of the GPS and InSAR data, a 3D heterogeneous numerical model was generated having realistic shear moduli structure, a 3D model of the subducting slab, and a spherical Earth. The results are similar to those from the half-space model, but there is significant refinement with a deeper rupture center at 60 km. Manuscript received 9 April 2003.

[1]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

[2]  T. Vanwagoner,et al.  Three-dimensional reflection image of the subducting Juan de Fuca plate , 2002 .

[3]  Kelin Wang,et al.  The Cascadia Subduction Zone and related subduction systems: seismic structure, intraslab earthquakes and processes, and earthquake hazards , 2002 .

[4]  Y. Okada.  Surface deformation due to shear and tensile faults in a half-space , 1985 .

[5]  Kelin Wang,et al.  Three‐dimensional dislocation model for great earthquakes of the Cascadia Subduction Zone , 1997 .

[6]  Stuart A. Spikin.  Estimation of earthquake source parameters by the inversion of waveform data: Global seismicity, 1981-1983 , 1986 .

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

[8]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

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

[10]  Howard A. Zebker,et al.  Decorrelation in interferometric radar echoes , 1992, IEEE Trans. Geosci. Remote. Sens..

[11]  Charles Werner,et al.  Accuracy of topographic maps derived from ERS-1 interferometric radar , 1994, IEEE Trans. Geosci. Remote. Sens..

[12]  G. Rogers,et al.  Moment Magnitude–Local Magnitude Calibration for Earthquakes off Canada's West Coast , 2003 .

[13]  Didier Massonnet,et al.  SATELLITE RADAR INTERFEROMETRY , 1997 .

[14]  Simon M. Peacock,et al.  Subduction factory 2. Are intermediate‐depth earthquakes in subducting slabs linked to metamorphic dehydration reactions? , 2003 .

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

[16]  Moment-tensor analysis using regional data: Application to the 25 March, 1993, Scotts Mills, Oregon, earthquake , 1995 .

[17]  Simon M. Peacock,et al.  Subduction factory 1. Theoretical mineralogy, densities, seismic wave speeds, and H 2 O contents , 2003 .

[18]  Charles A. Langston,et al.  Source inversion of seismic waveforms: The Koyna, India, earthquakes of 13 September 1967 , 1981 .

[19]  J. Nábělek,et al.  Seismotectonics of the Explorer region , 2002 .

[20]  Michael R. Craymer,et al.  Current tectonics of northern Cascadia from a decade of GPS measurements , 2003 .

[21]  K. Wang,et al.  Thermal structure and metamorphism of subducting oceanic crust: insight into Cascadia intraslab earthquakes , 2002 .