Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampled value”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Determination of the Raman polarizability tensor in the optically anisotropic crystal potassium dihydrogen phosphate and its deuterated analog

Abstract The Raman tensor of the dominant A 1 modes of the nonlinear optical crystalline material potassium dihydrogen phosphate and its 70% deuterated analog have been ascertained. Challenges in determining the A 1 mode tensor element values based on previous reports have been resolved using a specially designed experimental setup that makes use of spherical crystal samples. This novel experimental design enabled the determination of measurement artifacts, including polarization rotation of the pump and/or scattered light propagating through the sample and the contribution of additional overlapping phonon modes, which have hindered previous efforts. Results confirmed that the polarization tensor is diagonal, and matrix elements were determined with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fracture toughness of schist, amphibolite, and rhyolite from the Sanford Underground Research Facility (SURF), Lead, South Dakota

The Cracked Chevron Notched Brazilian Disc (CCNBD) method was selected for Mode I fracture toughness tests on Poorman schist, Yates amphibolite, and rhyolite dikes from the EGS Collab site at the SURF in Lead, South Dakota. The effects of lithology, anisotropy, and loading rate were investigated. Fracture toughness was greatest in amphibolite, with schist and rhyolite having similar toughness values ($K$ amphibolite > $K$ rhyolite ≈ $K$ schist ). The effects of anisotropy on fracture toughness were investigated in the foliated schist samples. Schist samples were prepared in three geometries (divider, arrester, and short transverse) which controlled how the fracture would propagate relative to foliations. The divider geometry was strongest and short transverse geometry was the weakest ($K$ divider > $K$ arrester > $K$ shorttransverse ). Fracture toughness was observed to decrease with decreasing loading rate. Optical and SEM microscopy revealed that for the short transverse geometry, fractures tended to propagate along grain boundaries, whereas in arrester and divider geometries fractures tended to propagate through grains. In foliated samples, the tortuosity of the fracture observed in thin section was greater in arrester and divider geometries than in short transverse geometries.

58 GEOSCIENCES↗
DOE OSTI2021

Determination of Seven Organic Impurities in FD&C Yellow No. 6 by Ultra-High-Performance Liquid Chromatography

The U.S. Food and Drug Administration (FDA) batch certifies FD&C Yellow No. 6 (Y6) to ensure that the color additive meets requirements published in the Code of Federal Regulations (CFR), including specifications for seven organic manufacturing impurities consisting of two intermediates, a reaction by-product, an impurity originating from an intermediate, and three subsidiary colors. An ultra-high-performance liquid chromatography (UHPLC) method was developed and validated for determining seven organic impurities in Y6. For comparison, the currently used high-performance liquid chromatography (HPLC) method was also validated. The new UHPLC method uses a 1.7 µm particle size biphenyl column with aqueous ammonium formate and methanol as eluants. Analytes are identified by comparing their retention times and UV-visible spectra to those of reference standards. Calibration is performed in the presence of the dye matrix and analyte levels are determined from their peak areas. UHPLC and HPLC validation studies obtained linear calibration curves and excellent values for limits of detection, limits of quantitation, recovery, repeatability, and intermediate precision for all analytes. Survey analyses of 30 samples using the UHPLC and HPLC methods yielded results that were consistent within experimental error; however, the UHPLC method provided more accurate results for two analytes that coelute using the HPLC method. The UHPLC method satisfies the accuracy and precision requisites for routine certification of Y6. The UHPLC method is faster and generates less waste than the HPLC method while delivering better sensitivity, separation, and accuracy.

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

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Zinc isotope constraints on the cycling of carbon in the Bermuda mantle source

Volatile recycling and storage in the mantle transition zone (MTZ) is important for the refertilization of the upper mantle and is associated with the generation of high-µ (HIMU, where µ is 238 U/ 204 Pb) mantle. One way to probe the MTZ and the processes associated with mantle convection is to sample lavas that originate from the shallow mantle and were contaminated by upwelling from the MTZ, such as at the previously proposed shallow plume of Bermuda. Here we present the first δ 66 Zn isotopic compositions of Bermuda silica-undersaturated and silica-saturated lavas to explore the origin of the carbon-rich lithologies and the genesis of the large seamount found in the western North Atlantic Ocean. Contrasting with global δ 66 Zn data sets, our results (δ 66 Zn between 0.24 ± 0.04 and 0.41 ± 0.04) do not support direct sampling of recycled marine carbonates in the Bermuda HIMU mantle. Instead, we show that δ 66 Zn fractionation toward higher values is associated with magmatic processes and incorporation of carbon sourced from deep fluids associated with the formation of carbonatites. These carbon-rich fluids are likely sourced from the metasomatic reactions between the subducted cold slab of the Iapetus oceanic lithosphere ca. 500 Ma and the thickened continental lithospheric mantle of Pangea. Melting of this metasomatized mantle was triggered by the arrival of the Farallon slab to the eastern North American margin in the late Cenozoic via shallow convection.

Mazza, Sarah E. [Smith College, Northampton, MA (U↗
DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

The large thermal gradients and high subsurface temperatures of the western region of the U.S. hold great potential for the implementation of enhanced geothermal systems (EGS). The development of these potential EGS resources requires stimulation of the reservoir to enhance permeability and it has been widely reported that a substantial amount of water will be required should conventional hydraulic stimulation be used. This presents a huge challenge and a high risk to the geothermal development because the water stress1 in these areas is already high or extremely high. The use of foam, a gas/liquid mixture predominantly composed of gas, in fracturing is considered and explored in this project as a potential approach to address water concerns with hydraulic stimulation in the development of EGS. This project, led by Oak Ridge National Laboratory (ORNL) in collaboration with Temple University, was awarded in an open lab call in 2018, and was part of the DOE GTO waterless stimulation initiative. The goal of the project was to demonstrate the feasibility of foam fracturing for EGS development through two primary tasks: Task 1: Laboratory study of the effectiveness of foam fracturing for representative geological materials, including cyclic pressurization using foam (led by ORNL) and Task 2: High temperature foam material selection and characterization (led by Temple University). In FY19, ORNL finished the critical review on serval issues associated with foam fracturing and the implementation of the proposed tasks in a lab study (Wang, et al., 2019), and completed the foam fracturing testing using cement as a model material (Wang, et al, 2020a). The work at ORNL was geared up to develop a brand-new foam testing system in FY20. The purchase of main components for the new system was finished in the first half of the FY20. The assembly of the foam testing system and foam fracturing testing were completed in the second half of the FY20 (Wang, et al., 2021a). Task 1 required the development of a test system which can be used to perform hydraulic fracturing of geological specimens with both water and foamed liquids at pressure up to 6,000 psi (41.4 MPa). The system possesses several capabilities that conventional injection systems lack for hydraulic fracturing. In addition to its ability to generate foam with controlled quality, it is capable of cycling pressure levels between specified values up to frequencies of 50 Hz. The latter capability was developed to evaluate the hypothesis that cyclic loading of samples would produce enhanced fracturing. The system consists of two sections: one for foam generation and another for foam injection. The foam is generated through separate control and pressurization of liquid and gas phases with controlled flow rates. The injection section is equipped with a low-flow Coriolis flow that monitors the density of foam to ensure the injection is in the range of target foam quality2. Experimental results of foam fracturing are reported for cylindrical granite specimens using water and aqueous N2 foam as the fracturing fluids. All experiments were performed for unconfined conditions. The effects of injection mode (i.e., monotonic vs cyclic pressurization) on breakdown pressure and failure response sample were investigated using water alone as a fracturing fluid and foams with a range of compositions. It was found that in the case of monotonic injection, the breakdown pressure of granite specimens tended to be slightly higher when fracturing with foam. Additionally, with a foam quality of 90%, the water use can be reduced by 50 to 84%, depending on hole size. On the other hand, it was observed that the breakdown pressure can be brought down to 70% of the monotonic breakdown pressure by using low cycle fatigue. Finally, discussions are presented regarding injectivity and water use reduction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

SMC2021 Where to go in atomic world

'Graphene_CrSi.npy' contains a scanning transmission electron microscopy (STEM) movie from graphene monolayer. The movie is a sequence of atom-resolved images from the same sample region that undergoes chemical and structural transformations due to interaction with electron beam (which is used to perform imaging). 'topo_defects.npy' contains coordinates of some of the objects of interest (topological defects in graphene). It is a dictionary, where keys are frame numbers and values are xy coordinates. These objects are usually localized in relatively small areas of the image and we are interested in identifying them without having to scan an entire grid (which leads to fast degradation of the sample).

36 MATERIALS SCIENCE↗
DOE OSTI2023

An attribution of the low single-scattering albedo of biomass burning aerosol over the southeastern Atlantic

Abstract. Aerosol over the remote southeastern Atlantic is some of the most sunlight-absorbing aerosol on the planet: the in situ free-tropospheric single-scattering albedo at the 530 nm wavelength (SSA530 nm) ranges from 0.83 to 0.89 within ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) aircraft flights from late August–September. Here we seek to explain the low SSA. The SSA depends strongly on the black carbon (BC) number fraction, which ranges from 0.15 to 0.4. Low organic aerosol (OA)-to-BC mass ratios of 8–14 and modified combustion efficiency values >0.975 point indirectly to the dry, flame-efficient combustion of primarily grass fuels, with back trajectories ending in the miombo woodlands of Angola. The youngest aerosol, aged 4–5 d since emission, occupied the top half of a 5 km thick plume sampled directly west of Angola with a vertically consistent BC:ΔCO (carbon monoxide) ratio, indicating a homogenization of the source emissions. The younger aerosol, transported more quickly off of the continent by stronger winds, overlaid older, slower-moving aerosol with a larger mean particle size and fraction of BC-containing particles. This is consistent with ongoing gas condensation and the coagulation of smaller non-BC particles upon the BC-containing particles. The particle volumes and OA:BC mass ratios of the older aerosol were smaller, attributed primarily to evaporation following fragmentation, instead of dilution or thermodynamics. The CLARIFY (CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017) aircraft campaign sampled aerosols that had traveled further to reach the more remote Ascension Island. CLARIFY reported higher BC number fractions, lower OA:BC mass ratios, and lower SSA yet larger mass absorption coefficients compared to this study's. Values from one ORACLES 2017 flight, held midway to Ascension Island, are intermediate, confirming the long-range changes. Overall the data are most consistent with continuing oxidation through fragmentation releasing aerosols that subsequently enter the gas phase, reducing the OA mass, rather than evaporation through dilution or thermodynamics. The data support the following best fit: SSA530nm=0.801+0055⋅(OA:BC) (r=0.84). The fires of southern Africa emit approximately one-third of the world's carbon; the emitted aerosols are distinct from other regional smoke emissions, and their composition needs to be represented appropriately to realistically depict regional aerosol radiative effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Utility-Scale Solar, 2022 Edition: Empirical Trends in Deployment, Technology, Cost, Performance, PPA Pricing, and Value in the United States [Slides]

Berkeley Lab’s “Utility-Scale Solar, 2022 Edition” provides an overview of key trends in the U.S. market, with a focus on 2021. Highlights of this year’s update include: -A record of nearly 12.5 GWAC of new utility-scale PV capacity came online in 2021, bringing cumulative installed capacity to more than 51.3 GWAC across 44 states. -90% of all new utility-scale PV capacity added in 2021 uses single-axis tracking. -Median installed project costs declined to $\$1.35$/WAC (or $\$1.02$/WDC) in 2021. -Project-level capacity factors vary widely, from 9% to 35% (on an AC basis), with a sample median of 24%. The report explores drivers of this variation. -Utility-scale PV’s LCOE fell to $\$33$/MWh in 2021 ($\$27$/MWh if factoring in the federal investment tax credit, or ITC). -PPA prices have largely followed the decline in solar’s LCOE over time, but have recently stagnated and even moved slightly higher. Prices from a sample of recent contracts average around $\$20$/MWh (levelized) in the West and $\$30-40$/MWh elsewhere in the continental US. In 2021, solar’s average market value (defined in the report to include only energy and capacity value) rose by 55% to $\$47$/MWh and exceeded average wholesale prices in 13 of the 17 balancing authorities analyzed. -Adding battery storage is one way to increase the value of solar. Our public data file tracks metadata and PPA prices from 67 PV+battery hybrid projects that are already online or that have secured offtake arrangements. -At the end of 2021, there were at least 674 GW of utility-scale solar power capacity within the interconnection queues across the nation, 284 GW of which include batteries. For more information, and to explore related interactive data visualizations, go to utilityscalesolar.lbl.gov.

14 SOLAR ENERGY↗
DOE OSTI2024

High-resolution ion mobility based on traveling wave structures for lossless ion manipulation resolves hidden lipid features

Abstract High-resolution ion mobility (resolving power > 200) coupled with mass spectrometry (MS) is a powerful analytical tool for resolving isobars and isomers in complex samples. High-resolution ion mobility is capable of discerning additional structurally distinct features, which are not observed with conventional resolving power ion mobility (IM, resolving power ~ 50) techniques such as traveling wave IM and drift tube ion mobility (DTIM). DTIM in particular is considered to be the “gold standard” IM technique since collision cross section (CCS) values are directly obtained through a first-principles relationship, whereas traveling wave IM techniques require an additional calibration strategy to determine accurate CCS values. In this study, we aim to evaluate the separation capabilities of a traveling wave ion mobility structures for lossless ion manipulation platform integrated with mass spectrometry analysis (SLIM IM-MS) for both lipid isomer standards and complex lipid samples. A cross-platform investigation of seven subclass-specific lipid extracts examined by both DTIM-MS and SLIM IM-MS showed additional features were observed for all lipid extracts when examined under high resolving power IM conditions, with the number of CCS-aligned features that resolve into additional peaks from DTIM-MS to SLIM IM-MS analysis varying between 5 and 50%, depending on the specific lipid sub-class investigated. Lipid CCS values are obtained from SLIM IM ( TW(SLIM) CCS) through a two-step calibration procedure to align these measurements to within 2% average bias to reference values obtained via DTIM ( DT CCS). A total of 225 lipid features from seven lipid extracts are subsequently identified in the high resolving power IM analysis by a combination of accurate mass-to-charge, CCS, retention time, and linear mobility-mass correlations to curate a high-resolution IM lipid structural atlas. These results emphasize the high isomeric complexity present in lipidomic samples and underscore the need for multiple analytical stages of separation operated at high resolution. Graphical abstract

Reardon, Allison R. (ORCID:0000000165830134)↗
DOE OSTI2023

Logging Crew Attributes by Region in the Southeast USA

Urbanization, shrinking markets, and reduced forestry investment may affect harvesting efficiency in regions of the US South. To monitor these conditions, logging businesses have been tracked by surveys conducted by universities and trade associations. This project used a sampling approach coordinated with FIA utilization studies to sample logging crews based on a harvesting location. The approach was used to develop relationships among firm attributes and site attributes in six southeastern states (AL, GA, FL, NC, SC, and VA) from 2011 to 2018. The data included harvest attributes (location, harvest size and stand type) and logging firm attributes (production, crew labor, crew number, the number of machines by type, and machine age). For crew capital value, an equation was developed for this study using machine number and average machine age. The data from logging crews on 419 harvests were analyzed by region, harvest size, and stand type. Mean values for crew labor ranged from 3.1 to 7.1 workers. The average capital value per crew ranged from $\$220,000$ to $\$524,000$ per crew in the Coastal Plain with a narrower range in the Piedmont. In the Coastal Plain, higher productivity was detected for larger harvests and pine versus hardwood and mixed stands; however, in the Piedmont those trends were less obvious. Ratio of feller-bunchers, skidders and loaders were mostly 1:1:1 or 1:2:1 with 41% and 24% of samples, respectively. There were notable trends among Coastal Plain loggers regarding capital value and productivity with evidence supported by a production function. The differences in Piedmont (e.g., ownership size, market access, terrain, population density, etc.) may combine to limit daily production and labor productivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

NOvA 2024 official data release (26.61E20 neutrino + 12.5E20 antineutrino)

This data release corresponds to the Bayesian 2024 analysis of NOvA $\nu_e$ appearance and $\nu_{\mu}$ disappearance data, corresponding to analysis described in https://arxiv.org/abs/2509.04361. File `NOvA_2024_data_histograms.root` contains data histograms for all the NOvA data samples. Exposure: * neutrino-enhanced beam: 26.61E20 protons on target * antineutrino-enhanced beam: 12.5E20 protons on target External constraints: * ss2th12=0.851, dm21=7.53e-5 are fixed at 2019 PDG values, with negligible effect on NOvA predictions. * ss2th13 & dm32: * RCDB1D: 1D constraint from Daya Bay for ss2th13, 0.0851+/-0.0024 * RCDB2D: Correlated 2D constraint from Daya Bay on ss2th13 & dm32, available in their official 2023 data release: https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.130.161802#supplemental An additional file containing the predictions (`NOvA_2024_prediction_with_systs_histograms.root`) for all channels’ signal and background components computed at the NOvA best-fit oscillation parameters and systematic pull terms. Numu samples include the no-oscillation case as well. Please note that any fits performed with these histograms are not expected to exactly reproduce the official NOvA results as parameterizations of the numerous systematic uncertainties considered in the official fits are not included in this release. The zip file also contains the 2D credible interval contours, with details in a README.md

Sztuc, Artur [University Coll. London] (ORCID:0000↗
DOE OSTI2026

Data for Yield from Iowa’s first commercial miscanthus fields: implications of spatial variability for productivity and sustainability beyond research plots

This dataset contains biomass yield measurements and associated vegetation index data collected from commercial Miscanthus × giganteus fields in eastern Iowa during the 2022–2023 growing seasons. The data support the analyses presented in the article: “Yield From Iowa's First Commercial Miscanthus Fields: Implications of Spatial Variability for Productivity and Sustainability Beyond Research Plots.” We collected 105 ground-truth biomass samples from four mature commercial fields (>4 years old) covering 92.81 ha. Samples were taken from 3 m² quadrats that were hand-harvested in alignment with commercial harvest timing. Stem biomass (excluding leaves) was weighed, moisture-corrected, and converted to dry-matter yield expressed in Mg DM ha⁻¹. Sampling locations were selected to capture spatial variability visible in aerial imagery and were recorded using RTK GPS. Each biomass observation was paired with vegetation indices derived from high-resolution PlanetScope satellite imagery (3 m resolution). Images were acquired throughout the growing season, and indices were calculated to evaluate their ability to predict end-of-season biomass yield. Statistical and machine learning approaches were used to identify key predictors, and a linear regression model based on end-of-July Green Normalized Difference Vegetation Index (GNDVI) was developed and evaluated. This repository includes the data used in that modeling workflow. Management practices, economic data, full imagery time series, and additional methodological details are described in the associated publication and are not included here. The dataset consists of three comma-separated value (CSV) files: 1. Combine_Groundtruth_Yield_VI_22_23.csv This file contains ground-truth biomass yield measurements and associated key vegetation index values collected during the 2022 and 2023 growing seasons. Rows: 105 observations Columns: Year — Year of observation (2022 or 2023) Field — Field location identifier Sample_number — Unique sample identifier GNDVI_End_Jul — Green Normalized Difference Vegetation Index calculated at end of July GNDVI_End_Aug — Green Normalized Difference Vegetation Index calculated at end of August NDRE_End_Aug — Normalized Difference Red Edge index calculated at end of August Biomass_Stem_Yield_MgDM/ha — Measured stem biomass yield (megagrams dry matter per hectare) 2. trainData_GNDVI.csv This file contains the subset of observations used to train the predictive relationship between July GNDVI and biomass yield. Rows: 76 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Stem_Yield_MgDM/ha — Observed stem biomass yield (Mg DM ha⁻¹) 3. testData_GNDVI.csv This file contains the test dataset used to evaluate model performance. Rows: 29 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Predicted_Yield_MgDM/ha — Model-predicted stem biomass yield (Mg DM ha⁻¹) Observed_Yield_MgDM/ha — Measured stem biomass yield (Mg DM ha⁻¹)

Potential yield, yield gap, in-field management, y↗
DOE OSTI2021

Exploring DFT+U parameter space with a Bayesian calibration assisted by Markov chain Monte Carlo sampling

Abstract The density-functional theory is widely used to predict the physical properties of materials. However, it usually fails for strongly correlated materials. A popular solution is to use the Hubbard correction to treat strongly correlated electronic states. Unfortunately, the values of the Hubbard U and J parameters are initially unknown, and they can vary from one material to another. In this semi-empirical study, we explore the U and J parameter space of a group of iron-based compounds to simultaneously improve the prediction of physical properties (volume, magnetic moment, and bandgap). We used a Bayesian calibration assisted by Markov chain Monte Carlo sampling for three different exchange-correlation functionals (LDA, PBE, and PBEsol). We found that LDA requires the largest U correction. PBE has the smallest standard deviation and its U and J parameters are the most transferable to other iron-based compounds. Lastly, PBE predicts lattice parameters reasonably well without the Hubbard correction.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A self-referenced in-situ arrival time monitor for X-ray free-electron lasers

We present a novel, highly versatile, and self-referenced arrival time monitor for measuring the femtosecond time delay between a hard X-ray pulse from a free-electron laser and an optical laser pulse, measured directly on the same sample used for pump-probe experiments. Two chirped and picosecond long optical supercontinuum pulses traverse the sample with a mutually fixed time delay of 970 fs, while a femtosecond X-ray pulse arrives at an instant in between both pulses. Behind the sample the supercontinuum pulses are temporally overlapped to yield near-perfect destructive interference in the absence of the X-ray pulse. Stimulation of the sample with an X-ray pulse delivers non-zero contributions at certain optical wavelengths, which serve as a measure of the relative arrival time of the X-ray pulse with an accuracy of better than 25 fs. We find an excellent agreement of our monitor with the existing timing diagnostics at the SACLA XFEL with a Pearson correlation value of 0.98. We demonstrate a high sensitivity to measure X-ray pulses with pulse energies as low as 30 μJ. Using a free-flowing liquid jet as interaction sample ensures the full replacement of the sample volume for each X-ray/optical event, thus enabling its utility even at MHz repetition rate XFEL sources.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

CHESS 2025: Leaf Area Index (LAI) for meadow, shrub, tree, and understory vegetation

This dataset contains Leaf Area Index (LAI) measurements made as part of the Colorado Headwaters Ecological Spectroscopy Study (CHESS) during June and July of 2025. Data were collected in the Upper Gunnison Basin, Colorado, across three study domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). Field observations of LAI were collected within 72 hours of airborne data collection by the National Ecological Observatory Network’s Aerial Observation Platform (NEON AOP). The NEON AOP collected waveform LiDAR (Light Detection and Ranging) and imaging spectrometer data in 426 spectral bands from the visible to shortwave infrared. LAI measurements were collected using the LICOR LAI-2200C Plant Canopy Analyzer following protocols outlined in the instrument manual (LI-COR 2019). Sampling targeted four distinct vegetation types: meadows, shrubs, trees, and aspen forest understory. We have archived data separately by site type because different field methods were used for each. At meadow sites, measurements were made at the four corners of 1m x 1m plots, with the instrument moving inward toward the center of the plot. At shrub sites, we measured the canopies of individual shrubs. At tree sites, we made measurements within a 10m x 10m subplot centered around a focal tree, with 30 observations taken on a regular grid. At aspen understory sites, we measured overstory trees following the tree protocol and understory herbaceous vegetation following the meadow protocol. All measurements included above-canopy (A) and below-canopy (B) readings, with specific protocols for scattering correction measurements in direct-sun conditions. Data were processed using the R package `rlai` (Worsham 2025). This package includes functions to calculate LAI, gap fraction, apparent clumping factor (Ω), scattering correction, and other canopy metrics. Package contents: Full file descriptions appear in ‘flmd.csv’. Files named according to the convention ‘lai_*_summary_data_cleaned.csv’ contain summary values of LAI, apparent clumping factor (Ωapp), and scattering correction factors for each site. These are the analysis-ready products that most data users will work with. Files named ‘lai_*_metadata_cleaned.csv’ contain additional site-level observations made during field collection. We have also archived intermediate and supplementary data for users who wish to check our processing approach or apply alternative methods. ‘raw_lai_2200C.zip’ contains the raw files as read from the LI-COR instrument, with no processing applied, in TXT format. The zip archive contains subdirectories by site type, which are further subdivided by sampling area. Filenames correspond to the sampling site number. ‘intermediate_results.zip’ contains detailed output from the processing routines, in JSON format. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘scattering_correction_logs.zip’ contains logfiles from the implementation of Kobayashi et al.'s (2013) scattering correction algorithm. The logfiles report values of several parameters at each iteration of the algorithm, as the model converges toward a stable solution. They are intended for users who want to verify scattering correction performance. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘spot_checks.csv’ reports LAI and other values for a small number of files processed with LI-COR FV2200 software (LI-COR 2013) using the same control parameters as in our R-based approach. Additional metadata are provided in a data dictionary describing column names and definitions (dd.csv), and in a file-level metadata file (flmd.csv). All zip files can be expanded with common archive utilities. TXT, CSV, and JSON files can be ingested into R or Python computing environments or read in common text editor utilities. Geospatial information: Geospatial data for mapping measurement site locations are in the files CHESS_polygons_lai_UTM.geojson, CHESS_polygons_shrub_UTM.geojson, and CHESS_polygons_meadow_UTM.geojson in the companion geospatial package for the 2025 CHESS campaign, ‘CHESS 2025: Location data for field observations and sampling’ (Henderson et al., 2026). CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. * Todorov and Worsham are co–first authors.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2022

RAPID UNSEPARATED RARE EARTH ELEMENT ANALYSES BY ISOTOPE DILUTION MULTICOLLECTOR INDUCTIVELY COUPLED MASS SPECTROMETRY (ID-MC-ICP-MS)

We present a method for analyzing rare earth elements (REEs) by isotope dilution without separation of individual elements. Utilizing a Nu Plasma 2 multicollector inductively coupled plasma mass spectrometer (MC-ICP-MS) equipped with an Aridus II desolvation nebulizer, sample introduction tuning parameters (torch position, gas flows) were investigated and used to minimize oxide formation (e.g., NdO/Nd <0.05%) and remove the need to analyze the REEs as separated “cuts.” Ratios for spike/sample isotopes in mixed synthetic standards containing varying amounts of Ce, Nd, Sm, Eu, Gd, Dy, Er, and Yb were measured accurately within 1% of IUPAC values. La and Lu in the standards were more than 2% off from expected values due to more demanding interference corrections. A USGS BCR-2 rock standard was also processed with flux fusion, iron coprecipitation and a single column pass bulk separation followed by analysis. The concentrations for Ce, Nd, Sm, Eu, Dy, Er, and Yb fell within 2% of consensus values. A United States Geological Survey (USGS) BIR-1A rock standard was also processed in triplicate by hydrofluoric and nitric acid digestion, yielding rare earth concentrations within 3% of consensus values with 0.3% relative standard deviation (n=3). This demonstrates the ability to produce high quality REE patterns of comparable quality to traditional methods involving labor and time intensive REE separation, with sample preparation and analytical times that are comparable to analyses against isotope dilution by single-collector ICP-MS.

Shen, Steve D.↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗