A Compact Weighing Lysimeter to Estimate the Water Inﬁltration Rate in Agricultural Soils

: Inﬁltration estimation is made by tests such as concentric cylinders, which are prone to errors, such as the lateral movement under the ring. Several possibilities have been developed over the last decades to compensate these errors, which are based on physical, electronic, and mathematical principles. In this research, two approaches are proposed to measure the water inﬁltration rate in a silty loam soil by means of the mass values of a lysimeter weighing under rainfall conditions and different moisture contents. Based on the fact that with the lysimeter it is possible to determine acting soil ﬂows very precisely, then with the help of mass conservation and assuming a downward vertical movement, 12 rain events were analyzed. In addition, it was possible to monitor the behavior of soil moisture and to establish the content at ﬁeld capacity from the values of the weighing lysimeter, from which both approach are based. The inﬁltration rate of these events showed a variable rate at the beginning of the rainfall until reaching a maximum, to descend to a stable or basic rate. This basic inﬁltration rate was 1.49 ± 0.36 mm/h, and this is because soils with ﬁne textures have reported low inﬁltration capacity. Four empirical or semi-empirical models of inﬁltration were calibrated with the values obtained with our approaches, showing a better ﬁt with the Horton’s model.


Introduction
To understand the behavior of the hydrological cycle in the soil remains a challenge for science, specifically the movement of water through the soil and its capacity to retain it. To estimate the water exchange between a well-defined portion of soil and other physical systems to which water is transferred through evaporation, transpiration, percolation or drainage, several methods as gravimetric, tensiometry, humidity sensing, and lysimetry have been proposed in the specialized literature [1][2][3].
Infiltration is the hydrological process that describes the entry of water into a soil, and the amount of water that enters the soil in a given time represents the infiltration rate [4,5]. It is variable in space and time, due to hydrodynamic and physical soil properties, flow conditions, and cultivation practices, among others [6][7][8][9]. As water enters the soil profile through the forces of adhesion, the interaction of water and soil particles, and gravity, it fills the pores of the soil, defining a wet front as it passes through the soil profile. If the soil has a good structure and large pores, its infiltration rate will be high (as in the case of coarse-textured soils) and it will reach a deeper wet front compared to a fine-textured soil whose pores are smaller. It is assumed that if the soil is at its lower limit of initial The soil texture was classified as silt loam, based on the soil percentages obtained with Bouyoucos hydrometer method and the USDA texture triangle figure [40]. Percentages of soil granulometry are 27.1% sand, 51.5% silt, and 21.4% clay; an organic matter of 2.77%, and a bulk density of 1.38 g/cm 3 . The hydrodynamic characteristics of the field capacity and permanent wilting point were obtained by the membrane and pressure plate apparatus, which resulted in 0.35 cm 3 /cm 3 and 0.16 cm 3 /cm 3 , respectively. The initial soil moisture content on February 8 were the experiment started was inferred knowing the bulk density and total volume with 0.23 cm 3 /cm 3 .

Materials
The compact weighing lysimeter, LISITN-AC-E1 model (teleNatura, Alicante, Spain), consists of two vessels. The upper vessel called soil vessel has dimensions of 0.56 m × 0.96 m and 0.35 m depth in the center, which stores a volume of the reconstituted soil from the plot trying to keep the same conditions from where it is extracted, with a resolution of 20 g or 0.033 mm ( Figure 2). The soil is isolated of its environment in such a manner that the lateral and capillary rise underground flows are zero [26], so that the water balance terms can be accurately obtained.  The soil texture was classified as silt loam, based on the soil percentages obtained with Bouyoucos hydrometer method and the USDA texture triangle figure [40]. Percentages of soil granulometry are 27.1% sand, 51.5% silt, and 21.4% clay; an organic matter of 2.77%, and a bulk density of 1.38 g/cm 3 . The hydrodynamic characteristics of the field capacity and permanent wilting point were obtained by the membrane and pressure plate apparatus, which resulted in 0.35 cm 3 /cm 3 and 0.16 cm 3 /cm 3 , respectively. The initial soil moisture content on February 8 were the experiment started was inferred knowing the bulk density and total volume with 0.23 cm 3 /cm 3 .

Materials
The compact weighing lysimeter, L ISI TN-AC-E 1 model (teleNatura, Alicante, Spain), consists of two vessels. The upper vessel called soil vessel has dimensions of 0.56 m × 0.96 m and 0.35 m depth in the center, which stores a volume of the reconstituted soil from the plot trying to keep the same conditions from where it is extracted, with a resolution of 20 g or 0.033 mm ( Figure 2). The soil is isolated of its environment in such a manner that the lateral and capillary rise underground flows are zero [26], so that the water balance terms can be accurately obtained.
The second vessel is located below the previous one. It is called drainage vessel and serves to store the excess of water drained by the soil profile with a resolution of 1 g (or 0.002 mm). The soil vessel has a hole in the lower central part that allows the water flows to an electrovalve, then goes to a funnel and finally reaches the drainage vessel, where it will be stored and weighed (Figure 2c). The drainage vessel was configured to store only 3000 g of water, so when this increase is reported, the electrovalve located under it is activated and starts emptying it, while this occurs the electrovalve of the soil vessel is closed to avoid losses of water quantification. Once the emptying process has been completed, the electrovalves return to their initial state, opening the soil vessel one and closing the drainage vessel one, allowing the water flow from the soil vessel to the drainage vessel.
The mass variations of both vessels are detected by load cells, four for the soil vessel and one for the drainage vessel. These load cells send an electrical signal to the datalogger (CR300 model, Campbell Scientific, Logan, UT, USA), according to the deformation suffered by the object that they are supporting. This electrical signal is converted to a mass value by the same datalogger. According to the Spanish norm UNE 500520 [41], the mass data of the two vessels was sampled every second and the average value of 60 s (one minute) was registered as an instantaneous value. This guarantees the maximum quality and representativeness of the data obtained. For further details of the compact weighing lysimeter see the work of Nicolás-Cuevas et al. [42]. The compact weighing lysimeter, LISITN-AC-E1 model (teleNatura, Alicante, Spain), consists of two vessels. The upper vessel called soil vessel has dimensions of 0.56 m × 0.96 m and 0.35 m depth in the center, which stores a volume of the reconstituted soil from the plot trying to keep the same conditions from where it is extracted, with a resolution of 20 g or 0.033 mm ( Figure 2). The soil is isolated of its environment in such a manner that the lateral and capillary rise underground flows are zero [26], so that the water balance terms can be accurately obtained.

Flow Calculation
The lysimeter contains an isolated soil profile, so that is possible to determine the water inputs and outputs in the soil based on the mass data records, through the mass balance equation, shown in Equation (1): where ± dS dt is the change of soil water storage (g), dP dt is the amount of water by precipitation (g), dIR dt is the amount of water by irrigation (g), dD dt is the amount of water drainage beyond the root zone (g), dET dt is the amount of water by the crop evapotranspiration (g), and dSr dt is the surface runoff (g) along interval of time [6].
In this work, only the rainfall events were analyzed to know the rate of entry of water into the soil. The rainfall has a longer duration of application and therefore gives a better accuracy of the results. So, variables as irrigation and evapotranspiration were not considered, except in determining θ FC , because they cannot occur when there is rain [32,43], so we reformulate Equation (1) as Equation (2): Surface runoff is neglected because it was not observed here. The only known variables are the accumulated masses of the soil vessel (S) and the drainage vessel (D), so the flow is determined by the increase in the mass (m) of the vessel over time [32], as the Equation (3): where dD dt is the flow of drainage vessel over time, dS dt is the flow of the soil vessel over time, t i is the time in the instant i and t i+1 is the time in the instant i + 1.
Rainfall is calculated based on literature reports [27,32,43,44], so that the lysimeter must satisfy the condition dD dt + dS dt > 0. Therefore, the following Equation (4) was used:

Calculation of the Infiltration Rate of Water into the Soil
The water infiltration has a downward vertical behavior in the lysimeter [26], and conditions of free drainage [7]. When the soil moisture was below field capacity water content, the increments of soil water stored in the soil vessel represent the soil infiltration rate dI dt without any drainage component dD dt . The infiltration rate was then represented with the following Equation (5) which was called Approach 1 [45]: However, once the soil reached the water content at field capacity, it is no longer able to store more water and thus, the incoming water will be drained [45]. In that way, the sum of soil water storage and drainage represents the water infiltration rate of the soil profile. The infiltration rate model is represented as the Equation (6) which was called Approach 2: Both approaches are valid if the water content at soil saturation is not reached due to the rainfall rate not overpassing the soil basic infiltration capacity. In this experiment the soil never overpass the saturation capacity, therefore there was no surface runoff. This complies with the principle of mass conservation and continuity equation.
The resulting flows and infiltration rate are representing the area of the lysimeter (0.54 m 2 ), so with the assumption that 1 L of water is equal to 1 kg of water mass and that 1 mm is equal to 1 L/m 2 [27], 1 kg/min of lysimeter is equal to 1.85 kg/min or 1.85 mm/min of 1 m 2 of surface.
The irrigation events were not analyzed because their time and volume of application were small, so they do not allow a visualization of the infiltration capacity curve. The infiltration rate was constant throughout the two or three hours that the irrigations lasted.

Estimation of Soil Moisture Content
The soil profile water content inside the soil vessel was determined by the gravimetric method, since the mass of the wet soil profile was known at any time. Once obtained, the gravimetric water content was converted to volumetric (θ v ) as in Equation (7): where θ v is the volumetric water content is of the soil at instant i (m 3 /m 3 ), ρ b is the bulk density of the soil, m i is the wet mass of the soil inside the soil vessel at instant i (g or kg) and m dry is the dry mass of the soil inside the soil vessel. This last value was determined from the apparent density of the soil and the total volume of the soil vessel (0.161 m 3 ) [13]. The moisture content at field capacity (θ FC ) was established two or three days after the rain events, when the soil stopped draining and the soil vessel value was constant, which indicates that the soil is already capable of retaining water [4]. No runoff was observed in any of the rain events. Evapotranspiration ( dET dt ) was calculated because it affects the redistribution of soil moisture, which was obtained using the mass balance method (Equation (1)), remembering that evapotranspiration and rain cannot occur at the same time. The mass at instant i represents soil-water aggregates, and the mass of the vessel structure is considered zero as it was tared to zero before inserting the soil profile into the soil vessel.

Model Calibrations
The results of the water infiltration rate dI dt into the soil were used to calibrate empirical and physics infiltration models from the previously mentioned authors.
One of the most used models to represent the infiltration rate is the Kostiakov's model, which can be generally represented as Equation (8): where t is the infiltration time and a and b are empirical constants [4,6,[46][47][48]. This equation has a boundary condition: lim t→∞ dI dt = 0 [49,50]. All empirical parameters of all infiltration equations were determined by an iterative fit of the experimental data obtained with our approaches and selected using the best root mean square error (RMSE) in the fitting [46].
Parameter b of Kostiakov's equation is limited to being greater than zero and less than one. Therefore, the initial value of the infiltration rate tends to be infinite and has a basic rate to zero in an infinite time. In long irrigation events, Lewis [47] observed that the infiltration rate tends to a constant value before the end of the irrigation. So, in these long events, Kostiakov's equation usually overestimates the infiltration [51,52]. Therefore, a constant was added to this equation to represent the basic infiltration rate, which resulted in the modified Kostiakov's equation: where c is an empirical constant, which represents the basic infiltration rate [6,53,54]. This equation also is known as the Kostiakov-Lewis or Mezencev model. Its boundary condition is: lim . Philip's equation for the homogeneous porous medium infiltration rate is a simplification of Richard's model. The first term of the equation establishes the absorption of water into the soil and the second term represents a truncated series of small flooding times (A), as the Equation (10): where S is the sorptivity [50,53] and lim t→∞ dI dt = A, where A is a constant value [49]. Horton establishes an equation of infiltration capacity based on the behavior of nature that follows a law of inverse exponential. The infiltration rate decreases to a constant value, with an exponential decrease of change, and proportional to the infiltrated water volume as follows: (11) where V f is the final infiltration rate or basic rate, V 0 is the initial infiltration rate, and β is an empirical constant of change from the initial rate to the basic rate [17]. The boundary condition is lim t→∞ dI dt ≥ 0 [50]. When comparing the four models, it can be seen that the parameters of Kostiakov's model have no physical meaning and for longer times of water infiltration, it loses its representativeness and accuracy of the rate of water infiltration into the soil. For t = 0, the infiltration rate remains undetermined. In the modified Kostiakov's model, a term was added to represent the basic rate observed at long times of water infiltration. In the case of the Horton model, the initial infiltration rate is defined by a finite value and the model describes the behavior of the soil infiltration capacity curve. Finally, Philip described the rate of infiltration into the soil surface when there is a small amount of ponding, so its terms do have physical significance for the behavior of the soil when water is allowed to enter.
The statistics used to evaluate the accuracy of the estimation of the water infiltration rate into the soil using the weighing lysimeter and using the above-mentioned models were: the standard deviation (s), the root mean square error (RMSE), and the coefficient of determination (R 2 ) [55]. Figure 3 shows the mass curves of the two vessels of the lysimeter in the analyzed time interval. The initial data reported by the soil vessel indicates the mass of the soil particles plus the water content. Consequently, the mass variations reported after this initial point indicate the variations of the water content within the soil. In the curve of the soil vessel, the inputs were the rain and the irrigation and they can be differentiated by the type of increase they present. Irrigation has a very vertiginous increase in a short time, while the rain has a discontinuous growth in time, on the other hand, the descents can represent the evapotranspiration losses, drainage losses or both.

Results and Discussion
of determination (R ) [55]. Figure 3 shows the mass curves of the two vessels of the lysimeter in the analyzed time interval. The initial data reported by the soil vessel indicates the mass of the soil particles plus the water content. Consequently, the mass variations reported after this initial point indicate the variations of the water content within the soil. In the curve of the soil vessel, the inputs were the rain and the irrigation and they can be differentiated by the type of increase they present. Irrigation has a very vertiginous increase in a short time, while the rain has a discontinuous growth in time, on the other hand, the descents can represent the evapotranspiration losses, drainage losses or both.

Results and Discussion
In the case of the data from the drainage vessel (red curve), the rise of the curve represents that the drainage has started and the sudden descents show us the emptying of the vessel. As shown in the Figure 3, only two drainage events were registered by the vessel on 13-14 March and 28-30 April. The drainage vessel was programmed to open its valve at 3000 g of storage and discharge (to avoid damaging the load cell). Therefore, the drainage vessel curve in the Figure 3 shows accumulated masses less than 4000 g. The emptying is almost instantaneous, it takes less than 3 min and, as mentioned, no drop is lost. The March event caused more drainage than the April event, in less time of rain. This indicates that the rain intensity of the March event was higher than in the April event, ass more water fell in less time. The seven rainfall events between February and March contributed to the first drainage event occurred in 14 March. Meanwhile, the five rainfall events occurred between April and May caused the second drainage event in 28-30 April.  In the case of the data from the drainage vessel (red curve), the rise of the curve represents that the drainage has started and the sudden descents show us the emptying of the vessel. As shown in the Figure 3, only two drainage events were registered by the vessel on 13-14 March and 28-30 April. The drainage vessel was programmed to open its valve at 3000 g of storage and discharge (to avoid damaging the load cell). Therefore, the drainage vessel curve in the Figure 3 shows accumulated masses less than 4000 g. The emptying is almost instantaneous, it takes less than 3 min and, as mentioned, no drop is lost. The March event caused more drainage than the April event, in less time of rain. This indicates that the rain intensity of the March event was higher than in the April event, ass more water fell in less time. The seven rainfall events between February and March contributed to the first drainage event occurred in 14 March. Meanwhile, the five rainfall events occurred between April and May caused the second drainage event in 28-30 April.

Estimation of Rain Inflow
The rainy days total of 12 were analyzed and were determined with the Equation (4), representing 12 events. The Table 1 shows the rainfall recorded by day in the weighing lysimeter and those observed in a rain gauge installed near the weighing lysimeter at 800 m. The rain calculated with the lysimeter for the first analysis day concerning the Instituto Técnico Agronómico Provincial (ITAP, Spanish acronym) station has a difference of 0.05 mm. The differences found can be explained by the space behavior of the rainfall variability in a region, which has a greater impact in semi-arid zones such as the study zone [56]. There is an absolute average difference of 1.59 mm between the two devices.

Estimation of Soil Moisture Content
The soil moisture content was obtained using the gravimetric method (Equation (7)). Two drainage events were caused with the rain events of 13-14 March, 28-30 April so only with these events was possible to determine the soil field capacity (θ FC ) based on its definition. The Figure 4 shows the rain behavior and its relation to the moisture content of the soil, which shows that when rain is falling, there is a mass increase, and when the soil reaches the field capacity and the rain continues, there is drainage; when there is a fall in the humidity curve, drainage increases.
By plotting the water content curve with that of the drainage vessel (red curve) it was possible to determine the moment when the moisture content at field capacity of the soil (θ FC ) was reached. In both months, the θ FC was established when the drainage vessel curve became constant, indicating that the soil had stopped draining and there was no variation in the weighing lysimeter.
A horizontal line was drawn on the water content at the instant the mass drainage vessel was not increased, indicating that the θ FC was reached and was established two days after the beginning of the rain in 0.33 cm 3 /cm 3 .
Since the θ FC was established days after the rain events, evapotranspiration dET dt was presents, so it was calculated as indicated by [27,32,43] in their work. Table 2 has the water balance of the days involved in the calculation of the θ FC . dET dt is a variable that subtracts moisture from the soil surface, this variable for the first days of each event only represents 0% of water losses and rainfall; it increases as the rain decreases over the days, so that in the second day dET dt represents 15.47% and 10.32% of water losses and 18.96% and 7.83% of rainfall. The total dET dt losses for all days of each event used for the θ FC calculation represent 11.04% for the first (Figure 4a) and 24.37% for the second (Figure 4b).  Since the was established days after the rain events, evapotranspiration ( ⁄ ) was presents, so it was calculated as indicated by [27,32,43] in their work. Table   2 has the water balance of the days involved in the calculation of the . ⁄ is a variable that subtracts moisture from the soil surface, this variable for the first days of each event only represents 0% of water losses and rainfall; it increases as the rain decreases over the days, so that in the second day ⁄ represents 15.47% and 10.32% of water losses and 18.96% and 7.83% of rainfall. The total ⁄ losses for all days of each event used for the calculation represent 11.04% for the first (Figure 4a) and 24.37% for the second (Figure 4b).
⁄ represents 8.28% of the rain of the first event ( Figure 4a) and 35.38% of the second event (Figure 4b). The Table 3 reports the moisture content at field capacity values obtained and some that report the literature for a silty loam soil. The values with less standard deviation with the result of the weighing lysimeter were the value obtained in the laboratory in this investigation and the maximum value reported by Assi et al. [57] with 0.014 cm 3 /cm 3 . These are followed by the values reported by FAO [58] in graphs and tables based on soil texture to obtain the hydrodynamic characteristics of the soil (θsaturation, θFC and θwilting point) with 0.02 cm 3 /cm 3 , which shows a good agreement with results obtained in laboratory and lysimeter.  The Table 3 reports the moisture content at field capacity values obtained and some that report the literature for a silty loam soil. The values with less standard deviation with the result of the weighing lysimeter were the value obtained in the laboratory in this investigation and the maximum value reported by Assi et al. [57] with 0.014 cm 3 /cm 3 . These are followed by the values reported by FAO [58] in graphs and tables based on soil texture to obtain the hydrodynamic characteristics of the soil (θ saturation , θ FC and θ wilting point ) with 0.02 cm 3 /cm 3 , which shows a good agreement with results obtained in laboratory and lysimeter.

Water Infiltration Rate
Seven rains were analyzed with the Approach 1. The soil was below field capacity so that the positive variations reported soil vessel indicated the water infiltration rate (13,18,19 and 24 February,27 April and 10 and 30 May). According to USDA [60] and Evanylo and McGuinn [61] the best way to determine the rate or speed of infiltration is when the soil is close to field capacity; four rain events had drainage due to soil was above field capacity; in these rains, approach 2 was applied (13-14 March and 28, 29 and 30 April). Figure 5 shows the 11 scatter soil water infiltration rate diagrams of the rain 12 events analyzed, most of the graphs have a duration of less than 24 h. The soil shows a very similar behavior in all the graphs, it starts with lower rates until it reached a maximum peak and then decreases to a stable or basic rate. The basic infiltration rate was obtained considering the eleven rainfall events in 1.49 mm/h with a 0.36 mm/h standard deviation. According to USDA [60], this rate occurs when the soil is almost saturated and it does not decrease or increase as more water is added.

Water Infiltration Rate
Seven rains were analyzed with the Approach 1. The soil was below field capacity so that the positive variations reported soil vessel indicated the water infiltration rate (13,18,19 and 24 February,27 April and 10 and 30 May). According to USDA [60] and Evanylo and McGuinn [61] the best way to determine the rate or speed of infiltration is when the soil is close to field capacity; four rain events had drainage due to soil was above field capacity; in these rains, approach 2 was applied (13)(14)29 and 30 April). Figure 5 shows the 11 scatter soil water infiltration rate diagrams of the rain 12 events analyzed, most of the graphs have a duration of less than 24 h. The soil shows a very similar behavior in all the graphs, it starts with lower rates until it reached a maximum peak and then decreases to a stable or basic rate. The basic infiltration rate was obtained considering the eleven rainfall events in 1.49 mm/h with a 0.36 mm/h standard deviation. According to USDA [60], this rate occurs when the soil is almost saturated and it does not decrease or increase as more water is added. The variability of the standard deviation could be due the experimental error. The infiltration rate or speed modifications reaching a 23%, could probably be due the flow paths alterations at forcing the water to pass through an orifice. Further research is needed in order to minimize this variability, i.e., a number of orifices placed in the horizontal plane The variability of the standard deviation could be due the experimental error. The infiltration rate or speed modifications reaching a 23%, could probably be due the flow paths alterations at forcing the water to pass through an orifice. Further research is needed in order to minimize this variability, i.e., a number of orifices placed in the horizontal plane in the soil vessel bottom could represent a variability of the measured infiltration rate/speed minimizing the measurement errors. Nevertheless, our methodology represents an approach to estimate the infiltration rate and could be compared with other methodologies.
The variability of the initial infiltration rate may be due to rainfall intensity and initial moisture content [62,63]. The initial moisture contents of the seven infiltration rate scatter diagrams below field capacity ranged from 0.23 to 0.30 cm 3 /cm 3 and the other four scatter diagrams were 0.30 to 0.33 cm 3 /cm 3 .
The basic infiltration found in both approaches (summarized in the Table 4) were low compared to the values reported in the literature. However, the values found in this investigation coincide with those reported by Li et al. [64] for homogeneous wetting and their value was obtained with an initial soil water content of 0.034 m 3 /m 3 , which was air-dried and placed in a 50 cm column. The lower value obtained in Li et al. [64] was with a slightly water-repellent silty loam soil (this is achieved by adding octadecylamine to the soil). Concentric cylinders [6] The basic infiltration rate obtained with the weighing lysimeter, was lower compared with the research works of Maldonado [65] and Ali [6] considering the same soil type. Some factors that could have caused the low value of the basic infiltration rate are organic matter content, vegetative cover, and the crop root [6,63]. Nevertheless, further investigation is needed to corroborate it.
The results obtained with the two proposed approaches are simple and representative of the current capacity of infiltration of the soil, where roots have been growing within the soil profile, fertilizers have been added, and the profile has been compacted by the same force exerted by rainwater and irrigation. All of this is not reflected in the values seen in Table 4.

Model Calibrations
The 11 infiltration curves were calibrated with the four models described above, from its highest point. In Figure 6, only the four models with longer rainfall duration are presented. The model that showed the best fit was Horton's, with a RMSE of 0.89 ± 0.40 and R 2 of 0.65 ± 0.30, followed by Kostiakov with a RMSE of 1.06 ± 0.40 and R 2 of 0.63 ± 0.23, modified Kostiakov with a RMSE of 1.03 ± 0.43 and R 2 of 0.62 ± 0.22, and Philip with a RMSE of 1.30 ± 0.42 and R 2 of 0.48 ± 0.17. The only model that established a lower initial infiltration rate was Horton.
The models establish the value for the basic infiltration rate in some of its variables in the equation. In the case of Kostiakov, it does not have an established value. Modified Kostiakov established the basic infiltration at 0.14 ± 0.14 mm/h (in its constant c) [53,54], Philip at 0.16 ± 0.11 mm/h (in A), and Horton at 0.38 ± 0.34 mm/h (in V f ) [17,50]. The models calibration showed a behavior as indicated by their boundary limits, as the infiltration rate tends to constant values close to zero when time tends to infinity [49,50]. In the research works of Rodríguez-Vásquez et al., [66] and Mirzaee et al., [67] measured infiltration rate with double ring infiltrometer apparatus for silt loam and clayey soil, their parameter values obtained for the different infiltration models are shown in the Table 5. In the first work, the water contents were between 24.3% and 31.1%, while in the second article this was not mentioned, and our range was of 0.23-0.33 cm 3 /cm 3 . In the research works of Rodríguez-Vásquez et al., [66] and Mirzaee et al., [67] measured infiltration rate with double ring infiltrometer apparatus for silt loam and clayey soil, their parameter values obtained for the different infiltration models are shown in the Table 5. In the first work, the water contents were between 24.3% and 31.1%, while in the second article this was not mentioned, and our range was of 0.23-0.33 cm 3 /cm 3 . The fit parameters of Kostiakov model in this investigation compared with the Rodríguez-Vásquez et al., [66], have a difference of 39% for the parameter a, while for the parameter b is of 20%. In the Philip model, there is a difference of 35% for S parameter and the value of A parameter of Rodríguez-Vásquez et al., [66] is 25 times higher than ours. This difference probably is due to that error of estimation is greater in the concentric cylinders.
There is a significant difference with respect to Mirzaee et al., [67] work in the four models varying from 68 to 245%. because their reported values are of a soil profile with siltyloam and clayey horizons, which the infiltration capacity is very low due to characteristics of the clays.

Conclusions
In this research, the mass values of a compact weighing lysimeter were used to measure soil moisture behavior and the water infiltration rate of a soil silt loam under unsteady rainfall conditions, using two approaches, the first part when the soil is below and the second when the soil is above field capacity. The results showed that with the two approaches it is possible to establish the soil infiltration curve and for short and long application times in real conditions. The soil is subject to a variable and faster infiltration rate at the beginning, but it reaches a basic or stable infiltration as the soil fills its pores.
These approaches are simple but satisfy the conservation of mass and continuity equation for incomprehensible flows. Although they are used for non-saturation conditions, they can be used for saturation conditions by adding the surface runoff variable that would subtract incoming flow to the soil.
The use of the weighing lysimeter with these approaches allowed to characterize the infiltration rate of a homogeneous and non-homogeneous soil at depths where agriculture regularly takes place. However, the limitation is that it is not possible to monitor the progress of the wet front, so further research is necessary in order to have greater measurement instruments (i.e., soil moisture sensors) to characterize the water movement in the soil and the forces acting within it (i.e., soil water potential sensor).