Application of GIS-Interval Rough AHP Methodology for Flood Hazard Mapping in Urban Areas

Floods are natural disasters with significant socio-economic consequences. Urban areas with uncontrolled urban development, rapid population growth, an unregulated municipal system and an unplanned change of land use belong to the highly sensitive areas where floods cause devastating economic and social losses. The aim of this paper is to present a reliable GIS multi-criteria methodology for hazard zones’ mapping of flood-prone areas in urban areas. The proposed methodology is based on the combined application of geographical information systems (GIS) and multi-criteria decision analysis (MCDA). The methodology considers six factors that are relevant to the hazard of flooding in urban areas: the height, slope, distance to the sewage network, the distance from the water surface, the water table and land use. The expert evaluation takes into account the nature and severity of observed criteria, and it is tested using three scenarios: the modalities of the analytic hierarchy process (AHP). The first of them uses a new approach to the exploitation of uncertainty in the application of the AHP technique, the interval rough numbers (IR’AHP). The second one uses the fuzzy technique for the exploitation of uncertainty with the AHP method (F’AHP), and the third scenario contemplates the use of the traditional (crisp) AHP method. The proposed methodology is demonstrated in Palilula Municipality, Belgrade, Serbia. In the last few decades, Palilula Municipality has been repeatedly devastated by extreme flood events. These floods severely affected the transportation networks and other infrastructure. Historical flood inundation data have been used in the validation process. The final urban flood hazard map proves a satisfactory agreement between the flood hazard zones and the spatial distribution of historical floods that happened in the last 58 years. The results indicate that the scenario in which the IR’AHP methodology is used provides the highest level of compatibility with historical data on floods. The produced map showed that the areas of very high flood hazard are located on the left Danube River bank. These areas are characterized by lowland morphology, gentle slope, sewage network, expansion of impermeable locations and intense urbanization. The proposed GIS-IR’AHP methodology and the results of this study provide a good basis for developing a system of flood hazard management in urban areas and can be successfully used for spatial city development policy.


Introduction
Floods fall into the most serious natural disasters in the world that endanger more lives and cause more property damage than any other natural phenomenon [1][2][3][4][5].Human activities, such as the growth of settlements and economic assets in flooding areas and the reduction of the area for natural water retention as a result of land use and climate change, contribute to increasing the probability of the adverse impacts of floods.Due to the absolute and relative population growth of urban areas, cities are physically wider, which means the expanding residential zone, the business, commercial and industrial areas and the road and railway network in the natural and agricultural land around cities. It usually involves the removal or reduction of vegetation and construction of impervious surfaces in the form of buildings, sidewalks, parking lots and roads.The negative effect of urbanization is two-fold: on the one hand, it accumulates objects threatened by flooding, while on the other hand, narrowing the watercourse profile and reducing the throughput of land.The spread of impervious surfaces is the main driver of hydrological change and leads to an increase and acceleration of runoff storm water [6][7][8][9][10].In underdeveloped areas, the development of cities and the changes that this development causes, usually are not accompanied by the construction of an adequate drainage system for the atmospheric water management.The classical approach based on the collection of rain water from the urban area sewage system and taking it through the fastest route to the closest recipient has been used instead.Most often, this approach results in inefficient solutions for the removal of atmospheric water, leading to increasingly frequent flooding.To overcome such a difficult situation, it is necessary to establish a systematic approach in order to define flood areas in urban zones as soon as possible [11].
Risk reduction and flood management in urban areas inhabited with a large population have multiple significances.In addition to reducing material and human losses, they reduce the uncontrolled contamination with pollutants of known or unknown origin in the urban environment, but also create conditions for optimal land use.Contemporary trends of integrated water management and planning in urban environments include adequate risk assessment of their appearance and implementation of a number of technical and preventive measures that enable the control of the movement of rainwater and waste water in all projected hydrological regimes and their subsequent draining through the drain and sewage system.Flood management cannot become technically controllable without a proper assessment of flood hazard mapping and flood hazard.
Flood mapping is a crucial element of flood risk management.Directive 2007/60/EC on the assessment and management of flood risks will require Member States to prepare two types of maps by 2013 [12]: Flood hazard maps, showing the extent and expected water depths/levels of an area flooded in three scenarios, a low probability scenario or extreme events, in a medium probability scenario (at least with a return period of 100 years) and, if appropriate, a high probability scenario.
Flood risk maps shall also be prepared for the areas flooded under these scenarios showing potential population, economic activities and the environment at potential risk from flooding and other information that Member States may find useful to include, for instance other sources of pollution.Accordingly, the first step is the hazard mapping, containing information on the flooding areas, which is the subject of the work.
In addition to the assessment, the flood hazard maps are useful tools for spatial planning and development of the city, especially in the area where it is necessary to identify infrastructure guidelines for the construction of a sewage drainage system and to determine the extent of the strict construction of new settlements [13].
Nowadays, the risk analysis of natural disasters is unimaginable without the digital backing and support of geographic information systems (GIS).Natural disasters are multidimensional phenomena with a spatial dimension, which makes GIS very applicative for such analysis [14][15][16][17][18][19][20][21][22][23].Geographic information systems are suitable for this type of study because the powerful geostatistical tool can effectively manage large volumes of spatial data, and many studies dealing with flood risks and flood modeling effectively use GIS [24][25][26].
However, the application of GIS in the process of decision-making has certain shortcomings, which are reflected in the inability of multiple criteria analysis, which affects the decision-making.
In order to enable multi-criteria decision-making and presentation of the results in the area using GIS, it is necessary to combine tools for multi-criteria decision analysis (MCDA) with GIS.The application of MCDA in the GIS environment is currently a fundamental tool to solve problems in many areas, because it allows greater flexibility and accuracy in decision-making [27][28][29][30][31][32].
Over the last few decades, the combination of GIS and MCDA has been successful, and in a number of studies, it is often used in the risk assessment and generation of flood risk maps [8,15,[33][34][35][36][37][38][39][40].The researchers have used different methods and models for mapping of natural hazards using the GIS techniques [41][42][43][44][45].A number of studies show that AHP in the GIS environment is the most popular and powerful method for generating flood hazard maps, with a good degree of accuracy, and it is suitable for other hazard studies.Many authors describe AHP as a user-friendly, cost-effective and convenient method for flood risk assessment [9,35,[46][47][48][49][50][51][52][53][54][55][56].
This paper investigates the use of the GIS-MCDA methodology for assessing flood hazard in urban areas.The methodology considers six criteria: height, slope, distance to the sewers, distance from the river, the water table and land use.The goal of implementing the GIS interval rough numbers AHP (IR'AHP) methodology is the mapping of flood hazards.This map is the first step in the development plans of flood risk management [57].
Criteria weights in GIS IR'AHP methodology were determined by the IR'AHP technique, which represents a new approach in this field.As one of the most popular methods of multi-criteria decision-making, AHP has been widely used in a variety of policy issues, particularly in the problem of solving criteria weights [58].It provides the ability to measure the consistency of decision-makers (DM) preferences in group decision-making and allows the manipulation of qualitative and quantitative criteria.The final decision when using the AHP method is made on the basis of subjective evaluations of DM [59].Because of the subjectivity and ambiguity that occur in group decision-making, this paper used the interval rough numbers, combined with the AHP method in order to exploit that subjectivity.After determining the weight coefficients applying criteria IR'AHP techniques, a comparison of the results (weights) obtained that way and using fuzzy AHP (F'AHP) and the traditional (crisp) AHP technique was made.The process of comparing results or the comparison of the weight criteria included three scenario modalities.The first scenario involved the use of a new approach, or IR'AHP technique.The second and third scenario implied the application of the F'AHP and traditional (crisp) AHP techniques, respectively.After the validation, using historical events with the generated flood maps of the flood hazard in scenarios, the F'AHP technique that uses interval rough numbers is proven to be the most accurate.
It turned out that the IR'AHP methodology objectively considers uncertainties that occur in group decision-making, which is one of the contributions of this work.Another contribution of this work is the fact that the IR'AHP GIS-MCDA methodology rationally exploits uncertainty in group decision-making and objectively brings out flood hazard maps.The third contribution of this paper is the improvement of the methodology for evaluating the hazard of flooding in urban areas through a new approach to the treatment of uncertainty.The proposed methodology allows the evaluation of alternatives in spite of uncertainties in the process of decision-making and the lack of quantitative information.The results show a significant advantage of the IR'AHP GIS-MCDA methodology in comparison to those who have used the traditional and fuzzy AHP methodology with GIS for the evaluation of the risk of flooding in urban areas [46][47][48][49][50][51][52][53][54][55][56].According to the authors' knowledge, the methodology presented in this paper represents a new approach to the GIS MCDA literature that contributes to the improvement of MCDA techniques.
The aim of this paper is to present a reliable GIS multi-criteria methodology for hazard zones' mapping of flood-prone areas in urban areas.The paper is organized into six sections.In the Introduction section, the importance of the problem of the increasingly frequent occurrence of floods in urban areas is defined and highlighted, thus giving an overview of the GIS and MCDA-AHP methodology relevant to the analysis and preparation of flood hazard maps.The second section contains a brief description of the geographical urban area of Palilula in Belgrade, Serbia.The methodology and the description of the GIS-IR'AHP methodology phases are presented in the third section of the paper.The fourth section covers the estimation of flood-prone areas in Palilula municipality.The fifth section covers the discussion and the results' validation.Concluding considerations are given in the last, sixth section.

Study Area
The study area includes the urban community of Palilula, located in the northern part of the city of Belgrade in Serbia.Palilula extends to 71 km 2 and lies on the alluvial banks of the Danube River, of which 11.9 km 2 is covered by water surfaces (Figure 1).A total of 32% of the municipality is highly urbanized.On the non-urban part of it, mostly on the left river bank, the land is flat, swampy and marshy.The area of Palilula is characterized by the highest groundwater level in Belgrade.According to the last census of 2011, the population of the Palilula urban community is 110,637 persons, and it is constantly increasing due to the present process of urbanization [60].
Water 2017, 9, 360 4 of 25 municipality.The fifth section covers the discussion and the results' validation.Concluding considerations are given in the last, sixth section.

Study Area
The study area includes the urban community of Palilula, located in the northern part of the city of Belgrade in Serbia.Palilula extends to 71 km 2 and lies on the alluvial banks of the Danube River, of which 11.9 km 2 is covered by water surfaces (Figure 1).A total of 32% of the municipality is highly urbanized.On the non-urban part of it, mostly on the left river bank, the land is flat, swampy and marshy.The area of Palilula is characterized by the highest groundwater level in Belgrade.According to the last census of 2011, the population of the Palilula urban community is 110,637 persons, and it is constantly increasing due to the present process of urbanization [60].These changes negatively affect the expansion of impermeable surfaces and cause an increase and acceleration of storm water runoff.The expansion of drains and the sewage network is not adequately monitored with the urbanization and population growth, but the system for collecting rainwater from the area and its transfer mostly relies on the existing sewer system that does not have sufficient capacity and is obsolete.Recipients of rainfall water are existing melioration channels, inadequately maintained and with suspicious functionality due to the urbanization.When the rain is of high intensity for a period of 24 h, the existing systems are unable to absorb rainfall, leading to their spillover phenomena, causing extensive damage of floods.A particular problem of Palilula, especially in areas on the left bank of the Danube, is unplanned construction, a high level of groundwater and unsatisfactory sanitary conditions of life of the population in most parts of the municipality.For these reasons, frequent floods, in addition to considerable material damage, contamination with pollutants, disgorging sewage and septic tanks cause environmental disaster and endanger the health of the population.Remediation of this problem is an urgent problem of the city administration.
The mean annual precipitation for the period of 1851-2015 is approximately 710 mm (Figure 2).The majority of rains fall in June and May.In June, there is 12-13% of the total annual precipitation sum [61].
Data from the Belgrade meteorological station show several average annual precipitations for different periods (Figure 2).It can be seen from the total annual precipitation that for a total of four These changes negatively affect the expansion of impermeable surfaces and cause an increase and acceleration of storm water runoff.The expansion of drains and the sewage network is not adequately monitored with the urbanization and population growth, but the system for collecting rainwater from the area and its transfer mostly relies on the existing sewer system that does not have sufficient capacity and is obsolete.Recipients of rainfall water are existing melioration channels, inadequately maintained and with suspicious functionality due to the urbanization.When the rain is of high intensity for a period of 24 h, the existing systems are unable to absorb rainfall, leading to their spillover phenomena, causing extensive damage of floods.A particular problem of Palilula, especially in areas on the left bank of the Danube, is unplanned construction, a high level of groundwater and unsatisfactory sanitary conditions of life of the population in most parts of the municipality.For these reasons, frequent floods, in addition to considerable material damage, contamination with pollutants, disgorging sewage and septic tanks cause environmental disaster and endanger the health of the population.Remediation of this problem is an urgent problem of the city administration.
The mean annual precipitation for the period of 1851-2015 is approximately 710 mm (Figure 2).The majority of rains fall in June and May.In June, there is 12-13% of the total annual precipitation sum [61].
Data from the Belgrade meteorological station show several average annual precipitations for different periods (Figure 2).It can be seen from the total annual precipitation that for a total of four out of the five analyzed periods, the total annual amount of precipitation was about 700 mm, while in the period of 1901-1930, there was a very strong minimum.However, when it comes to the trend, we cannot talk about some significant changes over long periods of time.
Water 2017, 9, 360 5 of 25 out of the five analyzed periods, the total annual amount of precipitation was about 700 mm, while in the period of 1901-1930, there was a very strong minimum.However, when it comes to the trend, we cannot talk about some significant changes over long periods of time.Lately, there have been more extreme rainfalls, thus influencing the character of fluvial floods in urban areas.The occurrence of fluvial floods in the Palilula area is in a causal relation with the occurrence of extreme rain events.Therefore, predicting such phenomena is of great importance.Till now, several authors have contributed to this matter.Đorđević [62] analyzed trends for the annual and seasonal rainfalls and temperatures along with the extremity indexes of these parameters for Belgrade in the period of 1888-2006.It is found that the annual rainfalls increase while the number of rainy days is lessened for all seasons, except summer.An increase in the frequency of extreme rainfall events is also found, which is of particular importance for the occurrence of floods.Unkašević et al. [63], in a statistical analysis of monthly rainfall and maximum daily precipitation for Belgrade in the period of 1888-1995, indicate the importance of climate changes.According to them, the most common maximum daily precipitation events are associated with more frequent cold fronts followed by torrential rain and stormy weather and cyclones over the Eastern Mediterranean with centers in the southwest coast of the Black Sea.Precisely as a result of the formation of the spatial cyclone, which is shifted over the Adriatic to the Balkans, in the period from 14 to 16 May 2014, the largest flooding in the Palilula area is recorded when in the two days, the measured precipitation was 152.3 mm [62].
Relying on those studies, the increase of the extreme precipitation frequency, and thus, the number of floods with catastrophic consequences, is expected in the future in the area of Belgrade and Palilula.In the prevention of damage caused by floods, it is necessary to know the maximum range of this extreme event.It is of particular importance in the shortest possible time to define the boundaries of flooded areas.The first step is the mapping of flood hazards for Palilula.

3.1.MethodologicalBackground
The methodological hierarchy in this paper is based on the GIS-multi-criteria decision analysis (MCDA) structure.This approach uses the capabilities of GIS in the management of geospatial data and the flexibility MCDA to combine factual information (e.g., land use, slope, drainage system, etc.) with value-based information (e.g., expert opinion, standards, surveys, etc.).From a methodological point of view, the proposed defined flood zones in the municipality of Palilula include the following main steps (Figure 3).Lately, there have been more extreme rainfalls, thus influencing the character of fluvial floods in urban areas.The occurrence of fluvial floods in the Palilula area is in a causal relation with the occurrence of extreme rain events.Therefore, predicting such phenomena is of great importance.Till now, several authors have contributed to this matter.Đorđević [62] analyzed trends for the annual and seasonal rainfalls and temperatures along with the extremity indexes of these parameters for Belgrade in the period of 1888-2006.It is found that the annual rainfalls increase while the number of rainy days is lessened for all seasons, except summer.An increase in the frequency of extreme rainfall events is also found, which is of particular importance for the occurrence of floods.Unkašević et al. [63], in a statistical analysis of monthly rainfall and maximum daily precipitation for Belgrade in the period of 1888-1995, indicate the importance of climate changes.According to them, the most common maximum daily precipitation events are associated with more frequent cold fronts followed by torrential rain and stormy weather and cyclones over the Eastern Mediterranean with centers in the southwest coast of the Black Sea.Precisely as a result of the formation of the spatial cyclone, which is shifted over the Adriatic to the Balkans, in the period from 14 to 16 May 2014, the largest flooding in the Palilula area is recorded when in the two days, the measured precipitation was 152.3 mm [62].
Relying on those studies, the increase of the extreme precipitation frequency, and thus, the number of floods with catastrophic consequences, is expected in the future in the area of Belgrade and Palilula.In the prevention of damage caused by floods, it is necessary to know the maximum range of this extreme event.It is of particular importance in the shortest possible time to define the boundaries of flooded areas.The first step is the mapping of flood hazards for Palilula.

Methodological Background
The methodological hierarchy in this paper is based on the GIS-multi-criteria decision analysis (MCDA) structure.This approach uses the capabilities of GIS in the management of geospatial data and the flexibility MCDA to combine factual information (e.g., land use, slope, drainage system, etc.) with value-based information (e.g., expert opinion, standards, surveys, etc.).From a methodological point of view, the proposed defined flood zones in the municipality of Palilula include the following main steps (Figure 3).

Interval Rough Numbers
Since the expansion of the rough numbers is a new methodology, the next section will illustrate the basic ideas and preferences of interval rough numbers (IRN).The process of group decision-making is accompanied by a large amount of uncertainty and subjectivity, and decision-makers often have a dilemma when allocating the value of certain attributes' decision.These uncertainties and ambiguities are commonly exploited with the application of interval numbers [64,65], fuzzy sets [66][67][68], rough numbers [58,59,69,70], gray theory [71][72][73] and by applying other approaches.The basic idea of applying algorithms to make decisions that are based on interval access (interval numbers, gray theory, and so on) involves the use of interval numbers presenting attribute value decisions.However, the boundaries of the interval in interval numbers are very difficult to determine and are based on the experience and intuition of decision-makers.
A large number of authors uses the process of the multi-criteria decision-making benefits of fuzzy sets to exploit ambiguity as their default [66] or through the different types of fuzzy theory extensions [74][75][76][77][78][79].In addition, fuzzy sets are a very powerful tool for representing imprecision; the choice of membership functions of fuzzy sets is based on subjectivity and shall be based on experience and intuition [80].In the process of decision-making, the intention of interval fuzzy techniques is that the crisp numbers were transformed into fuzzy numbers using the membership functions, which show uncertainties that exist in the real environment.In contrast to the theory of fuzzy sets whose application requires the definition of partial membership functions without clear boundaries set, rough set theory uses the border area set to express concerns.Unlike the fuzzy theory and probability theory in which the degree of uncertainty is defined based on assumptions, the theory of rough sets [81,82] determines the ambiguity on the basis of approximations, which is the basic concept of rough sets.

Interval Rough Numbers
Since the expansion of the rough numbers is a new methodology, the next section will illustrate the basic ideas and preferences of interval rough numbers (IRN).The process of group decision-making is accompanied by a large amount of uncertainty and subjectivity, and decision-makers often have a dilemma when allocating the value of certain attributes' decision.These uncertainties and ambiguities are commonly exploited with the application of interval numbers [64,65], fuzzy sets [66][67][68], rough numbers [58,59,69,70], gray theory [71][72][73] and by applying other approaches.The basic idea of applying algorithms to make decisions that are based on interval access (interval numbers, gray theory, and so on) involves the use of interval numbers presenting attribute value decisions.However, the boundaries of the interval in interval numbers are very difficult to determine and are based on the experience and intuition of decision-makers.
A large number of authors uses the process of the multi-criteria decision-making benefits of fuzzy sets to exploit ambiguity as their default [66] or through the different types of fuzzy theory extensions [74][75][76][77][78][79].In addition, fuzzy sets are a very powerful tool for representing imprecision; the choice of membership functions of fuzzy sets is based on subjectivity and shall be based on experience and intuition [80].In the process of decision-making, the intention of interval fuzzy techniques is that the crisp numbers were transformed into fuzzy numbers using the membership functions, which show uncertainties that exist in the real environment.In contrast to the theory of fuzzy sets whose application requires the definition of partial membership functions without clear boundaries set, rough set theory uses the border area set to express concerns.Unlike the fuzzy theory and probability theory in which the degree of uncertainty is defined based on assumptions, the theory of rough sets [81,82] determines the ambiguity on the basis of approximations, which is the basic concept of rough sets.
In theory, rough sets are used exclusively with in-house knowledge and operational data, and there is no need to rely on model assumptions.In other words, for the application of rough sets, instead of various additional/external parameters, the structure of the data provided is used exclusively [83].In rough sets, the measurement of uncertainty is based on the uncertainty that is already contained in the data [84].Thus, we come to objective indicators contained in the data.
In this paper, a new approach in the theory of rough sets is applied for the measurement of uncertainty that is already contained in the data, an approach that is based on intermittent rough numbers.Since this is a new approach, the authors found in the literature only a few papers that deal with the application of interval rough numbers in the field of multi-criteria decision-making [80].One of the goals of this paper is to initiate other authors for the wider application of IRN in MCDM, since the IRN benefits outlined in this paper represent a logical motive for their wider application.
Suppose there is a set of k classes that represents DM preferences, R = (J 1 , J 2 , ..., J k ), on the condition that they belong to a series, which satisfies the requirement that the J 1 < J 2 <, ..., < J k and second set of m classes that also represents DM preferences, R * = (I 1 , I 2 , ..., I k ).All objects are defined in the universe and connected with the DM preferences.In R * , each object class is presented in the interval I i = {I li , I ui }, where the next condition is satisfied I li ≤ I ui (1 ≤ i ≤ m), and also I li , I ui ∈ R.Then, I li denotes lower interval limit, while I ui denotes upper interval limit of the i-th object class.If both boundaries of a class of objects (upper and lower limit) are listed in a way that is we can define the lower approximations of I * li and I * ui as follows: The upper approximations of I * li and I * ui are defined in the next equations: Both object classes (upper and lower I * li and I * ui ) are defined with their lower limits Lim(I * li ) and Lim(I * ui ) and upper limits Lim(I * li ) and Lim(I * ui ), respectively.
where M L and M * L represent the sum of objects that are contained in the lower approximation of a class of objects I * li and I * ui , respectively.Upper limits Lim(I * li ) and Lim(I * ui ) are defined with Equations ( 7) and (8).
where M U and M * U represent the sum of objects that are contained in the upper approximation of a class of objects I * li and I * ui , respectively.
The rough boundary interval for the lower object class I * li is represented as RB(I * li ) and indicates the interval between the lower and upper limit: while for the upper class of objects, the rough boundary interval I * ui is calculated as: Then, uncertain object classes I * li and I * ui can be shown by their lower and upper limit.
As is seen, each class of objects is defined by its lower and upper limits, which make up an interval rough number, which is defined as: Interval rough numbers characterize specific arithmetic operations, which differ from the arithmetic operations with the classic rough numbers that are shown in Appendix A.

IR'AHP Mathematical Model
The next section describes the mathematical formulation of IR'AHP model.The IR'AHP model generally involves the implementation of the traditional AHP model [85,86].The IR'AHP model allows decision-makers, if there are fluctuations in the comparison in pairs criteria, to make the comparisons matrix, enrolling two values from the Saaty scale on which the decision-maker is ambiguous.Therefore, we get the comparisons matrix in pairs that contain uncertainty that existed for the expert evaluation criteria.In a further implementation of the method, the uncertainty of the IR'AHP model is used to create interval rough numbers, Expressions (1)- (13).The algorithm of the IR'AHP model is implemented through five steps, which are presented in the following: Step 1. Establishing a hierarchical structure of evaluation criteria: A group of experts is formed who carry out the selection criteria and define a problem hierarchy with global goals at the upper and at the lower level criteria.
Step 2. Filling a matrix to compare evaluation criteria in pairs: Members of the group of experts compare in pairs the evaluation criteria in order to define the criteria weights.Comparison in pairs is done using the Saaty ninth degree linguistic scale [86].Each e-th expert presents his/her comparisons by using matrix: where x k ij and x k ij are phrases from Saaty ninth degree linguistic scale used by expert k to present his/her comparison in criteria pairs.
If expert k finds uncertainty while comparing pairs of criteria (i, j), or actually expert k cannot decide between two values from the Saaty ninth degree linguistic scale, then both values are named (x k ij = x k ij ).If there is no uncertainty, then an expert unambiguously selects one value.Then, into the matrix of the comparison criteria (Z k ), the same value (i, j) is entered,x k ij = x k ij .For example, an expert Water 2017, 9, 360 9 of 26 during the criterion comparison at the position (1,2) cannot decide between two linguistic values (e.g., 5 and 6), then at the position (1,2) in the matrix Z k x e 12 = 5 and x e 12 = 6 are entered.Thus, the Z 1 ,Z 2 ,...,Z e matrix is obtained with e experts giving their comparison in criteria pairs.
Step 3. Determination of experts' weight coefficients: For each comparison, matrix Z k is determined by the consistency of experts' evaluation.Saaty [86] suggested the consistency ratio (CR) for the consistency check.The calculation of the degree of consistency is done in two steps.In the first step, the consistency index (CI) is calculated CI = (λ max − n)/(n − 1), where n is matrix rank and λ max the maximum Eigen value of the comparison matrix.
The second step calculates CR as the ratio between CI and the random index (RI).
The random index (RI) depends on matrix rank, and its values are calculated with the random generation of 500 matrices [80]; see Table 1.If the consistency ratio (CR) is less than or equal to 0.10, that proves that the expert was consistent and that there is no need to repeat the evaluations [85].If the CR is higher than 0.10, the decision-maker should repeat (or modify) his/her evaluation in order to improve his/her own consistency.
As every expert in the matrix Z k at (i,j) enters two values, two CRs are obtained for each expert, CRe and CRe'.Prior to determining the final value (CR k ), the values of CRe and CRe' must meet the following requirements CRe ≥ 0.1 and CRe' ≥ 0.1.The final CR is calculated as medium value Experts' weighted coefficients are obtained by normalizing the reciprocal value of the consistency ratio, Equations ( 16) and (17).
where CR k is C r for expert e and W ke the weight coefficient for expert e. Normalization of the experts' weight coefficient is carried out using the additive normalization.
where W i is the weight coefficient for expert e.
Step 4. Construction of the average interval rough comparison matrix: Using Matrix ( 14) from all k experts, two matrices of aggregated experts' sequences are obtained, X* L and X* U . where are sequences used to describe the relative significance of criterion i towards criterion j.Using Equations ( 1)-( 13), each sequence Such rough sequences are defined in the matrices ( 18) and (19).Thus, rough matrices X 1L , X 2L , . . ., X mL are obtained (where m is the number of experts) for the first rough sequence RN x kL ij and X 1 U , X 2'U , . . ., X m U for the second rough sequence RN x k U ij .For the first group of rough matrices X 1L , X 2L , . . ., X mL at (ij) rough sequence is obtained as follows: In the same way, for the second group of rough matrices 21) and ( 22), average rough sequences are calculated: where e is the e-th expert (e = 1, 2, ..., m) and RN(z L ij ) and RN(z U ij ) are rough sequences representing the lower and upper limit of . Therefore, we get the average interval rough comparison matrix in pairs of evaluation criteria (Z): Step 5. Calculation of priority criteria vectors: The priority criteria vector represents an interval rough weighted coefficient IRN(w j ), and it is determined for each n evaluation criteria.Interval rough weighted coefficient IRN(w j ) is calculated using Equations ( 24)- (27).Using Equation (24), matrix Z elements are summed by columns.
The final interval rough weight coefficients (IRN(w j )) for the evaluation criteria are calculated using Equation (27).
where n is the number of evaluation criteria.The values of the criteria weight coefficients are in the interval IRN(w j ) = w L j , w U j , w L j , w U j where the condition 0 ≤ w L j ≤ w L j ≤ w U j ≤ w U j ≤ 1 is fulfilled for each evaluation criteria x j ∈ X.However, it is generally necessary that a sum of criteria weight coefficients is equal to one.In our case, since it is an interval rough weight coefficients criterion, using the expression (27), we obtain the weights Thus, the condition is satisfied for the weight coefficients to be in the interval w i ∈ [0, 1], (i = 1, 2, ..., n).

Estimation of Flood-Prone Areas in Palilula Municipality
This section is divided into subsections.It provides a concise and precise description of the aim of this study.The aim is to propose a reliable GIS-MCDA methodology for the flood hazard mapping of urban areas, which could serve as a useful tool for preventing and reducing the hazard of flood damage for spatial planners to create spatial policies and systems for water management.In the criteria selection and zoning of flood hazard in the local community of Palilula, 10 experts participated.They have experience in the field of crisis management, hydrology, spatial planning and environmental protection.The experts' interviews are used to collect data that are further processed, and the aggregation of experts' opinions is executed.

Criteria Selection
The criteria selection for evaluating the hazard of floods is an important step of the analysis.It is essential to identify floods factors in order to create a reliable susceptibility to flood map [44].Based on previous studies [10,55], experts' opinions and longer observations from the field, this study adopted six criteria that are an important cause of flooding in the local community of Palilula.The selected criteria with a brief description are as follows: Elevation (C1): The study area is located on the alluvial plain of the Danube River and in the lowest part of Belgrade, in the relative range between 69 and 240 m above sea level.This criterion has a key role in controlling the movement of the overflow direction and in the depth of the water level [87].Height values were obtained using a digital elevation model (DEM 25) made at the Military Geographical Institute in Belgrade.
Slope (C2) is an important indicator of a surface zones, which are highly prone to flooding [88,89].Slope is a major factor in determining the rate and duration of water flow.On the flatter surface, water is moving more slowly, collects longer and accumulates so these areas are riskier with respect to the occurrence of floods in relation to the steeper surfaces.Map inclination is prepared as percentages using the DEM of the studied area.
Distance from the drainage network (C3) has the greatest significance for the flood mapping in urban areas.As a result of rainfall overflow, areas near these channels during flooding in urban areas are the most affected [9].In field studies, the sewerage network consists of channel wastewater systems, storm water systems and general system.Data for this criterion were obtained using the map from general solutions for the channeling of rainwater and wastewater from the General Urban Plan of Belgrade, which was established by the Town Planning Institute of Belgrade.
Distance from water surfaces (C4) has a significant impact on the spread and size of floods in the projected area [90].River overflows are one of the main causes for the initiation of flooding.Often, flooding starts from the river or canal riverbed and expands in the region.The areas near bodies of water are very high-hazard areas for the occurrence of floods, and the effect of this criterion is reduced with increasing distance.The study area is dominated by the Danube as the largest body of water to which many open channels and streams gravitate.
Water table (C5) directly affects the infiltration capacity of the soil.In areas with low levels of groundwater, land area with the appearance of water is rapidly saturated, and the water is accumulated and spread to the surrounding area [9].The data used in this study were obtained using the data perennial average groundwater level map from hydrology and geology "groundwater level", prepared by the Directorate for Building Land and Construction of Belgrade.
Land cover use (C6) is also one of the main factors that contributes to the occurrence of floods and has an important impact on runoff and the ability of soil to act as storage for water.Urban and industrial areas are mainly made of impervious surfaces (buildings, roads and parking lots), act as a barrier, reducing the infiltration capacity, retain water and are prone to floods.On the other side, the fields, pastures and forest vegetation are less prone to flooding [22,91].The layer criteria for land cover of the studied area was prepared after integration of the Database Online Map Street within the ArcGIS software package [92].For this study, land use is systematized into eight categories: urbanized areas, industrial areas, agricultural land, land covered with sparse vegetation, grassy areas and parks, forests, marshy areas and water surfaces.
The selected criteria in this paper are based on available data and are adapted to the objective natural and social conditions of the studied zone of Palilula.The applied methodology allows inclusion of other criteria that were not included in this study.This is primarily related to the flow accumulation, flow direction, wetness index, rainfall and the like.All newly-added criteria may have an additional role in the election of the best solution.
Once the criteria are defined, the next step of their evaluation is to build the spatial database and input this into GIS.In this way, each of the criteria has been converted in the form of spatially-defined layers of maps with the same grid cells size 15 m × 15 m, which are the units that are estimated.All of the processes of the transformation and GIS data modeling were performed using the integrated tool ESRI ArcGIS 10.2 software.Maps' criteria height and slope were obtained using the 3D Analyst algorithm based on DEM.Maps from the sewage network and water were obtained by the use of Radial Distance tools.Maps of the groundwater depth are obtained by georeferencing of the analogue map "groundwater level", drawn up by Directorate for building land and construction of Belgrade, while the land cover use map was obtained by importing OpenStreetMap (OSM) bases within the ArcGIS software environment.

GIS-MCDA
This phase involves the standardization, expert work, weighting, summary analysis, aggregation of all criteria to be considered in the decision-making process and their validation.
Given that the data were collected in different ways and have different formats, the first step is that all datasets have to be standardized and in units that can be compared.Based on the literature and the experience of experts, in this study, the fuzzy concept has been used to standardize the data criteria.The fuzzy logic concept is flexible and is suitable for data modeling in which there is no exact boundaries of the elements belonging to the set, determined as zero or one [64].In such cases, the elements belonging to the set are defined on the basis of the degree of affiliation to a function (sigmoidal, J-shaped, linear or user-defined).Which membership functions will be used depends on the nature of the data and on the basis of the decisions and experience of experts.
In this case, with the criteria set whose elements have categorical values (land cover use), we used the discrete classification, in which experts directly determined the values of the elements of fuzzy sets.For other criteria, which are the values of the gradual change from one location to another, elements of the set are standardized using the fuzzy concept based on the linear membership function.A scale ranging from 0 to 1 bytes was used for fuzzification, where zero is the least hazard and one is the most dangerous element of the set value in relation to the likelihood of the occurrence of flooding (Table 2, Figure 4).Given that the data were collected in different ways and have different formats, the first step is that all datasets have to be standardized and in units that can be compared.Based on the literature and the experience of experts, in this study, the fuzzy concept has been used to standardize the data criteria.The fuzzy logic concept is flexible and is suitable for data modeling in which there is no exact boundaries of the elements belonging to the set, determined as zero or one [64].In such cases, the elements belonging to the set are defined on the basis of the degree of affiliation to a function (sigmoidal, J-shaped, linear or user-defined).Which membership functions will be used depends on the nature of the data and on the basis of the decisions and experience of experts.
In this case, with the criteria set whose elements have categorical values (land cover use), we used the discrete classification, in which experts directly determined the values of the elements of fuzzy sets.For other criteria, which are the values of the gradual change from one location to another, elements of the set are standardized using the fuzzy concept based on the linear membership function.A scale ranging from 0 to 1 bytes was used for fuzzification, where zero is the least hazard and one is the most dangerous element of the set value in relation to the likelihood of the occurrence of flooding (Table 2, Figure 4).Water table(C5) Linear monotonically decreasing c = 100 cm; d = 5000 cm 0-100 cm equal to 1, 100-4000 cm between 1 and 0, more than 5000 cm equal to 0 Land cover use (C6) Discrete categorical data Water areas equal 1; wetlands equal 0.9; urbanized areas equal 0.8; industrial areas equal 0.7; agriculture equals 0.5; land covered with sparse vegetation equals 0.4; grass and parks equal 0.2; forests equal 0.1 After standardization, it is necessary for decision-makers to define the significance factors of particular criteria using the appropriate coefficient weights (weights) or the criteria weights.The next section explains a method of determining the criteria weights presented in the studied case.As previously mentioned, the study involved 10 experts who carried out the expert evaluation criteria; see Table 3.After filling in the matrix of comparisons in pairs (Table 3), we see that there is a definite dose of uncertainty when comparing the expert criteria.It is concluded on the basis of the values, which are given out by experts in the comparison matrices.Thus, at the position of C1-C5, in a comparison matrix DM1, the first expert entered two values (z 15 = 4; z 15 = 5).This means that the expert could not choose one of the values four or five.If there is no uncertainty, then the expert e will undoubtedly select one value.An example of this is for the position of C2-C5 in the matrix DM1.The first expert introduced two same value in the comparison matrix DM1 (z 25 = 9; z 25 = 9).
The consistency ratio for the comparison matrix is determined after the criteria comparison in pairs.Since the expert enters two values for each position in matrix Z k , (z e ij ; z e ij ), two consistency ratios are obtained for each expert CR e and CR e .After calculating the consistency ratio for the comparison matrix (Table 4), it can be concluded that the research is valid, because all values CR e < 0.1.Thus, e.g., for the first expert (Table 4), it is CR 1 = (0.022 + 0.088)/2 = 0.055.In order to obtain an interval rough averaged comparison matrix, based on the data in Table 3 and using Expressions (1)-( 13), the elements x e ij and x e ij of matrix X k are transformed into interval rough number IRN(z e ij ).Thus, ten interval rough matrices are calculated X k .Using Equations ( 22), ( 23) and experts' weight coefficients (Table 4), the average interval rough comparison matrix is calculated in evaluation criteria pairs; see Table 5.
Based on data from Table 5, and using Equations ( 25) and ( 26), the normalized matrix of weight coefficients W are calculated (Table 6.).Using (28), interval rough weight coefficients for evaluation criteria are obtained; see Table 7.The process of weight criteria determining in this study was carried out in three scenarios in which, under the AHP arithmetic methods, different modalities' numbers are used (interval rough numbers, triangular fuzzy numbers and crisp value), with the aim of highlighting the most reliable modality.
The paired comparison criterion in fuzzy AHP methods was performed by applying the fuzzy scale with triangular fuzzy numbers; see Table 8 [93][94][95][96].For the comparison in pairs of the traditional AHP method, the classic Saaty scale is used, in Table 8 [86].The most reliable modality has been confirmed by the validation process, which is presented in Section 5. Table 7 shows the interval rough weights and the values of coefficients, which are obtained by the traditional crisp AHP and fuzzy AHP (F'AHP) methods.When calculating the value weights using the F'AHP method, the symmetrical triangular fuzzy numbers form is used.From Table 7, it can be seen that all three methods of generating weight coefficient sequences are of the same rank (C5 > C2 > C7 > C6 > C1 > C4), but with different values.

Aggregation of Weighted Linear Combination
Weighted linear combination (WLC) method is used in the process of criteria map aggregation.In addition, it is compensatory, meaning that low scores in one criterion can be compensated by high scores in another one, which is desired for this particular decision problem.For these reasons, WLC was selected as the method of aggregation.The weighted linear combination (WLC) method multiplies each fuzzy standardized criteria map (i.e., each raster cell 15 m × 15 m) with criteria weights, obtaining different variations from the AHP method, and then sums the results.The mathematical expression for calculating the suitability index in the WLC is given as follows (28): where S is the suitability index, w i is the normalized value of the factor weight and x i is the criterion score of factor i.
On the basis of the adopted criteria and the determination of their weight, according to the three scenarios described in the previous part, the WLC is used to execute the aggregation map of the criteria in the final flood hazard map, which is presented in the same fuzzy value range from 0 to 1.Finally, the flood hazard index is calculated using the defuzzification algorithm by the Standard Deviation method from the Reclass Spatial Analyst Tools ArcGIS 10.2 software release.Based on this, each cell is classified into five categories and receives a new value from 1 to 5 representing the flood hazard index (FHI).
The results of the urban flood hazard assessment are given in the maps of Figure 5.The maps were made according to three different scenarios using the AHP method.The first scenario exploits uncertainty using interval rough numbers in the AHP method (IR'AHP).The second uses the fuzzy technique method for exploiting of the uncertainty in the AHP (F'AHP), and the third scenario considers the use of the traditional (crisp) AHP method.The study area of each map was classified into five sections corresponding to very high, high, moderate, low and very low flood hazard zones.The percentages of these classes, in relation to the entire study area, are given in Table 9. Figure 5 displays the final flood risk maps obtained by the WLC aggregating method for three scenarios.Statistical analysis of the results is obtained from the final flood risk maps (Table 9).According to Scenario 1, when applying the IR'AHP technique, very high surface areas (FHI 5) in Palilulaare 12.9 km 2 , which is about 18.5% of the territory of the community.Under Scenario 2, which uses the fuzzy technique in the AHP method (F'AHP), the FHI 5 area is 8.8 km 2 , while in Scenario 3, with the traditional (crisp) AHP method applied, the FHI 5 area is calculated to be 7.8 km 2 for the occurrence of floods.Spatially, all three maps shows that these are the parts on the left bank of the Danube, which gravitate towards the sewage network.At the same time, these are the most ecologically-sensitive parts in which uncontrolled floods can lead to catastrophic environmental and social consequences.

Results-Discussions
Validation of flood hazard maps was performed on the basis of historical flood events that have been recorded in the area of research.Data on many high waters and floods in the past have been gathered from all relevant stakeholders involved in flood protection: the Republic Hydro-meteorological Service, the Republic Water Directorate, public water companies and the competent authorities of Palilula.Some data, for flooding that occurred in 2010 and 2014, by the "Questionnaire on flood event", were collected directly in the field.In the period from 1938 to 2016, 31 major flood events were identified caused by spillover of the sewage network and the spill from  Statistical analysis of the results is obtained from the final flood risk maps (Table 9).According to Scenario 1, when applying the IR'AHP technique, very high surface areas (FHI 5) in Palilulaare 12.9 km 2 , which is about 18.5% of the territory of the community.Under Scenario 2, which uses the fuzzy technique in the AHP method (F'AHP), the FHI 5 area is 8.8 km 2 , while in Scenario 3, with the traditional (crisp) AHP method applied, the FHI 5 area is calculated to be 7.8 km 2 for the occurrence of floods.Spatially, all three maps shows that these are the parts on the left bank of the Danube, which gravitate towards the sewage network.At the same time, these are the most ecologically-sensitive parts in which uncontrolled floods can lead to catastrophic environmental and social consequences.

Results-Discussions
Validation of flood hazard maps was performed on the basis of historical flood events that have been recorded in the area of research.Data on many high waters and floods in the past have been gathered from all relevant stakeholders involved in flood protection: the Republic Hydro-meteorological Service, the Republic Water Directorate, public water companies and the competent authorities of Palilula.Some data, for flooding that occurred in 2010 and 2014, by the "Questionnaire on flood event", were collected directly in the field.In the period from 1938 to 2016, 31 major flood events were identified caused by spillover of the sewage network and the spill from the Water 2017, 9, 360 18 of 26 bed of the channel and the river, mainly on the sections along which there are systems built for flood protection, but also on protected areas due to spill-over protective structures.
In Figure 6, the green circles show the locations of historical flood events.The majority of the city's historic flooding is located left, next to the Danube River, and in the urban part on the left bank of the Danube.
For the validation of the final results obtained by applying the GIS MCDA methodology for flood hazard assessment, it is necessary to determine the spatial relationship between the historical flooding sites and the flood hazard maps made.For the purposes of this analysis, we used the algorithm "Extract by Mask" integrated within the software environment ArcGis10.2.Based on this tool, the cell grid that corresponds to historically-flooded areas, based on the spatial coincidence, is extracted in one of five different levels of flood hazard, according to different scenarios.The results of this analysis are presented in Table 10.
Water 2017, 9, 360 18 of 25 the bed of the channel and the river, mainly on the sections along which there are systems built for flood protection, but also on protected areas due to spill-over protective structures.In Figure 6, the green circles show the locations of historical flood events.The majority of the city's historic flooding is located left, next to the Danube River, and in the urban part on the left bank of the Danube.
For the validation of the final results obtained by applying the GIS MCDA methodology for flood hazard assessment, it is necessary to determine the spatial relationship between the historical flooding sites and the flood hazard maps made.For the purposes of this analysis, we used the algorithm "Extract by Mask" integrated within the software environment ArcGis10.2.Based on this tool, the cell grid that corresponds to historically-flooded areas, based on the spatial coincidence, is extracted in one of five different levels of flood hazard, according to different scenarios.The results of this analysis are presented in Table 10.Based on the validation results from Table 10, a relatively high consistency of the applied GIS-AHP methodology can be seen.The biggest consistency of historical flooding events of the final flood hazard map was drawn up on the basis of the IR'AHP method use, according to Scenario 1. Scenario 2 has been less consistent with the final flood risk map made using F'AHP.The least consistency has been shown using Scenario 3, in which the criteria weights were calculated using the traditional (crisp) AHP method.
Considering the flood hazard map in Scenario 1 (Figure 6), it can be seen that 27% or 87.1% of the recorded historical flood locations coincide with the zone of very high flooding hazard, while the four floods coincide with the zone of very high flooding hazard.In this way, the reliability of the proposed methodology is confirmed, thus ensuring the certainty of the results of this analysis.Therefore, it can be concluded that the IRN's use is justified in the process of determining the criteria weights, since the proposed IR'AHP methodology successfully exploited ambiguities and objectively reflects the perceptions of decision-makers.Based on the validation results from Table 10, a relatively high consistency of the applied GIS-AHP methodology can be seen.The biggest consistency of historical flooding events of the final flood hazard map was drawn up on the basis of the IR'AHP method use, according to Scenario 1. Scenario 2 has been less consistent with the final flood risk map made using F'AHP.The least consistency has been shown using Scenario 3, in which the criteria weights were calculated using the traditional (crisp) AHP method.
Considering the flood hazard map in Scenario 1 (Figure 6), it can be seen that 27% or 87.1% of the recorded historical flood locations coincide with the zone of very high flooding hazard, while the four floods coincide with the zone of very high flooding hazard.In this way, the reliability of the proposed methodology is confirmed, thus ensuring the certainty of the results of this analysis.Therefore, it can be concluded that the IRN's use is justified in the process of determining the criteria weights, since the proposed IR'AHP methodology successfully exploited ambiguities and objectively reflects the perceptions of decision-makers.
Spatial analysis obtained by overlapping flood hazard zone maps with different thematic content maps can give general conclusions about the potential risk of flooding in Palilula.Table 11 shows the data of the exposure of certain land use and infrastructure elements according to the flooding hazard zones.From Table 11, it can be seen that about 75 km of roads and 3 km of railways and 950 residential buildings with about 72,500 inhabitants are located in zones of high and very high flood hazard.Furthermore, this zone covers 15 km 2 of urban area.The general conclusion is that the risk of flood covers much of the demographic and infrastructural resources of Palilula municipality.

Conclusions
The paper tests the application of the GIS-IR'AHP modified methodology for zoning flood hazard in urban areas in the case of the Municipality of Palilula in Belgrade, Serbia.It is a methodology that can be applied to other sites with similar characteristics.The association of these tools is provided in the Advanced ArcGIS 10.2 software environment of the company ESRI.The methodology considers six factors that are relevant to the hazard of flooding in urban areas: the height, slope, distance to the sewage network, the distance from the water surface, the depth of groundwater and land use.The process of comparing and determining the weight criteria included three scenarios used under the AHP method: interval rough numbers, fuzzy numbers and crisp value.Final maps of areas of the hazard of flooding, according to different scenarios, were obtained using weighed linear combination (WLC).The validation of the results was carried out on the basis of comparison, an event of historic floods that were recorded in the area of research, with defined zones of hazard on the final maps.
Results of the validation show that Scenario 1 is the most consistent with the historical events of flooding.Based on these indicators, the authors propose an original approach for determining the weight criteria implementing the new methodology IR'AHP.The results of the GIS-IR'AHP methodology of Scenario 1 resulted in the final aggregation of risk maps of flood zones, the definition of the flood hazard index (FHI).They show that 12.9 km 2 or 18.5% of the Palilula municipality has a very high-hazard, the high-hazard area being 23.2% and 21.6% for moderate hazard with respect to the occurrence of flood events.On the other hand, 8.5 km 2 or 12.2% of the area has very low and 24.6% has low hazard in the occurrence of floods.From these results, it is evident that Palilula municipality is mostly located in areas of high hazard with the probability of the occurrence of high water and flooding.These parts are mostly on the left bank of the Danube and parts that gravitate toward the sewage canal network of the community.
Application of GIS-IR'AHP in defining flood hazard zones in urban areas proved to be justified because it is based on the adopted criteria, with a high degree of reliability of the differentiated parts on the premises of a high-hazard of the occurrence of floods.The result is a mapping of flood hazard, which is the first step in the development of flood hazard management plans and determining areas where there is or might appear significant hazard of flooding, with harmful consequences for human health, the environment, economic activities and cultural heritage.In addition to the guidelines in water management, the application of the methodology can serve spatial planners to create a space policy in the area of urban land use plans and insurance companies for the evaluation of potential flood damage.The results of the work, in front of spatial planners and city government, impose the necessity to take austerity measures regarding uncontrolled urbanization, particularly in high-hazard areas near bodies of water and clogging city canal water.In addition, it is necessary to define new urban plan construction rules that require the increase of the permeable surface and the increase of the value of the runoff coefficient to surface runoff.Furthermore, one of the potentially good measures related to the collection and storm water drainage in the municipality of Palilula is to partially convert sewage into the general sewer separation system.However, only in places where it is close to the recipient and where it is estimated that the water will be too polluted, it does not have to go to the treatment plant prior to discharging to the recipient.
This study presents evidence of the importance and reliability of the application of the mixed GIS and MCDA technique in the assessment of phenomena that require expert opinion and manipulate a large number of data from different sources.The current issue discusses the new GIS-IR'AHP methodology, which creates the basis for further theoretical and practical upgrade.The present methodology allows the inclusion of new criteria that were not considered in the work, which are relevant to this issue (flow accumulation, flow direction, wetness index, rainfall, soil type, etc.).Furthermore, it is possible to use digital urban projects and large-scale orthophotos and satellite images with high spatial resolution, to provide detailed mapping of hazard zones in urban areas, which can improve the results.All of these factors may have an additional role in the upgrading of the present methodology and choosing the best solutions.

Figure 1 .
Figure 1.Location of the study area.

Figure 1 .
Figure 1.Location of the study area.

Figure 2 .
Figure 2. Average annual value of rainfall for the Palilula area.

Figure 2 .
Figure 2. Average annual value of rainfall for the Palilula area.

Figure 3 .
Figure 3. Flowchart of the applied methodology.

Figure 3 .
Figure 3. Flowchart of the applied methodology.

Figure 4 .
Figure 4. Maps of standardized criteria for the studied area.Figure 4. Maps of standardized criteria for the studied area.

Figure 4 .
Figure 4. Maps of standardized criteria for the studied area.Figure 4. Maps of standardized criteria for the studied area.

Figure 5 .
Figure 5. Final flood hazard maps for the scenarios.

Figure 5 .
Figure 5. Final flood hazard maps for the scenarios.

Figure 5
Figure5displays the final flood risk maps obtained by the WLC aggregating method for three scenarios.Statistical analysis of the results is obtained from the final flood risk maps (Table9).According to Scenario 1, when applying the IR'AHP technique, very high surface areas (FHI 5) in Palilulaare 12.9 km 2 , which is about 18.5% of the territory of the community.Under Scenario 2, which uses the fuzzy technique in the AHP method (F'AHP), the FHI 5 area is 8.8 km 2 , while in Scenario 3, with the traditional (crisp) AHP method applied, the FHI 5 area is calculated to be 7.8 km 2 for the occurrence of floods.Spatially, all three maps shows that these are the parts on the left bank of the Danube, which gravitate towards the sewage network.At the same time, these are the most ecologically-sensitive parts in which uncontrolled floods can lead to catastrophic environmental and social consequences.

Figure 6 .
Figure 6.Validation of the final flood hazard map (IR'AHP) with historically-flooded points.

Figure 6 .
Figure 6.Validation of the final flood hazard map (IR'AHP) with historically-flooded points.
respectively, then, we can define two new sets containing the lower class of objects R * l = (I * l1 , I * l2 , ..., I * lj ) and the upper class of objects R * u = (I * u1 , I * u2 , ..., I * uk ), respectively.Then, for any class of objects

Table 1 .
Random index (RI) value depending on the matrix rank.
where Lim(x kL ij ) and Lim(x k U ij ) are upper limits, while Lim(x kL ij ) and Lim(x k U ij ) are upper limits of rough sequences RN x kL

Table 2 .
Fuzzy standardization of criteria.

Table 3 .
Matrices of experts' comparison in theevaluation criteria pairs.

Table 4 .
CR e comparison matrix and experts' weights.

Table 5 .
Interval rough average matrix.

Table 6 .
Normalized matrix of weight coefficients.

Table 7 .
Weight coefficients for evaluation criteria.

Table 9 .
The areas of classes from the final flood hazard map for the scenarios.IR, interval rough number; F, fuzzy.

Table 9 .
The areas of classes from the final flood hazard map for the scenarios.IR, interval rough number; F, fuzzy.

Table 10 .
Spatial relations of the historical floods locations and flood hazard zones according to the scenarios.

Table 10 .
Spatial relations of the historical floods locations and flood hazard zones according to the scenarios.

Table 11 .
The exposure to the flooding hazard in Palilula.