Different methods for modelling the areal infiltration of a grass field under heavy precipitation

The areal infiltration behaviour of a grass field is studied using a data set of 78 sprinkler infiltration experiments. The analysis of the experimental data shows a distinct event dependency: once runoff begins, the final infiltration rate increases with increasing rainfall intensity. This behaviour is attributed to the effects of small‐scale variability. Increasing rainfall intensity increases the ponded area and therefore the portion of the plot which infiltrates at maximum rate. To describe the areal infiltration behaviour of the grass field the study uses two different model structures and investigates different approaches for consideration of subgrid variability. It is found that the effective parameter approach is not suited for this purpose. A good representation of the observed behaviour is obtained by using a distribution function approach or a parameterization approach. However, it is not clear how the parameters can be derived for these two approaches without a large measurement campaign. The data analysis and the simulations show the great importance of considering the effects of spatial variability for the infiltration process. This may be significant even at a small scale for a comparatively homogeneous area. The consideration of heterogeneity seems to be more important than the choice of the model type. Furthermore, similar results may be obtained with different modelling approaches. Even the relatively detailed data set does not seem to permit a clear model choice. In view of these results it is questionable to use very complex and detailed simulation models given the approximate nature of the problem. Although the principle processes may be well understood there is a lack of models that represent these processes and, more importantly, there is a lack of techniques to measure and parameterize them. Copyright © 2002 John Wiley & Sons, Ltd.     

Robert E. Horton's perceptual model of infiltration processes

Robert E. Horton is best known as the originator of the infiltration excess overland flow concept for storm hydrograph analysis and prediction, which, in conjunction with the unit hydrograph concept, provided the foundation for engineering hydrology for several decades. Although these concepts, at least in their simplest form, have been largely superseded, a study of Horton's archived scientific papers reveals that his perceptual model of infiltration processes and appreciation of scale problems in modelling were far more sophisticated and complete than normally presented in hydrological texts. His understanding of surface controls on infiltration remain relevant today. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of infiltration models with regard to design of rectangular infiltration trenches

ABSTRACT

Guidelines for the design of infiltration trenches of rectangular cross-section are compared. Four concepts of infiltration into a native soil from linear soakaways of rectangular cross-section are discussed. The results of calculations using analytical equations are compared with numerical two-dimensional (2D) simulations using HYDRUS computer code and field measurements. Satisfactory agreement is achieved for Kozeny’s concept of infiltration and the free-surface approach, applied to a trench in a highly saturated sandy loam. Theoretical solutions, neglecting matric potential, suggest that the mean infiltration rate along the wetted perimeter of a rectangular trench varies in the range 1.0 to 1.5 of the value of saturated hydraulic conductivity.     

Simplified modelling of areal average infiltration at the hillslope scale

Two models for estimating expected areal‐average infiltration rate, ī, at the hillslope scale are presented. The first relies upon the condition of a negligible infiltration of surface water running downslope (run‐on process) into a previous heterogeneous soil. It is an adapted version of an earlier semi‐analytical model. The second incorporates the run‐on process and is based on a lumped approach that uses an effective saturated hydraulic conductivity. This latter was parameterized in terms of the main characteristics of rainfall and soil. Both the models were tested by comparison with the results carried out by Monte‐Carlo simulations over different soil types. It was found that the first model simulated ī with maximum errors in magnitude typically less than 10%. The second model provided similar errors in the total volume of overland flow, and the rising limb of the hydrograph experienced a distortion. Lastly, satisfactory results were obtained by comparing the model without run‐on with an empirical approach particularly accurate for fine‐textured soils. Copyright © 2002 John Wiley & Sons, Ltd.     

A combined rainfall infiltration model based on Green–Ampt and SCS‐curve number

The Green–Ampt infiltration equation is an incomplete governing equation for rainfall infiltration due to the absence of an inertia term. The estimation of the capillary pressure head at the wetting front is difficult to determine. Thus, a major limitation of the Green–Ampt model is the constant, non‐zero surface ponding depth. This paper proposes an integrated rainfall infiltration model based on the Green–Ampt model and the SCS‐CN model. It achieves a complete governing equation for rainfall infiltration by momentum balance and the water budget based on the Green–Ampt assumption, and uses the curve number from the SCS‐CN method to calculate the initial abstraction, which is used as a basic parameter for the governing equation of the intensity of rainfall loss during the runoff period. The integrated rainfall infiltration model resolves the dilemma for capillary pressure head estimation, overcomes the limitation of constant, non‐zero surface ponding depth, and facilitates the calculation of runoff for individual flood simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Modified models for better prediction of infiltration rates in trapezoidal permeable stormwater channels

ABSTRACT

In stormwater management, it is important to accurately quantify the infiltration rates to solve urban runoff-related problems. This study proposes a method to improve estimates of the infiltration rate in permeable stormwater channels. As part of the analysis, five infiltration models were evaluated: the Kostiakov, Horton, modified Kostiakov, Philip and SCS (Soil Conservation Service) models. Infiltration tests with various initial water levels were performed on channel models with differing base width and side slopes. The results show that the addition of three parameters that describe the trapezoidal cross-sectional area, i.e. the depth, side slope and base width, in the infiltration models yielded better estimates of the infiltration rate. A comparison of the infiltration capacity values obtained from the models after the three parameters were added with those that were experimentally obtained, shows that the improved modified Kostiakov model is the most suitable model to predict infiltration rates in trapezoidal permeable stormwater channels.     

Soil infiltration processes of different underlying surfaces in the permafrost region on the Tibetan Plateau

ABSTRACT

Soil infiltration processes were evaluated under field conditions by double-ring infiltrometers with different underlying surfaces in permafrost regions of the Tibetan Plateau. The results show that initial infiltration rates, stable soil infiltration rates and cumulative soil infiltration are strongly dependent on the underlying surface types, with the highest initial and stable soil infiltration rates in the alpine desert steppe, and the lowest in alpine meadow. The effects of soil moisture and texture on infiltration processes were also assessed. Within the same underlying surfaces, the values of infiltration parameters increased with the amount of vegetation cover, while soil moisture and soil infiltration rates displayed opposing trends, with fitting slopes of ?0.03 and ?0.01 for the initial and stable soil infiltration rates, respectively. The accuracies of the five models in simulating soil infiltration rates and seven models in predicting cumulative infiltration rates were evaluated against data generated from field experiments at four sites. Based on a comparative analysis, the Horton model provided the most complete understanding of the underlying surface effects on soil infiltration processes. Altogether, these findings show that different underlying surfaces can alter soil infiltration processes. This study provides a useful reference for understanding the parameterization of land surface processes for simulating changes in hydrological processes under global warming conditions in the permafrost region on the Tibetan Plateau.     

Effects of long-term wastewater irrigation on soil physical properties and performance of selected infiltration models in a semi-arid region

ABSTRACT

In many arid and semi-arid countries, wastewater irrigation is becoming a common practice in agriculture. In this study, the effect of long-term (40 years) wastewater irrigation on selected physical and hydraulic properties of soil in different parts of a landscape was investigated. The performance of some infiltration models, including Philip (Ph), Kostiakov (Kos), Kostiakov-Lewis (Kos-L), Horton (Ho), Huggins and Monke (Hug-M), and linear and nonlinear Smith-Parlange (S-P(L) and S-P(NL)), was compared. This study was performed in the Urmia region, Iran, where flooding wastewater irrigation has been practised for at least 40 years. Five paired sites, each of which contained a measurement location at the wastewater-irrigated (WWI) and adjacent control area were studied. Accuracy of the infiltration models was evaluated using several statistical criteria, including root mean square error (RMSE) and Akaike information criterion (AIC). The models were classified into groups using cluster analysis based on level of similarity in their performance. The cumulative water infiltration into soils after 1 h (I_1h) was calculated using the selected most accurate models and introduced so as to use only one term to compare the infiltration behaviour of soils. Based on RMSE and AIC, the performance of the Ph, Ho, Kos and Kos-L models was considerably better than that of Hug-M, S-P(L) and S-P(NL). The ranking of the models in terms of their AIC values was: Kos-L > Ho > Kos > Ph > S-P(L) > Hug-M > S-P(NL). The models were classified into two distinct groups. The similarity among Ph, Ho, Kos and Kos-L models was more than 80% and for Hug-M, S-P(L), and S-P(NL) models, it was more than 79%. However, the similarity between these two groups of models was less than 58%.

Editor M.C. Acreman; Associate editor not assigned     

Comparison of stormflow responses of surface‐mined and forested watersheds in the Appalachian Mountains,USA

The results of a hydrological analysis that was conducted as part of a larger, multifaceted, collaborative effort to quantify ecosystem functions in watersheds subjected to land‐use and land‐cover change are presented. The primary goal of the study was to determine whether a small watershed in the Appalachian region (USA) that was recently subjected to surface mining and reclamation practices produces stormflow responses to rain events that are different from those produced by a nearby reference watershed covered by young, second‐growth forest. Water balances indicated that runoff yields did not vary significantly between the two watersheds on an annual basis. Statistically significant differences (p?0·05) in runoff responses were observed on an event basis, however, with the mined/reclaimed watershed producing, on average (a) higher storm runoff coefficients (2·5×), (b) greater total storm runoff (3×), and (c) higher peak hourly runoff rates (2×) when compared with the reference watershed. Results of a unit hydrograph analysis also showed, unexpectedly, that the modelled unit responses of the two watersheds to effective rainfall pulses were similar, despite the noted differences in land cover. Differences in stormflow responses were thus largely explained by dramatic reductions in cumulative rates of rainfall abstraction (measured using infiltrometers) attributable to soil compaction during land reclamation. Additional field hydrological measurements on other mined watersheds will be needed to generalize our results, as well as to understand and predict the cumulative hydrological impacts of widespread surface mining in larger watersheds and river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

Performance of a conceptual and physically based model in simulating the response of a semi‐urbanized watershed in San Antonio,Texas

The need for accurate hydrologic analysis and rainfall–runoff modelling tools has been rapidly increasing because of the growing complexity of operational hydrologic and hydraulic problems associated with population growth, rapid urbanization and expansion of agricultural activities. Given the recent advances in remote sensing of physiographic features and the availability of near real‐time precipitation products, rainfall–runoff models are expected to predict runoff more accurately. In this study, we compare the performance and implementation requirements of two rainfall–runoff models for a semi‐urbanized watershed. One is a semi‐distributed conceptual model, the Hydrologic Engineering Center‐Hydrologic Modelling System (HEC‐HMS). The other is a physically based, distributed‐parameter hydrologic model, the Gridded Surface Subsurface Hydrologic Analysis (GSSHA). Four flood events that took place on the Leon Creek watershed, a sub‐watershed of the San Antonio River basin in Texas, were used in this study. The two models were driven by the Multisensor Precipitation Estimator radar products. One event (in 2007) was used for HEC‐HMS and GSSHA calibrations. Two events (in 2004 and 2007) were used for further calibration of HEC‐HMS. Three events (in 2002, 2004 and 2010) were used for model validation. In general, the physically based, distributed‐parameter model performed better than the conceptual model and required less calibration. The two models were prepared with the same minimum required input data, and the effort required to build the two models did not differ substantially. Copyright © 2012 John Wiley & Sons, Ltd.     

An evaluation of approaches for modelling hydrological processes in high‐elevation,glacierized Andean watersheds

We use two hydrological models of varying complexity to study the Juncal River Basin in the Central Andes of Chile with the aim to understand the degree of conceptualization and the spatial structure that are needed to model present and future streamflows. We use a conceptual semi‐distributed model based on elevation bands [Water Evaluation and Planning (WEAP)], frequently used for water management, and a physically oriented, fully distributed model [Topographic Kinematic Wave Approximation and Integration ETH Zurich (TOPKAPI‐ETH)] developed for research purposes mainly. We evaluate the ability of the two models to reproduce the key hydrological processes in the basin with emphasis on snow accumulation and melt, streamflow and the relationships between internal processes. Both models are capable of reproducing observed runoff and the evolution of Moderate‐resolution Imaging Spectroradiometer snow cover adequately. In spite of WEAP's simple and conceptual approach for modelling snowmelt and its lack of glacier representation and snow gravitational redistribution as well as a proper routing algorithm, this model can reproduce historical data with a similar goodness of fit as the more complex TOPKAPI‐ETH. We show that the performance of both models can be improved by using measured precipitation gradients of higher temporal resolution. In contrast to the good performance of the conceptual model for the present climate, however, we demonstrate that the simplifications in WEAP lead to error compensation, which results in different predictions in simulated melt and runoff for a potentially warmer future climate. TOPKAPI‐ETH, using a more physical representation of processes, depends less on calibration and thus is less subject to a compensation of errors through different model components. Our results show that data obtained locally in ad hoc short‐term field campaigns are needed to complement data extrapolated from long‐term records for simulating changes in the water cycle of high‐elevation catchments but that these data can only be efficiently used by a model applying a spatially distributed physical representation of hydrological processes. Copyright © 2013 John Wiley & Sons, Ltd.     

A probabilistic investigation of infiltration in the vadose zone: proposal for a new formula of infiltration rate

The infiltration rate in the unsaturated zone is analysed from a probabilistic point of view. It is shown that the empirical formulas of Horton and Kostiakov, without apparent physical basis, are explained in a probabilistic approach. Horton's formula reflects a Markovian process contrary to Kostiakov's formula. This approach made it possible to explain why Kostiakov's formula is more powerful than that of Horton. A new equation of infiltration is proposed. The three formulas were compared, for four types of soil, with the model of van Genuchten based on the Richards equation. Copyright © 2006 John Wiley & Sons, Ltd.     

Spatial aggregation of time‐variant stream water ages in urbanizing catchments

We calibrated an integrated flow–tracer model to simulate spatially distributed isotope time series in stream water in a 7.9‐km² catchment with an urban area of 13%. The model used flux tracking to estimate the time‐varying age of stream water at the outlet and both urbanized (1.7 km²) and non‐urban (4.5 km²) sub‐catchments over a 2.5‐year period. This included extended wet and dry spells where precipitation equated to >10‐year return periods. Modelling indicated that stream water draining the most urbanized tributary was youngest with a mean transit time (MTT) of 171 days compared with 456 days in the non‐urban tributary. For the larger catchment, the MTT was 280 days. Here, the response of urban contributing areas dominated smaller and more moderate runoff events, but rural contributions dominated during the wettest periods, giving a bi‐modal distribution of water ages. Whilst the approach needs refining for sub‐daily time steps, it provides a basis for projecting the effects of urbanization on stream water transit times and their spatial aggregation. This offers a novel approach for understanding the cumulative impacts of urbanization on stream water quantity and quality, which can contribute to more sustainable management. Copyright © 2015 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 new saturated/unsaturated model for stormwater infiltration systems

Infiltration systems are widely used as an effective urban stormwater control measure. Most design methods and models roughly approximate the complex physical flow processes in these systems using empirical equations and fixed infiltration rates to calculate emptying times from full. Sophisticated variably saturated flow models are available, but rarely applied owing to their complexity. This paper describes the development and testing of an integrated one‐dimensional model of flow through the porous storage of a typical infiltration system and surrounding soils. The model accounts for the depth in the storage, surrounding soil moisture conditions and the interaction between the storage and surrounding soil. It is a front‐tracking model that innovatively combines a soil‐moisture‐based solution of Richard's equation for unsaturated flow with piston flow through a saturated zone as well as a reservoir equation for flow through a porous storage. This allows the use of a simple non‐iterative numerical solution that can handle ponded infiltration into dry soils. The model is more rigorous than approximate stormwater infiltration system models and could therefore be valuable in everyday practice. A range of test cases commonly used to test soil water flow models for infiltration in unsaturated conditions, drainage from saturation and infiltration under ponded conditions were used to test the model along with an experiment with variable depth in a porous storage over saturated conditions. Results show that the model produces a good fit to the observed data, analytical solutions and Hydrus. Copyright © 2008 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.