skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Importance of Capturing Topographic Features for Modeling Groundwater Flow and Transport in Mountainous Watersheds
Abstract Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds.  more » « less
Award ID(s):
1830172
PAR ID:
10453285
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
54
Issue:
12
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. How does the physical and chemical structure of the Critical Zone (CZ), defined as the zone from treetops to the bottom of groundwater, govern its hydro-biogeochemical functioning? Multiple lines of evidence from past and newly emerging research have prompted the shallow and deep partitioning concentration-discharge (C-Q) hypothesis. The hypothesis states that in-stream C-Q relationships are shaped by distinct source waters from flow paths at different depths. Base flows are often dominated by deep groundwater and mostly reflect groundwater chemistry, whereas high flows are often dominated by shallow soil water and thus mostly reflect soil water chemistry. The contrasts between shallow soil water versus deeper groundwater chemistry shape stream solute export patterns. In this context, the vertical connectivity that regulates the shallow and deep flow partitioning is essential in determining chemical contrasts, biogeochemical reaction rates in soils and parent rocks, and ultimately solute export patterns. This talk will highlight insights gleaned from multiple lines of recent studies that include collation of water chemistry data from soils, rocks, and streams in intensively monitored watersheds, meta-analysis of stream chemistry data at the continental scale, and integrated reactive transport modeling at the hillslope and watershed scales. The hypothesis underscores the importance of subsurface vertical structure and connectivity relative to the extensively studied horizontal connectivity. It also alludes to the potential of using streams as mirrors for subsurface water chemistry, and the potential of using C-Q relationships to infer flow paths and biogeochemical reaction rates and the response of earth’s subsurface to climate and human perturbations. Broadly, this simple conceptual framework links CZ subsurface structure to its functioning under diverse climate, geology, and land cover conditions. 
    more » « less
  2. Abstract Mountain System Recharge processes are significant natural recharge pathways in many arid and semi‐arid mountainous regions. However, Mountain System Recharge processes are often poorly understood and characterized in hydrologic models. Mountains are the primary water supply source to valley aquifers via lateral groundwater flow from the mountain block (Mountain Block Recharge) and focused recharge from mountain streams contributing to focused Mountain Front Recharge at the piedmont zone. Here, we present a multi‐tool isogeochemical approach to characterize mountain flow paths and Mountain System Recharge in the northern Tulare Basin, California. We used groundwater chemistry data to delineate hydrochemical facies and explain the chemical evolution of groundwater from the Sierra Nevada to the Central Valley aquifer. Stable isotopes and radiogenic groundwater tracers validated Mountain System Recharge processes by differentiating focused from diffuse recharge, and estimating apparent groundwater age, respectively. Novel application of End‐Member Mixing Analysis using conservative chemical components revealed three Mountain System Recharge end‐members: (a) evaporated Ca‐HCO3water type associated with focused Mountain Front Recharge, (b) non‐evaporated Ca‐HCO3and Na‐HCO3water types with short residence times associated with shallow Mountain Block Recharge, and (c) Na‐HCO3groundwater type with long residence time associated with deep Mountain Block Recharge. We quantified the contribution of each Mountain System Recharge process to the valley aquifer by calculating mixing ratios. Our results show that deep Mountain Block Recharge is a significant recharge component, representing 31%–53% of the valley groundwater. Greater hydraulic connectivity between the Sierra Nevada and Central Valley has significant implications for parameterizing groundwater flow models. Our framework is useful for understanding Mountain System Recharge processes in other snow‐dominated mountain watersheds. 
    more » « less
  3. Abstract High latitude regions across the globe are undergoing severe modifications due to changing climate. A high latitude region of concern is the Gulf of Alaska (GoA), where these changes in hydroclimate undoubtedly affect the hydrogeochemistry of freshwater discharging to the nearshore ecosystems of the region. To fill the knowledge gap of our understanding of freshwater stream geochemistry with the GoA, we compile stream water chemistry data from 162 stream sites across the region. With an inverse model, we estimate fractional contributions to solute fluxes from weathering of silicate, carbonate, and sulfide minerals, and precipitation. We assess weathering rates across the region and compare them against global river yields. The median fractional contribution of carbonate weathering to total weathering products is 78% across all stream sites; however, there are several streams where silicate weathering is a dominant source of solutes. Weathering by sulfuric acid is elevated in glacierized watersheds. Finally, cation weathering rates are lower in GoA streams compared to the world's largest rivers; however, weathering rates are similar when compared to a global dataset of glacier fed streams. We suggest that hydrologic changes driven by glacier ice loss and increased precipitation will alter river water quality and chemical weathering regimes such that silicate weathering may become a more important source of solutes and sulfide oxidation may decrease. This contribution provides a platform to build from for future investigations into changes to stream water chemistry in the region and other high latitude watersheds. 
    more » « less
  4. High elevation mountain watersheds are undergoing rapid warming and declining snow fractions worldwide, causing earlier and quicker snowmelt. Understanding how this hydrologic shift affects subsurface flow paths, biogeochemical reactions, and solute export has been challenging due to the entanglement of hydrological and biogeochemical processes. Coal Creek, a high-elevation catchment (2,700 3,700 m, 53 km2) in Colorado, is experiencing a higher rate of warming than surrounding low-lying areas. This warming corresponds with dynamic and increased responses from biogenic solutes and dissolved organic carbon (DOC), whereas the behavior of geogenic solutes and dissolved inorganic carbon (DIC) has remained relatively unchanged. DOC has experienced the largest concentration increase (>3x), with annual average flow weighted concentrations positively correlated to average annual temperature. This suggests temperature is the main driver of increasing DOC levels. Although DOC and DIC response to warming is influenced by many drivers, the relative contribution of each remains unknown. DOC and DIC were analyzed to incorporate both carbon component products of soil respiration (DOC and CO2) and to represent high solute concentrations transported by shallow (DOC) versus deep (DIC) subsurface flow. The contrasting behavior of these carbon solutes indicates climate change and warming are driving changes in organic matter decomposition and soil respiration. Modeling results from the process-based model HBV-BioRT show increased temperatures cause earlier snowmelt and streamflow generation and lower peak discharge. As stream flow generation occurs earlier, so do DOC flushing and DIC dilution events. Additionally, post-snowmelt periods show greater DOC production and concentrations under warming scenarios. Results indicated increased production of DOC in post-snowmelt periods. DOC is then flushed out by earlier snowmelt partitioned through the shallow soil zone. Most process-based studies lack a watershed-scale understanding of carbon transformation and flow path alterations. This work demonstrates complex hydrologic and biogeochemical coupling at the watershed scale to illustrate how water flow paths and chemistry are responding to a changing climate in highelevation mountain watersheds. 
    more » « less
  5. Soil biota generate CO2 that can vertically export to the atmosphere, and dissolved organic and inorganic carbon (DOC and DIC) that can laterally export to streams and accelerate weathering. These processes are regulated by external hydroclimate forcing and internal structures (permeability distribution), the relative influences of which are rarely studied. Understanding these interactions is essential a hydrological extremes intensify in the future. Here we explore the question: How and to what extent do hydrological and permeability distribution conditions regulate soil carbon transformations and chemical weathering? We address the questions using a hillslope reactive transport model constrained by data from the Fitch Forest (Kansas, United States). Numerical experiments were used to mimic hydrological extremes and variable shallow-versus-deep permeability contrasts. Results demonstrate that under dry conditions (0.08 mm/day), long water transit times led to more mineralization of organic carbon (OC) into inorganic carbon (IC) form (>98\%). Of the IC produced, ~ 75\% was emitted upward as CO2 gas and ~ 25\% was exported laterally as DIC into the stream. Wet conditions (8.0 mm/day) resulted in less mineralization (~88\%), more DOC production (~12\%), and more lateral fluxes of IC (~50\% of produced IC). Carbonate precipitated under dry conditions and dissolved under wet conditions as the fast flow rapidly droves the reaction to disequilibrium. The results depict a conceptual hillslope model that prompts four hypotheses for our community to test. H1: Droughts enhance carbon mineralization and vertical upward carbon fluxes, whereas large hydrological events such as storms and flooding enhance subsurface vertical connectivity, reduce transit times, and promote lateral export. H2: The role of weathering as a net carbon sink or source to the atmosphere depends on the interaction between hydrologic flows and lithology: transition from droughts to storms can shift carbonate from a carbon sink (mineral precipitation) to carbon source (dissolution). H3: Permeability contrasts regulate the lateral flow partitioning via shallow flow paths versus deeper groundwater though this alter reaction rates negligibly. H4: Stream chemistry reflect flow paths and can potentially quantify water transit times: solutes enriched in shallow soils have a younger water signature; solutes abundant at depth carry older water signature. 
    more » « less