Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial variations”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An intermediate model for fitting triplet–triplet annihilation in phosphorescent organic light emitting diode materials

Triplet–triplet annihilation (TTA) is one of the primary contributors to efficiency roll-off and permanent material degradation in phosphorescent organic light-emitting diodes. Here, the two limiting case models typically used to quantify this quenching mechanism are multi-step Dexter and single-step Förster, which, respectively, assume ideal Fickian diffusion or perfect trapping of triplet excitons. For device-relevant guest doping levels (typically 5–12 vol. %), both significant diffusion of excitons and trapping due to spatial and energetic disorder exist, so neither conventional model fits experimental data well. We develop and validate an intermediate TTA model, which is a weighted average of the limiting cases of pure radiative decay (no TTA) and multi-step Dexter based TTA that returns an effective TTA rate constant and a parameter quantifying the portion of well-isolated excitons. Kinetic Monte-Carlo simulations and time-resolved photoluminescence measurements of an archetype host–guest system demonstrate that our intermediate model provides significantly improved fits with more realistic physical values, is more robust to variations in experimental conditions, and provides an analysis framework for the effects of trapping on TTA.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Global root traits (GRooT) database

Motivation: Trait data are fundamental to the quantitative description of plant form and function. Although root traits capture key dimensions related to plant responses to changing environmental conditions and effects on ecosystem processes, they have rarely been included in large-scale comparative studies and global models. For instance, root traits remain absent from nearly all studies that define the global spectrum of plant form and function. Thus, to overcome conceptual and methodological roadblocks preventing a widespread integration of root trait data into large-scale analyses we created the Global Root Trait (GRooT) Database. GRooT provides readyto- use data by combining the expertise of root ecologists with data mobilization and curation. Specifically, we (a) determined a set of core root traits relevant to the description of plant form and function based on an assessment by experts, (b) maximized species coverage through data standardization within and among traits, and (c) implemented data quality checks. Main types of variables contained: GRooT contains 114,222 trait records on 38 continuous root traits. Spatial location and grain: Global coverage with data from arid, continental, polar, temperate and tropical biomes. Data on root traits were derived from experimental studies and field studies. Time period and grain: Data were recorded between 1911 and 2019. Major taxa and level of measurement: GRooT includes root trait data for which taxonomic information is available. Trait records vary in their taxonomic resolution, with subspecies or varieties being the highest and genera the lowest taxonomic resolution available. It contains information for 184 subspecies or varieties, 6,214 species, 1,967 genera and 254 families. Owing to variation in data sources, trait records in the database include both individual observations and mean values. Software format: GRooT includes two csv files. A GitHub repository contains the csv files and a script in R to query the database.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Microtopographic and Depth Controls on Active Layer Chemistry in Arctic Polygonal Ground: Supporting Data

Polygonal ground is a signature characteristic of the Arctic, and permafrost thaw can potentially generate substantial feedbacks to Arctic ecosystems and climate. This study describes the first comprehensive spatial examination of active layer biogeochemistry that extends across high- and low-centered polygons and their features, including depth. Water chemistry measurements were made on active layer water samples collected near Barrow, Alaska during summer, 2012. Several significant differences in chemistry were observed between high- and low-centered polygons suggesting polygon types may be useful for landscape-scale geochemical classification. However, differences were found for polygon features (centers and troughs) for analytes that were not significant for type, suggesting that finer scale features control biogeochemistry in a different way than polygon type. Depth variations were also significant, demonstrating important multi-dimensional aspects of polygonal ground biogeochemistry. These results have major implications for understanding how polygonal ground ecosystems function, and how they may respond to future change. This data package includes one *.csv data file and one *.pdf user guide. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Seismic monitoring of CO 2 geosequestration using multi-well 4D DAS VSP: Stage 3 of the CO2CRC Otway project

An important part of any CO 2 geosequestration project is to ensure CO 2 containment and conformance in the subsurface. This is generally done by implementing a comprehensive, risk-based Measurement, Monitoring and Verification plan, a key element of which is active time-lapse seismic monitoring. However, high cost and environmental impact of the standard surface seismic monitoring dictate the need for a cost-effective and environmentally friendly alternative. An opportunity to develop such method emerges with advances in distributed acoustic sensing (DAS) technology, which turns an optical fibre into a seismic sensor with dense spatial sampling. DAS can be permanently deployed in multiple wells across the geosequestration site providing a robust and non-intrusive network of seismic receivers. This approach was developed and tested in the CO2CRC Otway project, where injection of 15 kt of CO 2 at 1.5 km depth was monitored with a 4D vertical seismic profiling (VSP) using five borehole DAS arrays and mobile vibroseis sources. The 4D DAS VSP in each of the five wells provides broadly consistent images of the CO 2 plume with some differences due to different illumination of the target horizon, lateral variation of velocities, and seismic anisotropy. When the newly injected CO 2 reaches a CO 2 plume created as a result of an earlier injection into the same formation ~600 m updip, 4D DAS VSP shows a change in reflectivity in that area and beyond. Furthermore, this shows a potential of 4D DAS VSP for monitoring gas injection into gas-saturated reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

Introduction to special issue: Colloidal quantum dots

The quantum confinement effect has driven the development and exploration of colloidal semiconductor nanocrystals over approximately the past three decades. These semiconductor nanostructures are typically synthesized as stable colloidal dispersions in solution-phase chemical batch reactions at relatively low temperatures. Researchers have developed a wide range of reaction conditions to produce a variety of shapes, compositions, and structures with well-controlled sizes and low polydispersities. The shape can be controlled to produce confinement in one, two, or all three spatial dimensions, the latter defining a “quantum dot.” Combining two or more semiconductors within the various nanostructures produces additional optical and electronic degrees of freedom. These mixed compositional nanostructures have properties not accessible in bulk semiconductor systems. Further, researchers have fabricated novel nanoheterostructures, such as “dot-in-rods,” quantum wells, core/shell structures, giant-core/shell tetrapods where the core is one material and the arm is another, Janus structures that consist of spherical structures with well-defined internal heterostructures, nano-dumbbells, hollow-core structures, and many other novel and interesting structures. The structural variations are seemingly only limited by the creativity of the various researchers engaged in these endeavors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Climate and air pollution impacts of generating biopower from forest management residues in California

California faces crisis conditions on its forested landscapes. A century of aggressive logging and fire suppression in combination with conditions exacerbated by climate change have created an ongoing ecological, economic, and public health emergency. Between commercial harvests on California’s working forestlands and the increasing number of acres the state treats each year for fire risk reduction and carbon sequestration, California forests generate millions of tons of woody residues annually—residues that are typically left or burned in the field. State policymakers have turned to biomass electricity generation as a key market for woody biomass in the hope that it can support sustainable forest management activities while also providing low-carbon renewable electricity. However, open questions surrounding the climate and air pollution performance of electricity generation from woody biomass have made it difficult to determine how best to manage the risks and opportunities posed by forest residues. The California Biomass Residue Emissions Characterization (C-BREC) model offers a spatially-explicit life cycle assessment framework to rigorously and transparently establish the climate and air pollution impacts of biopower from forest residues in California under current conditions. The C-BREC model characterizes the variable emissions from different biomass supply chains as well as the counterfactual emissions from prescribed burn, wildfire, and decay avoided by residue mobilization. We find that the life cycle ‘carbon footprint’ of biopower from woody residues generated by recent forest treatments in California ranges widely—from comparable with solar photovoltaic on the low end to comparable with natural gas on the high end. This variation stems largely from the heterogeneity in the fire and decay conditions these residues would encounter if left in the field, with utilization of residue that would otherwise have been burned in place offering the best climate and air quality performance. California’s energy and forest management policies should account for this variation to ensure desired climate benefits are achieved.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Energy Distribution of the Galactic Center Excess’s Sources

The Galactic Center Excess (GCE) may yet herald the discovery of annihilating dark matter. Weighing against that conclusion are analyses showing evidence for dim point sources within the spatial structure of the emission. Because of technical limitations these analyses are purely spatial with all spectral information that could disentangle the excess from astrophysical backgrounds discarded. Here, we demonstrate that a neural network simulation-based inference approach can jointly analyze the spatial and spectra data. The addition is profound: energy information drives the putative point sources to be significantly dimmer, indicating either the GCE is truly diffuse in nature or made of an exceptionally large number of sources. Quantitatively, for our best fit background model, the excess is essentially consistent with Poisson emission as predicted by dark matter. If due to point sources, our median prediction is O(10^{5}) sources, or more than 35 000 at 90% confidence-both orders of magnitude larger than the hundreds preferred by earlier point-source analyses of the GCE, although variations allowed by background systematics could reduce the required number of sources by roughly an order of magnitude.

List, Florian↗
DOE OSTI2023

Multiwavelength Bulge–Disk Decomposition for the Galaxy M81 (NGC 3031). I. Morphology

A panchromatic investigation of morphology for the early-type spiral galaxy M81 is presented in this paper. We perform bulge–disk decomposition in M81 images at a total of 20 wave bands from far-UV to near-IR (NIR) obtained with GALEX, Swift, Sloan Digital Sky Survey, WIYN, Two Micron All Sky Survey, Wide-field Infrared Survey Explorer, and Spitzer. Morphological parameters such as Sérsic index, effective radius, position angle, and axis ratio for the bulge and the disk are thus derived at all of the wave bands, which enables quantifying the morphological K -correction for M81 and makes it possible to reproduce images for the bulge and the disk in the galaxy at any wave band. The morphology as a function of wavelength appears as a variable-slope trend of the Sérsic index and the effective radius, in which the variations are steep at UV–optical and shallow at optical–NIR bands; the position angle and the axis ratio keep invariable at least at optical–NIR bands. It is worth noting that the Sérsic index for the bulge reaches ~4–5 at optical and NIR bands, but drops to ~1 at UV bands. This difference brings forward a caveat that a classical bulge is likely misidentified for a pseudobulge or no bulge at high redshifts where galaxies are observed through rest-frame UV channels with optical telescopes. The next work of this series is planned to study spatially resolved spectral energy distributions for the bulge and the disk, respectively, and thereby explore stellar population properties and star formation/quenching history for the galaxy composed of the subsystems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗
DOE OSTI2021

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Differential holography

Differential Holography technology measures the amplitude and/or phase of, e.g., an incident linearly polarized spatially coherent quasi-monochromatic optical field by optically computing the first derivative of the field and linearly mapping it to an irradiance signal detectable by an image sensor. This information recorded on the image sensor is then recovered by a simple algorithm. In some embodiments, an input field is split into two or more beams to independently compute the horizontal and vertical derivatives (using amplitude gradient filters in orthogonal orientations) for detection on one image sensor in separate regions of interest (ROIs) or on multiple image sensors. A third unfiltered beam recorded in a third ROI directly measures amplitude variations in the input field to numerically remove its contribution as noise before recovering the original wavefront using a numerical in algorithm. When combined, the measured amplitude and phase constitute a holographic recording of the incident optical field.

Ravizza, Frank↗
DOE OSTI2022

3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy

A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.

58 GEOSCIENCES↗
DOE OSTI2025

Wind Turbine Sound Setbacks and Supply Curves: Ordinances and Extrapolated Trends, 110 Hub Height, 130 Rotor Diameter

This dataset provides a comprehensive set of wind turbine sound setbacks from every residential structure in the contiguous United States (CONUS). A sound setback is defined as the minimum required distance between a residential structure and a hypothetical turbine installation site to ensure that modeled sound levels received at the residence do not exceed local sound ordinances, which are commonly expressed in A-weighted decibels (dBA). Therefore, sound setbacks are a local spatial assessment combining multiple factors, including the sound pressure curve as a function of the observer location (distance and direction) relative to the turbine, local sound regulations, and the geographical distribution of residential structures. The dataset is organized into multiple scenario-based products, detailed as follows: 1. Existing and extrapolated sound setbacks. An existing scenario characterizes sound setbacks only in states or counties that have implemented sound regulations as of 2022. The extrapolated scenarios extend a constant sound threshold to counties that lack explicit sound regulations, with thresholds ranging from 35 to 60 dBA, in 5-dBA increments reflecting the variation observed in current sound ordinances. 2. Sound setbacks in directional and worst scenarios. The directional scenario accounts for the distance and orientation of residential structures relative to a hypothetical turbine location, utilizing the turbine's sound emissions in that specific direction. In contrast, the worst scenario takes loudest sound level at each distance step from the turbine, irrespective of directional considerations, which aligns with current industry practice. 3. Supply curves for Open and Reference Access scenarios. This dataset includes supply curves generated by the reV model, which integrates each of the above sound setbacks into both Open and Reference siting scenarios. In addition, two Open and Reference baselines scenarios were included which do not consider sound setbacks for comparative analysis. All sound setback data are stored in TIF files, with partial maps of the data provided in PNG format. The values in the sound setback raster range from 0 to 1, representing the fraction of developable land within a 90 meter by 90 meter pixel due to sound ordinances. A value of 0 indicates areas where wind energy development is prohibited, while a value of 1 signifies areas fully permissible. The wind turbine parameters used in the sound modeling are based on the land-based turbine from International Energy Agency (IEA), featuring a rated electrical power of 3.4 MW, a rotor diameter of 130 meters, and a hub height of 110 meters. The atmospheric conditions, including wind speed/direction, turbulence, air temperature, relative humidity, and air pressure, that drive the sound generation are obtained from the WIND Toolkit dataset.

Array↗
DOE OSTI2025

Angularly resolved spectral reconstruction of x rays via filter pack attenuation

We have designed a new filter pack array to measure angular variations in x-ray spectra during a single shot. The filter pack was composed of repeating identical columns of aluminum and copper filters of varying thicknesses. These columns were located at different positions to measure the spectrum at each corresponding angle. This array was utilized in an experiment to measure the energy evolution of betatron x rays in a laser wakefield accelerator by curving the wakefield with a transverse density gradient, streaking the x rays across the array in front of an x-ray charge-coupled device (CCD) camera. After subtracting the background and “flattening” the image to remove spatial nonuniformities, a critical energy was calculated for each position that produced the best agreement with the measured signal. There was a clear change in critical energy with angle, shedding light on the dynamics of the electrons that traveled through the accelerator. Furthermore, these angles correspond to distinct emission times, covering a timescale of tens of picoseconds. The filter pack was capable of recovering these angular details without the impact of errors introduced by shot-to-shot variability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effects of boundary hydraulics, dissolved oxygen, and dissolved organic carbon on growth and death dynamics of aerobic microbes in riverbed dune-induced hyporheic zones

Surface and groundwater interact in the hyporheic zone beneath and adjacent to rivers in the presence of a diverse microbial community. Heterotrophic bacteria mediate a range of environmentally important reactions, yet few studies have quantified bacterial growth and death dynamics in the hyporheic zone, and none have systematically analyzed their response to variations in hydraulic or chemical conditions. Here we used MODFLOW and SEAM3D to simulate hydraulics; dissolved oxygen (DO) and dissolved organic carbon (DOC) transport; and aerobic microbial metabolism, growth, and death in hyporheic zones induced by riverbed dunes. We ran simulations both with and without growth/death processes, and varied hydraulic parameters and DO/DOC boundary concentrations. Microbial biomass reached steady state (t = 3 days) in every simulation, at which time there was greater biomass and DOC biodegradation rates in the hyporheic flowcell (300% and 85% higher for the base case, respectively) when accounting for microbial growth dynamics. This occurred as microbial biomass tailored its spatial distribution to the availability of DO and DOC, demonstrating the importance of simulating growth/death processes. Biomass generally increased with hyporheic flow cell area as upwelling groundwater decreased. When varying surface water DO and DOC source concentrations relative to the base case, the greatest effect on biomass occurred when increasing DOC and decreasing DO. We determined minimum DO and DOC steady-state concentrations required for microbial growth, but the minimums were not absolute or related by stoichiometry. Increasing DOC created a smaller area of microbes with higher concentrations relative to the base case. Increasing DO slightly increased the area occupied by microbes while keeping the total biomass nearly constant. Overall, microbial growth and death dynamics depend on DO and DOC availability in the hyporheic zone, which is dependent on DOC/DO boundary concentrations and hyporheic flow paths, and in turn the hydraulic interaction between surface water and groundwater.

54 ENVIRONMENTAL SCIENCES↗