Similarity and Error Intercomparison of the GPM and Its Predecessor-TRMM Multisatellite Precipitation Analysis Using the Best Available Hourly Gauge Network over the Tibetan Plateau

The performance of Day-1 Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM) mission (IMERG) and its predecessor, the Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis 3B42 Version 7 (3B42V7), was cross-evaluated using data from the best-available hourly gauge network over the Tibetan Plateau (TP). Analyses of three-hourly rainfall estimates in the warm season of 2014 reveal that IMERG shows appreciably better correlations and lower errors than 3B42V7, though with very similar spatial patterns for all assessment indicators. IMERG also appears to detect light rainfall better than 3B42V7. However, IMERG shows slightly lower POD than 3B42V7 for elevations above 4200 m. Both IMERG and 3B42V7 successfully capture the northward dynamic life cycle of the Indian monsoon reasonably well over the TP. In particular, the relatively light rain from early and end Indian monsoon moisture surge events often fails to be captured by the sparsely-distributed gauges. In spite of limited snowfall field observations, IMERG shows the potential of detecting solid precipitation, which cannot be retrieved from the 3B42V7 products.


Introduction
Reliable precipitation data and information are essential for understanding the water and energy cycles at both regional and global scales.However, accurate observation of precipitation is challenging in many remote regions of the world, such as the Tibetan Plateau (TP), due to sparse gauge networks and high spatial variability in precipitation [1,2].Satellite remote sensing has provided unprecedented precipitation information at a broader range of time and space scales, representing a significant contribution toward mapping global rainfall [3][4][5].
The Global Precipitation Measurement (GPM) mission is an important new program designed for global satellite precipitation estimation based on an international satellite constellation, which provides precipitation measurements from space at a spatial resolution of 0.1° × 0.1 ° and a half-hourly temporal resolution [6].As the successor of the highly successful TRMM, the GPM mission provides four levels of products based on various algorithms.Level 1 data consist of geolocated, calibrated Dual-frequency Precipitation Radar (DPR) power, GPM Microwave Imager (GMI) brightness temperatures and intercalibrated brightness temperatures from partner radiometers at the Instantaneous Field Of View (IFOV); Level 2 products include geolocated, geophysical data and DPR reflectivity at the IFOV; Level 3 products consist of gridded time-space sampled geophysical data (e.g., precipitation rates) from the GPM Core sensors and partner radiometers; and Level 4 products include merged remotely-sensed and model information [6].The Integrated Multi-satellitE Retrievals for GPM mission (IMERG) was designed to create Level 3 products, a consistently calibrated global precipitation product with appropriate error and metadata information.
The intercomparison of GPM with current satellite precipitation products and corresponding analysis across the TP is urgently needed.Efforts have been made to assess the suitability of satellite precipitation products over the TP [1,[7][8][9][10][11]. Previous studies found that the TMPA precipitation products outperform other satellite products (e.g., PERSIANN and CMORPH) across the plateau with lower errors and biases [9,11].The work in [9] showed that Global Precipitation Climatology Center (GPCC) gauges employed in the post-real-time TMPA system significantly improve the data accuracy.The work in [11] also confirmed that due to the monthly gauge adjustments, the research products of TMPA exhibit great hydrological potential at both monthly and daily scales over the TP.The work in [12] showed that the latest 3B42 Version 7 (3B42V7) clearly improves upon 3B42V6 over China in terms of daily mean precipitation.Here, the TMPA 3B42V7 product of the best performance is selected to compare to the GPM IMERG product.
Despite the aforementioned studies, the newly-launched GPM satellite-based precipitation estimates warrant evaluation over the TP urgently because: (1) the similarity and differences of GPM Day-1 IMERG and TMPA 3B42V7 rainfall products should be investigated; (2) error analysis at the sub-daily scale is seldom performed; and (3) impacts of elevation and latitude on satellite rainfall estimates should be thoroughly explored.
In the cold season (i.e., October to the next March), snowfall occurs frequently across the TP [13,14].Figures S1 and S2 also show snowfall regimes (including snowfall occurrence and amounts) retrieved from the IMERG product for the period October 2014-March 2015.Moreover, the automatic gauge network is generally switched off from October to the next March on the TP [15].Thus, this study focuses on the error analysis and comprehensive assessments of GPM-era Day-1 IMERG and TRMM-era 3B42V7 products only in the warm season of 2014 over the TP.In order to understand solid precipitation regimes from IMERG on the TP, we perform an initial comparison with the field snowfall experiment at a site in the northeastern TP.Results of this study are expected to promote the improvement of future versions of the IMERG algorithms.
This paper is organized as follows: Section 2 provides the details of the study area, dataset and method.Sections 3 and 4 present error analysis and the discussion, as well as constructive suggestions for the IMERG algorithm over regions of complex terrain and high elevations.The main conclusions and future work are provided in Section 5.

Study Region
The TP, with a mean elevation of more than 4000 m above sea level (a.s.l.), is the highest and largest plateau in the world.Topography plays a vital role in the formation of distinct local climates and environments ranging from forests to alpine grassland across the whole plateau [16].The TP is also subjected to the effect of multiple climate systems.Atmospheric circulation patterns over the TP are primarily characterized by the Indian monsoon in the warm season and the mid-latitude Westerlies in the cold season [17].Moreover, the TP is Asia's water tower, being the source region of many major rivers in the Asian continent, e.g., the Brahmaputra, Yellow, Yangtze, Mekong, Salween, and Indus Rivers.The three watersheds (i.e., Brahmaputra, Yellow and Yangtze Rivers) totally within China's territory are selected for basin-scale error analysis (Figure 1).

Satellite Retrievals and Gauge Data
High-resolution IMERG precipitation products have been released since mid-March 2014.It is necessary and urgent to conduct a comprehensive assessment all over the world, especially in areas with complex terrain and climates, such as the TP.A detailed description of IMERG is provided in [6,18].In this research, we use the very early Day-1 stage of the IMERG "Final run" product.The IMERG data are downloaded from the GPM website [19].The IMERG algorithm is designed to intercalibrate, merge and interpolate all possible satellite microwave precipitation estimates, accompanied with microwave-calibrated infrared satellite estimates, gauge precipitation analyses and other potential precipitation estimators at reliable temporal and spatial scales [18].In order to retrieve precipitation from the GPM constellation, the 2014 version of the Goddard Profiling Algorithm (GPROF2014) was used and then gridded, intercalibrated to the GPM Combined Instrument product with a 0.10 degree scale [18].The CMORPH Kalman Filter Lagrangian time interpolation scheme [20] and PERSIANN-CCS re-calibration scheme [21] were also provided in the computation routines.
The rain gauge network used for validation comprises 258 automatic gauges over the whole TP (Table S1).The used hourly gauge data are available from China's hourly gauge network, which consists of a total of ~2000 automatic weather stations distributed across China [15].Hourly precipitation is observed by siphon or tipping-bucket rain gauges and recorded automatically.Hourly precipitation reports from these stations are transferred to China Meteorological Administration (CMA).Because all the stations used are automatic observation systems, and we do not have the bias-adjusted precipitation procedure for automatic rain gauges, the precipitation undercatch is not considered in this study.However, all of the gauge data have undergone strict quality control in three levels by CMA, including: (1) the extreme values' check; (2) internal consistency check; and (3) spatial consistency check [26].The quality-controlled hourly gauge data are used widely for satellite precipitation estimation and high spatiotemporal gauge-satellite precipitation integration in mainland China [10,15,27].To our best understanding, this is the best hourly rainfall observational dataset over the study region.We have some other gridded rainfall datasets, for instance, China Gauge-based Daily Precipitation Analysis (CGDPA).However, this is based on a daily scale that does not meet our requirement for evaluation at the sub-daily scale.In terms of the independence of the intercomparison between gauge-and satellite-based data, there are 19 GPCC stations used for bias adjustment or the update of TRMM/GPM research level products over the TP, accounting for a very small portion of the total gauges (~7%) (Table S1).Therefore, the evaluation and comparison of GPM and TRMM are based on ~93% independent gauge stations.In addition, both TRMM and GPM products combine GPCC gauge data at the monthly scale [18,22], while the analysis of two satellite products are performed at the sub-daily scale.
Considering data availability and quality, this study uses the best hourly gauge data covering the warm period from 1 April-30 September 2014.Among all of the stations, there are 52 in the Yellow River basin, 66 in the Yangtze River basin and 69 in the Brahmaputra basin.The gauge densities in the Yellow River, Yangtze River and Brahmaputra basins are 2.01 and 1.34, 1.78 per ten thousand square kilometers, respectively.Moreover, a field snowfall experiment from April 2014-February 2015 at a site (99°52′54″ E, 38°16′06″ N, 2980 m) of the northeastern TP (Figure 1) was performed by [28], and the data are used for IMERG snowfall error analysis.Details of the field site information can be found from [28].

Evaluation Indicators
To better evaluate the satellite rainfall products over the TP, several statistical indices were selected, including Relative Bias (RB), Standard Deviation (SD), Root Mean Square Error (RMSE), and Correlation Coefficients (CC).Moreover, the Probability Of Detection (POD) and False Alarm Ratio (FAR) were calculated to check the appearance possibility of rainfall events from satellite products.The definitions of the indicators are as follows: (2) where is the number of samples; and stand for individual and averaged satellite rainfall estimates; and indicate individual and averaged gauge-based measurements, respectively.represents the observed rain by gauge, which is detected by satellite correctly, while is the observed rain not detected, and is not the observed rain, but detected falsely.The detailed information of the evaluation indices can be found in [29,30].

Similarity of Spatial Rainfall Patterns
Figure 2 describes spatial patterns of average three-hourly IMERG and 3B42V7 rainfall estimates in the warm season (from April-September) of 2014 over the TP.Both products show a decreasing trend from the southeast to the northwest over the plateau, which are generally consistent with other satellite products (e.g., 3B42RT, CMORPH and PERSIANN) [9][10][11].The average accumulated rainfall in the warm season of 2014 is 347 mm over the whole plateau, and the maximum rainfall primarily occurs in the Himalayan region with a value of 2609 mm.In contrast, most parts of the west and north show lower values, where the Westerlies do not prevail and the Indian monsoon is relatively weaker during the warm period [17].

Pixel-Scale Error Intercomparison
To ensure proper error analysis at the grid cell scale, we firstly check the specified location or pixel for each rain gauge and then evaluate the spatial-temporally coincident satellite pixel with corresponding gauges over the TP.Note that the zero-event records are not considered as the rainfall events during the calculation.
In general, the IMERG product showed a limited capability in estimating hourly rainfall in the warm season of 2014, with the averaged CC of 0.32 for all gauges across the TP (Table 1).Higher CC was observed in the northeast and southeast TP (Figure 3).Along the Himalayas (i.e., the southern part of the TP), CC was smaller with the mean value less than 0.40.The weaker correlation was attributed to the complex topography in the Himalayas, which might influence the retrieval of satellite precipitation [10].In spite of higher elevation (i.e., more than 4000 m a.s.l.) in the central TP, the mean CC was more than 0.60.At the three-hourly scale, IMERG and 3B42V7 showed similar RB, RMSE and CC patterns in space across the TP (Figure 3).Both of them showed overestimation in the south and west and underestimation in the northeast.However, several negative RB values of IMERG rainfall estimates were detected along the southern boundary of the TP, whereas 3B42V7 generally showed positive RB values.In addition, lower absolute RB values were mainly distributed in the eastern TP for both products at the three-hourly scale.Over the plateau, both IMERG and 3B42V7 showed overall overestimation of 28.1% and 48.6% at the three-hourly scale, respectively.As for RMSE, both products decreased gradually from the southeast to northwest across the TP, with averaged values of 3.26 mm (IMERG) and 3.47 mm (3B42V7), respectively (Table 1).3B42V7 showed slightly worse agreement with gauge data than IMERG, with the CC value of 0.42 at the three-hourly scale.Spatially, the CC values for both products were generally higher in the east than those in the west.

Basin-Scale Error Intercomparison
To facilitate direct comparison between IMERG and 3B42V7 products, the spatial resolutions of gauge-and satellite-based products should be consistent firstly.The gauge data available in each basin (including the Yellow, Yangtze and Brahmaputra basins) were initially interpolated to a 0.10° grid with the Inverse Distance Weighting (IDW) method, which had been shown to be reliable and efficient in many previous basin-scale studies [11,[31][32][33][34].In order to make the different data comparable, the 3B42V7 data were resampled to the 0.10° × 0.10° resolution, the same as IMERG and the interpolated gauge data, by using the standard bilinear interpolation method [35].We adopt a consistent common 0.10° grid cell for both satellite rainfall comparison and validation.Note that potential errors would possibly be introduced during the resampling procedure.
At the three-hourly scale, both IMERG and 3B42V7 showed similar bias patterns and bias ratio distribution in the Yellow and Yangtze River basins (Figure 4).In the Yellow River basin, the two products showed apparent underestimation in the north and south parts and overestimation in the other parts (e.g., central and southeast); whereas in the Yangtze River basin, both showed overestimation in the west and east and underestimation in the central part.Furthermore, IMERG showed lower biases than 3B42V7 for both the Yellow and Yangtze River basins.The averaged RB values in the Yellow River basin were 3.7% and 5.8% for IMERG and 3B42V7, respectively.In the Yangtze River basin, the RB values were 10.8% (IMERG) and 13.8% (3B42V7), respectively.In the Brahmaputra basin, IMERG showed obvious overestimation with an overall mean RB value of 35.6% and SD value of 49.5%, respectively.3B42V7 showed more significant overestimation around the whole basin, and the averaged RB and SD values were 92.6% and 59.6%, respectively (Table 2).In short, both satellite products showed more acceptable performance in the Yellow and Yangtze River basins than in the Brahmaputra basin.In addition, time series analyses of IMERG, 3B42V7 and gauge-based rainfall are shown for each major river basin during the warm season of 2014 (Figure 5).Although both satellite products showed large standard deviations, IMERG did not have extreme values and appeared to be a more acceptable satellite rainfall product in each basin.

Elevation Impact Analyses
All gauges in the TP were rearranged in Figures 6-8 according to their elevations ranging from 1200 m-5135 m.A comparison of CC, POD and FAR for IMERG and 3B42V7 along elevated sites is presented at the three-hourly scale in Figure 6.All of the metrics showed similar trends for both satellite products as elevations increased from 1200 m-5135 m.No significant CC differences occurred in the two satellite products along the elevated stations.Furthermore, there was no apparent difference for POD and FAR values between IMERG and 3B42V7 products with increasing elevations.3B42V7, however, showed slightly higher POD for the elevation above 4200 m and obviously higher FAR for the elevation below 3000 m.This implied that regions characterized by higher altitudes (e.g., >4200 m) might still be problematic for the Day-1 multisatellite retrieval in GPM over the TP.Note that potential errors originating from different places under various rainfall regimes might influence the statistical results.(a-c) the Correlation Coefficients (CC), Probability Of Detection (POD) and False Alarm Ratio (FAR) between satellite data and gauge observations as the elevation increases.Note that the grey square stands for 258 gauges recorded according to elevation; the red lines stand for IMERG estimates; the blue lines stand for 3B42V7 products and the comparison was calculated at 3-hourly scale.
Figure 7 shows daily rainfall regimes, i.e., maximum, mean and conditional mean values, with elevations ranging from 1200 m-5135 m in the warm season of 2014 over the TP.The conditional mean rainfall is defined as the average value among the detected rainfall days within the warm season of 2014 at each station.There was a gradually downward trend for the daily rainfall, including maximum, mean and conditional mean values, below the elevation of 3000 m.When the elevation exceeds 3000 m, no clear trend was detected for each daily rainfall regime.Regarding the daily mean rainfall, high consistency was observed among IMERG, 3B42V7 products and gauge measurements for various elevations.However, there were large biases between IMERG and 3B42V7 in terms of daily maximum and conditional mean rain rates at elevations higher than 3000 m.This phenomenon was more significant in the conditional mean value due to the removal of multiple zero rainfall events.Compared to 3B42V7, IMERG showed higher biases regarding maximum daily rainfall at elevations of approximately 4200 m and slightly higher positive errors in the conditional mean daily estimates at most sites.Rain probability is a description of the likelihood of rainfall within a specified period and location.We define the warm season of 2014 as the specified period and calculate the percentages of rainfall days during this period at each station.Along the elevated stations, there was no noticeable trend for the raining probability (Figure 8a).Significantly positive errors were observed from 3B42V7 estimates for raining probability, indicating that 3B42V7 overestimated rainfall events along the elevated gauges.The classification of raining intensities (i.e., 0-10 mm/d for light rainfall, 10-25 mm/d for moderate rainfall and >25 mm/d for heavy rainfall) is primarily from [9,36].Less significant biases were found in IMERG for light rainfall (Figure 8b).Compared to light and heavy rainfall, moderate rainfall derived from the IMERG and 3B42V7 products matched the gauge measurements well along the elevated sites (Figure 8c).In addition, light rainfall events were far more frequent than moderate and heavy rainfall events for both satellite products, where the proportion of heavy rainfall was relatively low, with the mean value lower than 5%.This means that the sample sizes to verify heavy rainfall events from a single warm season are possibly small, which may introduce a potential uncertainty for satellite rainfall estimates.

Latitude Impact Analyses
The 258 gauges over the TP were reorganized with increasing latitudes (northward), and their daily rainfall intensity was plotted according to time and latitude.Figure 9 presents a series of continuous rainfall events from the gauge, IMERG and 3B42V7 datasets in the warm season of 2014 when the Indian monsoon prevailed over the plateau.Before June, sporadic rainfall activities occurred from satellite products (i.e., IMERG and 3B42V7) and gauge measurements.Additionally, there was intermittent rainfall activity captured by both satellite products at the beginning of the Indian monsoon in early May 2014.It gradually weakened from latitudes of 26.5°-32.5° in the plateau.
Moisture surge events dominated the plateau from early June to mid-September in 2014, which were characterized by abundant rainfall activities from north to south for both satellite products and gauge observations.For instance, there were several latitudinal gradients of rainfall occurring from mid-June to the end of August, which were regarded as typical moisture surge events over the TP.Within the gauge-based domain, an insignificant latitudinal gradient existed during the mature phase of the Indian monsoon when rainfall should be relatively widespread across all latitudes.Over the plateau, rain gauges are often installed at low altitudes or sites with residents.It is worth noting that [37] proposed that moisture surge events tend to influence low-elevation stations more than middle or high elevations, and the gauge-based rainfall presented in Figure 9 might ignore such an elevation-specific response.
In addition, ground-based gauges failed to capture rainfall events during the end of the Indian monsoon period in September.We deduce that sparse gauges are the primary reason; however, the specific reason needs to be further investigated in the future.Overall, both satellite products could detect rainfall activities synchronously under the domination of Indian monsoon periods, and IMERG does not show obvious superiority to 3B42V7 during the comparison.

Snowfall Regime Analyses
Based on limited field snowfall data at a site of the northeastern TP (Figure 1), we initially analyzed IMERG-retrieved snowfall characteristics.During the sixteenth snowfall days, IMERG captured four times which primarily have larger snowfall amounts (Figure 10).Although the probability of detection on IMERG snowfall was not high (i.e., 25%) at the experimental site, the error values were not much reduced for the four right times in the high-elevation region (Table 3).

Discussion
Several studies [8][9][10][11] have been performed to evaluate satellite precipitation products over the TP.However, none of these studies were concerned with sub-daily scale examination in this region.It should also be emphasized that the performance of the GPM-era Day-1 IMERG product was examined over regions of high elevation and complex terrain for the first time.The differences between two satellite products (i.e., IMERG and 3B42V7) and gauge data showed clearly whether the IMERG datasets perform better.Note that the resampling and interpolation methods are used to make the different satellite (i.e., IMERG and 3B42V7) and gauge data comparable.More studies will be carried out to evaluate the possible errors introduced by the above two methods.
Satellite precipitation products can provide alternatives for estimating rainfall data and create new opportunities for hydrologists in understanding and applying remotely-sensed information.The work in [38] forced the Coupled Routing and Excess STorage (CREST) hydrological model with 3B42V7 products and yielded comparable skill scores with in situ gauges even without recalibration of the hydrological model by the satellite precipitation.TRMM-era research products also showed a satisfactory capability as the forcing of the Variable Infiltration Capacity (VIC) model to simulate streamflow in basins [11,39,40].GPM is regarded as the successor of TRMM for satellite precipitation estimates.Although no significant advantages were observed for the GPR-era IMERG product at the sub-daily scale over the TP, it seems to be acceptable for regional hydrological studies.Given the sparse distribution of rain gauges over the TP, especially in the west and central part, the IMERG datasets are extremely valuable in deriving temporal and spatial patterns of precipitation with high resolution.It can also be a substitute for the TRMM-era rainfall product, such as 3B42V7.
The results of IMERG and 3B42V7 data assessment in river basins, especially over the upper reaches of the Yellow and Yangtze Rivers, showed a new opportunity for basin-scale hydrological applications.However, a few extreme high positive or negative relative biases (i.e., >100% and <−50%) were detected in the Brahmaputra basin.The work in [41] indicated that regions characterized by complex terrain and a rigid climate would still be problematic for the Day-1 IMERG and 3B42V7 products under current observing skills.We found that the complex terrain and local strong convectional weather significantly influence the satellite rainfall retrieval and result in unexpected errors between satellite estimates and gauge observations in the Brahmaputra basin.For satellite rainfall products over land, the passive microwave algorithms depend primarily on scattering by ice, but the orographic rains might not produce much ice aloft [42].Thus, the noticeable IMERG and 3B42V7 rainfall errors in the Brahmaputra basin might be due to complex orography and associated warm-rain process [43,44].Furthermore, it is noted that discrepancies between satellite products and gauge measurements might be due to the shortcomings in the automatic gauge network.Although the gauge data have gone through three levels of quality control, systematic biases induced by wind and blown snow undercatch, wetting and evaporation errors may contribute to potential sources of uncertainty in mountainous regions over the TP [2].Moreover, nearly no rain gauges are installed over the western and central TP.More efforts are therefore urgently needed to build a denser rain gauge observation network in these regions.
The influence of the Indian monsoon is critical for the advection of heat and moisture and climate patterns in the TP region [17,45].The continuous domains of IMERG and 3B42V7 precipitation estimates help characterize rainfall space-time structures under the domination of the Indian monsoon in the warm season over the TP.The onset, maturation and dissipation of the Indian monsoon were successfully reflected from the satellite rainfall estimates, whereas some moisture surge activities were not successfully detected by the sparse gauge-based observations.Our efforts of examining the variability in satellite rainfall estimates with latitude are helpful for understanding the moisture flux transport across the TP.Because of limited rainfall records, the Westerlies were not included for understanding the whole climate system in the plateau.

Conclusions
This study provides an early and timely quantitative study of the error characteristics for GPM and its predecessor-TRMM rainfall products with high temporal resolution (i.e., three-hourly) across the TP.The intercomparison is based on multiple evaluations with the best available rain gauge network at the grid cell and basin scales.Furthermore, we have investigated the impact of elevation and latitude on rainfall analysis.In addition to the IMERG and 3B42V7 satellite rainfall error analysis, the snowfall error analysis from IMERG is carried out at a site of the TP.Primary conclusions are summarized as follows: (1) Overall, IMERG and 3B42V7 show similar rainfall patterns across the TP, with a decreasing trend from the southeast to the northwest.However, IMERG performs better, with slightly higher correlation and lower bias than 3B42V7 at the three-hourly scale, although both datasets show similar spatial patterns for all assessment indicators.(2) As elevations increase, there is no obvious difference in CC between IMERG and 3B42V7.
However, IMERG shows slightly lower POD for elevations above 4200 m and significantly smaller FAR for elevations below 3000 m.In addition, higher consistency is detected from the IMERG products for light rainfall than 3B42V7 estimates.(3) The latitudinal rainfall gradients of the Indian monsoon dynamics are successfully detected by the two satellite estimates.In particular, the relatively light rain from the early and end Indian monsoon moisture surge events is often not observed by the sparse gauges over the TP.(4) Over high and complex TP regions, this study provides the very first understanding of GPM at the hourly scale and confirms the desired results that the GPM-era Day-1 Level-3 products perform at least equally well, often times better for light and solid precipitation at high elevations, than the TRMM V7 mature products.
The GPM Core Observatory sensors are expected to have the ability to detect cold-season solid precipitation, but only one snowfall site in the northeastern TP is available for evaluating the IMERG snowfall product.Conducting snowfall observation and evaluation across the TP are our ongoing work.

Figure 1 .
Figure 1.Overview of the geographical and topographical location of the Tibetan Plateau (TP) and three main river basins, where the red, blue and black polygons indicate the basin of the Yellow River, Yangtze River and Brahmaputra, respectively; the black dots show rain gauges, and the red rectangle shows the field experimental site for snowfall observation.

Figure 2 .
Figure 2. Maps of the mean three-hourly (a) IMERG and (b) 3B42V7 rainfall estimates in the warm season of 2014 over the TP.

Figure 3 .
Figure 3. Maps of (a,b) Relative Bias (RB); (c,d) Root Mean Square Error (RMSE); (e,f) Correlation Coefficients (CC) and (g,h) scatter diagram between satellite and gauge-based three-hourly rainfall estimates in the warm season (April-September) of 2014 over the TP.Notes: the left panel, i.e., (a,c,e,g), stands for IMERG rainfall estimates, and the right panel, i.e., (b,d,f,h) stands for 3B42V7 products.

Figure 4 .
Figure 4. Spatial patterns of averaged three-hourly relative bias (RB) (left and middle column) and its percentage of RB basin area (%) (right column) derived from satellite products and gauge observations in the three main river basins, where the top panel, i.e., (a-c), is the Yellow River basin; the middle panel, i.e., (d-f), is the Yangtze River basin; and the bottom panel, i.e., (g-i), is the Brahmaputra basin.The left column stands for the IMERG rainfall estimates and the middle side is the 3B42V7 data.

Figure 5 .
Figure 5. Variation of area mean daily rainfall for each river basin, i.e., (a) Yellow River; (b) Yangtze River and (c) Brahmaputra, in the warm season of 2014.The black line denotes gauge observations; the red one stands for the IMERG estimates; and the blue one is the 3B42V7 products.

Figure 6 .
Figure 6.Elevation-related analyses of three-hourly rainfall indices in the warm season of 2014:(a-c) the Correlation Coefficients (CC), Probability Of Detection (POD) and False Alarm Ratio (FAR) between satellite data and gauge observations as the elevation increases.Note that the grey square stands for 258 gauges recorded according to elevation; the red lines stand for IMERG estimates; the blue lines stand for 3B42V7 products and the comparison was calculated at 3-hourly scale.

Figure 7 .
Figure 7.Comparison of IMERG and gauge-based daily rainfall in the warm season of 2014 over the TP: (a-c) the maximum, mean and conditional mean of daily rainfall in the warm season as the elevation increases.Note that the grey square stands for 258 gauges recorded according to elevation; the black line denotes gauge observations; the red one stands for the IMERG estimates; and the blue one is the 3B42V7 products.

Figure 8 .
Figure 8.Comparison of IMERG and gauge-based daily rainfall in the warm season of 2014 over the TP: (a-d) the probability of raining days, light rainfall (0-10 mm/d), moderate rainfall (10-25 mm/d) and heavy rainfall (>25 mm/d) in the warm season as the elevation increases.Note that the grey square stands for 258 gauges recorded according to elevation; the black line denotes gauge observations; the red one stands for the IMERG estimates; and the blue one is the 3B42V7 products.

Figure 9 .
Figure 9.Comparison of satellite products and gauge observations with latitudes in the warm season of 2014 over the TP: (a) the 258 gauges recorded according to their latitudes; (b-d) the surging of Indian monsoon illustrated from the gauge, IMERG and 3B42V7 datasets.Note that the unit in (b-d) is mm per day (mm/d).

Figure 10 .
Figure 10.IMERG retrieved total precipitation and snowfall amount comparison within the ground observed snowfall at the field experimental site of northeastern TP.

Table 1 .
Mean values of Correlation Coefficients (CC), Root Mean Square Error (RMSE) (mm) and Relative Bias (RB) (%) for IMERG and 3B42 Version 7 (3B42V7) rainfall estimates at hourly and 3-hourly timescales in the warm season of 2014 over the Tibetan Plateau (TP).

Table 2 .
List of Relative Bias (RB) features (i.e., minimum, mean and maximum) and corresponding Standard Deviation (SD) for IMERG and 3B42V7 rainfall estimates at hourly and 3-hourly timescales in the Yellow River, Yangtze River and Brahmaputra basins in the warm season of 2014 over the Tibetan Plateau (TP).

Table 3 .
The successfully detected snowfall amounts' comparison between IMERG retrieval and field snowfall observation at a site of northeastern TP.