Slow-slip, slow earthquakes, period-two cycles, full and partial ruptures, and deterministic chaos in a single asperity fault

• Published in: Tectonophysics Vol. 768; p. 228171
• Main Author: Barbot, Sylvain
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 05.10.2019; Elsevier BV
• Subjects: Asperity; Bifurcations; Chaos; Coefficient of friction; Computer simulation; Cycles; Dimensional analysis; Dynamics; Earthquake source; Earthquakes; Fault dynamics; Fault lines; Faults; Friction; Geological faults; Geophysics; Mathematical models; Modeling; Nucleation; Parameters; Physical properties; Physics; Rate-and-state friction; Seismic activity; Seismic cycle; Slip; Static friction; Upper bounds; Velocity; Fault dynamics; Seismic cycle; Rate-and-state friction; Modeling; Earthquake source
• ISSN: 0040-1951; 1879-3266
• DOI: 10.1016/j.tecto.2019.228171

The assimilation of geological, geophysical, and laboratory data in physics-based models of fault dynamics promises increasingly realistic simulations of the seismic cycle. To assist this effort, I explore the dynamics of a single velocity-weakening asperity to delineate the relationship between a fault physical properties and the style of faulting during the seismic cycle based on a micro-physical model of rate-and-state friction. From a dimensional analysis, four non-dimensional parameters control the characteristics of slip events and the behavior during the interseismic period, including the Dieterich-Ruina-Rice number Ru, which is proportional to the ratio of asperity to critical nucleation size W/h*, and Rb = (b − a)/b, where a and b control the direct and local evolutionary effects, respectively. For sufficiently large Rb, the Ru number controls the spectrum of fault slip, from periodic slow-slip events to aperiodic fast ruptures. In finite faults, but not in semi-infinite faults, the transition between these end-members includes a bifurcation that involves period-two, period-four, period-six cycles, and deterministic chaos with increasing Ru. Macroscopic fault slip is stabilized as Rb vanishes, even as Ru ≫ 1. For finite faults, but not for semi-infinite faults, the transition between slow and fast ruptures for increasing Rb includes another type of bifurcation with period-two and period-four cycles of slow and fast ruptures. For 0 <Rb ≪ 1 and Ru ≫ 1, slow ruptures become chaotic, characterized by aperiodic bursts of slow earthquakes within a longer slow-slip episode. The third non-dimensional parameter is a cut-off velocity that affects the maximum slip speed. The static friction coefficient μ0 strongly affects the rupture style, as fault strength provides an upper bound for stress drop. Many styles of instabilities and interseismic behaviors emerge depending on the coordinates in phase space, most importantly Ru and Rb. Combining complementary observations at different stages of the seismic cycle may offer an opportunity to infer these parameters in nature. •Four non-dimensional parameters control the dynamics of a single-asperity fault.•The Dieterich-Ruina-Rice number Ru controls the complexity of fast slip, with full and partial ruptures and aftershocks.•The Rb number measures the relative proximity to velocity-neutral. Slow slip with slow earthquakes can occur for low Rb numbers.•The transition from slow to fast ruptures forms a bifurcation in phase space involving period-two cycles and deterministic chaos.•Fault dynamics is more complex on weak faults (small coefficient of static friction), every other parameters being the same.