# Dynamic Vulnerability Analysis of Mountain Settlements Exposed to Geological Hazards: A Case Study of the Upper Min River, China

Authors: Mingtao Ding; Chuan Tang; Tao Huang; Zemin Gao
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

The upper reaches of Min River (The upper Min River) is located in Southwest China with significant mountain settlements, which are vulnerable to frequent geological hazards. Based on a field investigation, collation of yearbook data, and analysis through the use of SPSS statistical software, a vulnerability evaluation index system of geological hazards was devised. According to the actual field situation and the acquired data of the study area in 2006, 2009, and 2015, 16 indicators were selected as settlement vulnerability evaluation indexes of geological hazards. The indexes included population density, building coverage, and economic density. Based on the comprehensive evaluation model of entropy value, the dynamic change in the settlement vulnerability of geological hazards was analyzed. The results showed that population density, building coverage, economic density, and road density were the factors that affected the settlement vulnerability of geological hazards the most—Wenchuan earthquake caused considerable damage to the upper Min River, making the area the most vulnerable in 2009. However, its vulnerability decreased in 2015, which indicated that postearthquake reconstruction achieved significant results. Thus, the vulnerability has emerged as an important indicator reflecting the safety and healthy development of mountain settlements.


---

## Body

## 1. Introduction

The vulnerability assessment of geological hazards is a necessary process of regional risk analysis. Furthermore, the vulnerability of geological hazards fluctuates with time. The upper Min River, as a topographic uplift and economic valley area, is in the process of rapid construction of social structure after the settlement was devastated by the massive Wenchuan earthquake of 2008 [1]. Therefore, the dynamic assessment of mountain settlement vulnerability of geological hazards in the upper Min River is an essential part of the mechanism for regional disaster risk reduction, prevention, and capacity building for emergencies.In 1848, John Stuart Mill pointed out in his book, thePrinciples of political economy with some of their applications to social philosophy, that all signs of destruction caused by earthquakes, floods, hurricanes, and wars would disappear in a short time and that the country would recover rapidly from the state of disasters, indicating that human and social factors were important aspects affecting disaster risk [2, 3]. As a result, it was gradually realized that improving social disaster-bearing and emergency response capacity was more accessible than changing disaster-causing factors which further helped people understand the importance of vulnerability to disaster research. For natural disasters, the important thing is not to explain the characteristics, consequences, and causes of disasters by the conditions or the actions of disasters but to analyze the contemporary social order, assess the daily relations of disasters, and shape the historical environment of these characteristics [4]. In addition, socioeconomic factors, such as population characteristics, industrial structure, and spatial distribution, have an impact on natural disasters. Therefore, vulnerability is not only a central concept in the field of disaster research but also a means to solve the problems related to population, development, and environment [5, 6].At present, vulnerability research has developed rapidly—not only as qualitative research but also as quantitative research. Guo [7] analyzed the relationship between natural disasters and social vulnerability, believing that vulnerability includes material and spiritual aspects and that there is a close relationship between social vulnerability and disaster types, disaster nature, disaster intensity, and disaster occurrence time. Hu et al. [8] analyzed the impact of social factors on disasters, pointing out that social factors mainly affect disasters by influencing the source of disasters, the ability of disaster-bearing bodies to prevent and resist disasters, and the level of social and economic development. Liu et al. [9, 10] divided vulnerability evaluation factors into property index and population index. He put forward the conversion assignment function, which solved the problem of unified scaling and comprehensive expression of human, financial, and material aspects and discussed the vulnerability calculation method. Pang et al. [11] analyzed five parts of dynamic risk assessment of natural disasters: time constraints, disaster-causing factors, disaster-bearing bodies, coupling, and scenario output. They also gave a basic model of dynamic risk assessment of natural disasters in a narrow sense.Meanwhile, many scholars have made vulnerability analyses based on a large number of cases. Taking Kangding County as an example, Wang et al. [12] put forward a new method for evaluating the socioeconomic vulnerability of urban geological hazards and established a vulnerability evaluation system. Tang [13] selected population density, housing property, GDP, and cultivated land and highway distribution as vulnerability evaluation factors of landslide disasters in compiling the landslide risk map of the Red River Basin. Wu et al. [14] selected housing construction, road engineering, lifeline engineering, population distribution, and land use as vulnerability assessment factors in landslide hazard risk assessment in the new urban area of Badong County. Miao and Ding [15] selected eight factors as evaluation indexes. The factors were population density, building coverage, road density, number of hospitals per 10,000 persons, number of welfare institutes per 10,000 persons, number of village committees per 10,000 persons, mobile users' proportion at the end of the year, and number of schools per 10,000 persons. They used the entropy method to obtain the vulnerability zoning map of the Lushan earthquake area.Vulnerability research has become an indispensable part of natural disaster risk research [16]. However, the overall research still relates mainly to the analysis of the natural factors causing disasters, and there is still a lack of correct understanding and in-depth discussion of the social factors involved in disasters.This study chooses the mountain settlements in the upper Min River as the research object. It uses the entropy method to carry out the dynamic assessment of the vulnerability of mountain settlements to geological hazards under time constraints and obtains the zoning maps of the vulnerability in different periods. Thus, the study provides an essential scientific basis for the compilation of data on regional geological hazards, planning risk prevention, and the formulation of emergency measures.


## 2. Study Area

The upper Min River is situated between 30°45′N to 33°10′N and 102°35′E to 103°57′E. This area belongs to the uplift area of China’s topography and the active belt of the tectonic movement. The occurrence of several strong earthquakes in the region has resulted in the loosening of the surface material, damage to deep rocks and the soil mass, and the decrease in water saturation rate in the study area [17, 18]. All these factors have led to frequent secondary disasters in the area, such as topples, slides, and flows. The paper mainly studies the dynamic vulnerability analysis of mountain settlements and does not involve unmanned areas. Therefore, it only selects the geological disaster points that bring harm to human survival or damage the human living environment. The geoenvironmental information system of Sichuan Province [19] stores such data and does not involve geological phenomena such as slides and flows in unmanned areas. According to the geoenvironmental information system of Sichuan province, there have been 2,778 geological hazards in the study area. According to the new classification system proposed by Hunger et al. [20], these geological hazards include 803 slides, 726 topples, 770 flows, and 479 slopes deformation (Figure 1). The study area is inhabited by 22 ethnic groups, such as Tibet, Qiang, Hui, and Han.Figure 1
Location of the upper Min River (inset map shows Sichuan and China) and the distribution of geological disasters in the study area.Moreover, the area is home to the largest Qiang ethnic community in China. Its ethnic culture has distinctive characteristics of mixed cohabitation. Besides, under the influence of the farming culture of the Central Plains and the traditional Qiang-Tibetan culture, the region’s economy is still relatively backward [21].


### 2.1. Distribution and Evolution of Mountain Settlements

The spatial pattern of settlements in the upper Min River has changed significantly from 2006 to 2015. Its distribution is manifested mainly by the migration of the settlements with time from high altitude areas and large slope areas to the river valley platforms. The trend indicated a shift from decentralization to agglomeration. Consequently, the settlement area gradually decreased from 955 km2 in 2006 to 462 km2 in 2015 (Figure 2).Figure 2
Distribution of mountain settlements in the upper Min River for different years. (a) 2005. (b) 2009. (c) 2015.


## 2.1. Distribution and Evolution of Mountain Settlements

The spatial pattern of settlements in the upper Min River has changed significantly from 2006 to 2015. Its distribution is manifested mainly by the migration of the settlements with time from high altitude areas and large slope areas to the river valley platforms. The trend indicated a shift from decentralization to agglomeration. Consequently, the settlement area gradually decreased from 955 km2 in 2006 to 462 km2 in 2015 (Figure 2).Figure 2
Distribution of mountain settlements in the upper Min River for different years. (a) 2005. (b) 2009. (c) 2015.


## 3. Materials and Methods

### 3.1. Division of Cell Meshes

The primary data of this evaluation are based on the township data, and the operation and processing of the data are determined from the division of regional grid cells. This is done due to the apparent advantage of using grid cells in spatial data superposition calculation. The division of grid size directly affects the rationality of evaluation results. At present, the selection of grid size depends mainly on the resolution of original data and the experience and knowledge of the experts.“Technical Guide for Evaluation of Resources and Environment Carrying Capacity and Suitability of Land and Space Development” (for trial implementation) [22] proposes to use a grid with a spatial resolution of 50 × 50 m at the regional scale. Therefore, the vulnerability assessment unit of this study adopts a 50 × 50 m resolution grid unit.


### 3.2. Construction of the Evaluation Index System

The index system of vulnerability assessment includes population, economy, society, and other factors, which referred to the existing research results on disaster vulnerability and other related index systems in mountain areas [23–26]; 43 evaluation indicators were selected from the aspects of population characteristics, social security, social economy, structural characteristics, and natural characteristics. The evaluation index system for the settlement vulnerability of geological hazards in the upper Min River was devised at four levels: target level, criterion level, element layer, and index level (Table 1).Table 1
Vulnerability evaluation index system of the upper Min River.
Target layerCriteria layerElement layerIndex layerAttributeEvaluation index of geological disaster vulnerability of mountain settlementsSocial populationDisadvantagedFemale ratioPositiveNumber of the disabledPositiveNumber of minimum living security population for urban residentsPositiveThe proportion of agricultural and animal husbandry populationPositivePopulation ratio under 4 years old and over 65 years oldPositivePopulation pressurePopulation densityPositiveBuilding coveragePositiveVulnerable occupationNumber of people in the transportation industryPositiveNumber of mining workersPositiveConstruction practitionersPositiveSocioeconomicEconomic developmentEconomic densityPositivePer capita GDPNegativeLocal fiscal revenueNegativeSavings balance of urban and rural residentsNegativeGDP ratio of the primary industryPositivePer capita income of farmersNegativeTertiary industry incomeNegativeSocial capitalInvestment attraction assetsNegativePer capita investment in fixed assetsNegativeGeneral public budget expenditureNegativeSocial securitySocial securityRoad densityNegativeNumber of beds in social welfare institutionsNegativeThe proportion of social security expenditureNegativeNumber of village committees per 10,000 personsNegativeNumber of health institutionsNegativeNumber of medical technicians per 10,000 personsNegativeNumber of beds per 10,000 personsNegativeSocial cultureSocial civilizationNumber of people with high school background per 10,000 personsNegativeTeacher-student ratioNegativeRadio and television coverageNegativeMobile base station coverage rateNegativePer capita mobile phonesNegativeIlliteracy ratioPositiveCustoms and cultureThe ethnic minority population ratioPositiveEthnic culture categoriesPositiveHazard indexDanger of hazardHazard densityPositiveNumber of people threatened by geological disastersPositiveProperty threatened by geological disastersPositiveProtective engineeringProtection engineeringNegativeInvestment in geological hazard controlNegativeSurface protectionSoil and water loss intensityPositiveVegetation coverageNegativeTopographic undulationPositiveThe bold text represents the final selection indexes.The vulnerability evaluation indexes of geological hazards are different at different scales. The social and economic development of the upper Min River is reduced, and data collection is complicated. In the selection of indicators, we must not only follow the scientific principle but also fully consider the reality to ensure the availability of data. 27 indexes were excluded because they could not collect the complete time series data. 16 indexes are selected as vulnerability evaluation factors of the geological disasters in the upper Min River. These indexes include property threatened by geological disasters (D1), number of people threatened by geological disasters (D2), population density (D3), hazard density (D4), teacher-student ratio (D5), economic density (D6), GDP ratio of primary industry (D7), building coverage (D8), road density (D9), number of village committees per 10,000 persons (D10), per capita GDP (D11), per capital investment in fixed assets (D12), general public budget expenditure (D13), per capita mobile phones (D14), number of medical technicians per 10,000 persons (D15), and engineering protection (D16).The basic data of this study are obtained mainly from the following three aspects:(1)
Database on the potential points of geological hazards: According to the geoenvironmental information system of Sichuan province [19] (http://202.61.89.33:16003), information related to coordinates, geographic location, direction, volume, scale, property threatened by geological disasters, number of people threatened by geological disasters, number of households threatened by geological disasters, and prevention engineering of the potential geological hazards in the upper Min River since 2008 was obtained.(2)
Remote sensing images and vector data: Based on the National Earth System, science data-sharing infrastructure [27] (http://www.geodata.cn), the enhanced thematic mapper (ETM) image of the study area, and the essential data related to fields like topography, land use, and earthquakes were obtained.(3)
Data from yearbooks and statistical bulletins: The data on administrative areas, economic output, a permanent population, practitioners of secondary and tertiary industries, industrial output, built-up areas of towns and townships, and other related data were obtained from several yearbooks and statistical bulletins for the years 2006 to 2016. The sources included “China County Economic Statistical Yearbook,” “Sichuan Statistical Yearbook,” “Aba Statistical Yearbook,” “Chengdu Statistical Yearbook,” “Heishui Statistical Yearbook,” “Wenchuan Statistical Yearbook,” “Mao Statistical Yearbook,” “Songpan Statistical Yearbook,” and “Li Statistical Yearbook.” The data not included in yearbooks and statistical bulletins were obtained through field investigation. For the readjustment of the administrative division of Dujiangyan City, at the end of 2014, Zipingpu Town and Hongkou Town were merged into Longchi Town. However, since statistical data are still divided into the yearbooks for three townships, this article does not make any adjustment to it.The indices selected for vulnerability assessment were defined as follows:(1)
Property threatened by geological disasters (D1): geological disaster threatening property refers to the amount of property that may be lost after a geological disaster occurs. The more threatening property, the more property that may be damaged, and the higher the vulnerability, which is a positive index.(2)
The number of people threatened by geological disasters (D2): the number of people threatened by geological disasters is the number of people threatened by geological disasters in each township. The more people threatened by geological disasters, the more people may be damaged by geological disasters, and the higher the vulnerability, which is a positive index.(3)
Population density (D3): population density is the number of people per unit area. Vulnerability is first related to population density. The higher the population density, the higher the vulnerability of settlements, which is a positive index.(4)
Hazard density (D4): the more geological disaster points, the greater the potential threat to residents and the greater the vulnerability, which is a positive index.(5)
Teacher-student ratio (D5): it is the ratio of students to teachers. The more students assigned by each teacher, the weaker the ability to respond to the crisis. The role of the teacher is very important when disaster strikes. The more teachers are relative to the students in the evacuation and transfer project, the more conducive to ensuring the safety of students’ lives, the higher the antirisk ability, and the lower the vulnerability, which is a negative index.(6)
Economic density (D6): economic density refers to the ratio of GDP to the regional area, which characterizes the efficiency of economic activities per unit area of settlements and the intensity of land use. The greater the economic density, the higher the probability of loss in the event of a geological disaster, and the higher the vulnerability of the settlement, which is a positive index.(7)
GDP ratio of primary industry (D7): the GDP ratio of primary industry in the upper Min River mainly refers to agriculture, forestry, animal husbandry, and fisheries. The GDP ratio of the primary industry is an industry with more inputs and less income. Such industries were the leading industries in the upper Min River before 2006, reflecting the society. The higher the degree of closure of the structure, the higher the output value ratio, the higher the vulnerability, which is a positive index.(8)
Building coverage (D8): buildings are also sensitive to geological disasters and are the main disaster-bearing bodies of geological disasters in the upper reaches of the Minjiang River. The higher the building coverage, the higher the vulnerability of settlements, which is a positive index.(9)
Road density (D9): road density refers to the ratio of the length of the road in the area to the area of the area, and it reflects the density of roads. The greater the road density, the higher the probability of damage (89-90), and the higher the vulnerability, which is a positive index.(10)
The number of village committees per 10,000 persons (D10): this is the number of village committees owned by every 10,000 people in the study area. Village committees are grassroots mass organizations for self-service, self-management, and self-education of villagers. They are stable and developed in a small area. The primary organization is the grassroots force for the transmission and implementation of national disaster prevention and reduction policies. The more village committees there are, the better the national disaster prevention and mitigation policies are communicated and implemented, and the lower the vulnerability of settlements, which is a negative index.(11)
Per capita GDP (D11): per capita GDP is mainly used to assess the economic development level of a region. At present, there are differences between urban and rural areas in my country. In terms of economic conditions, the financial status of rural residents is worse than that of urban residents. And the more barren a place is, the lower the danger awareness of geological disasters, and the higher the degree of damage when a disaster occurs. And the more personal wealth, the stronger the ability to reconstruct after the disaster, which can reduce the vulnerability of the system to a certain extent. Therefore, the higher the per capita GDP of a region, the lower the vulnerability, which is a negative index.(12)
Per capita investment in fixed assets (D12): per capita investment in fixed assets is the economic activity of construction and purchase of fixed assets, that is, the reproduction of fixed assets. Investment in fixed assets is the primary means of social fixed asset reproduction. Per capita, fixed assets reflect the strength and ability of disaster reconstruction. The greater the per capita investment in fixed assets, the lower the vulnerability, which is a negative index.(13)
General public budget expenditure (D13): it is the expenditure that the region allocates and uses in a planned manner to the centralized budget revenue, and it is the embodiment of the local grassroots construction capacity. The higher the general public budget expenditure, the better the implementation of disaster prevention and mitigation policies, and the lower the vulnerability. It is a negative index.(14)
Per capita mobile phones (D14): the number of mobile phones is an essential index of regional communication capabilities. The more the number of telephones per capita, the faster the transmission of disaster information, which is more conducive to benefit and avoid harm, and the lower the vulnerability. It is a negative index.(15)
The number of medical technicians per 10,000 persons (D15): it indicates that there are medical technicians for every 10,000 people in the study area, and many places in the mountain regions cannot build a medical system. The number of medical technicians can reflect the local rescue ability. The greater the number of doctors, the better the rescue ability after the disaster and the lower the vulnerability of the settlement, which is a negative index.(16)
Engineering protection (D16): protection engineering refers to the number of construction protection projects to control geological disasters. Engineering protection is an effective measure to protect the safety of local residents' lives and property. Therefore, the greater the number of engineering protection, the lower the vulnerability, which is a negative index.


### 3.3. Evaluation Method

Entropy is an objective, comprehensive evaluation method whose concept originated from thermodynamics [28–31]. It is a measure of uncertainty—the smaller the information entropy of the indicators, the higher the difference coefficient and the weight of the indicators. This means that the information provided by the indicators is large. Thus, it plays an essential role in the comprehensive application. The entropy method is increasingly applied in engineering, management, and social sciences and has gradually become an indispensable practical method in scientific research. Miao and Ding [30] took Lushan and Ludian earthquake-stricken areas as research areas and used the entropy comprehensive evaluation method to determine the social vulnerability of the affected areas. This further proved that the entropy method is practical in social vulnerability assessment. According to the theory and method of information entropy, Ding et al. [32] put forward the concept of “geological hazard entropy” and applied it to the evaluation of geological hazards in the Ranwu-Dongjiu section of Sichuan-Tibet Highway and achieved excellent results.This study combines the entropy method with geographic information system (GIS) to objectively and intuitively analyze the settlement vulnerability of geological hazards in the upper Min River. The main steps of the entropy comprehensive evaluation method are as follows [31]:(a)
Normalization of indexes [31]:Due to the differences in the dimensions, orders of magnitude, units, and quantity changes of the selected indicators in the study, these differences may have an impact on the final analysis results. In order to eliminate these effects, it is necessary to uniformly eliminate the dimensions of the data to solve the problem of the homogenization of qualitative index values. Moreover, because the values of positive and negative indicators represent different meanings (the greater the value of the positive indicator, the greater the vulnerability of mountain settlements, and the smaller the value of the negative indicator, the greater the vulnerability of mountain settlements) Therefore, we use different algorithms to standardize data for positive and negative indicators.The mountain settlements in the upper Min River are divided inton grids, and p indicators are selected. xij is the value of the jth indicator of the ith grid unit (i = 1, 2, …, n; j = 1,2, …, p). The matrix is as follows:(1)X=x11x12⋯x1px21x22⋯x2p⋮⋮⋮xn1xn2⋯xnp=X1,X2,…,Xp.When Xij is a positive indicator [31],(2)Xij∗=Xij−minX1j,X2j,L,XnjmaxX1j,X2j,L,Xnj−minX1j,X2j,L,Xnj,i=1,2,L,n;j=1,2,L,p.When Xij is a negative indicator [31],(3)Xij∗=maxX1j,X2j,L,Xnj−XijmaxX1j,X2j,L,Xnj−minX1j,X2j,L,Xnj,i=1,2,n;Lj=1,2,Lp,wherei represents the sample; j represents the indicator, Xij′ represents the standardized value of the ith sample under the j indicator, minX1j,X2j,L,Xnj represents the minimum value in the j indicator value, and maxX1j,X2j,L,Xnj represents the maximum value in the j indicator value.Normalized data after calculation:(4)X∗=x11∗x12∗⋯x1p∗x21∗x22∗⋯x2p∗⋮⋮⋮xn1∗xn2∗⋯xnp∗=X1∗,X2∗,⋯,Xp∗.(b)
Computation of entropy value [31]:(5)ej=−k∑i=1npijlnpij,whereej represents the entropy value of index j; k > 0; pij=xij∗/∑ni=1xij∗ represents the weight of scheme i under index j; n represents the number of samples.(c)
Computation of the different coefficient [31]:(6)gj=1−ej,wheregj represents the different coefficient of index j and ej is the same as in (5) above.(d)
Computation of the index weight [31]:(8)aj=gj∑j=1ngj,where aj represents the weight of index j and gj is the same as in (6) above.Through ArcGIS software, establish the raster layer of each index, use the above entropy value calculation formula, and use the raster calculator to calculate the vulnerability to obtain the vulnerability value of the mountain settlement. At the same time, in order to conduct a comparative analysis of the vulnerability of mountain settlements in the upper Min River in different periods, the results of mountain settlement evaluations in 2006, 2009, and 2005 were subjected to dispersion standardization (Min-Max standardization). So that the vulnerability index of each period is changed to between [0,1]. The vulnerability level is divided into 5 categories using the equal division method, which are in turn, low (0.0–0.2), comparatively low (0.2–0.4), medium (0.4–0.6), comparatively high (0.6–0.8), and high (0.8–1.0) five mountain settlement vulnerable areas. The standardized transfer function for dispersion is as follows [31]:(9)X′=X−minmax−min,where X′ is the normalized value of the sample data, X is the original value of the sample data, max is the maximum value of the sample data, and min is the minimum value of the sample data.


## 3.1. Division of Cell Meshes

The primary data of this evaluation are based on the township data, and the operation and processing of the data are determined from the division of regional grid cells. This is done due to the apparent advantage of using grid cells in spatial data superposition calculation. The division of grid size directly affects the rationality of evaluation results. At present, the selection of grid size depends mainly on the resolution of original data and the experience and knowledge of the experts.“Technical Guide for Evaluation of Resources and Environment Carrying Capacity and Suitability of Land and Space Development” (for trial implementation) [22] proposes to use a grid with a spatial resolution of 50 × 50 m at the regional scale. Therefore, the vulnerability assessment unit of this study adopts a 50 × 50 m resolution grid unit.


## 3.2. Construction of the Evaluation Index System

The index system of vulnerability assessment includes population, economy, society, and other factors, which referred to the existing research results on disaster vulnerability and other related index systems in mountain areas [23–26]; 43 evaluation indicators were selected from the aspects of population characteristics, social security, social economy, structural characteristics, and natural characteristics. The evaluation index system for the settlement vulnerability of geological hazards in the upper Min River was devised at four levels: target level, criterion level, element layer, and index level (Table 1).Table 1
Vulnerability evaluation index system of the upper Min River.
Target layerCriteria layerElement layerIndex layerAttributeEvaluation index of geological disaster vulnerability of mountain settlementsSocial populationDisadvantagedFemale ratioPositiveNumber of the disabledPositiveNumber of minimum living security population for urban residentsPositiveThe proportion of agricultural and animal husbandry populationPositivePopulation ratio under 4 years old and over 65 years oldPositivePopulation pressurePopulation densityPositiveBuilding coveragePositiveVulnerable occupationNumber of people in the transportation industryPositiveNumber of mining workersPositiveConstruction practitionersPositiveSocioeconomicEconomic developmentEconomic densityPositivePer capita GDPNegativeLocal fiscal revenueNegativeSavings balance of urban and rural residentsNegativeGDP ratio of the primary industryPositivePer capita income of farmersNegativeTertiary industry incomeNegativeSocial capitalInvestment attraction assetsNegativePer capita investment in fixed assetsNegativeGeneral public budget expenditureNegativeSocial securitySocial securityRoad densityNegativeNumber of beds in social welfare institutionsNegativeThe proportion of social security expenditureNegativeNumber of village committees per 10,000 personsNegativeNumber of health institutionsNegativeNumber of medical technicians per 10,000 personsNegativeNumber of beds per 10,000 personsNegativeSocial cultureSocial civilizationNumber of people with high school background per 10,000 personsNegativeTeacher-student ratioNegativeRadio and television coverageNegativeMobile base station coverage rateNegativePer capita mobile phonesNegativeIlliteracy ratioPositiveCustoms and cultureThe ethnic minority population ratioPositiveEthnic culture categoriesPositiveHazard indexDanger of hazardHazard densityPositiveNumber of people threatened by geological disastersPositiveProperty threatened by geological disastersPositiveProtective engineeringProtection engineeringNegativeInvestment in geological hazard controlNegativeSurface protectionSoil and water loss intensityPositiveVegetation coverageNegativeTopographic undulationPositiveThe bold text represents the final selection indexes.The vulnerability evaluation indexes of geological hazards are different at different scales. The social and economic development of the upper Min River is reduced, and data collection is complicated. In the selection of indicators, we must not only follow the scientific principle but also fully consider the reality to ensure the availability of data. 27 indexes were excluded because they could not collect the complete time series data. 16 indexes are selected as vulnerability evaluation factors of the geological disasters in the upper Min River. These indexes include property threatened by geological disasters (D1), number of people threatened by geological disasters (D2), population density (D3), hazard density (D4), teacher-student ratio (D5), economic density (D6), GDP ratio of primary industry (D7), building coverage (D8), road density (D9), number of village committees per 10,000 persons (D10), per capita GDP (D11), per capital investment in fixed assets (D12), general public budget expenditure (D13), per capita mobile phones (D14), number of medical technicians per 10,000 persons (D15), and engineering protection (D16).The basic data of this study are obtained mainly from the following three aspects:(1)
Database on the potential points of geological hazards: According to the geoenvironmental information system of Sichuan province [19] (http://202.61.89.33:16003), information related to coordinates, geographic location, direction, volume, scale, property threatened by geological disasters, number of people threatened by geological disasters, number of households threatened by geological disasters, and prevention engineering of the potential geological hazards in the upper Min River since 2008 was obtained.(2)
Remote sensing images and vector data: Based on the National Earth System, science data-sharing infrastructure [27] (http://www.geodata.cn), the enhanced thematic mapper (ETM) image of the study area, and the essential data related to fields like topography, land use, and earthquakes were obtained.(3)
Data from yearbooks and statistical bulletins: The data on administrative areas, economic output, a permanent population, practitioners of secondary and tertiary industries, industrial output, built-up areas of towns and townships, and other related data were obtained from several yearbooks and statistical bulletins for the years 2006 to 2016. The sources included “China County Economic Statistical Yearbook,” “Sichuan Statistical Yearbook,” “Aba Statistical Yearbook,” “Chengdu Statistical Yearbook,” “Heishui Statistical Yearbook,” “Wenchuan Statistical Yearbook,” “Mao Statistical Yearbook,” “Songpan Statistical Yearbook,” and “Li Statistical Yearbook.” The data not included in yearbooks and statistical bulletins were obtained through field investigation. For the readjustment of the administrative division of Dujiangyan City, at the end of 2014, Zipingpu Town and Hongkou Town were merged into Longchi Town. However, since statistical data are still divided into the yearbooks for three townships, this article does not make any adjustment to it.The indices selected for vulnerability assessment were defined as follows:(1)
Property threatened by geological disasters (D1): geological disaster threatening property refers to the amount of property that may be lost after a geological disaster occurs. The more threatening property, the more property that may be damaged, and the higher the vulnerability, which is a positive index.(2)
The number of people threatened by geological disasters (D2): the number of people threatened by geological disasters is the number of people threatened by geological disasters in each township. The more people threatened by geological disasters, the more people may be damaged by geological disasters, and the higher the vulnerability, which is a positive index.(3)
Population density (D3): population density is the number of people per unit area. Vulnerability is first related to population density. The higher the population density, the higher the vulnerability of settlements, which is a positive index.(4)
Hazard density (D4): the more geological disaster points, the greater the potential threat to residents and the greater the vulnerability, which is a positive index.(5)
Teacher-student ratio (D5): it is the ratio of students to teachers. The more students assigned by each teacher, the weaker the ability to respond to the crisis. The role of the teacher is very important when disaster strikes. The more teachers are relative to the students in the evacuation and transfer project, the more conducive to ensuring the safety of students’ lives, the higher the antirisk ability, and the lower the vulnerability, which is a negative index.(6)
Economic density (D6): economic density refers to the ratio of GDP to the regional area, which characterizes the efficiency of economic activities per unit area of settlements and the intensity of land use. The greater the economic density, the higher the probability of loss in the event of a geological disaster, and the higher the vulnerability of the settlement, which is a positive index.(7)
GDP ratio of primary industry (D7): the GDP ratio of primary industry in the upper Min River mainly refers to agriculture, forestry, animal husbandry, and fisheries. The GDP ratio of the primary industry is an industry with more inputs and less income. Such industries were the leading industries in the upper Min River before 2006, reflecting the society. The higher the degree of closure of the structure, the higher the output value ratio, the higher the vulnerability, which is a positive index.(8)
Building coverage (D8): buildings are also sensitive to geological disasters and are the main disaster-bearing bodies of geological disasters in the upper reaches of the Minjiang River. The higher the building coverage, the higher the vulnerability of settlements, which is a positive index.(9)
Road density (D9): road density refers to the ratio of the length of the road in the area to the area of the area, and it reflects the density of roads. The greater the road density, the higher the probability of damage (89-90), and the higher the vulnerability, which is a positive index.(10)
The number of village committees per 10,000 persons (D10): this is the number of village committees owned by every 10,000 people in the study area. Village committees are grassroots mass organizations for self-service, self-management, and self-education of villagers. They are stable and developed in a small area. The primary organization is the grassroots force for the transmission and implementation of national disaster prevention and reduction policies. The more village committees there are, the better the national disaster prevention and mitigation policies are communicated and implemented, and the lower the vulnerability of settlements, which is a negative index.(11)
Per capita GDP (D11): per capita GDP is mainly used to assess the economic development level of a region. At present, there are differences between urban and rural areas in my country. In terms of economic conditions, the financial status of rural residents is worse than that of urban residents. And the more barren a place is, the lower the danger awareness of geological disasters, and the higher the degree of damage when a disaster occurs. And the more personal wealth, the stronger the ability to reconstruct after the disaster, which can reduce the vulnerability of the system to a certain extent. Therefore, the higher the per capita GDP of a region, the lower the vulnerability, which is a negative index.(12)
Per capita investment in fixed assets (D12): per capita investment in fixed assets is the economic activity of construction and purchase of fixed assets, that is, the reproduction of fixed assets. Investment in fixed assets is the primary means of social fixed asset reproduction. Per capita, fixed assets reflect the strength and ability of disaster reconstruction. The greater the per capita investment in fixed assets, the lower the vulnerability, which is a negative index.(13)
General public budget expenditure (D13): it is the expenditure that the region allocates and uses in a planned manner to the centralized budget revenue, and it is the embodiment of the local grassroots construction capacity. The higher the general public budget expenditure, the better the implementation of disaster prevention and mitigation policies, and the lower the vulnerability. It is a negative index.(14)
Per capita mobile phones (D14): the number of mobile phones is an essential index of regional communication capabilities. The more the number of telephones per capita, the faster the transmission of disaster information, which is more conducive to benefit and avoid harm, and the lower the vulnerability. It is a negative index.(15)
The number of medical technicians per 10,000 persons (D15): it indicates that there are medical technicians for every 10,000 people in the study area, and many places in the mountain regions cannot build a medical system. The number of medical technicians can reflect the local rescue ability. The greater the number of doctors, the better the rescue ability after the disaster and the lower the vulnerability of the settlement, which is a negative index.(16)
Engineering protection (D16): protection engineering refers to the number of construction protection projects to control geological disasters. Engineering protection is an effective measure to protect the safety of local residents' lives and property. Therefore, the greater the number of engineering protection, the lower the vulnerability, which is a negative index.


## 3.3. Evaluation Method

Entropy is an objective, comprehensive evaluation method whose concept originated from thermodynamics [28–31]. It is a measure of uncertainty—the smaller the information entropy of the indicators, the higher the difference coefficient and the weight of the indicators. This means that the information provided by the indicators is large. Thus, it plays an essential role in the comprehensive application. The entropy method is increasingly applied in engineering, management, and social sciences and has gradually become an indispensable practical method in scientific research. Miao and Ding [30] took Lushan and Ludian earthquake-stricken areas as research areas and used the entropy comprehensive evaluation method to determine the social vulnerability of the affected areas. This further proved that the entropy method is practical in social vulnerability assessment. According to the theory and method of information entropy, Ding et al. [32] put forward the concept of “geological hazard entropy” and applied it to the evaluation of geological hazards in the Ranwu-Dongjiu section of Sichuan-Tibet Highway and achieved excellent results.This study combines the entropy method with geographic information system (GIS) to objectively and intuitively analyze the settlement vulnerability of geological hazards in the upper Min River. The main steps of the entropy comprehensive evaluation method are as follows [31]:(a)
Normalization of indexes [31]:Due to the differences in the dimensions, orders of magnitude, units, and quantity changes of the selected indicators in the study, these differences may have an impact on the final analysis results. In order to eliminate these effects, it is necessary to uniformly eliminate the dimensions of the data to solve the problem of the homogenization of qualitative index values. Moreover, because the values of positive and negative indicators represent different meanings (the greater the value of the positive indicator, the greater the vulnerability of mountain settlements, and the smaller the value of the negative indicator, the greater the vulnerability of mountain settlements) Therefore, we use different algorithms to standardize data for positive and negative indicators.The mountain settlements in the upper Min River are divided inton grids, and p indicators are selected. xij is the value of the jth indicator of the ith grid unit (i = 1, 2, …, n; j = 1,2, …, p). The matrix is as follows:(1)X=x11x12⋯x1px21x22⋯x2p⋮⋮⋮xn1xn2⋯xnp=X1,X2,…,Xp.When Xij is a positive indicator [31],(2)Xij∗=Xij−minX1j,X2j,L,XnjmaxX1j,X2j,L,Xnj−minX1j,X2j,L,Xnj,i=1,2,L,n;j=1,2,L,p.When Xij is a negative indicator [31],(3)Xij∗=maxX1j,X2j,L,Xnj−XijmaxX1j,X2j,L,Xnj−minX1j,X2j,L,Xnj,i=1,2,n;Lj=1,2,Lp,wherei represents the sample; j represents the indicator, Xij′ represents the standardized value of the ith sample under the j indicator, minX1j,X2j,L,Xnj represents the minimum value in the j indicator value, and maxX1j,X2j,L,Xnj represents the maximum value in the j indicator value.Normalized data after calculation:(4)X∗=x11∗x12∗⋯x1p∗x21∗x22∗⋯x2p∗⋮⋮⋮xn1∗xn2∗⋯xnp∗=X1∗,X2∗,⋯,Xp∗.(b)
Computation of entropy value [31]:(5)ej=−k∑i=1npijlnpij,whereej represents the entropy value of index j; k > 0; pij=xij∗/∑ni=1xij∗ represents the weight of scheme i under index j; n represents the number of samples.(c)
Computation of the different coefficient [31]:(6)gj=1−ej,wheregj represents the different coefficient of index j and ej is the same as in (5) above.(d)
Computation of the index weight [31]:(8)aj=gj∑j=1ngj,where aj represents the weight of index j and gj is the same as in (6) above.Through ArcGIS software, establish the raster layer of each index, use the above entropy value calculation formula, and use the raster calculator to calculate the vulnerability to obtain the vulnerability value of the mountain settlement. At the same time, in order to conduct a comparative analysis of the vulnerability of mountain settlements in the upper Min River in different periods, the results of mountain settlement evaluations in 2006, 2009, and 2005 were subjected to dispersion standardization (Min-Max standardization). So that the vulnerability index of each period is changed to between [0,1]. The vulnerability level is divided into 5 categories using the equal division method, which are in turn, low (0.0–0.2), comparatively low (0.2–0.4), medium (0.4–0.6), comparatively high (0.6–0.8), and high (0.8–1.0) five mountain settlement vulnerable areas. The standardized transfer function for dispersion is as follows [31]:(9)X′=X−minmax−min,where X′ is the normalized value of the sample data, X is the original value of the sample data, max is the maximum value of the sample data, and min is the minimum value of the sample data.


## 4. Results and Analysis

### 4.1. Computation and Analysis of Vulnerability Indicators

The weights of various indexes changed greatly in 2006, 2009, and 2015 (Table2). In 2006, the weights of road density, economic density, building coverage, and population density were larger, and the weight of per capita GDP and general public budget expenditure was lower. In 2009, the weights of building coverage and road density were higher, and the weight of per capita GDP, general public budget expenditure, and number of medical technicians per 10,000 persons were lower. In 2015, the weights of building coverage and highway density were higher, while the weights of general public budget expenditure, per capita mobile phones, and the number of medical technicians per 10,000 persons were lower.Table 2
Entropy value, difference coefficient, and weight of vulnerability indicators of geological hazards in the upper Min River.
No.Parameter200620092015ejgjajejgjajejgjajD1Population density0.77020.22980.13220.84240.15760.11610.84220.15780.1028D2Building coverage0.76450.23550.13540.72870.27130.19990.60450.39550.2577D3Economic density0.66070.33930.19510.85900.14100.10390.77730.22270.1451D4Per capita GDP0.98500.01500.00860.99500.00500.00370.98430.01570.0102D5GDP ratio of primary industry0.93230.06770.03890.94060.05940.04380.87540.12460.0812D6Per capita investment in fixed assets0.95130.04870.02800.91890.08110.05980.98180.01820.0119D7General public budget expenditure0.99640.00360.00210.99460.00540.00400.99630.00370.0024D8Road density0.66710.33290.19140.76370.23630.17410.78960.21040.1371D9Number of village committees per 10,000 persons0.97640.02360.01350.98280.01720.01270.98010.01990.0130D10Number of medical technicians per 10,000 persons0.96800.03200.01840.99010.00990.00730.99350.00650.0042D11Teacher-student ratio0.96030.03970.02280.97440.02560.01890.95100.04900.0319D12Per capita mobile phone0.94390.05610.03230.95510.04490.03310.99190.00810.0053D13Hazard density0.89550.10450.06010.89320.10680.07870.89320.10680.0696D14Number of people threatened by geological disasters0.89040.10960.06300.90530.09470.06980.90530.09470.0617D15Property threatened by geological disasters0.90360.09640.05540.90390.09610.07080.90390.09610.0626D16Engineering protection0.99520.00480.00280.99520.00480.00360.99520.00480.0031In the evaluation process using the entropy method, the indicators with too small weights can be excluded according to the size of the entropy weight, in order to facilitate a more accurate and reliable evaluation [27]. According to Table 2, teacher-student ratio, number of village committees per 10,000 persons general public budget expenditure, number of medical technicians per 10,000 persons, and engineering protection with aj < 0.03 should be excluded. The 11 indicators finally selected property threatened by geological disasters (N1), the number of people threatened by geological disasters (N2), population density (N3), hazard density (N4), economic density (N5), GDP ratio of primary industry (N6), building coverage (N7), road density (N8), per capita GDP (N9), per capital investment in fixed assets (N10), and per capita mobile phones (N11). Finally, calculate the weight of each indicator in 2006, 2009, and 2015 again (Table 3).Table 3
Final weight of vulnerability indicators of geological hazards in the upper Min River.
N1N2N3N4N5N6N7N8N9N10N11aj20060.05900.06700.14050.06390.20740.04140.14400.20350.00920.02980.034320090.07430.07320.12180.08250.10900.04590.20960.18260.00380.06270.034720150.06630.06530.10880.07360.15350.08590.27260.14500.01080.01260.0056


### 4.2. Results of Vulnerability Assessment and Its Dynamic Changes

As shown in Table4, the medium and low vulnerability areas accounted for the highest proportion between 2006 and 2015. The vulnerable areas changed from medium to low vulnerable areas, which indicates that the implementation of risk-aversion relocation and land-change relocation policies changed the vulnerability of the upper Min River. According to the percentage of each vulnerable area, in the highly vulnerable area, the percentage trend of vulnerable areas was 2006 > 2015 > 2009; in the middle and comparatively high vulnerable areas, the trend was 2006 > 2009 > 2015; and in the comparatively low and low vulnerable areas, the trend of the proportion of lower vulnerability was 2015 > 2009 > 2006. The main reason for this variation was the evacuation and relocation projects of the Wenchuan disaster reconstruction plan that made many people move to dangerous areas.Table 4
Vulnerability zoning results of geological hazards in the upper Min River for 2006, 2009, and 2015.
Vulnerability zoning200620092015Area (km2)Area (%)Area (km2)Area (%)Area (km2)Area (%)High142.166314.8856.82836.9861.779413.35Comparatively high101.227810.59169.200920.7831.25956.75Middle440.935646.14325.55739.9972.763515.72Comparatively low177.231618.54139.737717.16234.849750.73Low94.14319.85122.788215.0862.273713.45In 2006, most settlements in Wenchuan, Mao, and Heishui counties had a high vulnerability to geological hazards—mainly due to their dense population, frequent economic activities, and more potential hazards. Comparatively, higher and middle vulnerable areas were distributed mainly on both sides of the Min River in Wenchuan, Mao, and Songpan counties. The vulnerability of geological hazards in most areas of Li County was relatively low, and the relatively high vulnerability areas were mainly near Zagunao Town. Wolong Town, and Gengda Town in Wenchuan County which are the natural reserves of the giant panda. Because of better vegetation, less human engineering activities, and sparse distribution of geological hazards and settlements, the vulnerability of geological hazards was relatively low (Figure3(a)).Figure 3
Vulnerability zoning of geological hazards in the upper Min River for several years. (a) 2006. (b) 2009. (c) 2015.In 2009, the areas with the high and comparatively high vulnerability of geological hazards in the upper Min River were concentrated mainly in Songpan County, Heishui County, and townships on both sides of the mainstream of Min River in Mao County. Except for the high and comparatively high vulnerability in Songpan County, the upper reaches of the Heishui River Basin and the villages and towns on the right side of the river course were moderately vulnerable. The vulnerability values of the Zagunao River Basin in Li County and the Yuzixi River Basin in Wenchuan County were low as a whole, while those of the Shouxi River Basin were lower (Figure3(b)).In 2015, the settlement vulnerability of geological hazards in the upper Min River was relatively low (Figure3(c)). Nevertheless, the highly vulnerable areas were mainly Fengyi Town in Mao County, Weizhou Town in Wenchuan County, and Jinan Town and Jinan Hui Township in Songpan County. These areas are mainly the government locations of the districts and counties, which are relatively low on both sides of the river and are conducive to settlement production. However, the upper Min River has deep mountains and valleys, which are prone to geological disasters and have a greater potential for vulnerability. Hongkou Township of Dujiangyan has a comparatively high vulnerability. Due to its large population and dense economy, the potential loss is enormous. Except for the high and comparatively high vulnerability in Songpan County, the upper reaches of the Heishui River Basin and the villages and towns on the right side of the river course were moderately vulnerable. The vulnerability values of the Zagunao River Basin and the Yuzixi River Basin were low as a whole.


## 4.1. Computation and Analysis of Vulnerability Indicators

The weights of various indexes changed greatly in 2006, 2009, and 2015 (Table2). In 2006, the weights of road density, economic density, building coverage, and population density were larger, and the weight of per capita GDP and general public budget expenditure was lower. In 2009, the weights of building coverage and road density were higher, and the weight of per capita GDP, general public budget expenditure, and number of medical technicians per 10,000 persons were lower. In 2015, the weights of building coverage and highway density were higher, while the weights of general public budget expenditure, per capita mobile phones, and the number of medical technicians per 10,000 persons were lower.Table 2
Entropy value, difference coefficient, and weight of vulnerability indicators of geological hazards in the upper Min River.
No.Parameter200620092015ejgjajejgjajejgjajD1Population density0.77020.22980.13220.84240.15760.11610.84220.15780.1028D2Building coverage0.76450.23550.13540.72870.27130.19990.60450.39550.2577D3Economic density0.66070.33930.19510.85900.14100.10390.77730.22270.1451D4Per capita GDP0.98500.01500.00860.99500.00500.00370.98430.01570.0102D5GDP ratio of primary industry0.93230.06770.03890.94060.05940.04380.87540.12460.0812D6Per capita investment in fixed assets0.95130.04870.02800.91890.08110.05980.98180.01820.0119D7General public budget expenditure0.99640.00360.00210.99460.00540.00400.99630.00370.0024D8Road density0.66710.33290.19140.76370.23630.17410.78960.21040.1371D9Number of village committees per 10,000 persons0.97640.02360.01350.98280.01720.01270.98010.01990.0130D10Number of medical technicians per 10,000 persons0.96800.03200.01840.99010.00990.00730.99350.00650.0042D11Teacher-student ratio0.96030.03970.02280.97440.02560.01890.95100.04900.0319D12Per capita mobile phone0.94390.05610.03230.95510.04490.03310.99190.00810.0053D13Hazard density0.89550.10450.06010.89320.10680.07870.89320.10680.0696D14Number of people threatened by geological disasters0.89040.10960.06300.90530.09470.06980.90530.09470.0617D15Property threatened by geological disasters0.90360.09640.05540.90390.09610.07080.90390.09610.0626D16Engineering protection0.99520.00480.00280.99520.00480.00360.99520.00480.0031In the evaluation process using the entropy method, the indicators with too small weights can be excluded according to the size of the entropy weight, in order to facilitate a more accurate and reliable evaluation [27]. According to Table 2, teacher-student ratio, number of village committees per 10,000 persons general public budget expenditure, number of medical technicians per 10,000 persons, and engineering protection with aj < 0.03 should be excluded. The 11 indicators finally selected property threatened by geological disasters (N1), the number of people threatened by geological disasters (N2), population density (N3), hazard density (N4), economic density (N5), GDP ratio of primary industry (N6), building coverage (N7), road density (N8), per capita GDP (N9), per capital investment in fixed assets (N10), and per capita mobile phones (N11). Finally, calculate the weight of each indicator in 2006, 2009, and 2015 again (Table 3).Table 3
Final weight of vulnerability indicators of geological hazards in the upper Min River.
N1N2N3N4N5N6N7N8N9N10N11aj20060.05900.06700.14050.06390.20740.04140.14400.20350.00920.02980.034320090.07430.07320.12180.08250.10900.04590.20960.18260.00380.06270.034720150.06630.06530.10880.07360.15350.08590.27260.14500.01080.01260.0056


## 4.2. Results of Vulnerability Assessment and Its Dynamic Changes

As shown in Table4, the medium and low vulnerability areas accounted for the highest proportion between 2006 and 2015. The vulnerable areas changed from medium to low vulnerable areas, which indicates that the implementation of risk-aversion relocation and land-change relocation policies changed the vulnerability of the upper Min River. According to the percentage of each vulnerable area, in the highly vulnerable area, the percentage trend of vulnerable areas was 2006 > 2015 > 2009; in the middle and comparatively high vulnerable areas, the trend was 2006 > 2009 > 2015; and in the comparatively low and low vulnerable areas, the trend of the proportion of lower vulnerability was 2015 > 2009 > 2006. The main reason for this variation was the evacuation and relocation projects of the Wenchuan disaster reconstruction plan that made many people move to dangerous areas.Table 4
Vulnerability zoning results of geological hazards in the upper Min River for 2006, 2009, and 2015.
Vulnerability zoning200620092015Area (km2)Area (%)Area (km2)Area (%)Area (km2)Area (%)High142.166314.8856.82836.9861.779413.35Comparatively high101.227810.59169.200920.7831.25956.75Middle440.935646.14325.55739.9972.763515.72Comparatively low177.231618.54139.737717.16234.849750.73Low94.14319.85122.788215.0862.273713.45In 2006, most settlements in Wenchuan, Mao, and Heishui counties had a high vulnerability to geological hazards—mainly due to their dense population, frequent economic activities, and more potential hazards. Comparatively, higher and middle vulnerable areas were distributed mainly on both sides of the Min River in Wenchuan, Mao, and Songpan counties. The vulnerability of geological hazards in most areas of Li County was relatively low, and the relatively high vulnerability areas were mainly near Zagunao Town. Wolong Town, and Gengda Town in Wenchuan County which are the natural reserves of the giant panda. Because of better vegetation, less human engineering activities, and sparse distribution of geological hazards and settlements, the vulnerability of geological hazards was relatively low (Figure3(a)).Figure 3
Vulnerability zoning of geological hazards in the upper Min River for several years. (a) 2006. (b) 2009. (c) 2015.In 2009, the areas with the high and comparatively high vulnerability of geological hazards in the upper Min River were concentrated mainly in Songpan County, Heishui County, and townships on both sides of the mainstream of Min River in Mao County. Except for the high and comparatively high vulnerability in Songpan County, the upper reaches of the Heishui River Basin and the villages and towns on the right side of the river course were moderately vulnerable. The vulnerability values of the Zagunao River Basin in Li County and the Yuzixi River Basin in Wenchuan County were low as a whole, while those of the Shouxi River Basin were lower (Figure3(b)).In 2015, the settlement vulnerability of geological hazards in the upper Min River was relatively low (Figure3(c)). Nevertheless, the highly vulnerable areas were mainly Fengyi Town in Mao County, Weizhou Town in Wenchuan County, and Jinan Town and Jinan Hui Township in Songpan County. These areas are mainly the government locations of the districts and counties, which are relatively low on both sides of the river and are conducive to settlement production. However, the upper Min River has deep mountains and valleys, which are prone to geological disasters and have a greater potential for vulnerability. Hongkou Township of Dujiangyan has a comparatively high vulnerability. Due to its large population and dense economy, the potential loss is enormous. Except for the high and comparatively high vulnerability in Songpan County, the upper reaches of the Heishui River Basin and the villages and towns on the right side of the river course were moderately vulnerable. The vulnerability values of the Zagunao River Basin and the Yuzixi River Basin were low as a whole.


## 5. Discussion

The main indicators of the vulnerability of geological hazards in the upper Min River were population density, building coverage, economic density, and road density. They indicate that population, buildings, and highways are the main disaster-bearing bodies of geological hazards in the upper Min River and that regional economic value is the most direct loss of disaster threat objects. Among the factors contributing less, GDP per capita and general public budget expenditure were relatively low in 2006. The reasons that accounted for this were the mountain areas in the upper Min River, poor economic foundation, lack of revenue sources for the government, low corresponding public budget expenditure, and internal differences. In 2009, per capita GDP, general public budget expenditure, and the number of medical technicians per 10,000 persons accounted for the lowest share. Compared with 2006, after the Wenchuan earthquake, there was an increase in the number of medical technicians per 10,000 persons, showing higher medical and health investment by the government in remote areas. As a result, the regional differences in the distribution of doctors and technicians narrowed. More importantly, although the share of general public budget expenditure was still small, it was quite different from the situation in 2006. The government still invested less in the area in 2006, but the state invested a lot of money after the Wenchuan earthquake, and with support from all parts of the country, the general public budget input in 2009 was very high compared with that in 2006.Compared with 2006 and 2009, in 2015, the weight of per capita mobile phones was relatively lower. The study area still had mainly fixed phones and functional computers in 2006 and 2009, and the mobile infrastructure was not perfect. With imperfection in mobile infrastructure, there was a big difference in the per capita telephones in each town. In 2015, due to the government’s construction of basic communication facilities in the region, coupled with the popularization of 3G and 4G signals, smartphones entered thousands of mountain families, and now almost every household has a mobile phone.The average vulnerability of villages and towns in the upper Min River was 0.2213 in 2006, 0.2290 in 2009, and 0.2251 in 2015. The order of the vulnerability value was 2009 > 2015 > 2006, and the overall vulnerability of the upper Min River was low. The main reason for the highest vulnerability value in the area in 2009 was earthquake-stricken areas of Wenchuan, where production and life were seriously affected by earthquakes and the secondary geological hazards caused by the earthquake.The vulnerability value in 2015 was much lower than that in 2009, but it was still higher than that in 2006, which indicates that the reconstruction effect of the Wenchuan earthquake was remarkable, but there was a need for strengthening disaster prevention and mitigation mechanism. Before 2015, five counties in the upper Min River were deeply poverty-stricken (among them Heishui County was a state-level poverty-stricken county). With the deepening implementation of the precise poverty alleviation policy, the economic development of the region has been rapid in the past 10 years. Especially after 2008, the national postdisaster reconstruction has invested a lot of manpower and material resources in the region, bringing Mao, Wenchuan, and Li counties out of poverty successively in 2017, which is consistent with the results of vulnerability study.According to the ascending order of vulnerability values of 84 townships in the upper Min River in 2006, 2009, and 2015 (Table5), the time-varying curve of settlement vulnerability of geological hazards can be obtained (Figure 4). From the trend of change, the vulnerability curve is 76–78 in villages and towns, which divides the whole area into two parts: high value and low-value areas. The vulnerability trend of high value area is 2009 > 2015 > 2006, and that of low value area is 2009 > 2006 > 2015. In each year, the proportion of towns in the low-value areas is higher, while that of the high-value area is lower. Vulnerability rises faster in the high-value area and slower in the low-value area, which is like the “Matthew effect” in economics. This special “Matthew effect” can better reflect the changing trend of vulnerability on a time scale. That is, the settlement vulnerability of geological hazards has polarized—the areas with high vulnerability value continue to rise and the areas with low vulnerability value continue to decrease. This also indirectly indicates that we should focus on building a disaster prevention and emergency response capacity in highly vulnerable areas.Table 5
Vulnerability of townships to geological hazards in the upper Min River for 2006, 2009, and 2015.
No.Vulnerability in 2006TownshipVulnerability in 2009TownshipVulnerability in 2015Township10.0000Jiabi0.0000Jiabi0.0000Jiabi20.0268Miyaluo0.0390Sanjiang0.0211Miyaluo30.0362Guergou0.0459Miyaluo0.0398Guergou40.0393Shangmeng0.0526Shangmeng0.0419Shangmeng50.0545Puxi0.0572Guergou0.0488Puxi60.0623Putou0.0597Gengda0.0570Weicheng70.0634Tonghua0.0649Yinxing0.0610Songpinggou80.0638Sanjiang0.0821Putou0.0664Putou90.0924Yinxing0.0829Weicheng0.0672Sanjiang100.0960Gengda0.0860Puxi0.0721Tonghua110.0971Taoping0.0955Tonghua0.0755Diexi120.1022Xuecheng0.1019Wolong0.0766Yinxing130.1042Xiameng0.1062Songpinggou0.0813Wolong140.1073Wolong0.1102Diexi0.0819Gengda150.1129Ganbao0.1163Taoping0.0850Wadi160.1172Muka0.1197Longxi0.0851Kalong170.1223Caopo0.1218Xiameng0.0879Taoping180.1319Xiabazhai0.1305Ganbao0.0951Taiping190.1350Yanyun0.1327Shidaguan0.0984Sanlong200.1368Weicheng0.1353Xuecheng0.0997Caopo210.1404Shangbazhai0.1377Heihu0.1007Muka220.1408Longxi0.1396Wadi0.1017Shidaguan230.1439Hongzha0.1412Kalong0.1017Xiameng240.1442Songpinggou0.1424Caopo0.1026Xuecheng250.1446Yanmen0.1424Sanlong0.1035Heihu260.1493Caoyuan0.1624Qugu0.1040Shashiduo270.1497Yingxiu0.1626Yadu0.1075Ganbao280.1500Dazhai0.1626Shashiduo0.1171Longxi290.1502Shuijing0.1628Taiping0.1185Yanmen300.1504Mouni0.1653Baixi0.1198Luoduo310.1609Shanba0.1737Yanmen0.1249Qugu320.1656Diexi0.1739Muka0.1265Yadu330.1685Taiping0.1764Luoduo0.1291Baixi340.1698Miansi0.1800Ciba0.1291Ciba350.1705Shashiduo0.1806Miansi0.1343Longba360.1709Hongkou0.1861Longba0.1407Huilong370.1712Daxing0.1905Yanyun0.1412Yingxiu380.1754Heihu0.1962Shangbazhai0.1421Zhimulin390.1783Hongtu0.1990Goukou0.1423Zhawo400.1790Wadi0.1994Hongzha0.1461Waboliangzi410.1819Shidaguan0.1994Xiabazhai0.1467Qinglang420.1823Sanlong0.2009Zhimulin0.1530Mawo430.1834Xiaoxing0.2023Waboliangzi0.1552Keku440.1888Kalong0.2033Shuijing0.1576Goukou450.1890Mawo0.2037Mouni0.1632Shidiaolou460.1915Luoduo0.2039Qinglang0.1764Shuangliusuo470.1953Ciba0.2063Shidiaolou0.1828Weigu480.1964Zhenping0.2112Nanxin0.1852Hongyan490.1968Baixi0.2114Zhawo0.1853Feihong500.1968Longba0.2116Caoyuan0.2044Weimen510.2042Keku0.2128Shanba0.2193Musu520.2064Yadu0.2148Dazhai0.2538Miansi530.2064Longchi0.2267Huilong0.2578Hongkou540.2084Qinglang0.2278Hongtu0.2672Xiabazhai550.2131Chuanzhusi0.2298Mawo0.2710Longchi560.2149Shuangliusuo0.2306Yingxiu0.2721Hongzha570.2162Zhawo0.2306Daxing0.2721Yanyun580.2162Zhimulin0.2322Feihong0.2750Dazhai590.2180Minjiang0.2354Keku0.2752Shangbazhai600.2241Qugu0.2401Xiaoxing0.2794Shuijing610.2245Hongyan0.2437Hongyan0.2817Mouni620.2299Shidiaolou0.2510Weigu0.2835Daxing630.2323Shili Hui0.2545Shuangliusuo0.2886Shanba640.2325Zhenjiangguan0.2547Seergu0.2897Caoyuan650.2328Goukou0.2597Weimen0.3040Hongtu660.2354Waboliangzi0.2636Minjiang0.3062Xiaoxing670.2368Anhong0.2660Zhenping0.3177Zhenping680.2439Huilong0.2769Chuanzhusi0.3200Minjiang690.2546Qingyun0.2797Anhong0.3293Nanxin700.2571Musu0.2807Musu0.3411Zhenjiangguan710.2716Weigu0.2809Zhenjiangguan0.3458Anhong720.2745Zagunao0.3028Shili Hui0.3467Chuanzhusi730.2825Feihong0.3087Hongkou0.3542Zipingpu740.2825Seergu0.3204Shuimo0.3609Qingyun750.3040Weimen0.3289Qingyun0.3616Xuankou760.3129Nanxin0.3461Luhua0.3699Shili Hui770.3930Shuimo0.3475Longchi0.3870Seergu780.5282Jinan Hui0.4193Xuankou0.3874Shuimo790.6186Luhua0.4935Zipingpu0.5389Luhua800.6978Jinan0.5829Jinan Hui0.6135Zagunao810.7608Fengyi0.6386Zagunao0.8388Fengyi820.7730Weizhou0.8669Jinan0.8412Jinan Hui830.8081Zipingpu0.9191Weizhou0.9766Jinan841.0000Xuankou1.0000Fengyi1.0000WeizhouThe bold text indicates the transition zone from low‐value vulnerability area to high-value area.Figure 4
Vulnerability change of geological hazards in the upper Min River.


## 6. Conclusions

Vulnerability, as a part of risk assessment, is an indispensable part of regional risk control. Wenchuan earthquake, geological hazards, settlement migration, socioeconomic conditions, and ethnic cultural changes have a profound impact on the mountain settlement vulnerability of geological hazards in the upper Min River. The migration of mountain settlements in this area varies greatly. Due to national policy, the area of settlements has been drastically reduced, and the population living in traditional farming and animal husbandry has gradually moved to cities. As main destinations of settlement migration, the county towns of each district and county should strengthen the capacity building of disaster prevention and emergency response.The contribution share of each index is different in different years. Population density, building coverage, economic density, and road density are the main factors affecting the settlement vulnerability of geological hazards in the upper Min River. Therefore, strengthening the construction of the regional population, economy, and infrastructure is conducive to reducing the vulnerability of the area.The settlement vulnerability of geological hazards in the upper Min River is witnessing a dynamic change, with the highest vulnerability in 2009, the second in 2015, and the lowest in 2006. Although the Wenchuan earthquake caused considerable damage in the area, significant achievements have also been made in the postearthquake reconstruction. The time fluctuation of vulnerability distribution in villages and towns is consistent. The trend of vulnerability distribution is basically in line with the “Matthew effect.” The areas with high vulnerability to settlement geological hazards will continue to increase. Therefore, we should focus on the prevention of geological hazards, the creation of an investment in reconstruction, and building emergency capacity in highly vulnerable areas.



---

*Source: 8887487-2020-08-28.xml*