Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical inference”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Star-spots and magnetism: testing the activity paradigm in the Pleiades and M67

We measure star-spot filling fractions for 240 stars in the Pleiades and M67 open star clusters using APOGEE high-resolution H-band spectra. For this work, we developed a modified spectroscopic pipeline which solves for star-spot filling fraction and star-spot temperature contrast. We exclude binary stars, finding that the large majority of binaries in these clusters (80 per cent) can be identified from Gaia DR3 and APOGEE criteria – important for field star applications. Our data agree well with independent activity proxies, indicating that this technique recovers real star-spot signals. In the Pleiades, filling fractions saturate at a mean level of 0.248 ± 0.005 for active stars with a decline at slower rotation; we present fitting functions as a function of Rossby number. In M67, we recover low mean filling fractions of 0.030 ± 0.008 and 0.003 ± 0.002 for main sequence GK stars and evolved red giants, respectively, confirming that the technique does not produce spurious spot signals in inactive stars. Star-spots also modify the derived spectroscopic effective temperatures and convective overturn time-scales. Effective temperatures for active stars are offset from inactive ones by −109 ± 11 K, in agreement with the Pecaut & Mamajek empirical scale. Star-spot filling fractions at the level measured in active stars changes their inferred overturn time-scale, which biases the derived threshold for saturation. Finally, we identify a population of stars statistically discrepant from mean activity–Rossby relations and present evidence that these are genuine departures from a Rossby scaling. Our technique is applicable to the full APOGEE catalogue, with broad applications to stellar, galactic, and exoplanetary astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Inferring subhalo effective density slopes from strong lensing observations with neural likelihood-ratio estimation

ABSTRACT Strong gravitational lensing has emerged as a promising approach for probing dark matter (DM) models on sub-galactic scales. Recent work has proposed the subhalo effective density slope as a more reliable observable than the commonly used subhalo mass function. The subhalo effective density slope is a measurement independent of assumptions about the underlying density profile and can be inferred for individual subhaloes through traditional sampling methods. To go beyond individual subhalo measurements, we leverage recent advances in machine learning and introduce a neural likelihood-ratio estimator to infer an effective density slope for populations of subhaloes. We demonstrate that our method is capable of harnessing the statistical power of multiple subhaloes (within and across multiple images) to distinguish between characteristics of different subhalo populations. The computational efficiency warranted by the neural likelihood-ratio estimator over traditional sampling enables statistical studies of DM perturbers and is particularly useful as we expect an influx of strong lensing systems from upcoming surveys.

Astronomy & Astrophysics↗
DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Nondeterministic data base for computerized visual perception

A description is given of the knowledge representation data base in the perception subsystem of the Mars robot vehicle prototype. Two types of information are stored. The first is generic information that represents general rules that are conformed to by structures in the expected environments. The second kind of information is a specific description of a structure, i.e., the properties and relations of objects in the specific case being analyzed. The generic knowledge is represented so that it can be applied to extract and infer the description of specific structures. The generic model of the rules is substantially a Bayesian representation of the statistics of the environment, which means it is geared to representation of nondeterministic rules relating properties of, and relations between, objects. The description of a specific structure is also nondeterministic in the sense that all properties and relations may take a range of values with an associated probability distribution.

Yakimovsky, Y.↗
DOE OSTI2023

Inferring Viral Transmission Time from Phylogenies for Known Transmission Pairs

When the time of an HIV transmission event is unknown, methods to identify it from virus genetic data can reveal the circumstances that enable transmission. We developed a single-parameter Markov model to infer transmission time from an HIV phylogeny constructed of multiple virus sequences from people in a transmission pair. Our method finds the statistical support for transmission occurring in different possible time slices. We compared our time-slice model results to previously described methods: a tree-based logical transmission interval, a simple parsimony-like rules-based method, and a more complex coalescent model. Across simulations with multiple transmitted lineages, different transmission times relative to the source’s infection, and different sampling times relative to transmission, we found that overall our time-slice model provided accurate and narrower estimates of the time of transmission. We also identified situations when transmission time or direction was difficult to estimate by any method, particularly when transmission occurred long after the source was infected and when sampling occurred long after transmission. Applying our model to real HIV transmission pairs showed some agreement with facts known from the case investigations. We also found, however, that uncertainty on the inferred transmission time was driven more by uncertainty from time calibration of the phylogeny than from the model inference itself. Encouragingly, comparable performance of the Markov time-slice model and the coalescent model—which make use of different information within a tree—suggests that a new method remains to be described that will make full use of the topology and node times for improved transmission time inference.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ionizing radiation effects in SONOS-based neuromorphic inference accelerators

Here, we evaluate the sensitivity of neuromorphic inference accelerators based on Silicon-Oxide-Nitride-Oxide-Silicon (SONOS) charge trap memory arrays to total ionizing dose (TID) effects. Data retention statistics were collected for 16 Mbit of 40 nm SONOS digital memory exposed to ionizing radiation from a Co-60 source, showing good retention of the bits up to the maximum dose of 500 krad(Si). Using this data, we formulate a rate-equation-based model for the TID response of trapped charge carriers in the ONO stack, and predict the effect of TID on intermediate device states between ‘program’ and ‘erase’. This model is then used to simulate arrays of low-power, analog SONOS devices that store 8-bit neural network weights and support in situ matrix-vector multiplication. We evaluate the accuracy of the irradiated SONOS-based inference accelerator on two image recognition tasks – CIFAR-10 and the challenging ImageNet dataset – using state-of-the-art convolutional neural networks, such as ResNet-50. We find that across the datasets and neural networks evaluated, the accelerator tolerates a maximum TID between 10 krad(Si) and 100 krad(Si), with deeper networks being more susceptible to accuracy losses due to TID.

43 PARTICLE ACCELERATORS↗
NASA NTRS1984

Crop moisture estimation over the southern Great Plains with dual polarization 1.66 centimeter passive microwave data from Nimbus 7

Moisture content of snow-free, unfrozen soil is inferred using passive microwave brightness temperatures from the scanning multichannel microwave radiometer (SMMR) on Nimbus-7. Investigation is restricted to the two polarizations of the 1.66 cm wavelength sensor. Passive microwave estimates of soil moisture are of two basic categories; those based upon soil emissivity and those based upon the polarization of soil emission. The two methods are compared and contrasted through the investigation of 54 potential functions of polarized brightness temperatures and, in some cases, ground-based temperature measurements. Of these indices, three are selected for the estimated emissivity, the difference between polarized brightness temperatures, and the normalized polarization difference. Each of these indices is about equally effective for monitoring soil moisture. Using an antecedent precipitation index (API) as ground control data, temporal and spatial analyses show that emissivity data consistently give slightly better soil moisture estimates than depolarization data. The difference, however, is not statistically significant. It is concluded that polarization data alone can provide estimates of soil moisture in areas where the emissivity cannot be inferred due to nonavailability of surface temperature data.

Mcfarland, M. J.↗
DOE OSTI2026

CAHS: Context-Aware Homology Search

Protein homology search is foundational to bioinformatics: it supports annotation transfer, structure/function inference, and evolutionary analysis over rapidly expanding sequence repositories (e.g., UniProtKB). Profile hidden Markov models (pHMMs), as implemented in HMMER, remain the most widely trusted approach because they provide statistically calibrated E-values; however, their gap behavior is fixed once a profile is trained, despite biological evidence that insertion/deletion tolerance varies across flexible loops and intrinsically disordered regions. We present CAHS (Context-Aware Homology Search), a lightweight query-time adapter for pHMM search that incorporates learned and biologically motivated signals without changing HMMER's downstream search pipeline or its calibrated E-value reporting. Given a query sequence, CAHS computes per-residue representations from a protein language model and a disorder predictor, maps these to profile coordinates, and modulates only match-state transition rows (gap-open and gap-extension probabilities) while preserving Plan7 constraints. We comprehensively evaluate CAHS across six structurally diverse protein families and multi-domain architectures against a 570k-sequence target corpus. CAHS expands detection capability, retrieving thousands of additional remote homologs at relaxed thresholds by maintaining alignment quality through flexible regions. For multi-domain proteins, context-aware modulation resolves 94% of fragmented alignments. Crucially, CAHS preserves hit-set invariance at stringent operating points (E<10-10), demonstrating increased statistical confidence without inflating false positives. Furthermore, sharper statistical distinction between homologs and background noise during early filter stages yields up to a 3.87× acceleration in end-to-end wall-clock time on high-performance computing clusters. Overall, CAHS illustrates a practical AI-for-science design pattern: augmenting a trusted probabilistic model with query-specific learned signals to improve interpretable, reproducible inference in data-rich biology.

Bhattaram, Swethasree [Georgia Institute of Techno↗
DOE OSTI2020

A Probabilistic Model-Based Diagnostic Framework for Nuclear Engineering Systems

A fault diagnostic framework was investigated in this study for applications in thermal–hydraulic systems of nuclear power plants. The proposed framework consists of quantitative model-based diagnosis, statistical change detection and probabilistic reasoning. The use of physics-based diagnostic models provides high detection sensitivity and allows noise and measurement uncertainty to be incorporated robustly. Performance-related parametric models for each component are constructed based on first principles. Numerical model residuals are generated using the concept of analytical redundancy. Statistical change detection methods are employed to detect non-zero residuals in the presence of uncertainty. The diagnosis task is performed using Bayesian inference to detect and localize possible faults. Application to a single-phase heat exchanger for demonstration showed that the proposed probabilistic framework can provide improved results in comparison with traditional approaches while remaining less sensitive to false alarms in the presence of measurement and modeling uncertainty.

Bayesian network↗
NASA NTRS2013

Spatial Statistical Data Fusion (SSDF)

As remote sensing for scientific purposes has transitioned from an experimental technology to an operational one, the selection of instruments has become more coordinated, so that the scientific community can exploit complementary measurements. However, tech nological and scientific heterogeneity across devices means that the statistical characteristics of the data they collect are different. The challenge addressed here is how to combine heterogeneous remote sensing data sets in a way that yields optimal statistical estimates of the underlying geophysical field, and provides rigorous uncertainty measures for those estimates. Different remote sensing data sets may have different spatial resolutions, different measurement error biases and variances, and other disparate characteristics. A state-of-the-art spatial statistical model was used to relate the true, but not directly observed, geophysical field to noisy, spatial aggregates observed by remote sensing instruments. The spatial covariances of the true field and the covariances of the true field with the observations were modeled. The observations are spatial averages of the true field values, over pixels, with different measurement noise superimposed. A kriging framework is used to infer optimal (minimum mean squared error and unbiased) estimates of the true field at point locations from pixel-level, noisy observations. A key feature of the spatial statistical model is the spatial mixed effects model that underlies it. The approach models the spatial covariance function of the underlying field using linear combinations of basis functions of fixed size. Approaches based on kriging require the inversion of very large spatial covariance matrices, and this is usually done by making simplifying assumptions about spatial covariance structure that simply do not hold for geophysical variables. In contrast, this method does not require these assumptions, and is also computationally much faster. This method is fundamentally different than other approaches to data fusion for remote sensing data because it is inferential rather than merely descriptive. All approaches combine data in a way that minimizes some specified loss function. Most of these are more or less ad hoc criteria based on what looks good to the eye, or some criteria that relate only to the data at hand.

Braverman, Amy J.↗
NASA NTRS2009

Statistical and Scaling Properties of Remotely-Sensed Soil Moisture in Two Contrasting Domains in the North American Monsoon Region

Characterizing soil moisture (theta) variability is important for inferring high-resolution information from coarse estimates provided by remote sensors. In this study, we analyze the spatial variability and scale invariance of high-resolution theta estimates collected in two contrasting semiarid areas, Arizona (AZ) and Sonora (SON), during the Soil Moisture Experiment - North American Monsoon in 2004 (SMEX04- NAME). Results reveal that as the mean theta condition ( ) becomes drier, the spatial standard deviation becomes smaller in both domains. The coefficient of variation of theta decreases with in SON, but does not display a clear tendency with in AZ. We also found the presence of scale invariance and multifractality in the range of support scales from 51.2 km to 0.8 km for all soil moisture fields in the two regions. The multifractal properties of theta are clearly linked to in SON, while the relation is affected by more dispersion in AZ. We argue this is due to differences in the dynamic (rainfall) and static (vegetation) controls on theta in the two domains.

Mascaro, Giuseppe↗
DOE OSTI2025

Dark Energy Survey Year 3 results: $w$CDM cosmology from simulation-based inference with persistent homology on the sphere

We present cosmological constraints from Dark Energy Survey Year 3 (DES Y3) weak lensing data using persistent homology, a topological data analysis technique that tracks how features like clusters and voids evolve across density thresholds. For the first time, we apply spherical persistent homology to galaxy survey data through the algorithm TopoS2, which is optimized for curved-sky analyses and HEALPix compatibility. Employing a simulation-based inference framework with the Gower Street simulation suite, specifically designed to mimic DES Y3 data properties, we extract topological summary statistics from convergence maps across multiple smoothing scales and redshift bins. After neural network compression of these statistics, we estimate the likelihood function and validate our analysis against baryonic feedback effects, finding minimal biases (under $0.3σ$) in the $Ω_\mathrm{m}-S_8$ plane. Assuming the $w$CDM model, our combined Betti numbers and second moments analysis yields $S_8 = 0.821 \pm 0.018$ and $Ω_\mathrm{m} = 0.304\pm0.037$-constraints 70% tighter than those from cosmic shear two-point statistics in the same parameter plane. Our results demonstrate that topological methods provide a powerful and robust framework for extracting cosmological information, with our spherical methodology readily applicable to upcoming Stage IV wide-field galaxy surveys.

Prat, J. [Nordita; Royal Inst. Tech., Sodertalje; ↗
DOE OSTI2020

Broadband quadrature-squeezed vacuum and nonclassical photon number correlations from a nanophotonic device

We report demonstrations of both quadrature-squeezed vacuum and photon number difference squeezing generated in an integrated nanophotonic device. Squeezed light is generated via strongly driven spontaneous four-wave mixing below threshold in silicon nitride microring resonators. The generated light is characterized with both homodyne detection and direct measurements of photon statistics using photon number–resolving transition-edge sensors. We measure 1.0(1) decibels of broadband quadrature squeezing (~4 decibels inferred on-chip) and 1.5(3) decibels of photon number difference squeezing (~7 decibels inferred on-chip). Nearly single temporal mode operation is achieved, with measured raw unheralded second-order correlations g (2) as high as 1.95(1). Multiphoton events of over 10 photons are directly detected with rates exceeding any previous quantum optical demonstration using integrated nanophotonics. These results will have an enabling impact on scaling continuous variable quantum technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Tensile adhesion testing methodology for thermally sprayed coatings

The structure of thermally sprayed coatings consists of lamellae which are oriented parallel to the substrate surface. The lamellae separate and fracture by distinctive mechanisms which are reflected in the failure morphology, and these may be described as adhesive (between the coating and substrate), cohesive (within the coating), or mixed mode. There is a large variability in the failure stress for any nominally identical group of coatings. A lower bound for the fracture toughness of alumina coatings can be calculated as 0.2 MNm exp -3/2. The coating strength values may also be treated as belonging to the statistical distribution of the Weibull function. The Weibull modulus of the coating strength varied from 1.4 to 3.8. This analysis infers that the flaw size within coatings is highly variable and that the flaws are nonuniformly dispersed. The present work focuses on the question of whether tensile adhesion tests are an appropriate testing method for thermally sprayed materials.

Berndt, Christopher C.↗
NASA NTRS1993

Comparative study of x ray and microwave emissions during solar flares

The work supported by the grant consisted of two projects. The first project involved making detailed case studies of two flares using SMM data in conjunction with ground based observations. The first flare occurred at 1454 UT on June 20, 1989 and involved the eruption of a prominence near the limb. In the study we used data from many wavelength regimes including the radio, H-alpha, hard X-rays, and soft X-rays. We used a full gyrosynchrotron code to model the apparent presence of a 1.4 GHz source early in the flare that was in the form of a large coronal loop. The model results lead us to conclude that the initial acceleration occurs in small, dense loops which also produced the flare's hard X-ray emission. We also found evidence that a source at 1.4 GHz later in the event was due to second harmonic plasma emission. This source was adjacent to a leg of the prominence and comes from a dense column of material in the magnetic structure supporting the prominence. Finally, we investigated a source of microwaves and soft X-rays, occurring approximately 10 min after the hard X-ray peak, and calculate a lower limit for the density of the source. The second flare that was studied occurred at 2156 UT on June 20, 1989 and was observed with the VLA and the Owens Valley Radio Observatory (OVRO) Frequency Agile Array. We have developed a gyrosynchrotron model of the sources at flare peak using a new gyrosynchrotron approximation which is valid at very low harmonics of the gyrofrequency. We found that the accelerated particle densities of the sources decreased much more with radius from the source center than had been supposed in previous work, while the magnetic field varied less. We also used the available data to analyze a highly polarized source which appeared late in the flare. The second project involved compiling a statistical base for the relative timing of the hard X-ray peak, the turbulent and blue-shift velocities inferred from soft X-ray line emissions observed by SMM and the microwave peak as determined from ground-based observations. This timing was then used to aid the testing of newly developed global models for flares that incorporate the global magnetic topology as well as the electron dynamics that are responsible for the hard X-rays and microwaves.

Winglee, Robert M.↗
NASA NTRS2018

MMS Examination of FTEs at the Earth's Subsolar Magnetopause

Determining the magnetic field structure, electric currents, and plasma distributions within flux transfer event (FTE)-type flux ropes is critical to the understanding of their origin, evolution, and dynamics. Here the Magnetospheric Multiscale mission's high-resolution magnetic field and plasma measurements are used to identify FTEs in the vicinity of the subsolar magnetopause. The constant-alpha flux rope model is used to identify quasi-force free flux ropes and to infer the size, the core magnetic field strength, the magnetic flux content, and the spacecraft trajectories through these structures. Our statistical analysis determines a mean diameter of 1,700 ± 400 km (~30 ± 9 d(sub i)) and an average magnetic flux content of 100 ± 30 kWb for the quasi-force free FTEs at the Earth's subsolar magnetopause which are smaller than values reported by Cluster at high latitudes. These observed nonlinear size and magnetic flux content distributions of FTEs appear consistent with the plasmoid instability theory, which relies on the merging of neighboring, small-scale FTEs to generate larger structures. The ratio of the perpendicular to parallel components of current density, R(sub J), indicates that our FTEs are magnetically force-free, defined as R(sub J) < 1, in their core regions (<0.6 R(sub flux rope)). Plasma density is shown to be larger in smaller, newly formed FTEs and dropping with increasing FTE size. It is also shown that parallel ion velocity dominates inside FTEs with largest plasma density. Field-aligned flow facilitates the evacuation of plasma inside newly formed FTEs, while their core magnetic field strengthens with increasing FTE size.

Akhavan-Tafti, M.↗
NASA NTRS2012

Inferred Lunar Boulder Distributions at Decimeter Scales

Block size distributions of impact deposits on the Moon are diagnostic of the impact process and environmental effects, such as target lithology and weathering. Block size distributions are also important factors in trafficability, habitability, and possibly the identification of indigenous resources. Lunar block sizes have been investigated for many years for many purposes [e.g., 1-3]. An unresolved issue is the extent to which lunar block size distributions can be extrapolated to scales smaller than limits of resolution of direct measurement. This would seem to be a straightforward statistical application, but it is complicated by two issues. First, the cumulative size frequency distribution of observable boulders rolls over due to resolution limitations at the small end. Second, statistical regression provides the best fit only around the centroid of the data [4]. Confidence and prediction limits splay away from the best fit at the endpoints resulting in inferences in the boulder density at the CPR scale that can differ by many orders of magnitude [4]. These issues were originally investigated by Cintala and McBride [2] using Surveyor data. The objective of this study was to determine whether the measured block size distributions from Lunar Reconnaissance Orbiter Camera - Narrow Angle Camera (LROC-NAC) images (m-scale resolution) can be used to infer the block size distribution at length scales comparable to Mini-RF Circular Polarization Ratio (CPR) scales, nominally taken as 10 cm. This would set the stage for assessing correlations of inferred block size distributions with CPR returns [6].

Baloga, S. M.↗
DOE OSTI2026

Hierarchical Gaussian Random Field Sampling for Multilevel Markov Chain Monte Carlo: Coupling Stochastic Partial Differential Equation and the Karhunen–Loève Decomposition

This work introduces structure preserving hierarchical decompositions for sampling Gaussian random fields (GRFs) within the context of multilevel Bayesian inference in high-dimensional space. Existing scalable hierarchical sampling methods, such as those based on stochastic partial differential equations (SPDEs), often reduce the dimensionality of the sample space at the cost of accuracy of inference. Other approaches, such that those based on Karhunen-Loève (KL) expansions, offer sample space dimensionality reduction but sacrifice GRF representation accuracy and ergodicity of the Markov chain Monte Carlo (MCMC) sampler and are computationally expensive for high-dimensional problems. The proposed method integrates the dimensionality reduction capabilities of KL expansions with the scalability of SPDE-based sampling, thereby providing a robust, unified framework for high-dimensional uncertainty quantification (UQ) that is scalable and accurate, preserves ergodicity, and offers dimensionality reduction of the sample space. The hierarchy in our multilevel algorithm is derived from the geometric multigrid hierarchy. By constructing a hierarchical decomposition that maintains the covariance structure across the levels in the hierarchy, the approach enables efficient coarse-to-fine sampling while ensuring that all samples are drawn from the desired distribution. The effectiveness of the proposed method is demonstrated on a benchmark subsurface flow problem, demonstrating its effectiveness in improving computational efficiency and statistical accuracy. Furthermore, our proposed technique is more efficient and accurate and displays better convergence properties than existing methods for high-dimensional Bayesian inference problems.

Gaussian random fields↗