1. Introduction
The irrigated agriculture sector is the prime user of freshwater resources around the world and consumes approximately 69% of the freshwater withdrawal [
1]. Asia has the largest consumption of around 56% of global fresh water for irrigation purposes [
2]. Freshwater is becoming increasingly scarce in general and in Asia in particular [
3]. The existing canal irrigation system in Pakistan is a supply-based system working on the principle of equitable distribution independent of the crop cover [
4]. Major constraints faced by irrigated agriculture in Pakistan are canal water scarcity and its uneven distribution. Irregular rainfall pattern and water logging also have severe impact on the crop productivity [
5,
6]. The water resources of Pakistan, both groundwater and surface water, have become inadequate to fulfill the growing demands of the irrigation-based agriculture sector [
7]. Canal water is not sufficient to solely satisfy the crop water needs as it fulfills only 37.5% of the crop demand in Punjab. Rainfall contributes 15.5% and groundwater contributes 47% of the remaining demand. Crops receive only 70% of their demand in Punjab despite consumptive use of all sources including rainfall [
8]. Punjab province, the largest irrigation area of Pakistan, is an intensely cultivated region covering an area of about 8.4 million ha [
9]. This area observes cultivation throughout the year in two crop seasons namely summer (kharif) and winter (rabi). The high flow irrigation period falls into the summer from June to August, while late winter season from February up to early June is observed as low flow season [
10]. A rotational 8-day irrigation scheme is implemented to equally satisfy the water need in all irrigation districts [
11]. There is a realization globally of the need to efficiently utilize water resources to avoid crop water stress and increase crop productivity [
12].
Crop water consumption, generally known as evapotranspiration (ET), can be calculated accurately at pixel level using remote sensing technology [
13]. ET is the key parameter used for improvement in agriculture water management. Several methods are being used for direct or indirect measurement of crop ET; these include weighing lysimeter-based measurements, eddy covariance method, Bowen ration surface energy balance (BREB), surface energy balance algorithm for land (SEBAL), and using reference ET by multiplying it with crop coefficients [
14]. Estimation of reference evapotranspiration based on crop coefficients is being used commonly for farm level agriculture management [
6]. Reference evapotranspiration (ETo) is mostly calculated using weather parameters by common equations known as the Hargreaves equation and the Penman–Monteith equation. The Hargreaves equation requires minimum and maximum temperature as input for ETo calculation, while the Penman–Monteith equation requires temperature, wind speed, humidity, and sun radiation as input weather parameters [
15]. ETo calculated through Penman–Monteith or Hargreaves equations reflect the values of a reference crop in ideally suited climate conditions and needs to be converted to crop specific consumption, i.e., actual evapotranspiration (ETa) using crop coefficients [
9]. Crop coefficients reflect variable water need as per stage of the crop from sowing up to harvesting and vary for each crop as per its leaf area density and climate conditions [
16].
The crop coefficient approach is considered as reliable and efficient in determination of the spatiotemporal growth of crop water needs and crop specific variations show strong correlations with the satellite derived spectral index, normalized difference vegetation index (NDVI) [
6,
17,
18,
19]. Through regression calculation between crop coefficients and the NDVI, a set of crop specific parameters is created that allows the calculation of reflectance crop coefficients (Kcr) [
17]. By applying Kcr to both, an ideal reference crop and the monitored crop, the difference between these two can be translated into crop water stress. This approach is quite significant as it provides crop coefficients (Kc) for each stage of the crop rather than providing a static estimate of Kc found in published reports or manuals. This method includes the data of land use for different dates throughout the crop growing season derived from the freely available satellite images of the whole year. This method is considered cost-effective in terms of data requirements, i.e., satellite images are available freely and climate data are available online in open access in comparison to the data collected through field surveys [
20].
The spectral crop coefficient approach was used to monitor crop water consumption and to test the accuracy of crop coefficients for three fields in the Texas High Plains, all planted with cotton. Comparison was made of crop water use calculated from the reflectance-based crop coefficients approach and eddy covariance measurement approach (measured from a flux tower). In investigations, researchers used Landsat data for NDVI calculation, and ETo was calculated after the Penman-Monteith equation. They concluded that the use of reflectance crop coefficients is effective under a variety of different irrigation conditions and also stressed the adaptive character of the approach accounting for specific field conditions as opposed to standard crop coefficients [
21]. NDVI-based Kc calculation approach is significant in crop water requirement estimation. The crop specific water requirements are expressed in terms of Kc and have been quantified for different crops. Water requirements can then be found through multiplying crop coefficients times the reference evapotranspiration [
22,
23,
24,
25,
26]. Using modeled potential evapotranspiration (PET) and the relationship between NDVI and Kc, NDVI differences between ideal and monitored crop can be translated into irrigation water needs. Successful applications of this technique or approach have been reported in many studies and projects [
25,
26,
27,
28].
The purpose of this study is to provide the spatial distribution pattern of crop water requirement by monitoring the crop health in the Lower Bari Doab Canal (LBDC) command of the Punjab province. The results from this study would provide useful information for allocating water in main and secondary canals in order to avoid crop water stress and provide equitable water across the canal command. In summary, the objectives of this study are: (i) identification and quantification of different cropping patterns in LBDC; (ii) developing a decision support model for crop water efficacy at monthly intervals by using hydro-meteorological, geographical, and solar parameters; and (iii) scrutinizing the potential water deficit along the canal by counting upon the contribution of irrigation supplies and rainfall for both rabi and kharif seasons.
3. Results
3.1. Cropping Pattern in LBDC
The vegetation index (NDVI) time series approach generates a total number of 12 land covers in the study area including six crop types. The classification scheme follows the crop calendar and defines summer (rabi) and winter (kharif) season crops separately.
The time series NDVI profiles of each class behave differently due to the intrinsic properties of crops. For example, wheat and rice are the perennial crops of rabi and kharif seasons substantiated in the study area, respectively. The NDVI profiles of both crops are closely related to each other, but the seasonal crop calendar (sowing and harvesting dates) is completely different. However, different weather conditions may lead to slight variations in the NDVI profile of same crops among different years. The results show that wheat, maize, and potato crops cover a major portion for rabi season in the LBDC region illustrated in
Figure 8.
On the other hand, maize and cotton are the major crops of the kharif season as shown in
Figure 9. The area comprises a combination of perennial and non-perennial crops as shown in
Table 2.
Spatial distribution shows that the northeast (NE) region of the LBDC command comprises maize; however, a huge concentration of cotton was found while moving towards the southern side of the area. The spatial distribution describing rice, fodder, and sugarcane classes covering small portions in the study area as shown in
Figure 8 and
Figure 9.
3.2. Crop Calendar and Coefficients of LBDC
The sowing and harvesting periods and stages are different for different spatial extents. However, the LBDC area comprises a large spatial extent and has variable stages of crops. The annual crop cycles for each crop in the study area were collected from the field survey and by considering the local expert knowledge. The crop calendar describes the temporal variation of individual crop growth over a year in the study area. Different stages of crops according to their sowing and harvesting time periods are presented in
Table 3. The reflection-based crop coefficients (Kcr) were generated for individual crops using Equation (2). The average Kcr values of individual crops for years 2014–2016 are shown below in
Figure 10.
The smaller NDVI values resulted in smaller Kc values and vice versa. The maximum value of Kc for each crop was observed between the middle to final stage whereas the initial stage reflected a low value. For instance, following the rice–wheat crop combination, the Kc value curve attains maximum value for wheat crop in the months of February (middle) and March (final stage). However, the same line goes to minimum in the months of May and June when wheat is harvested, and it again starts to climb when rice is cultivated. It again attains peak and goes to bottom after the harvesting of rice. The same is the case with other crop combinations corresponding to their sowing and harvesting periods. The case of sugarcane is different as it is considered a whole year crop, and its Kc curve attains peak value depending upon its sowing period.
The reflection-based Kc trend corresponded with the stages of crops (initial, middle, and final stages) tabulated in the crop calendar as shown in
Table 3, which was prepared using crop phonological profile. Crop calendar was also verified during field survey by interviewing farmers.
3.3. Spatial Distribution of Meteorological Parameters (ETo, T, and P)
Temperature, evapotranspiration, and precipitation are very crucial parameters in calculation of crop water requirement. The parameters for the quantification of water deficit are used on a monthly basis for the years 2014, 2015, and 2016.
Figure 11 shows the variation of evapotranspiration, temperature, and rainfall for the years 2014, 2015, and 2016.
The results show that high values of evapotranspiration, temperature, and rainfall lie between the months April-August, April–September, and July–September, respectively, for the years 2014, 2015 and 2016. Maximum values of evapotranspiration and rainfall were observed in the months of June and July, respectively. These two parameters are in inverse relation to each other by indicating a decreasing trend of ETo compared to the rise in rainfall in the month of July. However, with the increase in temperature, the values of evapotranspiration also increased as shown in
Figure 11, which indicates a positive relation between the two parameters.
Chronological behavior of evapotranspiration, temperature, and rainfall shows variability in the spatial distribution pattern in LBDC as shown in
Figure 12. It has been observed that temperature and evapotranspiration have an increasing trend from east towards west in the study area. On the contrary, precipitation has a decreasing trend moving from west towards east. In general, similar pattern and correlation between evapotranspiration and temperature parameters have been observed for all the months of three years in the study area.
3.4. Crop Water Deficit and Requirement
The incorporated meteorological parameters and crop cycles resulted in chronological crop water deficit assessment in the LBDC region. Seasonal mean water deficit for the years 2014, 2015, and 2016 has been mapped in the study area as shown in
Figure 13. It depicts the spatial distribution pattern of water deficit for rabi and kharif seasons in millimeters. Generally, the same spatial distribution pattern of water deficit for rabi and kharif seasons was observed in the study area. However, a deficit value in the kharif season is much higher than rabi season for all three years. This is because the kharif season lies in the hottest months of the year (i.e., June and July), which causes maximum evapotranspiration resulting in high water deficit.
Figure 13 illustrates extremely high deficit observed at the tail (southwest) of the study area due to less irrigation supply because of maximum conveyance losses occurring and being located far away from the canal head. Similarly, southeast parts show high deficit for both seasons of these three years. These areas also lie along the border of the canal command and at the tail of the secondary canals, whereas northeast and northwest parts face considerably less deficit due to the advantage of lying close to the canal head and near to the Ravi river.
Irrigation supply, groundwater, and rainfall are the available sources to fulfill the crop water needs. However, groundwater pumping is not feasible or available for all farmers on demand, and small land holders usually rely on canal supply and rainfall. Influence and correlation of irrigation supply and rainfall parameters with crop water deficit measured on a monthly basis in LBDC are shown in
Figure 14.
The results show that less than 10 mm mean deficit is observed in the winter months (October-February) for the years 2014, 2015, and 2016. The low values of water deficits are mainly due to the lower evapotranspiration rate because of temperature falling down in these months. Specifically, December and January show relatively the least deficit due to the teleconnection of water deficit with evapotranspiration and temperature. Irrigation supply is not available in the months of January due to annual canal closure, but rainfall has supplemented to fill the gap. The crop calendar and trend graph of Kc also complement the deficit result because rabi crops are at an initial stage of their growth period during these months and require a substantially lower amount of water as shown in
Figure 15.
The annual cycle of water budget for 2014 (
Figure 15) shows that more than 40 mm of irrigation water was supplied to the LBDC throughout the year except for the months of January and February. This constant volume of canal supply is not sufficient for the crops in the peak summer months of July and August. Although aid of precipitation, around 45 mm, was available in the month of July, it was not able to meet the crop demand. However, it is obvious that when the peak summer months were over and temperature started to fall in the month of September, rainfall was also available in a good proportion of 50 mm, the deficit fell considerably. The deficit values dropped to a range of 5–10 mm from the month of October onwards because of the sufficient decrease in Kc values, hence low water requirement of crops.
In the water budget graph of the year 2015 (
Figure 16), water deficit did not attain peaks in peak summer months of June to August, because the LBDC command area received high precipitation in these months. It was observed that a very high value of rainfall (>100 mm) for the month of July restricted the deficit to as low as 5 mm. However, the irrigation supply remained in the range of 40–45 mm for April-November. It is obvious from the water budget graph of the year 2015 that rainfall volume in the study area directly influences the crop water deficit. The water deficit pattern in the winter months is much more similar to the year 2014.
Similar rainfall and irrigation supply patterns were observed in the year 2016 as shown in
Figure 17. However, rabi crops are observed with a shortage of water in the peak demand month of March due to a lower amount of rainfall as compared to years 2014 and 2015. High volume of rainfall in June and July for the kharif season kept the deficit under control, similar to the year 2015.
From results of all three years, the fact is established that the study area observes relatively high water deficit in the kharif season due to higher temperature, which causes maximum evapotranspiration in the months of June, July, and August. Although irrigation supply remains relatively high during these kharif months, crop demand is not satisfied unless sufficient volume of precipitation is available or groundwater is pumped in ample proportions. For the rabi season, the months of March and April are critical to meet crop water demand because of wheat being the dominant crop at its peak growth. Moreover, increase in ET values and decline in precipitation during peak demand months of rabi and kharif seasons favor water deficit in crops. On the contrary, the months of December, January, and February, being the coolest months of the year, face minimal ET deficit in the study area because crop water demand is fulfilled by the precipitation and irrigation supply and groundwater abstraction.
4. Discussion
Canal water scarcity and inefficient management of irrigation water channels are resulting in maximization of the use of groundwater to fulfill the crop water needs. Groundwater pumping is at maximum during the driest months of the year because of the higher rate of evapotranspiration. This fact is also confirmed in different published research that shows that in south Asia, Pakistan is the country with the highest groundwater usage of around 53% [
38]. It is imperative to improve crop water allocation by identifying the particular cropland under stress in irrigation districts and disseminating it to canal managers and farmers in a timely manner. Visualizing the impact of spatial precipitation variations is necessary to rationalize canal water allocation, which is possible by establishing an approach that is capable of providing estimation of the spatiotemporal dissemination of crop water requirements [
39]. Availability of spatial and temporal crop distribution data is one of the major components in devising policies for the agriculture sector regarding sustainable irrigation water management. A decision support tool for water resources management should be developed by acquiring the understanding of agricultural practices and the availability of reliable assessments of crop water requirements [
40]. Particular challenges in the development of a crop water management system are the size of the Punjab irrigation area, spatial and temporal crop variability, the location specific stress situation of crops that cannot be compared to standard crops, and the developmental status of technically viable monitoring techniques/models. Spatial analyses at such scale and at such short repetitive steps are only feasible with remotely sensed images.
This pilot research study attempts to provide a cost- and time-effective approach for efficient monitoring of crop health over one of the largest canal command areas in Punjab, Pakistan. The major goal of this study was to develop a spatial decision support tool to integrate and standardize the various parameters for estimation of crop water needs. Multiple datasets including hydrometeorology, irrigation, and remote sensing satellite information collected from different sources including organizational and open-source platforms were employed in GIS and remote sensing environments. Reference evapotranspiration (ETo) and crop coefficients were of major concern for quantification of crop water requirements in the study area. The Hargreaves method, considered as simple and reliable, especially in data poor areas, has been employed for calculation of reference evapotranspiration (ETo) [
35,
41], while reflected crop coefficients (Kcr), which proved highly correlated with the NDVIs, were derived using MODIS NDVI product and FAO defined crop coefficients (Kc) for calculation of reference and actual crop coefficients [
42]. Crop type and pattern were identified using Landsat imagery for accurate crop water quantification.
The research results verified the worth of satellite data for proficient appraisal of crop health monitoring. They also highlighted the worth of open-source information for crop water requirement estimation. The multi-step scheme for near real-time irrigation water stress estimation is quite effective in nature for large scale studies with limited data availability and therefore could be employed on other regions of the same agro-climatic conditions.
5. Conclusions
This study presents a technique for crop water deficit modeling in data poor areas and concluded the following:
(i) NDVI-based crop classification revealed that wheat and cotton are the major cultivated crops for rabi and kharif seasons with an area of approximately 60% and 30%, respectively. However, maize crop is non-perennial in nature, showing a variable cultivated area in rabi and kharif seasons. Sugarcane and rice, being the perennial and non-perennial crops, respectively, are also cultivated in a small portion of the study area. It is clear from the crop statistics that the kharif season has more variation in crop cover than the rabi season, hence requiring more dynamic water allocation.
(ii) The reflection-based Kc values remained consistent with crop calendar of LBDC generated by the local experts and published reports. The Kc values remained high at middle and final stages for both rabi and kharif season crops due to the maximum NDVI value, hence maximum crop water requirement at these stages.
(iii) The LBDC command area receives a continuous canal water supply throughout the year, except the month of January, which is observed as the annual canal closure. A relatively low volume of rainfall, i.e., 50 mm, was observed in July 2014, which caused the deficit to rise up to 30 mm in July, one of the hottest months, while the year 2015 received higher rainfall, especially in July, of up to 120 mm, minimizing the deficit to as low as 5 mm. However, the study area received 80 mm rainfall in July for the year 2016, and the deficit climbed up to 40 mm because of relatively low rainfall accompanied with a reduction in irrigation supplies due to maintenance work going on in LBDC [
8].
(iv) The consumptive use of irrigation supply and groundwater is vital to meet crop water needs, though consumptive use is still not able to meet crop demand as observed in the results. However, the rainfall pattern directly influences the quantitative behavior of the deficit. The research results endorsed the water budget published by the On Farm Water Management (OFWM) Agriculture Department that stated that groundwater and irrigation supply are not sufficient to meet crop demand in the LBDC, and rainfall plays an important role especially in the summer to reduce deficit. Moreover, the deficit values remained high in March and April for rabi crops and in August and September for kharif crops due to the high crop water demand. Water deficit was observed to be consistent towards the tail end of the main canal and secondary canals, which needs to be addressed to ensure equitable distribution of irrigation water resources.
There is a dire need to use modern technology of remote sensing and GIS to efficiently utilize irrigation water resources and avoid crop water stress by monitoring crop health and climate conditions simultaneously. A spatial decision support system (SPSS) needs to be developed for better management of canal water and alternate sources available for crop water use. The current study has laid the foundations for using independent and open-source technology to monitor crop water stress on a near real-time basis and guide decision makers to take informed decisions.