Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequential sampling”

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 37 records · Page 2

Classifier-guided sampling for discrete variable, discontinuous design space exploration: Convergence and computational performance

A classifier-guided sampling (CGS) method is introduced for solving engineering design optimization problems with discrete and/or continuous variables and continuous and/or discontinuous responses. The method merges concepts from metamodel-guided sampling and population-based optimization algorithms. The CGS method uses a Bayesian network classifier for predicting the performance of new designs based on a set of known observations or training points. Unlike most metamodeling techniques, however, the classifier assigns a categorical class label to a new design, rather than predicting the resulting response in continuous space, and thereby accommodates nondifferentiable and discontinuous functions of discrete or categorical variables. The CGS method uses these classifiers to guide a population-based sampling process towards combinations of discrete and/or continuous variable values with a high probability of yielding preferred performance. Accordingly, the CGS method is appropriate for discrete/discontinuous design problems that are ill-suited for conventional metamodeling techniques and too computationally expensive to be solved by population-based algorithms alone. In addition, the rates of convergence and computational properties of the CGS method are investigated when applied to a set of discrete variable optimization problems. Results show that the CGS method significantly improves the rate of convergence towards known global optima, on average, when compared to genetic algorithms.

42 ENGINEERING↗

Scheduled Sequential Compressed Spectrum Sensing for Wideband Cognitive Radios

The support for high data rate applications with the cognitive radio technology necessitates wideband spectrum sensing. However, it is costly to apply long-term wideband sensing and is especially difficult in the presence of uncertainty, such as high noise, interference, outliers, and channel fading. In this paper, we propose scheduling of sequential compressed spectrum sensing which jointly exploits compressed sensing (CS) and sequential periodic detection techniques to achieve more accurate and timely wideband sensing. Instead of invoking CS to reconstruct the signal in each period, our proposed scheme performs backward grouped-compressed-data sequential probability ratio test (backward GCD-SPRT) using compressed data samples in sequential detection, while CS recovery is only pursued when needed. This method on one hand significantly reduces the CS recovery overhead, and on the other takes advantage of sequential detection to improve the sensing quality. Furthermore, we propose (a) an in-depth sensing scheme to accelerate sensing decision-making when a change in channel status is suspected, (b) a block-sparse CS reconstruction algorithm to exploit the block sparsity properties of wide spectrum, and (c) a set of schemes to fuse results from the recovered spectrum signals to further improve the overall sensing accuracy. Finally, extensive performance evaluation results show that our proposed schemes can significantly outperform peer schemes under sufficiently low SNR settings.

42 ENGINEERING↗

Automated Coupling of Nanodroplet Sample Preparation with Liquid Chromatography–Mass Spectrometry for High-Throughput Single-Cell Proteomics

Single-cell proteomics can provide critical biological insight into the cellular heterogeneity that is masked by bulk-scale analysis. Here, we have developed a nanoPOTS (nanodroplet processing in one pot for trace samples) platform and demonstrated its broad applicability for single-cell proteomics. However, because of nanoliter-scale sample volumes, the nanoPOTS platform is not compatible with automated LC-MS systems, which significantly limits sample throughput and robustness. To address this challenge, we have developed a nanoPOTS autosampler allowing fully automated sample injection from nanowells to LC-MS systems. We also developed a sample drying, extraction, and loading workflow to enable reproducible and reliable sample injection. The sequential analysis of 20 samples containing 10 ng tryptic peptides demonstrated high reproducibility with correlation coefficients of >0.995 between any two samples. The nanoPOTS autosampler can provide analysis throughput of 9.6, 16, and 24 single cells per day using 120, 60, and 30 min LC gradients, respectively. As a demonstration for single-cell proteomics, the autosampler was first applied to profiling protein expression in single MCF10A cells using a label-free approach. At a throughput of 24 single cells per day, an average of 256 proteins was identified from each cell and the number was increased to 731 when the Match Between Runs algorithm of MaxQuant was used. Using a multiplexed isobaric labeling approach (TMT-11plex), ~77 single cells could be analyzed per day. We analyzed 152 cells from three acute myeloid leukemia cell lines, resulting in a total of 2558 identified proteins with 1465 proteins quantifiable (70% valid values) across the 152 cells. These data showed quantitative single-cell proteomics can cluster cells to distinct groups and reveal functionally distinct differences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applications of Hyporheic Zone Sampling for Assessing VOC Plume Engineered Remediation and Natural Attenuation Effectiveness - 20414

Determining the efficacy of engineered remedial actions and/or natural attenuation as an acceptable EPA Record of Decision (ROD) must address long-term temporal and spatial variability. Since final plume downslope destinations are aquifers and related surface streams and rivers, hyporheic zone sampling is essential. Comprehensive hyporheic sampling documents plume groundwater remediation prior to upwelling into stream benthic strata and natural attenuation within this dynamic zone itself. Our intensive sampling of the Savannah River Site (SRS) Chemical, Metals, and Pesticides (CMP) Pits VOC plume-fringe from 2005 to 2018 documented PCE natural attenuation and effectiveness of In Situ Thermal Treatment (ISTT). Most plume hyporheic upwelling was delimited to a 20 m inflow reach along Pen Branch. Our study utilized corroborative sampling methodologies: ambient hole-water grab samples: passive diffusion bags (PDB); and hole-core sediment samples to confirm station to station spatial and temporal trends. VOC hot spot detections were confirmed from both groundwater and soil core samples. VOC pathways through the Pen Branch hyporheic zone were assessed at hole-installation using sequential soil samples from augered cores and grab-samples of hole-bottom water. Longer term, composite water samples were collected from PDBs deployed for a minimum two-week equilibration period. Soil samples were collected from the 15 and 65 cm core levels in December 2017 and June 2018 at six stations. Additional core depths of 40 and 85 cm were sampled in June 2018 at the highest VOC station (5DB80). All samples were analyzed in an EPA certified lab by purge and trap GC-MS. Effectiveness of previous ISTT and soil vapor extraction (SVE) for this CMP Pits PCE plume, was clearly documented by significant downward trends in VOC at hyporheic zone stations prior to December 2015. Achieving final stage conversions of PCE daughter products to ethene and CO{sub 2} may be limited depending upon suitable interactions between hyporheic zone sediment morphology, hydrology, and redox conditions. Because declines in VOC daughter products at our hyporheic stations may also be caused by physical actions of dilution and dispersion rather than true biochemical dechlorination, full quantification of natural attenuation requires the use of compound-specific isotope ratio analysis (CSIA) and microbial 16S rRNA with qPCR. To assist these next-phase efforts, a natural attenuation index (NAI) was used to help produce quicker, more cost-effective decision making for locating CSIA and microbial sampling efforts. The NAI used natural logarithm transformed GC-MS results for station VOC degradation compounds versus VOC parent compounds. Wide differences were observed between stations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data from a throughfall exclusion experiment: Fine root dynamics, morphology, chemistry, and AMF colonization across four lowland Panamanian forests

Fine roots regulate forest nutrient, carbon, and water cycling, yet their variation within and among tropical forests remains under-characterized. We quantified root productivity, disappearance, and stocks to 1 m using minirhizotron imaging, and we measured morphology, elemental composition [root carbon (C), root nitrogen (N), root phosphorus (P)], and arbuscular mycorrhizal fungi (AMF) colonization to 20 cm using ingrowth cores and sequential coring. Sampling took place in four distinct lowland Panamanian forests (32 plots; 8 per forest) from 2018 through 2022 under control and throughfall-exclusion (drought) treatments in the Panama Rainforest Changes with Experimental Drying (PARCHED) experiment.The dataset is presented as an Excel workbook with six tabs. The first tab is the data dictionary. Tab S1 contains ingrowth-core production and mortality, morphology and soil moisture. Tab S2 contains sequential-coring standing stocks with associated morphology and soil moisture. Tab S3 contains minirhizotron row data records to 1 m depth, including per-frame root length and diameter, normalized length metrics, and session timing. Tab S4 contains AMF colonization. Tab S5 contains fine-root chemistry at 0–10 cm, reporting %P, %C, %N, and C:N for samples collected via ingrowth cores and sequential-coring standing stocks. CSV mirrors for each tab are provided, and a KML file supplies coordinates for all 32 plots.Key variables span live and dead fine-root biomass (and coarse fractions where applicable), specific root length (SRL) and area (SRA), diameter, root tissue density (RTD), soil moisture, AMF colonization, root %N, %C, %P, and C:N, along with minirhizotron root length and diameter. Depth, season, treatment, and plot/site identifiers are included to support cross-tab integration and analysis from 0–100 cm (minirhizotron) and 0–20 cm (cores).Units are reported in-column and missing values are coded as NA. No special software is required to open or use the files (Excel, CSV, and KML compatible).

54 ENVIRONMENTAL SCIENCES↗

At-Line Sampling and Characterization of Pyrolytic Vapors from Biomass Feedstock Blends Using SPME-GC/MS-PCA: Influence of Char on Fast Pyrolysis

Solid-phase microextraction (SPME) coupled with GC-MS analysis was used for at-line sampling of pyrolytic vapors produced during fast pyrolysis of biomass. The pure and binary blends of switchgrass (SWG) and pine harvest residues (PT6) were used as feedstock. Sequential SPME sampling allowed for monitoring of changes in the pyrolysis vapors as char accumulated in the fluid bed. The concentration and composition of the vapors desorbed from the SPME fibers were investigated using GC-MS, and the data sets were then analyzed using principal component analysis (PCA) to compare the composition of the pyrolysis vapors over the course of the pyrolysis run. The chemical composition of both carbohydrate and lignin fragments varied as the char builds up in the reactor bed. Fragments derived from cellulose and xylan included anhydrosugars, furans and light oxygenated compounds. Lignin fragments included methoxyphenols, phenolic ketones and aldehydes, low molecular weight aromatics. The composition of the carbohydrate fragments changed more than the lignin fragments as the char build-up in the fluid bed. This combination of SPME-GC/MS-PCA were a novel, easy and effective method for measuring the composition and changes in the composition of pyrolysis vapors during the fast pyrolysis process. Here, this work also highlighted the effect of char build-up on the composition of the overall pyrolysis vapors.

09 BIOMASS FUELS↗

Persistent Sampling: Enhancing the Efficiency of Sequential Monte Carlo

Sequential Monte Carlo (SMC) samplers are powerful tools for Bayesian inference but suffer from high computational costs due to their reliance on large particle ensembles for accurate estimates. We introduce persistent sampling (PS), an extension of SMC that systematically retains and reuses particles from all prior iterations to construct a growing, weighted ensemble. By leveraging multiple importance sampling and resampling from a mixture of historical distributions, PS mitigates the need for excessively large particle counts, directly addressing key limitations of SMC such as particle impoverishment and mode collapse. Crucially, PS achieves this without additional likelihood evaluations-weights for persistent particles are computed using cached likelihood values. This framework not only yields more accurate posterior approximations but also produces marginal likelihood estimates with significantly lower variance, enhancing reliability in model comparison. Furthermore, the persistent ensemble enables efficient adaptation of transition kernels by leveraging a larger, decorrelated particle pool. Experiments on high-dimensional Gaussian mixtures, hierarchical models, and non-convex targets demonstrate that PS consistently outperforms standard SMC and related variants, including recycled and waste-free SMC, achieving substantial reductions in mean squared error for posterior expectations and evidence estimates, all at reduced computational cost. PS thus establishes itself as a robust, scalable, and efficient alternative for complex Bayesian inference tasks.

Karamanis, Minas↗

Toward a modular multi-material nanoparticle synthesis and assembly strategy via bionanocombinatorics: bifunctional peptides for linking Au and Ag nanomaterials

Materials-binding peptides represent a unique avenue towards controlling the shape and size of nanoparticles (NPs) grown under aqueous conditions. Here, employing a bionanocombinatorics approach, two such materials-binding peptides were linked at either end of a photoswitchable spacer, forming a multi-domain materials-binding molecule to control the in situ synthesis and organization of Ag and Au NPs under ambient conditions. These multi-domain molecules retained the peptides’ ability to nucleate, grow, and stabilize Ag and Au NPs in aqueous media. Disordered co-assemblies of the two nanomaterials were observed by TEM imaging of dried samples after sequential growth of the two metals, and showed a clustering behavior that was not observed without both metals and the linker molecules. While TEM evidence indicated the formation of AuNP/AgNP assemblies upon drying, SAXS analysis indicated that no extended assemblies existed in solution, suggesting that sample drying plays an important role in facilitating NP clustering. Molecular simulations and experimental data revealed tunable materials-binding based upon the isomerization state of the photoswitchable unit and metal employed. This work is a first step in generating externally actuated biomolecules with specific material-binding properties that could be used as the building blocks to achieve multi-material switchable NP assemblies.

36 MATERIALS SCIENCE↗

Evolution of LIBS technology to mobile instrumentation for expediting firearm-related investigations at the laboratory and the crime scene

Gunshot residue (GSR) is one of the few forensic disciplines that lack accurate screening techniques. This study proposes using a mobile LIBS instrument to detect inorganic GSR and compares performance to a previously validated laboratory instrument. The mobile LIBS is designed with advanced configurations specifically for on-site GSR analysis, including a CMOS detector and a sampling chamber that holds up to six typical GSR collection devices with separate gas flow ports to prevent cross-contamination. A significant novelty of the portable instrument is its image magnification, which allows quick searching and visualization of GSR particle morphology. The single-particle imaging and elemental composition capability is one of a kind and offers superior confirmatory features for GSR. The mobile LIBS performance was evaluated for residues collected from the hands of shooters (100 samples) and non-shooters (200 background samples), analyzed sequentially by the mobile instrument and then the laboratory instrument. Accuracies better than 98.8% were obtained by both instruments, demonstrating their suitability for trace IGSR detection from skin specimens. Implementation of this methodology is anticipated to drastically speed up response times (i.e., from several hours per sample by standard SEM-EDS practice to a few minutes by LIBS). The screening methods can be easily incorporated into workflows to improve decision-making processes at the crime scene and laboratory settings, reduce backlogs, and improve case management.

47 OTHER INSTRUMENTATION↗

Hydrogeological Model Selection Among Complex Spatial Priors

Hydrogeological field studies rely often on a single conceptual representation of the subsurface. This is problematic since the impact of a poorly chosen conceptual model on predictions might be significantly larger than the one caused by parameter uncertainty. Furthermore, conceptual models often need to incorporate geological concepts and patterns in order to provide meaningful uncertainty quantification and predictions. Consequently, several geologically realistic conceptual models should ideally be considered and evaluated in terms of their relative merits. Here, we propose a full Bayesian methodology based on Markov chain Monte Carlo to enable model selection among 2-D conceptual models that are sampled using training images and concepts from multiple-point statistics. More precisely, power posteriors for the different conceptual subsurface models are sampled using sequential geostatistical resampling and Graph Cuts. To demonstrate the methodology, we compare and rank five alternative conceptual geological models that have been proposed in the literature to describe aquifer heterogeneity at the MAcroDispersion Experiment site in Mississippi, USA. We consider a small-scale tracer test for which the spatial distribution of hydraulic conductivity impacts multilevel solute concentration data observed along a 2-D transect. The thermodynamic integration and the stepping-stone sampling methods were used to compute the evidence and associated Bayes factors using the computed power posteriors. We find that both methods are compatible with multiple-point statistics-based inversions and provide a consistent ranking of the competing conceptual models considered.

58 GEOSCIENCES↗

Rare earth element behavior during groundwater–seawater mixing along the Kona Coast of Hawaii

Groundwater and seawater samples were collected from nearshore wells and offshore along the Kona Coast of the Big Island of Hawaii to investigate rare earth element (REE) behavior in local subterranean estuaries. Previously we saw that submarine groundwater discharge (SGD) is the predominant flux of terrestrial waters to the coastal ocean along the arid Kona Coast of Hawaii. Groundwater and seawater samples were filtered through 0.45 μm and 0.02 μm pore-size filters to evaluate the importance of colloidal and soluble (i.e., truly dissolved ionic species and/or low molecular weight [LMW] colloids) fractions of the REEs in the local subterranean estuaries. Mixing experiments using groundwater collected immediately down gradient from a wastewater treatment facility (WWTF) proximal to the Kaloko-Hanokohau National Historic Park, and more “pristine” groundwater from a well constructed in a lava tube at Kiholo Bay, were conducted with local seawater to study the effect of solution composition (i.e., pH, salinity) on the concentrations and fractionation behavior of the REEs as groundwater mixes with seawater in Kona Coast subterranean estuaries. The mixed waters were also filtered through 0.45 or 0.02 μm filters to ascertain the behavior of colloidal and soluble fractions of the REEs across the salinity gradient in each mixing experiment. Concentrations of the REEs were statistically identical (two-tailed Student t-test, 95% confidence) between the sequentially filtered sample aliquots, indicating that the REEs occur as dissolved ionic species and/or LMW colloids in Kona Coast groundwaters. The mixing experiments revealed that the REEs are released to solution from suspended particles or colloids when Kona Coast groundwater waters mix with local seawater. The order of release that accompanies increasing pH and salinity follows light REE (LREE) > middle REE (MREE) > heavy REE (HREE). Release of REEs in the mixing experiments is driven by decreases in the free metal ion activity in solution and the concomitant increase in the amount of each REE that occurs in solution as dicarbonato complexes [i.e., Ln(CO 3 ) 2 - ] as pH increases across the salinity gradient. Input-normalized REE patterns of Kona Coast groundwater and coastal seawater are nearly identical and relatively flat compared to North Pacific seawater, indicating that SGD is the chief source of these trace elements to the ocean along the Kona Coast. Additionally, REE concentrations of the coastal seawater are between 10 and 50 times higher than previously reported open-ocean seawater values from the North Pacific, further demonstrating the importance of SGD fluxes of REEs to these coastal waters. Altogether, these observations indicate that large-scale removal of REEs, which characterizes the behavior of REEs in the low salinity reaches of many surface estuaries, is not a feature of the subterranean estuary along the Kona Coast. A large positive gadolinium (Gd) anomaly characterizes groundwater from the vicinity of the WWTF. The positive Gd anomaly can be traced to the coastal ocean, providing further evidence of the impact of SGD on the coastal waters. Estimates of the SGD fluxes of the REEs to the coastal ocean along the Kona Coast (i.e., 1.3–2.6 mmol Nd day -1 ) are similar to recent estimates of SGD fluxes of REEs along Florida’s east coast and to Rhode Island Sound, all of which points to the importance of SGD as significant flux of REEs to the coastal ocean.

58 GEOSCIENCES↗

A method for determining average damage depth of sawn crystalline silicon wafers

The depth of surface damage (or simply, damage) in crystalline silicon wafers, caused by wire sawing of ingots, is determined by performing a series of minority carrier lifetime (MCLT) measurements. Samples are sequentially etched to remove thin layers from each surface and MCLT is measured after each etch step. The thickness-removed (δt) at which the lifetime reaches a peak value corresponds to the damage depth. This technique also allows the damage to be quantified in terms of effective surface recombination velocity (S eff ). To accomplish this, the MCLT data are converted into an S eff vs δt plot, which represents a quantitative distribution of the degree of damage within the surface layer. We describe a wafer preparation procedure to attain reproducible etching and MCLT measurement results. We also describe important characteristics of an etchant used for controllably removing thin layers from the wafer surfaces. Lastly, some typical results showing changes in the MCLT vs δt plots for different cutting parameters are given.

14 SOLAR ENERGY↗

The laser shock station in the dynamic compression sector. I

The Laser Shock Station in the Dynamic Compression Sector [Advanced Photon Source (APS), Argonne National Laboratory] links a laser-driven shock compression platform with high energy x-ray pulses from the APS to achieve in-situ, time-resolved x-ray measurements (diffraction and imaging) in materials subjected to well-characterized, high stress, short duration planar shock waves. Additionally, the combination of this station and the other DCS experimental stations provides a unique and versatile experimental facility to study condensed state phenomena by undertaking real-time, multi-scale measurements in materials subjected to a wide range of shock amplitudes (to above ~350 GPa) and time-durations (~10 ns to 1 μs). The Laser Shock Station uses a 100 joule, 5-17 ns, 351 nm frequency tripled Nd:glass laser with programmable pulse shaping and focal profile smoothing for maximum precision. The laser is capable of a shot every 30 minutes. The interaction chamber contains multiple ports for simultaneous diagnostics, a sample holder that allows sequential exposure for 14 samples without breaking vacuum, the capability to vary the angle between the x-ray beam and laser beam by 135°, and the ability to translate to select one of two types of x-ray beams. X-ray data can be taken with a temporal resolution of ~90 psec. Highly reproducible, planar shocks can be propagated in samples to undertake precise and well-characterized experiments. In a series of ten shots, the absolute variation in the time the shock breaks out at the back of the sample is less than 500ps and the variation in peak particle velocity at the sample/window interface is 4.3%. Lastly, this paper describes the entire DCS Laser Shock Station, including sample fabrication and diagnostics, as well as representative experimental results from shock compressed tantalum to demonstrate the excellent quality of x-ray diffraction data in solid tantalum and to establish the stress for tantalum melting.

47 OTHER INSTRUMENTATION↗

Using Minority Carrier Lifetime Measurement to Determine Saw Damage Characteristics on Si Wafer Surfaces

The damage on the Si wafer surfaces, caused by ingot cutting, is determined from measurement of minority carrier lifetime (..tau..eff). Samples are sequentially etched to remove thin layers from each surface and lifetime is measured after each etch step. The thickness-removed at which the lifetime reaches a peak value corresponds to the damage depth. This technique also allows the depth distribution of the damage to be quantified in terms of surface recombination velocity (SRV). An accurate measurement of ..tau..eff requires corrections to optical reflection, and transmission to account for changes in the surface morphology and in the wafer thickness.

14 SOLAR ENERGY↗

CHESS 2025: Waveform LiDAR data from NEON AOP surveys

This dataset provides Level 1 (L1) full-waveform light detection and ranging (LiDAR) data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. Waveform LiDAR data can provide more detailed information about objects on the ground than discrete point clouds typically do, and they are often used for granular target segmentation and characterization of subcanopy vegetation. 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. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary waveform LiDAR data delivered by NEON and are provided per flightline in compressed Pulsewaves format, an open-source binary file standard. A Pulsewaves object comprises a two files: a pulse (.pls) file, which stores the geographic origin, outgoing vector, and metadata for every laser pulse emitted by the scanner, and a wave file (.wvs), which stores the sequential amplitude samples of the outgoing pulse and the returning signals. The files are published here in their compressed forms (.plz, .wvz). All waveform data were processed following the theoretical workflow described in the NEON L0-to-L1 Waveform LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022a); however, the Pulsewaves output format differs from a legacy format described in that document. Waveform amplitude samples are recorded at 1 nanosecond intervals. All coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Waveform data for the UPTA survey area were collected without incident and the published records are complete. However, both the ALMO and CRBU collections experienced issues that resulted in incomplete data for those areas. On collection day 2018-06-16 a hardware failure caused the waveform digitizer to lose data from the eastern edge of the ALMO site (Figure 22). The waveform data for flightlines 2–20 could not be extracted from the digitizer, and the data proved unrecoverable. As a result, a portion of the site does not have coverage with waveform data. Although no hardware failure was observed during collection over the CRBU area, final waveform files generated by vendor software contained only ~25% of the expected number of return pulses. After discovery, NEON initiated troubleshooting with the vendor. The root cause of the data ablation had not been identified at the time of publication. Additional data will be published in an update to this package if further recovery proves successful. 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↗

Geochemical characterization of lithium deposition in fossil energy wastes

Lithium is a critical mineral used in rechargeable batteries for electric vehicles (EVs) and future modernization of electric grids with sensitive supply chains that are subject to volatility. Therefore, methods to recover lithium from domestic unconventional sources are actively being pursued in the United States. Herein, we report on the lithium recovery potential from four fossil energy waste feedstocks: oil and gas drill cuttings (sample number n = 16); oil and gas produced waters (n > 200) from Marcellus, Bakken and Permian basins; coal byproducts (n>20), such as coal ash; and acid mine drainage treatment solids (AMD solids). Select solid samples underwent sequential extraction to explore lithium hosting phases and simple organic and inorganic acid extraction to explore lithium recovery potential. Our preliminary results found lithium concentrations up to 440 mg/kg in select AMD treatment solids and up to 300 mg/L in Marcellus Shale produced waters, demonstrating these sources have lithium yield potentials comparable to conventionally mined lithium ores and brines. Sequential extractions revealed that lithium is mostly associated with clay and silicate phases in shale drill cuttings, whereas lithium resides in Fe,Mn-oxide reducible phases in AMD solids. Further, the lithium content in produced waters has a strong linear relationship with total dissolved solid (TDS) levels in three separate oil and gas basins. At the same TDS level, Marcellus Shale produced waters contain more lithium compared to Bakken Shale and Permian Basin waters, with higher percentages of Ca and Mg, major cations that might impact lithium recovery efficiencies. Our study demonstrates the heterogeneity of lithium hosts from different fossil energy wastes. Characterization results will inform future lithium recovery from these different feedstocks.

Stuckman, Mengling↗

Tractable minor-free generalization of planar zero-field Ising models

In this work, we present a new family of zero-field Ising models over N binary variables/spins obtained by consecutive 'gluing' of planar and O(1)-sized components and subsets of at most three vertices into a tree. The polynomial time algorithm of the dynamic programming type for solving exact inference (computing partition function) and exact sampling (generating i.i.d. samples) consists of sequential application of an efficient (for planar) or brute-force (for O(1)-sized) inference and sampling to the components as a black box. To illustrate the utility of the new family of tractable graphical models, we first build a polynomial algorithm for inference and sampling of zero-field Ising models over K 33 -minor-free topologies and over K 5 -minor-free topologies—both of which are extensions of the planar zero-field Ising models—which are neither genus- nor treewidth-bounded. Second, we empirically demonstrate an improvement in the approximation quality of the NP-hard problem of inference over the square-grid Ising model in a node-dependent nonzero 'magnetic' field.

97 MATHEMATICS AND COMPUTING↗

Anomaly Detection in Dynamic Networks

Anomaly detection in dynamic communication networks has many important security applications. These networks can be extremely large and so detecting any changes in their structure can be computationally challenging; hence, computationally fast, parallelisable methods for monitoring the network are paramount. For this reason the methods presented here use independent node and edge based models to detect locally anomalous substructures within communication networks. As a first stage, the aim is to detect changes in the data streams arising from node or edge communications. Throughout the thesis simple, conjugate Bayesian models for counting processes are used to model these data streams. A second stage of analysis can then be performed on a much reduced subset of the network comprising nodes and edges which have been identified as potentially anomalous in the first stage. The first method assumes communications in a network arise from an inhomogeneous Poisson process with piecewise constant intensity. Anomaly detection is then treated as a changepoint problem on the intensities. The changepoint model is extended to incorporate seasonal behavior inherent in communication networks. This seasonal behavior is also viewed as a changepoint problem acting on a piecewise constant Poisson process. In a static time frame, inference is made on this extended model via a Gibbs sampling strategy. In a sequential time frame, where the data arrive as a stream, a novel, fast Sequential Monte Carlo (SMC) algorithm is introduced to sample from the sequence of posterior distributions of the change points over time. A second method is considered for monitoring communications in a large scale computer network. The usage patterns in these types of networks are very bursty in nature and don’t fit a Poisson process model. For tractable inference, discrete time models are considered, where the data are aggregated into discrete time periods and probability models are fitted to the communication counts. In a sequential analysis, anomalous behavior is then identified from outlying behavior with respect to the fitted predictive probability models. Seasonality is again incorporated into the model and is treated as a changepoint model on the transition probabilities of a discrete time Markov process. Second stage analytics are then developed which combine anomalous edges to identify anomalous substructures in the network.

97 MATHEMATICS AND COMPUTING↗