SEARCH · Search NASA
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.
This dataset consists of the ultra high-resolution (km-scale) global gravity wave momentum fluxes, computed from the ECMWF IFS Experimental Nature Run at 1-km (XNR1K) simulation data, using Helmholtz decomposition. The raw model output from the ECMWF XNR1K data is composed of simulations corresponding to Nov 2018 - Feb 2019.
The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard ΛCDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Lyα flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.
Liquid argon time projection chambers, (LArTPCs), are particle detectors used to collect ionization charge information from particle trajectories, facilitating detailed particle track analysis. They are currently used as particle detectors in major physics projects such as the deep underground neutrino experiment (DUNE), to detect and study the nature of the elusive neutrino particle. TinyTPC is a small scale LArTPC featuring a pixelated readout system (LArPix) that we will use to study liquid argon doping. We expect this doping to enhance the resolution of LAr detectors, especially for low energy particles below 10 MeV which would expand the capabilities of currently running experiments. Housed within a cryostat filled with liquid argon, the TinyTPC will rely on a high and a low voltage system to collect data. In preparation for deployment we found and resolved issues in the HV and LV systems and we determined optimal running conditions in a test vessel. This presentation will go over the detector commissioning process that enabled data taking with the TinyTPC.
We report that following the initial development of High-Speed Counter-Current Chromatography (HSCCC) in the 1960s, several studies have explored its applicability in the separation of rare earth elements (REEs). More recently, however, HSCCC publications have transitioned towards the separation of natural products or pharmaceuticals, leaving the application for REEs largely unexplored from a practical standpoint. Herein, we expand upon prior work in this field by evaluating the suitability of HSCCC to separation of a subset of non-radioactive REEs (Nd, Sm, Eu, Tb, and Y) at 10 -4 mol levels using di-(2-ethylhexyl)phosphoric acid (HDEHP) in n-heptane as the stationary phase and hydrochloric acid as the mobile phase. First, the effect of flow rate on the stationary phase volume retention ratio and resolution of Nd/Sm/Eu subgroup was evaluated followed by optimization of step-gradient elution profiles resulting in additional recovery of Tb and Y within a seven-hour window. The five REEs were separated at the baseline resolution level or above. Elution profiles obtained from multiple runs across two independently operated columns and across independent runs were cross analyzed. Reproducibility in elution profiles point to future applications in radioelement separation chemistry, where both chemical and radiochemical purity are of importance.
The thermal stratification phenomenon plays a crucial role in the safety of various nuclear systems, including the Gen-III + Light Water Reactors (LWR) and the Gen-IV reactors. The phenomenon is of particular importance for the pool-type Sodium-cooled Fast Reactors (SFRs) because it may cause neutronic and thermal-hydraulic in- stabilities in the reactor core, or lead to damages of both the reactor vessel and in-vessel components due to the growth of thermal fatigue cracking. More significantly, thermal stratification could impede the establishment of the natural circulation during accidental scenarios and introduce uncertainties to the core safety of SFRs. Efforts for modeling of the thermal stratification in SFRs have been made for decades to prevent or mitigate the damage caused by the phenomenon. This paper gives a review of the advances that have been made in recent 10 years on the computational modeling methods for thermal stratification phenomenon in SFRs. These methods can be generally drawn into two categories. The first one is the system-level methods which provide fast-running but approximate calculations, and the second one is the CFD methods which provide high-resolution calculations at high computational expense. Finally, after introducing the efforts that have been made to improve the one- dimensional (1-D) models, the paper envisioned the possible research directions that could be pursued to enhance the modeling of thermal stratification in the near future.
An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.
Focal Area: Focal Area 3: Insight gleaned from complex simulated data using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.
ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.
In this conducted research, optical fiber sensors are used to measure low strength alternating magnetic fields. Various fiber sensor configurations are tested and investigated to demonstrate sensing capabilities at different field magnitudes and frequencies. Distributed acoustic sensing fibers (DAS) have been largely studied and documented across a variety of applications and sensing systems. This research uses the DAS technology in tandem with magnetostrictive materials to create a distributed multi-material optical fiber magnetic sensor. Magnetic sensing has high demand across different fields and often runs into challenges of extreme environments including high temperature, corrosion, and areas with poor accessibility. The robust and distributed nature of optical fiber sensors which can be cheaply produced for long lengths is an attractive option over other single point magnetic sensors. In down hole applications specifically, having a distributed sensor able to be deployed easily and over long distances for magnetic sensing would be a large improvement to bulkier traditional magnetometers. In the conducted study, different magnetostrictive materials are implemented in distributed optical fiber sensors to analyze and compare the effective sensitivity and potential commercial viability. Nickel, galfenol alloy, and MetGlas alloy inclusions are drawn into fused silica optical fibers with Bragg gratings inscribed later on for DAS capability. Each was investigated for its response to varying AC magnetic fields to determine relative sensitivity and resolution for distributed magnetic field sensing.
Here, we present the near-Maximal Algorithm for Poisson-disk Sampling (nMAPS) to generate point distributions for variable resolution Delaunay triangular and tetrahedral meshes in two and three-dimensions, respectively. nMAPS consists of two principal stages. In the first stage, an initial point distribution is produced using a cell-based rejection algorithm. In the second stage, holes in the sample are detected using an efficient background grid and filled in to obtain a near-maximal covering. Extensive testing shows that nMAPS generates a variable resolution mesh in linear run time with the number of accepted points. We demonstrate nMAPS capabilities by meshing three-dimensional discrete fracture networks (DFN) and the surrounding volume. The discretized boundaries of the fractures, which are represented as planar polygons, are used as the seed of 2D-nMAPS to produce a conforming Delaunay triangulation. The combined mesh of the DFN is used as the seed for 3D-nMAPS, which produces conforming Delaunay tetrahedra surrounding the network. Under a set of conditions that naturally arise in maximal Poisson-disk samples and are satisfied by nMAPS, the two-dimensional Delaunay triangulations are guaranteed to only have well-behaved triangular faces. While nMAPS does not provide triangulation quality bounds in more than two dimensions, we found that low-quality tetrahedra in 3D are infrequent, can be readily detected and removed, and a high-quality balanced mesh is produced.
We use spin torque ferromagnetic resonance and ferromagnetic-resonance-driven spin pumping to detect spin-charge interconversion at room temperature in heterostructure devices that interface an archetypal Dirac semimetal, Cd 3 As 2 , with a metallic ferromagnet, Ni 0.80 Fe 0.20 (permalloy). Angle-resolved photoemission directly reveals the Dirac-semimetal nature of the samples prior to device fabrication and high-resolution transmission electron microscopy is used to characterize the crystalline structure and the relevant heterointerfaces. We find that the spin-charge interconversion efficiency in Cd 3 As 2 /permalloy heterostructures is comparable to that in heavy metals and that it is enhanced by the presence of an interfacial oxide. Spin torque ferromagnetic resonance measurements reveal an in-plane spin polarization regardless of an oxidized or pristine interface. We discuss the underlying mechanisms for spin-charge interconversion by comparing our results with first principles calculations and conclude that extrinsic mechanisms dominate the observed phenomena. Our results indicate a need for caution in interpretations of spin-transport and spin-charge conversion experiments in Cd 3 As 2 devices that seek to invoke the role of topological Dirac and Fermi arc states.
The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.
Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.
Nitrate concentrations in streams and rivers in the Midwestern United States are often elevated, reflecting the predominance of agriculture in the surrounding landscape. Recent advances in technology, including surface water drones and more precise sensors, provide opportunities to investigate nitrate dynamics with high spatial and temporal resolution. We deployed an aquatic drone, the AquaBOT, in a sixth-order, agriculturally influenced river to examine longitudinal patterns in water quality. Our goal was to measure the spatial and temporal heterogeneity in nitrate and nitrate removal processes and determine the influence of tributary inputs on main stem chemistry. We navigated the drone along a 12-km reach of the Des Moines River (Iowa, USA) nine times between June 2021 and August 2022. Across the deployments, mean nitrate concentration was positively related to discharge and was nearly two orders of magnitude higher in spring than summer. We observed contrasting patterns in main stem nitrate, which decreased downstream during some runs (e.g., 3.1–2.7 mg N L −1 in June 2021), demonstrating net nitrate uptake along the reach, and remained constant on other dates. Similarly, tributaries to the Des Moines had a varied influence on riverine nitrate. Tributaries either increased or decreased main stem nitrate concentrations depending on the tributary and the date. Nitrate removal rates were spatially and temporally variable but showed some consistency at the subreach (2 km) scale, with two subreaches often showing elevated rates of nitrate removal across dates. In conclusion, our study reveals nuanced heterogeneity in nitrate dynamics of the Des Moines River despite the homogeneity of agricultural land cover in the watershed.
Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.
We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0'.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.
This dataset provides multi-year, ground-truth visual observations of surface water flow/no-flow conditions within the East-Taylor Watershed, Colorado, collected during June and July of 2025 and 2026. In June and July 2025, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function Scientific Focus Area (SFA) and Rocky Mountain Biological Laboratory (RMBL) Colorado Headwaters Ecological Spectroscopy Study (CHESS) campaign (further details are provided within the CHESS Project Description). We obtained 377 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. These ground-truth observations were collected to validate classification maps from remote sensing data and model results within the East-Taylor Watershed. In 2025, flow/no-flow measurements were collected using a field-based app for the CHESS Campaign (Zerion iForm). Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included information about visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, and beaver activity. For 2025 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2025_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2025_Water_Observations_Locations.kmz); (3) photos (.jpg and .jpeg) of the water observation points, organized by location, contained within 2025_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2025_Water_Observation_Protocols.pdf). In June and July 2026, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function SFA project. We obtained 365 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. The 2026 observations focused on collecting repeat measurements at the 2025 flow/no-flow observation locations conducted as part of the CHESS campaign. These ground-truth observations were collected to understand differences in flow/no-flow in 2026, given the unprecedented 2026 drought in Colorado. In 2026, flow/no-flow measurements were collected using ArcGIS (Geographic Information System) Survey123. Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included repeat information from the 2025 water observation effort, including visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, beaver activity, and a new metadata component of estimated stream depth (for select locations). For 2026 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2026_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2026_Water_Observations_Locations.kmz); (3) photos (.jpg) of the water observation points, organized by location, contained within 2026_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2026_Water_Observation_Protocols.pdf). For 2025 and 2026 water observations, this dataset contains: (1) a location metadata file (locations.csv); (6) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (7) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. 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. This work was 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. 2026-09-02: This dataset was updated to include 2026 water observation measurements. The 2025 observation files were also updated to ensure a consistent file naming convention across water observation years.