^{1}

^{2}

^{*}

^{2}

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.

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 [

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 [

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 [

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 [

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 [

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 [

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. [

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

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 [

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) [

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.

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 [

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 [

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 [

Horton [

i = i c + ( i o − i c ) e − β 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 ( i − i c ) = ln ( i o − i c ) − β t (4)

Parameters i_{0}, i_{c} and k were obtained by plotting ln ( i − i c ) against t.

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 [_{c}) were obtained from soil hydrologic group built on values founded by [

M i = θ s − θ i (5)

The average moisture deficit (M_{i}) of the control depth, multiplied by the assessed depth gave available storage. Holtan and Lopez [

i p = G I a S A 1.4 + i c (6)

where: i_{p}—infiltration rate (Lt^{−1}), SA—available storage (L), GI—growth index of crop in percent of maturity, a—an index of surface connected porosity (in.hr.^{−1} per (in.)^{1.4} of storage), i_{c}—steady state infiltration rate [Lt^{−1}].

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

S A = ( θ s − θ i ) d (7)

S A = M i × d (8)

where: θ_{s}—saturated soil water content (L^{3}L^{−3}), θ_{i}—volumetric soil water

Model Name | Model Equation | Parameters | Estimated Parameter Value | ||
---|---|---|---|---|---|

Bontanga | Golinga | Libga | |||

Kostiakov | I = K t n ^{ } | K and n | K = 28.0027 n = 0.5902 | K = 17.0294 n = 0.4504 | K = 23.0356 n = 0.6339 |

Philip | I = S t 0.5 + A t | S and A | S = 49.778 A = −11.254 | S = 41.858 A = −14.163 | S = 33.25 A = −2.3953 |

Horton | I = i c t + i o − i c k [ 1 − e − k t ] | i_{o}, i_{c} and k | i_{0} = 73.4417 i_{c} = 2.66 k = 1.9794 | i_{0} = 45.6589 i_{c} = 0.4 k = 1.8578 | i_{0} = 46.2287 i_{c} = 5 k = 1.6332 |

Green and Ampt | I = X + Y / I | X and Y | X = 3.325 Y = 395.16 | X = −11.683 Y = 341.33 | X = −6.3031 Y = 210.08 |

Holtan | I = G I a S A 1.4 + i c | GI, a and i_{c}_{ } | GI = 0.3 a = 1.0 i_{c} = 7.6 | GI = 0.3 A = 1.0 i_{c} = 7.6 | GI = 0.3 A = 1.0 i_{c} = 7.6 |

SCS | I = c t d + 0.6985 | c and d | c = 40.0867 d = 0.5902 | c = 24.3837 d = 0.4504 | c = 32.9761 d = 0.6339 |

Content (L^{3}L^{−3}) and d—depth of the sampled soil from the surface layer (L).

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

I = S t 0.5 + A t (9)

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

A—parameter is known as transmissivity (LT^{−}^{1}).

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

d I / d t = i ( t ) = 0.5 S 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.

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

i = K ( 1 + μ S c I ) (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 ( 1 + y x I ) (12)

i = x + 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).

The USDA SCS model has been described by [

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)

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.

The model parameters as base values of the various infiltration models are presented in

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.

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).

[

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

[

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

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

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

Negative values, for assessed transmissivity, indicate extra infiltration took place and/or exudation of water from the ground was realized as noted by [

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

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

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^{−1} 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.

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

[

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.

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.

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

Salifu, A., Abagale, F.K. and Kranjac-Berisavljevic, G. (2021) Estimation of Infiltration Models’ Parameters Using Regression Analysis in Irrigation Fields of Northern Ghana. Open Journal of Soil Science, 11, 164-176. https://doi.org/10.4236/ojss.2021.113009