skip to main content


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
NSF-PAR ID:
10419717
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

    Uncertainty in the estimation of hydrologic export of solutes has never been fully evaluated at the scale of a small‐watershed ecosystem. We used data from the Gomadansan Experimental Forest, Japan, Hubbard Brook Experimental Forest, USA, and Coweeta Hydrologic Laboratory, USA, to evaluate many sources of uncertainty, including the precision and accuracy of measurements, selection of models, and spatial and temporal variation. Uncertainty in the analysis of stream chemistry samples was generally small but could be large in relative terms for solutes near detection limits, as is common for ammonium and phosphate in forested catchments. Instantaneous flow deviated from the theoretical curve relating height to discharge by up to 10% at Hubbard Brook, but the resulting corrections to the theoretical curve generally amounted to <0.5% of annual flows. Calibrations were limited to low flows; uncertainties at high flows were not evaluated because of the difficulties in performing calibrations during events. However, high flows likely contribute more uncertainty to annual flows because of the greater volume of water that is exported during these events. Uncertainty in catchment area was as much as 5%, based on a comparison of digital elevation maps with ground surveys. Three different interpolation methods are used at the three sites to combine periodic chemistry samples with streamflow to calculate fluxes. The three methods differed by <5% in annual export calculations for calcium, but up to 12% for nitrate exports, when applied to a stream at Hubbard Brook for 1997–2008; nitrate has higher weekly variation at this site. Natural variation was larger than most other sources of uncertainty. Specifically, coefficients of variation across streams or across years, within site, for runoff and weighted annual concentrations of calcium, magnesium, potassium, sodium, sulphate, chloride, and silicate ranged from 5 to 50% and were even higher for nitrate. Uncertainty analysis can be used to guide efforts to improve confidence in estimated stream fluxes and also to optimize design of monitoring programmes. © 2014 The Authors.Hydrological Processespublished John Wiley & Sons, Ltd.

     
    more » « less
  3. 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.

     
    more » « less
  4. 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.

     
    more » « less
  5. 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