skip to main content


This content will become publicly available on January 24, 2025

Title: Receiver Function Derived Structural Constraints on Dynamic Processes Associated with the Young Nazca Lithosphere Subducting beneath Colombia
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
Award ID(s):
2041631
NSF-PAR ID:
10487638
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
AGU Fall Meeting 2023
Date Published:
Journal Name:
AGU Fall Meeting 2023
Subject(s) / Keyword(s):
["subduction","Nazca plate","South America","geology","geophysics","seismology","receiver functions"]
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The structure of the lithosphere-asthenosphere boundary (LAB) beneath oceanic plates is key to understanding how plates interact with the underlying mantle. Prior contradictory geophysical observations have been used to argue for a thin, melt-rich boundary that decouples the plate from the rest of the mantle, or for a much broader anisotropic and thermally controlled boundary that indicates significant coupling with the rest of the mantle. The predictions of models based on these interpretations can be tested most easily in a subduction zone setting where the steady increase in pressure at the base of the subducting plate’s LAB will have differing effects on melt and anisotropy. Melt remains stable within the mantle to ~150-250 km (for carbonate melt) or to ~330 km (for silicate melt), while anisotropy induced by different processes should have no significant change until ~250 km to ~440 km depth. We calculate P-to-S receiver functions (PRFs) using varying frequency bands at broadband seismic stations with >4 years of data from the Servicio Geológico Colombiano’s Red Sismológica Nacional de Colombia to investigate the characteristics of the LAB of the subducting Nazca oceanic plate from the coast to the Andean foreland (corresponding to slab LAB depths of ~50 km to >400 km). The use of PRFs permits identification and analysis of anisotropy across the boundary while calculation at a range of frequency bands permits tuning of the PRFs to differing spatial scales to determine the size and abruptness of the boundary. We find that the P-to-S converted phase of the subducted Nazca plate’s LAB is detectable 4-5 seconds after the converted phase of the plate’s Moho to at least ~150 km depth. Assuming the slab has an average Vp/Vs of 1.75 to 1.78 and Vp of 8.2 km/s (+2.5% dVp), this corresponds to a plate thickness of ~50 km, matching the expected thickness given the Nazca plate’s age in the region (~10-20 Myrs). We find that the Nazca plate’s LAB is most consistently detectable in the <0.24 Hz band and largely undetectable in the <2.4 Hz band, indicating the LAB is gradational and between 10 and 30 km in thickness. Amplitude variations and complexities in the LAB converted phases further indicate that the boundary marks a change in anisotropy most consistent with the LAB representing a sheared zone between the plate and underlying mantle. 
    more » « less
  2. SUMMARY

    The Nazca-South America subduction zone in Ecuador is characterized by a complicated along-strike geometry as the slab transitions from flat slab subduction in the south, with the Peruvian flat slab, to what has been characterized as ‘normal’ dipping subduction beneath central Ecuador. Plate convergence additionally changes south to north as the trench takes on a convex shape. Highly heterogeneous bathymetry at the trench, including the aseismic oceanic Carnegie Ridge (CR), and sparse intermediate-depth seismicity has led many to speculate about the behaviour of the downgoing plate at depth. In this study, we present a finite-frequency teleseismic P-wave tomography model of the northern Andes beneath Ecuador and Colombia from 90 to 1200 km depth. Our model builds on prior tomography models in South America by adding relative traveltime residuals recorded at stations in Ecuador. The complete data set is comprised of 114 096 relative traveltime residuals from 1133 stations across South America, with the added data serving to refine the morphology of the Nazca slab in the mantle beneath the northern Andes. Our tomography model shows a Nazca slab with a fragmented along-strike geometry and the first teleseismic images of several proposed slab tears in this region. At the northern edge of the Peruvian flat slab in southern Ecuador, we image a shallow tear at 95–200 km depth that appears to connect mantle flow from beneath the flat slab to the Ecuadorian Arc. Beneath central Ecuador at the latitudes of the CR, the Nazca slab is continuous into the lower mantle. Beneath southern Colombia, the Malpelo Tear breaks the Nazca slab below ∼200 km depth.

     
    more » « less
  3. null (Ed.)
    Geologic processes at convergent plate margins control geochemical cycling, seismicity, and deep biosphere activity in subduction zones and suprasubduction zone lithosphere. International Ocean Discovery Program (IODP) Expedition 366 was designed to address the nature of these processes in the shallow to intermediate depth of the Mariana subduction channel. Although no technology is available to permit direct sampling of the subduction channel of an intraoceanic convergent margin at depths up to 18 km, the Mariana forearc region (between the trench and the active volcanic arc) provides a means to access this zone. Active conduits, resulting from fractures in the forearc, are prompted by along- and across-strike extension that allows slab-derived fluids and materials to ascend to the seafloor along associated faults, resulting in the formation of serpentinite mud volcanoes. Serpentinite mud volcanoes of the Mariana forearc are the largest mud volcanoes on Earth. Their positions adjacent to or atop fault scarps on the forearc are likely related to the regional extension and vertical tectonic deformation in the forearc. Serpentinite mudflows at these volcanoes include serpentinized forearc mantle clasts, crustal and subducted Pacific plate materials, a matrix of serpentinite muds, and deep-sourced formation fluid. Mud volcanism on the Mariana forearc occurs within 100 km of the trench, representing a range of depths and temperatures to the downgoing plate and the subduction channel. These processes have likely been active for tens of millions of years at this site and for billions of years on Earth. At least 10 active serpentinite mud volcanoes have been located in the Mariana forearc. Two of these mud volcanoes are Conical and South Chamorro Seamounts, which are the furthest from the Mariana Trench at 86 and 78 km, respectively. Both seamounts were cored during Ocean Drilling Program (ODP) Legs 125 and 195, respectively. Data from these two seamounts represent deeper, warmer examples of the continuum of slab-derived materials as the Pacific plate subducts, providing a snapshot of how slab subduction affects fluid release, the composition of ascending fluids, mantle hydration, and the metamorphic paragenesis of subducted oceanic lithosphere. Data from the study of these two mud volcanoes constrain the pressure, temperature, and composition of fluids and materials within the subduction channel at depths of about 18 to 19 km. Understanding such processes is necessary for elucidating factors that control seismicity in convergent margins, tectonic and magma genesis processes in the forearc and volcanic arc, fluid and material fluxes, and the nature and variability of environmental conditions that impact subseafloor microbial communities. Expedition 366 centered on data collection from cores recovered from three serpentinite mud volcanoes that define a continuum of subduction-channel processes defined by the two previously cored serpentinite mud volcanoes and the trench. Three serpentinite mud volcanoes (Yinazao, Fantangisña, and Asùt Tesoro) were chosen at distances 55 to 72 km from the Mariana Trench. Cores were recovered from active sites of eruption on their summit regions and on the flanks where ancient flows are overlain by more recent ones. Recovered materials show the effects of dynamic processes that are active at these sites, bringing a range of materials to the seafloor, including materials from the lithosphere of the Pacific plate and from subducted seamounts (including corals). Most of the recovered material consists of serpentinite mud containing lithic clasts, which are derived from the underlying forearc crust and mantle and the subducting Pacific plate. Cores from each of the three seamounts drilled during Expedition 366, as well as those from Legs 125 and 195, include material from the underlying Pacific plate. A thin cover of pelagic sediment was recovered at many Expedition 366 sites, and at Site U1498 we cored through serpentinite flows to the underlying pelagic sediment and volcanic ash deposits. Recovered serpentinites are largely uniform in major element composition, with serpentinized ultramafic rocks and serpentinite muds spanning a limited range in SiO2 , MgO, and Fe2 O3 compositions. However, variation in trace element composition reflects pore fluid composition, which differs as a function of the temperature and pressure of the underlying subduction channel. Dissolved gases H2 , CH4 , and C2 H6 are highest at the site furthest from the trench, which also has the most active fluid discharge of the Expedition 366 serpentinite mud volcanoes. These dissolved gases and their active discharge from depth likely support active microbial communities, which were the focus of in-depth subsampling and preservation for shore-based analytical and culturing procedures. The effects of fluid discharge were also registered in the porosity and GRA density data indicated by higher than expected values at some of the summit sites. These higher values are consistent with overpressured fluids that minimize compaction of serpentinite mud deposits. In contrast, flank sites have significantly greater decreases in porosity with depth, suggesting that processes in addition to compaction are required to achieve the observed data. Thermal measurements reveal higher heat flow values on the flanks (~31 mW/m2) than on the summits (~17 mW/m2) of the seamounts. The new 2G Enterprises superconducting rock magnetometer (liquid helium free) revealed relatively high values of both magnetization and bulk magnetic susceptibility of discrete samples related to ultramafic rocks, particularly in dunite. Magnetite, a product of serpentinization, and authigenic carbonates were observed in the mudflow matrix materials. In addition to coring operations, Expedition 366 focused on the deployment and remediation of borehole casings for future observatories and set the framework for in situ experimentation. Borehole work commenced at South Chamorro Seamount, where the original-style CORK was partially removed. Work then continued at each of the three summit sites following coring operations. Cased boreholes with at least three joints of screened casing were deployed, and a plug of cement was placed at the bottom of each hole. Water samples were collected from two of the three boreholes, revealing significant inputs of formation fluids. This suggests that each of the boreholes tapped a hydrologic zone, making these boreholes suitable for experimentation with the future deployment of a CORK-lite. An active education and outreach program connected with many classrooms on shore and with the general public through social media. 
    more » « less
  4. null (Ed.)
    Geologic processes at convergent plate margins control geochemical cycling, seismicity, and deep biosphere activity in subduction zones and suprasubduction zone lithosphere. International Ocean Discovery Program Expedition 366 was designed to address the nature of these processes in the shallow to intermediate depth of the Mariana subduction channel. Although no technology is available to permit direct sampling of the subduction channel of an intraoceanic convergent margin at depths up to 19 km, the Mariana forearc region (between the trench and the active volcanic arc) provides a means to access materials from this zone. Active conduits, resulting from fractures in the forearc, are prompted by along- and across-strike extension that allows slab-derived fluids and materials to ascend to the seafloor along associated faults, resulting in the formation of serpentinite mud volcanoes. Serpentinite mud volcanoes of the Mariana forearc are the largest mud volcanoes on Earth. Their positions adjacent to or atop fault scarps on the forearc are likely related to the regional extension and vertical tectonic deformation in the forearc. Serpentinite mudflows at these volcanoes include serpentinized forearc mantle clasts, crustal and subducted Pacific plate materials, a matrix of serpentinite muds, and deep-sourced formation fluid. Mud volcanism on the Mariana forearc occurs within 100 km of the trench, representing a range of depths and temperatures to the downgoing plate and the subduction channel. These processes have likely been active for tens of millions of years at the Mariana forearc and for billions of years on Earth. At least 19 active serpentinite mud volcanoes have been located in the Mariana forearc. Two of these mud volcanoes are Conical and South Chamorro Seamounts, which are the farthest from the Mariana Trench at 86 and 78 km, respectively. Both seamounts were cored during Ocean Drilling Program Legs 125 and 195, respectively. Data from these two seamounts represent deeper, warmer examples of the continuum of slab-derived materials as the Pacific plate subducts, providing a snapshot of how slab subduction affects fluid release, the composition of ascending fluids, mantle hydration, and the metamorphic paragenesis of subducted oceanic lithosphere. Data from the study of these two mud volcanoes constrain the pressure, temperature, and composition of fluids and materials within the subduction channel at depths of up to 19 km. Understanding such processes is necessary for elucidating factors that control seismicity in convergent margins, tectonic and magma genesis processes in the volcanic arc and backarc areas, fluid and material fluxes, and the nature and variability of environmental conditions that impact subseafloor microbial communities. Expedition 366 focused on data collection from cores recovered from three serpentinite mud volcanoes that define a continuum of subduction-channel processes to compare with results from drilling at the two previously cored serpentinite mud volcanoes and with previously collected gravity, piston, and remotely operated vehicle push cores across the trench-proximal forearc. Three serpentinite mud volcanoes (Yinazao, Fantangisña, and Asùt Tesoro) were chosen at distances 55 to 72 km from the Mariana Trench. Cores were recovered from active sites of eruption on their summit regions and on the flanks where ancient flows are overlain by more recent ones. Recovered materials show the effects of dynamic processes that are active at these sites, bringing a range of materials to the seafloor, including materials from the crust of the Pacific plate, most notably subducted seamounts (even corals). Most of the recovered material consists of serpentinite mud containing lithic clasts, which are derived from the underlying forearc crust and mantle and the subducting Pacific plate. A thin cover of pelagic sediment was recovered at many Expedition 366 sites, and at Site U1498 we cored through distal serpentinite mudflows and into the underlying pelagic sediment and volcanic ash deposits. Recovered serpentinized ultramafic rocks and mudflow matrix materials are largely uniform in major element composition, spanning a limited range in SiO2, MgO, and Fe2O3 compositions. However, variation in trace element composition reflects interstitial water composition, which differs as a function of the temperature and pressure of the underlying subduction channel. Dissolved gases H2, CH4, and C2H6 are highest at the site farthest from the trench, which also has the most active fluid discharge of the Expedition 366 serpentinite mud volcanoes. These dissolved gases and their active discharge from depth likely support active microbial communities, which were the focus of in-depth subsampling and preservation for shore-based analytical and culturing procedures. The effects of fluid discharge were also registered in the porosity and gamma ray attenuation density data indicated by higher than expected values at some of the summit sites. These higher values are consistent with overpressured fluids that slow compaction of serpentinite mud deposits. In contrast, flank sites have significantly greater decreases in porosity with depth, suggesting that processes in addition to compaction are required to achieve the observed data. Thermal measurements reveal higher heat flow values on the flanks (~31 mW/m2) than on the summits (~17 mW/m2) of the seamounts. The new 2G Enterprises superconducting rock magnetometer (liquid helium free) revealed relatively high values of both magnetization and bulk magnetic susceptibility of discrete samples related to ultramafic rocks, particularly dunite. Magnetite, a product of serpentinization, and authigenic carbonates were observed in the mudflow matrix materials. In addition to coring operations, Expedition 366 focused on the deployment and remediation of borehole casings for future observatories and set the framework for in situ experimentation. Borehole work commenced at South Chamorro Seamount, where the original-style CORK was partially removed. Work then continued at each of the three summit sites following coring operations. Cased boreholes with at least three joints of screened casing were deployed, and a plug of cement was placed at the bottom of each hole. Water samples were collected from two of the three boreholes, revealing significant inputs of formation fluids. This suggests that each of the boreholes tapped a hydrologic zone, making these boreholes suitable for experimentation with the future deployment of a CORK-Lite. 
    more » « less
  5. Abstract

    An apparent gap in the Andean volcanic arc in the Pampean section of the subduction zone in Chile (~28°–33°S) marks a section of flat slab subduction. In this tectonic environment, the fate of fluids released from the subducting Nazca slab remains uncertain and the degree of their interaction with the basal layer of the continental lithosphere is poorly understood. Results of a RF investigation and forward modeling effort at three long‐running stations of the Chilean National Seismic Network allow us to constrain the position of the subducting Nazca slab and to address the physical properties of the interplate contact zone in Central Chile. Our observations suggest a transition in seismic character, from a weakly anisotropic contact in the normally subducting section north of the flat slab region to a strongly anisotropic plate contact within the flat slab region. We attribute this change to a transition from sheared olivine to serpentinized peridotite generated as a result of fluid release across the flat slab. This interpretation is supported by forward modeling synthetic RFs at each of the stations. We propose that the identified layer extends across the flat slab region, acting as a mineral reservoir that captures and, possibly, transports fluids from the dehydrating Nazca Plate as it subducts below South America. We note that the Ps converted phase at the slab interface at southernmost station GO04 suffers a 2‐s discontinuity at 180° back azimuth, consistent with a 15–20‐km scarp or kink in the Nazca slab to the south of the station.

     
    more » « less