Research at long‐term catchment monitoring sites has generated a great volume, variety, and velocity of data for analysis of stream water chemistry dynamics. To harness the potential of these big data and extract patterns that are indicative of underlying functional relationships, machine learning tools have advantages over traditional statistical methods, and are increasingly being applied for dimension reduction, feature extraction, and trend identification. Still, as examples of complex systems, catchments are characterized by multivariate factor interactions and equifinality that are not easily identified by most machine‐learning methods. Using dissolved organic carbon (DOC) dynamics as an illustration, we applied a new evolutionary algorithm (EA) to extract geologic, topographic, meteorologic, hydrologic, and land use attributes that were correlated to mean stream DOC concentration in forested catchments distributed across the continental United States. The EA reduced dimensionality of our attribute dataset to identify the combination of factors, and their specific value ranges, that interacted to drive membership in High or Low mean DOC clusters. High mean DOC concentrations were associated with two distinct geographic locations of variable climatic and vegetative conditions, indicating equifinality. Our findings underscore the importance of critical zone structure in mediating hydrological and biogeochemical processes to govern DOC dynamics at the catchment scale. This multi‐scale, pattern‐to‐process approach is being applied to refine hypotheses for process‐based modeling of DOC dynamics in forested headwater streams at catchment to site scales.
more » « less- Award ID(s):
- 2012123
- PAR ID:
- 10419716
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 59
- Issue:
- 3
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
null (Ed.)Understanding and predicting catchment responses to a regional disturbance is difficult because catchments are spatially heterogeneous systems that exhibit unique moderating characteristics. Changes in precipitation composition in the Northeastern U.S. is one prominent example, where reduction in wet and dry deposition is hypothesized to have caused increased dissolved organic carbon (DOC) export from many northern hemisphere forested catchments; however, findings from different locations contradict each other. Using shifts in acid deposition as a test case, we illustrate an iterative “process and pattern” approach to investigate the role of catchment characteristics in modulating the steam DOC response. We use a novel dataset that integrates regional and catchment-scale atmospheric deposition data, catchment characteristics and co-located stream Q and stream chemistry data. We use these data to investigate opportunities and limitations of a pattern-to-process approach where we explore regional patterns of reduced acid deposition, catchment characteristics and stream DOC response and specific soil processes at select locations. For pattern investigation, we quantify long-term trends of flow-adjusted DOC concentrations in stream water, along with wet deposition trends in sulfate, for USGS headwater catchments using Seasonal Kendall tests and then compare trend results to catchment attributes. Our investigation of climatic, topographic, and hydrologic catchment attributes vs. directionality of DOC trends suggests soil depth and catchment connectivity as possible modulating factors for DOC concentrations. This informed our process-to-pattern investigation, in which we experimentally simulated increased and decreased acid deposition on soil cores from catchments of contrasting long-term DOC response [Sleepers River Research Watershed (SRRW) for long-term increases in DOC and the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO) for long-term decreases in DOC]. SRRW soils generally released more DOC than SSHCZO soils and losses into recovery solutions were higher. Scanning electron microscope imaging indicates a significant DOC contribution from destabilizing soil aggregates mostly from hydrologically disconnected landscape positions. Results from this work illustrate the value of an iterative process and pattern approach to understand catchment-scale response to regional disturbance and suggest opportunities for further investigations.more » « less
-
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.
-
Temporal patterns in stream chemistry provide integrated signals describing the hydrological and ecological state of whole catchments. However, stream chemistry integrates multi-scale signals of processes occurring in both the catchment and stream. Deconvoluting these signals could identify mechanisms of solute transport and transformation and provide a basis for monitoring ecosystem change. We applied trend analysis, wavelet decomposition, multivariate autoregressive state-space modeling, and analysis of concentration–discharge relationships to assess temporal patterns in high-frequency (15 min) stream chemistry from permafrost-influenced boreal catchments in Interior Alaska at diel, storm, and seasonal time scales. We compared catchments that varied in spatial extent of permafrost to identify characteristic biogeochemical signals. Catchments with higher spatial extents of permafrost were characterized by increasing nitrate concentration through the thaw season, an abrupt increase in nitrate and fluorescent dissolved organic matter (fDOM) and declining conductivity in late summer, and flushing of nitrate and fDOM during summer rainstorms. In contrast, these patterns were absent, of lower magnitude, or reversed in catchments with lower permafrost extent. Solute dynamics revealed a positive influence of permafrost on fDOM export and the role of shallow, seasonally dynamic flowpaths in delivering solutes from high-permafrost catchments to streams. Lower spatial extent of permafrost resulted in static delivery of nitrate and limited transport of fDOM to streams. Shifts in concentration–discharge relationships and seasonal trends in stream chemistry toward less temporally dynamic patterns might therefore indicate reorganized catchment hydrology and biogeochemistry due to permafrost thaw.more » « less
-
Abstract Stream and shallow groundwater responses to rainfall are characterized by high spatial variability, but hydrologic response variability across small, agro‐forested sub‐catchments remains poorly understood. Conceivably, improved understanding in this regard will result in agricultural practices that more effectively limit nutrient runoff, erosion, and pollutant transport. Terrestrial hydrologic response approaches can provide valuable information on stream‐aquifer connectivity in these mixed‐use watersheds. A study was implemented, including eight stream and co‐located shallow groundwater monitoring sites, in a small sub‐catchment of the Chesapeake Bay watershed in the Northeast, USA to advance this ongoing need. During the study period, 100 precipitation‐receiving days (i.e., 24‐hour periods, midnight to midnight) were observed. On average, the groundwater table responded more to precipitation than stream stage (level change of 0.03 vs. 0.01 m and rainfall‐normalized level change estimate of 3.81 vs. 3.37). Median stream stage responses, groundwater table responses, and response ratios were significantly different between sub‐catchments (
n = 8;p < 0.001). Study area average precipitation thresholds for runoff and shallow groundwater flow were 2.8 and 0.6 cm, respectively. Individual sub‐catchment thresholds ranged from 0.5 to 2.8 cm for runoff and 0.2 to 1.3 cm for shallow groundwater flow. Normalized response lag times between the stream and shallow groundwater ranged from −0.50 to 3.90 s·cm−1, indicating that stormflow in one stream section was regulated by groundwater flow during the period of study. The observed differences in hydrologic responses to precipitation advance future modelling efforts by providing examples of how terrestrial groundwater response methods can be used to investigate sub‐catchment spatial variability in stream‐aquifer gradients with co‐located shallow groundwater and stream stage data. Additionally, results demonstrate asynchronous stream and shallow groundwater responses on precipitation‐receiving days, which may hold important implications for modelling hydrologic and biogeochemical fate and transport processes in small, agro‐forested catchments. -
Abstract Tracer‐aided rainfall‐runoff modelling is a promising tool for understanding catchment hydrology, particularly when tracers provide information about coupled hydrological‐biogeochemical processes. Such models allow for predicting the quality and quantity of water under changing climatic and anthropogenic conditions. Here, we present the Spatially‐distributed Tracer‐Aided Rainfall‐Runoff model with a coupled biogeochemical reactive tracer module (STARR‐DOC) to simulate dissolved organic carbon (DOC) dynamics and sources. The STARR‐DOC model was developed and tested for a humid high Andean ecosystem (páramo) using high‐resolution hourly DOC and hydrometeorological data to simulate hourly discharge and DOC at a fine spatial (10 × 10 m) resolution. Overall, the model was able to acceptably reproduce discharge (KGE ~ 0.45) and stream DOC (KGE ~ 0.69) dynamics. Spatially distributed DOC simulations were independently compared using point DOC measurements for different soil types across the catchment, which allowed for identifying DOC production hot spots and hot moments. Results showed higher hydrological connectivity between slopes and valleys with increasing precipitation. Wetter conditions also favoured DOC production (wet month = 82 mg L−1, dry month = 5 mg L−1) and transport to the stream network (DOC concentrations: during events ~15 mg L−1, during baseflows ~4 mg L−1). Our results also suggest that minor changes in meteorological conditions directly affect páramo soil water dynamics and biogeochemistry. Knowledge of when and where DOC production in mountain catchments is greatest is important for water managers to understand when they make decisions about water security, especially considering climate change predictions for the Andean region.