A Coupled Hydrologic–Hydraulic Model (XAJ–HiPIMS) for Flood Simulation
To protect ecologies and the environment by preventing floods, analysis of the impact of climate change on water requires a tool capable of considering the rainfall-runoff processes on a small scale, for example, 10 m. As has been shown previously, hydrologic models are good at simulating rainfall-runoff processes on a large scale, e.g., over several hundred km2, while hydraulic models are more advantageous for applications on smaller scales. In order to take advantages of these two types of models, this paper coupled a hydrologic model, the Xinanjing model (XAJ), with a hydraulic model, the Graphics Processing Unit (GPU)-accelerated high-performance integrated hydraulic modelling system (HiPIMS). The study was completed in the Misai basin (797 km2), located in Zhejiang Province, China. The coupled XAJ–HiPIMS model was validated against observed flood events. The simulated results agree well with the data observed at the basin outlet. The study proves that a coupled hydrologic and hydraulic model is capable of providing flood information on a small scale for a large basin and shows the potential of the research.
Keywords: Xinanjiang model; high-performance integrated hydraulic modelling system (HiPIMS); coupled hydrologic–hydraulic model; two-dimensional surface confluence calculation over basin
1. Introduction
Flood, as the most common natural hazard, causes huge causalities, heavy economic losses, and ecological and environmental risks worldwide [[1]]. A warmer climate creates more extreme storm events, combined with rapid land use/cover change and a growing population, the risk of floods to urban and farming areas is increasing and posing a threat to ecologies and the environment [[2]]. Extreme rainfall is considered as the main cause of flooding [[7]], however, inadequate flood risk management also contributes to the increase of flood risk [[8]]. There is an urgent requirement for flood simulation/forecasting and integrated flood risk management in our changing environment.
Hydrologic models are commonly implemented for rainfall-runoff simulations to predict the basin-scale flow of water passing through the basin outlet, normally on a daily/hourly temporal scale, e.g., Soil and Water Assessment Tool (SWAT) [[9]], Xinanjing (XAJ) [[11]], Stanford Watershed Model (SWM) [[13]], Sacramento Model (SAC) [[14]]. However, high spatial heterogeneity of a realistic basin raises the requirement for fine flood management characterized by spatial information inside of the complex basin that affects its response to extreme storms. Numerous distributed grid-based hydrologic models have been developed to represent spatial variability [[15]]. River channel flow is computed by the Muskingum method, e.g., [[15]], or the Muskingum method is coupled with a one-dimensional (1D) hydraulic model to generate water depth along a river channel, e.g., [[18]]. The overland flow is barely considered as a fully dynamic flow, and is normally solved for by hydrologic concepts such as linear reservoir [[19]] or is assumed to be a two-dimensional (2D) kinematic/diffusive wave equation [[16]]. Hence, the flood characteristics (i.e., flood extent, water depth, and velocity) cannot be represented inside the two-dimensional domain.
Currently, the two-dimensional hydraulic model provides a good solution for flood routing processes for fine flood risk assessment and management [[20]]. The fully 2D hydraulic model is a robust numerical methodology to provide accurate simulation of runoff flow on a minutely temporal scale, based on high-resolution data sources and GPU-accelerated techniques [[23]]. However, runoff generation is still a big concern for the application of the fully 2D hydraulic model in the large-scale of a basin, especially for small and medium floods. The typical solution of infiltration is to adopt the Kostiakov equation (e.g., [[26]]), Green–Ampt model (e.g., [[27]]), Horton equation (e.g., [[28]]), or Richard model (e.g., [[30]]).
In recent years, how to jointly use the hydrologic and hydraulic models for flood risk assessment and reservoir operation has been discussed [[31]]. As has been shown previously, hydrologic models are good at simulating rainfall-runoff processes for a basin on a large scale, e.g., several hundred km2, while hydraulic models are more advantageous for applications on smaller scales. In order to take advantage of these two types of models, this study attempts to couple a hydrologic model with a hydraulic model.
In this study, the Misai river basin (797 km2), located in Zhejiang Province, China, is taken as the study area. For the hydrologic model, XAJ was selected as it has been widely used in China for a long time [[34]]. For the hydraulic model, the high-performance integrated hydraulic modelling system (HiPIMS) was chosen because of its high-speed calculations and wide application [[23]]. This paper introduces how the XAJ–HiPIMS is coupled from XAJ and HiPIMS and validates the new model against typical floods. Section 2 outlines the materials and methods, Section 3 presents results and discussion, and Section 4 gives conclusions.
2. Materials and Methods
2.1. XAJ
XAJ is a lumped rainfall-runoff model developed by the authors of [[11]]. Based on the watershed saturation-excess runoff theory, evapotranspiration is calculated in three soil layers. The soil's moisture storage capacity distribution curve is utilized to provide a non-uniform distribution of storage capacity in the hydrologic unit. The total runoff is composed of surface runoff, interflow, and groundwater runoff. The distribution of these three components can be decided by using a free water capacity distribution curve. The surface runoff is routed by the unit hydrograph. The interflow and groundwater flow are routed by the linear reservoir method. Hence, four modules compose the XAJ: evapotranspiration, runoff generation, runoff sources partition, and runoff routing ([[11]], Figure 1). The parameters of XAJ are described in Table 3. The output of XAJ is outflow discharge at the outlet section of the basin.
2.2. HiPIMS
The governing equations of HiPIMS are fully 2D shallow water equations based on the Cartesian uniform grid [[38]], as follows,
(1)
(2)
In Equation (1), q stands for the vector of conserved flow variables; f and g are the flux vectors in the x- and y-directions, respectively; and R, Sb, and Sf denote the source terms of rainfall, bed slope, and friction, respectively. In Equation (2), h is the water depth (h = η − zb, where η is the water surface elevation and zb is the bed elevation); u and v are the depth-averaged velocity components in the x- and y-directions, respectively; g is the acceleration of gravity; ρ is the water density; r is the generated surface runoff (which is the in Equation (2)); ∂b / ∂x and ∂b / ∂y denote the bed slope in the x- and y-directions; τbx and τby are the bed friction stresses, calculated by the following equations:
(3)
where Cf (= gn2 / h1/3) stands for the bed roughness coefficient; and n is the Manning's coefficient.
In this study, HiPIMS employs a first-order finite volume Godunov-type numerical scheme to solve the above governing equations. Equation (1) is discretized by Equation (4) and the flow variables can be updated as follows.
(4)
where the superscript m denotes the present time level; Δt is the present time step (which is the in Equation (5)); Ωi stands for the area of cell i; Fk(q) presents the fluxes that are normal to the cell edges; k is the index of the cell edges (k = 1–4); lk is the corresponding cell edge length; and n = (nx,ny) denotes the unit vector of the outward normal direction. In HiPIMS, the flux terms and bed slope term are calculated by an explicit scheme. The interface fluxes are calculated by the HLLC (Harten–Lax–van Leer-Contact) approximate Riemann solver. The local Riemann problems at the cell interfaces are solved by implementing the surface reconstruction method [[39]]. The friction term is solved by an implicit scheme [[40]], which is capable of maintaining the numerical stability when dealing with very small water depth. The adaptive time step method is implemented and controlled by the Courant–Friedrichs–Lewy (CFL) criterion. The high-performance of this CPU/GPU-integrating model is achieved by adopting the OpenCL programming framework [[38], [41]].
2.3. Coupling Framework
In the coupling framework, XAJ generates the surface runoff and produces interflow and underground flow while the hydraulic model is run for 2D confluence calculations of the surface runoff generated by XAJ. The generated surface runoff was presented in the governing equations of HiPIMS as the source term in the mass conservation equation. The temporal scale of XAJ is hourly while the computing time-step of HiPIMS is in seconds, and the spatial scale of XAJ is the whole basin, while computing grid-cell of HiPIMS is in meters. In order to deal with this temporal and spatial downscaling issue, functions (Equations (5) and (6)) were created as below.
(5)
(6)
In Equation (5), RXAJ denotes the surface runoff calculated by the XAJ; is the time-step of HiPIMS; n means the presently considered time-step; and stands for the downscaled surface runoff. In Equation (6), means the surface runoff in terms of the source term in the mass conservation equation of HiPIMS and m is the amount of the computing grid cell.
Then, the surface runoff generated by XAJ inputs to the hydraulic model to produce two-dimensional flow over the basin. Thus, the procedure is proposed as follows (shown in Figure 2): (1) Calculate the surface runoff by XAJ with the temporal scale of 1 h; (2) Allocate the surface runoff using Equation (5) in order to transform the temporal scale from hours to seconds; (3) Distribute the surface runoff into each computing grid cell of HiPIMS using Equation (6); (4) Compute the dynamic surface flow process and produce the spatial distribution of the inundated extent, water depth, and velocity in the whole basin by HiPIMS; (5) Overlay the surface flow process at the basin outlet produced by HiPIMS and the processes of interflow and underground flow given by XAJ to form the simulated flow process at the basin outlet and validate it against the flow observed at the outlet.
2.4. Statistical Method
In this work, the statistical methods absolute relative error of flood-peak discharge (ARED) (Equation (7)), difference of peak arrival time (DPAT) (Equation (8)), and the Nash–Sutcliffe efficiency coefficient (NSE) (Equations (9) and (10)) are employed to estimate the discrepancy of the discharge. Three types of data are considered herein, i.e., observation, simulation by XAJ, and simulation by XAJ–HiPIMS.
(7)
where MAX(Qs) stands for the simulated peak discharge and MAX(Qo) denotes the observed peak discharge. ARED∈ (−∞,+∞),
(8)
where TMAX(Qs) is the arrival time of simulated peak discharge by XAJ and XAJ–HiPIMS and TMAX(Qo) is the arrival time of the observed peak discharge.
(9)
(10)
where NSE1 is the Nash–Sutcliffe efficiency coefficient computed by the XAJ simulation and observation; NSE2 is the Nash–Sutcliffe efficiency coefficient computed by the XAJ–HiPIMS simulation and observation; Qo is the observed discharge; Q1 is the simulated discharge of XAJ; Q2 is the simulated discharge of XAJ-HiPIMS; and is the averaged observed discharge.
2.5. Study Area
The Misai basin is located in the Zhejiang province of China with an area of 797 km². The bed elevation differs from 105 to 1253 m. In this study, the 30-m ASTER GDEM V2 elevation data were adopted for computing, as shown in Figure 3. Six precipitation gauging stations are located in this basin. One is located at the outlet of the basin, and it also served as stream gauging station. The 30-m Landsat TM image data were used to provide the land use information. Considering the land characteristics, the land use types were classified into seven categories in this work, i.e., forest, heavy brush, cultivated land, grass land, pond/river, bare land, and urban land (distribution shown in Figure 4). The proportional area and Manning's coefficients of land use types are summarized in Table 1.
2.6. Flood Processes
Three typical flood processes (FPs) were chosen to represent big (B: peak discharge > 1000 m3/s), medium (M: 500 m3/s < peak discharge < 1000 m3/s), and small (S: peak discharge < 500 m3/s) floods, as shown in Table 2.
2.7. Modelling Set
XAJ: The rainfall input was given at 1-h intervals from the historic records from the six gauging stations in the basin. The 13 coefficients were calibrated and are shown in Table 3.
HiPIMS: The considered basin was divided into 909,824 grid cells upon a regular uniform grid of resolution 30 m. The Manning's n are given in Table 1. The CFL number was set to 0.35 in the present work.
3. Results and Discussion
In this study, the simulations of XAJ–HiPIMS and XAJ were assessed by comparing their results with the observed discharge at the outlet gauging station, as presented in Figure 5. Good agreement of the simulation with the observation was demonstrated for the three typical floods. Better performance was found in the bigger flood category.
Further assessment has been done by implementing the statistical methods of ARED, DPAT, and NSE, shown in Table 4. We found that (1) in terms of the peak discharge, the simulation of XAJ–HiPIMS presented good agreement with the measurement; (2) in terms of ARED, the XAJ–HiPIMS achieved a better flood peak simulation, especially for the bigger flood; (3) in terms of DPAT, the differences between simulations and measurements were no more than 1 h and smaller differences were found in the medium flood; (4) in terms of NSE, XAJ–HiPIMS obtained good flood simulations, which were greater than 0.85. The above results demonstrate XAJ–HiPIMS could provide satisfactory simulation at different flood scales.
According to the 2D distribution of flood information, inundated extent and maximum water depth were investigated. The threshold of inundated depth was set to 0.01 m in this study. The snapshots of inundated extent were produced by XAJ–HiPIMS at the peak moment (the moment of peak discharge) shown in Figure 6 (left column). The value of water depth in each grid cell was recorded during the calculation. Maximum water depth is defined as the maximum value of water depth in each grid cell in the entire flood event record. The maps of maximum water depth were analyzed and are presented in Figure 6 (right column).
We found that: (1) Inundated area was mostly located along the river network and in the low-lying urban areas and cultivated land. Inundated extent and maximum water depth were obviously increased due to the bigger rainfall height. (2) In three FPs, maximum water depth in low-lying areas of urban land exceed 1 meter. In FP(B), cultivated land was also inundated with the maximum water depth of more than 1 meter. The results suggest that spatial distribution of flood information could be adopted in flood risk assessment and also help with the planning of flooding mitigation and adaptation strategies.
These findings imply that the XAJ–HiPIMS (1) presents good agreement with the observations; (2) can provide 2D presentation of flood processes and more information of flood characteristics at the basin scale; (3) provides a new perspective and methodology for ecological and environmental protection, flood prevention and management, and urban water landscape design.
4. Conclusions
This study presented a coupled hydrologic–hydraulic model, XAJ–HiPIMS. The model was validated in the Misai basin, Zhejiang Province, China by comparing its performance with measurements and statistically assessing the absolute relative error of flood-peak discharge, the difference of peak arrival time, and the Nash–Sutcliffe efficiency coefficient.
The present study: (1) demonstrates that the XAJ–HiPIMS can provide a good performance in basin-scale flood simulation and deliver a 2D viewpoint of flood characteristics; (2) reveals the potential application of XAJ–HiPIMS in ecological and environmental protection, urban flood management, landscape design, and climate change analysis; and (3) brings forward the necessity and recommendations for further research.
Figures and Tables
Graph: Figure 1 The framework of the Xinanjing model (XAJ) [[37]].
Graph: Figure 2 Flowchart of the coupling framework.
MAP: Figure 3 Map of the Misai basin.
Graph: Figure 4 Land use of the Misai basin.
Graph: Figure 5 Time history of discharge at the outlet gauging station.
MAP: Figure 6 (a) FP(S): Snapshot of inundated extent at the peak moment (left) and map of maximum water depth (right); (b) FP(M): Snapshot of inundated extent at the peak moment (left) and map of maximum water depth (right); (c) FP(B): Snapshot of inundated extent at the peak moment (left) and map of maximum water depth (right).
Graph: water-12-01288-g006b.tif
Table 1 The proportional area and Manning's n of different land use types in the Misai basin.
No. | Land Use Type | Area Proportion (%) | Manning's n1 |
---|
L1 | Forest | 45.5 | 0.139 |
L2 | Heavy brush | 37.3 | 0.098 |
L3 | Cultivated land | 2.2 | 0.041 |
L4 | Grass land | 7.3 | 0.031 |
L5 | Pond/river | 5.8 | 0.021 |
L6 | Bare land | 0.7 | 0.026 |
L7 | Urban land | 1.2 | 0.013 |
1 summarized from [[42]].
Table 2 The three typical flood processes (FPs).
Name | Start Time | End Time | Rainfall Height | Peak Flow |
---|
FP(S) | 1987/5/26 8:00 | 1987/5/27 20:00 | 65 mm | 202 m3/s |
FP(M) | 1985/5/5 8:00 | 1985/5/7 12:00 | 106 mm | 708 m3/s |
FP(B) | 1983/5/29 8:00 | 1983/5/30 13:00 | 222 mm | 1820 m3/s |
Table 3 The parameters of XAJ.
Module | Parameters | Physical Meaning | Value |
---|
Evapotranspiration | WUM | Averaged soil moisture storage capacity of the upper layer | 14 |
| WLM | Averaged soil moisture storage capacity of the lower layer | 86 |
| WDM | Averaged soil moisture storage capacity of the deep layer | 33 |
| K | Conversion coefficient of evaporation | 1 |
| C | Coefficient of the deep layer | 0.126 |
Runoff generation | B | Exponential of the distribution to tension water capacity | 0.375 |
| IMP | Percentage of impervious and saturated areas in the catchment | 10 |
Runoff source partition | SM | Areal mean free water capacity of the surface soil layer | 97 |
| EX | Exponent of the free water capacity curve influencing the development of the saturated area | 1.03 |
| KG | Outflow coefficients of the free water storage to groundwater relationships | 0.459 |
| KSS | Outflow coefficients of the free water storage to interflow relationships | 0.07 |
Runoff routing | KKG | Recession constants of the groundwater storage | 0.997 |
| KKSS | Recession constants of the lower interflow storage | 0.747 |
Table 4 The statistical analysis of the measured and simulated discharge (Q).
FP No. | Peak discharge (m3/s) | ARED (%) | DPAT (hour) | NSE |
---|
O | H | C | H | C | H | C | NSE1 | NSE2 |
---|
FP(S) | 202 | 291 | 250 | 0.44 | 0.24 | 0 | 1 | 0.6392 | 0.8459 |
FP(M) | 708 | 900 | 822 | 0.27 | 0.16 | 0 | 0 | 0.7111 | 0.8524 |
FP(B) | 1820 | 1620 | 1817 | 0.11 | 0.002 | −1 | −1 | 0.9757 | 0.9422 |
O: Observed data; H: Simulated by XAJ; C: Simulated by XAJ–HiPIMS.
Author Contributions
All authors have read and agreed to the published version of the manuscript. Conceptualization, Y.W. and X.Y.; methodology, Y.W. and X.Y.; software, Y.W. and X.Y.; validation, Y.W.; formal analysis, Y.W.; investigation, Y.W. and X.Y.; resources, Y.W. and X.Y.; data curation, Y.W.; writing—original draft preparation, Y.W.; writing—review and editing, Y.W. and X.Y.; visualization, Y.W.; supervision, X.Y.; project administration, X.Y.; funding acquisition, Y.W.
Funding
This research was funded by the National Natural Science Foundation of China, grant numbers 41890823 and 51609231.
Conflicts of Interest
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Acknowledgments
The authors would like to acknowledge Liang Qiuhua of Loughborough University of UK for kindly providing the HiPIMS software used in this research, and Huili Chen of Loughborough University of UK for her interpretation of land use information, and Chaowei Xu of Peking University for his helping in XAJ use. Special thanks go to the three anonymous reviewers for their comments and suggestions, valuable and helpful in improving the quality of this paper.
[
References
1
Tsakiris G. Flood risk assessment: concepts, modelling, applications. Nat. Hazards Earth Syst. Sci. 2014; 14: 1361-1369. 10.5194/nhess-14-1361-2014
2
IPCC, IPCC Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate ChangeStocker T.F., Qin D.P., Allen S.K., Boschung J., Nauels A., Xia Y., Bex V., Midgley P.M. ; Cambridge University Press: Cambridge, UK: New York, NY, USA. 1535: 23
3
Kundzewicz Z.W., Kanae S., Seneviratne S.I., Handmer J., Nicholls N., Peduzzi P., Mechler R., Bouwer L.M., Arnell N., Mach K. Flood risk and climate change: global and regional perspectives. Hydrol. Sci. J. 2014; 59: 1-28. 10.1080/02626667.2013.857411
4
Hirabayashi Y., Mahendran R., Koirala S., Konoshima L., Yamazaki D., Watanabe S., Kim H., Kanae S. Global flood risk under climate change. Nat. Clim. Chang. 2013; 3: 816-821. 10.1038/nclimate1911
5
Wang Y., Yang X. Assessing the effect of land use/cover change on flood in Beijing, China. Front. Environ. Sci. Eng. 2013; 7: 769-776. 10.1007/s11783-013-0542-z
6
Yang X., Chen H., Wang Y., Xu C. Evaluation of the effect of land use/cover change on flood characteristics using an integrated approach coupling land and flood analysis. Hydrol. Res. 2016; 47: 1161-1171. 10.2166/nh.2016.108
7
Courty L.G., Rico-Ramirez M.A., Pedrozo-Acuña A. The Significance of the Spatial Variability of Rainfall on the Numerical Simulation of Urban Floods. Water. 2018; 10207. 10.3390/w10020207
8
Erena S.H., Worku H. Flood risk analysis: Causes and landscape based mitigation strategies in Dire Dawa city, Ethiopia. Geoenviron. Disasters. 2018; 5: 16. 10.1186/s40677-018-0110-8
9
Arnold J.G., Srinivasan R., Muttiah R.S., Williams J.R. Large area hydrologic modelling and assessment Part I: model development. J. Am. Water Resour. Assoc. 1998; 34: 73-89. 10.1111/j.1752-1688.1998.tb05961.x
Neitsch S.L., Arnold J.G., Williams J.R., Kiniry J.R., King K.W. Soil and Water Assessment Tool Theoretical Documentation: Version 2000. 2002Available online: https://swat.tamu.edu/media/1290/swat2000theory.pdf(accessed on 20 April 2020)
Zhao R.J., Zhang Y.L., Fang L.R., Liu X.R., Zhang Q.S. The Xinanjiang model. Proceedings of the Hydrological Forecasting Symposium. Wallingford, UK. 15–18 April 1980; IAHS: Wallingford, UK: 351-356
Zhao R., Liu X. The Xinanjiang model. Computer Models of Watershed HydrologySingh V.P. ; Water Resource Publications, LLC: Highlands Ranch, CO, USA. 1995
Crawford N.H., Linsley R.K.. Digital simulation in hydrology: Stanford Watershed Model IV; Stanford University: Stanford, CA, USA. 1966
Finnerty B.D., Smith M.B., Seo D.-J., Koren V., Moglen G.E. Space-time scale sensitivity of the Sacramento model to radar-gage precipitation inputs. J. Hydrol. 1997; 203: 21-38. 10.1016/S0022-1694(97)00083-8
Meng C., Zhou J., Zhong D., Wang C., Guo J. An Improved Grid-Xinanjiang Model and Its Application in the Jinshajiang Basin, China. Water. 2018; 101265. 10.3390/w10091265
Yao C., Li Z., Yu Z., Zhang K. A priori parameter estimates for a distributed, grid-based Xinanjiang model using geographically based information. J. Hydrol. 2012; 468: 47-62. 10.1016/j.jhydrol.2012.08.025
Refsgaard J.C., Storm B. MIKE SHE. Computer Models of Watershed HydrologySingh V.P. ; Water Resource Publications, LLC: Highlands Ranch, CO, USA. 1995
Waseem M., Kachholz F., Klehr W., Tränckner J. Suitability of a Coupled Hydrologic and Hydraulic Model to Simulate Surface Water and Groundwater Hydrology in a Typical North-Eastern Germany Lowland Catchment. Appl. Sci. 2020; 101281. 10.3390/app10041281
Xiong L., Guo S., Tian X. DEM-based distributed hydrological model and its application. Adv. Water Sci. 2004; 15: 517-520(In Chinese)
Wang Y., Liang Q., Kesserwani G., Hall J.W. A 2D shallow flow model for practical dam-break simulations. J. Hydraul. Res. 2011; 49: 307-316. 10.1080/00221686.2011.566248
Liang Q., Xia X., Hou J. Efficient urban flood simulation using a GPU-accelerated SPH model. Environ. Earth Sci. 2015; 74: 7285-7294. 10.1007/s12665-015-4753-4
Wang C., Hou J., Miller D., Brown I., Jiang Y. Flood risk management in sponge cities: The role of integrated simulation and 3D visualization. Int. J. Disaster Risk Reduct. 2019; 39: 101139. 10.1016/j.ijdrr.2019.101139
Smith L.S., Liang Q. Towards a generalised GPU/CPU shallow-flow modelling tool. Comput. Fluids. 2013; 88: 334-343. 10.1016/j.compfluid.2013.09.018
Liang Q., Smith L.S. A high-performance integrated hydrodynamic modelling system for urban flood simulations. J. Hydroinform. 2015; 17: 518-533. 10.2166/hydro.2015.029
Xia X., Liang Q., Ming X. A full-scale fluvial flood modelling framework based on a high-performance integrated hydrodynamic modelling system (HiPIMS). Adv. Water Resour. 2019; 132: 103392. 10.1016/j.advwatres.2019.103392
Bellos V., Tsakiris G. A hybrid method for flood simulation in small catchments combining hydrodynamic and hydrological techniques. J. Hydrol. 2016; 540: 331-339. 10.1016/j.jhydrol.2016.06.040
Yin J., Yu D., Yin Z., Liu M., He Q. Evaluating the impact and risk of pluvial flash flood on intra-urban road network: A case study in the city center of Shanghai, China. J. Hydrol. 2016; 537: 138-145. 10.1016/j.jhydrol.2016.03.037
Horton R.E. The Rôle of infiltration in the hydrologic cycle. Trans. Am. Geophys. Union. 1933; 14: 446. 10.1029/TR014i001p00446
Horton R.E. The Interpretation and Application of Runoff Plat Experiments with Reference to Soil Erosion Problems. Soil Sci. Soc. Am. J. 1939; 3: 340-349. 10.2136/sssaj1939.036159950003000C0066x
Mudd S. Investigation of the hydrodynamics of flash floods in ephemeral channels: Scaling analysis and simulation using a shock-capturing flow model incorporating the effects of transmission losses. J. Hydrol. 2006; 324: 65-79. 10.1016/j.jhydrol.2005.09.012
Anselmo V., Galeati G., Palmieri S., Rossi U., Todini E. Flood risk assessment using an integrated hydrological and hydraulic modelling approach: A case study. J. Hydrol. 1996; 175: 533-554. 10.1016/S0022-1694(96)80023-0
Seibert S., Skublics D., Ehret U. The potential of coordinated reservoir operation for flood mitigation in large basins—A case study on the Bavarian Danube using coupled hydrological—hydrodynamic models. J. Hydrol. 2014; 517: 1128-1144. 10.1016/j.jhydrol.2014.06.048
Li W., Lin K., Zhao T., Lan T., Chen X.-H., Du H., Chen H. Risk assessment and sensitivity analysis of flash floods in ungauged basins using coupled hydrologic and hydrodynamic models. J. Hydrol. 2019; 572: 108-120. 10.1016/j.jhydrol.2019.03.002
Ren-Jun Z. The Xinanjiang model applied in China. J. Hydrol. 1992; 135: 371-381. 10.1016/0022-1694(92)90096-E
Jayawardena A., Zhou M. A modified spatial soil moisture storage capacity distribution curve for the Xinanjiang model. J. Hydrol. 2000; 227: 93-113. 10.1016/S0022-1694(99)00173-0
Zhang D., Zhang L., Guan Y., Chen X., Chen X. Sensitivity analysis of Xinanjiang rainfall–runoff model parameters: A case study in Lianghui, Zhejiang province, China. Hydrol. Res. 2012; 43: 123-134. 10.2166/nh.2011.131
Zhao R.J.. Water Hydrological Modelling—Xinanjiang Model and Shanbei Model; China Water Resources and Hydropower Publishing House: Beijing, China. 1984: 125-129(in Chinese)
Liang Q., Smith L., Xia X. New prospects for computational hydraulics by leveraging high-performance heterogeneous computing techniques. J. Hydrodyn. 2016; 28: 977-985. 10.1016/S1001-6058(16)60699-6
Xia X., Liang Q., Ming X., Hou J. An efficient and stable hydrodynamic model with novel source term discretization schemes for overland flow and flood simulations. Water Resour. Res. 2017; 53: 3730-3759. 10.1002/2016WR020055
Jin S. Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: A review. Rivista di Matematica della Universita di Parma. 2010; 2: 177-216
Liang Q., Xia X., Hou J. Catchment-scale High-resolution Flash Flood Simulation Using the GPU-based Technology. Procedia Eng. 2016; 154: 975-981. 10.1016/j.proeng.2016.07.585
Zheng Z., He S., Wu F., Hu J. Relationship between surface roughness and Manning roughness. J. Mt. Sci. 2004; 22: 236-239(in Chinese). 10.16089/j.cnki.1008-2786.2004.02.018
Tang X. Characteristics of runoff and soil erosion of typical forest types in the sources of Qiantang River. Master Thesis; Zhejiang A&F University: Hangzhou, China22May. 2012(in Chinese)
Huang J. Evaluation of anti-water erosion function of main forests in the south hilly and mountain region of Jiangsu Province. PhD Thesis; Nanjing Forestry University: Nanjing, China16June. 2011(in Chinese)
]
By Yueling Wang and Xiaoliu Yang
Reported by Author; Author