Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, 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 181 records · Page 10

Automatic Generation of Algorithms for the Statistical Analysis of Planetary Nebulae Images

Analyzing data sets collected in experiments or by observations is a Core scientific activity. Typically, experimentd and observational data are &aught with uncertainty, and the analysis is based on a statistical model of the conjectured underlying processes, The large data volumes collected by modern instruments make computer support indispensible for this. Consequently, scientists spend significant amounts of their time with the development and refinement of the data analysis programs. AutoBayes [GF+02, FS03] is a fully automatic synthesis system for generating statistical data analysis programs. Externally, it looks like a compiler: it takes an abstract problem specification and translates it into executable code. Its input is a concise description of a data analysis problem in the form of a statistical model as shown in Figure 1; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. Internally, however, it is quite different: AutoBayes derives a customized algorithm implementing the given model using a schema-based process, and then further refines and optimizes the algorithm into code. A schema is a parameterized code template with associated semantic constraints which define and restrict the template s applicability. The schema parameters are instantiated in a problem-specific way during synthesis as AutoBayes checks the constraints against the original model or, recursively, against emerging sub-problems. AutoBayes schema library contains problem decomposition operators (which are justified by theorems in a formal logic in the domain of Bayesian networks) as well as machine learning algorithms (e.g., EM, k-Means) and nu- meric optimization methods (e.g., Nelder-Mead simplex, conjugate gradient). AutoBayes augments this schema-based approach by symbolic computation to derive closed-form solutions whenever possible. This is a major advantage over other statistical data analysis systems which use numerical approximations even in cases where closed-form solutions exist. AutoBayes is implemented in Prolog and comprises approximately 75.000 lines of code. In this paper, we take one typical scientific data analysis problem-analyzing planetary nebulae images taken by the Hubble Space Telescope-and show how AutoBayes can be used to automate the implementation of the necessary anal- ysis programs. We initially follow the analysis described by Knuth and Hajian [KHO2] and use AutoBayes to derive code for the published models. We show the details of the code derivation process, including the symbolic computations and automatic integration of library procedures, and compare the results of the automatically generated and manually implemented code. We then go beyond the original analysis and use AutoBayes to derive code for a simple image segmentation procedure based on a mixture model which can be used to automate a manual preproceesing step. Finally, we combine the original approach with the simple segmentation which yields a more detailed analysis. This also demonstrates that AutoBayes makes it easy to combine different aspects of data analysis.

Fischer, Bernd↗

Components of interannual ozone change based on Nimbus 7 TOMS data

A multiple regression statistical model is applied to estimate the latitude and seasonal dependences of the solar cycle, quasi-biennial oscillation (QBO), and anthropogenic trend components of stratospheric total ozone change using 13.2 years of Nimbus 7 TOMS data. The characteristics of the linear trend component are in agreement with earlier studies. The QBO regression coefficient is significantly different from zero at high southern latitudes in the Austral spring supporting earlier evidence that the Antarctic ozone depletion is modulated by the QBO. The existence of a solar cycle component is indicated by empirical studies of model residuals and by the approximate agreement of the derived global mean solar coefficient amplitude with photochemical calculations. Initial estimates for the latitude dependence of the solar coefficient suggest higher amplitudes with increasing latitude, especially in the Southern Hemisphere in spring. The statistical model predicts a return to more rapid ozone depletions during the next 4 years as solar minimum is approached.

Hood, Lon L.↗

Uncertainty in Model Predictions of Vibrio Vulnificus Response to Climate Variability and Change: A Chesapeake Bay Case Study

The effect that climate change and variability will have on waterborne bacteria is a topic of increasing concern for coastal ecosystems, including the Chesapeake Bay. Surface water temperature trends in the Bay indicate a warming pattern of roughly 0.3-0.4 C per decade over the past 30 years. It is unclear what impact future warming will have on pathogens currently found in the Bay, including Vibrio spp. Using historical environmental data, combined with three different statistical models of Vibrio vulnificus probability, we explore the relationship between environmental change and predicted Vibrio vulnificus presence in the upper Chesapeake Bay. We find that the predicted response of V. vulnificus probability to high temperatures in the Bay differs systematically between models of differing structure. As existing publicly available datasets are inadequate to determine which model structure is most appropriate, the impact of climatic change on the probability of V. vulnificus presence in the Chesapeake Bay remains uncertain. This result points to the challenge of characterizing climate sensitivity of ecological systems in which data are sparse and only statistical models of ecological sensitivity exist.

model predictions↗

Spatial Statistical Data Fusion (SSDF)

As remote sensing for scientific purposes has transitioned from an experimental technology to an operational one, the selection of instruments has become more coordinated, so that the scientific community can exploit complementary measurements. However, tech nological and scientific heterogeneity across devices means that the statistical characteristics of the data they collect are different. The challenge addressed here is how to combine heterogeneous remote sensing data sets in a way that yields optimal statistical estimates of the underlying geophysical field, and provides rigorous uncertainty measures for those estimates. Different remote sensing data sets may have different spatial resolutions, different measurement error biases and variances, and other disparate characteristics. A state-of-the-art spatial statistical model was used to relate the true, but not directly observed, geophysical field to noisy, spatial aggregates observed by remote sensing instruments. The spatial covariances of the true field and the covariances of the true field with the observations were modeled. The observations are spatial averages of the true field values, over pixels, with different measurement noise superimposed. A kriging framework is used to infer optimal (minimum mean squared error and unbiased) estimates of the true field at point locations from pixel-level, noisy observations. A key feature of the spatial statistical model is the spatial mixed effects model that underlies it. The approach models the spatial covariance function of the underlying field using linear combinations of basis functions of fixed size. Approaches based on kriging require the inversion of very large spatial covariance matrices, and this is usually done by making simplifying assumptions about spatial covariance structure that simply do not hold for geophysical variables. In contrast, this method does not require these assumptions, and is also computationally much faster. This method is fundamentally different than other approaches to data fusion for remote sensing data because it is inferential rather than merely descriptive. All approaches combine data in a way that minimizes some specified loss function. Most of these are more or less ad hoc criteria based on what looks good to the eye, or some criteria that relate only to the data at hand.

Braverman, Amy J.↗

Use of observational and model-derived fields and regime model output statistics in mesoscale forecasting

Various empirical and statistical weather-forecasting studies which utilize stratification by weather regime are described. Objective classification was used to determine weather regime in some studies. In other cases the weather pattern was determined on the basis of a parameter representing the physical and dynamical processes relevant to the anticipated mesoscale phenomena, such as low level moisture convergence and convective precipitation, or the Froude number and the occurrence of cold-air damming. For mesoscale phenomena already in existence, new forecasting techniques were developed. The use of cloud models in operational forecasting is discussed. Models to calculate the spatial scales of forcings and resultant response for mesoscale systems are presented. The use of these models to represent the climatologically most prevalent systems, and to perform case-by-case simulations is reviewed. Operational implementation of mesoscale data into weather forecasts, using both actual simulation output and method-output statistics is discussed.

Forbes, G. S.↗

Multi-source Estimates of Land / Ocean Moisture Transport Variability over the Satellite Era

It is widely appreciated that atmospheric transport of water from the world’s oceans is a process key to planetary energy balance as well as Earth's habitability. What is not yet clear is the extent of variability in moisture transports, the relative importance of interdecadal variability versus climate change signals, and importantly, our ability to quantify these changes. This work assesses variations in moisture transport variability during the satellite era (~1980 to present) by comparing several different estimates. (i) The most direct estimate is the vertically integrated flux convergence of moisture from reanalyses which use observed wind and moisture information. (ii) One alternative estimate comes from P-ET over land taken from global hydrologic models constrained with precipitation and near-surface meteorology. Here we use an ensemble of six models. An adjunct to this method is to employ satellite derive ET (e.g., GLEAM or DOLCE). (iii) Complementary to this is E-P over the global oceans derived from satellite estimates of P such as TRMM, GPM and GPCP and SeaFlux V3 or J-OFURO3 estimates of E, all relying heavily upon microwave measurements. Transport between land and oceans must essentially balance at monthly scales, i.e., vanish globally. (iv) a fourth perspective comes from estimate of terrestrial RO + storage rate, delta S. G-RUN Ensemble which uses observed streamflow and P measurements to calibrate a statistical model provides the former while GRACE, GRACE-FO provide total water storage anomalies used to calculate storage rate changes. GRACE REC uses GRACE data to train a precipitation-driven statistical model to extend storage estimates before the GRACE era. (All of these alternatives to reanalysis estimates also consider the small atmospheric column water vapor contribution.) We examine the transport changes from these three different methodologies, their relative accuracies and discuss the origin of their differences. Regional trends in moisture flux divergence and their role in multi-decadal trends are considered. Interannual variability arising in connection with ENSO variability is a dominant signal, driven largely by P changes. Trends since 1980 include reductions in moisture delivery to the western U.S., eastern Brazil, and central Africa with recovery of moisture convergence to the Sahel and parts of eastern North America.

Franklin Robertson↗

Statistical Projection of Solar Cycle 24 for the Exposure Estimates

A solar cycle statistical model has been developed based on the accumulating cycle sunspot data to estimate future levels of the solar cycle activity. Since the current solar cycle 24 has progressed about three years, the cycle activity levels are estimated with an accurately defined solar minimum 24. Then, solar cycle 24 is projected with the cycle activity levels using the statistical model. The projection of solar cycle 24 is then coupled to space related quantities of interest to radiation protection, because the interplanetary plasma and radiation fields are modulated by the degree of disturbance in the solar surface and the radiation doses received by astronauts in interplanetary space are likewise influenced. The resultant projection of solar cycle 24 provides a basis for estimating exposure in future space missions, and projection errors can be corrected as the cycle progresses and observations become available because this model is shown to be self-correcting.

Kim, Myung-Hee↗

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

The AutoBayes Program Synthesis System: System Description

AUTOBAYES is a fully automatic program synthesis system for the statistical data analysis domain. Its input is a concise description of a data analysis problem in the form of a statistical model; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. AUTOBAYES synthesizes code by a schema-guided deductive process. Schemas (i.e., code templates with associated semantic constraints) are applied to the original problem and recursively to emerging subproblems. AUTOBAYES complements this approach by symbolic computation to derive closed-form solutions whenever possible. In this paper, we concentrate on the interaction between the symbolic computations and the deductive synthesis process. A statistical model specifies for each problem variable (i.e., data or parameter) its properties and dependencies in the form of a probability distribution, A typical data analysis task is to estimate the best possible parameter values from the given observations or measurements. The following example models normal-distributed data but takes prior information (e.g., from previous experiments) on the data's mean value and variance into account.

Fischer, Bernd↗

Information retrieval from wide-band meteorological data - An example

The methods proposed by Smith and Adelfang (1981) and Smith et al. (1982) are used to calculate probabilities over rectangles and sectors of the gust magnitude-gust length plane; probabilities over the same regions are also calculated from the observed distributions and a comparison is also presented to demonstrate the accuracy of the statistical model. These and other statistical results are calculated from samples of Jimsphere wind profiles at Cape Canaveral. The results are presented for a variety of wavelength bands, altitudes, and seasons. It is shown that wind perturbations observed in Jimsphere wind profiles in various wavelength bands can be analyzed by using digital filters. The relationship between gust magnitude and gust length is modeled with the bivariate gamma distribution. It is pointed out that application of the model to calculate probabilities over specific areas of the gust magnitude-gust length plane can be useful in aerospace design.

Adelfang, S. I.↗

System statistical reliability model and analysis

A digital computer code was developed to simulate the time-dependent behavior of the 5-kwe reactor thermoelectric system. The code was used to determine lifetime sensitivity coefficients for a number of system design parameters, such as thermoelectric module efficiency and degradation rate, radiator absorptivity and emissivity, fuel element barrier defect constant, beginning-of-life reactivity, etc. A probability distribution (mean and standard deviation) was estimated for each of these design parameters. Then, error analysis was used to obtain a probability distribution for the system lifetime (mean = 7.7 years, standard deviation = 1.1 years). From this, the probability that the system will achieve the design goal of 5 years lifetime is 0.993. This value represents an estimate of the degradation reliability of the system.

Lekach, V. S.↗

Statistical-band-model analysis and integrated intensity for the 21.8 micrometer bands of HNO3 vapor

A large number of low-resolution quantitative spectra were obtained for a band-model analysis of 21.8-micrometer bands. The experimental investigations were conducted with pure nitric-acid vapor at 40 C, taking into account pressures in the range from a few torr to a pressure which is close to the saturation pressure. It is pointed out that the obtained data, in addition to those from the 11.3-micrometer band, can be used for an independent determination of the nitric acid vertical distribution in the stratosphere.

Goldman, A.↗

Evaluation of the radar response to land surfaces and volumes: Examination of theoretical models, target statistics, and applications

Four areas of L-band radar remote sensing of terrain were examined: (1) the behavior of the radar backscatter coefficient of distributed surface and volumes as a function of the targets' dielectric and geometric parameters and as a fnction of their physical parameters; (2) the correspondence of the angular behavior of the relative backscatter coefficient as extracted from SIR-B digital imagery and truck mounted L-band scatterometer measurements for about 100 fields; (3) the statistical behavior of SIR-B image density for targets that appear homogeneous on Thematic Mapper (TM) optical imagery and/or color IR photography; and (4) the applicability of SIR-B imagery both alone and in conjunction with TM imagery for the classification and monitoring of land cover and renewable resources.

Ulaby, F. T.↗

Mechanisms of crustal deformation in the western US

The deformation processes in the western United States were studied, considering both deterministic models and random or statistical models. The role of the intracrustal delamination and mechanisms of crustal thinning were also examined. The application of fractal techniques to understand how the crust is deforming was studied in complex regions. Work continued on the development of a fractal based model for deformation in the western United States. Fractal studies were also extended to the study of topography and the geoid.

Turcotte, Donald L.↗

Application of large eddy interaction model to a mixing layer

The large eddy interaction model (LEIM) is a statistical model of turbulence based on the interaction of selected eddies with the mean flow and all of the eddies in a turbulent shear flow. It can be utilized as the starting point for obtaining physical structures in the flow. The possible application of the LEIM to a mixing layer formed between two parallel, incompressible flows with a small temperature difference is developed by invoking a detailed similarity between the spectra of velocity and temperature.

Murthy, S. N. B.↗

Development of a two-dimensional zonally averaged statistical-dynamical model. III - The parameterization of the eddy fluxes of heat and moisture

A number of perpetual January simulations are carried out with a two-dimensional zonally averaged model employing various parameterizations of the eddy fluxes of heat (potential temperature) and moisture. The parameterizations are evaluated by comparing these results with the eddy fluxes calculated in a parallel simulation using a three-dimensional general circulation model with zonally symmetric forcing. The three-dimensional model's performance in turn is evaluated by comparing its results using realistic (nonsymmetric) boundary conditions with observations. Branscome's parameterization of the meridional eddy flux of heat and Leovy's parameterization of the meridional eddy flux of moisture simulate the seasonal and latitudinal variations of these fluxes reasonably well, while somewhat underestimating their magnitudes. New parameterizations of the vertical eddy fluxes are developed that take into account the enhancement of the eddy mixing slope in a growing baroclinic wave due to condensation, and also the effect of eddy fluctuations in relative humidity. The new parameterizations, when tested in the two-dimensional model, simulate the seasonal, latitudinal, and vertical variations of the vertical eddy fluxes quite well, when compared with the three-dimensional model, and only underestimate the magnitude of the fluxes by 10 to 20 percent.

Stone, Peter H.↗