Application and Research of Liuxihe Model in the Simulation of Inﬂow Flood at Zaoshi Reservoir

: Floods occur frequently in China, and watershed ﬂoods are caused mainly by intensive rainfall, but the spatial distribution of this rainfall is often very uneven. Thus, a watershed hydrological model that enables a consideration of a heterogeneous spatial distribution of rainfall is needed. In this study, a ﬂood forecasting scheme based on the Liuxihe model is established for the Zaoshi Reservoir. The particle swarm optimization (PSO) algorithm is used to optimize the model parameters for ﬂood simulation, and the model’s performance is assessed by a comparison with measured ﬂood data. The spatial distributions of rainfall selected for this study are non-uniform, with much greater rainfall in some areas than in others in some cases. Rainfall may be concentrated in the middle of the basin, in the reservoir area, or in the upstream portion of the basin. The Liuxihe-model-based ﬂood inﬂow forecasting scheme for the Zaoshi Reservoir demonstrates an excellent simulation effect, with an average peak simulation accuracy of 96.3%, an average peak time of 1.042 h early, and an average Nash–Sutcliffe coefﬁcient of 0.799. Under the condition of an uneven spatial distribution of rainfall, the Liuxihe model simulates ﬂoods well. The PSO algorithm signiﬁcantly improves the model’s simulation accuracy, and its practical application requires only the selection of a typical ﬂood for parameter optimization. Thus, the ﬂood simulation effect of the Liuxihe model is ideal for the watershed above the Zaoshi Reservoir, and the scheme developed in this study can be applied for operational ﬂood forecasting.


Introduction
Flooding has been one of the most widespread and destructive types of natural disasters since ancient times [1][2][3][4].With global climate change [5][6][7], the impact of human activities, and accelerating urbanization [8][9][10][11], flood disasters have become increasingly frequent, with non-negligible impacts on human welfare, sustainable economic development, and the natural environment [12].Reducing flood disaster risk has become an important component of disaster prevention and mitigation [13].In traditional flood forecasting, the prediction scheme is constructed mainly by artificial means.A typical rainstorm process in the region is observed and measured, and a formula to convert rainfall into discharge is constructed for artificial forecasting.However, the accuracy of traditional flood prediction models tends to decrease with increasing prediction cycle length [14], making it difficult to meet timeliness and accuracy requirements for practical application.For this reason, these traditional models have been gradually replaced with automatic forecasting based on hydrological models.In China, reservoirs are built mainly in mountainous and hilly areas with complex and diverse landforms, where rainfall intensity is high, and the flood burst is strong.For these reasons, rainfall and flood process data are difficult to obtain.The accurate, quick, and automatic forecasting of floods in such reservoir basins has become a focus of many Chinese and international scholars.
Watershed floods are caused mainly by torrential rain, and the spatial distribution of precipitation [15] has an important impact on flood warning and forecasting.Thus, the inflow flood forecasting model used for prediction must reflect the impact of the uneven spatial distribution of rainfall on floods.
At present, watershed flood forecasting models include lumped [16][17][18] and distributed [19] models.The lumped model considers the watershed as a whole and does not distinguish between slope confluence and channel confluence.The structure and parameters of the lumped hydrological model reflect the average situation across the whole basin, and input elements are spatially uniform, making it difficult to consider the heterogeneity of the underlying surface or other types of spatial heterogeneity.Additionally, a large amount of historical flood data is needed to optimize the model parameters, making this model unsuitable for mountainous watersheds with limited data.The distributed model is a hydrological model based on physical processes.It divides the whole basin into refined units based on their physical characteristics and calculates production and confluence at the unit level, enabling the description of the movement of a flood along the basin's slope and river channel and, thus, a consideration of the uneven spatial distributions of rainfall and the impact of reservoir storage on flood flows.The distributed hydrological model was first proposed by Freeze and Harlan [20], who noted areas that require further study.After years of development, a variety of distributed hydrological models have been proposed in China and other countries.They include the SHE model [21], VIC model [22], WetSpa model [23], Vflo model [24], LL model [25], Liuxihe model [26], etc.
In this paper, the Liuxihe model is used to formulate a flood forecasting scheme.The influence of the spatial distribution characteristics of rainfall events in the study area on the model simulation results is analyzed.The particle swarm optimization (PSO) algorithm [27][28][29][30] is used to optimize the model parameters before flood process simulation.Based on a comparison of the results, the flood simulation accuracy of the model under the condition of uneven rainfall in the basin is evaluated quantitatively to demonstrate its feasibility for flood forecasting for the Zaoshi Reservoir, providing data to support decision making related to flood forecasting and warning in this basin and the implementation of flood control measures in downstream areas.

Liuxihe Model
The Liuxihe model is a distributed physically based hydrological model developed by Yangbo Chen [26,31], which is used mainly for flood forecasting and simulation at the basin scale.The model has six sub-models for watershed division, evapotranspiration calculation, runoff calculation, confluence calculation, parameter determination, and simulation calculation.It has been used to achieve gridded, fully distributed watershed flood simulation and forecasting, and it meets the accuracy requirements of engineering applications.The Liuxihe model has been applied successfully to reservoir inflow flood forecasting, small and medium rivers, mountain flood disaster warning and forecasting, coupled hydro-meteorological flood forecasting for large watersheds, flood forecasting for urban watersheds, and hydrological process simulation [32,33].

Study Area
The Zaoshi Reservoir is located in the downstream portion of the Xieshui River, a first-order tributary of the Lishui River, and its dam is 19 km from Shimen County and 2 km from Zaoshi Town.The main function of the reservoir is flood control.It also supports functions such as power generation, irrigation, and shipping, and acts as a hub of water conservancy for general utilization.Figure 1 is a map of the Zaoshi Reservoir watershed.The Xieshui Basin, shaped like a crescent moon, is located at 29 • 36 -30 • 11 N and 110 • 27 -111 • 22 E. The river originates from Du Mountain, Wufeng County, Hubei Province, and flows mainly through Shimen County, Hunan Province.The Nishi, Suoshi, Moshi, and Zaoshi rivers merge into the main stream of Lishui at Sanjiangkou, located 2 km upstream of Shimen County.The overall area of the Xieshui Basin is 3201 km 2 , the river length is 171 km, and the total drop is 1401 m.The rainwater collection area above the dam is 3000 km 2 , accounting for 93.7% (3201 km 2 ) of the area of the Xieshui watershed.The river flows through a mountainous area.The terrain of the basin is high in the northwest and low in the southeast.The elevation is 2098.7 m at Huping Mountain in the northwest and decreases gradually toward the southeast, reaching 150-250 m in the downstream portion of the river.This river has a steep slope and rapid flow and, thus, is classified as a mountain river.The basin is located in a zone with a humid mid-subtropical monsoon climate, in an area of heavy rainfall in the Yangtze River basin.Rainfall is concentrated from May to August, with precipitation of about 700-1200 mm, accounting for 55-60% of the total annual precipitation.This region has characteristics of high precipitation variability, heavy rainstorms, continuous droughts in summer and autumn, huge differences in runoff characteristics, frequent floods [34], and an uneven spatial distribution of rainfall.
15, x FOR PEER REVIEW 3 of 14 Zaoshi rivers merge into the main stream of Lishui at Sanjiangkou, located 2 km upstream of Shimen County.The overall area of the Xieshui Basin is 3201 km 2 , the river length is 171 km, and the total drop is 1401 m.The rainwater collection area above the dam is 3000 km 2 , accounting for 93.7% (3201 km 2 ) of the area of the Xieshui watershed.The river flows through a mountainous area.The terrain of the basin is high in the northwest and low in the southeast.The elevation is 2098.7 m at Huping Mountain in the northwest and decreases gradually toward the southeast, reaching 150-250 m in the downstream portion of the river.This river has a steep slope and rapid flow and, thus, is classified as a mountain river.The basin is located in a zone with a humid mid-subtropical monsoon climate, in an area of heavy rainfall in the Yangtze River basin.Rainfall is concentrated from May to August, with precipitation of about 700-1200 mm, accounting for 55-60% of the total annual precipitation.This region has characteristics of high precipitation variability, heavy rainstorms, continuous droughts in summer and autumn, huge differences in runoff characteristics, frequent floods [34], and an uneven spatial distribution of rainfall.

Collation of Measured Flood Data
For this study, data on a total of 25 major flood processes occurring since 2009 were collected.The Zaoshi Reservoir area has four hydrological stations in Zaoshi, Suoshi, Weixinchang, and Jiangping and nine rainfall stations in Kangjiaping, Nishi, Nanzhen, Luoping, Nanping, Zhuxiangkou, Guanzhuang, Jinping, and Hekou (Figure 1).In this study, hourly rainfall data from the rain stations and the Thiessen polygon method [35][36][37] were used for spatial interpolation and watershed division according to the distribution of the rain gauges.The rainfall intensity was assigned to each grid unit based on the rainfall data collected at the rain stations within it.Detailed information about the floods is provided in Table 1.

Collation of Measured Flood Data
For this study, data on a total of 25 major flood processes occurring since 2009 were collected.The Zaoshi Reservoir area has four hydrological stations in Zaoshi, Suoshi, Weixinchang, and Jiangping and nine rainfall stations in Kangjiaping, Nishi, Nanzhen, Luoping, Nanping, Zhuxiangkou, Guanzhuang, Jinping, and Hekou (Figure 1).In this study, hourly rainfall data from the rain stations and the Thiessen polygon method [35][36][37] were used for spatial interpolation and watershed division according to the distribution of the rain gauges.The rainfall intensity was assigned to each grid unit based on the rainfall data collected at the rain stations within it.Detailed information about the floods is provided in Table 1.To facilitate analysis, the data for the three underlying surfaces used in this study were resampled to a spatial resolution of 90 m × 90 m.Five land cover types are present in the study area: evergreen coniferous forest, evergreen broad-leaved forest, shrub, thin forest, and cultivated land, which account for 17.15%, 12.39%, 69.20%, 0.14%, and 1.12%, respectively, of the basin area.Shrub occupies the largest proportion (>40%) of the basin area.In total, 15 soil types are present, among Five land cover types are present in the study area: evergreen coniferous forest, evergreen broad-leaved forest, shrub, thin forest, and cultivated land, which account for 17.15%, 12.39%, 69.20%, 0.14%, and 1.12%, respectively, of the basin area.Shrub occupies the largest proportion (>40%) of the basin area.In total, 15 soil types are present, among which low-fertility and strongly acidic soil (Haplic Acrisol) occupies the largest proportion (>20%) of the basin area.The spatial distributions of land use type and soil type data show little variation and are represented fully by the data used for this study.

Construction of Liuxihe Model 2.4.1. Watershed Division and Cell Type Determination
The Liuxihe model divides the watershed into grid cells with independent physical characteristics and rainfall, classified as channel cells, slope cells, and reservoir cells.The model determines the flow direction of each unit based on the D8 [38] flow direction method, and the generated flow diagram is used to calculate the cumulative flow for each cell.According to the cumulative flow threshold, channel cells are extracted, after which the normal water level of Zaoshi Reservoir is used as the threshold to distinguish reservoir cells from slope cells.In this study, the Strahler method [39] was employed to classify the extracted river channels, and flow accumulation (FA 0 ) values were used to divide the river channels into a maximum of six orders.The corresponding critical FA 0 values and number of river channel cells assigned to each order are listed in Table 2.As not all corresponding tributaries for first-and second-order river channels are included in the model, these levels should not be selected.For fourth-, fifth-, and sixthorder river channels, the river divisions are too closely spaced, and preliminary analysis of the river section using Google Earth remote sensing images showed that the estimation of channel section size is difficult.Additionally, the branches of fourth-, fifth-, and sixth-order river channels are very dense, and their deviation from the real river is large, making it inconvenient to measure river width from remote sensing images.Thus, third-order river channels were selected for analysis in this study.A range of FA 0 values for the generation of third-order river channels is assigned, within which the threshold can be adjusted.The SHP (shapefile) files generated for a third-order river system with the adjusted threshold are imported into the remote sensing image to determine the deviation from the real river system.In this study, 5000 was selected as a reasonable threshold.The river classification results are shown in Figure 3.To facilitate the estimation of river section size in the Liuxihe model according to the unit division and river classification results, as well as changes in the Google Earth remote sensing images and DEM of the basin above Zaoshi Reservoir, nodes are set, and virtual river reaches are divided into locations of large tributary confluences, river width changes, and river flow direction changes.Thus, the river bottom width, side slope, river roughness, and bottom slope properties are set.

Initial Model Parameters
Among the parameters characterizing the land use type, the initial value for the slope element roughness coefficient is determined according to the method of the Liuxihe model, and the evaporation coefficient is set uniformly to 0.7 [31] (Table 3).Soil type parameters include the thickness of soil layer, saturated water content, field moisture retention, wilting moisture content, and saturated hydraulic conductivity.Empirical thickness values for various soil types are shown in Table 4.The soil porosity characteristic, b, is set uniformly to 2.5 [31] for all cells, based on the literature.The saturated water content, field moisture retention, saturated hydraulic conductivity, and wilting moisture content are calculated using the soil hydraulic properties calculator proposed by Arya et al., which was developed by Saxton et al. and can be operated online ("http://www.bsyse.wsu.edu/saxton"(accessed on 15 November 2022)) [40].Potential evapotranspiration is determined empirically from the climatic conditions of the basin, with one value (0.23 mm/h) applied across the whole catchment.The underground water recession coefficient is set empirically to 0.995.The topographic parameters of flow direction and slope are non-adjustable model parameters calculated from DEM data (Figure 4).[40].Potential evapotranspiration is determined empirically from the climatic conditions of the basin, with one value (0.23 mm/h) applied across the whole catchment.The underground water recession coefficient is set empirically to 0.995.The topographic parameters of flow direction and slope are non-adjustable model parameters calculated from DEM data (Figure 4).

Model Parameter Optimization
The Liuxihe model consists of 13 parameters per unit, which can be classified into four types: meteorological parameters, topographical parameters, soil parameters, and vegetation parameters.Meteorological parameters include potential evapotranspiration, topographical parameters include flow direction and slope, vegetation parameters include evaporation coefficient and slope roughness, and soil parameters include the thickness of soil layer, soil porosity characteristics (b), saturated water content, field moisture retention, wilting moisture content, and saturated hydraulic conductivity.The flow direction and slope are determined using DEM calculation and remain unchanged; they are also known as non-adjustable parameters.Other adjustable parameters are optimized based on initial values using the particle swarm optimization (PSO) algorithm.
In 1995, when American psychologist James Kennedy and electrical engineer Russell Eberhart were studying the diversity of biological populations, they found that birds showed a particular social behavior while hunting.Inspired by this behavior, they proposed the framework of the PSO algorithm [41].
The idea of the PSO algorithm is as follows: Each particle in the algorithm represents a set of parameter solutions.The particle adjusts its velocity and direction by remembering and following the individual and global best positions, thus achieving the optimization process.The particle's velocity and position transformations are implemented using the following formulas: where V i,k is the velocity of the ith particle at time k, X i,k is the position of the ith particle at time k, X i,pBest is the optimal position for individual particle i at time k, X gBest is the global optimal position of all particles at time k, ω is the inertia acceleration speed, C 1 and C 2 are learning factors, and rand is a random number between 0 and 1.Meanwhile, the performance of the algorithm can be improved by dynamically adjusting the values of ω.For ω, the linearly decreasing inertia weight strategy (LDIW) proposed by Shi [42] is used for dynamic adjustment in this paper.The calculation formula is as follows: where i is the current evolution number, MaxN is the maximum evolution number, and ω max and ω min are the maximum and minimum values of the inertia acceleration speed ω, which are 0.9 and 0.1, respectively.
Sustainability 2023, 15, 9857 8 of 14 The Liuxihe model is a distributed model that employs the PSO algorithm for parameter optimization, and the advantage of this model is that data from only one flood are needed for this optimization.Flood 20090607 was selected for model parameter optimization, and the other 24 floods were used for model validation.For PSO algorithm parameter setting, the population size (particle number) of the particle swarm is set to 20.The number of iterations is set to 200, and the total number of calculations is 1000.The inertia factor value range is [0.1, 0.9], with a linear decrease within this range occurring during optimization.The learning acceleration factors C1 and C2 have a value range of [0.5, 2.5] and are optimized dynamically and iteratively within this range according to the anticosine acceleration algorithm.The parameter optimization process and its results are shown in Figure 5.
posed by Shi [42] is used for dynamic adjustment in this paper.The calculation formula is as follows: where i is the current evolution number,  is the maximum evolution number, and  and  are the maximum and minimum values of the inertia acceleration speed ω, which are 0.9 and 0.1, respectively.
The Liuxihe model is a distributed model that employs the PSO algorithm for parameter optimization, and the advantage of this model is that data from only one flood are needed for this optimization.Flood 20090607 was selected for model parameter optimization, and the other 24 floods were used for model validation.For PSO algorithm parameter setting, the population size (particle number) of the particle swarm is set to 20.The number of iterations is set to 200, and the total number of calculations is 1000.The inertia factor value range is [0.1, 0.9], with a linear decrease within this range occurring during optimization.The learning acceleration factors C1 and C2 have a value range of [0.5, 2.5] and are optimized dynamically and iteratively within this range according to the anticosine acceleration algorithm.The parameter optimization process and its results are shown in Figure 5.The results confirm that the optimization effect of parameters obtained using the PSO algorithm is good and that the simulated flood process is consistent with the measured flood process.The peak error for flood 20090607 was reduced from 0.443 to 0.028 after parameter optimization (Table 5).Other evaluation indices verify the high precision and good effect of optimization.Thus, the optimized flood model parameters can be considered to be close to the true values.The results confirm that the optimization effect of parameters obtained using the PSO algorithm is good and that the simulated flood process is consistent with the measured flood process.The peak error for flood 20090607 was reduced from 0.443 to 0.028 after parameter optimization (Table 5).Other evaluation indices verify the high precision and good effect of optimization.Thus, the optimized flood model parameters can be considered to be close to the true values.

Flood Simulation Results
The optimized model parameters were used to simulate the remaining 24 flood processes.Selected simulation results are shown in Figure 6.

Model Performance Evaluation
In order to quantify the coincidence between measured and simulated flows, six evaluation indices were used: the Nash-Sutcliffe coefficient, correlation coefficient, process relative error, water balance coefficient, peak error, peak time error, and average value.They were calculated using formulas obtained from the literature [43].The results are provided in Table 6.
With the application of the optimized parameters, the average peak simulation accuracy is 96.3%, the average peak time is 1.042 h early, and the average Nash-Sutcliffe coefficient value is 0.799.Additionally, the simulated flood process curve aligns well with the measured flood process curve.These results demonstrate that the Liuxihe model can simulate floods in a mountainous area accurately to support inflow flood forecasting.

Flood Simulation Results
The optimized model parameters were used to simulate the remaining 24 flood processes.Selected simulation results are shown in Figure 6.

Simulation to Flood with Uneven Precipitation Spatial Distribution
To assess the model's simulation performance, spatial distribution maps of accumulated rainfall for the 24 floods were drawn and analyzed.The uneven spatial distribution of rainfall (with much more rainfall in some areas than others) is universal and is exhibited by almost every rainfall event.According to their characteristics, the rainfall distributions are divided into four types: downstream, midstream, upstream, and irregular (Figure 7).Four typical floods (one of each type) with extremely uneven spatial rainfall distributions were selected, and their simulation was assessed to verify the applicability of the model as such floods are difficult to simulate using hydrological models.The simulation results for these four typical floods are shown in Figure 8.The peak accuracies are 99.7%, 98.9%, 99.4%, and 99.7%, respectively.According to the simulation results for the four floods, the Liuxihe model has an excellent simulation effect under the condition of an uneven rainfall distribution, demonstrating the applicability of this distributed hydrological model.

Analysis of Model Parameter Optimization Results
To explore the influence of the PSO method on the inflow flood simulation accuracy of the Liuxihe model, the optimized and initial parameter values were used for the comparison of 24 floods, and average indicator values for the flood simulation results were obtained.All indicators were significantly improved after parameter optimization (Table 7).Specifically, the average Nash-Sutcliffe coefficient increased by 25.1% from 0.548 to 0.799, the average correlation coefficient increased by 8%, the average process relative error was reduced by 160%, and the average peak error was reduced by 38.8%.Thus, the results show that the PSO algorithm significantly improved the simulation accuracy of the Liuxihe model.This method requires only the selection of a typical flood type for parameter optimization to support practical application, effectively improving the model's performance.
were selected, and their simulation was assessed to verify the applicability of the model as such floods are difficult to simulate using hydrological models.The simulation results for these four typical floods are shown in Figure 8.The peak accuracies are 99.7%, 98.9%, 99.4%, and 99.7%, respectively.According to the simulation results for the four floods, the Liuxihe model has an excellent simulation effect under the condition of an uneven rainfall distribution, demonstrating the applicability of this distributed hydrological model.were selected, and their simulation was assessed to verify the applicability of the model as such floods are difficult to simulate using hydrological models.The simulation results for these four typical floods are shown in Figure 8.The peak accuracies are 99.7%, 98.9%, 99.4%, and 99.7%, respectively.According to the simulation results for the four floods, the Liuxihe model has an excellent simulation effect under the condition of an uneven rainfall distribution, demonstrating the applicability of this distributed hydrological model.

Analysis of Model Parameter Optimization Results
To explore the influence of the PSO method on the inflow flood simulation accuracy of the Liuxihe model, the optimized and initial parameter values were used for the com-

Conclusions
For exploring the effectiveness of the Liuxihe model in real-time flood forecasting for the Zaoshi Reservoir, this study employs a distributed hydrological model, the Liuxihe model, which is based on the DEM, land use data, and soil type data with a resolution of 90 m × 90 m, and uses the PSO algorithm to optimize the model parameters for flood simulation.Using the uneven spatial distributions of flood rainfall, flood prediction for the Zaoshi Reservoir was evaluated.From this study, the following conclusions can be drawn: 1.
Each rainfall event analyzed in this paper exhibits an uneven spatial distribution, and the floods can be divided into four types: downstream, midstream, upstream, and irregular.2.
In this study, the Liuxihe model's flood inflow forecasting scheme for the Zaoshi Reservoir showed an excellent simulation effect, with an average peak simulation accuracy of 96.3%, average peak time of 1.042 h early, and average Nash-Sutcliffe coefficient of 0.799.Under the condition of an uneven spatial distribution of rainfall, the Liuxihe model simulates floods well.Thus, flood simulation with the Liuxihe model is ideal for the watershed above the Zaoshi Reservoir, and the model can be applied in operational flood forecasting.

3.
The PSO algorithm significantly improved the simulation precision of the Liuxihe model, and its practical application requires only the selection of a typical flood for parameter optimization, thereby effectively improving the model's performance.
The basic data needed to construct the Liuxihe model are digital elevation model (DEM) data with spatial resolution of 30 m × 30 m, soil type data, and land use type data with spatial resolution of 1000 m × 1000 m.Aster V3 data, jointly developed and published by METI of Japan and NASA of the United States, are employed for the DEM (Figure 2a).The lowest elevation in the basin is 130 m, the highest elevation is 2240 m, and the average elevation is 696 m.Land use type data were obtained from the United States Geological Survey's Global Land Cover database ("http://landcover.usgs.gov"(accessed on 15 November 2022)), and soil type data were obtained from the International Food and Agriculture Organization (http://www.isric.org/(accessed on 15 November 2022)).
The basic data needed to construct the Liuxihe model are digital elevation model (DEM) data with spatial resolution of 30 m × 30 m, soil type data, and land use type data with spatial resolution of 1000 m × 1000 m.Aster V3 data, jointly developed and published by METI of Japan and NASA of the United States, are employed for the DEM (Figure 2a).The lowest elevation in the basin is 130 m, the highest elevation is 2240 m, and the average elevation is 696 m.Land use type data were obtained from the United States Geological Survey's Global Land Cover database ("http://landcover.usgs.gov"(accessed on 15 November 2022)), and soil type data were obtained from the International Food and Agriculture Organization (http://www.isric.org/(accessed on 15 November 2022)).To facilitate analysis, the data for the three underlying surfaces used in this study were resampled to a spatial resolution of 90 m × 90 m.

Sustainability 2023 , 14 Figure 3 .
Figure 3. Setting of third-order channel segment points in the Liuxihe model.To facilitate the estimation of river section size in the Liuxihe model according to the unit division and river classification results, as well as changes in the Google Earth remote sensing images and DEM of the basin above Zaoshi Reservoir, nodes are set, and virtual

Figure 3 .
Figure 3. Setting of third-order channel segment points in the Liuxihe model.

Figure 4 .
Figure 4. Non-adjustable model parameters of the Zaoshi Reservoir watershed: (a) flow direction and (b) slope.

Figure 5 .
Figure 5. Parameter optimization results: (a) change in the fitness value; (b) changes in parameter values; (c) process curves before and after parameter optimization.

Figure 5 .
Figure 5. Parameter optimization results: (a) change in the fitness value; (b) changes in parameter values; (c) process curves before and after parameter optimization.

Table 2 .
Statistical table for channel order division.

Table 3 .
Parameters of land use type.

Table 4 .
Parameters of soil type.

Table 5 .
Simulation evaluation indices for flood 20090607 before and after parameter optimization.

Table 6 .
Flood simulation index values.