skip to main content

Title: Shallow Damage Zone Structure of the Wasatch Fault in Salt Lake City from Ambient-Noise Double Beamforming with a Temporary Linear Array
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
Award ID(s):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Seismological Research Letters
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY Ocean bottom distributed acoustic sensing (OBDAS) is emerging as a new measurement method providing dense, high-fidelity and broad-band seismic observations from fibre-optic cables deployed offshore. In this study, we focus on 35.7 km of a linear telecommunication cable located offshore the Sanriku region, Japan, and apply seismic interferometry to obtain a high-resolution 2-D shear wave velocity (VS) model below the cable. We first show that the processing steps applied to 13 d of continuous data prior to computing cross-correlation functions (CCFs) impact the modal content of surface waves. Continuous data pre-processed with 1-bit normalization allow us to retrieve dispersion images with high Scholte-wave energy between 0.5 and 5 Hz, whereas spatial aliasing dominates dispersion images above 3 Hz for non-1-bit CCFs. Moreover, the number of receiver channels considered to compute dispersion images also greatly affects the resolution of extracted surface-wave modes. To better understand the remarkably rich modal nature of OBDAS data (i.e. up to 30 higher modes in some regions), we simulate Scholte-wave dispersion curves for stepwise linear VS gradient media. For soft marine sediments, simulations confirm that a large number of modes can be generated in gradient media. Based on pre-processing and theoretical considerations, we extract surface wave dispersion curves from 1-bit CCFs spanning over 400 channels (i.e. ∼2 km) along the array and invert them to image the subsurface. The 2-D velocity profile generally exhibits slow shear wave velocities near the ocean floor that gradually increase with depth. Lateral variations are also observed. Flat bathymetry regions, where sediments tend to accumulate, reveal a larger number of Scholte-wave modes and lower shallow velocity layers than regions with steeper bathymetry. We also compare and discuss the velocity model with that from a previous study and finally discuss the combined effect of bathymetry and shallow VS layers on earthquake wavefields. Our results provide new constraints on the shallow submarine structure in the area and further demonstrate the potential of OBDAS for high-resolution offshore geophysical prospecting. 
    more » « less
  2. 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. 
    more » « less
  3. 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

    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
  5. null (Ed.)
    We present a database and analyze ground motions recorded during three events that occurred as part of the July 2019 Ridgecrest earthquake sequence: a moment magnitude (M) 6.5 foreshock on a left‐lateral cross fault in the Salt Wells Valley fault zone, an M 5.5 foreshock in the Paxton Ranch fault zone, and the M 7.1 mainshock, also occurring in the Paxton Ranch fault zone. We collected and uniformly processed 1483 three‐component recordings from an array of 824 sensors spanning 10 seismographic networks. We developed site metadata using available data and multiple models for the time‐averaged shear‐wave velocity in the upper 30 m (⁠VS30⁠) and for basin depth terms. We processed ground motions using Next Generation Attenuation (NGA) procedures and computed intensity measures including spectral acceleration at a number of oscillator periods and inelastic response spectra. We compared elastic and inelastic response spectra to seismic design spectra in building codes to evaluate the damage potential of the ground motions at spatially distributed sites. Residuals of the observed spectral accelerations relative to the NGA‐West2 ground‐motion models (GMMs) show good average agreement between observations and model predictions (event terms between about −0.3 and 0.5 for peak ground acceleration to 5 s). The average attenuation with distance is also well captured by the empirical NGA‐West2 GMMs, although azimuthal variations in attenuation were observed that are not captured by the GMMs. An analysis considering directivity and fault‐slip heterogeneity for the M 7.1 event demonstrates that the dispersion in the near‐source ground‐motion residuals can be reduced. 
    more » « less