Abstract. Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of three shale-underlain headwater catchments located in Pennsylvania, USA (the forested Shale Hills Critical Zone Observatory), and Wales, UK (the peatland-dominated Upper Hafren and forest-dominated Upper Hore catchments in the Plynlimon forest), dissimilar concentration–discharge (C–Q) behaviors are best explained by contrasting landscape distributions of soil solution chemistry – especially dissolved organic carbon (DOC) – that have been established by patterns of vegetation and soil organic matter (SOM). Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Furthermore, concentration–discharge relationships of non-chemostatic solutes changed following tree harvest in the Upper Hore catchment in Plynlimon, while no changes were observed for chemostatic solutes, underscoring the role of vegetation in regulating the concentrations of certain elements in the stream. These results indicate that differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where SOM is dominantly in lowlands (e.g., Shale Hills), we infer that non-chemostatic elements associated with organic matter are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), these non-chemostatic elements are released later during rainfall events. The distribution of SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.
Snowmelt‐fed springs and small (0.5 km2) upland catchments in alpine areas of the western United States contribute significantly to the quantity and inorganic chemistry of water delivered to downstream basins but have not been studied extensively. Mineral weathering, transit time, and hydrologic mixing control the solute chemistry of waters that drain the upland zone of Niwot Ridge, Colorado Front Range, and adjacent areas in the granitic core of the Southern Rocky Mountains. Water in 37 springs sampled in this study flows in generally short steep paths (~0.3 km) through shallow regolith with mean transit times (MTT) of weeks to months, producing solutions dominated by Si, Ca2+, Na+, and HCO3−, locally SO42−. Rock type is a significant control on spring, surface, and shallow groundwater chemistry, and plagioclase (oligoclase) is the major source of dissolved Na+and Si. Concentrations of Ca2+exceed stoichiometric predictions of oligoclase weathering by ~3.5×; excess Ca2+likely represents weathering of aeolian material, vein calcite, or trace minerals. Concentrations of base cations and Si increase slowly with estimated MTT of 0.2 years for Niwot Ridge spring waters, and several years for shallow groundwater sampled by wells. Chemical weathering of silicate minerals is slow with estimated rates of ~2.0 and 0.2 pmol·m−2·s−1for oligoclase and microcline, respectively; the most mineralized spring waters are saturated only with respect to kaolinite and montmorillonite. More than 50% of the dissolved base cations + Si measured in Boulder Creek at Orodell (~25 km downstream) accumulate before water emerges from alpine springs on Niwot Ridge. Warming global temperatures are shifting more high‐elevation precipitation to rain, potentially changing run‐off patterns, transit time, and solute loads. Acquisition of solutes by alpine waters thus has implications far beyond small upland catchments.
more » « less- Award ID(s):
- 1637686
- NSF-PAR ID:
- 10461670
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 33
- Issue:
- 12
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 1756-1774
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of headwater catchments underlain by shale in Pennsylvania, USA (Shale Hills) and Wales, UK (Plynlimon), dissimilar concentration-discharge behaviors are best explained by contrasting landscape distributions of soil solution chemistry – especially dissolved organic carbon (DOC) – that have been established by patterns of vegetation. Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic "bioactive" behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where soil organic matter (SOM) is dominantly in lowlands (e.g., Shale Hills), bioactive elements are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), bioactive elements are released later during rainfall events. The distribution of vegetation and SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.more » « less
-
more » « less
Abstract Basaltic watersheds such as those found in Iceland are thought to be important sites of CO2sequestration via silicate weathering. However, determining the magnitude of CO2uptake depends on accurately interpreting river chemistry. Here, we compile geochemical data from Iceland and use them to constrain weathering processes. Specifically, we use a newly developed inverse model to quantify solute supply from rain and hydrothermal fluids as well as allow for variable silicate end‐member compositions, solutes to be removed via secondary phase formation, and some Ca to be sourced from carbonate dissolution. While some of these processes have been considered previously, they have not been considered together allowing us to newly determine their relative contributions. We find that weathering in Iceland is incongruent in two ways. First, solute release from primary silicates is characterized by a higher proportion of Na than would be expected from bulk basalts, which may reflect preferential weathering or some contribution from rhyolites. This Na enrichment is further enhanced by preferential Mg and K uptake by secondary phases. No samples in our data set (
n = 537) require carbonate dissolution even if isotopic data (δ 26Mg,δ 30Si,δ 44Ca, and/or87Sr/86Sr) are included. While some carbonate weathering is allowable, silicate weathering likely dominates. The complexity we observe in Iceland underscores the need for inverse models to account for a wide range of processes and end‐members. Given that riverine fluxes from Iceland are more Na‐rich than expected for congruent basalt weathering, the characteristic timescale of CO2drawdown is likely affected. -
more » « less
In weathered bedrock aquifers, groundwater is stored in pores and fractures that open as rocks are exhumed and minerals interact with meteoric fluids. Little is known about this storage because geochemical and geophysical observations are limited to pits, boreholes, or outcrops or to inferences based on indirect measurements between these sites. We trained a rock physics model to borehole observations in a well-constrained ridge and valley landscape and then interpreted spatial variations in seismic refraction velocities. We discovered that P-wave velocities track where a porosity-generating reaction initiates in shale in three boreholes across the landscape. Specifically, velocities of 2.7 ± 0.2 km/s correspond with growth of porosity from dissolution of chlorite, the most reactive of the abundant minerals in the shale. In addition, sonic velocities are consistent with the presence of gas bubbles beneath the water table under valley and ridge. We attribute this gas largely to CO2produced by 1) microbial respiration in soils as meteoric waters recharge into the subsurface and 2) the coupled carbonate dissolution and pyrite oxidation at depth in the rock. Bubbles may nucleate below the water table because waters depressurize as they flow from ridge to valley and because pores have dilated as the deep rock has been exhumed by erosion. Many of these observations are likely to also describe the weathering and flow path patterns in other headwater landscapes. Such combined geophysical and geochemical observations will help constrain models predicting flow, storage, and reaction of groundwater in bedrock systems.
-
more » « less
Abstract Understanding relationships between stream chemistry and watershed factors: land use/land cover, climate, and lithology are crucial to improving our knowledge of critical zone processes that influence water quality. We compiled major ion data from >100 monitoring stations collected over 60 years (1958–2018) across the Colorado River Watershed in Texas (103,000 km2). We paired this river chemistry data with complementary lithology, land use, climate, and stream discharge information. Machine learning techniques were used to produce new insights on controls of stream water chemical behavior, which were validated using traditional multivariate analyses. Studies on stream flow and chemistry in the American west and globally have shown strong relationships between major ion chemical composition, climate, and lithology which hold true for the Colorado River basin in this study. Reactive minerals, including carbonates and evaporites, dominate major ion chemistry across the upper, low‐precipitation regions of the watershed. Upstream and middle reaches of the Colorado River showed shifts from Na‐Cl‐SO4dominated water from multiple sources including dissolution of gypsum and halite in shallow groundwater, and agricultural activities, to Ca‐HCO3water types controlled by carbonate dissolution. In the lower portion of the watershed multiple analyses demonstrate that stream chemistry is more influenced by greater precipitation and the presence of silicate minerals than the middle and upstream reaches. This study demonstrates the power of applying machine learning approaches to publicly available long term water chemistry data sets to improve the understanding of watershed interactions with surficial lithology, salinity sources, and anthropogenic influences of water quality.
