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.     

Temporal and spatial variation of infiltration in urban green infrastructure

Infiltration is the primary mechanism in green stormwater infrastructure (GSI) systems to reduce the runoff volume from urbanized areas. Soil hydraulic conductivity is most important in influencing GSI infiltration rates. Saturated hydraulic conductivity (Ksat) is a critical parameter for GSI design and post-construction performance. However, Ksat measurement in the field is problematic due to temporal and spatial variability and measurement errors. This review paper focuses on a comparison of methods for in-situ Ksat measurement and the causes of temporal and spatial variations of Ksat within GSI systems. Automated infiltration testing methods, such as the Modified Philip–Dunne (MPD) and SATURO infiltrometers, show promise for efficient Ksat measurements. Soil Ksat values can change over time and substantially vary throughout a GSI, which can be attributed to multiple factors, including but not limited to temperature changes, soil composition and properties, soil compaction level, plant root morphology and distribution, biological and macrofauna activities in the soil, inflow sediment characteristics, quality of infiltrating water, and measurement errors. There is evidence that infiltration rates in vegetated urban GSI systems are sustained given an appropriate GSI design, reasonable concentration of suspended sediments in the inflow runoff, and routine maintenance procedures. These observations indicate that clogging can be counteracted by processes that tend to increase the soil hydraulic conductivity (e.g., plant root and biological activities). This self-sustainability underlines that infiltration-based GSI systems are a reliable long-term stormwater management solution. Recommendations on how to incorporate the temporal changes of Ksat in GSI design and on obtaining a spatially-representative Ksat for the GSI design are presented.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

A modified subsurface stormflow model of hillsides in forest catchment

Subsurface flow plays an important role in forest catchment hydrological processes, for which a modified model is established in this paper. Firstly, by taking soil samples from Natural Preserve Forests of Changbai Mountain, two crucial parameters for subsurface flow, the saturated hydraulic conductivity and effective porosity, were measured in the laboratory. Secondly, submodels of the two parameters varying logarithmically with soil depth were established through regressive analysis. Then a modified subsurface stormflow model (Modified model) was founded by substituting the submodels into a storage–discharge model (Sloan's model), established by Sloan in 1983. Finally, to verify the Modified model, five rainfall events on a simulated hillside were carried out. The subsurface flow processes were simulated using the Modified model, Sloan's model and Robinson's model. The comparison of simulated subsurface stormflow processes using the three models respectively with measured ones showed that the Modified model obtained better accuracy for peak flow and total amount of subsurface stormflow than the other two models. Copyright © 2005 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.     

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.     

Tracer and hydrometric study of preferential flow in large undisturbed soil cores from the Georgia Piedmont,USA

We studied the temporal patterns of tracer throughput in the outflow of large (30 cm diameter by 38 cm long) undisturbed cores from the Panola Mountain Research Watershed, Georgia. Tracer breakthrough was affected by soil structure and rainfall intensity. Two rainfall intensities (20 and 40 mm hr⁻¹) for separate Cl⁻ and Br⁻ amended solutions were applied to two cores (one extracted from a hillslope soil and one extracted from a residual clay soil on the ridge). For both low and high rainfall intensity experiments, preferential flow occurred in the clay core, but not in the hillslope core. The preferential flow is attributed to well‐developed interpedal macrochannels that are commonly found in structured clay soils, characteristic of the ridge site. However, each rainfall intensity exceeded the matrix infiltration capacity at the top of the hillslope core, but did not exceed the matrix infiltration capacity at the middle and bottom of the hillslope core and at all levels in the clay core. Localized zones of saturation created when rainfall intensity exceeds the matrix infiltration capacity may cause water and tracer to overflow from the matrix into macrochannels, where preferential flow occurs to depth in otherwise unsaturated soil. Copyright © 1999 John Wiley & Sons, Ltd.