Flat-slab subduction occurs when the descending plate becomes horizontal at some depth before resuming its descent into the mantle. It is often proposed as a mechanism for the uplifting of deep crustal rocks ('thick-skinned' deformation) far from plate boundaries, and for causing unusual patterns of volcanism, as far back as the Proterozoic eon. For example, the formation of the expansive Rocky Mountains and the subsequent voluminous volcanism across much of the western USA has been attributed to a broad region of flat-slab subduction beneath North America that occurred during the Laramide orogeny (80-55 million years ago). Here we study the largest modern flat slab, located in Peru, to better understand the processes controlling the formation and extent of flat slabs. We present new data that indicate that the subducting Nazca Ridge is necessary for the development and continued support of the horizontal plate at a depth of about 90 kilometres. By combining constraints from Rayleigh wave phase velocities with improved earthquake locations, we find that the flat slab is shallowest along the ridge, while to the northwest of the ridge, the slab is sagging, tearing, and re-initiating normal subduction. On the basis of our observations, we propose a conceptual model for the temporal evolution of the Peruvian flat slab in which the flat slab forms because of the combined effects of trench retreat along the Peruvian plate boundary, suction, and ridge subduction. We find that while the ridge is necessary but not sufficient for the formation of the flat slab, its removal is sufficient for the flat slab to fail. This provides new constraints on our understanding of the processes controlling the beginning and end of the Laramide orogeny and other putative episodes of flat-slab subduction.

We have determined the Wadati-Benioff Zone seismicity and state of stress of the subducting Nazca slab beneath central and southern Peru using data from three recently deployed local seismic networks. Our relocated hypocenters are consistent with a flat slab geometry that is shallowest near the Nazca Ridge, and changes from steep to normal without tearing to the south. These locations also indicate numerous abrupt along-strike changes in seismicity, most notably an absence of seismicity along the projected location of subducting Nazca Ridge. This stands in stark contrast to the very high seismicity observed along the Juan Fernandez ridge beneath central Chile where, a similar flat slab geometry is observed. We interpret this as indicative of an absence of water in the mantle beneath the overthickened crust of the Nazca Ridge. This may provide important new constraints on the conditions required to produce intermediate depth seismicity. Our focal mechanisms and stress tensor inversions indicate dominantly down-dip extension, consistent with slab pull, with minor variations that are likely due to the variable slab geometry and stress from adjacent regions. We observe significantly greater variability in the P-axis orientations and maximum compressive stress directions. The along strike change in the orientation of maximum compressive stress are likely related to slab bending and unbending south of the Nazca Ridge.

Flat or near-horizontal subduction of oceanic lithosphere has been an important tectonic process both currently and in the geologic past. Subduction of the aseismic Nazca Ridge beneath South America has been associated with the onset of flat subduction and the termination of arc volcanism in Peru, making it an ideal place to study flat-slab subduction. Recently acquired seismic recordings for 144 broadband seismic stations in Peru permit us to image the Mohorovičić discontinuity (Moho) of the subducted oceanic Nazca plate, Nazca Ridge, and the overlying continental Moho of the South American crust in detail through the calculation of receiver functions. We find that the subducted over-thickened ridge crust is likely significantly eclogitized ~350 km from the trench, requiring that the inboard continuation of the flat slab be supported by mechanisms other than low-density crustal material. This continuation coincides with a low-velocity anomaly identified in prior tomography studies of the region immediately below the flat slab, and this anomaly may provide some support for the flat slab. The subduction of the Nazca Ridge has displaced most, if not the entire South American lithospheric mantle beneath the high Andes as well as up to 10 km of the lowermost continental crust. The lack of deep upper-plate seismicity suggests that the Andean crust has remained warm during flat subduction and is deforming ductilely around the subducted ridge. This deformation shows significant coupling between the subducting Nazca oceanic plate and overriding South American continental plate up to ~500 km from the trench. These results provide important modern constraints for interpreting the geological consequences of past and present flat-slab subduction locations globally.

The occurrence of intermediate depth seismicity (70-300 km) is commonly attributed to the dehydration of hydrous phases within the downgoing oceanic plate. While some water is incorporated into the oceanic crust at formation, a significant amount of water is introduced into the plate immediately before subduction along outer-rise faults. These faults have been shown to extend to depths of over 30 km and can channel water to depths of 20 km or more beneath the seafloor. However, the amount of water introduced into the oceanic mantle lithosphere, and the role of that water in the formation of intermediate depth seismicity, has been the topic of ongoing research.Here we compile evidence from areas where the subducted oceanic crust is likely thicker than the penetration depth of water into the downgoing plate. These regions comprise aseismic plateaus and ridges (hot spot tracks) that can be compared directly to adjacent segments of the oceanic plate where oceanic crust of normal thickness is subducted. Regions with thick oceanic crust show little to no seismicity at intermediate depths, whereas adjacent regions with normal oceanic crust (∼6-8 km thick) have significant seismicity at similar depths and distances from the trench. We hypothesize that intermediate depth earthquakes observed in regions with thinner oceanic crust are caused by mantle dehydration reactions that are not possible in regions where the oceanic mantle was never hydrated because the thickness of the oceanic crust exceeded the penetration depth of water into the plate. We compare our observations to phase diagrams of hydrous basalt and hydrated depleted peridotite to determine pressures and temperatures that would be consistent with our observations. These can provide valuable constraints, not only on the degree of hydration and dehydration in the downgoing plate, but also as ground-truth for thermal models of these regions, all of which have complex, three-dimensional, time-variant subduction geometries and thermal histories.

No abstract