Vulnerability Analysis of Land Instability Using Multi-Criteria Evaluation for Urban Sustainability: Methodological Overview and Case Study Assessment

There were many developed urban areas have been established without well studied from the hazards perspective. However, Mokattam Plateau is one of the most vulnerable areas of frequent landslides and rock-falls disasters (Moustafa et al., 1991). So, an integrated analysis method is highly needed. Spatial analysis of Multi-criteria evaluation (MCE) provided an assessment of the hazards’ factors (i.e. faults, joints, lithology, slope, old wadies/surface drainage, and quarries) at Mokattam area. The data have been analyzed by the MCE, Analytical Hierarchy Process (AHP)/Ordered Weighted Average (OWA)/Weighted Linear Combination (WLC). The research found that the geological factors (faults, joints, and lithology) were the highest contributors by about 65% of the hazardous driving forces, while the geomorphological processes (slope and old wadies/surface drainage) were contributed by about 30%. In addition, the impact of the human activities such as random urbanization, excess use of irrigation water and the transportation are critical hidden drivers that affect the land instability and accelerates the landslides and rock-falls (Amasha, 2009). Therefore, the decision makers and urban planners have to consider the four scenarios of low risk-high tradeoff (MIDAND), and high risk-some tradeoff (MIDOR) in their disaster risk management plans. While the risk-taking (OR) option is highly recommended for the new urban development projects to ensure the sustainability and risk resilience. While the risk-averse (AND) scenario is not recommended.


Introduction
Egypt has rapid population growth that reached 92.1 Million inhabitants (inside the country) in 2017 and 90.1 in 2016 compared with 71.3 in 2006 by increase rate about 20% through 10 years (CAPMAS, 2017). Such overpopulation growth requires more residential housing and urban development to accommodate population. The Egyptian government has encouraged the urban development projects outside the Nile Delta and flood plain to rescue the fertile land for agriculture and food security. According to this governmental plan, there are many new cities have been developed over the last two decades for residential and industrial purposes.
In the case of Greater Cairo, most of the urban expansion was normally slummy adjacent to the old urbanized areas. It is however, expanded to the desert with more modern urban planning. The critical issue is that, most of these new cities and urbanization have not thoroughly investigated in respect of determining the land stability and geo-environmental hazards. Unfortunately, there is potential disaster, related to the bedrock instability and compliance with urbanization in addition to flood hazard in other areas. Therefore, there is an urgent need to evaluate these areas before any development in order to reduce the number and the frequency of hazards and to mitigate their risk.
The GIS modeling is a potential approach for multiple analyses to estimate different scenarios of site hazard/suitability. Such approach enables decision makers to obtain alternative developmental strategies to be examined and compared before making approval of the implementation plans (Albrecht, 1996). GIS has been widely used in assessing geologic hazards, site selection, and geo-environmental evaluation (Lee & Pradhan, 2006;. The GIS techniques and statistical models are used in evaluating landslide hazards (Carrara et al., 1991), in addition to risk analysis (Pradhan et al., 2009;Pradhan & Lee, 2009a, 2009bCarver, 1991), integrated multi-criteria evaluation and GIS with geological information. It had been developed a method to combine mapping and engineering properties using a GIS to evaluate natural hazards (Gunawan et al., 1992).
The AHP is a multi-criteria decision-making approach introduced by (Saaty, 1977), in which a series of pair-wise comparisons between factors have developed to create a scaled set of preferences. It also well documented in (Saaty, 1996) and (Saaty, 2003) literatures. The added-value of this multi-criteria model is considering the weights and cross correlation of each factor relative to every other factor.
It is quite difficult to define precisely MCE Analysis. However, various definitions appear in literature. One common definition is that of (Roy, 1996), who

Area of Study
The area of study-Mokatam Plateau Area (MPA) is located eastern of the Cairo built up area, Figure 1. It characterizes by the informal urbanization over and down the plateau where it recorded the highest population density in Cairo, more than 25,000 inhabitants in the sq. km (Amasha, 2009). It is considered as one of the most vulnerable areas for natural hazards with major impact on so-

Vulnerability of the Area of Study to Natural Hazards
Mohatam plateau area (MPA) is considered a common vulnerable area for landslides and natural hazards. Few critical factors are driving its sensitivity to natural hazards, such as atmospheric conditions that critically play a key role in the exogenic driving force (Amasha, 2009). The system in equilibrium can be defined by real parameters that represent the limits of the balance condition, whereas the limits of equilibrium have critical conditions called thresholds. The threshold conditions develop in response to gradual changes within the system. In many cases, threshold represents a deterioration of resistance rather than an increase in the driving forces (Schumm, 1973;Ritter, 1988).
Regarding to the lithological aspects, there are five claystone beds, 5 to 10 meters thickness, existing between Upper Eocene Ain Musa Bed and Middle Eocene limestones inter-bedded with marl. The claystone, especially the uppermost beds, are characterized by flow structures. This claystone beds has the highest swelling potential due to the occurrence of montmorillonite clay mineral (Moustafa et al., 1991). Topographically, the relative high cliffs of the studied area (223 m a.m.s.l.), the slope stability plays a key factor that controlling the mass movements and slope failure. It is evaluated based on the other factors of lithology, geologic structure, differential physical weathering and the topography. Slope stability, therefore, represents some balance between driving forces (shear stresses) and resisting forces (shear strength), this can be expressed as a safety ratio (Ritter et al., 2001): Theoretically, failure occurs when F = 1. Beside the natural driving forces affecting the high potential risk on MPA, the impact of human activities, is an important factor which accelerates the land instability. The human activities include the urbanization on the highly-potential risk foundation bed rocks, especially of Upper Eocene rocks of the Upper Plateau. The leakage of sewage water, the irrigation of the urban green areas and the limestone quarrying pumps are considered as important factors in the land instability of the Mokattam Plateau area (Amasha, 2009).

Materials and Methodology
To apply the GIS modelling and analysis, an initial step was to prepare geographic layers. Six GIS layers that represent hazards factors were prepared and classified into four sub-groups as follows: 3) Geomorphology (slope as topography and old wad is/surface drainage) 4) Human activities (old quarries) in addition to the urbanization settlements.
ArcGIS Desktop package was used for spatial and layer preparation, IDRISI package was used for MCE Model and make statistical and mapping analysis via the following steps: • GIS functions were applied to create distance layers of faults, joints, drainage, slope and lithology. This followed by creating proximity layers for faults, joints, surface drainage (wadies) and quarries.
• Weighting scale was applied to lithology layer from 1 -9 scale according to the effect on hazards. However, slope was scaled as gradient in degrees.
Therefore, the vulnerability of slope as a hazard factor is clearly observed at the boundaries of the Plateaux where cliffs occur.
• All layers were aggregated to 0 -255 scaled before MCE analyses.
• Define the hazard assessment as single objective for the MCE model; this was functioned by overlaying of several thematic layers to identify zones that are common to given criteria. Furthermore sieve mapping is used to examine all layers to identify zones that encompass all desired criteria (Eastman, 2006).
To meet a specific objective, it is frequently the case that several criteria will need to be evaluated. Such a procedure is called Multi-Criteria Evaluation (Voogd, 1983;Carver, 1991).
Due to the existence of different scales of the criteria that might create uncertainty, it is necessary that all criteria factors are standardized before combination. A review of variety of procedures for standardization typically had been used the minimum and maximum values as scaling points (Carver, 1991). The simplest approach that adopted is a linear scaling (Equation (2) In IDRISI, the FUZZY module is provided for the standardization of factors using a whole range of fuzzy set membership functions (Eastman & Jiang, 1996).
The module provides the option of standardizing factors to either a 0 -1 real number scale or a 0 -255 byte scale. This latter option is recommended because the MCE module has been optimized for speed using a 0 -255 level standardization. Importantly, the higher value of the standardized scale must represent the case of being more likely to belong to the decision set.
The implementation of this standardization approach, for the used criteria, aggregations were as the following processes: • Fault-distance layer is processed as Sigmoidal function type and Monotonically decreasing function shape; • Joints-distance layer is processed as Sigmoidal function type and Monotonically decreasing function shape; • Lithology-distance layer is processed as Sigmoidal function type and Monotonically decreasing function shape; • Slope-distance layer is processed as J-shaped function type and Monotonically increasing function shape; • Wadies-distance layer is processed as J-shaped function type and Monotonically decreasing function shape; • Quarries-distance layer is processed as Linear function type and Monotonically decreasing function shape; • The AHP pair-wise comparison matrix between the standardized 6 factors is edited to produce Eigen vector values. These Eigen values considered as weights of factors.
Finally, based on the pair-wise comparison matrix that had developed by (Saaty, 1980), an estimated weight is calculated to derive a consistency ratio (CR) to produce eigenvector values in terms of weights. The calculated weights of the applied factors in our work are: Faults 0.0914, Joints 0.2227, Lithology 0.3343, Slope 0.1741, Wadies 0.1234 and Quarries 0.0541. The (CR) in our work is calculated as 0.08 which is less than 0.10 as recommended. So, it is accepted.

The Simple Additive Weighting (SAW) or Weighted Linear Combination
(WLC) is the most often used technique in multi-criteria decision making. Criteria here may include weighted factors and constraints. Calculating the product of weight and factor multiplied with all constraints at any location, and then summing up all products yields a total overall score. The score for each alternative A is shown as in Equation (3) (Eastman, 2006): If a constraint is part of the decision, wi = weight of factor i, xi = criterion score of factor i, cj = criterion score of constraint j. The OWA may also be incorporated into the criteria aggregation process. The order weights are applied pixel by pixel to the order of suitability scores. The effect of order weights is most easily understood in terms of levels of risk and tradeoff (Yager, 1988). Factor Weights; Specify weights for factors to produce an output of a decision map using the user-defined weights by applying the pair-wise comparison matrix (Saaty, 1977). This has finally created decision maps using the AHP weights and the ordered weighted averaging (OWA), Table   1 and Figure 3.
With a weighted linear combination (WLC), factors are combined by applying a weight to each, followed by summation of the results to yield a suitability map (Voogd, 1983), i.e.:

S wixi
where S = suitability, wi = weight of factor i, xi = criterion score of factor i.  In cases where Boolean constraints also apply; the procedure can be modified by multiplying the suitability calculated from the factors by the product of the constraints (Voogd, 1983), i.e.: where cj = criterion score of constraint j, Π = product.

MCE and Ordered Weighted Average (OWA)
The outcomes of the MCE are categorized into different scenarios according to the options have been taken between the risk and trade off. The order weights that generated these decision maps are scaled from (risk-averse, minimum, AND) to (risk-taking, maximum, OR)- (Figure 4).
This has resulted in kind of decision maps that support the decision makers and planners. Each decision map has categorized into equal interval-five classes based on the aggregated scale (0 -255). These four scenarios are: 1) The first scenario is based on applying the minimum risk throughout the risk-averse (AND) model where there is no tradeoff and minimum risk. It considers the minimum weights of the driving factors generating the hazards. . This is actually representing nearly 90% of the study area, Figure 9 and Table 2.
2) The second scenario is the maximum risk is applied in the risk-taking (OR) model. It considers the maximum weights of the driving factors contribute to the hazards. Figure 6 shows the decision map resulted from this scenario. It shows that the aggregated hazard weights are distributed from 173 to 255.
This scenario shows that the majority of the study area is dominated by class 5 (extremely hazards). This means that 99% of the study area is under high risk of hazards from this scenario (Figure 9, Table 2).
3) The third scenario is the alternatives between low risk and high tradeoff where low risk-high trade off (MIDAND) model is applied. It considers the median order weights of the no-and-full trade off and the risk-averse (AND) values of the factors affecting the hazards. Figure 7 shows the decision map resulted from this scenario, where the hazard weights are scaled from 26 to 250. This scenario has categorized the study area into various hazards categories. Nearly 50% of the study area is in class 4 (high hazards), nearly 20% is located in class 3 (moderate hazards), while about 15% is located in class 2 (low hazards), and 13% in class 1 (very low hazards), Figure 9 and Table 2.
This means that more stable areas of classes 1 and 2 are shown in this scenario.  Figure 8 shows the outcomes of this scenario, the decision map shows a scaled hazard weight values from 133 to 255. This scenario has also categorized the area of study into different hazards zones. The area of study is fully hazards area, the majority (75%) is considered extreme hazard (class 5), however only 25% is high risk (class 4) Figure 9 and Table 2. This means that the more stable areas of classes 1 and 2 are totally disappeared in this scenario while the more hazardous classes are dominated.

Conclusions
Six     This is followed by the high risk-low tradeoff (MIDOR) scenario with some conditions to be considered. The limitations and conditions will be increased when using the low risk-high tradeoff (MIDAND) scenario, and even further constrains for the risk-averse (AND) scenario, which are minimizing the contribution of the risk factors in the MCE-OWA model. Figure 10 shows disaster, increase preparedness for response and recovery, and thus strengthen resilience". Therefore, the decision-makers should be taken into consideration the factors accelerating the hazards and act to mitigate the risks by applying the disaster risk management plans at the most vulnerable areas. By applying this methodology to evaluate an existing built up area, it is not easy to decide which scenario to apply. When I made a mask of urban area polygons on the output maps of the four scenarios, I found that almost 100% of the urbanization is located in the class 5 in the scenarios risk-taking (OR) and high risk-low trade off (MIDOR). While about 75% of the urbanization is located in class 4 in the scenario (MIDAND) and about 93% in class 3 in the scenario (AND), Figure 9. Therefore, the decision here is very difficult to maintain the risky situation of such urban areas.
Finally, it is recommended to apply the MCE (AHP-OWA-WLC) methodology, for its integrated and comprehensive perspective, to achieve more resilient and sustainable development of the new urban development projects. The most suitable scenario to evaluate these projects is the risk-taking (OR) then the high risk-low trade off (MIDOR) models. But for stability evaluation of our study area (MPA), the decision makers need to interactively use this model with updated information to ensure optimum accuracy and safer scenario.