Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multivariate time series”

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 37 records · Page 2

Identification of multivariable high performance turbofan engine dynamics from closed loop data

The multivariable instrumental variable/approximate maximum likelihood (IV/AML) method of recursive time-series analysis is used to identify the multivariable (four inputs-three outputs) dynamics of the Pratt and Whitney F100 engine. A detailed nonlinear engine simulation is used to determine linear engine model structures and parameters at an operating point using open loop data. Also, the IV/AML method is used in a direct identification made to identify models from actual closed loop engine test data. Models identified from simulated and test data are compared to determine a final model structure and parameterization that can predict engine response for a wide class of inputs. The ability of the IV/AML algorithm to identify useful dynamic models from engine test data is assessed. Previously announced in STAR as N82-20339

Merrill, W.↗

Antarctic Ice-Sheet Mass Balance from Satellite Altimetry 1992 to 2001

A major uncertainty in understanding the causes of the current rate of sea level rise is the potential contributions from mass imbalances of the Greenland and Antarctic ice sheets. Estimates of the current mass balance of the Antarctic ice sheet are derived from surface- elevation changes obtained from 9 years of ERS - 1 & 2 radar altimeter data. Elevation time-series are created from altimeter crossovers among 90-day data periods on a 50 km grid to 81.5 S. The time series are fit with a multivariate linear/sinusoidal function to give the average rate of elevation change (dH/dt). On the major Rome-Filchner, Ross, and Amery ice shelves, the W d t are small or near zero. In contrast, the ice shelves of the Antarctic Peninsula and along the West Antarctic coast appear to be thinning significantly, with a 23 +/- 3 cm per year surface elevation decrease on the Larsen ice shelf and a 65 +/- 4 cm per year decrease on the Dotson ice shelf. On the grounded ice, significant elevation decreases are obtained over most of the drainage basins of the Pine Island and Thwaites glaciers in West Antarctica and inland of Law Dome in East Antarctica. Significant elevation increases are observed within about 200 km of the coast around much of the rest of the ice sheet. Farther inland, the changes are a mixed pattern of increases and decreases with increases of a few centimeters per year at the highest elevations of the East Antarctic plateau. The derived elevation changes are combined with estimates of the bedrock uplift from several models to provide maps of ice thickness change. The ice thickness changes enable estimates of the ice mass balances for the major drainage basins, the overall mass balance, and the current contribution of the ice sheet to global sea level change.

Zwally, H. Jay↗

Statistical Evaluation of Time Series Analysis Techniques

The performance of a modified version of NASA's multivariate spectrum analysis program is discussed. A multiple regression model was used to make the revisions. Performance improvements were documented and compared to the standard fast Fourier transform by Monte Carlo techniques.

Benignus, V. A.↗

Seasonal and Interannual Variations of Ice Sheet Surface Elevation at the Summit of Greenland: Observed and Modeled

Observed seasonal and interannual variations in the surface elevation over the summit of the Greenland ice sheet are modeled using a new temperature-dependent formulation of firn-densification and observed accumulation variations. The observed elevation variations are derived from ERS (European Remote Sensing)-1 and ERS-2 radar altimeter data for the period between April 1992 and April 1999. A multivariate linear/sine function is fitted to an elevation time series constructed from elevation differences measured by radar altimetry at orbital crossovers. The amplitude of the seasonal elevation cycle is 0.25 m peak-to-peak, with a maximum in winter and a minimum in summer. Inter-annually, the elevation decreases to a minimum in 1995, followed by an increase to 1999, with an overall average increase of 4.2 cm a(exp -1) for 1992 to 1999. Our densification formulation uses an initial field-density profile, the AWS (automatic weather station) surface temperature record, and a temperature-dependent constitutive relation for the densification that is based on laboratory measurements of crystal growth rates. The rate constant and the activation energy commonly used in the Arrhenius-type constitutive relation for firn densification are also temperature dependent, giving a stronger temperature and seasonal amplitudes about 10 times greater than previous densification formulations. Summer temperatures are most important, because of the strong non-linear dependence on temperature. Much of firn densification and consequent surface lowering occurs within about three months of the summer season, followed by a surface build-up from snow accumulation until spring. Modeled interannual changes of the surface elevation, using the AWS measurements of surface temperature and accumulation and results of atmospheric modeling of precipitation variations, are in good agreement with the altimeter observations. In the model, the surface elevation decreases about 20 cm over the seven years due to more compaction driven by increasing summer temperatures. The minimum elevation in 1995 is driven mainly by a temporary accumulation decrease and secondarily by warmer temperatures. However, the overall elevation increase over the seven years is dominated by the accumulation increase in the later years.

Zwally, H. Jay↗

MARGInS Model-Based Analysis of Realizable Goals in Systems

The high complexity of modern aircraft and spacecraft requires elaborate Verification and Validation (V&V) approaches to make sure that such complex systems work properly and reliably. MARGInS is a framework for the analysis, understanding, and prediction of the behavior of a complex, hybrid system. MARGInS contains a set of machine learning and statistical algorithms for multivariate clustering, treatment learning, critical factor determination, time-series analysis, event prediction, and safety-boundary detection and characterization. The framework supports system testing and can be configured to find novel features in test suites, determine classes of behavior, propose new experiments that can efficiently explore and characterize the boundaries between classes of system behavior, and to create visualizations and reports.

He, Yuning↗

Enhancements of Bayesian Blocks; Application to Large Light Curve Databases

Bayesian Blocks are optimal piecewise linear representations (step function fits) of light-curves. The simple algorithm implementing this idea, using dynamic programming, has been extended to include more data modes and fitness metrics, multivariate analysis, and data on the circle (Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations, Scargle, Norris, Jackson and Chiang 2013, ApJ, 764, 167), as well as new results on background subtraction and refinement of the procedure for precise timing of transient events in sparse data. Example demonstrations will include exploratory analysis of the Kepler light curve archive in a search for "star-tickling" signals from extraterrestrial civilizations. (The Cepheid Galactic Internet, Learned, Kudritzki, Pakvasa1, and Zee, 2008, arXiv: 0809.0339; Walkowicz et al., in progress).

Kepler light curve archive↗

Application of multivariate autoregressive spectrum estimation to ULF waves

The estimation of the power spectrum of a time series by fitting a finite autoregressive model to the data has recently found widespread application in the physical sciences. The extension of this method to the analysis of vector time series is presented here through its application to ULF waves observed in the magnetosphere by the ATS 6 synchronous satellite. Autoregressive spectral estimates of the power and cross-power spectra of these waves are computed with computer programs developed by the author and are compared with the corresponding Blackman-Tukey spectral estimates. The resulting spectral density matrices are then analyzed to determine the direction of propagation and polarization of the observed waves.

Ioannidis, G. A.↗

Stratospheric Ozone Trends and Variability as Seen by SCIAMACHY from 2002 to 2012

Vertical profiles of the rate of linear change (trend) in the altitude range 15-50 km are determined from decadal O3 time series obtained from SCIAMACHY/ENVISAT measurements in limb-viewing geometry. The trends are calculated by using a multivariate linear regression. Seasonal variations, the quasi-biennial oscillation, signatures of the solar cycle and the El Nino-Southern Oscillation are accounted for in the regression. The time range of trend calculation is August 2002-April 2012. A focus for analysis are the zonal bands of 20 deg N - 20 deg S (tropics), 60 - 50 deg N, and 50 - 60 deg S (midlatitudes). In the tropics, positive trends of up to 5% per decade between 20 and 30 km and negative trends of up to 10% per decade between 30 and 38 km are identified. Positive O3 trends of around 5% per decade are found in the upper stratosphere in the tropics and at midlatitudes. Comparisons between SCIAMACHY and EOS MLS show reasonable agreement both in the tropics and at midlatitudes for most altitudes. In the tropics, measurements from OSIRIS/Odin and SHADOZ are also analysed. These yield rates of linear change of O3 similar to those from SCIAMACHY. However, the trends from SCIAMACHY near 34 km in the tropics are larger than MLS and OSIRIS by a factor of around two.

stratosphere↗

Analysis models for the estimation of oceanic fields

A general model for statistically optimal estimates is presented for dealing with scalar, vector and multivariate datasets. The method deals with anisotropic fields and treats space and time dependence equivalently. Problems addressed include the analysis, or the production of synoptic time series of regularly gridded fields from irregular and gappy datasets, and the estimate of fields by compositing observations from several different instruments and sampling schemes. Technical issues are discussed, including the convergence of statistical estimates, the choice of representation of the correlations, the influential domain of an observation, and the efficiency of numerical computations.

Carter, E. F.↗

Surface Elevation Changes in West Antarctica from Satellite Altimetry: Mass Balance Implications

Time-series of surface elevation change, which are constructed from 7-years (1992-1999) of ERS-1 and 2 satellite radar altimeter data of Antarctica, show significant seasonal, inter-annual, and long-term changes. Elevation time-series are created from altimeter crossovers among 90-day data periods on a 50 km grid to 81.5 degrees S and fit with a multivariate linear/sinusoidal function to give the average rate of elevation change (dH/dt) and account for seasonal changes. On the major Ronne, Filchner, and Ronne ice shelves, the dH/dt are small or near zero. In contrast, the ice shelves of the Antarctic Peninsula and along the West Antarctic coast appear to be thinning significantly, with a 23 +/- 3 cm a(exp -1) surface elevation decrease on the Larsen ice shelf and a 65 +/- 4 cm a(exp -1) decrease on the Dotson ice shelf. Significant elevation decreases are obtained over most of the drainage basins of the Pine Island and Thwaites glaciers. Significant increases are obtained over most of the other grounded ice in Marie Byrd Land, the Antarctic Peninsula, and Coates Land. Over the sector from 85 degrees W to 115 degrees W, which includes the Pine Island and Thwaites basins, the average elevation is significantly decreasing by 8.1 cm a(exp -1). The corresponding ice thickness change is about -11 cm a(exp -1), with a corresponding mass loss of 82 Gt a(exp -1), and a 0.22 mm a(exp -1) contribution to global sea level rise. In terms of elevation change, the decrease in the Pine Island-Thwaites sector is largely balanced by the increase in the Marie Byrd Land, but only balanced by about 1/4 in terms of ice thickness change and contribution to sea level rise. The overall average elevation change for the grounded ice is + 1.2 cm a(exp -1). Using an average bedrock uplift of 2.5 cm a(exp -1), implies an average ice thickness decrease of 1.3 cm a(exp -1), a mass loss of 22 Gt a(exp -1), and a 0.06 mm a(exp -1) contribution to global sea level rise.

Zwally, H. Jay↗

On the Choice of Variable for Atmospheric Moisture Analysis

The implications of using different control variables for the analysis of moisture observations in a global atmospheric data assimilation system are investigated. A moisture analysis based on either mixing ratio or specific humidity is prone to large extrapolation errors, due to the high variability in space and time of these parameters and to the difficulties in modeling their error covariances. Using the logarithm of specific humidity does not alleviate these problems, and has the further disadvantage that very dry background estimates cannot be effectively corrected by observations. Relative humidity is a better choice from a statistical point of view, because this field is spatially and temporally more coherent and error statistics are therefore easier to obtain. If, however, the analysis is designed to preserve relative humidity in the absence of moisture observations, then the analyzed specific humidity field depends entirely on analyzed temperature changes. If the model has a cool bias in the stratosphere this will lead to an unstable accumulation of excess moisture there. A pseudo-relative humidity can be defined by scaling the mixing ratio by the background saturation mixing ratio. A univariate pseudo-relative humidity analysis will preserve the specific humidity field in the absence of moisture observations. A pseudorelative humidity analysis is shown to be equivalent to a mixing ratio analysis with flow-dependent covariances. In the presence of multivariate (temperature-moisture) observations it produces analyzed relative humidity values that are nearly identical to those produced by a relative humidity analysis. Based on a time series analysis of radiosonde observed-minus-background differences it appears to be more justifiable to neglect specific humidity-temperature correlations (in a univariate pseudo-relative humidity analysis) than to neglect relative humidity-temperature correlations (in a univariate relative humidity analysis). A pseudo-relative humidity analysis is easily implemented in an existing moisture analysis system, by simply scaling observed-minus background moisture residuals prior to solving the analysis equation, and rescaling the analyzed increments afterward.

Dee, Dick P.↗

Regional climate change predictions from the Goddard Institute for Space Studies high resolution GCM

A new diagnostic tool is developed for examining relationships between the synoptic scale circulation and regional temperature distributions in GCMs. The 4 x 5 deg GISS GCM is shown to produce accurate simulations of the variance in the synoptic scale sea level pressure distribution over the U.S. An analysis of the observational data set from the National Meteorological Center (NMC) also shows a strong relationship between the synoptic circulation and grid point temperatures. This relationship is demonstrated by deriving transfer functions between a time-series of circulation parameters and temperatures at individual grid points. The circulation parameters are derived using rotated principal components analysis, and the temperature transfer functions are based on multivariate polynomial regression models. The application of these transfer functions to the GCM circulation indicates that there is considerable spatial bias present in the GCM temperature distributions. The transfer functions are also used to indicate the possible changes in U.S. regional temperatures that could result from differences in synoptic scale circulation between a 1XCO2 and a 2xCO2 climate, using a doubled CO2 version of the same GISS GCM.

Crane, Robert G.↗

Using Concurrent Cardiovascular Information to Augment Survival Time Data from Orthostatic Tilt Tests

Orthostatic Intolerance (OI) is the propensity to develop symptoms of fainting during upright standing. OI is associated with changes in heart rate, blood pressure and other measures of cardiac function. Problem: NASA astronauts have shown increased susceptibility to OI on return from space missions. Current methods for counteracting OI in astronauts include fluid loading and the use of compression garments. Multivariate trajectory spread is greater as OI increases. Pairwise comparisons at the same time within subjects allows incorporation of pass/fail outcomes. Path length, convex hull area, and covariance matrix determinant do well as statistics to summarize this spread Missing data problems Time series analysis need many more time points per OTT session treatment of trend? how incorporate survival information?

Feiveson, Alan H.↗

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↗

X-Windows PVT Widget Class

The X-Windows Process Validation Table (PVT) Widget Class ( Class is used here in the object-oriented-programming sense of the word) was devised to simplify the task of implementing network registration services for Information Sharing Protocol (ISP) graphical-user-interface (GUI) computer programs. Heretofore, ISP PVT programming tasks have required many method calls to identify, query, and interpret the connections and messages exchanged between a client and a PVT server. Normally, programmers have utilized direct access to UNIX socket libraries to implement the PVT protocol queries, necessitating the use of many lines of source code to perform frequent tasks. Now, the X-Windows PVT Widget Class encapsulates ISP client server network registration management tasks within the framework of an X Windows widget. Use of the widget framework enables an X Windows GUI program to interact with PVT services in an abstract way and in the same manner as that of other graphical widgets, making it easier to program PVT clients. Wrapping the PVT services inside the widget framework enables a programmer to treat a PVT server interface as though it were a GUI. Moreover, an alternate subclass could implement another service in a widget of the same type. This program was written by Matthew R. Barry of United Space Alliance for Johnson Space Center. For further information, contact the Johnson Technology Transfer Office at (281) 483-3809. MSC-23582 Shuttle Data Center File- Processing Tool in Java A Java-language computer program has been written to facilitate mining of data in files in the Shuttle Data Center (SDC) archives. This program can be executed on a variety of workstations or via Web-browser programs. This program is partly similar to prior C-language programs used for the same purpose, while differing from those programs in that it exploits the platform-neutrality of Java in implementing several features that are important for analysis of large sets of time-series data. The program supports regular expression queries of SDC archive files, reads the files, interleaves the time-stamped samples according to a chosen output, then transforms the results into that format. A user can choose among a variety of output file formats that are useful for diverse purposes, including plotting, Markov modeling, multivariate density estimation, and wavelet multiresolution analysis, as well as for playback of data in support of simulation and testing.

Barry, Matthew R.↗

Shuttle Data Center File-Processing Tool in Java

A Java-language computer program has been written to facilitate mining of data in files in the Shuttle Data Center (SDC) archives. This program can be executed on a variety of workstations or via Web-browser programs. This program is partly similar to prior C-language programs used for the same purpose, while differing from those programs in that it exploits the platform-neutrality of Java in implementing several features that are important for analysis of large sets of time-series data. The program supports regular expression queries of SDC archive files, reads the files, interleaves the time-stamped samples according to a chosen output, then transforms the results into that format. A user can choose among a variety of output file formats that are useful for diverse purposes, including plotting, Markov modeling, multivariate density estimation, and wavelet multiresolution analysis, as well as for playback of data in support of simulation and testing.

Barry, Matthew R.↗

Modeling and control for closed environment plant production systems

A computer program was developed to study multiple crop production and control in controlled environment plant production systems. The program simulates crop growth and development under nominal and off-nominal environments. Time-series crop models for wheat (Triticum aestivum), soybean (Glycine max), and white potato (Solanum tuberosum) are integrated with a model-based predictive controller. The controller evaluates and compensates for effects of environmental disturbances on crop production scheduling. The crop models consist of a set of nonlinear polynomial equations, six for each crop, developed using multivariate polynomial regression (MPR). Simulated data from DSSAT crop models, previously modified for crop production in controlled environments with hydroponics under elevated atmospheric carbon dioxide concentration, were used for the MPR fitting. The model-based predictive controller adjusts light intensity, air temperature, and carbon dioxide concentration set points in response to environmental perturbations. Control signals are determined from minimization of a cost function, which is based on the weighted control effort and squared-error between the system response and desired reference signal.

NASA Discipline Life Support Systems↗