Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Label_Title

Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Label_Title

Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Use of L-moments in the Regional Frequency Analysis of Meteorological Drought Duration in the Transboundary Comoé River Basin

Ouattara Ismaïla*, Bamba Adama, Coulibaly Léréyaha, Dao Amidou, Kamagaté Bamory

Department of Mines and Reservoirs, Faculty of Geological and Mining Sciences, University of Man, Man, Côte d’Ivoire

Abstract: This study aims to estimate the duration of meteorological droughts in the Comoé River Basin. The sparse distribution of rainfall stations motivated the use of a regional frequency analysis based on L-moments. This method was applied to drought duration series derived from the 3-month Standardised Precipitation Index (SPI-3), using a drought threshold of -0.84. Classification and homogeneity tests identified four regions characterised by different drought patterns. Region R2 (Abidjan, Adiaké, Alépé, Bassam, Adzopé) represents the transitional equatorial climate in the southern part of the basin, while Region R4 (Tafiré, Tengréla, Ferkessédougou, Bobo-Dioulasso, Niangoloko, Kampti) is typical of the Sudanian climate in the north. The other two regions, R1 (Orodara, Ouangolo, Abengourou, Ouéllé, Dabakala, Bongouanou, Daoukro) and R3 (Bouna, Bondoukou, Agnibilekrou, M’Bahiakro, Aboisso), exhibit contrasting climatic characteristics because some rainfall stations that make up these groups are not located within the same climate zone. The Wakeby distribution was found to be the most suitable for modelling drought durations. Regional quantiles, estimated using the index flood approach, demonstrated greater accuracy than site-specific quantiles. These quantiles were subsequently interpolated spatially. The results indicate a relatively homogeneous pattern of drought durations across the basin, except in the far north. High probabilities of prolonged drought were observed in this northern area, which corresponds to a transitional tropical climate zone. For return periods (T) of 2, 5, and 10 years, the average drought duration was estimated at 2, 3, and 4 months across most of the basin, compared with 4, 6, and 8 months in the northern region. For T = 20 and 50 years, drought durations are expected to reach 7–8 months in the south, and up to 11 and 16 months, respectively, in the north.

Keywords: regional frequency analysis, L-moments, drought duration, west african monsoon, Comoé river basin

Received September 06, 2025; Revised October 08, 2025; Accepted October 15, 2025

Citation: Ouattara Ismaïla, Bamba Adama, Coulibaly Léréyaha, Dao Amidou, Kamagaté Bamory, Use of L-moments in the Regional Frequency Analysis of Meteorological Drought Duration in the Transboundary Comoé River Basin. American Journal of Water Resources, 13 (4). 158-169. January 2025. doi:10.12691/ajwr-13-4-5

1. Introduction

Drought is one of the most severe natural hazards, with far-reaching impacts on water resources, natural ecosystems, and agriculture. Frequent and intense droughts limit vegetation growth and increase soil vulnerability to erosion caused by heavy rainfall [3]. In developing countries, droughts are a major cause of famine, disease outbreaks, and land degradation, while in developed regions they result in substantial economic losses [4]. Unlike other natural hazards, drought is more difficult to detect, as it typically arises from multiple contributing factors that only become apparent after a prolonged period of deficit [5]. West Africa and Southern Africa are experiencing severe drought conditions that are disrupting agricultural and livestock production systems across nearly 14 countries, while agriculture in the region remains approximately 95% rain-fed. As such, it is highly vulnerable to fluctuations in rainfall patterns [6].

The Comoé river basin, which is the focus of this study, is a transboundary catchment shared by Côte d’Ivoire, Burkina Faso, Ghana, and Mali. Climatic variability in this part of Africa is marked by a 20% to 30% decline in rainfall and an even greater reduction in river discharge [7]. This watershed exhibits a pronounced asymmetry between the southern coastal zone, which is directly influenced by moist Atlantic monsoon flows and receives higher rainfall (ranging from 1,200 to 2,000 mm/year), and the northern part, which is dominated by the dry Harmattan winds [8]. The severe global drought episode of 1982–1983 highlighted the need for further research into the causes and impacts of drought in order to better anticipate future events [9]. In Côte d’Ivoire, these extreme climatic events had serious consequences between 1982 and 1983, including forest and plantation fires, a significant decline in agricultural output, and widespread power outages [10].

Drought is a stochastic natural hazard triggered by an intense and persistent lack of precipitation. Following an initial meteorological event, subsequent impacts occur on agriculture and hydrology. Among natural hazards, droughts exhibit unique characteristics: in addition to delayed effects, they vary across multiple dynamic dimensions such as severity and duration. These factors, combined with their pervasive and subjective impacts, make droughts difficult to characterise. Effective drought management therefore requires a clear characterisation of drought in order to facilitate both retrospective analyses and prospective planning [11]. However, drought research faces major challenges due to the scarcity of direct observations, the limitations of the drought indices employed, and geographical inconsistencies in observed drought trends. The complexity and variability of atmospheric mechanisms such as the El Niño–Southern Oscillation (ENSO) and sea surface temperature (SST) anomalies contribute significantly to the difficulty of estimating droughts [1,4]. An efficient drought characterisation requires the identification of contiguous regions with similar drought histories. The first step in risk management is thus to establish management zones by subdividing areas into subregions based on drought criteria. Building on this principle, the regional approach has been incorporated into modelling in recent studies [12,13]. The advantage of this approach lies in its ability to characterise drought evolution in areas with sparse rain gauge networks. Furthermore, the ultimate goal of drought studies remains the prediction of future episodes in order to mitigate their impacts. To this end, more robust stochastic models have been employed by various authors [14,15]. Among the methods used are Markov chains and frequency analyses employing three, four, or even five-parameter probability distributions.

Most studies on the Comoé basin address drought from the perspective of climate variability using the Nicholson index at the annual scale, combined with trend and homogeneity analyses [16,17]. However, meteorological droughts often manifest at finer temporal scales (monthly or seasonal). Focusing solely on the annual scale tends to obscure the actual duration of drought episodes within the year. Only the work of Ouattara et al. [18] offers a genuine characterisation of drought at the monthly scale, including a spatial analysis of the phenomenon. Nevertheless, drought duration is not modelled in terms of return periods. The study by Meledje et al. [19] on the neighbouring Bia basin provides a useful approach to describing the dynamics of drought state transitions using a Markov chain framework. Given the frequency of droughts and the magnitude of their impacts, public authorities should place greater emphasis on research related to this issue in order to build a robust information database. Identifying zones with similar drought manifestation patterns could also help policymakers design efficient response plans. The attempt to address this issue in a drought-vulnerable area specifically the Comoé River watershed, according to global climate models has motivated this study. The objective is to spatially predetermine the meteorological drought “duration” parameter across the Comoé watershed using Regional Frequency Analysis (RFA). The aim of this frequency analysis is to link the duration of drought episodes to their occurrence frequency through probability distributions.

2. Materials and Methods

2.1. Study Area

The Comoé River Basin is a transboundary watershed shared by Côte d’Ivoire, Burkina Faso, Ghana, and Mali (Figure 1). It is located between latitudes 5°02’N and 11°04’N, and longitudes 2°07’W and 5°81’W, within the WGS 84 geographic coordinate system (Zone 30N). The study area lies in the eastern part of Côte d’Ivoire. It is bordered to the north by the Niger River Basin, to the north-east by the Volta Basin, to the south-east by the Bia Basin, to the west by the Bandama Basin and the Agnéby Coastal Basin, and to the south by the Atlantic Ocean in the Gulf of Guinea (Figure 1). Administratively, the basin spans five administrative regions in Côte d’Ivoire, one region in Burkina Faso, and one region in Mali. In Côte d’Ivoire, these are the regions of Gontougo, Bounkani, Poro, Hambol, and Iffou. In Burkina Faso, it extends into the Banfora region, and in Mali, into the Sikasso region [20].

A large part of the Comoé River watershed lies within the eastern region of Côte d’Ivoire. The regional climatic phenomena prevailing in this area help to better characterise the basin’s climate. The eastern part of Côte d’Ivoire encompasses two main climatic zones that correspond to the two dominant landscapes: pre-forest savannah in the northeast and dense, humid forest in the southeast.The climate of the southeast is marked by two rainy seasons: the longer and more intense rainy season peaks in June, while the shorter rainy season occurs around October. These two seasons are separated by a short dry spell during August and September. The major dry season, lasting on average three months, extends from December to February. In contrast, the northeast experiences only one rainy season, with a peak in September. The single dry season lasts between six and seven months and becomes increasingly pronounced moving northwards [8].

Figure 1. Location of the Comoé river basin and the rainfall stations used in the study

2.2. Data and Sources

The data used in this study primarily consist of monthly rainfall totals derived from daily observations. The rain gauge stations are irregularly distributed across the watershed, with most of them located outside the Comoé River Basin. However, only those stations influencing the basin’s rainfall dynamics identified through Thiessen polygon analysis were selected and are shown in Figure 1. The dataset contains some missing values, which were imputed using a regression-based method implemented via the VIM package from the R CRAN (Comprehensive R Archive Network) repository [21]. To ensure a sufficient number of stations for spatial interpolation, the selected data cover the period from 1960 to 2000. Each rainfall time series included in the study spans at least 40 years. Rainfall data for Côte d’Ivoire were obtained from SODEXAM (Société d’Exploitation et de Développement Aéroportuaire, Aéronautique et Météorologique). In addition, several rain gauge stations located in Burkina Faso were incorporated into the analysis, with data sourced from the SIEREM database (Système d'Informations Environnementales sur les Ressources en Eau et leur Modélisation) [22].

2.3. Methods

2.3.1. Drought Index, Threshold and Duration

The Standardised Precipitation Index (SPI) was developed by McKee et al. [23] as a tool to define and monitor drought events. The SPI compares precipitation over a given timescale (typically from 1 to 24 months) to the long-term average precipitation observed at the same location. The calculation steps and associated formulas are detailed in Edwards and McKee [24]. Based on the monthly rainfall series collected over the Comoé River Basin and using the DrinC 1.5 software [25], SPI values were calculated at a 3-month timescale (referred to as SPI-3). The 3-month SPI was chosen in this study because it better reflects short- to medium-term wet or dry conditions. It is also widely used for analysing meteorological droughts [26].

The duration of a drought event is defined as the period between its onset and termination, characterised by a sustained rainfall deficit. Technically, a month is considered dry when the SPI value falls below the threshold of -0.84 (Figure 2). The duration of a dry episode, denoted as L, corresponds to the total number of consecutive months during which SPI remains below -0.84. A drought event ends when the SPI rises above this threshold [27].

Figure 2. Schematic representation of a drought episode with a threshold value of –0.84

2.3.2. Regional Frequency Analysis

- Rainfall gauges SPI-index patterning with Self- Organising Maps

The variables selected for grouping drought manifestation patterns include both geographical and statistical characteristics. These consist of the spatial coordinates X (longitude) and Y (latitude) of the various rainfall stations, as well as descriptive statistical parameters derived from the SPI-3 time series. Specifically, kurtosis, skewness, the first quartile, the median (second quartile), the third quartile, and the interquartile range are used. The classification model is based on Kohonen’s Self-Organising Maps (SOM) [28]. This type of neural network, inspired by sensory activation patterns in the human cerebral cortex, is trained in an unsupervised manner using a simple heuristic approach and is capable of identifying hidden nonlinear structures in high-dimensional data. Fundamentally, a SOM consists of a discrete two-dimensional grid of neurons (Figure 3). The hidden layer, often referred to as the Kohonen layer, is fed by an input layer [29].

Figure 3. Overview of the Self-Organising Map (SOM) neural network

Each neuron in the hidden layer has several “neighbouring” neurons. We define the “distance” between each neuron and its neighbours as the Euclidean distance between the weights connecting the input layer to the hidden layer of each neuron. The distance defined in this way becomes important later, as it encodes low-dimensional information about the original dataset. Each hidden-layer neuron therefore has a weight vector with n components (i.e., the input is densely connected to the hidden layer). How, then, is this network trained? The training procedure and a complete description of the algorithm are provided in Kohonen [28]. The Self-Organising Maps (SOM) functions are implemented in Matlab, enabled the entire classification procedure to be carried out. The validity of the classification (or regionalisation) was further tested using the discordance measure (Di) and the heterogeneity index (H), as recommended by Hosking and Wallis [12].

Discordance measure (Di)

The purpose of this measure is to identify sites (here, the sites are represented by rainfall stations) that are markedly discordant with the group (or class) as a whole in terms of drought manifestation. It is based on the L-moments of the data from the different sites. Consider a group of sites, and suppose that there are N sites in the group. Let us denote by a vector Ui (Eq. 1):

[picture] | (1)

This vector contains the values, and of site, and the superscript denotes the transpose of a matrix or vector. The group mean value is denoted by ū and is expressed as follows (Eq. 2):

[picture] | (2)

The sum-of-squares-and-cross-products matrix is defined as follows (Eq. 3):

[picture] | (3)

Thus, the discordance measure Di for site i is defined as follows (Eq. 4):

[picture] | (4)

A site i is declared discordant if its Di value is large. The notion of “large” depends on the number of sites in the group. A site is considered discordant if its Di value exceeds the critical value given in Table 1.

Table 1. Critical Values for the Discordance Statistic [12]

Heterogeneity measure (H)

The heterogeneity of a proposed region is based on the magnitude of the between-site variability in L-moment ratios compared to the level of variability expected in a homogeneous region [12]. The heterogeneity measures H₁, H₂, and H₃ correspond, respectively, to L-Cv, L-skewness, and L-kurtosis. In practice, H1 is the most informative measure, as it reflects whether the observed variability in L-Cv is significant. Conversely, high natural variability in L-skewness and L-kurtosis, measured by H₂ and H₃ respectively, has relatively low discriminating power. Consequently, the heterogeneity measure H1 for the variability in L-Cv effectively becomes the primary metric for assessing the relative heterogeneity of a given region (or group). The heterogeneity measure H₁ is calculated as follows:

• a regional weighted mean of L-Cv (L-CVR) is computed from the sample L-Cv values of the sites within a region, with weights proportional to the length of each series ;

• The standard deviation of the weighted mean for the L-Cv values across all sites is then calculated and denoted (V) ;

• a four-parameter Kappa probability distribution is fitted using the regional weighted means of the L-moment diagrams from the sites. Subsequently, 500 computer simulations are conducted using the four-parameter Kappa distribution, where each simulation has the same number of sites and series lengths as the proposed region. The mean (µv) and standard deviation (σv) are calculated from the 500 simulated samples ;

• finally, H₁ is computed as follows (Eq. 5):

[picture] | (5)

A value of H₁ equal to zero indicates that the between-site variability of the regional L-Cv values is identical to that expected for a homogeneous region, given the observed L-moment ratios fitted by a four-parameter Kappa distribution. Positive H₁ values indicate that the between-site variability of L-Cv at a site is greater than expected for a homogeneous region, with larger values suggesting potential or probable heterogeneity.

Conversely, negative H₁ values indicate that the between-site variability of L-Cv at a site is lower than expected, and the proposed region would be considered

homogeneous. Experience to date indicates that an H₁ value of 2 is a reasonable threshold for distinguishing regions that are likely heterogeneous from those that are likely homogeneous [12]. Table 2 provides guidelines for the acceptance of homogeneous regions, where external sources contribute to the variability in L-Cv beyond the statistical considerations of sampling variability [30].

Table 2. Guidelines for Accepting or Rejecting the Homogeneity of a Given Region [30]

2.3.3. Selection of a Regional Probability Distribution Using L-Moments

Techniques developed from order statistics, aimed at constructing estimators that better capture the tails of distributions, led to the concept of L-moments. These combinations of order statistics account for distribution tails more effectively than conventional moments. They allow fitting to generalised families of distributions (Generalised Logistic, Generalised Extreme Value, Generalised Pareto, Pearson Type III), Generalised Log-Normal with three parameters, or even four- and five-parameter distributions (Kappa and Wakeby), which are more robust and can have heavier tails than the normal distribution [31]. To select the best-fitting distribution for drought duration series, the L-moment ratio diagram was used, together with the ZDIST statistic, as recommended by Hosking and Wallis [12].

L-moment ratio diagram

L-moment ratio diagrams are increasingly used in the literature to select a probability distribution function for regional frequency analysis [12,28]. In particular, the diagram plots the L-kurtosis against the L-skewness [33]. Two graphical methods are often employed for distribution selection: the sample mean and the best-fit line through the L-moment ratios of the sample. Homogeneous and heterogeneous regional samples have been simulated to highlight the usefulness of these two methods. For homogeneous regional data, the best distribution choice is obtained from the sample mean. For highly heterogeneous regional data, exhibiting a wide range of shape parameters, the best-fit line proves more useful for selecting the distribution. These results emphasise the importance of using heterogeneity tests, such as the ZDIST statistic, alongside L-moment ratio diagrams [34].

Calculation of the ZDIST statistic

According to Ben-Zvi and Azmon [35], the L-moment diagram allows the identification of acceptable distributions, but not necessarily the most appropriate one. Statistical tests should be used in conjunction with the L-moment ratio diagram [36]. To this end, Hosking and Wallis [12] developed the Z statistic test, which aims to assess whether a given theoretical frequency model can adequately represent the statistical behaviour of the parameter under study (here, drought duration) across the different sites (stations) of a homogeneous region (group).

The ZDIST test measures the deviation of the L-kurtosis ratios between the theoretical value of the selected distribution and the mean value obtained from Monte Carlo simulations of the original region. The ZDIST statistic was calculated using an R programme provided by Hosking [37] in the lmomRFA package (Version 3.0-1). The goodness-of-fit was evaluated for five candidate distributions: Generalised Logistic (GLO), Generalised Extreme Value (GEV), Generalised Lognormal (LN3), Pearson Type III (PE3), and Generalised Pareto (GPA). Provided that the region is acceptably homogeneous, the fit can be considered acceptable at the 10% significance level if the absolute value of ZDIST is less than 1.645.

2.3.4. Regional Quantile Estimation

Regional estimation transfers information from gauged sites within a region to an ungauged or partially gauged site for which insufficient data are available. According to Ouarda et al. [38], the most widely used methods that provide adequate estimates are the Index Flood method [39] and the regression method. In this study, the flood index method was employed.

Originally proposed by Dalrymple [39], the basic assumption of the method is that data at different stations are independent and follow the same statistical distribution, up to a scale factor. The method comprises several steps, as described in [12].

2.3.5. Reliability of the Regional Quantile Estimation Method

The reliability of the regional quantile estimation method can be assessed by calculating the root mean square error (RMSE) associated with the regional estimates. This measure quantifies the overall deviation (or error) of the estimated quantiles relative to the true (empirical) quantiles. For each station within the homogeneous region, these values are computed for return periods of 2, 5, 10, 20, and 50 years [36]. Thus, it serves as an important criterion for evaluating the superiority of one model over another. After several iterations (Monte Carlo simulation), the regional relative mean RMSE of the estimated quantile is given by the following equation (Eq. 6):

Here et represent the return-period T quantiles estimated for site i from the regional and local distribution parameters, respectively, and N is the number of stations.

2.3.6. Spatial Interpolation Technique for Estimated Quantiles

The estimated drought duration quantiles were spatially interpolated using the Inverse Distance Weighting (IDW) method. This is a simple and effective interpolation technique based on the assumption that the values of a variable at a location to be estimated are similar to those at neighbouring observation points. The method assumes that each station exerts a local influence, which decreases with distance through the use of a power parameter [40]. In this study, a power value of p = 0.3 was applied.

3. Results and Discussion

3.1. Formation of Rainfall Station Groups Based on Drought Manifestation

Table 3. Matrix Selection Based on Quantisation (QE) and Topology Errors (TE)

The Kohonen self-organising map produced the results shown in Table 3 and Figure 4. The [3×8] matrix was found to be the most appropriate configuration, as it minimised topology and quantisation errors (Table 3).

Figure 4a shows a Kohonen self-organising map (SOM) illustrating the clustering of rainfall stations based on drought manifestation. The stations distributed across the different Kohonen cells were subsequently aggregated using an agglomerative hierarchical clustering (Figure 4b), which groups the stations on the basis of similarities in spatial coordinates and SPI-3 time series statistics. With a cut-off level of 0.7 on the Euclidean distance tree (Figure 4d), four distinct groups were obtained. However, when applying an euclidean distance of 1, three groups were identified (Figure 4c). The transition from three to four groups results from the subdivision of Group 2 in the classification of Figure 4c. For the remainder of the study, the classification yielding four groups or regions, denoted R1 (Orodara, Ouangolo, Abengourou, Ouéllé, Dabakala, Bongouanou, Daoukro), R2 (Abidjan, Adiaké, Alépé, Bassam, Adzopé), R3 (Bouna, Bondoukou, Agnibilekrou, M’Bahiakro, Aboisso) and R4 (Tafiré, Tengréla, Ferkessédougou, Bobo-Dioulasso, Niangoloko, Kampti), was retained (Figure 4d).

The distribution of the different regions identified within the Comoé River Basin roughly corresponds to climatic zoning. Group R2 is characteristic of the transitional equatorial climate, as all stations within this group are located in this zone. The stations in Region R4 (Tafiré, Tengréla, Ferkessédougou, Bobo-Dioulasso, Niangoloko, and Kampti) are typical of drought manifestation in the tropical transitional climate. Groups R1 and R3 appear to be heterogeneous in terms of climate typology: the stations of M’Bahiakro, Aboisso, and Agnibilékrou, which belong to R3, and those of Ouéllé, Daoukro, Bongouanou, and Abengourou, which belong to R1, are situated in the transitional equatorial climate zone, as is R2. The other two stations of R3, namely Bouna and Bondoukou, and the station of Dabakala in R1, are located in the attenuated equatorial climate. Furthermore, Region R1 appears more heterogeneous than R3, with stations such as Orodara and Ouangolodougou situated in the tropical transitional zone, while others are in the transitional equatorial climate.

3.2. Assessment of Cluster Validity Based on Homogeneity and Discordancy Tests

In all groups, the skewness coefficient (L-Skewness = t₃) is less than 0.23. The Hosking and Wallis type 1 (H₁) test is therefore recommended for assessing the heterogeneity of the different groups. The values of the H₁ heterogeneity index for the groups resulting from the classification indicate that they are all homogeneous (Table 4). Groups R1, R2, and R3 are homogeneous with H₁ values of 0.5, –0.27, and –0.56, respectively. They are more homogeneous compared with region R4, which is acceptably homogeneous with H₁ = 1.71. The values of the regional L-Cv (denoted by t) for all groups range between 0.33 and 0.36 in absolute terms, indicating a large regional variability in the distribution of SPI-3 values. These high L-Cv values are also associated with strong observed skewness, reflecting a marked asymmetry or dispersion of the different distributions at the regional scale. These conclusions are not sufficient to determine the homogeneity of a region within the framework of regional frequency analysis. The discordancy statistic (Di) must also be satisfied in order to assess whether a site (station) is discordant or inappropriate within a given region (group). The computed Di values are presented in Table 5. The sites belonging to the different regions are therefore correctly assigned. No manual reorganisation or adjustment of the clusters was carried out, as all groups resulting from the SOM procedure are homogeneous according to the recommended H₁ and Di indices. In the subsequent analyses, these groups will be considered for further investigation of the drought duration parameter.

Figure 4. Group formation: (a) synoptic view of the different Kohonen cells, (b) various cutting levels of the dendrogram, (c) and (d) clusters obtained with Euclidean distances of 1 and 0.7, respectively

Table 4. Regional L-Moments for the Drought Duration Groups

Table 5. L-Moment Statistics and Discordancy Index (Di) for the Sites (Stations) in the Previously Delineated Regions (Groups)

3.3. Regional Probability Distribution of Drought Duration

3.3.1. ZDIST Statistic

The interpretation of the ZDIST statistic results indicates that only the Pearson Type III and Generalised Pareto probability distributions provide the best fit to drought duration in region R2, as the absolute value of the ZDIST statistic is less than 1.645 (Table 6). For the other regions, none of the tested distributions appear suitable for the estimation of drought duration across the basin. This statistic should be combined with the L-moment ratio diagram to ensure an optimal distributional choice.

Table 6. L-Kurtosis Coefficient and ZDIST Statistic of the Drought Duration Series

3.3.2. L-Moment Ratio Diagram

L-moment diagrams (more precisely, diagrams of the L-skewness (t₃) and L-kurtosis (t₄) are similar to classical moment diagrams but have the advantage of being less sensitive to sampling fluctuations. They allow the assessment of the degree of homogeneity of a region and the selection of a theoretical distribution to construct a growth curve. Assuming that the group formed by the rainfall stations in region R1 is homogeneous, all stations would be represented by the same point in these diagrams (generally by the mean). This regional value is denoted by the symbol ‘+’. Based on this diagram, none of the distributions is formally appropriate for predicting the future evolution of drought duration at the regional scale (Figure 5). However, region R2 appears close to the Pearson Type III (PE3) model. Additionally, at the individual site level, some stations resemble the Generalised Logistic distribution in terms of drought duration. In such cases, it is recommended to use the Wakeby or 5-parameter Kappa distributions, which are more robust. Table 7 summarises the probability distributions selected for the subsequent analyses.

Figure 5. L-moment ratio diagrams of drought duration for the groups or regions

Table 7. Parameters of the Probability Distributions Selected for the Estimation of Drought Duration

3.4. Estimation of Regional Drought Quantiles Associated with Return Periods

The quantile here represents occurrence of a drought of a given duration (in months) over a return period of T years. Longer drought durations for a given return period indicate a higher risk of drought onset. Figure 6 shows the probable drought durations for several return periods. For a return period T of 5 years, the probable drought duration would be on average 3 months for groups R1, R3, and R4, whereas it would be 4 months for group R2, characteristic of the transitional equatorial climate. For T = 10 years, this duration increases to 4 months for groups R1, R3, and R4, while it would be approximately 5 months for group R2. For long-term estimation (T = 50 years), groups R2 and R1 would have a probable duration of 8 months, and groups R3 and R4 a duration of 6 months. However, this variability in estimated duration remains relatively small. For example, for return periods between 1 and 4 years, the estimated drought durations, ranging from 2 to 3 months, are relatively homogeneous across all groups. For a 50-year return period, the drought duration would be estimated at 6–7 months for groups R1 and R2, whereas it would be 5–6 months on average for groups R3 and R4. The calculation of the relative mean square error (RMSE) indicates that the selected distributions provide the best fit to the different components of drought at the regional scale. Indeed, the computed RMSE values are below 15% for the regional growth curves.

Figure 6. Regional quantiles of drought duration estimated for different return periods

Only the extreme right-hand tails exhibit errors approaching 50%, but these remain within acceptable limits. Regional estimation appears more reliable than site-specific (station-based) estimation, given the smaller RMSE values in the regional model compared with the site model. Consequently, site-specific quantiles will be estimated by normalising them using the Index Flood method. Following the estimation of station-based quantiles on the basis of the Index Flood, Section 3.5 addresses the interpolation of these estimated durations in order to provide insight into the manifestation of drought across the Comoé River basin.

3.5. Mapping of Station-Based Estimated Durations of Meteorological Drought

The catchment is predominantly characterised by short and relatively homogeneous drought durations over most of its area (Figure 7). In contrast, in the far north of the basin (the portion located in Burkina Faso), longer drought durations are observed compared with the rest of the basin. Drought duration increases with return period, ranging from a maximum of 4 months for T = 2 years to 20 months for a return period of 100 years. By way of illustration, the most probable maximum drought duration is around 2, 4, and 5 months for return periods of 2, 5, and 10 years, respectively (short-term prediction), across the majority of the basin. For these same return periods, drought durations of 4, 6, and 8 months are expected in the extreme north of the basin. For long-term prediction (T = 20 and 50 years), the maximum expected drought persistence ranges between 6, 8, and 10 months over most of the basin. In the northern part of the basin, rainfall deficiency may extend over a full hydrological year or even two, with maximum drought durations of 11 and 16 months corresponding to return periods of 20 and 50 years, respectively.

Figure 7. Spatial distribution of meteorological drought duration for return periods of 2, 5, 10, 20, and 50 years

3.6. Discussion

The classification produced four groups on the basis of the model’s input parameters. These regions appear well delineated due to the variation in kurtosis from one region to another, confirming the discriminating power of the classification. Furthermore, the Hosking and Wallis [12] heterogeneity tests applied to the pre-defined regions corroborate their homogeneity. Ward’s criterion remains a robust classification method, as highlighted by [41]. The correspondence between the regions and climatic zones is discussed in this section. Regions R2 and R4 are characteristic of the transitional equatorial climate and transitional tropical climate, respectively. By contrast, the remaining groups, R1 and R3, reflect a form of spatial and/or climatic heterogeneity, since stations in these regions belong to multiple climatic zones. Such heterogeneity in drought manifestation has also been highlighted in previous studies and is evident both at national and basin scale [38,39]. According to Vicente-Serrano [13], the spatial distribution of drought is generally complex, with marked inhomogeneity even between neighbouring regions. This variability may be attributed to the complexity of atmospheric circulation and to the fact that drought cannot be linked to a single type of atmospheric condition. Moreover, pronounced inhomogeneity in drought distribution is most frequently observed in transitional climatic regions [44]. The Comoé River basin, which is the focus of this study, spans three distinct climatic zones from north to south. Additionally, this transboundary basin lies within the West African monsoon domain, one of the most complex climatic systems.

In the process of stochastic modelling, uncertainties must be minimised as much as possible, and robust methods should be employed at each stage. The use of L-moments in this study is consistent with previous research aimed at improving the robustness of estimation methods. According to Lana et al. [45], the Euclidean distance defined by L-skewness and L-kurtosis is an effective indicator for achieving a better fit between empirical distributions and statistical models, especially when sample sizes are small and the 95% Kolmogorov–Smirnov confidence interval is not restrictive. Furthermore, Hosking and Wallis [12] demonstrated that the statistical differences between return periods derived from conventional moments and those derived from L-moments are non-negligible. The adoption of more robust estimation methods (L-moments), as in the present study, is therefore recommended. The combination of the L-moment ratio diagram and the ZDIST statistic proved decisive in selecting the final regional model. Based on ZDIST, the GNO, PE3, and GPA distributions provided a satisfactory fit at site (station) scale. However, at the regional scale, only the Wakeby distribution was found to be the most appropriate. According to Hosking and Wallis [12], the five-parameter Wakeby distribution has the advantage of encompassing the forms of several asymmetric distributions commonly used (GEV, GNO, PE3, etc.). The technical parameterisation of the Wakeby distribution reveals a generalised form of the Generalised Pareto probability law. In drought analyses by [13] and [46], the Wakeby distribution has been widely used to simulate both drought duration and severity using the SPI and SPEI indices. Hosking and Wallis [12] suggest that two- or three-parameter distributions are adequate for site-based frequency analysis, whereas regional analyses are better suited to three- to five-parameter distributions for improved estimation.

High probabilities of long drought durations are observed in the transitional tropical climate (Region R4) rather than in areas close to the equator. The risk of drought occurrence is relatively homogeneous across the basin, except in the extreme north. In this part of the basin, drought severity is exacerbated by long durations for return periods of 20 and 50 years. This northern zone is the most arid, being close to the Sudano-Sahelian climate. According to Lebel and Vischel [47], tropical regimes exhibit two features that distinguish them from mid-latitude regimes: (i) their seasonal cycle is well marked and more stable from year to year; and (ii) rainfall is below potential evapotranspiration (PET) for most of the year, particularly in semi-arid subtropical areas. These two features help to explain why, apart from areas close to the equator where the ITCZ influence is almost permanent, tropical regions are especially sensitive to drought, even though their annual rainfall totals are often comparable to, or greater than, those of temperate regions. In these tropical zones, water balances are controlled by rainfall amounts, whereas in temperate regions they are controlled primarily by temperature. With respect to future drought occurrence, the findings of this study are consistent with projections based on global climate models, which suggest that drought frequency and intensity will likely increase in the Mediterranean basin and West Africa, but decrease in central North America and north-western Australia since 1950 [1]. Africa is expected to be particularly vulnerable to extreme droughts affecting multiple regions simultaneously [42]. More localised studies further indicate that semi-arid and sub-humid regions in Africa are especially prone to drought [48]. The uncertainties estimated using RMSE remain within acceptable limits, lending credibility to the modelling results. However, it is important to note that considerable uncertainties and inconsistencies persist in drought projections for Africa, owing to model limitations and the complexity of multiple mechanisms influencing rainfall trends [42]. These uncertainties are further compounded by the scarcity of direct observations, the dependence of inferred trends on the choice of drought indices, and the geographic inconsistency of observed trends [1].

4. Conclusion

Current and future climatic conditions project the persistence of hydrological extremes (droughts and floods). At the same time, pressures on water resources are intensifying due to the diversification of uses and demographic growth. One of the main challenges for basin managers remains ensuring the availability of resources in both quality and quantity, while also establishing an early warning system to mitigate the adverse impacts of drought. This study was undertaken to construct a spatio-temporal perspective of the occurrence of meteorological drought at the scale of the Comoé River basin. A threshold of –0.84 in the Standardised Precipitation Index (SPI) was adopted as the limit below which drought conditions are declared. Regional analysis subsequently identified four homogeneous regions, one of which is representative of the transitional tropical climate (Region R4) in the north, another representative of the transitional equatorial climate (Region R2) in the south, and two others (Regions R1 and R3) with contrasting climatic characteristics. The reliability of the classification, assessed using the discordance (Di) and heterogeneity (H₁) indices, led to the implementation of regional frequency analysis. Accordingly, the five-parameter Wakeby probability distribution was found to be the most appropriate for stochastically reproducing the duration of meteorological drought across the basin. The estimation results indicate prolonged drought events in the far north of the basin (transitional tropical climate), in contrast with the central and southern zones, which are characterised respectively by attenuated transitional equatorial climate and transitional equatorial climate. This study concludes that the Comoé River basin is vulnerable to extended droughts. This vulnerability is further amplified by the basin’s strong dependence on rain-fed agriculture and by the lack of water resource management and drought disaster management infrastructure. These important preliminary findings provide a foundation for more detailed research aimed at establishing an integrated system and a chain of scientific approaches for drought early warning in the transboundary Comoé basin.

ACKNOWLEDGEMENTS

All this work could not be done without material and financial support. I would like to express my special thanks to the facilitators of the PAES (Programme d’Appui à l’Enseignement Supérieur) in the ECOWAS region for agreeing to finance my research project.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1]  
IPCC, “Summary for Policymakers,” in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex and P.M. Midgley (eds.)], Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press, 2013. [Online]. Available: https://www.ipcc.ch/site/assets/uploads/2018/03/WG1AR5_SummaryVolume_FINAL.pdf.

[2]   | V. Murray and K. L. Ebi, “IPCC Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (SREX),” J. Epidemiol. Community Health, vol. 66, no. 9, pp. 759–760, Sept. 2012.
[3]   | C. Lécuyer, “Evolution de la désertification en Afrique de l’ouest,” Institut de Recherche pour le Développement, Montpellier, Rapport de stage, 2012.
[4]   | S. E. Nicholson, “Climatic and environmental change in Africa during the last two centuries,” Clim. Res., vol. 17, pp. 123–144, 2001.
[5]   | D. A. Wilhite, Ed., Drought Assessment, Management, and Planning: Theory and Case Studies. Boston, MA: Springer US, 1993.

[6]  
AfDB, “Lutte contre la sécheresse en Afrique de l’Ouest et australe : la BAD octroie une enveloppe de 549 millions $,” Banque Africaine de Développement. Accessed: Sept. 16, 2025. [Online]. Available: https:// www.afdb.org/ fr/news-and-events/ afdb-group-announces-us-549-million-drought-response-package-for-eastern-and-southern-africa-15547.

[7]   | Y. L’Hôte, G. Mahé, B. Somé, and J. P. Triboulet, “Analysis of a Sahelian annual rainfall index from 1896 to 2000; the drought continues,” Hydrol. Sci. J., vol. 47, no. 4, pp. 563–572, Aug. 2002.
[8]   | D. Noufé, “Changements hydroclimatiques et transformations de l’agriculture : l’exemple des paysanneries de l’Est de la Côte d’Ivoire,” Carnets Géographes, no. 3, Dec. 2011.
[9]   | D. A. Wilhite and M. H. Glantz, “Understanding: the Drought Phenomenon: The Role of Definitions,” Water Int., vol. 10, no. 3, pp. 111–120, Jan. 1985.
[10]   | FAO, Situation des forêts dans le monde., FAO. Rome, 1999.
[11]   | A. Zargar, R. Sadiq, B. Naser, and F. I. Khan, “A review of drought indices,” Environ. Rev., vol. 19, pp. 333–349, Dec. 2011.
[12]   | J. R. M. Hosking and J. R. Wallis, Regional frequency analysis: an approach based on L-moments. Cambridge: Cambridge University Press, 1997.
[13]   | S. M. Vicente-Serrano, “Spatial and temporal analysis of droughts in the Iberian Peninsula (1910–2000),” Hydrol. Sci. J., vol. 51, no. 1, pp. 83–97, Feb. 2006.
[14]   | I. Abi-Zeid and B. Bobée, “La modélisation stochastique des étiages: une revue bibliographique,” Rev. Sci. Eau, vol. 12, no. 3, pp. 459–484, Apr. 2005.
[15]   | A. A. Paulo and L. S. Pereira, “Stochastic Prediction of Drought Class Transitions,” Water Resour. Manag., vol. 22, no. 9, pp. 1277–1296, Sept. 2008.
[16]   | E. A. Assemian et al., “Étude de l’impact des variabilités climatiques sur les ressources hydriques d’un milieu tropical humide: Cas du département de Bongouanou (Est de la Côte d’ivoire),” Rev. Sci. L’eau, vol. 26, no. 3, pp. 247–261, Oct. 2013.
[17]   | K. E. Kouakou, B. T. A. Goula, and I. Savane, “Impacts de la variabilité climatique sur les ressources en eau de surface en zone tropicale humide: cas du bassin versant transfrontalier de la comoé (côte D’ivoire - burkina faso),” Eur. J. Sci. Res., vol. 16, no. 1, pp. 31–43, 2007.
[18]   | I. Ouattara et al., “Assessment of Meteorological Drought in Climate Variability Context within the Comoe River Transboundary Watershed,” Atmospheric Clim. Sci., vol. 12, no. 01, pp. 1–17, 2022.
[19]   | N. H. Meledje, K. Lazare Kouassi, Y. A. N’Go, and I. Savane, “Characterization of dry period in the transfrontier Bia watershed in Côte d’Ivoire and Ghana: contribution of Markov chains,” Cah. Agric., vol. 24, no. 3, pp. 186–197, May 2015.
[20]   | JICA, “Plan directeur de gestion intégrée des ressources en eau en République de Côte d’Ivoire,” JICA, Côte d’Ivoire, Rapport principal, 2001.
[21]   | A. Kowarik and M. Templ, “IImputation with the R Package VIM,” J. Stat. Softw., vol. 74, no. 7, 2016.
[22]   | J. F. Boyer et al., “SIEREM: An environmental information system for water resources,” in Proceedings of the 5th FRIEND World Conference on Water resource variability: hydrological impacts, Havana, Cuba, 2006.
[23]   | T. B. McKee, N. J. Doesken, and J. Kleist, “The relationship of drought frequency and duration to time scales,” in Proceedings of the 8th conference on applied climatology, Anaheim: American Meteorological Society, 1993, pp. 179–184.
[24]   | D. C. Edwards and T. B. McKee, “Characteristics of 20th century drought in the United States at multiple time scales,” Colorado State University Dept. of Atmospheric Science, Fort Collins Colorado, Climatology Rep, 97-2, 1997.
[25]   | D. Tigkas, H. Vangelis, and G. Tsakiris, “DrinC: a software for drought analysis based on drought indices,” Earth Sci. Inform., vol. 8, no. 3, pp. 697–709, Sept. 2015.
[26]   | M. D. Svoboda and B. A. Fuchs, Handbook of drought indicators and indices. Geneva: World Meteorological Organization, 2016.
[27]   | C. T. Agnew, “Using the SPI to identify drought,” Drought Netw. News, vol. 12, no. 1, pp. 6–12, 2000.
[28]   | T. Kohonen, Self-organizing maps, 3rd ed. in Springer series in information sciences, no. 30. Berlin ; New York: Springer, 2001.
[29]   | S. S. Haykin, Neural networks and learning machines, 3rd ed. New York: Prentice Hall, 2009.
[30]   | M. Ribatet, E. Sauquet, J.-M. Grésillon, and T. B. M. J. Ouarda, “A regional Bayesian POT model for flood frequency analysis,” Stoch. Environ. Res. Risk Assess., vol. 21, no. 4, pp. 327–339, Mar. 2007.
[31]   | T. Öztekin, “Wakeby distribution for representing annual extreme and partial duration rainfall series,” Meteorol. Appl., vol. 14, no. 4, pp. 381–387, Dec. 2007.
[32]   | R. M. Vogel and N. M. Fennessey, “L-moment diagrams should replace product moment diagrams,” Water Resour. Res., vol. 29, no. 6, pp. 1745–1752, June 1993.
[33]   | Y. D. Chen, G. Huang, Q. Shao, and C.-Y. Xu, “Regional analysis of low flow using L-moments for Dongjiang basin, South China,” Hydrol. Sci. J., vol. 51, no. 6, pp. 1051–1064, Dec. 2006.
[34]   | M. C. Peel, Q. J. Wang, R. M. Vogel, and T. A. McMAHON, “The utility of L-moment ratio diagrams for selecting a regional probability distribution,” Hydrol. Sci. J., vol. 46, no. 1, pp. 147–155, Feb. 2001.
[35]   | A. Ben-Zvi and B. Azmon, “Joint use of L-moment diagram and goodness-of-fit test: a case study of diverse series,” J. Hydrol., vol. 198, no. 1–4, pp. 245–259, Nov. 1997.
[36]   | K. Benhattab, C. Bouvier, and M. Meddi, “Analyse fréquentielle régionale des précipitations journalières maximales annuelles dans le bassin hydrographique - Chéliff, Algérie,” Rev. Sci. L’eau, vol. 27, no. 3, pp. 189–203, Dec. 2014.
[37]   | J. R. M. Hosking, lmomRFA: Regional Frequency Analysis using L-Moments. (2015). CRAN. Accessed: Sept. 17, 2025. [Online]. Available: https://CRAN.R-project.org/package=lmomRFA.
[38]   | T. B. M. J. Ouarda, M. Lang, B. Bobée, J. Bernier, and P. Bois, “Synthèse de modèles régionaux d’estimation de crue utilisée en France et au Québec,” Rev. Sci. Eau, vol. 12, no. 1, pp. 155–182, Apr. 2005.
[39]   | T. Dalrymple, “Flood-frequency analyses,” USGS, Washington DC, Manual of Hydrology Part 3, 1960.
[40]   | S. M. Pingale, D. Khare, M. K. Jat, and J. Adamowski, “Spatial and temporal trends of mean and extreme rainfall and temperature for the 33 urban centers of the arid and semi-arid state of Rajasthan, India,” Atmospheric Res., vol. 138, pp. 73–90, Mar. 2014.
[41]   | Y. Unal, T. Kindap, and M. Karaca, “Redefining the climate zones of Turkey using cluster analysis,” Int. J. Climatol., vol. 23, no. 9, pp. 1045–1055, July 2003.
[42]   | I. Masih, S. Maskey, F. E. F. Mussá, and P. Trambauer, “A review of droughts on the African continent: a geospatial and long-term perspective,” Hydrol. Earth Syst. Sci., vol. 18, no. 9, pp. 3635–3649, Sept. 2014.
[43]   | F. E. F. Mussá, Y. Zhou, S. Maskey, I. Masih, and S. Uhlenbrook, “Groundwater as an emergency source for drought mitigation in the Crocodile River catchment, South Africa,” Hydrol. Earth Syst. Sci., vol. 19, no. 2, pp. 1093–1106, Feb. 2015.
[44]   | C. Rodriguez-Puebla, A. H. Encinas, S. Nieto, and J. Garmendia, “Spatial and temporal patterns of annual precipitation variability over the Iberian Peninsula,” Int. J. Climatol., vol. 18, no. 3, pp. 299–316, Mar. 1998.
[45]   | X. Lana, M. D. Martínez, A. Burgueño, C. Serra, J. Martín‐Vide, and L. Gómez, “Distributions of long dry spells in the iberian peninsula, years 1951–1990,” Int. J. Climatol., vol. 26, no. 14, pp. 1999–2021, Nov. 2006.
[46]   | M. B. Masud, M. N. Khaliq, and H. S. Wheater, “Analysis of meteorological droughts for the Saskatchewan River Basin using univariate and bivariate approaches,” J. Hydrol., vol. 522, pp. 452–466, Mar. 2015.
[47]   | T. Lebel and T. Vischel, “Climat et cycle de l’eau en zone tropicale : un problème d’échelle,” Comptes Rendus Géoscience, vol. 337, no. 1–2, pp. 29–38, Dec. 2004.
[48]   | World Bank and GFDRR, “Disaster Risk Profile Mozambique,” World Bank: Washington, DC, USA. Accessed: Sept. 05, 2024. [Online]. Available: https://www.worldbank.org.