Estimation of Infiltration Models’ Parameters Using Regression Analysis in Irrigation Fields of Northern Ghana

Abdallah Salifu¹, Felix. K. Abagale^2*, Gordana Kranjac-Berisavljevic²

¹Department of Agricultural Mechanisation and Irrigation Technology, Faculty of Agriculture, Food and Consumer Sciences, University for Development Studies, Tamale, Ghana.
²West African Centre for Water, Irrigation and Sustainable Agriculture, University for Development Studies, Tamale, Ghana.

DOI: 10.4236/ojss.2021.113009                

Abstract

The quantification of soil infiltration is necessary for the estimation of water accessibility in soils for plant growth and development. Field infiltration tests runs were conducted on agricultural soils in three irrigation sites of Northern Ghana. The field data were fitted into Green and Ampt, Kostiakov, Philip, Holtan, Soil Conservation Service and Horton infiltration models for the determination of the unknown model parameters. Regression analysis at 95% confidence level using GraphPad Prism 8. Laboratory and field data on infiltration were used for the model fitting and the unknown parameters were determined using the calibrated models. The k and n parameters of Kostiakov model at Bontanga study site were found to be k = 28.0027 and n = 0.5902, k = 17.0294 and n = 0.4504 for Golinga and k = 23.0356 and n = 0.6339 for Libga. For all the models, the coefficient of determination ranged from 0.7612 to 0.9876 indicating strong relationships. Only Holtan model gave different values at all the three study sites. The parameter GIa and i_c of the Holtan model drawn from hydrologic soil group were the same at all the study sites because of the same vegetative cover and surface conditions. The study observed the values of the parameters to have influenced the models’ performance.

Keywords

Parameter, Infiltration Model, Regression, Sorptivity, Transmissivity

1. Introduction

It is necessary to quantify infiltration in order to estimate the accessibility of water for plant growth and to determine the volume of extra water required for irrigation [1]. For this study, field infiltration data were fitted into some models for the determination of unknown model parameters. Six (6) infiltration models were used in the parameter estimation using data from the field and in the laboratory. These models are Green and Ampt (1911), Kostiakov (1932), Horton (1940), Philip (1954), Holtan (1961) and USDA (1972).

The infiltration models that can properly predict infiltration rates are paramount for determination of behaviour of water in soil and water management activities. The challenge in predicting infiltration is mainly because of the variation of infiltration-related soil physical properties from a different site in the field indirectly influencing the model’s parameters. In order to simplify infiltration predictions, several researchers [2] and [7] introduced simple time-dependent algebraic equations.

Numerous models that make the concept of infiltration an easy process have been established for field use. The data from the field experiment are usually processed and fitted into the models to determine numerical values of the model’s parameters and to appraise the applicability of the infiltration models in the soil or site. These parameters are influenced by soil texture, moisture content, soil bulk density, and other soil properties and are site-specific [8].

Truly based approximation equations use obtainable parameters which acquisition from soil water properties is likely and do not need estimated infiltration information.

Different equations are somewhat or completely exact and boundaries should be gotten from estimated infiltration information or generally assessed by different methods. Empirical equations, for example, Kostiakov and Horton are not much prohibitive to the method of water usage since they do not need the suppositions with respect to soil surface and profile conditions that the truly based conditions need [9].

Where soils are diverse, and features, for example, large pore size stream and ensnared air entangle the water penetration processes, empirical equations might possibly give better predictions, once they are utilized by considering comparative settings aided in their establishment. Fundamentally, the grounds that their underlying boundaries are resolved dependent on real field-estimated infiltration information [10] [11] [12].

All the infiltration models use soil boundaries, for example, sorptivity/initial infiltration, pressure driven conductivity, and volumetric water content. They likewise utilize the soil water retention-curve; which esteems are characterized by alignment fit or estimations. These boundaries have vulnerability, brought about by mistakes in information assortment from field or research facility, imprecision or mistaken alignment hardware and spatial and transient varieties.

Numerous penetration models have been created to foresee infiltration rate as a component of time and physical or empirical boundaries, for an alternate sort of soils. Every one of these models is appropriate to specific parameters and the parameters feeding these models have operational challenges. These challenges become significant relying upon the affectability of them to the rate of infiltration [13].

Kostiakov’s infiltration model’s parameters such as the values of empirical constant “c” were found in the range of 0.140 - 0.290 and the values of infiltration decay constants “β” were attained in the range of 0.307 - 0.433. The values of infiltration decay constant “β” were in harmony with the theory of infiltration that defined the values to be positive and always less than unity. It has also been reported that most values of these parameters lie between 0.2 - 0.9. The Philip’s infiltration model parameters such as sorptivity “S” values were found in the range of 0.167 - 0.288 cm/min and conductivity constant “A” values were obtained as −0.001 to −0.009 cm/min for all the six points. [14] Described the high variability of these parameters along the top sequence. Similarly, Horton’s model’s parameters such as constant “β” was determined and in the range of −1.619 to −1.238. The estimated values of these infiltration models’ parameters can be used for developing the infiltration equation for the study area after verification [8].

The study estimated parameters of six (6) soil infiltration models that are commonly used in the prediction of infiltration of agricultural soils.

2. Materials and Methods

2.1. Description of the Study Area

The study was conducted at Bontanga, Golinga and Libga irrigation sites in the Northern Region of Ghana.

Bontanga Irrigation Scheme is located at Bontanga in the Kumbungu District, 34 km Northwest of Tamale. It lies between latitude 9˚57'N and longitude 1˚02'W.

The Golinga Irrigation Scheme is situated in the Tolon District of the Northern Region of Ghana and it lies on latitude 9˚4'N and longitude 1˚0'W.

The Libga Irrigation Scheme is located in the Savelugu Municipality in the Northern Region of Ghana. It lies between latitude 9˚59'N and longitude 0˚85'W [15].

Crops usually cultivated at the irrigation schemes are Rice (Oryza sativa), Onion (Allium cepa), Roselle (Hibiscus sabdariffa), Maize (Zea mays), Okra (Hibiscus esculentus), Tomato (Lycopersicon esculentum) Pepper (Capsicum frutescens), Vegetable jute (Corchorus olitorius), Garden eggs (Solanum melongena), Groundnut (Arachis hypogea), Cowpea (Vigna unguiculata), Amaranth, “Alefu” (Amaranthus spinosus) [16].

2.2. Data Sources

Data was obtained from field and laboratory work. For data collection of the study, undisturbed samples of soil were collected from the field to determine the hydraulic conductivity and, gravimetric water content in the laboratory whilst soil infiltration test was carried out using double ring infiltrometer in the field.

2.3. Model Calibration and Parameter Determination

Calibration or adjustment is the iterative process of contrasting the model and genuine framework, changing the model if essential, re-examining it, until a model is acknowledged or accepted.

The need to estimate the models’ parameters was to create the basis upon which the sensitivity of the models’ parameters and their predictions of the infiltration rates could be determined and to also ease the comparisons between predicted and observed infiltration rates.

The infiltration curves from Libga, Bontanga and Golinga irrigation sites, were used to adjust the empirical Kostiakov, SCS and Horton Models. Infiltration curve parameters for Kostiakov, Horton and SCS models were obtained from Equations (2), (4) and (16) using logarithmic and semi-logarithmic methods.

Philip and Green-Ampt models’ parameters were found using regression analysis. Parameters for the Holtan equation were gotten from estimated starting and saturated moisture contents and from organized standards conforming to the soil state and the nature of vegetation characterizing the study locations [12].

2.4. Kostiakov Equation

The two (2) empirical constants k and n of the Kostiakov model were obtained using Equation (2).

The gradient and intercept of the linear equation were found with a plot of logI verse logt [17] and logk of the equation was the intercept and the slope was n.

$I=k{t}^n$ (1)

That is;

$\log I=\log k-n\log t$ (2)

If the value of n is higher, the slope is steeper and the larger the rate of decay of infiltration. The bigger the value of k, the bigger the initial infiltration value [12] [18].

2.5. Horton Equation

Horton [4] equation for infiltration rate (i) is:

$i=i_c+\left(i_o-i_c\right){\text{e}}^{-\beta t}$ (3)

where: i_c—the continuous infiltration capacity as time gets to infinity; i_o—infiltration capacity at the beginning of water penetration; and β—the positive persistent dependent on soil and introductory conditions of water.

The parameter i_c was found from the infiltration curve whilst i_o was deduced from the experimental values for i and the common log of the outcome plotted as time dependent. β was determined from the gradient of the line and i_o from the intercept. The base value for parameter β was equivalent to the negative gradient, and the start of infiltration capacity, i_o was determined using Equation (4).

$\ln \left(i-i_c\right)=\ln \left(i_o-i_c\right)-\beta t$ (4)

Parameters i₀, i_c and k were obtained by plotting $\ln \left(i-i_c\right)$ against t.

2.6. Holtan Equation

The Holtan equation parameters are; GI, a, and i_c. Based on the site location, the average values for a, were gotten from hydrologic group [12] [19] as presented in Table 1. The values of the steady infiltration (i_c) were obtained from soil hydrologic group built on values founded by [20]. GI was estimated by considering the age-long existence of ground cover and with available storage determined from the moisture deficit using Equation (5);

$M_i={\theta }_s-{\theta }_i$ (5)

The average moisture deficit (M_i) of the control depth, multiplied by the assessed depth gave available storage. Holtan and Lopez [21] indicated that after some number of adjustments, the resulting equation is Equation (6);

$i_p=GIaS{A}^{1.4}+i_c$ (6)

where: i_p—infiltration rate (Lt⁻¹), SA—available storage (L), GI—growth index of crop in percent of maturity, a—an index of surface connected porosity (in.hr.⁻¹ per (in.)^1.4 of storage), i_c—steady state infiltration rate [Lt⁻¹].

SA was calculated using Equations (7) or (8):

$SA=\left({\theta }_s-{\theta }_i\right)d$ (7)

$SA=M_i\times d$ (8)

where: θ_s—saturated soil water content (L³L⁻³), θ_i—volumetric soil water

Table 1. Estimated infiltration model parameters.

Content (L³L⁻³) and d—depth of the sampled soil from the surface layer (L).

2.7. Philip Equation

Cumulative infiltration, for Philip’s model is described using Equation (9):

$I=S{t}^{0.5}+At$ (9)

where; I—cumulative infiltration (L), S—sorptivity (LT^1/2), t—time of infiltration (T) and

A—parameter is known as transmissivity (LT⁻¹).

The infiltration rate was established by differentiating Equation (10).

$\text{d}I/\text{d}t=i\left(t\right)=0.5S{t}^{-0.5}+A$ (10)

Parameters A and S were estimated by plotting the graph of i against t^−0.5 and compared with straight line equation gives A as intercept and S as slope or gradient.

2.8. Green-Ampt Equation

Green and Ampt recommended an equation for infiltration rate developed from Darcy’s law as in Equation (11) [22].

$i=K\left(1+\frac{\mu S_c}{I}\right)$ (11)

where: μ—soil porosity, S_c—capillary suction at the wetting front, and K—hydraulic conductivity.

Equation (11) can be re-written as Equations (12) or (13):

$i=X\left(1+\frac{y}{xI}\right)$ (12)

$i=x+\frac{y}{I}$ (13)

where: x and y are parameter values for Green and Ampt model obtained by plotting a graph of i against 1/I and drawing the best fit lines, where y equals to gradient and x equals to intercept of i-axis. Cumulative (I) values were used to predict i values using Equation (13).

2.9. USDA Soil Conservation Service (SCS) Model

The USDA SCS model has been described by [23] and [24] as one of the empirically developed approaches to water infiltration and is written as Equation (14).

$I=c{t}^d+0.6985$ (14)

where: I—cumulative infiltration rate, c and d—constants.

The constants c andd were determined from observed infiltration data by taking logarithm to both sides of Equation (14) to give Equation (15), and log I was plotted against log t and the equation was compared to linear equation to determine the parameter c and d.

$\log I=\log C+\log 0.6985+d\log t$ (15)

2.10. Data Analysis

Linear regression analysis was undertaken at 95% confidence level whilst Graph Pad Prism 8 was used for the data analysis and Microsoft Excel 10 was used for plotting of graphs.

3. Results and Discussion

3.1. Estimated Infiltration Model Parameters

The model parameters as base values of the various infiltration models are presented in Table 1. Model parameters are very important in predicting the infiltration characteristics of a typical soil condition.

These base values were obtained using fittings of the models constructed from the experimental data with the exception of Holtan model parameters, whose values were obtained from soil hydrologic group. The parameters of the models determined were used to calculate the condition numbers and subsequently test the degree of sensitivity of each model to change in model parameters. The values of the parameters were also used for prediction of models’ infiltration rates.

3.2. Kostiakov Infiltration Model Parameters

The model parameters of the Kostiakov model at the Bontanga irrigation scheme were k = 28.0027 and n = 0.5902, k = 17.0294 and n = 0.4504 for Golinga and k = 23.0356 and n = 0.6339 for Libga irrigation scheme using fittings of Kostiakov Equation (15).

[25] at Hisar and Kurukshetra districts of Haryana, India reported Kostiakov parameters, k = 23.56 and 18.32; and n = 0.69 and 0.36 and noted to be lower than the k and n values obtained at Bontanga, but similar to values reported at Golinga and Libga study sites. Literature reveals that k is consistently ≥ 1 and n is consistently ≤ 1.

The fitting of the model to the infiltration data yielded coefficient of determination of 0.9702 at Bontanga, 0.9565 at Golinga and 0.9876 at Libga study sites as in Figure 1. These indicate good correlation between the observed infiltration and the modeled infiltration, since high values of the coefficient of determination > 95% indicates a good relationship.

[9] proved that the value of n parameter in the Kostiakov equation can be more than unity as oppose to the general opinion that n parameter is found between zero and one. [26], however, found analytically that the value of n was constantly less than one and this was found to be consistent with the findings in this study [12].

3.3. Horton Infiltration Model Parameters

Plots of in (i_o-i_c) against time using Equation (4), aided determination of the parameters of the model using constant infiltration rate of i_c = 2.66 mm/h. At

Figure 1. Kostiakov infiltration model fitting for the study sites. (a) Bontanga; (b) Golinga; (c) Libga.

Bontanga, the i_o parameter was 73.4417 mm/h and the infiltration decay rate (β) = 1.9794 mm/h; whilst at Golinga, the i_o was 45.6589 mm/h and the infiltration decay rate (β) = 1.8578 mm/h and i_c = 0.4 mm/h and with Libga recording i_o of 46.2287 mm/h, β = 1.6332 mm/h and i_c = 5 mm/h.

The model fitting to the infiltration data yielded coefficient of determination of 0.9602 at Bontanga, 0.9404 at Golinga and 0.8781 at Libga as presented in Figure 2. The coefficient of variation of the three (3) sites presents a good relationship between the measured and modeled infiltration.

3.4. Philip Infiltration Model Parameters

Plots of i mm/h against t^−0.5 (h) from which sorptivity (S) of 49.778 and transmissivity (A) of −11.254 mm/h were obtained at the Bontanga study site; S = 41.858 and A = −14.163 mm/h at Golinga study site and S = 33.25 and A = −2.3953 mm/h at Libga sites were obtained as the model parameters Figure 3. Fitting the model to the infiltration data yielded coefficient of determination of 0.9177, 0.9714 and 0.917 for the Bontanga, Golinga and Libga study sites respectively. The coefficient of determination indicates a very good relationship between field and modeled infiltration data of the study sites.

Negative values, for assessed transmissivity, indicate extra infiltration took place and/or exudation of water from the ground was realized as noted by [27] [28]. However, it has been noted to be impossible to interpret the negative value for transmissivity physically. Sites with high negative transmissivity e.g. Bontanga (−11.254 mm/h) and Golinga (−14.123 mm/h) recorded lower infiltration rates of 0.64 mm/h by the models’ prediction as compared to infiltration rates of 9.36 mm/h at Libga with −2.3953 mm/h transmissivity.

3.5. Green and Ampt Infiltration Model Parameters

Green-Ampt model parameters X and Y obtained for the various study sites

Figure 2. Horton infiltration model fitting for the study sites. (a) Bontanga; (b) Golinga; (c) Libga.

Figure 3. Philip infiltration model fitting for the study sites. (a) Bontanga; (b) Golinga; (c) Libga.

were X = 3.325 mm/h and Y = 395.16 mm/h at Bontanga; X = −11.683 mm/h and Y = 341.33 mm/h at Golinga and X = 6.303 mm/h and Y = 210.08 at Libga. The fitting of the model to the infiltration data yielded coefficient of determination of 0.7612 at Bontanga, 0.9007 at Golinga and 0.8068 at Libga study sites as presented in Figure 4. These values indicate a good correlation between the observed infiltration values and the model infiltration values. What these values tell

Figure 4. Green and ampt infiltration model fitting for the study sites. (a) Bontanga; (b) Golinga; (c) Libga.

us is that the model relates better to the infiltration rates data at Golinga than to the infiltration rates at Bontanga and Libga.

These values reflect the prediction ability of the model. Considering Libga and Golinga, the direct association between i, X, Y and I⁻¹ was observed. Comparing Equation (20) to Equation (21) indicates that Y represents physical and hydraulic parameters such as porosity, suction front and hydraulic conductivity. Therefore, increasing these parameters will increase the infiltration rates and vice versa. These parameters were used in Green and Ampt model prediction and determination of the sensitivity of the model to change in its parameter values at the Bontanga study site.

3.6. Soil Conservation Service (SCS) Infiltration Model Parameters

The SCS model parameters c and d obtained through plots of infiltration data were, c = 40.0867 mm/h and d = 0.5902 mm/h for Bontanga; c = 24.3837 mm/h and d = 0.4504 mm/h for Golinga and c = 32.9761 mm/h and d = 0.6339 mm/h for Libga study sites. These parameters were used to predict the models’ infiltration rates and also undertake the model sensitivity analysis. Fitting the model to the infiltration data yielded coefficient of determination of 0.9702 at Bontanga, 0.9565 at Golinga and 0.9876 at Libga as presented in Figure 5. These values indicate an excellent relationship between observed infiltration and modeled infiltration values.

[25] in a similar study obtained c = 48.39 and d = 0.62; c = 34.15 and d = 0.76 at Hisar and Kurukshetra districts of Haryana, India. These values conform with the values obtained in this study on the basis established by [28] that c > 0 and 0 < d < 1.

Figure 5. SCS infiltration model fitting for the study sites. (a) Bontanga; (b) Golinga; (c) Libga.

4. Conclusion

Each model’s parameters for all the six (6) models except Holtan model gave different values at all the three study sites. The parameter GIa and i_c, of Holtan model drawn from hydrologic soil group, were the same at all the study sites because of the similar vegetative cover and surface conditions. The values of the parameters influenced the models’ performance. For all the models, the coefficient of determination indicated good relationships and with the parameter values also influencing models’ performance.

Acknowledgements

The findings of this paper are part of a Ph.D. thesis and the authors are grateful to the staff and students of the Department of Agricultural Mechanisation and Irrigation Technology, Faculty of Agriculture, Food and Consumer Sciences, University for Development Studies, Ghana. Also, field data collection for this work was supported by Staff and Farmers of the Botanga, Goling and Libga Irrigation Projects. This support is therefore very much appreciated by the authors.

Conflicts of Interest

The authors declare that they have no competing interest regarding this publication.

References

[1]	Tuffour, O. (2015) Physically based Modeling of Water Infiltration with Soil Particle Phase. Ph.D. Thesis, Kwame Nkrumah University of Science and Technology, Kumasi.
[2]	Green, W.H. and Ampt, G.A. (1911) Studies in Soil Physics. The Journal of Agricultural Science, 4, 1-24. https://doi.org/10.1017/S0021859600001441
[3]	Kostiakov, A.N. (1932) On the Dynamics of the Confinement of Water Percolation in Soils and on the Necessity of Studying It from a Dynamic Point of View for Purpose of Amelioration. Society of Soil Science, 14, 17-21.
[4]	Horton, R.E. (1941) An Approach toward a Physical Interpretation of Infiltration Capacity. Soil Science Society of America Journal, 5, 399-417. https://doi.org/10.2136/sssaj1941.036159950005000C0075x
[5]	Philip, J.R. (1954) An Infiltration Equation with Physical Significance. Soil Science, 77, 153-158. https://doi.org/10.1097/00010694-195402000-00009
[6]	Holtan, H.N. (1961) A Concept for Infiltration Estimates in Watershed Engineering. Agricultural Research Service, United States Department of Agriculture, USA, 41-51.
[7]	USDA Soil Conservation Service (1972) National Engineering Handbook, Section 4, Hydrology. USDA Soil Conservation Service, Washington DC. USDA-ARS Tech. Bull. No. 1435, Agricultural Research Station, Beltsville, Md.
[8]	Farid, H.U., Mahmood-Khan, Z., Ahmad, I., Shakoor, A., Anjum, M.N., Iqbal, M. M., Mubeen, M. and Asghar, M. (2019) Estimation of Infiltration Models Parameters and their Comparison to Simulate the Onsite Soil Infiltration Characteristics. International Journal of Agricultural and Biological Engineering, 12, 84-91. https://doi.org/10.25165/j.ijabe.20191203.4015
[9]	Hillel, D. (1998) Environmental Soil Physics. Academic Press, San Diego, CA.
[10]	Rawls, W.J., Brakensiek, D.L. and Shirmohammadi, A. (1993) Infiltration and Soil Water Movement. In: Maidment, D.R., Ed., Handbook of Hydrology, McGraw-Hill, New York, USA.
[11]	Skaggs, R.W. and Khaleel, R. (1982) Hydrologic Modeling of Small Watersheds. In: Haan, C.T., Johnson, H.P. and Brakenstek, D.L., Eds., An ASAE Monograph Number 5 in a Series, American Society of Agricultural Engineers, USA.
[12]	Turner, E.R. (2006) Comparison of Infiltration Equations and Their Field Validation with Rainfall Simulation. M.Sc. Thesis, University of Maryland, College Park, Maryland.
[13]	Estefanía, M.H. and André, L.B.C. (2015) Sensitivity Analysis of One-Dimensional Infiltration Models. Electronic Journal of Geotechnical Engineering, 20, 4313-4324.
[14]	Azuka, C.V., Mbagwu, J.S.C. and Oyerinde, G.T. (2013) Infiltration Characteristics and Their Prediction on a Toposequence at Nsukka, South-Eastern Nigeria. International Journal of Science and Advanced Technology, 3, 1-7.
[15]	Adongo, T.A., Kyei-Baffour, N., Abagale F.K. and. Agyare, W.A. (2019) Assessment of Reservoir Sedimentation of Irrigation Dams in Northern Ghana. Lake and Reservoir Management, 36, 87-105. https://doi.org/10.1080/10402381.2019.1659461
[16]	Ghana Irrigation Development Authority (2017) Factsheet of Ghana Irrigation Authority. Ghana Irrigation Development Authority, Ghana.
[17]	Criddle, W.D., Davis, S., Pair, C.H. and Shockley, D.G. (1956) Methods for Evaluating Irrigation Systems. Soil Conservation Service, United States Department of Agriculture, USA, 24 p.
[18]	Naeth, M.A., Chanasyk, D.S. and Bailey, A.W. (1991) Applicability of the Kostiakov Equation to Mixed Prairie and Fescue Grasslands of Alberta. Journal of Range Management, 44, 18-21.
[19]	Frere, M.H., Onstad, C.A. and Holtan, H.N. (1975) ACTMO, An Agricultural Chemical Transport Model. ARS-H-3, Agricultural Research Service, US Department of Agriculture, Washington DC, 54.
[20]	Musgrave, G.W. (1955) How Much of the Rain Enters the Soil? In: The United States Department of Agriculture, Ed., Water: Year book of Agriculture, The United States Department of Agriculture, Washington DC, 151-159.
[21]	Holtan, H.N. and Lopez, N.C. (1971) USDAHL-70 Model of Watershed Hydrology. Technical Bulletins from United States Department of Agriculture, Economic Research Service, Article No. 171863.
[22]	Subramanya, K. (2008) Engineering Hydrology. 3rd Edition, Tata McGraw-Hill Publishing Company Limited, New Delhi.
[23]	Jury, W.A., Gardner, W.R. and Gardner, W.H. (1991) Soil Physics. 5th Edition. John Wiley & Sons, New York.
[24]	Parveen, S., Tiwari, N.K. and Subodh, R. (2017) Estimation and Inter-comparison of Infiltration Models. Water Science, 31, 34-43. https://doi.org/10.1016/j.wsj.2017.03.001
[25]	Parveen, S., Tiwari, N.K. and Subodh, R. (2018) Performance Evaluation of Infiltration Models. Journal of Indian Water Resources Society, 38.
[26]	Mbagwu, J.S.C. (1990) Mulch and Tillage Effects on Water Transmission Characteristics of an Ultisol and Maize Grain Yield in SE. Nigeria. Pédologie, 40, 155-168.
[27]	Machiwal, D., Madan, K.J. and Mal, B.C. (2006) Modelling Infiltration and Quantifying Spatial Soil Variability in a Watershed of Kharagpur, India. Biosystems Engineering, 95, 569-582. https://doi.org/10.1016/j.biosystemseng.2006.08.007
[28]	Haghighi-Fashi, F. Sharifi, F. and Kamali, K. (2014) Modelling Infiltration And Geostatistical Analysis of Spatial Variability of Sorptivity and Transmissivity in a Flood Spreading Area. Spanish Journal of Agricultural Research, 12, 277-288. https://doi.org/10.5424/sjar/2014121-4659