Evapotranspiration in the Nile Basin: Identifying Dynamics and Drivers, 2002–2011

Analysis of the relationship between evapotranspiration (ET) and its natural and anthropogenic drivers is critical in water-limited basins such as the Nile. The spatiotemporal relationships of ET with rainfall and vegetation dynamics in the Nile Basin during 2002–2011 were analyzed using satellite-derived data. Non-parametric statistics were used to quantify ET-rainfall interactions and trends across land cover types and subbasins. We found that 65% of the study area (2.5 million km 2) showed significant (p < 0.05) positive correlations between monthly ET and rainfall, whereas 7% showed significant negative correlations. As expected, positive ET-rainfall correlations were observed over natural vegetation, mixed croplands/natural vegetation, and croplands, with a few subbasin-specific exceptions. In particular, irrigated croplands, wetlands and some forests exhibited negative correlations. Trend tests revealed spatial clusters of statistically significant trends in ET (6% of study area was negative; 12% positive), vegetation greenness (24% negative; 12% positive) and rainfall (11% negative; 1% positive) during 2002–2011. The Nile Delta, Ethiopian highlands and central Uganda regions showed decline in ET while central parts of Sudan, South Sudan, southwestern Ethiopia and northeastern Uganda showed increases. Except for a decline in 4915 ET in central Uganda, the detected changes in ET (both positive and negative) were not associated with corresponding changes in rainfall. Detected declines in ET in the Nile delta and Ethiopian highlands were found to be attributable to anthropogenic land degradation, while the ET decline in central Uganda is likely caused by rainfall reduction.


Introduction
Terrestrial evapotranspiration (ET)-the water transferred from the soil-plant complex to the atmosphere-plays an important role in water and energy budgets [1].Next to precipitation, ET constitutes the largest component of the terrestrial water budget [2].About three-fifths of the net solar radiation absorbed by earth's surface is used in ET and the accompanying sensible heat flux [3].Characterizing the variation of terrestrial ET is also a key issue for understanding terrestrial components of the global water cycle given that about 40% of terrestrial precipitation originates from terrestrial ET; while about 57% of all terrestrial ET returns as terrestrial precipitation [4].Consequently, accurate estimates of ET are key environmental constraints for climate and hydrological simulations [5][6][7], and for ensuring sustainable use of water resources [5,8] and hence the ability of biological systems to support human needs.This is especially true in the Nile Basin (Figure 1), a highly populated region in northeastern Africa, where surface water primarily originates from rainfall over an area encompassing about 25% of the Basin, with substantial parts of the Basin (about 40%) remaining water-limited and dependent on the Nile River for irrigation.The vulnerability of the Nile Basin region to climate [9] and land use [10,11] change and the responses of the hydrological cycle to these changes remain poorly understood.
Several studies have investigated the links between vegetation, ET and climate in Africa [10][11][12][13][14][15][16][17][18][19].For instance, Jung et al. [12] analyzed long-term in situ, remote sensing and meteorological measurements, and attributed decreasing trends in ET in East Africa from 1998-2008 to decreasing trends in soil moisture.Marshall et al. [19] analyzed 31 years of trends in ET with trends in rainfall and surface temperature using remote sensing and surface climate reanalysis data; and suggested observed declines in ET in the Ethiopian highlands could be attributable to over-cultivation leading to land degradation.Pricope et al. [10] studied the spatial relationships between climate, vegetation change, degradation and population density in the East African Horn.Using the parametric statistical approach of linear regression trend tests with annual data, Alemu et al. [20] also detected declining ET trends in central portions of Uganda, eastern Sudan and Ethiopian highlands, but did not attempt to relate these trends directly to drivers.Other researchers found that both climatic and anthropogenic factors influence vegetation dynamics and thereby ET in east Africa [10,11].Nonetheless, the interactions between ET and climatic drivers in Africa remain relatively unexplored [19].
Access to sufficiently reliable and consistent in situ hydrological data in the Nile Basin remains a major challenge [21], despite the monumental need to accurately assess and monitor the limited water and land resources in the Basin, and the growing demand for water from an increasing population.However, the availability of long-term, basin-wide temporally and spatially consistent remote sensing data provides an opportunity to consistently and reliably provide geophysical and land resource data in a geographically comprehensive manner.Remote sensing provides the only feasible means to estimate ET over larger regions [2].In addition, the reliability of remote sensing data for accurately estimating spatially explicit hydrological and biophysical data has improved significantly [21] to the extent that they can be used with confidence for water accounting and hydrological modeling, with overall accuracies of 95% for ET, 82% for rainfall, and 85% for land use [22].Here we use publicly-available satellite-derived ET, rainfall, vegetation index, and land cover data to analyze ecohydrological dynamics and interrelationships in the Nile Basin during 2002 to 2011.
We analyzed the monthly relationship between ET and rainfall and examined how this relationship was influenced by land cover/land use in the Basin during the period from 2002 through 2011.For this analysis, we estimated per-pixel correlations at 1-km spatial resolutions between monthly time series of ET and rainfall and compared the results by land cover and subbasins.In order to evaluate the influence of rainfall and anthropogenic factors on changes in ET, we restricted our analysis to spatial associations of trends in rainfall and NDVI with corresponding trends in ET in the vegetated parts of the basin.For the spatiotemporal analysis we assessed pixel-scale temporal trends in ET, rainfall and vegetation greenness, spatially synthesized individual trends, and analyzed the relationships.Our objective was to address the following questions: (1) How are monthly ET dynamics related to rainfall in the Nile Basin?(2) How does this relationship vary with land cover?(3) What are the geographic extents and patterns of trends the Basin exhibited by ET, rainfall and vegetation dynamics?(4) Which regions exhibited change in ET as a result of rainfall-driven or anthropogenic factors?

Study Area
The Nile Basin (Figure 1) is a major river basin in northeastern Africa and is drained by the longest river in the world.Home to 238 million people and shared by eleven countries, the basin traverses an extremely wide band of latitude from 4° S to 32° N. Its geographic location and latitudinal extent lead to a pronounced north-south rainfall gradient that divides the basin into three main rainfall regimes: arid, tropical and equatorial [23] with a corresponding gradient of vegetation productivity (Figure 1b).
Total basin area is ~3.7 million km 2 .A water mask was used to exclude major water bodies of the basin (including Lakes Victoria, Albert, Tana, Kyoga, Edward) that constitute ~3% of the entire basin area.As this study is primarily a comparison of an operational ET product (SSEBop) with other datasets, a mask was generated for ET pixels with 0% frequency 8-day time series during the study period.Additionally, a mask was also generated for non-vegetated portions of the basin which corresponds to pixels with median of positive NDVI < 0.12 during 2002-2011 (cf. Figure 1b).The non-vegetated area mask covered nearly one-third (29%) of the entire basin area.The remaining unmasked 2.5 million km 2 of vegetated land surface comprised the study area.

Land Cover
The Moderate Resolution Imaging Spectroradiometer (MODIS) land cover type product (MCD12Q1) is a global map of land cover types derived using observations collected by the MODIS sensors onboard the Aqua and Terra satellites [24].It is produced at annual time steps with a spatial resolution of 500 m and comes in five different classification schemes.We used the 2007 land cover product with the International Geosphere Biosphere Programme (IGBP) classification scheme, which contains 17 land cover classes, as representative of the study period (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011).The 16 land cover units found in the Nile Basin (Figure 1a) were thematically aggregated into 9 major land cover classes, viz., (1) Croplands; (2) Croplands and Natural Vegetation mosaic (CNV); (3) Forest; (4) Savannas; (5) Woody Savannas; (6) Shrublands; (7) Grasslands; (8) Permanent wetlands; and (9) Bare Lands.We synthesized two new land cover classes (i) Forest, by aggregating 5 IGBP classes (Evergreen Needleleaf Forest, Evergreen Broadleaf forest, Deciduous Needleleaf Forest, Deciduous Broadleaf Forest, Mixed Forest) and (ii) Shrublands (Closed Shrublands and Open Shrublands).The resulting land cover map (Figure 1a) was resampled to 1-km to match the spatial resolution of ET and vegetation data using the majority class filter as described in Kaptué et al. [26].Our data encompassed 10 MODIS tiles (h20v05-09 and h21v05-09), which were acquired from Reverb, a metadata and service discovery tool [27].

Evapotranspiration Data
The monthly Operational Simplified Surface Energy Balance (SSEBop) ET product [30] at 1-km spatial resolution (Figure 2b shows mean annual ET average for 2002-2011) was acquired from the USGS Earth Resources Observation and Science (EROS) Center.The SSEBop is an operational parameterization of the Simplified Surface Energy Balance (SSEB) model [31].It uses MODIS-derived land surface temperature and model-assimilated meteorological data to generate a gridded ET product [30].

Precipitation Data
Daily rainfall data were obtained from the Tropical Rainfall Measuring Mission (TRMM) at 0.25° spatial resolution between 2002 and 2011 (product 3B42, version 7) [32].The daily data (in mm/day) were temporally aggregated to monthly time series and resampled to 1-km spatial resolution using bilinear interpolation technique to match the spatial resolution of ET and vegetation data.The data were acquired from NASA's TRMM webpage [33].Figure 2a presents rainfall distribution across the basin using mean annual rainfall during 2002-2011.Moreover, the difference between rainfall and ET showing key water source and sink areas is shown in Figure 2c.

Vegetation Data
As a proxy for vegetation productivity, the Normalized Difference Vegetation Index (NDVI) data were generated from the MODIS nadir bidirectional reflectance distribution function (BRDF) adjusted reflectance (NBAR) product (MCD43B4) [25].The NDVI was calculated for each 16-day NBAR composite (updated every 8 days) as the ratio of the difference between the near-infrared (0.841-0.876 μm) and red (0.620-0.670 μm) bands normalized by the sum of the two bands.A monthly time series of NDVI data was then compiled using the maximum composited NDVI for each month during 2002-2011.The 10 MODIS NBAR tiles (h20v05-09 and h21v05-09) were acquired from Reverb, a metadata and service discovery tool [27].

Correlation Analysis
To investigate the interactions between rainfall and ET from 2002 to 2011, nonparametric Kendall rank correlations [34,35] between monthly time series data of ET and rainfall were calculated and tested against the null hypothesis that the correlation values were not significantly different from zero.To analyze the relationship between land cover and ET-rainfall interactions, the median of ET-rainfall correlations were calculated for individual land cover classes in all subbasins.

Trend Analysis
The nonparametric Seasonal Mann-Kendall trend test (SMK) [34][35][36] was applied to each of the three time series (ET, rainfall, and NDVI).For each time series, the resulting SMK statistic and associated p-value maps were combined to produce a three-class map (Positive (p), Negative (n), Not Significant (o)) at three significance levels: p < 0.05, p < 0.01 and p < 0.001.
To attenuate noise and remove individual outlier pixels, a 3 × 3 median filter was applied to both correlation (Figure 3) and trend maps (Figure 5).The Kendall rank correlations and SMK tests were performed using Kendall library [34] in the R software environment [37].[30] and TRMM rainfall [32] at spatial resolution of 1-km.Numbered polygons are major subbasins of the Nile.Major water bodies are shown in grey.

Synthesis of Trends
Spatial association between the trend maps were used to evaluate how trends of ET were associated with trends of rainfall and NDVI.To analyze how changes in ET associate with significant changes and/or no-changes in rainfall and NDVI, we combined the two trend categories of ET (positive and negative, p < 0.05) with three trend categories of rainfall and NDVI (Positive, Negative, Not Significant classes, p < 0.05) and generated a geodatabase with 18 classes (Table 1).

Temporal Correlations across Space
Figure 3 illustrates the spatial distribution of the Kendal tau (τ) rank correlation coefficients between monthly time series of ET and rainfall from 2002 through 2011.The dominance of statistically significant (p < 0.05) positive τ-values, covering about two-thirds of the study area, strengthens the conclusions from earlier reports that showed that rainfall is the key driver of seasonal and interannual variability in ET in the region [12,20,38].About 28% of the study area (2.5 M•km 2 ) showed no statistically significant relationship.
Areas with relatively strong positive correlations (light blue to dark blue) were observed stretching across the Sahelian belt (10° N to 15° N) from western Sudan to Ethiopian highlands (subbasins 3-6).This pattern agrees with the consensus that rainfall is the dominant limiting factor on semi-arid vegetation growth and ET [13,[38][39][40].Farther south, a diagonal band of positive correlation pixels was identified stretching from southwestern South Sudan (subbasin 7) to northeastern Uganda (subbasin 10).The southern portions of the basin surrounding Lake Victoria (subbasins 10 and 11), a region characterized by wetter climate, also showed positive correlation pixels.
Areas with negative correlations (yellow to dark-red, Figure 3) constituted ~7% of the study area.In subbasin 1, the Nile delta and valley in Egypt exhibited relatively strong negative correlations (−0.5 ≤ τ < −0.2).This pattern occurs because the Nile delta and valley are intensively irrigated agricultural regions [41] in an arid climate with very limited rainfall [20,23].Crop productions in those regions are almost wholly dependent on irrigation from the Nile River [41,42].

Land Cover-ET-Rainfall Relationships
The color-coded matrix in Figure 4 shows ET-rainfall correlation coefficients with respect to land cover types for each of the eleven subbasins.Negative correlation coefficients were the results of the land cover decoupling the local effects of precipitation on ET, through the provision of an alternate source of water to supply ET demands.On the other hand, the interaction remained positive when there is little or no influence of land cover.Land cover types that exhibited positive ET-rainfall interactions included natural vegetation types (shrublands, savannas, grasslands and woody savannas), CNV, croplands, and bare lands.In subbasins that have wetter hydro-climatic characteristics such as subbasin 8, where moisture availability is supplied by sources other than rainfall (e.g., wetlands, soil moisture or base flow), woody savanna and grasslands exhibited negative ET-rainfall relationships.Negative correlations observed over CNV and croplands in subbasin 1 resulted from the subbasin being almost entirely irrigated as rainfall is negligible.

Temporal Dynamics
Seasonal Mann-Kendall (SMK) trend results for monthly time series of ET, rainfall and monthly-maximum NDVI data during 2002-2011 are presented in Figure 5 at three different significance levels (0.05, 0.01 and 0.001).

ET Patterns
The extent of statistically significant trends in ET detected in the basin was limited, with just 12% and 6% of the unmasked area showing positive and negative ET changes, respectively, at a p-value of 0.05 or less (Figure 5a).The central regions of the basin in parts of Sudan and South Sudan (subbasins 5,8,9) showed an overall positive change in ET.In addition, small pockets of positive ET changes were identified in eastern Ethiopian highlands (subbasins 3,4), northeastern Uganda (subbasin 10) and Rwanda/Burundi regions (subbasin 11).The ET trends in the inland portions of the Nile delta and Nile valley in subbasin 1 were negative (Figure 5a), in contrast with the positive ET trends detected on the fringes of the delta.In addition, smaller spots of negative ET trends were detected in eastern Sudan (subbasins 3,4), the western Ethiopian highlands (subbasins 3,4), in regions surrounding the Sudd Wetlands in South Sudan (subbasin 8), central Uganda (subbasin 10) and Tanzania (subbasin 11) regions.The declines in ET in the Ethiopian highlands were in agreement with the findings of Marshall et al. [19].As the significance levels are tightened, the corresponding areas of positive and negative trends contract.Positive trends occurred at 5% and 2% for p < 0.01 and p < 0.001, respectively; whereas, negative trends occurred at 3% and 1% for p < 0.01 and p < 0.001, respectively.
Significant declines in ET in the agricultural lands of the Nile delta and valley in Egypt have been observed in subbasin 1.This was likely caused by a combination of man-made land degradation (chemical and physical) that occurred in the delta.Several researchers have reported soil degradation through increased soil salinization as a result of saline drainage water reuse for irrigation [43][44][45].Moreover, some portions of the delta also exhibited physical land degradation from urban land encroachment on agricultural areas [46,47].However, the fringes of the Nile delta showed increases in ET, which were likely a consequence of expanding irrigation in newly reclaimed agricultural lands [48][49][50].

Rainfall Patterns
The results of the monthly rainfall trend analysis over 10 years in Figure 5b showed that substantial parts (88%) of the unmasked area experienced no statistically significant (p < 0.05) change.However, 11% and 3% of the area were characterized by declines in rainfall at the 0.05 and 0.01 significance levels respectively.The most conspicuous regions where rainfall declines were detected include along the Uganda-South Sudan border (subbasins 8 and 10), and in central Uganda, north of Lake Victoria (subbasin 10).In general, the rainfall dynamics in the basin during 2002-2011 showed no trends, but given high interannual variation in timing of rainfall [51], the length of the period examined was likely too short to reveal significant trends.
These observations in rainfall trends in the basin were also reported by Alemu et al. [20] who used linear regression to detect trends in annual rainfall using TRMM data (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011).They found large-scale reductions in rainfall in South Sudan-Uganda border (subbasins 8 and 10), central Uganda (subbasin 10) and eastern Sudan (subbasin 3 and 4).A recent study by Hoscilo et al. [11] who also used linear regression to detect trends in annual rainfall (FEWSNET RFE 2.0 data for 2001-2010) also confirmed declines in rainfall across large swaths of central Uganda (subbasin 10).However, our results largely diverge from the Hoscilo et al. [11] findings.Our results indicate no part of the unmasked area exhibited increase in rainfall with the exception of few pixels bordering Lake Tana and Lake Victoria.Other researchers using gauge data also found no signs of recent rainfall trends in eastern parts of subbasin 4 [52,53].

Vegetation Patterns
The NDVI trend map generated using the monthly maximum NDVI time series (Figure 5c) showed that a sizeable 36% of the study area exhibited statistically significant (p < 0.05) changes-12% positive and 24% negative in vegetative greenness.The affected area was 7% (4%) positive and 16% (10%) negative at p < 0.01 (p < 0.001).
The Ethiopian highlands covering subbasins 3 and 4 exhibited the most geographically contiguous decline in vegetative greenness.This pattern agrees with Pricope et al. [10] who found declining trends in NDVI in Ethiopia during 1999-2011.Considerable vegetative decline was detected in the Nile delta/valley regions in subbasin 1 while both sides of the fringes of the Nile delta exhibited increasing greenness and much of the inland delta exhibited declines.The changes in vegetative greenness in the Nile delta and valley followed similar patterns as the ET changes.Increases in greenness in the fringes of the delta are indicators of the expansion of irrigation activities in areas that were initially bare lands.The decline in greenness in the inland parts of the delta may result from a combination of anthropogenically-induced chemical (soil salinity) and physical (loss of agricultural fields) land degradation caused by irrigation by reuse of saline drainage water, population pressure leading to loss of agricultural fields to urban sprawl [43][44][45][46][47].
Regions of central Uganda (subbasin 10), south of Lake Victoria in Tanzania (subbasin 11) and western (subbasin 6) and eastern (subbasins 3 and 4) Sudan also exhibited vegetative declines.The regions with a significantly positive NDVI trend were mainly located in the central parts of the basin (subbasins 5,8,9), and also in smaller pockets in western Sudan (subbasin 6) and northeastern Uganda (subbasin 10).

Synthesis of Individual Trends
Using individual trends in ET, rainfall and NDVI, we synthesized 18 combinations (Table 1) of trends for regions where significant changes in ET were detected.The geographic extent of synthesized combination of trends spanned over an area of 0.5 M km 2 or nearly one-fifth of the unmasked area.Positive ET changes represented about two-thirds (~33.5 thousand km 2 ) of this area, while the remaining one-third (~17.7 thousand km 2 ) were negative.Only five of the 18 classes constituted at least 1% of the total study area (~20 thousand km 2 ) each.With respect to the trend direction in ET, these five classes can be placed into two groups: (1) where ET increases (poo & pop); and (2) where ET decreases (non, noo and nnn).Altogether, these two groups constituted 18% of the entire study area, and 88% of the area that exhibited changes in ET.The synthesis map of these 5 classes as Groups 1 and 2 revealed the distinct geographic distribution between the two groups across the basin (Figure 6).Group 1 includes areas in which ET exhibited increasing trends with no trend in rainfall, and either no significant change in NDVI (poo) or a significant positive trend in NDVI (pop).Group 1 covers 14% of the study area and 57% of the areas that exhibited change in ET.Occurrences of Group 1 tended to cluster together and are largely located in the center of the Basin (subbasins 5,8,9) and, to a lesser extent, in northeastern part of subbasin 10 and western part of subbasin 11.ET increases in those regions are not associated with corresponding increases in rainfall.
Group 2 includes areas that exhibited significant negative trends in ET and constitutes about 7% of the study area or 31% of the area that exhibited ET changes.Both classes non and noo represent no trend in rainfall, while non represents downward trending NDVI as well.The third class in Group 2 covers areas in which significant decreasing trends were observed in ET, rainfall, and the NDVI.Group 2 conditions occurred in the Nile delta in subbasin 1, in eastern part of the study area in subbasins 3 and 4, and in regions north and south of Lake Victoria in subbasins 10 and 11 (Figure 6).
Almost the entire regions of the inland Nile delta (subbasin 1) were classified as non.As the delta is highly populated and intensively irrigated [41], the decline in ET and NDVI is the likely consequence of soil salinization (as saline drainage water is used for irrigation) and/or agricultural land degradation (as urban encroachment on agricultural fields) [46,47].This is in contrast to the fringe regions of the delta that were classified as pop, indicating an increase in both ET and NDVI with no change in rainfall.This finding in the peripheries of the delta is in agreement with other researchers who identified expansion of irrigation activity in newly reclaimed land in the fringes of the delta [48,49].
The decline in NDVI in association with the decline in ET (class non) was mainly observed in subbasins 3 and 4 and is likely caused by non-rainfall-driven land degradation.This non-rainfall-driven degradation in the Ethiopian highlands is likely from increasing population pressure and has been reported by other researchers [10,11,54].On the other hand, class noo was mainly observed around the Sudd Wetlands in subbasin 8, and can be explained by the dynamics of wetland ecosystems where ET decline is independent of rainfall.Class nnn was observed in subbasin 10 north of Lake Victoria.A previous study also found linear trends of declining rainfall and ET in central parts of Uganda (subbasin 10) during 2002-2011 [20].

Conclusions
We have explored seasonal ET dynamics as well as multiyear trends, attempting to understand the relationships between ET, rainfall and land cover both in terms of monthly variability as well as the changes that occurred over the period of 2002-2011.Results from monthly ET-rainfall correlations indicated that 65% of the basin showed positive relationships while only 7% showed negative.The semi-arid regions of the basin (Sahelian belt), large parts of the Ethiopian highlands, the regions surrounding Lake Victoria exhibited positive ET-rainfall relationships, as expected.On the other hand, the Nile valley and delta, the Sudd Wetlands and eastern South Sudan and southwestern Ethiopia exhibited negative ET-rainfall relationships.Analysis of the monthly ET-rainfall relationships with land cover showed that the relationship is generally negative over forests and wetlands; while it is positive in most subbasins over woody savannas, savannas, shrublands, and grasslands as well as anthropogenic land covers such as croplands and CNVs.
Analysis and synthesis of the trends in ET, rainfall and NDVI revealed that the inland parts of the Nile delta exhibited declines in ET and NDVI, possibly indicating soil salinization as a result of use of saline drainage water for irrigation and/or agricultural land degradation as a result of urban sprawl in the delta.In contrast, we found increasing trends in both ET and NDVI indicating signs of agricultural expansion and land reclamation in both fringes of the Nile delta.We have also revealed downward trends in ET and NDVI in the Ethiopian highlands where there has been no significant trend in rainfall.Since this is highly populated land with mainly rainfed subsistence agriculture, the decline in both ET and greenness is likely an indication of land degradation.Central parts of Uganda exhibited declining trends in ET and NDVI as well as in rainfall.These declines in both ET and NDVI are likely driven by associated decreases in rainfall.Large areas in the central part of the basin where increasing trends in ET and NDVI were observed were not associated with corresponding increasing trends in rainfall.
It should be noted that the relatively short-time period of the study (10 years) most likely captured the influences of short-term climatic variations rather than slower evolving dynamics resulting from global climate change.As a result, the trends documented in this study should not be extrapolated in to the future, but instead should be viewed as significant changes of historical fluctuations over the decade in climatic conditions and human land use.The various model-assimilated and satellite-derived data products used in the study, including ET, rainfall estimates, vegetation indices, land cover classification, as well as methods employed such as downscaling of rainfall data, all introduce levels of uncertainty, sometimes not well-characterized, that likely increased the noise in the data and limited our ability to detect relatively-short-term trends.
Despite these limitations, the study provides evidence that land use/land cover influence the association between local rainfall and ET and is a likely driver of short-term trends in ET throughout many parts of the basin.Continued basin-wide monitoring of hydrological processes and land use and land cover characteristics using satellite remote sensing will provide a basis for characterizing future changes and continuing to improve our understanding of the drivers of ET in the Nile Basin.

Figure 1 .
Figure 1.Nile Basin.(a) MODIS land cover map (MCD12Q1, IGBP scheme) for the year 2007 [24]; (b) Median of positive NDVI values for the period of 2002-2011 derived from time series of MODIS NBAR (MCD43B4) [25] composite data.Black lines show subbasin boundaries; grey lines show country boundaries; numbers show subbasin ID.

Figure 3 .
Figure 3. Spatial distribution of Kendall rank correlation coefficients ( ) between monthly time series of ET and rainfall.Displayed correlation coefficients are statistically significant (p < 0.05).Data used are positive values of SSEBop ET [30] and TRMM rainfall [32] at spatial resolution of 1-km.Numbered polygons are major subbasins of the Nile.Major water bodies are shown in grey.

Figure 4 .
Figure 4. Color-coded matrix displays statistically significant (p < 0.05) Kendal tau (τ) rank correlation coefficients between ET and rainfall averaged by land cover type for the 11 major subbasins of the Nile basin.Colors range from dark-red (for strongly negative) to dark-blue (for strongly positive) correlations.A 3 × 3 majority filter was applied on the land cover map.Empty boxes indicate an absence of the land cover type in the subbasin.

Figure 5 .
Figure 5. Maps of statistically significant Seasonal Mann-Kendall trends of (a) ET; (b) Rainfall; and (c) NDVI for the period between 2002 and 2011 using monthly time series of monthly ET, monthly rainfall, and monthly maximum NDVI at 1-km spatial resolution.

Figure 6 .
Figure 6.Per-pixel synthesis map of significant trends (p < 0.05) in ET, rainfall and NDVI with 5 major combination classes of trends.Gray areas represent regions that exhibit no significant trend in ET.The colored 5 major classes encoding the combined trends in ET-rainfall-NDVI order are shown with their proportional area with respect to study area.

Table 1 .
Table showing synthesis of SMK trends (p < 0.05) of monthly ET, rainfall and monthly maximum NDVI.The 18-classes encode the combined trends in ET-rainfall-NDVI order.