Simultaneous body and surface wave retrieval from the seismic ambient field and discrimination from unavoidably arising spurious artifacts

Author(s):  
Ali Riahi ◽  
Zaher-Hossein Shomali ◽  
Anne Obermann ◽  
Ahmad Kamayestani

<p>We simultaneously extract both, direct P-waves and Rayleigh waves, from the seismic ambient noise field recorded by a dense seismic network in Iran. With synthetics, we show that the simultaneous retrieval of body and surface waves from seismic ambient noise leads to the unavoidable appearance of spurious arrivals that could lead to misinterpretations.</p><p>We work with 2 months of seismic ambient noise records from a dense deployment of 119 sensors with interstation distances of 2 km in Iran. To retrieve body and surface waves, we calculate the cross-coherency in low-frequency ranges, i.e. frequencies up to 1.2 Hz, to provide the empirical Green’s functions between each pair of stations. To separate the P and Rayleigh waves, we use the polarization method that also enhances the small amplitude body waves.</p><p>We observe both P and Rayleigh waves with an apparent velocity of 4.9±0.3 and 1.8±0.1 km/s in the studied area, respectively, as well as S or higher mode of Rayleigh waves, with an apparent velocity of 4.1±0.1 km/s. Besides these physical arrivals, we also observe two spurious arrivals with similar amplitudes before/after the P and/or Rayleigh waves that render the discrimination challenging.</p><p>To better understanding these arrivals, we perform synthetic tests. We show that simultaneously retrieving the body and surface waves from seismic ambient noise sources will unavoidably lead to the appearance of superior arrivals in the calculation of empirical Green’s functions.</p>

2021 ◽  
Author(s):  
Korbinian Sager ◽  
Victor Tsai ◽  
Yixiao Sheng ◽  
Florent Brenguier ◽  
Pierre Boué ◽  
...  

The theory of Green's function retrieval essentially requires homogeneously distributed noise sources. Even though these conditions are not fulfilled in nature, low-frequency (<1 Hz) surface waves generated by ocean-crust interactions have been used successfully to image the crust with unprecedented spatial resolution. In contrast to low-frequency surface waves, high-frequency (>1 Hz) body waves have a sharper, more localized sensitivity to velocity contrasts and temporal changes at depth. In general, their retrieval using seismic interferometry is challenging, and recent studies focus on powerful, localized noise sources. They have proven to be a promising alternative but break the assumptions of Green's function retrieval. In this study, we present an approach to model correlations between P waves for these scenarios and analyze their sensitivity to 3D Earth structure. We perform a series of numerical experiments to advance our understanding of these signals and prepare for an application to fault monitoring. In the considered cases, the character of the signals strongly diverges from Green's function retrieval, and the sensitivity to structure has significant contributions in the source direction. An accurate description of the underlying physics allows us to reproduce observations made in the context of monitoring the San Jacinto Fault in California using train-generated seismic waves. This approach provides new perspectives for detecting and localizing temporal velocity changes previously unnoticed by commonly exploited surface-wave reconstructions.


Geophysics ◽  
1946 ◽  
Vol 11 (1) ◽  
pp. 1-9 ◽  
Author(s):  
C. Y. Fu

By taking the apparent velocity along the boundary as the parameter instead of the angle of incidence, the equations for the different wave amplitudes may be put in more symmetrical forms. In this way, it is more convenient to discuss both the body waves and the Rayleigh waves at the same time. A difficulty in the plotting of the square root of the wave intensity against the angles is also discussed. When the reflection or refraction coefficient is not real, the meaning of the intensity, as obtained by squaring the absolute value of the latter quantity, needs clarification.


2021 ◽  
Author(s):  
Mehrdad Fotouhimehr ◽  
Elham Shabani

&lt;p&gt;Knowledge about seismic ambient noise wavefield through decomposition into different participant waves is of special importance in geophysical studies. In this study, WaveDec technique (Maran&amp;#242; et al., 2012) as an array statistical signal processing technique was used to decompose seismic ambient noise wavefield and to estimate wavefield parameters. In this method, the measurements from all components of stations and parameters of interest are modeled jointly which leads to significant improvement in extracting characteristics of surface waves. Considering the contribution of both Love and Rayleigh waves in the wavefield, the method estimates the desired parameters including amplitude, phase, azimuth, wave number and the ellipticity angle (for the Rayleigh wave only) based on the Maximum Likelihood Estimation method. One of the main characteristic of WaveDec is estimating the ellipticity angle of Rayleigh waves. This is very beneficial in determining retrograde and prograde particle motion and also in mode distinction.&lt;/p&gt;&lt;p&gt;In the WaveDec algorithm, the Truncated Newton method is used to optimize likelihood functions with respect to wavefield parameters. Furthermore, Bayesian Information Criterion (BIC) is used to select the best model and wave type determination (Rayleigh, Love, body wave or noise). Regarding a group of generated models for different wave types, the one with the smallest BIC is chosen.&lt;/p&gt;&lt;p&gt;We examined consistency of WaveDec algorithm by applying different numerical optimization methods; Truncated Newton, L-BFGS-B quasi-Newton and simplex-based Nelder-Mead methods. Furthermore, different model selection criteria; BIC, Akaike Information Criterion (AIC) and Hannan&amp;#8211;Quinn Information Criterion (HQC) were examined to study the quality of generated models. They possess different penalty terms to avoid overfitting the models on data. All possible pairs of optimization methods and model selection criteria were utilized and replaced in WaveDec algorithm. In order to compare the resultant dispersion curves of surface waves and ellipticity angle curves of Rayleigh waves, SESAME model M2.1 synthetic data and some seismic ambient noise measurements in Colfiorito basin in Italy (Array B) were analyzed.&lt;/p&gt;


2021 ◽  
Vol 8 ◽  
Author(s):  
Takashi Tonegawa ◽  
Toshinori Kimura ◽  
Eiichiro Araki

Ambient noise correlation is capable of retrieving waves propagating between two receivers. Although waves retrieved using this technique are primarily surface waves, the retrieval of body waves, including direct, refracted, and reflected waves, has also been reported from land-based observations. The difficulty of body wave extraction may be caused by large amplitudes and little attenuation of surface waves excited by microseisms, indicating that body wave extraction using seafloor records is very challenging because microseisms are generated in ocean areas and large amplitudes of surface waves are presumably observed at the seafloor. In this study, we used a unique dataset acquired by dense arrays deployed in the Nankai subduction zone, including a permanent cabled-network of 49 stations, a borehole sensor, and 150 temporary stations, to attempt to extract near-field body waves from ambient seafloor noise observed by multivariate sensors of broadband and short-period seismometers, differential pressure gauges (DPGs), and hydrophones. Our results show that P waves are extracted only in the DPG-record correlations at a frequency of 0.2–0.5 Hz, which can be seen up to a separation distance of two stations of 17 km with an apparent velocity of 3.2 km/s. At 1–3 Hz, P waves are observed only in the vertical-record correlations up to a separation distance of 11 km with an apparent velocity of 2.0 km/s. These velocity differences reflect the vertical velocity gradient of the accretionary prism, because the P waves at low frequencies propagate at relatively long distances and therefore the turning depth is greater. Moreover, the long-period and short-period P waves are observed at the slope and flat regions on the accretionary prism, respectively. To investigate the retrieved wavefield characteristics, we conducted a two-dimensional numerical simulation for wave propagations, where we located single sources at the sea surface above the flat and slope bathymetry regions. Based on our observations and simulations, we suggest that the retrieval of near-field body waves from ambient seafloor noises depends on the relative amplitudes of P and other surface waves in the ambient noise wavefield, and those are controlled by the subseafloor velocity structure, seafloor topography, and water depth.


Biology ◽  
2021 ◽  
Vol 10 (7) ◽  
pp. 672
Author(s):  
Bruce A. Young ◽  
Skye Greer ◽  
Michael Cramberg

In the viper boa (Candoia aspera), the cerebrospinal fluid (CSF) shows two stable overlapping patterns of pulsations: low-frequency (0.08 Hz) pulses with a mean amplitude of 4.1 mmHg that correspond to the ventilatory cycle, and higher-frequency (0.66 Hz) pulses with a mean amplitude of 1.2 mmHg that correspond to the cardiac cycle. Manual oscillations of anesthetized C. aspera induced propagating sinusoidal body waves. These waves resulted in a different pattern of CSF pulsations with frequencies corresponding to the displacement frequency of the body and with amplitudes greater than those of the cardiac or ventilatory cycles. After recovery from anesthesia, the snakes moved independently using lateral undulation and concertina locomotion. The episodes of lateral undulation produced similar influences on the CSF pressure as were observed during the manual oscillations, though the induced CSF pulsations were of lower amplitude during lateral undulation. No impact on the CSF was found while C. aspera was performing concertina locomotion. The relationship between the propagation of the body and the CSF pulsations suggests that the body movements produce an impulse on the spinal CSF.


2015 ◽  
Author(s):  
Rabah Bensalem* ◽  
Djamal Machane ◽  
Jean-Luc Chatelain ◽  
Mohamed Djeddi ◽  
Hakim Moulouel ◽  
...  

1967 ◽  
Vol 57 (1) ◽  
pp. 55-81
Author(s):  
E. J. Douze

abstract This report consists of a summary of the studies conducted on the subject of short-period (6.0-0.3 sec period) noise over a period of approximately three years. Information from deep-hole and surface arrays was used in an attempt to determine the types of waves of which the noise is composed. The theoretical behavior of higher-mode Rayleigh waves and of body waves as measured by surface and deep-hole arrays is described. Both surface and body waves are shown to exist in the noise. Surface waves generally predominate at the longer periods (of the period range discussed) while body waves appear at the shorter periods at quiet sites. Not all the data could be interpreted to define the wave types present.


1991 ◽  
Vol 81 (5) ◽  
pp. 1900-1922
Author(s):  
Arthur Frankel ◽  
Susan Hough ◽  
Paul Friberg ◽  
Robert Busby

Abstract A small aperture (≈300 m), four-station array was deployed in Sunnyvale, California for 5 days to record aftershocks of the Loma Prieta earthquake of October 1989. The purpose of the array was to study the seismic response of the alluvium-filled Santa Clara Valley and the role of surface waves in the seismic shaking of sedimentary basins. Strong-motion records of the Loma Prieta mainshock indicate that surface waves produced the peak velocities and displacements at some sites in the Santa Clara Valley. We use the recordings from the dense array to determine the apparent velocity and azimuth of propagation for various arrivals in the seismograms of four aftershocks with magnitudes between 3.6 and 4.4. Apparent velocities are generally observed to decrease with increasing time after the S wave in the seismograms. Phases arriving less than about 8 sec after the S wave have apparent velocities comparable to the S wave and appear to be body waves multiply reflected under the receiver site or reflected by crustal interfaces. For times 10 to 30 sec after the direct S wave, we observe long-period (1 to 6 sec) arrivals with apparent velocities decreasing from 2.5 to 0.8 km / sec. We interpret these arrivals to be surface waves and conclude that these surface waves produce the long duration of shaking observed on the aftershock records. Much of the energy in the 40 sec after the S-wave is coming approximately from the direction of the source, although some arrivals have backazimuths as much as 60° different from the backazimuths to the epicenters. Two of the aftershocks show arrivals coming from 30 to 40° more easterly than the epicenters. This energy may have been scattered from outcrops along the southeastern edge of the basin. In contrast, the deepest aftershock studied (d = 17 km) displays later arrivals with backazimuths 30 to 40° more westerly than the epicenter. A distinct arrival for one of the aftershocks propagates from the southwest, possibly scattered from the western edge of the basin. Synthetic seismograms derived from a plane-layered crustal model do not produce the long-period Love waves observed in the waveforms of the ML 4.4 aftershock. These Love waves may be generated by the conversion of incident S waves or Rayleigh waves near the edge of the basin.


1967 ◽  
Vol 57 (1) ◽  
pp. 83-90
Author(s):  
J. A. Hudson ◽  
L. Knopoff

abstract The two-dimensional problems of the scattering of harmonic body waves and Rayleigh waves by topographic irregularities in the surface of a simplified model of the earth are considered with especial reference to the processes of P-R, SV-R and R-R scattering. The topography is assumed to have certain statistical properties; the scattered surface waves also have describable statistical properties. The results obtained show that the maximum scattered seismic noise is in the range of wavelengths of the order of the lateral dimensions of the topography. The process SV-R is maximized over a broader band of wavelengths than the process P-R and thus the former may be more difficult to remove by selective filtering. An investigation of the process R-R shows that backscattering is much more important than forward scattering and hence topography beyond the array must be taken into account.


Sign in / Sign up

Export Citation Format

Share Document