skip to main content

Title: Denoising Surface Waves Extracted From Ambient Noise Recorded by 1‐D Linear Array Using Three‐Station Interferometry of Direct Waves

We develop an automatic workflow for enhancing surface wave signals in ambient noise cross correlations (ANCs) calculated for a one‐dimensional (1‐D) linear array. The proposed array‐based method is applied to a 1.6 km‐long dense linear nodal array crossing surface traces of the San Jacinto fault near Anza, California. Fundamental and higher modes of surface waves are observed in ANCs of the nodal array. After attenuating the surface wave overtones by applying a frequency‐dependent tapering window to the ANCs, signals dominated by the fundamental mode surface wave are then enhanced through a denoising process based on three‐station interferometry of direct waves. The signal‐to‐noise ratio is significantly increased at high frequencies (>2 Hz) after denoising. Phase travel times are extracted reliably in the frequency domain for the period ranges of 0.3–1.2 s and 0.3–1.6 s for Rayleigh and Love waves, respectively. The corresponding period‐dependent phase velocity profiles derived from the eikonal equation reveal high‐resolution details of fault zone internal structures beneath the array. A broad (500–1,000 m) low‐velocity zone that narrows with increasing period is observed, illuminating a flower‐shaped structure of the San Jacinto fault damage zone.

more » « less
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present observations and modeling of spatial eigen‐functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5‐Hz geophones that cross the SJFZ with ~10–30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after theSwave arrival are consistent for more than 50 near‐fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio show three main frequency peaks at ~1.3, ~2.0, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen‐functions of resonance waves in a low‐velocity layer (fault zone) between two quarter‐spaces. Using a grid search‐based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65%Swave velocity reduction compared to the host rock. The SW edge of the low‐velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together.

    more » « less
  2. Abstract

    We use Eikonal tomography to derive phase and group velocities of surface waves for the plate boundary region in Southern California. Seismic noise data in the period range 2 and 20 s recorded in year 2014 by 346 stations with ~1‐ to 30‐km station spacing are analyzed. Rayleigh and Love wave phase travel times are measured using vertical‐vertical and transverse‐transverse noise cross correlations, and group travel times are derived from the phase measurements. Using the Eikonal equation for each location and period, isotropic phase and group velocities and 2‐psi azimuthal anisotropy are determined statistically with measurements from different virtual sources. Starting with the SCEC Community Velocity Model, the observed 2.5‐ to 16‐s isotropic phase and group dispersion curves are jointly inverted on a 0.05° × 0.05° grid to obtain local 1‐D piecewise shear wave velocity (Vs) models. Compared to the starting model, the final results have generally lowerVsin the shallow crust (top 3–10 km), particularly in areas such as basins and fault zones. The results also show clear velocity contrasts across the San Andreas, San Jacinto, Elsinore, and Garlock Faults and suggest that the San Andreas Fault southeast of San Gorgonio Pass is dipping to the northeast. Investigation of the nonuniqueness of the 1‐DVsinversion suggests that imaging the top 3‐kmVsstructure requires either shorter period (≤2 s) surface wave dispersion measurements or other types of data set such as Rayleigh wave ellipticity.

    more » « less

    We study temporal changes of seismic velocities associated with the 10 June 2016 Mw 5.2 Borrego Springs earthquake in the San Jacinto fault zone, using nine component Green's function estimates reconstructed from daily cross correlations of ambient noise. The analysed data are recorded by stations in two dense linear arrays, at Dry Wash (DW) and Jackass Flat (JF), crossing the fault surface trace ∼3 km northwest and southeast of the event epicentre. The two arrays have 9 and 12 stations each with instrument spacing of 25–100 m. Relative velocity changes (δv/v) are estimated from arrival time changes in the daily correlation coda waveforms compared to a reference stack. The obtained array-average δv/v time-series exhibit changes associated with the Borrego Springs event, superposed with seasonal variations. The earthquake-related changes are characterized by a rapid coseismic velocity drop followed by a gradual recovery. This is consistently observed at both arrays using time- and frequency-domain δv/v analyses with data from different components in various frequency bands. Larger changes at lower frequencies imply that the variations are not limited to the near surface material. A decreasing coseismic velocity reduction with coda wave lapse time indicates larger coseismic structural perturbations in the fault zone and near-fault environment compared to the surrounding rock. Observed larger changes at the DW array compared to the JF array possibly reflect the northwestward rupture directivity of the Borrego Springs earthquake.

    more » « less
  4. null (Ed.)
    Abstract We image the shallow structure across the East Bench segment of the Wasatch fault system in Salt Lake City using ambient noise recorded by a month-long temporary linear seismic array of 32 stations. We first extract Rayleigh-wave signals between 0.4 and 1.1 s period using noise cross correlation. We then apply double beamforming to enhance coherent cross-correlation signals and at the same time measure frequency-dependent phase velocities across the array. For each location, based on available dispersion measurements, we perform an uncertainty-weighted least-squares inversion to obtain a 1D VS model from the surface to 400 m depth. We put all piece-wise continuous 1D models together to construct the final 2D VS model. The model reveals high velocities to the east of the Pleistocene Lake Bonneville shoreline reflecting thinner sediments and low velocities particularly in the top 200 m to the west corresponding to the Salt Lake basin sediments. In addition, there is an ∼400-m-wide low-velocity zone that narrows with depth adjacent to the surface trace of the East Bench fault, which we interpret as a fault-related damage zone. The damage zone is asymmetric, wider on the hanging wall (western) side and with greater velocity reduction. These results provide important constraints on normal-fault earthquake mechanics, Wasatch fault earthquake behavior, and urban seismic hazard in Salt Lake City. 
    more » « less
  5. Abstract

    We image the shallow seismic structure across the Southern San Andreas Fault (SSAF) using signals from freight trains and trucks recorded by a dense nodal array, with a linear component perpendicular to SSAF and two 2D subarrays centered on the Banning Fault and Mission Creek Fault (MCF). Particle motion analysis in the frequency band 2–5 Hz shows that the examined traffic sources can be approximated as moving single‐ or multi‐point sources that primarily induce Rayleigh waves. Using several techniques, we resolve strong lateral variations of Rayleigh wave velocities andQ‐values across the SSAF, including 35% velocity reduction across MCF toward the northeast and strong attenuation around the two fault strands. We further resolve 10% mass density reduction and 45% shear modulus decrease across the MCF. These findings suggest that the MCF is currently the main strand of the SSAF in the area with important implications for seismic hazard assessments.

    more » « less