skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Contribution of Viscoelastic Stress to the Synchronization of Earthquake Cycles on Oceanic Transform Faults
Abstract 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
PAR ID:
10420683
Author(s) / Creator(s):
 ;  ;  
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
  1. Abstract 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,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
  2. Recent seismic and geodetic observations indicate that interseismic creep rate varies in both time and space. The spatial extent of creep pinpoints locked asperities, while its temporary accelerations, known as slow-slip events, may trigger earthquakes. Although the conditions promoting fault creep are well-studied, the mechanisms for initiating episodic slow-slip events are enigmatic. Here we investigate surface deformation measured by radar interferometry along the central San Andreas Fault between 2003 and 2010 to constrain the temporal evolution of creep. We show that slow-slip events are ensembles of localized creep bursts that aseismically rupture isolated fault compartments. Using a rate-and-state friction model, we show that effective normal stress is temporally variable on the fault, and support this using seismic observations. We propose that compaction-driven elevated pore fluid pressure in the hydraulically isolated fault zone and subsequent frictional dilation cause the observed slow-slip episodes. We further suggest that the 2004 Mw 6 Parkfield earthquake might have been triggered by a slow-slip event, which increased the Coulomb failure stress by up to 0.45 bar per year. This implies that while creeping segments are suggested to act as seismic rupture barriers, slow-slip events on these zones might promote seismicity on adjacent locked segments. 
    more » « less
  3. SUMMARY Observations of historical tsunami earthquakes reveal that ruptures of these earthquakes propagate slowly at shallow depth with longer duration, depletion in high-frequency radiation and larger discrepancy of Mw–Ms than ordinary megathrust earthquakes. They can effectively generate tsunami and lead to huge damage to regional populated areas near the coast. In this study, we use a recently developed dynamic earthquake simulator to explore tsunami earthquake generation from a physics-based modelling point of view. We build a shallow-dipping subduction zone model in which locally locked, unstable patches (asperities) are distributed on a conditionally stable subduction interface at shallow depth. The dynamic earthquake simulator captures both quasi-static and dynamic processes of earthquake cycles. We find that earthquakes can nucleate on these asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A high normal stress asperity, representing a subducted seamount, can act as an asperity in some events but as a barrier in other events over multiple earthquake cycles. Low normal stress asperities typically act as asperities in tsunami earthquakes. The degree of velocity-weakening in the conditionally stable zone, which may sustain rupture at different speeds or stop rupture, is critical for tsunami earthquake generation and affects its recurrence interval. Distributed asperities may rupture in isolated events separated by tens of years, or in a sequence of events separated by hours to days, or in one large event in a cascade fashion, demonstrating complex interactions among them. The recurrence interval on a high normal stress asperity is much larger than that on low normal stress asperities. These modelling results shed lights on the observations from historical tsunami earthquakes, including the 1994 and 2006 Java tsunami earthquakes and 2010 Mentawai tsunami earthquake. 
    more » « less
  4. 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
  5. Abstract The rheology of the crust and mantle and the interaction of viscoelastic flow with seismic/aseismic slip on faults control the state of stress in the lithosphere over multiple seismic cycles. The rheological behavior of rocks is well constrained in a laboratory setting, but thein situproperties of the lithosphere and its lateral variations remain poorly known. Here, we access the lower‐crustal rheology in Southern California by exploiting 8 years of geodetic postseismic deformation following the 2010 El Mayor‐Cucapah earthquake. The data illuminate viscoelastic flow in the lower crust with lateral variations of effective viscosity correlated with the geological province. We show that a Burgers assembly with dashpots following a nonlinear constitutive law can approximate the temporal evolution of stress and strain rate, indicating the activation of nonlinear transient creep before steady‐state dislocation creep. The transient and background viscosities in the lower crust of the Salton Trough are on the order of ~1018and ~1019 Pa s, respectively, about an order of magnitude lower than those in the surrounding regions. We highlight the importance of transient creep, nonlinear flow laws, and lateral variations of rheological properties to capture the entire history of postseismic relaxation following the El Mayor‐Cucapah earthquake. 
    more » « less