Integrated Water Resource Assessment of Irrigation System of Haryana

Agriculture is the major activity in the state of Haryana and large volume of water is required to meet the irrigation demands of the crops grown. But, there is limited water availability in the state. Haryana receives water from Yamuna River and Bhakra system. Sowmelt, rainfall and groundwater are main sources of water in the catchment. It is essential to integrate the manmade canal system with hydrological system. This paper focuses on integrated hydrological modeling framework to conceptualize the system and to assess the Water Resources of the state. Snowmelt and Rainfall runoff modeling using GR4JSG model were combined to model the inflows to the irrigation system of Haryana. Irrigator canal model of eWater Source has been used to generate water demands from crops grown. The water balance and water use efficiency have been worked out for each district of Haryana. The hydro climate input data, stream flows, crop data and soil data have been used in the study. The flows modeled at Tuini (P), Yashwant Nagar, Bausan, Haripur, Poanta and HKB sites were compared with the observed flows. The objective function of NSE Daily and log Flow duration was used for model calibration and validation at various locations up to Mathura, the outlet of the study area. The value of the objective function at Mathura was 0.54, a fairly good value. The results of the Irrigator canal model have shown that all the Inflows, Outflows and the Utilizations of water have been properly balanced for each district. The water use efficiency of districts varies from 27% to 59%. The overall water use efficiency for Haryana canal system has been calculated as 39%. This is low value indicating excess water is being extracted to meet the water demands.


Introduction
Haryana state uses majority of water resources in agriculture sector.The state receives surface water from Yamuna, Sutlej, Ravi and Beas rivers as per various interstate water sharing agreements.The state does not have any perennial river.
The Western Jamuna Canal (WJC) system and Bhakra system are the two main canal systems irrigating 2.97 million hectare area.The intensity of canal irrigation is not uniform throughout the state due to highly skewed distribution of canal water in different commands.Some areas are getting more water and the others getting less because of capacity and other constraints.The water use efficiency in the state is 60% [1] and for WJC system the water use efficiency is 36% to 43% [2].The present water availability in the state is 14 Million Acre Feet (MAF) or 17 Billion Cubic Meter (BCM) and projected water demands is 36 MAF (45 BCM) in 2045 [1].There is a huge gap between demand and supply of water resources.Ground Water (GW) has been over exploited that have caused depletion of water table and deterioration of quality of water.The ground water has declined in parts of the state and also increased in other parts of the WJC command [3].The sweet water is going down and saline water is rising up.The per capita freshwater availability is declining and there is rapid deterioration of the available resources in the state and low flows in Yamuna causing environmental issues.
Snowmelt is a major contributor to the Hydrology of Yamuna and other Himalayan Rivers.Snow melt contributes about 30% -50% of annual stream flow in Himalayan Rivers [4].The flows in the Yamuna River are higher from June to September due to monsoon rainfall and snowmelts in the upper catchment.But during non monsoon the water in river is due to snowmelt and base flow.The variation in rainfall has visible effect on the flow.There is temporal as well as spatial variability in the amount of water supplied due to lack of storage across Yamuna.At HKB, the average high flow is 7079 cumecs and minimum flow is 70 cumecs [2].
Snow melt water, surface water of canals/streams/drainage system, ground water and crop water demands are interconnected and play an important role in the study area.It is essential to integrate the manmade canal system with hydrological system.Hence, there is an urgent need of an integrated hydrological modeling framework to conceptualize the system and to assess the Water Resources of the state.The system has been conceived to have the catchment models and schematic canal model.The hydrological modeling (snow and rainfall runoff) of the river basin is the usual water resource assessment methodology by means of considering basin as a natural unit and using land uses and hydrological response units.The methodology used in this study has been depicted in Figure 1.Rainfall runoff model with snow melt from hilly areas are conceptualized and the water flows are diverted to canal system of Haryana through WJC system.The flows are utilized to meet the crop water demand generated by Irrigator canal model.

Study Area
The present study covers the Yamuna River Basin up to Mathura and adjacent Irrigation Canal System of Haryana (Figure 2).The outlet of the study area is at Mathura on Yamuna.Yamuna originates from the Yamunotri Glacier near Bander Punch peak (38˚59'N, 78˚27'E) in the Mussourie range of lower Himalayas, at an elevation of 6320 meter above mean sea level (amsl) in Uttarkashi District of Uttarakhand.The study is focused on an area of 60,557 km 2 including 11,397.78km 2 hilly area upstream of HKB and 27,740 km 2 irrigation command area of the state.The elevation ranges from 150 m amsl to 6253 m amsl in the study area.The river water at HKB is diverted into Western Jamuna Canal (WJC) and Eastern Jamuna Canal (EJC) to supply water for irrigation and drinking purposes to the States of Haryana and Uttar Pradesh respectively and also drinking water supply of Delhi.
Yamuna basin is a heterogeneous and has high spatial variation of rain fall varying from 200 mm to 2350 mm.The annual average rainfall is 906 mm.The rainfall in the study area is from South West monsoon and major part of rain fall is received between June to September.The rainfall is scanty in winter.However, temperature and rainfall vary with altitude.The mean maximum and minimum temperature vary between 24˚C to 45˚C and −1˚C to 11˚C respectively.The annual rainfall is much lower i.e. about 1000 mm in the foot hills to less than 600 mm in Haryana and Delhi.Haryana is mostly arid or semiarid.There is limited rainfall in Haryana ranging from 300 mm in the South West to 1100 mm in the North East.Average annual rainfall is 573 mm.The soils in the Yamuna basin are predominately alluvial.

Methods and Materials
This section describes the GR4JSG model for snow melt and rainfall runoff, the Irrigator Demand model for generation of crop water demand and the input data used in the study.

Description of GR4JSG Model
GR4JSG model developed by [5] has been used for rainfall runoff, snowmelt and glacial melt simulation.GR4JSG model is based on GR4J model [6].The degree day factor has been use to describe the snow and ice melt processes [7].The GR4JSG model is used within the eWater Source modeling framework [8].The conceptual structure of the model is shown in Figure 3, with the original GR4J model developed by [6] and it is shown on the left side and the additional snow and ice stores shown on the right side.The GR4J model is a daily lumped rainfall-runoff model with four parameters (x 1 , x 2 , x 3 and x 4 ).GR4J consists of two main stores: the production store and the routing store as shown in the conceptual structure of Figure 3.The time series inputs to the model are rainfall depth (P) in mm/day and potential evapotranspiration (E) in mm/day.The parameter x 1 controls the size of the production store (mm), x 2 controls the flux to groundwater (mm/day), x 3 controls the size of the routing store (mm) and x 4 controls the recession of the unit hydrograph (days).
The snow and glacier processes are represented in GR4JSG model by two conceptual stores for snow and ice.The temperature is used to decide if precipitation is rainfall or snow.The snow accumulation and melt are represented by snow store and glacier melt processes are represented by the ice store.Snow is Figure 3. Conceptual structure of the GR4JSG (GR4J + Snow and Glacier) model [5].melted prior to the glacier melting and it is accumulated in the snow store.The snowmelt obtained from the snow area infiltrates the soil and enters the production store.Afterwards it is treated as though it was rainfall.The production store gets contributions from snowmelt and rainfall.The amount of snowmelt contribution to the runoff is calculated separately.For more details on GR4J and GR4JSG models please refer to [5] and [6].

Irrigator Demand Model (Crop Model)
A supply point and water use nodes have been used to model the water use in Source platform [9].The demand models are provided by the water user nodes and that can be configured to represent irrigation demand.Irrigator is applied through a water user node and connected to at least one supply point node to provide water to meet with the irrigation demand.Irrigator maintains a daily water balance for each cropping area during its planting season so as to calculate the daily soil water deficit and irrigation requirement.A water balance of the fallow crop is also used to initialize the soil moisture store when new crops are grown.The fallow cropping area does not order or receive irrigation water.The daily water balance includes effective rainfall, irrigation, runoff, evapotranspiration and deep percolation that affect changes in soil water storage.The irrigator uses the method described in FAO 56 [10] to represent the daily crop water balance.Irrigator is conceptualized to represent at irrigation district scale rather than an individual irrigator.The total requirement for the district is calculated by summing each of the individual cropping area requirements.The district requirement is then adjusted to allow for any district supply escapes and losses associated with delivering water to the district.For more details refer to FAO 56 [10] [11] [12].

Basin Area and Elevations
Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) has been used for delineation of the study area watershed [http://srtm.csi.cgiar.org/].The horizontal resolution is 90 × 90 m of used in this study.The projections are WGS 1984, 43N.ArcMap 10.2 is used to extract sub catchments and to calculate the area and the area-elevation distribution of the basin.The catchment has been divided into 13 sub catchments, including 7 sub catchments in hilly areas.The elevation of the study area varies from 150 m to 6253 m.The basin is subdivided in 11 elevation zones/bands each with a vertical difference of 600 m in hilly areas (Table 1).Each elevation zone acts as a functional unit.Source conceptualizes catchment as a combination of functional units [13].A functional unit (FU) is a part of catchment that functions in a similar way for hydrological processes, and which may be defined by different climate, land use, soil type or runoff processes.
Functional unit may have different model parameters or inputs to another FU in the catchment.For this study, temperature, precipitation and evapotranspiration were used for each FU.The runoff at all 11 FUs from GR4JSG models was summed to calculate runoff at a point.The areas enclosed by various elevation contours can be determined with Arc Map 10.2 spatial analyst tool by making use of the zone boundaries and other selected contours in the basin.The mean hypsometric elevation of each zone is used as the elevation to which base station temperatures are extrapolated for the calculation of number of zonal degree-days (Table 1).

Hydro Meteorological Input data
The time series data of daily precipitation, daily minimum and maximum temperatures, potential evapotranspiration (PET), stream flow and monthly snow cover has been used to calibrate and validate the model.The data of Indian Meteorology department (IMD) and Central Water Commission (CWC) has been used in this study.

1) Temperature
The IMD gridded maximum and minimum temperatures data from year 1995 to 2013 has been used in this study.As temperature varies with elevation, so, temperature lapse rate has been used to estimate the temperature at the unmeasured locations.The temperature lapse rate has been taken from the literature of nearby basins of Himalayan region due to lack of historical data of the temperature lapse rate in the Yamuna basin.The annual lapse rates 0.51˚C/100 m to 0.68˚C /100 m have been found by [14] for Snowmelt Runoff Model to Tamor catchment.The seasonal lapse rates of 0.55˚C/100 m for monsoon (June to September) and 0.6˚C/100 m for the dry season (October to May) has been derived by [15] in Dudh Koshi catchment.The lapse rate from 0.65˚C and 0.69˚C per 100 m were found by [16] in Beas basin during the years 1998 and 1999 respectively.The temperature lapse rate 0.65˚C/100 m was used by [17] in Alakhnanda and Bhagirathi Basin of Ganga basin.A temperature lapse rate of -0.65˚C per 100 m has been used in the study.

2) Precipitation
The IMD gridded precipitation data from year 1995 to 2013 has been used in this study.It is difficult to extrapolate precipitation in mountainous watersheds due to lack of weather stations in the watershed and local factors.The topography strongly influences the spatial distribution.However; the rainfall scaling factor of 10 mm/km has been used from literature.

3) Potential evapotranspiration
The global data from the http://hydrology.princeton.edu/data.measures_pet.php from year 1995-2013 has been used.The PET is corrected using PET decay rate of −0.0005 mm/1000 m for Himalayan area [18].

4) Discharge Data
The The calibration was done using the "Shuffled Complex Evolution then "Rosenbrock" option [11].The objective function used is NSE Daily and log flow durations.NSE is Nash Sutcliffe Efficiency [19].
The model has been calibrated at six different CWC Discharge sites.The Model is first calibrated at Tuini Pabra which covers sub catchment 1 and 2. Then it is calibrated at Yashwantnagar site for sub catchment 3. Next calibration site is Bausan for sub catchment no. 5, Haripur for Sub catchment no. 1, 2, 3, 4, Poanta sahib site for sub catchment no. 1, 2, 3, 4, 5 and 7.The model is then calibrated at HKB, the outlet of the study area covering all upstream sub catchments from 1 to 7. The results of Calibration and validation showing NSE daily and log flow duration and values of eight Metaparameters are tabulated in Table 2.During the calibration period the Mataparameters have been modified so as to improve the fit between the model values and the observed values.The objective function of NSE Daily & log Flow duration have been found between 0.69 to 0.85 for calibration and between 0.60 to 0.73 for validation period and the values are fairly good.

Water Balance at District Level
After running of both GR4JSG model and GR4J model, the Irrigator canal model (Crop model) is run.Then, water balance of all the districts in Haryana has been worked out.Water balance is evaluation of all the water inflows, outflows and utilizations of all canals in the district.The water balance of Karnal District is being explained in this paper to show how the water balance of each district has been worked out both for calibration and validation periods (Figure 4).all the four outflow links (Indri_2 link, WJC_MB upstream SYL link, NBK link to WJC-MB and the Default link 104) and also the downstream flow volume of all the above five supply points.It is important to add the downstream flow volume from supply points to the outflow.The reason we added the downstream flow volume from supply point to outflow is that all the water that comes to supply points is not extracted by the water user, so we need to add the downstream flow volume at the supply point to outflows.For Irrigator water user, the ordered water is calculated based on forecast rainfall and evaporation, so it is actually a forecast order.Water that will be released from upstream to meet the demand of the water user is based on this forecast order.But the extracted volume is based on real order which is re-calculated based on current time step's rainfall and evaporation.So, it could happen that water user orders more than actually needed, if it rains more and evaporates less on that time step.The total volume of water utilizations in Karnal district is the sum of the extracted volume of all the five supply points.The water balance of Karnal district has been worked to satisfy the following equation: ( )  4 and Table 3).The average inflows, outflows and utilizations in Karnal District for the period 1995-2013 come out to be 12,304,553.53 ML,12,222,694.84 ML and 81,858.69ML respectively.
It has been found that the above water balance equation has been properly satisfied for Karnal district.It proves that irrigator canal model has been properly set up.Similarly, the water balance has been worked out for the remaining districts of Haryana for calibration and validation period.The water balance of all districts of Haryana state has been shown in (Figure 5 and Figure 6).For all districts the water balance equation is satisfied.It shows that the crop model has been properly set up.The overall water balance of Haryana canal system for calibration, validation period and also overall period 1995-2013 is shown in Table 4.
Table 3. Water balance of Karnal district (Calibration and Validation).

Ground Water Fluctuations
The ground water levels are observed in June and October every year.The fluctuations in ground water level (depth) for the year 1974 to 2010 both for June and October have been shown in

Conclusions
In the paper a detailed study of hydrological models and irrigator canal model

Figure 1 .
Figure 1.Concept of the integrated water resources assessment model.

Figure 2 .
Figure 2. Yamuna basin up to Mathura and adjoining Haryana area (Source Google Map).

3 . 2 . 3 .
Discharge data of Central Water Commission of India (CWC) is used for CWC sites at Tuini (Pabra), Yashwantnagar, Bausan, Haripur, Poanta and Mathura.The flow data of various canals from year 1995-2015 has been collected from Haryana Irrigation and Water Resources Department.Crop and Soil Data Wheat, Barley, Rabi oil seeds, Gram and Sun Flower (Rabi crops) and Rice, Bajra, Maize, Cotton, Jowar and Sugar Cane (Kharif crops) are the main crops grown in Haryana state.The data of crop areas, crop factors, root zone depth, depletion factor, planting and harvesting dates, target soil depletion, field capacity, permanent wilting point, soil moisture capacity, soil types and ground water extractions etc have been obtained from Haryana Agriculture department, reference [1] [10] and Indian Agricultural Research Institute (IARI), New Delhi.The split-sample technique has been used in the study.The model has been calibrated from January 1, 1995 to December 31, 2005 and then validated from period January 1, 2006 to December 31, 2013.A warm-up period of one year (Jan-uary1, 1995 to December 31, 1995) was sufficient to remove any trends in the snow store states.The model was calibrated using calibration tool in the Source.
The water is supplied to Karnal district by five canals: Chautang feeder, Offtakes of WJC-MB, Nardak Gogripur offtakes of NBK link, other offtakes of NBK and Goli distributory to meet the irrigation water demands generated by the five water user node on each of the above canals.The total inflow volume in Karnal district is the sum of downstream flow volumes of all the four inflow links (In-dri_1 link, Default link 788, NBK link 53 and the Default link 103).The total outflow volume in Karnal district is the sum of downstream flow volumes of

Figure 4 .
Figure 4. Water balance of Karnal district (Calibration and Validation).
has been made to evaluate flows at different sites, water balance and water use efficiency at the district level.An integrated model of the state of Haryana and Yamuna catchment contributing flows in Haryana has been conceptualized.GR4J, GR4JSG and Irrigator models have been successfully applied to evaluate the modeled stream flow.The evaluation with observed stream flow was done using NSE daily and Log flow duration.The model calibration and validation are based on split-sample.The study has been successful in achieving the objectives of water balance and water use efficiency at district level for proper planning of water resources.The water use efficiency has been worked out to be low for Haryana under the case as usual.It is suggested to application of micro irrigation techniques such as drip and sprinkler irrigation.The water thus saved can further be utilized for increasing low flows in Yamuna downstream of HKB/Tajewala.

Table 1 .
S. S.Rawat et al.Area of elevation zones and hypsometric mean elevations.

Table 2 .
Calibration and validation results of at 6 CWC stations up to HKB.

Table 5 .
The overall fluctuations of Haryana state have been shown in Figure7.

Volume (ML) Year Calibration -District Mahendergarh
) Figure 6.Water balance of all districts of Haryana (validation).Inflows Outflows Utilization Figure 8. Validation at Mathura, the outlet of study area.S. S. Rawat et al.DOI: 10.4236/as.2018.94034509 Agricultural Sciences