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.