A modelling framework for evaluation of the hydrological impacts of nature‐based approaches to flood risk management,with application to in‐channel interventions across a 29‐km2 scale catchment in the United Kingdom

Nature‐based approaches to flood risk management are increasing in popularity. Evidence for the effectiveness at the catchment scale of such spatially distributed upstream measures is inconclusive. However, it also remains an open question whether, under certain conditions, the individual impacts of a collection of flood mitigation interventions could combine to produce a detrimental effect on runoff response. A modelling framework is presented for evaluation of the impacts of hillslope and in‐channel natural flood management interventions. It couples an existing semidistributed hydrological model with a new, spatially explicit, hydraulic channel network routing model. The model is applied to assess a potential flood mitigation scheme in an agricultural catchment in North Yorkshire, United Kingdom, comprising various configurations of a single variety of in‐channel feature. The hydrological model is used to generate subsurface and surface fluxes for a flood event in 2012. The network routing model is then applied to evaluate the response to the addition of up to 59 features. Additional channel and floodplain storage of approximately 70,000 m³ is seen with a reduction of around 11% in peak discharge. Although this might be sufficient to reduce flooding in moderate events, it is inadequate to prevent flooding in the double‐peaked storm of the magnitude that caused damage within the catchment in 2012. Some strategies using features specific to this catchment are suggested in order to improve the attenuation that could be achieved by applying a nature‐based approach.     

Hydraulic and geomorphic processes in an overbank flood along a meandering,gravel‐bed river: implications for chute formation

Hydraulic interactions between rivers and floodplains produce off‐channel chutes, the presence of which influences the routing of water and sediment and thus the planform evolution of meandering rivers. Detailed studies of the hydrologic exchanges between channels and floodplains are usually conducted in laboratory facilities, and studies documenting chute development are generally limited to qualitative observations. In this study, we use a reconstructed, gravel‐bedded, meandering river as a field laboratory for studying these mechanisms at a realistic scale. Using an integrated field and modeling approach, we quantified the flow exchanges between the river channel and its floodplain during an overbank flood, and identified locations where flow had the capacity to erode floodplain chutes. Hydraulic measurements and modeling indicated high rates of flow exchange between the channel and floodplain, with flow rapidly decelerating as water was decanted from the channel onto the floodplain due to the frictional drag provided by substrate and vegetation. Peak shear stresses were greatest downstream of the maxima in bend curvature, along the concave bank, where terrestrial LiDAR scans indicate initial floodplain chute formation. A second chute has developed across the convex bank of a meander bend, in a location where sediment accretion, point bar development and plant colonization have created divergent flow paths between the main channel and floodplain. In both cases, the off‐channel chutes are evolving slowly during infrequent floods due to the coarse nature of the floodplain, though rapid chute formation would be more likely in finer‐grained floodplains. The controls on chute formation at these locations include the flood magnitude, river curvature, floodplain gradient, erodibility of the floodplain sediment, and the flow resistance provided by riparian vegetation. Copyright © 2015 John Wiley & Sons, Ltd.     

Implications of heterogeneous flood‐frequency distributions on traditional stream‐discharge prediction techniques

Traditional flood‐frequency analysis involves the assumption of homogeneity of the flood distribution. However, floods are often generated by heterogeneous distributions composed of a mixture of two or more populations. Differences between the populations may be the result of a number of factors, including seasonal variations in the flood‐producing mechanisms, changes in weather patterns resulting from low‐frequency climate shifts and/or El Niño/La Nina oscillations, changes in channel routing owing to the dominance of within‐channel or floodplain flow, and basin variability resulting from changes in antecedent soil moisture. Not recognizing these physical processes in conventional flood‐frequency analysis probably is the main reason that many frequency distributions do not provide an acceptable fit to flood data. In this paper, we use long‐term hydroclimatic records from the Gila River basin of south‐east and central Arizona in the USA to explore the extent and significance of mixed populations. First, we discuss the probable causes of heterogeneity in the frequency distribution of annual flood and present evidence of its occurrence. Second, we investigate the implications of using various popular homogeneous distributions for predicting peak flows for basins that exhibit mixed population characteristics. Third, we demonstrate how alternative frequency models that explicitly account for floods generated by a mixture of two or more populations are both hydrologically and statistically more appropriate. We illustrate how the selection of the most plausible distribution for flood‐frequency analysis also should be based on hydrological reasoning as opposed to the sole application of the traditional statistical goodness‐of‐fit tests. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Modelling stream flow for use in ecological studies in a large,arid zone river,central Australia

Australian arid zone ephemeral rivers are typically unregulated and maintain a high level of biodiversity and ecological health. Understanding the ecosystem functions of these rivers requires an understanding of their hydrology. These rivers are typified by highly variable hydrological regimes and a paucity, often a complete absence, of hydrological data to describe these flow regimes. A daily time‐step, grid‐based, conceptual rainfall–runoff model was developed for the previously uninstrumented Neales River in the arid zone of northern South Australia. Hourly, logged stage data provided a record of stream‐flow events in the river system. In conjunction with opportunistic gaugings of stream‐flow events, these data were used in the calibration of the model. The poorly constrained spatial variability of rainfall distribution and catchment characteristics (e.g. storage depths) limited the accuracy of the model in replicating the absolute magnitudes and volumes of stream‐flow events. In particular, small but ecologically important flow events were poorly modelled. Model performance was improved by the application of catchment‐wide processes replicating quick runoff from high intensity rainfall and improving the area inundated versus discharge relationship in the channel sections of the model. Representing areas of high and low soil moisture storage depths in the hillslope areas of the catchment also improved the model performance. The need for some explicit representation of the spatial variability of catchment characteristics (e.g. channel/floodplain, low storage hillslope and high storage hillslope) to effectively model the range of stream‐flow events makes the development of relatively complex rainfall–runoff models necessary for multisite ecological studies in large, ungauged arid zone catchments. Grid‐based conceptual models provide a good balance between providing the capacity to easily define land types with differing rainfall–runoff responses, flexibility in defining data output points and a parsimonious water‐balance–routing model. Copyright © 2004 John Wiley & Sons, Ltd.     

Sediment routing analyses based on chronological changes in hillslope and riverbed morphologies

Recent studies of sediment delivery and budgets in the United States indicate that upland erosion rates at a given time may not explain contemporaneous sediment yields from a drainage basin. This suggests temporal discontinuities in sediment delivery associated with hillslope and channel storage processes. Integration of sediment production, storage and transport is essential to understand sediment routing in basins. We analysed each process chronologically using aerial photographs, monitoring data of sediment movement and annual tree-rings, and then compared estimated temporal changes in sediment production from hillslopes, floodplain disturbance areas and sediment transport in river channels. Toeslopes, floodplains and alluvial fans together contained 59 per cent of sediment eroded from uplands over the last 30 years. Monitoring results of riverbed changes showed that the volume of stored sediment on floodplains decreased exponentially with succeeding floods. The age distribution of floodplain deposits reflected the disturbance history of a river channel, and followed an exponential decrease with age. The results of this study may have important implications for sediment control plans for watersheds in steep regions.     

Simulating hydrologic and hydraulic processes throughout the Amazon River Basin

Presented here is a model framework based on a land surface topography that can be represented with various degrees of resolution and capable of providing representative channel/floodplain hydraulic characteristics on a daily to hourly scale. The framework integrates two models: (1) a water balance model (WBM) for the vertical fluxes and stores of water in and through the canopy and soil layers based on the conservation of mass and energy, and (2) a routing model for the horizontal routing of surface and subsurface runoff and channel and floodplain waters based on kinematic and diffusion wave methodologies. The WBM is driven by satellite‐derived precipitation (TRMM_3B42) and air temperature (MOD08_M3). The model's use of an irregular computational grid is intended to facilitate parallel processing for applications to continental and global scales. Results are presented for the Amazon Basin over the period Jan 2001 through Dec 2005. The model is shown to capture annual runoff totals, annual peaks, seasonal patterns, and daily fluctuations over a range of spatial scales (>1, 000 to < 4·7M km²). For the period of study, results suggest basin‐wide total water storage changes in the Amazon vary by approximately + /? 5 to 10 cm, and the fractional components accounting for these changes are: root zone soil moisture (20%), subsurface water being routed laterally to channels (40%) and channel/floodplain discharge (40%). Annual variability in monthly water storage changes by + /? 2·5 cm is likely due to 0·5 to 1 month variability in the arrival of significant rainfall periods throughout the basin. Copyright © 2009 John Wiley & Sons, Ltd.     

A hydro‐economic modelling framework for flood damage estimation and the role of riparian vegetation

A modelling framework for the quick estimate of flood inundation and the resultant damages is developed in this paper. The model, called the flood economic impact analysis system (FEIAS), can be applied to a river reach of any hydrogeological river basin. For the development of the integrated modelling framework, three models were employed: (1) a modelling scheme based on the Hydrological Simulation Program FORTRAN model that was developed for any geomorphological river basin, (2) a river flow/floodplain model, and (3) a flood loss estimation model. The first sub‐model of the flood economic impact analysis system simulates the hydrological processes for extended periods of time, and its output is used as input to a second component, the river/floodplain model. The hydraulic model MIKE 11 (quasi‐2D) is the river/floodplain model employed in this study. The simulated flood parameters from the hydraulic model MIKE 11 (quasi‐2D) are passed, at the end of each time step, to a third component, the flood loss model for the estimation of flood damage. In the present work, emphasis was given to the seasonal variation of Manning's coefficient (n), which is an important parameter for the determination of the flood inundation in hydraulic modelling. High values of Manning's coefficient for a channel indicate high flow resistance. The riparian vegetation can have a large impact on channel resistance. The modelling framework developed in this paper was used to investigate the role of riparian vegetation in reducing flood damage. Moreover, it was used to investigate the influence of cutting riparian vegetation scenarios on the flow characteristics. The proposed framework was applied to the downstream part of the Koiliaris River basin in Crete, Greece, and was tested and validated with historical data. Copyright © 2012 John Wiley & Sons, Ltd.     

Morphology of Red Creek,Wyoming, an arid-region anastomosing channel system

Red Creek, in the Red Desert area of the Great Divide Basin, Wyoming, is an arid-region anastomosing stream. The narrow, deep, and sinuous main channel is flanked by anastomosing flood channels, or anabranches. Most anabranches are initiated at meander bends. The primary mechanism of anabranch initiation is avulsion during overbank floods. Anabranch enlargement occurs by headward erosion. Anabranches act as distributary channels during floods, when water and sediment from overbank flows are transported to and deposited on the floodplain via the anabranches. During periods of low discharges, the anabranches act as tributaries to the main channel, transporting runoff from the floodplain and surrounding hillslopes to the main channel of Red Creek. Aggradation is occurring in the main channel and on the floodplain throughout the study reach. Infilling of the main channel occurs primarily by lateral accretion, while the floodplain accretes vertically through deposition of overbank sediment from the main channel and anabranches. Infilling of the main channel may cause avulsion of the main channel into an anabranch. The abandoned main channel segment may then fill completely or act as an anabranch. Because lateral migration of channels is inhibited by the high cohesion of the silt and clay channel sediment, periodic avulsion is the primary form of lateral mobility in the system.     

Ecological Conditions of Watercourses in the Middle Amur Lowland in the Areas of Drainage Reclamation

The effect of land reclamation on the migration and accumulation of heavy metals and some pollutants in small rivers is determined. The study was conducted in the Middle Amur lowland in different phases of hydrological regime and at different extent of floodplain inundation. The formation of the hydrological regime in small rivers in this area is governed by the irregular annual runoff distribution. The seasonal character of flood periods requires water samples to be taken before spring flood (April) and after floods have passed (September–October), and at various extent of floodplain inundation. The field studies of water-courses were carried out in 2009–2014. The water samples were analyzed to determine the concentrations of heavy metals (iron, manganese, copper, nickel, cobalt, lead and zinc), suspended substances, organic carbon (total, dissolved, and suspended), humic and fulvic acids, and volatile organic compounds. The mobility of heavy metals in surface watercourses was shown to reflect the effect of drainage and surface runoff from soil horizons, an increase in the concentrations of suspended and organic compounds, especially, fulvic acids, which enhance their geochemical mobility. During floods, an increase in heavy metals washout from floodplain soils and the dilution of their concentrations causes the formation of a single-type concentration series of heavy metals. In addition, the processes of pollutants migration show an effect of changes in the geomorphological characteristics of floodplain–channel complexes and a decrease in flow velocity in watercourses in the areas where land reclamation was applied. Thus, all factors mentioned above lead to a decrease in water pollution index in the river.     

Including prior information in the estimation of effective soil parameters in unsaturated zone modelling

In this paper we propose a methodology to include prior information in the estimation of effective soil parameters for modelling the soil moisture content in the unsaturated zone. Laboratory measurements on undisturbed soil cores were used to estimate the moisture retention curve and hydraulic conductivity curve parameters. The soil moisture content was measured at 25 locations along three transects and at three different depths (surface, 30 and 60 cm) on an 80×20 m hillslope for the year 2001. Soil cores were collected in 84 locations situated in three profile pits along the hillslope. For the estimation of the effective soil hydraulic parameters the joint probability distribution of measured parameter values was used as prior information. A two-horizon single column 1D MIKE SHE model based on Richards' equation was set-up for nine soil moisture measurement locations along the middle transect of the hillslope. The goal of the model is to simulate the soil moisture profile at each location. The shuffled complex evolution (SCE) algorithm has been applied to estimate effective model parameters using either wide parameter ranges, referred to as the ‘no-prior’ case, or the joint probability distribution of measured parameter values as prior information (‘prior’ case). When the prior information is incorporated in the SCE optimisation the goodness-of-fit of the model predictions is only slightly worse compared to when no-prior information is incorporated. However, the effective parameter estimates are more realistic when the prior information is incorporated. For both the no-prior and prior case the generalised likelihood uncertainty estimation procedure (GLUE) was subsequently used to estimate the uncertainty bounds (UB) on the model predictions. When incorporating the prior information more parameter sets were accepted for the estimation of the predictive uncertainty and the parameter values were more realistic. Moreover, UB better enclosed the observations. Thus, incorporating prior information in GLUE reduces the amount of model evaluations needed to obtain sufficient behavioural parameter sets. The results indicate the importance of prior information in the SCE and GLUE parameter estimation strategies.     

Modelling hyperconcentrated flow in the Yellow River

A large amount of the total sediment load in the Chinese Yellow River is transported during hyperconcentrated floods. These floods are characterized by very high suspended sediment concentrations and rapid morphological changes with alternating sedimentation and erosion in the main channel, and persistent sedimentation on the floodplain. However, the physical mechanisms driving these hyperconcentrated floods are still poorly understood. Numerical modelling experiments of these floods reveal that sedimentation is largely caused by large vertical concentration gradients, both in the channel during the rising stage of the flood, as well as on the floodplains, during a later stage of the flood. These vertical concentration gradients are large because the turbulent mixing rates are reduced by the increased sediment‐induced density gradients, resulting in a positive feedback mechanism that produces high deposition rates. Erosion prevails when the sediment is largely held in suspension due to hindered settling, and is strengthened by the reduced wetted cross‐section caused by massive sedimentation on the floodplain. Observed patterns of erosion and sedimentation during these floods can be qualitatively reproduced with a numerical model in which sediment‐induced density effects and hindered settling are included. Copyright © 2009 John Wiley & Sons, Ltd.     

An integral-balance nonlinear model to simulate changes in soil moisture,groundwater and surface runoff dynamics at the hillslope scale

We present a system of ordinary differential equations (ODEs) capable of reproducing simultaneously the aggregated behavior of changes in water storage in the hillslope surface, the unsaturated and the saturated soil layers and the channel that drains the hillslope. The system of equations can be viewed as a two-state integral-balance model for soil moisture and groundwater dynamics. Development of the model was motivated by the need for landscape representation through hillslopes and channels organized following stream drainage network topology. Such a representation, with the basic discretization unit of a hillslope, allows ODEs-based simulation of the water transport in a basin. This, in turn, admits the use of highly efficient numerical solvers that enable space–time scaling studies. The goal of this paper is to investigate whether a nonlinear ODE system can effectively replicate observations of water storage in the unsaturated and saturated layers of the soil. Our first finding is that a previously proposed ODE hillslope model, based on readily available data, is capable of reproducing streamflow fluctuations but fails to reproduce the interactions between the surface and subsurface components at the hillslope scale. However, the more complex ODE model that we present in this paper achieves this goal. In our model, fluxes in the soil are described using a Taylor expansion of the underlying storage flux relationship. We tested the model using data collected in the Shale Hills watershed, a 7.9-ha forested site in central Pennsylvania, during an artificial drainage experiment in August 1974 where soil moisture in the unsaturated zone, groundwater dynamics and surface runoff were monitored. The ODE model can be used as an alternative to spatially explicit hillslope models, based on systems of partial differential equations, which require more computational power to resolve fluxes at the hillslope scale. Therefore, it is appropriate to be coupled to runoff routing models to investigate the effect of runoff and its uncertainty propagation across scales. However, this improved performance comes at the expense of introducing two additional parameters that have no obvious physical interpretation. We discuss the implications of this for hydrologic studies across scales.     

The introduction of runoff routing into large-scale hydrological models

Two alternative schemes are presented that are appropriate for the representation of runoff routing in large-scale grid-based hydrological models and atmospheric general circulation models (AGCMs). The first scheme characterizes routing processes as a single conceptual store. The second scheme, developed by Naden (1992), uses the normalized network width function to characterize the channel network form and a linear solution to the convective diffusion equation of one-dimensional flow to characterize the routing effect of a single channel. Both schemes are applied to the Severn catchment at the daily time-scale for the period 1981 to 1990 using a grid resolution of 40 km. Comparable results were obtained using both schemes (efficiencies were of the order of 80% in both cases). A combined model using a conceptual reservoir to represent hillslope routing and the network-based scheme to represent channel routing was developed to investigate the relative roles of hillslope and channel routing at the catchment scale. The application of this model demonstrated the important role of hillslope routing in reproducing the low frequency component of the catchment response. However, in terms of goodness-of-fit there was little to choose between the three schemes. Consequently, it is recommended that additional a priori knowledge of the routing processes should be used to condition the choice of model structure. © 1997 John Wiley & Sons, Ltd.     

Temporal variability of contemporary floodplain sedimentation in the Rhine–Meuse delta,The Netherlands

It is often believed that extreme but infrequent events are most important in the development of landforms. When evaluating the overall effect of large floods on floodplain sedimentation, quantitative measurements of both high- and low-magnitude events should be considered. To analyse the role of flood magnitude on floodplain sedimentation, we measured overbank sedimentation during floods of different magnitude and duration. The measurements were carried out on two embanked floodplain sections along the rivers Rhine and Meuse in The Netherlands, using sediment traps made of artificial grass. The results showed an increase in total sediment accumulation with flood magnitude, mainly caused by enhanced accumulation of sand. At low floodplain sections the increase in sediment deposition was smaller than expected from the strong increase in suspended sediment transport in the river. Spatial variability in sediment accumulation was found to depend both on flood magnitude and duration. Deposition of sand on natural levees mainly takes place during high-magnitude floods, whilst low floods and slowly receding floods are important for the deposition of silt and clay in low-lying areas, at greater distance from the main channel. © 1998 John Wiley & Sons, Ltd.     

Floodplain inundation: Processes and relationships with channel discharge

The factors affecting the relationship between channel discharge and volume of inundation are discussed. For many floodplains this relationship is not simple, but involves hysteretic effects which vary according to the hydrological characteristics of individual floods and the way in which these interact with the surface form of the floodplain reach. Some methods of deriving or estimating the extent of the hysteresis on floodplain reaches are suggested, although an acute lack of available data prevents detailed examination at present. Further investigations are required before floodplain conductivity relations can be used as an aid to flood routing procedures or in floodplain management problems. A major priority must be the acquisition of sequential inundation data.     

Quantitative model of the growth of floodplains by vertical accretion

A simple one‐dimensional model is developed to quantitatively predict the change in elevation, over a period of decades, for vertically accreting floodplains. This unsteady model approximates the monotonic growth of a floodplain as an incremental but constant increase of net sediment deposition per flood for those floods of a partial duration series that exceed a threshold discharge corresponding to the elevation of the floodplain. Sediment deposition from each flood increases the elevation of the floodplain and consequently the magnitude of the threshold discharge resulting in a decrease in the number of floods and growth rate of the floodplain. Floodplain growth curves predicted by this model are compared to empirical growth curves based on dendrochronology and to direct field measurements at five floodplain sites. The model was used to predict the value of net sediment deposition per flood which best fits (in a least squares sense) the empirical and field measurements; these values fall within the range of independent estimates of the net sediment deposition per flood based on empirical equations. These empirical equations permit the application of the model to estimate of floodplain growth for other floodplains throughout the world which do not have detailed data of sediment deposition during individual floods. Copyright © 2000 John Wiley & Sons, Ltd.     

Influence of rainfall and soil properties spatial aggregation on extreme flash flood response modelling: An evaluation based on the Sesia river basin,North Western Italy

High resolution radar rainfall fields and a distributed hydrologic model are used to evaluate the sensitivity of flood and flash flood simulations to spatial aggregation of rainfall and soil properties at catchment scales ranging from 75 to 983 km². Hydrologic modeling is based on a Hortonian infiltration model and a network-based representation of hillslope and channel flow. The investigation focuses on three extreme flood and flash flood events occurred on the Sesia river basin, North Western Italy, which are analysed by using four aggregation lengths ranging from 1 to 16 km. The influence of rainfall spatial aggregation is examined by using the flow distance as a spatial coordinate, hence emphasising the role of river network in the averaging of space–time rainfall. The effects of reduced and distorted rainfall spatial variability on peak discharge have been found particularly severe for the flash flood events, with peak errors up to 35% for rainfall aggregation of 16 km and at 983 km² catchment size. Effects are particularly remarkable when significant structured rainfall variability combines with relatively important infiltration volumes due to dry initial conditions, as this emphasises the non-linear character of the rainfall–runoff relationship. In general, these results confirm that the correct estimate of rainfall volume is not enough for the accurate reproduction of flash flood events characterised by large and structured rainfall spatial variability, even at catchment scales around 250 km². However, accurate rainfall volume estimation may suffice for less spatially variable flood events. Increasing the soil properties aggregation length exerts similar effects on peak discharge errors as increasing the rainfall aggregation length, for the cases considered here and after rescaling to preserve the rainfall volume. Moreover, peak discharge errors are roughly proportional to runoff volume errors, which indicates that the shape of the flood wave is influenced in a limited way by modifying the detail of the soil property spatial representation. Conversely, rainfall aggregation may exert a pronounced influence on the discharge peak by reshaping the spatial organisation of the runoff volumes and without a comparable impact on the runoff volumes.     

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.     

Precipitation‐runoff processes in Shimen hillslope micro‐catchment of Taihang Mountain,north China

Precipitation runoff is a critical hillslope hydrological process for downslope streamflow and piedmont/floodplain recharge. Shimen hillslope micro‐catchment is strategically located in the central foothill region of Taihang Mountains, where runoff is crucial for water availability in the piedmont corridors and floodplains of north China. This study analyzes precipitation‐runoff processes in the Shimen hillslope micro‐catchment for 2006–2008 using locally designed runoff collection systems. The study shows that slope length is a critical factor, next only to precipitation, in terms of runoff yield. Regression analysis also shows that runoff is related positively to precipitation, and negatively to slope length. Soil mantle in the study area is generally thin and is therefore not as critical a runoff factor as slope length. The study shows a significant difference between overland and subsurface runoff. However, that between the 0–10 and 10–20 cm subsurfaces is insignificant. Runoff hardly occurs under light rains (<10 mm), but is clearly noticeable under moderate‐to‐rainstorm events. In the hillslope catchment, vertical infiltration (accounting for 42–84% of the precipitation) dominates runoff processes in subsurface soils and weathered granite gneiss bedrock. A weak lateral flow (at even the soil/bedrock interface) and the generally small runoff suggest strong infiltration loss via deep percolation. This is critical for groundwater recharge in the downslope piedmont corridors and floodplains. This may enhance water availability, ease water shortage, avert further environmental degradation, and reduce the risk of drought/flood in the event of extreme weather conditions in the catchment and the wider north China Plain. Copyright © 2011 John Wiley & Sons, Ltd.