Nonuniform and Unsteady Sliding of a Plate Boundary in a Great Earthquake Cycle: A Numerical Simulation Using a Laboratory-derived Friction Law

A numerical study is conducted to simulate complicated sliding behavior and earthquake activity on a subducting plate boundary. A 2-D model of a uniform elastic half-space with a semi-infinite thrust fault is set up, and the frictional stress prescribed by a rate- and state-dependent friction law is assumed to act on the plate boundary fault. Spatial nonuniformity of friction parameters representing rate-dependence of friction and of slip-dependence of friction are introduced in the model to obtain complicated sliding behavior in the numerical simulation. Analogs of great earthquakes that break the entire seismogenic plate boundary repeatedly occur at a constant time interval. Smaller events of seismic or aseismic sliding occur during a great earthquake cycle. Regions of rate-strengthening of friction and of a large characteristic distance in slip-dependence of friction behave as barriers or asperities. Rupture propagation is often arrested in such a region and a great earthquake occurs later when the region is broken. The variety of earthquake activity observed in many regions along real plate boundaries may be explained by similar nonuniformity in friction parameters. Conversely, the friction parameters on plate boundaries might be estimated from comparison of theoretical simulations with observations of earthquake activity. Simulation results indicate that spatiotemporal variation in stress due to aseismic sliding may play an important part in generating earthquakes.