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.
- Award ID(s):
- 1953710
- NSF-PAR ID:
- 10334031
- 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
-
SUMMARY -
We have conducted three-dimensional (3D) 0–7.5 Hz physics-based wave propagation simulations to model the seismic response of the Long Valley Dam (LVD), which has formed Lake Crowley in Central California, to estimate peak ground motions and settlement of the dam expected during maximum credible earthquake (MCE) scenarios on the nearby Hilton Creek Fault (HCF). We calibrated the velocity structure, anelastic attenuation model, and the overall elastic properties of the dam via linear simulations of a Mw3.7 event as well as the Mw6.2 Chalfant Valley earthquake of 1986, constrained by observed ground motions on and nearby the LVD. The Statewide California Earthquake Center (SCEC) Community Velocity Model CVM-S4.26.M01 superimposed with a geotechnical layer using [Formula: see text] information tapered from the surface to a 700-m depth was used in the simulations. We found optimal fit of simulated and observed ground motions at the LVD using frequency-independent attenuation of [Formula: see text] ([Formula: see text] in m/s). Using the calibrated model, we simulated 3D nonlinear ground motions at the LVD for Mw6.6 rupture scenarios on the HCF using an Iwan-type, multi-yield-surface technique. We use a two-step method where the computationally expensive nonlinear calculations were carried out in a small domain with the plane wave excitation along the bottom boundary obtained from a full-domain 3D linear finite-fault simulation. Our nonlinear MCE simulation results show that peak ground velocities (PGVs) and peak ground accelerations (PGAs) as high as 72 cm/s and 0.55 g, respectively, can be expected at the crest of the LVD. Compared with linear ground motion simulation results, our results show that Iwan nonlinear damping reduces PGAs on the dam crest by up to a factor of 8 and increasingly depletes the high-frequency content of the waves toward the dam crest. We find horizontal relative displacements of the material inside the dam of up to [Formula: see text] and up to [Formula: see text] of vertical subsidence, equivalent to 1% of the dam height.
-
Abstract The Húsavík‐Flatey Fault Zone (HFFZ) is the largest strike‐slip fault in Iceland and poses a high seismic risk to coastal communities. To investigate physics‐based constraints on earthquake hazards, we construct three fault system models of varying geometric complexity and model 79 3‐D multi‐fault dynamic rupture scenarios in the HFFZ. By assuming a simple regional prestress and varying hypocenter locations, we analyze the rupture dynamics, fault interactions, and the associated ground motions up to 2.5 Hz. All models account for regional seismotectonics, topo‐bathymetry, 3‐D subsurface velocity, viscoelastic attenuation, and off‐fault plasticity, and we explore the effect of fault roughness. The rupture scenarios obey earthquake scaling relations and predict magnitudes comparable to those of historical events. We show how fault system geometry and segmentation, hypocenter location, and prestress can affect the potential for rupture cascading, leading to varying slip distributions across different portions of the fault system. Our earthquake scenarios yield spatially heterogeneous near‐field ground motions modulated by geometric complexities, topography, and rupture directivity, particularly in the near‐field. The average ground motion attenuation characteristics of dynamic rupture scenarios of comparable magnitudes and mean stress drop are independent of variations in source complexity, magnitude‐consistent and in good agreement with the latest regional empirical ground motion models. However, physics‐based ground motion variability changes considerably with fault‐distance and increases for unilateral compared to bilateral ruptures. Systematic variations in physics‐based near‐fault ground motions provide important insights into the mechanics and potential earthquake hazard of large strike‐slip fault systems, such as the HFFZ.
-
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
-
SUMMARY The near-surface seismic structure (to a depth of about 1000 m), particularly the shear wave velocity (VS), can strongly affect the propagation of seismic waves and, therefore, must be accurately calibrated for ground motion simulations and seismic hazard assessment. The VS of the top (<300 m) crust is often well characterized from borehole studies, geotechnical measurements, and water and oil wells, while the velocities of the material deeper than about 1000 m are typically determined by tomography studies. However, in depth ranges lacking information on shallow lithological stratification, typically rock sites outside the sedimentary basins, the material parameters between these two regions are typically poorly characterized due to resolution limits of seismic tomography. When the alluded geological constraints are not available, models, such as the Southern California Earthquake Center (SCEC) Community Velocity Models (CVMs), default to regional tomographic estimates that do not resolve the uppermost VS values, and therefore deliver unrealistically high shallow VS estimates. The SCEC Unified Community Velocity Model (UCVM) software includes a method to incorporate the near-surface earth structure by applying a generic overlay based on measurements of time-averaged VS in top 30 m (VS30) to taper the upper part of the model to merge with tomography at a depth of 350 m, which can be applied to any of the velocity models accessible through UCVM. However, our 3-D simulations of the 2014 Mw 5.1 La Habra earthquake in the Los Angeles area using the CVM-S4.26.M01 model significantly underpredict low-frequency (<1 Hz) ground motions at sites where the material properties in the top 350 m are significantly modified by the generic overlay (‘taper’). On the other hand, extending the VS30-based taper of the shallow velocities down to a depth of about 1000 m improves the fit between our synthetics and seismic data at those sites, without compromising the fit at well-constrained sites. We explore various tapering depths, demonstrating increasing amplification as the tapering depth increases, and the model with 1000 m tapering depth yields overall favourable results. Effects of varying anelastic attenuation are small compared to effects of velocity tapering and do not significantly bias the estimated tapering depth. Although a uniform tapering depth is adopted in the models, we observe some spatial variabilities that may further improve our method.