3.2. Crop Water Requirement and Irrigation Schedule
Table 2 and
Table 3 show, respectively, the crop water requirement (
CWR) and the irrigation scheduling for different irrigation conditions for cotton. Similarly,
Tables S1–S8 (S stands for supplementary data) show the same for sugarcane, maize, sorghum, and sunflower, respectively. In
Table 3,
Tables S2, S4, S6 and S8, the water stress coefficient was estimated by the following set of equations [Equations (8a) and (8b)] [
14]:
where
Ks = water stress coefficient,
Dr = root zone depletion,
p = fraction of
TAW that is available at the root zone to be used by the crops,
RAW = readily available water, and
TAW = total available water.
The water stress coefficient (
Ks) was used to calculate the adjusted crop evapotranspiration under non-standard conditions for different crops while simulating the irrigation scheduling. In
Table 3,
Tables S2, S4, S6 and S8, the adjusted crop evapotranspiration under non-standard conditions was estimated as (Equation (9)) [
14]:
where
ETa = adjusted crop evapotranspiration under non-standard conditions and
ETc = crop evapotranspiration under standard conditions.
Table 2 indicates that cotton requires fewer irrigation operations for the period August–October, while a significant amount of water is required during November–January. The temporal trend of
Kc is uniform during the initial stage, has an increasing slope in the crop development stage, is still uniform in the mid-season stage, and has a decreasing slope at the time of harvesting. The planting date for cotton is 03/08/2005, but the irrigation is performed on 28/11/2005, depending on the timing of the irrigation and application of irrigation water. The water stress criteria are not significant here, and the crop can beneficially extract the water it needs. Effective rainfall satisfies almost 46% of the total water requirement of cotton, and no irrigation operation was done for the months of August and September. Frequent irrigation should be done during the mid- and late-season stages, mainly in the late-season stage as there was no rainfall. The maximum crop coefficient was achieved at the mid-season stage, resulting in a higher
ETc. Effective rainfall exceeds crop water demand during the initial and crop development stages of cotton, resulting in no irrigation operations during that period. Maximum irrigation was done at mid-season and late-season stages, with a total water depth of 346 mm constituting almost 90.7% of the total irrigation requirement.
Table 3 shows the irrigation scheduling for cotton, obtained at a fixed application depth of 100 mm and irrigation at critical depletion. The proper amount of water and timing of irrigation are determined by the irrigation schedule [
26,
27,
28]. The irrigation started 118 days after planting. The next watering is performed after 33 and 29 days, respectively, from the previous watering. The rooting depth of cotton at the initial growth stage is 0.3 m and is linearly increasing with increasing the growing period and will attain a constant value of 1.40 m. The critical depletion fraction is significant and at its maximum stage of harvesting. Therefore, drought stress occurs at the root zone, affecting crop production. Thus, a larger amount of water is needed to be supplied at a high flow rate of 0.5 L/sec/hectare. At this stage, the soil moisture deficit is also high. During the entire growing period of cotton, the root zone depletion is below the
RAW, resulting in no water stress. Root zone depletion was found to be comparatively higher at the mid-season and late-season stages due to less rainfall. Adjusted evapotranspiration almost satisfied crop evapotranspiration as there was no water stress in the soil moisture zone. Comparatively, less water was supplied at a flow rate of 0.14 L/sec/hectare during the early period of the mid-season stage due to a low water deficit. During the transmission of water from the source to the root zone, part of it is lost due to percolation, infiltration, surface runoff, and transmission losses. Therefore, irrigation efficiency must be calculated to account for the losses and to evaluate
NIR to increase the moisture content up to field capacity. The irrigation module is calculated by considering the highest
ETc value [
26].
CWR and irrigation scheduling were similarly analysed for other crops (
Tables S1–S8).
Sugarcane, being a perennial crop, requires water for the entire base period. Whereas sorghum requires a larger amount of water for cultivation as it is a rabi crop. Maize has a base period from May to September, and almost enough rain occurs during these months. Hence, frequent applications of water are not necessary.
For sunflowers, very little rainfall would occur. Hence, the cultivation of this crop involves a significant cost, as almost 91% of the total water requirement has to be satisfied by irrigation. Irrigation scheduling for sugarcane is obtained at refill soil at field capacity and irrigation at critical depletion. The irrigation operation is started after 124 days of planting, and there is no water stress problem throughout the entire crop period. Since, in this case, the soil is not saturated up to its field capacity, a significant amount of water needs to be supplied to avoid root zone depletion. Simultaneously, GIR is also high at each irrigation phase. Since this crop canopy expansion rate is very high, significant transpiration would take place, resulting in a high value of Kc throughout the entire growing period. Subsequently, ETc is also high.
Almost 94% of
ETc is fulfilled by irrigation for sorghum. The
Kc value is comparatively low as it is a medium-deep-rooted crop. It is cultivated at a 50 mm depth of water application and at critical depletion. The first watering is given to the crop 68 days after direct sowing. The entire
NIR is satisfied by three intervals of irrigation. At each phase of irrigation, a significant moisture deficit is observed, which ultimately affects crop productivity. Almost 54% of the moisture in the available water is depleted at harvesting for sorghum. Actual irrigation requirements and actual and potential water used by any crop can be estimated by the following equations (Equations (10)–(12)):
where
Ia = actual irrigation requirement (mm),
DMH = moisture deficit at harvest (mm),
AET = actual water used by the crop (mm),
Peff = effective rainfall (mm),
PET = potential water used by the crop (mm), and
ηi = efficiency irrigation scheduling (%).
All the parameters are simulated by CROPWAT 8.0. Equations (10)–(12) are valid only when irrigation is done at critical depletion and for a fixed application depth at field capacity.
Sunflower has a constant Kc value in the initial and mid-season stages. As it is a deep-rooted crop and it uses almost 77% of its base period to attain its maximum height, the canopy expansion rate is initially low, and at the mid-season stage, the transpiration rate is almost uniform. Thus, the Kc value remains the same throughout the initial and mid-season stages of crop development.
The following
NIR values of 200 mm, 1423.2 mm, 220 mm, 150 mm, and 150 mm for cotton, sugarcane, maize, sorghum, and sunflower, respectively, were obtained. The number of irrigation applications for cotton, sorghum, and sunflower was three times, for sugarcane it was five times, and for maize it was four times [
26]. The gross irrigation requirement was estimated by considering 70% application efficiency for the furrow irrigation system. It basically varied between 60 and 80% for surface irrigation systems [
22]. It is a user-defined parameter. Similar results were obtained for the remaining crops by using CROPWAT 8.0.
Figure 4 shows the soil water retention-growing length curve for the cotton.
Figures S1–S4 show the same for the other crops. Soil water retention basically indicates the depletion of soil moisture with reference to the Total Available Water (
TAW) and the existing water that the plants can extract for their growth. Growing length for a crop will give the duration of different growth stages from sowing to harvesting.
For cotton,
TAW is uniform after the crop development stage, resulting in constant crop height. However,
RAW is constant up to 55 days of crop growth, and after that, it will significantly fall as moisture depletion is very high.
Figure 4 shows the irrigation scheduling for cotton at a fixed application depth of 100 mm and irrigation at critical moisture depletion. From
Figure 4, it is clear that initially the moisture depletion level for cotton is low as a percentage of field capacity,
FC, as effective rainfall almost satisfies the water demand. It increased linearly after the crop development stage, and moisture depletion was at its maximum, i.e., almost 71.4% of total available moisture (
TAW), at 150 days from its showing due to less rainfall. The water deficit at the root zone increased by 7.14% within 32 days after the crop development stage, and it was again decreased by 3.6% within 30 days after the mid-season stage as the flowering stage had been reached at that time and the water requirement was comparatively low for attaining full maturity.
From
Figures S1–S4, it can be observed that the moisture depletion for sugarcane is not seen in the soil water retention—crop period curve after 124 days of planting. After 49, 38, 35, and 43 days of consecutive irrigation applications, no moisture deficit was found. Here,
TAW and
RAW are constants throughout the planting and harvesting of the crop. Initially, the moisture depletion did not occur until 65 days after planting. Therefore, it did not affect the crop growth, but at the time of harvesting, it was low, and the variation of moisture depletion showed a marginal difference throughout the initial to mid-season stage. For maize,
TAW was initially 55 mm and then increased linearly up to 180 mm, and after that, it was constant at 55 days after planting. The crop height for maize also varies between 0.3 m and a maximum of 1.0 m. Thus, a linear relationship (
Figure S5) was obtained between the
TAW and the crop height for maize.
where
h = crop height (m) and
TAW = total available water as soil water retention (mm)—as already remarked previously.
For sorghum, irrigation started at a 50% moisture depletion level. However, with the passage of time, significant irrigation helped to achieve
RAW as 80% of
TAW, and this condition occurred after 90 days from the crop planting (i.e., after the mid-season stage). Here also, a linear relationship was found (
Figure S6) between the crop height and
TAW:
For sunflowers, the maximum moisture depletion occurred 105 days after planting (i.e., 93% of
TAW). From
Figure S4,
RAW was 45.4% of
TAW, and it was then linearly varying up to the crop development stage. Then, it remained constant up to the end of the mid-season stage and was almost 50% of
TAW. Then
RAW achieved 78.7% of
TAW at the end of the mid-season stage and up to harvesting. The following best-fit curve between the moisture depletion percentage and
TAW was obtained for sunflower (
Figure S7):
where
MD = moisture depletion percentage. The above relationships (Equations (13)–(15)) were achieved based on the simulated output obtained from CROPWAT 8.0.
3.3. Crop Characteristics
Table 4 shows different conservative and non-conservative crop parameters for cotton.
Tables S9–S12 show the same for maize, sorghum, sugarcane, and sunflower, respectively, under no-water, salinity, and fertility stress.
The initial canopy cover for cotton was 32% of the maximum canopy cover, and the canopy senescence occurred after 32 days from when the maximum canopy had been reached. The time taken to reach maximum canopy was 62% of the total growing period, and senescence occurred 32 days after reaching maximum canopy cover as significant moisture depletion was seen at this phase. This indicates that significant transpiration can take place throughout the entire growing cycle. The temporal growth of the cotton canopy is given by Equation (16)—with a coefficient of determination of
R2 = 0.99—and shown in
Figure 5.
where
GCC = growth of canopy cover (%) and
B is the crop growth period in days.
GCC indicates the development of canopy structure, and crop growth period indicates the time taken by a crop to reach every stage of its development.
Similarly, the variation of the root zone development for cotton with the crop growth cycle (
Figure 6) can be expressed as follows (
R2 = 0.99):
where
Zr = root zone depth (m) and
B is the crop growth period in days, as previously stated. The developed equations (i.e., Equations (16) and (17)) from
Figure 5 and
Figure 6, respectively, are based on the simulated crop characteristics for cotton in AQUACROP 6.1, depending on the available climatic and soil parameters.
The development of the Harvest Index (
HI) with the growth period for cotton showed a gradually increasing slope with a minimum value at the start of the flowering stage and a constant value at the start of maturity (i.e.,
HI around 35%). Being a C3 type of crop, the slope of the curve between the biomass production and ∑
Tr/
ETo is less as compared to that of C4 crops. Here,
Tr is transpiration from the plant stomata. Thus, the overall yield will be low. For the given field condition, 38% of water stress was estimated during the stomatal closure and 14% at early senescence. It should be noted that we have neither used any sensors nor conducted any field experiments for the concerned crops for which water stress was calculated. The water stress coefficient was determined by using Equation (8), which is evaluated in CROPWAT 8.0 software. Moreover, we have used AQUACROP 6.1 software to calculate the yield of crops, and by default, this software would determine the water stress for different stages of crop development like canopy expansion, stomatal closure, early senescence, etc. [
20]. We have not directly measured this parameter. Basically, the water stress coefficient during early senescence (
Ks,sen) is responsible for the reduction of green canopy cover, whereas the same for stomatal closure (
Ks,sto) plays a major role in reducing crop transpiration and root zone expansion. The target model parameters for calculating
Ks,sen and
Ks,sto, are green canopy cover (
CC), crop transpiration coefficient (
KTR), harvest index (
HI), and depth of root zone (
Zr).
Therefore, the transpiration process was affected, resulting in a lower yield for the crop. The
HI is required to be adjusted by a factor that takes into account inadequate photosynthesis, pollination stress, and water stress [
8]. For cotton,
HI will be 2% higher in the vegetative period and 20% higher during yield formation to take into account the above-mentioned stresses. A significant amount of biomass can actually be produced as compared to potential biomass. AQUACROP 6.1 will give the value of biomass for each crop depending on crop characteristics, soil condition, cultivation operation, and climatic conditions. The biomass is calculated by using Equation (18) [
8]:
where
B is the biomass produced (tonne/ha),
WP* is the normalised water productivity for a reference CO
2 concentration of 369.41 ppm (g/m
2),
Tr is the daily transpiration from plant stomata (mm), and
ET0 is the daily reference crop evapotranspiration (mm). The normalised water productivity is the slope of the curve between
B and
.
For sugarcane, the initial canopy cover is very high, resulting in a higher leaf area index (
LAI). The maximum canopy cover was developed at the crop development stage, and the crop senescence occurred 266 days after the maximum canopy cover developed. The growth of canopy cover for sugarcane is related to the crop growth period according to Equation (19), as shown in
Figure S8:
where
GCC and
B are expressed in percentage and days, respectively.
The root zone expansion rate is very high for sugarcane, so the maximum crop height can be achieved within the minimum possible time after transplantation. Sugarcane is a C4-type crop. Hence, water productivity is also high. No water stress was seen during canopy expansion, stomatal closure, and early senescence, resulting in high crop production; thus, actual biomass production tends to be the same as the potential biomass. The root zone development for sugarcane shows a linear relationship with the growing period according to Equation (20), as shown in
Figure S9:
For maize, initial canopy cover is very high, and plant density is also very high, resulting in a higher
LAI. The maximum canopy cover occurred at the end of the development stage. The duration of flowering was only 13 days, so the production of dry yield declined. The root zone development for maize shows a linear relationship with the growing period according to Equation (21) as shown in
Figure S10:
The water productivity for maize was high and within the range of 30–35 gm/m2. The water stress was small at stomatal closure; hence, the transpiration process was not affected at the crop development stage; the water stress was comparatively high during canopy expansion. An adjusting factor of 1.01 was applied to the reference harvest index (HIo). No adjustment was required during the vegetative period, but only an adjustment of 1% of HIo was needed during yield formation to calculate adjusted HI. Here, the adjustment is positive. In a similar way, relationships between growth of canopy cover (GCC) versus growing days and root development versus growing period length were achieved for the other crops (e.g., sorghum and sunflower).
The above Equations (16)–(21) were achieved based on the simulated output obtained from AQUACROP 6.1.
3.4. Soil Surface Characteristics
Table 5,
Table 6 and
Table 7 show different soil surface characteristics for black clay, medium (loam), and red loamy soil, respectively, as derived from CROPWAT 8.0.
Black clay soil and medium (loam) soil exhibited the highest values of TAW and the maximum rain infiltration rate. Black clay soil had an initial moisture depletion level of 50% as a percentage of TAW. Hence, cotton and sugarcane faced significant water stress problems, which could be overcome by frequent irrigation to increase the moisture content up to FC. However, medium and red loamy soil initially had no moisture depletion. So, the entire available water could be utilised by other crops.
AQUACROP 6.1 deals with (and determines) more parameters than CROPWAT 8.0, which could be more convenient for further research.
Table 8,
Table 9 and
Table 10 show the different soil surface characteristics for the above-mentioned soils from AQUACROP 6.1.
AQUACROP 6.1 underestimates TAW for black clay soil, whereas it overestimates the same for medium and red loamy soil. Hence, it could be advantageous to adopt CROPWAT 8.0 for black clay soil and AQUACROP 6.1 for medium and red loamy soil in this regard. The saturated hydraulic conductivity for black clay soil was almost equal to the maximum infiltration rate. For black clay soil, FC almost met the saturation limit; hence, readily evaporable water was very high, resulting in a higher ETC.
From
Table 9 and
Table 10, a relationship can be observed between drainage characteristics (ζ) and saturated hydraulic conductivity, expressed by Equations (22) and (23)—as shown in
Figure 7 and
Figure 8—for medium and red loamy soil, respectively. The variation is basically a power law with a high value of
R2. The developed equations (i.e., Equations (22) and (23)) from
Figure 7 and
Figure 8, respectively, are based on the simulated soil characteristics for medium and red loamy soil in AQUACROP 6.1, depending on the crop phenology.
where
ζ expresses the fraction of drainable water day by day. For both equations,
R2 = 0.99. AQUACROP 6.1 also determines the Curve Number (
CN) and the hydrological soil group for black clay, medium (loam), and red loamy soils. They are D, A, and C, respectively. It can be observed that for medium soil, the infiltration rate was higher as compared to the other two soils, resulting in higher percolation losses. Hence, water availability at the field could be lower, affecting the duty of the cultivated crops; thus, less area could be irrigated, though significant water supply could be carried out.
3.5. Furrow Irrigation Design Parameters
Table 11 shows different furrow irrigation design parameters (e.g., application efficiency, storage efficiency, cut-off time, advance time, recession time, furrow length, surface runoff ratio, and deep percolation ratio) by taking into account different input parameters (e.g., flow rate, furrow cross-section) for cotton. The FURDEV model calculates the values of uniformity coefficient (
UC), distribution uniformity (
DU), deep percolation ratio, and surface runoff ratio for furrow irrigation based on input parameters like flow rate, furrow cross sections, and soil infiltration characteristics for the crop under consideration. These above-mentioned parameters have been determined by using Equations (24)–(27) [
19,
24].
where
Di is the depth of water in the furrow for
ith emitter (mm);
Davg is the average infiltrated depth (mm);
n is the number of emitters along the length of the furrow;
Dmin is the minimum infiltrated depth (mm);
Ddp is the deep percolation depth (mm); and
Dsr is the surface runoff depth (mm), which is the difference between the actual infiltrated depth (
Da) and the average infiltrated depth (mm) (i.e.,
Dsr =
Da −
Davg).
Tables S13–S15 show the same for maize, sorghum, and sunflower.
The programme FURDEV clearly shows the impact of different furrow irrigation design parameters on the yield of crops as well as the proper management of irrigation water with the help of AQUACROP 6.1. Among all the three forms of operation modes, the cutback flow method and tailwater reuse method were found to be more convenient because the application efficiency and storage efficiency were relatively higher as compared to the fixed flow method for the same flow rate, soil infiltration characteristics, and furrow cross sections [
19]. However, the under-irrigation depth and over-irrigation depth were almost the same in these three modes.
From
Table 11, it was observed that cotton has a less advanced ratio, requiring more uniformity in irrigation. A significant amount of water can be utilised from the tail end of the field as the recovery ratio is relatively high, which ensures an economical and controlled use of water. Cotton has a very short depletion phase and a long ponding phase, indicating that the rate of decrease in surface water storage was low and the availability of irrigation water in the field was high for a long period of time. Due to the long ponding phase, the depth of water in the field increased, resulting in a lower moisture deficit. Cutback flow mode was more convenient for cotton as surface runoff depth was low. More water can be used to increase the water content in the soil moisture zone. Application efficiency increases by 16.14% and 13.3% for cutback flow and tailwater reuse methods compared to fixed flow methods, respectively. Whereas distribution uniformity and uniformity coefficient were independent of furrow length and flow rate for cotton. Over irrigation depth was slightly higher as compared to under irrigation depth for cotton, maize, and sunflower, except for sorghum. Additionally, it was higher when the fixed flow method was used. Therefore, the other two modes of operation could be beneficial for obtaining the maximum possible efficiency. More controlled water application may be possible for maize with the lowest cutback ratio. However, for maize, surface runoff depths were comparatively higher, and they were higher for the tailwater reuse method for all the crops. Whereas the cutback flow method simulates a very small amount of surface runoff for all the crops with the same input parameters. A reverse curve was obtained between the relationship between flow rate and application efficiency for the crops for cutback flow and tailwater reuse methods.
From
Table 11 and
Tables S9–S12, one can observe: (i) spacing of furrows, gradient, surface roughness, and furrow cross sections had no influence on the furrow length, flow rate, cutoff time, application efficiency, storage efficiency, and many other design parameters; rather, they depended upon crop type; (ii) a larger furrow length would result in a larger distribution uniformity and uniformity coefficient. Hence, more uniformity in irrigation could be obtained. In this regard, sunflower had the minimum Distribution Uniformity (
DU) and Uniformity Coefficient (
UC) and maize had the maximum
UC and
DU; and (iii) surface runoff ratio increased with increasing furrow length except for sunflower. A lower flow rate will give a higher deep percolation ratio for all the crops as the applied water depth on the field will be less, resulting in the development of significant lengths of under irrigation, which clearly signifies the occurrence of water stress and lower application efficiency.
Application efficiency increases with increasing the furrow length for sunflower, cotton, and sorghum, but storage efficiency is almost independent of furrow length. Lower furrow length gives lower application efficiency, resulting in a more non-uniform water application. Additionally, improper water distribution does not satisfy the crop’s water demand. The distribution uniformity (DU) and uniformity coefficient (UC) are independent from the flow rate, furrow length, and operation mode for all the crops. The average infiltration opportunity time for cotton, maize, sunflower, and sorghum is calculated as 1513 min, 1093 min, 794 min, and 1027 min, respectively. So, more non-uniformity in irrigation operations has been seen for cotton, which has a large opportunity time.
Additionally, by using the results of FURDEV, the advance curve and the irrigation profile were prepared.
Figure 9,
Figure 10 and
Figure 11 show the advance curve and the irrigation profile for the cotton for three different modes of operation. In a similar way, the same profiles were obtained for the remaining crops and shown in
Figures S11–S19 under the same modes of operation. The advance curve shows the variation of different hydraulic phases of furrow irrigation with furrow length and gives the value of infiltration opportunity time, which is the difference between advance time and recession time. It will signify uniformity in irrigation operations, whereas the irrigation profile shows different water depths in the field along the furrow length, which further indicates whether there is any possibility of waterlogging due to over-irrigation, surface runoff loss, or the occurrence of water stress [
19].
From the advance curve, it can be observed that cotton has the highest opportunity time while sunflower has the lowest one. Maize and sorghum have almost the same opportunity time, indicating that more uniformity in the irrigation operation will ensure more infiltration depth through the soil. From the advance curve, it can be observed that the advance time gradually increases with increasing the furrow length for all the crops under all the modes of operation. In the case of cotton, a relationship (with
R2 = 0.997) can be established between the advance time and the furrow length as follows (
Figure 12). The relationship was established based on the simulated output obtained from FURDEV and the curve of best fit.
where
Ta is the advance time (min) and
L is the length of the furrow (m).
Analogously, the advance curve for sorghum (
Figure S20) can be expressed by Equation (29). This relationship was established based on the simulated output obtained from FURDEV and the curve of best fit:
Similar relationships can be obtained for the remaining crops. More advance time will ensure that the water takes a long time to reach the downstream end of the field until the inflow has been started; therefore, more infiltration loss and deep percolation loss will occur. Subsequently, the application’s efficiency will decrease. Whereas the furrow length has no influence on the recession curve, 83.4% of the total furrow length will consume maximum over-irrigation length for maize and 58.42% of the total furrow length will consume minimum over-irrigation length for sorghum. This clearly indicates that a significant portion of over-irrigation length would occur, resulting in deeper percolation loss and surface runoff loss for the crops.
Figure 13 shows that a polynomial relationship can be obtained between the under-irrigation depth and the furrow length for cotton (Equation (30)). The relationship was established based on the simulated output obtained from FURDEV and the curve of best fit (with
R2 = 0.993):
where
Du = under-irrigation depth (mm) and
L = length of the furrow (m). Similar relationships were achieved for the other crops.
Figure 14 shows the crop productivity depending on the irrigation water for cotton. Similarly,
Figures S21 and S22 show the same for maize and sunflower, respectively.
The crop productivity-irrigation water relationship (
Figure 14,
Figures S21 and S22) for cotton, maize, and sunflower can be established from the output obtained from AQUACROP 6.1.
Initially, for cotton, the water productivity (WP) was small, and it increased with decreasing water application. Similarly, the dry yield for cotton decreased from 2.6 tonne/ha to 1.68 tonne/ha at a decreasing water depth from 1000 mm to 200 mm. It was noted that there was a significant reduction in the irrigation water requirement of 80%, with a reduction in the ET water productivity of 10.6%.
Similarly, for maize, the optimum depth of water application was 275 mm and the maximum
WP was 2.4 kg/m
3; dry yield linearly increased with increasing the irrigation water. The highest yield was obtained at a water depth of 770 mm.
Figure S21 shows a reduction in the irrigation water requirement of 85.7% with a reduction in the
ET water productivity of 7.76%. Similarly,
Figure S22 shows a reduction in the irrigation water requirement of 50% with a reduction in the
ET water productivity of 1.45% for sunflower. These results clearly highlight the need for controlled water management for irrigation to optimize the beneficial use of water.