Title: Strong correlation between stress drop and peak ground acceleration for recent M1–4 earthquakes in the San Francisco Bay Area

Abstract

Theoretical and observational studies suggest that between-event variability in the median ground motions of larger ( M≥5 ) earthquakes is controlled primarily by the dynamic properties of the earthquake source, such as Brune-type stress drop. Analogous results remain equivocal for smaller events due to the lack of comprehensive and overlapping ground-motion and source-parameter datasets in this regime. Here in this paper, we investigate the relationship between peak ground acceleration (PGA) and dynamic stress drop for a new dataset of 5297 earthquakes that occurred in the San Francisco Bay area from 2002 through 2016. For each event, we measure PGA on horizontal-component channels of stations within 100 km and estimate stress drop from P-wave spectra recorded on vertical-component channels of the same stations. We then develop a nonparametric ground-motion prediction equation (GMPE) applicable for the moderate (M 1–4) earthquakes in our study region, using a mixed-effects generalization of the Random Forest algorithm. We use the Random Forest GMPE to model the joint influence of magnitude, distance, and near-site effects on observed PGA. We observe a strong correlation between dynamic stress drop and the residual PGA of each event, with the events with higher-than-expected PGA associated with higher values of stress drop.more » The strength of this correlation increases as a function of magnitude but remains significant even for smaller magnitude events with corner frequencies that approach the observable bandwidth of the acceleration records. Mainshock events are characterized by systematically higher stress drop and PGA than aftershocks of equivalent magnitude. Coherent local variations in the distribution of dynamic stress drop provide observational constraints to support the future development of nonergodic GMPEs that account for variations in median stress drop at different source locations.« less

@article{osti_1430024,
title = {Strong correlation between stress drop and peak ground acceleration for recent M1–4 earthquakes in the San Francisco Bay Area},
author = {Trugman, Daniel Taylor and Shearer, Peter M.},
abstractNote = {Theoretical and observational studies suggest that between-event variability in the median ground motions of larger ( M≥5 ) earthquakes is controlled primarily by the dynamic properties of the earthquake source, such as Brune-type stress drop. Analogous results remain equivocal for smaller events due to the lack of comprehensive and overlapping ground-motion and source-parameter datasets in this regime. Here in this paper, we investigate the relationship between peak ground acceleration (PGA) and dynamic stress drop for a new dataset of 5297 earthquakes that occurred in the San Francisco Bay area from 2002 through 2016. For each event, we measure PGA on horizontal-component channels of stations within 100 km and estimate stress drop from P-wave spectra recorded on vertical-component channels of the same stations. We then develop a nonparametric ground-motion prediction equation (GMPE) applicable for the moderate (M 1–4) earthquakes in our study region, using a mixed-effects generalization of the Random Forest algorithm. We use the Random Forest GMPE to model the joint influence of magnitude, distance, and near-site effects on observed PGA. We observe a strong correlation between dynamic stress drop and the residual PGA of each event, with the events with higher-than-expected PGA associated with higher values of stress drop. The strength of this correlation increases as a function of magnitude but remains significant even for smaller magnitude events with corner frequencies that approach the observable bandwidth of the acceleration records. Mainshock events are characterized by systematically higher stress drop and PGA than aftershocks of equivalent magnitude. Coherent local variations in the distribution of dynamic stress drop provide observational constraints to support the future development of nonergodic GMPEs that account for variations in median stress drop at different source locations.},
doi = {10.1785/0120170245},
journal = {Bulletin of the Seismological Society of America},
number = 2,
volume = 108,
place = {United States},
year = {2018},
month = {3}
}

Figures / Tables:

Figure 1: Map view of the San Francisco Bay area study region: M ≥ 1:0 seismicity from the relocated catalog of Waldhauser and Schaff (2008); focal mechanisms for M ≥ 3:5 events, Northern California Seismic Network station coverage (triangular symbols), and mapped fault structures (see Data and Resources) are shownmore » for reference. The color version of this figure is available only in the electronic edition.« less

The objectives of this study are to improve low-magnitude (concentrating on M2.5-5) regional seismic discrimination by performing a thorough investigation of earthquake source scaling using diverse, high-quality datasets from varied tectonic regions. Local-to-regional high-frequency discrimination requires an estimate of how earthquakes scale with size. Walter and Taylor (2002) developed the MDAC (Magnitude and Distance Amplitude Corrections) method to empirically account for these effects through regional calibration. The accuracy of these corrections has a direct impact on our ability to identify clandestine explosions in the broad regional areas characterized by low seismicity. Unfortunately our knowledge at small magnitudes (i.e., m{sub b}more » < {approx} 4.0) is poorly resolved, and source scaling remains a subject of on-going debate in the earthquake seismology community. Recently there have been a number of empirical studies suggesting scaling of micro-earthquakes is non-self-similar, yet there are an equal number of compelling studies that would suggest otherwise. It is not clear whether different studies obtain different results because they analyze different earthquakes, or because they use different methods. Even in regions that are well studied, such as test sites or areas of high seismicity, we still rely on empirical scaling relations derived from studies taken from half-way around the world at inter-plate regions. We investigate earthquake sources and scaling from different tectonic settings, comparing direct and coda wave analysis methods that both make use of empirical Green's function (EGF) earthquakes to remove path effects. Analysis of locally recorded, direct waves from events is intuitively the simplest way of obtaining accurate source parameters, as these waves have been least affected by travel through the earth. But finding well recorded earthquakes with 'perfect' EGF events for direct wave analysis is difficult, limits the number of earthquakes that can be studied. We begin with closely-located, well-correlated earthquakes. We use a multi-taper method to obtain time-domain source-time-functions by frequency division. We only accept an earthquake and EGF pair if they are able to produce a clear, time-domain source pulse. We fit the spectral ratios and perform a grid-search about the preferred parameters to ensure the fits are well constrained. We then model the spectral (amplitude) ratio to determine source parameters from both direct P and S waves. We analyze three clusters of aftershocks from the well-recorded sequence following the M5 Au Sable Forks, NY, earthquake to obtain some of the first accurate source parameters for small earthquakes in eastern North America. Each cluster contains a M{approx}2, and two contain M{approx}3, as well as smaller aftershocks. We find that the corner frequencies and stress drops are high (averaging 100 MPa) confirming previous work suggesting that intraplate continental earthquakes have higher stress drops than events at plate boundaries. We also demonstrate that a scaling breakdown suggested by earlier work is simply an artifact of their more band-limited data. We calculate radiated energy, and find that the ratio of Energy to seismic Moment is also high, around 10{sup -4}. We estimate source parameters for the M5 mainshock using similar methods, but our results are more doubtful because we do not have a EGF event that meets our preferred criteria. The stress drop and energy/moment ratio for the mainshock are slightly higher than for the aftershocks. Our improved, and simplified coda wave analysis method uses spectral ratios (as for the direct waves) but relies on the averaging nature of the coda waves to use EGF events that do not meet the strict criteria of similarity required for the direct wave analysis. We have applied the coda wave spectral ratio method to the 1999 Hector Mine mainshock (M{sub w} 7.0, Mojave Desert) and its larger aftershocks, and also to several sequences in Italy with M{approx}6 mainshocks. The Italian earthquakes have higher stress drops than the Hector Mine sequence, but lower than Au Sable Forks. These results show a departure from self-similarity, consistent with previous studies using similar regional datasets. The larger earthquakes have higher stress drops and energy/moment ratios. We perform a preliminary comparison of the two methods using the M5 Au Sable Forks earthquake. Both methods give very consistent results, and we are applying the comparison to further events.« less

The objectives of this study are to improve low-magnitude regional seismic discrimination by performing a thorough investigation of earthquake source scaling using diverse, high-quality datasets from varied tectonic regions. Local-to-regional high-frequency discrimination requires an estimate of how earthquakes scale with size. Walter and Taylor (2002) developed the MDAC (Magnitude and Distance Amplitude Corrections) method to empirically account for these effects through regional calibration. The accuracy of these corrections has a direct impact on our ability to identify clandestine explosions in the broad regional areas characterized by low seismicity. Unfortunately our knowledge of source scaling at small magnitudes (i.e., m{sub b}more » < {approx}4.0) is poorly resolved. It is not clear whether different studies obtain contradictory results because they analyze different earthquakes, or because they use different methods. Even in regions that are well studied, such as test sites or areas of high seismicity, we still rely on empirical scaling relations derived from studies taken from half-way around the world at inter-plate regions. We investigate earthquake sources and scaling from different tectonic settings, comparing direct and coda wave analysis methods. We begin by developing and improving the two different methods, and then in future years we will apply them both to each set of earthquakes. Analysis of locally recorded, direct waves from events is intuitively the simplest way of obtaining accurate source parameters, as these waves have been least affected by travel through the earth. But there are only a limited number of earthquakes that are recorded locally, by sufficient stations to give good azimuthal coverage, and have very closely located smaller earthquakes that can be used as an empirical Green's function (EGF) to remove path effects. In contrast, coda waves average radiation from all directions so single-station records should be adequate, and previous work suggests that the requirements for the EGF event are much less stringent. We can study more earthquakes using the coda-wave methods, while using direct wave methods for the best recorded subset of events so as to investigate any differences between the results of the two approaches. Finding 'perfect' EGF events for direct wave analysis is difficult, as is ascertaining the quality of a particular EGF event. We develop a multi-taper method to obtain time-domain source-time-functions by frequency division. If an earthquake and EGF event pair are able to produce a clear, time-domain source pulse then we accept the EGF event. We then model the spectral (amplitude) ratio to determine source parameters from both direct P and S waves. We use the well-recorded sequence of aftershocks of the M5 Au Sable Forks, NY, earthquake to test the method and also to obtain some of the first accurate source parameters for small earthquakes in eastern North America. We find that the stress drops are high, confirming previous work suggesting that intraplate continental earthquakes have higher stress drops than events at plate boundaries. We simplify and improve the coda wave analysis method by calculating spectral ratios between different sized earthquakes. We first compare spectral ratio performance between local and near-regional S and coda waves in the San Francisco Bay region for moderate-sized events. The average spectral ratio standard deviations using coda are {approx}0.05 to 0.12, roughly a factor of 3 smaller than direct S-waves for 0.2 < f < 15.0 Hz. Also, direct wave analysis requires collocated pairs of earthquakes whereas the event-pairs (Green's function and target events) can be separated by {approx}25 km for coda amplitudes without any appreciable degradation. We then apply coda spectral ratio method to the 1999 Hector Mine mainshock (M{sub w} 7.0, Mojave Desert) and its larger aftershocks. We observe a clear departure from self-similarity, consistent with previous studies using similar regional datasets.« less

According to the article, we numerically model broadband ground motion (up to 5–7.5 Hz⁠) from blind‐thrust scenario earthquakes matching the fault geometry of the 1994 Mw 6.7 Northridge earthquake. Several realizations are modeled (by varying the hypocenter location in the dynamic rupture simulation) in a 1D‐layered velocity profile. In addition, we include Q(f )⁠, nonlinear effects from Drucker–Prager plasticity, and superimpose small‐scale medium complexity in both a 1D‐layered and 3D velocity model within the subsequent wave propagation. We investigate characteristics of the ground motion and its variability up to 50 km from the fault by comparing them with ground‐motion prediction equations (GMPEs),more » simple proxy metrics, as well as strong ground motion records from the Northridge event. We find that median ground motion closely follows the trend predicted by GMPEs and that the intraevent standard deviation, although varying with hypocenter location, lies near that of GMPE models. Plasticity affects ground‐motion amplitudes in regions near the source, reducing intraevent variability above ∼0.5 Hz⁠. Heterogeneity in the velocity structure on both the regional and small scales is needed for the simulated data to match two proxy metrics: the period‐to‐period correlation of spectral acceleration (SA) and the ratio of maximum‐to‐median SA. Although small‐scale heterogeneity has a negligible effect on median SA for this style of rupture, it serves to significantly increase the cumulative absolute velocity, better agreeing with observations. When compared with strong‐motion data, we find that long‐wavelength velocity structure within our deterministic simulations reduces bias at both short and long periods. Finally, synthetic ground motion at both footwall and hanging‐wall sites has no clear dependence on the distance to rupture (at both short and long periods); directivity is likely overpowering any hanging‐wall effect.« less

In this paper, we present model simulations of ground motions caused by CO 2 -injection-induced fault reactivation and analyze the results in terms of the potential for damage to ground surface structures and nuisance to the local human population. It is an integrated analysis from cause to consequence, including the whole chain of processes starting from earthquake inception in the subsurface, wave propagation toward the ground surface, and assessment of the consequences of ground vibration. For a small magnitude (M w =3) event at a hypocenter depth of about 1000m, we first used the simulated ground-motion wave train in anmore » inverse analysis to estimate source parameters (moment magnitude, rupture dimensions and stress drop), achieving good agreement and thereby verifying the modeling of the chain of processes from earthquake inception to ground vibration. We then analyzed the ground vibration results in terms of peak ground acceleration (PGA), peak ground velocity (PGV) and frequency content, with comparison to U.S. Geological Survey's instrumental intensity scales for earthquakes and the U.S. Bureau of Mines' vibration criteria for cosmetic damage to buildings, as well as human-perception vibration limits. Our results confirm the appropriateness of using PGV (rather than PGA) and frequency for the evaluation of potential ground-vibration effects on structures and humans from shallow injection-induced seismic events. For the considered synthetic M w =3 event, our analysis showed that the short duration, high frequency ground motion may not cause any significant damage to surface structures, but would certainly be felt by the local population.« less

In the article, we model deterministic broadband (-7.5 Hz) ground motion from an M w 7.1 bilateral strike-slip earthquake scenario with dynamic rupture propagation along a rough-fault topography embedded in a medium including small-scale velocity and density perturbations. Spectral accelerations (SAs) at periods 0.2-3 s and Arias intensity durations show a similar distance decay (at the level of 1-2 interevent standard deviations above the median) when compared to Next Generation Attenuation-West2 (NGA)-West2 ground-motion prediction equations (GMPEs) using a Q(f) power-law exponent of 0.6-0.8 above 1 Hz in models with a minimum V S of 750 m/s. With a trade-off frommore » Q(f), the median ground motion is slightly increased by scattering from statistical models of small-scale heterogeneity with standard deviation (σ) of the perturbations at the lower end of the observed range (5%) but reduced by scattering attenuation at the upper end (10%) when using a realistic 3D background velocity model. The ground-motion variability is strongly affected by the addition of small-scale media heterogeneity, reducing otherwise large values of intraevent standard deviation closer to those of empirical observations. These simulations generally have intraevent standard deviations for SAs lower than the GMPEs for the modeled bandwidth, with an increasing trend with distance (most pronounced in low-to-moderate scattering media) near the level of observations at distances greater than 35 km from the fault. Durations for the models follow the same increasing trend with distance, in which σ~5% produces the best match to GMPE values. We find that a 3D background-velocity model reduces the pulse period into the expected range by breaking up coherent waves from directivity, generating a lognormal distribution of ground-motion residuals. These results indicate that a strongly heterogeneous medium is needed to produce realistic deterministic broadband ground motions. Finally, the addition of a thin surficial layer with low, frequency-independent Q in the model (with a minimum V S of 750 m/s) controls the high-frequency decay in energy, as measured by the parameter κ, that may be necessary to include as simulations continue to extend to higher frequencies.« less