Search

Ischia Island has been repeatedly affected by mass collapses, which are mainly caused by the steepness of the main peak (Mt Epomeo) and by phenomena related to its volcanic activity. The most relevant cases of mass failure studied in the literature and postulated to be tsunamigenic cover a wide spectrum of sizes, from sector collapse to small-volume mass transports. Tsunamis generated by landslides in Ischia may affect the coast of the Campania mainland, including the Gulf of Naples. The focus of this work is an evaluation of the pattern of the maximum tsunami energy. To this purpose, we perform a series of numerical simulations by moving the same landslide source in different hypothetical positions around the island. The landslide dynamics are computed through the code UBO-BLOCK, and the tsunami propagation by employing the code UBO-TSUFD, both developed in-house. The final goal is to characterize the coastal areas of the Campania mainland most exposed to tsunami attack from Ischia sources. It is found that the position of the landslide deeply influences the distribution of the tsunami elevation in the coastal stretch north of the Procida Mt, while, remarkably, it is irrelevant inside the Gulf of Naples where the bathymetric effect prevails.

In this research, we have conducted a numerical simulation to obtain the seismic source, coseismic deformation field, and the tsunami propagation of the great Lima-Peru 1940 earthquake and tsunami, based on macroseismic information and focal mechanism from the literature. The seismic dimensions of our preferred model were set at 162 km×71 km according to VIII isoseismal intensity. The slip distribution is homogeneous with a mean value of 2.7 m. The fault plane orientation was set at strike = 330∘, dip = 20∘, and rake = 90∘. The maximum simulated uplift was 1.27 m and the maximum subsidence was 0.36 m. Due to the proximity of the seismic source to the coast, the city of Huacho was uplifted around 11 cm. The fault plane of the 1940 earthquake was located in the updip near the coast. The maximum tsunami height simulated in Huacho was 1.80 m. We suggest that there is a high tsunamigenic potential near the trench to generate a shallow earthquake.

Geologically, Sulawesi Island is a complex area. It is caused by the convergence between three lithosphere plates, where the Australian Plate drifting to North, Pacific to West, and Eurasia from South to South-East respective. The future tsunami simulation was done by using SiTProS (Siam Tsunami Propagation Simulator) Version 1.5., which uses data such as earthquakes coordinate, magnitudes, focus depths, and propagation scheme, as the main input parameters in the simulations. Based on those future tsunami simulations results, the mitigation scheme chosen for Banggae District of Majene Region, West Sulawesi Province is the Vertical Tsunami Mitigation scheme. Due to its detailed parameters i.e. population details, evacuation time, distance and speed, the Vertical Tsunami Mitigation Scheme can determine the existence of shelters parameters and analyze the additional shelters needed in precise manner.

The 1783 Scilla landslide–tsunami (Calabria, southern Italy) is a well-studied event that caused more than 1500 fatalities on the beaches close to the town. This paper complements a previous work that was based on numerical simulations and was focused on the very local effects of the tsunami in Scilla. In this study we extend the computational domain to cover a wider portion of western Calabria and northeastern Sicily, including the western side of the Straits of Messina. This investigation focuses on Capo Peloro area (the easternmost cape of Sicily), where the highest tsunami effects outside Scilla were reported. Important tsunami observations, such as the wave height reaching 6 m at Torre degli Inglesi and flooding that reached over 600 m inland, have been successfully modeled but only by means of a high-resolution (10 m) topo-bathymetric grid, since coarser grids were inadequate for the purpose. Interestingly, the inundation of the small lake of Pantano Piccolo could not be reproduced by using today's coastal morphology, since a coastal dune now acts as a barrier against tsunamis. Historical analysis suggests that this dune was not in place at the time of the tsunami occurred and that a ground depression extending from the lake to the northern coast is a remnant of an ancient channel that was used as a pathway in Roman times. The removal of such an obstacle and the remodeling of the coeval morphology allows the simulations to reproduce the tsunami penetration up to the lake, thus supporting the hypothesis that the 1783 tsunami entered the lake following the Roman channel track. A further result of this study is that the computed regional tsunami propagation pattern provides a useful hint for assessing tsunami hazards in the Straits of Messina area, which is one of the most exposed areas to tsunami threats in Italy and in the Mediterranean Sea overall.

The frequent occurrence of landslide tsunami has aroused great social concern, hence an efficient approach for modeling the phenomena is required urgently. Herein, we present and supplement a meshless numerical model named ‘Tsunami Squares’ (TS). TS simulates both the landslide and its generated wave using an updated numerical scheme that has enhanced computational efficiency. TS also provides an artful way of approximating linear dispersive wave propagation and non-linear wave inundation that makes it quite capable of simulating real-scale events. A well-designed test that compares TS output with an analytical formula for a tsunami propagation calculation provides firm validation of the approach. We then applied the ‘Tsunami Squares’ to simulate the 1792 Unzen-Mayuyama mega slide and generated tsunami. This event remains the largest volcanic disaster in the history of Japan. Landslide dynamics simulated by TS coincide with historical reports and the known geometry of landslide deposits. The model output of tsunami runup heights well match the observed values recorded by the locations of historical ‘Tsunami Stones’. This well documented and data-rich case further validates the ‘Tsunami Squares’ model and the excellent reproduction of the 1792 event provides new insights for better understanding of landslide tsunami hazards.

Earthquake rupture dynamic models capture the variability of slip in space and time while accounting for the structural complexity which is characteristic for subduction zones. The use of a geodynamic subduction and seismic cycling (SC) model to initialize dynamic rupture (DR) ensures that initial conditions are self-consistent and reflect long-term behavior. We extend the 2D geodynamical subduction and SC model of van Zelst et al. (2019) and use it as input for realistic 3-dimensional DR megathrust earthquake models. We follow the subduction to tsunami run-up linking approach described in Madden et al. (2020), including a complex subduction setup along with their resulting tsunamis. The distinct variation of shear traction and friction coefficients with depth lead to realistic average rupture speeds and dynamic stress drop as well as efficient tsunami generation. We here analyze a total of 14 subduction-initialized 3D dynamic rupture-tsunami scenarios. By varying the hypocentral location along arc and depth, we generate 12 distinct unilateral and bilateral earthquakes with depth-variable slip distribution and directivity, bimaterial, and geometrical effects in the dynamic slip evolutions. While depth variations of the hypocenters barely influence the tsunami behavior, lateral varying nucleation locations lead to a shift in the on-fault slip which causes time delays of the wave arrival at the coast. The fault geometry of our DR model that arises during the SC model is non-planar and includes large-scale roughness. These features (topographic highs) trigger supershear rupture propagation in up-dip or down-dip direction, depending on the hypocentral depth.In two additional scenarios, we analyze variations in the energy budget of the DR scenarios. In the SC model, an incompressible medium is assumed (ν=0.5) which is valid only for small changes in pressure and temperature. Unlike in the DR model where the material is compressible and a different Poisson’s ratio (ν=0.25) has to be assigned. Poisson’s ratios between 0.1 and 0.4 stand for various compressible materials. To achieve a lower shear strength of all material on and off the megathrust fault and to facilitate slip, we increase the Poisson ratio in the DR model to ν=0.3 which is consistent with basaltic rocks. As a result, larger fault slip is concentrated at shallower depths and generates higher vertical seafloor displacement and earthquake moment magnitude respectively. Even though the tsunami amplitudes are much higher, the same dynamic behavior as in the twelve hypocenter-variable models can be observed. Lastly, we increase fracture energy by changing the critical slip distance in the linear slip-weakening frictional parameterization. This generates a tsunami earthquake (Kanamori, 1972) characterized by low rupture velocity (on average half the amount of s-wave speed) and low peak slip rate, but at the same time large shallow fault slip and moment magnitude. The shallow fault slip of this event causes the highest vertical seafloor uplift compared to all other simulations. This leads to the highest tsunami amplitude and inundation area while the wavefront hits the coast delayed compared to the other scenarios.

International audience; We propose a numerical method for a family of two-dimensional dispersive shallow water systems with topography. The considered models consist in shallow water approximations without the hydrostatic assumption-of the incompressible Euler system with free surface. Hence, the studied models appear as extensions of the classical shallow water system enriched with dispersive terms. The model formulation motivates to use a prediction-correction scheme for its numerical approximation. The prediction part leads to solving a classical shallow water system with topography while the correction part leads to solving an elliptic-type problem. The numerical approximation of the considered dispersive models in the two-dimensional case over unstructured meshes is described, it requires to combine finite volume and finite element techniques. A special emphasis is given to the formulation and the numerical resolution of the correction step (variational formulation, inf-sup condition, boundary conditions,.. .). The numerical procedure is confronted with analytical and experimental test cases. Finally, an application to a real tsunami case is given.

The paper deals with the analysis of the propagation, heights, and arrival times of a tsunami that may occur on the coastal areas of Cyprus and Eastern Mediterranean in the case of an earthquake in southern Cyprus. Following a review of the seismic risk and historical earthquakes which occurred in southern Cyprus, it was concluded that this region may be subject to high vulnerability if a tsunami occurs. A study was conducted on the numerical modelling of a possible tsunami generated by movement along the fault of the 1222 Paphos earthquake. The region where the earthquake occurred can be attributed to the Cyprian Arc in the southwest of Cyprus. This arc is one of the most active seismic zones in the Mediterranean, which has led to the occurrence of earthquakes from submerged seismogenic sources. A methodology is used to simulate tsunami wave propagation for the Eastern Mediterranean coastal areas which requires the initial wave due to fault parameters as well as the bathymetry data. The GEBCO30 bathymetry data are used which have a grid spacing of 0.30 arc min. The fault parameters are deduced from the maximum stress directions and source geometry of the region from the moment tensor solutions derived from analyzing of earthquake waveforms. The numerical tsunami propagation model was performed by using SWAN code. The simulated highest tsunami heights were 4.02 m in Kouklia (Cyprus); 2.85 m in Paphos Ktima (Cyprus); 2.58 in Episkopi (Cyprus); 2.06 in Peyia (Cyprus); 1.76 in Yennadhi, Rhodes (Greece); 1.53 in Burg Migheizil (Egypt); 1.46 m in Tarabulus (Lebanon); 1.39 m in Bur Said (Egypt); 1.28 in Al-Burj (Egypt); and 0.60 in Mugla-Aksaz (Turkey). The results of the model outline the extent of the tsunami waves of damaging size, but destructive event in the region.

Tsunami propagating into a coastal river can severely damage hydraulic structures including bridge piers as reported during past tsunami events. Hence, in the present study, laboratory experiments are performed initially by generating tsunami bores on to a wet bed in a 24m long, 0.8m wide wave flume and later in a 12mx0.5mx0.465m recirculating channel to maintain a time scale of 1:6 required as per Froude scaling to observe reliable scour depth at the pier model. Three-dimensional flow velocities (using rigid bed) and scour patterns, depths (using mobile bed) around a circular bridge pier model were obtained to understand the scour process. Finally, the scour depth due to tsunami bore is compared to empirical steady flow scour estimations. Therefore, the present study proposes new scour design standards for bridge piers or vertical evacuation towers which are already inundated by the first wave of tsunami.

This study aims to simulate tsunami wave heights and their travel time toward the WestKalimantan coastal area. A Nuclear Power Plant (NPP) will be developed in a coastal area in the WestKalimantan Province, Indonesia. The two scenarios of events that may trigger the tsunami waves are volcanicevents caused by Mount Krakatau located in the Sunda Strait and tectonic subduction located at the IndianOcean on the west side of the Sunda Strait. Modeling the tsunami propagation using a finite-element-basedhydrodynamic model developed by the US Army Corps of Engineer, namely the Surface-water ModelingSystem (SMS), is carried out. The model simulates tsunami wave propagation for six alternative locations atthe proposed site of NPP. The model domains consist of Sunda Strait and Karimata Strait domains. The SundaStrait domain model is validated by the observed historical tsunami heights reaching the 12 locations inLampung, Banten, West Java, and Jakarta of Indonesia. The model validations show a good agreement. Thevalidated hydrodynamic model results of the Sunda Strait model are used as the boundary conditions for theKarimata Strait domain model. The Karimata Strait model results show that the peaks of tsunami wave heightsthat reach the western Kalimantan coast are between 12 cm and 116 cm. The minimum wave height peak thatreaches the six prospective locations is 12 cm, which is located in Sambas Regency.

We propose a two-layer model with two different axes of integration and a well-balanced finite volume method. The purpose is to study submarine avalanches and generated tsunamis by a depth-averaged model with different averaged directions for the fluid and the granular layers. Two-layer shallow depth-averaged models usually consider either Cartesian or local coordinates for both layers. However, the motion characteristics of the granular layer and the water wave are different: the granular flow velocity is mainly oriented downslope while water motion related to tsunami wave propagation is mostly horizontal. As a result, the shallow approximation and depth-averaging have to be imposed (i) in the direction normal to the topography for the granular flow and (ii) in the vertical direction for the water layer. To deal with this problem, we define a reference plane related to topography variations and use the associated local coordinates to derive the granular layer equations whereas Cartesian coordinates are used for the fluid layer. Depthaveraging is done orthogonally to that reference plane for the granular layer equations and in the vertical direction for the fluid layer equations. Then, a finite volume method is defined based on an extension of the hydrostatic reconstruction. The proposed method is exactly well-balanced for two kind of stationary solutions: the classical one, when both water and granular masses are at rest; the second one, when only the granular mass is at rest. Several tests are presented to get insight into the sensitivity of the granular flow, deposit and generated water waves to the choice of the coordinate systems. Our results show that even for moderate slopes (up to 30◦), strong relative errors on the avalanche dynamics and deposit (up to 60%) and on the generated water waves (up to 120%) are made when using Cartesian coordinates for both layers instead of an appropriate local coordinate system as proposed here. Ministerio de Economía y Competitividad (MINECO). España European Commission (EC). Fondo Europeo de Desarrollo Regional (FEDER) Agence Nationale de la Recherche. France European Research Council (ERC)

A numerical model that simulates the propagation of tsunamis produced by submarine earthquakes was applied to the Caspian Sea. The model was first applied to the June 1990 (Rudbar) earthquake and our results proved to be coherent with previous simulations. It was also applied to the Krasnovodsk earthquake (1895), which had not been simulated before, and it was found that the coastal area was flooded by a tsunami, as reported in literature. Nevertheless, it was a rather local effect since the tsunami was generated in very shallow waters. Some worst-case hypothetical earthquakes were simulated in the most seismically active regions of the Caspian. In general, both the east and west coasts of the central Caspian Sea would be the most affected by these earthquake-induced tsunamis, with potential significant effects in some cases in cities such as Baku.

Tsunami overland flow may cause severe destruction to the affected coastal area. Several tsunamis had hit Indonesia. Studies have shown that there is a high probability of a massive earthquake in the Mentawai fault, Sumatra, Indonesia. A massive earthquake with a magnitude of over 8 Mw might occur in the next 50 years. The nearshore tsunami propagation and inundation in Painan City, West Sumatra, Indonesia, is assessed using HEC-RAS. Various land cover types are simulated to investigate their effects on the tsunami overland flow. The model showed satisfactory results. The model can show the known behaviors of tsunami propagation in urban areas and in rivers. It was found that there is a significant effect of land cover types on the inundation at the coast. Inundation depth increases with roughness along the coast. On the other hand, high roughness decreases the inundated area. Furthermore, the land cover has a significant effect on the velocity. The velocity is much lower in higher Manning’s value. Simulation with a high surface roughness shows a significant reduction in both inundation area and velocity. The high roughness causes a blockage effect for tsunami propagation. Thus, the Tsunami cannot penetrate further inland. The information can be used for land use control suggestions. These provide significant information for future mitigation plans, e.g., placement of evacuation buildings and evacuation routes.

Shallow water equations were analyzed numerically for simulation of tsunami propagation from the source area to coastal zone of Mentawai Islands along West Sumatra Island. A triangular mesh was set for spatial discretization and a system of partial differential equation is reduced to a system of ordinary differential equations. Initial displacement based on Okada model was applied. Our numerical techniques were demonstrated with a simulation of tsunami generated from an earthquake with epicenter near South Pagai Island, Mentawai Islands, Indonesia.

In the present study, the local tsunami amplification observed in Ryori Bay, located on the Sanriku coast of Japan, was investigated using numerical simulations. Large-scale tsunami propagation simulations and tsunami inundation simulations for the bay were systematically conducted to estimate and model the 2011, 1933, and 1896 tsunamis that occurred off the Sanriku coast and which resulted in large run-ups. The simulation results, which are moderately consistent with observations, presented larger run-up heights and inundations for the 1933 and 1896 tsunamis (which followed relatively small earthquakes) compared to those of the 2011 tsunami (which followed a larger earthquake). Furthermore, the frequency analysis indicated that the former two tsunamis comprised higher predominant components. A tsunami inundation simulation using parametrized synthetic waveforms was conducted to identify the contributing factors associated with the large amplification and run-ups. The results indicated that the predominant components are significantly amplified in the bay and the initial decrease in the water surface elevation prior to the primary waves of the two tsunamis leads to an increase in their run-up heights. Furthermore, the simulated waveforms of the tsunamis revealed that the 1933 and 1896 tsunamis had their wavefronts changed into a steep wavefront, i.e., a bore-like wave, during their wave developments in the bay, attributed to shoaling, narrowing bay width, and the nonlinear effect of the wave. These results, therefore, indicate that bores which are known to generate large run-up heights were generated in the bay during the two tsunamis.

Tsunami propagation over the entire Dayawan Bay was simulated using a 2D shallow water model to provide hydrodynamic boundary conditions for a 3D flow model, by which tsunami propagation inside the inlet of water-intake basin of the Dayawan Nuclear Power Plant (DNPP) was simulated. The 2D model covers the whole bay plus a little extension of the waters around the bay mouth. The sea side open boundary of the 2D model is generated from a superposition of an assumed tsunami process on a tide with exceedance probability of 10%. The event of encounter of a 10-min highest waters of the assumed tsunami process with such level of tide is estimated to be approximately a return period of 1490 years. Tsunami run-up in front of the revetment of the pump room was simulated using the 3D flow models. The rubble-mound breakwater was treated as porous medium in the simulation. Computational results of flow velocity indicate that tsunami propagation inside the basin is effectively prevented by the east breakwater. Computational results of waters in front of the revetment of the pump room indicate that the possible maximum water level is lower than the present elevation of the revetment, therefore the pump room has low risk of inundation.

The effect of giant tsunamis such as the Indian Ocean Tsunami in 2004 and the Great East Japan Earthquake Tsunami in 2011 has been devastating. In this study, a numerical simulation of the tsunami has been developed to estimate the physical characteristics of tsunamis and their effect on human society. Several laws and equations have been introduced for the simulation of tsunami propagation in the ocean, tsunami refraction, and tsunami run-up on land under a stable computational condition with acceptable accuracy. Our proposed method has been accepted as the world standard since 1997 and has been widely distributed to many countries through UNESCO.1) Computer graphic animations prepared by using the results of numerical simulation have been effectively used in public education and to increase the understanding of behaviors of the tsunami on the earth. When the numerical prediction of tsunami becomes possible with sufficient accuracy, then their results can be used to predict future damages and prevent the occurrence of a disaster. Data in the past were collected and expressed in terms of a newly introduced tsunami intensity which is related to the locally observed tsunami heights.

In this research work, tsunami wave propagation models and the underlying theory are examined through solving numerical wave-based partial differential equations using the leapfrog scheme. Following the wave physics theory, shallow water equations are explored for their basics and influential parameters; then, numerical methods are employed to solve said equations. As a result, the derivatives of shallow water equations for tsunami propagation modeling have been successfully formulated.

This paper aims at best describing the submarine landslide which induced partial submersion of the atolls of Mururoa and Fangataufa in 1979. More precisely, waves propagated along the south coast of Mururoa atoll and penetrated into its lagoon some minutes after the landslide triggering (t = 0 s), whereas a train of eight water waves reached the runway located on the north-east coast of Fangataufa (40 km south of Mururoa) between t = 7 min 30 s and t = 20 min. A numerical model based on shallow water equations is used to simulate the landslide as well as the associated tsunami. Saint-Venant equations are used to propagate the tsunami in coastal areas, whereas the offshore propagation is simulated by solving weakly nonlinear Boussinesq equations. Low- and high-resolution nested grids are used to simulate the tsunami propagation in deep sea and in shallow waters, respectively. Several scenarios have been tested to reproduce the observed water and run-up heights in the near and far fields. The best scenarios correspond to a landslide with a volume in the range (75–90 Mm3) (for a basal friction angle of 35°) and with a basal friction angle in the range (30°–40°) (for a volume of 80 Mm3). These results have been completed by a parametric study on the slide parameters.

Recent studies suggest that the tsunami speed can be slowed down by around 1% due to Earth elasticity, water compressibility and density stratification. Analytical solutions of wave dispersion relationship, accounting for such effects, were found in previous studies. In this paper, we investigate the additional effects of water viscosity, ocean stratification due to temperature/salinity and numerical dispersion. Theoretical solutions are derived and checked with known solutions. All the formulas are then simplified for long tsunamis waves so that the propagation speed can be calculated explicitly. The simplified solutions are evaluated using realistic geophysical parameters. For a typical tsunami wavelength of ∼200km, the viscous effect is found to be negligible; ocean stratification due to temperature/salinity causes significant speed reduction because of the high density change rate, which has been ignored before. We also evaluate the numerical dispersion of tsunami simulations, which is shown to be potentially comparable to physical dispersion.

We apply a genetic algorithm (GA) to find the optimized unit sources using dispersive tsunami synthetics to estimate the tsunami source of the 2012 Haida Gwaii earthquake. The optimal number and distribution of unit sources gives the sea surface elevation similar to that from our previous slip distribution on a fault using tsunami data, but different from that using seismic data. The difference is possibly due to submarine mass failure in the source region. Dispersion effects during tsunami propagation reduce the maximum amplitudes by up to 20% of conventional linear long wave propagation model. Dispersion effects also increase tsunami travel time by approximately 1 min per 1,300 km on average. The dispersion effects on amplitudes depend on the azimuth from the tsunami source reflecting the directivity of tsunami source, while the effects on travel times depend only on the distance from the source.

Early tsunami warning critically hinges on rapid determination of the tsunami hazard potential in real-time, before waves inundate critical coastlines. Tsunami energy can quickly characterize the destructive potential of generated waves. Traditional seismic analysis is inadequate to accurately predict a tsunami’s energy. Recently, two independent approaches have been proposed to determine tsunami source energy: one inverted from the Deep-ocean Assessment and Reporting of Tsunamis (DART) data during the tsunami propagation, and the other derived from the land-based coastal global positioning system (GPS) during tsunami generation. Here, we focus on assessing these two approaches with data from the March 11, 2011 Japanese tsunami. While the GPS approach takes into consideration the dynamic earthquake process, the DART inversion approach provides the actual tsunami energy estimation of the propagating tsunami waves; both approaches lead to consistent energy scales for previously studied tsunamis. Encouraged by these promising results, we examined a real-time approach to determine tsunami source energy by combining these two methods: first, determine the tsunami source from the globally expanding GPS network immediately after an earthquake for near-field early warnings; and then to refine the tsunami energy estimate from nearby DART measurements for improving forecast accuracy and early cancelations. The combination of these two real-time networks may offer an appealing opportunity for: early determination of the tsunami threat for the purpose of saving more lives, and early cancelation of tsunami warnings to avoid unnecessary false alarms.

Aoyama, Y.; Adityawan, M.B., Widiyanto, W., Mitobe, Y., Komori, D., and Tanaka, H., 2016. Numerical Study on Tsunami Propagation into a River. In: Vila-Concejo, A.; Bruce, E.; Kennedy, D.M., and McCarroll, R.J. (eds.), Proceedings of the 14th International Coastal Symposium (Sydney, Australia). Journal of Coastal Research, Special Issue, No. 75, pp. 1017 - 1021. Coconut Creek (Florida), ISSN 0749-0208. The Tohoku Earthquake and Tsunami in 2011 caused serious damage not only to coastal structures, but also to riverine infrastructure due to long distance of wave propagation into rivers located on the coast facing the Pacific Ocean. Although there have been numerous investigations on tsunami, limited number of studies have been made for tsunami propagation into a river channel. This study investigates tsunami propagation into a river numerically by comparing with laboratory experiment data. The present numerical simulation is based on shallow water equations, which are solved with the MacCormack sch...

Tsunamis are one of the most destructive natural phenomena (Pugh and Woodworth, 2014). For past decades, the risk of potential earthquake zones has been issued by Japanese researchers, which could occur tsunami, in Niigata and Tottori area located in the west coast of Japan and the Ryukyu trench located in Okinawa as well (Disaster Prevention Research Institute, Japan). Also, there is a Yamamoto Rise on the East Sea, where it strongly affects tsunami propagation. This topography causes high tsunami energy to concentrate on the east coast of Korea (Cho and Lee, 2013). For example, the 1983 Akita and 1993 Hokkaido earthquake induced tsunamis, Japan, respectively cause the property and life damages to the certain cities on the east coast of Korea. Therefore, it is important to propose the optimal deployment location of offshore tsunami observation instruments to contribute to the tsunami early warning system by increasing probability of tsunami detection with the minimal number of instruments considering a large number of potential tsunami scenarios. In this study, by considering various factors, the optimal location of instruments is suggested based on the numerical model results of possible tsunami scenarios.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/q3PLE6EgtrY

Tsunami modeling and simulation has changed in the past few years more than it had in decades, especially so with respect to coastal inundation. Among other things, this change is supported by the approaching era of exa-scale computing, whether via GPU or more likely forms of hybrid computing whose presence is growing across the geosciences. For reasons identified across this review, exa-scale computing efforts will impact the on-shore, highly turbulent r\'egime to higher degree than the 2D shallow water equations used to model tsunami propagation in the open ocean. This short review describes the different approaches to tsunami modeling from generation to impact and underlines the limits of each model based on the flow r\'egime. Moreover, from the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, we underline the current challenges associated with this approach and review the few efforts that are currently underway to achieve this goal. A table of existing tsunami software packages is provided along with an open Github repository to allow developers and model users to update the table with additional models as they are published and help with model discoverability.