The 15 April 2015
Most large tsunamis are generated by earthquakes on offshore plate boundary megathrusts. The primary factors influencing tsunami excitation are the seismic moment, faulting geometry, and depth of the faulting. Efforts to provide rapid tsunami warning have emphasized seismic and geodetic methods for quickly determining the event size and faulting geometry. It remains difficult to evaluate the updip extent of rupture, which has significant impact on tsunami excitation. Teleseismic
- Award ID(s):
- 1802364
- NSF-PAR ID:
- 10455667
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 124
- Issue:
- 10
- ISSN:
- 2169-9313
- Page Range / eLocation ID:
- p. 10583-10609
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract M w7.8 Nepal Gorkha earthquake occurred on a shallowly dipping portion of the Main Himalayan Thrust (MHT). Notable features of the event include (1) the dominance of a slip pulse of about 6‐s duration that unlocked the lower edge of the MHT and (2) the near‐horizontal fault geometry, which, combined with proximity of the free surface, allows surface‐reflected phases to break the across‐fault symmetries of the seismic wavefield. Our dynamic rupture simulations in an elastoplastic medium yield earthquake parameters comparable to those deduced from kinematic inversions, including seismic moment and rupture velocity. The simulations reproduce pulse‐like behavior predicting pulse widths in agreement with those kinematic studies and supporting an interpretation in which the pulse‐like time dependence of slip is principally controlled by rupture geometry. This inference is strongly supported by comparison of synthetic ground velocity with the near‐field high‐rate GPS recording at station KKN4, which shows close agreement in pulse width, amplitude, and pulse shape. That comparison also constrains the updip extent of rupture and disfavors significant coseismic slip on the shallow ramp segment. Over most of the rupture length, the simulated rupture propagates at a near‐constant maximum velocity (~90% of theS wave speed) that is controlled by the antiplane geometry and off‐fault plastic yielding. Simulations also reveal the role of reflected seismic waves from the free surface, which may have contributed ~30% elongation of the slip pulse, and show the potential for significant free‐surface interaction effects in shallow events of similar geometry. -
Abstract When slip at shallow depth occurs during large subduction zone thrust events,
P wave energy enters the water layer and establishespwP , the reverberating waves called “water bounces” that followpP . For water depths ≥5–6 km (i.e., near the trench) above the shallow slip,pwP manifests in a strong ~10‐s period ringing that can persist for minutes into the teleseismicP wave coda at all azimuths. Deeper slip can generate shorter‐periodpwP ringing at trenchward azimuths. At large distances,P codawindows have several‐minute‐long intervals free of secondary arrivals. We consider rmsP coda/rmsP amplitude ratios at distances from 80° to 120° as a potential proxy for occurrence of shallow slip for 39M W 7.5+ megathrust earthquakes from 1990 to 2016 with estimated slip distributions. Ratios for the 15‐ to 7‐s‐period band have a strong bimodal distribution, with higher averageP coda/P amplitudes observed for ruptures with slip extending to shallow depth. -
Abstract The Shumagin seismic gap along the Alaska Peninsula experienced a major,
M W 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 teleseismicP andSH 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. -
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. -
Abstract Recent development of dense strong‐motion networks and seismic array processing enables rapid tsunami predictions based on the back‐projection (BP) approach. We develop a multi‐array local BP method (MLBP) using seismic networks with epicentral distance from 0.7° to 3.5°. The local BPs using individual arrays are first calculated and are then merged into a single image of the rupture process. This multi‐array approach circumvents the issue of artifact in single‐array BP caused by the overlapping of multiple phases and coda waves. Based on the local BP approach, tsunami predictions are available 7 min after the origin time. Case studies of the 2003 Mw 8.1 Tokachi‐oki and the 2011 Mw 9.0 Tohoku earthquakes show that their rupture zones are well resolved and are comparable with principal slip areas inferred from tsunami observations. The amplitude and arrival time errors of the predicted tsunami waves are within −1.59 to 3.74 m and −10.0 to 10.0 min, which are sufficiently small for warning purposes.