geosciencesGeosciencesGeosciencesGeosciences2076-3263MDPI10.3390/geosciences11010005geosciences-11-00005ReviewModeling and Simulation of Tsunami Impact: A Short Review of Recent Advances and Future Challengeshttps://orcid.org/0000-0002-7498-049XMarrasSimone1*†https://orcid.org/0000-0002-8267-5989MandliKyle T.2†Department of Mechanical and Industrial Engineering Center for Applied Mathematics and Statistics, New Jersey Institute of Technology, Newark, NJ 07102, USADepartment of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027, USA; kyle.mandli@columbia.eduCorrespondence: smarras@njit.edu
Tsunami modeling and simulation has changed in the past few years more than it has in decades, especially 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 in this review, exa-scale computing efforts will impact the on-shore, highly turbulent régime to a 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égime. 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.
The failure of a multi-billion dollar wall designed to protect the To¯hoku coasts of Japan (Figure 1) from a level-2 tsunami (Level-2 tsunami: infrequent but highly destructive [1]) in 2011 triggered an important debate about alternative approaches to tsunami risk reduction. This debate is ongoing and has attracted broad media attention worldwide (Reuters [2], The Guardian [3], The Economist [4], The New York Times [5], Wired [6]). The question whether a wall is the best solution to tsunami mitigation lies in the significant expense required for a wall that does not guarantee protection from big tsunamis. Even when cost is not the main constraint—consider Japan whose GDP accounts for 4.22% of the world economy versus Indonesia’s (0.93%) or Chile’s (0.24%)—relying on traditional concrete-based solutions alone may not be desirable or sustainable, partly because of their potential long-term negative impact on the population [2], coastal ecosystems [7,8,9], and shoreline stability [10,11]. For these reasons, decision-makers and engineers are increasingly considering protection solutions that rely on green designs as sustainable and effective alternatives to seawalls. These designs, three examples of which are shown in Figure 2 along the Ring of Fire, are usually human-made hillscapes erected on the shoreline to protect the communities behind them by partially dissipating and partially reflecting the tsunami energy [12]. The question remains whether these alternatives truly protect the people and properties behind it. Robust computational simulation capabilities are required that can model such diverse types of infrastructure including the forces upon them.
While the numerical capabilities to model tsunamis in the deep ocean are well established and mature, as it will be further pointed out throughout this review, this is still not the case for the modeling of tsunami–shore interaction once the wave propagates inland. The limit lies in the level of detail necessary to correctly estimate the forces involved, especially on complex structures. In only five years since the review by Behrens and Dias [14], important steps towards high-fidelity tsunami modeling at all scales have been made, especially with respect to the fine-scale simulation of tsunami run-up and tsunami–shore interaction (see, e.g., in [15,16,17,18]), contributing to further advancing the forecast of tsunami impact as envisioned by Bernard et al. [19] ten years earlier.
In light of the development of a comprehensive, multi-scale tsunami modeling infrastructure as proposed in the latest 10-year science plan of the Natural Hazards Engineering Research Infrastructure [20], this article reviews models of tsunami propagation from generation to impact with respect to their régimes of validity and model-to-model interaction. As an aide to this, a basic derivation of the assumptions that lead to the basic two-dimensional shallow water equations are provided along with some of the most direct extensions and their limitations. Connecting these assumptions and their limitations to the much more complex three-dimensional and non-hydrostatic models follows and demonstrates the complexity the current state-of-the-art faces.
The article is organized as follows. Section 2 reviews the state of the field in tsunami modeling and simulation. The different mathematical models of classical use are described in Section 3, followed by a review of the different numerical methods to solve them given in Section 4. The idea of a future comprehensive and multi-scale simulation framework is discussed in Section 5. Conclusions are given in Section 6. A table of some available tsunami software packages is given in the Appendix A.
2. State of the Field in Tsunami Forward Modeling
The multi-scale nature of tsunami dynamics makes their study in a laboratory setting challenging [21,22]. For this reason, and supported by ever cheaper computing power and data storage, numerical modeling and simulation has become the most widely utilized tool for large-scale tsunami modeling and analysis [14]. The numerical simulation of tsunamis started in Japan sixty years ago with the work by Aida [23,24] and has shown to be effective to model their generation (see, e.g., in [25,26,27,28,29]), propagation (see, e.g., in [30,31,32]), and inundation [12,16,33,34,35], although the problem of inundation is still partially open when it comes to its numerical treatment (see, e.g., in [36,37,38] and references therein).
When it comes to tsunami–shore interaction, the effects of isolated components such as bathymetry and vegetation have been studied for idealized one-dimensional (1D) and two-dimensional (2D) settings using numerical simulations for years. For example, in [12,33,39,40,41] 1D and 2D shallow water models were used to demonstrate that the nonlinearity of tsunami waves has a significant effect on the on-shore flow velocities and propagation. Many inundation modeling efforts aim to obtain an accurate prediction of the water elevation level. While these models typically do a good job at this, their correct prediction of water level seldom translates into the accurate prediction of the forces at play [42,43], possibly indicating the limitations of 1D and 2D models in specific flow régimes.
A tsunami is a naturally multi-scale phenomenon whose characteristic length scales range from the planetary scale spanning oceans to the small scales of turbulence. On the one hand, in the open ocean it is effectively a fast moving two-dimensional long wave that travels undisturbed for thousands of kilometers and can be accurately described by the 2D nonlinear shallow water equations. On the other hand, as the flow approaches the shore and moves inland, its three-dimensionality and turbulent nature become important, with boundary layer dynamics that contribute to an important shear-driven erosion and sediment transport (see Figure 3).
As the tsunami advances further inland, sand, dirt, vegetation, structures, and other debris are scoured off of the sea bed, transported, and deposited along its path, potentially leading to important coastal morphological changes (see, e.g., in [44,45,46,47,48,49,50]). By construction, however, the shallow water equations are not apt to model such complex flows, although they have been used or enhanced in the context of erosion and sediment transport in [47,51,52,53,54], for example.
Towards an understanding of the possible limitations of the shallow water equations for a fine-grained tsunami modeling, Qin et al. [15] compared the 3D Reynolds-Averaged Navier–Stokes (RANS) solution of a turbulent propagating bore using OpenFOAM [55] against the 2D shallow water solution of the same problem using GeoClaw [56]; Qin et al. demonstrated that a 3D model for turbulent flows is necessary to correctly predict tsunami inundation and the fluid forces involved. Among the first numerical studies that attempted to model tsunami run-up as the solution of the Navier–Stokes equations, we find the 2D simulations by Hsiao and Lin [57] from 2011 and by Larsen and Fuhrman [17,18] who recently underlined the growing need for high-fidelity and multi-scale models to study tsunami risk and inland propagation. Furthermore, the RANS equations classically used in engineering were shown to overproduce turbulence beneath surface waves [58], therefore suggesting the need for high-fidelity turbulence models such as large eddy simulation (LES). LES was used in the context of plunging breakers by Christensen [59] who underlined the numerical challenges of models in reproducing experimental results on the surf zone dynamics. Much attention has been given to run-up; however, not enough work has addressed the analysis of the interaction of the tsunami with the built environment and vegetation and of the erosion caused by large tsunamis [18], all critical at improving the prediction of inundation. Modeling erosion with shallow water models can be difficult as they generally require some sort of reconstruction of vertical velocity profile to reconstruct the velocity shear in these models; however, shear is the driving mechanism of erosion and entrainment. Moreover, the modeling of the dynamics of the near-bed region and of erosion is still an open challenge in coastal engineering because it is governed by important particle–particle and particle–fluid interactions [60], which require grain-resolving simulations [61,62]. In a recent inter-comparison of models to study shear and separation driven coastal currents, Lynett et al. [63] concluded that “[...]In general, we find that models of increasing physical complexity provide better accuracy, and that low-order three-dimensional models are superior to high-order two-dimensional models[...]”.
Proper fine-grained modeling is important for studying erosion. In the context of tsunami impact mitigation analysis, the next frontier may lie in the detailed simulation of vegetation–tsunami and vegetation–morphology interactions. In light of the pioneering findings by the U.S. Army Corps of Engineers that vegetation may act as a bio-shield against flooding from storm surges [64,65], the designs of tsunami mitigation solutions around the world often incorporate vegetation. Significant experimental and numerical work has been done to analyze the effect of vegetation on wind waves, currents and, albeit less so, tsunamis (see, e.g., in [66,67,68,69,70,71,72,73,74,75,76]), particularly in the context of steady flows [77]. After the 2004 tsunami in Sumatra, Bayas et al. [78] estimated that vegetation along the west coast of Aceh may have reduced casualties by 5%. The benefits of coastal vegetation are not limited to attenuation. It is well known that bathymetry, topography, and coastal geomorphology (rivers, canals, and barrier islands) have a profound effect on run-up for tsunamis [12] and storm surge [79,80,81,82]. Furthermore, the presence of vegetation alters sediment-transport processes and landform evolution [83,84,85,86,87]. The modeling of the interaction between the tsunami and vegetation is numerically difficult because of the flexible nature of vegetation. Efforts have been made to include the effect of flexible vegetation on the flow. Examples of one- and two-way coupling of the flow with the dynamics of vegetation can be found in the recent work by Mattis et al. [88,89] and Mukherjee et al. [90], who, for the first time in the context of inundation, used fluid–structure interaction algorithms to study the effect of deformation of idealized trees on the flow (See Figure 4).
3. Mathematical Representations and Assumptions
While the most complete model to describe a tsunami as a highly turbulent free surface flow is given by the Navier–Stokes equations of incompressible flows, simpler models of wave propagation have been successfully adopted for decades to study tsunami propagation in the open ocean and run-up [14,91]. In this section, we introduce these models and underline their properties and limits of validity with the aim to identify the assumptions that go into them and where those assumptions may break down, especially in the near-shore. Based on this analysis, we justify the ever increasing use of the 3D Navier–Stokes equations to study the tsunami boundary layer in the near-shore and during run-up where 1D and 2D models were shown to lead to inaccurate solutions [15,63], leading to the idea suggested in Figure 5 of having both the 2D shallow water and 3D Navier–Stokes models working together to provide a full picture of a tsunami’s impact.
3.1. The 3D Navier–Stokes Equations
We first start with arguably the best model of water flow, the three-dimensional Navier–Stokes equations. During inundation, the interaction of the flow with the coastal features, erodible terrain, sediments, vegetation, and structures is such that the flow is fully turbulent, and thus three-dimensional and characterized by shear. The correct estimation of the hydrodynamic forces responsible for the damages during run-up requires a model that can capture these flow characteristics [17,63]. The Navier–Stokes equations are the most complete model that can capture these features. Omitting their derivation from first principles, the 3D Navier–Stokes equations of incompressible flows are written as
ux+vy+wz=0ut+(u2)x+(uv)y+(uw)z=−1ρPx+ν(uxx+uyy+uzz)vt+(vu)x+(v2)y+(vw)z=−1ρPy+ν(vxx+vyy+vzz)wt+(wu)x+(wv)y+(w2)z=−1ρPz+ν(wxx+wyy+wzz)−g
3.2. Depth-Averaged Models
Although (1) are the model equations for tsunamis that we would like to use, they are computationally too costly to use everywhere and more than likely not needed everywhere in the domain of interest for a tsunami. We therefore need to carefully identify where and when we can approximate (1). The primary type of approximation associated with tsunamis are characterized by averaging through the depth of the water column in the gravity centered coordinate system (see Figure 6). This precludes modifications that would handle steeper terrain in favor of reducing complexity but should be noted as an alternative for simple topographical features. We will also forgo analysis of layered depth-averaged models as though they have unique properties and largely follow the same derivation.
3.2.1. Scaled Equations
We first start with the inviscid version of the two-dimensional Navier–Stokes equations (i.e., the incompressible Euler equations) of incompressible flows, one vertical and horizontal dimension, where we will ignore one of the horizontal directions for simplicity of notation and without loss of generality. With this setup and assuming a free surface and non-flat bottom boundary we have the equations
ux+wz=0ut+(u2)x+(uw)z=−1ρPxwt+(uw)x+(w2)z=−1ρPz−g,
with the boundary conditions
w=ηt+uηxP=Paatz=ηw=ubxatz=b,
where the velocity is represented by u in the horizontal direction and w in the vertical, ρ is the density of water, P its pressure—which includes non-hydrostatic components at this juncture, and g is the gravitational acceleration. The underscores t,x, and y indicate the partial derivatives with respect to them.
The first step in defining and justifying the depth-averaged equations lies in the non-dimensionalization assumptions made during their definition. Because of this, we will be explicit about these assumptions as they will play an important role later on. The primary scalings are
xλ=x˜za=z˜tT=t˜
where λ is a length scale, often taken as the characteristic wave-length of the waves involved; a the characteristic depth that the wave is propagating through; and T the characteristic time-scale or period of the wave. With these values defined, we can also normalize the velocities and P such that
uU=u˜wW=w˜PP0=P˜
where P0 is the pressure normalization, usually taken to be atmospheric pressure, and the non-dimensional quantities are those with ·˜. Introducing then the traditional shallow water parameter ϵ=a/λ we can write the velocity and temporal normalizations as U=ga, W=ϵU, and T=λ/U.
Applying these scalings to (2) and simplifying we are lead to the system of equations
u˜x˜+w˜z˜=0u˜t˜+(u˜2)x˜+(u˜w˜)z˜=−P˜x˜ϵ2w˜t˜+(u˜w˜)x˜+(w˜2)z˜=−P˜z˜−g
where we have also assumed that P0/ρU2=1 without loss of generality. For the boundary conditions we introduce a new scaling, η/δ=η˜ and b/β=b˜ which, for the top boundary, results in
ϵλδw˜=η˜t˜+u˜η˜x˜ϵλβw˜=u˜b˜x˜.
where δ is surface amplitude scaling. For the upper boundary condition, if δ/λ=O(ϵ) the boundary condition is all of the same order. This is probably not the case as δ≪a unless in very shallow water but also should be noted for later.
The lower boundary condition is a similar situation, if β/λ=O(ϵ) the boundary condition is all of the same order. In the case of the bathymetry scale, this is probably the case as the bathymetry will vary on the scale of the depth scale a. The final non-dimensionalized equations are then (dropping the ·˜)
ux+wz=0ut+(u2)x+(uw)z=−Pxϵ2wt+(uw)x+(w2)z=−Pz−1
with the boundary conditions
ϵλδw=ηt+uηxP=Pa′z=ηϵλβw=ubxz=b
where Pa′ is the appropriately scaled atmospheric pressure condition.
3.2.2. Depth Integration
The next step to deriving the shallow water and similar equations is to depth integrate (7)–(9). Here, we will skip many of the details other than to point out two important assumptions that are often misrepresented. Continuing on the first step is to integrate the incompressibility Equation (7) such that we find
∫bηux+wzdz=0⇒ht+∂∂x∫bηudz=0.
The next equation to be integrated is the horizontal momentum Equation (8). For this we will first introduce the assumption that the total pressure P can be split up additively into an hydrostatic and non-hydrostatic component such that P(x,z,t)=Pa+(η−z)+p where the term (η−z) represents the depth from the surface η and therefore the non-dimensionalized hydrostatic component of the pressure and p(x,z,t) non-hydrostatic pressure component. Then, integrating the left-hand-side of (8), we find
∫bηut+(u2)x+(uw)zdz=∂∂t∫bηudz+∂∂x∫bηu2dz+∫bηpdz+12h2
where we have simplified the notation using our average values making sure to maintain the non-commutativity of the average and squaring operators. For the right-hand-side of (8) we similarly conclude that
−∫bη(Pa+(η−z)+p)xdz=−bx(h+p|z=b).
Finally, we also integrate the vertical momentum Equation (9) through the depth. This is very similar to the horizontal equations leading to
ϵ2∫bηwt+(uw)x+(w2)zdz=ϵ2∂∂t∫bηwdz+∂∂x∫bηuwdz
for the left-hand-side and
−∫bηPa+(η−z)+p)z+1dz=p|z=b
for the right-hand side of the equation.
Putting all this together leads us to the vertically integrated equations
ht+∂∂x∫bηudz=0∂∂t∫bηudz+∂∂x∫bηu2dz+∫bηpdz+12h2=−bx(h+p|z=b)ϵ2∂∂t∫bηwdz+∂∂x∫bηuwdz=p|z=b.
At this juncture it is useful to start to introduce the notation that indicates average quantities through the depths. We will denote these by
·¯=1h∫bη(·)dz.
Herein lies one of the critical misrepresentations of depth-averaged models, that they assume constant velocity, or even more general quantities, throughout their vertical profiles. For example, in general, u2¯≠u¯2. In fact we will make some assumptions about the commutativity of the averaging operator with others but this does not preclude the consideration of non-constant functions of the velocity for instance. This leads to the following system of equations,ht+(hu¯)x=0(hu¯)t+hu2¯+hp¯+12h2x=−(h+p|z=b)bxϵ2(hw¯)t+(uw¯)x=p|z=b
One last equation for the pressure will be useful as it will give us a means to calculate p¯ depending on the approximations that we will make in the next section. We can come to this equation by reconsidering the scaled vertical momentum equation and integrating from η to a vertical level b+αh where α∈[0,1] leading toϵ2∂∂t∫b+αhηwdz+∂∂x∫b+αhηuwdz+w[α(ηt+uηx)+(1+α)ubx−w]z=b+αh=p|z=b+αh.
3.2.3. Approximations
At this point we can use the equations from the previous section to derive a number of approximations commonly used in tsunami numerical modeling. This will not be exhaustive but rather suggestive to where this basic framework can be taken to derive and analyze the validity of these approximations.
The first of these approximations will assume that the averaging operator will commute with multiplication, e.g., u¯2≈u2¯. The terms that this ignores are often termed as differential advection; however, note that this approximation does not imply that the velocity profiles are constant but rather that the averages commute. This can be important when comparing boundary layer physics for instance. Moving forward, this allows us to rewrite (11) and simplify notation so that
ht+(hu)x=0(hu)t+hu2+hp+12h2x=−(h+p|z=b)bxϵ2(hw)t+(uw)x=p|z=b
where we have also now dropped the ·¯ notation.
3.2.4. Shallow Water
Finally for the shallow water equations we simply need to assume that the shallow water parameter ϵ≪1 implying that 0=p|z=b. As the non-hydrostatic pressure is defined such that p(x,η,t)=0 we can show that p¯=0 therefore implying that (11) become
ht+(hu)x=0(hu)t+hu2+12h2x=−hbx,
which are of course the traditional, non-dimensionalized shallow water equations. Note that (12) can be used here to show that p¯=0 is in fact true here.
3.2.5. Not So Shallow Equations
The next perhaps logical step is to maintain the assumption that differential advection is still small and therefore averages commute but that vertical momentum is not as ignorable as before, i.e., ϵ≥O(1). Instead we might make some assumptions about the depth profile of w and use that to derive a form of a closure expression for the non-hydrostatic pressure. Here, we will give one example of an analysis where we assume that w has a linear profile through the depth along with the equations that arise due to this assumption.
First, if w is linear in z we can use the incompressibility condition to find w(x,z,t)=w(x,t)−uxz. The crux of these calculations is then the computation of w¯ and p¯. Leaving out the tedious details of these calculations but reporting the results the final equations with these assumptions leads to
ht+(hu)x=0,(hu)t+hu2+12h2+hp¯x=−bx(h+p|z=b),and(hw)t+(huw)x−hux12h+bt+huux12h+bx=p|z=b,
where
p¯=ϵ2h∂∂t12hw+12uxhb+13uxh2+∂∂x12huw+14(u2)xhb+16(u2)xh2−32w2+32uwbx+3bwux−32b2(ux)2−34((ub)2)x+116hwux−32hb(ux)2−14h(bu2)x−23(hux)2−13hbx(ux)2
Note that this complexity is a result of the attempt to recreate the 3D structure that we have already integrated out. This may be counterintuitive as to why we might want to do this but in fact schemes of this nature have been quite successful, namely, the Boussinesq, Serre–Green–Naghdi, and other types of semi-2D types of equations.
3.3. Mathematical Conclusions
As we have seen, (1) can be greatly reduced in complexity by integration and appropriate assumptions. While this means that the 2D equations can easily recreate the extent of flooding, the recreation of the 3D velocity fields and therefore the appropriate forces on structures is however difficult even for schemes that have been successfully implemented such as the Boussinesq and Serre–Green–Naghdi models, leaving the need for 3D Navier–Stokes models as the gold standard for impact questions. Furthermore, as we have shown, the assumptions we have made are often not as stringent as is often thought in the literature. The importance of this is that the shallow water 2D field does not preclude there being a complex horizontal velocity field, rather that its average have the special property that it commutes as mentioned above. This critical property prescribes exactly the compatibility condition that we need to have between a 2D shallow water model and a 3D Navier–Stokes model whose field does not require this commutativity property.
4. Numerical Solution and Computational Considerations
The equations described across this paper have been solved for decades using a wide range of numerical approximations. From the classical finite difference method (FD) used in, e.g., MOST [92], to finite volumes (FV; see, e.g., [31,56,93,94,95]), to Galerkin methods such as finite elements (FE; e.g., [90,96,97]) and high order spectral elements and discontinuous Galerkin (SE and DG; e.g., [12,32,36,97]), to the less common smoothed particle hydrodynamics method (SPH) used in, e.g., [98,99].
Each numerical method comes with its advantages and disadvantages. FD are the least costly, but they lack the geometrical flexibility of Galerkin methods and FV. FE/DG and FV are optimal for grid refinement [32,56,100] and are geometrically flexible to handle complex coastal lines; furthermore, FE/SE/DG are inherently optimal for massive parallelism [101]. SPH is good at tracking the free surface of the three dimensional flow, but it is expensive and it is still unclear how to treat boundary layers with it. Regardless of the numerical method of approximation, the numerical handling of the wet–dry interface to model inundation has been improved throughout the years as demonstrated in, e.g., [36,37,38,102,103,104,105,106]. It should also be noted that new methods in high-performance computing continue to push the envelope forward and will consequently lead to advances in tsunami science. This advancement will be either hampered by the same issues encountered in general geophysical turbulence computational modeling, or enjoy the benefits of embarrassingly parallel ensemble modeling as is the case with probabilistic tsunami hazard assessment (PTHA).
While Equation (1) models turbulent incompressible flows, their solution on a computational grid using any of the methods above (except for SPH) requires some attention with respect to the treatment of the viscous stresses. If we could afford to have an infinitely fine computational grid, they could be solved directly; this approach is known as direct numerical simulation (DNS). However, even with the approaching era of exa-scale computing, DNS is much too costly to be viable for tsunami simulations (The number of grid points in a computational grid necessary for DNS scales as the (9/4)th power of the Reynolds number.) [107]. Instead, by means of either a Reynolds averaging approach that gives rise to the RANS equations, or scale separation filtering for LES, turbulence can be modeled with sufficient precision at a drastically reduced computational cost with respect to DNS. RANS was first used to simulate tsunami run-up in the software COBRAS [108,109]. However, it took approximately one more decade to become a popular choice to model tsunami generated turbulence. Examples of RANS-based tsunami models are described in [17,18,110,111,112].
In terms of solution cost, LES lies between DNS and RANS. Due to its computational cost, it is not popular among tsunami modelers, although it was recently considered as an optimal choice to study breaking waves and run-up in highly non-steady dynamics (see, e.g., in [90]). As the era of exa-scale computing on hybrid CPU–GPU architectures is fast approaching, we expect LES to become the model of choice among tsunami inundation modelers interested in the analysis of tsunami impact and coastal propagation, although it would be unreasonable to use either RANS or LES to model for off-shore propagation where, instead, the 2D shallow water equations do a perfect job.
The free surface in a model based on the solution of (1) is typically accounted for in either one of two ways: by means of a volume of fluid approach (in, e.g., [112,113]) or by a level-set tracking of the water–air interface in a two-fluid model (in, e.g., [90,114]). Both of these approaches make it natural to handle the wet–dry front as the flow reaches land; when high order methods are utilized, this may not be the case for shallow-water equations whose solution is challenged in wet–dry regions for certain numerical approximations, as demonstrated by the continuous work of, for example, [36,37,38,103,104,105].
In summary, all the numerical methods of solution as well as the mathematical models have their use case and there is no one that can rule them all. Similarly, the adage “using the right tool for the job” holds here as much as it does when doing carpentry. It is important therefore to consider what metrics may lead to the decisions to choose one model over another. The most prominent metrics are listed as follows.
Flow scale and régime. For example, are we modeling breaking waves in the vicinity of the shore or linear waves propagating in the ocean basin? Is turbulence an important factor?
Complexity of the physics needed. Here the difference between using a wall boundary condition at the shore and doing true wetting and drying may be significant as does the representation of true turbulent flow.
Performance on the computing architecture being considered.
Overall size of the problem considered. Is one interested in a single simulation or a large ensemble in order to account for uncertainty?
These are of course only some of the issues but they highlight the difficulties faced and the need for hybrid solver techniques to bridge the gaps between the advantages and disadvantages between solvers.
5. Towards a Multi-Scale Framework from Source to Impact
It seems plausible, given the ongoing research on all fronts from tsunami generation, propagation, to fine-scale run-up, that tsunami modelers may soon join forces to build a unified, all-scale framework to model a full tsunami event from source to impact (see, e.g., in [27,115]). The all-scale idea is naïvely represented in Figure 5 where the tsunami’s domain is subdivided in different zones as a function of the flow régime and complexity. The tsunami assumes different shapes during its lifespan from generation to impact; with time scales that range from seconds (100–1000 s at the source [50]) to up to several hours during propagation and run-up/run-down, the flow is driven by both two-dimensional and three-dimensional dynamics across spatial scales from submeters to hundreds of kilometers. Due to such enormous scale variability in time and space, not all régimes can or should be described by the same physical model. While it may be or become technically possible, it is not sensible to solve the 3D Navier–Stokes equations for turbulent flows in the far-field where the shallow water equations are appropriate and inexpensive to solve. The wave propagates in the open ocean with a small amplitude from a few centimeters up to several tens of centimeters (Ref. [50] and citations therein), which is always much smaller than the open ocean depth which is, in turn, much smaller than the wavelength. Under these circumstances, the tsunami motion is well represented by the shallow water equations [116,117,118] for linear, quasi-linear, and quasi-dispersive waves. While the wave approaches land it decelerates, its amplitude increases, and its wavelength is drastically reduced; furthermore, it loses its linear, non-dispersive nature. In these conditions a non-linear model that describes turbulent breaking waves becomes necessary. The most complete model of this is given by the 3D Navier–Stokes equations of incompressible flows. If the only metric of interest during run-up were the inundation level, the 2D shallow water equations would still provide relatively satisfactory results, but their limits become clear when the correct hydrodynamic forces on structures are to be evaluated.
To achieve a complete multi-scale model of the tsunami, different solvers may be coupled as illustrated in such a way that the 3D Navier–Stokes solver is forced at the boundaries of the shore proximity area by the shallow water model of the tsunami moving off-shore [119,120]. At the source the shallow water model is forced by an earthquake simulator or landslide model. Arguably the most difficult and still unclear part of a modeling infrastructure of this type lies in the coupling across models. Coupling across models is an open field of research of its own (see, e.g., in [121,122,123,124].) Parallel performance, data exchange, and time scale interactions are among the difficulties to be overcome in designing coupling algorithms across software packages. We envision a major effort to optimize software coupling and make it efficient for costly tsunami simulations towards real-time modeling of a full tsunami event.
6. Conclusions
This article gives a concise review of the different approaches to tsunami modeling from generation to impact. From the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, the article underlines the challenges associated with this approach and reviews the efforts that are underway to achieve this goal.
The current state of the field seems to indicate that the numerical tsunami modeling community is moving towards the high-fidelity modeling of tsunamis across scales, either by means of grid refinement or high order methods in the open ocean, or by means of a complete three-dimensional model of the flow in the near shore region. We expect an increasing effort to couple different models and allow the interaction of an earthquake simulator with the large scale dynamics in the open sea and fine scale dynamics in the coastal region, all within the same software infrastructure. While improved earthquake–tsunami modeling at the source, adaptive grid refinement, and ever faster and inexpensive hybrid high-performance computing (i.e., graphics processing units) have contributed to the advancement of tsunami modeling at large scales, we are now witnessing the very beginning of a new tsunami modeling era; the introduction of high-fidelity simulations of inundation and enhanced coupling of software packages across all scales is finally leading the way towards full tsunami forecast as envisioned by Synolakis et al. [91] fifteen years ago.
While it is difficult to estimate how long it will be before a full scale tsunami simulation can be run in real time (or faster) on a laptop, the pace at which the necessary tools are being developed is fast. At risk of being too optimistic, the advent of exa-scale computing within the next decade and the ever-increasing presence of general purpose GPUs mounted on personal computers by default led us to think that such a tool may be available within the next two decades.
Author Contributions
S.M. and K.T.M. contributed equally to the planning and writing of this review. All authors have read and agreed to the published version of the manuscript.
Funding
This material is based upon work supported by the National Science Foundation under Grant No. DMS-1720288 and OAC-1735609 and work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under Cooperative Agreement No. 1852977.
Acknowledgments
The authors would like to thank two anonymous reviewers whose comments and suggestions helped improve the final version of the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
Appendix A
Some of the most common tsunami simulators are listed in Table A1. The table is maintained in an open repository hosted at https://github.com/mandli/tsunami-models to allow developers and scientists to add their models.
geosciences-11-00005-t0A1_Table A1
Some available tsunami models. Acronyms used in this table: SW Shallow Water, B Boussinesq, SGN Serre–Green–Naghdi, N−S Navier–Stokes, SPH Smoothed Particle Hydrodynamics, LES Large Eddy Simulation, WM Wall Modeled LES, RANS Reynolds Averaged Navier–Stokes equations, FSI Fluid–Structure Interaction, MP Multi-Phase capabilities for erosion and sediment transport, and SD for spatial discretization. The use of 2D(1/2) indicates that the model can run in multi-layer modality [95]. Storm surge and ocean dynamics models are omitted.
Model
Space Dim.
Equations
Turbul.
Wave Break.
FSI
MP
SD
GeoCLAW [56]
1D/2D/2D12
SW
No
No
No
No
FV
NUMA2D [125,126]
1D/2D
SW
No
No
No
No
SE/DG
MOST [92]
1D/2D
SW
No
No
No
No
FD
Cliffs [92,127]
1D/2D
SW
No
No
No
No
FD
Tsunami-HySEA [128,129,130]
1D/2D
SW/B
No
Yes
No
No
FV
Multilayer-HySEA [131,132]
1D/2D(1/2)
SW/B
No
Yes
No
Yes
FV
TUNAMI [133,134]
1D/2D
SW
No
No
No
No
FD
NAMI-DANCE [135]
1D/2D
SW
No
No
No
No
FD
COMCOT [136]
1D/2D
SW
No
No
No
No
FD
SELFE [96]
1D/2D
SW
No
No
No
No
FE
TsunAWI [97]
1D/2D
SW
No
No
No
No
FE
TsunaFlash [137]
1D/2D
SW
No
No
No
No
FE/DG
VOLNA [94,138]
1D/2D
SW
No
No
No
No
FV
Delft3D [139]
1D/2D
SW
No
No
No
Yes
FD
Basilisk [140,141,142]
2D/3D
SGN
No
Yes
No
Yes
FV
BOSZ [143]
1D/2D
B
No
No
No
No
FV/FD
Celeris [144]
1D/2D
B
No
No
No
No
FV
FUNWAVE [145,146]
1D/2D
B
No
No
No
No
FV/FD
pCOULWAVE [147,148]
2D/3D
B
Yes
No
No
No
FV
NEOWAVE [149]
2D
B
No
No
No
No
FD
GPUSPH [98]
3D
SPH
No
Yes
No
No
SPH
SCHISM [111]
1D/2D/3D
N−S
RANS
Yes
No
No
FE/FV
COBRAS [108,109]
2D/3D
N−S
RANS
Yes
No
No
FD
TSUNAMI3D [110,150]
2D/3D
N−S
RANS
Yes
No
No
FD
waves2FOAM [17,18,112]
2D (tsunami)
N−S
RANS
Yes
No
No
FV
NHWAVE [151]
2D/3D
N−S
LES
Yes
Yes
Yes
FV/FD
Alya [90,152]
2D/3D
N−S
LES/WM/RANS
Yes
Yes
Yes
FE
ReferencesIshiwatariM.SagaraJ.ReutersSeven Years after Tsunami, Japanese Live Uneasily with Seawalls2018Available online: https://www.reuters.com/article/us-japan-disaster-seawalls/seven-years-after-tsunami-japanese-live-uneasily-with-seawalls-idUSKCN1GL0DK(accessed on 21 December 2020)GuardianT.After the Tsunami: Japan’s Sea Walls–in Pictures2018Available online: https://www.theguardian.com/world/gallery/2018/mar/09/after-the-tsunami-japan-sea-walls-in-pictures(accessed on 21 December 2020)EconomistT.The Great Wall of Japan2014Available online: https://www.economist.com/asia/2014/06/14/the-great-wall-of-japan(accessed on 21 December 2020)TimesN.Y.Seawalls Offered Little Protection against Tsunami’s Crushing Waves2011Available online: https://www.nytimes.com/2011/03/14/world/asia/14seawalls.html(accessed on 21 December 2020)WiredW.Ominous Views of Japan’s New Concrete Seawalls2018Available online: https://www.wired.com/story/photo-gallery-japan-seawalls/(accessed on 21 December 2020)PetersonM.LoweM.Implications of Cumulative Impacts to Estuarine and Marine Habitat Quality for Fish and Invertebrate ResourcesDuganJ.HubbardD.Ecological Effects of Coastal Armoring: A Summary of Recent Results for Exposed Sandy Beaches in Southern CaliforniaBulleriF.ChapmanM.The introduction of coastal infrastructure as a driver of change in marine environmentsDeanR.G.DalrympleR.A.KomarP.LunghinoB.Santiago TateA.MazereeuwM.MuhariA.GiraldoF.MarrasS.SuckaleJ.The protective benefits of tsunami mitigation parks and ramifications for their strategic designRenew Yamada2011Available online: http://renewyamada.org(accessed on 10 July 2020)BehrensJ.DiasF.New computational methods in tsunami scienceQinX.MotleyM.LeVequeR.GonzálezF.MuellerK.A comparison of a two-dimensional depth-averaged flow model and a three-dimensional RANS model for predicting tsunami inundation and fluid forcesPrasetyoA.YasudaT.MiyashitaT.MoriN.Physical Modeling and Numerical Analysis of Tsunami Inundation in a Coastal CityLarsenB.FuhrmanD.Full-scale CFD simulation of tsunamis. Part 1: Model validation and run-upLarsenB.FuhrmanD.Full-scale CFD simulation of tsunamis. Part 2: Boundary layers and bed shear stressesBernardE.MofjeldH.TitovV.SynolakisC.GonzálezF.Tsunami: Scientific frontiers, mitigation, forecasting and policy implicationsNHERIChenJ.HuangZ.JiangC.DengB.LongY.An experimental study of changes of beach profile and mean grain size caused by tsunami-like wavesJiangC.ChenJ.YaoY.LiuJ.DengY.Study on threshold motion of sediment and bedload transport by tsunami wavesAidaI.Numercal experiments for the tsunami propagation of the 1964 Niigata tsunami and 1968 Tokachi-Oki tsunamiAidaI.Numerical computational of a tsunami based on a fault origin model of an earthquakeBorreroJ.C.LeggM.R.SynolakisC.E.Tsunami sources in the southern California bightPeltiesC.de la PunteP.AmpueroJ.P.BrietzkeG.KäserM.Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructured tetrahedral meshesLópez-VenegasA.HorrilloJ.Pampell-ManisA.HuérfanoV.MercadoA.Advanced Tsunami Numerical Simulations and Energy Considerations by use of 3D-2D coupled Models: The October 11, 1918, Mona Passage TsunamiGálvezP.AmpueroJ.P.DalguerL.SomalaS.N.Nissen-MeyerT.Dynamic earthquake rupture modelled with an unstructured 3-D spectral element method applied to the 2011 M9 Tohoku earthquakeUlrichT.VaterS.MaddenE.BehrensJ.van DintherY.van ZelstI.FieldingJ.LiangC.GabrielA.A.Coupled, Physics-Based Modeling Reveals Earthquake Displacements are Critical to the 2018 Palu, Sulawesi TsunamiTitovV.GonzálezF.BernardE.EbleM.C.MofjeldH.NewmanJ.C.VenturatoA.Real-Time Tsunami Forecasting: Challenges and SolutionsLeVequeR.GeorgeD.BergerM.Tsunami modelling with adaptively refined finite volume methodsBonevB.HesthavenJ.S.GiraldoF.X.KoperaM.A.Discontinuous Galerkin scheme for the spherical shallow water equations with applications to tsunami modeling and predictionLynettP.J.Effect of a shallow water obstruction on long wave runup and overland flow velocityParkH.CoxD.T.LynettP.J.WiebeD.M.ShinS.Tsunami inundation modeling in constructed environments: A physical and numerical comparison of free-surface elevation, velocity, and momentum fluxOishiY.ImamuraF.SugawaraD.Near-field tsunami inundation forecast using the parallel TUNAMI-N2 model: Application to the 2011 Tohoku-Oki earthquake combined with source inversionsVaterS.BeisiegelN.BehrensJ.A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: One-dimensional caseVaterS.BeisiegelN.BehrensJ.A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: Triangular gridsAyogJ.KesserwaniG.ShawJ.SharifianM.BauD.Second-order discontinuous Galerkin flood model: Comparison with industry-standard finite volume modelsLinP.A numerical study of solitary wave interaction with rectangular obstaclesApotsosA.JaffeB.GelfenbaumG.Wave characteristic and morphologic effects on the onshore hydrodynamic response of tsunamisOkalE.SynolakisC.Sequencing of tsunami waves: Why the first wave is not always the largestLynettP.BorreroJ.WeissR.GreerD.RenteriaW.Observations and modeling of tsunami-induced currents in ports and harborsBorreroJ.LynettP.KalligerisN.Tsunami currents in portsBorreroJ.SynolakisC.FritzH.Northern Sumatra field survey after the December 2004 Great Sumatra Earthquake and Indian Ocean tsunamiParisR.WassmerP.SartohadiJ.LavigneF.BarthomeufB.DesgagesE.GrancherD.BaumertP.VautierF.BrunsteinD.Tsunamis as geomorphic crises: Lessons from the december 26, 2004 tsunami in Lhok Nga, West Banda Aceh (Sumatra, Indonesia)KatoF.NoguchiK.SuwaY.SakagamiT.SatoY.Field survey on tsunami induced topographical changeKuriyamaY.TakahashiK.YanagishimaS.TomitaT.Beach profile change at Hasaki, Japan caused by 5-m-high tsunami due to the 2011 off the Pacific coast of Tohoku earthquakeYamashitaK.SugawaraD.TakahashiT.ImamuraF.SaitoY.ImatoY.KaiT.UeharaH.KatoT.NakataK.Numerical simulations of large-scale sediment transport caused by the 2011 Tohoku earthquake tsunami in Hirota Bay. Southern Sanriku CoastUdoK.TakedaY.TanakaH.Coastal morphology change before and after 2011 off the Pacific coast of Tohoku earthquake tsunami at Rikuzen-Takata coastLevinB.NosovM.LiL.QiuQ.HuangZ.Numerical modeling of the morphological change in Lhok Nga, West Banda Aceh, during the 2004 Indian Ocean tsunami: Understanding tsunami deposits using a forward modeling methodOntowirjoB.ParisR.ManoA.Modeling of coastal erosion and sediment deposition during the 2004 Indian Ocean tsunami in Lhok Nga, Sumatra, IndonesiaSugawaraD.TakahashiT.ImamuraF.Sediment transport due to the 2011 Tohoku-Oki tsunami at Sendai: Results from numerical modelingGarres-DíazJ.Castro DíazM.KoellermeierJ.Morales de LunaT.Shallow Water Moment models for bedload transport problemsOpenFOAM v2.3.12014Available online: https://openfoam.org/release/2-3-1/(accessed on 21 December 2020)BergerM.GeorgeD.LeVequeR.MandliK.The GeoClaw software for depth-averaged flows with adaptive refinementHsiaoS.C.LinT.C.Tsunami-like solitary waves impinging and overtopping an impermeable seawall: Experiment and RANS modelingLarsenB.FuhrmanD.On the over-production of turbulence beneath surface waves in Reynold-averaged Navier-Stokes modelsChristensenE.Large eddy simulation of spilling and plunging breakersMeiburgE.RadhakrishnanS.Nasr-AzadaniM.Modeling Gravity and Turbidity Currents: Computational Approaches and ChallengesYuX.OzdemirC.E.HsuT.J.BalachandarS.Numerical Investigation of Turbulence Modulation by Sediment-Induced Stratification and Enhanced Viscosity in Oscillatory FlowsUhlmannM.An immersed boundary method with direct forcing for the simulation of particulate flowsLynettP.J.GatelyK.WilsonR.MontoyaL.ArcasD.AytoreB.BaiY.BrickerJ.D.CastroM.J.CheungK.F.Inter-model analysis of tsunami-induced coastal currentsUSACEFosbergF.R.ChapmanV.Mangroves v. tidal wavesMazdaY.MagiM.KogoM.HongP.N.Mangroves as a coastal protection from waves in the Tong King Delta, VietnamMöllerI.SpencerT.FrenchJ.R.LeggettD.J.DixonM.Wave transformation over salt marshes: A field and numerical modeling study from North Norfolk, EnglandHaradaK.ImamuraF.Effects of coastal forest on tsunami hazard mitigation—A preliminary investigationTanakaN.SasakiY.MowjoodM.JinadasaK.HomchuenS.Coastal vegetation structures and their functions in tsunami protection: Experience of the recent Indian Ocean tsunamiTanakaN.NandasenaN.JinadasaK.SasakiY.TanimotoK.MowjoodM.Developing effective vegetation bioshield for tsunami protectionStoesserT.Palau SalavadorG.RodiW.DiplasP.Large eddy simulation of flow through submerged vegetationIimuraK.TanakaN.Numerical simulation estimating effects of tree density distribution in coastal forest on tsunami mitigationTanakaN.JinadasaK.MowjoodM.FaslyM.Coastal vegetation planting projects for tsunami disaster mitigation: Effectiveness evaluation of new establishmentsTanakaN.Effectiveness and limitations of vegetation bioshield in coast for tsunami disaster mitigagationThuyN.B.TanakaN.TanimotoK.Tsunami mitigation by coastal vegetation considering the effect of tree breakingMuhammadR.TanakaN.Energy Reduction of a Tsunami Current through a Hybrid Defense System Comprising a Sea Embankment Followed by a Coastal ForestNepfH.Drag, turbulence, and diffusion in flow through emergent vegetationBayasJ.MarohnC.DerconG.DewiS.PiephoH.JoshiL.van NoordwijkM.CadischG.Influence of coastal vegetation on the 2004 tsunami wave impact in West AcehWesterinkJ.LuettichR.BaptistaA.ScheffnerN.FarrarP.Tide and storm surge predictions using finite element modelKânoğluU.SynolakisC.Long wave runup on piecewise linear topographiesWesterinkJ.LuettichR.FeyenJ.AtkinsonJ.DawsonC.RobertsH.PowellM.DunionJ.KubatkoE.PourtaheriH.A basin-to-channel-scale unstructured grid hurricane storm surge model applied to southern LouisianaSunF.CarsonR.Coastal wetlands reduce property damage during tropical cyclonesJacksonN.NordstromK.FeaginR.SmithW.Coastal geomorphology and restorationDayJ.BoeschD.ClairainE.E.A.Restoration of the Mississippi Delta: Lessons from Hurricanes Katrina and RitaGedanK.KirwanM.WolanskiE.BarbierE.SillimanB.The present and future role of coastal wetland vegetation in protecting shorelines: Answering recent challenges to the paradigmShephardC.CrainC.BeckM.The protective role of coastal marshes: A systematic review and meta-analysisHorstmanE.Dohmen-JanssenC.NarraP.van den BergN.SiemerinkM.HulscherS.Wave attenuation in mangroves: A quantitative approach to field observationsMattisS.DawsonC.KeesC.FarthingM.An immersed structure approach for fluid-vegetation interactionMattisS.KeesC.WeiM.DimakopoulosA.DawsonC.N.Computational Model for Wave Attenuation by Flexible VegetationMukherjeeA.CajasJ.HouzeauxG.LehmkuhlO.VázquezM.SuckaleJ.MarrasS.Using fluid-structure interaction to evaluate the energy dissipation of a tsunami run-up through idealized flexible treesSynolakisC.BernardE.Tsunami science before and beyond boxing dayTitovV.V.GonzálezF.GeorgeD.L.LeVequeR.L.Finite volume methods and adaptive refinement for global tsunami propagation and local inundationDutykhD.PoncetR.DiasF.The VOLNA code for the numerical modeling of tsunami waves: Generation, propagation and inundationMandliK.A numerical method for the two layer shallow water equations with dry statesZhangY.BaptistaA.An efficient and robust tsunami model on unstructured grids. Part I: Inundation benchmarksHarigS.ChaeroniX.PranowoW.BehrensJ.Tsunami simulations on several scales: Comparison of approaches with unstructured meshes and nested gridsWeiZ.DalrymplR.HéraultA.BilottaG.RusticoE.YehH.SPH modeling of dynamic impact of tsunami bore on bridge pierKlappJ.Areu-RangelO.CruchagaM.AranguizR.BonasiaR.GodoyM.Silva-CasarínR.Tsunami hydrodynamic force on a building using a SPH real-scale numerical simulationBehrensJ.MüllerA.KoperaM.MarrasS.WilcoxL.IsaacT.GiraldoF.Strong scaling for numerical weather prediction at petascale with the atmospheric model NUMACastroM.Ferreiro FerreiroA.García-RodriguezJ.González-VidaJ.MacíasJ.ParésC.Elena Vázquez-CendónM.The numerical treatment of wet/dry fronts in shallow flows: Application to one-layer and two-layer systemsBunyaS.KubatkoE.J.WesterinkJ.J.DawsonC.A wetting and drying treatment for the Runge–Kutta discontinuous Galerkin solution to the shallow water equationsKesserwaniG.LiangQ.Well-balanced RKDG2 solutions to the shallow water equations over irregular domains with wetting and dryingBerthonC.MarcheF.TurpaultR.An efficient scheme on wet/dry transitions for shallow water equations with frictionLeH.LambrechtsJ.OrtlebS.GratiotN.DeleersnijderE.AoS.F.An implicit wetting–drying algorithm for the discontinuous Galerkin method: Application to the Tonle Sap, Mekong River BasinPopeS.B.LinP.LiuP.A numerical study of breaking waves in the surf zoneLinP.LiuP.Turbulence transport, vorticity dynamics, and solute mixing under plunging breaking waves in surf zoneHorrilloJ.WoodA.KimG.ParambathA.A simplified 3-D/Navier–Stokes numerical model for landslide tsunami: Application to the Gulf of MexicoZhangY.YeF.StanevE.GrashornS.Seamless cross-scale modeling with SCHISM, Ocean ModellingJacobsenN.FuhrmanD.FredsoeJ.A wave generation toolbox for the open-source CFD library: OpenFOAM (R)GhoshD.MittalA.BhattacharyyaS.Multiphase modeling of tsunami impact on building with openingsOwenH.HouzeauxG.SamaniegoC.CucchiettiF.MarinG.TripianaC.HadrienC.VázquezM.Two Fluids Level Set: High Performance Simulation and Post ProcessingMaddenE.BaderM.BehrensJ.van DintherY.GabrielA.A.RannabauerL.UlrichT.UphoffC.VaterS.WollherrS.Linked 3D modeling of megathrust earthquake-tsunami events: From subduction to tsunami run upPelinovsckyE.N.Non-linear hyperbolic equations and run-up of Huge sea wavesRiveraP.Modeling the Asian tsunami evolution and propagation with a new generation mechanism and non-linear dispersive wave modelBabaT.AllgeyerS.HossenJ.CumminsP.R.TsushimaH.ImaiK.YamashitaK.KatoT.Accurate numerical simulation of the far-field tsunami caused by the 2011 Tohoku earthquake, including the effects of Boussinesq dispersion, seawater density stratification, elastic loading, and gravitational potential changeChenX.Dynamic coupling of a three-dimensional hydrodynamic model with a laterally averaged, two-dimensional hydrodynamic modelChenX.J.Coupling an Unstructured Grid three-dimensional Model with a Laterally Averaged two-dimensional Model for Shallow Water Hydrodynamics and Transport ProcessesBungartzH.LindnerF.GatzhammerB.MehlM.ScheufeleK.ShukaevA.UekermannB.preCICE–a fully parallel library for multi-physics surface couplingMehlM.UekermannB.BijlH.BlomD.GatzhammerB.Van ZuijlenA.Parallel coupling numerics for partitioned fluid–structure interaction simulationsLottoG.JeppsonT.DunhamE.Fully Coupled Simulations of Megathrust Earthquakes and Tsunamis in the Japan Trench, Nankai Trough, and Cascadia Subduction ZoneASCETEAdvanced Simulation of Coupled Earthquake and Tsunami EventsAvailable online: https://t3projects.cen.uni-hamburg.de/index.php?id=2099(accessed on 21 December 2020)MarrasS.KoperaM.GiraldoF.X.Simulation of Shallow Water Jets with a Unified Element-based Continuous/Discontinuous Galerkin Model with Grid Flexibility on the SphereMarrasS.KoperaM.ConstantinescuE.SuckaleJ.GiraldoF.A Residual-based Shock Capturing Scheme for the Continuous/Discontinuous Spectral Element Solution of the 2D Shallow Water EquationsTolkovaE.Land–water boundary treatment for a tsunami model with dimensional splittingMacíasJ.CastroM.OrtegaS.EscalanteC.González-VidaJ.Performance Benchmarking of Tsunami-HySEA Model for NTHMP’s Inundation Mapping ActivitiesMacíasJ.CastroM.OrtegaS.González-VidaJ.Performance assessment of Tsunami-HySEA model for NTHMP tsunami currents benchmarking. Field casesMacíasJ.CastroM.EscalanteC.Performance assessment of the Tsunami-HySEA model for NTHMP tsunami currents benchmarking. Laboratory dataMacíasJ.EscalanteC.CastroM.J.Multilayer-HySEA model validation for landslide generated tsunamis. Part I Rigid slidesMacíasJ.EscalanteC.CastroM.J.Multilayer-HySEA model validation for landslide generated tsunamis. Part II Granular slidesImamuraF.ImamuraF.YalcinerA.OzyurtG.YalcinerA.PelinovskyE.ZaytsevA.KurkinA.OzerC.KarakusH.WangX.PranowoW.BehrensJ.SchlichtJ.ZiemerC.Adaptive mesh refinement applied to tsunami modeling: TsunaFLASHRegulyI.Z.GilesD.GopinathanD.QuivyL.BeckJ.H.GilesM.B.GuillasS.DiasF.The VOLNA-OP2 tsunami code (version 1.5)RoelvinkJ.Van BanningG.Design and development of DELFT3D and application to coastal morphodynamicsPopinetS.Basilisk: Simple abstractions for octree-adaptive schemePopinetS.A quadtree-adaptive multigrid solver for the Serre–Green–Naghdi equationsPopinetS.A vertically-Lagrangian, non-hydrostatic, multilayer model for multiscale free-surface flowsRoeberV.CheungK.Boussinesq-type model for energetic breaking waves in fringing reef environmentTavakkolS.LynettP.Celeris: A GPU-accelerated open source software with a Boussinesq-type wave solver for real-time interactive simulation and visualizationKennedyA.ChenQ.KirbyJ.DalrympleR.Boussinesq modeling of wave transformation, breaking and runup, part I: 1DShiF.KirbyJ.HarrisJ.GeimanJ.GrilliS.A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundationLynettP.J.WuT.R.LiuP.L.F.Modeling wave runup with depth-integrated equationsKimD.LynettP.Turbulent mixing and passive scalar transport in shallow flowsYamazakiY.KowalikZ.CheungK.Depth-integrated, non-hydrostatic model for wave breaking and run-upHorrilloJ.GrilliS.NicolskyD.RoeberV.ZhangJ.Performance benchmarking tsunami models for NTHMP’s inundation mapping activitiesMaG.ShiF.KirbyJ.Shock-capturing non-hydrostatic model for fully dispersive surface wave processesVázquezM.HouzeauxG.Alya: Multiphysics Engineering Simulation Towards ExascaleFigures
Pieces of the To¯hoku multi-billion dollar tsunami protection wall after a level-2 tsunami his the coasts of Japan in 2011. While each individual concrete section of the wall did not suffer significant damage, the erosion at the foundations of the wall happened so quickly that the concrete barrier simply fell (Picture taken from in [13].).
Map of the Ring of Fire, a long coastal stretch that is most likely impacted by large tsunamis. Some tsunami mitigation parks are being constructed in South Java, Indonesia (Image: Indonesia Ministry of Marine Affairs and Fisheries), Miyagi Prefecture, Japan (Image: the Morino Project), and Constitución, Chile (Image: Architect Magazine). Adapted from the work in [12].
The off-shore tsunami is, effectively, a two-dimensional long-wave (left, Chilean tsunami, 2010). On-shore, however, it becomes turbulent and multi-phase—water, sand, dirt, and debris (right, To¯hoku tsunami, Japan, 2011). Left: adapted from NOAA Tsunami Warning System. Right: image credit: Reuter/Kyodo found at www.theatlantic.com/photo/2011/03/earthquake-in-japan/100022/.
Two examples of flexible vegetation in a three-dimensional flow at large Reynolds number. Left: bent idealized 1D vegetation submerged in a channel flow (adapted from the work in [88]). Right: two-way fluid–structure interaction model of 3D idealized trees in a dam break triggered flow (image courtesy of Abhishek Mukherjee).
Representation of a comprehensive off-to-on-shore flow framework. The 3D Navier–Stokes solver is forced at the boundaries by the shallow water model of the tsunami off-shore, which is forced by an earthquake simulator. In fact, the 2D shallow water solver can extend all the way into the shore depending on the needs of the solver and compatibility layer.
Gravity-centered, depth-averaged coordinate system for the shallow water equations. Here h is the depth of the water column; η the difference between a defined datum, commonly a given sea-level, and the sea-surface; and b the bathymetry surface as measured from the same datum.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.