Implementing a nonlinear groundwater module in the soil and water assessment tool (SWAT)

The Soil and Water Assessment Tool (SWAT) is widely used in modeling water quantity and quality. In the original SWAT, groundwater flow is calculated using a linear‐reservoir model, with outflow proportional to storage. However, observations show that this assumption is not always applicable; for example, macropores in Karst formations would seriously affect the groundwater behavior. A nonlinear groundwater algorithm was introduced in a new version of the SWAT model, called ISWAT. The Shenandoah Valley area in the Eastern U.S., which includes a number of geologic formations including Karst, was selected to test the modified ISWAT model. Parameter ESTimation (PEST) was coupled with ISWAT to auto‐calibrate the nonlinear parameter values. Ten years of record at 15 stream gauges were used to calibrate the model. The nonlinear ISWAT, statistically and visually, performed better in stream discharge estimation especially during baseflow recession and low‐flow periods. This indicated that the nonlinear algorithm can better represent groundwater behavior. The coupled ISWAT‐PEST approach can be used in future stream discharge simulation. Copyright © 2013 John Wiley & Sons, Ltd.     

Understanding salt mobilization from an irrigated catchment in south‐eastern Australia

Groundwater discharge from the Riverine Plains of the southern Murray‐Darling Basin is a major process contributing salt to the Murray River in Australia. In this study, data from an irrigated 60 000 ha catchment in the Riverine Plains were analysed to understand groundwater discharge into deeply incised drains, the process dominating salt mobilization from the catchment. We applied three integrated methodologies: classification and regression trees (CART), conceptual modelling and artificial neural networks (ANNs) to a comprehensive, spatially lumped, monthly data set from July 1975 to December 2004. Using CART analysis, it was shown that rainfall was the most important variable consistently explaining the salt load patterns at the catchment outlet. Using the conceptual model representing spatially lumped groundwater discharge into deeply incised drains, we demonstrated that salt mobilization from the study catchment can be well represented by a rainfall contribution, influenced by the hydraulic head in the deep regional aquifer and potential evapotranspiration. Using ANNs, it was confirmed that rainfall had a much higher impact on salt loads at the catchment outlet than irrigation water use. All these results demonstrate that under conditions similar to those experienced from 1975 to 2004, it is rainfall rather than irrigation water use that governs salt mobilization from the study catchment. Management of salt mobilization from irrigated catchments has traditionally focussed on the improvement of irrigation practices but it could be equally important to further understand the scope for management to control groundwater discharge in these irrigation areas. Copyright © 2010 John Wiley & Sons, Ltd.     

Response of hydrological processes to land‐cover and climate changes in Kejie watershed,south‐west China

Land‐cover/climate changes and their impacts on hydrological processes are of widespread concern and a great challenge to researchers and policy makers. Kejie Watershed in the Salween River Basin in Yunnan, south‐west China, has been reforested extensively during the past two decades. In terms of climate change, there has been a marked increase in temperature. The impact of these changes on hydrological processes required investigation: hence, this paper assesses aspects of changes in land cover and climate. The response of hydrological processes to land‐cover/climate changes was examined using the Soil and Water Assessment Tool (SWAT) and impacts of single factor, land‐use/climate change on hydrological processes were differentiated. Land‐cover maps revealed extensive reforestation at the expense of grassland, cropland, and barren land. A significant monotonic trend and noticeable changes had occurred in annual temperature over the long term. Long‐term changes in annual rainfall and streamflow were weak; and changes in monthly rainfall (May, June, July, and September) were apparent. Hydrological simulations showed that the impact of climate change on surface water, baseflow, and streamflow was offset by the impact of land‐cover change. Seasonal variation in streamflow was influenced by seasonal variation in rainfall. The earlier onset of monsoon and the variability of rainfall resulted in extreme monthly streamflow. Land‐cover change played a dominant role in mean annual values; seasonal variation in surface water and streamflow was influenced mainly by seasonal variation in rainfall; and land‐cover change played a regulating role in this. Surface water is more sensitive to land‐cover change and climate change: an increase in surface water in September and May due to increased rainfall was offset by a decrease in surface water due to land‐cover change. A decrease in baseflow caused by changes in rainfall and temperature was offset by an increase in baseflow due to land‐cover change. Copyright © 2009 John Wiley & Sons, Ltd.     

Multi-objective autocalibration of SWAT model for improved low flow performance for a small snowfed catchment

ABSTRACT

A reliable modelling framework needs to ensure that the model is simulating reality with limited uncertainty, thus enhancing its predictive ability. In the literature, hydrological model assessment using one or more metrics is reported to be inadequate when the river flow regime is required to be reproduced comprehensively. This research is aimed to: (a) calibrate the Soil and Water Assessment Tool (SWAT) based on the concept of multi-objective optimization by applying the Borg multi-objective evolutionary algorithm (MOEA); (b) apply hydrological signatures as objective functions; and (c) adopt a multi-metric approach for model evaluation. The SWAT model was coupled with a relatively newer and powerful Borg MOEA. The inclusion of hydrological signatures as objective functions along with the conventional statistical functions assisted in improving the performance for low flows by 135% in terms of volume efficiency and 65% for flow time series simulation.     

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.     

SWAT模型在斯里兰卡河流径流预测中的运用

本文运用SWAT模型和新安江模型对斯里兰卡卡鲁河流域上游地区日径流进行了预测.卡鲁河是斯里兰卡的第二大河,由于流域的降雨量很大,上游地区河流沿峡谷流下,中下游平原地区河床平坦.卡鲁河流域的洪水变的很正常.应用SWAT模型来对卡鲁河的日径流量进行预测,并同应用新安江模型所得到的结果做对比.研究表明,新安江模型要比SWAT (分布式水文模型)模型在卡鲁河日径流量预测上稍微好一些.实际上,或许数据质量不高或不恰当是部分原因,因为SWAT的输出成果严格取决于其输入的数据质量.此外,在斯里兰卡,许多人的日常用水是靠井水.当把流域看作一个整体,通常都是一个很大的范围,那样的话就不可能详尽的记录所有各个小规模的水利用,例如:小灌溉、小规模的家畜管理和工业水利用.这些水利用累积起来或许就很可观.这些数据的缺失对分布式水文模型在水平衡的应用有着独特的影响.但是概念水文模型(如新安江模型)可以根据实际情况在校正中调节它的参数,因为这些参数并没有实质的物理含义.因此,在流域特征和模型输入数据有限或不完整的情况下,概念水文模型比分布式水文模型更具优势.     

Morphological signatures of normal faulting in low‐gradient alluvial rivers in south‐eastern Louisiana,USA

We explore the fluvial response to faulting in three low‐gradient, sand‐bed rivers in south‐eastern Louisiana, USA, that flow across active normal faults from footwall (upstream) to hangingwall (downstream). We calculate sinuosity, migration rate and migration direction in order to identify anomalies spatially associated with fault scarps. In two of the rivers we model one‐dimensional steady water flow to identify anomalies in surface water slope, width‐to‐depth ratio, and shear stress. In each of these rivers there is one location where flow modeling suggests potential channel incision through the footwall, as indicated by relatively high surface water slopes and shear stress values. In one of these footwall locations, the river straightens and width‐to‐depth ratios decrease, likely contributing to higher surface water slopes and shear stress. This is in contrast to previous studies that have proposed increased sinuosity across fault footwalls and decreased sinuosity across hangingwalls. However, in two hangingwall locations we also observe relatively less sinuous channels. Other planform changes on the hangingwall include topographic steering of channels along and towards the fault and one example of an avulsion. The most notable anomaly in migration rate occurs on the hangingwall of a fault where a river has cut off a meander loop. Although fluvial response to faulting varies here, comparatively large and small channels exhibit similar responses. Further, Pleistocene fault slip rates are orders of magnitude lower than the channel migration rates, suggesting that faulting should not be a major influence on the fluvial evolution. Nonetheless, notable channel anomalies exist near faults, suggesting that recent fault slip rates are higher than Pleistocene rates, and/or that low‐gradient alluvial channels are more sensitive to faulting than previous studies have suggested. Rivers appear to be influenced by faulting in this setting, however background rates of meander loop cutoff may be just as influential as faulting. Copyright © 2015 John Wiley & Sons, Ltd.     

Development and application of a physically based landscape water balance in the SWAT model

Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamic morphology of small south‐eastern Mediterranean river mouths: a conceptual model

River mouths along the Israeli Mediterranean coast are characterized by a dynamic morphology as their channels migrate hundreds of meters along the coast. This study examines the dynamic morphology of seven such river mouths. It offers a conceptual model aimed at generalizing and describing their spatial and temporal morphological patterns, and the environmental factors influencing them. The study methodology comprised a detailed monitoring and mapping by GIS techniques, with quantitative data derived from historic aerial photographs, river discharge records, wave measurements, and a digital elevation model. These data were incorporated into a homogenous database and subsequently applied in the investigation of the morphological patterns of these mouths, and the analysis of their influencing factors. River mouths in this study occur in two distinctive topographic settings. In one setting (here termed barrier topography) the river mouth is deflected alongshore by a sandy barrier. In the second setting (termed funnel topography) the river mouth is confined to a funnel‐shaped topographic depression perpendicular to the coast. The behavior of river mouths in these two settings is quite distinctive. Barrier mouths usually migrate over larger distances, as they tend to deflect along a sand barrier and establish semi‐permanent channels along the dune toe. This enables the wide range migration of semi‐permanent channels over decades. Funnel topography mouths deflect over shorter distances and they rapidly migrate within the funnel boundaries. This study concludes that the topographic setting of the beach, a constant element in the temporal scale of this study, is the primary influencing factor on the morphology of the mouths studied. The influence of other factors on the morphology of these mouths differs in space and time and depends on the topographic settings.     

Application of semi‐distributed hydrological model for basin level water balance of the Ken basin of Central India

In the present study, a semi‐distributed hydrological model soil and water assessment tool (SWAT) has been employed for the Ken basin of Central India to predict the water balance. The entire basin was divided into ten sub basins comprising 107 hydrological response units on the basis of unique slope, soil and land cover classes using SWAT model. Sensitivity analysis of SWAT model was performed to examine the critical input variables of the study area. For Ken basin, curve number, available water capacity, soil depth, soil evaporation compensation factor and threshold depth of water in the shallow aquifer (GWQ_MN) were found to be the most sensitive parameters. Yearly and monthly calibration (1985–1996) and validation (1997–2009) were performed using the observed discharge data of the Banda site in the Ken basin. Performance evaluation of the model was carried out using coefficient of determination, Nash–Sutcliffe efficiency, root mean square error‐observations standard deviation ratio, percent bias and index of agreement criterion. It was found that SWAT model can be successfully applied for hydrological evaluation of the Ken basin, India. The water balance analysis was carried out to evaluate water balance of the Ken basin for 25 years (1985–2009). The water balance exhibited that the average annual rainfall in the Ken basin is about 1132 mm. In this, about 23% flows out as surface run‐off, 4% as groundwater flow and about 73% as evapotranspiration. Copyright © 2013 John Wiley & Sons, Ltd.     

Surface water quality in a sulfide mineral‐rich arid zone in North‐Central Chile: Learning from a complex past,addressing an uncertain future

This study presents an analysis of up to 30 years of hydrological variables and selected water quality parameters (pH, SO₄, Fe, Cu, and As) in the upper area of the Elqui River basin in North‐Central Chile. A correlation analysis determined statistically significant positive relationship for SO₄‐Cu, Fe‐As, and Fe‐Cu. In terms of historical behaviour, no statistically significant trends were detected for precipitation or temperature. In contrast, for flow, there is an overall decreasing pattern for the entire area of study, although only in one case this trend was statistically significant. Along with the aforementioned analysis, a characterization of the flow‐water quality relationships is considered for the time period analyzed. Although erratic behaviours were confirmed, a negative (i.e., inverse) flow‐concentration relationship was identified for SO₄, a positive (i.e., direct) relationship for Fe, and undefined relationships for As and Cu were obtained. From these analyses and based on previous studies on projections regarding climate change for the Andean region, and in particular for the upper Elqui zone, an estimation of the possible effects of the change in water regimes on water quality in the area of study is developed. It is likely that a decrease in surface flow, as a consequence of climate change could translate into improvements in water quality in terms of Fe and eventually As and Cu, but into an impairment in the case of SO₄. In any case, this is a complex situation that demands special attention in the face of industrial activities that could be developed in tributaries like the Claro River, which currently play an important role in depurating or diluting contaminants in the waters of the Elqui River. Finally, it should be noted that this study addresses an issue that goes beyond the local interest and could be used as a reference to compare other transitional environments containing sulphide ores or areas of hydrothermal alterations, which are considered to be highly vulnerable to climate change and variability.     

How to improve the representation of hydrological processes in SWAT for a lowland catchment – temporal analysis of parameter sensitivity and model performance

Model diagnostic analyses help to improve the understanding of hydrological processes and their representation in hydrological models. A detailed temporal analysis detects periods of poor model performance and model components with potential for model improvements, which cannot be found by analysing the whole discharge time series. In this study, we aim to improve the understanding of hydrological processes by investigating the temporal dynamics of parameter sensitivity and of model performance for the Soil and Water Assessment Tool model applied to the Treene lowland catchment in Northern Germany. The temporal analysis shows that the parameter sensitivity varies temporally with high sensitivity for three groundwater parameters (groundwater time delay, baseflow recession constant and aquifer fraction coefficient) and one evaporation parameter (soil evaporation compensation factor). Whereas the soil evaporation compensation factor dominates in baseflow and resaturation periods, groundwater time delay, baseflow recession constant and aquifer fraction coefficient are dominant in the peak and recession phases. The temporal analysis of model performance identifies three clusters with different model performances, which can be related to different phases of the hydrograph. The lowest performance, when comparing six performance measures, is detected for the baseflow cluster. A spatially distributed analysis for six hydrological stations within the Treene catchment shows similar results for all stations. The linkage of periods with poor model performance to the dominant model components in these phases and with the related hydrological processes shows that the groundwater module has the highest potential for improvement. This temporal diagnostic analysis enhances the understanding of the Soil and Water Assessment Tool model and of the dominant hydrological processes in the lowland catchment. Copyright © 2013 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.     

Assessment of surface water resources and evapotranspiration in the Haihe River basin of China using SWAT model

Quantitative assessment of surface water resources (SWRs) and evapotranspiration (ET) is essential and significant for reasonably planning and managing water resources in the Haihe River basin which is facing severe water shortage. In this study, a distributed hydrological model of the Haihe River basin was constructed using the Soil and Water Assessment Tool, well considering the reservoirs and agricultural management practices for reasonable simulation. The crop parameters were independently calibrated with the observed crop data at six experimental stations. Then, sensitivity ranks of hydrological parameters were analysed, which suggested the important parameters used for calibration. The model was successfully calibrated using the monthly observed data of discharge in around 1970–1991 and actual ET (ET_a) in 2002–2004 for the mountainous area and Haihe plain, respectively. Meanwhile, good agreements between the simulated and statistical crop yields in 1985–2005 further verified the model's appropriateness. Finally, the calibrated model was used to assess SWRs and ET_a in time and space during 1961–2005. Results showed that the average annual natural SWRs and the ET_a were about 17.5 billion cubic metre and 542 mm, respectively, both with a slight downward trend. The spatial distributions of both SWRs and ET_a were significantly impacted by variations of precipitation and land use. Moreover, the reservoir in operation was the main factor for the noticeable decline of actual SWRs. In the Haihe plain, the ET_a with irrigation was increased by 46% compared with that under rainfed conditions. In addition, this study identified the regions with potential to improve the irrigation effects on water use. Copyright © 2012 John Wiley & Sons, Ltd.     

Meta‐analysis of flow modeling performances—to build a matching system between catchment complexity and model types

Hydrological models play a significant role in modelling river flow for decision making support in water resource management. In the past decades, many researchers have made a great deal of efforts in calibrating and validating various models, with each study being focused on one or two models. As a result, there is a lack of comparative analysis on the performance of those models to guide hydrologists to choose appropriate models for the individual climate and physical conditions. This paper describes a two‐level meta‐analysis to develop a matching system between catchment complexity (based on catchment significant features (CSFs)) and model types. The intention is to use the available CSF information for choosing the most suitable model type for a given catchment. In this study, the CSFs include the elements of climate, soil type, land cover and catchment scale. Specific choices of model types in small and medium catchments are further explored with all CSF information obtained. In particular, it is interesting to find that semi‐distributed models are the most suitable model type for catchments with the area over 3000 km², regardless of other CSFs. The potential methodology for expanding the matching system between catchment complexity and model complexity is discussed. Copyright © 2014 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.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of AnnAGNPS and SWAT model simulation results in USDA‐CEAP agricultural watersheds in south‐central Kansas

This study was conducted under the USDA‐Conservation Effects Assessment Project (CEAP) in the Cheney Lake watershed in south‐central Kansas. The Cheney Lake watershed has been identified as ‘impaired waters’ under Section 303(d) of the Federal Clean Water Act for sediments and total phosphorus. The USDA‐CEAP seeks to quantify environmental benefits of conservation programmes on water quality by monitoring and modelling. Two of the most widely used USDA watershed‐scale models are Annualized AGricultural Non‐Point Source (AnnAGNPS) and Soil and Water Assessment Tool (SWAT). The objectives of this study were to compare hydrology, sediment, and total phosphorus simulation results from AnnAGNPS and SWAT in separate calibration and validation watersheds. Models were calibrated in Red Rock Creek watershed and validated in Goose Creek watershed, both sub‐watersheds of the Cheney Lake watershed. Forty‐five months (January 1997 to September 2000) of monthly measured flow and water quality data were used to evaluate the two models. Both models generally provided from fair to very good correlation and model efficiency for simulating surface runoff and sediment yield during calibration and validation (correlation coefficient; R², from 0·50 to 0·89, Nash Sutcliffe efficiency index, E, from 0·47 to 0·73, root mean square error, RMSE, from 0·25 to 0·45 m³ s^?1 for flow, from 158 to 312 Mg for sediment yield). Total phosphorus predictions from calibration and validation of SWAT indicated good correlation and model efficiency (R² from 0·60 to 0·70, E from 0·63 to 0·68) while total phosphorus predictions from validation of AnnAGNPS were from unsatisfactory to very good (R² from 0·60 to 0·77, E from ? 2·38 to 0·32). The root mean square error–observations standard deviation ratio (RSR) was estimated as excellent (from 0·08 to 0·25) for the all model simulated parameters during the calibration and validation study. The percentage bias (PBIAS) of the model simulated parameters varied from unsatisfactory to excellent (from 128 to 3). This study determined SWAT to be the most appropriate model for this watershed based on calibration and validation results. Copyright © 2008 John Wiley & Sons, Ltd.