scholarly journals Finite frequency inversion of cross-correlation amplitudes for ambient noise source directivity estimation

2019 ◽  
Author(s):  
Arjun Datta ◽  
Shravan Hanasoge ◽  
Jeroen Goudswaard
2020 ◽  
Author(s):  
Alexey Gokhberg ◽  
Laura Ermert ◽  
Jonas Igel ◽  
Andreas Fichtner

<p>The study of ambient seismic noise sources and their time- and space-dependent distribution is becoming a crucial component of the real-time monitoring of various geosystems, including active fault zones and volcanoes, as well as geothermal and hydrocarbon reservoirs. In this context, we have previously implemented a combined cloud - HPC infrastructure for production of ambient source maps with high temporal resolution. It covers the entire European continent and the North Atlantic, and is based on seismic data provided by the ORFEUS infrastructure. The solution is based on the Application-as-a-Service concept and includes (1) acquisition of data from distributed ORFEUS data archives, (2) noise source mapping, (3) workflow management, and (4) front-end Web interface to end users.</p><p>We present the new results of this ongoing project conducted with support of the Swiss National Supercomputing Centre (CSCS). Our recent goal has been transitioning from mapping the seismic noise sources towards modeling them based on our new method for near real-time finite-frequency ambient seismic noise source inversion. To invert for the power spectral density of the noise source distribution of the secondary microseisms we efficiently forward model global cross-correlation wavefields for any noise distribution. Subsequently, a gradient-based iterative inversion method employing finite-frequency sensitivity kernels is implemented to reduce the misfit between synthetic and observed cross correlations.</p><p>During this research we encountered substantial challenges related to the large data volumes and high computational complexity of involved algorithms. We handle these problems by using the CSCS massively parallel heterogeneous supercomputer "Piz Daint". We also apply various specialized numeric techniques which include: (1) using precomputed Green's functions databases generated offline with Axisem and efficiently extracted with Instaseis package and (2) our previously developed high performance package for massive cross correlation of seismograms using GPU accelerators. Furthermore, due to the inherent restrictions of supercomputers, some crucial components of the processing pipeline including the data acquisition and workflow management are deployed on the OpenStack cloud environment. The resulting solution combines the specific advantages of the supercomputer and cloud platforms thus providing a viable distributed platform for the large-scale modeling of seismic noise sources.</p>


2021 ◽  
Author(s):  
Changjiang Zhou ◽  
Jianghai Xia ◽  
Feng Cheng ◽  
Jingyin Pang ◽  
Xinhua Chen

<p>Abundant noise sources in urban area has been widely utilized for subsurface investigations based on the seiemic interferometry. Reliable dispersion extraction between two seismic stations is an essential basis of surface wave imaging. Noise source directivity has become an inescapable obstacle and a main concern for passive seismic surveys: it basically breaks the physics of Green’s function retrieval in travel-time tomography; Moreover, the azimuthal effect of ambient noise sources would inherently cause different levels of early arrival on cross-correlation functions, so that the apparent velocity of surface wave could be overestimated in multichannel slant stackings.</p><p> </p><p>Instead of the conventional frequency-time analysis, which aims to extract the apparent dispersions of phase/group velocity between seismic stations, we proposed a method to jointly invert noise source distributions and the corresponding unbiased surface wave velocities based on the theoretical framework of full waveform ambient noise inversion. Waveform itself could intrinsically contains the features of travel-time, energy and asymmetry of ambient noise cross correlation functions (NCF). And they could in return map the resulted NCF into the noise source distributions and velocity structures. The L2 norm of cross-correlating waveform misfits was taken as the objective function to conduct gradient based inversion (i.e. the L-BFGS algorithm). We parametrized the noise source distributions as a temporally ensemble averaged model, which was discretized as a spatially plane grid of normalized source strength. The surface wave velocity model was approximated as the straight-ray interstation velocity. The two kinds of variants were decoupled in waveform misfit function by their corresponding partial derivatives to iteratively update the model space.</p><p> </p><p>The effectiveness of source-velocity joint imaging using above full waveform inversion work flow was qualified by both the synthetic test and the applied research in Hangzhou urban area. The inverted noise source model was comparable with the urban traffic- and construction- noise distributions. And the truthful surface wave velocities were achieved considering the constraint of noise source distributions, they were also prior constrained and later verified by local borehole datasets.</p>


2020 ◽  
Author(s):  
Jonas Igel ◽  
Laura Ermert ◽  
Andreas Fichtner

<p>Common assumptions in ambient noise seismology such as Green’s function retrieval and equipartitioned wavefields are often not met in the Earth. Full waveform ambient noise tomography methods are free of such assumptions, as they implement knowledge of the time- and space-dependent ambient noise source distribution, whilst also taking finite-frequency effects into account. Such methods would greatly simplify near real-time monitoring of the sub-surface. Additionally, the distribution of the secondary microseisms could act as a new observable of the ocean state since its mechanism is well understood (e.g. Ardhuin et al., 2011).</p><p>To efficiently forward-model global noise cross-correlations we implement (1) pre-computed high-frequency wavefields obtained using, for example, AxiSEM (Nissen-Meyer et al., 2014), and (2) spatially variable grids, both of which greatly reduce the computational cost. Global cross-correlations for any source distribution can be computed within a few seconds in the microseismic frequency range (up to 0.2 Hz). Similarly, we can compute the finite-frequency sensitivity kernels which are then used to perform a gradient-based iterative inversion of the power-spectral density of the noise source distribution. We take a windowed logarithmic energy ratio of the causal and acausal branches of the cross-correlations as measurement, which is largely insensitive to unknown 3D Earth structures.</p><p>Due to its parallelisation on a cluster, our inversion tool is able to rapidly invert for the global microseismic noise source distribution with minimal required user interaction. Synthetic and real data inversions show promising results for noise sources in the North Atlantic with the structure and spatial distribution resolved at scales of a few hundred kilometres. Finally, daily noise sources maps could be computed by combining our inversion tool with a daily data download and processing toolkit.</p>


2020 ◽  
Author(s):  
Laura Ermert ◽  
Jonas Igel ◽  
Korbinian Sager ◽  
Eléonore Stutzmann ◽  
Tarje Nissen-Meyer ◽  
...  

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.


Sign in / Sign up

Export Citation Format

Share Document