Seasonal variations in infiltrability of moss‐dominated biocrusts on aeolian sand and loess soil in the Chinese Loess Plateau

Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (K_s) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased K_s both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P < .05) and on loess (by 42%, 92%, and 10%, respectively; P <.05). As compared with the bare soil, the decreased K_s in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in K_s, but their patterns were quite different. Specifically, the K_s of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the K_s of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in K_s of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of K_s were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.     

Variation in soil saturated hydraulic conductivity along the hillslope of collapsing granite gullies

Saturated hydraulic conductivity (K_s) affects the soil hydrological process and is influenced by many factors that exhibit strong spatial variations. To accurately measure K_s and its scale, spatial variability and relationship with collapsing gullies, we analysed four double-ring infiltrometer diameters in three soil layers during in situ experiments designed to measure K_s in two typical collapsing gullies (three slope sites) in Tongcheng County of China. The results showed that K_s increased with increasing inner ring diameter, but no significant difference existed between inner diameters of 30 and 40 cm. The K_s in red soil layers was higher than that in sandy soil layers, the transition layers had the lowest values. K_s also varied with slope position, gradually decreasing with distance from the gully head. The suggestion is that the spatial variation in K_s is affected not only by the intrinsic soil properties but also by the interaction with the collapsing gully.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

A modified multiple tension upward infiltration method to estimate the soil hydraulic properties

Determination of saturated hydraulic conductivity, K_s, and the shape parameters α and n of the water retention curve, θ(h), is of paramount importance to characterize the water flow in the vadose zone. This work presents a modified upward infiltration method to estimate K_s, α and n from numerical inverse analysis of the measured cumulative upward infiltration (CUI) at multiple constant tension lower boundary conditions. Using the HYDRUS‐2D software, a theoretical analysis on a synthetic loam soil under different soil tensions (0, 0–10, 0–50 and 0–100 cm), with and without an overpressure step of 10 cm high from the top boundary condition at the end of the upward infiltration process, was performed to check the uniqueness and the accuracy of the solutions. Using a tension sorptivimeter device, the method was validated in a laboratory experiment on five different soils: a coarse and a fine sand, and a 1‐mm sieved loam, clay loam and silt‐gypseous soils. The estimated α and n parameters were compared to the corresponding values measured with the TDR‐pressure cell method. The theoretical analysis demonstrates that K_s and θ(h) can be simultaneously estimated from measured upward cumulative infiltration when high (>50 cm) soil tensions are initially applied at the lower boundary. Alternatively, satisfactory results can be also obtained when medium tensions (<50 cm) and the K_s calculated from the overpressure step at the end of the experiment are considered. A consistent relationship was found between the α (R² = 0.86, p < 0.02) and n (R² = 0.97, p < 0.001) values measured with the TDR‐pressure cell and the corresponding values estimated with the tension sorptivimeter. The error between the α (in logarithm scale) and n values estimated with the inverse analysis and the corresponding values measured with pressure chamber were 3.1 and 6.1%, respectively. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimating field‐saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment

Field‐saturated soil hydraulic conductivity, K_fs, is highly variable. Therefore, interpreting and simulating hydrological processes, such as rainfall excess generation, need a large number of K_fs data even at the plot scale. Simple and reasonably rapid experiments should be carried out in the field. In this investigation, a simple infiltration experiment with a ring inserted shortly into the soil and the estimation of the so‐called α* parameter allowed to obtain an approximate measurement of K_fs. The theoretical approach was tested with reference to 149 sampling points established on Burundian soils. The estimated K_fs with the value of first approximation of α* for most agricultural field soils (α* = 0.012 mm^?1) differed by a practically negligible maximum factor of two from the saturated conductivity obtained by the complete Beerkan Estimation of Soil Transfer parameters (BEST) procedure for soil hydraulic characterization. The measured infiltration curve contained the necessary information to obtain a site‐specific prediction of α*. The empirically derived α* relationship gave similar results for K_fs (mean = 0.085 mm s^?1; coefficient of variation (CV) = 71%) to those obtained with BEST (mean = 0.086 mm s^?1; CV = 67%), and it was also successfully tested with reference to a few Sicilian sampling points, since it yielded a mean and a CV of K_fs (0.0094 mm s^?1 and 102%, respectively) close to the values obtained with BEST (mean = 0.0092 mm s^?1; CV = 113%). The developed method appears attractive due to the extreme simplicity of the experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Upward infiltration–evaporation method to estimate soil hydraulic properties

Determination of saturated hydraulic conductivity, K_s, and the van Genuchten water retention curve θ(h) parameters is crucial in evaluating unsaturated soil water flow. The aim of this work is to present a method to estimate K_s, α and n from numerical analysis of an upward infiltration process at saturation (Cap0), with (Cap0 + h) and without (Cap0) an overpressure step (h) at the end of the wetting phase, followed by an evaporation process (Evap). The HYDRUS model as well as a brute-force search method were used for theoretical loam soil parameter estimation. The uniqueness and the accuracy of solutions from the response surfaces, K_s–n, α–n and K_s–α, were evaluated for different scenarios. Numerical experiments showed that only the Cap0 + Evap and Cap0 + h + Evap scenarios were univocally able to estimate the hydraulic properties. The method gave reliable results in sand, loam and clay-loam soils.     

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.     

Predicting physical equations of soil detachment by simulated concentrated flow in Ultisols (subtropical China)

Soil detachment in concentrated flow is due to the dislodging of soil particles from the soil matrix by surface runoff. Both aggregate stability and shear strength of the topsoil reflect the erosion resistance of soil to concentrated runoff, and are important input parameters in predicting soil detachment models. This study was conducted to develop a formula to predict soil detachment rate in concentrated flow by using the aggregate stability index (As), root density (R_d) and saturated soil strength (σ_s) in the subtropical Ultisols region of China. The detachment rates of undisturbed topsoil samples collected from eight cultivated soil plots were measured in a 3.8 m long, 0.2 m wide hydraulic flume under five different flow shear stresses (τ = 4.54, 9.38, 15.01, 17.49 and 22.54 Pa). The results indicated that the stability index (As) was well related with soil detachment rate, particularly for results obtained with high flow shear stress (22.54 Pa), and the stability index (As) has a good linear relationship with concentrated flow erodibility factors (K_c). There was a positive linear relationship between saturated soil strength (σ_s) and critical flow shear stress (τ_c) for different soils. A significant negative exponential relationship between erodibility factors (K_c) and root density (R_d) was detected. This study yielded two prediction equations that allowed comparison of their efficiency in assessing soil detachment rate in concentrated flow. The equation including the root density (R_d) may have a better correlation coefficient (R² = 0.95). It was concluded that the formula based on the stability index (As), saturated soil strength (σ_s) and root density (R_d) has the potential to improve methodology for assessing soil detachment rate in concentrated flow for the subtropical Chinese Ultisols. Copyright © 2012 John Wiley & Sons, Ltd.     

Measuring streambank erosion due to ground water seepage: correlation to bank pore water pressure,precipitation and stream stage

Limited information exists on one of the mechanisms governing sediment input to streams: streambank erosion by ground water seepage. The objective of this research was to demonstrate the importance of streambank composition and stratigraphy in controlling seepage flow and to quantify correlation of seepage flow/erosion with precipitation, stream stage and soil pore water pressure. The streambank site was located in Northern Mississippi in the Goodwin Creek watershed. Soil samples from layers on the streambank face suggested less than an order of magnitude difference in vertical hydraulic conductivity (K_s) with depth, but differences between lateral K_s of a concretion layer and the vertical K_s of the underlying layers contributed to the propensity for lateral flow. Goodwin Creek seeps were not similar to other seeps reported in the literature, in that eroded sediment originated from layers underneath the primary seepage layer. Subsurface flow and sediment load, quantified using 50 cm wide collection pans, were dependent on the type of seep: intermittent low‐flow (LF) seeps (flow rates typically less than 0·05 L min^?1), persistent high‐flow (HF) seeps (average flow rate of 0·39 L min^?1) and buried seeps, which eroded unconsolidated bank material from previous bank failures. The timing of LF seeps correlated to river stage and precipitation. The HF seeps at Goodwin Creek began after rainfall events resulted in the adjacent streambank reaching near saturation (i.e. soil pore water pressures greater than ?5 kPa). Seep discharge from HF seeps reached a maximum of 1·0 L min^?1 and sediment concentrations commonly approached 100 g L^?1. Buried seeps were intermittent but exhibited the most significant erosion rates (738 g min^?1) and sediment concentrations (989 g L^?1). In cases where perched water table conditions exist and persistent HF seeps occur, seepage erosion and bank collapse of streambank sediment may be significant. Copyright © 2007 John Wiley & Sons, Ltd.     

Near‐saturated surface soil hydraulic properties under different land uses in the St Denis National Wildlife Area,Saskatchewan, Canada

Surface soil hydraulic properties are key factors controlling the partition of rainfall and snowmelt into runoff and soil water storage, and their knowledge is needed for sound land management. The objective of this study was to evaluate the effects of three land uses (native grass, brome grass and cultivated) on surface soil hydraulic properties under near‐saturated conditions at the St Denis National Wildlife Area, Saskatchewan, Canada. For each land use, water infiltration rates were measured using double‐ring and tension infiltrometers at ?0·3, ?0·7, ?1·5 and ?2·2 kPa pressure heads. Macroporosity and unsaturated hydraulic properties of the surface soil were estimated. Mean field‐saturated hydraulic conductivity (K_fs), unsaturated hydraulic conductivity at ?0·3 kPa pressure head, inverse capillary length scale (α) and water‐conducting macroporosity were compared for different land uses. These parameters of the native grass and brome grass sites were significantly (p < 0·1) higher than that of the cultivated sites. At the ?0·3 kPa pressure head, hydraulic conductivity of grasslands was two to three times greater than that of cultivated lands. Values of α were about two times and values of K_fs about four times greater in grasslands than in cultivated fields. Water‐conducting macroporosity of grasslands and cultivated fields were 0·04% and 0·01% of the total soil volume, respectively. Over 90% of the total water flux at ?0·06 kPa pressure head was transmitted through pores > 1·36 × 10^?4 m in diameter in the three land uses. Land use modified near‐saturated hydraulic properties of surface soil and consequently may alter the water balance of the area by changing the amount of surface runoff and soil water storage. Copyright © 2004 John Wiley & Sons, Ltd.     

Alternative analysis of transient infiltration experiment to estimate soil water repellency

The repellency index (RI) defined as the adjusted ratio between soil‐ethanol, S_e, and soil‐water, S_w, sorptivities estimated from minidisk infiltrometer experiments has been used instead of the widely used water drop penetration time and molarity of ethanol drop tests to assess soil water repellency. However, sorptivity calculated by the usual early‐time infiltration equation may be overestimated as the effects of gravity and lateral capillary are neglected. With the aim to establish the best applicative procedure to assess RI, different approaches to estimate S_e and S_w were compared that make use of both the early‐time infiltration equation (namely, the 1 min, S1, and the short‐time linearization approaches), and the two‐term axisymmetric infiltration equation, valid for early to intermediate times (namely, the cumulative linearization and differentiated linearization approaches). The dataset included 85 minidisk infiltrometer tests conducted in three sites in Italy and Spain under different vegetation habitats (forest of Pinus pinaster and Pinus halepensis, burned pine forest, and annual grasses), soil horizons (organic and mineral), postfire treatments, and initial soil water contents. The S1 approach was inapplicable in 42% of experiments as water infiltration did not start in the first minute. The short‐time linearization approach yielded a systematic overestimation of S_e and S_w that resulted in an overestimation of RI by a factor of 1.57 and 1.23 as compared with the cumulative linearization and differentiated linearization approaches. A new repellency index, RI_s, was proposed as the ratio between the slopes of the linearized data for the wettable and hydrophobic stages obtained by a single water infiltration test. For the experimental conditions considered, RI_s was significantly correlated with RI and WDPT. Compared with RI, RI_s includes information on both soil sorptivity and hydraulic conductivity and, therefore, it can be considered more physically linked to the hydrological processes affected by soil water repellency.     

A field‐scale infiltration model accounting for spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

This study first explores the role of spatial heterogeneity, in both the saturated hydraulic conductivity K_s and rainfall intensity r, on the integrated hydrological response of a natural slope. On this basis, a mathematical model for estimating the expected areal‐average infiltration is then formulated. Both K_s and r are considered as random variables with assessed probability density functions. The model relies upon a semi‐analytical component, which describes the directly infiltrated rainfall, and an empirical component, which accounts further for the infiltration of surface water running downslope into pervious soils (the run‐on effect). Monte Carlo simulations over a clay loam soil and a sandy loam soil were performed for constructing the ensemble averages of field‐scale infiltration used for model validation. The model produced very accurate estimates of the expected field‐scale infiltration rate, as well as of the outflow generated by significant rainfall events. Furthermore, the two model components were found to interact appropriately for different weights of the two infiltration mechanisms involved. Copyright © 2005 John Wiley & Sons, Ltd.     

Identification of optimal soil hydraulic functions and parameters for predicting soil moisture

Abstract

The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, K_s , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to K_s than to the soil pore-size distribution parameter.

Editor D. Koutsoyiannis; Associate editor D. Hughes

Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737.     

A semi-analytical model of expected areal-average infiltration under spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

A semi-analytical model for the estimate of expected areal-average infiltration rate at hillslope scale is presented. It accounts for spatial heterogeneity of the saturated hydraulic conductivity, K_s, and rainfall rate, r. The K_s field is characterized by a lognormal probability density function while the rainfall rate r is represented by a uniform distribution between two extreme values. The model formulation relies upon the use of cumulative infiltration as the independent variable which is then expressed as a function of an expected time for use in practical applications. The solution is applicable for those ranges of r and K_s that allow for neglecting the infiltration of surface water running downslope into pervious soils (run-on process). The model was tested by comparisons with Monte Carlo simulations carried out for a variety of coefficients of variation of r and K_s over a clay loam soil and a sandy loam soil. The model was found to be very reliable both with coupled spatial variability of r and K_s and when only one variable is characterized by spatial heterogeneity while the other is uniform.     

Geomorphic controls on the physical and hydrologic properties of soils in a valley floor

The distribution of soil hydraulic and physical properties strongly influences runoff processes in landscapes. Although much work has been done to quantify and predict the properties of hillslope soils, far less is known about the distribution of soil properties in valley floors. A technique that links the estimation and distribution of soil hydraulic properties in valleys, with easily identified geomorphic features, was developed along a 2 km length of a valley at Brooks Creek in New South Wales, Australia. Soil physical and hydraulic property data were collected across a set of floodplain and fan features within the valley and analysed statistically to determine if soil properties varied significantly between geomorphic features and stratigraphic layers. The results show that the depth‐averaged saturated hydraulic conductivity, K_s, of the soil varies significantly with landform: fan units have K_g values that are twice that of floodplains and colluvial toeslope deposits have K_s values four times higher than floodplains. Given the notorious variability of K_s values in space, the strong statistical separation of soil properties by landform, backed up by strong separation of soil particle size by landform, suggests a way forward in understanding the distribution of soil properties in valleys and their influence on catchment hydrology. Copyright © 2000 John Wiley & Sons, Ltd.     

An approach to analysing plot scale infiltration and runoff responses to rainfall of fluctuating intensity

Simulated rainfall of fluctuating intensity was applied to runoff plots on bare dryland soils in order to explore a new method for analysing the non‐steady‐state responses of infiltration and overland flow. The rainfall events all averaged 10 mm/h but included intensity bursts of up to 70 mm/h and lasting 5–15 min, as well as periods of low intensity and intermittency of up to 25 min. Results were compared with traditional steady‐state estimates of infiltrability made under simulated rainfall sustained at a fixed intensity of 10 mm/h. Mean event infiltration rate averaged 13.6% higher under fluctuating intensities, while runoff ratios averaged only 63% of those seen under constant intensity. In order to understand the changing soil infiltrability, up to three affine Horton infiltration equations were fitted to segments of each experiment. All equations had the same final infiltrability f_c, but adjusted values for coefficients f₀ (initial infiltrability) and K_f (exponential decay constant) were fitted for periods of rainfall that followed significant hiatuses in rainfall, during which subsurface redistribution allowed near‐surface soil suction to recover. According to the fitted Horton equations, soil infiltrability recovered by up 10–24 mm/h during intra‐event rainfall hiatuses of 15 to 20‐min duration, contributing to higher overall event infiltration rates and to reduced runoff ratios. The recovery of infiltrability also reduced the size of runoff peaks following periods of low intensity rainfall, compared with the predictions based on single Horton infiltration equations, and in some cases, no runoff at all was recorded from late intensity peaks. The principal finding of this study is that, using a set of affine equations, the intra‐event time variation of soil infiltrability can be tracked through multiple intensity bursts and hiatuses, despite the lack of steady‐state conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Using a transient infiltrometric technique for intensively sampling field‐saturated hydraulic conductivity of a clay soil in two runoff plots

The point measurement of soil properties allows to explain and simulate plot scale hydrological processes. An intensive sampling was carried out at the surface of an unsaturated clay soil to measure, on two adjacent plots of 4 × 11 m² and two different dates (May 2007 and February–March 2008), dry soil bulk density, ρ_b, and antecedent soil water content, θ_i, at 88 points. Field‐saturated soil hydraulic conductivity, K_fs, was also measured at 176 points by the transient Simplified Falling Head technique to determine the soil water permeability characteristics at the beginning of a possible rainfall event yielding measurable runoff. The ρ_b values did not differ significantly between the two dates, but wetter soil conditions (by 31%) and lower conductivities (1.95 times) were detected on the second date as compared with the first one. Significantly higher (by a factor of 1.8) K_fs values were obtained with the 0.30‐m‐diameter ring compared with the 0.15‐m‐diameter ring. A high K_fs (> 100 mm h^?1) was generally obtained for low θ_i values (< 0.3 m³m^?3), whereas a high θ_i yielded an increased percentage of low K_fs data (1–100 mm h^?1). The median of K_fs for each plot/sampling date combination was not lower than 600 mm h^?1, and rainfall intensities rarely exceeded 100 mm h^?1 at the site. The occurrence of runoff at the base of the plot needs a substantial reduction of the surface soil permeability characteristics during the event, probably promoted by a higher water content than the one of this investigation (saturation degree = 0.44–0.62) and some soil compaction due to rainfall impact. An intensive soil sampling reduces the risk of an erroneous interpretation of hydrological processes. In an unstable clay soil, changes in K_fs during the event seem to have a noticeable effect on runoff generation, and they should be considered for modeling hydrological processes. Copyright © 2012 John Wiley & Sons, Ltd.     

INVESTIGATION OF INFILTRATION CHARACTERISTICS AND DEBRIS FLOW INITIATION CONDITIONS IN DEBRIS FLOW SOURCE AREAS USING A RAINFALL SIMULATOR

Rainfall was simulated on unconfined plots on regolith in debris flow source areas using a portable simulator. In total, 351 simulations were carried out on steep slopes (27–54°) with rainfall intensities of 28–291 mm/h. From these rainfall simulations the infiltration parameters sorptivity (S) and steady-state infiltration capacity (K) of the regolith, and a threshold for the occurrence of micro-scale mass movements, were obtained. Two evaluation methods were used to obtain the infiltration parameters K and S. The ‘infiltration envelope’ method uses rainfall intensity and time to ponding from multiple tests and fits an infiltration envelope through the data from which K and S can be obtained. The ‘constant runoff’ method uses rainfall intensity and overland flow intensity to calculate K, after which S can be calculated in several ways by using time to ponding. The constant runoff method produced K values of 16.6–128 mm/h, which usually show a log-normal distribution. K values depend on the regolith parent material and rainfall intensity. Using this method, S values are 0.088–0.381 cm/min^1/2. The infiltration envelope method produced K values of 9.8–131 mm/h and S values of 0.14–0.32 cm/min^1/2. It can be argued that both methods overestimate K as well as S, but quantitative relations between measured/calculated and actual values of K and S have not yet been obtained. At high rainfall intensities, typically 100 mm/h or more, micro-scale mass movements sometimes occur. A lower threshold curve for the occurrence of these micro-scale mass movements has been constructed. It is a function of both slope angle and rainfall intensity. The micro-scale mass movements could play an important part in the initiation of debris flows in the study area, possibly by delivering sediment to overland flow. On the very steep slopes, the sediment-rich overland flow can easily mobilize coarse material.