skip to main content

Title: Spatiotemporal Drivers of Hydrochemical Variability in a Tropical Glacierized Watershed in the Andes

There is a critical knowledge gap about how glacier retreat in remote and rapidly warming tropical montane watersheds will impact solute export, which has implications for downstream geochemical cycling and ecological function. Because tropical glacierized watersheds are often uniquely characterized by year‐round ablation, upslope vegetation migration, and significant groundwater flow, baseline understanding is needed of how spatiotemporal variables within these watersheds control outlet hydrochemistry. We implemented a recently developed reactive transport watershed model, BioRT‐Flux‐PIHM, for a sub‐humid glacierized watershed in the Ecuadorian Andes with young volcanic soils and fractured bedrock. We found a unique simulated concentration and discharge (C‐Q) pattern that was mostly chemostatic but superimposed by dilution episodes. The chemostatic background was attributed to large simulated contributions of groundwater (subsurface lateral flow) to streamflow, of which a notable fraction (37%) comprised infiltrated ice‐melt. Relatively constant concentrations were further maintained in the model because times and locations of lower mineral surface wetting and dissolution were offset by concentrating effects of greater evapotranspiration. Ice‐melt did not all infiltrate in simulations, especially during large precipitation events, when high surface runoff contributions to discharge triggered dilution episodes. In a model scenario without ice‐melt, major ion concentrations, including Na+, Ca2+, and Mg2+, became more strongly chemostatic and higher, but weathering rates decreased, attenuating export by 23%. We expect this reduction to be exacerbated by higher evapotranspiration and drier conditions with expanded vegetation. This work brings to light the importance of subsurface meltwater flow, ecohydrological variability, and interactions between melt and precipitation for controlling hydrochemical processes in tropical watersheds with rapidly retreating glaciers.

more » « less
Award ID(s):
1759071 1758854
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Climate models predict amplified warming at high elevations in low latitudes,making tropical glacierized regions some of the most vulnerable hydrologicalsystems in the world. Observations reveal decreasing streamflow due toretreating glaciers in the Andes, which hold 99% of all tropicalglaciers. However, the timescales over which meltwater contributes tostreamflow and the pathways it takes – surface and subsurface – remainuncertain, hindering our ability to predict how shrinking glaciers willimpact water resources. Two major contributors to this uncertainty are thesparsity of hydrologic measurements in tropical glacierized watersheds andthe complication of hydrograph separation where there is year-round glaciermelt. We address these challenges using a multi-method approach that employsrepeat hydrochemical mixing model analysis, hydroclimatic time seriesanalysis, and integrated watershed modeling. Each of these approachesinterrogates distinct timescale relationships among meltwater, groundwater,and stream discharge. Our results challenge the commonly held conceptualmodel that glaciers buffer discharge variability. Instead, in a subhumidwatershed on Volcán Chimborazo, Ecuador, glacier melt drives nearly allthe variability in discharge (Pearson correlation coefficient of 0.89 insimulations), with glaciers contributing a broad range of 20%–60%or wider of discharge, mostly (86%) through surface runoff on hourlytimescales, but also through infiltration that increases annual groundwatercontributions by nearly 20%. We further found that rainfall may enhanceglacier melt contributions to discharge at timescales that complement glaciermelt production, possibly explaining why minimum discharge occurred at thestudy site during warm but dry El Niño conditions, which typicallyheighten melt in the Andes. Our findings caution against extrapolations fromisolated measurements: stream discharge and glacier melt contributions intropical glacierized systems can change substantially at hourly tointerannual timescales, due to climatic variability and surface to subsurfaceflow processes.

    more » « less
  2. Abstract

    As glaciers around the world rapidly lose mass, the tight coupling between glaciers and downstream ecosystems is resulting in widespread impacts on global hydrologic and biogeochemical cycling. However, a range of challenges make it difficult to conduct research in glacierized systems, and our knowledge of seasonally changing hydrologic processes and solute sources and signatures is limited. This in turn hampers our ability to make predictions on solute composition and flux. We conducted a broad water sampling campaign in order to understand the present‐day partitioning of water sources and associated solutes in Alaska's Wolverine Glacier watershed. We established a relationship between electrical conductivity and streamflow at the watershed outlet to divide the melt season into four hydroclimatic periods. Across hydroclimatic periods, we observed a shift in nonglacial source waters from snowmelt‐dominated overland and shallow subsurface flow paths to deeper groundwater flow paths. We also observed the shift from a low‐ to high‐efficiency subglacial drainage network and the associated flushing of water stored subglacially with higher solute loads. We used calcium, the dominant dissolved ion, from watershed outlet samples to estimate solute fluxes for each hydroclimatic period across two melt seasons. We found between 40% and 55% of Ca2+export occurred during the late season rainy period. This partitioning of the melt season coupled with a characterization of the chemical makeup and magnitude of solute export provides new insight into a rapidly changing watershed and creates a framework to quantify and predict changes to solute fluxes across a melt season.

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

    Threshold changes in rainfall‐runoff generation commonly represent shifts in runoff mechanisms and hydrologic connectivity controlling water and solute transport and transformation. In watersheds with limited human influence, threshold runoff responses reflect interaction between precipitation event and antecedent soil moisture. Similar analyses are lacking in intensively managed landscapes where installation of subsurface drainage tiles has altered connectivity between the land surface, groundwater, and streams, and where application of fertilizer has created significant stores of subsurface nitrogen. In this study, we identify threshold patterns of tile‐runoff generation for a drained agricultural field in Illinois and evaluate how antecedent conditions—including shallow soil moisture, groundwater table depth, and the presence or absence of crops—control tile response. We relate tile‐runoff thresholds to patterns of event nitrate load observed across multiple storm events and evaluate how antecedent conditions control within‐event nitrate concentration‐discharge relationships. Our results demonstrate that an event tile‐runoff threshold emerges relative to the sum of gross precipitation and indices of antecedent shallow soil moisture and antecedent below‐tile groundwater moisture deficit, indicating that both shallow soil and below‐tile storages must be filled to generate significant runoff. In turn, event nitrate load shows a linear dependence on runoff for most time periods, suggesting that subsurface nitrate export and storage can be estimated using runoff threshold relationships and long‐term average nitrate concentrations. Finally, within‐event nitrate concentration‐discharge relationships are controlled by event size and the antecedent tile flow state because these factors dictate the sequence of flow path activation and tile connectivity over a storm event.

    more » « less
  5. Abstract

    This study focuses on characterizing the contributions of key terrestrial pathways that deliver dissolved organic carbon (DOC) to streams during hydrological events and on elucidating factors governing variation in water and DOC fluxes from these pathways. We made high‐frequency measurements of discharge, specific conductance (SC), and fluorescent dissolved organic matter (FDOM) during 221 events recorded over 2 years within four Vermont (USA) watersheds that range in area from 0.4 to 139 km2. Using the SC measurements, together with statistical information on discharge, we separated the event hydrographs into contributions from three terrestrial pathways, which we refer to as riparian quickflow, subsurface quickflow, and slow‐flow groundwater. The pathway discharges were used as input to a mixing model that closely approximated sub‐hourly streamwater DOC concentrations as measured with the FDOM sensors. Subsurface quickflow, comprised of pre‐event water, was the leading contributor to streamwater DOC fluxes, while riparian quickflow, comprised of event water, was the second‐leading contributor to streamwater DOC fluxes, despite comprising the smallest proportion of streamflow yield among the three end‐member pathways. Fixed‐effects regression analysis revealed that the relationship between DOC fluxes from the end‐member pathways and event magnitude was consistent across the four watersheds. This analysis also showed that DOC fluxes from the quickflow pathways increased significantly with temperature and varied inversely, but weakly, with catchment antecedent wetness. We believe that our approach, which leverages in‐stream sensors that enable high‐frequency measurements over extended periods, may be applicable for evaluating controls on DOC export from other watersheds within and beyond our study region.

    more » « less