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.     

Water flux and sediment transport within a forested landscape: the role of connectivity,subsurface flow,and slope length scale on transport mechanism

The connectivity and upscaling of overland runoff and sediment transport are important issues in hillslope hydrology to identify water flux and sediment transport within landscape. These processes are highly variable in time and space with regard to their interactions with vegetation and soil surface conditions. The generation of overland runoff and its spatial connectivity were examined along a slope to determine the variations in the transport mechanism of runoff and soil particles by rain splash and overland runoff. Field experiments were conducted by erosion plots on a steep hillslope at lengths of 5, 10, and 15 m. The overland runoff connectivity and flow transport distance decreased with the slope length, while spatial variability of infiltration increased significantly with the slope length. Observation of subsurface flow revealed that surface soil and litter layer could have important role in water transport. However, the surface soil water content and water flux transport along the slope was highly variable for different storm events; the variability was related to the complexity of the system, mainly by way of the initial wetness conditions and infiltration characteristics. Only net rain‐splashed soil was measurable, but examination of the water flux, overland runoff and sediment transport connectivity, characteristics of sheetwash, and the variability in spatial infiltration indicated an increase in the contribution of the rain splash transport mechanism along the slope. Copyright © 2013 John Wiley & Sons, Ltd.     

Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Physically based modelling of sediment generation and transport under a large rainfall simulator

A series of large rainfall simulator experiments was conducted in 2002 and 2003 on a small plot located in an experimental catchment in the North Island of New Zealand. These experiments measured both runoff and sediment transport under carefully controlled conditions. A physically based hydrological modelling system (SHETRAN) was then applied to reproduce the observed hydrographs and sedigraphs. SHETRAN uses physically based equations to represent flow and sediment transport, and two erodibility coefficients to model detachment of soil particles by raindrop erosion and overland flow erosion. The rate of raindrop erosion also depended on the amount of bare ground under the simulator; this was estimated before each experiment. These erodibility coefficients were calibrated systematically for summer and winter experiments separately, and lower values were obtained for the summer experiments. Earlier studies using small rainfall simulators in the vicinity of the plot also found the soil to be less erodible in summer and autumn. Limited validation of model parameters was carried out using results from a series of autumn experiments. The modelled suspended sediment load was also sensitive to parameters controlling the generation of runoff from the rainfall simulator plot; therefore, we found that accurate runoff predictions were important for the sediment predictions, especially from the experiments where the pasture cover was good and overland flow erosion was the dominant mechanism. The rainfall simulator experiments showed that the mass of suspended sediment increased post‐grazing, and according to the model this was due to raindrop detachment. The results indicated that grazing cattle or sheep on steeply sloping hill‐country paddocks should be carefully managed, especially in winter, to limit the transport of suspended sediment into watercourses. Copyright © 2006 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 1. model development and event‐scale dynamics

This paper presents an erosion model, ARMOUR, which simulates time‐varying runoff, erosion, deposition and surface armour evolution down a hillslope either as a result of a single erosion event or as the cumulative impact of many events over periods up to decades. ARMOUR simulates sediment transport for both cohesive and non‐cohesive soil and dynamically differentiates between ‘transport‐limited’ and ‘source‐limited’ processes. A variety of feasible processes for entrainment of different size classes can be modelled and evaluated against data. The generalized likelihood of uncertainty estimation (GLUE) technique was used to calibrate and validate ARMOUR using data collected during rainfall simulator experiments at two contrasting sites: (1) non‐cohesive stony sediments at Ranger Uranium Mine, Northern Territory, Australia; and (2) cohesive silty sediments at Northparkes Gold Mine, NSW, Australia. The spatial and temporal variations of model predictions within the individual runoff events showed that some entrainment processes could not model the spikes in concentration and subsequent depletion, while the hiding model of Andrews and Parker best simulated the concentration trends for both calibrated and independent runoff events. ARMOUR also successfully captured the coarsening of the surface material, though small, over the duration of the rainfall simulator trials. This was driven by the depletion of the finest size class of the soil. For a constant discharge, ARMOUR simulated higher sediment flux at the start of the storm with the sediment flux and concentration diminishing with time. For natural rainfall a power law relationship between sediment flux and discharge was observed. The calibration exercise showed that sediment concentration and discharge alone are insufficient to calibrate all aspects of the physics, in particular the armour depth. This appears to be because the armouring during the short duration events is driven by depletion of the finest classes of the sediments (diameters less then 62·5 mm), which are not normally measured. Copyright © 2006 John Wiley & Sons, Ltd.     

Hortonian overland flow from Japanese forest plantations—an aberration,the real thing,or something in between?

There is a growing opinion that poorly managed plantation forests in Japan are contributing to increased storm runoff and erosion. Here we present evidence to the contrary from runoff plots at two scales (hillslope and 0·5 × 2 m plots) for several forest conditions in the Mie and Nariki catchments. Runoff coefficients from small plots in untended hinoki forests were variable but typically higher than from better managed or deciduous forests during small storms at Nariki; at Mie, runoff during small events was highly variable from all small plots but runoff coefficients were similar for hinoki plots with and without understory vegetation, while the deciduous plot had lower runoff coefficients. Storm runoff was less at the hillslope scale than the plot scale in Mie; these results were more evident at sites with better ground cover. During the largest storms at both sites, differences in runoff due to forest condition were not evident regardless of scale. Dynamic soil moisture tension measurements at Nariki indicated that during a large storm, flow in the upper organic‐rich and root‐permeated soil horizons was 3·2 times higher than measured overland runoff from a small hinoki plot with poor ground cover and 8·3 times higher than runoff from a deciduous forest plot. On the basis of field observations during storms, at least a portion of the monitored ‘Hortonian overland flow’ was actually occurring in this near‐surface ‘biomat’. Therefore our field measurements in both small and large plots potentially included biomat flow in addition to short‐lived Hortonian runoff. Because overland flow decreased with increasing scale, rill erosion did not occur on hillslopes. Additionally, runoff coefficients were not significantly different among cover conditions during large storms; thus, the ‘degraded’ forest conditions appear not to greatly enhance peak flows or erosion potential at larger scales, especially when biomat flow is significant. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling water erosion in the Sahel: application of a physically based soil erosion model in a gentle sloping environment

Water is a major limiting factor in arid and semi‐arid agriculture. In the Sahelian zone of Africa, it is not always the limited amount of annual rainfall that constrains crop production, but rather the proportion of rainfall that enters the root zone and becomes plant‐available soil moisture. Maximizing the rain‐use efficiency and therefore limiting overland flow is an important issue for farmers. The objectives of this research were to model the processes of infiltration, runoff and subsequent erosion in a Sahelian environment and to study the spatial distribution of overland flow and soil erosion. The wide variety of existing water erosion models are not developed for the Sahel and so do not include the unique Sahelian processes. The topography of the Sahelian agricultural lands in northern Burkina Faso is such that field slopes are generally low (0–5°) and overland flow mostly occurs in the form of sheet flow, which may transport large amounts of fine, nutrient‐rich particles despite its low sediment transport capacity. Furthermore, pool formation in a field limits overland flow and causes resettlement of sediment resulting in the development of a surface crust. The EUROSEM model was rewritten in the dynamic modelling code of PCRaster and extended to account for the pool formation and crust development. The modelling results were calibrated with field data from the 2001 rainy season in the Katacheri catchment in northern Burkina Faso. It is concluded that the modified version of EUROSEM for the Sahel is a fully dynamic erosion model, able to simulate infiltration, runoff routing, pool formation, sediment transport, and erosion and deposition by inter‐rill processes over the land surface in individual storms at the scale of both runoff plots and fields. A good agreement is obtained between simulated and measured amounts of runoff and sediment discharge. Incorporating crust development during the event may enhance model performance, since the process has a large influence on infiltration capacity and sediment detachment in the Sahel. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Interstorm surface preparation and sediment detachment by vehicle traffic on unpaved mountain roads

Road survey and field rainfall simulation experiments have shown that the erodibility of a road surface is dynamic. In the absence of extreme runoff events, dynamic erodibility results from the generation and removal of easily entrained surface material by human road surface maintenance activities, vehicular detachment and overland flow events. Maintenance activities introduce easily transportable material to the road surface where it can be entrained by overland flow. Traffic in dry conditions detaches material that is quickly removed during subsequent overland flow events. The pre‐storm erodibility of a road is therefore largely a function of maintenance and vehicle activity since the last overland flow event. During rainstorms, vehicle passes increase sediment production by detaching/redistributing surface material and creating efficient overland flow pathways for sediment transport. However, if incision of tracks by overland flow does not occur, post‐pass sediment transport quickly returns to pre‐pass rates. Field rainfall simulation data suggest that sediment transport resulting from during‐storm vehicle passes is greatly influenced by the presence of existing loose material, which again is a function of prior road usage and maintenance activities. Incorporation of vehicular effects into physically based road erosion models may be possible by parameterizing both during‐storm and inter‐storm changes in the supply of loose surface material as changes in surface erodibility. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of climate change on overland flow generation: a case study in central Germany

The impact of global climate change on runoff components, especially on the type of overland flow, is of utmost significance. High‐resolution temporal rainfall plays an important role in determining the hydrological response of quick runoff components. However, hydrological climate change scenario analyses with high temporal resolution are rare. This study investigates the impact of climate change on discharge peak events generated by rainfall, snowmelt, and soil‐frost induced runoff using high‐resolution hydrological modelling. The study area is Schäfertal catchment (1.44 km²) in the lower Harz Mountains in central Germany. The WaSiM‐ETH hydrological model is used to investigate the rainfall response of runoff components under near future (2021–2050) and far‐distant future (2071–2100) climatic conditions. Disaggregated daily climate variables of WETTREG2010 SRES scenario A1B are used on a temporal resolution of 10 min. Hydrological model parameter optimization and uncertainty analysis was conducted using the Differential Evolution Adaptive Metropolis (DREAM_(ZS)) uncertainty tool. The scenario results show that total runoff and interflow will increase by 3.8% and 3.5% in the near future and decrease by 32.85% and 31% in the far‐distant future compared to the baseline scenario. In contrast, overland flow and the number and size of peak runoff will decrease moderately for the near future and drastically for the far‐distant future compared to the baseline scenario. We found the strongest decrease for soil‐frost induced discharge peaks at 79.6% in the near future and at 98.2% in the far‐distant future scenario. It can be concluded that high‐resolution hydrological modelling can provide detailed predictions of future hydrological regimes and discharge peak events of the catchment. Copyright © 2014 John Wiley & Sons, Ltd.     

Uncertainties in DRAINMOD predictions of subsurface drain flow for an Indiana silt loam using the GLUE methodology

Good modelling practice requires the incorporation of uncertainty analysis into hydrologic/water quality models. The generalized likelihood uncertainty estimation procedure was used to evaluate the uncertainty in DRAINMOD predictions of daily, monthly, and yearly subsurface drain flow. A variance‐based sensitivity analysis technique, the extended Fourier amplitude sensitivity test, was used to identify the main sources of prediction uncertainty. The analysis was conducted for the experimental drainage field at the Southeast Purdue Agricultural Center in Indiana. Six years of data were used and the uncertainties in eight model parameters were considered to analyse how uncertainties in input parameters propagate to model outputs. The width of 90% confidence interval bands of drain flow ranged from 0 to 0·6 cm day^?1 for daily predictions, from 0 to 3·1 cm month^?1 for the monthly predictions, and from 7·6 to 12·4 cm year^?1 for yearly predictions. Annual drain flow predicted by DRAINMOD fell well within the 90% confidence bounds. Model results were most sensitive to the vertical saturated hydraulic conductivity of the restrictive layer and the lateral hydraulic conductivity of the deepest soil layer, followed by the lateral hydraulic conductivity of the top soil layer and surface micro‐storage. Parameter interactions also contributed to the prediction uncertainty. Copyright © 2006 John Wiley & Sons, Ltd.     

Prescribed‐fire effects on rill and interrill runoff and erosion in a mountainous sagebrush landscape

Changing fire regimes and prescribed‐fire use in invasive species management on rangelands require improved understanding of fire effects on runoff and erosion from steeply sloping sagebrush‐steppe. Small (0·5 m²) and large (32·5 m²) plot rainfall simulations (85 mm h^–1, 1 h) and concentrated flow methodologies were employed immediately following burning and 1 and 2 years post‐fire to investigate infiltration, runoff and erosion from interrill (rainsplash, sheetwash) and rill (concentrated flow) processes on unburned and burned areas of a steeply sloped sagebrush site on coarse‐textured soils. Soil water repellency and vegetation were assessed to infer relationships in soil and vegetation factors that influence runoff and erosion. Runoff and erosion from rainfall simulations and concentrated flow experiments increased immediately following burning. Runoff returned to near pre‐burn levels and sediment yield was greatly reduced with ground cover recovery to 40 per cent 1 year post‐fire. Erosion remained above pre‐burn levels on large rainfall simulation and concentrated flow plots until ground cover reached 60 per cent two growing seasons post‐fire. The greatest impact of the fire was the threefold reduction of ground cover. Removal of vegetation and ground cover and the influence of pre‐existing strong soil‐water repellency increased the spatial continuity of overland flow, reduced runoff and sediment filtering effects of vegetation and ground cover, and facilitated increased velocity and transport capacity of overland flow. Small plot rainfall simulations suggest ground cover recovery to 40 per cent probably protected the site from low‐return‐interval storms, large plot rainfall and concentrated flow experiments indicate the site remained susceptible to elevated erosion rates during high‐intensity or long duration events until ground cover levels reached 60 per cent. The data demonstrate that the persistence of fire effects on steeply‐sloped, sandy sagebrush sites depends on the time period required for ground cover to recover to near 60 per cent and on the strength and persistence of ‘background’ or fire‐induced soil water repellency. Published in 2009 by John Wiley & Sons, Ltd.     

Scale effects of Hortonian overland flow and rainfall‐runoff dynamics: laboratory validation of a process‐based model

Hortonian runoff was measured in the laboratory from uniform slopes of lengths of 1·5, 3·0, and 6·0 m for steady, high‐intensity rainstorms with durations of 1·0 to 7·5 min. A clear reduction in runoff per unit slope length was found as slope lengths were increased. This effect becomes more pronounced with decreasing storm duration. The runoff data were used to validate a simple process‐based model that combines the Philip‐two‐term infiltration equation with the kinematic wave overland flow principle. The predicted and experimental results agreed well. Laboratory findings were extrapolated with the aid of the model to slopes and rainfall durations similar to those found under West African conditions. The calculated reduction of runoff per unit length is similar to reported observations. Thus, this process‐based model can largely explain the phenomenon of runoff reduction with increasing slope length. Copyright © 2002 John Wiley & Sons, Ltd.     

Experimental study and mathematical modelling of soluble chemical transfer from unsaturated/saturated soil to surface runoff

A one‐dimensional, two‐layer solute transport model is developed to simulate chemical transport process in an initially unsaturated soil with ponding water on the soil surface before surface runoff starts. The developed mathematical model is tested against a laboratory experiment. The infiltration and diffusion processes are mathematically lumped together and described by incomplete mixing parameters. Based on mass conservation and water balance equations, the model is developed to describe solute transport in a two‐zone layer, a ponding runoff zone and a soil mixing zone. The two‐zone layer is treated as one system to avoid describing the complicated chemical transport processes near the soil surface in the mixing zone. The proposed model was analytically solved, and the solutions agreed well with the experimental data. The developed experimental method and mathematical model were used to study the effect of the soil initial moisture saturation on chemical concentration in surface runoff. The study results indicated that, when the soil was initially saturated, chemical concentration in surface runoff was significantly (two orders of magnitude) higher than that with initially unsaturated soil, while the initial chemical concentrations at the two cases were of the same magnitude. The soil mixing depth for the initially unsaturated soil was much larger than that for the initially saturated soil, and the incomplete runoff mixing parameter was larger for the initially unsaturated soil. The higher the infiltration rate of the soil, the greater the infiltration‐related incomplete mixing parameter. According to the quantitative analysis, the soil mixing depth was found to be sensitive for both initially unsaturated and saturated soils, and the incomplete runoff mixing parameter was sensitive for initially saturated soil but not for the initially unsaturated soil; the incomplete infiltration mixing parameter behaved just the opposite. Some suggestions are made for reducing chemical loss from runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling catchment‐scale erosion patterns in the East African Highlands

Prompt location of areas exposed to high erosion is of the utmost importance for soil and water conservation planning. Erosion models can be useful tools to locate sources of sediment and areas of deposition within a catchment, but the reliability of model predictions of spatial patterns of erosion at catchment scale has seldom been validated against observations. This study aimed to evaluate the performance of a simple empirical model (Morgan, Morgan and Finney model, MMF) in predicting spatial patterns of erosion at two small catchments in the East African Highlands: Kwalei (Tanzania) and Gikuuri (Kenya). Erosion maps predicted by the MMF model were compared with erosion maps obtained by direct survey. In Kwalei, erosion features were especially frequent in fields of annual crops. In Gikuuri, slope was the critical erosion factor, with estimated erosion rates >10 kg m^?2 a^?1 on slopes >18 per cent. Predicted erosion rates were mainly transport‐limited and ranged from <0·01 to 13·50 kg m^?2 a^?1 in Kwalei and 9·29 kg m^?2 a^?1 in Gikuuri. The performance of the MMF model in predicting the spatial patterns of erosion was acceptable in Kwalei, but poor in Gikuuri. However, by excluding the elements at the valley bottoms in Gikuuri Catchment, the performance of the model improved dramatically. The spatial pattern of erosion predicted by the MMF model was driven by the accumulation of surface runoff, which did not consider the possibility of re‐infiltration along the slope. As a result, the MMF erosion patterns predicted by the model increased invariably from the ridges to the valley bottoms, hampering the model suitability for locating areas subjected to high and very high erosion. It is concluded that the model predictions could be substantially improved by introducing a more realistic hydrological component for the prediction of surface runoff along the hill‐slope. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigations on the runoff generation at the profile and plot scales,Swiss Emmental

This article describes an investigation on runoff generation at different scales in the forested catchment of the Sperbelgraben in the Emmental region (Swiss Prealps) where studies in the field of forest hydrology have a history of 100 years. It focuses on the analysis of soil profiles and the subsequent sprinkling experiments above them (1 m²), as well as on surface runoff measurements on larger plots (50 to 110 m²). In the Sperbelgraben investigation area, two very distinct runoff reactions could be observed. On the one hand, very high production of saturation overland flow was registered on wet areas of gleyic soils, with runoff coefficients between 0·39 and 0·94 for profile irrigation. On the other hand, almost no surface runoff was measured on Cambisols, with the exception at some sites of a hydrophobic reaction detected at the beginning of storms after dry periods (coefficients for profile irrigation: 0·01–0·16). This pattern was observed during 1 m² soil plot irrigation and on surface runoff plots. Apart from a less distinctive signal of the water‐repellent litter layer on the larger surface runoff plots, the dominant hydrological processes at the two scales are the same. The determined runoff reaction at the two scales corresponds well with information from a forest site type map describing soil and vegetation characteristics and used as a substitute for a soil map in this study. Theoretical considerations describing forest influence on flood discharge are discussed and evaluated to be in good agreement with observations. These findings are a sound foundation for application in hydrological catchment modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Erosion prediction on unpaved mountain roads in northern Thailand: validation of dynamic erodibility modelling using KINEROS2

The event‐ and physics‐based KINEROS2 runoff/erosion model for predicting overland flow generation and sediment production was applied to unpaved mountain roads. Field rainfall simulations conducted in northern Thailand provided independent data for model calibration and validation. Validation shows that KINEROS2 can be parameterized to simulate total discharge, sediment transport and sediment concentration on small‐scale road plots, for a range of slopes, during simulated rainfall events. The KINEROS2 model, however, did not accurately predict time‐dependent changes in sediment output and concentration. In particular, early flush peaks and the temporal decay in sediment output were not predicted, owing to the inability of KINEROS2 to model removal of a surface sediment layer of finite depth. After 15–20 min, sediment transport declines as the supply of loose superficial material becomes depleted. Modelled erosion response was improved by allowing road erodibility to vary during an event. Changing the model values of erosion detachment parameters in response to changes in surface sediment availability improved model accuracy of predicted sediment transport by 30–40%. A predictive relationship between road erodibility ‘states’ and road surface sediment depth is presented. This relationship allows implementation of the dynamic erodibility (DE) method to events where pre‐storm sediment depth can be estimated (e.g., from traffic usage variables). Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluating and modelling the hydrological and erosive behaviour of an olive orchard microcatchment under no‐tillage with bare soil in Spain

Despite the high risk of erosion in olive orchards located in mountainous areas in Spain, little research has been carried out to account for the complexity and interaction of the natural processes of runoff and soil erosion on the catchment scale or small catchment scale. In this study, a microcatchment of 6·7 ha in a mountainous area under no‐tillage farming with bare soil was set up to record runoff and sediment. Soil erosion and runoff patterns were monitored over a two‐year period. Totally, 22 events were observed. The data were analysed, and then used to calibrate the AnnAGNPS model, which allowed us to complete the data period and describe the hydrological and erosive behaviour on a monthly and annual basis. A high variability in catchment responses was observed, due to differences in the storms and to the effect of the surface soil moisture content. Maximum intensities of 10 and 30 min determined the final runoff values while the total sediment loads were dependent on the rainfall depth. The impact of management on the reduction of porosity can explain the relationship between runoff and intensity in the microcatchment. However, the impact of the spatial scale meant that the transport of sediment required substantial rainfall depths to ensure a continuous flow from the hillslopes. The results of the calibration (E >0·60 and r >0·75) on the event and monthly scale confirmed the applicability of AnnAGNPS to predict runoff and erosion in the microcatchment. The predicted average runoff coefficient was 3·3% for the study period and the total average sediment loads, 1·3 Mg/ha/yr. Despite these low values, the model simulation showed that much larger runoff coefficients and soil losses can be expected for periods with several consecutive years in which the annual rainfall depth was over 500 mm. The use of cover is recommended to prevent the high levels of erosion associated with these conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Incremental distributed modelling investigation in a small agricultural catchment: 2. Erosion and phosphorus transport

Distributed physically based erosion and phosphorus (P) transport models, run by the overland flow model described in Taskinen and Bruen (2006. Hydrological Processes 20 : this issue), are described. In the erosion model, the additional components to the basic model were the outflow of the particles by infiltration and a new model component, i.e. deposition when rainfall stops. Two ways of calculating the shielding factor due to the flow depth were compared. The P transport model had both dissolved P (DP) and particulate P (PP) components. The processes included in the DP model were desorption from the soil surface, advection, storage in the overland flow and infiltration. The PP model accounted for advection, storage in the flow, infiltration, detachment from the soil surface by flow and rainfall and deposition both when transport capacity of suspended solids (SS) is exceeded and when rainfall ceases. When the models were developed and validated in small agricultural fields of cohesive soil types in southern Finland, comparisons were made between corresponding processes and the significance of added components were estimated in order to find out whether increased model complexity improves the model performance. The sedigraphs were found to follow the dynamics of rainfall, emphasizing the importance of the rainfall splash component. The basic model was too slow to react to changes in rainfall and flow rates, but infiltration and deposition that acts during the cessation in rainfall improved the model significantly by enabling the modelled SS to fall sharply enough. The shielding effect of flow depth from the splash detachment was found to play a significant role. Transport capacity should also be included in erosion models when they are applied to cohesive soils. In this study, the Yalin method worked well. A strong correlation was obvious between the measured SS and total P concentrations, indicating that the main form of P in runoff is PP. This emphasizes the importance of a good sediment transport model in P transport modelling. The submodel used for DP desorption from the soil surface produced plausible results without any calibration. Copyright © 2006 John Wiley & Sons, Ltd.