Floods are natural hazards aggravated by both climatic factors (i.e. climate variability and climate change) and non-climatic factors (e.g. changes in land cover, use, vegetation, etc.) , and result in the destruction and disruption of socio-economic activities, damage to property and infrastructure, loss of lives, and financial loss . In Nigeria—the case-study for this research, frequent and unprecedented levels of flooding and impacts has increased the concern from the public, government, and other stakeholders about the probability of flood recurrence, thus reinforcing the need to establish appropriate mitigation measures to minimize flood impacts .
Knowledge of flood frequency estimates is crucial to ensure socio-economic activities and infrastructural development are planned appropriately to improve resilience . Accurate estimates of flood intensities and frequencies are also important for the design of critical infrastructure required to flood risk reduction (dykes, levees, dams, etc.), construction of hydraulic structures (bridges, culverts, drainages), the development of floodplain and urban land-use regulations, emergency management, and disaster risk insurance . Under-estimating design flood could lead to increased flood risk with potentially catastrophic consequences, while over-estimation, on the other hand, could cause resource wastage and aggravate upstream and downstream flooding .
To accurately estimate expected flood magnitudes and return periods, networks of gauging stations are typically established to collect hydrological data over a long period. However, in many developing regions, establishing the optimal number of hydrological stations is usually hampered by challenges such as the high cost associated with gauging equipment. Therefore, several locations are usually left ungauged or contain short-length of record if the gauging stations are newly established, discontinued, or damaged. In many low/middle-income countries, many catchments are sparsely gauged due to factors that include lack of commitment by station operators; deteriorating conditions of observation equipment; insecurity/theft, and inaccessibility to remote locations  . The absence of high-quality and sufficient data results in poor flood predictions in these areas, and consequently, flawed flood risk management interventions . Therefore, it is essential to explore techniques capable of extracting maximum value from available data to develop robust flood frequency estimates in data-poor regions with reduced uncertainty.
Generally, the choice of a flood frequency estimation approach depends on the extent of availability of historical flood records at/or around the specific site of interest. When sufficient historical flood data are available, direct flood frequency analysis is performed by fitting a pre-defined probability distribution to the annual maximum flood or partial flood time series . Where data is insufficient, indirect flood estimation procedures are used, such as the adoption of hydro-meteorological data from other locations similar in characteristics to the site of interest   or the incorporation of data from other sources including remote sensing  . The present study adopts the former approach.
A major factor that affects future flood regimes and must be considered when estimating flood magnitudes is the changing climate, whether characterized by shorter-term cyclical climate variability or longer-term trends related to climate change. Climate variability refers to the Spatio-temporal deviations/anomalies of climatic statistics over a given period (e.g. a month, season or year), while climate change refers to the longer-term non-cyclical trends that are projected to lead to more intense and frequent flooding . Estimating flood intensities and frequencies under changing climatic conditions requires the incorporation of the effects of non-stationarity defined by statistically significant breakpoints and trends within historical time series  . Understanding the effects of climate on extreme hydro-meteorological events in a region requires the integration of measures of climate variability such as hydro-climatic indices  . Several studies have revealed that incorporating climatic variability into flood frequency estimation procedures result in more robust estimates   .
This study aims to tackle the problem of data sparsity and limited resources to estimate flood magnitude and frequency while taking into consideration the effects of climate variability using open-source tools. The study is presented as follows: Section (2): Description of the study area, datasets and sources; Section (3): Methodology-Data preparation, preliminary analysis and L-moment statistics-based regional flood frequency techniques, taking climate variability into account (climate indices); Section (4): Presentation of the results of preliminary analysis, direct and regional flood frequency estimates and incorporating climate variability; and Section (5): Conclusion and implication of the results on flood risk management.
2. Study Area, Datasets and Sources
The Ogun-Osun River Basin (OORB) is in Western Nigeria (6˚30' - 8˚20'N latitude and 3˚23' - 5˚10'E longitude), and encompasses four states including Ogun, Osun, Oyo and Lagos, within an area of 66,264 km2. The Basin is drained by two major tributaries, the Ogun and Osun Rivers, and other minor tributaries including the Yewa, Ibu, Ona, Sasa and Ofiki Rivers, into the Lagos Lagoon (See Figure 1).
Figure 1. Study area: The Ogun-Osun River Basin (OORB), Nigeria.
The climate of the OORB is influenced by tropical continental and maritime air masses  and experiences an annual rainfall of 1400 mm to 1500 mm; mean annual air temperature ranges from 25.7˚C to 30˚C, and relative humidity varies from 37% to 85% . OORB experiences recurring flooding caused by the increased frequency of intense precipitation events, poor urban planning and waste management practices, and failure of upstream hydraulic systems, which have led to significant socio-economic, infrastructural, ecological and environmental impacts . Also, recent evidence from studies in West Africa    and Nigeria  suggests the presence of strong correlations between climatic variability and hydro-meteorological events in the regions.
Hydrological data (discharge, water levels and rating curves) used for this study were provided by the Ogun-Osun River Basin Development Authority (OORBDA), the agency responsible for collecting and managing data within the Basin. Additional data sets for two hydrological stations, Yewa Mata and Ona River/Sala village were extracted from published research by Olukanni and Alatise  and Ewemoje and Ewemooje  respectively, using WebPlotDigitizer tool . The catchment area for each station was delineated from 30-m resolution Shuttle Radar Topography Mission (SRTM) digital elevation model data  using the Arc Hydro tool in ArcMap 10.2. The properties of the gauging stations for OORB are presented in Table 1, while the spatial distribution of gauging stations is presented in Figure 1, showing the spread and sparsity of the hydrological monitoring network.
Climate indices data used in this study originates from the data repository of the United States National Oceanographic and Atmospheric Administration (NOAA)  and are embedded within the open-source International Centre for Integrated Water Resources Management Regional Analysis of Frequency
Table 1. Properties of OORB gauging stations.
Tool (ICI-RAFT). This data includes multi-decadal hydro-climate indices such as the Pacific Decadal Oscillation, El Nino/Southern Oscillation, Madden Julian Oscillation, North Atlantic Oscillation and others. As a data limitation, it is important to note the variation in the period of data availability for the respective sites restricts the direct comparison of sites for analytical purposes.
3.1. Data Preparation and Preliminary Analysis
Prerequisite data preparation undertaken includes data formatting and conversion, infilling of missing data and statistical testing. First, historical river water level data were converted to discharge using available up-to-date rating curves provided by the OORBDA. Multiple imputation techniques based on coupled Markov Chain Monte Carlo and ordinary least squares regression is applied using Microsoft XLSTAT tool to fill the gaps to approximate missing peak annual discharge (flood) . FFA application typically assumes that available data satisfies conditions of randomness, serial non-correlation, the absence of outliers, and homogeneity, to reduce the inherent data uncertainty .
The randomness of hydrologic data points at each station is estimated using the Mann-Kendall (M-K) test , which assesses the increasing and decreasing trends in the time series . The presence of serial correlation within hydrological records at a particular site results in discrepancies in regional variance and increased data skewness , thus contributing to uncertainty in regional flood frequency estimates . To assess the magnitude of the serial correlation, 1-unit Lag correlation coefficients  are applied to derive values ranging from −1 (perfect non-correlation) to 1 (perfect correlation). The presence of outliers also affects data quality and can be attributed to gauge failure, sampling inconsistencies, typo errors, or gauge disruptions; they are not considered part of the real flood population . Outliers were identified by using the Grubbs and Beck test . Finally, breakpoint analysis  is applied to assess significant homogeneity within the hydrological time series.
3.2. Climate Indices—Assessing Climate Variability Effect on Flood Frequency Estimates
Climate variability depicts the pattern of climate dynamics on both temporal and spatial scales, identified as fluctuation above or below the average climate pattern for a short period, thus influencing the magnitude and frequency of extreme flood events  . While past hydrologic models have assumed stationarity, current climatic conditions suggest that the future is expected to differ despite what is known of the past and present . Ocean-atmosphere processes that influence precipitation, atmospheric pressure and temperature are defined by climatic indices and are useful in tracking long-term decadal hydrological changes  . Some key climate indices that characterize the frequency, intensity and duration of extreme climatic events include the Arctic Oscillation (AO), North Pacific Oscillation (NPO), North Atlantic Oscillation (NAO), Pacific Decadal Oscillation (PDO), Pacific/North American Index (PNA), El Nino/Southern Oscillation (ENSO), and Madden-Julian Oscillation (MJO) .
In this study, the correlation between the annual maximum discharge time-series and climatic indices are evaluated as an indicator of climate variability influence on flood magnitude and frequencies  . ICI-RAFT developed by Giovannettone and Wright,  contains a database of 30 hydro-climate indices (HCI), including those previously mentioned, to facilitate correlation analysis, taking into account a stipulated Lag, as well as lower and upper index limits. To consider climate variability influence on flood frequency analysis, the ICI-RAFT program recomputes flood magnitude for each return period using only non-zero data that satisfies the index limits set for the HCI with the highest correlation coefficient (R2) .
3.3. L-Moment—Index Flood Regional Flood Frequency Analysis (RFFA)
Regional flood frequency analysis is based on the agglomeration of hydrological data from homogeneous regions characterised by similar physiographical parameters (e.g. catchment area, catchment slope, stream length, precipitation, and elevation). Hydrological data available at sites within a defined region are used to estimate the regional flood quantile based on the assumption that they are defined by the same probability distribution and differ only by the index flood . This approach helps reduce inconsistencies associated with data shortage .
The Index flood technique developed by Dalrymple  has been widely applied in determining flood estimates for gauged and ungauged catchments of varying sizes at the global, regional and local scales   . The general assumption for this method is that the probability distribution of the annual maximum floods across sites in the region are similar and differ only by a site-specific scaling factor termed the index flood (mean or median flood)  . The flood quantile ( ) for a T-year return period at a site of interest (i), given a common regional probability distribution factor ( ), can be mathematically expressed as:
The index-flood ( ) for an ungauged site of interest is usually derived from establishing a relationship between available catchment characteristic (such as catchment area, elevation, annual precipitation, etc.) and the index-flood of gauged sites within a homogeneous region . The regional probability distribution is a dimensionless parameter determined using a best-fit statistical approach discussed later in Section 3.3.3.
L-moment based flood frequency analysis was undertaken using ICI-RAFT , and the procedure includes 1) data screening and site clustering to derive discordancy measure (D) based on the Wards hierarchical clustering approach, 2) regional homogeneity testing using the heterogeneity measure (H), and 3) selection of the appropriate distribution using the goodness-of-fit measure (Z) . The L-moments is a widely-preferred method for regional flood frequency analysis due to the robustness of linear moments in comparison to ordinary moments in handling extreme values over a wider range of probability distributions and its reduced susceptibility to bias. The components of L-moment analysis are detailed in Hosking and Wallis  and other studies  and are summarized below.
3.3.1. Data Screening
The discordancy measure based on L-Moments (L-Mean, L-Covariance, L-Kurtosis and L-Skewness) is applied to identify sites whose L-Moment ratio are discordant from that of the whole group, denoted by a critical value of (D ≥ 3). Where, L-mean is similar to conventional mean, defined as a measure of central tendency; L-Covariance is defined as a dimensionless measure of variability, ranging from low (0) to very high (0.4) variability; L-Skewness (t3) is a measure of the degree of symmetry of a sample. L-Skewness value typically lies between −1 and +1, t3 = 0 suggests symmetric distribution. L-Kurtosis (t4), is a measure of peakedness or the flatness of the frequency distribution curve near its centre. Formulas that define L-moment statistics are presented in Hosking and Wallis .
3.3.2. Homogeneity Testing
Heterogeneity measure (H) compares the variation between L-moments for a group of sites and what is expected of a homogeneous region to justify that a group of sites are defined by a similar probability distribution. The region is deemed acceptably homogeneous if H <1, possibly heterogeneous if 1 ≤ H < 2, and heterogeneous if H ≥ 2 .
3.3.3. Probability Distribution Selection
The Z-Statistic is a goodness-of-fit measure that assesses the probability distribution that best fits the weighted-average regional L-moment parameters of each site in a homogeneous region (L-Skewness and L-Kurtosis). A preliminary probability distribution can also be approximated and visualized using an L-moment diagram (L-Kurtosis vs. L-Skewness), with the best distribution defined as the distribution curve closest to the majority of the sample data points .
4. Results and Discussion
4.1. Data Characteristics and Preliminary Analysis
Data preparation results are presented in Table 2. The 1-unit Lag correlation results show that the serial correlation between data sets at each site varied from −0.002 to 0.516 (−1 = perfect inverse correlation; 1 = perfect correlation; and 0 = no correlation), suggesting the absence of a strong relationship among peak annual discharge at each site. Grubbs and Beck test detected no low outlier, but high
Table 2. Results of preliminary analysis.
Outliers and Homogeneity presented as p-values.
outliers observed at Oba/Oyo-Ogbomosho, Ofiki/Ofiki town, Ofiki/Iganna-Ilere road, Ofiki/Igangan, Ogun/Shepeteri, Ogun/Oyo-Iseyin road, and Ogun/Ibaragun gauging stations were consistent with recorded flood events reported in previous studies . Trend and breakpoint analysis (homogeneity test) results revealed significant upward trends at Ijaka-Oke, Oyan/Ilaji-Ile, and Oba/Oyo-Ogbomosho stations, while no significant trends were identified at the remaining sites. These trends were consistent with those of the neighbouring Oueme River Basin in the Benin Republic, which is influenced by similar climatic conditions . The time series plots presented in Figures 2(a)-(d) show the annual maximum discharge of selected stations for further analysis. These selections capture the varying spectrum of trends, displaying spikes and troughs that represent flood variability at Ijaka-Oke and Ofiki-Igangan (Figure 2(a) & Figure 2(b)). Changes in hydrologic regimes defined by the breakpoint analysis were found at Ofiki/Iganna-Ilere road and Oba/Oyo-Ogbomosho stations and are presented in Figure 2(c) & Figure 2(d). Changes in the hydrological regime are evident in the breakpoint analysis plots from 1965 to 1957 and 1979 to 1988, corresponding to documented years of changes in precipitation patterns in Nigeria and West Africa that suggest intense drought to rainfall zone transitions  .
Figure 2. Trend and breakpoint plots for selected gauging stations at (a) Ijaka-Oke, (b) Ofiki-Igangan, (c) Ofiki/Iganna-Ilere road, and (d) Oba/Oyo-Ogbomosho. Where: mu1 (---) and mu2 (---) represent the mean discharge of both breakpoints, representing the average difference in hydrological regimes.
4.2. Identification of Homogeneous Regions and Determination of Discordancy Measure
Regional discordancy (D) and heterogeneity (H) statistics are presented in Table 3, while site-specific results of the same statistics for sites within each sub-region are presented in Table 4. An H-statistic of 8.89 (i.e. H > 1) was realized for the entire catchment area, suggesting heterogeneity . Consequently, the region was divided into three sub-regions and tested for homogeneity (Table 3); the resulting L-moment statistics and discordancy of sites constituting each sub-region are presented in Table 4. The H-Statistics for sub-regions 2 and 3 showed homogeneity (H < 1), while sub-region 1 exhibited heterogeneity (H >> 1). In terms of discordancy, only Idogo was discordant (D = 4.2232) and, therefore,
Table 3. L-moments and Homogeneity statistics per sub-region.
L = Linear, CV = Covariance, Skew = Skewness, Kurt = Kurtosis, Dis = Discordancy, H = Heterogeneity statistic.
Table 4. L-moment and discordancy statistics for gauging stations per sub-region.
was removed from further analysis. The discordancy of other sites in the homogeneous sub-regions was within the prescribed critical limit of D < 3.
The agglomeration of historical data within the homogeneous sub-regions resulted in improved data availability for flood magnitude estimation, i.e. a combined 126 data points (sub-region 2) and 141 data points (sub-regions 3), thus exceeding the 30 years data length required for regional FFA in Nigeria, according to the Technical Guidelines on Soil Erosion, Flood and Coastal Zone Management .
4.3. Regional Distribution and Goodness of Fit Measures
The L-Moment ratio diagram (Figure 3) displays the relationship between
Figure 3. L-Moment ratio diagram for the three (3) sub-regions. (a) Region 1; (b) Region 2; (c) Region 3.
regional average L-skewness and L-kurtosis fitted to varying probability distributions for all three regions. The 3-parameter distribution line/curve closest to L-moment ratio points of sub-regional sites portrays an initial deduction concerning an optimal distribution  ; in this case, the Generalized Logistic (GLO) curve satisfies this approximation. A 3-parameter distribution is selected instead of its 2-parameter counterpart due to its robustness and ability to optimally represent the probability distribution parameters .
In addition to L-Moment ratio diagram, Z Statistics provides a viable statistical approach to identify the optimal distribution that best fits the data for each sub-region. Table 5 shows the Z Statistics for all distributions for each sub-region and reveals that GLO is the most significant at a 90% confidence interval (Z ≤ |1.64|) as prescribed by Hosking and Wallis  for regions 2 and 3. The result is consistent with deductions from the L-Moment Ration diagram.
The identified optimal probability distribution corresponds with those applied in previous single-site and regional studies undertaken for catchment areas close to our study area  . The insignificance of the probability distribution for all combined sites and region 1 (Z > 1.65) suggests that all individual sites within this region are not defined by the same particular distribution due to their apparent heterogeneity.
4.4. Regional Flood Frequency and Parameter Estimation
After identifying GLO as the optimal probability distribution for regions 2 and 3, a flood frequency relationship was established to derive flood magnitudes. The GLO probability distribution function is given by:
Table 5. Region clusters and Z statistics for a selection of an optimal probability distributions.
a = optimal distribution.
where ξ, α and k depict the location, scale and shape parameters, respectively .
The range of x is defined as if ; If ; if .
The location parameter (ξ) dictates the position of the distribution about a symmetric axis, the scale parameter (α) define the distribution spread, and the shape parameter (k) indicates the behaviour of the upper tail of the distribution. These parameters are derived from L-moments and applied to derive T-year flood exceedance based on the GLO (XT) defined by:
where: , T is the return period and ZT is the growth curve of T.
GLO distribution parameters estimated for each sub-region using L-moments were substituted into Equation (3) to estimate the sub-regional growth factors for ungauged and sparsely gauged basins and are presented in Table 6.
4.5. Climate Indices Correlation with Peak Annual Discharge
Ijaka-Oke, Oba/Oyo-Ogbomosho, Ofiki/Igangan-Ilere road and Ofiki-Igangan were identified by breakpoints and trend analysis to be heterogeneous and were further investigated to ascertain the influence of climate variability by correlating peak annual discharge and global climate indices. Regional and direct flood frequency estimates were then determined in ICI-RAFT using the highest correlated indices.
The Madden-Julian Oscillation (MJO) demonstrated the highest correlation with the annual maximum flood time series for the four sites (Figure 4) using an optimal lag time of 1 month (Climate Index Lag), considering that only single peak discharge values for each year were applied in the analysis. Correlation coefficients (R2) based on MJO phase 7 which originates from the Western
Table 6. Region clusters and parameters for the selected distribution functions.
Figure 4. Plots showing the relationship between climate-indices and peak annual discharge for selected gauging stations. (a) Ijaka-Oka; (b) Ofiki-Igangan; (c) Ofiki/Igangan-Ilere road; (d) Oba/Oyo-Ogbomosho.
Pacific (longitude 40W) were 0.27, 0.50 0.31 and 0.45 for Ijaka-Oke, Ofiki Igangan, Ofiki/Igangan-Ilere road and Oba/Oyo-Ogbomosho, respectively, suggesting that between 27 to 50 percent of the variability in the annual maximum flood series is attributable to climate dynamics. The correlation values derived in this study were consistent with those revealed in other studies . The remaining variability in peak annual discharge can be linked factors such as local catchment properties, land use/cover changes and hydrodynamics  , which is beyond the scope of this study.
MJO is known to be a strong driver of rainfall variability in tropical regions , governing atmospheric pressure and temperature around the equator. The MJO significantly influences regional rainfall   and was reported to have influenced the rainfall dynamics that triggered the unprecedented 2012 flood event in Nigeria . Arnold et al.,  and Caballero and Huber,  further suggested that, due to the dependence of MJO on Sea Surface Temperature (SST) and Outgoing Longwave Radiation (OLR), MJO activities may increase in response to global warming, resulting in more frequent MJO-influenced flood events.
4.6. Climate Variability Effect and Flood Quantile Estimation
Results presented in Table 7 and Figure 5 show flood frequency estimates derived when all available data points for the hydrological time series are used for regional and direct flood frequency estimation, as well as when data points that correlate and satisfies the limits set for the MJO climate variability index are used, thus accounting for the influence of climate variability. These results suggest
Table 7. Flood frequency estimates per return period (regional and Direct, and considering climate variability).
Unit of flood magnitude = m3/s.
Figure 5. Flood frequency plots of regional and direct flood frequency analysis taking climate variability into account. (a) Ijaka-Oke; (b) Ofiki-Igangan; (c) Ofiki/Igangan-Ilere road; (d) Oba/Oyo-Ogbomosho.
that climate variability accounted for increased flood magnitude for regional and direct flood estimates varying from 0% - 35%, and demonstrate that multi-decadal changes in ocean-atmosphere conditions can influence both small and larger floods . Also, the influence of climate variability was most evident at sites that exhibited a higher correlation with HCI (i.e. Ofiki Igangan, Ofiki/Igangan-Ilere road and Oba/Oyo-Ogbomosho). These results are generally consistent with those revealed by other studies where flood estimates that accounted for climate variability were higher than those estimated under the assumption of stationarity .
The criss-cross plot pattern observed at Ijaka-Oke for the inclusion of climate variability in regional flood frequency estimation suggests that caution must be taken when integrating climate variability into FFA , especially when the relationship between climate indices is low (R2 = 0.28). Additionally, the significance of the homogeneity rather than trends is identified as the key indicator of nonstationarity , as evident at the Ijaka-Oke gauging station (p-values: trend = 0.001, homogeneity = 0.081).
Furthermore, Figure 4 reveals that the historical maximum flood experienced at each site in the OORB is less than the 1-in-100 year flood guideline stipulated for flood management planning in Nigeria . Thereby reinforcing the need to implement flood management measures (both structural and non-structural) based on a 1-in-100 year to curtail recurring flood impacts.
The impact of flooding in Nigeria has increased over the last two decades, resulting in the displacement of persons, disruption of socio-economic activities, damage of infrastructures and loss of lives. Therefore, efficient flood risk management is urgently needed to reduce the vulnerability and exposure of the local population and assets. Flood frequency analysis is usually the first step towards flood risk management, aimed to determine flood magnitudes for varying return periods. This is, however, challenging in many developing and remote regions, including the Ogun-Osun River Basin of Nigeria due to the absence of sufficient historic hydrological data, due to financial, technical and organizational drawbacks.
We have presented a robust flood estimation approach based on L-moment regional flood frequency analysis that combines multiple short-length historical hydrological data to curb aleatoric uncertainty. Building on the evidence of climate influence on the changing hydrological regime in the region, this study accounts for climate variability effect on flood frequency estimates through climate indices.
Two homogeneous regions are identified based on clustering algorithm and statistical tests, while GLO three-parameter distribution is identified as the best fit distribution for flood frequency analysis in the 2 sub-regions based on L-moment ratio diagrams and goodness of fit test (Z-statistic). Also, Madden-Julian Oscillation (MJO) is identified as the most influential climate indices for the region, resulting in increased flood magnitude for direct and regional flood frequency estimates. This further reinforces the need to integrate climate variability into flood frequency analysis as more climate-driven events are expected due to global warming.
In conclusion, integrating climate variability into regional and direct flood frequency estimation results in more robust outcomes—a useful input into hydraulic modelling and flood risk mapping needed to inform resilient structural and non-structural flood risk management interventions. However, for the result of this study to be transferable to ungauged areas within the homogeneous regions, further analysis is needed to determine the physiographical parameters required to establish the best relationship with the index flood of sites in homogeneous sub-regions; for example, using approaches such as artificial neural networks (ANN)
The authors acknowledge the Niger Delta Development Commission (NDDC), Nigeria for funding Iguniwari Thomas Ekeu-wei’s Ph.D. at Lancaster University, United Kingdom (NDDC/DEHSS/2013PGFS/BY/5), from which this paper is a product; we also acknowledge The Ogun Osun River Basin Authority (OORBA) for providing the in-situ river hydrological data. We appreciate the US Army Corps of Engineers provided free ICI-RAFT software license used for flood frequency analysis.
 Tehrany, M.S., Pradhan, B. and Jebur, M.N. (2014) Flood Susceptibility Mapping Using a Novel Ensemble Weights-of-Evidence and Support Vector Machine Models in GIS. Journal of Hydrology, 512, 332-343.
 The Federal Government of Nigeria (2013) Post-Disaster Needs Assessment 2012 Floods.
 Agada, S. and Nirupama, N. (2015) A Serious Flooding Event in Nigeria in 2012 with Specific Focus on Benue State: A Brief Review. Natural Hazards, 77, 1405-1414.
 Kjeldsen, T.R., Smithers, J.C. and Schulze, R.E. (2002) Regional Flood Frequency Analysis in the KwaZulu-Natal Province, South Africa, Using the Index-Flood Method. Journal of Hydrology, 255, 194-211.
 Mishra, B., Takara, K., Yamashiki, Y. and Tachikawa, Y. (2009) Hydrologic Simulation-Aided Regional Flood Frequency Analysis of Nepalese River Basins. Journal of Flood Risk Management, 2, 243-253.
 Ampadu, B., Chappell, N.A. and Kasei, R.A. (2013) Rainfall-River Flow Modelling Approaches: Making a Choice of Data-Based Mechanistic Modelling Approach for Data Limited Catchments: A Review. Canadian Journal of Pure and Applied Sciences, 7, 2571-2580.
 Dano Umar, L., et al. (2011) Geographic Information System and Remote Sensing Applications in Flood Hazards Management: A Review. Research Journal of Applied Sciences, Engineering and Technology, 3, 933-947.
 O’Brien, N.L. and Burn, D.H. (2014) A Nonstationary Index-Flood Technique for Estimating Extreme Quantiles for Annual Maximum Streamflow. Journal of Hydrology, 519, 2040-2048.
 Hounkpè, J., Diekkrüger, B., Badou, D.F. and Afouda, A.A. (2015) Non-Stationary Flood Frequency Analysis in the Ouémé River Basin, Benin Republic. Hydrology, 2, 210-229.
 Li, J. and Tan, S. (2015) Nonstationary Flood Frequency Analysis for Annual Flood Peak Series, Adopting Climate Indices and Check Dam Index as Covariates. Water Resources Management, 29, 5533-5550.
 Adeaga, O., Oyebande, L. and Depraetere, C. (2006) Surface Runoff Simulation for Part of Yewa Basin. In: Predictions in Ungauged Basins: Promise and Progress, No. 303, IAHS Publ., Wallingford, 382.
 Adeleke, O.O., Makinde, V., Eruola, A.O., Dada, O.F., Ojo, A.O. and Aluko, T.J. (2015) Estimation of Groundwater Recharges Using Empirical Formulae in Odeda Local Government Area, Ogun State, Nigeria. Challenges, 6, 271-281.
 Mouhamed, L., Traore, S.B., Alhassane, A. and Sarr, B. (2013) Evolution of Some Observed Climate Extremes in the West African Sahel. Weather and Climate Extremes, 1, 19-25.
 Diatta, S. and Fink, A.H. (2014) Statistical Relationship between Remote Climate Indices and West African Monsoon Variability. International Journal of Climatology, 34, 3348-3367.
 Salau, O.R., Fasuba, A., Aduloju, K.A., Adesakin, G.E. and Fatigun, A.T. (2016) Effects of Changes in ENSO on Seasonal Mean Temperature and Rainfall in Nigeria. Climate, 4, 5.
 Ewemoje, T.A. and Ewemooje, O. (2011) Best Distribution and Plotting Positions of Daily Maximum Flood Estimation at Ona River in Ogun-Oshun River Basin, Nigeria. Agricultural Engineering International: CIGR Journal, 13, 1-10.
 GCOS-AOPC/PPOC. Download Climate Time Series.
 Ekeu-wei, I.T., Blackburn, G.A. and Pedruco, P. (2018) Infilling Missing Data in Hydrology: Solutions Using Satellite Radar Altimetry and Multiple Imputation for Data-Sparse Regions. Water, 10, 1483.
 Yue, S. and Wang, C. (2002) The Influence of Serial Correlation on the Mann– Whitney Test for Detecting a Shift in Median. Advances in Water Resources, 25, 325-333.
 Pedruco, P., Nielsen, C., Kuczera, G. and Rahman, A. (2014) Combining Regional Flood Frequency Estimates with an at Site Flood Frequency Analysis Using a Bayesian Framework: Practical Considerations. In: Hydrology and Water Resources Symposium, ACT: Engineers Australia, Barton, 766-773.
 Grubbs, F.E. and Beck, G. (1972) Extension of Sample Sizes and Percentage Points for Significance Tests of Outlying Observations. Technometrics, 14, 847-854.
 Kwon, H.H., Brown, C. and Lall, U. (2008) Climate Informed Flood Frequency Analysis and Prediction in Montana Using Hierarchical Bayesian Modeling. Geophysical Research Letters, 35, L05404.
 Gutiérrez, F. and Dracup, J.A. (2001) An Analysis of the Feasibility of Long-Range Streamflow Forecasting for Colombia Using El Ni?o-Southern Oscillation Indicators. Journal of Hydrology, 246, 181-196.
 Sayers, P., et al. (2015) Strategic Flood Management: Ten “Golden Rules” to Guide a Sound Approach. International Journal of River Basin Management, 13, 137-151.
 Machado, M.J., Botero, B.A., López, J., Francés, F., Díez-Herrero, A. and Benito, G. (2015) Flood Frequency Analysis of Historical Flood Data under Stationary and Non-Stationary Modelling. Hydrology and Earth System Sciences Discussions, 12, 525-568.
 Padi, P.T., Baldassarre, G.D. and Castellarin, A. (2011) Floodplain Management in Africa: Large Scale Analysis of Flood Data. Physics and Chemistry of the Earth, 36, 292-298.
 Stedinger, J.R. and Griffis, V.W. (2008) Flood Frequency Analysis in the United States: Time to Update. Journal of Hydrologic Engineering, 13, 199-204.
 Izinyon, O. and Ehiorobo, J. (2014) L-Moments Approach for Flood Frequency Analysis of River Okhuwan in Benin-Owena River Basin in Nigeria. Nigerian Journal of Technology, 33, 10-18.
 Komi, K., Amisigo, B.A., Diekkrüger, B. and Hountondji, F.C. (2016) Regional Flood Frequency Analysis in the Volta River Basin, West Africa. Hydrology, 3, 5.
 Peel, M., Wang, Q.J., Vogel, R. and McMahon, T. (2001) The Utility of L-Moment Ratio Diagrams for Selecting a Regional Probability Distribution. Hydrological Sciences Journal, 46, 147-155.
 Hailegeorgis, T.T. and Alfredsen, K. (2017) Regional Flood Frequency Analysis and Prediction in Ungauged Basins Including Estimation of Major Uncertainties for Mid-Norway. Journal of Hydrology: Regional Studies, 9, 104-126.
 Madden, R.A. and Julian, P.R. (1971) Detection of a 40-50 Day Oscillation in the Zonal Wind in the Tropical Pacific. Journal of the Atmospheric Sciences, 28, 702-708.
 Mohino, E., Janicot, S., Douville, H. and Li, L. (2012) Impact of the Indian Part of the Summer MJO on West Africa Using Nudged Climate Simulations. Climate Dynamics, 38, 2319-2334.
 ACMAD (2012) Flood Report over West Africa—September 2012. African Centre of Meteorological Applications for Development (ACMAD).
 Arnold, N., Branson, M., Kuang, Z., Randall, D. and Tziperman, E. (2015) MJO Intensification with Warming in the Superparameterized CESM. Journal of Climate, 28, 2706-2724.
 Sheng, H., et al. (2020) Frequency and Magnitude Variability of Yalu River Flooding: Numerical Analyses for the Last 1000 Years. Hydrology and Earth System Sciences.
 López, J. and Francés, F. (2013) Non-Stationary Flood Frequency Analysis in Continental Spanish Rivers, Using Climate and Reservoir Indices as External Covariates. Hydrology and Earth System Sciences, 17, 3189-3203.
 Serinaldi, F., Kilsby, C.G. and Lombardo, F. (2018) Untenable Nonstationarity: An Assessment of the Fitness for Purpose of Trend Tests in Hydrology. Advances in Water Resources, 111, 132-155.