Search NASA⌕ Search

SEARCH · Search NASA

Results for “covariance selection”

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 19 records

Stable and verifiable state estimation methods and systems with spacecraft applications

The stability of a recursive estimator process (e.g., a Kalman filter is assured for long time periods by periodically resetting an error covariance P(t.sub.n) of the system to a predetermined reset value P.sub.r. The recursive process is thus repetitively forced to start from a selected covariance and continue for a time period that is short compared to the system's total operational time period. The time period in which the process must maintain its numerical stability is significantly reduced as is the demand on the system's numerical stability. The process stability for an extended operational time period T.sub.o is verified by performing the resetting step at the end of at least one reset time period T.sub.r whose duration is less than the operational time period T.sub.o and then confirming stability of the process over the reset time period T.sub.r. Because the recursive process starts from a selected covariance at the beginning of each reset time period T.sub.r, confirming stability of the process over at least one reset time period substantially confirms stability over the longer operational time period T.sub.o.

Li, Rongsheng↗

Pneumatic distortion compensation for aircraft surface pressure sensing devices

In this paper a technique of compensating for pneumatic distortion in aircraft surface pressure sensing devices is developed. The compensation allows conventional pressure sensing technology to obtain improved unsteady pressure measurements. Pressure distortion caused by frictional attenuation and pneumatic resonance within the sensing system makes obtaining unsteady pressure measurements by conventional sensors difficult. Typically, most of the distortion occurs within the pneumatic tubing used to transmit pressure impulses from the surface of the aircraft to the measurement transducer. This paper develops a second-order distortion model that accurately describes the behavior of the primary wave harmonic of the pneumatic tubing. The model is expressed in state-variable form and is coupled with standard results from minimum-variance estimation theory to develop an algorithm to compensate for the effects of pneumatic distortion. Both postflight and real-time algorithms are developed and evaluated using simulated and flight data. Covariance selection and filter-tuning examples are presented. Results presented verify that, given appropriate covariance magnitudes, the algorithms accurately reconstruct surface pressure values from remotely sensed pressure measurements.

Whitmore, Stephen A.↗

Mariner/Jupiter/Saturn navigation in the presence of massive planetary satellites

Orbit determination accuracies attainable during a Mariner/Jupiter/Saturn encounter have been established by means of covariance analyses. Advanced earth-based multistation radiometric data, optical data consisting of star/satellite pictures provided by a narrow-angle TV camera on board the spacecraft, and batch sequential filtering methods were employed. An evaluation of sequential filter sensitivities to errors in modeling small nongravitational spacecraft accelerations is presented. Selected covariances developed in the orbit determination analysis are used in an investigation of the trajectory correction costs associated with a close Ganymede encounter. The 99 percentile Delta V contours in the Ganymede aiming plane are obtained for satellite encounters before or after Jupiter closed approach, and the sources of the dominant contributions to the Delta V costs are identified.

Hildebrand, C. E.↗

Generic Kalman Filter Software

The Generic Kalman Filter (GKF) software provides a standard basis for the development of application-specific Kalman-filter programs. Historically, Kalman filters have been implemented by customized programs that must be written, coded, and debugged anew for each unique application, then tested and tuned with simulated or actual measurement data. Total development times for typical Kalman-filter application programs have ranged from months to weeks. The GKF software can simplify the development process and reduce the development time by eliminating the need to re-create the fundamental implementation of the Kalman filter for each new application. The GKF software is written in the ANSI C programming language. It contains a generic Kalman-filter-development directory that, in turn, contains a code for a generic Kalman filter function; more specifically, it contains a generically designed and generically coded implementation of linear, linearized, and extended Kalman filtering algorithms, including algorithms for state- and covariance-update and -propagation functions. The mathematical theory that underlies the algorithms is well known and has been reported extensively in the open technical literature. Also contained in the directory are a header file that defines generic Kalman-filter data structures and prototype functions and template versions of application-specific subfunction and calling navigation/estimation routine code and headers. Once the user has provided a calling routine and the required application-specific subfunctions, the application-specific Kalman-filter software can be compiled and executed immediately. During execution, the generic Kalman-filter function is called from a higher-level navigation or estimation routine that preprocesses measurement data and post-processes output data. The generic Kalman-filter function uses the aforementioned data structures and five implementation- specific subfunctions, which have been developed by the user on the basis of the aforementioned templates. The GKF software can be used to develop many different types of unfactorized Kalman filters. A developer can choose to implement either a linearized or an extended Kalman filter algorithm, without having to modify the GKF software. Control dynamics can be taken into account or neglected in the filter-dynamics model. Filter programs developed by use of the GKF software can be made to propagate equations of motion for linear or nonlinear dynamical systems that are deterministic or stochastic. In addition, filter programs can be made to operate in user-selectable "covariance analysis" and "propagation-only" modes that are useful in design and development stages.

Lisano, Michael E., II↗

Sensitivity of Lower Stratospheric Assimilated Ozone on Error Covariance Modeling and Data Selection

Assimilated ozone is produced at the NASA/Goddard Data Assimilation Office by blending ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line transport model. The current system tends to overestimate the amount of lower stratospheric ozone. This is a region where ozone plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. Since these are all important applications of assimilated ozone products, effort is being directed towards reducing this bias. The SBUV ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone in the lower stratosphere, and its vertical distribution in particular, are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Accordingly, the sensitivity of the assimilated lower stratospheric ozone fields to changes in the TOMS error-covariance modeling and the SBUV data selection has been investigated. The use of a spatially correlated TOMS error covariance model led to improvements in the product. However, withholding the SBUV/2 data for the layer between 63 and 126 hPa typically degraded the product, a result which vindicates the use of this layer ozone product, despite its known errors. These efforts to improve the lower stratospheric distribution will be extended to include a more advanced forecast error covariance model, and by assimilating ozone products from new instruments on Envisat and EOS Aura.

Stajner, Ivanka↗

Least squares collocation applied to local gravimetric solutions from satellite gravity gradiometry data

An autonomous spaceborne gravity gradiometer mission is being considered as a post Geopotential Research Mission project. The introduction of satellite diometry data to geodesy is expected to improve solid earth gravity models. The possibility of utilizing gradiometer data for the determination of pertinent gravimetric quantities on a local basis is explored. The analytical technique of least squares collocation is investigated for its usefulness in local solutions of this type. It is assumed, in the error analysis, that the vertical gravity gradient component of the gradient tensor is used as the raw data signal from which the corresponding reference gradients are removed to create the centered observations required in the collocation solution. The reference gradients are computed from a high degree and order geopotential model. The solution can be made in terms of mean or point gravity anomalies, height anomalies, or other useful gravimetric quantities depending on the choice of covariance types. Selected for this study were 30 x 30 foot mean gravity and height anomalies. Existing software and new software are utilized to implement the collocation technique. It was determined that satellite gradiometry data at an altitude of 200 km can be used successfully for the determination of 30 x 30 foot mean gravity anomalies to an accuracy of 9.2 mgal from this algorithm. It is shown that the resulting accuracy estimates are sensitive to gravity model coefficient uncertainties, data reduction assumptions and satellite mission parameters.

Robbins, J. W.↗

Improved effectiveness of GPS RAIM through ridge regression signal processing

A measurement processing method has been developed which markedly improves the GPS Receiver Autonomous Integrity Monitoring (RAIM) software-based algorithm system's effectiveness in detecting satellite signal failures. Detection is via the consistency of a redundant set of pseudorange measurements. When five satellites are in view, five different subsolutions can be calculated; the integrity alarm is triggered on the basis of subsolution comparisons. Because a poor distribution of the satellites also causes RAIM subsolution scattering, a methodology for selecting the covariance matrix is presented which incorporates ridge regression into a Kalman filter.

Kelly, Robert J.↗

Estimation of Model Error Variances During Data Assimilation

Data assimilation is all about understanding the error characteristics of the data and models that are used in the assimilation process. Reliable error estimates are needed to implement observational quality control, bias correction of observations and model fields, and intelligent data selection. Meaningful covariance specifications are obviously required for the analysis as well, since the impact of any single observation strongly depends on the assumed structure of the background errors. Operational atmospheric data assimilation systems still rely primarily on climatological background error covariances. To obtain error estimates that reflect both the character of the flow and the current state of the observing system, it is necessary to solve three problems: (1) how to account for the short-term evolution of errors in the initial conditions; (2) how to estimate the additional component of error caused by model defects; and (3) how to compute the error reduction in the analysis due to observational information. Various approaches are now available that provide approximate solutions to the first and third of these problems. However, the useful accuracy of these solutions very much depends on the size and character of the model errors and the ability to account for them. Model errors represent the real-world forcing of the error evolution in a data assimilation system. Clearly, meaningful model error estimates and/or statistics must be based on information external to the model itself. The most obvious information source is observational, and since the volume of available geophysical data is growing rapidly, there is some hope that a purely statistical approach to model error estimation can be viable. This requires that the observation errors themselves are well understood and quantifiable. We will discuss some of these challenges and present a new sequential scheme for estimating model error variances from observations in the context of an atmospheric data assimilation system.

Dee, Dick↗

Techniques in Assessing Random Uncertainty of Wind Tunnel Replicate Test Data

The preferred method to determine the random uncertainty of variables of interest during wind tunnel tests is by direct analysis of replicate data. However, current information is limited to analysis of a single replicate set condition that does not consider data covariance. This paper provides additional information on considerations for quantifying random uncertainty including the evaluation and correction of replicate data covariance, the selection of replicate set conditions for the test, and the interpolation of random uncertainty at other than replicate set conditions. The analysis and techniques are supported by data from a conventional aircraft configuration test at the Oran Nicks Low-Speed Wind Tunnel and an ASME nozzle flow characterization at the Propulsion Systems Laboratory.

Uncertainty↗

Population Pharmacokinetics of Intranasal Scopolamine

Introduction: An intranasal gel dosage formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness (SMS).The bioavailability and pharmacokinetics (PK) was evaluated using data collected in Phase II IND protocols. We reported earlier statistically significant gender differences in PK parameters of INSCOP at a dose level of 0.4 mg. To identify covariates that influence PK parameters of INSCOP, we examined population covariates of INSCOP PK model for 0.4 mg dose. Methods: Plasma scopolamine concentrations versus time data were collected from 20 normal healthy human subjects (11 male/9 female) after a 0.4 mg dose. Phoenix NLME was employed for PK analysis of these data using gender, body weight and age as covariates for model selection. Model selection was based on a likelihood ratio test on the difference of criteria (-2LL). Statistical significance for base model building and individual covariate analysis was set at P less than 0.05{delta(-2LL)=3.84}. Results: A one-compartment pharmacokinetic model with first-order elimination best described INSCOP concentration ]time profiles. Inclusion of gender, body weight and age as covariates individually significantly reduced -2LL by the cut-off value of 3.84(P less than 0.05) when tested against the base model. After the forward stepwise selection and backward elimination steps, gender was selected to add to the final model which had significant influence on absorption rate constant (ka) and the volume of distribution (V) of INSCOP. Conclusion: A population pharmacokinetic model for INSCOP has been identified and gender was a significant contributing covariate for the final model. The volume of distribution and Ka were significantly higher in males than in females which confirm gender-dependent pharmacokinetics of scopolamine after administration of a 0.4 mg dose.

Wu, L.↗

Unsupervised segmentation of polarimetric SAR data using the covariance matrix

An unsupervised selection of polarimetric features useful for the segmentation and analysis of polarimetric synthetic aperture radar (SAR) data is presented. The technique is based on multidimensional clustering of the parameters composing the polarimetric covariance matrix of the data. Clustering is performed on the logarithm of these quantities. Once the polarimetric cluster centers have been determined, segmentation of the polarimetric data into regions is performed using a maximum likelihood polarimetric classifier. Segmentation maps are further improved using a Markov random field to describe the statistics of the regions and computing the maximum of the product of the local conditional densities. Examples with real polarimetric SAR imagery are given to illustrate the potential of this method.

Rignot, Eric↗

Assimilation of satellite microwave observations over the rainbands of tropical cyclones

A novel Bayesian Monte Carlo integration (BMCI) technique was developed to retrieve geophysical variables from satellite microwave radiometer data in the presence of tropical cyclones. The BMCI technique includes three steps: generating a stochastic database, simulating satellite brightness temperatures using a radiative transfer model, and retrieving geophysical variables such as profiles of temperature, relative humidity, and cloud liquid and ice water content from real observations. The technique also provides uncertainty estimates for each retrieval and can output the error covariance matrix of selected parameters. The measurements from the Advanced Technology Microwave Sounder (ATMS) on board Suomi National Polar-Orbiting Partnership (Suomi NPP) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI) were used as input. A new technique was developed to correct the ATMS and GMI observations for the beam-filling effect, which is due to small-scale variability of precipitation and clouds when compared with the instrument footprint and also the nonlinear relation between the brightness temperature and precipitation. In addition, the assimilation of the BMCI retrievals into the NASA GEOS model is discussed for Hurricane Maria. The results show that assimilating the BMCI retrievals can influence the dynamical features of the cyclone, including a stronger warm core, a symmetric eye, and vertically aligned wind columns. Two possible factors that may limit the impact of the BMCI retrievals include 1) the resolution of the model (about 25 km), which was too coarse to show the potential of the BMCI data in improving the representation of tropical storms in the model forecast, and 2) the data assimilation system not being able to consider vertically correlated observation errors.

Isaac Moradi↗

Evaluation of orbits with incomplete knowledge of the mathematical expectancy and the matrix of covariation of errors

The problem of selecting the optimal algorithm of filtration and the optimal composition of the measurements is examined assuming that the precise values of the mathematical expectancy and the matrix of covariation of errors are unknown. It is demonstrated that the optimal algorithm of filtration may be utilized for making some parameters more precise (for example, the parameters of the gravitational fields) after preliminary determination of the elements of the orbit by a simpler method of processing (for example, the method of least squares).

Bakhshiyan, B. T.↗

Specialization of the DNA-Cleaving Activity of a Group I Ribozyme Through In Vitro Evolution

In an earlier study, an in vitro evolution procedure was applied to a large population of variants of the Tetrahymena group 1 ribozyme to obtain individuals with a 10(exp 5)-fold improved ability to cleave a target single-stranded DNA substrate under simulated physiological conditions. The evolved ribozymes also showed a twofold improvement, compared to the wild-type, in their ability to cleave a single-stranded RNA substrate. Here, we report continuation of the in vitro evolution process using a new selection strategy to achieve both enhanced DNA and diminished RNA-cleavage activity. Our strategy combines a positive selection for DNA cleavage with a negative selection against RNA binding. After 36 "generations" of in vitro evolution, the evolved population showed an approx. 100-fold increase in the ratio of DNA to RNA-cleavage activity. Site-directed mutagenesis experiment confirmed the selective advantage of two covarying mutations within the catalytic core of ribozyme that are largely responsible for this modified behavior. The population of ribozymes has now undergone a total of 63 successive generations of evolution, resulting in an average 28 mutations relative to the wild-type that are responsible for the altered phenotype.

Tsang, Joyce↗

Can Assimilation of Satellite Ozone Data Contribute to the Understanding of the Lower Stratospheric Ozone?

We study the quality of lower stratospheric ozone fields from a three- dimensional global ozone assimilation system. Ozone in this region is important for the forcing of climate, but its global distribution is not fully known because of its large temporal and vertical variability. Modeled fields often have biases due to the inaccurate representation of transport processes in this region with strong gradients. Accurate ozonesonde or satellite occultation measurements have very limited coverage. Nadir measurements, such as those from the Solar Backscatter Ultraviolet/2 (SBUV/2) instrument that provide wide latitudinal coverage, lack the vertical resolution needed to represent sharp vertical features. Limb measurements, such as those from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), provide a finer vertical resolution. We show that assimilation of MIPAS data in addition to SBUV/2 data leads to better estimates of ozone in comparison with independent high quality satellite, aircraft, and ozone sonde measurements. Other modifications to the statistical analysis that have an impact on the lower stratospheric ozone will be mentioned: error covariance modeling and data selection. Direct and indirect impacts of transport and chemistry models will be discussed. Implications for multi-year analyses and short-tern prediction will be addressed.

Stajner, I.↗

Users manual for flight control design programs

Computer programs for the design of analog and digital flight control systems are documented. The program DIGADAPT uses linear-quadratic-gaussian synthesis algorithms in the design of command response controllers and state estimators, and it applies covariance propagation analysis to the selection of sampling intervals for digital systems. Program SCHED executes correlation and regression analyses for the development of gain and trim schedules to be used in open-loop explicit-adaptive control laws. A linear-time-varying simulation of aircraft motions is provided by the program TVHIS, which includes guidance and control logic, as well as models for control actuator dynamics. The programs are coded in FORTRAN and are compiled and executed on both IBM and CDC computers.

Nalbandian, J. Y.↗

Sensitivity Studies for Assimilated Ozone Profiles

An ozone data assimilation system at the NASA/Goddard Data Assimilation Office (DAO) produces three-dimensional global ozone fields. They are obtained by assimilating ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line parameterized chemistry and transport model. In this talk we focus on the quality of lower stratospheric assimilated ozone profiles. Ozone in the lower stratosphere plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. The SBUV/2 ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the lower stratosphere are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. We investigate the sensitivity of assimilated ozone profiles to changes in a variety of system inputs: TOMS and SBUV/2 data selection, forecast and observations error covariance models, inclusion or omission of a parameterized chemistry model, and different versions of DAO assimilated wind fields used to drive the transport model. Comparisons of assimilated ozone fields with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes.

Stajner, Ivanka↗

Assimilation of Satellite Ozone Observations

This talk will discuss assimilation of ozone data from satellite-borne instruments. Satellite observations of ozone total columns and profiles have been measured by a series of Total Ozone Mapping Spectrometer (TOMS), Solar Backscatter Ultraviolet (SBUV) instruments, and more recently by the Global Ozone Monitoring Experiment. Additional profile data are provided by instruments on NASA's Upper Atmosphere Research Satellite and by occultation instruments on other platforms. Instruments on Envisat' and future EOS Aura satellite will supply even more comprehensive data about the ozone distribution. Satellite data contain a wealth of information, but they do not provide synoptic global maps of ozone fields. These maps can be obtained through assimilation of satellite data into global chemistry and transport models. In the ozone system at NASA's Data Assimilation Office (DAO) any combination of TOMS, SBUV, and Microwave Limb sounder (MLS) data can be assimilated. We found that the addition of MLS to SBUV and TOMS data in the system helps to constrain the ozone distribution, especially in the polar night region and in the tropics. The assimilated ozone distribution in the troposphere and lower stratosphere is sensitive also to finer changes in the SBUV and TOMS data selection and to changes in error covariance models. All results are established by comparisons of assimilated ozone with independent profiles from ozone sondes and occultation instruments.

Stajner, I.↗