Subduction zones host some of Earth's most damaging natural hazards, including megathrust earthquakes and earthquake‐induced tsunamis. A major control on the initiation and rupture characteristics of subduction megathrust earthquakes is how the coupled zone along the subduction interface accumulates elastic strain between events. We present results from observations of slow slip events (SSEs) in Cascadia occurring during the interseismic period downdip of the fully coupled zone, which imply that the orientation of strain accumulation within the coupled zone can vary with depth. Interseismic GPS motions suggest that forces derived from relative plate motions across a shallow, offshore locked plate interface dominate over decadal timescales. Deeper on the plate interface, below the locked (seismogenic) patch, slip during SSEs dominantly occurs in the updip direction, reflecting a dip‐parallel force acting on the slab, such as slab pull. This implies that in subduction zones with obliquely convergent plate motions, the seismogenic zone of the megathrust is loaded by forces acting in two discrete directions, leading to a depth‐varying orientation of strain accumulation on the plate interface.
- Award ID(s):
- 1821853
- PAR ID:
- 10097700
- Date Published:
- Journal Name:
- Pure and Applied Geophysics
- ISSN:
- 0033-4553
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We use densely spaced campaign GPS observations and laboratory friction experiments on fault rocks from one of the world's most rapidly slipping low‐angle normal faults, the Mai'iu fault in Papua New Guinea, to investigate the nature of interseismic deformation on active low‐angle normal faults. GPS velocities reveal 8.3 ± 1.2 mm/year of horizontal extension across the Mai'iu fault, and are fit well by dislocation models with shallow fault locking (above 2 km depth), or by deeper locking (from ~5–16 km depth) together with shallower creep. Laboratory friction experiments show that gouges from the shallowest portion of the fault zone are predominantly weak and velocity‐strengthening, while fault rocks deformed at greater depths are stronger and velocity‐weakening. Evaluating the geodetic and friction results together with geophysical and microstructural evidence for mixed‐mode seismic and aseismic slip at depth, we find that the Mai'iu fault is most likely strongly locked at depths of ~5–16 km and creeping updip and downdip of this region. Our results suggest that the Mai'iu fault and other active low‐angle normal faults can slip in large (Mw > 7) earthquakes despite near‐surface interseismic creep on frictionally stable clay‐rich gouges.
-
more » « less
Abstract The Cascadia subduction zone (CSZ) is known to host M9 megathrust ruptures; however, no such event has occurred in historical observation. The distribution and characteristics of small‐ to moderate‐sized earthquakes can be used to determine the behavior of the megathrust fault but are notably absent offshore the CSZ due to the distance from onshore seismometers. We use automated subspace detection coupled with an onshore‐offshore seismic deployment to find small‐magnitude earthquakes in the offshore seismogenic zone and analyze their locations in the context of interseismic locking and seismogenic zone extent. We detected and located 5,282 earthquakes, 4,096 of which had been previously undetected. We find that the downdip extent of the seismogenic zone as defined by interplate seismicity agrees with the 20% locking contour of the Schmalzle et al. (2014,
https://doi.org/10.1002/2013GC005172 ) geodetic model and extends deeper than predicted by previous thermal models. We cannot determine the updip extent of the seismogenic zone; this may be due to a lack of templates for detection in the updip source area, stress shadows updip of asperity loading, and/or strong locking to the trench. We present a map of possible asperities determined by the small earthquakes in this study. Our asperity locations and extents show some, but not complete, agreement with the asperities modeled from the 1700 M9 rupture and geodetic locking models, and good agreement with the paleo‐rupture extents determined from offshore turbidites and forearc basin‐based asperity estimates. This highlights the need of continued offshore observations over time, and to elucidate fine‐scale variation in locking. -
more » « less
Abstract From California to British Columbia, the Pacific Northwest coast bears an omnipresent earthquake and tsunami hazard from the Cascadia subduction zone. Multiple lines of evidence suggests that magnitude eight and greater megathrust earthquakes have occurred ‐ the most recent being 321 years ago (i.e., 1700 A.D.). Outstanding questions for the next great megathrust event include where it will initiate, what conditions are favorable for rupture to span the convergent margin, and how much slip may be expected. We develop the first 3‐D fully dynamic rupture simulations for the Cascadia subduction zone that are driven by fault stress, strength and friction to address these questions. The initial dynamic stress drop distribution in our simulations is constrained by geodetic coupling models, with segment locations taken from geologic analyses. We document the sensitivity of nucleation location and stress drop to the final seismic moment and coseismic subsidence amplitudes. We find that the final earthquake size strongly depends on the amount of slip deficit in the central Cascadia region, which is inferred to be creeping interseismically, for a given initiation location in southern or northern Cascadia. Several simulations are also presented here that can closely approximate recorded coastal subsidence from the 1700 A.D. event without invoking localized high‐stress asperities along the down‐dip locked region of the megathrust. These results can be used to inform earthquake and tsunami hazards for not only Cascadia, but other subduction zones that have limited seismic observations but a wealth of geodetic inference.
-
more » « less
Abstract Slow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014
M w7.3 Guerrero, Mexico earthquake which was preceded by aM w7.6 SSE. However, the underlying physics connecting SSEs to earthquakes remains elusive. Here, we link 3D slow‐slip cycle models with dynamic rupture simulations across the geometrically complex flat‐slab Cocos plate boundary. Our physics‐based models reproduce key regional geodetic and teleseismic fault slip observations on timescales from decades to seconds. We find that accelerating SSE fronts transiently increase shear stress at the down‐dip end of the seismogenic zone, modulated by the complex geometry beneath the Guerrero segment. The shear stresses cast by the migrating fronts of the 2014M w7.6 SSE are significantly larger than those during the three previous episodic SSEs that occurred along the same portion of the megathrust. We show that the SSE transient stresses are large enough to nucleate earthquake dynamic rupture and affect rupture dynamics. However, additional frictional asperities in the seismogenic part of the megathrust are required to explain the observed complexities in the coseismic energy release and static surface displacements of the Guerrero earthquake. We conclude that it is crucial to jointly analyze the long‐ and short‐term interactions and complexities of SSEs and megathrust earthquakes across several (a)seismic cycles accounting for megathrust geometry. Our study has important implications for identifying earthquake precursors and understanding the link between transient and sudden megathrust faulting processes.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s00024-018-1991-x
