Search for: All records

Abstract Understanding controls on solute export to streams is challenging because heterogeneous catchments can respond uniquely to drivers of environmental change. To understand general solute export patterns, we used a large‐scale inductive approach to evaluate concentration–discharge (C–Q) metrics across catchments spanning a broad range of catchment attributes and hydroclimatic drivers. We leveraged paired C–Q data for 11 solutes from CAMELS‐Chem, a database built upon an existing dataset of catchment and hydroclimatic attributes from relatively undisturbed catchments across the contiguous USA. Because C–Q relationships with Q thresholds reflect a shift in solute export dynamics and are poorly characterized across solutes and diverse catchments, we analysed C–Q relationships using Bayesian segmented regression to quantify Q thresholds in the C–Q relationship. Threshold responses were rare, representing only 12% of C–Q relationships, 56% of which occurred for solutes predominantly sourced from bedrock. Further, solutes were dominated by one or two C–Q patterns that reflected vertical solute–source distributions. Specifically, solutes predominantly sourced from bedrock had diluting C–Q responses in 43%–70% of catchments, and solutes predominantly sourced from soils had more enrichment responses in 35%–51% of catchments. We also linked C–Q relationships to catchment and hydroclimatic attributes to understand controls on export patterns. The relationships were generally weak despite the diversity of solutes and attribute types considered. However, catchment and hydroclimatic attributes in the central USA typically drove the most divergent export behaviour for solutes. Further, we illustrate how our inductive approach generated new hypotheses that can be tested at discrete, representative catchments using deductive approaches to better understand the processes underlying solute export patterns. Finally, given these long‐term C–Q relationships are from minimally disturbed catchments, our findings can be used as benchmarks for change in more disturbed catchments. 

Abstract. Large sample datasets are transforming the catchment sciences, but there are few off-the-shelf stream water chemistry datasets with complementary atmospheric deposition, streamflow, meteorology, and catchment physiographic attributes. The existing CAMELS (Catchment Attributes and Meteorology for Large-sample Studies) dataset includes data on topography, climate, streamflow, land cover, soil, and geology across the continental US. With CAMELS-Chem, we pair these existing attribute data for 516 catchments with atmospheric deposition data from the National Atmospheric Deposition Program and water chemistry and instantaneous discharge data from the US Geological Survey over the period from 1980 through 2018 in a relational database and corresponding dataset. The data include 18 common stream water chemistry constituents: Al, Ca, Cl, dissolved organic carbon, total organic carbon, HCO3, K, Mg, Na, total dissolved N, total organic N, NO3, dissolved oxygen, pH (field and lab), Si, SO4, and water temperature. Annual deposition loads and concentrations include hydrogen, NH4, NO3, total inorganic N, Cl, SO4, Ca, K, Mg, and Na. We demonstrate that CAMELS-Chem water chemistry data are sampled effectively across climates, seasons, and discharges for trend analysis and highlight the coincident sampling of stream constituents for process-based understanding. To motivate their use by the larger scientific community across a variety of disciplines, we show examples of how these publicly available datasets can be applied to trend detection and attribution, biogeochemical process understanding, and new hypothesis generation via data-driven techniques. 

Key Points We re‐evaluate equations proposed by Francis Hall to assess concentration‐discharge ( C ‐ Q ) relationships using newly available long‐term and high‐frequency data sets Across time steps we find that log‐log and log‐linear models perform equally well to describe C ‐ Q relationships Parametrization of storage‐discharge relationships via recession analyses provides additional insight to C ‐ Q relationships 

The concurrent reduction in acid deposition and increase in precipitation impact stream solute dynamics in complex ways that make predictions of future water quality difficult. To understand how changes in acid deposition and precipitation have influenced dissolved organic carbon (DOC) and nitrogen (N) loading to streams, we investigated trends from 1991 to 2018 in stream concentrations (DOC, ~3,800 measurements), dissolved organic nitrogen (DON, ~1,160 measurements), and dissolved inorganic N (DIN, ~2,130 measurements) in a forested watershed in Vermont, USA. Our analysis included concentration-discharge (C-Q) relationships and Seasonal Mann-Kendall tests on long-term, flow-adjusted concentrations. To understand whether hydrologic flushing and changes in acid deposition influenced long-term patterns by liberating DOC and dissolved N from watershed soils, we measured their concentrations in the leachate of 108 topsoil cores of 5 cm diameter that we flushed with solutions simulating high and low acid deposition during four different seasons. Our results indicate that DOC and DON often co-varied in both the long-term stream dataset and the soil core experiment. Additionally, leachate from winter soil cores produced especially high concentrations of all three solutes. This seasonal signal was consistent with C-Q relation showing that organic materials (e.g., DOC and DON), which accumulate during winter, are flushed into streams during spring snowmelt. Acid deposition had opposite effects on DOC and DON compared to DIN in the soil core experiment. Low acid deposition solutions, which mimic present day precipitation, produced the highest DOC and DON leachate concentrations. Conversely, high acid deposition solutions generally produced the highest DIN leachate concentrations. These results are consistent with the increasing trend in stream DOC concentrations and generally decreasing trend in stream DIN we observed in the long-term data. These results suggest that the impact of acid deposition on the liberation of soil carbon (C) and N differed for DOC and DON vs. DIN, and these impacts were reflected in long-term stream chemistry patterns. As watersheds continue to recover from acid deposition, stream C:N ratios will likely continue to increase, with important consequences for stream metabolism and biogeochemical processes. 

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. 

Abstract Oxygen (O₂) regulates soil reduction‐oxidation processes and therefore modulates biogeochemical cycles. The difficulties associated with accurately characterizing soil O₂variability have prompted the use of soil moisture as a proxy for O₂, as O₂diffusion into soil water is much slower than in soil air. The use of soil moisture alone as a proxy measurement for O₂could result in inaccurate O₂estimations. For example, O₂may remain high during cool months when soil respiration rates are low. We analyzed high‐frequency sensor data (e.g., soil moisture, CO₂, gas‐phase soil pore O₂) with a machine learning technique, the Self‐Organizing Map, to pinpoint suites of soil conditions associated with contrasting O₂regimes. At two riparian sites in northern Vermont, we found that O₂levels varied seasonally, and with soil moisture. For example, 47% of low O₂levels were associated with wet and cool soil conditions, whereas 32% were associated with dry and warm conditions. Contrastingly, the majority (62%) of high O₂conditions occurred under dry and warm conditions. High soil moisture levels did not always lead to low O₂, as 38% of high O₂values occurred under wet and cool conditions. Our results highlight challenges with predicting soil O₂solely based on water content, as variable combinations of soil and hydrologic conditions can complicate the relationship between water content and O₂. This indicates that process‐based ecosystem and denitrification models that rely solely on soil moisture to estimate O₂may need to incorporate other site and climate‐specific drivers to accurately predict soil O₂. 

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. 

Abstract The shallow and deep hypothesis suggests that stream concentration‐discharge (CQ) relationships are shaped by distinct source waters from different depths. Under this hypothesis, baseflows are typically dominated by groundwater and mostly reflect groundwater chemistry, whereas high flows are typically dominated by shallow soil water and mostly reflect soil water chemistry. Aspects of this hypothesis draw on applications like end member mixing analyses and hydrograph separation, yet direct data support for the hypothesis remains scarce. This work tests the shallow and deep hypothesis using co‐located measurements of soil water, groundwater, and streamwater chemistry at two intensively monitored sites, the W‐9 catchment at Sleepers River (Vermont, United States) and the Hafren catchment at Plynlimon (Wales). At both sites, depth profiles of subsurface water chemistry and stream CQ relationships for the 10 solutes analyzed are broadly consistent with the hypothesis. Solutes that are more abundant at depth (e.g., calcium) exhibit dilution patterns (concentration decreases with increasing discharge). Conversely, solutes enriched in shallow soils (e.g., nitrate) generally exhibit flushing patterns (concentration increases with increasing discharge). The hypothesis may hold broadly true for catchments that share such biogeochemical stratifications in the subsurface. Soil water and groundwater chemistries were estimated from high‐ and low‐flow stream chemistries with average relative errors ranging from 24% to 82%. This indicates that streams mirror subsurface waters: stream chemistry can be used to infer scarcely measured subsurface water chemistry, especially where there are distinct shallow and deep end members.