Rainfall Comparison from Different Precipitation Estimates over West Africa ()
1. Introduction
Recent understanding of rainfall estimation from satellite sensors was achieved from past intercomparison programs, which were conducted during the last two decades, to investigate the various rainfall algorithms. For instance, the third Algorithm Intercomparison Program (AIP-3) coordinated by the GPCP in the framework of the World Climate Research Program (WCRP) analysed and compared different rainfall estimates from over 50 algorithms with ground- based radar data in the TOGA-COARE oceanic region from November 1992 to February 1993. While detailing the results from AIP-3, [1] indicated that there is a better correlation of the Microwave (MW) techniques with the validation data when compared with the Infra-Red (IR) or combined IR/MW techniques for instantaneous rainfall estimation, whereas, for monthly rainfall mean, the combined IR/MW techniques performed better than the other rainfall estimation techniques [1] .
Other studies, both global and regional, indicated the relevance of techniques that combine gauge and satellite data when compared to techniques not applying any calibration (e.g. [2] and [3] ). As regard the various rainfall products obtained from TRMM data detailed in [4] , [5] presented the results obtained after two years in orbit of TRMM satellite and the deviations among the Microwave (MW), Infra-Red (IR), and combined IR/MW.
Advanced algorithms have been proposed by several authors, which were developed to produce surface rainfall estimates averaged on short time scales to satisfy the larger scientific community. For example, [6] conducted a validation study using surface data from two different dense rain gauge networks located in Tropical Africa, at various temporal resolution, for surface rainfall products obtained by using different source data. Three surface rainfall satellite products were considered at the 1-day timescale; 1) The RACC product obtained from geo-IR data calibrated by TRMM Microwave Imager (TMI) Microwave (MW) data accumulated for 1-day and 1-degree grid boxes, 2) The TRMM 3B42 product, which combines geo-IR and TMI-MW data, and 3) The GPCP one degree-daily product combines geo-IR, SSMI MW and gauge analyses. A microwave-based satellite rainfall product was also considered: it was computed from the MW data delivered by all the available SSMI’s, but it is cumulated for 5-day.
[7] showed that the correlation coefficients, from [6] , were low for the three satellite products at the daily time scale, and coefficients of about 0.7 were reached with time integration for time accumulation of at least 5 days. The performance obtained for GPCP one degree daily is slightly better than for TRMM 3B42, showing the improvement introduced by the use of the gauge analysis, while the RACC product which is mainly from the IR data were less performing, according to [7] . The SSMI 5-day product showed poor performances, mainly due to the very bad time sampling of MW data. [7] demonstrated that the performances get better when all the source data (satellite geo-IR and MW data and gauge data) were used together.
Similar results regarding the performance of short timescale satellite rainfall products were given by [8] who compared the GPCP one degree daily precipitation field with the high-resolution precipitation data from dense surface networks of the BALTEX countries. For precipitation amounts, the comparison results indicated differences between gauge data and satellite data, which were higher on a daily time scale than on a monthly time scale, and indicated that the resolution of one degree and one day could be too high for quantitative applications. Nevertheless the GPCP one degree daily estimates showed well the time development of the large-scale spatial precipitation systems.
[9] performed an extensive evaluation of two categories of products. The first category had low spatial (2.5˚) and temporal (monthly) resolution and included the GPCP) the National Oceanographic and Atmospheric Administration Climate Prediction Center (NOAA-CPC) merged analysis (CMAP), and the TRMM-3B43. The second category comprised relatively high spatial (0.1˚ to 1˚) and temporal (3-hourly to 10-daily) resolution products. These included the NOAA-CPC African rainfall estimation algorithm, GPCP one-degree-daily (1DD), TRMM-3B42, Tropical Applications of Meteorology using SATellite and other data (TAMSAT) estimates, and the CPC morphing technique (CMORPH). These products were aggregated to a 10-day total and rescaled to spatial resolutions of 1˚, 0.5˚ and 0.25˚. [9] revealed that TRMM-3B43 and CMAP from the first category and CMORPH, TAMSAT and TRMM-3B42 from the second category performed reasonably well.
[10] described a detailed investigation of a new high-resolution gauge―satellite-based analysis of daily precipitation over continental South America during 2004. The methodology adopted by [10] is based on a combination of additive and multiplicative bias correction schemes to get the lowest bias when compared with the observed values. The technique by [10] appeared to be a suitable tool to generate high-resolution and real-time, gauge-based and satellite-based analyses of daily precipitation over land in various regions.
[11] performed a comparison of precipitation estimates from near-real-time of satellite observations and numerical models. The results from [11] confirmed that the reliability of satellite precipitation estimates and the model quantitative precipitation forecasts (QPFs) are both highly dependent on the rainfall regime, and essentially opposed to each other―numerical models are most accurate during winter and at higher latitudes, while satellite estimates of rainfall show the greatest skill during summer and at lower latitudes. [11] further showed that the more the precipitation regime tends toward deep convective system, the less (more) accurate the model (satellite) estimates will be.
[12] evaluated four satellite-derived rainfall products (CMORPH, PERSIANN, TMPA-RT, and TMPA-V6) with regards to their reliability in capturing extreme precipitation. The results of [12] showed that CMORPH and PERSIANN rainfall products indicated better precipitation estimates. However, the false alarm ratio (FAR) and volume of CMORPH and PERSIANN data sets are higher than those of TMPA-RT and TMPA-V6. According to [12] , no single rainfall product can be considered perfect for detecting extreme precipitation events.
[13] assessed the suitability of NASA’s TRMM 3B42 Version 6 rainfall product, and the satellite-based NOAA Climate Prediction Center (CPC) morphing technique (CMORPH), for assimilation applications in hydrological land data. [13] showed that at seasonal or annual time scales, CMORPH has much higher biases and RMS errors than TRMM 3B42. CMORPH shows season-dependent biases, which underestimate in winter and overestimate in summer. This leads to 50% higher RMS errors in CMORPH’s area-averaged daily precipitation than TRMM 3B42. At higher temporal resolution (5 days or less), CMORPH uncertainty is slightly less, and about 10% - 20% higher chances of detection of rainfall events than TRMM 3B42. In addition, the satellite estimates detect more high-intensity events, causing a remarkable shift in precipitation spectrum.
2. Data and Methodology
Three instantaneous (level 3) surface rainfall products (i.e. 3A12, 3A25, and 3B43 in TRMM data classification) from the TRMM precipitation radar (TPR) and GPCC Monthly regional rain gauge data (derived from the Global Precipitation Climatology Centre) were used in this research which span between 1998-2014. The 3A12 dataset is based on the TRMM 2A12 TMI retrieval algorithm. The 3A12 data sets are monthly averaged values of the surface rain rate (or intensity). The 3A25 dataset is based on the TRMM 2A25 PR retrieval algorithm. This dataset are monthly averaged values of the surface rain rate (or intensity). The 3B43 is TRMM merged product with other satellite product (3B42). Monthly regional rain gauge data are derived from the Global Precipitation Climatology Centre (GPCC). The analysis is performed on a 2.5˚ × 2.5˚ grid based on an objective analysis technique presented by [14] . This resolution makes it compatible with most blended and satellite products. The GPCC gauge-only dataset is based on a network of stations that varies from about 500 to 600 stations per month in 1988 to approximately 2500 stations in 1994. The data sources are quality-controlled globally via exchanged synoptic weather reports (‘‘SYNOP’’) and monthly climate reports (‘‘CLIMAT’’). The GPCC gauge dataset is interpolated to points that have no gauges [15] .
Although the PR (Precipitation Radar) and TMI (TRMM Microwave Imager) sensors are on board the same satellite, significant differences exist in inter-an- nual variability of regional mean rainfall between rain estimates from the PR and TMI. As a result, intensity and frequency of TRMM and multi-satellite/gauge merged climate rainfall products were analysed and compared over West Africa using GPCC (Global Precipitation Climatology Center) data as ground truthing.
Analysis and plots are restricted to two seasonal time scales of March-May (MAM) and June-September (JJAS), because of low rainfall in the other months of October-February. Two kinds of validation methodology were carried out. In the first, the rainfall intensity and frequency fields for MAM and JJAS season for various TRMM products were compared.
The second exercise/methodology is the estimation of error statistics for MAM and JJAS season. The bias is calculated as the difference between the satellite mean and the gauge mean, and the rms error is calculated after the bias is removed from the satellite estimate.
The statistical approaches used for the validation are as follows:
1) the Root Mean Square Error (RMSE), which measures the average error magnitude but gives greater weight to the larger errors, expressed as:
(1)
where xi indicates the estimated value at a point or grid box i, yi indicates the observed value, and n is the number of samples.
2) The mean absolute error (MAE) measures the average magnitude of the error. It is given as:
(2)
3) The linear correlation coefficient (CC). this is expressed as:
(3)
where xi indicates the estimated value at a point or grid box i, yi indicates the observed value, and n is the number of samples.
The map of the study area, showing the three zones under investigation, Figure 1 is as shown below.
3. Results and Discussion
3.1. Frequency and Intensity Analysis Using Cumulative Distribution Functions during Mam and JJAS
Figures 2(a)-(c) showed the frequency and intensity distribution for all
Figure 1. Showing the three zones in the study area.
(a)(b)(c)
Figure 2. (a)-(c) Showing Cumulative Distribution Functions (CDF) in Guinea, Savannah, and Sahel during MAM.
the four (4) rainfall products in each zone under investigation. In Figure 2(a) (Guinea), the 3A12 rainfall product indicated higher intensity than the other three (3) rainfall products. It is obvious that, below 50% of the distribution, frequency and intensity are almost the same for the other three rainfall products (3A25, 3B43 and GPCC), with a little spread, between 50% and 90% of the distribution. All the rainfall products reached peak intensities between 12 mm/day and 17 mm/day, but with very few occurrences (frequencies). In Figure 2(b) (Savannah) above, the same pattern is observed as in Figure 2(a) (Guinea). Figure 1(b) indicated that the 3A12 rainfall product estimated rainfall at a higher intensity than the other three (3) products. The 3A12 and 3A25 rainfall products, which are pure TRMM rainfall products from TRMM’s TMI and TRMM’s PR respectively, showed higher number of occurrences (frequencies) of high intensity rainfall than the gauge (GPCC) and the 3B43 rainfall products. A close observation shows that the overall rainfall intensity had decreased in Savannah when compared to Guinean zone. Figure 2(c) represents the CDF of the rainfall products in Sahel region. The Figure showed that the overall intensity had decreased further when compared to Savannah and Guinean Zones. About 70% of the rainfall, in this zone, is below 0.5 mm/day during MAM, 30% of the 70% is completely 0 mm/day during the MAM period. The maximum observed or estimated rainfall in this zone is below 4 mm/day as indicated by the 3A12 and 3A25 rainfall products. The 3B43 and gauge (GPCC) rainfall products indicated maximum intensities below 2 mm/day in the Sahel zone. With this, it can be concluded that rainfall intensity generally decreases with latitude during March to May (MAM) period in West Africa.
The same investigation is performed for June to September (JJAS) period/ season, which is mostly characterised by increased rainfall when compared to the MAM period. In Figure 3(a) (Guinean zone), the gauge (GPCC) and the 3B43 rainfall products had higher intensities compared to the other two (2) rainfall products (3A12 and 3A25). 10% - 95% (representing 85% of the total rainfall in this zone), had intensities between 2.5 mm/day and 12.5 mm/day. This observation indicates that rainfall is most frequent within these intensities range of 2.5 mm/day and 12.5 mm/day in the Guinean zone during JJAS. Only 5% of the total rainfall had intensities within the range of 12.5 mm/day and 20 mm/day, as observed by all the rainfall products in Guinea. In Figure 3(b) (Savannah zone), the highest observed or estimated intensities were indicated by the 3A25 and 3B43 rainfall products, but the intensities observed or estimated by the 3A12 and the 3B43, between 0% - 99%, were generally higher than the other two (2) rainfall products. Both Guinea and Savannah zones peak between 20 mm/day and 25 mm/day. The pattern of the distribution is slightly different for the Sahel zone, as shown in Figure 3(c) (Sahel) above. In this zone, about 20% of the total rainfall distribution indicated no rainfall (0 mm/day intensity). 20% - 99% (about 80% proportion) had rainfall intensities below 8 mm/day while only about 1% of rainfall occurrences had intensities between 8 mm/day and maximum of 14 mm/day, on the average, in the Sahel.
(a)(b)(c)
Figure 3. (a)-(c) Showing Cumulative Distribution Functions (CDF) in Sahel during JJAS.
3.2. Average Rainfall Intensity during MAM and JJAS
The comparison of mean rainfall intensity of all the four (4) rainfall products for each zone during MAM and JJAS are shown in Figure 4(a)-(c). The results, as shown in Figure 4(a) and Figure 4(c), indicated that, on the average, rainfall intensity decreases from Guinea to Sahel in MAM and JJAS, but in JJAS the rainfall intensity in Guinea and Savannah are almost the same. In Sahel, as shown in Figure 4(a) during MAM, all the products indicated very low rainfall intensities, on the average, which are almost the same, in Sahel, the 3A12 rainfall product had the highest mean intensity, while in Guinea, the 3A12 also indicated the highest rainfall mean intensity, compared to the other rainfall products. In Figure 4(c) during JJAS, all the rainfall products are almost the same for Sahel zone except for the 3A25 rainfall product. In Savannah, the mean 3B43 and 3A12 rainfall products are almost the same whereas the 3A25 and gauge (GPCC)
(a)(c)
Figure 4. (a) and (c) Showing averaged rainfall intensities during JJAS. Note that (b) and (d) are not shown.
rainfall products exhibited close mean value in Savannah. The Guinea zone during JJAS showed that the gauge (GPCC) rainfall product had the highest mean value, followed by 3B43 rainfall product, then 3A25 rainfall product and finally 3A12 rainfall product.
This had clearly shown that during peak rainfall episodes, as in JJAS, all the rainfall products under estimate rainfall in Guinea and Sahel, but overestimate rainfall in the Savannah zone, whereas during low rainfall episodes, as in MAM, all the rainfall products over estimate rainfall when compared to gauge (GPCC) rainfall product.
3.3. Correlation, Regression and Statistical Errors
The correlation among the rainfall products are tabulated in Table 1 below. The Gauge (GPCC)-3A12 correlation is the least observed correlation (0.707329) or relationship among the products in Guinea during MAM, while the Gauge (GPCC)-3B43 relationship has the highest correlation (0.874176) in Guinea during MAM. The Gauge (GPCC)-3A25 relationship has the least correlation (0.8186) in Savannah during MAM, also the Gauge (GPCC)-3B43 relationship has the highest correlation (0.919766) in Savannah during MAM. The table also showed that the Gauge (GPCC)-3A25 relationship had the least correlation in Sahel during MAM, and as inferred in Guinea and Savannah, the Gauge (GPCC)-3B43 relationship is nearly perfectly correlated in Sahel During MAM. With this tabulated information, it can be deduced that the 3A12 rainfall product/algorithm does not perform well in the Guinea zone which is closer to the Atlantic Ocean but improved in observation inland. The 3B43 rainfall product/algorithm is the best performing algorithm among the rainfall products and it also improved inland. The 3A25 rainfall product which had the least relationship with Gauge (GPCC) also improved inland. All the rainfall products showed
Table 1. Correlation between each of the Rainfall Products during MAM.
very strong positive correlation with Gauge (GPCC) in all the zones in West Africa during the March-May (MAM) period.
Table 2 below showed that the Root Mean Square Error (RMSE) is always larger than the Mean Absolute Error (MAE). The lower the value of the RMSE, the better is the relationship, and so the results tabulated in Table 2 revealed that the RMSE decreases from Guinea to Sahel for each of the rainfall products. It is also observed that the RMSE decreases from 3A12 to 3A25, and then from 3A25 to 3B43. Thus 3B43 has the least error among the rainfall products.
Because rainfall amount, frequency, and intensity are high during the JJAS season, especially in the Guinean zone, results Tabulated in Table 3 showed that correlation is moderate between Gauge (GPCC) and 3A12 rainfall product, with an improved relationship between Gauge (GPCC) and 3A25, which improved further between Gauge (GPCC) and 3B43 rainfall product. In Savannah, also in Table 3, the correlation also improved with rainfall product as were observed in Guinea zone, but in Savannah the correlation is much better than were observed in Guinea zone. The relationship in Sahel is a little different. In Sahel, Gauge
Table 2. Root Mean Square and Mean Absolute errors of each product with respect to GPCC.
Table 3. Correlation between each of the Rainfall Products during JJAS.
(GPCC)-3A12 had a slightly better correlation than the Gauge (GPCC)-3A25 correlation. The Gauge (GPCC)-3B43 correlation maintained the best relationship with Gauge (GPCC) among the rainfall products.
The Root Mean Square Error (RMSE) for JJAS, in Table 4, also decreases with increasing latitude in West Africa, and also decreases with the rainfall products from 3A12 to 3A25, to 3B43 rainfall product, which implied that the 3B43 rainfall product estimate rainfall better than the other two (2) rainfall products when compared with Gauge (GPCC) rainfall product.
4. Conclusions
This research work had compared observed frequencies and intensities of rainfall estimates from different measurement systems, comprising Tropical Rainfall Measurement mission (TRMM, 3A12 and 3A25 algorithms), the 3B43 rainfall product (The 3B43 combines the estimates generated by the TRMM and other satellites product (3B42) and the CAMS global gridded rain gauge data, produced by NOAA’s Climate Prediction Center (CPC) and/or the global rain gauge product produced by the Global Precipitation Climatology Center (GPCC)), and the GPCC as in-situ data for ground truthing.
Rainfall intensity generally decreases with latitude during March to May (MAM) period in West Africa, with the TRMM’s 3A12 rainfall product showing stronger intensity, in Guinea, Savannah and Sahel, than the other three rainfall products. Both Guinea and Savannah zones peak between 20 mm/day and 25 mm/day. The pattern of the distribution is slightly different for the Sahel zone, as about 20% of the total rainfall distribution indicated no rainfall (0 mm/day intensity). 20% - 99% (about 80% proportion) had rainfall intensities below 8 mm/day while only about 1% of rainfall occurrences had intensities between 8 mm/day and maximum of 14 mm/day, on the average, in the Sahel. The 3A12 rainfall product over estimate rainfall during MAM period, but under estimate rainfall in Guinea during the JJAS season. The 3A12 rainfall product also over estimate rainfall in Savannah during the JJAS season, but in Sahel both during MAM and JJAS, the 3A12 rainfall had close average values with the gauge
Table 4. Root Mean Square and Mean Absolute errors of each product with respect to GPCC.
(GPCC, which represented ground truth observation) rainfall product. The other rainfall products (3B43 and 3A25) also over estimate rainfall in Savannah both during MAM and JJAS, but under estimate rainfall in Guinea during JJAS.
This had clearly shown that during peak rainfall episodes, as in JJAS, all the rainfall products under estimate rainfall in Guinea and Sahel, but overestimate rainfall in the Savannah zone, whereas during low rainfall episodes, as in MAM, all the rainfall products over estimate rainfall when compared to gauge (GPCC) rainfall product.
A view of the intensities and frequencies of each of the rainfall products in each of the zones indicated that rainfall intensities are generally high in Guinea when compared to the other two (2) zones; also intensities are generally high in Savannah when compared to Sahel zone, as indicated by each of the rainfall product during MAM. During JJAS, 3A12 rainfall product estimated stronger intensities, in Savannah zone, than the other two zones. The 3A25 and 3B43 rainfall products exhibited almost the same intensities and frequencies for the Guinea and Savannah zones, whereas, the GPCC indicated stronger intensities and frequencies in the Guinea zone when compared to the Savannah and Sahel zones.
The Gauge (GPCC)-3B43 had the best relationship (highest correlation) in all the three zones during MAM. The 3B43 rainfall product/algorithm is the best performing algorithm among the rainfall products and it also improved inland. The 3A25 rainfall product which had the least relationship with Gauge (GPCC) also improved inland. All the rainfall products showed very strong positive correlation with Gauge (GPCC) in all the zones in West Africa during the March-May (MAM) period. The Gauge (GPCC)-3B43 correlation maintained the best relationship with Gauge (GPCC) among the rainfall products, during JJAS.