Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.     

Comparison of stormflow responses of surface‐mined and forested watersheds in the Appalachian Mountains,USA

The results of a hydrological analysis that was conducted as part of a larger, multifaceted, collaborative effort to quantify ecosystem functions in watersheds subjected to land‐use and land‐cover change are presented. The primary goal of the study was to determine whether a small watershed in the Appalachian region (USA) that was recently subjected to surface mining and reclamation practices produces stormflow responses to rain events that are different from those produced by a nearby reference watershed covered by young, second‐growth forest. Water balances indicated that runoff yields did not vary significantly between the two watersheds on an annual basis. Statistically significant differences (p?0·05) in runoff responses were observed on an event basis, however, with the mined/reclaimed watershed producing, on average (a) higher storm runoff coefficients (2·5×), (b) greater total storm runoff (3×), and (c) higher peak hourly runoff rates (2×) when compared with the reference watershed. Results of a unit hydrograph analysis also showed, unexpectedly, that the modelled unit responses of the two watersheds to effective rainfall pulses were similar, despite the noted differences in land cover. Differences in stormflow responses were thus largely explained by dramatic reductions in cumulative rates of rainfall abstraction (measured using infiltrometers) attributable to soil compaction during land reclamation. Additional field hydrological measurements on other mined watersheds will be needed to generalize our results, as well as to understand and predict the cumulative hydrological impacts of widespread surface mining in larger watersheds and river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 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.     

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.     

Plot and small‐watershed scale runoff from two reclaimed surface‐mined watersheds in Alberta

Infiltration excess overland flow has been identified as the dominant flow pathway in recently reclaimed surface mined watersheds as a result of compaction and sorting during the reclamation procedure. Therefore, there could be a fairly direct relationship between runoff generated from the hillslopes to that measured at the watershed outlet. A 3‐year study was initiated in 1993 to determine how well surface runoff at a watershed scale could be predicted from 1‐m² runoff frames placed on hillslopes in two reclaimed surface‐mined watersheds in central Alberta. Runoff from the hillslope frames suggests outlet discharge should be high from the 3\4‐ha Sandy Subsoil Watershed and much less for the 9\8‐ha West Watershed, but the opposite occurred. Most of the hillslope runoff from the Sandy Subsoil Watershed infiltrated once it reached the channel and depression storage played an insignificant role in determining runoff. In contrast, most of the runoff from the West Watershed originated from rain falling directly on the saturated channel (depression storage) or near‐channel saturated areas, rather than the hillslopes. Neither watershed runoff magnitude nor timing could be predicted from the same parameters for hillslope runoff frames for either reclaimed watershed. Copyright © 2000 John Wiley & Sons, Ltd.     

Understanding plot‐scale hydrology of Lesser Himalayan watershed—A field study and HYDRUS‐2D modelling approach

Soil moisture dynamics have a significant effect on overland flow generation. Catchment aspect is one of the major controlling factors of overland flow and soil moisture behaviour. A few experimental studies have been carried out in the uneven topography of the Himalayas. This study presents plot‐scale experiments using portable rainfall simulator at an altitude of 1,230 m above mean sea level and modelling of overland flow using observed datasets. Two plots were selected in 2 different aspects of Aglar watershed of Lesser Himalaya; the agro‐forested (AF) plot was positioned at the north aspect whereas the degraded (DE) plot was located at the south aspect of the hillslope. HS flumes and rain gauges were installed to measure the runoff at the outlet of the plot and the rainfall depth during rainfall simulation experiments. Moreover, 10 soil moisture sensors were installed at upslope and downslope locations of both the plots at 5, 15, 25, 35, and 45 cm depth from ground level to capture the soil moisture dynamics. The tests were conducted at intensities of 79.8 and 75 mm/hr in AF plot and 82.2 and 72 mm/hr in the DE plot during Test 1 and Test 2, respectively. The observed data indicate the presence of reinfiltration process only in the AF plot. The high water holding capacity and the presence of reinfiltration process results in less runoff volume in the AF plot compared with the DE plot. The Hortonian overland flow mechanism was found to be the dominant overland flow mechanism as only a few layers of top soil get saturated during all of the rainfall–runoff experiments. The runoff, rainfall, and soil moisture data were subsequently used to calibrate the parameters of HYDRUS‐2D overland flow module to simulate the runoff hydrograph and soil moisture. The components of hydrograph were evaluated in terms of peak discharge, runoff volume and time of concentration, the results were found to be within the satisfactory range. The goodness of fit of simulated hydrographs were more than 0.85 and 0.95 for AF and DE plot, respectively. The model produced satisfactory simulation results of soil moisture for all of the rainfall–runoff experiments. The HYDRUS‐2D overland flow module was found promising to simulate the runoff hydrograph and soil moisture in plot‐scale research.     

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.     

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.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. 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.     

Impact of spatial variations of land surface parameters on regional evaporation: a case study with remote sensing data

Most precipitation in watersheds is consumed by evaporation, thus techniques to appraise regional evaporation are important to assess the availability of water resources. Many algorithms to estimate evaporation from remotely sensed spectral data have been developed in the recent past. In addition to differences in the physical parameterization of surface fluxes, these algorithms have different solutions for describing spatial variations of the parameters in the soil–vegetation–atmosphere–transfer (SVAT) continuum. In this study, the necessity to spatially distinguish SVAT parameters for computing surface heat fluxes is analysed for the Naivasha watershed in the Kenyan Rift Valley. Landsat Thematic Mapper (TM) spectral data have been used to first delineate the watershed into 15 hydrological units using surface temperature, normalized difference vegetation index and surface albedo as attributes. Thereafter, semi‐empirical relationships between these TM‐based parameters and other SVAT parameters have been applied to compute the spatial variation of SVAT parameters and the associated evaporation from the different hydrological units. The impact of using watershed‐constant or watershed‐distributed SVAT parameters on the fluxes is analysed. The determination of watershed averaged evaporation with area‐aggregated SVAT parameters is feasible without significant loss of accuracy. Distributed evaporation in heterogeneous watersheds, however, can be investigated only with remote sensing flux algorithms that can account for spatially variable air temperature, surface roughness, surface albedo and the stability correction of the temperature profile due to buoyancy. Erroneous results can be expected if area‐aggregated SVAT parameters are used to calculate local evaporation. As most of the recently developed remote sensing flux algorithms are based on areal constant SVAT parameters, direct applications in watersheds are still limited. Copyright © 2001 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.     

Macropores and preferential flow—a love‐hate relationship

Preferential flow is of high relevance for runoff generation, transport of chemicals and nutrients, and the transit time distribution of water in the soil or watershed. However, preferential flow effects are generally ignored in lumped hydrological models. And even most physically‐based models ignore macropores and preferential flow features at the soil and hillslope scale. Keith Beven was never satisfied with this situation and he tried again and again to convince the scientific community to focus their research on the complex topic of macropore and preferential flow. Although he recognized how difficult it is to correctly include preferential flow in hydrological models, he made substantial progress defining and describing macropore flow and showing its relevance, developing models to simulate preferential flow, and in particular, the interaction between macropores and the soil matrix. In this short commentary, I reflect on these achievements and outline a vision for research in preferential flow experiments and modeling.     

A kinematic‐wave‐based distributed watershed model using FEM,GIS and remotely sensed data

For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of infiltration‐based stormwater management to restore hydrological processes in urban headwater streams

Urbanization threatens headwater stream ecosystems globally. Watershed restoration practices, such as infiltration‐based stormwater management, are implemented to mitigate the detrimental effects of urbanization on aquatic ecosystems. However, their effectiveness for restoring hydrologic processes and watershed storage remains poorly understood. Our study used a comparative hydrology approach to quantify the effects of urban watershed restoration on watershed hydrologic function in headwater streams within the Coastal Plain of Maryland, USA. We selected 11 headwater streams that spanned an urbanization–restoration gradient (4 forested, 4 urban‐degraded, and 3 urban‐degraded) to evaluate changes in watershed hydrologic function from both urbanization and watershed restoration. Discrete discharge and continuous, high‐frequency rainfall‐stage monitoring were conducted in each watershed. These datasets were used to develop 6 hydrologic metrics describing changes in watershed storage, flowpath connectivity, or the resultant stream flow regime. The hydrological effects of urbanization were clearly observed in all metrics, but only 1 of the 3 restored watersheds exhibited partially restored hydrologic function. At this site, a larger minimum runoff threshold was observed relative to the urban‐degraded watersheds, suggesting enhanced infiltration of stormwater runoff within the restoration structure. However, baseflow in the stream draining this watershed remained low compared to the forested reference streams, suggesting that enhanced infiltration of stormwater runoff did not recharge subsurface storage zones contributing to stream baseflow. The highly variable responses among the 3 restored watersheds were likely due to the spatial heterogeneity of urban development, including the level of impervious cover and extent of the storm sewer network. This study yielded important knowledge on how restoration strategies, such as infiltration‐based stormwater management, modulated—or failed to modulate—hydrological processes affected by urbanization, which will help improve the design of future urban watershed management strategies. More broadly, we highlighted a multimetric approach that can be used to monitor the restoration of headwater stream ecosystems in disturbed landscapes.     

Two‐dimensional physically based finite element runoff model for small agricultural watersheds: II. Model testing and field application

Mathematical models are useful analysis tools to understand problems in watersheds associated with runoff, and to find solutions through land use changes and best management practices. However, before a model is applied in the field, it must be tested and checked to ensure that the model represents the real world adequately. In this paper, a two‐dimensional physically based finite element runoff model ROMO2D has been verified and validated by comparing the model output with analytic solution under simplified conditions, published data, and field measurements. Calibration of the model was done manually through a multi‐objective calibration procedure, using observed field data. Before going for field validation/application of ROMO2D, analysis was carried out to determine the optimum number of finite elements into which the watershed should be discretized and the size of the time step. A sensitivity analysis of the model was performed using the observed values of watershed parameters. The model was applied to a 1·45 ha agricultural watershed located in the Shiwalik foothills (India) to simulate runoff. The results demonstrated the potential of the model to simulate runoff from small agricultural watersheds for individual storm events with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

不同尺度流域地表径流氮、磷浓度比较

选择太湖上游为研究对象,采集了1-400 km2不同尺度小流域产出径流TN、TP浓度实测数据,结合前期开展的地表坡面流人工暴雨实验监测结果,开展不同尺度流域水质监测对水体面源污染产出浓度估算影响的比较研究,探讨流域尺度之间入渗、汇流以及伴随的流域生态系统营养盐调节机制的差异．结果表明,流域监测尺度对土地利用面源污染产出浓度估算有较大影响．地表坡面流由于未经过流域汇流过程伴随的下渗滤过与吸附等过程,产出径流TN、TP浓度一般高于小流域．小流域林地生态系统具有较强的入渗机制、接近自然的生态沟谷汇流网络,对面源污染TN、TP有较强的削减作用．农业生态系统较弱的入渗机制、人工沟渠汇流网络对面源污染TN、TP的削减作用较弱．现代农业造成流域面源污染增加不仅仅是因为人类农业活动对流域局部土体及养分的改变,农业生态系统改变流域自然生态系统整体水文过程及营养盐调节机制也是面源污染增加的重要因素之一,恢复小尺度的生态沟谷网络系统对削减流域面源污染具有重要的意义．     

Effect of watershed subdivision on simulation of water balance components

The present effect of watershed subdivision on simulated water balance components using the thoroughly tested Soil and Water Assessment Tool (SWAT) model has been evaluated for the Nagwan watershed in eastern India. Observed meteorological and hydrological data (daily rainfall, temperature, relative humidity and runoff) for the years 1995 to 1998 were collected and used. The watershed and sub‐watershed boundaries, slope and soil texture maps were generated using a geographical information system. A supervised classification method was used for land‐use/cover classification from satellite imagery of 1996. In order to study the effect of watershed subdivision, the watershed was spatially defined into three decomposition schemes, namely a single watershed, and 12 and 22 sub‐watersheds. The simulation using the SWAT model was done for a period of 4 years (1995 to 1998). Results of the study showed a perfect water balance for the Nagwan watershed under all of the decomposition schemes. Results also revealed that the number and size of sub‐watersheds do not appreciably affect surface runoff. Except for runoff, there was a marked variation in the individual components of the water balance under the three decomposition schemes. Though the runoff component of the water balance showed negligible variation among the three cases, variations were noticed in the other components: evapotranspiration (5 to 48%), percolation (2 to 26%) and soil water content (0·30 to 22%). Thus, based on this study, it is concluded that watershed subdivision has a significant effect on the water balance components. Copyright © 2005 John Wiley & Sons, Ltd.     

Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.