Seasonal Effect on Spatial and Temporal Consistency of the New GPM-Based IMERG-v 5 and GSMaP-v 7 Satellite Precipitation Estimates in Brazil ’ s Central Plateau Region

This study assesses the performance of the new Global Precipitation Measurement (GPM)-based satellite precipitation estimates (SPEs) datasets in the Brazilian Central Plateau and compares it with the previous Tropical Rainfall Measurement Mission (TRMM)-era datasets. To do so, the Integrated Multi-satellitE Retrievals for GPM (IMERG)-v5 and the Global Satellite Mapping of Precipitation (GSMaP)-v7 were evaluated at their original 0.1◦ spatial resolution and for a 0.25◦ grid for comparison with TRMM Multi-satellite Precipitation Analysis (TMPA). The assessment was made on an annual, monthly, and daily basis for both wet and dry seasons. Overall, IMERG presents the best annual and monthly results. In both time steps, IMERG’s precipitation estimations present bias with lower magnitudes and smaller root-mean-square error. However, GSMaP performs slightly better for the daily time step based on categorical and quantitative statistical analysis. Both IMERG and GSMaP estimates are seasonally influenced, with the highest difficulty in estimating precipitation occurring during the dry season. Additionally, the study indicates that GPM-based SPEs products are capable of continuing TRMM-based precipitation monitoring with similar or even better accuracy than obtained previously with the widely used TMPA product.


Introduction
Precipitation is a key component of the water cycle, which is facing unprecedented changes related to both climate change and human growth population.Climate change redistributes precipitation's seasonality [1] and intensity [2], threatening the six-fold increase in water extraction observed in the 20th century in response to the increasing world population, food needs, and economics level (vision, water council, 2000).In this context, The World Meteorological Organization (WMO) defines precipitation as an Essential Climate Variable (ECV) to understand and adapt to these ongoing changes.In a general way, precipitation information obtained from point gauges measurement is spatialized to follow precipitation in time and space.However, in relation to (i)-the access difficulty (forest, desert, mountainous), and/or to (ii)-political, social, and economic issues, the gauges network is often too sparse and unevenly distributed to adequately follow precipitation in space and time [3,4].It could be even more marked over tropical regions due to high spatial and temporal precipitation variability.In this context, with an almost global coverage, remote sensing data provide relevant information to overcome these obstacles.Since the first meteorological satellite, Television and Infra-Red Observation Satellite (TIROS), launched in April 1960, a large number of sensors have been developed and used to monitor the precipitation from space [5].They fall into two main categories: visible and infrared (VIS/IR) sensors available from geostationary (GEO) and low-earth orbiting (LEO) satellites; and passive microwave sensors (PMW) available from LEO satellites [6].VIS/IR-based precipitation estimates use the relation between cloud-top temperature and the probability and intensity of precipitation, whereas PMW-based precipitation estimates use a more direct method related to the precipitation-sized particle.Although PMW-based precipitation estimates are a more accurate source of precipitation than VIS/IR-based estimates, the spatial and temporal coverage of LEO limits their application [7].Indeed, a short-term precipitation event can be lost by the PMWs sensors as a result of the irregular and limited overpass of LEO [8].VIS/IR sensors on GEO satellites can overcome this limitation since they provide high-temporal-resolution global data.
In this context, much research has been conducted seeking the development of algorithms to combine both PMW and VIS/IR observations to provide operational satellite precipitation estimates (SPEs) [5][6][7].
These efforts were initiated after the Tropical Rainfall Measurement Mission (TRMM) was launched, in 1997.TRMM is an international project of the National Aeronautics and Space Administration (NASA) and Japan Aerospace Exploration Agency (JAXA), from which a first generation of SPEs products was developed.TRMM-based SPEs include the TRMM Multi-Satellite Precipitation Analysis (TMPA) [9], the Climate Prediction Center Morphing technique (CMORPH) [10], the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) [11], and the Global Satellite Mapping of Precipitation (GSMaP) [12].
To ensure continuity of the TRMM mission, the Global Precipitation Measurement (GPM) mission's Core Observatory was deployed in February 2014, with new technologies and improved global coverage as a result of opportunities from existing satellite programs and new partners obtained from bilateral agreements with either NASA or JAXA [13].From the GPM project, a second generation of SPEs was made available to ensure continuity in precipitation monitoring started by the TRMM-based SPEs.GPM-based SPEs includes the integrated multi-satellite retrievals for GPM (IMERG) developed by NASA [14,15] and a new GSMaP version developed by JAXA [12].
If the TRMM-based SPEs have proved effective in precipitation estimation over the world, as reviewed by [16,17], GPM-based SPEs' reliability is being currently assessed and discussed by the scientific community.In the last 5 years, successive versions of the GPM-based SPEs were released.Up to date, three IMERG versions (IMERG-v3, -v4 and -v5) and two GSMaP versions (GSMaP-v6 and -v7) are available.In this context, timely studies reporting on the reliability of the successive released versions have been made available by the scientific community to contribute to the enhancement of the quality of the successively released versions.
First studies focused on the first GPM-based SPEs, IMERG version (IMERG-v3), to report on its potential improvement in precipitation estimate in comparison to its TRMM-based SPEs predecessors.A first attempt was made at the global scale by comparing IMERG-v3 with its predecessor TMPA-v7 at the monthly time step to highlight the difference in both datasets according to surfaces and precipitations rates [18].Since then, different studies have assessed IMERG-v3's reliability at more local scales in different regions.For example, in India, IMERG-v3 shows notable improvement in comparison to TMPA-v7 to capture heavy rainfall [19] and to represent the mean monsoon rainfall and its variability [20].In China, IMERG-v3 is more accurate than TMPA-v7 at monthly [21], sub-daily, and daily scales [22] and therefore represents a better option for hydrological applications [23].More generally, an overall improvement from TRMM-based SPEs to IMERG-v3 were also highlighted in Singapore [24], Iran [25], Bolivia [26], United Arab Emirates [27], Tibet [28], East Africa [29], Continental United State (CONUS) [30], Korea, and Japan [31].
The second released IMERG version (IMERG-v4) was compared to its previous version (IMERG-v3) to highlight potential improvement on precipitation estimates accuracy.Interestingly, IMERG-v4 did not exhibit the anticipated improvement when compared to IMERG-v3 in China [32,33] and Pakistan [34] with higher bias, lower Correlation Coefficient and Root Mean Square Error.Still IMERG-v4 provides better precipitation estimates than TRMM-based SPE in Malaysia [35] and Pakistan [36].
It is worth mentioning that even if IMERG-v5 globally provides more accurate precipitation estimates than the TRMM-based SPEs, the opposite can be observed at the very local scale.For example, in Pakistan, at the grid-cell scale, TMPA-v7 and CMORPH-BLD provided more accurate monthly precipitation estimates than IMERG-v5 [45].Similarly, the other GPM-based SPEs (GSMaP-v7) also provided better precipitation estimates than IMERG-v5.However, only a few studies have compared GSMaP-v7 to IMERG-v5.To our knowledge, GSMaP-v7's reliability was only assessed and compared to IMERG-v5 in Myanmar [42], Tibet [46], Pakistan [45], Brazil [43], and CONUS [44].Therefore, there is a need to report on the latest released GPM-based SPEs versions to understand their reliability and potential complementarity.
In the described context, this study aims to evaluate the performance of the latest GPM-based IMERG-v5 and GSMaP-v7 SPEs, for the first time, over the Brazilian Central Plateau region.The assessment was made by comparing their estimates with measured data from precipitation gauges at annual, monthly, and daily scales for both wet and dry seasons, respectively.The TRMM-based TMPA-v7 dataset for a common period of data availability is also considered to assess the potential improvement of IMERG-v5 and GSMaP-v7 compared to TMPA-v7.This timely study aims to accompany the transition from TRMM to GPM-based SPE products and is expected to provide relevant information to ongoing GPM precipitation estimate algorithms and to researchers when selecting the most reliable SPEs for their studies.

Study Area
The study area is situated in the Brazilian Central Plateau, located between latitudes 15.5 • S and 16.1 • S and longitudes 47.3 • S and 48.5 • S, where three South American river basins originate: São Francisco, Tocantins, and Paraná (Figure 1a).The region presents a mean slope value of 5.2 • (9.3%), with elevations ranging from 667 to 1356 m a.s.l.(Figure 1b).The study area is located in the Cerrado biome, which is exceeded in terms of area only by the Amazon and is composed of a mosaic of grass fields, savanna, and forest vegetation formations [47,48].Recently, the Cerrado is the Brazilian region with the most important beef production, and has experienced the main expansion in total cropland area over the last two decades [49].If the actual scenario of large-scale agriculture expansion remains, projections indicate that around 30% of the Cerrado's remaining natural vegetation is likely to be cleared in the next three decades as a consequence of limited legal protection and market pressure [50].The climate is tropical, with a dry winter season from May to September and a rainy summer season in the remaining months, which accounts for more than 90% of average annual precipitation (Figure 1c).observation time; and the final-run, which uses monthly precipitation data from the Global Precipitation Climatology Centre (GPCC) to perform a final adjustment [14].The IMERG datasets are currently at version 5.In this study, only IMERG final-run version 5 (called IMERG hereafter) is considered, as it is the latest version.Early-and late-run versions were not considered, as they consist of intermediate step with an overall low efficiency due to the absence of the gauge's adjustment.Data were downloaded from the internet (ftp://arthurhou.pps.eosdis.nasa.gov/gpmdata/).

Global Satellite Mapping of Precipitation (GSMaP)
The GSMaP project was established in 2002 by the Japan Science and Technology Agency (JST) under the Core Research for Evolutional Science and Technology (CREST).Today, the activities are implemented under the Japanese Precipitation Measuring Mission (PMM) Science Team, which is a joint science team between TRMM and its successor mission, GPM [51].
The GSMaP Project has produced different SPE products from PMW, and has merged IR and PMW datasets [52,12].This study used the 0.1° spatial-resolution GSMaP product available at onehour time interval, with coverage between 60° N and 60° S, and gauge-adjusted with the Climate Prediction Center (CPC) global rain gauge daily database.This dataset uses an approach similar to CMORPH, where a moving vector derived from two successive IR images is used to propagate rainy areas from microwave radiometry.GSMaP focuses on the relationship between precipitation strength and IR brightness temperature, and a Kalman filter is then applied to the precipitation rate to update the rate on the IR brightness temperature [12].Currently, the GSMaP gauge-adjusted product is in its

Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG)
IMERG is a 0.1 • spatial-resolution product available at 30 min time intervals, with spatial coverage between 60 • N and 60 • S. IMERG is produced by NASA using concepts and components from previous algorithms: TMPA, CMORPH, and PERSIANN.Precipitation estimates from available LEO satellites' PMW sensors are assembled, gridded, and intercalibrated, then combined with data from GEO IR sensors using the Morphing-Kalman Filter, CMORPH-KF, Lagrangian time interpolation scheme, and the PERSIANN Cloud Classification System (PERSIANN-CCS) to produce half-hourly estimates [14].
There are three stage level available IMERG datasets: the early-run, providing data approximately five hours after observation time; the late-run, with approximately 16 h delay from observation time; and the final-run, which uses monthly precipitation data from the Global Precipitation Climatology Centre (GPCC) to perform a final adjustment [14].The IMERG datasets are currently at version 5.In this study, only IMERG final-run version 5 (called IMERG hereafter) is considered, as it is the latest version.Early-and late-run versions were not considered, as they consist of intermediate step with an overall low efficiency due to the absence of the gauge's adjustment.Data were downloaded from the internet (ftp://arthurhou.pps.eosdis.nasa.gov/gpmdata/).

Global Satellite Mapping of Precipitation (GSMaP)
The GSMaP project was established in 2002 by the Japan Science and Technology Agency (JST) under the Core Research for Evolutional Science and Technology (CREST).Today, the activities are implemented under the Japanese Precipitation Measuring Mission (PMM) Science Team, which is a joint science team between TRMM and its successor mission, GPM [51].
The GSMaP Project has produced different SPE products from PMW, and has merged IR and PMW datasets [12,52].This study used the 0.1 • spatial-resolution GSMaP product available at one-hour time interval, with coverage between 60 • N and 60 • S, and gauge-adjusted with the Climate Prediction Center (CPC) global rain gauge daily database.This dataset uses an approach similar to CMORPH, where a moving vector derived from two successive IR images is used to propagate rainy areas from microwave radiometry.GSMaP focuses on the relationship between precipitation strength and IR brightness temperature, and a Kalman filter is then applied to the precipitation rate to update the rate on the IR brightness temperature [12].Currently, the GSMaP gauge-adjusted product is in its seventh version, and will be called GSMaP hereafter.GSMaP data were downloaded from the internet (https://www.gportal.jaxa.jp/gp).

Tropical Rainfall Measurement Mission Multi-Satellite Precipitation Analysis (TMPA)
The TMPA was developed by NASA [53,54] and has been widely assessed in many regions [35][36][37][38].TMPA was developed to assemble precipitation estimates from different satellite systems with land-surface precipitation gauges in order to provide the best satellite estimate [53].The TMPA products are produced in four stages: combination and calibration of PMW precipitation estimates using the Goddard Profiling Algorithm (GPROF); IR precipitation estimates are created using the merged PMW data; the PMW and IR are combined; and rain gauge data are re-scaled, calibrated, and used [24,53].
TMPA has a 0.25 • spatial resolution over 50 • N and 50 • S, and is available at three-hour and monthly time steps, which are called 3B42 and 3B43, respectively.TMPA 3B42 has a real-time and an adjusted dataset using gauge-based data from the GPCC and the Climate Assessment and Monitoring System (CAMS) developed by the CPC.This study adopted the TMPA 3B42-v7 adjusted, available on the internet (https://mirador.gsfc.nasa.gov/),hereafter called TMPA.

Rain Gauges
A total of 56 rain gauges with daily records for the period from 12 March 2014 to 30 November 2016 were made available by different institutions which are responsible for quality assessment before making it available to users (Agência Reguladora de Águas, Energia e Saneamento do Distrito Federal (ADASA), Agência Nacional de Águas (ANA), Companhia de Saneamento Ambiental do Distrito Federal (CAESB), Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), and Instituto Nacional de Meteorologia (INMET)).The rain gauge network extends over approximately 6000 km 2 and the stations are distributed on forty 0.1 • grid-cells and on twelve 0.25 • grid-cells (Figure 1d,e).

Pre-processing
In order to compare IMERG and GSMaP with TMPA on a common period of data availability (12 March 2014 to 29 November 2016), it was necessary to aggregate IMERG and GSMaP values from their initial 0.1 • grid to TMPA's 0.25 • grid scale [21].First, both IMERG and GSMaP were resampled from their 0.1 • grid resolution into a 0.05 • resolution dataset, using the nearest neighbor method.Then, from this 0.05 • grid, an average precipitation value was calculated for each set of 5 × 5 grid-cells which were within the overlapping 0.25 • grid [26].Together, these two steps provided 0.25 • grids for IMERG and GSMaP.
In order to compare the SPEs with the in situ precipitation dataset, it was also necessary for them all to have the same temporal resolution.The precipitation values from IMERG, GSMaP, and TMPA, which were taken every 30 min, one hour, and three hours, respectively, were converted into daily (from 12Z to 12Z), monthly, and yearly time steps.
Finally, since the number of precipitation gauges inside the same grid-cell varied from one to four on the 0.1 • grid (Figure 1d), and from one to ten on the 0.25 • grid (Figure 1e), the analysis also required processing the in situ precipitation dataset before comparing it with the SPEs.An average precipitation time-series was generated from each set of gauges within each grid-cell with more than Water 2019, 11, 668 6 of 18 one precipitation gauge inside.It is worth mentioning that the number of gauges can impact the results, due to the localized nature of rainfall measurements obtained by rain gauges and to the localized precipitation events that can occur.The influence of gauge density will be addressed in the discussion section.

Comparison Methodology
According to the considered time step (annual, monthly, daily), different methods were used to compare the SPEs with the observed gauge values and to analyze the difference between them.For the monthly and daily time step assessments, two types of datasets were generated to assess the SPEs potential.The first dataset (hereafter called Mean Regional) consists of an average daily and monthly time series obtained by taking the mean of the time series records of all the grid-cells including at least one gauge.The resulting series count was 994 daily and 32 monthly values for both the 0.1 • and 0.25 • grid scale.The second dataset (hereafter called Mean Spatial) consists of records available from all the grid-cells including at least one gauge in a single series for both daily and monthly time steps.The resulting series counts were 34,799 daily and 1065 monthly values for the 0.1 • grid scale and 11,782 daily and 375 monthly values for the 0.25 • grid scale.The consideration of both Mean Spatial and Mean Regional datasets aims to assess the potential of SPEs at different spatial scales.The Mean Regional series comparisons will provide insight into the potential of SPEs at the regional scale, whereas the Mean Spatial series comparisons will provide feedback on the spatial consistency of SPEs over the considered region.Finally, both datasets (Mean Regional and Mean Spatial) were split in two according to the dry and wet season to observe the response of the SPEs to the climatic seasonality.

Annual Comparison
For the annual time step assessment, the 2015 annual precipitation maps were derived from the gauges and SPEs for the 0.1 • and 0.25 • grid to visually compare precipitation spatial distribution from all considered datasets (gauges and SPEs).Then, the 2015 annual precipitation obtained from all considered SPEs at all the grid-cells including at least one gauge were compared to the values obtained from the gauges.The comparison used the correlation coefficient (CC), the percent root-mean-square error (%RMSE), and the percent bias (%B) (see Equations ( 1)-( 3)): where N is the number of values, P is the in situ observed precipitation, P is the average P for the considered period, SPE is the precipitation estimate of the considered SPEs, and SPE is the average of the values of the SPE parameters for the considered period.

Monthly Comparison
The monthly time step analysis used the CC, %RMSE, and %B Equations ( 1)-(3) for both the Mean Regional and the Mean Spatial datasets.Previous studies assessing SPEs (e.g., [26,[39][40][41]55]) have used the following reference values to define acceptable monthly results: CC ≥ 0.7%; RMSE ≤ 50%; and %B between −10% to 10%.In order to facilitate the SPEs potential inter-comparison, results are also presented in the form of a Taylor diagram [56].The Taylor diagram shows the standard  4) normalized (divided) by the reference STD.
In the Taylor diagram, the reference is presented by a black dot where the STD and the CC are equal to 1 and the RMSEc is equal to 0. Therefore, the closer the SPEs' statistics are to the reference values, the closer their estimates are to the observed precipitation values.Additionally, for all considered SPEs, the CC, %RMSE, and %B are computed for all the grid-cells including at least one gauge.This aims to provide more details in the statistical score obtained from the Mean Regional datasets.The monthly analysis is performed for each dataset (Mean Regional and Mean Spatial) considering the whole period and the dry and wet season separately.

Daily Comparison
The daily time step analysis also used the CC, %RMSE, and %B Equations ( 1)-( 3) for both the Mean Regional and the Mean Spatial datasets.Results are presented in the form of a Taylor diagram to facilitate the SPEs potential inter-comparison.Additionally, SPEs were also evaluated based on categorical statistics considering precipitation events as discrete values with only two possible cases (rainy day or not).Four cases are possible: (a) both SPEs and rain gauge register a precipitation event; (b) only the SPEs reports precipitation; (c) only the rain gauge registers precipitation; and (d) precipitation is not reported on either SPEs or rain gauge.After dividing data into these four classes, it was possible to calculate the probability of detection (POD), the false alarm ratio (FAR), the Critical Success Index (CSI), and the bias (B) using the following equations: POD, FAR, and CSI vary from 0 to 1.For POD, 1 is the perfect score and indicates the SPEs' capacity to correctly detect precipitation within the study area.CSI also has a perfect score of 1, which represents the ratio between rain events correctly detected by the SPEs and the total rain events reported by the SPEs and the rain gauges.FAR has 0 as a perfect score, and is an indicator of how much the SPEs wrongly report precipitation when it does not occur.The success ratio (SR = 1 − FAR) can be derived from FAR, having 1 as a perfect score.Different to POD, FAR, and CSI, B does not vary from 0 to 1. Rather, this index takes 1 as a reference value.If B is higher than 1, it indicates overestimation of the number of rain events within the SPEs product in the studied period, while values lower than 1 indicate underestimation.
A synthesis diagram was also prepared to facilitate the comparison between the categorical statistics of all considered SPEs.The performance diagram [57] builds a response surface from the geometric relationship between the POD, CSI, SR, and B, where the results of each assessed SPEs scenario are plotted.This type of analysis has been conducted in previous studies [26,58], and considerably helps to interpret the global performance of the results.The daily analysis is done for each dataset (Mean Regional and Mean Spatial) considering the whole period and the dry and wet season separately.

Annual Analysis
Figure 2a presents the 2015 annual precipitation maps registered by the SPEs and the gauges.According to the gauge observations, the study area presents a marked west-east gradient, with the highest (lowest) precipitation occurring over the western (eastern) part (Figure 2a).Annual precipitation ranges from 800 to 1500 mm.Despite an overall underestimation of precipitation, with %B close to −20%, GSMaP (0.1 • and 0.25 • ) captured the gradient (Figure 2a,b).On the other hand, IMERG had difficulties identifying the gradient, presenting smaller annual precipitation variability over the study area when compared to GSMaP, with annual precipitation values between 1000 and 1400 mm.TMPA also had difficulties identifying the local gradient and registered even smaller spatial variability, with annual values ranging between 1100 and 1300 mm, resulting in worse statistics when compared to the precipitation gauges.Moreover, the CC results presented statistical significance for all evaluated SPEs (p-value < 0.01), except for TMPA, in Figure 2.
Water 2019, 11, x FOR PEER REVIEW 8 of 18

Annual Analysis
Figure 2a presents the 2015 annual precipitation maps registered by the SPEs and the gauges.According to the gauge observations, the study area presents a marked west-east gradient, with the highest (lowest) precipitation occurring over the western (eastern) part (Figure 2a).Annual precipitation ranges from 800 to 1500 mm.Despite an overall underestimation of precipitation, with %B close to −20%, GSMaP (0.1° and 0.25°) captured the gradient (Figure 2a,b).On the other hand, IMERG had difficulties identifying the gradient, presenting smaller annual precipitation variability over the study area when compared to GSMaP, with annual precipitation values between 1000 and 1400 mm.TMPA also had difficulties identifying the local gradient and registered even smaller spatial variability, with annual values ranging between 1100 and 1300 mm, resulting in worse statistics when compared to the precipitation gauges.Moreover, the CC results presented statistical significance for all evaluated SPEs (p-value < 0.01), except for TMPA, in Figure 2. Despite the SPEs respective ability to represent the regional precipitation pattern, Figure 2 allows an inter-comparison of the potential of SPEs to estimate the annual precipitation amount.When comparing the GPM-based SPEs with the TRMM-based SPEs, both IMERG and GSMaP are more able to quantify the annual precipitation amount.Indeed, in the Taylor diagram (Figure 2c), GSMaP and IMERG for the 0.25° grid are closer to the reference dot than TMPA, which locates farther from the reference dot than all the considered SPEs and resolutions.Interestingly, passing from 0.1° GSMaP and IMERG for the 0.25 • grid are closer to the reference dot than TMPA, which locates farther from the reference dot than all the considered SPEs and resolutions.Interestingly, passing from 0.1 • to 0.25 • resolution enhances the quality of IMERG's and GSMaP's annual precipitation estimates with a slightly increase (decrease) of the CC (%RMSE).This could be because the lower resolution tends to smooth the annual precipitation and thus to produce better results.

Monthly Analysis
Figure 3 shows the potential of SPEs for reproducing monthly mean regional precipitation for the inter-seasonal and for both the dry and wet seasons.The SPEs are able to capture the seasonal cycle from April 2014 to November 2016 with marked dry (from May to September) and wet (from October to April) seasons, following the reference's monthly seasonal tendency (Figure 3a,b).Figure 4 presents the spatial variability of monthly statistics within the study area for the interseasonal temporal scale in the form of %RMSE, %B, and CC for IMERG, GSMaP, and TMPA.It is worth mentioning that only the 0.25° spatial resolution was considered, as similar results were obtained for the 0.1° and 0.25° resolutions for the previous annual and mean monthly regional assessments.
It is possible to notice that most of the assessed grid-cells present %RMSE < 50%, achieving the established quality threshold.TMPA and IMERG present the same %RMSE pattern, whereas GSMaP presents two additional grid-cells with %RMSE between 50 and 75%.Regarding the %B, TMPA (GSMaP) is the SPE with the highest (lowest) proportion of grid-cell, with %B within the established threshold (−10% < %B < 10%).IMERG and TMPA present a similar spatial distribution of %B, with some grid-cells overestimating the reference.On the other hand, GSMaP underestimated Figure 3c presents the mean regional monthly CC, %RMSE, and %B for all considered SPEs, and Figure 3d presents the Taylor diagram, which is intended to facilitate the inter-comparison of SPE potential.The CC values presented in Figure 3 were statistically significant (p-value < 0.01) for all considered SPEs and scenarios.
In the inter-seasonal analysis, all SPEs performed better than during both wet and dry seasons.Indeed, in the Taylor diagram, all SPEs' relative positions are closer to the reference dot for the inter-seasonal dataset.Concerning inter-seasonal evaluation, all values, except the %B for GSMap (0.1 • and 0.25 • grids), are within the established quality threshold (CC ≥ 0.7; %RMSE < 50% and −10% < %B < 10%).As observed for the annual time step, GSMaP tends to underestimate the precipitation amount, with a %B of approximately −10% for both 0.1 • and 0.25 • spatial resolution.Overall, the new GPM-based SPEs presented similar statistics to TRMM-based SPEs in the inter-seasonal analysis, especially IMERG.This can be visually checked in the Taylor diagram, in which all SPEs present very close relative position to the reference dot (Figure 3d).Aggregating IMERG and GSMap datasets from 0.1 • to 0.25 • slightly improved their accuracy, with lower %RMSE and %B for both SPEs.
During the dry season, all SPEs achieved the established quality threshold for CC (CC ≥ 0.7), with CC systematically superior to 0.9 and the highest value obtained for TMPA (CC = 0.96).Regarding the %RMSE, only GSMaP achieved the established quality threshold (%RMSE < 50).However, the %RMSE values of IMERG (0.1 • ) and TMPA (%RMSE = 56%) are very close.With a %B value of 3.2, only IMERG (0.1 • ) achieved the acceptable quality threshold (−10 < %B< 10%).TMPA and IMERG (0.25 • ) overestimate the regional monthly precipitation amount, whereas GSMaP (0.1 • and 0.25 • ) underestimates the regional monthly precipitation (Figure 3).Despite a significant precipitation underestimation for GSMaP 0.1 • and 0.25 • , with %B of −29.1% and −21%, respectively, GSMaP can capture monthly regional precipitation more accurately during the dry season.Indeed, as observed in the Taylor diagram, GSMaP is closer to the reference dot than all other considered SPEs.On the contrary, IMERG is the least efficient SPE, having the farthest position from the dot reference in the Taylor diagram.It is worth mentioning that aggregating from 0.1 • to 0.25 • increased (decreased) GSMaP's (IMERG's) %RMSE and %B.
During the wet season, all considered SPEs present CC, %RMSE, and %B values that are within the established quality threshold (CC ≥ 0.7; %RMSE < 50%; and −10% < %B < 10%) with CC values systematically superior to 0.92 and %RMSE values systematically inferior to 22%.The %B of GSMaP (0.25 • ) was the only exception (%B = −11.1%),confirming the tendency of GSMaP to underestimate monthly regional precipitation amounts, as observed for the inter-seasonal and the dry season.Despite the %B statistics, all SPEs performed similarly, as illustrated in the Taylor diagram, where their positions relative to the reference dot are very close (Figure 3d).
Finally, SPEs performed better during the wet season than during the dry season.As illustrated in Figure 3d, SPEs are closer to the reference dot for the wet season than for the dry season.
Figure 4 presents the spatial variability of monthly statistics within the study area for the inter-seasonal temporal scale in the form of %RMSE, %B, and CC for IMERG, GSMaP, and TMPA.It is worth mentioning that only the 0.25 • spatial resolution was considered, as similar results were obtained for the 0.1 • and 0.25 • resolutions for the previous annual and mean monthly regional assessments.
It is possible to notice that most of the assessed grid-cells present %RMSE < 50%, achieving the established quality threshold.TMPA and IMERG present the same %RMSE pattern, whereas GSMaP presents two additional grid-cells with %RMSE between 50 and 75%.Regarding the %B, TMPA (GSMaP) is the SPE with the highest (lowest) proportion of grid-cell, with %B within the established threshold (−10% < %B < 10%).IMERG and TMPA present a similar spatial distribution of %B, with some grid-cells overestimating the reference.On the other hand, GSMaP underestimated precipitation, with negative %B over most of the considered grid-cells.Regarding the CC, all SPEs systematically present values within the defined quality threshold (CC > 0.7).All CC values presented statistical significance (p-value < 0.01) in Figure 4c,d).IMERG and TMPA presented the same CC pattern over the region, with CC superior to 0.8 over all of the considered region.
Figure 4d shows the SPEs' statistics based on the Mean Spatial dataset in the form of the Taylor diagram.It aims to summarize observations from Figure 4a-c to quantify the SPEs' potential to reproduce monthly precipitation in space for the inter-seasonal, dry, and wet seasons.With relative position closer to the reference dot, all SPEs performed better at the inter-seasonal temporal scale than on a seasonal basis.At the inter-seasonal scale, the new GPM-based SPEs present similar statistics to the TRMM-based SPEs.All three SPEs presented satisfactory statistical scores, with %RMSE, %B, and CC within the defined quality threshold except for GSMaP, which slightly underestimates the precipitation amount (%B = −11.5).Considering the seasonal scale (wet and dry seasons), all SPEs are more efficient during the wet season than during the dry season.From the wet to the dry season, %RMSE is more than three times greater for all SPEs, %B is twice as large for GSMaP and four times as large for both IMERG and TMPA, while a similar decrease in CC is observed for all SPEs.Consequently, none of the considered SPEs presents an %RMSE or %B value within the established quality threshold.
For all the considered temporal scales (inter-seasonal, wet season, dry season), the new GPM-based SPEs presented similar statistics to TRMM-based SPEs, especially IMERG.Indeed, the positions of IMERG and TMPA in the Taylor diagram are very close, and they both present similar %B values.In the Taylor diagram, GSMaP is systematically slightly farther to the reference dot than IMERG and TMPA, highlighting an overall lower ability to represent mean monthly precipitation in space.The main difference between GSMaP and both IMERG and TMPA occurs in terms of %B, with GSMaP underestimating precipitation amount.
Figure 4d shows the SPEs' statistics based on the Mean Spatial dataset in the form of the Taylor diagram.It aims to summarize observations from Figure 4a-c to quantify the SPEs' potential to reproduce monthly precipitation in space for the inter-seasonal, dry, and wet seasons.With relative position closer to the reference dot, all SPEs performed better at the inter-seasonal temporal scale than on a seasonal basis.At the inter-seasonal scale, the new GPM-based SPEs present similar statistics to the TRMM-based SPEs.All three SPEs presented satisfactory statistical scores, with %RMSE, %B, and CC within the defined quality threshold except for GSMaP, which slightly underestimates the precipitation amount (%B = −11.5).Considering the seasonal scale (wet and dry seasons), all SPEs are more efficient during the wet season than during the dry season.From the wet to the dry season, %RMSE is more than three times greater for all SPEs, %B is twice as large for GSMaP and four times as large for both IMERG and TMPA, while a similar decrease in CC is observed for all SPEs.Consequently, none of the considered SPEs presents an %RMSE or %B value within the established quality threshold.
For all the considered temporal scales (inter-seasonal, wet season, dry season), the new GPMbased SPEs presented similar statistics to TRMM-based SPEs, especially IMERG.Indeed, the positions of IMERG and TMPA in the Taylor diagram are very close, and they both present similar %B values.In the Taylor diagram, GSMaP is systematically slightly farther to the reference dot than IMERG and TMPA, highlighting an overall lower ability to represent mean monthly precipitation in space.The main difference between GSMaP and both IMERG and TMPA occurs in terms of %B, with GSMaP underestimating precipitation amount.

Daily Analysis
The daily analysis results are illustrated in Figure 5.The daily quantitative results are summarized in the form of a Taylor diagram for both Mean Regional (Figure 5a) and Mean Spatial (Figure 5c) precipitation datasets.The daily contingency assessment is summarized in the form of a performance diagram for both Mean Regional (Figure 5b) and Mean Spatial (Figure 5d) precipitation datasets.
The daily analysis results are illustrated in Figure 5.The daily quantitative results are summarized in the form of a Taylor diagram for both Mean Regional (Figure 5a) and Mean Spatial (Figure 5c) precipitation datasets.The daily contingency assessment is summarized in the form of a performance diagram for both Mean Regional (Figure 5b) and Mean Spatial (Figure 5d) precipitation datasets.The Taylor diagram shows an overall low performance for all SPEs and both Mean Regional (Figure 5a) and Mean Spatial (Figure 5c) precipitation datasets.Indeed, all SPEs registered very low CC values, which were systematically inferior to 0.4 and 0.5 for the Mean Regional and Mean Spatial datasets, respectively.All CC values presented statistical significance (p-value < 0.01) in Figure 5a,c).Aggregation from 0.1° to 0.25° did not significantly improve the considered statistical metrics.Still, as observed previously in the monthly assessment, the wet season presents better statistics than the dry season and the best results are obtained for the inter-seasonal temporal scale.IMERG presents similar statistics to TMPA, with close positions on the Taylor diagram (Figure 5a,c), and GSMaP appears as the most effective SPE, with closer position to the reference dot on the Taylor diagram for all considered temporal scales (inter-seasonal, wet, and dry seasons).
Different to the daily quantitative statistical results, whose performance was unsatisfactory, as reported in the Taylor diagram (Figure 5a,c), all SPEs present a good ability to detect daily precipitation events (Figure 5b,d).The statistical scores are represented in the form of a performance diagram, in which geometric relationships between POD, SR, B, and CSI facilitate the intercomparisons between the SPEs.The closer the SPEs are to the right-hand top corner black dot, the higher is the SPEs' potential to detect daily precipitation events.Overall, IMERG presented the best statistical scores according to the performance diagram for the Mean Regional dataset assessment (Figure 5b), while GSMaP presented the best results considering the Mean Spatial dataset (Figure 5d).
The ability of SPEs for detecting precipitation events is seasonally dependent.Indeed, for both datasets (Mean Regional and Mean Spatial), the best scores in the performance diagram are reached during the wet period, whereas the lowest scores are reached during the dry season.On the other hand, intermediate scores are reached for the inter-seasonal scale.Interestingly, aggregating IMERG and GSMaP from the 0.1° grid to the 0.25° grid yielded better results considering the Mean Spatial The Taylor diagram shows an overall low performance for all SPEs and both Mean Regional (Figure 5a) and Mean Spatial (Figure 5c) precipitation datasets.Indeed, all SPEs registered very low CC values, which were systematically inferior to 0.4 and 0.5 for the Mean Regional and Mean Spatial datasets, respectively.All CC values presented statistical significance (p-value < 0.01) in Figure 5a,c).Aggregation from 0.1 • to 0.25 • did not significantly improve the considered statistical metrics.Still, as observed previously in the monthly assessment, the wet season presents better statistics than the dry season and the best results are obtained for the inter-seasonal temporal scale.IMERG presents similar statistics to TMPA, with close positions on the Taylor diagram (Figure 5a,c), and GSMaP appears as the most effective SPE, with closer position to the reference dot on the Taylor diagram for all considered temporal scales (inter-seasonal, wet, and dry seasons).
Different to the daily quantitative statistical results, whose performance was unsatisfactory, as reported in the Taylor diagram (Figure 5a,c), all SPEs present a good ability to detect daily precipitation events (Figure 5b,d).The statistical scores are represented in the form of a performance diagram, in which geometric relationships between POD, SR, B, and CSI facilitate the inter-comparisons between the SPEs.The closer the SPEs are to the right-hand top corner black dot, the higher is the SPEs' potential to detect daily precipitation events.Overall, IMERG presented the best statistical scores according to the performance diagram for the Mean Regional dataset assessment (Figure 5b), while GSMaP presented the best results considering the Mean Spatial dataset (Figure 5d).
The ability of SPEs for detecting precipitation events is seasonally dependent.Indeed, for both datasets (Mean Regional and Mean Spatial), the best scores in the performance diagram are reached during the wet period, whereas the lowest scores are reached during the dry season.On the other hand, intermediate scores are reached for the inter-seasonal scale.Interestingly, aggregating IMERG and GSMaP from the 0.1 • grid to the 0.25 • grid yielded better results considering the Mean Spatial dataset, whereas no significant changes were observed considering the Mean Regional dataset.Comparing IMERG and GSMaP with TMPA for the 0.25 • grid shows an improvement in the ability of SPEs to detect precipitation events from the TRMM-based SPEs to the GPM-based algorithms.Indeed, IMERG and GSMaP are systematically located closer to the reference dot on the performance diagram than TMPA.

Discussion
In the present study, the potential of the last-released GPM-based SPEs, IMERG and GSMaP, was compared to the widely used TRMM-based TMPA SPEs.The analysis was performed considering different time steps (annual, monthly, and daily) and considering the potential seasonal effect over their precipitation estimates.
For the monthly time step, the consideration of fixed quality-threshold values highlighted the ability of SPEs to correctly capture monthly regional precipitation pattern in space and time.However, high CC values with unsatisfactory %B can occur independently of the considered SPEs, which is in line with previous studies [24,26,59].According to our results, a bias correction of approximately 10% is recommended before using GSMaP to follow monthly precipitation amounts.The negative bias of GSMaP supports the results of a previous study over Bolivia where GSMaP-v6 presented a %B of approximately −22.4% [26].The lowest %B value in the present study (−10%) may be the result of the improvement of the gauge calibration applied to produce the considered GSMaP-v7.
However, the inter-comparison exercise between SPEs conducted over distinct regions remains challenging and should be considered with caution.Indeed, due to the difference in spatial scale between areal (SPEs grid-cell) and point (gauges) measurement, potential conclusions using rain gauges as reference data are influenced by the gauges' locations and distribution density [26,60].Precipitation events estimated by the SPEs may not be detected by the gauges, as it might rain in other location within the grid-cell area [26,55].Therefore, part of the variability in the ability of SPEs to predict precipitation from one study to another could be attributed to the distribution density and location of the gauges used as references.In this study, data from the 56 rain gauges were aggregated into twelve 0.25 • grid-cells (mean of 4.7 gauges per 0.25 • grid-cell) with one grid-cell using up to 10 stations to calculate the reference average precipitation, whereas the 56 gauges were aggregated into forty 0.1 • grid-cells joining a maximum of four gauges within a grid-cell (mean of 1.4 gauges per 0.1 • grid-cell) (Figure 1d,e).In this context, the gauges managed to better capture the precipitation variability for the coarser grid (0.25 • ) than for the finer grid (0.1 • ).Therefore, the aggregation from 0.1 • to 0.25 • slightly enhanced the quality of SPEs from IMERG and GSMaP for both monthly and daily time steps, confirming previous results [26].To provide more insight on the conclusion reliability regarding SPEs' sensitivity to the rain gauges density, Figure 6 shows the statistical score (RMSE normalized) obtained in function of the gauge number used as reference to compare with the SPEs at the 0.25 • grid-cell scale.For both daily and monthly time step, the SPEs accuracy increase with the number of gauges used as reference (Figure 6a,b).For example, at the daily time step, highest (1) and lowest (0) RMSE values are obtained for grid-cells counting with 1 and 10 gauges, respectively, for both IMERG and GSMaP (Figure 6a).Similarly, at the monthly time scale and for all SPEs, higher normalized RMSE is obtained when only one gauge is used as reference (Figure 6b).Therefore, different conclusion on SPEs' reliability should be drawn when using more or less dense rain gauges network as reference.
It is worth mentioning that part of the variability observed in Figure 6 could also be attributed (i) to the gauges distribution into the considered grid-cells and (ii) to the consideration of different grid-cells for the different gauge number threshold (Figure 1).Therefore, Figure 6 illustrates the potential influence of gauge density into the SPEs' reliability conclusion and more studies should be developed to clearly state on gauge density influence.
To avoid any influence of the gauge density and/or distribution on the conclusion regarding the SPEs' reliability, the SPEs could be used as forcing data in hydrological modeling.Indeed, comparing simulated streamflow using SPEs as forcing data with observed streamflow aims to enhance the quality of spatial representativeness between the assessed SPEs and the reference used [34].It is worth mentioning that part of the variability observed in Figure 6 could also be attributed (i) to the gauges distribution into the considered grid-cells and (ii) to the consideration of different grid-cells for the different gauge number threshold (Figure 1).Therefore, Figure 6 illustrates the potential influence of gauge density into the SPEs' reliability conclusion and more studies should be developed to clearly state on gauge density influence.
To avoid any influence of the gauge density and/or distribution on the conclusion regarding the SPEs' reliability, the SPEs could be used as forcing data in hydrological modeling.Indeed, comparing simulated streamflow using SPEs as forcing data with observed streamflow aims to enhance the quality of spatial representativeness between the assessed SPEs and the reference used [34].
Additionally, the different dataset, metrics, and time scale used from one study to another complicate the inter-comparison exercise between SPEs conducted over distinct regions.Actually, as shown in the present study, the conclusions regarding the potential of SPEs differ according to the considered dataset (Mean Regional or Mean Spatial).For both the monthly and daily time steps, the SPEs appear to be globally more accurate when considering the Mean Regional dataset compared to the Mean Spatial dataset.Secondly, as observed for the daily time step, the conclusions regarding the potential of SPEs vary according to the considered metrics: quantitative (Taylor diagram) or categorical (performance diagram).SPEs were able to correctly detect daily precipitation events (categorical metrics) but failed to retrieve daily precipitation amounts (quantitative metrics).Finally, the consideration of different time scales (inter-seasonal, wet, and dry season) reveals a strong seasonal influence on the potential of SPEs, with SPEs performing better during the wet season than during the dry season for both monthly and daily time steps.This finding is in line with previous studies reporting on the influence of season on the potential of SPEs [24,26,58,61].Evaluated independently, the time scale may influence the conclusion regarding the potential of the SPEs.For the above-detailed reasons, the use of variable datasets, metrics, and time scales should introduce important discrepancy on the conclusion regarding the potential of SPEs obtained from one study to another.Therefore, studies on the potential of SPEs to predict precipitation for a determined region should not be used to determine the performance of SPEs for another region, and the assessment of the potential of SPEs proves to be an essential step before their use in impact studies.
Finally, according to the presented results, the SPEs appear to be unsuitable for following the monthly precipitation dynamic during the dry season.However, this observation has to be weighted.Indeed, the mean monthly precipitation values for the dry season are very low (Figure 3).Therefore, a minimum error in precipitation estimates will automatically induce very large %RMSE and %B.Actually, even if large %RMSE and %B are observed for the dry season, these errors are insignificant regarding the total annual precipitation input over the Brazilian Central Plateau, where the dry months contribute less than 15% of the annual total precipitation.The relatively high and satisfactory CC values (>0.7) obtained by all the considered SPEs mean that the SPEs correctly capture the Additionally, the different dataset, metrics, and time scale used from one study to another complicate the inter-comparison exercise between SPEs conducted over distinct regions.Actually, as shown in the present study, the conclusions regarding the potential of SPEs differ according to the considered dataset (Mean Regional or Mean Spatial).For both the monthly and daily time steps, the SPEs appear to be globally more accurate when considering the Mean Regional dataset compared to the Mean Spatial dataset.Secondly, as observed for the daily time step, the conclusions regarding the potential of SPEs vary according to the considered metrics: quantitative (Taylor diagram) or categorical (performance diagram).SPEs were able to correctly detect daily precipitation events (categorical metrics) but failed to retrieve daily precipitation amounts (quantitative metrics).Finally, the consideration of different time scales (inter-seasonal, wet, and dry season) reveals a strong seasonal influence on the potential of SPEs, with SPEs performing better during the wet season than during the dry season for both monthly and daily time steps.This finding is in line with previous studies reporting on the influence of season on the potential of SPEs [24,26,58,61].Evaluated independently, the time scale may influence the conclusion regarding the potential of the SPEs.For the above-detailed reasons, the use of variable datasets, metrics, and time scales should introduce important discrepancy on the conclusion regarding the potential of SPEs obtained from one study to another.Therefore, studies on the potential of SPEs to predict precipitation for a determined region should not be used to determine the performance of SPEs for another region, and the assessment of the potential of SPEs proves to be an essential step before their use in impact studies.
Finally, according to the presented results, the SPEs appear to be unsuitable for following the monthly precipitation dynamic during the dry season.However, this observation has to be weighted.Indeed, the mean monthly precipitation values for the dry season are very low (Figure 3).Therefore, a minimum error in precipitation estimates will automatically induce very large %RMSE and %B.Actually, even if large %RMSE and %B are observed for the dry season, these errors are insignificant regarding the total annual precipitation input over the Brazilian Central Plateau, where the dry months contribute less than 15% of the annual total precipitation.The relatively high and satisfactory CC values (>0.7) obtained by all the considered SPEs mean that the SPEs correctly capture the monthly precipitation dynamics, even during the dry season.Interestingly, during the dry season, higher distribution is observed in both the Taylor and performance diagrams, meaning that SPEs' sensitivity to low precipitation amounts or low number of precipitation events (dry season) differs more than to high precipitation amounts or large number of precipitation events (wet season).According to this study, GSMaP is the most adapted SPE to follow daily and monthly precipitation dynamics during the dry seasons over the Brazil Central Plateau region.However, as discussed before, a bias correction is strongly recommended.

Conclusions
This study evaluated the performance of the new GPM-based SPEs in the Brazilian Central Plateau region, using data from precipitation gauges as reference.The assessment compared IMERG and GSMaP at their original 0.1 • spatial resolution and used an adjusted/refitted 0.25 • grid to compare them with TMPA on an annual, monthly, and daily basis.The major conclusions are:

•
Generally, GPM-based SPEs products (IMERG, GSMaP) are able to ensure precipitation monitoring with similar or even better accuracy than that obtained previously with TRMM-based TMPA products in the Brazilian Central Plateau.

•
IMERG presents the best annual and monthly results in nearly all metrics used in the analysis for all considered temporal scales (inter-seasonal, wet, and dry season).

•
GSMaP precipitation estimations presented negative bias for monthly and annual precipitation amounts for all considered time scales (inter-seasonal, wet, and dry season).A bias correction of approximately 10% is recommended for monthly time step.

•
For the daily time step, all SPEs correctly detected precipitation events but considerably failed in the quantification of daily precipitation amount.Among the considered SPEs, GSMaP presented the highest ability for detecting and quantification of daily precipitation events.

•
Despite the differences between IMERG and GSMaP, their results indicate an improvement in detecting precipitation events from TRMM-based SPEs to GPM-based SPEs.

•
The potential of all SPEs presented a strong seasonal variability.For both the daily and monthly time steps, all SPEs performed better during the wet season than during the dry season.The new GPM-based SPEs, IMERG and GSMaP, still have difficulty in estimating precipitation in the dry period, when precipitation events are generally less intense, of lower volume, and more sparsely distributed across the territory.Among the considered SPEs, GSMaP presents the highest ability for monthly and daily precipitation monitoring during the dry season.
More generally, this study also shows that the SPEs' reliability conclusion is sensitive to the density of the rain gauges used as reference.Indeed, more optimist/pessimist conclusions on SPEs' reliability should be drawn if considering denser/scarcer gauge network.Similarly, SPEs' reliability conclusions depend on the considered indicator.Different conclusions on SPEs' reliability can be drawn according to the considered assessment.In this line, future study should focus on the impacts of modeling sensitivity to various SPEs to assess their reliability in a more applicative context.However, for the time being, only four years of data are available for IMERG and GSMaP, which compromises robust calibration/validation steps for impact models.

Figure 1 .
Figure 1.Study area: location within four important hydrographic regions (a); rain gauges and elevation map (b); average monthly precipitation between 1963 and 2016 for Brasilia's rain gauge (c); number of gauges by grid-cell for the 0.1° grid (d); and number of gauges by grid-cell for the 0.25° grid (e).

Figure 1 .
Figure 1.Study area: location within four important hydrographic regions (a); rain gauges and elevation map (b); average monthly precipitation between 1963 and 2016 for Brasilia's rain gauge (c); number of gauges by grid-cell for the 0.1 • grid (d); and number of gauges by grid-cell for the 0.25 • grid (e).

Figure 2 .
Figure 2. 2015 mean spatial annual precipitation registered by satellite precipitation estimates (SPEs) and gauges (a) and accuracy expressed in statistical metrics; (b) summarized in the form of a Taylor diagram with values normalized by the reference standard deviation (STD); (c).IMERG: Integrated Multi-satellitE Retrievals for Global Precipitation Measurement; GSMaP: Global Satellite Mapping of Precipitation; TMPA: Tropical Rainfall Measurement Mission (TRMM) Multi-satellite Precipitation Analysis; RMSE: root-mean-square error.

Figure 2 .
Figure 2. 2015 mean spatial annual precipitation registered by satellite precipitation estimates (SPEs) and gauges (a) and accuracy expressed in statistical metrics; (b) summarized in the form of a Taylor diagram with values normalized by the reference standard deviation (STD); (c).IMERG: Integrated Multi-satellitE Retrievals for Global Precipitation Measurement; GSMaP: Global Satellite Mapping of Precipitation; TMPA: Tropical Rainfall Measurement Mission (TRMM) Multi-satellite Precipitation Analysis; RMSE: root-mean-square error.Despite the SPEs respective ability to represent the regional precipitation pattern, Figure2allows an inter-comparison of the potential of SPEs to estimate the annual precipitation amount.When comparing the GPM-based SPEs with the TRMM-based SPEs, both IMERG and GSMaP are more able to quantify the annual precipitation amount.Indeed, in the Taylor diagram (Figure2c),

Water 2019 ,
11, x FOR PEER REVIEW 10 of 18

Figure 6 .
Figure 6.Gauges number influence onto the SPEs accuracy expressed in term of normalized RMSE between 0 and 1.The analysis is applied for the 0.25° grid-cells scale at both daily (a) and monthly; (b) time step.Dotted line shows the linear trends.

Figure 6 .
Figure 6.Gauges number influence onto the SPEs accuracy expressed in term of normalized RMSE between 0 and 1.The analysis is applied for the 0.25 • grid-cells scale at both daily (a) and monthly; (b) time step.Dotted line shows the linear trends.