Critical Zone (CZ) Export to Streams as Indicator for CZ Structure and Function

Over the past two decades, headwater streams of the northern hemisphere have shown increased amounts of dissolved organic carbon (DOC), coinciding with decreased acid deposition. The exact nature of the mechanistic link between precipitation composition and stream water DOC is still widely debated in the literature. We hypothesize that soil aggregates are the main source of stream water DOC and that DOC release is greater in organic rich, riparian soils vs. hillslope soils. To test these hypotheses, we collected soils from two main landscape positions (hillslope and riparian zones) from the acid-impacted Sleepers River Research Watershed in northeastern Vermont. We performed aqueous soil extracts with solutions of different ionic strength (IS) and composition to simulate changes in soil solution. We monitored dynamic changes in soil particle size, aggregate architecture and composition, leachate DOC concentrations, dissolved organic matter (DOM) characteristics by fluorescence spectroscopy and trends in bioavailability. In low IS solutions, extractable DOC concentrations were significantly higher, particle size (by laser diffraction) was significantly smaller and organic material was separated from mineral particles in scanning electron microscope observations. Furthermore, higher DOC concentrations were found in Na⁺ compared to Ca²⁺ solutions of the same IS. These effects are attributed to aggregate dispersion due to expanding diffuse double layers in decreased IS solutions and to decreased bridging by divalent cations. Landscape position impacted quality but not quantity of released DOC. Overall, these results indicate that soil aggregates might be one important link between Critical Zone inputs (i.e., precipitation) and exports in streams.

Concentration-discharge relationships are a key tool for understanding the sources and transport of material from watersheds to fluvial networks. Storm events in particular provide insight into variability in the sources of solutes and sediment within watersheds, and the hydrologic pathways that connect hillslope to stream channel. Here we examine high-frequency sensor-based specific conductance and turbidity data from multiple storm events across two watersheds (Quebrada Sonadora and Rio Icacos) with different lithology in the Luquillo Mountains of Puerto Rico, a forested tropical ecosystem. Our analyses include Hurricane Maria, a category 5 hurricane. To analyze hysteresis, we used a recently developed set of metrics to describe and quantify storm events including the hysteresis index (HI), which describes the directionality of hysteresis loops, and the flushing index (FI), which can be used to infer whether the mobilization of material is source or transport limited. We also examine the role of antecedent discharge to predict hysteretic behavior during storms. Overall, specific conductance and turbidity showed contrasting responses to storms. The hysteretic behavior of specific conductance was similar across sites, displaying clockwise hysteresis and a negative FI indicating proximal sources of solutes and consistent source limitation. In contrast, the directionality of turbidity hysteresis was significantly different between watersheds, although both had strong flushing behavior indicative of transport limitation. Overall, models that included antecedent discharge did not perform any better than models with peak discharge alone, suggesting that the magnitude and trajectory of an individual event was the strongest driver of material flux and hysteretic behavior. Hurricane Maria produced unique hysteresis metrics within both watersheds, indicating a distinctive response to this major hydrological event. The similarity in response of specific conductance to storms suggests that solute sources and pathways are similar in the two watersheds. The divergence in behavior for turbidity suggests that sources and pathways of particulate matter vary between the two watersheds. The use of high-frequency sensor data allows the quantification of storm events while index-based metrics of hysteresis allow for the direct comparison of complex storm events across a heterogeneous landscape and variable flow conditions.