Abstract:The time-dependent viscoelastic deformation of host rocks is important when considering the dynamics of fault behavior, specifically in brittle–ductile transitional regions or shallow subduction zones, because it relaxes stress heterogeneity and affects loading to the fault. For a rate-and-state fault embedded in a Maxwell viscoelastic medium, a previous study discovered a transition from repeated earthquakes to the permanent stuck of a rate-weakening patch (EQ–ST transition) with decreased viscoelastic relaxation time . This transition differs from the well-known seismic–aseismic transition explained by the Hopf bifurcation at the critical stiffness of an elastic system. To better understand the EQ–ST transition, quantifying the effect of heterogeneous frictional strength is important, because this effect is characteristic to the viscoelastic medium and is absent in the elastic limit. Previous experimental studies suggest a potential contrast in frictional strength in such a way that a rate-weakening patch is stronger than a rate-strengthening region containing clay minerals. Here, we conducted two-dimensional, fully dynamic earthquake sequence simulations for a fault in a Maxwell viscoelastic medium; we investigated the EQ–ST transition in the two-dimensional parameter space of and . With , the EQ–ST transition occurred at about 1 order of magnitude larger than in the case with . We constructed a coarse-grained model with only two degrees of freedom based on the spatial average. Consequently, the coarse-grained model behaves remarkably similar to the continuum model, and the EQ–ST transition is associated with a homoclinic bifurcation. The EQ–ST boundary in the parameter space can be quantitatively explained by considering elastic loading due to creep in the rate-strengthening region and unloading by viscoelastic relaxation of the stress heterogeneity comparable to . A larger rate-weakening patch is anticipated to become aseismic earlier as decreases, because the elastic loading rate is inversely correlated with the patch size. This may be qualitatively consistent with the change in the size distribution of events around the brittle–ductile transition in observations and laboratory experiments; however, further investigations on, for example, the interactions of events and changes in frictional parameters with depth are required for quantitative discussion.

Steady‐State Effective Normal Stress in Subduction Zones Based on Hydraulic Models and Implications for Shallow Slow Earthquakes

Two sides of a fault: Grain-scale analysis of pore pressure control on fault slip

A weak subducting slab at intermediate depths below northeast Japan

Fault Valving and Pore Pressure Evolution in Simulations of Earthquake Sequences and Aseismic Slip

3-D thermal regime and dehydration processes around the regions of slow earthquakes along the Ryukyu Trench

2-D elastic FEM simulation on stress state in the deep part of a subducted slab

Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

Effective Bulk Rheology of a Two‐Phase Subduction Shear Zone: Insights From Micromechanics‐Based Modeling and Implications for Subduction Interface Slow Slip Events

Analysis of Near-Field Stresses in an Analogue Strike-Slip Fault Model

Stress change in southwest Japan due to the 1944–1946 Nankai megathrust rupture sequence based on a 3-D heterogeneous rheological model

Earthquake Nucleation and Slip Behavior Altered by Stochastic Normal Stress Heterogeneity

Homoclinic bifurcation of a rate-weakening patch in a viscoelastic medium and effect of strength contrast

Path and Slip Dependent Behavior of Shallow Subduction Shear Zones During Fluid Overpressure

Stress‐driven Relaxation of Heterogeneous Upper Mantle and Time‐dependent Afterslip Following the 2011 Tohoku Earthquake

Seismic Versus Aseismic Slip: Probing Mechanical Properties of the Northeast Japan Subduction Zone

Influence of Farallon Slab Loading on Intraplate Stress and Seismicity in Eastern North America in the Presence of Pre‐Existing Weak Zones

Viscoelastic simulation of deep tectonic deformation of the Longmenshan fault zone and its implication for strong earthquakes

Frictional instabilities in clay illuminate the origin of slow earthquakes

Spatiotemporal Evolution of Pore Pressure Changes and Coulomb Failure Stress in a Poroelastic Medium for Different Faulting Regimes

The Effect of Normal Stress Oscillations on Fault Slip Behavior Near the Stability Transition From Stable to Unstable Motion

Field‐based description of near‐surface crustal deformation in a high‐strain shear zone: A case study in southern Kyushu, Japan