Search NASA⌕ Search

SEARCH · Search NASA

Results for “Empirical orthogonal functions”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Long-Term, Non-Anthropogenic Groundwater Storage Changes Simulated by Three Global-Scale Hydrological Models

This study examined long-term, natural (i.e., excluding anthropogenic impacts) variability of groundwater storage worldwide. Groundwater storage changes were estimated by forcing three global-scale hydrological models with three 50+ year meteorological datasets. Evaluation using in situ groundwater observations from the U.S. and terrestrial water storage derived from the Gravity Recovery and Climate Experiment (GRACE) satellites showed that these models reasonably represented inter-annual variability of water storage, as indicated by correlations greater than 0.5 in most regions. Empirical orthogonal function analysis revealed influences of the El Niño Southern Oscillation (ENSO) on global groundwater storage. Simulated groundwater storage, including its global average, exhibited trends generally consistent with that of precipitation. Global total (natural) groundwater storage decreased over the past 5-7 decades with modeled rates ranging from 0.01 to 2.18 mm year-1. This large range can be attributed in part to groundwater’s low frequency (inter-decadal) variability, which complicates identification of real long-term trends even within a 50+ year time series. Results indicate that non-anthropogenic variability in groundwater storage is substantial, making knowledge of it fundamental to quantifying direct human impacts on groundwater storage.

Groundwater storage↗
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↗
NASA NTRS2014

TPSAS-NF1676L-18357-DND

High-spectral resolution measurements (in the IR region) from space born infrared sensors can provide sufficient spectral information about cloud properties like cloud top pressure, cloud particle size and cloud optical thickness. Empirical orthogonal function (EOF) based regression scheme was used to retrieve cloud properties along with atmospheric state profiles. A channel based physical inversion scheme was also tried on cloud-contaminated radiances to improve retrieval accuracy achieved by regression scheme. Recently developed principle component (PC) based physical retrieval algorithm utilized the EOF transformation to compress both hyperspectral radiance spectra and atmospherical profiles into PC scores to greatly boost the computational speed while maintaining the information content. The accuracy of the (PC) based radiation transfer model (PCRTM) to simulate top-of the atmosphere(TOA) high-spectral resolution radiance spectra for both clear and cloudy sky conditions has been well tested by comparing with the other line-by-line radiation transfer model and using the observation data. The application of the PCRTM based physical inversion method on simulated IASI radiance under 'cloudy' conditions will be a good integrity test for the new algorithm.

X Liu↗
NASA NTRS2019

Long-Term, Non-Anthropogenic Groundwater Storage Changes Simulated by Three Global-Scale Hydrological Models

This study examined long-term, natural (i.e., excluding anthropogenic impacts) variability of groundwater storage worldwide. Groundwater storage changes were estimated by forcing three global-scale hydrological models with three 50+ year meteorological datasets. evaluation using in situ groundwater observations from the U.S. and terrestrial water storage derived from the Gravity Recovery and Climate experiment (GRACe) satellites showed that these models reasonably represented inter-annual variability of water storage, as indicated by correlations greater than 0.5 in most regions. Empirical orthogonal function analysis revealed influences of the El Niño Southern Oscillation (ENSO) on global groundwater storage. simulated groundwater storage, including its global average, exhibited trends generally consistent with that of precipitation. Global total (natural) groundwater storage decreased over the past 5–7 decades with modeled rates ranging from 0.01 to 2.18 mm year−1. this large range can be attributed in part to groundwater’s low frequency (inter-decadal) variability, which complicates identification of real long-term trends even within a 50+ year time series. Results indicate that non-anthropogenic variability in groundwater storage is substantial, making knowledge of it fundamental to quantifying direct human impacts on groundwater storage.

Bailing Li↗
NASA NTRS2020

Interannual Variations of TOA Albedo over the Arctic, Antarctic and Tibetan Plateau in 2000–2019

Recent changes in Earth’s climate system have significantly affected the radiation budget and its year-to-year variations at top of the atmosphere (TOA). Observing high-latitude TOA fluxes is still challenging from space, because spatial inhomogeneity of surface/atmospheric radiative processes and spectral variability can reflect sunlight very differently. In this study we analyze the 20-year TOA flux and albedo data from CERES and MISR over the Arctic, the Antarctic, and Tibetan Plateau (TP), and found overall great consistency in the TOA albedo trend and interannual variations. The observations reveal a lagged correlation between the Arctic and subarctic albedo fluctuations. The observed year-to-year variations are further used to evaluate the reanalysis data, which exhibit substantial shortcomings in representing the polar TOA flux variability. The observed Arctic flux variations are highly correlated with cloud fraction (CF), except in the regions where CF > 90% or where the surface is covered by ice. An empirical orthogonal function (EOF) analysis shows that the first five EOFs can account for ~50% of the Arctic TOA variance, whereas the correlation with climate indices suggests that Sea Ice Extent (SIE), North Atlantic Oscillation (NAO) and 55°N–65°N cloudiness are the most influential processes in driving the TOA flux variabilities.

interannual variations↗
NASA NTRS2020

Decadal Changes in the Interannual Variability of Heat Waves in East Asia Caused by Atmospheric Teleconnection Changes

The heat wave in East Asia is examined by using empirical orthogonal function analysis to isolate dominant heat-wave patterns in the ground-based temperature observations over the Korean Peninsula and China and related large-scale atmospheric circulations obtained from the National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis 1 during 1973–2012. This study focuses particularly on the interannual variability of heat waves and its decadal change. The analysis identifies two major atmospheric teleconnection patterns playing an important role in developing typical heat-wave patterns in East Asia—the Scandinavian (SCAND) and the circumglobal teleconnection (CGT) patterns, which exhibit a significant decadal change in the interannual variability in the mid-1990s. Before the mid-1990s, heat-wave occurrence was closely related to the CGT pattern, whereas the SCAND pattern is more crucial to explain heat-wave variability in the recent period. The stationary wave model experiments suggest an intensification of the SCAND pattern in the recent period driven by an increase in land–atmosphere interaction over Eurasia and decadal change in the dominant heat-wave patterns in East Asia.

heat wave patterns↗
NASA NTRS2021

Reduced El NiñO Variability in the Mid-Pliocene According to the Pliomip2 Ensemble

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period during which atmospheric CO 2 levels were similar to recent historical values (∼400 ppm). Several proxy reconstructions for the mid-Pliocene show highly reduced zonal sea surface temperature (SST) gradients in the tropical Pacific Ocean, indicating an El Niño-like mean state. However, past modelling studies do not show these highly reduced gradients. Efforts to understand mid-Pliocene climate dynamics have led to the Pliocene Model Intercomparison Project (PlioMIP). Results from the first phase (PlioMIP1) showed clear El Niño variability (albeit significantly reduced) and did not show the greatly reduced time-mean zonal SST gradient suggested by some of the proxies. In this work, we study El Niño–Southern Oscillation (ENSO) variability in the PlioMIP2 ensemble, which consists of additional global coupled climate models and updated boundary conditions compared to PlioMIP1. We quantify ENSO amplitude, period, spatial structure and “flavour”, as well as the tropical Pacific annual mean state in mid-Pliocene and pre-industrial simulations. Results show a reduced ENSO amplitude in the model-ensemble mean (−24 %) with respect to the pre-industrial, with 15 out of 17 individual models showing such a reduction. Furthermore, the spectral power of this variability considerably decreases in the 3–4-year band. The spatial structure of the dominant empirical orthogonal function shows no particular change in the patterns of tropical Pacific variability in the model-ensemble mean, compared to the pre-industrial. Although the time-mean zonal SST gradient in the equatorial Pacific decreases for 14 out of 17 models (0.2 ∘C reduction in the ensemble mean), there does not seem to be a correlation with the decrease in ENSO amplitude. The models showing the most “El Niño-like” mean state changes show a similar ENSO amplitude to that in the pre-industrial reference, while models showing more “La Niña-like” mean state changes generally show a large reduction in ENSO variability. The PlioMIP2 results show a reasonable agreement with both time-mean proxies indicating a reduced zonal SST gradient and reconstructions indicating a reduced, or similar, ENSO variability.

mid-Pliocene warm period↗
NASA NTRS2023

Operational Forecasting Inundation Extents using REOF analysis (FIER) over Lower Mekong and its Economic Impact on Agriculture

In the Lower Mekong River Basin floodplains, rice cultivation is highly crucial for regional and global food security. However, prolonged flooding can pose damage to rice cultivation and other socio-economic aspects. Yet, there is no rapid operational inundation forecasting system that can help decision-makers proactively mitigate flood damages. Here, we integrated the so-called Forecasting Inundation Extents using Rotated empirical orthogonal function analysis (FIER) framework with an altimetry-based operational Mekong River level forecasting system and built an operational web application, FIER-Mekong, (https://fier-mekong.streamlit.app/) that generates daily skillful forecasted inundation extents (>70% of critical success index) and depths in about 3 and 30 s, respectively, with up to 18-day lead times. One of its applications, predicting flood-induced rice economic losses, is also presented. Had FIER-Mekong being adopted, we estimated that the rice damages, up to 87 and 53 million US dollars during the 2020 and 2021 harvest time, respectively, could have been avoided.

Chi-Hung Chang↗
NASA NTRS1998

Radiance Covariance and Climate Models

Spectral Empirical Orhtogonal Functions (EOFs) derived from the covariance of satellite radiance spectra may be interpreted in terms of the vertical distribution of the covariance of temperature, water vapor, and clouds. The purpose of the investigation is to demonstrate the important constraints that resolved spectral radiances can place upon climate models.

Spectral Empirical Orthogonal Functions EOFs radia↗
NASA NTRS1999

Teleconnection of the 1997 El Nino Observed by Spaceborne Sensors and the Dacadal Anomalies in the Northeast Pacific

Liu et al.[1998] (hereafter referred as LTH), superimposed wind velocity anomalies observed by the NASA Scatterometer (NSCAT) on the map of sea surface temperature (SST) anomalies observed by the Advanced Very High Resolution Radiometer (AVHRR) in the Pacific at the end of May 1997, and illustrated that the three regions of anomalous warming in the North Pacific Ocean are related to wind anomalies through different mechanisms.

global change Empirical Orthogonal Function (EOF) ↗
NASA NTRS1992

Wall functions for the kappa-epsilon turbulence model in generalized nonorthogonal curvilinear coordinates

A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.

Sondak, D. L.↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
NASA NTRS2003

On the Confidence Limit of Hilbert Spectrum

Confidence limit is a routine requirement for Fourier spectral analysis. But this confidence limit is established based on ergodic theory: For stationary process, temporal average equals the ensemble average. Therefore, one can divide the data into n-sections and treat each section as independent realization. Most natural processes in general, and climate data in particular, are not stationary; therefore, there is a need for the Hilbert Spectral analysis for such processes. Here ergodic theory is no longer applicable. We propose to use various adjustable parameters in the shifting processes of the Empirical Mode Decomposition (EMD) method to obtain an ensemble of Intrinsic Mode Function 0 sets. Based on such an ensemble, we introduce a statistical measure in. a form of confidence limits for the Intrinsic Mode Functions, and consequently, the Hilbert spectra. The criterion of selecting the various adjustable parameters is based on the orthogonality test of the resulting M F sets. Length-of-day data from 1962 to 2001 will be used to illustrate this new approach. Its implication in climate data analysis will also be discussed.

Huang, Norden↗
NASA NTRS2005

On the Hilbert-Huang Transform Theoretical Developments

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.

Kizhner, Semion↗
NASA NTRS2006

On Certain Theoretical Developments Underlying the Hilbert-Huang Transform

One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,

Kizhner, Semion↗
NASA NTRS1977

An empirical model for ocean radar backscatter and its application in inversion routine to eliminate wind speed and direction effects

Several regression models were tested to explain the wind direction dependence of the 1975 JONSWAP (Joint North Sea Wave Project) scatterometer data. The models consider the radar backscatter as a harmonic function of wind direction. The constant term accounts for the major effect of wind speed and the sinusoidal terms for the effects of direction. The fundamental accounts for the difference in upwind and downwind returns, while the second harmonic explains the upwind-crosswind difference. It is shown that a second harmonic model appears to adequately explain the angular variation. A simple inversion technique, which uses two orthogonal scattering measurements, is also described which eliminates the effect of wind speed and direction. Vertical polarization was shown to be more effective in determining both wind speed and direction than horizontal polarization.

Dome, G. J.↗
NASA NTRS2000

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.↗
NASA NTRSDate not supplied

Mapping Aerosol Lidar Ratios Over Ocean for CALIPSO Using Modis AOD-Constrained Retrievals and A Global Aerosol Model

The CALIPSO aerosol algorithms currently assign one lidar ratio value for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms by developing regional and seasonal lidar ratio climatologies. In this study, aerosol lidar ratios are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) datasets. The analysis is subsampled for only those profiles that are cloud-free and contain one aerosol type (based on CALIOP feature classification). In addition, the CALIOP profiles are collocated with aerosol volume fractions simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. In this talk, the twelve-year (2006-2017) mean spatial distributions of inferred aerosol lidar ratios for CALIOP-classified marine aerosols and the corresponding modeled sea salt volume fractions (SSVF) will be shown. Model-assisted climatological lidar ratio maps on a seasonal scale will also be presented, developed from the empirical relationship found between the modeled SSVF and lidar ratios. A comparison with past studies will be provided as well as results of a sensitivity study regarding the variability of retrieved lidar ratios as a function of horizontal averaging resolution. While the majority of this talk will focus on lidar ratios for CALIOP-classified marine aerosols, preliminary results will be shown for other aerosol types over ocean, such as dust and elevated smoke.

Travis Toth↗