Finite-fault Source Models Research Articles

Influence of Super-Shear Earthquake Rupture Models on Simulated Near-Source Ground Motion from the 1999 Izmit, Turkey, Earthquake

We numerically simulate seismic wave propagation from the 1999 Mw 7.4 earthquake in Izmit, Turkey, using a 3D finite difference method based on published finite-source models obtained by waveform inversions. This earthquake has been reported, based on observations at the near-fault station SKR, as an example of supershear rupture propagation toward the east. Although the modeled ground motion does show a characteristic Mach wave from the fault plane, it is difficult to identify any particular effects in terms of peak ground velocity (PGV), an important parameter in earthquake engineering. This is because the fault's spatial heterogeneity is strong enough to mask the properties of supershear rupture, which has been reported through several numerical simulations mostly based on homogeneous fault conditions. This article demonstrates the importance of studying ground motions for known earthquakes through numerical simulations based on finite-fault source models.

Rupture Process of the 1999 Mw 7.1 Duzce Earthquake from Joint Analysis of SPOT, GPS, InSAR, Strong-Motion, and Teleseismic Data: A Supershear Rupture with Variable Rupture Velocity

We analyze the rupture process of the 1999 M_w 7.1 Duzce earthquake using seismological, remote sensing, and geodetic data. Ground deformation measured from the subpixel cross correlation of Satellite Pour l'Observation de la Terre (SPOT) images reveals a 55 km long fault trace and smooth surface-slip distribution peaking at 3.5–4 m. The westernmost segment overlaps for over 10 km with ruptures from the M_w 7.4 Izmit earthquake. The 15 km long easternmost segment, which cuts across mountainous topography, had not been reported previously. We determine a well-constrained source model using a four-segment fault geometry using constraints on surface fault slip and inverting Global Positioning System and Interferometric Synthetic Aperture Radar data along with strong-motion records. Our results show that some variability of the rupture velocity and an eastward supershear velocity are required to fit the strong-motion data. The rise time, up to 6 sec, correlates with cumulative slip, suggesting a sliding velocity of about 1 m/sec. The source model predicts teleseismic waveforms well, although early by 2 sec. This time shift is probably due to the weak beginning of the earthquake that is not observable at teleseismic distances. Strong-motion records are relatively well predicted from a source model derived from the teleseismic data using the fault geometry derived from the satellite images. This study demonstrates the benefit of using accurate fault geometries to determine finite-fault source models.

Spatial correlation of aftershock locations and on‐fault main shock properties

We quantify the correlation between spatial patterns of aftershock hypocenter locations and the distribution of coseismic slip and stress drop on a main shock fault plane using two nonstandard statistical tests. Test T1evaluates if aftershock hypocenters are located in low‐slip regions (hypothesis H1), test T2evaluates if aftershock hypocenters occur in regions of increased shear stress (hypothesis H2). In the tests, we seek to reject the null hypotheses H0: Aftershock hypocenters are not correlated with (1) low‐slip regions or (2) regions of increased shear stress, respectively. We tested the hypotheses on four strike‐slip events for which multiple earthquake catalogs and multiple finite fault source models of varying accuracy exist. Because we want to retain earthquake clustering as the fundamental feature of aftershock seismicity, we generate slip distributions using a random spatial field model and derive the stress drop distributions instead of generating seismicity catalogs. We account for uncertainties in the aftershock locations by simulating them within their location error bounds. Our findings imply that aftershocks are preferentially located in regions of low‐slip (u≤umax) and of increased shear stress (Δσ < 0). In particular, the correlation is more significant for relocated than for general network aftershock catalogs. However, the results show that stress drop patterns provide less information content on aftershock locations. This implies that static shear stress change of the main shock may not be the governing process for aftershock genesis.

Prediction of near-field strong ground motions for scenario earthquakes on active fault

A method to predict near-field strong ground motions for scenario earthquakes on active faults is proposed. First, macro-source parameters characterizing the entire source area, i.e., global source parameters, including fault length, fault width, rupture area, average slip on the fault plane, etc., are estimated by seismogeology survey, seismicity and seismic scaling laws. Second, slip distributions characterizing heterogeneity or roughness on the fault plane, i.e., local source parameters, are reproduced/evaluated by the hybrid slip model. Finally, the finite fault source model, developed from both the global and local source parameters, is combined with the stochastically synthetic technique of ground motion using the dynamic corner frequency based on seismology. The proposed method is applied to simulate the acceleration time histories on three base-rock stations during the 1994 Northridge earthquake. Comparisons between the predicted and recorded acceleration time histories show that the method is feasible and practicable.