Coastal subsidence, dating of plant remains and tree rings, and evidence for tsunami inundation point to coseismic activity on a sizable portion of the Cascadia subduction zone around three centuries ago. A tsunami of remote origin in 1700 C.E., probably from Cascadia, caused flooding and damage in Japan. In previous modeling, this transpacific evidence was found most simply explained by one Cascadia rupture about 1,000 km long. Here I model tens of thousands of ruptures and simulate their subsidence and tsunami signals and show that it is possible that the earthquake was part of a sequence of several events. Partial rupture of ∼400 km offshore southern Oregon and northern California in one large M ≥ 8.7 earthquake can explain the tsunami in Japan without conflicting with the subsidence. As many as four more earthquakes with M ≤ 8.7 can complete the subsidence signal without their tsunamis being large enough to be recorded in Japan. The purpose of this study is not to find a single, most likely, scenario or disprove the single‐rupture hypothesis favored by alternative evidence such as turbidites. Rather, it demonstrates that a multiple rupture sequence may explain part of the available data, and therefore cannot be discounted. Given the gaps in the presently available estimates of subsidence it is also possible that segments of the megathrust, for example from Copalis to the Strait of Juan de Fuca, did not rupture in 1700. The findings have significant implications for Cascadia geodynamics and how earthquake and tsunami hazards in the region are quantified.
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- Award ID(s):
- 1663769
- NSF-PAR ID:
- 10446007
- 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
-
more » « less
Abstract -
more » « less
Abstract Given recent advances in geodetic data, interseismic locking models along the megathrust now become useful to qualitatively evaluate future earthquake potential. However, an individual earthquake's true rupture potential is challenging, as it depends on more than just a static image of prior locking. Here, we test the determinism of interseismic locking models using spontaneous rupture simulations and the well‐resolved processes associated with the 2012 moment magnitude (Mw) 7.6 Nicoya earthquake. To do so, we estimate initial megathrust stress from locking by assuming that the entire slip deficit will be released in the next megathrust earthquake. Then we initiate spontaneous ruptures at the hypocenter of the 2012 Nicoya earthquake. We find scenarios that approximate the same coseismic slip distribution and final earthquake moment magnitude as obtained from seismic and geodetic observations, demonstrating that deriving potential rupture scenarios from interseismic locking is feasible. We also find that spontaneous rupture scenarios from different locking models differ in moment rate duration and thus ground motion prediction, although the final slip distribution and moment magnitude were similar. The results highlight that quantifying rupture scenarios and ground motions from reliable locking models by dynamic rupture simulations can be an effective tool for seismic hazard assessment in subduction zones.
-
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. -
Abstract On 19 September 2022, a major earthquake struck the northwestern Michoacán segment along the Mexican subduction zone. A slip model is obtained that satisfactorily explains geodetic, teleseismic, and tsunami observations of the 2022 event. The preferred model has a compact large-slip patch that extends up-dip and northwestward from the hypocenter and directly overlaps a 1973 Mw 7.6 rupture. Slip is concentrated offshore and below the coast at depths from 10 to 30 km with a peak value of ∼2.9 m, and there is no detected coseismic slip near the trench. The total seismic moment is 3.1×1020 N·m (Mw 7.6), 72% of which is concentrated in the first 30 s. Most aftershocks are distributed in an up-dip area of the mainshock that has small coseismic slip, suggesting near-complete strain release in the large-slip patch. Teleseismic P waveforms of the 2022 and 1973 earthquakes are similar in duration and complexity with high cross-correlation coefficients of 0.68–0.98 for long P to PP signal time windows, indicating that the 2022 earthquake is a quasi-repeat of the 1973 earthquake, possibly indicating persistent frictional properties. Both the events produced more complex P waveforms than comparable size events along Guerrero and Oaxaca, reflecting differences in patchy locking of the Mexican megathrust.more » « less
-
more » « less
Abstract There is a strong need to model potential rupture behaviors for the next Cascadia megathrust earthquake. However, there exists significant uncertainty regarding the extent of downdip rupture and rupture speed. To address this problem, we study how the transition region (i.e., the gap), which separates the locked from slow‐slip regions, influences coseismic rupture propagation using 2‐D dynamic rupture simulations governed by a slip‐weakening friction law. We show that rupture propagation through the gap is strongly controlled by the amount of accumulated tectonic initial shear stress and gap friction level. A large amplitude negative dynamic stress drop is needed to arrest downdip rupture. We also observe downdip supershear rupture when the gradient in effective normal stress from the locked to slow‐slip regions is dramatic. Our results justify kinematic rupture models that extend below the gap and suggests the possibility of high‐frequency energy radiation during the next Cascadia megathrust earthquake.
