scholarly journals Elastic anisotropy estimation from laboratory measurements of velocity and polarization of quasi-P-waves using laser interferometry

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA83-WA89 ◽  
Author(s):  
Maxim Lebedev ◽  
Andrej Bóna ◽  
Roman Pevzner ◽  
Boris Gurevich
Keyword(s):  
Particle Velocity ◽  
Elastic Anisotropy ◽  
Laser Interferometry ◽  
Elasticity Tensor ◽  
Ultrasonic Waves ◽  
Body Waves ◽  
P Wave ◽  
Surface Element ◽  
Laboratory Measurements ◽  
P Waves

A new method for conducting laboratory measurements of the velocities and polarizations of compressional and shear waves in rock samples uses a laser Doppler interferometer (LDI). LDI can measure the particle velocity of a small (0.03 mm2) element of the surface of the sample along the direction of the laser beam. By measuring the particle velocity of the same surface element in three linearly independent directions and then transforming those velocities to Cartesian coordinates, three orthogonal components of the particle-velocity vector are obtained. Thus, LDI can be used as a localized three-component (3C) receiver of ultrasonic waves, and, together with a piezoelectric transducer as a source, it can simulate a 3C seismic experiment in the laboratory. Performing such 3C measurements at various locations on the surface of the sample produces a 3C seismogram, which can be used to separate the P-wave and two S-waves and to find the polarizations and traveltimes of those waves. Then, the elasticity tensor of the medium can be obtained by minimizing the misfit between measured and predicted polarizations and traveltimes. Computation of the polarizations and traveltimes of body waves inside a sample with a given elasticity tensor is based on the Christoffel equation. The predicted polarizations on the surface then are obtained using the anisotropic Zoeppritz equations. The type of velocity measured (phase or group velocity) depends on the acquisition geometry and the material properties. This is taken into account in the inversion procedure. A “walkaway” laboratory experiment demonstrates the high accuracy of this method.

Elastic velocity anisotropy in the presence of an anisotropic initial stress

1972 ◽  
Vol 62 (5) ◽  
pp. 1183-1193 ◽  
Author(s):  
F. A. Dahlen
Keyword(s):  
Initial Stress ◽  
Elastic Anisotropy ◽  
Body Wave ◽  
Homogeneous Medium ◽  
Body Waves ◽  
P Waves ◽  
S Waves ◽  
Velocity Anisotropy ◽  
First Order ◽  
Wave Velocities

Abstract The effect of a homogeneous anisotropic initial stress on the propagation of infinitesimal amplitude elastic body waves in a perfectly elastic, homogeneous medium is investigated. If the medium is inherently isotropic in the reference configuration and if the magnitude τ0 of the deviatoric part of the initial static stress is small compared to the rigidity μ of the medium, then the apparent body-wave velocities of P waves are unaffected by the initial stress to first order in τ0/μ. The apparent body-wave velocities of S waves are rendered anisotropic to first order, and this effect is described explicitly. It is concluded that the direct effect of an anisotropic initial stress cannot contribute appreciably to the observed velocity anisotropy of horizontally propagating P waves in the oceanic upper mantle. Those observations require an inherent elastic anisotropy of the oceanic uppermantle material.

scholarly journals Comparison of recent S-wave indicating methods

2018 ◽  
Vol 29 ◽  
pp. 00019
Author(s):  
Katarzyna Hubicka ◽  
Jakub Sokolowski
Keyword(s):  
Real Data ◽  
The Body ◽  
Body Waves ◽  
Identification Accuracy ◽  
P Wave ◽  
S Wave ◽  
P Waves ◽  
S Waves ◽  
Mode Decomposition ◽  
Wave Arrival

Seismic event consists of surface waves and body waves. Due to the fact that the body waves are faster (P-waves) and more energetic (S-waves) in literature the problem of their analysis is taken more often. The most universal information that is received from the recorded wave is its moment of arrival. When this information is obtained from at least four seismometers in different locations, the epicentre of the particular event can be estimated [1]. Since the recorded body waves may overlap in signal, the problem of wave onset moment is considered more often for faster P-wave than S-wave. This however does not mean that the issue of S-wave arrival time is not taken at all. As the process of manual picking is time-consuming, methods of automatic detection are recommended (these however may be less accurate). In this paper four recently developed methods estimating S-wave arrival are compared: the method operating on empirical mode decomposition and Teager-Kaiser operator [2], the modification of STA/LTA algorithm [3], the method using a nearest neighbour-based approach [4] and the algorithm operating on characteristic of signals’ second moments. The methods will be also compared to wellknown algorithm based on the autoregressive model [5]. The algorithms will be tested in terms of their S-wave arrival identification accuracy on real data originating from International Research Institutions for Seismology (IRIS) database.

Lattice preferred orientation and seismic anisotropy of chloritoid in subduction zone

2021 ◽  
Author(s):  
Jungjin Lee ◽  
Mainak Mookherjee ◽  
Taehwan Kim ◽  
Haemyeong Jung ◽  
Reiner Klemd
Keyword(s):  
Subduction Zone ◽  
Subduction Zones ◽  
Seismic Anisotropy ◽  
Elastic Anisotropy ◽  
P Wave ◽  
Stability Field ◽  
S Wave ◽  
P Waves ◽  
Hydrous Minerals ◽  
Natural Rock

<p>Subduction zones are often characterized by the presence of strong trench-parallel seismic anisotropy and large delay times. Hydrous minerals, owing to their large elastic anisotropy and strong lattice preferred orientations (LPOs) are often invoked to explain these observations. However, the elasticity and LPO of chloritoid, which is one such hydrous phases relevant in subduction zone settings, is poorly understood. In this study, we measured the LPO of polycrystalline chloritoid in natural rock samples and obtained the LPO-induced seismic anisotropy and evaluated the thermodynamic stability field of chloritoid in subduction zones. The LPO of chloritoid aggregates displayed a strong alignment of the [001] axes subnormal to the rock foliation, with a girdle distribution of the [100] axes and the (010) poles subparallel to the foliation. New elasticity data of single-crystal chloritoid showed a strong elastic anisotropy of chloritoid with 47% for S-waves (V<sub>S</sub>) and 22% for P-waves (V<sub>P</sub>), respectively. The combination of the LPO and the elastic anisotropy of the chloritoid aggregates produced a strong S-wave anisotropy of AV<sub>S</sub> = 18% and a P-wave anisotropy of AV<sub>P</sub> = 10%. Our results indicate that the strong LPO of chloritoid along the hydrated slab-mantle interface and in subducting slabs can influence trench-parallel seismic anisotropy in subduction zones with “cold” geotherms.</p>

Seismic signals and noise recorded on a seafloor vertical hydrophone array and a colocated OBS during the LFASE experiment

1999 ◽  
Vol 89 (2) ◽  
pp. 423-432
Author(s):  
Mark A. Riedesel ◽  
John A. Orcutt ◽  
J. Andrew Adams
Keyword(s):  
Signal To Noise Ratio ◽  
Seismic Refraction ◽  
Vertical Direction ◽  
Body Waves ◽  
P Wave ◽  
Ocean Bottom ◽  
Seismic Signals ◽  
Signal To Noise ◽  
P Waves ◽  
Hydrophone Array

Abstract During the LFASE experiment conducted in the summer of 1989, a vertical hydrophone array (VHA) was deployed at the site of DSDP Hole 534b in the Blake-Bahama basin. The VHA consisted of 16 hydrophones spaced 30 m apart rising vertically above the seafloor with the bottom of the array anchored to a digital recording package located on the seafloor. The main purpose of the experiment was to record seismic signals and noise above, on, and beneath the seafloor at the same site; the VHA recorded data from above the seafloor, an ocean-bottom seismograph (OBS) was on the seafloor, and a borehole geophone array (BHA) was in a borehole beneath the seafloor. Signal-to-noise measurements of P waves made for airgun shots recorded simultaneously by the VHA and the OBS show that the OBS has a signal-to-noise ratio (SNR) 5 to 10 dB greater than that of a single hydrophone, but P-wave stacks of the VHA data have SNRs for body waves 10 to 15 dB greater than the OBS. This implies that a vertical array is capable of recording distinct body waves out to significantly greater ranges than is a single receiver and so might be a useful tool in seismic refraction experiments. Other stacks of the VHA data revealed S waves not visible on the data from a single hydrophone or geophone. The VHA was also used to construct images of the wavefronts arriving from a distant source and to determine the vertical direction of the source. This array processing capability shows the potential of VHAs for use in 3D seismic reflection surveys, particularly in cases where it is not convenient to use a multi-channel streamer.

Simulation of thermoelastic waves based on the Lord-Shulman theory

Geophysics ◽  
2021 ◽  
Vol 86 (3) ◽  
pp. T155-T164
Author(s):  
Wanting Hou ◽  
Li-Yun Fu ◽  
José M. Carcione ◽  
Zhiwei Wang ◽  
Jia Wei
Keyword(s):  
Thermal Conductivity ◽  
Expansion Coefficient ◽  
Velocity Dispersion ◽  
Wave Attenuation ◽  
Splitting Method ◽  
Grazing Incidence ◽  
P Wave ◽  
Staggered Grid ◽  
Thermal Mode ◽  
P Waves

Thermoelasticity is important in seismic propagation due to the effects related to wave attenuation and velocity dispersion. We have applied a novel finite-difference (FD) solver of the Lord-Shulman thermoelasticity equations to compute synthetic seismograms that include the effects of the thermal properties (expansion coefficient, thermal conductivity, and specific heat) compared with the classic forward-modeling codes. We use a time splitting method because the presence of a slow quasistatic mode (the thermal mode) makes the differential equations stiff and unstable for explicit time-stepping methods. The spatial derivatives are computed with a rotated staggered-grid FD method, and an unsplit convolutional perfectly matched layer is used to absorb the waves at the boundaries, with an optimal performance at the grazing incidence. The stability condition of the modeling algorithm is examined. The numerical experiments illustrate the effects of the thermoelasticity properties on the attenuation of the fast P-wave (or E-wave) and the slow thermal P-wave (or T-wave). These propagation modes have characteristics similar to the fast and slow P-waves of poroelasticity, respectively. The thermal expansion coefficient has a significant effect on the velocity dispersion and attenuation of the elastic waves, and the thermal conductivity affects the relaxation time of the thermal diffusion process, with the T mode becoming wave-like at high thermal conductivities and high frequencies.

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.

Evidence of tectonic release from underground nuclear explosions in long-period P waves

1983 ◽  
Vol 73 (2) ◽  
pp. 593-613
Author(s):  
Terry C. Wallace ◽  
Donald V. Helmberger ◽  
Gladys R. Engen
Keyword(s):  
Upper Mantle ◽  
High Stress ◽  
Test Site ◽  
Strike Slip ◽  
Body Waves ◽  
Release Time ◽  
P Waves ◽  
Underground Nuclear Explosions ◽  
Nuclear Explosions ◽  
Long Period

abstract In this paper, we study the long-period body waves at regional and upper mantle distances from large underground nuclear explosions at Pahute Mesa, Nevada Test Site. A comparison of the seismic records from neighboring explosions shows that the more recent events have much simpler waveforms than those of the earlier events. In fact, many of the early events produced waveforms which are very similar to those produced by shallow, moderate-size, strike-slip earthquakes; the phase sP is particularly obvious. The waveforms of these explosions can be modeled by assuming that the explosion is accompanied by tectonic release represented by a double couple. A clear example of this phenomenon is provided by a comparison of GREELEY (1966) and KASSERI (1975). These events are of similar yields and were detonated within 2 km of each other. The GREELEY records can be matched by simply adding synthetic waveforms appropriate for a shallow strike-slip earthquake to the KASSERI observations. The tectonic release for GREELEY has a moment of 5 ՠ1024 dyne-cm and is striking approximately 340°. The identification of the sP phase at upper mantle distances indicates that the source depth is 4 km or less. The tectonic release time function has a short duration (less than 1 sec). A comparison of these results with well-studied strike-slip earthquakes on the west coast and eastern Nevada indicate that, if tectonic release is triggered fault motion, then the tectonic release is relatively high stress drop, on the order of several hundred bars. It is possible to reduce these stress drops by a factor of 2 if the tectonic release is a driven fault; i.e., rupturing with the P velocity. The region in which the stress is released for a megaton event has a radius of about 4 km. Pahute Mesa events which are detonated within this radius of a previous explosion have a substantially reduced tectonic release.

Source dynamics of the 1979 Imperial Valley earthquake from near-source observations (of ground acceleration and velocity)

1982 ◽  
Vol 72 (6A) ◽  
pp. 1957-1968
Author(s):  
Mansour Niazi
Keyword(s):  
Particle Acceleration ◽  
Particle Velocity ◽  
High Frequency ◽  
Observation Point ◽  
P Wave ◽  
Imperial Valley ◽  
Ground Acceleration ◽  
Data Set ◽  
Ground Shaking ◽  
Fault Surface

abstract Two sets of observations obtained during the 15 October 1979 Imperial Valley earthquake, MS 6.9, are presented. The data suggest different dynamic characteristics of the source when viewed in different frequency bands. The first data set consists of the observed residuals of the horizontal peak ground accelerations and particle velocity from predicted values within 50 km of the fault surface. The residuals are calculated from a nonlinear regression analysis of the data (Campbell, 1981) to the following empirical relationships, PGA = A 1 ( R + C 1 ) − d 1 , PGV = A 2 ( R + C 2 ) − d 2 in which R is the closest distance to the plane of rupture. The so-calculated residuals are correlated with a positive scalar factor signifying the focusing potential at each observation point. The focusing potential is determined on the basis of the geometrical relation of the station relative to the rupture front on the fault plane. The second data set consists of the acceleration directions derived from the windowed-time histories of the horizontal ground acceleration across the El Centro Differential Array (ECDA). The horizontal peak velocity residuals and the low-pass particle acceleration directions across ECDA require the fault rupture to propagate northwestward. The horizontal peak ground acceleration residuals and the high-frequency particle acceleration directions, however, are either inconclusive or suggest an opposite direction for rupture propagation. The inconsistency can best be explained to have resulted from the incoherence of the high-frequency radiation which contributes most effectively to the registration of PGA. A test for the sensitivity of the correlation procedure to the souce location is conducted by ascribing the observed strong ground shaking to a single asperity located 12 km northwest of the hypocenter. The resulting inconsistency between the peak acceleration and velocity observations in relation to the focusing potential is accentuated. The particle velocity of Delta Station, Mexico, in either case appears abnormally high and disagrees with other observations near the southeastern end of the fault trace. From the observation of a nearly continuous counterclockwise rotation of the plane of P-wave particle motion at ECDA, the average rupture velocity during the first several seconds of source activation is estimated to be 2.0 to 3.0 km/sec. A 3 km upper bound estimate of barrier dimensions is tentatively made on the basis of the observed quasiperiodic variation of the polarization angles.

Estimates of velocity structure and source depth using multiple P waves from aftershocks of the 1987 Elmore Ranch and Superstition Hills, California, earthquakes

1991 ◽  
Vol 81 (2) ◽  
pp. 508-523
Author(s):  
Jim Mori
Keyword(s):  
Velocity Structure ◽  
Velocity Model ◽  
P Wave ◽  
Sediment Layer ◽  
Source Depth ◽  
Imperial Valley ◽  
P Waves ◽  
Main Shocks ◽  
Aftershock Sequences ◽  
Event Record

Abstract Event record sections, which are constructed by plotting seismograms from many closely spaced earthquakes recorded on a few stations, show multiple free-surface reflections (PP, PPP, PPPP) of the P wave in the Imperial Valley, California. The relative timing of these arrivals is used to estimate the strength of the P-wave velocity gradient within the upper 5 km of the sediment layer. Consistent with previous studies, a velocity model with a value of 1.8 km/sec at the surface increasing linearly to 5.8 km/sec at a depth of 5.5 km fits the data well. The relative amplitudes of the P and PP arrivals are used to estimate the source depth for the aftershock distributions of the Elmore Ranch and Superstition Hills main shocks. Although the depth determination has large uncertainties, both the Elmore Ranch and Superstition Hills aftershock sequences appear to have similar depth distribution in the range of 4 to 10 km.

Reflection moveout approximations for P-waves in a moderately anisotropic homogeneous tilted transverse isotropy layer

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
3D Space ◽  
Axis Of Symmetry ◽  
Relative Errors ◽  
Symmetry Tests ◽  
Ti Media

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.

Sign in / Sign up

Export Citation Format

Share Document