Rainfall,run‐off,and sediment transport dynamics in a humid mountain badland area: Long‐term results from a small catchment

In this study, we investigated rainfall, run‐off, and sediment transport dynamics (414 run‐off events and 231 events with sediment information) of a humid mountain badland area—the Araguás catchment (Central Pyrenees, Spain)—from October 2005 to September 2016. Use of this long‐term database allows characterization of the hydrological response, which consist of low‐magnitude/high‐frequency events and high‐magnitude/low‐frequency events, and identification of seasonal dynamics and rainfall‐run‐off thresholds. Our results indicate that the Araguás catchment, similarly to other humid badlands, had high hydrological responsiveness (mean annual run‐off coefficient: 0.52), a non‐linear relationship of rainfall with run‐off (common in Mediterranean environments), and seasonal hydrological and sedimentological dynamics. We created and validated a multivariate regression model to characterize the hydrological variables (stormflow and peak discharge) and sedimentological variables (mean and maximum suspended sediment concentrations and total suspended sediment load). In summer and at the beginning of autumn, the response was mainly related to rainfall intensity, suggesting a predomination of Hortonian flows. In contrast, in spring and winter, the responses were mainly related to the antecedent conditions (previous rainfall and baseflow), suggesting the occurrence of saturated excess flow processes, and the contribution of neighbouring vegetated areas. The multivariate analysis also showed that total sediment load is better predicted by a multivariate regression model that integrates pre‐event, rainfall, and run‐off variables. In general, our models provided more accurate predictions of small‐magnitude/high‐frequency events than high‐magnitude/low‐frequency events. This study highlights the high inter‐ and intra‐annual variability response in humid badland areas and that long‐term records are needed to reduce the uncertainty of hydrological and sedimentological responses in Mediterranean badland areas.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Effects of rock fragment cover on hydrological processes under rainfall simulation in a semi‐arid region of China

Rock fragment cover has long been an important agricultural crop production technique on the Loess Plateau, China. Although this approach plays an important role in controlling hydrological processes and preventing soil erosion, inconsistent results have been recovered in this field. In this study, we investigated the effects of rock fragment cover on infiltration, run‐off, soil erosion, and hydraulic parameters using rainfall simulation in the field in a semi‐arid region of China. Two field plots encompassing 6 rock fragment coverages (0%, 10%, 20%, 25%, 30%, and 40%), as well as 2 rock fragment positions and sizes were exposed to rainfall at a particular intensity (60 mm h^?1). The results of this study showed that increasing the rock fragment coverage with rock fragments resting on the soil surface increased infiltration but decreased run‐off generation and sediment yield. A contrasting result was found, however, when rock fragments were partially embedded into the soil surface; in this case, a positive relationship between rock fragment coverage and run‐off rate as well as a nonmonotonic relationship with respect to soil loss rate was recovered. The size of rock fragments also exerted a positive effect on run‐off generation and sediment yield but had a negative effect on infiltration. At the same time, both mean flow velocity and Froude number decreased with increasing rock fragment coverage regardless of rock fragment position and size, whereas both Manning roughness and Darcy–Weisbach friction factor were positively correlated. Results show that stream power is the most sensitive hydraulic parameter affecting soil loss. Combined with variance analysis, we concluded that the order of significance of rock fragment cover variables was position followed by coverage and then size. We also quantitatively incorporated the effects of rock fragment cover on soil loss via the C and K factors in the Revised Universal Soil Loss Equation. Overall, this study will enable the development of more accurate modelling approaches and lead to a better understanding of hydrological processes under rock fragment cover conditions.     

Hillslope run‐off thresholds with shrink–swell clay soils

Irrigation experiments on 12 instrumented field plots were used to assess the impact of dynamic soil crack networks on infiltration and run‐off. During applications of intensity similar to a heavy rainstorm, water was seen being preferentially delivered within the soil profile. However, run‐off was not observed until soil water content of the profile reached field capacity, and the apertures of surface‐connected cracks had closed >60%. Electrical resistivity measurements suggested that subsurface cracks persisted and enhanced lateral transport, even in wet conditions. Likewise, single‐ring infiltration measurements taken before and after irrigation indicated that infiltration remained an important component of the water budget at high soil water content values, despite apparent surface sealing. Overall, although the wetting and sealing of the soil profile showed considerable complexity, an emergent property at the hillslope scale was observed: all of the plots demonstrated a strikingly similar threshold run‐off response to the cumulative precipitation amount. Copyright © 2014 John Wiley & Sons, Ltd.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Effective run‐off flow length over biological soil crusts on silty loam soils in drylands

Surface hydrological behaviour is important in drylands because it affects the distribution of soil moisture and vegetation and the hydrological functioning of slopes and catchments. Microplot scale run‐off can be relatively easily measured, i.e. by rainfall simulations. However, slope or catchment run‐off cannot be deduced from microplots, requiring long‐time monitoring, because run‐off coefficients decrease with increasing drainage area. Therefore, to determine the slope length covered by run‐off (run‐off length) is crucial to connect scales. Biological soil crusts (BSCs) are good model systems, and their hydrology at slope scale is insufficiently known. This study provides run‐off lengths from BSCs, by field factorial experiments using rainfall simulation, including two BSC types, three rain types, three antecedent soil moistures and four plot lengths. Data were analysed by generalized linear modelling, including vascular plant cover as covariates. Results were the following: (i) the real contributing area is almost always much smaller than the topographical contributing area; (ii) the BSC type is key to controlling run‐off; run‐off length reached 3 m on cyanobacterial crust, but hardly over 1 m on lichen crust; this pattern remained through rain type or soil moisture; (iii) run‐off decreased with BSC development because soil sealing disappears; porosity, biomass and roughness increase and some changes occur in the uppermost soil layer; and (iv) run‐off flow increased with both rain type and soil moisture but run‐off coefficient only with soil moisture (as larger rains increased both run‐off and infiltration); vascular plant cover had a slight effect on run‐off because it was low and random. Copyright © 2014 John Wiley & Sons, Ltd.     

Surface moisture area during rainfall–run‐off events to understand the hydrological dynamics of a basin in a plain region

The understanding of the hydrology of plain basins may be improved by the combined analysis of rainfall–run‐off records and remote sensed surface moisture data. Our work evaluates the surface moisture area (SMA) produced during rainfall–run‐off events in a plain watershed of the Argentine Pampas Region, and studies which hydrological variables are related to the generated SMA. The study area is located in the upper and middle basins of the Del Azul stream, characterized by the presence of small gently hilly areas surrounded by flat landscapes. Data from 9 rainfall–run‐off events were analysed. MODIS surface reflectance data were processed to calculate SMA subsequent to the peak discharge (post‐SMA), and previous to the rainfall events (prev‐SMA), to consider the antecedent wetness. Rainfall–run‐off data included total precipitation depth (P), maximum intensity of rainfall over 6 hr (I6max), surface run‐off registered between the beginning of the event and the day previous to the analysed MODIS scene (R), peak flow (Qp), and flood intensity (IF). In contrast with other works, post‐SMA showed a negative relationship with the R. Three groups of cases were identified: (a) Events of low I6max, high prev‐SMA, and low R were associated with slow and weakly channelized flow over plain areas, leading to saturated overland flow (SOF), with large SMA; (b) events of high I6max, low prev‐SMA, and medium to high R were rapidly transported along the gentle slopes of the basin, related to Hortonian overland flow (HOF) and low post‐SMA; and (c) events of medium to high I6max and prev‐SMA with medium R were related to heterogeneous input‐antecedent‐run‐off conditions combined: Local spatial conditions may have produced HOF or SOF, leading to an averaged response with medium SMA. The interactions between the geomorphology of the basin, the characteristics of the events, and the antecedent conditions may explain the obtained results. This analysis is relevant for the general knowledge of the hydrology of large plains, whose functioning studies are still in their early stages.     

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

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

Sources of uncertainty in the NRCS CN model: Recognition and solutions

The Natural Resources Conservation Service Curve Number model is one of the most recognizable procedures in the field of rainfall–run‐off estimation. It has been widely applied for different purposes in hydrological models. In spite of its widespread use, some uncertainties have not even clarified and must be examined for its proper application. Particularly, choosing the most representative rainfall–run‐off events, and the coefficient λ that relates the parameters of the model (curve number CN and initial abstraction I_a). In this research, an advanced analysis is developed to evaluate the influence of λ for a set of representative watersheds of the Agricultural Research Service of the United Stated Department of Agriculture. They are characterized by different soil properties, land uses, and climatic conditions. Finally, 2 novel methodologies for the selection of the most representative rainfall–run‐off events and for the adaptation of coefficient λ are included, based on the pattern of rainfall distribution.     

Quantifying rainfall controls on catchment‐scale landslide erosion in Taiwan

Landslide erosion is a dominant hillslope process and the main source of stream sediment in tropical, tectonically active mountain belts. In this study, we quantified landslide erosion triggered by 24 rainfall events from 2001 to 2009 in three mountainous watersheds in Taiwan and investigated relationships between landslide erosion and rainfall variables. The results show positive power‐law relations between landslide erosion and rainfall intensity and cumulative rainfall, with scaling exponents ranging from 2·94 to 5·03. Additionally, landslide erosion caused by Typhoon Morakot is of comparable magnitude to landslide erosion caused by the Chi‐Chi Earthquake (M_W = 7·6) or 22–24 years of basin‐averaged erosion. Comparison of the three watersheds indicates that deeper landslides that mobilize soil and bedrock are triggered by long‐duration rainfall, whereas shallow landslides are triggered by short‐duration rainfall. These results suggest that rainfall intensity and watershed characteristics are important controls on rainfall‐triggered landslide erosion and that severe typhoons, like high‐magnitude earthquakes, can generate high rates of landslide erosion in Taiwan. Copyright © 2012 John Wiley & Sons, Ltd.     

On the relationship between flood and contributing area

Although it is well known that the vast majority of the time only a portion of any watershed contributes run‐off to the outlet, this extent is rarely documented. Also, the power law form of the streamflow and contributing area (Q‐A_c) relationship has been known for a half century, but it is uncommon for it to be quantified, and time series of contributing area extensive enough to calculate its frequency distribution are almost non‐existent. Data from the Canadian Prairies, where there are extensive estimates of contributing area during the median annual flood, imply that the power law coefficient for any Q‐A_c curve is a function of flow magnitude and return period. These data also suggest that regional flood frequency curves are a construct of Q‐A_c curves from individual basins. This paper will discuss research that attempted to reproduce the Q‐A_c curves for the La Salle River Watershed with a semidistributed numerical hydrological model, MESH‐PDMROF. The model simulated streamflow reasonably well (Nash Sutcliffe values = 0.62) compared with published examples of comparable models applied in the region. Estimates of the coefficient and exponent of the Q‐A_c power law function ranged from 0.08–0.14 and 0.9–1.12, respectively. These exponent values were lower than those of regional flood frequency curves and support the theory that regional flood frequency curves are a construct of Q‐A_c curves. Simulations of the area contributing to the median annual flood were lower (0.3) than those derived from independent topographic analysis (0.9) described in earlier literature though there is uncertainty in both these estimates. This uncertainty was extended across the flood frequency distribution and may be too large to definitively verify the study hypothesis.     

Tracing variability of run‐off generation in mountainous permafrost of semi‐arid north‐eastern Mongolia

The headwaters of mountainous, discontinuous permafrost regions in north‐eastern Mongolia are important water resources for the semi‐arid country, but little is known about hydrological processes there. Run‐off generation on south‐facing slopes, which are devoid of permafrost, has so far been neglected and is totally unknown for areas that have been affected by recent forest fires. To fill this knowledge gap, the present study applied artificial tracers on a steppe‐vegetated south‐facing and on two north‐facing slopes, burned and unburned. Combined sprinkling and dye tracer experiments were used to visualize processes of infiltration and water fluxes in the unsaturated zone. On the unburned north‐facing slope, rapid and widespread infiltration through a wet organic layer was observed down to the permafrost. On the burned profile, rapid infiltration occurred through a combusted organic and underlying mineral layer. Stained water seeped out at the bottom of both profiles suggesting a general tendency to subsurface stormflow (SSF). Ongoing SSF could directly be studied 24 h after a high‐intensity rainfall event on a 55‐m hillslope section in the burned forest. Measurements of water temperature proved the role of the permafrost layer as a base horizon for SSF. Repeated tracer injections allowed direct insights into SSF dynamics: A first injection suggested rather slow dispersive subsurface flow paths; whereas 18 h later, a second injection traced a more preferential flow system with 20 times quicker flow velocities. We speculate that these pronounced SSF dynamics are limited to burned slopes where a thermally insulating organic layer is absent. On three south‐facing soil profiles, the applied tracer remained in the uppermost 5 cm of a silt‐rich mineral soil horizon. No signs of preferential infiltration could be found, which suggested reduced biological activity under a harsh, dry and cold climate. Instead, direct observations, distributed tracers and charcoal samples provided evidence for the occurrence of overland flow. Copyright © 2014 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

Estimating daily inflows of large lakes using a water‐balance‐based runoff coefficient scaling approach

Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall–runoff relationship at a reference catchment was analysed to select rainfall input and estimate run‐off coefficient firstly, and the run‐off coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Run‐off coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared with the commonly used drainage area ratio method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R² and Nash–Sutcliffe efficiency values of 0.69 and 0.6, respectively. This approach outperforms drainage area ratio method in terms of mean R² and Nash–Sutcliffe efficiency values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.     

Run‐off formation in a humid,temperate headwater catchment using a combined hydrological,hydrochemical and isotopic approach (Jizera Mountains,Czech Republic)

The formation of baseflow and stormflow was examined in the 1.18 km² part of the headwater catchment Uhlí?ská, Jizera Mountains, Czech Republic, over the period 2007–2011, by means of run‐off data and environmental tracers ¹⁸O and SiO₂. The baseflow, computed using the digital filter approach BFLOW, contributes 67% to total streamflow and has a mean residence time of 12.3 months. It is formed by groundwater discharge from the valley deluviofluvial granitic sediments, in combination with soil water in weathered layers on hillslopes during rainfall and snowmelt periods. The prevailing source of the groundwater is the infiltration of snowmelt water. Analysis of 20 run‐off events and their hysteretic patterns demonstrated that the stormflow water has a residence time of about 4 months and is generated by preferential flow on hillslopes combined by soil matrix drainage. Because of slower flow in the soil matrix, the enrichment of pore water in SiO₂ is more pronounced. The stormflow and snowmelt water flowing via preferential pathways of upslope minerals soils pushes the pre‐event groundwater through the pathways in wetlands to the stream, and the wetland can be therefore considered as groundwater supplied. This mechanism has been found to be typical for the groundwater‐supplied headwater catchments of the Jizera Mountains and can be also assumed in other mountainous headwaters of the granitic massif in Central Europe. The main methodological contribution of this study are the residence time calculations stratified by baseflow and event flow, identifying run‐off components of different travel times to streams and linking them with geochemical run‐off sources. This achievement was possible because of a comprehensive dataset on hydrology, stable isotopes and silica hydrochemistry in all relevant run‐off generation components. This concept indicates that a possible long‐term change in snowmelt may affect the run‐off regime of headwater catchments to climate or land‐use changes. Copyright © 2013 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Occurrence of water ponding on soil surfaces depending on infiltration rates on Mongolian rangeland

The occurrence of water ponding on soil surfaces during and after heavy rainfall produces surface run‐off or surface water accumulation in low‐lying areas, which might reduce the water supply to soils and result in a reduction of the soil water that plants can use, especially in arid climates. On Mongolian rangeland, we observed ponded water on the surface of a specific soil condition subjected to a heavy rainfall of 30 mm/hr. By contrast, ponded water was not observed for the same type of soil where livestock grazing had been removed for 6–8 years via a fence or for nearby soil containing less clay. We measured the infiltration rate (the saturated hydraulic conductivity of the surface soil, K_s) of the three sites by applying ponded water on the soil surface (an intake rate test). The results showed that K_s in the rangeland was lower than the rainfall intensity in the site where water ponded on the soil surface; however, K_s of the soil inside of the fence has recovered to 3 times that of the soil outside of the fence to exceed the rainfall intensity. Heavy rainfall that exceeds the infiltration rate occurs several times a year at the livestock grazing site where we observed ponded water. Slight water repellency of the soil reduces rain infiltration to increase the possibility of surface ponding for the soil.     

The roles of channels and hillslopes in rainfall/run‐off lag times during intense storms in a steep catchment

The response time (lag time) between rainfall input and run‐off output in headwater catchments is a key parameter for flood prediction. Lag times are expected to be controlled by run‐off processes, both on hillslopes and in channels. To demonstrate these effects on peak lag times within a 4.5‐km² catchment, we measured stream water levels at up to 16 channel locations at 1‐min intervals and compared the lag times with topographic indices describing the length and gradient of the hillslope and channel flow path. We captured storm events with a total precipitation of 38–198 mm and maximum hourly precipitation intensity of 9–90 mm/hr. There were positive relationships between lag time and flow path length as well as the ratio of the flow path length and the square root of the gradient of channels for the most intense storms, demonstrating that channel flow paths generally defined the variation in lag times. Topographic analysis showed that hillslope flow path lengths were similar among locations, whereas channel flow path length increased almost one order of magnitude with a 100‐fold increase in catchment area. Thus, the relative importance of hillslope flow path decreased with increasing catchment area. Our results indicate that the variation in lag times is small when hillslopes are sufficiently wet; thus, catchment‐scale variation in lag times can be explained almost entirely by channel processes. Detailed topographic channel information can improve prediction of flood peak timing, whereas hillslopes can be treated as homogeneous during large flood events.     

Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.