On 15 January 2022 at 04:15 UTC, the Hunga Tonga-Hunga Ha’apai (HTHH) volcano in Tonga produced a massive eruption that triggered a transoceanic tsunami generated by the coupled ocean and atmospheric shock wave produced during the explosion. The tsunami first reached the coast of Tonga and eventually reached many coasts around the world. This volcano previously underwent a massive eruption in 1100 AD, and an eruption occurs approximately every 1000 years. The 2022 HTHH event provides an opportunity to study a major volcanically generated tsunami that caused substantial damage. In this study, we present a numerical simulation of a tsunami with a state-of-the-art numerical model based on a submarine explosion scenario. We constrain the geometry and magnitude of the explosion energy source based on analyses of pre- and post-event satellite images, which demonstrate that the explosion magnitude varied from 1 to 90 megatons of trinitrotoluene (Mt). Estimated submarine explosion geometries result in a suitable explosion magnitude of approximately 25 Mt, as determined with the waveform from the tide gauge in the time and frequency domains. The tsunami wave first reached the northwestern part of Tonga’s Tongatapu within 10 min, with a maximum runup height of approximately 15 m, and covered the whole of Tongatapu within 30 min. Finally, the numerical simulation provides deep insights into the physical volcanic explosion processes and improves our understanding and forecasting capabilities of frequent and catastrophic tsunamis caused by submarine volcanic explosions.

We developed tsunami fragility functions using three sources of damage data from the 2018 Sulawesi tsunami at Palu Bay in Indonesia obtained from (i) field survey data (FS), (ii) a visual interpretation of optical satellite images (VI), and (iii) a machine learning and remote sensing approach utilized on multisensor and multitemporal satellite images (MLRS). Tsunami fragility functions are cumulative distribution functions that express the probability of a structure reaching or exceeding a particular damage state in response to a specific tsunami intensity measure, in this case obtained from the interpolation of multiple surveyed points of tsunami flow depth. We observed that the FS approach led to a more consistent function than that of the VI and MLRS methods. In particular, an initial damage probability observed at zero inundation depth in the latter two methods revealed the effects of misclassifications on tsunami fragility functions derived from VI data; however, it also highlighted the remarkable advantages of MLRS methods. The reasons and insights used to overcome such limitations are discussed together with the pros and cons of each method. The results show that the tsunami damage observed in the 2018 Sulawesi event in Indonesia, expressed in the fragility function developed herein, is similar in shape to the function developed after the 1993 Hokkaido Nansei-oki tsunami, albeit with a slightly lower damage probability between zero-to-five-meter inundation depths. On the other hand, in comparison with the fragility function developed after the 2004 Indian Ocean tsunami in Banda Aceh, the characteristics of Palu structures exhibit higher fragility in response to tsunamis. The two-meter inundation depth exhibited nearly 20% probability of damage in the case of Banda Aceh, while the probability of damage was close to 70% at the same depth in Palu.

Indonesia has experienced several tsunamis triggered by seismic and non-seismic (i.e., landslides) sources. These events damaged or destroyed coastal buildings and infrastructure and caused considerable loss of life. Based on the Global Earthquake Model (GEM) guidelines, this study assesses the empirical tsunami fragility to the buildings inventory of the 2018 Sunda Strait, 2018 Sulawesi–Palu, and 2004 Indian Ocean (Khao Lak–Phuket, Thailand) tsunamis. Fragility curves represent the impact of tsunami characteristics on structural components and express the likelihood of a structure reaching or exceeding a damage state in response to a tsunami intensity measure. The Sunda Strait and Sulawesi–Palu tsunamis are uncommon events still poorly understood compared to the Indian Ocean tsunami (IOT), and their post-tsunami databases include only flow depth values. Using the TUNAMI two-layer model, we thus reproduce the flow depth, the flow velocity, and the hydrodynamic force of these two tsunamis for the first time. The flow depth is found to be the best descriptor of tsunami damage for both events. Accordingly, the building fragility curves for complete damage reveal that (i) in Khao Lak–Phuket, the buildings affected by the IOT sustained more damage than the Sunda Strait tsunami, characterized by shorter wave periods, and (ii) the buildings performed better in Khao Lak–Phuket than in Banda Aceh (Indonesia). Although the IOT affected both locations, ground motions were recorded in the city of Banda Aceh, and buildings could have been seismically damaged prior to the tsunami's arrival, and (iii) the buildings of Palu City exposed to the Sulawesi–Palu tsunami were more susceptible to complete damage than the ones affected by the IOT, in Banda Aceh, between 0 and 2 m flow depth. Similar to the Banda Aceh case, the Sulawesi–Palu tsunami load may not be the only cause of structural destruction. The buildings' susceptibility to tsunami damage in the waterfront of Palu City could have been enhanced by liquefaction events triggered by the 2018 Sulawesi earthquake.