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 lower Vs in 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‐D Vs inversion suggests that imaging the top 3‐km Vs structure requires either shorter period (≤2 s) surface wave dispersion measurements or other types of data set such as Rayleigh wave ellipticity.

S U M M A R YThe internal structure of the Clark fault in the trifurcation area of the San Jacinto fault zone is imaged using seismograms recorded by a dense linear array (Jackass Flat, JF) crossing the surface trace of the fault and an adjacent array (TR) to the SW. Delay times between phase arrivals associated with ∼3500 local earthquakes and nine teleseismic events are used to estimate velocity variations within the arrays. The teleseismic P waves travel faster beneath the TR than the JF array, in contrast to larger scale tomographic results. Statistical analysis of local P-wave delay times indicates that the entire JF array, with an aperture of ∼400 m, is inside a low-velocity damage zone. This low-velocity zone is bounded on the NE side by a shallow bimaterial interface generating fault zone head waves, and it contains an inner zone of more intense damage generating fault zone trapped waves. The P-wave velocity contrast across the local bounding bimaterial interface is 10-15 per cent. The trapping structure is associated with a width of ∼200 m, S-wave velocity reduction of ∼35 per cent with respect to the surrounding rock, Q-value of ∼20 and depth of ∼3.5 km. The imaging results suggest that the main seismogenic fault is near the SW end of the JF array, in agreement with a prominent geomorphologic feature. The existence of intense local damage on the crustal block with faster larger scale velocity at depth is consistent with common propagation of earthquake ruptures in the area to the NW.

A self‐consistent regional‐scale seismic velocity model with resolution from seismogenic depth to the surface is crucial for seismic hazard assessment. Though Southern California is the most seismically imaged region in the world, techniques with high near‐surface sensitivity have been applied only in disparate local areas and have not been incorporated into a unified model with deeper resolution. In the present work, we obtain isotropic values for Rayleigh wave phase velocity and ellipticity in Southern California by cross‐correlating daily time series from the year 2015 across 315 regional stations in period ranges 6 to 18 s. Leveraging the complementary sensitivity of the two Rayleigh wave data sets, we combine H/V and phase velocity measurements to determine a new 3‐D shear velocity model in a Bayesian joint inversion framework. The new model has greatly improved shallow resolution compared to the Southern California Earthquake Center CVMS4.26 reference model. Well‐known large‐scale features common to previous studies are resolved, including velocity contrasts across the San Andreas, San Jacinto, Garlock, and Elsinore faults, midcrustal high‐velocity structure beneath the Mojave Desert, and shallow Moho beneath the Salton Trough. Other prominent features that have previously only been imaged in focused local studies include the correct sedimentary thickness of the southern Central Valley, fold structure of the Ventura and Oak Ridge Anticlines, and velocity contrast across the Newport‐Inglewood fault. The new shallow structure will greatly impact simulation‐based studies of seismic hazard, especially in the near‐surface low‐velocity zones beneath densely populated areas like the Los Angeles, San Bernardino, and Ventura Basins.

S U M M A R YLocal and teleseismic earthquake waveforms recorded by a 180-m-long linear array (BB) with seven seismometers crossing the Clark fault of the San Jacinto fault zone northwest of Anza are used to image a deep bimaterial interface and core damage structure of the fault. Delay times of P waves across the array indicate an increase in slowness from the southwest most (BB01) to the northeast most (BB07) station. Automatic algorithms combined with visual inspection and additional analyses are used to identify local events generating fault zone head and trapped waves. The observed fault zone head waves imply that the Clark fault in the area is a sharp bimaterial interface, with lower seismic velocity on the southwest side. The moveout between the head and direct P arrivals for events within ∼40 km epicentral distance indicates an average velocity contrast across the fault over that section and the top 20 km of 3.2 per cent. A constant moveout for events beyond ∼40 km to the southeast is due to off-fault locations of these events or because the imaged deep bimaterial interface is discontinuous or ends at that distance. The lack of head waves from events beyond ∼20 km to the northwest is associated with structural complexity near the Hemet stepover. Events located in a broad region generate fault zone trapped waves at stations BB04-BB07. Waveform inversions indicate that the most likely parameters of the trapping structure are width of ∼200 m, S velocity reduction of 30-40 per cent with respect to the bounding blocks, Q value of 10-20 and depth of ∼3.5 km. The trapping structure and zone with largest slowness are on the northeast side of the fault. The observed sense of velocity contrast and asymmetric damage across the fault suggest preferred rupture direction of earthquakes to the northwest. This inference is consistent with results of other geological and seismological studies.

in the southern part of the Walker Lane shear zone (Figure 1) were felt throughout southern California and produced a vigorous aftershock sequence. These events led to rapid deployments of seismic arrays across and around the Ridgecrest earthquake sequence (Catchings et al., 2020). Kinematic rupture processes of the Mw 6.4 and Mw 7.1 events, surface deformation, and properties of the aftershocks show complex patterns, with strong variations both along strike of the rupture zones and in depth (e.g.,

Summary 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–2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analyzed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km x 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 years. 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 second 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 multi-scale 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.

We image the internal structure of the San Jacinto fault zone (SJFZ) in the trifurcation area southeast of Anza, California, with seismic records from dense linear and rectangular arrays. The examined data include recordings from more than 20 000 local earthquakes and nine teleseismic events. Automatic detection algorithms and visual inspection are used to identify P and S body waves, along with P-and S-types fault zone trapped waves (FZTW). The location at depth of the main branch of the SJFZ, the Clark fault, is identified from systematic waveform changes across lines of sensors within the dense rectangular array. Delay times of P arrivals from teleseismic and local events indicate damage asymmetry across the fault, with higher damage to the NE, producing a local reversal of the velocity contrast in the shallow crust with respect to the large-scale structure. A portion of the damage zone between the main fault and a second mapped surface trace to the NE generates P-and S-types FZTW. Inversions of high-quality S-type FZTW indicate that the most likely parameters of the trapping structure are width of ∼70 m, S-wave velocity reduction of 60 per cent, Q value of 60 and depth of ∼2 km. The local reversal of the shallow velocity contrast across the fault with respect to large-scale structure is consistent with preferred propagation of earthquake ruptures in the area to the NW.

SUMMARY 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.