Next Article in Journal
Time-Series Analysis of Isotope Composition of Precipitation in Zagreb, Croatia
Previous Article in Journal
Projection of Future Water Resources Carrying Capacity in the Huang-Huai-Hai River Basin under the Impacts of Climate Change and Human Activities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantifying the UDS Hydraulic and Social Resilience to Flooding: An Index-Based Approach vs. a Parameter-Based MCDM Method

1
Department of Environmental Engineering, Faculty of Engineering, University of Tehran, Tehran 1417853111, Iran
2
Institute for Infrastructure Engineering, University of Innsbruck, 6020 Innsbruck, Austria
3
Department of Engineering, University of Messina, 98158 Messina, Italy
*
Authors to whom correspondence should be addressed.
Water 2022, 14(13), 2007; https://doi.org/10.3390/w14132007
Submission received: 28 April 2022 / Revised: 7 June 2022 / Accepted: 15 June 2022 / Published: 23 June 2022
(This article belongs to the Section Hydrology)

Abstract

:
Various approaches to quantifying resilience have sparked debate in past years. This study measures the social and hydraulic resilience of an Urban Drainage System (UDS) using two methodologies, i.e., a proposed index-based framework and a parametric-based Multi-Criteria Decision Making (MCDM) method. The former positions flexibility as a key property of the resilience concept, and after quantifying four characteristics of a flexible system, measures the flood resilience based on quantifying the two resilience components, i.e., flexibility and resistance; while. The latter calculates the resilience through a linear combination of the actual values and relative weights for different UDS sub-characteristics which best describe the resiliency of the system. The methodologies were applied to a basin, focusing on flooding as a disaster, to quantitatively evaluate the behavior of UDS under both single-event and long-term precipitation. The results of both methods were indicative of a relatively low level of flood resiliency for the UDS in the studied area. Among different scenarios examined for the studied region, using Best Management Practices (BMPs) was turned out to be an effective adaptive measure to enhance the UDS resilience to a higher level.

Graphical Abstract

1. Introduction

Building resilience into urban critical infrastructures is essential to have the urban daily functions continue as normal in the face of various disasters. Measuring resilience is a crucial step in gaining such an achievement. It is known that resilience is the ability of a system to revert to its previous state after being disturbed by a perturbation. Flexibility, on the other hand, is sometimes used in the same meaning as resilience [1], which means there is a lack of focus on this concept as a different notion from resilience. Although- in some of the past studies and particularly in the field of psychology- flexibility has been mentioned as a close concept to resilience and related to it, it has not been characterized as an element under resilience structure. In other words, the concept of flexibility has not been well placed within the resilience mechanism. Only a few studies have focused on the notion of flexibility in detail, while the quantification of this concept has been missing in most cases. Some studies believe that a flexible system has the main capabilities of resilience.
However, most studies agree that robustness (as an attribute of the resistance concept) could be at least a part of resilience. Restemeyer et al. [2] argue that resistance can be one important aspect of resilience. In line with this, Tempels and Hartmann [3] believe that prioritizing resistance or flexibility over each other does not seem to be true and state that there should be a balanced approach towards flexibility and robustness [4]. Asokan et al. [5] defined flexibility as the representative of two main resilience parameters, i.e., robustness and transformation, but not as a factor inside the resilience configuration. In the current study, we have a new definition of resilience based on system flexibility. In the proposed definition, flexibility is nearly the same as resilience, excluding one metric (the robustness), which is the only term that resilience has more than flexibility. In fact, flexibility is most of the resilience concept. Such an approach seems to agree with the idea represented by Husdal [6], who believes robustness and flexibility are essential tools for handling risks and uncertainties.
It could be said that despite the frequent use of the term flexibility, very little work has been accomplished to quantify this term, especially regarding the urban drainage system. Some studies focused on flexibility measurements concerning the other meanings of flexibility suitable for the other types of systems and provided with the quantification methods are not necessarily applicable to the UDS. For example, Stinner et al. [7] proposed a quantification method of flexibility specific to the Building Energy Systems with Heat Pumps. Fletcher [8] addressed flexibility in Water Supply Infrastructure Planning. However, the purpose was not to quantify the flexibility or resilience of the system. Sidqi et al. [9] defined flexibility in the household as the capacity to increase, decrease, or shift the consumption of domestic appliances over time and quantified the consumed power in 194 households, and computed the share of both heating needs, i.e., space heating and hot water heating.
Almost all previous studies agreed on the measuring method of UDS robustness in the face of extreme rainfall (which is accepted by this paper and is explained later). But the question arises: how is the system flexibility measured? This study mainly focuses on this part of quantifying the resilience concept. In the proposed approach, finding the answer to the seven following questions is necessary: (1) how much of the disturbance can be tolerated by the system? (2) how fast does the system degrade in facing the disturbance? (3) how big is the response by the system when hit by an external disturbance? (4) how long does it take for the system to be able to survive from the minimum failure state? (5) how fast is the process of regaining stability and restoration to an acceptable performance level? (6) how much improvement is obtained after the recovery process? And (7) how well has the system been adapted to the new condition? The index-based framework discussed in the following sections is proposed to answer the above questions.
The results of the proposed index-based method will be compared with results from a parametric-based methodology, which has an emphasis on the multi-dimensional nature of resilience and is developed based on a Multi-Criteria Decision Making (MCDM) method.
This study attempts to fill in the gaps in previous works by concentrating on an urban drainage system resilience and flexibility from various viewpoints, i.e., hydraulic and social aspects, and instant and long-term evaluation of the system. Tahmasebi Birgani and Yazdandoost [10] proposed a framework to quantify resilience in the case of a drainage system facing a single-event rainfall based on three indicators of “reaction magnitude”, “graduality”, and “recovery capacity”. They suggested that future studies evaluate the long-term performance of UDS while going under long-term rainfall events. This has been focused on as a part of the current study. Most other studies, such as the works done by Mugume and Butler [11] and Bakhshipour et al. [12], quantified the UDS functional/structural resilience; however, they didn’t consider the quantification of resilience indicators and, instead, defined direct formulas for calculating the system resiliency. This is while the current study introduces indicator-based methodologies for UDS resilience quantification in the face of extreme rainfall events.

2. Literature Review

To provide a framework for quantifying the resilience of a disturbed system based on the notions of robustness and flexibility, after reviewing some of the most well-known definitions provided in recent years, it should be clarified what is meant exactly by these three concepts to be able to clarify how they are related to each other.

2.1. Definition of Resilience Based on “Robustness” and “Flexibility”

Resilience is defined as the system’s capacity to cope, adapt, and improve in response to both chronic stresses and acute shocks applied to it [13,14]. Before explaining how such a definition could be expressed based on two elements of “robustness” and “flexibility”, we had better have a comparison between the two components.
Robustness points to the ability of a system to “stay on course” and absorb unforeseen and uncertain future disturbances [6]. A Robust system accommodates the perturbations and resists the disturbance so that it continues its pre-defined path without functional degradation. This is why flexibility means the capacity to accommodate unforeseen disturbances “by changing tracks and being open to temporarily deviate from the initial course to reach the desired goal” [6]. A flexible system can gracefully degrade while exposed to a disturbance and restructure itself rapidly with an adaptive approach and by learning from past experiences. It seems that concerning the flexibility concept, what is important is the final goal without being hurt, no matter how the path has to be paced during a disaster, whereas robustness focuses on the capacity of the system to endure such changes [15]. Totally, as Husdal [6] states, “robustness reduces uncertainty and flexibility adapts to it”.
According to resilience’s general definition and descriptions of system robustness and flexibility presented above, the resilience of a system is highly dependent on these two concepts. Thus, measuring the system resilience requires quantifying the subset components, i.e., robustness and flexibility. Walker et al. [16] define flexibility as “dynamic robustness”, which enables a plan to change in response to changing conditions over time.
A literature review indicates that flexibility has been viewed from two different angles based on where on the system functionality curve (Figure 1) it will be considered: in the response phase or within the recovery phase and after that. In the first meaning, flexibility refers to the degree of proportionality in response to distribution and means not having an abrupt failure due to a sudden shock. In this meaning, a system is flexible if the applied perturbation is mostly absorbed, and therefore, consequences or damages emerge little by little. Based on such a definition, with a flexible system, there is enough time to mobilize the equipment and resources and cope with the impacts of the shock while the damages are not too large. This meaning of flexibility has been expressed in the current study as “braking capacity”, which has a close relationship with the notion of “graduality” introduced by De Bruijn [17]. Graduality refers to the increase in damage with increased disturbance magnitude [18]. This is why braking capacity represents the trend of system degradation with time. (For more information on the difference between the two concepts of graduality and braking capacity, see Supplementary Materials, Figure S1).
In the second meaning, flexibility is related to uncertainties in future conditions and, thus, uncertainties in estimations for future planning. In this meaning, flexibility refers to the willingness and ability to adopt alternative strategies in response to changing circumstances or sudden crises [13]. In other words, despite the primary estimations for future conditions and planned implementations, a flexible system could still be stable and welcome such unforeseen changes. Heggarty et al. [19] considered this meaning of flexibility for power systems and defined it as “the ability to cope with variability and uncertainty in demand”. This translation of flexibility has been referred to as “adaptive capacity” in the current study.
However, some other descriptions of flexibility can be found in the literature covering both above-mentioned definitions, according to which flexibility is described as system change, evolution, and adaptation in response to changing circumstances. This definition matches the most to what the current study means by flexibility. The current definition encompasses a wide range of system behavior in the face of an external perturbation and involves both the meanings pointed out above. This study believes that flexibility is composed of four metrics of response, survivability, recovery, and adaptability, and defines flexibility as the ability of a disturbed system to have (1) a soft response (i.e., with a small response-magnitude and a graceful degradation), (2) a short stay in the failure mode and ability to immediate rise-up, (3) a rapid restoration to an acceptable level of performance or ideally, an improved one, and (4) adapt to the new conditions after changing the surrounding environment or the system condition.

2.2. System Behavior in Response to a Disturbance

Before quantifying the resilience components, it is necessary to explain the system performance and behavior as well as the critical thresholds in response to an external perturbation. The behavior of a disturbed system vs. time has been presented in Figure 1. Three types of critical thresholds can be recognized in such a performance curve. The first one is the “reaction threshold”, which is located at the end of the resistance domain. After that, the system is no longer able to endure the applied perturbation and thus, starts responding to it and degrades. Concerning urban flooding, the reaction threshold occurs when we have the largest load on the system (i.e., the largest amount of extreme rainfall), which is not supposed to create flooding [20,21,22]. At the end of the response process, when the system performance reaches its minimum level possible, the “tolerance threshold” occurs, which means a threshold where the system cannot bear the disturbance anymore. This is why before this threshold, a part of the disturbance could be absorbed by the system alongside the performance degradation. In fact, the “tolerance threshold” is located before the maximum failure of the system. Apart from the two mentioned thresholds, a third could be identified in the system behavior curve, which is called the “non-recovery” threshold and is not determined in Figure 1 since it occurs only in certain conditions. The “non-recovery” threshold happens when the system is no longer able to return to its previous status and enters a new attraction basin and new stability domain. A regime change happens in system performance and status, which may occur when the failure time duration and system stagnation is too long. If the magnitude of the system response to a disturbance is so large that it causes the system to spend a long time at its minimum level of performance, or if the disturbance magnitude is too large and goes beyond certain limits, the system is no longer able to return to its previous stability status, which is where the “non-recovery” threshold makes sense [23,24].
As shown in Figure 1, the system behavior curve experiences five different stages when exposed to an external disturbance (i.e., resistance, response, rise-up, recovery, and adaptability phases). Quantifying these phases means measuring the components of a disaster cycle, from resistance against a disaster to recovery, mitigation, and preparedness for future similar hazards. Figure 2 demonstrates the components of system resilience in the proposed framework.

3. Introduction to the Case Study

A semi-urbanized basin in the northern areas of Tehran, Iran, was considered the case study in which the proposed methodology (explained in the following section) would be applied. Three major rivers—Darakeh, Farahzad, and Hesarak—are located in this basin whose waters discharge to a constructed diversion channel west of Tehran, which ultimately joins the Kan River. In the current study, we call the basin the “West Flood Diversion (WFD)” catchment. The basin covers an area of more than 140 Km2 and involves both urban and suburban mountainous areas upstream of the urbanized region. Figure 3 depicts the situation of the catchment in Tehran and Iran. The blue and gray areas in the basin represent the urban and suburban areas, respectively. Black lines represent the sub-basins’ boundaries. In addition, Table 1 represents the physiographic characteristics of the sub-basins.

4. Materials and Methods

Figure 4 represents the flowchart of the works accomplished in the current study. This flowchart gives a total view of the process of measuring resilience and flexibility, comparing the results with an MCDM method and choosing the best solutions for improving the resilience index of the system, which are discussed in detail in the following sections.

4.1. Index-Based Measurement of Resilience

4.1.1. System Definition and Characterizing the Disturbance and the System Response

To use the framework, the system boundaries should first be defined, and the type of disturbance applied to the system and its main characteristics should be determined. The system under investigation in the current study is Urban Drainage Infrastructure (UDI), embedded in the surrounding urban complex environment and is in mutual dependency with it through social feedback. Since this study will examine the resilience of UDS to flooding, any cases of extreme rainfall, pluvial flooding, rapid urbanization, or climatic changes (as the source of extreme events) could be considered the disturbance. Herein, the precipitation was chosen as the applied perturbation for the drainage system, and therefore, the surcharges from the UDS would be the response to such a disturbance. The response magnitude would be the surcharged depth, while the response trend would be determined through the increase in inundation depth with time.
To determine the characteristics of the extreme rainfalls (as the perturbation) and the system surcharges (as a consequence), the Intensity-Duration-Frequency (IDF) curves alongside the Stormwater Management Model (SWMM) were utilized (Sensitivity analysis and calibration of the rainfall-runoff model are presented in Supplementary Material, Tables S1–S3). The total rainfall depth was extracted from IDF curves, and parameters such as the volume/depth of the surcharged runoff and flood duration were obtained from the simulation results in EPA SWMM (as well as PCSWMM). SWMM is a dynamic hydrology-hydraulic water quality simulation model often used for a single event or continuous simulation of runoff quantity and quality in urban areas [26,27].
It should be noted that what is meant by flood in the current study is any surcharging from the drainage network, which inundates the surrounding area. And what is meant by extreme rainfall is the precipitation that results in a surcharge from UDS and thus causes inundation in the surrounding area.

4.1.2. Developing IDF Curves for the Region

The Equations (1) and (2) have been proposed by Ghahreman [28] for the condition of Iran—based on the relation offered by Bell [29]—which calculate the rainfall with a duration of 15 min to 20 hrs and a return period of 2 to 100 years [30].
P T t = 0.4524 + 0.2471 ln T 0.6000 0.3710 + 0.6184 t 0.4484 P 10 60
In Equation (1), t represents the time (min), T is the return period (year), P 10 60 is the amount of 10-year, 1-h rainfall, which is obtained through Equation (2) [31,32]:
P 10 60 = e 0.291 × P 24 h 0.694
where P 24 h is the maximum average of daily rainfall (mm).
To develop IDF curves for the studied region, Equations (1) and (2) were modified and calibrated based on the IDF equation presented for Tehran by Mahab-e-Ghods Consulting Engineers [33] (Equation (3)).
P = ( C Alt . RP D 0.645 ) · D / 60
In which P is the rainfall value and D is the rainfall duration in minutes ( D 360   min ). C Alt . RP is a coefficient that is calculated from Equation (4) and is dependent on the two factors of the return period (T, year) and average ground level (H, MASL).
C Alt . RP = 108 × 0.0008 × H + 0.2 × 0.304 × Ln T + 0.78928
After calibration and verification, the IDF curve for the WFD catchment could be obtained as Equation (5).
P T t = 0.3 + 0.114 ln T 0.57 0.3 + 0.6 t 0.4692 P 10 60
where:
P 10 60 = e 0.2 × P 24 h 0.6

4.1.3. Calculating the System Robustness

Robustness could be evaluated as the resistance phase in the system’s behavior curve, during which the disturbance is endured by the system, and no reaction is shown. This characteristic represents the tolerance capacity of the system under perturbation. On the system behavior curve, resistance refers to the duration from when the disturbance is applied to the system until the system can fully resist it (where no performance degradation occurs anymore), which is representative of “full absorption” of the disturbance by the system. For instance, in a normal system undergoing precipitation less than the designed rainfall (for which the system is designed), no response or reaction is made by the system. In such a situation, it could be said that the system behavior is in the resistance phase. In most previous studies [10,23,34,35], the resistance capacity (RC) has been quantified as a quota of the magnitude of the disturbance, which could be absorbed by the system without performance degradation and with no essential change in the system structure.
According to what was mentioned, two terms of Hydraulic Resistance Capacity (HRC) and Social Resistance Capacity (SRC) can be extracted for the current study. Hydraulic Resistance Capacity for Single-event rainfalls (HRCS) can be measured as the largest rainfall depth that could be tolerated by the system without any failures (Equation (7)). The amount of such rainfall with a certain return period could be extracted from the IDF curves of the studied region. The focus of the current study is on a 6-h design rainfall, which is common in the study site [36], and therefore, the largest amount of such a rainfall depth could be obtained through IDF curves for different return periods, and the surcharge from the system was calculated using the rainfall-runoff modeling in SWMM. However, it is better to normalize the system resistance to be comparable with the other cases and, simultaneously, apply it to the current study. To do this, the maximum rainfall depth that would not result in a surcharge from the UDS ( M a x ( R D i ) ) , is divided by Probable Maximum Precipitation (PMP) for the region.
M a x ( R D i ) < R D S H R C S = M a x ( R D i ) P M P
In Equation (7), R D S is the smallest depth of rainfall which results in a surcharge from the system, M a x ( R D i ) is the largest rainfall which does not cause a surcharge from UDS, and i is counter of the rainfall-induced runoff depths during the simulation process.
For long-term events, the basis for the calculation of the HRC would be the same as single event rainfall in general. However, since the return periods may not be known in a long-term continuous precipitation time series, HRC against a long-term rainfall disturbance (HRCL) can be defined as the quotient of the largest rainfall depth in the time series, which causes no surcharge from the UDS to the amount of PMP for the basin.
In the social context, quantifying the UDS robustness may be somehow different. When we are focusing on social problems regarding floods, not every inundation depth is going to be considered a hazard; therefore, we need to, first, have a threshold for that. Such a threshold may depend on the community under investigation. Herein, a 60 mm depth was considered a critical depth in the social context, and it is assumed that the flowing water with a depth of less than 60 mm could not make serious problems for the residents and commuters in the region [37]. On this basis, the largest rainfall value, which does not result in surface runoff of more than 60 mm in depth, is considered the Social Resistance Capacity against single-event rainfall (SRCS). For long-term precipitation, Social Resistance Capacity (SRCL) could be obtained as the quotient of the largest rainfall depth in the time series for which the surcharges from the drainage system do not exceed 60 mm to the PMP minus the given inundation threshold causing societal problems (herein 60 mm).
M a x ( R D i 60 ) < R D = 60   mm H R C S = M a x ( R D i 60 )
H R C L = M a x ( R D i 60 ) P M P 60 mm

4.1.4. Quantifying the System Flexibility

System Response

What is meant by the response is the system’s reaction to the disturbance. The response of the system can be argued from two viewpoints: response magnitude (i.e., degradation capacity, D C ) and response trend (i.e., braking capacity, B r C ). Before addressing the two components of the response phase, the quantification of absorptive capacity must be explained.

Absorptive Capacity

Absorptive capacity refers to that part of the behavior of a disturbed system that is resisting transformation under a disturbance. According to this definition, absorptive capacity is determined as the sum of the resistance domain and the quota of the response phase during which a part of the disturbance is absorbed by the system. Looking at a system behavior curve (Figure 1), the system may absorb the disturbance for a certain time duration and remain unchanged (i.e., resistance domain), but after reaching the reaction threshold, it is no longer able to endure the disturbance and, therefore, its track change (i.e., performance degradation) commences. In this phase (response domain), only a part of the disturbance is absorbed until the system reaches its minimum performance regarding the applied perturbation. This partial absorption in the response phase plus the total absorption happening during the resistance domain forms the “absorptive capacity” of a system. To quantify the absorptive capacity, the following equation is suggested:
A b C = R C + B r C + 1 / D C
where R C is the resistance capacity.

Degradation Capacity

The degradation capacity of a system is indicative of the fact as to what degree the system is impacted by the disturbance. In some previous works [1,17], the Expected Annual Damage (EAD) has been utilized to measure the response magnitude of a disturbed system, which is a commonly used index representing the expenses related to damages resulting from the perturbation applied to a system [38]. However, in the current study, the system response is described as the inundation depth considering the extreme rainfall as the disturbance on the UDS. (See Figure S2 (Supplementary Material) for more available approaches towards a system response to a disturbance). Therefore, the simulation results by SWMM can be utilized to determine this indicator for the system. In the current study, the Hydraulic Degradation Capacity (HDC) could be measured as the average surcharged runoff depth from the UDS. Social Degradation Capacity (SDC) is approximately the same other than subtracting the 60 mm surface runoff depth, which does not create significant social problems. This can be described as the following equation:
S D C = H D C 60 ( If S D C < 0 S D C = 0 )
where SDC and HDC are the values of social and hydraulic degradation capacities in mm.
For long-term time series, the average of surcharged runoff depths could represent the HDC, and accordingly, SDC can be calculated using Equation (9). In addition, SDC could be measured by considering the number of casualties and deaths during floods, the average delay time in reaching the destination, and increased traffic jams due to inundation problems. However, such a measurement requires having enough reliable data across the region. It is not the case in the current study since most of this kind of information is qualitative and not reliably documented.

Braking Capacity

Braking capacity (i.e., response trend) is described as the extent covering the onset point of performance reduction to the minimum performance level of the system (Figure 1). The longer the time duration is within this extent, or the gentler the performance reduction, the higher the braking capacity of the system.
This indicator is quite dependent on the “absorptive capacity” of the system, as the more the system absorption is, the less the probability of abrupt failure is, and the system degrades more gracefully—due to the absorption of a part of the disturbance by the system. According to the behavior of a disturbed system (Figure 3), the UDS braking capacity (BrC) is the inverse slope of the functionality degradation line. It can be calculated as described in Equation (10).
B r C = 1 i = 1 N D C i ( T P i ) N = N i = 1 N D C i ( T P i )
In the above equation, DCi is the degradation capacity of the ith flooded node (practically the surcharged flow-depth for that node), N is the total number of flooded nodes, and i is the counter. TP is time to peak, calculated as ( T M T 0 ) where T0 and TM are representative of the flooding start-time and peak-flood time, respectively. The higher the value of BrC is, the higher the system braking capacity (and consequently, resilience) is.
For calculating social braking capacity, TP would be the difference between T M and the time the surcharged runoff depth exceeds 60 mm (T60). Equation (10), which calculates the inverse slope of the response curve in Figure 1, can be utilized for both hydraulic and social braking capacities for single-event rainfall. To have the index for the long-term time series, Equation (11) can be used instead.
B r C L = N i = 1 N j = 1 M D C i j ( T P i j )
where M is the number of single rainfalls in the precipitation time series. Supplementary Material (Figure S2) explains the difference in a system’s behavior (i.e., response and recovery phases) under a similar disturbance.

4.1.5. Rise-Up Capacity

This study introduces an indicator describing a part of the system’s behavior, hereafter called “survivability”, which means the ability of the system to rise after a minimum failure state. In other words, the degree of the system agility in restoration after reaching the minimum possible performance under a specific disturbance. This system property has an inverse proportionality to the time duration during which the system remains in minimum failure mode. Considering the important parameter of flood depth as the multiplying weight for each node, the rise-up capacity (i.e., survivability) can be defined as follows:
S R V = i = 1 N 1 T S t = i = 1 N D i i = 1 N ( T S 1 i T M i ) D i
where T S t is the stagnation time, i.e., the time duration the system remains in the minimum failure mode, which is the time between maximum inundation ( T M ) and the time of flood subsidence in at least one node of the drainage system ( T S 1 ) . For the social aspect, T S 1 represents the time the flooding subsides below the 60 mm threshold at least in one node of the UDS. D i indicates the flood depth above the ith node of UDS, and therefore, SRV for the social aspect is meaningful only when the inundation depth reaches 60 mm. For long-term precipitation time series, this indicator should be calculated by averaging the obtained SRV for all extreme rainfall events (resulting in surcharges from the UDS). SRV values have been expressed with the unit 1/hr. A high value of SRV represents a system highly resilient to the given disturbance.

4.1.6. Restorative Capacity

By investigating the recovery phase in the system behavioral curve (Figure 1), we know the speed and the trend of the system in bouncing back to its previous (or acceptable) level of performance. Thus, recovery (i.e., restorative capacity) encompasses two characteristics of the system in response to disturbances: (1) system rapidity and (2) system improvability.

Rapidity

The Rapidity (RPD) indicator explains how fast the disturbed system returns to normal (or acceptable) performance after its functionality has reached the minimum possible level.
The recovery time (TR) is defined as the time duration the largest value of runoff depth starts to subside in one node ( T S 1 ) until it is 100% finished in all UDS nodes ( T S S ). To formulate UDS rapidity, the inversely weighted average of the recovery time is calculated. Inundation depth is the multiplying weight for each flooded node as the flood depth well represents the severity of the disaster. Accordingly, rapidity is formulated as the inverse of the weighted average of flood duration above the drainage nodes (Equation (15)).
R P D = i = 1 N D i i = 1 N T R i D i = i = 1 N D i i = 1 N ( T S S i T S 1 i ) D i
Equation (15) can obtain both hydraulic and social rapidity. To calculate social rapidity, T S S represents the time the flood depth reaches less than 60 mm in all UDS nodes ( T S S ), and like the SRV indicator, it is meaningless if the surcharged runoff from the system is below 60 mm. For long-term precipitation time-series, the rapidity would be the average of calculated RPD for any of the rainfall events.

Improvability

Improvability refers to the increase in the performance level (after resuming the normal functioning) compared to the pre-disaster level of functionality. Improvability is important for moving towards sustainability, which means the same disturbance with the same characteristics would not be able to disrupt the system efficiency for a second time. It is in agreement with the idea of some other researchers [14], who believe that a recovered system is expected to reach a higher resilience than the pre-disaster condition, and a “new normal” state is required at the end of this process, so that the system may even enter a new stability state and attraction basin.
To measure the improvability of the system, we use a modified version of the relationships proposed in the study by Binesh et al. [39] and calculate this system property as the difference between the UDS performances before and after implementing the adaptive solutions. The whole system’s improvability is obtained by averaging the calculated index for every single conduit, according to Equation (16) (provided that the modeling time steps are equal for all conduits’ simulations before and after implementing the adaptive strategies).
I M P = j = 1 N s V s b j V s a j V s b j N c
In Equation (16), Vsbj and Vsaj are representative of the surcharge volume in jth surcharged conduit in the UDS before the disaster and after the recovery process of the disturbed system, respectively; Ns is the number of surcharged conduits in the simulating process, and Nc is the whole number of UDS conduits. The IMP indicator varies between zero and one.
Equation (16) applies to single-event precipitation and long-term rainfall time series. Regarding the social aspect, a surcharged runoff depth equal to 60 mm (or less) plays the role of “no runoff” for the hydraulic improvability of UDS, which means only those surcharges from the system are accounted for that exceed 60 mm depth.
Herein, the improvability index was defined from a “volumetric” viewpoint. It could also be described from a “temporal” or “number of occurrence” point of view, which means how much improvability has been gained in terms of reduction in inundation duration time or frequency of surcharges from the system.

4.1.7. Adaptive Capacity

The factor of adaptability is of the most important components of resiliency as some studies consider the whole resiliency of a system mainly within this single property and define resiliency as “the ability of the system to adapt” [40]. It could be said that adaptability is considered a basic property and a key measure of resiliency [41,42]. An adaptive system changes its behavior in response to the change in the environment, reducing uncertainty. For example, using some adaptive solutions may accommodate different volumes of rainfall-induced runoff, and therefore, up to a certain threshold, no inundation (or little inundation) occurs. To calculate the adaptability of UDS (after considering the adaptive solutions alongside the drainage system), the following equation is used herein:
A D P = R D L R D L R R D S R R D S R
where R D L R and R D S R represent the inundation depth following the largest and smallest rainfalls that occurred in a rainfall time series or a certain period, while D L R and R D S R are the rainfall depths for the largest and smallest rainfall events in the time series or a certain period, respectively. The ADP index is expected to be a number less than one. An ADP closer to zero represents a highly adaptive system. For addressing the social aspect for all indicators, 60 mm of inundation depth is ignored during the calculation, as expressed in the previous sections.

4.1.8. Calculating the System Flexibility and Resilience

After estimating the quantities mentioned above, the system flexibility is calculated by averaging the normalized values of its components. It is worth noting that among all flexibility components, the system degradation capacity is inversely proportional to flexibility; therefore, before using this component in calculating the final flexibility index, it should be modified as Equation (18).
D C m = 1 D C
where D C m is the modified degradation capacity.
The Flexibility Index (FLX) is then obtained using Equation (19) for both single-event and long-term events in both hydraulic and social dimensions.
F L X = D C m + B r C + S R V + R P D + I M P + A D P / 6
Afterward, having the flexibility value and the resistance capacity, the UDS resilience index (RI) can be obtained by averaging the two components (Equation (20)), which could be interpreted as Table S4 (Supplementary Material). The classification of resiliency levels is suggested with equal intervals as proposed in some similar previous studies [43].
R I = F L X + R C / 2

4.2. Evaluating Different Scenarios for Improving the UDS Resiliency

The index-based methodology was applied to the WFD catchment using four types of drainage scenarios, including the existing condition of UDS, installing the BMPs, increasing the system capacity, and adding an auxiliary channel to a section of UDS where most inundation problems usually occur.

4.2.1. Scenario No. 1: The Existing Condition of UDS

The first scenario is considered the current situation of UDS without any changes to the existing system or using any adaptive measures. This scenario is required to be included to compare the degree of improvement achieved by other scenarios to see if they are helpful in comparison with the existing situation.
A primary simulation of UDS and rainfall-runoff modeling represent that the system encounters some flooded nodes and surcharged conduits when it goes under a 10-year, 6-h design rainfall. Figure S3 (Supplementary Material) shows the simulated drainage network in PCSWMM.

4.2.2. Scenario No. 2: Increasing the UDS Capacity

A proposed scenario could be the enlargement of the drainage channels in critical points to make up for the lack of hydraulic capacity in conveying the stormwater runoff. To apply this scenario, the flow conveyance capacity of the channel was increased up to 30 percent.

4.2.3. Scenario No. 3: Using Auxiliary Channels (Parallel Channels)

Given the discharge of all drainage channels of the area to the WFD channel in the south of the catchment, the probability of emerging bottlenecks alongside the main WFD channel is high [44]. Therefore, a scenario of adding a rectangular closed-channel to the location indicated in Figure S3 (Supplementary Material) is proposed. The auxiliary channel’s depth, width, and length are 2 m, 3 m, and 5 km, respectively. By adding this channel, it is expected that the system redundancy will increase, and the inundation points and surcharges from the system will decrease as a result.

4.2.4. Scenario No. 4: BMP Implementation alongside the Drainage System

The BMP installation scenario includes two sub-scenarios of (1) installing BMPs merely in the urban sub-basin and (2) installing BMPs in urban and upstream suburban areas. Two types of BMPs were determined for the catchment- i.e., Porous Pavement (PP) and Detention Pond (DP) according to the previous studies performed for the basin and some areas in the neighborhood [39,45]. Tahmasebi Birgani et al. [45] ranked eight BMPs according to the managers’ preferences in four areas of hydraulic, environmental, social, and economic aspects for an urban part of northern Tehran. The results (Table S5, Supplementary Material) indicated that some BMPs, including PP and DP, are among the most suitable for the studied area in various priorities considering the four studied aspects. Both selected BMPs were simulated in urban parts of the catchment. However, since PP is mostly suitable for urbanized areas, only DP was considered for the suburban areas of the basin.

4.3. Parametric-Based Resilience Calculation Using an MCDM Method

The second proposed method for calculating the system resilience index is a parametric MCDM methodology, namely IMPRESS (i.e., an Integrated Methodology for Parametric-based quantification of RESilience). In this approach, different characteristics of the system which play a role in accelerating/alleviating the urban flooding and inundation problems are first identified (i.e., the number and severity of rainfalls and inundations, average annual precipitation, UDS materials’ deterioration and aging, the drainage system density per unit area, etc.). These characteristics were chosen considering both single-event and continuous rainfalls/flooding. In addition, some characteristics of resilience (some accepted from previous works and some are proposed here in this study) are considered, including robustness, amplitude (i.e., the magnitude of the system response to a disturbance), graduality (i.e., the response trend), redundancy (i.e., representing an excess capacity and the extent to which system components are substitutable), modularity (i.e., the ability of each part of the system to perform independently and on its own), survivability (i.e., the degree of agility for bouncing back towards the previous functionality after reaching the minimum failure state), and rapidity (i.e., the speed of recovery in a disturbed system).
In the next step, the actual values (vobs) for each system’s characteristics are acquired based on several databases, including historical data, official reports, and data recorded in the Iranian Meteorological Organization (IMO), Iranian Water Resources Management Company, Iranian Red Crescent and Disaster Management Organizations, Environmental Conservation agencies, Tehran Municipality as well as doing interviews and field observations. In addition, the corresponding minimum and maximum values (Vmin and Vmax) are furthermore obtained according to the available data and observations.
Since the identified system characteristics do not have the same importance regarding resilience to floods, the relative importance of each factor should be determined by a weight given to it by experts. The actual values and weights were normalized to be used as a non-dimensional number between 0 and 1.
In the next stage, a matrix approach is utilized, and the resilience matrix is developed. For this purpose, the system characteristics are categorized under relevant domains and resilience properties, as indicated in Figure 5. Finally, PROMETHEE (i.e., preference ranking organization method for Enrichment of Evaluations) was utilized to quantify the resilience index based on an MCDM approach. (For more information on the PROMETHEE method, refer to Brans et al. [46] and Brans and Vincke [47]).

5. Results and Discussion

5.1. Data Collection

The map of the studied region, as well as the UDS configuration and the conduits’ characteristics, were obtained from Yekom Consulting Engineers [48] and some past studies performed on the basin [49]. The data required for rainfall-runoff modeling—such as sub-basins’ curve number (CN), manning coefficient, UDS characteristics, and land properties—were extracted from relevant sources [33,50,51]. Wherever needed, information on the catchment was gained through recorded data in the Tehran Municipality Database and interviews with informed individuals and authorities.
Rainfall data was mainly extracted from two stations inside the WFD catchment, Mehrabad and Darakeh (Haft-Howz), whose characteristics are listed in Table 2. For long-term precipitation as the disturbance, a 30-year rainfall time series from 1985 to 2015 was utilized.

5.2. Estimating the Resilience Indicators

To estimate UDS robustness to a 6-h rainfall, it is first required that total rainfall depth is calculated through IDF curves obtained for the studied region. Figure S4 (Supplementary Material) shows the catchment’s IDF curves for different return periods, and Table S6 (Supplementary Material) represents the values of total rainfall depths for different return periods of a 6-h rainfall.
By introducing the total rainfalls for any return periods presented in Table S6 (Supplementary Material) to the rainfall-runoff model, one can have the stormwater runoff surcharged from the UDS. Table 3 has listed the average runoff depth surcharged from the whole drainage network in the face of rainfalls with different return periods while using various scenarios. Since stormwater runoff depth directly correlates with the consequent damages [10,52,53,54,55,56], it seems to be a good indicator representing the UDS response to the disturbance (i.e., extreme rainfall). As can be seen in Table 3, the Response Magnitude (RM) values for the current situation of UDS in the study region are significantly larger than RM for the proposed adaptive measures for both single-event and continuous rainfall. Among all scenarios, BMP implementation in both urban and non-urban parts of the basin reduces the average surcharged runoff from the UDS the most.
Robustness, which is representative of UDS resistance capacity, was estimated according to what was explained in the Methods section and is similar to the previous works [10]. In this viewpoint, the threshold to which the drainage network can bear the disturbance- no surcharge occurs- is considered as the reaction threshold, shown by a unique number, which is the rainfall value for the corresponding return period (presented in Table S6, Supplementary Material). For instance, using the first scenario (channel enlargement), a 2-year, 6-h rainfall can be tolerated, whose value is 32.62 mm according to Table S6 (Supplementary Material). For the second scenario, the resistance capacity is the same. BMP installation has the most impact on improving the system’s robustness, particularly if BMPs are installed in both urbanized and upstream suburban areas. This is while the resistance capacity for the current condition of UDS is nearly zero. On this basis, the values of UDS robustness to a 6-h rainfall could be represented as indicated in Table 4. Such values represent the capacity of adapted UDS to absorb the disturbance and withstand it.
A similar procedure was accomplished for a continuous rainfall time series in both social and hydraulic domains. The long-term system robustness was obtained by estimating the surcharges from the UDS and the rate of rainfalls in the time series the system can tolerate. It can be concluded that for both single-event and continuous rainfalls, the resistance domain would become larger under the given adaptive scenarios.
Based on the average degradation capacity obtained for UDS, the braking capacity is estimated using Equations (10) and (11). Table 5 indicates that the system braking capacity is the largest in the case of BMP installation and is the lowest for the existing condition of UDS, which leads to lower resilience.
To summarize the evaluation briefly and rapidly, the rest of the indicators are calculated only for the 100-year, 6-hr rainfall. A 100-year rainfall could represent an extreme event, and a 6-h rainfall is the typical duration of design rainfall in the study area. The calculation results have been presented in Table 6 and Table S7 (Supplementary Material) for the hydraulic and social domains, respectively. It can be observed that nearly in all cases, BMP installation in upstream non-urban areas as long as the urbanized parts of the basin would be the best adaptive scenario, both for single-event and continuous rainfall. Absorptive capacity, as well as survivability and rapidity, are the least for the current state of UDS while possessing the highest values for BMP implementation scenarios. Among two other management scenarios, channel enlargement is mostly more successful than using an auxiliary channel. This may be due to the local effectiveness of the auxiliary channel while inundated locations are expanded through the entire drainage system. Regarding social resilience, according to the obtained results, there is no serious problem even with the existing condition of the drainage system in the studied area. However, even in the social context, utilizing the management strategies could help improve system sustainability.
It should be noted that to estimate the total flexibility and resilience, the calculated indices would be normalized considering the maximum possible value for the given indicator.

5.3. Final Resilience of the System

Table 7 represents the results of estimating UDS future resilience calculated under various scenarios for hydraulic and social aspects and two types of disturbances, i.e., single-event and continuous rainfall. The table shows that in the system in both single event and long-term rainfall cases, social resistance values are higher than hydraulic ones, resulting from the methodology that considers a 60 mm surcharged runoff in the social domain the same as zero runoff in the hydraulic domain. In addition, for continuous rainfall time series, the resilience values—as well as the values of its components—are near the ones for single-event rainfall, though a little bit smaller, which again is a consequence of the proposed methodology and formula for gaining the value of such indicators.
According to the results and considering Table S4 (Supplementary Material), the UDS resilience for the studied area in its current condition is evaluated as “low” except for social resilience, which hardly enters the “medium” level. Using adaptive scenarios, the resilience level of the drainage system would improve significantly (Table 7). So that by implementing the scenarios of “adding an auxiliary channel to the system”, the resilience level falls into one level higher, i.e., “medium level”. For both scenarios of “enlarging the UDS channels” and “BMP installation only in the urban part of the catchment”, the UDS resilience would mostly reach the “acceptable level”; and the scenario of “BMP installation in both urban and upstream non-urban areas” would result in a system with “high” resilience. Such findings indicate the importance of utilizing various adaptive measures in urbanized catchments and would help urban planners find the best solutions for reducing urban inundation problems. The results are in agreement with some similar previous works, including the one performed by Tahmasebi Birgani and Yazdandoot [10], which was accomplished for a region near the WFD catchment, and emphasizes the fact that, considering most of the resilience indicators, using BMPs are more effective than channels’ enlargement for improving the UDS resiliency to inundation and flooding.

5.4. Comparison of the Resilience Calculated from the two Methodologies

Table 8 represents the resilience calculated for social and hydraulic aspects of the drainage system by the IMPRESS method for the existing condition of the UDS. Comparing the results with the index-based method shows that the calculated index from both methodologies is almost the same and relatively low (below 50%) for the studied region. Therefore, using adaptive measures such as BMPs (which, in the previous section, were proved to be the best solutions for the studied region) seems to be of great necessity in this semi-urbanized basin.
According to the results of the IMPRESS framework, the social resilience index shows a lower value than the hydraulic resilience index. This is why when using the index-based methodology the social aspects have an almost better condition compared to hydraulic issues of the drainage system. It can be related to the formulas used for calculating the social resilience, considering the water depths below 6 cm as no threat to social disorder and community safety. This is while the calculated index in the parameter-based method is based on previous events, observations, and interviews, which is more probable to be coincident with reality.

6. Conclusions

Unlike the traditional risk management of urban critical infrastructures, which makes an effort to resist the system in the face of all types of perturbations (as much as possible), the resilience-based approach emphasizes the system’s ability to restore from the unexpected disturbances given the fact that not all disasters can be predicted or even recognized. The flexibility of a system has recently gained attention as an important factor related to resilience. However, their relative position to each other is still a matter of discussion. The current study focused mainly on providing a framework for quantifying the UDS resilience to disasters (flooding, in particular) based on a main resilience component, i.e., flexibility, and aimed to cover the gap in previous research by accounting for recovery from long-term precipitation considering both hydraulic and social aspects.
The resilience of a drainage system in a semi-urbanized basin was calculated considering four different scenarios to find the best solution for improving the UDS functionality and resilience to flooding. Within this proposed framework, seven indicators representing the system behavior in the face of a disaster were quantified. The resilience indices were obtained and discussed from both hydraulic and social viewpoints under both single-event and continuous precipitation, and compared with the results from an MCDM approach. The findings of this study suggested a relatively low level of flood resiliency for the study site, and it was revealed that BMPs as structural adaptive solutions are among the best choices for improving the system’s resilience to flooding.
By focusing on resilience measurement to quantify its components, this study can be considered the first step to the sustainable management and planning of the drainage infrastructure in urban areas. Resilience quantification for UDSs is a crucial task for evaluating the system as it can provide a better view of the direction the system functionality is moving toward. This could help the urban managers and practitioners evaluate the performance of such critical infrastructures as a primary step to managing them more efficiently and making more realistic decisions.
In the current study, the social resilience indicators were mostly obtained considering the fact that a 6 cm inundation depth is tolerable for the residents and commuters in the region (based on the information from reliable sources). This amount of inundation depth may differ for other societies based on other sources of information. Besides, in future works, other social parameters could be added to such a methodology to evaluate the social resilience of UDS to flooding in the context of the complicated urban society. The proposed framework in the current study can be utilized in predicting the UDS resilience in future periods and, therefore, can be used as a method to estimate the design criteria of UDS in the face of future climatic changes and urban development.
However, performing this methodology has its limitations and difficulties. Gathering data from reliable sources is the most important challenge for having accurate results, and this can be more challenging when using this methodology for estimating infrastructure resilience in the future. Nevertheless, if reliable data are available, this proposed framework can act as a helpful tool for the primary assessment of UDS performance which covers the gaps in some similar tools developed in the past. The other issue, which is suggested to be considered in similar studies, is the need to examine if the final solutions to solve the UDS problems are economically and environmentally feasible and sustainable. In the current study, the focus was mainly on quantifying social and hydraulic aspects of UDS resilience. However, a cost-benefit analysis could be added to identify the best compromise between increasing resilience and investments based on the probability of rainfall events and the impacts on UDS. It should be emphasized that the proposed methodology and the indicators could be modified to be usable for resilience assessment from other viewpoints, such as economic and environmental aspects.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/xxx/s1:Figure S1: Graduality vs. braking capacity; Figure S2: The difference in systems’ behavior under a similar disturbance: (A) resistance capacity, (B) response trend and (C) system recovery; Figure S3: The WFD basin simulated by PCSWMM and the location of proposed auxiliary channel; Figure S4: Calculated IDF curves for the studied region; Table S1: Some key parameters of SWMM model (for Sensitivity Analysis) and their Permissible range of change; Table S2: The calibrated values of some important parameters of SWMM model; Table S3: The results of rainfall-runoff SWMM verification; Table S4: Evaluating the calculated adaptability index; Table S5: The selected BMPs’ ranking according to the manager’s preferences for the northern Tehran; Table S6: Total rainfall depth for different return periods; Table S7: The quantified (normalized) sub-factors of UDS flexibility for various scenarios for the social aspect [57,58,59].

Author Contributions

Conceptualization, methodology, simulation, validation, formal analysis, investigation, data gathering, writing and original draft preparation, N.B.; Supervision and project administration, M.H.N. and A.S.; Review and editing, W.R.; Review and funding acquisition, G.T.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tahmasebi Birgani, Y. Sustainable Management of Urban Drainage Based on Flexible Risk Management. Ph.D. Thesis, K.N. Toosi University of Technology, Tehran, Iran, 2014. [Google Scholar]
  2. Restemeyer, B.; Woltjer, J.; Van Den Brink, M. A strategy-based framework for assessing the flood resilience of cities—A Hamburg case study. Plan. Theory Pract. 2015, 16, 45–62. [Google Scholar] [CrossRef]
  3. Tempels, B.; Hartmann, T. A co-evolving frontier between land and water: Dilemmas of flexibility versus robustness in flood risk management. Water Int. 2014, 39, 872–883. [Google Scholar] [CrossRef]
  4. McClymont, K.; Morrison, D.; Beevers, L.; Carmen, E. Flood resilience: A systematic review. J. Environ. Plan. Manag. 2020, 63, 1151–1176. [Google Scholar] [CrossRef] [Green Version]
  5. Asokan, V.A.; Yarime, M.; Esteban, M. Introducing Flexibility to Complex, Resilient Socio-Ecological Systems: A Comparative Analysis of Economics, Flexible Manufacturing Systems, Evolutionary Biology, and Supply Chain Management. Sustainability 2017, 9, 1091. [Google Scholar] [CrossRef] [Green Version]
  6. Husdal, J. Flexibility and Robustness as Options to Reduce Risk and Uncertainty; Unpublished Draft Paper; Molde University College: Molde, Norway, 2004; Available online: http://husdal.com/2004/12/15/flexibility-and-robustness-as-options-to-reduce-risk-and-uncertainty/ (accessed on 6 August 2018).
  7. Stinner, S.; Huchtemann, K.; Muller, D. Flexibility Quantification for building energy systems with heat pumps. In Proceedings of the 15th IBPSA Conference (Building Simulation 2017), San Francisco, CA, USA, 7–9 August 2017. [Google Scholar]
  8. Fletcher, S.M. Learning and Flexibility for Water Supply Infrastructure Planning under Diverse Uncertainties. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2018. [Google Scholar]
  9. Sidqi, Y.; Ferrez, P.; Gabioud, D.; Roduit, P. Flexibility quantification in households: A swiss case study. Energy Inform. 2020, 3, 23. [Google Scholar] [CrossRef]
  10. Birgani, Y.T.; Yazdandoost, F. A framework for evaluating the persistence of urban drainage risk management systems. J. Hydro-Environ. Res. 2014, 8, 330–342. [Google Scholar] [CrossRef]
  11. Mugume, S.N.; Butler, D. Evaluation of functional resilience in urban drainage and flood management systems using a global analysis approach. Urban Water J. 2016, 14, 727–736. [Google Scholar] [CrossRef]
  12. Bakhshipour, A.E.; Dittmer, U.; Haghighi, A.; Nowak, W. Toward Sustainable Urban Drainage Infrastructure Planning: A Combined Multiobjective Optimization and Multicriteria Decision-Making Platform. J. Water Resour. Plan. Manag. 2021, 147, 04021049. [Google Scholar] [CrossRef]
  13. Arup (International Development). City Resilience Framework, The Rockerfeller Foundation, 2015. p. 5. Available online: https://assets.rockefellerfoundation.org/app/uploads/20160105134829/100RC-City-Resilience-Framework.pdf (accessed on 29 October 2018).
  14. Lioyd’s and Arup. Future Cities: Building infrastructure resilience, City Infrastructure Resilience–Designing the Future, Emerging Risk Report: Society and Security. 2017, p. 66. Available online: https://assets.lloyds.com/assets/pdf-arup-lloyds-arup-future-cities-2017/1/pdf-arup-lloyds-arup-future-cities-2017.pdf (accessed on 16 August 2018).
  15. Ku, A. Modelling Uncertainty in Electricity Capacity Planning. Ph.D. Thesis, London Business School, London, UK, 1995. [Google Scholar]
  16. Walker, W.E.; Haasnoot, M.; Kwakkel, J.H. Adapt or Perish: A Review of Planning Approaches for Adaptation under Deep Uncertainty. Sustainability 2013, 5, 955–979. [Google Scholar] [CrossRef] [Green Version]
  17. De Bruijn, K. Resilience and flood risk management. Water Policy 2004, 6, 53–66. [Google Scholar] [CrossRef]
  18. De Bruijn, K.M. Resilience indicators for flood risk management systems of lowland rivers. Int. J. River Basin Manag. 2004, 2, 199–210. [Google Scholar] [CrossRef]
  19. Heggarty, T.; Bourmaud, J.-Y.; Girard, R.; Kariniotakis, G. Quantifying power system flexibility provision. Appl. Energy 2020, 279, 115852. [Google Scholar] [CrossRef]
  20. De Bruijn, K.M. Resilience and flood risk management. A systems approach applied to lowland rivers. Ph.D. Thesis, Technical University Delft, Delft, The Netherlands, 2005. [Google Scholar]
  21. Gersonius, B. Can resilience support integrated approaches to urban drainage management? In Proceedings of the 11th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 31 August–5 September 2008. [Google Scholar]
  22. Zevenbergen, C. Flood resilience. In IRGC Resource Guide on Resilience; Resilience Guide Edition; Linkov, I., Florin, M.-V., Eds.; International Riks Governance Council (IRGC), EPFL International Risk Governance Center: Lausanne, France, 2016; Available online: https://www.irgc.org/risk-governance/resilience/ (accessed on 29 October 2018).
  23. Mens, M.J.P.; Klijn, F.; de Bruijn, K.M.; van Beek, E. The meaning of system robustness for flood risk management. Environ. Sci. Policy 2011, 14, 1121–1131. [Google Scholar] [CrossRef]
  24. Butler, D. Crowdsourcing goes mainstream in typhoon response. Nature 2013. [Google Scholar] [CrossRef]
  25. Moghimi, E.; Saffari, A. Geomorphologic Assessment of Urban Development in the Territory surface of Drainage Basins Case Study: Tehran Metropolis. Hum. Sci. MODARES J. (MJSP) 2010, 14, 1–31. [Google Scholar]
  26. Rossman, L.A. Stormwater Management Model (SWMM); USEPA: Washington, DC, USA, 2017. [Google Scholar]
  27. Simon, M. An Introduction to the Stormwater Management Model (SWMM); Surface Water Quality Modeling: Cincinnati, OH, USA, 2018. [Google Scholar]
  28. Ghahreman, B. Updated IDF equation for rainfall in Iran using 1-hour, 10-year rainfall. J. Agric. Sci. 1996, 6, 13–30. (In Persian) [Google Scholar]
  29. Bell, F.C. Generalized Rainfall-Duration-Frequency Relationship. J. Hydraul. Division. ASCE 1957, 95, 6357. [Google Scholar] [CrossRef]
  30. Alizadeh, A. Principles of Applied Hydrology, 31st ed.; University of Imam Reza (Peace be upon him): Mah’had, Iran, 2010; p. 735. (In Persian) [Google Scholar]
  31. Ghahreman, B.; Sepaskhah, A. Estimating IDF equation for rainfall in Iran using 1-hour, 10-year rainfall. In 3rd International Congress on Road and Structural Engineering; Faculty of Engineering, University of Shiraz: Shiraz, Iran, 1980. (In Persian) [Google Scholar]
  32. Bakhtiari, B.; Pour Mousavi, S.H.; Sayyari, N. Investigating the effect of climate change on IDF curves for rainfalls in Babolsar station during 2011–2030. Iran. J. Irrig. Drain. 2014, 4, 694–704. (In Persian) [Google Scholar]
  33. MGCE and Poyry. Tehran Stormwater Management Master Plan; MGCE: Tehran, Iran, 2011. [Google Scholar]
  34. Holling, C.; Meffe, G.K. Command and Control and the Pathology of Natural Resource Management. Conserv. Biol. 1996, 10, 328–337. [Google Scholar] [CrossRef] [Green Version]
  35. Furuta, K. Resilience engineering. In Reflections on the Fukushima Daiichi Nuclear Accident; Ahn, J., Carlson, C., Jensen, M., Juraku, K., Nagasaki, S., Tanaka, S., Eds.; Springer Open: Berlin/Heildelberg, Germany, 2014; pp. 435–454. [Google Scholar]
  36. Nazif, S. Developing an Algorithm for Climate Change Assessment on Urban Water Cycle; University of Tehran: Tehran, Iran, 2010. [Google Scholar]
  37. Iranian Red Crescent Society (IRCS). Actions during Flood and after That, 2019. Available online: http://rcs.ir (accessed on 29 October 2018).
  38. Zhou, Q.; Mikkelsen, P.S.; Halsnæs, K.; Arnbjerg-Nielsen, K. Framework for economic pluvial flood risk assessment considering climate change effects and adaptation benefits. J. Hydrol. 2012, 414–415, 539–549. [Google Scholar] [CrossRef]
  39. Binesh, N.; Niksokhan, M.H.; Sarang, A.; Rauch, W. Improving sustainability of urban drainage systems for climate change adaptation using best management practices: A case study of Tehran, Iran. Hydrol. Sci. J. 2019, 64, 381–404. [Google Scholar] [CrossRef]
  40. Rogers, L.K.; Lucchesi, P.A. Stress adaptation and the resilience of youth: Fact or Fiction? Physiology 2014, 29, 156. [Google Scholar] [CrossRef] [Green Version]
  41. Carmichael, D.G. Incorporating resilience through adaptability and flexibility. Civ. Eng. Environ. Syst. 2015, 32, 31–43. [Google Scholar] [CrossRef]
  42. Sideroff, S. The Path: Mastering the Nine Pillars of Resilience and Success; Third Wind Press: Los Angeles, CA, USA, 2015. [Google Scholar]
  43. Karamouz, M.; Zeynolabedin, A.; Olyaei, M.A. Regional Drought Resiliency and Vulnerability. J. Hydrol. Eng. 2016, 21. [Google Scholar] [CrossRef]
  44. Yazdi, J. Improving Urban Drainage Systems Resiliency Against Unexpected Blockages: A Probabilistic Approach. Water Resour. Manag. 2018, 32, 4561–4573. [Google Scholar] [CrossRef]
  45. Birgani, Y.T.; Yazdandoost, F.; Moghadam, M. Role of Resilience in Sustainable Urban Stormwater Management. J. Hydraul. Struct. Shahid Chamran Univ. 2013, 1, 44–53. [Google Scholar] [CrossRef]
  46. Brans, J.P.; Mareschal, B.; Vincke, P. PROMETHEE: A New Family of Outranking Methods in MCDM; IFORS 84: Amsterdam, The Netherlands, 1984; pp. 477–490. [Google Scholar]
  47. Brans, J.P.; Vincke, P. A Preference Ranking Organisation Method: (The PROMETHEE Method for Multiple Criteria Decision-Making). J. Manag. Sci. 1985, 31, 647–656. [Google Scholar] [CrossRef] [Green Version]
  48. YCE (Yekom Consulting Engineers). A Project on (Quantitative & Qualitative) Rehabilitation of Surface Waters in Southern Tehran; Yekom Consulting Engineering Company: Tehran, Iran, 2009; Volume 3. (In Persian) [Google Scholar]
  49. Moafi Rabori, A. Optimal Design of WFD (West Flood-Diversion) Dimensions based on Upland Catchment’ Characteristics. Master’s Thesis, University of Tehran, Tehran, Iran, 2012. [Google Scholar]
  50. McCuen, R.H. A Regional Approach to Urban Stormwater Detention. Geophys. Res. Lett. 1974, 1, 321–322. [Google Scholar] [CrossRef]
  51. Parsa, V.; Moti’ee, H. Modeling urban floods using Storm-Cad. In Proceedings of the 5th Conference on Iran Water Resources Management, Tehran, Iran, 17–18 February 2014. [Google Scholar]
  52. Dutta, D.; Herath, S.; Musiake, K. A mathematical model for flood loss estimation. J. Hydrol. 2003, 277, 24–49. [Google Scholar] [CrossRef]
  53. Penning-Rowsell, E.; Johnson, C.; Tunstall, S.; Tapsell, S.; Morris, J.; Chatterton, J.; Coker, A.; Green, C. The Benefits of Flood and Coastal Defence, Techniques and Data 2002; Flood Hazard Research Centre, Middlesex University: London, UK, 2003. [Google Scholar]
  54. FLOODsite. Flood Risk Assessment and Flood Risk Management. In An Introduction and Guidance Based on Experiences and Findings of FLOODsite (An EU-Funded Integrated Project); Deltares/Delft Hydraulics: Delft, The Netherlands, 2009; Available online: http://www.floodsite.net/html/partner_area/project_docs/T29_09_01_Guidance_Screen_Version_D29_1_v2_0_P02.pdf (accessed on 29 October 2018).
  55. Jonkman, S.; Bočkarjova, M.; Kok, M.; Bernardini, P. Integrated hydrodynamic and economic modelling of flood damage in the Netherlands. Ecol. Econ. 2008, 66, 77–90. [Google Scholar] [CrossRef]
  56. Kazama, S.; Sato, A.; Kawagoe, S. Evaluating the Cost of Flood Damage Based on Changes in Extreme Rainfall in Japan. Sustain. Sci. 2010, 6, 23–27. [Google Scholar] [CrossRef]
  57. Javaheri Seyyed, M.A. Robustness Assessment of Optimal Management Options for Improving the Performance of Urban Stormwater Drainage Systems. Master’s Thesis, University of Tehran, Tehran, Iran, 2015. [Google Scholar]
  58. Lee, M. Resilience Assessment in Geotechnical Engineering. Ph.D. Thesis, University of Waterloo, Waterloo, ON, Canada, 2016. [Google Scholar]
  59. Harte, D.; Williams, A.; Turner, J. Reciprocity and the hyper-local journalist. Journal. Pract. 2017, 11, 160–176. [Google Scholar]
Figure 1. The system behavior in response to a disturbance.
Figure 1. The system behavior in response to a disturbance.
Water 14 02007 g001
Figure 2. The proposed framework for quantifying system resiliency to disasters.
Figure 2. The proposed framework for quantifying system resiliency to disasters.
Water 14 02007 g002
Figure 3. The situation of WFD catchment in Tehran (province) and Iran.
Figure 3. The situation of WFD catchment in Tehran (province) and Iran.
Water 14 02007 g003
Figure 4. Flowchart of index-based measurement of flexibility and resilience.
Figure 4. Flowchart of index-based measurement of flexibility and resilience.
Water 14 02007 g004
Figure 5. The flowchart of resilience quantification in a parametric-based method.
Figure 5. The flowchart of resilience quantification in a parametric-based method.
Water 14 02007 g005
Table 1. The sub-basins’ physiographic characteristics in WFD catchment [25].
Table 1. The sub-basins’ physiographic characteristics in WFD catchment [25].
The ParameterArea (Km2)Perimeter (m)Length (m)Min height (m)Max height (m)Time of Concentration (min)Average Slope (%)
The Sub-Basin
Hesarak (Suburb)1417,5167052155027503550
Hesarak (Urban)Main Channel814,763.7684012431570565.3
Shahin & Shaghayegh Channels1724,942848412431985537.2
Farahzad (Suburb)2324,412.79909160034004557.4
Farahzad (Urban)Main Channel2027,337.6945012692224549.4
Moradabad channel39141.34027156525202036.9
Behroud channel410,833.54887162827902336.8
Darakeh (Suburb)2524,121.89376170038003959
Darakeh (Urban)Main channel1229,472.79384132022005513.2
Tappeh-Neizar channel1319,239.863713001800547.9
Khoshkeh channel29151.33750165027001846.7
Table 2. The characteristics of two main stations in the basin for rainfall data extraction.
Table 2. The characteristics of two main stations in the basin for rainfall data extraction.
Station NameStation TypeLatitudeLongitudeHeight (MASL)
Darakeh (Haft-howz)Hydrometer35°48′ N51°23′ E1700
MehrabadSynoptic35°41′ N51°19′ E1190.8
Table 3. Average surcharge (mm) from the UDS (i.e., degradation capacity) calculated for various scenarios.
Table 3. Average surcharge (mm) from the UDS (i.e., degradation capacity) calculated for various scenarios.
Return Period2-year5-year10-year25-year50-year100-yearContinuous Rainfall
Scenario
Channel enlargement0102339589219
Auxiliary channel01529476610323
BMP in urban areas0081625339
BMP in both urban and non-urban sub-basins000411193
Existing UDS1637496812933754
Table 4. The value of UDS resistance capacity (RC) under different scenarios.
Table 4. The value of UDS resistance capacity (RC) under different scenarios.
ScenarioSingle-Event RainfallContinuous Rainfall
Hydraulic RobustnessSocial RobustnessHydraulic RobustnessSocial Robustness
Channel enlargement32.6278.911.423.2
Auxiliary channel32.6263.597.3116.6
BMP in urban areas44.96(more than) 78.932.549.4
BMP in both urban and non-urban sub-basins53.2(more than) 78.944.876
Existing UDS053.21.889.34
Table 5. Braking capacity (BrC) for various scenarios (mm/min).
Table 5. Braking capacity (BrC) for various scenarios (mm/min).
Return Period2-year5-year10-year25-year50-year100-yearContinuous Rainfall
Scenario
Channel enlargement-4.82.30.940.780.344.4
Auxiliary channel-41.60.70.50.23.9
BMP in urban areas--5.34.93.32.85.7
BMP in both urban and non-urban sub-basins---7.46.95.66.9
Existing UDS82.91.10.550.360.112.6
Table 6. The quantified (normalized) sub-factors of UDS flexibility for various scenarios for the hydraulic aspect.
Table 6. The quantified (normalized) sub-factors of UDS flexibility for various scenarios for the hydraulic aspect.
ScenariosType of DisturbanceAbCSRVRPDIMPADP
Channel enlargementSingle-event rainfall32.970.510.160.540.54
Continuous rainfall15.850.480.150.490.5
auxiliary channelSingle-event rainfall32.830.380.110.510.53
Continuous rainfall11.250.370.090.450.47
BMP in urban areasSingle-event rainfall47.79304.10.790.77
Continuous rainfall38.31263.80.740.75
BMP in both urban and non-urban sub-basinsSingle-event rainfall58.85400070.830.84
Continuous rainfall52.03906.70.80.79
Existing condition of UDSSingle-event rainfall0.110.190.04--
Continuous rainfall4.500.170.03--
Table 7. The calculated resilience (percent) by different scenarios.
Table 7. The calculated resilience (percent) by different scenarios.
ResilienceSingle-Event RainfallContinuous Rainfall
Scenario Hydraulic ResilienceSocial ResilienceHydraulic ResilienceSocial Resilience
Enlarging the channels0.640.660.610.62
Adding an auxiliary channel0.530.590.50.55
BMPs in Urban sub-basin0.790.830.750.78
BMPs in both Urban and Suburban areas0.870.920.830.89
Existing UDS0.390.440.350.39
Table 8. The calculated resilience index by IMPRESS methodology for various resilience domains.
Table 8. The calculated resilience index by IMPRESS methodology for various resilience domains.
Resilience DomainsSingle-Event RainfallContinuous Rainfall
Scenario Hydraulic ResilienceSocial ResilienceHydraulic ResilienceSocial Resilience
Existing UDS0.410.390.370.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Binesh, N.; Sarang, A.; Niksokhan, M.H.; Rauch, W.; Aronica, G.T. Quantifying the UDS Hydraulic and Social Resilience to Flooding: An Index-Based Approach vs. a Parameter-Based MCDM Method. Water 2022, 14, 2007. https://doi.org/10.3390/w14132007

AMA Style

Binesh N, Sarang A, Niksokhan MH, Rauch W, Aronica GT. Quantifying the UDS Hydraulic and Social Resilience to Flooding: An Index-Based Approach vs. a Parameter-Based MCDM Method. Water. 2022; 14(13):2007. https://doi.org/10.3390/w14132007

Chicago/Turabian Style

Binesh, Negin, Amin Sarang, Mohammad Hossein Niksokhan, Wolfgang Rauch, and Giuseppe T. Aronica. 2022. "Quantifying the UDS Hydraulic and Social Resilience to Flooding: An Index-Based Approach vs. a Parameter-Based MCDM Method" Water 14, no. 13: 2007. https://doi.org/10.3390/w14132007

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