Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Topographic Response to Simulated Mw 6.5–7.0 Earthquakes on the Seattle Fault
ABSTRACT We explore the response of ground motions to topography during large crustal fault earthquakes by simulating several magnitude 6.5–7.0 rupture scenarios on the Seattle fault, Washington State. Kinematic simulations are run using a 3D spectral element code and a detailed seismic velocity model for the Puget Sound region. This model includes realistic surface topography and a near-surface low-velocity layer; a mesh spacing of ∼30 m at the surface allows modeling of ground motions up to 3 Hz. We simulate 20 earthquake scenarios using different slip distributions and hypocenter locations on a planar fault surface. Results indicate that average ground motions in simulations with and without topography are similar. However, shaking amplification is common at topographic highs, and more than a quarter of all sites experience short-period (≤2 s) ground-motion amplification greater than 25%–35%, compared with models without topography. Comparisons of peak ground velocity at the top and bottom of topographic features demonstrate that amplification is sensitive to period, with the greatest amplifications typically manifesting near a topographic feature’s estimated resonance frequency and along azimuths perpendicular to its primary axis of elongation. However, interevent variability in topographic response can be significant, particularly at shorter periods (<1 s). We do not observe a clear relationship between source centroid-to-site azimuths and the strength of topographic amplification. Overall, our results suggest that although topographic resonance does influence the average ground motions, other processes (e.g., localized focusing and scattering) also play a significant role in determining topographic response. However, the amount of consistent, significant amplification due to topography suggests that topographic effects should likely be considered in some capacity during seismic hazard studies.  more » « less
Award ID(s):
1953710
PAR ID:
10334031
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Bulletin of the Seismological Society of America
Volume:
112
Issue:
3
ISSN:
0037-1106
Page Range / eLocation ID:
1436 to 1462
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT We simulate shaking in Tacoma, Washington, and surrounding areas from Mw 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite-fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region-specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground-motion models (GMMs). However, long-period (T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM-derived short-period (T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different Mw 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short-period shaking adjacent to the fault, while reducing far-field shaking. Intermediate-period shaking is amplified within the Tacoma basin, likely due to enhanced surface-wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near-surface velocity model resolution as waveform simulations progress toward higher frequencies. 
    more » « less
  2. ; (, Seismological Research Letters)
    Abstract The ShakeOut scenario of an M 7.8 northwestward rupture on the southern San Andreas fault (SSAF) (Jones et al., 2008) predicted significant long-period ground-motion amplification in the greater Los Angeles, California, area, caused by a waveguide from interconnected sedimentary basins. However, the early ShakeOut ground-motion simulations omitted important model features with immature versions of the velocity structure and fault geometry. Here, we present 0–1 Hz 3D numerical wave propagation simulations for the ShakeOut scenario including surface topography, as well as updated high-resolution velocity structures and SSAF geometry. Spectral accelerations at 3 s are increased by the local high-resolution basin models (25%–45%) but decreased from complexity in velocity and density updates outside the basins (65%–100%) and inclusion of surface topography (∼30%). The updated model reduces the simulated long-period ground motions in the waveguide by 60%–70%, bringing the predictions significantly closer to the values by a leading Next Generation Attenuation-West2 ground-motion model. 
    more » « less
  3. ; ; (, Geophysical Journal International)
    SUMMARY 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. ; ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT Following the 2019 Mw 7.1 Ridgecrest, California, earthquake, 15 dense 1D and 2D arrays (461 sites) were deployed around the main ruptures, including four 1D arrays across the surface ruptures of the mainshock. The dense arrays captured numerous aftershocks and provided an unprecedented dataset for studying near-fault ground motions. We combine the near-fault ground motions measured from the dense arrays with a regional ground-motion dataset and develop a near-fault partially nonergodic ground-motion model (GMM). The partially nonergodic GMM results are used to investigate several aspects of the near-fault ground motions that are related to source and fault zone properties. We find that shallow events tend to produce weaker short-period energy due to smaller stress drop or stronger attenuation but stronger long-period energy from surface-wave excitation. Significant ground-motion variations across the fault zone are observed, which include amplifications within the damage zones due to both reduced seismic velocities and trapped waves at certain periods. We also observe double-couple radiation pattern effects only in long-period ground motions, agreeing with a period-dependent earthquake rupture process. Between 3 and 15 days following the mainshock, the ground motions of aftershocks decrease with time, likely reflecting decreasing stress drops. No temporal evolution of site response and wave propagation effects within the fault zone are observed between 10 and 35 days following the mainshock. Our study demonstrates the capability of dense array data to improve our understanding of earthquake processes, fault zone properties, and near-fault seismic hazard. 
    more » « less
  5. (, Bulletin of the Seismological Society of America)
    ABSTRACT A depletion of high-frequency ground motions on soil sites has been observed in recent large earthquakes and is often attributed to a nonlinear soil response. Here, I show that the reduced amplitudes of high-frequency horizontal-to-vertical spectral ratios (HVSRs) on soil can also be caused by a smooth crustal velocity model with low shear-wave velocities underneath soil sites. I calculate near-fault ground motions using both 2D dynamic rupture simulations and point-source models for both rock and soil sites. The 1D velocity models used in the simulations are derived from empirical relationships between seismic wave velocities and depths in northern California. The simulations for soil sites feature lower shear-wave velocities and thus larger Poisson’s ratios at shallow depths than those for rock sites. The lower shear-wave velocities cause slower shallow rupture and smaller shallow slip, but both soil and rock simulations have similar rupture speeds and slip for the rest of the fault. However, the simulated near-fault ground motions on soil and rock sites have distinct features. Compared to ground motions on rock, horizontal ground acceleration on soil is only amplified at low frequencies, whereas vertical ground acceleration is deamplified for the whole frequency range. Thus, the HVSRs on soil exhibit a depletion of high-frequency energy. The comparison between smooth and layered velocity models demonstrates that the smoothness of the velocity model plays a critical role in the contrasting behaviors of HVSRs on soil and rock for different rupture styles and velocity profiles. The results reveal the significant role of shallow crustal velocity structure in the generation of high-frequency ground motions on soil sites. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email