Evaluating the Spatio-Temporal Distribution of Irrigation Water Components for Water Resources Management Using Geo-Informatics Approach

: Spatio-temporal distribution of irrigation water components was evaluated at the canal command area in Indus Basin Irrigation System (IBIS) by using a remote sensing-based geo-informatics approach. Satellite-derived MODIS product-based Surface Energy Balance Algorithm for Land (SEBAL) was used for the estimation of the actual evapotranspiration (ETa). The ground data-based advection aridity method (AA) was used to calibrate and validate the model. Statistical analysis of the SEBAL based ETa and AA shows the mean values of 87.1 mm and 47.9 mm during Kharif season (May–November) and 100 mm and 77 mm during the Rabi Season (December–April). Mean NSEs of 0.72 and 0.85 and RMSEs 34.9 and 5.76 during the Kharif and the Rabi seasons were observed for ETa and AA, respectively. Rainfall data were calibrated with the point observatory data of the metrological stations. The average annual ETa was found 899 mm for deﬁned four cropping years (2011–2012 to 2014–2015) with the minimum average value of 63.3 mm in January and the maximum average value of 110.6 mm in August. Average of the sum of net canal water use (NCWU) and rainfall during the study period of four years was 548 mm (36% of ETa). Seasonal analysis revealed 39% and 61% of groundwater extraction proportion during Rabi and Kharif seasons, dependent upon the occurrence of rainfall and crop phenology. Overall, the results provide insight into the interrelationships between key water resources management components and the variation of these through time, offering information to improve the strategic planning and management of available water resources in this region.


Introduction
Management of water resources to feed the swift increasing population is the biggest challenge of the 21st century. Ref. [1] found an exponential rise in population in South Asian countries, including Pakistan. Pakistan is one of the ten largest populous countries in the world. The agriculture-based economy of Pakistan is dependent on irrigation waters supplied by the Indus River and its tributaries. The largest irrigation system in the world is

SEBAL for the Estimation of Actual Evapotranspiration
Surface Energy Balance Algorithm for Land (SEBAL) is widely used for the estimation of the actual evapotranspiration in hydrological studies. The input used in this energy balance was the satellite-derived remote sensing imagery. SEBAL was successfully used in the water balance studies at different temporal and spatial scales [16,21,[28][29][30] researched in Pakistan and provided the theory of the SEBAL algorithm. The methodology was successfully applied in the IBIS by [31]. The residual of the energy budget was actual evapotranspiration, which is demonstrated below (Equation (1)). Canal water is the main source of irrigation in the LCC system along with a significant proportion of groundwater. The lower Chenab canal system is distributed into seven main distributaries. All these distributaries are monitored by the Program Monitoring and Implementation Unit (PMIU) of the Punjab Irrigation Department. Irrigation efficiency is a product of field and irrigation network efficiency. Irrigation efficiency was incorporated from the study of [6] that defines the field efficiency of 75% and the irrigation network efficiency of 48%. According to the aforementioned definition of irrigation efficiency, 36% irrigation efficiency of the system was obtained.

SEBAL for the Estimation of Actual Evapotranspiration
Surface Energy Balance Algorithm for Land (SEBAL) is widely used for the estimation of the actual evapotranspiration in hydrological studies. The input used in this energy balance was the satellite-derived remote sensing imagery. SEBAL was successfully used in the water balance studies at different temporal and spatial scales [16,21,[28][29][30] researched in Pakistan and provided the theory of the SEBAL algorithm. The methodology was successfully applied in the IBIS by [31]. The residual of the energy budget was actual evapotranspiration, which is demonstrated below (Equation (1)).
where, R n , G o , H, and λE represent net radiation (W m −2 ), soil heat flux (W m −2 ), sensible heat flux (W m −2 ), and latent heat flux (W m −2 ), respectively. The concept of evaporative fraction was introduced by [16] and, based on this concept, latent heat flux can be represented (Equation (2)) by considering net available energy (R n − G o ) and evaporative fraction: where, ∆ = evaporative fraction and was formulated as: Net available energy can be calculated from instant, daily, or monthly timescales. For the calculation of the daily ET, soil heat flux can be ignored at the time scale of one day, and net radiation (R n ) is considered as the net available energy. Daily ET can be calculated as Equation (4): where, R n24 , λ, and ρ w describe the 24 h averaged net radiation, latent heat of vaporization, and water density, respectively. Ref. [32] described the important variables for ET estimation, including land cover, land surface albedo, and land surface temperature. Remote sensing estimated the variables possible using different satellite sensors [30]. MODIS product was used due to its superior spatiotemporal resolution and spectral bands. The MODIS standard products can be downloaded for free from the USGS website https://lpdaac.usgs.gov/get_data/data_pool (accessed on 3 March 2021). Table 1 describes the detail of the product used in this study.

Advection Aridity Approach
The advection aridity (AA) model was successfully used by many researchers under different climatic conditions. With the water balance and the eddy covariance methods [33][34][35] applied, the advection aridity method at the basin scale showed good agreement with SEBS.
Ref. [36] applied the AA model at a regional scale and found good agreement between SEBAL and AA models. Ref. [34] applied three different methods for mapping the ETa in China and found the AA model most accurate. Firstly, Ref. [37] proposed the AA model using [38] for the wet surface ET as formulated in Equation (5). For the potential evapotranspiration, Ref. [39]'s equation was used as given below (Equation (6)), and an empirical wind function was used, as given below in Equation (7).
where, ∝ e = 1.26 is the Priestly and Taylor coefficient, and variables (R n − G), ∆, γ, U 2 , e s and e a represents the available energy at the surface, slope of saturated water vapor pressure curve at the current time (kP • C −1 ), psychometric constant (kP • C −1 ), mean wind speed above the 2 m of the ground surface, saturated vapor pressure and actual vapor pressure respectively. Ref. [37] combined the potential ET (Equation (5)) and the wet surface ET (Equation (6)) following the relationship given by [40], as shown below (Equation (9)): In this study, due to the unavailability of the directly measured ET, the AA model was used for the evaluation of SEBAL estimated ETa. The climatic data used for the AA model were collected from the University of Agriculture, Faisalabad, and the Pakistan Meteorological Department.

Satellite Derived Rainfall
The Pakistan Meteorological Department (PMD) is responsible for the measurement of meteorological data at a regional scale. There are 97 total observatories in Pakistan (PMD online source). However, there are only two observatory stations in the study area. GIS-based interpolation was not found appropriate for the low-resolution dataset. Therefore, satellite-based raster data were acquired with a spatial resolution of 25 km along with monthly temporal resolution using Tropical Rainfall Measurement Mission (TRMM). Low-resolution rain gauge availability creates a limitation for the calibration of TRMM data. Ref. [41] provided the solution under this limited environment by averaging the measurement over time. Therefore, 3B43 product was accumulated into seasonal basis using the cropping pattern. The Kharif season includes mid-April to the middle of November. Following the experience of [36], local calibration and validation of this data were performed with point observatory data. Monthly rainfall data were summed up to obtain seasonal rainfall for all Kharif and Rabi seasons. Further, GIS was used to extract the data from each pixel, and the detailed data were interpolated using kriging to create the high-resolution map.

Groundwater Irrigation
Net groundwater irrigation at the spatial scale of distributaries and two temporal scales of LCC were estimated using geo-informatics techniques. This approach provides the advances above conservative and indirect approaches [22]. Satellite-derived actual evapotranspiration was used as a key component for the establishment of a water balance approach in the unsaturated zone. This approach was successfully applied by [16] (Equation (10)).
NGW I = ET a − P e − NCW I where, NCWI, P e , and NGWi describe net canal water irrigation, effective rainfall, and net groundwater irrigation, respectively. This approach provided the net groundwater irrigation that was the water used by the plants or evaporated from the area under consideration, while the pumping from the aquifer was more than the net groundwater irrigation due to the inefficient delivery system of irrigation. Therefore, groundwater abstraction was estimated using the field application efficiency. In this study, the results of [6] were used that showed field application efficiency of 75% and network efficiency of 90%. This resulted in 68% irrigation efficiency of the system. The framework of the study is shown in Figure 2.

Canal Water Availability
Regarding estimation of canal water availability, the maximum gross canal water (GCW) irrigation depth of 343 mm was found during the 2014-2015 cropping year, while in the cropping year of 2011-2012, the GCW was found at a minimum with the value of 305 mm for the LCC (Table 2). Water supply during Kharif season was 32% more than Rabi season due to canal closure during Rabi. Every year, for the removal of accumulated silt in the irrigation system, canals remain closed during December and January [42] and sometimes for the month of February. During the spatial analysis, the gross canal water distribution was found to be maximum at the Upper Gugera, while the minimum was the Sagar distributaries ( Figure 3). The distribution of gross canal water from the LCC system was minimal during the month of February and maximal during the months of July and August. This lower canal water availability in the month of February was due to the closure of the canal for the annual removal of silts from the irrigation system. The maximal levels in July and August were due to the maximum surface water availability in the Indus basin due to heavy snow melting and monsoon rainfall. Net canal water distribution (mm) in LCC from 2011-2012 to 2014-2015 is shown in Table 3.

Canal Water Availability
Regarding estimation of canal water availability, the maximum gross canal water (GCW) irrigation depth of 343 mm was found during the 2014-2015 cropping year, while in the cropping year of 2011-2012, the GCW was found at a minimum with the value of 305 mm for the LCC (Table 2). Water supply during Kharif season was 32% more than Rabi season due to canal closure during Rabi. Every year, for the removal of accumulated silt in the irrigation system, canals remain closed during December and January [42] and sometimes for the month of February. During the spatial analysis, the gross canal water distribution was found to be maximum at the Upper Gugera, while the minimum was the Sagar distributaries ( Figure 3). The distribution of gross canal water from the LCC system was minimal during the month of February and maximal during the months of July and August. This lower canal water availability in the month of February was due to the closure of the canal for the annual removal of silts from the irrigation system. The maximal levels in July and August were due to the maximum surface water availability in the Indus basin due to heavy snow melting and monsoon rainfall. Net canal water distribution (mm) in LCC from 2011-2012 to 2014-2015 is shown in Table 3.

Calibration and Validation of SEBAL with Advection Aridity Method
None of the previous researchers deny the validation of the SEBAL model. Validation was done with the available data source of ET measurements such as Bowen ration energy balance [43], meteorological data [44], advection aridity method [36], and eddy covariance [45]. Ref. [34] found the AA to be the most accurate method after applying three different techniques for mapping the ETa at a regional scale in China. Validation of SEBAL is important for further use in hydrological studies. The calibration process was carried out by using SEBAL pixels and meteorological data available at the nearby climate monitoring station and is shown in Figure 4. SEBAL was successfully validated with the AA model by [36] at the regional scale studies in the Indus Basin Irrigation System. The results obtained were satisfactory for Rabi and Kharif seasons. Nash suit cliff efficiency was found to be higher in Rabi season as compared to Kharif season. Similarly, low bias values were found in the Rabi season relative to the Kharif season. The results are supported by the outcomes of [21]. Similarly, [34,46] argued that AA models work well during cold seasons. A significant difference in the ET was found for both seasons. Ref. [36] argued that this difference was due to the difference in climatic condition, change in crop type, and water availability. The variation in the ET was found from the upper reaches to the lower reaches of the LCC system with the same crop type. For the Kharif season, the difference was higher than the Rabi season. This variation was 9.63 mm for the rice crop, 9.71 mm for cotton, and 10.2 mm for sugarcane. For the Rabi crop, wheat variation was 1.03%, and sugarcane variation was 1.05%. Usman et al. (2014) argued that this variation increases if the analysis is performed at a finer resolution. Figure 4 represents the calibration and the validation of SEBAL with AA. Statistical analysis of seasonal and annual ETa (mm) of SEBAL and AA are explained in Table 4. of the LCC system with the same crop type. For the Kharif season, the difference was higher than the Rabi season. This variation was 9.63 mm for the rice crop, 9.71 mm for cotton, and 10.2 mm for sugarcane. For the Rabi crop, wheat variation was 1.03%, and sugarcane variation was 1.05%. Usman et al. (2014) argued that this variation increases if the analysis is performed at a finer resolution. Figure 4 represents the calibration and the validation of SEBAL with AA. Statistical analysis of seasonal and annual ETa (mm) of SEBAL and AA are explained in Table 4.

Spatially Distributed ETa
Seasonal analysis of Eta revealed 47% more ETa during the Kharif season as compared to the Rabi season, while the variation in the annual scale analysis showed less variation (Table 5). High ETa during Kharif was due to more surface water availability in the irrigation system along with the growing period of the high water demanding crops such as cotton and rice as compared to wheat. The month-based analysis revealed the variation in ETa with the minimum average value of 23 mm in January and a maximum average value of 123 mm in August and a four year average annual value of 897 mm during the study period (2010-2011 to 2014-2015). In the case of temporal analysis, May to August showed the peak rates of ETa, while, in the case of spatial analysis, peak rates were observed in rice-cropped areas. Similarly, December to January showed the lowest ETa rates due to less water supply along with decreased crop water requirements [42]. During the design of the LCC irrigation system, the cropping intensity was considered 75% [47]. Though, it increased more than double its design due to the green revolution of the 1960s. Despite this significant increase in demand due to high cropping intensity, the deliveries from the canal water remain stagnant, which causes a big gap between canal water availability and its demand and leads towards the severe scarcity in the region.  Figures 5 and 6 show the four years average of Rabi and Kharif season ETas, respectively, at the LCC system. The ETa map shows the significant variation in the actual evapotranspiration; the actual evapotranspiration was high at the head reaches of the LCC system. This higher ETa was due to more rice cultivation at the head reaches than the tail of the irrigation system, while, in the Rabi season, ETa was found higher from the center to the lower part of the system. This higher ET was due to the cultivation of more Rabi fodder crops than Kharif crops, especially the cultivation of the berseem crop.   . There was no significant temporal variation in the ETa between different CAs. In the case of Sagar CCA, the average ETa was 903 mm with a standard deviation of 12.5 mm. Mian Ali canal command showed the maximum standard deviation of 20 mm and Jhang canal command showed the lowest standard deviation of 12 mm. Due to fixed water allowance, stagnant resources, and locally fixed cropping zone, there were no changes in supply of irrigation water, practices of irrigation method, and or cropping patterns. Rakh and Mian Ali canal command showed low average ETa values as compared to the CCAs in the system. Likewise, Mian Ali and Rakh canal command areas had poor quality groundwater. Maximum ET at the Sagar canal commands was due to its location upstream of the irrigation system and cultivation of rice crop.

Spatially Distributed Rainfall
The existing studies describe over-and underestimations of rainfall from the TRMM relative to the changes in topography. Ref. [20] described, in detail, over-and underestimations of rainfall from the TRMM in the Indus basin. Similarly, [48] described the similar trend of rainfall over-and underestimations in Bangladesh. Further, [49] described that that satellite overestimated the rainfall in areas with the rainfall <400 mm month −1 and underestimated in the areas of rainfall >400 mm month −1 . These results were further strengthened by the findings of [48]. They concluded that the TRMM has limitations in accurate rainfall detection at low or high rainfall rates. Thus, in the study area, TRMM overestimated the rainfall that was calibrated with the point observatory data. The spatial distributions of rainfall for Rabi and Kharif seasons are shown in Figures 8 and 9, respec-

Spatially Distributed Rainfall
The existing studies describe over-and underestimations of rainfall from the TRMM relative to the changes in topography. Ref. [20] described, in detail, over-and underestimations of rainfall from the TRMM in the Indus basin. Similarly, [48] described the similar trend of rainfall over-and underestimations in Bangladesh. Further, [49] described that that satellite overestimated the rainfall in areas with the rainfall <400 mm month −1 and underestimated in the areas of rainfall >400 mm month −1 . These results were further strengthened by the findings of [48]. They concluded that the TRMM has limitations in accurate rainfall detection at low or high rainfall rates. Thus, in the study area, TRMM overestimated the rainfall that was calibrated with the point observatory data. The spatial distributions of rainfall for Rabi and Kharif seasons are shown in Figures 8 and 9, respectively. These figures clearly show the decrease in rainfall towards the tail of the irrigation system. Tables 6 and 7 describe the temporal analysis of gross and effective rainfalls. Table 6 shows the temporal rainfall distribution. The maximum mean monthly of four years of rainfall was found at a maximum of 95.6 mm in September and a minimum of 2.6 mm in November. The average annual rainfall for the four years was 343.5 mm.    Figure 10 shows that there were significant temporal variations in the effective rainfall at each CCA, while there was no significant average annual rainfall between different CCAs, as shown in the figure. The maximum average four years effective rainfall was 302 mm at Sagar, and the minimum was 258 mm at the Burala branch.    Figure 10 shows that there were significant temporal variations in the effective rainfall at each CCA, while there was no significant average annual rainfall between different CCAs, as shown in the figure. The maximum average four years effective rainfall was 302 mm at Sagar, and the minimum was 258 mm at the Burala branch.   Rainfall results are shown in Table 7. The annual effective rainfall variation ranged between 232.6 mm year −1 and 350.8 mm year −1 , with an average of 274 mm year −1 during 2010-2011 to 2014-2015. Variation in the mean monthly analysis for effective rainfall showed variations of 2.1 mm (minimum) in November and 76.5 mm (maximum) in September. The rainfall that occurred during the months of monsoon (June-September) was more than 55 mm month −1 . The analysis of seasonal rainfall occurrence found averages of 23% during Rabi and 77% during Kharif for total rainfall amount. Figure 10 shows the effective rainfall in different CCAs during the last four years (2011-2012 to 2014-2015). The maximum four years average annual effective rainfall in the LCC irrigation scheme was 302 mm at the Sagar. Figure 10 shows that there were significant temporal variations in the effective rainfall at each CCA, while there was no significant average annual rainfall between different CCAs, as shown in the figure. The maximum average four years effective rainfall was 302 mm at Sagar, and the minimum was 258 mm at the Burala branch.

Groundwater Irrigation
In arid and semi-arid areas, groundwater is a complementary source of irrigation under high cropping intensity. Groundwater contribution is significant in the IBIS. Ref. [16] described 42% of the groundwater pumping at the head of the irrigation system. Groundwater irrigation results showed that groundwater has a significant share to fulfill the irrigation water demand. The groundwater irrigation map obtained from the water balance described the significant groundwater contribution at the head of the canal system. Irrigation water supply from the system was aggregated to monthly, seasonal, and annual irrigation water volumes. Based on spatial mapping, four years' average net groundwater irrigation range was found from 195 mm to 88 mm during Rabi season and 314 mm to 99 mm during the Kharif season. The four years average of gross groundwater abstraction (GWA) was a maximum of 750 mm at the upstream areas of the irrigation system. The minimum GWA was 275 mm at the tail of the LCC system. This described more use of canal water at the head of the irrigation system. Four years average share of the net groundwater irrigation (NGWI) was 283 mm, which represented 44% of the average annual supply of canal water (Figures 11 and 12). The contribution of the average annual groundwater in the ETa was 42%. The monthly maximum average NGWI was during the month of May, whereas the minimum was during the month of November.

Groundwater Irrigation
In arid and semi-arid areas, groundwater is a complementary source of irrigatio under high cropping intensity. Groundwater contribution is significant in the IBIS. Re [16] described 42% of the groundwater pumping at the head of the irrigation system Groundwater irrigation results showed that groundwater has a significant share to fulfi the irrigation water demand. The groundwater irrigation map obtained from the wate balance described the significant groundwater contribution at the head of the canal sy tem. Irrigation water supply from the system was aggregated to monthly, seasonal, an annual irrigation water volumes. Based on spatial mapping, four years' average n groundwater irrigation range was found from 195 mm to 88 mm during Rabi season an 314 mm to 99 mm during the Kharif season. The four years average of gross groundwate abstraction (GWA) was a maximum of 750 mm at the upstream areas of the irrigatio system. The minimum GWA was 275 mm at the tail of the LCC system. This describe more use of canal water at the head of the irrigation system. Four years average share o the net groundwater irrigation (NGWI) was 283 mm, which represented 44% of the ave age annual supply of canal water (Figures 11 and 12). The contribution of the averag annual groundwater in the ETa was 42%. The monthly maximum average NGWI wa during the month of May, whereas the minimum was during the month of November.    The remaining crop water requirement is highly dependent on groundwater use, as the sum of the rainfall plus the irrigation water supply is less than the demand. The sum of both surface water sources (net canal plus rainfall) four years average in the study period was found to be 548 mm (36% of ETa). Therefore, farmers depend on the use of groundwater to fulfill the irrigation water need. The average of the study period showed 739.6 mm of groundwater abstraction, while the annually based variation ranged between 632 mm and 780 mm. Seasonal analysis revealed 39% and 61% of groundwater abstraction proportions during Rabi and Kharif seasons, respectively ( Table 8). The fluctuations in four cropping years for monthly groundwater abstraction ranged between 28.7 mm and 120.3 mm. This variation was high in the 2011-2012 to 2012-2013 cropping year (0 mm to The remaining crop water requirement is highly dependent on groundwater use, as the sum of the rainfall plus the irrigation water supply is less than the demand. The sum of both surface water sources (net canal plus rainfall) four years average in the study period was found to be 548 mm (36% of ETa). Therefore, farmers depend on the use of groundwater to fulfill the irrigation water need. The average of the study period showed 739.6 mm of groundwater abstraction, while the annually based variation ranged between 632 mm and 780 mm. Seasonal analysis revealed 39% and 61% of groundwater abstraction proportions during Rabi and Kharif seasons, respectively ( Table 8). The fluctuations in four cropping years for monthly groundwater abstraction ranged between 28.7 mm and 120.3 mm. This variation was high in the 2011-2012 to 2012-2013 cropping year (0 mm to 148.7 mm); this variation is dependent upon the occurrence of rainfall and crop phenology. This lowest abstraction in September could be due to high monsoon rainfall and to less water demand. Net groundwater irrigation estimated after incorporating the efficiencies was 503 mm year −1 ( Table 9) on average for the four cropping years. Unlike surface water use from canals, the maximum groundwater use in the month of May was due to the cultivation of rice as well as insufficient rainfall and canal water availability.   Figure 13 shows the gross groundwater irrigation (GGWA) in LCC system CCAs during the last four years (2010-2011 to 2014-2015). The four years average annual GGWA was found to be 526 mm for the whole LCC system. There was significant temporal variation in GGWA in the CCAs of the LCC system ( Figure 13). The four years average gross groundwater abstraction was 887 mm at the Sagar, and the minimum was at the Upper Gugera with zero groundwater abstraction for this study period. This was due to the maximum canal water supply to the Upper Gugera. The maximum abstraction was in the Sagar CCA, while the minimum abstraction was in the Upper Gugera CCA. Higher GGWA in the Sagar CCAS was due to the intensive cultivation of rice crop.

Discussion
One of the largest contiguous irrigation systems in the world is only covered with 96 metrological stations, and only three metrological stations cover the Lower Chenab Canal System. Effective and reliable estimation of actual evapotranspiration is not possible through the empirical models for the very fine assessment of irrigation water management components. SEBAL is the most reliable and effective technique applied in different irrigation systems for the estimation of actual evapotranspiration. SEBAL was successfully applied in the different regions of the world for the estimation of the actual evapotranspiration [16,44,50,51]. The gross canal water availability in the irrigation system described the quantity of water available from the irrigation sources, while the net canal water availability described the actual amount of water available in the field after the exclusion of the system losses. Irrigation system efficiency was about 45% from the point of water diversion from the river to the point of application in the field. Estimation of effective rainfall described the amount of water available as a supplementary source of irrigation and re-

Discussion
One of the largest contiguous irrigation systems in the world is only covered with 96 metrological stations, and only three metrological stations cover the Lower Chenab Canal System. Effective and reliable estimation of actual evapotranspiration is not possible through the empirical models for the very fine assessment of irrigation water management components. SEBAL is the most reliable and effective technique applied in different irrigation systems for the estimation of actual evapotranspiration. SEBAL was successfully applied in the different regions of the world for the estimation of the actual evapotranspiration [16,44,50,51]. The gross canal water availability in the irrigation system described the quantity of water available from the irrigation sources, while the net canal water availability described the actual amount of water available in the field after the exclusion of the system losses. Irrigation system efficiency was about 45% from the point of water diversion from the river to the point of application in the field. Estimation of effective rainfall described the amount of water available as a supplementary source of irrigation and reduced the demand for irrigation. The gross groundwater recharge was the amount of water that percolated due to the excessive irrigation from the canal system and the rainfall. This study was very comprehensive about the irrigation water management components and can help the policymakers effectively utilize the irrigation sources by reducing losses and improving effective irrigation scheduling by raising the flexibility. Remotely sensed data are acquired instantaneously and can only provide the instantaneous two-dimensional spatial distribution of land surface variables such as surface albedo, surface vegetation fraction, surface temperature, surface net radiation, soil moisture, etc., which are indispensable variables to know for remote sensing estimates of land surface ET. The methodology developed for the evaluation of irrigation water management components cannot evaluate short-term dynamics. The satellite data retrieved from the sources had a temporal resolution of 8 days and considered the cumulative effect; the daily change in the crop phenology was not considered. The method was found suitable and was widely applied for irrigation system management to develop crop combination maps, rainfall mapping, actual evapotranspiration mapping, and groundwater mapping where data scarcity is the main concern [20,25,44,[52][53][54]. The method is not limited to the model framework chosen here; it also has application in the dry areas irrigation systems such as India, Egypt, and African countries. The current methods cover the issue of data scarcity and provide a suitable solution for mapping the irrigation water management components at a very fine resolution. The key parameters involved in this study were the retrieval of the MODIS data with 8 days temporal and 250 m spatial resolution as well as the daily discharge data from the irrigation system and the retrieval of rainfall at very fine temporal and spatial resolution of TRMM. The study provided a very fine estimation of the irrigation water management components. Rainfall was estimated based on the study performed by the [20] at basin scale. Actual evapotranspiration was estimated using SEBAL [16] for the equity assessment at the very fine spatial scale of canal command area.

Conclusions
Four year average annual based actual evapotranspiration of 899 mm described the actual water requirement for the crop production. The water availability from the rainfall and the canal irrigation system during the four year study period only was 548 mm, which was only 36% of actual evapotranspiration. The contributions of the most flexible water resources of groundwater were 39% and 61% during Rabi and Kharif seasons. The most important resources of surface water supply variation at different scales were assessed for effective management. The spatial analysis of the system at canal command areas of the distributaries revealed the highest canal water availability at the Upper Gogera canal command and the lowest at the Sagar canal command. Therefore, groundwater extraction was found at its maximum at Sagar command due to less supply and cultivation of the high irrigation demanding crop (rice). Similarly, Jhang canal command, Burala canal command, and lower Goger canal command showed significant groundwater irrigation. Based on the study period, this study focused seriously on water allocation plans to avoid water logging in some areas and higher chances of secondary salinization in the major portion of the study area. It directed the water managers for the effective management of canal water resources to sustain soil and water productivity of the system for contribution in the achievement of sustainable development goals, as the rainfall occurrence is the uncontrolled resource, and groundwater is the most flexible water resource of the system.