SEARCH · Engineering Papers

Results for “infrasonic methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

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 OSTI2025

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗

DOE OSTI2024

Finite-frequency modeling of regional tropospheric infrasound using realistic atmospheres and terrain

Infrasonic waves have been observed to propagate to regional (greater than 15 km) distances through the troposphere. Infrasound propagation in the geometric acoustics approximation has shown that realistic terrain can scatter acoustic energy from tropospheric ducts; however, ray methods cannot intrinsically capture finite-frequency behavior such as diffraction. A two-dimensional finite-difference time-domain (FDTD) method has been developed to solve linearized equations for infrasound propagation with realistic terrain. Acoustic wave propagation over 100 km with both flat terrain and a Gaussian hill was first simulated in order to compare finite-frequency propagation with ray predictions. The effects of realistic terrain and atmospheres on infrasound signals from a 2012 surface explosion at the Utah Testing and Training Range are then investigated. Propagation through the troposphere is suggested by array processing results, but eigenrays are not predicted due to weak to nonexistent ducting conditions. FDTD modeling suggests that the inclusion of terrain and finite frequency effects helps explain much of the observed signal in a realistic scenario. Furthermore, these results suggest that geometric acoustics may underestimate propagation through the troposphere, and that recorded waveforms at regional distances may be noticeably affected by terrain.

58 GEOSCIENCES↗

DOE OSTI2024

The Utility of Infrasound in Global Monitoring of Extraterrestrial Impacts: A Case Study of the 2008 July 23 Tajikistan Bolide

Abstract Among the various observational techniques used for the detection of large bolides on a global scale is a low-frequency sound known as infrasound. Infrasound, which is also one of the four sensing modalities used by the International Monitoring System, offers continuous global monitoring and can be leveraged for planetary defense. Infrasonic records can provide an additional dimension for event characterization and a distinct perspective that might not be available through any other observational method. This paper describes the infrasonic detection and characterization of the bolide that disintegrated over Tajikistan on 2008 July 23. This event was detected by two infrasound stations at distances of 1530 and 2130 km. Propagation paths to one of the stations were not predicted by the model, despite being clearly detected. The presence of the signal is attributed to the acoustic energy being trapped in a weak but leaky stratospheric AtmoSOFAR channel. The infrasound signal analysis indicates that the shock originated at the point of the main breakup at an altitude of 35 km. The primary mode of the shock production of the signal detected at the two stations was a spherical blast resulting from the main gross fragmentation episode. The energy estimate, based on the signal period, is 0.17–0.51 kt of TNT equivalent, suggesting a mass of 6.6–23.5 tons. The corresponding object radius, assuming the chondritic origin, was 0.78–1.18 m.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗

DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗

DOE OSTI2024

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

42 ENGINEERING↗

DOE OSTI2024

Modeling Partial Reflection Paths for Infrasound Analysis

Numerical methods enabling simulation of scattered and partially reflected infrasonic propagation paths produced by interaction with fine-scale structure in the middle atmosphere have been implemented in the infraGA ray tracing software. This capability enables simulation of ensonification in the classical stratospheric “shadow zone” that has been observed during the Humming Roadrunner and LSECE surface explosion campaigns as well as in other data sets. In the case of LSECE, a pair of stations roughly 140 kilometers east of the source location observed arrivals with celerities (horizontal group velocities) slightly slower than observed stratospheric paths at similar azimuths. The arrivals exhibited increasing trace velocity later in the wavetrain indicating a steepening of the arrival path for longer or slower propagation paths. Simulation of partially reflected paths using the updated infraGA software methods finds good agreement between observed and predicted infrasonic ensonification at these locations within the stratospheric shadow zone. Further development of the partial reflection physics and comparison with other data sets is needed to more robustly understand how such anomalous infrasonic signals can be predicted; however, the demonstration of this capability is a promising first step in such analyses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Seismoacoustic Signatures Observed During a Long-Term Deployment of Infrasound Sensors at the Nevada National Security Site

Earthquakes have repeatedly been shown to produce inaudible acoustic signals (<20 Hz), otherwise known as infrasound. These signals can propagate hundreds to thousands of kilometers and still be detected by ground-based infrasound arrays depending on the source strength, distance between source and receiver, and atmospheric conditions. Another type of signal arrival at infrasound arrays is the seismic induced motion of the sensor itself, or ground-motion-induced sensor noise. Measured acoustic and seismic waves produced by earthquakes can provide insight into properties of the earthquake such as magnitude, depth, and focal mechanism, as well as information about the local lithology and atmospheric conditions. Large earthquakes that produce strong acoustic signals detected at distances greater than 100 km are the most commonly studied; however, more recent studies have found that smaller magnitude earthquakes ($M$ w <2.0) can be detected at short ranges. In that vein, this study will investigate the ability for a long-term deployment of infrasound sensors (deployed as part of the Source Physics Experiments [SPE] from 2014 to 2020) to detect both seismic and infrasonic signals from earthquakes at local ranges (<50 km). Methods used include a combination of spectral analysis and automated array processing, supported by U.S. Geological Survey earthquake bulletins. Further, this investigation revealed no clear acoustic detections for short range earthquakes. However, secondary infrasound from an $M$ w 7.1 earthquake over 200 km away was detected. Important insights were also made regarding the performance of the SPE networks including detections of other acoustic sources such as bolides and rocket launches. Finally, evaluation of the infrasound arrays is performed to provide insight into optimal deployments for targeting earthquake infrasound.

58 GEOSCIENCES↗

DOE OSTI2025

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗

DOE OSTI2022

stochprop Documentation (Release 1.0)

Simulations of infrasonic propagation in the atmosphere typically utilize a single atmospheric specification describing the acoustic sound speed, ambient winds, and density as a function of altitude. Due to the dynamic and sparsely sampled nature of the atmosphere, there is a notable amount of uncertainty in the atmospheric state at a given location and time so that a more robust analysis of infrasonic propagation requires inclusion of this uncertainty. This Python library, stochprop, has been implemented using methods developed jointly by infrasound scientists at Los Alamos National Laboratory (LANL) and the University of Mississippi’s National Center for Physical Acoustics (NCPA). This software library includes methods to quantify variability in the atmospheric state, identify typical seasonal variability in the atmospheric state and generate suites of representative atmospheric states during a given season, as well as perform uncertainty analysis on a specified atmospheric state given some level of uncertainty. These methods have been designed to interface between propagation modeling capabilities such as InfraGA/GeoAc and NCPAprop and signal analysis methods in the LANL InfraPy tool.

58 GEOSCIENCES↗

DOE OSTI2021

Synthetic Evaluation of Infrasonic Multipole Waveform Inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies are able to help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

58 GEOSCIENCES↗

DOE OSTI2020

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗

DOE OSTI2024

Reverse Time Migration of Explosion-generated Infrasound and Ground-Coupled Airwaves

Accurately locating sources of seismic and infrasonic energy is integral to global monitoring of earthquakes and explosions. Infrasound arrivals times can be used to calculate the origins of events that generate acoustic energy. Picking times of emergent infrasound arrivals, however, can be difficult and prone to uncertainty. Reverse time migration (RTM) is a waveform based location method that does not rely on picked arrival times. Here we use RTM to locate a known chemical explosion that generated acoustic and acoustic-to-seismic signals on 26 and 108 receivers, respectively. All location predictions are less than 24 km from the known location with time errors of less than three minutes. We find strong overall agreement between our results and those of existing RTM and arrival time based methods. Our initial results suggest that RTM is a promising method of event location using acoustic arrivals recorded on both infrasound and seismic instrumentation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

Back-Azimuth Estimation of Air-to-Ground Coupled Infrasound from Transverse Coherence Minimization [Slides]

Infrasonic signals are generated by a variety of natural and anthropogenic sources and can propagate long distances through a windy, layered atmosphere. When these signals impinge on the ground, the resulting deformation can be recorded by buried seismometers, we call this signal a ground coupled airwave (GCA). GCAs are relatively common, and several methods now exist to leverage an infrasound microphone to separate GCAs from pure seismic motion.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Evaluation of a generalized least squares algorithm for infrasound beamforming with coherent background noise

Infrasonic signals of interest can occur during periods with persistent, coherent, background noise, which may be natural or anthropogenic. For high signal-to-noise (SNR) ratio transient signals, an ‘overprinting’ of the coherent background may occur, and the signal may still be detected. However, this approach fails for low SNR signals of interest, which may be obscured by coherent noise. An infrasound beamforming method based on generalized least squares (GLS) is investigated for detecting transient signals of interest in the presence of coherent and incoherent background noise. This approach relies on an estimate of the noise covariance, captured in a covariance matrix, to effectively null contributions to the array response from noisy directions of arrival. Synthetic array data is used to investigate the performance of the GLS beamformer compared to the Bartlett beamformer when coherent and incoherent backgrounds are present. Additionally, the effects of array element number and relative strength of the interfering signal on the GLS estimates is investigated. GLS empirical area under the curve estimates suggest that the beamformer can recover coherent power for a signal of interest lower in amplitude than the coherent background, but this effectiveness degrades more quickly with SNR for a four element array compared to a six or eight element infrasound array. Finally, infrasound from the Forensic Surface Experiment, a bolide signal observed at IMS array I37NO, and a volcanic signal recorded at the Alaska Volcano Observatory array ADKI are used to evaluate GLS performance on recorded data. A ten minute window was used to capture the background noise, and the coherent background signal was nulled in all three examples.

58 GEOSCIENCES↗

DOE OSTI2023

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗