Journal
EARTH PLANETS AND SPACE
Volume 71, Issue 1, Pages -Publisher
SPRINGEROPEN
DOI: 10.1186/s40623-019-1113-8
Keywords
Earthquake sequence simulation; Viscoelasticity; Spectral boundary integral equation method; Memory variable method
Categories
Funding
- MEXT KAKENHI [26109007]
- JSPS KAKENHI [19K04038]
- Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan
- Grants-in-Aid for Scientific Research [19K04038, 26109007] Funding Source: KAKEN
Ask authors/readers for more resources
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) 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) 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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available