Assessment of Irrigation Performance in Large River Basins under Data Scarce Environment—A Case of Kabul River Basin, Afghanistan

The Kabul River basin (KRB) of Afghanistan, a lifeline of around 10 million people, has multiplicity of governance, management, and development-related challenges leading to inequity, inadequacy, and unreliability of irrigation water distribution. Prior to any uplifting intervention, there is a need to evaluate the performance of irrigation system on a long term basis to identify the existing bottlenecks. Although there are several indicators available for the performance evaluation of the irrigation schemes, we used the coefficient of variation (CV) of actual evapotranspiration (ETa) in space (basin, sub-basin, and provincial level), relative evapotranspiration (RET), and temporal CV of RET, to assess the equity, adequacy, and reliability of water distribution, respectively, from 2003 to 2013. The ETa was estimated through a surface energy balance system (SEBS) algorithm and the ETa estimates were validated using the advection aridity (AA) method with a R2 value of 0.81 and 0.77 at Nawabad and Sultanpur stations, respectively. The global land data assimilation system (GLDAS) and moderate-resolution imaging spectroradiometer (MODIS) products were used as main inputs to the SEBS. Results show that the mean seasonal sub-based RET values during summer (May–September) (0.37 ± 0.06) and winter (October–April) (0.40 ± 0.08) are below the target values (RET ≥ 0.75) during 2003–2013. The CV of the mean ETa, within sub-basins and provinces for the entire study period, has an equitable distribution of water from October–January (0.09 ± 0.04), whereas the highest inequity (0.24 ± 0.08) in water distribution is during early summer. The range of the CV of the mean ETa (0.04–0.06) on a monthly and seasonal basis shows the unreliability of water supplies in several provinces or sub-basins. The analysis of the temporal CV of mean RET highlights the unreliable water supplies across the entire basin. The maximum ETa during the study period was estimated for the Shamal sub-basin (552 ± 43 mm) while among the provinces, Kunar experienced the highest ETa (544 ± 39 mm). This study highlights the dire need for interventions to improve the irrigation performance in time and space. The proposed methodology can be used as a framework for monitoring and implementing water distribution plans in future.


Introduction
The Kabul River Basin (KRB) is strategically very important for its transboundary nature and around 50% flow contribution to the overall outflow of the basin. Irrigated agriculture in the KRB irrigation performance, remote sensing data-based analyses are extremely useful [11], with which one can go back in time and carry out comparative analyses with the present situation. This is usually required for irrigation system uplifting in terms of efficiency and productivity. It therefore provides a low-cost remote sensing-based approach for making a quick analysis of the performance of an irrigation system [11]. Owing to the successful application and implementation of the SEBS in different ecological zones of the world for ET a estimation, it has many benefits in the cases of large temporal and spatial scales. However, there are some limitations to this study too, which affected the overall quality of the results targeted at the performance evaluation of the irrigation system. One of these limitations is that that it did not simulate the irrigation application at farms, and no direct in-situ measurements of the irrigated water use were made. This was driven by the limited resources for sampling across the large basin as well as the prevailed insecurity conditions that didn't allow for extensive measurements in different irrigated lands. Beside this, the SEBS estimations were in reality a model (SEBS) to model (AA) comparisons where AA estimations were used to validate a SEBS estimated ET a . This is because this study lacks the ability to provide detailed information on the actual field evaporative fraction calculated by means of other methods (i.e., Energy Balance Bowen Ratio, scintillometry, and Eddy Covariance systems). Thus, in the absence of actual field measurements of ET a , the AA method could be used to validate the results of the ET a estimated by energy balance algorithms as performed by the previous studies [12,13].
Equity, reliability, and adequacy are the main challenges in water distribution at the system level for water managers, and not only in the KRB but also in other large irrigation schemes around the world [14,15]. In the current study, we assessed the equity, reliability, and adequacy of KRB at the system level (e.g., provinces and sub-basin scale). We used metrics that require only remote sensing data and have been documented well in other studies. For example, RET was used to assess the adequacy of water distribution at the water user's association level by using satellite remote sensing [13,15]. Similarly, Ahmad, Bastiaanssen and Usman, et al. [9] used the CV of actual evapotranspiration (ET a ) for assessing equity, which is the seasonal evaporative fraction for assessing the adequacy and temporal CV of the evaporative fraction for describing reliability. Thus, the objective of the current study is to assess the irrigation performance of the data-scarce KRB through a set of criteria (equity, adequacy and reliability) across the constituent sub-basins and provinces of the KRB for the period of 2003-2013. The coefficient of variation (CV) of the mean ET a , RET, and CV of RET were the metrics used to assess the equity, adequacy, and reliability, respectively. This study not only provides guidelines for the water managers in the region to optimize the operational and strategical performance of the KRB, but also the proposed set of indicators can be used as part of the monitoring and evaluation criteria for irrigation performance by relevant platforms in the future. The proposed set of indicators could also be helpful for the performance evaluation of the irrigations systems in data-scarce basins where ground-based measurements are missing or otherwise difficult to perform.

Description of the Study Site
The Kabul River Basin (KRB) is one of the five major river basins of Afghanistan having around 72,646 km 2 of land area. It spans a wide swath of land starting from the central highlands of the country at 6000 m above mean sea level down to the valleys in the east at 400 m above mean sea level ( Figure 1).
In 2013, the KRB received an annual precipitation of 327 mm at the downstream with a usual fringe effect of the Indian monsoon coming from the South Asian Himalayas, as well as around 418 mm at the upstream. The mean annual temperature at the central upstream and downstream locations were 13 • C and 23 • C, respectively ( Figure 2).  Around 9% of the total land area of the KRB is cropland [16]. The main sources of irrigation of the cropland in the KRB are canals, streams, springs, karezes (underground water channels which exploit unconfined aquifers in alluvial fans which are recharged by snowmelt), as well as groundwater wells. Due to the relatively intensive conventional canal networks developed in the eastern provinces, mainly Nangarhar, Laghman, and Kunar, irrigated agriculture contributes the most to food production in the KRB. There is additional potential for the enhancement of irrigated agriculture [17] provided with modifications in the irrigation infrastructure and economical investment [18]. The existing irrigation system in the KRB is made of conventional schemes usually developed, constructed, operated, and maintained by farmers according to traditional customs and practices with some exceptions of recent interventions by the Afghan government. Most of the farmers are unaware of the actual crop water requirements and scheduled irrigations for the crops that are being grown locally. As a consequence, the delivery of water to the field is based on the rule  In 2013, the KRB received an annual precipitation of 327 mm at the downstream with a usual fringe effect of the Indian monsoon coming from the South Asian Himalayas, as well as around 418 mm at the upstream. The mean annual temperature at the central upstream and downstream locations were 13 °C and 23 °C, respectively ( Figure 2). Around 9% of the total land area of the KRB is cropland [16]. The main sources of irrigation of the cropland in the KRB are canals, streams, springs, karezes (underground water channels which exploit unconfined aquifers in alluvial fans which are recharged by snowmelt), as well as groundwater wells. Due to the relatively intensive conventional canal networks developed in the eastern provinces, mainly Nangarhar, Laghman, and Kunar, irrigated agriculture contributes the most to food production in the KRB. There is additional potential for the enhancement of irrigated agriculture [17] provided with modifications in the irrigation infrastructure and economical investment [18]. The existing irrigation system in the KRB is made of conventional schemes usually developed, constructed, operated, and maintained by farmers according to traditional customs and practices with some exceptions of recent interventions by the Afghan government. Most of the farmers are unaware of the actual crop water requirements and scheduled irrigations for the crops that are being grown locally. As a consequence, the delivery of water to the field is based on the rule Around 9% of the total land area of the KRB is cropland [16]. The main sources of irrigation of the cropland in the KRB are canals, streams, springs, karezes (underground water channels which exploit unconfined aquifers in alluvial fans which are recharged by snowmelt), as well as groundwater wells. Due to the relatively intensive conventional canal networks developed in the eastern provinces, mainly Nangarhar, Laghman, and Kunar, irrigated agriculture contributes the most to food production in the KRB. There is additional potential for the enhancement of irrigated agriculture [17] provided with modifications in the irrigation infrastructure and economical investment [18]. The existing irrigation system in the KRB is made of conventional schemes usually developed, constructed, operated, and maintained by farmers according to traditional customs and practices with some exceptions of recent interventions by the Afghan government. Most of the farmers are unaware of the actual crop water requirements and scheduled irrigations for the crops that are being grown locally. As a consequence, the delivery of water to the field is based on the rule of maximizing the amount captured, leading to imbalances and water losses at different reaches along the canals with a potential yield and biomass loss.

Diagnosis of the Irrigation System Performance of the KRB
For the diagnosis of irrigation system performance, researchers use different indicators including overall consumed ratio [19], crop water deficit [20], relative water supply [21], relative evapotranspiration [14,15], delivery performance ratio, drainage ratio and depleted fraction [15], and equity and reliability [8,14], etc. For this study, equity, adequacy, and reliability were used as the primary criteria for the assessment of irrigation performance at administratively important spatial units (sub-basin and provincial) and strategically viable time steps (annual, monthly, and seasonal). The choice of these indicators is mainly due to being required by the mainly remote sensing based ET a data rather than detailed secondary data or field measurements, which were hindered by prevailed insecurity and political instability in the rural irrigated regions of the KRB. Moreover, the selection is also based on the consideration of quantifiability, cost effectiveness, and accuracy required [22].
Equity is an extremely important aspect of irrigation management, which deals with the spatial distribution of water from the system manager's point of view for large supply-based irrigation systems [23]. But in areas like KRB, suffering from multilateral issues (e.g., infrastructure damage, poor hydro-meteorological network, and water scarcity, etc.), equity in water consumption is more relevant at the system level and therefore can be computed from the remote sensing-based estimates of ET a . The CV of ET a was used as the key indicator to assess this criterion across sub-basins and provinces both in the summer and winter seasons, from 2003 to 2013. The CV used here is a measure of the relative variability which that calculated as the ratio of the standard deviation of the mean monthly ET a (2003-2013) to the mean monthly ET a (Equation (1)): The high CV of ET a , across the spatial constituent units of the KRB (i.e., sub-basins and provinces), would reveal high inequity in the irrigation system and vice versa. Adequacy is used to assess the reduction in ET a as well as to evaluate the sufficiency of water delivery to a known command area [21]. For the evaluation of the adequacy, the relative evapotranspiration (RET) was used to detect the areas with water shortages over the time period of 2003-2013. The RET is very important in relation to the monitoring of agricultural drought as well as crop yield forecasting. Relative evapotranspiration is known to be a reliable measure of plant available soil water and is proportional to plant growth [24]. The RET was calculated as a ratio of ET a to ET c as given in Equation (2): The crop evapotranspiration (ET c ) in this case is determined by the crop coefficient approach, as explained in the FAO irrigation and drainage paper 56 [25], and was calculated using by using Equation (3): where K c is the crop coefficient and ET 0 is the reference evapotranspiration calculated through the FAO Penman-Monteith method [25].
The higher values of RET in a specific sub-basin shows that this sub-basin is under lesser stress, and vice versa. Usually RET ≥ 0.75 is acceptable for irrigated agriculture, although it is constantly over time [26,27], while RET between 0.75 and 1.0 shows adequate water supply [15]. The value of RET = 0.65 is considered to be the critical value [28] and below that level (RET < 0.65) is considered to be poor performance of the irrigation system, and eventually leads to a loss to the crop yield and biomass. Reliability is the measure that describes the sufficiency of water availability for crops' consumption throughout the season. Therefore, the temporal CV of RET across the constituent provinces of the KRB was used for the assessment of the water reliability.

Actual Evapotranspiration Estimated by Remote Sensing
There are different algorithms for mapping the ET a , but their selection depends on the level of accuracy required, the geographical conditions of the targeted area, model limitations, etc. [29]. Based on the acceptable estimates of the Surface Energy Balance System (SEBS) algorithm [30] used at the Indus Basin and similar regions [12,13], it was also chosen for this study to estimate ET a for the KRB. The GLDAS and MODIS data were used as inputs to the SEBS. Furthermore, it is a well-established model and has been validated under the different agro-climatic conditions of the world. Former studies show that SEBS results, when validated against the observed eddy covariance (EC) technique [31], ground-measured ET a [32], and the ET a estimated with advection aridity (AA) in the Indus Basin [13], yielded an acceptable correlation. In this study, we used the AA method, due to its fewest ground data requirements, high reliability, and successful implementation in the same region (Indus Basin) [13].
Therefore, a methodological framework ( Figure 3) was established to estimate the ET a covering the study period (2003-2013) by the SEBS model for administratively important spatial units (basin, sub-basin, and provincial) and strategically viable time steps (annual, monthly, and seasonal). Main Input Data Characteristics for the SEBS The Global Land Data Assimilation System (GLDAS) is a unique uncoupled land surface modeling system that drives multiple models and integrates a large quantity of observed data purposed to ingest satellite and ground-based measured data. It runs globally with a spatial resolution of 0.25° with 3 h of step information [38]. The use of GLDAS datasets is very helpful while dealing with areas suffering from data scarcity or the absence of ground meteorological information [39,40]. The meteorological variables (Table 1)   The SEBS algorithm is a single-source model used for the estimation of atmospheric turbulent fluxes and the surface evaporative fraction derived from the remote sensing data. The SEBS algorithm employs meteorological and satellite spectral reflectance and radiance data for the estimation of the turbulent heat fluxes. It is based entirely on the rationale of the basic equation (Equation (4)) used for the computation of surface energy balance and is given below: λE is the turbulent latent heat flux (Wm −2 ) with λ being the latent heat of vaporization (Jkg −1 ), and E is evapotranspiration. It is worth mentioning that most of the energy balance models differ from each other in the way, they calculate the sensible heat flux. SEBS is basically a single-source algorithm, which uses either the Monin-Obukov Similarity (MOS) theory [33] for atmospheric surface layer (ASL) scaling or the Bulk Atmospheric Similarity (BAS) theory [34] used for atmospheric boundary layer scaling, depending on the height at which the measurements were taken in order to simulate the sensible heat flux used in Energy budget Equation (4). In this study, all measurements were conducted within the ASL. The MOS theory was considered for the allocation of H from the available energy through an iteration procedure in this study [30]. SEBS is a widely used and published model, thus readers are referred to previous studies [13,30,31,35] for a detailed description of the algorithm. But it is important to say that during cloud-free sky situations the evaporative fraction (Λ) (known to be conservative in a diurnal cycle) is considered to be representative of daily energy partitioning. Based on this assumption, an instantaneous evaporative fraction (Λ) during the satellite overpass time was used for the computation of the instantaneous ET at different timescales. The instantaneous ET a (mm s −1 ) was scaled up to daily values, ET a24 (mm/day), by estimating the averaged 24 h net radiation (R N24 ) [35,36], and for this reason we summed up the instantaneous ET a values over 24 h. This adaptation was based on the assumption that the evaporative fraction almost remains constant, though the sensible and latent heat fluxes may vary strongly during a day [37].

Main Input Data Characteristics for the SEBS
The Global Land Data Assimilation System (GLDAS) is a unique uncoupled land surface modeling system that drives multiple models and integrates a large quantity of observed data purposed to ingest satellite and ground-based measured data. It runs globally with a spatial resolution of 0.25 • with 3 h of step information [38]. The use of GLDAS datasets is very helpful while dealing with areas suffering from data scarcity or the absence of ground meteorological information [39,40]. The meteorological variables (Table 1) extracted from the Goddard Earth Sciences Data and Information Services Center (GES DISC-http://disc.sci.gsfc.nasa.gov/hydrology) for tiles H 23 V 5 , which covers the study area, were used in the SEBS for the ET a estimation in the data-scarce KRB. The MODIS data sets (Table 2) used in the SEBS were downloaded from the Land processes Distributed Active Archive Center (LP DAAC) of the United States Geological Survey (USGS) (https://lpdaac.usgs.gov/products/modis_products_table). The downloaded meteorological variables were interpolated in a linear way to match the MODIS's temporal resolutions over-pass time over the KRB. While using the MODIS re-projection tool, the downloaded data sets were re-sampled by using the nearest neighbor interpolation method. The MODIS land surface temperature data is produced daily (instantaneous), while the leaf area index (LAI) is an 8-day composite dataset. Because land surface status, defined by LAI and surface reflectance or surface albedo (MCD43B3), does not alter significantly over short periods, the 8-day interval is enough to portray the land surface properties [35]. The land cover map (MCD12Q1) was used in the analysis of the evaporative behavior of different land cover types in the KRB.

Evaluation of the SEBS Actual Evapotranspiration through Advection-Aridity Model (AA)
The advection-aridity model (AA model) is an energy balance model [41] used in this study to estimate the actual evapotranspiration from the meteorological data of two stations, located at the Nangarhar and Kunar provinces. Both the stations are located at places where there are no major variations in cropping pattern for a few kilometers from the stations. The spatial difference in estimated ET of both methods is small and negligible. The results of the AA model were used to evaluate the ET a estimates of SEBS model. The AA model was chosen because of its suitability under arid and semi-arid conditions [13,42,43]. The main benefit of the AA complementary method is that it does not require any information on plant canopy resistance, stomatal resistance vegetation properties, soil moisture, or other measures of aridity, because it depends mainly on meteorological parameters [41,44]. The details of the AA model have been elaborated in several studies under various geographic and climatic conditions [12,45,46], which states that the AA method was able to produce accurate estimations of the ET a when compared to the eddy covariance measurements or the water balance approach. The AA model [41] for regional evapotranspiration estimation, which is based on Bouchet's complementary relationship [47], and expresses the ET a as a combination of the wet environment (ET w ) and potential evapotranspiration (ET p ) (Equation (5)): where ET a is the actual evapotranspiration, ET w is the evapotranspiration under a wet surface, and ET p is the potential evapotranspiration.

Mann-Kendall Test for Monotonic Trend in Temperature
The trend in the meteorological time-series is complicated because of the skewness, persistence, and seasonal characteristics in the data. Different researchers [48,49] have used the non-parametric Mann-Kendall [50,51] and Sen's [52] slope estimator. Several researchers [53] prefer the use of the Mann-Kendall test because it provides the best alternative when the observed data is skewed either positively or negatively, or otherwise correlated cyclically or serially [54]. Therefore, to evaluate if there Remote Sens. 2018, 10, 972 9 of 21 has been a significant monotonic upward or downward trend of the variable (temperature) over time that could trigger possible changes in the ET a trend during 2003-2013, the Mann Kendall (MK) test was used, which is applied to the monthly inputs of the targeted variable. The monthly mean temperature (2003-2013) was used for the MK trend analysis test, and the MK test statistics are calculated on a monthly basis by using Equation (6), while the variance (MK test) is computed by using Equation (7). Besides this, the MK statistic and its variance are used to calculate the cumulative MK statistics and variance by using Equations (8) and (9), respectively, for a specific season/year. Finally, the cumulative MK statistic and variance is used to compute the Z statistic by the set of Equation (10): where in Equation (6) When each time step is shown by a single data point (with sample size > 10), then the variance of MK statistics is calculated through the following equation: where g i = number of groups with tied data in month i; t ip = number of tied data points in group p of month i; and N i = number of data points (over years) for month i. After calculating S i and VAR(S i ) for individual months, the annual sum was calculated by using equations C and D: The standardized Z MK statistics were calculated using the following equation which follows a standard normal distribution and can be related to a P value for testing the significance of the trend [55]: , i f S < 0 (10)

Validation of the ET a Estimated through SEBS with AA Model Estimates
Based on the validation results shown in Figure 4, there has been a sound fitness between the ET a estimated through SEBS and the ET a estimated through AA at the Kunar and Nangarhar provinces. The coefficients of determination thus obtained were R 2 = 0.81 and R 2 = 0.77, respectively (Figure 4). For the cool months of the year, the ET a values of the AA model are lower than those of SEBS. The reason behind is that the AA model uses a form of the Penman equation that does not work well for those periods for which the available energy (R n ) is negative or otherwise very close to zero. A similar result has been obtained from another study [56] whereby they estimated ET a through Complementary Relationship Areal Evapotranspiration (CRAE) and AA models, which resulted in lower values by using the AA model against those of CRAE in the cool months of winter. Another study in the Indus Basin [13] also shows that the ET a calculated through AA was lower in the cool months of winter (October to March) compared to ET a estimates through SEBS for the same period. The AA model yields lower value of ET a under high precipitation conditions [43], which goes in line with the result of this study.

Validation of the ETa Estimated through SEBS with AA Model Estimates
Based on the validation results shown in Figure 4, there has been a sound fitness between the ETa estimated through SEBS and the ETa estimated through AA at the Kunar and Nangarhar provinces. The coefficients of determination thus obtained were R 2 = 0.81 and R 2 = 0.77, respectively ( Figure 4). For the cool months of the year, the ETa values of the AA model are lower than those of SEBS. The reason behind is that the AA model uses a form of the Penman equation that does not work well for those periods for which the available energy (Rn) is negative or otherwise very close to zero. A similar result has been obtained from another study [56] whereby they estimated ETa through Complementary Relationship Areal Evapotranspiration (CRAE) and AA models, which resulted in lower values by using the AA model against those of CRAE in the cool months of winter. Another study in the Indus Basin [13] also shows that the ETa calculated through AA was lower in the cool months of winter (October to March) compared to ETa estimates through SEBS for the same period. The AA model yields lower value of ETa under high precipitation conditions [43], which goes in line with the result of this study.

Annual Distribution of Actual Evapotranspiration across the KRB, Constituent Sub-Basins and Provinces
There is around 9% increment in the mean annual ET a of KRB from 2003 to 2013, with around 2% of agricultural land cover expansion across the KRB [16]. The LULC analysis, as well as the secondary data [57], shows that from 2003 until 2013 there was increase in the cultivation of wheat and rice across the country, which are the main staple foods in Afghanistan. The minimum mean annual ET a was estimated for the Chak aw Logar sub-basin (i.e., 420 ± 21 mm), where the relatively higher standard deviation has been driven by the drought conditions in 2004 that caused a drop of 19% in the ET a in the Shamal sub-basin in comparison to the annual ET a in the preceding year (i.e., 2003). In comparison to 2003, Kunar sub-basin experienced the highest increase (13%) in the mean annual ET a , while Alingar experienced the least increase (4%). The outliers in the box-plot (Figure 5a) for the Shamal, Kunar, and Gomal subbasins show the effects of drought in 2004, which were relatively dominant in the mentioned sub-basins. In reference to 2003, the Nangarhar province experienced the highest increase in the ET a (15%) in 2013, followed by the Khost province (i.e., 12%), while the Paktya province experienced approximately an 8% reduction. The reason for the increase in the Nangarhar and Khost provinces could be attributed to infrastructural development leading to irrigated land area enhancement as well as the resultant increased crop rotation compared to a decade earlier, while the reduction of ET a in the Paktya province is perhaps due to the reduced water supply for irrigation purpose.
for the Shamal, Kunar, and Gomal subbasins show the effects of drought in 2004, which were relatively dominant in the mentioned sub-basins. In reference to 2003, the Nangarhar province experienced the highest increase in the ETa (15%) in 2013, followed by the Khost province (i.e., 12%), while the Paktya province experienced approximately an 8% reduction. The reason for the increase in the Nangarhar and Khost provinces could be attributed to infrastructural development leading to irrigated land area enhancement as well as the resultant increased crop rotation compared to a decade earlier, while the reduction of ETa in the Paktya province is perhaps due to the reduced water supply for irrigation purpose.

Monthly Distribution of Actual Evapotranspiration Across the KRB, Constituent Sub-Basins and Provinces
The mean annual ET a shows that the usual highest values of ET a across the KRB are in the months of May-July, which are around 70 ± 1.6 mm. The seasonal variability in the ET a in the summer season (May-September) across the KRB stays almost consistent with a seasonal mean ET a of 333 ± 19 mm. Overall, the least monthly mean ET a for the month of January and December across all the sub-basins was 9 ± 3 mm) and (7 ± 1 mm), respectively (Figure 6a,b).
The Panjshir province experienced the lowest mean ET a in the winter season during the study period (Figure 6c). Similarly the highest ET a in summer season among provinces was observed by Kunar and Nangarhar provinces due to their relatively easier access to river and canal water supplies (Figure 6d). Contrary to these estimations, another study [58] estimated 570 mm as the seasonal ET a (May-September) being an average of 3 years (2003-2005) for which the estimate of the current study for the very period is about 259 mm. The high level of ET a estimations in the Kabul province by another study [58] may be due to the tendency of SSEB's overestimation of ET a , both at the local and regional scales, probably due to rainfall contributions and abundant soil moisture that naturally supplement crop water needs [59].
Kunar and Nangarhar provinces due to their relatively easier access to river and canal water supplies (Figure 6d). Contrary to these estimations, another study [58] estimated 570 mm as the seasonal ETa (May-September) being an average of 3 years (2003-2005) for which the estimate of the current study for the very period is about 259 mm. The high level of ETa estimations in the Kabul province by another study [58] may be due to the tendency of SSEB's overestimation of ETa, both at the local and regional scales, probably due to rainfall contributions and abundant soil moisture that naturally supplement crop water needs [59]. A routine inter-seasonal comparison of ETa in summer and winter shows a higher range of ETa in summer 2013 compared to that in the winters of 2012-2013 ( Figure 7). The reasons for this include the peak irrigation demand, favorable meteorological conditions for crop growth, diversity and abundance of crops, and the vegetables and fruit orchards mostly in the summer season while winter is limited to fewer crops and is mostly dominated by wheat and barley, etc.  A routine inter-seasonal comparison of ET a in summer and winter shows a higher range of ET a in summer 2013 compared to that in the winters of 2012-2013 ( Figure 7). The reasons for this include the peak irrigation demand, favorable meteorological conditions for crop growth, diversity and abundance of crops, and the vegetables and fruit orchards mostly in the summer season while winter is limited to fewer crops and is mostly dominated by wheat and barley, etc.

Land Cover Based Variation of ET a Across the KRB
The land cover-based ET a results (Figure 5b) show a high ET a in the areas with maximum precipitation, which usually falls in the months of November to March, and has been used in the peak irrigation period with various frequencies. Among crops, wheat is the highly consumed and cultivated crop across the KRB. It is normally cultivated in rotation with maize and rice mostly in the downstream of the KRB, dominantly on irrigated lands, while in the central KRB it is cultivated with fruit orchards, in contrast to the downstream regions. Therefore, from 2003-2013 the mean annual ET a of wheat-maize, wheat-rice, and wheat alone across the KRB was 377 ± 49, 472 ± 27, and 352 ± 54 mm, respectively. For wheat, the results from the neighboring Uzbekistan showed the ET a values for wheat to be 397 mm [60]. statistical details have been explained in the text).
A routine inter-seasonal comparison of ETa in summer and winter shows a higher range of ETa in summer 2013 compared to that in the winters of 2012-2013 (Figure 7). The reasons for this include the peak irrigation demand, favorable meteorological conditions for crop growth, diversity and abundance of crops, and the vegetables and fruit orchards mostly in the summer season while winter is limited to fewer crops and is mostly dominated by wheat and barley, etc.

Evaluation of the Results of Mann-Kendall Test for Monotonic Trend in Temperature
The Mann-Kendall test has been used for the statistical analyses of trend detection in temperature data from 2003 to 2013 (11 years). It was observed that there is no significance trend found in the temperature time series data of the KRB that could possibly influence the monthly ET a values in a trending way. The resultant mean P-value for the entire study period was 0.53, greater than the significance level (α = 0.05) that lead us to accept the null hypothesis (H 0 ) ( Table 3).

Analysis of Spatial Equity
A recent study in the KRB [16] shows that the ratio of cropped to non-cropped area varies substantially among different provinces and sub-basins. Thus, comparing the ET a for total area by considering cropped and non-cropped areas between the provinces and sub-basins will be misleading. We, therefore, delineated the cropped area from the non-cropped area and hence used the ET a from the cropped area among different provinces. (Figure 8). Moreover, different cropping patterns in different sub units would result in different RET values and hence the comparison of adequacy will be misleading. However, the irrigation system of KRB is supply-based rather than demand-based. The supply-based irrigation system provides a thin layer of an equal amount of water to different sub units regardless of the cropping pattern [15]. Therefore, the limited and equitable water amounts to all sub units restricting the farmers to grow similar crops. The results of [16] also show that the cropping patterns in KRB do not vary substantially in different sub-basins and provinces. The results in Figure 9 show that, generally, in winter seasons the CV of mean the ETa in the months of October to January are the least (0.09 ± 0.04) for all the sub-basins, which shows that greater equity exists among the sub-basins during the colder months of the year where there is no irrigation required. The highest CV was assessed for Chak aw Logar (0.41) in the month of February, followed by the second highest CV in the Shamal (0.35) and Kabul (0.33) sub-basins during the peak irrigation period (e.g., in the month of April) (Figure 9a). The reason behind the highest CV in the month of February is the earlier season warm up and raised temperatures that melt down the snow and contribute to the irrigation period. The high CV values for the Kabul and Shamal sub-basins are perhaps due to the steady melting of the snow packs on the tops and contributions to the peak irrigation season, which is usually from February to May.
Moreover, this sub-basin hosts the most irrigated area compared to the remaining sub-basins of the KRB. During summer seasons, there is high inequity in the Kabul sub-basin compared to the others (Figure 9b). Similarly, the spatial variation of CV in the Panjshir province among the provincial units shows the highest CV for the months of October, November, and April ( Figure 9c). Generally, the CV for areas with larger irrigation areas was higher than those with more rainfed agriculture, which shows high inequity among the spatial units of the KRB (Figure 9c,d). The results in Figure 9 show that, generally, in winter seasons the CV of mean the ET a in the months of October to January are the least (0.09 ± 0.04) for all the sub-basins, which shows that greater equity exists among the sub-basins during the colder months of the year where there is no irrigation required. The highest CV was assessed for Chak aw Logar (0.41) in the month of February, followed by the second highest CV in the Shamal (0.35) and Kabul (0.33) sub-basins during the peak irrigation period (e.g., in the month of April) (Figure 9a). The reason behind the highest CV in the month of February is the earlier season warm up and raised temperatures that melt down the snow and contribute to the irrigation period. The high CV values for the Kabul and Shamal sub-basins are perhaps due to the steady melting of the snow packs on the tops and contributions to the peak irrigation season, which is usually from February to May. Moreover, this sub-basin hosts the most irrigated area compared to the remaining sub-basins of the KRB. During summer seasons, there is high inequity in the Kabul sub-basin compared to the others (Figure 9b). Similarly, the spatial variation of CV in the Panjshir province among the provincial units shows the highest CV for the months of October, November, and April ( Figure 9c). Generally, the CV for areas with larger irrigation areas was higher than those with more rainfed agriculture, which shows high inequity among the spatial units of the KRB (Figure 9c,d).

Analysis of Seasonal Adequacy
Throughout the KRB, the mean maximum RET for summer (May-September) and winter (October-April) was 0.84 and 0.73, respectively ( Figure 10). In the irrigated parts of the KRB, the RET, estimated both for summer and winter, was below the critical value (i.e., RET ≤ 0.65). Similarly, in areas where ET a and ET c are closer to each other (RET > 0.65), they are elevated at higher altitudes with either snow cover or otherwise forested. For winter season, the highest RET for Ghorband aw Panjshir and Chak aw Logar sub-basins was 0.50 and 0.48, respectively, while the lowest RET was in the Gomal sub-basin (i.e., 0.28) (Figure 11). The reason for the higher RET in the Ghorband, Panjshir, and Chak aw Logar sub-basins might be the lower temperatures, and also low irrigation that is a consequence of the early stages of crop development with low water demand.
Overall, the RET for the winter season is higher than that of in summer. For the summer season the highest RET values were experienced by the Kunar and Alingar sub-basins, which were 0.46 and 0.43, respectively. While for the remaining sub-basins the RET was in the range of almost 50% of the critical value of RET and include the following: Chak aw Logar (0.33), Ghorband aw Panjshir (0.35), Gomal (0.31), Kabul (0.34), and Shamal (0.37). The comparative higher RET in Kunar and Alingar could be attributed to the presence of the less irrigated areas compared to the rest because the Kabul sub-basin hosts the highest irrigated land area compared to the rest and therefore stays under stress during the peak irrigation season, which takes place from month of April onwards.
the highest RET values were experienced by the Kunar and Alingar sub-basins, which were 0.46 and 0.43, respectively. While for the remaining sub-basins the RET was in the range of almost 50% of the critical value of RET and include the following: Chak aw Logar (0.33), Ghorband aw Panjshir (0.35), Gomal (0.31), Kabul (0.34), and Shamal (0.37). The comparative higher RET in Kunar and Alingar could be attributed to the presence of the less irrigated areas compared to the rest because the Kabul sub-basin hosts the highest irrigated land area compared to the rest and therefore stays under stress during the peak irrigation season, which takes place from month of April onwards.

Analysis of Temporal Reliability
The temporal reliability of water supplies for irrigation purposes in the KRB was evaluated by using the CV of RET at all the constituent provinces of the KRB, which shows that the CV of RET in winter (October-April) is usually higher than in the winter season. The higher CV in winter is strongly driven by the higher CV of RET, especially in the months of February to April (winter) and May (summer) (Figure 12) because it is the time where irrigation starts based on the irrigation water demand and temperatures on the rise, and therefore the unreliability of water supplies are experienced for the crops grown across the basin. The

Analysis of Temporal Reliability
The temporal reliability of water supplies for irrigation purposes in the KRB was evaluated by using the CV of RET at all the constituent provinces of the KRB, which shows that the CV of RET in winter (October-April) is usually higher than in the winter season. The higher CV in winter is strongly driven by the higher CV of RET, especially in the months of February to April (winter) and May (summer) (Figure 12) because it is the time where irrigation starts based on the irrigation water demand and temperatures on the rise, and therefore the unreliability of water supplies are experienced for the crops grown across the basin.
The Similarly, in the summer season (May-September), throughout the study period, the fluctuation in the CV of RET has largely been affected by the higher CV in the month of May. It is perhaps due to the fact that maize and some local vegetables are being grown in this month where fields are irrigated to bring them to a field capacity level for creating a germination-friendly crop-soil-moisture atmosphere.
The CV of RET varied within a range of 43-60% throughout 2003-2013 across all the constituent provinces of the KRB, which is quite narrow compared to that of winter season. Likewise, in winter season, the range of CV of RET in the summer was high in the Panjshir province whereas the highest CV of RET was 0.60 in the year 2009, which is considered to be a wet year due to receiving more precipitation than during a normal year. The Panjshir province is located at the upstream heights of the KRB, which receives a larger share of the annual precipitation that is a result of melt that is flown downstream and contributes to the irrigation network in the KRB. After the Panjshir province, Nangarhar has the highest range of the CV of RET during 2003-2013, as it receives snow-melt from upstream in the form of streamflow and therefore is the largest irrigated land-holding province in the KRB, and thus enjoys crop rotation (i.e., usually two crops per year). It shows that the largest irrigated part of the KRB is suffering from unreliable water supplies both in summer and winter. Nuristan province has the narrowest range of the CV of RET (0.49-0.52), as it can be justified with its mountainous terrain with dominant forests and limited cropped area, which receives rainfall (other than winter precipitation in the form of snow) that contribute to its canals for irrigation purposes each year due to the fringe effects of monsoon from Himalayas.  Similarly, in the summer season (May-September), throughout the study period, the fluctuation in the CV of RET has largely been affected by the higher CV in the month of May. It is perhaps due to the fact that maize and some local vegetables are being grown in this month where fields are irrigated to bring them to a field capacity level for creating a germination-friendly crop-soil-moisture atmosphere.
The CV of RET varied within a range of 43-60% throughout 2003-2013 across all the constituent provinces of the KRB, which is quite narrow compared to that of winter season. Likewise, in winter season, the range of CV of RET in the summer was high in the Panjshir province whereas the highest CV of RET was 0.60 in the year 2009, which is considered to be a wet year due to receiving more precipitation than during a normal year. The Panjshir province is located at the upstream heights of the KRB, which receives a larger share of the annual precipitation that is a result of melt that is flown downstream and contributes to the irrigation network in the KRB. After the Panjshir province, Nangarhar has the highest range of the CV of RET during 2003-2013, as it receives snow-melt from upstream in the form of streamflow and therefore is the largest irrigated land-holding province in the KRB, and thus enjoys crop rotation (i.e., usually two crops per year). It shows that the largest In general, the higher range of the CV of RET shows that the supply-based irrigation system in the KRB delivers unreliable water supplies to the users, which require management interventions for raising water productivity across the KRB which hosts the highest population of the country compared to any other basin in Afghanistan. The different ranges of the CV of RET at different provinces is driven by various factors, e.g., geographic location (upstream receives the snowfall as the winter precipitation while downstream is the recipient of snow-melt), topographic relief and canal network existence, and certainly management interventions by the relevant department and ministry plays a role too.

Conclusions and Recommendations
Decades of political instability and conflicts caused infrastructural damage that ultimately led to data scarcity on land and water resources in Afghanistan. The agriculture sector consumes around 98% of the total water withdrawn for irrigation purposes. Having consumed such a high amount of water for agriculture, the country's dependency on foreign food imports highlights the inefficiency of the irrigation system as well as the responsible institutions. In order to sustainably manage these water resources, it is extremely important to analyze the long-term operational and strategical performance of the irrigation system. To assess the irrigation performance for around a decade (2003-2013), we identified those sets of indicators that requires minimum secondary data, which is very effective in data-scarce basins like the KRB. The irrigation performance assessment in the KRB shows that the irrigation system in the KRB is experiencing inequity, inadequacy, and unreliability, especially during the peak demand period that results in inefficient irrigated agriculture and poor crop-water productivity. It therefore leads to larger voids between water demand and supply, and ultimately leads to failure to meet growing food demand across the country. Besides the water losses that occur during conveyance and application, it affects other water consuming sectors too that are deemed to have increased water needs in the case of industrial and municipal developments in the near future. The relevant platforms need to strategically invest in the development of irrigation infrastructure (e.g., canal lining, water storage structures, etc.) as well as the capacity development of the local manpower at the Ministries of Agriculture, Irrigation, and Livestock, as well as the Ministry of Energy and Water to safeguard the irrigation quota for the dominant crops, especially in the peak demand period. As shown in the study, the relative evapotranspiration across all the sub-basins shows the entire KRB, hosting around one-third of the country's population, to be under water stress, which needs to be addressed and managed in order to meet the growing water and food demand in the country. The ET a , derived under this research could also be used in comparisons to future studies over the effects/impacts of climate change on ET a . Furthermore, the results of this study will considerably contribute to the investment plans addressing water resources management issues in the Kabul River Basin.
Author Contributions: The correspondent author of this research paper, F.A., wrote the manuscript and analyzed the data sets used in this study; U.K.A., B.T. and U.W.L. reviewed, edited and provided valuable suggestions and inputs for the final manuscript.
Funding: This research received no external funding other than what is acknowledged below.