A magnetotelluric survey in the Barotse Basin of western Zambia shows clear evidence for thinned lithosphere beneath an orogenic belt. The uppermost asthenosphere, at a depth of 60–70 km, is highly conductive, suggestive of the presence of a small amount of partial melt, despite the fact that there is no surface expression of volcanism in the region. Although the data support the presence of thicker cratonic lithosphere to the southeast of the basin, the lithospheric thickness is not well resolved and models show variations ranging from ~80 to 150 km in this region. Similarly variable is the conductivity of the mantle beneath the basin and immediately beneath the cratonic lithosphere to the southeast, although the conductivity is required to be elevated compared to normal lithospheric mantle. In a general sense, two classes of model are compatible with the magnetotelluric data: one with a moderately conductive mantle and one with more elevated conductivities. This latter class would be consistent with the impingement of a stringer of plume‐fed melt beneath the cratonic lithosphere, with the melt migrating upslope to thermally erode lithosphere beneath the orogenic belt that is overlain by the Barotse Basin. Such processes are potentially important for intraplate volcanism and also for development or propagation of rifting as lithosphere is thinned and weakened by melt. Both models show clear evidence for thinning of the lithosphere beneath the orogenic belt, consistent with elevated heat flow data in the region.
This study probes the lithosphere‐asthenosphere system beneath 155 Ma Pacific seafloor using teleseismic S‐to‐p receiver functions at the Pacific Lithosphere Anisotropy and Thickness Experiment project ocean‐bottom‐seismometers. Within the lithosphere, a significant velocity decrease at 33–50 km depth is observed. This mid‐lithospheric discontinuity is consistent with the velocity contrast between the background mantle and thin, trapped layers of crystallized partial melt, in the form of either dolomite or garnet granulite. These melts possibly originated from deeper asthenospheric melting beneath the flanks of spreading centers, and were transported within the cooling lithosphere. A positive velocity increase of 3%–6% is observed at 130–155 km depth and is consistent with the base of a layer with partial melt in the asthenosphere. A shear velocity decrease associated with the lithosphere‐asthenosphere boundary at 95–115 km depth is permitted by the data, but is not required.
more » « less- PAR ID:
- 10498306
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 51
- Issue:
- 5
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract The thermal regime of continental lithosphere plays a fundamental role in controlling the behavior of tectonic plates. In this work, we assess the thermal state of the North American upper mantle by combining shear‐wave velocity models, calculated using data from the EarthScope facility, with empirically derived anelasticity models and basalt thermobarometry. We estimate the depth of the thermal lithosphere‐asthenosphere boundary (LAB), defined as the intersection of a geotherm with the 1300°C adiabat. Results show lithospheric thicknesses across the contiguous US vary between ∼40 km and >200 km. The thinnest thermal lithosphere is observed in the tectonically active western US and the thickest lithosphere in the midcontinent. By combining geotherm estimates with solidus curves for peridotite, we show that a pervasive partial melt zone is common within the western US upper mantle and that partial melt is absent in the eastern and central US without significant metasomatism.
-
more » « less
Abstract The Colorado Plateau and its surroundings serve as an archetypal case to investigate the interaction of mantle melting processes and lithospheric structure. It has been hypothesized that widespread Cenozoic volcanism indicates the encroachment of the convective upwelling of asthenosphere toward the Plateau center. In this study, we generate a Common Conversion Point (CCP) stack of S‐to‐p (Sp) receiver functions to image the locations of lithospheric discontinuities in the southwestern United States. Our results are broadly similar to prior work, showing a strong and continuous Negative Velocity Gradient (NVG) consistent with the Lithosphere‐Asthenosphere Boundary (LAB) over much of the study area. However, with several methodological improvements, we are able to obtain more reliable NVG depth picks below the Colorado Plateau where the LAB becomes weaker, deeper, and broader. We compare the inferred topography of NVGs with the locations of volcanoes, and find that the majority of recent volcanoes are co‐located with lithosphere that is ∼80 km thick. This appears to be the critical depth at which partial melt from upwelling asthenosphere pooling at the base of (or within) the lithosphere may percolate to the surface. We compare our CCP profiles with magma equilibration conditions determined from petrologic analysis and find good agreement between the depth of NVGs and depth of magma equilibration. This analysis provides insight into the progression of magmatism and lithospheric loss toward the center of the Colorado Plateau, and demonstrates how small‐scale processes like melting influence lithosphere‐asthenosphere interactions that persist over large temporal and spatial scales.
-
more » « less
Abstract The relative motion of the lithosphere with respect to the asthenosphere implies the existence of a boundary zone that accommodates shear between the rigid plates and flowing mantle. This shear zone is typically referred to as the lithosphere‐asthenosphere boundary (LAB). The width of this zone and the mechanisms accommodating shear across it have important implications for coupling between mantle convection and surface plate motion. Seismic observations have provided evidence for several physical mechanisms that might help enable relative plate motion, but how these mechanisms each contribute to the overall accommodation of shear remains unclear. Here we present receiver function constraints on the discontinuity structure of the oceanic upper mantle at the NoMelt site in the central Pacific, where local constraints on shear velocity, anisotropy, conductivity, and attenuation down to ∼300 km depth provide a comprehensive picture of upper mantle structure. We image a seismic discontinuity with a Vsv decrease of 4.5% or more over a 0–20 km thick gradient layer centered at a depth of ∼65 km. We associate this feature with the Gutenberg discontinuity (G), and interpret our observation of G as resulting from strain localization across a dehydration boundary based on the good agreement between the discontinuity depth and that of the dry solidus. Transitions in Vsv, azimuthal anisotropy, conductivity, and attenuation observed at roughly similar depths suggest that the G discontinuity represents a region of localized strain within a broader zone accommodating shear between the lithosphere and asthenosphere.
-
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.
