An official website of the United States government
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
Passive Source Reverse Time Migration Based on the Spectral Element Method
Abstract 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
- Award ID(s):
- 2042098
- 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
-
-
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 MohoP‐to‐sconversions 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 We develop a 3‐D isotropic shear velocity model for the Alaska subduction zone using data from seafloor and land‐based seismographs to investigate along‐strike variations in structure. By applying ambient noise and teleseismic Helmholtz tomography, we derive Rayleigh wave group and phase velocity dispersion maps, then invert them for shear velocity structure using a Bayesian Monte Carlo algorithm. For land‐based stations, we perform a joint inversion of receiver functions and dispersion curves. The forearc crust is relatively thick (35–42 km) and has reduced lower crustal velocities beneath the Kodiak and Semidi segments, which may promote higher seismic coupling. Bristol Bay Basin crust is relatively thin and has a high‐velocity lower layer, suggesting a dense mafic lower crust emplaced by the rifting processes. The incoming plate shows low uppermost mantle velocities, indicating serpentinization. This hydration is more pronounced in the Shumagin segment, with greater velocity reduction extending to 18 ± 3 km depth, compared to the Semidi segment, showing smaller reductions extending to 14 ± 3 km depth. Our estimates of percent serpentinization from VSreduction and VP/VSare larger than those determined using VPreduction in prior studies, likely due to water in cracks affecting VSmore than VP. Revised estimates of serpentinization show that more water subducts than previous studies, and that twice as much mantle water is subducted in the Shumagin segment compared to the Semidi segment. Together with estimates from other subduction zones, the results indicate a wide variation in subducted mantle water between different subduction segments.more » « less
-
null (Ed.)Abstract Observations of teleseismic earthquakes using broadband seismometers on the Ross Ice Shelf (RIS) must contend with environmental and structural processes that do not exist for land-sited seismometers. Important considerations are: (1) a broadband, multi-mode ambient wavefield excited by ocean gravity wave interactions with the ice shelf; (2) body wave reverberations produced by seismic impedance contrasts at the ice/water and water/seafloor interfaces and (3) decoupling of the solid Earth horizontal wavefield by the sub-shelf water column. We analyze seasonal and geographic variations in signal-to-noise ratios for teleseismic P-wave (0.5–2.0 s), S-wave (10–15 s) and surface wave (13–25 s) arrivals relative to the RIS noise field. We use ice and water layer reverberations generated by teleseismic P-waves to accurately estimate the sub-station thicknesses of these layers. We present observations consistent with the theoretically predicted transition of the water column from compressible to incompressible mechanics, relevant for vertically incident solid Earth waves with periods longer than 3 s. Finally, we observe symmetric-mode Lamb waves generated by teleseismic S-waves incident on the grounding zones. Despite their complexity, we conclude that teleseismic coda can be utilized for passive imaging of sub-shelf Earth structure, although longer deployments relative to conventional land-sited seismometers will be necessary to acquire adequate data.more » « less
