# Probabilistic Tsunami Hazard and Risk Analysis: A Review of Research Gaps

Jörn Behrens, Finn Løvholt, Fatemeh Jalayer, Stefano Lorito, Mario Salgado-Gálvez, Mathilde Sørensen, Stéphane Abadie, Ignacio Aguirre-Ayerbe, Iñigo Aniel-Quiroga, Andrey Babeyko, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Jörn Behrens, Finn Løvholt, Fatemeh Jalayer, Stefano Lorito, Mario Salgado-Gálvez, et al. Probabilistic Tsunami Hazard and Risk Analysis: A Review of Research Gaps. Frontiers in Earth Science, 2021, 9, 10.3389/feart.2021.628772. [https://uca.hal.science/hal-03227520](https://uca.hal.science/hal-03227520)

HAL Id: hal-03227520

https://uca.hal.science/hal-03227520

Submitted on 17 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for
updates

# Probabilistic Tsunami Hazard and Risk Analysis: A Review of Research Gaps

Jörn Behrens¹*, Finn Løvholt², Fatemeh Jalayer³, Stefano Lorito⁴,
Mario A. Salgado-Gálvez⁵,⁶, Mathilde Sørensen⁷, Stephane Abadie⁸,
Ignacio Aguirre-Ayerbe⁹, Iñigo Aniel-Quiroga⁹, Andrey Babeyko¹⁰, Marco Baiguera¹¹,
Roberto Basili⁴, Stefano Belliazzı³, Anita Grezio¹², Kendra Johnson¹³, Shane Murphy¹⁴,
Raphaël Paris¹⁵, Irina Rafliana¹⁶,¹⁷, Raffaele De Risi¹⁸, Tiziana Rossetto¹¹, Jacopo Selva¹²,
Matteo Taroni⁴, Marta Del Zoppo³, Alberto Armigliato¹⁹, Vladimir Bureš²⁰, Pavel Cech²⁰,
Claudia Cecioni²¹, Paul Christodoulides²², Gareth Davies²³, Frédéric Dias²⁴,
Hafize Başak Bayraktar³, Mauricio González⁹, Maria Gritsevich²⁵,²⁶,²⁷, Serge Guillas¹¹,
Carl Bonnevie Harbitz², Utko Kânoğlu²⁸, Jorge Macías²⁹, Gerassimos A. Papadopoulos³⁰,
Jascha Polet³¹, Fabrizio Romano⁴, Amos Salamon³², Antonio Scala³, Mislav Stepinac³³,
David R. Tappin¹¹,³⁴, Hong Kie Thio³⁵, Roberto Tonini⁴, Ioanna Triantafyllou³⁶,
Thomas Ulrich³⁷, Elisa Varini³⁸, Manuela Volpe⁴ and Eduardo Vyhmeister³⁹

OPEN ACCESS

**Edited by:**

Victoria Miller,
The University of the West Indies St.
Augustine, Trinidad and Tobago

**Reviewed by:**

Eric Geist,
United States Geological Survey
(USGS), United States
Patricio Andres Catalan,
Federico Santa María Technical
University, Chile

*Correspondence:

Jörn Behrens,
joern Behrens@uni-hamburg.de

**Specialty section:**

This article was submitted to
Geohazards and Georisks,
a section of the journal
Frontiers in Earth Science

**Received:** 12 November 2020

**Accepted:** 10 February 2021

**Published:** 29 April 2021

**Citation:**

Behrens J, Løvholt F, Jalayer F, Lorito S, Salgado-Gálvez MA, Sørensen M, Abadie S, Aguirre-Ayerbe I, Aniel-Quiroga I, Babeyko A, Baiguera M, Basili R, Belliazzı S, Grezio A, Johnson K, Murphy S, Paris R, Rafliana I, De Risi R, Rossetto T, Selva J, Taroni M, Del Zoppo M, Armigliato A, Bureš V, Cech P, Cecioni C, Christodoulides P, Davies G, Dias F, Bayraktar HB, González M, Gritsevich M, Guillas S, Harbitz CB, Kânoğlu U, Macías J, Papadopoulos GA, Polet J, Romano F, Salamon A, Scala A, Stepinac M, Tappin DR, Thio HK, Tonini R, Triantafyllou I, Ulrich T, Varini E, Volpe M and Vyhmeister E (2021) Probabilistic Tsunami Hazard and Risk Analysis: A Review of Research Gaps. *Front. Earth Sci.* 9:628772. doi: 10.3389/feart.2021.628772

¹Department of Mathematics/CEN, Universität Hamburg, Hamburg, Germany, ²NGI - Norwegian Geotechnical Institute, Oslo, Norway, ³University of Naples Federico II, Naples, Italy, ⁴Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy, ⁵ERN Internacional, Mexico City, Mexico, ⁶Centre Internacional de Metodes Numerics en Enginyeria (CIMNE), Barcelona, Spain, ⁷University of Bergen, Bergen, Norway, ⁸Université de Pau et des Pays de L'Adour, E2S UPPA, SIAME, France, ⁹ICantabria - Instituto de Hidráulica Ambiental de La Universidad de Cantabria, Santander, Spain, ¹⁰Deutsches GeoForschungsZentrum GFZ, Potsdam, Germany, ¹¹University College London, London, United Kingdom, ¹²Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Bologna, Italy, ¹³GEM Foundation, Pavia, Italy, ¹⁴Ifremer, Plouzané, France, ¹⁵Laboratoire Magmas and Volcans, Aubière, France, ¹⁶German Development Institute/Deutsches Institut für Entwicklungspolitik (DIE), Bonn, Germany, ¹⁷Indonesian Institute of Sciences (LIPI), Jakarta, Indonesia, ¹⁸University of Bristol, Bristol, United Kingdom, ¹⁹Alma Mater Studiorum - University of Bologna, Bologna, Italy, ²⁰University of Hradec Kralove, Hradec Kralove, Czech Republic, ²¹Università Degli Studi Roma Tre, Rome, Italy, ²²Cyprus University of Technology, Limassol, Cyprus, ²³Geoscience Australia, Canberra, ACT, Australia, ²⁴University College Dublin, Dublin, Ireland, ²⁵Finnish Geospatial Research Institute (FGI), Masala, Finland, ²⁶Department of Physics, University of Helsinki, Helsinki, Finland, ²⁷Institute of Physics and Technology, Ural Federal University, Ekaterinburg, Russian Federation, ²⁸Department of Engineering Sciences, Middle East Technical University, Ankara, Turkey, ²⁹Facultad de Ciencias, Universidad de Málaga, Málaga, Spain, ³⁰International Society for the Prevention and Mitigation of Natural Hazards, Athens, Greece, ³¹Cal Poly Pomona, Pomona, CA, United States, ³²Geological Survey of Israel, Jerusalem, Israel, ³³University of Zagreb, Zagreb, Croatia, ³⁴British Geological Survey, Nottingham, United Kingdom, ³⁵AECOM, Los Angeles, CA, United States, ³⁶Department of Geology and Geoenvironment, National and Kapodistrian University of Athens, Athens, Greece, ³⁷Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Munich, Germany, ³⁸C.N.R. - Institute for Applied Mathematics and Information Technologies, Milano, Italy, ³⁹Insight Research Centre, University College Cork, Cork, Ireland

Tsunamis are unpredictable and infrequent but potentially large impact natural disasters. To prepare, mitigate and prevent losses from tsunamis, probabilistic hazard and risk analysis methods have been developed and have proved useful. However, large gaps and uncertainties still exist and many steps in the assessment methods lack information, theoretical foundation, or commonly accepted methods. Moreover, applied methods have very different levels of maturity, from already advanced probabilistic tsunami hazard analysis for earthquake sources, to less mature probabilistic risk analysis. In this review we give an overview of the current state of probabilistic tsunami hazard and risk analysis. Identifying research gaps, we offer suggestions for future research directions. An extensive literature list allows for branching into diverse aspects of this scientific approach.

**Keywords:** tsunami, probabilistic method, hazard, risk, research gap
# INTRODUCTION

Tsunamis are rare but potentially devastating natural hazards. With often limited available data, a coherent framework that incorporates data, physical assumptions (i.e., the general model of the system), and statistical methods for hazard and risk analysis is necessary to assess consequences affecting different layers of societies. To further develop, standardize and document such a framework is the underlying objective of COST Action AGITHAR (Accelerating Global Science in Tsunami Hazard and Risk Analysis; AGITHAR, 2020) and this article forms one outcome of the Action.

Probabilistic tsunami hazard and risk analyses (PTHA and PTRA, respectively) offer structured and rigorous procedures that allow for tracing and weighting the key elements in understanding the potential tsunami hazard and risk in globally distributed applications (e.g., Basili et al., 2021). Because of this, PTHA are becoming a standard basis for tsunami risk assessment around the world. Significant challenges in this analysis method are 1) the choice of hypothetical events and assigning “correct” probabilities, and ii) the impact of source regions distributed throughout an ocean basin and, conceivably, unifying distinct types of sources in a homogeneous probabilistic framework with a comprehensive treatment of uncertainty. The great importance of PTHA is due to its practical implications for society providing information for long-term planning and coastal management in areas where potential tsunamis may occur. Conversely, PTRA are still less abundant and standardized than PTHA, as elaborated in this review.

Few mega-tsunamis have been observed in the instrumental period, a timeframe spanning from approximately the 1960s to today. Thus, it is challenging to confidently assess the rate at which consequential tsunamis will occur. Predominantly seismically triggered tsunamis comprise about 80% of all tsunamis worldwide (e.g., Harbitz et al., 2014) with the remainder caused by landslides, volcanoes, or meteorological phenomena.

The sparsity of background data and requirements in engineering applications have driven the development of probabilistic methods for assessing tsunami hazard and risk aiming for unbiased comparisons of different hazards (natural and anthropogenic) as well as their uncertainty quantification. In recent years, the probabilistic framework has been increasingly applied. However, broadly accepted approaches are not yet defined, and potentially incompatible implementations of probabilistic methods are used in different regions across the world, and different tsunami source types are often treated separately and are rarely combined.

In this study, we have documented current gaps and open research questions related to PTHA and PTRA. We have organized this review into two main sections, one focused on tsunami hazard and the other on risk. We preface these topics with a brief introduction to the probabilistic framework underlying both PTHA and PTRA. Note that we grouped the gaps in numerical modeling in the hazard analysis related section, even though modeling may also be considered a cross-cutting

topic. We believe, however, that the mentioned gaps are more related and addressed in a similar way as the other hazards related research gaps.

While PTHA and PTRA allow for including uncertainty in a consistent way, it is necessary to point out that it is not always simple to describe the knowledge gaps formally, for example through alternative models, and quantify their impact on hazard and risk models in terms of epistemic uncertainty (i.e., caused by lack of knowledge or data, Kiureghian and Ditlevsen, 2009). Overall, the research gaps identified in this study are “known unknowns” (e.g., Logan, 2009) and deserve more thorough research efforts, in order to determine their influence on the overall outcome of the PTHA or PTRA workflow.

This fact makes it hard to determine quantitatively the importance of each of the research gaps. Nevertheless, we tried to assess—in a qualitative way—the relative priority of research gaps and discuss this in the last section of this report.

# PROBABILISTIC FRAMEWORK

In this section, we present a structure for probabilistic hazard and risk analyses. An overview is given in **Figure 1**. More in-depth reviews of identified gaps related to the individual probabilistic framework components are discussed in sections “Probabilistic Tsunami Hazard Analysis” and “Probabilistic Tsunami Risk Assessment”.

The purpose of PTHA is to find the probability for a tsunami intensity measure (IM) to exceed a given threshold in a predefined time interval. Note that, in the PTHA framework, “Intensity Measure” is used with a meaning that differs from the “tsunami intensity scale” used, for instance, in tsunami catalogs to define the “size” of a tsunami or the effects it produces inland. In the PTHA context, an IM is a physical observable strictly connected to the physics of the process. Common IMs are wave amplitude, flow depth, current velocity, momentum flux, or maximum inundation height, depending on the problem setting (Grezio et al., 2017).

Different probabilistic framework alternatives for computational PTRA exist. One option, rooted in seismic risk analysis, is performance-based risk assessment, presented by PEER (Pacific Earthquake Engineering Center) in 2000. The term performance-based is often used in contraposition to traditional prescriptive assessment procedures for seismic-resistant building design (Fardis, 2009). The performance-based framework aims to provide a practical yet rigorous workflow and has also been used for risk assessment for hurricanes (e.g., van de Lindt and Dao, 2009; Barbato et al., 2013), floods (De Risi et al., 2013; Jalayer et al., 2016), and tsunamis (Chock et al., 2011; Chock, 2016; Attary et al., 2017). This framework can be organized in different modules; for example, hazard and vulnerability or hazard, fragility and consequence. Modules communicate with each other through intermediate variables and their conditional probabilities. Examples of intermediate variables are intensity measure (IM), damage measure (DM) and decision variable (DV). IM serves as an intermediate variable between hazard and vulnerability. DM
<figure><img src="image_4.png" /><figcaption>FIGURE 1 | Roadmap of PTHA and PTRA frameworks: The entire process of risk evaluation needs to interact with the (risk-informed) decision-making process. Composite multi-dimensional risk and vulnerability indicators ("Probabilistic Tsunami Hazard Analysis" section) are shown as defining the context for the complex tsunami risk evaluation. The exposure modeling ("Probabilistic Tsunami Risk Assessment" section) defines groups of individuals and assets at risk. The horizontal flowchart at the bottom of the figure shows the PEER-like workflow for risk assessment. Probabilistic hazard analysis ("Probabilistic Tsunami Hazard Analysis" section) discusses estimation of the mean annual frequency (rate) of exceedance of a given value (im) of an intensity measure (IM, Eq. 1) commonly visualized as a hazard curve. The IM can be a vector or a scalar that describes the intensity of a tsunami. Examples of IM's are flow depth, maximum tsunami inundation height, etc. M refers in a generic manner to the size of various tsunami sources (e.g., earthquake magnitude, landslide volume). The tsunami sources, probability and modeling (earthquake, landslide, volcanic and meteotsunami) are discussed in "Probabilistic Tsunami Hazard Analysis" section. s denotes the vector of source parameters. N denotes the number of tsunamigenic sources. $\lambda_{M_{min}}$ denotes the mean annual frequency of tsunamigenic events exceeding $M_{min}$ from source i. "Gaps in Hydrodynamic Tsunami Modeling, Generation, Propagation, and Run-up" section discusses hydrodynamic tsunami modeling, generation, propagation and run-up. The physical vulnerability ("Gaps in Physical Vulnerability" section) discusses the estimation of the probability distribution for a damage measure (DM, specific value dm) given IM (specific value im), known as the fragility function. The most common example of a DM is the physical damage state. The risk and resilience metrics section ("Gaps in Risk and Resilience Metrics" section) discusses the estimation of various decision variables (e.g., fatalities, repair costs, downtime) denoted as DV (specific value dv). More specifically, it discusses the probability distribution for DV given DM also known as the consequence function. The vulnerability function (Eq. 3) describes the (mean and standard deviation) of the probability distribution for DV given IM and is obtained by integrating over the entire domain of DM. One way to show the PTRA results is through visualizing the mean annual frequency of exceeding a specific value dv of DV (e.g., the loss exceedance curve (LEC) or the annual average loss (AAL)) shown in Eq. 2, referred to generically as the risk curve.</figcaption></figure>

connects vulnerability with fragility and describes physical damage. DV connects fragility with consequences and reaches out to decision-makers with numbers of casualties, repair costs, or downtimes. Interestingly, several risk-informed decision-making processes related to tsunamis are based on PTHA information only (e.g., hazard-based evacuation zones, hazard-based land-use planning). As an example, the criterion "flow depth (IM) larger than a threshold (im)" can be used as a basis for decision-making (e.g., assigning evacuation zones). In other words, an IM can act as an intermediate variable (intensity measure) as well as a decision variable.

To illustrate the framework, suppose a finite set of N hypothetical tsunamigenic sources representing all possible tsunami events affecting the site of interest. Each event occurs randomly in time and independently of all others (i.e., as a Poisson process). The tsunami hazard curve—the main outcome of PTHA—describes the mean annual rate of a

tsunami event affecting location x with an intensity measure IM(x) greater than some threshold im, denoted as $\lambda(\text{IM}(x) \ge \text{im})$. This can be expressed as:

$$ \lambda(\text{IM}(x) \ge \text{im}) = \sum_{i}^{N} \lambda_{M_{\text{min},i}} \int P(\text{IM}(x) \ge \text{im} | s, m) f_{s|M}(s|m) f_{M}(m) \, \text{ds}dm \quad (1) $$

where $\lambda_{M_{min,i}}$ is the mean annual rate of occurrence of tsunamigenic events from source i (e.g., earthquakes, landslides, etc.) having magnitudes M exceeding $M_{min}$, $f_M$ is the conditional probability density function for $M \ge M_{min,p}$ and $f_{s|M}$ is the probability density function of the set of source parameters S given magnitude M. The aleatoric uncertainty associated with variable source characteristics can be represented by probabilistic prediction models of the source
parameters. Finally, $p(\text{IM}(x)\ge im|s,m)$ is the complementary cumulative distribution function of IM given $S = s$ and $M = m$, and can be evaluated through tsunami simulations. Note that Eq. 1 can be used only if sources are independent; a counterexample being a landslide generated from the same earthquake that amplifies the ensuing tsunami's destruction.

Epistemic uncertainty in PTHA is often accounted for using logic trees or, more recently, ensemble modeling, which allow alternative hypotheses for uncertain parameters, each of which is assigned a weight reflecting confidence in the respective parameter value (e.g., Geist and Parsons, 2006; Selva et al., 2016; Grezio et al., 2017). Equation 1 is computed for each logic tree 'end branch'.

Building on tsunami hazard, the tsunami loss curve at any location is calculated by convolving vulnerability and hazard over the entire span of IM:

$$ \lambda(DV \ge dv) = \int_{im} G_{DV|IM}(dv|im) d\lambda (IM(x) \ge im) \quad (2) $$

where $\lambda(DV \ge dv)$ is the mean annual rate of occurrence of DV larger than a threshold $dv$. Vulnerability is expressed through the complementary probability distribution function denoted as $G_{DV|IM}(dv|im)$, for DV given IM, and is itself calculated by integrating fragility and consequence functions (see also Figure 1):

$$ G_{DV|IM}(dv|im) = \int_{dm} G_{DV|DM}(dv|dm) f_{DM|IM}(dm|im) ddm, \quad (3) $$

with

* $f_{DM|IM}$, the tsunami fragility function, predicts the probability of incurring a particular value (dm) of damage measure DM (e.g., damage states) for a given IM = im;

* $G_{DV|DM}(dv|dm)$, the tsunami consequence function (e.g., the damage-loss function), expressed as the complementary cumulative distribution function of DV given DM.

Strictly speaking, Eqs. 1 and 2 do not consider multi-hazard and multi-risk aspects such as cascading effects, combined damage due to tsunami loading and earthquake shaking. Assuming a Poisson process, the rate of exceedance $\lambda$ is often transformed the first excursion of a specific value $dv$ for a generic decision variable DV in the time $\Delta t$ (e.g., 1 year, 50 years):

$$ P(DV \ge dv; \Delta t) = 1 - \exp(-\lambda(DV \ge dv)\Delta t) \quad (4) $$

## PROBABILISTIC TSUNAMI HAZARD ANALYSIS

This section discusses gaps in PTHA, focusing on those in tsunami sources and hydrodynamic modeling. Each subsection includes a summary of the present state-of-the-art, followed by an in-depth discussion of the gaps.

## Gaps in Earthquake Source Representation Existing Methods

Seminal Seismic PTHA (SPTHA) was performed using crude source and tsunami representations (Lin and Tung, 1982; Rikitake and Aida, 1988; Tinti, 1991). Since then, the methodology has evolved dramatically (Geist and Parsons, 2006; Annaka et al., 2007; Power et al., 2007; Thio et al., 2007; Burbidge et al., 2008; González et al., 2009; Sørensen et al., 2012; Hoechner et al., 2016; Miyashita et al., 2020), also in the framework of large programs (e.g., Horspool et al., 2014; Davies et al., 2018; Davies and Griffin, 2018; Basili et al., 2021).

SPTHA methodology for spatio-temporal and kinematic source treatment and the basic uncertainty framework were mostly transcribed from Probabilistic Seismic Hazard Analysis (PSHA, Esteva, 1967; Cornell, 1968; a historical perspective: McGuire, 2008). Due to tsunami data scarcity, it is challenging to derive hazard estimates directly from historical records (Geist and Parsons, 2006; Grezio et al., 2017). Consequently, numerical modeling is a distinctive characteristic of SPTHA where seafloor displacement and tsunami evolution from generation to inundation are simulated for each scenario (Geist and Parsons, 2006; Geist and Lynett, 2014). Source parameters can be inferred from past seismicity or from balancing the seismic moment across a fault zone, potentially constrained by geodetic strain rates (Grezio et al., 2017). Often only major subduction zones are considered in SPTHA, assuming that they are the main hazard drivers (e.g., González et al., 2009; Davies et al., 2018). In this case, spatial characterization provides geometrical and kinematic constraints, such as the fault geometry, preferential slip direction, and other source zone properties. Crustal and general seismicity from unconstrained or unknown faults is treated with a larger uncertainty (e.g., Selva et al., 2016; Basili et al., 2021). Earthquakes are usually simplified to having either uniform (e.g., Horspool et al., 2014) or heterogeneous instantaneous slip (e.g., De Risi and Goda, 2017). Seafloor deformation is predominantly computed analytically assuming an elastic homogeneous half-space (Mansinha and Smylie, 1971; Okada, 1992; Meade, 2007; Nikkhoo and Walter, 2015).

State-of-the-art seismic source representation for tsunami simulations is reviewed by Geist and Oglesby (2014) and Geist et al. (2019). Additionally, we note some innovative efforts for complex, yet computationally affordable, approaches to source simulation (Melgar et al., 2016; Murphy et al., 2016; Sepúlveda et al., 2017; Scala et al., 2020), and methods for handling source modeling uncertainties and sensitivity including temporal aspects and recurrence (Grezio et al., 2010; Basili et al., 2013; Lorito et al., 2015; Selva et al., 2016; Lotto et al., 2017; Davies, 2019; Goda, 2019; Davies and Griffin, 2020).

## Identified Gaps

*Limited Past Events and Data to Inform Hazard Models (S1)* Completeness and quality of historical earthquake data, needed to constrain seismic source parameters, varies greatly depending on the history of the investigated geographical region (Stucchi et al., 2004; Albini et al., 2014). Event catalogs are generally too short to account for the source frequency needed to model large average return periods in PTHA. The description of earthquake
recurrence and of the tail of the frequency-magnitude distribution is highly uncertain (Kagan, 2002; Geist and Parsons, 2014; Rong et al., 2014; Bommer and Stafford, 2016). In the attempt of constraining this uncertainty, seismic source parameters have been estimated globally using seismic or geodetic data or both (e.g., Bird and Kagan, 2004; Bird et al., 2015; Bird and Kreemer, 2015); however, these types of input data are not always considered by PTHAs. Moreover, a framework for constraining PTHA directly from tsunami observations exists (Geist and Parsons, 2006; Grezio et al., 2017), while treatment of incomplete catalogs is described by Smit et al. (2017). Where possible, other data types should also be considered. Paleo-seismic and paleo-tsunami catalogs may help constrain or validate at least large event recurrence (e.g., Priest et al., 2017; Paris et al., 2020), while GPS-constrained strain accumulation can indicate the total seismic moment rate (e.g., Hayes et al., 2018). Care should be taken of potential biases coming from overweighting evidence of large tsunamis in the past (Geist and Parsons, 2006).

## Fault Identification, Fault and Source Zone Parameterization and Tsunamigenic Potential Characterization (S2)

Tsunami sources are often constrained from infrequent offshore geologic studies investigating very large areas. Therefore, geologic fault data are often incomplete, causing a wide range of source knowledge levels (Basili et al., 2013). Seismic source characterization for SPTHA generally refers to properties of pre-existing large faults, and often only to great subduction zone sources. All other-mostly crustal-faults are seldomly considered in PTHA, although non-subduction earthquakes may control tsunami hazard, especially when located near the target site (Selva et al., 2016). Despite overall good constraint of subduction interface geometries (e.g., Hayes et al., 2018), along-strike trench segmentation and its impact on rupture propagation remains uncertain, limiting rupture forecasts and hindering estimates of subduction earthquake maximum magnitude (e.g., Bilek, 2010; Kopp, 2013; Grezio et al., 2017). Whenever fault knowledge is incomplete, more randomized "background" seismicity modeling is needed, with less predictable geometry and seismic behavior compared to subduction interfaces (Sørensen et al., 2012; Selva et al., 2016). Fault slip rates can constrain seismicity recurrence parameters; these can vary both spatially (Zechar and Frankel, 2009) and temporally (e.g., Ota and Yamaguchi, 2004; Ramírez-Herrera et al., 2011; Tiberti et al., 2014) but usually only averages are considered due to a lack of information. Kagan and Jackson (2014) pointed out that more research would be needed for focal mechanism forecasting; identifying the prevailing faulting mechanism is a critical task particularly in tectonically complex environments. This is expected, in turn, to exert a strong influence on tsunami hazard.

## Variety, Complexity, and Dynamics of Fault Mechanics (S3)

Source simplification represents a dominant uncertainty in SPTHA (Geist and Oglesby, 2014). Its effect on seafloor deformation needs to be investigated better, concerning deformation models that incorporate complex material properties, geometrical complexity, varying depth-dependent

fault conditions, dynamic simulations including off-fault damage and near-surface amplification, which all may increase tsunami hazard (Masterlark, 2003; Ma, 2012; Kozdon and Dunham, 2013; Ryan et al., 2015; Murphy et al., 2016; Lotto et al., 2017; Murphy et al., 2018; Scala et al., 2019; Scala et al., 2020; Tonini et al., 2020). Secondary ruptures including splay faulting may happen as an independent source or as part of a large earthquake on the subduction interface (Wendt et al., 2009; Li et al., 2014; Hananto et al., 2020).

Tsunami earthquakes produce excessively large tsunami intensities compared to their moment magnitude (Polet and Kanamori, 2016), and their global and local frequency is unconstrained. A simplified characterization of tsunami earthquakes, which is sometimes adopted, assumes larger slip associated with less rigid materials at shallow depths to preserve the seismic moment (e.g., Bilek and Lay, 1999; Geist and Bilek, 2001). These and other very complex ruptures, potentially containing fault branching, rupture jumping, and mixed-mode slip (e.g., Ulrich et al., 2019a; Ulrich et al., 2019b), are not well represented in PTHA. On a larger scale, rupture area may be shared by more than one subduction interface, like in the case of triple junctions (e.g., Solomon event 2007, Lorito et al., 2016). Due to a lack of observations the likelihood of such events is uncertain and quantification of their relative contribution to SPTHA therefore difficult.

Due to all these uncertainties and the extreme computational demand for dynamic computation, numerical simulations are *de facto* replaced with heterogeneous stochastic slip modeling (e.g., Herrero and Bernard, 1994; Mai and Beroza, 2002; Davies et al., 2015; Sepúlveda et al., 2017), and less frequently with stochastic stress modeling (e.g., Wendt et al., 2009). Because source observations are relatively scarce, more statistical tests (Davies and Griffin, 2019) are needed for source model validation.

## Empirical Scaling Relations (S4)

Several different empirical scaling relations are used to define earthquake rupture properties, such as length, width, average slip, and earthquake magnitude (e.g., Wells and Coppersmith, 1994; Murotani et al., 2008; Blaser et al., 2010; Strasser et al., 2010; Murotani et al., 2013; Goda et al., 2016; Skarlatoudis et al., 2016; Allen and Hayes, 2017; Thingbaijam et al., 2017). These relationships quantify appreciable uncertainties that are seldomly accounted for in SPTHA. These relations imply stress drop and time-dependent rupture characteristics and self-similarity of earthquakes across scales, but this is apparently violated in some cases. For example, the 2011 Tohoku earthquake released a huge amount of slip in a relatively small portion of the subduction interface compared to the Sumatra 2004 or Chile 1960 earthquakes (Okal, 2015); scaling relations are not directly applicable to abnormally slow and unusually large shallow slip occurring in low-rigidity zones during tsunami earthquakes.

## Complex, Non-stationary Seismic Cycle (S5)

Even in the simplest subduction environment, the seismic cycle over co-seismic, inter-seismic and post-seismic phases is complex and non-stationary, for example due to visco-elastic rheology and
the role of fluids (Wang et al., 2012; Moreno et al., 2014; Melnick et al., 2017). Time-dependent models could potentially be used to estimate the stress transfer from one earthquake to the neighboring faults (King et al., 1994). Stress transfer from megathrust earthquakes triggering outer-rise ruptures or possibly even the opposite are such examples (e.g., Lorito et al., 2016). Based on seismic catalogs, it is possible to infer non-Poissonian earthquake recurrence, for example earthquake clustering (Kagan and Jackson, 1991). A time-dependent model, which could better describe the probability of earthquake occurrence for some specific applications or timeframes, is taken into account by only a few PTHAs (e.g., Goda et al., 2017; Goda, 2020).

## Other Constraints (S6)

It is reasonable to assume that high seismic coupling correlates with future slip location. Under simplifying assumptions, along-strike geodetic coupling variation can be inferred from geodetic strain (Métois et al., 2012). Large uncertainty remains, particularly regarding the near-trench region (Loveless and Meade, 2011). Recent developments in seafloor geodesy and modeling techniques are offering improved constraints (e.g., increasing offshore coupling resolution, Bürgmann and Chadwell, 2014; Foster et al., 2020), and slow slip events and consequently the stress evolution on the fault (e.g., Araki et al., 2017). High seismic coupling combined with stress accumulation in areas of seismic inactivity is described as a seismic gap. The possibility of using seismic gaps to identify zones of enhanced seismic hazard has long been debated (e.g., Bilek and Lay, 2018). Attempts to obtain physically motivated constraints on the maximum magnitude utilizing convergence rate, age of the oceanic crust and sediment thickness have been rather unsuccessful (Okal, 2015). Ongoing research explores these and other controlling factors of subduction zone seismicity, including small- and large-scale roughness of the subduction interface, static friction coefficient, upper plate strain and rigidity, dip angle and curvature (e.g., Heuret et al., 2012; Bletery et al., 2016; Sallarès and Ranero, 2019; Rijsingen, et al., 2019; Muldashev and Sobolev, 2020). Additionally, rupture cycles and supercycles over multiple segments controlled by geological asperities have been proposed (Philibosian and Meltzner, 2020). Similar to some of the previously discussed items in this section, no consensus has been reached on the statistical meaning of such information and on how to frame it within SPTHA.

## Gaps in Landslide Source Representation
Existing Methods

Landslide tsunami PTHA (LPTHA) was introduced less than a decade ago (Geist and Lynett, 2014). Its application is often similar to SPTHA (e.g., ten Brink et al., 2006; Lane et al., 2016), but can also be based on geotechnical interpretations with a strong emphasis on expert judgment (e.g., Grilli et al., 2009; Hermanns et al., 2013; Løvholt et al., 2020). Salamon and Di Manna (2019) derive empirical scaling relations for landslides triggered by onshore earthquakes. In LPTHA, the landslide volume is used analogously to the seismic moment in SPTHA

as a rate of occurrence. The slide volume is generally also the most influential factor on tsunami genesis (Snelling et al., 2020). Landslide motion has a strong influence too (Løvholt et al., 2015b; Yavari-Ramshe and Ataie-Ashtiani, 2016). LPTHA source models are coupled to numerical tsunami models in Monte Carlo simulations. Methods for simulating both the landslide dynamics and tsunami generation range from block models (Harbitz, 1992; Tinti et al., 1997; Watts, 2000; Grilli and Watts, 2005; Tinti et al., 2006; Løvholt et al., 2015b), depth-averaged rheological models of viscoplastic or granular nature (e.g., Kelfoun and Druitt, 2005; Jop et al., 2006; Løvholt et al., 2017; Kim et al., 2019), to computational fluid dynamics (CFD) based approaches with different landslide complexity (e.g., Crosta et al., 2016; Abadie et al., 2020). Submarine landslide tsunamis are mainly characterized by the Froude number (landslide velocity over wave celerity) measuring the degree of critical landslide velocity, the landslide acceleration, and the rate of landslide mass mobilization (e.g., Ward, 2001; Løvholt et al., 2015b). Subaerial landslides are characterized by the landslide frontal area, along with the Froude number, landslide density, and slope angle (e.g., Fritz et al., 2003; Heller and Hager, 2010; Bullard et al., 2019).

## Identified Gaps

### Lack of Understanding and Likelihoods for Tsunamigenic Landslide Volumes (L1)

For submarine landslides, we refer to the reviews of Huhn et al. (2019) and Harbitz et al. (2014). The challenge can be attributed to several factors:

* Limited or insufficient mapping of past landslide occurrences. Their characteristics and lack of dating prevent constraining the age of the sediments without excessive uncertainty ranges (e.g., Geist et al., 2013). The new global landslide database initiative (Clare et al., 2019) is a good starting point for standardizing, but not yet complete enough for feeding LPTHA. Good data coverage exists for certain regions such as the Mediterranean (Urgeles and Camerlenghi, 2013), Gulf of Mexico (Pampell-Manis et al., 2016) and the US East Coast (Chaytor et al., 2009, Geist et al., 2014).

* Limited understanding of how past landslide recurrence can be projected into the future hazard, including time and geological context dependency. For example, we cannot yet generally link climatically driven trends to past landslide frequency (Urlaub et al., 2013). However, it is concluded that the last ice age affect present landslide probability offshore US (Lee, 2009) and Norway (Bryn et al., 2005).

* Limited available geological and geotechnical data inhibit identification of failure-prone sediments and discrimination from stable areas, including weak failure zones, pore pressure conditions or fractures, as well as obstacles or structures. When data exist, they may be proprietary, and a challenge is related to the need for covering very large geographical and heterogeneous regions. A methodological gap exists in bridging geotechnical data and slope stability models (e.g., Carlton et al., 2019) to volume-frequency relationships.
* Limited data and knowledge on triggers of landslides, such as meteorological or seismic events, impedes the quantitative assessment of potential landslide magnitude.

## Difference of Onshore and Offshore Landslides (L2)

The specific character of subaerial and submarine landslides is often vastly different. Potential direct or indirect trigger mechanisms are sometimes not fully understood or difficult to embed into the probability of failure (e.g., precipitation-induced landslides, weak zones and fluid overpressure, range of failure propagation and cascading failure propagation spread). Understanding and estimating the annual probability of landslide failure in rock slopes with complex fracture patterns and stress conditions is associated with large uncertainty. Extensively monitored rock slopes in Norway (e.g., Blika et al., 2005) show large motion over decades before failure takes place, rendering assessment of failure probability difficult. Matching expert judgment (e.g., Hermanns et al., 2013) to observed landslide magnitude frequency statistics (e.g., Nes, 2018) will help aggregate understanding of landslide frequencies and help link knowledge on failure-prone areas to probability. While epistemic uncertainties in the described situations are large, current LPTHA models do not incorporate them.

## Limited Constraints on Landslide Dynamics and Material Behavior (L3)

The interplay of diverse tsunamigenic landslide parameters makes the generation complex, implying that much less voluminous landslides may be more effective tsunami generators than respectively larger ones. As an example, we note that the approximately 500 km³ Trænadjupet Slide that occurred 4,500 years BP likely produced a moderate coastal impact possibly of just a few meters (e.g., Løvholt et al., 2017), while the 100 times less voluminous 1998 Papua New Guinea landslide induced more than 10 m run-up locally (e.g., Tappin et al., 2008). Because tsunami genesis is tightly linked to landslide acceleration as well as rate of mobilization of the landslide volume (e.g., Løvholt et al., 2005), quantifying the rate and nature of the slope failure is important. Just a few studies discuss the effect of initial failure rate on tsunami generation (e.g., Trapper et al., 2015; Germanovich et al., 2016; Puzrin et al., 2016) and related aspects such as remoulding and cascading failures on the landslide tsunami generation (e.g., Løvholt et al., 2017; Kim et al., 2019; Zengaffinen et al., 2020). How to include these factors and their associated probabilities in PTHA is not resolved. While advanced numerical models for landslide dynamics exist (e.g., Tinti et al., 1997; Jop et al., 2006; Savage et al., 2014; Si et al., 2018a; Si et al., 2018b; Kim et al., 2019; Wang et al., 2019; Gallotti and Tinti, 2020), their complexity and variety obfuscate understanding on which models are most suitable to be used. Furthermore, some models (e.g., Savage et al., 2014; Si et al., 2018a; Si et al., 2018b) are presently too comprehensive to be used in PTHA. Procedures for linking them to measured material properties and geological settings are not in place. Finally, fluid resistance forces (pressure drag, skin friction, and added mass) may be as important as the landslide properties, in particular for

submarine landslides and further investigating physical understanding is necessary to constrain epistemic uncertainty.

## Limited Availability of Benchmarks (L4)

Suitable benchmarks have recently been made available (e.g., Huang and Garcia, 1998; NTHMP, 2018; Kim et al., 2019), but are arguably less mature and fewer than their hydrodynamic modeling counterparts (e.g., Pedersen, 2008; Synolakis et al., 2008). A challenge is a transition from simplified laboratory tests to real-world landslide representation. Moreover, while numerous empirical lab experiments exist, they are significantly influenced by scale effects (Heller, 2011). Neither complex rheological behavior nor real-world complexity is covered in the benchmarks. Complex laboratory experiments (e.g., Rondon et al., 2011) can be used for validating CFD models, but CFD models are presently too computationally expensive for tsunami hazard analysis modeling.

## Limited Past Events to Inform Hazard Models (L5)

Information about past landslides and tsunamis can be used to infer landslide dynamics uncertainty. This can be done using landslide run-out information alone (e.g., Salmanidou et al., 2017), which consequently yields broad epistemic uncertainties in LPTHA. By using tsunami information, such uncertainties can be drastically reduced (e.g., Gylfadóttir et al., 2017; Kim et al., 2019; Løvholt et al., 2020). In practice, however, very few landslide tsunami data are available.

# Gaps in Volcano Source Representation
Existing Methods

Volcanic PTHA, coined VPTHA here, is even less developed than LPTHA (Grezio et al., 2017). Among the few examples are the VPTHA framework developed in Ulvrova et al. (2016) and Paris et al. (2019) for underwater explosions at Campi Flegrei, and Grezio et al. (2020) for pyroclastic flows of Vesuvius. Given that risk reduction measures at volcanoes are often related to the identification of precursory patterns preceding eruptions or to recognizing unrest episodes with increased volcanic activity, the volcanic hazard is often computed conditional to eruptions or unrest, and without an explicit quantification of long-term probability. For example, in Paris et al. (2019), the hazard analysis (Campi Flegrei, Naples, Italy) is confined to conditional tsunami intensity probabilities, due to probabilistic realizations of eruptions with different vent size and location.

## Identified Gaps

### Variety of Potential Volcanic Sources (V1)

Tsunamigenic volcanic events are diverse and they include both eruptive and non-eruptive triggering phenomena, such as underwater explosions, pyroclastic flows, lahars, slope failures, volcanic earthquakes, shock waves from large explosions, and caldera subsidence (Latter, 1981; Kienle et al., 1987; Begét et al., 2005; Day, 2015; Paris, 2015; Grezio et al., 2017). A large range of wave characteristics is typical for volcano tsunamis, even if most such sources are localized and generate mainly short-period waves with greater dispersion and limited far-field effects compared to earthquake-generated tsunamis (e.g., Yokoyama,
1987; Nomanbhoy and Satake, 1995; Le Méhauté and Wang, 1996; Choi et al., 2003; Watts and Waythomas, 2003; Bellotti et al., 2009; Maeno and Imamura, 2011; Ulvrova et al., 2016; Selva et al., 2019, 2020). However, tsunamis are among the farthest propagating volcanic perils, often generating regional impact (e.g., Krakatau, Stromboli, Ischia, etc., see for example Paris et al., 2014; Rosi et al., 2018; Selva et al., 2019; Gallotti et al., 2020). Notably, some of the tsunamigenic volcanic events overlap with those recorded for seismic and landslide tsunami: flank collapse, slope failure and even pyroclastic flows are related to landslides (Løvholt et al., 2015b; Paris, 2015); volcano-tectonic earthquakes occur with high frequency in volcanic areas (Paris, 2015). Such frequency information as well as understanding material properties and transformation during flow should draw upon volcanological expertise. It is often difficult to define a single generation phenomenon since different potentially tsunamigenic processes can occur during the same volcanic episode, especially during large caldera-forming eruptions (Paris, 2015).

## Difficulties in Constraining Recurrence Rates (V2)

Since volcanic tsunami generation is so diverse, constraining recurrence rates for the different source types as eruptive phenomena (Walter et al., 2019), unrest episodes (Tinti et al., 1999; Selva et al., 2020), and triggered subaerial landslides (Selva et al., 2019) is difficult. The integration into a multi-source VPTHA is further complicated by the need for accounting for the complex interdependencies that may exist among the different source mechanisms. The hazard is often nonstationary through time (e.g., Bebbington, 2008; Bebbington, 2010), which represents another challenge.

## Gaps in Modeling Tsunami Generation and Propagation (V3)

Extensive reviews on existing strategies to model volcanic sources are found in Paris, (2015), Grezio et al. (2017) and Paris et al. (2019). Given the complexity, an important part of the hazard analysis is oriented toward understanding the physical mechanism of generation, and how to represent this probabilistically. Similar to landslide generated tsunamis, volcano tsunami modeling suffers from the difficulty of coupling the complex dynamics of the generating event and its interaction with wave propagation. For example, pyroclastic flows are complex, multi-phase phenomena involving the interaction of high-temperature gases and volcanic clasts covering a very large range of granulometric dimensions (Freundt, 2003; Bougouin et al., 2020). This difficulty leads to simplified modeling schemes (e.g., Bevilacqua et al., 2017; Sandri et al., 2018). These simplified strategies may be too reduced for an effective constraint of their tsunami potential (Grezio et al., 2020). Some phenomena may be represented by empirical models (for submarine explosions, see Paris et al., 2019, and for caldera collapse, see Ulvrova et al., 2016). Experimental and numerical simulations coupled with field data increased understanding of the physics and main parameters of volcanic tsunamis (Grezio et al., 2017).

## Lack of Data From the Geological Record (V4)

Tsunami is often not dealt with in the volcanological community, although it may be more fatal than other volcanic hazards such as lava flows or ash falls (Auker et al., 2013; Brown et al., 2017). Consequently, a systematic investigation of tsunami-related data in geological surveys at volcanoes is often missing. Because different volcanic phenomena may trigger tsunamis, even when tsunami data exist, attributing the observation to a specific mechanism is difficult (e.g., Krakatau 1883 eruption: Paris et al., 2014). Therefore, a systematic collection of available volcano-generated tsunami data and linking to potential volcanic generating processes is required. This will imply defining a strategy of tsunami-oriented monitoring around coastal volcanoes. It would be useful to combine such efforts with existing data collections such as the Global Volcanism Program (Global Volcanism Program, 2013).

## Limited Availability of Well Recorded Past Events or Benchmark Studies (V5)

Only a few past events are well constrained in terms of both the source and of the subsequent tsunami (e.g., Unzen 1792, Karymskoye Lake 1996; Montserrat 1997 and 2003, Anak Krakatau 2018; Stromboli 2002 and 2019). The lack of consensus in modeling procedures for each type of tsunamigenic volcanic event, along with the tendency to consider all sources as “unique”, complicates the task of defining benchmarks for volcano tsunamis.

## Gaps in Meteorological Source Representation

### Existing Methods

Meteotsunami PTHA, coined MPTHA here, was developed only recently (see Grezio et al., 2017). A framework for PTHA development is proposed by Geist et al. (2014). The dynamics of meteotsunamis are fairly well-known (e.g., Monserrat et al., 2006; Sibley et al., 2020), related to unusually strong and rapid atmospheric pressure fluctuations and resonance effects causing strong waves closely associated with the behavior of tsunamis. The source mechanisms of meteotsunamis are also well understood (Monserrat et al., 2006; Pattiaratchi and Wijeratne, 2015) with a major driver a Proudman resonance (Proudman, 1929). Because meteotsunamis are strongly linked to (un)favorable combinations of pressure fluctuations, shallow (shelf) bathymetry, and directivity of the weather system, they take place more frequently in specific geographical areas, such as in the Adriatic Sea (Vilibić and Šepić, 2009), the Baltic Sea (Pellikka et al., 2020), and the East Coast of the United States (Pasquet et al., 2013). The main input data for meteotsunamis include meteorological pressure data, preferably with full spatial and temporal characteristics of the pressure field for given meteorological events. Such data can be used to provide synthetic probabilistic source scenarios as input to an PTHA, where an example for the Northeast US coastline is given by Geist et al. (2014). While this field does
not share the data sparsity issues that are associated with volcanoes and landslides, large uncertainties persist, as briefly discussed below.

## Identified Gaps

### Lack of Understanding the Potential and Likelihood for Tsunamigenic Meteorological Patterns (M1)

A systematic assessment of potential source areas and exposed coastal regions is not available. Some regional studies can serve as a preliminary indication (e.g., Dusek et al., 2019; Šepić et al., 2012; Šepić et al., 2016), but a rigorous catalog is missing. Climatological information is likely available, but a systematic extraction of data concerning meteotsunami potential has not been performed. It is not clear whether the resolution of available climatological data (e.g., from reanalysis) is sufficiently fine to allow for the extraction of corresponding relevant meteotsunami source patterns.

### High Sensitivity to Several Parameters and Lack of Understanding of Local Amplification Factors (M2)

Whitmore and Knight (2014) demonstrate the high sensitivity of typical tsunami impact to source parameters and hence a large gap in knowledge on relevant localized parameters. The size, speed, amplitude, directivity, and duration of an atmospheric disturbance resonating with the water column in a specific topographic setting need to be known to assess the hazard. Therefore, such parameters need to be derived for all tsunamigenic regions, then applied to available climatological data sets, and finally fed into corresponding models for assessment of hazard. An assessment of amplifying tidal conditions in each of such regions is also missing.

### Limited Availability of Benchmark Studies (M3)

While there are many individual meteotsunami events described in the literature (e.g., Churchill et al., 1995; González et al., 2001; Pasquet et al., 2013; Vilibić et al., 2014), no truly validated benchmark data are available for meteotsunami benchmarking. In principle, a similar methodology as described in Synolakis et al. (2008) could be followed. However, only very little unification of source modeling has been achieved and except for preliminary simplified tests (as in Vilibić, 2008), there exists no widely accepted test suite. This applies in particular to verification and validation of the probabilistic workflow of MPTHA.

### Limited Past Events and Data to Inform Hazard Models (M4)

There is no consistent catalog of occurrences, although regional studies have been performed (e.g., Haslett et al., 2009; Woodruff et al., 2018). As stated before, there are no unified parameterizations of meteotsunami sources, which could be entered into such a catalog. Even though many individual events are described in the literature (see subsection above), these are by no means representative or complete to be used in hazard models. More rigorous collection of data with the special focus on meteotsunamis-background climatology, meteorological situation, ocean state, topo-bathymetry-for the diverse areas of interest would be desirable.

## Gaps in Hydrodynamic Tsunami Modeling, Generation, Propagation, and Run-up Existing Methods

Hydrodynamic tsunami modeling includes numerical simulation of tsunami generation, propagation as well as coastal and onshore impact. It is an essential part of any PTHA or PTRA analysis. Reviews of commonly applied methods are available (e.g., Pedersen, 2008; Synolakis et al., 2008; Behrens and Dias, 2015). The pre-eminent challenge is the need to bridge a broad range of scales. First, in the probabilistic regime, a comprehensive PTRA must consider a very large number of scenarios to cover all relevant tsunamigenic sources, explore wave physics, and quantify uncertainties. Second, for each individual scenario source, large-scale propagation and coastal inundation modeling (optimally at scales of 1–10 m) need to be represented to quantify tsunami-related on-shore damages and losses. However, the fastest HPC simulation workflows (e.g., de la Asunción et al., 2013; Oishi et al., 2015; Macías et al., 2017; Musa et al., 2018) still require typically 10–60 min to simulate tsunami inundation at a scale of tens of meters, rendering them unsuitable for extensive PTRA studies with up to millions of scenarios (Basili et al., 2021). To overcome this “challenge of scales”, modeling approximations are presently necessary for PTHA feasibility and can either involve 1) largely reducing the number of inundation scenarios (e.g., González et al., 2009; Lorito et al., 2015; Volpe et al., 2019; Williamson et al., 2020), 2) use of approximate models or statistics such as amplification factors (e.g., Løvholt et al., 2012; Kriebel et al., 2017; Gailler et al., 2018; Glimsdal et al., 2019), or 3) machine learning-based tsunami emulators (e.g., Sarri et al., 2012; Salmanidou et al., 2017; Giles et al., 2020).

## Identified Gaps

### PTHA Uncertainty Treatment for Tsunami Inundation Processes (H1)

At present, we lack well tested local PTHA benchmarks where the sources of uncertainties are effectively characterized, in a way that allows their formal propagation along the PTHA/PTRA assessment chain. Moreover, the effect of coseismic coastal displacement due to near field sources (e.g., Volpe et al., 2019), which affects tsunami inundation, should be investigated more deeply, especially when using techniques for reducing the number of scenarios. For this purpose, a large number of inundation scenarios are needed to quantify the epistemic uncertainty and bias caused by simplifications introduced through approximate methods. A local PTHA application using more than 40,000 earthquake sources (Gibbons et al., 2020) is only a start.

### Tsunami Generation (H2)

Unit source models (Kajiura, 1963; Nosov and Kolesov, 2007; Molinari et al., 2016) of varying computational cost and complexity approximate the volumetric deep-water source displacements. While Lotto et al. (2019) clarified that the horizontal momentum does not effectively contribute to tsunami generation in deep-water sources, an extensive sensitivity analysis of how such simplifications affect PTHA
has not been carried out. Incorporating time-dependent and moving sources, be it earthquakes (e.g., Ulrich et al., 2019a), landslides (e.g., Løvholt et al., 2015b) or volcanoes, will involve much higher computational burden. How to limit the number of source time steps for time-dependent source modeling is sparsely studied (e.g., Zengaffinen et al., 2020). For landslide tsunamis, closed-form models (e.g., Watts et al., 2003; Cecioni and Bellotti, 2010) represent a simple alternative but can introduce biases when conveyed to real geographical settings, due to oversimplification or inadequacy for the real situation. Subaerial landslides and volcanoes are often simplified because the required consideration of full 3D hydrodynamics (e.g., Abadie et al., 2020) into PTHA poses too high computational demand. Hence, more research is needed for developing simplified time-dependent sources compatible with PTHA demands, while quantifying the epistemic uncertainty and bias caused by the simplification. New methods may facilitate more detailed characterization of past inundation scenarios and their sources (e.g., Chagué-Goff et al., 2012; Sugawara et al., 2014; Paris et al., 2020).

## Uncertainty and Variability due to Numerical Model for Tsunami Propagation (H3)

Most non-linear shallow water (NLSW) simulation codes produce similar results in the propagation phase, in particular in controlled benchmark cases (e.g., Synolakis et al., 2008). However, clear model differences can appear due to varying components (applied numerical method, workflow, sources, setup etc.) in practical applications. Comparing different numerical forecast models in the Indian Ocean, Greenslade et al. (2014) found large variations, attributed to differences in the workflow and source representation rather than to the tsunami model itself. Testing how such kinds of uncertainty quantification relate to “heterogeneous modeling practices” has not been carried out systematically. Moreover, a rigorous investigation of the performance of far-field propagation is sparse (Dao and Tkalich, 2007; Davies and Griffin, 2020). Differences in numerical dissipation and discretization can also contribute to modeling deviations. As there is no standardized test case for far-field propagation that could reveal the differences in performance of different approaches, it is pressing to address these issues more systematically. Due to the computational burden, most PTHA applications today employ shallow water type models, neglecting frequency dispersion, which can lead to bias. Dispersion can be incorporated through conventional dispersive wave solvers (e.g., Bellotti et al., 2008; Løvholt et al., 2008; Kim et al., 2009; Shi et al., 2012), or through manipulating numerical schemes in NLSW codes (like MOST, e.g., Wei et al., 2008), although the general applicability of the latter is presently not clear. A systematic investigation quantifying the effect of dispersion (as in Glimsdal et al., 2013) on PTHA for practical source configurations would be desirable.

## Nonlinearity and Resonances (H4)

Most tsunami simulations to date start from an ocean at rest and assume that interaction of currents with gravity waves is

negligible. Androsov et al. (2011) demonstrated that significant alterations of the wave height can be attributed to tidal activity. A quantitative sensitivity analysis of this effect, its dependence on bathymetry, and its correlation to the choice of model (NLSW) is necessary. Huthnance (1975) described the phenomenon of trapped waves on continental shelves that may trigger edge waves and other amplified phenomena. Tsunami resonance effects in Chile and the Balearic Islands are studied in Aranguiz et al. (2019) and Vela et al. (2014). Pattiariatchi and Wijeratne (2015) describe the effect of such phenomena as amplifying factors for meteotsunamis. It is currently unclear how such amplifying phenomena can be represented in the numerical model, nor if the strength is captured adequately.

## Quantifying the Influence of Modeling Assumptions and Scaling (H5)

A hierarchy of modeling approaches, from shallow water assumption, over dispersive long wave solvers, to Navier Stokes type models, can be used to numerically treat tsunami hazard analysis in varying complexity. Due to ever-increasing computational resources, a trend toward more involved model equations can be observed. However, a clear quantitative assessment of the difference has only partly been performed. Lynett et al. (2017) use extensive benchmarking to study and compare modeling approaches to currents induced by tsunami waves. While this study is enlightening and provides very good benchmarking tools, further assessment is necessary to quantify the influence of higher fidelity modeling techniques. Generally, we note that current benchmarking (e.g., Synolakis et al., 2008) stays behind current high-fidelity modeling capabilities. Additionally, some benchmarks based on laboratory experiments have issues with scaling (see Heller, 2011; Pedersen et al., 2013), and related bias and accuracy have not been investigated systematically.

## Modeling Situations With Complex Tsunami Inundation (H6)

NLSW models are predominantly used to simulate tsunami inundation. However, real inundation situations involve features too complex for NLSW approximate modeling, such as urban structures, or damage and erosion due to debris transport. At present, these topics are only partly represented, often using heuristic model formulations. Examples include spatially variable friction mapping (e.g., Gayer et al., 2010; Kaiser et al., 2011), or porous body equivalent friction models representing buildings (e.g., Yamashita et al., 2018). Bottom friction parameterization is almost insensitive for offshore modeling (see Arcos and LeVeque, 2015). However, variable bottom friction parameterizations may pose a viable tool for simulating detailed inundation, but large uncertainties still prevail (e.g., Griffin et al., 2015; Macías et al., 2020). While small scale laboratory tests exist (Park et al., 2013), the heuristic nature of named models and the difficulty to perform controlled tests, implies potentially large epistemic uncertainties. Debris impact and transport are predominantly addressed through post-disaster surveys and experimental analysis of data so far (e.g., Nistor et al., 2017a; Nistor et al.,
2017b; Stolle et al., 2019), and is mostly embedded in only vulnerability analysis (see below), and not in hydrodynamic modeling or PTHA to our knowledge. Extending the modeling dimensions and physical complexity is desirable (e.g., Marras and Mandli, 2021). Open and related to this issue is the influence and potential bias of the accuracy of topo-bathymetric grids, including filtering of structures and vegetation, on the accuracy of inundation simulations (see Griffin et al., 2015; Goda and Song, 2019). Unphysical bias can also be introduced when coupling high resolution (nested) models to large-scale propagation models as shown in Harig et al. (2008).

# PROBABILISTIC TSUNAMI RISK ASSESSMENT

This section discusses identified gaps in PTRA. We go through current state for exposure modeling, physical vulnerability, and risk and resilience metrics, as they naturally follow each other in a consequence-based risk workflow (Figure 1). Methods characterizing the complex social, organizational, and economic context in a tsunami risk assessment are discussed subsequently.

## Gaps in Exposure Modeling
### Existing Methods

Exposure data provide information about the characteristics and location of people and assets at risk. There are several techniques for the acquisition of exposure data, with different degrees of resolution and precision (Pittore et al., 2017). Data from governmental agencies are most commonly used, as they are open and available in most developed countries. These data often provide coverage for the entire building inventory (e.g., physical assets) and are regularly updated for asset management (e.g., national technical maps) and fiscal reasons (e.g., cadastral data). Different exposure databases exist. The Global Exposure Database—GED (De Bono and Mora, 2014; De Bono and Chatenoux, 2015) developed for GAR13 and updated later for GAR15 (UNISDR, 2013; 2015) provides a global dataset at 5 km grid resolution at inland and 1 km at coastal locations, including data for buildings, their use, and exposed value. The 2013 and 2015 versions of the GED served as the exposure databases for the global risk model by the United Nations Office for Disaster Risk Reduction, which considered earthquakes, hurricanes, tsunamis and riverine floods as hazards. The DRMKC Risk Data Hub WebGIS tool (Antofie et al., 2019) has been developed to provide access and sharing of EU-wide information relevant for disaster risk management. Initiatives such as the Open Exposure Data (OED) with roots in proprietary catastrophe modeling and reinsurance industry, provide the opportunity to generate exposure data, including those relevant to tsunami risk, with interoperability between different modeling tools. These databases mainly contain data from census or remote sensing. A recent interview-based approach, relying on local practicing engineers with knowledge of building features, has been adopted for the compilation of building inventories at regional scales (Polese et al., 2020). Careful validation needs to address possible

heterogeneity in data. At present, the only guidelines and tools that exist for capturing and classifying exposure data for a tsunami are the multi-hazard exposure taxonomy, and associated tools, provided by GED4ALL (Silva et al., 2018b).

## Identified Gaps

### Lack of Detail (E1)

Most available exposure data have not been collected for the purpose of tsunami risk assessment and may be missing important information for modeling tsunami fragility or vulnerability. For instance, population cadastral data are often collected at the municipal, district or residential unit level, requiring extra assumptions to determine the geographical distribution. Tsunami hazard intensities can vary considerably between two nearby locations. Accurate geo-localization of the exposed assets and people is needed to obtain robust results, necessitating a minimum resolution level for the exposure databases. While main building construction characteristics are often known, tsunami relevant features (e.g., building lateral load resistance, foundation) are missing (Rivera et al., 2020). Exposure data for critical structures and infrastructure should include functionality information for the exposed asset. This would allow for proper modeling and hence assessment of community resilience, considering different services such as healthcare and education. In other cases, data gaps and uncertainties are associated with regulatory and privacy limitations or outdated sources.

### Lack of Exposure Data (E2)

In many developing countries, where cities have rapid urbanization processes and long-term planning is not consistently enforced, exposure data are not always available or updated. Such data may be inferred from satellite and aerial imagery, from freeware data made available from international projects (e.g., NASA's EOSDIS), from volunteered geographic information systems (e.g., Huyck et al., 2011; Huyck and Eguchi, 2017; OpenStreetMap, 2020), or through intergovernmental organizations (e.g., JRC Risk Data Hub, 2020).

### Lack of Tsunami Exposure Model and Taxonomy (E3)

Significant efforts have been made in the earthquake risk community to define a common exposure taxonomy (e.g., GED4GEM, Silva et al., 2018a; METEOR, Huyck et al., 2019). However, these taxonomies do not contain all the required structural attributes for estimating tsunami risk such as geomorphological, land use, and land cover datasets, or number and size of openings in buildings. A recent development is GED4ALL, a multi-hazard taxonomy (Silva et al., 2018b), which considers tsunami as a hazard. GED4ALL also discusses multiple asset types like buildings, people, infrastructure systems and crops. Common taxonomy and attributes are fundamental to avoid heterogeneity, especially when considering multiple asset types.

### Spatio-Temporal Variability (E4)

Most exposure models are static in time and do not consider the spatio-temporal variability of exposure components. This aspect
is critical when modeling human exposure since there can be daily and seasonal variations. For example, coastal regions often attract tourists, visitors and seasonal workers, leading to significant seasonal fluctuations in the population (Fraser et al., 2014). Spatio-temporal variation in exposure heavily influences the tsunami risk.

# Gaps in Physical Vulnerability

## Existing Methods

As tsunami losses are closely connected to damages to buildings and infrastructure, the vulnerability component is often cut into two parts: a tsunami-to-damage fragility function, and a damage-to-loss consequence function (Figure 1). Advancements in tsunami vulnerability models have significantly lagged behind those of tsunami hazard, with almost no studies found to precede the 2004 Indian Ocean Tsunami (Charvet et al., 2017). However, with the recent devastating tsunamis providing a large quantity of observed damage and loss data to develop and validate fragility and vulnerability models, this field of study has rapidly grown. Several empirical fragility functions for the assessment of buildings (Koshimura et al., 2009; Mas et al., 2012; Suppasri et al., 2014; Charvet et al., 2015; Chock et al., 2016) and infrastructure (Eguchi et al., 2014; Hatayama, 2014) have been derived from observed damage in the 2004 Indian Ocean, 2009 Samoa, 2010 Chile, and 2011 Tohoku tsunamis. Recently, analytical fragility functions were derived from numerical simulations of building response under tsunami inundation (Petrone et al., 2017; Alam et al., 2018; Karafagka et al., 2018; Páez-Ramírez et al., 2020), and under sequential earthquake and tsunami impact (Park et al., 2012; Attary et al., 2019; Petrone et al., 2020). Only a few studies exist that move from fragility to vulnerability modeling (De Risi et al., 2017). There is a lack of consensus on many aspects of physical fragility and vulnerability modeling.

## Identified Gaps

### Limitation in Asset Types and Geographical Scope (P1)

The vast majority of existing tsunami fragility and vulnerability models relate to buildings, few exist for bridges, fuel tanks, or other types of infrastructure. The main reason is that most fragility functions are empirical, and few observational damage or loss data are available for infrastructure components. Even for buildings, the geographical scope of existing vulnerability and fragility models is limited. Most empirical fragility functions are based on data from the 2004 Indian Ocean and 2011 Tohoku events, and hence represent non-engineered buildings in countries surrounding the Indian Ocean and engineered buildings typical of Japan. With analytical fragility functions only covering a small number of building types, large portions of the world's exposure remain unrepresented by current studies.

### Effect of Multiple Hazard on Empirical Tsunami Fragility Mode (P2)

Tsunamis are commonly triggered by large earthquakes. Near-source, observational data on asset damage and loss collected after the tsunami often include the combined effects of earthquake ground shaking and tsunami inundation. Hence,

empirical fragility and vulnerability models derived from such data inherently comprise the effects of both hazards. Therefore, corresponding empirical fragility models may be regarded as inappropriate for use in a tsunami-only risk assessment. Pure tsunami damage data is rare and currently limited to non-engineered structures (Charvet et al., 2017).

### Lack of Consensus Regarding Best Tsunami Intensity Measure (P3)

The intensity measure *IM* (Figure 1) links the hazard and vulnerability components within risk models. Traditionally, tsunami inundation maps are presented in terms of inundation depth. While the majority of fragility and vulnerability models adopt inundation depth as *IM*, other tsunami *IM* have also been used such as the flow velocity or momentum flux. The absence of inundation velocity measurements in field data requires running tsunami inundation simulations to use such *IM* (Koshimura et al., 2009; Song et al., 2017). More recently, force-based *IM* (e.g., flow velocity, momentum flux) were used in fragility functions for engineered buildings yielding better correlation to observed damage than inundation depth (Macabuag et al., 2016). However, no consensus on the most appropriate *IM* could be reached. As a consequence, mismatches between representations of hazard and vulnerability in risk modeling may exist.

### Gaps in Building Analysis and Assessment for Use in Analytical Tsunami Fragility (P4)

Buildings are often used as vertical evacuation shelters and an assessment of their structural fragility is therefore an important information in the risk assessment workflow. Tsunami engineering being a younger discipline than earthquake engineering has adopted approaches from the latter community. This was supported by the physical similarity of both hazards applying predominantly horizontal loads to structures. However, there are fundamental differences in how earthquake and tsunami loads are applied to buildings. For example, tsunami loads affect the lower floors of a high-rise building, whereas seismic loads are inertial forces usually causing increasing magnitude for higher floors (Baiguera et al., 2019). Thus, earthquakes induce large bending moments in structural elements, whereas tsunamis typically induce large shear. Since typical structural modeling approaches tend to prioritize flexural effects, the bias in tsunami fragility assessment may be large. Furthermore, seismic loads are dynamic, whereas loads from tsunami inundation can be considered quasi-static, and Rosetto et al. (2018) have shown that building ductility is often not crucial in the tsunami response of structures. Although no consensus has been reached in this regard, more fragility functions based on static rather than time-dependent non-linear approaches are derived now (Petrone et al., 2017; Rosetto et al., 2019). As a tsunami applies direct pressures to a structure, non-structural components like infill walls (and their openings) are seen to play an important role in determining tsunami forces (Del Zoppo et al., 2021). Furthermore, buoyancy, foundation scour and debris impact, which significantly affect building damage from tsunami inundation are rarely modeled (Del Zoppo et al., 2019). These
effects are still to be investigated; therefore, published analytical tsunami fragility functions are subjected to large modeling uncertainties. Progress towards more comprehensive and reliable analytical fragility and vulnerability models is needed.

# Gaps in Risk and Resilience Metrics
## Existing Methods

Tsunami risk assessments typically reflect the impact on the exposed population and infrastructure. The most commonly used decision variables (or metrics) are the number of fatalities, injuries, affected people, besides the direct and indirect economic losses. Direct economic losses represent the repairing/replacement costs of damaged assets, whereas indirect losses reflect costs as down-time, partial loss of functionality of buildings and infrastructure, loss or reduction in network connectivity, flow and/or capacity. These metrics can be used in alternative approaches such as worst-case scenarios, scenario-based for a prescribed return-period, and fully probabilistic. A review of early methods for tsunami risk assessment can be found in Jelínek and Krausmann (2008).

Fully probabilistic risk assessments require the integration of hazard estimates (PTHA) with vulnerability functions (see **Figure 1**, Løvholt et al., 2015a; 2019). Since the results of PTHA are not always available, tsunami risk assessments are often performed considering selected (worst-case) scenarios as hazard input (e.g., Triantafyllou et al., 2019), which sometimes represent past disasters (e.g., Daniell et al., 2017). Having the results of PTHA available, tsunami risk assessment can also be performed for a limited number of scenarios (e.g., Nadim and Glade, 2006; Okumura et al., 2017). When the PTHA results are available in the form of stochastic event sets, a fully probabilistic tsunami risk assessment (PTRA) can be performed (Ordaz, 2000; Strunz et al., 2011; Salgado-Gálvez et al., 2014), although these types of analyses usually demand an extensive computational effort (e.g., Løvholt et al., 2015a; Jaimes et al., 2016; Goda and Song, 2019; Ordaz et al., 2019).

In a fully probabilistic tsunami risk assessment workflow, risk results are obtained in terms of exceedance frequencies for the above-mentioned metrics (**Figure 1**). For instance, loss exceedance curves (LEC) provide the relationship between loss values and their exceedance frequencies (Løvholt et al., 2015a; Jaimes et al., 2016; Attary et al., 2017; Ordaz et al., 2019). The area under the LEC corresponds to the average annual loss (AAL), a metric that provides a long-term overview of risk and accounts for the contribution of large and infrequent events as well as small and more frequent ones. From the LEC, loss values associated with a given return period can be obtained, such as loss values estimated by Løvholt et al. (2015a) at a global level representing direct losses. The Hazus tsunami loss estimation methodology provides state-of-the-art decision-support software for estimating potential losses from tsunami events (FEMA, 2017a; FEMA, 2017b).

Risk assessment is not necessarily limited to quantifying the direct and indirect impact on exposed populations and infrastructures. The evaluation of safety and reliability of physical systems is of interest too and for this, fragility functions ("Gaps in Physical Vulnerability" section) can be integrated with hazard to obtain the frequency of exceeding a

given damage level (see **Figure 1**, e.g., Park et al., 2019; Fukutani et al., 2019). The risk metrics provide valuable data also for the assessment of quantitative resilience (also denoted as engineering resilience), which aims to estimate the resilience of a network, an infrastructure, or even an urban ecosystem to a specific natural hazard (see Mebarki et al., 2016 for industrial plants, Akiyama et al., 2020 for bridges). Quantitative resilience should not be confused with coastal community resilience which is discussed in detail in the following section.

## Identified Gaps

### Gaps Related to Characterization and Propagation of Uncertainties (R1)

Most existing PTRA models rely on a homogeneous Poisson process as the probabilistic backbone for the occurrence process (Eq. 4). The Poisson model, strictly speaking, should be used for propagating only those uncertainty sources that renew with the occurrence of each new event (Kiureghian, 2005). This means that propagation of other sources of uncertainties in a PTRA framework (i.e., those that lack renewal properties), such as the uncertainties in modeling, analysis method, and in general epistemic uncertainties, need more research (Goda and De Risi, 2018; Goda, 2020). One possible direction could point to Bayesian methods (Jalayer and Ebrahimian, 2020).

### Challenges in Characterizing Vulnerability Functions (R2)

PTRA lacks a clear distinction and definition of the different loss components that are quantified through the vulnerability functions. On the one hand, direct economic losses can be estimated with a good degree of confidence using existing methodologies (Pagnoni et al., 2019). Long-term direct (e.g., cost of maintenance) and indirect losses (e.g., down-time and reduced functionality including business interruption) typically represent a significant component of the total economic loss (direct + indirect) yet require better quantitative approaches.

### Lack of a Tsunami Consequences Database (R3)

There is a lack of tsunami-specific consequence databases accounting for casualties and losses (Yamao et al., 2015). These types of databases exist for disasters in general (e.g., EM-DAT) and more specifically for earthquakes (So et al., 2012; Cardona et al., 2018). They are useful not only to keep a consistent record of past events and the affected regions but to disaggregate the impacts of large events in terms of losses (direct and indirect) and casualties (fatalities and injured), besides assessing the consequences in particular sectors (e.g., road networks, heritage sites, etc.) at different resolution levels. The information included in the consequences databases provides valuable data to validate and calibrate different components of the models (e.g., fragility curves, vulnerability functions). Some data can be partially acquired from collections of documented eyewitness accounts (Santos and Koshimura, 2015), or other sources (e.g., ITIC, 2020).

### General Lack of Risk Studies for Networks and Lifelines (R4)

Current implementations of PTRA are mainly focused on residential buildings and emergency planning activities such as
the definition of evacuation routes. However, the resilience of coastal areas relies on conventional and strategic infrastructures (Akiyama et al., 2013; Pitilakis et al., 2019). Conventional infrastructure such as roads, bridges, power, water, sanitation and communication networks, underpin economic and social activities in most urban areas (Salgado-Gálvez et al., 2018). Schools and hospitals support the provision of education and health services, which are essential to recovery. Critical infrastructures in coastal areas include harbors (nuclear) power plants, gas and oil storage, and early warning infrastructure, such as tidal buoys and offshore bottom pressure gauges (De Risi et al., 2018). Such infrastructures are complex, often interconnected and geographically distributed systems involving multiple sectors (Duenas-Osorio and Vemuru, 2009; Argyroudis et al., 2019), where further research is needed to quantify their resilience to tsunamis.

## Assessing Tsunami Risk in a Multi-Hazard and Multi-Risk Framework (R5)

As triggered events, tsunamis fit naturally within a multi-hazard framework. Moreover, there can be several cascading consequences associated with the occurrence of tsunamis, such as technological disasters induced by natural hazards known as NATECH risks (e.g., the Fukushima Disaster), disruption to supply chains, and societal impacts. Therefore, management and decision-making for tsunami risk should be framed in a multi-risk context. To be able to make risk-informed decisions considering tsunamis, it is important to model the interaction of tsunamis with other phenomena at the level of hazards, vulnerabilities, and socio-economic consequences. An important gap related to risk assessment for tsunamis (and in general) is the lack of a streamlined and standard workflow for modeling the multi-hazard and multi-risk aspects. Currently, most studies consider the different hazards to be independent or “simultaneous” (e.g., earthquake and tsunami as independent events); whereas, few works consider interacting hazards such as coupled simulation of tsunami and earthquake (De Risi and Goda, 2016; Goda et al., 2017; Goda and De Risi, 2018; Ordaz et al., 2019; Park et al., 2019), the cumulation of tsunami and earthquake damages and losses (Ordaz, 2015; Attary et al., 2019; Park et al., 2019; Petrone et al., 2020), and interaction of tsunami and aging infrastructure (Akiyama et al., 2020).

## Lack of Understanding and Quantification of Mortality (R6)

Strikingly, the 2004 tsunami with more than 226,000 dead and missing people (EM-DAT, 2020) caused an order of magnitude higher fatalities than the 2011 Tohoku tsunami with 19,846 (EM-DAT, 2020). Hence, past major disasters indicate that the vulnerability to tsunami mortality of a population is much more sensitive to demographic factors (Løvholt et al., 2014) than to physical vulnerabilities (“Gaps in Physical Vulnerability” section). Correlations of tsunami flow depth and number of fatalities following the 2004 Indian Ocean, 2006 Java and 2011 Tohoku tsunamis reveal much larger scatter than those observed in physical vulnerability functions, even when derived from the same events (Reese et al., 2007; Koshimura et al., 2009; Suppasri et al., 2016). As human behavior influences mortality strongly

(Johnston et al., 2016; Blake et al., 2018), deriving simplified vulnerability charts based on single tsunami intensity measures may not be appropriate. Tsunami awareness and availability of tsunami early warning systems and infrastructure are important (Gregg et al., 2006; Fraser et al., 2014), as well as proximity to source areas. Our understanding and ability to quantify and assess the effect of all these factors on tsunami mortality is still very limited.

## The Weakness of Capturing Multi-Faceted Aspects of Vulnerability (R7)

Quantitative risk assessments typically address several socio-economic parameters (e.g., safety, downtime, direct and indirect economic losses, and even human behavior and response) as dimensions of consequences to disruptive tsunami events. However, PTRA falls short in modeling some dimensions of vulnerability that are part of a given context and not directly caused by a disruptive event (e.g., governance-related issues, adaptation and coping capacities, societal inequalities). There are no established methodologies, within the context of the PTRA framework (Equations 1–4), for characterizing context-based impacts of tsunami on the social, political and economic dimensions, leaving it unclear how to address these dimensions. Integrated and heuristic approaches such as "MOVE" (Birkmann et al., 2013) or holistic approaches as those proposed by Carreño et al. (2007) or Aguirre-Ayerbe et al. (2018), have strived to address the context-based and multi-dimensional nature of vulnerability and risk and could be adapted to be used as physical risk indicators in the outcomes of PTRA.

# Gaps in Social Vulnerability, Multi-Dimensional Vulnerability and Risk Indicators

## Existing Methods

Although not directly addressing tsunami risk, Jasanoff (1993) pointed out the urge to bridge the two cultures of quantitative and qualitative risk assessment, stressing the importance to view risk in a larger context of social justice (*who* should we protect, from *which* harm, at *what* cost, and by foregoing *what* other opportunities). The societal factors impacting vulnerability and risk are mainly rooted in a complex and diverse aggregate, which varies over time and space. Qualitative vulnerability investigations use models and frameworks considering several dimensions (e.g., economic, demographic, psychological, political or physical), summarized by composite vulnerability and risk indices. These indicators can be distinguished from the risk and resilience metrics discussed in the previous section (“Gaps in Risk and Resilience Metrics” section) since some of them cannot be directly integrated into a computational PTRA procedure. Examples of existing multi-dimensional vulnerability and risk indicators are: The *community resilience* (e.g., Lam et al., 2016; Saja et al., 2019), the *urban disaster risk index* (Carreño et al., 2007; Salgado-Gálvez et al., 2016), the *social vulnerability index* (Cutter et al., 2003; Flanagan et al., 2011), the *Coastal vulnerability index* (McLaughlin and Cooper, 2010), the *Metropolitan Tsunami Human Vulnerability Assessment* (Tufekci et al., 2018).
## Identified Gaps

### The Difficulty of “Quantifying” Social Vulnerability (II)

Social vulnerability describes combinations of social, cultural, economic, political, and institutional processes that determine differentials in the experience of hazards and recovery from dangerous events (Spielman et al., 2020). Experts may construct meaningful indicators to include a social component in hazard planning, preparation, and response. Integrating social vulnerability research into emergency and disaster risk management is essential, but caution is required to assign quantitative elements. Integration of social factors may allow planners and decision-makers to better identify problems in case of destructive events and provide insights into addressing recovery solutions (Cardona, 2001; Chakraborty et al., 2005; Schmidtlein et al., 2008). Social Vulnerability Index (SoVI) is a single quantitative indicator which was developed through a review of hazard case studies by Cutter et al. (2003) examining the spatial patterns of social vulnerability to natural hazards at the county level in the United States. Because of the complex and multidimensional nature of factors contributing to vulnerability, no variable has yet been identified to fully validate SoVI. An alternative approach to assess its reliability is to identify how the changes in the SoVI algorithm construction may lead to the changes in the outcome. Schmidtlein et al. (2008) investigated the sensitivity of quantitative features of the SoVI such as the scale of application, the set of used variables, and various geographic contexts.

### Ambiguities in Definition of Community Resilience (I2)

Resilience is a frequently used term to discuss the capacity of a society or ecosystem to recover quickly from a disaster. The United Nations Office for Disaster Risk Reduction has defined resilience as “the capacity of a system, community or society potentially exposed to hazards to adapt, by resisting or changing in order to reach and maintain an acceptable level of functioning and structure. This is determined by the degree to which the social system is capable of organizing itself to increase this capacity for learning from past disasters for better future protection and to improve risk reduction measures” (UNISDR, 2007). A comprehensive review of various definitions of resilience can be found in Davoudi et al. (2012) and Ayyub (2014). The definition of coastal resilience is hindered by varying definitions and non-unified terminology, difficulties in selecting and combining different resilience indicators, and lack of data for validation (Lam et al., 2016). In fact, resilience is still lacking rigorous measurement methods (Bozza et al., 2015), especially in the context of tsunami hazard (Genadt, 2019).

### Lack of Tsunami Vulnerability Index (I3)

A specific tsunami Disaster Risk Index (TDRI), similarly to the Disaster Risk Index (DRI) developed by the UN Development Program to compare disaster risk between countries exposed to hazards (UNDP, 2004) or the Urban Seismic Risk Index by Carreño et al. (2007) should be developed.

### Integrated Approaches to Consider the Multi-Dimensional Aspects of Tsunami Risk (I4)

Vulnerability and risk are multi-faceted concepts and encompass various assets, physical, organizational, and institutional

dimensions (e.g., Herslund et al., 2016). Vulnerability and risk assessment considering these different facets often requires different scientific backgrounds and approaches (Hufschmidt et al., 2005). A consequence-based approach to risk assessment (e.g., the PEER framework, or computational PTRA in general) has its roots in engineering. The approach follows a logical flow from causes associated with a disruptive event toward quantifying its direct and indirect socio-economic consequences. This approach focuses on the physical dimension of vulnerability, acting as a “container” of functions and services and thereby invokes–directly or indirectly–other dimensions of vulnerability such as social, economic and organizational vulnerability. On the other hand, the context-based approach (e.g., approaches based on integrated indicators) has its roots in the humanities and social disciplines. This approach deals with the context and the interactions between different actors, the respective territory, the different drivers (climate, societal changes) and how decisions can affect the overall context and the complex interplay between actors and drivers. Needless to say, the two approaches complement each other and have to be taken into account in policymaking for DRR in an integrated manner (O'Brien et al., 2007).

### Considering Community Response and Organizational Capacities (I5)

Recent tsunami events worldwide have highlighted the need to critically revisit how human behavior in tsunami evacuation, and more generally, the human dimension of preparedness for tsunamis is addressed within the risk assessments. The lessons from Japan 2011, Chile 2010 and Indonesia in 2010 and 2018 events highlight such needs. Questions arise on how and if the different and seemingly inconsistent human behavior can be addressed in tsunami risk assessments. Moreover, atypical events such as the Krakatoa, Indonesia 2018, do not allow for conventional prevention, warning and mitigation strategies. In most cases, aid and help arrive late due to limited organizational capacities, leaving the affected communities in even more vulnerable conditions, especially during the first critical hours and days after the event. Events with growing levels of complexity are likely to continue to occur in the future and this calls for a more in-depth consideration of how different communities respond and how those variations can be integrated within the risk assessment framework.

### Incorporating Risk Perception in the Formulation and Analysis of Complex Risks (I6)

Perceptions are dynamic and socially constructed. Perceptions can change abruptly or gradually, depending on the context. Understanding evacuation behavior requires an understanding of risk perceptions. This can help explain why the response to tsunami drills may be different than when responding to a real event. It is quite challenging for risk methodologies to consider the dynamic, complex and subjective aspects of risk perception. Only by understanding the subjective meanings of perceived risks allows risk communication to be designed and applied more effectively.
# CONCLUSIONS AND DIRECTIONS

In this review, we discuss a large number of research gaps in PTHA and PTRA. It becomes obvious that methods have substantially improved over the past decades, but also that open questions remain in the physical description, conceptualization, modeling, as well as the social and psychological dimensions of the topic.

The physics and geological complexity of tsunamigenic sources are still not captured nor understood adequately, leading to large uncertainties. For SPTHA, neither all earthquake faults nor their exact location, geometry, boundary and initial conditions (e.g., stress, friction) are fully constrained. Statistical models of recurrence constitute the largest uncertainties in large and rare events, including tsunami earthquakes. Uncertainty may become excessive for landslide tsunamis, where statistics on past events often are absent, and our understanding of slope failure probability is limited. The need for covering vast geographical scales, source diversity and related uncertainties, render LPTHA extremely challenging. For VPTHA additional difficulties arise due to the complexity of tsunamigenic volcano sources and triggers, but they are constrained spatially. MPTHA may benefit from a large meteorological data network allowing for (prototypical) forecasting as well as PTHA applications, but sensitivity to source parameters is still unconstrained.

While modeling and parameterization of individual phenomena are possible, they are often excessively computationally expensive or highly uncertain due to missing constraints on input parameters. The multiple scales involved in PTHA from far-field propagation over oceanic distances to the need to resolve small scale inundation features while capturing physics and resolving uncertainties still represent an open challenge. Yet, this solution is needed to convey PTHA information properly into risk analysis.

Even more challenging is the situation in PTRA, where gaps exist in the transformation of physical hazard to risk and quantifying the uncertainties in the assessment of risk and resilience. Key concepts, such as physical vulnerability and mortality and their related uncertainties, are less developed than the main PTHA elements. There are gaps regarding selection of IM, limited observed damage asset- and location-wise, and limited experimental validation.

Furthermore, tsunami science is immature concerning embedding issues with intrinsically multi-hazard and multi-risk aspects, such as the cascading events that are entangled with tsunami hazards. A weakly developed link between quantitative PTRA and the social sciences is a clear gap. At this point, it is worth noting that terms “vulnerability” and “resilience” are multi-dimensional concepts that are used both in the consequence-based-natural sciences inspired-as well as context-based approaches-motivated by social sciences. Therefore, they may have quite different interpretations depending on the analysis context.

The overarching issue is integrating all the above components and developing an overall consistent sensitivity and uncertainty quantification framework, to understand tsunami risk and

identify risk drivers, from the probability of the sources causing hazards to the probability of their physical consequences and societal impact. This understanding must be developed and prioritized in future research.

To guide such efforts, we have performed an expert judgment exercise that we discuss in the following subsection. It may help to identify most pressing research needs as well as prioritize research efforts.

## Prioritizing Research Gaps

A scientific sensitivity analysis of the impact of each research gap, as conducted for individual sources in Sepúlveda et al. (2017) or Davies and Griffin (2020), on the overall result of a PTHA or PTRA is certainly out of the scope of a single review paper. However, some guidance on prioritization of efforts is certainly desirable. Since we focused our description on research gaps, we suggest two important metrics for the prioritization: The susceptibility of PTHA and PTRA results on uncertainty due to the research gap (sensitivity) and the difficulty or amount of research effort needed to fill that respective gap (tractability).

In order to assess these two metrics, we conducted a first-pass expert judgment among the more than 50 co-authors of this article-all experts in one or more of the aspects of our review. A questionnaire was designed that asked three questions for each of the 47 research gap subsections that we have described before. The first two questions addressed the two metrics just mentioned. The third question asked if experts were of the opinion if the research gap existed because of a missing theoretical understanding, a lack of data, or both. While this somewhat ad hoc prioritization is not as solid as a rigorous expert elicitation (e.g., Cooke, 1991; Budnitz et al., 1997; Morgan, 2014; for tsunami hazard see an application in Basili et al., 2021, or the discussion in Grezio et al., 2017) and hence could be somehow biased, we believe it still provides a valuable starting point for future efforts. It is a qualitative broad-brush answer to the question, which research gap may be of highest importance. More details on this exercise are given in the **Supplementary Material**.

The result of this exercise is visualized in a priority matrix (Figure 2). It may appear natural to respond first to those research gaps that are located in the left upper quadrant of the matrix, since these gaps are considered less difficult to solve, while they are expected to influence the risk considerably. It can be noted that most of the research gaps are judged hard to solve but with a highly sensitive impact on the overall result. This seems natural, since high impact but simple problems would have been solved already.

Based on our qualitative assessment, we can therefore identify some overall trends. First, we see some common challenges related to establishing annualized source probability of occurrence, which tend to cluster in the upper right corner of Figure 2. This means that they are considered relatively most important, yet hardest to solve. Of these, obtaining landslide related annual source probabilities (L1) is considered both the largest yet most important obstacle, while a just slightly lower similar prioritization is evident for earthquake and volcano sources (S1 and V2). Another aspect that is considered important (and challenging) is the multi-hazard and cascading
<figure><img src="image_5.png" /></figure>

hazard aspect (R5). On the other hand, the research gaps that appear to be least sensitive and also easy to be filled are related to the numerical modeling of wave propagation (H3), as well as lack of joint intensity measures (I3) and gaps related to earthquake scaling relations (S4). Finally, we also note **Figure 2** allows us to analyze several instances of components with similar sensitivity but with clearly different tractability. For instance, the lack of tsunami exposure data (E2) is considered as important as modeling complicated aspects of inundation (H6), but the former is assumed by the authors of this paper to be more easily achieved. Several other similar examples can be analyzed from **Figure 2**.

It is noteworthy that most of the research gaps that most experts find consensus on are highly sensitive in their impact (all located at the upper margin of the point cloud). It is also worth

noting that most research gaps are considered to relate to data and theory gaps and that those that relate to only a missing theoretical understanding are considered of relatively low sensitivity. This may be related to the fact that when we don't understand a phenomenon, we cannot really judge whether it affects our results or not. In other words, this may be an "unknown". Whereas a data related research gap may already have proved to be sensitively influential by a specific example, but due to a lack of data cannot be involved concisely into the workflow.

This priority matrix is just a very first approach. Since tsunami research eventually aims at protecting life from natural hazard, one could also prioritize those research gaps with direct impact on this goal. These would be in particular those topics mentioned in sections "Gaps in Physical Vulnerability," "Gaps in Risk and
*Resilience Metrics,*" and “Gaps in Social Vulnerability, Multi-Dimensional Vulnerability and Risk Indicators" (marked with P, R, and I; respectively).

Final Considerations

We have described and prioritized a comprehensive list of
research gaps in PTHA and PTRA. While our approach to
prioritization and the metric used to do so are to some extent
subjective, it remains for the scientific community and further
investigation as well as future incentives to decide, which
directions to choose from. Nevertheless, our priority matrix
will serve as a first impression on the weight of each of the
identified research gaps.

An important part of the future puzzle will be exploring how
uncertainties propagate to risk across disciplines. While
uncertainties are more extensively explored in earthquake-
related hazard analysis, non-seismic hazard, vulnerability,
exposure and risk are lagging behind. On the other hand,
different levels of maturity of methods and understanding will
always exist. Hence, it is imperative to develop PTRA standards
and guidelines to appropriately merge all risk analysis
components considering their different uncertainty exploration
and maturity level.

While validation of individual components has been addressed
in several of the sections in our text, validating the PTHA and
PTRA workflow as a whole is still ongoing research. Marzocchi
and Jordan (2014) propose a methodology for a meaningful test
of general probabilistic hazard models and an example of a
successful application can be found in Meletti et al. (2021).

Certainly, research gaps exist also outside of the scope of PTHA and PTRA. New computational methods, like fuzzy methods, machine learning techniques and even advances in classical computational methods have to be considered. Rigorous, information theory inspired approaches to validation may also be explored.

Considering the goals of the Sendai Framework for Disaster Risk Reduction and acknowledging the vast number of challenges outlined in the sections before, a concerted interdisciplinary effort to close the most pressing gaps is required. Attempts to gather expertise, facilitate exchange and development, and coordinate community efforts are represented by the Global Tsunami Model (GTM, 2020) and the COST Action AGITHAR. A thorough consolidation of available sources of information in openly accessible databases, documentation of standard workflows, unification of terminology and metrics, as well as information hubs need to be established.

REFERENCES

Abadie, S., Paris, A., Ata, R., Le Roy, S., Arnaud, G., Poupardin, A., et al. (2020). La Palma landslide tsunami: calibrated wave source and assessment of impact on French territories. *Nat. Hazards Earth Syst. Sci.* 20, 1–20. doi:10.5194/nhess-20-1-202010.5194/nhess-20-3019-2020

AGITHAR. (2020). Accelerating global science in tsunami hazard and risk analysis. Available at: https://www.agithar.uni-hamburg.de (Accessed November 5, 2020).

AUTHOR CONTRIBUTIONS

All Authors contributed in early stages of the manuscript by
individual contributions from their respective research field.
Major contributions are listed as follows: Early compilation of
text: MS, FL, SL, JB, IA-Q, TR; Abstract and Conclusions: JB, FL,
SL, JS; Introduction: FL, FJ, SL, JB, JS, MS-G; Probabilistic
Framework: FJ, MS-G, IA-A; Earthquake Sources: SL, JS, FL,
KJ, MV, SM; Landslide Sources: FL, RP, SA, AS; Volcano Sources:
JS, FL, SA; Meteotsunami Sources: JB; Hydrodynamic Modeling:
JB, FL, JM, IA-Q, AB; Exposure: FJ, TR, MS-G; Physical
Vulnerability: TR, FJ, JS, MDZ; Risk/Resilience: IA-A, FJ;
Social Vulnerability: FJ, IA-A, IR; substantial revisions: AG,
KJ, SL, SM, RDR, MS, JS, RP, AB, RB, SA, MV, MDZ, AS, IA-
Q, CC; internal review: all authors.

FUNDING

This article is based upon work from COST Action CA18109 AGITHAR, supported by COST (European Cooperation in Science and Technology). VB and PC obtained support through the VES20 Inter-Cost LTC 20020 project. MS-G obtained support through the Severo Ochoa Centers of Excellence Program (Ref. CEX 2018-000797-S). TU acknowledges funding from the European Union's Horizon 2020 research and innovation program (ChEESE project, Grant Agreement No. 823844).

ACKNOWLEDGMENTS

This work would not have been possible without the help of many
researchers, who are too many to be named individually, but we
thank all participants of COST Action AGITHAR for their
valuable input in official or informal conversations, workshops
and interactions. We would like to thank the reviewers and editor
for very considerate and challenging suggestions that improved
this manuscript substantially.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at:
https://www.frontiersin.org/articles/10.3389/feart.2021.628772/
full#supplementary-material

Aguirre-Ayerbe, I., Martínez Sánchez, J., Aniel-Quiroga, Í., González-Riancho, P., Merino, M., Al-Yahyai, S., et al. (2018). From tsunami risk assessment to disaster risk reduction—the case of Oman. *Nat. Hazards Earth Syst. Sci.* 18, 2241–2260. doi:10.5194/nhess-18-2241-2018

Akiyama, M., Frangopol, D. M., Arai, M., and Koshimura, S. (2013). Reliability of bridges under tsunami hazards: emphasis on the 2011 Tohoku-oki earthquake. *Earthq. Spectra* 29 (S1), S295–S314. doi:10.1193/1.4000112

Akiyama, M., Frangopol, D. M., and Ishibashi, H. (2020). Toward life-cycle reliability-, risk- and resilience-based design and assessment of bridges and bridge networks under independent and interacting hazards: emphasis on
earthquake, tsunami and corrosion. *Struct. Infrastruct. Eng.* 16 (1), 26–50.
doi:10.1080/15732479.2019.1604770

Alam, M. S., Barbosa, A. R., Scott, M. H., Cox, D. T., and van de Lindt, J. W. (2018).
Development of physics-based tsunami fragility functions considering
structural member failures. *J. Struct. Eng.* 144 (3), 04017221. doi:10.1061/
(ASCE)ST.1943-541X.0001953

Albini, P., Musson, R. M. W., Rovida, A., Locati, M., Gomez Capera, A. A., and
Viganò, D. (2014). The global earthquake history. *Earthq. Spectra* 30 (2), 607.
doi:10.1193/122013EQS297

Allen, T. I., and Hayes, G. P. (2017). Alternative rupture-scaling relationships for subduction interface and other offshore environments. *Bull. Seismol. Soc. Amer.* 107, 1240–1253. doi:10.1785/0120160255

Androsov, A., Behrens, J., and Danilov, S. (2011). "Tsunami modelling with unstructured grids Interaction between tides and tsunami waves," in, *Computational science and high performance computing IV*. Editors E. Krause, Y. Shokin, M. Resch, D. Kröner, and N. Shokina (Heidelberg: Springer Berlin), 115, 191–206. doi:10.1007/978-3-642-17770-5_15

Annaka, T., Satake, K., Sakakiyama, T., Yanagisawa, K., and Shuto, N. (2007).
Logic-tree approach for probabilistic tsunami hazard analysis and its
applications to the Japanese coasts. *Pure Appl. Geophys.* 164, 577–592.
doi:10.1007/s00024-006-0174-3

Antofie, T., Luoni, S., Faiella, A., and Marin Ferrer, M. (2019). *Risk data hub—web platform to facilitate management of disaster risks*. Luxembourg: Publications Office of the European Union. doi:10.2760/68372

Araki, E., Saffer, D. M., Kopf, A. J., Wallace, L. M., Kimura, T., Machida, Y., et al.
(2017). Recurring and triggered slow-slip events near the trench at the Nankai
Troaugh subduction megathrust. *Science* 356 (6343), 1157–1160. doi:10.1126/
science.aan3120

Aranguiz, R., Catalán, P. A., Cecioni, C., Bellotti, G., Henriquez, P., and González, J.
(2019). Tsunami resonance and spatial pattern of natural oscillation modes
with multiple resonators. *J. Geophys. Res. Oceans* 124 (11), 7797–7816. doi:10.
1029/2019JC015206

Arcos, M. E. M., and LeVeque, R. J. (2015). Validating velocities in the GeoClaw tsunami model using observations near Hawaii from the 2011 Tohoku tsunami. *Pure Appl. Geophys.* 172, 849–868. doi:10.1007/s00024-014-0980-y

Argyroudis, S. A., MitoulisA., S., Winter, M. G., and Kaynia, A. M. (2019). Fragility of transport assets exposed to multiple hazards: state-of-the-art review toward infrastructural resilience. *Reliab. Eng. Syst. Saf.* 191. 106567. doi:10.1016/j.ress.2019.106567

Attary, N., Unnikrishnan, V. U., van de Lindt, J. W., Cox, D. T., and Barbosa, A. R.
(2017). Performance-based tsunami engineering methodology for risk assessment of
structures. *Eng. Struct.* 141, 676–686. doi:10.1016/j.engstruct.2017.03.071

Attary, N., Van De Lindt, J. W., Barbosa, A. R., Cox, D. T., and UnnikrisAhnan, V. U. (2019). Performance-based tsunami engineering for risk assessment of structures subjected to multi-hazards: tsunami following earthquake. *J. Earthq. Eng.* 7, 1–20. doi:10.1080/13632469.2019.1616335

Auker, M. R., Sparks, R. S. J., Siebert, L., Crosweller, H. S., and Ewert, J. (2013). A statistical analysis of the global historical volcanic fatalities record. *J. Appl. Volcanol.* 2, 2. doi:10.1186/2191-5040-2-2

Ayyub, B. M. (2014). Systems resilience for multihazard environments: definition, metrics, and valuation for decision making. *Risk Anal.* 34, 340–355. doi:10.1111/risa.12093

Baiguera, M., Rossetto, T., Robertson, I. N., and Petrone, C. (2019). "Towards a tsunami nonlinear static analysis procedure for the ASCE 7 standard," in 2nd International conference on natural hazards and infrastructure. Available at: https://discovery.ucl.ac.uk/id/eprint/10078809

Barbato, M., Petrini, F., Unnikrishnan, V. U., and Ciampoli, M. (2013).
Performance-based hurricane engineering (PBHE) framework. *Struct. Saf.*
45, 24–35. doi:10.1016/j.strusafe.2013.07.002

Basili, R., Brizuela, B., Herrero, A., Iqbal, S., Lorito, S., Maesano, F. E., et al. (2021).
The making of the NEAM tsunami hazard model 2018 (NEAMTHM18). *Front. Earth Sci.* 8, 616594. doi:10.3389/feart.2020.616594

Basili, R., Tiberti, M. M., Kastelic, V., Romano, F., Piatanesi, A., Selva, J., et al. (2013). Integrating geologic fault data into tsunami hazard studies. *Nat. Hazards Earth Syst. Sci.* 13, 1025–1050. doi:10.5194/nhess-13-1025-2013

Beauval, C., Yepes, H., Audin, L., Alvarado, A., Nocquet, J.-M., Monelli, D., et al. (2014). Probabilistic seismic-hazard assessment in quito, estimates and uncertainties. *Seismol. Res. Lett.* 85 (6), 1316–1327. doi:10.1785/0220140036

Bebbington, M. (2010). Trends and clustering in the camps of volcanic eruptions. *J Geophys. Res.* 115, B01203. doi:10.1029/2009JB006581

Bebbington, M. (2008). Incorporating the eruptive history in a stochastic model for volcanic eruptions. *J. Volcanol. Geothermal Res.* 175, 325–333. doi:10.1016/j.jvolgeores.2008.03.013

Begét, J. E., Cecioni, C., and Yepes, H. (2000). "Volcanic tsunamis," in *Encyclopedia of volcanoes*. Editor Sigurdsson, H. (London, UK: Academic Press), 1005–1013.

Behrens, J., and Dias, F. (2015). New computational methods in tsunami science. *Phil. Trans. R. Soc. A.* 373 (2053), 20140382. doi:10.1098/rsta.2014.0382

Bellotti, G., Cecioni, C., and de Girolamo, P. (2008). Simulation of small-amplitude frequency-dispersive transient waves by means of the mild-slope equation. *Coastal Eng.* 55, 447–458.

Bellotti, G., Di Risio, M., and De Girolamo, P. (2009). Feasibility of tsunami early warning systems for small volcanic islands. *Nat. Hazards Earth Syst. Sci.* 9, 1911–1919. doi:10.5194/nhess-9-1911-2009

Bevilacqua, A., Neri, A., Bisson, M., Esposti Ongaro, T., Flandoli, F., Isaia, R., et al. (2017). The effects of vent location, event scale, and time forecasts on pyroclastic density current hazard maps at Campi Flegrei caldera (Italy). *Front. Earth Sci.* 5, 112. doi:10.3389/feart.2017.00072

Bilek, S. L. (2010). Invited review paper: seismicity along the South American subduction zone: review of large earthquakes, tsunamis, and subduction zone complexity. *Tectonophysics* 495 (1–2), 2–14. doi:10.1016/j.tecto.2009.02.037

Bilek, S. L., and Lay, T. (1999). Rigidity variations with depth along interplate megathrust faults in subduction zones. *Nature* 400, 443–446. doi:10.1038/22739

Bilek, S. L., and Lay, T. (2018). Subduction zone megathrust earthquakes. *Geosphere* 14, 1468–1500. doi:10.1130/GES01608.1

Bird, P., Jackson, D. D., Kagan, Y. Y., Kreemer, C., and Stein, R. S. (2015). GEAR1: a global earthquake activity rate model constructed from geodetic strain rates and smoothed seismicity. *Bull. Seismol. Soc. Amer.* 105, 2538–2554. doi:10.1785/0120150058

Bird, P., and Kagan, Y. Y. (2004). Plate-tectonic analysis of shallow seismicity: apparent boundary width, beta, corner magnitude, coupled lithosphere thickness, and coupling in seven tectonic settings. *Bull. Seismol. Soc. Amer.* 94 (6), 2380–2399. doi:10.1785/0120030107

Bird, P., and Kreemer, C. (2015). Revised tectonic forecast of global shallow seismicity based on version 2.1 of the global strain rate map. *Bull. Seismol. Soc. Amer.* 105, 152–166. doi:10.1785/0120140129

Birkmann, J., Cardona, O. D., Carreño, M. L., Barbat, A. H., Pelling, M., Schneiderbauer, S., et al. (2013). Framing vulnerability, risk and societal responses: the MOVE framework. *Nat. Hazards* 67 (2), 193–211. doi:10.1007/s11069-013-0558-5

Blake, D., Johnston, D., Leonard, G., McLaren, L., and Becker, J. (2018). A citizen science initiative to understand community response to the kaikoura earthquake and tsunami warning in petone and eastbourne, wellington, aotearoa/New Zealand. *Bull. Seismol. Soc. Amer.* 108 (3B), 1807–1817. doi:10.1785/0120170292

Blaser, L., Krüger, F., Ohrnberger, M., and Scherbaum, F. (2010). Scaling relations of earthquake source parameter estimates with special focus on subduction environment. *Bull. Seismol. Soc. Amer.* 100, 2914–2926. doi:10.1785/0120100111

Bletery, Q., Thomas, A. M., Rempel, A. W., Karlstrom, L., Sladen, A., and De Barros, L. (2016). Mega-earthquakes rupture flat megathrusts. *Science* 354 (6315), 1027–1031. doi:10.1126/science.aag0482

Blikra, L. H., Longva, O., Harbitz, C., and Løvholt, F. (2005). Quantification of rock-avalanche and tsunami hazard in Storfjorden, western Norway. *Landslides Avalanches: ICFL Conf. Proc.* 12, 57–64. doi:10.1007/springerreference_5010

Bommer, J. J., and Stafford, P. J. (2016). *Seismic hazard and earthquake actions*, in *Seismic design of buildings to eurocode 8*. Editor A. Elghazouli (Boca Raton: CRC Press), 21–54. doi:10.1201/9781315368221-9

Bougouin, A., Paris, R., and Roche, O. (2020). Impact of fluidized granular flows into water: implications for tsunamis generated by pyroclastic flows. *J. Geophys. Res. Solid Earth* 7, e2019JB018954. doi:10.1002/asl2.471

Bozza, A., Asprone, D., and Manfredi, G. (2015). Developing an integrated framework to quantify resilience of urban systems against disasters. *Nat. Hazards* 78 (3), 1729–1748. doi:10.1007/s11069-015-1798-3

Brown, S. K., Jenkins, S. F., Sparks, R. S. J., Odbert, H., and Auker, M. R. (2017).
Volcanic fatalities database: analysis of volcanic threat with distance and victim
classification. *J. Appl. Volcanol.* 6, 15. doi:10.1186/s13617-017-0067-4
Bryn, P., Berg, K., Forsberg, C. F., Solheim, A., and Kvalstad, T. J. (2005). Explaining the storegga slide. *Mar. Pet. Geol.* 22 (1–2), 11–19. doi:10.1016/j.marpetgeo.2004.12.003

Budnitz, R. J., Apostolakis, G., and Boore, D. M. (1997). Nuclear regulatory commission, Washington: The Academy Press. Recommendations for probabilistic seismic hazard analysis: guidance on uncertainty and use of experts. Technical Report. doi:10.2172/479072

Bullard, G. K., Mulligan, R. P., Carreira, A., and Take, W. A. (2019). Experimental analysis of tsunamis generated by the impact of landslides with high mobility. *Coastal Eng.* 152, 103538. doi:10.1016/j.coastaleng.2019.103538

Burbidge, D., Cummins, P. R., Mleczko, R., and Thio, H. K. (2008). A probabilistic tsunami hazard assessment for western Australia. *Pure Appl. Geophys.* 165, 2059–2088. doi:10.1007/s00024-008-0421-x

Bürgmann, R., and Chadwell, D. (2014). Seafloor geodesy. *Annu. Rev. Earth Planet. Sci.* 42, 509–534. doi:10.1146/annurev-earth-060313-054953

Cardona, O. D. (2001). Estimación holística del riesgo sísmico utilizando sistemas dinámicos complejos. Ph.D. Thesis. Barcelona, Spain: Polytechnic University of Catalonia.

Cardona, O. D., Ordaz, M., Salgado-Gálvez, M. A., Barbat, A. H., and Carreño, M. L. (2018). Latin American and caribbean earthquakes in the GEM's earthquake consequences database (GEMECD). *Nat. Hazards* 93 (1), S113–S125. doi:10.1007/s11069-017-3087-9

Carlton, B., Vanneste, M., Forsberg, C. F., Knudsen, S., Løvholt, F., Kvalstad, T., et al. (2019). Geohazard assessment related to submarine instabilities in Bjørnafjorden, London. *Spec. Publ. Geol. Soc.* 477, 549–566. doi:10.1007/bf00305006

Carreño, M.-L., Cardona, O. D., and Barbat, A. H. (2007). Urban seismic risk evaluation: a holistic approach. *Nat. Hazards* 40 (1), 137–172. doi:10.1007/s11069-006-0008-8

Cecioni, C., and Bellotti, G. (2010). Inclusion of landslide tsunamis generation into a depth integrated wave model. *Nat. Hazards Earth Syst. Sci.* 10 (11), 2259–2268. doi:10.5194/nhess-10-2259-2010

Chagué-Goff, C., Andrew, A., Szczuciński, W., Goff, J., and Nishimura, Y. (2012). Geochemical signatures up to the maximum inundation of the 2011 Tohoku-oki tsunami—implications for the 869 AD Jogan and other palaeotsunamis. *Sediment. Geol.* 282, 65–77. doi:10.1016/j.sedgeo.2012.05.021

Chakraborty, J., Tobin, G. A., and Montz, B. E. (2005). Population evacuation: assessing spatial variability in geophysical risk and social vulnerability to natural hazards. *Nat. Hazards Rev.* 6 (1), 23–33. doi:10.1061/(asce)1527-6988(2005)6:1(23)

Charvet, I., Macabuag, J., and Rossetto, T. (2017). Estimating tsunami-induced building damage through fragility functions: critical review and research needs. *Front. Built Environ.* 3, 121. doi:10.3389/fbuil.2017.00036

Charvet, I., Suppasri, A., Kimura, H., Sugawara, D., and Imamura, F. (2015). A multivariate generalized linear tsunami fragility model for Kesennuma City based on maximum flow depths, velocities and debris impact, with evaluation of predictive accuracy. *Nat. Hazards* 79 (3), 2073–2099. doi:10.1007/s11069-015-1947-8

Chaytor, J. D., Uri, S., Solow, A. R., and Andrews, B. D. (2009). Size distribution of submarine landslides along the US Atlantic margin. *Mar. Geol.* 264 (1-2), 16–27. doi:10.1016/j.margeo.2008.08.007

Chock, G. (2016). Design for tsunami loads and effects in the ASCE 7-16 standard. *J. Struct. Eng.* 142 (11), 04016093. doi:10.1061/(ASCE)ST.1943-541X.0001565

Chock, G. Y., Robertson, I., and Riggs, H. R. (2011). “Tsunami structural design provisions for a new update of building codes and performance-based engineering,” in *Solutions to coastal disasters* 2011, Editors Wallendorf, L.A., Jones, C., and Ewing, L. (London: ASCE) 423–435. doi:10.1061/41185(417)38

Chock, G., Yu, G., Thio, H. K., and Lynett, P. J. (2016). Target structural reliability analysis for tsunami hydrodynamic loads of the ASCE 7 standard. *J. Struct. Eng.* 142, 04016092–4016112. doi:10.1061/(ASCE)ST10.1061/(asce)st.1943-541x.0001499

Choi, B. H., Pelinovsky, E., Kim, K. O., and Lee, J. S. (2003). Simulation of the trans-oceanic tsunami propagation due to the 1883 Krakatau volcanic eruption. *Nat. Hazards Earth Syst. Sci.* 3, 321–332. doi:10.5194/nhess-3-321-2003

Churchill, D. D., Houston, S. H., and Bond, N. A. (1995). The daytona beach wave of 3-4 July 1992: a shallow-water gravity wave forced by a propagating squall line. *Bull. Amer. Meteorol. Soc.* 76 (1), 21–32. doi:10.1175/1520-0477(1995)076<0021:tdbwoj>2.0.co;2

Clare, M., Chaytor, J., Dabson, O., Gamboa, D., Georgiopoulou, A., Eady, H., et al. (2019). A consistent global approach for the morphometric characterization of subaqueous landslides. *Geol. Soc. Lond. Spec. Publ.* 477 (1), 455–477. doi:10.1144/sp477.15

Cooke, R. M. (1991). *Experts in uncertainty: opinion and subjective probability in science*. Oxford: Oxford University Press.

Cornell, C. A. (1968). Engineering seismic risk analysis. *Bull. Seismol. Soc. Amer.* 58, 1583–1606.

Crosta, G. B., Imposimato, S., and Roddeman, D. (2016). Landslide spreading, impulse water waves and modelling of the Vajont rockslide. *Rock Mech. Rock Eng.* 49 (6), 2413–2436. doi:10.1007/s00603-015-0769-z

Cutter, S. L., Boruff, B. J., and Shirley, W. L. (2003). Social vulnerability to environmental hazards*. *Soc. Sci. Q.* 84 (2), 242–261. doi:10.1111/1540-6237.8402002

Daniell, J. E., Schaefer, A. M., and Wenzel, F. (2017). Losses associated with secondary effects in earthquakes. *Front. Built Environ.* 3, 133. doi:10.3389/fbuil.2017.00030

Dao, M. H., and Tkalich, P. (2007). Tsunami propagation modelling—a sensitivity study. *Nat. Hazards Earth Syst. Sci.* 7, 741–754. doi:10.5194/nhess-7-741-2007

Davies, G., and Griffin, J. (2018). The 2018 Australian probabilistic tsunami hazard assessment Record 2018/41. Canberra: Geoscience Australia. doi:10.1163/Record

Davies, G., Griffin, J., Løvholt, F., Glimsdal, S., Harbitz, C., Thio, H. K., et al. (2018). A global probabilistic tsunami hazard assessment from earthquake sources. *Geol. Soc. Lond. Spec. Publ.* 456, 219–244. doi:10.1144/SP456.5

Davies, G., and Griffin, J. (2019). Sensitivity of probabilistic tsunami hazard assessment to far-field earthquake slip complexity and rigidity depth-dependence: case study of Australia. *Pure Appl. Geophys.* 177, 1521–1548. doi:10.1007/s00024-019-02299-w

Davies, G., and Griffin, J. (2020). Sensitivity of probabilistic tsunami hazard assessment to far-field earthquake slip complexity and rigidity depth-dependence: case study of Australia. *Pure Appl. Geophys.* 177, 1521–1548. doi:10.1007/s00024-019-02299-w

Davies, G., Horspool, N., and Miller, V. (2015). Tsunami inundation from heterogeneous earthquake slip distributions: evaluation of synthetic source models. *J. Geophys. Res. Solid Earth* 120, 6431–6451. doi:10.1002/2015JB012272

Davies, G. (2019). Tsunami variability from uncalibrated stochastic earthquake models: tests against deep ocean observations 2006-2016. *Geophys. J. Int.* 218, 1939–1960. doi:10.1093/gji/ggz260

Davoudi, S., Shaw, K., Haider, L. J., Quinlan, A. E., Peterson, G. D., Wilkinson, C., et al. (2012). Resilience: a bridging concept or a dead end? “Reframing” resilience: challenges for planning theory and practice interacting traps: resilience assessment of a pasture management system in northern Afghanistan urban resilience: what does it mean in planning practice? Resilience as a useful concept for climate change adaptation? The politics of resilience for planning: a cautionary note. *Plann. Theor. Pract.* 13 (2), 299–333. doi:10.1080/14649357.2012.677124

Day, S. J., et al. (2015). Volcanic tsunamis, in *Encyclopedia of volcanoes*. Editor Sigurdsson, H. (London, UK: Academic Press), 993–1009.

De Bono, A., and Chatenoux, B. (2015). A global exposure model for GAR 2015. UNEP-GRID, GAR.

De Bono, A., and Mora, M. G. (2014). A global exposure model for disaster risk assessment. *Int. J. Disaster Risk Reduct.* 10, 442–451. doi:10.1016/j.ijdrr.2014.05.008

de la Asunción, M., Castro, M. J., Fernández-Nieto, E. D., Mantas, J. M., Acosta, S. O., and González-Vida, J. M. (2013). Efficient GPU implementation of a two waves TVD-WAF method for the two-dimensional one layer shallow water system on structured meshes. *Comput. Fluids* 80, 441–452. doi:10.1016/j.complfluid.2012.01.012

De Risi, R., and Goda, K. (2016). Probabilistic earthquake-tsunami multi-hazard analysis: application to the Tohoku region. *Jpn. Front. Built Environ.* 2, 25. doi:10.3389/fbuil.2016.00025

De Risi, R., Goda, K., Mori, N., and Yasuda, T. (2017). Bayesian tsunami fragility modeling considering input data uncertainty. *Stoch. Environ. Res. Risk Assess.* 31, 1253–1269. doi:10.1007/s00477-016-1230-x

De Risi, R., and Goda, K. (2017). Simulation-based probabilistic tsunami hazard analysis: empirical and robust hazard predictions. *Pure Appl. Geophys.* 174, 3083–3106. doi:10.1007/s00024-017-1588-9
De Risi, R., Jalayer, F., De Paola, F., Iervolino, I., Giugni, M., Topa, M. E., et al. (2013). Flood risk assessment for informal settlements. *Nat. Hazards* 69 (1), 1003–1032. doi:10.1007/s11069-013-0749-0

De Risi, R., Muhammad, A., and Goda, K. (2018). A way forward to resilient infrastructures against earthquake-tsunami multi-hazard in Routledge handbook of sustainable and resilient infrastructure. Editor P. Gardoni (Oxon, U.K.: Routledge), 180–196.

Del Zoppo, M., Di Ludovico, M., and Prota, A. (2021). Methodology for assessing the performance of RC structures with breakaway infill walls under tsunami inundation. *J Struct Eng.* 147 (2), 04020330. doi:10.1016/b978-0-12-107180-6. 50009-3

Del Zoppo, M., Rossetto, T., Di Ludovico, M., and Prota, A. (2019). "Assessing the effect of tsunami-induced vertical loads on RC frames," in Proceedings of the 1st fib Italy Symposium on Concrete and Concrete Structures. Parma, Italy.

Dueñas-Osorio, L., and Vemuru, S. M. (2009). Cascading failures in complex infrastructure systems. *Struct. Saf.* 31 (2), 157–167. doi:10.1016/j.strusafe.2008.06.007

Dusek, G., DiVeglio, C., Licate, L., Heilman, L., Kirk, K., Paternostro, C., et al. (2019). A meteotsunami climatology along the U.S. East coast. *Bull. Amer. Meteorol. Soc.* 100, 1329–1345. doi:10.1175/BAMS-D-18-0206.1

Global Volcanism Program. (2013). "Volcanoes of the world," Editor Venzke, E. (London: Smithsonian Institution) (Accessed Januray 5, 2021). doi:10.5479/si.GVP.VOTW4-2013

Eguchi, R. T., Eguchi, M. T., Bouabid, J., Koshimura, S., and Graf, W. P. (2014). "Benchmarking, validation and calibration of newly-developed HAZUS tsunami methodology." in 10th U.S. National conference on earthquake engineering: Frontiers of earthquake engineering. NCEE 2014. Anchorage: United States. doi:10.4231/D31834329

EM-DAT. (2020). The international disaster database CRED. Available at: https://www.emdat.be (Accessed November 4, 2020).

Esteva, L. (1967). "Criteria for the construction of spectra for seismic design," in Presented at third Panamerican symposium on structures. Caracas, Venezuela.

Fardis, M. N. (2009). *Seismic design, assessment and retrofitting of concrete buildings: based on EN-Eurocode 8*. Berlin: Springer. doi:10.1007/978-1-4020-9842-0

FEMA (2017b). Federal emergency management agency. Hazus tsunami model technical guidance. FEMA, fed. Emerg. Manag. Agency. Washington. DC. Available at: https://www.fema.gov/media-library-data/1511284000276-4f18206fb0c7bab3c5ecbbbdf504b9fd/Hazus_40_Tsunami_Tech_Manual.pdf

FEMA (2017a). Federal emergency management agency. Hazus tsunami model user guidance. FEMA, fed. Emerg. Manag. Agency. Washington. DC. Available at: https://www.fema.gov/media-library-data/1511284000292-4f18206fb0c7bab3c5ecbbbdf504b9fd/Hazus_40_Tsunami_User_Manual.pdf.

Flanagan, B. E., Gregory, E. W., Hallisey, E. J., Heitger, J. L., and Lewis, B. (2011). A social vulnerability index for disaster management. *J. Homeland Security Emerg. Manag.* 8 (1), 3. doi:10.2202/1547-7355.1792

Foster, J. H., Ericksen, T. L., and Bingham, B. (2020). Wave glider-enhanced vertical seafloor geodesy. *J. Atmos. Oceanic Technol.* 37, 417–427. doi:10.1175/JTECH-D-19-0095.1

Fraser, S. A., Wood, N. J., Johnston, D. M., Leonard, G. S., Greening, P. D., and Rossetto, T. (2014). Variable population exposure and distributed travel speeds in least-cost tsunami evacuation modelling. *Nat. Hazards Earth Syst. Sci.* 14 (11), 2975–2991. doi:10.5194/nhess-14-2975-2014

Freundt, A. (2003). Entrance of hot pyroclastic flows into the sea: experimental observations. *Bull. Volcanol.* 65 (2-3), 144–164. doi:10.1007/s00445-002-0250-1

Fritz, H. M., Hager, W. H., and Minor, H.-E. (2003). Landslide generated impulse waves. *Experiments in Fluids* 35 (6), 505–519. doi:10.1007/s00348-003-0659-0

Fukutani, Y., Moriguchi, S., Terada, K., Kotani, T., Otake, Y., and Kitano, T. (2019). Tsunami hazard and risk assessment for multiple buildings by considering the spatial correlation of wave height using copulas. *Nat. Hazards Earth Syst. Sci.* 19 (11), 2619–2634. doi:10.5194/nhess-19-2619-2019

Gailler, A., Hébert, H., Schindelé, F., and Reymond, D. (2018). Coastal amplification laws for the French tsunami warning center: numerical modeling and fast estimate of tsunami wave heights along the French riviera. *Pure Appl. Geophys.* 175 (4), 1429–1444. doi:10.1007/s00024-017-1713-9

Gallotti, G., Passaro, S., Armigliato, A., Zaniboni, F., Pagnoni, G., Wang, L., et al. (2020). Potential mass movements on the Palinuro volcanic chain (southern

Tyrrhenian Sea, Italy) and consequent tsunami generation. *J. Volcanol. Geothermal Res.* 404, 107025. doi:10.1016/j.jvolgeores.2020.107025

Gallotti, G., and Tinti, S. (2020). A new approach for landslide modeling: application to the Scilla 1783 tsunamigenic landslide, South Italy. *Pure Appl. Geophys.* 177, 3563–3576. doi:10.1007/s00024-020-02546-5

Gayer, G., Leschka, S., Nöhren, I., Larsen, O., and Günther, H. (2010). Tsunami inundation modelling based on detailed roughness maps of densely populated areas. *Nat. Hazards Earth Syst. Sci.* 10 (8), 1679–1687. doi:10.5194/nhess-10-1679-2010

Geist, E. L., and Bilek, S. L. (2001). Effect of depth-dependent shear modulus on tsunami generation along subduction zones. *Geophys. Res. Lett.* 28, 1315–1318. doi:10.1029/2000GL012385

Geist, E. L., Brink, U. S., and Parsons, U. S. (2019). Offshore landslide hazard curves from mapped landslide size distributions. *J. Geophys. Res. Solid Earth* 124 (4), 3320–3334. doi:10.1029/2018JB017236

Geist, E. L., Chaytor, J. D., Parsons, T., and Brink, U. T. (2013). Estimation of submarine mass failure probability from a sequence of deposits with age dates. *Geosphere* 9 (2), 287–298. doi:10.1130/ges00829.1

Geist, E. L., and Oglesby, D. D. (2014). Earthquake mechanism and seafloor deformation for tsunami generation, in *Encyclopedia of earthquake engineering*. Editors M. Beer, I. A. Kougïoumtzoglou, E. Patelli, and I. S.-K. Au (Berlin, Heidelberg: Springer), 1–17. doi:10.1007/978-3-642-36197-5_296-1

Geist, E. L., Oglesby, D. D., and Ryan, K. J. (2019). "Tsunamis: stochastic models of occurrence and generation mechanisms," in *Encyclopedia of complexity and systems science*, Editor R. A. Meyers (Berlin, Heidelberg: Springer), 1–30. doi:10.1007/978-3-642-27737-5_595-2

Geist, E. L., and Parsons, T. (2006). Probabilistic analysis of tsunami hazards*. *Nat. Hazards* 37, 277–314. doi:10.1007/s11069-005-4646-z

Geist, E. L., and Parsons, T. (2014). Undersampling power-law size distributions: effect on the assessment of extreme natural hazards. *Nat. Hazards* 72, 565–595. doi:10.1007/s11069-013-1024-0

Geist, E. L., ten Brink, U. S., and Gove, M. (2014). A framework for the probabilistic analysis of meteotsunamis. *Nat. Hazards* 74, 123–142. doi:10.1007/s11069-014-1294-1

Geist, E., and Lynett, P. (2014). Source processes for the probabilistic assessment of tsunami hazards. *Oceanog* 27 (2), 86–93. doi:10.5670/oceanog.2014.43

Genadt, A. (2019). Three lessons from Japan on architectural resilience. *Architectural Histories* 7 (1), 16. doi:10.5334/ah.393

Germanovich, L. N., Kim, S., and Puzrin, A. M. (2016). Dynamic growth of slip surfaces in catastrophic landslides. *Proc. R. Soc. A.* 472 (2185), 20150758. doi:10.1098/rspa.2015.0758

Gibbons, S. J., Lorito, S., Macias, J., Løvholt, F., Selva, J., Volpe, M., et al. (2020). Probabilistic tsunami hazard analysis: high performance computing for massive scale inundation simulations. *Front. Earth Sci.* 8, 623. doi:10.3389/feart.2020.591549

Giles, D., Gopinathan, D., Guillas, S., and Dias, F. (2020). Faster than real time tsunami warning with associated hazard uncertainties. *Front. Earth Sci.* 8, 66. doi:10.3389/feart.2020.597865

Glimsdal, S., Løvholt, F., Harbitz, C. B., Romano, F., Lorito, S., Orefice, S., et al. (2019). A new approximate method for quantifying tsunami maximum inundation height ProbabilityA new approximate method for quantifying tsunami maximum inundation height probability. Pure appl. *Pure Appl. Geophys.* 176, 3227–3246. doi:10.1007/s00024-019-02091-w

Glimsdal, S., Pedersen, G. K., Harbitz, C. B., and Løvholt, F. (2013). Dispersion of Tsunamis: does it really matter?. *Nat. Hazards Earth Syst. Sci.* 13, 1507–1526. doi:10.5194/nhess-13-1507-2013

Goda, K., and De Risi, R. (2018). Multi-hazard loss estimation for shaking and tsunami using stochastic rupture sources. *Int. J. Disaster Risk Reduct.* 28, 539–554. doi:10.1016/j.ijdrr.2018.01.002

Goda, K. (2020). Multi-hazard portfolio loss estimation for time-dependent shaking and tsunami hazards. *Front. Earth Sci.* 8, 512. doi:10.3389/feart.2020.592444

Goda, K., Petrone, C., De Risi, R., and Rossetto, T. (2017). Stochastic coupled simulation of strong motion and tsunami for the 2011 Tohoku, Japan earthquake. *Stoch. Environ. Res. Risk Assess.* 31 (9), 2337–2355. doi:10.1007/s00477-016-1352-1

Goda, K., and Song, J. (2019). Influence of elevation data resolution on tsunami loss estimation and insurance rate-making. *Front. Earth Sci.* 7, 246. doi:10.1016/
b978-0-12-107180-6_50007-x
Goda, K. (2019). Time-dependent probabilistic tsunami hazard analysis using stochastic rupture sources. *Stoch. Environ. Res. Risk Assess.* 33, 341–358. doi:10.1007/s00477-018-1634-x

Goda, K., Yasuda, T., Mori, N., and Maruyama, T. (2016). New scaling relationships of earthquake source parameters for stochastic tsunami simulation. *Coastal Eng.* J. 58, 1650010–1650011. doi:10.1142/S0578563416500108

González, F. I., Geist, E. L., Jaffe, B., Kânoğlu, U., Mofjeld, H., Synolakis, C. E., et al. (2009). Probabilistic tsunami hazard assessment at Seaside, Oregon, for near- and far-field seismic sources. *J. Geophys. Res.* 114, 37. doi:10.1029/2008JC005132

González, J. I., Farreras, S. F., and Ochoa, J. (2001). Seismic and meteorological tsunami contributions in the Manzanillo and Cabo San Lucas Seiches of September 14, 1995. *Mar. Geodesy* 24 (4), 219–227. doi:10.1080/014904101753227860

Greenslade, D. J. M., Annunziato, A., Babeyko, A. Y., Burbidge, D. R., Ellguth, E., Horspool, N., et al. (2014). An assessment of the diversity in scenario-based tsunami forecasts for the Indian Ocean. *Continental Shelf Res.* 79, 36–45. doi:10.1016/j.csr.2013.06.001

Gregg, C. E., Houghton, B. F., Paton, D., Lachman, R., Lachman, J., Johnston, D. M., et al. (2006). Natural warning signs of tsunamis: human sensory experience and response to the 2004 great Sumatra earthquake and tsunami in Thailand. *Thai. Earthq. Spectra* 22 (S3), 671–691. doi:10.1193/1.2206791

Grezio, A., Babeyko, A., Baptista, M. A., Behrens, J., Costa, A., Davies, G., et al. (2017). Probabilistic tsunami hazard analysis: multiple sources and global applications. *Rev. Geophys.* 55, 1158–1198. doi:10.1002/2017RG000579

Grezio, A., Cinti, F. R., Costa, A., Faenza, L., Perfetti, P., Pierdominici, S., et al. (2020). Multisource bayesian probabilistic tsunami hazard analysis for the Gulf of Naples (Italy). *J. Geophys. Res. Oceans* 125, 12. doi:10.1029/2019JC015373

Grezio, A., Marzocchi, W., Sandri, L., and Gasparini, P. (2010). A bayesian procedure for probabilistic tsunami hazard assessment. *Nat. Hazards* 53, 159–174. doi:10.1007/s11069-009-9418-8

Griffin, J., Latief, H., Kongko, W., Harig, S., Horspool, N., Hanung, R., et al. (2015). An evaluation of onshore digital elevation models for modeling tsunami inundation zones. *Front. Earth Sci.* 3, 32. doi:10.3389/feart.2015.00032

Grilli, S. T., Taylor, O. D. S., Baxter, C. D., and Maretzki, S. (2009). A probabilistic approach for determining submarine landslide tsunami hazard along the upper east coast of the United States. *Mar. Geol.* 264 (1–2), 74–97. doi:10.1016/j.margeo.2009.02.010

Grilli, S. T., and Watts, P. (2005). Tsunami generation by submarine mass failure: I. modeling, experimental validation, and sensitivity analyses. *J. Waterway, Port, Coastal, Ocean Eng.* 131 (6), 283–297. doi:10.1061/(asce)0733-950x(2005)131:6(283)

GTM (2020). Global Tsunami Model. http://edanya.uma.es/gtm/ (Accessed November 12, 2020).

Gylfadóttir, S. S., Kim, J., Helgason, J. K., Brynjólfsson, S., Höskuldsson, Á., Jóhannesson, T., et al. (2017). The 2014 Lake Askja rockslide-induced tsunami: optimization of numerical tsunami model using observed data. *J. Geophys. Res. Oceans* 122 (5), 4110–4122. doi:10.1002/2016JC012496

Hananto, N. D., Leclerc, F., Li, L., Etchebes, M., Carton, H., Tapponnier, P., et al. (2020). Tsunami earthquakes: vertical pop-up expulsion at the forefront of subduction megathrust. *Earth Planet. Sci. Lett.* 538, 116197. doi:10.1016/j.epsl.2020.116197

Harbitz, C. B. (1992). Model simulations of tsunamis generated by the Storegga slides. *Mar. Geol.* 105 (1-4), 1–21. doi:10.1016/0025-3227(92)90178-k

Harbitz, C. B., Lovholt, F., and Bungum, H. (2014). Submarine landslide tsunamis: how extreme and how likely?. *Nat. Hazards* 72 (3), 1341–1374. doi:10.1007/s11069-013-0681-3

Harig, S., Chaeroni, Pranowo, W. S., and Behrens, J. (2008). Tsunami simulations on several scales. *Ocean Dyn.* 58, 429–440. doi:10.1007/s10236-008-0162-5

Haslett, S. K., Mellor, H. E., and Bryant, E. A. (2009). Meteo-tsunami hazard associated with summer thunderstorms in the United Kingdom. *Phys. Chem. Earth, Parts A/B/C* 34, 1016–1022. doi:10.1016/j.pce.2009.10.005

Hatayama, K. (2014). “Damage to oil storage tanks due to tsunami of the Mw9.0 2011 off the Pacific coast of Tohoku, Japan.” in 10th U.S. National conference on earthquake engineering, frontiers of earthquake engineering. Anchorage, USA: Earthquake Engineering Research Institute. doi:10.4231/D3VX0639C

Hayes, G. P., Moore, G. L., Portner, D. E., Hearne, M., Flamme, H., Furtney, M., et al. (2018). Slab2, a comprehensive subduction zone geometry model. *Science* 362 (6410), 58–61. doi:10.1126/science.aat4723

Heller, V., and Hager, W. H. (2010). Impulse product parameter in landslide generated impulse waves. *J. Waterway, Port, Coastal, Ocean Eng.* 136 (3), 145–155. doi:10.1061/(asce)ww.1943-5460.0000037

Heller, V. (2011). Scale effects in physical hydraulic engineering models. *J. Hydraulic Res.* 49 (3), 293–306. doi:10.1080/00221686.2011.578914

Hermanns, R., Oppikofer, T., Anda, E., Blikra, L., Böhme, M., Bunkholt, H., et al. (2013). Hazard and risk classification for large unstable rock slopes in Norway. *It. J. Eng. Geol. Env.* 11, 28. doi:10.4408/IJEGE.2013.06-B.22

Herrero, A., and Bernard, P. (1994). A kinematic self-similar rupture process for earthquakes. *Bull. Seismol. Soc. Amer.* 84, 1216–1228.

Herslund, L. B., Jalayer, F., Jean-Baptiste, N., Jørgensen, G., Kabisch, S., Kombe, W., et al. (2016). A multi-dimensional assessment of urban vulnerability to climate change in Sub-Saharan Africa. *Nat. hazards* 82 (2), 149–172. doi:10.1007/s11069-015-1856-x

Heuret, A., Conrad, C. P., Funiciello, F., Lallemand, S., and Sandri, L. (2012). Relation between subduction megathrust earthquakes, trench sediment thickness and upper plate strain. *Geophys. Res. Lett.* 39, a. doi:10.1029/2011GL050712

Hoechner, A., Babeyko, A. Y., and Zamora, N. (2016). Probabilistic tsunami hazard assessment for the Makran region with focus on maximum magnitude assumption. *Nat. Hazards Earth Syst. Sci.* 16, 1339–1350. doi:10.5194/nhess-16-1339-2016

Horspool, N., Pranantyo, I., Griffin, J., Latief, H., Natawidjaja, D. H., Kongko, W., et al. (2014). A probabilistic tsunami hazard assessment for Indonesia. *Nat. Hazards Earth Syst. Sci.* 14, 3105–3122. doi:10.5194/nhess-14-3105-2014

Huang, X., and García, M. H. (1998). A Herschel-Bulkley model for mud flow down a slope. *J. Fluid Mech.* 374, 305–333. doi:10.1017/S0022112098002845

Hufschmidt, G., Crozier, M., and Glade, T. (2005). Evolution of natural risk: research framework and perspectives. *Nat. Hazards Earth Syst. Sci.* 5 (3), 375–387. doi:10.5194/nhess-5-375-2005

Huhn, K., Arroyo, M., Cattaneo, A., Clare, M. A., Gràcia, E., Harbitz, C. B., et al. (2019). “Modern submarine landslide complexes: a short review,” in *Submarine. Landslides: Subaqueous mass transport deposits from outcrops to seismic profiles*. Editor Ogata, K. (Hoboken, NJ, USA:Wiley), 181–200.

Huthnance, J. M. (1975). On trapped waves over a continental shelf. *J. Fluid Mech.* 69 (4), 689–704. doi:10.1017/S0022112075001632

Huyck, C., Hu, Z., Amyx, P., Esquivias, M. H., Eguchi, M., Contributors, K., et al. (2019). METEOR: exposure data classification. *Metadata Popul. Confidence Assess. Rep. Number* 3, 2. doi:10.3403/30126912u

Huyck, C. K., and Eguchi, M. T. (2017). GFDRR Africa disaster risk financing-result area 5 exposure development. Replace. Cost Refinements Expo. Data. Prepared for World Bank/GFDRR.

Huyck, C. K., Esquivias, G., Gamba, P., Hussain, M., Odhiambo, O., Jaiswal, K., et al. (2011). D2.2 survey of available input databases for GED, report produced in the context of the global exposure databases for the global earthquake model (IDCT). *Constraints Databases* 4, 139–161. doi:10.1007/978-1-4615-5515-5_7

ITIC (2020). International tsunami information center. Available at: http://itic.iocunesco.org/ (Accessed November 11, 2020).

Jaimes, M. A., Reinoso, E., Ordaz, M., Huerta, B., Silva, R., Mendoza, E., et al. (2016). A new approach to probabilistic earthquake-induced tsunami risk assessment. *Ocean Coastal Manage.* 119, 68–75. doi:10.1016/j.ocecoaman.2015.10.007

Jalayer, F., Carozza, S., De Risi, R., Manfredi, G., and Mbuya, E. (2016). Performance-based flood safety-checking for non-engineered masonry structures. *Eng. Struct.* 106, 109–123. doi:10.1016/j.engstruct.2015.10.007

Jalayer, F., and Ebrahimian, H. (2020). Seismic reliability assessment and the nonergodicity in the modelling parameter uncertainties. *Earthq. Eng. Struct. Dyn.* 49 (5), 434–457. doi:10.1002/eqe.3247

Jasanoff, S. (1993). Bridging the two cultures of risk Analysis1,2. *Risk Anal.* 13 (2), 123–129. doi:10.1111/j.1539-6924.1993.tb01057.x

Jelinek, R., and Krausmann, E. (2008). Approaches to tsunami risk assessment. *JRC Sci. Tech. Rep.* 48713, 112. doi:10.4324/9781351140843-3

Johnston, D., Tarrant, R., Tipker, K., Lambie, E., Crawford, M., Becker, J., et al. (2016). Towards tsunami-safer schools in the Wellington Region, New Zealand:
evaluating drills and awareness programmes. *Aust. J. Emerg. Manage.* 31, 59–66. doi:10.7717/peerj-cs.386/table-1

Jop, P., Forterre, Y., and Pouliquen, O. (2006). A constitutive law for dense granular flows. *Nature* 441 (7094), 727–730. doi:10.1038/nature04801

JRC Risk Data Hub (2020). Available at: https://ec.europa.eu/jrc/en (Accessed November 11, 2020).

Kagan, Y. Y., and Jackson, D. D. (1991). Long-term earthquake clustering. *Geophys. J. Int.* 104, 117–134. doi:10.1111/j.1365-246X.1991.tb02498.x

Kagan, Y. Y., and Jackson, D. D. (2014). Statistical earthquake focal mechanism forecasts. *Geophys. J. Int.* 197, 620–629. doi:10.1093/gji/ggu015

Kagan, Y. Y. (2002). Seismic moment distribution revisited: I. Statistical results. *Geophys. J. Int.* 148, 520–541. doi:10.1046/j.1365-246x.2002.01594.x

Kaiser, G., Scheele, L., Kortenhaus, A., Lovholt, F., Römer, H., and Leschka, S. (2011). The influence of land cover roughness on the results of high resolution tsunami inundation modeling. *Nat. Hazards Earth Syst. Sci.* 11 (9), 2521–2540. doi:10.5194/nhess-11-2521-2011

Kajijiura, K. (1963). The leading wave of a tsunami. *B. Earthq. Res. Inst.* 41, 535–571.

Karafagka, S., Fotopoulou, S., and Pitilakis, K. (2018). Analytical tsunami fragility curves for seaport RC buildings and steel light frame warehouses. *Soil Dyn. Earthq. Eng.* 112, 118–137. doi:10.1016/j.soildyn.2018.04.037

Kelfoun, K., and Druitt, T. H. (2005). Numerical modeling of the emplacement of socompa rock avalanche, Chile. *J. Geophys Res: solid Earth*. 110 (B12), 113. doi:10.3133/jfr01383

Kienle, J., Kowalik, Z., and Murty, T. S. (1987). Tsunamis generated by eruptions from mount st. Augustine volcano, Alaska. *Science* 236, 1442–1447. doi:10.1126/science.236.4807.1442

Kim, D. H., Lynett, P. J., and Socolofsky, S. A. (2009). A depth-integrated model for weakly dispersive, turbulent, and rotational fluid flows. *Ocean Model.* 27 (3–4), 198–214. doi:10.1016/j.ocemod.2009.01.005

Kim, J., Lovholt, F., Issler, D., and Forsberg, C. F. (2019). Landslide material control on tsunami genesis-the storegga slide and tsunami (8,100 Years BP). *J. Geophys. Res. Oceans* 124 (6), 3607–3627. doi:10.1029/2018jc014893

King, G. C. P., Stein, R. S., and Lin, J. (1994). Static stress changes and the triggering of earthquakes. *Bull. Seismol. Soc. Am.* 84, 935–953.

Kiureghian, A. D., and Ditlevsen, O. (2009). Aleatory or epistemic? Does it matter? *Struct. Saf.* 31 (2), 105–112. doi:10.1016/j.strusafe.2008.06.020

Kiureghian, A. D. (2005). Non-ergodicity and PEER's framework formula. *Earthq. Eng. Struct. Dyn.* 34 (13), 1643–1652. doi:10.1002/eqe.504

Kopp, H. (2013). Invited review paper: the control of subduction zone structural complexity and geometry on margin segmentation and seismicity. *Tectonophysics* 589, 1–16.

Koshimura, S., Oie, T., Yanagisawa, H., and Imamura, F. (2009). Developing fragility functions for tsunami damage estimation using numerical model and post-tsunami data from Banda Aceh, Indonesia. *Coastal Eng. J.* 51 (3), 243–273. doi:10.1142/s0578563409002004

Kozdon, J. E., and Dunham, E. M. (2013). Rupture to the trench: dynamic rupture simulations of the 11 march 2011 Tohoku earthquake. *Bull. Seismol. Soc. Amer.* 103, 1275–1289. doi:10.1785/0120120136

Kriebel, D. L., Lynett, P. J., Cox, D. T., Petroff, C. M., Robertson, I. N., and Chock, G. Y. (2017). Energy method for approximating overland tsunami flows. *J. Waterway Port Coastal Ocean Eng.* 143 (5),04017014. doi:10.1061/(ASCE) WW.1943-5460.0000393

Lam, N. S., Reams, M., Li, K., Li, C., and Mata, L. P. (2016). Measuring community resilience to coastal hazards along the northern Gulf of Mexico. *Nat. Hazards Rev.* 17 (1), 113. doi:10.1061/(ASCE)NH.1527-6996.0000193

Lane, E. M., Mountjoy, J. J., Power, W. L., and Mueller, C. (2016). “Probabilistic hazard of tsunamis generated by submarine landslides in the cook strait canyon (New Zealand), in *Global tsunami science Past and Future*. Editor Geist, H. (Basel, Switzerland: birkhäuser), Vol. I, 3757–3774.

Latter, J. H. (1981). Tsunamis of volcanic origin: summary of causes, with particular reference to Krakatoa, 1883. *Bull. Volcanol* 44, 467–490. doi:10.1007/bf02600578

Lee, H. J. (2009). Timing of occurrence of large submarine landslides on the Atlantic Ocean margin. *Mar. Geol.* 264 (1-2), 53–64. doi:10.1016/j.margeo. 2008.09.009

LeMéhauté, B. L., and Wang, S. (1996). *Water waves generated by underwater explosion. Advanced series on ocean engineering*. New Jersey: World Scientific.

Li, S., Moreno, M., Rosenau, M., Melnick, D., and Oncken, O. (2014). Splay fault triggering by great subduction earthquakes inferred from finite element models. *Geophys. Res. Lett.* 41, 385–391. doi:10.1002/2013GL058598

Lin, I.-C., and Tung, C. C. (1982). A preliminary investigation of tsunami hazard. *Bull. Seismol. Soc. Amer.* 72, 2323–2337.

Logan, D. C. (2009). Known knowns, known unknowns, unknown unknowns and the propagation of scientific enquiry. *J. Exp. Bot.* 60 (3), 712–714. doi:10.1093/ jxb/erp043

Lorito, S., Romano, F., and Lay, T. (2016). “Tsunamigenic major and great earthquakes (2004–2013): source processes inverted from seismic, geodetic, and sea-level data,” in *Encyclopedia of complexity and systems science*. Editor R. A. Meyers (Berlin, Heidelberg: Springer), 1–52. doi:10.1007/978-3-642- 27737-5_641-1

Lorito, S., Selva, J., Basili, R., Romano, F., Tiberti, M. M., and Piatanesi, A. (2015). Probabilistic hazard for seismically induced tsunamis: accuracy and feasibility of inundation maps. *Geophys. J. Int.* 200 (1), 574–588. doi:10.1093/gji/ggu408

Lotto, G. C., Dunham, E. M., Jeppson, T. N., and Tobin, H. J. (2017). The effect of compliant prisms on subduction zone earthquakes and tsunamis. *Earth Planet. Sci. Lett.* 458, 213–222. doi:10.1016/j.epsl.2016.10.050

Lotto, G. C., Jeppson, T. N., and Dunham, E. M. (2019). Fully coupled simulations of megathrust earthquakes and tsunamis in the Japan trench, nankai trough, and cascadia subduction zone. *Pure Appl. Geophys.* 176, 4009–4041. doi:10. 1007/s00024-018-1990-y

Loveless, J. P., and Meade, B. J. (2011). Spatial correlation of interseismic coupling and coseismic rupture extent of the 2011 MW= 9.0 Tohoku-oki earthquake. *Geophys. Res. Lett.* 38 (17), 39. doi:10.1029/2011GL048561

Lovholt, F., Bondevik, S., Laberg, J. S., Kim, J., and Boylan, N. (2017). Some giant submarine landslides do not produce large tsunamis. *Geophys. Res. Lett.* 44 (16), 8463–8472. doi:10.1002/2017GL074062

Lovholt, F., Fraser, S., Salgado-Galvez, M., Lorito, S., Selva, J., Romano, F., et al. (2019). Global trends in advancing tsunami science for improved hazard and risk understanding. Contributing Paper to GAR19, June. doi:10.1109/ urgenthpc49580.2019.00011

Lovholt, F., Glimsdal, S., and Harbitz, C. B. (2020). On the landslide tsunami uncertainty and hazard. *Landslides* 17, 2301–2315. doi:10.1007/s10346-020- 01429-z

Lovholt, F., Glimsdal, S., Harbitz, C. B., Zamora, N., Nadim, F., Peduzzi, P., et al. (2012). Tsunami hazard and exposure on the global scale. *Earth-Sci. Rev.* 110 (1-4), 58–73. doi:10.1007/978-3-642-27737-5_642-1

Lovholt, F., Griffin, J., and Salgado-Gálvez, M. (2015a). Tsunami hazard and risk assessment on the global scale, in *Encyclopedia of complexity and systems science*. Editor R. A. Meyers (Berlin, Heidelberg: Springer), 1–34. doi:10.1007/ 978-3-642-27737-5642-1

Lovholt, F., Harbitz, C. B., and Haugen, K. B. (2005). A parametric study of tsunamis generated by submarine slides in the Ormen Lange/Storegga area off western Norway. *Mar. Pet. Geol.* 22 (1–2), 219–231. doi:10.1016/b978-0-08- 044694-3.50023-8

Lovholt, F., Pedersen, G., and Gisler, G. (2008). Oceanic propagation of a potential tsunami from the La palma island. *J Geophys Res: Oceans.* 113 (C9), 38. doi:10. 1029/2007JC004603

Lovholt, F., Pedersen, G., Harbitz, C. B., Glimsdal, S., and Kim, J. (2015b). On the characteristics of landslide tsunamis. *Phil. Trans. R. Soc. A.* 373 (2053), 20140376. doi:10.1098/rsta.2014.0376

Lovholt, F., Setiadi, N. J., Birkmann, J., Harbitz, C. B., Bach, C., Fernando, N., et al. (2014). Tsunami risk reduction–are we better prepared today than in 2004? *Int. J. Disaster Risk reduction* 10, 127–142.

Lynett, P. J., Gately, K., Wilson, R., Montoya, L., Arcas, D., Aytore, B., et al. (2017). Inter-model analysis of tsunami-induced coastal currents. *Ocean Model.* 114, 14–32. doi:10.1016/j.ocemod.2017.04.003

Ma, S. (2012). A self-consistent mechanism for slow dynamic deformation and tsunami generation for earthquakes in the shallow subduction zone. *Geophys. Res. Lett.* 39 (11), 12. doi:10.1029/2012GL051854

Macabuag, J., Rossetto, T., Ioannou, I., Suppasri, A., Sugawara, D., Adriano, B., et al. (2016). A proposed methodology for deriving tsunami fragility functions for buildings using optimum intensity measures. *Nat. Hazards* 84 (2), 1257–1285. doi:10.1007/s11069-016-2485-8

Macías, J., Castro, M. J., and Escalante, C. (2020). Performance assessment of the Tsunami-HySEA model for NTHMP tsunami currents benchmarking.
Laboratory data. *Coastal Eng.* 158, 103667. doi:10.1016/j.coastaleng.2020.103667

Macías, J., Castro, M. J., Ortega, S., Escalante, C., and González-Vida, J. M. (2017). Performance benchmarking of tsunami-HySEA model for NTHMP's inundation mapping activities. *Pure Appl. Geophys.* 174, 3147–3183. doi:10.1007/s00024-017-1583-1

Maeno, F., and Imamura, F. (2011). Tsunami generation by a rapid entrance of pyroclastic flow into the sea during the 1883 Krakatau eruption, Indonesia. *J Geophys Res.* 116 (B09), 205. doi:10.1029/2011JB008253

Mai, P. M., and Beroza, G. C. (2002). A spatial random field model to characterize complexity in earthquake slip. *J. Geophys. Res.* 107, 10–11. doi:10.1029/2001JB000588

Mansinha, L., and Smylie, D. E. (1971). The displacement of the earthquake fault model. *Bull. Seismol. Soc. Amer.* 61, 1433–1440.

Marras, S., and Mandli, K. T. (2021). Modeling and simulation of tsunami impact: a short review of recent advances and future challenges. *Geosciences* 11 (1), 5. doi:10.3390/geosciences11010005

Marzocchi, W., and Jordan, T. H. (2014). Testing for ontological errors in probabilistic forecasting models of natural systems. *Proc. Natl. Acad. Sci.* 111 (33), 11973–11978. doi:10.1073/pnas.1410183111

Mas, E., Koshimura, S., Suppasri, A., Matsuoka, M., Matsuyama, M., Yoshii, T., et al. (2012). Developing Tsunami fragility curves using remote sensing and survey data of the 2010 Chilean Tsunami in Dichato. *Nat. Hazards Earth Syst. Sci.* 12, 2689–2697. doi:10.5194/nhess-12-2689-2012

Masterlark, T. (2003). Finite element model predictions of static deformation from dislocation sources in a subduction zone: sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions. *J. Geophys. Res.* 108, 12. doi:10.1029/2002JB002296

McGuire, R. K. (2008). Probabilistic seismic hazard analysis: early history. *Earthq. Eng. Struct. Dyn.* 37, 329–338. doi:10.1002/eqe.765

McLaughlin, S., and Cooper, J. A. G. (2010). A multi-scale coastal vulnerability index: a tool for coastal managers?. *Environ. Hazards* 9 (3), 233–248. doi:10.3763/ehaz.2010.0052

Meade, B. J. (2007). Algorithms for the calculation of exact displacements, strains, and stresses for triangular dislocation elements in a uniform elastic half space. *Comput. Geosci.* 33, 1064–1075. doi:10.1016/j.cageo.2006.12.003

Mebarki, A., Jerez, S., Prodhomme, G., and Reimeringer, M. (2016). Natural hazards, vulnerability and structural resilience: tsunamis and industrial tanks. *Geomat. Nat. Hazards Risk* 7 (Suppl. 1), 5–17. doi:10.1080/19475705.2016.1181458

Meletti, C., Marzocchi, W., D'Amico, V., Lanzano, G., Luzi, L., and Martinelli, F., et al.; MPS19 Working Group (2021). The new Italian seismic hazard model (MPS19). *Ann. Geophys.* 64 (1), SE112. doi:10.4401/ag-8579

Melgar, D., LeVeque, R. J., Dreger, D. S., and Allen, R. M. (2016). Kinematic rupture scenarios and synthetic displacement data: an example application to the Cascadia subduction zone. *J. Geophys. Res. Solid Earth* 121, 6658–6674. doi:10.1002/2016JB013314

Melnick, D., Moreno, M., Quinteros, J., Baez, J. C., Deng, Z., Li, S., et al. (2017). The super-interseismic phase of the megathrust earthquake cycle in Chile. *Geophys. Res. Lett.* 44, 784–791. doi:10.1002/2016GL071845

Métois, M., Socquet, A., and Vigny, C. (2012). Interseismic coupling, segmentation and mechanical behavior of the central Chile subduction zone. *J. Geophys. Res.* 117, 18. doi:10.1029/2011JB008736

Miyashita, T., Mori, N., and Goda, K. (2020). Uncertainty of probabilistic tsunami hazard assessment of Zihuatanejo (Mexico) due to the representation of tsunami variability. *Coastal Eng. J.* 62, 413–428. doi:10.1080/21664250.2020.1780676

Molinari, I., Tonini, R., Lorito, S., Piatanesi, A., Romano, F., Melini, D., et al. (2016). Fast evaluation of tsunami scenarios: uncertainty assessment for a Mediterranean Sea database. *Nat. Hazards Earth Syst. Sci.* 16 (12), 2593–2602. doi:10.5194/nhess-16-2593-2016

Monserrat, S., Vilibić, I., and Rabinovich, A. B. (2006). Meteotsunamis: atmospherically induced destructive ocean waves in the tsunami frequency band. *Nat. Hazards Earth Syst. Sci.* 6, 1035–1051. doi:10.5194/nhess-6-1035-2006

Moreno, M., Haberland, C., Oncken, O., Rietbrock, A., Angiboust, S., and Heidbach, O. (2014). Locking of the Chile subduction zone controlled by

fluid pressure before the 2010 earthquake. *Nat. Geosci.* 7, 292–296. doi:10.1038/ngeo2102

Morgan, M. G. (2014). Use (and abuse) of expert elicitation in support of decision making for public policy. *Proc. Natl. Acad. Sci.* 111 (20), 7176–7184. doi:10.1073/pnas.1319946111

Muldashev, I. A., and Sobolev, S. V. (2020). What controls maximum magnitudes of giant subduction earthquakes?. *Geochem. Geophys. Geosyst.* 21, e2020GC009145. doi:10.1029/2020GC009145

Murotani, S., Miyake, H., and Koketsu, K. (2008). Scaling of characterized slip models for plate-boundary earthquakes. *Earth Planet. Sp* 60, 987–991. doi:10.1186/BF03352855

Murotani, S., Satake, K., and Fujii, Y. (2013). Scaling relations of seismic moment, rupture area, average slip, and asperity size for M ~9 subduction-zone earthquakes. *Geophys. Res. Lett.* 40, 5070–5074. doi:10.1002/grl.50976

Murphy, S., Di Toro, G., Romano, F., Scala, A., Lorito, S., Spagnuolo, E., et al. (2018). Tsunamigenic earthquake simulations using experimentally derived friction laws. *Earth Planet. Sci. Lett.* 486, 155–165. doi:10.1016/j.epsl.2018.01.011

Murphy, S., Scala, A., Herrero, A., Lorito, S., Festa, G., Trasatti, E., et al. (2016). Shallow slip amplification and enhanced tsunami hazard unravelled by dynamic simulations of mega-thrust earthquakes. *Sci. Rep.* 6, 35007. doi:10.1038/srep35007

Musa, A., Watanabe, O., Matsuoka, H., Hokari, H., Inoue, T., Murashima, Y., et al. (2018). Real-time tsunami inundation forecast system for tsunami disaster prevention and mitigation. *J. Supercomput.* 74, 3093–3113. doi:10.1007/s11227-018-2363-0

Nadim, F., and Glade, T. (2006). "On tsunami risk assessment for the west coast of Thailand," in *Geohazards*. Available at: https://dc.engconfintl.org/geohazards/28/

Nes, N. (2018). Past slope collapse and current unstable slopes along a cliff at Hyllestad Kommune. MSc thesis, Norway: Technical University Science.

Nikkhoo, M., and Walter, T. R. (2015). Triangular dislocation: an analytical, artefact-free solution. *Geophys. J. Int.* 201, 1119–1141. doi:10.1093/gji/ggy035

Nistor, I., Goseberg, N., and Stolle, J. (2017b). Tsunami-driven debris motion and loads: a critical review. *Front. Built Environ.* 3, 2. doi:10.3389/fbuil.2017.00002

Nistor, I., Goseberg, N., Stolle, J., Mikami, T., Shibayama, T., Nakamura, R., et al. (2017a). Experimental investigations of debris dynamics over a horizontal plane. *J. Waterway Port Coastal Ocean Eng.* 143, 04016022. doi:10.2139/srn.2277688

Nomanbhoy, N., and Satake, K. (1995). Generation mechanism of tsunamis from the 1883 Krakatau eruption. *Geophys. Res. Lett.* 22, 509–512. doi:10.1029/94gl03219

Nosov, M. A., and Kolesov, S. V. (2007). Elastic oscillations of water column in the 2003 Tokachi-oki tsunami source: *in-situ* measurements and 3-D numerical modelling. *Nat. Hazards Earth Syst. Sci.* 7, 243–249. doi:10.5194/nhess-7-243-2007

NTHMP. (2018). National tsunami hazard mitigation program Landslide tsunami modelling benchmark workshop—Benchmark problems. http://www1.udel.edu/kirby/landslide/problems.html

O'Brien, K., Eriksen, S., Nygaard, L. P., and Schjolden, A. (2007). Why different interpretations of vulnerability matter in climate change discourses. *Clim. Pol.* 7 (1), 73–88. doi:10.1080/14693062.2007.9685639

Oishi, Y., Imamura, F., and Sugawara, D. (2015). Near-field tsunami inundation forecast using the parallel TUNAMI-N2 model: application to the 2011 Tohoku-Oki earthquake combined with source inversions. *Geophys. Res. Lett.* 42, 1083–1091. doi:10.1002/2014GL062577

Okada, Y. (1992). Internal deformation due to shear and tensile faults in a half-space. *Bull. Seismol. Soc. Am.* 82, 1018–1040.

Okal, E. A. (2015). The quest for wisdom: lessons from 17 tsunamis, 2004–2014. *Phil. Trans. R. Soc. A.* 373 (2053), 20140370. doi:10.1098/rsta.2014.0370

Okumura, N., Jonkman, S. N., Esteban, M., Hofland, B., and Shibayama, T. (2017). A method for tsunami risk assessment: a case study for Kamakura, Japan. *Nat. Hazards* 88 (3), 1451–1472. doi:10.1007/s11069-017-2928-x

OpenStreetMap. (2020). Available at: https://www.openstreetmap.org/ (Accessed November 1, 2020).

Ordaz, M. (2000). Metodología para la evaluación del riesgo sísmico enfocada a la gerencia de seguros por terremoto. Universidad Nacional Autónoma de México. México DF.
Ordaz, M. (2015). A simple probabilistic model to combine losses arising from the simultaneous occurrence of several hazards. *Nat. Hazards* 76, 389–396. doi:10.1007/s11069-014-1495-7

Ordaz, M., Salgado-Gálvez, M. A., Huerta, B., Rodríguez, J. C., and Avelar, C. (2019). Considering the impacts of simultaneous perils. *Dpm* 28 (6), 823–837. doi:10.1108/DPM-09-2019-0295

Ota, Y., and Yamaguchi, M. (2004). Holocene coastal uplift in the western Pacific Rim in the context of late Quaternary uplift. *Quat. Int.* 120 (1), 105–117. doi:10.1016/j.quaint.2004.01.010

Páez-Ramírez, J., Lizarazo-Marriaga, J., Medina, S., Estrada, M., Mas, E., and Koshimura, S. (2020). A comparative study of empirical and analytical fragility functions for the assessment of tsunami building damage in Tumaco, Colombia. *Coastal Eng.* J. 62 (3), 362–372. doi:10.1080/21664250.2020.1726558

Pagnoni, G., Armigliato, A., and Tinti, S. (2019). *Estimation of human damage and economic loss of buildings related to tsunami inundation in the city of Augusta, Italy. Special publications*, 501. London: Geological Society.

Pampell-Manis, A., Horillo, J., Shigihara, Y., and Parambath, L. (2016). Probabilistic assessment of landslide tsunami hazard for the northern Gulf of Mexico. *J. Geophys. Res. Oceans* 121 (1), 1009–1027. doi:10.1007/978-3-642-13296-4_8

Paris, R., Falvard, S., Chagué, C., Goff, J., Etienne, S., and Doumalin, P. (2020). Sedimentary fabric characterized by X-ray tomography: a case-study from tsunami deposits on the Marquesas Islands, French Polynesia. *Sedimentology* 67, 1207–1229. doi:10.1111/sed.12582

Paris, R. (2015). Source mechanisms of volcanic tsunamis. *Phil. Trans. R. Soc. A.* 373, 20140380. doi:10.1098/rsta.2014.0380

Paris, R., Switzer, A. D., Belousova, M., Belousov, A., Ontowirjo, B., Whelley, P. L., et al. (2014). Volcanic tsunami: a review of source mechanisms, past events and hazards in Southeast Asia (Indonesia, Philippines, Papua New Guinea). *Nat. Hazards* 70, 447–470. doi:10.1007/s11069-013-0822-8

Paris, R., Ulvrova, M., Selva, J., Brizuela, B., Costa, A., Grezio, A., et al. (2019). Probabilistic hazard analysis for tsunamis generated by subaqueous volcanic explosions in the Campi Flegrei caldera, Italy. *J. Volcanol. Geothermal Res.* 379, 106–116. doi:10.1016/j.yvolgeores.2019.05.010

Park, H., Alam, M. S., Cox, D. T., Barbosa, A. R., and van de Lindt, J. W. (2019). Probabilistic seismic and tsunami damage analysis (PSTDA) of the Cascadia Subduction Zone applied to Seaside, Oregon. *Int. J. Disaster Risk Reduct.* 35, 101076. doi:10.1016/j.ijdrr.2019.101076

Park, H., Cox, D. T., Lynett, P. J., Wiebe, D. M., and Shin, S. (2013). Tsunami inundation modeling in constructed environments: a physical and numerical comparison of free-surface elevation, velocity, and momentum flux. *Coastal Eng.* 79, 9–21. doi:10.1016/j.coastaleng.2013.04.002

Park, S., van de Lindt, J. W., Cox, D., Gupta, R., and Aguiniga, F. (2012). Successive earthquake-tsunami analysis to develop collapse fragilities. *J. Earthq. Eng.* 16 (6), 851–863. doi:10.1080/13632469.2012.685209

Pasquet, S., Vilićić, I., and Šepić, J. (2013). A survey of strong high-frequency sea level oscillations along the US East Coast between 2006 and 2011. *Nat. Hazards Earth Syst. Sci.* 13, 473–482. doi:10.5194/nhess-13-473-2013

Pattiaratchi, C. B., and Wijeratne, E. M. S. (2015). Are meteotsunamis an underrated hazard? *Phil. Trans. R. Soc. A.* 373, 20140377. doi:10.1098/rsta. 2014.0377

Pedersen, G. (2008). *Advanced numerical models for simulating tsunami waves and runup. Modeling runup with depth integrated equation models.* Singapore: World Scientific, 3–41.

Pedersen, G. K., Lindstrøm, E., Bertelsen, A. F., Jensen, A., Laskovski, D., and Sælevik, G. (2013). Runup and boundary layers on sloping beaches. *Phys. Fluids.* 25 (1), 012102. doi:10.1063/1.4773327

Pellikka, H., Laurila, T. K., Boman, H., Karjalainen, A., Björkqvist, J.-V., and Kahma, K. K. (2020). Meteotsunami occurrence in the Gulf of Finland over the past century. *Nat. Hazards Earth Syst. Sci.* 20, 2535–2546. doi:10.5194/nhess-20-2535-2020

Petrone, C., Rossetto, T., Baiguera, M., la Barra Bustamante, C. D., and Ioannou, I. (2020). Fragility functions for a reinforced concrete structure subjected to earthquake and tsunami in sequence. *Eng. Struct.* 205, 110120. doi:10.1016/j.engstruct.2019.110120

Petrone, C., Rossetto, T., and Goda, K. (2017). Fragility assessment of a RC structure under tsunami actions via nonlinear static and dynamic analyses. *Eng. Struct.* 136, 36–53. doi:10.1016/j.engstruct.2017.01.013

Philibos, B., and Meltzner, A. J. (2020). Segmentation and supercycles: a catalog of earthquake rupture patterns from the Sumatran Sunda Megathrust and other well-studied faults worldwide. *Quat. Sci. Rev.* 241, 106390. doi:10.1016/j. quascirev.2020.106390

Pitilakis, K., Argyroudis, S., Fotopoulou, S., Karafagka, S., Kakderi, K., and Selva, J. (2019). Application of stress test concepts for port infrastructures against natural hazards. The case of Thessaloniki port in Greece. *Reliab. Eng. Syst. Saf.* 184, 240–257. doi:10.1016/j.ress.2018.07.005

Pittore, M., Wieland, M., and Fleming, K. (2017). Perspectives on global dynamic exposure modelling for geo-risk assessment. *Nat. Hazards* 86 (1), 7–30. doi:10.1007/s11069-016-2437-3

Polese, M., Di Ludovico, M., Gaetani d'Aragona, M., Prota, A., and Manfredi, G. (2020). Regional vulnerability and risk assessment accounting for local building typologies. *Int. J. Disaster Risk Reduct.* 43, 101400. doi:10.1016/j.ijdrr.2019.101400

Polet, J., and Kanamori, H. (2016). "Tsunami earthquakes," in *Encyclopedia of complexity and systems science*. Editor R. A. Meyers (Berlin, Heidelberg: Springer), 1–22. doi:10.1007/978-3-642-27737-5567-2

Power, W., Downes, G., and Stirling, M. (2007). Estimation of tsunami hazard in New Zealand due to south American earthquakes. *Pure Appl. Geophys.* 164, 547–564. doi:10.1007/s00024-006-0166-3

Priest, G. R., Witter, R. C., Zhang, Y. J., Goldfinger, C., Wang, K., and Allan, J. C. (2017). New constraints on coseismic slip during southern Cascadia subduction zone earthquakes over the past 4600 years implied by tsunami deposits and marine turbidites. *Nat. Hazards* 88, 285–313. doi:10.1007/s11069-017-2864-9

Proudman, J. (1929). The effects on the sea of changes in atmospheric pressure. *International* 2, 197–209. doi:10.1111/j.1365-246x.1929.tb05408.x

Puzrin, A. M., Germanovich, L. N., and Friedli, B. (2016). Shear band propagation analysis of submarine slope stability. *Géotechnique* 66 (3), 188–201. doi:10.1680/jgeot.151.lm.002

Ramírez-Herrera, M. T., Kostoglodov, V., and Urrutia-Fucugauchi, J. (2011). Overview of recent coastal tectonic deformation in the Mexican subduction zone. *Pure Appl. Geophys.* 168, 1415–1433. doi:10.1007/s00024-010-0205-y

Reese, S., Cousins, W. J., Power, W. L., Palmer, N. G., Tejakusuma, I. G., and Nugrahadi, S. (2007). Tsunami vulnerability of buildings and people in South Java—field observations after the July 2006 Java tsunami. *Nat. Hazards Earth Syst. Sci.* 7, 573–589. doi:10.5194/nhess-7-573-2007

Rijsingen, E.V., Funiciello, F., Corbi, F., and Lallemand, S. (2019). Rough subducting seafloor reduces interseismic coupling and mega-earthquake occurrence: insights from analogue models. *Geophys. Res. Lett.* 46 (6), 3124–3132. doi:10.1029/2018GL081272

Rikitake, T., and Aida, I. (1988). Tsunami hazard probability in Japan. *Bull. Seismol. Soc. Amer.* 78, 1268–1278.

Rivera, F., Rossetto, T., and Twigg, J.(2020). An interdisciplinary study of the seismic exposure dynamics of Santiago de Chile. *Int. J. Disaster Risk Reduct.* 48, 101581. doi:10.1016/j.ijdrr.2020.101581

Rondon, L., Pouliquen, O., and Aussillous, P. (2011). Granular collapse in a fluid: role of the initial volume fraction. *Phys. Fluids* 23 (7), 073301. doi:10.1063/1.3594200

Rong, Y., Jackson, D. D., Magistrale, H., and Goldfinger, C. (2014). Magnitude limits of subduction zone earthquakes. *Bull. Seismol. Soc. Amer.* 104, 2359–2377. doi:10.1785/0120130287

Rosi, M., Di Traglia, F., Pistolesi, M., Esposti Ongaro, T., de' Michieli Vitturi, M., and Bonadonna, C. (2018). Dynamics of shallow hydrothermal eruptions: new insights from Vulcano's Breccia di Commenda eruption. *Bull. Volcanol.* 80, 83. doi:10.1007/s00445-018-1252-y

Rossetto, T., De la Barra, C., Petrone, C., De la Llera, J. C., Vásquez, J., and Baiguera, M. (2019). Comparative assessment of nonlinear static and dynamic methods for analysing building response under sequential earthquake and tsunami. *Earthq. Eng. Struct. Dyn.* 48, 867–887. doi:10.1002/eqe.3167

Rossetto, T., Petrone, C., Eames, I., De la Barra, C., Foster, A., and Macabuag, J. (2018). *Advances in the assessment of buildings subjected to earthquakes and tsunami. In European Conference on Earthquake Engineering Thessaloniki Greece* (Cham, Switzerland: Springer), 545–562.

Ryan, K. J., Geist, E. L., Barall, M., and Oglesby, D. D. (2015). Dynamic models of an earthquake and tsunami offshore Ventura, California. *Geophys. Res. Lett.* 42, 6599–6606. doi:10.1002/2015GL064507

Šepić, J., Vilićić, I., and Strelec Mahović, N. (2012). Northern Adriatic meteorological tsunamis: observations, link to the atmosphere, and predictability. *J. Geophys. Res.* 117, C02002. doi:10.1029/2011JC007608
Saja, A. M. A., Goonetilleke, A., Teo, M., and Ziyath, A. M. (2019). A critical review of social resilience assessment frameworks in disaster management. *Int. J. Disaster Risk Reduct.* 35, 101096. doi:10.1016/j.ijdrr.2019.101096

Salamon, A., and Di Manna, P. (2019). Empirical constraints on magnitude-distance relationships for seismically-induced submarine tsunamigenic landslides. *Earth-Sci. Rev.* 191, 66–92. doi:10.1016/j.earscirev.2019.02.001

Salgado-Gálvez, M. A., Zuloaga, D., Henao, S., Bernal, G. A., and Cardona, O. D. (2018). Probabilistic assessment of annual repair rates in pipelines and of direct economic losses in water and sewage networks. Application to Manizales, Colombia. *Nat. Hazards* 93 (1), S5–S24. doi:10.1007/s11069-017-2987-z

Salgado-Gálvez, M. A., Zuloaga Romero, D., Velásquez, C. A., Carreño, M. L., Cardona, O.-D., and Barbat, A. H. (2016). Urban seismic risk index for Medellín, Colombia, based on probabilistic loss and casualties estimations. *Nat. Hazards* 80 (3), 1995–2021. doi:10.1007/s11069-015-2056-4

Salgado-Gálvez, M. A., Zuloaga-Romero, D., Bernal, G. A., Mora, M. G., and Cardona, O.-D. (2014). Fully probabilistic seismic risk assessment considering local site effects for the portfolio of buildings in Medellín, Colombia. *Bull. Earthq. Eng.* 12 (2), 671–695. doi:10.1007/s10518-013-9550-4

Sallarès, V., and Ranero, C. R. (2019). Upper-plate rigidity determines depth-varying rupture behaviour of megathrust earthquakes. *Nature* 576, 96–101. doi:10.1038/s41586-019-1784-0

Salmanidou, D. M., Guillas, S., Georgiopoulou, A., and Dias, F. (2017). Statistical emulation of landslide-induced tsunamis at the rockall bank, NE atlantic. *Proc. R. Soc. A.* 473 (2200), 20170026. doi:10.1098/rspa.2017.0026

Sandri, L., Tierz, P., Costa, A., and Marzocchi, W. (2018). Probabilistic hazard from pyroclastic density currents in the Neapolitan area (Southern Italy). *J. Geophys. Res. Solid Earth* 123, 3474–3500. doi:10.1002/2017JB014890

Santos, A., and Koshimura, S. (2015). The historical review of the 1755 lisbon tsunami. *J. Geodesy Geomat. Eng.* 1, 38–52. doi:10.17265/2332-8223/2015/04.004

Sarri, A., Guillas, S., and Dias, F. (2012). Statistical emulation of a tsunami model for sensitivity analysis and uncertainty quantification. *Nat. Hazards Earth Syst. Sci.* 12 (6), 2003. doi:10.5194/nhess-12-2003-2012

Savage, S. B., Babaei, M. H., and Dabros, T. (2014). Modeling gravitational collapse of rectangular granular piles in air and water. *Mech. Res. Commun.* 56, 1–10. doi:10.1016/j.mechrescom.2013.11.001

Scala, A., Festa, G., Vilotte, J. P., Lorito, S., and Romano, F. (2019). Wave interaction of reverse-fault rupture with free surface: numerical analysis of the dynamic effects and fault opening induced by symmetry breaking. *J. Geophys. Res. Solid Earth* 124, 1743–1758. doi:10.1029/2018JB016512

Scala, A., Lorito, S., Romano, F., Murphy, S., Selva, J., Basili, R., et al. (2020). Effect of shallow slip amplification uncertainty on probabilistic tsunami hazard analysis in subduction zones: use of long-term balanced stochastic slip models. *Pure Appl. Geophys.* 177, 1497–1520. doi:10.1007/s00024-019-02260-x

Schmidtlein, M. C., Deutsch, R. C., Piegorsch, W. W., and Cutter, S. L. (2008). A sensitivity analysis of the social vulnerability index. *Int. J.* 28 (4), 1099–1114. doi:10.1111/j.1539-6924.2008.01072.x

Selva, J., Acocella, V., Bisson, M., Caliro, S., Costa, A., Della Seta, M., et al. (2019). Multiple natural hazards at volcanic islands: a review for the Ischia volcano (Italy). *J. Appl. Volcanol.* 8 (1), 5. doi:10.1186/s13617-019-0086-4

Selva, J., Bonadonna, C., Branca, S., De Astis, G., Gambino, S., Paonita, A., et al. (2020). Multiple hazards and paths to eruptions: a review of the volcanic system of Vulcano (Aeolian Islands, Italy). *Earth-Sci. Rev.* 207, 103186. doi:10.1016/j.earscirev.2020.103186

Selva, J., Tonini, R., Molinari, I., Tiberti, M. M., Romano, F., Grezio, A., et al. (2016). Quantification of source uncertainties in seismic probabilistic tsunami hazard analysis (SPTHA). *Geophys. J. Int.* 205, 1780–1803. doi:10.1093/gji/ggw107

Sepúlveda, I., Liu, P. L.-F., Grigoriu, M., and Pritchard, M. (2017). Tsunami hazard assessments with consideration of uncertain earthquake slip distribution and location. *J. Geophys. Res. Solid Earth* 122, 7252–7271. doi:10.1002/2017JB014430

Shi, F., Kirby, J. T., Harris, J. C., Geiman, J. D., and Grilli, S. T. (2012). A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundation. *Ocean Model.* 43, 36–51. doi:10.1016/j.ocemod.2011.12.004

Si, P., Shi, H., and Yu, X. (2018b). A general numerical model for surface waves generated by granular material intruding into a water body. *Coastal Eng.* 142, 42–51. doi:10.1016/j.coastaleng.2018.09.001

Si, P., Shi, H., and Yu, X. (2018a). Development of a mathematical model for submarine granular flows. *Phys. Fluids* 30, 083302. doi:10.1063/1.5030349

Sibley, A. M., Cox, D., and Tappin, D. R. (2020). Convective rear-flank downdraft as driver for meteotsunami along English Channel and North Sea coasts 28-29 May 2017. *Nat. Hazards* 14, 114. doi:10.1007/s11069-020-04328-7

Silva, V., Amo-Oduro, D., Calderon, A., Dabbeek, J., Despotaki, V., Martins, L., et al. (2018a). *Global Earthquake Model (GEM) Seismic Risk Map* (version 2018.1), Available at: https://maps.openquake.org/map/global-seismic-risk-map (Accessed November 10, 2020).

Silva, V., Yepes-Estrada, C., Dabbeek, J., Martins, L., and Brzev, S. (2018b). GED4ALL-Global exposure database for multi-hazard risk analysis-multi-hazard exposure taxonomy. GEM Technical Report 2018-01. Pavia: GEM Foundation.

Skarlatoudis, A. A., Somerville, P. G., and Thio, H. K. (2016). Source-scaling relations of interface subduction earthquakes for strong ground motion and tsunami simulation. *Bull. Seismol. Soc. Amer.* 106, 1652–1662. doi:10.1785/0120150320

Smit, A., Kijko, A., and Stein, A. (2017). Probabilistic tsunami hazard assessment from incomplete and uncertain historical catalogues with application to tsunamigenic regions in the pacific ocean. *Pure Appl. Geophys.* 174, 3065–3081. doi:10.1007/s00024-017-1564-4

Snelling, B., Neethling, S., Horsburgh, K., Collins, G., and Piggott, M. (2020). Uncertainty quantification of landslide generated waves using Gaussian process emulation and variance-based sensitivity analysis. *Water* 12 (2), 416. doi:10.3390/w12020416

So, E. K. M., Pomonis, A., Below, R., Cardona, O., King, A., Zulfikar, C., et al. (2012). “An introduction to the global earthquake consequences database (GEMECD),” in *Proceedings of the 15th world conference on earthquake engineering*. Lisbon: Springer.

Song, J., De Risi, R., and Goda, K. (2017). Influence of flow velocity on tsunami loss estimation. *Geosciences* 7 (4), 114. doi:10.3390/geosciences7040114

Sorensen, M. B., Spada, M., Babeyko, A., Wiemer, S., and Grünthal, G. (2012). Probabilistic tsunami hazard in the mediterranean sea. *J. Geophys. Res.* 117 (B1), 13. doi:10.1029/2010JB008169

Spielman, S. E., Tuccillo, J., Folch, D. C., Schweikert, A., Davies, R., Wood, N., et al. (2020). Evaluating social vulnerability indicators: criteria and their application to the Social Vulnerability Index. *Nat. Hazards* 100 (1), 417–436. doi:10.1007/s11069-019-03820-z

Stolle, J., Takabatake, T., Hamano, G., Ishii, H., Iimura, K., Shibayama, T., et al. (2019). Debris transport over a sloped surface in tsunami-like flow conditions. *Coastal Eng.* J. 61 (2), 241–255. doi:10.1080/21664250.2019.1586288

Strasser, F. O., Arango, M. C., and Bommer, J. J. (2010). Scaling of the source dimensions of interface and intraslab subduction-zone earthquakes with moment magnitude. *Seismol. Res. Lett.* 81, 941–950. doi:10.1785/gssrl81.6.9410.1785/gssrl81.6.941

Strunz, G., Post, J., Zosseder, K., Wegscheider, S., Mück, M., Riedlinger, T., et al. (2011). Tsunami risk assessment in Indonesia. *Nat. Hazards Earth Syst. Sci.* 11 (1), 67–82. doi:10.5194/nhess-11-67-2011

Stucchi, M., Albini, P., Mirto, M., and Rebez, A. (2004). Assessing the completeness of Italian historical earthquake data. *Ann. Geophys.* 47, 35. doi:10.4401/ag-3330

Sugawara, D., Goto, K., and Jaffe, B. E. (2014). Numerical models of tsunami sediment transport—current understanding and future directions. *Mar. Geol.* 352, 295–320. doi:10.1016/j.margeo.2014.02.007

Suppasri, A., Charvet, I., Imai, K., and Imamura, F. (2015). Fragility curves based on data from the 2011 tohoku-oki tsunami in ishinomaki city, with discussion of parameters influencing building damage. *Earthq. Spectra* 31, 841–868. doi:10.1193/053013EQS138M

Suppasri, A., Hasegawa, N., Makinoshima, F., Imamura, F., Latcharote, P., and Day, S. (2016). An analysis of fatality ratios and the factors that affected human fatalities in the 2011 Great East Japan tsunami. *Front. Built Environ.* 2, 32. doi:10.3389/fbuil.2016.00032

Synolakis, C. E., Bernard, E. N., Titov, V. V., Kânoğlu, U., and González, F. I. (2008). Validation and verification of tsunami numerical models. *Pure Appl. Geophys.* 165, 2197–2228. doi:10.1007/s00024-004-0427-y
Šepić, J., Vilićić, I., and Monserrat, S. (2016). Quantifying the probability of meteotsunami occurrence from synoptic atmospheric patterns. *Geophys. Res. Lett.* 43, 10377–10384. doi:10.1002/2016GL070754

Tappin, D. R., Watts, P., and Grilli, S. T. (2008). The Papua New Guinea tsunami of 17 July 1998: anatomy of a catastrophic event. *Nat. Hazards Earth Syst. Sci.* 8 (2), 243–266. doi:10.5194/nhess-8-243-2008

ten Brink, U., Geist, E. L., and Andrews, B. D. (2006). Size distribution of submarine landslides and its implication to tsunami hazard in Puerto Rico. *Geophys. Res. Lett.* 33 (11), 23. doi:10.1142/9789812774613_0004

Thingbaijam, K. K. S., Martin Mai, P., and Goda, K. (2017). New empirical earthquake source-scaling laws. *Bull. Seismol. Soc. Am.* 107, 2225–2246. doi:10.1785/0120170017

Thio, H. K., Somerville, P., and Ichinose, G. (2007). Probabilistic analysis of strong ground motion and tsunami hazards in southeast asia. *J. Earthq. Tsunami* 01, 119–137. doi:10.1142/S1793431107000080

Tiberti, M. M., Basili, R., and Vannoli, P. (2014). Ups and downs in western Crete (Hellenic subduction zone). *Sci. Rep.* 4, 5677. doi:10.1038/srep05677

Tinti, S. (1991). Assessment of tsunami hazard in the Italian seas. *Nat. Hazards* 4, 267–283. doi:10.1007/bf00162792

Tinti, S., Bortolucci, E., and Armigliato, A. (1999). Numerical simulation of the landslide-induced tsunami of 1988 on Vulcano Island, Italy. *Bull. Volcanol.* 61 (1), 121–137. doi:10.1007/s004450050267

Tinti, S., Bortolucci, E., and Vannini, C. (1997). A block-based theoretical model suited to gravitational sliding. *Nat. Hazards* 16, 1–28. doi:10.1023/A10.1023/a:1007934804464

Tinti, S., Pagnoni, G., and Zaniboni, F. (2006). The landslides and tsunamis of the 30th of December 2002 in Stromboli analysed through numerical simulations. *Bull. Volcanol.* 68, 462–479. doi:10.1007/s00445-005-0022-9

Tonini, R., Basili, R., Maesano, F. E., Tiberti, M. M., Lorito, S., Romano, F., et al. (2020). Importance of earthquake rupture geometry on tsunami modelling: the Calabrian Arc subduction interface (Italy) case study. *Geophys. J. Int.* 223, 1805–1819. doi:10.1093/gji/gga409

Trapper, P. A., Puzrin, A. M., and Germanovich, L. N. (2015). Effects of shear band propagation on early waves generated by initial breakoff of tsunamigenic landslides. *Mar. Geol.* 370, 99–112. doi:10.1016/j.margeo.2015.10.014

Triantafyllou, I., Novikova, T., Charalampakis, M., Fokaefs, A., and Papadopoulos, G. A. (2019). Quantitative tsunami risk assessment in terms of building replacement cost based on tsunami modelling and GIS methods: the case of crete isl., hellenic arc. *Pure Appl. Geophys.* 176 (7), 3207–3225. doi:10.1007/s00024-018-1984-9

Tufekci, D., Suzen, M. L., Yalciner, A. C., and Zaytsev, A. (2018). Revised MeTHuVA method for assessment of tsunami human vulnerability of Bakirkoy district, Istanbul. *Nat. Hazards* 90 (2), 943–974. doi:10.1007/s11069-017-3082-1

Ulrich, T., Gabriel, A.-A., Ampuero, J.-P., and Xu, W. (2019a). Dynamic viability of the 2016 Mw 7.8 Kaikoura earthquake cascade on weak crustal faults. *Nat. Commun.* 10 (1), 1213. doi:10.1038/s41467-019-09125-w

Ulrich, T., Vater, S., Madden, E. H., Behrens, J., van Dinther, Y., van Zelst, I., et al. (2019b). Coupled, physics-based modeling reveals earthquake displacements are critical to the 2018 palu, sulawesi tsunami. *Pure Appl. Geophys.* 176, 4069–4109. doi:10.1007/s00024-019-02290-5

Ulvrova, M., Paris, R., Nomikou, P., Kelfoun, K., Leibrandt, S., Tappin, D. R., et al. (2016). Source of the tsunami generated by the 1650 AD eruption of Kolumbo submarine volcano (Aegean Sea, Greece). *J. Volcanol. Geothermal Res.* 321, 125–139. doi:10.1016/j.jvolgeores.2016.04.034

UNISDR. (2013). *United Nations international Office for disaster risk reduction. Global assessment report on disaster risk reduction 2013*. Geneva: Springer.

UNISDR. (2015). *United Nations international Office for disaster risk reduction. Global assessment report on disaster risk reduction 2015*. Geneva: Switzerland.

UNISDR. (2007). United Nations international Office for disaster risk reduction. Hyogo framework for action 2005–2015: building the resilience of Nations and communities to disasters. Geneva. Switzerland. Available at: https://www.unisdr.org/files/1037_hyogoframeworkforactionenglish.pdf (Accessed October 22, 2020).

Urgeles, R., and Camerlenghi, A. (2013). Submarine landslides of the Mediterranean Sea: trigger mechanisms, dynamics, and frequency-

magnitude distribution. *J. Geophys. Res. Earth Surf.* 118 (4), 2600–2618. doi:10.1002/2013jf002720

Urlaub, M., Talling, P. J., and Masson, D. G. (2013). Timing and frequency of large submarine landslides: implications for understanding triggers and future geohazard. *Quat. Sci. Rev.* 72, 63–82. doi:10.1016/j.quascirev.2013.04.020

van de Lindt, J. W., and Dao, T. N. (2009). Performance-based wind engineering for wood-frame buildings. *J. Struct. Eng.* 135 (2), 169–177. doi:10.1061/(asce)0733-9445(2009)135:2(169)

Vela, J., Pérez, B., González, M., Otero, L., Olabarrieta, M., Canals, M., et al. (2014). Tsunami resonance in palma bay and harbor, majorca island, as induced by the 2003 western mediterranean earthquake. *J. Geol.* 122 (2), 165–182. doi:10.1086/675256

Vilićić, I. (2008). Numerical simulations of the Proudman resonance. *Continental Shelf Res.* 28 (4), 574–581. doi:10.1016/j.csr.2007.11.005

Vilićić, I., Horvath, K., Strelec Mahović, N., Monserrat, S., Marcos, M., Amores, A., et al. (2014). Atmospheric processes responsible for generation of the 2008 Boothbay meteotsunami. *Nat. Hazards* 74, 25–53. doi:10.1007/s11069-013-0811-y

Vilićić, I., and Šepić, J. (2009). Destructive meteotsunamis along the eastern adriatic coast: overview. Physics and chemistry of the Earth. *Parts A/B/C* 34 (17–18), 904–917. doi:10.1016/j.pce.2009.08.004

Volpe, M., Lorito, S., Selva, J., Tonini, R., Romano, F., and Brizuela, B. (2019). From regional to local SPTHA: efficient computation of probabilistic tsunami inundation maps addressing near-field sources. *Nat. Hazards Earth Syst. Sci.* 19, 455–469. doi:10.5194/nhess-19-455-2019

Walter, T. R., Haghshenas Haghighi, M., Schneider, F. M., Coppola, D., Motagh, M., Saul, J., et al. (2019). Complex hazard cascade culminating in the Anak Krakatau sector collapse. *Nat. Commun.* 10, 4339. doi:10.1038/s41467-019-12284-5

Wang, K., Hu, Y., and He, J. (2012). Deformation cycles of subduction earthquakes in a viscoelastic Earth. *Nature* 484, 327–332. doi:10.1038/nature11032

Wang, L., Zhang, X., and Tinti, S. (2019). Shallow landslides modeling using a particle finite element model with emphasis on landslide evolution. 11, 21. doi:10.5194/esurf-2019-17

Ward, S. N. (2001). Landslide tsunami. *J. Geophys. Res.* 106 (B6), 11201–11215. doi:10.1029/2000jb900450

Watts, P., Grilli, S. T., Kirby, J. T., Fryer, G. J., and Tappin, D. R. (2003). Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model. *Nat. Hazards Earth Syst. Sci.* 3, 391–402. doi:10.5194/nhess-3-391-2003

Watts, P. (2000). Tsunami features of solid block underwater landslides. *J. Waterway Port Coastal Ocean Eng.* 126 (3), 144–152. doi:10.1061/(asce)0733-950x(2000)126:3(144)

Watts, P., and Waythomas, C. F. (2003). Theoretical analysis of tsunami generation by pyroclastic flows. *J. Geophys. Res.* 108 (B12), 2563. doi:10.1029/2002JB002265

Wei, Y., Bernard, E. N., Tang, L., Weiss, R., Titov, V. V., Moore, C., et al. (2008). Real-time experimental forecast of the Peruvian tsunami of August 2007 for U.S. coastlines. *Geophys. Res. Lett.* 35 (4), 67. doi:10.1029/2007GL032250

Wells, D. L., and Coppersmith, K. J. (1994). New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. *Bull. Seismol. Soc. Am.* 84, 974–1002.

Wendt, J., Oglesby, D. D., and Geist, E. L. (2009). Tsunamis and splay fault dynamics. *Geophys. Res. Lett.* 36 (15), 34. doi:10.1029/2009GL038295

Whitmore, P., and Knight, B. (2014). “Meteotsunami forecasting: sensitivities demonstrated by the 2008 Boothbay, Maine, event,” in *Meteorological Tsunamis, U. S. The East Coast, and Other Coastal Regions*. Editors Vilićić, I., Monserrat, S., and Rabinovich, A. B. (Cham, Switzerland: Springer), doi:10.1007/978-3-319-12712-5_2

Williamson, A. L., Rim, D., Adams, L. M., LeVeque, R. J., Melgar, D., and González, F. I. (2020). A source clustering approach for efficient inundation modeling and regional scale probabilistic tsunami hazard assessment. *Front. Earth Sci.* 8, 442. doi:10.3389/feart.2020.591663

Woodruff, I., Kirby, J., Shi, F., and Grilli, S. (2018). Estimating meteo-tsunami occurrences for the US East coast. *Int. Conf. Coastal. Eng.* 1 (36) 66. doi:10.9753/icce.v36.currents.66
Yamao, S., Esteban, M., Yun, N. Y., Mikami, T., and Shibayama, T. (2015). "Estimation of the current risk to human damage life posed by future tsunamis in Japan," in *Handbook of coastal disaster mitigation for engineers and planners*. Editor Esteban, M. (Boston: Butterworth-Heinemann), 257-275.

Yamashita, K., Suppasri, A., Oishi, Y., and Imamura, F. (2018). Development of a tsunami inundation analysis model for urban areas using a porous body model. *Geosciences* 8 (1), 12. doi:10.3390/geosciences8010012

Yavari-Ramshe, S., and Ataie-Ashtiani, B. (2016). Numerical modeling of subaerial and submarine landslide-generated tsunami waves-recent advances and future challenges. *Landslides* 13 (6), 1325-1368. doi:10.1007/s10346-016-0734-2

Yokoyama, I. (1987). A scenario of the 1883 Krakatau tsunami. *J. Volcanol. Geothermal Res.* 34, 123–132. doi:10.1016/0377-0273(87)90097-7

Zechar, J. D., and Frankel, K. L. (2009). Incorporating and reporting uncertainties in fault slip rates. *J. Geophys. Res.* 114, 121. doi:10.1029/2009JB006325

Zengaffinen, T., Løvholt, F., Pedersen, G. K., and Harbitz, C. B. (2020). *Effects of rotational submarine slump dynamics on tsunami genesis: new insight from idealized models and the 1929 Grand Banks event Special publications*. London: Geological Society.

**Conflict of Interest:** HT was employed by AECOM.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Behrens, Løvholt, Jalayer, Lorito, Salgado-Gálvez, Sørensen, Abadie, Aguirre-Ayerbe, Aniel-Quiroga, Babeyko, Baiguera, Basili, Belliazzi, Grezio, Johnson, Murphy, Paris, Rafliana, De Risi, Rossetto, Selva, Taroni, Del Zoppo, Armigliato, Bureš, Cech, Cecioni, Christodoulides, Davies, Dias, Bayraktar, González, Gritsevich, Guillas, Harbitz, Kânoğlu, Macías, Papadopoulos, Polet, Romano, Salamon, Scala, Stepinac, Tappin, Thio, Tonini, Triantafyllou, Ulrich, Varini, Volpe and Vyhmeister. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.