One-dimensional hydrodynamic modeling approaches are useful for flood simulations; however, most studies often neglect intermediate discharges due to difficulties in obtaining the associated data. Herein, we produced the XAJ-H1DM model by coupling the Xinanjiang (XAJ) model, without the Muskingum module, with a one-dimensional hydrodynamic (H1DM) model, using regionalization approaches to test their practicality. Another model, named H1DM-XAJ, was also produced by orderly calibrating the H1DM and XAJ models to achieve improved flood simulations in poorly gauged catchments. The flood simulation capabilities of the four models (including the single XAJ and H1DM models) were investigated and compared at a daily time scale in the Three Gorges Reservoir Region, China. The results show that the regionalization approaches can be successfully used in the application of the integrated hydrologic and hydrodynamic model in ungauged intermediate catchments. Further, the coupled models produced markedly improved estimates of peak discharge and runoff volume compared to the single models. Moreover, the ability of the coupled models to simulate the peak water level and hydrograph, which hydrological models lack, is significantly better than that of the single H1DM model. The framework presented can be applied in other data-scarce catchments worldwide for better understanding of the hydrodynamic processes.
Keywords: Xinanjiang model; one-dimensional hydrodynamic model; flood simulation; Three Gorges Reservoir Region; ungauged catchments; integrated hydrologic and hydrodynamic model
Floods are one of the most serious natural disasters associated with water, causing severe property damage and loss of life. To reduce the impact of flooding, operational flood simulation systems have been developed around the world [[
Multi-model coupling integrates the selected model software to form a secondary processing integration process. The more thorough the secondary processing integration, the tighter the coupling simulation system integration. According to the ease of the secondary development in question, secondary processing integration can be divided into four categories: non-integrated, loosely coupled, tightly coupled, and fully integrated [[
Satisfactory hydrodynamic modeling of many receiving water bodies is limited because in situ observations of river discharge are scarce or unavailable [[
However, the simulation of runoff in ungauged basins is difficult (e.g., determining a suitable set of parameters without calibration) in hydrology [[
The first objective of this paper is, therefore, to assess the flood simulation capability of an integrated hydrologic and hydrodynamic model using regionalization approaches. The Xinanjiang model, which is the most popular conceptual rainfall-runoff model in China, has been extensively used for flood simulation and operational forecasting in gauged catchments since it was first developed in the 1970s. The flood events simulated by the Xinanjiang (XAJ) model are extremely sensitive to the parameters of the Muskingum runoff routing method [[
Second, we intend to address how to improve the flood simulation capability of the integrated hydrologic and hydrodynamic model for ungauged intermediate catchments. When the rainfall-runoff model is applied in a non-intermediate basin, where there is only one outlet in the entire basin, all runoff is generated by rainfall [[
Our specific objectives are: (i) to assess the flood simulation capability of the XAJ model and the H1DM model respectively; (ii) to evaluate the performance of the integrated hydrologic and hydrodynamic model using parameter replacement from interior basins; and (iii) to compare two coupled models to determine which is more appropriate in an ungagged intermediate basin.
The Zhutuo-Three Gorges Dam area, located in the middle reaches of the Yangtze River, China, is chosen as our study area (Figure 1). The Three Gorges Reservoir (Figure 1, SX), a typical River type reservoir with a modified end, is located approximately 40 km west of Yichang City in Hubei Province, central China (Figure 1, 106°53′–109°43′ E, 28°60′–30°56′ N). It is one of the largest reservoirs in the Yangtze River basin, with a water surface area of 1084 km
Figure 1 shows a schematic diagram of the ZT to Dam section of the Three Gorges Reservoir area. There are five hydrometric stations in study area, Zhutuo (ZT), Cuntan (CT), Qingxichang (QXC), Wanxian (WX), and Sanxia (SX), from upstream to downstream. The storage control stations of the Jialing River and Wujiang River, two major tributaries of the reservoir area, are the Beibei hydrometric station (BB) and the Wulong hydrometric station (WUL), respectively.
The landscape characteristics of all catchments were derived from a DEM with a spatial resolution of 90 m (3") provided by Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (
The river section data includes 368 sections from the ZT to the Three Gorges Dam, 30 sections from the BB to the Yangtze River mainstream, and 46 sections from the WUL Station to the Yangtze River mainstream. In addition, there are 12 tributaries with measured section data and a total of 144 sections.
There are several hydrometric stations in the study area. According to the location distribution of key hydrometric stations in mainstream, the study area is divided into four sub-regions, namely ZT-CT, CT-QXC, QXC-WX, and WX-SX. The dividing line is determined according to the topography of the natural watershed. The areas of these four sub-regions are: 14,555 km
The Xinanjiang model was first used to predict the storage of the Xinanjiang Reservoir and later became a general-purpose rainfall-runoff model. Its main feature is the formation of runoff as a dependent variable of the storage recharge, that is, no runoff is generated before the soil water content in the aeration zone reaches the field water holding capacity. Thereafter, the runoff is equal to the excess rainfall without further loss. The Xinanjiang model, with four runoff components, has been widely used in humid or semi-humid regions in China, as well as in many other countries. The overall structure of the Xinanjiang model is shown in Figure 2.
For each sub-catchment, the outflow is calculated by four major modules:
- Evapotranspiration module: the three-layer evapotranspiration mode—that is, the total evaporation is composed of three parts, namely, surface evaporation, shallow evaporation, and deep evaporation.
- Runoff generation module: runoff generation under saturated conditions—that is, the precipitation does not produce flow before the field water capacity is satisfied, and all precipitation is absorbed by the soil. After the precipitation meets the field water holding capacity, all precipitation (excluding the evaporation in the same period) produces flow.
- Runoff separation module: the division of the three water sources—that is, according to the free water capacity distribution curve, the total production flow is divided into the surface runoff, soil middle flow, and underground runoff.
- Runoff routing module: this module is divided into the river network confluence and the river channel confluence. The river network confluence is runoff that flows directly into the sub-catchment outlet; the river confluence is the runoff at the sub-catchment outlet that uses the Muskingum method to evolve the basin exit.
The 1-D hydrodynamic model (H1DM) was used to simulate the hydrodynamic processes of the study area. One-dimensional mathematical models have been widely used to solve several engineering problems and have obtained many satisfactory results in China and abroad [[
(
(
where Q is the flow rate, q is the lateral inflow, A is the wetted area, η is the water level, g is the gravity acceleration, x is the river mileage coordinate, S
The friction slope S
(
where R is the hydraulic radius, n is the Manning coefficient, and u is the water velocity.
Equation (
(
In this study, we developed the H1DM to simulate floods. The finite volume method and the finite element method were used to discretize the hydrodynamic model control equations, and the ELM (Eulerian–Lagrangian method) was used to solve the convection term in the momentum equation [[
According to different coupling methods, two hydrological-hydrodynamic coupled models are proposed. Under these two coupling modes, the Xinanjiang model is used without the Muskingum module. The coupled models share 15 of the 17 Xinanjiang model parameters. The KX and EX parameters of the Muskingum method in the original Xinanjiang model are not needed in the coupled model. Instead, the roughness coefficient of the H1DM is introduced in the coupled model.
According to the two assumptions of the Muskingum method, the flow in a section of a river and the water demand of the river section has a single linear relationship, and the water surface line of the river section is a straight line; however, in a natural river channel, these two linear relationships are difficult to satisfy [[
We used the Xinanjiang model without the Muskingum method. The confluence of the Xinanjiang model can be divided into the river network confluence and the river channel confluence. When the watershed area is less than 1000 km
Based on the H1DM, the XAJ-H1DM model is developed using the intermediate flow predicted by the Xinanjiang model as the intermediate boundary conditions.
This paper developed the XAJ-H1DM in a loosely coupled way. First, the Xinanjiang model of the sub-catchment in the sub-region are calibrated with observed data and the calibrated parameters were transposed to the catchment in the same sub-region. In the receiving catchment, the hydrology data were gauged. In addition, the one-dimensional hydraulic model was calibrated. Then, the prediction results of the Xinanjiang model are appropriately increased according to the area ratio of the sub-catchment and sub-region. Finally, the increased results are added to the one-dimension hydrodynamic model as the boundary condition of the one-dimensional hydrodynamic model. The flow chart of the XAJ-H1DM and H1DM-XAJ are shown in Figure 3.
In this coupled model, the upper boundary uses the measured flow from the ZT, BB, and WUL stations, the lower boundary uses the measured water levels in front of the dam during the same period, and the internal boundary uses the results of the Xinanjiang model forecast (Figure 4a).
The relative position relationship of sub-catchments and sub-regions is shown in the Figure 5.
In contrast to the previous coupling method, in which the Xinanjiang model and the H1DM model were calibrated independently, in the coupling mode described here, the calibration of the Xinanjiang model was premised on the results of the H1DM model.
First, the H1DM model was separately calibrated, and then the optimization algorithm was used to calibrate the Xinanjiang model, of which the results of the calibrated H1DM model were used as the basis for the coupled model. The results of the one-dimensional hydraulic model used here are fixed and the only the parameters of the Xinanjiang model were calibrated. Finally, the results of the Xinanjiang model are added to the one-dimensional hydrodynamic model as the intermediate flow of the one-dimensional hydrodynamic model. The flow chart of the H1DM-XAJ is shown in Figure 3.
In this coupled model, the Xinanjiang model is used to predict the discharge values of the four sub-regions, which are added to the H1DM model as the intermediate inflow values (Figure 4b).
To compare the performance of the Xinanjiang model to that of the coupled model, four sets of simulations were implemented at a daily time scale: (
Daily flood events from 2007 to 2008 for each catchment were used for calibration, while the remaining events (daily flood events from 2009) were used for validation. In the period of calibration and validation, the Xinanjiang model was applied to each of the six study catchments and calibrated separately in each individual catchment. All parameters changed within their "reasonable" range based on the literature and/or experience with the model implementation [[
The shuffled complex evolution (SCE-UA) global optimization algorithm [[
(
where N
In accordance with the accuracy standard established by the MWR (Ministry of Water Resources of the People's Republic of China, 2008) [[
(
(
(
where R
The simulation time of the hydrodynamic model is set to one year and the time step is set to 900 s. There are three inflow inputs in this model, the ZT hydrometric station flow, BB hydrometric station flow, and WUL hydrometric station flow. In the downstream of the model, the discharge of the model is controlled by the water level of the Fenghuangshan hydrometric station, which is the water level control station of the three gorges reservoir.
The roughness parameter plays a significant role in the hydrodynamic model. Although their values are priori requirements for the solution of the Saint-Venant equations, there have been less definitive measurements of these values. Accurate roughness values are not well established [[
A trial-and-error method was used to calibrate model parameters. First, the operator selected a set of parameters based on practical experience or existing literature, which was used to drive the model. Then, the simulated value of water level was compared with the measured value in the section that needed to be calibrated. If the measured value was higher than the simulated value, the roughness value was increased; otherwise, the roughness value was decreased. The hydrometric station divides the river into several sections whose roughness values are calibrated in order from downstream to upstream.
The H1DM was calibrated to minimize the mean error of the water level estimation, and the calibrated parameters are listed in Table 1. The simulation results obtained were close to those of previous studies [[
Hydrologic modeling and parameter replacement were combined to obtain runoff values for this poorly gauged basin. The hydrological similarity among the sub-basins was first analyzed using hierarchical cluster analysis (HCA) [[
In this study, the topographic index was chosen as the defining characteristic of the sub-basins, and in the upstream sub-basins this value was extracted from the DEM. The topographic index, also known as the topographic wetness index (TWI), was developed by Beven and Kirkby [[
The single flow direction (SFD) algorithm [[
Therefore, based on the results from HCA, the topographic index probability density distribution curves were also compared for different sub-basins. Finally, parameter replacement was performed for groups with hydrological similarity.
The performance of the XAJ-H1DM was investigated to determine whether improved simulations were obtained using the parameters transposed from the donor catchments, and the performances of H1DM-XAJ, XAJ-H1DM, and XAJ models were compared to identify whether improved simulations were obtained from the H1DM-XAJ for ungauged intermediate catchments.
The calibrated parameter values of the Xinanjiang module for the three models are listed in Table 2. There are 16 parameters in the Xinanjiang model, while the Muskingum parameters (KE and XE) not included in the XAJ-H1DM and H1DM-XAJ models. The outflow is finally routed by the H1DM instead of the Muskingum successive-reaches model in the latter two models to produce the flow at the outlet of the whole catchment. The XAJ-H1DM is applied to sub-catchments, while the remaining models are applied to the sub-regions.
Four evaluation indexes involving the relative runoff volume error (RRE, %), relative peak discharge error (RPE, %), peak time error (PTE, day), and Nash–Sutcliffe efficiency (NSE) are used to evaluate the simulation accuracy of the flow for the four models. The accuracy statistics of flow for the four models are summarized in Table 3, which shows the absolute means with respect to the RRE, RPE, PTE, and NSE values for the calibration and validation periods.
For the Xinanjiang model, the absolute means of the RRE, RPE, and PTE values in the simulations for the four study regions range from 1.01% to 4.8%, −4.75% to 12.24% and 0 to −1 day, respectively (Table 2). The PTE values of the Xinanjiang model simulations are mostly zero at a daily time scale except in one special case (Table 2 result of SX in 2008). The NSE values vary between 0.9776 and 0.9937. These results indicate that the Xinanjiang model is capable of reproducing both the magnitude and dynamics of the observed flood events at different catchment scales once its parameters are calibrated.
For the H1DM model, the absolute means of the RRE and RPE simulations for the four study regions range from −10.19% to −2.3% and −20.56% to 0.78%, while the PTE is always 0 days (Table 2). The NSE varies between 0.9212 and 0.9960 (Table 2), indicating that the H1DM model can reasonably simulate the general shapes of the storm hydrographs when calibrated singly. This suggests that the H1DM model is highly capable of simulating the total runoff volume and has an inferior capability to simulate the peak discharge when only the two main tributaries were taken into consideration in these catchments. However, the original Xinanjiang model shows relatively better performance in simulating the total runoff volume and the peak values of discharge in all catchments and periods, since both the main tributaries and rainfall-runoff values in the catchments are reasonably simulated in the original Xinanjiang model.
Three additional evaluation indexes, the absolute mean water level error (AMWE, m), relative maximum water level error (RMWE, %), and Nash–Sutcliffe efficiency (NSE), are used to evaluate the simulation accuracy of water level for the four models, as follows:
(
(
where
The accuracy statistics of the water level for the three models' performances are summarized in Table 4, which shows the absolute means with respect to the AMWE, RMWE, and NSE values for the calibration and validation periods. Meanwhile, the AMWE and RMWE simulations for the four study regions range from 0 m to −0.34 m and 0% to −2.19%, respectively. The NSE values of the water level vary between 0.9839 and 1 (Table 4). The water level simulations and observations are almost identical in SX because the water level was set as the downstream boundary of the H1DM model while SX was too close to the downstream boundary to avoid the influence.
Relative to the flow simulations of H1DM model, the water level simulations show a superior performance in NSE since the roughness, the only parameter calibrated in the H1DM model, only affects the water level simulations. In other words, the water level simulations can be adjusted to the best state but the flow simulations cannot be adjusted once the boundary conditions of H1DM are confirmed. The roughness values are listed in Table 1. Table 3 shows that the validated roughness values are basically consistent with those of previous studies. These results indicate that the H1DM model is capable of reproducing both the flow and water levels of observed flood events once its parameters are calibrated.
The physical characteristics among the 6 sub-basins and the 4 sub-regions upstream of the Three Gorges Reservoir Basin showed distinct differences. The topographic index probability density distribution curves were compared among the 10 regions (Figure 6a). Based on the topographic index probability density distribution and the locations, the 10 regions were divided into three categories.
Finally, three groups with hydrological similarities were obtained from the results of the topographic index probability density distribution curve, and parameter replacement could then be conducted within each group.
For the XAJ-H1DM model, the ungauged catchments were modeled using the calibrated parameter sets from similar catchments in the same sub-region.
The donor catchments for the receptor catchments listed in Table 5 were selected according to their locations (see Figure 4) and topographic index probability density distribution. For ZT-CT sub-region and CT-QXC sub-region, model parameters were transposed from a single donor catchment, i.e., WC. For QXC-WX sub-region, model parameters were transposed from SZ catchment, while multiple donor catchments were used for WX-SX sub-region. The inflows of the hydrodynamic model in the coupled model increased from 3 to 15, while the output controlled by the water level remained the same.
The calibrated parameter values of the Xinanjiang module of XAJ-H1DM for each sub-catchment are listed in Table 2. The accuracy statistics of the model performance are summarized in Table 5, which represents the annual NSE values and the average NSE values for the separate calibration and validation flood events. For all sub-catchments, since all watersheds are mesoscale, and from small watersheds and the Xinanjiang module was only used to simulate runoff, the Muskingum module in runoff routing was not included in the Xinanjiang module.
According to Table 5, the NSE values for SZ were 0.291, 0.223 and 0.235 from 2007 to 2009. There are two possible explanations for this result: the catchment area is too small, and the outflow of the catchment during this period is too small (less 350 m
In addition to the above conditions, the Xinanjiang module performs well, and NSE is higher than 0.5 for all remaining flood events. For all calibration and validation events, the average NSE values for the six study sub-catchments range from 0.249 to 0.781. In general, the simulations of the Xinanjiang module can be used to simulate rainfall-runoff accurately. According to the relative position and the nesting relationship of the basin, the parameters of the Xinanjiang module in sub-catchments are simply transposed to the sub-regions.
Table 3 and Table 4 provide the accuracy statistics of the XAJ-H1DM model results by coupling the singly calibrated Xinanjiang model with the singly calibrated H1DM model. The simulations by the XAJ-H1DM model show obvious improvement for the RRE and RPE metrics relative to those of the single H1DM model. The RRE and RPE tremendously improved at SX in 2007, form −10.19% to −2.28% and from −20.56% to 14.38% (Table 2). It is not surprising that the performance of XAJ-H1DM model has a remarkable improvement than H1DM model in all aspects. For the results of water level (Table 4), relative to H1DM model, both the simulations of AWME and the simulations of RMWE have decreased for XAJ-H1DM in varying degrees, while their NSE increased accordingly.
Collective analysis of the simulation results shown in Table 3 and Table 4 clarifies that the XAJ-H1DM model outperforms the H1DM model in most cases and can be compared with Xinanjiang model. Moreover, extra simulations, such as water level and water surface profile in any reach, can be obtained from the XAJ-H1DM model. Therefore, it is evident that the integrated hydrologic and hydrodynamic model has the ability to simulate floods using the regionalization approaches.
To address the issue of the simulation results of the XAJ-H1DM model not being satisfactory with regard to ungauged intermediate catchments, another coupled model (H1DM-XAJ model) is proposed. For the H1DM-XAJ model, four Xinanjiang models established to predict intermediate flow were applied to ZT-CT, CT-QXC, QXC-WX, and WX-SX sub-regions. The calibration of H1DM module in H1DM-XAJ model and the calibration of a single H1DM model are the same. The calibration of Xinanjiang module in H1DM-XAJ model is different from the calibration of Xinanjiang module in XAJ-H1DM model and the single Xinanjiang model. In other words, there are four different Xinanjiang parameter sets for all four models. The calibration method and the calibration target are the same as the previous model.
In this coupling mode, the calibration of the Xinanjiang model is based on the results of the H1DM model. The calibrated parameter values of the Xinanjiang module of H1DM for each sub-catchment are listed in Table 2. The accuracy statistics of the H1DM-XAJ model performance are summarized in Table 3 and Table 4. The absolute means of RRE and RPE simulations for the four study regions range from −1.51% to 5.63% and −3.52% to 6.97%, while the PTE is always 0 day (Table 3). The NSE varies between 0.9760 and 0.9976. From Table 4, it is evident that the AMWE and RMWE simulations for the four study regions range from 0 m to 0.237 m and 0% to 1.38%, respectively. These indicate that the H1DM-XAJ model shows a veracious simulation and can be compared with other models.
In addition, the water level and discharge curves of several stations are shown in Figure 7 that allow the comparison of the results of the coupled model. As can be seen from the figure, H1DM-XAJ performs well in most cases, and the simulation accuracy of water level and flow is higher than that for other models. Especially regarding the peak simulation, the H1DM-XAJ model is better than other models.
To compare the simulation results of the four models more intuitively, the contents of Table 3 and Table 4 are rearranged and are shown as Figure 8, Figure 9 and Figure 10. The absolute values of all metrics were taken and were averaged with respect to time. These were then plotted according to the sequence of the hydrological stations. As can be seen from Figure 8, RRE of H1DM-XAJ model is significantly smaller than that of other models in CT and QXC, while XAJ-H1DM model replaces the position of H1DM-XAJ in WX and SX. The difference between H1DM-XAJ and XAJ-H1DM is the simulations of rainfall-runoff in Xinanjiang module. Therefore, the simulations of rainfall-runoff were overrated by H1DM-XAJ model in WX and SX.
According to Figure 8, the performance of RPE in these four models shows a similar outcome. Meanwhile, the simulations in terms of NSE are larger than 0.98 for all sub-regions except for the single H1DM model (Figure 9). According to Figure 8b, H1DM-XAJ and XAJ-H1DM show worse results in AMWE. The reason for this may be that the simulations of H1DM model have been calibrated to the best status before the simulations of the Xinanjiang model were taken into account. Therefore, the addition of simulations of the Xinanjiang model made the average water level too high. On the contrary, with the addition of the Xinanjiang model, the ability of the coupled models to simulate the flood peak improved (Figure 8b), and the fitting degree of water level process line also increased (see Figure 9b). These show that the Xinanjiang model, the XAJ-H1DM model, and the H1DM-XAJ model, once calibrated with parameters, were able to reproduce the magnitude and dynamics of observed flood events on different catchment scales. Furthermore, the H1DM-XAJ model performed better than the other three in majority of the instances and showed obvious improvement compared with the other two.
Table 6 presents the ratio of the intermediate inflow calculated by the two coupled models to the measured flow at the downstream site. The intermediate inflow is relatively small, and the intermediate inflow has little effect on the simulation results, which explains why the NSE values of the two models do not change much. In addition, from upstream to downstream of the river basin, the proportion of intermediate inflow gradually increased, which may be related to the steep terrain of the downstream river basin and the amount of rainwater.
The above discussion is based on the results obtained from a simple comparison of the evaluation indexes. We further discuss the model performance by comparing the flood simulation capabilities of the four models for different catchments. In this approach, the correlation coefficient R
The simulation results of the four models are distributed around the measured values, and the simulation results of the coupled models are more densely distributed than the single Xinanjiang model and the H1DM model. In addition, relative to XAJ-H1DM model, simulation results of H1DM-XAJ model are more uniformly and intensively distributed (see Figure 11 and Figure 12). It is worth mentioning that the H1DM model shows several extreme anomalies in SX, which may be due to the fact that SX is too close to the Three Gorges Dam, and H1DM does not consider the huge inflows of tributaries contributing to the mainstream flow.
In addition, the simulations of the H1DM-XAJ model are much better than those of the other three models. The Xinanjiang model tends to overestimate the flow value, while the XAJ-H1DM model does the opposite; the latter always simulates a smaller simulation value than the real data. H1DM has a more balanced simulation capability, and simulation results are always evenly distributed around the observed values. These conclusions justify the results presented in Table 3 and Table 4. In terms of water level simulation results, compared with the single H1DM model, the coupled models also showed a significant improvement (see Figure 11 and Figure 12).
The coupled models perform better when simulating middle water levels relative to simulated high and low water levels, and the simulation results are closer to the real measurements when simulating middle water levels (Figure 12a). In addition, hydrodynamic models tend to simulate smaller results than the actual values at the CT station. The possible reason for this phenomenon is that CT is located near the end of the backwater fluctuation area of the Three Gorges Reservoir; the riverbed is silted, which results in the deviation between the topographic data and the measured data. Based on the results presented here and the aforementioned similarity analysis of flow and water level, it is suggested that the coupled model, developed by coupling Xinanjiang model with H1DM model, is more likely suitable for simulation of flood in the Three Gorges Reservoir, and has the potential to produce substantially improved flood simulations in this ungauged intermediate catchment.
The calculations of the XAJ-H1DM and H1DM-XAJ models are almost the same. However, the structures of the two models are different, with H1DM-XAJ being better than XAJ-H1DM in terms of model structure and simulation results. In addition, the calculation volume of the two integrated models is larger than that of the single Xinanjiang model, requiring more data and outputting more results. Coupling the hydrological model with the hydrodynamic model can overcome the shortcomings of the single hydrological model that cannot predict the water level and the lack of accuracy of the single hydrodynamic model due to the lack of intermediate inflow simulation. Moreover, it combines the advantages of both models and takes rainfall as well as terrain information into consideration.
The models proposed in this paper integrate hydrological and hydrodynamic models in the simplest coupling way. Compared with other fully integrated distributed models, such as the coupled model of MIKE11 and MIKE SHE [[
- in the simulation of rainfall-runoff, only the effects of rainfall and evaporation on the runoff generation and routing in the basin were considered, and no consideration was given to meteorological factors (such as humidity, wind direction), land use, and vegetation cover;
- compared with the distributed hydrological model, the Xinanjiang model is faster in the calculation and requires less data, but the physical meaning is not clear enough and there are more empirical parameters;
- the H1DM used in this paper has fast calculation speed and high accuracy but does not consider the different roughness of the riverbed and the flood plain.
However, the coupling way proposed in this paper has proven to be effective, and the integrated model also has reliable accuracy that can be applied to flood simulation.
One-dimensional hydrodynamic models are very useful in simulating flood routing in natural rivers, but most studies often do not take intermediate inflows into account because of the difficulty in obtaining the data regarding intermediate flows. In view of the deficiencies of the general one-dimension hydrodynamic model, we proposed its coupling with the Xinanjiang model and produce two coupled models, namely the XAJ-H1DM and H1DM-XAJ models, which differed in their manner of coupling.
For the XAJ-H1DM model, the Xinanjiang model and the H1DM model were separately calibrated, and the results of the Xinanjiang model in the coupled model were derived from the Xinanjiang model for the enlarged sub-catchment; the results were verified by a small amount of gauged data. In the H1DM-XAJ model, H1DM model was calibrated first, and then the Xinanjiang model was calibrated on the basis of the already calibrated H1DM model for each catchment. The four models were applied to a poorly gauged catchment, Three Gorges Reservoir Region, where the intermediate runoff cannot be obtained directly but can be predicted by the rainfall-runoff model.
In this study, we evaluated the flood simulation capability of the Xinanjiang model, H1DM model and two coupled models at the daily time scale for the purpose of improving the capability of flood simulation in ungauged intermediate catchments.
- Results indicate that the regionalization approaches can be successfully used in the application of the integrated hydrologic and hydrodynamic model.
- Our results show that both coupled models are capable of providing satisfactory and comparative simulations of runoff volume, peak discharge, peak time, and flood hydrograph for the data-scarce catchments once they are well-calibrated. Integrated models also have the ability to simulate floods using the regionalization approaches. The coupled models produced markedly improved estimates of peak discharge and runoff volume as compared to the single H1DM model, indicating that the intermediate flow is a major uncertainty source for the application of the H1DM model in large scale watersheds. Coupling Xinanjiang with H1DM has the potential to substantially improve the flood simulation capability of the H1DM model in poorly gauged catchments.
- The coupled models have shown improvements in peak discharge, runoff volume, peak time, and hydrograph as compared to the XAJ model. Moreover, the ability of the coupled models to simulate peak water level and hydrograph, which any hydrological model lack, is significantly better than that of the single H1DM model. This study demonstrates the importance of incorporating intermediate inflows, which can be obtained from rainfall-flow model predictions in data-deficient areas, in hydraulic models.
- The coupled way of the H1DM-XAJ model provides a realizable direction to improve the flood simulation capability of the integrated hydrologic and hydrodynamic models in ungauged intermediate catchments.
The presented framework can be applied to other data-scarce catchments worldwide to gain a better understanding of the hydrodynamic processes.
MAP: Figure 1 Map of the study area and locations of the river gauging stations. ZT, CT, QXC, WX, SX, BB and WUL stands for the Zhutuo, Cuntan, Qingxichang, Wanxian, Sanxia, Beibei and Wulong hydrometric station respectively.
Graph: Figure 2 Flow chart of the Xinanjiang model.
Graph: Figure 3 Flow chart of the Xinanjiang model with the one-dimension hydrodynamic model, XAJ-H1DM (left) and H1DM-XAJ (right).
Graph: Figure 4 Schematic diagram of XAJ-H1DM (a) and H1DM-XAJ (b).
Graph: Figure 5 Relative positioning of the sub-catchments and sub-regions. WC, SZ, WQ, CTA, WXI and XS stand for the Wucha, Shizhu, Wenquan, Changtan, Wuxi, and Xingshan hydrometric station respectively.
Graph: Figure 6 Topographic index probability density distribution curve for six sub-basins in different groups. The topographic index probability density distribution curves of the 10 regions (a), the first group (b), the second group (c), the third group (d).
Graph: Figure 7 Streamflow hydrographs for several sites. (a) CT-2009-Flow; (b) CT-2009-Water level; (c) QXC-2007-Flow; (d) QXC-2007-Water level; (e) QXC-2009-Flow; (f) QXC-2009-Water level.
Graph: water-770896-g007b.tif
Graph: water-770896-g007c.tif
Graph: water-770896-g007d.tif
Graph: Figure 8 Average RRE (a) and AWME (b) of the four sites for all periods.
Graph: Figure 9 Average RPE (a) and RWME (b) of the four sites during all periods.
Graph: Figure 10 Average NSE (a) and NSE-L (b) of the four sites for all periods.
Graph: Figure 11 Flow scatter diagram and correlation coefficient (R2) of the four models in CT (a), QXC (b), WX (c), SX (d) for all periods.
Graph: Figure 12 Water level scatter diagram and R2 of three models in CT (a), QXC (b), WX (c), SX (d) for all periods.
Table 1 Roughness of H1DM.
Subdomains Range of Roughness ZT-CT 0.025–0.028 CT-QXC 0.028–0.033 QXC-WX 0.033–0.035 WX-SX 0.035–0.056
Table 2 List of parameters of the models and the calibrated values for the study catchments.
Xinanjiang Model H1DM-XAJ Model XAJ-H1DM Model Module Para ZT-CT CT-QXC QXC-WX WX-SX ZT-CT CT-QXC QXC-WX WX-SX WC SZ WQ CTA WXI XS Evapotranspiration K 0.500 0.501 1.100 0.953 0.500 0.500 0.500 0.500 0.500 0.500 1.002 0.701 0.587 0.501 UM 65 75 74 74 65 73 70 71 13 13 26 21 29 16 LM 76 77 77 61 85 64 65 62 66 66 90 74 79 71 DM 31 26 37 18 60 15 16 15 15 15 49 15 60 40 C 0.16 0.11 0.15 0.10 0.14 0.15 0.12 0.13 0.13 0.13 0.11 0.09 0.12 0.16 Runoff generation IM 0.00 0.03 0.01 0.01 0.03 0.03 0.02 0.01 0.007 0.007 0.030 0.000 0.029 0.023 B 0.1 0.2 0.2 0.3 0.4 0.2 0.2 0.3 0.1 0.1 0.3 0.4 0.5 0.3 Runoff separation SM 44.00 16.50 40.81 10.04 22.02 16.72 10.00 10.01 4.55 4.55 1.00 6.93 1.00 8.63 EX 0.50 0.50 1.16 2.00 0.50 0.50 0.50 0.50 0.66 0.66 0.50 0.50 0.50 0.50 KI 0.376 0.449 0.350 0.390 0.425 0.380 0.429 0.426 0.450 0.450 0.406 0.450 0.350 0.450 KG 0.305 0.250 0.350 0.250 0.250 0.250 0.250 0.250 0.252 0.252 0.307 0.350 0.269 0.250 Runoff routing CS 0.500 0.412 0.164 0.500 0.408 0.500 0.495 0.469 0.410 0.410 0.269 0.273 0.545 0.467 CI 0.567 0.900 0.754 0.604 0.503 0.565 0.500 0.559 0.691 0.691 0.790 0.900 0.726 0.898 CG 0.992 0.990 0.998 0.998 0.994 0.990 0.990 0.990 0.998 0.998 0.998 0.996 0.998 0.998 KE 13.289 2.877 11.324 0.214 - - - - - - - - - - XE 0.282 0.000 0.122 0.003 - - - - - - - - - -
Table 3 Accuracy statistics of the four model flow simulations with calibrated parameters for both calibration and validations events.
2007 2008 2009 XAJ H1DM XAJ-H1DM H1DM-XAJ XAJ H1DM XAJ-H1DM H1DM-XAJ XAJ H1DM XAJ-H1DM H1DM-XAJ CT RRE (%) 1.70% −2.39% −0.20% −0.51% 1.76% −2.54% −0.83% −0.67% 2.34% −3.94% −2.19% −1.51% RPE (%) −4.75% 0.78% 1.17% 1.05% −3.08% −3.53% −3.40% −3.15% 8.77% −10.95% −1.92% −1.32% PTE (day) 0 0 0 0 0 0 0 0 0 0 0 0 NSE 0.9891 0.9960 0.9975 0.9976 0.9909 0.9949 0.9976 0.9976 0.9899 0.9925 0.9975 0.9978 QXC RRE (%) 3.06% −5.84% −3.92% −1.14% 3.42% −6.40% −4.80% −1.38% 2.82% −5.56% −3.68% −0.57% RPE (%) 2.12% −9.24% −8.09% −3.52% 3.12% −5.18% −4.39% −0.88% 1.46% −15.25% −5.98% −1.60% PTE (day) 0 0 0 0 0 0 0 0 0 0 0 0 NSE 0.9910 0.9810 0.9883 0.9955 0.9900 0.9774 0.9827 0.9896 0.9915 0.9776 0.9867 0.9913 WX RRE (%) 1.01% −4.68% −1.84% 2.16% 1.03% −3.37% −0.87% 3.92% 1.05% −3.18% −0.47% 3.88% RPE (%) −2.90% −10.77% −8.42% −0.72% 0.31% −2.40% 0.56% 6.94% −3.77% −13.69% −3.87% 2.09% PTE (day) 0 0 0 0 0 0 0 0 0 0 0 0 NSE 0.9930 0.9834 0.9921 0.9951 0.9936 0.9885 0.9936 0.9886 0.9937 0.9864 0.9935 0.9893 SX RRE (%) 4.80% −10.19% −2.28% 3.70% 4.27% −8.95% −2.09% 5.63% 3.58% −7.28% −0.85% 3.20% RPE (%) 6.42% −20.56% −14.38% −2.89% 12.24% −7.75% 5.27% 6.07% 0.20% −1.24% 1.76% 8.92% PTE (day) 0 0 0 0 −1 0 0 0 0 0 0 0 NSE 0.9776 0.9212 0.9766 0.9915 0.9808 0.9426 0.9871 0.9788 0.9895 0.9524 0.9803 0.9750
Table 4 Accuracy statistics of the four model water level simulations with calibrated parameters for both calibration and validations events.
2007 2008 2009 XAJ H1DM XAJ-H1DM H1DM-XAJ XAJ H1DM XAJ-H1DM H1DM-XAJ XAJ H1DM XAJ-H1DM H1DM-XAJ CT AWME (m) - −0.340 −0.204 −0.237 - −0.215 −0.112 −0.099 - −0.114 −0.028 0.003 RMWE (%) - −0.76% −0.74% −0.74% - −0.64% −0.60% −0.60% - −2.19% −1.37% −1.38% NSE - 0.9839 0.9885 0.9867 - 0.9922 0.9942 0.9943 - 0.9890 0.9938 0.9940 QXC AWME (m) - −0.046 0.046 0.132 - 0.017 0.083 0.180 - −0.044 0.015 0.083 RMWE (%) - −1.09% −0.88% −0.36% - −0.13% 0.03% 0.05% - −0.04% −0.03% −0.02% NSE - 0.9896 0.9934 0.9920 - 0.9993 0.9990 0.9977 - 0.9991 0.9995 0.9990 WX AWME (m) - −0.046 0.026 0.064 - 0.007 0.057 0.096 - −0.009 0.030 0.061 RMWE (%) - −0.02% −0.01% 0.02% - −0.01% 0.06% 0.07% - −0.01% 0.00% 0.01% NSE - 0.9966 0.9981 0.9974 - 0.9998 0.9996 0.9993 - 0.9998 0.9987 0.9995 SX AWME (m) - 0.000 0.000 0.000 - 0.000 0.000 0.000 - 0.000 0.000 0.000 RMWE (%) - 0.00% 0.00% 0.00% - 0.00% 0.00% 0.00% - 0.00% 0.00% 0.00% NSE - 1.0000 1.0000 1.0000 - 1.0000 1.0000 1.0000 - 1.0000 1.0000 1.0000
Table 5 Accuracy statistics of the Xinanjiang module in XAJ-H1DM.
Catchment NSE 2007 2008 2009 Average WC 0.514 0.799 0.826 0.713 SZ 0.291 0.223 0.235 0.249 WQ 0.848 0.662 0.832 0.781 CTA 0.741 0.741 0.590 0.690 WXI 0.830 0.330 0.394 0.518 XS 0.640 0.642 0.602 0.628
Table 6 Ratio of intermediate flow to the gaged inflow.
Catchment 2007 2008 2009 ZT-CT 2.22%/1.5% 1.83%/2.09% 1.88%/2.55% CT-QXC 2.08%/3.94% 1.97%/5.26% 2.12%/5.38% QXC-WX 3.15%/5.94% 2.85%/7.45% 2.92%/7.38% WX-SX 8.06%/12.02% 7.07%/14.47% 6.64%/13.85%
All listed authors have contributed substantially to the work, namely: Conceptualization, J.Z. and Y.Z.; Methodology, Y.Z.; Software, Y.Z.; Validation, Y.Z. and C.L.; Formal Analysis, Y.Z.; Investigation, Y.Z. and C.L.; Resources, J.Z.; Data Curation, Y.Z.; Writing—Original Draft Preparation, Y.Z. and C.L.; Writing—Review & Editing, Y.Z.; Visualization, Y.Z.; Supervision, J.Z.; Project Administration, J.Z.; Funding Acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.
This study was funded by the National Natural Science Foundation of China, grant number U1865202 and 91547208.
The authors declare no conflict of interest.
By Yulong Zhang; Jianzhong Zhou and Chengwei Lu
Reported by Author; Author; Author