Characterisation of Hydrological Drought and Implications for Sustainable Water Resources Management in the Sokoto-Rima River Basin (SRRB), Nigeria

Hydrological drought is usually characterised by water loss over time from both underground and surface supplies. Thus for this study, the assessment of hydrological drought was carried out by employing Cumulative Rainfall/Streamflow Anomaly as preliminary tools for the presence of drought signatures while detailed characterisation was via Streamflow Drought Index (SDI). The results revealed that hydrological drought was observed in all the stations; however, though in general, the stations could be classified as experiencing near normal drought conditions with mild drought signatures. The findings also revealed that the average streamflow deficit volume and durations of the hydrological drought severity were 1.780 Mm 3 and 192 months, 1.444 Mm 3 and 252 months, 3.148 Mm 3 and 252 months, and 0.159 Mm 3 and 372 months for Bakolori, Goronyo (pre dam construction era), Goronyo (post dam construction era) and Zobe stations, respectively. The results also revealed the relevance of flow duration curve and analysis of frequency of drought state transition for the development of scenario-based basin water resources management protocol. The coefficient of determination (R 2 ) statis-tic of the developed regression models indicate that 73.3% and 86.5% variation in streamflow dynamics across the Basin can be explained by climate change variables. However, for sustainable management of water resources in the Basin, it is hydrological risk. Hence, a robust framework that finds application both for planning mitigation actions which embody strategic, tactical and emergency components should be designed; to this end, analysis of persistence and recurrence of drought in time and determination of possible recurrent patterns are necessary.


Introduction
Drought is a stochastic natural phenomenon that has great impacts on the economy, society and environment of any nation. Thus its determination, monitoring and characterisation are of great significance in water resources planning and management. It is a complex natural hazard which is a composition of many factors such as hydrological, meteorological and agricultural; all are highly influenced by lack of precipitation amount and its frequency. But as reported by Ali et al. [1], hydrological drought occurs when dry weather patterns outweigh other climatic conditions. This calls for routine assessment of water resource availability to meet demands coupled with evaluations of water shortage risk through the analysis of different elements of the water supply system. Hydrological drought is determined accordingly by Van Loon and Laaha [2] in terms of the propagation of meteorological drought through the terrestrial hydrological cycle. This means that it is the significant decrease in the availability of water in all forms appearing in the land phase of the hydrological cycle [3]; by extension, a hydrological drought episode is related to stream flow deficit with respect to normal condition [3]. Since the definition of drought is dependent on the objective of a particular study, which is perhaps important in an attempt to quantify drought, it then becomes imperative to look at hydrological drought as a broad term related to negative anomalies in surface and sub-surface water. But if hydrological drought is framed as a natural hazard, terms for the hazard literature such as vulnerability and resilience are employed. On the other hand, hydrological drought can be seen as a water resources issue which focuses more on the imbalance between water availability and demand [4].
According to Mishra and Singh [5], Pozzi et al. [6], hydrological drought deserves more attention due to its crucial link with drought impacts. In many cases, hydrological drought period leads to shortfall in hydro-power generation. In Nigeria, shortfall of power generation from Kainji dam due to drought has led to power rationing with many cities in black out. Similarly, according to IPCC report as in Sceneviratne et al. [7], there is need to pay more attention to the space-time development of hydrological drought. Against this backdrop there-

Study Location and Hydrometeorology
The Sokoto-Rima Basin is located within Latitude 10˚N and 14˚N and Longitude 4˚E and 9˚5'E, respectively. The entire Basin is classified as belonging to Hydrological Area One (HA.1) according to Nigerian Hydrological Service Agency (NIHSA). It belongs to the semi-arid region of the country with pre-dominant Sudan vegetation. The average annual rainfall is between 364 -970 mm. The River Basin is majorly serviced by four dams (i.e. Goronyo, Bakolori, Zobe, and Jibia); largely for water supply, irrigation and domestic use by the riparian communities. Figure 1 shows the gauging stations that were considered for the study; this includes the network of the stations and the dam sites.

Determination of Drought Threshold
To properly appreciate the spatial and temporal pattern of the drought regime, the Variable Threshold Method was employed; that is, because of seasonality implication, seasonally varying truncation levels against the traditional constant threshold approach was adopted. This approach was adopted against the backdrop of the fact that it does not require a prior knowledge of probability distribution; by extension, this method directly produces drought characteristics (e.g. frequency, duration, and severity).
According to the Threshold Level Method (TLM), a drought is observed once the variable of interest X (e.g. Streamflow) is equal to or drops below a pre-defined threshold. For this study, this threshold was defined from the streamflow percentile statistics. Based on the recommendation of Huijgevoort et al. [8], the 20th percentile, also known as the 80th exceedance percentile was adopted. To do this, the entire streamflow series of monthly series was arranged in ascending order; that is, from the smallest to the largest. The monthly data points were converted into their respective percentiles:

100
where n = total number of observation, x i = the pith percentage of the data set.
Due to the possibility of seasonal variability in the streamflow series, thresholds were defined for each month; i.e. P 20 , T,

Hydrological Drought Characterisation
In line with the thoughts of Nalbantis [3], various hydrological variables are used to describe the availability of water in all its forms in the land phase of the hydrological cycle but the most significant variable from the view point of quantity of water is streamflow. This connotes that a hydrological drought episode is related to streamflow deficit with respect to normal conditions. This forms the basis for the employment of the Streamflow Drought Index (SDI) for the hydrological drought analysis in real-time. Sequel to the proceeding section, exclusive use of streamflow was adopted here. Towards this end, as in the classical approach in treating time, successive non-overlapping time intervals were used. In this study, time was treated as follows: April, the first was considered the begin- where i denotes the hydrological year and j the month within the hydrological year (j = 1 for April and j = 12 for March). That is, where V i,k is the cumulative Streamflow volume for the i-th hydrological year and the k-th reference period. Using the cumulative streamflow volumes (V i,k ) as the platform, the Streamflow Drought Index (SDI) was defined for each reference period k of the i-th hydrological year as Open Journal of Modern Hydrology , , where k V and S k are respectively the mean and the standard deviation of cumulative Streamflow volume of reference period k; these were estimated over a long period of time. This means that the truncation level was set to k V . As reported in Nalbantis [3], Streamflow may follow a skewed probability distribution, thus a two parameter log-normal distribution in a simple normalisation exercise was done; in this case, natural logarithms of the streamflow for the stations were taken accordingly. This was implemented according as are the natural logarithms of cumulative streamflow with mean k Y and standard deviation, yik S .

Frequency of State Transition
Based on the historical streamflow series, SDI was computed; this yields a series of drought states. The series , i k SDI of the SDI index was then calculated by employing Equation (4) and from this, the series of states where n m,k is the number of occurrences of state m in reference period k within the available sample of N year.
Using the features vectors of Table 1, the frequency of state transition , , m m k F ′ from state in reference period k to state m ' in reference period k + 1 was computed as in Equation (7)

Development of Flow Duration Curves and Regression Models
To account for the variability in the nature of the flow trend regime for the sta-

Rainfall/Streamflow Anomaly Analysis
In view of the fact that the rainfall anomaly directly measures the shortage of rainfall, and is basically quantified as the difference between the observation and the long-term climatological mean, both discrete and cumulative anomalies on annual and seasonal/periodic basis were done for both rainfall and streamflow.
Since the discrete anomaly approach is typically primitive, the cumulative anomaly was thus employed along to complement the short comings of the former.
In specific terms, the deviations of the periods or monthly measurements from long-term monthly average were tallied. But since the relative importance of the cumulative precipitation anomaly depends largely upon the magnitude of the anomaly in relation to normal conditions [11], to account for this effect, Foley (1957) approach as reported in Keyantash and Dracup [11] was adopted. This requires that each anomaly was normalised with the cumulative anomaly being expressed in thousandths of the annual rainfall or streamflow.

Hydrological Characterisation and Frequency of State Transition
As expressed by Nalbantis [3], each drought event is characterised through four attributes: 1) its severity expressed by a drought index, 2) its time of onset and duration, 3) its areal extent, and 4) its frequency of occurrence, the most significant variable from the view point of quantity of water is the streamflow. Hence, a hydrological drought episode is related to streamflow deficit with respect to normal conditions. In the application of the SDI approach, a set of seasonally varying values of the truncation level was adopted; in this case, the periodic means.
M. Y. Otache et al. Table 2 presents the computed values of the SDI for each station. From the results, it is obvious that non-drought conditions dominate with minor traces of mild drought events. Based on probability of each state condition, using the SDI approach as an exclusive methodology for assessing the hydrological drought condition may not really augur well as seen here. It is glaring that the entire basin exhibits non-drought conditions as evidenced by the results using the exemplar stations even though extensive periods of drought are noticeable in some stations. Figure 2 and Table 2, Table 3 show the frequency values of different drought conditions and the state transitions for the respective stations. It is imperative to note that in order to really appreciate the true state of hydrological drought incidences, longer accumulation reference periods were considered as noted in Table 2, Table 3. The results showed clearly that the occurrence probability of non-drought in all the stations is staggeringly higher relatively. However, traces of mild drought are discernible ( Figure 2 and Table 2, Table 3). This is more so considering the fact that in the overall, the state transition frequency does not show significant dispersions.
Within this general context, it suffices to submit that arbitrary classification of drought states as espoused by the SDI methodology creates problems as the nature of a particular flow regime could be intermittent, ephemeral and perhaps perennial. This becomes critical especially as noted in the findings of Nalbantis [3], some drought states according to the definitions of SDI could become irrelevant since cumulative streamflow will always definitely have positive value.
Despite this though, as noted, the matrix of frequency of state transition allows for robust predictive capacity of drought. Against this backdrop therefore, since hydrological data length and its associated integrity are vital elements, cognisance should be taken to ensure that monitoring of drought should be accomplished based on priorities of sensitive areas rather than simple drought occurrences.

Implications and Strategy for Sustainable Water Resources Management
From the available record at a particular site for a certain period, the trends in streamflow can be detected. Three devices used to study the variability of the flow that have direct applications in resource planning are frequency curve, mass curve, and duration curve. The primary objective of a flow duration curve is to ascertain how often flow of a given magnitude occurred during the period of record. In a statistical sense, a duration curve represents the cumulating of the frequency distribution curve; it can be of daily flows, mean monthly flows, or mean annual flows; thus employing either of these strategies provides an appraisal platform for flow value analysis at various dependability. In the light of this, Figure 3 and  50  75  100  125  150  175  200  225  250  275  300  325  350  375  400  425  450  475  500  525  550   0  10  20  30  40  50  60  70  80  90 1953 1957 1961 1965 1969 1973 1977 1981 1985 1989 1993 1997 2001 2005 2009 Anomaly (  The physical implications of the findings generally hinge on the type of indices employed and absence or otherwise of robust criteria vis-à-vis robust policy framework for routine assessment. In this context, concern is more on how to deal with suitable choice of index to work with; it is difficult to deal with a variety of distinct characteristics and impacts of hydrological drought. For instance, duration of zero-flow periods, streams never fall dry and even failure of models to capture appropriately low flows sufficiently. For this particular case under discourse, there is implication for defined policy adoption; there is no seeming physical policy for adoption to recover from extreme drought event associated with lack of proper definition of territorial unit to be considered for different levels of planning and relevant institutional interventions. Against the backdrop of the concerns in the preceding section, for sustainability of long-term management design, recourse should be to ensuring that: 1) ecosystem based management as well as ecological restoration be considered, 2) routine characterisation of drought regime encompasses an evolution of the temporal and spatial distribution of drought within the basin as well as 3) designing policy to reflect the understanding that support is required along the value chain not just for production support but towards reducing the water dependence of added value.

Conclusion
Based on the method employed in evaluating the hydrological drought and its characterisation, all the stations: Bakolori, Goronyo, Jibia, Zobe, Katsina and Gusau, Bakolori experienced near normal drought classification conditions on the average though with minor traces of mild hydrologic drought conditions. The average deficit volume (severity of drought) and drought durations stood at 1.780 Mm 3 and 192 months, 1.444 Mm 3 and 252 months, 3.148 Mm 3 and 252 months, and 0.159 Mm 3 and 372 months for Bakolori, Goronyo (pre dam construction era), Goronyo (post dam construction era) and Zobe stations, respectively. In general, it suffices to note that the application of FDC is a veritable construct for flow analysis, especially for reservoir management in this regard. The developed regression models provide for multiple scenario basin water resources development and management frameworks as typified by the explana-  1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 Normalised Anomaly Water Year tory power of the R 2 along with the frequency of drought state transitions. Worthy of note is the adoption of other forms such as Variable Threshold methods and its primary form or its variants; this, sometimes might lead to overestimation of the number of drought events. However, the assessment of hydrological drought plays an essential role in sustainable water resources management considering the fact that recovery from it is usually prolonged in view of recovery time of river systems. Thus, it is apt to routinely evaluate the conditions of drought; in this regard, hydrological drought characterisation should encompass strategic, tactical and energy levels by detailing spatial components such as synchronicity, clustering and breaking up of clusters of droughts. In this regard too, it suffices to note that due to the fluctuations of the drought behavior as a result of natural climate variability, the drought characteristics derived from a streamflow record fluctuate from year to year. Thus, it is therefore imperative to analyse the persistence and recurrence in time. Similarly, since droughts are regional events, routine statistical analysis is needed to assess the variability within affected area and thus by extension determine the possible recurrent patterns.