ALBERT

Description: 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 (super 2) for the 30 degrees -dipping fault (geometry similar to the Chi-Chi earthquake), whereas for the 60 degrees - and 75 degrees -dipping faults this velocity is exceeded over an area of 2700 km (super 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.

Description: 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 sub-shear 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.

Description: 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.

Description: 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 (sub 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.