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 235 records · Page 13

Architectures of Exoplanetary Systems: Towards a Multi-planet Model for Reproducing the Planet Radius Valley and Intra-system Size Similarity

Multi-planet systems provide crucial insights into the architectures and correlations within planetary systems, which in turn offer clues into their formation and evolution histories. Through the thousands of exoplanet candidates that it discovered, NASA's Kepler mission revealed that systems with multiple transiting planets are common. We describe a new model for inferring the exoplanet population observed by Kepler, which combines a semi-parametric description of planetary systems with a prescription for envelope mass-loss due to photoevaporation. This new model is a"hybrid" of the maximum AMD (angular momentum deficit) model in which the orbital excitations of the planets are drawn by distributing the system's critical AMD (He et al. 2020), and a photoevaporation model for sculpting the observed radius valley (Neil & Rogers 2020). This is an update to the SysSim suite of exoplanet population models derived from forward modeling the Kepler detection pipeline. The hybrid model is constrained by the observed statistics of multi-transiting systems, including additional constraints from the distribution of planet radii and their intra-system size patterns. This model is capable of producing a planet radius valley as well as a preference for planet size similarity. However, we find that the Kepler catalog exhibits features that are even more significant (i.e., it has a deeper observed radius valley and stronger size similarity) than the best-fit model. We use this model to infer the intrinsic planet radius valley as well as the distribution of primordial planet masses and core masses.

kepler↗

Dynamics and configurations of galaxy triplets

The purpose is to infer the probable dynamical states of galaxy triplets by the observed data on their configurations. Two methods are proposed for describing the distributions of the triplet configuration parameters characterizing a tendency to alignment and hierarchy: (1) obtaining a representative sample of configurations and determining its statistical parameters (moments and percentages); and (2) dividing the region of possible configurations of triple systems (Agekian and Anosova, 1967) into a set of segments and finding the probabilities for the configurations to find themselves in each of them. Both these methods allow representation of the data by numerical simulations as well as observations. The effect of projection was studied. It rather overestimates the alignment and hierarchy of the triple systems. Among the parameters of interest there are found some parameters that are least sensitive to projection effects. The samples consist of simulated galaxy triplets (with hidden mass) as well as of 46 probably physical triple galaxies (Karachentseva et al., 1979). The observed triples as well as numerical models show a tendency to alignment. The triple galaxies do not show any tendency to hierarchy (formation of the temporary binaries), but this tendency may be present for simulated triplets without significant dark matter. The significant hidden mass (of order ten times the total mass of a triplet) decreases the probability of forming a binary and so weakens the hierarchy. Small galaxy groups consisting of 3 to 7 members are probably the most prevalent types of galaxy aggregate (Gorbatsky, 1987). Galaxy triplets are the simplest groups, but dynamically nontrivial ones.

Anosova, Joanna P.↗

Patterns of seismic activity preceding large earthquakes

A mechanical model of seismic faults is employed to investigate the seismic activities that occur prior to major events. The block-and-spring model dynamically generates a statistical distribution of smaller slipping events that precede large events, and the results satisfy the Gutenberg-Richter law. The scaling behavior during a loading cycle suggests small but systematic variations in space and time with maximum activity acceleration near the future epicenter. Activity patterns inferred from data on seismicity in California demonstrate a regional aspect; increased activity in certain areas are found to precede major earthquake events. One example is given regarding the Loma Prieta earthquake of 1989 which is located near a fault section associated with increased activity levels.

Shaw, Bruce E.↗

Constraints of Melting, Sea-Level and the Paleoclimate from GRACE

To gauge the accuracy of the GRACE data, we have undertaken a study to compare deformations predicted by GRACE inferences of seasonal water loading to crustal position variations determined from GRACE data. Two manuscripts that resulted from this study are attached. We found a very high correlation between the GRACE and GPS determinations for South America [Duvis et al., 2004]. We also developed a statistical approach for choosing which Stokes coefficients to include. This approach proves to be somewhat more accurate than the traditional Gaussian filter [Duvis et al., 2005].

Davis, James L.↗

Using Historical Data to Automatically Identify Air-Traffic Control Behavior

This project seeks to develop statistical-based machine learning models to characterize the types of errors present when using current systems to predict future aircraft states. These models will be data-driven - based on large quantities of historical data. Once these models are developed, they will be used to infer situations in the historical data where an air-traffic controller intervened on an aircraft's route, even when there is no direct recording of this action.

trajectory generation↗

Cosmological Shocks Around Galaxy Clusters: A Coherent Investigation With DES, SPT & ACT

We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogues from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev–Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around 105 clusters with mass and redshift ranges 10 13.7 200m /M ⊙ <10 15.5 and 0.1 < z < 2, and the total sky coverage of the maps is ≈15000deg 2 ⁠. We find a clear pressure deficit at R/R 200m ≈ 1.1 in SZ profiles around both ACT and SPT clusters, estimated at 6σ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2–5 compared to predictions, and we discuss possible interpretations of this behaviour. An oriented stacking of clusters – where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogues – shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius.

galaxies↗

Pioneer-Venus radio occultation (ORO) data reduction: Profiles of 13 cm absorptivity

In order to characterize possible variations in the abundance and distribution of subcloud sulfuric acid vapor, 13 cm radio occultation signals from 23 orbits that occurred in late 1986 and 1987 (Season 10) and 7 orbits that occurred in 1979 (Season 1) were processed. The data were inverted via inverse Abel transform to produce 13 cm absorptivity profiles. Pressure and temperature profiles obtained with the Pioneer-Venus night probe and the northern probe were used along with the absorptivity profiles to infer upper limits for vertical profiles of the abundance of gaseous H2SO4. In addition to inverting the data, error bars were placed on the absorptivity profiles and H2SO4 abundance profiles using the standard propagation of errors. These error bars were developed by considering the effects of statistical errors only. The profiles show a distinct pattern with regard to latitude which is consistent with latitude variations observed in data obtained during the occultation seasons nos. 1 and 2. However, when compared with the earlier data, the recent occultation studies suggest that the amount of sulfuric acid vapor occurring at and below the main cloud layer may have decreased between early 1979 and late 1986.

Steffes, Paul G.↗

Robust Fuzzy Controllers Using FPGAs

Electro-mechanical device controllers typically come in one of three forms, proportional (P), Proportional Derivative (PD), and Proportional Integral Derivative (PID). Two methods of control are discussed in this paper; they are (1) the classical technique that requires an in-depth mathematical use of poles and zeros, and (2) the fuzzy logic (FL) technique that is similar to the way humans think and make decisions. FL controllers are used in multiple industries; examples include control engineering, computer vision, pattern recognition, statistics, and data analysis. Presented is a study on the development of a PD motor controller written in very high speed hardware description language (VHDL), and implemented in FL. Four distinct abstractions compose the FL controller, they are the fuzzifier, the rule-base, the fuzzy inference system (FIS), and the defuzzifier. FL is similar to, but different from, Boolean logic; where the output value may be equal to 0 or 1, but it could also be equal to any decimal value between them. This controller is unique because of its VHDL implementation, which uses integer mathematics. To compensate for VHDL's inability to synthesis floating point numbers, a scale factor equal to 10(sup (N/4) is utilized; where N is equal to data word size. The scaling factor shifts the decimal digits to the left of the decimal point for increased precision. PD controllers are ideal for use with servo motors, where position control is effective. This paper discusses control methods for motion-base platforms where a constant velocity equivalent to a spectral resolution of 0.25 cm(exp -1) is required; however, the control capability of this controller extends to various other platforms.

Monroe, Author Gene S., Jr.↗

Estimating the Length of the North Atlantic Basin Hurricane Season

For the interval 1945-2011, the length of the hurricane season in the North Atlantic basin averages about 130 +/- 42 days (the +/-1 standard deviation interval), having a range of 47 to 235 days. Runs-testing reveals that the annual length of season varies nonrandomly at the 5% level of significance. In particular, its trend, as described using 10-yr moving averages, generally has been upward since about 1979, increasing from about 113 to 157 days (in 2003). Based on annual values, one finds a highly statistically important inverse correlation at the 0.1% level of significance between the length of season and the occurrence of the first storm day of the season. For the 2012 hurricane season, based on the reported first storm day of May 19, 2012 (i.e., DOY = 140), the inferred preferential regression predicts that the length of the current season likely will be about 173 +/- 23 days, suggesting that it will end about November 8 +/- 23 days, with only about a 5% chance that it will end either before about September 23, 2012 or after about December 24, 2012.

Wilson, Robert M.↗

Bayesian Estimation of Earth’s Undiscovered Mineralogical Diversity Using Noninformative Priors

Recently, statistical distributions have been explored to provide estimates of the mineralogical diversity of Earth, and Earth-like planets. In this paper, a Bayesian approach is introduced to estimate Earth’s undiscovered mineralogical diversity. Samples are generated from a posterior distribution of the model parameters using Markov chain Monte Carlo simulations such that estimates and inference are directly obtained. It was previously shown that the mineral species frequency distribution conforms to a generalized inverse Gauss–Poisson (GIGP) large number of rare events model. Even though the model fit was good, the population size estimate obtained by using this model was found to be unreasonably low by mineralogists. In this paper, several zero-truncated, mixed Poisson distributions are fitted and compared, where the Poisson-lognormal distribution is found to provide the best fit. Subsequently, the population size estimates obtained by Bayesian methods are compared to the empirical Bayes estimates. Species accumulation curves are constructed and employed to estimate the population size as a function of sampling size. Finally, the relative abundances, and hence the occurrence probabilities of species in a random sample, are calculated numerically for all mineral species in Earth’s crust using the Poisson-lognormal distribution. These calculations are connected and compared to the calculations obtained in a previous paper using the GIGP model for which mineralogical criteria of an Earth-like planet were given.

Bayesian statistics↗

Earth rotation as a proxy for interannual variability in atmospheric circulation, 1860-present

Modern atmospheric and geodetic data sets have demonstrated that changes in the axial component of the atmosphere's angular momentum and in the rotation rate of the solid earth are closely coupled on time scales of up to several years. The feasibility of using a historical record of the earth's rotation as a proxy for year-to-year changes in the zonal wind field over the globe is examined. The bulk of the earth rotation series acquired for this purpose is based on telescopic observations of the occultation of stars by the moon; semiannual values of changes in the length of day derived from these observations have acceptably small errors from about 1860 onwards. These values are filtered to remove decade-scale fluctuations, which are driven primarily by nonatmospheric processes, and the resulting proxy series is examined to see if it contains a signal associated with one of the major modes of interannual variability in the atmosphere, namely that due to the El Nino/Southern Oscillation (ENSO). According to tests of statistical significance, such a signal is present in the historical earth rotation series, in that the day is typically longer during the year following an ENSO oceanic warm event than otherwise. Therefore other signals of interannual variability in the proxy series are considered. In particular, it is inferred that noteworthy trends in atmospheric interannual variability have occurred over the last century; for example, the decade of the 1920s was marked by much larger year-to-year changes in the zonal circulation over the globe than that of the 1940s. Based on modern atmospheric data, it is suggested that most of these circulation changes have resulted from anomalies in the region between 30 deg N and 30 deg S.

Salstein, David A.↗

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

Digital Lunar Exploration Sites (DLES) Terrain Crafting

Humans will soon be returning to the surface of the Moon with NASA’s Artemis program. The Artemis program is an international collaboration that will consist of a complex series of space systems and missions to explore the lunar surface and pave the way for the future exploration of Mars. NASA and its partners rely heavily on simulation for lighting and navigation studies as well as training astronauts, flight controllers, and mission support staff. The NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch (ER7) in the Engineering Directorate at NASA’s Johnson Space Center has built up many simulation products to support this effort, one of which is the Digital Lunar Exploration Sites (DLES). DLES is a collection of products used to simulate and render the lunar surface in a digital environment. We discussed and presented an overview of the DLES products at the 2022 IEEE Aerospace Conference in Big Sky, MT with a paper titled "Digital Lunar Exploration Sites". This “DLES Terrain Crafting” paper will expand on the information previously provided in “DLES” paper and dive deeper into the details of the terrain crafting process and the toolsets used to support this task. The best digital data currently available of the lunar surface is provided by the Lunar Reconnaissance Orbiter (LRO). Its Lunar Orbiter Laser Altimeter (LOLA) achieves an impressive resolution of 5m per pixel at the Lunar South Pole (LSP) and can generate datasets covering a large continuous region near the LSP. There are a few additional methods, such as Shape from Shading which can infer higher resolution data (up to 1m per pixel) from the LRO Narrow Angle Camera (NAC) images. However, surface-based simulations require higher-resolution data, and this paper will discuss the process of enhancing the terrain to meet that need. The process begins with capturing statistical data of craters in the regions of interest using images provided by the LRO NAC. This data is then used to scatter artificial features which are not captured in the truth data, resulting in an enhanced DEM with a much higher resolution of 20cm per pixel. Many tools were built up to assist in the creation of these artificial Digital Elevation Models (DEM), which this paper will discuss in detail. DEMs themselves are a very powerful representation of a planetary surface, and many operations and tools can utilize the data they contain. This paper includes a description of the rendering of the lunar surface in a graphics engine, generation of contact patches to simulate tire to ground interaction, and ray tracing utilities to model Line of Sight (LOS) interactions with the terrain. This paper will also explore some new tool sets currently under development which aim to utilize Machine Learning (ML) to assist in the identification of craters from LRO NAC imagery. While this is not a novel idea, the NExSyS team is developing a unique approach which may result in more robust identification of crater characteristics.

Artemis↗

LANDSAT menhaden and thread herring resources investigation

The author has identified the following significant results. The relationship between the distribution of menhaden and selected oceanographic parameters (water color, turbidity, and possibly chlorophyll concentrations) was established. Similar relationships for thread herring were not established nor were relationships relating to the abundance of either species. Use of aircraft and LANDSAT remote sensing instruments to measure or infer a set of basic oceanographic parameters was evaluated. Parameters which could be accurately inferred included surface water temperature, salinity, and color. Water turbidity (Secchi disk) was evaluated as marginally inferrable from the LANDSAT MSS data and chlorophyll-a concentrations as less than marginal. These evaluations considered the parameters only as experienced in the two test areas using available sensors and statistical techniques.

Kemmerer, A. J.↗

Soil moisture from temperature measurements at the Earth's surface, update

Soil moisture budgets at the Earth's surface were investigated based on soil and atmospheric temperature variations. A number of data sets were plotted and statistically analyzed in order to accentuate the existence and the characteristics of mesoscale soil temperature extrema variations and their relations to other parameters. The correlations between diurnal temperature extrema for air and soil in drought and non-drought periods appear to follow different characteristic patterns, allowing an inference of soil moisture content from temperature data. The recovery of temperature extrema after a precipitation event also follows a characteristic power curve rise between two limiting values which is an indicator of evaporation rates. If these indicators are applied universally to regional temperature data, soil moisture content or drought conditions can be inferred directly from temperature measurements.

Welker, J. E.↗

Vertical Scales of Turbulence at the Mt. Wilson Observatory

The vertical scales of turbulence at the Mt. Wilson observatory are inferred from data from the University of California at Berkeley Infrared Spatial Interferometer (ISI), by modeling path length fluctuations observed in the interferometric paths to celestial objects and those in instrumental ground- based paths. The correlations between the stellar and ground-based path length fluctuations and the temporal statistic of those fluctuations are modeled on various time scales to constrain the vertical scales.

Mt. Wilson Observatory↗

Vertical Scales of Turbulence at the Mount Wilson Observatory

The vertical scales of turbulence at the Mount Wilson Observatory are inferred from data from the University of California at Berkeley Infrared Spatial Interferometer (ISI), by modeling path length fluctuations observed in the interferometric paths to celestial objects and those in instrumental ground-based paths. The correlations between the stellar and ground-based path length fluctuations and the temporal statistics of those fluctuations are modeled on various timescales to constrain the vertical scales. A Kolmogorov-Taylor turbulence model with a finite outer scale was used to simulate ISI data. The simulation also included the white instrumental noise of the interferometer, aperture-filtering effects, and the data analysis algorithms. The simulations suggest that the path delay fluctuations observed in the 1992-1993 ISI data are largely consistent with being generated by refractivity fluctuations at two characteristic vertical scales: one extending to a height of 45 m above the ground, with a wind speed of about 1 m/ s, and another at a much higher altitude, with a wind speed of about 10 m/ s. The height of the lower layer is of the order of the dimensions of trees and other structures near the interferometer, which suggests that these objects, including elements of the interferometer, may play a role in generating the lower layer of turbulence. The modeling indicates that the high- attitude component contributes primarily to short-period (less than 10 s) fluctuations, while the lower component dominates the long-period (up to a few minutes) fluctuations. The lower component turbulent height, along with outer scales of the order of 10 m, suggest that the baseline dependence of long-term interferometric, atmospheric fluctuations should weaken for baselines greater than a few tens of meters. Simulations further show that there is the potential for improving the seeing or astrometric accuracy by about 30%-50% on average, if the path length fluctuations in the lower component are directly calibrated. Statistical and systematic effects induce an error of about 15 m in the estimate of the lower component turbulent altitude.

Treuhaft, Robert N.↗

Galaxy Cluster Gas Mass Fractions From Sunyaev-Zel'dovich Effect Measurement: Constraints on Omega_M

Using sensitive centimeter-wave receivers mounted on the Owens Valley Radio Observatory and Berkeley-Illinois-Maryland-Association millimeter arrays, we have obtained interferometric measurements of the Sunyaev-Zei'dovich (SZ) effect toward massive galaxy clusters. We use the SZ data to determine the pressure distribution of the cluster gas and, in combination with published X-ray temperatures, to infer the gas mass and total gravitational mass of 18 clusters. The gas mass fraction, f~g, is calculated for each cluster, and is extrapolated to the fiducial radius r_{500} using the results of numerical simulations. The mean f_g within r_{500} is 0.081 + 0.009 - 0.011/(h_{100} (statistical uncertainty at 68% confidence level, assuming OmegaM=0.3, OmegaL=0.7). We discuss possible sources of systematic errors in the mean f 9 measurement. We derive an upper limit for OmegaM from this sample under the assumption that the mass composition of clusters within r_{500} reflects the universal mass composition: Omega_M h <Omega_B/f_g.The gas mass f on cosmology through the angular diameter distance and the r_{500} correction factors. For a flat universe (Omegal, = 1 - OmegaM) and h=0.7, we find the measured gas mass fractions are consistent with OmegaM less than 0.40, at 68% confidence. Including estimates of the baryons contained in galaxies and the baryons which failed to become bound during the cluster formation process, we find OmegaM\approximately 0.25.

Grego, Laura↗