Abstract

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 (http://www.gscloud.cn). 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 (http://www.resdc.cn). 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).

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.

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.

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.

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)).

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.

Acknowledgments

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.