Ambient noise data selection based on the asymmetry of cross-correlation functions for near surface applications

2018 ◽  
Vol 159 ◽  
pp. 803-813 ◽  
Author(s):  
Changjiang Zhou ◽  
Chaoqiang Xi ◽  
Jingyin Pang ◽  
Ya Liu
Author(s):  
Jie Zhou ◽  
Xiaofei Chen

ABSTRACT Frequency–Bessel (F-J) transform method can obtain higher-mode Rayleigh dispersion curves by multistation ambient noise data superposition (Wang et al., 2019). Because the dispersion curves of the overtones can provide more information compared with the single fundamental mode, the nonuniqueness of surface-wave inversion can be reduced. Because of the limited number of receivers, the integral in the process of transformation cannot be calculated precisely and there exists a kind of crossed artifacts which cuts off the real dispersion curves and contaminates the spectrum. Forbriger (2003) proposed to use the Hankel function instead of the Bessel function to conduct the transformation to remove the crossed artifacts. However, this method can reduce the resolution of the spectrum from ambient noise data. In this article, we give a complete workflow to deal with ambient noises which can eliminate the crossed artifacts without reducing the resolution. The Kramers–Kronig relations are used to obtain complete cross-correlation functions and a modified F-J transform is conducted to finally acquire the spectrum without crossed artifacts.


2021 ◽  
Author(s):  
Martha Savage ◽  
FC Lin ◽  
John Townend

Measurement of basement seismic resonance frequencies can elucidate shallow velocity structure, an important factor in earthquake hazard estimation. Ambient noise cross correlation, which is well-suited to studying shallow earth structure, is commonly used to analyze fundamental-mode Rayleigh waves and, increasingly, Love waves. Here we show via multicomponent ambient noise cross correlation that the basement resonance frequency in the Canterbury region of New Zealand can be straightforwardly determined based on the horizontal to vertical amplitude ratio (H/V ratio) of the first higher-mode Rayleigh waves. At periods of 1-3 s, the first higher-mode is evident on the radial-radial cross-correlation functions but almost absent in the vertical-vertical cross-correlation functions, implying longitudinal motion and a high H/V ratio. A one-dimensional regional velocity model incorporating a ~ 1.5 km-thick sedimentary layer fits both the observed H/V ratio and Rayleigh wave group velocity. Similar analysis may enable resonance characteristics of other sedimentary basins to be determined. © 2013. American Geophysical Union. All Rights Reserved.


2014 ◽  
Vol 136 (4) ◽  
pp. 2156-2156
Author(s):  
Xiaoqin Zang ◽  
Michael G. Brown ◽  
Neil J. Williams ◽  
Oleg A. Godin ◽  
Nikolay A. Zabotin ◽  
...  

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. F1-F8
Author(s):  
Eileen R. Martin

Geoscientists and engineers are increasingly using denser arrays for continuous seismic monitoring, and they often turn to ambient seismic noise interferometry for low-cost near-surface imaging. Although ambient noise interferometry greatly reduces acquisition costs, the computational cost of pair-wise comparisons between all sensors can be prohibitively slow or expensive for applications in engineering and environmental geophysics. Double beamforming of noise correlation functions is a powerful technique to extract body waves from ambient noise, but it is typically performed via pair-wise comparisons between all sensors in two dense array patches (scaling as the product of the number of sensors in one patch with the number of sensors in the other patch). By rearranging the operations involved in the double beamforming transform, I have developed a new algorithm that scales as the sum of the number of sensors in two array patches. Compared to traditional double beamforming of noise correlation functions, the new method is more scalable, easily parallelized, and it does not require raw data to be exchanged between dense array patches.


2021 ◽  
Author(s):  
Takashi Hirose ◽  
Hideki Ueda ◽  
Eisuke Fujita

<p>    Estimating seismic scattering and intrinsic absorption parameters, which are measures of medium heterogeneity, is important for understanding the complex structure in shallow regions of volcanoes. In recent years, seismic ambient noise cross-correlation functions (CCFs) have been used instead of records of natural earthquakes or active seismic experiments to estimate those parameters (e.g., Hirose et al., 2019; Hirose et al., 2020; van Dinther et al., 2020). This passive approach possibly allows us to estimate scattering and intrinsic absorption parameters in previously unmeasured regions and frequency bands. In this study, we apply the passive estimation method proposed by Hirose et al. (2019) to 18 active volcanoes in Japan and measure those parameters of Rayleigh waves. We used three-component seismic ambient noise data in the frequency bands of 0.5-1 Hz, 1-2 Hz, and 2-4 Hz at seismic stations of NIED, JMA, HSRI, and MFRI. Before computing CCFs, the temporal flattening technique (Weaver, 2011) was applied to ambient noise data for reducing the effect of temporal fluctuations in noise levels with retaining relative amplitudes among the stations. Daily CCFs of three components (ZZ, ZR, ZT) were computed by stacking 10-minutes-CCFs. We stacked daily CCFs over 1 year and computed mean squared envelopes by smoothing squared amplitude with 4 s (0.5-1 Hz), 2 s (1-2 Hz), or 1 s (2-4 Hz) long time windows. Scattering and intrinsic absorption parameters were estimated by modeling the space-time distributions of energy densities calculated from CCFs with 2D radiative transfer theory. Best-fit values of scattering mean free path at the 18 active volcanoes range between 1.0-4.6 km at 0.5-1Hz band, 0.7-2.9 km at 1-2 Hz band, and 0.9-2.9 km at 2-4 Hz band, respectively. These values are 2 orders of magnitude shorter than those in non-volcanic regions (e.g., Sato et al., 2012). Those of intrinsic absorption parameter range between 0.05-0.26 s<sup>-1</sup> at the 0.5-1 Hz band, 0.06-0.24 s<sup>-1</sup> at the 1-2 Hz band, and 0.06-0.32 s<sup>-1 </sup>at the 2-4 Hz band, respectively. They are at most one order of magnitude larger than those in the non-volcanic regions. Especially strong intrinsic attenuations are estimated at volcanic islands. Water-bearing layers at a depth of several hundred meters below these islands may cause such strong intrinsic attenuations. The frequency dependence of scattering attenuations is also strong at these volcanic islands, suggesting non-uniform structures that largely fluctuate along depths. The results of this study suggest that the passive estimation method of scattering and intrinsic absorption parameters proposed by Hirose et al. (2019) is applicable to various volcanoes. Comparing estimated values of these parameters at various volcanoes will improve our understanding of complex structure at the shallow regions of volcanoes. Moreover, the parameters estimated in this study will boost locating spatial distributions of seismic velocity and/or scattering property changes associated with volcanic activities at the 18 volcanoes.</p><p>Acknowledgments: We used seismograms recorded by Japan Meteorological Agency (JMA), Hot Springs Research Institute (HSRI) of Kanagawa Prefecture, and Mount Fuji Research Institute (MFRI), Yamanashi Prefectural Government.</p>


2019 ◽  
Vol 220 (3) ◽  
pp. 1521-1535
Author(s):  
Loïc Viens ◽  
Chris Van Houtte

SUMMARY Seismic interferometry is an established method for monitoring the temporal evolution of the Earth’s physical properties. We introduce a new technique to improve the precision and temporal resolution of seismic monitoring studies based on deep learning. Our method uses a convolutional denoising autoencoder, called ConvDeNoise, to denoise ambient seismic field correlation functions. The technique can be applied to traditional two-station cross-correlation functions but this study focuses on single-station cross-correlation (SC) functions. SC functions are computed by cross correlating the different components of a single seismic station and can be used to monitor the temporal evolution of the Earth’s near surface. We train and apply our algorithm to SC functions computed with a time resolution of 20 min at seismic stations in the Tokyo metropolitan area, Japan. We show that the relative seismic velocity change [dv/v(t)] computed from SC functions denoised with ConvDeNoise has less variability than that calculated from raw SC functions. Compared to other denoising methods such as the SVD-based Wiener Filter method developed by Moreau et al., the dv/v results obtained after using our algorithm have similar precision. The advantage of our technique is that once the algorithm is trained, it can be apply to denoise near-real-time SC functions. The near-real-time aspect of our denoising algorithm may be useful for operational hazard forecasting models, for example when applying seismic interferometry at an active volcano.


2021 ◽  
Author(s):  
◽  
Andy McNab

<p>This thesis applies ambient noise tomography to investigate the shallow structure of the Whataroa Valley. Ambient noise techniques are applied to continuous seismic recordings acquired on 158 geophones deployed during the Whataroa Active Source Seismic Experiment. Despite only having four days of data, a robust shear-wave velocity model is calculated using a phase-weighted stacking approach to improve the cross-correlation functions' signal-to-noise ratios, allowing for robust velocity measurements to be obtained between periods of 0.3 and 1.8\,s. This yields a database of 12,500 vertical component cross correlation functions and the corresponding Rayleigh wave phase and group velocity dispersion curves. Linearised straight-ray tomography is applied to phase and group velocity dispersion measurements at periods ranging from periods of 0.3 to 1.8\,s. The tomography reveals a velocity that decreases from the vicinity of the DFDP-2B borehole to the centre of the valley. This is interpreted to be the geologic basement deepening towards the centre of the valley. A Monte-Carlo inversion technique is used to jointly invert Rayleigh-wave phase and group velocity dispersion curves constructed from phase and group velocity tomography maps of successively higher periods. Linear interpolation of the resulting 1D shear-wave velocity profiles produces a pseudo-3D velocity model of the uppermost 1,000\,m of the Whataroa Valley. Using sharp increases in velocity to represent lithological change, we interpret two velocity contours at 1,150 and 1,250\,m/s as potential sediment-basement contacts. Depth isocontours of these velocities reveal that the basement deepens towards the centre of the valley, reaching a maximum depth of 400 or 600\,m for the 1,150 and 1,250\,m/s velocity contours respectively. These depths indicate strong glacial over-deepening and have implications for future drilling projects in the Whataroa Valley. A sharp velocity increase of 200\,m/s also occurs at 100\,m depth at the DFDP-2B borehole. We interpret this to be a change in sedimentary rock lithology from fluvial gravels to lacustrine silty sands, related to a change in sedimentary depositional environment.</p>


Sign in / Sign up

Export Citation Format

Share Document