Seismic imaging of the global boundaries of the mantle transition zone provides important constraints regarding mantle composition and convection. Additionally, modern hypotheses about mantle convection and increasing seismic data resources motivate attempts to image more complex regional transition zone structures such as dipping and discontinuous interfaces. Here we extend a 3‐D prestack migration method to make it more applicable for transition zone imaging with Ps receiver functions. After validation with 1‐D synthetic data, two types of hypothetical structures are adopted to demonstrate the strengths and weaknesses of different practical imaging parameters with 2‐D synthetic tests. The results show that the method can resolve dipping anomalies and laterally discontinuous low velocity layers. However, receiver spacing of 0.5° is inadequate to accurately migrate scattering from all possible angles in the transition zone at ~0.3 Hz. Application of a slowness cutoff window about the Ps raypath can mitigate artifacts due to low receiver density, but the tradeoff is a decrease in the maximum resolvable dip angle. Further synthetic tests evaluate source wavelet effects, addition of observed seismic noise to synthetic data, and imaging grid dimensions. Cumulatively, the tests indicate that generic application of migration algorithms to the limited source‐receiver distributions relevant for the transition zone should be treated cautiously. Imaging parameters, such as the slowness cutoff and wavelet, should be chosen based on synthetic tests for hypothetical targets and realistic source‐receiver geometries. Finally, observational Ps receiver functions from the USArray are migrated with the new method.
This content will become publicly available on October 12, 2025
Increasing deployment of dense arrays has facilitated detailed structure imaging for tectonic investigation, hazard assessment and resource exploration. Strong velocity heterogeneity and topographic changes have to be considered during passive source imaging. However, it is quite challenging for ray‐based methods, such as Kirchhoff migration or the widely used teleseismic receiver function, to handle these problems. In this study, we propose a 3‐D passive source reverse time migration strategy based on the spectral element method. It is realized by decomposing the time reversal full elastic wavefield into amplitude‐preserved vector P and S wavefields by solving the corresponding weak‐form solutions, followed by a dot‐product imaging condition to get images for the subsurface structures. It enables us to use regional 3‐D migration velocity models and take topographic variations into account, helping us to locate reflectors at more accurate positions than traditional 1‐D model‐based methods, like teleseismic receiver functions. Two synthetic tests are used to demonstrate the advantages of the proposed method to handle topographic variations and complex velocity heterogeneities. Furthermore, applications to the Laramie array data using both teleseismic P and S waves enable us to identify several south‐dipping structures beneath the Laramie basin in southeast Wyoming, which are interpreted as the Cheyenne Belt suture zone and agree with, and improve upon previous geological interpretations.
more » « less- PAR ID:
- 10548494
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
SUMMARY Temporal changes in seismic velocities are an important tool for tracking structural changes within the crust during transient deformation. Although many geophysical processes span the crust, including volcanic unrest and large-magnitude earthquakes, existing methods for seismic monitoring are limited to the shallow subsurface. We present an approach for deep seismic monitoring based on teleseismic receiver functions, which illuminate the crustal velocity structure from the bottom-up. Using synthetic waveform modelling, we show that receiver functions are uniformly sensitive to velocity changes throughout the crust and can locate the depth of the perturbation. We introduce a novel method based on optimal transport for measuring the non-linear time–amplitude signal variations characteristic of receiver function monitoring. We show that optimal transport enables comparison of full waveform distributions rather than relying on representative stacked waveforms. We further study a linearized version of optimal transport that renders time-warping signal variations into simple Euclidean perturbations, and use this capability to perform blind source separation in the space of waveform variations. This disentangles the effects of changes in the source–receiver path from changes in subsurface velocities. Collectively, these methods extend the reach of seismic monitoring to deep geophysical processes, and provide a tool that can be used to study heterogeneous velocity changes with different spatial extents and temporal dynamics.
-
more » « less
Abstract Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical‐component GSC velocity as an approximation for the incident source wavelet. The Moho
P ‐to‐s conversions are clearly visible on DAS receiver functions. We then derive meter‐scale arrival time measurements along the entire 20‐km‐long array. We are also able to measure path‐averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future. -
more » « less
Abstract Urban basin investigation is crucial for seismic hazard assessment and mitigation. Recent advances in robust nodal‐type sensors facilitate the deployment of large‐N arrays in urban areas for high‐resolution basin imaging. However, arrays typically operate for only one month due to the instruments' battery life, and hence, only record a few teleseismic events. This limits the number of available teleseismic events for traditional receiver function (RF) analysis‐the primary method used in sediment‐basement interface imaging in passive source seismology. Insufficient stacking of RFs from a limited number of earthquakes could, however, introduce significant biases to the results. In this study, we present a novel Bayesian array‐based Coherent Receiver Function (CRF) method that can leverage datasets from short‐term dense arrays to constrain basin geometry. We cast the RF deconvolution as a sparsity‐promoted inverse problem, in which the deconvolution at a single‐station involves the constraints from neighboring stations and multiple events. We solve the inverse problem using a trans‐dimensional Markov chain Monte Carlo Bayesian algorithm to find an ensemble of RF solutions, which provides a quantitative way of deciding which features are well resolved and warrant geological interpretation. An application in the northern Los Angeles basin demonstrates the ability of our method to produce reliable and easy‐to‐interpret RF images. The use of dense seismic networks and the state‐of‐the‐art Bayesian array‐based CRF method can provide a robust approach for subsurface structure imaging.
-
more » « less
Abstract Passive seismic methods for imaging the discontinuity structure of Earth have primarily focused on differences in vertically and radially polarized energy in the coda of earthquake‐generated body waves (e.g., receiver functions). To convert the timing of scattered wave arrivals to depth, three parameters must be known or inferred: depth or layer thickness (
H ),P ‐wave velocity (V P), andS ‐wave velocity (V S). A common way to solve for these parameters is throughH ‐κ stacking analysis, in which layer thickness and the ratio betweenV PandV S(κ ) is calculated while holding one of the velocity parameters constant. However, this assumption biases estimates of layer properties and leads to uncertainties that are not appropriately quantified. As these results are commonly used as starting models for more complex seismic or geodynamic analyses, these assumptions can propagate much further than the initial study. In this study, we introduce independent observations from body‐wave autocorrelations that can help constrain this underdetermined problem.P ‐wave autocorrelation allows for the recovery of the Moho‐reflectedP ‐wave phase from teleseismic earthquakes, which is removed during deconvolution in the calculation of receiver functions. As the Moho‐reflectedP ‐wave is independent ofV S, this constraint allows us to create a system of equations that better quantifies the thickness,V P, andV Sof a layer and produces a more appropriate estimation of associated uncertainties. We apply this to 88 seismic stations that are spatially distributed throughout the United States to obtain a model of crustal variability that is unbiased by a priori assumptions of velocity structure.
