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

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.

2020 ◽  
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>


Geophysics ◽  
2021 ◽  
pp. 1-58
Author(s):  
Deepankar Dangwal ◽  
Michael Behm

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.


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.


1978 ◽  
Vol 68 (1) ◽  
pp. 1-29 ◽  
Author(s):  
Charles A. Langston

abstract Teleseismic P, SV, and SH waves recorded by the WWSS and Canadian networks from the 1971 San Fernando, California earthquake (ML = 6.6) are modeled in the time domain to determine detailed features of the source as a prelude to studying the near and local field strong-motion observations. Synthetic seismograms are computed from the model of a propagating finite dislocation line source embedded in layered elastic media. The effects of source geometry and directivity are shown to be important features of the long-period observations. The most dramatic feature of the model is the requirement that the fault, which initially ruptured at a depth of 13 km as determined from pP-P times, continuously propagated toward the free surface, first on a plane dipping 53°NE, then broke over to a 29°NE dipping fault segment. This effect is clearly shown in the azimuthal variation of both long period P- and SH-wave forms. Although attenuation and interference with radiation from the remainder of the fault are possible complications, comparison of long- and short-period P and short-period pP and P waves suggest that rupture was initially bilateral, or, possibly, strongly unilateral downward, propagating to about 15 km depth. The average rupture velocity of 1.8 km/sec is well constrained from the shape of the long-period wave forms. Total seismic moment is 0.86 × 1026 dyne-cm. Implications for near-field modeling are drawn from these results.


1976 ◽  
Vol 66 (5) ◽  
pp. 1485-1499 ◽  
Author(s):  
L. J. Burdick ◽  
George R. Mellman

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.


1980 ◽  
Vol 70 (3) ◽  
pp. 791-808
Author(s):  
Anne Suteau ◽  
Louis Martel

abstract The transmitted field due to surface waves incident on a local irregularity of a plane-layered medium has been studied. A perturbation method to the first order and the Born approximation can be used if the variations in the thickness of the layers are sufficiently smooth and the wavelengths are long when compared to the size of the irregularities. The spectrum of the perturbed part of the displacement field at the surface is a sum over the surface-wave modes for the regular medium, with an additional term involving the scattered body waves. Numerical computations have been performed for structures composed of a layer overlying a half-space. The contribution of the various modes to the transmitted Love or Rayleigh fields has been studied for several structures. A general method has been obtained to analyze the effect of a complex structure as the superposition of the fields due to simpler ones. When the layer thickness is kept unchanged, the incident mode is not perturbed to the first order. Synthetic seismograms, computed at stations sufficiently close to the irregular region, show how the perturbation of the signal depends on distance. A comparison has been made for Love waves with a finite element method. Both methods give very similar results when the stations are not too close to the irregularities so that the body-wave contribution is negligible. The local phase velocity shows departures from the curves for a regular model.


1978 ◽  
Vol 68 (6) ◽  
pp. 1663-1677
Author(s):  
Stephen H. Hartzell ◽  
James N. Brune ◽  
Jorge Prince

abstract The Acapulco earthquake of October 6, 1974 (mb = 5.0, Ms = 4.75) resulted in 0.5 g accelerations in Acapulco at an epicentral distance of about 35 km. Extrapolation of the peak acceleration to the source region gives a near source acceleration of at least 1.0 g. If the teleseismically estimated source depth of 51 km is assumed, the Acapulco accelerogram must be interpreted as composed of primarily body waves. This assumption yields a moment estimate of 3.3 ×1023 dyne-cm and a stress drop of 1.5 kbar. However, strong evidence indicates that the source depth is only about 1.0 km and that the record is composed mainly of high frequency (1.0 to 4.0 Hz) surface waves. The character of the record is that of a normally dispersed surface wave. The relatively simple form and high acceleration may be attributed to the high rigidity, crystalline rock types in the region. The three component record is fitted by summing the fundamental and first higher mode Rayleigh and Love waves using a model consisting of a single layer over a homogeneous half-space. The results are also checked using a direct wave-number integration program developed by Apsel and Luco. The moment estimate from the surface-wave synthetics is 2.0 ×1023 dyne-cm.


2020 ◽  
Author(s):  
Nima Nooshiri ◽  
Ivan Lokmer ◽  
Chris Bean ◽  
Andrew Bell ◽  
Martin Möllhoff ◽  
...  

<p>Sierra Negra is a basaltic shield volcano in the Galapagos Archipelago (Ecuador) and is the largest of the Galapagos volcanoes. The 2018 eruption was a complex event that included eruptive fissures opening on the northern rim and north-western flank. In this study, we report observations of seismic signals recorded on a temporary dense local network consisting of 14 seismometers and nearby permanent seismic stations, and analyze this data set to retrieve the source mechanisms of moderate pre- and co-eruptive seismic events (body-wave magnitude range of M3.5-5.3). Because of the shallow depths of the seismic events (<2 km) and short source-receiver distances (~1.5-10 km), that are comparable to or shorter than the wavelengths of radiated waves, the effect of near- and intermediate-field terms on dynamic displacements can be significant and hence the far-field assumption may not be well-suited for determining fault-plane solutions. Therefore, we pay special attention on the polarization properties of seismic waves excited at the near-field and intermediate-field ranges, and model and analyze complete displacement wave-fields to determine seismic sources. The source mechanism solutions are also interpreted in light of the volcanic unrest leading to the 2018 eruption, GPS observations, and reported regional centroid moment tensors.</p>


Geophysics ◽  
1979 ◽  
Vol 44 (6) ◽  
pp. 1097-1115 ◽  
Author(s):  
Alfred L. Liaw ◽  
T. V. McEvilly

Frequency(f)‐wavenumber(k) spectra of seismic noise in the bands 1 ⩽ f ⩽ 10 Hz in frequency and |k| ⩽ 35.7 cycles/km in wavenumber, measured at several places in Grass Valley, Nevada, exhibit numerous features which can be correlated with variations in surface geology and sources associated with hot spring activity. Exploration techniques for geothermal reservoirs, based upon the spatial distribution of the amplitude and frequency characteristics of short‐period seismic noise, are applied and evaluated in a field program at this potential geothermal area. A detailed investigation of the spatial and temporal characteristics of the noise field was made to guide subsequent data acquisition and processing. Contour maps of normalized noise level derived from judiciously sampled data are dominated by the hot spring noise source and the generally high noise levels outlining the regions of thick alluvium. Major faults are evident when they produce a shallow lateral contrast in rock properties. Conventional seismic noise mapping techniques cannot differentiate noise anomalies due to buried seismic sources from those due to shallow geologic effects. The noise radiating from a deep reservoir ought to be evident as body waves of high‐phase velocity with time‐invariant source azimuth. A small two‐dimensional (2-D) array was placed at 16 locations in the region to map propagation parameters. The f‐k spectra reveal shallow local sources, but no evidence for a significant body wave component in the noise field was found. With proper data sampling, array processing provides a powerful method for mapping the horizontal component of the vector wavenumber of the noise field. This information, along with the accurate velocity structure, will allow ray tracing to locate a source region of radiating microseisms. In Grass Valley, and probably in most areas of sedimentary cover, the 2–10 Hz microseismic field is predominantly fundamental‐mode Rayleigh waves controlled by the very shallow structure.


Sign in / Sign up

Export Citation Format

Share Document