Application of the GPM-IMERG Products in Flash Flood Warning: A Case Study in Yunnan, China

: NASA’s Integrated Multi-Satellite Retrievals for Global Precipitation Measurement (IMERG) is a major source of precipitation data, having a larger coverage, higher precision, and a higher spatiotemporal resolution than previous products, such as the Tropical Rainfall Measuring Mission (TRMM). However, there rarely has been an application of IMERG products in ﬂash ﬂood warnings. Taking Yunnan Province as the typical study area, this study ﬁrst evaluated the accuracy of the near-real-time IMERG Early run product (IMERG-E) and the post-real-time IMERG Final run product (IMERG-F) with a 6-hourly temporal resolution. Then the performance of the two products was analyzed with the improved Rainfall Triggering Index (RTI) in the ﬂash ﬂood warning. Results show that (1) IMERG-F presents acceptable accuracy over the study area, with a relatively high hourly correlation coe ﬃ cient of 0.46 and relative bias of 23.33% on the grid, which performs better than IMERG-E; and (2) when the RTI model is calibrated with the gauge data, the IMERG-F results matched well with the gauge data, indicating that it is viable to use MERG-F in ﬂash ﬂood warnings. However, as the ﬂash ﬂood occurrence increases, both gauge and IMERG-F data capture fewer ﬂash ﬂood events, and IMERG-F overestimates actual precipitation. Nevertheless, IMERG-F can capture more ﬂood events than IMERG-E and can contribute to improving the accuracy of the ﬂash ﬂood warnings in Yunnan Province and other ﬂood-prone areas.


Introduction
Flash floods, triggered by heavy rainfall (i.e., short duration, high intensity), are the rapid flooding of water within minutes up to several hours in small basins (hundred square kilometers or less) [1]. It is one of the most devastating floods in the world, which can cause great economic losses and casualties. and rainfall intensity in the prediction of flash floods. It focuses on antecedent conditions and has been put into practice for over 10 years in Taiwan [13]. However, this method relies too much on rainfall stations, fails to fully consider intermittent rainfall, and is greatly affected by rainfall field segmentation. Meanwhile, using the deduction coefficient of "t" days in the RTI model for the antecedent rainfall calculation will result in a higher false alert rate under some rainfall patterns. Therefore, Chen et al. (2017) proposed an improved RTI to calculate the antecedent rainfall and effective accumulated rainfall to solve the abovementioned problems, which achieved good practical application effects in Taiwan's 2017 disaster warning [16].
As satellite precipitation data is widely used in large-scale watershed hydrological simulation or land surface process simulation, it is still difficult to meet the needs of flash flood prediction in small and medium-sized watersheds. Most of the existing research focus on the corrected products, but there scarcely has been application verification for real-time products. Most of the satellite products are applied on a large scale, and little attention has been paid to their applicability in small-scale flash flood warnings. Additionally, flash floods generally occur in small and medium-sized river basins with poor economic conditions, low station network coverage density, and there is great difficulty in obtaining data. Given that IMERG products have high spatiotemporal precision, it is therefore of great practical significance to evaluate the applicability of these products to flash flood warnings. Moreover, regional studies are also very common and popular because precipitation exhibits strong spatial variations and different products show varying performance over different regions. Meanwhile, there are still few studies on the inter-comparison of the IMERG products, especially in China, where the superiority of its application in flash flood warnings still needs further exploration [17]. Therefore, taking Yunnan Province in China as the study area, based on two fifth-generation IMERG products (IMERG-E, IMERG-F) and China Meteorological Administration (CMA) data, this study first evaluated these two products with particular attention paid to their systematic and random errors. Then, the empirical RTI method was utilized to evaluate the applicability of the different IMERG products in flash flood warnings. This paper is further organized as follows: Section 2 describes the materials and data; Section 3 presents an assessment of IMERG-E and IMERG-F based on locally measured data and analyzes the effects of flash flood warnings based on satellite precipitation data; and the conclusions are summarized in Section 4.

Study Area
Yunnan Province is located at the low latitude plateau of southwest China, with an area of 390,000 km 2 , accounting for 4.1% of China's total area [18]. The terrain in Yunnan Province is very complex, coupled with the low economic level of remote mountain areas, resulting in the lack of early disaster data. For example, the main terrains are mountains (84%), plateaus and hills (10%), and basins and valleys (6%). The altitude is high in the northwest and low in the southeast, descending stepwise from north to south in this region, and the average altitude is 1980 m, and 87.21% of the land is located in middle-latitude areas (1000-3500 m). The area with a slope below 25 • accounts for 56.46% of the national total area. Besides, Yunnan Province is characterized by the most widely distributed karst topography, accounting for 28.9% of the province's area; this results in low water retention and frequent flood disasters. Yunnan is rich in water resources, but unevenly distributed in time and space, where more than 70% of the water resources are concentrated in remote mountainous areas. The total population is 47.1 million in 2014, accounting for 3.5% of the national total population. Yunnan has a typical monsoon climate, mainly dominated by the east Asian monsoon and southwest monsoon, and 90% of the precipitation is concentrated from May to October, especially in the flood season (June to August). For example, long-term studies have shown that the critical rainfall is 35-200 mm in northwest Yunnan, 50-200 mm in southwest Yunnan, and 100-300 mm in eastern Yunnan [19]. The above reasons led to frequent flash flood disasters, which has caused extremely serious disasters.
For example, 72 people died from flash floods in Yunnan Province in 2014, accounting for 22.2% of the national total deaths from flash floods. All the preliminary works have promoted the development of flash flood prevention in Yunnan Province [20], but further exploration is still needed, and therefore Yunnan Province was selected as the study area in this study. Figure 1 shows the map of the topography and flash flood disaster distribution in Yunnan Province.
Remote Sens. 2020, 12, x FOR PEER REVIEW  4 of 16 is still needed, and therefore Yunnan Province was selected as the study area in this study. Figure 1 shows the map of the topography and flash flood disaster distribution in Yunnan Province.

Satellite Data
TRMM Multi-Satellite Precipitation Analysis (TMPA) (50°N-50°S) can retrieve microwaveinfrared satellite precipitation estimates with gauge adjustments and can generate rational precipitation estimates at fine spatiotemporal scales (0.25° × 0.25° and 3-hourly temporal intervals) within the global scope (60° N-60° S) [21,22]. As the successor of TRMM, the GPM constellation consists of a core observation platform and 10 cooperative satellites to observe global precipitation through inter-satellite cooperation. The GPM program offers three different levels of data products, all of which are available on the NASA website (http://www.nasa.gov/mission_pages/GPM). This study focuses on the GPM Level 3 products based on the IMERG algorithm. IMERG not only has a high spatial resolution (0.1° × 0.1°, 0.5 h) and fully global coverage (60° S~60° N), but also predicts flood disasters by reducing the uncertainty associated with short-term precipitation accumulation. At present, IMERG data has been developed from the first version to the sixth version, and IMERG-V06B is the latest version of the satellite rainfall data. However, Mohammad et al. (2019) verified that IMERG V05 performs better than V06, and therefore this study has employed IMERG V05 to estimate the precipitation data [23].
IMERG has three different products: Early, Late, and Final. In the real-time phase, the IMERG data generation system generates Early products after running once and then generates Late products after running again. The main difference between them is that the Early product is generated only by the forward propagation algorithm in the cloud mobile vector propagation algorithm, while the Late product is added with the backpropagation algorithm. The time-lag product Final introduces more sensor data sources based on the Late product. Meanwhile, the time of the extracted satellite data is

Satellite Data
TRMM Multi-Satellite Precipitation Analysis (TMPA) (50 • N-50 • S) can retrieve microwave-infrared satellite precipitation estimates with gauge adjustments and can generate rational precipitation estimates at fine spatiotemporal scales (0.25 • × 0.25 • and 3-hourly temporal intervals) within the global scope (60 • N-60 • S) [21,22]. As the successor of TRMM, the GPM constellation consists of a core observation platform and 10 cooperative satellites to observe global precipitation through inter-satellite cooperation. The GPM program offers three different levels of data products, all of which are available on the NASA website (http://www.nasa.gov/mission_pages/GPM). This study focuses on the GPM Level 3 products based on the IMERG algorithm. IMERG not only has a high spatial resolution (0.1 • × 0.1 • , 0.5 h) and fully global coverage (60 • S~60 • N), but also predicts flood disasters by reducing the uncertainty associated with short-term precipitation accumulation. At present, IMERG data has been developed from the first version to the sixth version, and IMERG-V06B is the latest version of the satellite rainfall data. However, Mohammad et al. (2019) verified that IMERG V05 performs better than V06, and therefore this study has employed IMERG V05 to estimate the precipitation data [23].
IMERG has three different products: Early, Late, and Final. In the real-time phase, the IMERG data generation system generates Early products after running once and then generates Late products after running again. The main difference between them is that the Early product is generated only by the forward propagation algorithm in the cloud mobile vector propagation algorithm, while the Late product is added with the backpropagation algorithm. The time-lag product Final introduces more sensor data sources based on the Late product. Meanwhile, the time of the extracted satellite data is matched with the station's data (08:00 a.m.) [24]. Table 1 shows the detailed information of the three IMERG products, with a study period of 2015-2018. Since IMERG-E and IMERG-L products are NRT products, they are released after 4 h and 12 h after observation, respectively; in turn, the IMERG-F is a PRT product, which has been calibrated for the deviation of the ground rainfall station, so it has a high accuracy and is usually released after two months of observation. Therefore, this study has adopted CMA data as calibration data, pays attention to the accuracy of the IMERG rainfall data to capture flash flood disasters, and combines the precipitation data estimated by IMERG-E and IMERG-F to discuss the accuracy of capturing flash flood, and thereby contributes to obtaining the static early warning thresholds.

Ground Observation Precipitation Data
The regional hourly precipitation integration products of the CMA (http://cdc.cma.gov.cn) is a reference product, which is based on the hourly precipitation observed by more than 30,000 automatic weather stations nationwide and the Climate Prediction Center morphing (CMORPH) satellite. Using the probability density function and optimal interpolation method, the precipitation fusion products with a high spatiotemporal resolution (0.1 • × 0.1 • and 1-h interval) have been generated in China. Among them, the optimal interpolation method is to combine the semi-hourly infrared data of the Geostationary Earth Orbit (GEO) satellite to interpolate the Passive microwave (PMW) inversion data, and then obtain the relatively fine precipitation. The gauge data involved in the CMA are subject to extreme quality control, including consistency checks, both internal and spatiotemporal. Meanwhile, the complex terrain, the low economic level in remote mountainous areas, and the sparse distribution of the ground stations have led to large errors in the measured data. The CMA data has a wide coverage and high spatiotemporal resolution; the overall error of this product is within 10%, and the error for heavy precipitation and site sparse areas is within 20%, which is more accurate than other similar products in China. Therefore, CMA products are suitable as the calibration data for product evaluation, but their scale should be consistent with the IMERG product in the calibration process [25]. Besides, when using CMA data to estimate the rainfall distribution, we should consider that the dataset may have a delayed phenomenon since it is actually a fusion of measured data and CMORPH satellite rainfall data.
The rainfall level is defined by the CMA, for example, the intraday rainfall is categorized into the rain and light-moderate rain (0-25 mm), heavy rain (25-50 mm), and rainstorm (>50 mm). Meanwhile, the minimum amount of precipitation that the gauges can measure is 0.1 mm. The critical rainfall in Yunnan Province mainly comes from the China rainstorm parameters atlas, which describes the statistical characteristics and laws of China's rainstorms. Based on the study period 2015-2018, the flash flood data were mainly from officially published data. In this study, a total of 120 flash flood events were analyzed, all of which resulted in deaths or missing of people. Additionally, some other data were also employed; for example, the topographic data came from a 1:250,000 digital elevation model (DEM) of the Yunnan Province.

RTI Method
Since 2005, the RTI model has been widely applied in flash flood warnings, which effectively reduced the casualties caused by flash floods. Considering that flash floods are mainly caused by hourly peak rainfall, this method is used to predict flash flood by multiplying the effective cumulative rainfall (R t ) and rainfall intensity (I). Among them, the effective accumulated rainfall mainly uses the accumulated rainfall in the previous 7 days before the flash flood occurs. Based on historical rainfall, this method obtains the flash flood probability under different rainfall conditions by calculating the RTI, and then divides the critical rainfall map into three regions for early warning (low possible, medium, and high) [13]. Meanwhile, it also defines how to segment the rainfall events, i.e., the start time is defined as the rainfall per hour exceeding 4 mm, and the end time is the rainfall per hour falling below 4 mm for six consecutive hours. A single day refers to the period from yesterday at 08:00 a.m. to today at 08:00 a.m. Rainfall data were selected according to the flash flood disaster events. If the flash flood disaster occurs in an area without a monitoring station, calculations are performed using data from the nearest rainfall stations (within a radius of 50 mm). The RTI equations are as follows: where RTI t is the RTI calculated at time t; I is the rainfall intensity; R t is the effective cumulative rainfall; i means the antecedent day number from one to n; α is the rainfall attenuation coefficient, mainly taken from the measured value of 0.78 by Cui Peng in Jiangjiagou, Yunnan Province [26]; α i is the reduction factor of the previous i day; and R i is the 24-h cumulative rainfall of the previous i-day, where the initial cumulative rainfall R 0 is 50 mm, which is mainly obtained through actual statistical analysis. Since the previous rainfall was calculated using "t" days of accumulated rainfall, the false alarms rate is higher in certain rainfall patterns (e.g., long-term duration and low rainfall intensity). Besides, this method does not consider the effects of intermittent rainfall and rainfall segmentation, etc.; all of the above results in low accuracy. Therefore, based on the reduction period and the reduction coefficient being unchanged, Chen et al. (2018) proposed an improved RTI method and provides a detailed flowchart describing the method [16]. The specific equation of this method is where I t is the current rainfall intensity at time t (mm/h); and R t−1 is the effective cumulative rainfall one hour earlier. In Formula (2), each operation needs to calculate R i separately, multiply it by α i , and then accumulate; however, in Formula (3), it has only one rainfall intensity and one accumulated rainfall, which greatly reduces the calculation amount and contributes to the future subsequent large-scale grid operations.

Evaluation Metrics
Six indicators were employed to evaluate the accuracy of the satellite precipitation products. Table 2 shows the formulas and optimal values of these indexes. CC (correlation coefficient) represents the correlation between the satellite precipitation data and site precipitation data, the greater the better; BIAS (Relative Bias) and RMSE (Root Mean Square Error) are quantitative indicators, representing the deviation degree between the satellite precipitation data and ground station precipitation data. The method also includes three classification indicators: Probability Of Detection (POD), False Alarm Ratio (FAR), and Critical Success Index (CSI), which can comprehensively reflect the ability of the precipitation products to estimate the probability of precipitation event occurrence. Among them, POD is the correct forecast rate, while FAR is the wrong forecast rate. The critical success index (CSI) is the function of POD and FAR that provides a more balanced estimate of the satellite products [27]. Table 2. List of the formulas and optimal values of the indexes.

Systematic or Random Error
Random error is an avoidable error caused by measurement or calculation. Systematic error is the unavoidable error caused by the experimental instrument or accuracy. Systematic and random errors are determined by objective and subjective factors, respectively, both of which can be reduced but not eliminated [28]. The spatiotemporal variability of precipitation, measurement error, and the uncertainty of sampling affects the satellite precipitation data accuracy, where the uncertainties include systematic and random errors (hereafter Syst and Rand), mainly from (1) sensor observation; (2) the algorithm used to estimate rainfall; and (3) sampling error. Reference [29] developed a method for estimating Syst and Rand for satellite precipitation products. The system of mean square error (MSE) and the formula for random error are where x is the satellite precipitation; y is the CMA precipitation; and n is the time steps number (here days); the formula for calculatingx is as follows: where a and b are the slope and intercept parameters that need to be calibrated, respectively. The calculation formula is Remote Sens. 2020, 12, 1954 8 of 15 Figure 2 shows the spatial distributions of the daily average precipitation in Yunnan Province captured by IMERG-E, IMERG-F, and CMA from March 2015 to December 2018. Using CMA data as the reference calibration data, the rainfall distribution in Yunnan Province is shown in Figure 2a. In Yunnan Province, the precipitation increases from northeast to southwest, with the largest precipitation occurring in the western border region. Referring to Figure 1, Yunnan's terrain is complex and changeable, with a high northwest and low southeast, descending stepwise from north to south. Moreover, flash floods are mainly concentrated in the southwest region and densely distributed in parts of the northeast, reflecting that the flash flood disaster is mainly affected by multiple factors with heavy rainfall as the main trigger. Meanwhile, there is a clear precipitation zone in the western region, with the smallest precipitation in the northwest and more precipitation in the central-western region. Figure 2b,c presents the distribution of the estimated precipitation in Yunnan Province using IMERG-E and IMERG-F, respectively. The overall trend is consistent with the precipitation distribution measured by CMA; that is, the estimated precipitation is relatively large in the southwest area, while low in the northeast relatively. Among them, the estimated rainfall of IMERG-E is less in the northeast region and has a larger coverage area, which is lower than that of CMA. The estimated maximum rainfall area is consistent with the CMA, but its coverage area is much smaller than the CMA. Therefore, IMERG-E underestimates the precipitation, especially in the southwest region. IMERG-F is the opposite of IMERG-E; its estimated rainfall covers a larger area, and the minimum and maximum rainfall are higher than the CMA. IMERG-F overestimates precipitation, especially in the central and western regions.  Figure 2 shows the spatial distributions of the daily average precipitation in Yunnan Province captured by IMERG-E, IMERG-F, and CMA from March 2015 to December 2018. Using CMA data as the reference calibration data, the rainfall distribution in Yunnan Province is shown in Figure 2a. In Yunnan Province, the precipitation increases from northeast to southwest, with the largest precipitation occurring in the western border region. Referring to Figure 1, Yunnan's terrain is complex and changeable, with a high northwest and low southeast, descending stepwise from north to south. Moreover, flash floods are mainly concentrated in the southwest region and densely distributed in parts of the northeast, reflecting that the flash flood disaster is mainly affected by multiple factors with heavy rainfall as the main trigger. Meanwhile, there is a clear precipitation zone in the western region, with the smallest precipitation in the northwest and more precipitation in the central-western region. Figure 2b, c presents the distribution of the estimated precipitation in Yunnan Province using IMERG-E and IMERG-F, respectively. The overall trend is consistent with the precipitation distribution measured by CMA; that is, the estimated precipitation is relatively large in the southwest area, while low in the northeast relatively. Among them, the estimated rainfall of IMERG-E is less in the northeast region and has a larger coverage area, which is lower than that of CMA. The estimated maximum rainfall area is consistent with the CMA, but its coverage area is much smaller than the CMA. Therefore, IMERG-E underestimates the precipitation, especially in the southwest region. IMERG-F is the opposite of IMERG-E; its estimated rainfall covers a larger area, and the minimum and maximum rainfall are higher than the CMA. IMERG-F overestimates precipitation, especially in the central and western regions.  Figure 3 presents the spatial patterns of six indicators obtained from IMERG-E and the CMA using hourly data over the Yunnan Province. Generally, CC was relatively low, changing from 0.2 to 0.5, especially over the northwest regions (CC < 0.2) with the highest altitudes. As for the RMSE and BIAS, their spatial patterns are similar, with an increasing trend for the RMSE and BIAS and a decreasing trend for relative error from northwest to southeast. Many previous studies have confirmed that satellite precipitation estimates are usually low with large errors in mountainous areas. For example, compared with the highest regions in the Northwest, IMERG-E's BIAS is underestimated by about 50% [30]. As for satellite-based precipitation metrics, there is also a significant trend; that is, POD and CSI are increasing, while FAR is decreasing, which is in harmony with the overall trend of the IMERG-E data. Table 3 shows the evaluation metrics for IMERG calculated with the mean value at hourly and daily timescales; the evaluation index calculated by the  Figure 3 presents the spatial patterns of six indicators obtained from IMERG-E and the CMA using hourly data over the Yunnan Province. Generally, CC was relatively low, changing from 0.2 to 0.5, especially over the northwest regions (CC < 0.2) with the highest altitudes. As for the RMSE and BIAS, their spatial patterns are similar, with an increasing trend for the RMSE and BIAS and a decreasing trend for relative error from northwest to southeast. Many previous studies have confirmed that satellite precipitation estimates are usually low with large errors in mountainous areas. For example, compared with the highest regions in the Northwest, IMERG-E's BIAS is underestimated by about 50% [30]. As for satellite-based precipitation metrics, there is also a significant trend; that is, POD and CSI are increasing, while FAR is decreasing, which is in harmony with the overall trend of the IMERG-E data. Table 3 shows the evaluation metrics for IMERG calculated with the mean value at hourly and daily timescales; the evaluation index calculated by the mean is larger on a daily scale than that on an hourly scale. Given the relatively low accuracy of IMERG-E, the differences between the IMERG-E and CMA were relatively obvious.

Evaluation of IMERG-E and IMERG-F
Remote Sens. 2020, 12, x FOR PEER REVIEW 9 of 16 mean is larger on a daily scale than that on an hourly scale. Given the relatively low accuracy of IMERG-E, the differences between the IMERG-E and CMA were relatively obvious. , and CSI(f)) generated from IMERG-E and CMA using hourly data in Yunnan Province. The spatial patterns of the six metrics derived from IMERG-F and CMA ( Figure 4) exhibited similar trends as those from IMERG-E and CMA. Table 2 revealed the four indicators CC, BIAS, POD, and CSI of NRL IMERG-E products are lower than these for PRL IMERG-F, regardless of daily and hourly data, while the two indicators RMSE and FAR are higher than that for IMERG-F. Therefore, IMERG-F indicated a better performance in this region. For instance, the POD and CSI of IMERG-F and CMA were also overall higher than those of IMERG-F and CMA, especially for CSI. Besides, in both Figures 3 and 4, the POD and CSI in these two figures show significant differences inside and outside of the Yunnan Province, which is mainly due to the complex and changeable terrain that induces large systematic errors in satellite precipitation. Besides, there is no actual rainfall measurement site outside the Yunnan border, and the rainfall distribution obtained only by interpolation appears discontinuous.  The spatial patterns of the six metrics derived from IMERG-F and CMA ( Figure 4) exhibited similar trends as those from IMERG-E and CMA. Table 2 revealed the four indicators CC, BIAS, POD, and CSI of NRL IMERG-E products are lower than these for PRL IMERG-F, regardless of daily and hourly data, while the two indicators RMSE and FAR are higher than that for IMERG-F. Therefore, IMERG-F indicated a better performance in this region. For instance, the POD and CSI of IMERG-F and CMA were also overall higher than those of IMERG-F and CMA, especially for CSI. Besides, in both Figures 3 and 4, the POD and CSI in these two figures show significant differences inside and outside of the Yunnan Province, which is mainly due to the complex and changeable terrain that induces large systematic errors in satellite precipitation. Besides, there is no actual rainfall measurement site outside the Yunnan border, and the rainfall distribution obtained only by interpolation appears discontinuous.
The following is a further analysis of the error source (i.e., systematic error or random error) of EMERG-E. The spatial patterns of systematic and random errors in IMERG-E at the 1 and 24 h temporal scales in Yunnan Province, respectively, demonstrated almost the same spatial patterns ( Figure 5). Overall, the system error is mainly distributed in the northwest region of Yunnan Province with a high altitude (about 80%), and the coverage area is relatively small. Meanwhile, except for the northwestern part of Yunnan Province, the systematic error at 1 h is significantly higher than that at 24 h, especially the error at 24 h mostly disappeared. Random error mainly occurs in the southern region of Yunnan Province with relatively high rainfall (about 80%), but a relatively low in the northwest region. Besides, the random error is higher at 24 h than at 1 h. Therefore, for the IMERG-E data, the estimated error in the high-altitude and low-precipitation areas of Yunnan Province is mainly determined by the systematic error; the southern area of Yunnan Province with a low elevation and high precipitation is dominated by random errors. The following is a further analysis of the error source (i.e., systematic error or random error) of EMERG-E. The spatial patterns of systematic and random errors in IMERG-E at the 1 and 24 h temporal scales in Yunnan Province, respectively, demonstrated almost the same spatial patterns ( Figure 5). Overall, the system error is mainly distributed in the northwest region of Yunnan Province with a high altitude (about 80%), and the coverage area is relatively small. Meanwhile, except for the northwestern part of Yunnan Province, the systematic error at 1 h is significantly higher than that at 24 h, especially the error at 24 h mostly disappeared. Random error mainly occurs in the southern region of Yunnan Province with relatively high rainfall (about 80%), but a relatively low in the northwest region. Besides, the random error is higher at 24 h than at 1 h. Therefore, for the IMERG-E data, the estimated error in the high-altitude and low-precipitation areas of Yunnan Province is mainly determined by the systematic error; the southern area of Yunnan Province with a low elevation and high precipitation is dominated by random errors. As for IMERG-F, systematic errors accounted for less than 20% of the total errors across the whole study area at both the 1 and 24 h temporal scales. Compared with IMERG-E, IMERG-F has been significantly improved because the systematic errors were effectively reduced, especially in mountainous areas with the complex terrain. Figure 6 shows that the IMERG-F data error (≥80%) is primarily a random error, independent of the underlying surface features. Besides, Figures 6 and 7 behave differently inside and outside the Yunnan border. This phenomenon is mainly due to the complex and changeable terrain of Yunnan Province, which further triggers the discontinuous As for IMERG-F, systematic errors accounted for less than 20% of the total errors across the whole study area at both the 1 and 24 h temporal scales. Compared with IMERG-E, IMERG-F has been significantly improved because the systematic errors were effectively reduced, especially in mountainous areas with the complex terrain. Figure 6 shows that the IMERG-F data error (≥80%) is primarily a random error, independent of the underlying surface features. Besides, Figures 6 and 7 behave differently inside and outside the Yunnan border. This phenomenon is mainly due to the complex and changeable terrain of Yunnan Province, which further triggers the discontinuous changes in satellite precipitation errors. It is a discontinuous problem in the original satellite's precipitation error, rather than using different spatiotemporal resolution rainfall. Therefore, the overall IMERG-F data is better than the IMERG-E data.

Applicability Analysis of IMERG in Flash Flood Warning
Based on historical flash flood disaster events, combined with three types of rainfall products (IMERG-E and -F, and CMA), the multi-period rainfall (1 h, 3 h, 6 h, and 24 h) is obtained. The resultant effective accumulated multi-period rainfall is calculated by the improved RTI method. Meanwhile, combined with the flash floods' actual frequency and the rainstorm statistical parameter atlas of China, the multi-period critical rainfall (1 h, 3 h, 6 h, and 24 h) (hereafter, CR1, CR3, CR6, and CR24) is obtained; the G (x) early warning model is constructed for effective cumulative rainfall (Rt) and corresponding period critical rainfall (CRt). Since this model does not take into account the potential multiple influencing factors, such as slope, vegetation, and human activity, when issuing the flash flood warning, we should make a comprehensive analysis regarding the flash flood risk distribution to determine whether a flash flood event has been captured. Considering the flash flood risk map obtained by Ma et al. (2019), the obtained results are shown in Figures 7 and 8 [31].
Obviously, for the same product, the captured flash flood events has decreased over time, where the hit rate of both CR6 and CR24 is less than 60%. Among them, the flash floods events captured by these three rainfall products are concentrated in western Yunnan, while there are fewer flash flood events in the relatively flat southwest and central Yunnan. In general, if a flash flood event cannot be captured in a short period, it will not be anyhow captured. Moreover, the effect of IMERG-E products

Applicability Analysis of IMERG in Flash Flood Warning
Based on historical flash flood disaster events, combined with three types of rainfall products (IMERG-E and -F, and CMA), the multi-period rainfall (1 h, 3 h, 6 h, and 24 h) is obtained. The resultant effective accumulated multi-period rainfall is calculated by the improved RTI method. Meanwhile, combined with the flash floods' actual frequency and the rainstorm statistical parameter atlas of China, the multi-period critical rainfall (1 h, 3 h, 6 h, and 24 h) (hereafter, CR1, CR3, CR6, and CR24) is obtained; the G (x) early warning model is constructed for effective cumulative rainfall (R t ) and corresponding period critical rainfall (CRt). Since this model does not take into account the potential multiple influencing factors, such as slope, vegetation, and human activity, when issuing the flash flood warning, we should make a comprehensive analysis regarding the flash flood risk distribution to determine whether a flash flood event has been captured. Considering the flash flood risk map obtained by Ma      Obviously, for the same product, the captured flash flood events has decreased over time, where the hit rate of both CR6 and CR24 is less than 60%. Among them, the flash floods events captured by these three rainfall products are concentrated in western Yunnan, while there are fewer flash flood events in the relatively flat southwest and central Yunnan. In general, if a flash flood event cannot be captured in a short period, it will not be anyhow captured. Moreover, the effect of IMERG-E products on capturing disasters is significantly lower than that of IMERG-F, and the capture rate in each period is less than 50%. However, the capture accuracy of IMERG-F at CR1 and CR3 is almost comparable to the CMA-only a 1% difference-and the capture rate of CR1 is nearly 80%, with an extremely high accuracy. Nonetheless, with the time marching, the capture effect decreased significantly. Meanwhile, if the flash flood events cannot be captured by IMERG-F and CMA, the same phenomenon occurs in IMERG-F; but, as time goes on, there is an out-of-sync phenomenon in the catching flash flood events by IMERG-F and CMA. Besides, the CMA data is a fusion of measured data and CMORPH, leading to possible delays in data acquisition. Therefore, the hit rate is decreasing with an increasing temporal resolution or increasing the averaging time of the satellite images.

Conclusions
This study first quantitatively evaluated the IMERG-E and IMERG-F satellite precipitation products, and then combined the improved RTI method to analyze its application effects in flash flood warnings, which has three major contributions: (i) Flash flood warning aspects are integrated, for the first time, with satellite precipitation to account for the applicability of satellite data in flash flood warnings. The result shows that the early warning effect of IMERG-F products is better at the 1 h and 3 h scale than that at the daily scale; (ii) the study area has not been documented in previous studies. Yunnan Province is characterized by a low latitude but high altitude where satellite precipitation exhibited some new characteristics, including that precipitation in Yunnan Province has increased from northeast to southwest, where the largest precipitation occurred in the flood-prone area in the southwest part; (iii) this study reveals some interesting phenomena that were not reported in related research [32].
For example, the systematic error of IMERG-E is mainly distributed in areas with high altitudes and low precipitation, and the random error is mainly distributed in areas with low altitudes and high precipitation. The most important thing is that as for the same satellite rainfall product, the flash flood disaster events that can be captured decreases with time. For different satellite rainfall products, the flash flood events captured by the IMERG-E products are significantly lower than IMERG-F. Meanwhile, the capture rate of each period is less than 50%.
To sum up, the above research results will contribute to the application of satellite remote sensing data in flash flood warnings.