Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling variability”

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

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

At least 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Isotopic determination of uranium in soil by laser induced breakdown spectroscopy

Laser-induced breakdown spectroscopy (LIBS) operated under ambient pressure has been evaluated for isotopic analysis of uranium in real-world samples such as soil, with U concentrations in the single digit percentage levels. The study addresses the requirements for spectral decomposition of 235 U and 238 U atomic emission peaks that are only partially resolved. Although non-linear least-square fitting algorithms are typically able to locate the optimal combination of fitting parameters that best describes the experimental spectrum even when all fitting parameters are treated as free independent variables, the analytical results of such an unconstrained free-parameter approach are ambiguous. In this work, five spectral decomposition algorithms were examined, with different known physical properties (e.g., isotopic splitting, hyperfine structure) of the spectral lines sequentially incorporated into the candidate algorithms as constraints. It was found that incorporation of such spectral-line constraints into the decomposition algorithm is essential for the best isotopic analysis. The isotopic abundance of 235 U was determined from a simple two-component Lorentzian fit on the U II 424.437 nm spectral profile. For six replicate measurements, each with only fifteen laser shots, on a soil sample with U concentration at 1.1% w/w, the determined 235 U isotopic abundance was (64.6 ± 4.8)%, and agreed well with the certified value of 64.4%. Another studied U line - U I 682.691 nm possesses hyperfine structure that is comparatively broad and at a significant fraction as the isotopic shift. Thus, 235 U isotopic analysis with this U I line was performed with spectral decomposition involving individual hyperfine components. For the soil sample with 1.1% w/w U, the determined 235 U isotopic abundance was (60.9 ± 2.0)%, which exhibited a relative bias about 6% from the certified value. The bias was attributed to the spectral resolution of our measurement system - the measured line width for this U I line was larger than its isotopic splitting. In conclusion, although not the best emission line for isotopic analysis, this U I emission line is sensitive for element analysis with a detection limit of 500 ppm U in the soil matrix; the detection limit for the U II 424.437 nm line was 2000 ppm.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

A conceptual geochemical model of the geothermal system at Surprise Valley, CA

Characterizing the geothermal system at Surprise Valley (SV), northeastern California, is important for determining the sustainability of the energy resource, and mitigating hazards associated with hydrothermal eruptions that last occurred in 1951. Previous geochemical studies of the area attempted to reconcile different hot spring compositions on the western and eastern sides of the valley using scenarios of dilution, equilibration at low temperatures, surface evaporation, and differences in rock type along flow paths. These models were primarily supported using classical geothermometry methods, and generally assumed that fluids in the Lake City mud volcano area on the western side of the valley best reflect the composition of a deep geothermal fluid. In this contribution, we address controls on hot spring compositions using a different suite of geochemical tools, including optimized multicomponent geochemistry (GeoT) models, hot spring fluid major and trace element measurements, mineralogical observations, and stable isotope measurements of hot spring fluids and precipitated carbonates. We synthesize the results into a conceptual geochemical model of the Surprise Valley geothermal system, and show that high-temperature (quartz, Na/K, Na/K/Ca) classical geothermometers fail to predict maximum subsurface temperatures because fluids re-equilibrated at progressively lower temperatures during outflow, including in the Lake City area. We propose a model where hot spring fluids originate as a mixture between a deep thermal brine and modern meteoric fluids, with a seasonally variable mixing ratio. The deep brine has deuterium values at least 3 to 4‰ lighter than any known groundwater or high-elevation snow previously measured in and adjacent to SV, suggesting it was recharged during the Pleistocene when meteoric fluids had lower deuterium values. The deuterium values and compositional characteristics of the deep brine have only been identified in thermal springs and groundwater samples collected in proximity to structures that transmit thermal fluids, suggesting the brine may be thermal in nature. On the western side of the valley at the Lake City mud volcano, the deep brine-meteoric water mixture subsequently boils in the shallow subsurface, precipitates calcite, and re-equilibrates at about 130 °C. On the eastern side of the valley, meteoric fluid mixes to a greater extent with the deep brine, cools conductively without boiling, and the composition is modified as dissolved elements are sequestered by secondary minerals that form along the cooling and outflow path at temperatures <130 °C. Re-equilibration of geothermal fluids at lower temperatures during outflow explains why subsurface temperature estimates based on classical geothermometry methods are highly variable, and fail to agree with temperature estimates based on dissolved sulfate-oxygen isotopes and results of classical and multicomponent geothermometry applied to reconstructed deep well fluids. Furthermore, the proposed model is compatible with the idea suggested by others that thermal fluids on the western and eastern side of the valley have a common source, and supports the hypothesis that low temperature re-equilibration during west to east flow is the major control on hot spring fluid compositions, rather than dilution, evaporation, or differences in rock type.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Atmospheric CO Concentrations from the CSIRO GASLAB Flask Sampling Network

Individual measurements have been obtained from flask air samples returned to the CSIRO GASLAB. Typical sample storage times range from days to weeks for some sites (e.g. Cape Grim) to as much as 1 year for Macquarie Island and the Antarctic sites. Experiments carried out to test for changes in sample CO mixing ratio during storage have shown significant drifts in some flask types over test periods of several months to years (Cooper et al., 1999). Corrections derived from the test results are applied to network data according to flask type. (Data from the S160 flasks have been rejected due to large and variable drift.)An annual cycle of CO is evident, largely due to an increase in its destruction by the OH radical during the summer months. Additional influences include spatial and seasonal differences in source strength associated with varying trajectories of arriving air at different times of the year. Annual average CO mixing ratios at Alert in 2001 had decreased by about 13 percent of their 1993 value. A string of abnormally high monthly mixing ratios during the winter of 1998-1999, coupled with higher than normal mixing ratios during the previous summer, led to an anomalously high annual average for 1998. Elevated CO mixing ratios were observed on a global scale during 1998 (see Langenfelds et al., 2002).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of chemometric models to classify solid-state U materials by micro-Raman spectroscopy

Discerning uranium (U) particles found in environmental sampling is of interest for monitoring the peaceful use of nuclear material. In this study, a soft independent modeling of class analogy (SIMCA) library was successfully developed for the classification of a four-class system consisting of α-U 3 O 8 , UO 2 , UO 2 (NO 3 ) 2 ·6H 2 O (UNH), and UO 2 O 2 ·4H 2 O (studtite) by Raman spectroscopy in the presence of matrix particulates and additional outliers. Spectral variability between numerous particles of each type revealed appreciable differences as a function of particle size with respect to hydration state and potential oxide phase within each class. Interclass variability was accounted for using both unsupervised and supervised chemometric models. The supervised SIMCA model displayed reasonable sensitivity for each U class and a high degree of specificity by returning whether a spectrum belonged to one class or not. This work demonstrates how Raman spectral features and chemometrics can be used to distinguish U materials from one another and from matrix materials such as flint clay. Combining the outlined chemometric approach with Raman mapping sequences could provide a rapid, nondestructive technique to characterize the chemical composition of a diverse collection of U compounds amid background samples for environmental sampling, nuclear forensics, and industrial applications.

Actinide↗
DOE OSTI2024

A fast and sensitive size-exclusion chromatography method for plasma extracellular vesicle proteomic analysis

Extracellular vesicles (EVs) carry diverse biomolecules derived from their parental cells, making their components excellent biomarker candidates. However, purifying EVs is a major hurdle in biomarker discovery since current methods require large amounts of samples, are time-consuming and typically have poor reproducibility. Here we describe a simple, fast, and sensitive EV fractionation method using size exclusion chromatography (SEC) on a fast protein liquid chromatography (FPLC) system. Our method uses a Superose 6 Increase 5/150, which has a bed volume of 2.9 mL. The FPLC system and small column size enable reproducible separation of only 50 µL of human plasma in 15 min. To demonstrate the utility of our method, we used longitudinal samples from a group of individuals who underwent intense exercise. A total of 838 proteins were identified, of which, 261 were previously characterized as EV proteins, including classical markers, such as cluster of differentiation (CD)9 and CD81. Quantitative analysis showed low technical variability with correlation coefficients greater than 0.9 between replicates. The analysis captured differences in relevant EV proteins involved in response to physical activity. Our method enables fast and sensitive fractionation of plasma EVs with low variability, which will facilitate biomarker studies in large clinical cohorts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical fracture behavior of doped UO 2 using a ball-on-ring equibiaxial flexure test method

Metal oxide dopants, such as titanium and chromium oxides, have garnered considerable attention for their potential to increase grain size (≥ 30 µm) in UO 2 fuel, purportedly enhancing fission gas retention during reactor operation. Fuel performance is significantly impacted by fuel fracture behavior, so it is important to understand the effects of enhanced grain size and dopant content on UO 2 fuel fracture. UO 2 pellets were doped with 0.1 wt% TiO 2 and 0.3 wt% Cr 2 O 3 to alter density and grain size. Inductively coupled plasma mass spectroscopy measured dopant levels pre- and post-sintering. X-ray diffraction revealed lattice changes and microstrain via Rietveld refinement. Field emission scanning electron microscopy determined grain sizes of approximately 30 µm for TiO 2 doping and 7 µm for Cr 2 O 3 doping. Transverse rupture strength tests were performed on over 30 samples per dataset to obtain characteristic strength and Weibull modulus. Results indicate no statistical difference in fracture strength between 0.1 wt% TiO 2 doped UO 2 and undoped UO 2 , while 0.3 wt% Cr 2 O 3 doped UO 2 exhibited a 20% decrease in fracture strength. Doped UO 2 samples also showed reduced Weibull modulus compared to undoped UO 2 , suggesting increased scatter in fracture strength. This study's findings suggest that titanium and chromium oxide doping in UO 2 , regardless of grain size, induce residual stresses, decreasing fracture strength and increasing variability in fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Advanced Method Optimization with Categorical and Constrained Continuous Parameters

Traditional approaches to analytical method optimization (e.g., univariate and “guess-and-check”) can be time-consuming, costly, and often fail to identify true optima within the parameter space. Previous work defined and implemented a generalized technique for method optimization for continuous method parameters, but a knowledge gap remains for the incorporation of categorical variables into these advanced method optimization schemes. This work presents and validates a generalized optimization approach that incorporates both continuous and categorical variables while also utilizing a multivariate, multiobjective optimization scheme with Karush–Kuhn–Tucker conditions to bound the optimization space to solutions within the physical limitations of the parameter space. Method optimization from a case study using GC–MS for the analysis of 11 analytical standards with objectives to minimize peak width and maximize peak height resulted in a 3 orders of magnitude improvement in the average peak height and a 2 orders of magnitude improvement in the average peak width compared to the least optimal (but reasonable) instrumental parameters utilized in this study. This approach to optimization allows for a customizable method optimization in which users can include both continuous and categorical variables to achieve objectives specific to their analytical goals. This approach significantly reduces the labor and cost associated with traditional method development approaches and can be applied in a variety of scientific fields across a range of laboratory techniques (e.g., instrument method development, sample preparation, and extraction techniques).

Amorphous materials↗
DOE OSTI2022

In-Use Passenger Vessel Emission Rates of Black Carbon and Nitrogen Oxides

This study quantified emission factors of black carbon (BC) and nitrogen oxides (NO x ) from 21 engines on in-use excursion vessels and ferries operating in California's San Francisco Bay, including EPA uncertified and Tier 1-4 engines and across engine operating modes. On average, ~60 fuel-based emission factors per engine were measured using a novel combination of exhaust plume capture combined with GPS location and speed data that can be more readily deployed than common portable emissions measurement systems. Here, BC and NO x emission factors (g kg -1 ) were lowest and least variable during fast cruising and highest during maneuvering and docked operation. Selective catalytic reduction (SCR) reduced NO x emissions by ~80% when functional. However, elevated NO x emissions that exceeded corresponding exhaust standards were measured on most Tier 3 and Tier 4 engines sampled, which can be attributed to inactive SCR during frequent low engine load operation. In contrast, BC emissions exceeded the PM emission standard for only one engine, and SCR systems employed as a NO x reduction technology also reduced emitted BC. Using these measured emission factors to compare commuting options, we show that the CO 2 -equivalent emissions per passenger-kilometer are comparable when commuting by car and ferry, but BC and NO x emissions can be several to more than ten times larger when commuting by ferry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Heterogeneity in Hyporheic Flow, Pore Water Chemistry, and Microbial Community Composition in an Alpine Streambed

Abstract The hyporheic zone, where surface water and groundwater mix, is an important microbial habitat where biogeochemical reactions influence water quality. We show that spatial variability in hyporheic flow in the East River near Crested Butte, CO, drives heterogeneity in streambed geochemical conditions and microbial community assemblages, but the diversity of microbial assemblages remains nearly constant throughout the reach. In July 2018, we collected approximately 100 pore water samples at 20‐cm depth and analyzed them for anions, cations, dissolved organic carbon, dissolved organic matter (DOM) quality, and basic water quality parameters. Vertical hydraulic head gradients were also measured to assess the potential for upward or downward flow, and heat tracing was used to quantify vertical flux rates at a subset of locations. We found that regions of the streambed that are more groundwater‐dominated contain less dissolved oxygen, higher concentrations of reduced metals, and more microbially processed, recalcitrant DOM, while more surface water‐dominated locations contain higher dissolved oxygen concentrations and terrestrially derived, labile DOM. 16S rRNA gene sequencing of extracted DNA revealed that microbial community composition varies with geochemical gradients related to hyporheic flow. These findings provide a better understanding of hyporheic controls on streambed biogeochemistry during the baseflow season, which is expected to lengthen with climate change in alpine watersheds due to earlier snowmelt onset and reduced snowpack.

Nelson, Amelia R.↗
DOE OSTI2015

A NuSTAR SURVEY OF NEARBY ULTRALUMINOUS INFRARED GALAXIES

We present a Nuclear Spectroscopic Telescope Array (NuSTAR), Chandra, and XMM-Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120–5453) has a spectrum consistent with a Compton-thick active galactic nucleus (AGN), but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189–2524 and Mrk 273, altering the classification of these borderline sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high-energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2–10 keV to bolometric luminosity, and unabsorbed 2–10 keV to mid-IR [O iv] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Studying Aerosol, Clouds, and Air Quality in the Coastal Urban Environment of Southeastern Texas

A multi-agency succession of field campaigns was conducted in southeastern Texas during July 2021 through October 2022 to study the complex interactions of aerosols, clouds and air pollution in the coastal urban environment. As part of the Tracking Aerosol Convection interactions Experiment (TRACER), the TRACER- Air Quality (TAQ) campaign the Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) and the Convective Cloud Urban Boundary Layer Experiment (CUBE), a combination of ground-based supersites and mobile laboratories, shipborne measurements and aircraft-based instrumentation were deployed. These diverse platforms collected high-resolution data to characterize the aerosol microphysics and chemistry, cloud and precipitation micro- and macro-physical properties, environmental thermodynamics and air quality-relevant constituents that are being used in follow-on analysis and modeling activities. We present the overall deployment setups, a summary of the campaign conditions and a sampling of early research results related to: (a) aerosol precursors in the urban environment, (b) influences of local meteorology on air pollution, (c) detailed observations of the sea breeze circulation, (d) retrieved supersaturation in convective updrafts, (e) characterizing the convective updraft lifecycle, (f) variability in lightning characteristics of convective storms and (g) urban influences on surface energy fluxes. The work concludes with discussion of future research activities highlighted by the TRACER model-intercomparison project to explore the representation of aerosol-convective interactions in high-resolution simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Application of a convolutional neural network in permeability prediction: A case study in the Jacksonburg-Stringtown oil field, West Virginia, USA

Permeability is a critical parameter for understanding subsurface fluid flow behavior, managing reservoirs, enhancing hydrocarbon recovery, and sequestering carbon dioxide. In general, permeability is measured in the laboratory based on subsurface core samples, calculated from well logs or estimated from well tests. However, laboratory measurements and well tests are expensive, time-consuming, and usually limited to a few core samples or wells in a hydrocarbon field or carbon storage site. Machine-learning techniques are good options for generating a rapid, robust, and cost-effective permeability prediction model because of their strengths to recognize the potential interrelationships between input and output variables. Convolutional neural networks (CNN), as a good pattern recognition algorithm, are widely used in image processing, natural language processing, and speech recognition, but are rarely used with regression problems and even less often in reservoir characterization. We have developed a CNN regression model to estimate the permeability in the Jacksonburg-Stringtown oil field, West Virginia, which is a potential carbon storage site and enhanced oil recovery operations field. We also evaluate the concept of the geologic feature image, which is converted from geophysical well logs. Five variables, including two commonly available conventional well logs (the gamma rays [GRs] and bulk density) and three well-log-derived variables (the slopes of the GR and bulk density curves, and shale content), are used to generate a geologic feature image. The CNN treats the geologic feature image as the input and the permeability as the desired output. In addition, the permeability predicted using traditional backpropagation artificial neural networks, which are optimized by genetic algorithms and particle swarm optimization, is compared with the permeability estimated using our CNN. Our results indicate that the CNN regression model provides more accurate permeability predictions than the traditional neural network.

Geochemistry & Geophysics↗
DOE OSTI2014

Bridging the Divide: Linking Genomics to Ecosystem Responses to Climate Change: Final Report

Over the project period, we have addressed the following objectives: 1) assess the effects of altered precipitation patterns (i.e., increased variability in growing season precipitation) on genetic diversity of the dominant C4 grass species, Andropogon gerardii, and 2) experimentally assess the impacts of extreme climatic events (heat wave, drought) on responses of the dominant C4 grasses, A. gerardii and Sorghastrum nutans, and the consequences of these response for community and ecosystem structure and function. Below is a summary of how we have addressed these objectives. Objective 1 After ten years of altered precipitation, we found the number of genotypes of A. gerardii was significantly reduced compared to the ambient precipitation treatments (Avolio et al., 2013a). Although genotype number was reduced, the remaining genotypes were less related to one another indicating that the altered precipitation treatment was selecting for increasingly dissimilar genomes (based on mean pairwise Dice distance among individuals). For the four key genotypes that displayed differential abundances depending on the precipitation treatment (G1, G4, and G11 in the altered plots and G2 in the ambient plots), we identified phenotypic differences in the field that could account for ecological sorting (Avolio & Smith, 2013a). The three altered rainfall genotypes also have very different phenotypic traits in the greenhouse in response to different soil moisture availabilities (Avolio and Smith, 2013c). Two of the genotypes that increased in abundance in the altered precipitation plots had greater allocation to root biomass (G4 and G11), while G1 allocated more biomass aboveground. These phenotypic differences among genotypes suggests that changes in genotypic structure between the altered and the ambient treatments has likely occurred via niche differentiation, driven by changes in soil moisture dynamics (reduced mean, increased variability and changes in the depth distribution of soil moisture) under a more variable precipitation regime, rather than reduced population numbers (A. gerardii tiller densities did not differ between altered and ambient treatments; p = 0.505) or a priori differences in genotype richness (Avolio et al.2013a). This ecological sorting of genotypes, which accounts for 40% of all sampled individuals in the altered plots, is an important legacy of the press chronic climate changes in the RaMPs experiment. Objective 2 In May 2010, we established the Climate Extremes Experiment at the Konza Prairie Biological Station. For the experiment, a gradient of temperatures, ranging from ambient to extreme, were imposed in 2010 and 2011 as a mid-season heat wave under well-watered or severe drought conditions. This study allowed us for the first time to examine species-specific thresholds of responses to climate extremes and assess how these phenotypic responses may impact selection of particular genotypes, with the ultimate goal of linking alterations in individual performance and genetic diversity to ecosystem structure and functioning. We found that tallgrass prairie was resistant to heat waves, but it was not resistant to extreme drought, which reduced aboveground net primary productivity (ANPP) below the lowest level measured in this grassland in almost thirty years (Hoover et al. in press(a)). This extreme reduction in ecosystem function was a consequence of reduced productivity of both C4 grasses and C3 forbs. This reduction in biomass of the C4 grasses (Andropogon gerardii and Sorghastrum nutans) was, in part, due to significant reductions in photosynthesis, leaf water potential and productivity with drought in the dominant grasses species, with S. nutans was more sensitive than A. gerardii to drought (Hoover et al. in press(b)). However, the dominant forb was negatively impacted by the drought more than the dominant grasses, and this led to a reordering of species abundances within the plant community. Although this change in community composition persisted post-drought, ANPP recovered completely the year after drought due to rapid demographic responses by the dominant grass, compensating for loss of the dominant forb. Overall, our results show that an extreme reduction in ecosystem function attributable to a climate extreme (e.g., low resistance) does not preclude rapid ecosystem recovery. Given that dominance by a few species is characteristic of most ecosystems, knowledge of the traits of these species and their responses to climate extremes will be key for predicting future ecosystem dynamics and function. In addition, our research suggests that water stress will dominate photosynthetic and productivity responses caused by discrete drought and heat wave events, rather than direct or additive effects of heat stress, with differential sensitivity in these grasses altering future ecosystem function.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Novel Geometric Operations for Linear Programming

This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Entropy and its Relationship with Statistics

The purpose of our report is to discuss the notion of entropy and its relationship with statistics. Our goal is to provide a manner in which you can think about entropy, its central role within information theory and relationship with statistics. We review various relationships between information theory and statistics—nearly all are well-known but unfortunately are often not recognized. Entropy quantities the "average amount of surprise" in a random variable and lies at the heart of information theory, which studies the transmission, processing, extraction, and utilization of information. For us, data is information. What is the distinction between information theory and statistics? Information theorists work with probability distributions. Instead, statisticians work with samples. In so many words, information theory using samples is the practice of statistics.

97 MATHEMATICS AND COMPUTING↗