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.

The Geo-INQUIRE (Geosphere INfrastructure for QUestions into Integrated REsearch) project, supported by the Horizon Europe Programme, is aimed at enhancing services to make data and high-level products accessible to the broad Geoscience scientific community. Geo-INQUIRE’s goal is to encourage curiosity-driven studies into understanding the geosphere dynamics at the interface between the solid Earth, the oceans and the atmosphere using long data streams, high-performance computing and cutting-edge facilities.In the framework of Geo-INQUIRE, Transnational Access (TA, both virtual and on-site) will be provided at six test beds across Europe: the Bedretto Laboratory, Switzerland; the Ella-Link Geolab, Portugal; the Liguria-Nice-Monaco submarine infrastructure, Italy/France; the Irpinia Near-Fault Observatory, Italy; the Eastern Sicily facility, Italy; and the Corinth Rift Laboratory, Greece. These test beds are state-of-the-art research infrastructures, covering the Earth’s surface, subsurface, and marine environments over different spatial scales, from small-scale experiments in laboratories to kilometric submarine fibre cables. The TA will revolve around answering scientific key-questions on the comprehension of fundamental processes associated with geohazards and georesources such as: the preparatory phases of earthquakes, the role of the fluids within the Earth crust, the fluid-solid interaction at the seabed, and the impact of geothermal exploitation. TA will be also offered for software and workflows belonging to the EPOS-ERIC and the ChEESE Centre of Excellence for Exascale in Solid Earth, to develop awarded user’s projects. These are grounded on simulation of seismic waves and rupture dynamics in complex media, tsunamis, subaerial and submarine landslides. HPC-based Probabilistic Tsunami, Seismic and Volcanic Hazard workflows are offered to assess hazard at high-resolution with extensive uncertainty exploration. Support and collaboration will be offered to the awardees to facilitate the access and usage of HPC resources for tackling geoscience problems. Geo-INQUIRE will grant TA to researchers to develop their own lab or numerical experiments with the aim of advancing scientific knowledge of Earth processes while fostering cross-disciplinary research across Europe. To be granted, researchers submit a proposal to the yearly TA calls that will be issued three times during the project life. Calls will be advertised at the Geo-INQUIRE web page https://www.geo-inquire.eu/ and through the existing community channels.To encourage the cross-disciplinary research, Geo-INQUIRE will also organize a series of training and workshops, focused on data, data products and software delivered by research infrastructures, and useful for researchers. In addition, two summer schools will be organized, dedicated to cross-disciplinary interactions of solid earth and marine science.The proposals, for both transnational access and training, will be evaluated by a panel that reviews the technical and scientific feasibility of the project, ensuring equal opportunities and diversity in terms of gender, geographical distribution and career stage. The first call is expected to be issued by the end of Summer 2023. The data and products generated during the TAs will be made available to the scientific community via the project’s strict adherence to FAIR principles.

Probabilistic Tsunami Forecasting (PTF) combines early estimates of earthquake parameters with large ensembles of urgent shallow water tsunami propagation simulations using the GPU based Tsunami-HySEA model (Selva et al. 2021, Nature Commun.). The present version of the PTF is initialised by the earthquake information, but not updated further with new data. In the recently started Horizon Europe project DT-GEO, this PTF is presently being upgraded to a Digital Twin. The first and essential upgrade to realise the Digital Twin is continuous data assimilation enabling a close to real time synthesis of data products and a set of numerical models that allow an updating of the model forecast as new data are continuously assimilated into the model. In DT-GEO, an extended set of data sources, including improved earthquake solutions, sea level tsunami data, and GNSS, will be integrated. A second objective of the PTF is to implement a modularised Digital Twin Component that allows for the inclusion of improved wave and source physics through dispersion, non-hydrostatic tsunami generation, inundation, improved earthquake physics, and cascading earthquake triggered landslide tsunamis. The model will be tested at site demonstrators, in the Mediterranean Sea for eastern Sicily and Samos, and in the Pacific Ocean for Chile and eastern Japan. The presentation will explain how the PTF as it works today, followed by an outline of the design of the components in the Digital Twin, as well as briefly describing initial improvements and plans for further development, including potential integration into Destination Earth. This work is supported by the European Union’s Horizon Europe Research and Innovation Program under grant agreement No 101058129 (DT-GEO, https://dtgeo.eu/).

Modern scientific endeavours already have the capacity to call upon a vast variety of data, often in huge volumes. However, the challenge is not only how to make the most of such a resource, but also how to make it available to the wider scientific community, especially for encouraging curiosity-driven research. Fifty-one institutions from 13 countries are currently working together in the Geo-INQUIRE (Geosphere INfrastructure for QUestions into Integrated REsearch) project.The main goal of this new project is to enhance, give access to, and make interoperable, key datasets of the Geoscience community. This will include "big" data streams and high-performance computing codes which are critical to studying the temporal variation of the solid Earth, forecasting multi-hazards, evaluating Georesources and the analysis of the interface between the solid Earth as well as oceans and atmosphere. About 150 access points – both on-site and virtually are involved. Transnational Access (TA, both virtual and on-site) will be provided at six test beds across Europe: the Bedretto Laboratory, Switzerland; the Ella-Link Geolab, Portugal; the Liguria-Nice-Monaco submarine infrastructure, Italy/France; the Irpinia Near-Fault Observatory, Italy; the Eastern Sicily facility, Italy; and the Corinth Rift Laboratory, Greece.Several European Research Infrastructure Consortia take part, namely the European Plate Observing System (EPOS) for solid Earth and geodynamics observations, the European Multidisciplinary Seafloor and Water Column Observatory (EMSO) for deep-sea and coastal observations, and ECCSEL for CO2 capture, utilization, transport, and storage, and geoenergy. This 16 million Euro project started in October 2022, within the Horizon Europe Infrastructure program of the European Union.The presentation will briefly describe the project and give examples of curiosity-driven research topics which will be made possible through such a multi-disciplinary project. We will finally present the challenges and efforts made to comply with FAIR principles and accompany the dissemination of the data with innovative and cross-disciplinary training activities.

Tsunamis caused by large volcanic eruptions and flanks collapsing into the sea are major hazards for nearby coastal regions. They often occur with little precursory activity and are thus challenging to detect in a timely manner. This makes the pre-emptive identification of volcanoes prone to causing tsunamis particularly important, as it allows for better hazard assessment and denser monitoring in these areas. Here, we present a catalogue of potentially tsunamigenic volcanoes in Southeast Asia and rank these volcanoes by their tsunami hazard. The ranking is based on a multicriteria decision analysis (MCDA) composed of five individually weighted factors impacting flank stability and tsunami hazard. The data are sourced from geological databases, remote sensing data, historical volcano-induced tsunami records, and our topographic analyses, mainly considering the eruptive and tsunami history, elevation relative to the distance from the sea, flank steepness, hydrothermal alteration, and vegetation coverage. Out of 131 analysed volcanoes, we found 19 with particularly high tsunamigenic hazard potential in Indonesia (Anak Krakatau, Batu Tara, Iliwerung, Gamalama, Sangeang Api, Karangetang, Sirung, Wetar, Nila, Ruang, Serua) and Papua New Guinea (Kadovar, Ritter Island, Rabaul, Manam, Langila, Ulawun, Bam) but also in the Philippines (Didicas). While some of these volcanoes, such as Anak Krakatau, are well known for their deadly tsunamis, many others on this list are lesser known and monitored. We further performed tsunami travel time modelling on these high-hazard volcanoes, which indicates that future events could affect large coastal areas in a short time. This highlights the importance of individual tsunami hazard assessment for these volcanoes, the importance of dedicated volcanological monitoring, and the need for increased preparedness on the potentially affected coasts.

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)

Tsunamis caused by large volcanic eruptions and flanks collapsing into the sea are major hazards for nearby coastal regions. They often occur with little precursory activity, and are thus challenging to detect in a timely manner. This makes the pre-emptive identification of volcanoes prone to causing tsunamis particularly important, as it allows for better hazard assessment and denser monitoring in these areas. Here, we present a catalogue of potentially tsunamigenic volcanoes in Southeast Asia and rank these volcanoes by their tsunami hazard. The ranking is based on a Multicriteria Decision Analysis (MCDA) composed of five individually weighted factors impacting flank stability and tsunami hazard. The data is sourced from geological databases, remote sensing data, historical volcano induced tsunami records and our topographic analyses, mainly considering the eruptive and tsunami history, elevation relative to the distance from the sea, flank steepness, hydrothermal alteration as well as vegetation coverage. Out of 131 analysed volcanoes, we found 19 with particularly high tsunamigenic hazard potential in Indonesia (Anak Krakatau, Batu Tara, Iliwerung, Gamalama, Sangeang Api, Karangetang, Sirung, Wetar, Nila, Ruang, Serua) and Papua New Guinea (Kadovar, Ritter Island, Rabaul, Manam, Langila, Ulawun, Bam), but also in the Philippines (Didicas). While some of these volcanoes, such as Anak Krakatau, are well-known for their deadly tsunamis, many others on this list are lesser known and monitored. We further performed tsunami travel time modelling on these high-hazard volcanoes, which indicates that future events could affect large coastal areas in a short time. This highlights the importance of individual tsunami hazard assessment for these volcanoes, dedicated volcanological monitoring, and the need for increased preparedness on the potentially affected coasts.

Landslides begin with the detachment of a mass and end with its impact at lower altitude. To simultaneously model seismic signals produced by these processes, we consider a double single-force model. We applied this source model to the seismic signals generated by the landslide in Uttarakhand, India, on February 7th, 2021. We model the seismic recordings at 12 seismic stations located at less than 100 km epicentral distance. We perform the source inversion by fitting low-frequency (0.08-0.15 Hz) three component (vertical, radial, transversal) waveforms in the time and simultaneously their amplitude spectra in the frequency domain. We compare our results with those obtained for a single-force model applied to each pulse separately. Our results identify two energetic pulses separated by a time delay of ~1 minute. The amplitude of the second pulse, interpreted as the impact, is ~3 times larger than the first one, and of opposite sign. Together with visual observations of the landslide itself, our results confirm that the first pulse was produced by the detachment of the rock mass and the second one by the impact of the mass in the valley. The orientations of the single forces are consistent with the slope geometry and the direction of the debris flow. We discuss statistical measures of fit of the two different inversion approaches and the possible strengths and weaknesses of the new double single-force model.

Extreme climatic and geophysical events pose a threat to societies and have the capacity to cause significant damage and losses whenever they occur, both in their immediate aftermath and in the medium- to long-term. Their consequences can be amplified even further when more than one event affects the same geographical areas within a short time. Be it cascading hazards, in which one event triggers the next, or simply hazards that happen to occur simultaneously (“compounding” hazards), estimation of their cumulative consequences is challenging because the action of one event affects the exposure and vulnerability to the next one. While the efforts from the research community to develop multi-hazard perspectives have increased considerably in recent years, multiple remaining challenges require strongly-coordinated efforts across different disciplines and areas of expertise to tackle them with the most appropriate tools.With a multidisciplinary team of scientists from four different Helmholtz research centres in Germany, we have started working on the CASCO project (2022-2024), in which we will develop an integrated risk workflow for CAScading and COmpounding hazards in COastal urban areas by focusing on a series of events occurring around Mount Etna (Italy). The case-scenario starts with a strong earthquake that triggers a submarine collapse at the eastern flank of Mount Etna, an area already known to be unstable, and both the earthquake and the landslide trigger a tsunami that hits the coasts of Sicily and Calabria. Almost concomitantly, a heatwave or heavy rainfall happens to affect the same regions, further stressing the population that had been affected by the combined effects of the earthquake and tsunami.The project will be directed towards the modelling of the cascading earthquake, landslide and tsunami events, the compounding heatwave and rainfall, as well as their immediate impacts in terms of cumulative damage and casualties. Moreover, the medium- to long-term response in urban dynamics and the effect of these extreme events on the economic development of the affected populations will be explored. By focusing on a tangible scenario, CASCO will not only tackle the challenges associated with bringing together the whole risk chain (which will be valid beyond our case-study) but also produce outcomes that help increase awareness of such extreme events and the need for societies to develop suitable strategies to strengthen their resilience and improve their disaster response.

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.

Controls on the deformation pattern (shortening mode and tectonic style) of orogenic forelands during tectonic shortening remain poorly understood. Here, we use high-resolution 2D thermomechanical ...

When the earthquake rupture is complex and ruptures of multiple fault segments contribute to the total energy release, the produced wavefield is the superposition of individual signals produced by single subevents. Resolving source complexity for large, shallow earthquakes can be used to improve ground shaking and surface slip estimations, and thus tsunami models. The 2018 Mw 7.9 Alaska earthquake showed such complexity: according to previous studies, the rupture initiated as a right-lateral strike-slip fault on a N-S oriented fault plane, but then jumped onto a left-lateral strike-slip fault oriented westward. Rupture complexity and presence of multiple subevents may characterize a number of other earthquakes. However, even when individual subevents are spatially and/or temporally separated, it is very difficult to identify them from far field recordings. In order to model complex earthquakes we have implemented a multiple double couple inversion scheme within Grond, a tool devoted to the robust characterization of earthquake source parameters included in the Pyrocko software. Given the large magnitude of the target earthquake, we perform our source inversions using broadband body waves data (P and S phases) at teleseismic distances. Our approach starts with a standard moment tensor inversion, which allows to get more insights about the centroid location and overall moment release. These values can then be used to constrain the double source inversion. We discuss the performance of the inversion for the Alaska earthquake, using synthetic and real data. First, we generated realistic synthetic waveforms for a two-subevents source, assuming double couple sources with the strike-slip mechanisms proposed for the Alaska earthquake. We model the synthetic dataset both using a moment tensor and a double double couple source, and demonstrate the stability of the double double couple inversion, which is able to reconstruct the two focal mechanisms, the moment ratio and the relative centroid locations of the two subevents. Synthetic tests show that the inversion accuracy can be in some cases reduced, in presence of noisy data and when the interevent time between subevents is short. A larger noise addition affects the retrieval of the focal mechanism orientations only in some cases, but in general all the parameters were well retrieved. Then, we test our tool using real data for the earthquake. The single source inversion shows that the centroid is shifted 27 s in time and 40 km towards NE with respect to the original assumed location retrieved from the gCMT catalogue. The following double double couple source inversion resolves two subevents with right-lateral and left-lateral strike-slip focal mechanisms and Mw 7.6 and 7.8 respectively. The subevent centroids are separated by less than 40 km in space and less than 20 s in time.

39 pages, 6 figures, 2 tables, This paper presents major gaps and challenges for implementing the UN Decade of Ocean Science for Sustainable Development (2021-2030) in the Mediterranean region. The authors make recommendations on the scientific knowledge needs and co-design actions identified during two consultations, part of the Decade preparatory-phase, framing them in the Mediterranean Sea’s unique environmental and socio-economic perspectives. According to the ‘Mediterranean State of the Environment and Development Report 2020’ by the United Nations Environment Programme Mediterranean Action Plan and despite notable progress, the Mediterranean region is not on track to achieve and fully implement the Sustainable Development Goals of Agenda 2030. Key factors are the cumulative effect of multiple human-induced pressures that threaten the ecosystem resources and services in the global change scenario. The basin, identified as a climate change vulnerability hotspot, is exposed to pollution and rising impacts of climate change. This affects mainly the coastal zones, at increasing risk of extreme events and their negative effects of unsustainable management of key economic assets. Transitioning to a sustainable blue economy is the key for the marine environment’s health and the nourishment of future generations. This challenging context, offering the opportunity of enhancing the knowledge to define science-based measures as well as narrowing the gaps between the Northen and Southern shores, calls for a joint (re)action. The paper reviews the state of the art of Mediterranean Sea science knowledge, sets of trends, capacity development needs, specific challenges, and recommendations for each Decade’s societal outcome. In the conclusions, the proposal for a Mediterranean regional programme in the framework of the Ocean Decade is addressed. The core objective relies on integrating and improving the existing ocean-knowledge, Ocean Literacy, and ocean observing capacities building on international cooperation to reach the “Mediterranean Sea that we want”, With the institutional support of the ‘Severo OchoaCentre of Excellence’ accreditation (CEX2019-000928-S)

Probabilistic Tsunami Hazard Analysis (PTHA) quantifies the probability of exceeding a specified inundation intensity at a given location within a given time interval. PTHA provides scientific guidance for tsunami risk analysis and risk management, including coastal planning and early warning. Explicit computation of site-specific PTHA, with an adequate discretization of source scenarios combined with high-resolution numerical inundation modelling, has been out of reach with existing models and computing capabilities, with tens to hundreds of thousands of moderately intensive numerical simulations being required for exhaustive uncertainty quantification. In recent years, more efficient GPU-based High-Performance Computing (HPC) facilities, together with efficient GPU-optimized shallow water type models for simulating tsunami inundation, have now made local long-term hazard assessment feasible. A workflow has been developed with three main stages: 1) Site-specific source selection and discretization, 2) Efficient numerical inundation simulation for each scenario using the GPU-based Tsunami-HySEA numerical tsunami propagation and inundation model using a system of nested topo-bathymetric grids, and 3) Hazard aggregation. We apply this site-specific PTHA workflow here to Catania, Sicily, for tsunamigenic earthquake sources in the Mediterranean. We illustrate the workflows of the PTHA as implemented for High-Performance Computing applications, including preliminary simulations carried out on intermediate scale GPU clusters. We show how the local hazard analysis conducted here produces a more fine-grained assessment than is possible with a regional assessment. However, the new local PTHA indicates somewhat lower probabilities of exceedance for higher maximum inundation heights than the available regional PTHA. The local hazard analysis takes into account small-scale tsunami inundation features and non-linearity which the regional-scale assessment does not incorporate. However, the deterministic inundation simulations neglect some uncertainties stemming from the simplified source treatment and tsunami modelling that are embedded in the regional stochastic approach to inundation height estimation. Further research is needed to quantify the uncertainty associated with numerical inundation modelling and to properly propagate it onto the hazard results, to fully exploit the potential of site-specific hazard assessment based on massive simulations.

Forecasting and early warning systems are important investments to protect lives, properties and livelihood. While early warning systems are frequently used to predict the magnitude, location and timing of potentially damaging events, these systems rarely provide impact estimates, such as the expected amount and distribution of physical damage, human consequences, disruption of services or financial loss. Complementing early warning systems with impact forecasts has a two‐fold advantage: it would provide decision makers with richer information to take informed decisions about emergency measures, and focus the attention of different disciplines on a common target. This would allow capitalizing on synergies between different disciplines and boosting the development of multi‐hazard early warning systems. This review discusses the state‐of‐the‐art in impact forecasting for a wide range of natural hazards. We outline the added value of impact‐based warnings compared to hazard forecasting for the emergency phase, indicate challenges and pitfalls, and synthesize the review results across hazard types most relevant for Europe. Plain language summary Forecasting and early warning systems are important investments to protect lives, properties and livelihood. While such systems are frequently used to predict the magnitude, location and timing of potentially damaging events, they rarely provide impact estimates, such as the expected physical damage, human consequences, disruption of services or financial loss. Extending hazard forecast systems to include impact estimates promises many benefits for the emergency phase, for instance, for organising evacuations. We review and compare the state‐of‐the‐art of impact forcasting across a wide range of natural hazards, and outline opportunities and key challenges for research and development of impact forecasting.

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.

Subduction systems globally terminate, allowing plate-motion to be transferred from the oceanic megathrust onto continental and/or oceanic transform faults. The mechanism of this kinematic transition over earthquake timescales is, however, poorly understood due to a lack of relevant data. Here, we study the 2016 Mw 7.8 Kaikōura Earthquake in New Zealand, the first large instrumentally-recorded earthquake across a subduction-termination, to investigate this transfer mechanism in detail. We find that the Kaikōura Earthquake, unlike standard subduction earthquakes globally, involved a predominance (∼80%) of coseismic-slip on upper-plate faults and minor triggered-slip on the underlying oceanic subduction-thrust. In the months following the earthquake, the subduction-thrust accommodated most of the earthquake's afterslip down-dip of its co-seismic rupture zone. This top-down strain-release mechanism is in accord with local geological, geodetic and historical seismicity data which suggest that the bulk of plate-convergence (>75%) is accommodated in the upper-plate. We suggest, therefore, that this alternative strain-release mechanism, which is distinct from standard plate-boundary earthquakes, is characteristic of subduction-terminations and results in the majority of seismic/tsunami hazard being on steep near-surface faults.

The tsunami disasters of 2004 in the Indian Ocean and of 2011 along the Tohoku coast of Japan revealed severe gaps between the anticipated risk and consequences, with resulting loss of life and property. A similar observation is also relevant for the smaller, yet disastrous, tsunamis with unusual source characteristics such as the recent events in Palu Bay and Sunda Strait in 2018. The severe consequences were underestimated in part due to the lack of rigorous and accepted hazard analysis methods and large uncertainty in forecasting the tsunami sources. Population response to small recent tsunamis in the Mediterranean also revealed a lack of preparedness and awareness. While there is no absolute protection against large tsunamis, a more accurate analysis of the potential risk can help to minimize losses. The tsunami community has made significant progress in understanding tsunami hazard from seismic sources. However, this is only part of the inputs needed to effectively manage tsunami risk, which should be understood more holistically, including non-seismic sources, vulnerability in different dimensions and the overall societal effects, in addition to its interaction with other hazards and cascading effects. Moreover, higher standards need to be achieved to manage and quantify uncertainty, which govern our basis for tsunami risk decision making. Hence, a collective community effort is needed to effectively handle all these challenges across disciplines and trades, from researchers to stakeholders. To coordinate and streamline these activities and make progress towards implementing the Sendai Framework of Disaster Risk Reduction (SFDRR) the Global Tsunami Model network (GTM) was initiated in 2015 towards enhancing our understanding of tsunami hazard and risk from a local to global scale. Here, we focus on coordinated European efforts, sharing the same goals as GTM, towards improving standards and best practices for tsunami risk reduction. The networking initiative, AGITHAR (Accelerating Global science In Tsunami HAzard and Risk Analysis), is a European COST Action, aims to assess, benchmark, improve, and document methods to analyse tsunami hazard and risk, understand and communicate the uncertainty involved, and interact with stakeholders in order to understand the societal needs and thus contribute to their effort to minimize losses. In this presentation, we provide an overview of the suite of methodologies used for tsunami hazard and risk analysis, review state of the art in global tsunami hazard and risk analysis, dating back to results from the Global Risk Model in 2015, and highlight possible gaps and challenges. We further discuss how AGITHAR and GTM will address how to tackle these challenges, and finally, discuss how global and regional structures such as the European Plate Observing System (EPOS) and the UNDRR Global Risk Assessment Framework (GRAF) can facilitate and mutually benefit towards an integrated framework of services aiding improved understanding of multiple hazards.

Historical data indicate that the Middle America subduction zone represents the primary tsunamigenic source that affects the Central American coastal areas. In recent years, the tsunami potential in the region has mainly been assessed using maximum credible earthquakes or historical events showing moderate tsunami potential. However, such deterministic scenarios are not provided with their adequate probability of occurrence. In this study, earthquake rates have been combined with tsunami numerical modeling in order to assess probabilistic tsunami hazard posed by local and regional seismic sources. The common conceptual framework for the probabilistic seismic hazard assessment has been adapted to estimate the probabilities of exceeding certain tsunami amplitudes along the Central American Pacific coast. The study area encompasses seismic sources related to the Central America, Colombia and Ecuador subduction zones. In addition to the classical subduction inter-plate events, this study also incorporates sources at the outer rise, within the Caribbean crust as well as intraslab sources. The study yields conclusive remarks showing that the highest hazard is posed to northwestern Costa Rica, El Salvador and the Nicaraguan coast, southern Colombia and northern Ecuador. In most of the region it is 50 to 80% likely that the tsunami heights will exceed 2 m for the 500 year time exposure (T). The lowest hazard appears to be in the inner part of the Fonseca Gulf, Honduras. We also show the large dependence of PTHA on model assumptions. While the approach taken in this study represents a thorough step forward in tsunami hazard assessment in the region, we also highlight that the integration of all possible uncertainties will be necessary to generate rigorous hazard models required for risk planning.

Applying probabilistic methods to infrequent but devastating natural events is intrinsically challenging. For tsunami analyses, a suite of geophysical assessments should be in principle evaluated because of the different causes generating tsunamis (earthquakes, landslides, volcanic activity, meteorological events, and asteroid impacts) with varying mean recurrence rates. Probabilistic Tsunami Hazard Analyses (PTHAs) are conducted in different areas of the world at global, regional, and local scales with the aim of understanding tsunami hazard to inform tsunami risk reduction activities. PTHAs enhance knowledge of the potential tsunamigenic threat by estimating the probability of exceeding specific levels of tsunami intensity metrics (e.g., run-up or maximum inundation heights) within a certain period of time (exposure time) at given locations (target sites); these estimates can be summarized in hazard maps or hazard curves. This discussion presents a broad overview of PTHA, including (i) sources and mechanisms of tsunami generation, emphasizing the variety and complexity of the tsunami sources and their generation mechanisms, (ii) developments in modeling the propagation and impact of tsunami waves, and (iii) statistical procedures for tsunami hazard estimates that include the associated epistemic and aleatoric uncertainties. Key elements in understanding the potential tsunami hazard are discussed, in light of the rapid development of PTHA methods during the last decade and the globally distributed applications, including the importance of considering multiple sources, their relative intensities, probabilities of occurrence, and uncertainties in an integrated and consistent probabilistic framework.

This study reports on recent developments of the Indonesia Tsunami Early Warning System (InaTEWS), specifically the tsunami modeling components used in the system. It is a dual system: firstly, InaTEWS operates a high-resolution scenario database pre-computed with the finite element model TsunAWI; running in parallel, the system also contains a supra real-time modeling component based on the GPU-parallelized linear long-wave model easyWave, capable of dealing with events outside the database coverage. The evolution of the tsunami scenario database over time is covered in the first sections also touching on the involved capacity building efforts. Starting with a coverage of just the Sunda Arc region, the database now includes scenarios for 15 fault zones. The study is augmented by an investigation of warning products used for early warning; the estimated wave height (EWH) and the estimated time of arrival (ETA). These parameters are determined by easyWave and TsunAWI with model specific approaches. Since the numerical setup of the two models is very different, the extent of variations in warning products is investigated for a number of scenarios, where both pure database scenarios and applications to real events are considered. Finally, the performance of the system in past tsunami events is reviewed to point out major system updates.

Historical data indicate that the Middle America subduction zone represents the primary tsunamigenic source that affects the Central American coastal areas. In recent years, the tsunami potential in the region has mainly been assessed using maximum credible earthquakes or historical events showing moderate tsunami potential. However, such deterministic scenarios are not provided with their adequate probability of occurrence. In this study, earthquake rates have been combined with tsunami numerical modeling in order to assess probabilistic tsunami hazard posed by local and regional seismic sources. The common conceptual framework for the probabilistic seismic hazard assessment has been adapted to estimate the probabilities of exceeding certain tsunami amplitudes along the Central American Pacific coast. The study area encompasses seismic sources related to the Central America, Colombia and Ecuador subduction zones. In addition to the classical subduction inter-plate events, this study also incorporates sources at the outer rise, within the Caribbean crust as well as intra-slab sources. The study yields conclusive remarks showing that the highest hazard is posed to northwestern Costa Rica, El Salvador and the Nicaraguan coast, southern Colombia and northern Ecuador. In most of the region it is 50 to 80% likely that the tsunami heights will exceed 2 m for the 500 year time exposure (T). The lowest hazard appears to be in the inner part of the Fonseca Gulf, Honduras. We also show the large dependence of PTHA on model assumptions. While the approach taken in this study represents a thorough step forward in tsunami hazard assessment in the region, we also highlight that the integration of all possible uncertainties will be necessary to generate rigorous hazard models required for risk planning.

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-East Atlantic, the Mediterranean, and connected Seas (NEAM). In this online data product, the hazard results are provided by hazard curves calculated at 2,343 Points of Interest (POI), distributed in the North-East Atlantic (1,076 POIs), the Mediterranean Sea (1,130 POIs), and the Black Sea (137 POIs) at an average spacing of ~20 km. For each POI, hazard curves are given for the mean, 2nd, 16th, 50th, 84th, and 98th percentiles. Maps derived from hazard curves are Probability maps for Maximum Inundation Heights (MIH) of 1, 2, 5, 10, 20 meters; Hazard maps for Average Return Periods (ARP) of 500, 1,000, 2,500, 5,000, 10,000 years. For each map, precalculated displays are provided for the mean, the 16th percentile, and the 84th percentile. All data are also made accessible through an interactive web mapper and through Open Geospatial Consortium standard protocols. The model was prepared in the framework of the European Project TSUMAPS-NEAM (http://www.tsumaps-neam.eu/) funded by the mechanism of the European Civil Protection and Humanitarian Aid Operations (grant no. ECHO/SUB/2015/718568/PREV26)., Co-financed by the European-Union Civil Protection Mechanism Agreement Number: ECHO/SUB/2015/718568/PREV26

Flank instability and sector collapses, which pose major threats, are common on volcanic islands. On 22 Dec 2018, a sector collapse event occurred at Anak Krakatau volcano in the Sunda Strait, triggering a deadly tsunami. Here we use multiparametric ground-based and space-borne data to show that prior to its collapse, the volcano exhibited an elevated state of activity, including precursory thermal anomalies, an increase in the island’s surface area, and a gradual seaward motion of its southwestern flank on a dipping décollement. Two minutes after a small earthquake, seismic signals characterize the collapse of the volcano’s flank at 13:55 UTC. This sector collapse decapitated the cone-shaped edifice and triggered a tsunami that caused 430 fatalities. We discuss the nature of the precursor processes underpinning the collapse that culminated in a complex hazard cascade with important implications for the early detection of potential flank instability at other volcanoes., On 22 December 2018, the western flank of Anak Krakatau collapsed into the sea of the Sunda Strait triggering a tsunami which killed approximately 430 people and displaced 33,000. Here, the authors show that Anak Krakatau exhibited an elevated state of activity several months prior to the collapse, including precursory thermal anomalies, an increase in the island’s surface area, and a gradual seaward motion of the southwestern flank.

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.

The Middle America convergent margin is capable of generating tsunamigenic earthquakes as portrayed through historical data. However, historical tsunamigenic sources appear to be unevenly distributed along the Middle America convergent margin. The last significant tsunami generated along the southern segments of this convergent margin occurred in 1992 causing 10-m run-up in Nicaragua. In particular, the Costa Rican Pacific coast has not clear evidences of significant wave heights or local tsunami events. Partly this could be explained by the lack of historical data since paleotsunami information is not available; thus, no tsunami data beyond 500 years could be found. Nevertheless, large rupture areas or shallow ruptures offshore Costa Rica cannot be neglected. We hereby present the results of numerical simulations aimed to assess the tsunamigenic potential of seismic sources offshore southern Central America. Our approach follows the deterministic method where earthquake maximum credible scenarios are assessed. We assumed hundreds of fault ruptures to obtain surface displacement and tsunami wave heights at the shoreline to assess the tsunami threat along south Nicaragua and Costa Rica. These fault scenarios are based on historical earthquakes as well as on the newest geodetic studies that provided slip deficit. Shallow part of locked segments in this region could release accumulated stress in the near future and may generate tsunamis similar to the 1992 Nicaragua M w 7.6 and 2012 El Salvador M w 7.4 tsunami earthquakes. The numerical simulations presented in this study show that maximum wave heights up to 2–6 m could affect the southern Pacific coast of Central America even when moderate-to-large M w 7.0–8.0 interplate or outer rise earthquakes occur.

On 16 September 2015, the M-W = 8.2 Illapel megathrust earthquake ruptured the Central Chilean margin. Combining inversions of displacement measurements and seismic waveforms with high frequency (HF) teleseismic backprojection, we derive a comprehensive description of the rupture, which also predicts deep ocean tsunami wave heights. We further determine moment tensors and obtain accurate depth estimates for the aftershock sequence. The earthquake nucleated near the coast but then propagated to the north and updip, attaining a peak slip of 5-6 m. In contrast, HF seismic radiation is mostly emitted downdip of the region of intense slip and arrests earlier than the long period rupture, indicating smooth slip along the shallow plate interface in the final phase. A superficially similar earthquake in 1943 with a similar aftershock zone had a much shorter source time function, which matches the duration of HF seismic radiation in the recent event, indicating that the 1943 event lacked the shallow slip.

Global Positioning System (GPS) has been proved to be a powerful tool for measuring co-seismic ground displacements with an application to seismic source inversion. Whereas most of the tsunamis are triggered by large earthquakes, GPS can contribute to the tsunami early warning system (TEWS) by helping to obtain tsunami source parameters in near real-time. Toward the end of 2012, the second phase of the BeiDou Navigation Satellite System (BDS) constellation was accomplished, and BDS has been providing regional positioning service since then. Numerical results indicate that precision of BDS nowadays is equivalent to that of the GPS. Compared with a single Global Satellite Navigation System (GNSS), combined BDS/GPS real-time processing can improve accuracy and especially reliability of retrieved co-seismic displacements. In the present study, we investigate the potential of BDS to serve for the early warning system of tsunamis in the South China Sea region. To facilitate early warnings of tsunamis and forecasting capabilities in this region, we propose to distribute an array of BDS-stations along the Luzon Island (Philippines). By simulating an earthquake with Mw = 8 at the Manila trench as an example, we demonstrate that such an array will be able to detect earthquake parameters in real time with a high degree of accuracy and, hence, contribute to the fast and reliable tsunami early warning system in this region.

Currently, China is strengthening infrastructure construction in South China Sea (SCS) area, it is believed that in the near future, SCS territory will play an increasingly important role in political and economic exploitation as well as international cooperation issues. However, political and economical development plans must also consider possible impact of various natural hazards as well as strategies for damage minimization. In particular, the Manila subduction zone has been identified as a zone of potential tsunami hazard in the SCS region. Despite the fact, that no earthquake higher than Mw7.6 has occurred in this region during last years, events of higher magnitudes cannot be excluded which could expose SCS to the tsunami risk. Recently, development of a region-wide Tsunami Early Warning System (TEWS) was proposed. Traditional TEWS incorporates seismic network for detection of earthquakes and estimating their tsunamigenic potential combined with ocean-based observations for monitoring of wave propagation. After the striking Great Sumatra 2004 earthquake and tsunami, new technologies were introduced into the tsunami early warning. In particular, continuous GPS-arrays able to measure earthquake source parameters in almost real-time. Around the end of 2012, China has completed the constellation of the regional Beidou Navigation Satellite System (BDS), and this system is running operationally since that time. Nowadays, China government is reinforcing the application of BDS in scientific and engineering fields. Numerous studies have demonstrated that the precise positioning performance of BDS is equal in match against GPS in the Asian-Pacific region, which promisingly suggests BDS can be implemented to facilitate tsunami early warning in the South China Sea. In present study, we use numerical simulation technique to assess the potential of the BDS for the tsunami early warning in SCS. We suggest deployment of a real-time BDS coastal network at the Luzon Island. Such a network, if established, will effectively contribute to fast source inversion at the Manila trench and, hence, to more rapid and reliable tsunami early warning in the whole SCS region.

Nowadays, Global Positioning System (GPS) plays an increasingly important role in retrieving real-time precise co-seismic displacements for geo-hazard monitoring and early warning. Several real-time positioning approaches have been demonstrated for such purpose, such as real-time kinematic relative positioning, precise point positioning, etc., where dual-frequency geodetic receivers are applied for the removal of ionosphere delays by inter-frequency combination. At the same time, it would be also useful to develop efficient algorithms for estimating precise displacements with low-cost GPS receivers since they can make a denser network or multi-sensors combination without putting too much financial burden. In this contribution, we present a new method to retrieve precise co-seismic displacements in real-time using a standalone single-frequency receiver. In the new method, observations prior to an earthquake are utilized to establish a precise ionospheric delay prediction model, so that precise co-seismic displacements can be obtained without any convergence process. Our method was validated with an outdoor experiment as well as by re-processing of 1-Hz GPS data collected by the GEONET network during the 2011 Tohoku Mw 9.0 earthquake. For the latter, RMS against dual-frequency receivers constituted 2 cm for horizontal components and 3 cm for the vertical component. We specially address the observation biases and their impact on the accuracy of single frequency positioning. Our approach makes real-time GPS displacement monitoring with dense network much more affordable in terms of financial costs.