Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably‐dipping planar and variably‐curved nonplanar megathrusts using the volumetric, high‐order accurate code
Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.
more » « less- PAR ID:
- 10366827
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 3
- Issue:
- 2
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly‐dipping megathrusts to host larger earthquakes than steeply‐dipping ones. Under elevated pore pressure and less strongly velocity‐weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly‐dipping and sharply‐curved megathrusts host multi‐period supercycles of slow‐to‐fast, small‐to‐large slip events under higher effective stresses and more strongly velocity‐weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second‐order effects from structurally‐derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co‐ and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults. -
more » « less
Abstract Long‐term slow slip events have been observed at several subduction zones around the globe, where they play an integral part in strain release along megathrust faults. Nevertheless, evidence for long‐term slow slip has remained elusive in the Cascadia subduction zone. Here we conduct a systematic analysis of 13 years of GNSS time series data from 2006 to 2019 and present evidence of at least one low‐amplitude long‐term slow slip event on the Cascadia subduction zone, with the possibility of others that are less resolved. Starting in mid‐2012, a 1.5‐year transient is observed in southern Cascadia, with a group of coastal GNSS stations moving ∼2 mm to the west. The data are modeled as a Mw 6.4 slow slip event occurring at 15–35 km depth on the plate interface, just updip of previously recognized short‐term slow slip and tremor. The event shares many characteristics with similar long‐term transient events on the Nankai subduction zone. However, the total fault slip amplitude is an order‐of‐magnitude smaller in Cascadia when compared to large events elsewhere, making long‐term slow slip detection challenging in Cascadia. While there are other westward long‐duration transients in the refined data set, the surface displacements are below the level of the noise or are limited spatially to a few neighboring stations, making interpretation unclear.
-
more » « less
Abstract Within the fore‐arc of the Cascadia Subduction Zone, there are significant along‐strike differences in the orientation of splay faults, sediment consolidation, and fault roughness. Here, we use dynamic rupture simulations of megathrust earthquakes on different realizations of a fault system that incorporate fore‐arc properties representative of offshore Oregon and Washington to estimate how splay faults may behave in future megathrust earthquakes in Cascadia. While splay faults were activated in all of our simulations, splay orientation is a primary control on slip amplitude. Seaward vergent faults accommodate significant amounts of slip resulting in large seafloor uplift and significantly larger tsunami amplitudes. For example, our median tsunami heights including splay faults are about a factor of two larger than those that did not include splay fault deformation. We suggest that there is an urgent need to revisit existing approaches to tsunami hazard assessment in Cascadia to include the influence of splay faults.
-
more » « less
Abstract 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.
-
more » « less
Because splay faults branch at a steep dip angle from the plate-boundary décollement in an accretionary wedge, their coseismic displacement can potentially result in larger tsunamis with distinct characteristics compared to megathrust-only fault ruptures, posing an enhanced hazard to coastal communities. Elsewhere, there is evidence of coseismic slip on splay faults during many of the largest subduction zone earthquakes, but our understanding of potentially active splay faults and their hazards at the Cascadia subduction zone remains limited. To identify the most recently active splay faults at Cascadia, we conduct stratigraphic and structural interpretations of near-surface deformation in the outer accretionary wedge for the ~400 km along-strike length of the landward vergence zone. We analyze recently acquired high-frequency sparker seismic data and crustal-scale multi-channel seismic data to examine the record of deformation in shallow slope basins and the upper ~1 km of the surrounding accreted sediments and to investigate linkages to deeper décollement structure. We present a new fault map for widest, most completely locked portion of Cascadia from 45 to 48°N latitude, which documents the distribution of faults that show clear evidence of recent late Quaternary activity. We find widespread evidence for active splay faulting up to 30 km landward of the deformation front, in what we define as the active domain, and diminished fault activity landward outside of this zone. The abundance of surface-deforming splay faults in the active outer wedge domain suggests Cascadia megathrust events may commonly host distributed shallow rupture on multiple splay faults located within 30 km of the deformation front.
