Comprehensive study of the wet and dry spells and their extremes in the Mediterranean climate basin Northern Tunisia

In agriculture, the characterization of dry spells is essential whether it is to calibrate the water needs of crops or the flow rates of rivers. This study seeks to develop a discretization of dry and wet spells on a monthly scale while evaluating the risk of extremes using the renewal wet-dry spell model. This model consists of defining the wet spell according a negligible precipitation threshold. The structure of the model is that all parameters of the climate cycle, including its length, are random variables. To study the trend of the parameters we use the Mann–Kendall test, while the magnitude is evaluated by the Sen’s estimation method. The approach is applied to Ichkeul Lake basin in northern Tunisia to demonstrate its capacity. This region is of great agricultural and water importance, although it holds six large dams. The results show that the duration of the dry and wet spells reach’s, respectively, 49 days and 17 days. The maximum dry spell was 49 days in 1982. The Mann–Kendall test revealed three stations with significant positive trend of the monthly extreme dry spell length (at March) located in south and east of the basin. The trend analysis of the seasonal rainfall number showed one station with significant negative trend in east and one station with significant negative trend in the center of the basin. Results indicated that no significant changes in the start and end of rainy season have occurred over the past years. But a great relation with a subsequent length exists. The results of this research assist farmers and managers in establishing drought management plans. It allow, among other things, to calibrate simulation models for a more realistic management of water reservoirs. It also makes it possible to plan irrigations on a more different basis from that of observations made at regular time intervals. We analyze the trends of the drought in Ichkeul lake basin, Northern Tunisia, characterized from the daily rainfall data of five stations. The alternating wet-dry spell model and a precipitation threshold value are used to define the rainfall/dry event. The Mann–Kendall test and the Sen’s estimation method were used to analyze the possible trends and the magnitude of variables analyzed, respectively. The results show an increasing trend of maximum monthly dry spells. No significant changes in the start and end of rainy season have occurred over the past years. These analyses provide useful information for science and society and make it possible to minimize unexpected damage due to long dry spells and to have effective and efficient planning for various stakeholders. We analyze the trends of the drought in Ichkeul lake basin, Northern Tunisia, characterized from the daily rainfall data of five stations. The alternating wet-dry spell model and a precipitation threshold value are used to define the rainfall/dry event. The Mann–Kendall test and the Sen’s estimation method were used to analyze the possible trends and the magnitude of variables analyzed, respectively. The results show an increasing trend of maximum monthly dry spells. No significant changes in the start and end of rainy season have occurred over the past years. These analyses provide useful information for science and society and make it possible to minimize unexpected damage due to long dry spells and to have effective and efficient planning for various stakeholders.

• The results show an increasing trend of maximum monthly dry spells. No significant changes in the start and end of rainy season have occurred over the past years. • These analyses provide useful information for science and society and make it possible to minimize unexpected damage due to long dry spells and to have effective and efficient planning for various stakeholders.

Introduction
It is well known that short or prolonged drought events, intense floods are often associated with adverse impacts either on society or the environment [1]. In the scientific literature search, there are several definitions for drought. A dry event can be defined as a series of consecutive dry days with no precipitation or less than a certain threshold value. This threshold can be defined according to physical considerations or by performing a statistical rainfall analysis [2]. The planning, development and management of water resources are affected by climate change, so today the issue of minimizing the adverse effects of drought has become very important [3]. Reiser and Kutiel [4] have studied the annual rainfall regime in the Valencia and Larnaca stations of the Mediterranean basin. They looked at several parameters, such as: total precipitation, date of accumulated percentage, length of rainy season, number of rainfall spells, rainfall spell yield, ratio between average and median rainfall, the relationship between accumulated precipitation and the number of rainfall spells and the dry spells distribution. They concluded that precipitation is totally absent during the summer in Larnaca, which better meets the definition of a Mediterranean climate than that of Vamancia. Caloiero and Coscarelli [5] studied the mean and maximum lengths of dry/wet spells, the number of dry and wet spells with different lengths, and the contribution of dry/wet spells with different lengths to the total number of dry/wet days in the Calabria region (southern Italy), over a period of 56 years , considering both yearly and seasonal scales. The results showed a similar behavior of dry and wet spells in the spatial distribution of mean lengths, and may be useful for water resource management (concerning dry spells) and risk reduction (regarding the wet spells). Raymond et al. [6] investigated the evolution of the long dry spells in the Mediterranean Basin during the rainy season and for regional future climate projections of the Mediterranean and European Coordinated Downscaling Experiment initiatives. The main results are an increase in the number of very long dry spells and an increase of their mean duration and spatial extent. Rivoire et al. [7] compared the extreme dry spells, defined either with a fixed threshold for a dry day (1 mm), or with a time-varying threshold estimated from reference evapotranspiration (ET0), of the basin Mediterranean. Results indicated positive trends in ET0 during the summer months in particular. A seasonal behavior has also been observed: in winter, the distribution of extreme dry spells is similar if it considers a fixed threshold (1 mm) or a time-varying threshold defined with ET0. However, during the summer, the extreme dry spell durations estimated with a threshold of 1 mm are greatly underestimated compared to extreme dry spells calculated with ET0. By analyzing the decadal to centennial scale hydrological response of southeastern Spain drylands, Machado et al. [8] pointed out that the recurrence of dry phases has been higher since the beginning of the seventeenth century, compared to the frequency of continuous wet phases with marked annual variability. They also concluded that the trend for longer dry spells and increased interannual variability and irregular seasonal distribution of precipitation is considered essential in the modeling for this specific arid Mediterranean region. Interest marked by Salvati et al. [9] is to develop a simple drought severity index (DSI), obtained from only two input variables (rainfall and reference evapotranspiration), able to produce information on a detailed geographic and time scale. Using dry spell climatology they estimated duration, intensity and severity of drought on a daily basis over the test period 1988-2003 in Italy. They concluded that the most arid period in Italy was recorded during spring-summer 2003, when especially dry conditions occurred in northern and central Italy. DSI correctly indicates a significant increase of drought from March to July 2003 throughout Italy. It may be used to assess the agro-climatic conditions occurring in a certain location through only two input variables. It allows a synthetic description of drought episodes in terms of both dry spell length and water balance.
The study by Ben Mhenni et al. [10] investigated the characteristics of drought and its effect on vegetation and agricultural productivity in three areas of Tunisia between 1982 and 2011, by exploiting meteorological drought indices and soil moisture and identifying the most appropriate index for each area.
The results showed that the meteorological droughts were short and frequent, triggering soil moisture droughts that were long lasting and intense. The standardized precipitation index was identified as the best indicator of agricultural vegetation and droughts in the northern forest, while the Palmer drought severity index was the best in the central steppe and southern desert. The lagcorrelation analysis found that the response of greenery, vegetation and wheat productivity to droughts was more pronounced and had a significantly longer lag in the central steppe than in other regions.
The review of Loon [11] provides an overview of the current state of scientific knowledge on the characterization of hydrological drought. It attaches particular attention to the influence of climatic and terrestrial properties on the hydrological characteristics of drought. In addition, it highlights that the current use and usefulness of different drought indicators is highlighted and it discusses recent advances in drought monitoring and forecasting. It provides a summary of the research on hydrological drought projections for the future. In a succinct way, the link of the hydrological characteristics of the drought with the impacts and the problems related to its management was also dealt with. Finally, the author defined the challenges for future research on hydrological drought that concern international initiatives. During rainy seasons, success and failure in crop production dependent on frequency and length of dry spells [12,13]. Water deficit resulting from drought reduces crop yield because of its negative impacts on plant growth [13]. Changes in the magnitude and frequency of droughts due to climate change will have severe impacts on agriculture, especially crop production, cropping systems, and livestock [13]. Crop water needs differ according to species and climatic parameters [14]. It should be noted that the additional irrigation of crops during the rainy season depends on the distribution of the precipitation during the season and therefore the distribution and the occurrence of dry spells. The scientific literature is very rich with regard to the modeling, analysis and investigation of the drought phenomenon.. This study has three main objectives: (1) Identify and analyze the trend of events (dry/wet) under prior knowledge of the characteristics of dry/wet spells and demarcation of start and end of rainy season in northern Tunisia with Mediterranean climate; (2) Test analysis of the observed variability of regional drought occurrence at the catchment scale (3) Contribute for the depreciation of the minimization of damage due to long dry spells and to have effective and efficient planning for various stakeholders of this region.

Case study
Tunisia is one of the Mediterranean countries least endowed with water resources. Water resources in Tunisia are characterized by scarcity and a pronounced irregularity. As a result, water has become in Tunisia in recent decades, a decisive factor in spatial planning and economic and social development [15]. The most recent assessment of water resources is 4503 million m 3 of annual resources, of which 2700 million m 3 is surface water [16]. In order to respond to a growing demand, strategies and master plans for mobilization and use have been drawn by the water services of Tunisia. The main surface water mobilization works (large dams, collinear dams and collinear lakes) were built on the rivers of the basins located in the Extreme North of the country, which provides 81% of the surface water resources (Table 1). These basins are characterized by young reliefs, in geological terms, presenting valleys and steep ridges, dominated, on the one hand, by marl soils and soils vulnerable to erosion and, on the other hand, by a forest cover limited to the most watered parts and a seasonal vegetation cover.
Ichkeul lake basin is one among the tributaries of the Extreme North basins. The total length of the main River (wadi Sejenane) is about 60 km and its area is 2120 sq km. It bounded by Latitude 4,066,689 m to 4,108,664 m N and Longitude 510,000 m to 570,500 m E (Fig. 1). Five large dams location serving approximately 7000 ha of irrigated perimeters of fruit trees, vegetable crops, fodder and cereals is reported also in this map. These dams contribute also to water transfer from North to Center and Sahel of the Country, the floods protection of Mateur city and the hand of the ecological balance of Ichkeul Lake.
Ichkeul lake basin is more often under the Mediterranean influence, sub-humid and has mild winters [17]. According rainfall data of stations operated by General Directorate of Water Resources (Tunisian Ministry of Agriculture), on average the rainiest months for autumn season are November (96.9 mm); December (105.5 mm) for winter season; March (61.3 mm) for spring season, and June (12.1 mm) for summer season. The driest month is July (3.3 mm) (see Fig. 2). On a monthly scale, apart from the dry season, the highest average precipitation is recorded in December (105.5 mm) then in January (96.9 mm). On the other hand, the lowest precipitations are observed in May and September with, respectively, 27.1 mm and 44.2 mm. But the rains are much more complex. They appear quite suddenly in September, quickly reaching a maximum in October. From November to February a glut of winter, this can last until March exist. From May the precipitation descend rapidly, and preceded the summer shortage that will last until August (Fig. 2). About 43.6% precipitation occurs during December to February, 31.8% during September to November and 21% during March to May ( Table 2). The seasonal distribution of precipitation shows a relative predominance of winter season, but autumn and spring are also relatively rainy ( Table 2). Most of this basin Mean monthly temperatures range from 11 and 9.5 °C (in January) to 27, respectively; 27.5 °C (in July). Evaporation was measured using the Colorado evaporation pan installed at Ghèzala dam. The annual total evaporation is 1518 mm. It varies between 105 and 35 mm during September-October; 30-135 mm during March to May and 203-232 mm during June to August. Maximum evaporation (262 mm) is in July and minimum (37 mm) is in January. But since, precipitation is spread over almost 8 months; this area does not receive precipitation more than natural evaporation. Generally, annual evaporation is double the annual precipitation. During all months (except November to February), evaporation is more than the naturally occurring precipitation. On September and May, the deficit is quite appreciable from this count. The deficit will have to be partially made good by conserving water from ravine into collinear lakes at several places and by increasing mulching agricultural fields to prevent evaporation [12]. In particular, the building in Tunisia of small earth dams in underprivileged rural areas contributes to expansion of agriculture through  the development of irrigated areas and aquaculture, the watering of cattle. In addition, they contribute to the protection of large dams against premature siltation. Daily rainfall data of five stations (Table 3, Fig. 1) were collected from start activity until 2010. These stations were selected because they are the oldest in the region and are fairly evenly spread throughout the study region. These daily records were used in this subsequent analysis study.

Method
The methodology is based on the wet-dry spell approach. Dry spells consist of a series of dry days framed by wet events (Fig. 4). Therefore, the rainy season is defined as a succession of wet spells and dry spells [2]. The beginning of the first wet spell in autumn, in September, marks the beginning of the rainy season, whiles the end of the last wet spell in the spring (April) marks  The wet events themselves are defined as an uninterrupted series of rainy days comprising at least one day having received a precipitation above a predetermined threshold value (Fig. 4) [2]. The definition of event is associated with a rainfall threshold value which defines wet. The limit 1 mm day −1 has been selected because it corresponds approximately to the minimum daily evapotranspiration in Ichkeul lake Basin during the rainy season [18]. Thus, at least a daily precipitation will balance the minimum daily evapotranspiration.
H n,m is given by the following equation: hj stands for the daily rainfall totals in mm. Note that hj > 0 and that for at least one hj > 1 mm. The dry spell Z n,m represents the number of days between two subsequent wet spells. The length of the rainy season L n is defined as the time span between the onset of the first and the end of the last event of the given season; while the length of the climatic cycle is determined as the time elapsed between the onsets of two subsequent rainy seasons (Fig. 4).

Identification of wet and dry spells
As mentioned already a physic criterion has been developed and applied to identify wet and dry spells over the five rain gauge stations. Rainfall threshold used in the criterion is derived from local evaporation climatology that is daily mean evaporation (DME) over the area of interest [12]. Computational steps of the schemes are as follows: 1. Identification of dry spells (DS) and wet spells (WS) (Fig. 4). Rainy days with less than the threshold depth of precipitation are accounted for the wet spell if they occur in an uninterrupted sequence-Only isolated sub threshold rainfall will be discarded, and considered as part of a dry spell-A dry spell is assigned to month in which it started even if it continues in the following month or months; 2. Calculation of duration, and rainfall amount of wet spells; 3. Calculation of the number of wet spells per rainy season n (M n ); 4. Deduction of the number of dry spells per rainy season n (Z n = M n − 1); Deduction of lengths of rainy season and climatic cycle; and 5. Identification of starting and ending date of rainy season.
Besides these, characteristics of important parameters (duration, number) of dry spells and wet spells duration and total accumulated rainfall depth per wet event, and his extremes in respect of yearly number of events; starting and ending date of rainy season have been studied in order to understand if there is any change taking place in the hydrological cycle of the area. All the parameters of the two elements wet and dry spells of the rainfall time distribution have been examined in this study. Description of spatial variation of dry spells, and rainfall amount of individual wet spell is given across the region.

Trend analysis
Statistically, trend is a significant change over time which is detectable by parametric and nonparametric procedures while trend analysis of a time series consists of the Fig. 4 Representation of the climatic cycle. Legend: m: the index m refers to the mth event (spell) of the nth rainy season; D n,m : duration of mth wet (rainfall) spell in a nth rainy season symbolizing the number of subsequent rainy days; H n,m : total accumulated rainfall depth of D n,m rainy days in mm; hj: daily rainfall amounts; j: the index j refers to the jth day. Z n,m : duration of mth dry spell in a nth rainy season.
Climatic "annual" cycle (C n) magnitude of trend and its statistical significance [19]. In this study the Mann-Kendall test was used to detect trends in the series of variables analyzed, and then the magnitude was examined by the nonparametric Sen's estimator method.

Mann-Kendal test
The nonparametric test, usually known as Kendall's τ statistics, [20,21] is a widely used in hydrology and climatology to test randomness against trend of hydrologic time series [19,22]. This test is not affected by the actual distribution of the data and is less sensitive to outliers. Therefore, it is more suitable for detecting trends in climatic and hydrological data, which are usually skewed, and may contain outlier observations [22]. For any sample of n variables, x 1 ,…, x n , the null hypothesis states that the sample is independent and identically distributed. The alternative hypothesis of a two-sided test is that the distributions of x i and x j are not identical for all k, j ≤ n with i ≠ j [23].
The MK test is based on test statistic S defined as follows: , where the x j are the sequential data values, n is the length of the time series, and sign ( The mean E[S] and variance V[S] of the statistic S may be given as follows: where t p is the number of ties for the pth value and q is the number of tied values. The second term represents an adjustment for tied or censored data. The standardized test statistic (Z MK ) is computed by: A positive Z MK indicates an increasing trend, whereas a negative Z MK indicates a decreasing trend. To test for either increasing or decreasing monotonic trend at p significance level, the null hypothesis is rejected if the absolute value of Z is greater than Z 1−p/2 , where Z 1−p/2 is obtained from the standard normal cumulative distribution tables. In this work, the significance level of p = 0.05 is applied.

Sen's estimator method
We used Sen's nonparametric estimate of the slope to determine the magnitude of trends [24] of dry and wet spells time series data. The Sen's slope estimator method uses a linear model for the trend analysis. The slope (T i ) of all data pairs is calculated using Eq. (6) by: where x j and x k are data values at time j and k (j > k) separately.
The median of these n values of T i is represented by Sen's slope of estimation which is calculated using Eq. (7): Sen's estimator (Q med ) is calculated using above equation depending upon value of n is either odd or even and then (Q med ) is computed using 100 (1 − α) % confidence interval using nonparametric test depending upon normal distribution. A positive value of Q i indicate increasing (upward) trend while negative value of Q i represent downward or decreasing trend of time series data.

Dry event analysis
In planning studies the dry spells of maximum duration are of crucial importance. These values and number of dry spells recorded during rainy seasons were extracted from the event based analysis. Longest dry spell durations of five rain gauges for Ichkeul Lake Basin are given in Table 4. The number of dry spell ranges between 20 and 32. Hence, supplementary irrigation iterations are required. The values in Table 4 show that out of 29 dry spells, there was a 49-day occurrence at Fritissa rain gauge station during   Table 4, it was observed that maximum dry spell (49 days) at Fritissa in 1982 and minimum 38 days at Sidi Salem in 2008. In order to set the crop planting and harvest of the main crops in the basin, the crop calendar was developed based on field surveys and internal reports of Agricultural Extension Services ( Table 5). The knowledge of occurrence of the first wet spell in autumn is very important because it affects germination for autumn sowing. Already, the temporary interruption of the rainy season by a dry spell can damage crops, for example, if the dry spell occurs at the time of pollination or before the crops begin to ripen. Damage can occur even if the total seasonal precipitation is normal. Among others, dry spells have adverse consequences on the flow of the rivers. While consideration of past wet spells, the autumn season is postponed by one month (September to October) which is shown in crop calendar of Ichkeul lake basin ( Table 5). Shifting of the entire cropping season by one month is affecting crop yield for various crops.
Analysis of maximum monthly dry spells is summarized in Table 6 and showed in Fig. 5. There is no statistical evidence that the average duration of this phenomenon is different in February, January or March. The mean maximum monthly dry spell during the core of the rainy season varies between 10 and 15 days. It is interesting to note that the length of dry spell in October, January and February at Sidi Abdelbasset is greater than for the other rainfall stations in Table 6. In fact, this station is situated in semi-arid region of the basin. This phenomenon can come also from the bias due to a short observation period. It can be seen that the dry spell seem to show randomness throughout the rainy season from September to April.
The effects of climate on drought duration were studied by statistical analysis. Short droughts lasting less than 3 days have been discarded. Subsequently, the mean, maximum and standard deviation statistics were calculated for the variable duration of dry spells (Table 7).
By comparing the dry spells between stations, Table 7 shows a small difference between dry spells By examining the data in Table 6, it is concluded that the average drought duration are a first indication of the water deficit that plants may suffer from it, and that this drought duration might have different governing factors. We conducted a linear regression model analysis to identify the trends in the basin. Figure 6 shows the number of drought affected Ichkeul lake basin. The regression trend lines are also given in this figure. Overall, the number of dry spell in the study area was increasing except for Ghèzala-Dam (Fig. 6). The number of drought affected basin is reported more in Ghèzala Dam and few number of drought affected Sidi Salem station. Most number of droughts affected Ghèzala Dam is observed in 1982,1985,1987,1989 to 1991, and 1995 to 1997 over Ichkeul lake basin (Fig. 6).
According to data from the General Directorate of Water Resources in Tunisia, recently, severe drought has been experienced in this region during September 2015 to September 2018, and indeed throughout the country, where surface water resources were almost entirely committed to human consumptive uses, groundwater exploited intensively and the discharge to the sea almost nil. During the entire operation period of the stations, a significant increase in the severity of droughts was reported over Ichkeul lake basin during the period of 1960-1961, 1968-1969, 1994, 2015-2016 and 2016-2017. Analysis carried out on the dry spells number is conducted in order to reflect the dry spells severity in this area. A conditioned number-duration by splitting the dry spell duration into different classes was performed  1960 1963 1966 1969 1972 1975 1978 1981 1984 1987 1990 1993 1996 1999 [25]. In fact, the results analysis has shown that a short to medium duration year can have a significant number of dry spells, without denying that there can be dry spell of long duration. Thus, dry spells occur randomly during the rainy season Mathlouthi and Lebdi [25]. Overall, the annual number of dry spells in Ghèzala Dam significantly decreased, with a rate of 0.217/year (p = 0.001). In other stations (Fritissa, Wadi Joumine Antra, Sidi Abdelbasset, Sidi Salem) the trend is not significant at 5% (Table 8). Table 9 shows the Mann-Kendall statistics for monthly extreme dry spell length and related p-values. From this table, it is clear that three stations, Fritissa, Ghèzala Dam and Wadi Joumine Antra, have a positive trend at 5% significant level with, respectively, a rate of 3.53 days/10 years (p = 0.049), 1.67 days/10 years (p = 0.034) and 1.56 days/10 years (p = 0.034), in March included growing season. This may show that the monthly extreme dry spell length have changed in recent years over Ichkeul lake basin for three stations. It should also be noted that Fritissa and Wadi Joumine Antra stations are located, respectively, in South East semi-arid region and south west sub-humid regions of the basin.
In Table 9, two stations have significant trend, among which, one station shows a negative trend at 5% significant level with a rate of − 7.39 days/10 years (p = 0.005) in September (Fritissa), and one station shows a negative trend at 5% significant level with a rate of − 2.31 days/10 years (p = 0.048) in October (Sidi Abdelbasset).
This section presents the nonparametric Mann-Kendall test results for monthly dry spell length. Regarding Table 10, as already mentioned each dry spell is assigned to the month in which it started. From Table 10, it is clear that two stations, Fritissa and Sidi Salem, have a negative trend at 5% significant level on September. This may shows that the monthly dry spell length has not changed in September except for two stations. Sidi Abdelbasset station shows a positive trend at 5% significant level in October (autumn season) and February and March (growing crop season).

Wet spell analysis
The exploratory data analysis is based on the calculation of data descriptive statistics. Table 11 shows observation number, maximum, average, and standard deviation of the selected variables. The average rainfall duration of Ichkeul Lake basin is 2.19 days and the average rainfall amount is 21.01 mm (Sidi Salem). Maximum seasonal rainfall number are observed at Ghèzala dam (30 events) and Wadi Joumine Antra (29 events) stations in center region of the basin, whereas the minimum seasonal rainfall number (23 events) belong to Sidi Abdelbasset station in south semi-arid region of the basin.
In order to understand the adequacy of precipitation during the first to 39th meteorological week for rainfed cultivation, a weekly analysis of precipitation during the rainy season (September-April) was carried out (Table 12). Calculus show that 20 mm of rainfall in week will be able to meet 0.54-0.75 times the evaporation demand, at the onset or end of the growing season. Allen et al. [26] stipulate that during the early stages of growth, the crop water requirement will be about 0.7 of the evaporation demand and would increase during the reproductive stage of crop growth. Longest observed wet spells in Ichkeul lake basin are given in Table 13. The number of wet spells varies between 24 and 32. From Table 13, it is quite clear that the longest wet spells are related to the annual total precipitation. The longest rainfall event was recorded in Sidi Salem rain gauge in 1962; 17 days resulting in 122 mm. The shortest is observed in Sidi Abdelbasset in 2002 with low annual precipitation; i.e., 10 days with an accumulation of 63.5 mm. Table 14 reports the Mann-Kendall test results for the wet spells duration, rainfall amount and seasonal number during the rainy season. From this table, we can see that for the wet spells duration the trend is significant for Sidi Abdelbasset and Sidi Salem stations with a Sen's slope of 0 (p = 0.001 and < 0.0001, respectively). At Sidi Abdelbasset station, the trend of the rainfall amount per wet spell is significant, but the Sen's slope is null (p = 0.012). A negative trend for the seasonal rainfall number is observed in two stations (Ghèzala Dam and Sidi Salem stations) with respective slopes of − 0.217 and − 0.111. Table 15 shows the main statistics of rainy season and climatic cycle lengths. From this table, we conclude that the average length of the rainy season differs slightly among stations. About the maximum, Fritissa records 258 days  and Sidi Abdelbasset scores 241 days. It appears that onset or end of the rainy season can move several weeks. The arithmetic mean of the climatic cycle length (Table 15) confirms the annual character for this climatic parameter. However, extreme lengths exceeding 400 days can be observed (445 days at Sidi Salem station). Table 16 shows there is no tend for these two seasonal parameters at 5% significance level.

Start and end of rainy season
The variability of start of the rainy season is similar for pairs of stations Ghèzala dam, Wadi Joumine Antra and Sidi Abdelbasset, Sidi Salem (Table 17). There is no variation for the end of the rainy season at all stations. Concerning the average, the rainy season (or climatic cycle) starts at mid-September in Fritissa, Wadi Joumine Antra and Sidi Salem and ends earlier at the end of April (Table 17). At Fritissa, Wadi Joumine Antra and Sidi Salem the rainy season starts on 14 September and ends, respectively, on 30, 27 and 23 April (Table 17) (Table 18).
Being in the drier part of the transect, the wet season at Fritissa and Sidi Abdelbasset stations ends earlier than the other stations along the transect. This confirms the results from Mupangwa et al. [27] who showed that as aridity increases, the season starts late and ends early. However, the size of the observation series is much shorter for these two stations to be able to consolidate these results. The most brutal end of the season in Sidi Abdelbasset was recorded on April 9.
The latest end of the wet season was recorded in 1968/69 and occurred on May 5. This analysis shows, among other things, that the abrupt end of the rainy seasons observed is in agreement with the observations made in semi-arid southern Ichkeul lake basin during the 1982/83 rainy season. The 1982/83 rainy season ended on April 2, 1983 (in southern region). It is also noted that the start of the rainy season is more variable in Fritissa, Sidi Abdelbasset and Sidi Salem than its end (Table 17). In our study, the farmers of Fritissa are faced with the prospects of a shorter rainy season compared to those of the Ghèzala dam, as it happens that the wet season did not start until September 13 and ended early on March 24; either giving a duration of only 190 days.

Relationships of rainy season characteristics
The linear regression model was applied to establish a linear relationship between the different characteristics of the rainy season in Ichkeul lake basin. Correlation analysis indicated that there is no significant relationship between the start and end of the rainy season (Table 19). A great relationship between the start and the duration of the rainy season was found (Table 19). The gradient in the relationship between the start and the end of the rainy season is most pronounced at Sidi Abdelbasset station, which implies that a delayed start of the season is probably followed by an early cessation of the rainy season. This station is located in the southern watershed, where the climate is semi-arid.
In Sidi Abdelbasst, the linear relationship between the start of the rainy season and the length of the climatic cycle is significant. The linear relationship suggests that a longer climate cycle translate into an early rainy season. This result provides an additional tool for characterizing climatic zones.

Discussion and conclusions
A set of five rain gauges in Ichkeul lake basin Northern Tunisia with daily rainfall recording for at least 28 years has facilitated a detailed analysis of the chorological spatial patterns of trends of dry and wet spell magnitudes. The length of dry spell in October, January and February at Sidi Abdelbasset is greater than for the other rainfall stations. This station is situated in semi-arid region of the basin, and this phenomenon can come also from the bias due to a short observation period. By comparing the dry spells between stations, a small difference between dry spells number (around 5.3 per year) and average duration (around 1.2 days) of droughts is showed. The significant difference affected by the variability in precipitation is recorded between Sidi Abdelbasset and Ghèzala Dam. Sidi Abdelbaset is characterized by less dry spells but longer, while in Ghèzala dam, the highest number of short dry spells reflects greater rainfall variability. The longest drought occurs at Fritissa (49 days). The number of dry spell in the study area was increasing except for Ghèzala-Dam generated by a small variation in annual rainfall (CV = 0.25). The annual number of dry spells in Ghèzala Dam significantly decreased, with a rate of 0.217/year (p = 0.001). In other stations, the trend is not significant at 5%. Focusing on the monthly extreme dry spell duration, Fritissa, Ghèzala Dam and Wadi Joumine Antra stations have a positive trend at 5% significant level with,  respectively, a rate of 3.53 days/10 years (p = 0.049), 1.67 days/10 years (p = 0.034) and 1.56 days/10 years (p = 0.034), in March (growing season). Subsequently, the monthly extreme dry spell length has changed in recent years over Ichkeul lake basin for three stations. Also note that Fritissa is located in South East semi-arid region and Wadi Joumine Antra is located in south west sub-humid region of the basin. These results also indicated that for the analyzed time period, there was significant climate change in the sub-humid region of Tunisia. The average rainfall duration of Ichkeul Lake basin is 2.19 days and the average rainfall amount is 21.01 mm (Sidi Salem station). Maximum seasonal rainfall number are observed at Ghèzala dam (30 events) and Wadi Joumine Antra (29 events) stations with greater and more regular annual rainfall, whereas the minimum seasonal rainfall number (23 events) belong to Sidi Abdelbasset station with low annual rainfall. in south semi-arid region of the basin. The trend is significant for Sidi Abdelbasset and Sidi Salem stations with a Sen's slope of 0 (p = 0.001 and < 0.0001, respectively). At Sidi Abdelbasset station, the trend of the rainfall amount per wet spell is significant, but the Sen's slope is null (p = 0.012).
Weeks receiving rainfall about 20 mm are not able to meet 0.54-0.75 times the evaporation demand in September and April, respectively. During the early stages of growth, the crop water requirement will be about 0.7 of the evaporation demand Allen et al. [26].
The variability of start of the rainy season is similar to Ghèzala dam-Wadi Joumine Antra stations and Sidi Abdelbasset-Sidi Salem stations. On average, at Fritissa, Wadi Joumine Antra and Sidi Salem the rainy season starts on 14 September and ends, respectively, on 30, 27 and 23 April giving an average season length of 228, 225 and 221 days. It was observed that in Ghèzala dam the probability to have a rainfall event starting in the fifteen days of September is 0.48 with a return period of two years [25] and continues up to 39 th meteorological weeks. The wet season at Fritissa and Sidi Abdelbasset stations ends earlier than the other stations, given their location in the driest part of the transect, This confirms the results from Mupangwa et al. [27] and Edao et al. [28] who showed that as aridity increases, the season starts late and ends early.
Correlation analysis indicated that there is no significant relationship between the start and end of the rainy season. The gradient in the relationship between the start and the end of the rainy season is most pronounced at Sidi Abdelbasset station, which implies that a delayed start of the season is probably followed by an early cessation of the rainy season. This can be explained by the fact that this station is located in the southern watershed, where the climate is semi-arid. Here also, the linear relationship between the start of the rainy season and the length of the climatic cycle is significant. It suggests that a longer climate cycle translate into an early rainy season. This result provides an additional tool for characterizing climatic zones.
It was pointed out that these results depend on the threshold value of negligible precipitation [29]. These results reflect only the case of 1 mm per day.

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.