An official website of the United States government
Abstract We report Global Positioning System (GPS) measurements of postseismic deformation following the 2015 Mw7.8 Gorkha (Nepal) earthquake, including previously unpublished data from 13 continuous GPS stations installed in southern Tibet shortly after the earthquake. We use variational Bayesian Independent Component Analysis (vbICA) to extract the signal of postseismic deformation from the GPS time series, revealing a broad displacement field extending >150 km northward from the rupture. Kinematic inversions and dynamic forward models show that these displacements could have been produced solely by afterslip on the Main Himalayan Thrust (MHT) but would require a broad distribution of afterslip extending similarly far north. This would require the constitutive parameter(a − b)σto decrease northward on the MHT to ≤0.05 MPa (an extreme sensitivity of creep rate to stress change) and seems unlikely in light of the low interseismic coupling and high midcrustal temperatures beneath southern Tibet. We conclude that the northward reach of postseismic deformation more likely results from distributed viscoelastic relaxation, possibly in a midcrustal shear zone extending northward from the seismogenic MHT. Assuming a shear zone 5–20 km thick, we estimate an effective shear‐zone viscosity of ~3·1016–3·1017 Pa·s over the first 1.12 postseismic years. Near‐field deformation can be more plausibly explained by afterslip itself and implies(a − b)σ ~ 0.5–1 MPa, consistent with other afterslip studies. This near‐field afterslip by itself would have re‐increased the Coulomb stress by ≥0.05 MPa over >30% of the Gorkha rupture zone in the first postseismic year, and deformation further north would have compounded this reloading.
more »
« less
Exploring GPS Observations of Postseismic Deformation Following the 2012 M W 7.8 Haida Gwaii and 2013 M W 7.5 Craig, Alaska Earthquakes: Implications for Viscoelastic Earth Structure
Abstract The Queen Charlotte‐Fairweather Fault (QC‐FF) system off the coast of British Columbia and southeast Alaska is a highly active dextral strike‐slip plate boundary that accommodates ∼50 mm/yr of relative motion between the Pacific and North America plates. NineMW ≥ 6.7 earthquakes have occurred along the QC‐FF system since 1910, including aMS(G‐R)8.1 event in 1949. Two recent earthquakes, the October 28, 2012 Haida Gwaii (MW7.8) and January 5, 2013 Craig, Alaska (MW7.5) events, produced postseismic transient deformation that was recorded in the motions of 25 nearby continuous Global Positioning System (cGPS) stations. Here, we use 5+ yr of cGPS measurements to characterize the underlying mechanisms of postseismic deformation and to constrain the viscosity structure of the upper mantle surrounding the QC‐FF. We construct forward models of viscoelastic deformation driven by coseismic stress changes from these two earthquakes and explore a large set of laterally heterogeneous viscosity structures that incorporate a relatively weak back‐arc domain; we then evaluate each model based on its fit to the postseismic signals in our cGPS data. In determining best‐fit model structures, we additionally incorporate the effects of afterslip following the 2012 event. Our results indicate the occurrence of a combination of temporally decaying afterslip and vigorous viscoelastic relaxation of the mantle asthenosphere. In addition, our best‐fit viscosity structure (transient viscosity of 1.4–2.0 × 1018 Pa s; steady‐state viscosity of 1019 Pa s) is consistent with the range of upper mantle viscosities determined in previous studies of glacial isostatic rebound and postseismic deformation.
more »
« less
- Award ID(s):
- 1802364
- PAR ID:
- 10446151
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 7
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract On 18 November 2022, a large earthquake struck offshore southern Sumatra, generating a tsunami with 25 cm peak amplitude recorded at tide gauge station SBLT. OurW‐phase solution indicates a shallow dip of 6.2°, compatible with long‐period surface wave radiation patterns. Inversion of teleseismic body waves indicates a shallow slip distribution extending from about 10 km deep to near the trench with maximum slip of ∼4.1 m and seismic moment of Nm (MW7.3). Joint modeling of seismic and tsunami data indicates a shallow rigidity of ∼23 GPa. We find a low moment‐scaled radiated energy of , similar to that of the 2010MW7.8 Mentawai event () and other tsunami earthquakes. These characteristics indicate that the 2022 event should be designated as a smaller moment magnitude tsunami earthquake compared to the other 12 well‐documented global occurrences since 1896. The 2022 event ruptured up‐dip of the 2007MW8.4 Bengkulu earthquake, demonstrating shallow seismogenic capability of a megathrust that had experienced both a deeper seismic event and adjacent shallow aseismic afterslip. We consider seismogenic behavior of shallow megathrusts and concern for future tsunami earthquakes in subduction zones globally, noting a correlation between tsunami earthquake occurrence and subducting seafloor covered with siliceous pelagic sediments. We suggest that the combination of pelagic clay and siliceous sediments and rough seafloor topography near the trench play important roles in controlling the genesis of tsunami earthquakes along Sumatra and other regions, rather than the subduction tectonic framework of accretionary or erosive margin.more » « less
-
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 Stress‐based postseismic deformation modeling including afterslip and viscoelastic relaxation usually assumes the coseismic slip distribution and the associated stress perturbation as known. However, that assumption biases the postseismic modeling results by the assumptions that underlie the coseismic models. Importantly, this misses an opportunity to iteratively constrain the coseismic slip model with postseismic observations. We used a broad set of seismic and geodetic data to create multiple coseismic slip models that only differ in the down‐dip extent of the rupture plane and fit the coseismic observations for the July 29, Mw 8.2 Chignik earthquake equally well. We then evaluated the quality of those coseismic slip models based on how well each of them predicts postseismic GNSS displacements using a stress‐driven afterslip model. We find that coseismic slip models that generate afterslip too far down‐dip systematically fail to predict postseismic deformation. We find that the postseismic observations are best predicted by a narrower coseismic slip model that terminates abruptly at its deepest extent. The model predictions improve further if stress‐driven afterslip is combined with a superimposed viscoelastic relaxation response of a 50 km thick elastic lithosphere for the overriding plate and an elastic cold nose to the mantle wedge.more » « less
-
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
