skip to main content


Title: Postseismic Deformation Following the 2015 M w 7.8 Gorkha (Nepal) Earthquake: New GPS Data, Kinematic and Dynamic Models, and the Roles of Afterslip and Viscoelastic Relaxation
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
Award ID(s):
1821853
NSF-PAR ID:
10453919
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
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
  1. 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
  2. Abstract

    Postseismic deformation following subduction earthquakes includes the combined effects of afterslip surrounding the coseismic rupture areas and viscoelastic relaxation in the asthenosphere and provides unique and valuable information for understanding the rheological structure. Because the two postseismic mechanisms are usually spatiotemporally intertwined, we developed an integrated model combining their contributions, based on 5 years of observations following the 2016 Pedernales (Ecuador) earthquake. The results show that the early, near‐field postseismic deformation is dominated by afterslip, both updip and downdip of the coseismic rupture, and requires heterogeneous interface frictional properties. Viscoelastic relaxation contributes more to far‐field displacements at later time periods. The best‐fit integrated model favors a 45‐km thick lithosphere overlying a Burgers body viscoelastic asthenosphere with a Maxwell viscosity of 3 × 1019 Pa s (0.9–5 × 1019 Pa s at 95% confidence), assuming the Kelvin viscosity equal to 10% of that value. In addition to the postseismic afterslip, the coastal displacement of sites north and south of the rupture clearly require extra slip in the plate motion direction due to slow slip events that may be triggered by the coseismic stress changes (CSC) but are not purely driven by the CSC. Spatially variable afterslip following the Pedernales event, combined with the SSEs during the interseismic period, demonstrate that spatial frictional variability persists throughout the whole earthquake cycle. The interaction of adjacent fault patches with heterogeneous properties may contribute to the clustered large earthquakes in this area.

     
    more » « less
  3. Abstract

    The 15 April 2015Mw7.8 Nepal Gorkha earthquake occurred on a shallowly dipping portion of the Main Himalayan Thrust (MHT). Notable features of the event include (1) the dominance of a slip pulse of about 6‐s duration that unlocked the lower edge of the MHT and (2) the near‐horizontal fault geometry, which, combined with proximity of the free surface, allows surface‐reflected phases to break the across‐fault symmetries of the seismic wavefield. Our dynamic rupture simulations in an elastoplastic medium yield earthquake parameters comparable to those deduced from kinematic inversions, including seismic moment and rupture velocity. The simulations reproduce pulse‐like behavior predicting pulse widths in agreement with those kinematic studies and supporting an interpretation in which the pulse‐like time dependence of slip is principally controlled by rupture geometry. This inference is strongly supported by comparison of synthetic ground velocity with the near‐field high‐rate GPS recording at station KKN4, which shows close agreement in pulse width, amplitude, and pulse shape. That comparison also constrains the updip extent of rupture and disfavors significant coseismic slip on the shallow ramp segment. Over most of the rupture length, the simulated rupture propagates at a near‐constant maximum velocity (~90% of theSwave speed) that is controlled by the antiplane geometry and off‐fault plastic yielding. Simulations also reveal the role of reflected seismic waves from the free surface, which may have contributed ~30% elongation of the slip pulse, and show the potential for significant free‐surface interaction effects in shallow events of similar geometry.

     
    more » « less
  4. Abstract

    The 2019 Ridgecrest conjugate Mw6.4 and Mw7.1 events resulted in several meters of strike‐slip and dip‐slip along an intricate rupture, extending from the surface down to 15 km. Now with >2 years of post‐rupture observations, we utilize these results to better understand vertical postseismic deformation from the Ridgecrest sequence and illuminate the emerging significance of vertical earthquake cycle deformation data. We determine the cumulative vertical displacement observed by the continuous GNSS network since Ridgecrest, which requires additional time series analyses to adequately resolve vertical deformation compared to the horizontal. Using a Maxwell‐type viscoelastic relaxation model, with a best fit time‐averaged asthenosphere viscosity of 4e17 Pa·s and a laterally heterogeneous lithosphere, we find that viscoelastic relaxation accounts for a majority of the cumulative vertical deformation at Ridgecrest and strongly controls far‐field observations in all north‐east‐up components. The viscoelastic model alone generally underpredicts deformation from GNSS and the remaining nonviscoelastic displacement is most prominent in the horizontal near‐field (−16 to 19 mm), revealing a deformation pattern matching the coseismic observations. This suggests that multiple deformation mechanisms are contributing to Ridgecrest's postseismic displacement, where afterslip likely dominates the near‐field while viscoelastic relaxation controls the far‐field. Similar deformation at individual GNSS stations has been observed for past earthquakes and additionally reveals long‐term transient viscosity over several years. Moreover, the greater temporal and spatial resolution of the GNSS array for Ridgecrest will help resolve the evolution of deformation for the entire network of observations as regional postseismic deformation persists for the next several years.

     
    more » « less
  5. Abstract

    Nearly every type of slip behavior of which megathrusts are known to be capable of has been observed within a narrow ~75‐km‐wide strip of the Sunda megathrust—the Banyak Islands subsegment. The diverse list of recorded slip events starts with the great 1861MW ~ 8.5 earthquake and includes the 1907MW ~ 8.2 shallow tsunami earthquake, the 1966–1981 slow slip event, the great 2005MW8.6 Nias‐Simeulue earthquake and associated afterslip, many small (M < 6) to moderate (6 M < 7) earthquakes, and the large 6 April 2010MW7.8 Banyak Islands earthquake. In this paper we map out the spatial and temporal relationships between these slip events to search for a pattern and possible controls on slip behavior. We use GPS and coral geodetic data to derive a coseismic slip distribution for the relatively understudied 2010 event and find that this event fits like a puzzle piece into a low slip patch left by the 2005 event. We compare the recent slip sequence with the historical sequence associated with the 1861 event and find that the modern slip pattern could be a repeat of the historical slip pattern, possibly suggesting the importance of stationary fault properties in controlling slip. The diverse set of slip events mostly fit our expectations within depth‐dependent rupture zones; however, the overlaps between events of different slip types may complicate this conceptual model. Along strike, we suggest that a subducting fracture zone could act to diversify the slip behavior, and we explore possible mechanisms for this slip control.

     
    more » « less