Achieving Real-World Saturated Hydraulic Conductivity: Practical and Theoretical Findings from Using an Exponential One-Phase Decay Model

: Obtaining accurate values of saturated hydraulic conductivity ( K sat ) is very important for managing all natural or artiﬁcial processes involving water ﬂow into soils. Double-ring inﬁl-tration (DRI) is one of the easiest-to-work-with techniques commonly used for K sat determination. Unfortunately, when improperly used, it leads to important variations and inaccurate results. This study was designed to investigate the necessary conditions to reach the true-value or real-world saturated hydraulic conductivity ( K sat-real-world ) in the ﬁeld. For this purpose, the effects of two factors—namely, the measured inﬁltration data type (cumulative, instant rate, and average rate) and the related non-linear regression equation type—were analyzed. Measurements with DRI were performed with samples from 106 locations in three West African countries, namely, Burkina Faso, Mali, and Cote d’Ivoire. The soils were composed of loam, sandy loam, and sandy clay loam. The results show that when inﬁltration rates are used rather than cumulative inﬁltration non-linear regression curves, the variability between the measured K sat and the real-world saturated hydraulic conductivity ( K sat-real-world ) could reach from 2.2% to 58.8%. This variability was caused by the approximate ampliﬁcation—according to the procedure used—of time-increment measurement errors. Extending the test duration to more than 4 h, especially when clay soils were involved, and using the exponential one-phase decay non-linear regression of the cumulative inﬁltration data based on a clear measurement protocol provided the K sat values that were closest to K sat-real-world .


Introduction
Saturated hydraulic conductivity (K sat ) stands out as one of the most important characteristics of soil, and its accurate determination is essential to the success of in-soil cropgrowing projects.The saturated hydraulic conductivity has an impact on both the ante-and post-development stages of soil [1].For example, the choice of irrigation system type-e.g., surface irrigation, sprinkler irrigation, or localized irrigation-and the associated cropping pattern (rice, maize, vegetables) suitable for certain soils have to be preceded by measurement campaigns to map K sat values, which will help with decision making.However, the pre-existence of certain types of vegetation or crop also induces changes and impacts on the K sat by changing the soil structure, as observed when a forest is replaced with an agricultural area [2].
The top soil pounding infiltration technique has been extensively used for a long time to determine agricultural and even urban soil properties.In history, the principle of this technique has evolved into two types of equipment.The first was proposed to be a "single ring" (SRI) and is considered by some studies to be as accurate as the second and more common configuration, named "double ring" (DRI), as far as the determination of saturated hydraulic conductivity (K sat ) is concerned [3,4].The usefulness of the technique is now admitted.DRI has been used, for example, to evaluate the difference in infiltration that would result from the decision to apply organic matter in soils [5].In the Sudan savanna, Ikazaki et al. [6] found, by using DRI, that minimum tillage (MT, a minimum soil disturbance and soil cover technique) and sorghum residue mulching (a minimum soil disturbance technique) without intercropping effectively reduced the annual soil loss by 54%.The study was designed to assess the differential importance of each of the three components of conservative agriculture, defined by the Food and Agriculture Organization of the United Nations (FAO) as being (i) minimum soil disturbance, (ii) soil cover, and (iii) crop rotation/association.This result was explained by an improvement in the soil hydraulic conductivity through the boring of termites and wolf spiders found under the sorghum stover mulch.Aiming to assess the effects of agricultural wastes on the hydraulic properties of loamy sand cropland in Turkey, Gülser and Candemir [7] used DRI and found that a high carbon/nitrogen ratio in agricultural waste leads to a significant increase in the initial infiltration rate of loamy sand soil.Exploring the impact of irrigation water's chemical composition on saline-sodic and normal soils, a study implemented both in the laboratory and on fields in India using DRI concluded that saturated hydraulic conductivity could drastically change in soils with a light texture as they were irrigated by chemical water [8].None of these studies targeted the accuracy of DRI.
Other investigations examined the factors that may have an impact on the measurements carried out via DRI.Such factors included the diameter of the inner ring; the depth at which both rings are driven into the soil; the single-ring and double-ring comparative effect; the temperature; the water viscosity; and the time of measurement.In China, in an in situ study using different double-ring infiltrator diameters-20 cm, 40 cm, 80 cm, and 120 cm-and in three types of texture (loamy sand, silt, and silt loam) covering seven measurement sites, through performing numerical simulations using the software HYDRUS-2D model (Version 2.0), Lai and Ren [9] found that the size of the inner ring during DRI significantly alters the measured hydraulic conductivity of heterogeneous soil.The field test used a constant-water-level double-ring infiltrator.The results showed that the variability in the measured hydraulic conductivity was greater for smaller inner rings.The field tests reported that a steady flow state was reached within 90 min on average and that the criterion to terminate the test was that the volume infiltrated in 5 min "remains constant for a 30 min period".In another study performed through simulation using the HYDRUS-2D model, Lai and Ren [10] investigated how the depth at which the double-ring infiltrator is driven into the soil could impact the hydraulic conductivity measurements.The results showed that the inner ring depth had a greater impact on K sat than the outer ring and that a depth of 5-15 cm, as is traditionally used in infiltration measurements, was acceptable, as further pushing down would result in soil and hydraulic conductivity disturbances.The duration of infiltration was 120 min in all the simulations.
Several investigations were performed to compare the accuracy of the instant infiltration rate when using a single-ring infiltrometer (SRI) or a double-ring infiltrator.Wu et al. [11] used simulations based on numerically solving the axisymmetric form of the Richard equation.They considered one-dimensional infiltration as the referential value and compared this with results obtained with an SRI (diameter: 20 cm), a double-ring infiltrator (inner ring: 30 cm; outer ring: 30 cm), and other diameters.The purpose of the study was not specifically to obtain the value of saturated hydraulic conductivity (K sat ).The simulations yielded variations in the relative infiltration rate-expressed as the ratio of instant infiltration over K sat (i/K sat )-as a function of time.The values of K sat were read from the literature for the three types of soils used in the simulations.Among other things, the study concluded that the DRI rates were about 80% of the one-dimensional infiltration vertical infiltration rates, while the infiltration rates of the SRI were f times greater than the one-dimensional infiltration rates, with f being a correction factor dependent on soil initial boundary conditions and ring geometry.Lewis et al. [3] conducted a study to compare the results of a single-ring infiltrometer (SRI) and a DRI.The measurements were performed on four sites in situ, with 4 replicates of 13 different experimental conditions or factor levels [12], which yielded 52 pairs of experimental units.The results showed that SRI and DRI lead to similar saturated hydraulic conductivity K sat values for a wide range of soil conditions.The duration of the SRI tests was from 20 min to 90 min while DRI tests lasted from 30 min to 180 min.Interestingly, this study concluded that "the high variability of infiltration measurements by either method should be factored into study design when assessing infiltration rates".Therefore, one should search for the variability causes through the measurement protocol and data-processing method.
Temperature and the daytime-measurement impacts on DRI measurements for saturated hydraulic conductivity were investigated by Clancy and Alba [13].This study, implemented in the United States, rightly noted that the main issue with the in situ measurement of K sat was its prominent variability.The study performed 67 K sat measurements on sand and loamy sand soils across a temperature ranging from 5 • C to 35 • C. The references used for comparison were the values of K sat , expressed as a function of water viscosity.The results showed that K sat for sandy soils was from 2.0 to 2.9 times higher than the values predicted by the K sat (viscosity) formula.The criterion for terminating the test was that the infiltration remained constant after three time increments of 5 min.The duration of the tests was, on average, 75 min.One of the main findings was that the time of day chosen to perform the tests affects the values.Morning measurements were found to be two times higher than the K sat function of the viscosity, while afternoon measurements were nearly three times higher.
All the above investigations made substantial contributions to the use of the DRI technique to determine soil's saturated hydraulic conductivity K sat .Nevertheless, one cannot determine from the results whether the real-world K sat-real-world was reached, for at least three reasons.Firstly, the duration of the tests was relatively short, generally not exceeding 90 min, while in cohesive soils, especially when clay is involved, it may take days and sometimes weeks to reach saturation.Secondly, none of the mentioned studies discussed the infiltration physical process pathway to select a suitable equation that would yield, after a confrontation with field measurement results, the K sat-real-world .Finally, the measurement protocol (in terms of data collection and processing methods), although recognized to impact the reliability of the results, changes from author to author, and thus may hinder efforts to reach K sat-real-world .The current study, implemented in three West African countries and covering 106 measurement sites, was designed to contribute to filling these gaps.In short, the aim of this study was to build a methodology based on theory and practice in order to obtain the true in situ saturated hydraulic conductivity K sat-real-world.

Investigation Sites
Aiming to ensure potential usefulness for subsequent broader regional rural development activity, the sites of this study were selected for an accurate determination of their saturated hydraulic conductivity K sat based on three main criteria.These criteria were: (i) a regional scope, with at least three countries whose economy is tightly linked to agriculture, (ii) the availability of different soils within or between the sites, and (iii) the existence of an irrigation project.
The three countries were Burkina Faso, Mali, and Cote d'Ivoire (Figure 1).Burkina Faso is a landlocked country in West Africa with an area of 274,200 km 2 .Its economy is largely based on agriculture, which employs 80% of the workforce [14].Mali is also a landlocked country in West Africa, with an area of over 1,240,000 km 2 .The country has a traditional economic system in which the majority of the population engages in subsistence agriculture confined to the riverine area irrigated by the Niger River [15].Cote d'Ivoire is a coastal country in West Africa, covering 322,463 km 2 .Cote d'Ivoire is heavily dependent on agriculture and related activities [16].
with occurrences of sandy clay loam and loam [17,18].The three sites in Mali-namely BaguinedaUp, BaguinedaDwn, and Moria-received 61 measurement points.Their soils are mainly composed of sandy clay loam, with occurrences of loam and sandy clay.Finally, Sema was the unique site selected in Cote d'Ivoire, with 22 measurement points, covering soil made of sandy clay loam, with occurrences of loam.
The sites were either in use or destined for agricultural irrigation.The three sites of Burkina Faso were in use for localized irrigation, microsprinkler irrigation (Kamboinsé), and gravity irrigation (Goupana and Rakaye).In Mali, gravity irrigation was also either practiced (BaguinedaUp and BaguinedaDwn) [19] or planned for implementation (Moria).Gravity irrigation was planned for implementation in the Sema site in Cote d'Ivoire.The soil types are related to the sites selected in each country (Table 1).The three sites of Burkina Faso-namely Kamboinsé, Goupana, and Rakaye-had, respectively, 2, 1 and 20 measurement points.These sites have soils that are essentially made of sandy loam, with occurrences of sandy clay loam and loam [17,18].The three sites in Mali-namely BaguinedaUp, BaguinedaDwn, and Moria-received 61 measurement points.Their soils are mainly composed of sandy clay loam, with occurrences of loam and sandy clay.Finally, Sema was the unique site selected in Cote d'Ivoire, with 22 measurement points, covering soil made of sandy clay loam, with occurrences of loam.The sites were either in use or destined for agricultural irrigation.The three sites of Burkina Faso were in use for localized irrigation, microsprinkler irrigation (Kamboinsé), and gravity irrigation (Goupana and Rakaye).In Mali, gravity irrigation was also either practiced (BaguinedaUp and BaguinedaDwn) [19] or planned for implementation (Moria).Gravity irrigation was planned for implementation in the Sema site in Cote d'Ivoire.

Double-Ring Infiltrometer
The double-ring infiltrometer is a well-known device to determine the soil infiltration curve and the saturated hydraulic conductivity.It has been used for various purposes related to saturated hydraulic conductivity, such as the determination of soil matrix infiltration in forest soils [20], simulations aiming to predict soil hydraulic parameters [21], or in comparing various infiltration devices' biases [22].Although sometimes automated and expensive [23], the basic device is essentially made of an inner ring and an outer ring (Figure 2).For the current investigations, the material was made of an inner ring with a diameter of 20 cm on the PVC pipe segment, and an outer ring with a diameter of 40 cm, also in the PVC pipe segment, with both having a height of 30 cm (Figure 2).However, as shown by Boivin [24], the size of the equipment-if not too small-has very little impact on the measurements.Before starting the measurement, the sites must be selected, and the number of sites will depend on the study area's size.At least 9 measurement sites are required to study the soil parameters of 1.0 ha agricultural land.In this investigation, this figure was applied.Some 24 h before the real measurements started on a specific site, the soil was wetted by slowly pouring several bucket contents of water on the soil in order to raise the moisture toward saturation, thus reducing the testing time before saturation was reached.The wetting had another purpose: the inner and outer rings were pushed 5 cm into the soil using a plastic head hammer.This ensured that water would not circulate from the inner ring toward the outer ring, and from the outer ring to the surrounding soil surface; therefore, vertical flow would not occur at the soil surface.Afterward, a 40 cm height ruler was fastened to the wall of the inner ring so that water-level measurements could be carried out at the same location (Figure 2).The operator also used a stopwatch.
After the installation of the device, the measurements were performed according to a precise procedure.The procedure was as follows: (i) both the inner and the outer rings were filled with water to the edge; (ii) the stopwatch was started by the operator for a chosen measurement time increment.Based on pretests in the field, the time increments were generally set at 10 min for the first 3 measurements, 20 min for the following 3 measurements, 30 min for the third 3 measurements, and so forth.However, depending on how fast or slow the infiltration was occurring, the number of measurements at each time increment could be increased or reduced.After the preset time increment was reached, the stopwatch was stopped, and both the time increment ∆t mi and the related infiltrated water layer ∆h mi were measured (read) by the operator.Then, the water level in both rings was refilled to the edges, and the chronometer was initialized and run again.After a preset time increment elapsed, a new measurement of the quantities (∆h mi , ∆t mi ) was performed.In this study, the operator knew that the saturation zone was reached by calculating the infiltration rate ∆h mi /∆t mi of the measured water layer increment ∆h mi and the related measured time increment ∆t mi , and by noticing that this rate no longer changed for three consecutive measurements.The constancy of the rate denoted that the soil moisture had reached the saturation zone from which the saturated hydraulic conductivity K sat could be drawn.

Measured Variables
Among the field measurement data, three variables are particularly important because they are used to make a comparison between the expressions modeling the infiltration process and to determine the saturated infiltration rate Ksat of the soil.These three variables are the measured cumulative infiltration ht,cum-,meas (mm) (see Equation ( 2)), the measured instant infiltration rate iinst-meas (mm h −1 ) (see Equation ( 3)), and the measured average infiltration rate iavrg-meas (mm h −1 ) (see Equation ( 4)).Before these variables were determined, the primary data measured in the field were the infiltrated water layers and their related time values, as expressed in Equation (1): where is the cumulative water layer infiltrated at the previous time t−1.The index "meas" is related to the measurement data.
The measured instant infiltration rate was given by: where is the instantaneous infiltration rate measured at instant t, expressed as the ratio of Δ , t meas h as determined in Equation ( 1), with the time increment Δt being the

Measured Variables
Among the field measurement data, three variables are particularly important because they are used to make a comparison between the expressions modeling the infiltration process and to determine the saturated infiltration rate K sat of the soil.These three variables are the measured cumulative infiltration h t,cum-,meas (mm) (see Equation ( 2)), the measured instant infiltration rate i inst-meas (mm h −1 ) (see Equation ( 3)), and the measured average infiltration rate i avrg-meas (mm h −1 ) (see Equation ( 4)).Before these variables were determined, the primary data measured in the field were the infiltrated water layers and their related time values, as expressed in Equation (1): where ∆h t,meas (mm) is the incremental water layer, computed as the difference between the water level h t,meas (mm) measured at time t and the water level h t−1,meas (mm) measured at time t−1.The index "meas" is related to the measured data.
The measured cumulative infiltration is expressed by: where h t,cum−meas (mm) is the cumulative water layer infiltrated from the beginning of the test to the instant current time t.This time t is actually the sum (so a cumulation) of individual time increments ∆t meas corresponding to ∆h t,meas , the incremental infiltrated water layer determined in Equation ( 1).The quantity h t−1,cum−meas is the cumulative water layer infiltrated at the previous time t − 1.The index "meas" is related to the measurement data.
The measured instant infiltration rate was given by: where i inst−meas is the instantaneous infiltration rate measured at instant t, expressed as the ratio of ∆h t,meas as determined in Equation ( 1), with the time increment ∆t being the duration between the measurements of the water level h t,meas (mm) and the water level h t−1,meas .
The measured average infiltration rate was provided by: where i avrg−meas (mm h −1 ) is the average infiltration rate measured, expressed as the ratio of the cumulative infiltrated water layer h t,cum−meas (mm) determined in Equation ( 2) and the elapsed time t (h) (from the beginning of the test).
The couples of measured data-(h t,cum−meas , t); (i inst−meas , t); (i avrg−meas , t)-would be plotted against their non-linear regression [12,25] counterparts (cf.Section 3.2) to determine which model will yield a K sat result that is as close as possible to the real-world or true value.

Processed Variables
Similar to the measured variables, three processed variables are of special importance because they can be used to determine which among the measurement expressions and mathematical models would lead to the closest real-world saturated infiltration K sat values.Conversely to the measurement, the first variable to be considered here is the instant infiltration rate because it better analytically models the infiltration process.The evolution of infiltrated water layer speed follows a process that is regularly found in nature [26]: a certain quantity varies at a speed that is proportional to the remaining amount of that quantity.In the case of soil infiltration measurement with a double ring, the infiltration speed depends on the standing water layer in the inner ring (Figure 2).Analytically, this relationship is expressed by Equation (5): where i t (t) (mm h −1 ) is the instant infiltration rate; dh (mm) and dt(h) are, respectively, water-level variation and the related duration.
Many processes occurring in nature follow the pattern described in Equation (5).For example, the determination of soil infiltration rate, the decay of radioactive isotope atoms [27,28], and the turnover of blood cells in HiV-1 infection [29].
It can be shown that the exponential one-phase decay function is a solution of Equation (5).Therefore, the instant infiltration equation is rewritten, yielding the first processed variable: where t(h) is the elapsed time, k sat (mm h −1 ) is the saturated infiltration rate toward which the instant infiltration rate will tend as time increases; i 0 (mm h −1 ) is the initial infiltration rate at the beginning of the tests; b is a coefficient of decay that can be determined through non-linear regression applied to the infiltration data collected from a double-ring test.
In other words, the three parameters (k sat , i 0 and b) of Equation ( 6) can be completely determined by running a non-linear regression on the field data (∆h t,meas , ∆t) measured for Equation (3).The resulting equation is Regr i nst-meas , in which numerical values are assigned to k sat , i 0 , and b.
The cumulative form is obtained from Equation ( 6) by the following analytical summation: This computation yields the second processed variable: where the three parameters (k sat , i 0 and b) of Equation ( 6) can be completely determined by running a non-linear regression on the field data (h t,cum−meas , t), measured for Equation (2).Two things are worth noting when considering Equations ( 6) and ( 8).The first is that k sat , drawn from Equations ( 6) and ( 8), is exactly the same variable, since one equation is derived from the other.Hence, the values of k sat drawn from field data by either method (instant infiltration or cumulative infiltration) should be exactly the same, theoretically.Discrepancies observed between both would mean that errors occurred and must be elucidated.A second aspect that is worthy of notice is that both non-linear equations include a linear part in their curve: (i) for Equation (6), this linear portion is a horizontal line whose value is equal to k sat , while (ii) for Equation (8), this linear portion is an oblique line whose slope is equal to the same saturated hydraulic conductivity k sat .Therefore, being able to identify the moment when the linear portion of each of the curves-materialized by a straight line-appears will be crucial in the determination of the appropriate time to stop the measurements because k sat is reached.
The third processed variable is the average infiltration rate, obtained by running a non-linear regression analysis on the average data i avrg-meas determined in Equation ( 4), as follows: Regr It is worth noting that the three parameters (k sat , i 0 and b) obtained from the three regressions Equations ( 6), (8), and ( 9) are not necessarily expected to be the same.This point was analyzed by plotting the following variables alongside the remaining data: (Regr i inst−meas , t), (Regr h cum−meas , t), and (Regr i avrg−meas , t).

How False Saturated Infiltration Values Are Measured with Various Soil Types
The sample presented in Figure 3 provides insights into several important aspects of the infiltration process.In graphs a-d in the three countries, namely Burkina Faso, Mali, and Cote d'Ivoire, the plots of the cumulative infiltration measured in the field (1) and the related fitted non-linear regression curve (2), the rate of infiltration measured in the field (4), and the related fitted non-linear regression curve (5) are presented.Even based only on visual inspection (the fit parameters are reported in Section 3.2), these figures show, in general, a good fit for the cumulative infiltration data (Figure 3a-d(2)), while the fits of the instant rate of infiltration indicate important gaps compared to the field data (5a-d).For all measurements related to the instant infiltration rate (Figure 3a-d(5)) in the three countries, discrepancies persistently occur at the beginning of the field measurements.This observation points to the volatility of the instant infiltration rate quantity ∆h/∆t because of the higher soil water absorption rate at the beginning [30], causing errors in the measurement of ∆h, but with a greater impact from errors related to the denominator ∆t.This hectic start influences the measured values and the regression fits that are produced based on the instant infiltration rate.The K sat values obtained from the two processesthe one from the non-linear regression of the cumulative infiltration and the one from the instant rate of infiltration, respectively-are not equal.The results yield for Kamb-1, 6.2 mm h −1 versus 11.1 mm h −1 ; for Moria-P01, 7.7 mm h −1 versus 4.2 mm h −1 ; for Moria-P07, 8.2 mm h −1 versus 12.4 mm h −1 ; and for Sema-P05-RCI, 1.8 mm h −1 versus 2.0 mm h −1 .Although the K sat value for cumulative infiltration cannot be ascertained to be the best assessment (rightness) (see Section 3.2) of the real-world K sat-real-world [31], it appears more precise (accuracy) because of the systematic better goodness of fit.The comparison between the average infiltration rate and the measured instant infiltration rate also reveals gaps, whose analysis leads to important clarifications.The two processes are illustrated in Figure 4. Three observations can be made.First, the instant infiltration rate data (Figure 4(4a-d)) tend to be systematically smaller than the average infiltration rate data (Figure 4(6a-d)).Therefore, one would expect a lower saturated hydraulic conductivity Ksat from the instant infiltration rate measurements.It could be that the average infiltration equation (Equation ( 4)) denominator, deduced from measurements from the beginning of the test to a certain elapsed time t (the sum of the incremental Δt) is marked by errors that tend to minimize the real duration related to the absorption of the water height hcum-meas.This involuntary shortening of the duration would lead to higher average infiltration rates.In fact, the same type of error occurs in incremental Δt measurements related to the instant infiltration rate in Equation ( 3), but its magnitude is much smaller than the one related to measurements of the average infiltration rate after 3 h or 5 h.The second observation from the processes depicted in Figure 4 is that, from graphical inspection, it can be seen that the goodness of fit of both infiltration processes is rather bad, at least until the midpoint of the test.The midpoints of the presented sample are 25 h for Kamb-1, 9 h for Moria P01, 5.5 h for Moria P07, and 4 h for Sema P05.It is only near the end of the process that the fit is good.This observation supports the assumption that the saturated hydraulic conductivity is better measured and modeled by non-linear regression when the soil moisture tends toward saturation [26,32].The third and final observation is related to the fact that, as one would expect, instant infiltration rate data with shorter intervals of measurement are much more hectic than the average infiltration rate data.Hence, one would be prone to conclude that the latter would yield a more accurate Ksat than the former [33].However, for the reasons provided in the previous observation, that is not the case (see further in Section 3.2).The Ksat values yielded by the average versus instant infiltration rates for the four sample sites are, respectively, 12.2 mm h −1 versus 11.1 The comparison between the average infiltration rate and the measured instant infiltration rate also reveals gaps, whose analysis leads to important clarifications.The two processes are illustrated in Figure 4. Three observations can be made.First, the instant infiltration rate data (Figure 4a-d(4)) tend to be systematically smaller than the average infiltration rate data (Figure 4a-d(6)).Therefore, one would expect a lower saturated hydraulic conductivity K sat from the instant infiltration rate measurements.It could be that the average infiltration equation (Equation ( 4)) denominator, deduced from measurements from the beginning of the test to a certain elapsed time t (the sum of the incremental ∆t) is marked by errors that tend to minimize the real duration related to the absorption of the water height h cum−meas .This involuntary shortening of the duration would lead to higher average infiltration rates.In fact, the same type of error occurs in incremental ∆t measurements related to the instant infiltration rate in Equation ( 3), but its magnitude is much smaller than the one related to measurements of the average infiltration rate after 3 h or 5 h.The second observation from the processes depicted in Figure 4 is that, from graphical inspection, it can be seen that the goodness of fit of both infiltration processes is rather bad, at least until the midpoint of the test.The midpoints of the presented sample are 25 h for Kamb-1, 9 h for Moria P01, 5.5 h for Moria P07, and 4 h for Sema P05.It is only near the end of the process that the fit is good.This observation supports the assumption that the saturated hydraulic conductivity is better measured and modeled by non-linear regression when the soil moisture tends toward saturation [26,32].The third and final observation is related to the fact that, as one would expect, instant infiltration rate data with shorter intervals of measurement are much more hectic than the average infiltration rate data.Hence, one would be prone to conclude that the latter would yield a more accurate K sat than the former [33].However, for the reasons provided in the previous observation, that is not the case (see further in Section 3.2).The K sat values yielded by the average versus instant infiltration rates for the four sample sites are, respectively, 12.2 mm h −1 versus 11.1 mm h −1 for Kamb-1; 24.2 mm h −1 versus 14.2 mm h −1 for Moria-P01; 18.7 mm h −1 versus 12.4 mm h −1 for Moria-P01; and 6.3 mm h −1 versus 2.0 mm h −1 for Sema-P05.Finally, as explained in the comments relating to Figure 3, the instant infiltration rate does not lead to the real-world K sat-real-world value.
Hydrology 2023, 10, x FOR PEER REVIEW 10 of 21 mm h −1 for Kamb-1; 24.2 mm h −1 versus 14.2 mm h −1 for Moria-P01; 18.7 mm h −1 versus 12.4 mm h −1 for Moria-P01; and 6.3 mm h −1 versus 2.0 mm h −1 for Sema-P05.Finally, as explained in the comments relating to Figure 3, the instant infiltration rate does not lead to the real-world Ksat-real-world value.With 2.5% organic matter content, Saxton [34] made the following classification of soils on the texture triangle: loam for Ksat from 8.6 to 46.9 mm/h; clay loam for Ksat from 2.01 to 7.79 mm/h; and sandy clay loam for Ksat between 3.94 and 20.06 mm/h.It can be noted that these three types of soil have overlapping areas in terms of Ksat, and for each type of soil, a whole range of Ksat exists.For example, the 6.2 mm h −1 of the site Kamb-1 is compatible with a loam, a clay loam, or even a sandy clay loam.Only a grain size analysis [35] will allow for a precise classification of the soils, though this does not hinder the importance of obtaining a real-world Ksat that could be used, for example, to design sprinkler irrigation systems [36].

Tracking the Real-World Ksat Value
The graphs displayed in Figure 5 are related to data samples collected and processed from 12 different testing sites out of the 106 addressed by this investigation.As shown in Table 1, many of the 106 investigation sites were located on similar soils.In each of the 12 graphs, the data are displayed, with identifying numbers and graph markers provided in the legend at the bottom of the graph.These numbers are as follows: "1" for the cumulative infiltrated water layer-Equation (2); "4" for the measured instant infiltration rate-Equation (3); and "6" for the measured average infiltration rate-Equation (4).Based on these three series of data, a non-linear regression computation-Equations ( 6), (8), and (9)-was performed, and the results are presented in three other graphs, connected with lines.These graphs are identified by the following numbers: "2" for the non-linear regression related to the measured cumulative infiltration hcum-meas, "5" for the non-linear With 2.5% organic matter content, Saxton [34] made the following classification of soils on the texture triangle: loam for Ksat from 8.6 to 46.9 mm/h; clay loam for Ksat from 2.01 to 7.79 mm/h; and sandy clay loam for Ksat between 3.94 and 20.06 mm/h.It can be noted that these three types of soil have overlapping areas in terms of Ksat, and for each type of soil, a whole range of Ksat exists.For example, the 6.2 mm h −1 of the site Kamb-1 is compatible with a loam, a clay loam, or even a sandy clay loam.Only a grain size analysis [35] will allow for a precise classification of the soils, though this does not hinder the importance of obtaining a real-world Ksat that could be used, for example, to design sprinkler irrigation systems [36].

Tracking the Real-World K sat Value
The graphs displayed in Figure 5 are related to data samples collected and processed from 12 different testing sites out of the 106 addressed by this investigation.As shown in Table 1, many of the 106 investigation sites were located on similar soils.In each of the 12 graphs, the data are displayed, with identifying numbers and graph markers provided in the legend at the bottom of the graph.These numbers are as follows: "1" for the cumulative infiltrated water layer-Equation (2); "4" for the measured instant infiltration rate-Equation (3); and "6" for the measured average infiltration rate-Equation ( 4).Based on these three series of data, a non-linear regression computation-Equations ( 6), (8), and ( 9)was performed, and the results are presented in three other graphs, connected with lines.These graphs are identified by the following numbers: "2" for the non-linear regression related to the measured cumulative infiltration h cum−meas , "5" for the non-linear regression on the instant infiltration rate data, and "7" for the non-linear regression run on the average infiltration rate data.To show the mathematical expressions of the 36 regression curves in Figure 5, 12 exponential one-phase decay regression equations [32] are presented in Table 2.In this table, only one infiltration test site was selected from each country.
Hydrology 2023, 10, x FOR PEER REVIEW 11 of 21 regression performed on the instant infiltration rate data, and "7" for the non-linear regression run on the average infiltration rate data.To show the mathematical expressions of the 36 regression curves in Figure 5, 12 exponential one-phase decay regression equations [32] are presented in Table 2.In this table, only one infiltration test site was selected from each country.As can be clearly seen from the asymptotic lines of the non-linear regression curves shown in Figure 5 and the sample extracts in Table 2, for the same infiltration measurement site, the saturated hydraulic conductivity Ksat is different for the three regression  As can be clearly seen from the asymptotic lines of the non-linear regression curves shown in Figure 5 and the sample extracts in Table 2, for the same infiltration measurement site, the saturated hydraulic conductivity K sat is different for the three regression processes, and sometimes the difference is rather significant.Based on Table 2, for example, K sat is 6.2 mm h −1 at Kamb-1-Burkina Faso if cumulative infiltration measurement data are used, while it would be, respectively, 11.10 mm h −1 and 12.21 mm h −1 if instant and average infiltration measurement data were used.The results for the site of BaguinedaUp-Mali are very close, even though the value obtained from the cumulative infiltration curve is the smallest.The results for this site yield a K sat of 3.1 mm h −1 when cumulative infiltration measurement data are used, and 3.1 mm h −1 and 3.8.mm h −1 , respectively, if instant average infiltration measurement data are used.However, the difference becomes significant again for the site of Sema-P1-Cote d'Ivoire: 2.3 mm h −1 with cumulative infiltration data, and 12.4 mm h −1 and 8.0 mm h −1 , respectively, when instant and average infiltration measurement data are used.

The Best Saturated Hydraulic Conductivity Results
Considering the graphs shown in Figure 5, two important questions need to be addressed to derive useful conclusions about the saturated hydraulic conductivity measurement process: (i) What can be considered real-world K sat-real-world ?(ii) Which of the three procedures-deriving K sat from non-linear regressions of h cum−meas , i inst-meas , and i avrg-meas -is better in terms of reaching K sat-real-world ?In real field experiments, when the soil moisture is brought close to saturation before the double-ring infiltration test starts, the cumulative infiltration curve remains a straight line, as can be seen in Figure 5e.It was noticed that, under these conditions, the values of the three regression lines were generally close, as for BaguinedaUp-Mali.From the graphs in Figure 5, it can also be deduced that two factors are very important: (a) the measurement of time increments ∆t (particularly in the curvature zone), and (b) the number of available data for the asymptote branch of the curve.When ∆t is constant and ∆h is constant, measurements are generally inaccurate and K sat will not be valid.To follow the real-world process of infiltration, the operator must exponentially dilate either ∆t or decrease ∆h to capture the filling of the exponential decrease in infiltration.When the linear segment of the infiltration rate curves is reached, then the ratio ∆h/∆t must remain constant, regardless of the values of ∆h and ∆t.Therefore, when the ratio (∆h/∆t) meas -computed directly from the measured data in the field-remains constant for three consecutive measurements, the soil saturation is reached and this constant value of (∆h/∆t) meas can be considered the real-world saturated hydraulic conductivity K sat-real-world .Reaching this value can be time-consuming in practice, especially when dealing with clayey soils.This is why regression analysis is very relevant: its prediction allows for the determination of a K sat that is as close as possible to K sat-real-world.
In order to answer the question "Which of the three procedures-deriving K sat from non-linear regressions of h cum−meas , i inst-meas , and i avrg-meas -is better in terms of reaching K sat-real-world ?", comparative figures of K sat are shown in Table 3.In the table, three quantities assessing the goodness of fit of non-linear regression are shown in columns 6-8.These are as follows: (i) R 2 , the coefficient of determination (the closer to 1, the better), (ii) SSE, the error sum of squares (the smaller, the better) [37], and (iii) MSE, the mean of squared error (the smaller, the better) [38].In the ninth column, the values of the saturated hydraulic conductivity K sat -derived from their respective non-linear equations-are shown for each of the 12 individual sites and each of the three procedures.In the 10th column that follows K sat , the average of the last three measurements (∆h/∆t) meas is presented.These values, as explained in the previous paragraph, can be considered the real-world or true saturated hydraulic conductivity K sat-real-world .The gap or discrepancy between K sat and K sat-real-world is shown in the last column of the table.For all 12 sites, it appears that the procedure that best estimates the saturated hydraulic conductivity is the one using the cumulative infiltration measurement data (Table 3(1-12a)).The gaps in absolute values vary from 2.2% to 58.8%.The worst estimation procedure is the averaging one, with a maximum gap reached on the site of Sema of 444.4%.The instant infiltration procedure is located between the two previous ones.Therefore, it can be stated that the cumulative infiltration measurement data procedure is the best way to estimate the saturated hydraulic conductivity of soil K sat .Once that is stated, the question is: how is the cumulative infiltration procedure more accurate than the instant infiltration rate procedure?A discussion based on both analytical expressions and graphics will help in answering this question.It is helpful to use the cumulative infiltration curve as the basis to assess which of the two non-linear regressions, Equations ( 6) and ( 8), will yield the most accurate saturated hydraulic conductivity.In Figure 6, it can be seen that errors occur during both the waterheight interval measurements ∆h mi and the time interval measurements ∆t i .Therefore, each individual segment of the measurement curve can be computed using the Pythagorean theorem [39] for rectangular triangles, as shown in Equation (10): where ∆C i is the length of individual segments of the measurement curve; ∆h i is the infiltrated water increment; and ∆t i is the time increment.
height interval Δhmi and the time interval measurements Δti.Therefore, each individual segment of the measurement curve can be computed using the Pythagorean theorem [39] for rectangular triangles, as shown in Equation ( 10): where ΔCi is the length of individual segments of the measurement curve; Δhi is the infiltrated water increment; and Δti is the time increment.Time and water head are independent variables at the macroscopic scale, and, as such, can be measured simultaneously with an arbitrary accuracy [40].However, in practice, water head measurement errors are much more easily minimized by a good operator and an accurate limnimetric device (Figure 2).Therefore, it is assumed that most errors only occur during the cumulative time interval measurements i m t Δ .With these hypothe- ses, the positive error i c ε between each individual segment k c Δ of the measurement curve when plotted against the corresponding regression curve segment (Figure 6) can be calculated using Equation (11):  Time and water head are independent variables at the macroscopic scale, and, as such, can be measured simultaneously with an arbitrary accuracy [40].However, in practice, water head measurement errors are much more easily minimized by a good operator and an accurate limnimetric device (Figure 2).Therefore, it is assumed that most errors only occur during the cumulative time interval measurements ∆t m i .With these hypotheses, the positive error ε c i between each individual segment ∆c k of the measurement curve when plotted against the corresponding regression curve segment (Figure 6) can be calculated using Equation (11): where ε c i is the error in absolute value (hence, positive) or the gap between each individual segment of the measurement curve and the (theoretically adjusted) non-linear regression curve, while ε t 1 , ε t 2 , ε t 3 , are the curve segment errors induced by time errors during the measurements of time increments ∆t m1 , ∆t m2 , ∆t m3 , such as where t mi is the cumulative measured time ("m" stands for measured); t i is the theoretical (true) cumulative time; ε t i is the time measurement error.
From the construction graph in Figure 6, the following can be derived: By inserting Equation ( 13) into Equation ( 11): The structure of Equation ( 14) two important aspects that can help to ensure more accurate infiltration measurements.First, as in any time series [41][42][43], the error ε C i at time step i depends on the errors in current measurements t mi , but also on measurements at previous times steps t mi−1 and t mi−2 .As written by Montgomery [41], in time series, a future observation depends on the current and the past observations, i.e., adjacent observations are dependent.There is an autocorrelation in the data.Secondly, the second and the third terms on the right side of Equation ( 14) constitute the random part of the error [41], with the deterministic part being t i + t i−1 , and show that the greater the interval between two measurements, the larger the gap between the real-world measurement and the regression curves.This result was expected since the straight line forming the curve segment ∆C i will be longer and difficult to fit into a curved segment of the regression line.Therefore, ideally, to minimize the gaps between the two curves by altering ∆t i , one would (i) proceed with independent measurements of time intervals by bringing the chronometer (and refilling the infiltrometer to reset the water head) to zero at each new measurement and (ii) choosing smaller time intervals.This last point, i.e., shortening the time interval, leads to two important proposals.The first is related to when it is better to keep short time intervals during the measurements.Field data (Figure 5) indicate that the change in the direction of the cumulative curve is faster at the beginning of the measurements and during the 4-5 following hours (depending on the soil type) than later on.This means that to keep the measurement curve as close as possible to the non-linear regression curve, one needs to use short time intervals ∆t mi at the beginning.However, while the soil becomes more and more saturated, the change in the direction of the cumulative regression curve slows down and the curve moves toward a straight line, with a constant slope, providing the expression of the saturated hydraulic conductivity: where K sat is the saturated hydraulic conductivity; ∆h msat is the infiltrated water layer measured while the straight-line zone is reached; and ∆t msat is the time interval related to ∆h msat .Hence, when the soil gets closer to saturation [44], longer time intervals are suitable to ensure an accurate fitting between the measurement curve and the non-linear regression curves (Figure 6).
From the elements of analysis presented in Figure 6 and Equation ( 6), one can also see how drawing the saturated hydraulic conductivity from the instantaneous infiltration rate measurements would lead to inaccurate results.In effect, each point of the infiltration rate is obtained by the slope: in which i mi is the "measured" (in reality, "calculated") instant infiltration rate; ∆h mii is, the measured infiltrated water layer; and ∆t mi is the time interval related to ∆h mi .The related error ε i mi is calculated (rather than measured) by applying the differential computation to Equation ( 16), which is a quotient, yielding: As shown in the expression of Equation (17), the smaller the time interval of measurements ∆t mi , which is suitable at the beginning of the infiltration measurements, as previously stated, the greater the error ε i mi in the soil infiltration rate data, resulting in important changes in the direction of the measurement curve.Due to these changes, the fitting between the infiltration rate measurement data and the related non-linear regression curve is generally worse than the discrepancies observed with the cumulative infiltration curves.As a result, the saturated infiltration value K sat obtained from the measured data and the regression equation present important, irreducible gaps.Therefore, even though standardization institutions such as ASTM [45] suggest deriving the permeability (saturated hydraulic conductivity) from infiltration rate data, this method usually leads to erroneous values.Hence, it is not surprising that, in its standard sheet D 5093-02 related to double-ring infiltration measurements, the values of K sat for i are usually so scattered that ASTM did not provide any estimate of the errors.

A Protocol for Accurate In Situ Measurement of Soil Saturated Hydraulic Conductivity K sat
When making infiltration measurements, it is always critical to determine when the investigation should be terminated [9,13].As shown in Equation (15) and Figure 5, the saturated hydraulic conductivity K sat is reached when the ratio of the hydraulic head variation and the related measured time interval becomes constant, or when the cumulative curve starts to change into a straight line.However, based on the experiments implemented in the three countries in this study, the figures indicate that when three consecutive computations of this ratio provide a constant, same value, the measurements could be stopped and K sat can be derived from Equation (8).Therefore, a full procedure-considering the necessity of (i) having shorter time intervals at the beginning of the measurements and longer time intervals near the saturated zone, (ii) ending the trials at the right time when the infiltration process reaches saturation-is proposed in Tables 4 and 5.

Steps Operation
Step 1 Wet the floor for at least 2-3 h and 24 h in advance Step 2 Set up two concentric rings by driving them 5-10 cm on wet ground and keeping them horizontal at a mason's level.Water must not be able to circulate from one ring to another or leak outside the guard ring.
Step 3 Tracing with an indelible marker, mark the place where all the water-level descent measurements ∆h will be made with a small line.
Step 4 Fill the two rings with water to the brim Step 5 Start the stopwatch and make the first measurement of the series I by stopping the stopwatch after ∆t = 10 min and reporting the corresponding ∆h (mm) Step 6 Immediately refill the two rings to the brim (to reduce the dead time error between two measurements to almost zero) Step 7 Repeat sequences 5-6 for the second then the third 10 min of series I Step 8 Apply sequences 5-7 for series II (with ∆t = 20 min), then for series III (with ∆t = 30 min), and so on until series VI, VII, and VIII.One can increase the number of measurements in a series or skip a series depending on the rate at which the water infiltrates.For example, one can repeat 20 min four times in series II, or go from series I to series III Step 9 The tests are stopped when a quick calculation shows that three successive values of the ratio ∆h (mm)/∆t (h) are equal or constant.This means that the permeability K sat-real-world has been reached, which must be confirmed by the non-linear regression curve.
The way to proceed with the field operation is described in Table 4, while the data should be written as shown in Table 5.The time measurements are clustered into series from I to V (Table 5) allowing for a well-organized and accurate experiment.In each series, the time interval is constant.For example, in series I, this is 10 min, in series II, this is 20 min, in series III, this is 30 min, etc.In a real field experiment, only time intervals (3rd column) are measured and the values of the cumulative time T(h) in hours (4th column) are deduced from the intervals.In the 5th column, the infiltrated water head increments are measured and correspond to the time intervals of the 3rd column.The cumulative infiltrated water head values in the 6th column are computed from the values in the 5th column.Finally, only data from column 4 and column 6 are used to plot the measurement curve and compute, and, against those data, a non-linear regression curve is plotted using, for example, Minitab [46] or XLSTAT 2014.5.03 software [47].

Conclusions
Targeting accuracy in the measurement of saturated hydraulic conductivity K sat , this study used data collected from 106 measurement locations in three West African countries.This how inaccurately K sat is measured in many studies when (i) the test duration is too short, (ii) the infiltration rate and average infiltration rate data, rather than cumulative infiltration data, are used, (iii) non-linear regression equation of infiltration rates instead of cumulative infiltration data are used, and (iv) the criterion for the test termination is not clearly specified.The measurement duration varied in this study from 4 to 50 h, while most of the tests found in the literature would not last more than 3 h, not allowing for the infiltration rate to reach K sat , especially when clay is involved.When infiltration rates are used instead of cumulative infiltration non-linear regression curves, the gap between the measured K sat and the real-world saturated hydraulic conductivity K sat-real-world was found to vary from 2.2% to 58.8%, even reaching 444.4% when average infiltration rate data were used.Using differential mathematical analysis, the study explained how operator errors in time measurement affect the infiltration rate more than cumulative infiltration curves and led to the conclusion that the cumulative infiltration measurement data procedure is the best way to estimate the saturated hydraulic conductivity of a soil K sat .The study proposed a clear investigation protocol that would help to achieve real-world saturated hydraulic conductivity K sat-real-world during investigations.

Figure 1 .
Figure 1.The seven investigation locations within three countries (Burkina Faso, Mali, and Cote d'Ivoire).

Figure 1 .
Figure 1.The seven investigation locations within three countries (Burkina Faso, Mali, and Cote d'Ivoire).

Figure 2 .
Figure 2. Double ring for soil infiltration measurements.Note: both rings are pushed 5 cm into the soil.

Figure 2 .
Figure 2. Double ring for soil infiltration measurements.Note: both rings are pushed 5 cm into the soil.

Hydrology 2023 ,
10, x FOR PEER REVIEW 9 of 21(rightness) (see Section 3.2) of the real-world Ksat-real-world[31], it appears more precise (accuracy) because of the systematic better goodness of fit.

Figure 3 .
Figure 3. Cumulative infiltration (1 and 2) versus instant infiltration rates (4 and 5).Measured data have no connecting lines (1 and 4), while non-linear regression curves (2 and 5) have connected markers.Numbers (1)-(5) refer to the individual graphs in the figure.The letters a-d refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.− = d(Regr ) / cum m inst eas i h d t

Figure 3 .
Figure 3. Cumulative infiltration (1 and 2) versus instant infiltration rates (4 and 5).Measured data have no connecting lines (1 and 4), while non-linear regression curves (2 and 5) have connected markers.Numbers (1)-(5) refer to the individual graphs in the figure.The letters (a-d) refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.i inst = d(Regr h cum−meas )/dt.

Figure 4 .
Figure 4. Average infiltration (6 and 7) versus instant infiltration rates (4 and 5).Measured data have no connecting lines (4 and 6), while non-linear regression curves (5 and 7) have connected markers.Numbers (1) to (7) refer to the individual graphs in the figure.The letters a-d refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.− = d(Regr ) / cum m inst eas i h d t.

Figure 4 .
Figure 4. Average infiltration (6 and 7) versus instant infiltration rates (4 and 5).Measured data have no connecting lines (4 and 6), while non-linear regression curves (5 and 7) have connected markers.Numbers (1) to (7) refer to the individual graphs in the figure.The letters (a-d) refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.i inst = d(Regr h cum−meas )/dt.

Figure 5 .
Figure 5. Infiltration measurement results for 12 locations in Burkina Faso, Mali, and Cote d'Ivoire.Comparing (i) measured data (1, 4, 6) vs. their corresponding non-linear regression fits (2, 5, 7), (ii) the non-linear regression fits, and (iii) the asymptotic saturated hydraulic conductivity Ksat values.Measured data are not connected with drawing lines.Numbers (1)-(7) refer to the individual graphs in the figure.The letters a-l refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.

Figure 5 .
Figure 5. Infiltration measurement results for 12 locations in Burkina Faso, Mali, and Cote d'Ivoire.Comparing (i) measured data (1, 4, 6) vs. their corresponding non-linear regression fits (2, 5, 7), (ii) the non-linear regression fits, and (iii) the asymptotic saturated hydraulic conductivity K sat values.Measured data are not connected with drawing lines.Numbers (1)-(7) refer to the individual graphs in the figure.The letters (a-l) refer to a specific site.The equation and relative curve (3) are not used because it does not bring new insight.It expressed instant infiltration as the derivative of the cumulative infiltration curve.i inst = d(Regr h cum−meas )/dt.

Figure 6 .
Figure 6.Time errors' impact on the fittings of the measured and regressed cumulative infiltration curves.

ε 1 t ε , 2 t ε , 3 tε
is the error in absolute value (hence, positive) or the gap between each indi- vidual segment of the measurement curve and the (theoretically adjusted) non-linear regression curve, while , are the curve segment errors induced by time errors during the measurements of time increments

Figure 6 .
Figure 6.Time errors' impact on the fittings of the measured and regressed cumulative infiltration curves.

Table 1 .
Site locations and main soils.

Lat (DD) Project Area A (ha) Nbr of Meas. Locations HWSD/USDA-Soil Texture (30 cm Top Soil) Dominant Soil Drainage Class (Slope 0.0-0.5%) (Total: 106) Dominant Occurrence
(mm) is the cumulative water layer infiltrated from the beginning of the test to the instant current time t.This time t is actually the sum (so a cumulation) of indi- Δ , t meas h (mm) is the incremental water layer, computed as the difference between the water level , t meas h (mm) measured at time t and the water level −1, t meas h (mm) measured at time t−1.The index "meas" is related to the measured data.The measured cumulative infiltration is expressed by: vidual time increments Δ meas t corresponding to Δ , t meas h , the incremental infiltrated water layer determined in Equation (1).The quantity − − 1, t cum meas h

Table 2 .
Nine non-linear regression infiltration equations from three sites in Burkina Faso, Mali, and Cote d'Ivoire.

Table 3 .
Comparing three procedures of regression for real-world K sat from three countries.

Table 4 .
Double-ring infiltration protocol of accurate K sat measurement.

Table 5 .
Protocol of time for infiltration measurements with examples of data.