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).
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 the
- Award ID(s):
- 1848192
- NSF-PAR ID:
- 10380383
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 9
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We use Sentinel‐1 and ALOS‐2 Interferometric Synthetic Aperture Radar (InSAR), and Global Navigation Satellite System (GNSS) data to investigate the mechanisms of coseismic and postseismic deformation due to the 2021 M7.4 Maduo (China) earthquake. We present a refined coseismic slip model constrained by the rupture trace and precisely located aftershocks. The InSAR time series corrected for the atmospheric and decorrelation noise reveal postseismic line of sight displacements up to ∼0.1 m. The displacements are discontinuous along the fault trace, indicating shallow afterslip and velocity‐strengthening friction in the top 2–3 km of the upper crust. The magnitude of shallow afterslip is however insufficient to compensate for the coseismic slip deficit, implying substantial off‐fault yielding. The observed surface deformation does not exhibit obvious features that could be attributed to poroelastic effects. We developed a fully coupled model that accounts for both stress‐driven creep on a deep localized shear zone and viscoelastic relaxation in the bulk of the lower crust. The mid‐ to near‐field data can be reasonably well explained by deep afterslip and/or non‐Maxwellian visco‐elasticity. Our results suggest a power‐law stress exponent of ∼4–4.5 assuming a power‐law rheology, and transient and steady‐state viscosities of 1018and 1019 Pa s, respectively, assuming a bi‐viscous (Burgers) rheology. However, a good fit to the GNSS data cannot be achieved assuming the bulk viscoelastic relaxation alone, and requires a contribution of deep afterlip and/or a localized shear zone extending through much of the lower crust.
-
more » « less
Abstract Rheologic variations in the Earth's crust (like elastic plate thickness [EPT] or crustal rigidity) modulate the rate at which seismic moment accumulates for potentially hazardous faults of the San Andreas Fault System (SAFS). To quantify rates of seismic moment accumulation, Global Navigation Satellite Systems, and Interferometric Synthetic Aperture Radar data were used to constrain surface deformation rates of a four‐dimensional viscoelastic deformation model that incorporates rheological variations spanning a 900 km section of the SAFS. Lateral variations in EPT, estimated from surface heat flow and seismic depth to the lithosphere‐asthenosphere boundary, were converted to lateral variations in rigidity and then used to solve for seismic moment accumulation rates on 32 fault segments. We find a cluster of elevated seismic moment rates (11–20 × 1015 Nm year−1km−1) along the main SAFS trace spanning the historical
M w 7.9 1857 Fort Tejon earthquake rupture length; present‐day seismic moment magnitude on these segments ranges fromM w 7.2–7.6. We also find that the average plate thickness in the Salton Trough is reduced to only 60% of the regional average, which results in a ∼60% decrease in moment accumulation rate along the Imperial fault. Likewise, a 30% increase of average plate thickness results in at least a ∼30% increase in moment rate and even larger increases are identified in regions of complex plate heterogeneity. These results emphasize the importance of considering rheological variations when estimating seismic hazard, suggesting that meaningful changes in seismic moment accumulation are revealed when considering spatial variations in crustal rheology. -
more » « less
Abstract Flexure occurs on intermediate geologic timescales (∼1 Myr) due to volcanic‐island building at the Island of Hawaii, and the deformational response of the lithosphere is simultaneously elastic, plastic, and ductile. At shallow depths and low temperatures, elastic deformation transitions to frictional failure on faults where stresses exceed a threshold value, and this complex rheology controls the rate of deformation manifested by earthquakes. In this study, we estimate the seismic strain rate based on earthquakes recorded between 1960 and 2019 at Hawaii, and the estimated strain rate with 10−18–10−15s−1in magnitude exhibits a local minimum or neutral bending plane at 15 km depth within the lithosphere. In comparison, flexure and internal deformation of the lithosphere are modeled in 3D viscoelastic loading models where deformation at shallow depths is accommodated by frictional sliding on faults and limited by the frictional coefficient (
μ f), and at larger depths by low‐temperature plasticity and high‐temperature creep. Observations of flexure and the seismic strain rate are best‐reproduced by models withμ f = 0.3 ± 0.1 and modified laboratory‐derived low‐temperature plasticity. Results also suggest strong lateral variations in the frictional strength of faults beneath Hawaii. Our models predict a radial pattern of compressive stress axes relative to central Hawaii consistent with observations of earthquake pressure (P) axes. We demonstrate that the dip angle of this radial axis is essential to discerning a change in the curvature of flexure, and therefore has implications for constraining lateral variations in lithospheric strength. -
Abstract The rheology of the upper mantle impacts a variety of geodynamic processes, including postseismic deformation following great earthquakes and post-glacial rebound. The deformation of upper mantle rocks is controlled by the rheology of olivine, the most abundant upper mantle mineral. The mechanical properties of olivine at steady state are well constrained. However, the physical mechanism underlying transient creep, an evolutionary, hardening phase converging to steady state asymptotically, is still poorly understood. Here, we constrain a constitutive framework that captures transient creep and steady state creep consistently using the mechanical data from laboratory experiments on natural dunites containing at least 94% olivine under both hydrous and anhydrous conditions. The constitutive framework represents a Burgers assembly with a thermally activated nonlinear stress-versus-strain-rate relationship for the dashpots. Work hardening is obtained by the evolution of a state variable that represents internal stress. We determine the flow law parameters for dunites using a Markov chain Monte Carlo method. We find the activation energy $$430\pm 20$$ 430 ± 20 and $$250\pm 10$$ 250 ± 10 kJ/mol for dry and wet conditions, respectively, and the stress exponent $$2.0\pm 0.1$$ 2.0 ± 0.1 for both the dry and wet cases for transient creep, consistently lower than those of steady-state creep, suggesting a separate physical mechanism. For wet dunites in the grain-boundary sliding regime, the grain-size dependence is similar for transient creep and steady-state creep. The lower activation energy of transient creep could be due to a higher jog density of the corresponding soft-slip system. More experimental data are required to estimate the activation volume and water content exponent of transient creep. The constitutive relation used and its associated flow law parameters provide useful constraints for geodynamics applications. Graphical Abstractmore » « less
