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.     

Evaluation of an interrill soil erosion model using laboratory catchment data

Physically based soil erosion simulation models require input parameters of soil detachment and sediment transport owing to the action and interactions of both raindrops and overland flow. A simple interrill soil water transport model is applied to a laboratory catchment to investigate the application of raindrop detachment and transport in interrill areas explicitly. A controlled laboratory rainfall simulation study with slope length simulation by flow addition was used to assess the raindrop detachment and transport of detached soil by overland flow in interrill areas. Artificial rainfall of moderate to high intensity was used to simulate intense rain storms. However, experiments were restricted to conditions where rilling and channelling did not occur and where overland flow covered most of the surface. A simple equation with a rainfall intensity term for raindrop detachment, and a simple sediment transport equation with unit discharge and a slope term were found to be applicable to the situation where clear water is added at the upper end of a small plot to simulate increased slope length. The proposed generic relationships can be used to predict raindrop detachment and the sediment transport capacity of interrill flow and can therefore contribute to the development of physically‐based erosion models. Copyright © 1999 John Wiley & Sons, Ltd.     

On predicting upland erosion losses from rainfall depth

Point rainfall triggers the complex processes of overland flow and surface erosion. The probability density functions of rainfall duration and intensity are coupled with a physically based dynamic formulation of rainfall-runoff-sediment transport relationships for upland areas. When considering a single storm, rainfall depth alone is a poor predictor of sediment transport because of the dispersion introduced by the effect of rainfall intensity. On a long terms basis, however, the total amount of rainfall can be used to predict total erosion losses.     

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.     

LISEM: A SINGLE-EVENT,PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. II: SENSITIVITY ANALYSIS,VALIDATION AND APPLICATION

A new hydrological and soil erosion model has been developed and tested: LISEM, the Limburg soil erosion model. The model uses physically based equations to describe interception, infiltration and soil water transport, storage in surface depressions, splash and flow detachment, transport capacity and overland and channel flow. From the validation results it is clear that, although the model has several advantages over other models, the results of LISEM 1.0 are far from perfect. Based on the sensitivity analysis and field observations, the main reasons for these differences seems to be the spatial and temporal variability of the soil hydraulic conductivity and the initial pressure head at the basin scale. Another reason for the differences between measured and simulated results is our lack or understanding of the theory of hydrological and soil erosion processes.     

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.     

Simulation of erosion-deposition processes at basin scale by a physically-based mathematical model

The development and application of the physically-based and spatially-distributed mathematical model CTSS8-SED is presented. The model simulates hydrologic-hydraulic processes produced by storm events and related soil erosion and sediment transport processes at basin scale in lowland areas. The model simulates (i) storm runoff,(ii) soil detachment by raindrop impact and overland flow (gross sediment yield),(iii) sediment transport by overland flow and associated erosion-deposition processes and (iv) sediment transport by stream flow and riverbed erosion-deposition processes. A quasi two-dimensional representation of water flow and sediment transport routing is made by means of interconnected cells approach. The model is applied to simulate two flooding events in the Luduea Creek basin (Santa Fe, Argentina) occurred in April 1994 and March 2007 due to extraordinary rainfalls.     

Effects of rainfall,overland flow and their interactions on peatland interrill erosion processes

Interrill erosion processes on gentle slopes are affected by mechanisms of raindrop impact, overland flow and their interaction. However, limited experimental work has been conducted to understand how important each of the mechanisms are and how they interact, in particular for peat soil. Laboratory simulation experiments were conducted on peat blocks under two slopes (2.5° and 7.5°) and three treatments: Rainfall, where rainfall with an intensity of 12 mm h^?1 was simulated; Inflow, where upslope overland flow at a rate of 12 mm h^?1 was applied; and Rainfall + Inflow which combined both Rainfall and Inflow. Overland flow, sediment loss and overland flow velocity data were collected and splash cups were used to measure the mass of sediment detached by raindrops. Raindrop impact was found to reduce overland flow by 10 to 13%, due to increased infiltration, and reduce erosion by 47% on average for both slope gradients. Raindrop impact also reduced flow velocity (80–92%) and increased roughness (72–78%). The interaction between rainfall and flow was found to significantly reduce sediment concentrations (73–85%). Slope gradient had only a minor effect on overland flow and sediment yield. Significantly higher flow velocities and sediment yields were observed under the Rainfall + Inflow treatment compared to the Rainfall treatment. On average, upslope inflow was found to increase erosion by 36%. These results indicate that overland flow and erosion processes on peat hillslopes are affected by upslope inflow. There was no significant relationship between interrill erosion and overland flow, whereas stream power had a strong relationship with erosion. These findings help improve our understanding of the importance of interrill erosion processes on peat. Copyright © 2017 John Wiley & Sons, Ltd.     

Overland runoff erosion dynamics on steep slopes with forages under field simulated rainfall and inflow

Although numerous studies have acknowledged that vegetation can reduce erosion, few process-based studies have examined how vegetation cover affect runoff hydraulics and erosion processes. We present field observations of overland flow hydraulics using rainfall simulations in a typical semiarid area in China. Field plots (5 × 2 m²) were constructed on a loess hillslope (25°), including bare soil plot as control and three plots with planted forage species as treatments—Astragalus adsurgens, Medicago sativa and Cosmos bipinnatus. Both simulated rainfall and simulated rainfall + inflow were applied. Forages reduced soil loss by 55–85% and decreased overland flow rate by 12–37%. Forages significantly increased flow hydraulic resistance expressed by Darcy–Weisbach friction factor by 188–202% and expressed by Manning's friction factor by 66–75%; and decreased overland flow velocity by 28–30%. The upslope inflow significantly increased overland flow velocity by 67% and stream power by 449%, resulting in increased sediment yield rate by 108%. Erosion rate exhibited a significant linear relationship with stream power. M. sativa exhibited the best in reducing soil loss which probably resulted from its role in reducing stream power. Forages on the downslope performed better at reducing sediment yield than upslope due to decreased rill formation and stream power. The findings contribute to an improved understanding of using vegetation to control water and soil loss and land degradation in semiarid environments.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

MODELING HEADCUT DEVELOPMENT AND MIGRATION IN UPLAND CONCENTRATED FLOWS

On hillslopes and agricultural fields, discrete areas of intense, localized soil erosion commonly take place in the form of migrating headcuts. These erosional features significantly increase soil loss and landscape degradation, yet the unsteady, transient, and migratory habits of headcuts complicate their phenomenological and erosional characterization. Here a unique experimental facility was constructed to examine actively migrating headcuts typical of upland concentrated flows. Essential components of the facility include a deep soil cavity with external drainage, rainfall simulator, capacity for overland flow, and a video recording technique for data collection. Results from these experiments show that： （1） after a short period of adjustment, headcut migration attained a steady-state condition, where the rate of migration, scour hole geometry, and sediment discharge remain constant with time; （2） boundary conditions of higher rates of overland flow, steeper bed slopes, and larger initial headcut heights produced systematically larger scour holes with higher rates of soil erosion; and （3） during migration, the turbulent flow structure within the scour hole remained unchanged, consisting of an overfall nappe at the brink transitioning into a reattached wall jet with two recirculation eddies within the plunge pool. The systematic behavior of headcut development and migration enabled the application of modified jet impingement theory to predict with good success the characteristics of the impinging jet, the depth of maximum scour, the rate of headcut migration, and the rate of sediment erosion. These laboratory data and the analytical formulation can be used in conjunction with soil erosion prediction technology to improve the management of agricultural areas impacted by headcut development and ephemeral gully erosion.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Predicting sediment transport by interrill overland flow on rough surfaces

Modelling soil erosion requires an equation for predicting the sediment transport capacity by interrill overland flow on rough surfaces. The conventional practice of partitioning total shear stress into grain and form shear stress and predicting transport capacity using grain shear stress lacks rigour and is prone to underestimation. This study therefore explores the possibility that inasmuch as surface roughness affects flow hydraulic variables which, in turn, determine transport capacity, there may be one or more hydraulic variables which capture the effect of surface roughness on transport capacity suffciently well for good predictions of transport capacity to be achieved from data on these variables alone. To investigate this possibility, regression analyses were performed on data from 1506 flume experiments in which discharge, slope, water temperature, rainfall intensity, and roughness size, shape and concentration were varied. The analyses reveal that 89·8 per cent of the variance in transport capacity can be accounted for by excess flow power and flow depth. Including roughness size and concentration in the regression improves that explained variance by only 3·5 per cent. Evidently, flow depth, when used in combination with excess flow power, largely captures the effect of surface roughness on transport capacity. This finding promises to simplify greatly the task of developing a general sediment equation for interrill overland flow on rough surfaces. Copyright © 1998 John Wiley & Sons, Ltd.     

Sediment particle selectivity and its response to overland flow hydraulics within grass strips

Particle selectivity plays an important role in clarifying sediment transport processes in vegetative filter strips (VFS). 10-m long grass strips at slopes of 5° and 15° were subjected to a series of sediment-laden inflows experiments with different particle sizes to investigate the sediment transport and its response to overland flow hydraulics. The inflow sediments came from local soil, river-bed sand, and mixed, with median particle size d₅₀ of 39.9, 207.9 and 77.4 μm, respectively. Three independent repeated experiments were carried for each treatment. The results show that when the sediment trapping lasted for a certain length of time, the re-entrainment of some small-sized particles was greater than the deposition; that is, net loss occurred, which was not erosion of the original soil. Net loss of particles is mainly determined by the particle diameter. The coarser the inflow sediment particles and/or the steeper the slope, the coarser the particles can be net lost. Deposited sediment causes the VFS bed surface to become smooth and hydraulic resistance decrease exponentially. Unit stream power P is more suitable than shear stress τ of overland flow to be used to describe the process of sediment particle transport in VFS. The relationship between P and d₅₀ of outflow sediment is very consistent with the form of power function with a constant term. These results are helpful to understand the physical process of sediment transport on vegetation hillslopes.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

The effect of antecedent moisture conditions on sediment and phosphorus loss during overland flow: Mahantango Creek catchment,Pennsylvania, USA

The loss of P in overland flow from most cultivated soils is controlled by erosion, and in‐turn soil moisture. We evaluated the effect of soil moisture on erosion and P transport in overland flow by applying rainfall (7 cm h^?1) to packed soil boxes (1 m long and 0·15 m wide) and field plots (1 and 10 m long by 1 m wide) of silt loams in a central Pennsylvania (USA) catchment. Flow from packed soil boxes took longer to initiate as antecedent soil moisture decreased from field capacity (2 min) to air dried (8 to 9 min). Even in the more complex field plots (i.e. soil heterogeneity and topography), the wetter site (1 by 10 m plot; 70% field capacity) produced flow more quickly (3 min) and in greater volume (439 L) than the drier site (1 by 10 m plot; 40% field capacity, 15 min, and 214 L, respectively). However, less suspended sediment was transported from wetter soil boxes (1·6 to 2·5 g L^?1) and field plots (0·9 g L^?1) than drier boxes (2·9 to 4·2 g L^?1) and plots (1·2 g L^?1). Differences are attributed to their potential for soil aggregate breakdown, slaking and dispersion, which contribute to surface soil sealing and crusting, as dry soils are subject to rapid wetting (by rainfall). During flow, selective erosion and antecedent moisture conditions affected P transport. At field capacity, DRP and PP transport varied little during overland flow. Whereas P transport from previously dry soil decreased rapidly after the initiation of flow (6 to 1·5 mg TP L^?1), owing to the greater slaking and dispersion of P‐rich particles into flow at the beginning than end of the flow event. These results indicate that soil moisture fluctuations greatly effect erosion and P transport potential and that management to decrease the potential for loss should consider practices such as conservation tillage and cover crops, particularly on areas where high soil P and erosion coincide. Copyright © 2002 John Wiley & Sons, Ltd.     

Estimating overland flow erosion capacity using unit stream power

Soil erosion caused by water flow is a complex problem. Both empirical and physically based approaches were used for the estimation of surface erosion rates. Their applications are mainly limited to experimental areas or laboratory studies. The maximum sediment concentration overland flow can carry is not considered in most of the existing surface erosion models. The lack of erosion capacity limitation may cause over estimations of sediment concentration. A correlation analysis is used in this study to determine significant factors that impact surface erosion capacity. The result shows that the unit stream power is the most dominant factor for overland flow erosion which is consistent with experimental data. A bounded regression formula is used to reflect the limits that sediment concentration cannot be less than zero nor greater than a maximum value. The coefficients used in the model are calibrated using published laboratory data. The computed results agree with laboratory data very well. A one dimensional overland flow diffusive wave model is used in conjunction with the developed soil erosion equation to simulate field experimental results. This study concludes that the non-linear regression method using unit stream power as the dominant factor performs well for estimating overland flow erosion capacity.     

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.     

Impact of rainfall pattern on interrill erosion process

The impact of rainfall pattern on the interrill erosion process is not fully understood despite its importance. Systematic rainfall simulation experiments involving various rainfall intensities, stages, intensity sequences, and surface cover conditions were conducted in this study to investigate their effects on the interrill erosion process. Five rainfall patterns designed with the same total kinetic energy/precipitation (increasing, decreasing, rising–falling, falling–rising and constant patterns) were randomly delivered to a pre‐wet clay loam soil surface at a 10° slope gradient. Significant differences in soil losses were observed among the different rainfall patterns and stages, but there was no obvious difference in runoff. Kinetic energy flux (KE_r) was a governing factor for interrill erosion, and constant rainfall pattern (CST) produced nine times greater soil loss than runs with no KE_r. Varied‐intensity patterns had a profound effect on raindrop‐induced sediment transport processes; path analysis results indicated that said effect was complex, interactive and intensity‐dependent. Low hydraulic parameter thresholds further indicated that KE_r was the dominant factor in detaching soil particles, while overland flow mainly contributed to transporting the pre‐detached particles. This study not only sheds light on the mechanism of interrill sediment transport capacity and detachability, but also may provide a useful database for developing event‐based interrill erosion prediction models. Copyright © 2017 John Wiley & Sons, Ltd.     

LISEM: A SINGLE-EVENT PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. I: THEORY,INPUT AND OUTPUT

A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.