Award Abstract # 2041631
Characterization of the Oceanic Lithosphere-Asthenosphere Boundary in the Colombian Subduction Zone through Receiver Function Analysis

NSF Org: EAR
Division Of Earth Sciences
Recipient: SAINT LOUIS UNIVERSITY
Initial Amendment Date: July 25, 2021
Latest Amendment Date: July 25, 2021
Award Number: 2041631
Award Instrument: Standard Grant
Program Manager: Eva Zanzerkia
EAR
 Division Of Earth Sciences
GEO
 Directorate for Geosciences
Start Date: August 1, 2021
End Date: July 31, 2024 (Estimated)
Total Intended Award Amount: $243,553.00
Total Awarded Amount to Date: $243,553.00
Funds Obligated to Date: FY 2021 = $243,553.00
History of Investigator:
  • Brandon Bishop (Principal Investigator)
    brandon.t.bishop@slu.edu
  • Linda Warren (Co-Principal Investigator)
Recipient Sponsored Research Office: Saint Louis University
221 N GRAND BLVD
SAINT LOUIS
MO  US  63103-2006
(314)977-3925
Sponsor Congressional District: 01
Primary Place of Performance: Saint Louis University
221 N. Grand Blvd.
St Louis
MO  US  63103-2006
Primary Place of Performance
Congressional District:
01
Unique Entity Identifier (UEI): JNBLLTBTLLD8
Parent UEI: JNBLLTBTLLD8
NSF Program(s): Marine Geology and Geophysics,
XC-Crosscutting Activities Pro
Primary Program Source: 010V2122DB R&RA ARP Act DEFC V
Program Reference Code(s): 102Z
Program Element Code(s): 162000, 722200
Award Agency Code: 4900
Fund Agency Code: 4900
Assistance Listing Number(s): 47.050

ABSTRACT

This award is funded in whole or in part under the American Rescue Plan Act of 2021 (Public Law 117-2).

The goal of this project is to test two competing models of the boundary separating the Earth?s stiff outer layer (lithosphere) from the underlying, weaker layer (asthenosphere). One model proposes that this boundary is a narrow channel that decouples the lithosphere and asthenosphere. The other model proposes that this boundary forms a broad band of deformation that couples the lithosphere and asthenosphere. Each model has important implications for our understanding of plate tectonics?the process which drives the Earth?s long-term geologic evolution?and its relation to the Earth?s interior. Each model carries predictions about the boundary between the lithosphere and asthenosphere. These predictions can be tested using receiver functions. In collaboration with the Universidad Nacional de Colombia, this project uses recordings from seismic stations of the Colombian National Seismic Network located above the Colombian subduction zone. In this region, the characteristics of Nazca oceanic plate?s lower boundary can be examined over a range of depths to determine which of the models is most accurate. Understanding which model is most accurate will benefit society by helping to constrain the force balance acting on tectonic plates, especially around subduction zones where potentially destructive earthquakes driven by these forces occur. This project further benefits society by deepening international collaboration between institutions in the US and Colombia as well as furthering STEM education by providing training in seismology to a postdoctoral fellow.

Observations of the oceanic lithosphere-asthenosphere boundary (LAB) support two categories of contradictory models for coupling between the lithosphere and underlying asthenosphere, and correspondingly contradictory implications for plate tectonics. One category interprets the boundary to be a narrow, melt-rich channel that decouples lithospheric plates from the asthenosphere. Plate motion is thus driven largely by slab pull. The other category interprets the boundary as a broad zone which couples plates to the asthenosphere. Plate motion is thus strongly influenced by asthenospheric flow induced drag. These models predict differing structures and localizations of shear along the LAB. Channel models predict two closely spaced, highly sheared boundaries while coupled models predict either a single boundary with little shear or a single highly sheared boundary. The predicted features respond differently to increasing pressure due to subduction. Channels narrow with depth and coupled structures persist to greater, differing depths, so observing a plate as it subducts offers an ideal natural experiment to test these models. Seismic waves with different frequency contents interact with these structures such that high frequencies can better image channels while low frequencies can better image coupled structures. This project uses P-to-S-wave receiver function analysis of radial and transverse components recorded at broadband seismic stations to examine structural characteristics and shear localization along the boundary. The complete characterization of structure and anisotropy caused by shear requires long-term observations at seismic stations from the coast to the back arc at a range of frequency bands. This project uses data from 2010 to the present from 26 stations in the Colombian National Seismic Network to investigate and characterize the subducting Nazca oceanic plate?s LAB from depths of <100 km at the coastline to >400 km in the back arc. These observations will provide constraints on the depth, thickness, and anisotropic properties (if present) of the boundary and use these to test the proposed models.

This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.

PUBLICATIONS PRODUCED AS A RESULT OF THIS RESEARCH

Note:  When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

Bishop, Brandon and Cho, Sung-Won and Warren, Linda M and Pedraza, Patricia and Prieto, German A and Dionicio, Viviana "Receiver Function Derived Structural Constraints on Dynamic Processes Associated with the Young Nazca Lithosphere Subducting beneath Colombia" AGU Fall Meeting 2023 , 2024 Citation Details
Bishop, Brandon and Fernandez, Ademar and Warren, Linda M. and Pedraza, Patricia and Prieto, German A. and Dionicio, Viviana "A Thick, Sheared Zone at the Base of an Oceanic Plate: The Lithosphere-Asthenosphere Boundary of the Subducting Nazca Plate in Colombia as Revealed through P-to-S Receiver Functions" AGU Fall Meeting , 2022 Citation Details

PROJECT OUTCOMES REPORT

Disclaimer

This Project Outcomes Report for the General Public is displayed verbatim as submitted by the Principal Investigator (PI) for this award. Any opinions, findings, and conclusions or recommendations expressed in this Report are those of the PI and do not necessarily reflect the views of the National Science Foundation; NSF has not approved or endorsed its content.

Tectonic plates represent the rigid outermost ~100 kilometer layer of the Earth. This outer layer is called the lithosphere. The lithosphere lies above a weaker layer of the Earth called the asthenosphere. The characteristics of the boundary between the lithosphere and asthenosphere, called the LAB, controls whether the lithosphere moves independently of the asthenosphere or whether motion in one layer causes motion in the other. If the LAB is narrow and dominated by melt, the two layers will move independently. If the LAB is broad with limited or no melt, motion in the lithosphere will cause drag in the upper part of the asthenosphere, slowing both.

Oceanic lithosphere is much younger than continental lithosphere, making it comparatively simpler to study as it has been affected by fewer geologic events. Oceanic lithosphere is also pulled to increasing depths during subduction, increasing pressure at the LAB while largely leaving the temperature contrast between the lithosphere and asthenosphere unaffected. This creates a natural experiment where the LAB is subjected to steadily increasing pressures. If melt dominates the LAB, it will tend to be forced upward by its buoyancy, making the LAB narrower with increasing depth. If drag dominates the LAB, there should be little effect on its thickness with depth. Drag affects different minerals in different ways, and for the main mineral making up the lithosphere and asthenosphere above ~410 km depth, drag causes an alignment of mineral grains which causes seismic waves traveling in different directions to travel faster or slower than they normally would (this is called anisotropy). Below 410 km depth, the dominant mineral does not cause anisotropy. If drag is an important characteristic of the LAB, anisotropy should affect seismic waves crossing the LAB until ~410 km depth.

We used data from the Colombian National Seismic Network (the circles and diamonds in the accompanying figure mark the location of the network’s seismometers) to investigate the LAB of the oceanic Nazca plate’s lithosphere as it subducts from the surface to >400 km depth below Colombia. The colored lines in the accompanying figure mark the estimated depth to the top of the Nazca plate’s subducted lithosphere. The Nazca plate is very young, and the age of the plate is shown in millions of years along the left side of the figure. The plate’s estimated motion with respect to the material below the asthenosphere is marked by a white arrow in the accompanying figure. If drag at a large scale is important for the asthenosphere, we might expect seismic anisotropy to approximately align with the arrow.

Our data was used to calculate receiver functions, a technique that uses seismic waves from distant earthquakes to detect boundaries in the Earth directly below a seismometer. We found that the Nazca plate’s LAB below all but 6 seismometers in the network has clear evidence for anisotropy. We further found that there is no LAB detectable at the seismometers above the deepest part of the subducted plate. Both of these observations indicate drag is the key characteristic for the LAB of this young plate. We further found that the orientation of anisotropy below the slab in the central and southern part of Colombia is highly variable and does not fit a simple model where anisotropy is mostly caused by the dip of the subducting plate (the cyan and magenta rays pointing out from the station symbol in the map key show the pattern that should be present if the dip of the plate is the primary cause of anisotropy).

This grant also supported a Master’s student in developing new receiver function techniques to potentially expand the usable data set and a post-doctoral researcher who helped mentor the Master’s student and carried out the main part of the project. Results from this grant were presented at the American Geophysical Union’s Fall Meeting in 2021-2023, and three manuscripts are being prepared to further disseminate our results. The >35,000 receiver functions calculated for this study will be publicly released for use in other studies alongside these manuscripts. In addition, the grant supported collaboration between seismologists at Saint Louis University, the Colombian Geological Survey, and the National University of Colombia.


Last Modified: 11/19/2024
Modified by: Brandon T Bishop

Please report errors in award information by writing to: awardsearch@nsf.gov.

Print this page

Back to Top of page