Cross-correlations of ambient seismic noise are widely used for seismic velocity imaging, monitoring and ground motion analyses. A typical step in analysing noise cross-correlation functions (NCFs) is stacking short-term NCFs over longer time periods to increase the signal quality. Spurious NCFs could contaminate the stack, degrade its quality and limit its use. Many methods have been developed to improve the stacking of coherent waveforms, including earthquake waveforms, receiver functions and NCFs. This study systematically evaluates and compares the performance of eight stacking methods, including arithmetic mean or linear stacking, robust stacking, selective stacking, cluster stacking, phase-weighted stacking, time–frequency phase-weighted stacking, Nth-root stacking and averaging after applying an adaptive covariance filter. Our results demonstrate that, in most cases, all methods can retrieve clear ballistic or first arrivals. However, they yield significant differences in preserving the phase and amplitude information. This study provides a practical guide for choosing the optimal stacking method for specific research applications in ambient noise seismology. We evaluate the performance using multiple onshore and offshore seismic arrays in the Pacific Northwest region. We compare these stacking methods for NCFs calculated from raw ambient noise (referred to as Raw NCFs) and from ambient noise normalized using a one-bit clipping time normalization method (referred to as One-bit NCFs). We evaluate six metrics, including signal-to-noise ratios, phase dispersion images, convergence rate, temporal changes in the ballistic and coda waves, relative amplitude decays with distance and computational time. We show that robust stacking is the best choice for all applications (velocity tomography, monitoring and attenuation studies) using Raw NCFs. For applications using One-bit NCFs, all methods but phase-weighted and Nth-root stacking are good choices for seismic velocity tomography. Linear, robust and selective stacking methods are all equally appropriate choices when using One-bit NCFs for monitoring applications. For applications relying on accurate relative amplitudes, the linear, robust, selective and cluster stacking methods all perform well with One-bit NCFs. The evaluations in this study can be generalized to a broad range of time-series analysis that utilizes data coherence to perform ensemble stacking. Another contribution of this study is the accompanying open-source software package, StackMaster, which can be used for general purposes of time-series stacking.
This content will become publicly available on August 7, 2024
- Award ID(s):
- 2003619
- NSF-PAR ID:
- 10440192
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 5
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
more » « less
Abstract Typical use of ambient noise interferometry focuses on longer period (>1 s) waves for exploration of subsurface structure and other applications, while very shallow structure and some environmental seismology applications may benefit from use of shorter period (<1 s) waves. We explore the potential for short‐period ambient noise interferometry to determine shallow seismic velocity structures by comparing two methodologies, the conventional amplitude‐based cross‐correlation and linear stacking (TCC‐Lin) and a more recently developed phase cross‐correlation and time‐frequency phase‐weighted‐stacking (PCC‐PWS) method with both synthetic and real data collected in a heterogeneous karst aquifer system. Our results suggest that the PCC‐PWS method is more effective in extracting short‐period wave velocities than the TCC‐Lin method, especially when using data collected in regions containing complex shallow structures such as the karst aquifer system investigated here. In addition to the different methodologies for computing the cross correlation functions, we also examine the relative importance of signal‐to‐noise ratio and number of wavelengths propagating between station pairs to determine data/solution quality. We find that the lower number of wavelengths of 3 has the greatest impact on the network‐averaged group velocity curve. Lastly, we test the sensitivity of the number of stacks used to create the final empirical Green's function, and find that the PCC‐PWS method required about half the number of cross‐correlation functions to develop reliable velocity curves compared to the TCC‐Lin method. This is an important advantage of the PCC‐PWS method when available data collection time is limited.
-
Abstract The CO Mapping Array Project (COMAP) aims to use line-intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1–0) emission from z = 2.4–3.4 and a fainter contribution from CO(2–1) at z = 6–8, the Pathfinder is surveying 12 deg 2 in a 5 yr observing campaign to detect the CO signal from z ∼ 3. Using data from the first 13 months of observing, we estimate P CO ( k ) = −2.7 ± 1.7 × 10 4 μ K 2 Mpc 3 on scales k = 0.051 −0.62 Mpc −1 , the first direct three-dimensional constraint on the clustering component of the CO(1–0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature bias product) of Tb 2 < 49 μ K 2 , nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 yr with signal-to-noise ratio (S/N) 9–17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO–galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR.more » « less
-
more » « less
Seismic imaging and monitoring of the near-surface structure are crucial for the sustainable development of urban areas. However, standard seismic surveys based on cabled or autonomous geophone arrays are expensive and hard to adapt to noisy metropolitan environments. Distributed acoustic sensing (DAS) with pre-existing telecom fiber optic cables, together with seismic ambient noise interferometry, have the potential to fulfill this gap. However, a detailed noise wavefield characterization is needed before retrievingcoherent waves from chaotic noise sources. We analyze local seismic ambient noise by tracking five-month changes in signal-to-noise ratio (SNR) of Rayleigh surface wave estimated from traffic noise recorded by DAS along the straight university campus busy road. We apply the seismic interferometry method to the 800 m long part of the Penn State Fiber-Optic For Environment Sensing (FORESEE) array. We evaluate the 160 virtual shot gathers (VSGs) by determining the SNR using the slant-stack technique. We observe strong SNR variations in time and space. We notice higher SNR for virtual source points close to road obstacles. The spatial noise distribution confirms that noise energy focuses mainly on bumps and utility holes. We also see the destructive impact of precipitation, pedestrian traffic, and traffic along main intersections on VSGs. A similar processing workflow can be applied to various straight roadside fiber optic arrays in metropolitan areas.
-
more » « less
SUMMARY Ambient noise tomography is a well-established tomographic imaging technique but the effect that spatially variable noise sources have on the measurements remains challenging to account for. Full waveform ambient noise inversion has emerged recently as a promising solution but is computationally challenging since even distant noise sources can have an influence on the interstation correlation functions and therefore requires a prohibitively large numerical domain, beyond that of the tomographic region of interest. We investigate a new strategy that allows us to reduce the simulation domain while still being able to account for distant contributions. To allow nearby numerical sources to account for distant true sources, we introduce correlated sources and generate a time-dependent effective source distribution at the boundary of a small region of interest that excites the correlation wavefield of a larger domain. In a series of 2-D numerical simulations, we demonstrate that the proposed methodology with correlated sources is able to successfully represent a far-field source that is simultaneously present with nearby sources and the methodology also successfully results in a robustly estimated noise source distribution. Furthermore, we show how beamforming results can be used as prior information regarding the azimuthal variation of the ambient noise sources in helping determine the far-field noise distribution. These experiments provide insight into how to reduce the computational cost needed to perform full waveform ambient noise inversion, which is key to turning it into a viable tomographic technique. In addition, the presented experiments may help reduce source-induced bias in time-dependent monitoring applications.
