skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Structure of the Ecuadorian forearc from the joint inversion of receiver functions and ambient noise surface waves
SUMMARY The Ecuadorian forearc is a complex region of accreted terranes with a history of large megathrust earthquakes. Most recently, a Mw 7.8 megathrust earthquake ruptured the plate boundary offshore of Pedernales, Ecuador on 16 April 2016. Following this event, an international collaboration arranged by the Instituto Geofisico at the Escuela Politécnica Nacional mobilized a rapid deployment of 65 seismic instruments along the Ecuadorian forearc. We combine this new seismic data set with 14 permanent stations from the Ecuadorian national network to better understand how variations in crustal structure relate to regional seismic hazards along the margin. Here, we present receiver function adaptive common conversion point stacks and a shear velocity model derived from the joint inversion of receiver functions and surface wave dispersion data obtained through ambient noise cross-correlations for the upper 50 km of the forearc. Beneath the forearc crust, we observe an eastward dipping slow velocity anomaly we interpret as subducting oceanic crust, which shallows near the projected centre of the subducting Carnegie Ridge. We also observe a strong shallow positive conversion in 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 also observe significant north–south changes in shear wave velocity. The velocity changes indicate variations in the accreted terranes and may indicate an increased amount of hydration beneath the Manabí Basin. This change in structure also correlates geographically with the southern rupture limit of multiple high magnitude megathrust earthquakes. The earthquake record along the Ecuadorian trench shows that no event with a Mw >7.4 has ruptured south of ∼0.5°S in southern Ecuador or northern Peru. Our observations, along with previous studies, suggest that variations in the forearc crustal structure and subducting oceanic crust may influance the occurrence and spatial distribution of high magnitude seismicity in the region.  more » « less
Award ID(s):
1723042 1723065
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Page Range / eLocation ID:
1671 to 1685
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY The 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. 
    more » « less
  2. null (Ed.)
    Recent GPS studies show that the Indo-Burma subduction system is locked with the implication of a potential large-magnitude earthquake. To inform better seismic hazard models in the region, we need an improved understanding of the crustal structure and the dynamics of the Indo-Burma subduction system. The Bangladesh-India-Myanmar (BIMA) tripartite project deployed 60 broadband seismometers across the subduction system and have been continuously recording data for ~2 years. In this study, we computed receiver functions from 30 high-quality earthquakes (M≥5.9) with epicentral distances between 30º and 90º recorded by the array. The algorithm utilized ensures the uniqueness of the seismic model and provides an uncertainty estimate of every converted wave amplitude. We stacked all the receiver functions produced at each station along the entire transect to generate a cross-sectional model of the average crustal structure. The level of detail in the image is improved by computing higher frequency receiver functions up to 4 Hz. The results represent some of the strongest constraints on crustal structure across the subduction system. Beneath the Neogene accretionary prism's outer belt, we observe a primary conversion associated with the Ganges Brahmaputra Delta that ranges in depth from ~10 km near the deformation front up to ~12 km at the eastern boundary. From the eastern end of the Neogene accretionary prism to the Sagaing Fault, we image the Indian subducting slab and the Central Myanmar basin. The depth-extent of seismicity associated with the Wadati-Benioff zone is consistent with the locations of primary conversions from the subducting plate. We further verify the converted phases of the slab by analyzing azimuthal moveout variations. The Central Myanmar basin is roughly bowl-shaped in cross-section with a maximum thickness of ~15 km about halfway between the Kabaw and Sagaing faults. The average crustal thickness beneath the Ganges-Brahmaputra delta is ~20 km, most likely representing a transitional crust formed from thinning of the continental crust intruded and underplated by igneous rocks. In contrast, the average thickness of the continental crust beneath the Central Myanmar basin is ~40 km. Our results provide a baseline model for future geophysical investigations of the Indo-Burma subduction zone. 
    more » « less
  3. Subduction of the very young (<15 Myr old) oceanic lithosphere of the Nazca plate in central to southern Colombia is observationally related to an unusually high and unusually variable amount of intermediate (>50 km) depth seismicity. From 2010 through 2019 89% of central and southern Colombia’s 11,466 intermediate depth events occurred between 3.5°N and 5.5°N, highlighting these unusual characteristics of the young slab. In addition, morphologic complexity and possible tears characterize the Nazca slab in Colombia and complicate mantle flow in the region. Prior SKS-phase shear-wave splitting results indicate sub-slab anisotropy is dominated by plate motion parallel-to-subparallel orientations in the region, suggesting the young slab has entrained a relatively thick portion of the sub-slab mantle. These observations suggest the subduction of young lithosphere has significant effects on both the overlying and underlying asthenosphere in the Colombia subduction zone. Here we use more than 10 years of data to calculate receiver functions for the Red Sismológica Nacional de Colombia’s network of broadband seismometers. These receiver functions allow us to tie these prior observations of the Colombia subduction zone to distinct, structural features of the slab. We find that the region of high seismicity corresponds to a low seismic velocity feature along the top of the subducting plate between 3.5°N and 5.5°N that is not present to the south. Moderately elevated P-wave velocity to S-wave velocity ratios are also observed within the slab in the north. This feature likely represents hydrated slab mantle and/or uneclogitized oceanic crust extending to a deeper depth in the north of the region which may provide fluids to drive slab seismicity. We further find evidence for a thick layer of material along the slab’s lithosphere-asthenosphere boundary characterized by spatially variable anisotropy. This feature likely represents entrained asthenosphere at the base of the plate sheared by both the overlying plate and complex flow related to proposed slab tears just north and south of the study region. These observations highlight how structural observations provide key contextual constraints on short-term (seismogenic) and long-term (anisotropic fabric) dynamic processes in the Colombia subduction zone. Plain-language Summary The Nazca oceanic plate is very young (<15 million years old) where it is pulled or subducted beneath the South America plate in central and southern Colombia. Earthquakes occurring in the subducted Nazca plate at depths greater than 50 km are nearly 9x more common in central Colombia than in southern Colombia. The subducted Nazca plate also has a complex shape in this region and may have been torn both in northern Colombia and to the south near the Colombia-Ecuador border. The slow flow of mantle rock beneath the subducted plate is believed to be affected by this and earlier studies have inferred this flow is mostly in the same direction as the subducting plate's motion. We have used 10+ years of data to calculate receiver functions, which can detect changes in the velocity of seismic waves at the top and bottom of the subducted plate to investigate these features. We found that the Nazca plate is either hydrated or has rocks with lower seismic velocities at its top in the central part of Colombia where earthquakes are common. We also find that a thick layer of mantle rock at the base of the subducted plate has been sheared. 
    more » « less
  4. Abstract

    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 Mw7.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 1906 Mw8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw7.7 earthquake.

    more » « less
  5. Abstract

    The south central Chilean margin regularly produces many of the world's largest earthquakes and tsunami, including the 2010Mw8.8 Maule and 1960Mw9.5 Valdivia events. In 2017, we acquired seismic reflection data along ~1,000 km of the margin using the R/VLangseth's 15 km long receiver array and 108.2 l (6,600 in3) seismic source to image structures associated with these ruptures. We focus on the Valdivia segment with the largest coseismic slip (~40 m). The outer 40 km of the forearc is an accretionary wedge constructed primarily of stacked sedimentary packages with irregular lengths and thicknesses and little along‐strike continuity. Forearc structures indicate that the accretionary wedge grows primarily through basal accretion of the downgoing trench fill. The décollement propagates along a weak boundary near the top of the trench fill but occasionally branches downward into the underthrust sediment along bedding horizons, peeling off slices that are underplated to the forearc. The shallow décollement level and the rarity of underplating events allow most of the trench sediment to subduct. As a result, only ~30% of the incoming sediment has been accreted since the Early Pliocene. This implies that, on average, ~1 km of sediment must subduct beyond the outer forearc, an inference that is supported by our seismic images. We propose that the thickness and great downdip and along‐strike extent of the underthrust layer, which separates the megathrust from the underlying roughness of the igneous ocean crust, ensures a smooth broad zone of strong coupling that generates the world's largest earthquakes and tsunami.

    more » « less