Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability prediction”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Predicting responses to climate change using a joint species, spatially dependent physiologically guided abundance model

Abstract Predicting the effects of warming temperatures on the abundance and distribution of organisms under future climate scenarios often requires extrapolating species–environment correlations to climatic conditions not currently experienced by a species, which can result in unrealistic predictions. For poikilotherms, incorporating species' thermal physiology to inform extrapolations under novel thermal conditions can result in more realistic predictions. Furthermore, models that incorporate species and spatial dependencies may improve predictions by capturing correlations present in ecological data that are not accounted for by predictor variables. Here, we present a joint species, spatially dependent physiologically guided abundance (jsPGA) model for predicting multispecies responses to climate warming. The jsPGA model uses a basis function approach to capture both species and spatial dependencies. We apply the jsPGA model to predict the response of eight fish species to projected climate warming in thousands of lakes in Minnesota, USA. By the end of the century, the cold‐adapted species was predicted to have high probabilities of extirpation across its current range—with 10% of lakes currently inhabited by this species having an extirpation probability >0.90. The remaining species had varying levels of predicted changes in abundance, reflecting differences in their thermal physiology. Though the model did not identify many strong species dependencies, the variation in estimated spatial dependence across species suggested that accounting for both dependencies was important for predicting the abundance of these fishes. The jsPGA model provides a new tool for predicting changes in the abundance, distribution, and extirpation probability of poikilotherms under novel thermal conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiscale Simulation of Precipitation Over East Asia by Variable Resolution CAM-MPAS

This study evaluates the precipitation over East Asia simulated by the variable resolution CAM-MPAS featuring the Model for Prediction Across Scales (MPAS) atmospheric dynamical core coupled with the physics parameterizations of the Community Atmosphere Model (CAM) version 5.4. Two CAM-MPAS experiments, one with a global quasi-uniform resolution mesh of 120 km (MPAS-UR) and the other with a variable resolution mesh of 30 to 120 km refined over East Asia (MPAS-VR), are conducted from 1989 to 2005 following the Atmospheric Model Intercomparison Project protocol. Two regional climate model simulations at ~25 km resolution from the Coordinated Regional Downscaling Experiment East Asia second plan are also analyzed for comparison. Results show that CAM-MPAS performs better than the selected regional models in simulating the precipitation climatology over East Asia. Compared with MPAS-UR, MPAS-VR with refinement over East Asia better simulates the precipitation over and around the Tibetan Plateau (TP), the frequency distribution of light, moderate, and heavy precipitation in various regions of China, and the interannual precipitation pattern associated with the East Asian summer monsoon. The improvements in MPAS-VR relative to MPAS-UR largely come from the decreased grid spacing rather than the use of finer-scale terrain information. Particularly, decreased grid spacing is more essential for simulating the precipitation magnitude over TP and the location of peak precipitation south of TP, likely due to the better resolved physical and dynamical processes associated with orographic precipitation. The 30-120 km CAM-MPAS model is shown to be a promising tool for precipitation simulation over East Asia.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Can socio-economic indicators of vulnerability help predict spatial variations in the duration and severity of power outages due to tropical cyclones?

Abstract Tropical cyclones are the leading cause of major power outages in the U.S., and their effects can be devastating for communities. However, few studies have holistically examined the degree to which socio-economic variables can explain spatial variations in disruptions and reveal potential inequities thereof. Here, we apply machine learning techniques to analyze 20 tropical cyclones and predict county-level outage duration and percentage of customers losing power using a comprehensive set of weather, environmental, and socio-economic factors. Our models are able to accurately predict these outage response variables, but after controlling for the effects of weather conditions and environmental factors in the models, we find the effects of socio-economic variables to be largely immaterial. However, county-level data could be overlooking effects of socio-economic disparities taking place at more granular spatial scales, and we must remain aware of the fact that when faced with similar outage events, socio-economically vulnerable communities will still find it more difficult to cope with disruptions compared to less vulnerable ones.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Predicted range expansion of Prostephanus truncatus (Coleoptera: Bostrichidae) under projected climate change scenarios

Abstract The larger grain borer (Prostephanus truncatus [Horn] [Coleoptera: Bostrichidae]) is a wood-boring insect native to Central America and adapted to stored maize and cassava. It was accidentally introduced to Tanzania and became a pest across central Africa. Unlike many grain pests, P. truncatus populations can establish and move within forests. Consequently, novel infestations can occur without human influence. The objectives of our study were to (i) develop an updated current suitability projection for P. truncatus, (ii) assess its potential future distribution under different climate change scenarios, and (iii) identify climate variables that best inform the model. We used WALLACE and MaxEnt to predict potential global distribution by incorporating bioclimatic variables and occurrence records. Future models were projected for 2050 and 2070 with Representative Concentration Pathways (RCPs) 2.6 (low change) and 8.5 (high change). Distribution was most limited by high precipitation and cold temperatures. Globally, highly suitable areas (> 75%) primarily occurred along coastal and equatorial regions with novel areas in northern South America, India, southeastern Asia, Indonesia, and the Philippines, totaling 7% under current conditions. Highly suitable areas at RCPs 2.6 and 8.5 are estimated to increase to 12% and 15%, respectively, by 2050 and increase to 19% in 2070 under RCP 8.5. Centroids of highly suitable areas show distribution centers moving more inshore and away from the equator. Notably, the result is a range expansion, not a shift. Results can be used to decrease biosecurity risks through more spatially explicit and timely surveillance programs for targeting the exclusion of this pest.

Entomology↗
DOE OSTI2021

Results from a wake-steering experiment at a commercial wind plant: investigating the wind speed dependence of wake-steering performance

Wake steering is a wind farm control strategy in which upstream wind turbines are misaligned with the wind to redirect their wakes away from downstream turbines, thereby increasing the net wind plant power production and reducing fatigue loads generated by wake turbulence. In this paper, we present results from a wake-steering experiment at a commercial wind plant involving two wind turbines spaced 3.7 rotor diameters apart. During the 3-month experiment period, we estimate that wake steering reduced wake losses by 5.6% for the wind direction sector investigated. After applying a long-term correction based on the site wind rose, the reduction in wake losses increases to 9.3%. As a function of wind speed, we find large energy improvements near cut-in wind speed, where wake steering can prevent the downstream wind turbine from shutting down. Yet for wind speeds between 6–8 m/s, we observe little change in performance with wake steering. However, wake steering was found to improve energy production significantly for below-rated wind speeds from 8–12 m/s. By measuring the relationship between yaw misalignment and power production using a nacelle lidar, we attribute much of the improvement in wake-steering performance at higher wind speeds to a significant reduction in the power loss of the upstream turbine as wind speed increases. Additionally, we find higher wind direction variability at lower wind speeds, which contributes to poor performance in the 6–8 m/s wind speed bin because of slow yaw controller dynamics. Further, we compare the measured performance of wake steering to predictions using the FLORIS (FLOw Redirection and Induction in Steady State) wind farm control tool coupled with a wind direction variability model. Although the achieved yaw offsets at the upstream wind turbine fall short of the intended yaw offsets, we find that they are predicted well by the wind direction variability model. When incorporating the expected yaw offsets, estimates of the energy improvement from wake steering using FLORIS closely match the experimental results.

17 WIND ENERGY↗
DOE OSTI2021

Decadal Variability of Southeast US Rainfall in an Eddying Global Coupled Model

Abstract Ocean variability is a dominant source of remote rainfall predictability, but in many cases the physical mechanisms driving this predictability are not fully understood. This study examines how ocean mesoscales (i.e., the Gulf Stream SST front) affect decadal Southeast US (SEUS) rainfall, arguing that the local imprint of large‐scale teleconnections is sensitive to resolved mesoscale features. Based on global coupled model experiments with eddying and eddy‐parameterizing ocean, we find that a resolved Gulf Stream improves localized rainfall and remote circulation response in the SEUS. The eddying model generally improves the air‐sea interactions in the Gulf Stream and the North Atlantic Subtropical High that modulate SEUS rainfall over decadal timescales. The eddy‐parameterizing simulation fails to capture the sharp SST gradient associated with the Gulf Stream and overestimates the role of tropical Pacific SST anomalies in the SEUS rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2021

Spatial and temporal variability of tritium in precipitation within South Africa and it's bearing on hydrological studies

Tritium, the radioactive isotope of hydrogen, has been used to understand groundwater recharge processes for decades. The current variation of tritium in the atmosphere is largely attributed to stratospheric production and fall out rates as well as global circulation phenomena controlling the hydrological cycle. Global controls on the variability in atmospheric tritium activity are poorly suited to explain local variation and tritium activities in precipitation are often assumed to be uniform over both local and regional catchments and watersheds. This assumption can result in both over and under estimation of modern recharge within an aquifer when using tritium as the recharge proxy. In order to minimize the inherent prediction residuals associated with tritium based recharge investigations, the variability of tritium in precipitation was modelled from 127 spatial precipitation samples taken over a two year period, combined with a 76 precipitation sample group-set taken over a one year period in a single location. Precipitation events were traced backward in time, from the point of collection, using HYSPLIT modelling to ascertain the origins of moisture content as well as the altitudes of moisture origin reached along the particle track. Tritium activities, collected over a one year period in Paarl, range from 0.45 to 4.16 TU and have a mean of 1.59 TU. Spatial storm events in the Western Cape in 2017 and 2018 had a range from 0 to 2.2 and 0.37 to 3.27 TU, respectively, with mean activities of 1.18 (n = 34) and 1.25 TU (n = 32). Both storm events had similar tritium variability (σ = 0.5 n = 35 and 0.48 n = 32). Regional precipitation events had the largest range of tritium activities (0.55–12.2 TU). Although not all tritium activities can be explained by interrogating the water mass origin, this study suggests that approximately 90% of events can be completely or partially attributed to the origin of the water mass. The variability of tritium, both spatially and temporally, was higher than expected, confirming that when uniform tritium inputs are used, the groundwater system would provide inaccurate modern recharge estimates. Higher spatial resolution of tritium variation in precipitation for a particular region will improve our ability to relate tritium activities in groundwater to local precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Reduced ecosystem resilience quantifies fine-scale heterogeneity in tropical forest mortality responses to drought

Sensitivity of forest mortality to drought in carbon-dense tropical forests remains fraught with uncertainty, while extreme droughts are predicted to be more frequent and intense. Here, in this study, the potential of temporal autocorrelation of high-frequency variability in Landsat Enhanced Vegetation Index (EVI), an indicator of ecosystem resilience, to predict spatial and temporal variations of forest biomass mortality is evaluated against in situ census observations for 64 site-year combinations in Costa Rican tropical dry forests during the 2015 ENSO drought. Temporal autocorrelation, within the optimal moving window of 24 months, demonstrated robust predictive power for in situ mortality (leave-one- out cross-validation R 2 = 0.54), which allows for estimates of annual biomass mortality patterns at 30 m resolution. Subsequent spatial analysis showed substantial fine-scale heterogeneity of forest mortality patterns, largely driven by drought intensity and ecosystem properties related to plant water use such as forest deciduousness and topography. Highly deciduous forest patches demonstrated much lower mortality sensitivity to drought stress than less deciduous forest patches after elevation was controlled. Our results highlight the potential of high-resolution remote sensing to “fingerprint” forest mortality and the significant role of ecosystem heterogeneity in forest biomass resistance to drought.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Impact of volcanic eruptions on CMIP6 decadal predictions: a multi-model analysis

Abstract. In recent decades, three major volcanic eruptions of different intensity have occurred (Mount Agung in 1963, El Chichón in 1982 and Mount Pinatubo in 1991), with reported climate impacts on seasonal to decadal timescales that could have been potentially predicted with accurate and timely estimates of the associated stratospheric aerosol loads. The Decadal Climate Prediction Project component C (DCPP-C) includes a protocol to investigate the impact of volcanic aerosols on the climate experienced during the years that followed those eruptions through the use of decadal predictions. The interest of conducting this exercise with climate predictions is that, thanks to the initialisation, they start from the observed climate conditions at the time of the eruptions, which helps to disentangle the climatic changes due to the initial conditions and internal variability from the volcanic forcing. The protocol consists of repeating the retrospective predictions that are initialised just before the last three major volcanic eruptions but without the inclusion of their volcanic forcing, which are then compared with the baseline predictions to disentangle the simulated volcanic effects upon climate. We present the results from six Coupled Model Intercomparison Project Phase 6 (CMIP6) decadal prediction systems. These systems show strong agreement in predicting the well-known post-volcanic radiative effects following the three eruptions, which induce a long-lasting cooling in the ocean. Furthermore, the multi-model multi-eruption composite is consistent with previous work reporting an acceleration of the Northern Hemisphere polar vortex and the development of El Niño conditions the first year after the eruption, followed by a strengthening of the Atlantic Meridional Overturning Circulation the subsequent years. Our analysis reveals that all these dynamical responses are both model- and eruption-dependent. A novel aspect of this study is that we also assess whether the volcanic forcing improves the realism of the predictions. Comparing the predicted surface temperature anomalies in the two sets of hindcasts (with and without volcanic forcing) with observations we show that, overall, including the volcanic forcing results in better predictions. The volcanic forcing is found to be particularly relevant for reproducing the observed sea surface temperature (SST) variability in the North Atlantic Ocean following the 1991 eruption of Pinatubo.

Bilbao, Roberto (ORCID:0000000307294980)↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2021

Integrated End-to-end Performance Prediction and Diagnosis for Extreme Scientific Workflows (IPPD) (Final Report)

This report details the accomplishments from the ASCR funded project “Integrated End-to-end Performance Prediction and Diagnosis for Extreme Scientific Workflows” under the award numbers FWP-66406 and DE-SC0012630, with a focus on the UC San Diego (Award No. DE-SC0012630) part of the accomplishments. We refer to the project as IPPD. The main activities of IPPD were centered on the development and integration of provenance information to capture empirically workflow information and identify the sources of bottlenecks as well as variability, and modeling and simulation to provide insights and predictively explore multiple scenarios for development of advanced techniques for optimization of resources and workflow execution. After a Phase I of the project, a Phase II effort focused on three major aspects: a) observe how data is generated, distributed, and used; b) analyze how data is (repeatedly) consumed with a focus both on repeated patterns and anomalies; and c) explore how to optimize data motion. The project leveraged and extended our existing tools with new research and demonstrated our work on the Belle II workflow suite as well as on workflows from NSLS-II.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

HȮ 2 + HȮ 2 : High level theory and the role of singlet channels

At low temperatures and high pressures, the HȮ 2 + HȮ 2 reaction is an important reaction in combustion. There are significant unresolved discrepancies between prior theoretical and experimental studies of this reaction. It has generally been presumed to occur as an abstraction on the triplet surface to produce H 2 O 2 + 3 O 2 . Here, we employ a combination of high-level electronic structure theory (the ANL0 composite method or multi-reference methods as appropriate), sophisticated transition state theory (vibrationally adiabatic torsions, variational, and variable reaction coordinate), and master equation analyses to predict the thermal kinetics on the H 2 O 4 surface. Notably, this analysis suggests a significant branching to HȮ 3 + ȮH near 1000 K via reaction on the singlet surface. This channel does not appear to have been considered in prior combustion models. HȮ 3 itself is a metastable complex (bound by only 3 kcal mol -1 ) that rapidly dissociates to ȮH + 3 O 2 . Thus, the net reaction for this channel, HȮ 2 + HȮ 2 → ȮH + ȮH + O 2 , converts two low reactivity HȮ 2 radicals into two highly reactive ȮH radicals. The ramifications of these newly derived rate expressions are highlighted through kinetic modeling studies of H 2 , CH 3 OH, $\textit{n}$-heptane, and isooctane oxidation; all at the high pressures of relevance to combustion devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A new phenomenological model to describe root-soil interactions based on percolation theory

In his paper on net primary productivity of terrestrial communities predicted from climatological data, Rosenzweig (1968) argued that variability in productivity is well accounted for by (evapo)-transpiration, and that water from transpiration is, on global scales, the most variable component in the photosynthesis reaction. The goal of this paper is to investigate whether variability in plant growth on local scales and within species is primarily related to transpiration under several scenarios including different terrain curvature, slope aspect, soil characteristics, and climate ranges. Here, we test the hypothesis that this relationship exists because root growth into the surface soil layers (0–2 m) tends to follow paths with minima in resistance, which in turn maximizes water flow and nutrient delivery rates that regulate growth. The set of all connected paths with individual pore-to-pore flow resistances less than a critical, percolating, value forms a cluster with mass fractal dimensionality, d f . We propose that roots follow paths through the 2D percolation cluster, defining the set of all optimal flow paths, making the 2D value of d f from percolation relevant to root fractal dimensionality. The tortuosity of such optimal paths as defined in percolation theory should then relate root length to root radial extent, linking the parameters of root tortuosity and plant productivity. Our analysis of large data sets across species implies that root radial extent and tree height are both proportional to cumulative transpiration until trees approached maximum height, and their growth rates are proportional to the transpiration rate, not to the moisture content. Local variations in tree height as functions of the variables investigated appear generally consistent with deduced variations in transpiration. Here this correlation is investigated more closely in the context of studies addressing individual tree species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Role of pH and Eh in geothermal systems: Thermodynamic examples and impacts on scaling and corrosion

The global demand for sustainable and clean renewable energy sources is rapidly increasing to meet growing energy demands while curtailing the CO2 emissions from fossil fuels. Harnessing the natural heat or geothermal energy of the Earth's interior can simultaneously address both issues by offering a clean and renewable energy source. Enhanced geothermal systems require the movement and cycling of fluid through a porous rock media. Further, this process requires understanding the chemistry occurring at the fluid-rock and fluid-metal interfaces both in the deep subsurface and surface, where fluid is circulated for energy transfer or storage. In this review paper, we evaluate the current state of research on how two key variables, pH and Eh, control the fluid-rock reactions controlling mineral scaling and the steel corrosion reactions associated with geothermal systems. These variables are frequently used in geochemical modeling to predict chemical reactions in geothermal systems. Prediction of the fluid-rock behavior can be utilized to formulate new strategies for extending the lifespan and efficiency of heat extraction and utilization and maximizing the safety of geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Shifting groundwater fluxes in bedrock fractures: Evidence from stream water radon and water isotopes

Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and below the fractured hillslope (non-fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone dominated, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. Importantly, we interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The effects of mixed-species root zones on the resistance of soil bacteria and fungi to long-term experimental and natural reductions in soil moisture

Mixed forest stands tend to be more resistant to drought than species-specific stands partially due to complementarity in root ecology and physiology. For this work, we asked whether complementary differences in the drought resistance of soil microbiomes might contribute to this phenomenon. We experimented on the effects of reduced soil moisture on bacterial and fungal community composition in species-specific (single species) and mixed-species root zones of Norway spruce and European beech forests in a 5-year-old throughfall-exclusion experiment and across seasonal (spring-summer-fall) and latitudinal moisture gradients. Bacteria were most responsive to changes in soil moisture, especially members of Rhizobiales, while fungi were largely unaffected, including ectomycorrhizal fungi (EMF). Community resistance was higher in spruce relative to beech root zones, corresponding with the proportions of drought-favored (more in spruce) and drought-sensitive bacterial taxa (more in beech). The spruce soil microbiome also exhibited greater resistance to seasonal changes between spring (wettest) and fall (driest). Mixed-species root zones contained a hybrid of beech- and spruce-associated microbiomes. Several bacterial populations exhibited either enhanced resistance or greater susceptibility to drought in mixed root zones. Overall, patterns in the relative abundances of soil bacteria closely tracked moisture in seasonal and latitudinal precipitation gradients and were more predictive of soil water content than other environmental variables. We conclude that complementary differences in the drought resistance of soil microbiomes can occur and the likeliest form of complementarity in mixed-root zones coincides with the enrichment of drought-tolerant bacteria associated with spruce and the sustenance of EMF by beech.

54 ENVIRONMENTAL SCIENCES↗