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.
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.
more » « less- NSF-PAR ID:
- 10371999
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Abstract Lithology and microfossil biostratigraphy beneath the marshes of a central Oregon estuary limit geophysical models of Cascadia megathrust rupture during successive earthquakes by ruling out >0.5 m of coseismic coastal subsidence for the past 2000 yr. Although the stratigraphy in cores and outcrops includes as many as 12 peat-mud contacts, like those commonly inferred to record subsidence during megathrust earthquakes, mapping, qualitative diatom analysis, foraminiferal transfer function analysis, and 14C dating of the contacts failed to confirm that any contacts formed through subsidence during great earthquakes. Based on the youngest peat-mud contact’s distinctness, >400 m distribution, ∼0.6 m depth, and overlying probable tsunami deposit, we attribute it to the great 1700 CE Cascadia earthquake and(or) its accompanying tsunami. Minimal changes in diatom assemblages from below the contact to above its probable tsunami deposit suggest that the lower of several foraminiferal transfer function reconstructions of coseismic subsidence across the contact (0.1–0.5 m) is most accurate. The more limited stratigraphic extent and minimal changes in lithology, foraminifera, and(or) diatom assemblages across the other 11 peat-mud contacts are insufficient to distinguish them from contacts formed through small, gradual, or localized changes in tide levels during river floods, storm surges, and gradual sea-level rise. Although no data preclude any contacts from being synchronous with a megathrust earthquake, the evidence is equally consistent with all contacts recording relative sea-level changes below the ∼0.5 m detection threshold for distinguishing coseismic from nonseismic changes.more » « less
-
Abstract Stratigraphic, lithologic, foraminiferal, and radiocarbon analyses indicate that at least four abrupt mud-over-peat contacts are recorded across three sites (Jacoby Creek, McDaniel Creek, and Mad River Slough) in northern Humboldt Bay, California, USA (∼44.8°N, −124.2°W). The stratigraphy records subsidence during past megathrust earthquakes at the southern Cascadia subduction zone ∼40 km north of the Mendocino Triple Junction. Maximum and minimum radiocarbon ages on plant macrofossils from above and below laterally extensive (>6 km) contacts suggest regional synchroneity of subsidence. The shallowest contact has radiocarbon ages that are consistent with the most recent great earthquake at Cascadia, which occurred at 250 cal yr B.P. (1700 CE). Using Bchron and OxCal software, we model ages for the three older contacts of ca. 875 cal yr B.P., ca. 1120 cal yr B.P., and ca. 1620 cal yr B.P. For each of the four earthquakes, we analyze foraminifera across representative mud-over-peat contacts selected from McDaniel Creek. Changes in fossil foraminiferal assemblages across all four contacts reveal sudden relative sea-level (RSL) rise (land subsidence) with submergence lasting from decades to centuries. To estimate subsidence during each earthquake, we reconstructed RSL rise across the contacts using the fossil foraminiferal assemblages in a Bayesian transfer function. The coseismic subsidence estimates are 0.85 ± 0.46 m for the 1700 CE earthquake, 0.42 ± 0.37 m for the ca. 875 cal yr B.P. earthquake, 0.79 ± 0.47 m for the ca. 1120 cal yr B.P. earthquake, and ≥0.93 m for the ca. 1620 cal yr B.P. earthquake. The subsidence estimate for the ca. 1620 cal yr B.P. earthquake is a minimum because the pre-subsidence paleoenvironment likely was above the upper limit of foraminiferal habitation. The subsidence estimate for the ca. 875 cal yr B.P. earthquake is less than (<50%) the subsidence estimates for other contacts and suggests that subsidence magnitude varied over the past four earthquake cycles in southern Cascadia.more » « less
-
more » « less
Abstract Tsunamis from earthquakes of various magnitudes have affected Cascadia in the past. Simulations of
M w 7.5–9.2 earthquake constrained by earthquake rupture physics and geodetic locking models show thatM w ≥ 8.5 events initiating in the middle segments of the subduction zone can create coastal tsunami amplitudes comparable to those from the largest expected event. Our rupture and tsunami simulations reveal that the concave coastline geometry of the Pacific Northwest coastline focuses tsunami energy between latitudes 44° and 45° in Oregon. The possible coastal tsunami amplitudes are largely insensitive to the choice of slip model for a given magnitude. These results are useful for identifying the most hazardous segments of the subduction zone and demonstrate that a worst‐case rupture scenario does not uniquely yield the worst‐case tsunami scenario at a given location. -
more » « less
Abstract A great earthquake struck the Semidi segment of the plate boundary along the Alaska Peninsula on 29 July 2021, re‐rupturing part of the 1938 rupture zone. The 2021
M W 8.2 Chignik earthquake occurred just northeast of the 22 July 2020M W 7.8 Simeonof earthquake, with little slip overlap. Analysis of teleseismicP andSH waves, regional Global Navigation Satellite System (GNSS) displacements, and near‐field and far‐field tsunami observations provides a good resolution of the 2021 rupture process. During ∼60‐s long faulting, the slip was nonuniformly distributed along the megathrust over depths from 32 to 40 km, with up to ∼12.9‐m slip in an ∼170‐km‐long patch. The 40–45 km down‐dip limit of slip is well constrained by GNSS observations along the Alaska Peninsula. Tsunami observations preclude significant slip from extending to depths <25 km, confining all coseismic slip to beneath the shallow continental shelf. Most aftershocks locate seaward of the large‐slip zones, with a concentration of activity up‐dip of the deeper southwestern slip zone. Some localized aftershock patches locate beneath the continental slope. The surface‐wave magnitudeM S of 8.1 for the 2021 earthquake is smaller thanM S = 8.3–8.4 for the 1938 event. Seismic and tsunami data indicate that slip in 1938 was concentrated in the eastern region of its aftershock zone, extending beyond the Semidi Islands, where the 2021 event did not rupture.
