Fully dynamic earthquake sequence simulation of a fault in a viscoelastic medium using a spectral boundary integral equation method: does interseismic stress relaxation promote aseismic transients?

Abstract

Abstract Along subduction interfaces or some major faults, a seismogenic layer in the upper crust is underlain by a zone of slow-slip events (SSEs) and tremors, and seismicity disappears at greater depths. The transition between seismogenic and aseismic behavior may be caused by changes in the frictional properties of the fault or changes in viscoelastic properties of the surrounding medium. Although aseismic transients have been numerically generated in previous studies by changing the frictional properties and compared to SSEs, the effect of viscoelasticity on the transition remains to be studied. In this study, we implemented interseismic viscoelastic stress relaxation in a simulation code for two-dimensional antiplane fully dynamic earthquake sequences in a uniform elastic material based on a spectral boundary integral equation method. In the implementation, we developed a suitable algorithm in which the viscoelastic relaxation is calculated by evolution of an “effective slip,” which gives viscoelastic stress change on the fault if convolved with a static Green’s function. We conducted parametric studies for a fault with a rate-weakening patch on two parameters: viscoelastic relaxation time $$t_{c}$$ tc and characteristic length of the state evolution in the rate- and state-dependent friction law L. The behavior of the simulated fault can be classified into four classes, earthquakes (EQ), aseismic transients (AT), stable sliding (SS), and stuck (ST), in which the central part of the rate-weakening patch has a diminishingly small slip rate and is permanently locked. A phase diagram of the fault behavior shows that there are two different types of seismogenic–aseismic transition. As L increases, an EQ patch changes to an AT patch before becoming an SS patch, as has been reported in previous studies in an elastic limit. The boundary between AT and SS can be explained via linear stability analysis of a system composed of a spring, a dashpot, and sliders. As $$t_{c}$$ tc decreases, the recurrence interval of the earthquakes diverges, and an EQ patch changes to an ST patch unless L is within a narrow range. Therefore, the transition associated with SSEs and tremors is dominated by changes in frictional properties rather than changes in viscoelastic properties.

Keywords

• Earthquake sequence simulation
• Viscoelasticity
• Spectral boundary integral equation method
• Memory variable method