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
Seismic Azimuthal Anisotropy Beneath a Fast Moving Ancient Continent: Constraints From Shear Wave Splitting Analysis in Australia
Seismic azimuthal anisotropy beneath Australia is investigated using splitting of the teleseismic PKS, SKKS, and SKS phases to delineate asthenospheric flow and lithospheric deformation beneath one of the oldest and fast‐moving continents on Earth. In total 511 pairs of high‐quality splitting parameters were observed at 116 seismic stations. Unlike other stable continental areas in Africa, East Asia, and North America, where spatially consistent splitting parameters dominate, the fast orientations and splitting times observed in Australia show a complex pattern, with a slightly smaller than normal average splitting time of 0.85 ± 0.33 s. On the North Australian Craton, the fast orientations are mostly N‐S, which is parallel to the absolute plate motion (APM) direction in the hotspot frame. Those observed in the South Australian Craton are mostly NE‐SW and E‐W, which are perpendicular to the maximum lithospheric horizontal shortening direction. In east Australia, the observed azimuthal anisotropy can be attributed to either APM induced simple shear or lithospheric fabric parallel to the strike of the orogenic belts. The observed spatial variations of the seismic azimuthal anisotropy, when combined with results from depth estimation utilizing the spatial coherency of the splitting parameters and seismic tomography studies, suggest that the azimuthal anisotropy in Australia can mostly be related to simple shear in the rheologically transition layer between the lithosphere and asthenosphere. Non‐APM parallel anisotropy is attributable to modulations of the mantle flow system by undulations of the bottom of the lithosphere, with a spatially variable degree of contribution from lithospheric fabric.
more » « less- Award ID(s):
- 2149587
- NSF-PAR ID:
- 10395212
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 2
- ISSN:
- 2169-9313
- 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
Abstract The breakup of supercontinent Pangea occurred ∼200 Ma forming the Eastern North American Margin (ENAM). Yet, the precise timing and mechanics of breakup and onset of seafloor spreading remain poorly constrained. We investigate the relict lithosphere offshore eastern North America using ambient‐noise Rayleigh‐wave phase velocity (12–32 s) and azimuthal anisotropy (17–32 s) at the ENAM Community Seismic Experiment (CSE). Incorporating previous constraints on crustal structure, we construct a shear velocity model for the crust and upper ∼60 km of the mantle beneath the ENAM‐CSE. A low‐velocity lid (
V S of 4.4–4.55 km/s) is revealed in the upper 15–20 km of the mantle that extends ∼200 km from the margin, terminating at the Blake Spur Magnetic Anomaly (BSMA). East of the BSMA, velocities are fast (>4.6 km/s) and characteristic of typical oceanic mantle lithosphere. We interpret the low‐velocity lid as stretched continental mantle lithosphere embedded with up to ∼15% retained gabbro. This implies that the BSMA marks successful breakup and onset of seafloor spreading ∼170 Ma, consistent with ENAM‐CSE active‐source studies that argue for breakup ∼25 Myr later than previously thought. We observe margin‐parallel Rayleigh‐wave azimuthal anisotropy (2%–4% peak‐to‐peak) in the lithosphere that approximately correlates with absolute plate motion (APM) at the time of spreading. We hypothesize that lithosphere formed during ultra‐slow seafloor spreading records APM‐modified olivine fabric rather than spreading‐parallel fabric typical of higher spreading rates. This work highlights the importance of present‐day passive margins for improving understanding of the fundamental rift‐to‐drift transition. -
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. -
more » « less
Abstract A total of 431 well‐defined and 632 null shear–wave splitting measurements obtained from 115 broadband seismic stations located in the eastern Himalayan syntaxis and adjacent areas is largely inconsistent with predicted fast orientations by absolute plate motion models. Spatial coherency analysis of the splitting parameters suggests that the observed azimuthal anisotropy is mostly located in the upper asthenosphere or the transitional layer between the lithosphere and asthenosphere, and the disagreement between the fast orientations and regional tectonic fabrics suggests an insignificant lithospheric contribution to the observed anisotropy. The observations may be attributed to flow systems that are driven by the westward rollback of the Indian slab beneath the Indo‐Burma block and are modulated by a previously revealed gap between the northward and eastward subducting slabs of the Indian plate.
