Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


Title: Performance of Slab Geometry Constraints on Rapid Geodetic Slip Models, Tsunami Amplitude, and Inundation Estimates in Cascadia
Tsunamigenic megathrust earthquakes along the Cascadia subduction zone present a major hazard concern. We can better prepare to model the earthquake source in a rapid manner by imbuing fault geometry constraints based on prior knowledge and by evaluating the capabilities of using existing GNSS sensors. Near-field GNSS waveforms have shown promise in providing rapid coarse finite-fault model approximations of the earthquake rupture that can improve tsunami modeling and response time. In this study, we explore the performance of GNSS derived finite-fault inversions and tsunami forecasting predictions in Cascadia that highlights the impact and potential of geodetic techniques and data in operational earthquake and tsunami monitoring. We utilized 1300 Cascadia earthquake simulations (FakeQuakes) that provide realistic (M7.5-9.3) rupture scenarios to assess how feasibly finite-fault models can be obtained in a rapid earthquake early warning and tsunami response context. A series of fault models with rectangular dislocation patches spanning the Cascadia megathrust area is added to the GFAST inversion algorithm to calculate slip for each earthquake scenario. Another method used to constrain the finite-fault geometry is from the GNSS-derived CMT fault plane solution. For the Cascadia region, we show that fault discretization using two rectangular segments approximating the megathrust portion of the subduction zone leads to improvements in modeling magnitude, fault slip, tsunami amplitude, and inundation. In relation to tsunami forecasting capabilities, we compare coastal amplitude predictions spanning from Vancouver Island (Canada) to Northern California (USA). Generally, the coastal amplitudes derived using fault parameters from the CMT solutions show an overestimation bias compared to amplitudes derived from the fixed slab model. We also see improved prediction values of the run-up height and maximum amplitude at 10 tide gauge stations using the fixed slab model as well.  more » « less
Award ID(s):
2103713
PAR ID:
10579415
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Seismica
Date Published:
Journal Name:
Seismica
Volume:
2
Issue:
4
ISSN:
2816-9387
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The USGS Powell Center Cascadia earthquake hazards working group compiled published onshore and offshore paleoseismic data along the Cascadia subduction zone, spanning sites from Vancouver Island to the Mendocino triple junction. Evidence for megathrust rupture includes coastal land-level change, tsunami inundation, onshore shaking proxies such as landslides or liquefaction, and offshore shaking proxies such as marine turbidites. The quality of paleoseismic data for megathrust rupture along the Cascadia subduction zone collected over the past three decades varies because analytical capabilities and data collection methodologies have evolved. Thus, as part of the compilation, we also present a ranking scheme to assess the quality of age estimates and evidence for great megathrust rupture. With the age ranking scheme, we ask: "How well is a proposed paleoseismic event dated?" based on the materials and methods used. With the evidence ranking scheme, we ask: "How confident are we that a proposed event is, in fact, the result of a Cascadia megathrust rupture?" based on the sedimentological characteristics, correlation, and mapping. The evidence ranking scheme also helps to evaluate possible alternative mechanisms for creating paleoseismic evidence such as crustal fault, intraslab, or distant tsunamigenic earthquake." 
    more » « less
  2. Complex interactions between strain accumulation, fault slip, and fluid migration influence shallow subduction zone dynamics. Using a decade of continuous ambient seismic data from Cascadia’s seafloor observatories, we identified distinct regional variations in subduction dynamics. Northern Cascadia exhibits a fully locked megathrust with persistent strain accumulation, while central Cascadia displays a slow slip event on protothrusts and rapid fluid migration along fault systems in the overriding plate. Effective fluid transport through the décollement and the Alvin Canyon Fault likely modulates the earthquake behavior but does not cause slow slip events on the megathrust and likely stabilizes large earthquakes, promoting rupture arrest. 
    more » « less
  3. Abstract Great subduction earthquakes (Mw ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward‐ or landward‐vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward‐vergent splays host more slip than those with steeper, landward‐vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip‐ and vergence‐dependence of splay rupture: static (mis)alignment relative to far‐field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip‐propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward‐vergent splays and inhibit backward branching onto landward‐vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments. 
    more » « less
  4. The Mediterranean Hellenic Arc subduction zone (HASZ) has generatedseveral Mw>=8 earthquakes and tsunamis.Seismic-probabilistic tsunami hazard assessment typically utilizesuniform or stochastic earthquake models, which may not represent dynamicrupture and tsunami generation complexity. We present an ensemble of ten3D dynamic rupture earthquake scenarios for the HASZ, utilizing arealistic slab geometry. Our simplest models use uniform along-arcpre-stresses or a single circular initial stress asperity. We thenintroduce progressively more complex models varying initial shear stressalong-arc, multiple asperities based on scale-dependent critical slipweakening distance, and a most complex model blending all aforementionedheterogeneities. Thereby, regional initial conditions are constrainedwithout relying on detailed geodetic locking models. Varying hypocenterlocations in the simplest, homogeneous model leads to different rupturespeeds and moment magnitudes. We observe dynamic fault slip penetratingthe shallow slip-strengthening region and affecting seafloor uplift.Off-fault plastic deformation can double vertical seafloor uplift. Asingle-asperity model generates a Mw~8 scenarioresembling the 1303 Crete earthquake. Using along-strike varying initialstresses results in Mw~8.0-8.5 dynamic rupture scenarioswith diverse slip rates and uplift patterns. The model with the mostheterogeneous initial conditions yields a Mw~7.5scenario. Dynamic rupture complexity in prestress and fracture energytends to lower earthquake magnitude but enhances tsunamigenicdisplacements. Our results offer insights into the dynamics of potentiallarge Hellenic Arc megathrust earthquakes and may inform futurephysics-based joint seismic and tsunami hazard assessments. 
    more » « less
  5. Abstract Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half‐space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low‐rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints. 
    more » « less