scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Enhancement of Teleseismic Body Phases with a Polarization Filter

01 Nov 1970-Geophysical Journal International (Oxford University Press)-Vol. 21, Iss: 2, pp 119-129
TL;DR: In this paper, a time domain polarization filter was used to increase the signal to noise ratio of teleseismic body phases, and the direction of polarization was found from the eigenvector of the largest principle axis.
Abstract: Summary A time domain polarization filter, originally proposed by Flinn, has been modified and is used to increase the signal to noise ratio of teleseismic body phases. Both the rectilinearity and direction of particle motion is obtained from the co-variance matrix for three components of ground motion over a small time interval. The estimate of rectilinearity is found by diagonalizing the matrix and forming a function that involves the ratio of the largest and intermediate principle axes of the matrix. The direction of polarization is found from the eigenvector of the largest principle axis. A set of time-varying operators are then obtained which act as a point by point gain control to modulate the digital seismic records. This non-linear filter is useful for enhancing P, pP, sP, PP, S and other compressional or shear phases. When applied to an array of three-component stations it appears to be possible to identify multiple events in the source function of some earthquakes. 1. Iotroductioo As any seismograph record is always noise contaminated, the detection and interpretation of seismic events requires knowledge of the characteristics of both signal and noise. Elastic body waves, which may be generally separated into P (compressional) and S (shear) phases, can be considered as non-dispersive group arrivals with maximum power in the 0.3 to 10-second period range. Surface Rayleigh and Love waves may be described as dispersive group arrivals with maximum power in the 2-to 100-s period range for earthquakes of moderate magnitude, the observable periods extending up to 57 min for larger teleseismic events. Superimposed on these signals is microseismic background noise as well as signal generated noise. Signal generated noise is the result of multiple reflections and refractions of P and S body waves at crustal interfaces and inhomogeneities and local conversion of body waves to surface waves; this type of noise originates primarily under the recording station. Microseismic noise, which is considered to consist mainly of fundamental and higher mode Rayleigh waves, has been shown to exhibit a sharp peak in the 540 8-s period range (Brune & Oliver 1959). As a result of this sharp peak in the microseismic noise spectrum, frequency bandpass filtering is often employed to improve the signal to noise ratio. Although this type of filtering is very effective in removing microseismic background from both long and short period data, it often cannot distinguish between signal and signal generated noise. Difficulty may still arise in the attempted identification of phases whose frequency characteristics are similar.

Content maybe subject to copyright    Report

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, a technique is presented for polarization analysis of three-component seismic array data, where the polarization ellipse is computed within sliding time windows by solving the eigenproblem for the covariance matrix.
Abstract: A technique is presented for polarization analysis of three-component seismic array data. The process is applied to a large suite of regional events recorded on the three-component sensors in the NORESS array in southern Norway. Polarization properties of the regional seismic phases P,, S,, and Lg are examined in detail. The analysis technique is based on a time-domain algorithm originally proposed by Flinn (1965). The polarization ellipse is computed within sliding time windows by solving the eigenproblem for the covariance matrix. Various attributes characterizing the particle motions are extracted from the motion ellipse. This technique is extended to multiple three-component sensors in an array configuration by averaging covariance matrices for the different sensors. In this case a 1/M reduction in the estimation variance is obtained (M is the number of sensors), when the noise and local scattering effects are uncorrelated. An important feature of this approach is that the phase velocities of coherent wavefronts across the array are not required to a high degree of accuracy. Significant results of the data analysis are the well-defined polarization of P, and Sn waves across the entire short-period band, the source azimuth estimates obtained from P, and Lg motions, and the distinct polarization for S, and Lg waves allowing these phases to be distinguished in most cases.

502 citations

Journal ArticleDOI
TL;DR: While no algorithm was clearly optimal under all source, receiver, path, and noise conditions tested, an STA/LTA algorithm incorporating adaptive window lengths controlled by nonstationary seismogram spectral characteristics was found to provide an output that best met the requirements of a global correlated event-detection and location system.
Abstract: Digital algorithms for robust detection of phase arrivals in the presence of stationary and nonstationary noise have a long history in seismology and have been exploited primarily to reduce the amount of data recorded by data logging systems to manageable levels. In the present era of inexpensive digital storage, however, such algorithms are increasingly being used to flag signal segments in continuously recorded digital data streams for subsequent processing by automatic and/or expert interpretation systems. In the course of our development of an automated, near-real-time, waveform correlation event-detection and location system (WCEDS), we have surveyed the abilities of such algorithms to enhance seismic phase arrivals in teleseismic data streams. Specifically, we have considered envelopes generated by energy transient (STA/LTA), Z -statistic, frequency transient, and polarization algorithms. The WCEDS system requires a set of input data streams that have a smooth, low-amplitude response to background noise and seismic coda and that contain peaks at times corresponding to phase arrivals. The algorithm used to generate these input streams from raw seismograms must perform well under a wide range of source, path, receiver, and noise scenarios. Present computational capabilities allow the application of considerably more robust algorithms than have been historically used in real time. However, highly complex calculations can still be computationally prohibitive for current workstations when the number of data streams become large. While no algorithm was clearly optimal under all source, receiver, path, and noise conditions tested, an STA/LTA algorithm incorporating adaptive window lengths controlled by nonstationary seismogram spectral characteristics was found to provide an output that best met the requirements of a global correlation-based event-detection and location system.

478 citations


Cites background or methods from "Enhancement of Teleseismic Body Pha..."

  • ...To emphasize this feature of the data, the three-component data, (z, n, e), are prefiltered with a polarization filter (e.g., Montalbetti and Kanasewich, 1970; Aster et al., 1990)....

    [...]

  • ...Aster et aL (1990), Magotra et al. (1987, 1989), Montalbetti and Kanasewich (1970), and Flinn (1965), to name a few, constructed detectors and filters based on the eigenvalues and eigenvectors of the three-component covariance matrix....

    [...]

Journal ArticleDOI
TL;DR: In this article, the authors extended the Montalbetti and Kanasewich scheme to analytic three-component seismograms, where the imaginary part of the signal is the Hilbert transform of the real part.
Abstract: Knowledge of particle motion polarization aids in identifying phases on three-component seismograms. The scheme of Montalbetti and Kanasewich (1970) is extended to analytic three-component seismograms, where the imaginary part of the signal is the Hilbert transform of the real part. This scheme has only one free parameter, the length of the time window over which the polarization parameters are estimated, so it can be applied in a routine way to three-component data. The azimuth and dip of the direction of maximum polarization and the degree of elliptical polarization as a function of time for the seismograms are obtained. Polarization analysis of strong motion data from the 1971 San Fernando earthquake aids in the discrimination between wave types, which is important for the understanding of the complicated earthquake-induced shaking observed in basins. Most arrivals are incident on the receivers in the direction of the back-azimuth to the epicenter, which suggests that despite the complicated motions, two-dimensional finite difference methods are sufficient to understand the effect on seismic waves of the Los Angeles and San Fernando basins (Vidale and Helmberger, 1986b).

393 citations


Cites background or methods from "Enhancement of Teleseismic Body Pha..."

  • ...This scheme differs from that of Montalbetti and Kanasewich (1970) in which a running average is subtracted from the seismograms before computing the covariance ....

    [...]

  • ...The real counterpart to the filter described below was first presented by Montalbetti and Kanasewich (1970)....

    [...]

  • ...In addition, and this is particularly true for the scheme of Montalbetti and Kanasewich (1970), the polarization character of the higher frequency energy is more stably estimated, as there are more wavelengths within the time-averaging window....

    [...]

  • ...The RE- MODE filter (Mims and Sax, 1965) and the polarization filter of Montalbetti and Kanasewich (1970) are used to enhance highly polarized portions of the signal automatically....

    [...]

Journal ArticleDOI
TL;DR: In this article, the amplitude and duration of the strong motion velocities recorded across the Los Angeles and San Fernando basins for the period range of 2 to 10 seconds were analyzed.
Abstract: Finite-difference seismograms calculated for the 1971 San Fernando earthquake show strong effects due to lateral variation in sediment thickness in the San Fernando valley and the Los Angeles basin. Using basin structure derived mostly from well logs and teleseismically determined source parameters, two-dimensional SH and P-SV finite-difference calculations reproduce the amplitude and duration of the strong-motion velocities recorded across the Los Angeles and San Fernando basins for the period range of 2 to 10 sec. The edges of basins nearest the seismic source show ground motion amplification up to a factor of 3 over the case without the basin and tend to convert direct shear waves into Love and Rayleigh waves that travel within the basins.

198 citations

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the wave fields generated by Strombolian activity using data from small-aperture seismic arrays deployed on the north flank of Stromboli and data from seismic and pressure transducers set up near the summit crater.
Abstract: The wave fields generated by Strombolian activity are investigated using data from small-aperture seismic arrays deployed on the north flank of Stromboli and data from seismic and pressure transducers set up near the summit crater. Measurements of slowness and azimuth as a function of time clearly indicate that the sources of tremor and explosions are located beneath the summit crater at depths shallower than 200 m with occasional bursts of energy originating from sources extending to a depth of 3 km. Slowness, azimuth, and particle motion measurements reveal a complex composition of body and surface waves associated with topography, structure, and source properties. Body waves originating at depths shallower than 200 m dominate the wave field at frequencies of 0.5–2.5 Hz, and surface waves generated by the surficial part of the source and by scattering sources distributed around the island dominate at frequencies above 2.5 Hz. The records of tremor and explosions are both dominated by SH motion. Far-field records from explosions start with radial motion, and near-field records from those events show dominantly horizontal motion and often start with a low-frequency (1–2 Hz) precursor characterized by elliptical particle motion, followed within a few seconds by a high-frequency radial phase (1–10 Hz) accompanying the eruption of pyroclastics. The dominant component of the near- and far-field particle motions from explosions, and the timing of air and body wave phases observed in the near field, are consistent with a gas-piston mechanism operating on a shallow (<200 m deep), vertical crack-like conduit. Models of a degassing fluid column suggest that noise emissions originating in the collective oscillations of bubbles ascending in the magma conduit may provide an adequate self-excitation mechanism for sustained tremor generation at Stromboli.

187 citations

References
More filters
Book
01 Jan 1968
TL;DR: In this paper, Spectral Analysis and its Applications, the authors present a set of applications of spectral analysis and its application in the field of spectroscopy, including the following:
Abstract: (1970). Spectral Analysis and its Applications. Technometrics: Vol. 12, No. 1, pp. 174-175.

4,220 citations

Journal ArticleDOI
TL;DR: In this article, spectral amplitudes and travel times of seismic body waves are used to determine mantle velocity structures appropriate to distinct structural provinces within the western continental United States, including the Basin and Range.
Abstract: The spectral amplitudes and travel times of seismic body waves are used to determine mantle velocity structures appropriate to distinct structural provinces within the western continental United States. In addition to basic amplitude and time data, travel-time delays and Pn velocity data from other studies are used as constraints in the systematic inversion of the data for mantle structure. The regional structures for the upper mantle determined in this manner show collectively rather sharp zones of transition (high velocity gradients) near 150, 400, 650 km and possibly near 1000 km. Comparatively, the regional structures indicate strong lateral variations in the upper mantle structure down to 150 km and possibly as deep as 200 km. The structures appropriate to the Rocky Mountain and Colorado plateau physiographic provinces show low-velocity zones capped by high-velocity lid zones, with variability in both the lid and the low-velocity zone properties from province to province and within these provinces to a much lesser degree. The mantle properties obtained for the Basin and Range contrast sharply with the plateau and mountain structures, with the lid zone being very thin or absent and abnormally low velocities extending from, or very near, the base of a thin crust to 150 km. The velocity determinations are coupled with estimates of the variation of the intrinsic dissipation function (Q) as a function of depth and frequency. These results show a pronounced low-Q zone corresponding to the average low-velocity zone depth range for the velocity models. The data suggest a frequency-dependent Q, with Q increasing with frequency. In total the results of the study strongly suggest phase transitions in the mantle, including a partially melted region corresponding to the low-velocity zone, the latter being highly variable in its properties over the region studied and strongly correlated with tectonic activity.

295 citations

Journal ArticleDOI
01 Dec 1965

240 citations


"Enhancement of Teleseismic Body Pha..." refers methods in this paper

  • ...Eccentricity, major axis, and angle of inclination of the ellipse from the vertical are displayed for the frequency at which maximum power is arriving and used to provide criterion for the identification of P and S V type motion. Various workers have applied polarization filtering techniques to recorded seismic data for improvement of the signal to noise ratio. Lewis & Meyer (1968) apply a phase filter of the REMODE type as described by Archambeau & Flinn (1965) and originally developed by Mims & Sax (1965) with subsequent work by Griffin (1966a, b) to data recorded during the Early Rise experiment in the summer of 1966....

    [...]

  • ...Various events were recorded at the University of Alberta’s Edmonton, Geophysical Observatory and at two field stations operated 5 and 10 km for the observatory using a broadband digital recording system described by Burke et al. (1970). The three data channels (vertical, east-west horizontal and north-south horizontal) were first bandpass filtered using a recursive zero phase shift digital filter and then rotated into an R-T-2 co-ordinate system so that the radial direction corresponded to the computed great circle azimuth from source to receiver and the transverse component was orthogonal to this....

    [...]

  • ...Eccentricity, major axis, and angle of inclination of the ellipse from the vertical are displayed for the frequency at which maximum power is arriving and used to provide criterion for the identification of P and S V type motion. Various workers have applied polarization filtering techniques to recorded seismic data for improvement of the signal to noise ratio. Lewis & Meyer (1968) apply a phase filter of the REMODE type as described by Archambeau & Flinn (1965) and originally developed by Mims & Sax (1965) with subsequent work by Griffin (1966a, b) to data recorded during the Early Rise experiment in the summer of 1966. Archambeau, Flinn & Lambert (1969) used the same filter to study multiple P phases from NTS explosions. In another application, Basham & Ellis (1969) use a REMODE filter designated as a &Detection (P-D) filter to process P-wave codes of numerous seismic events recorded in western Alberta....

    [...]

  • ...Eccentricity, major axis, and angle of inclination of the ellipse from the vertical are displayed for the frequency at which maximum power is arriving and used to provide criterion for the identification of P and S V type motion. Various workers have applied polarization filtering techniques to recorded seismic data for improvement of the signal to noise ratio. Lewis & Meyer (1968) apply a phase filter of the REMODE type as described by Archambeau & Flinn (1965) and originally developed by Mims & Sax (1965) with subsequent work by Griffin (1966a, b) to data recorded during the Early Rise experiment in the summer of 1966. Archambeau, Flinn & Lambert (1969) used the same filter to study multiple P phases from NTS explosions....

    [...]

Journal ArticleDOI
TL;DR: In this article, the maximum, average, and minimum values of surface particle displacement, velocity, and acceleration of earth noise as a function of period are illustrated in graphical form for periods less than about 5 seconds.
Abstract: Maximum, average, and minimum values of surface particle displacement, velocity, and acceleration of earth noise as a function of period are illustrated in graphical form. For periods less than about 5 seconds the amplitude curves rise rapidly with increasing period. The most prominent feature of the illustration is the sharp peak in the 5- to 8-second period range. There are virtually no data on noise in the range of periods between 20 seconds and the earth tide periods. With the exception of the 10- to 40-second period range, the data used are taken from the existing literature.

98 citations

Journal ArticleDOI
TL;DR: In this article, the maximum, average, and minimum values of surface particle displacement, velocity, and acceleration of earth noise as a function of period are illustrated in graphical form for periods less than about 5 seconds.
Abstract: Maximum, average, and minimum values of surface particle displacement, velocity, and acceleration of earth noise as a function of period are illustrated in graphical form. For periods less than about 5 seconds the amplitude curves rise rapidly with increasing period. The most prominent feature of the illustration is the sharp peak in the 5- to 8-second period range. There are virtually no data on noise in the range of periods between 20 seconds and the earth tide periods. With the exception of the 10- to 40-second period range, the data used are taken from the existing literature.

87 citations