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.     

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.     

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.     

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.     

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.     

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.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil hydrological response under simulated rainfall in the Dehesa land system (Extremadura,SW Spain) under drought conditions

Soil hydrology was investigated in the Guadelperalón experimental watershed in order to determine the influence of land use and vegetation cover on runoff and infiltration within the Dehesa land system. Five soil–vegetation units were selected: (1) tree cover, (2) sheep trials, (3) shrub cover, (4) hillslope grass and (5) bottom grass. The results of the simulated rainfall experiments performed at an intensity of 56·6 mm h⁻¹ during one hour on plots of 0·25 m², and the water drop penetration time test indicate the importance of water repellency in the Dehesa land system under drought conditions. Low infiltration rates (c. 9–44 mm h⁻¹) were found everywhere except at shrub sites and in areas with low grazing pressure. Soil water repellency greatly reduced infiltration, especially beneath Quercus ilex canopies, where fast ponding and greater runoff rates were observed. The low vegetation cover as a consequence of a prolonged drought and grazing pressure, in conjunction with the soil water repellency, induces high runoff rates (15–70 per cent). In spite of this, macropore fluxes were found in different locations, beneath trees, on shrub-covered surfaces, as well as at sites with a dominance of herbaceous cover. Discontinuity of the runoff fluxes due to variations in hydrophobicity causes preferential flows and as a consequence deeper infiltration, especially where macropores are developed. © 1998 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.     

Development and parameterization of an infiltration model accounting for water depth and rainfall intensity

Experimental work has clearly shown that the effective hydraulic conductivity (K_e) or effective infiltration rate (f_e) on the local scale of a plot cannot be considered as constant but are dependent on water depth and rainfall intensity because non‐random microtopography‐related variations in hydraulic conductivity occur. Rainfall–runoff models generally do not account for this: models assume that excess water is uniformly spread over the soil surface and within‐plot variations are neglected. In the present study, we propose a model that is based on the concepts of microtopography‐related water depth‐dependent infiltration and partial contributing area. Expressions for the plot scale K_e and f_e were developed that depend on rainfall intensity and runon from upslope (and thus on water depth). To calibrate and validate the model, steady state infiltration experiments were conducted on maize fields on silt loam soils in Belgium, with different stages and combinations of rainfall intensity and inflow, simulating rainfall and runon. Water depth–discharge and depth–inundation relationships were established and used to estimate the effect of inundation on K_e. Although inflow‐only experiments were found to be unsuitable for calibration, the model was successfully calibrated and validated with the rainfall simulation data and combined rainfall–runon data (R²: 0.43–0.91). Calibrated and validated with steady state infiltration experiments, the model was combined with the Green–Ampt infiltration equation and can be applied within a two‐dimensional distributed rainfall–runoff model. The effect of water depth–dependency and rainfall intensity on infiltration was illustrated for a hillslope. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of rainfall movement on soil crust development

The impact of rainfall event movement on soil crust development was investigated via laboratory experiments and physics-based hydrological-response simulation.Runoff and soil water content data were analyzed in order to study different crustal development periods under a variety of precipitation conditions.The results showed that the formation of crust significantly changed soil saturated hydraulic conductivity(K_s) and runoff,while a steady state,which could be broken by an extended period of raindrop impact,was found to exist for both the crust and,subsequently,runoff characteristics.The decrease in K_S was as high as 86%for downstream-moving events,and 72%for upstream-moving events,illustrating that the direction of rainfall movement has a more significant impact on crust formation than rainfall intensity and duration.As downstream-moving rainfall events may have a greater effect than those moving upstream,the development of a depositional crust before raindrop impact may play an important role in determining the final K_S value.     

Comparing Beerkan infiltration tests with rainfall simulation experiments for hydraulic characterization of a sandy‐loam soil

Saturated soil hydraulic conductivity, K _s , data collected by ponding infiltrometer methods and usual experimental procedures could be unusable for interpreting field hydrological processes and particularly rainfall infiltration. The K _s values determined by an infiltrometer experiment carried out by applying water at a relatively large distance from the soil surface could however be more appropriate to explain surface runoff generation phenomena during intense rainfall events. In this study, a link between rainfall simulation and ponding infiltrometer experiments was established for a sandy‐loam soil. The height of water pouring for the infiltrometer run was chosen, establishing a similarity between the gravitational potential energy of the applied water, E _p , and the rainfall kinetic energy, E _k . To test the soundness of this procedure, the soil was sampled with the Beerkan estimation of soil transfer parameters procedure of soil hydraulic characterization and two heights of water pouring (0.03 m, i.e., usual procedure, and 0.34 m, yielding E _p  = E _k ). Then, a comparison between experimental steady‐state infiltration rates, i _sR , measured with rainfall simulation experiments determining runoff production and K _s values for the two water pouring heights was carried out in order to discriminate between theoretically possible (i _sR  ≥ K _s ) and impossible (i _sR  < K _s ) situations. Physically possible K _s values were only obtained by applying water at a relatively large distance from the soil surface, because i _sR was equal to 20.0 mm h^?1 and K _s values were 146.2–163.9 and 15.2–18.7 mm h^?1 for a height of water pouring of 0.03 and 0.34 m, respectively. This result suggested the consistency between Beerkan runs with a high height of water pouring and rainfall simulator experiments. Soil compaction and mechanical aggregate breakdown were the most plausible physical mechanisms determining reduction of K _s with height. This study demonstrated that the height from which water is poured onto the soil surface is a key parameter in infiltrometer experiments and can be adapted to mimic the effect of high intensity rain on soil hydraulic properties.     

Effects of rainfall intensity and compaction on water transport from opencast coal mine soils: An experimental study

The use of heavy machinery during opencast coal mining can result in soil compaction. Severe soil compaction has a negative impact on the transport of water and gas in the soil. In addition, rainfall intensity has traditionally been related to soil surface sealing affecting water transport. To assess the effects of rainfall intensity and compaction on water infiltration and surface runoff in an opencast coal mining area, the disturbed soils from the Antaibao opencast mine in Shanxi Province, China, were collected. Four soil columns with different bulk densities (i.e., 1.4 g cm^-3, 1.5 g cm^-3, 1.6 g cm^-3, and 1.7 g cm^-3) were designed, and each column received water five times at rainfall intensities of 23.12, 28.91, 38.54, 57.81, and 115.62 mm hr^-1. The total volume of runoff, the time to start runoff, and the volumetric water contents at the depths of 5 cm, 15 cm, 25 cm, 35 cm, 45 cm, 55 cm, and 65 cm were measured. Under the same soil bulk density, high rainfall intensity reduced infiltration, increased surface runoff, and decreased the magnitude of change in the volumetric water contents at different depths. Under the same rainfall intensity, the soil column with a high bulk density showed relatively low water infiltration. Treatments 3 (1.6 g cm^-3) and 4 (1.7 g cm^-3) had very small changes in volumetric water contents of the profiles even under a lower rainfall intensity. Severe soil compaction was highly prone to surface runoff after rainfall. Engineering and revegetation measures are available to improve compacted soil quality in dumps. Our results provide a theoretical basis for the management of land reclamation in opencast coal mine areas.     

Groundwater impact on methane emissions from flooded paddy fields

High methane (CH₄) fluxes emitted from paddy fields strongly contribute to the accumulation of greenhouse gases into the atmosphere, compromising the eco-compatibility of one of the most important world foods. A strong link exists between infiltration rates of irrigation water and CH₄ emissions. Since depth to the groundwater table affects infiltration rates, a relevant groundwater impact is expected on CH₄ emissions from paddy fields. In this work, a theoretical approach is adopted to investigate the aquifer effect on CH₄ dynamics in paddies. Infiltration rates are strongly affected by the development of different connection states between aquifer and irrigation ponded water. A strong reduction in infiltration rates results from a water table near to the soil surface, when the system is hydraulically connected. When the groundwater level increases, the infiltration rate reduction due to the switch from disconnected to connected state promotes a relevant increase of CH₄ emissions. This is due to a strong reduction of dissolved organic carbon (DOC) percolation, which leads to higher DOC availability for microbial CH₄ production and, consequently, higher CH₄ emissions. Our simulations show that CH₄ fluxes can be reduced by up to 24% when groundwater level is decreased and the aquifer is disconnected from ponding water. In paddies with shallow aquifers, lowering the water table with a drainage system could thus represent a promising CH₄ mitigation option.     

Effects of initial soil water content and saturated hydraulic conductivity variability on small watershed runoff simulation using LISEM

Abstract

Soil water content (θ) and saturated hydraulic conductivity (K_s) vary in space. The objective of this study was to examine the effects of initial soil water content (θ_i) and K_s variability on runoff simulations using the LImburg Soil Erosion Model (LISEM) in a small watershed in the Chinese Loess Plateau, based on model parameters derived from intensive measurements. The results showed that the total discharge (TD) and peak discharge (PD) were underestimated when the variability of θ_i and K_s was partially considered or completely ignored compared with those when the variability was fully considered. Time to peak (TP) was less affected by the spatial variability compared to TD and PD. Except for TP in some cases, significant differences were found in all hydrological variables (TD, PD and TP) between the cases in which spatial variability of θ_i or K_s was fully considered and those in which spatial variability was partially considered or completely ignored. Furthermore, runoff simulations were affected more strongly by K_s variability than by θ_i variability. The degree of spatial variability influences on runoff simulations was related to the rainfall pattern and θ_i. Greater rainfall depth and instantaneous rainfall intensity corresponded to a smaller influence of the spatial variability. Stronger effects of the θ_i variability on runoff simulation were found in wetter soils, while stronger effects of the K_s variability were found in drier soils. For accurate runoff simulation, the θ_i variability can be completely ignored in cases of a 1-h duration storm with a return period greater than 10 years, while K_s variability should be fully considered even in the case of a 1-h duration storm with a return period of 20 years.

Editor D. Koutsoyiannis; Associate editor A. Fiori     

STUDY OF INFILTRATION INTO A HETEROGENEOUS SOIL USING MAGNETIC RESONANCE IMAGING

In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.     

A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.     

Comparing hydraulic properties of different forest floors

The forest floor plays an important role in runoff rate, soil erosion and soil infiltration capacity by protecting mineral soils from the direct impact of falling raindrops. Forest floor consists of different kinds of litter with different hydraulic properties. In this study, the inverse method was used to estimate the hydraulic properties of three kinds of forest floor (broad‐leaved, needle‐leaved and mixed‐stand) at three replications in a completely random design. Forest floor samples were collected from the Gilan Province, Iran. The samples were piled up to make long columns 40.88 cm high with an inner diameter of 18.1 cm. Artificial rainfall experiments were conducted on top of the columns, and free drainage from the bottom of the columns was measured in the laboratory. Saturated hydraulic conductivity (K_s), saturated water content and water retention curve parameters (van Genuchten equation) were estimated by the inverse method. The results showed that the K_s of needle‐leaved samples differed significantly (p < 0.05) from those of broad‐leaved and mixed‐stand samples, whereas the latter two did not differ in this regard. No significant differences emerged in the water retention function parameters of van Genuchten (θ_r, β and α) in the three forest floor samples. The saturated water content of mixed‐stand samples was significantly different (p < 0.05) from that of broad‐leaved and needle‐leaved treatments with the latter two samples showing no significant difference. The good agreement between simulated and observed free drainage for all forest floor samples in the validation period indicates that the estimated hydraulic properties efficiently characterize the unsaturated water flow in the forest floor. Copyright © 2013 John Wiley & Sons, Ltd.