Is Climate or Direct Human Inﬂuence Responsible for Discharge Decrease in the Tunisian Merguellil Basin?

: Climate change and direct anthropogenic impact are recognized as two major factors affecting catchment runoff. This study investigated the separate effect of each of these factors for runoff from the important Tunisian Merguellil catchment. For this purpose, more than forty years of hydrological data were used. The methodology was based on hydrological characterization, NDVI index to monitor land use dynamics, and the Budyko approach to specify origin of change. The results show that hydrological change is much more important upstream than downstream. The last three decades display a 40% reduction in runoff. This is associated with the direct inﬂuence of humans, who are responsible for about 78% of the variation in ﬂow. It appears that climate change contributes to less than about 22%. The combination of increased cultivated land and decreased annual rainfall is the main reason for reduced catchment runoff. Consequently, these effects threaten the sustainable runoff, water in reservoirs, and future water supply in general. Ultimately, the available runoff remains an important parameter and a key indicator to guide the choices of decision-makers and practitioners in current and future climatic conditions. This contributes to supporting sustainable management of remaining water resources. of precipitation temperature (T), potential evapotranspiration (E and in the


Introduction
The hydrological cycle is the process upon which all organisms on earth depend for their survival. It is the engine of environmental and ecological conditions on earth [1]. Over the past few decades, several studies have revealed that river flow is linked to many environmental relationships, especially in arid and semiarid areas [2]. Its variations are mainly due to two aspects, namely global climate change and direct human interventions. This contributes to the resulting consequences, a more pronounced scarcity of water resources and more serious environmental deterioration [3][4][5][6].
The accelerated global change leaves little time to adjust. Although the need for adaptation to climate change is recognized, there is also a need for adaptation to change affecting land use [7]. Thus, it is essential to understand the causal chains of the processes of river flow variation. As a result, the scientific community has started to analyze the factors behind changes and the possible consequences of these changes [8]. Several factors influence the catchment runoff, including precipitation, temperature, potential evapotranspiration, topography, soil texture, and land use [9]. Human activities are recognized as a major factor affecting river runoff [10]. This appears in deforestation, land use change, soil and water conservation practices, and exploitation of surface and groundwater reservoirs and water supply [11,12]. Human activities have many negative effects on the

Study Area
The watershed of the Merguellil Wadi is located in central Tunisia. It covers an area of 1183 km 2 . The El Haouareb dam, built in 1989, has a natural outlet ( Figure 1) and retains all water from Merguellil Wadi [30]. The shape of the watershed is almost triangular, the base being downstream (Figure 1). The main stream of the Merguellil Wadi directly drains most of the basin, the tributaries of some importance being located on the right bank. The river bed of the Wadi is very sensitive to erosion [31]. The altitude varies between 200 and 1200 m with an average of 500 m (Figure 1).
The watershed is semiarid, characterized by very high spatial and temporal variability of precipitation. Figure 2 shows the annual precipitation of the nearby Kairouan station during the period 1950-2020, including the study period 1976-2017. The displayed longterm rainfall averages aim to assess the representativity of the study period in relation to long-term rainfall conditions. The average annual precipitation during the 70 years is 307 mm (standard deviation ± 108 mm). The average for the period 1976-2017 is equal to 290 mm (standard deviation ± 89 mm). Thus, it can be said that the study period is representative for longer periods. In Tunisia, the rainfall regime is often subject to a semiperiodicity of about 10 years [32]. In addition, since the 1950s, Tunisia has experienced a decrease of 5% per decade in annual rainfall [33]. The ten-year average precipitation ( Figure 2) shows that there are successive periods of sequential surplus, deficit, and periods close to the long-term average. Figures 3 and 4 show the annual and intra-annual variability of precipitation (P), temperature (T), potential evapotranspiration (E 0 ), and river flow (Q) in the Merguellil Wadi.
The watershed is semiarid, characterized by very high spatial and temporal variability of precipitation. Figure 2 shows the annual precipitation of the nearby Kairouan station during the period 1950-2020, including the study period 1976-2017. The displayed long-term rainfall averages aim to assess the representativity of the study period in relation to long-term rainfall conditions. The average annual precipitation during the 70 years is 307 mm (standard deviation ± 108 mm). The average for the period 1976-2017 is equal to 290 mm (standard deviation ± 89 mm). Thus, it can be said that the study period is representative for longer periods. In Tunisia, the rainfall regime is often subject to a semi-periodicity of about 10 years [32]. In addition, since the 1950s, Tunisia has experienced a decrease of 5% per decade in annual rainfall [33]. The ten-year average precipitation ( Figure 2) shows that there are successive periods of sequential surplus, deficit, and periods close to the long-term average. Figures 3 and 4 show the annual and intra-annual variability of precipitation (P), temperature (T), potential evapotranspiration (E0), and river flow (Q) in the Merguellil Wadi.  The watershed is semiarid, characterized by very high spatial and temporal variability of precipitation. Figure 2 shows the annual precipitation of the nearby Kairouan station during the period 1950-2020, including the study period 1976-2017. The displayed long-term rainfall averages aim to assess the representativity of the study period in relation to long-term rainfall conditions. The average annual precipitation during the 70 years is 307 mm (standard deviation ± 108 mm). The average for the period 1976-2017 is equal to 290 mm (standard deviation ± 89 mm). Thus, it can be said that the study period is representative for longer periods. In Tunisia, the rainfall regime is often subject to a semi-periodicity of about 10 years [32]. In addition, since the 1950s, Tunisia has experienced a decrease of 5% per decade in annual rainfall [33]. The ten-year average precipitation ( Figure 2) shows that there are successive periods of sequential surplus, deficit, and periods close to the long-term average. Figures 3 and 4 show the annual and intra-annual variability of precipitation (P), temperature (T), potential evapotranspiration (E0), and river flow (Q) in the Merguellil Wadi.   Rainfall is often intense, especially in spring and autumn, which cause flooding [34]. The estimated peak flow during an exceptional flood event in September 1969 reached 3000 m 3 /s at Haffouz [35]. Periodic river flow of the Merguellil Wadi supplies the El Haouareb Dam that itself feeds the downstream aquifer of the Kairouan plain. The aquifer is heavily exploited for irrigated agriculture and drinking water supply. There are two phreatic aquifer systems in the watershed, the Bouhefna-Haffouz-Chérichira (BHC) and the Ain Beidha ( Figure 1). Little is known about their recharge and exchange, and links between river flow and the aquifer [36]. The soil texture can be classified into four main parts, namely alternating limestone-marl formations giving silty-clay soils with active clays conditioning active processes of soil degradation [37,38]. The geology is complex and mainly composed of limestone, marl-limestone, and marly sedimentary formations of the Cretaceous and Eocene [39,40]. The soil and water conservation practices consist of grass strips, stone barriers, contour ridges, and 44 lakes and hill dams. These structures occupy 260.45 km 2 or 22% of the watershed area [41].  Rainfall is often intense, especially in spring and autumn, which cause flooding [34]. The estimated peak flow during an exceptional flood event in September 1969 reached 3000 m 3 /s at Haffouz [35]. Periodic river flow of the Merguellil Wadi supplies the El Haouareb Dam that itself feeds the downstream aquifer of the Kairouan plain. The aquifer is heavily exploited for irrigated agriculture and drinking water supply. There are two phreatic aquifer systems in the watershed, the Bouhefna-Haffouz-Chérichira (BHC) and the Ain Beidha ( Figure 1). Little is known about their recharge and exchange, and links between river flow and the aquifer [36]. The soil texture can be classified into four main parts, namely alternating limestone-marl formations giving silty-clay soils with active clays conditioning active processes of soil degradation [37,38]. The geology is complex and mainly composed of limestone, marl-limestone, and marly sedimentary formations of the Cretaceous and Eocene [39,40]. The soil and water conservation practices consist of grass strips, stone barriers, contour ridges, and 44 lakes and hill dams. These structures occupy 260.45 km² or 22% of the watershed area [41].

Data Used
Rainfall and runoff data were used from the Skhira and Haffouz stations ( Figure 1). They were collected by the General Directorate for Water Resources [42]. Temperature (T) data were taken from the Tunisian Meteorological Institute [43] database for the Kairouan station ( Figure 1 and Table 1). Data gaps were less than 3% and filling of missing data was performed by multiple regression [44]. The daily data covered a 42-year   Rainfall is often intense, especially in spring and autumn, which cause flooding [34]. The estimated peak flow during an exceptional flood event in September 1969 reached 3000 m 3 /s at Haffouz [35]. Periodic river flow of the Merguellil Wadi supplies the El Haouareb Dam that itself feeds the downstream aquifer of the Kairouan plain. The aquifer is heavily exploited for irrigated agriculture and drinking water supply. There are two phreatic aquifer systems in the watershed, the Bouhefna-Haffouz-Chérichira (BHC) and the Ain Beidha ( Figure 1). Little is known about their recharge and exchange, and links between river flow and the aquifer [36]. The soil texture can be classified into four main parts, namely alternating limestone-marl formations giving silty-clay soils with active clays conditioning active processes of soil degradation [37,38]. The geology is complex and mainly composed of limestone, marl-limestone, and marly sedimentary formations of the Cretaceous and Eocene [39,40]. The soil and water conservation practices consist of grass strips, stone barriers, contour ridges, and 44 lakes and hill dams. These structures occupy 260.45 km² or 22% of the watershed area [41].

Data Used
Rainfall and runoff data were used from the Skhira and Haffouz stations ( Figure 1). They were collected by the General Directorate for Water Resources [42]. Temperature (T) data were taken from the Tunisian Meteorological Institute [43] database for the Kairouan station ( Figure 1 and Table 1). Data gaps were less than 3% and filling of

Data Used
Rainfall and runoff data were used from the Skhira and Haffouz stations ( Figure 1). They were collected by the General Directorate for Water Resources [42]. Temperature (T) data were taken from the Tunisian Meteorological Institute [43] database for the Kairouan station ( Figure 1 and Table 1). Data gaps were less than 3% and filling of missing data was performed by multiple regression [44]. The daily data covered a 42-year period (1976-2017).
The hydrological year begins on 1st September and ends on 31st August.  Table 2 shows descriptive statistics of the river flow, rainfall, evapotranspiration, and temperature series. The rainfall upstream in the basin (Skhira station) is significantly larger compared to the downstream (Haffouz station) due to orographic effects [45]. Figure 5 shows double mass curves for annual precipitation and river flow at the Skhira and Haffouz stations and approximate homogeneity of the data [46][47][48].   Table 2 shows descriptive statistics of the river flow, rainfall, evapotranspiration, and temperature series. The rainfall upstream in the basin (Skhira station) is significantly larger compared to the downstream (Haffouz station) due to orographic effects [45]. Figure 5 shows double mass curves for annual precipitation and river flow at the Skhira and Haffouz stations and approximate homogeneity of the data [46][47][48].  To analyze land use dynamics, we used multispectral Landsat 1, 5, and 8 satellite data. These images were used for calculation of NDVI index [49,50]. Landsat 1 MSS (Multi-Spectral Sensor) image was acquired for 21 July 1976, Landsat 5 TM (Thematic Mapper), and Landsat 8 OLI-TIRS (Operational Land Imager-Thermal Infrared Sensor) were for 19 July 1996 and 21 July 2016, respectively. The images have a spatial resolution of 30 m. All images were from July. This date coincides with start of the dry season and end of the rainy season. This made it possible to avoid effects of annual vegetation and to To analyze land use dynamics, we used multispectral Landsat 1, 5, and 8 satellite data. These images were used for calculation of NDVI index [49,50]. Landsat 1 MSS (Multi-Spectral Sensor) image was acquired for 21 July 1976, Landsat 5 TM (Thematic Mapper), and Landsat 8 OLI-TIRS (Operational Land Imager-Thermal Infrared Sensor) were for 19 July 1996 and 21 July 2016, respectively. The images have a spatial resolution of 30 m. All images were from July. This date coincides with start of the dry season and end of the rainy season. This made it possible to avoid effects of annual vegetation and to improve the knowledge on behavior and distribution of vegetation cover. Landsat images are downloaded free of charge from the USGS website (https://earthexplorer.usgs.gov (accessed on 3 October 2019)).

Runoff Trend Analysis
Identification of the flow trend was performed by using the ITA method (Innovative Trend Analysis) proposed byŞen [24,25]. It offers advantages over traditional statistical tools and its graphics present more detail by characterizing extreme events. In the ITA method, recorded flow data are divided into two equal parts from start to end of the time series. The two subseries are sorted separately in ascending order. For estimating the trend rate, the slope S ITA and trend indicator T ITA are calculated by: where n is sample size, X i , X j , X i ,nd X j are values and averages of the first and second half of the time series, respectively, µ is average of the first half of the time series.
The ITA method detects trends over a time series across two divided subseries of equal lengths. In addition, the standard deviation of the trend slope deviation formula assumes that variation is constant. This assumption is not valid for asymmetric and dependent series [51,52]. Therefore, ITA detects trends on a main time series with equal sub-series length [53]. Thus, we divided the 42-year study period into three periods : 20 years (1976-1995; 1996-2017), 30 years (1976-2005), and 40 years (1976-2017), and we verified that the series are not autocorrelated.
To implement the ITA method, the open source "trendchange" library was developed in R language and made available via the CRAN repository.

Climate Elasticity of Annual Flow
The climate elasticity of annual flow ε is an indicator of sensitivity of long-term hydrological systems to climatic fluctuation. The use of ε is informative for estimation of the impact of climate change on the hydrological system. It is a non-parametric estimator that calculates the elasticity value directly from time series data of annual flow in relation to the variability of precipitation or potential evapotranspiration [54,55]. Its evaluation offers an opportunity for efficient management of water resources [56,57]. It is calculated by: where Xnd Q are precipitation or potential evapotranspiration and average annual flow for the entire study period, respectively; X i and Q i are annual precipitation or potential evapotranspiration and flow for the ith year, respectively. The climate elasticity of annual flow, Equation (3), can be expressed according to the studies of Sankarasubramanian et al. [55] and Zheng et al. [54] in a non-parametric form: Thus, Equation (4) can be rewritten as: where ρ X,Q is the correlation coefficient between C X and C Q . These are the coefficients of variation for X and Q, respectively. Thus, these three parameters are calculated as: Potential evapotranspiration was calculated by the Thornthwaite method [58,59]: where E 0 (mm/month) is the average evapotranspiration of the month m (m = 1 to 12); T is interannual average temperature of the month in • C; I is annual thermal index, and F(m, ϕ) is a corrective factor depending on month (m) and latitude (ϕ). The parameters a, I and i are calculated as: 1.514 (12) F(m, ϕ) = L 12 N 30 (13) where N is the number of days in the month and L is the average day length (hours) of the month.

Connection between Climate Change, Anthropogenic Activity, and Streamflow
At watershed scale, streamflow can be simulated as a function of climate change (C) and anthropogenic activity (H) as [60]: where Q is flow; C represents integrated effect of climatic variables; H represents integrated effects of human activities. Therefore, the flow variation can be approximated by [17]: where ∆Q is total variation of streamflow; ∆Q C is change in streamflow due to climate change, and ∆Q H is change in flow caused by anthropogenic activities. The total variation of river flow is determined by: where ∆Q obs,1 and ∆Q obs,2 are the annual discharge observed in period P1 and P2, respectively. Precipitation (P) and potential evapotranspiration (E 0 ) are the dominant controlling factors for average annual water balances [6,15,18]. The variation in annual flow due to variation in annual average precipitation and potential evapotranspiration is determined by: where ε p and ε E 0 are the elasticity of precipitation and potential evapotranspiration of runoff. The ∆P and ∆E 0 are the variation in precipitation and potential evapotranspiration. The ∆P and ∆E 0 are determined by: where P obs,1 , P obs,2 , E 0,obs,1 and E 0,obs,2 are precipitation and potential evapotranspiration for the period P1 and P2, respectively.

Identification of Land Use Dynamics by NDVI Index
The normalized vegetation index (NDVI) was used to assess the importance of biomass and chlorophyll activity. The NDVI ranges from −1 to +1. Negative values correspond mainly to open water. Values close to zero are formed mostly of rock, roads, and bare soil. Moderate values (0 to 0.2) represent bare soil up to the early stage of vegetation growth, while values 0.2 to 0.8 correspond to plant variability. Values close to +1 (0.8 to 1) indicate forest areas or dense vegetation [61,62]. The NDVI index was used to assess the dynamics of land use in the Merguellil watershed. It was calculated by: where PIR is the infrared and R the red spectral band.

Analysis of Trend in River Flow
The    (e) (f) Figure 6. Trend analysis using the ITA method for different time steps (season, monthly, annual, and max and min) according to the four study periods for Skhira (a,c,e) and Haffouz (b,d,f).

Climate Elasticity of Annual Flow
To assess the impacts of climate change on the river flow of the Merguellil Wadi, the climate elasticity of the flow in relation to annual precipitation and potential evapotranspiration was calculated for the two periods using Equation (5) for the two stations Haffouz and Skhira. The is summarized in Table 5. The results of the climate elasticity for the streamflow vary from upstream to downstream of the watershed (Table 5). This translates into that a 10% change in precipitation implies a change from −6.2 to 16% in average annual flow for the Haffouz station and from 9 to 12% for the Skhira station during the two periods P1 and P2. A 10% decrease in potential evapotranspiration is associated with a −6 to 16.5% increase in streamflow (Table 5).  Like annual flow, the hydrological response of seasonal flow is much greater at the upstream than the downstream. The trend for period P1, P3, and P4 is greater in winter than during the other three seasons. The fall and summer seasons are characterized by a significant increasing trend during the 30-and 40-year periods at the upstream of the watershed. This trend is weaker at the downstream with a reduction of over 60%. Regarding the period P2, the trend is significantly decreasing.
On a monthly scale, in the upstream, September characterizes the fall trend. December and February play a primary role in the winter trend. The summer trend is influenced by July. In the downstream, January and February influence the winter trend. In addition, May and July characterize the trend of spring and summer, respectively.

Climate Elasticity of Annual Flow
To assess the impacts of climate change on the river flow of the Merguellil Wadi, the climate elasticity of the flow in relation to annual precipitation and potential evapotranspiration was calculated for the two periods using Equation (5) for the two stations Haffouz and Skhira. The ε is summarized in Table 5. The results of the climate elasticity for the streamflow vary from upstream to downstream of the watershed (Table 5). This translates into that a 10% change in precipitation implies a change from −6.2 to 16% in average annual flow for the Haffouz station and from 9 to 12% for the Skhira station during the two periods P1 and P2. A 10% decrease in potential evapotranspiration is associated with a −6 to 16.5% increase in streamflow (Table 5).  Table 6 illustrates the impact of climate change and anthropogenic activities on river flow in the Merguellil watershed. Indeed, by using Equation (17), in the upstream of the watershed, a decrease of 38 mm in precipitation results in a decrease of 0.17 mm in streamflow (Skhira). Thus, an increase of 109 mm in potential evapotranspiration will decrease streamflow depth by 0.03 mm (Table 6). Similarly, for the downstream of the Merguellil watershed, increase in precipitation and evapotranspiration results in a small increase in streamflow (Table 6). The cumulative effect of changes in precipitation and potential evapotranspiration results in a −5.2 to 0.6 mm increase in streamflow, i.e., a variation of 21.6 to 3.9% compared to observed variation (Table 6) at Skhira and Haffouz stations, respectively. Thus, it is estimated that anthropogenic activities account for more than 78% of the variation in streamflow for the Merguellil watershed (Table 6).

Impact of the Dynamics of Land Use on Streamflow
A decreasing trend in inflow to the El Houarreb Dam has been observed since the early 1970s. This decrease implies that the flow could be affected by anthropogenic activities more than climate change. In the present study, we assume that primary anthropogenic activity affects the flow of Merguellil Wadi through water withdrawal, land use, and land use change. The results of the calculation of the NDVI index for the Merguellil watershed for 1976, 1996, and 2016 are grouped together in Figure 7 and Table 7  During the study period, the Merguellil watershed underwent a change in land use. Indeed, the NDVI index is between −1 and 0.8, with an absence of values close to 1. The NDVI class between 0 and 0.2 has experienced a significant increase, while water bodies (NDVI < 0) have been reduced (Table 7). During the period 1976-1996, vegetation cover (0.2 < NDVI < 0.8) and open water bodies (NDVI < 0) increased by 13 and 0.3%, respectively. However, uncovered or bare surfaces (0 < NDVI < 0.2) decreased by 13% (Table 7). This is a result of conversion of denuded areas into forests, crops, and arboriculture as well as the construction of the El Houarreb Dam in 1989. However, during the period 1996-2016, the NDVI class between 0.2 and 0.8 experienced a decrease of 27%. At the same time, land covered with vegetation increased by 27%.

Discussion
Analysis of flow trend by the ITA method in the Merguellil watershed over the periods of 20, 30, and 40 years shows that the impact of time series length on the magnitude of trends is substantial. In fact, on a monthly scale a positive trend for P1 and negative for P2 for the 20-year period is not found in the 30 and 40-year periods, respectively. Indeed, at longer observation periods (≥30 year), the trend is increasing for annual flow. However, it is decreasing for extreme values of flow rates. In addition, one can notice a concordance of the magnitude of the trends during the four periods. This concordance appears, increasing in the watershed, particularly in spring during March and April. Downstream in the watershed, the concordance of magnitudes of trends characterized autumn, summer, and September, October, March, April, and June. The ITA method has  During the study period, the Merguellil watershed underwent a change in land use. Indeed, the NDVI index is between −1 and 0.8, with an absence of values close to 1. The NDVI class between 0 and 0.2 has experienced a significant increase, while water bodies (NDVI < 0) have been reduced (Table 7). During the period 1976-1996, vegetation cover (0.2 < NDVI < 0.8) and open water bodies (NDVI < 0) increased by 13 and 0.3%, respectively. However, uncovered or bare surfaces (0 < NDVI < 0.2) decreased by 13% (Table 7). This is a result of conversion of denuded areas into forests, crops, and arboriculture as well as the construction of the El Houarreb Dam in 1989. However, during the period 1996-2016, the NDVI class between 0.2 and 0.8 experienced a decrease of 27%. At the same time, land covered with vegetation increased by 27%.

Discussion
Analysis of flow trend by the ITA method in the Merguellil watershed over the periods of 20, 30, and 40 years shows that the impact of time series length on the magnitude of trends is substantial. In fact, on a monthly scale a positive trend for P1 and negative for P2 for the 20-year period is not found in the 30 and 40-year periods, respectively. Indeed, at longer observation periods (≥30 year), the trend is increasing for annual flow. However, it is decreasing for extreme values of flow rates. In addition, one can notice a concordance of the magnitude of the trends during the four periods. This concordance appears, increasing in the watershed, particularly in spring during March and April. Downstream in the watershed, the concordance of magnitudes of trends characterized autumn, summer, and September, October, March, April, and June. The ITA method has certain advantages for flow series analyses. The approach allows for more detailed interpretation of trend detection, which is useful for identifying hidden trends over time. The ITA method identified statistically significant upward and downward trends over the 30-year period. In this framework, we observed a high degree of connection of magnitudes of the trends for the 20 and 40-year series.
The 20-year trend magnitude results show a negative sign that is statistically significant. This is mainly due to the deficit years experienced in the study area ( Figure 2). Consequently, these results reflect the hydrological variability of the Merguellil watershed and not the hydrological trend linked to climate change. The 30-year series analysis identified a clear trend. Similarly, the 40-year series showed rather moderate amplitudes ( Figure 6). In conclusion, the 30 and 40-year series showed pronounced and moderate impact of climate change, respectively.
Statistical results of the study show that the hydrological response of stream flow is much greater upstream than downstream. The last three decades illustrate a 40% reduction in inflow to the dam. This is in accordance with what has been published in the yearly hydrological reports by the Kairouan hydrological department. Indeed, the Merguellil Wadi recorded zero flow for the first time in 1979 at Haffouz. Thus, a clear downward trend is observed between 1976-1977 and 1988-1989. This follows a succession of years of deficit due to the reduction in rainfall events of more than 30 mm [62]. In addition, despite the surplus period over five consecutive years from 2002-2003 to 2006-2007, an absence of decrease in rainfall over the watershed and across Tunisia was recorded [13,[63][64][65][66]. Thus, the Merguellil Wadi remains in a vulnerable situation.
Climate elasticity of annual flow expresses the percentage change in flow as a function of percentage change in climatic variables. It is proportional to precipitation, temperature, and increases with temperature if precipitation decreases and decreases with temperature if precipitation increases [67]. The results of the present study show that the climate elasticity of flow is variable. This lack of stationarity contradicts the conclusion of Sankarasubramanian et al. [55]. In addition, it can take negative signs (Table 5). This means that the reduction in input is combined with an increase in temperature, which tends to increase evapotranspiration, followed by a decrease in precipitation [68,69]. Therefore, the availability of water in the study area will likely be more critical in the future with current projections of global warming. An increase of 2 • C is likely to be followed by a significant decrease in water availability by 2 to 15% [70].
Our results suggest that the effect of climatic factors on the average annual flow increases from upstream to downstream of the watershed. This agrees with the results of Tsai [71]. Indeed, it is quite natural that our estimates of the effects of precipitation and evapotranspiration on the average flow in the Haffouz station differ from those estimated in the north of the watershed, at the Skhira station. A forest cover and an altitude, which varies from 700 to 1200 m, characterize the north. In addition, the effect of temperature on low water flow in the Haffouz region, an area characterized by a permeable soil texture, is greater than that on the average flow in the Skhira station.
The elasticity of potential evapotranspiration throughout the study period for the upstream is higher than downstream of the watershed. This indicates that the response of the Merguellil Wadi is more sensitive to potential evapotranspiration than precipitation. Certainly, evapotranspiration through vegetation significantly affects the flow, especially at the end of the growing season [72]. In general, the effects of non-climatic factors include changes in land use and human activities have been shown to be influencing factors for water supply. Vegetation changes in the Merguellil watershed (Figure 7) come about because of government intervention in the Kairouan region. This intervention aims to improve the country's self-sufficiency in certain products (vegetables, fruits, and cereals) and the supply of agricultural products to feed the livestock [73]. Thus, an extension of the cultivated area following governmental investments in water saving equipment linked to the development of drip irrigation between 1999 and 2005 [74].
Although, the results of the combination of the effects of climate change and land use show a reduction in water supplies to the basin. Human contribution plays an important role for the Merguellil watershed [22]. The growing demand of water for agricultural, domestic, and industrial purposes has led to overuse of surface and groundwater. The withdrawal of surface water will reduce the flow directly. Thus, the overexploitation of groundwater leaves a void to be filled by surface runoff, thus reducing the flow of watercourses.
Soil texture plays a role in the evolution of runoff. Indeed, upstream of the watershed (Skhira station), the riverbed of the Merguellil Wadi is formed by hard rocks with a stable section. Downstream, at the Haffouz station, the Merguellil bed is wide, formed by unstable sandy formations. It is then exposed to erosion phenomena, which means that the section changes shape almost at each passage of a flood [34]. Water for domestic, irrigation, industrial, and tourist use is being taken from deep groundwater of the Merguellil watershed. In addition, the disturbance of the hydrological balance in the study area can be explained by feedback from local residents following the disappearance of the low flow of Wadi Merguellil and its drying up since 1981 [34,75]. This consequence appeared following a political strategy to draw down the groundwater system of Bouhefna, located in the center of the study area (Figure 1), since the mid-1970s. This strategy aims to reduce water losses by evapotranspiration and for hydrogeological purposes. It therefore induces excessive hydrodynamic exploitation of underground resources [34]. This is planned in the expansion of agricultural pumping along the Wadi, by abandoning traditional irrigation canals, installing motor pumps, and setting up deep boreholes in the alluvium of the Wadi or wells on the banks from the 1980s. The transfer of water from Bouhefna to the coastal regions presents another explanatory factor for the hydrological changes experienced in the Merguellil watershed that is widely contested by the local population [21].
In a semiarid region such as the Merguellil watershed, the hydrological processes are characterized by an extreme spatiotemporal variability (e.g., [76]). The high rainfall variability requires human intervention to mitigate the landscape change. We particularly mention soil and water conservation practices and careful land use/land cover management. Indeed, surface runoff is generated by exceptional events that are rare but complex [77,78]. However, only extreme events produce significant runoff.
The expansion of water and soil conservation management has led to a marked decrease in runoff in the Merguellil watershed [41]. Many investigations have focused on the variability of runoff coefficient (e.g., [79][80][81]). In this regard, by analyzing 114 comparable events over the period 1989-2010, Ogilvie et al. [82] estimated a 40% reduction in the runoff coefficient following the expansion of water and soil conservation practices since 1996 with a decrease of 25% in annual flow at the basin outlet. Lacombe et al. [83] estimated that the sharp drop in flow arriving at the El Haouareb Dam over the past 40 years is due to a reduction of runoff generated by events below 40 mm from 28 to 32% and a reduction of 41-50% for all events over the period 1989-2005.
In conclusion, this research shows the importance of the flow to identify the impact of climatic change and anthropogenic activities in a Tunisian semiarid context. Thus, it is recommended that the length of the hydrological series should not be less than 30 years when analyzing the impact of climate change on streamflow. This can therefore serve as a benchmark for future research. In addition, this study confirms the importance of using seasonal and monthly flow data in trend analysis, rather than simple annual averages which may sometimes not reveal the full spatiotemporal complexity of flow trends.
Ultimately, this study identifies the importance of the continuity of research, which has been initiated in the Merguellil watershed for more than 50 years. This continuity has made it possible to show that anthropogenic effects make water scarcity a permanent situation in a semiarid context. The purpose of this is to guide rural development agents to consider water management and saving models adapted to climate change and land use to avoid finding themselves in similar situations in the future.

Conclusions
This study assessed impacts of climate change and direct human activities on streamflow in the Merguellil watershed. Trend at monthly, seasonal, and annual time scales was studied using the ITA method over the period 1976-2017 divided into periods of about twenty years each (P1: 1976-1995; P2: 1996-2017), about thirty years (P3: 1976-2005), and about forty years (P4: 1976-2017). The purpose was, on the one hand, to provide a length of series capable of informing us about the impact of climate change and on the other hand, allowing us to compare the hydrological contexts between the different study periods. Analysis of the trend shows that the hydrological response of flows is much greater upstream than downstream. The average annual flow towards the outlet of the watershed, the El Houarreb Dam, has decreased significantly over the past thirty years by about 40%. Extreme low flows have increased in high and medium altitudes, thus constituting a constant pressure on the availability of water resources in the watershed where the residents depend on agricultural activities. Strong points of this study are the precision of the detection of long-term trends. A 20-year time series reflects the hydrological variability and the characteristics of the watershed. This observation time, however, does not reflect the hydrological trends linked to climate change. Indeed, only series of 30 years and more can clearly and explicitly identify the impact of climate change. From the 40-year data, we can identify climatic impact that is more objective for an optimal management of water resources. In addition, this study confirms the importance of using seasonal and monthly flow data in trend analysis. Indeed, simple annual averages cannot reveal the spatiotemporal complexity of flow trends. The calculation of the NDVI index by using the three Landsat images provided an overview for a global qualitative analysis of the contribution of human activities on streamflow variation in the watershed during the study period. The spatial distribution of land cover for the three dates considered, namely 1976, 1996, and 2016, indicates an increase in cultivated land of 14% during the 40 years of the study. The elasticities of precipitation and potential evapotranspiration throughout the study period, for the upstream of the watershed are 1.26 and −0.26, respectively. While in the downstream of the watershed, they are 1.35 and −0.35, respectively. This fact indicates that the response of the Merguellil Wadi is more sensitive to precipitation than to potential evapotranspiration. The direct impact of human activities on water supplies in the Merguellil watershed is responsible for more than about 55% of overall river flow change. Climate change contributes less than about 22%. Therefore, the study area is subject to intensive human impact due to social and economic development for more than 50 years. Direct human activities are the main factor in the reduction of river flow. Even more, these results show the impact of historical context in terms of anthropogenic effects.
Ultimately, the flow remains a variable and a key indicator for decision-makers to consider for its spatiotemporal variability. Therefore, it is recommended that the hydrological observation times should not be less than 30 years when analyzing flows. This allows restoring the impact of climate change and those of human activities in Tunisia. This infor-mation is essential for the sustainable management of water resources. In addition, it can give implications for analytical decisions on interpretations of the hydrological response and on planning and development both upstream and downstream of the watershed.