Dr. Aránguiz-Muñoz, Rafael

Central Chile is exposed to tsunami hazard, and large earthquakes and tsunamis have occurred over the last 500 years. Tsunami hazard analysis in Chile has been traditionally implemented by means of a deterministic approach, which is based on historical events and uniform slip distribution. The objective of the present study is to improve tsunami hazard analysis in central Chile (30°S to 38°S). To encompass the purpose, stochastic earthquake scenarios of magnitude Mw 8.8 to 9.2 were generated. Two different sets of stochastic tsunami scenarios were selected by means of the Stochastic Reduced Order Model (SROM), which were applied to Quintero bay to perform a Probabilistic Tsunami Hazard Analysis (PTHA). The results showed that PTHA of Quintero bay from stochastic tsunami scenarios agrees with paleotsunami records in the bay, while a deterministic tsunami scenario underestimated the hazard. Two sets (50 and 100 scenarios, respectively) give similar results when smaller return periods are analyzed. However, for larger return periods (Unknown node type: font 2000 yr) the set of 100 scenarios show better results consistent with previous paleoseismological findings. The methodology implemented here can be replicated in other seismic regions in Chile as well as in other active subduction zones, thus, both near field and far field events can be analyzed.

The southern Peru and northern Chile coastal region is an active subduction zone that contains one of the most signifcant seismic gaps in the eastern Pacifc basin (~17°S–~24°S). Although the gap was partially flled by the 2014 Mw 8.1 Iquique earthquake, there is still a high seismogenic potential to release a Mw ~9 earthquake in the near future; therefore, all the near-feld coastal cities in the region face a latent tsunami threat. In this article, we propose a hybrid deterministic–stochastic multi-scenario approach to assess the current tsunami hazard level in the city of Iquique, an important commercial and industrial center of northern Chile that is home to 184,000 inhabitants. In our approach, we defned 400 stochastic, 10 deterministic and 10 homogeneous tsunamigenic earthquake scenarios, covering the entire area of the seismic gap. Based on the regional distribution of gravity anomalies and published interseismic coupling distributions, we interpreted the occurrence of four major asperities in the subduction interface of the seismic gap. The asperity pattern was used to construct a group of deterministic slip-defcit earthquake sources with seismic magnitudes ranging between Mw 8.4 and Mw 8.9. Additionally, we constructed 10 homogeneous slip scenarios to generate an inundation baseline for the tsunami hazard. Subsequently, following a stochastic scheme, we implemented a Karhunen–Loève expansion to generate 400 stochastic earthquake scenarios within the same magnitude range as the deterministic slip-defcit sources. All sources were used as earthquake scenarios to simulate the tsunami propagation and inundation by means of a non-hydrostatic model (Neowave 2D) with a classical nesting scheme for the city of Iquique. We obtained high-resolution data for fow depth, coastal surface currents and sea level elevation. The results suggest that the peak slip location and shelf resonance play an important role in the calculated coastal fow depths. The analysis of the entire set of simulated stochastic earthquake scenarios indicates that the worst-case scenario for Iquique is a Mw 8.9 earthquake. This scenario presented a tsunami arrival time of ~12 min, which is critical for the evacuation process. In addition, the maximum wave height and tsunami fow depth were found to be ~10 m and ~24 m, respectively. The observed coastal resonance processes exhibit at least three destructive tsunami wave trains. Based on historical and instrumental catalog statistics, the recurrence time of the credible worst-case earthquake scenario for Iquique (Mw 8.9) is 395 years, with a probability of occurrence of ~11.86% in the next 50 years.

The September 28 2018 Palu tsunami surprised the scientific community, as neither the earthquake magnitude nor its strike-slip mechanism were deemed capable of producing the wave heights that were observed. However, recent research has shown that the earthquake generated several landslides inside Palu bay. The authors conducted a post-disaster field survey of the area affected to collect spatial data on tsunami inundation heights, nearshore and bay bathymetry, and carried out eyewitness interviews to collect testimonies of the event. In addition, numerical simulations of the tsunami generation and propagation mechanisms were carried out and validated with the inferred time series. Seven small submarine landslides were identified along the western shore of the bay, and one large one was reported on the eastern shore of Palu City. Most of these landslides occurred at river mouths and reclamation areas, where soft submarine sediments had accumulated. The numerical simulations support a scenario in which the tsunami waves that arrived at Palu city 4–10 min after the earthquake were caused by the co-seismic seafloor deformation, possibly coupled with secondary waves generated from several submarine landslides. These findings suggest that more comprehensive methodologies and tools need to be used when assessing probabilistic tsunami hazards in narrow bays.