Using SWAT-MODFLOW to simulate groundwater flow and groundwater-surface water interactions in an intensively irrigated stream-aquifer system

Agricultural irrigated watersheds in semi-arid regions face challenges such as waterlogging, high soil salinity, reduced crop yield, and leaching of chemical species due to extreme shallow water tables resulting from long-term intensive irrigation. Hydrologic models can be used to evaluate the impact of land management practices on water yields and groundwater-surface water interactions in such regions. In this study, the newly developed SWAT-MODFLOW, a coupled surface/subsurface hydrologic model, is applied to a 950 km2 watershed in the Lower Arkansas River Valley (southeastern Colorado). The model accounts for the influence of canal diversions, irrigation applications, groundwater pumping, and earth canal seepage losses. The model provides a detailed description of surface and subsurface flow processes, thereby enabling detailed description of watershed processes such as runoff, infiltration, in-streamflow, three-dimensional groundwater flow in a heterogeneous aquifer system with sources and sinks (e.g. pumping, seepage to subsurface drains), and spatially-variable surface and groundwater exchange. The model was calibrated and tested against stream discharge from 5 stream gauges in the Arkansas River and its tributaries, groundwater levels from 70 observation wells, and evapotranspiration (ET) data estimated from satellite (ReSET) data during the 1999 to 2007 period. Since the water-use patterns within the study area are typical of many other irrigated river valleys in the United States and elsewhere, this modeling approach is transferable to other regions.