Seismotectonics of the southern subduction Chilean margin revealed by recent aftershock sequences

Subduction margins, as in the case of south-central Chile, are active seismotectonic environments and locus of the world largest earthquakes. In this thesis, two segments of the south-central Chilean subduction margin are studied: (A) the southernmost portion, at the termination of the Nazca-South America convergence (~46ºS), and (B) the segment located between 34º-38ºS, where the Mw 8.8 Maule Earthquake took place in 2010. Analysis of data from a local seismic network deployed in 2004-2005 in area A, indicates low levels of background seismicity with magnitudes ranging 0-3.4 Ml. The seismicity corresponds to shallow crustal events, mostly occurring within the upper 10 km. A third of the seismicity is associated to volcanic activity present in the area, while scarce seismicity is associated with a large strike-slip fault, the Liquiñe-Ofqui Fault System (LOFS), that intersects the region along the arc in a N-S-trend. In 2007, this region was affected by a seismic sequence with a peak of activity associated with a Mw 6.2 earthquake in April that year. A local seismic network was deployed after this main event in order to study its sequence of aftershocks, which provided a unique opportunity to characterise seismotectonically this area that usually lacks intermediate magnitude seismicity, including the calculation of a new local velocity model, accurate aftershock locations and computation of focal mechanisms. The results show P-wave velocities of ~5 km/s for the upper 5 km in accordance with the geology of the area, and low S-wave velocities for the upper 3 km of crust due to rock fracturing and the presence of fluids. An average Vp/Vs ratio of 1.76 was calculated for the region. The alignment of most of the aftershocks within the LOFS plus obtained focal mechanisms, indicate that this sequence had tectonic origin related to the re-activation of the LOFS. Further, a maximum seismogenic depth of about 15 km was determined for the entire region. Regarding area B, affected by a large megathrust earthquake in 2010, the study of moment tensor solutions for the sequence of aftershocks provided new insight into the distribution of postseismic activity relative to co-seismic slip and the release of seismic afterslip. Thrust aftershocks dominate the postseismic activity, but also normal faulting was detected in the outer-rise area and in the overriding plate near the coastline. The largest seismically released afterslip is located between the two main patches of co-seismic slip. Large aftershocks (M>4) occur along the megathrust interface, in zones of intermediate co-seismic slip associated to stress introduced on dislocation tips with high co-seismic slip contrast. On the other hand, smaller events (

Seismic imaging of the Northern Andean subduction zone from teleseismic tomography: a torn and fragmented Nazca slab

The Nazca-South America subduction zone in Ecuador is characterized by a complicated along-strike geometry as the slab transitions from flat slab subduction in the south, with the Peruvian flat slab, to what has been characterized as ‘normal’ dipping subduction beneath central Ecuador. Plate convergence additionally changes south to north as the trench takes on a convex shape. Highly heterogeneous bathymetry at the trench, including the aseismic oceanic Carnegie Ridge (CR), and sparse intermediate-depth seismicity has led many to speculate about the behaviour of the downgoing plate at depth. In this study, we present a finite-frequency teleseismic P-wave tomography model of the northern Andes beneath Ecuador and Colombia from 90 to 1200 km depth. Our model builds on prior tomography models in South America by adding relative traveltime residuals recorded at stations in Ecuador. The complete data set is comprised of 114 096 relative traveltime residuals from 1133 stations across South America, with the added data serving to refine the morphology of the Nazca slab in the mantle beneath the northern Andes. Our tomography model shows a Nazca slab with a fragmented along-strike geometry and the first teleseismic images of several proposed slab tears in this region. At the northern edge of the Peruvian flat slab in southern Ecuador, we image a shallow tear at 95–200 km depth that appears to connect mantle flow from beneath the flat slab to the Ecuadorian Arc. Beneath central Ecuador at the latitudes of the CR, the Nazca slab is continuous into the lower mantle. Beneath southern Colombia, the Malpelo Tear breaks the Nazca slab below ∼200 km depth

The 2012–2013 Montes Claros earthquake series in the São Francisco\ud Craton, Brazil: new evidence for non-uniform intraplate stresses\ud in mid-plate South America

On 2012 May 19, an mb = 4 earthquake shook the town of Montes Claros, Brazil in the\ud middle of the S˜ao Francisco Craton. Because of the scarce seismicity in the area, an event\ud like this could provide valuable information to characterize the governing seismotectonics and\ud stress field for the region. Here, we present the results of more than 1 yr of local seismic\ud monitoring after the main shock. We found that the seismicity originated at approximately\ud 1-km depth in an NNW-oriented blind reverse fault, dipping to the E. The magnitude of\ud the main shock was 4mb, with aftershocks reaching up to 3.6mb. Focal mechanisms from\ud first motion polarities and waveform moment tensor inversions indicate a reverse faulting in\ud agreement with the orientation of the aftershock locations. In addition, we derived a new 1-D\ud local velocity model using a simultaneous inversion of hypocentres and velocity layers. The\ud results indicate P-wave velocities of 4.5 km s−1 for the upper layer of carbonate rocks and 5.23\ud and 5.69 km s−1 for the lower fractured and compact crystalline basement layers, respectively.\ud Higher Vp/Vs ratios were obtained for the upper two layers compared to the lowermost layer,\ud possibly indicating presence of rock fracturing and percolated water. The calculated stress drop\ud for the main event is 0.33 MPa, which is a relatively low value for an intraplate earthquake\ud but still within the observed range. The inversion of the main shock focal mechanism and\ud previously published focal mechanisms suggests a compressional stress regime in the central\ud part of the S˜ao Francisco Craton, which is different from the strike-slip regime in the southern\ud part, although both have an EW-oriented σ1. On the other hand, focal mechanisms of events\ud located to the west of the craton indicate an NW–SE oriented σ1 for central Brazil. This\ud variability highlights the importance of local sources of stresses (e.g. flexural stresses) in\ud mid-plate South America, unlike other mid-plate areas of the world, such as central and east\ud North America, where a more uniform stress field is observedPetrobras - BRASIS projec

The tailings dam failure of 5 November 2015 in SE Brazil and its preceding seismic sequence

The collapse of a mine tailings dam and subsequent flood in SE Brazil on 5 November 2015 was preceded by a small-magnitude seismic sequence. In this report, we explore the spatiotemporal associations between the seismic events and the accident and discuss their possible connection. We also analyze the signals generated by the turbulent mudflow, as recorded by the Brazilian Seismographic Network (RSBR). In light of our observations, we propose as possible contributing factor for the dam collapse either ground shaking and/or soil liquefaction triggered by the earthquakes. The possibility of such a small-magnitude earthquake contributing to the collapse of a tailings dam raises important concerns regarding safety and related legislation of dams in Brazil and the world. ©2016. American Geophysical Union.H.A.D. and M.A. acknowledge support from Sao Paulo Research Foundation FAPESP grant 2014/09455-3 and CNPq grant 30.6547/2013-9.Peer reviewe

Triggered crustal earthquake swarm across subduction segment boundary after the 2016 Pedernales, Ecuador megathrust earthquake

Megathrust ruptures and the ensuing postseismic deformation cause stress changes that may induce seismicity on upper plate crustal faults far from the coseismic rupture area. In this study, we analyze seismic swarms that occurred in the north Ecuador area of Esmeraldas, beginning two months after the 2016 M$_{w}$ 7.8 Pedernales, Ecuador megathrust earthquake. The Esmeraldas region is 70 km from the Pedernales rupture area in a separate segment of the subduction zone. We characterize the Esmeraldas sequence, relocating the events using manual arrival time picks and a local a-priori 3D velocity model. The earthquake locations from the Esmeraldas sequence outline an upper plate fault or shear zone. The sequence contains one major swarm and several smaller swarms. Moment tensor solutions of several events include normal and strike-slip motion and non-double-couple components. During the main swarm, earthquake hypocenters increase in distance from the first event over time, at a rate of a few hundred meters per day, consistent with fluid diffusion. Events with similar waveforms occur within the sequence, and a transient is seen in time series of nearby GPS stations concurrent with the seismicity. The events with similar waveforms and the transient in GPS time series suggest that slow aseismic slip took place along a crustal normal fault during the sequence. Coulomb stress calculations show a positive Coulomb stress change in the Esmeraldas region, consistent with seismicity being triggered by the Pedernales mainshock and large aftershocks. The characteristics of the seismicity indicate that postseismic deformation involving fluid flow and slow slip activated upper plate faults in the Esmeraldas area. These findings suggest the need for further investigation into the seismic hazard potential of shallow upper plate faults and the potential for megathrust earthquakes to trigger slow-slip and shallow seismicity across separate segments of subduction zones

Structural Control on Megathrust Rupture and Slip Behavior: Insights From the 2016 Mw 7.8 Pedernales Ecuador Earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 M$_{w}$ 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906M$_{w}$ 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a M$_{w}$ 7.7 earthquake

3D Local Earthquake Tomography of the Ecuadorian Margin in the Source Area of the 2016 Mw 7.8 Pedernales Earthquake

Based on manually analyzed waveforms recorded by the permanent Ecuadorian network and our large aftershock deployment installed after the Pedernales earthquake, we derive three-dimensional Vp and Vp/Vs structures and earthquake locations for central coastal Ecuador using local earthquake tomography. Images highlight the features in the subducting and overriding plates down to 35 km depth. Vp anomalies (∼4.5–7.5 km/s) show the roughness of the incoming oceanic crust (OC). Vp/Vs varies from ∼1.75 to ∼1.94, averaging a value of 1.82 consistent with terranes of oceanic nature. We identify a low Vp (∼5.5 km/s) region extending along strike, in the marine forearc. To the North, we relate this low Vp and Vp/Vs (1.85) which we interpret as deeply fractured, probably hydrated OC caused by the CR being subducted. These features play an important role in controlling the seismic behavior of the margin. While subducted seamounts might contribute to the nucleation of intermediate megathrust earthquakes in the northern segment, the CR seems to be the main feature controlling the seismicity in the region by promoting creeping and slow slip events offshore that can be linked to the updip limit of large megathrust earthquakes in the northern segment and the absence of them in the southern region over the instrumental period

1D-velocity structure and seismotectonics of the Ecuadorian margin inferred from the 2016 Mw7.8 Pedernales aftershock sequence

International audienceOn April 16th 2016 a Mw 7.8 earthquake ruptured the central coastal segment of the Ecuadorian subduction zone. Shortly after the earthquake, the Instituto Geofisico de la Escuela Politecnica Nacional of Ecuador, together with several international institutions deployed a dense, temporary seismic network to accurately categorize the post-seismic aftershock sequence. Instrumentation included short-period and broadband sensors, along with Ocean Bottom Seismometers. This deployment complemented the permanent Ecuadorian seismic network and recorded the developing aftershock sequence for a period of one year following the main-shock. A subset of 345 events with ML > 3.5, were manually picked in the period of May to August 2016, providing highly accurate P- and S-onset times. From this catalogue, a high-quality dataset of 227 events, with an azimuthal gap <200°, are simultaneously inverted for, obtaining the minimum 1D velocity model for the rupture region, along with hypocentral locations and station corrections. We observe an average Vp/Vs of 1.82 throughout the study region, with relatively higher Vp/Vs values of 1.95 and 2.18 observed for the shallowest layers down to 7.5 km. The high relative Vp/Vs ratio (1.93) of the deeper section, between 30 km and 40 km, is attributed to dehydration and serpentinization processes. For the relocated seismicity distribution, clusters of events align perpendicular to the trench, and crustal seismicity is also evidenced, along with earthquakes located close to the trench axis. We also compute Regional Moment Tensors to analyze the different sources of seismicity after the mainshock. Aside from thrust events related to the subduction process, normal and strike-slip mechanisms are detected. We suggest that the presence of subducting seamounts coming from the Carnegie Ridge act as erosional agents, helping to create a scenario which promotes locking and allows seismicity to extend up to the trench, along zones of weakness activated after large earthquakes