Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution

Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D _c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D _c / D _tot frequency–magnitude distribution was found. D _c approaches the final slip D _tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.     

New approach in the kinematic k−2 source model for generating physical slip velocity functions

In an attempt to improve the ground motion modelling, the characteristics of the slip velocity functions (SVF) generated using the kinematic k ⁻² source are investigated and compared to the dynamic solutions proposed in the literature. Several numerical simulations were performed to test the influence of the model parameters on the SVF modelling. Overall, the shapes of SVF are very complex and exhibit a large variability in time and space. However, we found out that the mean SVF is a simple boxcar with duration equal to the largest rise time value. In the areas of weak slip, the SVFs are characterized by the existence of negative values, whereas in large slip areas, the SVF is more impulsive. Overall, on the examples investigated, the SVFs modelled with this k ⁻² source model are different from a typical Kostrov's solution. The critical analysis of the kinematic k ⁻² source led us to identify the Fourier decomposition of the slip to be responsible for these difficulties, and to propose a new recombination scheme. It consists of adding a positive correction to the Fourier slip components. The slip is described as the sum of positive contributions at various scales. The SVFs modelled using this new scheme are greatly improved. Moreover, through several parametrical analyses performed to qualify this new approach, we show that the SVF are corrected while preserving the essential quality of the k ⁻² modelling, that is, the ω² spectral shape and C _d apparent directivity of the synthetic accelerograms. Strong ground motion modelling in the near-fault region was made and numerical ground motion parameters were compared to the empirical relationships. We show that predicted peak ground motion is consistent with near-source attenuation laws.     

Anatomy of a small intraplate earthquake: a dissection of its rupture characteristics using regional data

The magnitude m _bLg 5.0 Mont-Laurier earthquake of 1990 October 19, in Quebec, Canada, was one of the largest to have occurred in eastern North America during the past decade. High-frequency ground motions recorded on regional network instruments exceeded values anticipated for an event of its size by a factor of 3. A commonly favoured explanation for the discrepancy is that the source was a rare 'high-stress' event. In this paper, detailed fault-slip models are derived to fit waveform and spectral characteristics of the regional data. The results establish that the effective rupture stress was normal (about 100 bars), that the fault rupture developed asymmetrically, and that the average slip time for points inside the rupture area (approx. 0.1 s) was significantly less than that associated with the standard Brune (1970) source spectral model. The rupture area developed in at least four distinct episodes, each extending the previously ruptured area. Taken together with similar results for the m _bLg 6.5 Saguenay earthquake of 1988 November, the results indicate that a widely used assumption in hazard analyses, that earthquake spectra are adequately represented by the standard Brune spectral model, is unreliable for the interpretation and prediction of strong ground motion.     

High-resolution structures of the Landers fault zone inferred from aftershock waveform data

High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.     

Determination of the fault plane using a single near-field seismic station with a finite-dimension source model

We explore the possibility of determining the actual fault plane of an earthquake from the inversion of near-source displacement seismograms of one station when a finite-dimension source is used instead of a point source model and when the complete displacement is taken into account, including near-field waves. Tests on synthetic seismograms and real data recorded at local distances show that this is possible even with a single, three-component station. A single accelerogram available for the Erzincan, Turkey, 1992 March 13, M _s = 6.8 earthquake is inverted and the solution found is compatible with other seismological studies and with the mechanism expected for the North Anatolian Fault.