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.     

Wavefield smoothing and the effect of rough velocity perturbations on arrival times and amplitudes

Geometric ray theory is an extremely efficient tool for modelling wave propagation through heterogeneous media. Its use is, however, only justified when the inhomogeneity satisfies certain smoothness criteria. These criteria are often not satisfied, for example in wave propagation through turbulent media. In this paper, the effect of velocity perturbations on the phase and amplitude of transient wavefields is investigated for the situation that the velocity perturbation is not necessarily smooth enough to justify the use of ray theory. It is shown that the phase and amplitude perturbations of transient arrivals can to first order be written as weighted averages of the velocity perturbation over the first Fresnel zone. The resulting averaging integrals are derived for a homogeneous reference medium as well as for inhomogeneous reference media where the equations of dynamic ray tracing need to be invoked. The use of the averaging integrals is illustrated with a numerical example. This example also shows that the derived averaging integrals form a useful starting point for further approximations. The fact that the delay time due to the velocity perturbation can be expressed as a weighted average over the first Fresnel zone explains the success of tomographic inversions schemes that are based on ray theory in situations where ray theory is strictly not justified; in that situation one merely collapses the true sensitivity function over the first Fresnel zone to a line integral along a geometric ray.     

Whole mantle SH velocity model constrained by waveform inversion based on three-dimensional Born kernels

A whole mantle SH velocity model is obtained by using a unique data set and techniques. Body and surface waveforms including major and multi-orbit phases are used as a data set and are inverted by using 3-D Born kernels. The resultant model, SH18CE, reveals the different natures of the two major upwelling systems: the strong low velocity anomalies beneath Africa extend for more than 1000 km from the core–mantle boundary (CMB), whereas those beneath the Pacific are restricted to 300–400 km from the CMB. The results also show the variable natures of stagnant slabs on the 670 discontinuity around Japan: the depths of the strongest high velocity anomalies within the stagnant slabs are different region by region, which is consistent with the detailed delay time tomography model in this area.