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 the S wave 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% S wave 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.

We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times between P wave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and local P waves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low-velocity zones with much slower inner cores that amplify S waveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault-zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high-amplitude waves, identifies candidate fault-zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high-quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km, S wave velocity reduction of ∼20% with respect to the surrounding rock, Q-value of ∼30, and S wave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault-zone internal structures (velocity contrasts and low-velocity zones) that vary along fault strike.

Near-surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near-surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short-period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near-surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation.

Abstract: We discuss general structural features of the Banning and Mission Creek strands (BF and MCF) of the southern San Andreas fault (SSAF) in the Coachella Valley, based on ambient noise and earthquake wavefields recorded by a seismic array with >300 nodes. Earthquake P arrivals show rapid changes in waveform characteristics over 20–40 m zones that coincide with the surface BF and MCF. These variations indicate that the BF and MCF are high-impedance contrast interfaces—an observation supported by the presence of seismic reflections. Another prominent but more diffuse change in SSAF structure is found ∼1 km northeast of the BF. This feature has average-to-low arrival times (P and S) and ambient noise levels (at <30 Hz), and likely represents a relatively fast velocity block sandwiched between broader MCF and BF zones. The maximal arrival delays (P ∼0.1 s and S ∼0.25 s) and the highest ambient noise levels (>2 times median) are consistently observed southwest of the BF—a combined effect of Coachella Valley sediments and rock damage on that side. Immediately northeast of the MCF, large S minus P delays suggest a broad high VP/VS zone associated with asymmetric rock damage across the SSAF. This general overview shows the BF and MCF as mature but distinctly different fault zones. Future analyses will further clarify these and other SSAF features in greater detail.

We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September-October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26-31 in the teleseismic P arrivals compared with the SW-most station, and generates both P-and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.