Brad T. Aagaard

Implications for prediction and hazard assessment from the 2004 Parkfield earthquake

Obtaining high-quality measurements close to a large earthquake is not easy: one has to be in the right place at the right time with the right instruments. Such a convergence happened, for the first time, when the 28 September 2004 Parkfield, California, earthquake occurred on the San Andreas fault in the middle of a dense network of instruments designed to record it. The resulting data reveal aspects of the earthquake process never before seen. Here we show what these data, when combined with data from earlier Parkfield earthquakes, tell us about earthquake physics and earthquake prediction. The 2004 Parkfield earthquake, with its lack of obvious precursors, demonstrates that reliable short-term earthquake prediction still is not achievable. To reduce the societal impact of earthquakes now, we should focus on developing the next generation of models that can provide better predictions of the strength and location of damaging ground shaking.

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 M_w 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.

Dynamic Earthquake Ruptures in the Presence of Lithostatic Normal Stresses: Implications for Friction Models and Heat Production

We simulate dynamic ruptures on a strike-slip fault in homogeneous and layered half-spaces and on a thrust fault in a layered half-space. With traditional friction models, sliding friction exceeds 50% of the fault normal compressive stress, and unless the pore pressures approach the lithostatic stress, the rupture characteristics depend strongly on the depth, and sliding generates large amounts of heat. Under application of reasonable stress distributions with depth, variation of the effective coefficient of friction with the square root of the shear modulus and the inverse of the depth creates distributions of stress drop and fracture energy that produce realistic rupture behavior. This ad hoc friction model results in (1) low-sliding friction at all depths and (2) fracture energy that is relatively independent of depth. Additionally, friction models with rate-weakening behavior (which form pulselike ruptures) appear to generate heterogeneity in the distributions of final slip and shear stress more effectively than those without such behavior (which form cracklike ruptures). For surface rupture on a thrust fault, the simple slip-weakening friction model, which lacks rate-weakening behavior, accentuates the dynamic interactions between the seismic waves and the rupture and leads to excessively large ground motions on the hanging wall. Waveforms below the center of the fault (which are associated with waves radiated to teleseismic distances) indicate that source inversions of thrust events may slightly underestimate the slip at shallow depths.

Near-Source Ground Motions from Simulations of Sustained Intersonic and Supersonic Fault Ruptures

We examine the long-period near-source ground motions from simulations of M 7.4 events on a strike-slip fault using kinematic ruptures with rupture speeds that range from subshear speeds through intersonic speeds to supersonic speeds. The strong along-strike shear-wave directivity present in scenarios with subshear rupture speeds disappears in the scenarios with ruptures propagating faster than the shear-wave speed. Furthermore, the maximum horizontal displacements and velocities rotate from generally fault-perpendicular orientations at subshear rupture speeds to generally fault-parallel orientations at supersonic rupture speeds. For rupture speeds just above the shear-wave speed, the orientations are spatially heterogeneous as a result of the random nature of our assumed slip model. At locations within a few kilometers of the rupture, the time histories of the polarization of the horizontal motion provide a better diagnostic with which to gauge the rupture speed than the orientation of the peak motion. Subshear ruptures are associated with significant fault-perpendicular motion before fault-parallel motion close to the fault; supershear ruptures are associated with fault-perpendicular motion after significant fault-parallel motion. Consistent with previous studies, we do not find evidence for prolonged supershear rupture in the long-period (>2 sec) ground motions from the 1979 Imperial Valley earthquake. However, we are unable to resolve the issue of whether a limited portion of the rupture (approximately 10 km in length) propagated faster than the shear-wave speed. Additionally, a recording from the 2002 Denali fault earthquake does appear to be qualitatively consistent with locally supershear rupture. Stronger evidence for supershear rupture in earthquakes may require very dense station coverage in order to capture these potentially distinguishing traits.

Characterization of near-source ground motions with earthquake simulations

We examine the characteristics of long-period near-source ground motions by conducting a sensitivity study with variations in six earthquake source parameters for both a strike-slip fault (M 7.0-7.1) and a thrust fault (M 6.6-7.0). The directivity of the ruptures creates large displacement and velocity pulses in the forward direction. The dynamic displacements close to the fault are comparable to the average slip. The ground motions exhibit the greatest sensitivity to the fault depth with moderate sensitivity to the rupture speed, peak slip rate, and average slip. For strike-slip faults and thrust faults with surface rupture, the maximum ground displacements and velocities occur in the region where the near-source factor from the 1997 Uniform Building Code is the largest. However, for a buried thrust fault the peak ground motions can occur up-dip from this region.

Ground-Motion Modeling of the 1906 San Francisco Earthquake, Part I: Validation Using the 1989 Loma Prieta Earthquake

We compute ground motions for the Beroza (1991) and Wald et al. (1991) source models of the 1989 magnitude 6.9 Loma Prieta earthquake using four different wave-propagation codes and recently developed 3D geologic and seismic velocity models. In preparation for modeling the 1906 San Francisco earthquake, we use this well-recorded earthquake to characterize how well our ground-motion simulations reproduce the observed shaking intensities and amplitude and durations of recorded motions throughout the San Francisco Bay Area. All of the simulations generate ground motions consistent with the large-scale spatial variations in shaking associated with rupture directivity and the geologic structure. We attribute the small variations among the synthetics to the minimum shear-wave speed permitted in the simulations and how they accommodate topography. Our long-period simulations, on average, under predict shaking intensities by about one-half modified Mercalli inten- sity (MMI) units (25%-35% in peak velocity), while our broadband simulations, on average, under predict the shaking intensities by one-fourth MMI units (16% in peak velocity). Discrepancies with observations arise due to errors in the source models and geologic structure. The consistency in the synthetic waveforms across the wave- propagation codes for a given source model suggests the uncertainty in the source parameters tends to exceed the uncertainty in the seismic velocity structure. In agree- ment with earlier studies, we find that a source model with slip more evenly distributed northwest and southeast of the hypocenter would be preferable to both the Beroza and Wald source models. Although the new 3D seismic velocity model improves upon previous velocity models, we identify two areas needing improvement. Nevertheless, we find that the seismic velocity model and the wave-propagation codes are suitable for modeling the 1906 earthquake and scenario events in the San Francisco Bay Area. Online Material: Modified Mercalli intensities and velocity waveforms, and a movie of simulated wave propagation.

The ShakeOut Earthquake Source and Ground Motion Simulations

The ShakeOut Scenario is premised upon the detailed description of a hypothetical Mw 7.8 earthquake on the southern San Andreas Fault and the associated simulated ground motions. The main features of the scenario, such as its endpoints, magnitude, and gross slip distribution, were defined through expert opinion and incorporated information from many previous studies. Slip at smaller length scales, rupture speed, and rise time were constrained using empirical relationships and experience gained from previous strong-motion modeling. Using this rupture description and a 3-D model of the crust, broadband ground motions were computed over a large region of Southern California. The largest simulated peak ground acceleration (PGA) and peak ground velocity (PGV) generally range from 0.5 to 1.0 g and 100 to 250 cm/s, respectively, with the waveforms exhibiting strong directivity and basin effects. Use of a slip-predictable model results in a high static stress drop event and produces ground motions somewhat higher than median level predictions from NGA ground motion prediction equations (GMPEs).

Effects of fault dip and slip rake angles on near-source ground motions: Why rupture directivity was minimal in the 1999 Chi-Chi, Taiwan, earthquake

We study how the fault dip and slip rake angles affect near-source ground velocities and displacements as faulting transitions from strike-slip motion on a vertical fault to thrust motion on a shallow-dipping fault. Ground motions are computed for five fault geometries with different combinations of fault dip and rake angles and common values for the fault area and the average slip. The nature of the shear-wave directivity is the key factor in determining the size and distribution of the peak velocities and displacements. Strong shear-wave directivity requires that (1) the observer is located in the direction of rupture propagation and (2) the rupture propagates parallel to the direction of the fault slip vector. We show that predominantly along-strike rupture of a thrust fault (geometry similar in the Chi-Chi earthquake) minimizes the area subjected to large-amplitude velocity pulses associated with rupture directivity, because the rupture propagates perpendicular to the slip vector; that is, the rupture propagates in the direction of a node in the shear-wave radiation pattern. In our simulations with a shallow hypocenter, the maximum peak-to-peak horizontal velocities exceed 1.5 m/sec over an area of only 200 km^2 for the 30°-dipping fault (geometry similar to the Chi-Chi earthquake), whereas for the 60°- and 75°-dipping faults this velocity is exceeded over an area of 2700 km^2. These simulations indicate that the area subjected to large-amplitude long-period ground motions would be larger for events of the same size as Chi-Chi that have different styles of faulting or a deeper hypocenter.

Dynamic Rupture Modeling of the Transition from Thrust to Strike-Slip Motion in the 2002 Denali Fault Earthquake, Alaska

We use three-dimensional dynamic (spontaneous) rupture models to investigate the nearly simultaneous ruptures of the Susitna Glacier thrust fault and the Denali strike-slip fault. With the 1957 M_w 8.3 Gobi-Altay, Mongolia, earthquake as the only other well-documented case of significant, nearly simultaneous rupture of both thrust and strike-slip faults, this feature of the 2002 Denali fault earthquake provides a unique opportunity to investigate the mechanisms responsible for development of these large, complex events. We find that the geometry of the faults and the orientation of the regional stress field caused slip on the Susitna Glacier fault to load the Denali fault. Several different stress orientations with oblique right-lateral motion on the Susitna Glacier fault replicate the triggering of rupture on the Denali fault about 10 sec after the rupture nucleates on the Susitna Glacier fault. However, generating slip directions compatible with measured surface offsets and kinematic source inversions requires perturbing the stress orientation from that determined with focal mechanisms of regional events. Adjusting the vertical component of the principal stress tensor for the regional stress field so that it is more consistent with a mixture of strike-slip and reverse faulting significantly improves the fit of the slip-rake angles to the data. Rotating the maximum horizontal compressive stress direction westward appears to improve the fit even further.

Long-Period Building Response to Earthquakes in the San Francisco Bay Area

This article reports a study of modeled, long-period building responses to ground-motion simulations of earthquakes in the San Francisco Bay Area. The earthquakes include the 1989 magnitude 6.9 Loma Prieta earthquake, a magnitude 7.8 simulation of the 1906 San Francisco earthquake, and two hypothetical magnitude 7.8 northern San Andreas fault earthquakes with hypocenters north and south of San Francisco. We use the simulated ground motions to excite nonlinear models of 20-story, steel, welded moment-resisting frame (MRF) buildings. We consider MRF buildings designed with two different strengths and modeled with either ductile or brittle welds. Using peak interstory drift ratio (IDR) as a performance measure, the stiffer, higher strength building models outperform the equivalent more flexible, lower strength designs. The hypothetical magnitude 7.8 earthquake with hypocenter north of San Francisco produces the most severe ground motions. In this simulation, the responses of the more flexible, lower strength building model with brittle welds exceed an IDR of 2.5% (that is, threaten life safety) on 54% of the urban area, compared to 4.6% of the urban area for the stiffer, higher strength building with ductile welds. We also use the simulated ground motions to predict the maximum isolator displacement of base-isolated buildings with linear, single-degree-of-freedom (SDOF) models. For two existing 3-sec isolator systems near San Francisco, the design maximum displacement is 0.5 m, and our simulations predict isolator displacements for this type of system in excess of 0.5 m in many urban areas. This article demonstrates that a large, 1906-like earthquake could cause significant damage to long-period buildings in the San Francisco Bay Area.