Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. Copyright © 2014 John Wiley & Sons, Ltd.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

A GIS‐based variable source area hydrology model

Effective control of nonpoint source pollution from contaminants transported by runoff requires information about the source areas of surface runoff. Variable source hydrology is widely recognized by hydrologists, yet few methods exist for identifying the saturated areas that generate most runoff in humid regions. The Soil Moisture Routing model is a daily water balance model that simulates the hydrology for watersheds with shallow sloping soils. The model combines elevation, soil, and land use data within the geographic information system GRASS, and predicts the spatial distribution of soil moisture, evapotranspiration, saturation‐excess overland flow (i.e., surface runoff), and interflow throughout a watershed. The model was applied to a 170 hectare watershed in the Catskills region of New York State and observed stream flow hydrographs and soil moisture measurements were compared to model predictions. Stream flow prediction during non‐winter periods generally agreed with measured flow resulting in an average r² of 0·73, a standard error of 0·01 m³/s, and an average Nash‐Sutcliffe efficiency R² of 0·62. Soil moisture predictions showed trends similar to observations with errors on the order of the standard error of measurements. The model results were most accurate for non‐winter conditions. The model is currently used for making management decisions for reducing non‐point source pollution from manure spread fields in the Catskill watersheds which supply New York City's drinking water. Copyright © 1999 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.     

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.     

Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

Estimating daily inflows of large lakes using a water‐balance‐based runoff coefficient scaling approach

Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall–runoff relationship at a reference catchment was analysed to select rainfall input and estimate run‐off coefficient firstly, and the run‐off coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Run‐off coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared with the commonly used drainage area ratio method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R² and Nash–Sutcliffe efficiency values of 0.69 and 0.6, respectively. This approach outperforms drainage area ratio method in terms of mean R² and Nash–Sutcliffe efficiency values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.     

Internal catchment process simulation in a snow‐dominated basin: performance evaluation with spatiotemporally variable runoff generation and groundwater dynamics

Hydrologic models have increasingly been used in forest hydrology to overcome the limitations of paired watershed experiments, where vegetative recovery and natural variability obscure the inferences and conclusions that can be drawn from such studies. Models are also plagued by uncertainty, however, and parameter equifinality is a common concern. Physically‐based, spatially‐distributed hydrologic models must therefore be tested with high‐quality experimental data describing a multitude of concurrent internal catchment processes under a range of hydrologic regimes. This study takes a novel approach by not only examining the ability of a pre‐calibrated model to realistically simulate watershed outlet flows over a four year period, but a multitude of spatially‐extensive, internal catchment process observations not previously evaluated, including: continuous groundwater dynamics, instantaneous stream and road network flows, and accumulation and melt period spatial snow distributions. Many hydrologic model evaluations are only on the comparison of predicted and observed discharge at a catchment outlet and remain in the ‘infant stage’ in terms of model testing. This study, on the other hand, tests the internal spatial predictions of a distributed model with a range of field observations over a wide range of hydroclimatic conditions. Nash‐Sutcliffe model efficiency was improved over prior evaluations due to continuing efforts in improving the quality of meteorological data collection. Road and stream network flows were generally well simulated for a range of hydrologic conditions, and snowpack spatial distributions were well simulated for one of two years examined. The spatial variability of groundwater dynamics was effectively simulated, except at locations where strong stream–groundwater interactions exist. Model simulations overall were quite successful in realistically simulating the spatiotemporal variability of internal catchment processes in the watershed, but the premature onset of simulated snowmelt for one of the simulation years has prompted further work in model development. Copyright © 2011 John Wiley & Sons, Ltd.     

Effective run‐off flow length over biological soil crusts on silty loam soils in drylands

Surface hydrological behaviour is important in drylands because it affects the distribution of soil moisture and vegetation and the hydrological functioning of slopes and catchments. Microplot scale run‐off can be relatively easily measured, i.e. by rainfall simulations. However, slope or catchment run‐off cannot be deduced from microplots, requiring long‐time monitoring, because run‐off coefficients decrease with increasing drainage area. Therefore, to determine the slope length covered by run‐off (run‐off length) is crucial to connect scales. Biological soil crusts (BSCs) are good model systems, and their hydrology at slope scale is insufficiently known. This study provides run‐off lengths from BSCs, by field factorial experiments using rainfall simulation, including two BSC types, three rain types, three antecedent soil moistures and four plot lengths. Data were analysed by generalized linear modelling, including vascular plant cover as covariates. Results were the following: (i) the real contributing area is almost always much smaller than the topographical contributing area; (ii) the BSC type is key to controlling run‐off; run‐off length reached 3 m on cyanobacterial crust, but hardly over 1 m on lichen crust; this pattern remained through rain type or soil moisture; (iii) run‐off decreased with BSC development because soil sealing disappears; porosity, biomass and roughness increase and some changes occur in the uppermost soil layer; and (iv) run‐off flow increased with both rain type and soil moisture but run‐off coefficient only with soil moisture (as larger rains increased both run‐off and infiltration); vascular plant cover had a slight effect on run‐off because it was low and random. Copyright © 2014 John Wiley & Sons, Ltd.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Two‐dimensional continuous simulation of spatiotemporally varied hydrological processes using the time‐area method

Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.     

On the relationship between flood and contributing area

Although it is well known that the vast majority of the time only a portion of any watershed contributes run‐off to the outlet, this extent is rarely documented. Also, the power law form of the streamflow and contributing area (Q‐A_c) relationship has been known for a half century, but it is uncommon for it to be quantified, and time series of contributing area extensive enough to calculate its frequency distribution are almost non‐existent. Data from the Canadian Prairies, where there are extensive estimates of contributing area during the median annual flood, imply that the power law coefficient for any Q‐A_c curve is a function of flow magnitude and return period. These data also suggest that regional flood frequency curves are a construct of Q‐A_c curves from individual basins. This paper will discuss research that attempted to reproduce the Q‐A_c curves for the La Salle River Watershed with a semidistributed numerical hydrological model, MESH‐PDMROF. The model simulated streamflow reasonably well (Nash Sutcliffe values = 0.62) compared with published examples of comparable models applied in the region. Estimates of the coefficient and exponent of the Q‐A_c power law function ranged from 0.08–0.14 and 0.9–1.12, respectively. These exponent values were lower than those of regional flood frequency curves and support the theory that regional flood frequency curves are a construct of Q‐A_c curves. Simulations of the area contributing to the median annual flood were lower (0.3) than those derived from independent topographic analysis (0.9) described in earlier literature though there is uncertainty in both these estimates. This uncertainty was extended across the flood frequency distribution and may be too large to definitively verify the study hypothesis.     

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (K_s) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K_s, and the spatial correlation length scale associated with K_s. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale K_s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Tracking sediment redistribution in a small watershed: implications for agro‐landscape evolution

A new, multi‐tracer method is used to track erosion, translocation, and redeposition of sediment in a small watershed, thus allowing for the ?rst time a complete, spatially distributed, sediment balance to be made as a function of landscape position. A 0·68 ha watershed near Coshocton, Ohio, USA was divided into six morphological units, each tagged with one of six rare earth element oxides. Sediment translocation was evaluated by collecting run‐off and by spatially sampling the soil surface. Average measured erosion rate was 6·1 t ha^?1, but varied between 40·4 t ha^?1 loss from the lower channels to 24·1 t ha^?1 gain on the toeslope. With this technique it was possible for the ?rst time to itemize the sediment budget for landscape elements into three components: (1) the soil from the element that left the watershed with run‐off; (2) soil from the element that was redeposited on lower positions, with the spatial distribution of that deposition; and (3) soil originating from the upper positions and deposited on the element, with quanti?cation of relative source areas. The results are incongruous with the current morphology of the watershed, suggesting that diffusion‐type erosion must also play a major role in de?ning the evolution of this landscape. Copyright © 2004 John Wiley & Sons, Ltd.     

Multi‐site calibration using a grid‐based event rainfall–runoff model: a case study of the upstream areas of the Nakdong River basin in Korea

Most runoff analyses using a grid‐based distributed model use one parameter group calibrated at the outlet of a watershed, instead of dividing the watershed into subwatersheds. Significant differences between the observed value and the simulation result of the subwatersheds can occur if just one parameter group is used in all subwatersheds that have different hydrological characteristics from each other. Therefore, to improve the simulation results of the subwatersheds within a watershed, a model calibrated at every subwatershed needs to be used to reflect the characteristics of each subwatershed. In this study, different parameter groups were set up for one or two sites using a distributed model, the GRM (Grid based Rainfall‐runoff Model), and the evaluations were based on the results of rainfall–runoff analysis, which uses a multi‐site calibration (MSC) technique to calibrate the model at the outlet of each site. The Hyangseok watershed in Naeseong River, which is a tributary of Nakdong River in Korea, was chosen as the study area. The watershed was divided into five subwatersheds each with a subwatershed outlet that was applied to the calibration sites . The MSC was applied for five cases. When a site was added for calibration in a watershed, the runoff simulation showed better results than the calibration of only one site at the most downstream area of the watershed. The MSC approach could improve the simulation results on the calibrated sites and even on the non‐calibrated sites, and the effect of MSC was improved when the calibrated site was closer to the runoff site. Copyright © 2014 John Wiley & Sons, Ltd.     

The role of land surface versus drainage network characteristics in controlling water quality and quantity in a small urban watershed

The processes that control run‐off quantity and quality in urban watersheds are complex and not well understood. Although impervious surface coverage has traditionally been used to examine altered hydrologic response in urban watersheds, several studies suggest that other elements of the urban landscape, particularly those associated with urban infrastructure and the drainage system, play an equally important role. The relative importance of impervious surfaces, stormwater ponds, expansion of the drainage network, and drainage network structures in controlling hydrologic response was examined in the subwatersheds of the Kromma Kill, an urban watershed located in Albany County, NY. In this study, geographic information systems was used to compute geospatial land surface and drainage network properties of 5 Kromma Kill subwatersheds. In these same subwatersheds, water quantity (rainfall and run‐off) and quality (macroinvertebrates, nitrate, total nitrogen, dissolved oxygen, total dissolved solids, and nonpurgable organic carbon) parameters were measured. Strong and significant correlations were identified between land surface and drainage network properties and field observations. Causal relationships were then tested using the Environmental Protection Agency's Stormwater Management Model. Field and model analyses suggest that whereas percent imperviousness is a dominant control on water quality, drainage density and slope are equally important. However, for water quantity, whereas imperviousness is positively correlated with increased run‐off volumes, drainage network properties and slope are the dominant controls on run‐off volumes. Results have important implications for stormwater management plans, especially those aimed at reducing the effective impervious surface coverage of urban watersheds. Reducing the percentage of effective imperviousness in a watershed is not a “one size fits all” solution and can help to meet some management objectives, such as reducing nitrogen concentrations and improving water quality, but may not serve as the most effective, and therefore economical, solution for every management objective including reducing run‐off volumes.     

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

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

Surface moisture area during rainfall–run‐off events to understand the hydrological dynamics of a basin in a plain region

The understanding of the hydrology of plain basins may be improved by the combined analysis of rainfall–run‐off records and remote sensed surface moisture data. Our work evaluates the surface moisture area (SMA) produced during rainfall–run‐off events in a plain watershed of the Argentine Pampas Region, and studies which hydrological variables are related to the generated SMA. The study area is located in the upper and middle basins of the Del Azul stream, characterized by the presence of small gently hilly areas surrounded by flat landscapes. Data from 9 rainfall–run‐off events were analysed. MODIS surface reflectance data were processed to calculate SMA subsequent to the peak discharge (post‐SMA), and previous to the rainfall events (prev‐SMA), to consider the antecedent wetness. Rainfall–run‐off data included total precipitation depth (P), maximum intensity of rainfall over 6 hr (I6max), surface run‐off registered between the beginning of the event and the day previous to the analysed MODIS scene (R), peak flow (Qp), and flood intensity (IF). In contrast with other works, post‐SMA showed a negative relationship with the R. Three groups of cases were identified: (a) Events of low I6max, high prev‐SMA, and low R were associated with slow and weakly channelized flow over plain areas, leading to saturated overland flow (SOF), with large SMA; (b) events of high I6max, low prev‐SMA, and medium to high R were rapidly transported along the gentle slopes of the basin, related to Hortonian overland flow (HOF) and low post‐SMA; and (c) events of medium to high I6max and prev‐SMA with medium R were related to heterogeneous input‐antecedent‐run‐off conditions combined: Local spatial conditions may have produced HOF or SOF, leading to an averaged response with medium SMA. The interactions between the geomorphology of the basin, the characteristics of the events, and the antecedent conditions may explain the obtained results. This analysis is relevant for the general knowledge of the hydrology of large plains, whose functioning studies are still in their early stages.     

Investigating hydrologic connectivity and its association with threshold change in runoff response in a temperate forested watershed

Hydrological studies across varied climatic and physiographic regions have observed small changes in the ‘states of wetness’; based on average soil moisture, can lead to dramatic changes in the amount of water delivered to the stream channel. This non-linear behaviour of the storm response has been attributed to a critical switching in spatial organization of shallow soil moisture and hydrologic connectivity. However, much of the analysis of the role of soil moisture organization and connectivity has been performed in small rangeland catchments. Therefore, we examined the relationship between hydrologic connectivity and runoff response within a temperate forested watershed of moderate relief. We have undertaken spatial surveys of shallow soil moisture over a sequence of storms with varying antecedent moisture conditions. We analyse each survey for evidence of hydrologic connectivity and we monitor the storm response from the catchment outlet. Our results show evidence of a non-linear response in runoff generation over small changes in measures of antecedent moisture conditions; yet, unlike the previous studies of rangeland catchments, in this forested landscape we do not observe a significant change in geostatistical hydrologic connectivity with variations in antecedent moisture conditions. These results suggest that a priori spatial patterns in shallow soil moisture in forested terrains may not always be a good predictor of critical hydrologic connectivity that leads to threshold change in runoff generation, as has been the case in rangeland catchments. Copyright © 2007 John Wiley & Sons, Ltd.