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.     

Hydrological modelling of a small watershed using generated rainfall in the soil and water assessment tool model

An adequately tested soil and water assessment tool (SWAT) model was applied to the runoff and sediment yield of a small agricultural watershed in eastern India using generated rainfall. The capability of the model for generating rainfall was evaluated for a period of 18 years (1981–1998). The watershed and subwatershed boundaries, drainage networks, slope, soil series and texture maps were generated using a geographical information system (GIS). A supervised classification method was used for land‐use/cover classification from satellite imageries. Model simulated monthly rainfall for the period of 18 years was compared with observations. Simulated monthly rainfall, runoff and sediment yield values for the monsoon season of 8 years (1991–1998) were also compared with their observed values. In general monthly average rainfall predicted by the model was in close agreement with the observed monthly average values. Also, simulated monthly average values of surface runoff and sediment yield using generated rainfall compared well with observed values during the monsoon season of the years 1991–1998. Results of this study revealed that the SWAT model can generate monthly average rainfall satisfactorily and thereby can produce monthly average values of surface runoff and sediment yield close to the observed values. Therefore, it can be concluded that the SWAT model could be used for developing a multiple year management plan for the critical erosion prone areas of a small watershed. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Assessing impacts of riparian buffer zones on sediment and nutrient loadings into streams at watershed scale using an integrated REMM‐SWAT model

Riparian buffer zones in agriculture dominated watersheds play important roles in reducing nonpoint source pollution into aquatic ecosystems and are widely used as a Best Management Practice. Assessment of the effectiveness of riparian buffer zones by modeling method is widely used for watershed management as field measurement‐based assessment is difficult and expensive. The integration of Riparian Ecosystem Management Model (REMM) and Soil and Water Assessment Tool (SWAT) has been developed to simulate the effect of nonpoint source pollution reduction by riparian buffer zones at subbasin scale. However, there are problems in using the integrated model at subbasin scale, as the size of subbasin partition could affect the pollutant reduction rate by riparian buffers. In this study, we partitioned a large watershed with size of 1331 ha into sub‐watersheds with sizes of 666, 333, 166, 83, 51, and 29 ha, and then compared the different simulation results. We found that the modeling could yield more convergent results when the sub‐watersheds were partitioned into suitable size. In the studied area, the suitable sub‐watershed size was less than about 166 ha for runoff and nitrogen and 83 ha for sediment and phosphorus. Among the eight sub‐watersheds (partitioned based on the size of 166 ha), results showed that the effects of riparian buffers on runoff and nutrient loading varied drastically. The reduction rate varied from 0.26% to 30.13% for runoff, 29.4% to 74.07% for sediment, 9.61% to 57.85% for nitrogen, and 18.61% to 68.12% for phosphorus, respectively.     

Application of the WEPP model for prioritization and evaluation of best management practices in an Indian watershed

The pre‐calibrated and validated physically based watershed model, water erosion prediction project (WEPP) was used as a modelling tool for the identification of critical watersheds and evaluation of best management practices for a small hilly watershed (Karso) of India. The land use/cover of the study area was generated using IRS‐1C LISS‐III (linear imaging self scanner) satellite data. The watershed and sub‐watershed boundaries, drainage, slope and soil map of the study area were generated using ARC/INFO geographic information system (GIS). The WEPP model was finally applied to the Karso watershed which lies within Damodar Barakar catchment of India to identify the critical sub‐watersheds on the basis of their simulated average annual sediment yields. Priorities were fixed on the basis of ranks assigned to each critical sub‐watershed based on the susceptibility to erosion. The sub‐watershed having the highest sediment yield was assigned a priority number 1, the next highest value was assigned a priority number 2, and so on. Subsequently, the model was used for evaluating the effectiveness of best management practices (crop and tillage) for conservation of soil for all the sub‐watersheds. On the basis of this study, it is realized that cash crops like soyabean should be encouraged in the upland portion of the sub‐watersheds, and the existing tillage practice (country plough/mould board plough) may be replaced by a field cultivation system for conservation of soil and water in the sub‐watersheds. Copyright © 2009 John Wiley & Sons, Ltd.     

Multi‐site evaluation of hydrology component of SWAT in the coastal plain of southwest Georgia

Simulation of watershed scale hydrologic and water quality processes is important for watershed assessments. Proper characterization of the accuracy of these simulations, particularly in cases with limited observed data, is critical. The Soil & Water Assessment Tool (SWAT) is frequently used for watershed scale simulation. The accuracy of the model was assessed by extrapolating calibration results from a well studied Coastal Plain watershed in Southwest Georgia, USA, to watersheds within the same geographic region without further calibration. SWAT was calibrated and validated on a 16.7‐km² subwatershed within the Little River Experimental Watershed by varying six model parameters. The optimized parameter set was then applied to a watershed of similar land use and soils, a smaller watershed with different land use and soils and three larger watersheds within the same drainage system without further calibration. Simulation results with percent bias (PB) ±15% ≤ PB < ±25% and Nash–Sutcliffe efficiency (NSE) 0.50 < NSE ≤ 0.65 were considered to be satisfactory, whereas those with PB < ±10% and 0.75 < NSE ≤ 1.00 were considered very good. With these criteria, simulation results for the five non‐calibration watersheds were satisfactory to very good. Differences across watersheds were attributed to differences in soils, land use, and surficial aquifer characteristics. These results indicate that SWAT can be a useful tool for predicting streamflow for ungauged watersheds with similar physical characteristics to the calibration watershed studied here and provide an indication of the accuracy of hydrologic simulations for ungauged watersheds. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the impact of agroforestry on hydrology of Mara River Basin in East Africa

Land‐use change is one of the main drivers of watershed hydrology change. The effect of forestry related land‐use changes (e.g. afforestation, deforestation, agroforestry) on water fluxes depends on climate, watershed characteristics and spatial scale. The Soil and Water Assessment Tool (SWAT) model was calibrated, validated and used to simulate the impact of agroforestry on the water balance in the Mara River Basin (MRB) in East Africa. Model performance was assessed by Nash–Sutcliffe Efficiency (NSE) and Kling–Gupta Efficiency (KGE). The NSE (and KGE) values for calibration and validation were: 0.77 (0.88) and 0.74 (0.85) for the Nyangores sub‐watershed, and 0.78 (0.89) and 0.79 (0.63) for the entire MRB. It was found that agroforestry in the watershed would generally reduce surface runoff, mainly because of enhanced infiltration. However, it would also increase evapotranspiration and consequently reduce baseflow and overall water yield, which was attributed to increased water use by trees. Spatial scale was found to have a significant effect on water balance; the impact of agroforestry was higher at the smaller headwater catchment (Nyangores) than for the larger watershed (entire MRB). However, the rate of change in water yield with an increase in area under agroforestry was different for the two and could be attributed to the spatial variability of climate within the MRB. Our results suggest that direct extrapolation of the findings from a small sub‐catchment to a larger watershed may not always be accurate. These findings could guide watershed managers on the level of trade‐offs that might occur between reduced water yields and other benefits (e.g. soil erosion control, improved soil productivity) offered by agroforestry. Copyright © 2016 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.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 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.     

A simple concept for calibrating runoff thresholds in quasi‐distributed variable source area watershed models

Most semi‐distributed watershed water quality models divide the watershed into hydrologic response units (HRU) with no flow among them. This is problematic when watersheds are delineated to include variable source areas (VSAs) because it is the lateral flows from upslope areas to downslope areas that generate VSAs. Although hydrologic modellers have often successfully calibrated these types of models, there can still be considerable uncertainty in model results. In this paper, a topographic‐index‐based method is described and tested to distribute effective soil water holding capacity among HRUs, which can be subsequently adjusted using the watershed baseflow coefficient. The method is tested using a version of the Soil and Water Assessment Tool (SWAT) model that simulates VSA runoff and is applied to two watersheds: a New York State (NYS) watershed, and one in the head waters of the Blue Nile Basin (BNB) in Ethiopia. Daily streamflow predicted using effective soil water storage capacities based only on the topographic index were reassuringly accurate in both the NYS watershed (daily Nash Sutcliffe (E) = 0·73) and in the BNB (E = 0·70). Using the baseflow coefficient to adjust the effective soil water storage capacity only slightly improved streamflow predictions in NYS (E = 0·75) but substantially improved the BNB predictions (E = 0·80). By comparison, the standard SWAT model, which uses the traditional look‐up tables to determine a runoff curve number, performed considerably less accurately in un‐calibrated form (E = 0·51 for NYS and E = 0·45 for BNB), but improved substantially when explicitly calibrated to streamflow measurements (E = 0·76 for NYS and E = 0·67 for the BNB). The calibration method presented here provides a parsimonious, systematic approach to using established models in VSA watersheds that reduces the ambiguity inherent in model calibration. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulating daily soil water under foothills fescue grazing with the soil and water assessment tool model (Alberta,Canada)

Grazing is common in the foothills fescue grasslands and may influence the seasonal soil‐water patterns, which in turn determine range productivity. Hydrological modelling using the soil and water assessment tool (SWAT) is becoming widely adopted throughout North America especially for simulation of stream flow and runoff in small and large basins. Although applications of the SWAT model have been wide, little attention has been paid to the model's ability to simulate soil‐water patterns in small watersheds. Thus a daily profile of soil water was simulated with SWAT using data collected from the Stavely Range Sub‐station in the foothills of south‐western Alberta, Canada. Three small watersheds were established using a combination of natural and artificial barriers in 1996–97. The watersheds were subjected to no grazing (control), heavy grazing (2·4 animal unit months (AUM) per hectare) or very heavy grazing (4·8 AUM ha^?1). Soil‐water measurements were conducted at four slope positions within each watershed (upper, middle, lower and 5 m close to the collector drain), every 2 weeks annually from 1998 to 2000 using a downhole CPN 503 neutron moisture meter. Calibration of the model was conducted using 1998 soil‐water data and resulted in Nash–Sutcliffe coefficient (EF or R²) and regression coefficient of determination (r²) values of 0·77 and 0·85, respectively. Model graphical and statistical evaluation was conducted using the soil‐water data collected in 1999 and 2000. During the evaluation period, soil water was simulated reasonably with an overall EF of 0·70, r² of 0·72 and a root mean square error (RMSE) of 18·01. The model had a general tendency to overpredict soil water under relatively dry soil conditions, but to underpredict soil water under wet conditions. Sensitivity analysis indicated that absolute relative sensitivity indices of input parameters in soil‐water simulation were in the following order; available water capacity > bulk density > runoff curve number > fraction of field capacity (FFCB) > saturated hydraulic conductivity. Thus these data were critical inputs to ensure reasonable simulation of soil‐water patterns. Overall, the model performed satisfactorily in simulating soil‐water patterns in all three watersheds with a daily time‐step and indicates a great potential for monitoring soil‐water resources in small watersheds. Copyright © 2004 John Wiley & Sons, Ltd.     

Distributed discharge and sediment concentration predictions in the sub‐humid Ethiopian highlands: the Debre Mawi watershed

Experimental research in the Ethiopian highlands found that saturation excess induced runoff and erosion are common in the sub‐humid conditions. Because most erosion simulation models applied in the highlands are based on infiltration excess, we, as an alternative, developed the Parameter Efficient Distributed (PED) model, which can simulate water and sediment fluxes in landscapes with saturation excess runoff. The PED model has previously only been tested at the outlet of a watershed and not for distributed runoff and sediment concentration within the watershed. In this study, we compare the distributed storm runoff and sediment concentration of the PED model against collected data in the 95‐ha Debre Mawi watershed and three of its nested sub‐watersheds for the 2010 and 2011 rainy seasons. In the PED model framework, the hydrology of the watershed is divided between infiltrating and runoff zones, with erosion only taking place from two surface runoff zones. Daily storm runoff and sediment concentration values, ranging from 0.5 to over 30 mm and from 0.1 to 35 g l^?1, respectively, were well simulated. The Nash Sutcliffe efficiency values for the daily storm runoff for outlet and sub‐watersheds ranged from 0.66 to 0.82, and the Nash–Sutcliffe efficiency for daily sediment concentrations were greater than 0.78. Furthermore, the model uses realistic fractional areas for surface and subsurface flow contributions, for example between saturated areas (15%), degraded areas (30%) and permeable areas (55%) at the main outlet, while close similarity was found for the remaining hydrology and erosion parameter values. One exception occurred for the distinctly greater transport limited parameter at the actively gullying lower part of the watershed. The results suggest that the model based on saturation excess provides a good representation of the observed spatially distributed runoff and sediment concentrations within a watershed by modelling the bottom lands (as opposed to the uplands) as the dominant contributor of the runoff and sediment load. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the rainfall–runoff process of the Mara River basin using the Soil and Water Assessment Tool

The delicate balance between human utilization and sustaining its pristine biodiversity in the Mara River basin (MRB) is being threatened because of the expansion of agriculture, deforestation, human settlement, erosion and sedimentation and extreme flow events. This study assessed the applicability of the Soil and Water Assessment Tool (SWAT) model for long‐term rainfall–runoff simulation in MRB. The possibilities of combining/extending gage rainfall data with satellite rainfall estimates were investigated. Monthly satellite rainfall estimates not only overestimated but also lacked the variability of observed rainfall to substitute gage rainfall in model simulation. Uncertainties related to the quality and availability of input data were addressed. Sensitivity and uncertainty analysis was reported for alternative model components and hydrologic parameters used in SWAT. Mean sensitivity indices of SWAT parameters in MRB varied with and without observed discharge data. The manual assessment of individual parameters indicated heterogeneous response among sub‐basins of MRB. SWAT was calibrated and validated with 10 years of discharge data at Bomet (Nyangores River), Mulot (Amala River) and Mara Mines (Mara River) stations. Model performance varied from satisfactory at Mara Mines to fair at Bomet and weak at Mulot. The (Nash–Sutcliff efficiency, coefficient of determination) results of calibration and validation at Mara Mines were (0.68, 0.69) and (0.43, 0.44), respectively. Two years of moving time window and flow frequency analysis showed that SWAT performance in MRB heavily relied on quality and abundance of discharge data. Given the 5.5% area contribution of Amala sub‐basin as well as uncertainty and scarcity of input data, SWAT has the potential to simulate the rainfall runoff process in the MRB. Copyright © 2012 John Wiley & Sons, Ltd.     

SWAT application to estimate design runoff curve number for South Korean conditions

Runoff estimations based on the standard USDA–NRCS curve number (CN) table without calibration have a tendency to give inaccurate results when the CN values are applied in South Korea which has many high slope watersheds and that has a continental monsoon climate. Particularly for the design flood estimation, accurately calibrated CN values are required because the estimated peak flow is very sensitive to the selection of CN. However, the lack of flood data makes it difficult to calibrate and assign runoff CNs to Korean watersheds. Even if sufficient data are available to estimate CN values, it is also difficult to obtain the direct flows by separating base flows from total runoff hydrographs due to the temporal irregularity of rainfall events and the resulting complex pattern of runoff. Therefore, an alternative method for estimating CNs needs to be developed to overcome these issues. The purpose of this study is to present a method for estimating runoff CNs using the soil and water assessment tool (SWAT) model which can take into account watershed heterogeneities such as climate conditions, land use and soil types. The proposed CN estimation method uses the simulated flow data by SWAT instead of using measured flow data. This method has advantages in estimating CN values spatially for each subbasin division considering watershed characteristics. The use of daily data can reduce the sensitivity to the abnormality that is commonly involved in flow data with a small time scale. The SWAT‐based CN estimation method, combined with the asymptotic CN method, was applied to the Chungju dam watershed in South Korea. A regression equation was then developed from this approach, which was used to estimate CN values that decrease exponentially as rainfall amounts increase and that converge to 60·6 and 79·4 without and with considering subsurface lateral flow, respectively. Furthermore, the CN values for the antecedent moisture conditions were determined using the probabilistic approach. The CN associated with the 50% probability for the Chungju dam watershed is 87·8 which can be taken to be representative of antecedent moisture condition (AMC) II. The CN_I and CN_III associated with 90% and 10% probabilities are 78·9 and 94·1, respectively. The estimated CN_II = 87·8 differs markedly from the geographic information system (GIS)‐based CN 65·0, which implies that the standard USDA–NRCS CN method should be calibrated to the studied area of interest. Copyright © 2010 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.     

Land cover effects on runoff patterns in eastern Piedmont (USA) watersheds

Physiography and land cover determine the hydrologic response of watersheds to climatic events. However, vast differences in climate regimes and variation of landscape attributes among watersheds (including size) have prevented the establishment of general relationships between land cover and runoff patterns across broad scales. This paper addresses these difficulties by using power spectral analysis to characterize area‐normalized runoff patterns and then compare these patterns with landscape features among watersheds within the same physiographic region. We assembled long‐term precipitation and runoff data for 87 watersheds (first to seventh order) within the eastern Piedmont (USA) that contained a wide variety of land cover types, collected environmental data for each watershed, and compared the datasets using a variety of statistical measures. The effect of land cover on runoff patterns was confirmed. Urban‐dominated watersheds were flashier and had less hydrologic memory compared with forest‐dominated watersheds, whereas watersheds with high wetland coverage had greater hydrologic memory. We also detected a 10–15% urban threshold above which urban coverage became the dominant control on runoff patterns. When spectral properties of runoff were compared across stream orders, a threshold after the third order was detected at which watershed processes became dominant over precipitation regime in determining runoff patterns. Finally, we present a matrix that characterizes the hydrologic signatures of rivers based on precipitation versus landscape effects and low‐frequency versus high‐frequency events. The concepts and methods presented can be generally applied to all river systems to characterize multiscale patterns of watershed runoff. Copyright © 2013 John Wiley & Sons, Ltd.     

Climate change sensitivity assessment of streamflow and agricultural pollutant transport in California's Central Valley using Latin hypercube sampling

Bracketing the uncertainty of streamflow and agricultural runoff under climate change is critical for proper future water resource management in agricultural watersheds. This study used the Soil and Water Assessment Tool (SWAT) in conjunction with a Latin hypercube climate change sampling algorithm to construct a 95% confidence interval (95CI) around streamflow, sediment load, and nitrate load predictions under changes in climate for the Sacramento and San Joaquin River watersheds in California's Central Valley. The Latin hypercube algorithm sampled 2000 combinations of precipitation and temperature changes based on Intergovernmental Panel on Climate Change projections from multiple General Circulation Models. Average monthly percent changes of the upper and lower 95CI limits compared to the present‐day simulation and a statistic termed the “r‐factor” (average width of the 95CI band divided by the standard deviation of the 95CI bandwidth) were used to assess watershed sensitivities. 95CI results indicate that streamflow and sediment runoff in the Sacramento River watershed are more likely to decrease under climate change compared to present‐day conditions, whereas the increase and decrease for nitrate runoff were found to be equal. For the San Joaquin River watershed, streamflow slightly decreased under climate change, whereas sediment and nitrate runoff increased compared to present‐day climate. Comparisons of watershed sensitivities indicate that the San Joaquin River watershed is more sensitive to climate changes than the Sacramento River watershed, which is largely caused by the high density of agricultural land. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

A distributed monthly hydrological model for integrating spatial variations of basin topography and rainfall

Hydrological models at a monthly time‐scale are important tools for hydrological analysis, such as in impact assessment of climate change and regional water resources planning. Traditionally, monthly models adopt a conceptual, lumped‐parameter approach and cannot account for spatial variations of basin characteristics and climatic inputs. A large requirement for data often severely limits the utility of physically based, distributed‐parameter models. Based on the variable‐source‐area concept, we considered basin topography and rainfall to be two major factors whose spatial variations play a dominant role in runoff generation and developed a monthly model that is able to account for their influences in the spatial and temporal dynamics of water balance. As a hybrid of the Xinanjiang model and TOPMODEL, the new model is constructed by innovatively making use of the highly acclaimed simulation techniques in the two existing models. A major contribution of this model development study is to adopt the technique of implicit representation of soil moisture characteristics in the Xinanjiang model and use the TOPMODEL concept to integrate terrain variations into runoff simulation. Specifically, the TOPMODEL topographic index ln(a/tanβ) is converted into an index of relative difficulty in runoff generation (IRDG) and then the cumulative frequency distribution of IRDG is used to substitute the parabolic curve, which represents the spatial variation of soil storage capacity in the Xinanjiang model. Digital elevation model data play a key role in the modelling procedures on a geographical information system platform, including basin segmentation, estimation of rainfall for each sub‐basin and computation of terrain characteristics. Other monthly data for model calibration and validation are rainfall, pan evaporation and runoff. The new model has only three parameters to be estimated, i.e. watershed‐average field capacity WM, pan coefficient η and runoff generation coefficient α. Sensitivity analysis demonstrates that runoff is least sensitive to WM and, therefore, it can be determined by a prior estimation based on the climate and soil properties of the study basin. The other two parameters can be determined using optimization methods. Model testing was carried out in a number of nested sub‐basins of two watersheds (Yuanjiang River and Dongjiang River) in the humid region in central and southern China. Simulation results show that the model is capable of describing spatial and temporal variations of water balance components, including soil moisture content, evapotranspiration and runoff, over the watershed. With a minimal requirement for input data and parameterization, this terrain‐based distributed model is a valuable contribution to the ever‐advancing technology of hydrological modelling. Copyright © 2006 John Wiley & Sons, Ltd.