On the mechanics of earthquake afterslip

We propose a model for earthquake afterslip based on rate and state variable friction laws. In the model, afterslip is attributed to the interaction of a velocity-weakening region at depth (within which earthquakes nucleate) with an upper region of velocity-strengthening frictional behavior. The existence of this upper region is supported by independent seismologic observations and the results of laboratory friction experiments. In our model, afterslip is the result of relaxation of a stress perturbation within the velocity-strengthening region, which arises when an earthquake propagates into that region from below. We derive the stress perturbation and its decay from the friction constitutive law using a simple, 1 degree-of-freedom approximation for the elastic interaction between the fault and its surroundings. This approximation is based on thickness-averaged displacements and slip velocities within the velocity-strengthening region, which is assumed to slip as a rigid block. Coseismic and postseismic slip are coupled through the thickness-averaged stiffness k of the velocity-strengthening region. We assume k to be inversely proportional to the thickness of this region, which means that thicker velocity strengthening regions have a greater tendency to arrest coseismic slip. We model the afterslip-time histories of the 1966 Parkfield and 1987 Superstition Hills earthquakes and relate the model parameters to physical parameters which may govern the rheologic behavior of the faults. In accord with field observations, our model predicts (1) that afterslip on some faults scales with the thickness of the (unconsolidated) sedimentary cover and (2) that proportionally more afterslip occurs for earthquakes in which coseismic surface slip is small compared with coseismic slip at depth. Velocity-strengthening frictional behavior is to be expected for faults within poorly consolidated sediments and for those that contain significant gouge zones (about >500 m) within their shallow regions (<3–5 km). Combining our results with those of recent laboratory friction studies indicates that relatively young faults with little accumulated fault gouge should exhibit little afterslip.

[1]  Robert V. Sharp,et al.  Surface faulting in the central Imperial Valley: in The Imperial Valley, California Earthquake of October 15, 1979 , 1982 .

[2]  H. Kanamori,et al.  Depth of seismicity in the Imperial Valley Region (1977–1983) and its relationship to heat flow, crustal structure and the October 15, 1979, earthquake , 1986 .

[3]  A. Ruina,et al.  Slip motion and stability of a single degree of freedom elastic system with rate and state dependent friction , 1984 .

[4]  R. Wesson Dynamics of fault creep , 1988 .

[5]  C. Crook Geodetic measurements of horizontal crustal deformation associated with the October 15, 1979, Imperial Valley (California) earthquake , 1984 .

[6]  Ralph J. Archuleta,et al.  A faulting model for the 1979 Imperial Valley earthquake , 1984 .

[7]  R. V. Sharp,et al.  Three-dimensional records of surface displacement on the Superstition Hills fault zone associated with the earthquakes of 24 November 1987 , 1989 .

[8]  W. D. Stuart,et al.  Long‐term fault creep observations in central California , 1982 .

[9]  C. Scholz,et al.  Dynamic properties of stick‐slip friction of rock , 1976 .

[10]  James H. Dieterich,et al.  Effects of physical fault properties on frictional instabilities produced on simulated faults , 1984 .

[11]  M. Wyss,et al.  Seismic and aseismic slip on the San Andreas Fault , 1969 .

[12]  J. Gu,et al.  Earthquake aftereffects and triggered seismic phenomena , 1983 .

[13]  H. Quin Dynamic stress drop and rupture dynamics of the October 15, 1979 Imperial Valley, California, earthquake , 1990 .

[14]  R. V. Sharp,et al.  Fault movement (afterslip) following the Guatemala earthquake of February 4, 1976 , 1978 .

[15]  J. Weeks,et al.  Constitutive behavior and stability of frictional sliding of granite , 1986 .

[16]  R. Wesson Modelling aftershock migration and afterslip of the San Juan Bautista, California, earthquake of October 3, 1972 , 1987 .

[17]  W. Brace,et al.  California Earthquakes: Why Only Shallow Focus? , 1970, Science.

[18]  Patrick L. Williams,et al.  Slip Along the Superstition Hills Fault Associated with the 24 November 1987 Superstition Hills, California, Earthquake , 1989 .

[19]  M. Wyss,et al.  DISPLACEMENT ON THE SAN ANDREAS FAULT SUBSEQUENT TO THE 1966 PARKFIELD EARTHQUAKE , 1968 .

[20]  Zhao Yusheng,et al.  Effects of load point velocity on frictional instability behavior , 1990 .

[21]  J. Weertman,et al.  A dislocation theory analysis of fault creep events , 1973 .

[22]  R. Bilham Surface slip subsequent to the 24 November 1987 Superstition Hills, California, earthquake monitored by digital creepmeters , 1989 .

[23]  Chris Marone,et al.  Frictional behavior and constitutive modeling of simulated fault gouge , 1990 .

[24]  Stephen R. Brown,et al.  Roughness of natural fault surfaces , 1987 .

[25]  John R. Rice,et al.  Crustal Earthquake Instability in Relation to the Depth Variation of Frictional Slip Properties , 1986 .

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

[27]  Chris Marone,et al.  The depth of seismic faulting and the upper transition from stable to unstable slip regimes , 1988 .

[28]  R. V. Sharp,et al.  Surface faulting near Lake Oroville, California in August, 1975 , 1976 .

[29]  C. H. Scholz,et al.  The critical slip distance for seismic faulting , 1988, Nature.

[30]  George A. McMechan,et al.  Crustal structure of the Imperial Valley region. , 1982 .

[31]  C. R. Allen,et al.  Fault slip in southern California , 1985 .

[32]  Duncan Carr Agnew,et al.  The 1987 Superstition Hills earthquake sequence: Strains and tilts at Piñon Flat Observatory , 1989 .

[33]  C. Scholz The Mechanics of Earthquakes and Faulting , 1990 .

[34]  William D. Stuart,et al.  Forecast model for great earthquakes at the Nankai Trough subduction zone , 1988 .

[35]  James H. Dieterich,et al.  The frictional properties of a simulated gouge having a fractal particle distribution , 1989 .

[36]  J. Dieterich Modeling of rock friction: 1. Experimental results and constitutive equations , 1979 .

[37]  Robert V. Sharp,et al.  Surface faulting along the Superstition Hills fault zone and nearby faults associated with the earthquakes of 24 November 1987 , 1989 .

[38]  J. Rice,et al.  Stressing of locked patches along a creeping fault , 1985 .

[39]  David J. Wald,et al.  Rupture process of the 1987 Superstition Hills earthquake from the inversion of strong-motion data , 1990, Bulletin of the Seismological Society of America.

[40]  A. Ruina Slip instability and state variable friction laws , 1983 .

[41]  T. Tullis,et al.  Geodetic Predictions of a strike-slip fault model: Implications for intermediate- and short-term earthquake prediction , 1989 .