Predicting phosphorus dynamics in complex terrains using a variable source area hydrology model

Phosphorus (P) loss from agricultural watersheds has long been a critical water quality problem, the control of which has been the focus of considerable research and investment. Preventing P loss depends on accurately representing the hydrological and chemical processes governing P mobilization and transport. The Soil and Water Assessment Tool (SWAT) is a watershed model commonly used to predict run‐off and non‐point source pollution transport. SWAT simulates run‐off employing either the curve number (CN) or the Green and Ampt methods, both assume infiltration‐excess run‐off, although shallow soils underlain by a restricting layer commonly generate saturation‐excess run‐off from variable source areas (VSA). In this study, we compared traditional SWAT with a re‐conceptualized version, SWAT‐VSA, that represents VSA hydrology, in a complex agricultural watershed in east central Pennsylvania. The objectives of this research were to provide further evidence of SWAT‐VSA's integrated and distributed predictive capabilities against measured surface run‐off and stream P loads and to highlight the model's ability to drive sub‐field management of P. Thus, we relied on a detailed field management database to parameterize the models. SWAT and SWAT‐VSA predicted discharge similarly well (daily Nash–Sutcliffe efficiencies of 0.61 and 0.66, respectively), but SWAT‐VSA outperformed SWAT in predicting P export from the watershed. SWAT estimated lower P loss (0.0–0.25 kg ha^?1) from agricultural fields than SWAT‐VSA (0.0–1.0+ kg ha^?1), which also identified critical source areas – those areas generating large run‐off and P losses at the sub‐field level. These results support the use of SWAT‐VSA in predicting watershed‐scale P losses and identifying critical source areas of P loss in landscapes with VSA hydrology. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

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.     

Assessment of runoff and sediment yield in the Miyun Reservoir catchment by using SWAT model

The Soil and Water Assessment Tool (SWAT) was used to simulate the transport of runoff and sediment into the Miyun Reservoir, Beijing in this study. The main objective was to validate the performance of SWAT and the feasibility of using this model as a simulator of runoff and sediment transport processes at a catchment scale in arid and semi‐arid area in North China, and related processes affecting water quantity and soil erosion in the catchment were simulated. The investigation was conducted using a 6‐year historical streamflow and sediment record from 1986 to 1991; the data from 1986 to 1988 was used for calibration and that from 1989 to 1991 for validation. The SWAT generally performs well and could accurately simulate both daily and monthly runoff and sediment yield. The simulated daily and monthly runoff matched the observed values satisfactorily, with a Nash‐Sutcliffe coefficient of greater than 0·6, 0·9 and a coefficient of determination 0·75, 0·9 at two outlet stations (Xiahui and Zhangjiafen stations) during calibration. These values were 0·6, 0·85 and 0·6, 0·9 during validation. For sediment simulation, the efficiency is lower than that for runoff. Even so, the Nash‐Sutcliffe coefficient and coefficient of determination were greater than 0·48 and 0·6 for monthly sediment yield during calibration, and these values were greater than 0·84 and 0·95 during validation. Sensitivity analysis shows that sensitive parameters for the simulation of discharge and sediment yield include curve number, base flow alpha factor, soil evaporation compensation factor, soil available water capacity, soil profile depth, surface flow lag time and channel re‐entrained linear parameter, etc. Copyright © 2009 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Use of a parsimonious rainfall–run‐off model for predicting hydrological response in ungauged basins

A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.     

Temporally weighted average curve number method for daily runoff simulation

The modified Soil Conservation Service curve number (CN) method is widely used in long‐term continuous models to predict daily surface runoff. However, it has been shown that this method gives poor results in reproducing peak flows in high rainfall periods. This is because there is an inaccuracy stemming from the model algorithm as it adjusts the daily runoff curve number as a function of soil moisture content at the end of the previous day. This paper proposes an alternative daily based curve number technique that can provide better prediction of daily runoff during the high flow season. The proposed method uses the temporally weighted average curve number (TWA‐CN) to estimate daily surface runoff, while considering the effect of rainfall during a given day as well as the antecedent soil moisture condition. To test the applicability of the TWA‐CN method, it was incorporated with the long‐term, continuous simulation watershed models SWAT and SWAT‐G. Simulations were conducted for the Miho River watershed located in the middle of South Korea. The graphical displays and statistics of the determination coefficient (R²) and the Nash–Sutcliffe model efficiency (NSE) of the observed and simulated daily runoff indicated that the modified SWAT with the TWA‐CN method may provide better runoff prediction (R² = 0·837, NSE = 0·833) than the original SWAT (R² = 0·815, NSE = 0·824). Likewise, the determination coefficient (R² = 0·816) and the Nash–Sutcliffe efficiency (NSE = 0·834) for the modified SWAT‐G are also higher than the original version (R² = 0·782, NSE = 0·825). It is expected that the improved capability in predicting surface runoff using the suggested CN estimate method will provide a sound contribution to the accurate simulations of water yield. Copyright © 2008 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.     

Multi‐variable and multi‐site calibration and validation of SWAT in a large mountainous catchment with high spatial variability

Many methods developed for calibration and validation of physically based distributed hydrological models are time consuming and computationally intensive. Only a small set of input parameters can be optimized, and the optimization often results in unrealistic values. In this study we adopted a multi‐variable and multi‐site approach to calibration and validation of the Soil Water Assessment Tool (SWAT) model for the Motueka catchment, making use of extensive field measurements. Not only were a number of hydrological processes (model components) in a catchment evaluated, but also a number of subcatchments were used in the calibration. The internal variables used were PET, annual water yield, daily streamflow, baseflow, and soil moisture. The study was conducted using an 11‐year historical flow record (1990–2000); 1990–94 was used for calibration and 1995–2000 for validation. SWAT generally predicted well the PET, water yield and daily streamflow. The predicted daily streamflow matched the observed values, with a Nash–Sutcliffe coefficient of 0·78 during calibration and 0·72 during validation. However, values for subcatchments ranged from 0·31 to 0·67 during calibration, and 0·36 to 0·52 during validation. The predicted soil moisture remained wet compared with the measurement. About 50% of the extra soil water storage predicted by the model can be ascribed to overprediction of precipitation; the remaining 50% discrepancy was likely to be a result of poor representation of soil properties. Hydrological compensations in the modelling results are derived from water balances in the various pathways and storage (evaporation, streamflow, surface runoff, soil moisture and groundwater) and the contributions to streamflow from different geographic areas (hill slopes, variable source areas, sub‐basins, and subcatchments). The use of an integrated multi‐variable and multi‐site method improved the model calibration and validation and highlighted the areas and hydrological processes requiring greater calibration effort. Copyright © 2005 John Wiley & Sons, Ltd.     

Parameter identification and uncertainty analysis for simulating streamflow in a river basin of Eastern India

Uncertainty is inherent in modelling studies. However, the quantification of uncertainties associated with a model is a challenging task, and hence, such studies are somewhat limited. As distributed or semi‐distributed hydrological models are being increasingly used these days to simulate hydrological processes, it is vital that these models should be equipped with robust calibration and uncertainty analysis techniques. The goal of the present study was to calibrate and validate the Soil and Water Assessment Tool (SWAT) model for simulating streamflow in a river basin of Eastern India, and to evaluate the performance of salient optimization techniques in quantifying uncertainties. The SWAT model for the study basin was developed and calibrated using Parameter Solution (ParaSol), Sequential Uncertainty Fitting Algorithm (SUFI‐2) and Generalized Likelihood Uncertainty Estimation (GLUE) optimization techniques. The daily observed streamflow data from 1998 to 2003 were used for model calibration, and those for 2004–2005 were used for model validation. Modelling results indicated that all the three techniques invariably yield better results for the monthly time step than for the daily time step during both calibration and validation. The model performances for the daily streamflow simulation using ParaSol and SUFI‐2 during calibration are reasonably good with a Nash–Sutcliffe efficiency and mean absolute error (MAE) of 0.88 and 9.70 m³/s for ParaSol, and 0.86 and 10.07 m³/s for SUFI‐2, respectively. The simulation results of GLUE revealed that the model simulates daily streamflow during calibration with the highest accuracy in the case of GLUE (R² = 0.88, MAE = 9.56 m³/s and root mean square error = 19.70 m³/s). The results of uncertainty analyses by SUFI‐2 and GLUE were compared in terms of parameter uncertainty. It was found that SUFI‐2 is capable of estimating uncertainties in complex hydrological models like SWAT, but it warrants sound knowledge of the parameters and their effects on the model output. On the other hand, GLUE predicts more reliable uncertainty ranges (R‐factor = 0.52 for daily calibration and 0.48 for validation) compared to SUFI‐2 (R‐factor = 0.59 for daily calibration and 0.55 for validation), though it is computationally demanding. Although both SUFI‐2 and GLUE appear to be promising techniques for the uncertainty analysis of modelling results, more and more studies in this direction are required under varying agro‐climatic conditions for assessing their generic capability. Copyright © 2015 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.     

Implementation of an automatic calibration procedure for HYDROTEL based on prior OAT sensitivity and complementary identifiability analysis

Efficiency of hydrological models mostly depends on the quality of the calibration performed prior to use. In this paper, an automatic calibration framework for the distributed hydrological model HYDROTEL is proposed. The calibration procedure was performed for three watersheds characterized with different hydroclimatological conditions: the Sassandra located in Ivory Coast, Africa, and the Montmorency and Beaurivage watersheds located in Quebec (Canada). Results of one‐a‐time (OAT) sensitivity analysis showed that the order of the most sensitive parameters differs for each watershed. Thus, the sensitivity depends on the hydroclimatic and physiographic characteristics of the watersheds. Co‐linearity indices showed that all model parameters were identifiable, that is, none of the studied parameters could be explained by a combination of the other parameters. Following these findings, an automatic calibration was run. Results indicated there was good agreement between simulated and measured streamflows at the outlet of each watershed; Nash–Sutcliffe efficiency (NSE) ranging between 0.77 and 0.92 and R² ranging from 0.87 to 0.97. When comparing NSE and R² values obtained using a process‐oriented, multiple‐objective, manual calibration strategy, a slight increase in model efficiency was reached with the automatic calibration procedure (4.15% for NSE and 2.95% for R²) improving predictions of peak flows for the Montmorency and Beaurivage watersheds (temperate climate conditions) and flows beyond the rainfall season in the Sassandra watershed. The proposed automatic calibration procedure introduced in this paper may be applied to other distributed hydrological model. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluating the spatiotemporal variations of water budget across China over 1951–2006 using IBIS model

The Integrated Biosphere Simulator is used to evaluate the spatial and temporal patterns of the crucial hydrological variables [run‐off and actual evapotranspiration (AET)] of the water balance across China for the period 1951–2006 including a precipitation analysis. Results suggest three major findings. First, simulated run‐off captured 85% of the spatial variability and 80% of the temporal variability for 85 hydrological gauges across China. The mean relative errors were within 20% for 66% of the studied stations and within 30% for 86% of the stations. The Nash–Sutcliffe coefficients indicated that the quantity pattern of run‐off was also captured acceptably except for some watersheds in southwestern and northwestern China. The possible reasons for underestimation of run‐off in the Tibetan plateau include underestimation of precipitation and uncertainties in other meteorological data due to complex topography, and simplified representations of the soil depth attribute and snow processes in the model. Second, simulated AET matched reasonably with estimated values calculated as the residual of precipitation and run‐off for watersheds controlled by the hydrological gauges. Finally, trend analysis based on the Mann–Kendall method indicated that significant increasing and decreasing patterns in precipitation appeared in the northwest part of China and the Yellow River region, respectively. Significant increasing and decreasing trends in AET were detected in the Southwest region and the Yangtze River region, respectively. In addition, the Southwest region, northern China (including the Heilongjiang, Liaohe, and Haihe Basins), and the Yellow River Basin showed significant decreasing trends in run‐off, and the Zhemin hydrological region showed a significant increasing trend. Copyright © 2009 John Wiley & Sons, Ltd.     

Season‐based rainfall–runoff modelling using the probability‐distributed model (PDM) for large basins in southeastern Brazil

Southeastern Brazil is characterized by seasonal rainfall variability. This can have a great social, economic, and environmental impact due to both excessive and deficient water availability. During 2014 and 2015, the region experienced one of the most severe droughts since 1960. The resulting water crisis has seriously affected water supply to the metropolitan region of São Paulo and hydroelectric power generation throughout the entire country. This research considered the upstream basins of the southeastern Brazilian reservoirs Cantareira (2,279 km²; water supply) and Emborcação (29,076 km²), Três Marias (51,576 km²), Furnas (52,197 km²), and Mascarenhas (71,649 km²; hydropower) for hydrological modelling. It made the first attempt at configuring a season‐based probability‐distributed model (PDM‐CEMADEN) for simulating different hydrological processes during wet and dry seasons. The model successfully reproduced the intra‐annual and interannual variability of the upstream inflows during 1985–2015. The performance of the model was very satisfactory not only during the wet, dry, and transitional seasons separately but also during the whole period. The best performance was obtained for the upstream basin of Furnas, as it had the highest quality daily precipitation and potential evapotranspiration data. The Nash–Sutcliffe efficiency and logarithmic Nash–Sutcliffe efficiency were 0.92 and 0.93 for the calibration period 1984–2001, 0.87 and 0.88 for the validation period 2001–2010, and 0.93 and 0.90 for the validation period 2010–2015, respectively. Results indicated that during the wet season, the upstream basins have a larger capacity and variation of soil water storage, a larger soil water conductivity, and quicker surface water flow than during the dry season. The added complexity of configuring a season‐based PDM‐CEMADEN relative to the traditional model is well justified by its capacity to better reproduce initial conditions for hydrological forecasting and prediction. The PDM‐CEMADEN is a simple, efficient, and easy‐to‐use model, and it will facilitate early decision making and implement adaptation measures relating to disaster prevention for reservoirs with large‐sized upstream basins.     

SWAT model application and prediction uncertainty analysis in the Lake Tana Basin,Ethiopia

Lake Tana Basin is of significant importance to Ethiopia concerning water resources aspects and the ecological balance of the area. Many years of mismanagement, wetland losses due to urban encroachment and population growth, and droughts are causing its rapid deterioration. The main objective of this study was to assess the performance and applicability of the soil water assessment tool (SWAT) model for prediction of streamflow in the Lake Tana Basin, so that the influence of topography, land use, soil and climatic condition on the hydrology of Lake Tana Basin can be well examined. The physically based SWAT model was calibrated and validated for four tributaries of Lake Tana. Sequential uncertainty fitting (SUFI‐2), parameter solution (ParaSol) and generalized likelihood uncertainty estimation (GLUE) calibration and uncertainty analysis methods were compared and used for the set‐up of the SWAT model. The model evaluation statistics for streamflows prediction shows that there is a good agreement between the measured and simulated flows that was verified by coefficients of determination and Nash Sutcliffe efficiency greater than 0·5. The hydrological water balance analysis of the basin indicated that baseflow is an important component of the total discharge within the study area that contributes more than the surface runoff. More than 60% of losses in the watershed are through evapotranspiration. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling potential hydrological impact of abandoned underground mines in the Monday Creek Watershed,Ohio

Abandoned underground mines (AUM) have caused dramatic environmental effects that are closely linked to regional sustainability. This paper explores the potential hydrological impact of AUM in the Monday Creek Watershed, a typically mined area in Appalachian region, using the Soil and Water Assessment Tool (SWAT 2005) model and Sequential Uncertainty Fitting (SUFI‐2), calibrated at both the global and local scales. The locally calibrated model better incorporates those key parameters relevant to AUM for specific sub‐basins and hydrologic response units. Data from the years 2003–2004 were used for calibration and 2005–2006 for validation. The results were quite satisfactory; both the coefficient of determination (R²) and the Nash–Sutcliffe efficiency statistic were over 0.80. The potential influences of AUM were assessed by modelling an alternative scenario assuming no AUM for the period 2003–2009. Results show that the hydrological process of lateral subsurface flow plays a dominant role in linking AUM to overall watershed hydrology. The potential hydrological impact of AUM is an increased annual lateral flow of 82.1%, and a decrease in annual surface flow by 15%, leading to an increase of 16.9% in annual water yield for the Monday Creek Watershed. The seasonal fluctuation of water yield has a similar trend to lateral flow, decreasing from March to August and increasing from August to January. Higher volume, higher flow peaks and higher recession constants characterized the hydrograph of daily streamflow from AUM. The results indicate that more emphasis should be put on lateral flow for further study of acid mine drainage and flooding control in those watersheds with AUM. Copyright © 2012 John Wiley & Sons, Ltd.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Using high resolution tracer data to constrain water storage,flux and age estimates in a spatially distributed rainfall‐runoff model

Models simulating stream flow and conservative tracers can provide a representation of flow paths, storage distributions and mixing processes that is advantageous for many predictive purposes. Compared with models that only simulate stream flow, tracer data can be used to investigate the internal consistency of model behaviour and to gain insight into model performance. Here, we examine the strengths and weaknesses of a data‐driven, spatially distributed tracer‐aided rainfall‐runoff model. The model structure allowed us to assess the influence of landscape characteristics on the routing and mixing of water and tracers. The model was applied to a site in the Scottish Highlands with a unique tracer data set; ~4 years of daily isotope ratios in stream water and precipitation were available, as well as 2 years of weekly soil and ground water isotopes. The model structure was based on an empirically based, lumped tracer‐aided model previously developed for the catchment. The best model runs were selected from Monte Carlo simulations based on dual calibration criteria using objective functions for both stream isotopes and discharge at the outlet. Model performance for these criteria was reasonable (Nash–Sutcliffe efficiencies for discharge and isotope ratios were ~0.4–0.6). The model could generally reproduce the variable isotope signals in the soils of the steeper hill slopes where storage was low, and damped isotope responses in valley bottom cells with high storage. The model also allowed us to estimate the age distributions of internal stores, water fluxes and stream flow. Average stream water age was ~1.6 years, integrating older groundwater in the valley bottom and dynamic younger soil waters. By tracking water ages and simulating isotopes, the model captured the changes in connectivity driven by distributed storage dynamics. This has substantially improved the representation of spatio‐temporal process dynamics and gives a more robust framework for projecting environmental change impacts. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Spatial Quantification of Groundwater Abstraction in the Irrigated Indus Basin

Groundwater abstraction and depletion were assessed at a 1‐km resolution in the irrigated areas of the Indus Basin using remotely sensed evapotranspiration (ET) and precipitation; a process‐based hydrological model and spatial information on canal water supplies. A calibrated Soil and Water Assessment Tool (SWAT) model was used to derive total annual irrigation applied in the irrigated areas of the basin during the year 2007. The SWAT model was parameterized by station corrected precipitation data (R) from the Tropical Rainfall Monitoring Mission, land use, soil type, and outlet locations. The model was calibrated using a new approach based on spatially distributed ET fields derived from different satellite sensors. The calibration results were satisfactory and strong improvements were obtained in the Nash‐Sutcliffe criterion (0.52 to 0.93), bias (?17.3% to ?0.4%), and the Pearson correlation coefficient (0.78 to 0.93). Satellite information on R and ET was then combined with model results of surface runoff, drainage, and percolation to derive groundwater abstraction and depletion at a nominal resolution of 1 km. It was estimated that in 2007, 68 km³ (262 mm) of groundwater was abstracted in the Indus Basin while 31 km³ (121 mm) was depleted. The mean error was 41 mm/year and 62 mm/year at 50% and 70% probability of exceedance, respectively. Pakistani and Indian Punjab and Haryana were the most vulnerable areas to groundwater depletion and strong measures are required to maintain aquifer sustainability.