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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Comment on “Advanced Testing of Low, Medium, and High ECS CMIP6 GCM Simulations Versus ERA5-T2m” by N. Scafetta (2022)

Scafetta (2022, https://doi.org/10.1029/2022gl097716) purports to test Coupled Model Intercomparison Project Phase 6 (CMIP6) climate models through a comparison of temperature changes over three decades. Unfortunately, the paper contains numerous conceptual and statistical errors that undermine all of the conclusions. First, no uncertainty is given for the observational temperature difference, making it impossible to assess compatibility with any model result. Second, the CMIP6 data are the ensemble means for each model, but the metric being tested is sensitive to the internal variability and so the full ensemble for each model must be used. When this is corrected, the conclusion that “all models with ECS > 3.0°C overestimate the observed global surface warming” is not sustained. Third, the statistical test in Section 2 would reject all models even in a perfect model setup given sufficient ensemble members, thus the second conclusion “that spatial t-statistics rejects the data-model agreement” is also not sustainable.

CMIP6↗
DOE OSTI2025

Map-level baryonification: unified treatment of weak lensing two-point and higher-order statistics

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matter-only simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales ℓ < 2000 and tomographic bins for a DES-Y3-like redshift distribution ( z ≲ 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
NASA NTRS2013

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗
NASA NTRS2022

Development of a “Nature Run” for Observing System Simulation Experiments (OSSEs) for Snow Mission Development

Snow is a fundamental component of global and regional water budgets, particularly in mountainous areas and regions downstream that rely on snowmelt for water resources. Land surface models (LSMs) are commonly used to develop spatially distributed estimates of snow water equivalent (SWE) and runoff. However, LSMs are limited by uncertainties in model physics and parameters, among other factors. In this study, we describe the use of model calibration tools to improve snow simulations within the Noah-MP LSM as the first step in an observing system simulation experiment (OSSE). Noah-MP is calibrated against the University of Arizona (UA) SWE product over a western Colorado domain. With spatially varying calibrated parameters, we run calibrated and default Noah-MP simulations for water years 2010–20. By evaluating both simulations against the UA dataset, we show that calibration decreases domain averaged temporal RMSE and bias for snow depth from 0.15 to 0.13 m and from −0.036 to −0.0023 m, respectively, and improves the timing of snow ablation. Increased snow simulation performance also improves estimates of model-simulated runoff in four of six study basins, though only one has statistically significant improvement. Spatially distributed Noah-MP snow parameters perform better than default uniform values. We demonstrate that calibrating variables related to snow albedo calculations and rain–snow partitioning, among other processes, is a necessary step for creating a nature run that reasonably approximates true snow conditions for the OSSEs. Additionally, the inclusion of a snowfall scaling term can address biases in precipitation from meteorological forcing datasets, further improving the utility of LSMs for generating reliable spatiotemporal estimates of snow.

Melissa L. Wrzesien↗
NASA NTRS2016

Development and Evaluation of High-Resolution Climate Simulations Over the Mountainous Northeastern United States

The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.

climatology↗
NASA NTRS2018

The Use of Gridded Fossil Fuel CO2 Emissions (FFCO2) Inventory for Climate Mitigation Applications: Errors, Uncertainties, and Current and Future Challenges

Emission Inventory (EI) is a fundamental tool to monitor global compliance of greenhouse gases (GHGs) emissions reduction actions. Inventory guidelines provide a best practice to help EI compilers to make comparable national emission estimates, in spite of the differences in data availability across countries and regions. There are a variety of sources of errors and uncertainties, however, that originate beyond what the inventory guidelines define. For example, spatially-explicit EIs, which are a key product for atmospheric modeling applications, are often developed for research purposes, and there are no specific guidelines to disaggregate emission estimates from country scale. On top of that, EIs are fundamentally prone to systematic biases due to the simple calculation methodology and thus an objective evaluation (e.g. atmospheric top-down estimates) is needed to assure the accuracy of the estimates. ODIAC is a global high-resolution (1x1 km) fossil fuel carbon dioxide (CO2) gridded EI that is now often used in atmospheric CO2 modeling. ODIAC is based on disaggregation of national emission estimates made by CDIAC, which is the well accepted standard in the community. The ODIAC emission data product is updated on an annual basis using best available statistical data. Subnational spatial emission patterns are estimated using power plant profiles and satellite-observations of nighttime lights. In addition to the conventional CDIAC gridded data product, ODIAC carries international bunker emissions (shipping and aviation), which allows flux inversion modelers to accurately impose the global total fossil fuel emissions and their horizontal and vertical distribution. We have extensively evaluated ODIAC emissions using fine-grained EIs as well as a high-resolution atmospheric model simulation across different scales (national, subnational/regional, and urban policy relevant) with a focus on the uncertainties associated with the emission disaggregation. We have examined the use of NASA's Black Marble Suomi-NPP/VIIRS nightlight data.

Oda, Tomohiro↗
DOE OSTI2022

The complex structural and chemical nature of monolithic U-10Mo fuel and Zr barrier layer

Nanoscale microstructural characterization by advanced transmission electron microscopy techniques on a U-10Mo/Zr barrier layer monolithic fuel plate was performed to evaluate the microstructural evaluation after high burn-up. Gas bubble superlattice evolution, grain restructuring, and evolution of the Zr interaction layer is investigated through detailed electron microscopy characterization. The use of automated crystallographic orientation mapping to irradiated U-10Mo fuel highlights that the restructured ultra-fine grains are separated by high angle grain boundaries at a burn up of 4.42 × 10 21 fissions/cm 3 . Additionally, advanced chemical analysis and multi-variable statistical analysis shows spatial clustering of solid fission product precipitates. Finally, characterization of a newly observed porous nanocrystalline Zr region in the barrier layer is studied. Finally, this work provides insights into the grain subdivision and restructuring process while using advanced microscopy techniques to analyze fission products in neutron-irradiated U-10Mo fuel.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Scales of North Atlantic wind stress curl determined from the comprehensive ocean-atmosphere data set

Nineteen years of wind data over the North Atlantic are used to calculate a field of wind stress curl. An empirical orthogonal function (EOF) analysis is performed on this field, resulting in spatial patterns of wind stress curl and associated time series. A Monte Carlo technique is used to establish the statistical significance of each spatial pattern, and the associated time series are spectrally analyzed. The first four statistically significant EOF modes represent more than 50 percent of the curl variance, and the spatial patterns of curl associated with these modes exhibit the major elements of North Atlantic climatology. Most of the time series spectral variance is contained in annual and semiannual frequencies. The features observed include the individual annual variation of the subtropical high and the subpolar low, the annual oscillation of intensity between pressure centers, the influence of localized strong SST gradients and associated cyclogenesis regions, and the constant nature of the trades.

Ehret, Laura L.↗
DOE OSTI2024

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Light scattering by V 4 O 7 film across the metal–insulator transition

The experimental study of the angle-resolved hemispherical light scattering by V 4 O 7 film within a broad temperature range across metal–insulator transition reveals complex structural reorganization of the film deposited on the c-cut sapphire crystal. The bidirectional scattering distribution function and the surface autocorrelation function were obtained from scattering data to visualize statistics of the spatially resolved contributions of optical inhomogeneities in normal and lateral directions to the surface. The measurements reveal an anisotropic surface roughness distribution due to the twinned domain structure, with significant anisotropy changes across the phase transition. The V 4 O 7 film deposited on sapphire leads to a polydomain structure, minimizing elastic strain energy with distinct multiscale distributions of surface domains. Near T c , the material shows the lowest roughness but the highest lateral disorder of the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗
DOE OSTI2025

Automated segmentation of soft X-ray tomography: Native cellular structure with submicron resolution at high-throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at suboptical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based autosegmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline uses manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the three-dimensional (3D) whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single-cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Chen, Jianhua [Lawrence Berkeley National Laborato↗
DOE OSTI2023

Different higher order kinematics between star-forming and quiescent galaxies based on the SAMI, MAGPI, and LEGA-C surveys

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star \gt 10^{10.5} \, \mathrm{M_\odot }$) at redshifts z = 0.05, 0.3, and 0.8 from the SAMI, MAGPI, and LEGA-C surveys, to measure the excess kurtosis h 4 of the stellar velocity distribution, the latter parametrized as a Gauss–Hermite series. We find that at all redshifts where we have large enough samples, h 4 anticorrelates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anticorrelation with rotation-to-dispersion ratio, we also find a correlation between h 4 and M ⋆ , potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between h 4 and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, h 4 also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower h 4 compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Insights into distorted lamellar phases with small-angle scattering and machine learning

Lamellar phases are essential in various soft matter systems, with topological defects significantly influencing their mechanical properties. In this report, we present a machine-learning approach for quantitatively analyzing the structure and dynamics of distorted lamellar phases using scattering techniques. By leveraging the mathematical framework of Kolmogorov–Arnold networks, we demonstrate that the conformations of these distorted phases – expressed as superpositions of complex waves – can be reconstructed from small-angle scattering intensities. Through the contour analysis of wave field phase singularities, we obtain the statistics of the spatial distribution of topological defects. Furthermore, we establish that the temporal evolution of these defects can be derived from the time-dependent traveling wave field, informed by the dispersion relation of spectral components. This method opens new avenues for investigating the dynamics of distorted lamellar phases using various dynamic scattering techniques such as neutron spin echo and X-ray photon correlation spectroscopy. These findings enhance our microscopic understanding of how defects influence the physical properties of lamellar materials, with implications for both equilibrium and non-equilibrium states in general lamellar systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraining the Multiscale Structure of Geophysical Fields in Machine Learning: The Case of Precipitation

The use of deep-learning algorithms for estimating the value of geophysical variables from remotely sensed information is rapidly expanding. The typical objective function minimized in such algorithms is the mean square error (MSE), which is known to lead to smooth estimates with compressed dynamical range as compared to the true distribution of the variable of interest. Here, we introduce and evaluate alternative objective functions, focusing on the retrieval of precipitation rates from satellite passive microwave radiometric measurements using a deep convolutional neural network. For this testbed application, the results show that explicitly imposing the preservation of the statistical distribution and spatial wavelet power spectrum of the target variable allows to accurately reproduce extreme values and sharp gradients across multiple scales in this study.

58 GEOSCIENCES↗
DOE OSTI2025

Prospects for silvicultural enhancement of fire resistance in mesic westside forests of the Pacific Northwest

Increasing wildfire activity in mesic, temperate Pacific Northwest forests west of the Cascade Range crest has stimulated interest in understanding whether alternative forest management practices could reduce risk of stand-replacing fire. To explore how management can enhance fire resistance in these forests and assess tradeoffs among resistance enhancement, carbon sequestration and storage, and economic returns, we conducted 40-year simulations of stand development with BioSum, a framework for conducting landscape analysis with the Forest Vegetation Simulator (FVS), utilizing a statistically representative and spatially balanced sample of Forest Inventory and Analysis (FIA) plots. Simulation outcomes under business-as-usual silviculture were contrasted with fire-aware silviculture, and treatment optimization logic was developed and applied to represent landscape-scale outcomes under business-as-usual and fire-focused management scenarios. Simulation results indicate that fire-aware prescriptions and fire-focused management can meaningfully enhance stand- and landscape-scale fire resistance of westside forests under less than extreme fire weather, but at the cost of lower economic returns and reduced net carbon storage and sequestration over the 40-year analysis window. Shifting from business-as-usual regeneration harvests with short rotations to fire-aware, episodic selection harvest improved fire resistance the most, especially in young privately-owned forests, and with only modest tradeoffs in carbon and economic outcomes. While fire-aware treatments generally reduced net present value from forest operations over business-as-usual, most treatments still generated positive net present value and could be implemented without subsidy. Fire-aware prescriptions that removed and utilized non-merchantable harvest residues instead of burning them, via either pile or broadcast burning, partially mitigated carbon emissions associated with fire-aware treatments, with about the same improvement in fire resistance. Given the currently limited institutional and financial capacity to implement fire resistance enhancing treatments at scale, the insights from this analysis may aid managers seeking to elevate fire resistance to prioritize where and how to manage.

Science & Technology - Other Topics↗