Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency estimation”

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 703 records · Page 39

Binary collisions of dark matter blobs

We describe the model-independent mechanism by which dark matter and dark matter structures heavier than ~ 8 × 10 11 GeV form binary pairs in the early Universe that spin down and merge both in the present and throughout the Universe’s history, producing potentially observable signals. Sufficiently dense dark objects will dominantly collide through binary mergers instead of random collisions. We detail how one would estimate the merger rate accounting for finite size effects, multibody interactions, and friction with the thermal bath. We predict how mergers of dark dense objects could be detected through gravitational and electromagnetic signals, noting that such mergers could be a unique source of high frequency gravitational waves. We rule out objects whose presence would contradict observations of the CMB and diffuse gamma-rays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamical friction from scalar dark matter in the relativistic regime

Light bosonic scalars (e.g., axions) may form clouds around black holes via superradiant instabilities or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we have studied numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. Here we show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar, the rapid force oscillations could, in principle, leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Strong parallel evidence of selection during switchgrass sward establishment in hybrid and lowland ecotypes

Switchgrass sward establishment results in up to 90% seedling mortality. The degree of selection during sward establishment has not been reported using modern genetic methods. Pooled leaf samples were sequenced from replicated swards of 46 half-sib families from two breeding groups (lowland and hybrid) before and through 3 years of stand establishment. Pooled allele frequencies were then assessed using fixation indices (Fst) and an independent data set was used to predict the polygenic impact of establishment selection on two traits (heading date and winter survivorship). Last, the DNA pools were assigned survival rankings to predict the sward survival genomically estimated breeding values within the training data set. Strong and parallel selection occured in both breeding groups. Five genomic regions exceeded the significant threshold of 99.9% in >10 families, indicating consistent selection across families and breeding groups. Polygenic trait predictions determined that establishment selection was partially associated with winter survivorship but resulted in variable heading date alterations. The genomewide variation is consistent with selection for a small number of related parental lines. This study observed strong selection for a small number of hybrid and coastal ecotype individuals which are promising germplasm sources for improved sward survival. This confirms prior reports of sward selection during grassland establishment and highlights the strength of pooled DNA sequencing for survival traits.

54 ENVIRONMENTAL SCIENCES↗

Wildfires increasingly impact western US fluvial networks

Wildfires are increasing globally in frequency, severity, and extent, but their impact on fluvial networks, and the resources they provide, remains unclear. We combine remote sensing of burn perimeter and severity, in-situ water quality monitoring, and longitudinal modeling to create the first large-scale, long-term estimates of stream+river length impacted by wildfire for the western US. We find that wildfires directly impact ~6% of the total stream+river length between 1984 and 2014, increasing at a rate of 342 km/year. When longitudinal propagation of water quality impacts is included, we estimate that wildfires affect ~11% of the total stream+river length. Our results indicate that wildfire activity is one of the largest drivers of aquatic impairment, though it is not routinely reported by regulatory agencies, as wildfire impacts on fluvial networks remain unconstrained. We identify key actions to address this knowledge gap and better understand the growing threat to fluvial networks, water security, and public health risks.

54 ENVIRONMENTAL SCIENCES↗

High-Spectral Resolution Lidar (HSRL) Instrument Handbook

High-spectral-resolution lidar (HSRL) systems provide vertical profiles of optical depth, backscatter cross-section, depolarization, and backscatter phase function. All HSRL measurements are calibrated by reference to molecular scattering, which is measured at each point in the lidar profile. Like the Raman lidar (RL), but unlike simple backscatter lidars such as the micropulse lidar (MPL), this enables the HSRL to measure backscatter cross-sections and optical depths without prior assumptions about the scattering properties of the atmosphere. The depolarization observations allow robust discrimination between ice and water clouds. Rigorous error estimates can be computed for all measurements. A very narrow angular field of view reduces multiple scattering contributions. The small field of view, coupled with a narrow optical bandwidth, nearly eliminates noise due to scattered sunlight. The laser transmitter is a diode-pumped, frequency-doubled Nd:YAG laser. Narrow-band, single-frequency operation is provided by injection seeding with a single-frequency, cw-diode-pumped diode laser. The main laser cavity is maintained in resonance with the seed laser by adjusting the cavity length to minimize the time between the Q-switch trigger and the emission of the laser pulse. The emission wavelength is tuned via temperature control of the seed laser crystal and is locked to line #1109 of the iodine absorption spectra. Locking is accomplished by minimizing the transmission through a 2-cm-long iodine absorption cell. Use of a high-repetition-rate laser and expansion of the transmitted beam through a 400-mm telescope reduces the transmitted energy density to eye-safe levels. It is possible to look directly into the output beam without hazard. The receiver and transmitter use the same afocal telescope, simplifying the maintenance of stable alignment of the transmitter and receiver although the angular FOV is only 100 μrad. The small FOV and the 4-kHz repetition rate also limit the near-field signal strength, making it possible to record continuous profiles that start at an altitude of ~100 m and extend to 30 km using photon counting detectors. The small FOV also suppresses multiple scattering contributions.

54 ENVIRONMENTAL SCIENCES↗

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, ΛCDM). Using the latest state-of-the-art cosmological ΛCDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass M * > 10 10 M ⊙ at redshift z = 0. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these ΛCDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at ≥12.52σ (TNG50-1) and ≥14.82σ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the ΛCDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an 8 5 times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the 5.58σ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in ΛCDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗

Multifrequency radar observations of clouds and precipitation including the G-band

Abstract. Observations collected during the 25 February 2020 deployment of the Vapor In-Cloud Profiling Radar at the Stony Brook Radar Observatory clearly demonstrate the potential of G-band radars for cloud and precipitation research, something that until now was only discussed in theory. The field experiment, which coordinated an X-, Ka-, W- and G-band radar, revealed that the Ka–G pairing can generate differential reflectivity signal several decibels larger than the traditional Ka–W pairing underpinning an increased sensitivity to smaller amounts of liquid and ice water mass and sizes. The observations also showed that G-band signals experience non-Rayleigh scattering in regions where Ka- and W-band signal do not, thus demonstrating the potential of G-band radars for sizing sub-millimeter ice crystals and droplets. Observed peculiar radar reflectivity patterns also suggest that G-band radars could be used to gain insight into the melting behavior of small ice crystals. G-band signal interpretation is challenging, because attenuation and non-Rayleigh effects are typically intertwined. An ideal liquid-free period allowed us to use triple-frequency Ka–W–G observations to test existing ice scattering libraries, and the results raise questions on their comprehensiveness. Overall, this work reinforces the importance of deploying radars (1) with sensitivity sufficient enough to detect small Rayleigh scatters at cloud top in order to derive estimates of path-integrated hydrometeor attenuation, a key constraint for microphysical retrievals; (2) with sensitivity sufficient enough to overcome liquid attenuation, to reveal the larger differential signals generated from using the G-band as part of a multifrequency deployment; and (3) capable of monitoring atmospheric gases to reduce related uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗

An Asymptotic Preserving Discontinuous Galerkin Method for a Linear Boltzmann Semiconductor Model

A key property of the linear Boltzmann semiconductor model is that as the collision frequency tends to infinity, the phase space density $f$ = $f$ ($x, v, t$) converges to an isotropic function $M (v)$$ρ$$(x, t)$, called the drift-diffusion limit, where $M$ is a Maxwellian and the physical density $ρ$ satisfies a second-order parabolic PDE known as the drift-diffusion equation. Numerical approximations that mirror this property are said to be asymptotic preserving. In this paper we build a discontinuous Galerkin method to the semiconductor model, and we show this scheme is both uniformly stable in $ε$, where 1/$ε$ is the scale of the collision frequency, and asymptotic preserving. Here in particular, we discuss what properties the discrete Maxwellian must satisfy in order for the schemes to converge in $ε$ to an accurate $h$-approximation of the drift-diffusion limit. Discrete versions of the drift-diffusion equation and error estimates in several norms with respect to $ε$ and the spacial resolution are also included.

97 MATHEMATICS AND COMPUTING↗

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

ARPA-E PNDIODES "PN Junctions by Ion Implantation" Project (Final Report)

Electricity accounts for nearly 40% of the energy generated in the United States, of which 65% is lost to electrical system inefficiency before reaching end-use sectors. Power electronics will play a growing role in this distribution and end-use, with up to 80% of electricity estimated to pass through power electronics from generation to consumption by 2030. This gives technical improvements in power electronics enormous potential to drive efficiency gains throughout the U.S. economy. Current Si-based power architectures including MOSFETs and IGBTs are limited by high losses, low switching frequency, and poor high-temperature performance. Given the material properties of GaN, devices can exhibit a given turn-on resistance and breakdown voltage with a minute fraction of the area of a Si device—all with higher switching speeds. While n-type doping of III-nitrides has reached a high level of maturity and integration into commercialized products, p-type doping presents distinct challenges that requires continued research. The success of the next generation of GaN-based power devices depends in large part on the ability to form selective-area p-type regions, which can be accomplished by ion implantation of the prevailing acceptor dopant, Mg. This approach is relatively inexpensive and is a fab-ready processing step for creation of current-generation power electronic devices including BJTs, IGBTs, and diode rectifiers. However, ion implantation induces lattice damage and creates point defects within the material, which negatively impacts dopant activation annealing. These forms of damage may be removed by high-temperature annealing, however short annealing time and higher pressures are required to prevent decomposition of the material. We approach this challenge through an innovative combination of approaches toward obtaining high-conductivity p-type GaN. By the novel employment of a gyrotron, a microwave source that can be shaped into a highly localized beam to dissipate very high power in a semiconductor target, a wafer of GaN can be annealed in short (<2 s) ‘pulses’ at temperatures of 1450 °C (2640 °F) to achieve p-type activation without appreciable surface degradation. The p-type activation of Mg-implanted GaN can be further aided by co-implantation with N, believed to enhance the removal of point defects that compensate the acceptors dopant. Elevation to such high temperatures for even a short time induces degradation of GaN-to overcome, we have implemented additional measures to retain high-quality crystalline GaN while removing implant damage and activating Mg dopant. A bilayer cap of AlN is deposited by a combination of chemical and physical vapor deposition to protect the GaN surface. To prevent decomposition of GaN loss of N during annealing, a custom vessel—to contain pressurized N while allowing injection of microwave energy—is employed. By detailed study of the interaction of various measures that are implemented to protect the sample surface (GaN), activate the implanted dopants, and mitigate damage of the lattice heating damage during annealing, we are able to achieve p-type GaN and observe recovery of the crystal and removal of substantial density of defects created during the implantation process.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Abstract Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss .

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the Sources and Transport Processes During Extreme Ammonia Episodes in the U.S. Corn Belt

Abstract Atmospheric ammonia (NH 3 ) is the primary form of reactive nitrogen (N r ) and a precursor of ammonium (NH 4 + ) aerosols. Ammonia has been linked to adverse impacts on human health, the loss of ecosystem biodiversity, and plays a key role in aerosol radiative forcing. The midwestern United States is the major NH 3 source in North America because of dense livestock operations and the high use of synthetic nitrogen fertilizers. Here, we combine tall‐tower (100 m) observations in Minnesota and Weather Research and Forecasting model coupled with Chemistry (WRF‐Chem) modeling to investigate high and low NH 3 emission episodes within the U.S. Corn Belt to improve our understanding of the distribution of emission sources and transport processes. We examined observations and performed model simulations for cases in February through November of 2017 and 2018. The results showed the following: (1) Peak emissions in November 2017 were enhanced by above‐normal air temperatures, implying a Q 10 (i.e., the change in NH 3 emissions for a temperature increase of 10°C) of 2.5 for emissions. (2) The intensive livestock emissions rom northern Iowa, approximately 400 km away from the tall tower, accounted for 17.6% of the abundance in tall‐tower NH 3 mixing ratios. (3) Ammonia mixing ratios in the innermost domain 3 frequently (i.e., 336 hr, 48% of November 2017) exceeded 5.3 ppb, an important air quality health standard. (4) In November 2017, simulated NH 3 net ecosystem exchange (the difference between NH 3 emissions and dry deposition) accounted for 60–65% of gross NH 3 emissions for agricultural areas and was 2.8–3.1 times the emissions of forested areas. (5) We estimated a mean annual NH 3 net ecosystem exchange of 1.60 ± 0.06 nmol · m −2 · s −1 for agricultural lands and −0.07 ± 0.02 nmol · m −2 · s −1 for forested lands. These results imply that future warmer fall temperatures will enhance agricultural NH 3 emissions, increase the frequency of dangerous NH 3 episodes, and enhance dry NH 3 deposition in adjacent forested lands.

Hu, Cheng↗

Identifiability and characterization of transmon qutrits through Bayesian experimental design

Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively.

97 MATHEMATICS AND COMPUTING↗

Analysis of Power-Maximizing Region 2 Controllers for Wind and Marine Turbines

Wind and marine energy are rapidly growing and complementary technologies that share some techniques for simplified modeling and control, particularly in below-rated flow speeds. A turbine operator has several choices of controller for maximizing power in Region 2. The simple and ubiquitous KΩ 2 control law is often effective but limited in its flexibility. Alternative controllers use reference tracking to split the control objectives into a low-bandwidth optimal tip-speed ratio tracking loop to maximize steady-state power and a higher-bandwidth proportional-integral control loop to reject inflow turbulence. Several options exist for identifying the slowly varying optimal set point during operation, based on estimating the inflow velocity or filtering the power or torque signals. This study compares the trade-offs between performance and other design priorities for a few choices of reference-tracking controller in the literature for reference wind and marine turbines. Analysis is performed in the frequency domain using the linearization of each controller, and the impact of turbulent disturbances on the closed-loop system is described. The controllers are simulated in OpenFAST to analyze their performance with higher-order nonlinear turbine dynamics.

17 WIND ENERGY↗

Seismoacoustic analysis of underground explosions using the Rayleigh integral

SUMMARY A coupled seismoacoustic model is developed for the analysis of acoustic signals produced by underground explosive events with an aim to develop a means of improving estimated depth and yield for explosion monitoring. A ground spall model is used to predict surface motion characteristics produced by an underground explosion and the Rayleigh integral is applied to relate the surface motion to the acoustic signal some distance from surface ground zero. The low-frequency component of the ground motion associated with the prolonged free fall of lofted material during spall is found to dominate the acoustic signal propagating away from surface ground zero at shallow angles. The model is applied to observed ground motion and acoustic signals recorded during the Source Physics Experiment (SPE) with promising results. In addition to accurately predicting characteristics of the observed acoustic signals from several of the SPE events, the model provides a means of explaining the lack of signals observed during several events in the SPE due to the directionality of the higher frequency acoustic signals associated with the uplift and closure components as well as the lack of a spatially localized, longer duration dwell.

58 GEOSCIENCES↗

Machine Learning Enabled Position Detection for 6.78 MHz UAV Wireless Power Transfer System

This paper presents a novel supervised machine learning (SML) approach for accurate position detection of the receiver coil in wireless power transfer (WPT) systems using only secondary-side electrical measurements, with applications in autonomous unmanned aerial vehicle (UAV) charging. The proposed method trains a supervised learning model to map measured secondary-side voltage and current features to the receiver’s spatial position with high precision. This enables an autonomous UAV to determine its location relative to the primary coil center, the optimal position for maximizing wireless charging efficiency. The sensing method is fully integrated into a standard WPT system, utilizing the same primary and secondary coils for both power transfer and position detection, thereby eliminating additional sensing hardware. The use of a 6.78 MHz operating frequency enhances positional sensitivity, as high-frequency near-field electromagnetic fields respond strongly to small spatial variations. Experimental validation is performed on a 30 W scaled prototype featuring a 210 mm × 140 mm primary coil, a 50 mm × 80 mm receiver coil, and a 15 mm air gap. Results demonstrate reliable position estimation and a strong correlation between predicted position and optimal coil alignment. This integrated framework unifying position detection and wireless charging offers a promising foundation for future autonomous electric vertical takeoff and landing (eVTOL) systems, enabling compact, hardware-efficient, and high-accuracy charging solutions.

Colak, Kerim [New York University]↗

Full‐Waveform Simulation of Infrasound Propagation in the Atmosphere: A Case Study of the 2023 April 20 SpaceX Starship Explosion

Infrasound, low‐frequency sound below 20 Hz, has been a key technology to monitor explosion events in the atmosphere. The International Monitoring System (IMS) of the Comprehensive Nuclear‐Test‐Ban Treaty Organization provides the means for continuous monitoring of infrasonic events worldwide. Infrasonic techniques for event location and size estimation can also complement other observational techniques for the detection and characterization of the entry of asteroids or large meteoroids. In this study, we describe the detection capability of IMS infrasound stations for an explosive event in the middle of the atmosphere. Full‐waveform simulations are performed with the specification of atmospheric conditions and incorporated into the event location and explosion yield estimation. We applied it to the 2023 April 20 SpaceX Starship explosion at 29 km altitude. Starship is a super heavy‐lift space vehicle constructed by SpaceX and known as the largest and most powerful rocket ever built. The Starship explosion created huge pressure disturbances in the atmosphere, and its infrasound was detected by the IMS arrays in North America between 2000 and 4000 km. Independent observational data and available ground‐truth information provide a rare opportunity to evaluate the monitoring capability of the IMS network for elevated sources in the atmosphere. We also demonstrate the capability of full‐waveform simulation for infrasound wavefield characterization and prediction to improve event location and yield estimation.

Geosciences↗