skip to main content

Title: Internal Structure of the Central Garlock Fault Zone From Ridgecrest Aftershocks Recorded by Dense Linear Seismic Arrays

We provide high‐resolution seismic imaging of the central Garlock fault using data recorded by two dense seismic arrays that cross the Ridgecrest rupture zone (B4) and the Garlock fault (A5). Analyses of fault zone head waves andP‐wave delay times at array A5 show that the Garlock fault is a sharp bimaterial interface withPwaves traveling ∼5% faster in the northern crustal block. The across‐fault velocity contrast agrees with regional tomography models and generates clearP‐wave reflections in waveforms recorded by array B4. Kirchhoff migration of the reflected waves indicates a near‐vertical fault between 2 and 6 km depth. TheP‐wave delay times imply a ∼300‐m‐wide transition zone near the Garlock fault surface trace beneath array A5, offset to the side with faster velocities. The results provide important constraints for derivations of earthquake properties, simulations of ruptures and ground motion, and future imaging studies associated with the Garlock fault.

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

    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 betweenPwave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and localPwaves 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 amplifySwaveforms, 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,Swave velocity reduction of ∼20% with respect to the surrounding rock,Q‐value of ∼30, andSwave 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.

    more » « less
  2. 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
  3. 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
  4. Abstract

    Fault zone complexities contain important information about factors controlling earthquake dynamic rupture. High‐resolution fault zone imaging requires high‐quality data from dense arrays and new seismic imaging techniques that can utilize large portions of recorded waveforms. Recently, the emerging Distributed Acoustic Sensing (DAS) technique has enabled near‐surface imaging by utilizing existing telecommunication infrastructure and anthropogenic noise sources. With dense sensors at several meters' spacing, the unaliased wavefield can provide unprecedented details for fault zones. In this work, we use a DAS array converted from a 10‐km underground fiber‐optic cable across Ridgecrest City, California. We report clear acausal and coda signals in ambient noise cross‐correlations caused by surface‐to‐surface wave scattering. We use these scattering‐related waves to locate and characterize potential faults. The mapped fault locations are generally consistent with those in the United States Geological Survey Quaternary Fault database of the United States but are more accurate than the extrapolated ones. We also use waveform modeling to infer that a 35 m wide, 90 m deep fault with 30% velocity reduction can best fit the observed scattered coda waves for one of the identified fault zones. These findings demonstrate the potential of DAS for passive imaging of fine‐scale faults in an urban environment.

    more » « less
  5. Abstract

    To discern spatial and explore possible existence of temporal variations of upper crustal anisotropy in an ∼15 km section of the San Jacinto Fault Zone (SJFZ) that is composed of the Buck Ridge and Clark faults in southern California, we conduct a systematic shear wave splitting investigation using local S‐wave data recorded by three broadband seismic stations located near the surface expression of the SJFZ. An automatic data selection and splitting measurement procedure is first applied, and the resulting splitting measurements are then manually screened to ensure reliability of the results. Strong spatial variations in crustal anisotropy are revealed by 1,694 pairs of splitting parameters (fast polarization orientation and splitting delay time), as reflected by the dependence of the resulting splitting parameters on the location and geometry of the raypaths. For raypaths traveling through the fault zones, the fast orientations are dominantly WNW‐ESE which is parallel to the faults and may be attributed to fluid‐filled fractures in the fault zones. For non‐fault‐zone crossing raypaths, the fast orientations are dominantly N–S which are consistent with the orientation of the regional maximum compressive stress. A three‐dimensional model of upper crustal anisotropy is constructed based on the observations. An increase in the raypath length normalized splitting times is observed after the 03/11/2013 M4.7 earthquake, which is probably attributable to changes in the spatial distribution of earthquakes before and after the M4.7 earthquake rather than reflecting temporal changes of upper crustal anisotropy.

    more » « less