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 19 records

Nonlinear Unsteady Aerodynamic Modeling Using Empirical Orthogonal Functions

Empirical orthogonal function modeling is explained and applied to identify compact discrete-time nonlinear unsteady aerodynamic models from data generated by an unsteady three-dimensional compressible Navier-Stokes flow solver for an airfoil undergoing various pitching motions. Model structures, model parameter estimates, and model parameter uncertainty estimates for nondimensional lift, drag, and pitching moment coefficient models were determined autonomously and directly from the data. Prediction tests using data that were not used in the modeling process showed that the identified models exhibited excellent prediction capability, which is a strong indicator of an accurate model.

empirical↗

Sampling errors in the estimation of empirical orthogonal functions

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

North, G. R.↗

Challenges and alternatives to empirical orthogonal functions for earth system data

Empirical orthogonal functions (EOFs) applied to gridded Earth system data enables users to diagnose modes of variability with relative ease. Yet, many challenges to interpretation exist such that they must be used with awareness and intention when applied to gridded climate data, especially with large ensembles. Utilizing data from two different Earth system modelling large ensemble frameworks, the Energy Exoscale Earth System Model and the Community Earth System Model, as well as reanalysis data, common EOF pitfalls are summarized and discussed. Challenges include erroneous mode swapping, sign flipping, and the temporal variability of the centers of action. For modes of variability with similar contribution to variance, mode swapping is not uncommon. Sign flipping can occur with almost any mode where the pattern is correct, but the sign is arbitrary. Although the variability of the center of action is not necessarily problematic, it potentially complicates interpretation over multi-century timescales. A wide variety of alternative methods to EOFs exist, but fitness-for-purpose must be evaluated. Additionally, illustrations of alternative methods and examples of proper use are provided. Alternative methods fit into three categories: EOF variants, linear methods, and multilinear methods.

54 ENVIRONMENTAL SCIENCES↗

Empirical orthogonal functions and normal modes

An attempt to provide physical insight into the empirical orthogonal function (EOF) representation of data fields by the study of fields generated by linear stochastic models is presented in this paper. In a large class of these models, the EOFs at individual Fourier frequencies coincide with the orthogonal mechanical modes of the system - provided they exist. The precise mathematical criteria for this coincidence are derived and a physical interpretation is provided. A scheme possibly useful in forecasting is formally constructed for representing any stochastic field by a linear Hermitian model forced by noise.

North, G. R.↗

Empirical orthogonal functions and multiple flow regimes in the Southern Hemisphere winter

A new empirical orthogonal function (EOF) analysis of winter 500 mb geopotential height anomalies in the Southern Hemisphere is presented. An earlier EOF analysis prefiltered the anomalies to exclude wavenumbers 5 and higher; the present analysis does not. The different preprocessing of data affects the results. All three distinct planetary flow regimes identified in the winter circulation of the Southern Hemisphere by a pattern correlation method are captured by the new set of EOFs; only two of those regimes were captured by the earlier set. The new results, therefore, lend further support to the idea that EOFs point to distinct planetary flow regimes.

Farrara, John D.↗

Empirical orthogonal function analysis of cloud-containing coastal zone color scanner images of northeastern North American coastal waters

Empirical-orthogonal-function (EOF) analyses were carried out on 36 images of the Mid-Atlantic Bight and the Gulf of Maine, obtained by the CZCS aboard Nimbus 7 for the time period from February 28 through July 9, 1979, with the purpose of determining pigment concentrations in coastal waters. The EOF procedure was modified so as to include images with significant portions of data missing due to cloud obstruction, making it possible to estimate pigment values in areas beneath clouds. The results of image analyses explained observed variances in pigment concentrations and showed a south-to-north pattern corresponding to an April Mid-Atlantic Bight bloom and a June bloom over Nantucket Shoals and Platts Bank.

Eslinger, David L.↗

Investigation of seasonal variability of the wind stress curl over the North Atlantic Ocean by means of empirical orthogonal function analysis

The seasonal variability of the wind stress curl over the North Atlantic is investigated by means of empirical orthogonal function (EOF) analysis. The curl field is calculated from 1 year of First Global GARP Experiment wind data. It was found that 44 percent of the variability is contained in four significant eigenvectors. Their spatial patterns are characterized by basin-sized oscillations with larger amplitude to the north of 40 deg N. Their associated time series coefficients have the highest amplitude during the winter and show a tendency toward a white frequency spectrum which nevertheless exhibits noticeable peaks or gaps at certain frequencies. Physically, the first EOF is seen as the seasonal fluctuations of the mean wind stress curl pattern. Five other eigenvectors are also found to be above the noise level, but they account for only a smaller percentage of variability (19 percent). They are characterized by smaller spatial scales than the basin size. Their time series coefficients show a whiter frequency spectrum.

Barnier, B.↗

A High-resolution Model of Field-aligned Currents Through Empirical Orthogonal Functions Analysis (MFACE)

Ten years of CHAMP magnetic field measurements are integrated into MFACE, a model of field-aligned currents (FACs) using empirical orthogonal functions (EOFs). EOF1 gives the basic Region-1/Region-2 pattern varying mainly with the interplanetary magnetic field Bz component. EOF2 captures separately the cusp current signature and By-related variability. Compared to existing models, MFACE yields significantly better spatial resolution, reproduces typically observed FAC thickness and intensity, improves on the magnetic local time (MLT) distribution, and gives the seasonal dependence of FAC latitudes and the NBZ current signature. MFACE further reveals systematic dependences on By, including 1) Region-1/Region-2 topology modifications around noon; 2) imbalance between upward and downward maximum current density; 3) MLT location of the Harang discontinuity. Furthermore, our procedure allows quantifying response times of FACs to solar wind driving at the bow shock nose: we obtain 20 minutes and 35-40 minutes lags for the FAC density and latitude, respectively.

CHAMP↗

A Rotation Procedure to Improve Seasonally Varying Empirical Orthogonal Function Bases for MJO Indices

Various indices have been defined to characterize the phase and amplitude of the Madden-Julian oscillation (MJO). One widely used index is the Outgoing Longwave Radiation (OLR) based MJO index (OMI), which is calculated using the spatial pattern of 30-96-day eastward-filtered OLR. The EOFs used to calculate the OMI in observations are prone to degeneracy and exhibit oscillations on the order of 10-20 days, despite initial filtering of the OLR. We propose a simple modification to the OMI that involves aligning the EOFs between neighboring days while retaining the spatial pattern described by the EOFs. This rotation method is implemented as a postprocessing procedure of the current OMI calculation and cleanly removes the spurious oscillations and degeneracy issues seen in the standard method. A similar rotation procedure can be implemented in calculations of other MJO indices.

Sarah Weidman↗

The structure and evolution of seasonal wind anomalies over the near-equatorial eastern Indian and western Pacific Oceans

The longitude-height-time structure and evolution of near-equatorial wind variability over the eastern Indian and western Pacific Oceans are studied using data obtained from a network of eight radiosonde stations extending from southern India to the central Pacific Ocean. The seasonal zonal wind anomalies observed at the cross section beween Trivandrum and Majuro stations are analyzed using an empirical orthogonal function. The Walker Circulation fluctuations are described in terms of standing oscillations in the longitude-height plane, and it is determined that Southern Oscillation propagation anomaly best represents the wind fluctuations. A complex empirical orthogonal function (CEOF) analysis and an El Nino compositing methodology are applied to the seasonal zonal wind anomalies. It is determined that the composite El Nino anomalies correspond to the spatial structure and temporal evolution of anomalies implied by the CEOF analysis.

Gutzler, David S.↗

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗

Geospatial coherence of surface-atmosphere fluxes in the upper Great Lakes region

Surface-atmosphere fluxes are known to vary at multiple time scales, but uncertainty is high as to how fluxes change spatially within regions. With an increase in the number of eddy covariance towers, we are now able to examine the geospatial coherence of ecosystem fluxes, using time-series correlation. Eighteen sites from Michigan and Wisconsin were used in this study, ranging from 100 m to 600 km apart. Surface-atmosphere fluxes from a six-month period were used to quantify spatial coherence on a pair-wise basis. Using geospatial statistics, carbon and sensible heat (H) fluxes were found to be 95% correlated directly outside of their flux footprint and 56% correlated up to a distance of ~35 km. Latent (LE) and momentum (τ) fluxes were less correlated, 83% directly outside of their flux footprint and 40% at a distance of ~130 km albeit, at a much larger spatial distance than for the carbon and sensible heat fluxes. All fluxes showed strong spectral resonance at diel and seasonal timescales, with 1-, 2- and 3-month periods being common modes of variability among H, LE, and τ fluxes. Results based on Empirical Orthogonal Function show distinct transitions of net ecosystem exchange from fall to winter before photosynthesis or respiration while H and τ do not exhibit coherent trends. Here, this work demonstrates the potential of quantifying geospatial coherence of surface-atmosphere fluxes in the Midwestern United States, with the ability to predict fluxes beyond the spatial limit of a single flux tower footprint. Ultimately, expanding the flux measurements to larger scales would allow better spatial scaling of terrestrial surface-atmosphere fluxes between tower footprint and modeling or remote sensing scales.

54 ENVIRONMENTAL SCIENCES↗

An uncertainty-focused database approach to extract spatiotemporal trends from qualitative and discontinuous lake-status histories

Changes in lake status are often interpreted as palaeoclimate indicators due to their dependence on precipitation and evaporation. The Global Lake Status Database (GLSDB) has since long provided a standardised synopsis of qualitative lake status over the last 30,000 14C years. Potential sources of uncertainty however are not recorded in the GLSDB. Here we present an updated and improved relational-database framework that incorporates uncertainty in both chronology and the interpretation of palaeoenvironmental data. The database uses peer-reviewed palaeolimnological studies to produce a consensus on qualitative lake-status histories, whose chronologies are revised and standardized through the recalibration of radiocarbon dates and the application of Bayesian age-depth modelling for stratigraphic archives. Quantitative information on absolute water-level elevation is preserved if available from geomorphological sources. We also propose a new probabilistic analytical framework that accounts for these uncertainties to reconstruct synoptic, integrated environmental signals. The process is based on a Monte Carlo algorithm that iteratively samples individual lake-status histories within the limits of their uncertainties to produce many possible scenarios. We then use Recursively-Subtracted Empirical Orthogonal Function analysis to extract dominant patterns of lake-status variability from these scenarios. As a proof of concept, we apply this framework to 67 sites in eastern and southern Africa whose lake-status histories cover part of the late Pleistocene and/or Holocene. We show that, despite the sometimes large temporal and interpretation uncertainties, and the inclusion of highly discontinuous lake-status time series, identifying the major known millennial-scale climatic phases during the last 20,000 years is possible. Our framework was also able to identify an antiphased response between the lake basins in eastern and interior southern Africa to these changes. Here, we propose that our new database and methodology framework serves as a template for efficient lake-status data synthesis, encourages the incorporation of lake-status data in palaeoclimate syntheses, and expands the possibilities for the use of such data in the evaluation of climate models.

58 GEOSCIENCES↗

Benchmarking performance changes in the simulation of extratropical modes of variability across CMIP generations

In this work, we evaluate extratropical modes of variability in the three most recent phases of the Coupled Model Intercomparison Project (CMIP3, 5, and 6) to gauge improvement of climate models over time. A suite of high-level metrics is employed to objectively evaluate how well climate models simulate the qobserved Northern Annular Mode (NAM), North Atlantic Oscillation (NAO), Pacific North America pattern (PNA), Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO), North Pacific Oscillation (NPO), and North Pacific Gyre Oscillation (NPGO). We apply a Common Basis Function (CBF) approach that projects model anomalies onto observed empirical orthogonal functions (EOFs), together with the traditional EOF approach, to CMIP Historical and AMIP models. We find simulated spatial patterns of those modes have been significantly improved in the newer models, although the skill improvement is sensitive to the mode and season considered. We identify some potential contributions to the pattern improvement of certain modes (e.g., the Southern Hemisphere jet and High-top vertical coordinate), however the performance changes are likely attributed to gradual improvement of the base climate and multiple relevant processes. Less performance improvement is evident in the mode amplitude of these modes and systematic overestimation of the mode amplitude in spring remains in the newer climate models. We find that the post-dominant season amplitude errors in atmospheric modes are not limited to coupled runs but are often already evident in AMIP simulations. This suggests that rectifying the egregious post-dominant season amplitude errors found in many models can be addressed in an atmospheric-only framework making it more tractable to address in the model development process.

54 ENVIRONMENTAL SCIENCES↗

Tropical synoptic interpretation of inter-channel correlations from Tiros N

Synoptic variability contained in a set of Tiros N channel brightness temperatures was analyzed for a nine-day period during January 1979, between the dateline and 90 deg W and between 30 deg N and 20 deg S. In the procedure, the spatial coherence of each channel is first described in terms of functions of interchannel differences and geographic and synoptic variations. Second, horizontal and vertical pattern of individual channels are examined through empirical orthogonal functional decomposition. Finally, response surface regression models are developed to predict missing Tiros N channels from available Tiros N and NOAA-5 data. The percentage variance was predicted for each channel from the Tiros N microwave and stratospheric channels. Generally speaking, about 80 percent of the variance in each channel can be predicted for the Tiros model, excepting the two tropospheric moisture channels, where prediction was less successful.

Mcguirk, J. P.↗