Interactions between Temperature and Pore-Fluid Pressure during Earthquake Faulting and a Mechanism for Partial or Total Stress Relief

Sibson, Richard H.

doi:10.1038/physci243066a0

Cited by 401 publications

(389 citation statements)

References 7 publications

Supporting

Mentioning

380

Contrasting

Unclassified

Order By: Relevance

“…Shimamoto (2003, 2005) and Spray (2005) observed slip weakening during high slip rate sliding when melting occurred. Thermal pressurization was proposed by Sibson (1973) and modeled by Lachenbruch (1980), Mase and Smith (1987), Andrews (2002), Wibberley and Shimamoto (2005), Bizzarri and Cocco (2006) and Rice (2006) and many others. It is difficult to confirm in the laboratory whether thermal pressurization occurs or not because no footprints of thermal pressurization remain after high-speed slip; all the above studies are based on theoretical considerations.…”

Section: Slip-weakening Distancementioning

confidence: 99%

Constitutive parameters for earthquake rupture dynamics based on high-velocity friction tests with variable sliprate

Fukuyama

Mizoguchi

2009

Int J Fract

View full text Add to dashboard Cite

To reproduce the dynamic rupture process of earthquakes, the fault geometry, initial stress distribution and a frictional constitutive law on the fault are important parameters as initial and boundary conditions of the system. Here, we focus on the frictional constitutive relation on the fault. During a high-speed rupture, fault strength decreases as slip develops which can be described by a slip weakening equation. To understand the physical process of stress breakdown during the dynamic rupture of earthquakes, we investigated the friction behavior of rocks in the laboratory by direct measurements of traction evolution with slip in response to a given slip history. We employed a highspeed rotary shear apparatus introduced at National Research Institute for Earth Science and Disaster Prevention (NIED). This apparatus has a capability of sliding with predefined variable velocities using a servo-controlled system. We used a pair of granite cylindrical specimens with a diameter of 25 mm. As an input signal, we used a regularized Yoffe function to investigate the scale dependence of fracture energy and slip weakening distance (D c ). We observed a positive correlation between D c and total slip, keeping the maximum slip velocity constant. These conditions correspond to those for earthquakes with the same stress drop and varying magnitudes. Finally, we used a

show abstract

Section: Slip-weakening Distancementioning

confidence: 99%

Constitutive parameters for earthquake rupture dynamics based on high-velocity friction tests with variable sliprate

Fukuyama

Mizoguchi

2009

Int J Fract

View full text Add to dashboard Cite

show abstract

“…Thus 7 controls the relative amount of slipweakening versus velocity-weakening effects [Shaw, 1995]. This formulation is based on a physical idea which goes back to Sibson [1973], whereby frictional heating increases pore fluid temperature, and therefore pressure, thereby reducing the effective normal stress. The slow heat dissipation slipweakening limit was solved by Lachenbruch [1980] and later discussO by Rice [1994].…”

Section: The Frictionmentioning

confidence: 99%

Existence of continuum complexity in the elastodynamics of repeated fault ruptures

Shaw¹,

Rice²

2000

J. Geophys. Res.

View full text Add to dashboard Cite

Abstract. What are the origins of earthquake complexity? The possibility that some aspects of the complexity displayed by earthquakes might be explained by stress heterogeneities developed through the self-organization of repeated ruptures has been suggested by some simple self-organizing models. The question of whether or not even these simple self-organizing models require at least some degree of material heterogeneity to maintain complex sequences of events has been the subject of some controversy. In one class of elastodynamic models, previous work has described complexity as arising on a model fault with completely uniform material properties. Questions'were raised, however, regarding the role of discreteness, the relevance of the nucleation mechanism, and special parameter choices, in generating the complexity that has been reported. In this paper, we examine the question of whether or not continuum complexity is achieved under the stringent conditions of continuous loading, and whether the results are similar to previously claimed findings of continuum complexity or its absence. The elastodynamic model that we use consists of a 1-D fault boundary with friction, a steady slowly moving 1-D boundary parallel to the fault, and a 2-D scalar elastic media connecting the two boundaries. The constitutive law used involves a pair of sequential weakening processes, one occurring over a small slip (or velocity) and accomplishing a small fraction of the total strength drop, and the other at larger slip (or velocity) and providing the remaining strength drop. The large-scale process is motivated by a heat weakening instability. Our main results are as follows. (1) We generally find complexity of type I, a broad distribution of large event sizes with nonperiodic recurrence, when the modeled region is very long, along strike, compared to the layer thickness. (2) We find that complexity of type II, with numerous small events showing a power law distribution, is not a generic result but does definitely exist in a restricted range of parameter space. For that, in the slip weakening version of our model, the strength drop and nucleation size in the small slip process must be much smaller than in the large slip process, and the nucleation length associated with the latter must be comparable to layer thickness. This suggests a basis for reconciling different previously reported results. (3) Bulk dispersion appears to be relatively unimportant to the results. In particular, motions on the fault plane are seen to be relatively insensitive to a wide range of changes in the dispersion in the bulk away from the fault, both at long wavelengths and at short wavelengths. In contrast, the fault properties are seen to be very important to the results. (4) Nucleation from slip weakening and time-dependent weakening showed similar large-scale behavior. However, not all constitutive laws are insensitive to all nucleation approximations; those making a model "inherently discrete" and hence grid-dependent, in particular, can affect large sc...

show abstract

“…When fluids are present and temperature changes in faults, thermal pressurization will yield resistance on the fault plane and thus play a significant role in earthquake rupture (Sibson, 1973;Lachenbruch, 1980;Chester and Higgs, 1992;Fialko, 2004;Fialko and Khzan, 2005;Bizzari and Cocco, 2006a, b;Rice, 2006;Wang, 2000Wang, , 2006Wang, , 2009Wang, , 2011Wang, , 2013Wang, , 2016bWang, , 2017Bizzarri, 2010Bizzarri, , 2011a. Rice (2006) proposed two end-member models for thermal pressurization: the adiabatic-undrained-deformation (AUD) model and the slip-on-a-plane (SOP) model.…”

Section: Friction Caused By Thermal Pressurizationmentioning

confidence: 99%

Multistable slip of a one-degree-of-freedom spring-slider model in the presence of thermal-pressurized slip-weakening friction and viscosity

Wang

2017

Nonlin. Processes Geophys.

View full text Add to dashboard Cite

Abstract. This study is focused on multistable slip of earthquakes based on a one-degree-of-freedom spring-slider model in the presence of thermal-pressurized slip-weakening friction and viscosity by using the normalized equation of motion of the model. The major model parameters are the normalized characteristic displacement, U c , of the friction law and the normalized viscosity coefficient, η, between the slider and background plate. Analytic results at small slip suggest that there is a solution regime for η and γ (= 1/U c ) to make the slider slip steadily. Numerical simulations exhibit that the time variation in normalized velocity, V /V max (V max is the maximum velocity), obviously depends on U c and η. The effect on the amplitude is stronger due to η than due to U c . In the phase portrait of V /V max versus the normalized displacement, U/U max (U max is the maximum displacement), there are two fixed points. The one at large V /V max and large U/U max is not an attractor, while that at small V /V max and small U/U max can be an attractor for some values of η and U c . When U c <0.55, unstable slip does not exist. When U c ≥ 0.55, U c and η divide the solution domain into three regimes: stable, intermittent, and unstable (or chaotic) regimes. For a certain U c , the three regimes are controlled by a lower bound, η l , and an upper bound, η u , of η. The values of η l , η u , and η u − η l all decrease with increasing U c , thus suggesting that it is easier to yield unstable slip for larger U c than for smaller U c or for larger η than for smaller η. When U c <1, the Fourier spectra calculated from simulation velocity waveforms exhibit several peaks, thus suggesting the existence of nonlinear behavior of the system. When U c >1, the related Fourier spectra show only one peak, thus suggesting linear behavior of the system.

show abstract

Interactions between Temperature and Pore-Fluid Pressure during Earthquake Faulting and a Mechanism for Partial or Total Stress Relief

Cited by 401 publications

References 7 publications

Constitutive parameters for earthquake rupture dynamics based on high-velocity friction tests with variable sliprate

Constitutive parameters for earthquake rupture dynamics based on high-velocity friction tests with variable sliprate

Existence of continuum complexity in the elastodynamics of repeated fault ruptures

Multistable slip of a one-degree-of-freedom spring-slider model in the presence of thermal-pressurized slip-weakening friction and viscosity

Contact Info

Product

Resources

About