An official website of the United States government
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
Comparison of methods for coupled earthquake and tsunami modelling
SUMMARY Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments.
more »
« less
- PAR ID:
- 10447515
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 234
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 404 to 426
- 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
-
Abstract One of the essential characteristics of earthquakes associated with large tsunami generation is the depletion of high‐frequency radiation, which is not well understood by elastic dislocation theory and largely not accounted for in most rupture models of real events. We present fully coupled models of dynamic rupture, ocean acoustic waves, and tsunami for the 1896 Sanriku earthquake with wedge inelasticity. The inelastic wedge deformation due to thick sediment in the northern Japan Trench is shown to generate efficient short‐wavelength seafloor uplift (>5 m), which is several times larger than the uplift by elastic dislocation models and generates impulsive tsunami that can have a large impact on the rugged Sanriku coast. Seismic moment due to inelastic wedge deformation has a reverse faulting focal mechanism with a steep plunge (>75°) of T axis, reflecting high efficiency in tsunami generation. However, the inelastic deformation is a large energy sink, which causes slower rupture velocity, weaker radiation of ocean acoustic and seismic waves, and ∼10 times lower moment‐scaled radiated energy than those of elastic models, explaining nearly all the anomalous characteristics of this tsunami earthquake. The anti‐plane off‐fault shear stress in the mode III rupture direction, limited by yielding, plays an important role in the slow rupture velocity and energy radiation along strike. Ocean acoustic waves may not provide robust signals for tsunami early warning due to weak high‐frequency radiation. Additionally, large, long‐duration ground velocity pulses can naturally result from inelastic deformation.more » « less
-
Tsunamis generated by seafloor displacements accompanying large submarine earthquakes provide sensitivity to absolute slip position and distribution for offshore faulting analogous to that of geodetic observations for landward faulting. Tsunami recordings at deep‐water and near‐shore ocean bottom pressure sensors and tide gauges, along with runup and inundation measurements, can now be reliably modeled using detailed bathymetric structures and robust numerical codes. As a result, tsunami observations now play an important role in quantifying coseismic slip distributions for large submarine earthquakes in subduction zones and other tectonic environments. Applications of joint modeling or inversion of seismic, geodetic and tsunami observations for recent major earthquakes are described, highlighting the specific contributions of the tsunami observations to source model resolution. Tsunami observations provide unique information on the up‐dip extent of earthquake coseismic slip on subduction zone megathrust faults and occurrence of near‐trench slip, which are usually not well constrained by seismic and land‐based geodetic signals. Tsunami signals also help to detect offshore slow slip that is not evident in seismic or land‐based geodetic data and to balance geophysical constraints on ruptures that extend from on‐shore to off‐shore. Tsunami runup measurements and stratigraphic deposits further provide unique constraints on large earthquake ruptures that occurred prior to modern geophysical instrumentation.more » « less
-
Tsunami wave observations far from the coast remain challenging due tothe logistics and cost of deploying and operating offshoreinstrumentation on a long-term basis with sufficient spatial coverageand density. Distributed Acoustic Sensing (DAS) on submarine fiber opticcables now enables real-time seafloor strain observations over distancesexceeding 100 km at a relatively low cost. Here, we evaluate thepotential contribution of DAS to tsunami warning by assessingtheoretically the sensitivity required by a DAS instrument to recordtsunami waves. Our analysis includes signals due to two effects induced by thehydrostatic pressure perturbations arising from tsunami waves: thePoisson’s effect of the submarine cable and the compliance effect of theseafloor. It also includes the effect of seafloor shear stresses andtemperature transients induced by the horizontal fluid flow associatedwith tsunami waves. The analysis is supported by fully coupled 3-Dphysics-based simulations of earthquake rupture, seismo-acoustic wavesand tsunami wave propagation. The strains from seismo-acoustic waves andstatic deformation near the earthquake source are orders of magnitudelarger than the tsunami strain signal. We illustrate a data processingprocedure to discern the tsunami signal. With enhanced low-frequencysensitivity on DAS interrogators (strain sensitivity ≈2×10 at mHz frequencies), we find that, on seafloorcables located above or near the earthquake source area, tsunamis areexpected to be observable with a sufficient signal-to-noise ratio withina few minutes of the earthquake onset. These encouraging results pavethe way towards faster tsunami warning enabled by seafloor DAS.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1093/gji/ggad053
