A Method for Estimating Urban Flood-Carrying Capacity Using the VIS-W Underlying Surface Model: A Case Study from Wuhan, China

: In 2016, 192 cities, including Wuhan, Nanjing, and Jiujiang, su ﬀ ered from severe ﬂooding, which raised social and government concerns in China. This paper proposes a method based on the underlying surface to estimate the urban ﬂood-carrying capacity. First, water is extended into the underlying surface vegetation-impervious surface-soil (VIS) model to form the vegetation-impervious surface-soil and water (VIS-W) model. Second, the watershed is delimitated using a digital elevation model (DEM). Third, the natural water system’s storage capacity, including rivers, lakes, and puddles, is estimated. Using the rainfall–runo ﬀ hydrology model and the VIS-W model, the storing ability, the receiving ability, the discharge potential, and the emergency of pumping outward are assessed to derive the ﬁnal ﬂood-carrying capacity. Finally, the result is compared with data on the waterlogging points collected in July 2016 during a ﬂood and waterlogging event in Wuhan. It is found that 84% of waterlogging points are located in the “weak” or “normal” areas, and 16% are located in the “strong” areas. Additionally, 99% of total waterlogging points are located upstream, based on the stream extracted by the DEM. This phenomenon indicates that Wuhan can mitigate ﬂood disasters by fully utilizing the natural water system storage capacity to corporate the city’s drainage schedule. from itself. The low-lying area F low is the area with a slope ratio of less than 1%, except for the area for which the land cover type is water. The multiplied results of Len ( F Outlet ) , Count ( F Areas ) , and F Low , which are recorded as f ( ∗ ) , will be used to deﬁne the level linearly, as “not urgent”, “normal”, and “urgent”.


Introduction
In recent years, China has faced serious challenges in terms of flood and waterlogging disasters. According to the Global Disaster Annual Report, the World Health Organization (WHO) pointed out that in 2016 China experienced its fifth highest number of natural disasters in the last decade, and was affected by 16 floods and landslides in that year [1]. Additionally, the national report released by the Ministry of Civil Affairs of the People's Republic of China (MCA) reported that in 2016, 192 cities, including Wuhan, Nanjing, and Jiujiang, suffered from serious rainstorms, which directly led to a disaster-affected population of 3.89 million and economic losses of RMB 9.1 billion [2]. The nationwide flood and waterlogging disasters caught the attention of researchers in the field of flood control and management [3]. Extreme rainstorms, the increasing ratio of impervious surface [4,5], and urbanization generated changes to transform original suburbs to semi-urban or urban land areas [6], which might affect the local original hydrological conditions [7][8][9]. Taking Wuhan as an example, the imperious surface increased from 903.45 km 2 to 3989.49 km 2 from 2010 to 2016 [10]. This phenomenon caused affecting the water route path could be used for reference in the estimation of the urban floodcarrying capacity studies.
In the urban catchment, the waterlogging points located in the center of the city will have a more serious effect than in the suburban or rural land. This indicates that the flood-carrying capacity needs to consider the potential route of the rainwater, the waterlogging distribution, and the residents' social activities. Therefore, in terms of assessing the degree of safety over a certain object-carrying capacity, the urban flood-carrying capacity should be defined as "the distribution capacity to store or trap rainwater in the urbanized watershed during the flood or waterlogging events under the condition that the city's normal routines will not be disturbed". It is obvious that the balance between budget of rainwater and storage capacity will affect the distribution of the surface overflow. The main purpose of this study is to estimate urban flood-carrying capacity in terms of the terrain's hydrological features, based on the underlying surface features and the hydrological characteristics of the terrain. In this paper, we attempt to derive the flood-carrying capacity based on the storage of the natural water system (e.g., rivers, lakes, ponds, and other entities), relayed by the digital elevation model (DEM) data set and the latest underlying surface information. This might provide timely information about the local rainwater-carrying limitations during daily flood control and management works.

Study Area
The study was conducted in the downtown area of Wuhan. Wuhan is the capital city of Hubei Province, and it is a metropolis located in central China. Wuhan is located at 113°41' E to 115°05' E and 29°58' N to 31°22' N, and it has a humid subtropical monsoon climate with abundant rainfall. The annual precipitation is 1150-1450 mm, and the rainfall is concentrated from June to August each year, so precipitation occurs for almost 40% of the whole year [36]. There are 13 districts, among which the Jiang'an, Jianghan, Qiaokou, Hanyang, Wuchang, Hongshan, and Qingshan districts are downtown areas, and the Dongxihu, Caidian, Jiangxia, Huangpi, Xinzhou, and Hannan districts are new towns. We chose the seven downtown districts to represent the study area in this research as shown in Figure 1. The world's third longest river, the Yangtze River, and its tributary, the Hangjiang River, meet in the center of Wuhan. They divide the city into three parts: Wuchang, Hankou, and Hanyang. Wuhan is a typical city facing severe flood and waterlogging threats due to the fact that the terrain is The world's third longest river, the Yangtze River, and its tributary, the Hangjiang River, meet in the center of Wuhan. They divide the city into three parts: Wuchang, Hankou, and Hanyang. Wuhan is a typical city facing severe flood and waterlogging threats due to the fact that the terrain is dominated by low-lying hills and almost plain areas, and most areas in the city are between 20 m and 24 m in elevation. Additionally, the ordinary water level is generally 23 meters, which is higher than the elevation of most urban areas, and local rainstorms can only be stored in lakes or pumped into the Yangtze River [14]. In flood seasons, the flood prevention level of the Yangtze River is 25 meters, but the water level may reach 27-28 meters during heavy rainstorms. For example, the peak water level reached 29.73 m and 29.43 m in 1954 and 1998, respectively. Meanwhile, there are 166 lakes in Wuhan, of which 43 are located in the center of the city. Therefore, realizing the flood storage potential of numerous lakes, and identify the storage capacity of lakes in the natural water system as a whole, might be effective measures for relieving flood and waterlogging disasters in Wuhan.

Data Collection
For the purposes of estimating the carrying capacity for storing and trapping rainwater, underlying surface characteristics such as the terrain, land cover, natural water system, and potential drainage ability are very important. We collected various data sets, including those on districts, the digital elevation model (DEM), the land use type, meteorological data, and the waterlogging points, as follows: • District dataset: We adopted Hubei province administrative division data (2016 edition) to obtain the water system information and city's district boundary of Wuhan.

Assessment Method
We propose an integrated estimation method in terms of DEM and the underlying surface model to obtain the urban flood-carrying capacity with respect to four basic abilities: rainwater storage ability, rainwater receiving ability, self-discharging potential, and the emergency of pumping outward. The main assessment process is shown in Figure 2. for the impact in the flow route path, underlying surface characteristics such as slope, elevation, land used type, and stream distribution are widely used in the flood risk qualitative analysis [37,38]. In this paper, we chose 11 factors to estimate the flood-carrying related abilities. There are three basic maps: land-use type, watershed, and slope, which are listed in Table 1. The remaining eight are logical indicators which are related to the storage capacity ( ), land-use type (Ai), the terrain feature ( ), and the stream link following the flow path ( , , , ( ), and ( )), as listed in Table 2 Fourth, we estimated the urban flood-carrying capacity. The four basic abilities related to the urban flood-carrying capacity are the rainwater-storing ability, the rainwater-receiving ability, the self-discharging potential, and the emergency of pumping outward, respectively. Based on the four basic abilities, we can derive the city's flood-carrying capacity.
(1) Estimate the four basic abilities A. Rainwater-storing ability The rainwater storing ability is determined by the precipitation interval according to the rainfall intensity over a 6-hour duration. This ability is used to express the rainwater capacity of the current sub-watershed. Compared with rainfall intensity of a 12-hour duration, the rainwater that converges in a 6-hour duration will have a better chance of being collected from runoff rather than from the infiltrate or other sources. It is labeled as "very weak", "weak", "normal", "strong", and "very strong". If the ( ) value is "<1", the ability will be defined as "very weak". The states of "weak", "normal" "strong", and "very strong" correspond to the return period intervals of "1~5", "5~20", "20~100", and "~100" or ">100", respectively. The urban flood-carrying capacity assessment process is made up of four steps. First, we extended the underlying surface model; second, we identified the parameters; third, we constructed the basic estimating maps and indicators; and finally, we estimated the urban flood-carrying capacity.
First, we extended the underlying vegetation-imperious surface-soil (VIS) surface model. We used the water land type from the traditional underlying surface model VIS, resulting in the VIS-W (vegetation-imperious surface-soil and water) model. Second, we pre-defined parameters related to the study area. In this paper, the runoff coefficient and the design precipitation of the city needed to firstly be defined. The former is used to calculate the competitive runoff coefficient of each sub-watershed, based on the VIS-W underlying surface model. The precipitation capacity corresponding to the storage capacity of each sub-watershed can be calculated using the hydrological rainfall-runoff model. The design precipitation can be used to identify the rainwater-storing capacity by the interval precipitation of a certain hydrological frequency.
Third, we constructed the basic estimating maps and indicators. In flood hazard analysis studies, for the impact in the flow route path, underlying surface characteristics such as slope, elevation, land used type, and stream distribution are widely used in the flood risk qualitative analysis [37,38]. In this paper, we chose 11 factors to estimate the flood-carrying related abilities. There are three basic maps: land-use type, watershed, and slope, which are listed in Table 1. The remaining eight are logical indicators which are related to the storage capacity (V Capacity ), land-use type (A i ), the terrain feature (F Low ), and the stream link following the flow path (L River , F Terminaltype , F Nexttype , Len(F Outlet ), and Count(F Areas )), as listed in Table 2. Storage capacity including river, lake, and puddle or pond / "<1", "1-5", "5-20", "20-100", "~100" and ">100" Fourth, we estimated the urban flood-carrying capacity. The four basic abilities related to the urban flood-carrying capacity are the rainwater-storing ability, the rainwater-receiving ability, the self-discharging potential, and the emergency of pumping outward, respectively. Based on the four basic abilities, we can derive the city's flood-carrying capacity.
(1) Estimate the four basic abilities A. Rainwater-storing ability The rainwater storing ability is determined by the precipitation interval according to the rainfall intensity over a 6-hour duration. This ability is used to express the rainwater capacity of the current sub-watershed. Compared with rainfall intensity of a 12-hour duration, the rainwater that converges in a 6-hour duration will have a better chance of being collected from runoff rather than from the infiltrate or other sources. It is labeled as "very weak", "weak", "normal", "strong", and "very strong". If the R Int(6h) value is "<1", the ability will be defined as "very weak". The states of "weak", "normal" "strong", and "very strong" correspond to the return period intervals of "1~5", "5~20", "20~100", and "~100" or ">100", respectively.
weak, R Int(6h) = "1-5" normal, R Int(6h) = "5-20" strong, R Int(6h) = "20-100" very strong, R Int(6h) = " ∼ 100" OR " > 100" (1) B. Rainwater-receiving ability The rainwater-receiving ability is estimated by the return period interval of 12 h in duration and its stream level. Considering the distance flow route and the effect of infiltrating, we adopted a duration of 12 hours to estimate the rainwater receiving ability. First, corresponding to the precipitation of 12-h rainfall, the interval return periods of "<1", "1-5", "5-20", "20-100", and "~100" or ">100" are labeled as "very weak", "weak", "normal", "strong", and "very strong" respectively. Then, the level will be changed according to the stream link level. If the stream level is 1, the level state is unchanged. If the stream level is 2, the ability state will be moved down 1 grade, e.g., the level of "strong" will be changed as "normal". Stream level 3 will be moved down two grades in terms of the original level, until its state is "very weak".
(R Int(12h) = "1-5" AND L River = 1) OR (R Int(12h) = "5-20" AND L River = 2) OR (R Int(12h) = "20-100"AND L River = 3) OR (R Int(12h) IN [" ∼ 100", " > 100"]AND L River = 4) normal, (R Int(12h) = "5-20"AND L River = 1) OR (R Int(12h) = "20-100" AND L River = 2) OR ((R Int(12h) IN [" ∼ 100", " > 100"]) AND L River = 3) strong, (R Int(12h) = "20-100"AND L River = 1) OR ((R Int(12h) IN [" ∼ 10", " > 100"]) AND L River = 2) very strong, (R Int(12h) IN [" ∼ 100", " > 100"]) AND L River = 1 (2) C. Self-discharging potential The self-discharging potential is determined by the rainwater flow path and the outlet. The outlet refers to the outlet of the whole area, including the entrance to the sea, or the lowest downstream outlet along the flow path. In this paper, the Yangtze River is the main river inside Wuhan, and its branch is the Han River. Therefore, we defined the range of river both the Yangtze River and Han River as the main outlet. The sub-watershed that contains or is adjacent to the main outlet will be labeled as "very strong". For the sub-watershed for which the flow path is out away from the study area, it will be considered as in the opposite direction of the main outlet, and its level will be "very weak". The other three levels, "strong", "normal", and "weak", will be determined in term of the terminal flow type and the neighbor watershed type along the flow path.
D. The pumping outward emergency The pumping outward emergency expresses the balance between local storing and collecting ability in low-lying parts of each sub-catchment. We used the amounts of sub-watersheds near the outlet (len(F river )), amounts of sub-watersheds upstream (Count(F areas )), and the ratio of low-lying areas (F low ) to estimate the emergency for pumping outward, rather than adopting the numerical hydrology model to simulate the budget of water storage capacity. The value of len(F river ) will be valued as 0 if it contains the main outlet. Meanwhile, if its terminal flow type is "outer", its value will be adopted twice of the maximum len(F river ) value. The default value of Count(F areas ) is set as 1, which means current sub-watershed is located at the head of stream, so that it only receives rainwater from itself. The low-lying area F low is the area with a slope ratio of less than 1%, except for the area for which the land cover type is water. The multiplied results of Len(F Outlet ), Count(F Areas ), and F Low , which are recorded as f ( * ), will be used to define the level linearly, as "not urgent", "normal", and "urgent".
(2) Derive the urban flood-carrying capacity The storing ability is used to level the inner water capacity, the receiving ability expresses the upstream converging water compared to the storage capacity, the potential drainage ability describes the capacity downstream, and the emergency for pumping outward describes the local low-lying terrain which needs to be pumped. It can be concluded that the first three abilities are used to express the water capacity along the flow path, while the last ability is related to the pumping demand. The first three abilities are denoted the storing ability, the receiving ability, and the potential drainage ability and are categorized as: "very strong", "strong", "normal", "weak", and "very weak". Therefore, the previous three indicators shall be joined together by intersection operation, as terrain which needs to be pumped. It can be concluded that the first three abilities are used to express the water capacity along the flow path, while the last ability is related to the pumping demand. The first three abilities are denoted the storing ability, the receiving ability, and the potential drainage ability and are categorized as: "very strong", "strong", "normal", "weak", and "very weak". Therefore, the previous three indicators shall be joined together by intersection operation, as Based on the previous four aspects of basic capacity, the urban flood-carrying capacity can be leveled as: "weak", "normal", and "strong". The final urban flood-carrying capacity can be derived by: = ∩ , as shown in Equation (5).

Preparation for the estimation
(1) Extend the underlying surface model from VIS to VIS-W The underlying surface is a general term referring to the surface from the bottom of the atmospheric environment to the part in contact with the Earth's surface. In the research of the urban heat island effect field, the underlying surface model can generally be classified as forest land, grassland, cement land, and water surfaces, according to the sensitivity of heat reflection and temperature changes. In urban land use research, it is modeled using urban land, farmland, grassland, forest, water surface, and bare land. In urban hydrological research, it is generally divided into roof waterproofing material, sidewalk paving material, and grassland, according to the hydrological characteristics of the covering. However, the generic model is the VIS model, which was firstly proposed by Ridd et al. [23]. The VIS model is widely used in the urban hydrology field, such as in the surface water distribution and water exchange field. However, in urban flood-carrying capacity research, water type represents the basic element of storing entities, such as rivers, lakes, reservoirs, and ponds in the natural water system. Therefore, in order to estimate the capacity of storing rainwater properly, we added the water type into the VIS model, so as to extend the VIS model to the VIS-W (vegetation-impervious surface-soil and water) model.
(2) Pre-define the runoff coefficient Except for the rainwater lost in the process of stratum infiltration, vegetation interception, and surface land evapotranspiration, most of the rainwater will be surface runoff, which will need to be discharged into the water system. The hydrology model integrating the runoff coefficient method is widely adopted for rainwater calculating. The runoff coefficient is related to the local terrain and hydrology features, and it needs to be calibrated by many rainstorm and runoff observation experiments. Therefore, we suggest to adopt the official provided runoff coefficient values.
Based on the previous four aspects of basic capacity, the urban flood-carrying capacity can be leveled as: "weak", "normal", and "strong". The final urban flood-carrying capacity can be derived by: C carrying = A capacity ∩ A emergency , as shown in Equation (5).

Preparation for the Estimation
(1) Extend the underlying surface model from VIS to VIS-W The underlying surface is a general term referring to the surface from the bottom of the atmospheric environment to the part in contact with the Earth's surface. In the research of the urban heat island effect field, the underlying surface model can generally be classified as forest land, grassland, cement land, and water surfaces, according to the sensitivity of heat reflection and temperature changes. In urban land use research, it is modeled using urban land, farmland, grassland, forest, water surface, and bare land. In urban hydrological research, it is generally divided into roof waterproofing material, sidewalk paving material, and grassland, according to the hydrological characteristics of the covering. However, the generic model is the VIS model, which was firstly proposed by Ridd et al. [23]. The VIS model is widely used in the urban hydrology field, such as in the surface water distribution and water exchange field. However, in urban flood-carrying capacity research, water type represents the basic element of storing entities, such as rivers, lakes, reservoirs, and ponds in the natural water system. Therefore, in order to estimate the capacity of storing rainwater properly, we added the water type into the VIS model, so as to extend the VIS model to the VIS-W (vegetation-impervious surface-soil and water) model.
(2) Pre-define the runoff coefficient Except for the rainwater lost in the process of stratum infiltration, vegetation interception, and surface land evapotranspiration, most of the rainwater will be surface runoff, which will need to be discharged into the water system. The hydrology model integrating the runoff coefficient method is widely adopted for rainwater calculating. The runoff coefficient is related to the local terrain and hydrology features, and it needs to be calibrated by many rainstorm and runoff observation experiments. Therefore, we suggest to adopt the official provided runoff coefficient values.
The comprehensive runoff coefficient used in urban areas is usually between 0.5 and 0.8, and its range is (0.4, 0.6) in the suburbs [39]. Referring to the official definition from the Wuhan Municipal Flood Control and Drainage Regulations (version year 2013) [40], it was calculated for the VIS-W model. Compared to the VIS-W model, the type of vegetation (V) has a direct corresponding value, and we adopted 0.275, which is the average value of the proposed values. The type of imperious surface (I) adopts the average value of the corresponding resident land, business land, industry land, and traffic land, and its value is 0.775. The type of soil (S) adopts the corresponding value of the other land, and the value is 0.275. Furthermore, the type of water (W) always represents the area for storing rainwater. Therefore, its value shall be considered as 1, which means that all of the rainwater produced in the range of the water area will be collected completely by the water storage system. All of the final values we adopted in this paper are listed in Table 3. The quantity of runoff rainwater needs to be calculated by the rainfall-runoff hydrology model. We used the integrated runoff coefficient related to the underlying surface model to obtain the value. The area of land type in VIS-W was considered as the weight to obtain the integrated runoff coefficient, as shown in Equation (6).
where Ψ int represents the integrated runoff coefficient, which is the average of the runoff coefficient weight by the area ratio of land cover type in the VIS-W model, including the land cover types of vegetation, impervious surface, soil, and water; A i (unit: 10 4 m 2 ) represents the area of land type in the sub-catchment; and ψ i represents the runoff coefficient of land use types in the VIS-W model, including vegetation, impervious surface, soil, and water.
(3) Pre-define the rainfall precipitation The precipitation is the main input when calculating rainwater. In a city, the statistical precipitation according to the hydrological frequency analysis can reflect the city's common rainfall intensity. In China, the city's frequency adjustment of design storm intensity usually adopts an empirical frequency distribution curve model, such as the extreme value distribution curve, the negative exponential distribution curve, and the Pearson-III type distribution curve [39]. In this paper, we quoted Wuhan's design rainfall identity (the precipitation of four return periods: 1, 5, 20, and 100 years) proposed by Guoping Hong et al., based on the theory of the extreme value distribution curve [41], and the values are listed in Table 4. In order to qualitatively estimate the rainwater capacity of each sub-watershed, the four designed precipitation values are divided into six intervals, which can be described as interval set :{"<1", "1-5", "5-20", "20-100", "~100" and ">100"} as shown in Table 5. The middle four intervals are calculated by the mean value of the two neighboring sides of precipitation; the first interval "<1" has the minimum value of "1-5", and the last interval ">100" has a maximum value of "~100".
where R Int can utilize 6-h or 12-h duration precipitation to calculate the return period interval. Its interval values can be quantified by the interpolation of precipitation according to the previous four return periods. R Int(6h) represents the return period interval and its precipitation, which uses the duration of 6 hours as the basic calculating value. R Int(12h) represents the return period interval precipitation and is based on 12 hours of precipitation. The interval value "<1" adopts the value of 1-year precipitation as the maximum. The interval values "1-5", "5-20", and "20-100" use the left neighbor's maximum as the minimum, and adopt the average value of the two interval sides corresponding to precipitation as the maximum.
The fifth interval "~100" uses the maximum of the previous neighbor interval "20-100" as the minimum, and the value of the 100-year return period is accumulated as the average value of the precipitation values of the 20-year and 100-year intervals. The sixth interval ">100" adopts the maximum value of the interval "~100" as the minimum limitation.

(4) Construct basic maps
We mainly need to map the slope, land-use type, and the watershed, as shown in Figure 3. With the exception of the land-use type, which is extracted by ENVI 5.3, the other maps are produced using the related tools in ARCGIS 10.2.
3) Slope < 1%. The slope < 1% raster map is used to identify the low-lying area, and it will be used to analyze the areas which need to pump outside of the watershed. It is calculated by the DEM data set using the slope analysis tool and the raster calculation tool in ARCGIS 10.2.

(5). Construct basic indicators 1)
. Corresponding precipitation in each sub-catchment is used to express the storage capacity. It is composed of three parts: the river, the lake, and the puddle or pond.
where (unit: m 3 ) represents the rainwater storage, and the storage volume of the river, lake, and puddle or pond are represented by the _ (c) The area with slope less than 1%. It is colored in two classes: slope = 0 and 0 < slope < 1%.
(1) Land-use type. The ratio of land cover type is used to calculate the rainfall-runoff production based on the urban runoff coefficient in each sub-watershed. The land cover type is composed of the vegetation, impervious surface, soil, and water, based on the VIS-W model. The maximum likelihood method was used to extract the land cover type by ENVI 5.3 from the GF-1 and GF-2 images. We firstly extracted the vegetation, impervious surface, soil, water, and shadow. Then, according to the spatial adjacency between the shadow and the water, we could identify the part of shadow which should belong to water, and we regarded the left part of shadows as the impervious surface. This might bring a slight overestimation of the impervious surface, but we considered this as acceptable because the rainwater produced by the overestimated impervious surface is too small to disturb the collected water.
(2) Watershed. The watershed polygon is the defined basic evaluation unit. We used the hydrology tool of ARCGIS 10.2, which is based on the D8 algorithm, to delimitate watershed. Related to the previous research, we used a refined DEM by water body for watershed delimitation, since the performances will increase the certainty of flow direction and improve the accuracy of the delimited watershed compared to the original DEM [42]. The adopted water body was extracted from GF-1 and GF-2 by the previous land type map step. Based on the smallest area of the Sanjiaohu lake watershed [43], we used 11,111 (10 km 2 / 30 m × 30 m ≈ 11,111) as the limitation of flow accumulation to extract the stream network, and we got 38 sub-catchments, finally. Considering the condition that the inside area ratio must be higher than 40%, there are 15 sub-catchments that have been chosen. They are divided into two groups: the "very important" group, for which the ratio is more than 89% (including No. 8,No. 11,No. 15,No. 17,No. 18,No. 19,No. 20,and No. 26), and the "important" group, for which the ratio is between 40% and 89% (including No. 5,No. 6,No. 16,No. 22,No. 25,No. 27,and No. 28).
(3) Slope < 1%. The slope < 1% raster map is used to identify the low-lying area, and it will be used to analyze the areas which need to pump outside of the watershed. It is calculated by the DEM data set using the slope analysis tool and the raster calculation tool in ARCGIS 10.2.

(5) Construct basic indicators
(1) V Capacity . Corresponding precipitation in each sub-catchment is used to express the storage capacity. It is composed of three parts: the river, the lake, and the puddle or pond.
where V Capacity (unit: m 3 ) represents the rainwater storage, and the storage volume of the river, lake, and puddle or pond are represented by the R_V Capacity (unit: m 3 ), L_V Capacity (unit: m 3 ), and P_V Capacity (unit: m 3 ), respectively.

River storage
The potential storage of a river is determined by the limitation of the water level in flood seasons. In this paper, we take three meters as the water level increase limitation based on the flood control water level of the Yangtze River and Han River. Equation (9) is used to calculate this storage capacity: where R_V Capacityy (its unit is 10 6 m 3 ) represents the storage capacity of the river area in the sub-catchment, A i (unit 10 6 m 3 ) is the area of water belonging to the river, ρ represents the area coefficient of the calculated river, and h (unit is m) is the storing water depth. We default set the following: ρ = 1.0 and h = 3.0.

Lake storage
Zhou et al. [43] proposed the identification of the storage of lakes in terms of pre-regular pumping (the water level should be 0.3 m, and the area coefficient should be 0.8), normal storage (the water level should be 1 m, and the area coefficient should be 1), and extended water level (the water level should be 0.3 m, and the area coefficient should be 1.2).
Here, L_V Capacity (unit: 10 6 m 3 ) represents the capacity of lake storage, V normal (unit: 10 6 m 3 ) represents the storage of a normal storage capacity, V pre_pumping (unit: 10 6 m 3 ) represents the volume of the pre-regular pumping water level before flood seasons, and V expansion (unit: 10 6 m 3 ) represents the magnifying storage according to the surveyed average flood control water level. A i (unit: 10 6 m 2 ) represents the area of lakes in the sub-catchment, ρ represents the calculating area coefficient of lakes, and h (unit: m) represents the permission increasing water level in lakes. We default set ρ = 1.0 and h = 1.0 to calculate the normal storage, set ρ = 0.8 and h = 0.3 to calculate the pre-pumping storage, and set ρ = 1.2 and h = 0.3 to calculate the lakes' expansion storage.

Puddle or pond storage
Puddles or ponds can temporarily store rainwater in residential areas. In this study, we defined single areas of water bodies measuring between 900 square meters and 100,000 square meters as puddle or pond entities. This is because 900 square meters represents an area of about the horizon resolution of DEM (length 30 meters, and width 30 meters). Additionally, the area of 100,000 square meters is close to the area of Sanjiao Lake, which is the smallest lake in Wuhan. The distribution of pits or ponds is wide and relatively separate. We found the location of pits or ponds to be generally near residents, which demonstrated that there is no room to extend the water level higher than the shores. We default set the average increase in water level as 0.8 m.
Here, P_V Capacity (unit: 10 6 m 3 ) represents the storage, A i (unit: 10 6 m 2 ) represents the area of the puddle or pond in the sub-catchment, ρ represents the area coefficient for calculating the puddle or pond, and h (unit: m) represents the water level. We default set ρ = 1.0 and h = 0.8.

The storage capacity identify by the interval of rainfall intensity
According to the pre-defined runoff coefficient, the integrated runoff coefficient by land cover type area in each catchment can be estimated according to Equation (6). Finally, the storage capacity of each sub-catchment can be calculated as follows: where the rainfall intensity Precipitation (unit : mm ) is restricted by the rainwater storage V Capacity (unit: m 3 ); the area of sub-catchment is recorded as Area(unit: m 2 ); and the integrated runoff is the coefficient Ψ int . R Int(12h) or R Int(6h) represent the precipitation of a certain return period interval, which is calculated by Equation (7).
(2) Ai. This represents the area of vegetation, soil, impervious surface, and water. It can be calculated by the statistical analysis tool in ARCGIS 10.2. The sub-catchment polygon is used to extract the area of vegetation, soil, impervious surface, and water from the land used type raster by the geometry analysis function in ARCGIS 10.2.
(3) F Low . The ratio of low-lying areas represents the area for which the slope is less than 1%, except the area for which the land cover type is water.
(4) L River . The stream level represents the flow path inside the urbanized catchment. It adopts the river level function in the hydrology analysis tool in ArcGIS 10.2 to calculate its value, using the stream data extracted from the flow accumulation raster. The level of stream is named according to the Strahler level definition, which adopts the river link to the stream level from upstream to downstream. When the branches with the same grade of streams converge together, the level of current link will move up one grade, otherwise, the higher value is kept among its branches.
(5) F Terminaltype . The terminal flow type of the sub-watershed current is analyzed by the stream link and the flow path. It is labeled as "river" "inner", and "outer". The term "river" represents the downstream link finally pouring into river and the term "outer" represents the downstream link finally pouring outside of the catchment; otherwise, the terminal flow type can be labeled as "inner".
(6) F Nexttype . The sub-catchment of the next flow type represents the type of neighbor sub-catchment along the flow path. It is also labeled as "river", "inner", and "outer".
(7) Len(F Outlet ). The length is from current sub-catchment to the main outlet. In this paper, it is used to record the sub-watersheds to the river. If the terminal flow type is "outlet", the value of Len(F Outlet ) shall use a value twice the Len(F Outlet ) maximum values among the sub-watersheds for which the terminal flow types are labeled as "river" or "inner".
(8) Count(F Areas ). The count of sub-catchment follows the flow path from upstream. Its value is default set as 1, which means that the area only receives rainwater from itself.
Finally, all of the basic indicators in each sub-catchment are listed in Table 6.

Results
According to the proposed model, the urban flood-carrying capacity can be derived by the storing ability, the receiving ability, the drainage potential, and the pumping emergency. The flood-carrying capacity estimating result is shown in Table 7. We find that the previous three abilities are summarized strictly by the lowest level, and the pumping emergency has the highest prior affectation in the final estimation results. The final result will be moved down one grade by the state of emergency compared with the previous three abilities combined. The state of "not urgent" keeps the original level unchanged. Taking No.5 and No.25 for example, for both the pumping emergency is labeled as "normal". Meanwhile, the original levels of No.5 and No.25 are "very strong" and "strong", respectively. Therefore, the finally flood-carrying capacities are "strong" and "normal", respectively. As shown in Table 7, there are four, five, and six sub-catchments whose capacities are labeled as "strong", "weak", and "normal", respectively. Figure 4 shows the distribution of the flood-carrying capacity. We can determine that the "strong" sub-catchment (e.g., No.5,No.8,No.20,and No.28) contains a part of river area. We can also find that the flood-carrying capacity is determined by the distance to the outlet, its rain water storage advantage, the receiving demand along the stream flow, and its low-lying areas. For example, comparing the distance to the river, both sub-catchments No. 25 and No.27 neighbor the river, but the final state of No.25 is normal, and that of No.27 is weak. As we can see in Figure 4

Comparing the Carrying Capacity with Waterlogging Points
According to the rainfall intensity during 12-hour in Wuhan shown in Table 4, the return period value in 5-20 years is 107.35-171.85 mm, and the return period values for 5 years, 20 years, and 100 years are 138.1 mm, 205.6 mm, and 336.2 mm, respectively. In the flood disaster from June 30 to July 6, 2016, the cumulative precipitation value on severe days was 583.1 mm, according to the

Comparing the Carrying Capacity with Waterlogging Points
According to the rainfall intensity during 12-hour in Wuhan shown in Table 4, the return period value in 5-20 years is 107. .85 mm, and the return period values for 5 years, 20 years, and 100 years are 138.1 mm, 205.6 mm, and 336.2 mm, respectively. In the flood disaster from 30 June to 6 July 2016, the cumulative precipitation value on severe days was 583.1 mm, according to the precipitation record of Station No. 57,494. The cumulative precipitation passed the limitations of the 100-year return period.
Indeed, as per the waterlogging point distribution released by Wuhan officials on July 6, 2016 (http://www.whwater.gov.cn/water/fxkht/7958.html, accessed: November 5, 2017), this flood and waterlogging event caused 99 waterlogging places in central Wuhan city, and all the water depths are too high to be convenient for motor vehicles to pass through. Inside Third Ring Road of Wuhan, there were 40 waterlogging points not suitable for motor vehicles to pass through, as shown in Figure 5. We found that 30 of them were located in the range of the study area. These points were grouped as "large cars cannot pass" and "small cars cannot pass". Here, "large cars" refer to large vehicles like buses and trucks, as compared to small household cars. The points labeled as "large car cannot pass" indicate none of the cars can pass through, while in the "small car cannot pass" range, only the large cars can pass through. As shown in Figure 5, there are a total of eight sub-catchments located inside Third Ring Road of Wuhan. The flood-carrying capacity in two of them are "weak", four of them are "normal", and the remaining two are "strong". On comparing the waterlogging distribution to the flood-carrying capacity of sub-watersheds, one of the "weak" sub-catchments (No.6) suffering from waterlogging was labeled as "large car cannot pass" and "small cars cannot pass", and one of the "normal" sub-catchments (No.18) has the same situation as No.6. There are four "normal" sub-catchments (No. 19,No.20,No.25,and No.26) which have waterlogging points labeled as "small car cannot pass". Furthermore, there was one "strong" sub-catchment (No.8) and one "weak" sub-catchment (No.27) that did not have waterlogging point, respectively. The accumulative waterlogging points of three level are listed in Table 8. Table 8. The typical waterlogging points distribution inside Third Ring Road of Wuhan. As listed in the number of sub-catchments, there are eight sub-catchments inside Third Ring Road, and there are two sub-catchments that have no waterlogging points. The columns "Large car cannot pass" and "Small car cannot pass" indicate different depths of waterlogging points, where the depth in the former is higher than in the latter. The distribution of waterlogging section lists the quality of points in three classes of sub-catchment, respectively. cannot pass". Furthermore, there was one "strong" sub-catchment (No.8) and one "weak" subcatchment (No.27) that did not have waterlogging point, respectively. The accumulative waterlogging points of three level are listed in Table 8. Figure 5. The typical waterlogging points distributed as overlapping the flood-carrying capacity result based on the actual data from the flood and waterlogging event which took place from June 30 to July 6, 2016. In the 30 waterlogging points, there were five points were not suitable for large motor vehicles to pass through, and 25 points were not suitable for small motor vehicles to pass through. Table 8. The typical waterlogging points distribution inside Third Ring Road of Wuhan. As listed in the number of sub-catchments, there are eight sub-catchments inside Third Ring Road, and there are two sub-catchments that have no waterlogging points. The columns "Large car cannot pass" and "Small car cannot pass" indicate different depths of waterlogging points, where the depth in the former is higher than in the latter. The distribution of waterlogging section lists the quality of points in three classes of sub-catchment, respectively. According to the statistical waterlogging points in Table 7, we can find that firstly, in terms of the sub-watershed, 100% of the "large cars cannot pass" waterlogging points occurred in the "weak" or "normal" sub-catchments. We also noticed that in the two sub-catchments for which the floodcarrying capacity level is "weak", No.27 had no waterlogging points. Compared to the real drainage pumping station, the Tangxun Lake pumping station is located in No.27, and it fully runs during the flood days. This may be the reason that this "weak" area met no waterlogging points at all. We can also see that 80% of the waterlogging points labeled as "small cars cannot pass" took place in the "weak" and "normal" sub-catchments, but there are still some waterlogging points in the "strong" According to the statistical waterlogging points in Table 7, we can find that firstly, in terms of the sub-watershed, 100% of the "large cars cannot pass" waterlogging points occurred in the "weak" or "normal" sub-catchments. We also noticed that in the two sub-catchments for which the flood-carrying capacity level is "weak", No.27 had no waterlogging points. Compared to the real drainage pumping station, the Tangxun Lake pumping station is located in No.27, and it fully runs during the flood days. This may be the reason that this "weak" area met no waterlogging points at all. We can also see that 80% of the waterlogging points labeled as "small cars cannot pass" took place in the "weak" and "normal" sub-catchments, but there are still some waterlogging points in the "strong" sub-catchment (No. 20). Secondly, in terms of the numbers of waterlogging points, as shown in Figure 6, we can see that 84% of the points located in the "weak" and "normal" sub-watershed, and 100% of the points labeled "large cars cannot pass" were located in the "weak" or "normal" sub-watersheds. As the accumulative precipitation nearly passed the return period of 100 years during the flood and waterlogging event, the waterlogging points which take place in areas with flood-carrying capacities of "weak" or "normal" are understandable. sub-catchment (No. 20). Secondly, in terms of the numbers of waterlogging points, as shown in Figure  6, we can see that 84% of the points located in the "weak" and "normal" sub-watershed, and 100% of the points labeled "large cars cannot pass" were located in the "weak" or "normal" sub-watersheds. As the accumulative precipitation nearly passed the return period of 100 years during the flood and waterlogging event, the waterlogging points which take place in areas with flood-carrying capacities of "weak" or "normal" are understandable.

Amount of Sub-Catchments Distribution of Waterlogging
. Figure 6. The ratio of waterlogging points grouped by "small cars cannot pass" and "large cars cannot pass" according to the location in the sub-catchment with "weak", "normal", and "strong" capabilities.
Generally, over 80% of the typical waterlogging points are located in area whose flood-carrying capacity states are "weak" or "normal". This statistical result indicates the urban flood-carrying capacity estimated by the proposed model is reasonable, and the proposed model can be used to estimate the urban flood-carrying capacity.

Discussion on the Storing Ability of Each Sub-Catchment in the Waterlogging Event
We further analyzed the potential capacity of the natural water system according to the flood and waterlogging disaster event. Based on the hydrological observation data, from 00:00 h on June 30, 2016 to 00:00 h on July 12, 2016 (Station: Hankou Station), we found that the original water level of Yangtze River on June 30, 2016 was 25 m, and the maximum level reached 28.37 m at 04:00 h on July 7, 2016. This record verified that the storage rainwater ability of the Yangtze River can be calculated by the increase in the limitation water level of 3 m. Additionally, the increased water level in the lakes in the direct pre-regular plan is 1 m, the normal rising level is 0.3 m, and the extended flood control level is 0.3 m, according to the flood control for Wuhan proposed in the literature [43]. Figure 6. The ratio of waterlogging points grouped by "small cars cannot pass" and "large cars cannot pass" according to the location in the sub-catchment with "weak", "normal", and "strong" capabilities.
Generally, over 80% of the typical waterlogging points are located in area whose flood-carrying capacity states are "weak" or "normal". This statistical result indicates the urban flood-carrying capacity estimated by the proposed model is reasonable, and the proposed model can be used to estimate the urban flood-carrying capacity.

Discussion on the Storing Ability of Each Sub-Catchment in the Waterlogging Event
We further analyzed the potential capacity of the natural water system according to the flood and waterlogging disaster event. Based on the hydrological observation data, from 00:00 h on 30 June 2016 to 00:00 h on 12 July 2016 (Station: Hankou Station), we found that the original water level of Yangtze River on 30 June 2016 was 25 m, and the maximum level reached 28.37 m at 04:00 h on 7 July 2016. This record verified that the storage rainwater ability of the Yangtze River can be calculated by the increase in the limitation water level of 3 m. Additionally, the increased water level in the lakes in the direct pre-regular plan is 1 m, the normal rising level is 0.3 m, and the extended flood control level is 0.3 m, according to the flood control for Wuhan proposed in the literature [43].
We finally calculated the result of the quantity of runoff based on the rainfall-runoff hydrological model by integrating the runoff coefficient of the underlying surface, which was about 235 million m 3 in the flood and waterlogging disaster event which occurred from 30 June 2016 to 6 July 2016. Meanwhile, the whole storage capacity was about 298 million m 3 . Therefore, the total balance of the flood-carrying capacity is higher than that of the runoff rainwater. In theory, the whole city can avoid waterlogging if the storing capacity can be regulated as freely as possible.
As shown in Figure 7, the rainwater storage capacities in the waterlogging disaster in sub-catchments No.6 (in red color), No. 22,No.25,and No.27 (in light brown colors), and No. 11,No.18,No.16,No.17,and No.26 (in light green colors), are overloaded. The storage capacities of sub-catchments No. 15 and No.19 in light blue color are slightly greater than the quantity of rainwater. Furthermore, the storage capacities in No. 5,No.8,and No.20 (in light purple colors) are much higher than the quantity of rainwater. We finally calculated the result of the quantity of runoff based on the rainfall-runoff hydrological model by integrating the runoff coefficient of the underlying surface, which was about 235 million m 3 in the flood and waterlogging disaster event which occurred from June 30, 2016 to July 6, 2016. Meanwhile, the whole storage capacity was about 298 million m 3 . Therefore, the total balance of the flood-carrying capacity is higher than that of the runoff rainwater. In theory, the whole city can avoid waterlogging if the storing capacity can be regulated as freely as possible.
As shown in Figure 7, the rainwater storage capacities in the waterlogging disaster in subcatchments No.6 (in red color), No. 22,No.25,and No.27 (in light brown colors), and No. 11,No.18,No.16,No.17,and No.26 (in light green colors), are overloaded. The storage capacities of subcatchments No. 15 and No.19 in light blue color are slightly greater than the quantity of rainwater. Furthermore, the storage capacities in No. 5,No.8,and No.20 (in light purple colors) are much higher than the quantity of rainwater.
Compared to the self-container storage capacity in the sub-catchment, the seven waterlogged places located in No. 5,No.8,and No.20, and the three waterlogged places located in No. 19, are not reasonable. The waterlogged places located in No. 25,No.26,and No.18 can be understood in terms of the natural water system as having a lesser capacity compared to the quantity of rainwater. The five waterlogged places located in the No.6 sub-catchment are reasonable. This is because the local storage of the natural water entity is not plentiful. According to the stream network, the flow path of the sub-catchment is finally outside of the study areas whose direction is opposite to the main outlet of the Yangtze River or Han River.
In our opinion, the waterlogging disaster could have been mitigated by a combination of fully using the natural drainage capacity and a reasonable drainage schedule of the city. According to the standards of the Wuhan government, if there exists waterlogging of higher than 0.4 m for more than 2 hours, the level of waterlogging can be defined as serious. In the serious flood and waterlogging disaster used in this paper, we believe the ability to regulate the natural water system is not sufficient. Taking catchment No.26 as an example, the longest duration of waterlogging is located here. From the city's drainage network, the rainwater naturally flows to No.27 and is finally pumped into the Compared to the self-container storage capacity in the sub-catchment, the seven waterlogged places located in No. 5,No.8,and No.20, and the three waterlogged places located in No. 19, are not reasonable. The waterlogged places located in No. 25,No.26,and No.18 can be understood in terms of the natural water system as having a lesser capacity compared to the quantity of rainwater. The five waterlogged places located in the No.6 sub-catchment are reasonable. This is because the local storage of the natural water entity is not plentiful. According to the stream network, the flow path of the sub-catchment is finally outside of the study areas whose direction is opposite to the main outlet of the Yangtze River or Han River.
In our opinion, the waterlogging disaster could have been mitigated by a combination of fully using the natural drainage capacity and a reasonable drainage schedule of the city. According to the standards of the Wuhan government, if there exists waterlogging of higher than 0.4 m for more than 2 h, the level of waterlogging can be defined as serious. In the serious flood and waterlogging disaster used in this paper, we believe the ability to regulate the natural water system is not sufficient. Taking catchment No.26 as an example, the longest duration of waterlogging is located here. From the city's drainage network, the rainwater naturally flows to No.27 and is finally pumped into the Yangtze River at the Tangxun lake pump station. During the waterlogging event, sub-catchment No.26 produced about 22 million m 3 of runoff rainwater, and its maximum regulated storing capacity is 13 million m 3 . Therefore, there remained 9 million m 3 of rainwater that needed to be discharged. According to the official information, the designed drainage capacity of the Tangxun Lake pump station is 120 m 3 /s, which means that the pump system needs to work for about 21 h to pump rainwater into the city's drainage system. However, in the pump work report, the Tangxun Lake pump station ran for 90 h from July 1 to 6. Meanwhile, we have found that there are six waterlogging places located in sub-catchment No.26.

Analyzing the Carrying Capacity with Land Type and Terrian Factors
According to the estimation model expressed in Section 2.3, we can see that the urban flood-carrying capacity corresponds with the underlying surface. For example, the terrain characteristics will affect the delimitation of the watershed and influence the flow path. The land cover type will determine the ratio of runoff yield, and the potential of discharge will affect the degree of waterlogging risk. As shown in Figure 8, we adopted statistical graphs to discuss the relation between the estimated results and the underlying surface characteristics. Figure 8a,b used the main axis on the left to display comparison factors (land cover type ratio, storing ability transforming into precipitation measured by runoff yield, elevation range, and slope), and adopted "1", "2", and "3" to map the carrying capacity states of "weak", "normal", and "strong", respectively, by the secondary axis. Taking catchment No.26 as an example, the longest duration of waterlogging is located here. From the city's drainage network, the rainwater naturally flows to No.27 and is finally pumped into the Yangtze River at the Tangxun lake pump station. During the waterlogging event, sub-catchment No.26 produced about 22 million m 3 of runoff rainwater, and its maximum regulated storing capacity is 13 million m 3 . Therefore, there remained 9 million m 3 of rainwater that needed to be discharged. According to the official information, the designed drainage capacity of the Tangxun Lake pump station is 120 m 3 /s, which means that the pump system needs to work for about 21 hours to pump rainwater into the city's drainage system. However, in the pump work report, the Tangxun Lake pump station ran for 90 hours from July 1 to 6. Meanwhile, we have found that there are six waterlogging places located in sub-catchment No.26.

Analyzing the Carrying Capacity with Land Type and Terrian Factors
According to the estimation model expressed in section 2.3, we can see that the urban floodcarrying capacity corresponds with the underlying surface. For example, the terrain characteristics will affect the delimitation of the watershed and influence the flow path. The land cover type will determine the ratio of runoff yield, and the potential of discharge will affect the degree of waterlogging risk. As shown in Figure 8, we adopted statistical graphs to discuss the relation between the estimated results and the underlying surface characteristics. Figure 8(a) and Figure 8(b) used the main axis on the left to display comparison factors (land cover type ratio, storing ability transforming into precipitation measured by runoff yield, elevation range, and slope), and adopted "1", "2", and "3" to map the carrying capacity states of "weak", "normal", and "strong", respectively, by the secondary axis. Figure 8(a) uses the cumulative ratio of land use type to draw a histogram chart. In Figure 8(a), we can find that the "weak" areas always have a similar shape to areas with a low ratio of the water body, while most of the "strong" areas are similar to those with a high ratio of the water body. Furthermore, the "normal" areas are not very significant in terms of the proportion of water. These phenomena indicate that the flood-carrying capacity of each sub-catchment is positively correlated with the water ratio. Moreover, the water area is the land cover of open water extracted by the remote sensing images. This indicates that the carrying capacity can reflect the storing rainwater ability to regulate the natural water system. Figure 8(b) uses the logarithmic function with a base of 2 to pre-calculate the underlying surface factors, and finally displays the results as a line chart. In Figure 8(b), we used two kinds of terrain factors (ratio of slope < 1% areas and elevation range) and one kind of rainfall factor (runoff) for analysis. Both of the terrain factors reflect the change of elevation from different aspects. The slope is the average value of elevation, which shows the local land relief character, and can distinguish the ratio of low-lying areas. The elevation range is the regional topographic drop, and can reflect the flow path of rainwater. The line shape and the amplitude of terrain factors are different. From Figure 8 Figure 8a uses the cumulative ratio of land use type to draw a histogram chart. In Figure 8a, we can find that the "weak" areas always have a similar shape to areas with a low ratio of the water body, while most of the "strong" areas are similar to those with a high ratio of the water body. Furthermore, the "normal" areas are not very significant in terms of the proportion of water. These phenomena indicate that the flood-carrying capacity of each sub-catchment is positively correlated with the water ratio. Moreover, the water area is the land cover of open water extracted by the remote sensing images. This indicates that the carrying capacity can reflect the storing rainwater ability to regulate the natural water system. Figure 8b uses the logarithmic function with a base of 2 to pre-calculate the underlying surface factors, and finally displays the results as a line chart. In Figure 8b, we used two kinds of terrain factors (ratio of slope < 1% areas and elevation range) and one kind of rainfall factor (runoff) for analysis. Both of the terrain factors reflect the change of elevation from different aspects. The slope is the average value of elevation, which shows the local land relief character, and can distinguish the ratio of low-lying areas. The elevation range is the regional topographic drop, and can reflect the flow path of rainwater. The line shape and the amplitude of terrain factors are different. From Figure 8b, we can find that the storing runoff ability is similar to the shape of the slope to a certain degree. The reason for this is that the ratio of low-lying land is an indicator of the storing ability. Additionally, the area with a "high" carrying capacity always has a peak value of storing ability of runoff. Sub-catchment No.25 does not have not the local peak runoff ability, but its slope of less than 1% is low and the elevation range is very high, which shows that the area might flow into a neighboring area. Therefore, its flood-carrying capacity is "strong". The phenomenon indicates that the estimating model is synthesized considering the water areas, runoff coefficient, and flow path.

Conclusions
The study proposed a novel model based on the underlying surface VIS-W model to estimate the urban flood-carrying capacity in terms of the consistency with the land use type, the terrain features, the storage, receiving and discharge ability, and pumping outward emergency of the natural water system, through (1) extending the VIS underlying surfaces model to form the VIS-W model, (2) delimiting the sub-watershed of the urban catchment by DEM data set constrained by the water body, and (3) calculating the rainwater budget with the storage capacity of rivers, lakes, and puddles or ponds in each sub-catchment by a rainfall-runoff model; combined with the hydrological features of the terrain, the flood-carrying capacity could be comprehensively derived.
On comparing to the verified result with the flood-carrying capacity and typical waterlogging points inside Third Ring Road of Wuhan, we found 80% of sub-catchments with flood-carrying capacity categorized as "weak" or "normal" suffered waterlogging disasters, and 84% were located in sub-watersheds labeled as "weak" or "normal". This proves that the proposed model could be suitable for estimating the flood-carrying capacity in cities like Wuhan. Furthermore, on comparing the surplus state of the storing ability, the inconsistency between the waterlogging places and the flood-carrying capacity indicated that the storing ability might have some space for improvement. In our opinion, waterlogging disasters might be mitigated by fully utilizing the natural storage capacity to coordinate the normal flood drainage schedule.
In future research, the potential for regulating the natural water system in urbanized catchments should be analyzed. Additionally, rainfall scenes simulated by 1D or 2D hydrodynamic models should be overlapped to analyze the drainage performances.