Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ecosystem functioning”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Forest carbon sequestration on the west coast, USA: Role of species, productivity, and stockability

Forest ecosystems store large amounts of carbon and can be important sources, or sinks, of the atmospheric carbon dioxide that is contributing to global warming. Understanding the carbon storage potential of different forests and their response to management and disturbance events are fundamental to developing policies and scenarios to partially offset greenhouse gas emissions. Projections of live tree carbon accumulation are handled differently in different models, with inconsistent results. We developed growth-and-yield style models to predict stand-level live tree carbon density as a function of stand age in all vegetation types of the coastal Pacific region, US (California, Oregon, and Washington), from 7,523 national forest inventory plots. We incorporated site productivity and stockability within the Chapman-Richards equation and tested whether intensively managed private forests behaved differently from less managed public forests. We found that the best models incorporated stockability in the equation term controlling stand carrying capacity, and site productivity in the equation terms controlling the growth rate and shape of the curve. RMSEs ranged from 10 to 137 Mg C/ha for different vegetation types. There was not a significant effect of ownership over the standard industrial rotation length (~50 yrs) for the productive Douglas-fir/western hemlock zone, indicating that differences in stockability and productivity captured much of the variation attributed to management intensity. Our models suggest that doubling the rotation length on these intensively managed lands from 35 to 70 years would result in 2.35 times more live tree carbon stored on the landscape. These findings are at odds with some studies that have projected higher carbon densities with stand age for the same vegetation types, and have not found an increase in yields (on an annual basis) with longer rotations. We suspect that differences are primarily due to the application of yield curves developed from fully-stocked, undisturbed, single-species, “normal” stands without accounting for the substantial proportion of forests that don’t meet those assumptions. The carbon accumulation curves developed here can be applied directly in growth-and-yield style projection models, and used to validate the predictions of ecophysiological, cohort, or single-tree style models being used to project carbon futures for forests in the region. Our approach may prove useful for developing robust models in other forest types.

Chisholm, Paul J. (ORCID:0000000238784707)↗
DOE OSTI2022

Microbial extracellular enzyme activity with simulated climate change

It is critical to understand the consequences of environmental change for the microbial regulation of carbon and nutrient cycling. Specifically, understanding microbial community traits, such as extracellular enzyme activity, can help inform nutrient cycling models and address knowledge gaps. We analyzed data on extracellular enzyme activities and litter decomposition from an 18-month experiment in which microbial communities were reciprocally transplanted along a climate gradient in Southern California. Communities were from desert, scrubland, grassland, pine–oak, and subalpine ecosystems. We aimed to test how enzyme activities responded to climate change following transplantation and how those responses related to decomposition rates. We hypothesized that microbial communities would specialize on their native climate conditions, resulting in higher enzyme activities when transplanted back into their native climate. We investigated the relationship between extracellular enzyme Vmax values, substrate mass loss, and microbial biomass as well as variation in these variables across the climate gradient. We found little evidence for climate specialization, and there was rarely a reduction in enzyme functioning after microbial communities were transplanted into new climate conditions. Moreover, observed differences in decomposition were not related to changes in extracellular enzyme potential, although there were significant differences in enzyme activities and decomposition rates across sites. These results suggest that direct, physiological impacts of climate are likely to be important for enzyme-mediated decomposition, but climate specialization will not constrain the microbial response to climate change in our system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Disturbance-accelerated succession increases the production of a temperate forest

Many secondary deciduous forests of eastern North America are approaching a transition in which maturing early successional tree species are declining, resulting in an uncertain future for these forests’ century-long carbon (C) sink. We initiated the Forest Accelerated Succession Experiment (FASET) at the University of Michigan Biological Station to examine the patterns and mechanisms underlying forest C cycling following the stem-girdling induced mortality of >6,700 early successional Populus spp. (aspen) and Betula papyrifera (paper birch). Since 2008, meteorological flux tower-based C cycling observations from the 33-ha treatment forest have been paired with those from a nearby unmanipulated forest. Following over a decade of observations, herein we revisit our core hypothesis: that net ecosystem production (NEP) would increase following the transition to mid-late successional species dominance due to increased canopy structural complexity. Supporting our hypothesis, NEP was stable, briefly declined, and then increased relative to the control in the decade following disturbance; however, increasing NEP was not associated with rising structural complexity but rather the rapid 1-year recovery of total LAI as mid-late successional Acer, Quercus, and Pinus assumed canopy dominance. The transition to mid-late successional species dominance improved carbon-use efficiency (CUE = NEP/gross primary production) as ecosystem respiration declined. Similar soil respiration rates in control and treatment forests, along with species differences in leaf physiology and the rising relative growth rates of mid-late successional species in the treatment forest, suggest changes in aboveground plant respiration and growth, respectively, were primarily responsible for increases in NEP. We conclude that deciduous forests transitioning from early to middle succession are capable of sustained or increased NEP, even when experiencing extensive tree mortality. This adds to mounting evidence that aging deciduous forests in the region will function as C sinks for decades to come.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Utilizing plant synthetic biology to accelerate plant-microbe interactions research

Plant-microbe interactions are critical to ecosystem resilience and substantially influence crop production. From the perspective of plant science, two important focus areas concerning plant-microbe interactions include: 1) understanding plant molecular mechanisms involved in plant-microbe interfaces and 2) engineering plants for increasing plant disease resistance or enhancing beneficial interactions with microbes to increase their resilience to biotic and abiotic stress conditions. Molecular biology and genetics approaches have been used to investigate the molecular mechanisms underlying plant responses to various beneficial and pathogenic microbes. While these approaches are valuable for elucidating the functions of individual genes and pathways, they fall short of unraveling the complex cross-talk across pathways or systems that plants employ to respond and adapt to environmental stresses. Also, genetic engineering of plants to increase disease resistance or enhance symbiosis with microbes has mainly been attempted or conducted through targeted manipulation of single genes/pathways of plants. Recent advancements in synthetic biology tool development are paving the way for multi-gene characterization and engineering in plants in relation to plant-microbe interactions. Here, we briefly summarize the current understanding of plant molecular pathways involved in plant interactions with beneficial and pathogenic microorganisms. Then, we highlight the progress in applying plant synthetic biology to elucidate the molecular basis of plant responses to microbes, enhance plant disease resistance, engineer synthetic symbiosis, and conduct in situ microbiome engineering. Lastly, we discuss the challenges, opportunities, and future directions for advancing plant-microbe interactions research using the capabilities of plant synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Effects of fire and fire-induced changes in soil properties on post-burn soil respiration

Boreal forests cover vast areas of land in the northern hemisphere and store large amounts of carbon (C) both aboveground and belowground. Wildfires, which are a primary ecosystem disturbance of boreal forests, affect soil C via combustion and transformation of organic matter during the fire itself and via changes in plant growth and microbial activity post-fire. Wildfire regimes in many areas of the boreal forests of North America are shifting towards more frequent and severe fires driven by changing climate. As wildfire regimes shift and the effects of fire on belowground microbial community composition are becoming clearer, there is a need to link fire-induced changes in soil properties to changes in microbial functions, such as respiration, in order to better predict the impact of future fires on C cycling. We used laboratory burns to simulate boreal crown fires on both organic-rich and sandy soil cores collected from Wood Buffalo National Park, Alberta, Canada, to measure the effects of burning on soil properties including pH, total C, and total nitrogen (N). We used 70-day soil incubations and two-pool exponential decay models to characterize the impacts of burning and its resulting changes in soil properties on soil respiration. Laboratory burns successfully captured a range of soil temperatures that were realistic for natural wildfire events. We found that burning increased pH and caused small decreases in C:N in organic soil. Overall, respiration per gram total (post-burn) C in burned soil cores was 16% lower than in corresponding unburned control cores, indicating that soil C lost during a burn may be partially offset by burn-induced decreases in respiration rates. Simultaneously, burning altered how remaining C cycled, causing an increase in the proportion of C represented in the modeled slow-cycling vs. fast-cycling C pool as well as an increase in fast-cycling C decomposition rates. Together, our findings imply that C storage in boreal forests following wildfires will be driven by the combination of C losses during the fire itself as well as fire-induced changes to the soil C pool that modulate post-fire respiration rates. Moving forward, we will pair these results with soil microbial community data to understand how fire-induced changes in microbial community composition may influence respiration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Using mid-infrared spectroscopy to estimate soil microbial properties at the continental scale

Understanding microbial community properties is critical to improving the predictions of biogeochemical processes for enhancing soil carbon sequestration. Here, in this observational study, mid-infrared (MIR) spectroscopy and partial least squares regression was used to predict soil microbial and chemical properties from diverse ecosystems across the continental USA. Random calibration and validation demonstrated the prediction potential for soil properties using MIR spectra, with the strongest predictions for microbial respiration, followed by microbial biomass carbon and nitrogen, ß-glucosidase activity, as well as soil chemical properties including organic carbon and total nitrogen. Microbial properties were mainly positively correlated to spectral regions associated with aliphatic C-H groups and C=O stretches of polysaccharides and negatively correlated to quartz and silicate-associated regions. We conclude that MIR spectroscopy can characterize soil microbial functions and be useful for the improvement of continental-scale soil carbon modeling and prediction programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Intra-canopy leaf trait variation facilitates high leaf area index and compensatory growth in a clonal woody encroaching shrub

Leaf trait variation enables plants to utilize large gradients of light availability that exist across canopies of high leaf area index (LAI), allowing for greater net carbon gain while reducing light availability for understory competitors. While these canopy dynamics are well understood in forest ecosystems, studies of canopy structure of woody shrubs in grasslands are lacking. To evaluate the investment strategy used by these shrubs, we investigated the vertical distribution of leaf traits and physiology across canopies of Cornus drummondii, the predominant woody encroaching shrub in the Kansas tallgrass prairie. We also examined the impact of disturbance by browsing and grazing on these factors. Our results reveal that leaf mass per area (LMA) and leaf nitrogen per area (Na) varied approximately threefold across canopies of C. drummondii, resulting in major differences in the physiological functioning of leaves. High LMA leaves had high photosynthetic capacity, while low LMA leaves had a novel strategy for maintaining light compensation points below ambient light levels. The vertical allocation of leaf traits in C. drummondii canopies was also modified in response to browsing, which increased light availability at deeper canopy depths. As a result, LMA and Na increased at lower canopy depths, leading to a greater photosynthetic capacity deeper in browsed canopies compared to control canopies. This response, along with increased light availability, facilitated greater photosynthesis and resource-use efficiency deeper in browsed canopies compared to control canopies. Furthermore, our results illustrate how C. drummondii facilitates high LAI canopies and a compensatory growth response to browsing—both of which are key factors contributing to the success of C. drummondii and other species responsible for grassland woody encroachment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Conterminous USA-Scale Map of Relative Tidal Marsh Elevation

Abstract Tidal wetlands provide myriad ecosystem services across local to global scales. With their uncertain vulnerability or resilience to rising sea levels, there is a need for mapping flooding drivers and vulnerability proxies for these ecosystems at a national scale. However, tidal wetlands in the conterminous USA are diverse with differing elevation gradients, and tidal amplitudes, making broad geographic comparisons difficult. To address this, a national-scale map of relative tidal elevation ( Z * MHW ), a physical metric that normalizes elevation to tidal amplitude at mean high water (MHW), was constructed for the first time at 30 × 30-m resolution spanning the conterminous USA. Contrary to two study hypotheses, watershed-level median Z * MHW and its variability generally increased from north to south as a function of tidal amplitude and relative sea-level rise. These trends were also observed in a reanalysis of ground elevation data from the Pacific Coast by Janousek et al. (Estuaries and Coasts 42 (1): 85–98, 2019). Supporting a third hypothesis, propagated uncertainty in Z * MHW increased from north to south as light detection and ranging (LiDAR) errors had an outsized effect under narrowing tidal amplitudes. The drivers of Z * MHW and its variability are difficult to determine because several potential causal variables are correlated with latitude, but future studies could investigate highest astronomical tide and diurnal high tide inequality as drivers of median Z * MHW and Z * MHW variability, respectively. Watersheds of the Gulf Coast often had propagated Z * MHW uncertainty greater than the tidal amplitude itself emphasizing the diminished practicality of applying Z * MHW as a flooding proxy to microtidal wetlands. Future studies could focus on validating and improving these physical map products and using them for synoptic modeling of tidal wetland carbon dynamics and sea-level rise vulnerability analyses.

Holmquist, James R. (ORCID:0000000325466766)↗
DOE OSTI2026

Data for "Depth of nutrient uptake by deep-rooted plants is regulated by water availability"

The data set consists of strontium (Sr) isotope ratios (87Sr/86Sr), water isotopes, soil cation concentrations, soil water potential sensor data, and results of 87Sr/86Sr mixing model. The plant canopy size files include the dataset of canopy dimension of sagebrush, lupine, and sunflower. The soil and plant ICPMS (Inductively Coupled Plasma Mass Spectrometry) data file includes both of 87Sr/86Sr, and cation concentration dataset from soil exchangeable pool, apatite pool, silicate extract, atmospheric rain deposition, and plant leaf and stem tissues. The plant dendrochronology file includes the dendrochronogical ring width of several sagebrush, and dendrochemical sample data includes the 87Sr/86Sr for each separated growth ring. The modeling result gives the proportion of nutrient sources of each plants (based on their 87Sr/86Sr in leaf tissues and growth rings) from atmospheric deposition and mineral weathering. Soil water potential data includes continuous collection of soil water potential dataset at 2 depths (30 cm and 60 cm, from Nov 24 - Jun 25) of the sampling site. All the samples were collected from 2 sampling campaign June and July 2023, and rain water is a separate sampling from Aug - Sept 2023, at north-facing hillslope near pumphouse site. The data showed that the depth of cation nutrient acquisition is thus tightly coupled with, and likely determined by, water availability in soil, saprolite and bedrock. The enhanced uptake of cations and water from regions of mineral weathering could confer plant and ecosystem resilience during low water years and may impact the rate of bedrock weathering and watershed chemistry during drought. This dataset includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type; a location metadata file (locations.csv); and a samples metadata file (samples.csv). All files are provided as comma-separated values (CSV) files (.csv). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Resolving temperature limitation on spring productivity in an evergreen conifer forest using a model–data fusion framework

Abstract. The flow of carbon through terrestrial ecosystems and the response to climate are critical but highly uncertain processes in the global carbon cycle. However, with a rapidly expanding array of in situ and satellite data, there is an opportunity to improve our mechanistic understanding of the carbon (C) cycle's response to land use and climate change. Uncertainty in temperature limitation on productivity poses a significant challenge to predicting the response of ecosystem carbon fluxes to a changing climate. Here we diagnose and quantitatively resolve environmental limitations on the growing-season onset of gross primary production (GPP) using nearly 2 decades of meteorological and C flux data (2000–2018) at a subalpine evergreen forest in Colorado, USA. We implement the CARbon DAta-MOdel fraMework (CARDAMOM) model–data fusion network to resolve the temperature sensitivity of spring GPP. To capture a GPP temperature limitation – a critical component of the integrated sensitivity of GPP to temperature – we introduced a cold-temperature scaling function in CARDAMOM to regulate photosynthetic productivity. We found that GPP was gradually inhibited at temperatures below 6.0 ∘C (±2.6 ∘C) and completely inhibited below −7.1 ∘C (±1.1 ∘C). The addition of this scaling factor improved the model's ability to replicate spring GPP at interannual and decadal timescales (r=0.88), relative to the nominal CARDAMOM configuration (r=0.47), and improved spring GPP model predictability outside of the data assimilation training period (r=0.88). While cold-temperature limitation has an important influence on spring GPP, it does not have a significant impact on integrated growing-season GPP, revealing that other environmental controls, such as precipitation, play a more important role in annual productivity. This study highlights growing-season onset temperature as a key limiting factor for spring growth in winter-dormant evergreen forests, which is critical in understanding future responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Environmental matrix and moisture are key determinants of microbial phenotypes expressed in a reduced complexity soil-analog

Soil moisture and porosity regulate microbial metabolism by influencing factors such as redox conditions, substrate availability, and soil connectivity. However, the inherent biological, chemical, and physical heterogeneity of soil complicates laboratory investigations into microbial phenotypes that mediate community metabolism. This difficulty arises from challenges in accurately representing the soil environment and in establishing a tractable microbial community that limits confounding variables. To address these challenges in our investigation of community metabolism, we use a reduced-complexity microbial consortium grown in a soil analog using a glass-bead matrix amended with chitin. Long-read and short-read metagenomes, metatranscriptomes, metaproteomes, and metabolomes were analyzed to test the effects of soil structure and moisture on chitin degradation. Our soil structure analog system greatly altered microbial expression profiles compared to the liquid-only incubations, emphasizing the importance of incorporating environmental parameters, like pores and surfaces, for understanding microbial phenotypes relevant to soil ecosystems. These changes were mainly driven by differences in overall expression of chitin-degrading Streptomyces species and stress-tolerant Ensifer. Our findings suggest that the success of Ensifer in a structured environment is likely related to its ability to repurpose carbon via the glyoxylate shunt while potentially using polyhydroxyalkanoate granules as a C source. We also identified traits like motility, stress resistance, and biofilm formation that underlie the degradation of chitin across our treatments and inform how they may ultimately alter carbon use efficiency. Together our results demonstrate that community functions like decomposition are sensitive to environmental conditions and more complex than the multi-enzyme pathways involved in depolymerization.

Rodriguez-Ramos, Josue A↗
DOE OSTI2026

The Global Spectra-Trait Initiative: A database of paired leaf spectroscopy and functional traits associated with leaf photosynthetic capacity

Accurate assessment of leaf functional traits is crucial for a diverse range of applications from crop phenotyping to parameterizing global climate models. Leaf reflectance spectroscopy offers a promising avenue to advance ecological and agricultural research by complementing traditional, time-consuming gas exchange measurements. However, the development of robust hyperspectral models for predicting leaf photosynthetic capacity and associated traits from reflectance data has been hindered by limited data availability across species and environments. Here we introduce the Global Spectra-Trait Initiative (GSTI), a collaborative repository of paired leaf hyperspectral and gas exchange measurements from diverse ecosystems. The GSTI repository currently encompasses over 7500 observations from 397 species and 41 sites gathered from 36 published and unpublished studies, thereby offering a key resource for developing and validating hyperspectral models of leaf photosynthetic capacity. The GSTI database is developed on GitHub (https://github.com/plantphys/gsti, last access: 4 January 2026) and published to ESS-DIVE https://doi.org/10.15485/2530733, Lamour et al., 2025). It includes gas exchange data, derived photosynthetic parameters, and key leaf traits often associated with traditional gas exchange measurements such as leaf mass per area and leaf elemental composition. By providing a standardized repository for data sharing and analysis, we present a critical step towards creating hyperspectral models for predicting photosynthetic traits and associated leaf traits for terrestrial plants.

Lamour, Julien [Université of Toulouse (France); U↗
DOE OSTI2020

A Genomic Catalogue of Earth’s Microbiomes - Introductory KBase Narrative

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled unprecedented insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this powerful approach to over 10,000 metagenomes collected from diverse habitats covering all of Earth’s continents and oceans, human- and animal-host associated microbiomes, engineered environments, and natural and agricultural soils to capture extant microbial metabolic and functional potential. We present a comprehensive catalogue of 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units (OTUs), spanning 135 phyla, which expand the known phylogenetic diversity of Bacteria and Archaea by 44%. We also demonstrate the utility of this collection for secondary metabolite biosynthetic potential and predicting host-virus linkage, which can provide a view into the global distribution of lysogenic viruses. This resource underscores the value of leveraging genome-centric approaches to reveal genomic properties of uncultivated microbes that impact on ecosystem processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The Straightening of a River Meander Leads to Extensive Losses in Flow Complexity and Ecosystem Services

To assist river restoration efforts we need to slow the rate of river degradation. This study provides a detailed explanation of the hydraulic complexity loss when a meandering river is straightened in order to motivate the protection of river channel curvature. We used computational fluid dynamics (CFD) modeling to document the difference in flow dynamics in nine simulations with channel curvature (C) degrading from a well-established tight meander bend (C = 0.77) to a straight channel without curvature (C = 0). To control for covariates and slow the rate of loss to hydraulic complexity, each of the nine-channel realizations had equivalent bedform topography. The analyzed hydraulic variables included the flow surface elevation, streamwise and transverse unit discharge, flow velocity at streamwise, transverse, and vertical directions, bed shear stress, stream function, and the vertical hyporheic flux rates at the channel bed. The loss of hydraulic complexity occurred gradually when initially straightening the channel from C = 0.77 to C = 0.33 (i.e., the radius of the channel is three-times the channel width), and additional straightening incurred rapid losses to hydraulic complexity. Other studies have shown hydraulic complexity provides important riverine habitat and is positively correlated with biodiversity. This study demonstrates how hydraulic complexity can be gradually and then rapidly lost when unwinding a river, and hopefully will serve as a cautionary tale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗