An official website of the United States government
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
Strong Physical Contrasts Across Two Mid‐Lithosphere Discontinuities Beneath the Northwestern United States: Evidence for Cratonic Mantle Metasomatism
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
- Award ID(s):
- 2123529
- PAR ID:
- 10480126
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 4
- Issue:
- 6
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract High pressure and temperature experiments were carried out on the oxide mixtures corresponding to the bridgmanite stoichiometry under the hydrous shallow lower mantle conditions (24–25 GPa and 1673–1873 K with 5–10 wt. % of water in the starting material). Oxide mixtures investigated correspond to MgSiO3, (Mg, Fe)SiO3, (Mg, Al, Si)O3, and (Mg, Fe, Al, Si)O3. Melting was observed in all runs. Partitioning of various elements, including Mg, Fe, Si, and H is investigated. Melting under hydrous lower mantle conditions leads to increased (Mg + Fe)O/SiO2in the melt compared to the residual solids. The residual solids often contain a large amount of stishovite, and the melt contains higher (Mg,Fe)O/SiO2ratio than the initial material. (Mg + Fe)O‐rich hydrous melt could explain the low‐velocity anomalies observed in the shallow lower mantle and a large amount of stishovite in the residual solid may be responsible for the scattering of seismic waves in the mid‐lower mantle and may explain the “stishovite paradox. Since stishovite‐rich materials are formed only when silica‐rich source rock (MORB) is melted (not a typical peridotitic rock [bulk silicate Earth]), seismic scattering in the lower mantle provides a clue on the circulation of subducted MORB materials. To estimate hydrogen content, we use a new method of estimating the water content of unquenchable melts, and also propose a new interpretation of the significance of superhydrous phase B inclusions in bridgmanite. The results provide revised values of water partitioning between solid minerals and hydrous melts that are substantially higher than previous estimates.more » « less
-
The East African rift overlies one or more mantle upwellings and it traverses heterogeneous Archaean-Paleozoic lithosphere rifted in Mesozoic and Cenozoic time. We re-analyze XKS shear wave splitting at publicly available stations to evaluate models for rifting above mantle plumes. We use consistent criteria to compare and contrast both splitting direction and strength, infilling critical gaps with new data from the Turkana Depression and North Tanzania Divergence sectors of the East African rift system. Our results show large spatial variations in the amount of splitting (0.1–2.5 s), with fast axes predominantly sub-parallel to the orientation of Cenozoic rifts underlain by thinned lithosphere with and without surface magmatism. The amount of splitting increases with lithospheric thinning and magmatic modification. Nowhere are fast axes perpendicular to the rift, arguing against the development of extensional strain fabrics. Thick cratons are characterized by small amounts of splitting (≤0.5 s) with a variety of orientations that may characterize mantle plume flow. Splitting rotates to rift parallel and increases in strength over short distances into rift zones, implying a shallow depth range for the anisotropy in some places. The shallow source and correlation between splitting direction and the shape of upper mantle thin zones suggests that the combination of channel flow and oriented melt pockets contribute > 1 s to the observed splitting delays. Enhanced flow, metasomatism, and melt intrusion at the lithosphere-asthenosphere boundary suggest that fluid infiltration to the base of the lithosphere may facilitate rifting of cratonic lithosphere.more » « less
-
Abstract A comprehensive analysis of experimental data and theoretical simulations on the partial molar volume of water in silicate melt indicates that finite strain theory successfully describes the compression of the H2O component dissolved in silicate melt at high pressures and temperatures. However, because of the high compressibility of the water component, a fourth order equation of state fit is required to accurately simulate experimental results on water's volume in silicate melts at a deep upper mantle, transition zone, and lower mantle pressures. Data from previous shock compression experiments on hydrous minerals in which melting occurs along the Hugoniot are used to provide an experimental constraint on the partial molar volume of water in silicate melt at deep mantle temperatures and pressures. The equation of state of the water component indicates that, depending on elastic averaging technique, the amount of water that could be present in neutrally or negatively buoyant mafic/ultramafic melts above the 410 km seismic discontinuity is upper‐bounded at 5.6 wt%: smaller than previously inferred, and consistent with melt being confined to a narrow depth range above the 410 km discontinuity. If melt is predominantly distributed along grain boundaries in low aspect ratio films, extents of melting as low as 2% could produce observed seismic velocity reductions. The ability of the lowermost mantle to contain negatively buoyant hydrous liquids hinges on the trade‐off between iron content and hydration: at these depths, substantially higher degrees of hydration could be present within partial melts.more » « less
