Trend and Sensitivity Analysis of Reference Evapotranspiration in the Senegal River Basin Using NASA Meteorological Data

Understanding evapotranspiration and its long-term trends is essential for water cycle studies, modeling and for water uses. Spatial and temporal analysis of evapotranspiration is therefore important for the management of water resources, particularly in the context of climate change. The objective of this study is to analyze the trend of reference evapotranspiration (ET0) as well as its sensitivity to climatic variables in the Senegal River basin. Mann-Kendall’s test and Sen’s slope were used to detect trends and amplitude changes in ET0 and climatic variables that most influence ET0. Results show a significant increase in annual ET0 for 32% of the watershed area over the 1984–2017 period. A significant decrease in annual ET0 is observed for less than 1% of the basin area, mainly in the Sahelian zone. On a seasonal scale, ET0 increases significantly for 32% of the basin area during the dry season and decreases significantly for 4% of the basin during the rainy season. Annual maximum, minimum temperatures and relative humidity increase significantly for 68%, 81% and 37% of the basin, respectively. However, a significant decrease in wind speed is noted in the Sahelian part of the basin. The wind speed decrease and relative humidity increase lead to the decrease in ET0 and highlight a “paradox of evaporation” in the Sahelian part of the Senegal River basin. Sensitivity analysis reveals that, in the Senegal River basin, ET0 is more sensitive to relative humidity, maximum temperature and solar radiation.


Introduction
One of the current challenges in hydrological studies is the assessment of the impact of climate change on water resources [1]. Indeed, streamflow projections allow water resource managers as well as decision-makers to compare estimated future water availability with water demands and thus develop strategies to better manage water, particularity in the context of climate change. Forecasting future streamflow is highly dependent on evapotranspiration, which is influenced by several climatic variables [2]. In addition, evapotranspiration plays a crucial role in determining crop water requirements [3], essential when planning and programming irrigation. Evapotranspiration is a key variable that establishes a relationship between the energy and the water balances on the earth's surface [4]. Thus, it is one of the parameters which control the transfer of water and energy between the earth and the atmosphere and is considered as an indicator of climate change [5].
Evapotranspiration is influenced by several factors: temperature, wind speed, relative humidity, solar radiation, plant cover, land use, etc. [6][7][8]. Variation in any of these factors can affect the spatiotemporal distribution of evapotranspiration [9][10][11][12]. In recent years, several authors [4,5,[13][14][15][16] have shown that global warming leads to an increase in evapotranspiration. Indeed, temperature increases due to climate change and the associated higher amount of available energy leads to an increase in evapotranspiration [13]. On the other hand, a decrease in evapotranspiration has been noted in certain climatic regions of the world: in dry and humid climates of the United States [17,18], in tropical and subtropical climates in China [13,14,19], in semi-arid climate of Turkey [20,21], and in arid climate of Iran [22]. The contrast between the increase in temperature and the decrease in evapotranspiration is known as the "evaporation paradox" [13,23]. This "paradox of evaporation" could be explained by the decrease in solar radiation, following the increase in cloud cover and the concentration of aerosols and pollutants in the atmosphere [17,24,25], the increase in relative humidity [26] and the reduction in wind speed [25].
In this regard, several studies have been interested in characterizing the trend of evapotranspiration and the climatic variables which condition it [4,12,16,27,28]. The method generally used by these authors are: (i) the use of the Mann-Kendall test for detecting the trend of ET 0 and climatic variables [7,13,27,29] (ii) and an analysis of the sensitivity of evapotranspiration to climatic variables [1,[30][31][32][33][34][35][36][37]. Sensitivity analysis helps to understand the link between evapotranspiration and climatic variables [31], to choose the methods for evapotranspiration estimation [30], and to detect the impact of climate change on evapotranspiration [16]. For example, Darshana et al. [34] used Mann-Kendall's test and sensitivity analysis to detect the trend of evapotranspiration and climatic variables in the Tons basin (India). Their results showed a significant decrease in annual and seasonal evapotranspiration; sensitivity analysis reveals that evapotranspiration is more sensitive to maximum temperature and solar radiation. Gong et al. [31] found that evapotranspiration is more sensitive to relative humidity, solar radiation, average air temperature, and wind speed, respectively, in the Changjiang Basin in China. Tabari et al. [36] noted that mean temperature, solar radiation, relative humidity and wind speed, respectively, are the main parameters that influence evapotranspiration in western Iran. These different studies show that the trend and sensitivity of evapotranspiration to climatic variables depend on the study regions. In Africa, Obada et al. [8] analyzed recent and future variations in evapotranspiration in Benin over the period 1981-2010 using the Lamb index [38] and the Mann-Kendall test. They noted a significant downward trend (95% of studied stations) in annual evapotranspiration during this period with high inter-annual variability. However, they did not analyze the trend of climatic variables (temperature, wind speed, relative humidity, and solar radiation) and their influence on evapotranspiration. Abiye et al. [7] used monthly data from the Climate Research Unit (CRU) (version TS 4.00) to analyze the trend of evapotranspiration over the period 1906-2015 on the scale of West Africa. Their results showed a significant upward evapotranspiration trend in the Sahelian zone and a downward trend in the Sahelian and Guinean areas. Their study showed that temperature is the key climatic parameter that explains evapotranspiration variation in West Africa. Although the study by Abiye et al. [7] is relevant in many ways, they did not include wind speed, humidity and solar radiation in the sensitivity analysis. The works of Gong et al. [31] in China, Irmak [30] in the United States, Bakhtiari and Liaghat [39] in Iran, Ambas and Baltas [40] in Greece, have shown that these climatic variables have an influence on evapotranspiration.
Thus, this work has a double objective: (i) detect the trends of annual and seasonal evapotranspiration at the scale of the Senegal River watershed and (ii) analyze the sensitivity of ET 0 to climatic variables (maximum temperature, minimum temperature, wind speed, relative humidity and solar radiation). The Senegal River basin is the second largest transboundary basin in West Africa after the Niger basin. It is at the heart of the development issues and challenges of its riparian countries: Guinea, Mali, Senegal and Mauritania and is managed by the Senegal River basin Organization (in French, Organisation pour la Mise en Valeur du fleuve Sénégal, OMVS). Its surface area is estimated at more than 300,000 km 2 [41] and its population was 6.5 million in 2015 and will reach 9.8 million in 2025 [42]. Senegal River basin water resources are used for irrigation (2 billion m 3 year −1 ), water supply for the capital cities of Senegal and Mauritania, hydroelectricity, water supply for livestock, industries and mining, navigation and fishery. The incomes of almost 6 million people depend on the Senegal River basin water resources [42]. Therefore, the study of the reference evapotranspiration which is one the main elements of the hydrological cycle is crucial, particularly the investigation of its long-term trends.

Study Area
The Senegal River basin is divided into two large geographic units: the upper basin and the lower basin. The upper basin is almost essentially centered on the mountain ranges of Fouta-Djalon in Guinea, source of the main flows in West Africa, and on the hills and plateau as the river flows north into Mali. This area is the rainiest part of the Senegal river basin and the source of the main tributaries (Bafing, Bakoye, Baoulé). In addition, the main hydraulic infrastructures (Manantali, Felou dams) were built in this area by OMVS. These infrastructures play the main role in river flow regulation, the production of hydroelectricity and the development of irrigation. In order to regulate the uncontrolled tributaries and increase the energy productivity of the basin, several other dams are planned in this area by OMVS. The lower basin, which corresponds to the Valley and the Delta, concentrates the main agricultural activities in the basin. From a climatic point of view, Dione [43] identified four main climatic zones on the basis of the latitudinal distribution of precipitation: Guinean (average annual rainfall p > 1500 mm); South Sudanian (1000 < p < 1500 mm); North Sudanian (500 < p < 1000 mm) and Sahelian (p < 500 mm). In this study, the Sudanian zone (South Sudanian and North Sudanian) has been considered as one climate zone. Thus, three climatic zones are considered: Guinean (1% of the river basin area), Sudanian (62%), and Sahelian (37%).

Data
In West Africa, climate data managed by national meteorological services are difficult for researchers to access because of their high cost of acquisition [44]. In addition, the low density of the observation network poses a problem of representativeness of these data at the scale of the watershed. However, a large-scale study of evapotranspiration requires several measurement points due to the heterogeneity of the landscapes and the variation in the energy transfer processes [45]. Therefore, in this study, data reanalysis and observation of the POWER project of NASA Langley Research Center (LaRC) funded by the NASA Earth Science/Applied Science program (https://power.larc.nasa.gov/data-accessviewer, accessed on 20 December 2018) were used as an alternative to observed data that are inaccessible and generally scattered [44,46]. These data, already used and validated in previous studies [47][48][49], have the advantage of having spatial and temporal coverage on a global scale [46,[50][51][52] and provide the climatic variables necessary for the estimation of evapotranspiration [46,53]. These datasets which result from satellite and model-derived weather data contain uncertainties however, and where possible these should be evaluated and validated locally against available in situ measurements [54]. When extracting NASA data, the coordinates of 31 weather stations from the OMVS database [44] and 33 virtual stations were used ( Figure 1). The climatic variables extracted on a daily scale over the period 1984-2017 are: temperature (maximum and minimum), relative humidity, wind speed and solar radiation. The average daily values of the main climatic variables are given in Figure 2. The temperature varies from 19 • C to more than 38 • C depending on the climatic zones. The maximum values are observed in the Sahelian zone and those minimum in the Guinean zone. The latter area has the highest relative humidity rates above 60%. zones. The maximum values are observed in the Sahelian zone and those minimum in the Guinean zone. The latter area has the highest relative humidity rates above 60%.   zones. The maximum values are observed in the Sahelian zone and those minimum in the Guinean zone. The latter area has the highest relative humidity rates above 60%.

Reference Evapotranspiration Estimation
Reference evapotranspiration (ET 0 ) is the evapotranspiration of a reference surface (grass or alfalfa), of uniform height (0.08 to 0.15 m) and suitably supplied with water [55]. It is considered as a climatic parameter and therefore can be estimated with only climatic variables [55]. Among the various existing methods for estimating ET 0 , Penman-Monteith is recognized as a standard method in all climatic conditions of the globe [55]. Thus, this method (Equation 1) is used in this study for the estimation of ET 0 of grass on a daily scale. The annual and seasonal ET 0 are obtained by summing daily ET 0 values. Two seasons are considered: a dry season (November to May) and a rainy season (June to October) even though, the length of the seasons varies according to the climatic zones of the basin. The method of interpolation of the Inverse Distance Weighting (IDW) [56,57] is used to spatialize the ET 0 on an annual and seasonal scale. The IDW interpolation technique is widely used to analyze the spatial distribution of evapotranspiration [10,[56][57][58]. It is a simple method that provides accurate estimates of evapotranspiration values [10] compared to more resource-intensive methods like kriging [59]: where ET 0 is the reference evapotranspiration (mm d −1 ), Rn: net radiation on the crop surface (MJ·m −2 d −1 ), G is the density of heat flux of the soil (MJ·m −2 d −1 ) which is ignored on a daily scale, T is the average daily air temperature at a height of 2 m ( • C), u2 is the wind speed at a height of 2 m (m s −1 ), es is the pressure of saturated vapor (kPa), ea is actual vapor pressure (kPa), (es−ea) is the saturation deficit (kPa) ∆ is the slope curve of the vapor pressure (kPa • C −1 ) and γ is the psychrometric constant (kPa • C). Due to the unavailability of maximum and minimum relative humidity in the reanalysis data used, the actual vapor pressure (ea) is estimated from the average relative humidity according to the formula proposed by Allen et al. [55]: where ea is the actual vapor pressure (kPa), Rh mean relative humidity, e • (Tmax) saturated vapor pressure at a maximum temperature, e • (Tmin) saturated vapor pressure at a minimum temperature.

Reference Evapotranspiration Trend
The Mann-Kendall test [60,61], is often used to detect trends in hydroclimatic data series [44,62,63]. The advantage of this test is that it does not require the data to follow any statistical distribution [56]. It is also not sensitive to extreme values [22]. The test is based on two hypotheses: null hypothesis (noted Ho) which supposes that the test is stationary, no trend and the alternative hypothesis (H1) which rejects Ho and indicates the existence of a trend. Mann-Kendall's statistical S is given by the following formula: where x i is the value of the variable at time i and x j is the value of the variable j, n is the length of the series and sign () is a function which is calculated as follows: An upward or downward trend is given by the sign of the value Z (positive or negative) which is obtained from the variance of S which is obtained by this formula: For n > 10, Z follows approximately a normal distribution and can be calculated as follows: The trend obtained is measured according to its level of significance i.e., the probability associated with the rejection or not of the null hypothesis. In this study, the significance level of 0.05 is used. When |Z|> 1.96 the null hypothesis is rejected and the trend is significant at 5%. If a trend is noted in the data series, its amplitude can be determined by the slope of the trend (noted β). The method of Sen [64] is generally used to estimate the slope of the trend: where x i and x j are the data values at times i and j, respectively. A positive β value indicates an upward trend and a negative value indicates a downward trend. The application of the Mann-Kendall test requires a verification of the autocorrelation of the data. If the data are not autocorrelated, the Mann-Kendall test can be applied directly without modification. However, the existence of an autocorrelation can increase the probability of detecting a significant trend [65]. A positive autocorrelation increases the rejection of the null hypothesis when it is true because the variance of the Mann-Kendall test increases with the magnitude of the autocorrelation of the series [5]. Therefore, the influence of autocorrelation should be removed from the series [66]. The autocorrelation function (ACF) was used to construct the autocorrelation graphs [67]. Figure 3 shows the autocorrelation tests of some stations in the basin. Significant autocorrelation (5%) was not observed for 80% of the stations. Nevertheless, some stations show significant autocorrelation but not very representative. Therefore, given the limited number of stations with significant autocorrelation, the Mann-Kendall trend test is directly applied to the raw data.

Sensitivity of ET 0 to Climatic Variables
Sensitivity analysis is the technique generally used to determine the influence of climatic variables on evapotranspiration [1,30]. The daily average variables used in this study are: maximum temperature (Tmax), minimum temperature (Tmin), wind speed (u2), relative humidity (Rh) and solar radiation (Sr). ET 0 sensitivity to climate variables was estimated by varying each variable individually by ±5, 10, 15, 20 and 25%, leaving the other variables constant [36]. A sensitivity coefficient (8) was calculated to determine the degree of influence of each climate variable on the ET 0 [11]: where S vi the sensitivity coefficient, ∂ET 0 variation of the reference evapotranspiration (ET 0 ) caused by the change of a variable ∂ vi , vi is the considered variable. A positive (negative) sensitivity coefficient indicates that the variable increases (decreases) the reference evapotranspiration. The higher the absolute value of the sensitivity coefficient, the more the variable has influence on ET 0 [11].
modification. However, the existence of an autocorrelation can increase the probability of detecting a significant trend [65]. A positive autocorrelation increases the rejection of the null hypothesis when it is true because the variance of the Mann-Kendall test increases with the magnitude of the autocorrelation of the series [5]. Therefore, the influence of autocorrelation should be removed from the series [66]. The autocorrelation function (ACF) was used to construct the autocorrelation graphs [67].   Figure 4 shows the spatial distribution of seasonal and annual reference evapotranspiration in the Senegal river basin averaged over 1984-2017. Reference evapotranspiration varies spatially from 1130-2170 mm during the dry season and decreases during the rainy season with values between 400 to 1500 mm. Annual ET 0 varies from 1500 to more 3288 mm with an average value of 2524 mm (i.e., between 4.1 and 9.0 mm d −1 , mean 6.9 mm d −1 ). These values, even if they seem high, are in line with those of Aubin [68] who showed that in West Africa, annual evapotranspiration varies from 3500 mm in the north of the Sahel to 1000 mm in the equatorial zone. From a spatial point of view, the distribution of annual ET 0 follows the opposite direction of the latitudinal distribution of precipitation. Indeed, the Sahelian zone which receives less rain records the highest ET 0 values. On the other hand, in the Guinean zone, the high value of relative humidity and the low ones of temperatures favor the drop-in evapotranspiration values.

Trend of Annual and Seasonal Evapotranspiration
Annual and seasonal evapotranspiration trends and their amplitude are given in Figures 5 and 6. Figure 5 shows that there is no spatial organization of the trend of evapotranspiration in the Senegal River basin. At the annual scale ( Figure 5), 32% of the basin area shows a significant increase in ET 0 . On the other hand, less than 1% of watershed area exhibits a significant decrease in ET 0, mainly in the Sahelian zone of basin. The amplitude of the annual ET 0 variation ranges from −15.5 to 13.8 mm year −1 with an average value of 0.72 mm year −1 (Figure 6). From a spatial point of view, there has been a significant increase in evapotranspiration in the Sahelian and Sudanian zones at annual and dry season scales. However, a significant decrease in ET 0 is observed locally in the Sahelian and Guinean zones. The annual ET 0 increases by 1.63 and 1.53 mm year −1 in the Sahelian and Sudanian zones and decreases by 2.59 mm year −1 in the Guinean zone (Table 1). These results are in agreement with those of Abiye et al. [7] who noted an increase in ET 0 of 0.02 and 0.07 mm year −1 in the Sahelian and Sudanian zones and a decrease of 0.02 mm year −1 in the Guinean zone of West Africa over 1906-2015. However, even if the trend is identical the amplitudes of variation differ significantly, revealing the intensification of ET 0 changes in the past 30 years.
Sensitivity analysis is the technique generally used to determine the influence of climatic variables on evapotranspiration [1,30]. The daily average variables used in this study are: maximum temperature (Tmax), minimum temperature (Tmin), wind speed (u2), relative humidity (Rh) and solar radiation (Sr). ET0 sensitivity to climate variables was estimated by varying each variable individually by ± 5, 10, 15, 20 and 25%, leaving the other variables constant [36]. A sensitivity coefficient (8) was calculated to determine the degree of influence of each climate variable on the ET0 [11]: where the sensitivity coefficient, ∂ET0 variation of the reference evapotranspiration (ET0) caused by the change of a variable , vi is the considered variable. A positive (negative) sensitivity coefficient indicates that the variable increases (decreases) the reference evapotranspiration. The higher the absolute value of the sensitivity coefficient, the more the variable has influence on ET0 [11].     zones. The annual ET0 increases by 1.63 and 1.53 mm year −1 in the Sahelian and Sudanian zones and decreases by 2.59 mm year −1 in the Guinean zone (Table 1). These results are in agreement with those of Abiye et al. [7] who noted an increase in ET0 of 0.02 and 0.07 mm year −1 in the Sahelian and Sudanian zones and a decrease of 0.02 mm year −1 in the Guinean zone of West Africa over 1906-2015. However, even if the trend is identical the amplitudes of variation differ significantly, revealing the intensification of ET0 changes in the past 30 years.    On a seasonal basis, 32% of the basin shows a significant positive trend in ET 0 during the dry season ( Figure 6). During the rainy season, a significant decrease trend is observed for 4% of the basin. 97% of this decrease of annual ET 0 is localized in the Sahelian zone of basin. However, the ET 0 increases significantly in the Guinean zone during the rainy season. The magnitude of the trend in seasonal evapotranspiration varies slightly from −8.36 to −9.27 mm year −1 and from −6.56 to 5.32 mm year −1 during the dry and rainy seasons, respectively. Overall, the analysis of annual and seasonal trends in ET 0 shows contrasting results (decline and increase) depending on the climatic zones of the basin. These results are in line with those of Abiye et al. [7] who showed: (i) a significant positive trend in ET 0 in the Sahelian and Sudanian areas and (ii) a negative trend in the Guinean area of West Africa, but this study reveals how these trends have intensified over 1987-2017. The positive trend of ET 0 found in this study is in agreement with the studies of Azzizadeh and Javan [62], which showed a significant increase in annual and seasonal ET 0 in the northwest of Iran in a semi-arid climate. Figure 7 gives the spatial distribution of the Z values of the Mann-Kendall test of the climatic variables on the annual and seasonal scales. On an annual scale, temperature (max and min), solar radiation (Sr) and relative humidity (Rh) generally show an increasing trend (Figure 7). In fact, maximum and minimum temperatures increase significantly for 68% and 81% of the basin area and relative humidity increases significantly for 37% of the area. Solar radiation stays unchanged for 97% of the basin; while, the wind speed (u2) shows a significant negative trend for 8% of the basin. From a spatial point of view, Figure 8 gives the distribution of the Sen slopes of the different climatic variables (wind speed, maximum and minimum temperature, solar radiation and relative humidity) and Table 1  basin show a significant negative trend in wind speed. Minimum temperature and relative humidity continue to increase, but a significant rising trend in maximum temperature and solar radiation is noted for 6% and 5% of the basin, respectively. From a spatial point of view, a decrease in wind speed, maximum temperature and solar radiation is noted in the Sahelian zone during the rainy season. In the Guinean zone, a decrease in radiation of 0.01 MJ·m −2 d −1 is observed during the dry season (Table  1).  Overall, wind speed is the only variable that decreases significantly on an annual scale and during the rainy period. Similar results have been obtained by Zeng et al. [69] who observed a significant decrease in wind speed during the period 1960-2018 in southwest China. In addition, Oguntunde et al. [70] also found a significant decrease in wind speed over the 1973-2008 period in Nigeria. McVicar et al. [71], based on a review of 148 studies carried out around the world, noted an average decrease in wind speed of −0.017 m year −1 . The decrease in solar radiation found in this study has also been observed in Mongolia [72] and in the Yangsu and Lijiang basins in China [13,56]. During the dry season, maximum and minimum temperatures increase respectively for 88% and 52% of the basin. Relative humidity and solar radiation increase respectively for 25% and 43% of the basin. Wind speed decreases significantly for only 2% of the area. During the wet season, 16% of the basin show a significant negative trend in wind speed. Minimum temperature and relative humidity continue to increase, but a significant rising trend in maximum temperature and solar radiation is noted for 6% and 5% of the basin, respectively. From a spatial point of view, a decrease in wind speed, maximum temperature and solar radiation is noted in the Sahelian zone during the rainy season. In the Guinean zone, a decrease in radiation of 0.01 MJ·m −2 d −1 is observed during the dry season (Table 1).

Trend of Annual and Seasonal Climatic Variables
Overall, wind speed is the only variable that decreases significantly on an annual scale and during the rainy period. Similar results have been obtained by Zeng et al. [69] who observed a significant decrease in wind speed during the period 1960-2018 in southwest China. In addition, Oguntunde et al. [70] also found a significant decrease in wind speed over the 1973-2008 period in Nigeria. McVicar et al. [71], based on a review of 148 studies carried out around the world, noted an average decrease in wind speed of −0.017 m year −1 . The decrease in solar radiation found in this study has also been observed in Mongolia [72] and in the Yangsu and Lijiang basins in China [13,56].
Overall, the upward trend in annual evapotranspiration can be explained by the increase in temperatures (maximum and minimum). This confirms at a more detailed scale the results of Abiye et al. [7] who note that temperature is the key factor behind the variation in evapotranspiration in West Africa. Tabari et al. [37] also obtained similar results in Iran. These results confirmed a direct relationship between evapotranspiration and climate change. Indeed, rising temperatures increase evapotranspiration rates and can worsen the arid conditions of an environment and accelerate the process of desertification [73]. The increase in evapotranspiration also increases agricultural water needs, raising the cost of irrigation and potentially affecting agricultural production and yields [74].
Despite the overall upward trend, a localized decline in ET 0 is observed in some zones in the Sahelian and Guinean areas. This negative trend can be explained by the increase in relative humidity and the decrease in solar radiation and wind speed. Lu et al. [29] also obtained similar results in the province of Yangsi in China. For them, the decrease in the duration of sunshine and the wind speed are the main causes of the decrease in evapotranspiration. Bian et al. [72] also noted that wind speed and radiation are the main factors behind the decline in ET 0 in Mongolia. This phenomenon of declining evapotranspiration despite an increase in temperature is called the "evaporation paradox". After the hydrological paradox highlighted in the Sahel by Descroix et al. [75,76], we are also confirming the paradox of evaporation [1,7,29,77] which is not only localized in the Sahelian zone but extends to the Guinean zone of the Senegal River basin. This paradox is more visible in the Sahelian zone where the high temperatures should be accompanied by a generalized increase in evapotranspiration. This phenomenon has recently been observed by Bian et al. [72] in arid and semi-arid areas of Mongolia.
The decrease in ET 0 noted especially in some areas of the Sahelian zone highlights the notion of "paradox of evaporation" which has been debated by certain authors [1,13,29,76]. Two main phenomena explain this decline in ET 0 : (i) the decrease in the duration of sunshine and solar radiation (Global Dimming) which promotes the reduction of the energy available for evapotranspiration [77] and (ii) the reduction in wind speed (Wind Stilling) which causes a decrease in evapotranspiration [24]. Increased cloud cover, concentration of aerosols and pollutants are the primary factors in the global dimming of the atmosphere [77]. Surface roughness and increased plant cover may also be reasonable explanations for the decrease in wind speed [11]. In addition, anthropogenic activities may play a role in ET 0 [78]. Indeed, hydraulic infrastructures, land use by agriculture and irrigation practices can contribute to the variation of evapotranspiration [4]. The increased precipitation may also be a plausible explanation for the decrease in evapotranspiration [25]. Figure 9 shows the relationship between evapotranspiration and climatic variables. The results show that temperatures (maximum and minimum), wind speed and radiation have a positive influence on evapotranspiration. The increase in these variables leads to an increase in evapotranspiration. On the other hand, an increase in relative humidity leads to a decrease in ET 0 . In addition, the influence of climatic variables on ET 0 varies depending on the climatic zones of the basin. In fact, in the Sahelian and Sudanian zones, evapotranspiration is more sensitive to maximum temperature, wind speed and solar radiation, respectively. In contrast, in the Guinean zone, maximum temperature, relative humidity and solar radiation have more influence on ET 0 . For example, on an annual scale, an increase of 5% to 25% of the maximum temperature leads to an increase in ET 0 from 0.48 mm to 2.89 mm in the Sahelian zone, from 0.37 mm to 2.24 mm in the Sudanian zone and by 0.15 mm to 0.89 mm in Guinean climate. Conversely, increasing the relative humidity from 5% to 25% favors a decrease in evapotranspiration from 0.06 mm to 0.38 mm according to climatic zones. These results are in line with those of Li et al. [11] who showed that relative humidity is the only variable that has a negative influence on ET 0 in the Huai basin in China.

Sensitivity of Evapotranspiration to Climatic Variables
show that temperatures (maximum and minimum), wind speed and radiation have a positive influence on evapotranspiration. The increase in these variables leads to an increase in evapotranspiration. On the other hand, an increase in relative humidity leads to a decrease in ET0. In addition, the influence of climatic variables on ET0 varies depending on the climatic zones of the basin. In fact, in the Sahelian and Sudanian zones, evapotranspiration is more sensitive to maximum temperature, wind speed and solar radiation, respectively. In contrast, in the Guinean zone, maximum temperature, relative humidity and solar radiation have more influence on ET0. For example, on an annual scale, an increase of 5% to 25% of the maximum temperature leads to an increase in ET0 from 0.48 mm to 2.89 mm in the Sahelian zone, from 0.37 mm to 2.24 mm in the Sudanian zone and by 0.15 mm to 0.89 mm in Guinean climate. Conversely, increasing the relative humidity from 5% to 25% favors a decrease in evapotranspiration from 0.06 mm to 0.38 mm according to climatic zones. These results are in line with those of Li et al. [11] who showed that relative humidity is the only variable that has a negative influence on ET0 in the Huai basin in China.  Figure 10 shows the spatial distribution of the mean annual and seasonal sensitivity coefficients of evapotranspiration to climatic variables. Overall, relative humidity has the highest sensitivity Figure 9. Sensitivity of evapotranspiration to climatic variables. (u2 wind speed, Tmax maximum temperature, Tmin minimum temperature, Rh relative humidity, Sr Solar radiation). Figure 10 shows the spatial distribution of the mean annual and seasonal sensitivity coefficients of evapotranspiration to climatic variables. Overall, relative humidity has the highest sensitivity coefficients followed by maximum temperature, solar radiation, minimum temperature and wind speed. From a spatial point of view, the influence of maximum temperature is more significant in the Sahelian and Sudanian zone. The impact of wind speed on ET 0 decreases from Sahelian to Guinean zone with the coefficients varying from −0.03 to 0.22 (Figure 10). Solar radiation and minimum temperature have sensitivity coefficients between 0.09 and 1.54 and between 0.04 and 2.12, respectively. The influence of solar radiation on ET 0 is greater in the Guinean area and during the rainy season. These results are in agreement with works in the Senegal River basin [79] which showed that the methods which integrate the wind speed and the temperature are more efficient for the estimation of the ET 0 in the Sahelian and Sudanian zones of the Senegal River basin. And the methods integrating solar radiation and relative humidity are more robust in the wetter Guinean domain. temperature have sensitivity coefficients between 0.09 and 1.54 and between 0.04 and 2.12, respectively. The influence of solar radiation on ET0 is greater in the Guinean area and during the rainy season. These results are in agreement with works in the Senegal River basin [79] which showed that the methods which integrate the wind speed and the temperature are more efficient for the estimation of the ET0 in the Sahelian and Sudanian zones of the Senegal River basin. And the methods integrating solar radiation and relative humidity are more robust in the wetter Guinean domain.  Table 2 gives the average sensitivity coefficients of the different variables depending on the climatic zones. Overall, relative humidity, maximum temperature and solar radiation are the variables that have the greatest influence on evapotranspiration in the Senegal River basin.   Table 2 gives the average sensitivity coefficients of the different variables depending on the climatic zones. Overall, relative humidity, maximum temperature and solar radiation are the variables that have the greatest influence on evapotranspiration in the Senegal River basin. Relative humidity has sensitivity coefficients ranging from −0.40 to −4.93, the maximum temperature from 0.63 to 1.88 and solar radiation 0.35 to 1.28 depending on the climatic zones. Thus, relative humidity is the climate variable that most influences ET 0 in the Senegal river basin. This finding is not in line with those of Abiye et al. [7] which suggested that ET 0 is most influenced by maximum temperature in West Africa. The wind speed has the lowest coefficients, less than 0.16.

Spatial Distribution of Annual and Seasonal Sensitivity Coefficients
Generally, it should be noted that the influence of climate variables on ET 0 varies from one region to another [80]. Indeed, Gong et al. [31] showed that reference evapotranspiration is more sensitive to relative humidity, solar radiation, temperature and wind speed in the Yangzi Basin in China. Li et al. [11] also noted that relative humidity is the variable that has the most influence on ET 0 in the Huai Basin in China. These results are confirmed in this study where ET 0 is more sensitive to relative humidity, maximum temperature and solar radiation. In addition, Rohina et al. [81] found that solar radiation, maximum temperature and wind speed have more influence on ET 0 in Iran. Ndiaye et al. [82] have argued that ET 0 is more sensitive to changes in solar radiation, maximum temperature and wind speed in Burkina Faso. Estévez et al. [32] observed that evapotranspiration is more sensible to air temperature, relative humidity and solar radiation respectively, in the semi-arid climate of Andalusia (Southern Europe). Irmak et al. [30] noted that in arid and semi-arid regions, ET 0 is more sensitive to wind speed and to saturation deficit. This is not fully verified in this study where the wind speed has the lowest coefficients. Furthermore, mutual effects between climate variables can influence ET 0 . An increase in temperature combined with a drop in relative humidity increases the saturation deficit; which results in an increase in evapotranspiration [83]. This finding is confirmed by Espadafor et al. [27] who identified that increases in air temperature and solar radiation and decreases in relative humidity lead to significant increases in ET 0 (up to 3.5 mm year −1 ) in Southern Europe.

Conclusions
The objective of this study was to analyze the trend of seasonal and annual reference evapotranspiration (ET 0 ) over the period 1984-2017 and its sensitivity to climatic variables in the Senegal River basin. Results reveal an increase in the ET 0 at the annual and dry season scales for 32% of the Senegal River basin. This increase in the ET 0 is localized in the Sudanian zone where the main hydraulic infrastructures are located (Manantali and Felou dams). However, during the rainy season, a localized downward trend in the ET 0 is noted in some areas of the Sahelian zone with a rise in temperatures which illustrates the concept of the "evaporation paradox". The sensitivity analysis reveals that evapotranspiration is more sensitive to relative humidity, maximum temperature and solar radiation. Increase in the ET 0 will increase crop water requirements, accentuate the water losses by evaporation in surface reservoirs, reduce the hydroelectricity producible and may lead to increased competition between different water uses. Therefore, the countries sharing the Senegal River basin must look at ways to reduce consumption where possible and optimize the allocation of available resources.
This research is based on 33-year reanalysis dataset. Therefore, the short time series of climate data used may be a limitation of this study. Thus, it would be important to extend the time series when data are available for future researches. It would also be interesting in further investigation, to evaluate the impact of the topography, the dynamics of land use and hydraulic infrastructures, and climate change on the trend of evapotranspiration in the Senegal river basin