Dr. Aranguiz-Muñoz, Rafael

2016, Dr. Aranguiz-Muñoz, Rafael, González, Gabriel, González, Juan, Catalán, Patricio, Cienfuegos, Rodrigo, Yagi, Yuji, Okuwaki, Ryo, Urra, Luisa, Contreras, Karla, Del Rio, Ian, Rojas, Camilo

On September 16, 2015 a magnitude Mw 8.3 earthquake took place off the coast of the Coquimbo Region, Chile. Three tsunami survey teams covered approximately 700 km of the Pacific coast. The teams surveyed the area, recording 83 tsunami flow depth and runup measurements. The maximum runup was found to be 10.8 m at only one small bay, in front of the inferred tsunami source area. However, it was observed that runup in other locations rarely exceed 6 m. Tsunami runup was larger than those of the 2014 Pisagua event, despite the similar earthquake magnitude. Moreover, tsunami arrival times were found to be shorter than those of previous tsunamis along the Chilean subduction zone. Numerical simulations of the tsunami event showed a good agreement with field data, highlighting that tsunami arrival time and the spatial variation of the tsunami amplitudes were strongly influenced by the bathymetry, coastal morphology and the slip distribution of the causative earthquake.

2020, Dr. Aranguiz-Muñoz, Rafael, Dr. Benavente-Bravo, Roberto, Becerra, Ignacio, González, Juan

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.

2016, Dr. Aranguiz-Muñoz, Rafael, Okuwaki, Ryo, Yagi, Yuji, González, Juan, González, Gabriel

We constructed a seismic source model for the 2015 MW 8.3 Illapel, Chile earthquake, which was carried out with the kinematic waveform inversion method adopting a novel inversion formulation that takes into account the uncertainty in the Green’s function, together with the hybrid backprojection method enabling us to track the spatiotemporal distribution of high-frequency (0.3–2.0 Hz) sources at high resolution by using globally observed teleseismic P-waveforms. A maximum slip amounted to 10.4 m in the shallow part of the seismic source region centered 72 km northwest of the epicenter and generated a following tsunami inundated along the coast. In a gross sense, the rupture front propagated almost unilaterally to northward from the hypocenter at \2 km/s, however, in detail the spatiotemporal slip distribution also showed a complex rupture propagation pattern: two up-dip rupture propagation episodes, and a secondary rupture episode may have been triggered by the strong high-frequency radiation event at the down-dip edge of the seismic source region. High-frequency sources tends to be distributed at deeper parts of the slip area, a pattern also documented in other subduction zone megathrust earthquakes that may reflect the heterogeneous distribution of fracture energy or stress drop along the fault. The weak excitation of high-frequency radiation at the termination of rupture may represent the gradual deceleration of rupture velocity at the transition zone of frictional property or stress state between the megathrust rupture zone and the swarm area.

2020, González, Juan, González, Gabriel, Aranguiz-Muñoz, Rafael, Melgar, Diego, Zamora, Natalia, Shrivastava, Mahesh N., Das, Ranjit, Catalán, Patricio A., Cienfuegos, Rodrigo

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.

2015, Catalán, Patricio, Aranguiz-Muñoz, Rafael, González, Gabriel, Tomita, Takashi, Cienfuegos, Rodrigo, González, Juan, Shrivastava, Mahesh N., Kumagai, Kentaro, Mokrani, Cyril, Cortés, Pablo, Gubler, Alejandra

On 1 April 2014, an earthquake with moment magnitudeMw8.2 occurred off the coast ofnorthern Chile, generating a tsunami that prompted evacuation along the Chilean coast. Here tsunamicharacteristics are analyzed through a combination of field data and numerical modeling. Despite theearthquake magnitude, the tsunami was moderate, with a relatively uniform distribution of runup, whichpeaked at 4.6 m. This is explained by a concentrated maximal slip at intermediate depth on the megathrust,resulting in a rapid decay of tsunami energy. The tsunami temporal evolution varied, with locations showingsustained tsunami energy, while others showed increased tsunami energy at different times after theearthquake. These are the result of the interaction of long period standing oscillations and trapped edgewave activity controlled by inner shelf slopes. Understanding these processes is relevant for the region,which still posses a significant tsunamigenic potential

2020, Dr. Aranguiz-Muñoz, Rafael, Esteban, Miguel, Takagi, Hiroshi, Mikami, Takahito, Takabatake, Tomoyuki, Gómez, Matías, González, Juan, Shibayama, Tomoya, Okuwaki, Ryo, Yagi, Yuji, Shimizu, Kousuke, Achiari, Hendra, Stolle, Jacob, Robertson, Ian, Ohira, Koichiro, Nakamura, Ryota, Nishida, Yuta, Krautwald, Clemens, Goseberg, Nils, Nistor, Ioan

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.