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.

Earthquake Cycles in a Model of Interacting Fault Patches: Complex Behavior at Transition from Seismic to Aseismic Slip

Complex Earthquake Cycle Simulations Using a Two-Degree-of-Freedom Spring-Block Model with a Rate- and State-Friction Law

The Positions, Colors, and Photometric Variability of Pluto's Small Satellites from HST Observations 2005-2006

Simple Spring-mass Model Simulation of Earthquake Cycle along the Nankai Trough in Southwest Japan

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

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

Characteristics of Earthquake Cycles: A Cross-Dimensional Comparison of 0D to 3D Numerical Models

Towards a simple, comprehensive model of regular earthquakes and slow slip events, part I: one-dimensional model

Phase response curves for models of earthquake fault dynamics

Dual effect of roughness during earthquake rupture sequences on faults with strongly rate-weakening friction

Progress in modeling the Tohoku-oki megathrust earthquake cycle and associated crustal deformation processes

Intrinsic Properties of a Burridge-Knopoff Model of an Earthquake Fault.

Toward quantitative characterization of simulated earthquake-cycle complexities

Different Earthquake Patterns for Two Neighboring Fault Segments Within the Haiyuan Fault Zone

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

Nature of the high-speed rupture of the two-dimensional Burridge-Knopoff model of earthquakes

Slow-Slip Phenomena Represented by the One-Dimensional Burridge-Knopoff Model of Earthquakes

Asperity Size and Neighboring Segments Can Change the Frictional Response and Fault Slip Behavior: Insights From Laboratory Experiments and Numerical Simulations

Kinematic Earthquake Sequences on Geometrically Complex Faults

Earthquake Nucleation and Slip Behavior Altered by Stochastic Normal Stress Heterogeneity

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