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 109 records · Page 6

Distinguishing Unique Earthquakes with Overlapping Signals in Oklahoma

Abstract During routine operations monitoring Oklahoma earthquakes, we found that certain earthquakes occurred closely both in space and time and had overlapping phases at the recording stations. Through further scrutiny and analysis, we determined that rather than being distinctly different earthquakes, some of the earthquakes exhibited multiphase arrivals and longer than expected coda due to unique ray paths that encounter impedance contrasts such as at the sedimentary rock-basement. Of course, some of these events truly were distinct events, which we term overlapping earthquakes, for which perceived coda duration overlaps and obscures the phase arrivals of the second event due to the source proximity in both time and space. We detail our classification scheme to separate the local earthquakes in Oklahoma as single, overlapping earthquakes, or those associated with multiphase arrivals. We forward model seismic wave propagation in a 2D crustal model and develop a methodology that utilizes waveform correlation to distinguish phases from overlapping earthquakes to those from crustal reverberations. Duration analysis shows a more elongated duration, qualitatively similar to the duration produced by overlapping earthquakes, at the sites where multiphase arrivals are observed.

Geochemistry & Geophysics↗

IC reports (2021)

We requested HPC support to continue research on the seismic waves generated by impacts. We used the following codes: (1) the Hybrid Optimization Software Suite (HOSS), developed at LANL. HOSS is based on a combined Finite and Discrete Element Method (FDEM). New material models are developed for the sedimentary rocks. HOSS has been recently benchmarked to iSale and FLAG codes. SPECFEM3D is an open-source code developed since the last 90s. It won the Gordon Bell award for best performance in 2003, was finalist again in 2008 for a run at 0.16 petaflops on 149,784 cores on the ‘Jaguar’ Cray system at Oak Ridge National Laboratory. It also won the BULL Joseph Fourier supercomputing award in 2010.; SW4 is a 4th-order finite difference code developed at LLNL which is currently actively developed to handle complex 3D models and to be ported on future exascale platforms. We assessed our need to a total of 3.9M CPU-hrs for year 1 and and 3.1 M for year 2.

79 ASTRONOMY AND ASTROPHYSICS↗

Challenges and Potential of Waveform Modeling for Crustal Scale Predictions

Waveform modeling is crucial to improving our understanding of observed seismograms. Forward simulation of wavefields provides quantitative methods of testing interactions between complicated source functions and the propagation medium. Here, we discuss three experiments designed to improve under standing of high frequency seismic wave propagation. First, we compare observed and predicted travel times of crustal phases for a set of real observed earthquakes with calculations and synthetic seismograms. Second, we estimate the frequency content of a series of nearly co-located earthquakes of varying magnitude for which we have a relatively well- known 1D velocity model. Third, we apply stochastic perturbations on top of a 3D tomographic model and qualitatively assess how those variations map to differences in the seismograms. While different in scope and aim, these three vignettes illustrate the current state of crustal scale waveform modeling and the potential for future studies to better constrain the structure of the crust.

58 GEOSCIENCES↗

Toward Global Regional Seismic Moment Tensor Inversion with Three-Dimensional Earth Models for Nuclear Explosion Monitoring with Sparse Networks: Demonstration of Reciprocity for Strain Greens Tensor Database Simulation with Salvus

Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).

58 GEOSCIENCES↗

S&TR September 2025: Computing Grand Challenge Turns 20

Livermore’s Computing Grand Challenge Program enters its 20th year with more unclassified high-performance computing (HPC) power than ever before. This unique, peer-reviewed competition awards HPC allocations on top supercomputers to multidisciplinary teams with high-impact projects. The Grand Challenge encourages researchers to innovate, pushes scientific discovery to new heights, improves the Laboratory’s HPC capabilities, and extends HPC accessibility to collaborators. Awardees must adapt to successive generations of HPC hardware and learn to run simulations at scale. The feature article spotlights three Grand Challenge teams whose research broke new ground in key scientific pursuits—the essence of dark matter, explosion-generated seismic waves, and protein interactions linked to cancer—while underscoring the importance of academic partnerships and considering the program’s future.

07 ISOTOPE AND RADIATION SOURCES↗

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗

Iron Bonding with Light Elements: Implications for Planetary Cores Beyond the Binary System

Light element alloying in iron is required to explain density deficit and seismic wave velocities in Earth’s core. However, the light element composition of the Earth’s core seems hard to constrain as nearly all light element alloying would reduce the density and sound velocity (elastic moduli). The alloying light elements include oxidizing elements like oxygen and sulfur and reducing elements like hydrogen and carbon, yet their chemical effects in the alloy system are less discussed. Moreover, Fe-X-ray Absorption Near Edge Structure (Fe-XANES) fingerprints have been studied for silicate materials with ferrous and ferric ions, while not many X-ray absorption spectroscopy (XAS) studies have focused on iron alloys, especially at high pressures. To investigate the bonding nature of iron alloys in planetary interiors, we presented X-ray absorption spectroscopy of iron–nitrogen and iron–carbon alloys at high pressures up to 50 GPa. Together with existing literature on iron–carbon, –hydrogen alloys, we analyzed their edge positions and found no significant difference in the degree of oxidation among these alloys. Pressure effects on edge positions were also found negligible. Our theoretical simulation of the valence state of iron, alloyed with S, C, O, N, and P also showed nearly unchanged behavior under pressures up to 300 GPa. This finding indicates that the high pressure bonding of iron alloyed with light elements closely resembles bonding at the ambient conditions. We suggest that the chemical properties of light elements constrain which ones can coexist within iron alloys.

Advanced Photon Source↗

Wellbore integrity assessment with casing-based advanced sensing

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage 1 and waste disposal. After installation, well casing and cement are subject to mechanical stress due to near-well pressure changes and fluid induced corrosion. This is exacerbated for geothermal wells where produced fluid is at high temperature and corrosive. The current state-of-the-art technologies for wellbore integrity assessments are an array of cased hole logging tools. Wireline deployed acoustic, electromagnetic and mechanical tools are all available to inspect steel casing corrosion and casing-cement bond and these tools can provide high-resolution assessment of borehole conditions. They are intrusive, however, in terms of borehole preparation and interruption to the normal operation of the wells, and not suitable for high temperature or highly deviated well deployments. In addition, these measurements are performed infrequently due to high cost, and are therefore incapable of providing frequent data to better predict borehole degradation trajectory, which can help provide early warning of potential borehole failures. For this project we are developing a suite of novel, non-invasive, casing based tools for wellbore integrity assessment, combining fast/low cost screening with higher-precision investigation. Our approach is based on monitoring the response of the casing when energized at the wellhead, thereby interrogating the casing without well intervention. Lab, field and numerical approaches are used in our study. During the early stage of the research, we focus on numerical simulations, which have shown the sensitivity of the low frequency electromagnetic (EM) signals to changes in borehole depths and have successfully tested the concept at a field site with different length well casings. Initial seismic modeling efforts have also demonstrated our capability to simulate seismic tube wave and seismic field alterations due to borehole breakage and associated fluid leakage. Further numerical, laboratory and field experiments are underway for additional technology sensitivity analysis, particularly the transient EM/Seismic reflectometry methods, data acquisition optimization, and numerical simulation improvements.

Wilt, Michael↗

Predicting Seafloor Seismic Noise from Sea Surface Wave Heights and Vice Versa

Seismic hazard monitoring and global seismic tomography efforts are improved by recording seismic signals close to sources and at a variety of azimuths surrounding the source. While seismic networks have been deployed at many locations on land, seafloor stations are still sparse. Inclusion of ocean-based data would greatly improve the global coverage of seismic networks. However, the use of seafloor seismic data to complement land-based detection and characterization of events is complicated by the fact that ambient noise on the seafloor is in general much higher than noise levels on land. This noise is driven primarily by sea surface waves and tides. In this project, I define an empirical relationship between wave height and seafloor seismic displacement, and tie this relationship to the microseism noise theory of Longuet-Higgins and Ursell (1948). I find that while there appears to be a linear relationship between wave height and seismic displacement, wave height values can only be used with low confidence to predict seismic displacement and vice versa, and furthermore, factors other than wave height are significant in determining the level of seismic noise at a particular site.

58 GEOSCIENCES↗

Seismic monitoring of well integrity

Seismic tube waves, produced by flexure of the well boundary, pressure changes in the fluid in the well, and deformation of the material immediately surrounding the well, are particularly sensitive to variations in the state of the well. In this work, we evaluate a direct approach for generating and observing tube waves as a means of detecting well damage. While we find that it can be difficult to reliably excite observable tube waves without a very strong surface source, time-frequency techniques can be employed to increase the detectability of tube wave reflections. New technologies, particularly distributed acoustic sensing, hold great promise for evaluating well integrity by monitoring tube waves, temperature changes, and seismic noise due to well deformation and fluid leakage.

58 GEOSCIENCES↗

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗

The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling

Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.

58 GEOSCIENCES↗

Improved Parametric Models for Explosion Pressure Signals Derived From Large Datasets

Accurate recording and characterization of explosion-induced pressure signals are key components of the forensic analysis of explosion events in the atmosphere. Parametric overpressure models based on several key waveform features (peak overpressure, positive pulse duration, and impulse) are widely used to estimate explosion energy in terms of trinitrotoluene equivalent yield. However, those models are often developed by a limited dataset, including only a few events or recordings at relatively short propagation distances. Here, we develop empirical waveform-parameter models based on a regression analysis of a large set of data curated from four chemical explosion experiments including 16 detonations. We measured peak overpressure and impulse for positive and negative phases from 1000 pressure signals recorded at local ranges (<20 km) with scaled distance up to 8000 m/kg 1/3 . Additionally, the measured waveform parameters showed large variation with respect to observing distances indicating the effects of atmospheric propagation. In this study, a second-order polynomial model was used in a least-squares regression to account for those propagation effects and to improve data fitting. In addition to model parameters for waveform features, we also determined range-dependent model uncertainties based on data variance. The model uncertainty represents the prediction error of our models and can be critical to evaluating the uncertainty of yield estimate.

58 GEOSCIENCES↗

Monitoring Seismic Velocity Changes Across the San Jacinto Fault Using Train‐Generated Seismic Tremors

Abstract Microseismic noise has been used for seismic velocity monitoring. However, such signals are dominated by low‐frequency surface waves that are not ideal for detecting changes associated with small tectonic processes. Here we show that it is possible to extract stable, high‐frequency body waves using seismic tremors generated by freight trains. Such body waves allow us to focus on small velocity perturbations in the crust with high spatial resolution. We report on 10 years of seismic velocity temporal changes at the San Jacinto Fault. We observe and map a two‐month‐long episode of velocity changes with complex spatial distribution and interpret the velocity perturbation as produced by a previously undocumented slow‐slip event. We verify the hypothesis through numerical simulations and locate this event along a fault segment believed to be locked. Such a slow‐slip event stresses its surroundings and may trigger a major earthquake on a fault section approaching failure.

58 GEOSCIENCES↗

Seismic constraints from a Mars impact experiment using InSight and Perseverance

NASA’s InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has operated a sophisticated suite of seismology and geophysics instruments on the surface of Mars since its arrival in 2018. On 18 February 2021, we attempted to detect the seismic and acoustic waves produced by the entry, descent and landing of the Perseverance rover using the sensors onboard the InSight lander. Similar observations have been made on Earth using data from both crewed and uncrewed spacecraft, and on the Moon during the Apollo era, but never before on Mars or another planet. This was the only seismic event to occur on Mars since InSight began operations that had an a priori known and independently constrained timing and location. It therefore had the potential to be used as a calibration for other marsquakes recorded by InSight. Here we report that no signal from Perseverance’s entry, descent and landing is identifiable in the InSight data. Nonetheless, measurements made during the landing window enable us to place constraints on the distance–amplitude relationships used to predict the amplitude of seismic waves produced by planetary impacts and place in situ constraints on Martian impact seismic efficiency (the fraction of the impactor kinetic energy converted into seismic energy).

79 ASTRONOMY AND ASTROPHYSICS↗

A “Floatilla” of Airborne Seismometers for Venus

Abstract Barometers floating on high‐altitude balloons in the relatively clement cloud layer on Venus could detect and characterize acoustic waves generated by seismic activity, avoiding the need for high‐temperature electronics required for surface seismology. Garcia et al. (2022, https://doi.org/10.1029/2022GL098844 ) recently demonstrated the detection of low‐frequency sound (infrasound) caused by earthquakes of magnitudes 7.3 and 7.5 from stratospheric balloons nearly 3,000 km away from the epicenter. They provided a preliminary demonstration of earthquake magnitude and location inversion, and the determination of S‐ and Rayleigh wave velocities using only their acoustic signature. Large earthquakes produce low‐frequency seismic waves that penetrate the interiors of planets; their detection at continental‐scale distances from a high‐vantage point demonstrates the feasibility of balloon‐based investigations of Venus' interior. We contextualize these results within the effort to perform seismology on Venus from balloons, discuss its limitations, and share perspectives on open research questions in this area.

Krishnamoorthy, Siddharth↗