The First Comparisons of IMERG and the Downscaled Results Based on IMERG in Hydrological Utility over the Ganjiang River Basin

: Rainfall information is a prerequisite to and plays a vital role in driving hydrological models. However, limited by the observation methods, the obtained precipitation data, at present, are still too coarse. In this study, a new downscaling method was proposed to obtain high spatial resolution (~1 km/hourly) precipitation estimates based on Integrated Multi-satellitE Retrievals for GPM (IMERG) data at hourly scale. Compared with original IMERG data, the downscaled precipitation results showed the similar spatial patterns with those of original IMERG data, but with ﬁner spatial resolution. In addition, the downscaled precipitation estimates were further analyzed to quantify their improvements using the Coupled Routing and Excess STorage (CREST) model across Ganjiang River basin. Compared with the observed streamﬂow, the downscaled precipitation results showed satisfying hydrological performance, with Nash-Sutcliffe Coefﬁcient of Efﬁciency (NSCE), Root Mean Square Error (RMSE), Relative Bias (BIAS), and Correlation Coefﬁcient (CC). The improvement in terms of four statistic metrics in terms of streamﬂow simulation also indicated great potential of hydrological utility for the downscaled precipitation results.


Introduction
The frequency and intensity of extreme rainfall events and severe storms present an increasing trend during recent years, resulting in disastrous floods and landslides, which pose great threats to personal security and the economy [1,2]. Distributed hydrologic models play irreplaceable roles in studying various hydrological variables and important hydrologic process at various scales and thus provide helpful information for water resources management, and disaster prediction (e.g., floods) [3,4].
(0.01 • × 0.01 • , hourly) used to downscale the IMERG data; and (2) far few rainfall pixels used for traditional regressing relationships between precipitation and environmental variables. In this case, most traditional or commonly used regression/empirical models and the state-of-art machine learning algorithms (such as Regression, Machine learning, Decision tree, Bayesian, Artificial neural network, Deep learning, Ensemble, Support vectors) are not suitable.
Motivated by the critical requirement to downscale IMERG data at hourly scale to obtain downscaled precipitation estimates with good quality, this study, specifically, proposed a geographically moving window technique to obtain downscaled precipitation estimates at hourly scale based on the IMERG V4 data, and then, both the original IMERG data, and the downscaled product were used to drive the distributed hydrological model to explore their capability in hydrological modeling.

Study Area
As the seventh largest branch of the Yangtze River, Ganjiang River basin is located within 113 • 30 -116 • 40 E and 24 • 29 -29 • 21 N (Figure 1), Jiangxi province, China. The total drainage area of the Ganjiang River basin is 81,158 km 2 above the Waizhou hydrologic station, which is the outlet of the basin. The topography of this basin is complex with elevation varying from 12 to 1930 m. The northwestern and eastern Ganjiang River basin are bordered by Jiuling mountain and Wuyi mountain, respectively.
With moderate climate and sufficient rainfall, the Ganjiang River basin is one of the typical rainstorm regions in China. The average annual rainfall is around 1500 mm. Over 70% of the total rainfall is happening in the plum rain season, from April to June. The maximum precipitation usually happens in the northwest and east of the Ganjiang River basin during the rainy season. Affected by both of the complex terrain and the East Asian Monsoon, this region is a historically flood-prone region and receives frequently flood disasters.
Water 2018, 10, x FOR PEER REVIEW 3 of 15 variables (0.01° × 0.01°, hourly) used to downscale the IMERG data; and (2) far few rainfall pixels used for traditional regressing relationships between precipitation and environmental variables. In this case, most traditional or commonly used regression/empirical models and the state-of-art machine learning algorithms (such as Regression, Machine learning, Decision tree, Bayesian, Artificial neural network, Deep learning, Ensemble, Support vectors) are not suitable. Motivated by the critical requirement to downscale IMERG data at hourly scale to obtain downscaled precipitation estimates with good quality, this study, specifically, proposed a geographically moving window technique to obtain downscaled precipitation estimates at hourly scale based on the IMERG V4 data, and then, both the original IMERG data, and the downscaled product were used to drive the distributed hydrological model to explore their capability in hydrological modeling.

Study Area
As the seventh largest branch of the Yangtze River, Ganjiang River basin is located within 113°30′-116°40′ E and 24°29′-29°21′ N (Figure 1), Jiangxi province, China. The total drainage area of the Ganjiang River basin is 81,158 km 2 above the Waizhou hydrologic station, which is the outlet of the basin. The topography of this basin is complex with elevation varying from 12 to 1930 m. The northwestern and eastern Ganjiang River basin are bordered by Jiuling mountain and Wuyi mountain, respectively.
With moderate climate and sufficient rainfall, the Ganjiang River basin is one of the typical rainstorm regions in China. The average annual rainfall is around 1500 mm. Over 70% of the total rainfall is happening in the plum rain season, from April to June. The maximum precipitation usually happens in the northwest and east of the Ganjiang River basin during the rainy season. Affected by both of the complex terrain and the East Asian Monsoon, this region is a historically flood-prone region and receives frequently flood disasters.

The Hydrological Stations
As shown in Figure 1, hourly streamflow data are available at five hydrological stations, Waizhou, Baisha, Dongbei, Xintian, and Saitang, during the period from 1 May to 30 September 2014.

IMERG V4 Data
The GPM mission is an international network of satellites that provide the next-generation global observations of rain and snow, building upon the success of TRMM. The Integrated Multi-satellitE Retrievals for GPM (IMERG) product was created by inter-calibrating, merging and interpolating all satellite microwave precipitation estimates, together with microwave-calibrated infrared (IR) satellite estimates, precipitation gauge analyses, and potentially other precipitation estimators at finer scales from the TRMM and GPM eras over the entire globe. The IMERG system has three products: the "Early" run (4 h after observation), the "Late" run (12 h after observation) multi-satellite products, and the "Final" run (two months after the observation month) satellite-gauge product with best accuracy and quality starting from March 2014 [24]. The IMERG "Final" run product (0.1 • /half-hourly) was used in this study, which was then accumulated to hourly resolution.

Topographical Characteristic
We collected DEM data from the Shuttle Radar Topography Mission (SRTM), which is operated by the National Geospatial-Intelligence Agency (NGA) and NASA, and provides high resolution DEMs between 56 • S and 60 • N globally [39]. DEM data were downloaded over the Ganjiang River basin at an original spatial resolution of 90 m, which was aggregated into results with 1 km.

Geographically Moving Window Weight Disaggregation Analysis
Li et al. [8] pointed out that the maximum rainfall usually occurs over the northwestern (Jiuling mountian) and eastern (Wuyi mountain) Ganjiang River basin, which means topography affects the spatial distribution of rainfall (e.g., orographic rain) over this basin. As well as to cater the two great challenges, aforementioned above, of downscaling IMERG data at hourly scale, we proposed a new algorithm called Geographically Moving Window Weight Disaggregation Analysis (GMWWDA) to obtain the disaggregation weights based on DEM. In addition, the downscaled results would be obtained by disaggregating the IMERG data using the disaggregation weights (1 km).
The steps to obtain disaggregation weights using DEM based on the assumption that the volume of the precipitation is inversely proportional to DEM, which means there is a general trend that as the elevation rising, the volume of rainfall decreases [16,40]. Using the GMWWDA method, the size of the moving window was the first key parameter needed to be determined. In this study, we have tried a series of window sizes, for example, 3 by 3, 4 by 4, 5 by 5 etc. In addition, then the corresponding downscaled results were used to drive the hydrological model, while the downscaled results using other window sizes performed similar or even worse than the window size of 3 by 3. Therefore, the window size of 3 by 3 has been selected in this study, and illustrated in the downscaling procedures.
(i) Firstly, we obtained the weights of all the pixels in the current moving window, 3 by 3, using Equation (1), (ii) obtained the rank of central pixel according to the weights in the current window from high to low by Equation (1), shown as Equation (2), where the function sort () was adopted from the R language (http://www.r-project.org/). (iii) the disaggregation weight based on DEM at the central pixel of the current window was obtained by Equation (3), according to the general opposite trend between DEM values and the rainfall volumes, where Rank 5 represents the rank of the central pixel shown in Figure 2, while N is the window size, and W is the weight from Equation (1).
where 5 represents the rank of the central pixel shown in Figure 2, while N is the window size, and W is the weight from Equation (1).

Main Downscaling Steps by Introducing GMWWDA
To downscale the coarse satellite-based precipitation estimates, IMERG, at hourly scale, GMWWDA was introduced as the core algorithm in this study, and the main sequential steps in the downscaling process were shown as below.
(1) The 0.0083° arc DEM data were obtained and aggregated to 0.01° in raster format.
(2) The disaggregation weights based on DEM were calculated according to the algorithms described in Section 3.1 and were constant at the temporal series. (3) The disaggregation weights based on DEM were transformed from raster format (0.01°) into points, which were then used to extract precipitation information from IMERG data in raster (0.1°). Finally, the downscaled results based on only DEM were obtained by multiplying IMERG value at 0.01° spatial resolution with the corresponding disaggregation weights.

The CREST Hydrologic Model
The CREST hydrologic model is a hybrid modeling strategy, which was jointly developed by the Oklahoma and the NASA SERVIR Project Team [41]. The CREST model simulates the water cycle on a regular grid which could be defined by the user, and therefore it can be applied at various scales. CREST has been used in the Flooded Locations and Simulated Hydrographs Project (FLASH) and the Near Real-time Global hydrological Simulation and Flood Monitoring Demonstration System to provide online flood predictions and multi-scale basin applications. The CREST model also has many successful applications in many basins [40,42,43]. Kan et al. [28] developed an updated CREST model Value (2) Value (3) Value (5) Value (6) 0.01°V alue (1) Value (4) Value (8) Value(9) Value (7) 0.03°F igure 2. Moving window and pixels in and overlapping it, at hourly temporal scale.

Main Downscaling Steps by Introducing GMWWDA
To downscale the coarse satellite-based precipitation estimates, IMERG, at hourly scale, GMWWDA was introduced as the core algorithm in this study, and the main sequential steps in the downscaling process were shown as below.
(1) The 0.0083 • arc DEM data were obtained and aggregated to 0.01 • in raster format.
(2) The disaggregation weights based on DEM were calculated according to the algorithms described in Section 3.1 and were constant at the temporal series. (3) The disaggregation weights based on DEM were transformed from raster format (0.01 • ) into points, which were then used to extract precipitation information from IMERG data in raster (0.1 • ). Finally, the downscaled results based on only DEM were obtained by multiplying IMERG value at 0.01 • spatial resolution with the corresponding disaggregation weights.

The CREST Hydrologic Model
The CREST hydrologic model is a hybrid modeling strategy, which was jointly developed by the Oklahoma and the NASA SERVIR Project Team [41]. The CREST model simulates the water cycle on a regular grid which could be defined by the user, and therefore it can be applied at various scales. CREST has been used in the Flooded Locations and Simulated Hydrographs Project (FLASH) and the Near Real-time Global hydrological Simulation and Flood Monitoring Demonstration System to provide online flood predictions and multi-scale basin applications. The CREST model also has many successful applications in many basins [40,42,43]. Kan et al. [28] developed an updated CREST model version 3.0, which mainly differs version 2.0 in following aspects: (1) switching from one soil layer to 3 soil layers, and adding 3-soil-layer tension water soil moisture and evapotranspiration computation; (2) using free water distribution curve to describe the sub-grid variations of free water storage; and (3) switching the flow concentration module into a 4-mechanisim-based cell-to-cell routing scheme, including overland flow, interflow, ground water, and open channel flow. In this study, we used CREST v3.0 to evaluate the hydrological performance of the IMERG products and downscaled results based on the IMERG data. The CREST model simulation was conducted at Ganjiang River basin, and other forcing data is similar with those in the study by Kan et al. [7]. The hourly observed streamflow data at five hydrological stations, namely Waizhou, Baisha, Dongbei, Xintian, and Saitang, during the period from 1 May to 30 September 2014 were used to do hydrological analysis.

Hydrological Utility of Downscaled Results
Both the original IMERG data (10 km/1 hourly) the downscaled results (1 km/1 hourly) in the period from 1 May to 30 September 2014, were used to force the hydrological model, CREST, to simulate the streamflow at the five hydrological stations. In addition, the four diagnostic statistics were used to assess the performances of the original IMERG data and the downscaled results.

Diagnostic Statistics
Several popularly used diagnostic statistics, correlation coefficient (CC), root mean square error (RMSE), mean absolute error (MAE), and Nash-Sutcliffe coefficient efficiency (NSCE) [8], were adopted in this study to quantitatively evaluate the quality of the IMERG product and downscaled results based on IMERG, at hourly temporal scale, in simulating the streamflow at five hydrological stations.

The Relationship between the IMERG and DEM over Ganjiang
To downscale the IMERG data at basin scale providing precipitation inputs with finer spatial resolution as well as quality, the relationship between precipitation and auxiliary data was critical and needed to be investigated firstly in the specific regions. In this study, the spatial patterns of mean hourly precipitation estimates measured by IMERG and DEM, with the spatial resolution of 0.1 • , were shown in Figure 3. From southwest to northeast, there was an increasing trend of mean hourly precipitation estimates captured by IMERG, while the DEM demonstrated a decreasing trend. The spatial patterns of the precipitation and DEM exhibited an intuitively negative relationship.
There were 825 pixels of both IMERG and DEM at the resolution of 0.1 • , over the Ganjiang basin. We explored the linear relationship between IMERG and DEM at this spatial scale in Figure 4. Similarly with those findings in Figure 3, the quantitatively relationship between IMERG and DEM was also negative. Over Ganjing basin, as the elevation rose 1 km, the precipitation decreased~0.075 mm/hour, which also passed the significant test with p value less than 2 × 10 −16 . The relationship revealed from Figures 3 and 4 explain the reasonability of the downscaling assumption that the volume of the precipitation is inversely proportional to DEM. There were 825 pixels of both IMERG and DEM at the resolution of 0.1°, over the Ganjiang basin. We explored the linear relationship between IMERG and DEM at this spatial scale in Figure 4. Similarly with those findings in Figure 3, the quantitatively relationship between IMERG and DEM was also negative. Over Ganjing basin, as the elevation rose 1 km, the precipitation decreased ~0.075 mm/hour, which also passed the significant test with p value less than 2 × 10 −16 . The relationship revealed from Figures 3 and 4 explain the reasonability of the downscaling assumption that the volume of the precipitation is inversely proportional to DEM.   There were 825 pixels of both IMERG and DEM at the resolution of 0.1°, over the Ganjiang basin. We explored the linear relationship between IMERG and DEM at this spatial scale in Figure 4. Similarly with those findings in Figure 3, the quantitatively relationship between IMERG and DEM was also negative. Over Ganjing basin, as the elevation rose 1 km, the precipitation decreased ~0.075 mm/hour, which also passed the significant test with p value less than 2 × 10 −16 . The relationship revealed from Figures 3 and 4 explain the reasonability of the downscaling assumption that the volume of the precipitation is inversely proportional to DEM.

Disaggregation Weights Based on DEM
To obtain the downscaled results, the disaggregation weight is a perquisite. According to the proposed GMWWDA in this study, the disaggregation weights have been obtained, shown in Figure 5.
Overall, the spatial patterns demonstrated some heterogeneous characteristics, for example, the values were relatively larger in the central Ganjiang River basin, where the DEM values were relatively smaller. The values of disaggregation weights vary from~0.13 to~3.34, with the average value~0.96 and standard error~0.27.
To obtain the downscaled results, the disaggregation weight is a perquisite. According to the proposed GMWWDA in this study, the disaggregation weights have been obtained, shown in Figure  5. Overall, the spatial patterns demonstrated some heterogeneous characteristics, for example, the values were relatively larger in the central Ganjiang River basin, where the DEM values were relatively smaller. The values of disaggregation weights vary from ~0.13 to ~3.34, with the average value ~0.96 and standard error ~0.27.

Downscaled Precipitation Results
After the disaggregation weights being obtained, the downscaled results were generated based on the IMERG, from 1 May to 30 September 2014. In most hours, there were very limited rainfall pixels, or no rainfall. Therefore, a typical hour with relatively large rainfall was selected to demonstrate the downscaled results. The spatial patterns of (a) original IMERG, (b) downscaled results, at 10 a.m. 3 July 2014 over the Ganjiang, were shown in Figure 6. Using the proposed downscaling algorithm by this study, the general spatial pattern of downscaled result was similar to that of the original IMERG data. One of the most important characteristics of the downscaled results was that there was still no rainfall over the region where the IMERG did not capture the rainfall happening. While the traditional regression method could not achieve this. At this typical hour, the value of IMERG ranged from 0 to 20.42 mm/hour, with average value ~1.18 mm/hour, and standard deviation ~2.05 mm/hour; while the value of downscaled IMERG ranged from 0 to 29.66 mm/hour, with average value ~1.14 mm/hour, and standard deviation ~2.03 mm/hour. By comparing the distributions of the values from both IMERG and the downscaled results, the volume magnitude of these original IMERG and the corresponding downscaled results were similar. However, the downscaled results could provide more detailed information with finer spatial resolutions.

Downscaled Precipitation Results
After the disaggregation weights being obtained, the downscaled results were generated based on the IMERG, from 1 May to 30 September 2014. In most hours, there were very limited rainfall pixels, or no rainfall. Therefore, a typical hour with relatively large rainfall was selected to demonstrate the downscaled results. The spatial patterns of (a) original IMERG, (b) downscaled results, at 10 a.m. 3 July 2014 over the Ganjiang, were shown in Figure 6. Using the proposed downscaling algorithm by this study, the general spatial pattern of downscaled result was similar to that of the original IMERG data. One of the most important characteristics of the downscaled results was that there was still no rainfall over the region where the IMERG did not capture the rainfall happening. While the traditional regression method could not achieve this. At this typical hour, the value of IMERG ranged from 0 to 20.42 mm/hour, with average value~1.18 mm/hour, and standard deviation~2.05 mm/hour; while the value of downscaled IMERG ranged from 0 to 29.66 mm/hour, with average value~1.14 mm/hour, and standard deviation~2.03 mm/hour. By comparing the distributions of the values from both IMERG and the downscaled results, the volume magnitude of these original IMERG and the corresponding downscaled results were similar. However, the downscaled results could provide more detailed information with finer spatial resolutions.

Validations against Ground Observations over the Tibetan Plateau
To validate the accuracy of downscaled results, this proposed downscaling method was also applied over the Tibetan Plateau (TP) where we have collected 258 hourly rain gauge data at automatic weather stations (AWSs) (shown in Figure 7), which were provided by the Chinese Meteorological Data Sharing Service System, and under quality control [44]. At the typical hour (11 a.m. 2 July 2014) at local scale in southeastern TP, the downscaled result demonstrated similar spatial trends with original IMERG trends, but with finer spatial resolution (~1 km), shown in Figure 7a

Validations against Ground Observations over the Tibetan Plateau
To validate the accuracy of downscaled results, this proposed downscaling method was also applied over the Tibetan Plateau (TP) where we have collected 258 hourly rain gauge data at automatic weather stations (AWSs) (shown in Figure 7), which were provided by the Chinese Meteorological Data Sharing Service System, and under quality control [44]. At the typical hour (11 a.m. 2 July 2014) at local scale in southeastern TP, the downscaled result demonstrated similar spatial trends with original IMERG trends, but with finer spatial resolution (~1 km), shown in Figure 7a

Validations against Ground Observations over the Tibetan Plateau
To validate the accuracy of downscaled results, this proposed downscaling method was also applied over the Tibetan Plateau (TP) where we have collected 258 hourly rain gauge data at automatic weather stations (AWSs) (shown in Figure 7), which were provided by the Chinese Meteorological Data Sharing Service System, and under quality control [44]. At the typical hour (11 a.m. 2 July 2014) at local scale in southeastern TP, the downscaled result demonstrated similar spatial trends with original IMERG trends, but with finer spatial resolution (~1 km), shown in Figure 7a  Both the hourly IMERG and downscaled results (in the period of July 2014) over the TP were validated against rain gauges, shown in Figure 8. Four commonly used indices were applied to evaluate the performances of IMERG and the corresponding downscaled results [15]. The original IMERG data has reasonable accuracy against ground observations with R 2~0 .51, MAE~0.72 mm/hour, RMSE~0.88 mm/hour, and Bias~56.77%, while the accuracy of the corresponding downscaled results have been significantly improved, with R 2~0 .69, MAE~0.46 mm/hour, RMSE~0.63 mm/hour, and Bias~13.81%, especially in terms of the Bias. Therefore, we inferred this downscaling algorithm was also suitable for Ganjiang basin. Both the hourly IMERG and downscaled results (in the period of July 2014) over the TP were validated against rain gauges, shown in Figure 8. Four commonly used indices were applied to evaluate the performances of IMERG and the corresponding downscaled results [15]. The original IMERG data has reasonable accuracy against ground observations with R 2~0 .51, MAE~0.72 mm/hour, RMSE~0.88 mm/hour, and Bias ~ 56.77%, while the accuracy of the corresponding downscaled results have been significantly improved, with R 2~0 .69, MAE~0.46 mm/hour, RMSE~0.63 mm/hour, and Bias~13.81%, especially in terms of the Bias. Therefore, we inferred this downscaling algorithm was also suitable for Ganjiang basin.

Comparisons of the Performances of Both IMERG and the Downscaled Results in CREST
Both the original IMERG data and the downscaled results were used to force the hydrological model, CREST, in Ganjiang River basin, from 1 May to 30 September 2014, and four diagnostic statistics have been used to assess the performance, shown in Figure 9. As for CC, the agreement of IMERG at Dongbei and Saitang was better than that of downscaled results, while at the other stations, the downscaled results demonstrated better agreements. Additionally, due to the updated CREST model, the CC of both the IMERG and downscaled in driving CREST to simulate the streamflow at five stations were very large, ranging from 0.7 to 0.90, and the differences were not very significant. In terms of MAE, and RMSE, the values of downscaled results were all smaller than those of original IMERG data in driving CREST to simulate the streamflow at each station. While for NSCE, the values of NSCE of downscaled results were larger than those of the original IMERG in forcing CREST model.
After spatially examining the performances of the both IMERG and downscaled results, we also found that at stations where the streamflow was relatively larger, such as the outlet Waizhou, and Dongbei, the models performed better (e.g., CC ~ 0.90, NSCE ~ 0.75) than those at stations where the streamflow was relatively smaller, such as Xintian, Saitang, and Baisha (e.g., CC ~ 0.77, NSCE ~ 0.50), which might be relevant with the CREST model.

Comparisons of the Performances of Both IMERG and the Downscaled Results in CREST
Both the original IMERG data and the downscaled results were used to force the hydrological model, CREST, in Ganjiang River basin, from 1 May to 30 September 2014, and four diagnostic statistics have been used to assess the performance, shown in Figure 9. As for CC, the agreement of IMERG at Dongbei and Saitang was better than that of downscaled results, while at the other stations, the downscaled results demonstrated better agreements. Additionally, due to the updated CREST model, the CC of both the IMERG and downscaled in driving CREST to simulate the streamflow at five stations were very large, ranging from 0.7 to 0.90, and the differences were not very significant. In terms of MAE, and RMSE, the values of downscaled results were all smaller than those of original IMERG data in driving CREST to simulate the streamflow at each station. While for NSCE, the values of NSCE of downscaled results were larger than those of the original IMERG in forcing CREST model.
After spatially examining the performances of the both IMERG and downscaled results, we also found that at stations where the streamflow was relatively larger, such as the outlet Waizhou, and Dongbei, the models performed better (e.g., CC~0.90, NSCE~0.75) than those at stations where the streamflow was relatively smaller, such as Xintian, Saitang, and Baisha (e.g., CC~0.77, NSCE~0.50), which might be relevant with the CREST model.

The Accuracy of the Downscaled Results Still Needs to Be Validated Spatially, Not Only at Limited Hydrological Stations
To meet the two challenges aforementioned above, a new downscaling algorithm proposed for obtaining downscaled results at finer tempo-spatial resolutions, especially over regions with limited ground-based rain gauges or radars. Though the downscaled results overall outperformed the original IMERG data in CREST simulating the streamflow at five hydrological stations, they were still needed to be validated using ground observations [45].

The Uncertainties and Sensitivities of the Parameters in This CREST Model
As pointed out by Kan et al. [28], there are commonly three approaches adopted to build a distributed hydrological model: the conceptually based models, the partial differential equation-based models, and conceptually physically coupled models. Even though many methods can be used to construct complex hydrological models, definitely precise description of the nature of hydrological processes is still a hard job. This is because the hydrological models inherently have some limitations, such as uncertainty issue, when facing with complex natural scenario simulations. Though there are numerous hydrological models developed and applied worldwide with reasonable accuracies and significantly cost-effective computational efficiency, the uncertainties in these models still exist at some degrees. There are 15 parameters in CREST 3.0 [28], and each of them has its own physical meaning. These model parameters were calibrated by using the SCE-UA (Shuffled Complex Evolution developed in the University of Arizona) optimization method [46], the uncertainty and sensitivity issues of these parameters do have non-ignorable influences in the performance of distributed hydrological model.
Besides the model structure and parameter uncertainty, data uncertainty issue is also a problem. Satellite-based products treated as data sources reflecting hydrological state variables can be adopted to calibrate the CEREST 3.0, such as soil moisture (Soil Moisture Active Passive, SMAP; and Soil Moisture and Ocean Salinity, SMOS), potential and actual evapotranspiration (MODerate-resolution Imaging Spectroradiometer, MODIS; and Advanced Very High Resolution Radiometer, AVHRR), land surface temperature (MODIS), and total water storage (Gravity Recovery and Climate Experiment, GRACE). Though these satellite-based products provided useful information, the uncertainties of these data still exist, which might be one of the error sources influencing the hydrological model performance, when using the downscaled precipitation dataset with finer spatial resolution (~1 km). The uncertainty issue may be amplified especially when the model was forced by satellite-based products, such as precipitation. The above-mentioned scientific issues should be noticed by modelers and users, and is worthy to be furtherly investigated in future studies.

Future Direction
In this study, only DEM was investigated as the auxiliary variable to downscale the IMERG data. As we know various related environmental variables are related to both spatial patterns and volumes of the happening/occurring rainfall, for example, the cloud information, and brightness temperature [25]. Therefore, these rainfall related variables also should be investigated to downscale the original IMERGE data.

Conclusions
Rainfall information is a perquisite for and plays a vital role in driving hydrological models. However, limited by the observation methods, the obtained precipitation data, at present, are still too coarse. In this study, a new method was introduced to downscale IMERG data at hourly scale. The downscaled precipitation product was compared to the original IMERG data in the temporal-spatial pattern and basin-scale hydrological simulation to assess its optimality. The main conclusions are summarized as follows: (a) the new proposed downscaling method at hourly temporal scale could obtain downscaled results with similar spatial patterns with those of original IMERG data, but with finer spatial resolution; and (b) the downscaled results based on IMERG overall outperformed the IMERG data in hydrological simulations.
Overall, the results showed that spatiotemporal characteristics by IMERG were well captured by the downscaled precipitation across the Ganjiang River basin. The work shed light on the newly proposed high-temporal-scale downscaling method. However, more assessments are still needed to provide a more comprehensive analysis on the downscaled products. In the future, more algorithms and data are still needed to be proposed to obtain downscaled results with finer tempo-spatial resolutions, for example, by introducing cloud information as one of the auxiliary data.
Author Contributions: Z.M., Y.C. and Y.H. designed this Manuscript; X.C., X.J., H.L. and J.L. processed and provided the data; G.K. run the hydrological model derived by the satellite-based precipitation data; Z.M. and X.T. wrote this manuscript; and Y.Y. modified and polished this manuscript.