Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameter estimation using data”

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 901 records · Page 50

Improved Soundings and Error Estimates using AIRS/AMSU Data

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of 1 K, and layer precipitable water with an rms error of 20 percent, in cases with up to 80 percent effective cloud cover. The basic theory used to analyze AIRS/AMSU/HSB data in the presence of clouds, called the at-launch algorithm, and a post-launch algorithm which differed only in the minor details from the at-launch algorithm, have been described previously. The post-launch algorithm, referred to as AIRS Version 4.0, has been used by the Goddard DAAC to analyze and distribute AIRS retrieval products. In this paper we show progress made toward the AIRS Version 5.0 algorithm which will be used by the Goddard DAAC starting late in 2006. A new methodology has been developed to provide accurate case by case error estimates for retrieved geophysical parameters and for the channel by channel cloud cleared radiances used to derive the geophysical parameters from the AIRS/AMSU observations. These error estimates are in turn used for quality control of the derived geophysical parameters and clear column radiances. Improvements made to the retrieval algorithm since Version 4.0 are described as well as results comparing Version 5.0 retrieval accuracy and spatial coverage with those obtained using Version 4.0.

Susskind, Joel↗

In vivo kinetics of SARS-CoV-2 infection and its relationship with a person’s infectiousness

Significance Quantifying the kinetics of SARS-CoV-2 infection and individual infectiousness is important for understanding SARS-CoV-2 transmission and evaluating intervention strategies. Here, we developed within-host models of SARS-CoV-2 infection, and by fitting them to clinical data, we estimated key within-host viral dynamic parameters. We also developed a mechanistic model for viral transmission and show that the logarithm of the viral load in the upper respiratory tract serves as an appropriate surrogate for a person’s infectiousness. Using data on how viral load changes during infection, we further evaluated the effectiveness of PCR and antigen-based testing strategies for averting transmission and identifying infected individuals.

60 APPLIED LIFE SCIENCES↗

MFNets: Multifidelity data-driven networks for Bayesian learning and prediction

This paper presents a multifidelity uncertainty quantification framework called MFNets. We seek to address three existing challenges that arise when experimental and simulation data from different sources are used to enhance statistical estimation and prediction with quantified uncertainty. Specifically, we demonstrate that MFNets can (1) fuse heterogeneous data sources arising from simulations with different parameterizations, e.g simulation models with different uncertain parameters or data sets collected under different environmental conditions; (2) encode known relationships among data sources to reduce data requirements; and (3) improve the robustness of existing multi-fidelity approaches to corrupted data. MFNets construct a network of latent variables (LVs) to facilitate the fusion of data from an ensemble of sources of varying credibility and cost. These LVs are posited as explanatory variables that provide the source of correlation in the observed data. Furthermore, MFNets provide a way to encode prior physical knowledge to enable efficient estimation of statistics and/or construction of surrogates via conditional independence relations on the LVs. We highlight the utility of our framework with a number of theoretical results which assess the quality of the posterior mean as a frequentist estimator and compare it to standard sampling approaches that use single fidelity, multilevel, and control variate Monte Carlo estimators. We also use the proposed framework to derive the Monte Carlo-based control variate estimator entirely from the use of Bayes rule and linear-Gaussian models -- to our knowledge the first such derivation. Finally, we demonstrate the ability to work with different uncertain parameters across different models.

97 MATHEMATICS AND COMPUTING↗

Evaluating the Feasibility of Data Simulation for Radiation Research: Estimating RBE Values for HZE Radiation

As the drive for deep space exploration intensifies, a comprehensive understanding of the health effects of radiation exposure becomes paramount to the future of human space flight. However, epidemiological data for radiation exposure, particularly to high-energy (HZE) ions, is limited, partially due to the financial and logistical costs of radiation studies. As an alternative, this study aims to assess the viability of data simulation strategies to accurately model potential study parameters prior to utilizing laboratory conditions. This study estimates a relative biological effectiveness (RBE) factor based on the solid tumor data for outbred mice provided by Edmundson et al. 2020. Excess relative risk (ERR) models for HZE-Fe ions and gamma radiation were estimated using Poisson regression with Weibull models to represent the background solid tumor hazard without radiation. RBE values were calculated from the ratio of the heavy ion linear slope to the gamma linear slope. The parameters from these models were then applied to simulate iterations of 300 datasets across HZE-Fe doses of 0.05, 0.2, 0.4, and 0.75 Gy; gamma radiation doses of 0.75, 2, and 3 Gy, as well as an unirradiated control group. Sample size per dose varied from 100-500 mice across simulations (800-4000 mice total per dataset). 1500 datasets were generated total (300 each for sample sizes 800, 1600, 3200, and 4000). For each dataset, ERR per radiation type and RBE for HZE-Fe were calculated. The RBE from the Edmondson data was calculated to be 5.55. The RBEs from the simulated data converged around this value across the increasing sample sizes. For n = 4000, the mean RBE across the 300 datasets was 5.71 (95% CI: 5.47-5.95). These results suggest that RBEs calculated from simulated data are credible. Based on this exploration, data simulation is a viable method of testing radiation studies. It provides a method of testing study assumptions and refining research questions prior to embarking on costly laboratory experiments.

Sneha Mehta↗

Estimating Geophysical Parameters From Gravity Data

ORBSIM program developed for accurate extraction of parameters of geophysical models from Doppler-radio-tracking data acquired from orbiting planetary spacecraft. Model of proposed planetary structure used in numerical integration along simulated trajectories of spacecraft around primary body. Written in FORTRAN 77.

Sjogren, William L.↗

Estimating a Service-Life Distribution Based on Production Counts and a Failure Database

A manufacturer wanted to compare the service-life distributions of two similar products. These concern product lifetimes after installation (not manufacture). For each product, there were available production counts and an imperfect database providing information on failing units. In the real case, these units were expensive repairable units warrantied against repairs. Failure (of interest here) was relatively rare and driven by a different mode/mechanism than ordinary repair events (not of interest here). Approach: Data models for the service life based on a standard parametric lifetime distribution and a related limited failure population were developed. These models were used to develop expressions for the likelihood of the available data that properly accounts for information missing in the failure database. Results: A Bayesian approach was employed to obtain estimates of model parameters (with associated uncertainty) in order to investigate characteristics of the service-life distribution. Custom software was developed and is included as Supplemental Material to this case study. One part of a responsible approach to the original case was a simulation experiment used to validate the correctness of the software and the behavior of the statistical methodology before using its results in the application, and an example of such an experiment is included here. Because of confidentiality issues that prevent use of the original data, simulated data with characteristics like the manufacturer’s proprietary data are used to illustrate some aspects of our real analyses. Lastly, we also note that, although this case focuses on rare and complete product failure, the statistical methodology provided is directly applicable to more standard warranty data problems involving typically much larger warranty databases where entries are warranty claims (often for repairs) rather than reports of complete failures.

97 MATHEMATICS AND COMPUTING↗

The Community Radiative Transfer Model (CRTM): Community-Focused Collaborative Model Development Accelerating Research to Operations

The Joint Center for Satellite Data Assimilation (JCSDA) Community Radiative Transfer Model (CRTM) is a fast, 1-D radiative transfer model used in numerical weather prediction, calibration/validation, etc. across multiple federal agencies and universities. The key benefit of the CRTM is that it is a satellite simulator. It provides a highly accurate representation of satellite radiances by using the specific sensor response functions convolved with a line-by-line radiative transfer model (LBLRTM). CRTM covers the spectral ranges consistent with all present operational and most research satellites, from visible to microwave. The capability to simulate ultraviolet radiances and support space-based radar sensors is being added over the next two years in CRTM Version 3.0. In addition to simulated radiances, the CRTM also provides Jacobian outputs needed to interpret satellite observations for numerical weather prediction. The Jacobian estimates how changes in geophysical parameters affect simulated measurements from satellite sensors. Using the Jacobian in modeling and weather prediction improves the accuracy and efficiency of data analysis, leading to better weather predictions. The CRTM model's success and growth depend on community contributions and evaluation. To facilitate this, we have made the CRTM highly accessible through modular programming, clear documentation and tutorials, public domain licensing, unfettered public access via Github, and a clear path to operational implementation for innovative research. We encourage and welcome contributions from the community to help us continue to improve the CRTM.

Benjamin T Johnson↗

Theoretical branching ratios for the 5I7 to 5I7 levels of Ho(3+) in the garnets A3B2C3O12 (A = Y,La,Lu,Gd; B = Al,Lu,Sc,Ga; C = Al,Ga)

Results are reported from an experimental study investigating triply ionized holmium in 10 garnets using the point-change model to predict theoretical energy levels and temperature-dependent branching ratios for the 5I7 to 5I8 manifolds for temperatures between 50 and 400 K. Plots were made for the largest lines at 300 K. YScAG was plotted twice, once for each set of X-ray data available. Energy levels are predicted based on theoretical crystal-field parameters, and good agreement to experiment is found. It is suggested that the present set of theoretical crystal-field parameters provides good estimates of the energy levels for the other hosts on which there are no experimental optical data. X-ray and index-of-refraction data are used to evaluate the performance of 10 lasers via a quantum mechanical model to predict the position of the energy levels and the temperature-dependent branching rations of the 5I7 to 5I8 levels of holmium. The fractional population inversion required for threshold is also evaluated.

Filer, Elizabeth D.↗

Multi-area parameter error identification for large power systems

Power grid model parameters may contain errors due to various reasons. Detecting and correcting parameter errors typically requires significant computational effort due to the size and complexity of the parameter database. While the normalized Lagrange multiplier (NLM) method can effectively detect, identify and correct parameter errors, its computational burden could rapidly grow with increasing system size. This paper addresses this issue by proposing a multi-area parameter error identification method. Each area has its own outlier detection tool for detecting the incorrect parameters and measurements within the area. On the other hand, due to the reduced redundancy at area boundaries, parameter errors on branches incident to boundary buses may not be detected. Such errors are subsequently detected by a coordination level estimator completing the system-wide parameter detection procedure. In conclusion, performance of the developed method is demonstrated using the IEEE 118-bus and 2000-bus Texas synthetic systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Robust State Estimation Framework Considering Measurement Correlations and Imperfect Synchronization

This paper develops a robust power system state estimation framework with the consideration of measurement correlations and imperfect synchronization. In the framework, correlations of SCADA and PMU measurements are calculated separately through unscented transformation and a vector autoregression (VAR) model. In particular, PMU measurements during the waiting period of two SCADA measurement scans are buffered to develop the VAR model. The latter takes into account the temporal and spatial correlations of PMU measurements and provides redundant measurements to suppress bad data and mitigate imperfect synchronization. In case of synchronization issues between SCADA and PMU measurements, either forecasted PMU measurements or prior SCADA measurements from the last estimation run are leveraged to restore system observability. After that, robust generalized maximum-likelihood (GM)-estimator is developed to integrate full measurement correlations and handle bad data in SCADA and PMU measurements. Thanks to the robustness and flexibility of the proposed robust estimator, bad data as well as imperfect measurement synchronization are mitigated, yielding high statistical efficiency. Comprehensive comparison results with other alternatives under various conditions demonstrate the benefits of the proposed framework.

Zhao, Junbo (ORCID:0000000284989666)↗

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

Quantifying Model Form Uncertainty of Epidemic Forecasting Models from Incidence Data

Objective We present a mathematical framework for non-parametric estimation of the force of infection, together with statistical upper and lower confidence bands. The resulting estimates allow to assess how well simpler models, such as SEIR, fit the observed time series of incidence data.IntroductionUncertainty Quantification (UQ), the ability to quantify the impact of sample-to-sample variations and model misspecification on predictions and forecasts, is a critical aspect of disease surveillance. While quantifying the impact of stochastic uncertainty in the data is well understood, quantifying the impact of model misspecification is significantly harder. For the latter, one needs a "universal model" to which more restrictive parametric models are compared too.MethodsThis talk presents a useful modeling framework for time series of incidence data from contagious diseases that enables one to identify and quantify the impact of model form uncertainty. Specifically, we propose to focus on estimating the timedependent force of infection. The latter is a universal parameters for all contagious disease model. Using a machine learning technique for estimating monotone functions, i.e., isotonic regression and its variants, one can estimate the force ofinfection without addtional assumptions. We note that most contagious disease model do satisfy this monotonicity assumption, due to a combination of factors: depletion of susceptibles, implementation of mitigation strategies, behavior change, etc. Comparing the resulting "non-parametric" estimate with parametric estimates, obtained by fitting an SEIR for example, can reveal model deficiencies and help quantify model form uncertainties.Finally, we discuss how ideas from "strict bound theory" can be used to develop upper and lower uncertainty bands for force of infection that acknowledge the intrinsic stochasticity in the data.ResultsWe demonstrate the application of the methodology to weekly Influenza Like Illness (ILI) incidence data from France andcompare the results to fitted SIR and SEIR models. This comparison can be seen as a nonparametric goodness of fit test, providing one with tools to do simple model selection.ConclusionsWe present a novel and flexible model to statistically describe the force of infection as a function of time. Comparing the fit to incidence data of that model with the fit of simpler parametric models enables the quantification of model form uncertainty and associated model selection.

Hengartner, Nicholas↗

A geopotential model from satellite tracking, altimeter, and surface gravity data: GEM-T3

An improved model of Earth's gravitational field, Goddard Earth Model T-3 (GEM-T3), has been developed from a combination of satellite tracking, satellite altimeter, and surface gravimetric data. GEM-T3 provides a significant improvement in the modeling of the gravity field at half wavelengths of 400 km and longer. This model, complete to degree and order 50, yields more accurate satellite orbits and an improved geoid representation than previous Goddard Earth Models. GEM-T3 uses altimeter data from GEOS 3 (1975-1976), Seasat (1978) and Geosat (1986-1987). Tracking information used in the solution includes more than 1300 arcs of data encompassing 31 different satellites. The recovery of the long-wavelength components of the solution relies mostly on highly precise satellite laser ranging (SLR) data, but also includes Tracking Network (TRANET) Doppler, optical, and satellite-to-satellite tracking acquired between the ATS 6 and GEOS 3 satellites. The main advances over GEM-T2 (beyond the inclusion of altimeter and surface gravity information which is essential for the resolution of the shorter wavelength geoid) are some improved tracking data analysis approaches and additional SLR data. Although the use of altimeter data has greatly enhanced the modeling of the ocean geoid between 65 deg N and 60 deg S latitudes in GEM-T3, the lack of accurate detailed surface gravimetry leaves poor geoid resolution over many continental regions of great tectonic interest (e.g., Himalayas, Andes). Estimates of polar motion, tracking station coordinates, and long-wavelength ocean tidal terms were also made (accounting for 6330 parameters). GEM-T3 has undergone error calibration using a technique based on subset solutions to produce reliable error estimates. The calibration is based on the condition that the expected mean square deviation of a subset gravity solution from the full set values is predicted by the solutions' error covariances. Data weights are iteratively adjusted until this condition for the error calibration is satisfied. In addition, gravity field tests were performed on strong satellite data sets withheld from the solution (thereby ensuring their independence). In these tests, the performance of the subset models on the withheld observations is compared to error projections based on their calibrated error covariances. These results demonstrate that orbit accuracy projections are reliable for new satellites which were not included in GEM-T3.

Lerch, F. J.↗

Fatigue behavior and modeling for additive manufactured 304L stainless steel: The effect of surface roughness

The fatigue strength of additively manufactured metallic parts in their as-built surface condition is mainly dominated by the surface roughness. Post-processing is often inevitable to reduce surface roughness effects even though post-processing diminishes the main advantage of additive manufacturing, which is net-shaped direct-to-service production. This work investigates the underlying mechanisms responsible for fatigue failure of additively manufactured 304L stainless steel parts in as-built and machined/polished surface conditions. Both strain- and force-controlled, fully reversed fatigue tests were conducted to gain a comprehensive understanding of surface roughness effects on fatigue behavior. The sensitivity to surface roughness is shown to be dependent on the control mode, with stress-based fatigue tests showing greater sensitivity than strain-based fatigue tests. Moreover, the fatigue life estimation for as-built specimens was performed based on surface roughness parameters as well as the fatigue properties of the specimens in machined/polished surface condition of the material without using any fatigue data of specimens in as-built surface condition. Accordingly, the effect of surface roughness on the fatigue behavior could be estimated reasonably well in both strain-life and stress-life approaches.

36 MATERIALS SCIENCE↗

Wind field model-based estimation of Seasat scatterometer winds

A model-based approach to estimating near-surface wind fields over the ocean from Seasat scatterometer (SASS) measurements is presented. The approach is a direct assimilation technique in which wind field model parameters are estimated directly from the scatterometer measurements of the radar backscatter of the ocean's surface using maximum likelihood principles. The wind field estimate is then computed from the estimated model parameters. The wind field model used in this approach is based on geostrophic approximation and on simplistic assumptions about the wind field vorticity and divergence but includes ageostrophic winds. Nine days of SASS data were processed to obtain unique wind estimates. Comparisons in performance to the traditional two-step (point-wise wind retrieval followed by ambiguity removal) wind estimate method and the model-based method are provided using both simulated radar backscatter measurements and actual SASS measurements. In the latter case the results are compared to wind fields determined using subjective ambiguity removal. While the traditional approach results in missing measurements and reduced effective swath width due to fore/aft beam cell coregistration problems, the model-based approach uses all available measurements to increase the effective swath width and to reduce data gaps. The results reveal that the model-based wind estimates have accuracy comparable to traditionally estimated winds with less 'noise' in the directional estimates, particularly at low wind speeds.

Long, David G.↗

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗

CONUS and Global Carbon Consumption Database for Wildland Fire

Fire plays a significant role on both national and global scales, profoundly impacting landscapes shaped by human activity as well as wildlands. Even though fire can be devastating, wildland fire is a natural and integral force on our landscapes. Fires can also serve to reduce fuels to mitigate wildfire risk and maintain healthy ecosystem functions. However, the smoke produced by fires, regardless of their size or purpose, can pose adverse effects on human health when inhaled downwind. Understanding the influence of smoke on air quality and human well-being necessitates the quantification of emissions that fires release into the atmosphere. In response to this need, we have developed a carbon consumption database for the Continental United States (CONUS) and are developing a global carbon consumption database, both of which are directly link to distinct fuels within various fire danger categories. Fuel Characteristic Classification System (FCCS) fuels are used to parameterize biomass, and consumption is broken down into five Fire Danger categories (Low, Moderate, High, Very High, Extreme), for both ‘new’ and ‘residual’ burning scenarios. Residual burned area is defined as burning in areas that have recently burned. We implemented the FCCS 30-meter CONUS fuelbed dataset during the 2019 Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) campaign to estimate daily carbon emissions. Our emissions estimates were rigorously compared against in-situ measurements of CO2, CO, and black carbon aerosols, revealing a robust agreement between the two datasets. The 300-meter global database builds upon the foundations of the Pettinari, M. Lucrecia (2015) Global Fuelbed database, a global fuel map with standardized FCCS biomass parameters. These products can serve as valuable tools when used in conjunction with burned area data to rapidly and accurately estimate carbon consumed and released into the atmosphere.

FCCS↗

Multistage variable probability forest volume inventory

An inventory scheme based on the use of computer processed LANDSAT MSS data was developed. Output from the inventory scheme provides an estimate of the standing net saw timber volume of a major timber species on a selected forested area of the Navajo Nation. Such estimates are based on the values of parameters currently used for scaled sawlog conversion to mill output. The multistage variable probability sampling appears capable of producing estimates which compare favorably with those produced using conventional techniques. In addition, the reduction in time, manpower, and overall costs lend it to numerous applications.

Anderson, J. E.↗