Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fitting”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Revisiting the Optical Spectrum of the Plutonyl Ion (PuO 2 ) 2+ in 1 M HClO 4

The analysis of the solution absorption spectrum of the plutonyl ion in an aqueous environment was given by Eisenstein and Pryce (E&P) in 1968. In 2011 a new spectrum was published of the (PuO 2 ) 2+ ion in 1 M HClO 4 . We have been provided with the original data of this spectrum and have found in the data a previously unreported low-lying transition at 7385 cm –1 which we have assigned as a magnetic dipole transition. We have fit most of the near-infrared and optical transitions with Gaussian fits and tabulated a new energy level list up to 22,000 cm –1 which mostly agrees with the data of E&P. We assumed a crystal field of D ∞h (only axial symmetry) and utilized the intensity calculations published for the isoelectronic (NpO 2 ) 1+ ion using a complete basis set for the 5f 2 problem including the Coulombic, spin–orbit as well as the crystal field Hamiltonian. Our results differ substantially from those of E&P. Subsequently, we used a truncated Hamiltonian to try to establish the effects of assuming the σ antibonding orbitals are at such high energies that we can ignore their contributions to the lower lying φ and δ orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Nanosecond Structure of Radical Pair Intermediates from High-Frequency Quantum Oscillations: Insight into the Q A •– to Q B Electron Transfer Step in Purple Bacterial Photosynthesis

We demonstrate the validity of our approach to deduce, from the anisotropy of quantum oscillations, the geometry of short-lived radical pair intermediates in photosynthesis. A global fit of a two-dimensional W-band (94 GHz) electron paramagnetic resonance (EPR) experiment provides the same global minimum values for the geometry of the A-side radical pair P 700 •+ A 1A •− in photosystem I (PSI) as observed in a previous Q-band (34 GHz) EPR study, yet with a significantly increased convergence rate of 62%. This demonstrates that the global fit yields the correct radical pair geometry even at Q-band frequencies. With this information, we revisit our previous Q-band study of the cofactor arrangement of P 865 •+ Q A •− , the stabilized charge-separated state in purple bacterial reaction centers (RCs). Analysis of calculated two-dimensional data sets of P 865 •+ Q A •− reveals that the quantum oscillation technique is unaffected by a mirror ambiguity in disordered solids and thus can provide unambiguous solutions for all five Euler angles of the radical pair geometry. This enables us to elucidate the Q A •− to Q B electron transfer step in purple bacterial photosynthesis, the subject of controversial discussions for more than 25 years. Our results show that this electron transfer step involves a gating mechanism requiring a 60° rotation of the headgroup of Q A •− in its binding pocket.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

High Pressure X-ray Diffraction and Equation of State of Hydrazine

Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.

diffraction↗
DOE OSTI2025

Comparing Machine Learning and Physics-Based Nanoparticle Geometry Determinations Using Far-Field Spectral Properties

Anisotropic metal nanostructures exhibit polarization-dependent light scattering, a property which has been widely studied and exploited to determine orientations of subwavelength structures using far-field microscopy. Here we explore the use of variational autoencoders (VAEs) to determine the geometries of gold nanorods (NRs) such as in-plane orientation and aspect ratio under linearly polarized dark-field illumination in an optical microscope. We enforce a shared latent space to connect two VAEs trained separately with polarized dark-field scattering spectra and electron microscopy images and achieve image prediction (shape, orientation, and size) of Au NRs using only polarized dark-field scattering spectra. We determine the geometrical parameters of orientational angle and aspect ratio quantitatively via both our dual-VAE and physics-based analysis on the input scattering spectra. We show that orientational angle prediction by dual-VAE performs well with only a small (~300 particle) training set, yielding a mean absolute error (MAE) of 14.4° and a concordance correlation coefficient (CCC) of 0.95. This performance is only marginally worse than the physics-based cos(2?) fitting approach between the scattering intensity and the polarizing angle, which achieves MAE of 8.78° and CCC of 0.99. Aspect ratio determination is also comparable for the dual-VAE and physics-based fitting comparison (MAE of 0.21 vs. 0.23 and CCC of 0.53 vs. 0.68). Here, this dual encoder-decoder architecture effectively exploits the structure-property relationships of plasmonic nanostructures to construct a cross-modal machine learning (ML) approach, providing a pathway to employ ML approaches to address other structure-property relationships in materials science.

Dark-field scattering↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2025

Self-Diffusivity Measurement of Eutectic F 7 LiNaK with and without Additives Using Quasi-Elastic Neutron Scattering

The atomic scale relaxation dynamics of eutectic F 7 LiNaK (46.5 LiF–11.5 NaF–42 KF mol %, Li-7 enriched) were measured using quasi-elastic neutron scattering (QENS) over a temperature range of 500–750 °C. Here, the effect of adding 0.988 mol % cerium, 0.499 mol % cesium, and 1.21 mol % zirconium individually to the dynamics of F 7 LiNaK was also investigated. The relaxation process in both pure and doped F 7 LiNaK molten salts was fit with a stretched exponential function and the temperature dependence follows an Arrhenius behavior over a wavevector transfer range of 0.4 Å –1 < Q < 0.9 Å –1 . The measured activation energy for self-diffusion is E a = 0.77 ± 0.02 eV/atom for pure molten F 7 LiNaK. The QENS response with additives added to F 7 LiNaK was also fit with a stretched exponential and the associated Arrhenius behavior was characterized with activation energies of E a = 0.88 ± 0.01 eV/atom for zirconium (1.21 mol %), E a = 1.02 ± 0.02 eV/atom for cerium (0.988 mol %), and E a = 0.71 ± 0.03 eV/atom for cesium (0.499 mol %). The measured diffusivities are compared to those simulated with a neural network force field model by Lee et al. [Lee, S.-C. Comparative Studies of the Structural and Transport Properties of Molten Salt FLiNaK Using the Machine-Learned Neural Network and Reparametrized Classical Forcefields. J. Phys. Chem. B 2021, 125(37), 10562–10570].

FLiNaK↗
DOE OSTI2025

Predicting Pulsed-Laser Deposition SrTiO 3 Homoepitaxy Growth Dynamics Using High-Speed Reflection High-Energy Electron Diffraction

Pulsed-laser deposition (PLD) is a powerful technique for growing complex oxides with controlled stoichiometry. To understand growth dynamics therein, it is common to leverage in situ spectroscopies, such as reflection high-energy electron diffraction (RHEED), to monitor surface crystallinity. Most commercial systems rely on video-rate cameras operating at 60-120 Hz that lack sufficient temporal resolution to capture growth dynamics at practical deposition frequencies. Here, a high-speed platform to record in situ dynamics via RHEED at >500 Hz is implemented. An open-source analysis package is designed to fit diffraction spots to 2D Gaussians, allowing single-pulse surface reconstruction kinetics extraction. Using homoepitaxially deposited (001)-oriented SrTiO 3 as a model system, we demonstrate how high-speed RHEED can provide real-time insight into growth processes obscured by slower acquisition systems. By fitting the single-pulse intensity to a set of exponential functions, we observe changes in the characteristic decay time and mechanism correlated to the substrate step width and surface termination. We observe distinct surface effects, with diffraction intensity decaying on lower-energy TiO 2 -terminated surfaces and stabilizing on SrO- or mixed-terminated surfaces. Similarly, using an exponential model, the extracted characteristic time of adatom deposition decreases with increased density of bonding sites associated with mixed termination and narrower step widths. Ultimately, this work shows how increasing RHEED temporal resolution can uncover new insights into growth processes, with practical implications for the design and control of PLD processes. This experimental platform provides new capabilities to enable data-driven machine learning analysis and autonomous control systems to enhance the complexity and fecundity of PLD.

(SrO)↗
DOE OSTI2025

Revisiting the Role of Mobile Ion Size for the Activation Barrier of Charge Transport in Single Ion Conducting Polymers

The development of polymer electrolytes for energy storage and conversion technologies requires a fundamental understanding of the material parameters controlling the energy barriers for ion transport. In glassy polymers, these activation barriers are usually extracted by using Arrhenius procedures. However, our recent studies of single ion conducting polymers reveal that this traditional Arrhenius description provides anomalously small prefactors, an issue that is widely disregarded in the literature. Here, to address it, an alternative approach was introduced for extracting the effective activation barrier E by imposing a physically valid limit (∼10 –13 s) for the ion hopping attempt time. Two consequences of this recent approach are that in polymer electrolytes (i) E is significantly (∼30–40%) lower than the barrier estimated using traditional Arrhenius fits and (ii) E displays significant temperature dependence even in the glassy state. Under these circumstances, we are revisiting the role of ion size, glass transition temperature T g , and dielectric and elastic constants of the polymer matrix on the effective activation barriers for ion transport as extracted using the recently proposed approach and highlight the differences from the picture that previously emerged based on simplistic apparent Arrhenius analysis. To this end, we investigate cation transport in three families of single ion conducting trifluoromethane sulfonimide-based polyanions with varied T g . Our results indicate that E decreases with mobile cation size and is highly sensitive to changes in the dielectric permittivity of the matrix, even for large cations. These insights call for revisions of many earlier results based on apparent Arrhenius fits as the proposed approach can provide more accurate guidance for the design of polymer electrolytes with enhanced ionic conductivity.

energy barrier↗
DOE OSTI2025

A Transferable Force Field for Predicting Adsorption and Diffusion of Water in Cationic Zeolites with Coupled Cluster Accuracy

We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.

Adsorption↗
DOE OSTI2025

Engineered Membrane Vesicle Production via oprF or oprI Deletion Has Distinct Phenotypic Effects in Pseudomonas putida

Membrane vesicle (MV) production is a natural phenomenon in Gram-negative bacteria and represents an emerging synthetic biology tool for the secretion of biomolecules or bioproducts. Manipulation of membrane components has proven successful in enhancing MV production. However, the impact of membrane disruptions on strain fitness and protein composition warrants further investigation for the use of MVs in industrial bioprocesses. Here, we identify and characterize two genetic engineering strategies for inducing hypervesiculation-deletion of genes for the outer membrane porin OprF or the lipoprotein OprI-in the commonly used platform Pseudomonas putida KT2440. Deletion of oprI generated up to a 1.5-fold increase in MVs, larger MVs with a greater proportion of outer membrane proteins, and no significant impact on strain fitness compared to wild type. In contrast, deletion of oprF, relative to wild type, generated up to a 4-fold increase in MVs but diminished growth, permeabilized membranes, and increased cytosolic protein packaging. Both hypervesiculation phenotypes increased nontargeted and MV-targeted mNeonGreen extracellular signal by up to 6-fold, demonstrating vesiculation as a mechanism for protein secretion. Despite increased blebbing of MVs from gene deletions, proteins involved in membrane biosynthesis were not elevated relative to wild type. Overexpression of gpsA, which initiates glycerophospholipid biosynthesis, in the ΔoprF background improved the membrane integrity by 37% and maintained MV formation, highlighting the importance of membrane biosynthesis in restoring the membrane in hypervesiculating strains. Together, this study provides genetic engineering strategies with corresponding phenotypic outcomes toward providing a synthetic biology toolset for MV deployment in P. putida.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization↗
DOE OSTI2025

Confronting Large‐Eddy Simulations With Stereo Camera Data by Means of Reconstructed Hemispheric Cloud Size Distributions

High-resolution hemispheric camera images at a meteorological site in western Germany are used to analyze the multi-dimensional spatial characteristics of continental cumulus cloud fields, and to evaluate Large-Eddy Simulations on this aspect. Traditional non-hemispheric cloud-detecting instruments provide additional reference data. The main model-observation comparison focuses on cloud size distributions (CSDs), employing two methods: (a) directly using three-dimensional model fields, direct CSDs, and (b) using rendered hemispheric images of the model fields as produced by a camera simulator based on path-tracing. In the latter method, both the real and rendered images are used to three-dimensionally reconstruct the cloud fields, yielding hemispheric CSDs. Advantages of hemispheric comparisons over more classic approaches include (a) fair comparisons between model and data, and (b) full use of the enhanced resolutions and hemispheric spatial coverage of the camera imagery. Basic evaluation of the simulations demonstrates good agreement on thermodynamic structure and its diurnal cycle. Cloud heights and cloud cover are intercompared between the model, camera data and other instrumentation, providing insight into their structural differences. A consistent alignment is found between the hemispheric CSDs from both the model and the cameras. Power law fits reveal structurally lower exponents in hemispheric CSDs compared to non-hemispheric CSDs, which particularly caution against directly comparing hemispheric CSDs to non-hemispheric distributions. This result is robust for sample size and fitting method. These findings inform future use of hemispheric camera systems for studying cumulus cloud field morphology and model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impacts of Bulk Microphysics Scheme Structural Choices on Simulations of Rain Initiation Through Drop Coalescence

This study examines how different structural choices in bulk microphysics schemes impact the simulation of warm rain initiation. A single liquid category (SLC) approach prognosing up to four moments of a single drop size distribution (DSD) is compared to the traditional two-category, two-moment approach with separate DSDs for cloud and rain (four total prognostic variables). Different methods for calculating tendencies of the prognostic variables from drop collision-coalescence are also tested: a discretized numerical-integration approach, machine learning via neural networks, lookup tables, and traditional power law fits. Relative to simulations using a bin microphysics model, SLC gives smaller error overall than the two-category approach when numerical integration is used to calculate the collision-coalescence tendencies for both. Replacing the numerical integration with a pre-computed lookup table reduces computational cost with little loss of accuracy. However, using fitted power laws with SLC to represent the collision-coalescence tendencies substantially reduces accuracy and leads to an order of magnitude increase in error. It is also demonstrated that with SLC, reasonably accurate solutions are obtained using only three prognostic moments, while a two-moment SLC scheme leads to substantial error. Overall, both the choice of prognostic moments (e.g., SLC vs. two-category) and method to calculate the collision-coalescence tendencies are important to consider for minimizing errors in bulk schemes. SLC with a sufficiently detailed calculation of the collision-coalescence tendencies provides accurate solutions for a reasonable computational cost, providing a viable alternative to the traditional two-category, two-moment approach for bulk microphysics.

320 (cloud physics and chemistry)↗
DOE OSTI2026

From Points to Planes: A Workflow for Converting Three‐Dimensional Point Cloud Data Into Discrete Fracture Network Flow and Transport Models

We present the Point cLoud Algorithm for NEtwork Extraction of Discrete Fracture Networks (PLANE-DFN), a point cloud–based algorithm for automatic fracture network extraction designed to support discrete fracture network (DFN) modeling workflows. PLANE-DFN segments three-dimensional fracture planes from raw point cloud data using RANdom SAmple Consensus coupled with statistical outlier removal and density-based clustering to isolate individual fracture features. Each candidate plane is constrained against site-specific structural constraints based on strike and dip. After segmentation, each fracture is converted into a 2-D convex polygon suitable for meshing and simulation. The PLANE-DFN algorithm is validated by comparing geometric and flow and transport data against data from dfnWorks simulations with ensembles of plane-fit networks. We find that the flow and transport in plane-fit networks are comparable to dfnWorks-generated networks when realistic network geometry is maintained. The PLANE-DFN algorithm provides an automated and streamlined workflow to transform point clouds of data into DFN network geometry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Genome-scale analysis of interactions between genetic perturbations and natural variation

Interactions between genetic perturbations and segregating loci can cause perturbations to show different phenotypic effects across genetically distinct individuals. To study these interactions on a genome scale in many individuals, we used combinatorial DNA barcode sequencing to measure the fitness effects of 8046 CRISPRi perturbations targeting 1721 distinct genes in 169 yeast cross progeny (or segregants). We identified 460 genes whose perturbation has different effects across segregants. Several factors caused perturbations to show variable effects, including baseline segregant fitness, the mean effect of a perturbation across segregants, and interacting loci. We mapped 234 interacting loci and found four hub loci that interact with many different perturbations. Perturbations that interact with a given hub exhibit similar epistatic relationships with the hub and show enrichment for cellular processes that may mediate these interactions. These results suggest that an individual’s response to perturbations is shaped by a network of perturbation-locus interactions that cannot be measured by approaches that examine perturbations or natural variation alone.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Growth rate as a link between microbial diversity and soil biogeochemistry

The growth rate of a microorganism is a simple yet profound way to quantify its impact on the world. The absolute growth rate of a microbial population reflects rates of resource assimilation, biomass production, and element transformation, some of the many ways that organisms affect Earth’s ecosystems and climate. Microbial fitness in the environment depends on the ability to reproduce quickly when conditions are favorable and adopt a survival physiology when conditions worsen, which cells coordinate by adjusting their relative growth rate. At the population level, relative growth rate is a sensitive metric of fitness, linking survival and reproduction to the ecology and evolution of populations. Techniques combining ‘omics and stable isotope probing enable sensitive measurements of growth rates of microbial assemblages and individual taxa in soil. Microbial ecologists can explore how the growth rates of taxa with known traits and evolutionary histories respond to changes in resource availability, environmental conditions, and interactions with other organisms. We anticipate that quantitative and scalable data on the growth rates of soil microorganisms, coupled with measurements of biogeochemical fluxes, will allow scientists to test and refine ecological theory and advance process-based models of carbon flux, nutrient uptake, and ecosystem productivity. Finally, measurements of in situ microbial growth rates provide insights into the ecology of populations and can be used to quantitatively link microbial diversity to soil biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The Q 10 of in situ microbial soil respiration varies with mean annual temperature, precipitation, pH, and plant cover: a meta-analysis and spatial prediction of Q 10

The temperature sensitivity of soil microbial respiration, commonly quantified using the Q 10 coefficient, is a key parameter in carbon cycle models. Uncovering how environmental factors affect in situ Q 10 values can therefore provide critical insight into potential shifts in global carbon stocks under climate change. We collected data from previously published field experiments that measured soil microbial respiration across a range of temperatures. We hypothesized that the Q 10 coefficient of in situ soil microbial respiration would vary based on environmental factors including mean annual temperature (MAT), mean annual precipitation (MAP), plant cover type, pH, soil C:N, and latitude. Linear regression revealed that Q 10 correlates negatively with MAT and MAP and positively with pH and absolute latitude. Additionally, average Q 10 varied significantly across different plant cover types; it was highest in mountain grasslands and lowest in tropical moist forests. Variation in microbial Q 10 across environmental factors may arise from underlying mechanisms such as enzyme kinetics, substrate availability and complexity, and microbial adaptation. To capture patterns in Q 10 more comprehensively, we developed a multiple linear regression model of Q 10 based on the most individually significant environmental drivers and applied it to public datasets to generate a global map of predicted Q 10 . Q 10 was higher in high-latitude and high-altitude regions, where large permafrost carbon stores are vulnerable to thawing and decomposition. We also compared fits between the Q 10 equation and a model produced from macromolecular rate theory (MMRT). We found that the MMRT model had the superior fit and may be better suited to model temperature sensitivity of complex biological reactions. Overall, our results emphasize that relationships between microbial Q 10 and environmental variables should be accounted for in climate models. Incorporating these variations in the Q 10 parameter, rather than using a fixed value, will help predict whether CO 2 emissions will be buffered or exacerbated by soil microbial respiration under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

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↗