Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process for Flood-Hazard Assessment in the Megalo Rema River Catchment (East Attica, Greece)

: This study deals with the ﬂood-hazard assessment and mapping in the catchment of Megalo Rema (East Attica, Greece). Flood-hazard zones were identiﬁed utilizing Multi-Criteria Decision Analysis (MCDA) integrated with Geographic Information System (GIS). Five factors were considered as the most inﬂuential parameters for the water course when high storm-water runoff exceeds drainage system capacity and were taken into account. These factors include slope, elevation, distance from stream channels, geological formations in terms of their hydro-lithological behavior and land cover. To obtain the ﬁnal weights for each factor, rules of the Analytic Hierarchy Process (AHP) were applied. The ﬁnal ﬂood-hazard assessment and mapping of the study area were produced through Weighted Linear Combination (WLC) procedures. The ﬁnal map showed that approximately 26.3 km 2 , which corresponds to 22.7% of the total area of the catchment, belongs to the high ﬂood risk zone, while approximately 25 km 2 , corresponding to ~15% of the catchment, is of very high ﬂood risk. The highly and very highly prone to ﬂooding areas are located mostly at the southern and western parts of the catchment. Furthermore, the areas on both sides of the channel along the lower reaches of the main stream are of high and very high risk. The highly and very highly prone to ﬂooding areas are relatively low-lying, gently sloping and extensively urbanized, and host the densely populated settlements of Raﬁna-Pikermi, Penteli, Pallini, Peania, Spata, Glika Nera, Gerakas and Anthousa. The accuracy of the ﬂood-hazard map was veriﬁed by correlating ﬂood events of the last 30 years, the Hydrologic Engineering Center’s River Analysis System (HEC–RAS) simulation and quantitative geomorphological analysis with the ﬂood-hazard level. The results of our approach provide decision makers with important information for land-use planning at a regional scale, determining safe and unsafe areas for urban development.


Introduction
Floods are among the most frequent and dangerous natural hazards and the leading cause of natural disaster fatalities worldwide [1]. They represent approximately one-third of all global hazardous events, and the number of extreme flood incidences has significantly increased over the past few decades [2]. Every year this type of hazard causes a significant loss of life and property [3] and severely affects the natural and human environment as well as the development of an area [4].
Floods in urban and peri-urban areas are caused by the complex combination of both natural and human-induced factors [5,6]. Among the most important physical factors are extreme precipitation events and the "organization" of the drainage network, as well as the geological and geomorphological conditions of the catchment [7]. Moreover, human

Study Area
The Megalo Rema River, located in East Attica, Greece, is an ephemeral stream and drains an area of approximately 115 km 2 ( Figure 1). The drainage basin reaches a maximum elevation of 950 m at its northeastern border (Penteli Mt.). The catchment includes the residential areas of the Rafina, Pikermi, Ntaou Penteli, Pallini, Gerakas, Anthousa and Glika Nera regions with a total population coverage of~117,000 residents. The sixth order (according to Strahler's [41] ordering system) main stream channel has a length of 25 km and follows a WSW-ENE flow direction. It discharges into the South Evoikos Gulf where the town of Rafina is located, with a population of~13,000 inhabitants.
In the asymmetric Megalo Rema basin the drainage system is well developed north of the main channel, whereas to the south a few channels of low order exist ( Figure 1) [42]. The catchment can be divided geomorphologically into two areas: (a) the relatively mountainous northern part of rough relief consisting mainly of metamorphic rocks (schists and marbles) from the Mesozoic age belonging to the autochthonous Almyropotamos geotectonic unit, and (b) the southern area of lower elevations, gentle slopes and a generally smoother relief that is composed mainly of lacustrine marls, silts, marly limestones and conglomerates of the Upper Miocene age [43,44].
The climate of the study area is typically subtropical Mediterranean, with prolonged hot and dry summers succeeded by considerably mild and wet winters. The mean annual precipitation ranges from 362 mm at the mouth of the river to 473 mm in the SW part of the catchment at Mount Ymittos, while the mean annual air temperature is 17.2 • C. The rainy period begins in October and ends in March.
The Megalo Rema drainage network is part of the GR06RAK0003 potentially significant flood risk (PSFR) zone, which was defined by the Greek Special Secretariat for Water (SSW) [45] under the European Council (EC) Floods Directive 2007/60/EC [46]. The Megalo Rema catchment is an area greatly prone to flash flooding. It is estimated that over the last twenty-seven years, one flash-flood event occurred every 2.7 years [47]. In addition, more than thirty-two flood events have affected the area between 2004 and 2014, with eighteen of these characterized as significant flood episodes based on the hydrological response intensity of the study catchment. The application of a system based on the coupling between the Weather Research and Forecasting (WRF) numerical weather prediction model and its hydrological extension package (WRF-Hydro) by Giannaros et al. [48] showed that the majority of the flash flood events in Megalo Rema catchment took place during the Sustainability 2021, 13, 10232 4 of 24 wet period of the hydrological year and were associated with typical, for the study area, wet-season cyclonic activity [49]. This is due to the global atmospheric circulation during the autumn and winter, which interacts with the complex geomorphology and land-sea temperature contrast in the Eastern Mediterranean region, favoring the development of cyclonic atmospheric conditions [50]. In particular, the flood episodes of the study area were driven by low-pressure systems, which, in most cases, affected the catchment while moving from the west towards the east [49]. However, the recent extreme flood events in Europe, caused by an abnormality in cyclonic seasonal activity, show an increase in rainfall intensity and a rise in frequency of flood rates due to climate change [51]. In Greece, even though the precipitation levels are predicted to decline an average of 17% annually, the flood hazard is expected to rise, thus increasing the cost of direct damage from floods up to 10% annually [52]. This means that flood rates and the period of occurrence of past events may not be reliable in the future due to climate change. Moreover, large wildfire-burn territories, increased urbanization, mild topographic slopes and the absence of efficient inundation protection make the catchment greatly prone to flash-flooding [42,49].
response intensity of the study catchment. The application of a system based on the coupling between the Weather Research and Forecasting (WRF) numerical weather prediction model and its hydrological extension package (WRF-Hydro) by Giannaros et al. [48] showed that the majority of the flash flood events in Megalo Rema catchment took place during the wet period of the hydrological year and were associated with typical, for the study area, wet-season cyclonic activity [49]. This is due to the global atmospheric circulation during the autumn and winter, which interacts with the complex geomorphology and land-sea temperature contrast in the Eastern Mediterranean region, favoring the development of cyclonic atmospheric conditions [50]. In particular, the flood episodes of the study area were driven by low-pressure systems, which, in most cases, affected the catchment while moving from the west towards the east [49]. However, the recent extreme flood events in Europe, caused by an abnormality in cyclonic seasonal activity, show an increase in rainfall intensity and a rise in frequency of flood rates due to climate change [51]. In Greece, even though the precipitation levels are predicted to decline an average of 17% annually, the flood hazard is expected to rise, thus increasing the cost of direct damage from floods up to 10% annually [52]. This means that flood rates and the period of occurrence of past events may not be reliable in the future due to climate change. Moreover, large wildfire-burn territories, increased urbanization, mild topographic slopes and the absence of efficient inundation protection make the catchment greatly prone to flashflooding [42,49].

Flood-Hazard Criteria
For the flood-hazard assessment, a GIS-based MCDA was implemented ( Figure 2). The analysis was based on five factors relevant to flood-hazard as documented in the literature [53]. The selected factors include slope, elevation, distance from the channels of the higher order streams, hydro-lithological characteristics of the geological formations and land cover. These factors control the water route when drainage system capacity is

Flood-Hazard Criteria
For the flood-hazard assessment, a GIS-based MCDA was implemented ( Figure 2). The analysis was based on five factors relevant to flood-hazard as documented in the literature [53]. The selected factors include slope, elevation, distance from the channels of the higher order streams, hydro-lithological characteristics of the geological formations and land cover. These factors control the water route when drainage system capacity is exceeded by high runoff and have been proved effective when included in flood-hazard assessment studies and applications [8,15,17,38]. weathering blocks their fracture system [63]. Higher rating 3 has been assigned to impermeable geological formations due to their lower infiltration capacity.

Figure 2.
Flow diagram of the methodology for the present study.

Analytical Hierarchy Process (AHP)
The final weights of each factor were defined following the AHP [39,40]. The first step in the AHP was the computation of the pair-wise comparison matrix, where each entry represents the relative significance of each factor/criterion to the others. The method of comparisons per pair is the most commonly used process for the calculation of criteria weight coefficients. AHP is based on the allocation of weights to the criteria used, according to their importance. The relative importance between two factors was measured according to a numerical scale from 1 to 9. The correlation between the numerical values and the intensity of importance was as follows: 1 = equal, 2 = weak or slight, 3 = moderate, 4 = moderate plus, 5 = strong, 6 = strong plus, 7 = very strong, 8 = extremely strong and 9 = of extreme importance. Inversely, less important variables were rated between 1 and 1/9 [19,20]. Completion of the degree of significance between two factors requires field experience, knowledge of the subject and/or the opinion of the experts, who assign the value aij according to their judgment for the relative importance of one criterion over the other [65]. If a criterion has equal importance with the one it is being compared with, then the preference value (aij) is equal to one. If the relative importance of a criterion is higher than the criterion it is being compared with, then aij > 1 and the reciprocal property is aji = 1/aij. When comparing a criterion with itself, the assumption is that aii = 1 for all the n criteria. In this way, a preference table is created (Table 1). Table 1. Pairwise comparison matrix of the criteria (preference value: aij = wi/wj, where wi/wj is the relative importance of factor i to factor j).
After assigning the preference values (αij), the columns of Table 1 are summed to Sum_j and then the preference values of each column (αij) in Table 1 are divided by the sum of the corresponding column (Sum_j) and the values are summed again, in rows this A crucial step in the flood-hazard assessment methodology was the classification of all of the factors. Hence, each factor was divided into classes with specific boundary values. The classes of the involved factors have to be standardized to a uniform suitability rating scale. The standardization method used in the analysis was consistently based on a fiveand three-grade scale, according to the hazard level each criterion contributes to the total flood-hazard. Integer numbers, ranging from 0 to 4 and 1 to 3, were assigned to every class of the five-or three-grade scaled criteria, respectively. The class which was rated as 0 represented no hazard level whereas the one rated as 4 represented a very high hazard level. Thus, each of the five factors involved in the flood-hazard assessment is represented by a spatial distribution map of the classified values, reclassified to a regular scale, and then by a thematic map of the spatial distribution of the standardized values.

Slope
Surface runoff and the water accumulation process in any geomorphic setting relies upon its surface slope appropriation [54]. Water flows from higher to lower elevations, which therefore means that slope influences the amount of surface runoff and infiltration. Flat and gently sloping areas (0-18 • ) in low elevations are more prone to flooding [55]. In these areas, the surface water runoff is usually accumulated easier and, as a result, intense precipitation can flood these low-slope areas faster than areas in higher elevation with a steeper slope.
The slope of the study area was computed utilizing ArcGIS/ArcMap ESRI ® by inputting elevation data derived from the detailed 5 m resolution DEM of the catchment, obtained from the Hellenic Cadastre (Ktimatologio S.A.). The resulting slope map consists of a 5 m resolution raster layer which was then converted to a vector layer, which represents the steepness of the terrain of the study area in degrees. The slope values were classified into five classes, based on the correlation of the slope with previous flood events. The classes were assigned a value from 0 to 4. Therefore, the gently sloping parts (≤2 • , 2-6 • , and 6-12 • ) of the catchment are considered as prone to flooding areas and were assigned the rating of 4, 3 and 2, respectively (see Table 4). Steeper parts of the area with a slope between 12 and 20 • and >20 • were assigned the ratings of 1 and 0.

Distance from Stream Channels
River-overflows are crucial for the initiation of a flood event. According to Predick and Turner [56], proximity to the stream channel critically increases the possibility of experiencing a riverine flood event. During extreme rainfalls, high discharges often cause a Sustainability 2021, 13, 10232 6 of 24 rapid rise in river level, leading to an overflow of water from the riverbed and inundating the surroundings. Hence, regions situated near rivers endure frequent flooding [57]. In contrast, areas far from stream channels generally have a lower probability of flooding, since the likelihood of being within the flood-induced riverbed is lower.
Areas near the stream channels are considered to have a high flood-hazard, whereas the effect of this parameter decreases with distance from the channels. It is obvious that there is no universal agreement on the critical distance that can have higher susceptibility, and this distance changes from river to river. In this study, the distance from the stream channel was taken into account in an analogical manner, by comparing the distances which were used in similar studies in correlation with the catchment size (see [8,15]).
The drainage network was acquired from the 1:5000 scale topographic maps of the Hellenic Military Geographical Service. This factor was calculated by imposing buffer zones within the ArcMap interface, by creating polygons enclosing the area on either side of the higher order stream channels, for each given distance from the channel. Records of historical floods in the study area [49] as well as the results of the simulation of the river hydraulics' behavior-using the Hydrologic Engineering Center's River Analysis System (HEC-RAS) developed by U.S. Army Corps of Engineers Hydrologic Engineering Center in GIS Environment-for the last 4 km before the river mouth [47], were particularly useful in defining the classes of this criterion. For the third order stream channels of the drainage network, 10, 20, 30, 40 and >40 m buffer zones were generated. Around the fourth order stream channels, buffer zones at the distances of 15, 30, 45, 50 and >50 m were created. For the fifth order streams, buffer zones at distances of 20, 40, 60, 80 and >80 m from the channels were considered, whereas buffer zones at the distances of 25, 50, 70, 100 and >100 m from the sixth order main channel of Megalo Rema were imposed. The shorter the distance, the higher the hazard level and the assigned value, and the longer the distance, the lower the hazard level and the assigned value.

Land Cover
Several hydrological processes such as surface runoff, infiltration rate and evapotranspiration, as well as interrelationship between surface and groundwater, are being significantly controlled by the land cover pattern of an area. As a result, land cover is considered an important parameter in flood-hazard assessment [58].
An integrated land cover layer of the drainage basin was created by the information obtained from CORINE Land Cover (CLC 2018) of Copernicus Land Monitoring Service and Agricultural Blocks (ILOTS 2012) of the Ministry of Rural Development and Food. The land cover type was classified into five categories: (i) densely urbanized areas, (ii) sparsely urbanized areas, (iii) agricultural land, (iv) land covered by sparse vegetation and (v) forests. Densely and sparsely urbanized areas that support the overland flow of water were assigned the highest ratings of 4 and 3, respectively, whereas forests have been assigned the lower rating of 0 since they favor infiltration.

Elevation
According to several relevant research studies and applications, elevation is one of the dominating factors that control flood occurrence [59][60][61]. Because of the gravitational force, water moves rapidly from upland to lowland areas and the water gets spread over the lower elevated plains, causing inundation [62]. However, it is not possible for a widely agreed flood elevation threshold to be set due to the various morphological/morphometric characteristics of each catchment. The relationship between elevation and flood events in the Megalo Rema catchment shows that almost 82% of the past events occurred in elevations lower than 200 m, while areas of elevations higher than 500 m have not been affected by floods.
The elevation grid of the study area was produced by the 5 m resolution DEM, and elevation values were divided into five categories by evaluating the elevation of the historical flood events in the study area catchment. The low-lying areas, with elevation ≤50 m, were assigned the highest rating 4, since they were considered as more prone to flooding.

Geological Formations
The hydro-lithological characteristics of the geological formations of an area are an important criterion, because they may amplify or extenuate the magnitude of flood events. The presence of permeable formations favors water infiltration while impermeable rocks favor surface runoff.
The geological formations of the catchment were derived from the corresponding 1:50,000 geological map of the Institute of Geology and Mineral Exploration [43,44], and were classified according to their hydro-lithological behavior into three categories: permeable, semi-permeable and impermeable. Permeable geological formations include tectonically fractured and karstified limestone and marbles, sandstones and alluvial deposits [63,64]. Marls, conglomerates and silts are considered semi-permeable formations while schist is considered impermeable rock, since the fine-grained material produced by their weathering blocks their fracture system [63]. Higher rating 3 has been assigned to impermeable geological formations due to their lower infiltration capacity.

Analytical Hierarchy Process (AHP)
The final weights of each factor were defined following the AHP [39,40]. The first step in the AHP was the computation of the pair-wise comparison matrix, where each entry represents the relative significance of each factor/criterion to the others. The method of comparisons per pair is the most commonly used process for the calculation of criteria weight coefficients. AHP is based on the allocation of weights to the criteria used, according to their importance. The relative importance between two factors was measured according to a numerical scale from 1 to 9. The correlation between the numerical values and the intensity of importance was as follows: 1 = equal, 2 = weak or slight, 3 = moderate, 4 = moderate plus, 5 = strong, 6 = strong plus, 7 = very strong, 8 = extremely strong and 9 = of extreme importance. Inversely, less important variables were rated between 1 and 1/9 [19,20]. Completion of the degree of significance between two factors requires field experience, knowledge of the subject and/or the opinion of the experts, who assign the value a ij according to their judgment for the relative importance of one criterion over the other [65]. If a criterion has equal importance with the one it is being compared with, then the preference value (a ij ) is equal to one. If the relative importance of a criterion is higher than the criterion it is being compared with, then a ij > 1 and the reciprocal property is a ji = 1/a ij . When comparing a criterion with itself, the assumption is that a ii = 1 for all the n criteria. In this way, a preference table is created (Table 1). Table 1. Pairwise comparison matrix of the criteria (preference value: a ij = w i /w j , where w i /w j is the relative importance of factor i to factor j). After assigning the preference values (α ij ), the columns of Table 1 are summed to Sum_j and then the preference values of each column (α ij ) in Table 1 are divided by the sum of the corresponding column (Sum_j) and the values are summed again, in rows this time (Row_sumi) ( Table 2). Finally, the set of each line (Row_sumi) is divided by the number of variables (n) and the result equals the weighting coefficient (W i ) of each criterion. AHP requires normalization of all factor weights, which was achieved using the following equation: After the calculation of the weight values, it is important to verify the consistency of each table matrix. The implication of each one was checked with the Consistency Ratio (CR). This ratio is used to avoid the creation of any incidental judgment in the matrix. When CR < 0.1, an acceptable level of consistency has been achieved, while judgments are tolerated if 0.1 < CR < 0.2 and rejected if CR > 0.2. CR is given by the following equation: where RI is the Random Index: a constant which depends on the order of the matrix (see Saaty, 1987; p. 171) [67], and CI is the Consistency Index calculated by the equation: where n is the number of items compared in the matrix and λmax is the maximum value of the eigenvalue that is obtained by the equation: where n ∑ i = 1 (WV i /W i ) = vector coherence (C), W i = the weighting coefficient estimated according to Table 2 and WV i = weighted sum vector (calculated according to Table 3).

Weighted Linear Combination (WLC)
The total score of the basic flood-hazard assessment for the study area was calculated with the linear combination of the selected factors/parameters, taking into account the relative weights. This involves superimposing the thematic maps with different weights in a vector-based GIS environment. The inclusion of the estimated factors was performed using the Weighted Linear Combination method, according to the following mathematical formula: where H is hazard degree, n is the number of factors, W i is the weight of factor i and X i is the rating of factor i, as shown in Table 4. The flood-hazard level scores of the study area were then classified into five classes using the quantile classification method, and the final basic flood-hazard map was created. Class 1 of the lower values corresponds to areas of very low flood susceptibility, class 2 corresponds to parts of the catchment with low flood susceptibility, class 3 to moderate susceptibility, class 4 to high susceptibility and, finally, class 5 of the highest values corresponds to extremely flood-prone areas.
The influence of the uncertainty of the adopted factor weights on the flood-hazard assessment was estimated using the following formula [8,68]: where ∆S is the error produced by independent errors ∆W i in the weighting coefficient values and X i is the rating of factor i. Equation (6) was applied to calculate the error (∆S). The independent errors (∆W i ) were set as 20% of each original factor weight for all the factors at the same time [8,69]. After, it was multiplied by 1.96 in order to compute 95% confidence level in the suitability values S. The map created by this process was used to calculate the upper and lower S values at 95% confidence level, once by adding it and once by subtracting it from the basic flood-hazard map. The two resulting maps represent the scenarios of maximum and minimum S values for the catchment.

Verification of the Flood-Hazard Assessment
For verification of the flood-hazard assessment, the spatial distribution of the past 45 severe flood events that occurred in the Megalo Rema River catchment was used, as shown in Figures 7a,b and 8. The geographic distribution of the 45 flood events was produced by plotting, as a map layer, the points referring to the sites affected by severe flood incidences. These flood events, which occurred during the last 30 years, were retrieved from the archives of the Region of Attica, the Hellenic Fire Service, the Hellenic Agricultural Insurance Organization and the Ministry of Environment and Energy. The verification was performed by applying a frequency ratio statistical analysis. For this purpose, firstly, the frequency distribution of flood events was calculated for each hazard level zone of the basic flood-hazard map (see Table 9). Then, the ratio of each flood-hazard level zone area (from very low to very high) to the total area of the catchment was computed. The frequency ratio for each flood-hazard level zone was calculated by dividing the percentage of the events (out of the total 45 events) which appear in each hazard level zone by the area percentage of each hazard level zone area, which is presented in Table 7.
The accuracy of the basic flood-hazard map was also examined, taking into consideration the results of the delineation of the potentially flooded area on both sides of the main channel of the river performed by Andreou et al. [47], by simulating the river hydraulics behavior using HEC-RAS Model in GIS Environment for the main stream channel for about 4 km before the river mouth. Detection of the potentially flooded areas was performed for three different high discharge values that correspond to the peak discharges of the storm events on 22 February 2013, as well as to the discharge values over return periods of 25 and 50 years, respectively. Furthermore, the results of the quantitative geomorphological analysis regarding irregularities in the hierarchical drainage by stream order performed for the Megalo Rema (Rafina) drainage network by Karymbalis et al. [42] were considered.
Finally, to obtain a preliminary assessment of the impacts of the flash flood events on the socio-economic activities of the study area, land use of the catchment was identified utilizing the Urban Atlas 2018 data of Copernicus Land Monitoring Service (see Table 8). Initially, twenty-one land use types were recognized and further grouped into eight categories according to the density of the urban fabric, the sector of economy and the vegetation density. The percentage of each one of these socio-economically and environmentally important land use classes that occupy each flood-hazard level zone of the final map was estimated by overlaying the layer of land use to the basic flood-hazard map.

Factor Classification
Thematic maps in Figures 3, 4 and 5a illustrate the spatial distribution of the parameters' values of slope, distance from stream channels and land cover in the catchment after their classification into five groups according to their impact on flooding.
gently sloping zone, since it is characterized by low (≤2°) and very low (2°-6°) slope values (Figure 3). On the contrary, about 10% of the catchment area has slope values >20° and is characterized as having very low susceptibility to flooding. The classes were assigned these rating values as depicted in Table 4. The highest slopes are found at the northern, and southwestern parts of the study area where the mountainous landscapes of Penteli and Ymittos Mountain exist, as well as at the hilly area along the southeastern water divide. According to the rating adopted in the methodology, the five classes of slope were assigned the rating values as shown in Table 4.
According to the ranking followed for the criterion "distance from the stream channels", riparian areas in immediate proximity to the channel bed of the higher Strahler's steam order of the drainage network are highly and very highly susceptible to flooding. These areas are located along the main channels of Megalo Rema and its major tributaries Valanaris and Krioneri (Figure 4).
Dense and sparse urban areas, which are highly and very highly susceptible to flooding, primarily dominate the western and eastern parts of the catchment, while the south- According to the rating adopted in the methodology, the five classes of slope were assigned the rating values as shown in Table 4. The classification of the geological formations of the catchment into three groups according to their hydrogeological behavior showed that the prevailing geological formations belong to the semi-permeable group ( Figure 6). The northern part of the catchment is dominated by impermeable rocks, whereas the westernmost part, as well as the broader area of the river mouth, is covered by permeable geological formations.    Table 4.
Regarding "slope", the range of the slope values lie between 0 and 44 • . The thematic map of this criterion shows that a significant part of the study area (approximately 54% Sustainability 2021, 13, 10232 13 of 24 of the total catchment area) belongs to the highly and very highly susceptible to flooding gently sloping zone, since it is characterized by low (≤2 • ) and very low (2 • -6 • ) slope values ( Figure 3). On the contrary, about 10% of the catchment area has slope values >20 • and is characterized as having very low susceptibility to flooding. The classes were assigned these rating values as depicted in Table 4. The highest slopes are found at the northern, and southwestern parts of the study area where the mountainous landscapes of Penteli and Ymittos Mountain exist, as well as at the hilly area along the southeastern water divide.
According to the ranking followed for the criterion "distance from the stream channels", riparian areas in immediate proximity to the channel bed of the higher Strahler's steam order of the drainage network are highly and very highly susceptible to flooding. These areas are located along the main channels of Megalo Rema and its major tributaries Valanaris and Krioneri (Figure 4).
Dense and sparse urban areas, which are highly and very highly susceptible to flooding, primarily dominate the western and eastern parts of the catchment, while the southern part of the study area along the main stream channel mainly hosts agricultural activities (Figure 5a). The mountainous parts of the catchment are less susceptible to floods, in terms of the land cover parameter, since they are covered by forests and sparse natural vegetation.
The elevation of the catchment ranges from 0 to 950 m. High-elevation (between 200 m and 500 m, as well as higher than 500 m) appears in the northern and western mountainous parts of the catchment. On the contrary, the lower elevation zones more susceptible to flooding (≤50 m and 50-100 m) are located on both sides of the main stream channel at the eastern and southern portions of the catchment (Figure 5b).
The classification of the geological formations of the catchment into three groups according to their hydrogeological behavior showed that the prevailing geological formations belong to the semi-permeable group ( Figure 6). The northern part of the catchment is dominated by impermeable rocks, whereas the westernmost part, as well as the broader area of the river mouth, is covered by permeable geological formations.  Table 4.

Multi-Criteria Decision Analysis (MCDA)-Flood-Hazard Assessment
The thematic maps of the spatial distribution of the five factors (Figures 3-6) were used as inputs in the MCDA. All values of all the parameters involved are classified into five classes and assigned to a uniform suitability five-grade rating scale. Higher classified

Multi-Criteria Decision Analysis (MCDA)-Flood-Hazard Assessment
The thematic maps of the spatial distribution of the five factors (Figures 3-6) were used as inputs in the MCDA. All values of all the parameters involved are classified into five classes and assigned to a uniform suitability five-grade rating scale. Higher classified ranking values correspond to areas more susceptible to floods, while lower values correspond to less flood-prone areas. Table 4 includes the factors/parameters, their classes and their ratings. All the parameters are assigned with ratings from 0 to 4, whereas the geological formation classes are assigned with ratings of 1 to 3, since there are no ideally permeable or impermeable geological formations. For this reason, the ratings 0 and 4 were not selected.
The results of the AHP procedure, the extraction of criteria weights and the calculation of the CR are shown in Table 5. The calculated CR (Table 5) in this study is 0.019 (lower than the threshold 0.1) within the range for the acceptance of the consistency of judgments in the pairwise comparison matrix. Consequently, the weights' consistency is affirmed.
As described in the methodology, the flood-hazard level values were calculated following the WLC procedure, using the classified and normalized values of the five criteria involved and their weighting coefficients (Tables 4 and 5). The final basic floodhazard assessment map (Hbasic) is presented in Figure 7a. It was produced after the reclassification of the flood-hazard level values into five categories (from very low to very high) with the quantile classification method. In addition, two other maps characterizing the maximum (Hmax = Hbasic + ∆S) and minimum (Hmin = Hbasic + ∆S) flood-hazard level values were produced (Figures 7b and 8) after the examination of the influence of uncertainty of the adopted factor weights on the flood-hazard assessment. This was necessary since uncertainty plays an important role in natural hazard evaluation [70,71] and can bias the outcome of every hazard assessment. The uncertainties of the weighting coefficient for each factor (∆W i ) are shown in Table 6.   The resultant flood-hazard assessment maps indicate that the areas of very high and high flood-hazard are distributed mostly on the lower reaches of the Megalo Rema River, along the main stream channel and around the mouth of the River (Figures 7a,b and 8). Additionally, an extensive low-lying area at the southern and western parts of the catchment belongs to the zones of high and very high flood-hazard risk. The flood-hazard map shows that the urban areas on both sides of the main channels of the drainage network are more prone to flooding. It is evident that these parts of the catchment are relatively low-lying, of gentle slope and densely populated. They include residential areas of the Rafina-Pikermi, Penteli, Pallini, Peania, Spata, Glika Nera, Gerakas and Anthousa settlements, with a total population of ~85,000 residents. This high population density is directly associated with the urbanization of this area, which increases the surface runoff during extreme rainfall events and contributes to an increased human exposure and social vulnerability to the flood-hazard. The northern and the southwestern mountainous parts of the study area are low and very low flood-hazard areas. Limited regions at the central part of the catchment are classified as being of high and very high flood risk; those that are include the settlements of Dion, Kallitechnoupolis and Drafi.  Regarding the spatial development of the five flood-hazard zones, the percentages corresponding to the area of each hazard zone for the three flood-hazard assessment maps (Hbasic, Hmax and Hmin) of Figures 7 and 8 are included in Table 7. An area of 26.3 km 2 , which corresponds to the 22.7% of the total catchment, belongs to the high flood-hazard zone, while 24.7 km 2 (some 21.3% of the catchment) is very highly prone to flooding. The areas of the high flood-hazard zones of the maximum flood-hazard assessment value map (Hmax) are negligibly increased compared to the corresponding zones of the basic map (Hbasic), whereas the areas of the very high flood zone are negligibly lower. This similar comparison proves that the spatial extents of moderate hazard zones decrease, while the moderate, low and very low hazard zone increase. Comparing the zones of the map of the minimum value of the flood-hazard assessment (Hmin) to those of the basic map (Hbasic), the areas of high and very high flood-hazard zones decrease, whereas the moderate, low and very low hazard areas slightly increase.
The resultant flood-hazard assessment maps indicate that the areas of very high and high flood-hazard are distributed mostly on the lower reaches of the Megalo Rema River, along the main stream channel and around the mouth of the River (Figures 7a,b  and 8). Additionally, an extensive low-lying area at the southern and western parts of the catchment belongs to the zones of high and very high flood-hazard risk. The floodhazard map shows that the urban areas on both sides of the main channels of the drainage network are more prone to flooding. It is evident that these parts of the catchment are relatively low-lying, of gentle slope and densely populated. They include residential areas of the Rafina-Pikermi, Penteli, Pallini, Peania, Spata, Glika Nera, Gerakas and Anthousa settlements, with a total population of~85,000 residents. This high population density is directly associated with the urbanization of this area, which increases the surface runoff during extreme rainfall events and contributes to an increased human exposure and social vulnerability to the flood-hazard. The northern and the southwestern mountainous parts of the study area are low and very low flood-hazard areas. Limited regions at the central part of the catchment are classified as being of high and very high flood risk; those that are include the settlements of Dion, Kallitechnoupolis and Drafi. The distribution of land use in the susceptible to flooding zones of the Hbasic flood hazard map is illustrated in the bar charts of Figure 9. Accordingly, 30.6% and 37.0% of the high and very high flood-hazard zones are occupied by agricultural land and related activities. Significant parts of the highly prone to flooding areas correspond to continuous and discontinuous dense urban fabric (12.6%), discontinuous medium-and low-density urban fabric (15.1%) and industrial, commercial, public and private units and construction sites (7.2%). Similarly, dense urban fabric, medium-to low-density urban fabric and industrial, commercial, public and private units and construction sites constitute 9.4%, 13.5% and 10.1% of the very high flood risk zone, respectively. It is worth noting that some 9.9% and 9.5% of the high and very high flood-prone areas, respectively, consist of road network and associated land. Very low to moderate prone areas appear mainly at forests and areas covered by natural vegetation (see Table 8). Discontinuous very low density urban fabric (S.L.: <10%), isolated structures, green urban areas 2.52 Discontinuous medium-and low-density urban fabric (S.L.: 10-50%) 13.57 Continuous and discontinuous dense urban fabric (S.L.: >50%) 9.69  The analysis performed by Giannaros et al. [48] in an attempt to introduce an operational impact-based warning system in the area, coupling a state-of-the-art numerical weather prediction model with an advanced, spatially explicit hydrological model, provided some preliminary discharge thresholds in terms of flash flooding's socioeconomic impacts. Minimal impact was reported when the peak discharge in the mouth of Megalo Rema was lower than 20 m 3 /s, whereas maximum stream discharges that ranged from 20 m 3 /s to 40 m 3 /s were associated with major impact. Significant impact was induced by events that were characterized by peak stream flows higher than 40-60 m 3 /s. These thresholds are considered as preliminary since their reliability is mainly moderated due to the lack of discharge data along the entire length of the Megalo Rema stream [48].
The flood-hazard assessment map provides valuable information for land-use planning at a regional scale, leading to the determination of safe and unsafe areas for urban development [72,73]. Even though hazard maps represent a snapshot situation regarding flood-hazard level variation within the catchment, they could still support policy makers with knowledge for future planning. This is very important, particularly for the study area, since the land cover properties are constantly changing. This can be attributed primarily to two factors: forest fires that devastated a significant part of the forested land of the area and the increased urbanization rate of the area, especially during the last 40 years. The increase in the population of the Municipality of Rafina-Pikermi and Penteli for the period between 1981 and 2011 was estimated up to 70% and 60%, respectively (Source: Hellenic Statistical Authority). One of the main reasons for this positive population growth is the construction of several public works in this area (i.e., the new international airport of Athens in Spata, the Attiki Odos motorway and the developing Rafina port), which contribute to a significant increase in private building activity in the neighboring settlements and Municipalities [74].

Verification of the Flood-Hazard Map's Accuracy
The accuracy of the basic flood-hazard map was verified by means of the frequency ratio and flood events which affected the study area over the past 30 years. The findings established that the vast majority of the flood events, almost 89% (40 out of 45 incidents), occurred within the limits of the high and very high flood-hazard zones ( Table 9). The frequency ratio is higher than one in the high and very high flood-hazard zones, which is indicative of the strong positive correlation between the hazard zones and the flood occurrences, since the flood-hazard intensity raises as the frequency ratio values increase (Table 9). Another way of validating the results of the final flood-hazard map is by overlaying the outcome results of the application of the one-dimensional steady flow hydraulic model HEC-RAS and comparing the potentially flooded areas with the borders of the high and very high flood-hazard zone. The flood-hazard mapping is in accordance with the results of the delineation of the potentially flooded area on both sides of the main stream channel, by simulating the river hydraulics' behavior using the HEC-RAS Model in a GIS environment for the last 4 km of its course before the river mouth ( Figure 10) [47]. The detection of potentially flooded areas was performed for three different high discharge values that correspond to the peak discharges of the storm events on 22 February 2013, as well as to the discharge values over return periods of 25 and 50 years, respectively. According to the results of the hydraulic model, the potentially flooded areas coincide with the highly and very highly flood-hazard prone areas, as mapped in this study. The results of a HEC-RAS simulation, taking into account higher discharges (i.e., a 100-year return period), would show that more of the high and very high flood-hazard zone could be included in the flood's extent, possibly including the flood event that lies in the moderate flood-hazard zone and outside the flooding extent of the simulation results in Figure 10.
The resultant flood-hazard map is also in line with the quantitative geomorphological analysis performed for the Megalo Rema drainage network by Karymbalis et al. [42]. The results of the aforementioned geomorphological analysis show irregularities in the hierarchical drainage by stream order that enhance flash floods. The most significant among them concerns the fourth order streams that drain directly into the sixth order main stream channel of the Megalo Rema (Rafina). The reason for this is the Neos Voutzas stream which joins the main channel at its lower reaches, enhancing its discharge when high precipitation events occur at the upper reaches of this tributary. This assumption is supported by the flood-hazard map of the present study, in which the broader area of the confluence of these two streams belongs to the very highly prone to flooding zone.
According to the results of the hydraulic model, the potentially flooded areas coincide with the highly and very highly flood-hazard prone areas, as mapped in this study. The results of a HEC-RAS simulation, taking into account higher discharges (i.e., a 100-year return period), would show that more of the high and very high flood-hazard zone could be included in the flood's extent, possibly including the flood event that lies in the moderate flood-hazard zone and outside the flooding extent of the simulation results in Figure  10. The resultant flood-hazard map is also in line with the quantitative geomorphological analysis performed for the Megalo Rema drainage network by Karymbalis et al. [42]. The results of the aforementioned geomorphological analysis show irregularities in the hierarchical drainage by stream order that enhance flash floods. The most significant among them concerns the fourth order streams that drain directly into the sixth order main stream channel of the Megalo Rema (Rafina). The reason for this is the Neos Voutzas stream which joins the main channel at its lower reaches, enhancing its discharge when high precipitation events occur at the upper reaches of this tributary. This assumption is supported by the flood-hazard map of the present study, in which the broader area of the confluence of these two streams belongs to the very highly prone to flooding zone.

Conclusions
The aim of the present paper was the development of a catchment scale methodology to identify and map the zones prone to flooding in the fast-growing urban catchment of Megalo Rema in East Attica, Greece. The proposed methodology was based on MCDA integrated with GIS (GIS-based MCDA) using AHP and took into account various critical factors such as slope, distance from stream channels, land cover, elevation and geology.

Conclusions
The aim of the present paper was the development of a catchment scale methodology to identify and map the zones prone to flooding in the fast-growing urban catchment of Megalo Rema in East Attica, Greece. The proposed methodology was based on MCDA integrated with GIS (GIS-based MCDA) using AHP and took into account various critical factors such as slope, distance from stream channels, land cover, elevation and geology.
The flood-hazard assessment map produced showed that 51 km 2 , which corresponds to 44% of the total catchment area, belong to the high and very high flood-hazard classes. The areas most prone to flooding are the low-lying, gently sloping and densely urbanized southern and western parts of the catchment, as well as the area of the river mouth. The areas highly prone to flooding host socioeconomically significant land uses such as continuous and discontinuous dense urban fabric (12.6%), discontinuous medium-and low-density urban fabric (15.1%) and industrial, commercial, public and private units and construction sites (7.2%). Similarly, the zone of very high flood risk is occupied by dense urban fabric (9.4%), medium-to low-density urban fabric (13.5%) and industrial, commercial, public and private units and construction sites (10.1%). In addition, 9.9% and 9.5% of the high and very high flood-prone zones of the catchment consist of the road network and its associated land, while 30.6% and 37.0% of these zones are occupied by agricultural land and related activities.
A significant correlation between the flood-hazard zones and the spatial distribution of past flood phenomena was detected, calculating the frequency ratio of the flood events as shown in Table 9. The flood-hazard map is also in accordance with the results of the quantitative geomorphological analysis of the Megalo Rema drainage network, as well as with the potentially flooded areas delineated by simulating the hydraulic behavior of the lower reaches of the river using the HEC-RAS Model. These verifications demonstrated the reliable results and high accuracy achieved by the created flood-hazard map.
The proposed methodology is simple and provides a tool for the improvement of flood risk management strategies and action plans, while it can also be useful in land use planning projects at a catchment scale.