Plants as river system engineers

Plants growing within river corridors both affect and respond to fluvial processes. Their above‐ground biomass modifies the flow field and retains sediment, whereas their below‐ground biomass affects the hydraulic and mechanical properties of the substrate and consequently the moisture regime and erosion susceptibility of the land surface. This paper reviews research that dates back to the 1950s on the geomorphological influence of vegetation within fluvial systems. During the late twentieth century this research was largely pursued through field observations, but during the early years of the twenty‐first century, complementary field, flume and theoretical/modelling investigations have contributed to major advances in understanding the influence of plants on fluvial systems. Flume experiments have demonstrated the fundamental role of vegetation in determining river planform, particularly transitions from multi‐ to single‐thread forms, and have provided insights into flow–vegetation–sediment feedbacks and landform building, including processes such as channel blockage and avulsion. At the same time, modellers have incorporated factors such as moisture‐dependent plant growth, canopy and root architecture and their influence on flow resistance and sediment/bank reinforcement into morphodynamic models. Meanwhile, field investigations have revealed that vegetation has a far more important and complex influence on fluvial systems than previously realized. It is now apparent that the influence of plants on river systems is significant across space scales from individual plants to entire forested river corridors. Small plant‐scale phenomena structure patch‐scale geomorphological forms and processes, and interactions between patches are almost certainly crucial to larger‐scale and longer‐term geomorphological phenomena. The influence of plants also varies continuously through time as above‐ and below‐ground biomass change within the annual growth cycle, over longer‐term growth trajectories, and in response to external drivers of change such as climatic, hydrological and fluvial fluctuations and extremes. Copyright © 2013 John Wiley & Sons, Ltd.     

From the blade scale to the reach scale: A characterization of aquatic vegetative drag

Previous studies have considered vegetative drag at different scales, the blade scale, the patch scale, and the reach scale, but few studies have considered the connection between these scales. In this paper, we develop simple, physically-based models that connect processes affecting the drag generated by aquatic vegetation at the blade and patch scale to the hydraulic resistance produced by vegetation at the reach scale. For fully developed flows through submerged patches of vegetation, velocities can be successfully predicted using a two-layer model in which momentum transfer from the unobstructed flow to the vegetation patch is characterized using a constant friction factor. To account for vegetation flexibility in this two-layer model, we develop an iterative procedure that calculates the reduction in plant height and drag for a given flow speed based on the plant material properties, and feeds this information back into the momentum balance. This simple iteration accurately predicts vegetation heights and velocities for submerged flexible vegetation. Finally, we consider the effect of varying vegetation distribution patterns by extending the two-layer model to account for more complex channel and patch geometries. The total hydraulic resistance produced by vegetation depends primarily on the blockage factor, i.e. the fraction of the channel cross-section blocked by vegetation. For a constant blockage factor, the specific distribution of vegetation can also play a role, with a large number of small patches generating more resistance than a single large patch. By considering models with different levels of complexity, we offer suggestions for what field measurements are needed to advance the prediction of channel resistance.     

The accuracy of drainage network delineation as a function of environmental factors: A case study in Central and Northern Sweden

Drainage networks delineated from Digital Elevation Models (DEMs), are the basis for the modelling of geomorphological and hydrological processes, biogeochemical cycling, and water resources management. Besides providing effective models of water flows, automatically extracted drainage networks based on topography can diverge from reality to varying degrees. The variability of such disagreement within catchments has rarely been examined as a function of the heterogeneity of land cover, soil type, and slope in the catchment of interest. This research gap might not only substantially limit our knowledge of the uncertainty of hydrological prediction, but can also cause problems for users attempting to use the data at a local scale. Using 1:100000 scale land cover maps, Quaternary deposits maps, and 2 m resolution DEMs, it is found that the accuracy of delineated drainage networks tends to be lower in areas with denser vegetation, lower hydraulic conductivity, and higher erodibility. The findings of this study could serve as a guide for the more thoughtful usage of delineated drainage networks in environmental planning, and in the uncertainty analysis of hydrological and biochemical predictions. Therefore, this study makes a first attempt at filling the knowledge gap described above.     

Sensitivity analysis of hydraulic conductivity and Manning's n parameters lead to new method to scale effective hydraulic conductivity across model resolutions

Hydrological modelling is an important tool for research, policy, and management, but uncertainty remains about parameters transferability from field observations made at small scale to models at the catchment scale and larger. This uncertainty compels the need to develop parameter relationships that are translatable across scale. In this study, we compare the changes to modelled processes as resolution is coarsened from 100‐m to 1‐km in a topographically complex, 255‐km² Colorado River headwater catchment. We conducted a sensitivity analysis for hydraulic conductivity (K) and Manning's n parameters across four orders of magnitude. Results showed that K acts as a moderator between surface and subsurface contributions to streamflow, whereas n moderates the duration of high intensity, infiltration‐excess flow. The parametric sensitivity analysis informed development of a new method to scale effective hydraulic conductivity across modelling resolutions in order to compensate for the loss of topographic gradients as resolution is coarsened. A similar mathematical relationship between n and lateral resolution changes was not found, possibly because n is also sensitive to time discretization. This research provides an approach to translate hydraulic conductivity parameters from a calibrated coarse model to higher resolutions where the number of simulations are limited by computational demand.     

Vegetation disturbances and flood energy during an extreme flood on a sub-tropical river

The advent of 2D hydraulic modelling has improved our understanding of flood hydraulics, thresholds, and dynamic effects on floodplain geomorphology and riparian vegetation at the morphological-unit scale. Hydraulic concepts of bed shear stress, stream power maxima, and energy (cumulative stream power) have been used to characterize floods and define their geomorphic effectiveness. These hydraulic concepts were developed in the context of reach-averaged, 1D hydraulic analyses, but their application to 2D model results is problematic due to differences in the treatment of energy losses in 1D and 2D analyses. Here we present methods for estimating total and boundary resistance from 2D modelling of an extreme flood on a subtropical river. Hydraulic model results are correlated with observations of the flood impacts on floodplain geomorphology and the riparian vegetation to identify thresholds and compute variants of flood energy. Comparison of LiDAR data in 2011 and 2014 shows that the 2011 flood produced 2–4 m of erosion on floodplain bars that were previously forested or grass-covered. Deposition on flood levees, dunes, and chute bars was up to 3.4 m thick. Various hydraulic metrics were trialled as candidates for thresholds of vegetation disturbance. The accuracy of thresholds using metrics extracted at the flood peak (i.e. boundary resistance and stream power maxima) was similar to that using energy as a threshold. Disturbance to forest and grass on vegetated bars was associated with stream powers of >834 W/m² and unit flows of >26 m²/s, respectively. Correlation of the hydraulic metrics with erosion and deposition depths showed no substantial improvement in using flood energy compared to metrics extracted at the flood peak for describing erosion and deposition. The extent of vegetation disturbances and morphological adjustments was limited for this extreme flood, and further 2D studies are needed to compare disturbance thresholds across different environments.     

Hydraulic conductivity in upland blanket peat: measurement and variability

A key parameter used in wetland hydrological and landform development models is hydraulic conductivity. Head recovery tests are often used to measure hydraulic conductivity, but the calculation techniques are usually confined to rigid soil theory. This is despite reports demonstrating the misapplication of rigid soil theory to non‐rigid soils such as peats. Although values of hydraulic conductivity calculated using compressible techniques have been presented for fenland peats, these data have never, to the authors' knowledge, been compared with such calculations in other peat types. Head recovery tests (slug withdrawal) were performed on piezometers at depths ranging from 10 to 80 cm from the surface on north Pennines blanket peats. Results were obtained using both rigid and compressible soil theories, thus allowing comparison of the two techniques. Compressible soil theory gives values for hydraulic conductivity that are typically a factor of five times less than rigid soil calculations. Hydraulic conductivity is often assumed to decrease with depth in upland peats, but at the study site in the northern Pennines it was not found to vary significantly with depth within the range of peat depths sampled. The variance within depth categories was not significantly different to the variance between depth categories showing that individual peat layers did not have characteristic hydraulic conductivity values. Thus, large lateral and vertical differences in hydraulic conductivity over short distances create problems for modelling but may help account for the high frequency of preferential flow pathways within what is otherwise a low matrix hydraulic conductivity peat. Hydraulic conductivity was found to vary significantly between sampling sites, demonstrating that hillslope‐ or catchment‐scale variability may be more important than plot‐scale variability. Values for compressibility of the peats are also reported. These generally decline with depth, and they also vary significantly between sampling sites. There are implications for the way in which measurements of hydraulic conductivity and other properties of blanket peat are interpreted, as the effects of environmental change in one part of a peat catchment may be very different to those in another. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of Vegetation Restoration on Soil Physical Properties in the Wind–Water Erosion Region of the Northern Loess Plateau of China

In the northern Loess Plateau that has been severely affected by wind–water erosion, shifts from arable land to forest or grasslands have been promoted since 1998, using both native and introduced vegetation. However, there is little knowledge of the ecological consequences and effectiveness of the vegetation restoration in the region. Therefore, relationships between watershed‐scale soil physical properties and plant recovery processes were analyzed. The results show that soil physical properties such as bulk density, hydraulic conductivity, mean weight diameter, and the stability of >1 mm macro‐aggregates have been significantly ameliorated in the 0–20 cm soil layer under secondary natural grasslands. In contrast, re‐vegetation with introduced species such as Caragana korshinskii or Medicago sativa had adversely affected the soil physical properties, probably due to the deterioration of soil water conditions and lower organic matter inputs resulting from severe erosion. Reductions in bulk density and increases in saturated hydraulic conductivity could be used as indicators of soil structure amelioration since they are closely related to most other measured properties. Practical considerations for future re‐vegetation projects are suggested, particularly that native species with lower water consumption rates than the introduced species should be used to avoid further soil degradation.     

Floodplain friction parameterization in two‐dimensional river flood models using vegetation heights derived from airborne scanning laser altimetry

Two‐dimensional (2‐D) hydraulic models are currently at the forefront of research into river flood inundation prediction. Airborne scanning laser altimetry is an important new data source that can provide such models with spatially distributed floodplain topography together with vegetation heights for parameterization of model friction. The paper investigates how vegetation height data can be used to realize the currently unexploited potential of 2‐D flood models to specify a friction factor at each node of the finite element model mesh. The only vegetation attribute required in the estimation of floodplain node friction factors is vegetation height. Different sets of flow resistance equations are used to model channel sediment, short vegetation, and tall and intermediate vegetation. The scheme was tested in a modelling study of a flood event that occurred on the River Severn, UK, in October 1998. A synthetic aperture radar image acquired during the flood provided an observed flood extent against which to validate the predicted extent. The modelled flood extent using variable friction was found to agree with the observed extent almost everywhere within the model domain. The variable‐friction model has the considerable advantage that it makes unnecessary the unphysical fitting of floodplain and channel friction factors required in the traditional approach to model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Runoff production processes in small alpine catchments within the unconsolidated Pleistocene sediments of the Lainbach area (Upper Bavaria)

In the Lainbach catchment, unconsolidated Pleistocene moraine sediments are widely distributed. Because of the great natural risk of floods, together with extreme loads of sediments, investigations of runoff production processes have been conducted in this area. At hillslope scale three test sites with different states of soil development and vegetation cover were instrumented with V‐shaped weirs, precipitation gauges and measurement devices for electrical conductivity (EC) of discharge water. The EC has been used as a geochemical tracer for hydrograph separation, since the statistical relationship between content of dissolved Ca²⁺, Mg²⁺ cations and EC is highly significant for different stages of runoff. This method allows hydrograph separation at high temporal resolution for both the rising and falling limb of the hydrograph. The following results of the investigations can be resumed. If relief conditions are similar, the effectiveness of runoff production decreases with an increasing density of vegetation cover. The runoff delivery ratio decreases as well as the peaks of runoff. In contrast, concentration times of hillslope catchments are equal, even if vegetation cover is of great density and soils are well developed. As a reason for the short reaction times, different runoff production processes have been detected. On bare ground, infiltration excess overland flow intensified by surface sealing processes is the main source for quick runoff. On hillslopes well covered by vegetation, translatory flow processes indicated by soil water with high solute contents force a rapid runoff reaction only a few minutes after rainfall has begun. It is to be assumed that translatory flow is a runoff production process typical for hillslopes covered by vegetation in a steep alpine relief. By means of the areal distribution of the topographic index, concentration of runoff production on a small part of the catchment has been demonstrated for hillslopes densely covered by vegetation. The investigations have shown that there is a lack of studies on runoff production processes in steep alpine relief, as well as a deficit of methods to quantify hydraulic properties of coarse‐grained soils with a wide grain size distribution. Copyright © 2003 John Wiley & Sons, Ltd.     

Solifluction meets vegetation: the role of biogeomorphic feedbacks for turf‐banked solifluction lobe development

Vegetation is an important factor influencing solifluction processes, while at the same time, solifluction processes and landforms influence species composition, fine‐scale distribution and corresponding ecosystem functioning. However, how feedbacks between plants and solifluction processes influence the development of turf‐banked solifluction lobes (TBLs) and their geomorphic and vegetation patterns is still poorly understood. We addressed this knowledge gap in a detailed biogeomorphic investigation in the Turtmann glacier foreland (Switzerland). Methods employed include geomorphic and vegetation mapping, terrain assessment with unmanned aerial vehicle (UAV) and temperature logging. Results were subsequently integrated with knowledge from previous geomorphic and ecologic studies into a conceptual model. Our results show that geomorphic and vegetation patterns at TBLs are closely linked through the lobe elements tread, risers and ridge. A conceptual four‐stage biogeomorphic model of TBL development with ecosystem engineering by the dwarf shrub Dryas octopetala as the dominant process can explain these interlinked patterns. Based on this model, we demonstrate that TBLs are biogeomorphic structures and follow a cyclic development, during which the role of their components for engineer and non‐engineer species changes. Our study presents the first biogeomorphic model of TBL development and highlights the applicability and necessity of biogeomorphic approaches and research in periglacial environments. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of leaf and canopy characteristics on rainfall interception and urban hydrology

Considering the rapid expansion of urban populations and the corresponding urbanization of landscapes, a dearth of knowledge exists regarding the role of urban vegetation in modulating urban ecosystem functioning. In response to the need for the development of new approaches to quantify ecohydrological processes along urban-to-rural gradients at alternate scales, this study explores the relationship between individual plant selection choices in landscaping and changes in urban hydrological functioning. This research examines differences in the variation of rainfall interception, leaf hydrophobicity, canopy structure, and water storage, between 13 species in an urban, semi-arid location. The species studied were selected based on resident preferences, and hence this research considers the role that urban residents play, through individual choices, in modifying the ecohydrology of an urban watershed. Rainfall interception, canopy surface storage, leaf hydrophobicity, and water droplet retention were significantly different between species. Results indicate that individual choice in plant selection for landscaping may influence urban hydrology.     

Dominant processes concept,model simplification and classification framework in catchment hydrology

Technological and methodological advances have facilitated tremendous growth in hydrology during the last century; however, there are also concerns that these advances indirectly contribute to additional problems in our research. An insight into hydrologic literature reveals our tendency to develop more complex models than perhaps needed, and our increasing emphasis on individual mathematical techniques rather than general hydrologic issues. Some recent studies of diverse forms have suggested that simplification in modeling and development of a common framework may help alleviate these problems. The present study is intended to bring such studies together towards a more coherent approach to research in catchment hydrology. This is done by highlighting the need for model simplification and generalization and proposing some potential directions for achieving such. Through a discussion of difficulties in data measurements, the need for moving beyond the notion of “modeling everything” to the notion of “capturing the essential features” is explained; the concept of dominant processes in model simplification and the utility of integration of concepts for modeling improvement are discussed. Formulation of a catchment classification framework is advocated as a possible means for a common framework in hydrology, and the role of dominant processes in this formulation is presented; the problems due to adoption of different modeling terminologies are highlighted and potential ways to overcome such are also discussed.     

CORRELATIONS BETWEEN DEM-DERIVED TOPOGRAPHIC INDICES AND REMOTELY-SENSED VEGETATION COVER IN RANGELANDS

The paper describes an attempt to relate patterns of vegetation cover with topography and a set of biological and grazing intensity variables in a mountain and piedmont area of arid central Australia. Vegetation cover, as measured by an index based on data from the Landsat satellite, can also be used as an erosion/deposition surrogate so the results have implications for distributed erosion models. A simple, analytically based erosion model derived from the continuity equation does not reproduce observed patterns of vegetation cover, and neither do various topographically based moisture indices. A regression approach shows that patterns of vegetation cover are related to topography but the most important predictors are biological ones, with percentage of bare ground upslope being the strongest. Tests with variable drainage area show that relationships between cover and topography, bare area upslope and grazing effects change systematically with basin size and that scale effects are present. Distributed erosion models are not yet capable of handling biological processes very well, yet these processes must be incorporated if erosion prediction is to be successful.     

The influence of riparian vegetation on near‐bank turbulence: a flume experiment

Measurements from a fixed‐bed, Froude‐scaled hydraulic model of a stream in northeastern Vermont demonstrate the importance of forested riparian vegetation effects on near‐bank turbulence during overbank flows. Sections of the prototype stream, a tributary to Sleepers River, have increased in channel width within the last 40 years in response to passive reforestation of its riparian zone. Previous research found that reaches of small streams with forested riparian zones are commonly wider than adjacent reaches with non‐forested, or grassy, vegetation; however, driving mechanisms for this morphologic difference are not fully explained. Flume experiments were performed with a 1:5 scale, simplified model of half a channel and its floodplain, mimicking the typical non‐forested channel size. Two types of riparian vegetation were placed on the constructed floodplain: non‐forested, with synthetic grass carpeting; and forested, where rigid, randomly distributed, wooden dowels were added. Three‐dimensional velocities were measured with an acoustic Doppler velocimeter at 41 locations within the channel and floodplain at near‐bed and 0·6‐depth elevations. Observations of velocity components and calculations of turbulent kinetic energy (TKE), Reynolds shear stress and boundary shear stress showed significant differences between forested and non‐forested runs. Generally, forested runs exhibited a narrow band of high turbulence between the floodplain and main channel, where TKE was roughly two times greater than TKE in non‐forested runs. Compared to non‐forested runs, the hydraulic characteristics of forested runs appear to create an environment with higher erosion potential. Given that sediment entrainment and transport can be amplified in flows with high turbulence intensity and given that mature forested stream reaches are wider than comparable non‐forested reaches, our results demonstrated a possible driving mechanism for channel widening during overbank flow events in stream reaches with recently reforested riparian zones. Copyright © 2007 John Wiley & Sons, Ltd.     

Towards estimation of extreme floods: examination of the roles of runoff process changes and floodplain flows

This paper presents the development and application of a distributed rainfall-runoff model for extreme flood estimation, and its use to investigate potential changes in runoff processes, including changes to the ‘rating curve’ due to effects of over-bank flows, during the transition from ‘normal’ floods to ‘extreme’ floods. The model has two components: a hillslope runoff generation model based on a configuration of soil moisture stores in parallel and series, and a distributed flood routing model based on non-linear storage-discharge relationships for individual river reaches that includes the effects of floodplain geometries and roughnesses. The hillslope water balance model contains a number of parameters, which are measured or derived a priori from climate, soil and vegetation data or streamflow recession analyses. For reliable estimation of extreme discharges that may extend beyond recorded data, the parameters of the flood routing model are estimated from hydraulic properties, topographic data and vegetation cover of compound channels (main channel and floodplains). This includes the effects of the interactions between the main channel and floodplain sections, which tend to cause a change to the rating curve. The model is applied to the Collie River Basin, 2545 km², in Western Australia and used to estimate the probable maximum flood (PMF) from probable maximum precipitation estimates for this region. When moving from normal floods to the PMFs, application of the model demonstrates that the runoff generation process changes with a substantial increase of saturation excess overland flow through the expansion of saturated areas, and the dominant runoff process in the stream channel changes from in-bank to over-bank flows. The effects of floodplain inundation and floodplain vegetation can significantly reduce the magnitude of the estimated PMFs. This study has highlighted the need for the estimation of a number of critical parameters (e.g. cross-sectional geometry, floodplain vegetation, soil depths) through concerted field measurements or surveys, and targeted laboratory experiments.     

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.     

Investigation of spatially distributed braided river flows using a two‐dimensional hydraulic model

This paper reports the application of a two‐dimensional hydraulic model to a braided reach of the Avoca River, New Zealand. Field measurements of water surface elevation, depth and velocity obtained at low flow were used to validate the model and to optimize the parameterization of bed friction. The main systematic trends in the measured flow variables are reproduced by the model. However, field data are characterized by greater spatial variability than model output reflecting differences in the scale of processes measured in the field and represented by the model. Additional model runs were conducted to simulate flow patterns within the study reach at five higher discharges. The purpose of these simulations was to evaluate the potential for using two‐dimensional hydraulic models to quantify the reach‐scale hydraulic characteristics of braided rivers and their dependence on discharge. Changes in flow depth and velocity with increasing discharge exhibit trends that are consistent with the results of previous field investigations, although the tendency for the wetted area of the braidplain within particular depth and velocity categories to remain fixed as discharge rises, as has been noted for several braided rivers in New Zealand, was not observed. Modelled shear stress frequency distributions fit gamma functions that incorporate a distribution shape parameter, the value of which follows clear systematic trends with rising discharge. These results illustrate both the problems of, and potential for, using two‐dimensional hydraulic models in braided river applications. This leads to something of a paradox in that while such models provide a means of generating hydraulic information that would be difficult to obtain in the field at an equivalent spatial resolution, they are, due to the problems inherent to data collection, difficult to validate conclusively. Despite this limitation, the application of spatially distributed models to investigate relationships between discharge and reach‐scale form and process variables appears to have considerable potential. Copyright © 2003 John Wiley & Sons, Ltd.     

Swash-zone morphodynamics

Hydrodynamic forcing with respect to sediment transport and morphologic change, paying particular attention to relevant swash asymmetries, is reviewed. The hydrodynamics are categorized into their individual effects: high- and low-frequency motions, bores and turbulence, in/exfiltration, shear stresses and friction coefficients. Individual effects are then related to their potential for driving or influencing sediment transport and morphological change. Additional concepts such as settling/scour lag and sediment advection that have largely been ignored are also discussed. A simple framework is presented for the morphological response of the beachface under swash zone hydrodynamic processes. The framework acknowledges that the beachface cannot be considered in isolation from the surf zone and that the two zones are strongly linked through feedback processes. The swash zone itself is also a morphodynamic system and morphological response occurs as a result of disequilibrium between the beachface gradient and asymmetries in the swash hydrodynamics. Any beachface morphological development in response to such equilibrium will have direct and indirect effects on swash hydrodynamics and sediment transport processes. It is concluded that the two issues requiring most urgent research attention with regards to swash zone sediment transport processes are the roles of sediment advection and longshore swash motion.     

Formation of asymmetrical loess gullies in the northeastern loess plateau of China

Loess gullies are the most active and changeable landform unit on the Loess Plateau of China. Under the influence of inhomogeneous internal and external forces, various gully morphologies have been identified as specific forms of asymmetrical loess gullies in the northeastern Loess Plateau. Thus, the formation mechanisms of asymmetrical gullies should be examined to better understand the gully evolution processes in this area. In this study, a typical asymmetrical gully area and its geological background in the northeastern Loess Plateau are investigated. Then, the asymmetrical gullies are extracted and ordered under different watershed hierarchies using 5 m horizontal resolution digital elevation models. The asymmetrical gullies are characterized using the gully deviation index and gully asymmetrical coefficient to quantitatively and qualitatively describe the gully formation from the perspective of gully morphology. Subsequently, environmental factors, such as the bedrock, climate, vegetation and interactions with neighbouring watersheds, are combined to achieve an in-depth understanding of the mechanisms of asymmetrical gully formation. The results show that most watersheds shift to the right side of the watershed geometric centre line, thereby forming a specific asymmetrical gully morphology. The phenomenon in which the asymmetrical degree characteristics decrease with the increase in drainage area suggests evident morphological differences on both sides of the main channel on a small scale, and relatively weak morphological differences on both sides of the main channel on a large scale. The degree of loess gully asymmetry appears higher in the area where only the windward slope is covered by loess than in areas where all slopes are covered by loess. The interaction between adjacent watersheds also influences the formation of asymmetrical gullies. These results support the understanding of asymmetrical gully formation in relation to the underlying bedrock structure and gully reorganization, thereby contributing to the development of process-based gully evolution models.     

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.