Search NASA⌕ Search

SEARCH · Search NASA

Results for “Functional Principal Components”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

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.↗
DOE OSTI2019

A Novel Statistical Method for Measuring the Temperature–Density Relation in the IGM Using the b – N Hi Distribution of Absorbers in the Ly α Forest

We propose a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Lyα forest. The distribution of Doppler parameter and column density (b–N H i distribution) is sensitive to the temperature–density relation T = T 0 (ρ/ρ 0 ) γ-1 , and previous work has inferred T 0 and γ by fitting its low-b cutoff. This approach discards the majority of available data and is susceptible to systematics related to cutoff determination. We introduce an approach that exploits all information encoded in the b – N H i distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, and then we use principal component decomposition to create an emulator that can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov Chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements, we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D ( T 0 , γ) measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at z = 2. Using 200 absorbers, equivalent in path length to a single Lya forest spectrum, we measure $\mathrm{log}{T}_{0}={4.092}_{-0.055}^{+0.050}$ and $\gamma ={1.49}_{-0.074}^{+0.073}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of log T 0 and γ with precision on $\mathrm{log}{T}_{0}$ (γ) nearly two (three) times higher for current data set sizes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

A unified development of several techniques for the representation of random vectors and data sets

Linear vector space theory is used to develop a general representation of a set of data vectors or random vectors by linear combinations of orthonormal vectors such that the mean squared error of the representation is minimized. The orthonormal vectors are shown to be the eigenvectors of an operator. The general representation is applied to several specific problems involving the use of the Karhunen-Loeve expansion, principal component analysis, and empirical orthogonal functions; and the common properties of these representations are developed.

Bundick, W. T.↗
NASA NTRS2011

TPSAS-NF1676L-12839-DND

Hyper-spectral remote sensors, such as Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI), provides top of atmospheric radiances with high information content on atmospheric and surface properties. In order to analyze these data in real time, fast, and accurate forward and inverse models are needed. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The method has been successfully applied to AIRS, IASI, and NAST-I data. By applying the same method to the current and future hyperspectral sounders, the derived products will have less errors due to biases introduced by different retrieval methods. We have compared the retrieved products to radiosondes, aircraft measurements and other validation dataset. Simulations have done to quantify retrieval errors associated with the retrievals. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a subset of the channels. The PCRTM forward model has also been used to performance end-to-end sensor performance simulations of the Climate Absolute Radiance and Refractivity Observatory (CLARREO).

Xu Liu↗
DOE OSTI2019

Determining moisture content of Traditional Chinese Medicines using a near-infrared LED-based moisture content sensor with spectrum analysis

In Traditional Chinese Medicine (TCM), moisture content (MC) inspection using voltage level or a lock-in amplifier, with a features-based near-infrared light-emitting diode sensor is possible. Abundant absorption response in a wide frequency range and information in both the time-domain and frequency-domain can be obtained, with the square excitation method. However, planeness, heterogeneous particle sizes, and measurement noise negatively affect the time-domain feature during MC inspection. In this work, spectrum analysis (SA) for feature extraction is investigated, to precisely measure the MC of Dutong pills. The principal components for low and high frequencies, are aimed at counteracting the dimensional methods, which are used for supplying the lower dimensional features. The SA experimental results demonstrate that the low frequency components, ranging from the base to the 3rd harmonic components, have a high linear correlation with the MC. The calibration and validation experimental results reveal that the MC can be a monotonic function about the SA amplitudes or first principal component (PC1), with the proposed sensor. This method has the potential to be used for gauging automatic, in situ MCs, particularly for quality control and Process Analytical Technology in TCM. Future research opportunities include improving accuracy and application of the method under field conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

PCAfold 2.0—Novel tools and algorithms for low-dimensional manifold assessment and optimization

We describe an update to our open-source Python package, PCAfold, designed to help researchers generate, analyze and improve low-dimensional data manifolds. In the current version, PCAfold 2.0, we introduce novel tools and algorithms for assessing and optimizing low-dimensional manifolds. This includes a method that generates a “map” of local feature sizes that can help pinpoint researchers to problematic regions on a manifold. We introduce a novel cost function that characterizes the quality of a manifold topology with a single number. We develop two algorithms for feature selection based on principal component analysis (PCA) that use the cost function as an objective function to minimize. We introduce a quantity of interest (QoI)-aware dimensionality reduction strategy where data projections are computed using an artificial neural network and are directly optimized towards representing various projection-independent and projection-dependent QoIs. We also introduce an implementation of partition of unity networks (POUnets) for efficient reconstruction of QoIs from low-dimensional manifolds based on combining neural network classification with localized polynomial regression. Our software can be broadly applicable in all domains of science and engineering that aim to reduce data dimensionality, as well as in the fundamental research on representation learning.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗
NASA NTRS1980

Executive and communications services to support the IPAD environment

The principal purposes of the prototype executive software are to provide a system independent interface to the underlying host system and to allow for extension to full IPAD executive services as described in the preliminary design. A basic set of functions is included in the prototype to meet the requirements of the other components of the prototype, principally IPID, the IPAD data management system. The functions were chosen so that they would be readily built on any of the proposed host systems with minimal redesign and execution overhead. The functions fall into five categories: access to host data, access to data files, access to communication services, data transformation, and instrumentation for performance measurement. Communication services provide message delivery between processes in a network of heterogeneous computers. Data transformation services and communication services ensure data type validity and data integrity of messages exchanged between processes.

Tanner, J. G.↗
DOE OSTI2018

Autoencoder node saliency: Selecting relevant latent representations

The autoencoder is an artificial neural network that performs nonlinear dimension reduction and learns hidden representations of unlabeled data. With a linear transfer function it is similar to the principal component analysis (PCA). While both methods use weight vectors for linear transformations, the autoencoder does not come with any indication similar to the eigenvalues in PCA that are paired with eigenvectors. In this work, we propose a novel autoencoder node saliency method that examines whether the features constructed by autoencoders exhibit properties related to known class labels. The supervised node saliency ranks the nodes based on their capability of performing a learning task. It is coupled with the normalized entropy difference (NED). We establish a property for NED values to verify classifying behaviors among the top ranked nodes. Lastly, by applying our methods to real datasets, we demonstrate their ability to provide indications on the performing nodes and explain the learned tasks in autoencoders.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Intraseasonal oscillations in the global atmosphere. II - Southern Hemisphere

A systematic examination of oscillatory modes in the global atmosphere is completed by studying 12 years of 500 mb geopotential heights in the Southern Hemisphere. The data were band-pass filtered to focus on intraseasonal (IS) phenomena, and spatial empirical orthogonal functions (EOFs) were obtained. The leading principal components were subjected to singular spectrum analysis in order to identify nonlinear IS oscillations with high statistical confidence. In the Southern Hemisphere, the dominant mode has a period of 23 days, with spatial patterns carried by the second and third winter EOF of the IS band. It has a zonal wavenumber-four structure. The 40-day mode is second, and is dominated by wavenumbers three and four, while a 16-day mode is too weak to separate its spatial behavior from the previous two. The IS dynamics in the Southern Hemisphere is more complex and dominated by shorter wavenumbers than the Northern Hemisphere. No statistically significant correlations between the Southern Hemisphere and the tropics or the Northern Hemisphere are apparent in the IS band.

Ghil, Michael↗
NASA NTRS2000

Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.

Liu, Antony K.↗
NASA NTRS2005

The North Alabama Lightning Mapping Array (LMA): A Network Overview

The North Alabama Lightning Mapping Array (LMA) is s a 3-D VHF regional lightning detection system that provides on-orbit algorithm validation and instrument performance assessments for the NASA Lightning Imaging Sensor, as well as information on storm kinematics and updraft evolution that offers the potential to improve severe storm warning lead time by up t o 50% and decrease te false alarm r a t e ( for non-tornado producing storms). In support of this latter function, the LMA serves as a principal component of a severe weather test bed to infuse new science and technology into the short-term forecasting of severe and hazardous weather, principally within nearby National Weather Service forecast offices. The LMA, which became operational i n November 2001, consists of VHF receivers deployed across northern Alabama and a base station located at the National Space Science and Technology Center (NSSTC), which is on t h e campus of the University of Alabama in Huntsville. The LMA system locates the sources of impulsive VHF radio signals s from lightning by accurately measuring the time that the signals aririve at the different receiving stations. Each station's records the magnitude and time of the peak lightning radiation signal in successive 80 ms intervals within a local unused television channel (channel 5, 76-82 MHz in our case ) . Typically hundreds of sources per flash can be reconstructed, which i n t u r n produces accurate 3-dimensional lightning image maps (nominally <50 m error within 150 la. range). The data are transmitted back t o a base station using 2.4 GHz wireless Ethernet data links and directional parabolic grid antennas. There are four repeaters in the network topology and the links have an effective data throughput rate ranging from 600 kbits s -1 t o 1.5 %its s -1. This presentation provides an overview of t h e North Alabama network, the data processing (both real-time and post processing) and network statistics.

Blakeslee, R. J.↗
NASA NTRS2013

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data: 1. an EOF Approach to the Spatial-Temporal Variability of Aerosol Optical Depth Using Multiple Remote Sensing Data Sets

Many remote sensing techniques and passive sensors have been developed to measure global aerosol properties. While instantaneous comparisons between pixel-level data often reveal quantitative differences, here we use Empirical Orthogonal Function (EOF) analysis, also known as Principal Component Analysis, to demonstrate that satellite-derived aerosol optical depth (AOD) data sets exhibit essentially the same spatial and temporal variability and are thus suitable for large-scale studies. Analysis results show that the first four EOF modes of AOD account for the bulk of the variance and agree well across the four data sets used in this study (i.e., Aqua MODIS, Terra MODIS, MISR, and SeaWiFS). Only SeaWiFS data over land have slightly different EOF patterns. Globally, the first two EOF modes show annual cycles and are mainly related to Sahara dust in the northern hemisphere and biomass burning in the southern hemisphere, respectively. After removing the mean seasonal cycle from the data, major aerosol sources, including biomass burning in South America and dust in West Africa, are revealed in the dominant modes due to the different interannual variability of aerosol emissions. The enhancement of biomass burning associated with El Niño over Indonesia and central South America is also captured with the EOF technique.

variability↗
NASA NTRS2026

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)↗
DOE OSTI2025

Birefringence axis rotation and nonlinear response of a side-hole fiber pressure sensor

The pressure-induced birefringence of side-hole fibers has been widely explored for hydrostatic pressure measurement. Here, this paper studies the effects of the initial birefringence of a side-hole fiber on the pressure response of polarimetric sensors based on side-hole fibers. A mathematical model, where the stress field induced by the pressure is assumed to be uniform in the core region of the fiber, was established to calculate the birefringence and optical axis rotation as a function of the differences between the two principal stress components of the stress field in the core region. Simulation results show that when there is misalignment between the initial birefringence axes and the principal axes of the stress field, pressure changes cause rotations of the birefringence axis, leading to complicated and nonlinear changes of birefringence as a function of pressure. Experiments were conducted using fiber-Bragg-grating-based resonators on side-hole fibers with different initial birefringence, and the results agree well with the simulation. The results underscore the importance of the control of magnitude and direction of the initial birefringence of the side-hole fiber and provide guidelines for the design and fabrication of side-hole fiber pressure sensors.

42 ENGINEERING↗
NASA NTRS1976

Three-axis attitude control for a geostationary satellite

The Applications Technology Satellite-6 was successfully launched into a geostationary orbit in May, 1974. Since that time, its attitude control subsystem has routinely demonstrated its ability to perform precision pointing and precision closed loop maneuvers. This paper describes the ATS-6 attitude control subsystem. The requirements imposed by the spacecraft mission objectives are outlined, and functional descriptions of the principal attitude control subsystem components and control laws are presented. The paper concludes with a summary comparison of in-orbit performance relative to its specified performance requirements.

Freeman, H. R.↗
NASA NTRS1984

Applications of array processors in the analysis of remote sensing images

The architectures, programming characteristics, and ranges of application of past, present, and planned array processors for the digital processing of remote-sensing images are compared. Such functions as radiometric and geometric corrections, principal-components analysis, cluster coding, histogram generation, grey-level mapping, convolution, classification, and mensuration and modeling operations are considered, and both pipeline-type and single-instruction/multiple-data-stream (SIMD) arrays are evaluated. Numerical results are presented in a table, and it is found that the pipeline-type arrays normally used with minicomputers increase their speed significantly at low cost, while even further gains are provided by the more expensive SIMD arrays. Most image-processing operations become I/O-limited when SIMD arrays are used with current I/O devices.

Ramapriyan, H. K.↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗