Search NASA⌕ Search

SEARCH · Search NASA

Results for “autoregressive model”

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.

126 records · Page 7

Enhancing Short-Range Weather Forecasts through Temporal Variation Encoding: A Multiperiod Embedding Approach

Machine learning (ML) techniques have emerged as promising approaches to improve regional weather forecast accuracy and reliability through data-driven methods. We propose a novel ML-based weather forecasting model, the Multiperiod Embed Net (MPENet). A key distinguishing feature of MPENet is its explicit utilization of the inherent cyclic nature in weather dynamics, unlike the autoregressive strategies commonly used in other ML weather forecasting approaches. Critical cyclic structures are identified via Fourier analyses of dynamic time series. Cyclicity in the convolutional representation is achieved by transforming one-dimensional time series of meteorological variables into two-dimensional tensors based on identified periods. This approach enables the model to leverage intrinsic weather patterns, enhancing regional forecast performance. To demonstrate the effectiveness of MPENet, we conduct a comparative analysis with Nvidia’s FourCastNet. Both models are trained on High-Resolution Rapid Refresh (HRRR) data from 2015 to 2022, over a 192 km × 192 km region in Tennessee. The comparisons are performed locally at two specific locations known to have different weather dynamics due to orographic effects: Crossville, on the relatively flat Cumberland Plateau with fewer topographic airflow disruptions, and Oak Ridge, in the ridge-and-valley region, where airflow is heavily influenced by surrounding valleys and mountains. Our results indicate that FourCastNet achieves strong accuracy at very short lead times, while MPENet maintains competitive skill and shows advantages in capturing temporal evolution over longer periods. Cross-correlation analyses of MPENet and FourCastNet predictions with the HRRR data suggest that encoding critical cyclicity into the network architecture leads to improvements in the forecasting skill.

Artificial intelligence↗

Intercomparison of Deep Learning Model Architectures for Atmospheric River Prediction

With a rapid surge in the application of machine learning (ML) for a diverse range of tasks in climate science, the present study addresses a challenge for climate scientists when selecting the optimal ML or deep learning (DL) architecture for a given application. In particular, a DL intercomparison study was performed with a focus on forecasting the position of atmospheric rivers (ARs) on short-range time scales (up to 5-day lead times). AR predictions from multiple DL architectures, including various types of convolutional autoencoders and a vision transformer (ViT), were compared against ECMWF ERA5 reanalysis and hindcasts from a global climate model. DL models with similar trainable parameters were trained on ERA5 reanalysis data and AR positions derived from a thresholding algorithm to ensure a fair comparison among the DL models. Each model’s performance and accuracy in forecasting AR location and key input fields within a 5-day window were assessed using metrics of root-mean-square error, anomaly correlation, and mean intersection over union. The ViT architecture outperformed other autoencoder models in most of the metrics. Incorporating additional meteorological fields only yielded slight improvements in forecasting certain fields at longer lead times. The results also suggest that a smaller number of input time steps or smaller number of autoregressive steps can achieve better prediction skills, while also improving the overall computational efficiency. This research offers valuable insights into the strengths and weaknesses of different DL techniques for AR forecasting, hopefully guiding the development of improved models for forecasting this phenomenon.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Multivariate Time Series Prediction Techniques for Emulating Noah-LSM Soil Moisture Outputs

Land surface models are crucial tools for many earth science applications including numerical weather prediction, water resource and crop monitoring, and climatological analysis. Given a set of atmospheric forcings, seasonal data, and static parameters, models like Noah-LSM solve for land surface quantities including skin temperature, sensible heat flux, and soil moisture. While these calculations are theoretically robust, they are often computationally expensive. Since artificial neural networks (ANNs) are universal function approximators, they can learn to emulate the output of a deterministic numerical model given a time series of input forcings, with the learned ANN having substantially shorter execution time. The ANN could efficiently parameterize other models, generate ensembles, and provide first-guess inputs for retrievals. As such, with the goal of developing a model that efficiently mimics the output of Noah-LSM given NLDAS2 forcings on a region covering much of the central US, we examine and compare several neural network architectures for the multi-horizon multivariate time series forecasting problem. Recent literature includes a diverse set of approaches including autoregressive architectures like LSTM and GRU, parametric and non-parametric statistical predictors (ForecastNet and MQRNN), self-attention (LSTM-attention-LSTM), and temporal convovlution (DeepTCN). We implement several of these models for the Noah-LSM prediction task, highlighting the features and challenges for each and providing practical insight on the training process.

Mitchell Dodson↗

Generalized Predictive Control for Active Stability Augmentation and Vibration Reduction on an Aeroelastic Tiltrotor Model

Tiltrotor aircraft are defining the state-of-the-art in vertical lift technology as they have the potential to greatly expand rotary-wing operational boundaries. However, they are often limited in forward flight speed due to complex coupled rotor and wing dynamic instabilities. The U.S. Army and NASA have been developing a new wind tunnel model, the TiltRotor Aeroelastic Stability Testbed(TRAST), to test proprotors in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) to investigate aeroelastic stability in cruise. The test is intended to provide high-quality research data for analytical tool development and validation. In addition, the TRAST model will support, develop, and mature new technologies for the design of advanced proprotor aircraft. Stability augmentation and vibration reduction during testing is planned with the use of an active control methodology known as Generalized Predictive Control(GPC). GPC is an autoregressive control law that experimentally acquires a system identification to derive the input-output relation of controls and corresponding sensors. This type of control law is especially useful for complex dynamic interactions that are difficult to explicitly model such as proprotor pylon instability, often referred to as whirl flutter. GPC has been successfully employed on other tiltrotor vehicles to suppress whirl flutter instabilities and vibrations. To aid in the characterization of the wind-tunnel model and in tool development, an analytical representation of the wind-tunnel model was developed using the rotorcraft comprehensive analysis system (RCAS) that simulates structural dynamics and aerodynamics. RCAS was used to derive state-space estimates of the physical plant at various flight conditions to test control law effectiveness. This paper will present an overview of the test article development, a description of RCAS, an explanation of the GPC methodology, and results of GPC being applied to state-space plant estimates of the TRAST model. In these simulations, GPC was effective at stabilizing the aircraft beyond the whirl-flutter boundary while simultaneously reducing vibrations across the flight regime. Additionally, a modern advancement to GPC, termed advanced GPC (AGPC), is introduced that enables a self-adapting system identification. Preliminary results show that AGPC is successful at self-correction as the plant changes from what was used for system identification.

tiltrotor↗

Representation of high frequency Space Shuttle data by ARMA algorithms and random response spectra

High frequency Space Shuttle lift-off data are treated by autoregressive (AR) and autoregressive-moving-average (ARMA) digital algorithms. These algorithms provide useful information on the spectral densities of the data. Further, they yield spectral models which lend themselves to incorporation to the concept of the random response spectrum. This concept yields a reasonably smooth power spectrum for the design of structural and mechanical systems when the available data bank is limited. Due to the non-stationarity of the lift-off event, the pertinent data are split into three slices. Each of the slices is associated with a rather distinguishable phase of the lift-off event, where stationarity can be expected. The presented results are rather preliminary in nature; it is aimed to call attention to the availability of the discussed digital algorithms and to the need to augment the Space Shuttle data bank as more flights are completed.

Spanos, P. D.↗

Spectral representation of high-frequency Space Shuttle data

High frequency Space Shuttle liftoff data are treated by autoregressive (AR) and autoregressive-moving-average (ARMA) digital algorithms. These algorithms provide useful information on the spectral densities of the data. Further, they yield spectral models, which lend themeselves to incorporation into the concept of the random response spectrum. This concept yields a reasonably smooth power spectrum for the design of structural and mechanical systems when the available data bank is limited. Due to the nonstationary of the liftoff event, the pertinent data are split into three slices. Each of the slices is associated with a rather distinguished phase of the liftoff event, in which stationarity can be expected. The presented results are preliminary in nature; they aim to call attention to the availability of the discussed concepts and to the need to augment the Space Shuttle data bank as more flights are completed.

Spanos, P. D.↗

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗

Search for deterministic pulse trends in gamma ray burst temporal profiles

Most cosmic gamma-ray burst temporal profiles appear to be comprised of several individual pulses, many of which overlap. It is advantageous to deconvolve the temporal structures into their constituent pulses, and thereby investigate the shape, intensity and temporal distributions of the pulses as a function of energy. Such fundamental pulse descriptors would provide constraints for theoretical modeling of the burst emission process, such as indications of source size, optical depth and geometry, as a function of time. We have developed a deconvolution algorithm which treats sequences of pulse shapes that change deterministically. The algorithm, a generalization of autoregressive techniques, has been applied to a few bright bursts observed by BATSE. Results indicate that, even within short intervals, constituent pulses are not self-similarly shaped, nor do pulse shapes evolve in a simple manner throughout a burst. Hence, the direction of our future work on pulse deconvolution will focus on analysis methods which allow pulse shape to vary.

Norris, J. P.↗

Construction of Covariance Functions with Variable Length Fields

This article focuses on construction, directly in physical space, of three-dimensional covariance functions parametrized by a tunable length field, and on an application of this theory to reproduce the Quasi-Biennial Oscillation (QBO) in the Goddard Earth Observing System, Version 4 (GEOS-4) data assimilation system. These Covariance models are referred to as multi-level or nonseparable, to associate them with the application where a multi-level covariance with a large troposphere to stratosphere length field gradient is used to reproduce the QBO from sparse radiosonde observations in the tropical lower stratosphere. The multi-level covariance functions extend well-known single level covariance functions depending only on a length scale. Generalizations of the first- and third-order autoregressive covariances in three dimensions are given, providing multi-level covariances with zero and three derivatives at zero separation, respectively. Multi-level piecewise rational covariances with two continuous derivatives at zero separation are also provided. Multi-level powerlaw covariances are constructed with continuous derivatives of all orders. Additional multi-level covariance functions are constructed using the Schur product of single and multi-level covariance functions. A multi-level powerlaw covariance used to reproduce the QBO in GEOS-4 is described along with details of the assimilation experiments. The new covariance model is shown to represent the vertical wind shear associated with the QBO much more effectively than in the baseline GEOS-4 system.

Gaspari, Gregory↗

Implementation of a modal filter on a five meter truss structure

Modal filtering is a spatial filtering technique which uses a weighted sum of a number of response measurements to extract the modal coordinates of the system from the physical response coordinates. No moving average or autoregressive calculations are required to implement the modal filter thus the modal coordinates may be calculated in real time. For practical implementation of the modal filter, the number and location of response locations must be chosen carefully. A modal filter is implemented on a five meter model space truss as a case study. The modal coordinates are extracted in real time using Hewlett Packard 3565 data acquisition and processig hardware. The effect of the number and location of response measurements on the performance of the modal filter is investigated. Applications of the modal filter to modal control and fast parameter identification are also discussed.

Shelly, S.↗

A two dimensional power spectral estimate for some nonstationary processes

A two dimensional estimate for the power spectral density of a nonstationary process is being developed. The estimate will be applied to helicopter noise data which is clearly nonstationary. The acoustic pressure from the isolated main rotor and isolated tail rotor is known to be periodically correlated (PC) and the combined noise from the main and tail rotors is assumed to be correlation autoregressive (CAR). The results of this nonstationary analysis will be compared with the current method of assuming that the data is stationary and analyzing it as such. Another method of analysis is to introduce a random phase shift into the data as shown by Papoulis to produce a time history which can then be accurately modeled as stationary. This method will also be investigated for the helicopter data. A method used to determine the period of a PC process when the period is not know is discussed. The period of a PC process must be known in order to produce an accurate spectral representation for the process. The spectral estimate is developed. The bias and variability of the estimate are also discussed. Finally, the current method for analyzing nonstationary data is compared to that of using a two dimensional spectral representation. In addition, the method of phase shifting the data is examined.

Smith, Gregory L.↗

Autoregressive harmonic analysis of the earth's polar motion using homogeneous International Latitude Service data

The homogeneous set of 80-year-long (1900-1979) International Latitude Service (ILS) polar motion data is analyzed using the autoregressive method (Chao and Gilbert, 1980), which resolves and produces estimates for the complex frequency (or frequency and Q) and complex amplitude (or amplitude and phase) of each harmonic component in the data. The ILS data support the multiple-component hypothesis of the Chandler wobble. It is found that the Chandler wobble can be adequately modeled as a linear combination of four (coherent) harmonic components, each of which represents a steady, nearly circular, prograde motion. The four-component Chandler wobble model 'explains' the apparent phase reversal during 1920-1940 and the pre-1950 empirical period-amplitude relation. The annual wobble is shown to be rather stationary over the years both in amplitude and in phase, and no evidence is found to support the large variations reported by earlier investigations. The Markowitz wobble is found to be marginally retrograde and appears to have a complicated behavior which cannot be resolved because of the shortness of the data set.

Chao, B. F.↗

Predictability of Malaria Transmission Intensity in the Mpumalanga Province, South Africa, Using Land Surface Climatology and Autoregressive Analysis

There has been increasing effort in recent years to employ satellite remotely sensed data to identify and map vector habitat and malaria transmission risk in data sparse environments. In the current investigation, available satellite and other land surface climatology data products are employed in short-term forecasting of infection rates in the Mpumalanga Province of South Africa, using a multivariate autoregressive approach. The climatology variables include precipitation, air temperature and other land surface states computed by the Off-line Land-Surface Global Assimilation System (OLGA) including soil moisture and surface evaporation. Satellite data products include the Normalized Difference Vegetation Index (NDVI) and other forcing data used in the Goddard Earth Observing System (GEOS-1) model. Predictions are compared to long- term monthly records of clinical and microscopic diagnoses. The approach addresses the high degree of short-term autocorrelation in the disease and weather time series. The resulting model is able to predict 11 of the 13 months that were classified as high risk during the validation period, indicating the utility of applying antecedent climatic variables to the prediction of malaria incidence for the Mpumalanga Province.

Grass, David↗

Sea spikes at moderate incidence and their relation to position on the waves

Most models of radar backscatter from the sea ignore sea spikes, the nonlinear effects that result in large excursions above the local mean signal. They are strong enough to change the mean scattering level significantly, and they may cause streaks in synthetic aperture radar (SAR) ocean images. A threshold based on short signal excursions above the local mean is used to identify spikes. With this method, spikes can be found in regions of low signal level. An autoregressive spectral method is used to identify the locations of the spikes on the dominant waves. Sample results from Ka-band measurements in the North Sea made during the SAXON-FPN experiment are given. It is shown that spikes can occur anywhere on the dominant wave, although they are most prevalent on the front face.

Salam, A.↗

Short-term fluctuations in the eddy heat flux and baroclinic stability of the atmosphere

National Meteorological Center data from midlatitudes for three Januaries is used in calculating time series of the zonal mean meridional eddy heat flux and the zonal mean baroclinic stability, as measured by the difference between the zonal wind shear and the critical value of the shear in two-level models. Time-lagged correlations between the two series reveal a highly significant negative correlation for short time lags, peaking at approximately -0.4 when the stability parameter lags one half day behind the eddy flux. They also reveal that strongly unstable conditions are not followed by significant increases in the eddy flux. These results are seen as indicating that the synoptic variations of the zonal mean eddy flux are not closely related to the degree of baroclinic instability of the zonal mean flow. The autocorrelation of the eddy flux is then compared with those expected for autoregressive processes. A Bayesian information criterion suggests that the behavior is represented best by a damped oscillation, with a damping time of 0.8 day and a period of five days.

Stone, P. H.↗

Simulation of ocean SAR images via phase history generation

A method for simulating an ocean synthetic-aperture radar (SAR) image is illustrated for a simple internal wave current pattern. The method calculates both the amplitude image and the radar signal history. The simulation model used consists of six stages. In the first stage, the full wave spectrum in two spatial coordinates is calculated from the wind speed and direction using the action spectral density equation. In the second stage, the pixel size is selected and used as the basis for dividing the spectrum into large and small scale motions. Realizations for the large-scale ocean surface height and velocity are then calculated. In the third stage, the sensor wavelength and geometry are used to calculate the small scale statistics (radar cross section, coherence time, root-mean square (RMS), radial velocity and RMS slope). In the fourth stage, an autoregressive method is used to generate a realization of the surface reflectivity history that is consistent with the radial velocity and the radial velocity variance. In the fifth stage, the signal history is generated by summing the reflectivities at the proper times with the appropriate antenna weighting function for the SAR. As a consequence of this process, speckle is automatically included in the signal. In the final stage, the SAR image can be created by using any of the traditional ways to process the signal history, including variable focusing and multilook processing.

Bennett, John R.↗

Respiratory sinus arrhythmia: time domain characterization using autoregressive moving average analysis

Fourier-based techniques are mathematically noncausal and are therefore limited in their application to feedback-containing systems, such as the cardiovascular system. In this study, a mathematically causal time domain technique, autoregressive moving average (ARMA) analysis, was used to parameterize the relations of respiration and arterial blood pressure to heart rate in eight humans before and during total cardiac autonomic blockade. Impulse-response curves thus generated showed the relation of respiration to heart rate to be characterized by an immediate increase in heart rate of 9.1 +/- 1.8 beats.min-1.l-1, followed by a transient mild decrease in heart rate to -1.2 +/- 0.5 beats.min-1.l-1 below baseline. The relation of blood pressure to heart rate was characterized by a slower decrease in heart rate of -0.5 +/- 0.1 beats.min-1.mmHg-1, followed by a gradual return to baseline. Both of these relations nearly disappeared after autonomic blockade, indicating autonomic mediation. Maximum values obtained from the respiration to heart rate impulse responses were also well correlated with frequency domain measures of high-frequency "vagal" heart rate control (r = 0.88). ARMA analysis may be useful as a time domain representation of autonomic heart rate control for cardiovascular modeling.

NASA Discipline Regulatory Physiology↗