Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deviated”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Autonomous Anomaly Detection For Continuous Streams

The code implements the Isolation Forest (IFML) algorithm within the digital twin (DT) of the AGN-201 nuclear reactor. The DT captures real-time operational data including control rod positions, reactor power, and temperature. The IFML model isolates anomalies by detecting patterns that deviate from expected operational behavior. The algorithm recursively partitions the data and assigns anomaly scores based on the isolation of rare and different events. By tuning parameters specific to the reactor’s operational data, the IFML identifies deviations such as unauthorized material insertions or reactor reactivity shifts. The system streams data using LabView and integrates with the DeepLynx data warehouse for anomaly processing.

Trevino, Eduardo↗
DOE OSTI2021

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How Time Inefficient and Uncertain are Paratransit Trips Compared to Car Trips

The Americans with Disabilities Act (ADA) stipulates that paratransit mobility be comparable to public transit. However, with just under 5% of the population relying on public transportation, transit is not a representative benchmark. To highlight this inequality, we compare trip times by paratransit to counterfactual car-based travel. Using 2.5 years of demand data from the Denver Regional Transportation District's Access-a-Ride service, we compare paratransit trip times to counterfactual car-based trips for the same origin-destination pairs, controlling for various trip characteristics. Compared to car travel, the findings indicate high variability and uncertainty associated with paratransit trip times. For the same origin-destination pairs, the mean paratransit trip time is almost twice that of the mean car travel time, with a standard deviation for paratransit trip time fourfold that of the standard deviation for car trip time. For perspective, traveling an average of 10 mi during the 7-8 a.m. morning peak can take about 16 min by car with almost no variability, while that same trip can take an average of 25 min by paratransit, with 5% of trips being outside the 95% confidence interval and thus unpredictable. Paratransit trip time inefficiency tends to be particularly worse for females; older adults; those making trips between 9 and 11 a.m.; cash-paying customers; those making shorter trips; and those traveling during inclement weather, including cold temperatures. These findings suggest a need to re-assess using public transit as a benchmark for paratransit supply as regulated by the ADA.

accessible transportation and mobility↗
DOE OSTI2025

Round Robin Analysis of Uranium Isotopics on Cotton Swipes Measured Using Microextraction-Based Sampling Versus Conventional Bulk Analysis

A comparison of microextraction sampling methods to directly analyze uranium isotopics on cotton swipes was performed concurrently with traditional bulk-processing methods. For the microextraction sampling approach, two different detection platforms were evaluated, a quadrupole-based inductively coupled plasma mass spectrometer (ICP-MS) and the liquid sampling-atmospheric pressure glow discharge (LS-APGD) coupled to an Orbitrap mass spectrometer. Results presented from this innovative sampling approach (i.e., microextraction) are compared with a more traditional approach employed for analysis of cotton-based environmental swipes, namely bulk ashing/digestion, separation, and subsequent analysis by high-precision multi-collector ICP-MS. Overall, the microextraction approach proved to be a reliable and accurate means to determine isotopic ratios of uranium collected on cotton swipes. The ICP-MS-based detection had relative standard deviations of <0.65% for the major isotopic determinations, whereas the LS-APGD-Orbitrap method had relative standard deviations of <3%. The percent relative difference for the 235 U/ 238 U ratios, in comparison to the expected values, was <1% for ICP-MS and <4% for the microplasma-Orbitrap method. Additionally, the microextraction ICP-MS accurately (<2%) and precisely (<5%) determined the minor isotopic compositions (i.e., 234 U/ 238 U and 236 U/ 238 U).

ICP-MS↗
DOE OSTI2020

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multigroup Examination of Nickel-Reflected HEU System

Critical configurations of HEU-MET-FAST-003 (HMF-003) from the International Criticality Safety Benchmark Evaluation Project Handbook were modeled for addition into the Verified, Archived Library of Inputs and Data maintained at Oak Ridge National Laboratory. HMF-003 contains models of a highly enriched uranium sphere with a spherical metal reflector of natural uranium, tungsten carbide, or nickel. For these simple models, good agreement is expected between the k eff calculated with continuous-energy (CE) and multigroup (MG) transport. However, using the CSAS5 sequence of the SCALE 6.2.4 package, the 8 in. thick nickel reflector model resulted in a difference of around 1.2 % Δ k . Results presented a clear indication of poor performance for the SCALE 252-group based on ENDF/B-VII.1 cross section library for this model. The bias was investigated over a range of nickel thickness and then compared with other MG libraries available in SCALE 6.2.4. As the nickel thickness increases, the MG k eff deviates from the CE result, confirming a performance issue with the MG calculation. Additionally, the reactions were compared between the two libraries to help determine the cause of k eff bias. Finally, additional MG libraries within SCALE 6.3 and based on ENDF/B-VII.1 and ENDF/B-VIII were reviewed. The 302-group and 1597-group result in significant improvements compared with the 252-group library. Nickel cross sections were modified from ENDF/B-VII.1 to ENDF/B-VIII. For the ENDF/B-VIII library, the deviation between the CE libraries and MG libraries is less pronounced than for ENDF/VII.1, except for the 1597-group. This examination clearly indicates the importance of validating calculations with applicable benchmark experiments and caution when using MG libraries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

36 Hz integral linewidth laser based on a photonic integrated 4.0 m coil resonator

Laser stabilization sits at the heart of many precision scientific experiments and applications, including quantum information science, metrology, and atomic timekeeping. Many of these systems narrow the laser linewidth and stabilize the carrier by use of Pound–Drever–Hall (PDH) locking to a table-scale, ultrahigh quality factor (Q), vacuum spaced Fabry–Perot reference cavity. Integrating these cavities to bring characteristics of PDH stabilization to the chip scale is critical to reducing their size, cost, and weight, and enabling a wide range of portable and system-on-chip applications. We report a significant advance in integrated laser linewidth narrowing, stabilization, and noise reduction by use of a photonic integrated 4.0 m long coil resonator to stabilize a semiconductor laser. We achieve a 36 Hz 1 / π <#comment/> -integral linewidth, Allan deviation of 1.8 × <#comment/> 10 − <#comment/> 13 at 10 ms measurement time, and a 2.3 kHz/s drift—to the best of our knowledge, the lowest integral linewidth and highest stability demonstrated for an integrated waveguide reference cavity. This performance represents over an order of magnitude improvement in integral linewidth and frequency noise over previous integrated waveguide PDH stabilized reference cavities and bulk-optic and integrated injection locked approaches, and over two orders of magnitude improvement in frequency and phase noise than integrated injection locked approaches. Two different wavelength coil designs are demonstrated, stabilizing lasers at 1550 nm and 1319 nm. The resonator is bus-coupled to a 4.0 m long coil, with a 49 MHz free spectral range, mode volume of 1.0 × <#comment/> 10 10 µ <#comment/> m 3 , and 142 million intrinsic Q , fabricated in a CMOS compatible, ultralow loss silicon nitride waveguide platform. Our measurements and simulations show that the thermorefractive noise floor for this particular cavity is reached for frequencies down to 20 Hz in an ambient environment with simple passive vibration isolation and without vacuum or thermal isolation. The thermorefractive noise limited performance is estimated to yield an 8 Hz 1 / π <#comment/> -integral linewidth and Allan deviation of 5 × <#comment/> 10 − <#comment/> 14 at 10 ms, opening a stability regime that heretofore has been available only in fundamentally non-integrated systems. These results demonstrate the potential to bring the characteristics of laboratory-scale stabilized lasers to the integrated, wafer-scale compatible chip scale, and are of interest for a number of applications in quantum technologies and atomic, molecular, and optical physics, and with further developments below 10 Hz linewidth, can be highly relevant to ultralow noise microwave generation.

Liu, Kaikai (ORCID:0000000213609646)↗
DOE OSTI2023

Structural aging of human neurons is opposite of the changes in schizophrenia

Human mentality develops with age and is altered in psychiatric disorders, though their underlying mechanism is unknown. In this study, we analyzed nanometer-scale three-dimensional structures of brain tissues of the anterior cingulate cortex from eight schizophrenia and eight control cases. The distribution profiles of neurite curvature of the control cases showed a trend depending on their age, resulting in an age-correlated decrease in the standard deviation of neurite curvature (Pearson’s r = -0.80, p = 0.018). In contrast to the control cases, the schizophrenia cases deviate upward from this correlation, exhibiting a 60% higher neurite curvature compared with the controls (p = 7.8 × 10 -4 ). The neurite curvature also showed a correlation with a hallucination score (Pearson’s r = 0.80, p = 1.8 × 10 -4 ), indicating that neurite structure is relevant to brain function. This report is based on our 3D analysis of human brain tissues over a decade and is unprecedented in terms of the number of cases. We suggest that neurite curvature plays a pivotal role in brain aging and can be used as a hallmark to exploit a novel treatment of schizophrenia.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗
DOE OSTI2026

Terrestrial laser scanning data (Levels 0 and 1) for Pasoh, Malaysia, Sep 2024

This data package contains data from terrestrial laser scanning (TLS) at the Pasoh Forest Reserve, Malaysia. The Pasoh Forest Reserve is a facility of the Forest Research Institute Malaysia, and contains evergreen lowland dipterocarp forest. The Next-Generation Ecosystem Experiments Tropics (NGEE-Tropics) study areas at Pasoh were established to study how different species respond to climatic variation and soil water availability. Two study areas were chosen representing different topography and species. The TLS data archived here were collected to provide detailed, three-dimensional information about forest structure. Specifically, data were collected to allow tree-level characterization of woody structure and leaf area for 12 focal trees with FloraPulse and sap flux sensors, facilitating estimation of woody biomass and leaf area to allow upscaling of water content and transpiration data to the tree-level. Scan positions were not selected to provide consistent data for non-focal trees with the study areas. This data package contains the following data: - High-level files document further details of the campaign and data package: 1_CampaignSummary.csv provides details about the campaign and study site, 2_ScanAreasDetail.csv provides details about each separate scan area (groups of scans post-processed into a single point cloud), 3_TerrestrialLidarSensor.csv provides further technical details about the Riegl VZ-400i TLS sensor, TLS_CSV_dd.csv is a CSV Data Dictionary providing information about the fields in CSV files following the ESS-DIVE CSV File Formatting Guidelines Reporting Format, TLS_flmd.csv is a File Level Metadata file providing information about each file in the data package following the ESS-DIVE File Level Metadata Reporting Format, and README.txt is a text file describing the overall project and file structure. - Level 0 data are the raw data (.PROJ folders) as recorded by the Riegl VZ-400i TLS instrument before scan co-registration and post-processing with the Riegl's proprietary RiSCAN PRO software, which requires a license. - Level 1 data contain post-processed, co-registered data from each scan area. The "PointClouds" folder for each scan area contains a .las file with 1 cm resolution point cloud data exported from RiSCAN PRO. These are the main files likely to be of interest to most users and can be further processed with any software capable of manipulating .las files (e.g. Python, R CloudCompare). The "Project Information" folder contains log files from post-processing in RiSCAN PRO that may be of interest to users who want to see detailed records of post-processing, including all PDF reports generated by RiSCAN PRO. The "ScanPositions" folder contains information about the final position of all TLS scans, after post-processing, in multiple formats. The file ScanPositions_*.csv provides final geo-referenced scan positions, and the file SOP_backup_*.csv can be used in RiSCAN PRO to restore the co-registered scan positions if users wish to re-process raw data (Level 0 .PROJ folders) with RiSCAN PRO software (e.g., subsample to a different resolution, exclude a certain scan position, or apply different filters on reflectance or deviation values) without redoing time-consuming co-registration steps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Terrestrial laser scanning data (Levels 0 and 1) from Urban Biogeochemistry Pilot Project sites, Knoxville, Tennessee, Jul 2024 - Jul 2025

This data package contains data from terrestrial laser scanning (TLS) at five urban park sites in Knoxville, Tennessee, USA. All parks include open-grown and/or closed-canopy trees and mixed nearby land use. These study sites were established as part of the Urban Biogeochemistry Pilot Project, which has an overall goal of better understanding how hydrobiogeochemical cycling is altered within the human environment. These five sites represent a gradient of urbanization, and were instrumented to understand hydrological and biogeochemical cycling (e.g., soil moisture, soil physical properties and biogeochemistry, tree transpiration, species type). The TLS data archived here were collected to provide detailed, three-dimensional information about forest structure. Specifically, data were collected to allow tree- and stand-level characterization of woody structure and leaf area. TLS scans were placed to capture the area around trees with sap flow sensors, and as much of a 50 m radius area around the meteorological station as possible given site property limits. Derived products will allow upscaling of water content and transpiration data. This data package contains the following data: - High-level files document further details of the campaign and data package: 1_CampaignSummary.csv provides details about the campaign and study site, 2_ScanAreasDetail.csv provides details about each separate scan area (groups of scans post-processed into a single point cloud), 3_TerrestrialLidarSensor.csv provides further technical details about the Riegl VZ-400i TLS sensor, TLS_CSV_dd.csv is a CSV Data Dictionary providing information about the fields in CSV files following the ESS-DIVE CSV File Formatting Guidelines Reporting Format, TLS_flmd.csv is a File Level Metadata file providing information about each file in the data package following the ESS-DIVE File Level Metadata Reporting Format, and README.txt is a text file describing the overall project and file structure. - Level 0 data are the raw data (.PROJ folders) as recorded by the Riegl VZ-400i TLS instrument before scan co-registration and post-processing with the Riegl's proprietary RiSCAN PRO software, which requires a license. - Level 1 data contain post-processed, co-registered data from each scan area. The "PointClouds" folder for each scan area contains a .las file with 1 cm resolution point cloud data exported from RiSCAN PRO. These are the main files likely to be of interest to most users and can be further processed with any software capable of manipulating .las files (e.g. Python, R CloudCompare). The "Project Information" folder contains log files from post-processing in RiSCAN PRO that may be of interest to users who want to see detailed records of post-processing, including all PDF reports generated by RiSCAN PRO. The "ScanPositions" folder contains information about the final position of all TLS scans, after post-processing, in multiple formats. The file ScanPositions_*.csv provides final geo-referenced scan positions, and the file SOP_backup_*.csv can be used in RiSCAN PRO to restore the co-registered scan positions if users wish to re-process raw data (Level 0 .PROJ folders) with RiSCAN PRO software (e.g., subsample to a different resolution, exclude a certain scan position, or apply different filters on reflectance or deviation values) without redoing time-consuming co-registration steps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Attribute-Guided Adversarial Training for Robustness to Natural Perturbations

We report while existing work in robust deep learning has focused on small pixel-level norm-based perturbations, this may not account for perturbations encountered in several real world settings. In many such cases although test data might not be available, broad specifications about the types of perturbations (such as an unknown degree of rotation) may be known. We consider a setup where robustness is expected over an unseen test domain that is not i.i.d. but deviates from the training domain. While this deviation may not be exactly known, its broad characterization is specified a priori, in terms of attributes. We propose an adversarial training approach which learns to generate new samples so as to maximize exposure of the classifier to the attributes-space, without having access to the data from the test domain. Our adversarial training solves a min-max optimization problem, with the inner maximization generating adversarial perturbations, and the outer minimization finding model parameters by optimizing the loss on adversarial perturbations generated from the inner maximization. We demonstrate the applicability of our approach on three types of naturally occurring perturbations --- object-related shifts, geometric transformations, and common image corruptions. Our approach enables deep neural networks to be robust against a wide range of naturally occurring perturbations. We demonstrate the usefulness of the proposed approach by showing the robustness gains of deep neural networks trained using our adversarial training on MNIST, CIFAR-10, and a new variant of the CLEVR dataset.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2020

A machine learning-based fast frequency response control for a VSC-HVDC system

An HVDC system can realize a very fast frequency response to the disturbed system under a contingency because its active power control is decoupled from the frequency deviation. However, most of existing HVDC frequency control strategies are coupled with system primary frequency control and secondary frequency control. Since the traditional system frequency control is dominated by the thermal generators, the advantage of the fast response of the HVDC system is not made fully used. The development of a frequency response estimation based on a machine learning algorithm provides another approach to improve the frequency response capability of the HVDC system. Different from other frequency deviation tracking strategies, a machine learning based HVDC frequency response control can directly increase the power flow of a HVDC system by estimation of the system generator or load lost. In this paper, a fast frequency response control using a HVDC system for a large power system disturbance based on the multivariate random forest regression (MRFR) algorithm is proposed. The simulation is carried out with an integrated power system model based on the North American interconnections. The simulation results indicate that the proposed MRFR based frequency response control can significantly improve the frequency low point during an event, while stabilizing the frequency in advance.

42 ENGINEERING↗
DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗
DOE OSTI2025

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3

Chang, Hilary↗
DOE OSTI2025

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗
DOE OSTI2023

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

We present the transverse coherence minimization method (TCM)—an approach to estimate the back–azimuth of infrasound signals that are recorded on an infrasound microphone and a colocated three–component seismometer. Accurate back–azimuth information is important for a variety of monitoring efforts, but it is currently only available for infrasound arrays and for seismoacoustic sensor pairs separated by 10 s of meters. Our TCM method allows for the analysis of colocated sensor pairs, sensors located within a few meters of each other, which may extend the capabilities of existing seismoacoustic networks and supplement operating infrasound arrays. This approach minimizes the coherence of the transverse component of seismic displacement with the infrasound wave to estimate the infrasound back–azimuth. After developing an analytical model, we investigate seismoacoustic signals from the August 2012 Humming Roadrunner experiment and the 26 May 2021 eruption of Great Sitkin Volcano, Alaska, U.S.A., at the ranges of 6.5–185 km from the source. We discuss back–azimuth estimates and potential sources of deviation (1°–15°), such as local terrain effects or deviation from common analytical models. This practical method complements existing seismoacoustic tools and may be suitable for routine application to signals of interest.

58 GEOSCIENCES↗