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
Upper-plate structure in Ecuador coincident with the subduction of the Carnegie Ridge and the southern extent of large mega-thrust earthquakes
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
- NSF-PAR ID:
- 10184035
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 220
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 1965 to 1977
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Through the Alaska Transportable Array deployment of over 200 stations, we create a 3‐D tomographic model of Alaska with sensitivity ranging from the near surface (<1 km) into the upper mantle (~140 km). We perform a Markov chain Monte Carlo joint inversion of Rayleigh wave ellipticity and phase velocities, from both ambient noise and earthquake measurements, along with receiver functions to create a shear wave velocity model. We also use a follow‐up phase velocity inversion to resolve interstation structure. By comparing our results to previous tomography, geology, and geophysical studies we are able to validate our findings and connect localized near‐surface studies with deeper, regional models. Specifically, we are able to resolve shallow basins, including the Copper River, Cook Inlet, Yukon Flats, Nenana, and a variety of other shallower basins. Additionally, we gain insight on the interaction between the upper mantle wedge, asthenosphere, and active and nonactive volcanism along the Aleutians and Denali volcanic gap, respectively. We observe thicker crust beneath the Brooks Range and south of the Denali fault within the Wrangellia Composite Terrane and thinner crust in the Yukon Composite Terrane in interior Alaska. We also gain new perspective on the Wrangell Volcanic Field and its interaction between surrounding asthenosphere and the Yakutat Terrane.
-
Palin, Richard (Ed.)more » « less
A 3D crustal model for the central Cascadia continental shelf and Coast Range between 44°N and 45°N shows that the crystalline crust of the forearc wedge beneath the coastline is characterized by a NW-trending, vertical slab of high-velocity rock interpreted to represent the dike complex that fed the Yachats Basalt, which was intruded into the forearc approximately 37 million years ago. A spatial correlation is observed between downward deflection of the crust of the subducting Juan de Fuca plate, inferred from inversion of PmP arrivals to image the Moho surface, and the high velocity (and consequently high density) anomaly underlying the Yachats Basalt. Apparent subsequent rebound of the subducting plate at greater depth suggests a primarily elastic response of the subducting plate to this load. Calculations for a range of plausible values for the magnitude of the load and the width and depth of the depression indicate that the effective elastic thickness of the subducted Juan de Fuca plate is < 6 km. Although our simple analytical models do not include partial support of the load of the slab by the adjacent upper plate crust or time dependence to account for the motion of the slab beneath the load, incorporation of those effects should decrease rather than increase the apparent strength of the subducted plate. We conclude that the subducted Juan de Fuca plate beneath the central Oregon margin is elastically thin and has the potential to store elastic strain energy before rupturing. Our model of a well-defined, focused and static upper plate load that locally deforms the subducted plate within the nominally seismogenic or transitional part of the Cascadia plate boundary may be unique in providing a relatively straightforward scenario for estimating the mechanical properties of the subducted Juan de Fuca plate. We extrapolate from these results to speculate that elastic deformation of the subducting plate may contribute to the low level of seismicity throughout much of the Cascadia forearc in the inter-seismic period between great earthquakes but note that our local results do not preclude faulting or elasto-plastic deformation of a thin and weak plate as it subducts. These results also suggest that the subducting plate should deform in response to larger scale variations in upper plate thickness and density.
-
null (Ed.)International Ocean Discovery Program (IODP) Expedition 372 combines two research topics, slow slip events (SSEs) on subduction faults (IODP Proposal 781A-Full) and actively deforming gas hydrate–bearing landslides (Proposal 841-APL). Our study area on the Hikurangi margin east of New Zealand provides unique locations for addressing both research topics. Gas hydrates have long been suspected of being involved in seafloor failure; not much evidence, however, has been found to date for gas hydrate–related submarine landslides. Solid, icelike gas hydrate in sediment pores is generally thought to increase seafloor strength, as confirmed by a number of laboratory measurements. Dissociation of gas hydrate to water and overpressured gas, on the other hand, may destabilize the seafloor, potentially causing submarine landslides. The Tuaheni Landslide Complex on the Hikurangi margin shows evidence for active, creeping deformation. Intriguingly, the landward edge of creeping coincides with the pinchout of the base of gas hydrate stability (BGHS) on the seafloor. We therefore hypothesize that gas hydrate may be linked to creeping by (1) repeated small-scale sliding at the BGHS, in a variation of the conventional model linking gas hydrates and seafloor failure; (2) overpressure at the BGHS due to a permeability reduction linked to gas hydrates, which may lead to hydrofracturing, weakening the seafloor and allowing transmission of pressure into the gas hydrate stability zone; or (3) icelike viscous deformation of gas hydrates in sediment pores, similar to onshore rock glaciers. The latter two processes imply that gas hydrate itself is involved in creeping, constituting a paradigm shift in relating gas hydrates to submarine slope failure. Alternatively, creeping may not be related to gas hydrates but instead be caused by repeated pressure pulses or linked to earthquake-related liquefaction. We have devised a coring and logging program to test our hypotheses. SSEs at subduction zones are an enigmatic form of creeping fault behavior. At the northern Hikurangi subduction margin (HSM), they are among the best-documented and shallowest on Earth. They recur about every 2 y and may extend close to the trench, where clastic and pelagic sediments about 1.0–1.5 km thick overlie the subducting, seamount-studded Hikurangi Plateau. The northern HSM thus provides an excellent setting to use IODP capabilities to discern the mechanisms behind slow slip fault behavior, as proposed in IODP Proposal 781A-Full. The objectives of Proposal 781A-Full will be implemented across two related IODP expeditions, 372 and 375. Expedition 372 will undertake logging while drilling (LWD) at three sites targeting the upper plate (midslope basin, proposed Site HSM-01A), the frontal thrust (proposed Site HSM-18A), and the subducting section in the trench (proposed Site HSM-05A). Expedition 375 will undertake coring at the same sites, as well as an additional seamount site on the subducting plate, and implement the borehole observatory objectives. The data from each expedition will be shared between both scientific parties. Collectively, the LWD and coring data will be used to (1) characterize the compositional, structural, thermal, and diagenetic state of the incoming plate and the shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock associated with SSEs at greater depth, and (2) characterize the material properties, thermal regime, and stress conditions in the upper plate above the SSE source region. These data will be used during Expedition 375 to guide the installation of CORK observatories at the frontal thrust and in the upper plate above the SSE source to monitor temporal variations in deformation, fluid flow, seismicity, and physical and chemical properties throughout the SSE cycle (Saffer et al., 2017). Together, these data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface.more » « less
-
Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100 km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggz558
