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

Postseismic deformation: the Esmeraldas, Ecuador Seismic sequence Following the 2016 Pedernales Megathrust Earthquake

International audienc

Structure of the Ecuadorian forearc from the joint inversion of receiver functions and ambient noise surface waves

International audienceThe Ecuadorian forearc is a complex region of accreted terranes with a history of largemegathrust earthquakes. Most recently, aMw7.8 megathrust earthquake ruptured the plateboundary offshore of Pedernales, Ecuador on 16 April 2016. Following this event, an inter-national collaboration arranged by the Instituto Geofisico at the Escuela Polit ́ecnica Nacionalmobilized a rapid deployment of 65 seismic instruments along the Ecuadorian forearc. Wecombine this new seismic data set with 14 permanent stations from the Ecuadorian nationalnetwork to better understand how variations in crustal structure relate to regional seismichazards along the margin. Here, we present receiver function adaptive common conversionpoint stacks and a shear velocity model derived from the joint inversion of receiver functionsand surface wave dispersion data obtained through ambient noise cross-correlations for theupper 50 km of the forearc. Beneath the forearc crust, we observe an eastward dipping slowvelocity anomaly we interpret as subducting oceanic crust, which shallows near the projectedcentre of the subducting Carnegie Ridge. We also observe a strong shallow positive conversionin the Ecuadorian forearc near the Borbon Basin indicating a major discontinuity at a depth of∼7 km. This conversion is not ubiquitous and may be the top of the accreted terranes. We alsoobserve significant north–south changes in shear wave velocity. The velocity changes indicatevariations in the accreted terranes and may indicate an increased amount of hydration beneaththe Manab ́ı Basin. This change in structure also correlates geographically with the southernrupture limit of multiple high magnitude megathrust earthquakes. The earthquake record alongthe Ecuadorian trench shows that no event with aMw>7.4 has ruptured south of∼0.5◦Sinsouthern Ecuador or northern Peru. Our observations, along with previous studies, suggestthat variations in the forearc crustal structure and subducting oceanic crust may influance theoccurrence and spatial distribution of high magnitude seismicity in the region

Structural Controls on Postseismic Deformation Following the Mw 7.8 Pedernales, Ecuador Megathrust Earthquake: Insights from Joint Tomographic Inversion and Aftershock Relocation

International audienc

Upper-plate structure in Ecuador coincident with the subduction of the Carnegie Ridge and the southern extent of large mega-thrust earthquakes

International audienceThe Ecuadorian convergent margin has experienced many large mega-thrust earthquakes in the past century, beginning with a 1906 event that propagated along as much as 500 km of the plate interface. Many subsections of the 1906 rupture area have subsequently produced Mw ≥ 7.7 events, culminating in the 16 April 2016, Mw 7.8 Pedernales earthquake. Interestingly, no large historic events Mw ≥ 7.7 appear to have propagated southward of ∼1◦S, which coincides with the subduction of the Carnegie Ridge. We combine data from temporary seismic stations deployed following the Pedernales earthquake with data recorded by the permanent stations of the Ecuadorian national seismic network to discern the velocity structure of the Ecuadorian forearc and Cordillera using ambient noise tomography. Ambient noise tomography extracts Vsv information from the ambient noise wavefield and provides detailed constraints on velocity structures in the crust and upper mantle. In the upper 10 km of the Ecuadorian forearc, we see evidence of the deepest portions of the sedimentary basins in the region, the Progreso and Manabı´ basins. At depths below 30 km, we observe a sharp delineation between accreted fast forearc terranes and the thick crust of the Ecuadorian Andes. At depths ∼20 km, we see a strong fast velocity anomaly that coincides with the subducting Carnegie Ridge as well as the southern boundary of large mega-thrust earthquakes. Our observations raise the possibility that upper-plate structure, in addition to the subducting Carnegie Ridge, plays a role in the large event segmentation seen along the Ecuadorian margin

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

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

International audienceMegathrust 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 Mw7.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

Repeating earthquakes cluster at the edges of the afterslip in the aftermath of the 16th April 2016 M7.8 Pedernales earthquake in Ecuador

International audienceThe study of repeating earthquakes can provide us with a deeper understanding of the interaction between seismic and aseismic slip. Repeating earthquakes are the result of repeated loading and failure of single fault asperities, likely driven by surrounding aseismic slip. As such, they are especially useful in the study of transient or stable slow slip, such as afterslip, and can provide information on slip properties and fault mechanics, with possible applications to earthquake triggering, loading rates and earthquake forecasting.Off the coast of Ecuador, seismicity during both the inter-seismic and post-seismic periods has been previously linked to aseismic slip. We analyze one year of aftershocks following the 16th April 2016 Mw 7.8 Pedernales earthquake to find repeating families, using data recorded by permanent and temporary seismological stations. We calculate waveform cross-correlation coefficients (CC) on all available catalogue events, which we use to sort events into preliminary families, using a minimum CC of 0.95. These events were then stacked and used to perform template-matching on the continuous data. In total, 376 earthquakes were classified into 62 families of 4 to 15 earthquakes, including 8 from the one-year period before the mainshock. We later relocated these earthquakes using a double-difference method, which confirmed that most of them did have overlapping sources. Finally, repeaters’ magnitudes and corner frequencies were determined within clusters using the spectral ratio method. With these we were able to calculate stress drop variations in time and space.Repeating earthquakes seem to concentrate largely around the areas of largest afterslip release, where afterslip gradient is the highest. We also find an increase in the recurrence time of repeating events with time after the mainshock, over the first year of the postseismic period, which highlights a possible timeframe for the afterslip’s deceleration. Our results suggest that while most repeating aftershocks are linked to afterslip release, the afterslip gradient may play a bigger role in determining their location than previously thought