Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance 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 145 records · Page 8

Calibration of the COBE FIRAS instrument

The Far-Infrared Absolute Spectrophotometer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite was designed to accurately measure the spectrum of the cosmic microwave background radiation (CMBR) in the frequency range 1-95/cm with an angular resolution of 7 deg. We describe the calibration of this instrument, including the method of obtaining calibration data, reduction of data, the instrument model, fitting the model to the calibration data, and application of the resulting model solution to sky observations. The instrument model fits well for calibration data that resemble sky condition. The method of propagating detector noise through the calibration process to yield a covariance matrix of the calibrated sky data is described. The final uncertainties are variable both in frequency and position, but for a typical calibrated sky 2.6 deg square pixel and 0.7/cm spectral element the random detector noise limit is of order of a few times 10(exp -7) ergs/sq cm/s/sr cm for 2-20/cm, and the difference between the sky and the best-fit cosmic blackbody can be measured with a gain uncertainty of less than 3%.

Fixsen, D. J.↗

Sequential estimation and satellite data assimilation in meteorology and oceanography

The central theme of this review article is the role that dynamics plays in estimating the state of the atmosphere and of the ocean from incomplete and noisy data. Objective analysis and inverse methods represent an attempt at relying mostly on the data and minimizing the role of dynamics in the estimation. Four-dimensional data assimilation tries to balance properly the roles of dynamical and observational information. Sequential estimation is presented as the proper framework for understanding this balance, and the Kalman filter as the ideal, optimal procedure for data assimilation. The optimal filter computes forecast error covariances of a given atmospheric or oceanic model exactly, and hence data assimilation should be closely connected with predictability studies. This connection is described, and consequences drawn for currently active areas of the atmospheric and oceanic sciences, namely, mesoscale meteorology, medium and long-range forecasting, and upper-ocean dynamics.

Ghil, M.↗

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation that is based on a coupling between inertial navigation and flush air data sensing pressure measurements. In this approach, the full navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to directly estimate atmospheric winds and density using a nonlinear weighted least-squares algorithm. The approach uses a high fidelity model of atmosphere stored in table-look-up form, along with simplified models of that are propagated along the trajectory within the algorithm to provide prior estimates and covariances to aid the air data state solution. Thus, the method is essentially a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere and winds are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the discrete-time observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content to the system. The algorithm is then applied to the design of the pressure measurement system for the Mars 2020 mission. The pressure port layout is optimized to maximize the observability of atmospheric states along the trajectory. Linear covariance analysis is performed to assess estimator performance for a given pressure measurement uncertainty. The results indicate that the new tightly-coupled estimator can produce enhanced estimates of atmospheric states when compared with existing algorithms.

Karlgaard, Christopher D.↗

Investigating the Causes and Consequences of Controlled Rest on the Flight Deck

INTRODUCTION: Long and irregular working hours can lead to fatigue in aviation operations. In some regions, a short nap taken on the flight deck (known as controlled rest) can be used as a countermeasure to unexpected in-flight sleepiness. We aimed to investigate the impact of taking controlled rest on self-reported sleepiness at top-of-descent. METHODS: Data from 120 long-haul (> 6 h flight duration), unaugmented flights were analyzed (n = 31 pilots). Pilots wore actigraphs and completed sleep logs before and during trips. At pre-flight and top-of-descent, pilots completed a Karolinska Sleepiness Scale (KSS). A mixed-effects model was used to assess the impact of controlled rest on KSS at top-of-descent. Sleep duration in the 48 hours prior to departure, timing of the flight (day vs. night), and pre-flight KSS scores were included as covariates. RESULTS: Due to missing data, complete data from 83 flights (n = 29 participants) were available in the analyses of the KSS. There were no differences by controlled rest status for KSS scores at top-of-descent (estimated marginal means ± SEM with controlled rest: 5.10 ± 0.22, without controlled rest: 5.37 ± 0.29; p = .45, η2p = 0.01). DISCUSSION: Our results suggest that there is no difference in self-reported sleepiness at top-of-descent on flights in which controlled rest was taken compared to flights without controlled rest. Further research is necessary to determine the impact of controlled rest on objective measures of performance at top-of-descent.

aviation↗

The Impact of Controlled Rest on Self-Reported Sleepiness at Top-of-Descent

INTRODUCTION: Long and irregular working hours can lead to fatigue in aviation operations. In some regions, a short nap taken on the flight deck (known as controlled rest) can be used as a countermeasure to unexpected in-flight sleepiness. We aimed to investigate the impact of taking controlled rest on self-reported sleepiness at top-of-descent. METHODS: Data from 120 long-haul (> 6 h flight duration), unaugmented flights were analyzed (n = 31 pilots). Pilots wore actigraphs and completed sleep logs before and during trips. At pre-flight and top-of-descent, pilots completed a Karolinska Sleepiness Scale (KSS). A mixed-effects model was used to assess the impact of controlled rest on KSS at top-of-descent. Sleep duration in the 48 hours prior to departure, timing of the flight (day vs. night), and pre-flight KSS scores were included as covariates. RESULTS: Due to missing data, complete data from 83 flights (n = 29 participants) were available in the analyses of the KSS. There were no differences by controlled rest status for KSS scores at top-of-descent (estimated marginal means ± SEM with controlled rest: 5.10 ± 0.22, without controlled rest: 5.37 ± 0.29; p = .45, η2p = 0.01). DISCUSSION: Our results suggest that there is no difference in self-reported sleepiness at top-of-descent on flights in which controlled rest was taken compared to flights without controlled rest. Further research is necessary to determine the impact of controlled rest on objective measures of performance at top-of-descent.

aviation↗

The Impact of Controlled Rest on Neurobehavioral Outcomes at Top-of-Descent

INTRODUCTION: Long and irregular working hours can lead to fatigue in aviation operations. In some regions, a short nap taken on the flight deck (known as controlled rest) can be used as a countermeasure to unexpected in-flight sleepiness. We aimed to investigate the impact of taking controlled rest on neurobehavioral measures at top-of-descent. METHODS: Data from 120 long-haul (> 6 h flight duration), unaugmented flights were analyzed (n = 31 pilots). Pilots wore actigraphs and completed sleep logs before and during trips. At pre-flight and top-of-descent, pilots completed a 5-minute psychomotor vigilance task (PVT) and Karolinska Sleepiness Scale (KSS). A series of mixed-effects models were conducted to assess the impact of controlled rest on outcome measures at top-of-descent. Sleep duration in the prior 48 hours, timing of the flight, and pre-flight scores for each measure were included as covariates. RESULTS: Due to missing data, complete data from 76 flights (n = 28 participants) were available in the models examining the PVT metrics, and data from 83 flights (n = 29 participants) were available in the analyses of the KSS. Pilots who took controlled rest had faster response speeds at top-of-descent (p = .03, η2p = 0.07; estimated marginal mean [EMM] = 4.19, standard error [SE] = 0.07, 95% CI [4.08, 4.29]) than those who did not take controlled rest (EMM = 4.00, SE = 0.05, 95% CI [3.86, 4.14]). There were no differences by controlled rest status for KSS scores and PVT lapses (p values > .05, η2p values ≤ 0.01). DISCUSSION: Our results suggest that taking controlled rest may improve vigilant attention at critical phases of flight and, thus, may be a useful fatigue management tool during unaugmented flights. Further research is necessary to determine the impact of factors on the decision to take controlled rest (e.g., airline culture, personal preference) and how the controlled rest policy is applied in practice.

aviation↗

Orbit Determination of Spacecraft in Earth-Moon L1 and L2 Libration Point Orbits

The ARTEMIS mission, part of the THEMIS extended mission, is the first to fly spacecraft in the Earth-Moon Lissajous regions. In 2009, two of the five THEMIS spacecraft were redeployed from Earth-centered orbits to arrive in Earth-Moon Lissajous orbits in late 2010. Starting in August 2010, the ARTEMIS P1 spacecraft executed numerous stationkeeping maneuvers, initially maintaining a lunar L2 Lissajous orbit before transitioning into a lunar L1 orbit. The ARTEMIS P2 spacecraft entered a L1 Lissajous orbit in October 2010. In April 2011, both ARTEMIS spacecraft will suspend Lissajous stationkeeping and will be maneuvered into lunar orbits. The success of the ARTEMIS mission has allowed the science team to gather unprecedented magnetospheric measurements in the lunar Lissajous regions. In order to effectively perform lunar Lissajous stationkeeping maneuvers, the ARTEMIS operations team has provided orbit determination solutions with typical accuracies on the order of 0.1 km in position and 0.1 cm/s in velocity. The ARTEMIS team utilizes the Goddard Trajectory Determination System (GTDS), using a batch least squares method, to process range and Doppler tracking measurements from the NASA Deep Space Network (DSN), Berkeley Ground Station (BGS), Merritt Island (MILA) station, and United Space Network (USN). The team has also investigated processing of the same tracking data measurements using the Orbit Determination Tool Kit (ODTK) software, which uses an extended Kalman filter and recursive smoother to estimate the orbit. The orbit determination results from each of these methods will be presented and we will discuss the advantages and disadvantages associated with using each method in the lunar Lissajous regions. Orbit determination accuracy is dependent on both the quality and quantity of tracking measurements, fidelity of the orbit force models, and the estimation techniques used. Prior to Lissajous operations, the team determined the appropriate quantity of tracking measurements that would be needed to meet the required orbit determination accuracies. Analysts used the Orbit Determination Error Analysis System (ODEAS) to perform covariance analyses using various tracking data schedules. From this analysis, it was determined that 3.5 hours of DSN TRK-2-34 range and Doppler tracking data every other day would suffice to meet the predictive orbit knowledge accuracies in the Lissajous region. The results of this analysis are presented. Both GTDS and ODTK have high-fidelity environmental orbit force models that allow for very accurate orbit estimation in the lunar Lissajous regime. These models include solar radiation pressure, Earth and Moon gravity models, third body gravitational effects from the Sun, and to a lesser extent third body gravitational effects from Jupiter, Venus, Saturn, and Mars. Increased position and velocity uncertainties following each maneuver, due to small execution performance errors, requires that several days of post-maneuver tracking data be processed to converge on an accurate post-maneuver orbit solution. The effects of maneuvers on orbit determination accuracy will be presented, including a comparison of the batch least squares technique to the extended Kalman filter/smoother technique. We will present the maneuver calibration results derived from processing post-maneuver tracking data. A dominant error in the orbit estimation process is the uncertainty in solar radiation pressure and the resultant force on the spacecraft. An estimation of this value can include many related factors, such as the uncertainty in spacecraft reflectivity and surface area which is a function of spacecraft orientation (spin-axis attitude), uncertainty in spacecraft wet mass, and potential seasonal variability due to the changing direction of the Sun line relative to the Earth-Moon Lissajous reference frame. In addition, each spacecraft occasionally enters into Earth or Moon penumbra or umbra and these shadow crossings reduche solar radiation force for several hours. The effects of these events on orbit determination accuracy will be presented. In order to plan for upcoming stationkeeping maneuvers, the maneuver planning team must take the current orbit estimate, propagate it forward to the planned maneuver time, and determine the optimal maneuver to maintain the Lissajous orbit for one or more revolutions. The propagation is performed using a Runge-Kutta 7/8 integrator and typically the position and velocity uncertainty increases with propagation time, increasing the overall uncertainty of the orbit state at the maneuver execution time. The effect of orbit knowledge uncertainty on stationkeeping operations will be presented.

Woodard, Mark↗

Version 2 of the global catalogue of large anthropogenic and volcanic SO 2 sources and emissions derived from satellite measurements

Sulfur dioxide (SO 2 ) measurements from the Ozone Monitoring Instrument (OMI), Ozone Mapping and Profiler Suite (OMPS), and TROPOspheric Monitoring Instrument (TROPOMI) satellite spectrometers were used to update and extend the previously developed global catalogue of large SO 2 emission sources. This version 2 of the global catalogue covers the period of 2005–2021 and includes a total of 759 continuously emitting point sources releasing from about 10 kt yr −1 to more than 4000 kt yr −1 of SO 2 , that have been identified and grouped by country and primary source origin: volcanoes (106 sources); power plants (477); smelters (74); and sources related to the oil and gas industry (102). There are several major improvements compared to the original catalogue: it combines emissions estimates from three satellite instruments instead of just OMI, uses a new version 2 of the OMI and OMPS SO 2 dataset, and updated consistent site-specific air mass factors (AMFs) are used to calculate SO 2 vertical column densities (VCDs). The newest TROPOMI SO 2 data processed with the Covariance-Based Retrieval Algorithm (COBRA), used in the catalogue, can detect sources with emissions as low as 8 kt yr −1 (in 2018–2021) compared to the 30 kt yr −1 limit for OMI. In general, there is an overall agreement within ±12 % in total emissions estimated from the three satellite instruments for large regions. For individual emission sources, the spread is larger: the annual emissions estimated from OMI and TROPOMI agree within ±13 % in 50 % of cases and within ±28 % in 90 % of cases. The version 2 catalogue emissions were calculated as a weighted average of emission estimates from the three satellite instruments using an inverse-variance weighting method. OMI, OMPS, and TROPOMI data contribute 7 %, 5 %, and 88 % to the average, respectively, for small (<30 kt yr −1 ) sources and 33 %, 20 %, and 47 %, respectively, for large (>300 kt yr −1 ) sources. The catalogue data show an approximate 50 % decline in global SO 2 emissions between 2005 and 2021, although emissions were relatively stable during the last 3 years. The version 2 of the global catalogue has been posted at the NASA global SO 2 monitoring website (https://doi.org/10.5067/MEASURES/SO2/DATA406, Fioletov et al., 2022).

volcanic↗

The processing of the Viking Orbiter range data and its contribution to Mars gravity solutions

The processing of Doppler data has been the primary method for deriving models of the Mars gravity field. Since the Mariner 9 and Viking spacecraft were placed in orbit about Mars, many models from degree and order 6 to degree and order 50 have been developed. However, during the Viking mission, some 26,000 range measurements to the two Viking Orbiters were also obtained. These data have not previously been used in the derivation of Mars gravity models. A portion of these range data have been processed simultaneously with the Doppler data. Normal equations were generated for both sets of data and were used to create two solutions complete to degree and order 30: a nominal solution including both the range and the Doppler data (MGM-R100), and another solution including only the Doppler data (MGM-R101). Tests with the covariances of these solutions, as well as with orbit overlap tests indicate that the interplanetary range data can be used to improve the modeling of the Mars gravity field.

Lemoine, Frank G.↗

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗

On the computation and updating of the modified Cholesky decomposition of a covariance matrix

Methods for obtaining and updating the modified Cholesky decomposition (MCD) for the particular case of a covariance matrix when one is given only the original data are described. These methods are the standard method of forming the covariance matrix K then solving for the MCD, L and D (where K=LDLT); a method based on Householder reflections; and lastly, a method employing the composite-t algorithm. For many cases in the analysis of remotely sensed data, the composite-t method is the superior method despite the fact that it is the slowest one, since (1) the relative amount of time computing MCD's is often quite small, (2) the stability properties of it are the best of the three, and (3) it affords an efficient and numerically stable procedure for updating the MCD. The properties of these methods are discussed and FORTRAN programs implementing these algorithms are listed.

Vanrooy, D. L.↗

Systems Engineering Design Via Experimental Operation Research: Complex Organizational Metric for Programmatic Risk Environments (COMPRE)

Unique and innovative graph theory, neural network, organizational modeling, and genetic algorithms are applied to the design and evolution of programmatic and organizational architectures. Graph theory representations of programs and organizations increase modeling capabilities and flexibility, while illuminating preferable programmatic/organizational design features. Treating programs and organizations as neural networks results in better system synthesis, and more robust data modeling. Organizational modeling using covariance structures enhances the determination of organizational risk factors. Genetic algorithms improve programmatic evolution characteristics, while shedding light on rulebase requirements for achieving specified technological readiness levels, given budget and schedule resources. This program of research improves the robustness and verifiability of systems synthesis tools, including the Complex Organizational Metric for Programmatic Risk Environments (COMPRE).

Mog, Robert A.↗

Contribution of X/Ka VLBI to Multi-Wavelength Celestial Frame Studies

This paper is an update of Sotuela et al. (2011) which improves their simulated Gaia frame tie precision by approximately 10% by adding three additional VLBI observing sessions. Astrometry at X/Ka-band (8.4/32 GHz) using NASAs Deep Space Network has detected 466 quasars with accuracies of 200-300 micro-arc seconds. A program is underway to reduce errors by a factor of 2-3. From our sample, 245 sources have optical magnitudes V less than 20 and should also be detectable by Gaia. A covariance study using existing X/Ka data and simulated Gaia uncertainties for the 345 objects yields a frame tie precision of 10-15 micro-arc seconds (1 - sigma). The characterization of wavelength dependent systematic from extended source morphology and core shift should benefit greatly from adding X/Ka-band measurements to S/X-band (2.3/8.4 GHz) measurements thus helping to constrain astrophysical models of the wavelength dependence of positions.

galaxies↗

Nearest-Neighbor Machine Learning Feature Selection for Interpretation of Microbial Molecular Signatures from Isotope Ratio Mass Spectrometry Data

Mass spectrometry (MS) promises to be a powerful tool for potential biosignature detection during astrobiological missions on ocean worlds in our solar system. Accurate and generalizable machine learning methods could enhance science return on investment by predicting seawater chemistry and classifying isotopic biosignatures, either as a signature consistent with microbial life (biotic) or as a novelty (unclassified/unique). However, machine learning models are likely to be complex and involve interactions between MS features, making biosignatures difficult to interpret. Feature selection methods provide biological and chemical context that help interpret the mechanisms of machine learning models, but these methods also need the ability to detect complex interactions. Previously, we developed a machine learning feature selection algorithm called nearest-neighbor projected distance regression (NPDR) that has the ability to identify important model features that involve complex interactions and automatically reduce correlation and the dimensionality in a high-dimensional variable space. The standard distance metrics used in NPDR – Manhattan and Euclidean – assume the multivariate data are isotropic, which is often violated in real data due to differences in the covariance between variables. Thus, we extend NPDR to include a random forest distance, and other anisotropic distance metrics, for computing nearest neighbors. We also augment the isotope-ratio MS data with time-series features from the raw MS signal to improve biotic classification. We test NPDR on our novel experimental ocean world seawater analog MS data. We measure isotope fractionations of volatile CO 2 that could be measured in exospheres or plumes. Samples include baseline abiotic conditions using a range of possible seawater chemistry consistent with Europa and Enceladus, and biotic samples that include microbes in these seawaters. We use penalized NPDR with random forest proximity to identify interpretable microbial molecular signatures. We compare features with random forest importance, and we train a classifier that discriminates between biotic and abiotic samples with high accuracy. These ML-trained ocean-world analog MS data could be used to assist in identifying biosignatures during future missions.

geochemistry↗

Computation of the factorized error covariance of the difference between correlated estimators

A state estimation problem where some of the measurements may be common to two or more data sets is considered. Two approaches for computing the error covariance of the difference between filtered estimates (for each data set) are discussed. The first algorithm is based on postprocessing of the Kalman gain profiles of two correlated estimators. It uses UD factors of the covariance of the relative error. The second algorithm uses a square root information filter applied to relative error analysis. In the absence of process noise, the square root information filter is computationally more efficient and more flexible than the Kalman gain (covariance update) method. Both the algorithms (covariance and information matrix based) are applied to a Venus orbiter simulation, and their performances are compared.

Wolff, Peter J.↗

Polarimetric remote sensing of sea ice in the Beaufort Sea

Data were collected in March 1988 with the Jet Propulsion Laboratory (JPL) polarimetric airborne synthetic aperture radar (SAR). Full covariance matrices were obtained from the scattering data with proper consideration of the polarimetric calibration. Images were processed to form the sea ice scenes in the Beaufort sea. These scenes contain various ice types including thin ice, thick first-year ice, and multiyear ice. Sea ice was modeled as a layer random medium and the polarimetric backscattering coefficients were calculated under the distorted Born approximation with effective permittivities obtained from the strong fluctuation theory. Sea ice signatures remotely observed by the JPL SAR over the Beaufort sea ice scenes were interpreted by comparing their behavior at different incident angles with results obtained from the theoretical model.

Nghiem, S. V.↗

On estimating the Venus spin vector

The improvement in spin vector and probe position estimates one may reasonably expect from the processing of such data is indicated. This was done by duplicating the ensemble calculations associated with a weighed least squares with a priori estimation technique applied to range rate data that were assumed to be unbiased and uncorrelated. The weighting matrix was assumed to be the inverse of the covariance matrix of the noise on the data. Attention is focused primarily on the spin vector estimation.

Argentiero, P. D.↗

Improving Reanalyses Using TRMM and SSM/I-Derived Precipitation and Total Precipitable Water Observations

Global reanalyses currently contain significant errors in the primary fields of the hydrological cycle such as precipitation, evaporation, moisture, and the related cloud fields, especially in the tropics. The Data Assimilation Office (DAO) at the NASA Goddard Space Flight Center has been exploring the use of rainfall and total precipitable water (TPW) observations from the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) and the Special Sensor Microwave/ Imager (SSM/I) instruments to improve these fields in reanalyses. The DAO has developed a "1+1"D procedure to assimilate 6-hr averaged rainfall and TPW into the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). The algorithm is based on a 6-hr time integration of a column version of the GEOS DAS. The "1+1" designation refers to one spatial dimension plus one temporal dimension. The scheme minimizes the least-square differences between the satellite-retrieved rain rates and those produced by the column model over the 6-hr analysis window. The control variables are analysis increments of moisture within the Incremental Analysis Update (IAU) framework of the GEOS DAS. This 1+1D scheme, in its generalization to four dimensions, is related to the standard 4D variational assimilation but differs in its choice of the control variable. Instead of estimating the initial condition at the beginning of the assimilation cycle, it estimates the constant IAU forcing applied over a 6-hr assimilation cycle. In doing so, it imposes the forecast model as a weak constraint in a manner similar to the variational continuous assimilation techniques. We present results from an experiment in which the observed rain rate and TPW are assumed to be "perfect". They show that assimilating the TMI and SSM/I-derived surface precipitation and TPW observations improves not only the precipitation and moisture fields but also key climate parameters directly linked to convective activities such as clouds, the outgoing longwave radiation, and the large-scale circulation in the tropics. In particular, assimilating these data types reduce the state-dependent systematic errors in the assimilated products. The improved analysis also leads to a better short-range forecast, but the impact is modest compared with improvements in the time-averaged fields. These results suggest that, in the presence of biases and other errors of the forecast model, it is possible to improve the time-averaged "climate content" in the assimilated data without comparable improvements in the short-range forecast skill. Results of this experiment provide a useful benchmark for evaluating error covariance models for optimal use of these data types.

Hou, Arthur Y.↗