Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Determining the tilt of the Raman laser beam using an optical method for atom gravimeters

The tilt of a Raman laser beam is a major systematic error in precision gravity measurement using atom interferometry. The conventional approach to evaluating this tilt error involves modulating the direction of the Raman laser beam and conducting time-consuming gravity measurements to identify the error minimum. In this work, we demonstrate a method to expediently determine the tilt of the Raman laser beam by transforming the tilt angle measurement into characterization of parallelism, which integrates the optical method of aligning the laser direction, commonly used in freely falling corner-cube gravimeters, into an atom gravimeter. A position-sensing detector (PSD) is utilized to quantitatively characterize the parallelism between the test beam and the reference beam, thus measuring the tilt precisely and rapidly. After carefully positioning the PSD and calibrating the relationship between the distance measured by the PSD and the tilt angle measured by the tiltmeter, we achieved a statistical uncertainty of less than 30 µrad in the tilt measurement. Furthermore, we compared the results obtained through this optical method with those from the conventional tilt modulation method for gravity measurement. The comparison validates that our optical method can achieve tilt determination with an accuracy level of better than 200 µrad, corresponding to a systematic error of 20 µGal in g measurement. This work has practical implications for real-world applications of atom gravimeters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Accuracy, Bias, and Improvements in Mapping Crops and Cropland across the United States Using the USDA Cropland Data Layer

The U.S. Department of Agriculture’s (USDA) Cropland Data Layer (CDL) is a 30 m resolution crop-specific land cover map produced annually to assess crops and cropland area across the conterminous United States. Despite its prominent use and value for monitoring agricultural land use/land cover (LULC), there remains substantial uncertainty surrounding the CDLs’ performance, particularly in applications measuring LULC at national scales, within aggregated classes, or changes across years. To fill this gap, we used state- and land cover class-specific accuracy statistics from the USDA from 2008 to 2016 to comprehensively characterize the performance of the CDL across space and time. We estimated nationwide area-weighted accuracies for the CDL for specific crops as well as for the aggregated classes of cropland and non-cropland. We also derived and reported new metrics of superclass accuracy and within-domain error rates, which help to quantify and differentiate the efficacy of mapping aggregated land use classes (e.g., cropland) among constituent subclasses (i.e., specific crops). We show that aggregate classes embody drastically higher accuracies, such that the CDL correctly identifies cropland from the user’s perspective 97% of the time or greater for all years since nationwide coverage began in 2008. We also quantified the mapping biases of specific crops throughout time and used these data to generate independent bias-adjusted crop area estimates, which may complement other USDA survey- and census-based crop statistics. Our overall findings demonstrate that the CDLs provide highly accurate annual measures of crops and cropland areas, and when used appropriately, are an indispensable tool for monitoring changes to agricultural landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗
DOE OSTI2017

Geometry and growth contributions to cosmic shear observables

We explore the sensitivity of weak-lensing observables to the expansion history of the Universe and to the growth of cosmic structures, as well as the relative contribution of both effects to constraining cosmological parameters. We utilize ray-tracing dark-matter-only N-body simulations and validate our technique by comparing our results for the convergence power spectrum with analytic results from past studies. We then extend our analysis to non-Gaussian observables which cannot be easily treated analytically. We study the convergence (equilateral) bispectrum and two topological observables, lensing peaks and Minkowski functionals, focusing on their sensitivity to the matter density Ω m and the dark energy equation of state w. We find that a cancellation between the geometry and growth effects is a common feature for all observables and exists at the map level. It weakens the overall sensitivity by factors of up to 3 and 1.5 for w and Ω m , respectively, with the bispectrum worst affected. However, combining geometry and growth information alleviates the degeneracy between Ω m and w from either effect alone. As a result, the magnitudes of marginalized errors remain similar to those obtained from growth-only effects, but with the correlation between the two parameters switching sign. Furthermore, these results shed light on the origin of the cosmology sensitivity of non-Gaussian statistics and should be useful in optimizing combinations of observables.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Framework for the Evolution of Heuristics in Advanced Manufacturing

Here, in this study, we work toward addressing a knowledge gap in understanding how heuristics are developed, retrieved, employed, and modified by designers. Having a better awareness of one’s own set of heuristics can be beneficial for relaying to other team members, improving a team’s training processes, and aiding others on their path to design expertise. The ability to understand and justify the use of a heuristic should lead to more effective decision-making in systems design. To do this, the heuristics and their characteristics must be extracted using a repeatable scientific research methodology. This study describes a unique extraction and characterization process compared to prior literature. It includes some of the first work towards documenting heuristics for both designers and operators in a hybrid manufacturing setting. Eight participants performed a series of two design journals, two interviews, and one survey. Heuristics were extracted and refined between each method and then verified by participants in the survey. The surveys produced novel statistically significant findings in regard to heuristic characterizations, impacting how participants view how often a heuristic is used, the reliability of the heuristic, and the evolution of the heuristic. Lastly, an alternate perspective of heuristics as an error management bias is highlighted and discussed.

42 ENGINEERING↗
DOE OSTI2012

In aqua vivo EPID dosimetry

Purpose: At the Netherlands Cancer Institute-Antoni van Leeuwenhoek Hospital in vivo dosimetry using an electronic portal imaging device (EPID) has been implemented for almost all high-energy photon treatments of cancer with curative intent. Lung cancer treatments were initially excluded, because the original back-projection dose-reconstruction algorithm uses water-based scatter-correction kernels and therefore does not account for tissue inhomogeneities accurately. The aim of this study was to test a new method, in aqua vivo EPID dosimetry, for fast dose verification of lung cancer irradiations during actual patient treatment. Methods: The key feature of our method is the dose reconstruction in the patient from EPID images, obtained during the actual treatment, whereby the images have been converted to a situation as if the patient consisted entirely of water; hence, the method is termed in aqua vivo. This is done by multiplying the measured in vivo EPID image with the ratio of two digitally reconstructed transmission images for the unit-density and inhomogeneous tissue situation. For dose verification, a comparison is made with the calculated dose distribution with the inhomogeneity correction switched off. IMRT treatment verification is performed for each beam in 2D using a 2D {gamma} evaluation, while for the verification of volumetric-modulated arc therapy (VMAT) treatments in 3D a 3D {gamma} evaluation is applied using the same parameters (3%, 3 mm). The method was tested using two inhomogeneous phantoms simulating a tumor in lung and measuring its sensitivity for patient positioning errors. Subsequently five IMRT and five VMAT clinical lung cancer treatments were investigated, using both the conventional back-projection algorithm and the in aqua vivo method. The verification results of the in aqua vivo method were statistically analyzed for 751 lung cancer patients treated with IMRT and 50 lung cancer patients treated with VMAT. Results: The improvements by applying the in aqua vivo approach are considerable. The percentage of {gamma} values {<=}1 increased on average from 66.2% to 93.1% and from 43.6% to 97.5% for the IMRT and VMAT cases, respectively. The corresponding mean {gamma} value decreased from 0.99 to 0.43 for the IMRT cases and from 1.71 to 0.40 for the VMAT cases, which is similar to the accepted clinical values for the verification of IMRT treatments of prostate, rectum, and head-and-neck cancers. The deviation between the reconstructed and planned dose at the isocenter diminished on average from 5.3% to 0.5% for the VMAT patients and was almost the same, within 1%, for the IMRT cases. The in aqua vivo verification results for IMRT and VMAT treatments of a large group of patients had a mean {gamma} of approximately 0.5, a percentage of {gamma} values {<=}1 larger than 89%, and a difference of the isocenter dose value less than 1%. Conclusions: With the in aqua vivo approach for the verification of lung cancer treatments (IMRT and VMAT), we can achieve results with the same accuracy as obtained during in vivo EPID dosimetry of sites without large inhomogeneities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Mammographic quantitative image analysis and biologic image composition for breast lesion characterization and classification

Purpose: To investigate whether biologic image composition of mammographic lesions can improve upon existing mammographic quantitative image analysis (QIA) in estimating the probability of malignancy. Methods: The study population consisted of 45 breast lesions imaged with dual-energy mammography prior to breast biopsy with final diagnosis resulting in 10 invasive ductal carcinomas, 5 ductal carcinomain situ, 11 fibroadenomas, and 19 other benign diagnoses. Analysis was threefold: (1) The raw low-energy mammographic images were analyzed with an established in-house QIA method, “QIA alone,” (2) the three-compartment breast (3CB) composition measure—derived from the dual-energy mammography—of water, lipid, and protein thickness were assessed, “3CB alone”, and (3) information from QIA and 3CB was combined, “QIA + 3CB.” Analysis was initiated from radiologist-indicated lesion centers and was otherwise fully automated. Steps of the QIA and 3CB methods were lesion segmentation, characterization, and subsequent classification for malignancy in leave-one-case-out cross-validation. Performance assessment included box plots, Bland–Altman plots, and Receiver Operating Characteristic (ROC) analysis. Results: The area under the ROC curve (AUC) for distinguishing between benign and malignant lesions (invasive and DCIS) was 0.81 (standard error 0.07) for the “QIA alone” method, 0.72 (0.07) for “3CB alone” method, and 0.86 (0.04) for “QIA+3CB” combined. The difference in AUC was 0.043 between “QIA + 3CB” and “QIA alone” but failed to reach statistical significance (95% confidence interval [–0.17 to + 0.26]). Conclusions: In this pilot study analyzing the new 3CB imaging modality, knowledge of the composition of breast lesions and their periphery appeared additive in combination with existing mammographic QIA methods for the distinction between different benign and malignant lesion types.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

TH-E-17A-10: Markerless Lung Tumor Tracking Based On Beams Eye View EPID Images

Purpose: Dynamic tumor tracking or motion compensation techniques have proposed to modify beam delivery following lung tumor motion on the flight. Conventional treatment plan QA could be performed in advance since every delivery may be different. Markerless lung tumor tracking using beams eye view EPID images provides a best treatment evaluation mechanism. The purpose of this study is to improve the accuracy of the online markerless lung tumor motion tracking method. Methods: The lung tumor could be located on every frame of MV images during radiation therapy treatment by comparing with corresponding digitally reconstructed radiograph (DRR). A kV-MV CT corresponding curve is applied on planning kV CT to generate MV CT images for patients in order to enhance the similarity between DRRs and MV treatment images. This kV-MV CT corresponding curve was obtained by scanning a same CT electron density phantom by a kV CT scanner and MV scanner (Tomotherapy) or MV CBCT. Two sets of MV DRRs were then generated for tumor and anatomy without tumor as the references to tracking the tumor on beams eye view EPID images. Results: Phantom studies were performed on a Varian TrueBeam linac. MV treatment images were acquired continuously during each treatment beam delivery at 12 gantry angles by iTools. Markerless tumor tracking was applied with DRRs generated from simulated MVCT. Tumors were tracked on every frame of images and compared with expected positions based on programed phantom motion. It was found that the average tracking error were 2.3 mm. Conclusion: This algorithm is capable of detecting lung tumors at complicated environment without implanting markers. It should be noted that the CT data has a slice thickness of 3 mm. This shows the statistical accuracy is better than the spatial accuracy. This project has been supported by a Varian Research Grant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Assessing North Atlantic Tropical Cyclone Rainfall Hazard Using Engineered-Synthetic Storms and a Physics-Based Tropical Cyclone Rainfall Model

In this study, we design a statistical method to couple observations with a physics-based tropical cyclone (TC) rainfall model (TCR) and engineered-synthetic storms for assessing TC rainfall hazard. We first propose a bias-correction method to minimize the errors induced by TCR via matching the probability distribution of TCR-simulated historical TC rainfall with gauge observations. Then we assign occurrence probabilities to engineered-synthetic storms to reflect local climatology, through a resampling method that matches the probability distribution of a newly proposed storm parameter named rainfall potential (POT) in the synthetic dataset with that in the observation. POT is constructed to include several important storm parameters for TC rainfall such as TC intensity, duration, and distance and environmental humidity near landfall, and it is shown to be correlated with TCR-simulated rainfall. The proposed method has a satisfactory performance in reproducing the rainfall hazard curve in various locations in the continental United States; it is an improvement over the traditional joint probability method (JPM) for TC rainfall hazard assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Chandra Detection of Intracluster X-Ray sources in Virgo

We present a survey of X-ray point sources in the nearest and dynamically young galaxy cluster, Virgo, using archival Chandra observations that sample the vicinity of 80 early-type member galaxies. The X-ray source populations at the outskirts of these galaxies are of particular interest. We detect a total of 1046 point sources (excluding galactic nuclei) out to a projected galactocentric radius of ∼40 kpc and down to a limiting 0.5–8 keV luminosity of ∼2×10{sup 38} erg s{sup −1}. Based on the cumulative spatial and flux distributions of these sources, we statistically identify ∼120 excess sources that are not associated with the main stellar content of the individual galaxies, nor with the cosmic X-ray background. This excess is significant at a 3.5σ level, when Poisson error and cosmic variance are taken into account. On the other hand, no significant excess sources are found at the outskirts of a control sample of field galaxies, suggesting that at least some fraction of the excess sources around the Virgo galaxies are truly intracluster X-ray sources. Assisted with ground-based and HST optical imaging of Virgo, we discuss the origins of these intracluster X-ray sources, in terms of supernova-kicked low-mass X-ray binaries (LMXBs), globular clusters, LMXBs associated with the diffuse intracluster light, stripped nucleated dwarf galaxies and free-floating massive black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Determination of the critical manifold tangent space and curvature with Monte Carlo renormalization group

We show that the critical manifold of a statistical mechanical system in the vicinity of a critical point is locally accessible through correlation functions at that point. A practical numerical method is presented to determine the tangent space and the curvature to the critical manifold with variational Monte Carlo renormalization group. Because of the use of a variational bias potential of the coarse-grained variables, critical slowing down is greatly alleviated in the Monte Carlo simulation. In addition, this method is free of truncation error. Furthermore, we study the isotropic Ising model on square and cubic lattices, the anisotropic Ising model, and the tricritical Ising model on square lattices to illustrate the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Advanced Signal Decomposition Analysis and Anomaly Detection in Photovoltaic Systems

With the rapid expansion of large-scale photovoltaic (PV) plants, it is paramount for solar stakeholders to understand the reliability and efficiency of their plants to inform maintenance decisions, increase production, and understand the design factors that impact performance. Diagnosing underperformance in PV plants is challenging due to the relatively few monitoring points with respect to the large geographic footprint of the plant. This work introduces a cutting-edge method that transforms the analysis and management of key factors influencing PV plant performance, including performance loss rate (PLR), recoverable soiling, and major system changes. Identifying these factors is critical for deriving actionable insights. Leveraging advanced analytical techniques such as wavelet transformation, robust regression, and extreme point analysis, this approach provides a nuanced understanding of these factors. This method has been tested across two synthetic datasets and one real dataset, consistently surpassing existing benchmarks by achieving a lower median mean absolute error and reduced error variability across all comparable components.

14 SOLAR ENERGY↗
DOE OSTI2015

Statistical precision of the intensities retrieved from constrained fitting of overlapping peaks in high-resolution mass spectra

Least-squares fitting of overlapping peaks is often needed to separately quantify ions in high-resolution mass spectrometer data. A statistical simulation approach is used to assess the statistical precision of the retrieved peak intensities. The sensitivity of the fitted peak intensities to statistical noise due to ion counting is probed for synthetic data systems consisting of two overlapping ion peaks whose positions are pre-defined and fixed in the fitting procedure. The fitted intensities are sensitive to imperfections in the m/Q calibration. These propagate as a limiting precision in the fitted intensities that may greatly exceed the precision arising from counting statistics. The precision on the fitted peak intensity falls into one of three regimes. In the "counting-limited regime" (regime I), above a peak separation &chi; ~ 2 to 3 half-widths at half-maximum (HWHM), the intensity precision is similar to that due to counting error for an isolated ion. For smaller &chi; and higher ion counts (~ 1000 and higher), the intensity precision rapidly degrades as the peak separation is reduced ("calibration-limited regime", regime II). Alternatively for &chi; < 1.6 but lower ion counts (e.g. 10–100) the intensity precision is dominated by the additional ion count noise from the overlapping ion and is not affected by the imprecision in the m/Q calibration ("overlapping-limited regime", regime III). The transition between the counting and m/Q calibration-limited regimes is shown to be weakly dependent on resolving power and data spacing and can thus be approximated by a simple parameterisation based only on peak intensity ratios and separation. A simple equation can be used to find potentially problematic ion pairs when evaluating results from fitted spectra containing many ions. Longer integration times can improve the precision in regimes I and III, but a given ion pair can only be moved out of regime II through increased spectrometer resolving power. As a result, studies presenting data obtained from least-squares fitting procedures applied to mass spectral peaks should explicitly consider these limits on statistical precision.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

A controlled statistical study to assess measurement variability as a function of test object position and configuration for automated surveillance in a multicenter longitudinal COPD study (SPIROMICS)

Purpose: A test object (phantom) is an important tool to evaluate comparability and stability of CT scanners used in multicenter and longitudinal studies. However, there are many sources of error that can interfere with the test object-derived quantitative measurements. Here the authors investigated three major possible sources of operator error in the use of a test object employed to assess pulmonary density-related as well as airway-related metrics. Methods: Two kinds of experiments were carried out to assess measurement variability caused by imperfect scanning status. The first one consisted of three experiments. A COPDGene test object was scanned using a dual source multidetector computed tomographic scanner (Siemens Somatom Flash) with the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) inspiration protocol (120 kV, 110 mAs, pitch = 1, slice thickness = 0.75 mm, slice spacing = 0.5 mm) to evaluate the effects of tilt angle, water bottle offset, and air bubble size. After analysis of these results, a guideline was reached in order to achieve more reliable results for this test object. Next the authors applied the above findings to 2272 test object scans collected over 4 years as part of the SPIROMICS study. The authors compared changes of the data consistency before and after excluding the scans that failed to pass the guideline. Results: This study established the following limits for the test object: tilt index ≤0.3, water bottle offset limits of [−6.6 mm, 7.4 mm], and no air bubble within the water bottle, where tilt index is a measure incorporating two tilt angles around x- and y-axis. With 95% confidence, the density measurement variation for all five interested materials in the test object (acrylic, water, lung, inside air, and outside air) resulting from all three error sources can be limited to ±0.9 HU (summed in quadrature), when all the requirements are satisfied. The authors applied these criteria to 2272 SPIROMICS scans and demonstrated a significant reduction in measurement variation associated with the test object. Conclusions: Three operator errors were identified which significantly affected the usability of the acquired scan images of the test object used for monitoring scanner stability in a multicenter study. The authors’ results demonstrated that at the time of test object scan receipt at a radiology core laboratory, quality control procedures should include an assessment of tilt index, water bottle offset, and air bubble size within the water bottle. Application of this methodology to 2272 SPIROMICS scans indicated that their findings were not limited to the scanner make and model used for the initial test but was generalizable to both Siemens and GE scanners which comprise the scanner types used within the SPIROMICS study.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

An exploration of heterogeneity in supernova type Ia samples

We examine three SNe Type Ia datasets: Union2.1, JLA and Panstarrs to check their consistency using cosmology blind statistical analyses as well as cosmological parameter fitting. We find that the Panstarrs dataset is the most stable of the three to changes in the data, although it does not, at the moment, go to high enough redshifts to tightly constrain the equation of state of dark energy, w . The Union2.1, drawn from several different sources, appears to be somewhat susceptible to changes within the dataset. The JLA reconstructs well for a smaller number of cosmological parameters. At higher degrees of freedom, the dependence of its errors on redshift can lead to varying results between subsets. Panstarrs is inconsistent with the other two datasets at about 2σ confidence level, and JLA and Union2.1 are about 1σ away from each other. For the Ω{sub 0} {sub m} − w cosmological reconstruction, with no additional data, the 1σ range of values in w for selected subsets of each dataset is two times larger for JLA and Union2.1 as compared to Panstarrs. The range in Ω{sub 0} {sub m} for the same subsets remains approximately similar for all three datasets. We find that although there are differences in the fitting and correction techniques used in the different samples, the most important criterion is the selection of the SNe, a slightly different SNe selection can lead to noticeably different results both in the purely statistical analysis and in cosmological reconstruction. We note that a single, high quality low redshift sample could help decrease the uncertainties in the result. We also note that lack of homogeneity in the magnitude errors may bias the results and should either be modeled, or its effect neutralized by using other, complementary datasets. A supernova sample with high quality data at both high and low redshifts, constructed from a few surveys to avoid heterogeneity in the sample, and with homogeneous errors, would result in a more robust cosmological reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

On the surface energy balance closure at different temporal scales

Measurements of the surface energy fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made in the Columbia River Basin (Oregon) in an area of complex terrain during a 10-month long portion of the second Wind Forecast Improvement Project (WFIP 2) field campaign are used to study the surface energy budget (SEB) and surface fluxes over different temporal scales. This study analyzes and discusses SEB closure based on half-hourly, daily, monthly, seasonal, and sub-annual (~10-month) temporal averages. The data were collected over all four seasons for different states of the underlying ground surface (dry, wet, and frozen). Our half-hourly direct measurements of energy balance show that the sum of the turbulent sensible and latent heat fluxes systematically underestimate positive net radiation by around 20–30% during daytime and overestimate negative net radiation at night. This imbalance of the surface energy budget is comparable to other terrestrial sites. However, on average, the residual energy imbalance is significantly reduced at daily, weekly, and monthly averaging timescales, and moreover, the SEB can be closed for this site within reasonable limits on seasonal and sub-annual timescales (311-day averaging for the entire field campaign dataset). Increasing the averaging time to daily and longer time intervals substantially reduces the ground heat flux and storage terms, because energy locally entering the soil, air column, and vegetation in the morning is released in the afternoon and evening. Averaging on daily to sub-annual timescales also reduces random instrumental measurement errors and other uncertainties as well as smooths out a hysteresis effect (phase lag) in the SEB relationship between different components. This study shows that SEB closure is better for dry soils compared to wet soils and the statistical dependence of the turbulent fluxes and net radiation for freezing soil surfaces appears weak, if not non-existent, apparently due to lack of the latent heat of fusion term in the traditional SEB equation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cluster-Graph Fingerprinting: A Framework for Quantitative Analysis of Machine-Learned Interatomic Model Training and Simulation Data

Machine-learned interatomic models represent a significant advancement in simulation methods, extending the predictive ability of first-principles methods to previously inaccessible length and time scales. However, the data-driven nature of these models can lead to difficult-to-detect errors that can compromise prediction accuracy. To address this challenge, we introduce a novel fingerprinting approach based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES) ML-IAM graph-based descriptor. Our strategy enables efficient and statistically rigorous analysis of system configurations used in ML-IAM training and those generated by their application, e.g., in molecular dynamics simulations. We demonstrate that these fingerprints can effectively assess novelty of a configuration relative to an existing data set and determine dissimilarity among individual configurations, which are two key tasks in workflows for active learning-based ML-IAM training, data set curation, and on-the-fly uncertainty quantification.

36 MATERIALS SCIENCE↗