Automatic slowness vector measurements of seismic arrivals with uncertainty estimates using bootstrap sampling, array methods and unsupervised learning

Author:

Ward J1ORCID,Thorne M2,Nowacki A1,Rost S1ORCID

Affiliation:

1. School of Earth and Environment, University of Leeds, Woodhouse, Leeds LS2 9JT, United Kingdom

2. Department of Geology and Geophysics, University of Utah, 115 S 1460 E, Salt Lake City, UT 84112, USA

Abstract

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.

Funder

NERC

NSF

Publisher

Oxford University Press (OUP)

Subject

Geochemistry and Petrology,Geophysics

Reference69 articles.

1. Azimuthal anisotropy from array analysis of Rayleigh waves in Southern California;Alvizuri;Geophys. J. Int.,2011

2. OPTICS: ordering points to identify the clustering structure;Ankerst;ACM . Record,1999

3. Estimation of slowness vectors and their uncertainties using multi-wavelet seismic array processing;Bear;Bull. seism. Soc. Am.,1997

4. Source directionality of ambient seismic noise inferred from three-component beamforming;Behr;J. geophys. Res.,2013

5. Scattering beneath Western Pacific subduction zones: evidence for oceanic crust in the mid-mantle;Bentham;Geophys. J. Int.,2014

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3