An official website of the United States government
Abstract Mid‐lithosphere discontinuities are seismic interfaces likely located within the lithospheric mantle of stable cratons, which typically represent velocities decreasing with depth. The origins of these interfaces are poorly understood due to the difficulties in both characterizing them seismically and reconciling the observations with thermal‐chemical models of cratons. Metasomatism of the cratonic lithosphere has been reported by numerous geochemical and petrological studies worldwide, yet its seismic signature remains elusive. Here, we identify two distinct mid‐lithosphere discontinuities at ∼87 and ∼117 km depth beneath the eastern Wyoming craton and the southwestern Superior craton by analyzing seismic data recorded by two longstanding stations. Our waveform modeling shows that the shallow and deep interfaces represent isotropic velocity drops of 2%–8% and 4%–9%, respectively, depending on the contributions from changes in radial anisotropy and density. By building a thermal‐chemical model including the regional xenolith thermobarometry constraints and the experimental phase‐equilibrium data of mantle metasomatism, we show that the shallow interface probably represents the metasomatic front, below which hydrous minerals such as amphibole and phlogopite are present, whereas the deep interface may be caused by the onset of carbonated partial melting. The hydrous minerals and melts are products of mantle metasomatism, with CO2‐H2O‐rich siliceous melt as a probable metasomatic reagent. Our results suggest that mantle metasomatism is probably an important cause of mid‐lithosphere discontinuities worldwide, especially near craton boundaries, where the mantle lithosphere may be intensely metasomatized by fluids and melts released by subducting slabs.
more »
« less
Global Patterns in Cratonic Mid‐Lithospheric Discontinuities From Sp Receiver Functions
Abstract The goal of this study is to constrain the origins of layering in the seismic velocity structure within the cratonic mantle lithosphere (i.e. mid‐lithospheric discontinuities [MLDs]). For long‐lived stations in cratons worldwide, we calculated S‐to‐P converted phase receiver function stacks using time domain deconvolution and a k‐means algorithm to select robust, consistent receiver functions. Negative MLDs appear in only 50% of the receiver function stacks, indicating that negative MLDs are common but intermittent. The negative MLDs correspond to shear velocity drops of 1%–4%, which could be caused by layers of minerals created by metasomatism, although vertical layering in seismic anisotropy cannot be ruled out. In craton interiors, negative MLDs have a lower amplitude (<3% velocity drops) and can be explained by metasomatism of the original Archean mantle. Negative MLD amplitudes increase with decreasing upper mantle shear velocity (toward the outer margins of the cratons), but do not depend on the age of the craton. Thus, negative MLD amplitudes are not dominated by age‐related variations in the cratonic mantle composition, and, instead, are more strongly correlated with proximity to tectonic and metasomatic activity that occurred long after craton formation. Negative MLDs are less numerous among stations that have Paleoproterozoic and Archean thermotectonic ages, consistent with the view that shallow release of slab‐derived fluids during early “warm” subduction was less favorable for negative MLD formation. We also observe velocity gradients below 150 km at stations in craton boundaries and interiors, indicating the presence of seismic velocity changes at the cratonic lithosphere‐asthenosphere boundary and/or Lehmann discontinuity.
more »
« less
- Award ID(s):
- 1763243
- PAR ID:
- 10449336
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 22
- Issue:
- 6
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Layering within the cratonic lithosphere has been explored and reported in different cratons using a range of techniques. However, whether there exists a common feature in the lithosphere for all the cratons is not clear yet. In this study, we carry out a comparison study between the Yilgarn craton in Western Australia and the Superior craton in North America that have never been in direct contact throughout their tectonic history. To have a detailed description of the lithospheric layering in both cratons, we employ receiver function analysis with harmonic decomposition to characterize the anisotropic seismic structure beneath 4 long‐operating sites in each craton. We can identify multiple unique anisotropic boundaries above 170 km at all sites in both cratons. Properties of the anisotropic boundaries are distinct both within and across the cratons. Our observation agrees with a commonly accepted view of the cratonic lithosphere consisting of at least two layers. Moreover, it adds new details to the previous view and reveals lateral variations of the anisotropic properties over distances of a few hundreds of kilometers. Such variations in anisotropic properties likely reflect the tectonic history predating the final assembly of cratons, and suggest horizontal movements are necessary for the formation of cratonic lithosphere.more » « less
-
Abstract The lithospheric structure of the contiguous US and surrounding regions offers clues into the tectonic history, including interactions between subducting slabs and cratons. In this paper, we present a new radially anisotropic shear wave speed model of the upper mantle (70–410 km) of the contiguous US and surrounding regions, constrained by seismic full‐waveform inversion. The new model (named CUSRA2021) utilizes frequency‐dependent travel time measurements, from 160 earthquake events recorded by 5,280 stations. The data coverage in eastern US is improved by incorporating more intraplate earthquakes. The final model exhibits clear and detailed shear wave speed anomalies correlating well with tectonic units such as North America Craton (high‐Vs), Cascadia subduction zones (high‐Vs), Columbia Plateau (low‐Vs), Basin and Range (low‐Vs), etc. In particular, the detailed structure of the North America Craton beneath Illinois basin is revealed. The depth of high‐Vs anomaly beneath the North America Craton correlates well with S‐to‐P receiver function and SH reflection results. Besides, the radial anisotropy in the Craton lithosphere shows a layering structure, which may relate to the process of lithospheric accretion and the origin of mid‐lithosphere discontinuities.more » « less
-
The cratonic lithosphere is a vast host for deep recycled carbon, trapping up to several weight percent CO2 at depths overlapping the seismic mid-lithospheric discontinuities (MLDs). However, the role of carbonates, especially for the latest discovered amorphous calcium carbonate (CaCO3), is underestimated in the formation of MLDs. Using the pulse-echo-overlap method in a Paris-Edinburgh press coupled with synchrotron X-ray diffraction, we explored the acoustic velocities of CaCO3 under high pressure-temperature (P-T) conditions relevant to the cratonic lithosphere. Two anomalous velocity drops were observed associated with the phase transition from aragonite to amorphous phase and with the pressure-induced velocity drop in the amorphous phase around 3 GPa, respectively. Both drops are comparable with approximately 35% and 52% reductions for compressional (VP) and shear (VS) wave velocities, respectively. The VP and VS values of the amorphous CaCO3 above 3 GPa are about 1/2 and 1/3 of those of the major upper-mantle minerals, respectively. These velocity reductions caused by the presence of CaCO3 would readily cause MLDs at depths of 70–120 km dependent on the geotherm even if only 1–2 vol.% CaCO3 is present in the cratonic lithosphere.more » « less
-
Abstract Lithospheric layering contains critical information about continental formation and evolution. However, discrepancies on the depth distributions of lithospheric layers have significantly limited our understanding of possible tectonic connections among the layers. Here, we construct a high‐resolution shear velocity model of eastern North America using full‐wave ambient noise simulation and inversion by integrating onshore and offshore seismic datasets. Our new model reveals large lateral variations of lithosphere thickness approximately across the major tectonic boundaries, strong low‐velocity anomalies underlying the thinner lithosphere, and multiple low‐velocity layers within the continental lithosphere. We suggest that the present mantle lithosphere beneath eastern North America was formed and modified through multiple stages of tectonic processes, among which metasomatism may have significantly contributed to the observed intralithospheric low‐velocity layers. The sharp thickness variation of lithosphere promoted edge‐driven mantle convection, which has been consequently modifying the overlying mantle lithosphere and further sharpening the gradient of lithosphere thicknessmore » « less
