Earthquake cycle, fault zones, and seismicity patterns in a rheologically layered lithosphere

We study the coupled evolution of earthquakes and faults in a model consisting of a seismogenic upper crust governed by damage rheology over a viscoelastic substrate. The damage rheology has two types of functional coefficients: (1) a “generalized internal friction” separating states associated with material degradation and healing and (2) damage rate coefficients for positive (degradation) and negative (healing) changes. The evolving damage modifies the effective elastic properties of material in the upper crust as a function of the ongoing deformation. This simulates the creation and healing of fault systems in the upper seismogenic zone. In addition to the vertically averaged thin sheet approximation we introduce a Green function for three-dimensional elastic half-space for the instantaneous component of deformation. The formulation accounts in an internally consistent manner for evolving deformation fields, evolving fault structures, aseismic energy release, and spatiotemporal seismicity patterns. These developments allow us to simulate long histories of crustal deformation and to study the simultaneous evolution of regional earthquakes and faults for various model realizations. To focus on basic features of a large strike-slip fault system, we first consider a simplified geometry of the seismogenic crust by prescribing initial conditions consisting of a narrow damage zone in an otherwise damage-free plate. For this configuration, the model generates an earthquake cycle with distinct interseismic, preseismic, coseismic, and postseismic periods. Model evolution during each period is controlled by a subset of physical properties, which may be constrained by geophysical, geodetic, rock mechanics, and seismological data. In the more generic case with a random initial damage distribution, the model generates large crustal faults and subsidiary branches with complex geometries. The simulated statistics depend on the space-time window of the observational domain. The results indicate that long healing timescale, τh, describing systems with relatively long memory, leads to the development of geometrically regular fault systems and the characteristic frequency-size earthquake distribution. Conversely, short τh (relatively short memory) leads to the development of a network of disordered fault systems and the Gutenberg-Richter earthquake statistics. For intermediate values of τh the results exhibit alternating overall switching of response from periods of intense seismic activity and the characteristic earthquake distribution to periods of low seismic activity and Gutenberg-Richter statistics.

[1]  Yehuda Ben-Zion,et al.  Wrinkle-like slip pulse on a fault between different , 1997 .

[2]  S. Wdowinski,et al.  On the choice of boundary conditions in continuum models of continental deformation , 1990 .

[3]  Russell Robinson,et al.  Synthetic seismicity models of multiple interacting faults , 1995 .

[4]  R. Sibson.  Fault zone models, heat flow, and the depth distribution of earthquakes in the continental crust of the United States , 1982 .

[5]  Ari Ben-Menahem,et al.  Seismic waves and sources , 1981 .

[6]  J. Rice.  Spatio‐temporal complexity of slip on a fault , 1993 .

[7]  P. Okubo.  Dynamic rupture modeling with laboratory‐derived constitutive relations , 1989 .

[8]  G. King,et al.  Stress Triggering of the 1994 M = 6.7 Northridge, California, Earthquake by Its Predecessors , 1994, Science.

[9]  Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: a physical model of earthquake generation processes. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Y. Ben‐Zion.  Stress, slip, and earthquakes in models of complex single-fault systems incorporating brittle and creep deformations , 1996 .

[11]  Yehuda Ben-Zion,et al.  Transitions between Gutenberg-Richter and Characteristic Earthquake behavior in Simple Mean-Field Models of Heterogeneous Faults , 1998 .

[12]  Gabi Laske,et al.  CRUST 5.1: A global crustal model at 5° × 5° , 1998 .

[13]  R. Dorsey,et al.  EARTHQUAKE CLUSTERING INFERRED FROM PLIOCENE GILBERT-TYPE FAN DELTAS IN THE LORETO BASIN, BAJA CALIFORNIA SUR, MEXICO , 1997 .

[14]  J. C. Savage,et al.  The velocity field along the San Andreas Fault in central and southern California , 1991 .

[15]  H. Kanamori,et al.  State of stress before and after the 1994 Northridge Earthquake , 1997 .

[16]  Peggy A. Johnson,et al.  Rupture mechanism and interface separation in foam rubber models of earthquakes: a possible solution to the heat flow paradox and the paradox of large overthrusts , 1993 .

[17]  M. Heimpel,et al.  A seismodynamical model of lithosphere deformation: Development of continental and oceanic rift networks , 1996 .

[18]  M. A. Chinnery.  The deformation of the ground around surface faults , 1961 .

[19]  Wayne Thatcher,et al.  Nonlinear strain buildup and the earthquake cycle on the San Andreas Fault , 1983 .

[20]  James N. Brune,et al.  Wrinkle-like Weertman pulse at the interface between two blocks of foam rubber with different velocities , 1999 .

[21]  James C. Savage,et al.  Postseismic deformation associated with the 1992 M ω=7.3 Landers earthquake, southern California , 1997 .

[22]  Peter Mora,et al.  Simulation of the frictional stick-slip instability , 1994 .

[23]  Yehuda Ben-Zion,et al.  Earthquake Failure Sequences Along a Cellular Fault Zone in a Three-Dimensional Elastic Solid Containing Asperity and Nonasperity Regions , 1993 .

[24]  J. C. Savage,et al.  Is internal friction friction , 1996 .

[25]  D. J. Andrews,et al.  Coupling of energy between tectonic processes and earthquakes , 1978 .

[26]  D. J. Andrews,et al.  Mechanics of fault junctions , 1989 .

[27]  Z. Garfunkel,et al.  Plate kinematics of the circum Red Sea—a re-evaluation , 1987 .

[28]  C. Scholz,et al.  Fault growth and fault scaling laws: Preliminary results , 1993 .

[29]  Thomas H. Heaton,et al.  Static deformations from point forces and force couples located in welded elastic poissonian half-spaces: Implications for seismic moment tensors , 1989, Bulletin of the Seismological Society of America.

[30]  P. England,et al.  A thin viscous sheet model for continental deformation , 1982 .

[31]  Yehuda Bock,et al.  Current plate motion across the Dead Sea Fault from three years of continuous GPS monitoring , 2002 .

[32]  Thomas K. Rockwell,et al.  Paleoseismology of the Johnson Valley, Kickapoo, and Homestead Valley Faults: Clustering of Earthquakes in the Eastern California Shear Zone , 2000 .

[33]  Didier Sornette,et al.  Statistical physics of fault patterns self-organized by repeated earthquakes , 1994 .

[34]  Patience A. Cowie,et al.  A healing–reloading feedback control on the growth rate of seismogenic faults , 1998 .

[35]  D. J. Andrews,et al.  Properties and implications of dynamic rupture along a material interface , 1998, Bulletin of the Seismological Society of America.

[36]  Steven N. Ward,et al.  A synthetic seismicity model for southern California: Cycles, probabilities, and hazard , 1996 .

[37]  R. Weinberger,et al.  Non-linear elastic behaviour of damaged rocks , 1997 .

[38]  Y. Ben‐Zion,et al.  Distributed damage, faulting, and friction , 1997 .

[39]  D. Lockner,et al.  The role of microcracking in shear-fracture propagation in granite , 1995 .

[40]  B. Shibazaki,et al.  Spontaneous processes for nucleation, dynamic propagation, and stop of earthquake rupture , 1992 .

[41]  C. Frohlich,et al.  Teleseismic b values; Or, much ado about 1.0 , 1993 .

[42]  M. Zoback,et al.  Evidence for near-frictionless faulting in the 1989 (M 6.9) Loma Prieta, California, earthquake and its aftershocks , 1993 .

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

[44]  James H. Dieterich,et al.  Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering , 1997 .

[45]  F. Amelung,et al.  Earthquake scaling laws for creeping and non‐creeping faults , 1997 .

[46]  G. Schubert,et al.  Modeling of periodic great earthquakes on the San Andreas Fault: Effects of nonlinear crustal rheology , 1994 .

[47]  J. Dieterich.  Earthquake nucleation on faults with rate-and state-dependent strength , 1992 .

[48]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[49]  T. Rockwell,et al.  Crusader castle torn apart by earthquake at dawn, 20 May 1202 , 1998 .

[50]  J. Rice,et al.  Stress diffusion along rupturing plate boundaries , 1981 .

[51]  R. Simpson,et al.  In the shadow of 1857‐the effect of the Great Ft. Tejon Earthquake on subsequent earthquakes in southern California , 1996 .

[52]  Y. Ben‐Zion.  The response of two joined quarter spaces to SH line sources located at the material discontinuity interface , 1989 .

[53]  Yehuda Ben-Zion,et al.  Self-driven mode switching of earthquake activity on a fault system , 1999 .

[54]  R. E. Wallace.  Grouping and migration of surface faulting and variations in slip rates on faults in the Great Basin province , 1987 .

[55]  J. Brune,et al.  Heat Flow, Stress, and Rate of Slip along the San Andreas Fault, , 1969 .

[56]  J. Rice,et al.  Slip patterns and earthquake populations along different classes of faults in elastic solids , 1995 .

[57]  Steven G. Wesnousky,et al.  The Gutenberg-Richter or characteristic earthquake distribution, which is it? , 1994 .

[58]  Yehuda Ben-Zion,et al.  Interaction of the San Andreas Fault Creeping Segment with Adjacent great rupture zones and earthquake recurrence at Parkfield , 1993 .

[59]  J. Rice,et al.  Crustal deformation in great California earthquake cycles , 1986 .

[60]  Yehuda Ben-Zion,et al.  Dynamic simulations of slip on a smooth fault in an elastic solid , 1997 .

[61]  Geoffrey C. P. King,et al.  The accommodation of large strains in the upper lithosphere of the earth and other solids by self-similar fault systems: the geometrical origin of b-Value , 1983 .

[62]  R. Meissner,et al.  Limits of stresses in continental crusts and their relation to the depth-frequency distribution of shallow earthquakes , 1982 .

[63]  Kunihiko Shimazaki,et al.  FAULT TRACE COMPLEXITY, CUMULATIVE SLIP, AND THE SHAPE OF THE MAGNITUDE-FREQUENCY DISTRIBUTION FOR STRIKE-SLIP FAULTS : A GLOBAL SURVEY , 1996 .

[64]  T. Heaton.  Evidence for and implications of self-healing pulses of slip in earthquake rupture , 1990 .

[65]  Raul Madariaga,et al.  On the relation between seismic moment and stress drop in the presence of stress and strength heterogeneity , 1979 .

[66]  D. L. Anderson,et al.  Theoretical Basis of Some Empirical Relations in Seismology by Hiroo Kanamori And , 1975 .

[67]  M. Stein,et al.  Long-term earthquake clustering: A 50,000-year paleoseismic record in the Dead Sea Graben , 1996 .