An official website of the United States government
Abstract Continent‐scale observations of seismic phenomena have provided multi‐scale constraints of the Earth's interior. Of those analyzed, array‐based observations of slowness vector properties (backazimuth and horizontal slowness) and multipathing have yet to be made on a continental scale. Slowness vector measurements give inferences on mantle heterogeneity properties such as velocity perturbation and velocity gradient strength and quantify their effect on the wavefield. Multipathing is a consequence of waves interacting with strong velocity gradients resulting in two arrivals with different slowness vector properties and times. The mantle structure beneath the contiguous Unites States has been thoroughly analyzed by previous seismic studies and is data‐rich, making it an excellent testing ground to both analyze mantle structure with our approach and compare with other imaging techniques. We apply an automated array‐analysis technique to an SKS data set to create the first continent‐scale data set of multipathing and slowness vector measurements. We analyze the divergence of the slowness vector deviation field to highlight seismically slow and fast regions. Our results resolve several slow mantle anomalies beneath Yellowstone, the Appalachian mountains and fast anomalies throughout the mantle. Many of the anomalies cause multipathing in frequency bands 0.15–0.30 and 0.20–0.40 Hz which suggests velocity transitions over at most 500 km exist. Comparing our observations to synthetics created from tomography models, we find model NA13 (Bedle et al., 2021,https://doi.org/10.1029/2021GC009674) fits our data best but differences still remain. We therefore suggest slowness vector measurements should be used as an additional constraint in tomographic inversions and will lead to better resolved models of the mantle.
more »
« less
Automatic slowness vector measurements of seismic arrivals with uncertainty estimates using bootstrap sampling, array methods and unsupervised learning
SUMMARY Horizontal slowness vector measurements using array techniques have been used to analyse many Earth phenomena from lower mantle heterogeneity to meteorological event location. While providing observations essential for studying much of the Earth, slowness vector analysis is limited by the necessary and subjective visual inspection of observations. Furthermore, it is challenging to determine the uncertainties caused by limitations of array processing such as array geometry, local structure, noise and their effect on slowness vector measurements. To address these issues, we present a method to automatically identify seismic arrivals and measure their slowness vector properties with uncertainty bounds. We do this by bootstrap sampling waveforms, therefore also creating random sub arrays, then use linear beamforming to measure the coherent power at a range of slowness vectors. For each bootstrap sample, we take the top N peaks from each power distribution as the slowness vectors of possible arrivals. The slowness vectors of all bootstrap samples are gathered and the clustering algorithm DBSCAN (Density-Based Spatial Clustering of Applications with Noise) is used to identify arrivals as clusters of slowness vectors. The mean of slowness vectors in each cluster gives the slowness vector measurement for that arrival and the distribution of slowness vectors in each cluster gives the uncertainty estimate. We tuned the parameters of DBSCAN using a data set of 2489 SKS and SKKS observations at a range of frequency bands from 0.1 to 1 Hz. We then present examples at higher frequencies (0.5–2.0 Hz) than the tuning data set, identifying PKP precursors, and lower frequency by identifying multipathing in surface waves (0.04–0.06 Hz). While we use a linear beamforming process, this method can be implemented with any beamforming process such as cross correlation beamforming or phase weighted stacking. This method allows for much larger data sets to be analysed without visual inspection of data. Phenomena such as multipathing, reflections or scattering can be identified automatically in body or surface waves and their properties analysed with uncertainties.
more »
« less
- Award ID(s):
- 1723081
- PAR ID:
- 10262163
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 226
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 1847 to 1857
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Seismic array processing is routinely used to infer detailed earthquakeproperties of intermediate and large events, however, the sourceproperties of microseismicity often remain elusive. In this study, weuse high signal-to-noise ratio seismograms of 204 earthquakes induced bythe 6 km deep 2018 Espoo/Helsinki geothermal stimulation to evaluate thecapabilities of beamforming and back-projection array methods. We showthat mini array beamforming is sensitive to medium heterogeneity andrequires calibration to mitigate systematic slowness biases.A combinedand wave back-projection approach significantly improves depthresolution, reducing offsets to catalogue locations from km to m.Supported by numerical experiments, we demonstrate that back-projectionswimming patterns can constrain focal mechanisms. Our results imply thatback-projection of data collected over a wide azimuthal range can beused to monitor and characterize local-scale microseismicity, whereasbeamforming calibration requires independently obtained referenceobservations.more » « less
-
Clustering is a fundamental task in machine learning. One of the most successful and broadly used algorithms is DBSCAN, a density-based clustering algorithm. DBSCAN requires ϵ-nearest neighbor graphs of the input dataset, which are computed with range-search algorithms and spatial data structures like KD-trees. Despite many efforts to design scalable implementations for DBSCAN, existing work is limited to low-dimensional datasets, as constructing ϵ-nearest neighbor graphs can be expensive in high-dimensions. This article introduces a modified DBSCAN, usingk-nearest neighbor (kNN) graphs to improve efficiency. We outline conditions forkNN-DBSCAN to match DBSCAN’s results and present a parallel implementation using OpenMP and MPI for shared and distributed memory systems. Testing on datasets up to 32 dimensions, we achieve remarkable scalability. Our implementation clusters one billion 3D points in under one second on 28K cores at TACC’s Frontera system. In a larger run, we cluster 65 billion points in 20 dimensions in under 40 seconds using 114,688 cores. Our method is up to 37× faster than state-of-the-art parallel DBSCAN on a 20-dimensional dataset with 4 million points. Code is available athttps://github.com/ut-padas/knndbscan.more » « less
-
Abstract Seismic arrays constrain local wave propagation that can be used to infer earthquake source characteristics. Array processing is routinely used to infer detailed earthquake properties of intermediate and large events. However, the source properties of microseismicity often remain elusive. In this study, we use high signal-to-noise ratio seismograms of 204 ML 0.0–1.8 earthquakes induced by the 6 km deep 2018 Espoo/Helsinki geothermal stimulation to evaluate the performance and capabilities of beamforming and backprojection array methods. Using accurate travel-time-based event locations as a reference, we first show that miniarray beamforming is sensitive to medium heterogeneities and requires calibration to mitigate local systematic slowness biases. A catalog-based calibration significantly improves our multiarray beam raytracing estimates of source locations. Second, the application of the backprojection technique using P-wave signals with sufficient azimuthal coverage yields hypocenter estimates with generally good horizontal but poor vertical resolution. The short local source–receiver distances result in incomplete separation of P- and S-wave arrivals during backprojection. Numerical tests show that the relatively large S-wave amplitudes can influence coherent P-wave stacks, resulting in large location errors. Our combined P- and S-wave backprojection approach mitigates the influence of the large S-wave amplitude and improves the depth resolution significantly. The average depth offset to the reference catalog locations reduces from ≥1.4 km to ∼91 m. Third, 3D numerical simulations demonstrate that backprojection swimming patterns are not merely processing or configuration artifacts. We show that the swimming patterns correlate with and can resolve the source focal mechanism when the azimuthal wavefield sampling is sufficiently complete. Our work demonstrates that the backprojection techniques can help to better constrain important properties of local-scale microseismicity.more » « less
-
Abstract Shear‐wave splitting measurements are commonly used to resolve seismic anisotropy in both the upper and lowermost mantle. Typically, such techniques are applied to SmKS phases that have reflected (m‐1) times off the underside of the core‐mantle boundary before being recorded. Practical constraints for shear‐wave splitting studies include the limited number of suitable phases as well as the large fraction of available data discarded because of poor signal‐to‐noise ratios (SNRs) or large measurement uncertainties. Array techniques such as beamforming are commonly used in observational seismology to enhance SNRs, but have not been applied before to improve SmKS signal strength and coherency for shear wave splitting studies. Here, we investigate how a beamforming methodology, based on slowness and backazimuth vespagrams to determine the most coherent incoming wave direction, can improve shear‐wave splitting measurement confidence intervals. Through the analysis of real and synthetic seismograms, we show that (a) the splitting measurements obtained from the beamformed seismograms (beams) reflect an average of the single‐station splitting parameters that contribute to the beam; (b) the beams have (on average) more than twice as large SNRs than the single‐station seismograms that contribute to the beam; (c) the increased SNRs allow the reliable measurement of shear wave splitting parameters from beams down to average single‐station SNRs of 1.3. Beamforming may thus be helpful to more reliably measure splitting due to upper mantle anisotropy. Moreover, we show that beamforming holds potential to greatly improve detection of lowermost mantle anisotropy by demonstrating differential SKS–SKKS splitting analysis using beamformed USArray data.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggab196
