Inter-Comparison and Evaluation of Remote Sensing Precipitation Products over China from 2005 to 2013

Precipitation is a key aspect of the climate system. In this paper, the dependability of five satellite precipitation products (TRMM [Tropical Rainfall Measuring Mission] 3BV42, PERSIANN [Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks] CDR, GSMaP [Global Satellite Mapping of Precipitation] RENALYSIS, CMORPH [Climate Prediction Center’s morphing technique] BLD and CMORPH_RAW) were compared with in situ measurements over China for the period of 2005 to 2013. To completely evaluate these precipitation products, the annual, seasonal and monthly precipitation averages were calculated. Overall, the Huaihe River and Qinlin mountains are shown to have heavy precipitation to the southeast and lighter precipitation to the northwest. The comparison results indicate that Gauge correction (CMORPH_BLD) improves the quality of the original satellite products (CMORPH_RAW), resulting in the higher correlation coefficient (CC), the low relative bias (BIAS) and root mean square error (RMSE). Over China, the GSMaP_RENALYSIS outperforms other products and shows the highest CC (0.91) and lowest RMSE (0.85 mm/day) and all products except for PERSIANN_CDR exhibit underestimation. GSMaP_RENALYSIS gives the highest of probability of detection (81%), critical success index (63%) and lowest false alarm ratio (36%) while TRMM3BV42 gives the highest of frequency bias index (1.00). Over Tibetan Plateau, CMORPH_RAW demonstrates the poorest performance with the biggest BIAS (4.2 mm/month) and lowest CC (0.22) in December 2013. GSMaP_RENALYSIS displays quite consistent with in situ measurements in summer. However, GSMaP_RENALYSIS and CMORPH_RAW underestimate precipitation over South China. CMORPH_BLD and TRMM3BV42 show consistent with high CC (>0.8) but relatively large RMSE in summer.


Introduction
Precipitation is one of important factors in the global water cycle and plays a key role in the global energy system.Accurate precipitation data is helpful for short-term and long-term weather forecasts, which allow proactive measures to prevent natural hazards, such as floods, landslide and debris flow, using rainfall models.Measurements of precipitation can be made using rain gauge stations, ground radar, or remote sensing technology.Rain gauge stations are traditional point measurements and have higher accuracy but it is difficult to obtain such precipitation measurements over oceanic and high-altitude areas where only sparse rain gauge stations are available.Ground radar is considered to more accurate for retrieving regional rainfall estimates but ground radar is aimed at monitoring extreme events over limited time spans [1].Additionally, ground radar is limited over relatively complex terrains, where radar signals maybe be affect by beam blockages and vertical reflectivity profiles [2].The development of satellite remote sensing has provided several advantages, including its broad spatial coverage, automatic data acquisition and long-term continuity but the signals themselves are affected by the atmospheric and the surface, which will lead to retrieval uncertainties [3].However, remote sensing is an important method for retrieving high spatial and temporal resolution rainfall measurements over complex terrains and mountainous areas, where both rain gauge stations and ground radar are very limited or unavailable.The satellite retrieval methodologies fall primarily into three categories, which there are the visible and infrared (VIS/IR) methods, microwave (MW) methods and multi-sensor methods.
The VIS precipitation retrieval method mainly utilizes the empirical relationship between cloud brightness and rainfall because precipitation is considered to occur under highly bright clouds, which can be associated with cloud types, cloud areas and cloud growth rates [4][5][6].IR imagery is associated with cloud top temperatures and cloud growth rates can be obtained via the thermal emissions during both night and day; generally, heavier rainfall tends to be associated with larger, taller clouds with colder cloud tops [5].Many researchers have developed different methods to retrieve precipitation based on VIR/IR data from geosynchronous earth orbit (GEO) and low-earth orbit (LEO) satellites, including the 3-hourly and monthly mean rainfalls [6][7][8][9][10].MW sensors can detect rain clouds directly and can provide information about the atmospheric constituents and hydro-meteorological profiles, which are more directly related to the ground precipitation rate [11,12].Many previous works have developed a number of MW methods to retrieve precipitation [13][14][15].However, precipitation is not always associated with clouds and the VIS/IR methods cannot detect information below clouds [16].Meanwhile, MW instruments are restricted to polar-orbiting platforms, meaning that they obtain only a small amount of data and have the disadvantages of low spatial and temporal resolutions [17].The idea of combining VIS/IR and MW observations to retrieve precipitation from multiple satellite sensors was proposed to compensate for the weaknesses of the two methods [18].Many products were developed by blending the MW and IR data, such as CMORPH (Climate Prediction Center's morphing technique), PERSIANN (Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks), TRMM (Tropical Rainfall Measuring Mission), GSMaP (Global Satellite Mapping of Precipitation).These precipitation products can help climate and hydrological research.Many researchers have compared different satellite precipitation products with in situ measurements over different areas and have proved that satellite estimates showed non-negligible biases.For example, Hirpa et al. evaluated the accuracies of three products (CMORPH, PERSIANN and TRMM3BV42) and found that both CMORPH and TRMM underestimated precipitation over higher elevations [19].Guo Hao et al. evaluated four products (TRMM, CMORPH, PERSIANN, GSMaP) using in situ measurements over Central Asia from 2004 to 2006 and most of the products overestimated the precipitation [20].Awange et al. used a "three-cornered-hat" method to assess six precipitation products and indicated that the RG-merged products had higher accuracies than the satellite-only products [21].It is necessary to discuss the performances of satellite precipitation products with the aim of determining whether a certain product is appropriate for a specific region.However, there are few researches among multi-satellite precipitation products evaluation based on long time series over China.
In this study, considering the study region and the availability of period for all satellite precipitation products, we will mainly assess the strengths and weaknesses of a few widely used rainfall products (CMORPH_CDR, CMORPH_RAW, TRMM3BV42, GSMaP_gauge_renalysis and PERSIANN_CDR) over China from 2005 to 2013.The spatial distributions, temporal variations and rainfall intensities were analyzed and compared with in situ measurements.Precipitation has a particular spatiotemporal distribution over China and thus, we discuss the precipitation over the Tibetan Plateau and South China in 2013.This paper aims to lead the reader to select more accurate satellite precipitation products over China and this research can help to improve the corresponding precipitation retrieval algorithms.This paper is structured as follows: the first part is an introduction.The rainfall products will be present in Section 2. Section 3 expounds on an integrated analytical approach and debates the performances of the space-time variations and errors of different datasets.Finally, the comparative result will be discussed.

Satellite-Based Precipitation Products
The TRMM satellite was launched in 1997 and was first used to measure precipitation.The National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) developed the Multisatellite Precipitation analysis (TMPA) algorithm.The TMPA algorithm combined multiple pass microwave (PMW) datasets, including the Microwave Imager (TMI), Special Sensor Microwave Imager (SSM/I), Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E), Advanced Microwave sounding Unit-B (AMSU-B) and IR data being obtained by the international constellation of GEOs satellites [22].The algorithm first calibrated the PMW data, which were used to create infrared precipitation data.Then, both datasets were combined and the rain gauge data were incorporated.TRMM has different spatial and temporal resolutions to satisfy the demands of its users.The real-time and post-real-time versions have the same spatial resolutions of 0.25 • but the data in the post-real-time version is not corrected by the rain gauge data and has a relatively narrow coverage area.The TRMM products have three different temporal resolutions, including three hourly, daily and monthly.This study obtained a TRMM3BV42 with a daily/0.25 • resolution, which was gauge-corrected by the Global Precipitation Climatology Center (GPCC) and the Climate Assessment and Monitoring System (CAMS).
CMORPH is produced based on PMW and IR measurements from the National Oceanic and Atmospheric Administration (NOAA)/Climate Prediction Center (CPC) and has a relatively high spatiotemporal resolution [11,23].IR images (including Meteosat-5/7, Geostationary Meteorological Satellite-5 [GMS-5] and Geostationary Operational Environmental Satellite [GOES-8/10]) are obtained every 30 min at large zenith angles [23] and PMW (including TMI, SSM/I and AMSU-B) can be used to retrieve rain rate relations to the IR images.Joyce et al. exploited IR data to deduce motion fields and spread the PMW rain domains across space and time via a time-weighting interpolation algorithm [11].There are three different spatial and temporal resolution CMORPH products: 8 km-30 min, 0.25 • -3 hourly and 0.25 • -daily.These products can be freely downloaded from NOAA for all dates since 2002.This paper employed two types of CMORPH data (CMORPH_RAW and CMORPH_BLD) to analyze the precipitation characteristics from 2005 to 2013 with a spatiotemporal resolution of 0.25 • /h.CMORPH_BLD indicates that the satellite data have been calibrated by CPC unified daily gauge analysis methods over land [24].
PERSIANN is a pure satellite precipitation product, which is produced based on an artificial neural network (ANN) model by the Center for Hydrometeorology and Remote Sensing (CHRS) of the University of California [25,26].The PERSIANN algorithm combines IR images of the Geostationary Environmental satellite and TMI on TRMM.The algorithm mainly depends on the statistical relationship between IR and the precipitation rate and the retrieved precipitation displays considerable uncertainty [26].Therefore, when PMW data are available, studies adopt a training technique to revise the PERSIANN products with the PMW data.PERSIANN products have three different temporal resolutions (3 h, 6 h and daily) and a spatial resolution of 0.25 • and the all products can be freely downloaded for dates since 2000.This paper employed the PERSIANN_CDR products with a spatiotemporal resolution of 0.25 • /daily; these are calibrated by the Global Precipitation Climatology Project (GPCP) monthly gauge analysis [27].
GSMaP is a multi-satellite precipitation data source developed by the Japan Science and Technology Agency (JST) and Japan Aerospace Exploration Agency (JAXA) in November 2002 [28,29].
The algorithms adopt a Kalman filtering technique to estimate the hourly global precipitation, wherein the highest spatiotemporal resolution of GSMaP is 30 min/0.0365• over the equator and 0.1 • -daily other area, as well as merged PMW (including TMI, AMSR-E and SSM/I) data and IR images from GEO satellite [28].GSMaP is researched to achieve three mainly goals: first, the project uses the current PMW data to produce high spatiotemporal resolution global precipitation maps.Second, the algorithm of the physical model of the precipitation and rainfall rate can constantly be improved.Third, evaluating the accuracy of the precipitation products can aid the Global Precipitation Measurement (GPM) projection, which launched in 2014.The GSMaP team has achieved different standard products according to diverse algorithms, including the GSMaP_NRT, GSMaP_MVK, GSMaP_Gauge and GSMaP_RENALYSIS.The study selected the GSMaP_RENALYSIS merged with the CPC global rain gauge data to analyze the precipitation over China after considering the features of the study period [30].Information concerning all the products is listed in Table 1.

Gauge Stations
This study adopted 830 rain gauge stations from the China Meteorological Data Sharing Service System (http://data.cma.cn/) and nine years of daily data were obtained from 2005 to 2013 to assess the distribution of precipitation over a long time.Meanwhile, we used data from 2280 rain gauge stations in 2013 to completely analyze the accuracies of the satellite precipitation products over the Tibetan Plateau and South China areas.Rain gauge stations are sparsely distributed across the Tibetan Plateau because of its extreme terrain and climatic conditions and some rain gauge stations are closed in the winter to protect their instruments, so some data is unavailable.There are more stations over the Yangtze River, South China and North China, which are used to accurately analyze the distributions of the rainfall.The western Northwest region also has relatively few stations.The elevation data of the study area were obtained from Shuttle Radar Topographic Mission (SRTM on http://srtm.csi.cgiar.org/), of which the spatial resolution is about 90 m.The rain gauge station distributions and the topography of China are shown Figure 1.
Remote Sens. 2018, 10, x FOR PEER REVIEW 4 of 18 and IR images from GEO satellite [28].GSMaP is researched to achieve three mainly goals: first, the project uses the current PMW data to produce high spatiotemporal resolution global precipitation maps.Second, the algorithm of the physical model of the precipitation and rainfall rate can constantly be improved.Third, evaluating the accuracy of the precipitation products can aid the Global Precipitation Measurement (GPM) projection, which launched in 2014.The GSMaP team has achieved different standard products according to diverse algorithms, including the GSMaP_NRT, GSMaP_MVK, GSMaP_Gauge and GSMaP_RENALYSIS.The study selected the GSMaP_RENALYSIS merged with the CPC global rain gauge data to analyze the precipitation over China after considering the features of the study period [30].Information concerning all the products is listed in Table 1.

Gauge Stations
This study adopted 830 rain gauge stations from the China Meteorological Data Sharing Service System (http://data.cma.cn/) and nine years of daily data were obtained from 2005 to 2013 to assess the distribution of precipitation over a long time.Meanwhile, we used data from 2280 rain gauge stations in 2013 to completely analyze the accuracies of the satellite precipitation products over the Tibetan Plateau and South China areas.Rain gauge stations are sparsely distributed across the Tibetan Plateau because of its extreme terrain and climatic conditions and some rain gauge stations are closed in the winter to protect their instruments, so some data is unavailable.There are more stations over the Yangtze River, South China and North China, which are used to accurately analyze the distributions of the rainfall.The western Northwest region also has relatively few stations.The elevation data of the study area were obtained from Shuttle Radar Topographic Mission (SRTM on http://srtm.csi.cgiar.org/), of which the spatial resolution is about 90 m.The rain gauge station distributions and the topography of China are shown Figure 1.

Statistical Evaluation Methods
In this study, five products are compared with the data from rain gauge stations, which are often viewed as true values.We will discuss the characteristics of the precipitation on the yearly, seasonally, monthly and daily time scales.We use statistical analysis methods, including the linear correlation coefficient (CC), relative bias (BIAS), root mean square error (RMSE), the frequency bias index (FBI), probability of detection (POD), false alarm ratio (FAR) and critical success index (CSI) to compare the differences from pixel to point [31,32].There is a high correlation between the satellite products and in situ measurements if the CC is greater than 0.7 [33].The result is considered underestimation (overestimation) of the precipitation when the BIAS is less than (more than) zero [34].POD and FAR express the ratios of the rainfall events, which are accurately measured and provide error warnings.The CSI denotes the overall ratio of the rainfall events that are accurately estimated by the satellite.The best comparison result occurs when the FBI, POD and CSI are equal to one and when FAR equals zero [35].All indices are defined as follows.
where S represents the satellite-retrieved precipitation and G represents the in situ-measured precipitation.S and G indicate the mean satellite precipitation and in situ-measured precipitation, respectively.N is the number of rainfall pairs between the gauge and satellite data.σ s and σ v are the standard deviations of the satellite and gauge.H, F and M denote the amounts of hits, false alarms and misses, which are calculated as described in Table 2 [32].The threshold was set 0.1 mm/day considering the minimum of the situ measurements.

Results and Discussion
4.1.Evaluation on the Regional Scale

Nine-Year Daily Mean Precipitation
China has special topographic and geomorphological features and its precipitation shows different distributions.Overall, all satellite-based products and in situ measurements have the same trends, wherein the precipitation distribution decreases from the southeast coast to the northwest inland.Figure 2 shows the nine-year daily mean precipitation of 2005 to 2013.The southeast has a higher mean daily precipitation than that in other areas, which some areas are as high as to 12 mm/day.The northwest has a typical semi-arid and arid climate with little precipitation, which most of areas are less than 2 mm/day.The phenomenon is mainly determined by climate and topography.CMORPH_BLD, CMORPH_RAW and TRMM3BV42 have clearly underestimated precipitation compared with the in-situ measurements in the south.This can be explained that there is more shallow and warm rains in the south, while the microwave-based algorithms can more easily obtain data for heavy rainfall events than weak rainfall events [36,37].PERSIANN-CDR show a remarkable overestimation of precipitation in most areas, such as southern Tibet and regions of South China (such as Guangzhou and Fujian), wherein PERSIANN mainly adopts the IR brightness temperature data from geostationary satellites to estimate rainfall rates.MW is used only to adjust the model parameters, which differs from other satellite retrieval algorithms [25,27].CMORPH_BLD, CMORPH_RAW, TRMM3BV42 and PERSIANN underestimated in north China, in where the precipitation mainly occurs in July and August.The precipitation is characterized by strong rainfall but short duration, so it is possible that sensors missed precipitation during the period.GSMaP_RENALYSIS is most consistent with the in-situ measurements but also overestimated values in some regions in the south (such as the south of Sichuan, Guangdong, Jiangxi) and northeast of China.All satellite products showed overestimations over the south of Tibet, in where are no-rain gauge stations (Figure 1).The precipitation was obtained by interpolation technique based on surrounding stations located semi-humid region.The south of Tibet belongs to humid region and has heavy precipitation.Therefore, the precipitation by interpolation is less than satellite retrieval.However, all satellite products underestimated precipitation over the most of arid region, the passive microwave retrieval of precipitation might be hampered by the snow and ice surface [38,39].The satellite bias-corrected product (CMORPH_BLD) is more consistent with the in-situ measurements than the original satellite product (CMORPH_RAW) according to the precipitation distribution.Therefore, the bias-correction using the in-situ measurements greatly improved the product accuracy.
The density-colored scatter plots are shown in Figure 3, which provides a quantitative comparison among the satellite products and in situ measurements for the nine-year mean daily precipitation values over China (including 830 stations).The bias-corrected products, except for that of PERSIANN_CDR, have greater R and lower RMSE than CMORPH_RAW.The R of CMORPH_BLD is 0.8 higher than CMORPH_RAW (0.73) and the slope of CMORPH_BLD and CMORPH_RAW are 0.58 and 0.32, respectively.The result suggest that the bias-corrected products have improved the accuracy greatly (Figure 3a,b).TRMM3BV42 and CMORPH_BLD have similar distributions of density scatter and the differences of R and RMSE are small (Figure 3b,c).The PERSIANN_CDR has a high R (0.78) but also yields high RMSE (2.93 mm/day) and the discrete degree of point is larger than other products, which PERSIANN uses only IR data to estimate rainfall rate (Figure 3d).The highest R and lowest RMSE of GSMaP_RENALYSIS are 0.91 and 0.85 mm/day and these points are uniform distributed, because GSMaP_RENALYSIS has inherited CMORPH's morphing algorithm and employs a new Kalman filter approach to assimilate IR-derived rain rates, which it can help to reduce the total errors (Figure 3e) [11,39,40].
To analyze the continuous variations, the CC, RMSE and BIAS are calculated for each year from 2005 to 2013, as is shown Table 3. GSMaP_RENALYSIS has the highest CC (more than 0.8) and the lowest RMSE and BIAS (about ±0.1 mm/day).GSMaP_RENALYSIS underestimated the precipitation in both 2005 and 2006 according to its BIAS, which is less than zero but it slightly overestimated the precipitation during 2007 to 2013.The PERSIANN_CDR overestimated the precipitation in each year and CMORPH_BLD, CMORPH_RAW and TRMM3BV42 underestimated the precipitation, which is similar to the nine-year daily mean precipitation.Note that all products have greater CC in 2007 than in 2008 but the RMSE are also higher, which is attributed to the occurrence of the decadal maximum of precipitation, such that the national mean precipitation was approximately 654.8 mm.In 2012, the CC of CMORPH_RAW and CMORPH_BLD were 0.37 and 0.76, respectively, which shows that the gauge correction greatly improved the precipitation estimation accuracies.Overall, GSMaP_RENALYSIS's performance is comparable to other satellite-based products over China.
Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 18 lowest RMSE and BIAS (about ±0.1 mm/day).GSMaP_RENALYSIS underestimated the precipitation in both 2005 and 2006 according to its BIAS, which is less than zero but it slightly overestimated the precipitation during 2007 to 2013.The PERSIANN_CDR overestimated the precipitation in each year and CMORPH_BLD, CMORPH_RAW and TRMM3BV42 underestimated the precipitation, which is similar to the nine-year daily mean precipitation.Note that all products have greater CC in 2007 than in 2008 but the RMSE are also higher, which is attributed to the occurrence of the decadal maximum of precipitation, such that the national mean precipitation was approximately 654.8 mm.In 2012, the CC of CMORPH_RAW and CMORPH_BLD were 0.37 and 0.76, respectively, which shows that the gauge correction greatly improved the precipitation estimation accuracies.Overall, GSMaP_RENALYSIS's performance is comparable to other satellite-based products over China.lowest RMSE and BIAS (about ±0.1 mm/day).GSMaP_RENALYSIS underestimated the precipitation in both 2005 and 2006 according to its BIAS, which is less than zero but it slightly overestimated the precipitation during 2007 to 2013.The PERSIANN_CDR overestimated the precipitation in each year and CMORPH_BLD, CMORPH_RAW and TRMM3BV42 underestimated the precipitation, which is similar to the nine-year daily mean precipitation.Note that all products have greater CC in 2007 than in 2008 but the RMSE are also higher, which is attributed to the occurrence of the decadal maximum of precipitation, such that the national mean precipitation was approximately 654.8 mm.In 2012, the CC of CMORPH_RAW and CMORPH_BLD were 0.37 and 0.76, respectively, which shows that the gauge correction greatly improved the precipitation estimation accuracies.Overall, GSMaP_RENALYSIS's performance is comparable to other satellite-based products over China.Figure 4 shows the distribution of the seasonal daily mean precipitation values interpolated from the rain gauge stations and the multiple satellite-based precipitation products over China.The CC, BIAS and RMSE statistical values are plotted in Figure 5.
As shown in the seasonal precipitation distribution interpolated from the rain gauge stations (GAUGE_STATION), the summer and spring experience more precipitation (the maximum amounts are 15 mm/day and 10 mm/day, respectively) because both are part of the "rainy" season in China influenced by the southwest monsoon and southeast monsoon.The winter and autumn experienced less precipitation (with maxima of 5 mm/day and 8 mm/day, respectively).This phenomenon is related to the passage of rain over China.It should be noted that there have more precipitation in summer and spring than in autumn and winter because the rainfall events begin in southern China from February to May every year.The precipitation then moves north with the heat of the seasons.Rainfall begins to fall rapidly when the cold air comes from the north at the end of August or the beginning of September.There is greater rainfall in the summer than in the autumn over the south because the period of the "plum rains" season occurs from April to June in China.The northwest experiences little precipitation in all four seasons as it belongs to the arid and semi-arid region of China, in where is far from the sea and is blocked by mountains.
Generally, the distributions of the satellite precipitation products are similar to that of the rain gauge stations, except for CMORPH_RAW, which is not bias corrected.The GSMaP_RENALYSIS are the most similar to the rain gauge data (CC more than 0.8 and BIAS less than ±0.1 mm/day, respectively).All satellite products underestimated the precipitation in the winter over some regions, such as northeast of Neimeng and northwest of Tibet.This might be attributed to the following reasons: The little precipitation in this area is mainly produced from the Atlantic Ocean through the western wind belt and it is difficult to capture low-level precipitation by PMW sensor.Meanwhile, it might be caused the more snow and ice cover on the ground in winter.All products overestimated the precipitation in the summer over the south of Tibet, in where have heavy precipitation.The phenomenon might be related to the sparsely and unevenly rain gauges.The interpolation technique is used to spread nonzero precipitation into zero and lead to obtain less precipitation.CMORPH_RAW is poorly consistent with all other products and with the in-situ measurements over all four seasons.CMORPH_RAW and CMORPH_BLD show some abnormal zero precipitation over Tibetan Plateau in winter.It is possible that the ice and snow surface tend to influence PWM-based accurate retrieval over land [37].PERSIANN_CDR exhibited serious overestimations of precipitation over southern China.This result may be attributed to the lack of training of the ANN parameters over China because PERSIANN_CDR is only adequately trained over the United States [41].We calculated the quality analysis of the satellite products and in situ measurements, as shown in Figure 5.All products provide the highest CC and lowest RMSE (greater than 0.8 and less than 1 mm/day, respectively) in the winter but the highest BIAS also occurs in the winter.This might be related to the following reasons: the rainfall events are less than other seasons due to climate and the rainfall intensity is usually small in winter.Next, the ice and snow surface might to influence satellite retrievals over land and the gauge stations might measure abnormal data over the high-latitude areas and high-altitude areas, which the data were considered abnormal data during dealing.This should be partially responsible for the high CC, BIAS and low RMSE in winter.There are more strong convective precipitation events well detected by microwave-base algorithms during the summer [36,37].However, the satellite might miss some precipitation events because the duration of precipitation is short, which may cause some errors.All products show the lowest BIAS (less than 0.5 mm/day) and highest RMSE (about 2 mm/day) in summer, meanwhile the CC (more than 0.7) is relative high.
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 18 than 0.8 and less than 1 mm/day, respectively) in the winter but the highest BIAS also occurs in the winter.This might be related to the following reasons: the rainfall events are less than other seasons due to climate and the rainfall intensity is usually small in winter.Next, the ice and snow surface might to influence satellite retrievals over land and the gauge stations might measure abnormal data over the high-latitude areas and high-altitude areas, which the data were considered abnormal data during dealing.This should be partially responsible for the high CC, BIAS and low RMSE in winter.
There are more strong convective precipitation events well detected by microwave-base algorithms during the summer [36,37].However, the satellite might miss some precipitation events because the duration of precipitation is short, which may cause some errors.All products show the lowest BIAS (less than 0.5 mm/day) and highest RMSE (about 2 mm/day) in summer, meanwhile the CC (more than 0.7) is relative high.

Monthly Daily Mean Precipitation
Statistical information (e.g., CC, BIAS and RMSE) is computed to study the monthly daily mean precipitation from 2005 to 2013, as shown in Figure 6.Some previous works have suggested that the biases and errors could be magnified over regions with no gauge stations [42].Therefore, this study calculated only the monthly mean daily data to guarantee more than one gauge station per pixel.Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 18 than 0.8 and less than 1 mm/day, respectively) in the winter but the highest BIAS also occurs in the winter.This might be related to the following reasons: the rainfall events are less than other seasons due to climate and the rainfall intensity is usually small in winter.Next, the ice and snow surface might to influence satellite retrievals over land and the gauge stations might measure abnormal data over the high-latitude areas and high-altitude areas, which the data were considered abnormal data during dealing.This should be partially responsible for the high CC, BIAS and low RMSE in winter.
There are more strong convective precipitation events well detected by microwave-base algorithms during the summer [36,37].However, the satellite might miss some precipitation events because the duration of precipitation is short, which may cause some errors.All products show the lowest BIAS (less than 0.5 mm/day) and highest RMSE (about 2 mm/day) in summer, meanwhile the CC (more than 0.7) is relative high.

Monthly Daily Mean Precipitation
Statistical information (e.g., CC, BIAS and RMSE) is computed to study the monthly daily mean precipitation from 2005 to 2013, as shown in Figure 6.Some previous works have suggested that the biases and errors could be magnified over regions with no gauge stations [42].Therefore, this study calculated only the monthly mean daily data to guarantee more than one gauge station per pixel.In this paragraph, the satellite precipitation products are compared with in situ measurements for detecting the precipitation events and intensities over China during the nine years from 2005 to 2013.The PDF is the proportion of the number of times that rainfall events from each bin occurs divided by the total number of rainfall events.To considering the rainfall intensity (R), this study has divided rainfall into eight bins [20]: (1) no rain (R = 0), (2) 0 < R ≤ 0.5 mm/day, (3) 0.5 mm/day < R ≤ 1 mm/day, (4) 1 mm/day < R ≤ 2 mm/day, (5) 2 mm/day < R ≤5 mm/day, (6) 5 mm/day < R ≤ 10 mm/day, (7) 10 mm/day< R ≤ 20 mm/day and (8) R > 20 mm/day.The pixels nearest the gauge stations were selected to calculate the rainfall intensities.Approximately 60% of the values computed using the rain gauge stations show no-rain events, as shown in Figures 7 and 8. CMORPH_RAW, CMORPH_BLD and TRMM3BV42 can be inclined to detect more no-rain events than the gauge stations (around 79%, 65.4% and 64.3%, respectively), which suggests that satellite measurements may miss some rain events.CMORPH_RAW detected more no-rain events than CMORPH_BLD, which showed that the bias-corrected satellite data tend to record more rainfall events.GSMaP_RENALYSIS and CMORPH_BLD record more light rainfall events (0 < R ≤ 0.5 mm/day, around 8.5% and 8.8%, respectively) than the other datasets.CMORPH_BLD, GSMaP_RENALYSIS and PERSIANN_CDR can record more rainfall events than the other products when the rainfall intensities are between 0.5 mm/day and 5 mm/day, where they may have an R value as high as 2 < R ≤ 5 mm/day (around 7.7%, 9.3% and 8.1%, respectively).These results can be explained follows: (1) the satellites can easily detect strong, convective rainfall events but are inclined to miss shallow and warm rains; (2) the bias-corrected methods can boost the amplitudes of the detected events to compensate for the missed events [37].This study analyzed FBI, POD, CSI and FAR to evaluate the qualities of the errors between the satellite products and the gauge station data.The results of the statistical analyses are shown in Table 4. CMORPH_RAW FBI is better than CMORPH_BLD, such that the POD and CSI are low and FAR is high.TRMM3BV42 FBI is the closest to the gauge data of all the products but the POD and CSI are lower and the FAR is higher than for GSMaP_RENALYSIS.Therefore, GSMaP_RENALYSIS can more accurately detect rain events than the other products.CMORPH_RAW, CMORPH_BLD and TRMM3BV42 can be inclined to detect more no-rain events than the gauge stations (around 79%, 65.4% and 64.3%, respectively), which suggests that satellite measurements may miss some rain events.CMORPH_RAW detected more no-rain events than CMORPH_BLD, which showed that the bias-corrected satellite data tend to record more rainfall events.GSMaP_RENALYSIS and CMORPH_BLD record more light rainfall events (0 < R ≤ 0.5 mm/day, around 8.5% and 8.8%, respectively) than the other datasets.CMORPH_BLD, GSMaP_RENALYSIS and PERSIANN_CDR can record more rainfall events than the other products when the rainfall intensities are between 0.5 mm/day and 5 mm/day, where they may have an R value as high as 2 < R ≤ 5 mm/day (around 7.7%, 9.3% and 8.1%, respectively).These results can be explained follows: (1) the satellites can easily detect strong, convective rainfall events but are inclined to miss shallow and warm rains; (2) the bias-corrected methods can boost the amplitudes of the detected events to compensate for the missed events [37].This study analyzed FBI, POD, CSI and FAR to evaluate the qualities of the errors between the satellite products and the gauge station data.
The results of the statistical analyses are shown in Table 4. CMORPH_RAW FBI is better than CMORPH_BLD, such that the POD and CSI are low and FAR is high.TRMM3BV42 FBI is the closest to the gauge data of all the products but the POD and CSI are lower and the FAR is higher than for GSMaP_RENALYSIS.Therefore, GSMaP_RENALYSIS can more accurately detect rain events than the other products.

Typical Regional Analysis
In this study, we collected data from 2280 rain gauge stations to analyze the precipitation in 2013 over the Tibetan Plateau and South China.The total precipitation was mapped in all four seasons and the CC, RMSE and BIAS were calculated in every month, which are shown in Figures 9-12.Tibetan Plateau is a typical of rainfall distribution including humid regions, semi-humid regions and semi-arid regions, which the distribution is influenced by atmospheric circulation and topography.Rain gauge stations are sparsely distributed over the Tibetan Plateau and the elevation ranges from 86 m to 8755 m and there are exceedingly few gauge stations built in areas with high elevation (Figure 9a). Figure 9b shows the spatial distribution of the total precipitation in four seasons.The total precipitation is relatively low and the spatial distribution trend of precipitation is decreasing from southeast to northwest.The satellite precipitation distribution is similar to the gauge stations but some regions show considerable differences.All products exhibited overestimation in the winter over the western regions, which the result related to the no-gauge stations and interpolated technique.All products overestimated the precipitation in the summer, excepting CMORPH_BLD.CMORPH_RAW overestimated the precipitation and its RMSE and BIAS showed large fluctuating patterns, demonstrating that uncorrected satellite products perform poorly over semi-arid and arid areas.CMORPH_RAW and CMORPH_BLD show some abnormal zero precipitation over Tibetan Plateau.It is possible that the ice and snow surface tend to influence PWM-based accurate retrieval over land [37].The RMSE of GSMaP_RENALYSIS showed large differences in four seasons and GSMaP_RENALYSIS produces underestimation in the spring and overestimation in the summer and autumn.
To evaluate the performances of the satellite-based precipitation measurements, we calculated the CC, RMSE and BIAS values for each month, as shown in Figure 10.The CC of all the products are high from April to October (more than 0.6) and are relatively low in those months when the climate is dry.This result can be explained by the following reasons: (1) compared with weak precipitation, satellite sensors are more able to obtain heavy rainfall events; (2) Over the higher elevation mountainous area of eastern China, it is difficult to obtain accurate precipitation values via remote sensing methods due to the snow and ice on the surface in the winter; (3) Because of the complicated land surface characteristics and their effects on the upward microwave radiation, impacting SSM/I, AMSR-E and AMSU-B, it is difficult for the satellites to resolve rainfall over areas with low rainfall amounts [44].PERSIANN_CDR has the worst CC among the satellite products and its RMSE and BIAS values show the largest errors, which can be explained by the poor input IR data over Tibetan Plateau.Overall, GSMaP_RENALYSIS has the best performance of all the products over the Tibetan Plateau and PERSIANN_CDR has the worst performance, both seasonally and monthly.There are more rain gauge stations and precipitation events in southern China than over the Tibetan Plateau, so we researched the regions including Fujian, Guangdong, Jiangxi, Hunan and Guangxi.The area of study and the elevation data are shown in Figure 11.The seasonal precipitation events detected by the satellites are similar to the in-situ measurements and the max annual precipitation is approximately 1500 mm in the spring and summer over Guangdong as well as over the south of Guangxi.This might be related to the period of precipitation: the study area is mainly affected by the southeast monsoon from the Pacific Ocean and the southwest monsoon from the India Ocean, which the period of precipitation is from May to September.The north of Hunan and Guangdong had more precipitation events than the other regions in all seasons except winter.Jiangxi and Fujian had more precipitation in the winter.This regional precipitation is mainly affected by the cold and warm air in the south.Note that little precipitation was measured over Hunan because of the effects of summer winds leading the precipitation northward in the summer.The CC, RMSE and BIAS values were calculated for each month and were compared with in situ measurements, as shown in Figure 12.The CC of all products are low in January, October and December and the RMSEs are high in the summer.The BIAS is near zero, except in CMORPH_RAW.CMORPH_RAW shows poor CC values, with a low near zero occurring in January.CMORPH_BLD and TRMM3BV42 show highly consistent precipitation distributions and both have higher CC than the other products.However, the CC of CMORPH_BDL is higher than that of TRMM3BV42 in January and November and its BIAS is lower.According to the precipitation distribution and the calculated statistical correlations, the CMORPH_BDL product is consistent with the in-situ measurements.There are more rain gauge stations and precipitation events in southern China than over the Tibetan Plateau, so we researched the regions including Fujian, Guangdong, Jiangxi, Hunan and Guangxi.The area of study and the elevation data are shown in Figure 11.The seasonal precipitation events detected by the satellites are similar to the in-situ measurements and the max annual precipitation is approximately 1500 mm in the spring and summer over Guangdong as well as over the south of Guangxi.This might be related to the period of precipitation: the study area is mainly affected by the southeast monsoon from the Pacific Ocean and the southwest monsoon from the India Ocean, which the period of precipitation is from May to September.The north of Hunan and Guangdong had more precipitation events than the other regions in all seasons except winter.Jiangxi and Fujian had more precipitation in the winter.This regional precipitation is mainly affected by the cold and warm air in the south.Note that little precipitation was measured over Hunan because of the effects of summer winds leading the precipitation northward in the summer.The CC, RMSE and BIAS values were calculated for each month and were compared with in situ measurements, as shown in Figure 12.The CC of all products are low in January, October and December and the RMSEs are high in the summer.The BIAS is near zero, except in CMORPH_RAW.CMORPH_RAW shows poor CC values, with a low near zero occurring in January.CMORPH_BLD and TRMM3BV42 show highly consistent precipitation distributions and both have higher CC than the other products.However, the CC of CMORPH_BDL is higher than that of TRMM3BV42 in January and November and its BIAS is lower.According to the precipitation distribution and the calculated statistical correlations, the CMORPH_BDL product is consistent with the in-situ measurements.There are more rain gauge stations and precipitation events in southern China than over the Tibetan Plateau, so we researched the regions including Fujian, Guangdong, Jiangxi, Hunan and Guangxi.The area of study and the elevation data are shown in Figure 11.The seasonal precipitation events detected by the satellites are similar to the in-situ measurements and the max annual precipitation is approximately 1500 mm in the spring and summer over Guangdong as well as over the south of Guangxi.This might be related to the period of precipitation: the study area is mainly affected by the southeast monsoon from the Pacific Ocean and the southwest monsoon from the India Ocean, which the period of precipitation is from May to September.The north of Hunan and Guangdong had more precipitation events than the other regions in all seasons except winter.Jiangxi and Fujian had more precipitation in the winter.This regional precipitation is mainly affected by the cold and warm air in the south.Note that little precipitation was measured over Hunan because of the effects of summer winds leading the precipitation northward in the summer.The CC, RMSE and BIAS values were calculated for each month and were compared with in situ measurements, as shown in Figure 12.The CC of all products are low in January, October and December and the RMSEs are high in the summer.The BIAS is near zero, except in CMORPH_RAW.CMORPH_RAW shows poor CC values, with a low near zero occurring in January.CMORPH_BLD and TRMM3BV42 show highly consistent precipitation distributions and both have higher CC than the other products.However, the CC of CMORPH_BDL is higher than that of TRMM3BV42 in January and November and its BIAS is lower.According to the precipitation distribution and the calculated statistical correlations, the CMORPH_BDL product is consistent with the in-situ measurements.

Conclusions
This paper mainly discusses the CMORPH_BLD, CMORPH_RAW, TRMM3BV42, PERSIANN_CDR and GSMaP_RENALYSIS precipitation products and their spatiotemporal precipitation distributions based on daily, monthly, seasonal and annual values.We compared the satellite products with data from rain gauge stations using statistical tools, such as their CC, BIAS and RMSE, for the period of 2005 to 2013.The five satellite rainfall products based on four different methods were analyzed via a quantitative study of their error features and the following points were concluded: (1) According to the nine-year daily mean precipitation over China (Figures 2 and 3), the characteristics of precipitation is clearly gradually increasing in the south and decreasing in the northern regions of the Qinling Mountain and Huaihe River.There is more than 10 mm/day in the south and less than 5 mm/day in the northwest.Among the five satellite products, GSMaP_RENALYSIS suggests the best performance with the highest of R (0.91) and the lowest RMSE (0.85 mm/day).CMORPH_RAW demonstrates the poorest capability with the lowest of R (0.73) and higher RMSE (2.79 mm/day).For the CC, RMSE and BIAS each annual daily mean precipitation from 2005 to 2013 (

Conclusions
This paper mainly discusses the CMORPH_BLD, CMORPH_RAW, TRMM3BV42, PERSIANN_CDR and GSMaP_RENALYSIS precipitation products and their spatiotemporal precipitation distributions based on daily, monthly, seasonal and annual values.We compared the satellite products with data from rain gauge stations using statistical tools, such as their CC, BIAS and RMSE, for the period of 2005 to 2013.The five satellite rainfall products based on four different methods were analyzed via a quantitative study of their error features and the following points were concluded: (1) According to the nine-year daily mean precipitation over China (Figures 2 and 3), the characteristics of precipitation is clearly gradually increasing in the south and decreasing in the northern regions of the Qinling Mountain and Huaihe River.There is more than 10 mm/day in the south and less than 5 mm/day in the northwest.Among the five satellite products, GSMaP_RENALYSIS suggests the best performance with the highest of R (0.91) and the lowest RMSE (0.85 mm/day).CMORPH_RAW demonstrates the poorest capability with the lowest of R (0.73) and higher RMSE (2.79 mm/day).For the CC, RMSE and BIAS each annual daily mean

Conclusions
This paper mainly discusses the CMORPH_BLD, CMORPH_RAW, TRMM3BV42, PERSIANN_CDR and GSMaP_RENALYSIS precipitation products and their spatiotemporal precipitation distributions based on daily, monthly, seasonal and annual values.We compared the satellite products with data from rain gauge stations using statistical tools, such as their CC, BIAS and RMSE, for the period of 2005 to 2013.The five satellite rainfall products based on four different methods were analyzed via a quantitative study of their error features and the following points were concluded: (1) According to the nine-year daily mean precipitation over China (Figures 2 and 3), the characteristics of precipitation is clearly gradually increasing in the south and decreasing in the northern regions of the Qinling Mountain and Huaihe River.There is more than 10 mm/day in the south and less than 5 mm/day in the northwest.Among the five satellite products, Shuttle Radar Topographic Mission (SRTM on http://srtm.csi.cgiar.org/)and NOAA (https://rda.ucar.edu/)scientific team for the provision of satellite data and situ measurement data utilized in this study.The authors acknowledge the China Meteorological Service Network and European Center, which made gauge station and satellite data available to us.

Figure 5 .
Figure 5. (a-c) denote the CC, BIAS and RMSE values of the satellite products and rain gauge stations for the seasonal daily mean precipitation values, respectively.

Figure 7 .
Figure 7. Probability density function of daily rainfall for the no-rain case.Figure 7. Probability density function of daily rainfall for the no-rain case.

Figure 7 .
Figure 7. Probability density function of daily rainfall for the no-rain case.Figure 7. Probability density function of daily rainfall for the no-rain case.

Figure 7 .
Figure 7. Probability density function of daily rainfall for the no-rain case.

Figure 8 .
Figure 8. Probability density function of the daily rainfall events with different intensities.Figure 8. Probability density function of the daily rainfall events with different intensities.

Figure 8 .
Figure 8. Probability density function of the daily rainfall events with different intensities.Figure 8. Probability density function of the daily rainfall events with different intensities.

Figure 9 .
Figure 9.The precipitation distributions of the rain gauge stations (a) and multiple satellite (b) in 2013 over southern China.

Figure 10 .
Figure 10.(a-c) denote the CC, RMSE and BIAS among a variety of products for each month in 2013 over the Tibetan Plateau, respectively.

Figure 9 .Figure 9 .
Figure 9.The precipitation distributions of the rain gauge stations (a) and multiple satellite (b) in 2013 over southern China.

Figure 10 .
Figure 10.(a-c) denote the CC, RMSE and BIAS among a variety of products for each month in 2013 over the Tibetan Plateau, respectively.

Figure 10 .
Figure 10.(a-c) denote the CC, RMSE and BIAS among a variety of products for each month in 2013 over the Tibetan Plateau, respectively.

Figure 11 .
Figure 11.The precipitation distributions from the rain gauge stations (a) and multiple satellites (b) in 2013 over southern China.

Figure 12 .
Figure 12. (a-c) denote the CC, RMSE and BIAS values of the studied products for each month in 2013 over southern China, respectively.

Figure 11 .Figure 11 .
Figure 11.The precipitation distributions from the rain gauge stations (a) and multiple satellites (b) in 2013 over southern China.

Figure 12 .
Figure 12. (a-c) denote the CC, RMSE and BIAS values of the studied products for each month in 2013 over southern China, respectively.

Figure 12 .
Figure 12. (a-c) denote the CC, RMSE and BIAS values of the studied products for each month in 2013 over southern China, respectively.

Table 1 .
Description of the precipitation datasets utilized in this study.

Table 1 .
Description of the precipitation datasets utilized in this study.

Table 2 .
Contingency table comparing rainfall retrievals by satellites and gauges.

Table 3 .
The CC, BIAS and RMSE values for a variety products and the gauge station data for a year of the daily mean precipitation values over China.

Table 4 .
The results of the statistics of FBI, POD, CSI and FAR.

Table 3
), GSMaP_RENALYSIS gives the highest CC (>0.8) and smallest RMSE (~1.2 mm/day) and the low BIAS (±0.1) in every year.GSMaP_RENALYSIS underestimates precipitation from 2005 to 2006 but it overestimates precipitation from 2007 to