Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Clutter Distributions for Tomographic Image Standardization in Ground-Penetrating Radar

Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce 3-D distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude and vary across different arrays, which makes a direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or log-normal distributions. Based on the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter- to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Lastly, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development that needs to be robust across multiple different arrays.

42 ENGINEERING↗
DOE OSTI2025

Historical and Future Windstorms in the Northeastern United States

Large-scale windstorms represent an important atmospheric hazard in the Northeastern US (NE) and are associated with substantial socioeconomic losses. Regional simulations performed with the Weather Research and Forecasting (WRF) model using lateral boundary conditions from three Earth System Models (ESMs: Geophysical Fluid Dynamics Laboratory (GFDL), Hadley Centre Global Environment Model (HadGEM) and Max Planck Institute (MPI)) are used to quantify possible future changes in windstorm characteristics and/or changes in the parent cyclone types responsible for windstorms. WRF nested within MPI ESM best represents important aspects of historical windstorms and the cyclone types responsible for generating windstorms compared with a reference simulation performed with the ERA-Interim reanalysis for the historical climate. The spatial scale and frequency of the largest windstorms in each simulation defined using the greatest extent of exceedance of local 99.9th percentile wind speeds (U > U999) plus 50-year return period wind speeds (U50,RP) do not exhibit secular trends. Projections of extreme wind speeds and windstorm intensity/frequency/geolocation and dominant parent cyclone type associated with windstorms vary markedly across the simulations. Only the MPI nested simulations indicate statistically significant differences in windstorm spatial scale, frequency and intensity over the NE in the future and historical periods. This model chain, which also exhibits the highest fidelity in the historical climate, yields evidence of future increases in 99.9th percentile 10 m height wind speeds, the frequency of simultaneous U > U999 over a substantial fraction (5–25%) of the NE and the frequency of maximum wind speeds above 22.5 ms−1. These geophysical changes, coupled with a projected doubling of population, leads to a projected tripling of a socioeconomic loss index, and hence risk to human systems, from future windstorms.

Pryor, Sara C. (ORCID:0000000348473440)↗
DOE OSTI2020

The Value of Information From Horizontal Distributed Acoustic Sensing Compared to Multicomponent Geophones Via Machine Learning

Abstract Faults play an important role in recharging many geothermal reservoirs, and seismic information can image the locations of these faults. The value of information (VOI) metric is used to objectively quantify and compare the value of two types of seismic receiver data via a machine learning approach. The demonstrated VOI methodology is novel by including spatial models from seismic data and obtaining the information statistics from machine learning. Our two-dimensional numerical experiments compare images created from sparsely spaced (80 m), two-component geophone sampling to high spatial resolution (1 m), single-component DAS. We used a three-fold cross validation of a U-Net convolutional neural networks to achieve average classification statistics. The results suggest that when horizontal sources are utilized, geophones and DAS identify reflectors and non-reflectors at roughly the same rate. The average F1 score for horizontal DAS is 0.939 and 0.931 for geophones. For images created from a vertical source, DAS performed marginally better (F1 = 0.919) than geophones (F1 = 0.877). Our transferrable methodology can provide guidance on which acquisition scenarios can improve images of important structures in the subsurface and present an efficient method for obtaining reliability statistics from high-dimensional, spatial data.

Energy & Fuels↗
DOE OSTI2020

How “hot” are hotspots: Statistically localizing the high-activity areas on soil and rhizosphere images

The topic of microbial hotspots in soil requires not only visualizing their spatial distribution and biochemical analyses, but also statistical approaches to identify these hotspots and separate them from the surrounding activities (background). We hypothesized that each hotspot type (e.g. enzyme activities in the rhizosphere, root exudation, localization of herbicide accumulation) is a result of local process driven by biotic and/or abiotic factors, and the process rates in the hotspots are much faster than those in the soil background. We further hypothesized that the background and hotspot activities in soil belong to different statistical distributions. Consequently, hotspot determination should be based on statistical separation of activities significantly higher than the background. We analyzed for the statistical distributions of grey values on three groups of published images: 1) 14 C images of carbon input by roots into the rhizosphere, 2) 14 C glyphosate accumulation in the plant, and 3) zymogram of leucine aminopeptidase activity in rooted soil. The two Gaussian distributions were fit (the first representing the background, the second the hotspots) to the distribution of grey values in the images, the parameters (means and standard deviations, SD) of the fitted distributions were calculated, and the background was removed. Thus, we identified hotspots as areas outside of the Mean+2SD image intensity (corresponding to the upper ~ 2.5% of activity, being over 97.5% of background values) and finally, visualized images of solely hotspot locations. Finally, these results were compared with previously used decisions on hotspot intensity thresholding (i.e. Top-25% and 17 standard thresholding approaches in ImageJ) and discussed the advantages of the Mean+2SD as well as Mean+3SD approaches. These advantages include: i) simple unification of the thresholding approach for several imaging methods with various principles of activity distribution, ii) identification of hotspots with various activity levels, iii) analysis of “time-specific” hotspots in temporal sequences of images. Compared with 17 standard thresholding methods, we concluded that objectively elucidating and separating the hotspots should be based on statistical distribution analysis, e.g. using the Mean+2SD or Mean+3SD approaches. Furthermore, this simple Mean+2SD approach delivered suitable results for three groups of images and so, helps to understand the processes responsible for the highest activities and elucidate hotspots.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Exploiting stochastic locality in lattice QCD: hadronic observables and their uncertainties

Abstract Because of the mass gap, lattice QCD simulations exhibit stochastic locality: distant regions of the lattice fluctuate independently. There is a long history of exploiting this to increase statistics by obtaining multiple spatially-separated samples from each gauge field; in the extreme case, we arrive at the master-field approach in which a single gauge field is used. Here we develop techniques for studying hadronic observables using position-space correlators, which are more localized, and compare with the standard time-momentum representation. We also adapt methods for estimating the variance of an observable from autocorrelated Monte Carlo samples to the case of correlated spatially-separated samples.

Physics↗
DOE OSTI2022

A practical guide to characterizing irradiated nuclear fuels using FIB tomography

We report focused ion beam (FIB) tomography with combined electron backscatter diffraction (EBSD) and energy dispersive x-ray spectroscopy (EDS) is a technique capable of statistically characterizing the microstructure and spatial compositional variation of nuclear fuel in three-dimensions (3D). The 3D visualization from FIB tomography provides a comprehensive picture of the interconnected microstructural and compositional features that can impact fuel performance. While these features are often characterized with surface examination, the complexity and relationship of fission products and grain boundary networks may not fully be captured by these 2D methods. This work presents a practical guide to FIB tomography that is tailored to nuclear fuel characterization. The steps used to collect and process the data are provided along with the scripts used to process the data. Additionally, suggestions for future characterization efforts utilizing this approach are given.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A Remote Sensing Technique to Upscale Methane Emission Flux in a Subtropical Peatland

Abstract Quantification of methane (CH 4 ) gas emission from peat is critical to understand CH 4 budget from natural wetlands under a climate warming scenario. Previous studies have focused on prediction and mapping of CH 4 emission flux using process‐based models, while application of statistical‐empirical models for upscaling spatially sparse in situ measurements is scarce. In this study, we developed an empirical remote sensing upscaling approach to estimate CH 4 emission flux in the Everglades using limited in situ point‐based CH 4 emission flux measurements and Landsat data during 2013–2018. We spatially and temporally linked in situ data with Landsat surface reflectance based on temporally composite data sets and developed an object‐based machine learning framework to model and map CH 4 emission flux. An ensemble analysis of two machine learning models, k ‐Nearest Neighbor ( k ‐NN) and Support Vector Machine (SVM), shows that the upscaling approach is promising for predicting CH 4 emission flux with a R 2 of 0.65 and 0.87 based on a fivefold cross‐validation for a dry season and wet season estimation, respectively. We generated emission flux map products that successfully revealed the spatial and temporal heterogeneity of CH 4 emission within the dominant freshwater marsh ecosystem in the Everglades. We conclude that Landsat is promising for upscaling and monitoring CH 4 emission flux and reducing the uncertainty in emission estimates from wetlands.

Zhang, Caiyun↗
DOE OSTI2023

pnnl/frequency_sensitivity

This research code base includes functions to compute statistics of image distributions in (spatial DFT) frequency space, train deep learning image classifiers on these distributions with variable depth and weight decay, and finally measure the sensitivity of the trained models to perturbations along (spatial DFT) frequency components.

Central, PNNL Developer↗
DOE OSTI2024

The BIRDIES Experiment: Measuring Beryllium Isotopes to Resolve Dynamics in the Stratosphere

Cosmogenic beryllium-10 and beryllium-7, and the ratio of the two (10Be/7Be), are powerful atmospheric tracers of stratosphere–troposphere exchange (STE) processes; however, measurements are sparse for altitudes well above the tropopause. We present a novel high-altitude balloon campaign aimed to measure these isotopes in the mid-stratosphere called Beryllium Isotopes for Resolving Dynamics in the Stratosphere (BIRDIES). BIRDIES targeted gravity waves produced by tropopause-overshooting convection to study their propagation and impact on STE dynamics, including the production of turbulence in the stratosphere. Two custom-designed payloads called FiSH and GASP were flown at altitudes approaching 30 km to measure in situ turbulence and beryllium isotopes (on aerosols), respectively. These were flown on nine high-altitude balloon flights over Kansas, USA, in summer 2022. The atmospheric samples were augmented with a ground-based rainfall collection targeting isotopic signatures of deep convection overshooting. Our GASP samples yielded mostly negligible amounts of both 10Be and 7Be collected in the mid-stratosphere but led to design improvements to increase aerosol capture in low-pressure environments. Observations from FiSH and the precipitation collection were more fruitful. FiSH showed the presence of turbulent velocity, temperature, and acoustic fluctuations in the stratosphere, including length scales in the infra-sonic range and inertial subrange that indicated times of elevated turbulence. The precipitation collection, and subsequent statistical analysis, showed that large spatial datasets of 10Be/7Be can be measured in individual rainfall events with minimum terrestrial contamination. While the spatial patterns in rainfall suggested some evidence for overshooting convection, inter-event temporal variability was clearly observed and predicted with good agreement using the 3D chemical transport model GEOS-CHEM.

58 GEOSCIENCES↗
DOE OSTI2025

PaleoSTeHM v1.0: a modern, scalable spatiotemporal hierarchical modeling framework for paleo-environmental data

Abstract. Geological records of past environmental change provide crucial insights into long-term climate variability, trends, non-stationarity, and nonlinear feedback mechanisms. However, reconstructing spatiotemporal fields from these records is statistically challenging due to their sparse, indirect, and noisy nature. Here, we present PaleoSTeHM, a scalable and modern framework for spatiotemporal hierarchical modeling of paleo-environmental data. This framework enables the implementation of flexible statistical models that rigorously quantify spatial and temporal variability from geological data while clearly distinguishing measurement and inferential uncertainty from process variability. We illustrate its application by reconstructing temporal and spatiotemporal paleo-sea-level changes across multiple locations. Using various modeling and analysis choices, PaleoSTeHM demonstrates the impact of different methods on inference results and computational efficiency. Our results highlight the critical role of model selection in addressing specific paleo-environmental questions, showcasing the PaleoSTeHM framework's potential to enhance the robustness and transparency of paleo-environmental reconstructions.

58 GEOSCIENCES↗
DOE OSTI2026

Associated Particle Imaging of Neutron Inelastic Scatter: 3-D Reconstruction, Capabilities, and Challenges

Associated particle imaging (API) offers unique advantages for 3-D imaging of neutron inelastic scatter, including single-view tomographic imaging and data acquisition when access is limited to only one side of the interrogated object. However, widespread adoption of neutron inelastic scatter imaging has been impeded by several inherent challenges, most prominently spatial resolution, self-attenuation, and statistical noise. Here, in this work, the capabilities and challenges of neutron inelastic scatter imaging are investigated. Instead of focusing on a single imaging application, we identify imaging principles that hold for various neutron inelastic scatter imaging techniques. The primary challenges for 3-D imaging are characterized. The inherent spatial resolution in the time-of-flight (TOF) dimension is derived based on the known system timing resolution and scan geometry. Three reconstruction algorithms are described and demonstrated, including the introduction of modern iterative reconstruction incorporating a physics-based system model. Simulation is leveraged to demonstrate imaging capability with varying coincidence count levels and system timing resolution. An example of measured data with both back-scatter and forward-scatter detector positioning is presented. System design characteristics and their effects on image quality are identified. The imaging framework presented in this article has the potential to facilitate growth of 3-D neutron inelastic scatter API by identifying applications that are a good match for the technique and by targeting system development resources toward the requirements of a specific imaging task.

Associated particle imaging (API)↗
DOE OSTI2022

Superior Daily and Sub-Daily Precipitation Statistics for Intense and Long-Lived Storms in Global Storm-Resolving Models

Here, daily and sub-daily precipitation statistics are investigated from three global model ensembles: (a) global storm-resolving models (GSRMs) with typical horizontal resolutions of ~4 km, (b) “high”-resolution global models with typical resolutions of ~50 km and (c) “standard”-resolution global models with typical resolutions of ~100 km. Compared to two satellite rainfall datasets, GSRMs convincingly exhibit superior performance for statistics of heavier rain rate events including their diurnal cycle, spatial propagation and the amount contributed by intense precipitation, but not for statistics of weaker or shorter duration precipitation. Both high- and standard-resolution models fail to simulate the correct phase and amplitude of diurnal cycle of precipitation and the propagating convection in the Central US, but high-resolution models show relative improvement in the distribution of precipitation frequency and amount, especially for intense precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Good and proper: self-similarity of N -body simulations with proper force softening

ABSTRACT Analysis of self-similarity in scale-free N-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordinates. We find that a proper softening that reaches roughly 1/30th of the inter-particle spacing by the end of the simulation resolves the same spatial and temporal scales as a comoving softening of the same length while using a third fewer time-steps, for a range of scale factors typical to Lambda cold dark matter (ΛCDM) simulations. We additionally infer an inherent resolution limit, set by the particle mass and scaling as a−1/2, beyond which reducing the softening does not improve the resolution. We postulate a mapping of these results with spectral index n = −2 to ΛCDM simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence of galaxy assembly bias in SDSS DR7 galaxy samples from count statistics

We present observational constraints on the galaxy–halo connection, focusing particularly on galaxy assembly bias from a novel combination of counts-in-cylinders statistics, P(N CIC ), with the standard measurements of the projected two-point correlation function w p (r p ), and number density n gal of galaxies. We measure n gal , w p (r p ), and P(N CIC ) for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the M r < -20.0 and M r <-19.5 samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the M r < -20.5 and M r < -19.0 samples. We find no evidence of central assembly bias in the M r < -21.0 sample. We observe only a marginal preference for negative satellite assembly bias in the M r < -20.0 and M r < -19.0 samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy–halo connection studies and cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Statistical upscaling of ecosystem CO 2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties

Abstract The regional variability in tundra and boreal carbon dioxide (CO 2 ) fluxes can be high, complicating efforts to quantify sink‐source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO 2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO 2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990–2015 from 148 terrestrial high‐latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO 2 fluxes and test the accuracy and uncertainty of different statistical models. CO 2 fluxes were upscaled at relatively high spatial resolution (1 km 2 ) across the high‐latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE‐focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO 2 sink strength was larger in the boreal biome (observed and predicted average annual NEE −46 and −29 g C m −2 yr −1 , respectively) compared to tundra (average annual NEE +10 and −2 g C m −2 yr −1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO 2 budget, estimated using the annual NEE ensemble prediction, suggests the high‐latitude region was on average an annual CO 2 sink during 1990–2015, although uncertainty remains high.

Virkkala, Anna‐Maria↗
DOE OSTI2024

Understanding and Predicting the Spatially Resolved Adsorption Properties of Nanoporous Materials

Using knowledge from statistical thermodynamics and crystallography, we develop an image–image translation model, called SorbIIT, that uses three-dimensional grids of adsorbate–adsorbent interaction energies as input to predict the spatially resolved loading surface of nanoporous materials over a broad range of temperatures and pressures. SorbIIT consists of a closed-form differential model for loading-surface prediction and a U-Net to generate spatial differential distributions from the energy grids. SorbIIT is trained using the energy grids and adsorbate distributions (obtained from high-throughput simulations) of 50 synthesized and 70 hypothetical zeolites and applied for predicting the adsorption of carbon dioxide, hydrogen sulfide, n-butane, 2-methylpropane, krypton, and xenon in other zeolites from 256 to 400 K. In conclusion, employing a quadratic isotherm model for the local differentiation, SorbIIT yields mean R 2 values of 0.998 for total adsorption and 0.6904 for local adsorption with a resolution of 0.2 Å, and a value of 0.721 for the structural similarity of the local loading distribution.

Sun, Yangzesheng↗
DOE OSTI2021

The first Hubble diagram and cosmological constraints using superluminous supernovae

In this work, we present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, Ω M , and the dark energy equation-of-state parameter, w(≡p/ρ). We build a sample of 20 cosmologically useful SLSNe I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak–decline SLSN I standardization relation with a larger data set and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardization via minimization of the χ2 computed from a covariance matrix that includes statistical and systematic uncertainties. For a spatially flat Λ cold dark matter (ΛCDM) cosmological model, we find $\Omega _{\rm M}=0.38^{+0.24}_{-0.19}$, with an rms of 0.27 mag for the residuals of the distance moduli. For a w 0 w a CDM cosmological model, the addition of SLSNe I to a ‘baseline’ measurement consisting of Planck temperature together with Type Ia supernovae, results in a small improvement in the constraints of w 0 and wa of 4 per cent. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat ΛCDM model with the same precision (considering only statistical uncertainties) as current surveys that use Type Ia supernovae, while providing a factor of 2–3 improvement in the precision of the constraints on the time variation of dark energy, w 0 and w a . This paper represents the proof of concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift (z > 1) universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Processing-dependent chemical ordering in Cu 3 Au characterized via non-destructive Bragg coherent diffraction imaging

Of current importance for alloy design is controlling chemical ordering through processing routes to optimize an alloy's mechanical properties for a desired application. However, characterization of chemical ordering remains an ongoing challenge, particularly when nondestructive characterization is needed. Here, in this study, Bragg coherent diffraction imaging is used to reconstruct morphology and lattice displacement in model Cu 3 Au nanocrystals that have undergone different heat treatments to produce variation in chemical ordering. The magnitudes and distributions of the scattering amplitudes (proportional to electron density) and lattice strains within these crystals are then analyzed to correlate them to the expected amount of chemical ordering present. Nanocrystals with increased amounts of ordering are found to generally have less extreme strains present and reduced strain distribution widths. In addition, statistical correlations are found between the spatial arrangement of scattering amplitude and lattice strains.

Warren, Nathaniel [Pennsylvania State Univ., Unive↗