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
This content will become publicly available on December 1, 2024
Insights From Layered Anisotropy Beneath Southern New England: From Ancient Tectonism to Present‐Day Mantle Flow
Seismic anisotropy beneath eastern North America, as expressed in shear wave splitting observations, has been attributed to plate motion‐parallel shear in the asthenosphere, resulting in fast axes aligned with the plate motion. However, deviations of fast axes from plate motion directions are observed near major tectonic boundaries of the Appalachians, indicating contributions from lithospheric anisotropy associated with past tectonic processes. In this study, we conduct anisotropic receiver function (RF) analysis using data from a dense seismic array traversing the New England Appalachians in Connecticut to examine anisotropic layers in the crust and upper mantle and correlate them with past tectonic processes as well as present‐day mantle flow. We use the harmonic decomposition method to separate directionally‐dependent variations of RFs and focus on features with the same harmonic signals observed across multiple stations. Within the crust, there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen‐parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra‐crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion‐parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area.
more » « less- Award ID(s):
- 1800923
- NSF-PAR ID:
- 10488118
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 24
- Issue:
- 12
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Seismic anisotropy beneath eastern North America likely reflects both the remnant lithospheric fabrics and the present‐day deformation of the asthenosphere. We report new observations of splitting in core‐refracted shear phases observed over 3–5 years at 33 sites in New Jersey, New York, and states in the New England region and also include data from eight previously studied locations. Our data set emphasizes back azimuthal coverage necessary to capture the directional variation of splitting parameters expected from vertically varying anisotropy. We report single‐phase splitting parameters as well as station‐averaged values based on splitting intensity technique that incorporates all observed records regardless of whether they showed evidence of splitting or not. Trends of averaged fast shear wave polarizations appear coherent and are approximately aligned with absolute plate motion direction. The general disparity between the fast axes and the trend of surface tectonic features suggests a dominant asthenosphere contribution for the observed seismic anisotropy. Averaged delay values systematically increase from ~0.5 s in New Jersey to ~1.4 s in Maine. Splitting parameters vary at all sites, and neighboring stations often show similar patterns of directional variation. We developed criteria to group stations based on their splitting patterns and identified four domains with distinct anisotropic properties. Splitting patterns of three domains suggest a layered anisotropic structure that is geographically variable, outlining distinct regions in the continental mantle, for example, the Proterozoic lithosphere of the Adirondack Mountains. A domain coincident with the North Appalachian Anomaly displays virtually no splitting, implying that the lithospheric fabric was locally erased.
-
more » « less
SUMMARY We present a new, 3-D model of seismic velocity and anisotropy in the Pacific upper mantle, PAC13E. We invert a data set of single-station surface-wave phase-anomaly measurements sensitive only to Pacific structure for the full set of 13 anisotropic parameters that describe surface-wave anisotropy. Realistic scaling relationships for surface-wave azimuthal anisotropy are calculated from petrological information about the oceanic upper mantle and are used to help constrain the model. The strong age dependence in the oceanic velocities associated with plate cooling is also used as a priori information to constrain the model. We find strong radial anisotropy with vSH > vSV in the upper mantle; the signal peaks at depths of 100–160 km. We observe an age dependence in the depth of peak anisotropy and the thickness of the anisotropic layer, which both increase with seafloor age, but see little age dependence in the depth to the top of the radially anisotropic layer. We also find strong azimuthal anisotropy, which typically peaks in the asthenosphere. The azimuthal anisotropy at asthenospheric depths aligns better with absolute-plate-motion directions while the anisotropy within the lithosphere aligns better with palaeospreading directions. The relative strengths of radial and azimuthal anisotropy are consistent with A-type olivine fabric. Our findings are generally consistent with an explanation in which corner flow at the ridge leads to the development and freezing-in of anisotropy in the lithosphere, and shear between the lithosphere and underlying asthenosphere leads to anisotropy beneath the plate. We also observe large regions within the Pacific basin where the orientation of anisotropy and the absolute-plate-motion direction differ; this disagreement suggests the presence of shear in the asthenosphere that is not aligned with absolute-plate-motion directions. Azimuthal-anisotropy orientation rotates with depth; the depth of the maximum vertical gradient in the fast-axis orientation tends to be age dependent and agrees well with a thermally controlled lithosphere–asthenosphere boundary. We observe that azimuthal-anisotropy strength at shallow depths depends on half-spreading rate, with higher spreading rates associated with stronger anisotropy. Our model implies that corner flow is more efficient at aligning olivine to form lattice-preferred orientation anisotropy fabrics in the asthenosphere when the spreading rate at the ridge is higher.
-
null (Ed.)Abstract The eastern margin of North America has been affected by a range of fundamental tectonic processes in the geologic past. Major events include the Paleozoic Appalachian orogeny, which culminated in the formation of the supercontinent Pangea, and the breakup of Pangea during the Mesozoic. The southern New England Appalachians exhibit a particularly rich set of geologic and tectonic structures that reflect multiple episodes of subduction and terrane accretion, as well as subsequent continental breakup. It remains poorly known, however, to what extent structures at depth in the crust and lithospheric mantle reflect these processes, and how they relate to the geological architecture at the surface. The Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn) was a deployment of 15 broadband seismometers in a dense linear array across northern Connecticut. The array traversed a number of major tectonic boundaries, sampling across the Laurentian margin in its western portion to the Avalonian terrane at its eastern end. It also crossed the Hartford rift basin in the central portion of the state. The SEISConn stations operated between 2015 and 2019; data from the experiment are archived at the Incorporated Research Institutions for Seismology Data Management Center and will be publicly available beginning in 2021. A suite of imaging techniques is being applied to SEISConn data, with the goal of providing a detailed view of the crust and mantle lithosphere (including discontinuities, seismic velocities, and seismic anisotropy) beneath the southern New England Appalachians. Results from these analyses will inform a host of fundamental scientific questions about the structural evolution of orogens, the processes involved in continental rifting, and the nature of crustal and mantle lithospheric deformation during subduction, terrane accretion, and continental breakup.more » « less
-
more » « less
Abstract To systematically investigate seismic azimuthal anisotropy in the Sumatra subduction zone and probe mantle dynamics associated with the subduction of the Australian Plate beneath the Sunda Plate, a total of 169 pairs of teleseismic XKS (including PKS, SKKS, SKS) and 115 pairs of local
S splitting parameters are obtained using broadband seismic data recorded at ~70 stations. Additionally, crustal anisotropy in the overriding Sunda Plate is measured by analyzing the moveout ofP ‐to‐S conversions from the Moho using a sinusoidal function. Comparison between the three sets of anisotropy measurements obtained using shear waves with different depths of origin suggests that (1) the crust of the Sunda Plate is anisotropic with mostly trench‐parallel fast orientations and a mean splitting time of 0.28 ± 0.05 s; (2) the mantle wedge is azimuthally anisotropic with dominantly trench‐parallel fast orientations and splitting times ranging from 0.22 to 0.81 s, which generally increase with the focal depth; and (3) subslab anisotropy is mostly trench‐normal beneath the fore‐arc region with an averaged splitting time of 1.48 ± 0.06 s, and becomes trench‐parallel beneath the arc and back‐arc areas with a mean splitting time of 0.33 ± 0.04 s. The resulting lateral and vertical distributions of anisotropy obtained using splitting of three types of shear waves advocate the presence of an entrained subslab flow that is deflected by the mantle transition zone. The flow enters the mantle wedge through a slab window and flows horizontally parallel to the trench.
