Tsunami Waveforms Research Articles

Finite-fault Slip Distribution and Initial Water Elevation of the 2024 Noto Peninsula, Japan Earthquake Estimated from Tsunami Waveforms

Finite-fault Slip Distribution and Initial Water Elevation of the 2024 Noto Peninsula, Japan Earthquake Estimated from Tsunami Waveforms

Tsunami Early Warning Using High-Frequency Ocean Radar System in the Kii Channel, Japan

Abstract The tsunami of the 2011 Tohoku earthquake was monitored by the high-frequency (HF) ocean radar system in the Kii Channel, Japan. In this study, we conducted a retrospective analysis of tsunami early warning based on HF radar observations. First, we synthesized the surface current velocity using observation data of two radar land stations: Minato and Saikazaki. After vector synthesis, the time series of current velocity showed high consistency with the results of forward simulation. Then, we used the time series starting from 08:00 (UTC), 11 March 2011, as the input for data assimilation, reconstructing the tsunami wavefield in the Kii Channel. Comparison with actual observations at Kobe showed that our approach can accurately forecast coastal tsunami waveforms at least 50 min prior to tsunami arrival. The normalized prediction error decreased over time, indicating a progressively improving forecast accuracy. This is the first study on tsunami data assimilation using HF radar for the early warning of an actual seismogenic tsunami. We also discussed the limitations of the current HF radar system in the Kii Channel and made suggestions regarding the deployment of future radar land stations.

Nearshore Propagation and Amplification of the Tsunami Following the 2024 Noto Peninsula Earthquake, Japan

AbstractOn 1 January 2024, a massive earthquake struck the Noto Peninsula in Japan, and the resulting tsunami significantly damaged the nearby coasts. This study focuses on Iida Bay, located on the peninsula that suffered severe coastal damage. We observed a tsunami waveform in Iida Bay from video analysis owing to the lack of instrument data. This, combined with the waveform observed outside the bay, was used to correct the initial sea‐surface displacements produced by a prior source model. We successfully reproduced the observed waveforms with corrected displacements and investigated how the tsunami propagated and amplified in the bay. Alongshore wave propagation and superposition were identified as crucial factors in this tsunami event.

Ocean Surface Gravity Waves Excited by the 2022 Eruption of Hunga Tonga‐Hunga Ha'apai Volcano

AbstractOn 15 January 2022, a massive underwater eruption occurred at the Hunga Tonga‐Hunga Ha'apai Volcano. The plume reached the mesosphere, and the eruption excited a significant atmospheric Lamb wave, which forced the tsunami. The complicated tsunami waveforms due to ocean‐atmosphere coupling prevented inferring the force history of the excitation. To address this, we analyze ocean surface gravity waves (OSWs) from 15 to 40 mHz, which are decoupled from the Lamb wave due to their slower phase velocities. Modeling these OSWs, we infer that the excitation started at 4:00 UTC with an amplitude of N and lasted for 5 hr, followed by a sub‐event at 8:40 UTC. The observations suggest an initial blowout of seawater above the summit and a subsequent outflow that excited a tsunami below 5 mHz. The 2 hr delayed OSW excitation from 6 to 15 mHz may indicate seawater inflow into the crater.

Estimation of Tsunami Waveform at Wajima Port by the 2024 Noto Peninsula Earthquake using Video Images

Estimation of Tsunami Waveform at Wajima Port by the 2024 Noto Peninsula Earthquake using Video Images

Adjoint Synthesis for Trans‐Oceanic Tsunami Waveforms and Simultaneous Inversion of Fault Geometry and Slip Distribution

AbstractTsunamis propagate over long distances and can cause widespread damage even after crossing ocean basins. Prediction of tsunamis in distant areas based on observations near their sources is critical to mitigating damage. In recent years, the accuracy of numerical models of trans‐oceanic tsunami propagation has improved significantly due to the incorporation of effects such as the solid earth response to tsunami loading and wave dispersion. However, these models are computationally expensive and have not been fully utilized for real‐time prediction. Here, we derive the adjoint operator for the linear set of equations describing deep‐ocean tsunami propagation and show how a pre‐computed database of adjoint states can achieve rapid synthesis of tsunami waveforms at target sites from nonpoint arbitrary tsunami sources. The adjoint synthesis method allows for an exhaustive parameter search for tsunami source estimation. A method for simultaneous inversion of fault geometry and slip distribution using adjoint synthesis with Sequential Monte Carlo method was proposed and applied to the 2012 Haida Gwaii earthquake tsunami. The influence of model accuracy and the amount of observed data on the estimation of tsunami sources and waveforms was examined. It was found that with a highly accurate propagation model, using only a limited amount of observed data produced source and waveform estimates very similar to the final models obtained with much larger data sets. The final inferred fault model involved megathrust slip distributed between the Haida Gwaii trench and the Queen Charlotte fault. The proposed method can also quantify the uncertainty of the waveform forecasts.

Tsunami waveform inversion using Green’s functions with advection effects: application to the 2003 Tokachi–Oki earthquake

We explored nonlinear effects within the context of tsunami waveform inversion, wherein Green's functions were linearly superimposed to estimate earthquake slips. We focused on these effects while developing a source model for the 2003 Tokachi–Oki earthquake off Hokkaido, Japan. A source model for this earthquake was developed based on linear tsunami waveform inversion using Green’s functions and tsunami waveforms observed at tide gauge stations. Subsequently, tsunami waveforms from the source were simulated at the stations using nonlinear long-wave theory and compared with those estimated by inversion. The comparisons demonstrated that the waveforms had a non-negligible discrepancy that was attributed to advection effects, even for the primary wave used in the inversion at the two stations. This result strongly suggests that advection effects should be considered in the source modeling of the 2003 earthquake based on tsunami waveforms observed by tide gauges. Based on these results, a new tsunami waveform inversion technique that incorporates linearly approximated advection effects and maintain the framework of linear tsunami waveform inversion using Green’s functions is proposed and applied. The proposed method successfully mimicked the advection effects during the 2003 tsunami, reproduced better tsunami waveforms, and developed a source model for the 2003 earthquake using these effects. The peak slip amount and seismic moment were greater in the source model with advection effects than those without the effects. This finding suggests that the values in the source models developed for other earthquake events without considering these effects may have been underestimated.Graphical abstract

Physical Modelling Of Boulder Transport Under The Influence Of Tsunami Waves

Tsunami events are traditionally represented in the geological record by a sequence of fine-grained sediments, but increasingly coastal boulder deposits are being used as indicators of past tsunami events. The emplacement mechanism of many boulder deposits, however, is heavily debated and determining whether the inundation event was a tsunami or storm remains an unresolved challenge (Cox et al., 2020). Using physical experiments, we aim to achieve a better understanding of how tsunamis move coastal boulders. This knowledge will aid field geomorphologists in the identification of the emplacement mechanism for coastal boulder deposits and allow for the determination of wave parameters. In January 2023, physical experiments using the HR Wallingford Tsunami Simulator were completed as part of the MAKEWAVES collaboration. These experiments investigated the movement of a cuboid and irregular shaped boulder model when impacted by different tsunami waveforms on a plane beach. We propose new empirical formulae to describe relationships between transport distance and different tsunami waves.

Modeling the 2022 Tonga Eruption Tsunami Recorded on Ocean Bottom Pressure and Tide Gauges Around the Pacific

AbstractTsunamis generated by the Hunga Tonga–Hunga Ha’apai volcanic eruption on January 15, 2022 were recorded on ocean bottom pressure and tide gauges around the Pacific Ocean, earlier than the expected arrival times calculated by tsunami propagation speed. Atmospheric waves from the eruption were also recorded globally with propagation speeds of ~ 310 m/s (Lamb wave) and 200–250 m/s (Pekeris wave). Previous studies have suggested that these propagating atmospheric waves caused at least the initial part of the observed tsunami. We simulated the tsunamis generated by the propagation of the Lamb and Pekeris waves by adding concentric atmospheric pressure changes. The concentric sources are parameterized by their propagation speeds, initial atmospheric wave amplitudes that decay with the distance, and a rise time. For the Lamb wave, inversions of the observed tsunami waveforms at 14 U.S. and nine New Zealand DART stations indicate the start of the positive rise at 4:16 UTC, the peak amplitude of 383 hPa, and the propagation speed of 310 m/s, assuming a rise time of 10 min. The later phases of the observed tsunami waveforms can be better reproduced by adding another propagating concentric wave (Pekeris wave) with a negative amplitude (− 50 hPa) and propagation speeds of 200–250 m/s. The DART records around the Pacific indicate that the Pekeris wave speed is faster toward the northwest and slightly slower toward the northeast. The synthetic waveforms roughly reproduced the far-field tsunami waveforms recorded at tide gauge stations, including the later phases, suggesting that the large amplitude in the later phase may be due to the coupling of the Pekeris wave and the tsunami, as well as resonance around tide gauge stations.

Slip distribution of the 2024 Noto Peninsula earthquake (MJMA 7.6) estimated from tsunami waveforms and GNSS data

AbstractThe 1 January 2024 Noto-Hanto (Noto Peninsula) earthquake (MJMA 7.6) generated strong ground motion, large crustal deformation and tsunamis that caused significant damage in the region. Around Noto Peninsula, both offshore submarine and partially inland active faults have been identified by previous projects: Ministry of Land, Infrastructure, Transport and Tourism (MLIT) and Japan Sea Earthquake and Tsunami Research Project (JSPJ). We inverted the tsunami waveforms recorded on 6 wave gauges and 12 tide gauges around Sea of Japan and the GNSS data recorded at 53 stations in Noto Peninsula to estimate the slip amount and seismic moment on each of active faults. The results show that the 2024 coseismic slips were 3.5 m, 3.2 m, and 3.2 m on subfaults NT4, NT5 and NT6 of the JSPJ model, located on the northern coast of Noto Peninsula and dipping toward southeast. A smaller slip, 1.0 m, estimated on NT8 on the southwestern end of the 2024 rupture, may be attributed to its previous rupture during the 2007 Noto earthquake. The total length of these four faults is ~ 100 km, and the seismic moment is 1.90 × 1020 Nm (Mw = 7.5). Almost no slip was estimated on the northeastern subfaults NT2 and NT3, which dip northwestward, opposite to NT4–NT5–NT6, and western subfault NT8. Aftershocks including the MJMA 6.1 event occurred in the NT2–NT3 region, but they are smaller than the potential magnitude (Mw 7.1) those faults can release in a tsunamigenic earthquake. Similar features are also found for the MLIT model; the 2024 slip was only on F43 along the northern coast of Noto Peninsula, and northeastern F42 did not rupture, leaving potential for future event. Graphical Abstract

Volcano generated tsunami recorded in the near source

Volcano sector collapse and pyroclastic density currents are common phenomena on active volcanoes and potentially a fatal source of tsunami waves which constitute a serious hazard for local as well as distant coastal population. Several examples in recent history, warn us on the urgent need to improve our mitigation counter-actions when tsunamis have volcanic origin. However, instrumental record of tsunami generated by mass movement along a volcano flank are still rare and not well understood yet. Small tsunamis (≤1 m) induced by pyroclastic density currents associated to violent explosions of Stromboli volcano were recorded in near-source conditions (<1.6 km). We show how tsunami waveform remains unaltered regardless of the two orders of variability in the landslide volume and dynamics. This unprecedented record is also providing the lesson to develop unconventional warning strategies necessary when the tsunamigenic source is expected to be very close (<10 minutes) to densely populated coasts and with a limited time to issue an alert based on simulation of wave propagation and inundation.

Linked and fully coupled 3D earthquake dynamic rupture and tsunami modeling for the Húsavík–Flatey Fault Zone in North Iceland

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.

Enigmatic Tsunami Waves Amplified by Repetitive Source Events Near Sofugan Volcano, Japan

AbstractOn 8 October 2023, mysterious tsunamis with a maximum wave height of 60 cm were observed in Izu Islands and southwestern Japan, although only seismic events with body‐wave magnitudes mb 4–5 have been documented to the west of Sofugan volcano. To investigate the source process, we analyze tsunami waveforms recorded by an array network of ocean bottom pressure gauges. Stacked waveforms of pressure gauge records suggest recurrent arrivals of multiple wave trains. Deconvolution of the stacked waveforms by tsunami waveforms from an earlier event revealed over 10 source events that intermittently generated tsunamis for ∼1.5 hr. The temporal history of this sequence corresponds to the origin times of T‐phases estimated by an ocean bottom seismometer and of the seismic swarm, implying a common origin. Larger events later in the sequence occurred at intervals comparable to the tsunami wave period, causing amplification of later phases of the tsunami waves.

Tsunami Waveform Characteristics of the Korean Peninsula

Tsunami Waveform Characteristics of the Korean Peninsula

Quantifying Magma Overpressure Beneath a Submarine Caldera: A Mechanical Modeling Approach to Tsunamigenic Trapdoor Faulting Near Kita‐Ioto Island, Japan

AbstractSubmarine volcano monitoring is vital for assessing volcanic hazards but challenging in remote and inaccessible environments. In the vicinity of Kita‐Ioto Island, south of Japan, unusual M ∼ 5 non‐double‐couple volcanic earthquakes exhibited quasi‐regular recurrence near a submarine caldera. Following the earthquakes in 2008 and 2015, a distant ocean bottom pressure sensor recorded distinct tsunami signals. In this study, we aim to find a source model of the tsunami‐generating earthquake and quantify the pre‐seismic magma overpressure within the caldera's magma reservoir. Based on the earthquake's characteristic focal mechanism and efficient tsunami generation, we hypothesize that submarine trapdoor faulting occurred due to highly pressurized magma. To investigate this hypothesis, we establish mechanical earthquake models that link pre‐seismic magma overpressure to the size of the resulting trapdoor faulting, by considering stress interaction between a ring‐fault system and a reservoir of the caldera. The trapdoor faulting with large fault slip due to magma‐induced shear stress in the submarine caldera reproduces well the observed tsunami waveform. Due to limited data, uncertainties in the fault geometry persist, leading to variations of magma overpressure estimation: the pre‐seismic magma overpressure ranging approximately from 5 to 20 MPa, and the co‐seismic pressure drop ratio from 10% to 40%. Although better constraints on the fault geometry are required for robust magma pressure quantification, this study shows that magmatic systems beneath calderas are influenced significantly by intra‐caldera fault systems and that tsunamigenic trapdoor faulting provides rare opportunities to obtain quantitative insights into remote submarine volcanism hidden under the ocean.

Estimation of the Seismic Source of the 1974 Lima Peru Earthquake and Tsunami (Mw 8.1)

In this investigation, we have conducted a long period teleseismic and tsunami waveform inversion to obtain the slip distribution of the 1974 Lima-Perú earthquake occurred in the central region of Peru. According to teleseismic inversion, the rupture process was complex with a duration of 90 s approximately and the main asperity was located in the northern side of the rupture geometry, offshore Lima and Callao. According to tsunami waveform inversion the main asperity was located in the southern side of the fault plane, offshore Cañete. However, a joint inversion of teleseismic and tsunami waveforms averaged the contribution of the two datasets and could give a better result. The maximum slip from the joint inversion was 7.25 m, and the moment magnitude was Mw 8.1. Despite the occurrence of this earthquake (almost 50 years ago), there is a high seismic potential to trigger a tsunamigenic earthquake in the central region of Peru. According to the interseismic coupling the next earthquake in the central region of Peru is ready to occur and it would be of 8.5–8.8 Mw.

Source Characteristics of the 18 November 2022 Earthquake ($${M}_{W}$$ 6.7), Offshore Southwest Sumatra, Indonesia, Revealed by Tsunami Waveform Analysis: Implications for Tsunami Hazard Assessment

Source Characteristics of the 18 November 2022 Earthquake ($${M}_{W}$$ 6.7), Offshore Southwest Sumatra, Indonesia, Revealed by Tsunami Waveform Analysis: Implications for Tsunami Hazard Assessment

Origins of the Tsunami Following the 2023 Turkey–Syria Earthquake

AbstractOn 6 February 2023, a local tsunami was recorded in the southeastern Mediterranean Sea following the Mw 7.7 Turkey–Syria inland strike‐slip earthquake. Due to the lack of underwater observation, the tsunami generation mechanism remains mysterious. To understand the source mechanisms, we analyzed the tsunami waveforms of four nearby tide gauges and located possible sources using a backward tsunami ray tracing approach. We then conducted forward numerical modelings for a range of possible source parameters. We show that there were probably two tsunami sources, inside and outside Iskenderun Bay, which may be related to thick coastal sediments. A source inside the Bay with a characteristic length of 7 km produced dominant periods of 10–30 min with negative initial motion, possibly generated by a landslide. Another source of 6 km length outside the Bay produced dominant periods of 2–10 min with positive initial motions, possibly related with liquefaction.

Multi-scale Simulation of Subsequent Tsunami Waves in Japan Excited by Air Pressure Waves Due to the 2022 Tonga Volcanic Eruption

The 2022 Hunga Tonga-Hunga Ha’apai eruption generated tsunamis that propagated across the Pacific Ocean. Along the coast of Japan, nearshore amplification led to amplitudes of nearly 1 m at some locations, with varying peak tsunami occurrence times. The leading tsunami wave can generally be reproduced by Lamb waves, which are a type of air-pressure wave generated by an eruption. However, subsequent tsunamis that occurred several hours after the leading wave tended to be larger for unknown reasons. This study performs multi-scale numerical simulations to investigate subsequent tsunami waves in the vicinity of Japan induced by air pressure waves caused by the eruption. The atmospheric pressure field was created using a dispersion relation of atmospheric gravity wave and tuned by physical parameters based on observational records. The tsunami simulations used the adaptive mesh refinement method, incorporating detailed bathymetry and topography to solve the tsunami at various spatial scales. The simulations effectively reproduced the tsunami waveforms observed at numerous coastal locations, and results indicate that the factors contributing to the maximum tsunami amplitude differ by region. In particular, bay resonance plays a major role in determining the maximum amplitude at many sites along the east coast of Japan. However, large tsunami amplification at some west coast locations was not replicated, probably because it was caused by amplification during oceanic wave propagation rather than meteorological factors. These findings enhance our understanding of meteotsunami complexity and help distinguish tsunami amplification factors.

Sensitivity of slip distribution on tsunami trace heights and geological evidences: a case study of the 2011 Tohoku-Oki earthquake

We examined whether it is possible to estimate the tsunami source model of the 2011 Tohoku-Oki earthquake from a comparison of numerical simulations of tsunami propagation and sediment transport, the measured trace heights, and the sediment thickness of tsunami deposits. Twelve models with different subfault numbers were prepared based on a reference model inferred from tsunami waveform inversion. The reference model with 55 subfaults considering rupture propagation and the model with instantaneous slip successfully reproduced both the tsunami trace heights and sediment thickness distribution of tsunami deposits in the Idagawa Lowland and Sendai Plain. Other models with the same moment magnitude but fewer subfaults could not reproduce the observed trace heights, and the reproducibility of sediment thickness distribution strongly depended on the slip distribution. Models with increased slip amounts and moment magnitude could reproduce the trace heights; however, the simulated sediment thickness was underestimated for the Idagawa Lowland while overestimated for the Sendai Plain. Our results indicate that the combination of trace heights, sediment thickness of tsunami deposits, and numerical simulations of tsunami propagation and sediment transport can be used to estimate historical earthquakes and tsunamis. Efforts should be made to increase the number of subfaults for the study of historical events, although the obtained solutions may not be unique because of fewer trace heights or tsunami deposit data.