P and S-wave displacements from kinematic dislocation models

1992 ◽  
Vol 82 (4) ◽  
pp. 1910-1926 ◽  
Author(s):  
R. R. Castro ◽  
J. G. Anderson ◽  
J. N. Brune

Abstract Numerical dislocation models based on Haskells (1969) formulation are used to estimate the amount of normal motion necessary to produce the high P/S spectral ratios observed from strong-motion records of the Guerrero, Mexico, subduction zone (Castro et al., 1991). While this depends on the nature of the assumed dislocations, a normal motion with amplitude of less than 10% of the amplitude of shear slip is sufficient to produce P/S values comparable with the observations. We model a planar fault with random patches distributed on the fault plane and a nonplanar fault in which a dilatational jog connects en-echelon fault segments (a structural system proposed by Sibson, 1985, 1989, which introduces normal motions on the fault that depend mainly on the geometry of the fault). For the planar fault model the magnitude of the normal motion is prescribed. Both models introduce complexity in both body-wave displacement time histories, although for the P waves this complexity is accentuated at higher frequencies (f > 1 Hz). For the nonplanar model, a jog with an angle of 10° introduces a normal component of 18% of the slip on the fault.

1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2<τ2(r^)>−1/2, satisfies ωc(P) > ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.


1974 ◽  
Vol 64 (6) ◽  
pp. 1621-1627 ◽  
Author(s):  
J. C. Savage

abstract A comprehensive set of body-wave spectra has been calculated for the Haskell fault model generalized to a circular fault surface. These spectra are used to show that in practice the P-wave corner frequency (ƒp) may exceed the S-wave corner frequency (ƒs) when near-sonic or transonic rupture propagation obtains. The explanation appears to be that in such cases ƒs is so large that it is not identified within the recorded band, but rather a secondary corner is mistaken for ƒs. As a consequence of failing to detect the true asymptotic trend, the high-frequency falloff of the spectrum with frequency is substantially less for S waves than for P waves. This explanation appears to be consistent with the demonstration by Molnar, Tucker, and Brune (1973) that ƒp may exceed ƒs.


1973 ◽  
Vol 63 (6-1) ◽  
pp. 2091-2104 ◽  
Author(s):  
Peter Molnar ◽  
Brian E. Tucker ◽  
James N. Brune

Abstract P- and S-wave spectra of 144 aftershocks 12≦M≦412 of the February 9, 1971 San Fernando earthquake corroborate previous work showing that the corner frequencies for P waves in general are greater than those for S waves. This observation is consistent not only with models that treat earthquakes as volume sources, but also with physically reasonable dislocation models for which (1) the source is approximately equidimensional, (2) both the duration of slip at each point on the fault and the time for the ruptured area to develop are not long compared with the time for seismic waves to cross the ruptured area, and (3) much of the source radiates essentially simultaneously. There may be other physically reasonable dislocation models compatible with the observations. Savage's calculations indicate that models that involve propagating dislocations on long thin faults are not adequate for describing most moderate and small earthquakes studied.


1984 ◽  
Vol 74 (5) ◽  
pp. 1615-1621
Author(s):  
David M. Boore ◽  
John Boatwright

Abstract Averages of P- and S-wave radiation patterns over all azimuths and various ranges of takeoff angles (corresponding to observations at teleseismic, regional, and near distances) have been computed for use in seismological applications requiring average radiation coefficients. Various fault orientations and averages of the squared, absolute, and logarithmic radiation patterns have been considered. Effective radiation patterns combining high-frequency direct and surfacere-flected waves from shallow faults have also been derived and used in the computation of average radiation coefficients at teleseismic distances. In most cases, the radiation coefficients are within a factor of 1.6 of the commonly used values of 0.52 and 0.63 for the rms of P- and S-wave radiation patterns, respectively, averaged over the whole focal sphere. The main exceptions to this conclusion are the coefficients for P waves at teleseismic distances from vertical strike-slip faults, which are at least a factor of 2.8 smaller than the commonly used value.


Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad

Compensation for geometrical spreading is important in prestack Kirchhoff migration and in amplitude versus offset/amplitude versus angle (AVO/AVA) analysis of seismic data. We present equations for the relative geometrical spreading of reflected and transmitted P‐ and S‐wave in horizontally layered transversely isotropic media with vertical symmetry axis (VTI). We show that relatively simple expressions are obtained when the geometrical spreading is expressed in terms of group velocities. In weakly anisotropic media, we obtain simple expressions also in terms of phase velocities. Also, we derive analytical equations for geometrical spreading based on the nonhyperbolic traveltime formula of Tsvankin and Thomsen, such that the geometrical spreading can be expressed in terms of the parameters used in time processing of seismic data. Comparison with numerical ray tracing demonstrates that the weak anisotropy approximation to geometrical spreading is accurate for P‐waves. It is less accurate for SV‐waves, but has qualitatively the correct form. For P waves, the nonhyperbolic equation for geometrical spreading compares favorably with ray‐tracing results for offset‐depth ratios less than five. For SV‐waves, the analytical approximation is accurate only at small offsets, and breaks down at offset‐depth ratios less than unity. The numerical results are in agreement with the range of validity for the nonhyperbolic traveltime equations.


1988 ◽  
Vol 78 (5) ◽  
pp. 1707-1724
Author(s):  
Masayuki Kikuchi ◽  
Yoshio Fukao

Abstract The seismic wave energy is evaluated for 35 large earthquakes by inverting far-field long-period P waves into the multiple-shock sequence. The results show that the seismic wave energy thus obtained is systematically less than that inferred from the Gutenberg-Richter's formula with the seismic magnitude. The difference amounts to one order of magnitude. The results also show that the energy-moment ratio is well confined to a narrow range: 10−6 < ES/Mo < 10−5 with the average of ∼5 × 10−6. This average value is exactly one order of magnitude as small as the energy-moment ratio inferred from the Gutenberg-Richter's formula using the moment magnitude. Comparing the energy-moment ratio with Δσo/2μ, where Δσo and μ are the stress drop and the rigidity, we obtain an empirical relation: ES/Mo ∼ 0.1 × Δσ0/2μ. Such a relation can be interpreted in terms of a subsonic rupture where the energy loss due to cohesion is not negligible to the seismic wave energy.


1995 ◽  
Vol 85 (6) ◽  
pp. 1821-1834
Author(s):  
Toshimi Satoh ◽  
Toshiaki Sato ◽  
Hiroshi Kawase

Abstract We evaluate the nonlinear behavior of soil sediments during strong ground shaking based on the identification of their S-wave velocities and damping factors for both the weak and strong motions observed on the surface and in a borehole at Kuno in the Ashigara Valley, Japan. First we calculate spectral ratios between the surface station KS2 and the borehole station KD2 at 97.6 m below the surface for the main part of weak and strong motions. The predominant period for the strong motion is apparently longer than those for the weak motions. This fact suggests the nonlinearity of soil during the strong ground shaking. To quantify the nonlinear behavior of soil sediments, we identify their S-wave velocities and damping factors by minimizing the residual between the observed spectral ratio and the theoretical amplification factor calculated from the one-dimensional wave propagation theory. The S-wave velocity and the damping factor h (≈(2Q)−1) of the surface alluvial layer identified from the main part of the strong motion are about 10% smaller and 50% greater, respectively, than those identified from weak motions. The relationships between the effective shear strain (=65% of the maximum shear strain) calculated from the one-dimensional wave propagation theory and the shear modulus reduction ratios or the damping factors estimated by the identification method agree well with the laboratory test results. We also confirm that the soil model identified from a weak motion overestimates the observed strong motion at KS2, while that identified from the strong motion reproduces the observed. Thus, we conclude that the main part of the strong motion, whose maximum acceleration at KS2 is 220 cm/sec2 and whose duration is 3 sec, has the potential of making the surface soil nonlinear at an effective shear strain on the order of 0.1%. The S-wave velocity in the surface alluvial layer identified from the part just after the main part of the strong motion is close to that identified from weak motions. This result suggests that the shear modulus recovers quickly as the shear strain level decreases.


1989 ◽  
Vol 166 (1-3) ◽  
pp. 189-203 ◽  
Author(s):  
S.J. Duda ◽  
T.B. Yanovskaya ◽  
E.N. Its ◽  
R. Nortmann
Keyword(s):  

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.


Sign in / Sign up

Export Citation Format

Share Document