Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and downstream of the fractured hillslope (non- fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone were high, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. We interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.
more »
« less
A Coupled SWAT-AEM Modelling Framework for a Comprehensive Hydrologic Assessment
This study attempts to integrate a Surface Water (SW) model Soil and Water Assessment Tool (SWAT) with an existing steady-state, single layer, unconfined heterogeneous aquifer Analytic Element Method (AEM) based Ground Water (GW) model, named Bluebird AEM engine, for a comprehensive assessment of SW and GW resources and its management. The main reason for integrating SWAT with the GW model is that the SWAT model does not simulate the distribution and dynamics of GW levels and recharge rates. To overcome this issue, often the SWAT model is coupled with the numerical GW model (either using MODFLOW or FEFLOW), wherein the spatial and temporal patterns of the interactions are better captured and assessed. However, the major drawback in integrating the two models (SWAT with—MODFLOW/FEM) is its conversion from Hydrological Response Unit’s (HRU)/sub-basins to grid/elements. To couple them, a spatial translation system is necessary to move the inputs and outputs back and forth between the two models due to the difference in discretization. Hence, for effective coupling of SW and GW models, it may be desirable to have both models with a similar spatial discretization and reduce the need for rigorous numerical techniques for solving the PDEs. The objective of this paper is to test the proof of concept of integrating a distributed hydrologic model with an AEM model at the same spatial units, primarily focused on surface water and groundwater interaction with a shallow unconfined aquifer. Analytic Element Method (AEM) based GW models seem to be ideal for coupling with SWAT due to their innate character to consider the HRU, sub-basin, River, and lake boundaries as individual analytic elements directly without the need for any further discretization or modeling units. This study explores the spatio-temporal patterns of groundwater (GW) discharge rates to a river system in a moist-sub humid region with SWAT-AEM applied to the San Jacinto River basin (SJRB) in Texas. The SW-GW interactions are explored throughout the watershed from 2000–2017 using the integrated SWAT-AEM model, which is tested against stream flow and GW levels. The integrated SWAT-AEM model results show good improvement in predicting the stream flow (R2 = 0.65–0.80) and GW levels as compared to the standalone SWAT model. Further, the integrated model predicted the low flows better compared to the standalone SWAT model, thus accounting for the SW-GW interactions. Almost 80% of the stream network experiences an increase in groundwater discharge rate between 2000 and 2017 with an annual average GW discharge rate of 1853 Mm3/year. The result from the study seems promising for potential applications of SWAT-AEM coupling in regions with considerable SW-GW interactions.
more »
« less
- Award ID(s):
- 1739977
- PAR ID:
- 10463744
- Date Published:
- Journal Name:
- Water
- Volume:
- 14
- Issue:
- 17
- ISSN:
- 2073-4441
- Page Range / eLocation ID:
- 2753
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Floodplains are essential ecosystems that provide a variety of economic, hydrologic, and ecologic services. Within floodplains, surface water‐groundwater exchange plays an important role in facilitating biogeochemical processes and can have a strong influence on stream hydrology through infiltration or discharge of water. These functions can be difficult to assess due to the heterogeneity of floodplains and monitoring constraints, so numerical models are useful tools to estimate fluxes, especially at large spatial extents. In this study, we use the SWAT+ (Soil and Water Assessment Tool) ecohydrological model to quantify magnitudes and spatiotemporal patterns of floodplain surface water‐groundwater exchange in a mountainous watershed using an updated version of thegwflowmodule that directly calculates floodplain‐aquifer exchange rates during periods of floodplain inundation. Thegwflowmodule is a spatially distributed groundwater modelling subroutine within the SWAT+ code that uses a gridded network and physically based equations to predict groundwater storage, groundwater head, and groundwater fluxes. We used SWAT+ to model the 7516 km2Colorado River headwaters watershed and streamflow data from USGS gages for calibration and testing. Models that included floodplain‐groundwater interactions outperformed those without such interactions and provided valuable information about floodplain exchange rates and volumes. Our analyses on the location of floodplain fluxes in the watershed also show that wider areas of floodplains, “beads” (e.g., like beads on a necklace), exchanged a higher net and per area volume of water, as well as higher rates of exchange, compared to narrower areas, “strings.” Study results show that floodplain channel‐groundwater exchange is a valuable process to include in hydrologic models, and model outputs could inform land conservation practises by indicating priority locations, such as beads, where substantial hydrologic exchange occurs.more » « less
-
Abstract Groundwater discharge zones connect aquifers to surface water, generating baseflow and serving as ecosystem control points across aquatic ecosystems. The influence of groundwater discharge on surface flow connectivity, fate and transport of contaminants and nutrients, and thermal habitat depends strongly on hydrologic characteristics such as the spatial distribution, age, and depth of source groundwater flow paths. Groundwater models have the potential to predict spatial discharge characteristics within river networks, but models are often not evaluated against these critical characteristics and model equifinality with respect to discharge processes is a known challenge. We quantify discharge characteristics across a suite of groundwater models with commonly used frameworks and calibration data. We developed a base model (MODFLOW‐NWT) for a 1,570‐km2watershed in the northeastern United States and varied the calibration data, control of river‐aquifer exchange directionality, and resolution. Most models (n = 11 of 12) fit similarly to calibration metrics, but patterns in discharge location, flow path depth, and subsurface travel time varied substantially. We found (1) a 15% difference in the percent of discharge going to first‐order streams, (2) threefold variations in flow path depth, and (3) sevenfold variations in the subsurface travel times among the models. We recalibrated three models using a synthetic discharge location data set. Calibration with discharge location data reduced differences in simulated discharge characteristics, suggesting an approach to improved equifinality based on widespread field‐based mapping of discharge zones. Our work quantifying variation across common modeling approaches is an important step toward characterizing and improving predictions of groundwater discharge characteristics.more » « less
-
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
-
Groundwater use for irrigation has a major influence on agricultural productivity and local water resources. This study evaluated the groundwater irrigation schemes, SWAT auto-irrigation scheduling based on plant water stress (Auto-Irr), and prescribed irrigation based on well pumping rates in MODFLOW (Well-Irr), in the U.S. Northern High Plains (NHP) aquifer using coupled SWAT-MODFLOW model simulations for the period 1982–2008. Auto-Irr generally performed better than Well-Irr in simulating groundwater irrigation volume (reducing the mean bias from 86 to −30%) and groundwater level (reducing the normalized root-mean-square-error from 13.55 to 12.47%) across the NHP, as well as streamflow interannual variations at two stations (increasing NSE from 0.51, 0.51 to 0.55, 0.53). We also examined the effects of groundwater irrigation on the water cycle. Based on simulation results from Auto-Irr, historical irrigation led to significant recharge along the Elkhorn and Platte rivers. On average over the entire NHP, irrigation increased surface runoff, evapotranspiration, soil moisture and groundwater recharge by 21.3%, 4.0%, 2.5% and 1.5%, respectively. Irrigation improved crop water productivity by nearly 27.2% for corn and 23.8% for soybean. Therefore, designing sustainable irrigation practices to enhance crop productivity must consider both regional landscape characteristics and downstream hydrological consequences.more » « less
An official website of the United States government

