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.     

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.     

Sodicity effects on hydrological processes of sodic soil slopes under simulated rainfall

Soil salinization can occur in many regions of the world. Soil sodicity affects rainfall‐runoff relationships and related erosion processes considerably. We investigated sodicity effects on infiltration, runoff and erosion processes on sodic soil slopes for two soils from China under simulated rainfall. Five sodicity levels were established in a silt loam and a silty clay with clay contents of 8.5% and 46.0%, respectively. The soils, packed in 50 cm × 30 cm × 15 cm flumes at two slope gradients (22° and 35°), were exposed to 60 min of simulated rainfall (deionized water) at a constant intensity of 125 mm h^?1. Results showed that, for both soils, increasing soil sodicity had some significant effects on hydrological processes, reducing the infiltration coefficient (pr = ?0.69, P  < 0.01) and the quasi‐steady final infiltration rate (pr = ?0.80, P  < 0.01), and increasing the mean sediment loss (pr = 0.39, P  < 0.05); however, it did not significantly affect the cumulative rainfall to ponding (P  > 0.05). Moreover, increasing sodicity significantly increased the Reynolds number and the stream power (pr = 0.78 and 0.66, P  < 0.01, respectively) of the runoff, decreased Manning roughness and Darcy–Weisbach coefficient (pr = ?0.52 and ?0.52, P  < 0.05, respectively), but did not significantly affect the mean flow velocity, mean flow depth, Froude number and hydraulic shear stress. Stream power was shown to be the most sensitive hydraulic variable affecting sediment loss for both soils. Furthermore, as sodicity increased, the values of critical stream power decreased for both the silt loam (R ² = 0.29, P  < 0.05) and the silty clay (R ² = 0.49, P  < 0.05). The findings of this study were applied to a real situation and identified some negative effects that can occur with increasing sodicity levels. This emphasized the importance of addressing the influences of soil sodicity in particularly high risk situations and when predicting soil and water losses.     

Determining short-term changes in the hydraulic properties of a sandy-loam soil by a three-run infiltration experiment

ABSTRACT

Soil structure-dependent parameters can vary rapidly as a consequence of perturbing events such as intense rainfall. Investigating their short-term changes is therefore essential to understand the general behaviour of a porous medium. The aim of this study is to gain insight into the effects of wetting, perturbation and recovery processes through different sequences of Beerkan infiltration experiments performed on a sandy-loam soil. Two different three-run infiltration experiments (LHL and LLL) were carried out by pouring water at low (L, non-perturbing) and high (H, perturbing) heights above the soil surface and at short time intervals (hours, days). The results demonstrate that the proposed method allows one to capture short-term variations in soil structure-dependent parameters. The developed methodology is expected to simplify the parameterization of hydrological models with temporally variable soil hydraulic properties.     

Evaluating water table response to rainfall events in a shallow aquifer and canal system

Shallow aquifers typically have greater hydrologic connectivity and response to recharge and changes in surface water management practices than deeper aquifers and are therefore often managed to reduce the risk of flooding. Quantification of the water table elevation response under different management scenarios provides valuable information in shallow aquifer systems to assess indirect influences of such modifications. The episodic master recession method was applied to the 15‐min water table elevation and NEXRAD rainfall data for 6 wells to identify water table response and individual rainfall events. The objectives of this study were to evaluate the effects of rainfall, water table elevation, canal stage, site‐specific characteristics, and canal structure modification/water management practice on the fluctuations in water table elevations using multiple/stepwise multiple linear regression techniques. With the modification of canal structure and operation adjustment, significant difference existed in water table response in the southern wells due to its relative downstream position regarding the general groundwater flow direction and the structural modification locations. On average, water table response height and flood risk were lower after than before the structure modification to canals. The effect of rainfall event size on the height of water table response was greater than the effect of antecedent water table elevation and canal stage on the height of water table response. Other factors including leakance of the canal bed sediment, specific yield, and rainfall on i  ? 1 day had significant effects on the height of water table response as well. Antecedent water table elevation and canal stage had greater and more linear effects on the height of water table response after the management changes to canals. Variation in water table response height/rainfall event size ratio was attributed to difference in S _y , antecedent soil water content, hydraulic gradient, rainfall size, and run‐off ratio. After the structure modification, water table response height/rainfall event size ratio demonstrated more linear and proportional relationship with antecedent water table elevation and canal stage.     

Beerkan multi-runs for characterizing water infiltration and spatial variability of soil hydraulic properties across scales

A method is presented for characterizing the spatial variability of water infiltration and soil hydraulic properties at the transect and field scales. The method involves monitoring a set of 10 Beerkan runs distributed over a 1-m length of soil, and running BEST (Beerkan estimation of soil transfer parameters) methods to derive hydraulic parameters. The Beerkan multi-runs (BMR) method provides a significant amount of data at the transect scale, allowing the determination of correlations between water infiltration variables and hydraulic parameters, and the detection of specific runs affected by preferential flow or water repellence. The realization of several BMRs at several transects on the same site allows comparison of the variation between locations (spatial variability at the field scale) and at the transect scale (spatial variability at the metre scale), using analysis of variance. From the results, we determined the spatial variability of water infiltration and hydraulic parameters as well as its characteristic scale (transect versus field).     

A comparison between the single ring pressure infiltrometer and simplified falling head techniques

Testing the relative performances of the single ring pressure infiltrometer (PI) and simplified falling head (SFH) techniques to determine the field saturated soil hydraulic conductivity, K_fs, at the near point scale may help to better establish the usability of these techniques for interpreting and simulating hydrological processes. A sampling of 10 Sicilian sites showed that the measured K_fs was generally higher with the SFH technique than the PI one, with statistically significant differences by a factor varying from 3 to 192, depending on the site. A short experiment with the SFH technique yielded higher K_fs values because a longer experiment with the PI probably promoted short‐term swelling phenomena reducing macroporosity. Moreover, the PI device likely altered the infiltration surface at the beginning of the run, particularly in the less stable soils, where soil particle arrangement may be expected to vary upon wetting. This interpretation was supported by a soil structure stability index, SSI, and also by the hydraulic conductivity data obtained with the tension infiltrometer, i.e. with a practically negligible disturbance of the sampled soil surface. In particular, a statistically significant, increasing relationship with SSI and an unsaturated conductivity greater than the saturated one were only detected for the K_fs data obtained with the PI. The SFH and PI techniques should be expected to yield more similar results in relatively rigid porous media (low percentages of fine particles and structurally stable soils) than in soils that modify appreciably their particle arrangement upon wetting. The simultaneous use of the two techniques may allow to improve K_fs determination in soils that change their hydrodynamic behaviour during a runoff producing rainfall event. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimation of soil hydraulic properties and their uncertainty through the Beerkan infiltration experiment

The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.     

Effects of biological crust coverage on soil hydraulic properties for the Loess Plateau of China

Biological soil crusts (BSCs) are ubiquitous living covers that have been allowed to grow on abandoned farmlands over the Loess Plateau because the “Grain for Green” project was implemented in 1999 to control serious soil erosion. However, few studies have been conducted to quantify the effects of BSC coverage on soil hydraulic properties. This study was performed to assess the effects of BSC coverage on soil hydraulic properties, which are reflected by the soil sorptivity under an applied pressure of 0 (S ₀ ) and ?3 (S ₃ ) cm, saturated hydraulic conductivity (K _s ), wetting front depth (WFD ), and mean pore radius (λ _m ), for the Loess Plateau of China. Five classes of BSC coverage (i.e., 1–20%, 20–40%, 40–60%, 60–80%, and 80–100%) and a bare control were selected at both cyanobacteria‐ and moss‐covered sites to measure soil hydraulic properties using a disc infiltrometer under 2 consecutive pressure heads of 0 and ?3 cm, allowing the direct calculation of S ₀ , S ₃ , K _s , and λ _m . The WFD was measured onsite using a ruler immediately after the experiments of infiltration. The results indicated that both cyanobacteria and moss crusts were effective in changing the soil properties and impeding soil infiltration. The effects of moss were greater than those of cyanobacteria. Compared to those of the control, the S ₀ , S ₃ , K _s , WFD , and λ _m values of cyanobacteria‐covered soils were reduced by 13.7%, 11.0%, 13.3%, 10.6%, and 12.6% on average, and those of moss‐covered soils were reduced by 27.6%, 22.1%, 29.5%, 22.2%, and 25.9%, respectively. The relative soil sorptivity under pressures of 0 (RS ₀ ) and ?3 (RS ₃ ) cm, the relative saturated hydraulic conductivity (RK _s ), the relative wetting front depth (RWFD ), and the relative mean pore radius (Rλ _m ) decreased exponentially with coverage for both cyanobacteria‐ and moss‐covered soils. The rates of decrease in RS ₀ , RS ₃ , RK _s , RWFD , and Rλ _m of cyanobacteria were significantly slower than those of moss, especially for the coverage of 0–40%, with smaller ranges. The variations of soil hydraulic properties with BSC coverage were closely related to the change in soil clay content driven by the BSC coverage on the Loess Plateau. The results are useful for simulating the hydraulic parameters of BSC‐covered soils in arid and semiarid areas.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Mapping potential surface ponding in agriculture using UAV-SfM

Among the environmental problems that could affect agriculture, one of the most critical is ponding. This may be defined as water storage on the surface in concavities and depressions due to soil saturation. Stagnant water can seriously affect crops and the management of agricultural landscapes. It is mainly caused by prolonged rainfall events, soil type, or wrong mechanization practices, which cause soil compaction. To better understand this problem and thus provide adequate solutions to reduce the related risk, high-resolution topographic information could be strategically important because it offers an accurate representation of the surface morphology. In the last decades, new remote sensing techniques provide interesting opportunities to understand the processes on the Earth's surface based on geomorphic signatures. Among these, Uncrewed Aerial Vehicles (UAVs), combined with the structure-from-motion (SfM) photogrammetry technique, represent a solid, low-cost, rapid, and flexible solution for geomorphological analysis. This study aims to present a new approach to detect the potential areas exposed to water stagnation at the farm scale. The high-resolution digital elevation model (DEM) from UAV-SfM data is used to do this. The potential water depth was calculated in the DEM using the relative elevation attribute algorithm. The detection of more pronounced concavities and convexities allowed an estimation and mapping of the potential ponding conditions. The results were assessed by observations and field measurements and are promising, showing a Cohen's k(X) accuracy of 0.683 for the planimetric extent of the ponding phenomena and a Pearson's r_xy coefficient of 0.971 for the estimation of pond water depth. The proposed workflow provides a useful indication to stakeholders for better agricultural management in lowland landscapes.     

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.     

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.     

FIELD ESTIMATION OF MACROPORE FUNCTIONING AND SURFACE HYDRAULIC CONDUCTIVITY IN A FEN PEAT

A limitation of existing models of water and solute movement in fen peats is that they fail to represent processes in the unsaturated zone. This limitation is largely due to a lack of data on the hydraulic properties of unsaturated peat, in particular the relationship between hydraulic conductivity (K) and pressure head (ψ). A tension infiltrometer was used to measure K(ψ) of a fen peat in Somerset, England. It was found that macropores could be important in water and solute movement in this soil type. It was also found that (i) variability of K in this peat was less than that reported for other peats and mineral soils, and (ii) the K data were better described by a log-normal distribution than a normal distribution in accord with findings from other peat and mineral soils. Recommendations on improving the understanding of water and solute movement in the unsaturated zone of this soil type are made. © 1997 by John Wiley & Sons, Ltd.     

Sediment and solute transport on soil slope under simultaneous influence of rainfall impact and scouring flow

Soil erosion and nutrient losses with surface runoff in the loess plateau in China cause severe soil quality degradation and water pollution. It is driven by both rainfall impact and runoff flow that usually take place simultaneously during a rainfall event. However, the interactive effect of these two processes on soil erosion has received limited attention. The objectives of this study were to better understand the mechanism of soil erosion, solute transport in runoff, and hydraulic characteristics of flow under the simultaneous influence of rainfall and shallow clear‐water flow scouring. Laboratory flume experiments with three rainfall intensities (0, 60, and 120 mm h⁻¹) and four scouring inflow rates (10, 20, 30, and 40 l min⁻¹) were conducted to evaluate their interactive effect on runoff. Results indicate that both rainfall intensity and scouring inflow rate play important roles on runoff formation, soil erosion, and solute transport in the surface runoff. A rainfall splash and water scouring interactive effect on the transport of sediment and solute in runoff were observed at the rainfall intensity of 60 mm h⁻¹ and scouring inflow rates of 20 l min⁻¹. Cumulative sediment mass loss (Ms) was found to be a linear function of cumulative runoff volume (Wr) for each treatment. Solute transport was also affected by both rainfall intensity and scouring inflow rate, and the decrease in bromide concentration in the runoff with time fitted to a power function well. Reynolds number (Re) was a key hydraulic parameter to determine erodability on loess slopes. The Darcy–Weisbach friction coefficients (f) decreased with the Reynolds numbers (Re), and the average soil and water loss rate (M_l) increased with the Reynolds numbers (Re) on loess slope for both scenarios with or without rainfall impact. Copyright © 2010 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytical solutions of infiltration process under ponding irrigation

The objective of this paper is to simulate the progress of the soil water content distribution in the soil profile with a water table at the bottom of the soil profile during ponding irrigation. This simulation can be done by solving the two‐dimensional Richards's equation for the assimilation of the advancing water jet, which uses the conditions of the two exponential functional forms k = k_s e^αψ and θ = θ_r + (θ_s − θ_r) e^αψ to represent the hydraulic conductivity and volumetric water content, with ψ the pressure as the third variable. We assume that the ground surface becomes ponded and saturated as soon as the water flux passes the dry ground surface. By the technique of transformation, the analytical solution of these two‐dimensional Richards' equations has enabled figures of volumetric water content distribution to be obtained in successive time periods after irrigation. For the example of loam soil, it can simulate the variation of volumetric water content during and after irrigation in the soil profile. The analytical solutions of this paper reflect the real situation simulated, and can be applied to verify those complicated solutions from other analytical models. Copyright © 2005 John Wiley & Sons, Ltd.     

Percolation losses in paddy fields with a dynamic soil structure: model development and applications

The hydraulic characteristics of the plough pan of paddy fields provide continuous ponding conditions during the growing season and control the water use efficiency in wet rice production. Its saturated hydraulic conductivity K_s, however, exhibits a large spatiotemporal variability as a consequence of a highly dynamic soil structure involving temporary shrinkage cracks. Water flow through the earthen bunds surrounding the fields further contributes to the uncertainty in water flux calculations. The objective of this study was to develop a simple deterministic model with stochastic elements (‘PADDY‐FLUX’) for depiction of deep percolation, and to assess the effect of different water management scenarios on percolation in two channel command areas. Darcy's law is used as the fundamental equation for water flow calculations with the ponding water depth h as a time‐dependent variable. Flux uncertainty is estimated by a Monte‐Carlo‐type implementation. K_s is treated as a random variable of a bimodal probability density function (PDF), which is the weighted sum of two Gaussian PDFs (accounting for a matrix and a preferential flow domain). The weighing factor α is a function of h, reflecting an increasing risk for preferential flow situations after desiccation and the development of shrinkage cracks. Under‐bund percolation is calculated using transfer functions. The results demonstrate that percolation losses increase in the following order: continuous soil saturation < continuous flooding (CF) < mid‐season drainage and intermittent irrigation (MD + II) < mid‐season drainage and continuous flooding. The bunds contribute up to 54 and 17% to total fluxes under CF and MD + II, respectively. Preferential water fluxes are responsible for the major part of water losses as soon as desiccation causes the formation of shrinkage cracks. As a conclusion, continuous soil saturation should be promoted as the least water‐intensive irrigation regime, while intermittent irrigation is recommended only in case that irreversible shrinkage cracks have already developed. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing the impacts of microtopography on soil erosion under simulated rainfall,using a multifractal approach

This study aimed to investigate the changing characteristics of microrelief of purple soil and its erosional response during successive stages of water erosion, including splash erosion, sheet erosion, and rill erosion. Methods employed included a rainfall simulator and the use of a laser scanner to generate a digital elevation model. Three artificial tillage practices, including conventional tillage (CT), artificial digging (AD), and ridge tillage (RT), were used to simulate different microrelief patterns. Eighteen artificial rainfall experiments were conducted using three 2 × 1 m boxes with a rainfall intensity of 1.5 mm min^?1 on a 15° slope. The results showed that the soil roughness (SR) index values for the tillage slopes were RT > AD > CT. The combined effects of detachment by raindrop impact and transport by run‐off decreased the SR index, whereas rill erosion increased the SR index during rainfall event. Microtopography and drainage networks have strong multifractal behaviours. The multifractal parameters of microtopography reflect the overall characteristics as well as the characteristics of the local soil surface. Within a certain range of threshold values, higher microrelief causes less soil erosion. However, when the parameters of spatial heterogeneity of microtopography exceed the threshold values, a higher degree of microrelief can increase soil erosion. These results help clarify the effect of microtopography on soil erosion and provide a theoretical foundation to guide future tillage practices on sloping farmland of purple soil.