Next Article in Journal
A Study of the Performance of Dielectric Barrier Discharge under Different Conditions for Nitrobenzene Degradation
Next Article in Special Issue
Simulation of Long-Term Soil Hydrological Conditions at Three Agricultural Experimental Field Plots Compared with Measurements
Previous Article in Journal
Assessing the Effectiveness of the WFD as a Tool to Address Different Levels of Water Scarcity Based on Two Case Studies of the Mediterranean Region
Previous Article in Special Issue
Estimates of Tillage and Rainfall Effects on Unsaturated Hydraulic Conductivity in a Small Central European Agricultural Catchment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying Optimal Irrigation Water Needs at District Scale by Using a Physically Based Agro-Hydrological Model

1
School of Agricultural, Forestry, Food and Environmental Sciences (SAFE), University of Basilicata, Viale dell’Ateneo Lucano, 10, 85100 Potenza PZ, Italy
2
Mediterranean Agronomic Institute, Land and Water Division, IAMB, 70010 Valenzano BA, Italy
3
Institute of Mediterranean Agriculture and Forestry Systems (ISAFOM), National Council of Research (CNR), 80056 Naples NA, Italy
4
Consorzio di Bonifica della Capitanata, 71121 Foggia FG, Italy
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 841; https://doi.org/10.3390/w11040841
Submission received: 10 March 2019 / Revised: 13 April 2019 / Accepted: 15 April 2019 / Published: 21 April 2019
(This article belongs to the Special Issue Soil Hydrology in Agriculture)

Abstract

:
This paper mainly aims to illustrate an irrigation management tool to simulate scheduling of district-level water needs over the course of an irrigation season. The tool is mostly based on a daily model for simulating flow of water (and solutes) in heterogeneous agri-environmental systems (called FLOWS-HAGES). The model produces information on the daily evolution of: soil water contents and pressure potentials in the soil profile; water uptake and actual evapotranspiration; stress periods for each crop; return fluxes to the groundwater and their quality in terms of solute concentrations (e.g., nitrates). FLOWS-HAGES provides a daily list of hydrants to be operated according to water or crop-based criteria. The daily optimal sequence of hydrant use may thus be established by passing the volumes to be delivered on to the model for simulating the hydraulics of the irrigation network, in order to ensure that the discharges flowing inside the network of distribution pipes are delivered under optimal pressure head distribution in the system. All the above evaluations can be carried out in a stochastic framework to account for soil heterogeneity and climate changes. To illustrate the potential of FLOWS-HAGES, a case study was considered for a selected sector of the Irrigation District 10 in the “Sinistra Ofanto” irrigation system (southern Italy, Apulia region). In a 139 ha area (Sector 6 of the Irrigation District), soil profiles were analyzed for characterization of hydraulic properties variability. Hydraulic properties were determined by a combination of field and laboratory measurements. Model simulations were validated by comparing soil water storage simulated and measured by a sensor based on electromagnetic induction technique. Irrigation water volumes and frequency calculated by the model were compared to the volumes actually supplied by the farmers. Compared to the farmers behavior, the model simulates more frequent irrigations with lower irrigation volumes. Finally, some indexes of irrigation performance were calculated for each farm under study. The resulting maps provide useful information on the spatial distribution of farmer behavior, indicating the abuse or underuse of water as well as the fraction of the water lost by drainage following the irrigation method applied.

1. Introduction

Water scarcity and quality degradation are set to become the main environmental problems for all countries in the Mediterranean region in the near future. Due mainly to population growth, pressure on agriculture and demand for water to irrigate food crops have also increased. Indeed, agriculture in Mediterranean countries (apart from those in the eastern Adriatic and France) accounted for 64% of overall demand in the period 2005–2010 (49% in the North, 74% and 81% in the South and East) [1]. In Mediterranean countries the irrigated area has more than doubled in 40 years, totaling 24,200,000 ha in 2009 (17.8 million in the Mediterranean Europe and 6.4 million in Northern Africa). The growing demand for irrigation water will be mainly met by over-exploiting both renewable and fossil groundwater, prompting seawater intrusion which leads to soil salinization problems in the frequent cases of ineffective irrigation management. Indeed, irrigation is often wasteful and highly inefficient [2]. In an attempt to tackle the problem in Italy and other Mediterranean regions, open-channel irrigation systems have been (or are going to be) converted into pressurized pipeline networks. In most cases, the "on-demand" method for water distribution has been adopted to replace the old rotational schedule.
However, experience shows that even for on-demand pressurized irrigation systems the performance of distribution networks is still frequently far from acceptable in terms of efficiency for both (1) economic/political and (2) technical reasons.
First, even in on-demand irrigation distribution networks, the problem of unmeasured and generally excessive application of irrigation water frequently occurs. This is mostly due to the fact that the price of the water does not reflect the true cost of provision but is kept deliberately low, sometimes for political reasons, by government subsidy. In such cases, the price of irrigation water may be lower than the cost of the labor or equipment needed to avoid over-irrigation, both for the farmer and the distribution network manager. Where water is delivered to the consumer only at fixed times and/or charges are set either per delivery or per hectare regardless of the actual water amount used, farmers tend to take as much water as they can while they can. This generally results in over-irrigation, which wastes both water and nutrients. In some Mediterranean countries the excessive water applied may lead to other serious problems, such as rising water tables, waterlogging and upward movement of salts.
Secondly, in the on-demand framework the manager of the distribution network has little control over water scheduling and hence on the performance of the distribution network. According to the particular combination of farmers who decide to irrigate at a given time (which is not known a priori), the discharges flowing inside the distribution network may change. Under such conditions the actual availability of water resources and the ability of the conveyance hydraulic system may not be able to provide the discharge which the farmers require at a given time. Moreover, under changing (and unknown a priori) discharges, the pressure head in the system can be subject to high fluctuation, with consequent impacts on the efficiency of on-farm irrigation systems [3] also in terms of energy to be spent on water distribution.
The twofold—economic and technical—nature of the problem calls for different kinds of solutions. On the one hand, solutions should come from a change of attitude in managing irrigation water by quantifying the profitability of irrigation water in much the same way as other economic products. In other words, the agricultural sector should seek to maximize the economic return per unit of water (water productivity). On the other hand, irrigation managers should be able to independently assess water demand in relation to the actual availability of water resources and the capability of the conveyance system (a top-down approach). The irrigation managers, provided with appropriate support tools, should be able to determine actual crop water requirements, deciding the optimal timing, quantity and quality of irrigation water to be supplied to the soil in each farm in order to control crop yield, while minimizing deep percolation fluxes of water and nutrients, hence water losses and groundwater degradation.
In this context, the use of decision support systems (DSS) may significantly enhance the management of on-demand irrigation systems. Advanced technical tools are nowadays available for monitoring and simulating the various physical processes involved in an irrigation system. Such tools can be used in an integrated way for estimating irrigation water demand and scheduling water use, and thus for simulating the operation of on-demand irrigation networks, under varying soil and weather conditions, crop management options and irrigation technologies.
In general, by integrating the above technical tools with economic instruments, such systems may facilitate the decision-making process on the quantities of water to be allocated to agricultural consumers. Such tools may take into account the most profitable cropping patterns given any water restrictions imposed by existing hydrological conditions, as well as the potential yields in each irrigation district according to its production characteristics, irrigation efficiency, economic scenario and external factors such as agricultural policies.
Against this background, this paper aims to illustrate the application at district scale of a physically-based agro-hydrological model for identifying optimal irrigation water needs in a context of on-demand irrigation systems. The model, called FLOWS-HAGES, simulates district-level water use and scheduling over the course of an irrigation season. The methodology will be applied for sector 6 of Irrigation District 10 in the “Sinistra Ofanto” irrigation system.
The paper is mostly devoted to illustrating the features of the FLOWS-HAGES model. The interaction of this model with a tool simulating the operation of on-demand irrigation networks by the irrigation manager will be the subject of a follow-up paper where, from a technical point of view, the following aspects of irrigation will be solved: spatial and temporal variability of soil water deficit, spatial and temporal variability of farmers’ water demand, water resources availability and structural capability of the irrigation conveyance and distribution network.

2. Materials and Methods

2.1. The FLOWS-HAGES Model

The model simulates water flow and solute transport in heterogeneous agri-environmental systems (termed FLOWS-HAGES). It produces information on the evolution of:
  • Soil water contents and pressure potentials in the soil profile;
  • Water uptake and actual evapotranspiration (actual water needs);
  • Solute (e.g., nitrates, pesticides, salts, heavy metals) concentrations in the soil profile;
  • Deep percolation water fluxes (return flow to the groundwater) and their quality in terms of solute (solute fluxes);
  • Stress periods for each crop.
The numerical tool is shaped to be run in different configurations and for different purposes:
  • To provide effective water requirement data to be used for optimal management of the irrigation network;
  • To facilitate the decision-making process on the quantities of water to be allocated to agricultural users;
  • To consider the most profitable cropping patterns given water availability restrictions imposed by the existing hydrological systems, and the potential yields reached in each irrigation district according to its productive characteristics, irrigation efficiency, economic scenario and external factors such as agricultural policies;
  • To predict the impact of anticipated climate changes on the irrigation system under the current land use and vegetation cover;
  • To predict the impact of human-driven changes in the land use on the irrigation system under current climate conditions;
  • To predict the impacts on the irrigation system under mixed conditions 4–5.
All the evaluations above can also be carried out in a stochastic framework, thus providing uncertainty bands for each of the main outputs.
The numerical code was written in Matlab and is based on a standard finite difference scheme [4,5,6]. In the model, vertical transient flow is simulated by numerically solving the 1D form of the Richards equation using an implicit, backward, finite differences scheme with explicit linearization.
The model discretizes the spatial flow field in a prescribed number of nodes (usually 100) of constant width (Δz). Time discretization starts with a prescribed initial time increment (Δt). This time increment is automatically adjusted at each time level according to the criteria proposed by Vogel [7].
Flow rates and pressure heads, whether constant or variable over time, can be assumed as the upper boundary condition. Gradients of different value, pressure heads or flow rates, again whether constant or variable, can be assumed at the bottom of the soil profile. In this study, daily potential evapotranspiration and rain were assumed at the upper boundary whereas free drainage was assumed at the bottom boundary.
Details on the water and solute transport equations and the sink term calculations are given in the Appendix. Details of the criteria adopted by the model to control irrigation are reported below.
The criterion used by the model to calculate the time for irrigation and the irrigation volume is schematically summarized in Figure 1a,b.
If hav remains higher than hcrit, this means that on average the pressure head h lies above the stress condition and no irrigation is required (Figure 1a). Irrigation starts any time hav becomes lower than hcrit (Figure 1b). In other words, it is assumed that stress starts when the average pressure head in the root zone becomes lower than the threshold for water stress.
The model evaluates the gross irrigation volume to be supplied at day t to the soil profile, Irrgross, as follows:
I r r g r o s s = 0 D r o o t ( θ f c θ ( z ) ) d z   ( cm   of   water )  
Again, irrigation aims to bring the actual water content at each depth to the water content at field capacity. Figure 2 provides a schematic view of this calculation. Net irrigation, Irrnet, to be supplied on day t is calculated as the difference between Irrgross and the eventual rain falling at the time when the irrigation starts.
I r r n e t = I r r g r o s s r a i n f a l l   ( cm   of   water )  
If the rainfall exceeds Irrgross, irrigation water is not supplied.
At the irrigation network management level, the actual irrigation amount (cm), Irract, may be obtained by dividing Irrnet by the on-farm irrigation efficiency, IE, of the irrigation system used for the crop considered, to give:
I r r a c t = I r r n e t / I E   ( cm   of   water )  
The duration of a hydrant opening, tirr, to satisfy a given Irract at day t is
t i r r = I r r a c t   A   k / q
q being the nominal discharge of the hydrant considered (l/s), A (ha) the area to be irrigated and k a factor to account for the different dimensions in the equation.
In practice, by repeating the same reasoning for each of the fields included in the irrigation sector, the FLOWS-HAGES model identifies the number, time and position of hydrants to be in operation at day t. By aggregating the irrigation demand (the Irract) for each field in the sector for day t, the model calculates the overall discharge to be delivered in that day (the discharge to be delivered at the sector head).

2.2. Sector 6 of the Capitanata Irrigation Network

The study area belongs to the Capitanata Irrigation Network (in Italian, Consorzio per la Bonifica della Capitanata CBC), Apulia, southern Italy) (Figure 3). The CBC is currently managed by a local water users’ organization (WUO), called Consorzio per la Bonifica della Capitanata and covers an area of 22,500 ha. The area is operated with a pressurized on-demand delivery schedule and with a “low zone,” where water is supplied to farms by gravity, and a “high zone” where farms are at higher elevations than the water source and irrigation water is conveyed and supplied by a lifting plant. The CBC-low zone consists of seven districts, each sub-divided into smaller operational units (“sectors”) consisting of several grouped farms.
District 10 of the low zone was specifically considered in this study, covering an overall area of about 2000 ha, of which the total irrigable area is 1679 ha and the area currently irrigated is 1423 ha. Irrigation water comes from a storage and compensation reservoir, with a total capacity of about 50,000 m3, which in turn is supplied with water by the Capacciotti dam. The hydraulic scheme of the low zone is presented in Figure 3b.
In this paper, the model was applied to Sector 6 (Figure 3c). The latter has an area of 139 ha with an average farm size ranging from 0.5 to 2.5 ha. Agriculture in the area is intensive and highly market-oriented. However, the average yearly precipitation in the area is about 550 mm, which is poorly distributed throughout the year. Consequently, intensive farming mostly depends on irrigation.

2.3. Soil Characterization in Sector 6

In order to define the pedological and hydraulic characterization of soils, five sites were selected. A pit was dug (Figure 4) for each profile and a pedological description was carried out. For each soil horizon, disturbed soil samples were collected to determine soil texture. Particle-size analysis was carried out using the hydrometer method combined with sieve analysis to characterize the range of particle diameter from 2 up to 2000 μm [8]. Sand, silt, and clay contents were expressed as percentages by mass of the fine-soil fraction (<2 mm) and soil texture was identified according to the United States Department of Agriculture classification (USDA).
Pedological descriptions and textural data for each soil profile are given in Table 1.
Hydraulic properties of each horizon for each profile were determined by both field and laboratory methods. As for the shallower horizon, field infiltration measurements were carried out using a tension disc infiltrometer from Soil Measurement Systems (Tucson, AZ) [9,10,11] with a 20 cm-diameter disc, covered with a highly permeable nylon membrane.
Each sequence of infiltration runs at a given location started with the lowest supply potential, increasing stepwise to the highest. Care was taken when setting up the instrument not to disturb the surface on which a thin layer of moist fine sand was applied at each measurement location as a patch with a diameter equal to that of the infiltrometer disc. This smoothed out any irregularities of the soil surface and assured good contact between infiltrometer and soil.
The infiltrometer supply pressure heads were −10, −7 and −4 cm (water column). For each pressure head, the infiltration test lasted about 13–15 minutes. Before each infiltration measurement, soil samples were taken with steel cylinders near the measurement location to obtain the initial water content. Immediately after each measurement, another sample was collected directly below the infiltration disc to quantify the final water content, as well as bulk density.
In order to deduce simultaneously the water retention and hydraulic conductivity curves from the same infiltration experiment, the transient experimental infiltration data, with initial and final water content, were used for a parameter estimation method based on a numerical solution of the Richards equation [12]. Specifically, in this study, the infiltrometer data were elaborated by the DISC code [13].
At the deeper soil horizons, laboratory hydraulic characterization was carried out on undisturbed soil samples collected by using stainless steel cylinders (inner diameter of 7.6 cm and height of 7 cm). Before performing the hydraulic tests in the laboratory, the height of each core was reduced to the value of 4.0 cm. Each soil core was slowly wetted from below at different steps to obtain the saturation conditions. We then measured water content θs at h = 0 by the gravimetric method and hydraulic conductivity Ks at θs by the falling-head method [14]. The θ(h) data points were measured using a sand-kaolin suction table. Water retentions were obtained at the following pressure heads: 1.0, 3.0, 7.0, 10.0, 15.0, 30.0, 70.0, 90.0, 130.0 and 180.0 cm (referring to the middle of the sample). We assumed that the average water contents actually corresponded to the indicated pressure head values [15].
The disturbed uppermost part of each soil core was also used to determine the water content θ at h values of 30, 60 and 120 m by the pressure membrane apparatus. The hydraulic properties measured in the five soil profiles were extended to all the fields in sector 6 by using a simple Thiessen polygon technique.

2.4. Crop Distribution and Actual Irrigation Volumes (Year 2016)

The most important species grown in sector 6 are grapevines, peaches and cereals. Figure 5 shows the cultivated area and relative spatial distribution for the year 2016 (data were obtained from CBC). Vineyards cover almost half (46%) the total study area; fruit trees with mainly early peach account for one fifth of the total area (21%); table grapes, apricots and olives are present in a smaller fraction (2%).
From the recorded volumes in Sector 6 at hydrant level for the year 2016, the total water consumption of farmers was about 125,000 m3. The distribution network of Sector 6 supplies water to 31 delivery hydrants, and covers a total cultivated area of 129 ha, of which 64 ha of vineyards consumed, according to CBC, more than 73,000 m3 in 2016 (Table 2).
However, the 30 ha of peach trees consumed more than 46,000 m3, amounting to 38% of total water consumption. N/A means Not Available

2.5. Evaluating Model Simulations by Direct Water Content Measurements

Model simulation reliability was evaluated by comparing simulation results with water storage monitored in the root zone in selected fields of the district by a proximal sensor based on the electromagnetic induction technique (EMI). EMI proximal sensors allow extensive surveys [16]. The EMI sensors measure depth-weighted electrical conductivity (the so-called apparent electrical conductivity, ECa), integrated up to a soil depth depending on both the type of sensor used and the configuration of the sensor coils. ECa depends on the soil profile distribution of water electrical conductivity, σw, volumetric water content, θ, tortuosity, τ, of the soil pore system, as well as on other factors related to the solid phase such as bulk density, clay content and mineralogy. Consequently, separating the effect of single soil properties (e.g., soil water) on ECa is no simple task (see for example, [17,18]). However, in two separate studies, [19,20] found spatial variation in soil water stored within the top 0.5 and 1.7 m to be highly correlated with the spatial variation in ECa measured with EMI sensors. Similar results were obtained by other authors (see for example, [21]. These studies demonstrated that in soils with low salinity, measurements of soil ECa can be used for estimating soil water stored in the soil profile. Estimating soil water content using ECa readings requires a site-specific calibration relating EMI measurements to simultaneous water contents measured by a different standard reference method. In this sense, the time-dynamics of water storage in the selected fields was measured by a Diviner 2000 (Sentek Pty Ltd, Stepney, South Australia) capacitance sensor. A Diviner access tube was installed in three sites of each field up to a depth of 1 m. The Diviner 2000 is a capacitance device consisting of a single sensor fitted to a square rod which can be lowered in a PVC plastic access tube (5.1-cm “i.d.”, 5.6-cm “o.d.”) [22].
In the same three calibration sites, the ECa was measured in both horizontal coplanar (HCP) and vertical coplanar (VCP) orientation by a CMD-Miniexplorer (GF Instruments, Brno, Czech Republic). The CMD operates at 30 kHz frequency and has three receiver coils with 0.32, 0.71 and 1.18 distance from the transmitter coil. The sensor has to be rotated 90° to change the orientation from VCP to HCP. According to the manufacturer, the sensor has an effective depth range of 0.5, 1.0 and 1.8 m in the HCP model, which is halved (0.25, 0.5 and 0.9 m) in the VCP mode. Thus, a single measurement allows six different ECa readings, each corresponding to a different depth exploration. Both the water content and ECa measurements were carried out at three times campaigns.
We tested which of the six ECa readings yields the best correlation with the total soil water storage in the 0–1 m depth soil layer. It was found that the 0.9 for VCP and 1.0 for HCP performed best. A simple linear regression accurately predicted the water stored in the profile (average of the three measurement sites in each field) at any time. Once calibrated, the CMD was used to monitor the evolution of the water storage in the root zone (0–80 cm) averaged over the whole field in a total of 12 monitoring campaigns.

3. Results and Discussion

In this paper, the FLOWS-HAGES model was applied twice for each farm: (1) to simulate the dynamics and distribution of water in the soil profile by assuming the irrigation time and volumes actually supplied by the farmer; (2) to simulate the dynamics and distribution of water in the soil profile by assuming that irrigation time and volumes were not decided by the farmer but calculated by the model. Water volume requirements from simulations—assumed to be the optimum—were compared to volumes actually delivered by the consortium, thus allowing an estimation of the irrigation efficiency both at field and district scale. Hereafter, labels MOD and ON-DEM will be used to identify respectively the results of the two simulations. Both simulations were carried out by assuming as natural top boundary conditions the rainfall and evapotranspiration data for the year 2016 recorded by a meteorological station located within the sector considered (data not shown here).

3.1. Irrigation Volumes Measured and Calculated by the Model

Figure 6a,b shows, respectively, the cumulative irrigation volumes for both grapevines and peaches, the most representative crops in the area. In the figure, the curves represent the irrigation volumes actually supplied by the farmers. The bold red line represents the curve calculated by the model. Changes in the slope of each curve indicate irrigation events.
What is immediately clear from the figures is the large variability in farmer behavior in irrigating the same crop. Besides the lowest volumes, which would suggest a likely use of private farm wells, the large range of cumulative irrigation values indicates that the farmer’s experience may be a rather misleading and erratic criterion to establish both the time of irrigation and the irrigation volumes. It is worth noting that the calculated irrigation volumes approximately represent the average behavior of farmers, assuming no variability in soil (or hydraulic behavior) in the irrigation sector.
Further details will be given below for some fields representing both crops, for the case of either lower or higher cumulative irrigation volumes than those calculated by the model. For these selected fields, water storage in the root zone (the 0–80 cm layer) simulated by the model were compared to those measured by using the EMI sensor.
Of course, if it is assumed that the model calculates irrigation times and volumes correctly, when the irrigation calculated by the model is lower than that applied by the farmer, this means that the farmer is using more water than necessary. By contrast, if the model estimates the need for more irrigation water than that applied by the farmer, then the farmer is inducing some stress, whose intensity will depend on the water volume actually applied by the farmer.

3.2. Irrigation Volumes, Pressure Heads in the Root Zone and Deep Percolation Fluxes

Regarding the grape crop, the farms identified as FARM 2-3-G and FARM 52-1-G were considered in our analysis of the cases of cumulative irrigation volumes lower and higher than those calculated.
The graphs in Figure 7a–d illustrate the irrigation volumes (both per single day and cumulative) calculated by the model (a), the irrigation volumes actually supplied by the farmer (b), the pressure head at 40 cm depth (c) and the percolation fluxes at 80 cm (d) for FARM 2-3-G. The depth of 80 cm was considered as the lowest depth for root uptake. In all the graphs, red and blue lines are used to identify respectively the MOD and ON-DEM simulation results. To demonstrate the goodness of the simulation results, the graph in Figure 7e compares the time evolution of the water storage predicted by the model and that estimated from the EMI sensor ECa readings.
Figure 8a–e shows the same graphs for FARM 52-1-G.
As for the peach crop, the farms identified as FARM 54-4-P and as FARM 49-5-P will be considered to analyze the cases of cumulative irrigation volumes lower and larger than the calculated ones. Figure 9a–e and Figure 10a–e show the graphs for the FARM 54-4-P and the FARM 49-5-P, respectively.
Irrespective of whether the farmer supplied more or less than the irrigation volumes calculated by the model, compared to the farmer, the model “would irrigate” more frequently with lower volumes.
In the graph of pressure heads, the horizontal line indicates the critical pressure head (hcrit) selected for the grape crop. It is interesting to see how the model allows the pressure head to remain regularly above the critical value, with frequent irrigation and low volumes applied. The low volumes have major implications in terms of overall discharge to be supplied at the head of the irrigation sector.
The MOD deep percolation fluxes may be either higher or lower than the ON-DEM, depending on the value of the cumulative irrigation volumes calculated by the model or actually supplied by the farmer. However, no general rule exists. For example, in the case of FARM 52-1-G, MOD irrigation volumes are higher than ON-DEM volumes, and yet deep percolation fluxes are comparable to those of ON-DEM. This may be explained by the fact that the farmer irrigates less frequently with higher volumes. This induces stress periods alternating with relatively large irrigation supplies, the latter leading to higher water contents (higher than the field capacity) and hence higher downward water fluxes due to the higher hydraulic conductivity.

3.3. List of Hydrants to be Opened and Efficiency Indices

The above discussion concerning some representative fields may be extended to the whole set of fields considered in the irrigation sector. The final result is the list of hydrants to be opened day by day during the growing season of a crop, along with irrigation volumes. Table 3 reports this list for the grape crop.
Based on the modeling results, the following indices may be mapped (see Figure 11a–d).
1. Ratio of irrigation volumes model/irrigation volumes on-demand: A value >1 indicates that the irrigation volumes modeled exceed those supplied by the farmer. If it is assumed that the model calculates irrigation volumes correctly, this means that the farmer is causing some stress, which may have an impact on crop yield. By contrast, a value <1 means higher volumes applied by the farmer, thus indicating that the farmer is using more water than necessary and is causing undue water (and nutrient) losses.
2. Ratio of drainage volumes model/drainage volumes on-demand: this index gives a measure of the deep percolation fluxes in the case of MOD irrigation volumes compared to ON-DEM. The index tends to follow the same pattern as the first. However, this cannot be generalized as there are fields where considerable MOD irrigation volumes result in MOD deep percolation fluxes which are somewhat greater than those of ON-DEM, due to the higher irrigation frequency and lower volumes applied in the MOD case. In other words, this means that when the model calculates the irrigation volumes, the water content in the root zone is maintained low enough to keep the hydraulic conductivity low, thus reducing the fluxes downward. However, the water content is kept stably above the critical pressure head so that the crop can stably uptake the optimal amount of water without stress.
The last two indices give information on the fraction of the water lost by percolation fluxes below the root zone in the MOD and ON-DEM cases, respectively. Interestingly, in general with the MOD irrigation system, a lower fraction of water is lost to drainage. In other words, of the water applied, a larger part remains in the root zone available for root uptake. This is what is commonly called application efficiency, which is evidently related to the frequency of irrigation and volumes applied for each irrigation, an issue extensively discussed above.

4. Conclusions

The main objective of this paper was to describe the potential of an irrigation management tool to simulate district-level water use and scheduling over the course of an irrigation season. The tool is largely based on the use of mathematical models for the simulation of soil water movement in the unsaturated zone. This is a precise approach to rational irrigation and the key to saving irrigation water. Compared to other similar existing approaches, the soil water balance is calculated by solving the partial differential equation known as the Richards equation. Additionally, the tool takes into account the spatial variability of the soil hydraulic parameters.
The combination of a detailed mathematical model, which solves the Richards equation, with some method of handling the spatial variability is highly desired. In this sense, the model provides a way to establish the right frequency and amount of irrigation to keep the pressure head above the critical value by minimizing both the stress periods and the water (and nutrients) lost by deep percolation fluxes.
The tool, due to its increased requirements in specialized data, is not aimed at the farmer but rather at an organization that could provide irrigation advisory services to farmers. It is designed to support farmers in improving irrigation management at the farm scale.
The tool responds to the question of balancing the pressure from growing water demands and the need to increase agricultural resilience to climate change through irrigation. Reducing water losses in the agricultural sector, by far the biggest water consumer, is probably the single most important step toward a more efficient management of water resources. To achieve this goal, the first step is undoubtedly to empower managers and through them the farmers.

Author Contributions

Conceptualization, A.C.; Data curation, G.D., N.N. and L.N.; Formal analysis, A.C. and A.B.; Investigation, G.D., A.S. and A.C.; Methodology, A.C. and N.L.; Software, A.C.; Validation, G.D., A.S. and A.B.; Visualization, A.S.; Writing—review & editing, A.C., A.B., A.S. and G.D.

Funding

This research received no external funding.

Acknowledgments

This research received no external funding. It was carried out thanks to soil and meteorological data, as well as data on irrigation volumes, provided by the Consorzio di Bonifica della Capitanata, Foggia, Italy. We are especially grateful to Pandi Zdruli of the IAMB, as well as to Osvaldo Moro and his collaborators of the CBC for their work on the field pedological and hydrological experiments and investigations.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. FLOWS-HAGES: Short Model Description

Appendix A.1. Water Flow

One-dimensional vertical transient water flow in this model is simulated by numerically solving the 1D form of the Richards equation (Equation (A1)) using an implicit, backward, finite differences scheme with explicit linearization similar to that adopted in the Soil-Water-Atmosphere-Plant (SWAP) model [23]:
C ( h ) h t = z ( K ( h ) h z K ( h ) ) S w ( h )
where C(h) = /dh (L−1) is the soil water capacity, h (L) is the soil water pressure head, t (T) is time, z (L) is the vertical coordinate being positive upward, K (h) (L·T−1) the hydraulic conductivity and Sw(h) (T−1) is a sink term describing water uptake by plant roots.

Appendix A.2. Solute Transport

The advection–dispersion equation (ADE) is used for prediction of solute transport (agrochemicals, salts, heavy metals, …):
θ C t + ρ b Q t = q C z + ( θ D C z ) z S s
Including non-linear adsorption, linear decay and proportional root uptake in unsaturated/saturated soil. In the equation, C(M·T−3) and Q (M·M−1) are the amount of solute in the liquid and adsorbed phases, respectively, q (L T−1) is the darcian flux, ρb (M L−3) is the bulk density, D (L2 T−1) the hydrodynamic dispersion coefficient. Ss is a source-sink term for solutes. Hydrodynamic dispersion is related to the molecular diffusion constant of the substance in bulk water, D0 (L2 T−1), and the pore water velocity, v = q/θ, as:
D = λ v + η ( θ ) D 0
where λ (L) is the dispersivity and η a tortuosity coefficient. However, as discussed by several researchers ([24,25], among others) the contribution of diffusion to the hydrodynamic dispersion D is often very small. Accordingly, in this study λ was simply derived from the ratio D/v.

Appendix A.3. Hydraulic Properties

Several water retention and hydraulic conductivity functions can be selected in the model. Water retention: Unimodal van Genuchten model [26]:
{ S e = θ θ r θ s θ r = [ 1 + | α V G h | n ] m ,   h < 0 θ = θ s ,   h 0
where h is the pressure head (h ≤ 0), Se is effective saturation and θ is the water content (θs and θr are the water content at h = 0 and for h→∞, respectively). αVG (cm−1), n and m = 1 − 1/n are shape parameters. Note the effective saturation, Se, is to be considered as a cumulative distribution function of pore size with a density function f(h) which may be expressed by:
f ( h ) = d S e d h
In natural soils, the presence of aggregates frequently results in a retention function curve having at least two points of inflection. To represent such behavior, a double porosity approach can be used which assumes that the pore space from θr to θs consists of two fi(h) distributions obtained by Equation (A5), each occupying a fraction ϕi of that pore space. Bimodal Durner model [27]:
S e = φ i [ 1 1 + ( α i h ) n i ] m i ,   0 < φ i < 1   and   φ i = 1     i = 1 , 2
in which ϕi is the weighting of the total pore space fraction to be attributed to the ith subcurve, and αi, ni and mi still represent the fitting parameters for each of the partial curves. Bimodal Ross and Smettem model [28]:
S e = φ 1 ( 1 + α 1 h ) e x p ( α 1 h ) + φ 2 [ 1 1 + ( α 2 h ) n 2 ] m 2 , 0 < φ i < 1   and   φ i = 1       i = 1 , 2
where the following simple one-parameter function:
S e = ( 1 + α h ) e x p ( α h )
describes macroporosity in aggregated soils.
Hydraulic conductivity: Mualem’s model is based on the capillary bundle theory [29] and relates relative hydraulic conductivity, Kr, to f(h) by the equation:
K r ( h ) = K ( h ) K 0 = S e τ [ η ( h ) / η ( 0 ) ] 2 ,   η ( h ) = h 1 f ( h ) d h
in which K0 is the hydraulic conductivity at h = 0, and τ is a parameter accounting for the dependence of the tortuosity and the correlation factors on the water content.
For the case of the van Genuchten unimodal soil water retention model and assuming m = 1 − 1/n, Equation (A9) becomes [26]:
K r ( S e ) = K ( S e ) K 0 = S e τ [ 1 ( 1 S e 1 m ) m ] 2
For the bimodal case, Equation (A10) becomes:
K r ( h ) = K ( h ) K 0 = S e τ [ i = 1 2 φ i η i ( h ) / i = 1 2 φ i η i ( 0 ) ] 2
with parameters as defined earlier.
In the case of bimodal Durner water retention, the Equation (A11) becomes [30]:
K r ( S e ) = K ( S e ) K 0 = ( i = 1 2 φ i S e , i ) τ { i = 1 2 φ i α i [ 1 ( 1 S e , i 1 / m i ) m i ] i = 1 2 φ i α i } 2
For the Ross and Smettem model, the η(h)functions in Equation (A9) are [28]:
η 1 ( h ) = α R S e x p ( α R S h ) ,         η 2 ( h ) = α V G n B ( S e 1 / m ,   m + 1 / n ,   1 1 / n )
where B( ) is the incomplete beta function.
Equation (A11) assumes that all pores in the material can interact [28].
In this case, an instantaneous equilibrium (i.e., instantaneous exchange of water) between the two pore systems is implicitly assumed. Nevertheless, for some materials part of the pore space may bypass the rest. This is the case when, for example, the cutans on the walls of the aggregates reduces the water transport between intra and inter-aggregate spaces. The conductivities of the two porous spaces are in parallel and, for such independent distributions, the following equations apply:
k ( h ) = i = 1 2 k 0 i k r i ( h ) ,       k r i ( h ) = S e τ [ p i ( h ) / p i ( 0 ) ] 2

Appendix A.4. Root Uptake

According to the approach of Feddes [31], the (i) root system distribution and (ii) the activity of the roots along the root zone during the growth season, are embedded in the macroscopic sink term as follows:
(i) The potential root water uptake flux per unit depth at a specific depth, Sp (d−1), is simulated macroscopically by distributing potential transpiration, Tp (cm·d−1), over the root zone depth, Dr (cm), on the basis of a normalized root density distribution, g(z) (cm−1), with depth z. The normalized root density distribution may be obtained by normalizing the root length density distribution Rld (cm·cm−1) by its integral across the rooting depth, Dr:
g ( z ) = R l d ( z ) 0 D r R l d ( z ) d z
Rld may be calculated as the ratio of the total length, Lr(z), of roots in a sample to the sample volume. Hereafter, we will define the denominator of Equation (A2) as root length density integral, RInt. g(z)dz is the fraction of roots located between z and z + dz. In the macroscopic approaches to root uptake, the function g(z) distributes the potential transpiration rate, Tp, through the root zone in proportion to the root distribution [31,32]:
S p ( z ) = g ( z ) T p
with
0 D r g ( z ) d z = 1
and thus
T p = 0 D r S p ( z ) d z
(ii) Low water contents and/or the presence of soluble salts in the soil lower the total hydraulic head and may reduce the water fluxes to the roots, thus reducing root activity and water uptake. Reduction coefficients to decrease the maximum water uptake according to the water and osmotic stresses may be calculated independently and multiplied to calculate the actual root uptake as:
S a = α r w ( h ) α r s ( h o ) S p = α r w ( h ) α r s ( h o ) g ( z ) T p
with αrw and αrs being reduction factors depending on the local (at a given z) water pressure head, h (cm) and osmotic head, ho (cm), respectively.
Accordingly,
T a = 0 D r S a ( z ) d z
and
T a T p = 0 D r g ( z ) α r w ( h ) α r s ( h o ) d z = β ( h , h o )
with Ta being the actual transpiration rate and β a dimensionless water stress index integrated over the whole rooted profile [33,34], providing a measure of total plant stress. A value of β equal to 1 indicates that there is no stress in the soil root zone and that the actual transpiration rate Ta is equal to the potential transpiration rate Tp.
When the soil is irrigated by keeping soil water content under optimal conditions, the eventual reduction in root uptake may only be induced by osmotic stress. Under only osmotic stresses, αrw = 1 and root uptake parameterization reduces to finding the factor αrs, depending on the osmotic potential (ho) induced by salts in the soil water.

Appendix A.5. Calculating the Water Stress Reduction Factor αrw

Feddes and Raats [32] reviewed various functional forms for αrw that have been proposed over the years. We mention here two general model types that have been used most often: piecewise linear functions and continuous smooth functions. To describe water stress, Feddes [31] proposed a piecewise linear reduction function parameterized by four critical values of the water pressure head, h4 < h3 < h2 < h1:
α r w ( h ) = { h h 4 h 3 h 4 ,   h 3 > h > h 4 1 ,       h 2 h h 3 h h 1 h 2 h 1 ,   h 1 > h > h 2 0 ,       h h 4   or   h h 1
Uptake is at the potential rate when the pressure head is h 3 h h 2 , drops off linearly when h > h2 or h < h3, and becomes zero when h h 4   o r   h h 1 . In general, the value of h3 is expected to be a function of evaporative demand.
van Genuchten [35] proposed an alternative smooth, S-shaped reduction function to account for water stress:
α r w ( h ) = 1 1 + ( h h 50 ) p 1
where h50 and p1 are adjustable parameters. Equation (A23) was motivated by a study [36] that found that an S-shaped function described salt tolerance yield reduction data better that other functions, including the threshold-slope model of Maas and Hoffman [37].
van Genuchten [35] reasoned that if water and salinity stress are assumed to have similar effects on yield, and hence on transpiration, then the S-shaped function should also be applicable to water stress response data.

Appendix A.6. Calculating the Salinity Stress Reduction Factor αrs

Despite the large body of studies on the topic (see amongst others [32,36,38,39,40,41]), proper modeling and parameterization of root water uptake as a function of salinity stresses remains a major challenge. The process is mainly determined by the specific soil-water-salt conditions and distributions which become established in the root zone according to the local soil physical-hydrological characteristics. This is why the well-known salt tolerance database presented by the U.S. Salinity Laboratory [37] should be handled with prudence in numerical modeling as most of those data were collected under well-controlled experimental conditions in which uniform salt distribution over the root zone was established. As confirmed by several researchers [34,38,39,40,41], the reduction function to be used in numerical modeling should be determined by analyzing the dynamics of water uptake under transient salinity conditions.
Several functional forms have been proposed for uptake reduction due to salinity (see amongst others [32,36]). The response function can be written in terms of concentration, or electrical conductivity of either the soil water or the soil saturation extract, or osmotic pressure head [36,37,42,43,44].
In terms of osmotic pressure heads, ho, consistent with the Mass and Hoffman [37] model for crop salt tolerance (Equation (A6)), the effects of salinity stress on root water uptake can be described using the piecewise linear (threshold-slope) function:
α ( h o ) = { 1 ,         a h o 0 1 + b ( h o a ) ,   a > h o > a 1 b 0 ,         h o a 1 b
where a and b are the adjustable parameters, often referred to as the salinity threshold and slope, respectively, thus mirroring the terminology used for the Maas–Hoffman parameters A and B. Note, however, that the parameter sets are not the same: A and B parameterize total yield reductions as a function of average root zone salinity, whereas a and b parameterize local reductions in the root water uptake rate as a function of osmotic head.
Similarly, in correspondence to Equation (A4), a van Genuchten-Hoffman type model) for the effects of salinity stress on root water uptake can be described as an S-shaped salinity-stress function:
α r s ( h o ) = 1 1 + ( h o h o , 50 ) p 2
where p2 and ho,50 are the adjustable parameters, the latter being the osmotic pressure head where uptake is halved.
The osmotic potential, expressed as osmotic head ho, is assumed to be a linear function of soil solution salinity ECw according to U. S. Salinity Laboratory Staff [45]:
h o = 360 E C
Here, 360 is a factor to convert the salinity-based values to cm osmotic head.

Appendix A.7. Combined Water and Salinity Stress

A major challenge is how to combine the effects of water and salinity stress. Uptake reductions due to a combination of water and salinity stresses could be modeled by assuming that the stresses are somehow additive or multiplicative. A generalized model for additive stresses is obtained by assuming that water uptake occurs in response to some weighted sum of the soil water pressure and osmotic heads:
α r c ( h , h o ) = 1 1 + [ α r w h + α r s h o h o 50 ] p 2
with simple additivity resulting when αrw = αrs = 1. Although additivity has been inferred from a number of laboratory and field experiments (e.g., [46,47,48,49,50]), its general applicability remains uncertain [51], especially for field conditions subject to relatively wide ranges in pressure heads (wetting/drying cycles).
Linear or weighted additivity is only one of several possibilities for combining the effects of water and salinity stress. Another approach presumes that water and salinity effects are multiplicative. In the general case this leads to:
α r c = α r w α r s ,
where αrw and αrs can be calculated as described above. In this case, Equation (A27) becomes:
α r c ( h , h o ) = 1 1 + ( h h 50 ) p 1 × 1 1 + ( h o h o 50 ) p 2
where h50, p50, p1 and p2 are the presumably crop, soil and climate-specific parameters.

Appendix A.8. Root Density Distribution

Several root density distributions, g(z), may be dealt with in the model for simulating the sink term in Equation (A1), assuming root distributions to be either homogeneous [31] or variable with depth [52,53,54], the latter accounting for the fact that in a moist soil the roots can mainly extract water from the upper root zone layers.
A constant distribution over the whole root zone would reduce simply to [31]:
g ( z ) = 1 D r
The model may also apply a modified version of the constant distribution in the case of layered root systems (see [55]).
A Prasad-type triangular g(z) function [52], with root water uptake at the bottom of the root zone (Dr) equal to zero:
g ( z ) = 2 ( D r z ) D r 2
A logistic root density distribution may also be used in the code
g ( z ) = a b c   e x p ( c z ) b + e x p ( c z ) 2
with a, b and c being empirical parameters to be estimated by fitting the logistic to measured values of g(z).

Appendix A.9. Solute Sink Term

The solute sink term Ss is:
S s =   μ ( θ C + ρ b Q ) R p f S w C
where μ is the first order rate coefficient of transformation (T−1), Rpf is the root uptake preference factor (dimensionless) accounting for positive or negative selection of solute ions relative to the amount of soil water that is extracted [23].

References

  1. Blinda, M.; Thivet, G. Ressources et demandes en eau en Méditerranée: Situation et perspectives. Sécheresse 2009, 20, 9–16. [Google Scholar] [CrossRef]
  2. Hsiao, T.; Steduto, P.; Fereres, E. A systematic and quantitative approach to improve water use efficiency in agriculture. Irrig. Sci. 2007, 25, 209–231. [Google Scholar] [CrossRef]
  3. Lamaddalena, N.; Sagardoy, J. Performance Analysis of On-Demand Pressurized Irrigation Systems; Food & Agriculture Organization of the United Nations: Rome, Italy, 2000. [Google Scholar] [CrossRef]
  4. Coppola, A.; Comegna, V.; Basile, A.; Lamaddalena, N.; Severino, G. Darcian preferential water flow and solute transport through bimodal porous systems: Experiments and modelling. J. Contam. Hydrol. 2009, 104, 74–83. [Google Scholar] [CrossRef]
  5. Coppola, A.; Randazzo, L.A. MatLab code for the transport of water and solutes in unsaturated soils with vegetation. In Tech. Rep. Soil and Contaminant Hydrology Laboratory; University of Basilicata: Potenza, Italy, 2006; Unpublished work. [Google Scholar]
  6. Coppola, A.; Randazzo, L.; Basile, A.; Fenu, C. FLOWS-HAGES (FLOw of Water and Solutes in Heterogeneous AGri-Environmental Systems): A MatLab Code for the Transport of Water and Solutes in Unsaturated Soils with Vegetation in Heterogenous Systems; D. S. U. o. B. Tech. Rep. Soil and Contaminant Hydrology Laboratory: Potenza, Basilicata, Italy, 2018; in press. [Google Scholar]
  7. Vogel, T. SWM II-Numerical Model of Two-Dimensional Flow in a Variably Saturated Porous Medium; Wageningen Agricultural University: Wageningen, The Netherlands, 1988. [Google Scholar]
  8. Gee, G.W.; Or, D. Particle-size analysis. In Methods of soil analysis, Part4. Physical methods. Soil Sci. Soc. Am. Book Ser. 2002, 5, 255–293. [Google Scholar]
  9. Ankeny, M.; Ahmed, M.; Kaspar, T.; Horton, R. Simple Field Method for Determining Unsaturated Hydraulic Conductivity. Soil Sci. Soc. Am. J. 1991, 55, 467–470. [Google Scholar] [CrossRef]
  10. Coppola, A.; Basile, A.; Wang, X.; Comegna, V.; Tedeschi, A.; Mele, G.; Comegna, A. Hydrological behaviour of microbiotic crusts on sand dunes: Example from NW China comparing infiltration in crusted and crust-removed soil. Soil Till. Res. 2011, 117, 34–43. [Google Scholar] [CrossRef]
  11. Basile, A.; Buttafuoco, G.; Mele, G.; Tedeschi, A. Complementary techniques to assess physical properties of a fine soil irrigated with saline water. Environ. Earth Sci. 2012, 66, 1797–1807. [Google Scholar] [CrossRef]
  12. Lazarovitch, N.; Ben-Gal, A.; Šimůnek, J.; Shani, U. Uniqueness of soil hydraulic parameters determined by a combined Wooding inverse approach. Soil Sci. Soc. Am. J. 2007, 71, 860–865. [Google Scholar] [CrossRef]
  13. Šimůnek, J.; van Genuchten, M.T.; Sejna, M. The DISC Computer Software for Analyzing Tension Disc Infiltrometer Data by Parameter Estimation; Version 1.0, US Salinity Laboratory: Riverside, CA, USA, 2000. [Google Scholar]
  14. Reynolds, D.; Elrick, D. Saturated and field-saturated water flow parameters: Falling head soil core (tank) method. Methods of Soil Analysis, Part 4. Phys. Methods 2002, 5, 809–812. [Google Scholar]
  15. Dane, J.H.; Hopman, J.W. Water retention and storage. In Methods of Soil Analysis: Part 4—Physical Methods; Dane, J.H., Topp, G.C., Eds.; Soil Science Society of America: Madison, WI, USA, 2002; pp. 721–738. [Google Scholar] [CrossRef]
  16. Rhoades, J.D.; Corwin, D.L. Soil electrical conductivity: Effects of soil properties and application to soil salinity appraisal. Commun. Soil Sci. Plant Anal. 1990, 21, 837–860. [Google Scholar] [CrossRef]
  17. Coppola, A.; Smettem, K.R.; Ajeel, A.; Saeed, A.; Dragonetti, G.; Comegna, A.; Lamaddalena, N.; Vacca, A. Calibration of an electromagnetic induction sensor with time-domain reflectometry data to monitor rootzone electrical conductivity under saline water irrigation. Eur. J. Soil Sci. 2016, 67, 737–748. [Google Scholar] [CrossRef]
  18. Dragonetti, G.; Comegna, A.; Ajeel, A.; Deidda, G.P.; Lamaddalena, N.; Rodriguez, G.; Vignoli, G.; Coppola, A. Calibrating electromagnetic induction conductivities with time-domain reflectometry measurements. Hydrol. Earth Syst. Sci. 2018, 22, 1509–1523. [Google Scholar] [CrossRef]
  19. Kachanoski, R.G.; Wesenbeeck, I.V.; Gregorich, E.G. Estimating spatial variations of soil water content using noncontacting electromagnetic inductive methods. Can. J. Soil Sci. 1988, 68, 715–722. [Google Scholar] [CrossRef]
  20. Kachanoski, R.G.; Wesenbeeck, I.V.; Jong, E.D. Field scale patterns of soil water storage from non-contacting measurements of bulk electrical conductivity. Can. J. Soil Sci. 1990, 70, 537–542. [Google Scholar] [CrossRef]
  21. Sheets, K.R.; Hendrickx, J.M. Noninvasive soil water content measurement using electromagnetic induction. Water Resour. Res. 1995, 31, 2401–2409. [Google Scholar] [CrossRef]
  22. Evett, S.R.; Tolk, J.A.; Howell, T.A. Soil profile water content determination. Vadose Zone J. 2006, 5, 894–907. [Google Scholar] [CrossRef]
  23. Van Dam, J.; Huygen, J.; Wesseling, J.G.; Feddes, R.A.; Kabat, P.; Van Walsum, P.E.; Groenendijk, P.; Van Diepen, C.A. Theory of SWAP, Version 2.0; Simulation of Water Flow, Solute Transport and Plant Growth in the Soil-Water-Atmosphere-Plant Environment; Technical Document; Wageningen Agricultural University: Wageningen, The Netherlands, 1998. [Google Scholar]
  24. Comegna, V.; Coppola, A.; Sommella, A. Nonreactive solute transport in variously structured soil materials as determined by laboratory-based time domain reflectometry (TDR). Geoderma 1999, 92, 167–184. [Google Scholar] [CrossRef]
  25. Vanderborght, J.; Vereecken, H. Review of Dispersivities for Transport Modeling in Soils. Vadose Zone J. 2007, 6, 29–52. [Google Scholar] [CrossRef]
  26. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef]
  27. Durner, W. Predicting the unsaturated hydraulic conductivity using multi-porosity water retention curves. Indirect Methods Estim. Hydraul. Prop. Unsaturated Soils 1992, 185–202. [Google Scholar]
  28. Ross, P.J.; Smettem, K.R. Describing soil hydraulic properties with sums of simple functions. Soil Sci. Soc. Am. J. 1993, 57, 26–29. [Google Scholar] [CrossRef]
  29. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef]
  30. Priesack, E.; Durner, W. Closed-form expression for the multi-modal unsaturated conductivity function. Vadose Zone J. 2006, 5, 121–124. [Google Scholar] [CrossRef]
  31. Feddes, R.A.; Kowalik, P.J.; Zaradny, J. Simulation of Field Water Use and Crop Yield; Centre for Agriculture Publishing and Documentation: Wageningen, The Netherland, 1978. [Google Scholar]
  32. Feddes, R.A.; Raats, P.A.C. Parameterizing the soil–water–plant root system. Unsaturated Zone Model. Prog. Chall. Appl. 2004, 6, 95–141. [Google Scholar]
  33. Jarvis, N.J. A simple empirical model of root water uptake. J. Hydrol. 1989, 107, 57–72. [Google Scholar] [CrossRef]
  34. Shouse, P.J.; Ayars, J.E.; Šimůnek, J. Simulating root water uptake from a shallow saline groundwater resource. Agric. Water Manag. 2011, 98, 784–790. [Google Scholar] [CrossRef]
  35. Van Genuchten, M.T. A Numerical Model for Water and Solute Movement in and Below the Root Zone; United States Department of Agriculture, Agricultural Research Service U.S. Salinity Laboratory: Riverside, CA, USA, 1987.
  36. Van Genuchten, M.T.; Hoffman, G.J. Analysis of crop salt tolerance data. In Soil Salinity under Irrigation; Shainberg, I., Shalhevet, J., Eds.; Ecological Studies; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1984; pp. 258–271. [Google Scholar]
  37. Maas, E.V.; Hoffman, G.J. Crop salt tolerance—Current assessment. J. Irrig. Drain. Div. 1977, 103, 115–134. [Google Scholar]
  38. Skaggs, T.H.; van Genuchten, M.T.; Shouse, P.J.; Poss, J.A. Macroscopic approaches to root water uptake as a function of water and salinity stress. Agric. Water Manag. 2006, 86, 140–149. [Google Scholar] [CrossRef]
  39. Homaee, M.; Dirksen, C.; Feddes, R.A. Simulation of root water uptake: I. Non-uniform transient salinity using different macroscopic reduction functions. Agric. Water Manag. 2002, 57, 89–109. [Google Scholar] [CrossRef]
  40. Homaee, M.; Feddes, R.A.; Dirksen, C. Simulation of root water uptake: II. Non-uniform transient water stress using different reduction functions. Agric. Water Manag. 2002, 57, 111–126. [Google Scholar] [CrossRef]
  41. Homaee, M.; Feddes, R.A.; Dirksen, C. Simulation of root water uptake: III. Non-uniform transient combined salinity and water stress. Agric. Water Manag. 2002, 57, 127–144. [Google Scholar] [CrossRef]
  42. Maas, E.V. Salt tolerance of plants. Appl. Agric. Res. 1986, 1, 12–26. [Google Scholar]
  43. Maas, E.V. Crop salt tolerance. In Agriculture Salinity Assessment and Management; Manuals and Reports on Engineering Practice; American Society Civil Engineers: New York, NY, USA, 1990; pp. 262–334. [Google Scholar]
  44. Maas, E.V.; Grattan, S.R. Crop yields as affected by salinity. Agronomy 1999, 38, 55–110. [Google Scholar]
  45. U.S. Salinity Laboratory Staff. Diagnosis and Improvement of Saline and Alkali Soils. USDA Handb.60; U.S. Govenment Printing Office: Washingtin, DC, USA, 1954.
  46. Wadleigh, C.H. The integrated soil moisture stress upon a root system in a large container of saline soil. Soil Sci. 1946, 61, 219–238. [Google Scholar] [CrossRef]
  47. Meiri, A.; Shalhevet, J. Crop growth under saline conditions. In Arid Zone Irrigation; Springer: Berlin/Heidelberg, Germany, 1973; pp. 277–290. [Google Scholar]
  48. Childs, S.W.; Hanks, R.J. Model of Soil Salinity Effects on Crop Growth 1. Soil Sci. Soc. Am. J. 1975, 39, 617–622. [Google Scholar] [CrossRef]
  49. Du Plessis, H.M. Evapotranspiration of citrus as affected by soil water deficit and soil salinity. Irrig. Sci. 1985, 6, 51–61. [Google Scholar] [CrossRef]
  50. Bresler, E.; Hoffman, G.J. Irrigation Management for Soil Salinity Control: Theories and Tests 1. Soil Sci. Soc. Am. J. 1986, 50, 1552–1560. [Google Scholar] [CrossRef]
  51. Shalhevet, J.; Hsiao, T.C. Salinity and drought. Irrig. Sci. 1986, 7, 249–264. [Google Scholar] [CrossRef]
  52. Prasad, R. A linear root water uptake model. J. Hydrol. 1988, 99, 297–306. [Google Scholar] [CrossRef]
  53. Raats, P.A.C. Steady Flows of Water and Salt in Uniform Soil Profiles with Plant Roots 1. Soil Sci. Soc. Am. J. 1974, 38, 717–722. [Google Scholar] [CrossRef]
  54. Vrugt, J.A.; Hopmans, J.W.; Šimunek, J. Calibration of a two-dimensional root water uptake model. Soil Sci. Soc. Am. J. 2001, 65, 1027–1037. [Google Scholar] [CrossRef]
  55. Coppola, A.; Comegna, A.; Dragonetti, G.; Gerke, H.H.; Basile, A. Simulated preferential water flow and solute transport in shrinking soils. Vadose Zone J. 2015, 14. [Google Scholar] [CrossRef]
Figure 1. Graphical view of the criterion used by the model to calculate the time for irrigation and the irrigation volume. (a) hav higher than hcrit, no irrigation is required; (b) hav lower than hcrit, irrigation is required to bring the pressure head at the field capacity, hfc.
Figure 1. Graphical view of the criterion used by the model to calculate the time for irrigation and the irrigation volume. (a) hav higher than hcrit, no irrigation is required; (b) hav lower than hcrit, irrigation is required to bring the pressure head at the field capacity, hfc.
Water 11 00841 g001
Figure 2. Graphical view of the irrigation volume calculation according to the criterion shown in Figure 1.
Figure 2. Graphical view of the irrigation volume calculation according to the criterion shown in Figure 1.
Water 11 00841 g002
Figure 3. The study area (a). The hydraulic scheme of the low-zone (b). Sector 6 of district 10 in the “Sinistra Ofanto” irrigation system (c).
Figure 3. The study area (a). The hydraulic scheme of the low-zone (b). Sector 6 of district 10 in the “Sinistra Ofanto” irrigation system (c).
Water 11 00841 g003
Figure 4. (a) Soil profile distribution; (b) Digging of a soil pit; (c) details of the soil profile vegetation context.
Figure 4. (a) Soil profile distribution; (b) Digging of a soil pit; (c) details of the soil profile vegetation context.
Water 11 00841 g004
Figure 5. Cropping pattern map of Sector 6 of District 10 in the “Sinistra Ofanto”.
Figure 5. Cropping pattern map of Sector 6 of District 10 in the “Sinistra Ofanto”.
Water 11 00841 g005
Figure 6. Cumulative irrigation actually supplied by farmers and calculated by the model (the bold red line) for grape (a) and peach (b).
Figure 6. Cumulative irrigation actually supplied by farmers and calculated by the model (the bold red line) for grape (a) and peach (b).
Water 11 00841 g006
Figure 7. FARM G 2-3. (a) Irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) irrigation volumes actually supplied by the farmer (ON-DEM); (c) pressure head at 40 cm depth; (d) percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the electromagnetic induction technique (EMI) sensor apparent electrical conductivity (ECa) readings (MEAS).
Figure 7. FARM G 2-3. (a) Irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) irrigation volumes actually supplied by the farmer (ON-DEM); (c) pressure head at 40 cm depth; (d) percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the electromagnetic induction technique (EMI) sensor apparent electrical conductivity (ECa) readings (MEAS).
Water 11 00841 g007
Figure 8. FARM G 52-1. (a) The irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) the irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Figure 8. FARM G 52-1. (a) The irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) the irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Water 11 00841 g008
Figure 9. FARM P 54-4. (a) The irrigation volumes (both per single day and cumulative) calculated by the model (MOD.); (b) the irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Figure 9. FARM P 54-4. (a) The irrigation volumes (both per single day and cumulative) calculated by the model (MOD.); (b) the irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Water 11 00841 g009
Figure 10. FARM P 49-5. (a) Irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Figure 10. FARM P 49-5. (a) Irrigation volumes (both per single day and cumulative) calculated by the model (MOD); (b) irrigation volumes actually supplied by the farmer (ON-DEM); (c) the pressure head at 40 cm depth; (d) the percolation fluxes at 80 cm; (e) comparison between the time evolution of the water storage predicted by the model (MOD) and that estimated starting from the EMI sensor ECa readings (MEAS).
Water 11 00841 g010
Figure 11. Maps of four indices: (a) Ratio of irrigation volumes model/irrigation volumes on-demand (IR MOD/IR ON-DEM); (b) ratio of drainage volumes model/drainage volumes on-demand (DR MOD/DR ON-DEM); (c) ratio of drainage volumes model/irrigation volumes on-demand (DR MOD/IR ON-DEM); (d) ratio of drainage volumes on-demand/irrigation volumes on-demand (DR ON-DEM/IR ON-DEM).
Figure 11. Maps of four indices: (a) Ratio of irrigation volumes model/irrigation volumes on-demand (IR MOD/IR ON-DEM); (b) ratio of drainage volumes model/drainage volumes on-demand (DR MOD/DR ON-DEM); (c) ratio of drainage volumes model/irrigation volumes on-demand (DR MOD/IR ON-DEM); (d) ratio of drainage volumes on-demand/irrigation volumes on-demand (DR ON-DEM/IR ON-DEM).
Water 11 00841 g011
Table 1. Pedological descriptions and textural data for each soil profile.
Table 1. Pedological descriptions and textural data for each soil profile.
ProfileHorizonTexture (%)Texture Class (USDA)
SandClaySilt
P1Ap17.9533.7558.30Silty-Clay-Loam
Bw15.2546.2548.50Silty-Clay
Bw25.7538.7555.50Silty-Clay-Loam
Bw36.6341.2552.13Silty-Clay
P2Ap127.0018.7554.25Silty-Loam
Ap233.4818.7547.78Loam
Bw134.7021.2544.05Loam
Bw234.9021.2543.85Loam
P3Ap20.4526.2553.30Silty-Loam
Bw120.0026.2553.75Silty-Loam
Bw219.8031.2548.95Silty-Clay-Loam
Bw314.5031.2554.25Silty-Clay-Loam
P4Ap32.1521.2546.60Loam
C135.5518.7545.70Loam
C240.3018.7540.95Loam
CK339.8518.7541.40Loam
P5Ap25.0033.7541.25Clay-Loam
B/C47.2523.7529.00Loam
Bw128.2533.7538.00Clay-Loam
Bw229.7033.7536.55Clay-Loam
Table 2. Crop water consumption for the year 2016 (Source: CBC).
Table 2. Crop water consumption for the year 2016 (Source: CBC).
CropArea (ha)Water Consumption
(m3) (Year 2016)
Water Consumption
(%)
Vineyard6473,640.360.4
Early peach3046,320.738
Autumn winter cereals20N/AN/A
Vegetable6.5N/AN/A
Fallow5N/AN/A
Table Grape1732.30.6
Apricot11163.71
Olive1N/AN/A
Table 3. List of hydrants to be opened day by day and volumes (m3) to be supplied during the growing season for grape (year 2016).
Table 3. List of hydrants to be opened day by day and volumes (m3) to be supplied during the growing season for grape (year 2016).
Water 11 00841 i001

Share and Cite

MDPI and ACS Style

Coppola, A.; Dragonetti, G.; Sengouga, A.; Lamaddalena, N.; Comegna, A.; Basile, A.; Noviello, N.; Nardella, L. Identifying Optimal Irrigation Water Needs at District Scale by Using a Physically Based Agro-Hydrological Model. Water 2019, 11, 841. https://doi.org/10.3390/w11040841

AMA Style

Coppola A, Dragonetti G, Sengouga A, Lamaddalena N, Comegna A, Basile A, Noviello N, Nardella L. Identifying Optimal Irrigation Water Needs at District Scale by Using a Physically Based Agro-Hydrological Model. Water. 2019; 11(4):841. https://doi.org/10.3390/w11040841

Chicago/Turabian Style

Coppola, Antonio, Giovanna Dragonetti, Asma Sengouga, Nicola Lamaddalena, Alessandro Comegna, Angelo Basile, Nicoletta Noviello, and Luigi Nardella. 2019. "Identifying Optimal Irrigation Water Needs at District Scale by Using a Physically Based Agro-Hydrological Model" Water 11, no. 4: 841. https://doi.org/10.3390/w11040841

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