Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparison of two modeling approaches for groundwater–surface water interactions

An assessment of interactions between groundwater and surface water was carried out by applying two different modeling approaches to a small‐scale study area in the municipality of Havelock, Quebec. The first approach involved a commonly used sequential procedure that consists in determining the daily recharge rate using a quasi 2D infiltration model (HELP), applied in the next step as an imposed flux to a 3D finite‐element groundwater flow model. The flow model was calibrated under steady‐state and transient conditions against measured water levels. The second approach was based on a recently developed physically based, 3D fully coupled groundwater–surface water flow model (CATHY) applied to the entire flow domain in an integrated manner. Implementation, calibration, and results of the simulations for both approaches are presented and discussed. For equal annual precipitation (1038 mm/y) and evapotranspiration (556 mm/y), the second approach computed a recharge rate of 233 mm/y (8.9% higher than the first approach) and a net upward flow from the fractured aquifer (the first approach predicted a net downward flow to the rock). The simulated annual discharge was similar for the two approaches (9.6% difference). Both approaches were found to be useful in understanding the interactions between groundwater and surface water, although limitations are apparent in the sequential procedure's inability to account for surface–subsurface feedbacks, for instance near stream reaches where groundwater discharge is prevalent. The decoupled, two‐model approach provides disaggregated surface, vadose, and aquifer flows, and a simple aperçu at the different components of total discharge. The fully coupled model accounts for continuous water exchanges between the land surface, subsurface, and stream channel in a more complex manner, and produces a better match against observed data. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved simulation of river water and groundwater exchange in an alluvial plain using the SWAT model

Hydrological interaction between surface and subsurface water systems has a significant impact on water quality, ecosystems and biogeochemistry cycling of both systems. Distributed models have been developed to simulate this function, but they require detailed spatial inputs and extensive computation time. The soil and water assessment tool (SWAT) model is a semi‐distributed model that has been successfully applied around the world. However, it has not been able to simulate the two‐way exchanges between surface water and groundwater. In this study, the SWAT‐landscape unit (LU) model – based on a catena method that routes flow across three LUs (the divide, the hillslope and the valley) – was modified and applied in the floodplain of the Garonne River. The modified model was called SWAT‐LUD. Darcy's equation was applied to simulate groundwater flow. The algorithm for surface water‐level simulation during flooding periods was modified, and the influence of flooding on groundwater levels was added to the model. Chloride was chosen as a conservative tracer to test simulated water exchanges. The simulated water exchange quantity from SWAT‐LUD was compared with the output of a two‐dimensional distributed model, surface–subsurface water exchange model. The results showed that simulated groundwater levels in the LU adjoining the river matched the observed data very well. Additionally, SWAT‐LUD model was able to reflect the actual water exchange between the river and the aquifer. It showed that river water discharge has a significant influence on the surface–groundwater exchanges. The main water flow direction in the river/groundwater interface was from groundwater to river; water that flowed in this direction accounted for 65% of the total exchanged water volume. The water mixing occurs mainly during high hydraulic periods. Flooded water was important for the surface–subsurface water exchange process; it accounted for 69% of total water that flowed from the river to the aquifer. The new module also provides the option of simulating pollution transfer occurring at the river/groundwater interface at the catchment scale. Copyright © 2015 John Wiley & Sons, Ltd.     

Integrated modeling as a decision‐aiding tool for groundwater management in a Mediterranean agricultural watershed

A decision‐aiding methodology for agricultural groundwater management is presented; it is based on the combination of a watershed model, a groundwater flow model, and an optimization model. This methodology was applied to an agricultural watershed in northeastern Greece. The watershed model used was the Soil and Water Assessment Tool (SWAT), which provided recharge rates for the aquifers. These recharge rates were imported in the well‐known MODFLOW groundwater flow model. Both models were calibrated and verified using field data. Then, the nonlinear optimization problem was solved by a piecewise linearization process, in which the Simplex algorithm was applied sequentially. Apart from several pumping and climate change sensitivity scenarios, a land use change scenario and a climate change scenario, combining the three models, were tested, showing the ability of this methodology to be used in the decision‐making process. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of stream flow components in a mountainous catchment in northern Thailand with SWAT,using the ANSELM calibration approach

Highland agriculture is intensifying rapidly in South‐East Asia, leading to alarmingly high applications of agrochemicals. Understanding the fate of these contaminants requires carefully planned monitoring programmes and, in most cases, accurate simulation of hydrological pathways into and through water bodies. We simulate run‐off in a steep mountainous catchment in tropical South‐East Asia. To overcome calibration difficulties related to the mountainous topography, we introduce a new calibration method, named A Nash–Sutcliffe Efficiency Likelihood Match (ANSELM), that allows the assignment of optimal parameters to different hydrological response units in simulations of stream discharge with the Soil and Water Assessment Tool (SWAT) hydrological model. ANSELM performed better than the Parasol calibration tool built into SWAT in terms of model efficiency and computation time. In our simulation, the most sensitive model parameters were those related to base flow generation, surface run‐off generation, flow routing and soil moisture change. The coupling of SWAT with ANSELM yielded reasonable simulations of both wet‐season and dry‐season storm hydrographs. Nash–Sutcliffe model efficiencies for daily stream flow during two validation years were 0.77 and 0.87. These values are in the upper range or even higher than those reported for other SWAT model applications in temperate or tropical regions. The different flow components were realistically simulated by SWAT, and showed a similar behaviour in all the study years, despite inter‐annual climatic differences. The realistic partitioning of total stream flow into its contributing components will be an important factor for using this hydrological model to simulate solute transport in the future. Copyright © 2014 John Wiley & Sons, Ltd.     

A multi‐storage groundwater concept for the SWAT model to emphasize nonlinear groundwater dynamics in lowland catchments

Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.     

Predicting saturation‐excess runoff distribution with a lumped hillslope model: SWAT‐HS

Saturation‐excess runoff is the major runoff mechanism in humid well‐vegetated areas where infiltration rates often exceed rainfall intensity. Although the Soil and Water Assessment Tool (SWAT) is one of the most widely used models, it predicts runoff based mainly on soil and land use characteristics, and is implicitly an infiltration‐excess runoff type of model. Previous attempts to incorporate the saturation‐excess runoff mechanism in SWAT fell short due to the inability to distribute water from one hydrological response unit to another. This paper introduces a modified version of SWAT, referred to as SWAT‐Hillslope (SWAT‐HS). This modification improves the simulation of saturation‐excess runoff by redefining hydrological response units based on wetness classes and by introducing a surface aquifer with the ability to route interflow from “drier” to “wetter” wetness classes. Mathematically, the surface aquifer is a nonlinear reservoir that generates rapid subsurface stormflow as the water table in the surface aquifer rises. The SWAT‐HS model was tested in the Town Brook watershed in the upper reaches of the West Branch Delaware River in the Catskill region of New York, USA. SWAT‐HS predicted discharge well with a Nash‐Sutcliffe Efficiency of 0.68 and 0.87 for daily and monthly time steps. Compared to the original SWAT model, SWAT‐HS predicted less surface runoff and groundwater flow and more lateral flow. The saturated areas predicted by SWAT‐HS were concentrated in locations with a high topographic index and were in agreement with field observations. With the incorporation of topographic characteristics and the addition of the surface aquifer, SWAT‐HS improved streamflow simulation and gave a good representation of saturated areas on the dates that measurements were available. SWAT‐HS is expected to improve water quality model predictions where the location of the surface runoff matters.     

Quantifying groundwater–surface water interactions in a proglacial valley,Cordillera Blanca,Peru

A myriad of downstream communities and industries rely on streams fed by both groundwater discharge and glacier meltwater draining the Cordillera Blanca, Northern Peruvian Andes, which contains the highest density of glaciers in the tropics. During the dry season, approximately half the discharge in the region's proglacial streams comes from groundwater. However, because of the remote and difficult access to the region, there are few field methods that are effective at the reach scale to identify the spatial distribution of groundwater discharge. An energy balance model, Rhodamine WT dye tracing, and high‐definition kite‐borne imagery were used to determine gross and net groundwater inputs to a 4‐km reach of the Quilcay River in Huascaran National Park, Peru. The HFLUX computer programme ( http://hydrology.syr.edu/hflux.html ) was used to simulate the Quilcay River's energy balance using stream temperature observations, meteorological measurements, and kite‐borne areal photography. Inference from the model indicates 29% of stream discharge at the reach outlet was contributed by groundwater discharge over the study section. Rhodamine WT dye tracing results, coupled with the energy balance, show that approximately 49% of stream water is exchanged (no net gain) with the subsurface as gross gains and losses. The results suggest that gross gains from groundwater are largest in a moraine subreach but because of large gross losses, net gains are larger in the meadow subreaches. These insights into pathways of groundwater–surface water interaction can be applied to improve hydrological modelling in proglacial catchments throughout South America. Copyright © 2016 John Wiley & Sons, Ltd.     

On the Panday and Huyakorn surface–subsurface hydrology test case: analysis of internal flow dynamics

A number of previous studies using models of integrated surface‐subsurface hydrology have adopted the Panday and Huyakorn (P&H) tilted V‐catchment test case (Panday S, Huyakorn PS. 2004. A fully coupled spatially distributed model for evaluating surface/subsurface flow. Advances in Water Resources 27: 361–382) to show inter‐code comparability. The P&H test case is used to evaluate models that simulate a broad range of hydrological processes, and yet only the catchment outflow hydrograph has been presented as verification of the consistency between codes. Therefore, a more comprehensive evaluation of the surface‐subsurface hydrology of the P&H case is needed. This study explores the internal catchment functioning of the P&H case, using the popular catchment simulator MODHMS. The processes leading to streamflow generation in the model are illustrated, including separation of overland flow (OLF) and groundwater discharge to the stream. The results identify non‐physical flow processes due to the problem set‐up, and modifications to the P&H case are suggested that include changes to stream roughness and incision of the stream channel to overcome these shortcomings. A modified P&H case produced more plausible transfers between OLF and the stream, and an increased groundwater discharge to the stream (6·5% of streamflow in the modified case compared to 0·5% in the original case). Despite changes to internal flow processes, near‐identical outflow hydrographs were obtained, showing the importance of considering and comparing internal flow processes when using surface‐subsurface hydrology test cases to evaluate integrated hydrological simulators. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of a grid‐based version of the SWAT landscape model

Integrated river basin models should provide a spatially distributed representation of basin hydrology and transport processes to allow for spatially implementing specific management and conservation measures. To accomplish this, the Soil and Water Assessment Tool (SWAT) was modified by integrating a landscape routing model to simulate water flow across discretized routing units. This paper presents a grid‐based version of the SWAT landscape model that has been developed to enhance the spatial representation of hydrology and transport processes. The modified model uses a new flow separation index that considers topographic features and soil properties to capture channel and landscape flow processes related to specific landscape positions. The resulting model is spatially fully distributed and includes surface, lateral and groundwater fluxes in each grid cell of the watershed. Furthermore, it more closely represents the spatially heterogeneous distributed flow and transport processes in a watershed. The model was calibrated and validated for the Little River Watershed (LRW) near Tifton, Georgia (USA). Water balance simulations as well as the spatial distribution of surface runoff, subsurface flow and evapotranspiration are examined. Model results indicate that groundwater flow is the dominant landscape process in the LRW. Results are promising, and satisfactory output was obtained with the presented grid‐based SWAT landscape model. Nash–Sutcliffe model efficiencies for daily stream flow were 0.59 and 0.63 for calibration and validation periods, and the model reasonably simulates the impact of the landscape position on surface runoff, subsurface flow and evapotranspiration. Additional revision of the model will likely be necessary to adequately represent temporal variations of transport and flow processes in a watershed. Copyright © 2014 John Wiley & Sons, Ltd.     

Hydrologic comparison between a forested and a wetland/lake dominated watershed using SWAT

The Soil and Water Assessment Tool (SWAT) is a physically‐based hydrologic model developed for agricultural watersheds, which has been infrequently validated for forested watersheds, particularly those with deep overwinter snow accumulation and abundant lakes and wetlands. The goal of this study was to determine the applicability of SWAT for modelling streamflow in two watersheds of the Ontonagon River basin of northern Michigan which differ in proportion of wetland and lake area. The forest‐dominated East Branch watershed contains 17% wetland and lake area, whereas the wetland/lake‐dominated Middle Branch watershed contains 26% wetland and lake area. The specific objectives were to: (1) calibrate and validate SWAT models for the East Branch and Middle Branch watersheds to simulate monthly stream flow, and (2) compare the effects of wetland and lake abundance on the magnitude and timing of streamflow. Model calibration and validation was satisfactory, as determined by deviation of discharge D and Nash and Sutcliffe coefficient values E that compared simulated monthly mean discharge versus measured monthly mean discharge. Streamflow simulation discrepancies occurred during summer and fall months and dry years. Several snow melting parameters were found to be critical for the SWAT simulation: TIMP (snow temperature lag factor) and SMFMX and SMFMN (melting factors). Snow melting parameters were not transferable between adjacent watersheds. Differences in seasonal pattern of long‐term monthly streamflow were found, with the forest‐dominated watershed having a higher peak flow during April but a lower flow during the remainder of the year in comparison to the wetland and lake‐dominated watershed. The results suggested that a greater proportion of wetland and lake area increases the capacity of a watershed to impound surface runoff and to delay storm and snow melting events. Representation of wetlands and lakes in a watershed model is required to simulate monthly stream flow in a wetland/lake‐dominated watershed. Copyright © 2007 John Wiley & Sons, Ltd.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved process-based representation of stream solute transport in the soil and water assessment tools

Hydrological models have long been used to study the interactions between land, surface and groundwater systems, and to predict and manage water quantity and quality. The soil and water assessment tool (SWAT), a widely used hydrological model, can simulate various ecohydrological processes on land and subsequently route the water quality constituents through surface and subsurface waters. So far, in-stream solute transport algorithms of the SWAT model have only been minimally revised, even though it has been acknowledged that an improvement of in-stream process representation can contribute to better model performance with respect to water quality. In this study, we aim to incorporate a new and improved solute transport model into the SWAT model framework. The new process-based model was developed using in-stream process equations from two well established models—the One-dimensional Transport with Inflow and Storage model and the Enhanced Stream Water Quality Model. The modified SWAT model (Mir-SWAT) was tested for water quality predictions in a study watershed in Germany. Compared to the standard SWAT model, Mir-SWAT improved dissolved oxygen (DO) predictions by removing extreme low values of DO (<6 mg/L) simulated by SWAT. Phosphate concentration peaks were reduced during high flows and a better match of daily predicted and measured values was attained using the Mir-SWAT model (R² = 0.17, NSE = −0.65, RSR = 1.29 with SWAT; R² = 0.28, NSE = −0.04, RSR = 1.02 with Mir-SWAT). In addition, Mir-SWAT performed better than the SWAT model in terms of Chlorophyll-a content particularly during winter months, improving the NSE and RSR for monthly average Chl-a by 74 and 42%, respectively. With the new model improvements, we aim to increase confidence in the stream solute transport component of the model, improve the understanding of nutrient dynamics in the stream, and to extend the applicability of SWAT for reach-scale analysis and management.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the effects of land‐use change on water and salt delivery from a catchment affected by dryland salinity in south‐east Australia

A comprehensive framework for the assessment of water and salt balance for large catchments affected by dryland salinity is applied to the Boorowa River catchment (1550 km²), located in south‐eastern Australia. The framework comprised two models, each focusing on a different aspect and operating on a different scale. A quasi‐physical semi‐distributed model CATSALT was used to estimate runoff and salt fluxes from different source areas within the catchment. The effects of land use, climate, topography, soils and geology are included. A groundwater model FLOWTUBE was used to estimate the long‐term effects of land‐use change on groundwater discharge. Unlike conventional salinity studies that focus on groundwater alone, this study makes use of a new approach to explore surface and groundwater interactions with salt stores and the stream. Land‐use change scenarios based on increased perennial pasture and tree‐cover content of the vegetation, aimed at high leakage and saline discharge areas, are investigated. Likely downstream impacts of the reduction in flow and salt export are estimated. The water balance model was able to simulate both the daily observed stream flow and salt load at the catchment outlet for high and low flow conditions satisfactorily. Mean leakage rate of about 23·2 mm year^?1 under current land use for the Boorowa catchment was estimated. The corresponding mean runoff and salt export from the catchment were 89 382 ML year^?1 and 38 938 t year^?1, respectively. Investigation of various land‐use change scenarios indicates that changing annual pastures and cropping areas to perennial pastures is not likely to result in substantial improvement of water quality in the Boorowa River. A land‐use change of about 20% tree‐cover, specifically targeting high recharge and the saline discharge areas, would be needed to decrease stream salinity by 150 µS cm^?1 from its current level. Stream salinity reductions of about 20 µS cm^?1 in the main Lachlan River downstream of the confluence of the Boorowa River is predicted. The FLOWTUBE modelling within the Boorowa River catchment indicated that discharge areas under increased recharge conditions could re‐equilibrate in around 20 years for the catchment, and around 15 years for individual hillslopes. Copyright © 2004 John Wiley & Sons, Ltd.     

Modeling groundwater–surface water interactions using the Dupuit approximation

Global errors in head and/or discharge may be introduced when groundwater flow to a stream is modeled using the Dupuit approximation. We consider a simple case of steady groundwater flow in the vertical plane to a horizontal stream bed in direct connection with the aquifer, and compare solutions to the exact problem with Dupuit solutions where common representations of the stream are chosen. In all cases considered, adopting the Dupuit approximation introduces global errors into the mathematical model, and the magnitude of the errors depends on the regional flow conditions. This behavior makes calibration of a model difficult and limits the predictive abilities of the model under conditions of changed regional flow. The global errors and their dependence on flow conditions can be minimized, but not eliminated by treating the resistance of a fictitious leaky stream bed as an effective parameter.We propose an alternate Dupuit model of groundwater–surface water interaction and demonstrate, for the case considered, that adding a second effective parameter allows us to eliminate global errors in head and discharge, and eliminate the dependence of the effective values on the flow field. Explicit expressions are provided to evaluate the two effective properties. We propose that the results be used as a general guideline for modeling groundwater–surface water interaction at streams.     

Advantages of the analytic element method for the solution of groundwater management problems

In the simulation‐optimization approach, a coupled optimization and groundwater flow/transport model is used to solve groundwater management problems. The efficiency of the numerical method, which is used to simulate the groundwater flow, is one the major reason to obtain the best solution for a management problem. This study was carried out to examine the advantages of the analytic element method (AEM) in the simulation‐optimization approach, for the solution of groundwater management problems. For this study, the AEM and finite difference method (FDM) based flow models were developed and coupled with the particle swarm optimization (PSO)‐based optimization model. Furthermore, the AEM‐PSO and FDM‐PSO models developed were applied in hypothetical as well as real field conditions to address groundwater management problems and the results were compared. For the real field situation, the models developed were applied to the Dore River basin in France to minimize the installation and operational cost of new pumping wells taking the location and discharge of the pumping wells as decision variables. The constraints of the problem were identified with the help of stakeholders and water authority officials. The AEM flow model was developed to facilitate the management model particularly when at each iteration, the optimization model calls for a simulation model to calculate the values of groundwater heads. The results show that, at some points, the AEM‐PSO model is efficient in identifying the optimal location of wells and consequently results in optimal costs, sometimes difficult when using the FDM. Copyright © 2011 John Wiley & Sons, Ltd.     

Improving hydrologic predictions of a catchment model via assimilation of surface soil moisture

This paper examines the potential for improving Soil and Water Assessment Tool (SWAT) hydrologic predictions of root-zone soil moisture, evapotranspiration, and stream flow within the 341 km² Cobb Creek Watershed in southwestern Oklahoma through the assimilation of surface soil moisture observations using an Ensemble Kalman filter (EnKF). In a series of synthetic twin experiments assimilating surface soil moisture is shown to effectively update SWAT upper-layer soil moisture predictions and provide moderate improvement to lower layer soil moisture and evapotranspiration estimates. However, insufficient SWAT-predicted vertical coupling results in limited updating of deep soil moisture, regardless of the SWAT parameterization chosen for root-water extraction. Likewise, a real data assimilation experiment using ground-based soil moisture observations has only limited success in updating upper-layer soil moisture and is generally unsuccessful in enhancing SWAT stream flow predictions. Comparisons against ground-based observations suggest that SWAT significantly under-predicts the magnitude of vertical soil water coupling at the site, and this lack of coupling impedes the ability of the EnKF to effectively update deep soil moisture, groundwater flow and surface runoff. The failed attempt to improve stream flow prediction is also attributed to the inability of the EnKF to correct for existing biases in SWAT-predicted stream flow components.