skip to main content


Title: Regional Drivers of Stream Chemical Behavior: Leveraging Lithology, Land Use, and Climate Gradients Across the Colorado River, Texas USA
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.

 
more » « less
Award ID(s):
1933259 1933261
NSF-PAR ID:
10382378
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
58
Issue:
11
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Quantifying the routing of snowmelt to surface water is critical for predicting the impacts of atmospheric deposition and changing land use on water quality in montane catchments. To investigate solute sources and streamflow in the montane Provo River watershed (Utah, USA), we used time‐series87Sr/86Sr ratios sampled at three sites (Soapstone, Woodland and Hailstone) across a gradient of bedrock types. Soils are influenced by aeolian dust contributions, with distinct87Sr/86Sr ratios relative to siliciclastic bedrock, providing an opportunity to investigate shallow versus deeper flow paths for controlling water chemistry. At the most upstream site (Soapstone), Sr concentrations averaged ~17 μg/L with minimal dilution during snowmelt suggesting subsurface flow paths dominated streamflow. However, a decrease in87Sr/86Sr ratios from ~0.717 during baseflow to as low as ~0.713 during snowmelt indicated the activation of shallow flow paths through dust‐derived soils. In contrast, downstream sites receiving water inputs from Sr‐rich carbonate bedrock (Woodland and Hailstone) exhibited strong dilution of Sr from ~120 to 20 μg/L and an increase in87Sr/86Sr ratios from ~0.7095 to ~0.712 during snowmelt. A three‐component mixing model using87Sr/86Sr ratios and Sr concentrations at Soapstone showed water inputs were dominated by direct snowmelt and flushed soil water during runoff and groundwater during baseflow. At Woodland and Hailstone, a two‐component mixing model showed that the river was a mixture of groundwater and up to 75% upstream channel water during snowmelt. Our findings highlight the importance of flushed soil water for controlling stream water discharge and chemistry during snowmelt, with the signal from the upstream site propagating downstream in a nested catchment. Further, aeolian dust contributes to the solute chemistry of montane streams with potential impacts on water quality along shallow flow paths. Potential contaminants in these surface soils (e.g., Pb deposition in dust) may have significant impacts on water quality during snowmelt runoff.

     
    more » « less
  2. Meteoric waters move along pathways in the subsurface that differ as a function of lithology because of the effects of chemical and physical weathering. To explore how this affects stream chemistry, we investigated watersheds around an igneous intrusion in the Luquillo Mountains (Puerto Rico). We analyzed streams on 1) unmetamorphosed country rock (volcaniclastic sedimentary strata, VC) surrounding an igneous intrusion, 2) the quartz-diorite intrusion (QD), and 3) the metamorphosed aureole rock (hornfels-facies volcaniclastics, HF). These lithologies differ physically and chemically but weather under the same tropical rain forest conditions. The sedimentary VC lithology is pervasively fractured while the massive QD and HF lithologies are relatively unfractured. However, the QD fractures during weathering to produce spheroidally-weathered corestones surrounded by cm-thick rindlets of increasingly weathered rock. Meteoric waters flow pervasively through the network of already-fractured VC rock and the spheroidally weathered rindlets on the QD, but only access a limited fraction of the HF, explaining why streams draining HF are the most dilute in the mountains. This results in various thicknesses of regolith from thick (VC) to moderate (QD) to thin or nonexistent (HF). The pervasive fractures allow groundwater to flow deeply through the VC and then return to the mainstem river (Río Mameyes) at lower elevations. These “rock waters” drive concentrations of rock-derived solutes (silica, base cations, sulfate, phosphate) higher in the lower reaches of the stream. Water also flows through weathering-induced fractures on the QD at high elevations where rindletted corestones are present in stacks, and this water flux dissolves plagioclase and hornblende and oxidizes biotite. This “QD rock water” is not generated at lower elevations in the Río Icacos watershed, where stacks of corestones are absent, and contributions to stream solutes derive from weathering of feldspar- and hornblende-depleted saprolite. The stream chemistry in the QD-dominated watershed (Río Icacos) thus varies from concentrated QD-rock water at channel heads below steep ridgelines toward more diluted “saprolite water” downstream. These observations emphasize the importance of lithology and fracture patterns in dictating water flowpaths, stream chemistry, and regolith development in headwater catchments. 
    more » « less
  3. Biogeochemical properties of soils play a crucial role in soil and stream chemistry throughout a watershed. How water interacts with soils during subsurface flow can have impacts on water quality, thus, it is fundamental to understand where and how certain soil water chemical processes occur within a catchment. In this study, ~200 soil samples were evaluated throughout a small catchment in the Front Range of Colorado, USA to examine spatial and vertical patterns in major soil solutes among different landscape units: riparian areas, alluvial/colluvial fans, and steep hillslopes. Solutes were extracted from the soil samples in the laboratory and analyzed for major cation (Li, K, Mg, Br, and Ca) and anion (F, Cl, NO 2 , NO 3 , PO 4 , and SO 4 ) concentrations using ion chromatography. Concentrations of most solutes were greater in near surface soils (10 cm) than in deeper soils (100 cm) across all landscape units, except for F which increased with depth, suggestive of surface accumulation processes such as dust deposition or enrichment due to biotic cycling. Potassium had the highest variation between depths, ranging from 1.04 mg/l (100 cm) to 3.13 mg/l (10 cm) sampled from riparian landscape units. Nearly every solute was found to be enriched in riparian areas where vegetation was visibly denser, with higher mean concentrations than the hillslopes and fans, except for NO 3 which had higher concentrations in the fans. Br, NO 2 , and PO 4 concentrations were often below the detectable limit, and Li and Na were not variable between depths or landscape units. Ratioed stream water concentrations (K:Na, Ca:Mg, and NO 3 :Cl) vs. discharge relationships compared to the soil solute ratios indicated a hydraulic disconnection between the shallow soils (<100 cm) and the stream. Based on the comparisons among depths and landscape units, our findings suggest that K, Ca, F, and NO 3 solutes may serve as valuable tracers to identify subsurface flowpaths as they are distinct among landscape units and depth within this catchment. However, interflow and/or shallow groundwater flow likely have little direct connection to streamflow generation. 
    more » « less
  4. Abstract. Salinity is one of the most common water quality threats in riverbasins and irrigated regions worldwide. However, no available numericalmodels simulate all major processes affecting salt ion fate and transport at the watershed scale. This study presents a new salinity module for the SWAT model that simulates the fate and transport of eight major salt ions(SO42-, Ca2+, Mg2+, Na+, K+, Cl−,CO32-, HCO3-) in a watershed system. The module accountsfor salt transport in surface runoff, soil percolation, lateral flow,groundwater, and streams, and equilibrium chemistry reactions in soil layersand the aquifer. The module consists of several new subroutines that areimbedded within the SWAT modelling code and one input file containing soilsalinity and aquifer salinity data for the watershed. The model is appliedto a 732 km2 salinity-impaired irrigated region within the ArkansasRiver Valley in southeastern Colorado and tested against root zone soilsalinity, groundwater salt ion concentration, groundwater salt loadings tothe river network, and in-stream salt ion concentration. The model can be auseful tool in simulating baseline salinity transport and investigatingsalinity best management practices in watersheds of varying spatial scales. 
    more » « less
  5. In the Baltimore urban long-term ecological research (LTER) project, (Baltimore Ecosystem Study, BES) we use the watershed approach to evaluate integrated ecosystem function. The LTER research is centered on the Gwynns Falls watershed, a 17,150 ha catchment that traverses a gradient from the urban core of Baltimore, through older urban residential (1900 - 1950) and suburban (1950- 1980) zones, rapidly suburbanizing areas and a rural/suburban fringe. Our long-term sampling network includes four longitudinal sampling sites along the Gwynns Falls as well as several small (40 - 100 ha) watersheds located within or near to the Gwynns Falls. The longitudinal sites provide data on water and nutrient fluxes in the different land use zones of the watershed (rural/suburban, rapidly suburbanizing, old suburban, urban core) and the small watersheds provide more focused data on specific land use areas (forest, agriculture, rural/suburban, urban). Each of the gaging sites is continuously monitored for discharge and is sampled weekly for chemistry. Additional chemical sampling is carried out in a supplemental set of sites to provide a greater range of land use. Weekly analyses includes nitrate, phosphate, total nitrogen, total phosphorus, chloride and sulfate, total suspended solids, turbidity, fecal coliforms, temperature, dissolved oxygen and pH. Cations, dissolved organic carbon and nitrogen and metals are measured on selected samples. This dataset presents stream chemistry from the Upper Gwynns Falls tributaries. From April 1999 to August 2000 Johns Hopkins University graduate student Mark Colosimo sampled a group of sites in the Upper Gwynns Falls (Red Run, Horsehead Branch, Scotts Level Branch, Holly Branch). There were two sites in the Red Run drainage. This watershed drains approximately 19 km2 and has been rapidly suburbanizing since the early 1990s. Percent impervious surface was approximately 10% as of 2002. Sampling station Red Run 1 (RR1) was approximately 35 m upstream of the crossing of Painters Mill Bridge Road, and 350 m upstream of the confluence with the Gwynns Falls. Sampling station Red Run 2 (RR2) was farther upstream, between the Pleasant Hill and Dolfield road crossings. There were two sites along Scotts Level Branch, an older suburban watershed which was approximately 25% impervious surface in 1970. Site SL1 drains approximately 11 km2 and is located at the outlet of the sub-watershed, just above the confluence with Gwynns Falls. Site SL2 is at the McDonogh Rd. bridge crossing. The Horsehead Branch (HH) sampling site was located at the McDonogh Road crossing. It drains approximately 5 km2 that has undergone rapid urbanization since the mid 1980s. As of 1997 percent impervious surface was approximately 12%. The Holly Bank (HB) sampling site was located just upstream of Gwynnbrook Ave. Seventy percent of land in this drainage is classified residential. The Gwynns Falls at McDonogh (GF5) site was located at the McDonogh school / McDonogh road crossing of the Gwynns Falls and samples a drainage area of approximately 51 km2, with approximately 20% impervious surface. 
    more » « less