Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Upper regularity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

A nonsmooth nonconvex optimization algorithm for two-stage optimization problems

An optimization algorithm for a group of nonsmooth nonconvex problems inspired by two-stage stochastic programming problems is proposed. The main challenges for these problems include (1) the problems lack the popular lower-type properties such as prox-regularity assumed in many nonsmooth nonconvex optimization algorithms, (2) the objective can not be analytically expressed and (3) the evaluation of function values and subgradients are computationally expensive. To address these challenges, this report first examines the properties that exist in many two-stage problems, specifically upper-C 2 objectives. Then, we show that quadratic penalty method for securityconstrained alternating current optimal power flow (SCACOPF) contingency problems can make the contingency solution functions upper-C 2 . Based on these observations, a simplified bundle algorithm that bears similarity to sequential quadratic programming (SQP) method is proposed. It is more efficient in implementation and computation compared to conventional bundle methods. Global convergence analysis of the algorithm is presented under novel and reasonable assumptions. The proposed algorithm therefore fills the gap of theoretical convergence for smoothed SCACOPF problems. The inconsistency that might arise in our treatment of the constraints are addressed through a penalty algorithm whose convergence analysis is also provided. Finally, theoretical capabilities and numerical performance of the algorithm are demonstrated through numerical examples.

97 MATHEMATICS AND COMPUTING↗

Impacts, processes and projections of the quasi-biennial oscillation

In the tropical stratosphere, deep layers of eastward and westward winds encircle the globe and descend regularly from the upper stratosphere to the tropical tropopause. With a complete cycle typically lasting almost 2.5 years, this quasi-biennial oscillation (QBO) is arguably the most predictable mode of atmospheric variability that is not linked to the changing seasons. Here, the QBO affects climate phenomena outside the tropical stratosphere, including ozone transport, the North Atlantic Oscillation and the Madden–Julian Oscillation, and its high predictability could enable better forecasts of these phenomena if models can accurately represent the coupling processes. Climate and forecasting models are increasingly able to simulate stratospheric oscillations resembling the QBO, but exhibit common systematic errors such as weak amplitude in the lowermost tropical stratosphere. Uncertainties about the waves that force the oscillation, particularly the momentum fluxes from small-scale gravity waves excited by deep convection, make its simulation challenging. Improved representation of the processes governing the QBO is expected to lead to better forecasts of the oscillation and its impacts, increased understanding of unusual events such as the two QBO disruptions observed since 2016, and more reliable future projections of QBO behaviour under climate change.

54 ENVIRONMENTAL SCIENCES↗

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

New particle formation in the remote marine boundary layer

Marine low clouds play an important role in the climate system, and their properties are sensitive to cloud condensation nuclei concentrations. While new particle formation represents a major source of cloud condensation nuclei globally, the prevailing view is that new particle formation rarely occurs in remote marine boundary layer over open oceans. Here we present evidence of the regular and frequent occurrence of new particle formation in the upper part of remote marine boundary layer following cold front passages. The new particle formation is facilitated by a combination of efficient removal of existing particles by precipitation, cold air temperatures, vertical transport of reactive gases from the ocean surface, and high actinic fluxes in a broken cloud field. The newly formed particles subsequently grow and contribute substantially to cloud condensation nuclei in the remote marine boundary layer and thereby impact marine low clouds.

54 ENVIRONMENTAL SCIENCES↗

The Short Life of Upvalley Wind in a High‐Altitude Valley in the Colorado Rocky Mountains

Thermally driven upvalley (UV) wind in the upper East River Valley in the Colorado Rocky Mountains often unexpectedly stops in midmorning and reverses back to downvalley (DV) wind. We use a comprehensive observational data set for a nearly two‐year long period to analyze the wind system and boundary layer evolution in this high‐altitude valley and determine the reason for this early wind reversal. Days with short UV wind predominantly occur during the warm season when the valley floor is free of snow and the convective boundary layer (CBL) grows well above the height of the surrounding ridges. UV wind persists throughout the day only on a few days during the warm season. We link differences in valley wind evolution to wind direction at upper levels at and above ridge height and propose forced channeling mechanisms to describe coupling between valley and upper‐level wind when the CBL grows above ridge height. The frequency distribution of upper‐level wind direction is such that channeling in the DV direction is favored, which explains the predominance of days with short UV wind. The deep CBL is supported by the presence of a deep weakly stably stratified residual layer with high aerosol content, which is regularly present over the mountain range during the warm season. On days when the CBL does not grow above ridge height, for example, when the valley floor is covered by snow, thermally driven UV wind is able to persist throughout the day independent of upper‐level wind direction.

54 ENVIRONMENTAL SCIENCES↗

Spectral Threshold for Extremal Cyclic Edge-Connectivity

In this report, the cyclic edge-connectivity of a graph G is the least k such that there exists a set of k edges whose removal disconnects G into components where every component contains a cycle. We show that for graphs of minimum degree at least 3 and girth g at least 4, the cyclic edge-connectivity is bounded above by (Δ-2)g where Δ is the maximum degree. We then prove that if the second eigenvalue of the adjacency matrix of a d-regular graph of girth g ≥ 4 is sufficiently small, then the cyclic edge-connectivity is (d-2)g, providing a spectral condition for when this upper bound on cyclic edge-connectivity is tight.

97 MATHEMATICS AND COMPUTING↗

Hydrokinetic energy conversion using flow induced oscillations of single-cylinder with large passive turbulence control

Various types of flow-induced oscillations (FIOs) have been implemented in development of marine hydrokinetic (MHK) energy converters. With passive turbulence control (PTC), energy harvesting starts at a flow speed of about 0.5 m/s. However, there is worldwide MHK energy available in even slower currents. In the present study, the effect of damping on FIO and power extraction is investigated for a converter with large turbulence stimulation (PTC) consisting of straight strips with a height of 15% of the cylinder diameter and placed symmetrically on the cylinder surface. The oscillating amplitude decreases, as the damping ratio increases, with unchanged sinusoidal pattern of the displacement time-history. The frequency ratio is also affected by damping especially in the VIV initial branch and transition region between VIV and galloping. An important flow characteristic of the large-PTC cylinder is that a recirculation region is formed behind the PTC, causing appreciable disturbance to the flow past the cylinder. Power can be harvested in the whole FIO range and the harnessed power maximum appears at the largest inflow velocity tested. However, the optimum of harnessing efficiency is located at the beginning of the VIV upper branch. The gap between VIV and galloping is bridged when large PTC is used, eliminating the drop in power and efficiency even at higher damping, which would be a weakness of regular-PTC cylinder for energy harvesting. Finally, the mechanism behind the variation of harnessing efficiency with inflow velocity and damping ratio is revealed, and the optimality criterion for the converter design is discussed.

16 TIDAL AND WAVE POWER↗

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

This dataset provides Level 1 (L1) and Level 3 (L3) orthorectified Red-Green-Blue (RGB) imagery collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). This high-resolution imagery is a photographic record of red, green, and blue visible light from sunlight reflected off of the Earth’s surface. The data comprise full-color images of the ground surface and are primarily intended to provide context to imaging spectroscopy and light detection and ranging (LiDAR) data. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. RGB images were acquired using the PhaseOne IXM-RS150F high-resolution digital camera onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). The package data include both an L1 product comprising one camera frame per file and an L3 mosaic aligned to the Universal Transverse Mercator (UTM) Zone 13N grid and the World Geodetic System (WGS) 84 projection. Both products are provided in geotif (.tif) format at 0.1 m ground resolution. The bulk of the imagery was collected during the main CHESS field campaign from June 13 to July 15, 2025. Additional images of a portion of the Upper Taylor (UPTA) domain were collected on September 18, 2025, to fill gaps in imagery identified after the main campaign was complete. RGB camera imagery is not radiometrically calibrated, and therefore pixel values should not be exploited for scientific analysis. Pixel values have undergone a manual adjustment to enhance feature identification. The imagery is rigorously geolocated which does allow for reliable geometric information to be retrieved. To generate the orthorectified imagery, the NEON AOP camera captured visible spectrum in red, green, and blue bands. The raw images were then processed using NEON’s camera orthorectification workflow. A boresight calibration flight was made to build a complete camera, distortion, and alignment model. Color balance/white balance and exposure correction were applied to the raw RGB images. The corrected images were orthorectified by ray-tracing image pixels to a lidar-derived digital surface model (DSM) mesh using the refined camera model, outputting orthorectified raster pixels on a regular grid. Flightline-level data were mosaicked by selecting per-pixel contributions from overlapping orthorectified images using line-of-sight (LOS) zenith angle minimization to reduce edge distortions. 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.

2018 NEON and 2025 CHESS Campaigns↗

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗

Infinite quantum signal processing

Quantum signal processing (QSP) represents a real scalar polynomial of degree d using a product of unitary matrices of size 2 × 2 , parameterized by ( d + 1 ) real numbers called the phase factors. This innovative representation of polynomials has a wide range of applications in quantum computation. When the polynomial of interest is obtained by truncating an infinite polynomial series, a natural question is whether the phase factors have a well defined limit as the degree d → ∞ . While the phase factors are generally not unique, we find that there exists a consistent choice of parameterization so that the limit is well defined in the ℓ 1 space. This generalization of QSP, called the infinite quantum signal processing, can be used to represent a large class of non-polynomial functions. Our analysis reveals a surprising connection between the regularity of the target function and the decay properties of the phase factors. Our analysis also inspires a very simple and efficient algorithm to approximately compute the phase factors in the ℓ 1 space. The algorithm uses only double precision arithmetic operations, and provably converges when the ℓ 1 norm of the Chebyshev coefficients of the target function is upper bounded by a constant that is independent of d . This is also the first numerically stable algorithm for finding phase factors with provable performance guarantees in the limit d → ∞ .

Dong, Yulong [Department of Mathematics, Universit↗

Energy Efficiency and Performance Evaluation of Self-Contained, Medium Temperature, Reach-In Refrigerated Display Cases

This project was part of an effort by Commonwealth Edison Company (ComEd) to evaluate the energy and peak demand saving potential of emerging technologies in the Chicago area. This document focuses on the assessment of energy-efficient, medium-temperature, self-contained refrigerated display cases utilizing environmentally friendly refrigerants. The results of this evaluation will be considered by CLEAResult to develop a new energy efficiency rebate measure for ComEd's incentive programs. This rebate measure will become an addition to the Technical Reference Manual (TRM). In 2016, the United Nations passed the Kigali Montreal Protocol Amendment which placed restrictions on certain types of refrigerants. In compliance with this amendment, the US Environmental Protection Agency (EPA) placed a ban on the manufacture of refrigeration systems using hydrofluorocarbons including R134a starting in January 2020. Although the ban has halted manufacture, the EPA continues to allow the use of these refrigerants. Therefore, it is critical to provide incentives for replacing these refrigerants with other environmentally friendly and energy-efficient alternatives. The energy-efficient refrigerator cases evaluated here (referred to as EE Case A and B) contain environmentally-friendly refrigerants in compliance with the EPA hydrofluorocarbon ban. These consist of natural refrigerant propane (R290), and HFC drop-in hydrofluoroolefin R513a, respectively. These cases also contain other energy-efficient components including efficient lighting, more robust evaporator and condenser fans and different-sized heat exchangers. EE Case A is also built with materials that better insulate the case, which improves energy efficiency. To ascertain the energy efficiency contribution of these design components, the consumption of the evaporator and condenser fan motors, compressor, and lighting/controller were evaluated individually. The medium-temperature, self-contained reach-in refrigerated display case was selected due to its widespread use in convenience stores and small supermarkets. Self-contained refrigeration has also seen increased use in restaurants due to curbside pickup during the COVID-19 pandemic. For this study, the refrigerated display cases' performance was evaluated in a controlled environmental chamber at representative indoor dry-bulb and humidity conditions found in supermarkets within ComEd's service territory climate zone. The test method used in this project was foundationally inspired by the ANSI/ASHRAE 72-2018 method of testing. However, modifications to the ANSI/ASHRAE test methodology were implemented to better represent customer operation of the units. In addition to the indoor supermarket conditions, the cases were also evaluated at the "upper target," or environmental conditions used in ANSI/ASHRAE 72-2018. The case total power and case components were metered to obtain their daily power (kW) and energy consumption(kWh). The cases were filled with thermal filler mass to replicate thermal mass of product loading. Additionally, product simulators were used to provide product temperature information. Door actuators were mounted to each of the cases' three doors to replicate regular door openings and effects of shopper traffic.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Search for the Anomalous Events Detected by ANITA Using the Pierre Auger Observatory

A dedicated search for upward-going air showers at zenith angles exceeding 110° and energies E > 0.1 EeV has been performed using the Fluorescence Detector of the Pierre Auger Observatory. The search is motivated by two “anomalous” radio pulses observed by the ANITA flights I and III that appear inconsistent with the standard model of particle physics. Using simulations of both regular cosmic-ray showers and upward-going events, a selection procedure has been defined to separate potential upward-going candidate events and the corresponding exposure has been calculated in the energy range [0.1–33] EeV. One event has been found in the search period between January 1, 2004, and December 31, 2018, consistent with an expected background of 0.27 ± 0.12 events from misreconstructed cosmic-ray showers. This translates to an upper bound on the integral flux of ( 7.2 ± 0.2 ) × 10 − 21 cm − 2 sr − 1 y − 1 and ( 3.6 ± 0.2 ) × 10 − 20 cm − 2 sr − 1 y − 1 for an E − 1 and E − 2 spectrum, respectively. An upward-going flux of showers normalized to the ANITA observations is shown to predict over 34 events for an E − 3 spectrum and over 8.1 events for a conservative E − 5 spectrum, in strong disagreement with the interpretation of the anomalous events as upward-going showers. Published by the American Physical Society 2025

Abdul Halim, A.↗

Free energy and specific heat near a quantum critical point of a metal

We analyze the free energy and the specific heat for fermions interacting with a gapless boson at a quantum-critical point (QCP) in a metal. We use the Luttinger-Ward-Eliashberg formula for the free energy in the normal state, which includes contributions from bosons, fermions, and their interaction, all expressed via fully dressed fermionic and bosonic propagators. The sum of the last two contributions is the free energy F γ of an effective low-energy model of fermions with boson-mediated dynamical 4-fermion interaction V(Ω m ) ∝ 1/|Ω m | γ (the γ model). This purely electronic model has been used to analyze the interplay between non-Fermi liquid behavior and pairing near a QCP, which are both independent of the upper energy cutoff Λ. However, the specific heat C γ (T), obtained from F γ , does depend on Λ. We argue that this dependence is spurious and cancels out, once we include the contribution from bosons. We further argue that the full C(T) is the sum of the contribution from free fermions and the one from a critical boson, with the fully dressed propagator, other terms cancel out. We compare the full C(T) with the C γ (T), obtained using recently proposed regularization of F γ . We argue that for γ < 1, the full C(T) and the regularized C γ (T) differ by a γ-dependent prefactor, while for γ > 1, the full C(T) and C γ (T) differ by the positive contribution from free massless fermions (a positive constant for the electron-phonon case γ = 2). For these γ, C γ (T) is negative, but the full C(T) is positive. We argue that only the full C(T) matters as the positive and the negative contributions originate from the term in C(T) which contains the fully dressed bosonic propagator. Here, we then argue that the normal state remains stable until the pairing instability develops.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Resolving extreme jet substructure

We study the effectiveness of theoretically-motivated high-level jet observables in the extreme context of jets with a large number of hard sub-jets (up to N = 8). Previous studies indicate that high-level observables are powerful, interpretable tools to probe jet substructure for N ≤ 3 hard sub-jets, but that deep neural networks trained on low-level jet constituents match or slightly exceed their performance. We extend this work for up to N = 8 hard sub-jets, using deep particle-flow networks (PFNs) and Transformer based networks to estimate a loose upper bound on the classification performance. A fully-connected neural network operating on a standard set of high-level jet observables, 135 N-subjetiness observables and jet mass, reach classification accuracy of 86.90%, but fall short of the PFN and Transformer models, which reach classification accuracies of 89.19% and 91.27% respectively, suggesting that the constituent networks utilize information not captured by the set of high-level observables. We then identify additional high-level observables which are able to narrow this gap, and utilize LASSO regularization for feature selection to identify and rank the most relevant observables and provide further insights into the learning strategies used by the constituent-based neural networks. The final model contains only 31 high-level observables and is able to match the performance of the PFN and approximate the performance of the Transformer model to within 2%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and testing of residential micro-CHP powered by opposed piston engine (CRADA Final Report)

A micro–combined heat and power (mCHP) prototype powered by innovative opposed piston engine technology was developed to simultaneously provide electricity and heat to residential or light commercial buildings. The mCHP prototype targeted at residential applications includes an opposed-piston four-stroke (OP4S) engine, generator, rectifier, inverter, battery energy storage system, 52 gal water tank, and application accessories for hot water supply and space heating. The OP4S engine can use renewable or regular natural gas, as well as hydrogen, to generate mechanical power and waste heat in form of hot coolant and exhaust gas simultaneously. The waste heat is recovered and stored in the water tank and can be used as a regular hot water supply and/or for space heating application. The tests show that the mCHP prototype enabled power outputs in the range of 3.2 –7.4 kW with up to 26.4% of AC electricity efficiency and up to 93.1% of the overall mCHP efficiency under stoichiometric combustion modes λ=~1.0. The mCHP was also run under lean combustion mode conditions at λ=~1.3. The lean mode operation enables more than 30% improvement in electrical energy efficiency. The maximum AC efficiency of the lean combustion mode attained was 35.2%, with the engine efficiency is approaching 40%. The exceptional electrical efficiency breaks the typical upper boundary of 30% for ICE-based mCHP. The engine exhaust temperatures in the lean modes are substantially less than in the stoichiometric modes. Moreover, the lean cases achieve high overall mCHP efficiencies: the overall mCHP efficiencies are all greater than 93%. Considering the mCHP prototype can achieve low-cost, flexible matching of thermal and electrical loads through reducing the complexity of distribution and installation, and high efficiency the novel technology will promote mCHP acceptance in the US residential and light commercial markets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Multi-year assessment of the impact of ship-borne radiosonde observations on polar WRF forecasts in the Arctic

Abstract To compensate for the lack of conventional observations over the Arctic Ocean, ship-borne radiosonde observations have been regularly carried out during summer Arctic expeditions and the observed data have been broadcast via the global telecommunication system since 2017. With these data obtained over the data-sparse Arctic Ocean, observing system experiments were carried out using a polar-optimized version of the Weather Research and Forecasting (WRF) model and the WRF Data Assimilation (WRFDA) system to investigate their effects on analyses and forecasts over the Arctic. The results of verification against reanalysis data reveal: (1) DA effects on analyses and forecasts; (2) the reason for the year-to-year variability of DA effects; and (3) the possible role of upper-level potential vorticity in delayed DA effects. The overall assimilation effects of the extra data on the analyses and forecasts over the Arctic are positive. Initially, the DA effects are the most apparent in the temperature variables in the middle/lower troposphere, which spread to the wind variables in the upper troposphere. The effects decrease with time but reappear after approximately 120 h, even in the 240-h forecasts. The effects on forecasts vary depending on the proximity of the radiosonde observation locations to the high synoptic variability. The upper-level potential vorticity is known to play an important role in the development of Arctic cyclones, and it is suggested as a possible explanation for the delayed DA effects after about 120 h.

Geology↗

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↗

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗