Single-station polarization analysis applied to seismic wavefields: A tutorial

Abstract The polarization properties of seismic wavefields recorded by triaxial (three-component) stations can been exploited for event detection, seismic direction finding, and wavefield filtering. Analysis tools have been derived by extending conventional electromagnetic polarization theory to transient seismic wavelets. The basis of this class of polarization processing is a (complex) covariance matrix formed over a window of data to characterize single polarized events within random noise. An eigendecomposition of the (complex) covariance matrix provides eigenvalues and eigenvectors to describe the degree and properties of polarization as well as to formulate polarization filters. Different modes of polarized events require different analysis domains. Rectilinearly polarized events can be examined using a real-valued covariance matrix. In contrast, elliptically polarized events (Rayleigh waves) require the analysis of the analytic signal and a complex-valued covariance matrix to describe the elliptical motion. The success of the polarization measures depends on the signal-to-noise ratio. Furthermore, if multiple events interfere, the basic assumptions used for single-station polarization analysis are violated and the method breaks down. The choice of analysis window length controls the resolution and success of the analysis and depends on the dataset being examined. Ideally, the window should be as long as possible but must not include two (or more) separate events. Synthetic datasets, constructed for different modes of events and different noise levels, as well as real data from a seismic reflection and vertical seismic profiling (VSP) survey together with recordings of an earthquake are used to illustrate the benefits and limitations of polarization processing.

[1]  Mark A. Meadows,et al.  Shear-wave polarizations and subsurface stress directions at Lost Hills field , 1991 .

[2]  David B. Harris,et al.  Comparison of the direction estimation performance of high-frequency seismic arrays and three-component stations , 1990, Bulletin of the Seismological Society of America.

[3]  S. P. Jarpe,et al.  Performance of high-frequency three-component stations for azimuth estimation from regional seismic phases , 1991 .

[4]  John E. Vidale,et al.  Complex polarization analysis of particle motion , 1986 .

[5]  Stewart Greenhalgh,et al.  Controlled Direction Reception Filtering of P- and S-Waves In τ-P Space , 1990 .

[6]  Bing Zhou,et al.  An analytical treatment of single station triaxial seismic direction finding , 2005 .

[7]  J. C. Samson,et al.  Pure states, polarized waves, and principal components in the spectra of multiple, geophysical time-series , 1983 .

[8]  S. Ronen,et al.  Ground Roll Detection and Attenuation by 3C Polarization Analysis , 2005 .

[9]  E. S. Husebye,et al.  Wavefield decomposition using ML-probabilities in modelling single-site 3-component records , 1988 .

[10]  S. Greenhalgh,et al.  Effect of coherent noise on single-station direction of arrival estimation , 2008 .

[11]  E. A. Flinn.  Signal analysis using rectilinearity and direction of particle motion , 1965 .

[12]  J. Achenbach.  Wave propagation in elastic solids , 1962 .

[13]  Thomas J. Owens,et al.  The TauP Toolkit: Flexible Seismic Travel-Time and Raypath Utilities , 1999 .

[14]  Kelly B. Knowlton,et al.  Polarization measurement uncertainty on three-component VSP , 1996 .

[15]  Lion Krischer,et al.  ObsPy: a bridge for seismology into the scientific Python ecosystem , 2015 .

[16]  A. Brink,et al.  Coda polarization properties of high-frequency microseismic events , 1988 .

[17]  E. I. Gal'perin,et al.  The Polarization Method of Seismic Exploration , 1983 .

[18]  A. Jurkevics.  Polarization analysis of three-component array data , 1988 .

[19]  Ram Sriram,et al.  Recording and processing vector wave field data: A review of the 1989 SEG Summer Research Workshop , 1990 .

[20]  W. Banerdt,et al.  Verifying single-station seismic approaches using Earth-based data: Preparation for data return from the InSight mission to Mars , 2015 .

[21]  Emil Wolf,et al.  Principles of Optics: Contents , 1999 .

[22]  F. Ardhuin,et al.  Polarized Earth's ambient microseismic noise , 2011 .

[23]  N. Magotra,et al.  Seismic event detection and source location using single-station (three-component) data , 1987 .

[24]  Gemma Musacchio,et al.  Polarization filter with singular value decomposition , 2001 .

[25]  J. Chiu,et al.  Polarization analysis of high-frequency, three-component seismic data , 1991, Bulletin of the Seismological Society of America.

[26]  Stewart Greenhalgh,et al.  6-C polarization analysis using point measurements of translational and rotational ground-motion: theory and applications , 2018 .

[27]  E. R. Kanasewich,et al.  Time sequence analysis in geophysics , 1973 .

[28]  A. Benhama,et al.  STUDY AND APPLICATIONS OF SPATIAL DIRECTIONAL FILTERING IN THREE‐COMPONENT RECORDINGS1 , 1988 .

[29]  James Gaiser.  3C Seismic and VSP: Converted waves and vector wavefield applications , 2016 .

[30]  Robert B. Herrmann,et al.  Ground roll: Rejection using polarization filters , 1990 .

[31]  E. Kragh,et al.  Ground Roll and Polarization , 1995 .

[32]  Ernest R. Kanasewich,et al.  Enhancement of Teleseismic Body Phases with a Polarization Filter , 1970 .

[33]  G. M. Jackson,et al.  Principal component transforms of triaxial recordings by singular value decomposition , 1991 .

[34]  Stewart Greenhalgh,et al.  Orientation of a downhole triaxial geophone , 1995 .

[35]  N. Magotra,et al.  Single-station seismic event detection and location , 1989 .

[36]  O. Nuttli,et al.  The effect of the earth's surface on the S wave particle motion , 1961 .

[37]  M. J. Rutty,et al.  14. Correlation Using Triaxial Data from Multiple Stations in the Presence of Coherent Noise , 1999 .

[38]  Johan O. A. Robertsson,et al.  Automated, six-component, single-station ground-roll identification and suppression by combined processing of translational and rotational ground motion , 2017 .

[39]  F. Scherbaum,et al.  Polarization analysis in the wavelet domain based on the adaptive covariance method , 2007 .

[40]  Stefanie Donner,et al.  Advances in 6C seismology: Applications of combined translational and rotational motion measurements in global and exploration seismology , 2018 .