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.
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.
more » « less- NSF-PAR ID:
- 10379491
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 232
- Issue:
- 3
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 1600-1618
- Size(s):
- ["p. 1600-1618"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography.
This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography.
Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function.
The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.
-
null (Ed.)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
-
The structure of the Antarctic crust is important to our understanding of processes occurring within the Antarctic cryosphere as well as to the Earth’s response to ice mass loss. With the increase in geophysical studies of Antarctica, crustal structure has become much better defined beneath many regions. Several crustal models have been created from seismic-derived and/or gravity-derived data, and some of these models incorporate sets of crustal receiver functions either as a priori constraints or to validate model results. However, receiver function constraints do not exist throughout large regions of Antarctica due to a lack of seismic coverage; given this, we search for additional metrics by which we can compare and contrast Earth models. One approach that has been utilized for other continents is to forward model accurate synthetic waveforms through existing seismic velocity models to identify which models most accurately reproduce seismic waveform datasets. Such waveform datasets may come from accurately determined seismic events or from ambient seismic noise. In an effort to assess existing Antarctic crustal models using a different metric to identify regions where crustal structure is still most uncertain, we have collected a suite of available seismic- and gravity-derived Antarctic crustal models. In the absence of accurately determined ‘ground-truthed’ seismic events in Antarctica, we use a frequency-time normalization approach to extract Rayleigh waves from ambient seismic noise, with periods of 15-55 seconds that are sensitive to crustal structure. We split the observations into two separate validation datasets. The first dataset includes all station-station cross-correlations, with at least one seismic station in each pair that has not been previously used to constrain prior tomographic inversions (a true validation dataset), and the second dataset includes all available station-station cross-correlations, including those that may have been used to constrain some of the models we are testing. We construct sets of Earth models from the available crustal models underlain by two different upper mantle models. We forward model synthetic waveforms using a finite difference approach through each of the Earth models and measure the phase delays between the synthetic waveforms and the ambient seismic noise dataset. Results from our waveform validation study and identification of the poorly characterized regions of Antarctic crust are forthcoming and will be presented.more » « less
-
null (Ed.)The origins of tectonic structures in East Antarctica, such as the Gamburtsev Subglacial Mountains (GSMs), the Wilkes Subglacial Basin (WSB), the Aurora Subglacial Basin (ASB), and the Transantarctic Mountains (TAMs), are not clearly understood. Previous investigations have proposed multiple origin models to explain the formation of these structures; however, existing tomographic images lack resolution and consistency given the sparse seismic coverage in East Antarctica. We use full-waveform ambient noise tomography to model the shear-wave velocity structure beneath East Antarctica to further investigate these features. We extract Rayleigh-wave Empirical Green’s Func-tions (EGFs) between periods of 15 and 340 secs from ambient seismic noise using a frequency time normalization technique. Synthetic waveforms are simulated through a 3-D heterogenous Earth model with a lateral grid spacing of 0.025º (~2.25 km) using a finite-difference wave propagation method. The synthetic seismograms are cross-correlated with the EGFs to measure the phase delays. The fi-nite-frequency sensitivity kernels are calculated using the scattering-integral approach and the shear-wave velocity model is computed by inverting the phase delays using a sparse damped least-square inversion method. Preliminary results show fast seismic velocities beneath the WSB, which may be associated with thick and stable lithosphere, and slow velocities beneath the ASB, possibly reflecting old rift systems or other inherited tectonic structures. Slow upper mantle velocities are also observed beneath the TAMs, possibly associated with a thermal load that contributes to the uplift of the moun-tain range. Slow shear-wave velocities in the vicinity of the GSMs may be associated with rifting along the extended Lambert Rift System. Our final tomographic model and associated tectonic inter-pretations will be shared.more » « less
