SEARCH · Search NASA

Results for “trace analysis”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Advanced Simulation of ITER Core X-ray Crystal Spectroscopy

X-Ray Simulation Analysis (XRSA) is an analytical ray-tracing mixed code developed specifically for the ITER Core X-Ray Crystal Spectroscopy (XRCS-Core) diagnostic, which employs a dual-reflection configuration incorporating multiple pre-reflectors made of Highly Oriented Pyrolytic Graphite (HOPG) and spherically curved analyzing crystals. The ITER XRCS-Core is designed for high spectral resolution measurement in specific wavelength ranges, including narrow bands around 1.354 Å for W 64+ , 2.19 Å for Xe 51+ , and 2.555 Å for Xe 44+ and Xe 47+ , enabling diagnostic capability across a broad electron temperature range in the ITER plasma. XRSA facilitates efficient simulation of the spectral performance of this complex X-ray spectroscopic system. Recent updates to the XRSA code have incorporated two critical effects: auto-focusing, which specifically applies to HOPG, and polarization. These two effects are particularly important in the dual-reflection configuration used in the ITER XRCS-Core system to provide more accurate modeling results. Here, simulations conducted with the updated code demonstrate that polarization has a substantial impact on the performance of the dual-reflection system. Additionally, the combined influence of polarization and system layout introduces performance variations across channels through the same crystal.

Computer graphics↗

DOE OSTI2026

Elevating SolTrace's Capabilities for the Next Generation of Concentrating Solar Analysis

SolTrace is an open-source Monte Carlo ray tracing software developed at NREL. SolTrace can characterize concentrating solar thermal (CST) collector optical performance and is CST technology agnostic. Shown in Fig. 1, SolTrace is a foundational tool in NREL's CST system and component modeling suite. SolTrace's generic surface elements can flexibly model novel collector and receiver designs to predict spatial and temporal flux distributions - critical to understand for CST component design, performance prediction, and system integration. Since its initial development, SolTrace has over 1,650 references on Google Scholar, over 9,800 downloads since 2017, and has served the CST research and development community as a benchmark of 3rd party verification. SolTrace provides users with many options for defining surface shape and boundaries. However, SolTrace provides limited documentation which can result in a steep learning curve for new users. Additionally, SolTrace lacks the computational performance required to evaluate optical performance of a CST system over the course of a year and/or iteratively over design parameters in a timely manner. To address this, we are working towards a new release of SolTrace that enables increased computational throughput by implementing ray tracing acceleration structures and enabling GPU parallelization. Additionally, we are working to improve SolTrace's usability, accessibility, and maintainability by (1) automating solar position time-dependent simulation processes, (2) creating general CST collector templates of grouped elements, (3) updating the user interface to better visualize model inputs and outputs, and (4) creating a user support network through forums, "how to" videos, and documentation.

14 SOLAR ENERGY↗

DOE OSTI2024

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Microscale Thermophoresis (MST) as a Tool to Study Binding Interactions of Oxygen-Sensitive Biohybrids

Microscale thermophoresis (MST) is a technique used to measure the strength of molecular interactions. MST is a thermophoretic-based technique that monitors the change in fluorescence associated with the movement of fluorescent-labeled molecules in response to a temperature gradient triggered by an IR LASER. MST has advantages over other approaches for examining molecular interactions, such as isothermal titration calorimetry, nuclear magnetic resonance, biolayer interferometry, and surface plasmon resonance, requiring a small sample size that does not need to be immobilized and a high-sensitivity fluorescence detection. In addition, since the approach involves the loading of samples into capillaries that can be easily sealed, it can be adapted to analyze oxygen-sensitive samples. In this Bio-protocol, we describe the troubleshooting and optimization we have done to enable the use of MST to examine protein–protein interactions, protein–ligand interactions, and protein–nanocrystal interactions. The salient elements in the developed procedures include 1) loading and sealing capabilities in an anaerobic chamber for analysis using a NanoTemper MST located on the benchtop in air, 2) identification of the optimal reducing agents compatible with data acquisition with effective protection against trace oxygen, and 3) the optimization of data acquisition and analysis procedures. The procedures lay the groundwork to define the determinants of molecular interactions in these technically demanding systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Utilization of the LS-APGD microplasma/orbitrap-FTMS booster system for detection and isotopic analysis of neodymium nanoparticles

Detection and isotopic analysis of particle populations has seen rapid growth across several application areas, including environmental analysis, nuclear forensics, and food safety. The ability to characterize the particles' unique elemental and isotopic fingerprints could provide information related to formation, processing history, and transport. Regarding nuclear forensics, isotopic analysis of particles derived from diverse materials is often used as a tool to trace the origin and processing history. Mass spectrometric-based techniques currently used for particle population analysis often suffer from limited mass resolution, particularly when dealing with real-world samples that are affected by isobaric and polyatomic interferences from the matrix. To address these analytical challenges, we propose a novel method utilizing the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an ultrahigh resolution Orbitrap mass spectrometer, further enhanced with the FTMS X2T Booster data acquisition and processing unit. The FTMS Booster enables acquisition of extended transient times of up to 3 s, significantly improving mass resolution, thereby reducing or even eliminating the need for prior separation of isobaric or polyatomic interferences. Additionally, the detection of low-abundance isotopes was improved by increasing the signal-to-noise (S/N) ratio. As proof of concept, this study demonstrates the feasibility of the LS-APGD/Orbitrap-FTMS X2T Booster platform for direct analysis using a suspension of well-characterized ∼120 nm neodymium particles. The quality of the isotope ratios values obtained from a few hundred particles were in good agreement with those obtained from homogeneous ionic solutions. These results highlight the potential of the LS-APGD/Orbitrap platform for rapid, accurate, and interference-resilient isotope ratio analysis of particle populations without the need for dissolution and subsequent chemical separations, offering significant advantages for nuclear forensics, safeguards, and environmental applications. The effort here also points to further paths forward, hopefully towards single particle (SP) analysis using microplasma ionization and the ultrahigh resolving power of the Orbitrap mass analyzer.

FTMS X2T booster↗

DOE OSTI2024

Generic and ML Workloads in an HPC Datacenter: Node Energy, Job Failures, and Node-Job Analysis

HPC datacenters offer a backbone to the modern digital society. Increasingly, they run Machine Learning (ML) jobs next to generic, compute-intensive workloads, supporting science, business, and other decision-making processes. However, understanding how ML jobs impact the operation of HPC datacenters, relative to generic jobs, remains desirable but understudied. In this work, we leverage long-term operational data, collected from a national-scale production HPC datacenter, and statistically compare how ML and generic jobs can impact the performance, failures, resource utilization, and energy consumption of HPC datacenters. Our study provides key insights, e.g., ML-related power usage causes GPU nodes to run into temperature limitations, median/mean runtime and failure rates are higher for ML jobs than for generic jobs, both ML and generic jobs exhibit highly variable arrival processes and resource demands, significant amounts of energy are spent on unsuccessfully terminating jobs, and concurrent jobs tend to terminate in the same state. We open-source our cleaned-up data traces on Zenodo (https://doi. org/10.5281/zenodo.13685426), and provide our analysis toolkit as software hosted on GitHub (https://github.com/atlarge-research/2024-icpads-hpc-workload-characterization). This study offers multiple benefits for data center administrators, who can improve operational efficiency, and for researchers, who can further improve system designs, scheduling techniques, etc.

crossanalysis↗

DOE OSTI2026

Engineered plants for the production of the antioxidants arbutin and gallate

The shikimate pathway is a crucial metabolic route for the biosynthesis of numerous valuable chemicals. In this study, we engineered the shikimate pathway in plants via expression of microbial enzymes to produce the two important antioxidants gallate and arbutin. The engineered pathways utilize the aromatics protocatechuate and 4-hydroxybenzoate as metabolic intermediates. Through transient expression in Nicotiana benthamiana leaves, we first identified biosynthetic routes for the production of gallate from either chorismate or 3-dehydroshikimate. Gallate production was then achieved in Arabidopsis using a genetic background that overproduces protocatechuate and via expression of a mutated version of the 4-hydroxybenzoate hydroxylase PobA from Pseudomonas sp. Arbutin production was obtained in Arabidopsis using a genetic background that overproduces 4-hydroxybenzoate and via expression of the monooxygenase MNX1 from Candida parapsilosis. The best Arabidopsis transgenic lines accumulated gallate and arbutin in the range of 0.25 and 0.93 dry weight % (dwt%), respectively. Using sorghum for large-scale in planta production, the titers of gallate and arbutin produced from the intermediate 4-hydroxybenzoate reached 0.58 dwt% and 0.50 dwt%, respectively, in mature transgenic plants, surpassing levels typically observed in plants that naturally produce these compounds. Gallate and arbutin were readily extracted from plant tissues using methanol solvent. Analysis of extractive-free biomass showed only trace amounts of gallate and its precursors 4-hydroxybenzoate and protocatechuate crosslinked to cell walls, suggesting that they mainly occur as soluble conjugated forms stored in the vacuole. This study presents alternative synthesis routes using plant hosts for the eco-friendly production of gallate and arbutin.

Antioxidant↗

DOE OSTI2026

Novel Cell-Type-Specific Drought-Responsive Proteins in Root Tips of Field-Grown Perennial Switchgrass

The root-tip region of plants, including the root cap, forms the most basal terminal of the root and exhibits a high degree of cellular complexity in terms of morphology, cytological function, and interaction with environmental cues in the soil. Cells in this region follow a developmental trajectory, transitioning from stem cells to meristematic cells, and ultimately to fully differentiated cell types. However, our understanding of root-tip cell-type specific proteomic responses to abiotic stresses, such as drought, particularly under field conditions, remains limited. This study aimed to identify spatially resolved, cell type-specific proteomes in switchgrass (Panicum virgatum) root tips under drought stress. Root tips were collected from seven-year-old, field-grown switchgrass ‘Alamo’ plants excavated under both well-watered and long-term drought conditions. Cell type-specific proteins were identified using laser capture microdissection (LCM) coupled with nanoPOTS (Nanodroplet Processing in One Pot for Trace Samples) and nano-LC-MS proteomics analysis. Five distinct cell types were targeted: (1) cells in the quiescent center and stem cell niche (QuC), (2) protodermal epidermal cells (PEC) in the meristematic zone, (3) epidermal cells in the transition and elongation zones above the root cap (Epi), (4) peripheral root cap cells (PRC), forming 2–3 layers below the PEC and 1–2 layers above the root border cells, and (5) columella root cap cells (Col) comprising of the columella initials and a single underlying layer of cells undergoing active growth. Principal component analysis (PCA) revealed clear separation among the five targeted cell types, confirming distinct proteomic profiles. Proteins predominantly enriched in each cell type were linked to distinct cellular functions, with QuC cells showing involvement in chromosomal behavior, DNA replication, and mitosis—key processes for stem cell niche regulation. Drought stress resulted in alterations of proteostasis, as evidenced by significant decreases in ribosomal proteins and increases in protein synthesis inhibitors. Moreover, drought stress induced unique cell-type–specific proteins involved in phytohormone biosynthesis and signaling pathways, including auxin, cytokinin, and jasmonic acid. In particular, QuC cells were more highly enriched in proteins associated with DNA repair and mitotic processes. Metabolic pathways related to amino acids, carbohydrates, and lipids were differentially affected in a cell-type–dependent manner, whereas general stress-responsive proteins exhibited consistent changes across all five cell types. Overall, this study provides unique spatially resolved, cell-type-specific proteomic profiles in root tips, representing a significant advancement in our understanding of the cellular mechanisms underlying plant responses to drought stress in natural field conditions.

perennial grass↗

DOE OSTI2025

Preliminary Study on Fine-Grained Power and Energy Measurements on Grace Hopper GH200 with Open-Source Performance Tools

The increasing adoption of tightly integrated, heterogeneous architectures, combined with the slowdown of Moore’s law, has made application power and energy-driven optimizations critical to efficiently use high-performance computing systems. This paper introduces a newly developed open-source toolkit that seamlessly integrates the Linux real-time hardware monitoring program hwmon with the Performance Application Programming Interface and the Score-P performance measurement system, thereby enabling fine-grained power and energy measurements for high-performance computing applications. Our primary target platform is the Wombat test bed, which is a system based on the NVIDIA GH200 superchip. The toolkit can capture transient power peaks with high temporal resolution (50 ms) and, thanks to Score-P integration, can map power metrics to specific code regions, thereby providing actionable information on power-intensive operations and inefficiencies. The toolkit also provides a holistic view of both the power and the energy consumption of the entire GH200 superchip by covering all major components: the Grace CPU, the Hopper GPU, and the I/O subsystem. Experiments that use Locally Self-consistent Multiple Scattering, which is an application for first-principles calculations of materials developed at Oak Ridge National Laboratory, have demonstrated the tool’s ability to identify transient power spikes and uncover opportunities for energy-aware optimizations. Additionally, we introduce a Python-based utility for converting Open Trace Format 2 traces to Parquet format, thus enabling advanced data analysis for numerical integration methods applied to power data for accurate energy profiling.

Hernandez Mendoza, Oscar [ORNL] (ORCID:00000002538↗

DOE OSTI2024

Measurement of Stark-split beam and carbon charge exchange emissions for simultaneous B-field and temperature/rotation analysis at DIII-D

A set of two newly designed, single-channel Czerny–Turner spectrometers has been deployed at the DIII-D tokamak for measurements of the motional Stark effect (MSE) split beam emission and the C6+ (CVI) carbon charge exchange recombination (CER) emission at high spectral (δλ = 0.13 nm) and temporal (1–5 kHz) resolution. High throughput optics (f/# = 2.8) allow for good signal-to-noise at high time resolution using fast EMCCD detectors. The MSE emission allows for spectral fitting of the magnitude and direction of the local B-field, while the carbon emission yields local ion temperature and toroidal rotation information. To reduce so-called Doppler broadening of the MSE emission, a new channel-specific variable lens-masking approach has been developed. Experimental data collected from the 2023 DIII-D experimental campaign demonstrate the signal quality and instrument fidelity for both diagnostic measurements. Moreover, initial CER data analysis shows a clear evolution of the toroidal rotation during edge localized modes. Initial progress on the advanced MSE model, including a new validated ray-trace model of the DIII-D collection optics, is shown via sensitivity analysis.

Instruments & Instrumentation↗

DOE OSTI2024

Strontium isoscape of sub-Saharan Africa allows tracing origins of victims of the transatlantic slave trade

Abstract Strontium isotope ( 87 Sr/ 86 Sr) analysis with reference to strontium isotope landscapes (Sr isoscapes) allows reconstructing mobility and migration in archaeology, ecology, and forensics. However, despite the vast potential of research involving 87 Sr/ 86 Sr analysis particularly in Africa, Sr isoscapes remain unavailable for the largest parts of the continent. Here, we measure the 87 Sr/ 86 Sr ratios in 778 environmental samples from 24 African countries and combine this data with published data to model a bioavailable Sr isoscape for sub-Saharan Africa using random forest regression. We demonstrate the efficacy of this Sr isoscape, in combination with other lines of evidence, to trace the African roots of individuals from historic slavery contexts, particularly those with highly radiogenic 87 Sr/ 86 Sr ratios uncommon in the African Diaspora. Our study provides an extensive African 87 Sr/ 86 Sr dataset which includes scientifically marginalized regions of Africa, with significant implications for the archaeology of the transatlantic slave trade, wildlife ecology, conservation, and forensics.

Science & Technology - Other Topics↗

DOE OSTI2024

Reversible Electron-Beam Patterning of Colloidal Nanoparticles at Fluid Interfaces

The directed self-assembly of colloidal nanoparticles (NPs) using external fields guides the formation of sophisticated hierarchical materials but becomes less effective with decreasing particle size. As an alternative, electron-beam-driven assembly offers a potential avenue for targeted nanoscale manipulation, yet remains poorly controlled due to the variety and complexity of beam interaction mechanisms. Here, we investigate the beam–particle interaction of silica NPs pinned to the fluid–vacuum interface of ionic liquid droplets. In these experiments, scanning electron microscopy of the droplet surface resolves NP trajectories over space and time while simultaneously driving their reorganization. With this platform, we demonstrate the ability to direct particle transport and create transient, reversible colloidal patterns on the droplet surface. By tuning the beam voltage, we achieve precise control over both the strength and sign of the beam–particle interaction, with low voltages repelling particles and high voltages attracting them. This response stems from the formation of well-defined solvent flow fields generated from trace radiolysis of the ionic liquid, as determined through statistical analysis of single-particle trajectories under varying solvent composition. Altogether, electron-beam-guided assembly introduces a versatile strategy for nanoscale colloidal manipulation, offering new possibilities for the design of dynamic, reconfigurable systems with applications in adaptive photonics and catalysis.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

AnnualReportingWind [SWR-26-026]

AnnualReportingWind is a version-controlled repository used to generate the cost and performance results published in major wind energy reports such as the Cost of Wind Energy Review (COWER) and the Annual Technology Baseline (ATB). The repository consolidates the model configurations, assumptions, and scripts used in each reporting year so that published results can be transparently reproduced and traced back to the exact analytical setup used in the analysis. The workflows rely on several techno-economic modeling tools developed at the National Laboratory of the Rockies (NLR), including WAVES, ORBIT, LandBOSSE, WOMBAT, and FLORIS.

Mulas Hernando, Daniel [National Laboratory of the↗

DOE OSTI2026

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Predominant Proton Insertion during Electrochemical Cycling of ε-VOPO 4 in a Nonaqueous Ca-Ion Electrolyte

Calcium-ion batteries (CIBs) are promising for next-generation energy storage systems. However, many reported Ca-ion storage mechanisms lack direct structural evidence, raising concerns that the observed electrochemistry may instead originate from proton insertion. Herein, we investigate the electrochemical cycling of a VOPO 4 electrode in an aprotic Ca electrolyte with trace water to clarify the charge storage mechanism. While electrochemistry, elemental analysis, and X-ray absorption spectroscopy seem to support Ca-ion intercalation in the VOPO 4 electrode, X-ray diffraction and electron microscopy unequivocally show proton insertion as the predominant charge storage mechanism. The apparent variation of Ca content in the electrode with the state of charge is attributed to the formation and dissolution of Ca-containing precipitates on the surface of VOPO 4 particles. Our work demonstrates the propensity of proton participation in the intercalation reaction in Ca-ion batteries and the necessity of structural evidence in understanding reactions in multivalent electrochemical systems.

Wang, Jiwei [Binghamton University, NY (United Sta↗

DOE OSTI2026

Geometric optics model of Thomson scattering enhanced by parametric coupling

We formulate a theory of Thomson scattering (TS) from inhomogeneous plasmas using a ray-tracing approach. The scattered radiation is described within the eikonal approximation, including convective gain arising from scattering instabilities driven by the TS probe laser. In the limit of a weak probe and vanishing gain, the theory recovers the standard TS cross section and scattered power due to volumetric particle noise. The ray-tracing formalism is implemented in a TS fitting framework, enabling accurate analysis of inhomogeneous plasmas whose background profiles can be provided by large-scale hydrodynamic simulations or experimental diagnostics. Here, several examples of synthetic TS spectra are presented, and further extensions of the theory are discussed.

Carleton, D. E. [University of Alberta, Edmonton, ↗

DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗