More Like this

1. Receiver Function Derived Structural Constraints on Dynamic Processes Associated with the Young Nazca Lithosphere Subducting beneath Colombia

   Bishop, Brandon; Cho, Sung-Won; Warren, Linda M; Pedraza, Patricia; Prieto, German A; Dionicio, Viviana (January 2024, AGU Fall Meeting 2023)

   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
2. Raising the Resurrection plate from an unfolded-slab plate tectonic reconstruction of northwestern North America since early Cenozoic time

   https://doi.org/10.1130/B35677.1  

   Fuston, Spencer; Wu, Jonny (October 2020, GSA Bulletin)

   null (Ed.)

   The configuration of mid-ocean ridges subducted below North America prior to Oligocene time is unconstrained by seafloor isochrons and has been primarily inferred from upper-plate geology, including near-trench magmatism. However, many tectonic models are permitted from these constraints. We present a fully kinematic, plate tectonic reconstruction of the NW Cordillera since 60 Ma built by structurally unfolding subducted slabs, imaged by mantle tomography, back to Earth’s surface. We map in three-dimensions the attached Alaska and Cascadia slabs, and a detached slab below western Yukon (Canada) at 400−600 km depth that we call the “Yukon Slab.” Our restoration of these lower plates within a global plate model indicates the Alaska slab accounts for Pacific-Kula subduction since ca. 60 Ma below the Aleutian Islands whereas the Cascadia slab accounts for Farallon subduction since at least ca. 75 Ma below southern California, USA. However, intermediate areas show two reconstruction gaps that persist until 40 Ma. We show that these reconstruction gaps correlate spatiotemporally to published NW Cordillera near-trench magmatism, even considering possible terrane translation. We attribute these gaps to thermal erosion related to ridge subduction and model mid-ocean ridges within these reconstruction gap mid-points. Our reconstructions show two coeval ridge-trench intersections that bound an additional “Resurrection”-like plate along the NW Cordillera prior to 40 Ma. In this model, the Yukon slab represents a thermally eroded remnant of the Resurrection plate. Our reconstructions support a “northern option” Farallon ridge geometry and allow up to ∼1200 km Chugach terrane translation since Paleocene time, providing a new “tomographic piercing point” for the Baja-British Columbia debate. 

   more » « less
3. Evidence for Deeply Subducted Lower‐Plate Seamounts at the Hikurangi Subduction Margin: Implications for Seismic and Aseismic Behavior

   https://doi.org/10.1029/2021JB022866  

   Chow, Bryant; Kaneko, Yoshihiro; Townend, John (January 2022, Journal of Geophysical Research: Solid Earth)

   Abstract Seamounts are found at many subduction zones and act as seafloor heterogeneities that affect slip behavior on megathrusts. At the Hikurangi subduction zone offshore the North Island, New Zealand, seamounts have been identified on the incoming Pacific plate and below the accretionary prism, but there is little concrete evidence for seamounts subducted beyond the present‐day coastline. Using a high‐resolution, adjoint tomography‐derived velocity model of the North Island, we identify two high‐velocity anomalies below the East Coast and an intraslab low‐velocity zone up‐dip of one of these anomalies. We interpret the high‐velocity anomalies as previously unidentified, deeply subducted seamounts, and the low‐velocity zone as fluid in the subducting slab. The seamounts are inferred to be 10–30 km wide and on the plate interface at 12–15 km depth. Resolution analysis using point spread functions confirms that these are well‐resolved features. The locations of the two seamounts coincide with bathymetric features whose geometries are consistent with those predicted from analog experiments and numerical simulations of seamount subduction. The spatial characteristics of seismicity and slow slip events near the inferred seamounts agree well with previous numerical modeling predictions of the effects of seamount subduction on megathrust stress and slip. Anomalous geophysical signatures, magnetic anomalies, and swarm seismicity have also been observed previously at one or both seamount locations. We propose that permanent fracturing of the northern Hikurangi upper plate by repeated seamount subduction may be responsible for the dichotomous slow slip behavior observed geodetically, and partly responsible for along‐strike variations in plate coupling on the Hikurangi subduction interface. 

   more » « less
4. Thermal modeling of subduction zones with prescribed and evolving 2D and 3D slab geometries

   https://doi.org/10.1186/s40645-024-00611-4  

   Sime, Nathan; Wilson, Cian R.; van Keken, Peter E. (March 2024, Progress in Earth and Planetary Science)

   Abstract The determination of the temperature in and above the slab in subduction zones, using models where the top of the slab is precisely known, is important to test hypotheses regarding the causes of arc volcanism and intermediate-depth seismicity. While 2D and 3D models can predict the thermal structure with high precision for fixed slab geometries, a number of regions are characterized by relatively large geometrical changes over time. Examples include the flat slab segments in South America that evolved from more steeply dipping geometries to the present day flat slab geometry. We devise, implement, and test a numerical approach to model the thermal evolution of a subduction zone with prescribed changes in slab geometry over time. Our numerical model approximates the subduction zone geometry by employing time dependent deformation of a Bézier spline that is used as the slab interface in a finite element discretization of the Stokes and heat equations. We implement the numerical model using the FEniCS open source finite element suite and describe the means by which we compute approximations of the subduction zone velocity, temperature, and pressure fields. We compute and compare the 3D time evolving numerical model with its 2D analogy at cross-sections for slabs that evolve to the present-day structure of a flat segment of the subducting Nazca plate. 

   more » « less
5. Assessing the Role of Water in Alaskan Flat‐Slab Subduction

   https://doi.org/10.1029/2021GC009734  

   Petersen, Sarah E.; Hoisch, Thomas D.; Porter, Ryan C. (May 2021, Geochemistry, Geophysics, Geosystems)

   Abstract Low‐angle subduction has been shown to have a profound impact on subduction processes. However, the mechanisms that initiate, drive, and sustain flat‐slab subduction are debated. Within all subduction zone systems, metamorphic dehydration reactions within the down‐going slab have been hypothesized to produce seismicity, and to produce water that fluxes melting of the asthenospheric wedge leading to arc magmatism. In this work, we examine the role hydration plays in influencing slab buoyancy and the geometry of the downgoing oceanic plate. When water is introduced to the oceanic lithosphere, it is incorporated into hydrous phases, which results in lowered rock densities. The net effect of this process is an increase in the buoyancy of the downgoing oceanic lithosphere. To better understand the role of water in low‐angle subduction settings, we model flat‐slab subduction in Alaska, where the thickened oceanic lithosphere of the Yakutat oceanic plateau is subducting beneath the continental lithosphere. In this work, we calculate the thermal conditions and stable mineral assemblages in the slab crust and mantle in order to assess the role that water plays in altering the density of the subducting slab. Our slab density results show that a moderate amount of hydration (1–1.5 wt% H₂O) in the subducting crust and upper lithospheric mantle reduces slab density by 0.5%–0.8% relative to an anhydrous slab, and is sufficient to maintain slab buoyancy to 300–400 km from the trench. These models show that water is a viable factor in influencing the subduction geometry in Alaska, and is likely important globally. 

   more » « less