Search NASA⌕ Search

SEARCH · Search NASA

Results for “covariance analysis”

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 613 records · Page 34

Conjugate point determination for multitemporal data overlay

The machine processing is discussed of spatially variant multitemporal data such as imagery obtained at different times which requires that these data be in geometrical registration so that the analysis processor may obtain the datum for a specified ground resolution element in each of the sets of imagery being utilized for analysis. Misregistration between corresponding subsets of imagery contains both a displacement and a geometrical distortion component. Search techniques utilizing the moduli of the Fourier Transforms are developed for estimating the coefficients of geometrical distortion components of this model. Following the correction of these distortion components, the displacement is located by the crosscorrelation of a template obtained from one set of data, termed the reference, with the second, or background data. This template, derived for the optimum discrimination of the reference data embedded in the background, is determined by the solution of a system of equations involving the reference data and the covariance matrix of these data.

Emmert, R. A.↗

Scale covariance and G-varying cosmology. III - The /m, z/, /theta sub m, z/, /theta sub i, z/, and /N/m/, m/ tests

The Einstein gravitational equations in atomic units are fully solved for a matter-dominated universe in the context of a recently proposed scale-covariant cosmology. The magnitude-redshift relation for elliptical galaxies is studied, the evolutionary parameter used in such a study is derived, and the relation between isophotal angular diameters and redshifts is investigated, along with the relation between metric angular diameters and redshifts, the N(m)-magnitude relation for QSOs, and the magnitude-redshift relation for QSOs. Results are presented for four gauges (i.e., relations between G and the scale function beta (t)), and no contradictions are found between the proposed theory and the observational data. It is shown that only an open universe can fit the data if certain gauges suggested by a recent analysis of the time variation of the moon's period are selected and that observations made with atomic instruments do not necessarily yield geometrical parameters unless specific assumptions are made regarding the relation between atomic and gravitational dynamics.

Canuto, V. M.↗

Comparison of TOPEX/Poseidon orbit determination solutions obtained by the Goddard Space Flight Center Flight Dynamics Division and Precision Orbit Determination Teams

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using the Goddard Trajectory Determination System (GTDS) and a real-time extended Kalman filter estimation system to process Tracking Data and Relay Satellite (TDRS) System (TDRSS) measurements in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. GTDS is the operational orbit determination system used by the FDD, and the extended Kalman fliter was implemented in an analysis prototype system, the Real-Time Orbit Determination System/Enhanced (RTOD/E). The Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generates an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the Geodynamics (GEODYN) orbit determination system with laser ranging tracking data. The TOPEX/Poseidon trajectories were estimated for the October 22 - November 1, 1992, timeframe, for which the latest preliminary POD results were available. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch cases were assessed using overlap comparisons, while the sequential cases were assessed with covariances and the first measurement residuals. The batch least-squares and forward-filtered RTOD/E orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 10 meters (m) for the batch least squares and less than 18 m for the sequential estimation solutions. The differences among the POD, GTDS, and RTOD/E solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C.↗

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.

Ramanathan, Anand K.↗

Information content of a single pass of phase-delay data from a short baseline connected element interferometer

An analytic development of the information array obtained with a single tracking pass of phase-delay measurements made from a short baseline interferometer is presented. Phase-delay observations can be made with great precision from two antennas using a single, common distributed frequency standard, hence the name connected element. With the information array, closed-form expressions are developed for the error covariance in declination and right ascension. These equations serve as useful tools for analyzing the relative merits of candidate station locations for connected element interferometry (CEI). The navigation performance of a short baseline interferometer located at the Deep Space Network's (DSN's) Goldstone intracomplex is compared with that which is presently achievable using Very Long Baseline Interferometry (VLBI) over intercontinental baselines. The performance of an intracomplex pair of short baselines formed by three stations is also investigated, along with the use of a single baseline in conjunction with conventional two-way Doppler data. The phase-delay measurement accuracy and data rate used in the analysis are based on the expected performance of an experimental connected element system presently under construction at Goldstone. The results indicate that the VLBI system that will be used during the Galileo mission can determine the declination and right ascension of a distant spacecraft to an accuracy of 20 to 25 nrad, while the CEI triad system and the combination of CEI-Doppler system are both capable of 30 to 70 nrad performance.

Thurman, S. W.↗

CDGPS-Based Relative Navigation for Multiple Spacecraft

This thesis investigates the use of Carrier-phase Differential GPS (CDGPS) in relative navigation filters for formation flying spacecraft. This work analyzes the relationship between the Extended Kalman Filter (EKF) design parameters and the resulting estimation accuracies, and in particular, the effect of the process and measurement noises on the semimajor axis error. This analysis clearly demonstrates that CDGPS-based relative navigation Kalman filters yield good estimation performance without satisfying the strong correlation property that previous work had associated with "good" navigation filters. Several examples are presented to show that the Kalman filter can be forced to create solutions with stronger correlations, but these always result in larger semimajor axis errors. These linear and nonlinear simulations also demonstrated the crucial role of the process noise in determining the semimajor axis knowledge. More sophisticated nonlinear models were included to reduce the propagation error in the estimator, but for long time steps and large separations, the EKF, which only uses a linearized covariance propagation, yielded very poor performance. In contrast, the CDGPS-based Unscented Kalman relative navigation Filter (UKF) handled the dynamic and measurement nonlinearities much better and yielded far superior performance than the EKF. The UKF produced good estimates for scenarios with long baselines and time steps for which the EKF would diverge rapidly. A hardware-in-the-loop testbed that is compatible with the Spirent Simulator at NASA GSFC was developed to provide a very flexible and robust capability for demonstrating CDGPS technologies in closed-loop. This extended previous work to implement the decentralized relative navigation algorithms in real time.

Mitchell, Megan Leigh↗

Dynamic Creep Buckling: Analysis of Shell Structures Subjected to Time-dependent Mechanical and Thermal Loading

The objective of the present research is to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. A complete, true ab-initio rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformations, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of metric and curvature in two coordinate systems: fixed (spatial) and convected (material). The relations between the time derivative and the covariant derivative (gradient) were developed for curved space and motion, so the velocity components supply the connection between the equations of motion and the time rates of change of the metric and curvature tensors.

Simitses, G. J.↗

Evaluating A Priori Ozone Profile Information Used in TEMPO Tropospheric Ozone Retrievals

Ozone (O3) is a greenhouse gas and toxic pollutant which plays a major role in air quality. Typically, monitoring of surface air quality and O3 mixing ratios is primarily conducted using in situ measurement networks. This is partially due to high-quality information related to air quality being limited from space-borne platforms due to coarse spatial resolution, limited temporal frequency, and minimal sensitivity to lower tropospheric and surface-level O3. The Tropospheric Emissions: Monitoring of Pollution (TEMPO) satellite is designed to address these limitations of current space-based platforms and to improve our ability to monitor North American air quality. TEMPO will provide hourly data of total column and vertical profiles of O3 with high spatial resolution to be used as a near-real-time air quality product. TEMPO O3 retrievals will apply the Smithsonian Astrophysical Observatory profile algorithm developed based on work from GOME, GOME-2, and OMI. This algorithm uses a priori O3 profile information from a climatological data-base developed from long-term ozone-sonde measurements (tropopause-based (TB) O3 climatology). It has been shown that satellite O3 retrievals are sensitive to a priori O3 profiles and covariance matrices. During this work we investigate the climatological data to be used in TEMPO algorithms (TB O3) and simulated data from the NASA GMAO Goddard Earth Observing System (GEOS-5) Forward Processing (FP) near-real-time (NRT) model products. These two data products will be evaluated with ground-based lidar data from the Tropospheric Ozone Lidar Network (TOLNet) at various locations of the US. This study evaluates the TB climatology, GEOS-5 climatology, and 3-hourly GEOS-5 data compared to lower tropospheric observations to demonstrate the accuracy of a priori information to potentially be used in TEMPO O3 algorithms. Here we present our initial analysis and the theoretical impact on TEMPO retrievals in the lower troposphere.

ozone↗

Evaluating a Priori Ozone Profile Information Used in TEMPO Tropospheric Ozone Retrievals

Ozone (O3) is a greenhouse gas and toxic pollutant which plays a major role in air quality. Typically, monitoring of surface air quality and O3 mixing ratios is primarily conducted using in situ measurement networks. This is partially due to high-quality information related to air quality being limited from space-borne platforms due to coarse spatial resolution, limited temporal frequency, and minimal sensitivity to lower tropospheric and surface-level O3. The Tropospheric Emissions: Monitoring of Pollution (TEMPO) satellite is designed to address these limitations of current space-based platforms and to improve our ability to monitor North American air quality. TEMPO will provide hourly data of total column and vertical profiles of O3 with high spatial resolution to be used as a near-real-time air quality product.TEMPO O3 retrievals will apply the Smithsonian Astrophysical Observatory profile algorithm developed based on work from GOME, GOME-2, and OMI. This algorithm uses a priori O3 profile information from a climatological data-base developed from long-term ozone-sonde measurements (tropopause-based (TB) O3 climatology). It has been shown that satellite O3 retrievals are sensitive to a priori O3 profiles and covariance matrices. During this work we investigate the climatological data to be used in TEMPO algorithms (TB O3) and simulated data from the NASA GMAO Goddard Earth Observing System (GEOS-5) Forward Processing (FP) near-real-time (NRT) model products. These two data products will be evaluated with ground-based lidar data from the Tropospheric Ozone Lidar Network (TOLNet) at various locations of the US. This study evaluates the TB climatology, GEOS-5 climatology, and 3-hourly GEOS-5 data compared to lower tropospheric observations to demonstrate the accuracy of a priori information to potentially be used in TEMPO O3 algorithms. Here we present our initial analysis and the theoretical impact on TEMPO retrievals in the lower troposphere.

TEMPO↗

Shortwave feedbacks and El Nino-Southern Oscillation: Forced ocean coupled ocean-atmosphere experiments

Changes in tropical sea surface temperature (SST) can produce changes in cloudiness that modify incoming solar shortwave (SW) radiation, which in turn affects SST. The effects of this negative feedback on Pacific interannual variability are examined in forced ocean model and hybrid coupled ocean-atmosphere model simulations. Two empirical schemes are used to model the large-scale, low-frequency response of surface SW to SST anomalies. The first scheme attempts to account for the nonlocal nature of the atmospheric response to SST based patterns of covariability analyzed through singular value decomposition. In the observations the primary coupled mode of variability is composed of a SW anomaly in the central Pacific that covaries with anomalous SST in the eastern Pacific. This is applied in the model as a nonlocal feedback. The second scheme examines the effects of a purely local feedback with a spatially varying coefficient of magnitude chosen similar to the first scheme. In almost all cases the second scheme behaved similarly to the first, presumably because the correlation scale of SST is large enough for El Nino-Southern Oscillation (ENSO) dynamics that there is little sensitivity to the local approximation in the SW feedback. In simulations forced by time series of observed wind stress the SW feedback induced very minor SST damping. Results for a simpified heat budget analysis showed that while the SW feedback increased the local heat flux damping on SST, it also induced a mean shallowing of the mixed layer. The resulting changes in both the local mean vertical temperature gradient and the zonal velocity response to the wind stress acted to oppose the local heat flux damping effects. When the observed SW anomalies were applied to forced simulations, the simulated SST anomalies were modified as expected, and agreement with observed SST improved. In coupled simulations the SW feedbacks had greater impact than in the case of specified stress. The main effects were to decrease the magnitude of the warm and cold SST anomalies in the central Pacific, while leaving the pattern and evolution of ENSO anomalies essentially unchanged elsewhere. The SW feedbacks thus produce a modest improvement of the model ENSO SST pattern compared with observations, although they tended to shorten the period of the model ENSO cycle. Overall the results suggest that large-scale SW feedbacks are of quantitative importance to simulating some aspects of the ENSO cycle but are not critical to the overall occurrence of the phenomenon.

Waliser, Duane E.↗

An analysis of the processing requirements of a complex perceptual-motor task

Current concerns in the assessment of mental workload are discussed, and the event-related brain potential (ERP) is introduced as a promising mental-workload index. Subjects participated in a series of studies in which they were required to perform a target acquisition task while also covertly counting either auditory or visual probes. The effects of several task-difficulty manipulations on the P300 component of the ERP elicited by the counted stimulus probes were investigated. With sufficiently practiced subjects the amplitude of the P300 was found to decrease with increases in task difficulty. The second experiment also provided evidence that the P300 is selectively sensitive to task-relevant attributes. A third experiment demonstrated a convergence in the amplitude of the P300s elicited in the simple and difficult versions of the tracking task. The amplitude of the P300 was also found to covary with the measures of tracking performance. The results of the series of three experiments illustrate the sensitivity of the P300 to the processing requirements of a complex target acquisition task. The findings are discussed in terms of the multidimensional nature of processing resources.

Kramer, A. F.↗

Estimating the instabilities of N clocks by means of comparison measurements

The estimation of individual instabilities of N clocks, compared by measuring the differences of their readings, is considered without assuming a priori any hypotheses on their uncorrelation. Instabilities of the N clocks are described by a complete (non-diagonal) N x N covariance matrix R. Only differences of clock readings are available in order to estimate R. Statistical processing of these data allows one to calculate the (N-1)x(N-l) covariance matrix S of the differences relative to the N-th(reference) clock. By analyzing the relationships tying R and S, several pieces of information can be inferred and, in particular, the conditions for the validity of the uncorrelation hypothesis are established. The estimation of R from S is not unique: in any case R must be positive definite. A theorem states that R is positive definite if and only if its determinant is positive. Nevertheless infinitely many acceptable choices of R still fulfill the condition of positive definiteness. This paper shows that, by increasing the number N of compared clocks, the amount of arbitrariness in estimating R is reduced. The analysis of some experimental data illustrates the capability of the method.

Premoli, Amedeo↗

A finite element technique for non-deterministic thermal deformation analyses including temperature dependent material properties

A technique utilizing the finite element displacement method is developed for the static analysis of structures subjected to non-deterministic thermal loading in which the material properties, assumed isotropic, are temperature dependent. Matrix equations are developed for the first two statistical moments of the displacements using a third order series expansion for the displacements in terms of the random temperatures. Sample problems are included to demonstrate the range of applicability of the third order series solutions. These solutions are compared with results from Monte Carlo analyses and also, for some problems, with solutions obtained by numerically integrating equations for the statistical properties of the displacements. In general, it is shown that the effect of temperature dependent material properties can have a significant effect on the covariances of the displacements.

Case, W. R., Jr.↗

LANDSAT-4 MSS and TM Spectral Class Comparison and Coherent Noise Analysis

A detailed spectral analysis is conducted of thematic mapper and MSS data for an area near Des Moines, Iowa. Data are utilized from 7 blocks distributed throughout the area which included agricultural, forest, suburan, urban, and water scene types. The blocks are processed using a clustering algorithm to produce up to 18 cluster groupings for each block. Each cluster class is then identified with a ground-cover class using aerial photography and maps of the area. The clusters from each of the 7 blocks are inspected with regard to separability, mean, and variances. The separability measure used in the transformed divergence function or processor measures the statistical distance between classes based on class means and covariance matrices. The measure has a maximum value of 2,000 and the minimum of 0. Spectrally, very close classes will typically have values as low as 50 to 500.

Anuta, P. E.↗

Geometrical approach to neural net control of movements and posture

In one approach to modeling brain function, sensorimotor integration is described as geometrical mapping among coordinates of non-orthogonal frames that are intrinsic to the system; in such a case sensors represent (covariant) afferents and motor effectors represent (contravariant) motor efferents. The neuronal networks that perform such a function are viewed as general tensor transformations among different expressions and metric tensors determining the geometry of neural functional spaces. Although the non-orthogonality of a coordinate system does not impose a specific geometry on the space, this "Tensor Network Theory of brain function" allows for the possibility that the geometry is non-Euclidean. It is suggested that investigation of the non-Euclidean nature of the geometry is the key to understanding brain function and to interpreting neuronal network function. This paper outlines three contemporary applications of such a theoretical modeling approach. The first is the analysis and interpretation of multi-electrode recordings. The internal geometries of neural networks controlling external behavior of the skeletomuscle system is experimentally determinable using such multi-unit recordings. The second application of this geometrical approach to brain theory is modeling the control of posture and movement. A preliminary simulation study has been conducted with the aim of understanding the control of balance in a standing human. The model appears to unify postural control strategies that have previously been considered to be independent of each other. Third, this paper emphasizes the importance of the geometrical approach for the design and fabrication of neurocomputers that could be used in functional neuromuscular stimulation (FNS) for replacing lost motor control.

Review↗

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗

Vertical Structure of a Springtime Smoky and Humid Troposphere Over the Southeast Atlantic From Aircraft and Reanalysis

The springtime atmosphere over the southeast Atlantic Ocean (SEA) is subjected to a consistent layer of biomass burning (BB) smoke from widespread fires on the African continent. An elevated humidity signal is coincident with this layer, consistently proportional to the amount of smoke present. The combined humidity and BB aerosol has potentially significant radiative and dynamic impacts. Here, we use aircraft-based observations from the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) deployments in conjunction with reanalyses to characterize covariations in humidity and BB smoke across the SEA. The observed plume–vapor relationship, and its agreement with the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) and Copernicus Atmosphere Monitoring Service (CAMS) reanalysis, persists across all observations, although the magnitude of the relationship varies as the season progresses. Water vapor is well represented by the reanalyses, while CAMS tends to underestimate carbon monoxide especially under high BB. While CAMS aerosol optical depth (AOD) is generally overestimated relative to ORACLES AOD, the observations show a consistent relationship between carbon monoxide (CO) and aerosol extinction, demonstrating the utility of the CO tracer to understanding vertical aerosol distribution. We next use k-means clustering of the reanalyses to examine multi-year seasonal patterns and distributions. We identify canonical profile types of humidity and of CO, allowing us to characterize changes in vapor and BB atmospheric structures, and their impacts as they covary. While the humidity profiles show a range in both total water vapor concentration and in vertical structure, the CO profiles primarily vary in terms of maximum concentration, with similar vertical structures in each. The distribution of profile types varies spatiotemporally across the SEA region and through the season, ranging from largely one type in the northeast and southwest to more evenly distributed between multiple types where air masses meet in the middle of the SEA. These distributions follow patterns of transport from the humid, smoky source region (greatest influence in the northeast of the SEA) and the seasonal changes in both humidity and smoke (increasing and decreasing through the season, respectively). With this work, we establish a framework for a more complete analysis of the broader radiative and dynamical effects of humid aerosols over the SEA.

Kristina Pistone↗

Improved solution accuracy for TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the operational orbit determination system used by the FDD in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. The extended Kalman filter was implemented in an orbit determination analysis prototype system, closely related to the Real-Time Orbit Determination System/Enhanced (RTOD/E) system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the geodynamics (GEODYN) orbit determination system with laser ranging and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) tracking measurements. The TOPEX/Poseidon trajectories were estimated for November 7 through November 11, 1992, the timeframe under study. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch-least-squares solutions were assessed based on the solution residuals, while the sequential solutions were assessed based on primarily the estimated covariances. The batch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 2 meters for the batch-least-squares and less than 13 meters for the sequential estimation solutions. After the sequential estimation solutions were processed with a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were obtained. The differences among the POD, GTDS, and filter/smoother solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C. E.↗