Journal of Irrigation and Drainage Engineering

Significant improvements in the profitability and sustainability of irrigated areas can be obtained by the application of new technologies. In this work, a model for the simulation of water flows in irrigation districts is presented. The model is based on the combination of a number of modules specialized on surface irrigation, open channel distribution networks, crop growth modeling, irrigation decision making, and hydrosaline balances. These modules are executed in parallel, and are connected by a series of variables. The surface irrigation module is based on a numerical hydrodynamic routine solving the Saint Venant equations, including the heterogeneity of soil physical properties. The simulation of water conveyance is performed on the basis of the capacity of the elements of the conveyance network. Crop growth is simulated using a scheme derived from the well-known model CropWat. The irrigation decision making module satisfies water orders considering water stress, yield sensitivity to stress, multiple water sources, and the network capacity. Finally, the hydrosaline module is based on a steady state approach, and provides estimations of the volume and salinity of the irrigation return flows for the whole irrigation season. The application of the model to district irrigation management and modernization studies may be limited by the volume of data required. In a companion paper, the model is calibrated, validated, and applied to a real irrigation district.

In a companion paper a model for the simulation of water flows in irrigation districts was formulated. The model combines a series of modules specialized in surface irrigation, open channel distribution networks, crop growth modeling, irrigation decision making, and hydrosaline balance. The objective of this paper is to calibrate, validate, and apply the model, using the Irrigation District Five of Bardenas (Spain) as a study area. Two years of study were used for the analysis, which could be classified as normal (2000) and dry (2001) from the point of view of crop water requirements. Model calibration was performed in one of the 11 hydrological sectors in which the district is divided. The control variable was the monthly water demand, while the calibration variables were related to irrigation operation and scheduling. The seasonal differences in observed and simulated water demand amounted to 0.9 and 1.9% for 2000 and 2001, respectively. Model validation was performed in the rest of the sectors, and the regression line of observed versus simulated monthly water demand could not be distinguished from a 1:1 line in both years. Model application explored scenarios based on management improvement (controlling the irrigation time) and structural improvement (increasing drainage water reuse for irrigation). These scenarios permitted one to sharply reduce water demand, halve the irrigation return flows, and reduce the daily irrigation period from $24\text{to}16\mathrm{h}$ .

Optimum land and water allocation to different crops grown in different regions of an irrigation scheme is a complex process, especially when these irrigation schemes are characterized by different soils and environment and by a large network of canals. At the same time if the water supply in the irrigation schemes is limited, there is a need to allocate water both efficiently and equitably. This paper describes the approach to include both productivity (efficiency) and equity in the allocation process and to develop the allocation plans for optimum productivity and/or maximum equity for such irrigation schemes. The approach presented in this paper considers the different dimensions of equity such as water distribution over the season, water distribution during each irrigation, and benefits generated. It also includes distribution and conveyance losses while allocating water equitably to different allocation units. This paper explains the approach with the help of the area and water allocation model which uses the simulation–optimization technique for optimum allocation of land and water resources to different crops grown in different allocation units of the irrigation scheme.

The area and water allocation model which uses simulation–optimization technique for optimum allocation of land and water resources to different crops cultivated in different allocation units of the irrigation scheme was modified to include both productivity and equity in the process of developing the allocation plans for optimum productivity and/or maximum equity. This paper illustrates the potential of this approach with the help of a case study on Nazare medium irrigation scheme in India. The allocation plans were developed for optimization of different performance parameters (productivity and equity) for different management strategies based on irrigation amount and irrigation interval and cropping distribution strategies of free and fixed cropping. The results indicated that the two performance objectives productivity and equity conflict with each other and in this case, equitable water distribution may be preferred over free water distribution at the cost of a small loss in productivity. Though these results relate to one case study, they show the value of the approach of incorporating productivity and equity in the allocation process with the help of the simulation-optimization model described in the companion paper.

Some irrigation schemes with limited water supply in Central and Southern India follow the area proportionate water distribution based on assumed uniform characteristics of the command area (planned schedule). However in most cases, this planned schedule is overridden by the practice in which users at head draw more than their share of water (actual schedule) due to human factors and technical limitations of the planned schedule. This practice is highly inequitable as users at tail end do not get any water. This paper considers alternative schedules based on full irrigation or deficit irrigation within the framework of area proportionate water distribution in such irrigation schemes and presents the simulation–optimization technique to develop the corresponding land area and water allocation plan for different allocation units by considering the heterogeneity of the irrigation scheme. This paper further demonstrates the utility of proposed alternative schedules by comparing the productivity and equity of these schedules with planned and actual schedules for one irrigation scheme in Central India. The results show that the actual schedule reduces both productivity and equity greatly and the productivity and equity with the alternative schedules are higher than with the planned schedule. The results also show that deficit irrigation has great potential to increase both productivity and equity of irrigation schemes.

The paper proposes a new method to tune robust distant downstream proportional integral (PI) controllers for an irrigation canal pool. The method emphasizes the role of gain and phase margins in the controller design, by linking the selection of these robustness indicators to the time domain specifications. This leads to link the frequency domain approach used by automatic control engineers to the time domain approach used by hydraulic engineers. The maximum error corresponding to an unpredicted perturbation is shown to be directly linked to the gain margin and the settling time to the phase margin of the controlled system. The tuning method gives analytical expressions for the controller parameters as function of physical parameters of the canal pool in order to satisfy desired performance requirements. The model is first expressed in terms of dimensionless variables, in order to get generic tuning formulas. The dimensionless PI coefficients are then expressed as functions of time-domain performance requirements. The PI tuning method is evaluated by simulation on a full nonlinear model for a canal pool taken from the ASCE test cases.

Based on the companion paper results, this paper presents a robust tuning method to tune a proportional integral (PI) controller that fulfills the design requirements for a single pool with different hydraulic conditions. The robust PI controller parameters are obtained analytically as a function of the physical parameters of the canal pool. Important implementation issues are also considered. Rules are provided for the sampling time selection in order to recover the continuous-time performance. When the sampling time is imposed, it has to be included in the controller design, by modifying the delay of the system. A simple way is proposed to take account of the gate opening as control action variable, instead of the upstream discharge. This robust PI tuning method is evaluated by simulation on a full nonlinear model of two different canal pools for different flow conditions and different implementations: continuous-time control, discrete-time control, using the discharge or the gate opening as control action variable. Simulation results show that the method leads to efficient realistic PI controllers for a canal pool.

Physically based modeling of the coupled one-dimensional surface and two-dimensional (2D) subsurface flow during furrow irrigation advance often causes numerical instabilities and nonconvergence problems. This is particularly the case for low irrigation advance rates when infiltration consumes a predominant part of the inflow volume. The proposed furrow advance phase model (FAPS) further develops the concepts of a previous study. An analytical zero-inertia surface flow model is iteratively coupled with the 2D subsurface water transport model HYDRUS-2. In contrast to the previous study, the flow domain is discretized using fixed space increments and the resulting set of nonlinear flow equations is solved using the Newton method. The complexity of the model was reduced by process adequate simplifications. FAPS exhibited better convergence, numerical stability, and less computational time than the original fixed time interval solution. The new solution converged rapidly for a number of model tests with various inflow rates including runs with very slow irrigation advance. Simulation model predictions agree very well with advance times measured in laboratory and field tests.

A one-dimensional water and solute transport UNSATCHEM model is calibrated and validated with a saline water use experiment for wheat and cotton crops. The model is further employed for regional scale salinity modeling with distributed data on soil, irrigation water supply, and its quality from six representative locations from the Kheri command of the Bhakra irrigation system. The wheat–cotton crop rotation, the main rotation in the command, is considered during long-term simulations. The CROPWAT model is used to determine the evapotranspiration requirements of different wheat and cotton crops, while soil water retention parameters are estimated by the RETC model. Atmospheric water and solute boundary conditions are assumed at the top boundary, while free drainage is considered for the lower boundary, as the watertable in the command is sufficiently deep. Simulated salinity and yield values are compared with observed values for regional validation of the model. Critical areas in the command are identified using regional scale modeling results, and applying irrigation water availability and root zone salinity criteria. Guidelines for sustainable conjunctive water use planning are for the Kheri command to get optimum agricultural production despite the use of saline water for irrigation under prevailing scenarios of water availability and its quality.

Irrigation application of effluent from lagoon based systems is an important aspect of animal manure management in swine rearing facilities. This paper describes a procedure for obtaining least-cost components and operational parameters for an irrigation system used to apply liquid manure. The paper employs propositions from hydraulics—affinity laws and Hazen–Williams formula—to derive statistically estimable relationships between pressure, liquid output, power requirement, and impeller speed of an irrigation system. Coefficient estimates are then used in a spreadsheet based decision support system that employs an optimization program to select the least costly irrigation system. A practical application is showcased through simulation of three hypothetical swine farms. Main findings are that the use of this optimization procedure can result in cost savings for swine growers applying manure by irrigation and that the cost savings potentials increase with the farm size.

Analytical expressions for rate and volume of flow depletion of semipervious streams due to sinusoidal variation in pumping rate are obtained. An analytical but approximate method is developed for obtaining the rate and volume of stream flow depletion due to arbitrary unsteady pumping discharge. The method uses the ramp kernel and convolution. The use of ramp kernels permits linear interpolation between two consecutive discretized discharge values. The analytical equations for the ramp kernels for the rate and volume of stream flow depletion are derived. The proposed method is applicable for homogeneous and isotropic aquifers that are hydraulically connected to streams.

This paper presents a theoretical derivation of a dimensionless unit hydrograph using the kinematic wave approach. The dimensionless kinematic wave unit hydrograph (KWUH) is normalized by the parameters associated with the equilibrium condition and can be generated according to the selected rainfall duration. For comparison, a KWUH whose rainfall duration is equal to the time of equilibrium produces good agreement with the Soil Conservation Service dimensionless unit graph (SCSUH). Both the KWUH and SCSUH are kinematic in nature and have a 37.5% runoff volume that occurred before the time to peak. To verify the applicability of the KWUH derived in this study, a set of hypothetical rectangular watersheds ranging from 20 to 200 acres was investigated. The predicted 100-year peak flow rates from these five hypothetical rectangular watersheds are compared with the SCSUH method and the Colorado Unit Hydrograph Procedures (CUHP). In general, the KWUH produces good agreement with the CUHP when the catchment has an area between 60 and 120 acres. Based on the study of hypothetical rectangular watersheds, it is recommended that the KWUH be applicable to urban catchments with a drainage area of up to 120 acres.

Knowledge of the sensitivity of various soil hydraulic properties is beneficial for model development and application purposes. It can lead to better estimated values, better understanding, and thus reduced uncertainty. In the present study, an extensive sensitivity analysis was performed to investigate the effects that various soil hydraulic properties have on subsurface water flow below furrows during two successive irrigation events to see which irrigation event was more sensitive and to analyze the effect of spatial variations in the initial soil water contents within the soil profile. Testing the sensitivity of the various soil hydraulic parameters in the van Genuchten-Mualem expression was carried out using the HYDRUS-2D model for two irrigation events $10\text{days}$ apart. Results showed that the first irrigation event was clearly more sensitive than the second one. The latter event was mainly associated with the nonuniformity of the initial soil water contents within the soil profile. Pressure heads in the soil profile were more sensitive than cumulative outlet fluxes and soil water contents. Sensitivity analysis results for pressure heads, cumulative fluxes, and water contents indicated that in every case the most sensitive parameter was the hydraulic property shape factor $\left(n\right)$ followed by the saturated water content $\left({\theta }_s\right)$ , the saturated hydraulic conductivity $\left(K_s\right)$ , the residual water content $\left({\theta }_r\right)$ , and the shape factor in the soil water retention curve $\left(\alpha \right)$ , with the pore-connectivity parameter $\left(l\right)$ the least sensitive parameter during both irrigation events. Pressure head sensitivity analysis for all parameters studied showed that the least sensitivity was linked with the wetting front as it gradually moved deeper with time, and the highest sensitivity was observed in those regions where the initial soil water contents were lower. Similarly, for water contents, higher sensitivity occurred in the drier regions during the first irrigation event and near the moisture front in the second irrigation event. Both pressure heads and water contents showed some sensitivity near the soil surface during both irrigation events, suggesting the importance of evaporation from the soil surface.

In the past, solutions to open channel flow problems involving free surfaces were generally found on the basis of experimental data or through the development of theoretical expressions using simplified assumptions. The volume of fluid (VOF) turbulence model is applied to obtain characteristics of three-dimensional open channel flows involving free surfaces. In particular, the VOF model is used to determine the pressure head distributions, velocity distributions, and water surface profiles for the free overfall in a trapezoidal open channel. The predictions of the proposed model are validated using existing experimental data.