Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
DOE OSTI2020

Mobility, Energy, and Electric Vehicle Typology for New York State

Mobility patterns, technology adoption and associated energy outcomes vary tremendously across settlement types. This paper explores how a highly geographically resolved exploration of the social, economic, techno-infrastructural and environmental domains of New York State is key to understanding observed variations in transportation technology adoption and associated mobility and energy outcomes. Current socioeconomic and mobility data sets at the census block group level are integrated in a hierarchical clustering approach to show how variations in mobility and energy outcomes are shaped by these domains - which may enhance progress on investments, or effectively inform planning and other decisions for state-wide goals. The clustering produces four settlement types to predict dependent variables of electric vehicle (EV) adoption rates, commute mode, vehicle fuel economy, and vehicles per household. This typology shows EV adoption rates among the core urban population, which is wealthy and highly-educated, are high - 3 EVs/1000 vehicles versus 1 EV/1000 vehicles among the urban working class. Commuting mode is closely linked with population and employment density - more than 90% of core urbanites use transit or active modes, compared with only 22% of suburbanites and 17% of rural residents. Household vehicle ownership also varies, with 2 vehicles per household in rural areas and only 0.5 in core urban settings. Important findings on differences among the rural, suburban, urban, and urban core settlement types suggest a need to explore how best to manage and anticipate very different types of services that may be supportive in achieving energy-efficient and affordable mobility systems state-wide.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Bootstrapping outperforms community‐weighted approaches for estimating the shapes of phenotypic distributions

Abstract Estimating phenotypic distributions of populations and communities is central to many questions in ecology and evolution. These distributions can be characterized by their moments (mean, variance, skewness and kurtosis) or diversity metrics (e.g. functional richness). Typically, such moments and metrics are calculated using community‐weighted approaches (e.g. abundance‐weighted mean). We propose an alternative bootstrapping approach that allows flexibility in trait sampling and explicit incorporation of intraspecific variation, and show that this approach significantly improves estimation while allowing us to quantify uncertainty. We assess the performance of different approaches for estimating the moments of trait distributions across various sampling scenarios, taxa and datasets by comparing estimates derived from simulated samples with the true values calculated from full datasets. Simulations differ in sampling intensity (individuals per species), sampling biases (abundance, size), trait data source (local vs. global) and estimation method (two types of community‐weighting, two types of bootstrapping). We introduce the traitstrap R package, which contains a modular and extensible set of bootstrapping and weighted‐averaging functions that use community composition and trait data to estimate the moments of community trait distributions with their uncertainty. Importantly, the first function in the workflow, trait_fill , allows the user to specify hierarchical structures (e.g. plot within site, experiment vs. control, species within genus) to assign trait values to each taxon in each community sample. Across all taxa, simulations and metrics, bootstrapping approaches were more accurate and less biased than community‐weighted approaches. With bootstrapping, a sample size of 9 or more measurements per species per trait generally included the true mean within the 95% CI. It reduced average percent errors by 26%–74% relative to community‐weighting. Random sampling across all species outperformed both size‐ and abundance‐biased sampling. Our results suggest randomly sampling ~9 individuals per sampling unit and species, covering all species in the community and analysing the data using nonparametric bootstrapping generally enable reliable inference on trait distributions, including the central moments, of communities. By providing better estimates of community trait distributions, bootstrapping approaches can improve our ability to link traits to both the processes that generate them and their effects on ecosystems.

Maitner, Brian S.↗
NASA NTRS1999

Global Geochemical Variation on the Lunar Surface: A Three-Element Approach

We present a method for displaying the relative abundances of three important elements (Th, Fe, and Ti) on the same map projection of the lunar surface. Using Th-, Fe-, and Ti-elemental abundances from orbital geochemical data and assigning each element a primary color, a false-color map of the lunar surface was created. This approach is similar to the ternary diagram approach presented by Davis and Spudis with some important differences, discussed later. For the present maps, Th abundances were measured by the Lunar Prospector (LP) Gamma-Ray Spectrometer(GRS).The new LPGRS low-altitude dataset was used in this analysis. Iron and Ti weight percentages were based on Clementine spectral reflectance data smoothed to the LP low altitude footprint. This method of presentation was designed to aid in the location and recognition of three principal lunar compositions: ferroan anorthosite (FAN), mare basalts (MB), and the Mg suite/ KREEP-rich rocks on the lunar surface, with special emphasis on the highlands and specific impact basins. In addition to the recognition of these endmember rock compositions, this method is an attempt to examine the relationship between elemental compositions that do not conform readily to previously accepted or observed endmember rocks in various specific regions of interest, including eastern highlands regions centered on 150 deg longitude, and a northern highlands Th-rich region observed. The LP low-altitude data has full width at half-maximum spatial resolution of about 40 km. The Clementine spectral reflectance datasets were adapted using an equal-area, gaussian smoothing routine to this footprint. In addition, these datasets, reported in weight percent of FeO and of Ti02, were adjusted to Fe and Ti weight percentages. Each dataset was then assigned one of the three primary colors: blue for Th, red for Fe, and green for Ti. For each element, the data range was normalized to represent the ratio of each point to the maximum in the dataset. (To view the color table, go to http://cass.jsc.nasa.gov/meetings/moon99/pdf/8033.pdf.) The full range of lunar longitudes is represented, but due to the lack of coverage of the Clementine data for latitudes > 70 deg and <-70 deg, the data for these regions is excluded. The differences between this approach and the ternary diagram approach of Davis and Spudis eliminate some of the uncertainty and ambiguity of the ternary diagram approach. Rather than using a ratio of Th to Ti normalized to CI chondritic ratios, and a ternary diagram with ternary apexes located at specific endmember compositional values, elemental compositions were used independently, eliminating the errors resulting from dividing numbers that can have high uncertainties, especially at low concentration. The three elements used in this method of presentation were chosen for several reasons. One reason for the inclusion of Th in this study is that it is an accurate indicator of KREEP. Iron and Ti concentrations are both low in highland regolith, causing any small fluctuations in Th to stand out very well. In addition, Fe and Ti are good compositional indicators of different mare basalts. Mixed with red for Fe, the green for Ti produces a yellow signal in high-Ti basalts. While universally high in Fe relative to the surrounding highlands, mare basalts have a diverse range of Ti values, making Ti concentration a valuable asset to the classification and identification of different basalt types. Finally, an important constraint in element selection is the availability of the global data, both from LP and Clementine results. Data for Th, Fe, and Ti are among the highest quality of existing lunar remote-sensing data. In addition, LP data for Fe and Ti will become available, enabling these data to be incorporated into the analysis. Using upper-limit values for end member rock compositions calculated from Korotev et al., attempts were made to locate the different endmember compositions of terranes on this diagram. Most strikingly, ferroan anorthosite (Th < and = 0.37 micro g/g; Fe (wt%)< and =2.29; Ti (wt%) < and = 0.22), which should appear as an almost black, reddish color, does not appear on the diagram at any noticeable frequency. Based on this analysis, the suggestion of extensive FAN regions on the lunar surface is not strong, especially at the presently accepted values for Fe and Th. However, to make sure this effect is not due to systematic errors, a thorough investigation of the precision, accuracy, and uncertainties of the Fe, Ti, and Th abundances needs to be carried out, especially at low concentrations. A particular region of interest is an area of high Th concentrations relative to Fe and Ti content north and east of Humboldtianum Crater. First observed by Lawrence et al., this region does not coincide with any visible impact structure and comprises one of the closest approximations to pure blue (high Th, very low Ti and Fe) on the lunar surface. Such an elemental composition does not lend itself readily to classification, and presents something of an anomaly. More detailed analysis of this region is needed to understand its structure and origin. There seems to be a longitudinal asymmetry in the Th concentrations of the highlands regolith. High-Th, low-Ti, and Fe regions are located between 135 deg and 180 deg longitude and between -30 deg and +30 deg latitude. While the Th levels are not high enough to attract attention in a single elemental display, the variation in the abundance of Th relative to Fe and Ti abundances can be clearly seen. The composition that these data suggest is not well represented in the sample return suite. In addition, these regions were largely missed by the Apollo orbital ground tracks, which only covered the outer edge of the areas of interest. The LP orbital Th data represent the first information about the Th concentrations in these regions of the highlands. Additional information contained in original.

Thomsen, D. R.↗
NASA NTRS2004

Variation in sulfide tolerance of photosystem II in phylogenetically diverse cyanobacteria from sulfidic habitats

Physiological and molecular phylogenetic approaches were used to investigate variation among 12 cyanobacterial strains in their tolerance of sulfide, an inhibitor of oxygenic photosynthesis. Cyanobacteria from sulfidic habitats were found to be phylogenetically diverse and exhibited an approximately 50-fold variation in photosystem II performance in the presence of sulfide. Whereas the degree of tolerance was positively correlated with sulfide levels in the environment, a strain's phenotype could not be predicted from the tolerance of its closest relatives. These observations suggest that sulfide tolerance is a dynamic trait primarily shaped by environmental variation. Despite differences in absolute tolerance, similarities among strains in the effects of sulfide on chlorophyll fluorescence induction indicated a common mode of toxicity. Based on similarities with treatments known to disrupt the oxygen-evolving complex, it was concluded that sulfide toxicity resulted from inhibition of the donor side of photosystem II.

Phylogeny↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Limited-angle computed tomography with deep image and physics priors

Computed tomography is a well-established x-ray imaging technique to reconstruct the three-dimensional structure of objects. It has been used extensively in a variety of fields, from diagnostic imaging to materials and biological sciences. One major challenge in some applications, such as in electron or x-ray tomography systems, is that the projections cannot be gathered over all the angles due to the sample holder setup or shape of the sample. This results in an ill-posed problem called the limited angle reconstruction problem. Typical image reconstruction in this setup leads to distortion and artifacts, thereby hindering a quantitative evaluation of the results. To address this challenge, we use a generative model to effectively constrain the solution of a physics-based approach. Our approach is self-training that can iteratively learn the nonlinear mapping from partial projections to the scanned object. Because our approach combines the data likelihood and image prior terms into a single deep network, it is computationally tractable and improves performance through an end-to-end training. We also complement our approach with total-variation regularization to handle high-frequency noise in reconstructions and implement a solver based on alternating direction method of multipliers. We present numerical results for various degrees of missing angle range and noise levels, which demonstrate the effectiveness of the proposed approach.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Toward Consistent High-Fidelity Quantum Learning on Unstable Devices via Efficient In-Situ Calibration

In the near-term noisy intermediate-scale quantum (NISQ) era, high noise will significantly reduce the fidelity of quantum computing. What's worse, recent works reveal that the noise on quantum devices is not stable, that is, the noise is dynamically changing over time. This leads to an imminent challenging problem: At run-time, is there a way to efficiently achieve a consistent high-fidelity quantum system on unstable devices? To study this problem, we take quantum learning (a.k.a., variational quantum algorithm) as a vehicle, which has a wide range of applications, such as combinatorial optimization and machine learning. A straightforward approach is to optimize a variational quantum circuit (VQC) with a parameter-shift approach on the target quantum device before using it; however, the optimization has an extremely high time cost, which is not practical at run-time. To address the pressing issue, in this paper, we proposed a novel quantum pulse-based noise adaptation framework, namely QuPAD. In the proposed framework, first, we identify that the CNOT gate is the fidelity bottleneck of the conventional VQC, and we employ a more robust parameterized multi-qubit gate (i.e., Rzx gate) to replace CNOT gate. Second, by benchmarking Rzx gate with different parameters, we build a fitting function for each coupling qubit pair, such that the deviation between the theoretic output of Rzx gate and its on-device output under a given pulse amplitude and duration can be efficiently predicted. On top of this, an evolutionary algorithm is devised to identify the pulse amplitude and duration of each Rzx gate (i.e., calibration) and find the quantum circuits with high fidelity. Experiments show that the runtime on quantum devices of QuPAD with 8–10 qubits is less than 15 minutes, which is up to 270 x faster than the parameter-shift approach. In addition, compared to the vanilla VQC as a baseline, QuPAD can achieve 59.33% accuracy gain on a classification task, and average 66.34% closer to ground state energy for molecular simulation.

Hu, Zhirui↗
NASA NTRS2022

Along-Orbit Analysis of GRACE Follow-On Inter-Satellite Laser Ranging Measurements for Sub-Monthly Surface Mass Variations

We examined the sensitivity of GRACE Follow-On (GRACE-FO) laser ranging interferometer (LRI) measurements to sub-monthly time-variable gravity (TVG) signals caused by transient, high-frequency mass changes in the Earth system. GRACE-FO LRI provides complementary inter-satellite ranging measurements with higher precision over a wider range of frequencies than the baseline K-band microwave ranging system. The common approach for studying mass variation relies on the inverted TVG or mascon solutions over a period of, for example, one month or 10 days which are adversely affected by temporal aliasing and/or smoothing. In this article, we present the alternative along-orbit analysis methodology in terms of line-of-sight gravity difference (LGD) to fully exploit the higher precision LRI measurements for examination of sub-monthly mass changes. The discrepancy between “instantaneous” LGD LRI observations and monthly mean LGD (from Level-2 data) at satellite altitude indicates the sub-monthly gravitational variability not captured by monthly-mean solutions. In conjunction with the satellite ocean altimetry observations, high frequency non-tidal atmosphere and ocean models, and hydrology models, we show that the LGD LRI observations detect the high-frequency oceanic mass variability in the Argentine Basin and the Gulf of Carpentaria, and sub-monthly variations in surface (river) water in the Amazon Basin. We demonstrate the benefits gained from repeat ground track analysis of GRACE-FO LRI data in the case of the Amazon surface water flow. The along-orbit analysis methodology based on LGD LRI time series presented here is especially suitable for quantifying temporal and spatial evolution of extreme, rapidly changing mass variations.

Khosro Ghobadi-Far↗
DOE OSTI2023

Short-term variation in leaf-level water use efficiency in a tropical forest

Here we report the representation of stomatal regulation of transpiration and CO 2 assimilation is key to forecasting terrestrial ecosystem responses to global change. Given its importance in determining the relationship between forest productivity and climate, accurate and mechanistic model representation of the relationship between stomatal conductance (g s ) and assimilation is crucial. We assess possible physiological and mechanistic controls on the estimation of the g 1 (stomatal slope, inversely proportional to water use efficiency) and g 0 (stomatal intercept) parameters, using diurnal gas exchange surveys and leaf level response curves of six tropical broadleaf evergreen tree species. g 1 estimated from ex-situ response curves averaged 50% less than g 1 estimated from survey data. While g 0 and g 1 varied between leaves of different phenological stages, the trend was not consistent among species. We identified a diurnal trend associated with g 1 and g 0 that significantly improved model projections of diurnal trends in transpiration. The accuracy of modelled g s can be improved by accounting for variation in stomatal behavior across diurnal periods, and between measurement approaches, rather than focusing on phenological variation in stomatal behavior. Additional investigation into the primary mechanisms responsible for diurnal variation in g 1 will be required to account for this phenomena in land surface models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Shifts in bacterial traits under chronic nitrogen deposition align with soil processes in arbuscular, but not ectomycorrhizal-associated trees

Nitrogen (N) deposition increases soil carbon (C) storage by reducing microbial activity. These effects vary in soil beneath trees that associate with arbuscular (AM) and ectomycorrhizal (ECM) fungi. Variation in carbon C and N uptake traits among microbes may explain differences in soil nutrient cycling between mycorrhizal associations in response to high N loads, a mechanism not previously examined due to methodological limitations. Here, we used quantitative Stable Isotope Probing (qSIP) to measure bacterial C and N assimilation rates from an added organic compound, which we conceptualize as functional traits. As such, we applied a trait-based approach to explore whether variation in assimilation rates of bacterial taxa can inform shifts in soil function under chronic N deposition. We show taxon-specific and community-wide declines of bacterial C and N uptake under chronic N deposition in both AM and ECM soils. N deposition-induced reductions in microbial activity were mirrored by declines in soil organic matter mineralization rates in AM but not ECM soils. Our findings suggest C and N uptake traits of bacterial communities can predict C cycling feedbacks to N deposition in AM soils, but additional data, for instance on the traits of fungi, may be needed to connect microbial traits with soil C and N cycling in ECM systems. In conclusion, our study also highlights the potential of employing qSIP in conjunction with trait-based analytical approaches to inform how ecological processes of microbial communities influence soil functioning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Automated label flows for excited states of correlation functions in lattice gauge theory

Extracting excited states from lattice gauge theory correlation functions can be achieved through $χ^2$ minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. To correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. This is a black-box approach in the sense that it gives a sensible set of labels without user guidance. Further, as an example, we pair one black box method with another, analyzing a lattice correlation function from real data using automated label flows in the context of Prony's method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

Multi-Variate LSTM Prediction of Alaska Magnetometer Chain Utilizing a Coupled Model Approach

During periods of rapidly changing geomagnetic conditions electric fields form within the Earth’s surface and induce currents known as geomagnetically induced currents(GICs), which interact with unprotected electrical systems our society relies on. In this study, we train multi-variate Long-Short Term Memory neural networks to predict magnitude of north-south component of the geomagnetic field (|BN|) at multiple ground magnetometer stations across Alaska provided by the SuperMAG database with a future goal of predicting geomagnetic field disturbances. Each neural network is driven by solar wind and interplanetary magnetic field inputs from the NASA OMNI database spanning from 2000–2015 and is fine tuned for each station to maximize the effectiveness in predicting |BN|. The neural networks are then compared against multivariate linear regression models driven with the same inputs at each station using Heidke skill scores with thresholds at the 50, 75, 85, and 99 percentiles for |BN|. The neural network models show significant increases over the linear regression models for |BN| thresholds. We also calculate the Heidke skill scores for d|BN|/dt by deriving d|BN|/dt from |BN| predictions. However, neural network models do not show clear outperformance compared to the linear regression models. To retain the sign information and thus predict BN instead of |BN|, a secondary so-called polarity model is utilized. The polarity model is run in tandem with the neural networks predicting geomagnetic field in a coupled model approach and results in a high correlation between predicted and observed values for all stations. We find this model a promising starting point for a machine learned geomagnetic field model to be expanded upon through increased output time history and fast turnaround times.

Matthew Blandin↗
DOE OSTI2025

Characterizing Seasonal Variation of the Atmospheric Mixing Layer Height Using Machine Learning Approaches

As machine learning becomes more integrated into atmospheric science, XGBoost has gained popularity for its ability to assess the relative contributions of influencing factors in the atmospheric boundary layer height. To examine how these factors vary across seasons, a seasonal analysis is necessary. However, dividing data by season reduces the sample size, which can affect result reliability and complicate factor comparisons. To address these challenges, this study replaces default parameters with grid search optimization and incorporates cross-validation to mitigate dataset limitations. Using XGBoost with four years of data from the atmospheric radiation measurement (ARM) (Southern Great Plains (SGP) C1 site, cross-validation stabilizes correlation coefficient fluctuations from 0.3 to within 0.1. With optimized parameters, the R value can reach 0.81. Analysis of the C1 site reveals that the relative importance of different factors changes across seasons. Lower tropospheric stability (LTS, ~0.53) is the dominant factor at C1 throughout the year. However, during DJF, latent heat flux (LHF, 0.44) surpasses LTS (0.22). In SON, LTS (0.58) becomes more influential than LHF (0.18). Further comparisons among the four long-term SGP sites (C1, E32, E37, and E39) show seasonal variations in relative importance. Notably, during JJA, the differences in the relative importance of the three factors across all sites are lower than in other seasons. This suggests that boundary layer development in the summer is not dominated by a single factor, reflecting a more intricate process likely influenced by seasonal conditions such as enhanced convective activity, higher temperatures, and humidity, which collectively contribute to a balanced distribution of parameter impacts. Furthermore, the relative importance of LTS gradually increases from morning to noon, indicating that LTS becomes more significant as the boundary layer approaches its maximum height. Consequently, the LTS in the early morning in autumn exhibits greater relative importance compared to other seasons. This reflects a faster development of the mixing layer height (MLH) in autumn, suggesting that it is easier to retrieve the MLH from the previous day during this period. The findings enhance understanding of boundary layer evolution and contribute to improved boundary layer parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Single Episode Policy Transfer in Reinforcement Learning

Transfer and adaptation to new unknown environmental dynamics is a key challenge for reinforcement learning (RL). An even greater challenge is performing near-optimally in a single attempt at test time, possibly without access to dense rewards, which is not addressed by current methods that require multiple experience rollouts for adaptation. To achieve single episode transfer in a family of environments with related dynamics, we propose a general algorithm that optimizes a probe and an inference model to rapidly estimate underlying latent variables of test dynamics, which are then immediately used as input to a universal control policy. This modular approach enables integration of state-of-the-art algorithms for variational inference or RL. Moreover, our approach does not require access to rewards at test time, allowing it to perform in settings where existing adaptive approaches cannot. In diverse experimental domains with a single episode test constraint, we report our method significantly outperforms existing adaptive approaches and shows favorable performance against baselines for robust transfer.

42 ENGINEERING↗