Tsunami risk management for crustal earthquakes and non-seismic sources in Italy

Destructive tsunamis are most often generated by large earthquakes occurring at subduction interfaces, but also other “atypical” sources—defined as crustal earthquakes and non-seismic sources altogether—may cause significant tsunami threats. Tsunamis may indeed be generated by different sources, such as earthquakes, submarine or coastal landslides, volcano-related phenomena, and atmospheric perturbations. The consideration of atypical sources is important worldwide, but it is especially prominent in complex tectonic settings such as the Mediterranean, the Caribbean, or the Indonesian archipelago. The recent disasters in Indonesia in 2018, caused by the Palu-Sulawesi magnitude Mw 7.5 crustal earthquake and by the collapse of the Anak-Krakatau volcano, recall the importance of such sources. Dealing with atypical sources represents a scientific, technical, and computational challenge, which depends on the capability of quantifying and managing uncertainty efficiently and of reducing it with accurate physical modelling. Here, we first introduce the general framework in which tsunami threats are treated, and then we review the current status and the expected future development of tsunami hazard quantifications and of the tsunami warning systems in Italy, with a specific focus on the treatment of atypical sources. In Italy, where the memory of historical atypical events like the 1908 Messina earthquake or the relatively recent 2002 Stromboli tsunami is still vivid, specific attention has been indeed dedicated to the progressive development of innovative strategies to deal with such atypical sources. More specifically, we review the (national) hazard analyses and their application for coastal planning, as well as the two operating tsunami warning systems: the national warning system for seismically generated tsunamis (SiAM), whose upstream component—the CAT-INGV—is also a Tsunami Service Provider of the North-eastern Atlantic, the Mediterranean and connected seas Tsunami Warning System (NEAMTWS) coordinated by the Intergovernmental Coordination Group established by the Intergovernmental Oceanographic Commission (IOC) of UNESCO, and the local warning system for tsunamis generated by volcanic slides along the Sciara del Fuoco of Stromboli volcano. Finally, we review the state of knowledge about other potential tsunami sources that may generate significant tsunamis for the Italian coasts, but that are not presently considered in existing tsunami warning systems. This may be considered the first step towards their inclusion in the national tsunami hazard and warning programs.

del Fuoco of Stromboli volcano. Finally, we review the state of knowledge about other potential tsunami sources that may generate significant tsunamis for the Italian coasts, but that are not presently considered in existing tsunami warning systems. This may be considered the first step towards their inclusion in the national tsunami hazard and warning programs.

Introduction
A tsunami is a series of gravity waves propagating in the water and generated by its sudden large-scale displacement due to one external source [1]. Tsunami waves are usually characterized by characteristic wavelengths longer than the sea depth. Most of the giant basin-wide tsunamis that are well known to the general public, such as the recent 2004 Indian ocean and the 2011 Tohoku tsunamis (e.g., [2]), are mostly caused by very large megathrust earthquakes (M > 8.0) occurring at subduction interfaces, which is where elastic energy is continuously stored and suddenly released by seismic events because of the friction between the incoming oceanic crust and the overlying continental crust, forced by the convergent motion between the two (e.g., [3]).
Some recent tsunamis, such as the Palu-Sulawesi and Sunda Strait events occurred in 2018, raised the attention toward tsunami sources different from the megathrust earthquakes occurring on the subduction interface. Indeed, the Palu-Sulawesi tsunami was generated by a magnitude Mw 7.5 crustal earthquake with a strike-slip source mechanism (that is, an earthquake occurring on a nearly vertical fault where two adjacent blocks move mostly horizontally with respect to each other) and with the possible contribution of landslides mobilized by the seismic shaking inside the Palu bay [4][5][6]; the Sunda Strait tsunami was instead generated by the collapse of the volcano edifice of Anak Krakatau, a volcanic island between Sumatra and Java islands [7][8][9][10]. Indeed, the initial water displacement developing in tsunami may be due to several causes, including the coseismic displacement caused by many kind of earthquakes, including non-subduction-interface earthquakes occurring inside or outside subduction zones, as well as by other causes such as landslides, volcanic phenomena, meteorological perturbations, or even rarer events like the impact of asteroids [11]. As usually not covered by tsunami warning systems (TWSs), such sources were recently defined as "atypical" by the IOC/UNESCO TOWS-WG Inter-ICG TTTWO (Task Team on Tsunami Watch Operations of the Working Group on Tsunamis and Other Hazards Related to Sea-Level Warning and Mitigation Systems connecting all Intergovernmental Coordination Group on Tsunami Warning and Mitigation Systems) [12]).
Atypical sources, as said, include both seismic and non-seismic sources. Earthquakes are generated by the stress build-up determined by tectonic plate movements that induces sudden displacement on faults affecting the Earth's brittle crust. Fault orientation in space and reciprocal movement of the two sides depends mainly on the orientation and sign of the stress field (compressional and dilatational), on the characteristics of the medium (e.g., elastic properties, rheology, friction), as well as on local geological features (e.g., pre-existing discontinuities) [1,13]. In subduction zones, large earthquakes tend to occur in the well-known interface that separates the two tectonic plates and that plunges with a relatively low angle into the Earth's mantle, typically with thrust earthquakes (that is, with contractional movement on the fault rupture). Crustal seismicity occurs instead within the crust, at plate boundaries and even in plate interiors. Within the crust, different earthquake mechanisms may coexist even within small areas. Atypical sources include all these crustal earthquakes capable of generating tsunamis outside the subduction zones, which are not always equally well mapped [14]. Noteworthy, the definition of atypical sources also includes the subduction-zone earthquakes generated outside the subduction interface, such as splay and back-thrust faults branching from the interface in the upper plate and accretionary wedge, and outer-rise earthquakes, mostly occurring seaward from the trench on the oceanic crust as a result of stress redistribution following subduction earthquakes. Among all tsunamigenic earthquakes, the "tsunami earthquakes" deserve particular attention since they generate larger tsunamis than expected given their seismic moment magnitude, and both their mechanical process and frequency of occurrence are not completely understood [15]. Atypical sources also include all the other non-seismic tsunamigenic sources like sub-aerial near-coast and submarine landslides, a wide range of volcanic phenomena occurring both during eruptions or in non-eruptive periods (pyroclastic flows-fast currents of hot gas and volcanic materials, submarine explosions, caldera collapses, and flank instability are the most important ones), and meteotsunamis-tsunamis generated by atmospheric pressure disturbances (e.g., [16][17][18][19][20][21][22]).
In complex tectonic settings and small basins, such as, for example, the Mediterranean region, the Caribbean [23,24], or around the Indonesian and Philippine archipelago north and east of the Sunda trench [19,25], atypical sources may even represent the majority of potential tsunamigenic sources. For example, in the Mediterranean Sea, numerous historical significant tsunamis were caused by non-subduction earthquakes or non-seismic sources [14,[26][27][28][29][30][31][32][33][34][35], and such sources generated the majority of the events causing damaging tsunamis in the available historical record [36]. All the largest recent tsunamis in the Mediterranean were caused by crustal seismicity: the 2003 M6.8 Zemmouri-Boumerdes, the 2017 M6.8 Kos-Bodrum, and the recent 2020 M7.0 Samos-Izmir earthquakes, all generating tsunami whose maximum run-up was larger than 1 m [37][38][39]. Noteworthy, also non-seismic or mixed seismic/non-seismic sources generated in the recent past significant tsunamis, like the one generated at Stromboli in 2002 [40][41][42][43] as well as the most destructive tsunami in the record, related to the 1908 Mw 7.1 Messina Straits earthquake [44]. This non-subduction earthquake caused at least 80,000 deaths, and a few hundred were added by the tsunami [45], whose seismic generation was possibly enhanced by an earthquake-induced landslide [46][47][48][49][50]. More in general, hazard disaggregation results [31,[51][52][53], which account for both source probability and tsunami potential and related uncertainty, show that the contribution to the total tsunami hazard of crustal sources is significant also along the coasts exposed to subduction tsunamis (e.g. eastern Sicily), and they may even represent the dominant contribution in many areas in the Mediterranean.
Tsunami forecast for Tsunami Warning Systems (TWSs) and Probabilistic Tsunami Hazard Analyses (PTHAs) are the main scientific inputs to the tsunami risk management [11,54,[54][55][56][57][58][59]. For TWSs, it is required to timely produce a forecast of the tsunami intensity at different sites due to an event (for example, an earthquake) that has just occurred, or it is forecasted to occur in the relatively near future (for example, for a volcanic eruption or a subaerial landslide) based on the information registered by a monitoring system. The forecast is then generally used to define alert levels associated with specific short-term risk reduction measures, like population evacuation [55,[60][61][62][63]. On the other hand, PTHA is a forecast of the potential future activity and it quantifies the probability of exceeding different tsunami intensities due to whatever possible tsunami in a target location within a relatively long time windows (e.g., 50 years). PTHA is based on the past behaviour of the system, and it is the main scientific input to long-term tsunami risk mitigation programs, like the definition of building codes or the design of evacuation areas themselves for TWSs [58,[60][61][62].
Tsunami forecasts should be in principle based on the definition of a set of individual sources sampling the whole natural variability of relevant physical generation mechanisms, the quantification of the probability of occurrence of each source, and the evaluation of the potential consequent tsunamis. Tsunami forecast typically is based on the numerical simulation of the tsunami generation and propagation for each individual source, usually considering non-linear models in the shallow water approximation, to model the tsunami impact on the coast [11,54,64,65]. This procedure makes these assessments potentially very challenging from both the modelling and the computational points of view, especially if atypical sources have to be taken into account.
The computational challenge arises from the fact that typically large source variability has to be modelled, implying the requirement of running a very large number of individual simulations. Indeed, non-subduction seismicity takes place over large source areas, potentially covering the entire sea-floor and the near-coast areas and, as discussed above, includes a broad variety of source mechanisms, that is, a broad variety of fault orientations and fault movements. In comparison, large subduction zone earthquake behaviour is indeed more predictable, as it is largely dominated by thrust events on the plate interface. To model this variability a very large set of scenarios is needed (e.g., > 10 7 , considering only seismic sources in a small basin like the Mediterranean, see [31]) and thus their integration into tsunami hazard and warning depends on the development of specific methodological and computational strategies [31,54,[66][67][68][69].
In addition, the identification and characterization of atypical sources is also difficult. Atypical earthquake faults and physics are usually less known than subduction sources, especially offshore, making it difficult to fully characterize their geometry and dynamics [14,31,70,71]. The characterization of the non-seismic sources is even more challenging due to their specific local geophysical (e.g. tectonic and geological) and geomechanical conditions (e.g. soil and rock properties), and the characterization and modelling of these sources and of the associated tsunamigenesis is still a matter of intense scientific research [16-18, 20, 72, 73].
The modelling challenge consists of capturing the fundamentals of the dynamics or the kinematics of the source process (e.g., the generation and the propagation of a submarine landslide or the seafloor deformation induced by an earthquake) and of the solid-fluid coupling that determines the input to tsunami propagation model, as well as the fundamental physics that the propagation model should include.
The earthquake induced deformations are determined by the dynamic rupture, controlled in turn by the friction, the medium properties and the pre-existing stress field (e.g. [13]). The seismic rupture results in a time-dependent slip distribution on the fault, which produces in turn the deformation of the seafloor that, once transmitted to the sea water column above it, induces a sea level perturbation then propagating under the action of gravity until the equilibrium is restored (e.g. [74,75]). The slower the rupture propagation and the local displacement on the fault, and the longer the earthquake duration, the larger are the departures from the widely used approximation of instantaneous sea floor displacement and hydrostatic transmission. The landslide rheology and related kinematics has direct consequence on the tsunami generation [76]. Advanced rheological and numerical models for the landslide dynamics normally distinguish between granular rheology (e.g. [77][78][79][80]) and viscoplastic rheology (e.g. [76,81]). Volcanic sources are especially difficult to model in this sense. Volcanic mass movements include subaerial pyroclastic currents, i.e., rapid granular flows at high (> 400°C) temperature, characterized by broad and dominantly sub-millimetric grain-size spectrum, enhancing the transport of interstitial pressure and consequent friction drop. Their tsunamigenic capability is still a matter of fundamental research but claims for urgent assessment [19,20].
Constraining fluid resistance forces (pressure drag, skin friction, and added mass) may be as important as source properties (e.g. [82,83]). There are indeed still open questions also about the physics of the solid/fluid (earth, ocean, and atmosphere) coupling controlling the tsunami generation mechanisms. This includes, for instance, whether or not horizontal momentum transferred to the water body by sudden seafloor motion plays a substantial role ( [84] and references therein), as well as the role played by water compressibility in generating hydro-acoustic precursors of gravity waves [85], or the generation mechanism of meteotsunamis [86,87]. Other important modelling challenges are related to tsunami shoaling inundation modelling, and thus in common to all tsunami sources. This includes, for example, phenomena that are sometimes observed, such as formation of near shore undular bores (e.g. [88]). More generally, conventional tsunami simulation typically simplifies the inundation modelling process through depth averaging the physics related to friction and boundary layer condition [64]. A recent investigation on this aspect was carried out by Qin et al. [89]. A deep understanding of these effects on tsunami hazards is still lacking.
To efficiently contribute to forecast procedures like TWS and PTHA, all the advancements in modelling techniques should also be accompanied by a synergic parallel effort dedicated to the management of existing uncertainty, both in forecasting the occurrence of the source process and in quantifying the potential consequent tsunami, in order to increase precision without loosing in accuracy. Also this field is indeed a field of intense research activity [11,31,54,63,[90][91][92], and it is particularly challenging when extending to the large variety of atypical tsunami sources (e.g. [11,17,32,54]).
As said, present TWSs primarily deal with large earthquakes occurring on subduction interfaces and were indeed created in response to tsunamis generated by this kind of seismic events [56]. This is often the case also for PTHA, because the subduction earthquakes may be thought to dominate the hazard either at the global (e.g., [93]) or at the local scale at specific places (e.g., [94]). Tsunamis generated by atypical sources are usually more local and probably less frequent than subduction zone events worldwide, and their effective hazard is more difficult to constrain, and thus it is often neglected [11,17]. Also, their forecasting in near-real-time for early warning purposes is rather challenging and not yet a well-established practice, particularly as far as non-seismic sources are concerned.
More specifically, TWSs and PTHA for seismically generated tsunamis are definitely the most advanced ones. Tsunami Service Providers (TSPs) of the IOC/UNESCO tsunami programme cover most of the seas worldwide. In the Northeastern Atlantic, the Mediterranean and connected seas Tsunami Early Warning and Mitigation System (NEAMTWS), to deal with crustal earthquakes and to deliver timely alerts also for local targets, all the TSPs operate based on decision matrices that define the alert levels solely based on the earthquake parameters (magnitude, location, and depth), the source-coast and the source-target distances. This simplifies extremely the forecasting procedure, as it ignores the specific characteristics of the source event and the asymmetric and complex tsunami propagation features while only considering attenuation with distance. For the future, more sophisticated forecasting strategies may be desired, which also explicitly deal with uncertainties (e.g., [95][96][97]).
PTHA for seismically induced tsunamis is, in most cases, limited to megathrust sources or based on the study of a few pre-selected scenarios ( [11,54] and reference therein). In the NEAM region, a recent regional-wide hazard model (NEAMTHM18; [51][52][53]) has been developed, which includes crustal earthquakes making large use of High-Performance Computing (HPC) resources [31,53,66].
Tsunami early warning and hazard assessment from non-seismic sources such as landslides and volcano flank instabilities, which are closely related phenomena, are less developed than their seismic equivalents. Consequently, despite few prototypal examples of PTHA including non-seismic sources exist [11,17,32,98], a standardized approach to model the hazard and manage the warnings for tsunamis caused by such non-seismic sources does not exist (e.g. [11,17,54,65]). As discussed above, this is related both to the very large variability of potential sources, the difficulties in efficiently modelling them, and the difficulty in reasonably constraining recurrences. In addition, with reference to seismically generated tsunamis, the smaller scale of the sources and their complex dynamics require the adoption of more advanced tsunami modelling, including, for example, wave dispersion, increasing the difficulty in developing widely accepted methods. However, it is noted that important steps forward have been developed in very recent years [16-22, 32, 65].
For submarine landslides, which in principle could occur at any margin proximal to coastlines (as indirectly evident, for example, from the Mediterranean landslide database; [99]), there is no regional tsunami early warning system equivalent to those from earthquakes in place worldwide. We are not aware of local early warning systems either. A key challenge is the rapid quantification of a characteristic observable of the landslide, playing a role similar to magnitude in the case of earthquakes. Proposed methodologies are also scarce, although one proposed approach is the use of GPS monitoring, including measuring the flexural rebound of the lithosphere when a very large landslide is released [100]. Moreover, bottom pressure sensors should, in principle, be able to also detect the propagating wave, although not associated with similar source information as for earthquakes. Similar issues remain for regional subaerial landslides. Yet, early warning systems are in place for site-specific local subaerial landslide tsunamis. Examples of such systems can be found in Norwegian fjords (e.g. [101]), and these are implemented in locations where there are significant slope instabilities and ongoing slope motion. These early warning systems in Norway differ significantly from conventional TWSs, as they are based on slope stability monitoring. Typically, this monitoring is extensive and consists of a range of different types of measurements, from remote characterization, microseismicity, to more local measurements using, for instance, strain and aperture motion. The early warning is based on slope acceleration, and the alert system relies on warning the population days ahead of the expected incident. Part of the reason for this is the short travel times (can be less than 5 min) and the huge amplitudes of the tsunamis (can be even larger than 100 m run-up heights). Examples of instrumentation used for the early warning and slope characterization are reported in several publications (e.g. [102][103][104]). The localization of the unstable slopes makes use of extensive and continuous remote sensing techniques and is connected to a rudimentary (semiquantitative) hazard analysis method based on the slope motion [105]. However, a firm link between tsunami hazard levels and evacuation zones are not yet in place (e.g. [17]).
Tsunami early warning in volcanic systems or for meteotsunamis face difficulties similar to the ones already discussed for landslides. Again, the main challenge is the early detection of the volcanic and atmospheric conditions that may favour the occurrence of tsunamigenic events (e.g. [10,73]). The focus on the triggering conditions is important to manage extended source areas, as for meteotsunamis, and to maximize the lead time to inundation, at the cost of increasing the uncertainty of the warning. Indeed, relying only on the tsunami signal drastically reduces the lead time for the warning, challenging the response to the warnings, at least in the near field. For example, in Norway, it has been decided that, because this lead time is too short, TWSs should not be based on the tsunami measurement but should rather be issued well in advance before the slope failure. Obviously, this has implications for possible high numbers of false warnings and evacuations, as evident from experience with the Mannen slope failure in Norway, where the local population has been evacuated many times before the slope failure actually happened. To monitor extended source areas, the alternative approach of direct detection of the tsunami waves is also challenging, since it might need very densely distributed measure points (DART-Deep-ocean Assessment and Reporting of Tsunamis, tide-gauges, elastic beacons, IDSL-Inexpensive Device for Sea Level Monitoring, etc.), with sufficiently high rates of sampling to measure the higher frequencies typically generated by non-seismic sources. Measure points must be densely distributed since an already propagating tsunami may reach inhabited areas in a few minutes; thus, very early detection is required. In this, it is interesting to note that a parallel can be drawn with near field first alert from seismic early warning, which is only based on earthquake detection, and for which both near-field monitoring and automatizations are very important [106][107][108]. For tsunamis, however, even with high quality and spatially dense data, the characterization of such sources from these signals could be excessively expensive and technically challenging too. At local level and for well-defined tsunamigenic sources, a combined approach based on both tsunami detection and source monitoring may instead be adopted, as in the case of the Stromboli TWS, which will be discussed in detail below.
The central position of Italy in the Mediterranean and the evident importance of atypical sources for Italian coasts stimulated the scientific community and decisionmakers to develop innovative procedures that progressively extended tsunami hazard analyses and warning systems to such atypical sources. The Italian national warning system for seismically generated tsunamis (Sistema di Allertamento nazionale per i Maremoti-SiAM, which includes the Italian National Civil Protection -DPC, the Istituto Superiore per la Protezione e la Ricerca Ambientale -ISPRA, and the Istituto Nazionale di Geofisica e Vulcanologia -INGV [109]) deals with all potential seismic sources in the Mediterranean basin, including both subduction and crustal seismicity. The upstream component of the SiAM, the CAT-INGV (Centro Allerta Tsunami of INGV; http://www.ingv.it/cat/it/, [110,111]) is also a TSP (Tsunami Service Provider) of the North-eastern Atlantic, the Mediterranean and connected seas Tsunami Warning System (NEAMTWS), coordinated by the Intergovernmental Coordination Group established by the Intergovernmental Oceanographic Commission of UNESCO. Also, coastal planning for SiAM is defined considering non-subduction zone earthquakes, being based on a regional PTHA including them [31, 51-53, 62, 63]. Another local warning system is dedicated to tsunamis generated by collapses occurring within the volcanic system of Stromboli. This system represents, to our knowledge, the only fully operational tsunami warning system for volcano sources. It has been operational since 2003, and it combines activations related to the early detection of the volcanic conditions that may lead to a tsunami with activations triggered by the early detection of the tsunami wave signal in the very near-field. Two recent small tsunamis, generated in 2019, challenged the system, providing the first indications for its further developments.
Here, we first introduce the general framework in which tsunami threats are treated by briefly reviewing the tsunami forecasting methods and the tsunami generation and propagation modelling techniques. Then we present the tectonic settings for the Mediterranean region in some details (Sect. 2). Next, we discuss the state of the art and the planned developments of tsunami monitoring and tsunami warning systems in Italy (Sects. 3 and 4), with specific reference to the treatment of atypical sources and related uncertainty. Starting from the existing records of historical tsunamis, we also offer a brief overview of the present state of knowledge about the most important sources for tsunamis that are not yet considered in present-day tsunami warning systems in Italy (Sect. 5). Finally, we outline the potential scientific paths toward the improvement of the efficiency and the source coverage of the tsunami warning service (Sect. 6).

Probabilistic tsunami hazard analysis and tsunami forecast for warning
TWSs and PTHAs are both based on forecasting techniques, where PTHA is focused on quantifying the long-term probability of exceedance of different tsunami intensities in the next years (usually tens of years) due to whatever tsunami that may occur, while tsunami forecasting for TWSs is instead oriented to forecasting the impact of an ongoing event.
Historically, tsunami forecasting procedures were based on numerical modelling of the tsunami generation and propagation process of one or a few scenarios. In the 1980s, some pioneering studies started better considering the source variability (e.g. [112][113][114]), as well as their probability of occurrence in the future. The occurrence of the 2004 Indian ocean tsunami, as well as the 2011 Tohoku tsunami, gave momentum to the development of more robust tsunami forecasting procedures, both for the management of the tsunami warnings and the development of appropriate planning strategies based on tsunami hazard, leading to the full establishment of probabilistic forecast approaches, progressively better and better incorporating uncertainty (e.g., [31,52,53,93,94,[115][116][117][118][119][120][121][122][123][124][125][126][127][128][129]).
Taking as a reference the equivalent development in the seismological field (e.g., from [130]; a recent review can be found in [131]), the attention was initially focused on the long-term forecasting procedure, namely on the Probabilistic Tsunami Hazard Analysis (PTHA; a recent review can be found in [11]). These studies mainly concentrated on seismically generated tsunamis (often named SPTHA) since they represent the most common source for tsunamis (> 75%, [36,132]) and also because SPTHA could be rooted in techniques already developed for seismic hazard (e.g., [115,129]). Progressively, SPTHA specialized in treating the peculiarity of tsunami sources and managing its inherent uncertainty, like techniques to account for the sensitivity of tsunamis to source magnitude, depth, and geometry or procedures to reduce the computational effort of inundation modelling (e.g., [31,93,94,133,134]).
The long-term forecast of PTHA consists of quantifying at each point x of the target domain the hazard curves, reporting the probability of exceeding different levels of tsunami intensity measure (using specific intensity measures like, for instance, wave height or momentum flux) at least once in a time window T , called exposure time. The hazard curves are usually evaluated in two steps [11]. First, the mean annual rate of exceedance is evaluated by combining a "source factor" with a "propagation factor", that is: where the source factor λ(σ ) represents the mean annual rate of the source scenario σ , the propagation factor P(z > Z ; x|σ i ) represents the exceedance probability at the target position x of the intensity Z, conditional upon the scenario σ , and the integration is made over the space containing all possible scenarios. The computation is usually performed over a discrete set of scenarios, as illustrated with the approximation by a discrete series in the right-hand side of Eq. (1). This type of approach, which combines source observations with numerical tsunami simulations, is most often preferred for direct quantification of tsunami hazard based tsunami observations, as they are typically few and sparse and hardly can represent the full variability of potential sources [11,54,65,115,135]. The evaluation of the source factor is based on seismic historical catalogues, geological constraints, and existing statistical laws (e.g., the Gutenberg-Richter distribution controlling the relative abundance of small vs. large earthquakes; [136,137]). This term is in many ways similar to the equivalent term in seismic hazard, with the difference that more parameters should be considered to describe tsunami generation (e.g., fault geometry and mechanism). The propagation factor is instead mainly based on the numerical modelling of tsunami generation, propagation, and inundation (see next section). This is very different from the seismic hazard, where this task is usually performed adopting empirical laws, the Ground Motion Prediction Equations (e.g., [138]). This is required to account for the tsunami impact patterns that are largely path-dependent. Yet, for the inundation phase, a stochastic approach to coastal evolution is sometimes introduced, based on tsunami observations and on numerical simulations (e.g. [93,139]), as discussed in Sect. 3.1.
The probability of exceedance in T is then quantified assuming that exceedance events are independent and occur randomly in time with a constant annual rate (stationary Poisson process) so that the probability of having one or more exceedances in T can be evaluated as: that is, it is evaluated as one minus the probability of no exceedances in T . As we will discuss in the next sections, intensity values to design long-term risk reduction measures (e.g., building codes, evacuation zones, etc.; [58,[60][61][62]) are defined either in terms of probability of exceedance in T or equivalently in terms of the mean return periods (often called average return periods-ARPs, reciprocal of λ(z > Z ; x) of Eq. 1). For example, seismic hazard design intensity usually corresponds to an exceedance probability of 2% or 5% in 50 years, corresponding to average return periods of 2475 and 475 years, respectively [131]. As discussed in Sect. 3.3, a common choice for the definition of evacuation areas for TWS is an average return period of 2500 years [60][61][62].
Given that hazard quantifications may be related to significant regulatory concerns, usually also the uncertainty on the hazard model is explicitly quantified [140]. This uncertainty is evaluated producing a number of alternative but "scientifically acceptable" models, that is a number of alternative estimations of λ(z > Z ; x). Such alternatives are usually not considered completely equivalent; rather, their different "credibility" is expressed by assigning weights to them in some specific framework. Different approaches have been developed to deal with this problem, ranging from logic trees to Bayesian and ensemble approaches [140][141][142][143][144]. This uncertainty is usually called epistemic uncertainty, as it arises from incomplete knowledge, giving rise to alternative models of nature, in contraposition to the so-called aleatory uncertainty arising from the natural variability of the phenomena [145]. Typically, epistemic uncertainty emerges from alternative models for the tsunami generation and propagation, alternative input and boundary conditions for source and propagation models, as well as alternative statistical distributions to describe the natural variability, for example, of the source (e.g. [11,53,54,115]). In principle, the epistemic uncertainty is reducible as new data and more accurate models become available, while aleatory uncertainty is not reducible, being related to the natural unpredictability. For this reason, epistemic uncertainty is usually kept separated, and represented as families of hazard curves and related statistics (mean hazard, percentiles, etc.; [11,131,142,145]).
Tsunami warning systems base their operational procedures on quasi-real-time tsunami forecasting of a possibly already ongoing tsunami. Nowadays, tsunami warning systems targeted to trans-oceanic seismic tsunamis (with travel times on the order of hours) are mostly based on the evaluation of the potential impact of one or few scenarios, selected inverting seismic and/or off-shore tsunami measures (e.g., [56]). When treated, near-source tsunami warnings are usually issued either based solely on source parameters or on near-source tsunami detection, possibly connecting to pre-computed tsunami scenario databases [59,[146][147][148][149][150]. Source parameters are often used as input to decision matrices that link magnitude, epicentre, and depth to tsunami intensity and alert levels. Such alert levels are usually also linked to specific actions in response to the tsunami (see, for example, Sects. 3.3 and 4.3). Quantitative and explicit management of the uncertainty affecting the short-term tsunami forecasting is still lacking, even if its importance has been already clearly acknowledged for a long time [95,96]. To this end, two main directions are being explored. On the one side, technologically-driven improvements of the monitoring network may and likely will substantially improve the time and accuracy of source parameter determination and tsunami detection [151,152]. On the other hand, formal procedures to quantitatively manage forecast uncertainty have been developed (e.g., [67,68,[153][154][155]) to be appropriately shaped for real-time use and accurately and systematically tested against real-data [68,156]. These probabilistic methods, similarly to long-term hazard, may be transparently connected to decision making for tsunami early warning, allowing for the separation between the scientific forecast regarding the hazard intensity and the political choices regarding safety levels and management of existing uncertainty [68,157,158].

Tsunami generation and propagation models
Numerical hydrodynamic modelling of tsunami evolution is an essential step of tsunami forecast procedures. Tsunami modelling consists of three fundamental steps: (i) the tsunami generation process, (ii) the propagation of the gravity waves generated by the source in the open sea/ocean, and (iii) the land inundation. A general review of modelling approaches and present-day challenges in the context of hazard and risk analysis can be found in [11,54] and in the references reported therein.

Tsunami generation models
The first step in tsunami modelling is the definition of the source and generation mechanism. As discussed in the Introduction, the most common sources for tsunamis are earthquakes, landslides, volcanic eruptions, and meteorological perturbations of the sea surface [11].
Earthquake sources generate tsunamis mainly by sudden deformations of the seafloor, inducing a disturbance of the overlying water body, which can be approximated as a sudden change of its potential energy. Tsunamigenic strength is primarily controlled by the fault dislocation (slip) and rupture areal extent (proportional to the seismic moment and hence linked to the magnitude), the geometry and the dynamics of the fault displacement, as well as the depth of the fault and the distance from the coast. Also, more complex features of the seismic source may have primary effects, like, for example, the spatial variability of dislocations (causing local slip concentrations) as well as spatiotemporal evolution of the seismic rupture, especially for large earthquakes featuring slow ruptures [159]. Seafloor deformation is modelled based on principles of continuum mechanics. Analytical expressions exist to calculate crustal deformation due to buried planar faults with instantaneous uniform slip in an elastic half-space exist [160,161], still representing a standard in many applications as they can also be extended to incorporate heterogeneous time-dependent slip in a simplified way. More advanced modelling approaches account for the details of Earth rheology, for example, considering layered elastic properties or even 3D variations through finite element models or models defining heterogeneous slip ruptures over non planar faults [92,[162][163][164][165][166][167][168]. While until recently tsunami forecasts were mainly based on the most simplistic approaches (uniform slip over planar faults), the recent scientific and technological advances allowed for the gradual inclusion of more complex modelling strategies [66,67,163,167,169], which made it imperative to cover complex slip distribution that can have first-order effects on the hazard. Even with these recent developments, many issues about seismic source modelling remain open [54], mainly related to the scarce knowledge about faulting, especially offshore, and to the variety and complexity of fault mechanics, with difficult-to-constrain features like slip distributions, potential surficial amplifications, slow events generating unusually large tsunamis (so-called tsunami-earthquakes, [15]), characterized by relatively slow dislocations.
Landslide sources vary from large submarine slides to more localized submarine landslides and subaerial landslides that enter the sea. As for the seismic case, the tsunami is primarily generated by the volumetric water displacement. Contrary to seismic sources, the time-dependent motion of the landslide plays a fundamental role in tsunami generation. Consequently, advanced models of tsunami generation from landslides should consider the landslide dynamics [16,18]. The main governing factors of landslide sources are acceleration, speed, volume, water depth, and geometrical configuration [16,170,171]. However, the generation mechanisms for submarine and subaerial landslides are distinctly different (e.g., the discussion in [17]). Submarine landslides are particularly sensitive to water depth and acceleration, hence not only the landslide dynamics is important, but also the slope failure rate and the mass mobilization rate, including remolding (e.g. [81,[170][171][172][173]). Subaerial landslides may impact the water at high speed, causing a different tsunami generation mechanism compared to submarine slides. The generation mechanisms scale non-linearly to the frontal area, Froude number, slope angle, and material density (e.g. [174][175][176][177]). Advanced tsunami generation models include the modelling of the landslide internal dynamics, with terms describing landslide rheology and density of different kinds, from viscoplastic models to frictional collisional models of different complexity and spatial averaging (e.g. [76-79, 81, 178]).
Volcanic sources are diverse and include source mechanisms potentially occurring on all volcanic phases, that is, in quiet periods, unrest episodes, as well as during eruptions. Tsunamigenic phenomena include slope failures and landslides, volcanic earthquakes, lahars-mudflow or debris flow composed of pyroclastic material-entering the water, as well as underwater explosions, pyroclastic flows, shock waves, or caldera subsidence [11,20,[179][180][181][182]. Notably, different potentially tsunamigenic processes can even concur in the tsunami generating phase, especially during large caldera-forming eruptions [20]. Similar to landslides, most of these source types require the joint modelling of source events and tsunami generation. However, some empirical models for some of these sources (e.g., explosions for which the water cavity scales with the energy forming the volcanic craters; e.g. [32,[183][184][185][186]), provide reasonable approximations also in the near field. Notably, some of these source processes overlap with seismic and landslide sources, like volcanic earthquakes, or flank collapses, slope failures and even pyroclastic flows [16,20], even though the peculiar material properties of volcanic systems and their potential transformations during the flow may force toward the development of ad hoc modelling strategies.
Strong, high-energy atmospheric phenomena may produce significant sea-level responses over a broad frequency band, which may generate significant oscillation due to resonance in the range between few minutes to hours [21,86,187,188] with dynamics significantly similar to tsunamis and different from wind-related long periodic storm surges. A very useful review of the processes underlying the generation of meteotsunamis can be found, for instance, in [86,87]. A moving atmospheric pressure disturbance and the related inverse barometer effect are the main driving source mechanisms. For a meteotsunami to be generated, the moving atmospheric pressure disturbance must attain suitable values in terms of magnitude, velocity (and hence of Froude number), duration and wavelength of the disturbance itself to be able trigger resonance phenomena, the most important of which is the Proudman resonance. Other, additional amplification mechanisms can occur when the pressure disturbance acts in the alongshore direction amplifying specific edge wave modes (Greenspan resonance) or due to topographic effects. In turn, the parent moving atmospheric pressure disturbance can be associated with very different atmospheric processes, such as squalls, thunderstorms, tropical and extratropical storms, atmospheric gravity waves. As each of these phenomena are often more frequent in specific geographical areas than in others, the processes that lead to the generation of a meteotsunami tend to be unique to a particular region, making them more frequent in favourable areas like Adriatic and the Baltic Seas, or the East Coast of the United States [22,189,190]. The tsunami generation may be modelled with coupled atmosphere-ocean models or coupling numerical wave propagation models with some source and parameter combination (e.g., [191,192]).
Other less frequent tsunami generating events exist. Among them, it is worth noting tsunamis generated by asteroid impacts, for which the source process scale from large impacts causing dry seabed to smaller events causing a surface cavity due to the direct impact of the asteroid or to its explosion near to the sea surface [11,[193][194][195][196][197][198]. Notably, especially for large-scale events, specific modelling approaches, including multi-phase (air, water, and solid) dynamics, are required to handle the thermodynamics and turbulence of the tsunami generation.
An important aspect of the deep-water tsunami generation process is that the water column here effectively works as a low-pass filter, suppressing short wavelengths (shorter than several times the water depth; [199]). For large earthquakes, this is only a relative problem, being dominated by long wavelengths, except perhaps for slip reaching at the trench in subduction zones (see anyway [200] and references therein for analytical treatment of the generation of tsunamis by underwater earthquakes).
On the contrary, smaller earthquakes, but particularly in deep-water landslide sources, need the introduction of low pass filters, with an efficiency that increases with depth and wave-number [41,75,199,201], to be applied at each step of the source evolution to avoid the introduction of short wavelengths determined by rapid accelerations in source dynamics. Recently, some authors performed 2D fully coupled earthquake and tsunami simulations of the subduction zones, documenting the smoothing effects induced by the nonhydrostatic response of the ocean at short wavelengths [202]. For landslides, unless a fully coupled multiphase Navier-Stokes model is used, filtering sources needs to be applied sequentially, which adds an additional computational burden (e.g. [16,170]).

Tsunami propagation and inundation models
The first numerical tsunami modelling procedures were developed in the early 80s [203,204], with the development of the first models based on a time-stepping scheme to solve the shallow water equations. To date, a large number of models have been developed, along with benchmarking and validation procedures against exact solutions for simplified problems, laboratory experiments, and/or past tsunami events (e.g., [205,206]).
The full three-dimensional Navier-Stokes equations are only rarely used for tsunami simulations: typically, they are needed only when non-hydrostatic and nonlinear effects may be relevant [207]. The effect of the curvature of the earth and of the Coriolis acceleration have only a minor impact on tsunami propagation [208] and they are usually neglected, but in very large domains to model trans-oceanic propagation where they play a more significant role. More commonly, different approximations are adopted based on the physics of the problem, with the purpose of making it computationally more accessible to enable a faster exploration of source variability. For most tsunamis, the vertical structure can be neglected, and two-dimensional depth-averaged systems of equations can generally be used, such as dispersive Boussinesq type or non-dispersive shallow-water approximation. In the deep sea, such equations may be linearized, being wave amplitude much smaller than sea depth, while the nonlinear shallow water equations should be solved approaching the shore (at depth < smaller than few tens of meters). This approximation is appropriate for sources dominated by long wavelength (e.g., large earthquakes), while strategies to include the vertical structure and to keep additional non-linear terms in asymptotic expansions are required for smaller and non-instantaneous sources like landslides and volcanic sources, as well as when nearshore coastal interaction and wave-induced currents in restricted environments (such as harbours or fjords) are significant. Many different sets of depthaveraged dispersive equations and layered non-hydrostatic models have been proposed through time to address these sources, often adopting dispersive Boussinesq solvers with more sophisticated numerical solution techniques (e.g., [209][210][211][212][213][214][215][216][217][218]).
Shoaling and inundation models are based on high-resolution bathymetric and topographic data, either including features such as buildings or vegetation or using Manning's roughness coefficients for accounting for their presence and, more in general, for the bottom friction. At the same time, inundation details may be sensitive to characteristics of the source, especially but not only in the near-field (e.g. [167,219]), thus, an appropriate level of details in the source process should be adopted depending on the distance of the source and characteristics of the target (see also [65]). Given the impossibility of extending high-resolution modelling techniques to entire basins, telescopic grids and different levels of approximations are usually adopted in the different phases of the propagation process. Non-linear shallow water models are most frequently used for this purpose (e.g., the review [64]); Boussinesq model applications are used as well, despite that they are more prone to instabilities (see discussion in [220]). Despite the relatively high accuracy of these modelling techniques, large uncertainty still exists (e.g., [54,63,69,92,127,221]), and the potential impact of uncertainty in initial/boundary conditions or of different modelling techniques still have to be fully understood.
Source characteristics are the primary driver for the selection of the appropriate modelling strategy. Tsunamis generated by large earthquakes are dominated by long wavelengths. These waves propagate more efficiently, reaching transoceanic distances, and are, in first approximation, non dispersive. On the contrary, landslide and volcanic sources, even if they encompass a broad range of source size, mostly generate localized tsunamis dominated by short-period waves with greater dispersion and limited far-field effects compared to earthquake-generated tsunamis (e.g., [19,91,184,[222][223][224][225][226][227][228][229][230][231][232]). Large volcanic explosions (e.g., caldera forming eruptions) and oceanic impacts of large asteroids may produce huge waves for which non-linearity may play a significant role for hundreds or even thousand of kilometres [197], with significant dispersive behaviours [196].
Boundary conditions also play a role. A typical computational domain used in tsunami simulations involves both offshore and coastal boundaries. In principle, tsunami waves should be fully transmitted through an offshore (open sea) boundary. In practice, the numerical implementation of these boundary conditions can lead to undesired effects, such as spurious reflections. To overcome this problem, several techniques have been proposed, including radiation conditions and absorbing layers (such as sponge layers and the "perfectly matched layer" method; see for instance [233,234] for a detailed discussion). The applicability and effectiveness of either methods is often dictated by the equations to be solved (shallow water, Boussinesq), by the specific numerical method adopted to solve the equations and by the geographical extension of the computational domain. The coastal boundary is treated in different ways depending on whether shoaling and inundation must be accounted for or not. In the first case, the complexity of the problem can vary significantly depending on whether the inundation process is simulated over a simplified topography or more realistic effects such as turbulence and the interaction with built environment must be taken into account: for the scope of this paper, it will be sufficient to address the interested reader to the review by Qin et al. [89]. To avoid unnecessary modelling of shoaling and inundation, reflecting boundaries are often used at shore-lines, with nearshore bathymetry corrections to manage, for example, minimum water depths. Since these conditions may introduce spurious reflections, a number of different approaches have been developed to manage such boundary conditions (e.g. [235]).
In regional studies, where inundation modelling is generally avoided, propagation modelling is stopped at the limit of linear shallow water approximations, that is several tens of meters depth (usually between 10 and 100 m, [11]). Then, amplifications models are applied to estimate the potential tsunami intensity (e.g., run-up) inland, ranging from a Green's law accounting for shoaling up to a bathymetric reference level, models amplifying solitary waves (e.g. [236,237]) or N-waves (e.g. [238][239][240]), or more advanced amplification models that account for the local bathymetric profiles and main characteristics (e.g. dominant period and polarity) of the incoming tsunami waves [139,241].
The availability of topographic and bathymetric data with suitable resolution and the capability of implementing them on properly built computational grids is a key yet often neglected aspect of any propagation and inundation modelling strategy. While this appears to be obvious when dealing with inland inundation, it is important to stress the importance of this aspect also for the propagation phase, where typical wave-related phenomena, like refraction and diffraction induced by the sea-bottom morphology, can be properly captured only by proper-resolution bathymetric data and computational grids.

Tectonic setting of the Mediterranean region
The Mediterranean region is characterized by basins and orogens resulting from a complex tectonic history started with the opening of oceans during the breakup of Pangea in the Triassic (∼ 250 million years ago), then followed by subduction and collision accommodating the convergence between the African (or Nubian) and Eurasian plates since the Jurassic (∼ 200 million years ago) [242]. Subduction also generated back-arc basins in the upper plate, leading to the formation of extended continental crust or even new oceans, such as the Gulf of Lion Basin, the Algerian Basin, and the Tyrrhenian Basin in the western Mediterranean [242]. As of today, the tectonic activity within the Mediterranean region is mainly driven by the still active convergence between Africa and Europe, in which the deformation pattern is further complicated by the presence of smaller plates (like the Aegean and Anatolian plates) or promontories (like the Adriatic block). A sketch map of the Mediterranean seismotectonic settings is reported in Fig. 1.
The boundary between the African and Eurasian plates takes a different character from place to place. From west to east, in southern Spain and northern Morocco across the Gibraltar Straits, one finds the Betic-Rif orogen surrounding the Alboran Basin. Then the Maghrebides orogen characterizes the western part of northern Africa from Morocco to Sicily and the Messina Straits. Next to the Messina Straits, southern Italy, the Calabrian Arc subduction accommodates the NW-directed subduction of the Ionian Basin lithosphere. Then the plate boundary bends sharply by wrapping around the Adria microplate, a promontory of the African plate and common foreland of the collisional belts of the East-directed Apennines, characterized by an extensional and contractional pair, the Alps, and the West-directed Dinarides, Albanides, and Hellenides. Here the plate boundary changes again its character into subduction with the Hellenic and the Cyprus Arcs, where the African plate is overridden by the Aegean and Anatolian plates (comprising Greece, Cyprus, and Turkey). The rates of convergence across the three subduction zones of the Mediterranean (Calabrian, Hellenic, and Cyprus Arcs) are in the order of few millimetres per year in the Calabrian Arc and Fig. 1 Seismotectonic map of the Mediterranean region. Earthquake locations were taken from the regional centroid moment tensors catalogue (RCMT; [389]) with mechanisms classified based on the Zoback [390] rules (TF thrust faulting, TS thrust strike, SS strike-slip, NF normal faulting, NS normal strike, UN undetermined). Tectonic structures were modified and simplified from [391] and various other regional geologic maps. Volcanoes from [392]. Topo-bathymetry from [393] in the eastern part of the Cyprus Arc [243][244][245][246], and of few tens of millimetres per year in the Hellenic Arc and in the western part of the Cyprus Arc [245,[247][248][249]. All other orogenic belts, in general, accommodate only a few millimetres per year of contraction. Extension in the order of a few millimetres per year characterizes the axis of peninsular Italy (inner Apennines chain) and the Aegean region, with faster extension rates in the Gulf of Corinth.
This complex geological setting is associated with frequent and intense seismicity throughout the Mediterranean basin. Offshore earthquakes in the region include different types of events: from subduction-related earthquakes, in the Calabrian, Hellenic, and Cyprus arcs, to crustal events, on thrust faults belonging to the Maghrebides, Apennines, Dinarides, Albanides, and Hellenides fronts, on normal faults in the Messina Straits, the Aegean Sea, and Corinth Gulf, and on strike-slip faults such as the Kefalonia-Lefkada (Ionian Sea), Paphos (Cyprus), and North Anatolia fault (Aegean and Marmara Seas). In association with this complex tectonic setting of oceanic and continental subduction/collision, several volcanoes have developed throughout the Mediterranean region, including syn-and post-orogenic volcanoes, many of which lie near the coasts or directly on the seafloor. Most of the active volcanoes are located in the back-arc of the Calabrian and Hellenic Arcs, in the Tyrrhenian and in the Aegenan seas, along with few other important volcanoes such as Mt. Etna, which is one of the most active volcanoes in the World, Mt. Vesuvius, Campi Flegrei, and Ischia in the Tyrrhenian coasts of southern Italy, as well as several volcanic edifices in the Sicily channel.
For the subduction-related events, either on slab interface or splay faults, magnitudes as high as 8 + have been hypothesized based on inferences from historical and geological data [250,251]. For crustal events occurring outside subduction zones, the maximum recorded magnitudes are slightly below magnitude 8, but with high tsunamigenic potential due to shallow faulting and dip-slip motion (when faulting mainly induces vertical displacement), as for the 1908 Messina-Reggio Calabria M7.1 or the 1956 Amorgos M7.8 earthquakes. Indeed, for both events, a tsunami run-up of up to 13 m and 25 m was reported, respectively [252][253][254][255]. It is possible that some of these events have triggered landslides, as hypothesized for both the 1908 case [50,253,256] and for the 1956 case [255,257], that could have increased the tsunami impact. More information about the past seismicity and the tsunami records in the Mediterranean area have been systematically collected in several historical and instrumental catalogues [36,45,[258][259][260][261].

Tsunami warning for earthquake-generated tsunamis
In the following sections, we briefly discuss hazard quantification, coastal planning, and tsunami warning in Italy for what it concerns seismic sources. Particular attention is given to the specific management of crustal sources, which constitute the large majority of potential tsunami sources in the Mediterranean.

Hazard quantification
The NEAM region (Northeastern Atlantic, the Mediterranean and connected seas) has its own recent SPTHA model, called NEAMTHM18 (NEAM Tsunami Hazard Model 2018; [51][52][53]), produced by the TSUMAPS-NEAM project (http://www.tsumapsneam.eu/). NEAMTHM18 is a hazard model obtained through the collaboration of a large scientific community, with the goal of building the input for a common and homogeneous strategy for tsunami risk management in the region. The methodology employed and the results are documented by Basili et al. [52,53].
NEAMTHM18 considers subduction zones and both inter-plate crustal seismicity (including oceanic ridges) and intra-plate diffuse seismicity. Subduction zone seismicity in the accretionary wedge, the continental crust, including splay faults, and outer-rise earthquakes in the oceanic crust are also included. As a general rule, NEAMTHM18 assumes that crustal earthquakes can take place everywhere. Only in few areas, when the regions are very stable and hence the seismicity rates are very low, the intraplate seismicity is disregarded, as in the oceanic crust in the Atlantic Ocean far from the ridges and other principal faults. No such areas are present in the Mediterranean.
The knowledge of the potential earthquake sources is somehow limited in many areas, and certain faults are known better than others [14]. To optimally deal with this heterogeneous degree of knowledge while maximizing the use of all the available infor-mation, the seismicity is subdivided into two categories, background seismicity, used for treating crustal earthquakes, and predominant seismicity, used for dealing mainly with subductions. Either type of seismicity adopts a different modelling approach for one or more seismic source parameters. Some values of the earthquake parameters can even be set constant when their variability is considered negligible with respect to other parameters. One example is fault geometry on subduction interfaces, which can be considered relatively well constrained with respect to other source parameters. This approach to seismicity types in probabilistic calculations was introduced into SPTHA by [31], in line with other approaches introducing faults into Probabilistic Seismic Hazard Analyses [262][263][264].
The sources considered, when focussing on the coasts of Italy, may be restricted to those located within the Mediterranean basin. Three subduction zones are taken into account (the Hellenic Arc, the Calabrian Arc, and the Cyprus Arc, the latter probably representing a relatively minor threat for Italian coasts). The subduction zone characterization and modelling in NEAMTHM18 are detailed in several papers [52,53,163,265,266]. We only note that a magnitude range from 6.0 to 9.1 over the 3D subduction geometries is used, and different alternatives are used for building slip distributions, including depth-dependent rigidity hence potentially enhanced shallow slip. The background crustal seismicity is instead modelled everywhere on a regular grid covering the entire Mediterranean, considering a magnitude range from 6 to 8.1, a depth range from the surface to the Moho, and 144 potential focal mechanisms exploring the variability of the strike, dip, and rake, for a total of more than 10 3 scenarios in each grid point. At each node of the grid, the distribution of faulting mechanisms is constrained by moment tensor catalogues and mapped faults [31]. Collectively, the number of subduction and background scenarios explored is very large (> 10 7 ), the majority of which are background seismicity. Hazard disaggregation results show that for the largest part of the Italian coasts, the background seismicity dominates over subduction seismicity, especially at the largest tsunami intensities [31,52,53], demonstrating the importance of including atypical sources in areas like the Mediterranean.
Tsunami propagation from each individual source is modelled up to a set of Points Of Interest (POIs). Considering the challenging number of individual scenarios, an ad hoc computational strategy has been adopted [66], mainly based on the tsunami wave reconstruction in the deep sea with linear combinations of elementary Gaussian sources or in some cases rectangular subfaults, modelled with the Tsunami-HySEA code (a finite volume solver for the non-linear shallow water equations developed for multi-Graphics Processing Unit-architectures; [267]) into a High-Performance Computing (HPC) environment. The POIs have been selected offshore approximately on the 50 m bathymetric line at a distance of approximately 20 km from each other. POI locations can be considered representative of a roughly 20 km long stretch of the coast behind them. Tsunami inundation is evaluated by adopting an amplification model that considers the tsunami characteristics (polarity, dominant period, and wave height) at the POI and the local bathymetric profile. The significant uncertainty introduced by this simplified model for inundation is propagated into the PTHA results.
The tsunami intensity is measured in terms of the Maximum Inundation Height (MIH; [93,139]). MIH is the maximum height reached by the water with respect to the sea level at rest (water + topography) on profiles orthogonal to the coast, and thus it can be higher than the run-up (maximum topographic altitude reached by the flood). MIH then represents the maximum inundation height averaged along the~20 km long stretch of coast represented by each POI. It follows that the locally maximum MIH and run-up values can be (up to 3-4 times) larger than those indicated by the hazard model, for example, due to tsunami energy focusing by topographic features [53,62,63,139].
In agreement with the PTHA approach, NEAMTHM18 evaluates results in terms of a suite of hazard curves, reporting the probability of exceeding different MIH values in 50 years in each POI (Fig. 2). More specifically, NEAMTHM18 provides mean hazard curves, and the model uncertainty is reported through different curves at different percentiles (Fig. 2b).
Being the most recent available hazard model including all Italian coasts, NEAMTHM18 was adopted as a starting point for the Italian national coastal planning for seismically-generated tsunamis [62], while working on a national hazard model (see Sects. 3.3 and 3.4). NEAMTHM18 is indeed the only one today that uses a homogeneous approach over the entire Italian coast, and which deals with all types of seismic sources. Cautious usage of the model is suggested when used for local planning due to the relatively coarse spatial scale of tsunami probability output and the relatively coarse sampling of the source parameter variability. Both scales of discretization are suitable for a large (NEAM region) scale model but certainly less well-calibrated for local (national) scale hazard analysis. For this reason, the practical use of NEAMTHM18 in coastal planning has been integrated with a specific strategy to account for its spatial resolution (Sect. 3.3). CAT-INGV operates as NTWC (National Tsunami Warning Centre) for the whole Mediterranean region (Fig. 3a), from Gibraltar to Eastern Mediterranean, and serves as NTWC and as TSP for the DPC (the Italian civil protection) and for several member States of UNESCO and European Institutions. CAT-INGV started monitoring the Mediterranean Sea in autumn 2014 as a candidate TSP. Monitoring at a global scale is also performed as a permanent training activity. Dedicated personnel are present on a 24/7 basis at the Seismic Surveillance Room of the INGV in Rome, assisted by a supervisor officer on call H24.

Italian Tsunami Warning System for seismically-induced tsunamis
After some years of pre-operational mode, CAT-INGV has been accredited as TSP by IOC/UNESCO in 2016, including the entire Mediterranean Sea in its monitoring area for tsunami sources. Since then, it has delivered seven tsunami alerts, including three Advisory messages and four local Watch messages, and several information messages for earthquakes that occurred in the Mediterranean region with magnitudes up to 7.0 [110,111]. As all the five accredited Tsunami Service Providers (TSPs) of the NEAMTWS, CAT-INGV is in charge of providing tsunami forecasting when earthquakes of magnitude greater or equal to 5.5 occur in their monitoring area. TSPs use combinations of global, regional, and national seismic and tide-gauge networks.
In Italy, the national seismic tsunami warning system for the Italian coasts (SiAM-Sistema di Allertamento da Maremoti) operates since January 2017 [109] (http://www.protezionecivile.gov.it/en/risk-activities/tsunami-risk/activities). Seismic monitoring, first alerting, and tide-gauge analysis are carried out by CAT-INGV, based on predefined and agreed rules. As in all NEAMTWS TSPs, CAT-INGV's Zemmouri-Boumerdes earthquake: red/orange/green triangles indicated watch/advisory/information levels in the warning message; the dotted circles indicated local/regional/basin spatial domain, as defined by the decision matrix (Table 1); the red star indicates the epicentral location. c As b, for the M6.8 2017 Kos-Bodrum earthquake tsunami forecasting relies on a decision matrix that converts the main earthquake parameters (hypocentral location and magnitude) into alert levels at predefined forecast points along the potentially threatened coastal locations. Forecast points generally correspond to the locations of sea-level gauges and/or to highly exposed/vulnerable locations. The national mareographic network (Rete Mareografica Nazionale, RMN, https://www.mareografico.it/) is managed by ISPRA. CAT-INGV communicates alerts to DPC, which is in charge of the dissemination to the operational structures and components of the National Civil Protection Service, aiming at reaching, in the shortest possible time, the population who would be potentially affected.
In NEAMTWS and in SiAM, two alert levels are foreseen: • Advisory/Orange: the detected seismic event could produce a tsunami wave height ≤ 0.5 m in front of the coast and run-up ≤ 1 m above the sea level; this may represent a marine and near-coast tsunami threat; • Watch/Red: the tsunami wave height is expected to be > 0.5 m in front of the coast and the run-up > 1 m above the sea level; hence, it may cause a significant inundation.
Moreover, information messages are sent when significant earthquakes (M ≥ 5.5) occur in the competence area, but it is unlikely that a tsunami will impact the coast because, according to the decision matrix, the seismic event is evaluated as not able to produce a tsunami wave that represents a relevant threat for the exposed coasts.
In NEAMTWS, there are two different decision matrices, for the Mediterranean and for the North-East Atlantic. Moreover, slightly different decision matrices, especially for the chosen earthquake magnitude thresholds, are used by the different TSPs. Each TSP has presented its own decision matrix during the accreditation procedure. The one currently used at CAT-INGV is presented in Table 1.
Since these decision matrices are based on the analysis of historical events, tsunami modelling and expert judgement, one can argue that they implicitly account for all types of seismic sources (including crustal non-subduction earthquakes); however, despite some degree of intended conservatism, as discussed below, decision matrices were certainly not conceived for including secondary events (e.g., seismically induced landslides), as present in the historical record (see Sect. 4). Those secondary events could have been initially accommodated in this scheme by operating on the thresholds for distance and magnitude.
More specifically, the decision matrices in the NEAMTWS are distance-based, so that the potential threat decreases with the distance from the earthquake source. Tsunami ranges (local, regional, basin-wide) are referred to as the distances between the earthquake epicentre and the forecast points. Each range is associated with a circle of fixed radius, as reported in Table 1. They do not take into account the source orientation (described by strike and dip angles) and mechanism (described by the rake angle) of the earthquake. Also, for this reason, while considering only tsunamis generated directly by seismic sources, decision matrices are often quite conservative on average, in the sense that they are "worst-case" oriented. Figure 3b, c reports the alert levels obtained by applying the decision matrix of Table 1 [110,111]. The colours red/orange/green correspond to watch/advisory/information levels in the warning message, respectively.
Based on the real-time analysis of seismic data, TSPs strive to issue initial alert messages within a few minutes (< 10 ) from any strong earthquake occurring at sea or near the coasts. At CAT-INGV, seismic parameters are obtained through the software Early-est [268][269][270][271][272], which elaborates in real-time data from hundreds of seismic stations worldwide and calculates a series of progressive localisations, with magnitude and hypocentre determination. While the first estimate is produced about 2 min after the earthquake, a good compromise between speed of calculation and accuracy was found with the Early-est 5th solution, generated about 7-8 min after the earthquake occurred. If no specific issues occur, this solution is used as input to the decision matrix and the consequent alert level is delivered within 10 min from origin time (in SiAM Informazione or Allerta Iniziale; in NEAMTWS Information or Advisory/Watch). In recent events, however, initial alert messages have been issued based on the Early-est 2nd solution, thus gaining 1-3 min; indeed, for both May 2020 (Crete) and October 2020 (Samos) events, initial messages were issued 8 min after the earthquakes' origin times. If seismic parameters significantly change afterward, an updating message (in SiAM Aggiornamento/Update; in NEAMTWS Ongoing) is delivered.
After the first alert is issued, eventual sea level anomalies around the earthquake epicentre are observed through the RMN (https://www.mareografico.it/) as well as the other national mareographic networks whose data are collected and distributed by IOC/UNESCO tide-gauges (http://www.ioc-sealevelmonitoring.org), by the JRC (Joint Research Centre of the European Union; [273]), or exchanged as a result of bilateral agreements.
The RMN and the other national mareographic networks operating in the Mediterranean consist of telemetered tide gauges only. The Italian network is composed of 36 tide gauges located in the harbours and/or close to the coast and equipped with a radar sensor coupled with a second floating back-up sensor based on a "shaft-encoder" technology. UMTS public channels are used for data transmission, furthermore signal redundancy is ensured through the IRIDIUM network. In the period 2014-2016, JRC started an installation campaign of a total of 20 IDSL (Inexpensive Device for Sea Level Measurements) devices in the NEAMS area integrating the national monitoring networks [273]. A floating GPS buoy suitable for the offshore installation has been tested in 2018 by JRC and ISPRA in the La Spezia Gulf. The testing returned good results but the device is not yet in production. In addition, recently a thematic board between ISPRA, INGV and MISE (Ministry of Economic Development) started analysing the possibility to put instruments in the dismissed submarine cables present in the Mediterranean Sea floor with bottom pressure sensors. Despite all these valid initiatives, unfortunately, to date, bottom pressure sensors like DART buoys or other sea-bottom real-time instruments (e.g. [151,152]) are not operational in the Mediterranean.
In case of significant sea-level anomalies in one (or more) of the tide gauges located near the epicentre, a confirmation message (in SiAM Conferma; in NEAMTWS Ongoing) is delivered, reporting also sea level observations, and an ending message (in Siam Fine Evento; in NEAMTWS Ending) is sent when all sea-level measures go back to background values. A cancellation message (in SiAM Revoca; in NEAMTWS Cancellation) is instead delivered when no sea level anomalies are detected.
All the messages produced by CAT-INGV are delivered in Italian to DPC and in English to all other countries/organizations which have subscribed to its service, based on predefined formats. DPC automatically disseminates these messages to all the interested emergency management system operators, including those at the local level. At present, the national tsunami warning system, SiAM, does not directly reach the population with its alerts.

Evacuation zones and civil protection coastal planning
As PTHA (and especially the Seismic-PTHA, SPTHA) is emerging as a standard method for tsunami hazard quantification [11,65], coastal planning is currently based on the most recent available SPTHA model that homogeneously covers all Italian coasts. As discussed in Sect. 3.1, this is NEAMTHM18 [51][52][53].
In SiAM, both Advisory/Orange and Watch/Red levels are connected to evacuation, corresponding to warning zones of type 1 and 2, respectively. Warning zones of type 1 are the inundation zones corresponding to a maximum run-up of 1 m or a wave height at the coast of 0.5 m above sea level. These thresholds directly derive from the definition of the alert level (Sect. 3.2). To define the tsunami inundation zones of type 2 to be used for evacuation in case of a Watch message, it is necessary to define an upper-bound since the Watch message corresponds to a possible run-up greater than 1 m. For this reason, the DPC chose as a reference the MIH value corresponding to the 2500 years average return period on the 84th percentile curve of the epistemic uncertainty in the reference hazard model NEAMTHM18 (Fig. 2b). This choice is consistent with what is suggested for evacuation mapping in New Zealand and for tsunami building codes in the US [58,60,61]. It is worth noting that the type 2 inundation zone is also used for coastal planning, which is discussed later in this Section.
The full procedure describing the conversion between MIH to maximum run-up and inundation areas is described in detail in [62,63]. Nevertheless, we provide here a short description of this procedure for the sake of completeness.
Through the inclusion of several safety factors, this procedure is meant to account for the limitations in the spatial resolution of the input hazard model (Sect. 3.1) and the limited resolution of digital elevation models available for all the Italian coasts. In short, to take into account the MIH fluctuations over the stretch of the coast behind the POI, a multiplicative safety factor of 3 is applied to translate MIH into maximum run-up along the coast. This factor corresponds to a very high percentile (> 95th) of the distribution of the ratio between the maximum run-up and the average MIH on stretches of coast, as modelled through nonlinear shallow water inundation for numerous tsunamis with different polarities and periods at different locations [52,53,63,139].
The limited spatial resolution and the relatively coarse sampling of the source parameter space (for example, the fault centres are discretized with a step of about 25 km, and their strike angles-indication fault orientation-with a 45°step) of NEAMTHM18 may hide un-modelled inter-POI finer scale run-up fluctuations. To account for this, the design run-up values were chosen by taking the maximum between several adjacent maximum run-up values within a search radius of 40 km. A specific procedure was reserved for some areas that had peculiar geographical and morphobathymetric features: the Northern Adriatic Sea, the southern Tyrrhenian Calabria, and several smaller islands.
To evaluate the inundation zones from maximum run-up values in a large coastal region, the so-called "bathtub" approach is the easiest and fast method. This hydrostatic approach, in which any coastal location with an altitude lower than the design run-up would be flooded, does not consider the dissipation of waves for the inland propagation and may provide overestimations, especially for the coastal areas with a flat morphology. To account for dissipation, a simple GIS-based approach was used, introducing an empirical dissipation factor in a rule linking the maximum run-up values with the maximum distance of inundation, and considering nominal Digital Elevation Model (DEM) uncertainty. The maximum inundation distance is calculated in 200 m per meter of the design run-up value. This rule is applied if the design runup is not already reached on the local topographic model (Fig. 4). The dissipation is halved (hence the maximum allowed inundation distance is doubled, 400 m per run-up meter) in the presence of rivers. These empirical relationships are very similar to those proposed for New Zealand and estimated on the basis of the numerous observations conducted following recent and historical tsunami events, which mainly occurred in the Pacific area [55,60,61,274]. This approach should, in principle, produce more realistic results than the simple hydrostatic model, even if it certainly does not take into account the actual physical behaviour in a complex coastline. More complex effects like interferences, interactions with territory and infrastructures, and between waves are to some extent incorporated when applying the safety factor and through the adoption of the maxima within search radii.
The procedure described above is applied to a predefined set of reference maximum run-up values (2, 5, 10, 15, and 20 m) for all the Italian coasts. In each area, to the design run-up value obtained from the hazard analysis, it is summed 1 m to account for the uncertainty in the Digital Terrain Model, and the reference maximum run-up value just larger than the obtained value is selected to define the evacuation areas. The obtained evacuation areas, available through the Tsunami Map Viewer (http:// sgi2.isprambiente.it/tsunamimap), finally were tested against a few data available in the Mediterranean, for example, by comparing the obtained inundation zones to the inundation values observed for the historical Messina 1908 tsunami. More specific testing is discussed in [63].
In 2018, Operational Guidelines were issued by the Head of DPC to support civil protection planning for tsunami risk [62]. These Guidelines, issued as an implementing act of the Directive that established the national warning system for seismically generated tsunamis (SiAM), refer to civil protection planning for the management of the risk arising from possible tsunami waves generated by earthquakes. The document provides indications to all levels of the public administrations and operational bodies of the National Service of Civil Protection for the protection of the coastal population from tsunami events.
It is useful to clarify the meaning of "civil protection planning" in the Italian system. Planning does not only mean organizing the actions of civil protection operators to intervene when an emergency occurs. A plan is also intended as a tool to increase risk awareness in "ordinary times" and encompasses pooling of resources, capacity-building activities among professionals, guaranteeing the link between different administrations and authorities. Civil protection planning is a consequent system-wide activity, to be carried out jointly by all local authorities and bodies involved in emergency preparedness and management, whether public or private. They provide a consistent and unique set of reference procedures and terminology. The plans must be constantly updated to accommodate changes in land use planning and flexible enough to be used in all emergencies, including unexpected ones.
In the Guidelines to support civil protection planning for tsunami risk in Italy, alert levels and corresponding expected tsunami inundation areas are defined for each part of the Italian coastline. In particular, national guidelines include: • a detailed explanation of the tsunami early warning system at the national level; • the tsunami inundation areas for each part of the Italian coastline for the two considered alert levels, "advisory" and "warning", corresponding to warning zones 1 and 2, respectively; • indications on the required contents of the civil protection plans, with specific sections dedicated to the municipal, provincial, regional level and one section dedicated to civil protection plans of operational bodies and lifeline providers; • reference material for the tsunami signs, which are recommended to be installed, and best practices on tsunami public alert communication measures.
The Guidelines foresee that plans at the territorial level must be defined according to the different SiAM alert messages produced by DPC on the basis of CAT/INGV messages: Informazione (information), Allerta (alert), Aggiornamento (update), Conferma (confirmation), Revoca (cancellation) and Fine Evento (end of the event). The various institutional levels have to translate the information received with the alert message in procedures to address the specific needs of the territory under their jurisdiction. According to the current legal framework, the public warning is the responsibility of Mayors; therefore, the population must be alerted through the procedures defined in the municipal management plans, with the support of public administrations at territorial levels (regions and prefectures).
The general strategy adopted is the preventive evacuation of the population present in pre-defined coastal zones at risk (zones 1 and 2). This evacuation can be either vertical, moving upward to the highest floors of buildings suitable to withstand the tsunami impact or to higher topographical heights, or horizontal, moving inland. The operational procedures should be inspired by those foreseen for the seismic risk, with the necessary adaptations, starting-for example-with verifying the location of the headquarters of the operational coordination centres with respect to the expected inundation areas and corresponding warning zones.
The plans at regional and provincial levels (Regions, Prefectures, and Maritime Authority) must both ensure the necessary support to the activities of coastal munic-ipalities during an emergency, and promote the organization of exercises and local hazard and vulnerability studies, to provide a higher level of detail for the tsunami risk management initiatives in their territory. The largest part of the Guidelines is the one addressed to Municipalities, which should prepare the municipal plans in four main sections: local hazard and corresponding warning zones; public warning procedures; operational response model and connected activities (e.g., installation of emergency signs for tsunami risk); communication plan and related activities (e.g., to increase community awareness and preparedness).
In the case of tsunamigenic earthquakes very close to the Italian coast, messages will most probably not reach near-field areas in time to activate preventive measures. In those cases, it is up to the citizens' ability to autonomously recognize potential precursor phenomena and immediately implement self-protection behaviours.
In this respect, periodic risk information campaigns and exercises become a crucial activity to engage, inform and develop a self-protection culture among the coastal population, such as the yearly nation-wide public information campaign "I don't take risks" (http://iononrischio.protezionecivile.it/en/homepage/), "I don't take risks-Tsunami", in particular, started in 2013 within the international exercise Twist-Tidal Wave In Southern Tyrrhenian Sea, funded by the European Commission.
In terms of exercises, SiAM (or its components) has participated in all the regional NEAMWave exercises organized by the Intergovernmental Coordination Group (ICG) of the NEAMTWS (NEAMWave12, NEAMWave14, NEAMWave17) and is going to participate in NEAMWave2021 with an Ionian Sea earthquake-tsunami scenario.
Starting from 2020, Italy has started to pursue the Tsunami Ready initiative, following the examples of the US National Weather Service (NWS) of the National Ocean and Atmospheric Administration (NOAA) that first launched this program in 2001, then followed by all TWS worldwide, including NEAMTWS in the last year. The goal of Tsunami Ready is to assist coastal communities to be better prepared to save lives through better planning, education, and awareness, involving all community levels.

Ongoing developments
Two main developments are being carried out to improve the forecasting procedures: the revision of the tsunami hazard model at the national scale and the development of a new real-time probabilistic tsunami forecasting. In addition, several activities are ongoing to improve coastal planning.

Toward a national S-PTHA
NEAMTHM18 is designed to produce a homogeneous assessment over the entire target area, extended to all the NEAM region (Fig. 2a). Limiting the target area to the Italian coasts would allow refinement of the hazard model; in fact, a specific SPTHA model for the Italian coasts is under development. This model is largely based on NEAMTHM18, but the following updates are being evaluated: • The use of updated local databases, like for example, the ones in development within the new European PSHA model [275] and the consequent re-evaluation of catalogue completeness (defining the minimum magnitude above which all the events are present in the catalogue in a given temporal window); • A revision of the regionalization (the process of defining homogeneous areas within the source domain) and all statistical models; • The improvement of the sampling of aleatory uncertainty, updating the sampling procedure for the sources most relevant for the Italian coasts; • The use of nonlinear shallow water simulations to update the tsunami intensity estimations of the scenarios relevant for Italian coasts. • The use of a set of target POIs more densely distributed along Italian coasts (from 20 to~5 km of average distance among POIs), also considering propagation in water shallower than 50 m (e.g., up to 10 m); • The update and refinement of the amplification model for the Italian coasts (similarly to [139]).
Once the national S-PTHA model is finalized, it may also be required to review the inundation zones, using a more accurate DEM, if available. Moreover, specific approaches could be adopted in the presence of, for example, high population density or critical infrastructures, and the GIS-based approach might be integrated with highresolution numerical modelling [63].

Toward a tsunami warning based on real-time probabilistic tsunami forecast
Considering that the alert levels are connected to the expected local tsunami intensity at the forecast points, decision matrices provide a basic tsunami forecasting procedure. As already discussed earlier, the main advantage of using decision matrices is the fact that forecasts and alert levels may be timely produced as soon as basic information about the source earthquake is available. Location, depth, and magnitude are typically available and sufficiently stable just a few minutes after the origin time [272]. More sophisticated source parameters (such as source geometry), which usually are produced later in time with respect to localization and magnitude, are not required for decision matrices.
There are several known drawbacks when decision matrices are adopted as the basis for tsunami forecasting. Some examples are shown in Fig. 5, where we report the results of the application of the decision matrix to the Norcia 2016 (M6.5) event, occurred at the core of the Apennines chain, and to the Peloponnese scenario (M8.5), adopted in the NEAMWave17 exercise [276].
In the case of the Norcia earthquake, both the Adriatic and Tyrrhenian Seas were within the areas with Advisory and Watch alert levels, even if the likelihood of tsunamis from an earthquake located in that area (usually faulting oriented NW-SE along the mountain chain) is clearly quite low. In the NEAMWave17 scenario, instead, the whole basin would be in a "Watch" status, while simulations show that many areas of the Mediterranean would be only slightly affected by a similar scenario. These issues are due to the fact that the decision matrix, as discussed in the previous section, does not consider the following: (1) any preferential earthquake fault orientation and mechanism, which in turn determines the extent of induced seafloor displacement and the directivity of the tsunami energy; (2) the tsunami propagation, not even at the first order. For example, in the specific case of the NEAMWave17 Peloponnese scenario,  Fig. 3, red/orange/green triangles indicate watch/advisory/information levels in the warning message; the dotted circles indicate local/regional/basin spatial domain, as defined by the decision matrix (Table 1); the red star indicates the epicentral location tsunami propagation would be limited through the Messina Straits and, to some extent, the Sicily channel, as well as within the Aegean Sea because of the chain of islands, and surely through the Dardanelles towards the Marmara Sea.
However, this cautionary approach not only leads to a large number of false alarms (when the tsunami is smaller than expected), but it may also be misleading, since it does not necessarily completely avoid the occurrence of missed alarms (when a tsunami larger than expected occurs) at specific locations, where energy focus may occur. Complex and asymmetrical propagation patterns of tsunamis, like the directivity of tsunami energy flux and the azimuthal anisotropy of the tsunami propagation on a complex bathymetry, may lead to important tsunami amplifications also far from the seismic source. Hence, depending on the relative source-forecast point positions, despite decision matrices are in average conservative, under-and over-estimation of the alert level may be expected (e.g. [68,277]). For example, for the 2003 Zemmouri-Boumerdes M6.8 tsunami (Fig. 4), according to the decision matrix, the Balearic islands are inside the advisory alert level (forecasting a tsunami with a run-up < 1 m or wave height < 0.5 m above sea level), but local amplifications lead to a tsunami intensity larger than expected [37,278].
To account for propagation is, however, a complex scientific challenge, considering the very vast variability of the potential sources and the need to provide robust estimations within the short times required in relatively small basins like the Mediterranean. To this end, a new probabilistic method coined "Probabilistic Tsunami Forecasting (PTF)" is being developed [68]. The PTF is based on the propagation of the uncertainty from the source parameters, as estimated in near real-time from the monitoring room, to the potential impact zone, through pre-calculated tsunami simulations and relative uncertainty. The method allows updating the input data through time, as new information about the seismic source and/or the propagating tsunami is available. The PTF provides as output the probability distribution of the tsunami inundation height at predefined target points.
To overcome computational issues and assure timely estimations, the PTF makes use of the pre-computed database of tsunami simulations derived from the regional hazard model NEAMTHM18, which includes, as said, both subduction and crustal seismicity [51][52][53].
In relatively small area like the Mediterranean, a very rapid assessment of the PTF is required, to deal with near-field tsunamis (with warnings to be delivered within 10 /15 from the earthquake occurrence; see Sect. 4.2 and references therein). To this end, the prototypal implementation of the method is based on the evaluation of the uncertainty on magnitude and epicentre from the seismic location algorithm adopted in the monitoring room (in CAT-INGV, Early-Est [272]). On the contrary, the estimation of the focal mechanism (defining geometry and relative movements on the fault through three angles: strike, dip, rake) is typically produced too late to be included. To compensate this lack, estimations derived from long-term hazard quantifications can be used. For example, Selva et al. [68] adopt the results of NEAMTHM18, which extend to all the source area the method developed in [31], in which the probability distribution describing the uncertainty on the source mechanisms in each source area is set based on local geological (e.g., mapped faults), historical and instrumental (e.g., earthquake focal mechanisms) information. New rapid source inversion techniques are rapidly progressing [268,[279][280][281][282][283][284][285] and will probably reduce this uncertainty in the future. In any case, in the Mediterranean region, the possibility to quantify this uncertainty is particularly important since it allows dealing with crustal seismicity, which is generally less constrained than subduction seismicity, for which it is customary in several TWS worldwide to assume a pure thrust faulting occurring at a depth of the subduction interface corresponding to the epicentral location.
In Fig. 6a, we report an example of the results for the M6.8 2003 Zemmouri-Boumerdes earthquake and tsunami. The figure shows a map of the probability of exceeding 0.1 m (map in the background). It can be noted that the results follow the expected tsunami focusing patterns in the area (e.g., toward SW Sardinia). The uncertainty stemming from limited source parameter knowledge just after the earthquake occurrence (strike, dip, rake, slip distribution) is propagated to the forecast points The passage between PTF results and alert levels can be automated by mapping PTF statistics (e.g., the mean) into the reference intensity intervals (in this case, wave height above sea level). Here, the three alert levels defined within the NEAMTWS are considered and expressed with probability distributions (instead of single tsunami intensity values; Fig. 6b). Similar results can be obtained for any earthquake of any magnitude occurring in the Mediterranean Sea.
The PTF results can be automatically connected to alert levels (Fig. 6b). For example, if we define a reference interval of tsunami intensities for each alert level, the level of each target point can be assigned by verifying in which interval fall some PTF statistics, like, for example, the mean or a pre-defined percentile [68].
Considering the three NEAMTWS levels for the tsunami warning (information, advisory, and watch) and assigning them reference tsunami intensity intervals (e.g. wave heights above sea level < 0.1 m, between 0.1 m and 0.5 m, and > 0.5 m respectively), the performance of alert levels adopting different PTF statistics (mean, median, different percentiles) can be tested for past events, quantifying the implication of different choices for example in terms of rates of missed-and false-alarms [68]. Being connected to risk reduction actions, the definition of the translation rule is not a task connected to tsunami forecasting, but it is connected to specific choices of decisionmakers (e.g., [286]), following a hazard-risk separation principle [157]. Before a PTF-based warning system can become operational, the Italian DPC, supported by the expert advice of its scientific community, would be called upon to define standardized rules to translate PTF's probability distributions into alert levels. Specific rules have not been defined so far.
PTF results seem promising but, before becoming operational, we think that an innovative method like PTF needs to pass through a thorough scientific revision and statistical validation both in hindcasting (retrospective) and forecasting modes, including data from large earthquakes and tsunamis and accurate tsunami records like the ones available from DART data [156]. To this end, to deal with the computational difficulties of extending PTF estimation worldwide, a specific workflow for High-Performance Computing (HPC) machines is under development, to extensively test PTF performances worldwide through HPC urgent computing [67].

Improving civil protection coastal planning
Planned steps for the future include the development of proactive initiatives in coastal areas, following a bottom-up approach to encourage Regions and Municipalities to apply the national Guidelines, drafting/updating their respective emergency plans for tsunami events. This process could take advantage of the promotion, at the national level, of the "International Tsunami Ready Guidelines". These guidelines, whose drafting is being coordinated by IOC/UNESCO in collaboration with experts from the Intergovernmental Coordination Groups of all tsunami warning systems at the global level, foresee a recognition process of the local community who manages to fulfil the 12 key indicators addressed to strengthen local capacities in case of a tsunami.
Given the importance of reaching the population in the shortest possible time, DPC is also working on the implementation of a multi-channel emergency messaging system, which will send messages via cell-broadcast technology, mobile apps, and web services. The system, called It-Alert, will be able to send also SiAM alert messages directly to the population. This channel will be implemented in parallel to the one already operating, which is targeted to all emergency management operators.

Tsunami warning at Stromboli volcano
Stromboli Island is an active volcano located in the southern Tyrrhenian Sea (Fig. 7a) and characterized by persistent Strombolian activity, with persistent low-energy explosions and episodic lava flows and larger explosions. Flank eruptions, deformations during effusive eruptions and paroxysmal explosions may trigger large mass failures along the Sciara del Fuoco (SdF), the most unstable sector of the volcanic edifice. Massive landslides (e.g., in the order 10 6 -10 7 m 3 ) at Stromboli can potentially trigger tsunamis, possibly impacting not only Stromboli but also the other Aeolian islands (in particular the closest, Panarea, 20 km SSW of Stromboli) and even mainland Italy [34,35,40]. From the early twentieth century, six tsunamis generated by landslides within the SdF have been documented [36], and the latest relevant one occurred on 30  [41,287,288]. The waves, which reached an altitude of about 10 m above sea level (run-up), caused significant damage to the buildings located near the beaches of Stromboli, and also reached the near Panarea Island and the coasts of the Southern Tyrrhenian Sea [40,41,43]. Fortunately, no fatalities occurred, mainly due to the fact that most of the affected houses were empty in winter-time and the few locals were alerted by the noise produced by the arriving tsunamis [42].
Landslides and tsunamis triggered by paroxysms usually occur a few minutes after the explosion. Such explosions typically do not have precursors in the hours or days preceding the event in terms of changes in the parameters recorded by the monitoring networks, but deformation of the edifice in the very short-term (minutes) [34,[289][290][291][292][293][294]. Tsunamis may be triggered also during effusive phases, for which ground deformation is crucial for the early detection of magma intrusions and consequent flank instability, as part of a slow and progressive dynamics that could be detected even days before it occurs [290,291,295,296].

Hazard quantification
Proper probabilistic hazard quantification for the tsunamis generated at the SdF of Stromboli is nowadays impossible, since (i) the tsunami source conditions are difficult to predict, observe, and model; (ii) the source is close to the shoreline; (iii) the bathymetry is very steep and irregular [288,297,298]. Therefore, in the past years, research has aimed to identify the most effective modelling approach (in terms of accuracy and computational efficiency) to inform the early warning system and to build a solid basis for future tsunami hazard assessments.
To this aim, several benchmarked models, at different levels of approximation and accuracy, have been used: the NHWAVE three-dimensional non-hydrostatic model in sigma-coordinates [215] and the HySEA family of geophysical codes [299] based on either single layer, two-layer [300] stratified systems or multilayer non-hydrostatic formulations [301] of the wave model.
Models results have been compared for the maximum run-up, inundation maps at the Stromboli village, and the waveform sampled at four proximal sites (two of them corresponding to the locations of the monitoring gauges offshore the SdF, see Sect. 4.2). Both rigid and deformable submarine landslide models (with granular flow rheology described by the Savage and Hutter model [302] and the Pouliquen and Forterre friction law [303]), with volumes ranging from 6 to 20 million cubic meters, have been used to trigger the water waves.
The comparisons between hydrostatic and non-hydrostatic models show that the simulated inundation maps at the Stromboli village are relatively similar, despite differences likely associated with differences in wave-breaking effects. On the other hand, as expected (i.e., due to frequency dispersion and other model differences) preliminary results indicate strong differences between the proximal waveforms produced by the considered models. Such an aspect is particularly critical in the perspective of using proximal sea-level gauges installed at Stromboli to characterize the tsunami source in an early-warning system (discussed in Sect. 4.2). In addition, we stress that subaerial landslides almost invariably generate dispersive waves (e.g. [16]).
This suggests that the use of non-hydrostatic models allows a better description of the proximal waveforms, but the source description remains the most sensitive and uncertain aspect of the modelling. On the other hand, hydrostatic models are computationally less expensive when it comes to modelling larger domains and hence the propagation of the waves far from the source [214]. However, the use of High-Performance Computing (HPC) techniques and of Graphics Processing Unit (GPU) accelerators for the numerical solution [304] has allowed much faster simulations even with complex non-hydrostatic models, which have opened new possibilities for probabilistic tsunami hazard quantification and tsunami early warning [67].
Based on benchmark study results, waveforms obtained with a non-hydrostatic model (Fig. 7b, c) have been selected to preliminary calibrate the tsunami alert system on the Stromboli island, as discussed below. In addition, a set of synthetic scenarios has been created to investigate the potential impact of a landslide-induced tsunami along the shoreline of Stromboli Island, with varying landslide volume, the initial height of the centre of mass, and the landslide-water density ratio.

Tsunami early-warning system at Stromboli
The tsunami early-warning system at Stromboli is based on three main components (outlined in separate subsections below) [150, 305]: 1. The multi-parametric monitoring network, designed to identify anomalies indicating the transition among different eruptive phases, e.g. from explosive (equilibrium) to effusive phase (disequilibrium); 2. Two elastic beacons located offshore in front of the SdF (Punta Labronzo and Punta dei Corvi, see Fig. 8), which are able to identify tsunami waves in real-time, and the data transmission infrastructure to the civil protection's Advanced Operational Centre (COA); 3. An acoustic warning system made of eight sirens distributed in Stromboli, Ginostra (village of Stromboli), Panarea Islands and one beeper inside the Coast Guard premises in Milazzo (Sicily island), interconnected by a dedicated VHF radio network [305].
The system has been operative since 2003, but several significant changes have been implemented through time, as discussed in detail below.

The multi-parametric monitoring network and tsunami warning
Although the mechanism of triggering tsunami waves at Stromboli is not well known, the literature shows that these phenomena occurred more frequently during periods of intense volcanic activity characterized by effusive phases, as well as a few minutes after the occurrence of paroxysmal explosions.
The higher supply of magma to the shallow reservoir is recorded by the multiparametric monitoring network (seismo-acoustic stations, thermal cameras, tiltmeters) months before the effusive onset and is responsible for the progressive transition towards a higher explosive regime with respect to the usual Strombolian activity [291,295,296].
When the dike reaches the surface, a new effusive vent opens, and lava is drained out of the shallow conduit system. During this stage, there is the highest probability of having a flank instability and potential generation of tsunamis (e.g., 2002). Ground deformation based on tiltmeter [289,291] and SAR interferometry (GB-InSAR) [290] of the SdF flank dynamics are therefore crucial to identify early stages of magma intrusions and assess instability scenarios.
According to the civil protection procedures implemented on the island and the overall hazard assessment, these slow and progressive dynamics could allow the civil protection authorities to issue a pre-warning, which could include preparation for the  20 19). c-e Location of the MEDA and the three networks that send data to the Advanced Operational Centre (COA) evacuation of people living in the tsunami inundation zone (e.g., 2007 eruptive crisis described in [295,306]). Furthermore, the activation of the acoustic warning system is automatic and triggered by the elastic beacons after the detection of the tsunami.
The direct consequence of the transition to the effusive regime is the progressive collapse of the crater terrace induced by a large amount of drained magma from the shallow portion of the conduit system [295,296]. Once the effusive phase has started and is efficient, hazards related to flank instability and to potential tsunamis decrease significantly [291].
The increasing hazard posed by volcanic activity (transition from explosive to effusive phase) is assessed through the adoption of a Volcanic Alert Level system based on four-color codes (from green to red). Depending on the alert level and on the overall risk assessment, local and national civil protection authorities implement appropriate actions to ensure public and private safety.

The elastic beacons (MEDA)
The University of Firenze-Laboratorio di Geofisica Sperimentale deployed two elastic beacons (MEDA) used to detect the tsunami proximal to the SdF. The first one (PDC) was installed in 2008, 260 m off the coast of Punta dei Corvi, while the second one (PLB) was set in 2017, at a distance of 350 m from Punta Labronzo (Fig. 8a).
The MEDA is a semi-rigid structure made of a 30 m long metallic pole, anchored to the shoal with a 20 Tons deadweight. The deadweight connects to the upper structure through an anti-torsion steel cable (Fig. 8b). The MEDA reaches 8 m of height above the sea surface and is equipped with sensors, power supply equipment, and a radio transmission system [150]. Sensors include hydrostatic pressure (125 Hz), temperature (1 Hz), hydroacoustic (40 Hz), GPS, and two tiltmeters (4 Hz), combining measures of interest for the tsunami detection and measures to monitor operability and performance. The MEDAs are installed 260 m off Punta dei Corvi and 350 m off Punta Labronzo (Fig. 8), have a length of 30 m, a height above sea level of 8 m, and a deadweight of 20 Tons.
Three different networks send data to the COA, as shown in Fig. 8c-e, in order to guarantee data transmission and redundancy.

Tsunami detection
Tsunami automatic detection algorithm [150,307] has been tested adopting two main criteria: (1) surface waves dispersion and (2) STA/LTA (Short Time Average over Long Time Average) ratio analysis.
For the surface wave dispersion, underwater hydrostatic pressure sensors (one at a depth of 50 m and the second on the seafloor) are used to make the sea-noise effect negligible and to increase tsunami wave detection capabilities. For the STA/LTA method, which is largely used in seismology for the automatic detection of earthquakes, we calculate the average absolute value of the signal's amplitude based on two-time windows.
Sensitivity analysis of the algorithm and of the threshold efficiency has been performed for at least five years of recorded data from PDC MEDA, using the database of simulations performed for the hazard quantification (Sect. 4.1) and the real-time measures (for non-tsunami signals). As a result, the algorithm is theoretically able to automatically detect tsunami waves in any sea condition since tsunami signals produce STA/LTA values well above the identified threshold [307]. The first real-time testing occurred during the small tsunamis that occurred in July and August 2019, as discussed below.

The acoustic warning system
The acoustic warning system is made of 8 sirens and one beeper (Fig. 9), interconnected by a dedicated radio network (VHF band), and three main base stations (radio links) located in Antennamare (in the province of Messina), Stromboli, and Panarea, that allow the simultaneous activation of the above-mentioned sirens [305].
Management software that runs on a workstation located inside the COA gives the opportunity to select three possible combinations: • Single: activation of one siren (out of the nine available); • Group: activation of a pre-selected group of sirens (for example, only Stromboli group or only Panarea group); • Global: activation of all the sirens.

Evacuation zones and civil protection coastal planning
The national civil protection emergency plan for volcanic events of Stromboli was first issued in 2003, after the 30 December 2002 tsunami caused by the sudden failure of the SdF. The plan was drawn up with the contribution of the Regional Civil Protection, the Prefecture-UTG of Messina and the Municipality of Lipari. In August 2015, a new updated version was released (http://www.protezionecivile.gov.it/attivita-rischi/ rischio-vulcanico/vulcani-italia/stromboli/pianificazione, [305]), taking into account new scenarios and the introduction of the national volcanic alert level system.
The plan is based on a tsunami inundation zone drawn upon the in-situ observation and measurements of the 30 December 2002 waves [40][41][42]. This choice is determined by the fact that this event is, to date, the best documented in terms of impact on the island, as well as the largest occurred in recent times. In terms of emergency management, this area is densely inhabited and it represents the exposed area that will be evacuated in case of (i) evidence of magma intrusions and deformations in the SdF, potentially leading to a collapse and subsequent tsunami (e.g., 2007 eruptive crisis described in [306]); (ii) automatic activation of the acoustic warning system (sirens) due to the detection of a tsunami. Noteworthy, also an automatic early warning system for paroxysmal explosions is in place [294]. This system triggers sirens locally at the villages of the island of Stromboli (Stromboli and Ginostra) [305]. The tsunami early warning system is instead automatically triggered only when a tsunami is detected by the MEDAs, with a global activation. For these two cases, the sirens have two different acoustic signals: (i) bi-tonal for paroxysmal explosions and (ii) mono-tonal for tsunamis.
Evacuation routes have been identified to guide people to the safe areas and reach the gathering points in case of tsunami, identified by emergency signals explaining the correct path to follow. More details can be found in [305].

Modelling of the 2019 events
On July the 3rd 2019, a paroxysmal eruption at Stromboli produced a sequence of highenergy pyroclastic flows, which generated several tsunami wave trains detected by the sea gauges (see Sect. 4.4.2). In response to the event, a preliminary assessment of the physical parameters characterizing the tsunami source was done by first comparing the recorded waveforms with the pre-computed synthetic waveforms and then running additional simulations with the multilayer-HySEA non-hydrostatic model.
Either one or two simultaneous pyroclastic flows were modelled starting from four different initial positions and using three different volumes. Simulated flows average speeds were compatible with that observed by webcams. Two possible scenarios were therefore taken into consideration: in one scenario, a single pyroclastic flow dominated the wave generation; in the other one, two pyroclastic flows contributed in a comparable way to wave generation.
The preliminary results of the numerical modelling allowed us to estimate, based on the wave heights observed at the wave gauges (see Sect. 2.2.2), the rock avalanche volumes and generation height. The simulated waveforms and arrival times appear to have characteristics very similar to those measured in the event.
The study of this event allowed comparing, for the first time in Stromboli, the wave generation and propagation model with the signals measured in the presence of a tsunamigenic event. The model has proven to be capable of reproducing the main characteristics of the waves in the two gauges, with satisfactory accuracy, at least in the first minutes of evolution.

Warning system
The acoustic warning system for tsunamis has been working in manual mode since its installation in 2003, following the 30 December 2002 landslide and tsunami. It was first activated on 27 February 2007 during the phases that preceded the lateral vent opening on the SdF, fearing a new flank collapse and subsequent tsunami [306]. Fortunately, the vent opened, and just a minor landslide occurred with no tsunami observed.
After the deployment of the two MEDAs, the DPC and the University of Firenze-Laboratorio of Geofisica Sperimentale started working on an automated procedure to activate the sirens through a radio-modem (working on the same frequency of the acoustic sirens) connected to the server running the detection algorithm [150].
On 3 April 2019, the first successful test was performed, and one single siren was remotely activated. The testing phase was still on-going when, on 3 July 2019, a paroxysmal explosion occurred, and the two elastic beacons recorded a small tsunami (1 m wave peak-to-through, Fig. 10a) generated by pyroclastic flows. This tsunami had no significant impact on the coast. The event was sudden, and the monitoring equipment located on the volcano was not able to detect any precursors, so there was no chance of activating the sirens manually.
On 28 August 2019, a new paroxysmal explosion was recorded, but this time local civil protection personnel deployed on the island were able to activate the sirens in manual mode as soon as the pyroclastic density currents reached the sea. In this case, a minor tsunami was recorded (60 cm peak-to-through, Fig. 10b), but with negligible impact on the coasts. In this event, the pyroclastic flow entered the sea approximately 35 after the volcanic eruption. The tsunami was generated a few seconds later, approximately 60 after the occurrence of the volcanic explosion (after 60 ).
In Fig. 10c, we report, as an example, the automatic tsunami detection for this event, based on the STA/LTA ratio [307]. The automatic trigger overcame the predefined threshold approximately 15 second after the generation of the tsunami, that is, 1 and 15 after the volcanic explosion. As far as we know, this is the first time a warning is delivered before the tsunami wave is fully developed.
The testing phase of the automated procedure ended on 9 September 2019 when the Mayor of Lipari, Regional and National civil protection authorities tested the entire chain of the procedure, simulating a trigger signal from the MEDAs.
Therefore, at this stage, if the MEDAs record a tsunami, the acoustic warning system is automatically activated (global mode) by the server that runs the tsunami detection algorithm (Fig. 9b) [150]. Should the automated system be down for any reason, there is the chance to manually activate the sirens from the workstation inside the COA or from three portable radios, which have been given to the Municipality of Lipari. Fig. 11 a Portion of the Euro-Mediterranean Tsunami Catalogue [36] involving historical events generated by offshore or coastal/inland earthquakes. The tsunami intensity is expressed in the Sieberg-Ambraseys scale. b Portion of the Euro-Mediterranean Tsunami Catalogue [36] involving historical events generated by non-seismic sources. The tsunami intensity is expressed in the Sieberg-Ambraseys scale

Other non-seismic tsunami sources in Italy
In this section, non-seismic tsunami sources are considered as a whole set of known and potential tsunamigenic sources, including landslides (in its widest meaning), volcanic activity, and atmospheric pressure disturbances. The portion of the Euro-Mediterranean Tsunami Catalogue (EMTC; [36]) describing the sources of the historical tsunamis hitting the Italian coastlines (Fig. 11) indicates that, although tectonic sources are the primary triggering mechanism (Fig. 11a, 49 entries), non-seismic sources are a non-negligible part of the reported events (Fig. 11b, 23 entries).
Looking at Fig. 11b, it is easily recognizable that historical tsunamis generated by sources other than earthquakes concentrate in volcanic areas (Gulf of Naples, Aeolian Archipelago, offshore Mount Etna) and along steep coastal slopes where pure gravitational instability of seismic shaking can trigger mass transport phenomena. We stress again that these can both take place as submarine (fully submerged) landslides or subaerial landslides plunging into the sea (e.g., the 1783 Scilla landslide; [308,309]). No mention is made in the EMTC to meteotsunami events, although some of these may fall under the events reported with "unknown" cause.
Regarding earthquake sources, the EMTC is considered complete starting approximately from the seventeenth century. However, it is very reasonable to argue that the EMTC is significantly incomplete as regards non-seismic sources due to their relatively more local impact and the fact that such sources were previously not equally well understood as potential triggers. Some of the most destructive historical tsunamis hitting the Italian coasts, such as the 11th January 1693 and the 28 December 1908 events, are traditionally associated with the large magnitude earthquakes occurred a few minutes before the tsunami, but studies published in the last decade have introduced and In the perspective of tsunami hazard assessment and warning related to non-seismic generation mechanisms, improved and deeper knowledge of historical events must be accompanied by systematic and detailed characterization of all potential sources. Significant progress in this sense was achieved by the MAGIC project, financed by DPC in the period 2007-2012 (http://www.protezionecivile.gov.it/media-communication/ dossier/detail/-/asset_publisher/default/content/progetto-magic). Many of the results matured in the framework of MAGIC were published in the scientific literature, allowing to get a significantly improved picture of the coastal areas where careful investigation on slope instability and possible future tsunamigenesis must be carried out.
In Figs. 12 and 13, with no claims to be exhaustive, we summarize the main geographical areas around Italy where tsunami generation by non-seismic sources occurred in the past, likely representing a threat also in the future. The main identified potential sources are discussed below.

Ligurian Sea
Only two historical tsunamis are present in EMTC for the Ligurian Sea (LS) with a non-seismic source generation. They occurred respectively on July 2nd, 1703, and on July 3rd, 1809, and both were moderate-intensity event (intensity 2 in the Sieberg-Ambraseys scale). Both events are labelled with an "unknown" source, which is an indicator of the lack of a thorough analysis of the tsunami hazard from non-seismic sources in this area.
Studies on offshore slope instability and on mass transport evidence along the Ligurian margin indeed exist (e.g., [310][311][312]), but a systematic approach to tsunami hazard and warning related to this domain appears to be still lacking, especially as regards the Italian coasts. It is important to recall that landslide-generated tsunamis can pose a serious threat for the coasts facing the Ligurian Sea, as demonstrated by the 1979 event in Nice (e.g., [313,314]).

Ischia and the Gulf of Naples (IGN)
The Ischia island is characterized by an intense volcanic activity whose main expression is the resurgence of Mt. Epomeo, causing instability phenomena along the slopes of the main relief due to both seismic shaking and increased slope gradients [232,315,316]. Several bathymetric campaigns highlighted the existence of slide deposits offshore the northern and western coasts of Ischia (e.g., [317]). Moreover, a debris avalanche offshore the southern coast of the island was mapped at about 1000 m depth, with a volume of 1.5-3.0 km 3 , composed of blocks (hummocks) and mixed and finer facies, up to 45-km distance from the coast [28,318]. An extreme scenario of tsunami generation by a single collapse involving the island top (about 800 m a.s.l.) down to 700 m depth b.s.l., with a total volume of more than 3 km 3 , was simulated by Tinti et al. [27]. The resulting simulated tsunami is catastrophic along the coasts of Ischia and Capri, with waves as high as 40 m; several meter-high waves are also expected inside the Gulf of Naples and in several locations to the north (Latium) and to the south (southern Campania and northern Calabria).
Another potential source of catastrophic collapse and consequent tsunami generation at Ischia is represented by the Monte Nuovo, which is a deep-seated block found at 400 m a.s.l. on the north-western flank of Mt. Epomeo (e.g., [232,319]). Zaniboni et al. [29] investigated the potential tsunami generation by a 160 × 10 6 m 3 rock collapse and its propagation in the Gulf of Naples area. Paparo and Tinti [320] performed a stability analysis for Monte Nuovo and found that it could be mobilized if an earthquake with magnitude comparable to the 1883 Casamicciola earthquake occurred on a nearby fault: the consequence of the destabilization would be a catastrophic mass failure with the potential to generate a disastrous tsunami. However, this possibility is still debated in the literature [232].
The activity of Mt. Vesuvius is responsible for the largest part of the historical tsunamis reported in EMTC: the two most famous and catastrophic eruptions occurred on 24 August 79 AD and on 17 December 1631 are both reported to have been accompanied by significant sea movements. A short review of the numerical simulations of pyroclastic flow impact and consequent tsunami generation in the Gulf of Naples can be found in [46,98,321].
Recently, the potential for generating tsunamis due to underwater explosions at the Campi Flegrei caldera has also been explored by Paris et al. [32]. Even in the absence of a past record of tsunamis generated by this mechanism, underwater vent opening cannot be excluded as possible tsunami trigger in the area (e.g., [322,323]), and they have the potential for large tsunamis impacting mainly within the Gulf of Naples. Submarine explosions appear to be the main potential generating mechanisms at Campi Flegrei because pyroclastic flows are highly dilute and, consequently, they have a low potential for generating tsunamis [20,324].

Southern Tyrrhenian Sea Mounts (STSM)
The complex tectonic setting related to the Tyrrhenian Sea evolution (Ionian slab moving towards NW and subducting underneath the Tyrrhenian crust) finds expression in a series of volcanic features that characterize its morphology. While the Aeolian Archipelago represents a typical volcanic-arc configuration, the back-arc expansion zone shows submarine volcanic structures, which were generated in different periods but that exhibit, in some cases, similar peculiarities.
The western one, named Vavilov seamount, is the oldest and presents an asymmetric morphology [325], suggesting the possible collapse of its western flank. The eastern slope average dip is about 15°, while the western one is higher (25°) and with a typical amphitheatre-like shape, supporting the hypothesis of a singular ancient catastrophic event [326], with volume ranging 40 km 3 . A first attempt to assess the consequences of such an event in terms of tsunami generation has been conducted by Arcangeli [327], with the use of numerical codes. The supposedly generated wave would impact the whole Tyrrhenian coasts with relevant waves.
The eastern seamount characterizing the area, named Marsili, is the highest volcanic edifice in the whole Mediterranean Sea, with dimensions 70 km × 20 km and height from the sea bottom exceeding 3000 m [328]. Even though no significant pieces of evidence of instability exist [329], there are some elements that justify the study concerning the potential collapse of portions of the Marsili volcanic edifice: • the presence of tectonic structures along its slopes, possibly contributing to volcano flank instability; • similarity, both in genesis and morphology, with the Vavilov seamount, which already was affected by a gigantic lateral collapse; • the presence, on the abyssal plain at the basis of the volcano, of deposits associated with mass-wasting processes originating from the volcano flanks.
The generation of tsunamis from mass collapses has been scarcely investigated yet. One first attempt to quantify the waves induced by landslide scenarios of different dimensions and depths has been conducted by Gallotti et al. [330].
Moving further to the east, close to the Campanian coast, recent studies have described the Palinuro volcanic chain, which displays a W-E trend [331] and is characterized by several volcanic cones and intense activity that can induce flank instability. The involved volumes range from 0.5 to a few km 3 [332]. This fact, combined with the relatively shallow water of the area and the steep slope of the southern margin, contributes to enhancing considerably the potential hazard associated with this area.
The other existing emerged (Ustica) or submerged (Enarete, Eolo, Alcione, Lamentini, Glabro, Enotrio, Ovidio) volcanic edifices in the southern sector of the Tyrrhenian Sea, in spite of some record of submarine gravitational fluxes, do not seem to present the conditions for a significant tsunami threat [333].

Aeolian Archipelago (AA)
Apart from the activity on the SdF, thoroughly discussed in Sect. 4, other tsunamigenic sources at Stromboli cannot be excluded, like sector collapses in the south-eastern zone (Sciara Vecchia, Schicciole, Rina Grande) or submarine eruptions, especially in the area between Punta Lena and Strombolicchio [333]. Also, other volcanic phenomena generated above the sea level and propagating into the sea in sectors different from the SdF cannot be excluded. However, the potential tsunami hazard due to these phenomena has never been quantitatively assessed.
The island of Volcano is another island of the Aeolian Archipelago where historical tsunamis occurred and where a potential for future events exists [26,231]. On April 20th, 1988, a 2 × 10 5 m 3 body [26,334,335] detached from the La Fossa crater and collapsed into the sea, generating a small local tsunami hitting Porto di Levante with a 1 m amplitude wave. The landslide dynamics and the ensuing tsunami have been simulated numerically by Tinti et al. [26]. It is important to stress here that the La Fossa crater is characterized by frequent sliding phenomena (e.g. [336,337]), and hence it represents a potential repeating tsunamigenic source. Apart from gravitational sources, even though no records of eruption-generated tsunamis are available (either submarine eruptions or volcanic phenomena entering the sea), this cause cannot be ruled out [231].
In addition to Stromboli and Vulcano, other gravitational instabilities potentially generating significant tsunamis have been tracked in Lipari, especially in its eastern sector and both in the emerged and the submerged parts of the volcanic edifice, and, with relatively smaller tsunamigenic potential, in Panarea [333].

Western Calabria Margin (WCM)
Tyrrhenian Calabria has been affected by several large historical earthquakes, often accompanied by moderate to destructive tsunamis [308]. Examples are the events that occurred in 1184 (M w 6.8), 1638 (M w 7.1), 1783 (a long sequence with three M w > 6.7 earthquakes in two months, M w 7.1 being the largest), 1905 (M w 7.0), and 1908 (M w 7.1). For some of them (mainly 1783, but also 1908), a significant contribution may be ascribed to seismically induced landslides (for the 1783 event see [338][339][340][341]; for the 1905 event see [342]).
Moderate-to-large landslides have also been mapped by offshore surveys (e.g., the late-Quaternary landslide mapped offshore Cape Licosa and in the Paola basin) and represent potential tsunamigenic sources (e.g. [343]); moreover, the Scilla landslide, triggered by the 6th February 1783 earthquake, is a historical example of a landslidegenerated tsunami hitting the Tyrrhenian coasts of Calabria violently. On that date, a rockslide detached from Mt. Campallà, entering the sea and generating a tsunami which impacted the nearby village of Scilla with waves up to 9 m high, killing around 1500 people. The simulation of the rockslide dynamics and of the ensuing tsunami has been addressed by several authors [338][339][340][341].
The Gioia Tauro basin is another sector that deserves special attention. It is characterized by a canyon [344] where past mass transport phenomena have been mapped and described. On 12 July 1977, a tsunami was observed with an initial withdrawal of 2-3 m, subsequently attacking the western dock of the Gioia Tauro harbour with

Eastern Calabria Margin (ECM)
The Assi landslide has been studied and characterized in the frame of the MAGIC project [346]. The main mass transport evidence is found offshore the Monasterace town and has been dated at about 2850 years ago, possibly resulting from two distinct events. The first of them involves a 1.85 km 3 volume, with an average thickness of 27 m. The landslide supposedly detached from a depth of 750 m b.s.l. and the deposit is found at 1400-1600 m depth, resulting in a run-out of the order of 20 km. The landslide has been simulated numerically [347]: the main result is that a landslide of this type could be responsible for waves with maximum amplitudes in the order of 1 m along a limited stretch of coast between Monasterace and Roccella Jonica. Although not catastrophic, such a tsunami might cause significant damage to small harbours along the coast and could also result in a threat to the local population, especially during the touristic seasons. In addition, this type of landslide poses a significant challenge from the warning point of view, as the first wave is likely to impact the closest coastlines nearly 7 min after the slide onset.
A second highly debated area is the so-called Crotone Swell offshore the Crotone peninsula (e.g. [348]), which is a morphological high approximately 16 km long and 30 km wide, exhibiting a very complicated structure. Zecchin et al. [348] interpret the Crotone Swell as a mega-landslide, up to 1.6 km thick, which started moving between Late Zanclean (5.3-3.6 million years ago) and Early Piacenzian (3.6-2.6 million years ago). According to the same authors, the seaward migration of such a mega-landslide significantly slowed down in the last 0.5 million years and currently its behaviour is not expected to pose any significant threat as a whole. Nonetheless, the complex and the not fully-understood structure of the landslide deserves a careful evaluation of the possibility that some portions, particularly the shallower parts of the Crotone Swell, might be set in motion by the seismic activity of the area, which is characterized by historical earthquakes with magnitude larger than 6.5 (08/06/1638 and 08/03/1832). For example, a seismic profile analysed in the framework of the recent (2017-2018) SPOT Project [349], financed by the Italian Ministry of Economic Development with the technical support of the Italian DPC, has highlighted the possible presence of a huge submarine landslide in the southern portion of the swell (offshore Capo Rizzuto). A possible reconstruction involves a 25 km 3 slump, covering an area of 85 km 2 with an average thickness of 300 m. Even with a limited run-out to its mainly rotational movement, it is reasonable to expect that this source could have been responsible for a very significant tsunami hitting the entire Ionian coasts of Calabria.
We have commented on only two main areas that have been extensively described in the literature. Other areas along the eastern Calabria margin deserve further investigation and careful assessment of the tsunami hazard related to coastal and offshore slope instabilities.

Eastern Sicily (ES): Messina Straits
The Messina Straits is an extremely complex area with very peculiar morphological features. As described recently by Ridente et al. [350], the main morpho-structural feature is represented by the Messina Canyon, striking mainly NNE-SSW and intercepting many tributary slope canyons that incise the steep Calabrian and Sicilian continental slopes. The landslides triggered by a cloudburst on 1 October 2009 in the area of Giampilieri and Roccalumera (Sicilian side of the Straits), resulting in 28 ascertained victims, and the intense debate on the source of the 1908 tsunami following the already mentioned devastating M 7.1 earthquake are the demonstration of the importance of a proper assessment of the susceptibility of the area to coastal and submarine mass movements. Regarding the debate on the 1908 tsunami, we will recall here that in the lack of any recognized primary faulting evidence ( [351], with references), many hypotheses on possible underwater slide positions and mechanisms have been proposed, starting from [47], with further contributions including the studies by Favalli et al. [50] and Schambach et al. [256]. In particular, in this very recent study, a key role in the generation of the 1908 tsunami is assigned to a ∼ 2 km 3 submarine mass failure mapped at the foot of the Fiumefreddo Valley, NE of Mount Etna. For the scope of the present discussion, we stress in particular what was proposed by Ridente et al. [350]: based on the results of high-resolution swath bathymetry, they pointed out that several relatively small-scale landslides are observed on both the Calabrian and Sicilian margins of the Straits, possibly related to the shaking associated with past earthquakes such as the relatively recent 1908 event. From the tsunami hazard assessment point of view, this must be carefully kept in mind as even small-scale landslides might determine local amplification of tsunami inundation effects along the coasts of the Straits.

Eastern Sicily (ES): Mount Etna
The area offshore Mount Etna has been intensively investigated by marine geology campaigns in the past, and several descriptions and interpretations about mass transport deposits can be found in the literature. A full review of the available studies is beyond the scope of this paper, but it is important to recall that Mount Etna is affected by sliding of its eastern flank, whose nature is the subject of intense debate (e.g. [352]). The recognized instability of this flank ( [352][353][354][355][356], among many others) and its ongoing deformation represent a critical subject in the frame of tsunami hazard and warning, at least for the eastern Sicily coastlines.
No unique view on the history of mass transport phenomena exists: the spectrum of hypotheses ranges from the small-scale deposits described in [354], reasonably triggered by volcano induced seismicity and flank deformation connected to the eruptive activity, to the catastrophic debris-avalanche deposit discussed in [357], related to the Valle del Bove scar and characterized by a 16-21 km 3 volume and up to 20 km offshore run-outs. A careful revision and a deeper understanding of the ongoing processes that might lead to future offshore mass movements are needed for a credible tsunami hazard assessment in this area.

Eastern Sicily (ES): Hyblean-Malta Escarpment
The Hyblean-Malta Escarpment (HME) is a geomorphological structure characterized by large bathymetric gradients. In principle, this is an important prerequisite to consider the HME as a place where mass collapses may occur, due either to pure gravitational instability or to earthquake shaking. Nonetheless, very few examples exist in the literature of documented slide deposits, or of scarps suggesting the occurrence of past slides, or of potentially unstable mass bodies.
The only clear evidence is reported in [358], where a 5 km 3 deposit mapped at 2000 m depth, offshore just north of the bay of Augusta, is described. It is argued that the main movement of the mapped mass was rotational, hence describing a slump, with an inferred 200 m dislocation over a 40 km 2 area, involving an average thickness of 100 m. Tonini et al. [359] and Argnani et al. [358] discussed the possible triggering of the slump by the 11th January 1693 earthquake, which was followed by a violent tsunami whose most dramatic effects were experienced by the town of Augusta, with observed wave amplitudes of more than 8 m [36,360]. Based on the results of numerical tsunami simulations, Tonini et al. [359] and Argnani et al. [358] agree that that a slump could not be the main responsible for the observed tsunami, hence leaving open the problem of the source of the 1693 tsunami.
More recently, Paparo et al. [361] studied the slope stability of the HME under seismic loading on a number of transects, trying to identify those sectors that may move under the effect of seismic shaking. One of the results of the study was the estimation of landslide volumes that might be set in motion by seismic shaking, representing a fundamental step toward quantitative tsunami hazard assessment for the nearby coastal areas.

Gela Margin (GM) and Sicily Channel (SC)
The Gela Basin Margin (GBM) is located in the Sicily Channel, south of central Sicily. This structure connects the shallow water characterizing the Malta Plateau (around 200 m) to the east, with a depression in deeper water (800 m depth b.s.l.). The whole area was the object of geophysical surveys, describing many features possibly related to mass-wasting processes: scars and canyons along the slope, huge amounts of sediments at its toe [362].
The most noticeable feature, in terms of involved volume, is the Gela slide [363], related to a prehistoric event (~600 thousand yeats before present) with a volume estimated at 630 km 3 . Such submarine collapse has been adopted as one of the scenarios for the study of the tsunami hazard involving the Malta archipelago [364] and has been found to be the major source of potentially catastrophic waves impacting the coasts.
More recent surveys have characterized further episodes along the escarpment. This is the case of the so-called Twin Slides, related to a double event placed in the northern portion of the GBM. Two evident scars develop just from the continental platform at 200 m depth b.s.l., 30 km offshore the city of Gela, with the corresponding deposit laying at their toes [365]. These two failures are the final phase of a complex series of collapses, tentatively reconstructed by Kuhlmann et al. [366], that repeatedly occurred in the last 10 thousand years. The Twin Slides volumes are estimated at about 0.5 km 3 each; the Northernmost one has been selected for scenario reconstruction and tsunami simulation [367]. The generated tsunami affects 60 km of southern Sicily coast with waves of at least 1 m maximum, with peaks of 3-4 m.
A larger scar, with the correspondent deposit at its toe, has been discovered on the southern part of GBM, 40 km offshore the Maltan island of Gozo. This landslide, whose volume has been estimated more than 3 km 3 , has been simulated numerically and further been coupled to a tsunami model obtaining large tsunami intensities (maximum wave height reaching 6 m [367]).
No historical tsunami generated by non-seismic sources is present in the Sicily Channel, according to EMTC. Despite this fact, a careful approach to tsunami hazard in the area must take into consideration at least the submarine volcanoes mapped in the NW sector of the Channel. For the Graham, Terrible, and Nerita banks (e.g. [368,369]), an important historical record of submarine eruption exists (e.g., the formation of the Ferdinandea island in 1831). While the potential for tsunamigenic gravitational instability seems rather low, eruptive activity potentially generating tsunamis cannot be excluded [333]. Other volcanic edifices are the Pelagie islands [370,371]. While tsunamigenic activity for these volcanoes has never been quantitatively studied, in Pantelleria and Linosa, the generation of tsunamigenic gravitational collapses and of potentially tsunamigenic submarine eruptions cannot be ruled out [333].

Adriatic Sea
The main potential tsunami source in the Adriatic Sea is the South Adriatic Margin (SAM). The SAM is a peculiar structure, interrupting the gentle slopes and shallow water characterizing the remaining part of the Adriatic Sea and presenting many features that can be linked to mass mobilization and rapid failures. In particular, the western flank, offshore Apulia, is characterized by extensive evidence of slope instability for at least 150 km, especially on the shelf edge and upper slope [372,373].
The main features of the Western SAM are: -the Bari Canyon System, characterized by intense bottom currents causing intense sediment deposition and erosion [374]; -the Vieste slide, on the north, involving low-gradient areas [375]; -the Gondola slide, which is the main body characterizing the margin, probably related to an event occurring simultaneously to that of the Vieste slide, due to the shaking provoked by a big earthquake during the Last Glacial Maximum (approximately 62 thousand years before present; [376]). Its failure mechanism is quite complex, mainly divided into three phases, and involves around 4.5 km 3 of material, with a run-out of more than 50 km for the looser portion. A pronounced scar (10 km wide, 250 m deep) cuts the platform at about 150 m depth, clearly marking the landslide detachment zone.
An attempt to assess the tsunami generation of the Gondola slide through numerical codes has been done by Tinti et al. [378]: in the simulation, the sea level has been lowered by more than 100 m, accounting for the likely coeval water depth. The ensuing waves reached 10-20 m on the Apulian coasts.
A historical tsunami in this area is also documented in association with the Mw 6.7 1627 Capitanata earthquake (CFTI5Med; [45]). The epicentre and the causative fault of this earthquake are located onshore [44,377]; therefore, the source of the related tsunami is still unknown.

Meteotsunamis
Several cases of meteotsunamis are documented throughout the Mediterranean, as in the cases described in Fig. 13a, which is taken from [379]. By investigating the resonance conditions in terms of Froude number (Fr), the authors showed that the most favourable conditions for meteotsunami generation (0.9 < Fr < 1.1) exist along most of the Italian coasts, especially in the Adriatic Sea and in the Strait of Sicily (Fig. 13a).
For example, a recent publication by Šepić et al. [380] describes the occurrence of a meteotsunami in the city of Mazara del Vallo in 2014, producing inundation of about 1 m, channelling upstream along the local river and some damages on boats. Another recent example of meteotsunami along the Italian coasts is reported in Fig. 13b, showing the record of a meteo-tsunami in the Adriatic Sea occurred on 1 August 2019. The inundation in some places may be quite sudden, as in the case of Ancona, where a sea-level rise of almost 1 m has been recorded by the tide gauge located in the harbour.
Attempts to quantify the probability of meteotsunami occurrence exist, based on probabilistic synoptic meteotsunami index [73] that can be potentially extended to integrate nonstationarities and time-dependencies induced by climatic factors (http:// www.savemedcoasts.eu/). These forecasts require the monitoring of synoptic conditions (highlighting the favourable condition for their generation), as well as the need to track small-scale air pressure disturbances (meteotsunami source) and atmosphere-sea interaction (meteotsunami generation and propagation). As pointed out above for seismic sources, these probabilistic methods may represent the future for science-informed tsunami warnings, but currently extensive testing against data is still needed.

Outlook
The critical challenge in common to all atypical (seismic and non-seismic) sources is the need to deal with significant uncertainty due to a large number of different and difficult-to-model potential sources and source mechanisms. The quantification of this uncertainty, which is even larger when there is a need to act fast to increase the lead time for local tsunami warning, is fundamental for tsunami risk management. This is a challenge from the scientific, computational, and decision-making points of view. Significant scientific research is still required to characterize many of the atypical sources in a simple and unified way, provided it is possible, and in many cases it is still needed to develop ad hoc modelling procedures, as well as monitoring (e.g., of the triggering conditions) and warning strategies. Moving away from deterministic previsions to probabilistic forecasts that manage source uncertainty (e.g. [96]), while reducing the uncertainty through the improvement of the source and tsunami monitoring systems and processing algorithms (e.g. [151]), are tasks that should be pursued synergically. Indeed, only within a coherent uncertainty treatment framework, the technologically-and modelling-driven uncertainty-reduction strategies can be exploited without creating dangerous biases due to overconfidence. In fact, uncertainty will never be eliminated, and a certain degree of unpredictability will always need to be managed.
The increasing number and quality of probabilistic hazard and risk analyses show that there are a growing capability and intent to study and characterize complex sources such as the atypical ones, especially the seismic ones. Recent hazard studies on seismically generated tsunamis have already started including non-megathrust events [11,31,52,53], and many efforts are also ongoing to improve the computational efficiency [66,67,69] and to extend toward non-seismic sources [11,17,32,98,191]. These hazard studies may provide important information also for tsunami warning, such as (i) the description of the full variability of the source mechanisms (e.g., [31]); (ii) the prioritization of sources based on the capability of generating tsunamis and the long-term frequency of tsunamigenic events (e.g., based on hazard disaggregation; [31,381]); (iii) the development of databases of pre-defined sources and related tsunami propagation scenarios that may be used in real-time tsunami warning (e.g. [52,53,66,68]). Important contributions may also arise from the development of optimized workflows making use of high-performance computing, like the ones in development within the ChEESE Centre of Excellence project (https://cheese-coe.eu/), providing the opportunity to refine the uncertainty quantification, as well as to improve the management of disasters caused by tsunamis through the post-event characterization evaluated by means of urgent computing [67].
Hazard quantifications for non-seismic sources are still very limited [11,54]. While a few prototypal probabilistic hazard analyses have been discussed in the literature (e.g., [17,30,32,65,98]), most present-day analyses are limited to the identification of reference past scenarios and the search of appropriate modelling tools to define datasets of possible scenarios, like in the discussed case of Stromboli. The lack of recognized modelling procedures and of probabilistic source models prevented until now the development of proper hazard characterization, and a strong research effort in this direction is still needed.
The main pathway for implementing tsunami warnings for atypical sources is the development of suitable source and tsunami monitoring systems, possibly keeping in mind some lessons from the present-day systems for seismic tsunamis. For the volcanic sources, the situation is complicated by the large variety of possible tsunami-generating events. In the case of Stromboli, the frequency of occurrence of tsunami generating events from one specific source area (the SdF slope) pushed toward the development of a specialized system completely devoted to this single source. Specialized systems focusing on only one or a few potential source areas may not be appropriate for many volcanoes. Even assuming that all the relevant source areas can be identified, due to their number, their direct monitoring may be practically impossible. In either case, the only viable solution is the implementation of prioritization efforts, based on a profound uncertainty analysis to allow the decision-makers to define which volcanoes and/or areas of gravitational instability are the most important ones to concentrate the efforts, based on both hazard frequency and intensity, as well as the potential impact due to high vulnerability and exposure. Ideally, it should be first clarified what is the level of residual risk is being accepted in order to establish to which extent different source mechanisms should be considered, even if risk is very difficult to be achieved at present given the large uncertainty associated with these phenomena.
The development of new tsunami monitoring technologies will bring new opportunities for tsunami warning [151,382], reinforcing the source monitoring that is foreseen to be subject as well to dramatic progress, including specific progress related to specific atypical sources such as for example tsunami earthquakes (e.g. [284,383]). The Stromboli warning system shows that source and tsunami monitoring should coexist, as indeed is also the case for the existing warning systems for tsunamigenic landslides in Norway. The current system at Stromboli Island is based on the direct detection of the tsunami that has been optimized by focusing on an exact location as the most probable source (SdF). Nevertheless, it is evident that a more constrained forecasting of the conditions for the tsunami generation can significantly increase the lead time and allow for a more effective evacuation.
The need for joint improvements in the forecasting/monitoring procedures for the triggering event is indeed a general remark valid for all sources, in which any tsunami detection should be at least accompanied by a specific monitoring system focused on the early detection of the triggering phenomena that may facilitate the tsunami generation [10,17,72,150,269,279]. This principle can also be extended to secondary sources like induced landslides, for which statistical methods are able to forecast the events based on the observation of triggering events, like rain or earthquakes. In the future, such forecasts may be coupled with local susceptibility analyses to develop probabilistic forecasts of the tsunami impact. While perhaps they will never reach high probability values in their forecasts, such methods may be used to provide significant probability gains that may help decision-making (e.g., [384]).
Probabilistic forecasting methods for tsunami warning may provide the background to merge the different strategies discussed above [68,95,96,153,154,385]. Indeed, they are based on the explicit consideration of ensembles of sources that virtually cover the whole natural variability, as well as of the uncertainty of tsunami generation and propagation. The definition and the weights of the ensembles can be based on direct monitoring of the tsunami-generating event, as well as by the monitoring of potential triggering events and/or by direct tsunami measures (e.g., through data assimilation techniques; [386][387][388]). The recent developments on the management of uncertainty in the short-term forecast of tsunamis caused by crustal sources show that the shift from empirical decision matrix or deterministic scenarios to computationally based probabilistic tsunami forecast is possible in the near future. The development of similar methods is ongoing also for meteotsunamis (e.g., [73,379]), based on the monitoring of the atmospheric conditions leading to these phenomena. In either case, thorough testing of existing procedures against future data will be required before this kind of system may become operational. coordinated by the Intergovernmental Oceanographic Commission of UNESCO. The Italian contribution to this working group initiated this review.
Funding Open access funding provided by Istituto Nazionale di Geofisica e Vulcanologia within the CRUI-CARE Agreement. This work benefited of the agreement between Istituto Nazionale di Geofisica e Vulcanologia and the Italian Presidenza del Consiglio dei Ministri, Dipartimento della Protezione Civile (DPC). This paper does not necessarily represent DPC official opinion and policies.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.