Search

Tsunami urgent computing procedures quantify the potential hazard due to an earthquake right after its occurrence, that is within a few hours. The hazard is quantified by simulating the propagation of the tsunami waves in the sea, accounting for the uncertainty due to the scarce knowledge of the source parameters and wave modelling uncertainty.In the context of the European project eflows4HPC, a workflow is currently in development for tsunamis hazard urgent computing, which consists of the following steps: 1) Retrieval of information about the tsunamigenic seismic event (magnitude, hypocentre and their uncertainties); 2) Definition of an ensemble of seismic sources; 3) Simulation of seismic/tsunamigenic waves propagation for each scenario in the ensemble; 4) Results aggregation to produce an estimate of seismic and tsunami hazard, which also incorporates a basic treatment of modelling uncertainty. The ensembles cover the uncertainty on source characteristics and may consequently be very large (generally 10,000 to 100,000 of scenarios; Selva et al., Nat. Comm.), requiring very high computational resources for the urgent computing context. It is thus necessary to reduce the size of these ensembles to limit the number of simulations and to converge faster towards stable results of hazard calculation.We developed and tested several sampling procedures aiming to reduce the number of scenarios in the ensemble and, at the same time, to integrate the new incoming information as they become available (e.g. solutions for focal mechanisms, seismic or tsunami records). When applied to several past earthquakes and tsunamis (e.g., the 2003 Boumerdes and the 2017 Kos-Bodrum earthquakes), our novel sampling strategies yielded a reduction of 1 or 2 order of magnitudes of the ensemble size, allowing a drastic reduction of the computational effort. Also, the update of the ensemble based on the incoming of new data, which strongly reduce the uncertainty, yields to an update of the probabilistic forecasts without compromising its accuracy. This may result very important for mitigating the risk far from the seismic source, as well as improving the risk management by better informing decision making in a frame of urgency.

The 2023 February 6 Mw 7.8 earthquake was the first one of a doublet which shook Türkiye and Syria causing, as per the estimates at the time of writing of this abstract, more than 45,000 casualties.The current standard operating procedures of the NEAMTWS (Tsunami Warning System in the North-Eastern Atlantic, the Mediterranean and connected seas, coordinated by UNESCO/IOC) for the initial tsunami warning message following an earthquake are based on a Decision Matrix (DM), whose input parameters are hypocentre and magnitude of the earthquake. Since the epicentre of this earthquake was located at a depth between 15-35 km at almost 100 km from the coast, both KOERI (Türkiye) and INGV (Italy) Tsunami Service Providers (TSPs) of the NEAMTWS issued a Tsunami Watch message (i.e., runup expected to exceed 1 m) for the whole Mediterranean Sea. NOA (Greece) did not issue any alert, because its initial location was more than 100 km from the coast.In response to the tsunami warning, trains were stopped in different locations in Southern Italy for several hours, and evacuation of some coastal areas was enforced. However, only a relatively small tsunami was recorded by Turkish close-by tide-gauges in the Eastern Mediterranean, with a maximum recorded amplitude of less than 50 cm. Based on these measurements and on others showing little to no tsunami at increasing distances, the alert was then ended after 5 and 9 hours by INGV and KOERI, respectively, based on the available tide-gauge recordings and interaction with Civil Protection Officers.This event has highlighted that NEAMTWS is an asset for the coastal communities. It can provide rapid alerts, which can save lives if the last-mile of the procedures is in place and the communities are “Tsunami Ready”, that is aware and prepared to respond with evacuations and other appropriate countermeasures. On the other hand, while it is reasonable – and dutiful based on current standard operation procedures – to issue a basin-wide, or at least a local alert, for an inland earthquake of unknown mechanism and of such a large magnitude, it is perhaps possible to improve the DM, which is totally heuristic and characterized by hard-thresholds, with consideration of numerical tsunami simulations and quantitative uncertainty treatment with more continuous variations. Moreover, there is no procedure currently in place to differentiate among locations where the expected time of arrival differs by many hours across the Mediterranean basin, nor a sufficient instrumental coverage that could make cancellation/ending faster due to a more solid observational basis.We will discuss some of the scientific and operational aspects with the aim of identifying which lessons can be learned to improve the NEAMTWS efficiency. We will also compare the DM-based alerts with those that would be produced with the recently introduced Probabilistic Tsunami Forecasting (PTF, Selva et al., 2021, Nature Communications), presently in pre-operational testing at INGV.

The destructive tsunami on 22 December 2018 due to the flank collapse of the Anak Krakatau volcano was a bitter reminder of large tsunami risks and of the shortcomings of the existing tsunami warning systems for atypical (non-seismic) sources. In the Mediterranean, several tsunamis were generated by landslides associated with volcanic systems in the past.The volcanic unrest experienced in 2011-2012 on the Santorini volcanic island in the Southern Aegean Sea pointed out the need to identify and quantify tsunami hazard and risk due to possible flank instability which may be triggered as a result of volcanic unrest or nearby seismotectonic activities. Inspired from this need, in this study we examined three possible landslide scenarios in Santorini Island with tsunamigenic potential. The results show that the scenarios considered in our study are able to generate significant local tsunamis impacting Santorini and the nearby islands, as well as producing significant impact along the coasts of the Southern Aegean Sea. While maximum tsunami amplitudes/arrival time ranges are 1.2m/30-90min for locations in the Greek-Turkish coasts in the far field, they are in the order of ≈60m/1-2min for some locations at the Santorini Island. The extreme tsunami amplitudes and short arrival times for locations inside the Santorini Island is a major challenge in terms of tsunami hazard warning and mitigation. As an effort to address this challenge, a discussion on the requirements for local tsunami warning system addressing atypical sources in the context of multi-hazard disaster risk reduction is also provided., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

The creation of hazard maps relative to volcanic phenomena requires taking into account the intrinsic complexity and variability of eruptions. Here we present an example of how HPC can allow producing a high resolution multi-source probabilistic hazard assessment due to tephra fallout over a domain covering Southern Italy. The three active volcanoes in the Neapolitan area, Somma-Vesuvius, Campi Flegrei and Ischia, were considered as volcanic sources. For each one, we explored three explosive size classes (Small, Medium and Large) for Somma Vesuvius and Campi Flegrei, and one explosive size class (Large) for Ischia. For each size class, we performed 1500 numerical simulations of ash dispersion (total of 10500) using the Fall3D (V8.0) model over a computational domain covering Southern Italy with a 0.03° ⨉ 0.03° (~3 km ⨉ 3 km) resolution. Within each size class, the eruptive parameters have been randomly sampled from well-suited probability distributions and with different meteorological conditions, obtained by randomly sampling a day between 1990 and 2020 and retrieving the relative data from the ECMWF ERA5 database. This allowed exploring the intra-class variability and to quantify aleatoric uncertainty. The results of these simulations have been post-processed with a statistical approach by assigning a weight to each eruption (based on its eruption magnitude) and the annual eruption rate of each size class. For the case of Campi Flegrei, the variability in the eruptive vent position has also been explored by constructing a grid of possible vent locations with different spatial probability. By merging the results obtained for each source and size class we produced a portfolio of hazard maps showing the expected mean annual frequency of overcoming selected thresholds in ground tephra load. A disaggregation analysis has also been performed in order to understand which particular source and/or size class had the greater impact on a particular area. The completion of this work, considering both numerical simulations and the statistical elaboration of the results has required a total of more than 5000 core hours and the processing of more than 2TB of data, an effort that wouldn’t have been possible without the access to high level HPC resources.

Tsunami urgent computing procedures quantify the potential hazard due to a seismically-induced tsunami right after an earthquake, that is from minutes to a few hours. The hazard is quantified by simulating the tsunami from source to shore, taking into account the uncertainty in the source parameters and the uncertainty associated with the wave generation, propagation, and inundation. In the European eFlows4HPC project, an HPC workflow for urgent computing of tsunami hazard assessment is currently being developed, consisting of the following steps: 1) retrieval of parameters for the tsunamigenic earthquake (magnitude, hypocentre and associated uncertainties), 2) definition of a seismic source ensemble, 3) simulation of the tsunami generated by each scenario in the ensemble, 4) aggregation of the results to produce an estimate of tsunami hazard, which also incorporates a basic treatment of uncertainty modelling and 5) update of the ensemble based on incoming data. Initially implemented on the Power-9 machine at BSC, the workflow has been fully embedded into a PyCOMPSs framework that enables parallel task execution and integrates full tsunami simulations for the first time. The tsunami numerical model (Tsunami-HySEA) computes the tsunami from the source to coastal impact using nested grids with resolution from kilometres to meters. To limit the number of simulations and converge faster towards stable hazard estimates, new methods for defining the seismic source ensembles have been developed. When applied to several past earthquakes and tsunamis (e.g., the 2003 Boumerdes and the 2017 Kos-Bodrum earthquakes), our new sampling strategy yielded a reduction of 1 or 2 orders of magnitude for ensemble size, allowing a drastic reduction in the computational effort. This reduction may be exploited to improve tsunami simulation accuracy, increasing the computational effort available for each simulation for the same overall cost. The workflow also allows the integration of new incoming data (focal mechanism, seismic or tsunami records) for an “on the fly” update of the PTF based on this new information. The improvement of the workflow through a well-defined ensemble of scenarios, realistic simulations and integration of incoming data, strongly reduces the uncertainty and yields to an update of the probabilistic forecasts without compromising theiraccuracy. This can be crucial in mitigating the risk far from the seismic source, and in improving risk management by better informing decision-making in an emergency framework.

Nowadays, tephra fallout hazard is based on coupling the physical modeling of the tephra dispersion processes with a probabilistic analysis that takes into account the natural variability of the volcanic phenomena in terms of eruption probability, eruption sizes, vent position and meteorological conditions. In this framework, we present a prototypal methodology to carry out a multi-volcano long-term tephra fallout hazard assessment in Southern Italy from the active Neapolitan volcanoes: Somma-Vesuvius, Campi Flegrei, and Ischia. FALL3D model (v.8.0) has been used to run thousands of numerical simulations (1,500 per eruption size class), considering the ECMWF ERA5 meteorological dataset over the last 30 years. The output in terms of tephra ground load has been processed within a new workflow for large-scale, high-resolution volcanic hazard assessment, in order to quantify the mean annual frequency with which the tephra load at the ground exceeds given critical thresholds at a target site within a 50-years exposure time, and the relative epistemic uncertainty. This work provides, for the first time, a multi-volcano probabilistic hazard analysis for tephra fallout, fully comparable with those used for seismic phenomena and other natural disasters in which multiple sources are integrated together, and it accounts for potential changes in regimes of each single considered volcano. This allows us to discuss also how the full information can be traced back to provide specific information about the prevalence of different volcanoes and eruptive style in the different target areas, based on hazard disaggregation. The methodology is applicable to any other volcanic areas or over different exposure times.

The main purpose of this article is to emphasize the importance of clarifying the probabilistic framework adopted for volcanic hazard and eruption forecasting. Eruption forecasting and volcanic hazard analysis seeks to quantify the deep uncertainties that pervade the modeling of pre-, sin- and post-eruptive processes. These uncertainties can be differentiated into three fundamental types: (1) the natural variability of volcanic systems, usually represented as stochastic processes with parameterized distributions (aleatory variability); (2) the uncertainty in our knowledge of how volcanic systems operate and evolve, often represented as subjective probabilities based on expert opinion (epistemic uncertainty); and (3) the possibility that our forecasts are wrong owing to behaviors of volcanic processes about which we are completely ignorant and, hence, cannot quantify in terms of probabilities (ontological error). Here we put forward a probabilistic framework for hazard analysis recently proposed by Marzocchi & Jordan (2014), which unifies the treatment of all three types of uncertainty. Within this framework, an eruption forecasting or a volcanic hazard model is said to be complete only if it (a) fully characterizes the epistemic uncertainties in the model's representation of aleatory variability and (b) can be unconditionally tested (in principle) against observations to identify ontological errors. Unconditional testability, which is the key to model validation, hinges on an experimental concept that characterizes hazard events in terms of exchangeable data sequences with well-defined frequencies. We illustrate the application of this unified probabilistic framework by describing experimental concepts for the forecasting of tephra fall from Campi Flegrei. Eventually, this example may serve as a guide for the application of the same probabilistic framework to other natural hazards.

Volcanic activity typically switches between high-activity states with many eruptions and low-activity states with few or no eruptions. We present a simple two-regime physics-informed statistical model that allows interpreting temporal modulations in eruptive activity. The model enhances comprehension and comparison of different volcanic systems and enables homogeneous integration into multivolcano hazard assessments that account for potential changes in volcanic regimes. The model satisfactorily fits the eruptive history of the three active volcanoes in the Neapolitan area, Italy (Mt. Vesuvius, Campi Flegrei, and Ischia) which encompass a wide range of volcanic behaviors. We find that these volcanoes have appreciably different processes for triggering and ending high-activity periods connected to different dominant volcanic processes controlling their eruptive activity, with different characteristic times and activity rates (expressed as number of eruptions per time interval). Presently, all three volcanoes are judged to be in a low-activity state, with decreasing probability of eruptions for Mt. Vesuvius, Ischia, and Campi Flegrei, respectively.

Special issue EPOS a Research Infrastructure in solid Earth: open science and innovation .-- 21 pages, 8 figures, Tsunamis constitute a significant hazard for European coastal populations, and the impact of tsunami events worldwide can extend well beyond the coastal regions directly affected. Understanding the complex mechanisms of tsunami generation, propagation, and inundation, as well as managing the tsunami risk, requires multidisciplinary research and infrastructures that cross national boundaries. Recent decades have seen both great advances in tsunami science and consolidation of the European tsunami research community. A recurring theme has been the need for a sustainable platform for coordinated tsunami community activities and a hub for tsunami services. Following about three years of preparation, in July 2021, the European tsunami community attained the status of Candidate Thematic Core Service (cTCS) within the European Plate Observing System (EPOS) Research Infrastructure. Within a transition period of three years, the Tsunami candidate TCS is anticipated to develop into a fully operational EPOS TCS. We here outline the path taken to reach this point, and the envisaged form of the future EPOS TCS Tsunami. Our cTCS is planned to be organised within four thematic pillars: (1) Support to Tsunami Service Providers, (2) Tsunami Data, (3) Numerical Models, and (4) Hazard and Risk Products. We outline how identified needs in tsunami science and tsunami risk mitigation will be addressed within this structure and how participation within EPOS will become an integration point for community development, With the institutional support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S)

Probabilistic volcanic hazard assessment is a standard methodology based on running a deterministic hazard quantification tool multiple times to explore the full range of uncertainty in the input parameters and boundary conditions, in order to probabilistically quantify the variability of outputs accounting for such uncertainties. Nowadays, different volcanic hazards are quantified by means of this approach. Among these, volcanic gas emission is particularly relevant given the threat posed to human health if concentrations and exposure times exceed certain thresholds. There are different types of gas emissions but two main scenarios can be recognized: hot buoyant gas emissions from fumaroles and the ground and dense gas emissions feeding density currents that can occur, e.g., in limnic eruptions. Simulation tools are available to model the evolution of critical gas concentrations over an area of interest. Moreover, in order to perform probabilistic hazard assessments of volcanic gases, simulations should account for the natural variability associated to aspects such as seasonal and daily wind conditions, localized or diffuse source locations, and gas fluxes. Here we present VIGIL (automatized probabilistic VolcanIc Gas dIspersion modeLling), a new Python tool designed for managing the entire simulation workflow involved in single and probabilistic applications of gas dispersion modelling. VIGIL is able to manage the whole process from meteorological data processing, needed to run gas dispersion in both the dilute and dense gas flow scenarios, to the post processing of models' outputs. Two application examples are presented to show some of the modelling capabilities offered by VIGIL., Fabio Dioguardi (FD) has been supported from UK National Capability funding (BGS Innovation Flexible Fund). This funding supported the code preparation, testing and publishing. Fabio Dioguardi (FD), Silvia Massaro (SM), Laura Sandri (LS), Jacopo Selva (JS), Giovanni Macedonio (GM) and Antonio Costa (AC) have also been supported by the European Union’s Horizon 2020 project EUROVOLC (grant agreement no 731070). This funding supported the initial idea and the interaction among the co‑authors. We thank John A. Stevenson (British Geological Survey) for his technical support in the development of VIGIL and its repository. This work is published with permission of the Executive Director of British Geological Survey (UKRI)

The b-value of the Gutenberg-Richter law is one of the most widely studied parameters regarding the distribution of earthquakes’ magnitude. The estimation of such a parameter and its uncertainty is critical, and it may become complex in the case of seismic catalogs with a non-uniform magnitude of completeness, or when different shapes of the Gutenberg-Richter relation (e.g. the tapered one) are adopted. Here we review recent results on the b-value estimation, also in the case of catalogs with multiple completeness levels, including the application of the weighted likelihood methodology, a method particularly suitable for spatial b-value mapping and for studying its temporal variations. These new techniques are applied to both global and regional seismic catalogs, in order to unveil the peculiarities of the b-value.

The effects of volcanic hazards can be quantified by applying new methods to provide support for rational decision-making. Mt Etna is one the most active volcanoes in the world, producing both effusive and explosive eruptions together with a very intense seismic activity, which significantly affect the territory and human society. We present the preliminary results obtained in the framework of the PANACEA project (INGV’s project “Pianeta Dinamico”, funded by the Italian Ministero dell’Università e la Ricerca, MUR) regarding the multi-hazard assessment around Mt Etna. These include: (i) the production of an updated spatio-temporal probability map of vent opening at Etna, using a procedure exploiting different Kernel functions (e.g. the exponential, Cauchy, and Gaussian functions), and testing volcanic deformation patterns to explore possible dynamic, structural conditioning on the vent opening process; (ii) the identification of a set of cascading effects scenarios that account for volcanic phenomena (i.e., volcanic unrest, seismicity, volcanic explosions, volcanic effusive events, lava flows, tephra/ballistic fall, and PDC), as well as other external hazards potentially linked in such chains (e.g., flooding, forest fires, etc.); and (iii) the identification of scenarios involving systemic impacts (e.g., impacts on the functionality or connectivity of networks).

The coasts of the Mediterranean Sea are densely populated and exposed to tsunami inundations as reported by historical evidence. Measures to mitigate the tsunami risk in this region are based on Probabilistic Tsunami Hazard Assessments (PTHA) computed considering present coastal morphologies. However, mean sea level projections for the 21st century indicated a general sea level rise which can be substantially modified if uplift or subsidence may occur locally due to other geological factors. In order to reduce the potential impact of tsunamis all factors (climatic or not) should be included in the tsunami hazard analysis. In this study we focus on the Eastern Sicily and we examine how the PTHA can significantly change when the general trend of sea level rise, based on AR-5 and AR-6 IPCC climate scenarios and rates of Vertical Land Movements, are included in the region. Moreover, we take into account associated epistemic uncertainties related to the future sea level rise under different conditions of low- and high-emission representative concentrations.

We present a refined methodological procedure for computationally efficient local SPTHA based on regional SPTHA. The adopted procedure extracts from the regional SPTHA the most impacting tsunami sources at the investigated site, and reconstructs hazard curves on high-resolution topobathymetric models based on a reduced set of inundation simulations. This procedure enhances the original workflow for local SPTHA quantification described by Volpe et al. (2019), applying some significant upgrades to simplify its application and improve the accuracy of the results. In particular, the description of local sources has been refined through a more detailed discretization of the natural variability (aleatory uncertainty), eventually reducing the epistemic uncertainty. Then, a more efficient filtering procedure, based on the strategy proposed by Williamson et al. (2020), is adopted to select a subset of scenarios to be modelled at high resolution, eventually reducing the epistemic uncertainty introduced by this selection. This allows to perform only coarse-grid simulations after the regional source filtering and local source refinement, and then combine coarse-grid results with fine-grid topography. Overall, the resulting method simplifies the original one, improving accuracy and decreasing uncertainty. The newly developed procedure is applied to an illustrative case study for the harbour of Ravenna (Northern Adriatic Sea, Italy).

Supplementary data: The following is the Supplementary data to this article: Acrobat PDF file (2MB) available at: https://ars.els-cdn.com/content/image/1-s2.0-S2212420921007329-mmc1.pdf Copyright © 2022 The Authors. Very large tsunamis are associated with low probabilities of occurrence. In many parts of the world, these events have usually occurred in a distant time in the past. As a result, there is low risk perception and a lack of collective memories, making tsunami risk communication both challenging and complex. Furthermore, immense challenges lie ahead as population and risk exposure continue to increase in coastal areas. Through the last decades, tsunamis have caught coastal populations off-guard, providing evidence of lack of preparedness. Recent tsunamis, such as the Indian Ocean Tsunami in 2004, 2011 Tohoku and 2018 Palu, have shaped the way tsunami risk is perceived and acted upon. Based on lessons learned from a selection of past tsunami events, this paper aims to review the existing body of knowledge and the current challenges in tsunami risk communication, and to identify the gaps in the tsunami risk management methodologies. The important lessons provided by the past events call for strengthening community resilience and improvement in risk-informed actions and policy measures. This paper shows that research efforts related to tsunami risk communication remain fragmented. The analysis of tsunami risk together with a thorough understanding of risk communication gaps and challenges is indispensable towards developing and deploying comprehensive disaster risk reduction measures. Moving from a broad and interdisciplinary perspective, the paper suggests that probabilistic hazard and risk assessments could potentially contribute towards better science communication and improved planning and implementation of risk mitigation measures. COST (European Cooperation in Science and Technology); Royal Society, UK (grant number CHL\R1\180173); Severo Ochoa Centers of Excellence Program (CEX 2018-000797-S) funded by MCIN/ AEI /10.13039/501100011033; Lloyd's Tercentenary Research Foundation, the Lighthill Risk Network, and the Lloyd's Register Foundation-Data Centric Engineering Programme of the Alan Turing Institute.

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.

Recent natural disasters that seriously affected critical infrastructure (CI) with significant socio-economic losses and impact revealed the need for the development of reliable methodologies for vulnerability and risk assessment. In this paper, a risk-based multi-level stress test method that has been recently proposed, aimed at enhancing procedures for evaluation of the risk of critical non-nuclear infrastructure systems against natural hazards, is specified and applied to six key representative CIs in Europe, exposed to variant hazards. The following CIs are considered: an oil refinery and petrochemical plant in Milazzo, Italy, a conceptual alpine earth-fill dam in Switzerland, the Baku–Tbilisi–Ceyhan pipeline in Turkey, part of the Gasunie national gas storage and distribution network in the Netherlands, the port infrastructure of Thessaloniki, Greece, and an industrial district in the region of Tuscany, Italy. The six case studies are presented following the workflow of the stress test framework comprised of four phases: pre-assessment phase, assessment phase, decision phase and report phase. First, the goals, the method, the time frame and the appropriate stress test level to apply are defined. Then, the stress test is performed at component and system levels and the outcomes are checked and compared to risk acceptance criteria. A stress test grade is assigned, and the global outcome is determined by employing a grading system. Finally, critical components and events and risk mitigation strategies are formulated and reported to stakeholders and authorities. © 2019, Springer Nature B.V.

Site-specific seismic probabilistic tsunami hazard analysis (SPTHA) is a computationally demanding task, as it requires, in principle, a huge number of high-resolution numerical simulations for producing probabilistic inundation maps. We implemented an efficient and robust methodology using a filtering procedure to reduce the number of numerical simulations needed while still allowing for a full treatment of aleatory and epistemic uncertainty. Moreover, to avoid biases in tsunami hazard assessment, we developed a strategy to identify and separately treat tsunamis generated by near-field earthquakes. Indeed, the coseismic deformation produced by local earthquakes necessarily affects tsunami intensity, depending on the scenario size, mechanism and position, as coastal uplift or subsidence tends to diminish or increase the tsunami hazard, respectively. Therefore, we proposed two parallel filtering schemes in the far- and the near-field, based on the similarity of offshore tsunamis and hazard curves and on the similarity of the coseismic fields, respectively. This becomes mandatory as offshore tsunami amplitudes can not represent a proxy for the coastal inundation in the case of near-field sources. We applied the method to an illustrative use case at the Milazzo oil refinery (Sicily, Italy). We demonstrate that a blind filtering procedure can not properly account for local sources and would lead to a nonrepresentative selection of important scenarios. For the specific source–target configuration, this results in an overestimation of the tsunami hazard, which turns out to be correlated to dominant coastal uplift. Different settings could produce either the opposite or a mixed behavior along the coastline. However, we show that the effects of the coseismic deformation due to local sources can not be neglected and a suitable correction has to be employed when assessing local-scale SPTHA, irrespective of the specific signs of coastal displacement.

Regional and global tsunami hazard analysis requires simplified and efficient methods for estimating the tsunami inundation height and its related uncertainty. One such approach is the amplification factor (AF) method. Amplification factors describe the relation between offshore wave height and the maximum inundation height, as predicted by linearized plane wave models employed for incident waves with different wave characteristics. In this study, a new amplification factor method is developed that takes into account the offshore bathymetry proximal to the coastal site. The present AFs cover the North-Eastern Atlantic and Mediterranean (NEAM) region. The model is the first general approximate model that quantifies inundation height uncertainty. Uncertainty quantification is carried out by analyzing the inundation height variability in more than 500 high-resolution inundation simulations at six different coastal sites. The inundation simulations are undertaken with different earthquake sources in order to produce different wave period and polarity. We show that the probability density of the maximum inundation height can be modeled with a log-normal distribution, whose median is quite well predicted by the AF. It is further demonstrated that the associated maximum inundation height uncertainties are significant and must be accounted for in tsunami hazard analysis. The application to the recently developed TSUMAPS-NEAM probabilistic tsunami hazard analysis (PTHA) is presented as a use case.

Inundation maps are a fundamental tool for coastal risk management and in particular for designing evacuation maps and evacuation planning. These in turn are a necessary component of the tsunami warning systems’ last-mile. In Italy inundation maps are informed by a probabilistic tsunami hazard model. Based on a given level of acceptable risk, Italian authorities in charge for this task recommended to consider, as design hazard intensity, the average return period of 2500 years and the 84th percentile of the hazard model uncertainty. An available, regional-scale tsunami hazard model was used that covers the entire Italian coastline. Safety factors based on analysis of run-up variability and an empirical coastal dissipation law on a digital terrain model (DTM) were applied to convert the regional hazard into the design run-up and the corresponding evacuation maps with a GIS-based approach. Since the regional hazard cannot fully capture the local-scale variability, this simplified and conservative approach is considered a viable and feasible practice to inform local coastal risk management in the absence of high-resolution hazard models. The present work is a first attempt to quantify the uncertainty stemming from such procedure. We compare the GIS-based inundation maps informed by a regional model with those obtained from a local high-resolution hazard model. Two locations on the coast of eastern Sicily were considered, and the local hazard was addressed with the same seismic model as the regional one, but using a higher-resolution DTM and massive numerical inundation calculations with the GPU-based Tsunami-HySEA nonlinear shallow water code. This study shows that the GIS-based inundation maps used for planning deal conservatively with potential hazard underestimation at the local scale, stemming from typically unmodeled uncertainties in the numerical source and tsunami evolution models. The GIS-based maps used for planning fall within the estimated “error-bar” due to such uncertainties. The analysis also demonstrates the need to develop local assessments to serve very specific risk mitigation actions to reduce the uncertainty. More in general, the presented case-studies highlight the importance to explore ways of dealing with uncertainty hidden within the high-resolution numerical inundation models, e.g., related to the crude parameterization of the bottom friction, or the inaccuracy of the DTM.

Advances in GPU-based High-Performance Computing (HPC) facilities, combined with improvements in GPU-optimized shallow water models for tsunami inundation, allow us to perform large numbers of numerical simulations of earthquake-generated tsunamis on high-resolution numerical grids. Large numbers of simulations are necessary to investigate the multi-dimensional parameter space that defines the tsunami hazard, including situations where the tsunami is generated outside major tectonic structures, where fault geometry is uncertain and can take widely different orientations. With over 1500 numerical simulations, we perform suites of systematic parameter searches to investigate the sensitivity of inundation at the towns of Catania and Siracusa on Sicily to changes both in the earthquake source parameters and in the specification of the Manning friction coefficient. The inundation is modelled using the GPU-based Tsunami-HySEA code on a system of nested topo-bathymetric grids with a finest spatial resolution of 10 meters. We consider tsunamigenesis by large earthquakes with uniform slip where the location, focal depth, fault dimensions and slip, together with the angles of strike, dip, and rake, are defined by the standard Okada parameters. We consider sources both close to the shore, in which significant co-seismic deformation occurs, and offshore, where co-seismic deformation is negligible. For the offshore earthquake sources, we see systematic and intuitive changes in the inundation with changes in strike, dip, rake, and depth. For the near-shore sources, the dependency is far more complicated and co-seismic deformation becomes significant in determining the inundation. The sensitivity studies provide clear guidelines as to the necessary resolution for source discretization for Probabilistic Tsunami Hazard Analysis, with a need for a far finer discretization of local sources than for more distant sources. For a small number of earthquake sources, we study systematically the inundation as a function of the Manning Friction Coefficient. The sensitivity of the inundation to this parameter varies greatly for different earthquake sources and topo-bathymetry at the coastline of interest. An understanding of all these dependencies is needed to better understand the consequences of tsunamigenic earthquake models with more complex geometries, and in quantifying the epistemic uncertainty in the tsunami hazard. This work is partially funded by the European Union’s Horizon 2020 Research and Innovation Program under grant agreement No 823844 (ChEESE Center of Excellence, www.cheese-coe.eu). Computational resources made available through Sigma2/UNINETT on Saga at NTNU, Trondheim, Norway (in project nn5008k) and through PRACE on Marconi-100 at CINECA, Rome, Italy (through PRACE grant Pra21_5386/TsuHazAP).