skip to main content

This content will become publicly available on October 31, 2024

Title: Spatially distributed tracer‐aided modelling to explore DOC dynamics, hot spots and hot moments in a tropical mountain catchment

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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Hydrological Processes
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Quantitative estimations of ecohydrological water partitioning into evaporation and transpiration remains mostly based on plot‐scale investigations that use well‐instrumented, small‐scale experimental catchments in temperate regions. Here, we attempted to upscale and adapt the conceptual tracer‐aided ecohydrology model STARRtropics to simulate water partitioning, tracer, and storage dynamics over daily time steps and a 1‐km grid larger‐scale (2565 km2) in a sparsely instrumented tropical catchment in Costa Rica. The model was driven by bias‐corrected regional climate model outputs and was simultaneously calibrated against daily discharge observations from 2 to 30 years at four discharge gauging stations and a 1‐year, monthly streamwater isotope record of 46 streams. The overall model performance for the best discharge simulations ranged in KGE values from 0.4 to 0.6 and correlation coefficients for streamflow isotopes from 0.3 to 0.45. More importantly, independent model‐derived transpiration estimates, point‐scale residence time estimates, and measured groundwater isotopes showed reasonable model performance and simulated spatial and temporal patterns pointing towards an overall model realism at the catchment scale over reduced performance in the headwaters. The simulated catchment system was dominated by low‐seasonality and high precipitation inputs and a marked topographical gradient. Climatic drivers overrode smaller, landcover‐dependent transpiration fluxes giving a seemingly homogeneous rainfall‐runoff dominance likely related to model input bias of rainfall isotopes, oversimplistic Potential Evapotranspiration (PET) estimates and averaged Leaf Area Index (LAI). Topographic influences resulted in more dynamic water and tracer fluxes in the headwaters that averaged further downstream at aggregated catchment scales. Modelled headwaters showed greater storage capacity by nearly an order of magnitude compared to the lowlands, which also favoured slightly longer residence times (>250 days) compared to superficially well‐connected groundwater contributing to shorter streamflow residence times (<150 days) in the lowlands. Our findings confirm that tracer‐aided ecohydrological modelling, even in the data‐scarce Tropics, can help gain a first, but crucial approximation of spatio‐temporal dynamics of how water is partitioned, stored and transported beyond the experimental catchment scale of only a few km2.

    more » « less
  2. Tracer-aided hydrological models (TAHMs) are one of the most powerful tools to identify new (event) and old (pre-event) water fractions contributing to stormflow because they account both for streamflow and tracer mixing dynamics in model calibration. Nevertheless, their representativeness of hydrograph dynamics is often limited due to the unavailability of high-resolution conservative tracer data (e.g., water stable isotopes or chloride). Hence, there is a need to identify alternative tracers yielding similar flow partitioning results than “ideal” ones while requiring fewer financial resources for high-frequency monitoring (e.g., sub-hourly). Here, wecompare flow partitioning results of a TAHM calibrated using high-frequency electrical conductivity (EC) and water stable isotope (18O) data collected during 37 rainfall-runoff events monitored during variable hydrometeorological conditions in the Zhurucay Ecohydrological Observatory, a tropical alpine catchment located in southern Ecuador. When the model was calibrated using the sampling resolution of stables isotopes (6-hours to 1-hour), no statistically significant differences of pre-event water fractions (PEWFs) using both tracers for model calibration were found. PEWF differences between both tracers for 89% of the events were < 20% regardless of the events’ antecedent moisture and rainfall conditions. Model transfer functions were also similar suggesting that catchment internal processes inferred using both tracers are comparable. Events presenting larger differences (n = 4; up to 27% PEWF difference) had no samples collected during peak flow. Calibration of the model using EC data collected at sub-hourly intervals (every 5-minutes) showed a significant increase in model performance as compared to the frequency of collection of isotopic data. Similarity in flow partitioning results can be attributed to a quasi-conservative nature of EC due to the presence of organic-rich riparian soils (peat-type) overlying compact bedrock across the catchment. Findings also highlight the importance of capturing rapidly occurring catchment mixing processes though high-temporal frequency monitoring of tracer data. Our study encourages the value of assessing the use of alternative tracers, such as EC, to identify fast occurring rainfall runoff processes, while lowering the costs needed to implement and sustain tracer data collection for long time periods. 
    more » « less
  3. Arctic landscapes are warming and becoming wetter due to changes in precipitation and the timing of snowmelt which consequently alters seasonal runoff and river discharge patterns. These changes in hydrology lead to increased mobilization and transport of terrestrial dissolved organic matter (DOM) to Arctic coastal seas where significant impacts on biogeochemical cycling can occur. Here, we present measurements of dissolved organic carbon (DOC) and chromophoric DOM (CDOM) in the Yukon River-to-Bering Sea system and two river plumes on the Alaska North Slope which flow into the Beaufort Sea. Our sampling characterized optical and biogeochemical properties of DOM during high and low river discharge periods for the Yukon River-Bering Sea system. The average DOC concentration at the multiple Yukon River mouths ranged from a high of 10.36 mg C L -1 during the ascending limb of the 2019 freshet (late May), 6.4 mg C L -1 during the descending limb of the 2019 freshet (late June), and a low of 3.86 mg C L -1 during low river discharge in August 2018. CDOM absorption coefficient at 412 nm ( a CDOM (412)) averaged 8.23 m -1 , 5.07 m -1 , and 1.9 m -1 , respectively. Several approaches to model DOC concentration based on its relationship with CDOM properties demonstrated cross-system seasonal and spatial robustness for these Arctic coastal systems despite spanning an order of magnitude decrease in DOC concentration from the lower Yukon River to the Northern Bering Sea as well as the North Slope systems. “Snapshot” fluxes of DOC and CDOM across the Yukon River Delta to Norton Sound were calculated from our measurements and modeled water fluxes forced with upstream USGS river gauge data. Our findings suggest that during high river flow, DOM reaches the delta largely unaltered by inputs or physical and biogeochemical processing and that the transformations of Yukon River DOM largely occur in the plume. However, during low summer discharge, multiple processes including local precipitation events, microbial decomposition, photochemistry, and likely others can alter the DOM properties within the lower Yukon River and Delta prior to flowing into Norton Sound. 
    more » « less
  4. Abstract

    Non‐perennial headwaters experience extremes in flow conditions that likely influence carbon fate. As surface waters contract through dry periods, reconnect during storms, and re‐expand or dry again, there is a great deal of variability in carbon emissions and export. We measured discharge, dissolved oxygen, carbon dioxide (CO2), and dissolved organic carbon (DOC) continuously in a persistent pool at the base of a non‐perennial, forested headwater stream in the southeastern United States to characterize how flow changes affect carbon emissions and export as the stream expands and shrinks. We also compared carbon concentrations and export during different stream flow categories before and after fall wet‐up. CO2concentrations were high when discharge was lowest (median = 10.2 mg L−1) and low during high flows (3.2 mg L−1) and storms (1.1 mg L−1). High CO2concentrations led to high emissions on a per area basis during low flow times, but whole‐channel stream CO2emissions were limited by the small surface area of the stream during periods of surface water disconnection. DOC concentration varied by season (range = 0.1–16.2 mg L−1) with large pulses during smaller summer storms. We found that CO2and DOC concentrations differed among binned stages of stream flow. As non‐perennial streams become more prevalent across the southeastern United States due to shifts in climate, the relationships between flow and carbon movement into and out of stream networks will become increasingly critical to understanding stream carbon biogeochemistry.

    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