An official website of the United States government
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
Heterogeneous Power‐Law Flow With Transient Creep in Southern California Following the 2010 El Mayor‐Cucapah Earthquake
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
- Award ID(s):
- 1848192
- 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
-
-
Abstract Following large earthquakes, viscoelastic stress relaxation may contribute to postseismic deformation observed at Earth's surface. Mechanical representations of viscoelastic deformation require a constitutive relationship for the lower crust/upper mantle material where stresses are diffused and, for non‐linear rheologies, knowledge of absolute stress level. Here, we describe a kinematic approach to representing geodetically observed postseismic motions that does not require an assumed viscoelastic rheology. The core idea is to use observed surface motions to constrain time‐dependent displacement boundary conditions applied at the base of the elastic upper crust by viscoelastic motions in the lower crust/upper mantle, approximating these displacements as slip on a set of dislocation elements. Using three‐dimensional forward models of viscoelastically modulated postseismic deformation in a thrust fault setting, we show how this approach can accurately represent surface motions and recover predicted displacements at the base of the elastic layer. Applied to the 1999 Chi‐Chi (Taiwan) earthquake, this kinematic approach can reproduce geodetically observed displacements and estimates of the partitioning between correlated postseismic deformation mechanisms. Specifically, we simultaneously estimate afterslip on the earthquake source fault that is similar to previous estimates, along with slip on dislocations at the base of the elastic layer that mimic predictions from viscous stress dissipation models in which viscosity is inferred to vary three‐dimensionally. A use case for the dislocation approach to modeling viscoelastic deformation is the estimation of spatiotemporally variable fault slip processes, including across sequential interseismic phases of the earthquake cycle, without assuming a lower crust/upper mantle 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.more » « less
-
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
-
Abstract The 2010MW7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events (M≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus onPnand its depth phasespPnandsPn, which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stackPnwaveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of thePn, directP, and directSphases recorded for event pairs. Focal depths of 93 out of 122M≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.more » « less
