Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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

Advancing Neutrino Simulation Modeling with MARLEY: Insights from the NNSA-MSIIP Internship

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. However, there is no experimental data to validate the current event generator MARLEY (Model of Argon Reaction Low Energy Yields) that can model supernova neutrinos. To validate the model, I developed a new muon capture feature within the MARLEY simulation framework by coding key functions in C++, Python, and ROOT. I generated and analyzed one million simulated muon capture events and comparing the results to experimental data. To improve the model, I worked on optimizing the model’s parameters to improve its precision using statistical methods.

Wong, Baker [Fermilab]↗
DOE OSTI2019

Pore Formation in the Laminar Structure of Wear-Resistant Modified Powder Coatings

The open porosity of wear-resistant modified coatings is investigated, and their macrostructure is analyzed. A statistical method is proposed for describing the formation of porosity in the laminar structure of a powder coating. The matrix of transition probabilities is determined for a uniform Markov chain. The porosity of a coating on a smooth substrate is calculated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing the Contribution of ECa and NDVI in the Delineation of Management Zones in a Vineyard

Precision fertilization implies the need to identify the variability of soil fertility, which is costly and time-consuming. Remotely measured data can be a solution. Using this strategy, a study was conducted, in a vineyard, to delineate different management zones using two indicators: apparent soil electrical conductivity (ECa) and normalized difference vegetation index (NDVI). To understand the contribution of each indicator, three scenarios were used for zone definition: (1) using only NDVI, (2) only ECa, or (3) using a combination of the two. Then the differences in soil fertility between these zones were assessed using simple statistical methods. The results indicate that the most beneficial strategy is the combined use of the two indicators, as it allowed the definition of three distinct zones regarding important soil variables and crop nutrients, such as soil total nitrogen, Mg2+ cation, exchange acidity, and effective cation exchange capacity, and some relevant cation ratios. This strategy also allowed the identification of an ionic unbalance in the soil chemistry, due to an excess of Mg2+, that was harming crop health, as reported by NDVI. This also impacted ECa and NDVI relationship, which was negative in this study. Overall, the results demonstrate the advantages of using remotely sensed data, mainly more than one type of sensing data, and suggest a high potential for differential crop fertilization and soil management in the study area.

Esteves, Catarina (ORCID:0000000327205906)↗
DOE OSTI2021

Spatial Interpolation of Air Pollutant and Meteorological Variables in Central Amazonia

The Amazon Rainforest is highlighted by the global community both for its extensive vegetation cover that constantly suffers the effects of anthropic action and for its substantial biodiversity. This dataset presents data of meteorological variables from the Amazon Rainforest region with a spatial resolution of 0.001° in latitude and longitude, resulting from an interpolation process. The original data were obtained from the GoAmazon 2014/5 project, in the Atmospheric Radiation Measurement (ARM) repository, and then processed through mathematical and statistical methods. The dataset presented here can be used in experiments in the field of Data Science, such as training models for predicting climate variables or modeling the distribution of species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Information Length Concept Applied to Plasma Turbulence

A methodology to study statistical properties of anomalous transport in fusion plasma is investigated. Three time traces generated by the full-f gyrokinetic code GKNET are analyzed for this purpose. The time traces consist of heat flux as a function of the radial position, which is studied in a novel manner using statistical methods. The simulation data exhibit transport processes with both medium and long correlation length along the radius. A typical example of a phenomenon with long correlation length is avalanches. In order to investigate the evolution of the turbulent state, two basic configurations are studied, one flux-driven and one gradient-driven with decaying turbulence. The information length concept in tandem with Boltzmann–Gibbs and Tsallis entropy is used in the investigation. It is found that the dynamical states in both flux-driven and gradient-driven cases are surprisingly similar, but the Tsallis entropy reveals differences between them. This indicates that the types of probability distribution function are nevertheless quite different since the higher moments are significantly different.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

THE THIRD FERMI GBM GAMMA-RAY BURST CATALOG: THE FIRST SIX YEARS

Since its launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has triggered and located on average approximately two γ -ray bursts (GRBs) every three days. Here, we present the third of a series of catalogs of GRBs detected by GBM, extending the second catalog by two more years through the middle of 2014 July. The resulting list includes 1405 triggers identified as GRBs. The intention of the GBM GRB catalog is to provide information to the community on the most important observables of the GBM-detected GRBs. For each GRB, the location and main characteristics of the prompt emission, the duration, peak flux, and fluence are derived. The latter two quantities are calculated for the 50–300 keV energy band where the maximum energy release of GRBs in the instrument reference system is observed, and also for a broader energy band from 10 to 1000 keV, exploiting the full energy range of GBM's low-energy [Nai[Tl)] detectors. Using statistical methods to assess clustering, we find that the hardness and duration of GRBs are better fit by a two-component model with short-hard and long-soft bursts than by a model with three components. Furthermore, information is provided on the settings and modifications of the triggering criteria and exceptional operational conditions during years five and six in the mission. This third catalog is an official product of the Fermi GBM science team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A NEW REFERENCE CHEMICAL COMPOSITION FOR TMC-1

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi and Kaifu and the values reported in the review from Agúndez and Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Sub-band Image Reconstruction Using Differential Chromatic Refraction

Refraction by the atmosphere causes the positions of sources to depend on the airmass through which an observation was taken. This shift is dependent on the underlying spectral energy of the source and the filter or bandpass through which it is observed. Wavelength-dependent refraction within a single passband is often referred to as differential chromatic refraction (DCR). With a new generation of astronomical surveys undertaking repeated observations of the same part of the sky over a range of different airmasses and parallactic angles, DCR should be a detectable and measurable astrometric signal. In this paper we introduce a novel procedure that takes this astrometric signal and uses it to infer the underlying spectral energy distribution of a source; we solve for multiple latent images at specific wavelengths via a generalized deconvolution procedure built on robust statistics. We demonstrate the utility of such an approach for estimating a partially deconvolved image, at higher spectral resolution than the input images, for surveys such as the Large Synoptic Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Photometric and Kinematic Study of the Open Clusters SAI 44 and SAI 45

We carry out a detailed photometric and kinematic study of the poorly studied sparse open clusters SAI 44 and SAI 45 using ground-based BVR {sub c} I {sub c} data supplemented by archival data from Gaia eDR3 and Pan-STARRS. The stellar memberships are determined using a statistical method based on Gaia eDR3 kinematic data, and we found 204 members in SAI 44 while only 74 members are identified in SAI 45. The average distances to SAI 44 and SAI 45 are calculated to be 3670 ± 184 and 1668 ± 47 pc. The logarithmic age of the clusters are determined to be 8.82 ± 0.10 and 9.07 ± 0.10 yr for SAI 44 and SAI 45, respectively. The color–magnitude diagram of SAI 45 hosts an extended main-sequence turnoff (eMSTO). The apparent age spread is found to be similar to the apparent age spread predicted on the basis of the age spread and cluster age relation predicted by rotation models. This indicates that eMSTO is a stellar evolution rather than star formation phenomenon in SAI 45. We conclude that eMSTO in SAI 45 is mainly caused by the different rotation rates of stars as the SYCLIST synthetic population with different rotation rates was able to reproduce the observed eMSTO, and stars in the red part of the eMSTO were preferentially concentrated in the inner region, which again hints at different rotations being the reason for the extension in the upper MS. This finding supports the theory attributing the origin of eMSTO to the different rotations of eMSTO stars. The mass function slopes are obtained as -2.24 ± 0.66 and -2.58 ± 3.20 in the mass rages 2.426–0.990 M{sub ⊙} and 2.167–1.202 M{sub ⊙} for SAI 44 and SAI 45, respectively. SAI 44 exhibits the signature of mass segregation while we found weak evidence of mass segregation in SAI 45 possibly due to tidal stripping. The dynamical relaxation times of these clusters indicate that both clusters are in a dynamically relaxed state. Using the AD-diagram method, the apex coordinates are found to be ([[Formula]], -[[Formula]]) for SAI 44 and (-[[Formula]], -[[Formula]]) for SAI 45. The average space velocity components of the clusters SAI 44 and SAI 45 are calculated in units of km s{sup -1} as (-15.14 ± 3.90, -19.43 ± 4.41, -20.85 ± 4.57) and (28.13 ± 5.30, -9.78 ± 3.13, -19.59 ± 4.43), respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

From Neutron Star Observables to the Equation of State. II. Bayesian Inference of Equation of State Pressures

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. Here, we present a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parameterized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, ρ{sub sat}, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ∼30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass–radius relation can lead to biases of nearly 1 km in radius, toward larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Determination of Dark Matter Halo Mass from Dynamics of Satellite Galaxies

We propose that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using nine dark matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo, but has an intrinsic dependence on the halo formation history. We build this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of nine tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim 40 \% $. The accuracy can be improved to within ~25% if 30 tracers are used. Yet, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim 20 \% $ that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase-space information of other kinematic tracers is used.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Quantitative Constraints on the Reionization History from the IGM Damping Wing Signature in Two Quasars at z > 7

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral universe, and enables measurements of the hydrogen neutral fraction x HI (z) at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Lyα. Here, we combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two quasars with the highest redshifts known, ULAS J1120+0641 and ULAS J1342+0928, and measure volume-averaged neutral fractions {x HI }(Z = 7.09) = 0.48 $+0.26\atop{-0.26}$and {x HI }(Z = 7.54) = 0.60 $+0.20\atop{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of 10 3 ≤ t q ≤ 10 8 yr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

MicroRNAs as diagnostic and prognostic biomarkers of age-related macular degeneration: advances and limitations

A main cause of vision loss in the elderly is age-related macular degeneration (AMD). Among the cellular, biochemical, and molecular changes linked to this disease, inflammation and angiogenesis appear as being crucial in AMD pathogenesis and progression. There are two forms of the disease: dry AMD, accounting for 80–90% of cases, and wet AMD. The disease usually begins as dry AMD associated with retinal pigment epithelium and photoreceptor degeneration, whereas wet AMD is associated with choroidal neovascularization resulting in severe vision impairment. The new vessels are largely malformed, leading to blood and fluid leakage within the disrupted tissue, which provokes inflammation and scar formation and results in retinal damage and detachment. MicroRNAs are dysregulated in AMD and may facilitate the early detection of the disease and monitoring disease progression. Two recent reviews of microRNAs in AMD had indicated weaknesses or limitations in four earlier investigations. Studies in the last three years have shown considerable progress in overcoming some of these concerns and identifying specific microRNAs as biomarkers for AMD. Further large-scale studies are warranted using appropriate statistical methods to take into account gender and age disparity in the study populations and confounding factors such as smoking status.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A Comprehensive Open-Source R Software For Statistical Metrology Calculations: From Uncertainty Evaluation To Risk Analysis

Whether calibrating equipment or inspecting products on the factory floor, metrology requires many complicated statistical calculations to achieve a full understanding and evaluation of measurement uncertainty and quality. In order to assist its workforce in performing these calculations in a consistent and rigorous way, the Primary Standards Lab at Sandia National Laboratories (SNL) has developed a free and open-source software package for computing various metrology calculations from uncertainty propagation to risk analysis. In addition to propagating uncertainty through a measurement model using the well-known Guide to Expression of Uncertainty in Measurement or Monte Carlo approaches, evaluating the individual Type A and Type B uncertainty components that go into the measurement model often requires other statistical methods such as analysis of variance or determining uncertainty in a fitted curve. Once the uncertainty in a measurement has been calculated, it is usually evaluated from a risk perspective to ensure the measurement is suitable for making a particular conformance decision. Finally, SNL’s software can perform all these calculations in a single application via an easy-to-use graphical interface, where the different functions are integrated so the results of one calculation can be used as inputs to another calculation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Describing rainfall in northern Australia using multiple climate indices

Savanna landscapes are globally extensive and highly sensitive to climate change, yet the physical processes and climate phenomena which affect them remain poorly understood and therefore poorly represented in climate models. Both human populations and natural ecosystems are highly susceptible to precipitation variation in these regions due to the effects on water and food availability and atmosphere–biosphere energy fluxes. Here we quantify the relationship between climate phenomena and historical rainfall variability in Australian savannas and, in particular, how these relationships changed across a strong rainfall gradient, namely the North Australian Tropical Transect (NATT). Climate phenomena were described by 16 relevant climate indices and correlated against precipitation from 1900 to 2010 to determine the relative importance of each climate index on seasonal, annual and decadal timescales. Precipitation trends, climate index trends and wet season characteristics have also been investigated using linear statistical methods. In general, climate index–rainfall correlations were stronger in the north of the NATT where annual rainfall variability was lower and a high proportion of rainfall fell during the wet season. This is consistent with a decreased influence of the Indian–Australian monsoon from the north to the south. Seasonal variation was most strongly correlated with the Australian Monsoon Index, whereas yearly variability was related to a greater number of climate indices, predominately the Tasman Sea and Indonesian sea surface temperature indices (both of which experienced a linear increase over the duration of the study) and the El Niño–Southern Oscillation indices. These findings highlight the importance of understanding the climatic processes driving variability and, subsequently, the importance of understanding the relationships between rainfall and climatic phenomena in the Northern Territory in order to project future rainfall patterns in the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.

54 ENVIRONMENTAL SCIENCES↗