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.     

EVALUATION OF VERTICAL AND LATERAL FLOW THROUGH AGRICULTURAL LOESSIAL HILLSLOPES USING A TWO-DIMENSIONAL COMPUTER SIMULATION MODEL

On four hillslopes in the loess region of The Netherlands pressure heads were monitored using stand alone measuring devices. During rain events pressure heads were measured regularly with time intervals of five minutes. Rainfall was recorded with a tipping bucket rainfall gauge. Two-dimensional simulations were executed to simulate water flow through these hillslopes during erosive rain events in summer and winter periods. From these simulations it appeared that vertical flow is dominant during rain events, and that lateral water transport is of minor importance despite distinct layer stratifications. During the selected rain showers, average lateral water movement varied between 1.6 and 4.7% of the total water displacement through the hillslopes. Therefore, it was decided to incorporate a one-dimensional water flow module into the event-based distributed soil erosion and hydrological model LISEM.     

Evaluation of the LISEM soil erosion model in two catchments in the East African Highlands

Under increasing population pressure, soil erosion has become a threat in the East African Highlands, and erosion modelling can be useful to quantify this threat. To test its applicability for this region, the LISEM soil erosion model was applied to two small catchments, one in the Usumbara Mountains, Tanzania, and the other on the slopes of Mount Kenya. Input data for the model were collected in both catchments, as were data on runoff and erosion that were used for calibration and validation of the model. LISEM was first calibrated on catchment outlet data, and afterwards simulated spatial patterns of erosion were compared to available erosion data. The results showed that LISEM can, after calibration, give good discharge predictions for some events, but not for all. However, LISEM generally overpredicted soil loss from the catchments. Comparison with observed erosion patterns did not show overprediction, but according to the model, erosion was more widespread than was observed. There are several reasons for these discrepancies. First, it is difficult to obtain enough accurate data to run the model, such as accurate maps, rainfall data and soil and plant characteristics. Second, it is also difficult to obtain accurate data to evaluate the performance of the model, either for the catchment outlet or spatially, therefore observed erosion rates are also uncertain. Third, the model could not deal correctly with complex events, i.e. those having double rainfall peaks, and might also have difficulties with catchment characteristics such as soil type and the complexity of land use. Finally, LISEM could not deal with events in which throughflow or baseflow played a role, which was to be expected since those processes are not simulated by LISEM. Nevertheless, LISEM could be calibrated to give good discharge predictions for some events, and also gave reasonable results when compared to data obtained from erosion plots. Furthermore, only complex, distributed, storm‐based models such as LISEM can give spatial predictions for single storms. Therefore, it is concluded that if the aim is spatial prediction on an event basis, there is no alternative to complex erosion models such as LISEM, but if the aim is to predict average annual erosion, the data‐demanding, physically based LISEM erosion model may not be the most appropriate model. Copyright © 2005 John Wiley & Sons, Ltd.     

RUNOFF GENERATION AND SOIL EROSION IN SMALL AGRICULTURAL CATCHMENTS WITH LOESS-DERIVED SOILS

Over a two-year period, rainfall, runoff and sediment output were measured in six small agricultural catchments (3–10 ha) in south Limburg (The Netherlands). These measurements were needed for validation of an erosion model for South Limburg (LISEM). In this paper, results of the measurements are presented and processes that determine surface runoff and sediment yield during winter and summer rainfall are identified. Before the start of the measurement programme, surface slaking and crust formation on the erodible loess soils were thought to be the main cause of overland flow and soil erosion in South Limburg. This was the starting point for soil conservation measures in the area. The measurement results discussed in this paper show that in some catchments much runoff occurred in winter and that soil moisture storage capacity may be just as important for runoff generation as infiltration capacity. Therefore, when modelling soil erosion and optimizing erosion control measures for South Limburg, runoff generation through Hortonian as well as through saturation overland flow must be considered.     

Effects of grid cell size and time step length on simulation results of the Limburg soil erosion model (LISEM)

With increasing computer power, process‐based models that use grids to discretize space have become increasingly popular. For such models, the simulation results might depend on both grid cell size and, in the case of dynamic models, on the time step length used in the model. In this study, the dynamic Limburg soil erosion model (LISEM) was applied to a small catchment on the Chinese Loess Plateau. To study the effect of grid cell size and time step length, simulations were performed for grid cell sizes ranging from 5 to 100 m for a single time step length, and for time step lengths ranging from 2 to 120 s for a single grid cell size. The results show that the LISEM results vary considerably as a function of both grid cell size and time step length. For both increase in cell size and increasing time step length, the trend was a decrease in predicted discharge and predicted soil loss. For discharge, the most important causes are likely to be a decrease in slope with increasing grid cell size, rainfall averaging for longer time step lengths, and numerical dispersion of the kinematic wave solution. For soil loss, the cause is less clear, reflecting the complexity of soil loss prediction, which depends on available water, transport capacity and sediment redistribution, all of which change in time and space. These results show that a choice for a certain grid cell size and a certain time step length should be made before calibration of the model. Similar erosion models are likely to have similar dependencies on grid size and time step length. Copyright © 2005 John Wiley & Sons, Ltd.     

Erosion characteristics of loess tunnels on the Loess Plateau: A field investigation and experimental study

Loess tunnels are a common geo-hazard in the Loess Plateau and not only cause considerable soil and water loss, but also aggravate and even induce the occurrence and development of other disasters such as ground fissures, mudflows, collapses, and landslides. To date, research on the hydrological characteristics and erosion behaviour of loess tunnel systems has focused on field investigation data and limited river basin observation data, whereas field test information and data are very scarce. In this study, field surveys, observations, field scouring experiments, and laboratory-based sediment percentage tests were conducted to analyse the erosion characteristics, spatial distribution, and hydrological characteristics of a large-scale loess tunnel system in the hilly Loess Plateau southeast region of northern China. The results showed that the loess tunnel erosion exhibited periodicity. Tunnel erosion in each period shows a similar erosion process, that is, thin-layer water flow erosion and lateral expansion, tunnel wall collapse and deposit due to the loss of support from the lateral erosion, and erosion and transport of deposits by water flow. Waterfall erosion, lateral erosion, headward erosion, and the resulting collapses were the main forms of tunnel erosion. Besides this, the base level of erosion significantly affects the erosion characteristics of the loess tunnel. The hydrological characteristics during field scouring experiments exhibited three different stages: a lag effect, attributed to the temporary loss of water velocity in the first stage; small water flow fluctuations in the second stage; and an increase in total seepage loss with increased water injection flow in the third stage. The erosion rate was positively correlated with the flow quantity. The results of this study not only provide valuable reference data for research on the mechanism and velocity of erosion events in loess, but also provide a theoretical basis for the prevention of loess tunnel disasters in engineering construction. © 2020 John Wiley & Sons, Ltd.     

MATHEMATICAL MODEL OF OVERLAND FLOW AND MECHANISM OF SOIL CONSERVATION FOR FORESTED STEEP HILLSLOPE－－（II）Mechanism of Soil Conservation Over Steep Hillslope

ＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＯＶＥＲＬＡＮＤＦＬＯＷＡＮＤＭＥＣＨＡＮＩＳＭＯＦＳＯＩＬＣＯＮＳＥＲＶＡＴＩＯＮＦＯＲＦＯＲＥＳＴＥＤＳＴＥＥＰＨＩＬＬＳＬＯＰＥ－－（ＩＩ）ＭｅｃｈａｎｉｓｍｏｆＳｏｉｌＣｏｎｓｅｒｖａｔｉｏｎＯｖｅｒＳｔ...     

An experimental investigation of the effects of thawed soil depth on rill flow velocity

Water flow velocity is an important hydraulic variable in hydrological and soil erosion models, and is greatly affected by freezing and thawing of the surface soil layer in cold high-altitude regions. The accurate measurement of rill flow velocity when impacted by the thawing process is critical to simulate runoff and sediment transport processes. In this study, an electrolyte tracer modelling method was used to measure rill flow velocity along a meadow soil slope at different thaw depths under simulated rainfall. Rill flow velocity was measured using four thawed soil depths (0, 1, 2 and 10 cm), four slope gradients (5°, 10°, 15° and 20°) and four rainfall intensities (30, 60, 90 and 120 mm·h⁻¹). The results showed that the increase in thawed soil depth caused a decrease in rill flow velocity, whereby the rate of this decrease was also diminishing. Whilst the rill flow velocity was positively correlated with slope gradient and rainfall intensity, the response of rill flow velocity to these influencing factors varied with thawed soil depth. The mechanism by which thawed soil depth influenced rill flow velocity was attributed to the consumption of runoff energy, slope surface roughness, and the headcut effect. Rill flow velocity was modelled by thawed soil depth, slope gradient and rainfall intensity using an empirical function. This function predicted values that were in good agreement with the measured data. These results provide the foundation for a better understanding of the effect of thawed soil depth on slope hydrology, erosion and the parameterization scheme for hydrological and soil erosion models.     

A GIS-BASED DISTRIBUTED SOIL EROSION MODEL： A CASE STUDY OF TYPICAL WATERSHED, SICHUAN BASIN

Based on the measuring data and Digital Elevation Data (DEM) in a typical watershed--Hemingguan Watershed, Nanbu County, Sichuan Province of China, a GIS-based distributed soil erosion model was developed particularly for the purple soil type. It takes 20 m × 20 m grid as calculating unit and operates at 10-minute time interval. The required input data to the model include DEM, soil, land use, and time-series of precipitation and evaporation loss. The model enables one to estimate runoff, erosion and sediment yield for each grid cell and route the flow along its flow path to the watershed outlet. Furthermore, the model is capable of calculating the total runoff; erosion and sediment yield for the entire watershed by recursion algorithm. The validation of the model demonstrated that it could quantitatively simulate the spatial distribution of hydrological variables in a watershed, such as runoff, vegetation entrapment, soil erosion, the degree of soil and water loss. Moreover, it can evaluate the effect of land use change on the runoff generation and soil erosion with an accuracy of 80% and 75% respectively. The application of this model to a neighboring watershed with similar conditions indicates that this distributed model could be extended to other similar regions in China.     

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.     

SPATIAL AND TEMPORAL VARIABILITY OF SOIL SURFACE ROUGHNESS AND THE APPLICATION IN HYDROLOGICAL AND SOIL EROSION MODELLING

Accurate estimations of water retention and detention are needed to simulate surface runoff and soil erosion following a rainfall event in a catchment. Several equations to estimate the amount of surface depressional storage, the fraction of the soil surface covered by water and the amount of rainfall excess needed to start surface runoff have been developed by Onstad (1984). The random roughness and slope gradient are needed for those estimations. Surface micro-elevation data have been gathered by a photographic method. The random roughness was determined from those elevation measurements. Several factors which have an impact on the soil surface roughness were taken into account. The main sources of influence are the type of land use, the crop stage within the growing period and tillage direction. Analyses of variance indicated that the variation in the RR-index could be explained mainly by type of land use, orientation and field type. The temporal variation was relatively small. Gradient data have been determined from a digital elevation model, constructed by digitizing contours. Combining the random roughness and the steepness of slope, the amounts of surface water retention and detention could be estimated. Knowledge of water retention and detention will improve the estimations of runoff and soil erosion modelling in catchments, such as those made with the LISEM model. The agricultural systems examined in this study have similar random roughness values in summer. Different soil erosion rates for several types of land use can not therefore be explained by the random roughness.     

Hydrological simulation of a conservation bench terrace system in a subhumid climate

Abstract

A finite element model to simulate runoff and soil erosion from agricultural lands has been developed. The sequential solutions of the governing differential equations were found: Richards' equation with a sink term for infiltration and soil water dynamics under cropped conditions; St Venant equation with kinematic wave approximation for overland and channel flow; and sediment continuity equation, for soil erosion. The model developed earlier has been improved to simulate erosion/deposition in impoundments and predicted and observed soil loss values were in reasonably good agreement when the model was tested for a conservation bench terrace (CBT) system. The finite element model was extensively applied to study the hydrological behaviour of a CBT system vis-à-vis the conventional system of sloping borders. The model estimates runoff and soil loss reasonably well, under varying conditions of rainfall and at different crop growth stages. The probable reasons for discrepancies between observation and simulation are reported and discussed. Sensitivity analysis was carried out to study the effect of various hydrological, soil and topographical parameters, such as ratio of contributing to receiving areas, weir length, depth of impoundment, slope of contributing area, etc. on the flow behaviour in a CBT system.     

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.     

Rapid assessment of field erosion and sediment transport pathways in cultivated catchments after heavy rainfall events

This paper explores a scale‐adapted erosion mapping method which aims at a rapid assessment of field erosion and sediment transport pathways in catchments up to several square kilometres and compares the results with the output of a well‐known erosion model (LISEM). The mapping method is based on an event‐defined classification scheme of erosion intensity (zero, weak, moderate and strong) that is applied to arable fields, in combination with incision measurements of erosion features for each erosion intensity class on a small sample of fields. Sediment deposition is classified on the basis of quantity indicators and abundance. In addition, relevant conditions and erosion factors are determined for each field. The method was applied to an agricultural catchment (4·2 km²) in the Sundgau (Alsace), after a short but violent thunderstorm in May 2001, to illustrate its potential use and its limitations. The rainfall event led to strong erosion on the arable fields and a muddy flow that caused significant damage in the built‐up area. On the basis of the analyses of the incision measurements in combination with the mapping of erosion intensity classes, total erosion for the catchment was estimated as 15 000 t (an average of about 36 t[sol ]ha). Sediment deposition was found to occur in three major locations: (1) in thalwegs at the interface between maize and downslope winter wheat fields, (2) in downslope headlands where the flow direction suddenly changed due to oriented tillage structures in the perpendicular direction, and (3) the lowest corners of fields which collect all the runoff from the field. Preliminary data analyses suggest that erosion intensity is related to field size and[sol ]or tillage direction and to slope morphology. Model output (LISEM) appeared to depend more strongly on slope gradient than the results obtained with the mapping method. The method yields a database, which can be used as a foundation for conservation strategies in small regions with similar land use and geomorphology. The mapping and modelling methods are compared, and their complementary aspects are highlighted. Copyright © 2005 John Wiley & Sons, Ltd.     

Negligible soil erosion in a burned mountain watershed,Canadian Rockies: field and modelling investigations considering the role of duff

Increased soil erosion in immediate post‐wildfire years has been well documented in the literature, but many unanswered questions remain about the factors controlling erosional responses in different regional settings. The field site for the present study was located in a closed canopy, subalpine forest in Kootenay National Park, British Columbia that was subjected to a high‐intensity crown fire in the summer of 2003. Low soil erosion values were documented at the study site in the years immediately following the 2003 wildfire, with estimates ranging from approximately 10^‐1 up to 10⁰ t ha^‐1. Following the wildfire, notable duff coverage (the duff layer is the combined fermentation and humus soil organic layers) remained above the mineral soil. This finding supports earlier studies documenting only partial duff consumption by high‐intensity wildfires in the boreal forest of Canada. It is postulated that remnant duff coverage after many high‐intensity wildfires impacts the hydrological and soil erosional response to rainstorm events in post‐wildfire years. In particular, duff provides detention storage for infiltrating rainfall and, therefore, may inhibit the generation of overland flow. Furthermore, duff also provides a physical barrier to soil erosion. The Green–Ampt model of rainfall infiltration is employed to better assess how interactions between rainfall duration/intensity and soil/duff properties affect hydrological response and the generation of overland flow. Model results show that duff provides an effective zone for detention storage and that duff accommodates all rainfall intensities to which it was subjected without the occurrence of surface ponding. In addition, the penetration of the wetting front is relatively slow in duff due to its high porosity and water storage potential. Copyright © 2011 John Wiley & Sons, Ltd.     

PROCESSES OF HEADCUT GROWTH AND MIGRATION IN RILLS AND GULLIES

lINTRoDUCTIONHeadcuterosionwithinrills,ephemeralguIlies,classicgullies,andstreamscausesseriousenvironmentalproblems.Headcuterosionacceleratesthelossoftopsoilanddecreasestheproductivityofagriculturallands.Erodedsedimentsoftenendupinreceivingstreams,causingwaterqualityproblemsandnegativelyimpactingbioIogicalprocesses.Inadditiontolandscapedegradation,gulliesarethedominantformofdamagetoearthspillways.lfagullycanmovethroughanearthspillwayandbreachthecrestofadam,thentheimpoundedfloodwaterswiI1…     

Modelling increased soil cohesion due to roots with EUROSEM

As organic root exudates cause soil particles to adhere firmly to root surfaces, roots significantly increase soil strength and therefore also increase the resistance of the topsoil to erosion by concentrated flow. This paper aims at contributing to a better prediction of the root effects on soil erosion rates in the EUROSEM model, as the input values accounting for roots, presented in the user manual, do not account for differences in root density or root architecture. Recent research indicates that small changes in root density or differences in root architecture considerably influence soil erosion rates during concentrated flow. The approach for incorporating the root effects into this model is based on a comparison of measured soil detachment rates for bare and for root‐permeated topsoil samples with predicted erosion rates under the same flow conditions using the erosion equation of EUROSEM. Through backwards calculation, transport capacity efficiencies and corresponding soil cohesion values can be assessed for bare and root‐permeated topsoils respectively. The results are promising and present soil cohesion values that are in accordance with reported values in the literature for the same soil type (silt loam). The results show that grass roots provide a larger increase in soil cohesion as compared with tap‐rooted species and that the increase in soil cohesion is not significantly different under wet and dry soil conditions, either for fibrous root systems or for tap root systems. Power and exponential relationships are established between measured root density values and the corresponding calculated soil cohesion values, reflecting the effects of roots on the resistance of the topsoil to concentrated flow incision. These relationships enable one to incorporate the root effect into the soil erosion model EUROSEM, through adapting the soil cohesion input value. A scenario analysis shows that the contribution of roots to soil cohesion is very important for preventing soil loss and reducing runoff volume. The increase in soil shear strength due to the binding effect of roots on soil particles is two orders of magnitude lower as compared with soil reinforcement achieved when roots mobilize their tensile strength during soil shearing and root breakage. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling the process of soil detachment by rill flow on steep loessial hillslopes

Soil detachment by rill flow is a key process of rill erosion, modelling this process can help in understanding rill erosion mechanisms. However, many soil detachment models are established on conceptual assumptions rather than experimental data. The objectives of this study were to establish a model of soil detachment by rill flow based on flume experimental data and to quantitatively verify the model. We simulated the process of soil detachment by rill flow in flume experiments with a soil-feeding hopper using loessial soil on steep slopes. Seven flow discharges, six slopes and five sediment loads were combined. Soil detachment capacity, sediment transport capacity, and soil detachment rate by rill flow under different sediment loads were measured. The process of soil detachment by rill flow can be modelled by a dual power function based on soil detachment capacity and transport capacity deficit as variables. The established model exhibited high credibility (NSE=0.97; R²=0.97). The contributions of soil detachment capacity and transport capacity deficit to soil detachment rate by rill flow reached 60% and 36%, respectively. Soil detachment capacity exerted more influence on soil detachment rate than did transport capacity deficit. The performance of the WEPP rill erosion equation is also favourable (NSE=0.95; R²=0.97). The two power exponents in the model we established strengthen the role of soil detachment capacity in soil detachment rate and weaken that for transport capacity deficit. Soil detachment capacity and transport capacity deficit played important roles in the determination of soil detachment rate by rill flow. The results can be applied to implement the numerical modeling and prediction of rill erosion processes on steep loessial hillslopes. © 2019 John Wiley & Sons, Ltd.     

SOIL EROSION PROCESS RESEARCH AND ITS POTENTIAL IMPACT ON EROSION PREDICTION MODEL DEVELOPMENT

During the past 50 years, many research efforts have been invested in understanding soil erosion process and development of erosion prediction models at various scales. This paper briefly introduces the erosion process and prediction model development in the USA. Especially, this paper focuses on discussing potential impacts of the erosion process on erosion model development, and future directions of the soil erosion process research and process- based model development. 1 DEVELOPMENT O…     

Morphology,structure and flow phases in soil pipes developing in forested hillslopes underlain by a Quaternary sand–gravel formation,Hokkaido, northern main island in Japan

Preferential flow pathways, such as soil pipes, are usually present in the soil of slopes. Subsurface flow through the soil pipes affects the subsurface drainage system and is responsible for sediment removal from slopes. However, a record of the inner structure of soil pipes has rarely been reported for slopes. A fibrescope examination of the morphology and flow phases in soil pipes in hillslopes underlain by a Quaternary sand–gravel formation provided the following information: the main pores of the soil pipes ran mostly parallel with the slope gradient; the cross‐sections of the soil pipes were approximately circular; and occurred on a few occasions; with some triple junctions being present. In addition, both full flow and partly full‐depth conditions occurred simultaneously in the soil pipe. The full flow condition has long been used in hydrological studies to model the pipe flow mechanism. Both the full flow condition and the partly full‐depth condition, however, must be examined closely in order to evaluate the subsurface hydrology in heterogeneous soil and the hydrogeomorphological processes of subsurface hydraulic erosion. Copyright © 2000 John Wiley & Sons, Ltd.