A surface wave particle motion discrimination process

Estimation of near-surface shear-wave velocities and quality factors using multichannel analysis of surface-wave methods

Using a colocated ring laser and an STS‐2 seismograph, we estimate the ratio of Rayleigh‐to‐Love waves in the secondary microseism at Wettzell, Germany, for frequencies between 0.13 and 0.30 Hz. Rayleigh wave surface acceleration was derived from the vertical component of STS‐2, and Love wave surface acceleration was derived from the ring laser. Surface wave amplitudes are comparable; near the spectral peak about 0.22 Hz, Rayleigh wave amplitudes are about 20% higher than Love wave amplitudes, but outside this range, Love wave amplitudes become higher. In terms of the kinetic energy, Rayleigh wave energy is about 20–35% smaller on average than Love wave energy. The observed secondary microseism at Wettzell thus consists of comparable Rayleigh and Love waves but contributions from Love waves are larger. This is surprising as the only known excitation mechanism for the secondary microseism, described by Longuet‐Higgins (1950), is equivalent to a vertical force and should mostly excite Rayleigh waves.

Comparative analysis of Rayleigh and Love waves detected propagating in the Nobi and Kanto basins during the 2004-, 2007- Chuetsu and 2011 Tohoku earthquakes

SUMMARY A 2-D orthogonal distributed acoustic sensing (DAS) array designed for seismic experiments was buried horizontally beneath the Kafadar Commons Geophysical Laboratory on the Colorado School of Mines campus at Golden, Colorado. The DAS system using straight fibre-optic cables is a cost-efficient technology that enables dense seismic array deployment for long-term seismic monitoring, favouring both earthquake-based and ambient-noise-based surface wave analysis for subsurface characterization. In our study, the horizontally orthogonal DAS array records ambient noise data for a period of about two months from November 2018 to January 2019. During this time, the array also detected seismic signals from an ML3.6 earthquake at Glenwood Springs, Colorado, which exhibit opposite signal polarities in the orthogonal DAS section recordings. We derive the transformation matrix for DAS strain measurements in horizontally orthogonal cables to retrieve both Rayleigh and Love wave dispersion information from the single-component DAS signals using the 2-D multichannel analysis of surface waves method. In addition, ambient noise interferometry is applied to long-term DAS noise recordings. Our theoretical derivation demonstrates that Rayleigh and Love wave Green's functions are coupled in the noise cross-correlation functions (NCFs) of DAS receiver pairs. Stacking NCFs over the horizontally orthogonal DAS array can constructively recover the radial Rayleigh wave component but destructively suppress the Love wave component. The multimodal Monte Carlo inversion of the earthquake-based Rayleigh wave and Love wave dispersion measurements and the noise-based Rayleigh wave measurement reveals a 1-D layered structure that agrees qualitatively with geological surveys of the site. Our study demonstrates that although straight fibre-optic cables lack broadside sensitivity, using appropriate DAS array configuration and seismic array methods can extend the seismic acquisition ability of DAS and enable its application to a broad range of scenarios.

We present a shear velocity model of the crust and uppermost mantle under the Aegean region by simultaneous inversion of Rayleigh and Love waves. The database consists of regional earthquakes recorded by portable broadband three-component digital stations that were installed for a period of 6 months in the broader Aegean region. For each epicentre–station ray path group velocity dispersion curves are measured using appropriate frequency time analysis (FTAN). The dispersion measurements for more than 600 Love wave paths have been used. We have also incorporated previous results for c . 700 Rayleigh wave paths for the study area. The single-path dispersion curves of both waves were inverted to regional group velocity maps for different values of period (6–32 s) via a tomographic method. The local dispersion curves of discrete grid points for both surface waves were inverted nonlinearly to construct 1D models of shear-wave velocity v. depth. In most cases the joint inversion of Rayleigh and Love waves resulted in a single model (from the multiple models compatible with the data) that could interpret both Rayleigh and Love wave data. Around 60 local dispersion curves for both Rayleigh and Love waves were finally jointly inverted. As expected, because of the complex tectonic environment of the Aegean region the results show strong lateral variations of the S-wave velocities for the crust and uppermost mantle. Our results confirm the presence of a thin crust typically less than 28–30 km in the whole Aegean Sea, which in some parts of the southern and central Aegean Sea becomes significantly thinner (20–22 km). In contrast, a large crustal thickness of about 40–45 km exists in western Greece, and the remaining part of continental Greece is characterized by a mean crustal thickness of about 35 km. A significant sub-Moho upper mantle low-velocity zone (LVLmantle) with velocities as low as 3.7 km s −1 , is clearly identified in the southern and central Aegean Sea, correlated with the high heat flow in the mantle wedge above the subducted slab and the related active volcanism in the region. The results obtained results are compared with independent body-wave tomographic information on the velocity structure of the study area and exhibit a generally good agreement, although significant small-scale differences are also identified.

Using ambient noise data from 10 to 50 s period across Alaska, we confirm previous estimates of Rayleigh wave 2 azimuthal anisotropy and present the first estimates of Rayleigh wave 4 and Love wave 2 and 4 azimuthal anisotropy, where is the angle of propagation. As in earlier studies, the fast orientations of Rayleigh wave 2 are mainly parallel to major faults in Alaska at all periods. We also find that on average the fast orientations of Love wave 4 are rotated 45° relative to Rayleigh wave 2, the fast orientation of Rayleigh wave 4 aligns with Love wave 2, and the fast orientation differences of Rayleigh and Love wave 2 range between 0° and 90° with many between 40° to 60°. These observations are consistent with non‐elliptical anisotropy with the ellipticity parameters and considerably smaller than 1. Observations of Love wave 2 and Rayleigh wave 4 reflect strong Rayleigh‐Love coupling, which causes the observed complementary trends with period of the amplitudes of Rayleigh and Love wave 2 and Rayleigh and Love wave 4 Recent theories of Rayleigh‐Love coupling based on a quasi‐degenerate theory allow these observations to be understood and to be used in the future to improve models of the elastic tensor in the crust and mantle.

In this paper, we investigated the propagation of surface waves in a nonhomogeneous rotating fibre-reinforced viscoelastic anisotropic media of higher order of nth order including time rate of strain. The general surface wave speed is derived to study the effect of rotation on surface waves. Particular cases for Stoneley, Love and Rayleigh waves are discussed. The results obtained in this investigation are more general in the sense that some earlier published results are obtained from our result as special cases. Also results for homogeneous media can be deduced from this investigation. For order zero our results are well agreed to fibre-reinforced materials. Also by neglecting the reinforced elastic parameters, the results reduce to well known isotropic medium. It is also observed that, surface waves cannot propagate in a fast rotating medium. Comparison was made with the results obtained in the presence and absence of rotation and parameters for fibre-reinforced of the material medium Numerical results are given and illustrated graphically. The results indicate that the effect of rotation and parameters for fibre-reinforced of the material are very pronounced.

SUMMARY In this paper, we report on the seismogram synthesis of Rayleigh and Love waves in horizontally double-layered elastic half-space. For purpose of analytically extracting surface wave components from full-wavefield, we first derive a compact expression of Green’s function in this medium, from which the excitation formulae of Rayleigh and Love waves are then obtained. Secondly, adopting a set of medium parameters close to those of the Earth, and with Rayleigh poles and Love poles solved in complex domain, we manage to compute synthetic seismograms of the two surface waves separately. Their phase-velocity dispersion and particle motion are also explored. For Rayleigh wave, the senior modes behave rather differently from the fundamental mode. While different Love modes share similar patterns, which are generally in conformity with the classic picture.

This study applies array methods to measure the relative proportions of Love and Rayleigh waves in the ambient vibration wavefield. Information on these properties is of special relevance for frequencies around the horizontal-to-vertical (H/V) spectral amplitude ratio peak. The analysis of H/V curves, a popular technique in site characterisation, commonly assumes that the curves represent the frequency-dependent Rayleigh wave ellipticity. For the detailed interpretation of amplitudes or the inversion of the curves, it is therefore necessary to estimate and correct for the contribution of other wave types to the ambient vibration wavefield. I use available ambient vibration array measurements to determine the relative amount of Love and Rayleigh waves on the horizontal components by frequency-dependent analysis of the main propagation and polarisation directions, with a special emphasis on the H/V peak frequency as determined from the same recordings. Tests with synthetic data demonstrate the feasibility of this approach, at least in the presence of dominant source regions. Analysis of the data from 12 measurements at nine European sites, which include shallow as well as deep locations that span a wide range of impedance contrasts at the sediment-bedrock interface, indicates that the relative contribution of Rayleigh waves varies widely with frequency, from close to 0% to more than 70%. While most data sets show relative Rayleigh wave contributions between 40% and 50% around the H/V peak, there are also examples where Love waves clearly dominate the wavefield at the H/V peak, even for a site with a low impedance contrast. Longer-term measurements at one site indicate temporal variations in the relative Rayleigh wave content between day- and nighttime. Results calculated with the method introduced herein generally compare well with results of modified spatial autocorrelation analysis. These two methods might be used in a complimentary fashion, as both rely on different properties of the ambient vibration wavefield. This study illustrates that it is possible to measure the relative Rayleigh wave content of the noise wavefield from array data. Furthermore, the examples presented herein indicate it is important to estimate this property, as the assumption that there are an equal proportion of Love and Rayleigh waves is not always correct.

Love wave dispersion and associated analytic partial derivatives are theoretically derived using the reflection∕transmission (R∕T) method. The numerical implementation is then applied to linearised inversion of synthetic and field data. For simple cases dominated by the fundamental mode, Love wave sensitivity and inversion stability is higher than the Rayleigh wave dispersion. However, in general, Rayleigh wave inversion converges more rapidly than Love waves, although similar low misfits can be achieved. When assumed interfaces are used, the inversion of fundamental-mode Love wave dispersion of the normally dispersive profile provides a more accurate result, because the Love wave dispersion is independent of the Poisson’s ratio. In more realistic, irregularly dispersive profiles, deep structure is only interpreted using higher-mode Love and Rayleigh wave dispersion. A field test over a shallow geologic fault with coincident Love and Rayleigh wavefields shows the fundamental—mode Love wave dispersion above [Formula: see text] to have at least 10% higher phase velocities than Rayleigh waves. However, at lower frequencies, Love wave dispersion is at least 10% slower. The resulting inverted models show up to 25% difference in shear-wave velocity. This is attributed to transverse isotropy of [Formula: see text] and [Formula: see text] in shallow fluvial sediments and allows improved geological horizon interpretation and soil discrimination over conventional, single-component surface wave inversion. A second test at a seismograph station site shows the Love wavefield less scattered and mode identification is simple, unlike Rayleigh waves over the same line, and the inversion results correlate well to downhole shear-wave velocity logs.

: Seismic surface waves are usually composed of overlapping wave trains representing multi-path propagation. A first task in the analysis of such waves is to identify and separate the various component wave trains so that each can be analyzed separately. Phase-matched filters are a class of linear filters in which the Fourier phase of the filter is made equal to that of a given signal. The authors previously described an iterative technique which can be used to find a phase-matched filter for a particular component of a seismic signal. Application of the filters to digital records of Rayleigh waves allowed multiple arrivals to be identified and removed, and allowed recovery of the complex spectrum of the primary wave train along with its apparent group velocity dispersion curve. A comparable analysis of Love waves presents additional complications. Love waves are contaminated by both Love and Rayleigh multipathing and by primary off-axis Rayleigh energy. In the case of explosions, there is much less energy generated as Love waves than as Rayleigh waves. The applicability of phase-matched filtering to Love waves is demonstrated by its use on earthquakes occurring in the Norwegian Sea and near Iceland and on a nuclear explosion in Novaya Zemlya. Despite severe multipathing in two of the three events, the amplitude and phase of each of the primary Love waves were recovered without significant distortion.

Approximately 2 yr (2010–2011) of continuous seismic records from a subset of the Antarctic component of the Polar Earth Observing Network (POLENET-ANET) seismic network deployed in West Antarctica are used to compute the nine components of the correlation tensor between each pair of stations in the network. Rayleigh wave velocity information from the vertical and radial components was extracted in the form of group and phase velocity dispersion curves, whereas the transverse component provided complimentary Love wave velocity information. The multicomponent Rayleigh wave measurements (ZZ, RR, ZR, and RZ) were averaged and used to infer the measurement uncertainties. The Rayleigh and Love wave group and phase velocities were then regionalized in space using a 2D deterministic tomography. A transect that spans the West Antarctic rift system was extracted from the tomography at individual periods between 7 and 60 s for the four types of surface wave velocities (i.e., Rayleigh and Love phase and group velocities). A transdimensional Bayesian joint inversion algorithm was used to invert these four datasets for a 1D model of isotropic shear-wave velocity versus depth at each point along the transect. In this way, surface wave dispersion curves from multicomponent noise correlations were used to build a 2D isotropic shear-wave velocity model down to ∼55 km depth. In this model, the top of the large low-velocity zone beneath Marie Byrd Land was imaged (up to a 5% decrease in velocity at ∼50 km depth), which provides further evidence for a mantle hot spot beneath the crust that supports the high topography in this region. We also observed a large velocity contrast in the lower crust beneath Marie Byrd Land at a depth where previous long-period seismicity has been observed. This strong contrast occurs more shallow than in previous crustal models, which compared to our model identify a deeper Moho (∼5–10 km deeper) beneath Marie Byrd Land. This new model has implications for interpreting earthquake locations in this region and perhaps necessitates that we revisit past hypocenter estimation studies using updated velocity models for the region.

Surface waves have been extensively studied in earthquake seismology. Surface waves are trapped near an infinitely large surface. The displacements decay exponentially with depth. These waves are also named Rayleigh and Love waves. Surface waves are also used for nondestructive testing of surface defects. Similar waves exist in finite width three-dimensional plates. In this case, displacements are no longer constant in the direction perpendicular to the wave propagation plane. Wave energy could still be trapped near the edge of the three-dimensional plate, and hence the term edge waves. These waves are thus different to the two-dimensional Rayleigh and Love waves. This paper presents a numerical model to study dispersion properties of edge waves in plates. A two-dimensional semi-analytical finite element method is developed, and the problem is closed by a perfectly matched layer adjacent to the edge. The numerical model is validated by comparing with available analytical and numerical solutions in the literature. On this basis, higher order edge waves and mode shapes are presented for a three-dimensional plate. The characteristics of the presented edge wave modes could be used in nondestructive testing applications.

This review covers the propagation of acoustic surface waves (ASW) in periodic structures. The authors discuss the main types of ASW: Rayleigh waves, Gulyaev?Bleustein waves, Love waves. Shear surface waves which can propagate along a periodically uneven surface of an elastic solid are described in detail. The authors rigorously treat the Bragg reflection of ASW from periodic arrays of grooves on the substrate for the cases of normal, oblique, and side incidence. Reciprocal conversion effects between bulk and surface waves in an array whose period is close to the ASW wavelength are discussed, together with such related effects as the interaction between Rayleigh waves and Lamb modes, acoustic wave transmission through a gap in a piezoelectric material, etc. Laser excitation of ASW and magnetostatic wave propagation in periodic structures are also reviewed.

SUMMARY This study deals with surface waves extracted from microseismic noise in the (0.1–0.2 Hz) frequency band with passive seismic-correlation techniques. For directive noise, we explore the concept of passive seismic-noise tomography performed on three-component sensors from a dense seismic network. From the nine-component correlation tensor, a rotation algorithm is introduced that forces each station pair to re-align in the noise direction, a necessary condition to extract unbiased traveltime from passive seismic processing. After rotation is performed, the new correlation tensor exhibits a surface wave tensor from which Rayleigh and Love waves can be separately extracted for tomography inversion. Methodological aspects are presented and illustrated with group-speed maps for Rayleigh and Love waves and ellipticity measurements made on the San Andreas Fault in the Parkfield area, California, USA.