Next Article in Journal
Minimal Building Flood Fragility and Loss Function Portfolio for Resilience Analysis at the Community Level
Next Article in Special Issue
Sustainable Intensification of Agriculture in the Context of the COVID-19 Pandemic: Prospects for the Future
Previous Article in Journal
Microclimate Characteristics and Evapotranspiration Estimates of Cucumber Plants in a Newly Developed Sunken Solar Greenhouse
Previous Article in Special Issue
Effect of Distribution Uniformity and Fertigation Volume on the Bio-Productivity of the Greenhouse Zucchini Crop
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Pulse Drip Irrigation Duration on Water Distribution Uniformity

1
Andalusian Institute of Agricultural and Fisheries Research and Training (IFAPA), Junta de Andalucía, “Alameda del Obispo”, 14004 Córdoba, Spain
2
Andalusian Institute of Agricultural and Fisheries Research and Training (IFAPA), Junta de Andalucía, “La Mojonera”, 04745 Almería, Spain
*
Author to whom correspondence should be addressed.
Water 2020, 12(8), 2276; https://doi.org/10.3390/w12082276
Submission received: 9 July 2020 / Revised: 2 August 2020 / Accepted: 9 August 2020 / Published: 13 August 2020
(This article belongs to the Special Issue Sustainable Irrigation Management in Agriculture)

Abstract

:
Developing an appropriate irrigation schedule is essential in order to save water while at the same time maintaining high crop yields. The standard procedures of the field evaluation of distribution uniformity do not take into account the effects of the filling and emptying phases of the irrigation system. We hypothesized that, in sloping sandy soils, when short drip irrigation pulses are applied it is important to take into account the total water applied from the beginning of irrigation until the emptying of the irrigation system. To compute distribution uniformity, we sought to characterize the filling, stable pressure, and emptying phases of a standard strawberry irrigation system. We found that the shorter the time of the irrigation pulse, the worse the distribution uniformity and the potential application efficiency or zero deficit are. This effect occurs because as the volume of water applied during filling and emptying phases increases, the values of the irrigation performance indicators decrease. Including filling and emptying phases as causes of non-uniformity has practical implications for the management of drip irrigation systems in sloping sandy soils.

1. Introduction

The relentless growth of the human population, the rise in general living standards, and the need for water to maintain the long-term sustainability of the environment and ecosystems are straining water resources all over the world [1]. There is a growing trend in water withdrawal for agriculture, which has a critical importance in semi-arid Mediterranean environments, such as Spain [2].
In Spain, horticultural crops with a high added value are highly dependent on irrigation. Berries from the province of Huelva are a clear example of this case. The annual export of berries is around one billion euros, which generates approximately 80,000 work positions. Seventy-five percent of the berry production in Huelva is situated near the Doñana National Park, which is catalogued as a UNESCO World Heritage site [3]. The main source of water for berry fields and for the Doñana National Park system is the same, the “Almonte-Marismas” aquifer, which has been declared to be in a “poor” state by water authorities [4]. Due to this, berry production in Huelva is criticized for negatively affecting the Doñana National Park because of poor irrigation management.
Berry production in Huelva has some specific drawbacks that hinder efficient irrigation management. One of the main irrigation drawbacks is having very sandy soils (more than 90% sand) in most of the berry production area [5]. Sandy soils not only have a lower water-holding capacity but also significant downward water movement [6,7]. Farmers seek to maximize berry production, which requires an irrigation strategy without water deficient crops. This approach may result in water and nutrient losses and groundwater pollution due to over irrigation in sandy soils [8,9]. In fact, irrigation frequency was found to affect nitrate leaching beneath the root zone [10,11,12].
A significant number of farmers in Huelva apply irrigation pulses of 30–40 min. [13] recommended different irrigation pulse durations for the strawberry producers of Huelva depending on the time of the season. During the first part of the season (October to January), [13] suggested that the pulse duration should range from 10 to 20 min. In the second half of the season, from February to June, irrigation pulses should be increased from 20 up to 60 min, depending on the root depth and water needs. On the other hand, some works revealed that pulse drip irrigation can improve strawberry production and marketable yields relative to the no pulse irrigation [14,15]. Therefore, developing appropriate irrigation pulse scheduling is essential to reduce current irrigation doses and reach high quantitative and qualitative production levels.
An additional drawback in high-frequency irrigation (short pulses) is water distribution on sloping plots. The effect of the slope on irrigation uniformity has also been addressed. [16] developed a method to design microirrigation laterals in uniformly and non-uniformly sloped fields. [17] developed analytical expressions relating water distribution indices to design variables that define trapezoidal drip irrigation plots. [18] evaluated the effect of hydraulic head and slope on water distribution uniformity. [19] analyzed the influence of the slopes on the water distribution of a single pipe during the stable pressure phase of irrigation. [20] developed an analytical method to design evenly sloped drip irrigation systems. It must be noted that all these studies assume stable pressures throughout the irrigation pulse, which may not be entirely true for short irrigation pulses on sloped fields.
On sloping plots, an uneven discharge of the irrigation tapes considering filling, stable pressure, and the emptying of the pipe flow process have been described by some authors [21,22]. In fact, these studies have analyzed the effects of the slope on the water emitter application at lateral levels according to the emitter position. Their results have been a starting point in understanding the effect of filling and emptying in a single tape in short pulses. However, these effects have not been studied at the irrigation unit scale for commercial fields.
The most popular methodologies for the field evaluation of distribution uniformity were proposed by [23,24,25]. These three methodologies differ in the number and location of the emitters selected to represent the entire irrigation unit.
Merriam and Keller’s methodology (adopted by FAO) is based on the measurements of the water discharge located in 16 locations. These 16 locations are selected from four laterals and four emitters per lateral located at the inlet, at a third and at two thirds of the distance to the end of the manifold and the lateral, and also at the end of the manifold and the lateral.
ASAE’s method is based on the works [26,27,28]. These works presented a statistical approach that uses water measured at 18 locations. The locations are selected randomly to measure both the discharge and pressure.
Burt’s method indicated that discharge and pressure should be measured from a total of 66 locations: 16 locations close to the water source, 16 locations near the middle of the field, and 28 locations at the end of the lateral most distantly to the manifold.
Merriam and Keller’s, ASAE’s, and Burt’s methodologies calculate the emitter’s discharges assuming that the pressures of the irrigation system are stable throughout the irrigation pulse. Merriam and Keller’s and Burt’s methods record water flow for an integer number of minutes (from 3 to 6 min). Merriam and Keller’s method uses a water volume between 100 and 250 mL for each emission point tested. On the other hand, ASAE’s method records the time required to fill a 100 mL container. Therefore, in these methods the distribution uniformity value is only representative of the stable pressure of the pipe flow process.
Burt’s methodology includes the emptying phase of the irrigation system as another component of the distribution uniformity, but indicates that a longer data collection time is not guaranteed because of the relatively small impact of the emptying phase of the distribution system on the overall distribution uniformity.
Standard irrigation units without anti-drain emitters have three pipe flow phases per irrigation pulse. The first phase is the filling phase. This starts at the beginning of the irrigation pulse and ends when the entire irrigation unit has stable pressure in all the emitters. The stable pressure phase starts once there is stable pressure and ends at the end of the irrigation pulse. The emptying phase starts at the end of the irrigation pulse and ends after complete water emptying of the irrigation system.
Our hypothesis is that, in high-frequency and short-duration irrigation and sloped fields, it is important to compute the water volume applied during the filling, stable pressure, and emptying phase of the pipe flow process to know the real water distribution uniformity. To the best of our knowledge, there are no studies that characterize the different phases of the pipe flow process in a drip irrigation unit of a real field.
The authors are aware that the use of different tapes or emitters could change the quantitative response of the filling, stable pressure, and emptying phase of the pipe flow process. This manuscript shows the influence of these three pipe flow phases on short-term irrigation in a specific case study. In this way, the authors understand the limitations of this work. However, we understand that this manuscript provides a step forward to manage short irrigation pulses on drip irrigation systems under more arbitrary irrigation material, soil, and slope conditions.
Therefore, the objectives of this work are: (i) to characterize the filling, stable pressure, and emptying phases of the pipe flow process in high-frequency and short-duration irrigation on a sloping field, and (ii) to evaluate the effect of duration of irrigation pulses on the water distribution uniformity and the application efficiency of a drip irrigation unit.

2. Materials and Methods

2.1. Description of the Study Irrigation Unit

The study was carried out during the 2015–2016 irrigation season in a commercial strawberry production farm located near Lepe, Huelva, Spain. The climate in Huelva is subhumid Mediterranean, classified as warm temperate with a maritime influence. Huelva has more than 3000 light hours per year and temperatures that only exceptionally fall below 0 °C. Therefore, frosts are unusual. In Huelva, the annual average rainfall is 553 mm, and the average, maximum and minimum temperatures are 17.4, 28.8, and 8 °C, respectively. The soils of berry farms located in Huelva are sandy (quite a few soils are more than 90% sand).
The study irrigation unit was composed of 12 strawberry greenhouses. Each greenhouse has an area of 290.42 m2 and consists of six rows 1.1 m separated. The area of the irrigation unit was 3485 m2 (79.2 m wide and 44 m long). The irrigation unit has an average ascending transverse slope of 0.3% and an average descending longitudinal slope of 3% (Figure 1). The slopes of the irrigation units were obtained by a topographic survey.
The irrigation unit had a manifold polyethylene pipe and 72 built in emitter laterals. The manifold pipe and built in emitter laterals are new every season. Thus, the operating conditions are similar throughout the entire irrigation season in the case study of this work. The manifold polyethylene pipe was 63 mm in diameter and 79.2 m in length. A total of 72 laterals (irrigation tapes) were inserted into the manifold pipe. Each lateral had a 16.2 mm inner diameter and a 44 m length. The laterals used were Netafim StreamlineTM 16,080, with a nominal flow and pressure of 5.25 l h−1 m−1 and 1 bar, respectively. The drip emitter separation was 0.2 m. The total number of drip emitters in the irrigation unit was 15,912.
The water applied and pressure at the head of the unit were measured rigorously and controlled during this study. The volume of irrigation applied was measured through a flow-meter located at the head of the irrigation unit. The inlet water pressure was regulated at the entrance of the irrigation unit using a pressure regulator. The pressure regulator was located after the flow-meter to avoid the head loss caused by the flow-meter. The inlet water pressure was fixed to 1 bar, which is the nominal pressure of the drip emitters.

2.2. Location of Measuring Points

Six laterals were selected for the irrigation unit studied (Figure 1). Four of them (L1, L2, L3, L4) correspond to the criteria proposed by Merriam and Keller’s methodology: Laterals L1, L2, L3, and L4 are located at the beginning, at a third, at two thirds of the distance to the end of the manifold pipe, and at the end of the manifold pipe, respectively. The location of the other two irrigation points was at a quarter (L1A) and at a half (L1B) of the L1 to L2 distance. This decision was based on a total applied volume gradient that was studied previously for a single pipe [21].
Sixty measuring points were used for studying the water distribution in the study irrigation unit. Ten emitters were used for each lateral. The location of the emitters was also based on a previous study of water distribution in a single sloped lateral [21]. From the 60 volume measuring points selected, 16 correspond to the criteria proposed by Merriam and Keller’s methodology. From now on, these 16 locations will be referred to as the 16 MK locations (Figure 1).
The pressure was measured at the 16 MK locations. To measure the water pressure, 16 manometer intakes were installed in the laterals. Sixteen glycerine manometers were located inside the manometer intakes close to the 16 MK selected locations. The manometers had a measuring range from 0 to 2.5 bar and a 0.05 bar accuracy (EN 837-1/6, weinmann-schanz©).

2.3. Characterization of the Pipe Flow Process

We characterized the three phases of the pipe flow process (filling, stable pressure, and emptying) in the irrigation unit studied. The water distribution uniformity of the study irrigation unit was evaluated at the beginning of the season (November). Firstly, the water distribution uniformity was evaluated by measuring the pressure during the stable pressure process and transforming it into volume at each of the 16 MK measuring points. Therefore, through this form of calculation, the water distribution uniformity is independent of the irrigation duration. Secondly, four different irrigation pulses of 5, 10, 15, and 20 min were evaluated considering the volume applied by the tested emitters from the beginning of the irrigation to the emptying of the pipe. In this work, we do not evaluate the effect of the same total irrigation time with a different number of irrigation pulses. The objective was to evaluate each irrigation pulse duration independently. For that purpose, each irrigation pulse (5, 10, 15, 20 min) was evaluated on different days to assure the same initial conditions.

2.3.1. Filling Phase

The filling phase was characterized by measuring the continuous evolution of the pressure at 16 MK locations. The 5 min and 10 min irrigation pulses were used to measure the evolution of pressure. Five minutes of irrigation pulses was enough to characterize the entire filling phase. To measure the pressure, we used 16 manometers simultaneously (with a measuring range from 0 to 2.5 bar). The manometer accuracy was 0.05 bar. The evolution of the pressure in each of the 16 manometers was recorded with 16 Gopro© video cameras operating simultaneously. The pressure was measured every 10 s and recorded. The pressure values used each 10 s for the characterization of the filling phase were the average values of the 5 min and 10 min irrigation pulses.
The emitter flow rate during the filling phase was calculated using the equation supplied by the drip emitter manufacturer. The manufacturer’s equation was tested before use (data not shown). The manufacturer’s equation was used in each one of the 16 MK points:
Qf ej = 0.3824 Pf ej 0.4384
where Qfej (l h−1) is the flow rate (liters per hour) of the emitter e at the instant j during the filling phase. Pfej is the pressure (meter’s water column) measured in the manometer intake e at the instant j during the filling phase. In this equation, the values of e range from 1 to 16, corresponding to the 16 MK measuring points (Figure 1). The values of j range from 1 to 30 because there are six pressure values per minute (1 each 10 s), and because the duration of the pressure recording of the irrigation pulse was 5 min.
During the filling phase, the volume of water applied to each of the 16 MK emitters was estimated using the following equation:
Vf e = j = 1 j = T f Qf ej Δ T f
where Vfe (l) is the water volume applied by the emitter e (1 to 16 MK) during the filling phase. ΔTf (h) is the time interval of the pressure recorded (10 s). Tf is the number of the first time interval measured once the pressure is stable in the 16 MK emitters.
From the water volume estimated during the filling phase (Vf1,… Vf16), an interpolation was performed to assign the volumes applied during the filling phase to the 15,912 emitters of the irrigation unit. Third degree polynomial equations were obtained based on the applied volume by 16 MK emitters during the emptying phase.

2.3.2. Stable Pressure Phase

The stable pressure phase was characterized by measuring the pressure twice at each of the 16 MK measuring points for each irrigation pulse (5, 10, 15, and 20 min). The pressure at the entrance of the irrigation unit was set at 1 bar for all the irrigation pulses. At each of the 16 MK points, the pressure was measured twice for each irrigation pulse. The value of the pressure assigned to each measurement point was the average of the two values. To measure the pressure, intake manometers were installed on the 16 MK measuring points (Figure 1). The flow rate was calculated using the equation proposed by the drip emitter manufacturer:
Qsp e = 0.3824 Psp eavg 0.4384
where Qspe (l h−1) is the flow rate of the emitter e during the stable pressure phase, and P s p e a v g is the average pressure measured (meter’s water column) in the intake manometer e during the stable pressure phase in the four irrigation pulses.
The volume of water applied during the stable pressure phase in each of the 16 MK emitters was estimated using the following equation:
Vsp e = Qsp e T sp ,
where Vspe (l) is the water volume applied by the emitter e (1 to 16 MK) during the stable pressure phase, and Tsp (h) is the duration of the stable pressure phase.

2.3.3. Emptying Phase

The emptying phase was characterized by measuring the water volume applied by 60 selected drip emitters (Figure 1) after the end of the irrigation pulse duration. Based on the experience of the authors, to ensure the complete emptying of the laterals the emptying phase was considered to be 2 h. After 2 h, the authors can guarantee that no dripper emits water. The data of the emptying volume was the average value of four independent irrigations on different days. The water of 60 drip emitters was stored in 0.7 l plastic containers. These plastic containers were placed under the 60 drip emitters after the end of the irrigation pulse. The emptying water volumes were measured using a graduated cylinder of 0.5 l capacity with a 5 mL accuracy. Third degree polynomial equations were developed based on the relationships between the volume applied by 60 emitters that were controlled during the emptying phase and the distance from the beginning of laterals.

2.4. Validation of the Water Distribution Model

An empirical water distribution model was developed to simulate the flow emitted by each one of the 15,912 emitters for the study irrigation unit. The model was based on the previous characterization of the filling, stable pressure, and emptying phases of the pipe flow process. To validate the water distribution model, four irrigation pulses were applied on four different days in May. The duration of the irrigation pulses was 5, 10, 15, and 20 min. The pressure at the entrance of the irrigation unit was also set at 1 bar in all the irrigation pulses.
To validate the model during the three irrigation phases (filling, stable pressure, and emptying), the water applied by the 60 selected emitters (Figure 1) was measured. The filling and stable pressure volumes were measured together as a unique stage called the volume of irrigation pulses. The water volumes at the emptying phase were measured separately. The pressure measurement methodology to validate the stable pressure phase was the same used during the characterization stage in the 16 MK measuring points. The values measured were compared to the values predicted by the model.
The water applied by emitters was stored in plastic containers, both from the beginning until the end of the irrigation pulse and from the end of the irrigation pulse until 2 h after the end of the irrigation pulse. Subsequently, the irrigation pulse water volumes were measured using a graduated cylinder of 0.5 l capacity with a 5 mL accuracy.
In the validation stage, 120 plastic containers of 0.7 l were used; 60 of them were used for measuring the volume of irrigation pulse (filling + stable pressure phases), and 60 were used for measuring the water volume of the emptying phase. The water volume was measured using a graduated cylinder of 0.5 l capacity with a 5 mL accuracy.

2.5. Irrigation Performance Indicators

2.5.1. Distribution Uniformity

Distribution uniformity is usually defined as the ratio of the smallest accumulated discharges in the irrigation water distribution to the average discharges of the whole distribution [29]. The distribution uniformity of the low quarter (UDlq) is an irrigation performance indicator developed by the United States Department of Agriculture–Natural Resources Conservation Service (USDA-NRCS). UDlq has been widely accepted by the international irrigation research community for decades [30,31], and there is a high correlation between distribution uniformity and other water distribution indicators [32,33]. UDlq is calculated using Equation (5):
UD lq = 100 Q 25 % Q avg = 100 V 25 % V avg
where Q25% is a 25% average flow (l h−1) with a lower flow of water, Qavg is the average flow (l h−1), V25% is the 25% average water applied volume (l) with a lower volume of water, and Vavg is the average water applied volume (l).
The performance indicator UDlq was calculated both with values measured in the validation stage or with values predicted using the water distribution model generated in this study.
  • Standard distribution uniformity:
Distribution uniformity in the stable pressure phase (UDlq_sp) is the water distribution uniformity indicator obtained using classical evaluation methodologies (i.e., Merriam and Keller’s, ASAE’s, and Burt’s). It is important to highlight that UDlq_sp is independent of the duration of the irrigation pulse. UDlq_sp was calculated using Equation (3), and Merriam and Keller’s methodology using the following equation:
UD lq sp = 100 Q 25 % sp Q avg sp
where Q25%_sp (l h−1) is the average flow of 4 MK emitters with a lower water flow during the stable pressure phase, and Qavg_sp (l h−1) is the average of the flow of the 16 MK emitters during the stable pressure phase.
  • Distribution uniformity of the total water applied in the irrigation unit studied:
The distribution uniformity of the total water applied in the field (UDlq_field_3phases) depends on the filling, stable pressure, and emptying phases of drip irrigation. UDlq_field_3phases was calculated for the water volume measured in irrigation pulses of 5, 10, 15, and 20 min in the validation stage. The volume of water applied during these three drip irrigation phases was measured in the 16 MK emitters from the beginning of the irrigation pulse until the emptying of the irrigation system. The duration of the irrigation tests was 2 h of the emptying phase plus the time of the irrigation pulses (5, 10, 15, and 20 min). The UDlq_field_3phases was calculated using the following equation:
UD lq field _ 3 phases = 100 V 25 % field _ 3 phases V avg field _ 3 phases
where V25%_field_3phases (l) is the average of 4 MK emitters with a lower volume of total water applied by the 16 MK emitters from the beginning of the irrigation pulse until the emptying of the irrigation system, and Vavg_field_3phases (l) is the average volume of the 16 MK water volumes measured from the beginning of the irrigation pulse until the emptying of the irrigation system.
  • Distribution uniformity of the total water applied by each emitter:
The distribution uniformity obtained with the empiric model (UDlq_model) was calculated for seven simulations of irrigation pulses. Simulations were made every 5 min, from 5 to 30 min, plus an additional simulation of 60 min. The UDlq_model was calculated using the following equation:
UD lq model = 100 V 25 % model V avg model
where V25%_model (l) is the average of 25% with a lower volume of total water applied (including three phases) by 15,912 drip emitters of the irrigation unit, and Vavg_model (l) is the average of the total water applied by the 15,912 drip emitters of the irrigation unit during the three drip irrigation phases.

2.5.2. Potential Application Efficiency for Zero Deficit

Application efficiency provides a general indication of how well an irrigation system performs in delivering water from the irrigation system to the crop. Therefore, application efficiency is a measure of the ratio between the crop water requirements (ETc) and the total volume of water delivered to the field by the irrigation system (Vwater delivered). ETc was calculated based on FAO methodology, using the measured reference evapotranspiration (ETo) and local crop coefficients (Kc). ETo was obtained from a meteorological station belonging to the Agrometeorological Network of Andalusia [34]. Kc was calculated based on crop coverage using the relationship obtained by [35].
Potential application efficiency without deficit irrigation at any location (PAEzerodeficit) occurs when there is no irrigation deficit, and no over-irrigation, at the point where the irrigation water applied was at a minimum. This irrigation performance indicator is especially important for the irrigation strategies of horticultural crops of high added value, such as berries. For these kinds of crops, the farmer’s production strategy is to maximize yield [14,35], which requires the absence of water deficit throughout the plot. PAEzerodeficit can be defined as:
PAE zerodeficit = 100 neV ETc V waterdelivered
where ne is the number of emitters of the irrigation unit evaluated and VETc is the minimum irrigation volume applied by one emitter in the irrigation unit, assuming that the minimum irrigation volume is equal to the crop water requirements.

2.6. Statistical Methods Used to Compare the Model and Data Fit

In order to carry out the subsequent statistical analysis for validating the irrigation distribution model, the following parameters were calculated.

2.6.1. Coefficient of Determination

The coefficient of determination (R2) is a key output of regression analysis. R2 is interpreted as the proportion of the variance in the dependent variable that is predictable from the independent variable. An R2 between 0 and 1 indicates the extent to which the dependent variable is predictable.

2.6.2. Theil’s Inequality Coefficient

Theil’s inequality coefficient (U) was calculated as [36]:
U = i = 1 n Ob i Pr i 2 i = 1 n Ob i 2
where Obi and Pri are the observed (measured) and predicted values and n is the number of data pairs. Parameter U could be 0 or greater. U = 0 means a perfect fit between the model results and observations. A larger U value means a poorer model performance [36].

2.6.3. Modeling Efficiency

Modeling efficiency (ME) was calculated as [37]:
ME = 1 i = 1 n Ob i Pr i 2 i = 1 n Ob i Pr avg 2
where Pr avg is the average of predicted values. ME = 1 indicates a perfect fit. ME = 0 reveals that the model is no better than a simple average, and negative values indicate a poor performance.

2.6.4. Standardized Residual

Residual methods calculate the difference between the observed and modelled data points. In our case, the residual plot was a plot of the Standardized Residual (SR) as the dependent variable and the location of the measured points as the descriptor variable [38]. The standardized residual was calculated as:
SR i = Ob i Pr i Ob i

3. Results

3.1. Modelling of the Filling Phase

The evolution of pressure measured during the filling phase depended on the location of the emitters. One hundred and forty seconds was the time it took for the entire system to fill and pressurize. The transition from zero pressure to the value of stable pressure is sharper as the distance of the pressure locations increases relative to the water entrance of the irrigation unit (Figure 2a). In this sense, 20 s was the amount of time it took to fully pressurize the system’s most distant extremities. In this regard, the pressure stabilized 160 s after the beginning of irrigation. Therefore, the duration of the filling phase in the irrigation unit in our study was 160 s. In the case of the farthest emitter (L4-G4), there was no pressure in the first 140 s, which implies that the application of different water volumes during the filling phase depended on the location of the emitters in the irrigation unit (Figure 2b).

3.2. Modelling of the Stable Pressure Phase

The decrease in ground level in the direction of the laterals (positive slope, Figure 1) caused a positive gradient of pressure. This positive gradient occurs because the decrease in ground level is higher than the head loss in the laterals. In contrast, in the manifold pipe there was a negative gradient of pressure, which was caused by the increment in ground level (negative slope, Figure 1) and the head loss in the manifold pipe (Table 1).
The increases and decreases in pressure were adjusted to second-degree equations in the laterals (Figure 3a) and in the manifold (Figure 3b). Second-degree equations were used to assign a predicted pressure to each of the 15,912 emitters of the irrigation unit during the stable pressure phase.

3.3. Modelling of the Emptying Phase

The effect of the emptying phase on sloped fields was much higher in the last third of the laterals (Figure 4). The analysis of the water volume measured during the emptying phase of the calibration stage showed that lower elevation laterals are more affected by the emptying phase of the manifold pipe (Table 2). To address the effect of emptying the manifold pipe, the study irrigation unit was divided into two areas: the irrigation area of lower elevation from L1 to L1B (see Figure 1), and the irrigation area of higher elevation from L1B to L4 (see Figure 1). In both areas, the third degree polynomial equations of water volume measured during the emptying phase versus the distance from the lateral beginning were obtained (Figure 4). The two equations shown in Figure 4 allow the simulation of the volume of water applied by the 159,212 emitters of the irrigation unit during the emptying phase.

3.4. Validation of the Water Distribution Model of the Study Irrigation Unit

3.4.1. Filling Phase

The results obtained in the validation phase show that the model is a good predictor of the filling phase. There was an excellent correlation (R2 = 0.9953) between the measured volume (Vfmeasured) and the emptying volumes predicted by the model (Vfmodel) during the emptying phase (Figure 5).

3.4.2. Stable Pressure Phase

The validation data set of the average measured pressure during the stable pressure phase on each of the 16 MK locations shows a clear negative trend in the pressure values of the manifold pipe and a slightly positive trend in the laterals (Table 3).
There was also a very good correlation (R2 = 0.9658) between the measured average pressure (Pspavg measured) and the pressure predicted by the model (Pspmodel) (Figure 6).

3.4.3. Emptying Phase

The validation data set of the average measured volume during the emptying phase on each 16 MK location shows a clear increase in the emptying volume values applied as the distance to the lateral origin increases. This increasing volume trend along the lateral in the validation stage (Table 4) is similar to the one we found with the calibration data set (Table 3).
There was an excellent correlation (R2 = 0.9942) between the measured average volume and the emptying volumes predicted by the irrigation distribution model in both L1-L1B and L1B-L4 (Figure 7).

3.4.4. Model Performance Evaluation

The performance of the irrigation distribution model developed in this study was evaluated by four criteria: the coefficient of determination (R2), the Theil’s index (U), the modelling efficiency (ME), and the standardized residuals (SR). The values of the first three indices are shown in Table 5. The fourth index were plotted against the location of measuring points (Figure 8).
The coefficient of determination (R2) was higher than 0.965 for all the irrigation phases, which implies an excellent goodness of fit. The values of U near to zero for the filling and stable pressure phases indicated an excellent performance in simulating the volume of water applied and the pressure. The emptying phase has a very good performance, but it was worse than the other two irrigation phases. In the case of modelling efficiencies, the ME values were near to 1, which means that the model performed very well in simulating volumes and pressures. Although different indices showed slightly different results, the irrigation distribution model is argued to be reliable in simulating the applied water volume by the emitter in the filling, stable pressure, and emptying phases.
A sensitivity analysis for quantifying the model behavior by the three drip irrigation phases according to the location of the measuring points was carried out. For that purpose, we used a standardized residual plot using the standardized residual (SR, Equation (12)) as the dependent variable and the location of the measuring point as the descriptor variable (Figure 8).
Analyzing the standardized residuals by the drip irrigation phases (Figure 8), the stable pressure phase has the closest values to zero, and therefore the best model fit. The filling phase has similar values to the pressure phase in the first 10 emitter locations (Lateral L1, Figure 1). However, on lateral L3 and especially on lateral L4, the values are worse. In the case of L4, the values of the volumes applied are less than 5 mL, which causes greater experimental errors in volume collection and measurement. The emptying phase shows worse values than the other two phases. The area between L1 and L2, the first third of the study unit (Figure 1), shows the worst values (Figure 8). In this area, there was a transition of equations used in the empirical model to predict the volumes of the emptying phase (Figure 4), which seems to generate greater standardized residuals.

3.5. Influence of the Drip Irrigation Phases on Water Applied According to the Irrigation Distribution Model

We found that, as the duration of the irrigation pulse increased, the influence of the filling and emptying phases decreased (Figure 9). The water volumes applied at each irrigation phase by each of the 15,912 emitters of the irrigation unit, were simulated using the empirical distribution model. The emptying phase had a greater influence than the filling phase on the irrigation water distribution. In any case, the phase of stable pressure had a greater influence than the filling and emptying phases together; however, the volume of water applied during the filling and emptying phases at an irrigation pulse of 5 min could reach half of the total. For irrigation pulses exceeding 20 and 60 min, the filling and emptying phases represent less than 10% and 4% of the total volume of water applied, respectively.

3.6. Sensitive Analyisis of Irrigation Performance Indicators Regarding the Duration of the Irrigation Pulse

UDlq_sp overestimated the distribution uniformity values in comparison with the methodologies that include the three phases of the pipe flow process. According to the duration of the irrigation pulse, the UDlq values change depending on the drip irrigation phases that are included in the calculation of the UDlq. The shorter the irrigation pulse, the greater the difference between the different UDlq was. UDlq _sp not only had the highest value but also was independent of the duration of the irrigation pulse. Both UDlq_field_3phases and UDlq_model had lower values than UDlq_sp. Except for UDlq_sp, the other two UDlq values (UDlq_field_3phases and UDlq_model) varied with the duration of the irrigation pulse. UDlq_field_3phases and UDlq_model had a similar distribution uniformity evolution depending on the duration of the irrigation pulse, even if UDlq_field_3phases had lower values than UDlq_model. The discrepancy between UDlq_field_3phases and UDlq_model was due to the difference in the size of the sample of emitters—16 and 15,912, respectively. Furthermore, the difference between the UDlq values of the field and model, both including three phases, were less than 3% when the duration of irrigation pulse exceeded 10 min. Therefore, the selection of the 16 MK emitters may be sufficient to obtain a good approximation of the water distribution uniformity for practical purposes.
The PAEzerodeficit value increased with the duration of the irrigation pulse (Figure 10). The value of the PAEzerodeficit indicator was 91%, assuming a potential UDlq_sp of 93%. However, in irrigation pulses below 20 min, the PAEzerodeficit was under 85%, while in the 5 min irrigation pulses the PAEzerodeficit was 57%. This deviation from the application efficiency to obtain a zero deficit in the irrigation schedule can have serious repercussions in terms of yield and economic performance.

4. Discussion

We wanted to know how important it was to include the volume applied during the filling and emptying phases of the pipe flow process on the water distribution and application efficiency without water deficit. Typically, distribution uniformity and application efficiency are calculated assuming constant pressures throughout the irrigation pulse [23,24,25,39,40]; however, we hypothesized that it was important to compute the water volume applied during the filling and emptying phases of the pipe flow process for short-duration irrigation pulses, particularly in sloping fields. The study found that the water distribution uniformity and potential application efficiency for zero deficit were greatly affected by the duration of the irrigation pulse for irrigation pulses shorter than 20 min (Figure 10). The distribution of the volume of water applied during the filling and emptying phases of the pipe flow process contributed decisively to this effect (Figure 9).
Moreover, the duration of the irrigation pulse has great relevance in sandy soils. In general, short-duration irrigation pulses are associated with lower water loss [41,42]. Even in a field experiment, high irrigation frequency led to an equivalent or improved yield and increased water productivity [8,43]. Previous studies also indicated the positive role of pulse drip irrigation with increasing the irrigation application efficiency. However, the irrigation events were between 15 and 25 min [44,45], which is consistent with the conclusions of this study.
However, it is important to take into account that the quantitative results obtained in this work are conditioned by the type of tapes and emitters of the studied irrigation unit. The used emitter constitutes a particular case of microirrigation systems. In fact, [46] concluded that the use of pressure-compensating emitters instead of built-in emitters changed the discharge results. In addition, laterals are thin-walled, collapsible polyethylene tapes with built-in emitters. The use of more rigid polyethylene tubes could also change the response of each phase. Despite this constraint, we think the conclusions are sufficiently robust and valid for the purpose of this work. Based on these results, new steps must be taken to generalize the results in any irrigation unit with any kind of irrigation system. Including a rational model in a decision-support tool will help to manage short irrigation pulses under more arbitrary irrigation material, soil, and slope conditions.
The distribution uniformity and potential application efficiency for zero deficit were overestimated for the short irrigation duration if some causes of non-uniformity are ignored. The greater the percentage of volume of water applied during the filling and emptying phases, the worse the performance of the water distribution uniformity and potential application efficiency for zero deficit are (Figure 9 and Figure 10). The distribution uniformity value was 93%, which is classified as “excellent” according to Merrian and Keller’s classification, using only stable pressure phase in a 5 min irrigation pulse. However, the value estimated by the water distribution model was 65% when the filling and emptying phases are included, which is classified as “poor” (Figure 10). Therefore, including the filling and emptying phases has profound practical implications for optimal irrigation scheduling in sloping sandy soils with very short irrigation durations. Optimal irrigation schedules are especially important in sandy soil because sandy soils require a balance in water application efficiency and deep percolation to save water without reducing yields.

5. Conclusions

The filling and emptying phases of drip irrigation both have an impact on accurately calculating water distribution uniformity and potential application efficiency with zero deficit. The influence of the filling and emptying phases varies according to the duration of the irrigation pulse. The shorter the duration of the irrigation pulse, the greater the influence of the filling and emptying phases on the irrigation performance indicators. For the better design and management of the drip irrigation system, it is necessary to know the upper threshold values of the real distribution uniformity and potential application efficiency for zero deficit depending on the duration of the irrigation pulse.
The irrigation performance indicators of very high frequency drip irrigation in sloping fields must include the total volume applied by the drip emitters. In order to include the total volume applied by the drip emitters, it is necessary to characterize the filling and emptying phases of the irrigation unit. This characterization can contribute significantly to improve irrigation scheduling, especially in sloping fields, which constitutes a significant breakthrough. The next step should be to make a rational model that allows a generalized quantitative analysis of any drip irrigation unit. Even so, we believe that the identification of this issue in drip irrigation systems in order to calculate irrigation performance indicators more accurately is an original contribution.

Author Contributions

Software and writing—original draft preparation, D.L.; Conceptualization, methodology, investigation, formal analysis, resources, funding acquisition, and writing—review and editing, D.L. and P.G.; data acquisition, validation, and supervision, D.L., N.R., J.I.C., R.B., and P.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially funded by the contract between IFAPA and “SAT Las Palmeritas” for the execution of the technical assistance “Experimentación y transferencia tecnológica para la mejora del riego de la fresa en la SAT ‘Las Palmeritas’ en el t.m. de Lepe (Huelva)” (CAICEM 134/2015) and for the project of “Gestión sostenible del regadío en la agricultura intensiva de Andalucía” (RTA2015-00029-C02). Financed by INIA and FEDER 2014-2020 (“Programa Operativo de Crecimiento Inteligente”).

Acknowledgments

The authors would like to thank the invaluable collaboration of the company “SAT Las Palmeritas” and especially the field staff who have collaborated intensively in this work. We thank the anonymous reviewers who provided constructive comments that improved this manuscript. David Lozano is grateful for the financing of his contract to the “Subprograma de Incorporación del Programa Estatal de Promoción del Talento y su Empleabilidad del Plan Estatal de Investigación Científica y Técnica y de Innovación 2013-2016 (DOC-INIA)”, co-financed by INIA, and for the European Social Fund (FSE).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

4 MK4 emitters with a lower water volume measured during the three phases (filling, stable pressures and emptying) using the Merriam and Keller’s methodology
16 MK16 emitters selected by the Merriam and Keller’s methodology to calculate water distribution uniformity in a drip irrigation system
ΔTfTime interval of pressure recorded (10 s)
ΔPLPressure variation depending on lateral distance (bar)
ΔPmPressure variation depending on manifold distance (bar)
DLDistance from lateral beginning (m)
DmDistance from manifold beginning (m)
e: emitter. Range from 1 to 16, corresponding to 16 MK (Figure 1), during the filling and the pressure stable phases. Range from 1 to 60 during the emptying phase
EToReference evapotranspiration
ETcCrop water requirements
G1Measuring points at the beginning of the laterals
G2Measuring points at a third of the distance to the end of the laterals
G3Measuring points at two third of the distance to the end of the laterals
G4Measuring points at the end of the laterals
iIndex of data pairs of observed and predicted values. 1 to 16 for emptying and stable pressure phases. 1 to 60 for irrigation pulse stage and emptying phase
jTime interval for measuring pressure. Range from 1 to 30. One value each 10 s during 5 min
KcCrop coefficient
L1Irrigation lateral located at the beginning of the manifold pipe
L1AIrrigation lateral at a quarter of the L1 to L2 distance
L1BIrrigation lateral at a half of the L1 to L2 distance
L2Irrigation lateral located at a third of the distance to the end of the manifold pipe
L3Irrigation lateral located at two third of the distance to the end of the manifold pipe
L4Irrigation lateral located at the end of the manifold pipe
MEModelling efficiency
nNumber of measuring points. 16 for filling and stable pressure phases and 60 for emptying phase
neTotal number of emitters of the irrigation unit evaluated
ObiObserved (measured) values
PfejPressure measured in the manometer intake e at the instant j during the filling phase (meter water column)
Pr avg Average of predicted values
PriPredicted values by a model
Pspavg measuredAverage stable phase pressure measured from the validation data set (bar)
PspmodelStable phase pressure predicted by the irrigation distribution model (bar)
PspavgAverage pressure measured in the manometer intakes e during the stable pressure phase (metre water column)
PAEzerodeficitPotential application efficiency without deficit irrigation at any location
Q25%25% average flow with a lower water flow (l h−1)
Q25%-sp25% average flow of 16 MK emitters with a lower water flow during the stable pressure phase (l h−1)
QavgAverage flow (l h−1)
Qavg-spAverage flow of 16 MK emitters during the stable pressure phase (l h−1)
QfejFlow rate of the emitter e at the instant j during the filling phase (l h−1)
QspeFlow rate of the emitter e during the stable pressure phase (l h−1)
SRStandardized residual
SDStandard deviation
TfNumber of the first time interval measured once the pressure is stable
TspDuration of the stable pressure phase (h)
UTheil’s inequality coefficient
UDlqDistribution uniformity of the low quarter (%)
UDlq_field_3phasesDistribution uniformity of the low quarter during filling, stable pressures and emptying phases measuring in 16 MK emitters
UDlq_modelDistribution uniformity of the low quarter during the three irrigation phases (filling, pressure stables and emptying) calculated with the irrigation distribution model for the 15,912 drip emitters of the study unit
UDlq_spDistribution uniformity of the low quarter on stable pressure phase measuring in the 16 MK emitters (%)
V25%25% average water applied volume with a lower volume of water (l)
V25%_model25% average of a lower volume of total water applied by the 15,9212 drip emitters of the study irrigation unit calculated with the empirical model
V25%_field_3phasesAverage of 4 MK emitters with a lower water volume measured during the three irrigation phases (filling, stable pressures and emptying)
VavgAverage water applied volume (l)
Vavg_field_3phasesAverage measured volume of 16 MK emitters during the three irrigation phases (filling, stable pressures and emptying)
Vavg_modelAverage of total water applied by the 15,912 drip emitters of the study irrigation unit calculated with the empirical model
Vepavg measured Average water volume measured during the emptying phase from validation data set (ml)
VepmodelWater volume predicted by the irrigation distribution model during the emptying phase (ml)
VETcVolume of crop water requirements (l)
VfeWater volume applied by the emitter e during the filling phase (l)
VfmodelPredicted water volume applied by emitters during the filling phase using the irrigation distribution model (ml)
VfmeasuredMeasured water volume applied by emitters during the filling phase (ml)
VG1Water volume applied by emitters in position G1 (Figure 1) during the filling phase (ml)
VG2G3G4Water volume applied by emitters in position G2, G3 and G4 (Figure 1) during the filling phase (ml)
VL1-L1BWater volume applied by emitters between the laterals L1 and L1B during the emptying phase (m)
VL1B-L4Water volume applied by emitters between the laterals L1B and L4 during the emptying phase (m)
VspeWater volume applied by the emitter e during the stable pressure phase (l)
Vwater deliveredTotal volume of water delivered to the field by the irrigation system

References

  1. Hsiao, T.C.; Steduto, P.; Fereres, E. A quantitative framework for the systematic analysis of potential water savings in agriculture. In Book Water Saving in Mediterranean Agriculture and Future Research Needs; Options Méditerranéennes: Série B; Lamaddalena, N., Bogliotti, C., Todorovic, M., Scardigno, A., Eds.; Etudes et Recherches 56: Valenzano, Italy, 2007; Volume 1, pp. 37–47. [Google Scholar]
  2. Lorite, I.J.; Ruiz-Ramos, M.; Gabaldón-Leal, C.; Cruz-Blanco, M.; Porras, R.; Santos, C. Water management and climate change in semiarid environments. In Water Scarcity and Sustainable Agriculture in Semiarid Environment, 1st ed.; García-Tejero, I., Duran, V.H., Eds.; Academic Press: Cambridge, MA, USA, 2018; pp. 3–40. [Google Scholar]
  3. Oñate, J.J.; Pereira, D.; Suarez, F. Strategic environmental assessment of the effects of European Union’s regional development plans in Donana National Park (Spain). Environ. Manag. 2003, 31, 0642–0655. [Google Scholar] [CrossRef] [PubMed]
  4. BOE. Apéndice 9. Dotaciones y Eficiencias. In Plan Hidrológico de la Demarcación Hidrográfica del Guadalquivir (2015–2021); Agencia Estatal Boletín Oficial del Estado: Madrid, Spain, 2016; Volume 16, pp. 3686–3809. [Google Scholar]
  5. García Morillo, J.; Rodríguez Díaz, J.A.; Camacho, E.; Montesinos, P. Drip Irrigation Scheduling Using Hydrus 2-D Numerical Model Application for Strawberry Production in South-West Spain. Irrig. Drain. 2017, 66, 797–807. [Google Scholar] [CrossRef]
  6. Clark, G.A.; Smajstrla, A.G. Water distributions in soils as influenced by irrigation depths and intensities. Proc. Soil Crop Sci. Soc. Florida 1983, 42, 157–165. [Google Scholar]
  7. Clark, G.A.; Stanley, C.D.; Zazueta, F.S. Qualitative sensing of water movement from a point-source emitter on a sandy soil. Appl. Eng. Agric. 1993, 9, 299–303. [Google Scholar] [CrossRef]
  8. Dukes, M.D.; Simonne, E.H.; Davis, W.E.; Studstill, D.W.; Hochmuth, R. Effect of sensor-based high frequency irrigation on bell pepper yield and water use. In Proceedings of the 2nd International Conference on Irrigation and Drainage, Phoenix, AZ, USA, 12–15 May 2003; pp. 12–15. [Google Scholar]
  9. Skaggs, T.H.; Trout, T.J.; Rothfuss, Y. Drip irrigation water distribution patterns: Effects of emitter rate, pulsing, and antecedent water. Soil Sci. Soc. Am. J. 2010, 74, 1886–1896. [Google Scholar] [CrossRef] [Green Version]
  10. Cote, C.M.; Bristow, K.L.; Charlesworth, P.B.; Cook, F.J.; Thorburn, P.J. Analysis of soil wetting and solute transport in subsurface trickle irrigation. Irrig. Sci. 2003, 22, 143–156. [Google Scholar] [CrossRef]
  11. Lembke, W.D.; Thorne, M.D. Nitrate leaching and irrigated corn production with organic and inorganic fertilizers on sandy soil. Trans. ASAE 1980, 23, 1153–1156. [Google Scholar] [CrossRef]
  12. Watts, D.G.; Martin, D.L. Effects of water and nitrogen management on nitrate leaching loss from sands. Trans. ASAE 1981, 24, 911–916. [Google Scholar] [CrossRef]
  13. García-Morillo, J.; Martín, M.; Camacho, E.; Rodríguez-Díaz, J.A.; Montesinos, P. Toward precision irrigation for intensive strawberry cultivation. Agric. Water Manag. 2015, 151, 43–51. [Google Scholar] [CrossRef]
  14. Gendron, L.; Létourneau, G.; Cormier, J.; Depardieu, C.; Boily, C.; Levallois, R.; Caron, J. Using Pulsed Water Applications and Automation Technology to Improve Irrigation Practices in Strawberry Production. HortTechnology 2018, 28, 642–650. [Google Scholar] [CrossRef]
  15. Létourneau, G.; Caron, J. Irrigation Management Scale and Water Application Method to Improve Yield and Water Productivity of Field-Grown Strawberries. Agronomy 2019, 9, 286. [Google Scholar] [CrossRef] [Green Version]
  16. Kang, Y.; Yuan, B.Z.; Nishiyama, S. Design of microirrigation laterals at minimum cost. Irrig. Sci. 1999, 18, 125–133. [Google Scholar] [CrossRef]
  17. Juana, L.; Rodríguez-Sinobas, L.; Sánchez, R.; Losada, A. Analytical expressions for hydraulic calculation of trapezoidal drip irrigation units. J. Irrig. Drain. Eng. 2005, 131, 420–432. [Google Scholar] [CrossRef]
  18. Ella, V.B.; Reyes, M.R.; Yoder, R. Effect of hydraulic head and slope on water distribution uniformity of a low-cost drip irrigation system. Appl. Eng. Agric. 2009, 25, 349–356. [Google Scholar] [CrossRef] [Green Version]
  19. Jiang, S.; Kang, Y. Evaluation of microirrigation uniformity on laterals considering field slopes. J. Irrig. Drain. Eng. 2010, 136, 429–434. [Google Scholar] [CrossRef]
  20. Zayani, K.; Hammami, M.; Alouini, A.; Souissi, A. Design of non-zero uniformly sloping laterals in trickle irrigation systems. J. Irrig. Drain. Eng. 2013, 139, 419–425. [Google Scholar]
  21. Baeza, R.; Gavilán, P.; Vargas, A.; Contreras, J.I. Influencia de la pendiente del terreno en la uniformidad de distribución de caudal en cintas de riego localizado. In Proceedings of the XXXII Congreso Nacional de Riegos, Madrid, Spain, 10–12 June 2014. [Google Scholar]
  22. Zapata, A.J.; López-Segura, J.G.; Alonso-López, F.; Cánovas-Fernández, G.; Baeza Cano, R.; Lozano Pérez, D.; Contreras Paris, J.I. Uniformidad de distribución en ramales de riego instalados en pendiente: Influencia del periodo de transición de vaciado. In Proceedings of the XXXV Congreso Nacional de Riegos, Tarragona, Spain, 6–8 June 2017. [Google Scholar]
  23. Merriam, J.L.; Keller, J. Farm Irrigation System Evaluation: A Guide to Management; Utah State University: Logan, UT, USA, 1978. [Google Scholar]
  24. ASAE. Design and installation of microirrigation systems. ASAE, EP405.1. In ASAE Standards 1998; American Society of Agricultural Engineers: St. Joseph, MI, USA, 1998. [Google Scholar]
  25. Burt, C.M. Rapid field evaluation of drip and microspray distribution uniformity. Irrig. Drain. Syst. 2004, 18, 275–297. [Google Scholar] [CrossRef]
  26. Bralts, V.F.; Kesner, D. Drip irrigation field uniformity estimation. Trans. ASAE 1983, 26, 1367–1374. [Google Scholar]
  27. Bralts, V.F.; Edwards, D.M. Field evaluation of drip irrigation submain units. Trans. ASAE 1986, 29, 1659–1664. [Google Scholar] [CrossRef]
  28. Bralts, V.F.; Edwards, D.M.; Wu, I.P. Drip irrigation design and evaluation based on the statistical uniformity concept. In Book Advances in Irrigation; Hillel, D., Ed.; Academic Press: New York, NY, USA, 1987; Volume 4, pp. 67–117. [Google Scholar]
  29. Kruse, E.G. Describing irrigation efficiency and uniformity. J. Irrig. Drain. Eng. Div. 1978, 104, 35–41. [Google Scholar]
  30. Burt, C.M.; Clemmens, A.J.; Strelkoff, T.S.; Solomon, K.H.; Bliesner, R.D.; Hardy, L.A.; Eisenhauer, D.E. Irrigation performance measures: Efficiency and uniformity. J. Irrig. Drain. Eng. 1997, 123, 423–442. [Google Scholar] [CrossRef] [Green Version]
  31. Clemmens, A.J.; Solomon, K.H. Estimation of global irrigation distribution uniformity. J. Irrig. Drain. Eng. 1997, 123, 454–461. [Google Scholar] [CrossRef]
  32. Wu, I.P.; Barragán, J. Design criteria for microirrigation systems. Trans. ASAE 2000, 43, 1145–1154. [Google Scholar]
  33. Barragan, J.; Bralts, V.; Wu, I.P. Assessment of emission uniformity for micro-irrigation design. Biosyst. Eng. 2006, 93, 89–97. [Google Scholar] [CrossRef]
  34. Gavilán, P.; Estévez, J.; Berengena, J. Comparison of standardized reference evapotranspiration equations in Southern Spain. J. Irrig. Drain. Eng. 2008, 134(1), 1–12. [Google Scholar] [CrossRef]
  35. Lozano, D.; Ruiz, N.; Gavilán, P. Consumptive water use and irrigation performance of strawberries. Agric. Water Manag. 2016, 169, 44–51. [Google Scholar] [CrossRef]
  36. Tian, H.; Xu, X.; Miao, S.; Sindhoj, E.; Beltran, B.J.; Pan, Z. Modelling ecosystem responses to prescribed fires in a phosphorus-enriched Everglades wetland: I. Phosphorus dynamics and cattail recovery. Ecol. Model. 2010, 221, 1252–1266. [Google Scholar] [CrossRef]
  37. Vanclay, J.K.; Skovsgaard, J.P. Evaluating forest growth models. Ecol. Model. 1997, 98, 1–12. [Google Scholar] [CrossRef] [Green Version]
  38. Bennett, N.D.; Croke, B.F.; Guariso, G.; Guillaume, J.H.; Hamilton, S.H.; Jakeman, A.J.; Marsili-Libelli, S.; Newham, L.T.H.; Norton, J.P.; Perrin, C.; et al. Characterising performance of environmental models. Environ. Modell. Softw. 2013, 40, 1–20. [Google Scholar] [CrossRef]
  39. Juana, L.; Rodríguez-Sinobas, L.; Sánchez, R.; Losada, A. Evaluation of drip irrigation: Selection of emitters and hydraulic characterization of trapezoidal units. Agric. Water Manag. 2007, 90, 13–26. [Google Scholar] [CrossRef]
  40. Xu, X.; Tian, H.; Pan, Z.; Thomas, C.R. Modelling ecosystem responses to prescribed fires in a phosphorus-enriched Everglades wetland: II. Phosphorus dynamics and community shift in response to hydrological and seasonal scenarios. Ecol. Model. 2011, 222, 3942–3956. [Google Scholar] [CrossRef]
  41. Coolong, T.; Surendran, S.; Warner, R. Evaluation of irrigation threshold and duration for tomato grown in a silt loam soil. HortTechnology 2011, 21, 466–473. [Google Scholar] [CrossRef] [Green Version]
  42. Eid, A.R.; Bakry, B.A.; Taha, M.H. Effect of pulse drip irrigation and mulching systems on yield, quality traits and irrigation water use efficiency of soybean under sandy soil conditions. Agric. Sci. 2013, 5, 249–261. [Google Scholar] [CrossRef] [Green Version]
  43. Segal, E.; Ben-Gal, A.; Shani, U. Root Water Uptake Efficiency under Ultra-High Irrigation Frequency. Plant Soil 2006, 282, 333–341. [Google Scholar] [CrossRef]
  44. Vázquez, N.; Pardo, A.; Suso, M.L.; Quemada, M. Drainage and nitrate leaching under processing tomato growth with drip irrigation and plastic mulching. Agric. Ecosyst. Environ. 2006, 112, 313–323. [Google Scholar] [CrossRef]
  45. Zotarelli, L.; Scholberg, J.M.; Dukes, M.D.; Muñoz-Carpena, R.; Icerman, J. Tomato yield, biomass accumulation, root distribution and irrigation water use efficiency on a sandy soil, as affected by nitrogen rate and irrigation scheduling. Agric. Water Manag. 2009, 96, 23–34. [Google Scholar] [CrossRef]
  46. Contreras, J.I.; Alonso, F.; Cánovas, G.; Baeza, R. Determinación de la ecuación que define la curva de descarga de una cinta de riego en función del tipo de emisor y la pendiente del terreno. Acta nº 71. In Proceedings of the XIV Congreso Nacional de Ciencias Hortícolas, Orihuela, Alicante, Spain, 3–5 June 2015. [Google Scholar]
Figure 1. Scheme of the irrigation unit evaluated. Measurements were made at 60 points, including locations that correspond to Merriam and Keller’s measuring points (16 Merriam and Keller (MK) locations).
Figure 1. Scheme of the irrigation unit evaluated. Measurements were made at 60 points, including locations that correspond to Merriam and Keller’s measuring points (16 Merriam and Keller (MK) locations).
Water 12 02276 g001
Figure 2. Modelling of the filling phase (160 s) in the 16 MK locations (L1-G1,…, L4-G4) using (a) the evolution of pressure from the beginning of irrigation or (b) the volume of water applied by emitters as a function of the distance from the beginning of the manifold (Dm).
Figure 2. Modelling of the filling phase (160 s) in the 16 MK locations (L1-G1,…, L4-G4) using (a) the evolution of pressure from the beginning of irrigation or (b) the volume of water applied by emitters as a function of the distance from the beginning of the manifold (Dm).
Water 12 02276 g002
Figure 3. Pressure variation (ΔP) with respect to the point closest to the water inlet of the irrigation unit (L1-G1) during the stable pressure phase (a) from the lateral beginning (DL) and (b) from the manifold beginning (Dm).
Figure 3. Pressure variation (ΔP) with respect to the point closest to the water inlet of the irrigation unit (L1-G1) during the stable pressure phase (a) from the lateral beginning (DL) and (b) from the manifold beginning (Dm).
Water 12 02276 g003
Figure 4. Average and standard deviation of the volume of water applied by the emitters in the emptying phases in four irrigation pulses (5, 10, 15, and 20 min) from the lateral beginning (DL). The two equations were adjusted for the areas between the laterals L1 and L1-B (L1-L1B) and L1-B and L4 (L1B-L4) using the calibration data set.
Figure 4. Average and standard deviation of the volume of water applied by the emitters in the emptying phases in four irrigation pulses (5, 10, 15, and 20 min) from the lateral beginning (DL). The two equations were adjusted for the areas between the laterals L1 and L1-B (L1-L1B) and L1-B and L4 (L1B-L4) using the calibration data set.
Water 12 02276 g004
Figure 5. Evolution of the measured volume and the volume estimated by the irrigation distribution model during the filling phase (160 s) in 16 MK locations during 5 min irrigation pulses in May.
Figure 5. Evolution of the measured volume and the volume estimated by the irrigation distribution model during the filling phase (160 s) in 16 MK locations during 5 min irrigation pulses in May.
Water 12 02276 g005
Figure 6. Correlation between the stable pressure phase from the validation data set (Pspavg measured) and the stable pressure phase predicted by the irrigation distribution model (Pspmodel) using the calibration data set.
Figure 6. Correlation between the stable pressure phase from the validation data set (Pspavg measured) and the stable pressure phase predicted by the irrigation distribution model (Pspmodel) using the calibration data set.
Water 12 02276 g006
Figure 7. Correlation between the volumes measured during the emptying phase from the validation data set (Vepavg measured) and the volume predicted by the irrigation distribution model during the emptying phase (Vepmodel).
Figure 7. Correlation between the volumes measured during the emptying phase from the validation data set (Vepavg measured) and the volume predicted by the irrigation distribution model during the emptying phase (Vepmodel).
Water 12 02276 g007
Figure 8. Standardized residuals plot regarding the location of the measuring points for the filling, stable pressure, and emptying phases.
Figure 8. Standardized residuals plot regarding the location of the measuring points for the filling, stable pressure, and emptying phases.
Water 12 02276 g008
Figure 9. Percentage of water applied per irrigation phase using the irrigation distribution model.
Figure 9. Percentage of water applied per irrigation phase using the irrigation distribution model.
Water 12 02276 g009
Figure 10. Distribution uniformity (%) according to the duration of the irrigation pulse taking into account the average pressure at 16 MK from the validation data set (UDlq_sp), the total water applied by the 16 MK from the beginning of the irrigation pulse until the end of the emptying phase from the validation data set (UDlq_field_3phases), and the volume of water applied by the 15,912 emitters of the irrigation unit using the irrigation distribution model (UDlq_model). Potential application efficiency for zero deficit (PAEzerodeficit) using the irrigation distribution model.
Figure 10. Distribution uniformity (%) according to the duration of the irrigation pulse taking into account the average pressure at 16 MK from the validation data set (UDlq_sp), the total water applied by the 16 MK from the beginning of the irrigation pulse until the end of the emptying phase from the validation data set (UDlq_field_3phases), and the volume of water applied by the 15,912 emitters of the irrigation unit using the irrigation distribution model (UDlq_model). Potential application efficiency for zero deficit (PAEzerodeficit) using the irrigation distribution model.
Water 12 02276 g010
Table 1. Calibration data set of average pressure (bar) measured during the stable pressure phase on each 16 MK location for 4 different irrigation pulses (5, 10, 15, and 20 min) in November.
Table 1. Calibration data set of average pressure (bar) measured during the stable pressure phase on each 16 MK location for 4 different irrigation pulses (5, 10, 15, and 20 min) in November.
Irrigation Pulse DurationL1-G1L2-G1L3-G1L4-G1L1-G2L2-G2L3-G2L4-G2
(bar)(bar)(bar)(bar)(bar)(bar)(bar)(bar)
5 min0.990.880.800.721.020.900.800.70
10 min0.930.810.720.660.950.840.700.65
15 min0.930.770.670.670.960.770.700.69
20 min0.900.700.650.660.850.720.670.67
Average0.940.790.710.680.950.810.720.68
SD0.030.070.060.020.060.070.050.02
L1-G3L2-G3L3-G3L4-G3L1-G4L2-G4L3-G4L4-G4
(bar)(bar)(bar)(bar)(bar)(bar)(bar)(bar)
5 min1.000.880.820.681.030.820.820.80
10 min1.020.850.700.71.070.870.730.75
15 min0.990.830.740.731.070.880.760.77
20 min0.910.760.70.710.950.820.750.77
Average0.980.830.740.711.030.850.770.77
SD0.040.040.050.020.050.030.030.02
Table 2. Calibration data set of the applied water volume (mL) during the emptying phase of the irrigation pulse in 60 selected locations for 4 different irrigation pulses (5, 10, 15, and 20 min) in November.
Table 2. Calibration data set of the applied water volume (mL) during the emptying phase of the irrigation pulse in 60 selected locations for 4 different irrigation pulses (5, 10, 15, and 20 min) in November.
Dm (m)L1L1AL1BL2L3L4L1L1BL1BL4L1L1AL1BL2L3L4L1L1BL1BL4
Average Volume (m)Standard Deviation (ml)
0.0000000000.000.000.000.000.000.000.000.00
7.3955221725.594.426.022.922.050.505.002.87
14.76011454361.500.008.010.000.710.500.752.30
22.01651191141091.640.872.872.595.723.271.253.61
29.329252614192027200.876.125.362.380.871.093.492.42
32.341354027273438322.493.273.966.304.664.442.884.84
35.253515136443852422.282.862.356.702.123.392.573.64
38.178676346485372531.805.022.458.733.277.353.415.45
41.192868061667389701.797.125.437.957.016.164.456.64
44.0126111111697093119863.4611.8019.735.895.395.107.639.03
Table 3. Validation data set of the average pressure (bar) measured during the stable pressure phase on each 16 MK location for 4 different irrigation pulses (5, 10, 15, and 20 min) in May.
Table 3. Validation data set of the average pressure (bar) measured during the stable pressure phase on each 16 MK location for 4 different irrigation pulses (5, 10, 15, and 20 min) in May.
Irrigation Pulse DurationL1-G1L2-G1L3-G1L4-G1L1-G2L2-G2L3-G2L4-G2
(bar)(bar)(bar)(bar)(bar)(bar)(bar)(bar)
5 min1.000.890.800.711.000.910.800.72
10 min0.980.830.750.711.000.850.750.72
15 min0.990.830.700.651.000.840.700.68
20 min1.010.810.700.681.020.810.700.68
Average1.000.840.740.691.010.850.740.70
SD0.010.030.040.020.010.040.040.02
L1-G3L2-G3L3-G3L4-G3L1-G4L2-G4L3-G4L4-G4
(bar)(bar)(bar)(bar)(bar)(bar)(bar)(bar)
5 min1.020.900.800.701.050.880.850.80
10 min1.020.880.800.701.090.880.850.89
15 min1.000.830.770.701.030.860.800.74
20 min1.040.860.740.701.080.850.800.72
Average1.020.870.780.701.060.870.830.79
SD0.010.030.020.000.020.010.030.07
Table 4. Validation data set of the applied water volume (ml) during the emptying phase of the irrigation pulse in 60 selected locations for 4 different irrigation pulses (5, 10, 15, and 20 min) in May.
Table 4. Validation data set of the applied water volume (ml) during the emptying phase of the irrigation pulse in 60 selected locations for 4 different irrigation pulses (5, 10, 15, and 20 min) in May.
L1L1AL1BL2L3L4L1L1BL1BL4L1L1AL1BL2L3L4L1L1BL1BL4
Dm (m)Average Volume (ml)Standard Deviation (ml)
0.0000000000.000.000.000.000.000.000.000.00
7.31166300924.152.594.330.000.004.154.531.73
14.74011120230.002.871.000.870.430.000.651.29
22.01912189104151011.263.037.533.563.7711.266.284.47
29.326282616201927203.462.871.920.501.923.462.671.80
32.343344332282838334.393.467.362.122.054.393.223.75
35.250505123473850400.872.177.713.422.120.870.983.85
38.193716545535882552.873.278.417.3813.992.875.558.26
41.110710078607370103706.385.726.388.0714.676.385.248.71
44.01231299678871071269214.7214.977.508.954.9514.7212.409.09
Table 5. Statistical comparison of the overall fit of the irrigation distribution model outputs against the measured data.
Table 5. Statistical comparison of the overall fit of the irrigation distribution model outputs against the measured data.
Irrigation Phases/IndicesR2UME
Filling phase0.99530.0570.99
Stable pressure phase0.96580.0400.92
Emptying phase0.99420.1410.96
Note: R2, coefficient of determination; U, Theil’s coefficient; ME, modelling efficiency.

Share and Cite

MDPI and ACS Style

Lozano, D.; Ruiz, N.; Baeza, R.; Contreras, J.I.; Gavilán, P. Effect of Pulse Drip Irrigation Duration on Water Distribution Uniformity. Water 2020, 12, 2276. https://doi.org/10.3390/w12082276

AMA Style

Lozano D, Ruiz N, Baeza R, Contreras JI, Gavilán P. Effect of Pulse Drip Irrigation Duration on Water Distribution Uniformity. Water. 2020; 12(8):2276. https://doi.org/10.3390/w12082276

Chicago/Turabian Style

Lozano, David, Natividad Ruiz, Rafael Baeza, Juana I. Contreras, and Pedro Gavilán. 2020. "Effect of Pulse Drip Irrigation Duration on Water Distribution Uniformity" Water 12, no. 8: 2276. https://doi.org/10.3390/w12082276

APA Style

Lozano, D., Ruiz, N., Baeza, R., Contreras, J. I., & Gavilán, P. (2020). Effect of Pulse Drip Irrigation Duration on Water Distribution Uniformity. Water, 12(8), 2276. https://doi.org/10.3390/w12082276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop