skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Rayleigh-wave multicomponent crosscorrelation-based source strength distribution inversions. Part 2: a workflow for field seismic data
SUMMARY Estimation of ambient seismic source distributions (e.g. location and strength) can aid studies of seismic source mechanisms and subsurface structure investigations. One can invert for the ambient seismic (noise) source distribution by applying full-waveform inversion (FWI) theory to seismic (noise) crosscorrelations. This estimation method is especially applicable for seismic recordings without obvious body-wave arrivals. Data pre-processing procedures are needed before the inversion, but some pre-processing procedures commonly used in ambient noise tomography can bias the ambient (noise) source distribution estimation and should not be used in FWI. Taking this into account, we propose a complete workflow from the raw seismic noise recording through pre-processing procedures to the inversion. We present the workflow with a field data example in Hartoušov, Czech Republic, where the seismic sources are CO2 degassing areas at Earth’s surface (i.e. a fumarole or mofette). We discuss factors in the processing and inversion that can bias the estimations, such as inaccurate velocity model, anelasticity and array sensitivity. The proposed workflow can work for multicomponent data across different scales of field data.  more » « less
Award ID(s):
1643795
PAR ID:
10203459
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
222
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
2084 to 2101
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Elastic full-waveform inversion (EFWI) is a state-of-the-art seismic tomographic method. Recent advances in technology and instrumentation, combining crosstalk-free source-encoded FWI (SE-FWI) with multicomponent marine data acquisition using ocean-bottom nodes (OBNs), enable full-physics wave propagation and parameter inversion without the computational burden of traditional FWI. With OBN acquisition, P waves, S waves, and P-to-S conversions are recorded. It is not well understood to what extent adding horizontal components to SE-FWI improves the resolution of subsurface modeling. We assess their potential for the reconstruction of shear and compressional wave speeds (VPand VS) by using a synthetic data set modeled after a recently acquired OBN survey in the North Sea. We perform synthetic inversion tests to design suitable strategies that leverage the information recorded in the horizontal components of the data to improve the reconstructed model resolution laterally and in depth. We advocate for a hierarchical inversion approach to recover the elastic parameters. We exploit the P and P-to-S converted waves recorded on the horizontal components to robustly reconstruct both VPand VS. Adding horizontal components to the SE-FWI modeling workflow results in improved spatial resolution, enhanced depth coverage, and more accurate elastic wave speed estimates. 
    more » « less
  3. Understanding subsurface conditions is critical to creating and maintaining resilient infrastructure systems, such as dams and levees. Seismic geophysical tools can be very effective for site characterization of these structures as they directly measure the elastic moduli and can provide insight into both the soil properties and groundwater conditions. Full waveform inversion (FWI) is one processing option for seismic geophysics that seeks to overcome some of the limitations in the traditional approaches by using the full time-domain recording of the wavefield to develop 2D or 3D profiles of shear wave velocity. In addition to providing characterization data, FWI can also potentially be used as a monitoring tool for dams and levees to assess how elastic moduli are changing with time and to infer how these changes might relate to changes in the hydromechanical properties of the soil. This study seeks to explore the use of seismic FWI as both a characterization and monitoring tool through numerical simulations of seismic surveys on a hypothetical levee with a low velocity anomaly in the foundation. The simulations are used to assess both the spatial resolution and the ability of the simulations to detect changes in properties that might be related to softening of the foundation or development of internal erosion failure modes. The findings from the study will be used to highlight potential benefits and challenges to using seismic FWI for characterization and monitoring of dams and levees. 
    more » « less
  4. We exploit the different but complementary data sensitivities of ground-penetrating radar (GPR) and electrical resistivity (ER) by applying a multiphysics, multiparameter, simultaneous 2.5D joint inversion without invoking petrophysical relationships. Our method joins full-waveform inversion (FWI) GPR with adjoint derived ER sensitivities on the same computational domain. We incorporate a stable source estimation routine into the FWI-GPR. We apply our method in a controlled alluvial aquifer using only surface-acquired data. The site exhibits a shallow groundwater boundary and unconsolidated heterogeneous alluvial deposits. We compare our recovered parameters to individual FWI-GPR and ER results, and we compare them to log measurements of capacitive conductivity and neutron-derived porosity. Our joint inversion provides a more representative depiction of subsurface structures because it incorporates multiple intrinsic parameters, and it is therefore superior to an interpretation based on log data, FWI-GPR, or ER alone. 
    more » « less
  5. SUMMARY The uneven distribution of earthquakes and stations in seismic tomography leads to slower convergence of nonlinear inversions and spatial bias in inversion results. Including dense regional arrays, such as USArray or Hi-Net, in global tomography causes severe convergence and spatial bias problems, against which conventional pre-conditioning schemes are ineffective. To save computational cost and reduce model bias, we propose a new strategy based on a geographical weighting of sources and receivers. Unlike approaches based on ray density or the Voronoi tessellation, this method scales to large full-waveform inversion problems and avoids instabilities at the edges of dense receiver or source clusters. We validate our strategy using a 2-D global waveform inversion test and show that the new weighting scheme leads to a nearly twofold reduction in model error and much faster convergence relative to a conventionally pre-conditioned inversion. We implement this geographical weighting strategy for global adjoint tomography. 
    more » « less