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 739 records · Page 41

Statistical chronometry of meteorites. I. A Test of 26 Al homogeneity and the Pb-Pb age of the solar system’s $t=0$

Here, we use rapidly cooled achondrites to test the assumption of 26 Al homogeneity in the solar nebula, by checking if there is a single value of t SS , the absolute “Pb-Pb” age of the Solar System’s t=0, that makes concordant their ages from the Al-Mg and Pb-Pb systems. We find that values t SS =4568.42 ± 0.24 Myr do make these ages concordant, and therefore the hypothesis of homogeneous 26 Al is not falsified. This age, defined to be when the solar nebula had ( 26 Al/ 27 Al) = 5.23 × 10 –5 , is significantly older than the ≈ 4567.3 Myr inferred from direct measurements of Pb-Pb ages in CAIs. Discrepancies between the Al-Mg and Pb-Pb chronometers in chondrules and CAIs have previously been interpreted as arising from heterogeneities in 26 Al, under the presumption that the Al-Mg and Pb-Pb systems in CAIs closed simultaneously. We examine this assumption and show that resetting is to be expected in CAIs. In particular, we quantitatively demonstrate that it is plausible that Pb-Pb ages of CAIs were reset at late times, without resetting the earlier Al-Mg ages, if they were transiently heated in the same manner as chondrules. We critically examine Pb-Pb isochrons, refining data and suggesting best practices for their calculation and reporting. We advocate reporting chronometry as times of formation after t=0 rather than absolute ages, as only the former is useful for astrophysical models of the solar nebula. We advocate averaging of multiple samples, rather than anchoring to individual meteorites, to improve precision.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗

Inference of the linear matter power spectrum at z = 0 using DESI DR1 Full-Shape data

Measurements of galaxy distributions at large cosmic distances capture clustering from the past. In this study, we use a cosmological model to translate these observations into the present-day galaxy distribution. Specifically, we reconstruct the 3D linear matter power spectrum at redshift z = 0 using Dark Energy Spectroscopic Instrument (DESI) Year 1 (DR1) galaxy clustering data and Cosmic Microwave Background (CMB) observations, assuming the ΛCDM model, and compare it to the result assuming the w 0 w a CDM model. Building on previous state-of-the-art methods, we apply Effective Field Theory (EFT) modelling of the galaxy power spectrum to account for small-scale effects in the 2-point statistics of galaxy data. Implementation of the EFT approach improves the modelling of the galaxy power spectrum, providing a more robust consistency test of the assumed cosmological model. By casting both CMB and galaxy clustering observations, spanning distinct redshift regimes, into k-space, we can identify discrepancies between the datasets of different redshifts, which would indicate potential inaccuracies in the assumed expansion history. While previous studies have shown consistency with ΛCDM, this work extends the analysis with higher-quality data to further test the expansion histories of both ΛCDM and w 0 w a CDM. Our findings show that both ΛCDM and w 0 w a CDM provide consistent fits to the linear matter power spectrum recovered from DESI DR1 data.

cosmological parameters from LSS↗

Inferring transmission heterogeneity using virus genealogies: estimation and targeted prevention

Spread of HIV typically involves uneven transmission patterns where some individuals spread to a large number of individuals while others to only a few or none. Such transmission heterogeneity can impact how fast and how much an epidemic spreads. Further, more efficient interventions may be achieved by taking such transmission heterogeneity into account. To address these issues, we developed two phylogenetic methods based on virus sequence data: 1) to generally detect if significant transmission heterogeneity is present, and 2) to pinpoint where in a phylogeny high-level spread is occurring. We derive inference procedures to estimate model parameters, including the amount of transmission heterogeneity, in a sampled epidemic. We show that it is possible to detect transmission heterogeneity under a wide range of simulated situations, including incomplete sampling, varying levels of heterogeneity, and including within-host genetic diversity. When evaluating real HIV-1 data from different epidemic scenarios, we found a lower level of transmission heterogeneity in slowly spreading situations and a higher level of heterogeneity in data that included a rapid outbreak, while R0 and Sackin’s index (overall tree shape statistic) were similar in the two scenarios, suggesting that our new method is able to detect transmission heterogeneity in real data. We then show by simulations that targeted prevention, where we pinpoint high-level spread using a coalescence measurement, is efficient when sequence data are collected in an ongoing surveillance system. Such phylogeny-guided prevention is efficient under both single-step contact tracing as well as iterative contact tracing as compared to random intervention.

59 BASIC BIOLOGICAL SCIENCES↗

A 5 per cent measurement of the Hubble–Lemaître constant from Type II supernovae

ABSTRACT The most stringent local measurement of the Hubble–Lemaître constant from Cepheid-calibrated Type Ia supernovae (SNe Ia) differs from the value inferred via the cosmic microwave background radiation (Planck+ΛCDM) by ∼5σ. This so-called Hubble tension has been confirmed by other independent methods, and thus does not appear to be a possible consequence of systematic errors. Here, we continue upon our prior work of using Type II supernovae to provide another, largely independent method to measure the Hubble–Lemaître constant. From 13 SNe II with geometric, Cepheid, or tip of the red giant branch (TRGB) host-galaxy distance measurements, we derive H$_0= 75.4^{+3.8}_{-3.7}$ km s−1 Mpc−1 (statistical errors only), consistent with the local measurement but in disagreement by ∼2.0σ with the Planck+ΛCDM value. Using only Cepheids (N = 7), we find H$_0 = 77.6^{+5.2}_{-4.8}$ km s−1 Mpc−1, while using only TRGB (N = 5), we derive H$_0 = 73.1^{+5.7}_{-5.3}$ km s−1 Mpc−1. Via 13 variants of our data set, we derive a systematic uncertainty estimate of 1.5 km s−1 Mpc−1. The median value derived from these variants differs by just 0.3 km s−1 Mpc−1 from that produced by our fiducial model. Because we only replace SNe Ia with SNe II – and we do not find statistically significant difference between the Cepheid and TRGB H0 measurements – our work reveals no indication that SNe Ia or Cepheids could be the sources of the ‘H0 tension.’ We caution, however, that our conclusions rest upon a modest calibrator sample; as this sample grows in the future, our results should be verified.

79 ASTRONOMY AND ASTROPHYSICS↗

Supernova-driven Turbulent Metal Mixing in High-redshift Galactic Disks: Metallicity Fluctuations in the Interstellar Medium and its Imprints on Metal-poor Stars in the Milky Way

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion among metal-poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this Letter is to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at parsec scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of α elements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal-poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low-star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical Assessment of Cetacean Stranding Events in Cape Cod Area

Cetacean (whales, dolphins and porpoises) mass strandings are a longstanding mystery in the field of marine biology that continue to be recorded in coastal environments around the world. For each of these events, anywhere from a few to several hundred otherwise healthy animals strand in onshore environments, often for no apparent reason. While the causes of these events remain unclear, anthropogenic and naturogenic mechanisms have been suggested. We present results of an inter-disciplinary study that draws expertise from space weather, marine mammal biology and ecology, and marine mammal stranding response. This study assessed 16 years of cetacean stranding events in the Cape Cod (Massachusetts, USA) area concurrently with a large dataset of meteorological, geophysical, biological, oceanographic and space weather data to produce inferences about possible causes for these unexplained events.

Dolphi↗

A causal data fusion method for the general exposure and outcome

Abstract With the advent of the big data era, the need to combine multiple individual data sets to draw causal effects arises naturally in many medical and biological applications. Especially each data set cannot measure enough confounders to infer the causal effect of an exposure on an outcome. In this article, we extend the method proposed by a previous study to causal data fusion of more than two data sets without external validation and to a more general (continuous or discrete) exposure and outcome. Theoretically, we obtain the condition for identifiability of exposure effects using multiple individual data sources for the continuous or discrete exposure and outcome. The simulation results show that our proposed causal data fusion method has unbiased causal effect estimate and higher precision than traditional regression, meta‐analysis and statistical matching methods. We further apply our method to study the causal effect of BMI on glucose level in individuals with diabetes by combining two data sets. Our method is essential for causal data fusion and provides important insights into the ongoing discourse on the empirical analysis of merging multiple individual data sources.

Li, Hongkai↗

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗

Small volcanoes in Tempe Terra, Mars: Their detailed morphometry and inferred geologic significance

The Tempe Terra province contains a variety of volcanic landforms that range in size from small vents (less than 10 km in diameter) to moderately sized volcanoes (150 km in diameter). The volcanoes are aligned along the dominant northerly and northeasterly trends of the faults in this region, and many of the volcanoes occur on grabens. Some workers have speculated on the nature of some of the volcanoes on the basis of their general morphology, shadow measurement height, lateral dimension, and geologic setting. As part of a larger study, we have obtained detailed photoclinometric profiles across five of the more conspicuous small volcanoes in the Tempe Terra region. For these data, we extracted for each volcano its flank width and edifice height and the width and depth of its summit crater. We statistically compared these dimensions for each volcano with a set of average dimensions for each type of terrestrial volcanic feature listed in Pike and Clow (1981). These comparisons indicate that the morphometries of the Martian volcanoes 1, 2, and 3 most closely match Earth's cratered basaltic lava shields, and the morphometries of volcanoes 4 and 5 are similar to those of terrestrial basaltic tuff rings.

Davis, P. A.↗

A review of some aspects of inferring the ozone profile by inversion of ultraviolet radiance measurements

The mathematical inversion of light-scattering observations to obtain the atmospheric ozone profile is discussed in terms of the filtering properties of the physical and mathematical processes for different spatial scales. Within this context, it is shown that the physical process of scattering acts as a low pass filter; which transfers large-scale profile information efficiently to radiance observations but very strongly attenuates small- or fine-scale profile information. To avoid domination of the mathematical inversion by the random error of radiance measurements, an equivalent spatial-scale filtering in the inversion procedure is essential. The available evidence suggests that mathematical inversion to obtain the low level ozone profile below 25 to 30 km is either inferior to or no better than the statistical estimation of the profile using total ozone as predictor. However, inversion profiles for high level ozone above 25 to 30 km appear to have moderately good accuracy.

Mateer, C. L.↗

A quantitative comparison of the fingerprint of twinned microstructures through surface and three-dimensional techniques

Assessing the fingerprint of a material’s microstructure is key for supporting materials design. With the emergence of a wide range of 3D characterization techniques, it is critical to understand the main differences in fingerprints reconstructed from 2D and 3D datasets. To this end, we introduce a graph-based microstructure reconstruction framework that enables structural comparisons of twin domain networks in high purity Ti using 3D and 2D electron backscatter diffraction. Insights into the structure of the twin networks are facilitated by combining statistical analysis of twin crystallography with visual and graphical analysis of the novel graph abstractions of the twins. We demonstrate that compared to 3D reconstructions, conventional 2D views of twinning miss key aspects of the microstructure including the high interconnectivity of domains into networks that span the full reconstruction volume. The reduced cross-grain and in-grain twin connectivity typically observed in 2D has notable implications on our understanding of how twinning mediates the plastic response of microstructures and how twin networks evolve. It is thus clear that 3D characterization is critical for accurately inferring both twin network morphologies as well as the key unit processes facilitating network formation.

36 MATERIALS SCIENCE↗

Robust Statistical Approach for Determination of Graphite Nitridation Using Bayesian Model Comparison

A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].

Engineering↗

Blinding multiprobe cosmological experiments

ABSTRACT The goal of blinding is to hide an experiment’s critical results – here the inferred cosmological parameters – until all decisions affecting its analysis have been finalized. This is especially important in the current era of precision cosmology, when the results of any new experiment are closely scrutinized for consistency or tension with previous results. In analyses that combine multiple observational probes, like the combination of galaxy clustering and weak lensing in the Dark Energy Survey (DES), it is challenging to blind the results while retaining the ability to check for (in)consistency between different parts of the data. We propose a simple new blinding transformation, which works by modifying the summary statistics that are input to parameter estimation, such as two-point correlation functions. The transformation shifts the measured statistics to new values that are consistent with (blindly) shifted cosmological parameters while preserving internal (in)consistency. We apply the blinding transformation to simulated data for the projected DES Year 3 galaxy clustering and weak lensing analysis, demonstrating that practical blinding is achieved without significant perturbation of internal-consistency checks, as measured here by degradation of the χ2 between the data and best-fitting model. Our blinding method’s performance is expected to improve as experiments evolve to higher precision and accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Calibration of the DSCOVR EPIC Visible and NIR Channels using Multiple LEO Radiometers

The first five years of operation of the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) at the Lagrange one point have produced results that uniquely complement the data from currently operating low orbit Earth-observing instruments. In this paper we describe an updated unified approach to EPIC calibration. In this approach, calibration coefficients and their trends were obtained by comparing EPIC observations to the measurements from polar orbiting radiometers. In this study L1B reflectances from Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua and Terra satellites, Multi-angle Imaging Spectroradiometer (MISR) onboard Terra and Visible Infrared Imaging Radiometer (VIIRS) onboard the Suomi National Polar-orbiting Partnership (Suomi NPP) spacecraft were used to infer calibration coefficients for four EPIC visible and near-infrared channels: 443 nm, 551 nm, 680 nm, and 780 nm. EPIC Version three measurements made between June 2015 and August 2020 were used for comparison. The calibration procedure identifies the most homogeneous low Earth orbit radiometer scenes matching scattering angles that are temporarily and spatially collocated with EPIC observations. These scenes are used to determine reflectance to count (R/C) ratios in spectrally analogous channels. Seasonal average R/C ratios were analyzed to obtain EPIC calibration gains and trends. The trends for the full dataset period are not statistically significant except in the 443 nm channel. No significant changes in calibration were found after the instrument’s exit from safe hold in March 2020. The R/C ratios were also used to determine the differences in EPIC gains resulting from separate calibrations: against MODIS Aqua or Terra, as well as against forward or aftward MISR cameras. Statistical tests indicate that the differences between the two datasets are not significant except in the 780 nm channels where Aqua-derived coefficients may be around 2% lower compared to Terra. The dependence of EPIC calibration gains on the instrument scattering angle and on DSCOVR-Earth distance were investigated. Lastly, model Low Earth Orbit (LEO) reflectances calculated to match the EPIC viewing geometry were employed to study how EPIC calibration coefficients depend on EPIC-LEO viewing geometry differences. The effect of LEO and EPIC angular mismatch on calibration was shown to be small.

deep space climate observatory/earth polychromatic↗

EFIT-Prime: Probabilistic and physics-constrained reduced-order neural network model for equilibrium reconstruction in DIII-D

We introduce EFIT-Prime, a novel machine learning surrogate model for EFIT (Equilibrium FIT) that integrates probabilistic and physics-informed methodologies to overcome typical limitations associated with deterministic and ad hoc neural network architectures. EFIT-Prime utilizes a neural architecture search-based deep ensemble for robust uncertainty quantification, providing scalable and efficient neural architectures that comprehensively quantify both data and model uncertainties. Physically informed by the Grad–Shafranov equation, EFIT-Prime applies a constraint on the current density J tor and a smoothness constraint on the first derivative of the poloidal flux, ensuring physically plausible solutions. Furthermore, the spatial location of the diagnostics is explicitly incorporated in the inputs to account for their spatial correlation. Extensive evaluations demonstrate EFIT-Prime's accuracy and robustness across diverse scenarios, most notably showing good generalization on negative-triangularity discharges that were excluded from training. Timing studies indicate an ensemble inference time of 15 ms for predicting a new equilibrium, offering the possibility of plasma control in real-time, if the model is optimized for speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constraints on the Cosmic Expansion History from GWTC–3

We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H(z), including its current value, the Hubble constant H0. Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H(z). The source mass distribution displays a peak around 34 M ⊙ , followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H(z) measurement, yielding H 0 = $68$$^{+12}_{–8}$km s –1 Mpc –1 (68% credible interval) when combined with the H 0 measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H 0 estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of H 0 = $68$$^{+8}_{–6}$km s –1 Mpc –1 with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H 0 studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H 0 ) is the well-localized event GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗

Adaptation of Mesoscale Weather Models to Local Forecasting

Methodologies have been developed for (1) configuring mesoscale numerical weather-prediction models for execution on high-performance computer workstations to make short-range weather forecasts for the vicinity of the Kennedy Space Center (KSC) and the Cape Canaveral Air Force Station (CCAFS) and (2) evaluating the performances of the models as configured. These methodologies have been implemented as part of a continuing effort to improve weather forecasting in support of operations of the U.S. space program. The models, methodologies, and results of the evaluations also have potential value for commercial users who could benefit from tailoring their operations and/or marketing strategies based on accurate predictions of local weather. More specifically, the purpose of developing the methodologies for configuring the models to run on computers at KSC and CCAFS is to provide accurate forecasts of winds, temperature, and such specific thunderstorm-related phenomena as lightning and precipitation. The purpose of developing the evaluation methodologies is to maximize the utility of the models by providing users with assessments of the capabilities and limitations of the models. The models used in this effort thus far include the Mesoscale Atmospheric Simulation System (MASS), the Regional Atmospheric Modeling System (RAMS), and the National Centers for Environmental Prediction Eta Model ( Eta for short). The configuration of the MASS and RAMS is designed to run the models at very high spatial resolution and incorporate local data to resolve fine-scale weather features. Model preprocessors were modified to incorporate surface, ship, buoy, and rawinsonde data as well as data from local wind towers, wind profilers, and conventional or Doppler radars. The overall evaluation of the MASS, Eta, and RAMS was designed to assess the utility of these mesoscale models for satisfying the weather-forecasting needs of the U.S. space program. The evaluation methodology includes objective and subjective verification methodologies. Objective (e.g., statistical) verification of point forecasts is a stringent measure of model performance, but when used alone, it is not usually sufficient for quantifying the value of the overall contribution of the model to the weather-forecasting process. This is especially true for mesoscale models with enhanced spatial and temporal resolution that may be capable of predicting meteorologically consistent, though not necessarily accurate, fine-scale weather phenomena. Therefore, subjective (phenomenological) evaluation, focusing on selected case studies and specific weather features, such as sea breezes and precipitation, has been performed to help quantify the added value that cannot be inferred solely from objective evaluation.

Manobianco, John T.↗