Probabilistic tsunami hazard analysis (PTHA) relies on computationally demanding numerical simulations of tsunami generation, propagation, and non-linear inundation on high-resolution topo-bathymetric models. Here we focus on tsunamis generated by co-seismic sea floor displacement, that is, on Seismic PTHA (SPTHA). A very large number of tsunami simulations are typically needed to incorporate in SPTHA the full expected variability of seismic sources (the aleatory uncertainty). We propose an approach for reducing their number. To this end, we (i) introduce a simplified event tree to achieve an effective and consistent exploration of the seismic source parameter space; (ii) use the computationally inexpensive linear approximation for tsunami propagation to construct a preliminary SPTHA that calculates the probability of maximum offshore tsunami wave height (HMax) at a given target site; (iii) apply a two-stage filtering procedure to these ‘linear’ SPTHA results, for selecting a reduced set of sources and (iv) calculate ‘non-linear’ probabilistic inundation maps at the target site, using only the selected sources. We find that the selection of the important sources needed for approximating probabilistic inundation maps can be obtained based on the offshore HMax values only. The filtering procedure is semi-automatic and can be easily repeated for any target sites. We describe and test the performances of our approach with a case study in the Mediterranean that considers potential subduction earthquakes on a section of the Hellenic Arc, three target sites on the coast of eastern Sicily and one site on the coast of southern Crete. The comparison between the filtered SPTHA results and those obtained for the full set of sources indicates that our approach allows for a 75–80 per cent reduction of the number of the numerical simulations needed, while preserving the accuracy of probabilistic inundation maps to a reasonable degree.

Tsunami warning centres face the challenging task of rapidly forecasting tsunami threat immediately after an earthquake, when there is high uncertainty due to data deficiency. Here we introduce Probabilistic Tsunami Forecasting (PTF) for tsunami early warning. PTF explicitly treats data- and forecast-uncertainties, enabling alert level definitions according to any predefined level of conservatism, which is connected to the average balance of missed-vs-false-alarms. Impact forecasts and resulting recommendations become progressively less uncertain as new data become available. Here we report an implementation for near-source early warning and test it systematically by hindcasting the great 2010 M8.8 Maule (Chile) and the well-studied 2003 M6.8 Zemmouri-Boumerdes (Algeria) tsunamis, as well as all the Mediterranean earthquakes that triggered alert messages at the Italian Tsunami Warning Centre since its inception in 2015, demonstrating forecasting accuracy over a wide range of magnitudes and earthquake types.

We propose a procedure for uncertainty quantification in Probabilistic Tsunami Hazard Analysis (PTHA), with a special emphasis on the uncertainty related to statistical modelling of the earthquake source in Seismic PTHA (SPTHA), and on the separate treatment of subduction and crustal earthquakes (treated as background seismicity). An event tree approach and ensemble modelling are used in spite of more classical approaches, such as the hazard integral and the logic tree. This procedure consists of four steps: (1) exploration of aleatory uncertainty through an event tree, with alternative implementations for exploring epistemic uncertainty; (2) numerical computation of tsunami generation and propagation up to a given offshore isobath; (3) (optional) site-specific quantification of inundation; (4) simultaneous quantification of aleatory and epistemic uncertainty through ensemble modelling. The proposed procedure is general and independent of the kind of tsunami source considered; however, we implement step 1, the event tree, specifically for SPTHA, focusing on seismic source uncertainty. To exemplify the procedure, we develop a case study considering seismic sources in the Ionian Sea (central-eastern Mediterranean Sea), using the coasts of Southern Italy as a target zone. The results show that an efficient and complete quantification of all the uncertainties is feasible even when treating a large number of potential sources and a large set of alternative model formulations. We also find that (i) treating separately subduction and background (crustal) earthquakes allows for optimal use of available information and for avoiding significant biases; (ii) both subduction interface and crustal faults contribute to the SPTHA, with different proportions that depend on source-target position and tsunami intensity; (iii) the proposed framework allows sensitivity and deaggregation analyses, demonstrating the applicability of the method for operational assessments.

The 2011 Tohoku earthquake (M_w = 9.1) highlighted previously unobserved features for megathrust events, such as the large slip in a relatively limited area and the shallow rupture propagation. We use a Finite Element Model (FEM), taking into account the 3D geometrical and structural complexities up to the trench zone, and perform a joint inversion of tsunami and geodetic data to retrieve the earthquake slip distribution. We obtain a close spatial correlation between the main deep slip patch and the local seismic velocity anomalies, and large shallow slip extending also to the North coherently with a seismically observed low-frequency radiation. These observations suggest that the friction controlled the rupture, initially confining the deeper rupture and then driving its propagation up to the trench, where it spreads laterally. These findings are relevant to earthquake and tsunami hazard assessment because they may help to detect regions likely prone to rupture along the megathrust, and to constrain the probability of high slip near the trench. Our estimate of ~40 m slip value around the JFAST (Japan Trench Fast Drilling Project) drilling zone contributes to constrain the dynamic shear stress and friction coefficient of the fault obtained by temperature measurements to ~0.68 MPa and ~0.10, respectively.