Effect of a compliant fault zone on the inferred earthquake slip distribution

[1] We present a new semi-analytic method to evaluate the deformation due to a screw dislocation in arbitrarily heterogeneous and/or anisotropic elastic half plane. The method employs integral transformations to reduce the governing partial differential equations to the integral Fredholm equation of the second kind. Dislocation sources, as well as spatial perturbations in the elastic properties are modeled using equivalent body forces. The solution to the Fredholm equation is obtained in the Fourier domain using a method of successive over-relaxation, and is mapped into the spatial domain using the inverse Fast Fourier Transform. We apply this method to investigate the effect of a soft damage zone around an earthquake fault on the co-seismic displacement field, and on the earthquake slip distribution inferred from inversions of geodetic data. In the presence of a kilometer-wide damage zone with a reduction of the effective shear modulus of a factor of 2, inversions that assume a laterally homogeneous model tend to underestimate the amount of slip in the middle of the seismogenic layer by as much as 20%. This bias may accentuate the inferred maxima in the seismic moment release at depth between 3–6 km suggested by previous studies of large strike-slip earthquakes.

[1]  G. Mavko.  Mechanics of Motion on Major Faults , 1981 .

[2]  Leon Knopoff,et al.  Body Force Equivalents for Seismic Dislocations , 1964 .

[3]  Yehuda Ben-Zion,et al.  Pulverized rocks in the Mojave section of the San Andreas Fault Zone , 2006 .

[4]  N. N. Ambraseys,et al.  Some characteristic features of the Anatolian fault zone , 1970 .

[5]  K. Rybicki.  The elastic residual field of a very long strike-slip fault in the presence of a discontinuity , 1971, Bulletin of the Seismological Society of America.

[6]  Yuri Fialko,et al.  Evidence of fluid-filled upper crust from observations of postseismic deformation due to the 1992 Mw7.3 Landers earthquake , 2004 .

[7]  C. Thurber,et al.  Earthquake locations and three‐dimensional fault zone structure along the creeping section of the San Andreas fault near Parkfield, CA: Preparing for SAFOD , 2003 .

[8]  F. Chester,et al.  Fracture surface energy of the Punchbowl fault, San Andreas system , 2005, Nature.

[9]  R. Parker.  Geophysical Inverse Theory , 1994 .

[10]  K. Kasahara,et al.  A strike-slip fault in a laterally inhomogeneous medium , 1977 .

[11]  D. Barnett.  On the screw dislocation in an inhomogeneous elastic medium: the case of continuously varying elastic moduli , 1972 .

[12]  Huajian Gao,et al.  Dislocations in inhomogeneous media via a moduli perturbation approach: General formulation and two‐dimensional solutions , 1994 .

[13]  J. C. Savage.  Equivalent strike‐slip earthquake cycles in half‐space and lithosphere‐asthenosphere earth models , 1990 .

[14]  R. Lebensohn.  N-site modeling of a 3D viscoplastic polycrystal using Fast Fourier Transform , 2001 .

[15]  Carl W. Gable,et al.  Convection in three dimensions with surface plates: Generation of toroidal flow , 1991 .

[16]  C. M. Budwine,et al.  Institute of geophysics and planetary physics , 1991 .

[17]  Hervé Moulinec,et al.  A numerical method for computing the overall response of nonlinear composites with complex microstructure , 1998, ArXiv.

[18]  J. Brune,et al.  Particle size and energetics of gouge from earthquake rupture zones , 2005, Nature.

[19]  Yehuda Ben-Zion,et al.  Earthquake cycle, fault zones, and seismicity patterns in a rheologically layered lithosphere , 2001 .

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

[21]  James R. Rice,et al.  Off-Fault Secondary Failure Induced by a Dynamic Slip Pulse , 2005 .

[22]  David S. Watkins.  Iterative Methods for Linear Systems , 2005 .

[23]  Y. Fialko.  Fracture and Frictional Mechanics: Theory , 2007 .

[24]  K. W. Cattermole.  The Fourier Transform and its Applications , 1965 .

[25]  Julia R. Weertman,et al.  Elementary Dislocation Theory , 1992 .

[26]  D. Sandwell.  Thermomechanical evolution of oceanic fracture zones , 1984 .

[27]  D. Turcotte,et al.  Micro and macroscopic models of rock fracture , 2003 .

[28]  L. Delves,et al.  Computational methods for integral equations: Frontmatter , 1985 .

[29]  William H. K. Lee,et al.  Fine Structure of the Landers Fault Zone: Segmentation and the Rupture Process , 1994, Science.

[30]  C. Scholz,et al.  The process zone: A microstructural view of fault growth , 1998 .

[31]  Mark Simons,et al.  Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy , 2001 .

[32]  D. Lieberman,et al.  Particle size and energetics of gouge from earthquake rupture zones , 2022 .

[33]  J. Steketee,et al.  SOME GEOPHYSICAL APPLICATIONS OF THE ELASTICITY THEORY OF DISLOCATIONS , 1958 .

[34]  S. Nemat-Nasser,et al.  Micromechanics: Overall Properties of Heterogeneous Materials , 1993 .

[35]  J. D. Eshelby.  The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

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

[37]  J. C. Savage,et al.  Geodetic determination of relative plate motion in central California , 1973 .

[38]  V. Alshits,et al.  Elasticity of multilayers I. Basic equations and solutions , 1995 .

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

[40]  R. Bürgmann,et al.  The Effect of Elastic Layering on Inversions of GPS Data for Coseismic Slip and Resulting Stress Changes: Strike-Slip Earthquakes , 2005 .

[41]  Fred F. Pollitz,et al.  Gravitational viscoelastic postseismic relaxation on a layered spherical Earth , 1997 .

[42]  D. Sandwell,et al.  A three-dimensional semianalytic viscoelastic model for time-dependent analyses of the earthquake cycle , 2004 .

[43]  Rongjiang Wang,et al.  PSGRN/PSCMP - a new code for calculating co- and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory , 2006, Comput. Geosci..

[44]  V. Alshits,et al.  Elasticity of multilayers II. Strips, coatings and sandwiches , 1995 .

[45]  C. Doubre,et al.  Slip accumulation and lateral propagation of active normal faults in Afar , 2001 .

[46]  Bernard Minster,et al.  Deformation on Nearby Faults Induced by the 1999 Hector Mine Earthquake , 2002, Science.

[47]  Sylvain Barbot,et al.  Space geodetic investigation of the coseismic and postseismic deformation due to the 2003 Mw7.2 Altai earthquake: Implications for the local lithospheric rheology , 2008 .

[48]  D. Krajcinovic,et al.  Introduction to continuum damage mechanics , 1986 .

[49]  Yuri Fialko,et al.  Structure and mechanical properties of faults in the North Anatolian Fault system from InSAR observations of coseismic deformation due to the 1999 Izmit (Turkey) earthquake , 2007 .

[50]  Y. Fialko.  Probing the mechanical properties of seismically active crust with space geodesy: Study of the coseismic deformation due to the 1992 Mw7.3 Landers (southern California) earthquake , 2004 .

[51]  Steven G. Johnson,et al.  FFTW: an adaptive software architecture for the FFT , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[52]  L. E. Malvern.  Introduction to the mechanics of a continuous medium , 1969 .

[53]  Amos Nur,et al.  Postseismic Viscoelastic Rebound , 1974, Science.

[54]  P. Shearer,et al.  Seismic Imaging of the Damage Zone Around the Calico Fault , 2006 .