Next Article in Journal
The Assessment of Indoor Formaldehyde and Bioaerosol Removal by Using Negative Discharge Electrostatic Air Cleaners
Previous Article in Journal
Analysis of the Effect of Human-Machine Co-Driving Vehicle on Pedestrian Crossing Speed at Uncontrolled Mid-Block Road Sections: A VR-Based Case Study
Previous Article in Special Issue
The Impact of Sustainability Goals on Productivity Growth: The Moderating Role of Global Warming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urban Inundation under Different Rainstorm Scenarios in Lin’an City, China

1
Institute of Remote Sensing and Geosciences, Hangzhou Normal University, Hangzhou 311121, China
2
Zhejiang Provincial Key Laboratory of Urban Wetlands and Regional Change, Hangzhou 311121, China
3
Faculty of Geo-Information Science and Earth Observation (ITC), University of Twente, 7500 AE Enschede, The Netherlands
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(12), 7210; https://doi.org/10.3390/ijerph19127210
Submission received: 9 April 2022 / Revised: 6 June 2022 / Accepted: 10 June 2022 / Published: 12 June 2022
(This article belongs to the Special Issue Global Warming & Water)

Abstract

:
Under the circumstances of global warming and rapid urbanization, damage caused by urban inundation are becoming increasingly severe, attracting the attention of both researchers and governors. The accurate simulation of urban inundation is essential for the prevention of inundation hazards. In this study, a 1D pipe network and a 2D urban inundation coupling model constructed by InfoWorks ICM was used to simulate the inundation conditions in the typical urbanized area in the north of Lin’an. Two historical rainfall events in 2020 were utilized to verify the modeling results. The spatial–temporal variation and the causes of urban inundation under different designed rainfalls were studied. The results were as follows: (1) The constructed model had a good simulation accuracy, the Nash–Sutcliffe efficiency coefficient was higher than 0.82, R2 was higher than 0.87, and the relative error was ±20%. (2) The simulation results of different designed rainfall scenarios indicated that the maximum inundation depth and inundation extent increased with the increase in the return period, rainfall peak position coefficient, and rainfall duration. According to the analysis results, the urban inundation in Lin’an is mainly affected by topography, drainage network (spatial distribution and pipe diameter), and rainfall patterns. The results are supposed to provide technical support and a decision-making reference for the urban management department of Lin’an to design inundation prevention measures.

1. Introduction

The spatial and temporal distributions, intensities, and frequencies of rainfall are changing because of global warming [1]. Rainfall is becoming more uneven and intense, and these variations in rainfall have aggravated the occurrence of urban inundation, soil erosion, and other hazards [2,3]. Meanwhile, the urban expansion in China is accelerating, increasing the impervious surface of the city, changing the urban hydrological process and increasing the possibility of urban inundation [4,5,6]. The preferential flow process that occurs in cases of heavy rainstorms cannot be neglected, in which soil material and soil moisture have an important influence [7,8,9]. Urban inundation is a common issue that can occasionally become destructive and cause considerable losses in human life and economy [10,11]. For example, affected by the typhoon in July 2021, Yuyao, Zhoushan, and many other cities in China were flooded, affecting a total of 2.711 million people and causing direct economic losses of CNY 3.35 billion [12]. Therefore, in recent years, urban inundation has attracted increasing social attention [13,14]. The accurate simulation of urban inundation scenarios is important for the prevention of inundation hazards and the reduction in hazard losses.
Today, urban inundation models are widely used for studying urban inundation [15]. Scholars have developed many urban inundation models, including SSCM (Urban Storm Pipeline Calculation Model) [16], CSYJM (Urban Storm Runoff Model) [17], SWMM (Storm Water Management Model) [18], InfoWorks ICM (InfoWorks Integrated Catchment Management) [19,20], and MIKE series [21,22]. Models such as SSCM and CSYJM lack good pre- and post-processing programs, and most of them are based on applying research in a specific region [23]. These models are usually limited to the internal use of the development teams and are not widely used owing to a lack of promotion [24,25]. On the other hand, urban inundation models such as SWMM, the MIKE series, and InfoWorks ICM have a great number of users. Yu et al. constructed an urban inundation model of Jinan based on SWMM and used it to analyze the flood control measures of Jinan [26]. Patro et al. constructed a one-dimensional (1D) river model and a two-dimensional (2D) hydrodynamic model for the Delta region of the Mahanadi River Basin in India using MIKE 11 and MIKE 21, respectively [27]. Then they used MIKE FLOOD coupling to simulate the inundation extents and depths of the study area under different grades of flood. Sidek et al. used the InfoWorks ICM to construct a 1D river and 2D surface inundation coupling model in the Damansara catchment, Malaysia [28]. The model was verified by monitoring data, and the results showed that the model has a good accuracy (R2 > 0.8). These efforts by previous scholars proved the suitability of using inundation models to study urban inundation issues.
SWMM is one of the most widely used models because its code is open source [18]. However, as a 1D model, SWMM cannot suitably model surface flooding [29]. Commercial software such as the MIKE series models and InfoWorks ICM have good pre- and post-processing abilities to model surface flooding [30]. However, the MIKE series models (such as MIKE 11/MIKE 21) do not allow complete coupling. The MIKE FLOOD model is needed to couple the 1D model (MIKE 11/MIKE Urban) and the 2D model (MIKE 21) [31]. As opposed to the MIKE series models, InfoWorks ICM is the first model in the world to combine a 1D hydraulic model of urban drainage networks and channels with a 2D flood inundation model of the urban catchment in a single simulation engine [32]. Furthermore, the operation of InfoWorks ICM is easier than that of the MIKE series models, and the simulation accuracy is better than that of SWMM [25]. In this context, InfoWorks ICM is now widely used in urban flood management and assessment [33,34].
In this study, a hydrological model was constructed by InfoWorks ICM to detect the spatial patterns and affecting rules of urban inundation on the northern part of Lin’an City, Hangzhou. The main objectives of this study are: (1) to analyze the 2D surface inundation in the study area through the constructed fine hydrological model; (2) to analyze the spatial-temporal variation and the causes of urban inundation in the study area under designed rainfalls. This study is supposed to provide scientific support for making urban inundation prevention and control plans in Lin’an City.

2. Materials and Methods

2.1. Study Area

Lin’an City, located in the northwest of Zhejiang Province, China, is a county-level city governed by Hangzhou. The city experiences a mid-subtropical monsoon climate, which is warm and humid with abundant rainfall. This area is easily affected by monsoon and typhoon weather, with an average annual precipitation of 1613.9 mm. Continuous rainfall and short-term heavy rainfall occur frequently in this area. In recent years, Lin’an has experienced rainstorms many times, resulting in inundation in large parts of the urban area. These rainstorms have caused traffic jams and paralysis, which resulted in huge losses to people’s property and social economy. In June 2015, there was a sudden heavy rainfall in Lin’an, causing the collapse of 91 houses and a direct economic loss of CNY 180 million [35]. Similarly, in August 2019, a short-term rainstorm in Lin’an caused four deaths and a direct economic loss of CNY 169 million [36].
The selected study area is located in the north of the main urban area of Lin’an (30°14′55″ N–30°15′48″ N, 119°42′7″ E–119°43′27″ E) (Figure 1). It is a relatively independent catchment unit formed by the intersection of Agricultural and Forestry Road (AFR), Maxi River, Baini Road (BNR), and Shimen Ridge, with a vertical range of 31.34–135.14 m and an area of 1.89 km2. The study area is a rapid urbanization region, with construction land accounting for 87.83% of the total area. The population and residential buildings are relatively concentrated, and the roads (such as Maxi Road and Xishu Street) have been inundated many times, threating people’s lives and urban transportation. Therefore, it is a suitable site to simulate urban inundation and analyze its causes.

2.2. Data Collection and Pre-Processing

2.2.1. Geographic Data

The geographic data applied in this study mainly included: (1) A digital elevation model (DEM) provided by the surveying and mapping department of Lin’an with a 2-m spatial resolution. The elevation range of the DEM is from 31.34 to 135.14 m. (2) A UAV image with 0.5 m spatial resolution, provided by the surveying and mapping department of Lin’an. The image was captured in December 2019. (3) Land-use data consisting of 21 land-use types (including residential land, industrial land, construction land, road, rivers, and other types). These data came from the surveying and mapping department of Lin’an. According to previous researches, the land-use types were classified into buildings, roads, green space, water, and others [25,37]. (4) Two historical rainfall events (29 May 2020 and 2 July 2020) were selected for modeling and derived from the meteorological department of Lin’an. (5) Drainage system data provided by the urban management department of Lin’an were used, which contained pipe diameter, pipe materials, and other information (Table 1).

2.2.2. Designed Rainfall Events

In order to study the inundation in the study area under different return periods, rainfall peak position coefficients, and rainfall duration, a variety of rainfall scenarios were set up for the analysis. According to the Intensity–Duration–Frequency (IDF) formula of a Lin’an rainstorm proposed by the Zhejiang Provincial Department of Housing and Urban–Rural Development, the rainfall intensity of designed rainfall is calculated with the following formula [38]:
q = 2763.132 × ( 1 + 0.399 l g P ) ( t + 10.870 ) 0.753
where q is the designed rainfall intensity (L/(s·hm2)); t is the rainfall duration (min); and P is the return period (a). Rainfall scenarios of 60, 120, and 180 min and 1-, 2-, 5-, and 10-year return periods were used to simulate the urban inundation depths and extents. Related studies have shown that the rainfall peak position coefficient, r, in most areas ranges from 0.3 to 0.5 [39]. Thus, rainfall peak position coefficients were set to 0.3, 0.4, and 0.5. In this study, a total of 36 rainfall scenarios were designed for analysis. Figure 2 presents a section of the hydrograph of the designed rainfall events (taking the rainfall duration of 120 min as an example).

2.3. Overall Workflow

This study included three main parts: the construction of InfoWorks ICM model, model validation, and simulation of different designed rainfalls (Figure 3). First, we input the processed geographic data (DEM data, drainage system data, land-use data, and historical rainfalls) into the InfoWorks ICM model. Then, the model was verified by the observed flow data and investigation results of inundation depths. Finally, we used the constructed model to simulate inundation depths and extents under different designed rainfall scenarios. According to the results, the spatial-temporal variation and causes of urban inundation under historical rainfalls and designed rainfalls were analyzed.

2.4. Model Description

As one of the best urban inundation models, InfoWorks ICM can achieve integrated 1D and 2D hydrological–hydraulic modeling. The model has many advantages. For example, it can shorten the modeling time with its parallel computation based on GPUs, especially 2D hydrodynamic modeling [25]. In addition, data in ArcGIS, AutoCAD, and Excel formats can be imported into InfoWorks ICM, which makes data processing easier [40]. Furthermore, the model has many modules, including hydrological, 1D pipe network flow hydraulic, 2D urban inundation, and real-time control modules [41]. The main modules involved in this study were the hydrological, 1D pipe network flow hydraulic, and 2D urban inundation modules. The hydrological module serves to calculate urban surface rainfall runoff according to the spatial division of the sub-catchment area and the surface composition of different runoff characteristics. Its main calculation modules include a runoff model and a routine model. In the hydraulic calculation of 1D pipe network, the flow velocity and water depth of the pipeline are simulated by solving the Saint-Venant equations (Equations (2) and (3)). The 2D urban inundation process is implemented by solving shallow water equations with the 2D finite volume method (Equations (4)–(6)):
A t + Q x = 0 ,
Q t + x ( Q 2 A ) + g A ( cos θ h x S 0 + Q | Q | K 2 ) = 0 ,
h t + ( h u ) x + ( h v ) y = q 1 D ,
( h u ) t + x ( h u 2 + g h 2 2 ) + ( h u v ) y = S 0 , x S f , x + q 1 D u 1 D ,
( h v ) t + y ( h v 2 + g h 2 2 ) + ( h u v ) x = S 0 , y S f , y + q 1 D v 1 D
where Q is the flow (m3/s); A is the cross-sectional area of the pipeline (m2); t is the time (s); x is the length along the x direction (m); h is the water depth (m); g is the gravity acceleration (m/s2); θ is the horizontal angle (degree); K is the conveyance and is calculated using the Colebrook–White formula or the Manning formula; S 0 is the bed slope; u and v are the velocity components in the x and y directions (m/s), respectively; S 0 , x and S 0 , y are the bed slope components in the x and y direction (m2/s2), respectively; S f , x and S f , y are the friction components in the x and y directions (m2/s2), respectively; q 1 D is the outflow per unit area (m3/s); u 1 D and v 1 D are the velocity components of q 1 D in the x and y directions (m/s), respectively.

2.5. Model Setup

2.5.1. Setup of Hydrological Model

The process of constructing the hydrological model mainly included: (1) the division of the sub-catchment area and (2) the determination of the runoff model and routine model. First, the sub-catchment areas were divided by the Thiessen polygon method based on the distribution of nodes. Then, the sub-catchment areas were modified according to the wall structure, building, slope, and aspect. The nearby sub-catchment without nodes flowed into the adjacent sub-catchment according to the aspect. In this study, a total of 247 sub-catchment areas were obtained. The nature of the underlying surface can affect the processes of rain water infiltration and confluence, so it was necessary to extract the runoff surface of the sub-catchment according to the land-use types. The Area Take Off (ATO) tool in InfoWorks ICM was used to extract the runoff surface type and its proportion in each sub-catchment. The parameters of each runoff-producing surface are shown in Table 2.

2.5.2. Setup of 1D Pipe Network Hydraulic Model

The construction of the 1D pipe network hydraulic model involved the construction of the drainage system, which was needed to deal with the original drainage system data. The main process was deletion of the rain grate and its connecting pipeline. Because the rain grate was connected to the drainage system through a thin and short pipeline, its deletion had little impact on the final simulation results. Finally, there were 245 pipelines, 244 manholes, and 5 outfalls in the study area, and the pipe diameter ranged between DN200 and DN2000. The Manning coefficient of the pipeline is an important input parameter of InfoWorks ICM, and different materials have different values. The pipelines in the study area were mainly made of concrete and plastic. Referring to the SWMM user manual [42], it was determined that the Manning coefficient of the concrete pipeline was 0.013, and that of the plastic pipeline was 0.012.

2.5.3. 1D–2D Model Coupling

In InfoWorks ICM, the 1D model was mainly used to simulate the drainage network flow, and determine the overflow nodes and overflow volume. The 2D model was used to simulate the diffusion of overflow generated by the drainage network on the complex surface. Through this process, the inundation depths and extents were obtained. The first step was building a 2D urban inundation model for the ground elevation model. The main process was extracting the elevation points from the DEM data, and then importing them into InfoWorks ICM to establish the TIN (Triangulated Irregular Network) model. Then, the 2D interval mesh was divided in the study area. The minimum mesh resolution was 25 m2, and the maximum was 100 m2. In this process, buildings were regarded as non-flooding zones and were not meshed [43]. In addition, considering that the road is generally lower than the sidewalks on both sides, its elevation was reduced by 15 cm and meshed. Finally, we set the flood type of manholes as “2D”, which achieved 1D and 2D model coupling.

2.6. Model Validation

In order to ensure the accuracy and applicability of the model, the observed data were used to verify the model. The two historical rainfalls on 29 May 2020 and 2 July 2020 were selected as the rainfall input data. The model was verified by comparing the observed discharge of flow monitoring points and the simulation results of the model. The coefficient of determination (R2), root mean square error (RMSE), and Nash–Sutcliffe efficiency coefficient (NSE) were used to evaluate the simulation results. The formula of NSE is as follows:
N S E = 1 t = 1 T ( Q o t Q m t ) 2 t = 1 T ( Q o t Q ¯ o ) 2
where Q o t is the observed flow at t time (m3/s), Q m t is the simulated flow at t time (m3/s), and Q ¯ o is the average value of observed flow (m3/s).

3. Results

3.1. Simulation Results of Historical Rainfalls

Figure 4 shows the validation results of the two historical rainfalls. Table 3 summarizes the simulated and observed flows in the validation processes. The NSE of the two historical rainfalls was more than 0.82, the R2 was above 0.87, and the relative error was ±20%, as shown in Table 3. In addition, after investigating the inundation depths and extents during the historical rainfalls in the study area, some areas were detected in inundation. Table 4 presents the recorded and simulated maximum inundation depths during the 2 July 2020 rainfall event. The table shows that the simulated maximum depth of the model was similar to the recorded depth. After calculation, the relative errors of MXR, XSS, and SLS were 15.5%, –17%, and –13%, respectively. Although there were some differences between the simulation results and the recorded results, the differences were in an acceptable range. The constructed model was able to simulate the inundation process in Lin’an and could be used for the following studies.

3.2. Simulation Results of Designed Rainfalls

3.2.1. Different Return Periods

The designed rainfalls of different return periods (1-, 2-, 5-, and 10-year), different rainfall peak position coefficients (r = 0.3/0.4/0.5), and different rainfall durations (t = 60 min/120 min/180 min) were input into the constructed model for simulation. After simulation, statistical analyses of the inundation depths and extents were accomplished. Statistics included the number of overflow nodes, maximum overflow volume of nodes, maximum inundation depth, and inundation extent under different designed rainfall scenarios. Table 5 and Figure 5 present the inundation depths and extents of the study area under different designed rainfall scenarios (t = 120 min).
Figure 5 shows the spatial distribution of inundation under different designed rainfall scenarios. The results indicated that the inundation of the study area became more serious as the return periods increased, and the inundation extents in the south of the study area were increasing. Table 5 presents the inundation depths and extents under different designed rainfall scenarios. When the rainfall peak position coefficient was 0.4, the number of overflow nodes was 13, 18, 31, and 37, respectively, in the return periods of 1, 2, 5, and 10 years. Similarly, the maximum overflow volumes of the nodes were 0.5 m3/s, 0.6 m3/s, 0.8 m3/s, and 0.9 m3/s, respectively. The maximum inundation depths were 0.847 m, 0.978 m, 1.126 m, and 1.222 m, respectively. In correspondence, the inundation extents were 84,453 m2, 103,769 m2, 148,234 m2, and 195,118 m2. Thus, with the same rainfall peak position coefficient, the number of overflow nodes, maximum overflow volumes of nodes, maximum inundation depths, and inundation extents increased as the return periods increased. Furthermore, through comparing the contribution areas of inundation depths under different return periods, the area with a depth of less than 0.2 m increased the most. The statistical results showed that the inundation area with a depth of less than 0.2 m accounted for about 90% of the increased area from 1 year to 10 years.

3.2.2. Different Rainfall Peak Position Coefficients

Table 5 shows that when the return period was 5 years, the numbers of overflow nodes under different rainfall peak position coefficients were 28, 31, and 32, respectively. The maximum overflow volumes of the nodes were 0.7 m3/s, 0.8 m3/s, and 0.8 m3/s. The maximum inundation depths were 1.096 m, 1.126 m, and 1.139 m, and the inundation extents were 137,329 m2, 148,234 m2, and 158,139 m2. Different return periods had similar results. The results showed that the number of overflow nodes, the maximum overflow volume of nodes, maximum inundation depths, and inundation extents increased with the increase in rainfall peak position coefficients in the same return period. By comparing the inundation conditions under different scenarios (see Table 5), we found that the inundation area with a depth of less than 0.8 m increased gradually, while the area with a depth of over 0.8 m remained almost unchanged. The lag of the rainfall peak increased the inundation area, and 99% of the increased area came from the area with a depth of less than 0.8 m. The area with a depth of less than 0.2 m increased most, accounting for 87% of the increased area (Figure 6).

3.2.3. Different Rainfall Durations

Table 6 presents the inundation depths and extents under different rainfall durations when the rainfall peak position coefficient was 0.4. With the increase in rainfall duration, the cumulative rainfall, maximum inundation depth, and inundation extents all increased in the same return period. However, the maximum inundation depth increased little as the rainfall duration increased. When the rainfall duration increased from 60 to 180 min, the maximum inundation depth increased by 0.06 m on average. Compared with the maximum inundation depth, the cumulative rainfall and inundation extent changed more significantly. Taking the 10-year rainfall events as an example, when the rainfall duration increased from 60 to 180 min, the cumulative rainfall increased by 42.20%, from 56.4 to 80.2 mm. Meanwhile the inundation area increasing by 10.92%, from 181,229 to 201,027 m2. These results indicated that the rainstorm with a longer rainfall duration brought more rainfall, which resulted in a larger inundation extent (Figure 7).

4. Discussion

4.1. Urban Inundation under Historical Rainfall

In this study, an InfoWorks ICM model was employed to build a 1D pipe network and 2D urban inundation coupling model. The values of the NSE, R2, and RMSE in the validation results of the coupled model were similar to the results of Zhang et al. [41] and Liew et al. [44]. This showed that the 1D–2D coupled model constructed in this study was reliable and able to accurately simulate the inundation depths and extents in the study area. Figure 8 shows that in the rainfall event on 2 July 2020, there were three obvious inundation areas in the study area, namely, MXR, XSS, and SLS. MXR was the most serious inundation area, with a maximum inundation depth of 20 cm. MXR is located in the west of the study area, whose topography in the south is about 80 cm higher than that in the north. Moreover, there is almost no pipeline distribution and only one manhole in the south. During the rainstorm, this manhole could not be discharged in time after gathering too much rainwater and thus caused overflowing. Rainwater flowed from south to north to form an inundation in low-lying areas, according to the topography (Figure 9). As a result, the characteristics of inundation in this area were deep in the north and shallow in the south, and the inundation area was wide. The inundation analysis results showed that the two important factors affecting the formation of inundation were topography and pipe network distribution. In order to avoid serious inundation in MXR, the low-lying areas should be filled. Moreover, additional drainage outlets should be added to speed up the discharge of rainwater. The simulation results also showed that the InfoWorks ICM model could simulate inundation depths and extents according to the local micro-topography of the study area.

4.2. Urban Inundation Projections under Different Designed Rainfalls

The inundations in the study area showed a similar pattern under different designed rainfall scenarios. With an increase in the return period, rainfall peak position coefficient, and rainfall duration, the maximum inundation depth and inundation extent increased. As a result, the urban inundation became more serious. This result was consistent with previous studies of urban areas that considered the influence of different designed rainfalls on urban inundation [45,46,47]. From the perspective of spatial distribution, the inundation areas are mainly concentrated in the south of the study area (Figure 10) because the topography in the south of the study area is lower than that in the north. During the rainstorm, most rainwater flowed southward and into the drainage network in the south. Compared with the drainage network in the north, the drainage network in the south was under higher drainage pressure. Therefore, the southern areas were more likely to be inundated than the northern ones. This was similar to the spatial distributions of inundation found in the Chebei Watershed [48], which were mainly concentrated in the lower terrain areas.
In order to further analyze the change in inundation of the study area, four areas (Figure 10A–D) were selected for detailed analysis. XJSR (Figure 10A) is located in the middle of the study area with an elevation ranging from 36.4 m to 38.8 m. The diameter of the drainage network in this area is generally 2000 mm. Under the designed rainfall scenarios, there was no inundation in this area. This result showed that the areas with higher terrain and larger pipe diameter are not easily inundated. NLS (Figure 10B) is located in the northwest of the study area and is higher than XJSR. In contrast to XJSR, the diameter of the drainage network in NLS is 600 mm. This area was inundated under the 1-year rainfall scenario. With the increase in return periods, the inundation depth increased. Under the 10-year rainfall scenario, the inundation depth in some areas was over 0.6 m. This happened because, during the rainstorm, the drainage network could not accommodate lots of rainwater and thus overflowed. Although some rainwater flowed to the lower places under the influence of the terrain, most of the rainwater still flowed into the drainage network. This result indicated that areas with higher terrain and smaller pipe diameter can cause inundation. AFR (Figure 10C) is located in the south of the study area and the terrain is flat (the elevation range is 35.3–35.6 m). The diameters of most of the pipes in the drainage network are 1500 and 1800 mm. This area began to be inundated under the 5-year rainfall scenario. The maximum inundation depth was 0.104 m under the 10-year rainfall scenario. Compared with XJSR, the result of the AFR showed that in the case of a similar pipe diameter, inundation was more likely to occur in lower terrain areas. In addition, AFR was inundated only under the 5-year rainfall scenario, which indicated that rainfall pattern was an important factor affecting the inundation. WSS (Figure 10D) is located in the south of the study area, and the terrain is flat. The diameter of the drainage network is small, only 300 mm. Thus, the drainage network could not accommodate much rainwater under the designed rainstorms, forming inundation in the surrounding area. The maximum inundation depth was 0.442 m here. Compared with AFR, the result showed that areas with lower terrain and smaller pipe diameter will be seriously inundated. Overall, urban inundation is mainly affected by topography. However, from a local point of view, urban inundation is not only affected by topography but also affected by the drainage network and rainfall patterns. In summary, inundation not only occurs in areas with low terrain but also occurs in areas where the pipe diameter is too narrow. Moreover, a short-term rainstorm is more likely to inundate cities. Li et al. [49] and Cheng et al. [50] reached a similar conclusion in their study of urban inundation. However, compared with them, we clearly and intuitively showed the impact of different factors on inundation through the presentation of data and details. For the urban management department of Lin’an City, detailed modifications including increasing are small pipe diameters and filling low-lying areas is recommended. Moreover, LID (Low-impact development) techniques, such as permeable pavements, green roofs, and rain barrels, can be implement in the study area to reduce overflows. In addition, we suggest installing flow control devices (gates) in overflowing manholes to reduce inundation and arranging drainage pumps in inundated areas to facilitate the timely removal of accumulated water. Now, we only propose ways to reduce overflows, but how to implement it needs further research with InfoWorks ICM because it could simulate pipe diameter expansion scenarios and different LID scenarios.

4.3. Limitations

Several limitations need to be addressed for a comprehensive understanding in this study. First, the parameters of the 1D–2D coupled model are not universal and only applicable to the current study area. The parameters must be adjusted when the model is extended to other study areas. Parameters can be selected by consulting the InfoWorks ICM user manual and the SWMM user manual. Second, the influence of the rivers was not considered in the modeling process. When there is a serious rainstorm, the excessive water level of the rivers will affect the discharge of rainwater and then affect the inundation depths and extents. Therefore, it is be necessary to collect information such as riverbed elevation and river cross-sectional coordinates to construct the river model in future studies, to improve the accuracy of the model. Third, the land-use data in this study were extracted from historical remote sensing images, which do not reflect the latest land use in the study area and may create uncertainty in the model. Finally, soil moisture will affect the simulation results. The widely used Horton’s model is used in the model construction. However, the model does not account for initial moisture conditions adequately [51]. In future studies, we can try to use the Green–Ampt infiltration model or other runoff models for simulation and compare the simulation accuracy.

5. Conclusions

In this study, a 1D–2D coupled model constructed with InfoWorks ICM was used to simulate the inundation scene in a typical inundation area in the north of Lin’an’s main urban area. The model was verified by two historical rainfalls. Different designed rainfalls were input into the coupled model to analyze the spatial-temporal variation in urban inundation in the study area. The main conclusions are listed below:
(1)
The simulation results of the 1D–2D coupled model were acceptable. The validation results of two historical rainfalls showed that the NSE values were all above 0.82, the R2 values were all above 0.87, and the relative errors were all ±20%. Meanwhile, the simulating depth was consistent with the observed depth. These results indicated that the coupled model could accurately simulate the inundation situation in the study area.
(2)
Scenario analysis demonstrated that with the increase in the return period, rainfall peak position coefficient, and rainfall duration, the maximum inundation depth and inundation extent increased. The increased inundation area was mainly concentrated in the south of the study area. Combined with the change in the inundation of the local area, inundation mainly occurred in areas with low terrain (such as AFR), lack of drainage network (such as MXR), and narrow pipe diameter (such as WSS). Rainfall patterns also affected inundation, as in the case of AFR.
The inundation analysis under the designed rainfall scenarios will help us better understand the main causes and changing characteristics of inundation in Lin’an. Furthermore, it will help policy makers in Lin’an City to formulate appropriate urban inundation prevention measures, such as increasing small pipe diameters or filling low-lying areas, to reduce the risk of urban inundation. The results of this study can also be used as a reference for solving urban inundation problems in other cities with a similar geological and climatic condition.

Author Contributions

Conceptualization, Y.C. and T.H.; methodology, Y.C., Y.L. and T.H.; formal analysis, Y.C. and L.W.; investigation, Y.C., L.W. and J.F.; resources, T.H. and B.W.; data curation, Y.C., T.H. and B.W.; writing—original draft preparation, Y.C.; writing—review and editing, Y.C., T.H., H.H. and Y.L.; supervision, T.H. and H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Zhejiang Provincial Natural Science Foundation of China, grant number: LY19D010004, and the Science and Technology Program of Hangzhou, grant number: 20191203B19.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to express gratitude to Pin Wang (Hangzhou Normal University) and Yuzhou Zhang (Zhejiang University of Water Resources and Electric Power) for their valuable support throughout this study. Furthermore, authors thank the main representative officers of Lin’an institutions involved in the data collection process.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC. Climate Change 2013–The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014; ISBN 9781107057999. [Google Scholar]
  2. Karamouz, M.; Hosseinpour, A.; Nazif, S. Improvement of Urban Drainage System Performance under Climate Change Impact: Case Study. J. Hydrol. Eng. 2011, 16, 395–412. [Google Scholar] [CrossRef]
  3. Kong, F.; Sun, S. Better Understanding Insurance Mechanism in Dealing with Climate Change Risk, with Special Reference to China. Int. J. Environ. Res. Public Health 2021, 18, 2996. [Google Scholar] [CrossRef] [PubMed]
  4. Yin, J.; Yu, D.; Yin, Z.; Wang, J.; Xu, S. Multiple Scenario Analyses of Huangpu River Flooding Using a 1D/2D Coupled Flood Inundation Model. Nat. Hazards 2013, 66, 577–589. [Google Scholar] [CrossRef]
  5. Chen, W.; Huang, G.; Zhang, H. Urban Stormwater Inundation Simulation Based on SWMM and Diffusive Overland-Flow Model. Water Sci. Technol. 2017, 76, 3392–3403. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, M.; Chen, W.; Cai, K.; Gao, X.; Zhang, X.; Liu, J.; Wang, Z.; Li, D. Analysis of the Spatial Distribution Characteristics of Urban Resilience and Its Influencing Factors: A Case Study of 56 Cities in China. Int. J. Environ. Res. Public Health 2019, 16, 4442. [Google Scholar] [CrossRef]
  7. Nimmo, J.R. The Processes of Preferential Flow in the Unsaturated Zone. Soil Sci. Soc. Am. J. 2021, 85, 1–27. [Google Scholar] [CrossRef]
  8. De Carlo, L.; Perkins, K.; Caputo, M.C. Evidence of Preferential Flow Activation in the Vadose Zone via Geophysical Monitoring. Sensors 2021, 21, 1358. [Google Scholar] [CrossRef]
  9. Scaini, A.; Hissler, C.; Fenicia, F.; Juilleret, J.; Iffly, J.F.; Pfister, L.; Beven, K. Hillslope Response to Sprinkling and Natural Rainfall Using Velocity and Celerity Estimates in a Slate-Bedrock Catchment. J. Hydrol. 2018, 558, 366–379. [Google Scholar] [CrossRef]
  10. Hammond, M.J.; Chen, A.S.; Djordjević, S.; Butler, D.; Mark, O. Urban Flood Impact Assessment: A State-of-the-Art Review. Urban Water J. 2015, 12, 14–29. [Google Scholar] [CrossRef]
  11. Zheng, Z.; Gao, J.; Ma, Z.; Wang, Z.; Yang, X.; Luo, X.; Jacquet, T.; Fu, G. Urban Flooding in China: Main Causes and Policy Recommendations. Hydrol. Process. 2016, 30, 1149–1152. [Google Scholar] [CrossRef]
  12. YNET. Typhoon “Fireworks” Caused Direct Economic Losses of 3.35 Billion Yuan to 2.711 Million People in 4 Provinces and Cities. Available online: https://t.ynet.cn/baijia/31182698.html (accessed on 28 February 2022).
  13. Yin, J.; Yu, D.; Wilby, R. Modelling the Impact of Land Subsidence on Urban Pluvial Flooding: A Case Study of Downtown Shanghai, China. Sci. Total Environ. 2016, 544, 744–753. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, Q.; Wu, Z.; Zhang, H.; Dalla Fontana, G.; Tarolli, P. Identifying Dominant Factors of Waterlogging Events in Metropolitan Coastal Cities: The Case Study of Guangzhou, China. J. Environ. Manag. 2020, 271, 110951. [Google Scholar] [CrossRef] [PubMed]
  15. Han, H.; Hou, J.; Xu, Z.; Jing, H.; Gong, J.; Zuo, D.; Li, B.; Yang, S.; Kang, Y.; Wang, R. A GPU-Accelerated Hydrodynamic Model for Urban Rainstorm Inundation Simulation: A Case Study in China. KSCE J. Civ. Eng. 2022, 26, 1494–1504. [Google Scholar] [CrossRef]
  16. Cen, G.P.; Zhan, D.J.; Hong, J.N. Urban Storm Pipeline Calculation Model. China Water Wastewater 1993, 9, 37–40. [Google Scholar] [CrossRef]
  17. Zhou, Y.W.; Zhao, H.B. A Study On Urban Storm Runoff Model. China Water Wastewater 1997, 13, 4–6+2. [Google Scholar] [CrossRef]
  18. Gironás, J.; Roesner, L.A.; Rossman, L.A.; Davis, J. A New Applications Manual for the Storm Water Management Model (SWMM). Environ. Model. Softw. 2010, 25, 813–814. [Google Scholar] [CrossRef]
  19. Gong, Y.; Li, X.; Zhai, D.; Yin, D.; Song, R.; Li, J.; Fang, X.; Yuan, D. Influence of Rainfall, Model Parameters and Routing Methods on Stormwater Modelling. Water Resour. Manag. 2018, 32, 735–750. [Google Scholar] [CrossRef]
  20. Innovyze. InforWorks ICM Product Information-Overview. Available online: https://www.innovyze.com/en-us/products/infoworks-icm (accessed on 26 January 2022).
  21. Liu, Z.; Cai, Y.; Wang, S.; Lan, F.; Wu, X. Small and Medium-Scale River Flood Controls in Highly Urbanized Areas: A Whole Region Perspective. Water 2020, 12, 182. [Google Scholar] [CrossRef]
  22. Samantaray, D.; Chatterjee, C.; Singh, R.; Gupta, P.K.; Panigrahy, S. Flood Risk Modeling for Optimal Rice Planning for Delta Region of Mahanadi River Basin in India. Nat. Hazards 2015, 76, 347–372. [Google Scholar] [CrossRef]
  23. Xia, J.; Zhang, Y.; Liang, C.M.; Liu, J. Review on Urban Storm Water Models. Eng. J. Wuhan Univ. 2018, 51, 95–105. [Google Scholar] [CrossRef]
  24. Wang, H.C.; Chen, J.G.; Zhang, S.H.; Lai, H.L.; Su, D.B. Application Status and Comparative Analysis of Urban Storm Water Models. Water Resour. Hydropower Eng. 2011, 42, 10–13. [Google Scholar] [CrossRef]
  25. Cheng, T.; Xu, Z.; Hong, S.; Song, S. Flood Risk Zoning by Using 2D Hydrodynamic Modeling: A Case Study in Jinan City. Math. Probl. Eng. 2017, 2017, 5659197. [Google Scholar] [CrossRef]
  26. Yu, H.; Huang, G.; Wu, C. Application of the Stormwater Management Model to a Piedmont City: A Case Study of Jinan City, China. Water Sci. Technol. 2014, 70, 858–864. [Google Scholar] [CrossRef] [PubMed]
  27. Patro, S.; Chatterjee, C.; Mohanty, S.; Singh, R.; Raghuwanshi, N.S. Flood Inundation Modeling Using MIKE FLOOD and Remote Sensing Data. J. Indian Soc. Remote Sens. 2009, 37, 107–118. [Google Scholar] [CrossRef]
  28. Sidek, L.M.; Jaafar, A.S.; Majid, W.H.A.W.A.; Basri, H.; Marufuzzaman, M.; Fared, M.M.; Moon, W.C. High-resolution Hydrological-hydraulic Modeling of Urban Floods Using Infoworks Icm. Sustainability 2021, 13, 10259. [Google Scholar] [CrossRef]
  29. Yin, D.; Evans, B.; Wang, Q.; Chen, Z.; Jia, H.; Chen, A.S.; Fu, G.; Ahmad, S.; Leng, L. Integrated 1D and 2D Model for Better Assessing Runoff Quantity Control of Low Impact Development Facilities on Community Scale. Sci. Total Environ. 2020, 720, 137630. [Google Scholar] [CrossRef]
  30. Salvadore, E.; Bronders, J.; Batelaan, O. Hydrological Modelling of Urbanized Catchments: A Review and Future Directions. J. Hydrol. 2015, 529, 62–81. [Google Scholar] [CrossRef]
  31. Papaioannou, G.; Loukas, A.; Vasiliades, L.; Aronica, G.T. Flood Inundation Mapping Sensitivity to Riverine Spatial Resolution and Modelling Approach. Nat. Hazards 2016, 83, 117–132. [Google Scholar] [CrossRef]
  32. Peng, H.; Liu, Y.; Wang, H.; Gao, X.; Chen, Y.; Ma, L. Urban Stormwater Forecasting Model and Drainage Optimization Based on Water Environmental Capacity. Environ. Earth Sci. 2016, 75, 1094. [Google Scholar] [CrossRef]
  33. Pyatkova, K.; Chen, A.S.; Butler, D.; Vojinović, Z.; Djordjević, S. Assessing the Knock-on Effects of Flooding on Road Transportation. J. Environ. Manag. 2019, 244, 48–60. [Google Scholar] [CrossRef]
  34. Zhou, Y.; Shen, D.; Huang, N.; Guo, Y.; Zhang, T.; Zhang, Y. Urban Flood Risk Assessment Using Storm Characteristic Parameters Sensitive to Catchment-Specific Drainage System. Sci. Total Environ. 2019, 659, 1362–1369. [Google Scholar] [CrossRef] [PubMed]
  35. Lin’an Release. Flood Bulletin. Available online: https://weibo.com/5181771928/Co1e0AD7W?refer_flag=1001030103_ (accessed on 28 February 2022).
  36. Hangzhou Daily. Lin’an Reports the Situation of Disaster Relief. Available online: http://www.hangzhou.gov.cn/art/2019/8/12/art_812262_36845347.html (accessed on 12 August 2019).
  37. Peng, H.Q.; Liu, Y.; Wang, H.W.; Ma, L.M. Assessment of the Service Performance of Drainage System and Transformation of Pipeline Network Based on Urban Combined Sewer System Model. Environ. Sci. Pollut. Res. 2015, 22, 15712–15721. [Google Scholar] [CrossRef] [PubMed]
  38. Zhejiang Provincial Department of Housing and Urban-Rural Development. Announcement on the Publication of the Engineering Construction Standard “Rainstorm Intensity Calculation Standard” in Zhejiang Province. Available online: http://jst.zj.gov.cn/art/2020/3/6/art_1228990170_208.html (accessed on 6 October 2020).
  39. Cen, G.P.; Shen, J.; Fan, R.S. Research on Rainfall Pattern of Urban Design Storm. Adv. Water Sci. 1998, 9, 42–47. [Google Scholar]
  40. Song, J.; Wang, J.; Xi, G.; Lin, H. Evaluation of Stormwater Runoff Quantity Integral Management via Sponge City Construction: A Pilot Case Study of Jinan. Urban Water J. 2021, 18, 151–162. [Google Scholar] [CrossRef]
  41. Zhang, H.; Yang, Z.; Cai, Y.; Qiu, J.; Huang, B. Impacts of Climate Change on Urban Drainage Systems by Future Short-Duration Design Rainstorms. Water 2021, 13, 2718. [Google Scholar] [CrossRef]
  42. EPA. Storm Water Management Model (SWMM) Version 5.1 User’s Manual. Available online: https://www.epa.gov/water-research/storm-water-management-model-swmm-version-51-users-manual (accessed on 10 March 2022).
  43. Chen, W.; Huang, G.; Zhang, H.; Wang, W. Urban Inundation Response to Rainstorm Patterns with a Coupled Hydrodynamic Model: A Case Study in Haidian Island, China. J. Hydrol. 2018, 564, 1022–1035. [Google Scholar] [CrossRef]
  44. Liew, Y.S.; Desa, S.M.; Md. Noh, M.N.; Tan, M.L.; Zakaria, N.A.; Chang, C.K. Assessing the Effectiveness of Mitigation Strategies for Flood Risk Reduction in the Segamat River Basin, Malaysia. Sustainability 2021, 13, 3286. [Google Scholar] [CrossRef]
  45. Wu, X.; Wang, Z.; Guo, S.; Liao, W.; Zeng, Z.; Chen, X. Scenario-Based Projections of Future Urban Inundation within a Coupled Hydrodynamic Model Framework: A Case Study in Dongguan City, China. J. Hydrol. 2017, 547, 428–442. [Google Scholar] [CrossRef]
  46. Mei, C.; Liu, J.H.; Wang, H.; Li, Z.J.; Yang, Z.Y.; Shao, W.W.; Ding, X.Y.; Weng, B.S.; Yu, Y.D.; Yan, D.Y. Urban Flood Inundation and Damage Assessment Based on Numerical Simulations of Design Rainstorms with Different Characteristics. Sci. China Technol. Sci. 2020, 63, 2292–2304. [Google Scholar] [CrossRef]
  47. Jiang, W.; Yu, J. Impact of Rainstorm Patterns on the Urban Flood Process Superimposed by Flash Floods and Urban Waterlogging Based on a Coupled Hydrologic–Hydraulic Model: A Case Study in a Coastal Mountainous River Basin within Southeastern China. Nat. Hazards 2022, 112, 301–326. [Google Scholar] [CrossRef]
  48. Zhang, M.; Xu, M.; Wang, Z.; Lai, C. Assessment of the Vulnerability of Road Networks to Urban Waterlogging Based on a Coupled Hydrodynamic Model. J. Hydrol. 2021, 603, 127105. [Google Scholar] [CrossRef]
  49. Li, Y.; Hu, T.; Zheng, G.; Shen, L.; Fan, J.; Zhang, D. An Improved Simplified Urban Storm Inundation Model Based on Urban Terrain and Catchment Modification. Water 2019, 11, 2335. [Google Scholar] [CrossRef]
  50. Cheng, T.; Xu, Z.; Yang, H.; Hong, S.; Leitao, J.P. Analysis of Effect of Rainfall Patterns on Urban Flood Process by Coupled Hydrological and Hydrodynamic Modeling. J. Hydrol. Eng. 2020, 25, 04019061. [Google Scholar] [CrossRef]
  51. Mallari, K.J.B.; Arguelles, A.C.C.; Kim, H.; Aksoy, H.; Kavvas, M.L.; Yoon, J. Comparative Analysis of Two Infiltration Models for Application in a Physically Based Overland Flow Model. Environ. Earth Sci. 2015, 74, 1579–1587. [Google Scholar] [CrossRef]
Figure 1. Location of the study area. AFR, Agricultural and Forestry Road; SLS, Shuanglin Street; LSR, Linshui Road; NLS, Nonglin Street; WSS, Wusu Street; XSS, Xishu Street; MXR, Maxi Road; XJSR, Xijingshan Road; BNR, Baini Road. (a) Zhejiang Province (b) Study area.
Figure 1. Location of the study area. AFR, Agricultural and Forestry Road; SLS, Shuanglin Street; LSR, Linshui Road; NLS, Nonglin Street; WSS, Wusu Street; XSS, Xishu Street; MXR, Maxi Road; XJSR, Xijingshan Road; BNR, Baini Road. (a) Zhejiang Province (b) Study area.
Ijerph 19 07210 g001
Figure 2. Hydrographs of the designed rainfall events (t = 120 min).
Figure 2. Hydrographs of the designed rainfall events (t = 120 min).
Ijerph 19 07210 g002
Figure 3. Technical flowchart.
Figure 3. Technical flowchart.
Ijerph 19 07210 g003
Figure 4. Comparison between simulated and observed flows.
Figure 4. Comparison between simulated and observed flows.
Ijerph 19 07210 g004
Figure 5. Simulation results under designed rainfall scenarios, where “1a-0.3” means P = 1a and r = 0.3 (t = 120 min).
Figure 5. Simulation results under designed rainfall scenarios, where “1a-0.3” means P = 1a and r = 0.3 (t = 120 min).
Ijerph 19 07210 g005
Figure 6. The inundation scene with different rainfall peak position coefficients (t = 120 min and P = 5a).
Figure 6. The inundation scene with different rainfall peak position coefficients (t = 120 min and P = 5a).
Ijerph 19 07210 g006
Figure 7. The inundation scene at different durations (P = 5a and r = 0.4).
Figure 7. The inundation scene at different durations (P = 5a and r = 0.4).
Ijerph 19 07210 g007
Figure 8. The inundation result from rainfall on 2 July 2020 ((A) is MXR, (B) is XSS, and (C) is SLS).
Figure 8. The inundation result from rainfall on 2 July 2020 ((A) is MXR, (B) is XSS, and (C) is SLS).
Ijerph 19 07210 g008
Figure 9. The inundation situation in MXR (the image came from aerial photography).
Figure 9. The inundation situation in MXR (the image came from aerial photography).
Ijerph 19 07210 g009
Figure 10. Inundation of the study area during different return periods (t = 120 min and r = 0.4): (A) XJSR, (B) NLS, (C) AFR, (D) WSS.
Figure 10. Inundation of the study area during different return periods (t = 120 min and r = 0.4): (A) XJSR, (B) NLS, (C) AFR, (D) WSS.
Ijerph 19 07210 g010
Table 1. List of geographic data in the study area.
Table 1. List of geographic data in the study area.
TypeFormatResolution (m)Main AttributesData Source
DEMGeoTIFF2ElevationSurveying and mapping department of Lin’an
UAV imageGeoTIFF0.5
Land-use dataShapefileLand-use typeMeteorological department of Lin’an
Rainfall dataExcelTime and rainfall
Drainage system dataShapefilePipe diameter, pipe materialsUrban management department of Lin’an
Table 2. Parameter attributes for five kinds of surfaces.
Table 2. Parameter attributes for five kinds of surfaces.
Surface TypeRoutine
Model
Routine
Parameter
Surface
Type
Runoff
Model
Runoff
Coefficient
Initial
Loss (m)
RoadSWMM0.02ImperviousFixed0.90.0015
BuildingSWMM0.02ImperviousFixed0.80.001
OthersSWMM0.025ImperviousFixed0.50.005
WaterSWMM0.03ImperviousFixed10
Green spaceSWMM0.2PerviousHorton0.005
Table 3. Model validation results.
Table 3. Model validation results.
EventsPeak Flow (m3/s)NSER2RMSERelative Error (%)
RecordSimulation
29 May 20200.5690.6030.820.940.076.0
2 July 20201.3511.1500.850.870.15–14.9
Table 4. Statistics of the maximum inundation depth of the 2 July 2020 rainfall event.
Table 4. Statistics of the maximum inundation depth of the 2 July 2020 rainfall event.
NumberPositionRecorded Depth (cm)Simulated Depth (cm)Errors (cm)
1MXR2016.93.1
2XSS1011.7–1.7
3SLS1011.3–1.3
Table 5. The inundation depths and extents in the study area under different designed rainfall scenarios (t = 120 min).
Table 5. The inundation depths and extents in the study area under different designed rainfall scenarios (t = 120 min).
Return Period (a)rNumber of Overflow NodesMaximum Overflow Volume of Nodes (m3/s)Maximum Inundation Depth (m)Inundation Extent (m2)Contribution Area of Different Inundation Depths (m2)
≤0.2 m0.2–0.4 m0.4–0.6 m0.6–0.8 m0.8–1 m>1 m
10.3120.50.81380,78175,77235151161309240
0.4130.50.84784,45378,71042281115376240
0.5130.50.8687,48381,45745111115376240
20.3160.60.946100,98193,34758091190564710
0.4180.60.978103,76995,64562991120620850
0.5200.70.99106,93298,46866031156620850
50.3280.71.096137,329126,2598708114771015748
0.4310.81.126148,234136,6528957171071013471
0.5320.81.139158,139146,1469368157984113471
100.3360.81.203188,477173,37711,171270292022285
0.4370.91.222195,118178,88411,8433084100022285
0.5440.91.235201,744185,23511,97432181010155152
Table 6. The simulation results of different rainfall durations (r = 0.4).
Table 6. The simulation results of different rainfall durations (r = 0.4).
Rainfall Duration (min)Return Periods (a)Cumulative Rainfall (mm)Maximum Inundation Depth (m)Inundation
Extent (m2)
Proportion of Inundation Area (%)
60140.30.80378,1154.13%
245.20.86499,0165.24%
551.61.093134,2007.11%
1056.41.199181,2299.59%
120150.80.84784,4534.47%
256.90.978103,7695.49%
565.01.126148,2347.84%
1071.11.222195,11810.33%
180157.30.85787,2364.62%
264.20.979106,4555.63%
573.31.135155,8238.25%
1080.21.229201,02710.64%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Hou, H.; Li, Y.; Wang, L.; Fan, J.; Wang, B.; Hu, T. Urban Inundation under Different Rainstorm Scenarios in Lin’an City, China. Int. J. Environ. Res. Public Health 2022, 19, 7210. https://doi.org/10.3390/ijerph19127210

AMA Style

Chen Y, Hou H, Li Y, Wang L, Fan J, Wang B, Hu T. Urban Inundation under Different Rainstorm Scenarios in Lin’an City, China. International Journal of Environmental Research and Public Health. 2022; 19(12):7210. https://doi.org/10.3390/ijerph19127210

Chicago/Turabian Style

Chen, Yan, Hao Hou, Yao Li, Luoyang Wang, Jinjin Fan, Ben Wang, and Tangao Hu. 2022. "Urban Inundation under Different Rainstorm Scenarios in Lin’an City, China" International Journal of Environmental Research and Public Health 19, no. 12: 7210. https://doi.org/10.3390/ijerph19127210

APA Style

Chen, Y., Hou, H., Li, Y., Wang, L., Fan, J., Wang, B., & Hu, T. (2022). Urban Inundation under Different Rainstorm Scenarios in Lin’an City, China. International Journal of Environmental Research and Public Health, 19(12), 7210. https://doi.org/10.3390/ijerph19127210

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop