Water Pricing and Implementation Strategies for the Sustainability of an Irrigation System: A Case Study within the Command Area of the Rakh Branch Canal

: The command area of the Rakh branch canal grows wheat, sugarcane, and rice crops in abundance. The canal water, which is trivial for irrigating these crops, is conveyed to the farms through the network of canals and distributaries. For the maintenance of this vast infrastructure; the end users are charged on a seasonal basis. The present water charges are severely criticized for not being adequate to properly manage the entire infrastructure. We use the residual value to determine the value of the irrigation water and then based on the quantity of irrigation water supplied to farm land coupled with the infrastructure maintenance cost, full cost recovery ﬁgures are executed for the study area, and policy recommendations are made for the implementation of the full cost recovery system. The approach is unique in the sense that the pricings are based on the actual quantity of water conveyed to the ﬁeld for irrigating crops. The results of our analysis showed that the canal water is severely under charged in the culturable command area of selected distributaries, thus negating the plan of having a self-sustainable irrigation system.


Introduction
Irrigation in Pakistan is indispensable for agriculture, as is true for most of the World's arid and semiarid environments. The area of irrigated land in Pakistan has grown from about 11.6 million hectares in 1947 to nearly 22.6 million ha in 1997 [1]. Therefore, irrigated farming accounts for 75% of total water withdrawals [2]. "Cost recovered" generally refers to the percentage of the full supply cost collected from the farmers. The full supply cost can be defined readily, owing to the ease in the estimation of its parameters (i.e., O&M cost and Establishment). However, the "full economic cost" and "full cost" are difficult to estimate due to the increased complexity involved in the estimation of the opportunity cost and other parameters [8].
In the case of Pakistan, water charges are based on either area or season-based flat rates. The present water pricing system does not encourage water use efficiency as it fails to relate the water price to the actual water volume supplied to the field [9]. Similarly, other major water pricing issues in Pakistan are the absence of linkage between water fee collection and service provision (lack of financial autonomy) and the lack of incentives for service providers to collect charges [10,11].
Several environmental laws, which support the use of economics as a key discipline to reach the objective of a self-sustainable irrigation system, have been introduced (e.g., [12,13]). The Water Framework Directive (WFD) approach hinges on the notion of full cost recovery for water services including environmental and resource costs from the farmers, which is precisely one of the most novel and interesting aspects in the norm. The WFD further strengthens the call for inculcating economic analysis to find the economic value of water for different alternatives to improve the water supply system for improved environmental status. [13] fully describes the implementation of WFD and explains the idea of full water cost coverage based on the concept of "water services".
Such analyses are urgent in regions like the Punjab province of Pakistan, where water scarcity is a critical issue. The province currently recovers only 20% of its irrigation water supply cost from consumers and intends to ensure full cost recovery as part of a new reform process [14]. The use of economic instruments such as full cost recovery and the analysis of measures applied to the irrigation sector need to take into account the level of demand for irrigation water in order to quantify the response of water use against the proposed measures (e.g., water saving as a response to increased water price). For this purpose, a detailed economic valuation of irrigation water might be useful to support the water pricing mechanism at the local scale.
Heal et al. [15] thoroughly discussed the methods available in literature to calculate the economic value of goods and services. The methods were later reviewed by Young [16] to estimate the economic value of water. Heal et al. [15] argues that in order to find the economic value of an unpriced input such as water, its contribution to the total value of the produced output should be isolated from the other inputs by subtracting their value addition to the product. Based on the data extracted from field experiments on crops in the 1970s and 1980s, the demand of irrigation water was estimated by statistically utilizing crop-water production functions [17,18]. The relationship between the output price and varying cost of water was studied to develop the demand functions.
Perhaps the most frequently used tool for estimating the value of irrigation water is linear programming (LP), which is based on the use of mathematical programing. The method presents an optimal solution by developing linearity between the constraints and objective function. The former is defined by the restrictions of the cropping area and water availability, whereas the latter is specified as the maximized return from several crops. The concept has been widely used for the estimation of the value of irrigation water; therefore, a large amount of literature can be found related to its application in the field of agriculture [19]. Chaudhry et al. [20] estimated the net return from water in Punjab province of Pakistan by using LP models. The study concluded that the marginal product value of irrigation water was high owing to its limited supply. However, the objective functions and constraints may not be linearly associated; hence the resulting regression equation may not be the true representative of the issue at hand [21]. Some researchers argued about the rigidity of LP in addressing several objectives of conflicting nature, simultaneously [22]. In this regard, Multicriteria Criteria Decision Making (MCDM) analysis is classified as a better alternative to LP (see [23] for a review). Researchers have proposed a number of multicriteria approaches to solve the problems in farming systems related to decision making (see [24] for a review of multicriteria methods in agricultural economics). A complete review of multicriteria methods was done by [25]. The method is unique in the sense that the most decisive factor in its implementation is not the sensitivity of the algorithms but the ways in which the process is organized, including stakeholder involvement and facilitation, and the ways in which the results are presented [26]. Berbel et al. [27] evaluated sustainability of the irrigation sector of Europe under various policy scenarios by generating farm-level water demand curves through MCDM. The authors claim that the methodology well grasped the complexities related to the European agriculture. However, the MCDM methods are criticized by a number of researchers because the methods are: (1) prone to manipulation, resulting in a false sense of accuracy [28] and (2) require the definition of a utility function, representing the preferences of decision-makers, through a complicated procedure [29]. A complete criticism of MCDM methods can be found in Mutikanga [30].
Occasionally, a hedonic method for estimating the value of irrigation water has also been employed and can be found in the literature to a limited extent. For example, Berbel and Mesa [31] used this technique in Guadalquivir Basin of Spain to estimate the value of irrigation water. Besides its numerous advantages, the major drawback of the hedonic method is its inability to translate the implicit costs into a meaningful financial sum.
Most recently, a technique called the Residual Value Method (RVM) has been employed by several researchers. The technique considers water as an intermediate input for the production of crops. The underlying assumption for the application of the model is that the farmer will use water to an extent where the net revenue generated per additional unit of water is equal to the marginal cost of the water [32]. Only a limited number of studies can be found in the literature making use of this technique [33,34] with almost no recently published studies. A complete review of this technique is provided by [16]. Using farm budget surveys data from 1987 to 1991, Bate and Dubourg [35] calculated the residual value of irrigation water applied to five crops in the eastern region of Anglia. Lange [32] used the same technique in the Orange River basin of Namibia. Moran and Dann [36] used secondary data sources in RVM to evaluate the economic value of water for the implementation of WFD. They further argued that the valuation of irrigation water will help in implementing the WFD. Speelman et al. [37] and Esmaeili and Vazirzadeh [38] applied the RVM technique for the valuation of irrigation water in South Africa and Iran, respectively.
In our work, we plan to use the RVM to execute the economic value of water for different crops. The estimated economic values define the income of farmers per unit of irrigation water applied to the crops. Afterwards, the full cost recovery figures are proposed for the operation and maintenance (O&M) of irrigation infrastructure by linking the water prices with the volume of water supplied to the field. Eventually, the implementation strategies are defined to ensure the same economic returns for the farmers during pre-and post-implementation phases of the full cost recovery prices.

Improved Water Pricing Mechanism: Why Now?
The province of Punjab, which is considered to be the agricultural hub of the country, is entitled to carry 55 million acre feet (MAF) of water according to Water Apportionment Accord 1991 [39]. Out of total Punjab's share, only 26 MAF (47.3%) reaches the farm gate while the rest is lost during conveyance in unlined water channels. These losses can be reduced from 60 to 80% by lining the water channels [40]. Over the last few years, a lot of investment has been made to reduce the post-farm gate losses by lining the water courses through successful cost sharing agreement between the government and farmers. However, no such investment is ever made at the national level to line the irrigation canals and distributaries, apart from few individual canals and distributaries which are lined recently (e.g., Lower Chenab Canal).
Unfortunately, due to the sociopolitical uncertainty in the country, no new water reservoirs have been constructed since the formulation of the Water Accord of 1991. With no new water-storage infrastructure and rising water demand, the water supply to each province is stagnant. The delivery of the allocated volume of water to the farmland takes place through the networks of canals and distributaries. These canals and distributaries require a huge financial capital per annum for their operation and maintenance (O&M). To recover the cost of O&M, a flat rate system of water pricing is currently in place in Punjab. Under the flat rate system, water charges are fixed per hectare of area sown in Rabi (winter) and Kharif (summer) seasons regardless of the type of crop being irrigated. The nominal rates of 123.55 PKR (Pakistani Rupees)/hectare and 210.035 PKR/hectare ($1 = 104.5 PKR) are charged for the winter (Rabi) and summer (Kharif) crops, respectively. The water charges are collected from the farmers by either the head of the village (locally called "Lambardar") or by the chairman of farmer's committee nominated by the Punjab Irrigation Department (PID). The flat rate water pricing system significantly increases the cropping intensities and production, but the system needs to be modified as it fails to ensure the sustainability of the irrigation system [41]. Our hypothesis is that irrigation water in Punjab is severely under charged and there is a need to revise the pricing system to ensure the sustainability of the irrigation system.
The long-forgotten issue of irrigation water pricing took center stage when the Asian Development Bank (ADB) and the Government of Punjab got involved in a process to finalize the loan of 2 million dollars to upgrade the Agriculture sector of the province. One of the conditions which ADB reiterated from its previous suggestions is that the Government of Punjab needs to improve the irrigation water pricing mechanism of the province [42]. Therefore, the aim of our research work is to find the value of irrigation water and then test the suitability of our newly developed pricing model which is based on the actual quantity of water supplied for irrigation. The model will suggest the appropriate prices of irrigation water for the study area to recover 100% of the O&M cost of irrigation system. The project will interpret the full cost recovery figures for three major crops (wheat, rice and sugarcane) grown in the area. It will also aid in understanding the issues and complexities related to the implementation of a full cost recovery mechanism. In the past, every effort to put in place such a mechanism was fully resisted by the average farmers of the area due to their delicate financial conditions, although the farmers were assured that the mechanism will be complimented with tax evasion in order to provide financial relief to them [14]. Therefore, the project will also assist in framing a mechanism to increase water prices in such a way that the average farmers are not financially over-burdened.

Study Area
The major sources of irrigation water in Punjab are: (1) canal water which contributes 40% of the total irrigation water (26 MAF); (2) tubewell water supplies 50% of the total irrigation requirements (33 MAF); and (3) rainfall accounts for 10% of water supplies (7 MAF) for irrigation [43,44]. Since our ultimate objective is to find the value of canal water, we therefore had to select a site where canal water is a major means of irrigation. The culturable command area (CCA) of the Rakh branch canal is selected for the implementation of the proposed methodology. The selected area is predominately irrigated by canal water due to better infrastructure which includes a lined branch canal, distributaries, and the majority of water courses. The study area is irrigated by four lined distributaries (Ratti, Gilluana, Khair ali, and Gatti distributaries) of the Rakh branch canal. The lining of irrigation channels ensures the delivery of irrigation water according to the designed criteria. The Rakh branch canal receives perennial canal supplies from the Lower Chenab Canal (LCC) system during the entire year of canal operation except December when the canal is closed for maintenance. The selected distributaries are located at the head, middle, and tail sections of the Rakh branch. A schematic diagram of the LCC illustrating the arrangement of distributaries is given in Figure 2. The area is defined as a mixed cropping zone, and analogous to rest of the country, the farmers in the study area are authorized to use the entire volume of water diverted to the water course by the distributary outlet, on a 7-day rotation period under the system called "Warabandi" [45]. Rice, wheat, and sugarcane are the main crops of the area while cotton and maize are not commonly cultivated in the area. The area is considered to be the most fertile area of the province with average output per hectare surpassing the national average. However, the area can still be more productive as the biological potential of the area is fairly higher than the present yield. The gap between average yield and biological potential for various crops in the study area is translated in the following Table 1. The details of the four distributaries irrigating a total of 10,615.136 hectares of the CCA are summed up in the following Table 2.  The Ratti distributary, being the smallest in our selection, has the lowest water allowance of 0.039 cumecs/1000 hectares while the other three distributaries have equal water allowance of 0.198 cumecs/1000 hectares. The allocated water is systematically delivered to the field through the outlets, engraved along the entire length of the distributaries, designed to carry only a specific amount of water. The magnitudes of discharge of each outlet in the selected distributaries are presented in Figure 3.

Value of Irrigation Water
The economic value of a service or good is measured according to its contribution towards the well-being of a given population [47]. The economic theory states that, analogous to land, water can also be treated as an input in the crop production [48].
This research uses the residual value of water for agriculture in the selected study area to study the allocation of the resource [16]. The application of the proposed model begins with the identification of variables that are defining features of agricultural property. Among other features, the geographical location of the study area is the core factor in defining all the parameters needed to implement our model. It should also be noted that increased size of the study area will certainly increase the number of dominating features that might affect the revenue generated by a farmer from agriculture productivity (e.g., local taxes on land and other inputs). The model conceptualizes around the fact that the agriculture production is the function of input variables. Mathematically, the model can be represented by the following equation: By considering technology available to the farmers in study area as constant and all other features as variables except water, the total revenue generated by a farmer from agricultural output can be written as: where (P Y ·Y) represents revenue generated from product Y for unit surface (e.g., hectares) and (V MP i ) is the value of the marginal product of each feature (X i ). The VMP defines the incremental change in output value attributed to a unit change in any single input item. For details on VMP, the readers are referred to [49]. The above equation can be made operational by assuming the following three points: (i) total revenue generated by the farmer is the function of capital invested on inputs except water; (ii) the inputs will maximize the output (maximizing behavior i.e., V MP i = P i ) and (iii) the costs of agricultural inputs are provided by the local markets. Through simple rearrangement of Equation (2), the value of irrigation water can be calculated as: where X W represents the actual quantity of water consumed by the crops. This variable is very important for the areas where irrigation water is not present in adequate quantity. It can be expressed by the definition of maximum irrigation needs, which is obtained by subtracting Evapotranspiration (ET) from usable (or effective) rainfall computed from an agronomic method based on FAO publications and the actual quantity of irrigation water given to crop. The farmers of areas with a surplus amount of water supply usually supply irrigation water greater than or equal to ET requirements of crop.
The amount of water supplied in excess to ET demands of a crop does not increase the yield and so it has no value beyond this point. The collected datasets both from the field and government departments are used to estimate the numerators of Equation (3), whereas the denominator, which represents the actual quantity of irrigation water supplied to the crop, is estimated by using the irrigation applied to each crop from sowing to harvesting and the depth of irrigation. While the major difference is observed in the values of X M , the values of other variables of Equation (3) are similar for the entire study area, predominately due to the similar farm practices and technology available to the farmers. The reason for the disagreement in the value of X M is the variable quantity of subsurface water extraction by using diesel motor pump tubewells which significantly increases the production cost associated with X M . The diesel driven tubewells in the study area consume 3 L (85.33 PKR/L) of diesel per hour to discharge 0.021 m 3 /s of water and require 12.355 h to irrigate one hectare of land. It is worth mentioning that due to the recent energy crisis in the country, the electric tubewells have become obsolete and are not used for irrigation in the study area [50].
The economic data for the application of the proposed methodology are based on the value and cost incurred by the farmers for the production of each crop. The field surveys are carried out at the head, middle and tail sections of the selected distributaries to collect the data on farming practices, production per hectare and cost incurred for the production of crops. Altogether, 120 farmers are interviewed. The sample is carefully selected to replicate the provincial average-to-progressive farmers proportion of 90 to 10% [51] because the full pricing mechanism is intended to be put in place in the study-area, evaluated, and extended from the pilot area to the entire province as per the directions of ADB. The progressive farmers, by definition, are considered to be the efficient growers producing relatively higher yield per hectare by eagerly adopting scientific methods of growing, whereas the average farmers adopt traditional inefficient methods of farming, resulting in a loss of precious resources and lower output. The output per hectare, in kilogram, against the resources spent on procuring inputs (in Pakistani rupees) is shown in Figure 4. It can be interpreted from Figure 4 that the farmers of the study area can be broadly categorized as: (1) average farmers and (2) progressive farmers. The classification is based on the farm size with the latter having a larger size (above 10.12 hectares) than the average farms (of less than 10.12 hectares) [52]. The progressive farmers, on average, are yielding 62.5%, 37.5% and 10.77% more yield per hectare (for wheat, sugarcane and rice, respectively) as compared to the average farmers. The significant increase in the yields are due to a nominal increase in production costs of 13.680, 10.382 and 12.888% for wheat, sugarcane and rice, respectively. The progressive farmers spend the difference in cost on: (1) sowing only certified and graded seeds; (2) soil testing to evaluate the fertility status in order to rationalize the optimum requirement of fertilizers for good yield; (3) consistent soil moisture monitoring for timely application of irrigation water; (4) usage of weedicides to prevent pre-harvest losses; and (5) effective utilization of machinery to prevent post-harvest losses.
Although progressive farmers are producing higher yields, even the present production levels are notably lower than the identified biological potential of the region (see, Table 1). The major constraints that bar farmers from reaching the biological potential are a lack of appropriate technology, inadequate investment in irrigation infrastructure, and lack of quality inputs and resources.
Besides production cost and revenue generated per hectare, the information related to crop rates is trivial for the implementation of Residual Value method (RVM). It was acquired from the official data of Ministry of agriculture and was further cross-validated during the field surveys, whereas the land rent statistics necessary for the execution of RVM are gathered from the "Lambardar" of the area falling within the CCA of the selected distributary and the farmers acquiring land on a rent basis.
However, the fundamental quantity required in the pricing of irrigation water is the volume of water that ultimately irrigates the farmland. The distributaries are the last point where the live flow data is recorded. The PID records data at the head and tail sections of each distributary. The data have recently been made available online at the official website of PID (irrigation.punjab.gov.pk/Search.aspx). Each distributary and its corresponding outlets acquire a fair share of water in case the main perennial canal is operational.
We conducted field surveys in the study area to determine the amount of irrigation that farmers are applying to each crop-from sowing to full maturity. The information regarding maximum evapotranspiration requirements of each crop and average dose of water provided to each crop are summed up in Table 3. Table 3. Irrigation water supply and relative use by crops in Rakh Branch Canal.

Crop
Average Dose (m 3 /hectare)  Table 3 shows the availability of water resources in the Rakh branch canal command area. We defined the ratio of water applied and water requirement of a crop as the Relative Irrigation Supply (RIS). For example, the average dose of wheat is 4065.04 m 3 against the maximum ET requirement of 4575.64 m 3 which results in an average RIS value of 0.88 (i.e., 88% of ET requirements of wheat are fulfilled). It can also be observed that the sugarcane crop (RIS = 1.408) is over irrigated and there is still some room to save water by improving irrigation practices. Besides sugarcane (ETP max = 108,040 m 3 ), the ET requirement of rice crop is very high (16,089.4 m 3 ) although the duration of the rice crop is only 100-120 days. The average water doses in the study area, as stated above in Table 3, are obtained from the multiple sources, viz. canal water, ground water and rainfall. Due to limited surface water supply in the face of an ever-increasing demand, the increased magnitude of rainfall will minimize the extraction of groundwater but will have no impact on the surface water supplies.

ETP-Max (m 3 /Hectare) RIS-Ratio
To determine the contribution of rainfall in irrigation, the effective rainfall depth for the study area is calculated by using the agronomic method outlined in FAO publications [53,54]. The effective rainfall depths for different crops in the study area along with the actual rainfall amount and the length of growing season are represented in Figure 5.
The horizontal lines in the Figure 5 represent the length of the growing season for each crop while the vertical lines quantify the contribution of effective rainfall depth to irrigation. Sugarcane is generally sown in February and is harvested in January. It therefore has the longest growing season and the greatest contribution of effective rainfall depth to irrigation. On the contrary, paddy crop has a span of only 100-120 days in the study area, falling in the summer monsoon season, and therefore has a reasonable effective rainfall depth as compared to the wheat crop which spans over 6 months__ 4 months in winter and 2 months in spring. Since the winter season receives considerably less rainfall, the effective rainfall depth is the lowest for wheat among the entire selection of crops. Assuming that the effective rainfall depth is constant over the entire study area, the percentage contributions of canal, tubewell and rainfall in irrigating wheat, sugarcane and rice crops cultivated in the CCA of the selected distributaries are displayed in Figure 6.

Pricing Irrigation Water
Within the context of the full supply cost, the Punjab Irrigation and Drainage Authority Act (PIDA) 1997 clearly states that the full supply cost of irrigation systems is to be recovered from the end users within the next 7 to 10 years [41]. Table 4 lists the total full supply cost, including its components (i.e., establishment cost and actual O&M cost), starting from 2003-2004 to 2016-2017. It can clearly be seen from Table 4 that the Irrigation Department collected only a fraction (about 20%) of the full supply cost while the difference is shared by the provincial government. Strictly speaking, a reasonable amount of 111,468 millions PKR is paid by the provincial government in the last 14 years (2003-2016) at an annual average of 7962 millions. In the given span of time (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), the full supply cost is increased from 4263 to 11,398 millions to the corresponding decrease of 1205 to 1017 millions for the revenue collected. Conclusively, the irrigation related infrastructure is persistently underfunded which speaks volume for its financial unsustainability [9]. It is also a matter of great concern that almost 70% of total full supply cost is consumed in operational needs (establishment cost) compared to only 30% which is actually spent on the maintenance of infrastructure.
For the sustainability of the irrigation system, it is inevitable to collect the full supply cost from the end users. To this end, we linked the irrigation pricing per hectare (P i/W ) with the full supply cost by using the following newly developed equation; where; The above Equation (4) consists of the product of the total volume of water in m 3 applied during irrigation (Q d * T * 60), number of irrigation events ( L GS W T ), and the maintenance cost (in PKR) required per unit volume of water (η C ). The assessed allocated irrigation time (T) in Equation (4) defines the actual time for which the CCA of a specific outlet is irrigated during each irrigation event in Warabandi. Figure 7 shows the values of "T" for each outlet of the selected distributaries, while the CCA of each outlet is shown in Figure 8.   Table 4, Figures 7 and 8, whereas, the availability of 26 MAF of water at the farmgate of Punjab and the 7-day water rotation cycle in the study area set the values of V F/G and W T , respectively.

Economic Characterization of Irrigated Crops in Rakh Branch Command Area
Equation (3) is solved for the CCA of Rakh branch for three major crops (Wheat, Sugarcane and Rice) and the results for average and progressive farmers are presented in Table 5. Table 5 describes the residual value of water for the considered crops in the study area for 2016-2017 data. The "apparent water productivity" is executed by dividing "total sales of crop per hectare" or "total income generated per hectare" by water usage, whereas, the residual value of water is calculated by excluding all the other paid factors by the farmer (e.g., Land rent).
It can be clearly noticed that (out of three considered crops) wheat has maximum apparent water productivity and residual value as compared to sugarcane and rice. From this it should never be concluded that irrigating other crops is irrational because crops cannot be analyzed in isolation [55]. For example, if irrigating sugarcane is considered in isolation, it may seem as if due to low value of water, sugarcane should not be irrigated. However, the whole sugar industry, which generates thousands of direct and indirect jobs, is dependent on this crop. Moreover, different crops are needed to optimize the agronomic system as a whole.

Irrigation Water Pricing for Full Cost Recovery
Having set the values of parameters in Equation (4), we calculate the annual price of irrigation water per hectare for the CCA falling under each outlet of the selected distributaries. The revised prices to recover the full cost of O&M are recommended in Figure 9A and are compared with the case when recovery cost includes infrastructure maintenance cost only (see Figure 9B).
It can clearly be interpreted from Figure 9A that the irrigation water in the command area of Rakh branch canal is severely under charged ever since the introduction of the flat-rate irrigation pricing system. The maximum deficit is observed for the Khairali distributary, whereas the Ratti distributary showed the least deficit. It can be further noted from Figure 9B that the prices for each distributary are reduced to almost one-third if only the infrastructure maintenance cost is recovered from the farmers. Apart from the annual water prices for the CCA falling under each outlet discussed in Figure 9, the crop-wise rates considering the volume of water each crop consumes from sowing to harvesting in the study area are also analyzed and summed in Table 6. The above table reports that the seasonal crops are consuming water far above the set seasonal prices in the CCA of all the selected distributaries. The situation is even worse for Gilluana, Khairali and Gatti distributaries as present water charges are not even recovering the O&M only, let alone the full cost recovery. On the contrary, due to comparatively lower water allowance, the Ratti distributary has not only the lowest deficit between cost and recovery but also recovers the O&M cost.
The advised rates are very high as compared to the current water charges. However, the proposed prices are nominal when compared to the pumping cost of groundwater using tube-wells. It is estimated for the case study area that 1000 m 3 of canal water even with modified rates costs 374.3 PKR against 3113.6 PKR needed to extract the same volume from sub-surface resources.
The seemingly obvious drawback of shifting from flat-rates to a full cost recovery model is a noticeable dip in the farmers' income. Therefore, the resistance of the average farmers against the full cost recovery model seems logical. The resistance can be minimized effectively by facilitating the farmers and ensuring the same profitability of the farmers even after the introduction of the full cost recovery system. It can be ensured by either one or a combination of the following options: (1) tax evasion; (2) decreasing conveyance losses or increasing availability of canal water; and (3) increasing average productivity per hectare. The response of water value against these options is mapped in the following panels of Figure 10.  The green shade in Figure 10 represents the extent to which these compensations can be extended in the present circumstances related to the agricultural sector of the province. Against the traditional perception, tax evasion has the least impact on the improvement of the value of water, whereas increased productivity exponentially improves the value of water for three crops. Moreover, although tax evasion can be an instantaneous relief, it will directly put pressure on public money, whereas, the increasing conveyance efficiency or constructing new reservoirs can take a longer time because the option requires lined water channels or public consensus on the issue (of constructing new water reservoirs). Therefore, the most feasible option is to increase productivity of the average farmers. The increased productivity is not only possible; physical examples of progressive farmers can be found in the region. Bringing production to the level of progressive farmers will not only increase the water value but will also increase the farmers' revenue even after increased water pricings. During the field survey, it was noticed with a grave concern that the lack of awareness and lack of intent to change are the core reasons preventing the average farmers from adopting a more functional frame of reference. For example, most of the farmers acquire the seeds from private dealers instead of buying from certified dealers of provincial seed corporations, thus compromising the productivity. Moreover, the average farmers seldom get their soil quality checked and continue using inputs based on their judgement. The continuous crop production squeezes the nutrients out of the soil and thus with the passage of time, the demand for fertilizer increases [56]. These small factors might not be seen as significant when observed in isolation but they impact greatly on the productivity when analyzed in totality. Therefore, following these rules can reduce the difference between average and progressive farmers. Alongside the aforementioned reasons, another reason for lower productivity is the continuous shrinking of average farm sizes in Punjab, which decreases efficiency of agricultural farming and ultimately the profitability [57].

Discussions
The application of this method has the advantage of simplicity and the easiness of replications once the values of variables defined in pricing equation are known for the outlets of a selected distributary. Evidently, the method defines the price of irrigation water as the function of volumetric quantity of water applied to the farm land each year. Therefore, the areas with lower water supply are charged less. For example, the Ratti distributary has a lower water allowance than the other three distributaries; therefore, the proposed prices for Ratti are lower, thus signifying how effectively the logical connection between water pricing and quantity supplied is explained by the proposed model.
The proposed water prices are 5.0%, 7.9% and 7.5% of the value of water in Gilluana, Khairali and Gatti distributaries and 6.0%, 9.7% and 8.0% in Ratti distributary for wheat, sugarcane and rice crops, respectively. The recommended prices are 3 to 9 times higher than the existing water prices against the value of water. Moreover, the deficit between full cost recovery and collection is increasing with the passage of time (see, Table 4). This difference should be minimized to achieve sustainability, keeping in mind the thin financial conditions of the small farmers which constitute the majority in Punjab's agriculture and are mostly going to be affected by an increase in water price. One recommended solution is a phasal increase in water cost with initial emphasis on either full cost recovery from the progressive farmers owning large farms or recovering 100% of infrastructure maintenance cost only. However, with the involvement of the majority of the rural work force in the field of agriculture, there are also political reasons for the significant under-charging of canal water. Therefore, the adoption of proposed water prices will be challenging and the implementation needs to be systematic so that the sales and incomes of the farmers are well in line with their logical expectations. There should be a long-term political consensus across the board to implement the plan of full recovery prices. The plan should then be rigorously followed by the elected political administration of the province.
Along with the efforts to put in place a full pricing mechanism, there also needs to be a policy-support mechanism to reduce the O&M cost of irrigation system. For example, the "establishment cost" which is 70% of the total O&M cost (see Table 4), can be reduced by empowering the farmers through decentralization and involving them in irrigation infrastructure development and maintenance. The most relevant example in this regard is the California State Water Project in which Farmer Organizations (Irrigation Districts) are responsible for overseeing a number of field tasks (e.g., collecting the O&M cost of the irrigation systems) [6]. The same system, if followed in Punjab, can significantly reduce the "establishment cost". Another pertinent solution is to create a consensus at the national level on the construction of new reservoirs and developmental projects to line the irrigation channels because increased availability of canal water at the farm-gate will bring down the prices per unit volume of water. Also, the improvement in per hectare productivity, which is presently on the lower side, will reduce the resistance against the full pricing mechanism. The major reason for lower productivity per hectare of average farmers is the shrinking farm size which can be properly addressed by encouraging corporate farming [58].
Finally, we urge that implementing any newer policy without winning hearts, minds and trust of farmers will be counterproductive. Farming, in general, is becoming more and more challenging with the unpredictable climate and out of all the businesses which are based on agriculture, only farming remains most vulnerable. Farmers know the importance of water in farming business and fully understand the importance of its sustainability. For example, lining of water course project was a massive hit only because of its visible benefits to the farmers. We believe that the farmers should be taught the importance of pricing the irrigation water for the sustainability of an irrigation system before increasing the prices of water. It can be concluded that any system developed for the sustainability of an irrigation system will have the full support of the farming community.

Conclusions
Water is the life line of agriculture, yet the sustainability of the Irrigation system of Punjab, Pakistan cannot be guaranteed. We believe that instead of the presently adopted flat rate system, the pricing of irrigation water should be linked with the pressure which the irrigation system has to sustain each year to ensure a proper pricing model that hinges on the concept of per unit water usage [59,60]. Increasing the price of water will encourage farmers to invest more into water saving technologies (HEIS), thus significantly improving the water productivity. However, the farmers of Punjab are already overcommitted and any unplanned increase in water cost would further aggravate their problems. We therefore suggested a number of solutions for the implementation of proposed methodology.
The results of this project will give insight into the practical implementation of the much delayed PIDA Act 1997. Apart from the sustainability of irrigation networks, the method can be helpful in the allocation of surface water to the intra-sector (tree irrigation versus herbaceous crops) by observing value per unit of water and allocating more water to high water-value crops. However, the estimation of the value of water for a larger study area would require the proper knowledge of area-based exceptions which cause disproportionality in the measurement of the value of irrigation water. These exceptional highlights also play a vital role in setting up cost-efficiency analysis which fundamentally needs to be based on the economic values of water.
We recommend that the proposed analysis should also be done for the other parts of the country to obtain a holistic picture of sustainability of the irrigation system. Moreover, future research works should also specifically focus on devising such policies or strategies which can rationalize the establishment cost resulting in a significant reduction in the total O&M cost of the irrigation system.