Synchronization behavior of large earthquakes (rupture of nearby faults close in time for many cycles) has been reported in many fault systems. The general idea is that the faults in the system have similar repeating intervals and are positively coupled through stress interaction. However, many details of such synchronization remain unknown. Here, we built a numerical model in the framework of rate‐and‐state friction to simulate earthquake cycles on the west Gofar fault, East Pacific Rise. Our model consists of two seismic patches separated by a barrier patch, which are constrained by seismic observations. We varied the parameters in the barrier to understand its role on earthquake synchronization. First, we found that when the barrier is relatively weak, synchronization can be achieved by afterslip or post‐seismic creep in the barrier patch. Second, static stress transfer can lead to synchronization, opposite to the suggestion by Scholz (2010,
Earthquake clustering can be promoted by local, regional, and remote triggering. The interaction between faults by static and dynamic stress transfer has received much attention. However, the role of quasi‐static stress interaction mediated by viscoelastic flow is still poorly understood. Here, we investigate whether the tight synchronization of moment‐magnitude 6 earthquakes every about 6 years on distant asperities in the Gofar‐Discovery fault system of the East Pacific Rise may be caused by mechanical coupling within the lithosphere‐asthenosphere system. We build a three‐dimensional numerical model of seismic cycles in the framework of rate‐ and state‐dependent friction with a brittle layer overlaying a viscoelastic mantle with nonlinear rheology to simulate earthquake cycles on separate asperities. The brittle section of the West Gofar fault consists of two frictionally unstable 20 km‐long by 5 km‐wide asperities separated by a velocity‐strengthening barrier, consistent with seismic observations, allowing stress transfer by afterslip and viscoelastic relaxation. We find that viscoelastic stress transfer can promote the synchronization of earthquakes. Even if the asperities are separated by as far as 30 km, synchronization is still possible for a viscosity of the underlying mantle of 1017 Pa s, which can be attained by dislocation creep or transient creep during the postseismic period. Considering the similarities in tectonic and structural settings, viscoelastic stress transfer and earthquake synchronization may also occur at 15’20 (Mid‐Atlantic Ridge), George V (Southeast Indian Ridge), Menard and Heezen transform fault (Pacific‐Antarctic Ridge).
more » « less- Award ID(s):
- 1848192
- NSF-PAR ID:
- 10420683
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 8
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract https://doi.org/10.1785/0120090309 ), which was based on results from a spring‐slider model using rate‐and‐state friction. Third, the width of the barrier is more important than its strength. When the barrier is narrow enough (no more than half the width of the seismic patch in our model), the system can achieve synchronization even with a very strong barrier. Fourth, for certain simulations, the interaction between the two seismic patches promotes partial rupture in the seismic patches and leads to complex behavior: the system switches from synchronized to unsynchronized over 10–20 cycles. -
more » « less
Oceanic transform faults display a unique combination of seismic and aseismic slip behavior, including a large globally averaged seismic deficit, and the local occurrence of repeating magnitude (M)
earthquakes with abundant foreshocks and seismic swarms, as on the Gofar transform of the East Pacific Rise and the Blanco Ridge in the northeast Pacific Ocean. However, the underlying mechanisms that govern the partitioning between seismic and aseismic slip and their interaction remain unclear. Here we present a numerical modeling study of earthquake sequences and aseismic transient slip on oceanic transform faults. In the model, strong dilatancy strengthening, supported by seismic imaging that indicates enhanced fluid-filled porosity and possible hydrothermal circulation down to the brittle–ductile transition, effectively stabilizes along-strike seismic rupture propagation and results in rupture barriers where aseismic transients arise episodically. The modeled slow slip migrates along the barrier zones at speeds ∼10 to 600 m/h, spatiotemporally correlated with the observed migration of seismic swarms on the Gofar transform. Our model thus suggests the possible prevalence of episodic aseismic transients in M rupture barrier zones that host active swarms on oceanic transform faults and provides candidates for future seafloor geodesy experiments to verify the relation between aseismic fault slip, earthquake swarms, and fault zone hydromechanical properties. -
more » « less
Abstract The constitutive behavior of faults is central to many interconnected aspects of earthquake science, from fault dynamics to induced seismicity, to seismic hazards characterization. Yet, a friction law applicable to the range of temperatures found in the brittle crust and upper mantle is still missing. In particular, rocks often exhibit a transition from steady‐state velocity‐strengthening at room temperature to velocity‐weakening in warmer conditions that is poorly understood. Here, we investigate the effect of competing healing mechanisms on the evolution of frictional resistance in a physical model of rate‐, state‐, and temperature‐dependent friction. The yield strength for fault slip depends on the real area of contact, which is modulated by the competition between the growth and erosion of interfacial micro‐asperities. Incorporating multiple healing mechanisms and rock‐forming minerals with different thermodynamic properties allows a transition of the velocity‐ and temperature‐dependence of friction at steady‐state with varying temperatures. We explain the mechanical data for granite, pyroxene, amphibole, shale, and natural fault gouges with activation energies and stress power exponent for weakening of 10–50 kJ/mol and 55–150, respectively, compatible with subcritical crack growth and inter‐granular flow in the active slip zone. Activation energies for the time‐dependent healing process in the range 90–130 kJ/mol in dry conditions and 20–65 kJ/mol in wet conditions indicate the prominence of viscoelastic collapse of microasperities in the absence of water and of pressure‐solution creep, crack healing, and cementation when assisted by pore fluids.
-
more » « less
Abstract Foreshocks are not uncommon prior to large earthquakes, but their physical mechanism remains controversial. Two interpretations have been advanced: (1) foreshocks are driven by aseismic nucleation and (2) foreshocks are cascades, with each event triggered by earlier ones. Here, we study seismic cycles on faults with fractal roughness at wavelengths exceeding the nucleation length. We perform 2‐D quasi‐dynamic, elastic simulations of frictionally uniform rate‐state faults. Roughness leads to a range of slip behavior between system‐size ruptures, including widespread creep, localized slow slip, and microseismicity. These processes are explained by spatial variations in normal stress (
σ ) caused by roughness: regions with lowσ tend to creep, while highσ regions remain locked until they break seismically. Foreshocks and mainshocks both initiate from the rupture of locked asperities, but mainshocks preferentially start on stronger asperities. The preseismic phase is characterized by feedback between creep and foreshocks: episodic seismic bursts break groups of nearby asperities, causing creep to accelerate, which in turns loads other asperities leading to further foreshocks. A simple analytical treatment of this mutual stress transfer, confirmed by simulations, predicts slip velocities and seismicity rates increase as 1/t , wheret is the time to the mainshock. The model reproduces the observed migration of foreshocks toward the mainshock hypocenter, foreshock locations consistent with static stress changes, and the 1/t acceleration in stacked catalogs. Instead of interpreting foreshocks as either driven by coseismic stress changes or by creep, we propose that earthquake nucleation on rough faults is driven by the feedback between the two. -
more » « less
Abstract Small repeating earthquakes are thought to represent rupture of isolated asperities loaded by surrounding creep. The observed scaling between recurrence interval and seismic moment,
T r ∼M 1/6, contrasts with expectation assuming constant stress drop and no aseismic slip (T r ∼M 1/3). Here we demonstrate that simple crack models of velocity‐weakening asperities in a velocity‐strengthening fault predict theM 1/6scaling; however, the mechanism depends on asperity radius,R . For small asperities (, where is the nucleation radius) numerical simulations with rate‐state friction show interseismic creep penetrating inward from the edge, and earthquakes nucleate in the center and rupture the entire asperity. Creep penetration accounts for ∼25% of the slip budget, the nucleation phase takes up a larger fraction of slip. Stress drop increases with increasing R ; the lack of self‐similarity being due to the finite nucleation dimension. Forsimulations exhibit simple cycles with ruptures nucleating from the edge. Asperities with exhibit complex cycles of partial and full ruptures. Here T r is explained by an energy criterion: full rupture requires that the energy release rate everywhere on the asperity at least equals the fracture energy, leading to the scalingT r ∼M 1/6. Remarkably, in spite of the variability in behavior with source dimension, the scaling ofT r with stress drop Δτ , nucleation length and creep ratev plis the same across all regimes:. This supports the use of repeating earthquakes as creepmeters and provides a physical interpretation for the scaling observed in nature.
