Tsunamigenic earthquakes pose considerable risks, both economically and socially, yet earthquake and tsunami hazard assessments are typically conducted separately. Earthquakes associated with unexpected tsunamis, such as the 2018 Mw 7.5 strike-slip Sulawesi earthquake, emphasize the need to study the tsunami potential of active submarine faults in different tectonic settings. Here, we investigate physics-based scenarios combining simulations of 3D earthquake dynamic rupture and seismic wave propagation with tsunami generation and propagation. We present time-dependent modeling of one-way linked and 3D fully coupled earthquakes and tsunamis for the ∼ 100 km long Húsavík–Flatey Fault Zone (HFFZ) in North Iceland. Our analysis shows that the HFFZ has the potential to generate sizable tsunamis. The six dynamic rupture models sourcing our tsunami scenarios vary regarding hypocenter location, spatiotemporal evolution, fault slip, and fault structure complexity but coincide with historical earthquake magnitudes. Earthquake dynamic rupture scenarios on a less segmented fault system, particularly with a hypocenter location in the eastern part of the fault system, have a larger potential for local tsunami generation. Here, dynamically evolving large shallow fault slip (∼ 8 m), near-surface rake rotation (± 20∘), and significant coseismic vertical displacements of the local bathymetry (± 1 m) facilitate strike-slip faulting tsunami generation. We model tsunami crest to trough differences (total wave heights) of up to ∼ 0.9 m near the town Ólafsfjörður. In contrast, none of our scenarios endanger the town of Akureyri, which is shielded by multiple reflections within the narrow Eyjafjörður bay and by Hrísey island. We compare the modeled one-way linked tsunami waveforms with simulation results using a 3D fully coupled approach. We find good agreement in the tsunami arrival times and location of maximum tsunami heights. While seismic waves result in transient motions of the sea surface and affect the ocean response, they do not appear to contribute to tsunami generation. However, complex source effects arise in the fully coupled simulations, such as tsunami dispersion effects and the complex superposition of seismic and acoustic waves within the shallow continental shelf of North Iceland. We find that the vertical velocity amplitudes of near-source acoustic waves are unexpectedly high – larger than those corresponding to the actual tsunami – which may serve as a rapid indicator of surface dynamic rupture. Our results have important implications for understanding the tsunamigenic potential of strike-slip fault systems worldwide and the coseismic acoustic wave excitation during tsunami generation and may help to inform future tsunami early warning systems.
more »
« less
Seafloor Observations Eliminate a Landslide as the Source of the 1918 Puerto Rico Tsunami
The 11 October 1918 devastating tsunami in northwest Puerto Rico had been used as an example for earthquake-induced landslide tsunami hazard. Three pieces of evidence pointed to a landslide as the origin of the tsunami: the discovery of a large submarine landslide scar from bathymetry data collected by shipboard high-resolution multibeam sonar, reported breaks of submarine cable within the scar, and the fit of tsunami models to flooding observations. Newly processed seafloor imagery collected by remotely operated vehicle (ROV) show, however, pervasive Fe–Mn crust (patina) on the landslide walls and floor, indicating that the landslide scar is at least several hundred years old. C14 dates of sediment covering the landslide floor verify this interpretation. Although we have not searched the region systematically for an alternative tsunami source, we propose a possible source—a two-segment normal-fault rupture along the eastern wall of Mona rift. The proposed fault location matches the published normal faults with steep bathymetry and is close to the International Seismological Center–Global Earthquake Model catalog locations of the 1918 mainshock and aftershocks. The ROV observations further show fresh vertical slickensides and rock exposure along the proposed fault trace. Hydrodynamic models from an Mw 7.2 earthquake rupture along the eastern wall of the rift faithfully reproduce the reported tsunami amplitudes, polarities, and arrival times. Our analysis emphasizes the value of close-up observations and physical samples to augment remote sensing data in natural hazard studies
more »
« less
- Award ID(s):
- 2103713
- PAR ID:
- 10579332
- Publisher / Repository:
- Geoscience World
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- Volume:
- 113
- Issue:
- 1
- ISSN:
- 0037-1106
- Page Range / eLocation ID:
- 268 to 280
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Tsunamigenic earthquakes pose considerable risks, both economically and socially, yet earthquake and tsunami hazard assessments are typically conducted separately. Earthquakes associated with unexpected tsunamis, such as the 2018 Mw 7.5 strike-slip Sulawesi earthquake, emphasize the need to study the tsunami potential of active submarine faults in different tectonic settings. Here, we investigate physics-based scenarios combining simulations of 3D earthquake dynamic rupture and seismic wave propagation with tsunami generation and propagation. We present time-dependent modeling of one-way linked and 3D fully coupled earthquakes and tsunamis for the ∼ 100 km long Húsavík–Flatey Fault Zone (HFFZ) in North Iceland. Our analysis shows that the HFFZ has the potential to generate sizable tsunamis. The six dynamic rupture models sourcing our tsunami scenarios vary regarding hypocenter location, spatiotemporal evolution, fault slip, and fault structure complexity but coincide with historical earthquake magnitudes. Earthquake dynamic rupture scenarios on a less segmented fault system, particularly with a hypocenter location in the eastern part of the fault system, have a larger potential for local tsunami generation. Here, dynamically evolving large shallow fault slip (∼ 8 m), near-surface rake rotation (± 20∘), and significant coseismic vertical displacements of the local bathymetry (± 1 m) facilitate strike-slip faulting tsunami generation. We model tsunami crest to trough differences (total wave heights) of up to ∼ 0.9 m near the town Ólafsfjörður. In contrast, none of our scenarios endanger the town of Akureyri, which is shielded by multiple reflections within the narrow Eyjafjörður bay and by Hrísey island. We compare the modeled one-way linked tsunami waveforms with simulation results using a 3D fully coupled approach. We find good agreement in the tsunami arrival times and location of maximum tsunami heights. While seismic waves result in transient motions of the sea surface and affect the ocean response, they do not appear to contribute to tsunami generation. However, complex source effects arise in the fully coupled simulations, such as tsunami dispersion effects and the complex superposition of seismic and acoustic waves within the shallow continental shelf of North Iceland. We find that the vertical velocity amplitudes of near-source acoustic waves are unexpectedly high – larger than those corresponding to the actual tsunami – which may serve as a rapid indicator of surface dynamic rupture. Our results have important implications for understanding the tsunamigenic potential of strike-slip fault systems worldwide and the coseismic acoustic wave excitation during tsunami generation and may help to inform future tsunami early warning systems.more » « less
-
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
-
On 29 July 2025, an Mw 8.8 earthquake struck off Kamchatka, Russia, generating a Pacific-wide tsunami and marking the largest earthquake since the launch of the surface water and ocean topography (SWOT) satellite in 2022. We analyze tsunami observations from SWOT together with three nearby deep-ocean assessment and reporting of tsunamis (DART) buoys to resolve the source of the event. SWOT provided the first high-resolution spaceborne track of a great subduction-zone tsunami, capturing waveforms that reveal complex propagation, dispersion, and scattering. Inversion of the DART time series using Gaussian unit sources shows that the rupture extended ∼400 km along strike, with peak uplift of ∼4 m, significantly different from the published finite-fault model. A blended source that combines the DART-inverted uplift with subsidence from the seismic–geodetic model best matches both datasets and reproduces the SWOT observations. Comparison with reconstructions of the 1952 Mw 9.0 Kamchatka earthquake indicates that the 2025 rupture likely reactivated significant portions of the megathrust that broke in 1952 but occurred farther down-dip and with little to no near-trench slip, consistent with its smaller tsunami impact. These findings highlight the hazard implications of short recurrence intervals of great earthquakes and show how rupture style governs tsunami severity. They also demonstrate the value of satellite altimetry for improving tsunami source characterization, post-event forecasting, and understanding of hydrodynamic processes.more » « less
-
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.more » « less
An official website of the United States government

