Integrating diverse expert opinions in hazard and risk projects is essential to managing subjective decisions and quantifying uncertainty to produce stable and trustworthy results. A structured procedure is necessary to organize the gathering of experts' opinions while ensuring transparency, accountability, and independence in judgements. We propose a novel Multiple-Expert management Protocol (MEP) to address this challenge, providing procedural guidelines for conducting single to multi-hazard risk analyses. MEP establishes a workflow to manage subjectivity rooted in (i) moderated and staged group interactions, (ii) trackable blind advice through written elicitations with mathematical aggregation, (iii) participatory independent review, (iv) close cooperation between scientific and managerial coordination, and (v) proper and comprehensive documentation. Originally developed for stress testing critical infrastructure, MEP is designed as a single, flexible, technology-neutral procedural workflow applicable to various sectors. Moreover, its scalability allows it to adapt from high to low-budget projects and from complex probabilistic multi-hazard risk assessments to standard single-hazard analyses, with different experts' degree and type of involvement depending on available funding and emerging controversies. We present two compelling case studies to showcase MEP's practical applicability: a multi-hazard risk analysis for a port infrastructure and a single-hazard regional tsunami hazard assessment.

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

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.

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k −2 stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model 150,000 earthquakes with M_w from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

Despite having been rather seismically quiescent for the last decades, the Makran subduction zone is capable of hosting destructive earthquakes and tsunami. In particular, the well-known thrust event in 1945 (Balochistan earthquake) led to about 4000 casualties. Nowadays, the coastal regions are more densely populated and vulnerable to similar events. Furthermore, some recent publications discuss rare but significantly larger events at the Makran subduction zone as possible scenarios. We analyze the instrumental and historical seismicity at the subduction plate interface and generate various synthetic earthquake catalogs spanning 300 000 years with varying magnitude-frequency relations. For every event in the catalogs we compute estimated tsunami heights and present the resulting tsunami hazard along the coasts of Pakistan, Iran and Oman in the form of probabilistic tsunami hazard curves. We show how the hazard results depend on variation of the Gutenberg–Richter parameters and especially maximum magnitude assumption.

We present a database of pre-calculated tsunami waveforms for the entire Mediterranean Sea, obtained by numerical propagation of uniformly spaced Gaussian-shaped elementary sources for the sea level elevation. Based on any initial sea surface displacement, the database allows the fast calculation of full waveforms at the 50 m isobath offshore of coastal sites of interest by linear superposition. A computationally inexpensive procedure is set to estimate the coefficients for the linear superposition based on the potential energy of the initial elevation field. The elementary sources size and spacing is fine enough to satisfactorily reproduce the effects of M> = 6.0 earthquakes. Tsunami propagation is modelled by using the Tsunami-HySEA code, a GPU finite volume solver for the non-linear shallow water equations. Like other existing methods based on the initial sea level elevation, the database is independent on the faulting geometry and mechanism, which makes it applicable in any tectonic environment. We model a large set of synthetic tsunami test scenarios, selected to explore the uncertainty introduced when approximating tsunami waveforms and their maxima by fast and simplified linear combination. This is the first time to our knowledge that the uncertainty associated to such a procedure is systematically analysed and that relatively small earthquakes are considered, which may be relevant in the near-field of the source in a complex tectonic setting. We find that non-linearity of tsunami evolution affects the reconstruction of the waveforms and of their maxima by introducing an almost unbiased (centred at zero) error distribution of relatively modest extent. The uncertainty introduced by our approximation can be in principle propagated to forecast results. The resulting product then is suitable for different applications such as probabilistic tsunami hazard analysis, tsunami source inversions and tsunami warning systems.

Real‐time GPS is nowadays considered as a valuable component of next generation near‐field tsunami early warning systems able to provide fast and reliable source parameters. Looking for optimal methodologies and assessing corresponding uncertainties becomes an important task. We take the opportunity and consider the 2014 Pisagua event as a case study to explore tsunami forecast uncertainty related to the GPS‐based source inversion. We intentionally neglect all other sources of uncertainty (observation set, signal processing, wave simulation, etc.) and exclusively assess the effect of inversion technique. In particular, we compare three end‐member methods: (1) point‐source fastCMT (centroid moment tensor), (2) distributed slip along predefined plate interface, and (3) unconstrained inversion into a single uniform slip finite fault. The three methods provide significantly different far‐field tsunami forecast but show surprisingly similar tsunami predictions in the near field.

A long foreshock series unlocked the South American plate boundary until eventually initiating the M 8.1 Iquique, Chile, earthquake. Two groups publishing in this issue of Nature analyse the seismic context of the Iquique earthquake that occurred off the coast of northern Chile on 1 April 2014 in a seismic zone that had been quiescent since a significant event in 1877. Gavin Hayes et al. identify areas of remaining or elevated earthquake hazard along the megathrust fault in the region, and conclude that the 2014 Iquique event was not the earthquake that had been anticipated. Given that significant sections of the northern Chile subduction zone have not ruptured in almost 150 years, they suggest that it is likely that future megathrust earthquakes will occur south and potentially north of the 2014 Iquique sequence. Bernd Schurr et al. show that the April 2014 earthquake broke a central fraction of the 'northern Chile seismic gap', the last major segment of the South American plate boundary that had yet to rupture in the past century. From July 2013 up to the April earthquake they identify three seismic clusters along this part of the plate boundary, each lasting a few weeks, with earthquakes of increasing peak magnitudes. They conclude that these seismic clusters and their slip transients reflect a gradual weakening of the central part of the seismic gap that was instrumental in initiating the final failure. On 1 April 2014, Northern Chile was struck by a magnitude 8.1 earthquake following a protracted series of foreshocks. The Integrated Plate Boundary Observatory Chile monitored the entire sequence of events, providing unprecedented resolution of the build-up to the main event and its rupture evolution. Here we show that the Iquique earthquake broke a central fraction of the so-called northern Chile seismic gap, the last major segment of the South American plate boundary that had not ruptured in the past century1,2. Since July 2013 three seismic clusters, each lasting a few weeks, hit this part of the plate boundary with earthquakes of increasing peak magnitudes. Starting with the second cluster, geodetic observations show surface displacements that can be associated with slip on the plate interface. These seismic clusters and their slip transients occupied a part of the plate interface that was transitional between a fully locked and a creeping portion. Leading up to this earthquake, the b value of the foreshocks gradually decreased during the years before the earthquake, reversing its trend a few days before the Iquique earthquake. The mainshock finally nucleated at the northern end of the foreshock area, which skirted a locked patch, and ruptured mainly downdip towards higher locking. Peak slip was attained immediately downdip of the foreshock region and at the margin of the locked patch. We conclude that gradual weakening of the central part of the seismic gap accentuated by the foreshock activity in a zone of intermediate seismic coupling was instrumental in causing final failure, distinguishing the Iquique earthquake from most great earthquakes. Finally, only one-third of the gap was broken and the remaining locked segments now pose a significant, increased seismic hazard with the potential to host an earthquake with a magnitude of >8.5.

Currently, inversion of geodetic data for earthquake fault ruptures is most based on a uniform half-space earth model because of its closed-form Green’s functions. However, the layered structure of the crust can significantly affect the inversion results. The other effect, which is often neglected, is related to the curved fault geometry. Especially, fault planes of most mega thrust earthquakes vary their dip angle with depth from a few to several tens of degrees. Also the strike directions of many large earthquakes are variable. For simplicity, such curved fault geometry is usually approximated to several connected rectangular segments, leading to an artificial loss of the slip resolution and data fit. In this presentation, we introduce a free FORTRAN code incorporating with the layered crust structure and curved fault geometry in a user-friendly way. The name SDM stands for Steepest Descent Method, an iterative algorithm used for the constrained least-squares optimization. The new code can be used for joint inversion of different datasets, which may include systematic offsets, as most geodetic data are obtained from relative measurements. These offsets are treated as unknowns to be determined simultaneously with the slip unknowns. In addition, a-priori and physical constraints are considered. The a-priori constraint includes the upper limit of the slip amplitude and the variation range of the slip direction (rake angle) defined by the user. The physical constraint is needed to obtain a smooth slip model, which is realized through a smoothing term to be minimized with the misfit to data. In difference to most previous inversion codes, the smoothing can be optionally applied to slip or stress-drop. The code works with an input file, a well-documented example of which is provided with the source code. Application examples are demonstrated.

Taking the 2011 Tohoku earthquake as an example, we demonstrate the ability of real-time GPS to provide qualified tsunami early warning within minutes. While in earlier studies we demonstrated the power of the so-called GPS shield concept based on synthetic data, we here present a complete processing chain starting from actual GPS raw data and fully simulate the situation as it would be in a warning center. The procedure includes processing of GPS observations with predicted high precision orbits, inversion for slip and computation of the tsunami propagation and coastal warning levels. We show that in case of the Tohoku earthquake, it would be feasible to provide accurate tsunami warning as soon as 3 min after the beginning of the earthquake.