Assessment of Geohazards and Preventative Countermeasures Using AHP Incorporated with GIS in Lanzhou, China

: This research investigates geohazards and preventative countermeasures for Lanzhou City, China. To investigate the factors related to the development of geohazards in Lanzhou, the regional geological conditions around Lanzhou were investigated. The geomorphology of the region is comprised of a loess landform underlying quaternary loess deposits. A large number of faults induced by strong neotectonic movements are present in the area. Therefore, earthquakes frequently occur around Lanzhou. Earthquakes cause numerous rock falls and landslides, with landslide masses found scattered on the upper middle level of the area’s mountains. When intense rainfall occurs, a lot of loose deposits are brought together along steep gullies, forming debris ﬂows; hence, a disaster chain of earthquake–landslide–debris ﬂow is formed. To evaluate the georisks around Lanzhou, the Analytic Hierarchy Process (AHP) was employed to assess geohazards. The spatial distribution of the evaluated georisk levels was mapped using a Geographic Information System (GIS). Based on the assessed results, about 55% of the urban area and 44% of Gaolan County have high or very high-risk levels. The ratio of relatively high disaster risk levels is up to 31% of the total area. To mitigate these geohazards, both strategic decision making and technical countermeasures should be implemented.


Introduction
In recent years, many cities in China have been subjected to extreme weather hazards due to climatic changes [1][2][3][4][5]. Recent extreme weather events have included flooding hazards in Guangzhou City, Guangdong Province [1] and tornado hazards in Yancheng City, Jiangsu Province [2]. Large-scale geohazards in the northwestern region of China may also occur due to its special geographical location and semi-arid climatic conditions. In the northwestern region of China, geohazards (for example, earthquake-landslide-debris flow and flooding) frequently occur [6][7][8][9][10][11][12]. Figure 1 shows the geographical environment whereby uplift of the Tibetan Plateau induces strong neotectonic movement and the development of a large number of faults in the northwestern region. In the northwestern parts of the region, three mountain ranges (the Altay, Tianshan, and Kunlun Mountains) can be found. Two low-lying desert basins (the Junggar and Tarim Basins) are also located between these three mountain ranges. The Tibetan Plateau, the largest plateau in the world, is located at the southwestern region. The monsoon areas of eastern China, which contain plains and hills at lower elevations, are located in the southeastern region. As shown in Figure 1, Lanzhou City is located in a complex geological environment, surrounded by active faults, normal faults, and reverse faults. hills at lower elevations, are located in the southeastern region. As shown in Figure 1, Lanzhou City is located in a complex geological environment, surrounded by active faults, normal faults, and reverse faults. With economic development, the urban scale of Lanzhou city is rapidly expanding. Numerous construction activities such as tunnel construction and excavation with ground water pumping have severely disturbed the natural geological environment [6]. To mitigate the influence caused by anthropic activities,   [13] proposed a serious of ground improvement techniques for implementation in soft soil deposits [14][15][16].   [17,18] also proposed a new method for predicting geohazards in the process of tunnel construction. Geohazards caused by construction activities [19][20][21] and other anthropic activities, such as deforestation for agricultural use and artificial-irrigation-induced landslides, are increasing in frequency [22,23]. Anthropic activities have caused numerous geological disasters over the years [24,25]. According to statistics, since 1949, geological disasters have caused more than 670 deaths and economic losses of 776 million Chinese RMB ($100 M USD) in Lanzhou [26].
Lanzhou, the capital city of Gansu Province, is located inland of the northwest region, which is the geometric center of China. Lanzhou City has frequently experienced strong tectonic movements and earthquakes [27,28]. The occurrence of an earthquake can induce landslides accompanied by debris flows [10,11]. Lanzhou also sits on loess deposits, which result in loess landslides under heavy rainfall, since loess is collapsible and loose in structure [8,28]. Intense short-duration rainfall is concentrated from July to September, and can easily induce debris flow and landslides because of the area's loose loess deposits. Moreover, mountains with rocky formations and geological faults are found in Lanzhou City. These mountains are in a weathered state and are covered by sedimentary rocks, also providing material for debris flows and landslides [26]. Geohazards readily occur in this region under the combined effects of complex geological formations, steep landforms, seasonal heavy rainfall, and anthropogenic activities. Therefore, it is urgent to perform a geohazard risk assessment for future sustainable development in and around Lanzhou.
There are three typically used methods of risk assessment: statistical analysis based on historical data [29][30][31]; multi-criteria indexing [32,33]; and geographic information system (GIS) and remote sensing technology [34,35]. Black and Burns (2002) [29] have presented an overview of changes in the estimation of Scottish river flood risks over time by reassessing flood records. Werritty (2013) [30] has applied trend analysis and current climatic change scenarios to identify issues in water resourcing. Li et al. (2014) [31] have argued that long-term historical flood materials are one of the most critical references for flood risk evaluation. Yet though assessment results can reflect the risks associated with a given study area, this approach requires a large amount of data and cannot accurately reflect spatial variation. Jiang et al. (2009) [32] and Kazakis (2015) [33]  With economic development, the urban scale of Lanzhou city is rapidly expanding. Numerous construction activities such as tunnel construction and excavation with ground water pumping have severely disturbed the natural geological environment [6]. To mitigate the influence caused by anthropic activities,   [13] proposed a serious of ground improvement techniques for implementation in soft soil deposits [14][15][16].   [17,18] also proposed a new method for predicting geohazards in the process of tunnel construction. Geohazards caused by construction activities [19][20][21] and other anthropic activities, such as deforestation for agricultural use and artificial-irrigation-induced landslides, are increasing in frequency [22,23]. Anthropic activities have caused numerous geological disasters over the years [24,25]. According to statistics, since 1949, geological disasters have caused more than 670 deaths and economic losses of 776 million Chinese RMB ($100 M USD) in Lanzhou [26].
Lanzhou, the capital city of Gansu Province, is located inland of the northwest region, which is the geometric center of China. Lanzhou City has frequently experienced strong tectonic movements and earthquakes [27,28]. The occurrence of an earthquake can induce landslides accompanied by debris flows [10,11]. Lanzhou also sits on loess deposits, which result in loess landslides under heavy rainfall, since loess is collapsible and loose in structure [8,28]. Intense short-duration rainfall is concentrated from July to September, and can easily induce debris flow and landslides because of the area's loose loess deposits. Moreover, mountains with rocky formations and geological faults are found in Lanzhou City. These mountains are in a weathered state and are covered by sedimentary rocks, also providing material for debris flows and landslides [26]. Geohazards readily occur in this region under the combined effects of complex geological formations, steep landforms, seasonal heavy rainfall, and anthropogenic activities. Therefore, it is urgent to perform a geohazard risk assessment for future sustainable development in and around Lanzhou.
There are three typically used methods of risk assessment: statistical analysis based on historical data [29][30][31]; multi-criteria indexing [32,33]; and geographic information system (GIS) and remote sensing technology [34,35]. Black and Burns (2002) [29] have presented an overview of changes in the estimation of Scottish river flood risks over time by reassessing flood records. Werritty (2013) [30] has applied trend analysis and current climatic change scenarios to identify issues in water resourcing. Li et al. (2014) [31] have argued that long-term historical flood materials are one of the most critical references for flood risk evaluation. Yet though assessment results can reflect the risks associated with a given study area, this approach requires a large amount of data and cannot accurately reflect spatial variation. Jiang et al. (2009) [32] and Kazakis (2015) [33] evaluated regional flood risk instead by establishing a multi-criteria indexing system. However, this method has a few deficits in determining subjective indices. GIS and remote sensing method have been adopted for flood risk assessment in urban areas; these two technologies also provide technological support for inundation Sustainability 2018, 10, 304 3 of 21 risk evaluation [34,35]. Among these presented approaches, Analytic Hierarchy Process (AHP) is a useful tool for providing qualitative and quantitative assessment to decision makers. GIS is a powerful information integration system, which can be used to map the evaluated results of geohazards. Yi et al. (2010) [36] and Ayalew et al. (2005) [37] have developed a GIS-based framework for disaster prevention in Korea using the AHP method. Arisido et al. (2017) [38] and Zhang et al. (2009) [39] used Bayesian networks combined with GIS technology to assess flood risks. For our study, the AHP method combined with GIS technology was used to evaluate geohazard risks.
The objectives of this paper are to investigate potential geohazards in Lanzhou, to evaluate geohazard risk levels for the various districts in Lanzhou City, and to propose appropriate countermeasures to mitigate future catastrophic disasters. Figure 2 shows the administrative map of Lanzhou, which includes five districts and three counties. The five districts constitute the urban area of Lanzhou: Chengguan, Qilihe, Xigu, Anning, and Honggu. The three counties are Yongdeng, Yuzhong, and Gaolan. The land area of Lanzhou is 13,085.6 km 2 and contains a population of 3.6 million people, and the urban area is located within Lanzhou Basin, one of the largest basins in the middle of the Yellow River, located on the northeast edge of the uplifting Tibetan Plateau. The basin extends 40 km from northwest to southeast and is 10 km wide. It is divided into two parts, an eastern basin and a western basin, by a folding fault between Gaolan Mountain to the south and Baita Mountain to the north. In the rising Tibetan Plateau, river terraces are usually regarded as topographic evidence of tectonic movement.

Study Area
Sustainability 2018, 10, x FOR PEER REVIEW 3 of 21 evaluated regional flood risk instead by establishing a multi-criteria indexing system. However, this method has a few deficits in determining subjective indices. GIS and remote sensing method have been adopted for flood risk assessment in urban areas; these two technologies also provide technological support for inundation risk evaluation [34,35]. Among these presented approaches, Analytic Hierarchy Process (AHP) is a useful tool for providing qualitative and quantitative assessment to decision makers. GIS is a powerful information integration system, which can be used to map the evaluated results of geohazards. Yi et al. (2010) [36] and Ayalew et al. (2005) [37] have developed a GIS-based framework for disaster prevention in Korea using the AHP method. Arisido et al. (2017) [38] and Zhang et al. (2009) [39] used Bayesian networks combined with GIS technology to assess flood risks. For our study, the AHP method combined with GIS technology was used to evaluate geohazard risks. The objectives of this paper are to investigate potential geohazards in Lanzhou, to evaluate geohazard risk levels for the various districts in Lanzhou City, and to propose appropriate countermeasures to mitigate future catastrophic disasters. Figure 2 shows the administrative map of Lanzhou, which includes five districts and three counties. The five districts constitute the urban area of Lanzhou: Chengguan, Qilihe, Xigu, Anning, and Honggu. The three counties are Yongdeng, Yuzhong, and Gaolan. The land area of Lanzhou is 13,085.6 km 2 and contains a population of 3.6 million people, and the urban area is located within Lanzhou Basin, one of the largest basins in the middle of the Yellow River, located on the northeast edge of the uplifting Tibetan Plateau. The basin extends 40 km from northwest to southeast and is 10 km wide. It is divided into two parts, an eastern basin and a western basin, by a folding fault between Gaolan Mountain to the south and Baita Mountain to the north. In the rising Tibetan Plateau, river terraces are usually regarded as topographic evidence of tectonic movement.

Topography
Topography has a critical influence on the development of geohazards. For example, the distribution of flooding hazards is significantly affected by variations in spatial topography [40]. Flat areas in particular are prone to flood. River networks play an important role in the reduction of floodwaters during flood events. The average altitude of China's northwestern region ranges from 2000 to 3000 m. The Yellow River, which starts on the Tibetan Plateau in Qinghai Province, passes through the provinces of Qinghai, Sichuan, Gansu, Ningxia, Inner Mongolia, Shanxi, and Henan. Figure 3 illustrates the altitude and distribution of the river system in Lanzhou. The altitude of Lanzhou ranges from 3796 to 1423 m with an average of approximately 2500 m. The Lanzhou urban area is located at a lower altitude, while Yongdeng County and Yuzhong County are located at higher altitudes. Furthermore, the urban area is located in the centre of Lanzhou, which tends to collect rainwater from the northern and southern areas, as can be seen in Figure 3; hence the Lanzhou urban area is prone to flooding during heavy rainstorm events. Since flooding always induces landslides accompanied with debris flows, the Lanzhou urban area is highly vulnerable to geohazards.

Topography
Topography has a critical influence on the development of geohazards. For example, the distribution of flooding hazards is significantly affected by variations in spatial topography [40]. Flat areas in particular are prone to flood. River networks play an important role in the reduction of floodwaters during flood events. The average altitude of China's northwestern region ranges from 2000 to 3000 m. The Yellow River, which starts on the Tibetan Plateau in Qinghai Province, passes through the provinces of Qinghai, Sichuan, Gansu, Ningxia, Inner Mongolia, Shanxi, and Henan. Figure 3 illustrates the altitude and distribution of the river system in Lanzhou. The altitude of Lanzhou ranges from 3796 to 1423 m with an average of approximately 2500 m. The Lanzhou urban area is located at a lower altitude, while Yongdeng County and Yuzhong County are located at higher altitudes. Furthermore, the urban area is located in the centre of Lanzhou, which tends to collect rainwater from the northern and southern areas, as can be seen in Figure 3; hence the Lanzhou urban area is prone to flooding during heavy rainstorm events. Since flooding always induces landslides accompanied with debris flows, the Lanzhou urban area is highly vulnerable to geohazards.

Loess Landforms
Longzhong Plain, known as the Chinese Loess Plateau, is covered by deep loess deposits at the northwestern side. These loess deposits are predominantly quaternary deposits, and can be divided into the Wuchen Loess, Lishi Loess, and Malan Loess deposits. Its loess landforms make the terrain of Lanzhou unique. Figure 4 shows the distribution of loess landforms across northwestern China [41]. As shown, these loess landforms are widely distributed, occurring in Lanzhou City, Gansu Province, Shanxi Province, Xining Province, and Yinchuan City. According to existing statistical data, quaternary loess covers 66% of the land surface of Lanzhou, with the thickness of loess

Loess Landforms
Longzhong Plain, known as the Chinese Loess Plateau, is covered by deep loess deposits at the northwestern side. These loess deposits are predominantly quaternary deposits, and can be divided into the Wuchen Loess, Lishi Loess, and Malan Loess deposits. Its loess landforms make the terrain of Lanzhou unique. Figure 4 shows the distribution of loess landforms across northwestern China [41]. As shown, these loess landforms are widely distributed, occurring in Lanzhou City, Gansu Province, Shanxi Province, Xining Province, and Yinchuan City. According to existing statistical data, quaternary loess covers 66% of the land surface of Lanzhou, with the thickness of loess deposits ranging from 100 to 150 m [27]. The loess landforms have many vertical joint fissures, which are always present in the loess plateau [42]. The characteristics of these sediments provide source material for the development of geohazards. deposits ranging from 100 to 150 m [27]. The loess landforms have many vertical joint fissures, which are always present in the loess plateau [42]. The characteristics of these sediments provide source material for the development of geohazards.

Geological Structures
Topographically, the Lanzhou urban area is a vulnerable place. Many geologic faults have developed in northwestern China, forming its special geographical environment. Figure 5 shows the distribution of faults around Lanzhou City. Because the Yellow River passes through the city, the area comprises four faulted-erosion basins. The Yellow River lies in the geometric centre of the rhomboidal western basin, which is demarcated by the Leitanhe Fault to the east, the Xinglong Fault to the south, the Shengou Fault to the west and the Jinchengguan Fault to the north [43]. Cretaceous and tertiary strata are deposited between Gaolan Mountain and Baita Mountain. Tectonic movements play an important role in the height of the uplifted strata, leading to a 2151 m height difference between Xinglong Mountain and the river valley region. Neotectonic movements are extremely strong in the Lanzhou urban area. Complex geological conditions and active faults around Lanzhou result in frequent earthquakes [44].

Geological Structures
Topographically, the Lanzhou urban area is a vulnerable place. Many geologic faults have developed in northwestern China, forming its special geographical environment. Figure 5 shows the distribution of faults around Lanzhou City. Because the Yellow River passes through the city, the area comprises four faulted-erosion basins. The Yellow River lies in the geometric centre of the rhomboidal western basin, which is demarcated by the Leitanhe Fault to the east, the Xinglong Fault to the south, the Shengou Fault to the west and the Jinchengguan Fault to the north [43]. Cretaceous and tertiary strata are deposited between Gaolan Mountain and Baita Mountain. Tectonic movements play an important role in the height of the uplifted strata, leading to a 2151 m height difference between Xinglong Mountain and the river valley region. Neotectonic movements are extremely strong in the Lanzhou urban area. Complex geological conditions and active faults around Lanzhou result in frequent earthquakes [44].
Sustainability 2018, 10, x FOR PEER REVIEW 5 of 21 deposits ranging from 100 to 150 m [27]. The loess landforms have many vertical joint fissures, which are always present in the loess plateau [42]. The characteristics of these sediments provide source material for the development of geohazards.

Geological Structures
Topographically, the Lanzhou urban area is a vulnerable place. Many geologic faults have developed in northwestern China, forming its special geographical environment. Figure 5 shows the distribution of faults around Lanzhou City. Because the Yellow River passes through the city, the area comprises four faulted-erosion basins. The Yellow River lies in the geometric centre of the rhomboidal western basin, which is demarcated by the Leitanhe Fault to the east, the Xinglong Fault to the south, the Shengou Fault to the west and the Jinchengguan Fault to the north [43]. Cretaceous and tertiary strata are deposited between Gaolan Mountain and Baita Mountain. Tectonic movements play an important role in the height of the uplifted strata, leading to a 2151 m height difference between Xinglong Mountain and the river valley region. Neotectonic movements are extremely strong in the Lanzhou urban area. Complex geological conditions and active faults around Lanzhou result in frequent earthquakes [44].

Hydrology
Lanzhou has a semi-arid climate and the loess is sensitive to climatic changes. Variations in altitude create temperature and climatic extremes. Figure 6 shows the hydrologic characteristics of Lanzhou City. As is shown, recorded data indicate significant temperature differences between the highest temperature (22.1 • C) and the lowest temperature (−6.1 • C). The rainfall in Lanzhou is relatively low, averaging 35.4 mm, and the heaviest rainfall is always accompanied by the greatest average temperature. Statistics show that almost 70% of rainfall is comprised of high-intensity and short-duration rainstorm events [26]. Continuous rainfalls of more than seven consecutive days occur from July to September [26].

Hydrology
Lanzhou has a semi-arid climate and the loess is sensitive to climatic changes. Variations in altitude create temperature and climatic extremes. Figure 6 shows the hydrologic characteristics of Lanzhou City. As is shown, recorded data indicate significant temperature differences between the highest temperature (22.1 °C) and the lowest temperature (−6.1 °C). The rainfall in Lanzhou is relatively low, averaging 35.4 mm, and the heaviest rainfall is always accompanied by the greatest average temperature. Statistics show that almost 70% of rainfall is comprised of high-intensity and short-duration rainstorm events [26]. Continuous rainfalls of more than seven consecutive days occur from July to September [26]. Rainfall plays an important role in the development of geohazards. Soils generally have high shear stress when unsaturated, but relatively low values when fully saturated [42,45]. Rainfall infiltration reduces soil shear strength, leading to the occurrence of geohazards [46,47]. The concentrated rainfall from July to September results in numerous geohazards during those months.

Methodology
Risk assessment is a quantitative or qualitative measure of risk related to a well-defined situation and an acceptable threat. Quantitative risk assessment needs to calculate two components of risk: the magnitude of potential loss and the probability of loss occurrence. An acceptable risk is one in which the cost or difficulty of implementing an effective countermeasure exceeds the expectation of loss.

Assessment Method
According to the literature [48,49], "disaster risk is a function of three components: hazard, exposure and vulnerability. Vulnerability is the reaction of the assets when exposed to the spatially variable forces produced by a hazard event. Vulnerability reveals the situation in future". This study does not consider the reaction after disaster; for our purposes, disaster risk is composed of hazard and exposure.
A hazard is an agent which may cause harm to a vulnerable target. Often the terms hazard and risk are used interchangeably. However, in the field of risk assessment, these two terms are different. A hazard is an agent that can induce harm or damage to humans, property, or the environment. Risk is the probability that exposure to a hazard will result in a negative consequence. Thus, a hazard poses no risk if there is no exposure to that hazard [49]. Therefore, in this study, Rainfall plays an important role in the development of geohazards. Soils generally have high shear stress when unsaturated, but relatively low values when fully saturated [42,45]. Rainfall infiltration reduces soil shear strength, leading to the occurrence of geohazards [46,47]. The concentrated rainfall from July to September results in numerous geohazards during those months.

Methodology
Risk assessment is a quantitative or qualitative measure of risk related to a well-defined situation and an acceptable threat. Quantitative risk assessment needs to calculate two components of risk: the magnitude of potential loss and the probability of loss occurrence. An acceptable risk is one in which the cost or difficulty of implementing an effective countermeasure exceeds the expectation of loss.

Assessment Method
According to the literature [48,49], "disaster risk is a function of three components: hazard, exposure and vulnerability. Vulnerability is the reaction of the assets when exposed to the spatially variable forces produced by a hazard event. Vulnerability reveals the situation in future". This study does not consider the reaction after disaster; for our purposes, disaster risk is composed of hazard and exposure.
A hazard is an agent which may cause harm to a vulnerable target. Often the terms hazard and risk are used interchangeably. However, in the field of risk assessment, these two terms are different. A hazard is an agent that can induce harm or damage to humans, property, or the environment. Risk is Sustainability 2018, 10, 304 7 of 21 the probability that exposure to a hazard will result in a negative consequence. Thus, a hazard poses no risk if there is no exposure to that hazard [49]. Therefore, in this study, geohazard risk is considered as the combination of hazard and exposure [50,51]. The disaster risk model used can be described by the following equation: where ⊗ is not a mathematical operation symbol. This equation is used only to express the definition of risk. Geohazard risk means that the development of disaster, which is associated with the interactions between hazard and exposure, is uncertain. Thus, the geohazard risk assessment model can be expressed by Equation (2): where R is risk and k is the index classification number. When k = 1, it is the hazard index; and when k = 2, it is the vulnerability index. S k refers to each factor in the assessment model, and can be calculated using Equation (3): where F k,i is the normalized value of index i of classification k, f k,i is the weight of index i of classification k, and Therefore, the assessment model can be redefined as Equation (4): where R is risk, H i , and E k are the hazard factor and exposure factor, and h i , and e k are the weight coefficients of each factor. To evaluate the geohazard risks of Lanzhou City, the analytic hierarchy process (AHP) was incorporated into the GIS. AHP is a useful method of evaluating both qualitative and quantitative risk [52][53][54]. The assessment procedure was implemented to build a risk assessment index system based on the factors affecting geohazards; to determine the weight of each assessment factor using AHP; to normalize each factor in the assessment system; to map the distribution of risk level of each assessment index; and to visualize the spatial distribution of geohazard risks in the regional scale using GIS.

AHP Assessment Structure
In this study, the risk of geohazards was evaluated by AHP in which the effects from both hazard and exposure were expressed by indices with different weights. Figure 7 shows the assessment structure for geohazard risk. This structure includes three parts: object layer, index layer, and factor layer. In the assessment structure, the object layer is geohazards risk (U). The index layer includes hazard index (U 1 ) and exposure index (U 2 ) the hazard index includes five factor layers-annual average rainfall (U 11 ), annual average rainy day (U 12 ), geological environment (U 13 ), earthquake magnitude (U 14 ), and historical disaster (U 15 )-while the exposure index includes four factor layers: topographical slope (U 21 ), topographical elevation (U 22 ), river density (U 23 ), and river proximity (U 24 ).

Weight Calculation by AHP
Based on the assessment structure, the judgment matrix of each assessment factor can be expressed by Equation (5): where Bp is the judgment matrix and bij is the element of judgment matrix, which represents the relative importance of factor bi to factor bj. When bi is significantly more important than bj, bij is set to 9 and bji is set to 1/9. Then, the weight coefficient of each factor can be calculated by Equation (6): where max ( )/( 1 ) CI n n λ = − − and λmax is the largest eigenvalue of the judgment matrix, which can be calculated based on Equation (8). RI is the average random consistency index.
If the value of CR is less than 0.1, then the matrix is reasonable, and should be reestablished. Using AHP incorporated with GIS technology, the level of geohazards risk can be evaluated and the spatial distribution of geohazards risk can be mapped.

Weight Calculation by AHP
Based on the assessment structure, the judgment matrix of each assessment factor can be expressed by Equation (5): where B p is the judgment matrix and b ij is the element of judgment matrix, which represents the relative importance of factor b i to factor b j . When b i is significantly more important than b j , b ij is set to 9 and b ji is set to 1/9. Then, the weight coefficient of each factor can be calculated by Equation (6): where M i = n n ∏ j=1 b ij ; b ij is the relative value of i factor to j factor, which ranges from 1 to 9 and their reciprocals. The consistency of judgment matrix can be validated by the value of consistency ratio (CR). It can be calculated by the following equation: where CI = (λ max − n)/(n − 1) and λ max is the largest eigenvalue of the judgment matrix, which can be calculated based on Equation (8). RI is the average random consistency index.
If the value of CR is less than 0.1, then the matrix is reasonable, and should be reestablished. Using AHP incorporated with GIS technology, the level of geohazards risk can be evaluated and the spatial distribution of geohazards risk can be mapped.

Normalization
To compare differences between each factor, all factors should be normalized before overlay analysis. After normalization, the value of each factor ranged from 0 to 1. For all assessment factors, risk will increase with greater topographic elevation and slope, and decrease with less river density. Thus, the assessment factors can be divided into two types: negative factors and positive factors. Positive factor and negative factor are normalized using Equations (9) and (10), respectively.
where ∨ i m and ∧ i m are the normalized values of the assessment factors, i m is the original value, i max is the maximum value of the factor, and i min is the minimum value. To facilitate the accuracy of the assessment results, all factors were based on a resolution of 30 m.

Extraction of Topography
Topography has critical effects in the development and redistribution of geohazards. Elevation and slope are widely used to reflect topographic characteristics. Natural hazards (for example, landslides accompanied with debris flows and floods) can cause considerably more catastrophic damage in flat regions than in steep regions. An area with a flat slope has a high risk of disaster [50,51]. Moreover, the river system also has an important role in inducing disasters, for example floods. The occurrence of flooding is related to the distribution of the drainage system. The river system not only can assist in draining rainwater but also can induce flooding [50,51]. Figure 8a-d show the area's topographic characteristics. For this paper, topographic elevation and slope were extracted from a digital elevation model (DEM) with a 30 m resolution in GIS environment. The original DEM data can be obtained from the geospatial data cloud (http://www.gscloud.cn/). River density refers to the length of the water channel per unit and was derived using a source radius of 1 km. River proximity refers to the distance to the closest river channels. The critical levels of proximity were set at 200 m, 400 m, 600 m, 800 m, and 1000 m, respectively.

Normalization
To compare differences between each factor, all factors should be normalized before overlay analysis. After normalization, the value of each factor ranged from 0 to 1. For all assessment factors, risk will increase with greater topographic elevation and slope, and decrease with less river density. Thus, the assessment factors can be divided into two types: negative factors and positive factors. Positive factor and negative factor are normalized using Equations (9) and (10), respectively.

Extraction of Topography
Topography has critical effects in the development and redistribution of geohazards. Elevation and slope are widely used to reflect topographic characteristics. Natural hazards (for example, landslides accompanied with debris flows and floods) can cause considerably more catastrophic damage in flat regions than in steep regions. An area with a flat slope has a high risk of disaster [50,51]. Moreover, the river system also has an important role in inducing disasters, for example floods. The occurrence of flooding is related to the distribution of the drainage system. The river system not only can assist in draining rainwater but also can induce flooding [50,51]. Figure 8a-d show the area's topographic characteristics. For this paper, topographic elevation and slope were extracted from a digital elevation model (DEM) with a 30 m resolution in GIS environment. The original DEM data can be obtained from the geospatial data cloud (http://www.gscloud.cn/). River density refers to the length of the water channel per unit and was derived using a source radius of 1 km. River proximity refers to the distance to the closest river channels. The critical levels of proximity were set at 200 m, 400 m, 600 m, 800 m, and 1000 m, respectively.

Investigation of Geohazards
Geohazards including debris flow, landslide, earthquake, and flooding are widely distributed and frequent in Lanzhou. Long-term historical geohazards are one of the most critical references for flood risk evaluation [29,30]. Figure 9a,b show the historical distribution of earthquakes and other geohazards in Lanzhou. Comparing Figure 9a to Figure 5 shows that the distribution of earthquake epicenters corresponds to the distribution of faults. As shown in Figure 9b, landslides are often accompanied by debris flows. The predominant geohazards in the region are debris flow; landslides accompanied by debris flow; intense rainstorms which induce flooding; earthquake; and flooding, which has significant effects on the development of debris flow and landslides.

Investigation of Geohazards
Geohazards including debris flow, landslide, earthquake, and flooding are widely distributed and frequent in Lanzhou. Long-term historical geohazards are one of the most critical references for flood risk evaluation [29,30]. Figure 9a

Investigation of Geohazards
Geohazards including debris flow, landslide, earthquake, and flooding are widely distributed and frequent in Lanzhou. Long-term historical geohazards are one of the most critical references for flood risk evaluation [29,30]. Figure 9a,b show the historical distribution of earthquakes and other geohazards in Lanzhou. Comparing Figure 9a to Figure 5 shows that the distribution of earthquake epicenters corresponds to the distribution of faults. As shown in Figure 9b, landslides are often accompanied by debris flows. The predominant geohazards in the region are debris flow; landslides accompanied by debris flow; intense rainstorms which induce flooding; earthquake; and flooding, which has significant effects on the development of debris flow and landslides. Figure 9. Distribution of (a) historical earthquakes (Data from Yuan [43]); (b) different geohazards collected from 1949 to 2000 (Data from Mei [26] and Chen [55]). Figure 9. Distribution of (a) historical earthquakes (Data from Yuan [43]); (b) different geohazards collected from 1949 to 2000 (Data from Mei [26] and Chen [55]).

Weights of Assessment Factors
As mentioned above, assessment of disaster risk using AHP incorporated into GIS is helpful in the assessment of geohazard risk [56][57][58][59]. In our assessment model, the hazard index includes annual average rainfall, annual average number of rainy days, geological conditions, earthquake magnitude, and distribution of historic disaster. The exposure index includes elevation, slope, river density, and river proximity. During the establishment of the judgment matrix, the influence of each factor on geohazard risk was investigated by a questionnaire, shown in the Appendix A. The questionnaire was given to more than 20 geohazard experts. Based on the questionnaire responses, the assessment matrixes of hazard index (B Hazard ) and exposure index (B Exposure ) were constructed as Equations (11) and (12). These two judgment matrixes are reasonable since the values of CR are less than 0.1.
By using Equation (6), the weight coefficient of each factor can be calculated. For example, the weight (F 1 ) for the hazard index can be calculated as follows: Similarly, the other weights can also be calculated. Table 1 lists the calculated weights of each assessment factor.

Assessment Results
To obtain the spatial distribution of geohazard risks, the typical factors were mapped and normalized as shown in Figure 10a-f. In these factors, elevation, slope, and river proximities are negative factors: with increased elevation and slope, the risk level will decrease. These negative factors are normalized using Equation (10). The positive factors can be normalized using Equation (9). We suppose that regions where past historical disasters are located have a higher risk than other areas. After normalization, the value of each factor ranged from 0 to 1 and reflected the effects of each disaster risk factor. Figure 10a-

Assessment Results
To obtain the spatial distribution of geohazard risks, the typical factors were mapped and normalized as shown in Figure 10a-f. In these factors, elevation, slope, and river proximities are negative factors: with increased elevation and slope, the risk level will decrease. These negative factors are normalized using Equation (10). The positive factors can be normalized using Equation (9). We suppose that regions where past historical disasters are located have a higher risk than other areas. After normalization, the value of each factor ranged from 0 to 1 and reflected the effects of each disaster risk factor. Figure 10a-e show the spatial distributions of normalized factors. The normalized factors adequately reflect the influence of geohazard risks. Using the normalized grid cell data and the weight coefficient for each factor, a spatial distribution of geohazard risk levels can be obtained. Figure 11 shows the spatial distribution of assessed geohazard risk levels in Lanzhou. As is shown, Lanzhou's urban area has a high risk of geohazards, particularly where there are lower elevations, greater river density, and a history of intensive disaster distribution. Gaolan County, which has a large river density, also has a high risk level. The results reveal that high geohazard risks always appear near rivers. Moreover, historically, geohazards occur at locations classified as having high risk levels. Therefore, long-term historical geohazards are also useful information in assessing risk levels.  Using the normalized grid cell data and the weight coefficient for each factor, a spatial distribution of geohazard risk levels can be obtained. Figure 11 shows the spatial distribution of assessed geohazard risk levels in Lanzhou. As is shown, Lanzhou's urban area has a high risk of geohazards, particularly where there are lower elevations, greater river density, and a history of intensive disaster distribution. Gaolan County, which has a large river density, also has a high risk level. The results reveal that high geohazard risks always appear near rivers. Moreover, historically, geohazards occur at locations classified as having high risk levels. Therefore, long-term historical geohazards are also useful information in assessing risk levels. Using the normalized grid cell data and the weight coefficient for each factor, a spatial distribution of geohazard risk levels can be obtained. Figure 11 shows the spatial distribution of assessed geohazard risk levels in Lanzhou. As is shown, Lanzhou's urban area has a high risk of geohazards, particularly where there are lower elevations, greater river density, and a history of intensive disaster distribution. Gaolan County, which has a large river density, also has a high risk level. The results reveal that high geohazard risks always appear near rivers. Moreover, historically, geohazards occur at locations classified as having high risk levels. Therefore, long-term historical geohazards are also useful information in assessing risk levels. Figure 11. Spatial distribution of geohazard risk levels in Lanzhou City. Figure 11. Spatial distribution of geohazard risk levels in Lanzhou City. Using the GIS system, the ratio of geohazard risk levels in different districts can be accounted. Table 2 tabulates the ratios of geohazard risk levels in different districts. As is shown, up to 32% of the total area is at high or very high risk. Over 36% of the total area is exposed to a medium level of disaster risk. Only 31% of the area is at low or very low risk. In terms of administrative division, the areas of relative high risk are mainly urban: more than 55% of the total high-risk area. The second greatest risk is to Gaolan Couty, with 44% at high or very high risk. Local governments can take measures to prevent geohazards in the regions at high risk based on these assessment results, which it is hoped will be helpful to local governments for sustainable development in urban planning.

Analysis
As is shown in Figure 9b, landslides are the most frequent geohazards in Lanzhou City. Thus, the following analyses will focus on landslides, discussing the relationships among topography, geological materials, and landslides.

Relationship between Topography and Landslides
The steep loess landforms that pass through the Yellow River result in a large number of geohazards, particularly loess landslides. Topographical characteristics often determine geohazard types. For example, Chongqing City is situated to the southeast of the Sichuan Basin, at the confluence of the Jialing River and Yangtze River, and has a similar topography within the Lanzhou region. Chongqing is surrounded by a dissected mountainous topography, exacerbated by deep weathering of Jurassic sandstone and mudstone. The city experiences at least 500 mostly rock landslides, 89 rock falls, and 185 unstable overhanging rocks per year [60,61].
Hong Kong, on the other hand, is a mountainous island with similar slope topography to Lanzhou. More than 65% of Hong Kong's land area has a slope gradient greater than 15 • , and 30% has a slope gradient greater than 30 • . Hong Kong, too, frequently suffers from catastrophic landslides [61]. The landslides in Hong Kong are the result of natural degradation of steep slopes. During intensive heavy rainstorms, landslides in soils and weathered rock mostly occur on steep slopes, and are generally characterized by rapid motions and fluid-like flows. Investigations have shown that the risk of landslides on the island is low where slopes are gentle and high where they are steep [62][63][64]. Landslides are prone to occur in steep topography, which may lead to severe landslide damages [55,65]. Figure 12 depicts the relationship between the frequency of landslides and slope steepness in Lanzhou City: about 40% of landslides occur in regions with slopes between 35 • and 50 • .

Relationship between Geological Materials and Landslides
Material characteristics and hydraulic factors both contribute to the initiation of geohazards. Loess soil is subject to collapse with rainfall infiltration. Soil collapsibility plays a critical role in slope stability [55,65]. Slopes with loess soils lose significant stability when encountering concentrated rainfall, which eventually leads to landslides. Figure 13 shows the number of landslides in Lanzhou, by type with their attributed triggers, from 1990 to 2009 [55]. As is shown, most landslides in Lanzhou are loess landslides triggered by rainfall. Moreover, quaternary loese deposits at steep topographies develop debris flow during rainstorms. Debris flow can be initiated from landslides and saturated slope failure, causing rapid sediment entrainment. Geological materials have significant effects on landslide velocity [66][67][68]. Loess and erosion settlements easily create the homogeneous fluids of debris flow, which can travel with high flow velocity, impact force, and long runout in Lanzhou City.

Relationship between Geological Materials and Landslides
Material characteristics and hydraulic factors both contribute to the initiation of geohazards. Loess soil is subject to collapse with rainfall infiltration. Soil collapsibility plays a critical role in slope stability [55,65]. Slopes with loess soils lose significant stability when encountering concentrated rainfall, which eventually leads to landslides. Figure 13 shows the number of landslides in Lanzhou, by type with their attributed triggers, from 1990 to 2009 [55]. As is shown, most landslides in Lanzhou are loess landslides triggered by rainfall. Moreover, quaternary loese deposits at steep topographies develop debris flow during rainstorms. Debris flow can be initiated from landslides and saturated slope failure, causing rapid sediment entrainment. Geological materials have significant effects on landslide velocity [66][67][68]. Loess and erosion settlements easily create the homogeneous fluids of debris flow, which can travel with high flow velocity, impact force, and long runout in Lanzhou City.

Relationship between Geological Materials and Landslides
Material characteristics and hydraulic factors both contribute to the initiation of geohazards. Loess soil is subject to collapse with rainfall infiltration. Soil collapsibility plays a critical role in slope stability [55,65]. Slopes with loess soils lose significant stability when encountering concentrated rainfall, which eventually leads to landslides. Figure 13 shows the number of landslides in Lanzhou, by type with their attributed triggers, from 1990 to 2009 [55]. As is shown, most landslides in Lanzhou are loess landslides triggered by rainfall. Moreover, quaternary loese deposits at steep topographies develop debris flow during rainstorms. Debris flow can be initiated from landslides and saturated slope failure, causing rapid sediment entrainment. Geological materials have significant effects on landslide velocity [66][67][68]. Loess and erosion settlements easily create the homogeneous fluids of debris flow, which can travel with high flow velocity, impact force, and long runout in Lanzhou City.

Strategic Decision Making Suggestions
Strategic decision making regarding geohazard mitigation should be conducted based on the comprehensive technical assessment of geohazard risk. The geographical position and environmental conditions of Lanzhou require specific strategies to prevent and mitigate geohazards. Strategies adopted by local residents and government should follow the procedures proposed by the United Nations Office for Disaster Risk Reduction (UNISDR) [69]: governments should organize professionals and researchers to develop a comprehensive unified plan; communities should improve popular awareness of environmental protection and prohibit human activities that may trigger geological disasters; and local governments or communities should carry out comprehensive surveys and establish instructive prevention plans against geohazards.

Countermeasures for Landslide-Accompanied Debris Flow
Landslides accompanied by debris flows have a significant relationship to rainfall events [70,71]. Along with the aforementioned strategic decision making, it is also important to adopt technical countermeasures to prevent landslides. To avoid deterioration, regular slope maintenance is essential. Effective strategies to guarantee slope stability against landslides are: revising slopes to reduce the inclination of the slope from steep to gentle; constructing retaining walls (for example, reinforced earth or concrete walls); providing anchors (for example, soils and pre-stressed anchors) to reduce slop failures; providing effective drainage systems; and compacting loose fill slopes [72,73]. The rest of this section discusses detailed countermeasures to mitigate landslides accompanied debris flows, since these geohazards are common in Lanzhou City.
Flexible barriers are common countermeasures to restrain debris flow and protect from rock fall hazards [74][75][76]. These barriers are made of a light structure in which a certain number of posts keep in position a steel mesh and a series of steel cables. The mesh stops the moving blocks and transfers all the forces to the anchors in the ground. Figure 14a,b show this type of flexible barrier: Figure 14a shows protection against rock hazard and Figure 14b shows the application of flexible barriers for debris flows in Zhouqu City, Gansu Province, which suffered a strong debris flow in 2014 [76]. As mentioned before, Lanzhou city is surrounded by several mountains (Gaolan Mountain, Xinlong Mountain, and Baita Mountain: see Figure 5). Structural countermeasures such as rigid and flexible barriers are commonly installed in mountainous regions to intercept mass-wasting processes, and are recommended to prevent debris flows in Lanzhou. In addition, to mitigate the damage of debris flow, local governments that monitor debris flow warning systems should adhere to operating standards in the monitoring process.

Strategic Decision Making Suggestions
Strategic decision making regarding geohazard mitigation should be conducted based on the comprehensive technical assessment of geohazard risk. The geographical position and environmental conditions of Lanzhou require specific strategies to prevent and mitigate geohazards. Strategies adopted by local residents and government should follow the procedures proposed by the United Nations Office for Disaster Risk Reduction (UNISDR) [69]: governments should organize professionals and researchers to develop a comprehensive unified plan; communities should improve popular awareness of environmental protection and prohibit human activities that may trigger geological disasters; and local governments or communities should carry out comprehensive surveys and establish instructive prevention plans against geohazards.

Countermeasures for Landslide-Accompanied Debris Flow
Landslides accompanied by debris flows have a significant relationship to rainfall events [70,71]. Along with the aforementioned strategic decision making, it is also important to adopt technical countermeasures to prevent landslides. To avoid deterioration, regular slope maintenance is essential. Effective strategies to guarantee slope stability against landslides are: revising slopes to reduce the inclination of the slope from steep to gentle; constructing retaining walls (for example, reinforced earth or concrete walls); providing anchors (for example, soils and pre-stressed anchors) to reduce slop failures; providing effective drainage systems; and compacting loose fill slopes [72,73]. The rest of this section discusses detailed countermeasures to mitigate landslides accompanied debris flows, since these geohazards are common in Lanzhou City.
Flexible barriers are common countermeasures to restrain debris flow and protect from rock fall hazards [74][75][76]. These barriers are made of a light structure in which a certain number of posts keep in position a steel mesh and a series of steel cables. The mesh stops the moving blocks and transfers all the forces to the anchors in the ground. Figure 14a,b show this type of flexible barrier: Figure 14a shows protection against rock hazard and Figure 14b shows the application of flexible barriers for debris flows in Zhouqu City, Gansu Province, which suffered a strong debris flow in 2014 [76]. As mentioned before, Lanzhou city is surrounded by several mountains (Gaolan Mountain, Xinlong Mountain, and Baita Mountain: see Figure 5). Structural countermeasures such as rigid and flexible barriers are commonly installed in mountainous regions to intercept mass-wasting processes, and are recommended to prevent debris flows in Lanzhou. In addition, to mitigate the damage of debris flow, local governments that monitor debris flow warning systems should adhere to operating standards in the monitoring process.

Early Monitoring System for Earthquake and Flood
Earthquakes can cause frequent rock falls and landslides in mountain ridges and watershed slopes. Landslide masses are scattered on the upper middle parts of the mountains in Lanzhou. When subjected to intense rainfall, loose deposits accumulate along gullies, forming devastating high-velocity debris flows: a typical earthquake-landslide-debris flow disaster chain [77]. Earthquakes and floods are the inducing factors for landslides and debris flows. Therefore, to mitigate geohazard ripple effects, it is important to monitor earthquakes and rainstorms. The China Earthquake Administration has authorized the China Earthquake Background Field Exploration Project and the National Earthquake Social Service Project, and constructed two demonstrative earthquake early warning systems in Beijing and Lanzhou [44]. The construction and operation of these early warning systems have provided excellent results for the development of earthquake early warning in China.
Since Lanzhou's rainfall is concentrated from July to September, potential flooding disaster may develop in these months. Due to lower rainfall in Lanzhou City, flood risk there is relatively lower than elsewhere in the province, but rainfall can induce other geohazards (for example, landslides accompanied with debris flows). Therefore, it is also necessary to monitor the rainfall in Lanzhou.

Conclusions
This work has investigated geohazards affecting Lanzhou. Local landforms, topography, geology, and hydrogeology were analyzed and risk assessment was applied to evaluate geohazard risk levels in different regions. Based on our investigations and risk assessment, countermeasures for prevention of geohazards have been discussed. The flowing conclusions can be drawn: 1.
The geohazards investigation shows that intensive seasonal short-duration rainstorm and special loess landforms with many vertical fissures contribute to frequent landslides, which are often accompanied by debris flows. Faults and strong neotectonic movements provide geological conditions which are conducive to earthquakes. The disaster chain of earthquake-landslide-debris flow is common in Lanzhou City.

2.
Geohazard risk assessment results show that 32% of the total area is at high or very high risk. About 55% of the urban area and 44% of Gaolan county are at high or very high risk. During the monitoring process, local governments should pay particular attention to these high-risk regions.

3.
Appropriate countermeasures include the establishment of relationships between geohazards (for example, landslides accompanied with debris flows), topography, and geological materials. Regular slope maintenance is an effective way to further mitigate landslide events, while the installation of flexible barriers is useful to restrain debris flows.