Effects of rainfall intensity fluctuations on infiltration and runoff: rainfall simulation on dryland soils,Fowlers Gap,Australia

Small plots and a dripper rainfall simulator were used to explore the significance of the intensity fluctuations (‘event profile’) within simulated rainfall events on infiltration and runoff from bare, crusted dryland soils. Rainfall was applied at mean rain rates of 10 mm/h. Fourteen simulated rainfall events each involved more than 5000 changes of intensity and included multipeak events with a 25‐mm/h peak of intensity early in the event or late in the event and an event that included a temporary cessation of rain. These are all event profiles commonly seen in natural rain but rarely addressed in rainfall simulation. A rectangular event profile of constant intensity, as commonly used in rainfall simulation experiments, was also adopted for comparative purposes. Results demonstrate that event profile exerts an important effect on infiltration and runoff for these soils and rainfall event profiles. ‘Uniform’ events of unvarying intensity yielded the lowest total runoff, the lowest peak runoff rate and the lowest runoff ratio (0.13). These parameters increased for ‘early peak’ profiles (runoff ratio 0.24) and reached maxima for ‘late peak’ profiles (runoff ratio 0.50). Differences in runoff ratio and peak runoff rate between the ‘uniform’ event profile and those of varying intensity were all statistically significant at p ≤ 0.01. Compared with ‘uniform’ runs, the varying intensity runs yielded larger runoff ratios and peak runoff rates, exceeding those of the ‘uniform’ events by 85%–570%. These results suggest that for small‐plot studies of infiltration and erosion, the continued use of constant rainfall intensity simulations may be sacrificing important information and misrepresenting the mechanisms involved in runoff generation. The implications of these findings for the ecohydrology of the research site, an area of contour‐aligned banded vegetation in which runoff and runon are of critical importance, are highlighted. Copyright © 2011 John Wiley & Sons, Ltd.     

Incorporating rainfall intensity into daily rainfall records for simulating runoff and infiltration into soil profiles

Hydrological modelling often implies the use of rainfall data. Its quality and resolution directly affect the accuracy of the simulation results. This study illustrates that a simple approach of incorporating rainfall intensity information in daily rainfall records significantly improves the simulation of surface runoff and rainfall infiltration into soil profiles. The procedure is developed using a frequency analysis on rainfall data of the Royal Meteorological Institute of Belgium, collected with a resolution of 10 min and for a consecutive period of 61 years. The frequency analysis of the data allowed the incorporation of rainfall intensity information into daily rainfall records. To test the effect of this approach the surface runoff and water flow into three different soil types was simulated using the HYDRUS‐1D model for a typical dry, normal and wet year. The simulation results whereby the observed 10‐min rainfall data was used as input were considered as the reference. Comparative analysis revealed that the simulations using the 10 min rainfall data deducted from the incorporation of rainfall intensity into daily rainfall records, deviate a maximum 1·2% from the reference and produce much better results than the Soil Conservation Service (SCS) runoff curve‐number method because rainfall intensity is considered in the procedure presented. The SCS curve‐number method typical overestimates surface runoff during periods of low rainfall intensity (winter) and underestimate runoff during periods of high rainfall intensities (summer). Copyright © 2002 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.     

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.     

Stemflow production and intrastorm rainfall intensity variation: an experimental analysis using laboratory rainfall simulation

Laboratory rainfall simulation experiments using a small artificial olive tree are used to show that the fraction of rain falling at a constant intensity that becomes stemflow rises from 9% at 2.5 mm/h to 36.1% at 35 mm/h. Natural rainfall events commonly exhibit wide fluctuations of intensity. Simulated rainfall events each having a mean intensity of 10 mm/h, but containing short intensity peaks of 20 to 100 mm/h at varying intra‐event positions, were used to explore the effect of varying intensity profiles. Results demonstrate that changes in rainfall event profile are associated with wide variation in stemflow flux, stemflow volume and stemflow fraction. When applied to an initially dry plant, rainfall events with a late intensity peak yielded an average peak stemflow flux up to 188% larger than events of contrasting profile, such as early peak events. The increase was smaller, up to 141%, when rain was applied to plants that were already partially wet, but was again found in events with a late intensity peak. Moreover, such events yielded a peak stemflow flux up to approximately seven times larger than comparable events of uniform intensity. Likewise, changing event profile with no change in rainfall depth was associated with a maximum stemflow fraction that was 31% larger than theminimum stemflow fraction, and a maximum stemflow volume that was nearly 37% larger than the minimum stemflow volume. These results suggest that rainfall event profile exerts a significant effect on all of the studied stemflow parameters. It is hypothesized that this is a consequence of the way in which intensity profile affects the rate of wetting‐up of trickle pathways on the plant, and variation in the time taken for these pathways to become fully connected. Event profile must therefore be considered along with plant architecture in seeking to understand stemflow. Copyright © 2014 John Wiley & Sons, Ltd.     

The role of rainfall temporal and spatial averaging in seasonal simulations of the terrestrial water balance

The partitioning of rainfall into surface runoff and infiltration influences many other aspects of the hydrologic cycle including evapotranspiration, deep drainage and soil moisture. This partitioning is an instantaneous non-linear process that is strongly dependent on rainfall rate, soil moisture and soil hydraulic properties. Though all rainfall datasets involve some degree of spatial or temporal averaging, it is not understood how this averaging affects simulated partitioning and the land surface water balance across a wide range of soil and climate types. We used a one-dimensional physics-based model of the near-surface unsaturated zone to compare the effects of different rainfall discretization (5-min point-scale; hourly point-scale; hourly 0.125° gridded) on the simulated partitioning of rainfall for many locations across the United States. Coarser temporal resolution rainfall data underpredicted seasonal surface runoff for all soil types except those with very high infiltration capacities (i.e., sand, loamy sand). Soils with intermediate infiltration capacities (i.e., loam, sandy loam) were the most affected, with less than half of the expected surface runoff produced in most soil types when the gridded rainfall dataset was used as input. The impact of averaging on the water balance was less extreme but non-negligible, with the hourly point-scale predictions exhibiting median evapotranspiration, drainage and soil moisture values within 10% of those predicted using the higher resolution 5-min rainfall. Water balance impacts were greater using the gridded hourly dataset, with average underpredictions of ET up to 27% in fine-grained soils. The results suggest that “hyperresolution” modelling at continental to global scales may produce inaccurate predictions if there is not parallel effort to produce higher resolution precipitation inputs or sub-grid precipitation parameterizations.     

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.     

Field studies on the influence of rainfall intensity,vegetation cover and slope length on soil moisture infiltration on typical watersheds of the Loess Plateau,China

Soil moisture is a key process in the hydrological cycle. During ecological restoration of the Loess Plateau, soil moisture status has undergone important changes, and infiltration of soil moisture during precipitation events is a key link affecting water distribution. Our study aims to quantify the effects of vegetation cover, rainfall intensity and slope length on total infiltration and the spatial variation of water flow. Infiltration data from the upper, middle and lower slopes of a bare slope, a natural grassland and an artificial shrub grassland were obtained using a simulated rainfall experiment. The angle of the study slope was 15° and rainfall intensity was set at 60, 90, 120, 150, and 180 mm/hr. The effect these factors have on soil moisture infiltration was quantified using main effect analysis. Our results indicate that the average infiltration depth (ID) of a bare slope, a grassland slope and an artificial shrub grassland slope was 46.7–73.3, 60–80, and 60–93.3 cm, respectively, and average soil moisture storage increment was 3.5–5.7, 5.0–9.4, and 5.7–10.2 mm under different rainfall intensities, respectively. Heavy rainfall intensity and vegetation cover reduced the difference of soil infiltration in the 0–40 cm soil layer, and rainfall intensity increased surface infiltration differences on the bare slope, the grassland slope and the artificial shrub grassland slope. Infiltration was dominated by rainfall intensity, accounting for 63.03–88.92%. As rainfall continued, the contribution of rainfall intensity to infiltration gradually decreased, and the contribution of vegetation cover and slope length to infiltration increased. The interactive contribution was: rainfall intensity * vegetation cover > vegetation cover * slope length > rainfall * slope length. In the grass and shrub grass slopes, lateral flow was found at a depth of 23–37 cm when the slope length was 5–10 m, this being related to the difference in soil infiltration capacity between different soil layers formed by the spatial cross-connection of roots.     

Reliability of rainfall kinetic power–intensity relationships

The rainfall erosivity plays a fundamental role in water soil erosion processes and it can be expressed by its kinetic power. At first in this paper, the raindrop‐size distributions measured, in the period June 2006–March 2014, by an optical disdrometer installed at the Department of Agricultural and Forestry Sciences of University of Palermo are aggregated into rainfall intensity classes, having different ranges, and the measured kinetic power values are determined. Measured kinetic power values are initially used for testing the applicability of the kinetic power‐rainfall intensity relationships proposed by Wischmeier and Smith ( 1978 ), used in Universal Soil Loss Equation (USLE), Brown and Foster ( 1987 ) (RUSLE), and McGregor et al. ( 1995 ) (RUSLE2). Then, the reliability of a theoretical relationship for estimating the kinetic power by rainfall intensity and median volume diameter is verified. Finally, using the literature available datasets, corresponding to measurements carried out by different techniques and in different geographical sites, the analysis demonstrated that the rainfall intensity is not sufficient to determine the rainfall kinetic power. On the contrary, the theoretically deduced relationship allows to reproduce adequately the kinetic power of all available datasets, demonstrating that the knowledge of both rainfall intensity and median volume diameter allows a reliable estimate of the rainfall erosivity.     

Automated rainfall simulator for variable rainfall on urban green areas

Rainfall simulators can enhance our understanding of the hydrologic processes affecting the total runoff to urban drainage systems. This knowledge can be used to improve urban drainage designs. In this study, a rainfall simulator is developed to simulate rainfall on urban green surfaces. The rainfall simulator is controlled by a microcomputer programmed to replicate the temporal variations in rainfall intensity of both historical and synthetic rainfall events with constant rainfall intensity on an area of 1 m². The performance of the rainfall simulator is tested under laboratory conditions with regard to spatial uniformity of the rainfall, the kinetic energy of the raindrops, and the ability to replicate historical and synthetic rainfall events with temporally varying intensity. The rainfall simulator is applied in the field to evaluate its functionality under field conditions and the influence of wind on simulated rainfall. Finally, a field study is carried out on the relationship between runoff, soil volumetric water content, and surface slope. Performance and field tests show that the simulated rainfall has a uniform spatial distribution, whereas the kinetic energy of the raindrops is slightly higher than that of other comparable rainfall simulators. The rainfall simulator performs best in low wind speed conditions. The simulator performs well in replicating historical and synthetic rainfall events by matching both intensity variations and accumulated rainfall depth. The field study shows good correlation between rainfall, runoff, infiltration, soil water content, and surface slope.     

Modeling scale-dependent runoff generation in a small semi-arid watershed accounting for rainfall intensity and water depth

Observed scale effects of runoff on hillslopes and small watersheds derive from complex interactions of time-varying rainfall rates with runoff, infiltration and macro- and microtopographic structures. A little studied aspect of scale effects is the concept of water depth-dependent infiltration. For semi-arid rangeland it has been demonstrated that mounds underneath shrubs have a high infiltrability and lower lying compacted or stony inter-shrub areas have a lower infiltrability. It is hypothesized that runoff accumulation further downslope leads to increased water depth, inundating high infiltrability areas, which increases the area-averaged infiltration rate. A model was developed that combines the concepts of water depth-dependent infiltration, partial contributing area under variable rainfall intensity, and the Green–Ampt theory for point-scale infiltration. The model was applied to rainfall simulation data and natural rainfall–runoff data from a small sub-watershed (0.4 ha) of the Walnut Gulch Experimental Watershed in the semi-arid US Southwest. Its performance to reproduce observed hydrographs was compared to that of a conventional Green–Ampt model assuming complete inundation sheet flow, with runon infiltration, which is infiltration of runoff onto pervious downstream areas. Parameters were derived from rainfall simulations and from watershed-scale calibration directly from the rainfall–runoff events. The performance of the water depth-dependent model was better than that of the conventional model on the scale of a rainfall simulator plot, but on the scale of a small watershed the performance of both model types was similar. We believe that the proposed model contributes to a less scale-dependent way of modeling runoff and erosion on the hillslope-scale.     

A 2D hydrodynamic model for shallow water flows with significant infiltration losses

Infiltration losses may be significant and warrant proper incorporation into mathematical models for river floods in arid and semi-arid areas, rainfall-induced surface runoffs in watersheds and swashes on beaches. Here, a depth-averaged two-dimensional hydrodynamic model is presented for such processes based on the cell-centred finite volume method on unstructured meshes, with the full Green-Ampt equation evaluating the infiltration rate. A local time stepping strategy is employed along with thread parallelization with Open Multi-processing and high-performance computing to reduce model run time and therefore facilitate applications for large-scale processes. The numerical solutions generally agree with the experimental and field-measured data for typical cases with significant infiltration losses. The case study shows that neglecting infiltration leads to an overestimated discharge hydrograph, which cannot be compensated by means of varied bed resistance as estimated by Manning roughness, and the infiltration parameters play disparate roles in modifying shallow flows compared with Manning roughness. In addition, infiltration affects bed shear stress, which in turn modifies the critical bed sediment size that could be initiated for incipient motion by the flow and therefore needs to be properly accounted for when sediment transport and morphological evolution are to be resolved.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial variability of rainfall on a sub‐kilometre scale

The variability of rainfall in space and time is an essential driver of many processes in nature but little is known about its extent on the sub‐kilometre scale, despite many agricultural and environmental experiments on this scale. A network of 13 tipping‐bucket rain gauges was operated on a 1·4 km² test site in southern Germany for four years to quantify spatial trends in rainfall depth, intensity, erosivity, and predicted runoff. The random measuring error ranged from 10% to 0·1% in case of 1 mm and 100 mm rainfall, respectively. The wind effects could be well described by the mean slope of the horizon at the stations. Except for one station, which was excluded from further analysis, the relative differences due to wind were in maximum ±5%. Gradients in rainfall depth representing the 1‐km² scale derived by linear regressions were much larger and ranged from 1·0 to 15·7 mm km^?1 with a mean of 4·2 mm km^?1 (median 3·3 mm km^?1). They mainly developed during short bursts of rain and thus gradients were even larger for rain intensities and caused a variation in rain erosivity of up to 255% for an individual event. The trends did not have a single primary direction and thus level out on the long term, but for short‐time periods or for single events the assumption of spatially uniform rainfall is invalid on the sub‐kilometre scale. The strength of the spatial trend increased with rain intensity. This has important implications for any hydrological or geomorphologic process sensitive to maximum rain intensities, especially when focusing on large, rare events. These sub‐kilometre scale differences are hence highly relevant for environmental processes acting on short‐time scales like flooding or erosion. They should be considered during establishing, validating and application of any event‐based runoff or erosion model. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Hydrology of stockpiled industrial coal exposed to rainfall

The movement of water and particulate matter through a simulated coal pile using a rainfall–coal‐pile simulator is investigated. In conjunction with coal‐pile management practices, this study indicates that coal‐pile systems change their internal and surficial structure with successive rainstorm events. This change in coal‐pile structure alters the proportions, rates and suspended solid concentrations within surface runoff and infiltration during subsequent storm events. Results demonstrate the need to determine the optimal coal‐pile slope that minimizes water contact time, but at the same time minimizes the potential for mass wasting and rill erosion. Promotion of a compacted surface is beneficial for the minimization of solids removal, and a stormwater detention pond with potential flocculant additive is suggested as a viable remedial measure for the containment and treatment of coal‐pile runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

Development of a portable field rainfall simulator for use in hillside vineyard runoff and erosion studies

An inexpensive, mobile field rainfall simulator and runoff plot frame were developed for use on hillside vineyards. The simulator framework and components were lightweight, readily available and easily manageable such that they can be handled by one person during transport, set–up and operation. The vineyard rainfall simulator was simpler than many of the machines in recent use for similar studies, yet offered equal or improved performance for small‐plot studies. The system developed consistent sized 2·58 mm raindrops at intensities ranging from 20 to 90 mm/h. The average distribution uniformity coefficient at an intensity of 60 mm/h was 91·7%, with a deviation of only 2·2%. This coefficient was similar to the range reported for a more complex rotating disk simulator, and was notably greater than that obtained for other similar devices. The system water capacity of 40 l allowed for 1‐h storm durations at 60 mm/h, usually sufficient time for commencement of erosion and runoff. The runoff plot frame was designed to be quickly installed, and to discourage sediment deposition in the routing of runoff to collect containers. Copyright © 2000 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.     

Effects of water quality on infiltration,runoff and interrill erosion processes during simulated rainfall

This paper discusses the effects of water quality on the hydrological and erosion response of non‐saline, non‐sodic soils during simulated rain experiments. It is well known that rain water quality affects the behaviour of saline soils. In particular, rain simulation experiments cannot be run using tap water if realistic values of infiltration rates and soil erosion are to be found. This paper reports on similar effects for non‐saline, non‐sodic soils. Two soils – a well‐aggregated clay‐rich soil developed on marine silty clay deposits and a soil developed on silt loam – were selected and subjected to a series of simulated rainstorms using demineralized water and tap water. The experiments were conducted in two different laboratories in order to obtain results independent of the tap water quality or the rainfall simulator characteristics. The results indicate that time‐to‐ponding is largely delayed by solute‐rich water (tap water). When tap water is used, infiltration rates are significantly overestimated, i.e. by more than 100 per cent. Interrill erosion rates increase by a factor of 2·5–3 when demineralized water is used. The silty clay soil was more affected by the water quality than the silt loam soil, with respect to infiltration and runoff production. Regarding interrill erosion rates, the two tested soils were similarly affected by the water quality. Therefore, it can be concluded that rainfall simulation experiments with non‐dispersive soils (e.g. non‐saline, non‐sodic) must also be conducted using water with very low electrical conductivity (i.e. less than 30–50 µS cm⁻¹), close to that of distilled water. The use of tap water certainly hampers comparisons and the relative ranking of the hydrological and erosion response of different soils, while parameter values, such as final infiltration rate or time‐to‐ponding, cannot be extrapolated and extended to natural situations. Therefore, the majority of hydrological and erosion models and parameter values measured during rainfall simulations in the past should be used with caution for all types of soils. Copyright © 2001 John Wiley & Sons, Ltd.