Engineering Papers⌕ Search

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.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multi-Event Explosive Seismic Source for the 2022 M w 6.3 Hunga Tonga Submarine Volcanic Eruption

The eruption of the Hunga Tonga–Hunga Ha’apai submarine volcano on 15 January 2022 produced a variety of geophysical responses, including a significant seismic signal. We study the seismic source process of this event by inverting for moment tensors (MTs) using regional surface waves (Rayleigh, Love). By comparing inversion results for the eruption with eight nearby earthquakes, we show that it is possible to discriminate MT source types. Our inversion yields a shallow explosive source for the eruption and reveals the importance of trade-offs among depth, magnitude, and source type. We illustrate these trade-offs by representing the misfit variations over the eigenvalue lune. Finally, we invert for the source-time function of the sequence of explosions that occurred in the first minutes of the eruption. The multi-event source-time function comprises four subevents spanning ~270 s, with a total magnitude estimate of M w 6.34 ± 0.10.

58 GEOSCIENCES↗
DOE OSTI2024

Seismoacoustic measurements of the OSIRIS-REx re-entry with an off-grid Raspberry PiShake

Hypersonic re-entries of spacecraft are valuable analogues for the identification and tracking of natural meteoroids re-entering the Earth's atmosphere. We report on the detection of seismic and acoustic signals from the OSIRIS-REx landing sequence, acquired near the point of peak capsule heating and recorded using a fully off-grid Raspberry PiShake sensor. This simple setup is able to record all the salient features of both the seismic and acoustic wavefields; including the primary shockwave, later reverberations, and possible locally induced surface waves. Peak overpressures of 0.7 Pa and ground velocities of 2x10 -6 m/s yield lower bound on the air-to-ground coupling factor between 3 and 44 Hz of 1.4x10 -6 m/s/Pa, comparable to results from other re-entries.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

Geotechnical investigation report for proposed array of six 40-meter diameter antennas, Pioneer site, DSS 11, Goldstone, California tracking complex

The geotechnical investigation was conducted in three disciplines: (1) geological field reconnaissance of the general area of proposed construction; (2) geophysical seismic refraction survey of the localized area surrounding the six proposed antenna sites, including shear wave velocity determination; and (3) detailed foundation engineering investigation of each of the six sites. The investigations indicate that the six sites selected are relatively free from geologic hazards which would inhibit the proposed construction or future antenna operations.

Sweitzer, J. S.↗
NASA NTRS2004

Near Space Tracking of the EM Phenomena Associated with the Main Earthquakes

Searching for electromagnetic (EM) phenomena originating in the Earth's crust prior to major earthquakes (M>5) are the object of this exploratory study. We present the idea of a possible relationship between: (1) electro-chemical and thermodynamic processes in the Earth's crust and (2) ionic enhancement of the atmosphere/ionosphere with tectonic stress and earthquake activity. The major source of these signals are proposed to originate from electromagnetic phenomenon which are responsible for these observed pre-seismic processes, such as, enhanced IR emission, also born as thermal anomalies, generation of long wave radiation, light emission caused by ground-to-air electric discharges, Total Electron Content (TEC) ionospheric anomalies and ionospheric plasma variations. The source of these data will include: (i) ionospheric plasma perturbations data from the recently launched DEMETER mission and currently available TEC/GPS network data; (ii) geomagnetic data from ORSTED and CHAMP; (iii) Thermal infra-red (TIR) transients mapped by the polar orbiting (NOAA/AVHRR, MODIS) and (iv) geosynchronous weather satellites measurements of GOES, METEOSAT. This approach requires continues observations and data collecting, in addition to both ground and space based monitoring over selected regions in order to investigate the various techniques for recording possible anomalies. During the space campaign emphasis will be on IR emission, obtained from TIR (thermal infrared) satellites, that records land/sea surface temperature anomalies and changes in the plasma and total electron content (TEC) of the ionosphere that occur over areas of potential earthquake activity.

Ouzounov, Dimitar↗
NASA NTRS2008

Verification of Atmospheric Signals Associated with Major Seismicity by Space and Terrestrial Observations

Observations from the last twenty years suggest the existence of electromagnetic (EM) phenomena during or preceding some earthquakes [Haykawa et a!, 2004; Pulinets at al, 1999,2004, 2006, Ouzounov et all 2007 and Liu et all 20041. Both our previous studies [Pulinets at al, 2005, 2006, Ouzounov et al, 2006, 20071 and the latest review by the Earthquake Remote Precursor Sensing panel [ERPS; 2003- 20051; have shown that there were precursory atmospheric TIR signals observed on the ground and in space associated with several recent earthquakes. [Tramutoli at al, 2005, 2006, Cervone et al, 2006, Ouzounov et all 2004,2006JT.o study these signals, we applied both multi parameter statistical analysis and data mining methods that require systematic measurements from an Integrated Sensor Web of observations of several physical and environmental parameters. These include long wave earth infra-red radiation, ionospheric electrical and magnetic parameters, temperature and humidity of the boundary layer, seismicity and may be associated with major earthquakes. Our goal is to verify the earthquake atmospheric correlation in two cases: (i) backward analysis - 2000-2008 hindcast monitoring of multi atmospheric parameters over the Kamchatka region, Russia ; and (ii) forward real-time alert analysis over different seismo-tectonic regions for California, Turkey, Taiwan and Japan. Our latest results, from several post-earthquake independent analyses of more then 100 major earthquakes, show that joint satellite and some ground measurements, using an integrated web, could provide a capability for observing pre-earthquake atmospheric signals by combining the information from multiple sensors into a common framework. Using our methodology, we evaluated and compared the observed signals preceding the latest M7.9 Sichuan earthquake (0511212008), M8.0 earthquake in Peru (0811512007), M7.6 Kashmir earthquake (1010812005) and M9.0 Sumatra earthquake (1212812004). We found evidence of the systematic appearance of both atmospheric and ionospheric anomalies preceding most of the major events during the period of our analysis 2001 -2008.

Taylor, Patrick↗
NASA NTRS1994

Seismic evidence for volatiles at large depth in the Earth

High resolution tomographic images that have been obtained of the subduction zones in the west Pacific do not only show very low P wave velocities above the leading edge of the slab at shallow depth, but also below 300 km. The two zones are generally, but not always, separated by a zone of normal shear velocity near 200 km depth. These findings confirm some earlier inferences from local seismic data in Japan, and similar findings of low Vp and low Q zones beneath South America. Surprisingly, such deep seismic low velocity zones have recently also been discovered beneath the locations of ancient subduction zones. A study of upper mantle shear velocity in Central Europe shows a similar distribution of low velocity zones beneath the Tornquist-Teisseyre line, the former west coast of the old continent of Baltica, and the site of the subducting Tornquist ocean in the early Palezoic. Preliminary results from a Russian-French experiment shows low P velocities below 250 km under the Urals, while an older tomographic study shows such low P velocities beneath the northern Appalachians.

Nolet, Guust↗
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
NASA NTRS2003

Control of the Laser Interferometer Space Antenna

The detection of gravity waves will open a new window of observation on the universe. Unlike typical observatories, which detect electromagnetic waves traveling through space-time, the Laser Interferometer Space Antenna (LISA) will detect ripples in space-time itself. Science targets include galactic binaries, merging supermassive black holes, intermediate-mass/seed black holes, and cosmological backgrounds. Gravity waves are detected by measuring the strain in space, i.e. the change in distance between a set of masses (test masses or proof masses) separated by a great distance. Ground based detection of gravity waves by Laser Interferometer Gravitational Wave Observatory (LIGO) and other observatories are possible with laser interferometry; hut the relatively short arm length (4 km) and seismic noise limit the measurement band to above 10 Hz on Earth. LISA also uses laser interferometric measurement of the change in distance between test masses, but does it in space. Each LISA spacecraft embodies two test masses. Space allows very long arm lengths (5 million km for LISA) and a very quiet acceleration environment (3.5x10(exp -15) meters per second squared/Hertz (sup 0.5) for LISA), which allows for the detection of gravity wave strains to a best sensitivity of 3x10(exp -24) strain/Hertz (sup 0.5) over the measurement band of 10(exp -4) to 10(exp -1) Hertz for a one-year observation. The LISA mission consists of three spacecraft in heliocentric orbit. The orbits are chosen so that the three spacecraft form a roughly equilateral triangle with its center located at a radius of 1 AU and 20 degrees behind the Earth, as shown. Requirements are placed on the rotational and translational dynamics of each spacecraft to ensure that the proper sensitivity for science measurements can be achieved.

Maghami, Peiman↗
DOE OSTI2026

Lunar accelerometer network gravitational observatory (LANGO)

With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Underground hydrogen storage leakage detection and characterization based on machine learning of sparse seismic data

Underground hydrogen storage (UHS) is considered as a scalable approach for massive storage and seasonal extraction of hydrogen (H 2 ). Although conventional leakage detection and characterization methods based on time-lapse seismic imaging and inversion generally apply to H 2 leakage detection problem, a high-fidelity yet cost effective geophysics approach is still missing to reliably inform leakage location and properties based on very sparse data. In response, we develop a novel supervised machine learning method to detect and characterize H 2 leakage from UHS. The input to our neural network are sparse time-lapse seismic waveforms, while the output from the neural network includes the spatial location and physical properties of a H 2 leakage. Here, we generate high-quality time-lapse waveforms using the elastic-wave equations to train the neural network. We train and validate our machine learning model and find that it attains high accuracy in using extremely sparse time-lapse seismic data to detect and characterize H 2 leakage. Our investigation is the first systematic study that focuses on applying machine learning to subsurface H 2 leakage detection and characterization and could potentially serve as a cost-effective geophysical tool for underground hydrogen leakage detection and characterization with high fidelity.

08 HYDROGEN↗
DOE OSTI2022

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
NASA NTRS2019

The Seismic Signatures of Recently Formed Impact Craters on Mars

We investigated the seismic signatures of recent impact crater clusters on Mars that would be recorded by the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight)seismometers. We used a database of 77 measured and dated impact sites, with craters with diameters between 2.1 and 33.8 m, along with inferred impact angle, bolide trajectory, and varying target material properties to empirically scale for the momentum, expected seismic source function, and radiation pattern of impacts. The impact source is simulated in a local 3‐Dfinite difference wave propagation code and coupled to teleseismic distances by scaling the spectra of 1‐D global synthetic seismograms. We use the InSight seismometer noise floors to estimate detectability of impact(s) across azimuth and distance. Our experiments reveal that impact clusters have a higher peak corner frequency resulting from energy contributed by smaller craters to the power spectrum. We also find that the time separation between individual impacts in a cluster is small (<10–15 ms) and requires a seismometer closely situated to the source (<10 km) and a high sampling rate (>100 Hz) to resolve individual impacts within the cluster. Two of the clusters in our database(>20‐m effective diameter) would have been detectable by InSight, with the assumptions that the Martian background noise and seismic attenuation are both low. Joint detection of surface changes from newly formed crater(s) in images and by InSight will provide precise source locations that are crucial for constraining the internal structure of Mars

N. C. Schmerra↗
NASA NTRS1984

A Crustal Structure Study of South America

Variations in the average crustal thickness of the South American platform were determined by measuring fundamental mode Rayleigh wave dispersion between nine pairs of South American WWSSN seismograph stations and inverting the measurements to obtain average shear wave velocity-depth models across the platform. Additional models were derived from the Rayleigh wave dispersion data previously measured over eastern South America. Both sets of models were interpreted in terms of average crustal thickness and average crustal shear wave velocity. The results obtained were depicted in a contour map of crustal thickness. Average values for seismic velocities in the crust and upper mantle were calculated and compared with maps of regional gravity models, MAGSAT magnetic anomalies, and major tectonic features. The correlations and relationships discovered are discussed.

Renbarger, K. S.↗
DOE OSTI2023

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

MASTODON: An Open-Source Software for Seismic Analysis and Risk Assessment of Critical Infrastructure

Seismic analysis and risk assessment of safety-critical infrastructure like hospitals, nuclear power plants, dams, and facilities handling radioactive materials involve computationally intensive numerical models and coupled multiphysics scenarios. They are also performed in a strict regulatory environment that requires high software quality assurance standards, and in the case of safety-related nuclear facilities, a conformance to the American Society of Mechanical Engineers Nuclear Quality Assurance (NQA-1) standard. This paper introduces the open-source finite-element software, MASTODON (Multi-hazard Analysis of Stochastic Time-Domain Phenomena), which implements state-of-the-art seismic analysis and risk assessment tools in a quality-controlled environment. MASTODON is built on MOOSE (Multi-physics Object-Oriented Simulation Environment), which is a highly parallelizable, NQA-1 conforming, coupled multiphysics, finite-element framework developed at Idaho National Laboratory. MASTODON is capable of fault rupture and source-to-site wave propagation using the domain reduction method, nonlinear site response, and soil-structure interaction analysis, implicit and explicit time integration, automated stochastic simulations, and seismic probabilistic risk assessment. When coupled with other MOOSE applications, MASTODON can also solve strongly and weakly coupled multiphysics problems. This paper presents a summary of the capabilities of MASTODON and some demonstrative examples.

42 ENGINEERING↗
DOE OSTI2023

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗