Interferometric body-wave retrieval from ambient noise after polarization filtering: application to shallow reflectivity imaging

Geophysics ◽  
2021 ◽  
pp. 1-58
Author(s):  
Deepankar Dangwal ◽  
Michael Behm
Keyword(s):  
Wave Energy ◽  
Ambient Noise ◽  
Body Wave ◽  
Body Waves ◽  
P Wave ◽  
Receiver Array ◽  
Short Period ◽  
Wave Imaging ◽  
Polarization Filtering ◽  
The Impact

Interferometric retrieval of body waves from ambient noise recorded at surface stations is usually challenged by the dominance of surface-wave energy, in particular in settings dominated by anthropogenic activities (e.g., natural resource exploitation, traffic, infrastructure construction). As a consequence, ambient noise imaging of shallow structures such as sedimentary layers remains a difficult task for sparse and irregularly distributed receiver networks. We demonstrate how polarization filtering can be used to automatically extract steeply inclined P-waves from continuous three-component recordings and in turn improves passive body-wave imaging. Being a single-station approach, the technique does not rely on a dense receiver array and is therefore well suited for data collected during surveillance monitoring for tasks such as reservoir hydraulic stimulation, CO_2 sequestration, and wastewater disposal injection. We apply the method on a continuous dataset acquired in the Wellington oilfield (Kansas, US), where local and regional seismicity, and other forms of ambient noise provide an abundant source of both surface- and body-wave energy recorded at 15 short-period receivers. We use autocorrelation to derive the shallow (lt; 1 km) reflectivity structure below the receiver array and validate our workflow and results with well logs and active seismic data. Raytracing analysis and waveform modeling indicates that converted shear waves need to be taken into account for realistic ambient noise body-wave source distributions, as they can be projected on the vertical component and might lead to misinterpretation of the P-wave reflectivity structure. Overall, our study suggests that polarization filtering significantly improves passive body-wave imaging on both autocorrelation and interstation crosscorrelation. It reduces the impact of time-varying noise source distributions and is therefore also potentially useful for time-lapse ambient noise interferometry.

The partition of radiated energy between P and S waves

1984 ◽  
Vol 74 (2) ◽  
pp. 361-376
Author(s):  
John Boatwright ◽  
Jon B. Fletcher
Keyword(s):  
Wave Energy ◽  
Body Wave ◽  
Focal Mechanisms ◽  
The Body ◽  
Body Waves ◽  
P Wave ◽  
Corner Frequency ◽  
S Wave ◽  
S Waves ◽  
Radiated Energy

Abstract Seventy-three digitally recorded body waves from nine multiply recorded small earthquakes in Monticello, South Carolina, are analyzed to estimate the energy radiated in P and S waves. Assuming Qα = Qβ = 300, the body-wave spectra are corrected for attenuation in the frequency domain, and the velocity power spectra are integrated over frequency to estimate the radiated energy flux. Focal mechanisms determined for the events by fitting the observed displacement pulse areas are used to correct for the radiation patterns. Averaging the results from the nine events gives 27.3 ± 3.3 for the ratio of the S-wave energy to the P-wave energy using 0.5 〈Fi〉 as a lower bound for the radiation pattern corrections, and 23.7 ± 3.0 using no correction for the focal mechanisms. The average shift between the P-wave corner frequency and the S-wave corner frequency, 1.24 ± 0.22, gives the ratio 13.7 ± 7.3. The substantially higher values obtained from the integral technique implies that the P waves in this data set are depleted in energy relative to the S waves. Cursory inspection of the body-wave arrivals suggests that this enervation results from an anomalous site response at two of the stations. Using the ratio of the P-wave moments to the S-wave moments to correct the two integral estimates gives 16.7 and 14.4 for the ratio of the S-wave energy to the P-wave energy.

scholarly journals Automatic 3D illumination-diagnosis method for large-N arrays: Robust data scanner and machine-learning feature provider

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. Q13-Q25 ◽  
Author(s):  
Michał Chamarczuk ◽  
Michał Malinowski ◽  
Yohei Nishitsuji ◽  
Jan Thorbecke ◽  
Emilia Koivisto ◽  
...  
Keyword(s):  
Ambient Noise ◽  
Body Wave ◽  
Synthetic Data ◽  
Large Data ◽  
Body Waves ◽  
Support Vector ◽  
Full Data ◽  
Noise Sources ◽  
Data Set ◽  
Receiver Array

The main issues related to passive-source reflection imaging with seismic interferometry (SI) are inadequate acquisition parameters for sufficient spatial wavefield sampling and vulnerability of surface arrays to the dominant influence of the omnipresent surface-wave sources. Additionally, long recordings provide large data volumes that require robust and efficient processing methods. We address these problems by developing a two-step wavefield evaluation and event detection (TWEED) method of body waves in recorded ambient noise. TWEED evaluates the spatiotemporal characteristics of noise recordings by simultaneous analysis of adjacent receiver lines. We test our method on synthetic data representing transient ambient-noise sources at the surface and in the deeper subsurface. We discriminate between basic types of seismic events by using three adjacent receiver lines. Subsequently, we apply TWEED to 600 h of ambient noise acquired with an approximately 1000-receiver array deployed over an active underground mine in Eastern Finland. We develop the detection of body-wave events related to mine blasts and other routine mining activities using a representative 1 h noise panel. Using TWEED, we successfully detect 1093 body-wave events in the full data set. To increase the computational efficiency, we use slowness parameters derived from the first step of TWEED as input to a support vector machine (SVM) algorithm. Using this approach, we detect 94% of the TWEED-evaluated body-wave events indicating the possibility to limit the illumination analysis to only one step, and therefore increase the time efficiency at the price of lower detection rate. However, TWEED on a small volume of the recorded data followed by SVM on the rest of the data could be efficiently used for a quick and robust (real-time) scanning for body-wave energy in large data volumes for subsequent application of SI for retrieval of reflections.

scholarly journals Near-Field Body-Wave Extraction From Ambient Seafloor Noise in the Nankai Subduction Zone

2021 ◽  
Vol 8 ◽  
Author(s):  
Takashi Tonegawa ◽  
Toshinori Kimura ◽  
Eiichiro Araki
Keyword(s):  
Surface Waves ◽  
Ambient Noise ◽  
Body Wave ◽  
Near Field ◽  
Accretionary Prism ◽  
Separation Distance ◽  
Body Waves ◽  
Apparent Velocity ◽  
Short Period ◽  
Nankai Subduction Zone

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.

scholarly journals t* for S waves with a continental ray path

1978 ◽  
Vol 68 (4) ◽  
pp. 1013-1030
Author(s):  
L. J. Burdick
Keyword(s):  
Body Wave ◽  
Body Waves ◽  
P Wave ◽  
Lateral Asymmetry ◽  
Deep Earthquake ◽  
Data Set ◽  
S Waves ◽  
Short Period ◽  
Wave Forms ◽  
Ray Path

abstract The purpose of this study was to determine t* for S waves with ray paths under the continental United States. The data set consists of long- and short-period body waves from the Borrego Mountain earthquake as observed in the northeastern U.S. The P wave forms are dominated by the sP phase and the SH wave forms by the sS. It is assumed that there are no losses in pure compression so that the relative attenuation rate of P and S waves is known. The initial source radiation is determined from the sP phase and the value of tβ* from the spectral content of the S wave. The results indicate that tβ* is 5.2 ± 0.7 sec along this ray path. Long- and short-period body waves from some deep South American events are used to test for lateral asymmetry of the Q distribution under the U.S. No lateral amplitude variation exists in this data, but this result is difficult to correlate with many previous results. The t°* value for a 600-km deep earthquake appears to be about 3 sec. A comparison of these values with values computed from current models of the Earth's Q distribution indicates that the models are slightly too high in Q overall and that more of the total body-wave attenuation occurs above 600 km than is indicated by the models.

A teleseismic body-wave analysis of the May 1980 Mammoth Lakes, California, earthquakes

1983 ◽  
Vol 73 (2) ◽  
pp. 419-434
Author(s):  
Jeffery S. Barker ◽  
Charles A. Langston
Keyword(s):  
Body Wave ◽  
Moment Tensor ◽  
Normal Fault ◽  
Strike Slip ◽  
Body Waves ◽  
P Wave ◽  
Moment Tensors ◽  
Double Couple ◽  
The Moment ◽  
Mammoth Lakes

abstract Teleseismic P-wave first motions for the M ≧ 6 earthquakes near Mammoth Lakes, California, are inconsistent with the vertical strike-slip mechanisms determined from local and regional P-wave first motions. Combining these data sets allows three possible mechanisms: a north-striking, east-dipping strike-slip fault; a NE-striking oblique fault; and a NNW-striking normal fault. Inversion of long-period teleseismic P and SH waves for the events of 25 May 1980 (1633 UTC) and 27 May 1980 (1450 UTC) yields moment tensors with large non-double-couple components. The moment tensor for the first event may be decomposed into a major double couple with strike = 18°, dip = 61°, and rake = −15°, and a minor double couple with strike = 303°, dip = 43°, and rake = 224°. A similar decomposition for the last event yields strike = 25°, dip = 65°, rake = −6°, and strike = 312°, dip = 37°, and rake = 232°. Although the inversions were performed on only a few teleseismic body waves, the radiation patterns of the moment tensors are consistent with most of the P-wave first motion polarities at local, regional, and teleseismic distances. The stress axes inferred from the moment tensors are consistent with N65°E extension determined by geodetic measurements by Savage et al. (1981). Seismic moments computed from the moment tensors are 1.87 × 1025 dyne-cm for the 25 May 1980 (1633 UTC) event and 1.03 × 1025 dyne-cm for the 27 May 1980 (1450 UTC) event. The non-double-couple aspect of the moment tensors and the inability to obtain a convergent solution for the 25 May 1980 (1944 UTC) event may indicate that the assumptions of a point source and plane-layered structure implicit in the moment tensor inversion are not entirely valid for the Mammoth Lakes earthquakes.

scholarly journals Triplicated P-wave measurements for waveform tomography of the mantle transition zone

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 339-354 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch ◽  
T. Nissen-Meyer
Keyword(s):  
Transition Zone ◽  
Epicentral Distance ◽  
Body Wave ◽  
Source Parameters ◽  
Body Waves ◽  
P Wave ◽  
Theoretical Modelling ◽  
Finite Frequency ◽  
Mantle Transition Zone ◽  
Band Limited

Abstract. Triplicated body waves sample the mantle transition zone more extensively than any other wave type, and interact strongly with the discontinuities at 410 km and 660 km. Since the seismograms bear a strong imprint of these geodynamically interesting features, it is highly desirable to invert them for structure of the transition zone. This has rarely been attempted, due to a mismatch between the complex and band-limited data and the (ray-theoretical) modelling methods. Here we present a data processing and modelling strategy to harness such broadband seismograms for finite-frequency tomography. We include triplicated P-waves (epicentral distance range between 14 and 30°) across their entire broadband frequency range, for both deep and shallow sources. We show that is it possible to predict the complex sequence of arrivals in these seismograms, but only after a careful effort to estimate source time functions and other source parameters from data, variables that strongly influence the waveforms. Modelled and observed waveforms then yield decent cross-correlation fits, from which we measure finite-frequency traveltime anomalies. We discuss two such data sets, for North America and Europe, and conclude that their signal quality and azimuthal coverage should be adequate for tomographic inversion. In order to compute sensitivity kernels at the pertinent high body wave frequencies, we use fully numerical forward modelling of the seismic wavefield through a spherically symmetric Earth.

Upper-mantle discontinuity from amplitude data of P′P′ and its precursors

1973 ◽  
Vol 63 (2) ◽  
pp. 587-597
Author(s):  
Ta-Liang Teng ◽  
James P. Tung
Keyword(s):  
Upper Mantle ◽  
Body Waves ◽  
P Wave ◽  
Specific Reference ◽  
Amplitude Data ◽  
Surface Displacements ◽  
Short Period ◽  
Mantle Velocity ◽  
Mantle Discontinuity ◽  
Upper Mantle Discontinuity

abstract Recent observations of P′P′ and its precursors, identified as reflections from within the Earth's upper mantle, are used to examine the structure of the uppermantle discontinuities with specific reference to the density, the S velocity, and the Q variations. The Haskell-Thomson matrix method is used to generate the complex reflection spectrum, which is then Fourier synthesized for a variety of upper-mantle velocity-density and Q models. Surface displacements are obtained for the appropriate recording instrument, permitting a direct comparison with the actual seismograms. If the identifications of the P′P′ precursors are correct, our proposed method yields the following: (1) a structure of Gutenberg-Bullen A type is not likely to produce observable P′P′ upper-mantle reflections, (2) in order that a P′P′ upper-mantle reflection is strong enough to be observed, first-order density and S-velocity discontinuities together with a P-wave discontinuity are needed at a depth of about 650 km, and (3) corresponding to a given uppermantle velocity-density model, an estimate can be made of the Q in the upper mantle for short-period seismic body waves.

WIDE‐BAND EXTRACTION OF MANTLE P WAVES FROM AMBIENT NOISE

Geophysics ◽  
1964 ◽  
Vol 29 (5) ◽  
pp. 672-692 ◽  
Author(s):  
Milo Backus ◽  
John Burg ◽  
Dick Baldwin ◽  
Ed Bryan
Keyword(s):  
Ambient Noise ◽  
Wide Band ◽  
P Wave ◽  
Cumberland Plateau ◽  
P Waves ◽  
Single Output ◽  
Short Period ◽  
Multichannel Filtering ◽  
On Line ◽  
Element Array

The spatial correlation characteristics of ambient short‐period (0.5 to 5 cps) noise at Ft. Sill, Oklahoma, and on the Cumberland Plateau in Tennessee were investigated on “permanent” arrays with 3–4 kilometer diameter. Dominant ambient noise at the two locations is spatially organized, and to first order may be treated as a combination of seismic propagating wave trains. At the Tennessee location noise energy above one cps is dominantly propagating with velocities from 3.5 to 4.5 km/sec, and must be carried in deeply trapped, high‐order modes. Generalized multichannel filtering (Burg) can be used to preserve a large class of mantle P‐wave signals, wide‐band, in a single output trace, while at the same time specifically rejecting ambient noise on the basis of its organization. Results of generalized multichannel filtering applied on‐line at the nineteen‐element array in Tennessee and applied off‐line are discussed.

scholarly journals Inversion of the body waves from the Borrego Mountain earthquake to the source mechanism

1976 ◽  
Vol 66 (5) ◽  
pp. 1485-1499 ◽  
Author(s):  
L. J. Burdick ◽  
George R. Mellman
Keyword(s):  
Body Wave ◽  
High Stress ◽  
Time Function ◽  
The Body ◽  
Body Waves ◽  
Polarization Angle ◽  
Long Period ◽  
Amplitude Data ◽  
Short Period ◽  
Wide Range

abstract The generalized linear inverse technique has been adapted to the problem of determining an earthquake source model from body-wave data. The technique has been successfully applied to the Borrego Mountain earthquake of April 9, 1968. Synthetic seismograms computed from the resulting model match in close detail the first 25 sec of long-period seismograms from a wide range of azimuths. The main shock source-time function has been determined by a new simultaneous short period-long period deconvolution technique as well as by the inversion technique. The duration and shape of this time function indicate that most of the body-wave energy was radiated from a surface with effective radius of only 8 km. This is much smaller than the total surface rupture length or the length of the aftershock zone. Along with the moment determination of Mo = 11.2 ×1025 dyne-cm, this radius implies a high stress drop of about 96 bars. Evidence in the amplitude data indicates that the polarization angle of shear waves is very sensitive to lateral structure.

Crustal structure in the Gangetic Plains of the Indian sub-continent from body waves

1974 ◽  
Vol 64 (3-1) ◽  
pp. 571-579
Author(s):  
R. K. Dube ◽  
J. C. Bhayana
Keyword(s):  
Crustal Structure ◽  
Body Wave ◽  
Bottom Layer ◽  
Body Waves ◽  
P Wave ◽  
Gangetic Plains ◽  
S Waves ◽  
Shear Wave Velocities ◽  
P Wave Velocity ◽  
Peninsular Shield

abstract Crustal structure in the Gangetic Plains of India has been investigated using body-wave data of earthquakes. A three-layered crust has been interpreted, consisting of a top layer of 2.7-km/sec P-wave velocity and 3.7-km thickness, an intermediate layer of 5.64-km/sec velocity and 15.2-km thickness and a bottom layer of 6.49-km/sec velocity and 21.4-km thickness. The average depth of the M discontinuity obtained is 40.3 km. The shear-wave velocities for the Sg, S*, Sn phases are 3.45, 3.85 and 4.61 km/sec, respectively. The velocities of both P and S waves are lower than those obtained for the Peninsular Shield of India.

Sign in / Sign up

Export Citation Format

Share Document