Search

Strong tsunami excitation from slow rupture of shallow subduction zone faults is recognized as a key concern for tsunami hazard assessment. Three months after the 22 July 2020 magnitude 7.8 thrust earthquake struck the plate boundary below the Shumagin Islands, Alaska, a magnitude 7.6 aftershock ruptured with complex intraplate faulting. Despite the smaller size and predominantly strike-slip faulting mechanism inferred from seismic waves for the aftershock, it generated much larger tsunami waves than the mainshock. Here we show through detailed analysis of seismic, geodetic, and tsunami observations of the aftershock that the event implicated unprecedented source complexity, involving weakly tsunamigenic fast rupture of two intraplate faults located below and most likely above the plate boundary, along with induced strongly tsunamigenic slow thrust slip on a third fault near the shelf break likely striking nearly perpendicular to the trench. The thrust slip took over 5 min, giving no clear expression in seismic or geodetic observations while producing the sizeable far-field tsunami. The 19 Oct 2020 MW 7.6 Shumagin earthquake involved unprecedented source complexity with two fast ruptures straddling the megathrust and strong tsunami excitation from a long-duration upper plate thrust rupture undetected by seismic and geodetic data. [ABSTRACT FROM AUTHOR]

The Shumagin seismic gap along the Alaska Peninsula experienced a major, MW 7.8, interplate thrust earthquake on 22 July 2020. Several available finite‐fault inversions indicate patchy slip of up to 4 m at 8–48 km depth. There are differences among the models in peak slip and absolute placement of slip on the plate boundary, resulting from differences in data distributions, model parameterizations, and inversion algorithms. Two representative slip models obtained from inversions of large seismic and geodetic data sets produce very different tsunami predictions at tide gauges and deep‐water pressure sensors (DART stations), despite having only secondary differences in slip distribution. This is found to be the result of the acute sensitivity of the tsunami excitation for rupture below the continental shelf in proximity to an abrupt shelf break. Iteratively perturbing seismic and geodetic inversions by constraining fault model extent along dip and strike, we obtain an optimal rupture model compatible with teleseismic P and SH waves, regional three‐component broadband and strong‐motion seismic recordings, hr‐GNSS time series and static offsets, as well as tsunami recordings at DART stations and regional and remote tide gauges. Slip is tightly bounded between 25 and 40 km depth, the up‐dip limit of slip in the earthquake is resolved to be well‐inland of the shelf break, and the rupture extent along strike is well‐constrained. The coseismic slip increased Coulomb stress on the shallow plate boundary extending to the trench, but the frictional behavior of the megathrust below the continental slope remains uncertain. Plain Language Summary: Several studies on the 22 July 2020 MW 7.8 thrust event have included finite‐fault inversions, yielding generally consistent slip models. However, none of these studies have fully considered the seismic and geodetic data together with the tsunami measurements at tide gauges and deep‐water pressure sensors. We select two representative slip models obtained from prior seismic and geodetic inversions for tsunami modeling and find that they give distinct tsunami predictions despite having only a minor shift of slip by ∼20 km relative to shelf break. Large differences in tsunami excitation occur when seafloor deformation extends to the continental slope under deeper water. We perform iterative inversion of seismic and geodetic observations and forward modeling of tsunami signals to obtain an optimal slip model compatible with teleseismic waves, regional seismic recordings, hr‐GNSS time series and static offsets, along with tsunami recordings at DARTs and tide gauges. Precise placement of the slip distribution beneath the continental shelf controls the strength and timing of seaward tsunamis, with waveforms varying nonlinearly with seafloor deformation extending across the shelf break. The earthquake increased Coulomb stress on the shallow megathrust, but uncertainty in the shallow frictional behavior leaves it unclear whether a future large event can occur there. Key Points: Published finite‐fault models for the 2020 MW 7.8 earthquake predict varying tsunami signals despite having similar slip patternsThe tsunami excitation is very sensitive to the absolute placement of slip relative to the shelf break and along strikeIteration of slip inversion parameters and tsunami predictions optimizes the model to fit all data well with slip only beneath the shelf [ABSTRACT FROM AUTHOR]