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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Principal component analysis of solar flares in the soft X-ray flux

The paper considers principal component analysis of solar flares in the soft X-ray flux, a technique for extracting the salient features from a mass of data. The method applies particularly to the analysis of nonstationary ensembles, and its computations require the evaluation of eigenvalues of matrices. The Eispack matrix eigen system routines were used to analyze full-disk proportional-counter data from the X-ray event analyzer which was part of the Skylab experiment. Empirical orthogonal functions were derived for events in the soft X-ray spectrum between 2.5 and 20 A during different time periods, indicating that about 90% of the cumulative power of each analyzed flare is contained in the largest eigenvector. The first two largest eigenvectors are sufficient for an empirical curve fit through the raw data and a characterization of solar flares in the soft X-ray flux, and power spectra of two largest eigenvectors reveal a reported periodicity of about 5 min.

Teuber, D. L.↗
NASA NTRS2007

An EOF Iteration Approach for Obtaining Homogeneous Radiative Fluxes from Satellites Observations

Conventional observations of climate parameters are sparse in space and/or in time and the representativeness of such information needs to be optimized. Observations from satellites provide improved spatial coverage than point observations however they pose new challenges for obtaining homogeneous coverage. Surface radiative fluxes, the forcing functions of the hydrologic cycle and biogeophysical processes, are now becoming available from global scale satellite observations. They are derived from independent satellite platforms and sensors that differ in temporal and spatial resolution and in the size of the footprint from which information is derived. Data gaps, degraded spatial resolution near boundaries of geostationary satellites, and different viewing geometries in areas of satellite overlap, could result in biased estimates of radiative fluxes. In this study, discussed will be issues related to the sources of inhomogeneity in surface radiative fluxes as derived from satellites; development of an approach to obtain homogeneous data sets; and application of the methodology to the widely used International Satellite Cloud Climatology Project (ISCCP) data that currently serve as a source of information for deriving estimates of surface and top of the atmosphere radiative fluxes. Introduced is an Empirical Orthogonal Function (EOF) iteration scheme for homogenizing the fluxes. The scheme is evaluated in several ways including comparison of the inferred radiative fluxes against ground observations, both before and after the EOF approach is applied. On the average, the latter reduces the rms error by about 2-3 W/m2.

Zhang, Banglin↗
DOE OSTI2024

Concurrent Inter-Model Spread of Boreal Winter Westerly Jet Meridional Positions Between the Northern and Southern Hemispheres in CMIP6 Models

Here, this study investigates the inter-model spread of climatological extratropical westerly jets in boreal winter, using the historical simulation of 52 Coupled Model Intercomparison Project phase 6 (CMIP6) models from 1851 to 2014. The results show that there is a substantial spread in the latitude of the upper-tropospheric westerly jet across models, characterised by large inter-model standard deviations to both the poleward and equatorward sides of the jet axis, although the multi-model ensemble mean (MME) performs well in simulating meridional position of westerly jets. Furthermore, we detect the consistency of inter-model jet position spread between the Northern and Southern Hemispheres, based on the inter-model empirical orthogonal function (EOF) decomposition and correlation of regional-averaged zonal winds. Specifically, the models that simulate the westerly jets poleward/equatorward relative to the MME position in one hemisphere also tend to simulate the jets poleward/equatorward in the other hemisphere. Accordingly, we define a global jet spread index to depict the concurrence of jet shift in the two hemispheres. The results of inter-model regression analyses based on this index indicate that the models positioning the jets poleward than the MME tend to simulate a wider Hadley Cell, a poleward-shifted Ferrel Cell in the Southern Hemisphere, enhanced precipitation in the subtropics and suppressed precipitation in the tropics, and warmer sea surface temperatures in the subtropics and mid-latitudes. The present results suggest that improving the simulation of jet positions in climate models requires a comprehensive consideration of thermal states in the tropics and subtropics/mid latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Global teleconnections influencing large-scale drought in the United States using SVDI

Understanding recent large-scale drought patterns and the mechanisms producing extreme drought events is vital for future drought forecasts and understanding future drought risks. Increasingly, vapor pressure deficit (VPD) has been used as an important measure of evaporative demand and proxy for drought detection. In this study, VPD is used to calculate the new Standardized VPD Drought Index (SVDI) with NASA North American Land Data Assimilation System (NLDAS) data. Previous studies have shown that SVDI accurately identifies the timing and magnitude short-term droughts in the United States (U.S). In the present study, SVDI is now used to identify large-scale drought patterns between 1980 and 2021 and drought variability driven by selected global teleconnections originating in the Pacific and Atlantic Oceans. Spatial drought characteristics were extracted from SVDI using empirical orthogonal function (EOF) analysis. Then a k-means clustering algorithm was applied to both EOF principal components and primary teleconnections, including the El Nino-Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO) to identify drought events driven by the Pacific Ocean. Results show that the SVDI is useful in evaluating large-scale drought variability in the U.S. related to global teleconnections, and that mechanisms influencing summer drought patterns in the Western and Southwestern U.S. are driven by a tropical-extratropical interactions originating in the equatorial Pacific Ocean related to ENSO dynamics with interdecadal variability modulated by PDO. The large-scale droughts in the Central and Southern U.S., like those in 2011 and 2012, on the other hand, are driven by the North Pacific Ocean warm pool during a strong negative PDO, which subsequently influenced variability in the Bermuda-Azores High in the Atlantic Ocean. In summer 2011, the Bermuda-Azores High weakened, reducing the onshore winds and moisture transport along the eastern Gulf of Mexico and contributing to ongoing drought in the region. The Northern Pacific and Atlantic Ocean sea surface temperatures (SSTs) have increased between 1980 and 2021. In conclusion, as SSTs continue to rise in the Northern Pacific Ocean, one consequence of the coupled North Pacific warm pool and atmospheric dynamics, is to increase summer drought variability over a large region in the southern and midwestern U.S. under global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Spatial Patterns and Trends of Summertime Low Cloudiness for the Pacific Northwest, 1996–2017

Summertime low clouds are common in the Pacific Northwest (PNW), but spatiotemporal patterns have not been characterized. We show the first maps of low cloudiness for the western PNW and North Pacific Ocean using a 22-year satellite-derived record of monthly mean low cloudiness frequency for May through September and supplemented by airport cloud base height observations. Domain-wide cloudiness peaks in midsummer and is strongest over the Pacific. Empirical orthogonal function (EOF) analysis identified four distinct PNW spatiotemporal modes: oceanic, terrestrial highlands, coastal, and northern coastal. There is a statistically significant trend over the 22-year record toward reduced low cloudiness in the terrestrial highlands mode, with strongest declines in May and June; however, this decline is not matched in the corresponding airport records. The coastal mode is partly constrained from moving inland by topographic relief and migrates southward in late summer, retaining higher late-season low cloud frequency than the other areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Impact of Inherited Morphology on Sandbar Migration During Mild Wave Seasons

Sandbars are ubiquitous morphologic features found in the nearshore environment throughout the world, yet predictive capabilities of their evolution remain limited. In order to provide new insights on the relevant processes controlling sandbar morphodynamics, this study uses a 41-year record of 637 monthly cross-shore profiles from Duck, North Carolina, USA, to derive complex empirical orthogonal functions representative of the two dominant modes of sandbar migrations: offshore and onshore propagation. Interference of these two modes produces commonly observed sandbar states. While mild wave energies are traditionally assumed to drive onshore sandbar migration, the offshore mode is repeatedly seen to dominate sandbar migration in mild wave seasons following anomalously high late-winter wave energy and when the inherited morphology is composed of a single bar or terrace, as opposed to a more common two-bar state. A data-derived conceptual model is presented synthesizing the effect of antecedent morphology and wave climate chronology on interannual trends in net offshore sandbar migration.

58 GEOSCIENCES↗
DOE OSTI2022

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Propagation and Periodicity of Mars's Northern Annular Mode Modulates the Dust Cycle

Abstract We document the propagation of annular modes—zonally symmetric patterns of variability—in Mars's atmosphere using a reanalysis dataset. Mars's Northern Annular Mode (MNAM) sees anomalies of zonal‐mean zonal wind emerge near the subtropics and migrate poleward with a period of 150 days, similarly to Earth's Southern Annular Mode. The mechanism of propagation involves the interaction of the two leading empirical orthogonal functions that define the MNAM. Moreover, the propagation encourages alternating bands of surface wind stress to migrate polewards with a 150‐day period. In addition, a 150‐day periodicity in anomalous column dust optical depth most likely emerges in response to extrema of the MNAM. The combination of the impact of the MNAM's internally forced periodicity on the surface wind stress and the seasonal cycle may contribute to the inter‐annual variability of global dust events, as suggested by a Monte Carlo estimate that correctly approximates the observed incidence of global dust events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗
DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

OpenCRUMS USA: An Open Machine Learning Framework for Characterizing Variability in Aerosol Reanalysis Data

Advances in artificial intelligence (AI) have called for exploring how these techniques can be used for exploring patterns in large climate datasets. To that regard, the U.S. Department of Energy AI for Earth System Predictability (AI4ESP) supported a pilot initiative called the Open Classification of Regimes in the Southeast USA (OpenCRUMS USA) project to explore how AI can be used to characterize modes of spatial variability in large climate datasets. For this study, we focus on comparing two methods for characterizing the modes of spatial variability of surface aerosol concentration over the Houston region: empirical orthogonal functions (EOFs) and layerwise relevance propagation (LRP) applied to a convolutional neural network (CNN) classifier. We show that EOF analysis typically attributes spatial variability modes that span all of southeast Texas, prohibiting the attribution of spatial variability to localized regions. However, using LRP on the CNN classifier resolves the explanatory parameters at a finer spatial resolution than EOFs. This allows for the attribution of the spatial variability of surface aerosols to local regions of organic carbon which was not possible using EOFs. In addition, the LRP analysis also suggests that synoptic-scale transport of dust is most prevalent during anticyclonic and pretrough synoptic conditions as categorized by self-organizing maps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of African Easterly Waves in a Global Climate Model

Abstract African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Beyond PCA: Additional Dimension Reduction Techniques to Consider in the Development of Climate Fingerprints

Abstract Dimension reduction techniques are an essential part of the climate analyst’s toolkit. Due to the enormous scale of climate data, dimension reduction methods are used to identify major patterns of variability within climate dynamics, to create compelling and informative visualizations, and to quantify major named modes such as El Niño–Southern Oscillation. Principal components analysis (PCA), also known as the method of empirical orthogonal functions (EOFs), is the most commonly used form of dimension reduction, characterized by a remarkable confluence of attractive mathematical, statistical, and computational properties. Despite its ubiquity, PCA suffers from several difficulties relevant to climate science: high computational burden with large datasets, decreased statistical accuracy in high dimensions, and difficulties comparing across multiple datasets. In this paper, we introduce several variants of PCA that are likely to be of use in climate sciences and address these problems. Specifically, we introduce non-negative , sparse , and tensor PCA and demonstrate how each approach provides superior pattern recognition in climate data. We also discuss approaches to comparing PCA-family results within and across datasets in a domain-relevant manner. We demonstrate these approaches through an analysis of several runs of the E3SM climate model from 1991 to 1995, focusing on the simulated response to the Mt. Pinatubo eruption; our findings are consistent with a recently identified stratospheric warming fingerprint associated with this type of stratospheric aerosol injection.

Weylandt, Michael↗
DOE OSTI2024

Investigating the Relative Contribution from Tropical Indo-Pacific SST to Asian Monsoon Precipitation Variability Using LIM

Abstract A critical issue is determining the factors that control the year-to-year variability in precipitation over southern Asia. In this study, we employ a cyclostationary linear inverse model (CS-LIM) to quantify the relative contribution of tropical Pacific and Indian Ocean sea surface temperature anomalies (SSTAs) to the interannual variability of the Asian monsoon, especially Indian summer monsoon rainfall (ISMR). Through a series of CS-LIM experiments, we isolate the impacts of the direct forcing from Pacific SSTAs, Indian Ocean SSTAs, and their interaction on Asian monsoon rainfall variability. Our results reveal distinct patterns of influence with the direct forcing from the Pacific (Indian) Ocean tending to enhance (reduce) the magnitude of precipitation variability, while the Indo-Pacific interaction acts to strongly damp the variability of Asian monsoon precipitation, especially over India. We further investigate these specific impacts on ISMR by analyzing the relationship between tropical Indo-Pacific SSTAs and the leading three empirical orthogonal functions (EOFs) of ISMR. The results from our CS-LIM experiments indicate that the direct forcing from El Niño–Southern Oscillation (ENSO) enhances the variability of the first and third EOFs, while the Indian Ocean SSTA opposes ENSO’s effects, which is consistent with previous studies. Our new results show that the tropical Indo-Pacific interaction strongly damps ISMR variability, which is due to the ENSO-induced Indian Ocean dipole (IOD) opposing the direct impacts from ENSO on ISMR. Additionally, reduced ENSO amplitude and duration associated with the Indo-Pacific interaction may also contribute to the damping effect on ISMR, but this requires further study to understand the relevant mechanisms.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

A double peak in the seasonality of California's photosynthesis as observed from space

Solar-induced chlorophyll fluorescence (SIF) has been shown to be a powerful proxy for photosynthesis and gross primary productivity (GPP). The recently launched TROPOspheric Monitoring Instrument (TROPOMI) features the required spectral resolution and signal-to-noise ratio to retrieve SIF from space. Here, we present a downscaling method to obtain 500 m spatial resolution SIF over California. We report daily values based on a 14 d window. TROPOMI SIF data show a strong correspondence with daily GPP estimates at AmeriFlux sites across multiple ecosystems in California. We find a linear relationship between SIF and GPP that is largely invariant across ecosystems with an intercept that is not significantly different from zero. Measurements of SIF from TROPOMI agree with MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation indices – the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation index (NIR$_v$) – at annual timescales but indicate different temporal dynamics at monthly and daily timescales. TROPOMI SIF data show a double peak in the seasonality of photosynthesis, a feature that is not present in the MODIS vegetation indices. The different seasonality in the vegetation indices may be due to a clear-sky bias in the vegetation indices, whereas previous work has shown SIF to have a low sensitivity to clouds and to detect the downregulation of photosynthesis even when plants appear green. We further decompose the spatiotemporal patterns in the SIF data based on land cover. The double peak in the seasonality of California's photosynthesis is due to two processes that are out of phase: grasses, chaparral, and oak savanna ecosystems show an April maximum, while evergreen forests peak in June. An empirical orthogonal function (EOF) analysis corroborates the phase offset and spatial patterns driving the double peak. The EOF analysis further indicates that two spatiotemporal patterns explain 84 % of the variability in the SIF data. Results shown here are promising for obtaining global GPP at sub-kilometer spatial scales and identifying the processes driving carbon uptake.

54 ENVIRONMENTAL SCIENCES↗
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 NTRS1980

Eigenvector analysis of some observed and model-generated Climatological fields

An empirical orthogonal function analysis of monthly mean synoptic fields computed with a general circulation model and the corresponding observed climatological patterns is presented. A computer program was used to compute the principal components of the data matrices representing global maps of various observed and model-generated climatological fields. The variables analyzed were: 500mb geopotential height, 850 mb temperature, and sea level pressure. It is concluded that eigenvector analysis of observed and model-generated synoptic climatological fields is a useful diagnostic aid in the evaluation of global model simulation. The maps of any climatological variables are readily decomposed into a few orthogonal components which account for a large fraction of the total spatial variance, and the observed and model-generated patterns may then be compared in terms of these principal components of the original data matrices.

Dennis, M.↗