Geology topics

R.W. Graves

Publications and source records attributed to R.W. Graves.

8 recordsLinked to original sources

USGS2012

Validation of ground-motion simulations for historical events using SDoF systems

The study presented in this paper is among the first in a series of studies toward the engineering validation of the hybrid broadband ground‐motion simulation methodology by Graves and Pitarka (2010). This paper provides a statistical comparison between seismic demands of single degree of freedom (SDoF) systems subjected to past events using simulations and actual recordings. A number of SDoF systems are selected considering the following: (1) 16 oscillation periods between 0.1 and 6 s; (2) elastic case and four nonlinearity levels, from mildly inelastic to severely inelastic systems; and (3) two hysteretic behaviors, in particular, nondegrading–nonevolutionary and degrading–evolutionary. Demand spectra are derived in terms of peak and cyclic response, as well as their statistics for four historical earthquakes: 1979 M w 6.5 Imperial Valley, 1989 M w 6.8 Loma Prieta, 1992 M w 7.2 Landers, and 1994 M w 6.7 Northridge.

Bulletin of the Seismological Society of America↗

USGS2011

Comparison of ground motions from hybrid simulations to nga prediction equations

We compare simulated motions for a Mw 7.8 rupture scenario on the San Andreas Fault known as the ShakeOut event, two permutations with different hypocenter locations, and a Mw 7.15 Puente Hills blind thrust scenario, to median and dispersion predictions from empirical NGA ground motion prediction equations. We find the simulated motions attenuate faster with distance than is predicted by the NGA models for periods less than about 5.0 s After removing this distance attenuation bias, the average residuals of the simulated events (i.e., event terms) are generally within the scatter of empirical event terms, although the ShakeOut simulation appears to be a high static stress drop event. The intraevent dispersion in the simulations is lower than NGA values at short periods and abruptly increases at 1.0 s due to different simulation procedures at short and long periods. The simulated motions have a depth-dependent basin response similar to the NGA models, and also show complex effects in which stronger basin response occurs when the fault rupture transmits energy into a basin at low angle, which is not predicted by the NGA models. Rupture directivity effects are found to scale with the isochrone parameter ?? 2011, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS2010

Broadband ground-motion simulation using a hybrid approach

This paper describes refinements to the hybrid broadband ground-motion simulation methodology of Graves and Pitarka (2004), which combines a deterministic approach at low frequencies (f< 1 Hz) with a semistochastic approach at high frequencies (f> 1 Hz). In our approach, fault rupture is represented kinematically and incorporates spatial heterogeneity in slip, rupture speed, and rise time. The prescribed slip distribution is constrained to follow an inverse wavenumber-squared fall-off and the average rupture speed is set at 80% of the local shear-wave velocity, which is then adjusted such that the rupture propagates faster in regions of high slip and slower in regions of low slip. We use a Kostrov-like slip-rate function having a rise time proportional to the square root of slip, with the average rise time across the entire fault constrained empirically. Recent observations from large surface rupturing earthquakes indicate a reduction of rupture propagation speed and lengthening of rise time in the near surface, which we model by applying a 70% reduction of the rupture speed and increasing the rise time by a factor of 2 in a zone extending from the surface to a depth of 5 km. We demonstrate the fidelity of the technique by modeling the strong-motion recordings from the Imperial Valley, Loma Prieta, Landers, and Northridge earthquakes.

Bulletin of the Seismological Society of America↗

USGS2008

Broadband simulations for Mw 7.8 southern san andreas earthquakes: Ground motion sensitivity to rupture speed

Using the high-performance computing resources of the Southern California Earthquake Center, we simulate broadband (0-10 Hz) ground motions for three Mw 7.8 rupture scenarios of the southern San Andreas fault. The scenarios incorporate a kinematic rupture description with the average rupture speed along the large slip portions of the fault set at 0.96, 0.89, and 0.84 times the local shear wave velocity. Consistent with previous simulations, a southern hypocenter efficiently channels energy into the Los Angeles region along the string of basins south of the San Gabriel Mountains. However, we find the basin ground motion levels are quite sensitive to the prescribed rupture speed, with peak ground velocities at some sites varying by over a factor of two for variations in average rupture speed of about 15%. These results have important implications for estimating seismic hazards in Southern California and emphasize the need for improved understanding of earthquake rupture processes. Copyright 2008 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2004

Observed and simulated ground motions in the San Bernardino basin region for the Hector Mine, California, earthquake

During the MW 7.1 Hector Mine earthquake, peak ground velocities recorded at sites in the central San Bernardino basin region were up to 2 times larger and had significantly longer durations of strong shaking than sites just outside the basin. To better understand the effects of 3D structure on the long-period ground-motion response in this region, we have performed finite-difference simulations for this earthquake. The simulations are numerically accurate for periods of 2 sec and longer and incorporate the detailed spatial and temporal heterogeneity of source rupture, as well as complex 3D basin structure. Here, we analyze three models of the San Bernardino basin: model A (with structural constraints from gravity and seismic reflection data), model F (water well and seismic refraction data), and the Southern California Earthquake Center version 3 model (hydrologic and seismic refraction data). Models A and F are characterized by a gradual increase in sediment thickness toward the south with an abrupt step-up in the basement surface across the San Jacinto fault. The basin structure in the SCEC version 3 model has a nearly uniform sediment thickness of 1 km with little basement topography along the San Jacinto fault. In models A and F, we impose a layered velocity structure within the sediments based on the seismic refraction data and an assumed depth-dependent Vp/Vs ratio. Sediment velocities within the SCEC version 3 model are given by a smoothly varying rule-based function that is calibrated to the seismic refraction measurements. Due to computational limitations, the minimum shear-wave velocity is fixed at 600 m/sec in all of the models. Ground-motion simulations for both models A and F provide a reasonably good match to the amplitude and waveform characteristics of the recorded motions. In these models, surface waves are generated as energy enters the basin through the gradually sloping northern margin. Due to the basement step along the San Jacinto fault, the surface wave energy is confined to the region north of this structure, consistent with the observations. The SCEC version 3 model, lacking the basin geometry complexity present in the other two models, fails to provide a satisfactory match to the characteristics of the observed motions. Our study demonstrates the importance of using detailed and accurate basin geometry for predicting ground motions and also highlights the utility of integrating geological, geophysical, and seismological observations in the development and validation of 3D velocity models.

California↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

The seismic response of the Los Angeles basin, California

Using strong-motion data recorded in the Los Angeles region from the 1992 ( M w 7.3) Landers earthquake, we have tested the accuracy of existing three-dimensional (3D) velocity models on the simulation of long-period (≧2 sec) ground motions in the Los Angeles basin and surrounding San Fernando and San Gabriel Valleys. First, the overall pattern and degree of long-period excitation of the basins were identified in the observations. Within the Los Angeles basin, the recorded amplitudes are about three to four times larger than at sites outside the basins; amplitudes within the San Fernando and San Gabriel Valleys are nearly a factor of 3 greater than surrounding bedrock sites. Then, using a 3D finite-difference numerical modeling approach, we analyzed how variations in 3D earth structure affect simulated waveforms, amplitudes, and the fit to the observed patterns of amplification. Significant differences exist in the 3D velocity models of southern California that we tested (Magistrale et al. , 1996; Graves, 1996a; Hauksson and Haase, 1997). Major differences in the models include the velocity of the assumed background models; the depth of the Los Angeles basin; and the depth, location, and geometry of smaller basins. The largest disparities in the response of the models are seen for the San Fernando Valley and the deepest portion of the Los Angeles basin. These arise in large part from variations in the structure of the basins, particularly the effective depth extent, which is mainly due to alternative assumptions about the nature of the basin sediment fill. The general ground-motion characteristics are matched by the 3D model simulations, validating the use of 3D modeling with geologically based velocity-structure models. However, significant shortcomings exist in the overall patterns of amplification and the duration of the long-period response. The successes and limitations of the models for reproducing the recorded ground motions as discussed provide the basis and direction for necessary improvements to earth structure models, whether geologically or tomographically derived. The differences in the response of the earth models tested also translate to variable success in the ability to successfully model the data and add uncertainty to estimates of the basin response given input “scenario” earthquake source models.

California↗