Subduction zones are home to the most seismically active faults on the planet. The shallow megathrust interfaces of subduction zones host Earth's largest earthquakes and are likely the only faults capable of magnitude 9+ ruptures. Despite these facts, our knowledge of subduction zone geometry-which likely plays a key role in determining the spatial extent and ultimately the size of subduction zone earthquakes-is incomplete. We calculated the three-dimensional geometries of all seismically active global subduction zones. The resulting model, called Slab2, provides a uniform geometrical analysis of all currently subducting slabs.

Using finite-frequency teleseismic P-wave tomography, we developed a new three-dimensional (3-D) velocity model of the mantle beneath Anatolia down to 900 km depth that reveals the structure and behavior of the subducting African lithosphere beneath three convergent domains of Anatolia: the Aegean, Cyprean, and Bitlis-Zagros domains. The Aegean slab has a relatively simple structure and extends into the lower mantle; the Cyprean slab has a more complex structure, with a western section that extends to the lower mantle with a consistent dip and an eastern section that is broken up into several pieces; and the Bitlis slab appears severely deformed, with only fragments visible in the mantle transition zone and uppermost lower mantle. In addition to the subducting slabs, high-amplitude slow velocity anomalies are imaged in the shallow mantle beneath recently active volcanic centers, and a prominent fast velocity anomaly dominates the shallow mantle beneath northern Anatolia and the southern Black Sea. As a whole, our model confirms the presence of well-established slow and fast velocity anomalies in the upper mantle beneath Anatolia and motivates two major findings about Eastern Mediterranean subduction: (1) Each of the slabs penetrates into the lower mantle, making the Eastern Mediterranean unique within the Mediterranean system, and (2) the distinct character of each slab segment represents different stages of subduction termination through progressive slab deformation. Our findings on the destructive processes of subduction termination and slab detachment have significant implications for understanding of the post detachment behavior of subducted lithosphere.

Slow seismic velocity anomalies are commonly imaged beneath subducting slabs in tomographic studies, yet a unifying explanation for their distribution has not been agreed upon. In South America two such anomalies have been imaged associated with subduction of the Nazca Ridge in Peru and the Juan Fernández Ridge in Chile. Here we present new seismic images of the subslab slow velocity anomaly beneath Chile, which give a unique view of the nature of such anomalies. Slow seismic velocities within a large hole in the subducted Nazca slab connect with a subslab slow anomaly that appears correlated with the extent of the subducted Juan Fernández Ridge. The hole in the slab may allow the subslab material to rise into the mantle wedge, revealing the positive buoyancy of the slow material. We propose a new model for subslab slow velocity anomalies beneath the Nazca slab related to the entrainment of hot spot material.

Nazca subduction beneath South America is one of our best modern examples of long‐lived ocean‐continent subduction on the planet, serving as a foundation for our understanding of subduction processes. Within that framework, persistent heterogeneities at a range of scales in both the South America and Nazca plates is difficult to reconcile without detailed knowledge of the subducted Nazca slab structure. Here we use teleseismic travel time residuals from >1,000 broadband and short‐period seismic stations across South America in a single tomographic inversion to produce the highest‐resolution contiguous P wave tomography model of the subducting slab and surrounding mantle beneath South America to date. Our model reveals a continuous trench‐parallel fast seismic velocity anomaly across the majority of South America that is consistent with the subducting Nazca slab. The imaged anomaly indicates a number of robust features of the subducted slab, including variable slab dip, extensive lower mantle penetration, slab stagnation in the lower mantle, and variable slab amplitude, that are incorporated into a new, comprehensive model of the geometry of the Nazca slab surface to ~1,100 km depth. Lower mantle slab penetration along the entire margin suggests that lower mantle slab anchoring is insufficient to explain along strike upper plate variability while slab stagnation in the lower mantle indicates that the 1,000 km discontinuity is dominant beneath South America.

Summary The Andean Subduction Zone is one of the longest continuous subduction zones on Earth. The relative simplicity of the two-plate system has makes it an ideal natural laboratory to study the dynamics in subduction zones. We measure teleseismic S and SKS travel-time residuals at > 1,000 seismic stations that have been deployed across South America over the last 30 years to produce a finite-frequency teleseismic S-wave tomography model of the mantle beneath the Andean Subduction Zone related to the Nazca Plate, spanning from ∼5° N to 45° S and from depths of ∼130 km to 1,200 km. Within our model, the subducted Nazca slab is imaged as a fast velocity seismic anomaly. The geometry and amplitude of the Nazca slab anomaly varies along the margin while the slab anomaly continues into the lower mantle along the entirety of the subduction margin. Beneath northern Brazil, the Nazca slab appears to stagnate at ∼1,000 km depth and extend eastward sub-horizontally for > 2,000 km. South of 25° S the slab anomaly in the lower mantle extends offshore of eastern Argentina, hence we do not image if a similar stagnation occurs. We image several distinct features surrounding the slab including two vertically-oriented slow seismic velocity anomalies: one beneath the Peruvian flat slab and the other beneath the Paraná Basin of Brazil. The presence of the latter anomaly directly adjacent to the stagnant Nazca slab suggests that the plume, known as the Paraná Plume, may be a focused upwelling formed in response to slab stagnation in the lower mantle. Additionally, we image a high amplitude fast seismic velocity anomaly beneath the Chile trench at the latitude of the Sierras Pampeanas which extends from ∼400 km to ∼1000 km depth. This anomaly may be the remnants of an older, detached slab, however its relationship with the Nazca-South America subduction zone remains enigmatic.

The April 2016, Pedernales Earthquake ruptured a 100 km by 40 km segment of the subduction zone along the coast of Ecuador in an Mw 7.8 megathrust event east of the intersection of the Carnegie Ridge with the trench. This portion of the subduction zone has ruptured on decadal time scales in similar size and larger earthquakes, and exhibits a range of slip behaviors, variations in segmentation, and degree of plate coupling along strike. Immediately after the earthquake, an international rapid response effort coordinated by the Instituto Geofísico at the Escuela Politécnica Nacional in Quito deployed 55 seismometers and 10 OBS above the rupture zone and adjacent areas to record aftershocks. In this article we describe the details of the US portion of the rapid response and present an earthquake catalog from May 2016-May 2017 produced using data recorded by these stations. Aftershocks focus in distinct clusters within and around the rupture area and match spatial patterns observed in long term seismicity. For the first two and a half months, aftershocks exhibit a relatively sharp cutoff to the north of the mainshock rupture. In early July an earthquake swarm occurred ~100 km to the northeast of the mainshock in the epicentral region of a Mw 7.8 earthquake in 1958. In December, an increase in seismicity occurred ~70 km to the northeast of the mainshock in the epicentral region of the 1906 earthquake. Data from the Pedernales earthquake and aftershock sequence recorded by permanent seismic and geodetic networks in Ecuador and the dense aftershock deployment provide an opportunity to examine the persistence of asperities for large to great earthquakes over multiple seismic cycles, the role of asperities and slow slip in subduction zone Meltzer et al., SRL Data Mine: Pedernales RAPID Deployment 3 megathrust rupture, and the relationship between locked and creeping parts of the subduction interface.

A tear in the subducting Nazca slab is located between the end of the Pampean flat slab and normally subducting oceanic lithosphere. Tomographic studies suggest mantle material flows through this opening. The best way to probe this hypothesis is through observations of seismic anisotropy, such as shear wave splitting. We examine patterns of shear wave splitting using data from two seismic deployments in Argentina that lay updip of the slab tear. We observe a simple pattern of plate‐motion‐parallel fast splitting directions, indicative of plate‐motion‐parallel mantle flow, beneath the majority of the stations. Our observed splitting contrasts previous observations to the north and south of the flat slab region. Since plate‐motion‐parallel splitting occurs only coincidentally with the slab tear, we propose mantle material flows through the opening resulting in Nazca plate‐motion‐parallel flow in both the subslab mantle and mantle wedge.

No abstract