skip to main content


Title: Supershear shock front contribution to the tsunami from the 2018 M w 7.5 Palu, Indonesia earthquake
SUMMARY Hazardous tsunamis are known to be generated predominantly at subduction zones. However, the 2018 Mw 7.5 Palu (Indonesia) earthquake on a strike-slip fault generated a tsunami that devastated the city of Palu. The mechanism by which this tsunami originated from such an earthquake is being debated. Here we present near-field ground motion (GPS) data confirming that the earthquake attained supershear speed, i.e. a rupture speed greater than the shear wave speed of the host medium. We subsequently study the effect of this supershear rupture on tsunami generation by coupling the ground motion to a 1-D non-linear shallow-water wave model accounting for both time-dependent bathymetric displacement and velocity. With the local bathymetric profile of Palu bay around a tidal station, our simulations reproduce the tsunami arrival and motions observed by CCTV cameras. We conclude that Mach (shock) fronts, generated by the supershear speed, interacted with the bathymetry and contributed to the tsunami.  more » « less
Award ID(s):
1753249
NSF-PAR ID:
10332935
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
230
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
2089 to 2097
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY In the Gulf of California, Mexico, the relative motion across the North America–Pacific boundary is accommodated by a series of marine transform faults and spreading centres. About 40 M> 6 earthquakes have occurred in the region since 1960. On 2009 August 3, an Mw 6.9 earthquake occurred near Canal de Ballenas in the region. The earthquake was a strike-slip event with a shallow hypocentre that is likely close to the seafloor. In contrast to an adjacent M7 earthquake, this earthquake triggered a ground-motion-based earthquake early warning algorithm being tested in southern California (∼600 km away). This observation suggests that the abnormally large ground motions and dynamic strains observed for this earthquake relate to its rupture properties. To investigate this possibility, we image the rupture process and resolve the slip distribution of the event using a P-wave backprojection approach and a teleseismic, finite-fault inversion method. Results from these two independent analyses indicate a relatively simple, unilateral rupture propagation directed along-strike in the northward direction. However, the average rupture speed is estimated around 4 km s−1, suggesting a possible supershear rupture. The supershear speed is also supported by a Rayleigh wave Mach cone analysis, although uncertainties in local velocity structure preclude a definitive conclusion. The Canal de Ballenas earthquake dynamically triggered seismicity at multiple sites in California, with triggering response characteristics varying from location-to-location. For instance, some of the triggered earthquakes in California occurred up to 24 hr later, suggesting that nonlinear triggering mechanisms likely have modulated their occurrence. 
    more » « less
  2. Tsunami generation from earthquake-induced seafloor deformations has long been recognized as a major hazard to coastal areas. Strike-slip faulting has generally been considered insufficient for triggering large tsunamis, except through the generation of submarine landslides. Herein, we demonstrate that ground motions due to strike-slip earthquakes can contribute to the generation of large tsunamis (>1 m), under rather generic conditions. To this end, we developed a computational framework that integrates models for earthquake rupture dynamics with models of tsunami generation and propagation. The three-dimensional time-dependent vertical and horizontal ground motions from spontaneous dynamic rupture models are used to drive boundary motions in the tsunami model. Our results suggest that supershear ruptures propagating along strike-slip faults, traversing narrow and shallow bays, are prime candidates for tsunami generation. We show that dynamic focusing and the large horizontal displacements, characteristic of strike-slip earthquakes on long faults, are critical drivers for the tsunami hazard. These findings point to intrinsic mechanisms for sizable tsunami generation by strike-slip faulting, which do not require complex seismic sources, landslides, or complicated bathymetry. Furthermore, our model identifies three distinct phases in the tsunamic motion, an instantaneous dynamic phase, a lagging coseismic phase, and a postseismic phase, each of which may affect coastal areas differently. We conclude that near-source tsunami hazards and risk from strike-slip faulting need to be re-evaluated.

     
    more » « less
  3. Abstract

    Supershear earthquakes, propagating faster than the Earth's shear wave velocity, can generate strong ground motion at distances far from the ruptured fault. Despite the hazards associated with these earthquakes, the exact fault properties that govern their occurrence are not well constrained. Although early studies associated supershear ruptures with simple fault geometries, recent dynamic rupture models have revealed a supershear transition mechanism over complex fault geometries such as fault stepovers. Here we present the first observation of a supershear transition on a fault stepover system during the 2017 Mw7.7 Komandorsky Islands earthquake. Using a high‐resolution back‐projection technique, we find that the earthquake's rupture velocity accelerates from 2.1 to 5.0 km/s between two offset faults, demonstrating the viability of a new supershear transition mechanism occurring in nature. Given the fault complexity of the Earth's transform plate boundaries, this result may improve our understanding of supershear rupture processes and their associated hazards.

     
    more » « less
  4. Abstract

    Megathrust earthquakes exhibit a ubiquitous seismic radiation style: low‐frequency (LF) seismic energy is efficiently emitted from the shallowest portion of the fault, whereas high‐frequency (HF) seismic energy is efficiently emitted from the deepest part of the fault. Although this is observed in many case‐specific studies, we show that it is ubiquitous in global megathrust earthquakes between 1995 and 2021. Previous studies have interpreted this as an effect of systematic depth variation in either the plate interface frictional properties (Lay et al., 2012) or the P wavespeeds (Sallarès & Ranero, 2019). This work suggests an alternative hypothesis: the interaction between waves and ruptures due to the Earth's free surface is the leading mechanism that generates this behavior. Two‐dimensional dynamic rupture simulations of subduction zone earthquakes support this hypothesis. Our simulations show that the interaction between the seismic waves reflected at the Earth's free surface and the updip propagating rupture results in LF radiation at the source. In contrast, the downdip propagation of rupture is less affected by the free surface and is thus dominated by HF radiation typical of buried faults. To a second degree, the presence of a realistic Earth structure derived from P‐wave velocity (VP) tomographic images and realistic VP/VSratio estimated in boreholes further enhances the contrast in source radiation. We conclude that the Earth's free surface is necessary to explain the observed megathrust earthquake radiation style, and the realistic structure of subduction zone is necessary to better predict earthquake ground motion and tsunami potential.

     
    more » « less
  5. null (Ed.)
    Abstract The Mw 7.1 47 km deep earthquake that occurred on 30 November 2018 had deep societal impacts across southcentral Alaska and exhibited phenomena of broad scientific interest. We document observations that point to future directions of research and hazard mitigation. The rupture mechanism, aftershocks, and deformation of the mainshock are consistent with extension inside the Pacific plate near the down‐dip limit of flat‐slab subduction. Peak ground motions >25%g were observed across more than 8000  km2, though the most violent near‐fault shaking was avoided because the hypocenter was nearly 50 km below the surface. The ground motions show substantial variation, highlighting the influence of regional geology and near‐surface soil conditions. Aftershock activity was vigorous with roughly 300 felt events in the first six months, including two dozen aftershocks exceeding M 4.5. Broad subsidence of up to 5 cm across the region is consistent with the rupture mechanism. The passage of seismic waves and possibly the coseismic subsidence mobilized ground waters, resulting in temporary increases in stream flow. Although there were many failures of natural slopes and soils, the shaking was insufficient to reactivate many of the failures observed during the 1964 M 9.2 earthquake. This is explained by the much shorter duration of shaking as well as the lower amplitude long‐period motions in 2018. The majority of observed soil failures were in anthropogenically placed fill soils. Structural damage is attributed to both the failure of these emplaced soils as well as to the ground motion, which shows some spatial correlation to damage. However, the paucity of instrumental ground‐motion recordings outside of downtown Anchorage makes these comparisons challenging. The earthquake demonstrated the challenge of issuing tsunami warnings in complex coastal geographies and highlights the need for a targeted tsunami hazard evaluation of the region. The event also demonstrates the challenge of estimating the probabilistic hazard posed by intraslab earthquakes. 
    more » « less