Search NASA⌕ Search

SEARCH · Search NASA

Results for “a priori estimates”

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 109 records · Page 6

Optimizing GRACE/GRACE-FO Data and A priori Hydrological Knowledge for Improved Global Terrestial Water Storage Component Estimates

The comprehensive information of global terrestrial water storage (TWS) components (soil moisture, groundwater, snow, surface water) is essential for effective assessment of water resource availability, climate variation, and disaster mitigation measures. Observational data provided by the Gravity Recovery And Climate Experiment (GRACE) and GRACE Follow-On satellite missions offer global TWS variation (ΔTWS) in terms of an integrated water column. However, GRACE spatial resolution is relatively coarse (i.e., 3°), and the vertically integrated value cannot be separated into ΔTWS components directly. This study demonstrates the feasibility to estimate ΔTWS components at any desired spatial-vertical resolution by effectively maintaining the native resolution of the employed hydrological knowledge. It utilizes a least-squares with constraints (LSC) approach to rigorously incorporate GRACE and GRACE-FO data and a priori hydrological knowledge, with the aim to improve global ΔTWS components’ accuracy and spatial resolution. The 3°×3° GRACE mascon derived ΔTWS data is disaggregated into the 0.5°×0.5° anomalous soil moisture storage (ΔSMS), groundwater storage (ΔGWS), snow water equivalent (ΔSWE), and surface water storage (ΔSWS) based on the covariance information obtained from the Community Atmosphere Biosphere Land Exchange (CABLE) and the PCRaster Global Water Balance (PCR-GLOBWB) models. Evaluation with different ground measurements and satellite products between 2002 and 2019 exhibits significantly improved accuracy in all individual ΔTWS components. This improvement is of particular note in ΔGWS and ΔSWS, where the LSC approach increases the globally averaged correlation values by approximately 0.13 and 0.05, respectively. Reliable prior knowledge leads to a more accurate ΔTWS component estimate, and the use of ensemble-mean knowledge yields the best result.

GRACE↗

Towards a Satellite Formaldehyde – in situ Hybrid Estimate for Organic Aerosol Abundance

Organic aerosol (OA) is one of the main components of the global particulate burden and intimately links natural and anthropogenic emissions with air quality and climate. It is challenging to accurately represent OA in global models. Direct quantification of global OA abundance is not possible with current remote sensing technology; however, it may be possible to exploit correlations of OA with remotely observable quantities to infer OA spatiotemporal distributions. In particular, formaldehyde (HCHO) and OA share common sources via both primary emissions and secondary production from oxidation of volatile organic compounds (VOCs). Here, we examine OA–HCHO correlations using data from summertime airborne campaigns investigating biogenic (NASA SEAC4RS and DC3), biomass burning (NASA SEAC4RS), and anthropogenic conditions (NOAA CalNex and NASA KORUS-AQ). In situ OA correlates well with HCHO (r=0.59–0.97), and the slope and intercept of this relationship depend on the chemical regime. For biogenic and anthropogenic regions, the OA–HCHO slopes are higher in low NOx conditions, because HCHO yields are lower and aerosol yields are likely higher. The OA–HCHO slope of wildfires is over 9 times higher than that for biogenic and anthropogenic sources. The OA–HCHO slope is higher for highly polluted anthropogenic sources (e.g., KORUS-AQ) than less polluted (e.g., CalNex) anthropogenic sources. Near-surface OAs over the continental US are estimated by combining the observed in situ relationships with HCHO column retrievals from NASA's Ozone Monitoring Instrument (OMI). HCHO vertical profiles used in OA estimates are from climatology a priori profiles in the OMI HCHO retrieval or output of specific period from a newer version of GEOS-Chem. Our OA estimates compare well with US EPA IMPROVE data obtained over summer months (e.g., slope =0.60–0.62, r=0.56 for August 2013), with correlation performance comparable to intensively validated GEOS-Chem (e.g., slope =0.57, r=0.56) with IMPROVE OA and superior to the satellite-derived total aerosol extinction (r=0.41) with IMPROVE OA. This indicates that OA estimates are not very sensitive to these HCHO vertical profiles and that a priori profiles from OMI HCHO retrieval have a similar performance to that of the newer model version in estimating OA. Improving the detection limit of satellite HCHO and expanding in situ airborne HCHO and OA coverage in future missions will improve the quality and spatiotemporal coverage of our OA estimates, potentially enabling constraints on global OA distribution.

anthropogenic emissions↗

Attitude determination error analysis - General model and specific application

This paper presents a comprehensive approach to filter and dynamics modeling for attitude determination error analysis. The discussion includes models for both batch least-squares and sequential estimators, a specific dynamic model for attitude determination error analysis of a three-axis stabilized spacecraft equipped with strapdown gyros, and the incorporation of general attitude sensor observations. An analyst using this approach to perform an error analysis chooses a subset of the spacecraft parameters to be 'solve-for' parameters, which are to be estimated, and another subset to be 'consider' parameters, which are assumed to have errors but not to be estimated. The result of the error analysis is an indication of overall uncertainties in the 'solve-for' parameters, as well as the contributions of the various error sources to these uncertainties, including those of errors in the a priori 'solve-for' estimates, of measurement noise, of dynamic noise (also known as process noise or plant noise), and of 'consider' parameter uncertainties. The analysis of attitude, star tracker alignment, and gyro bias uncertainties for the Gamma Ray Observatory spacecraft provide a specific example of the use of a general-purpose software package incorporating these models.

Markley, F. Landis↗

Velocimeter LIDAR Based Bulk Velocity Estimation for Terrain Relative Navigation Applications

A novel batch state estimation approach to estimate the translational and angular velocity states of a vehicle using a state-of-the-art velocimeter Light Detection and Ranging (LIDAR) sensor for use in Terrain and Hazard Relative Navigation (TRN/HRN) applications is presented in this paper. The velocimeter LIDAR is capable of measuring three dimensional position and line-of-sight (LOS) velocity associated with every pixel in the field-of-view (FOV). This new batch state estimation methodology is shown to provide accurate and statistically consistent velocity estimates with no a priori information. In contrast to traditional computer vision approaches for TRN, the proposed technique is not dependent on image features. This feature alleviates the need for accurate feature detection and correspondence to a predefined map making it suitable for unknown operating environments. Following a detailed development of the mathematical details associated with the batch state estimation methodology, the efficacy and utility of the proposed algorithms are evaluated through emulation robotics experiments performed at Texas A&M’s Land, Air, and Space Robotics (LASR) laboratory.

Skulsky, Eli↗

Asteroid orbital error analysis: Theory and application

We present a rigorous Bayesian theory for asteroid orbital error estimation in which the probability density of the orbital elements is derived from the noise statistics of the observations. For Gaussian noise in a linearized approximation the probability density is also Gaussian, and the errors of the orbital elements at a given epoch are fully described by the covariance matrix. The law of error propagation can then be applied to calculate past and future positional uncertainty ellipsoids (Cappellari et al. 1976, Yeomans et al. 1987, Whipple et al. 1991). To our knowledge, this is the first time a Bayesian approach has been formulated for orbital element estimation. In contrast to the classical Fisherian school of statistics, the Bayesian school allows a priori information to be formally present in the final estimation. However, Bayesian estimation does give the same results as Fisherian estimation when no priori information is assumed (Lehtinen 1988, and reference therein).

Muinonen, K.↗

LFSPMC: Linear feature selection program using the probability of misclassification

The computational procedure and associated computer program for a linear feature selection technique are presented. The technique assumes that: a finite number, m, of classes exists; each class is described by an n-dimensional multivariate normal density function of its measurement vectors; the mean vector and covariance matrix for each density function are known (or can be estimated); and the a priori probability for each class is known. The technique produces a single linear combination of the original measurements which minimizes the one-dimensional probability of misclassification defined by the transformed densities.

Guseman, L. F., Jr.↗

Knowledge-based vision for space station object motion detection, recognition, and tracking

Computer vision, especially color image analysis and understanding, has much to offer in the area of the automation of Space Station tasks such as construction, satellite servicing, rendezvous and proximity operations, inspection, experiment monitoring, data management and training. Knowledge-based techniques improve the performance of vision algorithms for unstructured environments because of their ability to deal with imprecise a priori information or inaccurately estimated feature data and still produce useful results. Conventional techniques using statistical and purely model-based approaches lack flexibility in dealing with the variabilities anticipated in the unstructured viewing environment of space. Algorithms developed under NASA sponsorship for Space Station applications to demonstrate the value of a hypothesized architecture for a Video Image Processor (VIP) are presented. Approaches to the enhancement of the performance of these algorithms with knowledge-based techniques and the potential for deployment of highly-parallel multi-processor systems for these algorithms are discussed.

Symosek, P.↗

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗

Synergistic Retrievals of Ice in High Clouds From Elastic Backscatter Lidar, Ku-band Radar and Submillimeter Wave Radiometer Observations

In this study, we investigate the synergy of elastic backscatter lidar, Ku-band radar, and sub-millimeter-wave radiometer measurements in the retrieval of ice from satellite observations. The synergy is analyzed through the generation of a large dataset of IceWater Content (IWC) profiles and simulated lidar, radar and radiometer observations. The characteristics of the instruments e.g. frequencies, sensitivities, etc. are set based on the expected characteristics of instruments of the Atmosphere Observing System (AOS) mission. A hold-out validation methodology is used to assess the accuracy of the IWC profiles retrieved from various combinations of observations from the three instruments. Specifically, the IWC and associated observations are randomly divided into two datasets, one for training and the other for evaluation. The training dataset is used to train the retrieval algorithm, while the evaluation dataset is used to assess the retrieval performance. The dataset of IWC profiles is derived from CloudSat reflectivity and CALIOP lidar observations. The retrieval of the ice water content IWC profiles from the computed observations is achieved in two steps. In the first step, a class, out of 18 potential classes characterized by different vertical distribution of IWC, is estimated from the observations. The 18 classes are predetermined based on the k-Means clustering algorithm. In the second step, the IWC profile is estimated using an Ensemble Kalman Smoother (EKS) algorithm that uses the estimated class as a priori information. The results of the study show that the synergy of lidar, radar, and radiometer observations is significant in the retrieval of the IWC profiles. Nevertheless, it should be mentioned that this synergy was found under idealized conditions, and additional work might be required to materialize it in practice. The inclusion of the lidar backscatter observations in the retrieval process has a larger impact on the retrieval performance than the inclusion of the radar observations. As ice clouds have a significant impact on atmospheric radiative processes, this work is relevant to ongoing efforts to reduce uncertainties in climate analyses and projections.

Mircea Grecu↗

The Impact of Prior Biosphere Models in the Inversion of Global Terrestrial CO2 Fluxes by Assimilating OCO-2 Retrievals

Atmospheric mixing ratios of carbon dioxide (CO2) are largely controlled by anthropogenic emissions and biospheric fluxes. The processes controlling terrestrial biosphere-atmosphere carbon exchange are currently not fully understood, resulting in terrestrial biospheric models having significant differences in the quantification of biospheric CO2 fluxes. Atmospheric transport models assimilating measured (in situ or space-borne) CO2 concentrations to estimate "top-down" fluxes, generally use these biospheric CO2 fluxes as a priori information. Most of the flux inversion estimates result in substantially different spatio-temporal posteriori estimates of regional and global biospheric CO2 fluxes. The Orbiting Carbon Observatory 2 (OCO-2) satellite mission dedicated to accurately measure column CO2 (XCO2) allows for an improved understanding of global biospheric CO2 fluxes. OCO-2 provides much-needed CO2 observations in data-limited regions facilitating better global and regional estimates of "top-down" CO2 fluxes through inversion model simulations. The specific objectives of our research are to: 1) conduct GEOS-Chem 4D-Var assimilation of OCO-2 observations, using several state-of-the-science biospheric CO2 flux models as a priori information, to better constrain terrestrial CO2 fluxes, and 2) quantify the impact of different biospheric model prior fluxes on OCO-2-assimilated a posteriori CO2 flux estimates. Here we present our assessment of the importance of these a priori fluxes by conducting Observing System Simulation Experiments (OSSE) using simulated OCO-2 observations with known "true" fluxes.

Biosphere↗

Polarized-interferometer feasibility study

The feasibility of using a polarized-interferometer system as a rendezvous and docking sensor for two cooperating spacecraft was studied. The polarized interferometer is a radio frequency system for long range, real time determination of relative position and attitude. Range is determined by round trip signal timing. Direction is determined by radio interferometry. Relative roll is determined from signal polarization. Each spacecraft is equipped with a transponder and an antenna array. The antenna arrays consist of four crossed dipoles that can transmit or receive either circularly or linearly polarized signals. The active spacecraft is equipped with a sophisticated transponder and makes all measurements. The transponder on the passive spacecraft is a relatively simple repeater. An initialization algorithm is developed to estimate position and attitude without any a priori information. A tracking algorithm based upon minimum variance linear estimators is also developed. Techniques to simplify the transponder on the passive spacecraft are investigated and a suitable configuration is determined. A multiple carrier CW signal format is selected. The dependence of range accuracy and ambiguity resolution error probability are derived and used to design a candidate system. The validity of the design and the feasibility of the polarized interferometer concept are verified by simulation.

Raab, F. H.↗

Estimation of the probability of error without ground truth and known a priori probabilities

The probability of error or, alternatively, the probability of correct classification (PCC) is an important criterion in analyzing the performance of a classifier. Labeled samples (those with ground truth) are usually employed to evaluate the performance of a classifier. Occasionally, the numbers of labeled samples are inadequate, or no labeled samples are available to evaluate a classifier's performance; for example, when crop signatures from one area from which ground truth is available are used to classify another area from which no ground truth is available. This paper reports the results of an experiment to estimate the probability of error using unlabeled test samples (i.e., without the aid of ground truth).

Havens, K. A.↗

Venus gravity field - Pioneer Venus Orbiter navigation results

The gravity field of Venus has been modeled by a spherical harmonic expansion of the potential to degree and order seven. The estimates of these coeficients were obtained by combining information from 43 short arcs (4 hr) of line-of-sight Doppler data centered at periapsis. The data arcs were distributed in longitude and time over more than two circulations of Venus by the Pioneer Venus Orbiter subperiapsis point which was confined to the band of latitudes from 14 deg N to 17 deg N. Convergence of the solution has been assured by iterating upon the initial estimate. All estimates were performed with zero a priori information on the gravity coefficients. Since the altitude of periapsis for most of the orbits was within the sensible Venusian atmosphere, drag effects on the estimated harmonics have been removed using an exponential atmosphere density model. Estimates of the mass parameter (GM) of Venus using this dataset are also evaluated.

Williams, B. G.↗

Trajectory Reconstruction and Uncertainty Analysis Using Mars Science Laboratory Pre-Flight Scale Model Aeroballistic Testing

As part of the Mars Science Laboratory (MSL) trajectory reconstruction effort at NASA Langley Research Center, free-flight aeroballistic experiments of instrumented MSL scale models was conducted at Aberdeen Proving Ground in Maryland. The models carried an inertial measurement unit (IMU) and a flush air data system (FADS) similar to the MSL Entry Atmospheric Data System (MEADS) that provided data types similar to those from the MSL entry. Multiple sources of redundant data were available, including tracking radar and on-board magnetometers. These experimental data enabled the testing and validation of the various tools and methodologies that will be used for MSL trajectory reconstruction. The aerodynamic parameters Mach number, angle of attack, and sideslip angle were estimated using minimum variance with a priori to combine the pressure data and pre-flight computational fluid dynamics (CFD) data. Both linear and non-linear pressure model terms were also estimated for each pressure transducer as a measure of the errors introduced by CFD and transducer calibration. Parameter uncertainties were estimated using a "consider parameters" approach.

Lugo, Rafael A.↗

Bounding filter - A simple solution to lack of exact a priori statistics.

Wiener and Kalman-Bucy estimation problems assume that models describing the signal and noise stochastic processes are exactly known. When this modeling information, i.e., the signal and noise spectral densities for Wiener filter and the signal and noise dynamic system and disturbing noise representations for Kalman-Bucy filtering, is inexactly known, then the filter's performance is suboptimal and may even exhibit apparent divergence. In this paper a system is designed whereby the actual estimation error covariance is bounded by the covariance calculated by the estimator. Therefore, the estimator obtains a bound on the actual error covariance which is not available, and also prevents its apparent divergence.

Nahi, N. E.↗

Stress distribution and topography of Tellus Regio, Venus

The Tellus Regio area of Venus represents a subset of a narrow latitude band where Pioneer Venus Orbiter (PVO) altimetry data, line-of-sight (LOS) gravity data, and Venera 15/16 radar images have all been obtained with good resolution. Tellus Regio also has a wide variety of surface morphologic features, elevations ranging up to 2.5 km, and a relatively low LOS gravity anomaly. This area was therefore chosen in order to examine the theoretical stress distributions resulting from various models of compensation of the observed topography. These surface stress distributions are then compared with the surface morphology revealed in the Venera 15/16 radar images. Conclusions drawn from these comparisons will enable constraints to be put on various tectonic parameters relevant to Tellus Regio. The stress distribution is calculated as a function of the topography, the equipotential anomaly, and the assumed model parameters. The topography data is obtained from the PVO altimetry. The equipotential anomaly is estimated from the PVO LOS gravity data. The PVO LOS gravity represents the spacecraft accelerations due to mass anomalies within the planet. These accelerations are measured at various altitudes and angles to the local vertical and therefore do not lend themselves to a straightforward conversion. A minimum variance estimator of the LOS gravity data is calculated, taking into account the various spacecraft altitudes and LOS angles and using the measured PVO topography as an a priori constraint. This results in an estimated equivalent surface mass distribution, from which the equipotential anomaly is determined.

Williams, David R.↗

A fast and accurate domain decomposition nonlinear manifold reduced order model

Here, this paper integrates nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD). NM ROMs approximate the full order model (FOM) state in a nonlinear-manifold by training a shallow, sparse autoencoder using FOM snapshot data. These NM-ROMs can be advantageous over linear-subspace ROMs (LS-ROMs) for problems with slowly decaying Kolmogorov n-width. However, the number of NM-ROM parameters that need to be trained scales with the size of the FOM. Moreover, for “extreme-scale” problems, the storage of high-dimensional FOM snapshots alone can make ROM training expensive. To alleviate the training cost, this paper applies DD to the FOM, computes NM-ROMs on each subdomain, and couples them to obtain a global NM-ROM. This approach has several advantages: Subdomain NM-ROMs can be trained in parallel, involve fewer parameters to be trained than global NM-ROMs, require smaller subdomain FOM dimensional training data, and can be tailored to subdomain specific features of the FOM. The shallow, sparse architecture of the autoencoder used in each subdomain NM-ROM allows application of hyper-reduction (HR), reducing the complexity caused by nonlinearity and yielding computational speedup of the NM-ROM. This paper provides the first application of NM-ROM (with HR) to a DD problem. In particular, this paper details an algebraic DD reformulation of the FOM, training a NM-ROM with HR for each sub domain, and a sequential quadratic programming (SQP) solver to evaluate the coupled global NM-ROM. Theoretical convergence results for the SQP method and a priori and a posteriori error estimates for the DD NM-ROM with HR are provided. The proposed DD NM-ROM with HR approach is numerically compared to a DD LS-ROM with HR on the 2D steady-state Burgers’ equation, showing an order of magnitude improvement in accuracy of the proposed DD NM-ROM over the DD LS-ROM.

97 MATHEMATICS AND COMPUTING↗