Engineering topics

Williams, Kenneth H.

Publications and source records attributed to Williams, Kenneth H..

At least 19 records

DOE OSTI2024

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Stream water sourcing from high-elevation snowpack inferred from stable isotopes of water: a novel application of d-excess values

About 80 % of the precipitation at the Colorado River's headwaters is snow, and the resulting snowmelt-driven hydrograph is a crucial water source for about 40 million people. Snowmelt from alpine and subalpine snowpack contributes substantially to groundwater recharge and river flow. However, the dynamics of snowmelt progression are not well understood because observations of the high-elevation snowpack are difficult due to challenging access in complex mountainous terrain as well as the cost and labor intensity of currently available methods. We present a novel approach to infer the processes and dynamics of high-elevation snowmelt contributions predicated upon stable hydrogen and oxygen isotope ratios observed in streamflow. We show that deuterium-excess (d-excess) values of stream water could serve as a comparatively cost-effective proxy for a catchment-integrated signal of high-elevation snowmelt contributions to catchment runoff. We sampled stable hydrogen and oxygen isotope ratios of the precipitation, snowpack, and stream water in the East River, a headwater catchment of the Colorado River, and the stream water of larger catchments at sites on the Gunnison River and Colorado River. The d-excess of snowpack increased with elevation; the upper subalpine and alpine snowpack (> 3200 m) had substantially higher d-excess compared to lower elevations (< 3200 m) in the study area. The d-excess values of stream water reflected this because d-excess values increased as the higher-elevation snowpack contributed more to stream water generation later in the snowmelt/runoff season. End-member mixing analyses based on the d-excess data showed that the share of high-elevation snowmelt contributions within the snowmelt hydrograph was on average 44 % and generally increased during melt period progression, up to 70 %. The observed pattern was consistent during 6 years for the East River, and a similar relation was found for the larger catchments on the Gunnison and Colorado rivers. High-elevation snowpack contributions were found to be higher for years with lower snowpack and warmer spring temperatures. Thus, we conclude that the d-excess of stream water is a viable proxy to observe changes in high-elevation snowmelt contributions in catchments at various scales. Inter-catchment comparisons and temporal trends of the d-excess of stream water could therefore serve as a catchment-integrated measure to monitor if mountain systems rely on high-elevation water inputs more during snow drought compared to years of average snowpack depths.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Old-Aged groundwater contributes to mountain hillslope hydrologic dynamics

Understanding connectivity between the soil and deeper bedrock groundwater is needed to accurately predict a watershed’s response to perturbation, such as drought. Yet, the bedrock groundwater dynamics in mountainous environments are typically under-constrained and excluded from watershed hydrologic models. Here, we investigate the role of groundwater characterized with decadal and longer water ages on the hydrologic and mass-transport processes within a steep snow-dominated mountain hillslope in the Central Rocky Mountains (USA). We quantify subsurface and surface water mass-balance, groundwater flowpaths, and age distributions using the ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models, which are compared to hydrometric and environmental tracer observations. An ensemble of models with varied soil and hydrogeologic parameters reproduces observed groundwater levels and century-scale mean ages inferred from environmental tracers. The numerical models suggest soil water near the toe of the hillslope contains considerable (>60 % of the mass-flux) contributions from bedrock flowpaths characterized with water ages >10 years. Flowpath connectivity between the deeper bedrock and soil systems is present throughout the year, highlighting the potentially critical role of groundwater with old ages on processes such as evapotranspiration and streamflow generation. The coupled numerical model and groundwater age observations show the bedrock groundwater system influences the hillslope hydrodynamics and should be considered in mountain watershed conceptual and numerical models.

54 ENVIRONMENTAL SCIENCES↗

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

Leveraging Groundwater Dynamics to Improve Predictions of Summer Low-Flow Discharges

Summer streamflow predictions are critical for managing water resources; however, warming-induced shifts from snow to rain regimes impact low-flow predictive models. Additionally, reductions in snowpack drive earlier peak flows and lower summer flows across the western United States increasing reliance on groundwater for maintaining summer streamflow. However, it remains poorly understood how groundwater contributions vary interannually. Here, we quantify recession limb groundwater (RLGW), defined as the proportional groundwater contribution to the stream during the period between peak stream flow and low flow, to predict summer low flows across three diverse western US watersheds. We ask (a) how do snow and rain dynamics influence interannual variations of RLGW contributions and summer low flows?; (b) which watershed attributes impact the effectiveness of RLGW as a predictor of summer low flows? Linear models reveal that RLGW is a strong predictor of low flows across all sites and drastically improves low-flow prediction compared to snow metrics at a rain-dominated site. Results suggest that strength of RLGW control on summer low flows may be mediated by subsurface storage. Subsurface storage can be divided into dynamic (i.e., variability saturated) and deep (i.e., permanently saturated) components, and we hypothesize that interannual variability in dynamic storage contribution to streamflow drives RLGW variability. In systems with a higher proportion of dynamic storage, RLGW is a better predictor of summer low flow because the stream is more responsive to dynamic storage contributions compared to deep-storage-dominated systems. Overall, including RLGW improved low-flow prediction across diverse watersheds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Advanced monitoring of soil-vegetation co-dynamics reveals the successive controls of snowmelt on soil moisture and on plant seasonal dynamics in a mountainous watershed

Evaluating the interactions between above- and below-ground processes is important to understand and quantify how ecosystems respond differently to atmospheric forcings and/or perturbations and how this depends on their intrinsic characteristics and heterogeneity. Improving such understanding is particularly needed in snow-impacted mountainous systems where the complexity in water and carbon storage and release arises from strong heterogeneity in meteorological forcing and terrain, vegetation and soil characteristics. This study investigates spatial and temporal interactions between terrain, soil moisture, and plant seasonal dynamics at the intra- and inter-annual scale along a 160 m long mountainous, non-forested hillslope-to-floodplain system in the upper East River Watershed in the upper Colorado River Basin. To this end, repeated UAV-based multi-spectral aerial imaging, ground-based soil electrical resistivity imaging, and soil moisture sensors were used to quantify the interactions between above and below-ground compartments. Results reveal significant soil-plant co-dynamics. The spatial variation and dynamics of soil water content and electrical conductivity, driven by topographic and soil intrinsic characteristics, correspond to distinct plant types, with highest plant productivity in convergent areas. Plant productivity in heavy snow years benefited from more water infiltration as well as a shallow groundwater table depth. In comparison, low snowpack years with an early first bare-ground date, which are linked to an early increase in plant greenness, imply a short period of saturated conditions that leads to lower average and maximum greenness values during the growing season. Overall, these results emphasize the strong impact of snowpack dynamics, and terrain and subsurface characteristics on the heterogeneity in plant type and seasonal dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

River thorium concentrations can record bedrock fracture processes including some triggered by distant seismic events

Fractures are integral to the hydrology and geochemistry of watersheds, but our understanding of fracture dynamics is very limited because of the challenge of monitoring the subsurface. Here we provide evidence that long-term, high-frequency measurements of the river concentration of the ultra-trace element thorium (Th) can provide a signature of bedrock fracture processes spanning neighboring watersheds in Colorado. River Th concentrations show abrupt (subdaily) excursions and biexponential decay with approximately 1-day and 1-week time constants, concentration patterns that are distinct from all other solutes except beryllium and arsenic. The patterns are uncorrelated with daily precipitation records or seasonal trends in atmospheric deposition. Groundwater Th analyses are consistent with bedrock release and dilution upon mixing with river water. Most Th excursions have no seismic signatures that are detectable 50 km from the site, suggesting the Th concentrations can reveal aseismic fracture or fault events. We find, however, a weak statistical correlation between Th and seismic motion caused by distant earthquakes, possibly the first chemical signature of dynamic earthquake triggering, a phenomenon previously identified only through geophysical methods.

58 GEOSCIENCES↗

DOE OSTI2023

Surface Water Disinfection Byproducts and Organic Matter Characterization Data Associated with: “Disinfection byproducts formed during drinking water treatment reveal an export control point for dissolved organic matter in a subalpine headwater stream”

This dataset is associated with the publication “Disinfection byproducts formed during drinking water treatment reveal an export control point for dissolved organic matter in a subalpine headwater stream” published in Water Research X (Leonard et al. 2022; https://doi.org/10.1016/j.wroa.2022.100144). The associated study analyzed temporal trends from the Town of Crested Butte water treatment facility and synoptic sampling at Coal Creek in Crested Butte, Colorado, US. This work demonstrates how drinking water quality archives combined with synoptic sampling and targeted analyses can be used to identify and understand export control points for dissolved organic matter. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) metadata and international geo-sample number (IGSN) mapping file; (4) dissolved organic carbon (DOC), ultraviolet absorbance at 254 nanometers (UV254), total nitrogen (TN), and specific ultraviolet absorbance (SUVA) data; (5) disinfection bioproduct formation potential (DBP-FP) data; (6) readme; (7) methods codes; (8) water collection protocol; (9) folder of high resolution characterization of organic matter via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory); and (10) folder of excitation emissions matrix (EEM) spectra. The FTICR folder contains a file of DOC (measured as non-purgeable organic carbon; NPOC) used for FTICR sample preparation. The FTICR folder also contains three subfolders: one subfolder containing the raw .xml data files, one containing the processed data, and the other containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .dat, .R, .ref, or .xml

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗

DOE OSTI2022

Borgs are giant genetic elements with potential to expand metabolic capacity

Anaerobic methane oxidation exerts a key control on greenhouse gas emissions1, yet factors that modulate the activity of microorganisms performing this function remain poorly understood. Here we discovered extraordinarily large, diverse DNA sequences that primarily encode hypothetical proteins through studying groundwater, sediments and wetland soil where methane production and oxidation occur. Four curated, complete genomes are linear, up to approximately 1 Mb in length and share genome organization, including replichore structure, long inverted terminal repeats and genome-wide unique perfect tandem direct repeats that are intergenic or generate amino acid repeats. We infer that these are highly divergent archaeal extrachromosomal elements with a distinct evolutionary origin. Gene sequence similarity, phylogeny and local divergence of sequence composition indicate that many of their genes were assimilated from methane-oxidizing Methanoperedens archaea. We refer to these elements as ‘Borgs’. We identified at least 19 different Borg types coexisting with Methanoperedens spp. in four distinct ecosystems. Borgs provide methane-oxidizing Methanoperedens archaea access to genes encoding proteins involved in redox reactions and energy conservation (for example, clusters of multihaem cytochromes and methyl coenzyme M reductase). These data suggest that Borgs might have previously unrecognized roles in the metabolism of this group of archaea, which are known to modulate greenhouse gas emissions, but further studies are now needed to establish their functional relevance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Modeling Snow Dynamics and Stable Water Isotopes Across Mountain Landscapes

Abstract A coupled hydrologic and snowpack stable water isotope model assesses controls on isotopic inputs across a mountainous basin. Annually, the most depleted isotope conditions occur in the upper subalpine where snow accumulation is high, and rainfall is low. Snowmelt isotopic evolution over time indicates fractionation processes account for <25% snowmelt enrichment. Meltwater isotopic inputs are largely determined by controls on the amount, phase and isotopic mass of precipitation coincident with the ablation period. Effect of vapor loss from the snowpack on d‐excess in snowmelt is a balance between energy and snow‐availability. It is highest above treeline, and in the grass and aspen‐dominated portions of the upper montane where vegetation shading is low. Deep snowpack in conifer forests limit the influence of vapor loss in the subalpine. Wet years reduce the effects of vapor loss on snowmelt across the basin, except in the lower montane where added snowfall bolsters snow‐limited conditions.

58 GEOSCIENCES↗

DOE OSTI2022

Evaluation of hydrograph separation techniques with uncertain end-member composition

Hydrograph separation is one of many approaches used to analyse shifts in source water contributions to stream flow resulting from climate change in remote watersheds. Understanding these shifts is vital, as shifts in source water contributions to a stream can shape water management decisions. Because remote watersheds are often inaccessible and have poorly characterized contributing water sources, or end-members, it is critical to understand the implications of using different hydrograph separation techniques in these data-limited environments. To explore the uncertainty associated with different techniques, results from two hydrograph separation techniques, mass balance and principle component analysis, were compared using 3 years of aqueous geochemical data from the East River watershed located in the Elk Mountains of Central Colorado. Solute concentrations of the end-members were characterized by both a limited set of direct chemical measurements of different sources and detailed seasonal instream chemistry to examine the influences of uncertain end-member compositions in a data-limited environment. Annual volumetric end-member contributions to stream flow had relatively good agreement across separation techniques. Large variations in time were observed in the hydrograph separations, depending on the end-member type, and estimated flow contributions varied between the selected solutes. End-member concentrations characterized by stream chemistry showed several limitations including a reduced number of distinguishable end-members and differences in timing of flow contributions. Here the results highlight the benefits of using multiple hydrograph separation techniques by providing a ‘weight-of-evidence’ approach to environments with limited end-member concentration data.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Variability in observed stable water isotopes in snowpack across a mountainous watershed in Colorado

In this study, isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Variability of Snow and Rainfall Partitioning Into Evapotranspiration and Summer Runoff Across Nine Mountainous Catchments

Abstract Understanding the partitioning of snow and rain contributing to either catchment streamflow or evapotranspiration (ET) is of critical relevance for water management in response to climate change. To investigate this partitioning, we use endmember splitting and mixing analyses based on stable isotope ( 18 O) data from nine headwater catchments in the East River, Colorado. Our results show that one third of the snow partitions to ET and 13% of the snowmelt sustains summer streamflow. Only 8% of the rainfall contributes to the summer streamflow, because most of the rain (67%) partitions to ET. The spatial variability of precipitation partitioning is mainly driven by aspect and tree cover across the sub‐catchments. Catchments with higher tree cover have a higher share of snow becoming ET, resulting in less snow in summer streamflow. Summer streamflow did not contain more rain with higher rainfall sums, but more rain was taken up in ET.

58 GEOSCIENCES↗

DOE OSTI2022

Sulfur Biogeochemical Cycling and Redox Dynamics in a Shale-Dominated Mountainous Watershed

Sulfur (S) is an essential macronutrient and important component of the earth’s crust, and its cycling has critical impacts on trace metal mobility, water quality, and human health. Pyrite weathering is the primary pathway by which sulfur enters surface waters. However, biogeochemical cycling of sulfur in soils and the river corridor mediates sulfate exports. In this study, we identified the major forms of sulfur across multiple compartments and scales in a pristine mountainous watershed, including shale bedrock weathering profiles, hillslope soils, and alluvial floodplain sediments, in order to provide insight into biogeochemical sulfur cycling in a hydrologically variable alpine system. X-ray absorption near-edge spectroscopy (XANES) analysis of shale weathering profiles showed clear evidence of pyrite oxidation to sulfate, with large accumulations of intermediate S(0) (20%–53%). Micro-scale XANES showed evidence of reprecipitation of pyrite at fracture surfaces within the permanently saturated zone. Organic sulfur dominated S speciation in shallow hillslope soil and floodplain sediment, with little evidence of reduced inorganic S. However, mackinawite formation, representing active sulfate reduction, was observed in saturated oxbow sediments and saturated weathered shale underlying floodplain sediments. Further evidence of sulfate reduction from aqueous sulfur isotopic analysis was observed in shallow groundwater transects across an Fe-reducing meander, whereas increases in pore water sulfate concentrations implied sulfur oxidation at other locations. Furthermore, the data present an integrated picture of sulfur cycling in a shale-dominated watershed, where riverine sulfate exports are mediated by biological cycling, particularly in redox-stratified and temporally dynamic hyporheic zone sediments.

58 GEOSCIENCES↗

DOE OSTI2022

Disinfection byproducts formed during drinking water treatment reveal an export control point for dissolved organic matter in a subalpine headwater stream

Changes in climate, season, and vegetation can alter organic export from watersheds. While an accepted tradeoff to protect public health, disinfection processes during drinking water treatment can adversely react with organic compounds to form disinfection byproducts (DBPs). By extension, DBP monitoring can yield insights into hydrobiogeochemical dynamics within watersheds and their implications for water resource management. In this study, we analyzed temporal trends from a water treatment facility that sources water from Coal Creek in Crested Butte, Colorado. These trends revealed a long-term increase in haloacetic acid and trihalomethane formation over the period of 2005-2020. Additionally, disproportionate export of dissolved organic carbon and formation of DBPs that exceeded regulatory guidelines were consistently recorded in association with late spring freshet. Synoptic sampling of the creek in 2020 and 2021 identified a biogeochemical hotspot for organic carbon export in the upper domain of the watershed that contained a prominent fulvic acid-like fluorescent signature. DBP formation potential analyses from this domain yielded similar ratios of dominant DBP species to those formed at the drinking water facility. Spectrometric qualitative analyses of pre and post-reacted waters with hypochlorite indicated ligninlike and condensed hydrocarbon-like molecules were the major reactive chemical classes during chlorine-based disinfection. This study demonstrates how drinking water quality archives can be coupled with synoptic sampling to identify and understand export control points for dissolved organic matter. Resultant understanding can be applied in water treatment and watershed management plans to improve drinking water quality sourced from analogous watersheds.

47 OTHER INSTRUMENTATION↗