Passive seismic methods in highly populated urban areas have gained much attention from geophysics and civil engineering communities because traditional seismic surveys, especially in complex urbanized environments, might be improperly applied. In passive seismic methods, directional noise sources will inevitably bring azimuthal effects and spatial aliasing to dispersion measurements due to the fact that true randomness of ambient noise cannot be achieved in reality. To solve these problems, multichannel analysis of passive surface (MAPS) waves based on long noise sequence crosscorrelations is proposed. We have introduced a hybrid method of seismic interferometry and the roadside passive multichannel analysis of surface waves (MASW) using crosscorrelation to produce common virtual source gathers from 1 h multichannel noise records. Common virtual source gathers are then used to do dispersion analysis with an active scheme based on phase-shift measurement. Synthetic tests demonstrated the advantages of this method with azimuthal adjustment and dispersion imaging for directional noise source distribution. Two field applications were conducted, and results from the roadside passive MASW, MAPS, and spatial autocorrelation method were compared. Our study indicated the superiority of MAPS over the roadside passive MASW on the validity of azimuth detection, feasibility of combining the active MASW and MAPS, and accuracy in determining dispersion energy trends, especially at a relative low-frequency range ($1\mathrm{Hz}\le f\le 10\mathrm{Hz}$) in urban areas.

The multichannel analysis of surface wave (MASW) method has been effectively and widely used to determine near-surface shear-wave velocity. Horizontal resolution of the MASW method represents the minimum horizontal length of recognizable geologic anomalous bodies on a pseudo-2D S-wave velocity [ ] section. Accurately assessing the achievable lateral resolution is one of the main issues in lateral variation reconstruction using the MASW method. It is difficult to quantitatively estimate the horizontal resolution of the MASW method because of the many influencing factors, such as parameters of the observation system, the depth of an anomalous body, and the velocity contrast between the anomalous body and the surrounding rocks. We first analyzed the horizontal resolution of the MASW method based on numerical simulation experiments. According to different influencing factors of the horizontal resolution, we established different laterally heterogeneous models and observation systems and then simulated several synthetic multichannel records with a finite-difference method along a linear survey line using the roll-along acquisition mode. After the extraction of dispersion curves of Rayleigh waves and inversion for S-wave velocity profiles for each synthetic shot gather, a pseudo-2D S-wave velocity section can be generated by aligning the 1D S-wave velocity models. Ultimately, we evaluated the horizontal resolution capability of the MASW method on pseudo-2D [ ] maps. Our numerical investigation results and field data analysis indicate that [ ] values on the maps are not the same as the true [ ] values for structures whose lateral dimension is shorter than a receiver spread length and that anomalous bodies, which are larger and have high velocity contrast, are easier to distinguish on [ ] maps with a shorter receiver spread length. The horizontal resolution decreases with the increasing depth and is approximately one-half of the shortest Rayleigh wavelength that can penetrate to the depth.

SUMMARY: Linear arrays are usually deployed for passive surface-wave investigations because of their high efficiency and convenience. In populated urban areas, it is almost impossible to set up a 2-D array in terms of the restriction from the existing infrastructures. The limited azimuthal coverage, however, lacks the ability to attenuate velocity overestimation caused by directional noise sources. We came up with a novel idea to compensate the azimuthal coverage by adding two more offline receivers to a conventional linear array, which is called pseudo-linear-array analysis of passive surface waves (PLAS). We used a beamforming algorithm to capture noise sources distribution and extract accurate dispersion curves. We used array response function to explain the superiority of the pseudo-linear array over the linear array and present the basic workflow of PLAS. Synthetic tests and field examples demonstrated the feasibility of PLAS to measure unbiased dispersion image. Comparison with mostly used passive surface wave methods (refraction microtremor, multichannel analysis of passive surface waves, spatial autocorrelation method, frequency–wavenumber analysis) suggested that PLAS can serve as an alternative passive surface wave method, especially in urban areas with restricted land accessibility and short-time acquisition demands.

Tomographic imaging based on long-term ambient seismic noise measurements, mainly the phase information from surface waves, has been shown to be a powerful tool for geothermal reservoir imaging and monitoring. In this study, we utilize seismic noise data from a dense nodal array (192 3C nodes within 20 ) over a ultra-short observation period (4.7 days) to reconstruct surface waves and determine the high-resolution (0.2 km) three-dimensional (3-D) S wave velocity structure beneath a rural town in Zhejiang, China. We report the advantage of cross-coherence over cross-correlation in suppressing pseudo-arrivals caused by persistent sources. We use ambient noise interferometry to retrieve high quality Rayleigh waves and Love waves. Body waves are also observed on the R-R component interferograms. We apply phase velocity dispersion measurements on both Rayleigh waves and Love waves and automatically pick more than 23,000 dispersion curves by using a Machine Learning technique. 3-D surface wave tomographic results after depth inversion indicate low-velocity anomalies (between −1% and −4%) from the surface to 2 km depth in the central area. Combined with the conductive characteristics observed on resistivity profile, the low-velocity anomalies are inferred to be a fluid saturated zone of highly fractured rock. Joint interpretation based on horizontal-to-vertical spectral ratio (HVSR) measurements, and existing temperature and fluid resistivity records observed in a nearby well, suggests the existence of the high-temperature geothermal field through the fracture channel. Strong correlation between HVSR measurements and the S wave velocity model highlights the potential of extraction of both amplitude and phase information from ambient noise.

The Imperial Valley, CA, is a tectonically active transtensional basin located south of the Salton Sea; the area hosts numerous geothermal fields, including significant hidden hydrothermal resources without surface manifestations. Development of inexpensive, rugged, and highly sensitive exploration techniques for undiscovered geothermal systems is critical for accelerating geothermal power deployment as well as unlocking a low‐carbon energy future. We present a case study utilizing distributed acoustic sensing (DAS) and ambient noise interferometry for geothermal reservoir imaging, utilizing unlit fiber‐optic telecommunication infrastructure (dark fiber). The study exploits two days of passive DAS data acquired in early November 2020 over a ∼28‐km section of fiber from Calipatria, CA to Imperial, CA. We apply ambient noise interferometry to retrieve coherent signals from DAS records and develop a bin stacking technique to attenuate the effects from persistent localized noise sources and to enhance retrieval of coherent surface waves. As a result, we are able to obtain high‐resolution two‐dimensional (2D) S wave velocity (V_s) structure to 3 km depth, based on joint inversion of both the fundamental and higher overtones. We observe a previously unmapped high V_s and low V_p/V_s ratio feature beneath the Brawley geothermal system, which we interpret to be a zone of hydrothermal mineralization and lower porosity. This interpretation is consistent with a host of other measurements including surface heat flow, gravity anomalies, and available borehole wireline data. These results demonstrate the potential utility of DAS deployed on dark fiber for geothermal system exploration and characterization in the appropriate geological settings.