Physical mechanisms of seismic-wave attenuation
Seismic wave attenuation mechanisms emphasizing partial melting, grain boundary relaxation and high temperature internal friction background
Anderson, D. L.↗
SEARCH · Engineering Papers
Results for “seismic wave”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Modeling of Seismic Waves Through Geologic Metamaterials
This project conducted a modeling study on seismic invisibility cloaks that render geologic targets invisible to seismic waves, using the concept of seismic metamaterials. We present a parametric numerical study on the behaviors of seismic waves through cloaks with different design parameters as well as degrees of geologic heterogeneity. In addition, a seismic cloaking strategy is proposed for a future field-scale experiment at a real-world test bed. This feasibility study will guide future field experiment designs and ultimately allow us to conduct systematic field-scale tests employing Sandia’s existing resources and field expertise. The ultimate goal is to develop methods and design parameters of seismic invisibility cloaks to protect against natural and man-made seismic waves. Seismic cloaking has potential applications in several areas of national security, energy, and natural hazard reduction.
58 GEOSCIENCES↗
The Effect of Redox Conditions on Seismic Waves in Iron-Bearing Olivine: Implications for Understanding Planetary Interiros Through Seismilogy
Seismic data, inclusive of velocities and attenuation, can be utilized to elucidate the physical state of planetary interiors]. However, numerous micromechanical factors have been either experimentally demonstrated, or theoretically considered, to affect the propagation and dissipation of seismic energy within crystalline solids - including, but not limited to, changes in grain size, temperature, melt fraction, pressure and dislocation density. Thus, observed variations in seismic wave speeds and attenuation may be used to ultimately ‘map’ variations in physical properties, such as those listed above, within planetary bodies. But, in order to complete a successful inversion of seismic data into representations of physical properties, a first requirement is to obtain a fundamental laboratory based understanding of how each of these possible factors individually influences seismic waves. Here we conduct an experimental study with the initial objective to further understand one of the most commonly invoked, yet least studied, mechanisms that could alter intrinsic seismic wave attenuation: water content (occurring as chemically-bound hydroxyl). The historical basis for determining the effect of water on seismic properties was established predominantly through analogy with large-strain creep experiments conducted on olivine under water-saturated conditions. While these deformation experiments routinely demonstrate a weakening of olivine in the presence of water, they represent a fundamentally different deformation regime in comparison to the microstrains experienced due to a passing seismic wave. Thus, in order to directly assess the effects of water on seismic properties, small-strain experiments are required. Substantially modified seismic properties in the presence of water have been observed previously at low strains and low frequencies, but only in a single exploratory study conducted under water-saturated conditions. Thus, to properly test the theoretical predictions we conducted a systematic study of the seismic properties of olivine using low-frequency torsional oscillation on aggregates containing varying concentrations of bound hydroxyl, for the first time at under saturated conditions.
Cline, C. J., II↗
The damping of seismic waves and its determination from reflection seismograms
The damping in theoretical waveforms is described phenomenologically and a classification is proposed. A method for studying the Earth's crust was developed which includes this damping as derived from reflection seismograms. Seismic wave propagation by absorption, attenuation of seismic waves by scattering, and dispersion relations are considered. Absorption of seismic waves within the Earth as well as reflection and transmission of elastic waves seen through boundary layer absorption are also discussed.
Engelhard, L.↗
Attenuation of seismic waves by grain boundary relaxation
Elastic seismic waves attenuation in polycrystalline ceramics and rocks by grain boundary relaxation, noting strong frequency dependence of earth mantle internal friction
Jackson, D. D.↗
A modular methodology for time-domain stochastic seismic wave propagation
Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.
58 GEOSCIENCES↗
Air-coupled seismic waves at long range from Apollo launchings.
Microphones and seismographs were co-located in arrays on Skidaway Island, Georgia, for the launchings of Apollo 13 and 14, 374 km to the south. Simultaneous acoustic and seismic waves were recorded for both events at times appropriate to the arrival of the acoustic waves from the source. The acoustic signal is relatively broadband compared to the nearly monochromatic seismic signal; the seismic signal is much more continuous than the more pulse-like acoustic signal; ground loading from the pressure variations of the acoustic waves is shown to be too small to account for the seismic waves; and the measured phase velocities of both acoustic and seismic waves across the local instrument arrays differ by less than 6 per cent and possibly 3 per cent if experimental error is included. It is concluded that the seismic waves are generated by resonant coupling to the acoustic waves along some 10 km of path on Skidaway Island.
Donn, W. L.↗
Investigations of acoustic-seismic effects at long range - Early-arriving seismic waves from Apollo 16
A reasonably comprehensive technical effort is described dealing with the investigations of acoustically generated seismic waves of Apollo 16 and Apollo 17 origin along the eastern seabord of the United States. This expanded effort is a continuation of earlier, rather successful detections of rocket-generated seismic disturbances on Skidaway Island, Georgia. The more recent effort has yielded few positive results other than a recording of an early-arriving seismic wave from Apollo 16 that was detected in Jacksonville. Evaluation of the negative results obtained in the Fort Monmouth area, with earlier studies of infrasound, local weather conditions, and geology, could be advantageous in the process of trying to gain a better insight into the acoustic-seismic resonance mechanism requiring phase-velocity matching at the atmosphere-ground interface.
Dalins, I.↗
Seismic-wave velocity patterns in some pahoehoe basalt flows
Seismic wave velocity measurements in pahoehoe basalt flows in lava beds, comparing with laboratory dilatational velocities
Greeley, R.↗
Seismic waves generated by sonic booms - A geoacoustical problem.
Air coupled seismic waves generated by jet fighters flying at high altitudes and supersonic speeds
Espinosa, A. F.↗
Analysis of Seismic Wave Amplification in Sedimentary Basins Using 3D Wavefield Simulations: Nenana Basin, Central Alaska
We investigate seismic wave amplification in Nenana basin, central Alaska, using 3D seismic wavefield simulations. We quantify frequency-dependent amplification by comparing synthetic seismograms among four different seismic velocity models: (a) a tomographic model, (b) a tomographic model with the uppermost 6.5 km replaced by a 1D basin profile, (c) a tomographic model with an embedded ellipsoidal basin, and (d) a tomographic model with an embedded realistic basin. For each model we perform wavefield simulations accurate up to 2 Hz for a set of 10 different earthquakes, which provide a range of frequencies and incident angles for waves interacting with the basin. Average amplification ratios are four on the horizontal components and seven on the vertical component. The dominant influence of the amplification is the 3D geometry of the basin, not the slow-velocity profile near the surface. Our synthetic amplification ratios are in general agreement with amplification estimates obtained from 14 stations that recorded the same set of earthquakes. Our approach offers a general strategy for documenting frequency-dependent basin amplification for a region with realistic basin structures and local earthquakes.
Geosciences↗
Analysis of long-period seismic waves excited by the May 18, 1980, eruption of Mount St. Helens - A terrestrial monopole
The eruption of Mount St. Helens on May 18, 1980, excited long-period seismic waves, and high-quality digital seismograms were recorded. The present investigation is concerned with the results of detailed analyses of Rayleigh and Love waves excited by this eruption. Since the elastic response of the earth is very accurately known, it is possible to retrieve the source parameters of this unique event from observations at far-field. It is shown that the source can be represented by a nearly horizontal single force. The conducted analysis is concerned with only long-period characteristics of the source. The short-period behavior of the source is difficult to determine from surface waves because the available knowledge of the earth's response is less accurate than at long periods.
Kanamori, H.↗
Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling
SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.
58 GEOSCIENCES↗
Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth
SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.
Geochemistry & Geophysics↗
Solving Seismic Wave Equations on Variable Velocity Models With Fourier Neural Operator
Here, in the study of subsurface seismic imaging, solving the acoustic wave equation is a pivotal component in existing models. The advancement of deep learning (DL) enables solving partial differential equations (PDEs), including wave equations, by applying neural networks to identify the mapping between the inputs and the solution. This approach can be faster than traditional numerical methods when numerous instances are to be solved. Previous works that concentrate on solving the wave equation by neural networks consider either a single velocity model or multiple simple velocity models, which is restricted in practice. Instead, inspired by the idea of operator learning, this work leverages the Fourier neural operator (FNO) to effectively learn the frequency domain seismic wavefields under the context of variable velocity models. We also propose a new framework paralleled FNO (PFNO) for efficiently training the FNO-based solver given multiple source locations and frequencies. Numerical experiments demonstrate the high accuracy of both FNO and PFNO with complicated velocity models in the OpenFWI datasets. Furthermore, the cross-dataset generalization test verifies that PFNO adapts to out-of-distribution velocity models. Finally, PFNO admits higher computational efficiency on large-scale testing datasets than the traditional finite-difference method. The aforementioned advantages endow the FNO-based solver with the potential to build powerful models for research on seismic waves.
58 GEOSCIENCES↗
Laboratory measurements of p-wave seismic Q on lunar and analog rocks
The longitudinal-wave, internal friction quality factor (Qp) of out-gassed rock subjected to hydrostatic confining pressure was measured by a technique of aerodynamically shaping the sample ends to remove most of the air drag. Q values of about 1000 and 100 were obtained at 100 MPa and in laboratory air, respectively. The temperature dependence of Q in lunar rock 70215.85 and an analog of lunar basalt was investigated over the range -100 to 450 deg C. A systematic increase in Q was observed at 50 Hz, 5 kHz and 50 kHz. When the temperature was lowered to -100 deg C moderate decreases in Q were observed from 100 to 250 deg C, and high Q values were obtained at 450 deg C. The temperature and pressure dependences of Q suggest that high lunar seismic Q values imply a very dry crust to depths of at least 50 km. These results further indicate that elastic waves are damped in volatile-rich rock by an absorption mechanism which involves changes in the bonding, structure and coverage of the molecular, physisorbed and chemisorbed H2O at crack and grain boundaries.
Tittmann, B. R.↗
SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle
SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.
58 GEOSCIENCES↗
Finite difference simulations of seismic scattering Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity
Synthetic seismographs that were obtained by the finite difference method are presently applied to the study of elastic and acoustic wave scattering in two-dimensional media with random spatial variations in seismic velocity. The seismograms are analyzed to determine the variation in travel times and waveforms across arrays of receivers. The random media with Gaussian and exponential correlation functions considered differ in the spectral falloff of their velocity fluctuations at wavelengths smaller than 2pi times the correlation distance. It is found that alternative models of crustal heterogeneity can be tested by improved measurements of the frequency dependence of the crustal Q at frequencies greater than about 1 Hz, assuming that scattering is responsible for most of the attenuation at such frequencies.
Frankel, A.↗