The Andean Foreland Basin in the Chaco-Pampean Plain of North-West Argentina is thought to have been tectonically inactive during the Cenozoic. However, re-interpreted industry seismic-reflection data and borehole information document a complex tectonic history in the subsurface at least until Palaeogene times. Data synopsis and re-analysis reveal two regionally extensive and approximately NW-SE-oriented basement highs beneath the flat present-day surface, the Quirquincho (or Rincón Caburé) High and the Pampeano-Chaqueño High. These large geological structures were described previously, but the mechanism that elevated these features relative to the surrounding stratigraphy and the timing of uplift has remained elusive.

This study documents and describes of the Quirquincho and Pampeano-Chaqueño Highs and their relationship to the depocenters around and the sedimentary successions of the Chacoparanaense, Salta Rift and Andean Foreland Basins. We studied palaeo-basin morphology, stratigraphy, stratal terminations and distance to the Andes to unravel viable mechanisms that influenced the genesis of the two structural highs. The work presented is based on the re-interpretation of a large set of subsurface information (2D seismic-reflection profiles and well reports) using a regional approach that considers the Chacoparanaense, Salta Rift and Andean Foreland Basins as a complex lateral arrangement of basins varying in activity through time, depending on their relative location concerning orogens and rifted ocean margins.

Our research reveals that the Quirquincho and Pampeano-Chaqueño Highs were elevated features from the Late Palaeozoic to the Palaeogene, strongly influencing the deposition of Mesozoic and Palaeogene sediments. The tectonic mechanism controlling the rise of the Quirquincho and Pampeano-Chaqueño Highs was initially flexural deformation in the foreland of the Gondwanide orogen in the Permian, subsequently influenced in the Mesozoic by the opening of the Southern Atlantic Ocean.

How to cite: Cortassa, V., Ondrak, R., Back, S., del Papa, C., and Rossello, E.: Quirquincho and Pampeano-Chaqueño Highs: forebulge or Palaeozoic structures?, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-990, https://doi.org/10.5194/egusphere-egu24-990, 2024.

The tectonic setting of the Northern Andes is delineated fundamentally by a western oceanic terrane that was juxtaposed to the continental margin along the now fossilized interandean Romeral suture since the Early Cretaceous. This constellation and the connection to the Caribbean Large Igneous province have been attributed to a far-travelled and now partially subducted, formerly coherent terrane with a trailing edge represented by the Panama-Choco block. A former disconnection between oceanic terrane and South American plate may, however, be contended by considering continental provenance data of siliciclastic and volcanic rock units and a widely distributed geochemical arc signature of the effusive rock series. Moreover, the emplacement of the basic igneous sequences was strongly controlled by extensional tectonics and subduction correlates in its lifetime with the production of oceanic crust, suggesting a coupling between intrusive activity and convergence. In this contribution, we examine apparently conflicting structural deformations that may be reconciled, however, with the opening of a forearc basin and a deformational imprint that affected extensively the continental margin, supposing the existence of a subduction system composed of two continentward dipping slabs.

How to cite: Kammer, A., Betancur, C. A., and Conde, C.: Lower plate retreat and opening of a Cretaceous forearc basin, Northern Andes, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12834, https://doi.org/10.5194/egusphere-egu24-12834, 2024.

The Central and South-Central Andes form a “two-sided” mountain belt bounded by distinct zones of convergence in the western forearc and eastern foreland flanks. Previous geodetic studies of interseismic deformation in the Bolivian Subandes and the Argentine Precordillera found that the forearc to foreland velocity field decayed too slowly to be explained purely by elastic shortening driven by locking of the Nazca megathrust. The velocity field is more precisely explained if elastic deformation is augmented by eastward displacement of the entire Andes. Here, we extend the earlier interpretation of interseismic motion and argue that foreland décollements can participate in the co- and postseismic phases of the earthquake deformation cycle associated with the Nazca megathrust. These findings have direct implications in estimating recurrence interval, slip rate, and probabilistic seismic hazard analysis on both sides of the orogen.

How to cite: Figueroa, M. A., Sobrero, F. S., Gómez, D. D., Smalley Jr., R., Bevis, M. G., Caccamise II, D. J., and Kendrick, E.: Microplate behaviour of the Andes during co- and early postseismic phases of the seismic cycle, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12909, https://doi.org/10.5194/egusphere-egu24-12909, 2024.

The South American active margin is one of the most important and well-studied subduction zones on the Earth. In the last decade, our knowledge about the geometry of its respective constituent parts in Northern Chile has been significantly expanded thanks to seismicity data from a large network of permanent seismic stations. That calls for an effort to summarize these diverse constraints in a single 3D model, which can be validated and optimized using satellite gravity data.

Integrated geophysical modelling is primarily based on gravity data, which bring information about different crustal density inhomogeneities and their sources. As an inverse geophysical problem, gravity modelling is ambiguous, and therefore it is necessary to include geometry constraints from other geophysical data as well as geological information. Different types of seismic data offer the most commonly used constraints due to their depth range and the relatively well-described relation between seismic velocities and densities.

We aim to compile a 3D integrated geophysical model for Northern Chile in the IGMAS+ software, based on gravity data of the global gravitational (or geopotential) model EIGEN-6C4, which include terrestrial, satellite and altimetry data to a high degree and order of spherical harmonic expansion. As the main geometry constraints of the model, we use the newest available seismicity catalogs in the study area together with crustal thickness values from receiver functions. Starting with the geometries of previously published density models in the area of the Central Andes, especially those located at least partly in Northern Chile, we will modify these models guided by the geometry constraints from seismic data and will also use regional and global models of crustal and lithospheric interfaces, such as the top of basement in sedimentary areas, plate interface geometry, depth to continental Moho, and the lithosphere-asthenosphere boundary (LAB). Densities of the modeled bodies will be selected based on previously published models or estimated from seismic velocities. We also plan to study the gravity effect of the different geometries deduced from different generations of seismic data.

Our contribution provides an overview of evidence compiled in previous studies and adds new information on the deep lithospheric structure of the North Chilean margin by integrating them into a single model.

How to cite: Godová, D. and Sippl, C.: Towards a 3D integrated geophysical model of Northern Chile, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8006, https://doi.org/10.5194/egusphere-egu24-8006, 2024.

Repeating earthquakes, corresponding to the rupture of the same asperity over time at more or less regular intervals, can be used to estimate the slip rate on a subduction plate interface. The purpose of this work is to build a catalog of repeaters for the central part of the Chilean subduction zone, extending from 24°S to 33°S latitude and centered on the Copiapo seismic gap. As a basis for our study, we used the seismicity catalog from the Centro Sismológico Nacional (CSN).

The similarity between waveforms gives a good criterion to assign earthquakes to a similar asperity. To measure it, we calculated for each pair of events the coherency, correlation and associated time lag between the vertical components of their P waves, on a 5 s window starting 1 s before the P arrival. We tested different frequency bands (1-4 Hz, 3-12 Hz, 5-20 Hz, ...) and kept each time the one with the best coherency value. We selected all pairs of earthquakes with coherency higher than 0.95, at three or more stations. To ensure a stable measurement, we imposed that the time lags from cross-correlation and coherency differ by less than 0.01 s. To verify that the earthquakes of a repeaters family take place at the same location, we relocated the events using a double-difference method and created clusters based on both coherency and location similarity. As coherency values are calculated on a 5 s window, we relocate the centroids of the events, i.e. the center of mass of the rupture. To estimate the surface rupture, we calculate the rupture radius based on the seismic moment and stress-drop values (Eshelby 1957), estimating the stress-drop and seismic moment with SourceSpec (Satriano 2023).          

As preliminary results, we found 347 families, mostly located between 30-60 km deep, and between 29-33.5°S. Almost no repeaters were found before 2015 due to the lack of available stations. Obtained families contain a few events, with 11 earthquakes in the biggest one. We compared the obtained repeaters with the coupling along the plate interface.  Most repeaters are located at the transition between strong and low coupling zones in the Illapel area, making a circle shape around the deep part of the Illapel coseismic slip. Furthermore, we investigated the evolution of the number of repeaters with time in different areas and found potential aseismic slip marked by repeaters' activity consistent with previous observations, such as before the 2017 Valparaiso sequence (Ruiz et al., 2017) or after the Atacama sequence (Klein et al., 2021).

In order to obtain more complete repeaters families, we created a new machine learning based earthquake catalog for the study area with SeisBench (Woollam et al., 2022), using data from permanent and temporary networks in Northern and Central Chile. We are currently applying our analysis to this new catalog.

How to cite: Chouli, A., Costes, L., Marsan, D., Münchmeyer, J., Giffard-Roisin, S., and Socquet, A.: Search for repeaters in the central part of the Chilean subduction zone , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17042, https://doi.org/10.5194/egusphere-egu24-17042, 2024.

The use of machine-learning based tools for phase picking and association is in the process of revolutionizing the field of seismicity analysis, leading to the simplified creation of “deep” seismicity catalogs often containing 10s or 100s of thousands of events. Having such catalogs as an available resource opens the field for novel ways of combining seismicity information with other types of datasets. At the same time, the sheer amount of data poses challenges for the visualization as well as joint analysis with other constraints.

In this contribution, I want to explore different ways of using and visualizing large seismicity catalogs, using a range of different recently compiled earthquake catalogs from the Chilean margin as showcase examples. Moreover, I attempt to find efficient ways of cross-plotting seismicity data from “deep” catalogs with other datasets such as interplate coupling maps, seismic velocity distributions, temperature models or inferred mineralogy maps.

How to cite: Sippl, C.: Using large microseismicity catalogs to constrain subduction zone processes – examples from the Chilean margin, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7484, https://doi.org/10.5194/egusphere-egu24-7484, 2024.

Northern Chile is one of the most seismogenic regions on the planet, and has been monitored by a permanent network of seismic stations since 2007. We here present a first step towards a new, more complete seismicity catalog for this region, leveraging modern deep-learning based algorithms for phase picking and association.

We first assessed the performance of EQTransformer, a deep learning based phase picker, in detecting and phase picking seismic data from the Northern Chile Subduction Zone by comparison with a large, meticulously handpicked dataset. We found that the "INSTANCE" model within SeisBench yielded the best performance for our study area. Through systematic threshold variations, we determined the optimal values using Precision-Recall curves (0.4 for event detection, 0.1 for P and S picks). Subsequently, we applied GaMMA, identified as the best performing phase associator in synthetic tests, coupled with NonLinLoc for initial event location. One of GaMMA's key operational criteria is the association threshold, where we required a minimum of five seismic phases to define an event, which yielded a high reliability in the phase association process. Moreover, we refined the catalog by automatically identifying and removing duplicate events. All associated events were consecutively relocated in a 1D and a 2D velocity model, using the VELEST and simul2000 algorithms. Events with disproportionally high RMS residuals as well as single picks with high residuals were removed in the process. In a final step, events were relocated with a double-difference approach.

A first application of this combined approach for the year 2020 yielded 2,838,080 P and S picks in the picking stage, with a total of 83,194 events after association and relocation. This is a nearly tenfold increase in event numbers compared to the IPOC catalog of Sippl et al. (2023), which contains 8,716 events for the same time interval.

In this contribution, we present results from a larger-scale application of our procedure to several years of IPOC data, and compare retrieved geometries as well as event numbers to the previously published IPOC catalog. Our findings demonstrate the potential of modern deep-learning algorithms in the creation of larger and more complete earthquake catalogs. Moving forward, our goal is to extend this preliminary catalog to span the entire 15 years of IPOC operation, facilitating in-depth analysis of regional processes.

How to cite: Najafipour, N., Puente Huerta, J. A., Sippl, C., Kasravi, J., Folesky, J., and Schurr, B.: The IPOC catalog goes deep: preliminary results, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9019, https://doi.org/10.5194/egusphere-egu24-9019, 2024.

Asperities inside the seismogenic zone of subduction megathrust are regions where specific frictional properties allow a stick-slip behavior characterized by the accumulation of slip deficit over decades to centuries and its sudden release during earthquakes. Despite its major role on the occurrence of the most devastating earthquakes and tsunamis on the planet, the physical nature of asperities and their limiting barriers is still unclear. This is partially due to an, often, ambiguous interpretation of individual geophysical proxies that are theoretically connected with the frictional structure of the megathrust at quite different time scales, ranging from 10⁰-10² yrs (seismicity patterns, geodetic locking, Vp/Vs and MT anomalies) to 10⁵-10⁷ yrs (coastal geomorphology, forearc wedge geometry and associated basal friction, magnetic and gravity anomalies). If transient phenomena, like slow slip events (SSEs) or stress shadows created by previous earthquakes, do not dominate the seismogenic behavior of the megathrust, then short- and long-term frictional proxies should coincidently illuminate the location of asperities and barriers. Moreover, this would imply that the nature of this features must be connected to the geology structure of both converging plates, with strong implications to seismic hazard assessment. A number of previous studies have explored a combination of several geophysical proxies for megathrust frictional structure, most of them along the Chilean margin. Here we expand over the work of Molina et al. (2021) and Sippl et al. (2021) by performing an integrated analysis of gravity anomalies, friction from critical wedge theory, geodetic locking and seismicity patterns for the entire 4000-km long Chilean megathrust. Particularly, we use available (micro)seismicity catalogues to compute maps of the b-value of the frequency-magnitude relationship. This parameter contributes with an independent short-term proxy for the stress state of the megathrust and we treat it as an additional continuous field into a principal component analysis (PCA) similar to Molina et al. (2021) that aims to quantify the main spatial correlation between the proxies. We will also test other techniques to measure the degree of spatial correlation, like AI-based pattern recognition methods. This integrated analysis will also consider rupture length of historical earthquakes over the last 500 yrs and slip distribution of instrumental earthquakes and SSEs. This will allow us to test contrasting hypothesis about the nature of seismic asperities and barriers along the Chilean megathrust and elsewhere.

How to cite: Tassara, A., Sippl, C., Riedel, M., Castro, C., and Carcamo, F.: Nature of asperities and barriers along the Chilean megathrust unveiled by an integrated analysis of seismicity, gravity, geodetic locking and wedge geometry, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14050, https://doi.org/10.5194/egusphere-egu24-14050, 2024.

The characterization of the spatial distribution of historical earthquake ruptures in a seismic segment plays a fundamental role in our understanding of the seismic cycle of significant earthquakes and in assessing the potential hazards associated with future events of this nature.

Due to its tectonic behavior, Chile has been impacted by megathrust earthquakes of considerable magnitude, such as the Valdivia 1960, Maule 2010, and more recently, the Illapel 2015 events. However, there are certain areas where no large earthquakes have occurred and are thus considered to be in a seismic gap. Despite experiencing some significant events, they do not manifest the required energy release properties and depth to compensate the accumulated friction. All these earthquakes, which represent varying stages of the seismic cycle, interact with different geological characteristics of the segment. This is evident in the central zone of Chile, specifically in the Valparaíso region, which has been in a seismic gap since the last major surface-rupturing earthquake of 1730.

During the Maule 2010 and Illapel 2015 earthquakes, rupture occurred only in the southern and northern segments in the mentioned area. Despite seismic activity in 1822, 1906, 1985, and 2017, and the presence of the Marga-Marga crustal fault in Viña del Mar, the energy release has not been sufficient to trigger the expected seismic sequence. It is worth noting that the fault is dangerously located in the most densely populated and frequented area of the city of Viña del Mar, presenting a threat to the surrounding population greater than what could be expected from a subduction earthquake itself.

This research aims to identify and quantify the interaction between the subducting and the Marga-Marga faults in order to assess the potential seismic activity in the area, considering that the crustal earthquakes caused by faults such as Marga-Marga are potentially more destructive than subduction earthquakes of equal magnitude. A relevant precedent is the interaction between the rupture of the Maule 2010 earthquake and the active fault segment of Pichilemu, which triggered a seismic swarm in 2011.

To achieve this, a study was conducted to characterize the slip associated with tsunamigenic events that occurred in the Central Chile segment in 1730, 1906, and 1985. The study revealed deformation patterns, indicating that the last shallow movement occurred in 1730, followed by deep patterns along the coast for subsequent events. Historical data was collected, and a stochastic modeling methodology was applied to comprehensively reconstruct the events. The Coulomb stress transmission between the Marga-Marga fault and subduction events, such as the one in 1906, was then characterized using the newly acquired information from historical deformation to identify potential activation zones of the crustal fault. Currently, efforts are underway to implement a methodology that uses computational simulation tools to visualize the impact of a coseismic event, such as the one in 1730, on the crustal fault and the surrounding region. The aim is to understand the past behavior of the region to be prepared for potential future activations.

How to cite: Álvarez, J., Calisto, I., Crempien, J., Cortés, J., Faccenna, C., and Araya, R.: Interaction between historical earthquakes in the seismic gap of central Chile and the Marga-Marga crustal Fault: The seismic potential of the Valparaiso region., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12800, https://doi.org/10.5194/egusphere-egu24-12800, 2024.

We derive a co-seismic slip model of the 2015 Mw 8.3 Illapel, Chile earthquake constrained by line-of-sight displacements from Sentinel-1 interferograms. Greens functions are calculated with 3D finite element models (FEMs). The FEMs simulate a non-uniform distribution of elastic material properties and a precise geometric configuration of the irregular topographical surface. The rupturing fault follows the curvilinear Peru-Chile Trench and Slab1.0. The optimal model that inherits heterogeneous material properties, provides a significantly better solution than that in a homogenous domain at the 95% confidence interval. The best-fit solution for the domain having a non-uniform distribution of material properties reveals a triangular slip zone. Slip is concentrated near the trench with a dip-slip up to 7.75 m, giving rise to a moment magnitude of Mw8.22 in general agreement with the seismological estimate. This methodology allows us to integrate multiple datasets of geodetic observations with seismic tomography, to achieve a better understanding of seismic ruptures within crustal heterogeneity and fault curvature.

How to cite: Tung, S. and Masterlark, T.: Revisit the coseismic slip model of the 2015 Mw 8.3 Illapel, Chile earthquake with curvilinear fault rupture, finite element model, and InSAR observations, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8443, https://doi.org/10.5194/egusphere-egu24-8443, 2024.

Characterizing the spatial distribution of ruptures from historical and recent earthquakes is key to understanding the seismic cycle of large earthquakes in subduction zones, and thus to assessing the potentialrisks associated with future earthquakes. Central Chile (35°S - 38°S) has been continuously affected by large earthquakes, such as the 2010 Maule (Mw 8.8) and the 1835 earthquakes witnessed by Robert Fitzroy (HMS Beagle captain). Here, we identify the rupture pattern and tsunami propagation of the 1751, 1835, and 2010 mega-earthquakes, events that overlapped in central Chile, by  compiling historical records and applying robust statistical tools. We used an adaptation of a logic tree methodology to generate random sources of slip distribution for each event, constrained by tsunami and deformation data. We find that the three events studied have different slip peaks. The 1751 earthquake has the largest slip with a maximum patch of ∼ 26 m, while the 2010 and 1835 earthquakes reach slips of ∼ 16 m and ∼ 10 m, respectively. Our results show that a part of the segment between 36◦S and 37◦S was consistently affected by large earthquakes, but with different slip and depth. The northern part of the segment accumulated energy for at least 300 years and was released by the 2010 earthquake. This work provides important information for identifying rupture patterns between historical and recent earthquakes, and highlights the importance of extending the time scale of earthquake slip distribution analyses to multiple cycles to describe both earthquake characteristics and their spatial relationship, and thus gain a better understanding of seismic hazard.

How to cite: Calisto, I., Cifuentes, R., San Martín, J., Alvarez, J., Ely, L., MacInnes, B., Quezada, J., and Stewart, D.: Characterizing past earthquakes through historical observations and logic tree approximation, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12727, https://doi.org/10.5194/egusphere-egu24-12727, 2024.

In recent years, several historical earthquakes have been studied in Chile to understand the seismotectonic context and anticipate the ground motion of these natural phenomena. One of the first great earthquakes documented by the Global Digital Seismographic Network (GDSN) occurred on March 3, 1985, off the coast of Valparaiso, with a moment magnitude (Mw) 8.0.

Several researchers have modeled the slip distribution at the seismic source, obtaining satisfactory results and fits, mainly at low frequencies and in far field. However, a discrepancy has been observed between the areas of maximum slip and the accelerations recorded in the near field.

In this study, the code proposed by Ruiz and Otarola (2016) was employed to stochastically generate synthetic accelerograms capable of accurately replicating the accelerations observed during near-field ground motion. This approach provides a realistic simulation of earthquake characteristics, source, path, and site.

The importance of generating synthetic accelerograms extends to critical sectors such as civil engineering, geophysics, construction, and urban planning. These simulations play a critical role in understanding ground behavior, predicting large seismic movements, and improving the development of earthquake-resistant structures. Furthermore, in the fields of construction and urban planning, synthetic accelerograms are essential for assessing the vulnerability of specific areas, diversifying applications in industry, and facilitating a more resilient design approach for future seismic events.

The results obtained by generating synthetic accelerograms can replicate the spectral and temporal shape, in agreement with the records provided by the National Seismological Center (CSN) network. Stochastic simulations have been run both in rocky environments and in areas with site effects.

How to cite: Torres, R. and Ruiz, S.: Stochastic Strong-Motion Simulation of Valparaiso 1985 Mw 8.0 Chile Earthquake, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13893, https://doi.org/10.5194/egusphere-egu24-13893, 2024.

The Andean Margin hosts alternating regions of “flat” and “normal” subduction, which includes the Pampean flat slab that extends from central Chile to Argentina. The discovery of an unusual travel time anomaly beneath the high Andes above the flat slab motivated a study to investigate the lithosphere in this region. Leveraging extensive archived seismic data from both Chile and Argentina, we performed a large-scale joint inversion of P and S body wave arrival times from earthquakes, and surface wave dispersion measurements from earthquakes and ambient noise. We created 3D Vp, Vs and Vp/Vs models using at least an order of magnitude more data than previous studies with about an 80% reduction in grid spacing. Our models corroborate results from previous studies: (1) a high velocity, high Vp/Vs region associated with a cool, slightly hydrated and depleted mantle above the flat slab, and (2) a low velocity structure beneath the high Andes interpreted as an overthickened crustal root, with our results showing that the root extends to just above the flat slab. Curiously, our models also reveal two low velocity zones within and below the flat slab seismic zone that have not been previously reported. Notably, the decrease in velocity is more pronounced in Vp than Vs. We postulate that the eastern low velocity anomaly is likely due to hot asthenosphere heating the slab, although no melting is occurring as the Vs is not significantly reduced. The western low velocity anomaly, which spatially correlates with the Juan Fernandez Ridge (JFR), we postulate is either due to the presence of supercritical fluids trapped within the JFR or an increase in silica content possibly linked to petit spot volcanism.

How to cite: Maharaj, A., Roecker, S., Comte, D., Saez, M., Trad, S., Ortiz, G., and Fernadez, M.: Joint Tomographic Inversion of the Pampean Flat Slab: Insights from Extensive Archived Seismic Data, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9376, https://doi.org/10.5194/egusphere-egu24-9376, 2024.

In subduction zones, intermediate-depth earthquakes typically occur in two discrete layers, delineating an upper and a lower seismicity plane with a separating aseismic or minimally seismic region, a phenomenon named Double Seismic Zone (DSZ). However, the seismicity pattern in Northern Chile features two parallel planes of seismic activity only in the shallower section of the slab. Around depths of 80–90 km, the seismicity undergoes a transition to a significantly seismogenic zone approximately 25–30 km thick, effectively connecting the initial seismicity planes. This variation presents a distinct form of intraslab seismicity that deviates from the traditional DSZ structure and prompts further investigation into its underlying mechanisms and implications for regional seismic hazard assessment. Insights derived from this region's seismicity could provide pivotal constraints and enhance our understanding of the complex interplay between geological processes, mineral transformations, and fluid migrations in shaping subduction zone seismicity.

The attenuation of seismic waves in a rock volume is a property that is highly sensitive to the presence and concentration of fluids as well as spatial variations of temperature. As intermediate-depth seismicity is thought to originate from dehydration processes in the downgoing slab, along-strike or along-dip changes in slab seismicity should have a signature in seismic attenuation of the slab as well as the overlying mantle wedge. We hence aim at better understanding the aforementioned seismicity configuration in Northern Chile by acquiring a 3D image of its attenuation signature. 

The primary dataset for our analysis comes from the seismic stations of the Integrated Plate boundary Observatory Chile (IPOC) network in Northern Chile's forearc, augmented by additional data from different temporary deployments. Using the extensive seismicity catalog of Sippl et al. (2023), we have about 180,000 events at over 50 seismic stations at our disposal from the period 2007 to 2021; we select only the high quality traces for the analysis. The rays are traced in a 3D velocity model. We invert the spectral ratios obtained with the coda normalization method to obtain total-Q values. We present images of the 3D attenuation structure of the Northern Chile Forearc between 21ºS and 23ºS, which are obtained with measurements of the coda normalization method using the Multi-Resolution Attenuation Tomography algorithm (Sketsiou et al., 2021).

How to cite: Castro-Melgar, I. and Sippl, C.: Towards 3D attenuation tomography of the Northern Chile forearc, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7958, https://doi.org/10.5194/egusphere-egu24-7958, 2024.

Intermediate-depth seismicity in Northern Chile shows a pattern which is distinct from a conventional double seismic zone (DSZ) setting. While two distinct seismicity planes are present in the updip part of the slab, there is a sharp transition to a highly seismogenic cluster of 25–30 km thickness at ∼80-90 km depth, extinguishing the gap between the two seismicity planes. As seismic velocities can be used to constrain mineralogy and fluid content, characterizing seismic wavespeeds of this subduction zone segment using local earthquake tomography can provide important constraints on the mineralogical processes that produce the seismicity pattern seen here.

We used the catalog of Sippl et al. (2018), which contains arrival time data from permanent stations of IPOC (Integrated Plate Boundary Observatory Chile), complemented by several temporary deployments spanning shorter time sequences. The catalog contains more than 100,000 earthquakes and 1,200,404 P- and 688,904 S-phase picks for the years 2007 to 2014. In order to use the best available picks for tomography, we limit our analysis to events that have more than 14 P-arrivals as well as more than 7 S-arrivals, leading to a total of 8883 events with 213,908 P- and 99466 S-arrivals.

Parallelly, we also attempt to obtain an estimate of the possible mineral compositions at the depths and P-T conditions relevant to our study in the DSZ setting. For this, we assume a simple model where the upper plane of the DSZ is considered to be evolving from MORB-like composition and the lower plane of the DSZ from depleted-mantle composition. These global average compositions are then fed into Perple X (Connolly & Kerrick, 1987) as starting compositions and pseudosections of possible mineral assemblages are constructed for P-T conditions significant to this study. We calculate the theoretical Vp and Vp/Vs values for those P-T conditions using the same software.

We present 3D models of P- and S-wavespeeds from the Northern Chile forearc between about 20.4° S and 22.5° S, as well as images of ray coverage, relocated seismicity, and synthetic resolution tests. Tomography models for different choices of grid spacing and damping-smoothing parameters are compiled and compared to derive the optimal settings for the inversion. The seismic velocity distribution obtained through tomography is compared with the aforementioned theoretical wavespeeds to narrow down the range of possible reactions occurring at depth.

How to cite: Hassan, N., Henry, S., and Sippl, C.: Combining local earthquake tomography and petrological models to constrain wavespeeds in the subducting Nazca Plate, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5730, https://doi.org/10.5194/egusphere-egu24-5730, 2024.

We collected earthquake waveform data recorded by permanent seismic stations in northern Chile from 2014 to 2019 to construct a new earthquake catalog, and integrated them with the previous catalog data. In total, the new catalog consisted of 536342 P and 453920 S arrival times from 52165 earthquakes and 245 stations. We resolved Vp, Vs, and Vp/Vs models and seismic locations for northern Chile by using a new version of double-difference seismic tomography method based on Vp/Vs model consistency constraint (Guo et al., 2018). The new velocity models provide a refined structure of the subducting slab down to 350 km.

The earthquake relocations reveal a distinct double seismic zone in northern Chile, but the gap between the two seismic planes disappears at a depth of approximately 100 km and replaced by a concentration of seismic cluster. Under this intermediate-depth seismic cluster, several isolated small seismic clusters remain. The tomography results indicate a strong correlation between seismicity distribution and high-velocity anomalies. The subducting Nazca Plate presents stripe-like high-velocity anomalies with clear segmentations, potentially related to the weakening at the outer-rise of the trench. Furthermore, our Vp/Vs model indicates that the upper seismic plane exhibits high Vp/Vs anomalies, which may indicate the presence of fluids released from dehydration reactions of various hydrous minerals. In contrast, lower seismic plane and deep seismic clusters are associated with low Vp/Vs anomalies, which could be related to supercritical fluids. Additionally, the enhanced seismicity and velocity anomalies in the region of 21-22ºS along the strike suggest a potential influence of the subduction of the Iquique Ridge of the Nazca Plate.

How to cite: Chen, Z., Zhang, H., and Gao, L.: Joint inversion for Vp, Vs, and Vp/Vs of subduction zone in northern Chile, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8641, https://doi.org/10.5194/egusphere-egu24-8641, 2024.

Between 2020 and 2022 the ANILLO+DEEPtrigger (Y6+XZ) Seismic Network, comprising 108 seismic stations, operated for eighteen months in Northern-Central Chile (24.5°S - 29°S). Employing Deep Learning (EQTransformer, Mousavi et al., (2022)) and Phase Association (GaMMA, Zhu et al. (2021)) algorithms, we identified over 30,000 seismic events in an area with a notable absence of moderate-to-large events in the past century, since the 1922 M8.5 Atacama earthquake.  
From the initial catalog, we selected a well-distributed subcatalog of 1000 earthquakes, consisting of 26,570 P- and 22,109 S-wave arrival times, by selecting for events with an optimal spatial distribution, small residuals, and abundant P- and S-arrivals. These selected events served as input for VELEST (Kissling et al., 1994) to compute a new 1-D velocity model representative of this region by minimizing the subset residuals. To reduce both residuals and location errors associated with the seismicity, we relocated the entire catalog using staggered tomographic inversions based on SIMUL2000 (Thurber & Eberhart-Phillips, 1999), simultaneously inverting for seismic velocity models and hypocentral parameters within the iterative damped least squares method. Following the proposed method, we gradually increased model complexity, transitioning from 2-D Vp and Vp/Vs to ultimately a 3-D fine Vp and Vp/Vs solution with low node separation.
Next, synthetic resolution tests were conducted to assess the reliability of the spatial limits and boundaries within the solutions. In this context, distinctive patterns were identified for each profile of the three-dimensional model, revealing enhanced horizontal and vertical resolution in the central region beneath the network. Conversely, a decline in resolution was noted at the peripheries, primarily attributable to reduced station coverage causing poorer seismic event relocations.
Our results reveal both regional and local patterns. We observed a mantle wedge with vertical thicknesses ranging from ~35km in the southernmost profiles less than 25 km in the northern region, consistent with previous seismic tomography observations in northern Chile (Pastén-Araya et al., 2021). The Vp/Vs ratio and Vp values allow us to discern the distribution of the hydrated slab, which, spatial correlated with seismicity, provides evidence of irregular dehydratation processes along both dip and strike directions.
Relocated seismicity exhibits some noteworthy features. Shallower crustal sesmicity is predominantly related to high rates of mining activity. In the subduction areas, the most prominent cluster is located at depths of 20-50 km, delineating the seismogenic zone. At greater depths, double and even triple seismic bands add structural complexities to the observations.
From 26.5°S to 29.5°S, between 20 km and ~75 km depth, seismicity predominantly aligns with the interplate contact defined by SLAB2 (Hayes et al., 2018). In contrast, northward from 26.5°S, our deepest seismicity, situated between 75 and 125 km depth, diverges from SLAB2, depicting a steeper dip angle.
Lastly, we recommend integrating OBS and back-arc stations, whose data would improve off-shore and back-arc resolution, contributing to a more comprehensive understanding of seismotectonic environments. Non-supervised Deep Learning results can provide exceptional databases for tomographic studies, yielding residuals similar to human-picked databases but within shorter timeframes.

How to cite: Hernández-Soto, N., Miller, M., Moreno, M., Lange, D., Socquet, A., Sippl, C., and González-Vidal, D.: From regional to local structures imaged by seismic tomography at the Atacama seismic gap, Central-Northern Chile (24.5-29°S), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-20700, https://doi.org/10.5194/egusphere-egu24-20700, 2024.