Dynamics of earthquake nucleation process represented by the Burridge-Knopoff model

Dynamics of earthquake nucleation process is studied on the basis of the one-dimensional Burridge-Knopoff (BK) model obeying the rate- and state-dependent friction (RSF) law. We investigate the properties of the model at each stage of the nucleation process, including the quasi-static initial phase, the unstable acceleration phase and the high-speed rupture phase or a mainshock. Two kinds of nucleation lengths L_sc and L_c are identified and investigated. The nucleation length L_sc and the initial phase exist only for a weak frictional instability regime, while the nucleation length L_c and the acceleration phase exist for both weak and strong instability regimes. Both L_sc and L_c are found to be determined by the model parameters, the frictional weakening parameter and the elastic stiffness parameter, hardly dependent on the size of an ensuing mainshock. The sliding velocity is extremely slow in the initial phase up to L_sc, of order the pulling speed of the plate, while it reaches a detectable level at a certain stage of the acceleration phase. The continuum limits of the results are discussed. The continuum limit of the BK model lies in the weak frictional instability regime so that a mature homogeneous fault under the RSF law always accompanies the quasi-static nucleation process. Duration times of each stage of the nucleation process are examined. The relation to the elastic continuum model and implications to real seismicity are discussed.Comment: Title changed. Changes mainly in abstract and in section 1. To appear in European Physical Journal

The Effects of Fault Roughness on the Earthquake Nucleation Process

We study numerically the effects of fault roughness on the nucleation process during earthquake sequences. The faults are governed by a rate and state friction law. The roughness introduces local barriers that complicate the nucleation process and result in asymmetric expansion of the rupture, nonmonotonic increase in the slip rates on the fault, and the generation of multiple slip pulses. These complexities are reflected as irregular fluctuations in the moment rate. There is a large difference between first slip events in the sequences and later events. In the first events, for roughness amplitude b_r ≤ 0.002, there is a large increase in the nucleation length with increasing br. For larger values of b_r, slip is mostly aseismic. For the later events there is a trade-off between the effects of the finite fault length and the fault roughness. For b_r ≤ 0.002, the finite length is a more dominant factor and the nucleation length barely changes with br. For larger values of b_r, the roughness plays a larger role and the nucleation length increases significantly with b_r. Using an energy balance approach, where the roughness is accounted for in the fault stiffness, we derive an approximate solution for the nucleation length on rough faults. The solution agrees well with the main trends observed in the simulations for the later events and provides an estimate of the frictional and roughness properties under which faults experience a transition between seismic and aseismic slip

The Effects of Fault Roughness on the Earthquake Nucleation Process

We study numerically the effects of fault roughness on the nucleation process during earthquake sequences. The faults are governed by a rate and state friction law. The roughness introduces local barriers that complicate the nucleation process and result in asymmetric expansion of the rupture, nonmonotonic increase in the slip rates on the fault, and the generation of multiple slip pulses. These complexities are reflected as irregular fluctuations in the moment rate. There is a large difference between first slip events in the sequences and later events. In the first events, for roughness amplitude b_r ≤ 0.002, there is a large increase in the nucleation length with increasing br. For larger values of b_r, slip is mostly aseismic. For the later events there is a trade-off between the effects of the finite fault length and the fault roughness. For b_r ≤ 0.002, the finite length is a more dominant factor and the nucleation length barely changes with br. For larger values of b_r, the roughness plays a larger role and the nucleation length increases significantly with b_r. Using an energy balance approach, where the roughness is accounted for in the fault stiffness, we derive an approximate solution for the nucleation length on rough faults. The solution agrees well with the main trends observed in the simulations for the later events and provides an estimate of the frictional and roughness properties under which faults experience a transition between seismic and aseismic slip

Nucleation process in the Burridge-Knopoff model of earthquakes

Nucleation process of the one-dimensional Burridge-Knopoff model of earthquakes obeying the rate- and state-dependent friction law is studied both analytically and numerically. The properties of the nucleation dynamics, the nucleation lengths and the duration times are examined together with their continuum limits.Comment: Title changed, one figure (previous Fig.2) omitted, several references (new Refs.5-8,21,24-28) added, and two eqs. (new eqs.1,5) added. Text expanded considerably, especially the part explaining the relation of the model to the elastic continuum model. To appear in Europhys. Letter

Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate- and state-dependent friction

We present an efficient and rigorous numerical procedure for calculating the elastodynamic response of a fault subjected to slow tectonic loading processes of long duration within which there are episodes of rapid earthquake failure. This is done for a general class of rate- and state-dependent friction laws with positive direct velocity effect. The algorithm allows us to treat accurately, within a single computational procedure, loading intervals of thousands of years and to calculate, for each earthquake episode, initially aseismic accelerating slip prior to dynamic rupture, the rupture propagation itself, rapid post seismic deformation which follows, and also ongoing creep slippage throughout the loading period in velocity-strengthening fault regions. The methodology is presented using the two-dimensional (2-D) antiplane spectral formulation and can be readily extended to the 2-D in-plane and 3-D spectral formulations and, with certain modifications, to the space-time boundary integral formulations as well as to their discretized development using finite difference or finite element methods. The methodology can be used to address a number of important issues, such as fault operation under low overall stress, interaction of dynamic rupture propagation with pore pressure development, patterns of rupture propagation in events nucleated naturally as a part of a sequence, the earthquake nucleation process, earthquake sequences on faults with heterogeneous frictional properties and/or normal stress, and others. The procedure is illustrated for a 2-D crustal strike-slip fault model with depth-variable properties. For lower values of the state-evolution distance of the friction law, small events appear. The nucleation phases of the small and large events are very similar, suggesting that the size of an event is determined by the conditions on the fault segments the event is propagating into rather than by the nucleation process itself. We demonstrate the importance of incorporating slow tectonic loading with elastodynamics by evaluating two simplified approaches, one with the slow tectonic loading but no wave effects and the other with all dynamic effects included but much higher loading rate

Foreshock properties illuminate nucleation processes of slow and fast laboratory earthquakes

Laboratory experiments demonstrate that prior to fast laboratory earthquakes the fault begins to unlock and creep, causing foreshocks to coalesce in both space and time. This demonstrates that the evolution of foreshocks is closely connected to the fault slip velocity.Understanding the connection between seismic activity and the earthquake nucleation process is a fundamental goal in earthquake seismology with important implications for earthquake early warning systems and forecasting. We use high-resolution acoustic emission (AE) waveform measurements from laboratory stick-slip experiments that span a spectrum of slow to fast slip rates to probe spatiotemporal properties of laboratory foreshocks and nucleation processes. We measure waveform similarity and pairwise differential travel-times (DTT) between AEs throughout the seismic cycle. AEs broadcasted prior to slow labquakes have small DTT and high waveform similarity relative to fast labquakes. We show that during slow stick-slip, the fault never fully locks, and waveform similarity and pairwise differential travel times do not evolve throughout the seismic cycle. In contrast, fast laboratory earthquakes are preceded by a rapid increase in waveform similarity late in the seismic cycle and a reduction in differential travel times, indicating that AEs begin to coalesce as the fault slip velocity increases leading up to failure. These observations point to key differences in the nucleation process of slow and fast labquakes and suggest that the spatiotemporal evolution of laboratory foreshocks is linked to fault slip velocity

Modeling afterslip and aftershocks following the 1992 Landers earthquake

One way to probe the rheology of the lithosphere and fault zones is to analyze the temporal evolution of deformation following a large earthquake. In such a case, the lithosphere responds to a known stress change that can be assessed from earthquake slip models constrained from seismology and geodesy. Here, we model the postseismic response of a fault zone that is assumed to obey a rate-strengthening rheology, where the frictional stress varies as aσ ln(ε), ε being the deformation rate and aσ > 0 a rheological parameter. The model is simple enough that these parameters can be estimated by inversion of postseismic geodetic data. We apply this approach to the analysis of geodetic displacements following the M_w 7.3, 1992, Landers earthquake. The model adjusts well the measured displacements and implies aσ ≈ 0.47–0.53 MPa. In addition, we show that aftershocks and afterslip follow the same temporal evolution and that the spatiotemporal distribution of aftershocks is consistent with the idea that they are driven by reloading of the seismogenic zone resulting from frictional afterslip

Theory of time-dependent rupture in the Earth

Fracture mechanics is used to develop a theory of earthquake mechanism which includes the phenomenon of subcritical crack growth. The following phenomena are predicted: slow earthquakes, multiple events, delayed multiple events (doublets), postseismic rupture growth and afterslip, foreshocks, and aftershocks. The theory predicts a nucleation stage prior to an earthquake, and suggests a physical mechanism by which one earthquake may 'trigger' another

The SCEC/USGS Dynamic Earthquake Rupture Code Verification Exercise

Numerical simulations of earthquake rupture dynamics are now common, yet it has been difficult to test the validity of thesesimulations because there have been few field observations and no analytic solutions with which to compare the results. This paper describes the Southern California Earthquake Center/U.S. Geological Surve(SCEC/USGS) Dynamic Earthquake Rupture Code Verification Exercise, where codes that simulate spontaneous rupture dynamics in three dimensions are evaluated and the results produced by these codes are compared using Web-based tools. This is the first time that a broad and rigorous examination of numerous spontaneous rupture codes has been performed—a significant advance in this science. The automated process developed to attain this achievement provides for a future where testing of codes is easily accomplished. Scientists who use computer simulations to understand earthquakes utilize a range of techniques. Most of these assume that earthquakes are caused by slip at depth on faults in the Earth, but hereafter the strategies vary. Among the methods used in earthquake mechanics studies are kinematic approaches and dynamic approaches. The kinematic approach uses a computer code that prescribes the spatial and temporal evolution of slip on the causative fault (or faults). These types of simulations are very helpful, especially since they can be used in seismic data inversions to relate the ground motions recorded in the field to slip on the fault(s) at depth. However, these kinematic solutions generally provide no insight into the physics driving the fault slip or information about why the involved fault(s) slipped that much (or that little). In other words, these kinematic solutions may lack information about the physical dynamics of earthquake rupture that will be most helpful in forecasting future events. To help address this issue, some researchers use computer codes to numerically simulate earthquakes and construct dynamic, spontaneous rupture (hereafter called “spontaneous rupture”) solutions. For these types of numerical simulations, rather than prescribing the slip function at each location on the fault(s), just the friction constitutive properties and initial stress conditions are prescribed. The subsequent stresses and fault slip spontaneously evolve over time as part of the elasto-dynamic solution. Therefore, spontaneous rupture computer simulations of earthquakes allow us to include everything that we know, or think that we know, about earthquake dynamics and to test these ideas against earthquake observations