skip to main content


Title: 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
NSF-PAR ID:
10463744
Author(s) / Creator(s):
; ;
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
  1. Abstract

    Several studies have focused on the importance of river bathymetry (channel geometry) in hydrodynamic routing along individual reaches. However, its effect on other watershed processes such as infiltration and surface water (SW)‐groundwater (GW) interactions has not been explored across large river networks. Surface and sbsurface processes are interdependent, therefore, errors due to inaccurate representation of one watershed process can cascade across other hydraulic or hydrologic processes. This study hypothesizes that accurate bathymetric representation is not only essential for simulating channel hydrodynamics but also affects subsurface processes by impacting SW‐GW interactions. Moreover, quantifying the effect of bathymetry on surface and subsurface hydrological processes across a river network can facilitate an improved understanding of how bathymetric characteristics affect these processes across large spatial domains. The study tests this hypothesis by developing physically based distributed models capable of bidirectional coupling (SW‐GW) with four configurations with progressively reduced levels of bathymetric representation. A comparison of hydrologic and hydrodynamic outputs shows that changes in channel geometry across the four configurations has a considerable effect on infiltration, lateral seepage, and location of water table across the entire river network. For example, when using bathymetry with inaccurate channel conveyance capacity but accurate channel depth, peak lateral seepage rate exhibited 58% error. The results from this study provide insights into the level of bathymetric detail required for accurately simulating flooding‐related physical processes while also highlighting potential issues with ignoring bathymetry across lower order streams such as spurious backwater flow, inaccurate water table elevations, and incorrect inundation extents.

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

    In mountainous river basins of the Pacific Northwest, climate models predict that winter warming will result in increased precipitation falling as rain and decreased snowpack. A detailed understanding of the spatial and temporal dynamics of water sources across river networks will help illuminate climate change impacts on river flow regimes. Because the stable isotopic composition of precipitation varies geographically, variation in surface water isotope ratios indicates the volume‐weighted integration of upstream source water. We measured the stable isotope ratios of surface water samples collected in the Snoqualmie River basin in western Washington over June and September 2017 and the 2018 water year. We used ordinary least squares regression and geostatistical Spatial Stream Network models to relate surface water isotope ratios to mean watershed elevation (MWE) across seasons. Geologic and discharge data was integrated with water isotopes to create a conceptual model of streamflow generation for the Snoqualmie River. We found that surface water stable isotope ratios were lowest in the spring and highest in the dry, Mediterranean summer, but related strongly to MWE throughout the year. Low isotope ratios in spring reflect the input of snowmelt into high elevation tributaries. High summer isotope ratios suggest that groundwater is sourced from low elevation areas and recharged by winter precipitation. Overall, our results suggest that baseflow in the Snoqualmie River may be relatively resilient to predicted warming and subsequent changes to snowpack in the Pacific Northwest.

     
    more » « less
  5. Abstract Non-perennial rivers and streams make up over half the global river network and are becoming more widespread. Transitions from perennial to non-perennial flow are a threshold-type change that can lead to alternative stable states in aquatic ecosystems, but it is unknown whether streamflow itself is stable in either wet (flowing) or dry (no-flow) conditions. Here, we investigated drivers and feedbacks associated with regime shifts between wet and dry conditions in an intermittent reach of the Arkansas River (USA) over the past 23 years. Multiple lines of evidence suggested that these regimes represent alternative stable states, including (a) significant jumps in discharge time series that were not accompanied by jumps in flow drivers such as precipitation and groundwater pumping; (b) a multi-modal state distribution with 92% of months experiencing no-flow conditions for <10% or >90% of days, despite unimodal distributions of precipitation and pumping; and (c) a hysteretic relationship between climate and flow state. Groundwater levels appear to be the primary control over the hydrological regime, as groundwater levels in the alluvial aquifer were higher than the stream stage during wet regimes and lower than the streambed during dry regimes. Groundwater level variation, in turn, was driven by processes occurring at both the regional scale (surface water inflows from upstream, groundwater pumping) and the reach scale (stream–aquifer exchange, diffuse recharge through the soil column). Historical regime shifts were associated with diverse pressures including network disconnection caused by upstream water use, increased flow stability potentially associated with reservoir operations, and anomalous wet and dry climate conditions. In sum, stabilizing feedbacks among upstream inflows, stream–aquifer interactions, climate, vegetation, and pumping appear to create alternative wet and dry stable states at this site. These stabilizing feedbacks suggest that widespread observed shifts from perennial to non-perennial flow will be difficult to reverse. 
    more » « less