The aim of the study is to present a combination of techniques for (a) the spatiotemporal analysis of mean monthly gridded precipitation datasets and (b) the evaluation of the relative position of the existing rain-gauge network. The mean monthly precipitation (P) patterns of Nigeria using ∼1 km2 grids for the period 1950–2000 were analyzed and the position of existing rain-gauges was evaluated. The analysis was performed through: (a) correlations of P versus elevation (H), latitude (Lat) and longitude (Lon); (b) principal component analysis (PCA); (c) Iso-Cluster and maximum likelihood classification (MLC) analysis for terrain segmentation to regions with similar temporal variability of mean monthly P; (d) use of MLC to create reliability classes of grid locations based on the mean clusters’ characteristics; and (e) analysis to evaluate the relative position of 33 rain-gauges based on the clusters and their reliability classes. The correlations of mean monthly P versus H, Lat, Lon, and PCA highlighted the spatiotemporal effects of the Inter Tropical Discontinuity phenomenon. The cluster analysis revealed 47 clusters, of which 22 do not have a rain-gauge while eight clusters have more than one rain-gauge. Thus, more rain-gauges and a better distribution are required to describe the spatiotemporal variability of P in Nigeria.

Over the last years, high resolution grids of various climatic parameters have been developed at global scale (e.g., Hijmans et al. 2005; Mitchell & Jones 2005; Sheffield et al. 2006) providing significant information for regions with no meteorological stations. The availability of such data stimulates interest for the implementation of various statistical techniques, which have already been used to capture the spatial and seasonal patterns of precipitation or other climatic parameters using data directly from meteorological stations (Wackernagel 2003; Chen et al. 2009; Buttafuoco et al. 2011; Marchetti et al. 2015).

The development of gridded climatic data is usually performed either using climatic models (e.g., general circulation models) (Sheffield et al. 2006; Shaffrey et al. 2009; Watanabe et al. 2010) in combination with downscaling techniques (Wilby & Wigley 1997), or using interpolation techniques (Ninyerola et al. 2000; Boer et al. 2001; Hijmans et al. 2005). Advanced interpolation techniques take into account not only the site-specific observations from stations and their location but also the effects of topography and other parameters (Ninyerola et al. 2000; Hijmans et al. 2005), while climatic models include also the effects of ocean and air masses’ circulation (Sheffield et al. 2006; Shaffrey et al. 2009; Watanabe et al. 2010).

Gridded multi-temporal climatic datasets can be used in order to define regions of similar spatiotemporal variability based on one or more climatic parameters using techniques such as cluster analysis in GIS environment. Such techniques have been used on land surface temperature (LST) gridded data in an attempt to define sub-regions with different seasonal LST variability, to assess its sensitivity to climatic change and to support environmental analysis (Miliaresis 2009, 2012; Miliaresis & Partsinevelos 2010; Maeda & Hurskainen 2014). In addition, Miliaresis (2013, 2014) proposed a method for standardizing multi-temporal LST imagery in order to reveal and describe thermal anomalies using elevation, latitude and longitude dependencies using correlation and principal component analysis (PCA). These techniques have also been adopted in order to analyze the spatiotemporal patterns of precipitation and reference crop evapotranspiration allowing the terrain segmentation of a territory based on the spatiotemporal variation of the two climatic parameters (Demertzi et al. 2014; Aschonitis et al. 2016).

The climatic regionalization/segmentation of a territory based on the spatiotemporal variation of a gridded climatic parameter and its dependency on spatial features (e.g., elevation, latitude, longitude, etc.) can further be used to evaluate the representativeness of existing meteorological stations to describe the climatic variability of adjacent territories with no observed data. This information could be used as an additional tool for the evaluation of existing networks of meteorological stations, especially for rain-gauge stations, since precipitation is the most variable and one of the most important climatic parameters. The design of the spatial structure/distribution of rain-gauge networks is very complex and efforts regarding its evaluation have been performed using satellite images (Caselles & Melia 1985), radar observations (Abdella & Alfredsen 2010), geostatistical methods and optimization algorithms (Filippini et al. 1994; Papamichail & Metaxa 1996; Tsintikidis et al. 2002; Nour et al. 2006; Adhikary et al. 2015), probabilistic approaches (Shafiei et al. 2014), ‘cost–benefit’ approaches using a step-by-step reduction of stations’ density (Vose & Menne 2004; Horel & Dong 2010), entropy approaches (Krstanovic & Singh 1992; Al-Zahrani & Husain 1998; Kawachi et al. 2001; Chen et al. 2008; Yoo et al. 2008), approaches based on intensity–duration–frequency (Chebbi et al. 2013), and combinations of the aforementioned approaches (Pardo-Igúzquiza 1998; Barca et al. 2008; Cheng et al. 2008). Considering the above, terrain segmentation/clustering approaches based on available gridded data could assist the aforementioned procedures for the evaluation and design of rain-gauge networks.

The aim of the study is to present a combination of techniques for spatiotemporal analysis of multi-temporal gridded climatic datasets, which can further be used to evaluate the relative position of existing meteorological stations. The techniques were applied for the analysis of gridded mean monthly precipitation data for the period 1950–2000 and for the evaluation of the existing rain-gauge network of Nigeria.

Study area and climate

The study area is the country of Nigeria which is confined between 4 ° to 14 ° north and 2 ° to 15 ° east covering a total area of ∼923,770 km2. The selection of the specific study area was based on the fact that the climate of the country presents extremely high latitudinal variability, which ranges from south to north as follows: tropical monsoon climate (Am), tropical savanna climate (Aw), hot semi-arid climate (BSh), hot desert climate (BWh) (Peel et al. 2007) (Figure S1(a) in Supplementary material, available with the online version of this paper). Nigeria has only two seasons, a dry and a rainy/monsoon season. The dry season is influenced by dry north-easterly winds, which regulate the flow of dust-laden air masses from the Sahara Desert, known as Continental Tropical (cT) air mass (locally known as Harmattan). The monsoon season is influenced by moist south-westerly winds that regulate the flow of air masses known as the Maritime Tropical (mT) air mass which reaches the country from the Gulf of Guinea coast (tropical Atlantic). The moving boundary between the two wind systems creates a zone of moisture discontinuity known as the Inter Tropical Discontinuity (ITD) zone. mT is responsible for Nigeria's rainy season, invades the country usually after February from the south and reaches the northern part of Nigeria during June. mT invasion is followed by the northward retreat of cT and both of them control the northward shift of ITD. Precipitation variability in Nigeria is not only regulated by the ITD mechanism but also by other tropical and extra-tropical factors, such as the tropical easterly jet (TEJ), the sea surface temperature (SST) anomaly, the stratospheric Quasi-Biennial Oscillation (QBO), the El Nino Southern Oscillation (ENSO), Southern Oscillation Index (SOI) and North Atlantic Oscillation (NAO), and the Atlantic Multi-Decadal Oscillation (AMO) (Olaniran 1983, 1988a, 1988b, 1990, 1991a, 1991b; Adefolalu & Oguntoyinbo 1985; Adefolalu 1988; Olaniran & Sumner 1989a, 1989b, 1990; Bello 1996, 1998; Odekunle 2010; Oguntunde et al. 2011, 2014; Alli et al. 2012; Ogungbenro & Morakinyo 2014).

Data

Gridded precipitation data for Nigeria were obtained from the WorldClim database (http://www.worldclim.org/). The database provides mean monthly values of precipitation at 30 arc-sec spatial resolution (∼1 × 1 km) for the period 1950–2000 at global scale (Hijmans et al. 2005). The dataset was produced using the thin-plate smoothing spline algorithm implemented in the ANUSPLIN package for interpolation, using latitude, longitude, and elevation as independent variables and mean monthly precipitation data from rain-gauge stations (Hijmans et al. 2005). The authors performed: (a) uncertainty quantification arising from the input data and the interpolation by mapping weather station density; (b) analysis of elevation bias in the weather stations and of elevation variation within grid cells; and (c) data partitioning and cross validation. A revised version of the GTOPO30 DEM based on SRTM DEM at 30 arc-sec spatial resolution was also obtained by the WorldClim database.

Location data of the current network of 33 rain-gauges in Nigeria were obtained by NOAA (www.nws.noaa.gov/tg/siteloc.php). The position of rain-gauge stations in comparison to (a) Köppen climate classification (Peel et al. 2007), (b) different ecoregions (Olson et al. 2001) and hydrologic network, and (c) elevation are given in Figure S1(a)–S1(c), respectively.

Methods of imagery analysis

A combination of three statistical methods consisting of correlation, principal components (PCs), and cluster analysis was used in GIS environment (Arc-Info 9.3, ESRI-GIS) for imagery analysis. Correlation analysis and PCA were used in order to investigate the spatiotemporal dependencies of P according to the proposed procedures by Miliaresis (2012) while cluster analysis was used to perform terrain segmentation based on the spatiotemporal variation of mean monthly precipitation.

The correlation analysis derived the Pearson correlation coefficient R of mean monthly precipitation P versus elevation H, latitude Lat, and longitude Lon. The analysis was performed using the following toolbox in Arc-Map (Spatial Analyst>Multivariate>Band Collection Statistics) which produces a text file with the correlations’ table.

PCA is a linear transformation technique that produces a set of images known as PCs that are uncorrelated with one another while they are ordered in terms of the amount of variance they explain from the original set of images (Maaten & Hinton 2008). PCs are computed from the linear combination of eigenvectors and the corresponding pixel values of the initial images (Mather 2004). PCA has traditionally been used in remote sensing as a means of data compaction since it is common to find that the first two or three components are able to explain the majority of the variability in data values, while later components tend to be dominated by noise effects. The rejection of these later components reduces the volume of data with no appreciable loss of information (Siljestrom et al. 1997; Miliaresis 2013). Standardized PCA (Eastman & Fulk 1993) was applied in this study (data of P per month is centered with mean 0.0 and standard deviation 1.0) and so each image is not weighed according to its variance. The standardization was performed with the following formula:
formula
(1)
where Psi is the raster of standardized precipitation values for the month i, Pi is the raster of real precipitation values for the month i, Pimax is the maximum real value of precipitation in the Pi raster, and Pimin is the minimum real value of precipitation in the Pi raster. The standardization of the rasters was performed in the Raster Calculator while the PCA analysis was performed using the 12 monthly standardized rasters of monthly precipitation and the following toolbox in Arc-Map (Spatial Analyst>Multivariate>PCs). The specific PCA analysis produces 12 PC rasters and a text file with their monthly factor loadings and the % of variance they explain.

Cluster analysis was used to segment Nigeria into a number of distinct territories (clusters) with similar intensity and temporal variability of mean monthly precipitation. Cluster signatures were derived by the Iso Cluster algorithm (Ball & Hall 1965), which uses a modified iterative optimization procedure known as the migrating means technique. The algorithm separates all cells into a predefined number of distinct unimodal groups (clusters) in the multidimensional space of the input bands. The algorithm iteratively computes the minimum Euclidean distance when assigning each candidate cell to a cluster. The process starts with arbitrary means assigned by the software for each cluster (the number of clusters is given by the user). Every cell is assigned to the closest of these means (all in the multidimensional attribute space). New means are iteratively recalculated for each cluster based on the attribute distances of the cells that belong to the cluster (200 iterations with sampling interval of 10 pixels were used in this study). The optimal number of clusters to specify is usually unknown. For this reason, an initial trial is performed using a high number in order to allow the algorithm to give its maximum number of stable clusters, which is smaller than the specified one by the user (Ball & Hall 1965; Richards 1986). The user can choose the maximum number of stable clusters or to reduce their number by repeating the procedure. In our case, the maximum number was used while the value of 50 km2 was assigned as the minimum threshold of total cluster area coverage. The aforementioned procedure was performed using the following toolbox in Arc-Map (Spatial Analyst>Multivariate>Iso Cluster) which produces a signatures file of variances and covariances of each cluster. The maximum likelihood classification (MLC) was used for terrain segmentation using the cluster signatures from the Iso Cluster. The MLC algorithm is based on Bayes’ theorem of decision (Gelman et al. 2014) considering that the cells in each cluster sample in the multidimensional space are normally distributed. The MLC considers both the variances and co-variances of the cluster signatures when assigning each cell to one of them. Each cluster is characterized by a mean vector and a covariance matrix. Taking into account these two characteristics, statistical probabilities are computed for the pixels to determine their proximity to the mean characteristics of the clusters to which they belong. Thus, MLC produced in our case a raster which shows the levels of % confidence of classification reliability for each classified pixel (eight reliability classes were used). A 0% reject fraction was used in order not to allow unclassified pixels. The aforementioned procedure was performed using the following toolbox in Arc-Map (Spatial Analyst>Multivariate>MLC) which produces two rasters: the raster of clusters classification and the raster of reliability classification classes.

Finally, the distribution of the 33 existing rain-gauges’ network was evaluated based on (a) the number of rain-gauges which exist in each cluster and (b) the reliability class values of the pixels in the positions of rain-gauges.

Spatial variation of mean annual precipitation

The 30 arc-sec resolution (∼1 × 1 km) maps of elevation H and mean annual precipitation P for the period 1950–2000 are given in Figure 1(a) and 1(b), respectively, while the respective cumulative area (based on the number of pixels) for the observations of elevation H and P in Figure 1(a) and 1(b) are given, respectively, in Figure 1(c) and 1(d). The respective overall mean annual value of precipitation and elevation over Nigeria are 1,174 ± 547 mm yr−1 and 329 ± 222 m above the mean sea level (m.a.s.l).

Figure 1

Maps of (a) elevation H (meters above sea level), (b) mean annual precipitation P (mm year−1) for the period 1950–2000, (c) cumulative area (based on number of pixels) for the observations of elevation H (raster values of (a)), and (d) cumulative area for the observations of mean annual precipitation P for the period 1950–2000 (raster values of (b)).

Figure 1

Maps of (a) elevation H (meters above sea level), (b) mean annual precipitation P (mm year−1) for the period 1950–2000, (c) cumulative area (based on number of pixels) for the observations of elevation H (raster values of (a)), and (d) cumulative area for the observations of mean annual precipitation P for the period 1950–2000 (raster values of (b)).

Close modal

The rainfall map (Figure 1(b)) shows that the highest annual P values are observed along the southern coastline while the lowest annual P values were observed in the northern part of the country. These trends indicate a clear decrease in rainfall amount with the increase of latitude verifying the results of previous studies (Bello 1998; Nicholson et al. 2000; Oguntunde et al. 2011). Taking into account the 1,077,730 pixel values of mean annual P raster, a simple regression was performed with latitude values in order to assess the latitudinal effect on annual P variability. The best fit was performed by a reciprocal-Lat function of P (P= 16,285.1/Lat ‒ 641.9, p < 0.001, R2 = 0.89) (Figure S2(a) in Supplementary material, available with the online version of this paper). The regression showed that the slope, which describes the decrease of annual precipitation per latitude unit increase, is steeper in the low latitude areas. The mean latitudinal effect for the entire country was found equal to 212 mm of P reduction per 1 ° of latitude increase. Splitting Nigeria into two regions between 4 °–9 ° N and 9 °–14 ° N the respective rates were found to be 468 and 130 mm of P reduction per 1 ° of latitude increase. The graphs between mean annual P versus longitude Lon and elevation H are also given in Figure S2(b) and S2(c), respectively, without providing fitting analysis due to weak correlation coefficients.

The 30 arc-sec resolution (∼1 × 1 km) maps of mean annual precipitation P for the period 1950–2000 showed that the overall mean annual value of precipitation over Nigeria is 1,174 ± 547 mm yr−1. Oguntunde et al. (2011), who also analyzed other gridded precipitation data of coarser resolution (0.5 °, CRU TS 2.1 dataset) (Mitchell & Jones 2005) for the period 1901–2000, found that the mean ± st. dev. precipitation is equal to 1,170 ± 109 mm yr−1. Although the mean annual precipitations given by the two different datasets are in proximity, the difference in their standard deviations is extremely high. The significantly lower standard deviation of the second dataset suggests a significant loss of detail in local precipitation patterns due to the lower resolution. This loss of detail in the general climatology of the country has also been identified by Hassan et al. (2009).

Spatial and seasonal variation of mean monthly precipitation

The minimum, maximum, and mean ± st. dev. values of monthly precipitation for the studied period are also presented in Figure 2(a). According to Figure 2(a), the minimum and mean monthly precipitations show a unimodal response while the maximum monthly precipitation shows a bimodal response with two maximum extremes during June (641 mm) and September (527 mm).

Figure 2

(a) Minimum, maximum, and mean ± st. dev. values of monthly precipitation in Nigeria for the period 1950–2000 and (b) correlations of the mean monthly rainfall P versus elevation H, latitude Lat, and longitude Lon for the period 1950–2000.

Figure 2

(a) Minimum, maximum, and mean ± st. dev. values of monthly precipitation in Nigeria for the period 1950–2000 and (b) correlations of the mean monthly rainfall P versus elevation H, latitude Lat, and longitude Lon for the period 1950–2000.

Close modal

The correlations of the mean monthly values of precipitation P versus elevation H, latitude Lat, and longitude Lon are given in Figure 2(b). It is indicative that the monthly variation of correlations of P versus H, Lat, and Lon follow similar patterns. All correlations are negative showing a relative stability during the period of October–June. The correlations’ magnitude of P versus these three parameters shows the following order Lat > H > Lon during this period. On the other hand, the period of July–September indicates a climatic disturbance, mainly due to ITD northward movement, which forces P-H, P-Lat, P-Lon relations to shift towards positive values. This disturbance succeeds in inverting the effects of H and Lon on monthly P variation in August (positive R) but not the effect of Lat, although it shows more intense change in comparison to H and Lon. This shift is also responsible for the bimodal precipitation patterns in the southern regions during the wet season. It is indicative that the monthly variation of absolute correlation for P-Lat is maximized during the onset (April–May) and cessation (October) of rainfall reaching values |R| > 0.9. The negative correlation between monthly P and elevation H is probably related to the fact that the higher altitude areas are located in the drier northern regions.

In order to further explore the monthly correlations observed in Figure 2(b), PCA was performed and the corresponding eigenvalues and factor loadings are presented in Table S1 (available with the online version of this paper). The first two PCs (PC-1 and PC-2) account for 91.7% of the total variance in the mean monthly P and they are given in Figure 3. The PC-1 map (Figure 3(a)) explains 80.3% of the total variance in the mean monthly P and amplifies the difference between August and the rest of the period of September–July (Table S1 and Figure 3(c)). The factor loadings of PC-1 (Figure 3(c)) are relatively stable during the period of September–July while during August they are significantly reduced. PC-1 is spatially maximized in the southern tropical zone while it is minimized in northern Nigeria, providing a first indication of which regions are more likely to present higher amounts of rainfall. PC-2 (Figure 3(b)) explains 11.4% of the total variance in the mean monthly P and amplifies the difference between the periods of October–May and June–September. The difference between these two periods is spatially maximized in the central mountainous regions of Plateau and Kaduna provinces and in the southeastern regions of the coastline. PC-2 is minimized in the southwestern region (rain forests) of Nigeria and in the extreme northeastern region (Figure 3(b)). PC-1 and PC-2 (Figure 3) provide a good visualization of the northward movement effects of IDT on precipitation patterns and the period of occurrence (period of June–September with a peak in August according to Figure 3(c)).

Figure 3

The first two PCs (a) PC-1 (80.3%), (b) PC-2 (11.4%) that account for 91.7% of the P variance, and (c) the monthly variation of the PC-1 and PC-2 factor loadings according to Table S1 (available with the online version of this paper).

Figure 3

The first two PCs (a) PC-1 (80.3%), (b) PC-2 (11.4%) that account for 91.7% of the P variance, and (c) the monthly variation of the PC-1 and PC-2 factor loadings according to Table S1 (available with the online version of this paper).

Close modal

Terrain segmentation based on the spatial and seasonal variation of mean monthly rainfall

Cluster analysis revealed 47 clusters and their characteristics such as percent area coverage, mean elevation and mean monthly P values are given in Table 1 while their spatial distribution is given in Figure 4(a). The levels of % confidence of classification reliability for each classified pixel are also given in Figure 4(b).

Table 1

Characteristics of the 47 precipitation clusters such as % coverage, mean monthly and mean annual values of P, and classification based on unimodal or bimodal precipitation patterns

Cluster no.% coverageMean elevation (m)JanFebMarAprMayJunJulAugSepOctNovDecAnnualModality subgroup
3.7% 313 11 31 100 145 49 344 1.1 
4.3% 325 20 56 137 177 74 10 479 1.1a 
5.2% 343 31 74 159 205 95 13 584 1.1b 
5.0% 379 10 41 89 178 233 114 17 682 1.1a 
4.7% 420 14 55 105 200 259 127 20 781 1.1b 
5.0% 361 23 69 116 189 243 157 31 837 1.1 
3.7% 505 22 77 126 224 275 154 28 910 1.1 
2.6% 664 37 106 151 245 283 188 41 1,058 1.1a 
4.7% 396 32 87 134 196 254 188 44 943 1.1a 
10 4.4% 242 49 111 152 193 230 203 66 1,013 1.1b 
11 3.8% 465 15 66 150 178 264 293 230 65 1,267 1.1a 
12 1.9% 952 42 105 150 205 271 232 61 1,078 1.1a 
13 2.5% 569 26 78 154 186 256 314 289 118 1,438 1.1 
14 5.3% 242 18 69 135 170 191 213 234 101 1,142 1.2 
15 2.9% 161 33 88 162 175 190 211 249 138 12 1,269 1.2 
16 4.0% 373 11 55 127 168 212 261 267 96 1,204 1.2b 
17 3.3% 376 27 83 154 180 223 253 259 125 1,320 1.2a 
18 1.9% 814 11 61 126 181 208 256 260 287 183 24 1,603 1.2 
19 0.8% 89 13 31 92 134 177 221 174 86 159 162 48 14 1,312 2.3 
20 1.7% 239 23 80 124 153 167 144 104 189 160 35 12 1,199 2.2b 
21 1.7% 310 16 60 106 145 161 154 128 206 138 20 1,152 2.2 
22 1.6% 145 28 80 144 172 190 184 235 111 1,164 2.2a 
23 2.9% 299 11 46 96 152 171 167 156 231 133 15 1,192 2.2b 
24 2.1% 399 22 76 118 159 173 167 138 234 164 31 15 1,306 2.2 
25 0.9% 130 24 78 149 175 166 198 227 132 1,161 2.1a 
26 1.1% 145 35 95 176 194 181 209 251 161 14 1,327 2.1a 
27 1.4% 199 17 55 108 178 189 186 175 270 174 27 1,394 2.2 
28 1.1% 162 45 106 193 214 191 212 269 197 23 1,471 2.1 
29 1.1% 176 10 22 78 125 200 216 222 187 291 205 36 1,602 2.2 
30 0.9% 177 13 56 117 208 237 209 226 286 235 33 1,636 2.1 
31 1.0% 217 32 96 135 169 205 217 147 239 182 49 17 1,498 2.2a 
32 0.7% 94 10 38 105 146 189 245 275 179 280 200 57 20 1,745 2.2 
33 0.5% 19 21 38 93 129 227 329 222 77 162 174 63 18 1,553 2.3b 
34 0.7% 34 16 38 97 133 210 294 271 115 225 201 69 18 1,686 2.3 
35 0.9% 132 13 27 89 142 218 245 261 209 309 231 44 10 1,799 2.2b 
36 0.8% 140 13 17 64 127 228 265 223 242 302 281 47 1,818 2.1a 
37 1.0% 126 17 30 89 155 240 265 291 255 308 264 52 10 1,978 2.2 
38 0.7% 74 15 42 107 164 218 281 322 229 345 243 62 18 2,046 2.2a 
39 0.6% 206 16 27 82 145 243 289 268 276 331 316 57 12 2,063 2.1 
40 0.9% 97 22 46 111 172 246 291 322 284 333 273 75 17 2,193 2.2 
41 0.9% 20 22 49 115 180 249 330 381 254 386 277 79 22 2,345 2.2 
42 0.8% 105 20 45 108 168 246 319 330 317 357 304 71 18 2,302 2.2a 
43 1.0% 86 28 62 136 183 254 320 367 320 389 306 110 27 2,502 2.2a 
44 0.7% 50 33 70 153 207 279 365 430 368 415 323 145 41 2,829 2.3a 
45 0.9% 14 32 64 130 198 278 378 431 293 443 317 109 32 2,704 2.2a 
46 0.7% 13 42 78 151 219 312 453 493 334 478 354 147 48 3,108 2.3 
47 0.8% 10 54 92 165 242 357 541 517 316 497 389 174 64 3,409 2.3 
Cluster no.% coverageMean elevation (m)JanFebMarAprMayJunJulAugSepOctNovDecAnnualModality subgroup
3.7% 313 11 31 100 145 49 344 1.1 
4.3% 325 20 56 137 177 74 10 479 1.1a 
5.2% 343 31 74 159 205 95 13 584 1.1b 
5.0% 379 10 41 89 178 233 114 17 682 1.1a 
4.7% 420 14 55 105 200 259 127 20 781 1.1b 
5.0% 361 23 69 116 189 243 157 31 837 1.1 
3.7% 505 22 77 126 224 275 154 28 910 1.1 
2.6% 664 37 106 151 245 283 188 41 1,058 1.1a 
4.7% 396 32 87 134 196 254 188 44 943 1.1a 
10 4.4% 242 49 111 152 193 230 203 66 1,013 1.1b 
11 3.8% 465 15 66 150 178 264 293 230 65 1,267 1.1a 
12 1.9% 952 42 105 150 205 271 232 61 1,078 1.1a 
13 2.5% 569 26 78 154 186 256 314 289 118 1,438 1.1 
14 5.3% 242 18 69 135 170 191 213 234 101 1,142 1.2 
15 2.9% 161 33 88 162 175 190 211 249 138 12 1,269 1.2 
16 4.0% 373 11 55 127 168 212 261 267 96 1,204 1.2b 
17 3.3% 376 27 83 154 180 223 253 259 125 1,320 1.2a 
18 1.9% 814 11 61 126 181 208 256 260 287 183 24 1,603 1.2 
19 0.8% 89 13 31 92 134 177 221 174 86 159 162 48 14 1,312 2.3 
20 1.7% 239 23 80 124 153 167 144 104 189 160 35 12 1,199 2.2b 
21 1.7% 310 16 60 106 145 161 154 128 206 138 20 1,152 2.2 
22 1.6% 145 28 80 144 172 190 184 235 111 1,164 2.2a 
23 2.9% 299 11 46 96 152 171 167 156 231 133 15 1,192 2.2b 
24 2.1% 399 22 76 118 159 173 167 138 234 164 31 15 1,306 2.2 
25 0.9% 130 24 78 149 175 166 198 227 132 1,161 2.1a 
26 1.1% 145 35 95 176 194 181 209 251 161 14 1,327 2.1a 
27 1.4% 199 17 55 108 178 189 186 175 270 174 27 1,394 2.2 
28 1.1% 162 45 106 193 214 191 212 269 197 23 1,471 2.1 
29 1.1% 176 10 22 78 125 200 216 222 187 291 205 36 1,602 2.2 
30 0.9% 177 13 56 117 208 237 209 226 286 235 33 1,636 2.1 
31 1.0% 217 32 96 135 169 205 217 147 239 182 49 17 1,498 2.2a 
32 0.7% 94 10 38 105 146 189 245 275 179 280 200 57 20 1,745 2.2 
33 0.5% 19 21 38 93 129 227 329 222 77 162 174 63 18 1,553 2.3b 
34 0.7% 34 16 38 97 133 210 294 271 115 225 201 69 18 1,686 2.3 
35 0.9% 132 13 27 89 142 218 245 261 209 309 231 44 10 1,799 2.2b 
36 0.8% 140 13 17 64 127 228 265 223 242 302 281 47 1,818 2.1a 
37 1.0% 126 17 30 89 155 240 265 291 255 308 264 52 10 1,978 2.2 
38 0.7% 74 15 42 107 164 218 281 322 229 345 243 62 18 2,046 2.2a 
39 0.6% 206 16 27 82 145 243 289 268 276 331 316 57 12 2,063 2.1 
40 0.9% 97 22 46 111 172 246 291 322 284 333 273 75 17 2,193 2.2 
41 0.9% 20 22 49 115 180 249 330 381 254 386 277 79 22 2,345 2.2 
42 0.8% 105 20 45 108 168 246 319 330 317 357 304 71 18 2,302 2.2a 
43 1.0% 86 28 62 136 183 254 320 367 320 389 306 110 27 2,502 2.2a 
44 0.7% 50 33 70 153 207 279 365 430 368 415 323 145 41 2,829 2.3a 
45 0.9% 14 32 64 130 198 278 378 431 293 443 317 109 32 2,704 2.2a 
46 0.7% 13 42 78 151 219 312 453 493 334 478 354 147 48 3,108 2.3 
47 0.8% 10 54 92 165 242 357 541 517 316 497 389 174 64 3,409 2.3 

Clusters without asterisks do not have stations inside their boundaries.

aClusters with one rain-gauge.

bClusters with two rain-gauges.

Figure 4

(a) The spatial distribution of the 47 clusters of P and (b) levels of % confidence of classification reliability according to MLC.

Figure 4

(a) The spatial distribution of the 47 clusters of P and (b) levels of % confidence of classification reliability according to MLC.

Close modal

The clusters were further grouped based on P modality. The clusters 1 up to 18 present unimodal while the clusters 19 up to 47 present bimodal monthly variation of P. The clusters belonging in the bimodal group present higher precipitation in comparison to the unimodal group (Table 1) justifying the bimodal response of the maximum monthly P variation for the whole country (Figure 2(a)). The spatial threshold between the unimodal and bimodal responses of P is given in Figure 5. The % coverage of the unimodal and bimodal clusters was found equal to 69% and 31%, respectively (Table 1), justifying the dominance of the first group in the unimodal response of the mean monthly P variation for the whole country (Figure 2(a)). The general characteristics of the modality groups are the following:

  • The unimodal group is divided into two subgroups where the first (code 1.1) presents a unique P maximum extreme during August and the second one (code 1.2) during September (Table 1, Figure 5).

  • The bimodal group is divided into three subgroups where: (a) for the first subgroup (code 2.1), the minimum extreme, which is responsible for the bimodal rainy season, appears during July while September is the most rainy month; (b) for the second subgroup (code 2.2), the minimum extreme, which is responsible for the bimodal rainy season, appears during August while September is the most rainy month; and (c) for the third subgroup (code 2.3), the minimum extreme, which is responsible for the bimodal rainy season, appears during August while June or July is the most rainy month (Table 1, Figure 5).

  • The clusters 6, 9, 10, 11, 14, 16, and 42 present a distinct geographical dispersion in two different geographical regions (e.g., see black arrows for cluster 14 in Figure 4(a)). These clusters (except 42) are on the west and the east side of the upland regions of Bauchi, Kaduna, Kano, and Plateau states (located in the clusters 7, 8, 12, 13). This dispersion probably occurs because these upland regions force part of the incoming vapor masses from the south through mT to follow two different pathways through Niger and Benue River valleys (Figure S1(a)), respectively.

  • The mean monthly values of P for each cluster (Table 1) also provide the length of the rainy season and the total rainfall during the rainy season over the different areas of the country. The clusters 46 and 47 present the higher annual precipitation exceeding the threshold of 3,000 mm/yr while clusters 1 and 2 present the lower annual precipitation (below the threshold of 500 mm/yr). The clusters 1–13 and 16 present at least one month with no precipitation while the highest mean monthly precipitation appears in cluster 47 during June with a value of 541 mm/month.

Figure 5

The network of 33 rain-gauges and the groups and subgroups of precipitation patterns (unimodal, bimodal) based on the variation of mean monthly values for the period 1950–2000.

Figure 5

The network of 33 rain-gauges and the groups and subgroups of precipitation patterns (unimodal, bimodal) based on the variation of mean monthly values for the period 1950–2000.

Close modal

Evaluating the position of rain-gauges based on the cluster characteristics

The derived terrain segmentation of Figure 4(a) was used to evaluate the position of rain-gauges based on the cluster characteristics. Analyzing the position of the 33 existing stations (Figure 5), it was found that 22 out of 47 clusters do not contain a rain-gauge, while 17 and 8 clusters contain one and two stations, respectively (Table 1). The indication of reliability class <5% showed that 13 stations are located in regions of transitional precipitation patterns. In the case of clusters with more than one station, the reliability class can provide a good criterion for selecting the most representative one.

A coarser clustering, as is described by the zones of Figure 5, can also be used to analyze stations’ distribution. In this case, it was observed that the number of stations that are located in the modality subgroup zones with codes 1.1, 1.2, 2.1, 2.2, and 2.3 were 12, 3, 3, 12, and 3, respectively (Table 2). Their distribution is generally homogeneous in the three zones with codes 1.1, 2.1, and 2.2. In the case of 1.2 zone, the number of stations per unit area is the lowest and the stations are located close to the center of the zone and not at its west and east side where Niger and Benue River valleys are located. For the case of zone 2.3, it was observed that the southeast side contains only one station (Figure 5). In zone 2.3 are included the clusters 46 and 47 (Figure 4(a)) which are described by the highest annual precipitation (Table 1) without having a station within their boundaries.

Table 2

Main rain-gauges of Nigeria and their relative position in comparison to the terrain segmentation from cluster analysis (cluster number, reliability class within the cluster, and modality subgroup)

Station no.NameSt.elevation .Ha (m)Lat (dec.deg.)Lon (dec.deg.)Cluster numberReliability classModality subgroup
Abuza 344 9.2500 7.0000 17 75–95% 1.2 
Bauchi 609 10.2833 9.8167 25–50% 1.1 
Benin 79 6.3167 5.6000 38 <1% 2.2 
Bida 137 9.1000 6.0167 22 1–5% 2.2 
Calabar 63 4.9667 8.3500 44 <1% 2.3 
Enugu 137 6.4667 7.5500 35 5–25% 2.2 
Gusau 463 12.1667 6.7000 25–50% 1.1 
Ibadan 234 7.4333 3.9000 20 <1% 2.2 
Ibi 111 8.1833 9.7500 25 5–25% 2.1 
10 Ikom 93 5.9667 8.7167 42 75–95% 2.2 
11 Ilorin 308 8.4833 4.5833 23 1–5% 2.2 
12 Jos 1295 9.8667 8.9000 12 <1% 1.1 
13 Kaduna 645 10.6000 7.4500 16 25–50% 1.2 
14 Kano 476 12.0500 8.5333 1–5% 1.1 
15 Katsina 427 13.0167 7.6833 50–75% 1.1 
16 Lagos/Ikeja 40 6.5833 3.3333 33 <1% 2.3 
17 Lagos/Oshodi 19 6.5500 3.3500 33 <1% 2.3 
18 Lokoja 41 7.8000 6.7333 23 <1% 2.2 
19 Maiduguri 354 11.8500 13.0833 50–75% 1.1 
20 Makurdi 113 7.6833 8.6167 26 <1% 2.1 
21 Minna 262 9.6167 6.5333 16 25–50% 1.2 
22 Nguru 343 12.8833 10.4667 75–95% 1.1 
23 Ogoya 117 6.6667 8.8000 36 25–50% 2.1 
24 Ondo 287 7.1000 4.8333 31 5–25% 2.2 
25 Onitsha 86 6.1500 6.7833 35 5–25% 2.2 
26 Oshogbo 305 7.7833 4.4833 20 5–25% 2.2 
27 Port Harcourt 18 4.8500 7.0167 43 <1% 2.2 
28 Potiskum 488 11.7000 11.0333 5–25% 1.1 
29 Sokoto 351 13.0167 5.2500 5–25% 1.1 
30 Warri 5.5167 5.7333 45 <1% 2.2 
31 Yelwa 24 10.8833 4.7500 10 25–50% 1.1 
32 Yola 191 9.2333 12.4667 10 25–50% 1.1 
33 Zaria 664 11.1333 7.6833 11 5–25% 1.1 
Station no.NameSt.elevation .Ha (m)Lat (dec.deg.)Lon (dec.deg.)Cluster numberReliability classModality subgroup
Abuza 344 9.2500 7.0000 17 75–95% 1.2 
Bauchi 609 10.2833 9.8167 25–50% 1.1 
Benin 79 6.3167 5.6000 38 <1% 2.2 
Bida 137 9.1000 6.0167 22 1–5% 2.2 
Calabar 63 4.9667 8.3500 44 <1% 2.3 
Enugu 137 6.4667 7.5500 35 5–25% 2.2 
Gusau 463 12.1667 6.7000 25–50% 1.1 
Ibadan 234 7.4333 3.9000 20 <1% 2.2 
Ibi 111 8.1833 9.7500 25 5–25% 2.1 
10 Ikom 93 5.9667 8.7167 42 75–95% 2.2 
11 Ilorin 308 8.4833 4.5833 23 1–5% 2.2 
12 Jos 1295 9.8667 8.9000 12 <1% 1.1 
13 Kaduna 645 10.6000 7.4500 16 25–50% 1.2 
14 Kano 476 12.0500 8.5333 1–5% 1.1 
15 Katsina 427 13.0167 7.6833 50–75% 1.1 
16 Lagos/Ikeja 40 6.5833 3.3333 33 <1% 2.3 
17 Lagos/Oshodi 19 6.5500 3.3500 33 <1% 2.3 
18 Lokoja 41 7.8000 6.7333 23 <1% 2.2 
19 Maiduguri 354 11.8500 13.0833 50–75% 1.1 
20 Makurdi 113 7.6833 8.6167 26 <1% 2.1 
21 Minna 262 9.6167 6.5333 16 25–50% 1.2 
22 Nguru 343 12.8833 10.4667 75–95% 1.1 
23 Ogoya 117 6.6667 8.8000 36 25–50% 2.1 
24 Ondo 287 7.1000 4.8333 31 5–25% 2.2 
25 Onitsha 86 6.1500 6.7833 35 5–25% 2.2 
26 Oshogbo 305 7.7833 4.4833 20 5–25% 2.2 
27 Port Harcourt 18 4.8500 7.0167 43 <1% 2.2 
28 Potiskum 488 11.7000 11.0333 5–25% 1.1 
29 Sokoto 351 13.0167 5.2500 5–25% 1.1 
30 Warri 5.5167 5.7333 45 <1% 2.2 
31 Yelwa 24 10.8833 4.7500 10 25–50% 1.1 
32 Yola 191 9.2333 12.4667 10 25–50% 1.1 
33 Zaria 664 11.1333 7.6833 11 5–25% 1.1 

Limitations of the proposed methods and recommendations for optimum use

The application of the proposed methods using gridded climatic data for the evaluation of meteorological station networks presents the following limitations:

  • (a) The minimum period of temporal resolution of available gridded climatic data is usually a month. This attribute is a limiting factor, especially in the case of precipitation, since the proposed methods cannot capture differences in frequency, duration, and intensity of rainfall. The aforementioned limitations in the case of rainfall can lead to geographical dispersion of clusters (as in the case of clusters 6, 9, 10, 11, 14, 16, and 42). Clear dispersion of clusters (e.g., cluster 14, see black arrows in Figure 4(a)) can further be used to create sub-clusters. In the case of other climatic parameters (e.g., temperature, solar radiation, etc.), the proposed methodology is expected to be more robust.

  • (b) The results of the proposed methodology are highly dependent on the methods and data used for grid development. One of the most important problems in climatic grid development is the high uncertainty of grid values in high altitude regions or vast inaccessible regions where meteorological stations usually do not exist (Hijmans et al. 2005).

Due to the aforementioned limitations of the proposed methodology for evaluating the station network, it is recommended to be used as a supporting tool in combination with other methods such as those mentioned in the Introduction section. Some suggestions for future applications are:

  • (a) use of the proposed methodology as a primary analysis for setting the initial position and number of stations for the design of a new station network and then use of other methods (see Introduction) to optimize the network;

  • (b) use of the proposed methodology as an additional tool/criterion in order to evaluate the position of an existing network which has already been designed by other methods;

  • (c) use of the proposed methodology as a basic tool to investigate the changes on grid attributes and cluster formation due to changes in the position and number of stations. This case is very important since it can optimize not only the representativeness of the station network but also the quality of the grids, which is strongly dependent on the data of existing stations.

The application of the proposed statistical techniques captured the spatiotemporal variation of precipitation in Nigeria based on gridded precipitation data. Cross correlation revealed an intense change in monthly latitude, longitude, and elevation dependencies of precipitation. These changes were further investigated using PCA, which provided an excellent visualization of the northward movement effects of IDT on rainfall patterns, identifying also the period of its occurrence (June–September with a peak in August). The cluster analysis segmented the country of Nigeria into territories based on different precipitation patterns, which were further used to evaluate the representativeness of the existing network of rain-gauges to describe the spatiotemporal variation of precipitation. The cluster analysis revealed 47 clusters of which 22 do not have a rain-gauge, while eight clusters have more than one rain-gauge. Thus, more rain-gauges and a better distribution are required to describe the spatiotemporal variability of P in Nigeria. The proposed approach for the evaluation of rain-gauge networks (or meteorological stations in general) is mainly proposed as a complementary methodology and not as a complete method to substitute for other techniques since it is based on a coarse temporal resolution (mean monthly step), which is a limiting factor for capturing climatic peculiarities that may appear at smaller time steps. The combination of the specific methods can be implemented for all climatic variables (e.g., temperature, solar radiation, wind speed, etc.) and can support not only the evaluation of meteorological station networks but also other applications related to decision support systems for various management initiatives which are dependent on climate (e.g., water management in agriculture, flood risk assessment, etc.).

Adefolalu
,
D. O.
1988
Precipitation trends, evapotranspiration and the ecological zones of Nigeria
.
Theoretical and Applied Climatology
39
,
81
89
.
Adefolalu
,
D. O.
&
Oguntoyinbo
,
J. S.
1985
On rainfall distribution and agricultural planning in Nigeria
.
Malaysian Journal of Tropical Geography
11
,
1
11
.
Adhikary
,
S. K.
,
Yilmaz
,
A. G.
&
Muttil
,
N.
2015
Optimal design of rain gauge network in the Middle Yarra River catchment, Australia
.
Hydrological Processes
29
,
2582
2599
.
Alli
,
A. A.
,
Oguntunde
,
P. G.
,
Olufayo
,
A. A.
&
Fasinmirin
,
J. T.
2012
Implications of Trends and Cycles of Rainfall on Agriculture and Water Resource in the Tropical Climate of Nigeria
.
Hydrology for Disaster Management, Special Publication of the Nigerian Association of Hydrological Sciences
, pp.
188
200
.
Aschonitis
,
V. G.
,
Miliaresis
,
G.
,
Demertzi
,
K.
&
Papamichail
,
D.
2016
Terrain segmentation of Greece using the spatial and seasonal variation of reference crop evapotranspiration
.
Advances in Meteorology
,
Art.ID 3092671
, pp.
14
.
Ball
,
G. H.
&
Hall
,
D. J.
1965
A Novel Method of Data Analysis and Pattern Classification
.
Stanford Research Institute
,
Menlo Park, CA
,
USA
.
Barca
,
E.
,
Passarella
,
G.
&
Uricchio
,
V.
2008
Optimal extension of the rain gauge monitoring network of the Apulian Regional Consortium for Crop Protection
.
Environmental Monitoring and Assessment
145
,
375
386
.
Bello
,
N. J.
1998
A study of evidence of climate change based on rainfall seasonality and the reliability of rainfall regime in Nigeria
.
Proc. Sustained Africa
4
,
30
32
.
Boer
,
E. P. J.
,
de Beurs
,
K. M.
&
Hartkamp
,
A. D.
2001
Kriging and thin plate splines for mapping climate variables
.
International Journal of Applied Earth Observation and Geoinformation
3
,
146
154
.
Buttafuoco
,
G.
,
Caloiero
,
T.
&
Coscarelli
,
R.
2011
Spatial and temporal patterns of the mean annual precipitation at decadal time scale in southern Italy (Calabria region)
.
Theoretical and Applied Climatology
105
,
431
444
.
Caselles
,
V.
&
Melia
,
J.
1985
Use of Satellite Images in the Design of a Standard Meteorological Network
.
European Space Agency (Special Publication), ESA SP
, pp.
37
40
.
Chebbi
,
A.
,
Bargaoui
,
Z. K.
&
da Conceição Cunha
,
M.
2013
Development of a method of robust rain gauge network optimization based on intensity-duration-frequency
.
Hydrology and Earth System Sciences
17
,
4259
4268
.
Chen
,
Y. C.
,
Wei
,
C.
&
Yeh
,
H. C.
2008
Rainfall network design using kriging and entropy
.
Hydrological Processes
22
,
340
346
.
Chen
,
L.-J.
,
Chen
,
D.-L.
,
Wang
,
H.-J.
&
Yan
,
J.-H.
2009
Regionalization of precipitation regimes in China
.
Atmospheric and Oceanic Science Letters
2
,
301
307
.
Cheng
,
K. S.
,
Wei
,
C.
,
Cheng
,
Y. B.
,
Yeh
,
H. C.
&
Liou
,
J. J.
2008
Rain-gauge network evaluation and augmentation using geostatistics
.
Hydrological Processes
22
,
2555
2564
.
Demertzi
,
K.
,
Papamichail
,
D.
,
Aschonitis
,
V.
&
Miliaresis
,
G.
2014
Spatial and seasonal patterns of precipitation in Greece: the terrain segmentation approach
.
Global Nest Journal
16
,
988
997
.
Eastman
,
J. R.
&
Fulk
,
M.
1993
Long sequence time series evaluation using standardized principal components
.
Photogrammetric Engineering & Remote Sensing
59
,
1307
1312
.
Filippini
,
F.
,
Galliani
,
G.
&
Screpanti
,
F.
1994
Comparing optimization methods in the configuring of a network of meteorological stations
.
International Conference on Air Pollution – Proceedings
1
,
507
514
.
Gelman
,
A.
,
Carlin
,
J. B.
,
Stern
,
H. S.
,
Dunson
,
D. B.
,
Vehtari
,
A.
&
Rubin
,
D. B.
2014
Bayesian Data Analysis
, 3rd edn.
Chapman & Hall/CRC Press, Taylor and Francis Group
,
Boca Raton, FL
,
USA
.
Hassan
,
S. M.
,
Adefolalu
,
D. O.
&
Sani
,
M.
2009
Recent rainfall trends in the FCT, Abuja
.
Transactions of the Institute of Indian Geographers
31
,
49
56
.
Hijmans
,
R. J.
,
Cameron
,
S. E.
,
Parra
,
J. L.
,
Jones
,
P. G.
&
Jarvis
,
Α.
2005
Very high resolution interpolated climate surfaces for global land areas
.
International Journal of Climatology
25
,
1965
1978
.
Horel
,
J. D.
&
Dong
,
X.
2010
An evaluation of the distribution of remote automated weather stations (RAWS)
.
Journal of Applied Meteorology and Climatology
49
,
1563
1578
.
Kawachi
,
T.
,
Maruyama
,
T.
&
Singh
,
V. P.
2001
Rainfall entropy for delineation of water resources zones in Japan
.
Journal of Hydrology
246
,
36
44
.
Krstanovic
,
P. F.
&
Singh
,
V. P.
1992
Evaluation of rainfall networks using entropy. II: application
. Water Resources Management
6
,
295
314
.
Maaten
,
L.
&
Hinton
,
G.
2008
Visualizing high-dimensional data using t-SNE
.
Journal of Machine Learning Research
9
,
2579
2605
.
Marchetti
,
M.
,
Chapman
,
L.
,
Khalifa
,
A.
&
Buès
,
M.
2015
New role of thermal mapping in winter maintenance with principal components analysis
.
Advances in Meteorology
,
Art.ID 254795
, pp.
11
.
Mather
,
P. M.
2004
Computer Processing of Remotely-Sensed Images
, 3rd edn.
John Wiley and Sons
,
New York
,
USA
.
Miliaresis
,
G. C.
2013
Thermal anomaly mapping from night MODIS imagery of USA, a tool for environmental assessment
.
Environmental Monitoring and Assessment
185
,
1601
1612
.
Miliaresis
,
G. C.
2014
Spatiotemporal patterns of land surface temperature of Antarctica from MODIS monthly LST (MYD11C3) data
.
Journal of Spatial Science
59
,
157
166
.
Ninyerola
,
M.
,
Pons
,
X.
&
Roure
,
J. M.
2000
A methodological approach of climatological modelling of air temperature and precipitation through GIS techniques
.
International Journal of Climatology
20
,
1823
1841
.
Nour
,
M. H.
,
Smith
,
D. W.
&
Gamal El-Din
,
M.
2006
Geostatistical mapping of precipitation: implications for rain gauge network design
.
Water Science and Technology
53
(
10
),
101
110
.
Ogungbenro
,
S. B.
&
Morakinyo
,
T. E.
2014
Rainfall distribution and change detection across climatic zones in Nigeria
.
Weather and Climate Extremes
5
,
1
6
.
Oguntunde
,
P. G.
,
Abiodun
,
B. J.
&
Lischeid
,
G.
2011
Rainfall trends in Nigeria, 1901-2000
.
Journal of Hydrology
411
,
207
218
.
Oguntunde
,
P. G.
,
Lischeid
,
G.
,
Abiodun
,
B. J.
&
Dietrich
,
O.
2014
Analysis of spatial and temporal patterns in onset, cessation and length of growing season in Nigeria
.
Agricultural and Forest Meteorology
194
,
77
87
.
Olaniran
,
O. J.
1983
The monsoon factor and the seasonality of rainfall distribution in Nigeria
.
Malaysian Journal of Tropical Geography
7
,
38
45
.
Olaniran
,
O. J.
1988b
The July–August rainfall anomaly in Nigeria
.
Climatological Bulletin
22
,
26
38
.
Olaniran
,
O. J.
1990
Changing patterns of rain-days in Nigeria
.
GeoJournal
22
,
99
107
.
Olaniran
,
O. J.
1991b
Rainfall anomaly patterns in dry and wet years over Nigeria
.
International Journal of Climatology
11
,
177
204
.
Olaniran
,
O. J.
&
Sumner
,
G. N.
1989b
A study of climatic variability in Nigeria based on the onset, retreat, and length of the rainy season
.
International Journal of Climatology
9
,
253
269
.
Olaniran
,
O. J.
&
Sumner
,
G. N.
1990
Long-term variations of annual and growing season rainfalls in Nigeria
.
Theoretical and Applied Climatology
41
,
41
53
.
Olson
,
D. M.
,
Dinerstein
,
E.
,
Wikramanayake
,
E. D.
,
Burgess
,
N. D.
,
Powell
,
G. V. N.
,
Underwood
,
E. C.
,
D'Amico
,
J. A.
,
Itoua
,
I.
,
Strand
,
H. E.
,
Morrison
,
J. C.
,
Loucks
,
C. J.
,
Allnutt
,
T. F.
,
Ricketts
,
T. H.
,
Kura
,
Y.
,
Lamoreux
,
J. F.
,
Wettengel
,
W. W.
,
Hedao
,
P.
&
Kassem
,
K. R.
2001
Terrestrial ecoregions of the world: a new map of life on Earth
.
BioScience
51
,
933
938
.
Peel
,
M. C.
,
Finlayson
,
B. L.
&
McMahon
,
T. A.
2007
Updated world map of the Köppen–Geiger climate classification
.
Hydrology and Earth System Sciences
11
,
1633
1644
.
Richards
,
J. A.
1986
Remote Sensing Digital Image Analysis: An Introduction
.
Springer-Verlag
,
Berlin
,
Germany
.
Shaffrey
,
L. C.
,
Stevens
,
I.
,
Norton
,
W. A.
,
Roberts
,
M. J.
,
Vidale
,
P. L.
,
Harle
,
J. D.
,
Jrrar
,
A.
,
Stevens
,
D. P.
,
Woodage
,
M. J.
,
Demory
,
M. E.
,
Donners
,
J.
,
Clark
,
D. B.
,
Clayton
,
A.
,
Cole
,
J. W.
,
Wilson
,
S. S.
,
Connolley
,
W. M.
,
Davies
,
T. M.
,
Iwi
,
A. M.
,
Johns
,
T. C.
,
King
,
J. C.
,
New
,
A. L.
,
Slingo
,
J. M.
,
Slingo
,
A.
,
Steenman-Clark
,
L.
&
Martin
,
G. M.
2009
U.K. HiGEM: the new U.K. high-resolution global environment model – model description and basic evaluation
.
Journal of Climate
22
,
1861
1896
.
Shafiei
,
M.
,
Ghahraman
,
B.
,
Saghafian
,
B.
,
Pande
,
S.
,
Gharari
,
S.
&
Davary
,
K.
2014
Assessment of rain-gauge networks using a probabilistic GIS based approach
.
Hydrology Research
45
,
551
562
.
Siljestrom
,
P. A.
,
Moreno
,
A.
,
Vikgren
,
G.
&
Caceres
,
L. M.
1997
The application of selective principal components analysis to a Thematic Mapper image for the recognition of geomorphologic features configuration
.
International Journal of Remote Sensing
18
,
3843
3852
.
Tsintikidis
,
D.
,
Georgakakos
,
K. P.
,
Sperfslage
,
J. A.
,
Smith
,
D. E.
&
Carpenter
,
T. M.
2002
Precipitation uncertainty and raingauge network design within Folsom Lake watershed
.
Journal of Hydrologic Engineering
72
,
175
184
.
Wackernagel
,
H.
2003
Multivariate Geostatistics: An Introduction with Applications
.
Springer-Verlag
,
Berlin
,
Germany
.
Watanabe
,
M.
,
Suzuki
,
T.
,
O'Ishi
,
R.
,
Komuro
,
Y.
,
Watanabe
,
S.
,
Emori
,
S.
,
Takemura
,
T.
,
Chikira
,
M.
,
Ogura
,
T.
,
Sekiguchi
,
M.
,
Takata
,
K.
,
Yamazaki
,
D.
,
Yokohata
,
T.
,
Nozawa
,
T.
,
Hasumi
,
H.
,
Tatebe
,
H.
&
Kimoto
,
M.
2010
Improved climate simulation by MIROC5: mean states, variability, and climate sensitivity
.
Journal of Climate
23
,
6312
6335
.
Wilby
,
R. L.
&
Wigley
,
T. M. L.
1997
Downscaling general circulation model output: a review of methods and limitations
.
Progress in Physical Geography
21
,
530
548
.

Supplementary data