Distinct Source Water Chemistry Shapes Contrasting Concentration‐Discharge Patterns. Issue 5 (24th May 2019)

• Record Type: Journal Article
• Authors: Zhi, Wei; Li, Li; Dong, Wenming; Brown, Wendy; Kaye, Jason; Steefel, Carl; Williams, Kenneth H.
• Is Part Of: Water resources research. Volume 55:Issue 5(2019)
• Page Start: 4233
• Page End: 4251
• Publication Date: 2019-05-24
• Subjects: Hydrology -- Periodicals; 333.91
• DOI: 10.1029/2018WR024257
• Languages: English
• ISSNs: 0043-1397

Abstract

Abstract: Understanding concentration‐discharge (C‐Q) relationships are essential for predicting chemical weathering and biogeochemical cycling under changing climate and anthropogenic conditions. Contrasting C‐Q relationships have been observed widely, yet a mechanistic framework that can interpret diverse patterns remains elusive. This work hypothesizes that seemingly disparate C‐Q patterns are driven by switching dominance of end‐member source waters and their chemical contrasts arising from subsurface biogeochemical heterogeneity. We use data from Coal Creek, a high‐elevation mountainous catchment in Colorado, and a recently developed watershed reactive transport model (BioRT‐Flux‐PIHM). Sensitivity analysis and Monte‐Carlo simulations (500 cases) show that reaction kinetics and thermodynamics and distribution of source materials across depths govern the chemistry gradients of shallow soil water and deeper groundwater entering the stream. The alternating dominance of organic‐poor yet geo‐solute‐rich groundwater under dry conditions and organic‐rich yet geo‐solute‐poor soil water during spring melt leads to the flushing pattern of dissolved organic carbon and the dilution pattern of geogenic solutes (e.g., Na, Ca, and Mg). In addition, the extent of concentration contrasts regulates the power law slopes ( b ) of C‐Q patterns via a general equation b = δ b C ratio C ratio, 1 / 2 + C ratio + b min . At low ratios of soil water versus groundwater concentrations (Cratio = Csw /Cgw < 0.6), dilution occurs; at high ratios (Cratio > 1.8), flushing arises; chemostasis occurs in between. This equation quantitatively interprets b values of 11 solutes (dissolved organic carbon, dissolved P, NO3 −, K, Si, Ca, Mg, Na, Al, Mn, and Fe) from three catchments (Coal Creek, Shale Hills, and Plynlimon) of differing climate, geologic, and land cover conditions. This indicates potentially broad regulation of subsurface biogeochemical heterogeneity in determining C‐Q patterns and wide applications of this equation in quantifying b values, which can have broad implications for predicting chemical weathering and biogeochemical transformation at the watershed scale. Key Points: Contrasting C‐Q patterns are shaped by alternating dominance of soil water and groundwater with distinct water chemistries Dilution occurs for geogenic solutes concentrated in reactive rocks at depth, while flushing arises for biogenic solutes enriched in shallow soils A relationship between C‐Q slope and soil‐groundwater concentration ratios is established, explaining C‐Q patterns of 11 solutes from three catchments