Complexity / 2020 / Article
Special Issue

New Models, New Technologies, New Data and Applications of Urban Complexity from Spatio-temporal Perspectives

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 3484568 |

Li Wu, Binggeng Xie, Xiao Xiao, Bing Xue, Jingzhong Li, "Classification Method and Determination of Mountainous Area Types at Township Scales: A Case Study of Yuxi City, Yunnan Province", Complexity, vol. 2020, Article ID 3484568, 13 pages, 2020.

Classification Method and Determination of Mountainous Area Types at Township Scales: A Case Study of Yuxi City, Yunnan Province

Academic Editor: Jianhong (Cecilia) Xia
Received07 Jun 2020
Revised07 Aug 2020
Accepted23 Aug 2020
Published03 Sep 2020


The high-resolution regional division of mountainous area types has important scientific significance for promoting precise management of land space and regional sustainable development. At present, the classification method of mountainous area types is mainly at the county level and above, while classifications for towns and villages are nearly nonexistent, which poses a technical challenge for rural revitalization and the construction of ecological civilization in mountainous areas. We used Yuxi city, Yunnan Province, as the target area of this research, which was based on GIS technology and Digital Elevation Model (DEM) data and socioeconomic environmental monitoring data. The most appropriate statistical unit (e.g., 2.8224 km2) for topographic relief was defined, and the study area was divided into six mountain types: flatlands, hills, low mountains, medium-low mountains, midmountains, and subhigh mountains. Based on the township scale, an index system and classification method dominated by the plain comprehensive index were established to carry out mountain area classifications at township scales. The 75 towns of Yuxi city can be classified into 27 plain towns, 23 mountain-plain towns, and 25 mountain towns from an empirical application perspective, which can provide strong data support and a reference basis for studying the evolution characteristics of land use in different geographical spaces and their interrelationships as well as differentiated land space planning and governance.

1. Introduction

As important regions that support water, energy, and biological diversity on earth [1], mountains occupy 24% of the global land surface area, provide homes for 12% of the world’s population, indirectly offer goods and services for more than 40% of the global population, and play a very important role in the survival and development of humans [2, 3]. In 1992, chapter 13 of Agenda 21, which was titled “Managing fragile ecosystems: sustainable mountain development,” was released by the United Nations Conference on Environment and Development in Rio de Janeiro, while in November 1998, the United Nations General Assembly designated the year 2002 as the International Year of Mountains [3, 4]. From then on, mountains began to receive increasing attention, and development and construction in mountainous areas, as a common topic for mankind, has been highly valued by scholars at home and abroad [58]. As the natural foundation and core connotation of mountainous areas, mountains have some similarities and some differences with mountainous areas. Mountains are complex combinations of various geomorphological elements which focus on differences of landforms and geomorphology in nature [9, 10]. However, mountainous areas emphasize regional combinations of humanity and geomorphology [11, 12] and interactions between humans and the geographical environment which is the research carrier for exploring the relationship between humans, land, and sustainable development [13]. However, with the rapid development of cities, which are the main study areas of human-earth relations [14, 15], mountainous areas have gradually become marginal areas and problem areas for economic development [16, 17] due to the fragility of mountain ecosystems and excessive human activity, soil erosion, deterioration of the ecological environment, economic backwardness, population poverty, and other problems [1820]. The scientific definition of mountains and mountainous area types can provide a basis for classified development and classified policy, which are conducive to promoting the construction of ecological civilization and sustainable development in mountainous areas.

At present, many scholars have conducted research on the extraction of mountain information and classification methods of mountainous area types. The classification method of mountainous area types mainly uses the proportion of mountain area as the classification index, so the first task is to accurately define mountains types and their scopes [21]. Classification of mountain types is generally based on their geomorphologic elements such as topographic relief, elevation, and slope. In specific applications, some classifications use a single index to divide, while some classifications use a combination of multiple indicators to divide; compared with the use of a single index to divide mountain types, a comprehensive division based on multiple indicators can truly reflect mountain characteristics [22]. In the process of defining mountain types, remote sensing, geographic information technology, and DEM are mainly used to extract terrain factor information such as topographic relief and slope [23, 24]. Topographic relief, as the most important topographic factor for classifying mountain types [11, 25], is mainly influenced by determination of its optimal statistical unit [26, 27], which is generally determined by using mean-point analysis methods with the support of GIS software [28, 29]. However, results from the same region by using the same method are also different due to different levels of precision. For example, the optimal statistical unit of relief in Sichuan Province is calculated as 2.34 km2 [30] or 0.2025 km2 [31]; for the Hengduan Mountainous Region, this parameter is calculated as 14.98 km2 [32] or 2.62 km2 [33], and for the Tibetan Plateau, it is calculated as 1.17 km2 [34], 5.06 km2 [35], or 2.25 km2 [33]. It can be seen that the results of calculating the optimal statistical unit area of topographic relief in the same region are quite different; one of the reasons is the selection of the number of experimental units, which is further limited by the workload of manual extraction. Based on the Python programming language, the automatic extraction of experimental windows can be achieved, and the number of experimental units can be greatly increased, which is expected to improve the accuracy of the optimal statistical unit so that mountain types can be scientifically and reasonably defined based on the geomorphic characteristics of the study area.

Yuxi city in Yunnan Province, which has more mountains and less flatland, is located in the mountainous area of southwestern China, and these types of land resources are also a significant factor that affect the spatial differences in regional poverty [36]. As intermontane basins, river terraces, or relatively gentle piedmont valleys in the southwest mountainous area [37], plain areas have flat terrain and are conducive to construction and reclamation when compared with mountainous areas. Therefore, both productivity and land bearing capacity are relatively high in plain areas and are not only important agricultural production areas but are also gathering places for populations, economies, and cultures. However, in mountainous areas, natural conditions are harsh, loss of sand and soil is serious, and the level of social and economic development is low so the differences between these two geographical spaces are very clear. To achieve coordinated development of two different geographic spaces, scientific classification of mountainous areas is the prerequisite and foundation. The existing research on identification and classification of mountainous area types has created a relatively unified system at the macro scale, namely, the area proportions of mountain types are used [27, 38]. Scholars have divided mountainous and plain areas at the scale of provinces and counties, respectively, which has provided a certain theoretical basis for land use and protection of two different geographical spaces [39, 40]. However, with more detailed research, further refinement of spatial governance scales and further strengthening of the operability of local policies and with the implementation of the national “rural revitalization” policy, differentiated territorial space planning should proceed with more detail into the micro township scales. Hence, there is an urgent need to further subdivide mountainous area types at the microscale. In addition to dividing by the area proportions of mountain types as an important dividing index, superposition of human activities in land exploitation and utilization should also be considered to establish a comprehensive index system to conduct our high-precision classifications.

Therefore, on the basis of previous studies and by the combination of topographic and geomorphic data of Yuxi city, this study intends to increase the number of experimental units to 198 based on Python programming with the use of ArcGIS10.6 software to achieve automatic extraction of topographic relief. Then, the optimal statistical unit obtained is compared with relevant studies in an attempt to explore formation rules. At the same time, a comprehensive index system is constructed based on the proportion of mountain area, sum of flatland cultivated land and construction land area, and the comprehensive standardized value of flatland numbers and areas to divide the types of mountainous and plain areas at the scale of towns and villages. This classification in mountainous area types can provide data support for follow-up discussions on social and economic relationships between mountainous areas and plain areas and provide a scientific reference basis for the implementation of differentiated territorial space planning and precise governance.

2. Materials and Methods

2.1. Regional Overview

Yuxi city is located in the central region of Yunnan province. Its geographical coordinates are 23°19′∼24°53′N and 101°16′∼103°09′E. It is adjacent to the provincial capital, Kunming, which is to the northeast; Chuxiong Autonomous Prefecture in the north; Pu’er city in the southwest and Honghe Autonomous Prefecture in the southeast, as is shown in Figure 1. The city has a maximum distance of 172 km from east to the west and maximum distance of 163.5 km from north to south with a total land area of 1.49 × 104 km2. It has jurisdiction over 7 counties, 2 districts, and 75 townships.

Yuxi city is located on the Yunnan-Guizhou Plateau with high terrain in the Northwest and low terrain in the Southeast. The terrain is complex with staggered distributions of mountains, valleys, plateaus, and basins. Ailao Mountain to the West is a great barrier with mountains and valleys, and there are also some large fault collapse basins in the East and North. Due to the deep cutting of the Yuanjiang River Valley along the Yuanjiang fault zone on the east side of the Ailao Mountains, elevation differences greater than 2,000 m from the river to the top of the mountain and a mountain valley area were formed. Except for the Yuanjiang River, most of the research region is 1,500–1,800 m above sea level. The climate type is mainly a subtropical subhumid plateau monsoon climate without hot summers and cold winters with average annual temperatures of 17.4–23.8°C and average annual precipitation of 670–2,412 mm. There are three plateau faulted lakes with abundant rainfall and high levels of farming which enjoy the reputation of “granary in central Yunnan” and “hometown of Yunnan tobacco” due to high grain yields and good quality of tobacco leaves.

2.2. Data Source and Processing

The DEM data used in the study are ASTER GDEMV2 30M, and these data were obtained from the Geospatial Data Cloud website ( ArcGIS 10.6 was applied to splice projections, and the administrative boundary of Yuxi city was then selected to obtain the DEM data for Yuxi city. The land vector data were obtained from the Resource and Environment Data Center of the Chinese Academy of Sciences ( Grid transformation was carried out with ArcGIS 10.6, and the land vector data of each town in Yuxi city were obtained by selecting the boundaries of towns and villages in Yuxi city.

2.3. Methods
2.3.1. Mountain Type Definition Method

Considering the mountain classification index system proposed by the United Nations Environment Programme and relevant scholars [12, 22, 41], high-altitude areas are generally directly classified as mountains, while middle- and low-altitude areas should be defined as mountain types in combination with slopes, local elevation changes, and relief. According to the DEM used in the experiment, Yuxi city is located in the southwest mountainous area, and its elevation range is between 266 and 3,139 m. Statistically, areas up to 1,500–2,000 m sea level account for 49.15% of the city and form nearly half of the entire research area. Through slope analysis, it was found that the average slope is 21° and that mountains are widely distributed. The highest relief degree of the land surface (RDLS) is 1,335 m, with little overall undulation, but local undulation is great. Therefore, in the classification of mountainous area types at the township scale, RDLS should be considered as an important index and be combined with altitude and slope to classify mountain types. When the altitude is below 1,500 m, mountain types should be divided according to the degree of relief. The influence of slope factor should also be considered when the altitude is between 1,500 m and 3,500 m. Therefore, the mountains in Yuxi city are divided into hills, low mountains, medium-low mountains, midmountains, and subhigh mountains (Table 1).

ElevationRDLSSlopeMountain types

0–500>100 mHill
500–1000>300 m within the scope of 2.82 km2Low mountain
1000–1500>300 m within the scope of 2.82 km2Medium-low mountain
2500–3500>5°Subhigh mountain

2.3.2. Topographic Relief Calculation

At present, most scholars extract topographic relief by adopting window analysis methods based on the spatial analysis module in ArcGIS10.6. The level of topographic relief in a window can be obtained by calculating the difference between the maximum and minimum values of focus statistics for different window pixels n × n (n = 2, 3, …, 198, 199). Based on a window movement step distance of 1, automatic extraction of terrain relief for different windows can be achieved by programming Python modules.

For calculations of topographic relief, selection of the optimal size of statistical unit (the optimal moving window) is the key. In this study, mean variable point analysis was used to determine the optimal statistical unit. The basic principle of this method is based on the existence of only one inflection point which changes from steep to gradual. The mathematical expression is that there are data sequences , and among them, i = 1, 2, 3, …, N(1, 2, …, 198), where N is the number of samples and the sample is divided into two segments bounded by the point xi and then the arithmetic mean value of each section , and sample mean are calculated separately:where Si and S are intermediate calculated values t1 = 1, 2, …, i − 1 and t2 = i, i + 1, …, N. Then, the expected value can be calculated and the window size corresponding to the maximum expected value is the optimal statistical unit area.

2.3.3. Classification Methods for Mountainous and Plain Areas

Since there is no reference in existing studies to the classification method of mountainous area and plain area types based on the township scale, this study refers to the method of mountain area proportion based on the county scale, but there are improvements to two indicators of flatland-farmland area proportions, and number of flatlands (area ≥ 100 km2) [39]. With the rapid development of the social economy and constant increases of urbanization levels, large amounts of flatland cultivated land have been occupied by construction [42, 43]. Therefore, there is a certain deficiency when using only the proportion of flatland cultivated land for evaluations. In this study, the superposition of flatland construction land area proportions and flatland cultivated land area proportions is taken as an index, namely, the proportion of flatland cultivated land and construction land. At the same time, it is difficult for the area of a single flatland to reach 100 km2 at the township scale, so flatlands that are greater than or equal to 1 km2 will be chosen as the standard according to the definition of Agricultural Geography of Yunnan [44], field research, and data analysis. Meanwhile, it is unfair to evaluate only using the number of flatlands as the size of the flatland area has more practical significance for the division. Therefore, the standardized value of flatland (area ≥ 1 km2) in every village or town will be obtained by combining the number of flatlands and the area of flatlands. The formulas are as follows:where is the plain composite index, is the flatland farmland area, is the flatland construction land area, is the total farmland area, is the total construction land area, is the proportion of flatland farmland and construction land, is the flatland land area, is the total land area, is the proportion of flatland area, is the maximum proportion of flatland farmland and construction land, is the maximum proportion of flatland land area, is the standardized value of flatland (area ≥1 km2) in every villages or towns, is the standardized value of every flatland (area ≥1 km2), and m is the number of flatlands (area ≥ 1 km2); , , and , are the weights of , , and , respectively; and to ensure that the result of is between 0 and 100, with values of 50, 35, and 15, respectively.

Combined with the leading index and auxiliary index, the specific dividing standard for plain areas and mountainous areas is shown in Table 2.

Classification indicatorClassification criteria
Plain townsMountain-plain townsMountain towns

Leading index≥5020–50<20

Auxiliary index (%)≥2010–20<10

3. Results

3.1. Optimum Statistical Unit

By use of Python module programming and 2 × 2, 3 × 3, …, 198 × 198, to 199 × 199 windows, the topographic relief was automatically extracted, and different window sizes and their corresponding average topographic relief values were obtained, as shown in Table 3.

Window sizeArea (km2)Average relief (m)

2 × 20.003613.3291
3 × 30.008126.3569
4 × 40.014439.0506
5 × 50.022551.3448
6 × 60.032463.2301
7 × 70.044174.6776
8 × 80.057685.7005
9 × 90.072996.3601
10 × 100.0900106.5712
11 × 110.1089116.5744
53 × 532.5281387.8541
54 × 542.6244391.2082
55 × 552.7225396.9445
56 × 562.8224400.1518
57 × 572.9241406.2713
150 × 15020.25699.0040
151 × 15120.5209699.7820
152 × 15220.7936706.0452
194 × 19433.8724789.0474
195 × 19534.2225789.3074
196 × 19634.5744802.8152
197 × 19734.9281802.5625
198 × 19835.2836801.0867
199 × 19935.6409805.1577

Origin 2019 software was used to perform logarithmic equation fitting for the corresponding statistical areas and average topographic relief for each window shown in Table 3. The fitting curve of relief degree in Yuxi city was obtained. The equation used is , and the decision coefficient is (Figure 2). The fitting effect was good and passed the statistical test. Figure 2 shows that the average topographic relief in the research area increased with increases in window area. Below a window area of 2 km2, the average topographic relief increased greatly and the curve was steeper. However, between 2 and 4 km2, the curve changed from steep to gradual, and above 4 km2, the increase of the average topographic relief decreased and the curve became flatter. According to formulas (1) and (2), the Si, S and S − Si values and differential fitting curves of Si and S can be calculated (Table 4, Figure 3). Table 4 and Figure 3 show that the difference between Si and S reached a maximum value at the 55th point, i.e., the turning point where the curve changed from steep to gradual. Thus, a 56 × 56 window (2.8224 km2) is the optimal statistical unit area for the topographic relief of Yuxi city.

iS iS − Si


3.2. Definition of Mountain Types

Based on ASTER GDEM data, the elevations (Figure 4(a)) and slopes (Figure 4(b)) were extracted, and the topographic relief was extracted by using a 56 × 56 window (Figure 4(c)). The classification of mountain types in Yuxi city in ArcGIS10.6 was accomplished by using the previously mentioned mountain type definition indicators, and the results are shown in Figure 4(d). According to Figure 4(d), Yuxi city is mainly dominated by mountains, and the flatland is very small. In the defining results of various mountain types, the hilly area is 13,294.32 hm2, the low mountain area is 121,727.41 hm2, the medium-low mountain area is 250,393.05 hm2, the midmountain area is 820,116.80 hm2, and the subhigh mountain area is 20,460.24 hm2. The highest percentage of each mountain type in the total land area is midmountains (54.96%), which is followed by medium-low mountains (16.78%), low mountains (8.16%), subhigh mountains (1.37%), and hills (0.89%). Flatlands are mainly distributed in some large fault collapse basins in the east, namely, the Hongta District, Tonghai County, Jiangchuan County, and Chengjiang County. The elevations are between 1,500 and 2,500 m, but the slopes are relatively low (below 8°) and the topographic relief is less than 300 m (Figures 4(a)4(c)). Hills are mainly distributed in the part of Yuanjiang County that is located to the west of Yuxi city, in a valley with lower elevations. Xinping County also has a small amount of hills distribution, which is staggered with some flatlands. Because the elevations are between 500 and 1,000 m, the low mountains are mainly distributed around the southwest river valley in Xinping and Yuanjiang but tend to spread to the periphery of E’shan and Yimen. In addition to Xinping, Yuanjiang, and E’shan, medium-low mountains are widely distributed and extend from the west to the east along the southern periphery of each county, but there are no distributions in the north. With elevations of 1,500–2,500 m and as the largest distribution type, midmountains are distributed in all counties and districts, especially in the east and central regions. Subhigh mountain areas are less common and are mainly distributed in the Ailao Mountain area of Xinping but other counties such as Yimen, Chengjiang, and Huaning also have sporadic distributions.

3.3. Classification of Mountainous Areas and Plain Areas

With the help of formulas (3)–(6), the proportion of total flatland cultivated land and construction land area, proportion of flatland land area, and number and area of flatlands (area ≥1 km2) for each town in Yuxi city were calculated, and the plain comprehensive index was also calculated (Table 5, Figures 5(a)5(c)). According to the standard for dividing plain areas and mountainous areas in Yuxi city, the plain comprehensive index was used as the leading index and the other three indexes were used as auxiliary indexes; the 75 villages and towns were divided into plain areas, mountain-plain areas, and mountainous areas. First, by using the leading indicators, 14 plain regions, 34 mountain-plain regions, and 27 mountainous regions were defined. Plain regional villages and towns include Yudai, Xingmeng, Xiushan, Yousuo, Sijie, Qianwei, Jiangcheng, Longjie, Liqi, Dajie, Yuxing, Luju, Panxi, and Chunhe. Second, when combined with the field topographic conditions and by means of the auxiliary index “proportion of flatland farmland and construction land” and “proportion of flatland land area,” there were 13 villages and towns with a plain comprehensive index greater than 40, the proportions of flatland cultivated land and construction land were all greater than 50%, and the proportions of flatlands were all greater than 20% among the 32 divided mountain-plain villages and towns. Therefore, these 13 villages and towns were further divided into plain regions, which include Beicheng, Fenghuang, Jiujie, Xiongguan, Huanian, Guishan, and Longquan. Finally, among the 27 mountain regional villages or towns that were divided, Qinglong and Pingdian, whose plain comprehensive index was greater than 15 and the standardized value of flatland (area ≥ 1 km2) was greater than 0.2, were adjusted as mountain-plain regions. Therefore, 27 plain towns, 23 mountain-plain towns and 25 mountainous towns were classified in Yuxi city (Table 6 and Figure 5(d)).

Villages or townsProportion of flatland farmland and construction landProportion of flatland areaStandardized value of flatlandsPlain comprehensive index

Xiaojie (1)0.11240.07880.00008.60
Xiaojie (2)0.40760.21600.050629.32

ClassificationNumber of townsName of towns

Plain towns27Yudai, Xingmeng, Xiushan, Yousuo, Sijie, Qianwei, Jiangcheng, Longjie, Liqi, Dajie, Yuxing, Luju, Panxi, Chunhe, Beicheng, Nagu, Fenglu, Dayingjie, Yanhe, Gaocang, Yangguang, Fenghuang, Jiujie, Xiongguan, Huanian, Guishan, Longquan

Mountain-plain towns23Gaoda, Dong’e, Jiuxi, Tadian, Yangzong, Hexi, Haikou, Anhua, Yangwu, Xiaojie (2), Lishan, Xiaoshiqiao, Liujie, Dianzhong, Ningzhou, Fuliangpeng, Qinglong, Jiucun, Mosha, Pubei, Shuangjiang, Lijiang, Pingdian

Mountain towns25Jianxing, Luohe, Huaxi, Yinyuan, Dalongtan, Gasa, Xinhua, Shijie, Laochang, Chahe, Tongchang, Xiaojie (1), Qinglongchang, Yangjie, Yangchajie, Tonghongdian, Yaojie, Mili, Zhelong, Pingzhang, Lvzhi, Nanuo, Wazhi, Longtan, Shuitang

The spatial distributions of mountainous area types in Yuxi city are shown in Figure 5(d). Plain towns are mainly distributed in the eastern area, i.e., the alluvial plain areas around three lakes (e.g., Fuxian Lake, Xingyun Lake, and Qilu Lake). Only two of 27 plain towns are located in E’shan and Xinping in the central and western parts of Yuxi city while one is in the northern part of Yimen, nine are in Hongta, six are in Tonghai, five are in Jiangchuan, three are in Chengjiang, and one is in Huaning. The mountain-plain towns are mainly distributed on the periphery of plain regions and are located between plain and mountainous regions; the terrain consists of mountains and flatlands alternating with each other. Mountainous towns are mainly concentrated in Xinping in the western Yuxi area, Yimen in the northern area, and E’shan in the central area while Huaning County also has two mountainous regions. The main terrain is mountains with sporadic flatlands.

4. Conclusion and Discussion

(1)Through Python module programming using different windows for automatic extraction of topographic relief and the use of Origin 2019 software for logarithmic equation fitting, it was found that the average topographic relief and window area in Yuxi city had the highest fitting degree under the logarithmic equation and that the determination coefficient was 0.9967. However, the logarithmic fitting equation used by other scholars in the research process was , and the determination coefficients were close to 0.95 [45], 0.96 [34], 0.97, or 0.98 [30, 33]. Therefore, the log 3P1 logarithmic equation was adopted to conduct higher degrees of fitting between the window area and average topographic relief. At the same time, the optimal statistical unit in this study is 56 × 56, or 2.8224 km2, which is similar to the optimal statistical units of 2.43 km2 and 2.62 km2 in southwest China that were calculated by Zhong and Lu [33] and further verified the accuracy of the optimal statistical unit used in this study. However, Zhang and Zhao [32] calculated the area of the optimal statistical unit of topographic relief in the Hengduan Mountain area as 14.98 km2 by using the mean point variation method, which was quite different from the result of this study. First, the DEM data resolution adopted by Zhang et al. was 90 m but was 30 m in this study, which showed great differences in accuracy. Second, the number of experimental units studied by Zhang et al. was only 34 which was quite different from the 198 mobile units used in this study. In addition, Zhang et al. [46] calculated the optimal statistical unit area for Qihe River Basin topographic fluctuations as 0.1296 km2, and Yu et al. [25] calculated the optimal statistical unit area for the Three Gorges Area as 0.15 km2 by using the mean variation analysis method. Based on the above research results, it was determined that the more complex the terrain, the larger is the area of the optimal statistical unit; the flatter the surface, the smaller the area of the optimal statistical unit. Additionally, the larger the research scope, the smaller the area of the optimal statistical unit.(2)In the process of defining mountain types, it was found that it is most scientific to divide by using elevation, topographic relief, and slope through several experiments. Yuxi city is located in Mountains of Southwest China, with an average elevation of approximately 1,800 m. If the mountain type is defined only by using altitude, there are only mountains in our study area. However, during the field investigations, in addition to the low elevation valley areas along the Yuanjiang River, it was found that, in the alluvial plains in the eastern region, especially around the three large lakes, midelevation plain was present, which has elevations from 1,500 to 2,500 m and topographic relief of less than 100 m. In most studies, a slope of 5° is considered to be the critical slope for middle elevation plateaus and middle elevation mountains, but this value is not applicable to the definition of the critical mountain slope in Yuxi city. By comparing the results of multiple experiments with field and remote sensing images and by combining the definition of flatland slope by <Yunnan agricultural geography> and other scholars, it is most reasonable to define the mountain critical slope as 8°, which not only can meet the definition of the mountain type but can also conform to the difference between plain areas and mountainous areas in the field.(3)According to the classification method of Yang and Zhao [39], the classification index for mountainous areas and plain areas was adopted but two improvements were made. First, with increased urbanization levels, flatland farmland with low topographic relief and low slopes is constantly swallowed up by construction land. If only the proportion of flatland cultivated land is used as the measurement index, the results will be biased. The proportion of flatland cultivated land will be reduced to or even less than 50% due to occupation by construction land, which is inconsistent with reality. Second, the number of flatlands as a measurement index is not sufficiently accurate. When dividing plain areas and mountainous areas, it is more important to combine the area and number of flatlands, especially the area. If only the number of flatlands is used as the measurement, there will be large areas of flatlands with small numbers in some villages and towns but the comprehensive value is very low. For example, one flatland in Longjie has an area of 74.26 km2, while six flatlands in Gasa have a total area of only 14.18 km2. According to Yang’s calculation method, the index value of Longjie is only 2.5 while that of Gasa is 15, which is far from the situation. Therefore, it is more reasonable to standardize the flatland area (≥1 km2) and evaluate it in combination with the number of flatlands owned by each township. After this improvement, most of the results are consistent with the field research, but some villages and towns still show exceptions. For example, the proportion of flatland cultivated land in Dong’e town is greater than 35%, the proportion of flatland land is greater than 25%, and there are 6 flatlands whose areas are each greater than 1 km2. Meanwhile, the total area of flatlands is 49.83 km2, but the plain comprehensive index is only 36.21, which can only be classified into a mountain-plain town. Therefore, the accuracy of the dividing index between mountainous areas and plain areas needs to be further improved and discussed.

Data Availability

The dataset used to support the findings of this study is available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This research was financially supported by the National Natural Science Foundation of China under Grant nos. 41971116 and 41807081, the Natural Scientific Research Fund Project of Yunan Province under Grant no. 2017FH001-045, and the National Key Research and Development Project under Grant no. 2016YFC0502406.


  1. M. F. Price, “Forests in sustainable mountain development,” in IUFRO Division 8 Conference on Environmental Forest Science, vol. 54, pp. 39–49, Kyoto University, Kyoto, Japan, 1998. View at: Google Scholar
  2. L. Jansky, “UNU and sustainable mountain development,” Mountain Research and Development, vol. 20, no. 2, pp. 188-189, 2000. View at: Publisher Site | Google Scholar
  3. A. Schild, “The case of the Hindu Kush-Himalayas ICIMOD’s position on climate change and mountain systems,” Mountain Research and Development, vol. 28, no. 3/4, pp. 328–331, 2008. View at: Publisher Site | Google Scholar
  4. L. Jansky, J. D. Ives, K. Furuyashiki, and T. Watanabe, “Global mountain research for sustainable development,” Global Environmental Change, vol. 12, no. 3, pp. 231–239, 2002. View at: Publisher Site | Google Scholar
  5. T. Breu, D. Maselli, and H. Hurni, “Knowledge for sustainable development in the Tajik Pamir mountains,” Mountain Research and Development, vol. 25, no. 2, pp. 139–146, 2005. View at: Publisher Site | Google Scholar
  6. T. Scheurer and S. Jost, “Sustainable mountain development beyond Rio + 20: regional and national perspectives,” Mountain Research and Development, vol. 32, no. 1, pp. 92–95, 2012. View at: Publisher Site | Google Scholar
  7. W. Jin, Y. F. Cui, S. N. Wu, and D. Q. Cheng, “Ecological risk resonance of urbanization and its effect on geohazard disaster: the case of Freetown, Sierra Leone,” Urban Ecosystems Early Access, pp. 1–12, 2020. View at: Publisher Site | Google Scholar
  8. C. S. Zhan, R. X. Zhao, and S. Hu, “Emergy-based sustainability assessment of forest ecosystem with the aid of mountain eco-hydrological model in Huanjiang County, China,” Journal of Cleaner Production, vol. 251, Article ID 119638, p. 14, 2020. View at: Publisher Site | Google Scholar
  9. M. Meybeck, P. Green, and C. Vörösmarty, “A New typology for mountains and other relief classes: an application to global continental water resources and population distribution,” Mountain Research and Development, vol. 21, no. 1, pp. 34–45, 2001. View at: Publisher Site | Google Scholar
  10. Y. Li, X. Yang, H. Cai, L. Xiao, X. Xu, and L. Liu, “Topographical characteristics of agricultural potential productivity during cropland transformation in China,” Sustainability, vol. 7, no. 1, pp. 96–110, 2014. View at: Publisher Site | Google Scholar
  11. G. A. Tang, J. Strobl, J. Y. Gong, M. D. Zhao, and Z. J. Chen, “Evaluation on the accuracy of digital elevation models,” Journal of Geographical Sciences, vol. 11, no. 2, pp. 209–216, 2001. View at: Publisher Site | Google Scholar
  12. R. Hassan, R. Scholes, and N. Ash, Ecosystems and Human Well-Being: Current State and Trends: Findings of the Condition and Trends Working Group, Island Press, Washington, DC, USA, 2005.
  13. B. Xue, J. Z. Li, X. Xiao et al., “Overview of man-land relationship research based on POI data: theory, method and application,” Geography and Geo-Information Science, vol. 35, no. 6, pp. 51–60, 2019, in Chinese. View at: Google Scholar
  14. B. Xue, X. Xiao, J. Z. Li, and X. Xie, “Analysis of spatial economic structure of Northeast China cities based on points of interest big data,” Scientia Geographica Sinica, vol. 40, no. 5, pp. 691–700, 2020, in Chinese. View at: Publisher Site | Google Scholar
  15. J. Yang, W. L. Liu, Y. H. Li, Q. S. Ge, and M. Santos, “Simulating intraurban land use dynamics under multiple scenarios based on fuzzy cellular automata: a case study of Jinzhou District, Dalian,” Complexity, vol. 2018, Article ID 7202985, 17 pages, 2018. View at: Publisher Site | Google Scholar
  16. A. Schild and E. Sharma, “Sustainable mountain development revisited,” Mountain Research and Development, vol. 31, no. 3, pp. 237–241, 2011. View at: Publisher Site | Google Scholar
  17. Y. Yang, Y. Y. Zhao, G. D. Ding, and G. L. Gao, “Effects of urbanization on landscape patterns in a mountainous area: a case study in the Mentougou District, Beijing, China,” Sustainability, vol. 8, no. 11, p. 14, 2016. View at: Google Scholar
  18. T. Kohler and S. W. von Dach, “CDE links regional research and global efforts for sustainable mountain development,” Mountain Research and Development, vol. 35, no. 2, pp. 203–207, 2015. View at: Publisher Site | Google Scholar
  19. M. Kulikov, U. Schickhoff, A. Gröngröft, and P. Borchardt, “Modelling soil erodibility in mountain rangelands of southern Kyrgyzstan,” Pedosphere, vol. 30, no. 4, pp. 443–456, 2020. View at: Publisher Site | Google Scholar
  20. Y. X. Zhang, P. Halder, X. N. Zhang, and M. Qu, “Analyzing the deviation between farmers’ land transfer intention and behavior in China’s impoverished mountainous area: a Logistic-ISM model approach,” Land Use Policy, vol. 94, Article ID 104534, 2020. View at: Publisher Site | Google Scholar
  21. X. Huang, X. Huang, M. Liu, B. Wang, and Y. Zhao, “Spatial-temporal dynamics and driving forces of land development intensity in the Western China from 2000 to 2015,” Chinese Geographical Science, vol. 30, no. 1, pp. 16–29, 2020. View at: Publisher Site | Google Scholar
  22. X. H. Zhong and S. Z. Liu, “Research on the mountain classification in China,” Mountain Research, vol. 32, no. 2, pp. 129–140, 2014. View at: Google Scholar
  23. D. Karátson, T. Telbisz, and B. S. Singer, “Late-stage volcano geomorphic evolution of the Pleistocene San Francisco Mountain, Arizona (USA), based on high-resolution DEM analysis and 40Ar/39Ar chronology,” Bulletin of Volcanology, vol. 72, no. 7, pp. 833–846, 2010. View at: Publisher Site | Google Scholar
  24. W. Cheng, N. Wang, M. Zhao, and S. Zhao, “Relative tectonics and debris flow hazards in the Beijing mountain area from DEM-derived geomorphic indices and drainage analysis,” Geomorphology, vol. 257, pp. 134–142, 2016. View at: Publisher Site | Google Scholar
  25. H. Yu, Y. Luo, S. Liu, Y. Wang, Y. Yang, and W. D. Liu, “The influences of topographic relief on spatial distribution of mountain settlements in Three Gorges Area,” Environmental Earth Sciences, vol. 74, no. 5, pp. 4335–4344, 2015. View at: Publisher Site | Google Scholar
  26. C. Liu, W. Sun, and H. Wu, “Determination of complexity factor and its relationship with accuracy of representation for DEM terrain,” Geo-spatial Information Science, vol. 13, no. 4, pp. 249–256, 2010. View at: Publisher Site | Google Scholar
  27. O. D. A. Prima, A. Echigo, R. Yokoyama, and T. Yoshida, “Supervised landform classification of Northeast Honshu from DEM-derived thematic maps,” Geomorphology, vol. 78, no. 3-4, pp. 373–386, 2006. View at: Publisher Site | Google Scholar
  28. L. Wang and X. J. Tong, “Analysis on relief amplitude based on change point method,” Geography and Geo-Information Science, vol. 23, no. 6, pp. 65–67, 2007, in Chinese. View at: Google Scholar
  29. S. B. Liu, Y. P. Cui, F. X. Lu, and X. Y. Qin, “The analysis of relief amplitude in Qinghai Province based on GIS,” in Proceedings of the International Conference on Geoinformatics, IEEE, Kaifeng, China, June 2013. View at: Publisher Site | Google Scholar
  30. J. J. Zhang, Y. H. Li, and L. Q. Zhu, “Relief amplitude based on county units in West Henan mountain area and its correlation with distribution of population and economic activities,” Areal Research and Development, vol. 38, no. 2, pp. 55–60, 2019, in Chinese. View at: Google Scholar
  31. J. C. Zhang and W. Z. Zhou, “Spatial autocorrelation between topographic relief and population/economy in Sichuan Province,” Bulletin of Soil and Water Conservation, vol. 39, no. 1, pp. 250–257, 2019, in Chinese. View at: Publisher Site | Google Scholar
  32. Y. Zhang and Y. L. Zhao, “Mountainous area divisions at county level in Hengduan Mountains based on DEM,” Journal of Guizhou Normal University (Natural Sciences), vol. 34, no. 6, pp. 8–14, 2016, in Chinese. View at: Google Scholar
  33. J. Zhong and T. Lu, “Optimal statistical unit for relief amplitude in Southwestern China,” Bulletin of Soil and Water Conservation, vol. 38, no. 1, pp. 175–181, 2018, in Chinese. View at: Publisher Site | Google Scholar
  34. H. H. Han, T. Gao, H. Yi et al., “Extraction of relief amplitude based on change point method: a case study on the Tibetan Plateau,” Scientia Geographica Sinica, vol. 32, no. 1, pp. 101–104, 2012, in Chinese. View at: Google Scholar
  35. X. R. Zhang and K. Dong, “Neighborhood analysis-based calculation and analysis of multi-scales relief amplitude,” Advanced Materials Research, vol. 468–471, pp. 2086–2089, 2012. View at: Publisher Site | Google Scholar
  36. Z. Ma, X. Chen, and H. Chen, “Multi-scale spatial patterns and influencing factors of rural poverty: a case study in the Liupan Mountain Region, Gansu Province, China,” Chinese Geographical Science, vol. 28, no. 2, pp. 296–312, 2018. View at: Publisher Site | Google Scholar
  37. Y. B. Li, H. Chen, and G. J. Luo, “Spatial distribution of different scale flatland in Guizhou Province, China,” Scientia Geographica Sinica, vol. 39, no. 11, pp. 1830–1840, 2014, in Chinese. View at: Google Scholar
  38. Y.-p. Fang and B. Ying, “Spatial distribution of mountainous regions and classifications of economic development in China,” Journal of Mountain Science, vol. 13, no. 6, pp. 1120–1138, 2016. View at: Publisher Site | Google Scholar
  39. Z. S. Yang and Q. G. Zhao, “Study on dividing flatland county, semi-mountainous & semiflatland county and mountainous county in Yunnan Province based on the second national land survey,” Journal of Natural Resources, vol. 29, no. 4, pp. 564–574, 2014, in Chinese. View at: Google Scholar
  40. Y. L. Zhao and X. B. Li, “Spatial correlation between type of mountain area and land use degree in Guizhou Province, China,” Sustainability, vol. 8, no. 9, p. 849, 2016. View at: Publisher Site | Google Scholar
  41. M. Schuler, E. Stucki, and O. Roque, Mountain Areas in Europe: Analysis of Mountain Areas in EU Member States, Acceding and Other European Countries, European Commission, Brussels, Belgium, 2004.
  42. B. Xue, L. M. Zhang, Y. Geng, B. Mitchell, and W. X. Ren, “Extended land-use coding system and its application in urban brownfield redevelopment: case study of Tiexi District in Shenyang, China,” Journal of Urban Planning and Development, vol. 142, no. 3, Article ID 05015014, 2016, in Chinese. View at: Publisher Site | Google Scholar
  43. A. D. Guo, Y. Q. Zhang, Q. Hao, and J. H. C. Xia, “Monitoring and simulation of dynamic spatiotemporal land use/cover changes,” Complexity, vol. 2020, Article ID 3547323, 12 pages, 2020. View at: Publisher Site | Google Scholar
  44. Editing Group of Yunnan Agricultural Geography, Yunnan Agricultural Geography, Yunnan People’s Publishing House, Kunming, China, 1981, in Chinese.
  45. R. H. Wang, J. C. Yang, S. W. Zhang, L. P. Chang, L. M. Pu, and K. Bu, “Analysis on the relief amplitude in Northeast China based on ASTER GDEM and mean change point method,” Journal of Arid Land Resources and Environment, vol. 30, no. 6, pp. 49–54, 2016, in Chinese. View at: Publisher Site | Google Scholar
  46. J. Zhang, W. Zhu, F. Zhao et al., “Spatial variations of terrain and their impacts on landscape patterns in the transition zone from mountains to plains-a case study of Qihe River Basin in the Taihang Mountains,” Science China Earth Sciences, vol. 61, no. 4, pp. 450–461, 2018. View at: Publisher Site | Google Scholar

Copyright © 2020 Li Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.