Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bayesian estimation”

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 289 records · Page 16

Markov Chain Monte Carlo Bayesian Learning for Neural Networks

Conventional training methods for neural networks involve starting al a random location in the solution space of the network weights, navigating an error hyper surface to reach a minimum, and sometime stochastic based techniques (e.g., genetic algorithms) to avoid entrapment in a local minimum. It is further typically necessary to preprocess the data (e.g., normalization) to keep the training algorithm on course. Conversely, Bayesian based learning is an epistemological approach concerned with formally updating the plausibility of competing candidate hypotheses thereby obtaining a posterior distribution for the network weights conditioned on the available data and a prior distribution. In this paper, we developed a powerful methodology for estimating the full residual uncertainty in network weights and therefore network predictions by using a modified Jeffery's prior combined with a Metropolis Markov Chain Monte Carlo method.

Goodrich, Michael S.↗

Data Processing for High-Contrast Imaging with the James Webb Space Telescope

The James Webb Space Telescope (JWST) will probe circumstellar environments at an unprecedented sensitivity.However, the performance of high-contrast imaging instruments is limited by the residual light from the star at close separations (<2-3”), where the incidence of exoplanets increases rapidly. There is currently no solution to get rid of the residual light down to the photon noise level at those separations, which may prevent some crucial discoveries. JWST’s debut is planned for October 2021 for a baseline mission lifetime of only five years. Thus,it is crucial to find a solution to this problem before its launch. We are further developing and implementing an innovative approach of post-processing built on a Bayesian framework that provides a more robust determination of faint astrophysical structures around a bright source. This approach uses a model of high-contrast imaging instrument that takes advantage of prior information, such as data from wavefront sensing (WFS) operations on JWST, to estimate simultaneously instrumental aberrations and the circumstellar environment. With this approach, our goal is to further improve the contrast gain over the contrast that can be achieved with JWST instruments, starting with NIRCam imaging and coronagraphic imaging. This work will pave the way for the future space-based high-contrast imaging instruments such as the Nancy Grace Roman Space Telescope Coronagraph Instrument (Roman CGI). This technique will be crucial to make the best use of the telemetry data that will be collected during the CGI operations.

Jewell, Jeffrey B.↗

Data Applicability of Heritage and New Hardware For Launch Vehicle Reliability Models

Bayesian reliability requires the development of a prior distribution to represent degree of belief about the value of a parameter (such as a component's failure rate) before system specific data become available from testing or operations. Generic failure data are often provided in reliability databases as point estimates (mean or median). A component's failure rate is considered a random variable where all possible values are represented by a probability distribution. The applicability of the generic data source is a significant source of uncertainty that affects the spread of the distribution. This presentation discusses heuristic guidelines for quantifying uncertainty due to generic data applicability when developing prior distributions mainly from reliability predictions.

Al Hassan, Mohammad↗

Continuous Habitable Zones: Pairing a GCM and Bayesian Framework to Predict Habitable Zone Evolution

In the near-future, new space telescopes like JWST, LUVOIR, and HabEx will begin attempting to explore the properties of atmospheres of potentially habitable planets. This will require a significant amount of time and resources for even a single planet, which makes it essential to prioritize observations by those most-likely to have detectable life. Here we present a statistical method to estimate the probabilities that specific exoplanets have been continuously in the habitable zone of their host stars for more than 2 billion years, the approximate time it took life on Earth to significantly increase the oxygen content of the atmosphere. We introduce the use of statistics of an ensemble of 3D planetary general circulation models to estimate these probabilities, replacing prior 1D model estimates.

habitable planets↗

Pairing a GCM and Bayesian Framework to Predict Habitable Zone Evolution

In the near-future, new space telescopes like JWST, LUVOIR, and HabEx will begin attempting to explore the properties of atmospheres of potentially habitable planets. This will require a significant amount of time and resources for even a single planet, which makes it essential to prioritize observations by those most-likely to have detectable life. Here we present a statistical method to estimate the probabilities that specific exoplanets have been continuously in the habitable zone of their host stars for more than 2 billion years, the approximate time it took life on Earth to significantly increase the oxygen content of the atmosphere. We introduce the use of statistics of an ensemble of 3D planetary general circulation models to estimate these probabilities, replacing prior 1D model estimates.

habitable planets↗

A new method for the detection of a periodic signal of unknown shape and period

A method was developed for the detection and measurement of a periodic signal with unknown characteristics in a data set where the existence of such signal was not known. The method detects a signal by using Bayesian probability theory to compare a constant model for the signal to members of a class of models with periodic structure. The method was applied to simulated data generated with both stepwise and sinusoidal light curves, demonstrating that such signals can be sensitively detected and the signal frequency and its shape can be accurately estimated.

Gregory, P. C.↗

OGLE-2018-BLG-1011Lb,c: Microlensing Planetary System with Two Giant Planets Orbiting a Low-mass Star

We report a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (three lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios ∼9.5 X 10(exp −3) and ∼15 X 10(eap −3) are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses 1.8 (sup +3.4, sub -1.1) M(J) and 2.8 (sup +5.8, sub -1.7) M(J) around a host with a mass 0.18 (sup +0.33,sub - 0.10) M(ʘ) and located at a distance 7.1(sup +1.1, sub -1.5) kpc. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet–host separations are a(⊥,2) = 1.8 (sup +2.1, sub -1.5)au (0.8 (sup +0.9,sub -0.6 au) and a(⊥,3) = 0.8 (sup +0.9, sub -0.6)au , where the values of a(⊥,2) inside and outside the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

Cheongho Han↗

Model-Based Method for Sensor Validation

Fault detection, diagnosis, and prognosis are essential tasks in the operation of autonomous spacecraft, instruments, and in situ platforms. One of NASA s key mission requirements is robust state estimation. Sensing, using a wide range of sensors and sensor fusion approaches, plays a central role in robust state estimation, and there is a need to diagnose sensor failure as well as component failure. Sensor validation can be considered to be part of the larger effort of improving reliability and safety. The standard methods for solving the sensor validation problem are based on probabilistic analysis of the system, from which the method based on Bayesian networks is most popular. Therefore, these methods can only predict the most probable faulty sensors, which are subject to the initial probabilities defined for the failures. The method developed in this work is based on a model-based approach and provides the faulty sensors (if any), which can be logically inferred from the model of the system and the sensor readings (observations). The method is also more suitable for the systems when it is hard, or even impossible, to find the probability functions of the system. The method starts by a new mathematical description of the problem and develops a very efficient and systematic algorithm for its solution. The method builds on the concepts of analytical redundant relations (ARRs).

Vatan, Farrokh↗

Reducing the Risk and Improving Mission Success for NASA's Human Mission to a Near-Earth Asteroid: How Many Robotic Surveyors?

NASA's recent attention and interest in sending a human mission to land on a Near-Earth asteroid raises the question of need for a robotic surveyor. This paper describes a Bayesian approach for comparing the productivity and cost-risk tradeoffs of sending (versus not sending) one or more robotic surveyor missions prior to a human mission to land on an asteroid. The probability of finding an asteroid suitable for landing was derived from an analysis of more than 1200 asteroids in order to define a quantitative estimate of suitability. The low cost of the surveyors relative to the human mission underlined the multi-surveyor strategy as relatively inexpensive insurance against the risks of encountering an unsuitable asteroid for landing on arrival by a human mission.

decision analysis↗

Health Management and Prognostics for Electric Aircraft Powertrain

W and c Any air borne vehicle needs incorporating safety as key parameter of measure, and inclusion of autonomy raises the critical need for safety under autonomous operations. Management of faults and component degradation is key as complexity in autonomous operations grow over the period of time. Therefore, in addition to basic operational requirements, an autonomous electric vehicle should be able to make accurate estimates of its current system health and take the correct decisions to complete its mission successfully. Real-time safety and state-awareness tools are therefore essential for the vehicle to be able to reach its destination in a safe and successful manner. The need for safety assurance and health management capabilities is particularly relevant for aircraft electric propulsion systems, which are relatively new and with limited historical to learn. They are critical systems requiring high power density along with reliability, resilience, efficient management of weight, and operational costs. A model- based fault diagnosis and prognostics approach of complex critical systems can successfully accomplish the safety and state awareness goal for such electric propulsion systems, enabling autonomous decision making capability for safe and efficient operation. To identify critical components in the system a Qualitative Bayesian approach using FMECA is implemented. This requires the assessment of some quantities representing the state of the electric unmanned aerial systems (e-UAS), as well as look-ahead forecasts of such states during the entire flight, presented in form of safety metrics (SM). In-service data and performance data gathered from degraded components sup- ports diagnostic and prognostic methods for these systems, but this data can be difficult to obtain as weight and packaging restrictions reduce redundancy and instrumentation on-board the vehicle. Therefore, an model-based framework should be capable or operating with limited data. In addition to data scarcity, the variability of such complex critical systems re- quires the model-based framework to reason in the presence of uncertainty, such as sensor noise, and modeling imperfections. Quantification of errors and uncertainties in the measured states and quantities is therefore a fundamental step for a precise estimation of such SMs; un-modeled uncertainty may result in erroneous state assessment and un- reliable predictions of future states of e-UAVs. Typical, centralized model-based schemes suffer from inherent disadvantages such as computational complexity, single point of failure, and scalability issues, and therefore may fail in such a complex scenario. This paper presents a methodology for developing a system level diagnostics and prognostics approach using a Qualitative Bayesian FMECA approach along with a formal uncertainty management framework for an e-UAS. In this work we demonstrate the efficacy of the framework to predict effects of sub-system level degradation on vehicle operation incorporating uncertainty management to predict future behavior under different operating conditions.

Kulkarni, Chetan↗

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Overview - Part 1

Knowledge of the amount of water in the soil is of great importance to many earth science disciplines. Soil moisture is a key variable in controlling the exchange of water and energy between the land surface and the atmosphere. Thus, soil moisture information is valuable in a wide range of applications including weather and climate, runoff potential and flood control, early warning of droughts, irrigation, crop yield forecasting, soil erosion, reservoir management, geotechnical engineering, and water quality. Despite the importance of soil moisture information, widespread and continuous measurements of soil moisture are not possible today. Although many earth surface conditions can be measured from satellites, we still cannot adequately measure soil moisture from space. Research in soil moisture remote sensing began in the mid 1970s shortly after the surge in satellite development. Recent advances in remote sensing have shown that soil moisture can be measured, at least qualitatively, by several methods. Quantitative measurements of moisture in the soil surface layer have been most successful using both passive and active microwave remote sensing, although complications arise from surface roughness and vegetation type and density. Early attempts to measure soil moisture from space-borne microwave instruments were hindered by what is now considered sub-optimal wavelengths (shorter than 5 cm) and the coarse spatial resolution of the measurements. L-band frequencies between 1 and 3 GHz (10-30 cm) have been deemed optimal for detection of soil moisture in the upper few centimeters of soil. The Electronically Steered Thinned Array Radiometer (ESTAR), an aircraft-based instrument operating a 1,4 GHz, has shown great promise for soil moisture determination. Initiatives are underway to develop a similar instrument for space. Existing space-borne synthetic aperture radars (SARS) operating at C- and L-band have also shown some potential to detect surface wetness. The advantage of radar is its much higher resolution than passive microwave systems, but it is currently hampered by surface roughness effects and the lack of a good algorithm based on a single frequency and single polarization. In addition, its repeat frequency is generally low (about 40 days). In the meantime, two new radiometers offer some hope for remote sensing of soil moisture from space. The Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), launched in November 1997, possesses a 10.65 GHz channel and the Advanced Microwave Scanning Radiometer (AMSR) on both the ADEOS-11 and Earth Observing System AM-1 platforms to be launched in 1999 possesses a 6.9 GHz channel. Aside from issues about interference from vegetation, the coarse resolution of these data will provide considerable challenges pertaining to their application. The resolution of TMI is about 45 km and that of AMSR is about 70 km. These resolutions are grossly inconsistent with the scale of soil moisture processes and the spatial variability of factors that control soil moisture. Scale disparities such as these are forcing us to rethink how we assimilate data of various scales in hydrologic models. Of particular interest is how to assimilate soil moisture data by reconciling the scale disparity between what we can expect from present and future remote sensing measurements of soil moisture and modeling soil moisture processes. It is because of this disparity between the resolution of space-based sensors and the scale of data needed for capturing the spatial variability of soil moisture and related properties that remote sensing of soil moisture has not met with more widespread success. Within a single footprint of current sensors at the wavelengths optimal for this application, in most cases there is enormous heterogeneity in soil moisture created by differences in landcover, soils and topography, as well as variability in antecedent precipitation. It is difficult to interpret the meaning of 'mean' soil moisture under such conditions and even more difficult to apply such a value. Because of the non-linear relationships between near-surface soil moisture and other variables of interest, such as surface energy fluxes and runoff, mean soil moisture has little applicability at such large scales. It is for these reasons that the use of remote sensing in conjunction with a hydrologic model appears to be of benefit in capturing the complete spatial and temporal structure of soil moisture. This paper is Part I of a four-part series describing a method for intermittently assimilating remotely-sensed soil moisture information to improve performance of a distributed land surface hydrology model. The method, summarized in section II, involves the following components, each of which is detailed in the indicated section of the paper or subsequent papers in this series: Forward radiative transfer model methods (section II and Part IV); Use of a Kalman filter to assimilate remotely-sensed soil moisture estimates with the model profile (section II and Part IV); Application of a soil hydrology model to capture the continuous evolution of the soil moisture profile within and below the root zone (section III); Statistical aggregation techniques (section IV and Part II); Disaggregation techniques using a neural network approach (section IV and Part III); and Maximum likelihood and Bayesian algorithms for inversely solving for the soil moisture profile in the upper few cm (Part IV).

Crosson, William L.↗

Interpolating Fields of Carbon Monoxide Data Using a Hybrid Statistical-Physical Model

Atmospheric Carbon Monoxide (CO) is a pollutant gas of which the US congress has mandated regular monitoring, and satellite sensors can be used to retrieve regional concentrations of CO over several vertical layers. However, CO at cloudy locations cannot be observed and have to be estimated from the observed data set, resulting in an interpolation problem. The current state-of-the-art solution is to combine prior information, computed by a deterministic physical model, with observations. However, the deterministic model may introduce uncertainties that do not derive from the data. While sharing certain features with the physical model, this paper presents a Bayesian hierarchical model for interpolating CO on a 3-dimensional spatial grid, across time. To our knowledge such a model has not been considered before. The model is applied to a hypothetical air-quality monitoring scenario, and is compared to existing interpolation methods. The results provide motivation for the use of the statistical model for regional to local applications.

Arellano, A. A.↗

A Bayesian approach to nonlinear inversion

Powerful methods are now available for solving linear parametric inverse problems. However, many inverse problems which arise in geohysics are nonlinear. Fortunately, it is possible to treat most of these with the air of linear perturbation theory and liner inversion. But a convenient method is needed for assessing the importance of nonlinearity in these quasi-linear problems. The present paper provides such a method. Matsu'ura and Jackson (1984) have presented a simple algorithm for evaluating the asymptotic covariance matrix fo estimation errors. In the present investigation, aspects of linear inversion are discussed, taking into account linear parametric inverse problems, nonuniqueness, prior information, confidence limits, conditional and marginal statistics, the relative importance of the prior and observational data, and standardized variables. Attention is also given to nonlinear inversion, and the application of the considered approaches to a number of examples.

Jackson, D. D.↗

A Method for Retrieving Ground Flash Fraction from Satellite Lightning Imager Data

A general theory for retrieving the fraction of ground flashes in N lightning observed by a satellite-based lightning imager is provided. An "exponential model" is applied as a physically reasonable constraint to describe the measured optical parameter distributions, and population statistics (i.e., mean, variance) are invoked to add additional constraints to the retrieval process. The retrieval itself is expressed in terms of a Bayesian inference, and the Maximum A Posteriori (MAP) solution is obtained. The approach is tested by performing simulated retrievals, and retrieval error statistics are provided. The ability to retrieve ground flash fraction has important benefits to the atmospheric chemistry community. For example, using the method to partition the existing satellite global lightning climatology into separate ground and cloud flash climatologies will improve estimates of lightning nitrogen oxides (NOx) production; this in turn will improve both regional air quality and global chemistry/climate model predictions.

Koshak, William J.↗

Bayesian retrieval of optically relevant properties from hyperspectral water-leaving reflectances

Current methods to retrieve optically relevant properties from ocean color observations do not explicitly make use of prior knowledge about property distributions. Here we implement a simplified Bayesian approach that takes into account prior probability distributions on two sets of five optically relevant parameters, and conduct a retrieval of these parameters using hyperspectral simulated water-leaving reflectances. We focus specifically on the ability of the model to distinguish between two optically similar phytoplankton taxa, diatoms and Noctiluca scintillans. The inversion retrieval gives most-likely concentrations and uncertainty estimates, and we find that the model is able to probabilistically predict the occurrence of Noctiluca scintillans blooms using these metrics. We discuss how this method can be expanded to include a priori covariances between different parameters, and show the effect of varying measurement uncertainty and spectral resolution on Noctiluca scintillans bloom predictions.

Zachary K Erickson↗

Non-Targeted Effects and the Dose Response for Heavy Ion Tumorigenesis

BACKGROUND: There is no human epidemiology data available to estimate the heavy ion cancer risks experienced by astronauts in space. Studies of tumor induction in mice are a necessary step to estimate risks to astronauts. Previous experimental data can be better utilized to model dose response for heavy ion tumorigenesis and plan future low dose studies. DOSE RESPONSE MODELS: The Harderian Gland data of Alpen et al.[1-3] was re-analyzed [4] using non-linear least square regression. The data set measured the induction of Harderian gland tumors in mice by high-energy protons, helium, neon, iron, niobium and lanthanum with LET s ranging from 0.4 to 950 keV/micron. We were able to strengthen the individual ion models by combining data for all ions into a model that relates both radiation dose and LET for the ion to tumor prevalence. We compared models based on Targeted Effects (TE) to one motivated by Non-targeted Effects (NTE) that included a bystander term that increased tumor induction at low doses non-linearly. When comparing fitted models to the experimental data, we considered the adjusted R2, the Akaike Information Criteria (AIC), and the Bayesian Information Criteria (BIC) to test for Goodness of fit.In the adjusted R2test, the model with the highest R2values provides a better fit to the available data. In the AIC and BIC tests, the model with the smaller values of the summary value provides the better fit. The non-linear NTE models fit the combined data better than the TE models that are linear at low doses. We evaluated the differences in the relative biological effectiveness (RBE) and found the NTE model provides a higher RBE at low dose compared to the TE model. POWER ANALYSIS: The final NTE model estimates were used to simulate example data to consider the design of new experiments to detect NTE at low dose for validation. Power and sample sizes were calculated for a variety of radiation qualities including some not considered in the Harderian Gland data set and with different background tumor incidences. We considered different experimental designs with varying number of doses and varying low doses dependant on the LET of the radiation. The optimal design to detect a NTE for an individual ion had 4 doses equally spaced below a maximal dose where bending due to cell sterilization was < 2%. For example at 100 keV/micron we would irradiate at 0.03 Gy, 0.065 Gy, 0.13 Gy, and 0.26 Gy and require 850 mice including a control dose for a sensitivity to detect NTE with 80% power. Sample sizes could be improved by combining ions similar to the methods used with the Harderian Gland data.

Chappelli, Lori J.↗

On the design of optimal input signals in system identification

The problem of designing optimal inputs in the identification of linear systems with unknown random parameters is considered using a Bayesian approach. The information matrix, which is positive definite for the class of systems analyzed, gives a measure of performance for the system inputs. The computation of the optimal closed-loop input mappings is shown to be a nontrivial exercise in adaptive control. Deterministic optimal inputs are shown to be easily computable. Numerical examples are given. A Kalman filter is used to estimate the parameters. A necessary condition for the Kalman filter not to diverge when applying linear feedback is also given.

Lopez-Toledo, A. A.↗

The New NASA Orbital Debris Engineering Model ORDEM 3.0

The NASA Orbital Debris Program Office (ODPO) has released its latest Orbital Debris Engineering Model, ORDEM 3.0. It supersedes ORDEM 2.0. This newer model encompasses the Earth satellite and debris flux environment from altitudes of low Earth orbit (LEO) through geosynchronous orbit (GEO). Debris sizes of 10 microns through 1 m in non-GEO and 10 cm through 1 m in GEO are modeled. The inclusive years are 2010 through 2035. The ORDEM model series has always been data driven. ORDEM 3.0 has the benefit of many more hours from existing data sources and from new sources that weren't available to past versions. Returned surfaces, ground tests, and remote sensors all contribute data. The returned surface and ground test data reveal material characteristics of small particles. Densities of fragmentation debris particles smaller than 10 cm are grouped in ORDEM 3.0 in terms of high-, medium-, and lowdensities, along with RORSAT sodium-potassium droplets. Supporting models have advanced significantly. The LEO-to-GEO ENvironment Debris model (LEGEND) includes an historical and a future projection component with yearly populations that include launched and maneuvered intacts, mission related debris (MRD), and explosion and collision fragments. LEGEND propagates objects with ephemerides and physical characteristics down to 1 mm in size. The full LEGEND yearly population acts as an a priori condition for a Bayesian statistical model. Specific, well defined populations are added like the Radar Ocean Reconnaissance Satellite (RORSAT) sodium-potassium (NaK) droplets, recent major accidental and deliberate collision fragments, and known anomalous debris event fragments. For microdebris of sizes 10 microns to 1 mm the ODPO uses an in-house Degradation/Ejecta model in which a MLE technique is used with returned surface data to estimate populations. This paper elaborates on the upgrades of this model over previous versions highlighting the material density splits and consequences of that to the penetration risk to spacecraft.

Krisko, P. H.↗