Research Article Volume 2 Issue 3
Department of Civil Engineering, Dr. Subhash Technical Campus, India
Correspondence: Manoj Gundalia, Department of Civil Engineering, Dr. Subhash Technical Campus, India
Received: December 05, 2017  Published: May 4, 2018
Citation: Gundalia M. Estimation of infiltration rate based on complementary error function peak for Ozat watershed in Gujarat (India). Int J Hydro. 2018;2(3):289294. DOI: 10.15406/ijh.2018.02.00083
Infiltration refers to the downward movement of water into the soil from the surface, driven by force of gravity and capillary attraction. Infiltration of water in unsaturated soils has long been an important impact in soil science, hydrology, and geotechnical engineering.Fieldmeasurement of infiltration is very tedious and time consumingand, hence, it is often estimated from different conceptual models. Thus, in the present study, a model based on complementary error function peak (ERFC) is proposed to determine infiltration rate for the Ozat Watershed of Gujarat (India). The two quantitative standard statistical performance evaluation measures, refined Willmott’s index (d_{r}) and mean absolute error (MAE) are employed in comparingand evaluating the performance of theproposed model with existing Horton’s, Modified Kostiakov, and GreenAmpt infiltration models. The results of the models are validated with the Doublering infiltrometer field data of different soil types. The proposed ERFC based model was judged to be more consistent with d_{r}=(0.87 to 0.90) and MAE=(0.47 to 0.27cm/h_{r}). The results presented in this work are quite encouraging. Further the findings indicate that the use of ERFC based proposed model appears to be the most suitable and appropriate for estimation of infiltration rate in Ozat watershed.
Keywords: infiltration, infiltration rate, infiltration models, ERFC, Ozat watershed
Infiltration is the intake of water into the soil profile and it is one of the major components of the hydrologic cycle. According to^{1} Infiltration is the methodology by which water moves descending at the soil surface pass in the soil. Itseparates water into two major hydrologic components  surface runoff and subsurface recharge and it is usually measured by the depth of the water layer that can enter the soil in unit time.Infiltration is a complex process with many factors contributing to the infiltration rate. Water that falls as precipitation may run over land eventually reaching streams, lakes, rivers and oceans or infiltrate through the soil surface, into the soil profile. Water that runs off over land causes erosion, flooding and degradation of water quality. Infiltration, on the other hand, constitutes the sole source of water to sustain the growth of vegetation, is filtered by the soil which removes many contaminants through physical, chemical and biological processes, and replenishes the ground water supply to wells, springs and streams.^{2,3} Infiltration is critical because it supports life on land on our planet. The ability to quantify infiltration is of great importance in many hydrological problems such as runoff estimation, watershed management, soil moisture budgeting, irrigation planning, prediction of flooding, erosion and pollutant transport.
Infiltration rate characterizes by time from onset of rain or irrigation, initial water content of the soil, hydraulic conductivity, surface conditions, and profile depth and layering.^{4} Several wellknown empirical models such as Kostiakov, Horton, and Holtan, and approximate physically based models like those of Green and Ampt and Philip equations are used to estimate infiltration rate for field applications. Unfortunately, single model does not provide consistently better predictions for infiltration in different types of soil. Therefore, it is difficult for researchers to select the most suitable model to estimate infiltration rate. In the present study, a complementary error function peak (ERFC) based model is proposed to determine infiltration rate for the Ozat Watershed of Gujarat (India).The constant infiltration ratesof different soils of the study area were calculated by double ring infiltrometer method, and compared with calculated values from Horton’s, Modified Kostiakov, and GreenAmpt infiltration models along with proposed model. Infiltration is the process by which water enters the soil. It separates water into two major hydrologic components  surface runoff and subsurface recharge. The assessment of runoff risk has assumed an increased importance because of concerns about the associated pollution hazards. Accurate determination of infiltration rates is essential for reliable prediction of surface runoff. As environmental impact assessments are concerned with longterm effects, it is essential that the infiltration data on which they are based should be reasonably stable over decades. For planning purposes it is essential to know the stability of infiltration data for the infiltration capacity of individual soils is adequate to cope with the anticipated hydrologic loads. A high infiltration rate is generally desirable for plant growth and the environment. In some cases, soils that have unrestricted water movement through their profile can contribute to environmental concerns if misapplied nutrients and chemicals reach groundwater and surface water resources via subsurface flow. In India also, very few studies have been reported that focused on infiltration based rainfall simulator experiment.^{5}
A few studies had been carried out^{6‒9} for estimating the infiltration rates in various basins in different parts of India using double ring infiltrometer. The appearance of layered soil is more common than uniform soil in nature with the latter being the exception, vertical infiltration in layered soils has drawn much attention and been studied by many authors. The constant rate approached after a sufficiently large time is referred to as the steadyinfiltration rate. The process is described by the equations of^{10,11} which show a decreasing infiltration rate as a function of time. Many researchers have evaluated the models by comparing the computed and observed infiltration rates.^{12‒15} Most recent reports^{16‒21 }have shown that the ability to quantify infiltration is of great importance in soil management, especially, in irrigation and drainage designs. Hence, a robust infiltration model, that can correctly predict the actual infiltration, can be quite effective in planning and designing of water resources systems.^{19} The main objective of this study was to develop a model for determination of infiltration rate with specific focus on:
Ozat is a river flowing in western India in Gujarat state whose origin is near Visavadar and meets in Arebian Sea. Ozat is third largest river of Saurashtra region after Bhadar and Shetrunji rivers. Ozat watershed considered in this study geographically located within the latitudes 21^{0} 19^{’}N to 21^{0} 33^{’}N and the longitudes 70^{0} 37^{’}E to 70^{0} 59^{’}E respectively. Satellite Imageries: The Indian Remote Sensing satellite with Linear Imaging Self Scanning sensors (IRSLISS III) satellite data of scale 1:50000 were collected from Bhaskaracharya Institute for Space Applications and Geo Informatics (BISAG), Gandhinagar and used to prepare soil map of the study area. Soil taxonomy developed by United States Department of Agriculture and the National Cooperative Soil Survey is widely used system for classifying soils. The prime objective of soil taxonomy is to develop a hierarchical classification that reflects the relationships between different soils, and between soils and the factors responsible for their character. According to^{22} many countries have developed soil classificationsystems for national use, but Soil Taxonomy^{23} is used worldwide. There are 12 recognized soil orders in the world. The three soil orders, Entisols, Inceptisols, and Vertisols, are mapped in the study region.Locations for measurement of infiltration rate were selected based on different types of soil in the study area (Figure 1). The classification of the soils of Ozat watershed based on soil taxonomy is presented in Table 1.The total geographical area of the Ozat watershed is 351.0633 Km^{2} andthe major crops grown are Ground nut, wheat and Cotton. Periodic insufficient rainfall pattern, limited water storage capacity of aquifer and natural water conservation are vital issues for this region. Water availability is a critical factor in this area. Infiltration is the important parameter and significantly assessment of soil moisture, determination of surface runoff, to determination of efficiency of water harvesting structure in ground water recharging, and therefore accurate estimation of infiltration is needed.
Area 
Land without soil 
Entisols 
Inceptisols 
Vertisols 
Total 
Km^{2} 
38.9145 
7.4891 
229.9799 
74.6798 
351.0633 
% 
11.08 
2.13 
65.51 
21.27 
100.00 
Table 1 Spatial variation of soil properties in the Ozatwatershed
The following infiltration models were assessed for finding best fitting model to observed field infiltration rate data.
Horton’s model
The Horton model of infiltration^{11,24} is one of the bestknown models in hydrology. Horton recognized that infiltration capacity (f) decreased with time until it approached a minimum constant rate (f_{c}). He derived his equation for infiltration, which describes a pattern of exponential decay of infiltration rate from this basic relationship. The final form of the Horton equation is expressed as
$\text{f}={\text{f}}_{\text{c}}+\left({\text{f}}_{\text{o}}{\text{f}}_{\text{c}}\right){\text{e}}^{{k}_{h}\text{t}}$ (1)
Where
f is infiltration capacity at any time t.
f_{c} is final steady state infiltration capacity.
f_{o} is initial infiltration capacity
k_{h }Horton’s constant representing rate of decrease in infiltrationcapacity
t is time in hours.
He attributed this decrease in infiltration primarily to factors operating at the soil surface rather than to flow processes within the soil.^{25} Horton’s equation has been widely used because it generally provides a good fit to data. Although the Horton equation is empirical in that k_{h}, f_{c} and f_{o} must be calculated from experimental data, rather than measured in the laboratory, it does reflect the laws and basic equations of soil physics.^{26} However, the Horton equation is cumbersome in practice since it contains three constants that must be evaluated experimentally.^{4} A further limitation is that it is applicable only when rainfall intensity exceeds fc.^{2,27} criticized on its assumptions that hydraulic conductivity is independent of the soil water content. Horton’s approach has also been criticized because he neglects the role of capillary potential gradients in the decline of infiltration capacity over time and attributes control almost entirely on surface conditions. He also discovered that Horton’s perceptual model was far more sophisticated and complete than normally presented in hydrological texts.
Modified kostiakov model
Kostiakov^{10} and independently^{28} proposed a simple empirical infiltration equation based on curve fitting from field data. The major limitation of the KostiakovLewis variant, also known as the "Modified Kostiakov" is its reliance on the zero final intake rates. In most cases the infiltration rate instead approaches a finite steady value, which in some cases may occur after short periods of time.^{29} proposed a modification to Kostiakov’s equation by adding a constant to the equation that represents the final infiltration rate reached when the soil becomes saturated after prolonged infiltration.
$\text{f}={\text{k}}_{\text{k}}{\text{at}}^{\left(\text{a}1\right)}+{\text{f}}_{\text{c}}$ (2)
Where
k_{k} and a are constants whose values depends on soil type.^{30} showed that the ‘k’ and ‘a’ terms of the equation do have physical meaning even though several authors have described it as purely empirical.^{31} however, found empirically that the value of ‘a’ was consistently less than one, he recommended the modified Kostiakov equation for routine modeling of the infiltration process on soils with rapid water intake rates.
Greenamptmodel
Green, Ampt^{32} proposed an approximate model that directly applies Darcy’s law. The original equation was derived for infiltration from a ponded surface into a deep homogeneous soil with uniform initial water content. The Green and Ampt (GA) model has been found to apply best to infiltration into uniform, initially dry, coarse textured soils which exhibit a sharply defined wetting front.^{33} GA applies the physical principles governing infiltration for simplified boundary and initial conditions. It imply ponded surface conditions from time zero on,^{4} and is based on assumptions of uniform movement of water from the surface down through deep homogenous soil with a welldefined wetting front. These assumptions are more valid for sandy soils than for clay soils and reduce the amount of physical soil data needed from that of numerical solutions, but also limit their applicability under changing initial and boundary conditions.^{34} This pattern is often called a piston displacement profile or plug flow.
$\text{f}=\text{b}+\frac{\text{c}}{\text{F}}$ (3)
Where
F is cumulative infiltration.
b and c are GreenAmpt parameters of infiltration.
The GA equation produced reasonably good predictions for nonuniform soil profiles that become denser with depth,^{35} for profiles where hydraulic conductivity decreases with depth^{36} or increases with depth,^{37} and for soils with partially sealed surfaces.^{33,36} described a tabular procedure for calculating the GA relationship between cumulative infiltration and time for soils with nonuniform initial water contents and hydraulic conductivities. GA model parameters can be estimated from soil water properties and do not require measured infiltration data. Thus it should be able to produce estimates at lower cost than empirical equations.
In spite of the wider variety of models no generally accepted method to estimate infiltration rate is applicable to all field conditions or problems. On the other hand, accuracy and reliability of these methods vary widely according to regional soil types.^{38} stated that the various models can estimate different values of the final soil infiltration rate which seems incorrect as the parameter is soildependent. Therefore, it is difficult for many scientists to select the most suitable equation to use for a given study. The empirical equations, such as Kostiakov and Horton equations provide infiltration rates based on measured field data and therefore provide more realistic estimates with calibrated parameters for the same site or site having similar conditions. The model formulation was considered as infiltration rate has a tendency to decrease with the increase of time. The complementary error function is important special function which appears in the solutions of diffusion problems in heat, mass and momentum transfer, probability theory, the theory of errors and various branches of mathematical physics. The error function (erf) is obtained by integrating the normalized Gaussian distribution. The complementary error function (erfc) equals one minus the error function and is defined as:
$\text{erfcx}=1\frac{2}{\sqrt{\text{\pi}}}\underset{\text{x}}{\overset{\infty}{{\displaystyle \int}}}{\text{e}}^{{\text{t}}^{2}}$ (4)
Where
$t=\frac{x}{\sqrt{2\sigma}}$In the present work, an alternative strategy was followed; the ft relationship is expressed in terms of the complementary error function erfc with argument (β*t) as:
$\text{f}=\text{\alpha *ERFC}\left(\text{\beta *t}\right)+{\text{f}}_{\text{c}}$ (5)
Where
α and β are constants
Double ring infiltrometer method was used formeasurement of infiltration rates at selected sites. 30cm diameter inner ring and 60cm diameter outer ring with 25cm deep were driven at about 15cm deep in soil by using falling weight type hammer strikingon a wooden plank placed on top of ring uniformlywithout or undue disturbance to soil surface. Measurements were taken at regular time interval of 5, 10, 15, 20, 25 min. up to getting a constant infiltration rate. The measured steady state infiltration rate and other important properties of the soil at selected sites are shown in Table 2. According to soil taxonomy, the Ozat watershed is characterized by mainly three types of soil, i.e., Entisols, Inceptisols and Vertisols (Table1). Entisols are immature soils with little evidence of soil formation and indicated by “ent”. By texture, they are sandyclay, loam or clayloam to clay. Structurally these soils are weak, mainly subangular, blocky and sometimes crumblike, calcareous and alkaline in nature. Inceptisols are designated by “ept”. They are young soils with weakly developedsubsurface horizons but more developed than Entisols. These soils are texturally siltyloam to clay and neutral to alkaline in reaction. Vertisols are fairly deep, heavy clay soils, and have no definite structure. Because of the montmonllonitic nature of the clay minerals they shrink and crack when dry and expandwhen moist. The soils are saline and texturally sandy loam with silty clay loam. Vertisol are indicated by “ert” and characterized by low noncapillary pore spacewhich prevents drainage of excess water. The parameters of selected models were optimized by using solver tools of Microsoft excel to improve their performance. Optimized values of the parameters of different models for the selected sites of the study areaare presented in Table 3. It is evident from the Table 3 that the values of parameters of infiltration models vary for different soil types.The values of the parameters estimated were then incorporated into the respective models and the capability of each model to simulate infiltration rate for site was evaluated by comparing the models simulated data with field measured data.
Selected Sites 
Soil type 
Date of exp. 
Temp. (^{o}C) 
Dry density (Kg/m^{3}) 
Moist. cont. (%) 
Steady state infilt. rate (cm/hr) 
Lalpur 
Loamy, mixed, hyperthermic, calcareous, lithicustorthents 
20^{th} Feb., 2017 
26 
1109 
15 
2.39 
Kankachiyala 
Fine, montmorrilonitic, hyperthermic, calcareous, verticustochrepts 
1^{st} Feb., 2017 
22 
1606 
3.45 
3.6 
Chhelanka 
Clayey, montmorrilonitic, hyperthermic, lithic, ustochrepts 
25^{th} Mar., 2017 
26 
1572 
17.39 
4.32 
Vajadi 
Fine, montmorrilonitic, hyperthermic, calcareous, typicchromusterts 
12^{th} Jan., 2017 
18 
1089 
14 
8.16 
Table 2 Various soil properties and final rate of infiltration at selected sites
Selected sites/ Soil orders 
Horton’s model 
Modified kostiakov model 
Green – ampt model 
Proposed model 


K_{h} 
K_{k} 
a 
b 
c 
α 
β 
Lalpur (Entisols) 
2 
1748.15 
0.57 
2.44 
3.74 
6.23 
1.49 
Kankachiyala 
1.26 
1921.06 
0.69 
4.28 
12.51 
6.7 
0.44 
(Inceptisols) 

Chhelanka (Inceptisols) 
1.27 
1552.28 
0.73 
4.71 
8.59 
5.97 
0.58 
Vajadi (Vertisols) 
6.22 
3245.36 
0.66 
5.76 
34.49 
64.04 
8.41 
Table 3 The values of optimized parameters of infiltration models for selected sites
Selected sites/ Soil types 
Horton’s model 
Modified kostiakov model 
Green – ampt model 
Proposed model 

MAE (cm/hr) 
dr 
MAE (cm/hr) 
dr 
MAE (cm/hr) 
dr 
MAE (cm/hr) 
dr 

Lalpur (Entisols) 
0.69 
0.63 
0.68 
0.71 
1 
0.47 
0.32 
0.89 
Kankachiyala (Inceptisols) 
0.69 
0.83 
0.48 
0.82 
0.74 
0.65 
0.42 
0.87 
Chhelanka (Inceptisols) 
0.35 
0.88 
0.45 
0.81 
0.86 
0.51 
0.27 
0.9 
Vajadi (Vertisols) 
0.9 
0.75 
1.6 
0.66 
1.37 
0.67 
0.58 
0.86 
Table 4 Performance of selected models and proposed model for selected sites
Models evaluation was carried out with two popular statistical criteria refined Willmott’s index (d_{r}) (Dimensionless statistic) and mean absolute error (MAE) (error index statistic). MAE is the most natural and unambiguous measure of average error magnitude. It provides more robust measure of average model error, since it is not influenced by extreme outliers and estimates model error in the units of the variable.^{39} A higher MAE value indicates poor model performance and vice versa. MAE=0 indicates a perfect fit. The d_{r} is applied to quantify the degree to which values of measured infiltration rate are captured by the selected models. The range of d_{r} is from 1.0 to 1.0. A d_{r} of 1.0 indicates perfect agreement between model and observation and a d_{r} of 1. 0 indicates either lack of agreement between the model and observation or insufficient variation in observations to adequately test the model. The results of the statistical analysis of all the models versus observed measure values for selected sites are presented in Tables 4 respectively. From analysis it was found that for all soil type proposed model was best fitting with lower MAE and higher d_{r} values. The proposed model with a MAE of 0.27cm/hr and d_{r} value 0.90 most closely predicted the measured infiltration rate for Chhelanka, followed by Horton’s model for Chhelanka and Vajadi with MAE values of 0.35cm/hr and 0.90cm/hr and d_{r} values of 0.88 and 0.75respectively. Kostiakov Model provided good results with MAE of 0.68cm/hr and d_{r }value 0.71for Lalpur site. GreenAmpt Model provided less accurate estimates for all sites.Performances of all modelsfor selected sites are shown in Figures 25.
This paper investigates the ability of the novel infiltration model based on based on Complementary Error Function Peak in estimating the infiltration rate from actual field data in comparison to other three popular models in like situations. From the research work the following conclusions are drawn:
None.
The author declares there is no conflict of interest.
©2018 Gundalia. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work noncommercially.