Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayesian clustering”

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.

26 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

The X-Ray Halo Scaling Relations of Supermassive Black Holes

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, and masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being −, followed by −. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the −) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anticorrelation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as subdominant over most of cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion, the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev–Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Supermassive black holes↗
NASA NTRS2006

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗
NASA NTRS2021

Computational Inference of Vibratory System with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

Inverse analysis of vibratory system is an important subject in fault identification, model updating, and robust design and control. It is challenging subject because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield results that are similar with respect to actual response measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic framework using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate it into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element inverse analysis. The validity of this framework is systematically demonstrated through case studies.

K Zhou↗
NASA NTRS1993

On evaluating the accuracy of SAR sea-ice classification using multifrequency polarimetric AIRSAR data

We investigate how multifrequency and polarimetric synthetic aperture radar (SAR) imagery enhances present capability to discriminate different ice conditions in single-frequency, single-polarization satellite SAR data. Frequencies considered are C- (lambda = 5.6cm), L- (lambda = 24cm) and P- (lambda = 68cm) band. Radar backscatter characteristics of six radiometrically and polarimetrically distinct ice types are selected from a cluster analysis of the multifrequency polarimetric SAR data and used to classify SAR images. Validation of these ice conditions is based on information provided by aerial photos, weather and ice surface measurements acquired at an ice camp, together with airborne passive microwave imagery, and visual analysis of the SAR data. The six identified sea-ice types are: (1) multiyear sea-ice; (2) compressed first year ice; (3) first year rubble and ridges; (4) first year rough ice; (5) first year smooth ice; and (6) thin ice. Open water is absent in all analyzed data. Classification of the SAR imagery into those six ice types is performed using a Bayesian Maximum A Posteriori classifier. Two complete scenes acquired at different dates in different locations are classified. The scenes were chosen such that they are representative of typical ice conditions in the Beaufort Sea in March 1988 and because ancillary information is available for validating the segmentation of various ice surface conditions.

Drinkwater, Mark R.↗
NASA NTRS2023

Using Regionalized Air Quality Model Performance and Bayesian Maximum Entropy Data Fusion to Map Global Surface Ozone Concentration and Associated Uncertainty

Estimates of ground-level ozone concentrations have been improved through data fusion of observations and atmospheric chemistry models. Our previous global ozone estimates for the Global Burden of Disease study corrected for bias uniformly across continents and then corrected near monitoring stations using the Bayesian Maximum Entropy (BME) framework for data fusion. Here, we use the Regionalized Air Quality Model Performance (RAMP) framework to correct model bias over a much larger spatial range than BME can, accounting for the spatial inhomogeneity of bias and nonlinearity as a function of modeled ozone. RAMP bias correction is applied to a composite of 9 global chemistry-climate models, based on the nearest set of monitors. These estimates are then fused with observations using BME, which matches observations at measurement stations, with the influence of observations declining with distance in space and time. We create global ozone maps for each year from 1990 to 2017 at fine spatial resolution. RAMP is shown to create unrealistic discontinuities due to the spatial clustering of ozone monitors, which we overcome by applying a weighting for RAMP based on the number of monitors nearby. Incorporating RAMP before BME has little effect on model performance near stations, but strongly increases R 2 by 0.15 at locations farther from stations, shown through a checkerboard cross-validation. Corrections to estimates differ based on location in space and time, confirming heterogeneity. We quantify the likelihood of exceeding selected ozone levels, finding that parts of the Middle East, India, and China are most likely to exceed 55 parts per billion (ppb) in 2017. About 96% of the global population was exposed to ozone levels above the World Health Organization guideline of 60 µg m −3 (30 ppb) in 2017. Our annual fine-resolution ozone estimates may be useful for several applications including epidemiology and assessments of impacts on health, agriculture, and ecosystems.

Ozone↗
NASA NTRS2010

Initialization and Restart in Stochastic Local Search: Computing a Most Probable Explanation in Bayesian Networks

For hard computational problems, stochastic local search has proven to be a competitive approach to finding optimal or approximately optimal problem solutions. Two key research questions for stochastic local search algorithms are: Which algorithms are effective for initialization? When should the search process be restarted? In the present work we investigate these research questions in the context of approximate computation of most probable explanations (MPEs) in Bayesian networks (BNs). We introduce a novel approach, based on the Viterbi algorithm, to explanation initialization in BNs. While the Viterbi algorithm works on sequences and trees, our approach works on BNs with arbitrary topologies. We also give a novel formalization of stochastic local search, with focus on initialization and restart, using probability theory and mixture models. Experimentally, we apply our methods to the problem of MPE computation, using a stochastic local search algorithm known as Stochastic Greedy Search. By carefully optimizing both initialization and restart, we reduce the MPE search time for application BNs by several orders of magnitude compared to using uniform at random initialization without restart. On several BNs from applications, the performance of Stochastic Greedy Search is competitive with clique tree clustering, a state-of-the-art exact algorithm used for MPE computation in BNs.

Mengshoel, Ole J.↗
NASA NTRS1995

Inductive Approaches to Improving Diagnosis and Design for Diagnosability

The first research area under this grant addresses the problem of classifying time series according to their morphological features in the time domain. A supervised learning system called CALCHAS, which induces a classification procedure for signatures from preclassified examples, was developed. For each of several signature classes, the system infers a model that captures the class's morphological features using Bayesian model induction and the minimum message length approach to assign priors. After induction, a time series (signature) is classified in one of the classes when there is enough evidence to support that decision. Time series with sufficiently novel features, belonging to classes not present in the training set, are recognized as such. A second area of research assumes two sources of information about a system: a model or domain theory that encodes aspects of the system under study and data from actual system operations over time. A model, when it exists, represents strong prior expectations about how a system will perform. Our work with a diagnostic model of the RCS (Reaction Control System) of the Space Shuttle motivated the development of SIG, a system which combines information from a model (or domain theory) and data. As it tracks RCS behavior, the model computes quantitative and qualitative values. Induction is then performed over the data represented by both the 'raw' features and the model-computed high-level features. Finally, work on clustering for operating mode discovery motivated some important extensions to the clustering strategy we had used. One modification appends an iterative optimization technique onto the clustering system; this optimization strategy appears to be novel in the clustering literature. A second modification improves the noise tolerance of the clustering system. In particular, we adapt resampling-based pruning strategies used by supervised learning systems to the task of simplifying hierarchical clusterings, thus making post-clustering analysis easier.

Fisher, Douglas H.↗
NASA NTRS2022

On the Effects of Cloud Water Content on Passive Microwave Snowfall Retrievals

The Bayesian passive microwave retrievals of snowfall often rely on mathematical matching of the observed vectors of brightness temperature with an a priori database of precipitation profiles and their corresponding brightness temperatures. Mathematical proximity does not necessarily lead to consistent retrievals due to limited information content of passive microwave observations. This paper defines imposter (genuine) vectors of brightness temperature as those that are mathematically close but physically inconsistent (consistent) and characterizes them through the Silhouette Coefficient (SC) analysis. The Neyman–Pearson (NP) hypothesis testing is used to separate the imposter and genuine brightness temperatures based on their associated values of cloud ice (IWP) and liquid water path (LWP), given by coincidences of CloudSat Profiling Radar (CPR) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI). The study determines thresholds for IWP and LWP that allow optimal identification of imposter brightness temperatures of non-snowing and snowing clouds, which can mislead the passive microwave retrieval algorithms to falsely detect or miss the snowfall events. It is demonstrated that emission signal of supercooled liquid water in snowing clouds can lead to improved passive microwave retrieval of snowfall and conditioning the retrievals to the cloud IWP and LWP can result in marginal correction of the snowfall detection probability; however, reduce the probability of false alarm by 6%–8% over sea ice and open oceans.

Snowfall↗