Preliminary reference calibrating functions for body-wave magnitudes: refracted P-waves

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


2020 ◽  
Author(s):  
Pei-Ru Jian ◽  
Ban-Yuan Kuo

<p>Seismic attenuation accompanying the velocity structures demonstrates the variations of the physical and chemical properties of the earth. The t* measurement using the seismic body wave spectrum, however, typically encounters the trade-off of corner frequency, t*, and site effect. Ko et al, [2012] proposed the cluster event method (CEM) that reduced the model parameter numbers by grouping the spatial-closed enough events for those traveling to each station along the adjacent paths and sharing one t*. Yet, the site effects among different stations collected in the same cluster bring the challenges on fitting all spectrum. We adapt the cluster strategy to group multiple nearby events recorded by one station only. Moreover, the new iterative CEM algorithm includes both the spectrum and spectral ratio data which provide constraints on seismic moments and corner frequencies of each earthquake inside the cluster, respectively. The final t* and corner frequencies are determined again by including the side effects which are averaging from spectrum residuals in the initial CEM stage. We applied the iterative CEM for earthquakes recorded at dense deployed F-net and Hi-net by NIED in the Tohoku area, Japan. The multitaper spectrums are retrieved from direct P waves with coda wavetrains tapered. Combining the spectral ratio and spectrum data with proper weightings, our new approach increases the stability of t* measurements contributed from better constrains on the corner frequency estimations.</p>


1972 ◽  
Vol 62 (5) ◽  
pp. 1183-1193 ◽  
Author(s):  
F. A. Dahlen

Abstract The effect of a homogeneous anisotropic initial stress on the propagation of infinitesimal amplitude elastic body waves in a perfectly elastic, homogeneous medium is investigated. If the medium is inherently isotropic in the reference configuration and if the magnitude τ0 of the deviatoric part of the initial static stress is small compared to the rigidity μ of the medium, then the apparent body-wave velocities of P waves are unaffected by the initial stress to first order in τ0/μ. The apparent body-wave velocities of S waves are rendered anisotropic to first order, and this effect is described explicitly. It is concluded that the direct effect of an anisotropic initial stress cannot contribute appreciably to the observed velocity anisotropy of horizontally propagating P waves in the oceanic upper mantle. Those observations require an inherent elastic anisotropy of the oceanic uppermantle material.


1982 ◽  
Vol 72 (4) ◽  
pp. 1111-1128
Author(s):  
A. Deschamps ◽  
Y. Gaudemer ◽  
A. Cisternas

abstract We present a study of the El Asnam, Algeria, earthquake of 10 October 1980 from a large collection of long-period surface and body wave records. The focal mechanism of the main shock is well constrained by the P-wave first motions at teleseismic distances and field observations: it was a thrust event on a plane striking N45°E, and a dip angle of 54° to the north. It had a small component of left-lateral motion (λ = 83°) (Ouyed et al., 1981; Gaudemer et al., 1981). This earthquake was very well recorded on WWSSN stations and on GDSN and IDA digital stations, with a good azimuthal distribution. From these records, we confirm the focal solution and obtain a seismic moment Mo = 5 × 1026 dyne-cm. The P-wave seismograms indicate a complex source. We show that it is not sufficient to model the source by a multiple event, but it is also necessary to include a propagation effect in order to explain accurately the waveform. With assumptions based on field observation of the surface breaks, we model the P waves, including two discontinuities of the propagation, along the fault plane and obtain a good agreement of the waveform.


1975 ◽  
Vol 65 (3) ◽  
pp. 693-709 ◽  
Author(s):  
Otto W. Nuttli ◽  
So Gu Kim

abstract Body-wave magnitudes, mb, and surface-wave magnitudes, MS, were determined for approximately 100 Eurasian events which occurred during the interval August through December 1971. Body-wave magnitudes were determined from 1-sec P waves recorded by WWSSN short-period, vertical-component seismographs at epicentral distances greater than 25°. Surface-wave magnitudes were determined from 20-sec Rayleigh waves recorded by long-period, vertical-component WWSSN and VLPE seismographs. The earthquakes had mb values ranging from 3.6 to 5.7. Of 96 presumed earthquakes studied, 6 lie in or near the explosion portion of an mb:MS plot. The explosion mb:MS curve was obtained from seven Eurasian events which had mb values ranging from 5.0 to 6.2 and MS values from 3.2 to 5.1. All six anomalous earthquakes were located in the interior of Asia, in Tibet, and in Szechwan and Sinkiang provinces of China. In general, oceanmargin earthquakes were found to have more earthquake-like mb:MS values than those occurring in the continental interior. Neither focal depth nor focal mechanism can explain the anomalous events.


1964 ◽  
Vol 54 (6A) ◽  
pp. 1981-1996 ◽  
Author(s):  
John Dowling ◽  
Otto Nuttli

abstract Velocities within the earth can be determined from body wave time-distance (T-D) data by the Herglotz-Wiechert method provided the velocity does not decrease too rapidly with depth. Until the present time, the properties of T-D curves for rapid decreases of velocity with depth have been considered only qualitatively. This paper presents a technique for calculating a T-D curve for any velocity distribution, including continuous and discontinuous increases and decreases of velocity with depth. Some properties of T-D curves are quantitatively studied by systematically varying the characteristics of a single model and noting the corresponding variations in the calculated T-D curves. From this it is concluded that a significant low-velocity channel may not be evidenced by a shadow zone but rather by an overlapping of two distinct branches of the T-D curve. It is further concluded that the presence of a shadow zone implies a very gentle velocity gradient below the low-velocity channel. By fitting a calculated T-D curve to observed data one can determine velocity as a function of depth even when the velocity decreases rapidly with depth, when a low-velocity channel exists. Observed T-D data for two underground nuclear explosions (gnome and bilby) measured in four different azimuths were fitted with T-D curves calculated for assumed velocity distributions. It is concluded that these data can be satisfied by a low-velocity channel for P waves in the upper mantle. The character of this channel (depth, thickness and velocity) was determined in each azimuth. The depth to its top was shallow (70 ± km) in the western U.S. and deep (125 ± km) in the eastern U.S. The velocity gradient below the channel is sharp enough to produce no prominent shadow zones. There are significant lateral changes in upper mantle velocities in the western U. S.


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.


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.


Geophysics ◽  
1966 ◽  
Vol 31 (3) ◽  
pp. 524-535 ◽  
Author(s):  
Karl Fuchs

A numerical method for the synthesis of seismograms for body wave propagation in solid wedges is presented. The method is based on the superposition of multiple reflections arising from the entrance of a plane primary wave. Therefore the method is restricted to that part of the time domain where no diffracted waves from the wedge axis occur. In spite of this restriction, dispersion of body waves in wedges can well be studied by this method. Seismograms have been synthesized which show the dispersion of a primary p‐signal propagating in a solid 10‐degree and a 5‐degree wedge with free boundaries. For wedge angles less than 10 degrees the signal front (to be distinguished from the wavefront) suddenly decreases its velocity from that in the infinite medium to about that of the plate wave as the signal approaches the wedge axis. Simultaneously in this transition zone a decrease of the dominant period of the interference signal occurs. These observations are concordant with previous model studies. Particle motion diagrams disclose elliptical polarization of the interference signal in the neighborhood of the wedge axis; the polarization changes its sense from prograde to retrograde on passing through the transition zone.


Author(s):  
Tormod Kværna ◽  
Steven J. Gibbons ◽  
Sven Peter Näsholm

AbstractThe detection and location capability of the International Monitoring System for small seismic events in the continental and oceanic regions surrounding the Sea of Japan is determined mainly by three primary seismic arrays: USRK, KSRS, and MJAR. Body wave arrivals are coherent on USRK and KSRS up to frequencies of around 4 Hz and classical array processing methods can detect and extract features for most regional signals on these stations. We demonstrate how empirical matched field processing (EMFP), a generalization of frequency-wavenumber or f-k analysis, can contribute to calibrated direction estimates which mitigate bias resulting from near-station geological structure. It does this by comparing the narrowband phase shifts between the signals on different sensors, observed at a given time, with corresponding measurements on signals from historical seismic events. The EMFP detection statistic is usually evaluated as a function of source location rather than slowness space and the size of the geographical footprint valid for EMFP templates is affected by array geometry, the availablesignal bandwidth, and Earth structure over the propagation path. The MJAR arrayhas similar dimensions to KSRS but is sited in far more complex geology which results in poor parameter estimates with classical f-k analysis for all signals lacking energy at 1 Hz or below. EMFP mitigates the signal incoherence to some degree but the geographical footprint valid for a given matched field template on MJAR is very small. Spectrogram beamforming provides a robust detection algorithm for high-frequency signals at MJAR. The array aperture is large enough that f-k analysis performed on continuous AR-AIC functions, calculated from optimally bandpass-filtered signals at the different sites, can provide robust slowness estimates for regional P-waves. Given a significantly higher SNR for regional S-phases on the horizontal components of the 3-component site of MJAR, we would expect incoherent detection and estimation of S-phases to improve with 3-component sensors at all sites. Given the diversity of the IMS stations, and the diversity of the methods which provide optimal results for a given station, we advocate the development of seismic processing pipelines which can process highly heterogeneous inputs to help associate characteristics of the incoming signals with physical events.


Sign in / Sign up

Export Citation Format

Share Document