skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Machine‐Learning Reveals Equifinality in Drivers of Stream DOC Concentration at Continental Scales
Abstract 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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  
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
  1. 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
  2. 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 (CQ) 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. 
    more » « less
  3. 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
  4. Abstract Climate change is thawing and potentially mobilizing vast quantities of organic carbon (OC) previously stored for millennia in permafrost soils of northern circumpolar landscapes. Climate‐driven increases in fire and thermokarst may play a key role in OC mobilization by thawing permafrost and promoting transport of OC. Yet, the extent of OC mobilization and mechanisms controlling terrestrial‐aquatic transfer are unclear. We demonstrate that hydrologic transport of soil dissolved OC (DOC) from the active layer and thawing permafrost to headwater streams is extremely heterogeneous and regulated by the interactions of soils, seasonal thaw, fire, and thermokarst. Repeated sampling of streams in eight headwater catchments of interior Alaska showed that the mean age of DOC for each stream ranges widely from modern to ∼2,000 years B.P. Together, an endmember mixing model and nonlinear, generalized additive models demonstrated that Δ14C‐DOC signature (and mean age) increased from spring to fall, and was proportional to hydrologic contributions from a solute‐rich water source, related to presumed deeper flow paths found predominantly in silty catchments. This relationship was correlated with and mediated by catchment properties. Mean DOC ages were older in catchments with >50% burned area, indicating that fire is also an important explanatory variable. These observations underscore the high heterogeneity in aged C export and difficulty of extrapolating estimates of permafrost‐derived DOC export from watersheds to larger scales. Our results provide the foundation for developing a conceptual model of permafrost DOC export necessary for advancing understanding and prediction of land‐water C exchange in changing boreal landscapes. 
    more » « less
  5. Abstract Headwater catchments have strong impacts on downstream waterways, near‐shore ecosystems, and the quality of water available for growing human populations. Thus, understanding how water and solutes are exported through these upland landscapes is critically important. A growing body of literature highlights the interaction of topography, climate, and the critical zone structure as a key control on streamflow and chemical export. However, more focused work is needed to pinpoint how variability in subsurface structure across lithologically complex regions impacts streamflow and chemical signals at catchment outlets. Here, we aim to better understand how lithology and subsurface critical zones modulate streamflow response and solute export patterns in two central coastal California headwater catchments that are similar in topography, vegetation, and climate but have different lithologies. We monitored streamflow and collected surface water samples at the catchment outlets for dissolved major ions and organic carbon (DOC) for two consecutive water years. The catchment with mélange bedrock displayed much flashier hydrologic behavior with 7.8 times higher peak flow values and 1.9 times higher mean event concentrations of DOC, suggesting shorter and shallower hydrologic flow paths that likely arise from regions of shallower bedrock. Despite distinct hydrologic behavior and DOC export, dissolved major ion concentrations were broadly similar and chemostatic, which may be driven by rapid chemical reactions in the critical zone of both catchments. Our work contributes to building an integrated understanding of how subtle differences in catchment structure can have profound impacts on how water and solutes are routed through headwater catchments. 
    more » « less