The 2018 Palu tsunami contributed significantly to the devastation caused by the associated
- Award ID(s):
- 2026099
- Publication Date:
- NSF-PAR ID:
- 10269145
- Journal Name:
- Bulletin of the Seismological Society of America
- Volume:
- 111
- Issue:
- 2
- Page Range or eLocation-ID:
- 1032 to 1042
- ISSN:
- 0037-1106
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract 7.5 earthquake. This began a debate about how the moderate size earthquake triggered such a large tsunami within Palu Bay, with runups of more than 10 m. The possibility of a large component of vertical coseismic deformation and submarine landslides have been considered as potential explanations. However, scarce instrumental data have made it difficult to resolve the potential contributions from either type of source. We use tsunami waveforms derived from social media videos in Palu Bay to model the possible sources of the tsunami. We invert InSAR data with different fault geometries and use the resulting seafloor displacements to simulate tsunamis. The coseismic sources alone cannot match both the video‐derived time histories and surveyed runups. Then we conduct a tsunami source inversion using the video‐derived time histories and a tide gauge record as inputs. We specify hypothetical landslide locations and solve for initial tsunami elevation. Our results, validated with surveyed runups, show that a limited number of landslides in southern Palu Bay are sufficient to explain the tsunami data. The Palu tsunami highlights the difficulty in accurately capturing with tide gauges the amplitude and timing ofmore » -
Abstract Large earthquakes can construct mountainous topography by inducing rock uplift but also erode mountains by causing landslides. Observations following the 2008 Wenchuan earthquake show that landslide volumes in some cases match seismically induced uplift, raising questions about how the actions of individual earthquakes accumulate to build topography. Here we model the two‐dimensional surface displacement field generated over a full earthquake cycle accounting for coseismic deformation, postseismic relaxation, landslide erosion, and erosion‐induced isostatic compensation. We explore the related volume balance across different seismotectonic and topographic conditions and revisit the Wenchuan case in this context. The ratio (Ω) between landslide erosion and uplift is most sensitive to parameters determining landslide volumes (particularly earthquake magnitude
M w , seismic energy source depth, and failure susceptibility, as well as the seismological factor responsible for triggering landslides), and is moderately sensitive to the effective elastic thickness of lithosphere,T e . For a specified magnitude, more erosive events (higher Ω) tend to occur at shallower depth, in thicker‐T e lithosphere, and in steeper, more landslide‐prone landscapes. For given landscape and seismotectonic conditions, the volumes of both landslides and uplift to first order positively scale withM w and seismic momentM o . However, higherM w earthquakes generate lower landslide and uplift volumes per unitM o , suggesting lower efficiency in themore » -
Abstract Interferometric Synthetic Aperture Radar is an important tool for imaging surface deformation from large continental earthquakes. Here, we present maps of coseismic displacement and strain from the 2019 Ridgecrest earthquakes using multiple Sentinel-1 images. We provide three types of interferometric products. (1) Standard interferograms from two look directions provide an overview of the deformation and can be used for modeling coseismic slip. (2) Phase gradient maps from stacks of coseismic interferograms provide high-resolution (∼30 m) images of strain concentration and surface fracturing that can be used to guide field surveys. (3) High-pass filtered, stacked, unwrapped phase is decomposed into east–west and up–down, south–north components and is used to determine the sense of fault slip. The resulting phase gradient maps reveal over 300 surface fractures, including triggered slip on the Garlock fault. The east–west component of high-pass filtered phase reveals the polarity of the strike-slip offset (right lateral or left lateral) for many of the fractures. We find a small number of fractures that have slip polarity that is retrograde to the background tectonic stress. This is similar to observations of retrograde slip observed near the 1999 Mw 7.1 Hector Mine rupture, but the Ridgecrest observations are more completely imaged by the frequent and high-qualitymore »
-
Abstract Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene M w > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.
-
Abstract We develop finite element models of the coseismic displacement field accounting for the 3D elastic structures surrounding the epicentral area of the 2019 Ridgecrest earthquake sequence containing two major events of Mw7.1 and Mw6.4. The coseismic slip distribution is inferred from the surface displacement field recorded by interferometric synthetic aperture radar. The rupture dip geometry is further optimized using a novel nonlinear‐crossover‐linear inversion approach. It is found that accounting for elastic heterogeneity and fault along‐strike curvilinearity improves the fit to the observed displacement field and yields a more accurate estimate of geodetic moment and Coulomb stress changes. We observe spatial correlations among the locations of aftershocks and patches of high slip, and rock anomalous elastic properties, suggesting that the shallow crust's elastic structures possibly controlled the Ridgecrest earthquake sequence. Most of the coseismic slip with a peak slip of 7.4 m at 3.6 km depth occurred above a zone of reduced
S ‐wave velocity and significant post‐Mw7.1 afterslip. This implies that viscous materials or fluid presence might have contributed to the low rupture velocity of the mainshock. Moreover, the zone of high slip on the northwest‐trending fault segment is laterally bounded by two aftershock clusters, whose location is characterized by intermediate rock rigidity.more »