Search

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.

Subduction zones are the most seismically active regions in the world and hosted many great tsunamigenic earthquakes in the past, often with destructive coastal consequences. Hence, an accurate estimate of the tsunami forecast is crucial in Tsunami Early Warning Systems (TEWS) framework. However, the inherent uncertainties associated with the tsunami source estimation in real-time make tsunami forecasting challenging. In this study, we consider the South American subduction zone, where in the last 15 years occurred, three M8+ tsunamigenic earthquakes; in particular, we focus on the 2014 Mw 8.1 Iquique event.Here, we evaluate the variability of the tsunami forecasting for the Chilean coast as resulting i) from the coseismic slip model obtained by geophysical data inversion and ii) from an expeditious method for the tsunami source estimation, based on an extension of the well-known spectral approach. In the former method, we estimate the slip distribution of the 2014 Iquique earthquake by jointly inverting tsunami (DARTs and tide-gauges) and GPS data; we adopt a 3D fault geometry and Green’s functions approach.On the other hand, a set of stochastic slip models in the latter is generated through a Phase Variation Method (PVM), where realizations are obtained from both the wavenumber and phase spectra of the source.In the analysis, we also evaluate how the different physics complexity included in the tsunami modelling (e.g. by including dispersion or not) can be mapped into the tsunami forecasting uncertainty. Finally, as an independent check, we compare the predicted deformation field from the slip models (inverted or by PVM) with the RADARSAT-2 InSAR data.

We present a source solution for the tsunami generated by the Mw 6.6 earthquake that occurred on 2 May 2020, about 80 km offshore south of Crete, in the Cretan Passage, on the shallow portion of the Hellenic Arc subduction zone (HASZ). The tide gauges recorded this local tsunami on the southern coast of Crete and Kasos island. We used Crete tsunami observations to constrain the geometry and orientation of the causative fault, the rupture mechanism, and the slip amount. We first modelled an ensemble of synthetic tsunami waveforms at the tide gauge locations, produced for a range of earthquake parameter values as constrained by some of the available moment tensor solutions. We allow for both a splay and a back-thrust fault, corresponding to the two nodal planes of the moment tensor solution. We then measured the misfit between the synthetic and the Ierapetra observed marigram for each source parameter set. Our results identify the shallow, steeply dipping back-thrust fault as the one producing the lowest misfit to the tsunami data. However, a rupture on a lower angle fault, possibly a splay fault, with a sinistral component due to the oblique convergence on this segment of the HASZ, cannot be completely ruled out. This earthquake reminds us that the uncertainty regarding potential earthquake mechanisms at a specific location remains quite significant. In this case, for example, it is not possible to anticipate if the next event will be one occurring on the subduction interface, on a splay fault, or on a back-thrust, which seems the most likely for the event under investigation. This circumstance bears important consequences because back-thrust and splay faults might enhance the tsunamigenic potential with respect to the subduction interface due to their steeper dip. Then, these results are relevant for tsunami forecasting in the framework of both the long-term hazard assessment and the early warning systems.

We present a source solution for the tsunami generated by the Mw 6.6 earthquake that occurred on May 2, 2020, about 807thinsp;km offshore south of Crete, in the Cretan Passage, on the shallow portion of the Hellenic Arc Subduction Zone (HASZ). The tide-gauges recorded this local tsunami on the southern coast of Crete island and Kasos island. We used these tsunami observations to constrain the geometry and orientation of the causative fault, the rupture mechanism and the slip amount. We first modelled an ensemble of synthetic tsunami waveforms at the tide-gauge locations, produced for a range of earthquake parameter values as constrained by some of the available moment tensor solutions. We allow for both a splay and a back-thrust fault, corresponding to the two nodal planes of the moment tensor solution. We then measured the misfit between the synthetic and the observed marigrams for each source parameter set. Our results identify the shallow steeply-dipping back-thrust fault as the one producing the lowest misfit to the tsunami data. However, a rupture on a lower angle fault, possibly a splay fault, with a sinistral component due to the oblique convergence on this segment of the HASZ, cannot be completely ruled out. This earthquake reminds us that the uncertainty regarding potential earthquake mechanisms at a specific location remains quite significant. In this case, for example, it is not possible to anticipate if the next event will be one occurring on the subduction interface, on a splay fault, or on a back-thrust which seems the most likely for the event under investigation. This circumstance bears important consequences because back-thrust and splay faults might enhance the tsunamigenic potential with respect to the subduction interface due to their steeper dip. Then, these results are relevant for tsunami forecasting both in the framework of the long-term hazard assessment and of the early warning systems.

On 28 September 2018 a Mw 7.5 strike-slip earthquake occurred on the Palu-Koro fault system in the Sulawesi Island. Immediately after the earthquake a powerful tsunami hit the Palu Bay causing large damages and numerous fatalities.Several works, inverting seismic or geodetic data, clearly estimated the slip distribution of this event, but the causative source of the tsunami is still not completely understood; indeed, the strike-slip mechanism of the seismic source alone might not be sufficient to explain the large runups observed (> 6 m) along the coast of the Palu Bay, and thus one or more additional non-seismic sources like a landslide could have contributed to generate the big tsunami. An insight of that can be found in an extraordinary collection of amateur videos, and on the only available tide gauge in the Bay, at Pantoloan, that showed evidence for a short period wave of at least 2-3 minutes, compatible with a landslide.In this study, we attempt to discriminate the contribution in the tsunami generation of both the seismic source and some supposed landslides distributed along the coast of the Bay.In particular, we attempt to estimate the causative source of the tsunami by means of a nonlinear joint inversion of geodetic (InSAR) and runup data. We use a fault geometry consistent with the Sentinel-2 optical analysis results and analytically compute the geodetic Green’s functions. The same fault model is used to compute the initial condition for the seismic tsunami Green’s functions, including the contribution of the horizontal deformation due to the gradient of the bathymetry (10 m spatial resolution); the landslide tsunami Green’s functions are computed the software BingClaw by placing several hypothetical sources in the Bay. In both the cases the tsunami propagation is modelled by numerically solving the nonlinear shallow water equations.In this work we also attempt to address the validity of Green’s functions approach (linearity) for earthquake and landslide sources as well as the wave amplitude offshore as predictor of nearby runup.

We present the realization of a fault-source data set designed to become the starting point in regional-scale tsunami hazard studies. Our approach focuses on the parametric fault characterization in terms of geometry, kinematics, and assessment of activity rates, and includes a systematic classification in six justification levels of epistemic uncertainty related with the existence and behaviour of fault sources. We set up a case study in the central Mediterranean Sea, an area at the intersection of the European, African, and Aegean plates, characterized by a complex and debated tectonic structure and where several tsunamis occurred in the past. Using tsunami scenarios of maximum wave height due to crustal earthquakes (Mw=7) and subduction earthquakes (Mw=7 and Mw=8), we illustrate first-order consequences of critical choices in addressing the seismogenic and tsunamigenic potentials of fault sources. Although tsunamis generated by Mw=8 earthquakes predictably affect the entire basin, the impact of tsunamis generated by Mw=7 earthquakes on either crustal or subduction fault sources can still be strong at many locales. Such scenarios show how the relative location/orientation of faults with respect to target coastlines coupled with bathymetric features suggest avoiding the preselection of fault sources without addressing their possible impact onto hazard analysis results.

On 6 February 2013 an Mw = 8.0 subduction earthquake occurred close to Santa Cruz Islands at the transition between the Solomon and the New Hebrides Trench. The ensuing tsunami caused significant inundation on the closest Nendo Island. The seismic source was studied with teleseismic broadband P-wave inversion optimized with tsunami forward modelling at DART buoys (Lay et al., 2013) and with inversion of teleseismic body and surface waves (Hayes et al., 2014a). The two studies also use different hypocentres and different planar fault models and found quite different slip models. In particular, Hayes et al. (2014a) argued for an aseismic slip patch SE from the hypocentre. We here develop a 3-D model of the fault surface from seismicity analysis and retrieve the tsunami source by inverting DART and tide-gauge data. Our tsunami source model features a main slip patch (peak value of ~ 11 m) SE of the hypocentre and reaching the trench. The rake direction is consistent with the progressively more oblique plate convergence towards the Solomon trench. The tsunami source partially overlaps the hypothesized aseismic slip area, which then might have slipped coseismically.

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 coastal sites by linear superposition. A computationally inexpensive procedure is set to estimate the coefficients for the linear superposition. 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 volumes solver for the shallow water equations. Like other existing methods based on the initial sea level elevation, this database is independent on the faulting geometry and mechanism, which makes it applicable in any tectonic environment. However, this is the first time to our knowledge that the uncertainty associated to such a procedure is systematically analysed, and that relatively small earthquakes, which may be relevant in the near field of the source in a complex tectonic setting, are considered. 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. The resulting product is suitable for different applications such as probabilistic tsunami hazard analysis, tsunami source inversions and tsunami warning systems.

In this paper we present and discuss the performance of the procedure for earthquake location and characterization implemented in the Italian Candidate Tsunami Service Provider at the Istituto Nazionale di Geofisica e Vulcanologia (INGV) in Rome. Following the ICG/NEAMTWS guidelines, the first tsunami warning messages are based only on seismic information, i.e., epicenter location, hypocenter depth, and magnitude, which are automatically computed by the software Early-est. Early-est is a package for rapid location and seismic/tsunamigenic characterization of earthquakes. The Early-est software package operates using offline-event or continuous-real-time seismic waveform data to perform trace processing and picking, and, at a regular report interval, phase association, event detection, hypocenter location, and event characterization. Early-est also provides mb, Mwp, and Mwpd magnitude estimations. mb magnitudes are preferred for events with Mwp ≲ 5.8, while Mwpd estimations are valid for events with Mwp ≳ 7.2. In this paper we present the earthquake parameters computed by Early-est between the beginning of March 2012 and the end of December 2014 on a global scale for events with magnitude M ≥ 5.5, and we also present the detection timeline. We compare the earthquake parameters automatically computed by Early-est with the same parameters listed in reference catalogs. Such reference catalogs are manually revised/verified by scientists. The goal of this work is to test the accuracy and reliability of the fully automatic locations provided by Early-est. In our analysis, the epicenter location, hypocenter depth and magnitude parameters do not differ significantly from the values in the reference catalogs. Both mb and Mwp magnitudes show differences to the reference catalogs. We thus derived correction functions in order to minimize the differences and correct biases between our values and the ones from the reference catalogs. Correction of the Mwp distance dependency is particularly relevant, since this magnitude refers to the larger and probably tsunamigenic earthquakes. Mwp values at stations with epicentral distance Δ ≲ 30° are significantly overestimated with respect to the CMT-global solutions, whereas Mwp values at stations with epicentral distance Δ ≳ 90° are slightly underestimated. After applying such distance correction the Mwp provided by Early-est differs from CMT-global catalog values of about δ Mwp ≈ 0.0 ∓ 0.2. Early-est continuously acquires time-series data and updates the earthquake source parameters. Our analysis shows that the epicenter coordinates and the magnitude values converge within less than 10 min (5 min in the Mediterranean region) toward the stable values. Our analysis shows that we can compute Mwp magnitudes that do not display short epicentral distance dependency overestimation, and we can provide robust and reliable earthquake source parameters to compile tsunami warning messages within less than 15 min after the event origin time.

On 6 February 2013 an Mw 8.0 subduction earthquake occurred close to Santa Cruz Islands at the transition between the Solomon and the New Hebrides Trench. The ensuing tsunami caused significant inundation on the closest Nendo Island. The seismic source was studied with teleseismic broadband P waves inversion optimized with tsunami forward modeling at DART buoys (Lay et al., 2013), and with inversion of teleseismic body and surface waves (Hayes et al., 2014). The two studies also use different hypocenters and different planar fault models, and found quite different slip models. In particular, Hayes et al. (2014) argued for an aseismic slip patch SE from the hypocenter. We here develop a 3-D model of the fault surface from seismicity analysis and retrieve the tsunami source by inverting DART and tide-gauge data. Our tsunami source model features a main slip patch (peak value of ~11 m) SE of the hypocentre, and reaching to the trench. The rake direction is consistent with the progressively more oblique plate convergence towards the Solomon trench. The tsunami source partially overlaps the hypothesized aseismic slip area, which then might have slipped coseismically.