An official website of the United States government
SUMMARY Seismic interrogation of the upper mantle from the base of the crust to the top of the mantle transition zone has revealed discontinuities that are variable in space, depth, lateral extent, amplitude and lack a unified explanation for their origin. Improved constraints on the detectability and properties of mantle discontinuities can be obtained with P-to-S receiver function (Ps-RF) where energy scatters from P to S as seismic waves propagate across discontinuities of interest. However, due to the interference of crustal multiples, uppermost mantle discontinuities are more commonly imaged with lower resolution S-to-P receiver function (Sp-RF). In this study, a new method called CRISP-RF (Clean Receiver-function Imaging using SParse Radon Filters) is proposed, which incorporates ideas from compressive sensing and model-based image reconstruction. The central idea involves applying a sparse Radon transform to effectively decompose the Ps-RF into its underlying wavefield contributions, that is direct conversions, multiples, and noise, based on the phase moveout and coherence. A masking filter is then designed and applied to create a multiple-free and denoised Ps-RF. We demonstrate, using synthetic experiment, that our implementation of the Radon transform using a sparsity-promoting regularization outperforms the conventional least-squares methods and can effectively isolate direct Ps conversions. We further apply the CRISP-RF workflow on real data, including single station data on cratons, common-conversion-point stack at continental margins and seismic data from ocean islands. The application of CRISP-RF to global data sets will advance our understanding of the enigmatic origins of the upper mantle discontinuities like the ubiquitous mid-lithospheric discontinuity and the elusive X-discontinuity.
more »
« less
A Taxonomy of Upper‐Mantle Stratification in the US
Abstract The investigation of upper mantle structure beneath the US has revealed a growing diversity of discontinuities within, across, and underneath the sub‐continental lithosphere. As the complexity and variability of these detected discontinuities increase—for example, velocity increase/decrease, number of layers and depth—it is hard to judge which constraints are robust and which explanatory models generalize to the largest set of constraints. Much work has been done to image discontinuities of interest using S‐waves that convert to P‐waves (or top‐side reflected SS waves). A higher resolution method using P‐to‐S scattered waves is preferred but often obscured by multiply reflected waves trapped in a shallower layer, limiting the visibility of deeper boundaries. Here, we address the interference problem and re‐evaluate upper mantle stratification using filtered P‐to‐S receiver functions (Ps‐RFs) interpreted using unsupervised machine‐learning. Robust insight into upper mantle layering is facilitated with CRISP‐RF: Clean Receiver‐Function Imaging using Sparse Radon Filters. Subsequent sequencing and clustering organizes the polarity‐filtered Ps‐RFs into distinct depth‐based clusters. We find three types of upper mantle stratification beneath the old and stable continental US: (a) intra‐lithosphere discontinuities (paired or single boundary), (b) transitional discontinuities (single boundary or with a top layer), and (c) sub‐lithosphere discontinuities. Our findings contribute a more nuanced understanding of mantle discontinuities, offering new perspectives on the nature of upper mantle layering beneath continents.
more »
« less
- Award ID(s):
- 1818654
- PAR ID:
- 10539901
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 5
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Summary To provide constraints on a number of significant controversial issues related to the structure and dynamics of the Australian continent, we utilize P-to-S receiver functions (RFs) recorded by 182 stations to map the 410 and 660 km discontinuities (d410 and d660, respectively) bordering the mantle transition zone (MTZ). The RFs are stacked in successive circular bins with a radius of 1o under a non-plane wavefront assumption. The d410 and d660 depths obtained using the 1-D IASP91 Earth model show a systematic apparent uplifting of about 15 km for both discontinuities in central and western Australia relative to eastern Australia, as the result of higher seismic wavespeeds in the upper mantle beneath the former area. After correcting the apparent depths using the Australian Seismological Reference Model, the d410 depths beneath the West Australia Craton are depressed by ∼10 km on average relative to the normal depth of 410 km, indicating a positive thermal anomaly of 100 K at the top of the MTZ which could represent a transition from a thinner than normal MTZ beneath the Indian ocean and the normal MTZ beneath central Australia. The abnormally thick MTZ beneath eastern Australia can be adequately explained by subducted cold slabs in the MTZ. A localized normal thickness of the MTZ beneath the Newer Volcanics Province provides supporting evidence of non-mantle-plume mechanism for intraplate volcanic activities in the Australian continent.more » « less
-
Abstract In this study, we use data from the SEISConn seismic experiment to calculate Sp receiver functions in order to characterize the geometry of upper-mantle structure beneath southern New England (northeastern United States). We image robust negative-velocity-gradient discontinuities beneath southern New England that we interpret as corresponding to the lithosphere-asthenosphere boundary (LAB) and identify a well-defined step of 15 km in LAB depth at a longitude of 73°W, which we interpret to be the boundary between Laurentian and Appalachian lithosphere, although the offset may be larger if the putative LAB phase is reinterpreted to be a mid-lithospheric discontinuity. We infer that the lithosphere throughout the region is substantially thinner than elsewhere in the continental interior, consistent with regional tomographic studies and previously published Sp receiver function results. The presence of thinned lithosphere suggests that the low-velocity Northern Appalachian Anomaly (NAA) in the upper mantle may extend as far south as coastal Connecticut. The presence of regionally thinned lithosphere and a step in lithospheric thickness suggests that inherited structure may be preserved in present-day lithosphere, even in the presence of more recent dynamic processes associated with the NAA.more » « less
-
Abstract Abrupt velocity gradients in the upper mantle, detectable by receiver functions (RF) techniques, have been known to exist down to the depths of ∼110 km beneath northeastern North America. Comparisons with the surface wave velocity models have designated some negative velocity gradients (NVGs) as the lithosphere‐asthenosphere boundary (LAB), delineating a relatively thin lithosphere beneath this region. This work presents a systematic survey of upper mantle layering in seismic properties using P‐S RF analysis at 62 long‐running sites with dense lateral sampling. We examine both radial and transverse component RF for indicators of seismic anisotropy and adopt the notion of seismic attributes, utilized in active‐source seismology, to characterize the spatial distribution of directionally variant and invariant signal components. We confirm a widespread presence of NVGs at depths 60–100 km throughout the region, consistent with previous studies using mode‐converted body waves. We also find numerous converting boundaries that reflect changes in directional variation (anisotropy) of seismic velocity, indicating complexity of rock texture in the upper mantle. Some of these boundaries appear as deep as 185 km, implying that the lithosphere extends much deeper than the widespread NVGs would suggest. In this, our results agree with recent estimates of the lithospheric thickness in thermodynamically consistent models combining seismic, gravity, and heat flow constraints.more » « less
-
Abstract Lithospheric discontinuities, including the lithosphere‐asthenosphere boundary (LAB) and the enigmatic mid‐lithospheric discontinuities (MLDs), hold important clues about the structure and evolution of tectonic plates. However, P‐ and S‐receiver‐function (PRF and SRF) techniques, two traditional techniques to image Earth's deep discontinuities, have some shortcomings in imaging lithosphere discontinuities. Here, we propose a new method using reflections generated by teleseismic S waves (hereafter S‐reflections) to image lithospheric discontinuities, which are less affected by multiple phases than PRFs and have better depth resolution than SRFs. We apply this method to the data collected by the Transportable Array and other regional seismic networks and obtain new high‐resolution images of the lithosphere below the contiguous US. Beneath the tectonically active Western US, we observe a negative polarity reflector (NPR) in the depth range of 60–110 km, with greatly varying amplitude and depth, which correlates with active tectonic processes. We interpret this feature as the LAB below the Western US. Beneath the tectonically stable Central and Eastern US, we observe two NPRs in the depth ranges of 60–100 km and 100–150 km, whose amplitude and depth also vary significantly, and which appear to correlate with past tectonic processes. We interpret these features as MLDs below the Central and Eastern US. Our results show reasonable agreement with results from PRFs, which have similar depth resolution, suggesting the possibility of joint inversion of S‐reflections and PRFs to constrain the properties of lithospheric discontinuities.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2024JB028781
