#### Abstract

Scarce historical flood data in ungauged basins make it difficult to establish empirical and conceptual model forecast in these areas. The physical-based distributed model TOPKAPI is introduced for flood prediction in an ungauged basin by parameter transplant. Five main parameters are selected, and the sensitivity is analyzed by the GLUE method. The Xixian basin and Huangchuan basin in the upper Huaihe basin in China are chosen as study areas. The Xixian basin is regarded as a gauged basin for parameter calibration, and the Huangchuan basin is regarded as an ungauged basin by ignoring the historical discharge data. The model is calibrated in gauged Xixian basin, and then parameters are directly transplanted to adjacent “ungauged” Huangchuan basin to simulate flood forecast in an ungauged basin. The sensitivity analysis shows that soil thickness and soil saturated water content are the most sensitive parameters, and the Manning coefficient of main channel with high Strahler also significantly affects forecast results. According to the simulation results, the TOPKAPI model exhibits good performance in building and the prediction of the ungauged basin, in which the qualified rate of volume and peaks reaches 69.23%, and the average NSE criterion is over 0.67, which is acceptable forecast accuracy and has positive implication for the hydrological forecasting research.

#### 1. Introduction

The hydrological model is a vital tool for the current flood forecasting work by simulating the hydrological process. To estimate the outflow process at the exit section, it generates a certain flood runoff process in the basin under the circumstance of a certain structure and parameters [1]. The hydrological model basing on physical theory, according to the description method and generalization degree of the hydrological process, can be split into conceptual hydrological model and physical-based model [2]. The Xin’anjiang model proposed by Zhao of Hohai University is one of the conceptual hydrological models [3]. It is characterized by less model parameters and high forecasting accuracy, and it has been widely used in flood forecasting in various river basins in China [4]. However, as hydrogeological conditions of basins are highly generalized [5], the parameters should be calibrated using historical flood data [6], making it difficult to build a forecast model in the ungauged basins or basins using partial missing observation data [7]. Though the world’s construction of hydrological stations continues to develop, there are still many small basins without monitoring stations; these small basins are often the focus areas of flood disaster research. For this reason, the flood forecasting of ungauged basins has become the focus and hot topic of existing research, as well as a challenge for all hydrologists [8].

The study on ungauged basin flood forecasting aroused the attention from hydrologists in the end of last century. Vandewiele et al. [9] used the method of regionalization of model parameters to forecast monthly discharge in ungauged region in 1991. Despite using the conventional regression equations to get model parameters, Burn and Boorman [10] employed the basin classification method to find similar basins. Since usually there are insufficient historical discharge data in ungauged basins, the commonest approach to build the hydrological model in ungauged is parameter transferring [11–13], suggesting that the parameters required for models were calibrated in nearby basins with sufficient discharge data and then transplanted to ungauged basins. This approach, however, is based on the similarity of these two basins [14]. Günter [15] summarized several methods for ungauged basin prediction using both event models and explicit soil moisture accounting (ESMA) models. To solve the description of characteristics of ungauged basins, Murugesu [16] proposed the approach of physical controls (basin form and function) of the instantaneous unit hydrograph (IUH), which has been latterly achieved by Ellouze-Gargouri [17], using geomorphological instantaneous unit hydrograph (GIUH) and Copulas. In 2012, Moore et al. [18] made a prediction in ungauged basins in British Columbia (BC) by transposing model parameter from adjacent basin; they also achieved good results. Patil and Stieglitz [19], using a model developed themselves, analyzed selectively transferring sensitive versus insensitive parameters on flood forecasting in ungauged basins. Athira et al. [20] used SWAT model to explore the appropriate functional relationship between the parameters and basin characteristics. Waseem et al. [21] compared different interpolation schemes on parameter spatial distribution in ungauged basins. Hyun and Choi [22] analyzed the relationship between flooding index and rainfall pattern to forecast flood severity in ungauged basins. Sahoo et al. [23] adopted a novel routing method VPMD at ungauged rivers to plot rating curves. Cánovas et al. [24], using the 2D hydraulic model in Spain, conducted the research on flood discharge in ungauged basins.

With the advancement of geographic information technology, remote sensing technology, and computer science, even in the absence of hydrological observatories, the geography of the watershed can be obtained from remote sensing images, including digital elevation model (DEM) [5] and land use and soil classification maps [25]. To fully exploit the advantage of such development, the physical-based distributed hydrological models were created [2, 26]. They describe the underlying surface conditions of basins using physical parameters from remote senses, making it possible to simulate ungauged basin and become a new topic in the field of hydrological forecasting [27]. Among them, the TOPKAPI rainfall-runoff model refers to a distributed hydrological model with grid-based computational unit based on physical basis [28].

Professor Ezio Todini from the University of Bologna in Italy proposed the TOPKAPI (Topographic Kinematic Approximation and Integration) hydrological model in 1995 [29]. This model refers to a physics-based distributed hydrological model based on the study of rainfall-runoff relationship to explore the potential of hydrological model prediction based on physical theory in mountain flood forecast [29]. The model consists of several modules (e.g., evapotranspiration, snowmelt, soil flow, surface runoff, river runoff, and groundwater) [30]. In this model, the whole basin is split into cell grids, and each grid is considered a single calculation unit. Each calculation unit reflects the whole physical hydrological process [31]. This model, compared with conventional conceptual hydrological models (e.g., Xin’anjiang Model [3] and Sacramento Model [32]), fully considers the spatial heterogeneity of the underlying surface conditions, exhibiting the advantage of relatively simple structure, clear parameter meaning, large spatial scale elasticity, as well as wide application fields [33]. Flood simulation and forecasting by this model have been widely performed in many river basins in Italy, Germany, Spain, and other countries; good results have been achieved as well [30].

In this study, the capability to simulate flood in ungauged basin of TOPKAPI model is tested to give suggestions on flood forecast in ungauged basins. The sensitivity of the main parameters of the TOPKAPI model was analyzed, laying the theoretical basis for the parameter transplantation to build forecast model in ungauged basins. The Xixian basin located in the upstream of Huaihe River in China and the adjacent Huangchuan basin were taken as the research basins. The Huangchuan basin was considered the ungauged basin by ignoring the historical flood data. Parameters were calibrated in the Xixian basin and transplanted to the adjacent Huangchuan basin for flood simulation to test the application of TOPKAPI model in ungauged basins without parameter calibration.

#### 2. The TOPKAPI Model

##### 2.1. Model Introduction

The TOPKAPI model refers to a distributed hydrological model based on physical characteristics. It splits the study basin into detailed grids, and thus the difference of physical characteristic in spatial distribution can be fully considered under the calculation unit subdividing [34]. It can also give the specific condition of soil saturation, runoff yield, and confluence depth at each calculation unit point in the basin, which is of high reference implication in real-time flood forecasting, land use and environmental impact assessment, and hydrological process simulation in ungauged basins [35].

According to this model, the hydrological process in each calculation unit is generalized into three nonlinear reservoir equations, which are similar in structure, representing drainage in soil, surface runoff, and channel runoff on saturated soil and impermeable surface, respectively. The finite difference method is used in the calculation, so the correlated four surrounding grids are considered during the calculation of each grid unit [36]. The parameters used to describe the underlying surface conditions (e.g., soil thickness, vertical and horizontal saturated hydraulic conductivity, coefficient of nonlinear reservoir equation, and Manning coefficient of surface and evaporation coefficient of vegetation) could be extracted directly from the soil type classification map and land surface utilization map. Besides, the elevation and gradient distribution of the basin could be reflected by DEM (digital elevation model).

In terms of specific mechanism, TOPKAPI model employs Thornthwaite evaporation formula to calculate the potential evaporation of different vegetation covers in different growth stages; it also calculates the actual evaporation according to the actual wettability of upper soil [28]. According to the change of soil water content in the upper unsaturated region, the runoff in the soil and the surface runoff are first calculated and then aggregated to the total runoff. In the confluence part, the nonlinear reservoir equation is adopted to merge surface runoff and underground runoff for confluence calculation [28, 37]. The main parts and approaches of TOPKAPI model are listed in Figure 1.

To describe the movement of water in surface, soil, and river channel, TOPKAPI model adopts three nonlinear reservoir equations with a similar structure. In the earliest version, the variable step five-stage Runge–Kutta numerical method is employed to solve the differential equations. Later, it was found that under appropriate approximation conditions, these equations could be solved using analytical methods [28].

##### 2.2. Comparison with Conceptual Model

From the aspect of data required for the model building, the conceptual hydrological model should be based on the previous historical hydrological data. Taking Xin’anjiang model as an example, the parameters, including water storage capacity (WM), evapotranspiration coefficient (*K*), outflow and regression coefficient of the basin (KI and KG), should be calibrated through historical floods [3]. In other words, conceptual model is based on high-degree generalization of those physical parameters to reflect the characteristics of basins. This model has the advantage of relatively flexible adjustable parameters that can be fully adjusted to give optimized forecast results and avoid some relatively complex mechanisms in the hydrological process [38]. However, the disadvantage lies in the reliance on historical flood data, which is particularly evident in some areas in the absence of historical rainfall and runoff data where it is almost likely to calibrate parameters [39]. The existing feasible method is to search for data basins with similar geological and hydrological conditions for parameter transplantation, which will inevitably face systematic error of the model. In contrast, most of the data required for describing hydrological conditions and underlying surface conditions of watershed based on physical basis can be obtained through actual measurement (e.g., DEM, soil classification and soil thickness, hydraulic conductivity, and evaporation coefficient of land classification) [34]. These parameters are relatively easy to obtain online, including the Shuttle Radar Topography Mission (SRTM) launched by the Consultative Organization for International Agricultural Research (CGIAR), which can obtain the global 90 m resolution digital elevation model free of charge from its official website; the accessible website of the University of Maryland (UMD) in the United States. On this point, the hydrological model based on physical basis has obvious advantages over the conventional conceptual hydrological model when it is relatively easy to obtain the parameters required for model building in ungauged areas [28, 40].

#### 3. Study Area and Dataset

##### 3.1. Study Area

The neighboring Xixian and Huangchuan river basins (Figure 2) above the Wangjiaba basin in the upper reaches of the Huaihe River were taken as research area. The Wangjiaba basin is located in the southern part of Henan Province of China. It is a monsoon climate with an average annual rainfall of 1060 mm [30]. About 50% of the rainfall concentrates in the flood season, and the Xixian and Huangchuan river basins are the origins of the upstream river. The control area of Xixian Hydrological Station is 10 190 km^{2}. There are two large reservoirs in Nanwan (water catchment area 1058 km^{2}) and Shishankou (water catchment area 306 km^{2}) [41]. After deducting two large reservoirs, the catchment area is 8400 km^{2} [42]. The elevation range of the Xixian County is 37–963 m [43]. The terrain in the northeast is relatively flat, being mountainous in the west, southwest, and south. The soil type is mostly clay soil. It has soil types (e.g., loam, sandy loam, and clay). Most of the land is farmland, with a small portion of forest land and mixed forest land. The Huangchuan basin outlet station has a controlled area of 1989 km^{2} with Pohe reservoir built inside [44]. The elevation range of the area is 36–984 m. The watershed is mostly mountainous, the north is relatively flat, most of the soil types are clay and sandy loam, and the surface vegetation is primarily farmland [45, 46].

**(a)**

**(b)**

The digital elevation data of the basins could be obtained from the SRTM global 90 m digital elevation model provided by the Consultative Group on International Agricultural Research (Figure 3). GIS software was employed to extract the simulated stream nets and basin boundaries in the Xixian and Huangchuan basins. The rainfall data of 13 rainfall stations and the discharge data of Xixian hydrological station and that of the two reservoirs were collected from 1980 to 2003 according to the distribution of stations in the basin. The discharge data, as well as the rainfall data of 5 rainfall stations in the Huangchuan basin and the discharge data of the Huangchuan hydrological station from 2003 to 2005, were processed into 1 h interval. The meteorological data of Xinyang, Zhumadian, and Guangshan meteorological stations (including daily maximum and minimum temperatures, wind speed, humidity, and sunshine hours) were collected from 1980 to 2005.

##### 3.2. Underlying Surface Analysis

The underlying surface parameters of the TOPKAPI model were transplanted on the premise that the soil types and land cover of underlying surface of the two basins should be of the same type or have similar physical properties, so the distribution of soil (Figure 4 and Table 1) and surface in the two basins (Figure 5 and Table 2) significantly affected the ungauged area prediction [47]. According to the collected data, the Xixian and Huangchuan basins were located close to each other. Their soil types were similar to clay, and most of the land surface was the same type of farmland. The soil parameters and land use parameters calibrated by the Xixian basin in the TOPKAPI model were applicable to the adjacent Huangchuan basin. The soil data are derived from the FAO soil map of the world, Global soil profile databases, which is available online at http://www.ngdc.noaa.gov/seg/eco/cdroms/reynolds.htm. The land cover map is derived from National Administration of Surveying Mapping and Geoinformation of China (NASG), which is free downloadable for all Internet at http://www.globallandcover.com/GLC30Download/index.aspx and the characteristic data for each type of land cover are derived from Corinne Land Cover 2006 raster data by European Environment Agency, which is accessible at https://www.eea.europa.eu/data-and-maps/data/clc-2006-raster.

**(a)**

**(b)**

**(a)**

**(b)**

#### 4. Parameter Sensitivity Analysis

##### 4.1. Parameter Selection

In TOPKAPI model, the nonlinear reservoir equation can be written aswhere the variable *y* denotes average soil water content, river storage, surface water depth, etc. Variables *a*, *b*, and *c* keep constant in a single time step.

In the calculation of groundwater, the variables *a*, *b*, and *c* are expressed as follows:where denotes the precipitation, *X* is the size of grid unit, and refer to surface and subsurface inflow of a grid unit, respectively, and is the local hydraulic conductivity.

In the calculation of surface water, the variables a, *b*, and *c* are expressed as follows:where denotes the width of grid without river channel, represents the net precipitation, and is index in Manning equation with a constant value of 5/3.

In the calculation of channel water propagation, the variables *a*, *b*, and *c* are expressed as follows:where denotes the lateral inflow, including soil outflow to surface and channel.

The calculation formulas show that the main influence in the calculation of soil water originates from the upstream inflow, grid size, and hydraulic conductivity, in which hydraulic conductivity as the unknown coefficient of the nonlinear equation will directly affect the solution of the unknown term. Thus, hydraulic conductivity is considered the key parameter to be investigated. The calculation of surface water on slope is affected by slope, precipitation and surface Manning coefficient, and surface runoff after deducting infiltration is associated with soil conditions on underlying surface. Accordingly, in the part of surface flow calculation, the Manning coefficient of surface, the saturated water content associated with soil characteristics, and the soil thickness describing soil volume were taken as the key research parameters. In the calculation of channel flow using nonlinear reservoir method, the Manning coefficient of river roughness was taken as the main calibration parameter. All the selected parameters are listed in Table 3.

##### 4.2. Sensitivity Analysis

The GLUE method was first proposed by Beven and Binley in 1992 based on the RSA (regionalized sensitivity analysis) method proposed by Hornberger and Spear in 1981 [48]. Beven and Binley proposed that this method is not the numerical value of a single parameter, but the value of a group of parameters that significantly affect the prediction results of the model. In this method, Monte Carlo method is adopted to sample a series of parameter combinations in prior parameter range and substitute them into the model calculation. Likelihood function is selected as the evaluation and instruction of the parameter sets. The parameter set whose likelihood function value is higher than a certain threshold is considered to be able to represent the characteristic of study basin (which is called “behavioral”), whereas the parameter set with likelihood function value lower than the threshold is considered fail at reflecting it and be abandoned [49].

Several criterions could be used to describe the prediction performance of a model, among which the Nash–Sutcliffe model efficiency coefficient (NSE) [50] has been most widely used in model prediction performance judgment. In this study, the NSE is chosen as the likelihood function. The Nash–Sutcliffe model efficiency coefficient is expressed aswhere and denote the modeled and observed discharge at time *t*, respectively, and is the average value of observed discharge over the whole time series.

In this study, 3000 sets of parameters including those selected in the previous part were generated by random Monte Carlo sampling methodology within the prior range in Table 4. A threshold of 0.7 NSE is chosen to check whether these parameter sets are behavioral. Parameter uncertainty was analyzed by 8 floods in Xixian basin from 1980 to 2002 in TOPKAPI model. A total of 283 sets of parameters were found to be behavioral achieving 0.7 NSE likelihood.

Posterior parameter distribution was derived from 283 behavioral parameter sets. To get visual impression of these parameter sets, the histograms of all 18 parameters distribution are plotted in Figure 6. Histograms with sharp peaks or steep slopes means parameter are well identified, whereas parameters with relative flat histograms are associated with more uncertainty [51]. Because all the behavioral sets have performance over threshold in common, narrow and sharp peaks reflect relative strong relationship with model results and show that this parameter has more sensitivity, which means variation in this parameter may significantly affect model results. A statistic of parameter frequency was conducted in Table 5 by calculating standard deviation of frequencies of each parameter to get a numerical comparison.

**(a)**

**(b)**

**(c)**

In Table 5, the standard deviations of parameters frequency distribution are listed. It can be clearly seen that horizontal hydraulic conductivity has the least values and the histograms are all relatively flat, meaning that this parameter has low sensitivity and high uncertainty. This sort of parameters are difficult to be calibrated due to their weak relationship with model outputs, but they can only slightly affect the model performance. However, emphasis in calibration should be put on those parameters with high standard deviation and sensitivity such as soil thickness, saturated water content, and riverbed Manning coefficient. Soil type with high coverage percentage such as clay (with 71.48% coverage) has significantly high standard deviation value for both thickness and saturated water content. This is the same for the land use Manning coefficient of farmland (with 67.70% coverage) and the riverbed Manning coefficient of V level river channel.

#### 5. Application in Ungauged Basins

To verify the effectiveness of the model in ungauged basins, experiments were performed. The Xixian basin and adjacent Huangchuan basin in Wangjiaba were taken as study areas. The parameters of the model were calibrated in Xixian basin and subsequently transferred to Huangchuan basin, and flood simulations were carried out to verify the application effect of the model in ungauged basins.

The calibration and validation results are evaluated by “qualified rates” [52]. A flood forecast with flood volume and flood peak error less than 20% is considered to be “qualified,” and the “qualified rate” means the percentage of qualified flood simulations over the total number. The whole forecast scenario accuracy will be evaluated to a certain level according to Table 6.

##### 5.1. Parameter Calibration

Most parameters in physical-based hydrological model could be obtained by actual measurement; however, the measured values acquired on the point measurement might not be sufficiently representative at the calculation unit scale and might not reflect the temporal and spatial scale change within the unit. Thus, the parameters obtained by measurement or reference to the literature should be fine-tuned by trial before the simulation to compensate for the error caused by the spatial generalization of the parameters and the time generalization of the input parameters, making the results of the simulation closer to the actual situation.

The calibration parameters included soil horizontal water conductivity, saturated water content, land surface, and riverbed Manning coefficient for Muskingum algorithm. In the previous GLUE assessment, 3000 sets of parameters have been generated. Several “best” performed sets of parameters with reasonable values and high NSE are chosen for manually adjustment. These parameters were calibrated by the rainfall and discharge data of 16 floods from 1991 to 2003 in Xixian basin. Run the TOPKAPI model with the meteorological data (including precipitation and temperature) of 13 floods in Huangchuan watershed from 2003 to 2005 together with the calibrated parameters and adjusted these parameter according to the outcome flood results. Furthermore, to verify the accuracy of the simulation and to demonstrate the application of the model in ungauged basins, the flood simulation results were compared with the measured results.

The flood peak and NSE were taken as objective functions to calibrate the parameters. The percentage of initial soil water content in the basin would significantly affect the simulation results of the early floods in the basin. In this calibration, to determine the initial soil moisture, the “warm-up” method [53] was adopted to simulate the initial state through natural situation, which was initiated 30 days before the first flood in advance.

Figure 7 shows that, in the calibration, most of the stream flow series are well reproduced. Floods with large water volume are simulated better than those with low volume. From the results of model calibration in Table 7, the simulation results of the model in Xixian basin were good, and the average NSE were above 0.8. Among the 16 floods used in the calibration, 14 floods had NSE over 0.8, taking up 87.5% and 6 floods over 0.9. In terms of Nash–Sutcliffe model efficiency coefficient, the model exhibited certain applicability in this basin. According to the simulated flood volume, the qualified rate of TOPKAPI model flood volume was 81.3% in the 16 floods used for calibration and 93.8% in the flood peak simulation. From the results, it is revealed that under the underlying elevation data, soil and land use distribution parameters and basin hydrometeorological data meet the demand, the model achieved relatively good simulation in the study area. The average NSE were all above 0.7, achieving the accuracy of Class B forecasting [52], and the average error between flood volume and flood peak was controlled within 20%. This is associated with the wet area where the research basin is located. The full storage and runoff yield models used in the model were applicable to the wet area of Huaihe River, and the better the hydrogeological conditions of the basin, the smaller the error of the model confluence calculation will be. Thus, TOPKAPI model exhibits good applicability in this basin. All the calibrated parameters are listed in Tables 8–10.

From the process of model building, TOPKAPI model did not adjust too much parameters in the calibration process, and the prediction results were good. Accordingly, TOPKAPI model based on physical basis exhibited better applicability when building model in adjacent areas without historical flood data. However, the prediction accuracy of TOPKAPI model depended on the high-precision partition of computational grids, whereas high-resolution grids would increase the requirement of computational ability. Thus, to obtain higher prediction accuracy, the model requires better operation equipment and long calculation time.

##### 5.2. Parameter Validation

The parameters obtained by the calibration in adjacent basin were directly applied in Huangchuan basin. Besides, the flood discharge process of the outlet section was calculated using the rainfall data of 13 floods in 2003–2005 as input to simulate the application of the model in the ungauged basin. From the statistical Table 11, Figure 8, and the box plot Figure 9 of calculation and results, it is revealed that the overall NSE was relatively lower because of parameter transplantation. The average NSE was 0.67, and the qualified rate of flood volume and peak simulation was 69.23%. This was probably because Xixian is in the warm temperate zone and belongs to the semihumid basin, and the overall basin is dominated by excess of storage mode, which is consistent with the basic assumption of TOPKAPI model runoff generation. The accuracy forecast in ungauged basins, compared with the calibration results, decreased in varying degrees. This was because the parameter transplantation inevitably brings about inconsistencies with the actual conditions, which is also the most important problem facing the prediction of ungauged basins.

**(a)**

**(b)**

**(c)**

It was found that the flood process of TOPKAPI model had steep rise and fall just, as shown in the hydrographs. Accordingly, following the flood peak priority calibration strategy, the systematic error of overall smaller flood volume appeared. The water confluence process was faster than the actual situation, which is also one of the parts of the model to be further optimized. Nevertheless, TOPKAPI model still reflects good transplantability in flood peak prediction in ungauged basins. Since remote sensing data could reflect actual topography, vegetation, and soil types as much as possible, the eligibility rate of flood peak prediction reached about 70%. This is the advantage brought by detailed hydrological process based on the physical basis of the distributed model, fully displaying the physical conditions of underlying surface.

#### 6. Discussion and Conclusion

In this study, a physically based distributed hydrological model TOPKAPI was applied in upper Huaihe, Xixian, and adjacent Huangchuan basin. The model uncertainty caused by spatial and temporal generalization of model parameters was evaluated using GLUE method. Through uncertainty analysis and comparison of calibration simulation results with measured data and prediction results of the model in adjacent similar ungauged basins, parameter transplantation was tested. Based on the results, the following conclusions were drawn:(1)The description of basin characteristics in TOPKAPI model can be characterized by underlying surface data from land use and soil type obtained using remote sensing technology, which is easy to obtain and use. In this study, it was assumed that Huangchuan basin is an ungauged basin without any hydrological stations, whereas underlying surface data can still be acquired from satellite or remote sensing equipment, making it likely to describe the basin in hydrological model.(2)According to the results of the sensitivity analysis of the model parameters, the soil thickness and surface Manning coefficient, with a large proportion of coverage, significantly affected the prediction results of the model. The horizontal hydraulic conductivity has the least sensitivity. Since the classification of river depends on the area of the whole river basin and other factors, some differences would exist among different river basins for the same order. In this study, this is also a reason for the reduction of prediction accuracy of ungauged river basins. The rules of parameter transplantation in prediction research should be elucidated further.(3)The model exhibits good applicability in the upper Huaihe basin. The base assumptions of TOPKAPI model made it feasible in humid semihumid area, and thus the flood volume flood peak as long as NSE all achieved high accuracy in calibration.(4)In the prediction of ungauged basins, it was found that under premise of sufficient remote sensing data of underlying surface, the transplantation of parameters among similar basins could achieve good forecasting results.

Besides, the high prediction accuracy of TOPKAPI model should be based on detailed meshing of computing units. For instance, the grid resolution adopted in this study was 200 m, and the prediction results were relatively good. However, the calculation period took a long time. If the model is applied in a larger watershed, the problem will be more prominent. Thus, how to make an optimal trade-off between the calculation time of the model and the prediction accuracy needs further studies.

#### Data Availability

The digital elevation model data can be acquired from CGIAR (http://srtm.csi.cgiar.org/). The soil data are derived from the FAO soil map of the world, Global soil profile databases, which is available online at http://www.ngdc.noaa.gov/seg/eco/cdroms/reynolds.htm. The land cover map is derived from National Administration of Surveying Mapping and Geoinformation of China (NASG), which is freely downloadable for all Internet at http://www.globallandcover.com/GLC30Download/index.aspx and the characteristic data for each type of land cover are derived from Corine Land Cover 2006 raster data by European Environment Agency, which is accessible at https://www.eea.europa.eu/data-and-maps/data/clc-2006-raster.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by the National Key Research and Development Program of China (grant no. 2018YFC1508103) and the Fundamental Research Funds for the Central Universities (grant no. 2016B04714).

#### Supplementary Materials

The following data are available in supplementary materials: the precipitation data, temperature, and discharge data used for parameter calibration in Xixian basin and the coordinate of meteorological stations and precipitation gauges used in this research.* (Supplementary Materials)*