Assessment of Community Vulnerability to Different Types of Urban Floods: A Case for Lishui City, China

Urban flooding is a severe and pervasive hazard caused by climate change, urbanization, and limitations of municipal drainage systems. Cities face risks from different types of floods, depending on various geographical, environmental, and hydrometeorological conditions. In response to the growing threat of urban flooding, a better understanding of urban flood vulnerability is needed. In this study, a comprehensive method was developed to evaluate the vulnerability of different types of urban floods. First, a coupled urban flood model was built to obtain the extent of influence of various flood scenarios caused by rainfall and river levee overtopping. Second, an assessment framework for urban flood vulnerability based on an indicator method was used to evaluate the vulnerability in different flood hazard scenarios. Finally, the method was applied to Lishui City, China, and the distribution and pattern of urban flood vulnerability were studied. The results highlight the spatial variability of flooding and the vulnerability distributions of different types of urban floods. Compound floods were identified to cause more severe effects in the urban areas.


Introduction
Floods, which can cause widespread social and economic damage and even threaten human life, are one of the most severe and frequent natural hazards in urban areas [1,2]. Under rapid urbanization and climate change, the effects of the flooding are becoming worse in urban areas [2][3][4]. Reducing the effects of urban flooding on communities has become an essential requirement in urban disaster management [5,6]. Many cities and communities have tried to mitigate the effects of urban flooding through structural measures (e.g., increasing urban drainage capacity and building dams) and non-structural measures (e.g., enhancing flood risk awareness) [7][8][9]. Both approaches rely on the knowledge and resources of multiple sectors and communities in terms of risk reduction strategies, experiences, and lessons; knowledge sharing can be especially valuable among communities with similar vulnerability characteristics [10]. Therefore, vulnerability assessment has become an essential tool for communities to cope with the adverse effects of urban flooding.
Urban communities in complex hydrometeorological environments can face multiple threats from various types of flood, such as urban pluvial floods and fluvial floods [11]. Many studies have considered the risks and vulnerabilities of urban flooding; these research works have primarily focused on fluvial floods, which are considered to result in severe consequences. However, the impacts of other types of floods can be comparable to those of fluvial floods, such as pluvial flooding when rainfall reaches a certain level [12]. Different types of urban flooding are gradually being given more attention. Furthermore, flood risk and vulnerability for different types of floods have been addressed in several studies [13,14]; however, in each study, the different types of floods were considered sensitivity, and adaptive capacity was built using the weighted linear combination (WLC) method. Lishui City, China was selected for a case study for scenario analyses.

Study Area
The Liandu District, the main urban area of Lishui City, China, was selected as the study area. Lishui City is located in southeast Zhejiang Province. Liandu District is a unique municipal district, and 25 communities were selected for analysis ( Figure 1). The community, like a village, is the smallest administrative and demographic unit in China and is divided by local administrators. The study area covered approximately 45 km 2 and is characterized by low-lying topography with hills and mountains to the north and west and the Oujiang River to the south and east. indicators, and fully use the geospatial data. Finally, an integrated flood vulnerability index composed of exposure, sensitivity, and adaptive capacity was built using the weighted linear combination (WLC) method. Lishui City, China was selected for a case study for scenario analyses.

Study Area
The Liandu District, the main urban area of Lishui City, China, was selected as the study area. Lishui City is located in southeast Zhejiang Province. Liandu District is a unique municipal district, and 25 communities were selected for analysis ( Figure 1). The community, like a village, is the smallest administrative and demographic unit in China and is divided by local administrators. The study area covered approximately 45 km 2 and is characterized by low-lying topography with hills and mountains to the north and west and the Oujiang River to the south and east. From a hydrometeorological perspective, the area is vulnerable to different types of floods. First, the study area is a substantial fluvial-flood-affected area located in the confluence of the Haoxi River and the Daxi River in the middle and lower reaches of the Oujiang River basin. Upstream, the area is mountainous, and the river discharge is often very high during the rainy season. Liandu District features a subtropical monsoon climate and is high in annual precipitation (~1733 mm). June and July mark the rainy season, during which frequent and continuous rainfall occurs. During this season, local plums mature; thus, the period is called the plum rain season or East Asian rainy season, as it is the common climatic phenomenon in the eastern coastal regions of China. Although the city is not located in a coastal region, it is often affected by typhoons owing to its proximity to the East China Sea. The most recent flood event in the study area occurred on 20 August 2014, when From a hydrometeorological perspective, the area is vulnerable to different types of floods. First, the study area is a substantial fluvial-flood-affected area located in the confluence of the Haoxi River and the Daxi River in the middle and lower reaches of the Oujiang River basin. Upstream, the area is mountainous, and the river discharge is often very high during the rainy season. Liandu District features a subtropical monsoon climate and is high in annual precipitation (~1733 mm). June and July mark the rainy season, during which frequent and continuous rainfall occurs. During this season, local plums mature; thus, the period is called the plum rain season or East Asian rainy season, as it is the common climatic phenomenon in the eastern coastal regions of China. Although the city is not located in a coastal region, it is often affected by typhoons owing to its proximity to the East China Sea.
The most recent flood event in the study area occurred on 20 August 2014, when the entire Oujiang basin was hit by heavy and durative rainfall that induced a 50-year flood and caused considerable loss of economic property in the urban area.

Study Data
The vulnerability assessment was implemented using an indicator-based method, where indicators were derived from the results of flood modeling and existing physical, social, and economic data. Datasets were mainly used to create models with different flood scenarios and physical, social, and economic indicators for the vulnerability assessment.
A large amount of geographic and hydrometeorological data were used for the construction of the flood models. Terrain data and municipal drainage system data formed the essential framework of the flood models in this study. Terrain data, including digital elevation model data with 5 m resolution and river bathymetry data, were interpolated from river cross-section survey data provided by the Lishui Geographic Information Center (LSIC). Municipal drainage data, including central lines of rivers and channels, rainwater pipelines, and other drainage facilities, such as pumps and manholes, were also provided by the LSIC. Moreover, hydrometeorological data, which mainly include rainfall and river discharge, are important input parameters for flood modeling. Rainfall data were calculated using the local rainstorm intensity formula provided by the Lishui Meteorological Bureau: where q (mm/min) is the rainfall intensity, t (min) is the duration of rainfall, and P (years) is the return period of rainfall. The time series of river discharge and stage, including upstream and downstream, were obtained from the records of the local hydrological observatory during the latest 50-year flood on 20 August 2014. Except for flood modeling, the indicators used to build indices were derived from abundant physical and socioeconomic data. The physical data were used to evaluate receptors that can be affected by floods, including buildings and various infrastructures, and most of these data were spatial data with locations and relevant attribute information. An example is the building data including the footprints and attributes, such as floor number, building area, and usage. Moreover, POI data provided by the LSIC were used as an important data source for increased access to urban infrastructure data. Thus, sufficient infrastructure information, such as the locations and attributes of retail stores, schools, nursing homes, plants, and basements, were extracted from POI data used to produce relevant indicators. Much of the social and economic information was collected from local statistical yearbooks and interviews with local community staff.

Methodological Outline
According to a report of the Intergovernmental Panel on Climate Change (IPCC), vulnerability can be regarded as a function of exposure, sensitivity, and adaptive capacity-expressed as Equation (2) [22,34,35]-and we use this as the framework of this study: where E represents exposure, S represents sensitivity, and AC represents adaptive capacity. For flood disasters, vulnerability is the tendency towards an adverse effect of the flood. Exposure refers to the elements exposed to flood hazards in a specific region, such as people, buildings, infrastructure, and other personal or public assets. Many studies regard that all receptors can be affected by a disaster in the study area; that is, all receptors should be used for exposure assessment, no matter where they are located [10,36]. In our view, exposure should depend on a specific disaster. For example, the former approach is suitable for disasters with widespread effects, such as earthquakes.
However, for disasters where the scope of the effect can be determined, such as floods, the exposure of a population or infrastructure is determined by the extent of flooding and can vary under different flood hazards. Sensitivity is the intrinsic predisposition of an element to be affected by or susceptible to damage with the occurrence of a flood hazard. Adaptive capacity represents the ability of communities, institutions, or people to adjust to the flood, take advantage of opportunities, or respond to consequences. Based on this understanding of vulnerability to urban flooding, the research outline for this study was developed ( Figure 2). Different types of urban flooding were simulated using a coupled flood model, which was built using LISFLOOD-FP [37] and the stormwater management model (SWMM) [38]. Flood simulation results were combined with variables to create a set of vulnerability indicators. Components of vulnerability, including exposure, sensibility, and adaptive capacity, were combined from the corresponding indicators using the AHP method. Finally, the vulnerability index was integrated from its components by using the WLC method. In this approach, GIS techniques and methods are frequently used in flood simulation and assessment as a spatial data processing method. affected by or susceptible to damage with the occurrence of a flood hazard. Adaptive capacity represents the ability of communities, institutions, or people to adjust to the flood, take advantage of opportunities, or respond to consequences. Based on this understanding of vulnerability to urban flooding, the research outline for this study was developed ( Figure 2). Different types of urban flooding were simulated using a coupled flood model, which was built using LISFLOOD-FP [37] and the stormwater management model (SWMM) [38]. Flood simulation results were combined with variables to create a set of vulnerability indicators. Components of vulnerability, including exposure, sensibility, and adaptive capacity, were combined from the corresponding indicators using the AHP method. Finally, the vulnerability index was integrated from its components by using the WLC method. In this approach, GIS techniques and methods are frequently used in flood simulation and assessment as a spatial data processing method.

Coupled Urban Flood Model
The SWMM is a dynamic urban flood simulation model widely used for the simulation of flooding and water quality in urban areas. The SWMM model includes a hydrological rain-runoff model and a 1D hydrodynamic model for simulation of municipal drainage networks, and is frequently used for urban planning, analysis, and design relevant to urban stormwater management. Although the SWMM can handle the interaction between surface runoff and sewer network runoff, the model cannot intuitively reveal the condition of flood inundation on the surface caused by overflow when the runoff exceeds the capacity of the rainwater drainage networks. A 2D hydrodynamic model, LISFLOOD-FP, was coupled with the SWMM model to solve this problem. LISFLOOD-FP is a highly efficient 2D hydrodynamic model used to simulate surface land inundation and fluvial floods; it can adapt to various complex terrains suited for the urban environment. Moreover, LISFLOOD-FP is a raster-bashed model; the regular grid is used as a computational unit to simulate 2D flood inundation and the sub-grid channel method is used to simulate the 1D river. Details of the coupling technique and implementation can be found in the literature [31][32][33]. To ensure model accuracy, our model coupling LISFLOOD-FP and SWMM was validated using the observation records of the 2014 50-year flood. Meanwhile, field observations were carried out in communities 1 and 13, in which the urban area is ungauged. The verification process is described in detail in the literature [39,40]. The coupled urban flood model, which

Coupled Urban Flood Model
The SWMM is a dynamic urban flood simulation model widely used for the simulation of flooding and water quality in urban areas. The SWMM model includes a hydrological rain-runoff model and a 1D hydrodynamic model for simulation of municipal drainage networks, and is frequently used for urban planning, analysis, and design relevant to urban stormwater management. Although the SWMM can handle the interaction between surface runoff and sewer network runoff, the model cannot intuitively reveal the condition of flood inundation on the surface caused by overflow when the runoff exceeds the capacity of the rainwater drainage networks. A 2D hydrodynamic model, LISFLOOD-FP, was coupled with the SWMM model to solve this problem. LISFLOOD-FP is a highly efficient 2D hydrodynamic model used to simulate surface land inundation and fluvial floods; it can adapt to various complex terrains suited for the urban environment. Moreover, LISFLOOD-FP is a raster-bashed model; the regular grid is used as a computational unit to simulate 2D flood inundation and the sub-grid channel method is used to simulate the 1D river. Details of the coupling technique and implementation can be found in the literature [31][32][33]. To ensure model accuracy, our model coupling LISFLOOD-FP and SWMM was validated using the observation records of the 2014 50-year flood. Meanwhile, field observations were carried out in communities 1 and 13, in which the urban area is ungauged. The verification process is described in detail in the literature [39,40]. The coupled urban flood model, which integrates characteristics of different models, not only adapts to the urban environment, but can also express urban floods intuitively and visually, which is important for urban risk and vulnerability analysis. In this study, the model was applied to urban flooding under different scenarios.

Urban Flood Scenarios
Considering the local physical environment and hydrometeorological characteristics, three types of floods were considered (pluvial, fluvial, and compound). Pluvial floods were induced by storm rainfall and included seven scenarios based on Equation (1), namely 1-, 5-, 10-, 20-, 30-, 50-, and 100-year rainfall. Fluvial floods were induced by river levee overtopping and were simulated based on the 50-year flood event of 20 August 2014. Compound flooding was induced by rainfall and river overtopping. The compound flood scenarios were simulated with different rainfall events and the river flood event of 20 August 2014. Meanwhile, a method based on GIS was used to directly combine the simulation results of urban pluvial flooding and fluvial flooding to allow for comparison with the flood modeling method. Therefore, compound flooding was obtained by the flood modeling and GIS methods.

Vulnerability Assessment
The indicator-based method was used for vulnerability assessment. The indicators were selected and integrated based on the vulnerability framework, as shown in Equation (2). The selected vulnerability indicators were divided into three, namely exposure, sensitivity, and adaptive capacity. Finally, the sub-indices were integrated into a comprehensive vulnerability index.
Using a method of combining the flood model and GIS techniques, the level of exposure was determined by the receptors (i.e., the factors that suffer adverse effects from urban flooding within the community, such as populations and buildings). In this method, the results of flood modeling were simulated for the given flood scenario, including water depth and extent of flooding, and then overlaid on other spatial data layers, such as populations and buildings, to determine the affected elements. Then, the exposure elements were identified and counted in each neighborhood based on GIS methods and selected as exposure variates.
Sensitivity is related to physical, social, and economic weaknesses and disadvantages. For urban flooding, sensitivity indicators should indicate vulnerable places and socioeconomic receptors for the flood disaster. Moreover, the sensitivity indicator was divided into two categories, namely physical and socioeconomic sensitivity. For physical sensitivity, indicators were used to reveal areas and places vulnerable to flooding. Several flood-related indicators were selected, including elevation, slope, distance to channels, and reservoirs, to identify sensitive areas. Then, places vulnerable to urban flooding in each community, such as elementary schools, nursing homes, basements, and chemical plants, were counted. For socioeconomic sensitivity, the sensitive population, including the disabled, aged, children, and those on a low-income (i.e., those who receive social assistance), was selected as an indicator; the data on these groups were obtained from the Lishui City Statistics Bureau.
Adaptive capacity is the ability to adapt to a disaster. In this study, the convenience of obtaining rescue and help was used to assess the adaptive capacity when facing the threat of floods. Therefore, rescue and medical services, including hospitals, shelters, and police and fire stations, were selected, and these data were converted into indicators based on a GIS method of Euclidean distance analysis.
The selected indicators were normalized to process them under a unified framework and make the evaluation results comparable. In this study, the United Nations Development Program's (UNDP) Human Development Index normalization method was used to process data and obtain indicators free of units and scales [41]. Through normalization, each indicator is assigned a value between 0 and 1. Based on the relationship between indicators and the corresponding object of evaluation, different normalization methods were selected. First, indicators were classified as having a positive or negative relationship with the corresponding sub-index, including exposure, sensitivity, and adaptive capacity. In the assessment of exposure and adaptive capacity, all indicators were positive and enhanced the exposure and adaptive capacity by increasing the indicator value; the normalization was processed using Equation (3). However, the sensitivity of several indicators can decrease as the increasing indicator value increases; we used Equation (4) for normalization. Other positive indicators were also processed using Equation (4).
where x i is the x i value of the indicator i after normalization, and x max and x min are the maximum and minimum values of the corresponding indicator, respectively.
where V sub is the sub-index of vulnerability, including exposure, sensitivity, and adaptive capacity, x i is the value of indicator i, and w i is its weight.
After the development of indicator sets for each component of vulnerability, the AHP method, which is a multi-criteria decision-making (MCDM) technique initially introduced by Saaty [42], was used to assign weights to each indicator. AHP works well with GIS and was the main method for indicator integration in this study. Although AHP is a relatively subjective method, it can handle different insights from different perspectives, and allows decision-makers to change a complex problem into a hierarchical structure by identifying relationships between indicators. Therefore, AHP-GIS is frequently used in flood risk and vulnerability research [43]. Other methods to determine index weights are available, such as the entropy weight method, which is an objective method because the weight is determined by the corresponding indicator data. However, flood experience cannot be considered, and the relationships between indicators are ignored. Therefore, the AHP is considered a reasonable method for this study.
In this study, AHP was used as a weighted linear summation method in which the weight of an indicator was gained through pairwise comparison of the variable in a level of decision-making, and the sub-index was obtained from Equation (5). The final selected indicators and their corresponding weights are shown in Table 1. Finally, the WLC method was used to construct the sub-index and obtain the integrated vulnerability index, as shown in Equation (6): where V, E, S, and AC are the vulnerability, exposure, sensitivity, and adaptive capacity, respectively; w E , w s , and w AC are the corresponding weights, with values 0.5, 0.25, and 0.25 in this study, respectively [44].

Flood Modeling
Flooding maps, which were obtained from the results of urban flood modeling, display areas affected by different types of urban floods. This study extracted the maximum modeled floodwater depth at each location as the hazard factor; this was, then, used to identify the population and infrastructure in the hazardous areas to identify the worst-case flood scenario.
Pluvial and fluvial flooding have different causes. In general, urban pluvial flooding is caused by rainfall exceeding the drainage capacity of the city and always occurs in low-lying areas (e.g., community 13 of Lishui City), whereas fluvial flooding occurs due to river levee overtopping. In Lishui City, pluvial flooding is widely distributed in a discontinuous fashion across many communities (Figure 3a-g), including the urban center and suburban areas, and it continues to extend with increasing rainfall. In contrast, the 2014 50-year fluvial flood was continuously distributed in the low-lying area along the river (Figure 3h). As pluvial and fluvial flooding have different distribution characteristics, combining rainfall with river overtopping (i.e., compound flooding; Figures 4 and 5) has a significant effect on the spatial extent of the inundation area.   In Figure 6, the modeling output shows that pluvial flooding following a one-year rainfall had minimal effects on communities; only one community was slightly affected. However, as the rainfall intensity reached that of five year rainfall, the number of affected communities increased to 22, although the water depth remained shallow ( Figure 6a). As the rainfall increased, the water depth and the range of pluvial flooding gradually increased. For 50-year rainfall, almost all communities were threatened by urban pluvial flooding. However, the increase in flooding area was relatively slow compared to that in the flood-affected communities. During the 2014 50-year fluvial flood, nine communities were affected, and are represented by a horizontal line in Figure 6 because only one pluvial flood scenario was designed in this study.   For compound flooding, we compared two methods: GIS and flood modeling. Flooding caused by a combination of rainfall and river overtopping represented the most severe urban flood scenario. Inundation occurred in almost every community following the 20-year rainfall combined with fluvial flooding. Overall, the results of pluvial and compound flooding displayed the same trends under the influence of serial rainfall.
Given the combined effects of rainfall and fluvial flooding, compound floods are more serious than the pluvial flood. Nine communities were affected by more serious flooding than pluvial flooding, with fluvial flooding being the dominant contributor. However, compound and pluvial flooding results were similar when rainfall was 20-year rainfall or less (Figure 6b), with the curve from the modeling method located above that of the GIS method. Thus, the modeling of compound flooding shows a greater flood impact than the GIS method.  In Figure 6, the modeling output shows that pluvial flooding following a one-year rainfall had minimal effects on communities; only one community was slightly affected. However, as the rainfall intensity reached that of five year rainfall, the number of affected communities increased to 22, although the water depth remained shallow ( Figure 6a). As the rainfall increased, the water depth and the range of pluvial flooding gradually increased. For 50-year rainfall, almost all communities were threatened by urban pluvial flooding. However, the increase in flooding area was relatively slow compared to that in the flood-affected communities. During the 2014 50-year fluvial flood, nine communities were affected, and are represented by a horizontal line in Figure 6 because only one pluvial flood scenario was designed in this study.
For compound flooding, we compared two methods: GIS and flood modeling. Flooding caused by a combination of rainfall and river overtopping represented the most severe urban flood scenario. Inundation occurred in almost every community following the 20-year rainfall combined with fluvial flooding. Overall, the results of pluvial and compound flooding displayed the same trends under the influence of serial rainfall. Given the combined effects of rainfall and fluvial flooding, compound floods are more serious than the pluvial flood. Nine communities were affected by more serious flooding than pluvial flooding, with fluvial flooding being the dominant contributor. However, compound and pluvial flooding results were similar when rainfall was 20-year rainfall or less (Figure 6b), with the curve from the modeling method located above that of the GIS method. Thus, the modeling of compound flooding shows a greater flood impact than the GIS method.

Exposure
Based on the distribution results, exposure was divided into five levels: very low, low, medium, high, and very high. Community exposures to different types of floods are shown in Tables 2-4.   Comparing communities with different exposure levels, exposure shows obvious differences among different flood types. For pluvial flooding (Table 2), when the rainfall intensity was that of 20-year rainfall or less, all communities were at very low and low exposure levels. When rainfall intensity increased to that of 30-year rainfall, community 13 shifted to medium exposure. As rainfall intensity increased to that of 100-year rainfall, three communities were at a medium exposure level. For fluvial flooding, there were three communities at a medium exposure level. For compound flooding, under the combined effects of the two flood factors, the number of communities at medium or higher exposure levels increased significantly. When rainfall intensity was that of 20-year rainfall or less in the compound floods, the number of communities with medium exposure was similar to that in the fluvial flood scenario. More than 20% of the communities were at medium or higher exposure levels for the worst compound flood scenario. There were also differences among the compound floods obtained using the different methods. In particular, for rainfall intensity of 50-or 100-year rainfall, communities with high and very high exposure appeared in the results of flood modeling.
Considering the spatial locations of communities, the spatial distributions of exposure also differed for different types of floods. Communities with medium or higher exposure levels (communities 9, 13, 23) appeared in urban centers when the rainfall intensity was that of 30-year rainfall and in the suburban areas when rainfall intensity reached that of 100-year rainfall. Overall, only a few communities had medium or high exposure as a result of pluvial flooding, and these were mainly distributed in local low-lying areas. In terms of fluvial flooding, the three communities with medium exposure were distributed along the river. For compound flooding, communities with medium or higher exposure appeared in the city center and suburban areas simultaneously.

Sensitivity and Adaptive Capacity
The results of sensitivity and adaptive capacity were divided into four levels using the natural break method (Jenks) [45], as shown in Figure 7. For sensitivity, six communities were at a high level, seven communities were at medium levels, and the remaining were at low and very low levels. All high-sensitivity communities were located in the central urban area, except for community 15, whereas the communities with low and very low sensitivity were located on the edges of the urban area. For adaptive capacity, five communities with high adaptive capacity formed two high-value areas within the city; these communities had easy access to health services and shelters. Communities with low and very low adaptive capacity were predominantly located in the suburban and riverine areas, far from shelters. formed two high-value areas within the city; these communities had easy access to health services and shelters. Communities with low and very low adaptive capacity were predominantly located in the suburban and riverine areas, far from shelters.

Vulnerability Assessment
Based on the assessment of exposure, sensitivity, and adaptive capacity, the vulnerability was integrated using WLC. The vulnerability was also divided into five levels, namely very low, low, medium, high, and very high. The results are shown in Tables 5-7. Community vulnerability varied for different urban flood scenarios. For pluvial flooding (Table 5), when rainfall intensity was that of 20-year rainfall or less, all communities had low or very low vulnerability levels. We found that 20% of communities had a medium vulnerability in the worst compound flood scenario. No communities had high or very high vulnerability for the pluvial flood. For fluvial flooding (Table 5), four communities had a medium vulnerability, and there were no high-vulnerability communities. For compounding flood according to the GIS method (Table 6), when rainfall intensity was that of 20-year rainfall or less, the number of communities with medium vulnerability was similar to that for fluvial flooding. The number of medium-vulnerability communities significantly exceeded that for fluvial or pluvial flooding when the rainfall intensity reached that of 50-or 100-year rainfall. Compound flood vulnerability with flood modeling showed an increase in the number of communities with high or very high vulnerability levels.

Vulnerability Assessment
Based on the assessment of exposure, sensitivity, and adaptive capacity, the vulnerability was integrated using WLC. The vulnerability was also divided into five levels, namely very low, low, medium, high, and very high. The results are shown in Tables 5-7.   Community vulnerability varied for different urban flood scenarios. For pluvial flooding (Table 5), when rainfall intensity was that of 20-year rainfall or less, all communities had low or very low vulnerability levels. We found that 20% of communities had a medium vulnerability in the worst compound flood scenario. No communities had high or very high vulnerability for the pluvial flood. For fluvial flooding (Table 5), four communities had a medium vulnerability, and there were no high-vulnerability communities. For compounding flood according to the GIS method (Table 6), when rainfall intensity was that of 20-year rainfall or less, the number of communities with medium vulnerability was similar to that for fluvial flooding. The number of medium-vulnerability communities significantly exceeded that for fluvial or pluvial flooding when the rainfall intensity reached that of 50or 100-year rainfall. Compound flood vulnerability with flood modeling showed an increase in the number of communities with high or very high vulnerability levels.
According to the spatial distribution of vulnerability, for pluvial flooding, initially, communities with medium vulnerability were almost distributed in suburban areas, except for community 13 (city center). For fluvial flooding, all medium-vulnerability communities were distributed along the river. For compound flooding according to the GIS method, the distribution was the same as that for the superposition of pluvial and fluvial flooding. Most communities with medium or higher vulnerability were located in suburban areas or along the river. The vulnerabilities to compound flooding using the modeling method had a similar distribution to those using the GIS method.

Discussion
Assessing community vulnerability to urban flooding, which is important for flood hazard mitigation planning, requires the identification of vulnerable communities and factors that can increase the vulnerability. The results of this study suggest that the vulnerabilities of different communities within the same city vary considerably owing to the physical and socioeconomic gaps; the different vulnerability levels can be distinguished by analyzing the results of flood modeling and vulnerability assessment. Flood and vulnerability maps display the exposure and vulnerability of communities to different types of floods, can help city management departments and decision-makers in developing flood mitigation and adaptation strategies on a community scale, and can reduce the adverse effects of urban flooding in an efficient and target-specific manner.
The results of this study have implications for the methods and theory of vulnerability studies. First, a coupled flood model combined with GIS techniques was applied to map different types of floods and the relationships with physical and social indicators; from this, it was possible to map exposure, sensitivity, and adaptive capacity, and then to create a composite vulnerability index using the AHP and WLC methods. In particular, the assessment of exposure is linked to the flood hazard, and the results were consistent with reality; that is, the flood was more severe when a large population and more infrastructure were affected. Therefore, the exposure and vulnerability under different flood scenarios, which can help in flood management according to actual local conditions, are critical. This study addressed a gap in urban flood vulnerability research concerning different types of floods, including compound flooding induced by rainfall and river flooding. Currently, most flood vulnerability studies are focused on floods caused by single factors, such as river overtopping or storm surges. There are numerous vulnerability studies of river flooding with even the term "flood" commonly being defined as a river-related phenomenon. This study substantially differs from others because it simultaneously considered pluvial and fluvial floods (induced by rainfall and river overtopping, respectively) and analyzed the effects of these different flood types on vulnerabilities of urban communities by reexamining the theoretical basis of the "place-based" vulnerability and GIS methods [21,22,46,47].
Community vulnerability to different types of floods has different characteristics. The results of exposure show that fluvial flooding has a more substantial effect on riverine communities compared to pluvial flooding induced by 30-year or less rainfall. Therefore, fluvial flooding is frequently studied in vulnerability research. However, pluvial flooding can have an effect comparable to that of fluvial flooding when the rainfall intensity reaches that of 50-or 100-year rainfall, which is similar to the findings of existing studies [12]. Pluvial flooding and urban pluvial flooding can have intersections and interactions, such as compound flooding induced by rainfall and river flooding. Compound flooding has the most severe effect on the three types of communities because it has a wider distribution, reflecting the fact that fluvial and pluvial floods can have substantially different impact areas. Moreover, compound flooding results in more medium-and high-vulnerability communities. The vulnerability of communities to flooding can be underestimated when compound flooding is not considered. Multiple or compound hazards are increasing owing to climate change [18,20], and the assessment of compound flood vulnerability is necessary to aid communities to cope with possible future flood disasters.
The results of this study suggest that urbanization contributes to increased flood vulnerability. In the study area, increased infrastructure, such as hospitals and shelters in the urban center area, can enhance the adaptive capacity of a community, as nearby communities can easily access these services. Communities in suburban areas are less likely to be adaptive owing to the greater distance from infrastructure facilities. Thus, infrastructure development should meet the needs of urban communities; otherwise, communities within newly developed areas may become more susceptible to disasters. Central urban areas can be low-value aggregation areas of urban social vulnerability, whereas surrounding areas have a high-value aggregation of social vulnerability [48]. Moreover, this study provides suggestions for more effective urban flood management and urban planning. Our spatial vulnerability assessment approach can identify areas and communities vulnerable to urban flooding. Different response measures should be implemented according to the level of vulnerability and flood type. For example, dams along the river in communities 11, 14, 15, and 24 need to be strengthened, low-impact development measures should be carried out in community 13, and more medical facilities and shelters are needed in suburban areas.
The method used in this study is limited because of its high dependency on the flood modeling and the limited availability of socioeconomic data. Although the coupled model is a well-established technique that is widely used in flood risk and vulnerability studies, the calibration and validation are still characterized by uncertainty. Therefore, more validation should be done to verify the model further. In flood modeling, the influence of river water level on urban drainage is important; however, in this study, only one fluvial flood scenario was used (i.e., the worst fluvial flood scenario). This was chosen because flood modeling and flood hazards were not the focus of this research; however, the complexity and influencing factor combinations of more compound flood scenarios should be studied in the future. Data are a crucial constraint in vulnerability studies. For this reason, we selected POI data, whose collection requires manual data evaluation. However, vulnerability research requires current and real-life data, especially for communities with small amounts of demographic and socioeconomic statistics, and POI data may have a timeliness issue; for example, errors caused by POI data are not updated on time.

Conclusions
This study outlined an integrated method for the assessment of community vulnerability to urban flooding based on flood modeling and GIS. Modeling of various urban flood scenarios, including fluvial, pluvial, and compound flooding, was carried out to reveal the level of vulnerability to different types of urban floods. The method was applied to Lishui, Zhejiang, China; however, it could be used in other cities owing to its strong replicability. The flood modeling results show that different types of flood can occupy different areas and have adverse effects: Pluvial urban floods can occur in many parts of a city, fluvial flooding mainly affects areas along rivers, and compound flooding combines rainfall with river overtopping; it has a wider range of effects and can be more serious than the simple accumulation of fluvial and pluvial flooding. Exposure assessment based on flood modeling quantified the affected population and infrastructure. Sensitivity and adaptive capacity reflect the effect of a community's socioeconomic status on its flood vulnerability. The spatial distributions of the two were found to be similar; that is, urban centers tend to have higher sensitivity and adaptive capacity, and the opposite is true for suburban areas. Flood vulnerability assessment is a composite of multiple effects, including disaster, socioeconomic factors, and vulnerability, indicating that different flood vulnerabilities have substantially varied effects. Compound flood vulnerability is higher than single-factor flood vulnerability, especially in regions where multiple factors are combined.
The results of this study can guide urban managers and policymakers on community flood management. These results revisit the gaps in current flood measures and provide guidance to reduce vulnerability. Further validation of the assessment results is required in future studies, including those on flood simulation and vulnerability.