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.

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. Disproportionate export of dissolved organic carbon and formation of DBPs that exceeded maximum contaminant levels 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 regulated DBP classes to those formed at the drinking water facility. Spectrometric qualitative analyses of pre and post-reacted waters with hypochlorite indicated lignin-like and condensed hydrocarbon-like molecules were the major reactive chemical classes during chlorine-based disinfection. This study 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 approach could be applied to identify and characterize analogous watersheds where seasonal or climate-associated organic matter export challenge water treatment disinfection and by extension inform watershed management and drinking water treatment.

Although bedrock weathering strongly influences water quality and global carbon and nitrogen budgets, the weathering depths and rates within subsurface are not well understood nor predictable. Determination of both porewater chemistry and subsurface water flow are needed in order to develop more complete understanding and obtain weathering rates. In a long-term field study, we applied a multiphase approach along a mountainous watershed hillslope transect underlain by marine shale. Here we report three findings. First, the deepest extent of the water table determines the weathering front, and the range of annually water table oscillations determines the thickness of the weathering zone. Below the lowest water table, permanently water-saturated bedrock remains reducing, preventing deeper pyrite oxidation. Secondly, carbonate minerals and potentially rock organic matter share the same weathering front depth with pyrite, contrary to models where weathering fronts are stratified. Thirdly, the measurements-based weathering rates from subsurface shale are high, amounting to base cation exports of about 70 kmol_c ha^-1 y^-1, yet consistent with weathering of marine shale. Finally, by integrating geochemical and hydrological data we present a new conceptual model that can be applied in other settings to predict weathering and water quality responses to climate change.

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that, at this site, bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification probably occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadose zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that, under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle.