Performance Evaluation of Satellite-Based Rainfall Products over Nigeria

: Understanding the variability of rainfall is important for sustaining rain-dependent agriculture and driving the local economy of Nigeria. Paucity and inadequate rain gauge network across Nigeria has made satellite-based rainfall products (SRPs), which o ﬀ er a complete spatial and consistent temporal coverage, a better alternative. However, the accuracy of these products must be ascertained before use in water resource developments and planning. In this study, the performances of Climate Hazards Group Infrared Precipitation with Station data (CHIRPS), Precipitation estimation from Remotely Sensed Information using Artiﬁcial Neural Networks–Climate Data Record (PERSIANN-CDR), and Tropical Applications of Meteorology using SATellite data and ground-based observations (TAMSAT), were evaluated to investigate their ability to reproduce long term (1983–2013) observed rainfall characteristics derived from twenty-four (24) gauges in Nigeria. Results show that all products performed well in terms of capturing the observed annual cycle and spatial trends in all selected stations. Statistical evaluation of the SRPs performance show that CHIRPS agree more with observations in all climatic zones by reproducing the local rainfall characteristics. The performance of PERSIANN and TAMSAT, however, varies with season and across the climatic zones. Findings from this study highlight the beneﬁts of using SRPs to augment or ﬁll gaps in the distribution of local rainfall data, which is critical for water resources planning, agricultural development, and policy making. = Precipitation Estimation from Remotely Sensed Information using Artiﬁcial Neural Networks–Climate Data Record, TAMSAT = Tropical Applications of Meteorology using SATellite data and ground based observations.


Introduction
Over the last decades, occurrences of hydrologic extremes, such as flooding and droughts, have increased, as a result of human-induced climate change in West Africa. These impacts have led to tremendous socioeconomic losses in already vulnerable communities, and, most often, resulted in the deaths of human beings and livestock [1,2]. The economies of local agrarian communities have largely been affected because of their dependence on rain-fed agriculture [3,4]. Sustaining water resources development for improved agricultural production under varying climatic conditions and extreme climatic events has proved more challenging due to paucity of recorded climate data [5]. Consequently, sparse gauge networks do not allow for realistic temporal and spatial climate characteristics, as is evident in scientific literature from this region. Funding issues and lack of serious efforts have limited topography and multi-climate regions of Nigeria are still lacking in literature. Local evaluation of remotely sensed-data is necessary before utilization by government agencies because of their inherent uncertainties and limitations, especially in developing countries riddled with paucity of ground-truth data for adequate calibration and bias reductions [7,14]. Few studies [4,7,17] have attempted to validate the capability of some satellite-based rainfall products in Nigeria. However, these studies are region-specific with no study in literature showcasing satellite rainfall evaluation over the whole of Nigeria stretching over different climatic settings. In the northeastern part of Nigeria, the ability of satellite rainfall products to reproduce rainfall trends from 1981-2015 showed satisfactory results at decadal, monthly, and seasonal time scale [4]. In southwestern Nigeria, comparison of satellite and observed rainfall datasets from 1998 to 2016 also gave encouraging results [7]. However, owing to the inadequate state of climate gauge networks in Nigeria, assessing the performances of many of the existing satellite rainfall products for Nigerian conditions are lacking grossly in literature.
This study attempts to complement to the aforementioned few studies on satellite rainfall products evaluation in Nigeria by evaluating the performance of three products using twenty-four (24) ground-based stations located all over Nigeria. The selected products were based on rainfall datasets, which use a combination of MW and IR methods, and bias-corrected with rain gauge data. The objective of this study is to evaluate the capability of three satellite rainfall products (CHIRPS, PERSIANN-CDR, and TAMSAT) to reproduce local rainfall characteristic (seasonal and annual climatology) in Nigeria. We also aim to assess the utility of the gridded datasets in reproducing inter-annual variability of rainfall for the period 1983-2013. The Pearson coefficient of correlation (r), root mean square error (RMSE), and percent bias (PBIAS) were used to evaluate the performance of these satellite products over Nigeria. The remaining part of this study is structured as follows: Section 2 describes the data and methods, the findings are presented and discussed in Section 3, and the conclusions are made in Section 4.

Materials and Methods
Nigeria, the most populous country in West Africa has a land area of about 923,770 km 2 and is situated between Latitudes 4 • -14 • N and Longitudes 4 • -14 • E, as shown in Figure 1. The widely varying climatic pattern experienced in Nigeria is partly influenced by the presence of the Atlantic Ocean to the southern part, and the Sahara Desert to the northern part [20]. The climate pattern is also affected by distinct relief systems, such as lowlands, highlands, and plateau, as depicted in Figure 1. The major climatic zones divided latitudinally are Guinea (4 • -8 • N), Savannah (8 • -11 • N), and Sahel (11 • -14 • N) [21]. The spatial variability of Nigeria's climate is greatly influenced by the movement of the Inter-tropical Discontinuity (ITD), a narrow zone of trade-wind confluence between the southwest trade wind from the Atlantic Ocean and northeast trade wind from the Sahara Desert [22].
The Guinea region experiences a mean annual rainfall of about 1575-2533 mm, and the Savannah region is characterized by a mean annual rainfall of about 897-1535 mm while the Sahel region receives a mean annual rainfall of about 434-969 mm [23]. These authors noted that the Guinea and Savannah regions are characterized by a bimodal rainy season as a result of the abrupt non-linear latitudinal shift of rainfall band from a quasi-stationary position of 5 • N to about 10 • N. This process paves the way for the unimodal rainy season in Sahel during the period June-October, with the climatological peak in August. The Guinea and Savannah regions experience rainy seasons during March-May (MAM), June-August (JJA), and September-November (SON) seasons while the Sahel is majorly characterized by a peak rainfall during the July-August (JAS) season. However, all of the regions experience wide spread rainfall events during the period June-September (JJAS) as a result of more active convective activities that accompany the deep monsoon flow defined by the northward migration and surface position of the ITD.

Gauge Rainfall
This study used twenty-four (24) synoptic rainfall stations distributed across Nigeria as shown in Figure 1 and Table 1. Daily observed rainfall data for the period 1983-2013 for these synoptic stations, located across the three climatic regions, were obtained from the database of the Nigeria Meteorological Agency (NiMet) and used for evaluating the performance of the satellite rainfall products. However, the climatic stations are not evenly distributed within the three climatic zones. Five stations are found within the Sahel, eight stations within the savannah, while eleven stations are located within the Guinea climatic zone. The altitudes of the stations vary from about 40 m to 1160 m above sea level [23].

Gauge Rainfall
This study used twenty-four (24) synoptic rainfall stations distributed across Nigeria as shown in Figure 1 and Table 1. Daily observed rainfall data for the period 1983-2013 for these synoptic stations, located across the three climatic regions, were obtained from the database of the Nigeria Meteorological Agency (NiMet) and used for evaluating the performance of the satellite rainfall products. However, the climatic stations are not evenly distributed within the three climatic zones. Five stations are found within the Sahel, eight stations within the savannah, while eleven stations are located within the Guinea climatic zone. The altitudes of the stations vary from about 40 m to 1160 m above sea level [23]. The advent of satellite observing platforms in the late twentieth century has brought lots of benefits by providing more spatially complete datasets that helps to advance our knowledge in atmospheric science and other environmental-related disciplines. Improved spatiotemporal resolution of satellite rainfall products is useful as input in hydrologic modeling, especially in developing nations and data-scarce regions such as Nigeria. Table 2 shows an overview of the three gridded satellite rainfall products, which are freely available on the internet, and were assessed against gauged rainfall data in this study. These three (3) products under consideration in this study showed satisfactory results out of a total of ten (10) products, which were evaluated in a large-scale regional study in West Africa [24]. Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) data [10] was developed by scientists from the Climate Hazard Group at the University of California, Sant Barbara, in conjunction with the United States Geological Survey (USGS), specifically for monitoring droughts or to analyze shifts in rainfall in the data sparse African continent. Information on the data inputs used in developing CHIRPS is extensively reported on in literature [10].
Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record (PERSIANN-CDR; hereafter PERSIANN) was developed by the Center for Hydrometeorology and Remote Sensing group of the University of California, Irvine, in conjunction with National Oceanic and Atmospheric Administration (NOAA) [1].
Tropical Applications of Meteorology using SATellite data and ground based observations (TAMSAT) was developed at the University of Reading, and integrates about four thousand (4000) stations across Africa [24]. Information on its development is well documented in several studies [12,13].

Quality Control
Quality control and homogeneity tests were performed on all observed datasets using RClimDex and RHtestssoftware packages in a previous study by Gbode et al. [23]. Quality control was performed in other to remove erroneous values in the data series, such as negative rainfall values or days with daily rainfall values greater than 200 mm, and replaced with −999. A homogeneity test is important in order to correct anomalies that might result from a change in station location and faulty gauging equipment [25,26]. Quality control checks were also performed on all of the satellite rainfall datasets using RClimDex software. The permission to use this software package was obtained upon request from the model developers (http//:www.etccdi.pacificclimate.org). Missing and/or erroneous values within these data sets were replaced with −99 before use in this study.

Methodology
Daily rainfall data (1983-2013) were extracted from three satellite rainfall products using geographic coordinates of the twenty-four (24) synoptic stations (see Table 1) and processed for comparison with observed data. The satellite data were all extracted at grid points closest to the location of each synoptic station. The rainfall amounts were characterized into mean seasonal, annual, and inter-annual variations in order to evaluate how well the selected remotely sensed rainfall datasets reproduced observed rainfall characteristics from 1983 to 2013. Daily values were aggregated to monthly values for all datasets under study. Aggregated rainfall values for March-May (early rainfall season), June-August (mid rainfall season), and September-November (late rainfall season) were processed, and their mean determined to represent mean seasonal rainfall for the different seasons. Consequently, rainfall amounts for June-September were also aggregated and the mean for this period obtained to represent the JJAS season. The JJAS season is the only period in Nigeria when all climatic zones receive rainfall and it is highly influenced by the Intertropical Convergence Zone (ITCZ) oscillations.
Climatological and statistical evaluations of CHIRPS, PERSIANN, and TAMSAT were carried out against observed datasets. Spatial rainfall pattern representing seasonal climatology (MAM, JJA, and SON) of satellite products were assessed against observed datasets. Annual cycle and inter-annual variability for all locations ( Figure 1) were also evaluated against observed (gauged) data. The ability of the satellite rainfall products to reproduce observed extreme events (wet and dry spells) at the selected twenty-four climate stations were evaluated using the standardized precipitation index (SPI) [27].The SPI is based on the probability of rainfall for a given time period and applied in hydro-meteorological studies for monitoring of drought conditions. The SPI is calculated by fitting a rainfall time series to a probability distribution, which is then normalized so that the mean SPI for that location is zero [27]. These authors further stated that positive SPI values implies greater than median rainfall while negative values signifies lower than median rainfall.
Comparison of monthly rainfall cumulative distribution frequency (CDF) for all satellite product datasets against gauged datasets were performed at all stations to evaluate their deviations from observed/gauged patterns. Furthermore, trend analysis was performed using the Mann-Kendall (MK) test statistic [28] to detect monotonic changes in the rainfall time series at 5% significance level. This statistic tests whether to accept that there is a monotonic trend (alternative hypothesis, Ha) or to reject the null hypothesis (Ho), which states that no trend is present in the time series data [29]. The MK statistic has been used extensively in hydro-meteorological studies [29][30][31] for analyzing monotonic changes in time series data. The MK statistic is stated as [30]: where, x j and x i are sequential data values; n is length of datasets. The Signum (Sgn) function is given as: The statistics (S), the mean E(S), and the variance V(S) can be estimated as [29]: where, t i is the extent of any given tie. The standardized test statistic, Z is stated as: The remotely sensed rainfall datasets were validated for mean seasonal (for JJAS regime) and annual rainfall using the Pearson correlation coefficient (r), root mean square error (RMSE), and percent bias (PBIAS), and were presented by showing their spatial patterns. The equations representing these statistical models are shown below, as reported in other studies [1,24,[32][33][34].
where, O and M are rain gauge and model values, respectively, O and M are mean rain gauge and model values respectively, n = number of data pairs, α = measure of flow variability error, β = bias term. The spatial distribution of climate variables is important for understanding hydrological processes [35]. The inverse weighting distance (IDW) interpolation method was used in presenting spatial distribution of rainfall characteristics in this study. This method was adopted due to its popularity and wide applications in hydrology [36,37]. The IDW interpolation technique assumes that the weight between an observed point and unobserved point decreases exponentially as their distance increases. The characteristics of the interpolated cells were controlled by applying the variable search radius and adopting the default number of input points in the ArcMap window. The IDW function in ArcMap version 10.3 was used in this study and is implemented as follows [36]: where, Y(X o ) = interpolated value at point X o ; Y(X i ) = observed value at point X i ; n = number of observations; λ = weight; P = power; d io = distance between unknown point and known point.

Seasonal Climatology
The mean seasonal climatology results for March to May (MAM), June to August (JJA), and from September to November (SON), as well as the June to September (JJAS) period, during which the monsoon and associated rainfall widely dominates the West African region, are presented in this section. Evaluation of the JJAS period is important because about 60% of the West African population depends on rain-fed agriculture for their source of livelihood [18]. These classifications were computed using datasets from 1983 to 2013 for gauge data in comparison to CHIRPS, PERSIANN, and TAMSAT datasets for locations, as shown in Figure 1. These comparisons were performed in order to assess the difference between observed and remotely-sensed rainfall datasets.
The spatial distribution of mean MAM, JJA, and SON seasons for observed (gauged), CHIRPS, PERSIANN, and TAMSAT data are presented in Figure 2a-c. As seen in Figure 2a (MAM season), all of the products captured the high rainfall in the south and the low rainfall in the northeast of Nigeria. The lowest gauged seasonal rainfall for this period was about 5 mm at Nguru, while the highest amount of 232 mm was in Calabar, which is expected, as the latter station is located in the Guinea coast closer to the Atlantic Ocean, while the former is situated in the Sahel. Though the products performed reasonably well in capturing the seasonal amount of rainfall within the Savannah and Sahel regions, they poorly reproduced observed amounts in the Guinea coast region. All products overestimated MAM seasonal amounts by 1-23 mm for Guinea coast (Ikom), Savannah (Bida and Kaduna), and Sahel (Nguru) while underestimations occurred within the range of 2-52 mm in Guinea (Benin, Calabar, Enugu, Ibadan, Iseyin), Savannah (Bauchi), and Sahel (Kano, Katsina, Lokoja).

Seasonal Climatology
The mean seasonal climatology results for March to May (MAM), June to August (JJA), and from September to November (SON), as well as the June to September (JJAS) period, during which the monsoon and associated rainfall widely dominates the West African region, are presented in this section. Evaluation of the JJAS period is important because about 60% of the West African population depends on rain-fed agriculture for their source of livelihood [18]. These classifications were computed using datasets from 1983 to 2013 for gauge data in comparison to CHIRPS, PERSIANN, and TAMSAT datasets for locations, as shown in Figure 1. These comparisons were performed in order to assess the difference between observed and remotely-sensed rainfall datasets.
The spatial distribution of mean MAM, JJA, and SON seasons for observed (gauged), CHIRPS, PERSIANN, and TAMSAT data are presented in Figure 2a-c. As seen in Figure 2a (MAM season), all of the products captured the high rainfall in the south and the low rainfall in the northeast of Nigeria. The lowest gauged seasonal rainfall for this period was about 5 mm at Nguru, while the highest amount of 232 mm was in Calabar, which is expected, as the latter station is located in the Guinea coast closer to the Atlantic Ocean, while the former is situated in the Sahel. Though the products performed reasonably well in capturing the seasonal amount of rainfall within the Savannah and Sahel regions, they poorly reproduced observed amounts in the Guinea coast region. All products  During the JJA (Figure 2b), the seasonal spatial pattern of rainfall amount for all products was consistent with observed datasets. All of the products recorded low seasonal rainfall amounts in Nguru and high amount in Calabar. Specifically, all of the models showed underestimations in Guinea (Asaba, Awka, Benin, Ibadan, and Ijebu), Savannah (Bauchi and Bida), and Sahel (Gombe and Kano). The CHIRPS product agrees more with the observed data than PERSIANN and TAMSAT datasets in terms of capturing the spatial distribution pattern during the JJA season.
During the SON season (Figure 2c), all of the products were also able to reproduce the spatial and seasonal pattern of rainfall over Nigeria with the highest rainfall recorded in the south and lowest in the north. The observed seasonal amount was overestimated in many of the locations across the three climatic zones.
Generally, the spatial and seasonal patterns of all the rainfall products were consistent with the gauge (observed) dataset for all considered seasons. This revealed that all of the products captured the seasonal south-north rainfall oscillations, which are tightly coupled with the ITCZ latitudinal migration.

Annual Cycle of Mean Monthly Rainfall
The annual cycles (1983-2013) of mean monthly rainfall at twenty-four (24) point-based synoptic station scale were compared to corresponding point-based datasets for the satellite-based products in the Sahel, Savannah, and Guinea coast climate zones as shown in Figure 3a-c, respectively. The magnitude of the errors of the satellite-based products from observed data for all locations are presented in Table 3. Results showed that satellite product datasets were able to capture trends and peaks at all synoptic locations, with error deviations ranging from 3 to 60 mm. This shows that the latitudinal oscillations of the ITCZ from southern latitudes to northern latitudes, which results in convective processes, were well captured by all satellite-based products.
In the Sahel region (e.g., Nguru station, Figure 3a), the products underestimated observed rainfall peak at Gusau and Kano synoptic stations, with CHIRPS producing the lowest RMSE (Table 3) values of 13.5 mm and 35.6 mm, respectively. All gridded products produced good model fits with minimal residuals ranging from 3 mm to 11 mm for Maiduguri and Nguru synoptic stations. At Katsina, rainfall peak was over predicted by all products with RMSE values of 12 mm, 15 mm, and 17 mm for TAMSAT, CHIRPS, and PERSIANN, respectively.
The seasonal cycle was reasonably captured by all the entire products in the Savanna zone. Th mean seasonal cycle at Jos climate station is depicted in Figure 3b. CHIRPS and TAMSAT datasets represent the peaks better than PERSIANN, especially in Ilorin and Minna. Overestimations of rainfall peaks by all satellite products were observed at Bida, Ibi, Ilorin, Kaduna, and Jos with RMSE values ranging from 7-26 mm. Generally, CHIRPS and TAMSAT presented lower RMSE values ranging from 7.2-20.8 mm in comparison to higher values by the PERSIANN, which ranges from 12.3-23.2 mm.
Contrary to results from Sahel and Savannah regions, larger biases were detected for all of the products in all the locations within the Guinea region, except at Iseyin (Figure 3c) and Lokoja, as shown in Table 3. This notable difference can be partly explained by the heterogeneous land use land cover of the region, and the presence of the Atlantic Ocean at the south border of the country, which greatly influence rainfall processes. All of the satellite rainfall datasets showed a wide range of error magnitudes from 9.2 to 60.1 mm. High values of RMSE > 10 mm were recorded by all products at all synoptic stations, except for minimal values of RMSE < 10 mm recorded by CHIRPS at Lokoja. All satellite products overestimated (underestimated) rainfall peak at Ikom (Benin), showing large RMSE value > 29 mm.
It is important to also note that the entire satellite-based products showed strong agreement (r > 0.9) with gauged datasets at all locations for mean monthly data in Figure 3. Generally, the CHIRPS and TAMSAT datasets performed better in capturing the unimodal and bimodal annual rainfall pattern as influenced by the ITCZ for all locations studied. Unimodal and bimodal rainfall patterns as a result of oscillations of the ITCZ from 15 • S to 15 • N [34] were well represented by all satellite products in all climatic zones of Nigeria. In a previous study [7], the performance of CHIRPS rainfall estimates in comparison with gauge records for six (6) stations in southwest Nigeria showed strong relationship with high correlation (r) values greater than 0.70. The CHIRPS also showed strong relationships with gauge data at monthly and seasonal time resolutions in a study [4] conducted in the Sudano-Sahelian zone of Nigeria. In another study [5] in East Africa (Ethiopia, Kenya, and Tanzania), CHIRPS was reported as the preferential data source for climate change and hydrological studies in ungauged locations. These authors noted that CHIRPS was more accurate than African Rainfall Climatology, version 2.0, (ARC), Observational-Reanalysis Hybrid (ORH), and Regional Climate Models (RCMs) in reproducing mean monthly rainfall amount over East Africa. Dinku et al., in their study [8] over East Africa (Ethiopia, Kenya, Somalia, Uganda, Rwanda, and Tanzania), reported strong relationship with CHIRPS, TAMSAT, and gauge data with correlation values greater than 0.9 at monthly time steps. estimates in comparison with gauge records for six (6) stations in southwest Nigeria showed strong relationship with high correlation (r) values greater than 0.70. The CHIRPS also showed strong relationships with gauge data at monthly and seasonal time resolutions in a study [4] conducted in the Sudano-Sahelian zone of Nigeria. In another study [5] in East Africa (Ethiopia, Kenya, and Tanzania), CHIRPS was reported as the preferential data source for climate change and hydrological studies in ungauged locations. These authors noted that CHIRPS was more accurate than African Rainfall Climatology, version 2.0, (ARC), Observational-Reanalysis Hybrid (ORH), and Regional Climate Models (RCMs) in reproducing mean monthly rainfall amount over East Africa. Dinku et al., in a study [8] in East Africa (Ethiopia, Kenya, Somalia, Uganda, Rwanda, and Tanzania), reported strong relationship with CHIRPS, TAMSAT, and gauge data with correlation values greater than 0.9 at monthly time steps.

Inter-Annual Rainfall Anomaly
The capabilities of CHIRPS, PERSIANN, and TAMSAT products to reproduce observed year-toyear rainfall anomalies using standardized precipitation index (SPI) within the three climatic zones are shown in Figure 4a-c by using a single station in each zone for illustrations. The SPI at annual resolutions were estimated differently using rainfall records for gauge and satellite-based rainfall

Inter-Annual Rainfall Anomaly
The capabilities of CHIRPS, PERSIANN, and TAMSAT products to reproduce observed year-to-year rainfall anomalies using standardized precipitation index (SPI) within the three climatic zones are shown in Figure 4a-c by using a single station in each zone for illustrations. The SPI at annual resolutions were estimated differently using rainfall records for gauge and satellite-based rainfall products at point-based location scale for the period 1983-2013. Table 4 presents the correlation scores between each satellite product and gauge data for all locations within the study area at the annual time step. All products exhibited moderate agreement with gauged datasets with correlation values greater than 0.5 (r > 0.5), except at locations within the Sahel region (Katsina; TAMSAT in Gusau), Savannah region (Bida, Ilorin, Jos and Kaduna; PERSIANN at Bauchi; PERSIANN and TAMSAT at Ibi and Gombe), and the Guinea coast region (Calabar and Ikom; CHIRPS at Benin; PERSIANN at Asaba and Ibadan), as shown in Table 4. SPI results obtained for Kastina, Bida, and Ikom stations are shown in Figure 4a, 4b, and 4c, respectively. For the thirty-one (31) years of rainfall series considered, and in the majority of the locations, the CHIRPS dataset exhibited satisfactory performance over PERSIANN and TAMSAT in reproducing the year-year variations of rainfall anomalies. Unsatisfactory correlation found in most of the locations may be attributed to the presence of large-scale forcings on local climatic conditions [31].  Table 4 presents the correlation scores between each satellite product and gauge data for all locations within the study area at the annual time step. All products exhibited moderate agreement with gauged datasets with correlation values greater than 0.5 (r > 0.5), except at locations within the Sahel region (Katsina; TAMSAT in Gusau), Savannah region (Bida, Ilorin, Jos and Kaduna; PERSIANN at Bauchi; PERSIANN and TAMSAT at Ibi and Gombe), and the Guinea coast region (Calabar and Ikom; CHIRPS at Benin; PERSIANN at Asaba and Ibadan), as shown in Table 4. SPI results obtained for Kastina, Bida, and Ikom stations are shown in Figure 4a, 4b, and 4c, respectively. For the thirty-one (31) years of rainfall series considered, and in the majority of the locations, the CHIRPS dataset exhibited satisfactory performance over PERSIANN and TAMSAT in reproducing the year-year variations of rainfall anomalies. Unsatisfactory correlation found in most of the locations may be attributed to the presence of largescale forcings on local climatic conditions [31].

Empirical Cumulative Distribution Frequency
The ability of satellite-based rainfall products to reproduce the frequency of gauged monthly rainfall amounts from 1983 to 2013 is evaluated using the empirical cumulative distribution function (ECDF) plots. Comparisons of the ECDFs of mean monthly rainfall at point-based scale for gauged and satellite products within Sahel, Savannah, and Guinea Coast climatic zones are shown in Figure

Empirical Cumulative Distribution Frequency
The ability of satellite-based rainfall products to reproduce the frequency of gauged monthly rainfall amounts from 1983 to 2013 is evaluated using the empirical cumulative distribution function (ECDF) plots. Comparisons of the ECDFs of mean monthly rainfall at point-based scale for gauged and satellite products within Sahel, Savannah, and Guinea Coast climatic zones are shown in Figure 5a-c using one station in each climatic zone. Generally, similar patterns of gauged monthly rainfall distributions in all locations were captured by the satellite products. The frequency of occurrence of monthly rainfall within the Sahel climatic zone was significantly overestimated by all products at the Kano synoptic station ( Figure 5a) in the range of 200-400 mm/month. TAMSAT and PERSIANN exhibited overestimation (100-200 mm/month) at Gusau and Katsina, respectively. However, all products showed close monthly probabilities of rainfall to gauged datasets at Maiduguri and Nguru.
In the Savannah region, both CHIRPS and TAMSAT captured overall frequencies of observed values at all locations, although with slight margins, as shown in Figure 5b. The PERSIANN showed consistent frequencies with gauged datasets, but underestimated mean monthly rainfall frequency at Ibi, Ilorin, Jos, Kaduna, and Minna in the range of 150-300 mm/month.
In the Guinea coast climatic zone, the entire product showed consistent frequencies, except at Ikom (Figure 5c), where CHIRPS and TAMSAT significantly underestimated rainfall frequency in the range of 200-400 mm/month. Conversely, all products overestimated rainfall at more than 200 mm/month (at Asaba, Benin, Awka, Ibadan, Ijebu, and Iseyin), and more than 400 mm/month (at Calabar). 5a-c using one station in each climatic zone. Generally, similar patterns of gauged monthly rainfall distributions in all locations were captured by the satellite products. The frequency of occurrence of monthly rainfall within the Sahel climatic zone was significantly overestimated by all products at the Kano synoptic station (Figure 5a) in the range of 200-400 mm/month. TAMSAT and PERSIANN exhibited overestimation (100-200 mm/month) at Gusau and Katsina, respectively. However, all products showed close monthly probabilities of rainfall to gauged datasets at Maiduguri and Nguru.
In the Savannah region, both CHIRPS and TAMSAT captured overall frequencies of observed values at all locations, although with slight margins, as shown in Figure 5b. The PERSIANN showed consistent frequencies with gauged datasets, but underestimated mean monthly rainfall frequency at Ibi, Ilorin, Jos, Kaduna, and Minna in the range of 150-300 mm/month.
In the Guinea coast climatic zone, the entire product showed consistent frequencies, except at Ikom (Figure 5c), where CHIRPS and TAMSAT significantly underestimated rainfall frequency in the range of 200-400 mm/month. Conversely, all products overestimated rainfall at more than 200 mm/month (at Asaba, Benin, Awka, Ibadan, Ijebu, and Iseyin), and more than 400 mm/month (at Calabar).   Table 5 shows the p-values result of the non-parametric Mann-Kendall method used to assess the presence of significant trends at annual and seasonal (MAM, JJAS, and SON) time scales. The test was performed for all datasets at all synoptic stations under consideration in this study. In the Sahel, the entire satellite products and gauged datasets exhibited significant positive trend in the annual rainfall series, except at Gusau, where only TAMSAT and PERSIANN showed similar consistency. MAM rainfall features a positive significant trend in the local observed rainfall in Katsina, which is accurately replicated by PERSIANN. This significant positive trend was also observed at Gusau (PERSIANN and TAMSAT), Maiduguri (CHIRPS), and Nguru (CHIRPS and PERSIANN). The JJAS observed rainfall showed an increasing significant trend, which is in-line with all satellite products dataset in this zone, except for Gusau (gauge, CHIRPS, and TAMSAT were not significant). During the SON season, observed data from Kano and Maiduguri exhibited a significant trend that was replicated by PERSIANN and TAMSAT, as seen in Table 5.

Trend Analysis
In the Savannah, all satellite-based products and local rainfall annual series showed a significant positive trend at Bauchi. TAMSAT exhibited a significant increasing trend at all locations within this zone. The same was observed for PERSIANN at Bida, Jos, Kaduna, and Minna, while CHIRPS showed the same trend in Jos. At the seasonal scale, both CHIRPS and TAMSAT indicated a significant increasing trend, the same as seen in the gauged rainfall at Bauchi during the JJAS period. During the SON period, significant positive trends exhibited by the gauge time series at Ibi and Ilorin were reproduced by all satellite products, while in Minna, only PERSIANN and TAMSAT were able to replicate the gauged increasing significant trend successfully. Many of the products replicated seasonal observed trends in most locations, but with different levels of accuracy as presented in Table  5.
Over the Guinea coast, TAMSAT was consistent in exhibiting significant positive trends in all locations at the annual time step (Table 5), and was consistent with gauged data recorded at Benin, Calabar, Enugu, and Ijebu. From Table 5, significant positive trends shown in the gauge time series were also replicated by PERSIANN (at Ijebu), and CHIRPS (at Calabar). At the seasonal scale during the MAM period, (Table 5), TAMSAT exhibited positive significant trends in all stations within the Guinea coast, except in Iseyin. This significant trend was only exhibited in the gauged time series at Benin. For the high monsoon period (JJAS season), gauged data at Calabar exhibited same positive significant trends with CHIRPS and TAMSAT. Moreover, observed data at Ibadan, Ijebu, and Iseyin showed the same significant positive trends with all satellite products during the SON season. Similar  Table 5 shows the p-values result of the non-parametric Mann-Kendall method used to assess the presence of significant trends at annual and seasonal (MAM, JJAS, and SON) time scales. The test was performed for all datasets at all synoptic stations under consideration in this study. In the Sahel, the entire satellite products and gauged datasets exhibited significant positive trend in the annual rainfall series, except at Gusau, where only TAMSAT and PERSIANN showed similar consistency. MAM rainfall features a positive significant trend in the local observed rainfall in Katsina, which is accurately replicated by PERSIANN. This significant positive trend was also observed at Gusau (PERSIANN and TAMSAT), Maiduguri (CHIRPS), and Nguru (CHIRPS and PERSIANN). The JJAS observed rainfall showed an increasing significant trend, which is in-line with all satellite products dataset in this zone, except for Gusau (gauge, CHIRPS, and TAMSAT were not significant). During the SON season, observed data from Kano and Maiduguri exhibited a significant trend that was replicated by PERSIANN and TAMSAT, as seen in Table 5.

Trend Analysis
In the Savannah, all satellite-based products and local rainfall annual series showed a significant positive trend at Bauchi. TAMSAT exhibited a significant increasing trend at all locations within this zone. The same was observed for PERSIANN at Bida, Jos, Kaduna, and Minna, while CHIRPS showed the same trend in Jos. At the seasonal scale, both CHIRPS and TAMSAT indicated a significant increasing trend, the same as seen in the gauged rainfall at Bauchi during the JJAS period. During the SON period, significant positive trends exhibited by the gauge time series at Ibi and Ilorin were reproduced by all satellite products, while in Minna, only PERSIANN and TAMSAT were able to replicate the gauged increasing significant trend successfully. Many of the products replicated seasonal observed trends in most locations, but with different levels of accuracy as presented in Table 5.
Over the Guinea coast, TAMSAT was consistent in exhibiting significant positive trends in all locations at the annual time step (Table 5), and was consistent with gauged data recorded at Benin, Calabar, Enugu, and Ijebu. From Table 5, significant positive trends shown in the gauge time series were also replicated by PERSIANN (at Ijebu), and CHIRPS (at Calabar). At the seasonal scale during the MAM period, (Table 5), TAMSAT exhibited positive significant trends in all stations within the Guinea coast, except in Iseyin. This significant trend was only exhibited in the gauged time series at Benin. For the high monsoon period (JJAS season), gauged data at Calabar exhibited same positive significant trends with CHIRPS and TAMSAT. Moreover, observed data at Ibadan, Ijebu, and Iseyin showed the same significant positive trends with all satellite products during the SON season. Similar significant positive trends were also exhibited between gauge, PERSIANN, and TAMSAT datasets in Enugu and Ikeja, and between gauge data and TAMSAT at Benin in this season.

Evaluation of Satellite Rainfall Products
The capability of the satellite products in reproducing local rainfall characteristics at synoptic station-level was assessed using statistical methods. Results of these statistics at mean seasonal and annual time resolution for the period 1985-2013 are presented as spatial plots produced using the IDW interpolation method. Suitability of satellite products is dependent, not only on low RMSE and PBIAS values, but also on strong correlation coefficient between gauge and model values [31].

Inter-Annual Variation of Mean Seasonal Rainfall
The JJAS period was selected for evaluating the performance of the satellite products because, during this period, all three climatic zones in Nigeria experience rainfall [38]. Each satellite product was compared against gauge datasets using statistics: r, RMSE, and PBIAS, and are presented in Figures 6-8, respectively. Within the Sahel (Figure 6), CHIRPS and PERSIANN performed satisfactory exhibiting good correlations (r > 0.5) against gauge data, except in Kastina (where r < 0.4). However, TAMSAT presented weak correlations (in Kastina and Gusau), with an r value of -0.01 at Gusau. Results also indicate statistical significance at 95% confidence level for all products, except in Kastina. The spatial plot of the satellite product bias (Figure 7) exhibits low bias in the range -33.6 to 23.1% and RMSE (Figure 8) of 86-104 mm in Kano.
Results within the Savannah zone ( Figure 6) show that the gridded products generally were not suitable in simulating seasonal JJAS rainfall within this zone. Only the CHIRPS product reliably showed strong correlation (r > 0.5) in Bauchi, Gusau, Ilorin, and Ibi, but was not suitable (r < 0.5) at Bida, Kaduna, and Minna. PERSIAN and TAMSAT exhibited weak correlations (r < 0.4) at all locations, except that TAMSAT performed well at Bauchi (r = 0.6). However, low bias ( Figure 7) and RMSE ( Figure 8) were mostly prevalent over the Savanna zone.
In the Guinea coast zone (Figure 6), CHIRPS was more suitable and shows reasonable correlation with gauge datasets except at Ikom, where it exhibited weak correlation (r = 0.07) in comparison to PERSIANN(r = 0.3) and TAMSAT (0.1). Consequently, the correlation between CHIRPS and gauge were found to be significant at 95% confidence level, except at Ikom and Calabar. Furthermore, CHIRPS and TAMSAT exhibit bias of 26.3% and 20.3%, respectively. The RMSE values of above 78 mm were observed in Benin, Calabar, and Ikom were consistent for all products. Generally, the CHIRPS products showed good results during the JJAS rainfall season as compared to other products under study.
The capability of the satellite products in reproducing local rainfall characteristics at synoptic station-level was assessed using statistical methods. Results of these statistics at mean seasonal and annual time resolution for the period 1985-2013 are presented as spatial plots produced using the IDW interpolation method. Suitability of satellite products is dependent, not only on low RMSE and PBIAS values, but also on strong correlation coefficient between gauge and model values [31].

Inter-Annual Variation of Mean Seasonal Rainfall
The JJAS period was selected for evaluating the performance of the satellite products because, during this period, all three climatic zones in Nigeria experience rainfall [38]. Each satellite product was compared against gauge datasets using statistics: r, RMSE, and PBIAS, and are presented in Figures 6, 7, and 8, respectively. Within the Sahel (Figure 6), CHIRPS and PERSIANN performed satisfactory exhibiting good correlations (r > 0.5) against gauge data, except in Kastina (where r < 0.4). However, TAMSAT presented weak correlations (in Kastina and Gusau), with an r value of -0.01 at Gusau. Results also indicate statistical significance at 95% confidence level for all products, except in Kastina. The spatial plot of the satellite product bias (Figure 7) exhibits low bias in the range -33.6 to 23.1% and RMSE (Figure 8) of 86-104 mm in Kano.
Results within the Savannah zone ( Figure 6) show that the gridded products generally were not suitable in simulating seasonal JJAS rainfall within this zone. Only the CHIRPS product reliably showed strong correlation (r > 0.5) in Bauchi, Gusau, Ilorin, and Ibi, but was not suitable (r < 0.5) at Bida, Kaduna, and Minna. PERSIAN and TAMSAT exhibited weak correlations (r < 0.4) at all locations, except that TAMSAT performed well at Bauchi (r = 0.6). However, low bias ( Figure 7) and RMSE ( Figure 8) were mostly prevalent over the Savanna zone.
In the Guinea coast zone (Figure 6), CHIRPS was more suitable and shows reasonable correlation with gauge datasets except at Ikom, where it exhibited weak correlation (r = 0.07) in comparison to PERSIANN(r = 0.3) and TAMSAT (0.1). Consequently, the correlation between CHIRPS and gauge were found to be significant at 95% confidence level, except at Ikom and Calabar. Furthermore, CHIRPS and TAMSAT exhibit bias of 26.3% and 20.3%, respectively. The RMSE values of above 78 mm were observed in Benin, Calabar, and Ikom were consistent for all products. Generally, the CHIRPS products showed good results during the JJAS rainfall season as compared to other products under study.

Inter-Annual Variation of Mean Annual Rainfall
The results in this section reveal the inter-annual variation of annual rainfall (mm/year) of all products. Findings from Figure 9 show that CHIRPS was able to reproduce observed annual rainfall with good correlation values (r > 0.4) at Bida, Ikom, Kaduna, and Katsina. A greater portion of the Savannah zone was not well captured by the PERSIANN and TAMSAT data. The northwest region of the Sahel zone was also poorly represented by the TAMSAT product, exhibiting low correlation of less than 0.4. The spatial plot of satellite products bias from gauge datasets is depicted in Figure 10. Both CHIRPS and TAMSAT show the same pattern of bias with model overestimations of gauge datasets

Inter-Annual Variation of Mean Annual Rainfall
The results in this section reveal the inter-annual variation of annual rainfall (mm/year) of all products. Findings from Figure 9 show that CHIRPS was able to reproduce observed annual rainfall with good correlation values (r > 0.4) at Bida, Ikom, Kaduna, and Katsina. A greater portion of the Savannah zone was not well captured by the PERSIANN and TAMSAT data. The northwest region of the Sahel zone was also poorly represented by the TAMSAT product, exhibiting low correlation of less than 0.4.

Inter-Annual Variation of Mean Annual Rainfall
The results in this section reveal the inter-annual variation of annual rainfall (mm/year) of all products. Findings from Figure 9 show that CHIRPS was able to reproduce observed annual rainfall with good correlation values (r > 0.4) at Bida, Ikom, Kaduna, and Katsina. A greater portion of the Savannah zone was not well captured by the PERSIANN and TAMSAT data. The northwest region of the Sahel zone was also poorly represented by the TAMSAT product, exhibiting low correlation of less than 0.4. The spatial plot of satellite products bias from gauge datasets is depicted in Figure 10. Both CHIRPS and TAMSAT show the same pattern of bias with model overestimations of gauge datasets The spatial plot of satellite products bias from gauge datasets is depicted in Figure 10. Both CHIRPS and TAMSAT show the same pattern of bias with model overestimations of gauge datasets visible in southwestern and northeastern Nigeria. Both products overestimated gauge record in Savannah and the northwestern region of Nigeria.  The pattern of the RMSE (Figure 11) for all satellite products show a gradient from south to north, with RMSE values of above 300 mm within the Guinea zone (Benin, Calabar, Ikom) and Sahel (Kano). The pattern of the RMSE (Figure 11) for all satellite products show a gradient from south to north, with RMSE values of above 300 mm within the Guinea zone (Benin, Calabar, Ikom) and Sahel (Kano). The pattern of the RMSE (Figure 11) for all satellite products show a gradient from south to north, with RMSE values of above 300 mm within the Guinea zone (Benin, Calabar, Ikom) and Sahel (Kano).

Conclusions
In recent times, the easy accessibility of high-quality satellite rainfall products over Africa has stirred up research activities in rainfall-dependent sectors. Most applications at regional scales showed promising results and underscored the utility of these products in data-scarce regions. In this study, the capabilities of the CHIRPS, PERSIANN, and TAMSAT to reproduced local rainfall characteristics in Nigeria from 1983 to 2013 were evaluated at point-location scale. Satellite products performances were analyzed at seasonal, inter-seasonal, and inter-annual time scales using statistical approaches.
Findings from the spatial distribution of seasonal rainfall showed that all of the products followed the same pattern, depicting the south (high) to north (low) gradient of rainfall amount. All products visibly captured this trend during the early (MAM) and late (SON) season, but poorly reproduced it during the heavy rainfall period (JJA and JJAS). The satellite products captured the unimodal and bimodal trend of rainfall peaks, with high RMSE values observed in locations within the Guinea coast climatic zone. The products with varying magnitude of accuracy simulated the interannual variations of local rainfall anomaly. The pattern of frequency of monthly rainfall occurrence was also in agreement with gauge datasets, but with varying accuracy. It is worthy to note that

Conclusions
In recent times, the easy accessibility of high-quality satellite rainfall products over Africa has stirred up research activities in rainfall-dependent sectors. Most applications at regional scales showed promising results and underscored the utility of these products in data-scarce regions. In this study, the capabilities of the CHIRPS, PERSIANN, and TAMSAT to reproduced local rainfall characteristics in Nigeria from 1983 to 2013 were evaluated at point-location scale. Satellite products performances were analyzed at seasonal, inter-seasonal, and inter-annual time scales using statistical approaches.
Findings from the spatial distribution of seasonal rainfall showed that all of the products followed the same pattern, depicting the south (high) to north (low) gradient of rainfall amount. All products visibly captured this trend during the early (MAM) and late (SON) season, but poorly reproduced it during the heavy rainfall period (JJA and JJAS). The satellite products captured the unimodal and bimodal trend of rainfall peaks, with high RMSE values observed in locations within the Guinea coast climatic zone. The products with varying magnitude of accuracy simulated the inter-annual variations of local rainfall anomaly. The pattern of frequency of monthly rainfall occurrence was also in agreement with gauge datasets, but with varying accuracy. It is worthy to note that CHIRPS performed better than other products in reproducing local rainfall climatology inmost locations.
Within the Sahel, a significant increasing trend was accurately reproduced by all products at the annual time series except at Gusau, where only CHIRPS and gauge datasets exhibited a similar trend. At the seasonal time series, variations in capturing significant trends exist at different locations and for different products. Generally, product evaluation statistics show that CHIRPS scored high correlation (r > 0.5) in many of the locations within all zones. However, poor agreement between all products against observed data are seen around the southern and northwestern axis of Nigeria.