skip to main content

Title: Peak Ground Velocity Spatial Variability Revealed by Dense Seismic Array in Southern California

Understanding and modeling variability of ground motion is essential for building accurate and precise ground motion prediction equations, which can net site‐specific characterization and reduced hazard levels. Here, we explore the spatial variability in peak ground velocity (PGV) at Sage Brush Flats along the San Jacinto Fault in Southern California. We use data from a dense array (0.6 × 0.6 km2, 1,108 geophones, station spacings 10–30 m) deployed in 2014 for ~1 month. These data offer an opportunity to study small‐scale variability in this region. We examine 38 earthquakes (2 ≤ ML ≤ 4.2) within 200 km of the array. Fault strands and a small basin impact the ground motions, producing PGV variations up to 22% of the mean and a 40% reduction inPandSwave near‐surface velocities. We find along‐fault rupture directivity, source, and path effects can increase PGVs by 167%. Surface PGV measurements exceed the colocated borehole station (depth at 148 m) PGV by factors of 3–10, confirming the impact on PGV from near‐surface fault structures, basins, topography, and amplifications from soft sediments. Consistently, we find high PGVs within the basin structure. A pair of colocated GaML2.6 events produce repeatable PGV values with similar spatial patterns. The average corner frequencies of these two events are 11–16 Hz, and viable measurements of stress drop can differ by 6.45 MPa. Within this small array, the PGV values are variable implying spatial extrapolation of PGV to regions of known faults and basins, even across a small area, should be done with caution.

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

    The 14 September 2015 Hildale, Utah, storm resulted in 20 flash flood fatalities, making it the most deadly natural disaster in Utah history; it is the quintessential example of the “paroxysmal precipitation of the desert”. The measured peak discharge from Maxwell Canyon at a drainage area of 5.3 km2was 266 m3/s, a value that exceeds envelope curve peaks for Utah. The 14 September 2015 flash flood reflects features common to other major flash flood events in the region, as well as unique features. The flood was produced by a hailstorm that was moving rapidly from southwest to northeast and intensified as it interacted with complex terrain. Polarimetric radar observations show that the storm exhibited striking temporal variability, with the Maxwell Canyon tributary of Short Creek and a small portion of the East Fork Virgin River basin experiencing extreme precipitation. Periods of extreme rainfall rates for the 14 September 2015 storm are characterized byKDPsignatures of extreme rainfall in polarimetric radar measurements. SimilarKDPsignatures characterized multiple storms that have produced record and near‐record flood peaks in Colorado Plateau watersheds. The climatology of monsoon thunderstorms that produce flash floods exhibits striking spatial heterogeneities in storm occurrence and motion. The hydroclimatology of flash flooding in arid/semiarid watersheds of the southwestern United States exhibits relatively weak dependence on drainage basin area. Large flood peaks over a broad range of basin scales can be produced by small thunderstorms like the 14 September 2015 Hildale Storm, which pass close to the outlet.

    more » « less
  2. Abstract

    The Raton Basin has been an area of injection induced seismicity for the past two decades. Previously, the reactivated fault zone structures and spatiotemporal response of seismicity to evolving injection have been poorly constrained due to sparse publicly available seismic monitoring. The application of a machine‐learning phase picker to 4 years of continuous seismic data from a local array enables the detection and location of ∼38,000 earthquakes. The events from 2016 to 2020 are ∼2.5–6 km below sea level and range from ML < −1 to 4.2. Most earthquakes occur within previously identified ∼N‐S zones of seismicity, however our new catalog illuminates that these zones are composed of many short faults with variable orientations. The two most active zones, the Vermejo Park and Tercio zones, are potentially linked by small intermediate faults. In total, we find ∼60 short (<3 km long) basement faults with strikes from WNW to NNE. Faulting mechanisms are predominantly normal but some variability, including reverse dip‐slip and oblique‐slip, is observed. The Trinidad fault zone, which previously hosted a Mw5.3 earthquake in 2011, is quiescent during 2016–2020, likely in response to both slow accumulation of tectonic strain after the 2011 sequence, and the significant decrease (80% reduction) in nearby wastewater injection from 2012 to 2016. Unlike some other regions, where induced seismicity was triggered in response to higher injection rates, the Raton Basin's frequency‐magnitude and spatiotemporal statistics are not distinguishable from tectonic seismicity. The similarity suggests that seismicity in the Raton Basin is predominantly releasing tectonic stress.

    more » « less

    We have simulated 0–5 Hz deterministic wave propagation for a suite of 17 models of the 2014 Mw 5.1 La Habra, CA, earthquake with the Southern California Earthquake Center Community Velocity Model Version S4.26-M01 using a finite-fault source. Strong motion data at 259 sites within a 148 km × 140 km area are used to validate our simulations. Our simulations quantify the effects of statistical distributions of small-scale crustal heterogeneities (SSHs), frequency-dependent attenuation Q(f), surface topography and near-surface low-velocity material (via a 1-D approximation) on the resulting ground motion synthetics. The shear wave quality factor QS(f) is parametrized as QS, 0 and QS, 0fγ for frequencies less than and higher than 1 Hz, respectively. We find the most favourable fit to data for models using ratios of QS, 0 to shear wave velocity VS of 0.075–1.0 and γ values less than 0.6, with the best-fitting amplitude drop-off for the higher frequencies obtained for γ values of 0.2–0.4. Models including topography and a realistic near-surface weathering layer tend to increase peak velocities at mountain peaks and ridges, with a corresponding decrease behind the peaks and ridges in the direction of wave propagation. We find a clear negative correlation between the effects on peak ground velocity amplification and duration lengthening, suggesting that topography redistributes seismic energy from the large-amplitude first arrivals to the adjacent coda waves. A weathering layer with realistic near-surface low velocities is found to enhance the amplification at mountain peaks and ridges, and may partly explain the underprediction of the effects of topography on ground motions found in models. Our models including topography tend to improve the fit to data, as compared to models with a flat free surface, while our distributions of SSHs with constraints from borehole data fail to significantly improve the fit. Accuracy of the velocity model, particularly the near-surface low velocities, as well as the source description, controls the resolution with which the anelastic attenuation can be determined. Our results demonstrate that it is feasible to use fully deterministic physics-based simulations to estimate ground motions for seismic hazard analysis up to 5 Hz. Here, the effects of, and trade-offs with, near-surface low-velocity material, topography, SSHs and Q(f) become increasingly important as frequencies increase towards 5 Hz, and should be included in the calculations. Future improvement in community velocity models, wider access to computational resources, more efficient numerical codes and guidance from this study are bound to further constrain the ground motion models, leading to more accurate seismic hazard analysis.

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

    We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a possible period of ≈1000 days and an amplitude of Δm≈ 0.6 mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IRJandKsbands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-supernova outbursts in this time period. The IR brightness (MKs=10.7mag) and color (JKs= 1.6 mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature (Teff=35001400+800K) and luminosity (logL/L=5.1±0.2). This places the candidate among the most luminous Type II supernova progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass ofMinit= 17 ± 4M. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling atṀ3×105to 3 × 10−4Myr−1for an assumed wind velocity ofvw= 10 km s−1, perhaps pointing to enhanced mass loss in a pulsation-driven wind.

    more » « less