**GIS-Based Synthetic Measurement of Sustainable Development in Loess Plateau Ecologically Fragile Area—Case of Qingyang, China**

#### **Chenyu Lu, Chunjuan Wang, Weili Zhu, Hengji Li, Yongjin Li and Chengpeng Lu**

**Abstract:** Synthetic measurement of regional sustainable development has been one of the key issues in the research field of sustainability. In this paper, Qingyang City located in the Loess Plateau ecologically fragile area of Northwest China is used for a case study, and the present study aims to investigate the degree of sustainable development by conducting temporal- and spatial-scale based analysis, with the assessment index system, assessment model and GIS approach well integrated. The results show that the development pattern of Qingyang generally fits the mode of unsustainable development, even in the presence of certain levels of spatial differences. The sustainable development state in ecologically fragile area of China's Loess Plateau is non-optimistic, which is an uncoordinated status among subsystems of regional sustainable development. Although the level and tendency of regional sustainable development keeps increasing, such enhancement is abnormal. With the rapid deterioration of environmental and natural resources, their inhibitory effect on the economy and society would expand, eventually leading to the slow development rate or the recession of the entire system. The only solution is to change the traditional mode of economic development, to follow the guide of ecological economic conception so that the goal of achieving regional sustainable development strategies could be met ultimately. Meanwhile, the characteristics of different regions should be taken into account in order to achieve optimal spatial structure.

Reprinted from *Sustainability*. Cite as: Lu, C.; Wang, C.; Zhu, W.; Li, H.; Li, Y.; Lu, C. GIS-Based Synthetic Measurement of Sustainable Development in Loess Plateau Ecologically Fragile Areaü Case of Qingyang, China. *Sustainability* **2015**, *7*, 1576-1594.

#### **1. Introduction**

Ever since the Industrial Revolution took off in the 18th century, humanity's infinite material and spiritual needs have been satisfied through continuous expansion of traditional economic scale, characterized by "mass production, mass consumption and mass abandonment" [1]. Under such circumstances, a lot of problems have appeared to exert negative feedback on the normal functioning of social and economic activities [2]. Therefore, how to achieve sustainable development has become a common life aspiration of all human beings [3]. In China, environmental problems have become increasingly prominent, and the worsening conditions of the environment haven't been well controlled [4]. Therefore, sustainable development has become the basic guiding philosophy and national development strategy in China for coordinated development in social, economic and environmental aspects [5].

Studies of sustainable development in other countries can be divided into two phases: phase I (1962–1992) represents the phase when the concept of sustainable development was proposed and its associated theories explored. United Nations Conference on Environment and Development was held in Rio de Janeiro in 1992, after which the implementation phase of sustainable development strategy was officially initiated. Since the concept of sustainable development has been put forward, how to define this concept and how to link this abstract concept with cognitive and practical methods become two major challenges [6]. Meanwhile, the word sustainability with quantitative value, together with its related theoretical and practical studies, have become key components of study on sustainable development [7]. Therefore, the assessment of sustainable development has become a hot and central research topic in the field of sustainable development [8], and a series of theoretical and practical studies related to sustainability assessment have been done. In general, two types of assessment methods have been developed. The first type uses a single or composite index, whereas the second type applies multi-index system. Also, both types are still in the exploratory stage. Overall, the specific relationship between indexes and sustainable development goals, the determination of index weight and threshold, and the development of an integrated assessment method remains the core and challenge of any assessment system [9]. In particular, how to establish a reasonable sustainability assessment system with an integrated assessment model, and how to study the contributions of human activities to the realization of the goal of sustainable development represent key and hot topics of current sustainable development research [10].

In China, studies on sustainable development have started to track the trajectory of research on global sustainable development since the mid-1980s [11,12], which could be further divided into two parts, namely, theoretical studies and practical applications. Theoretical studies mostly focus on aspects such as the establishment of sustainability assessment system [13–15], the analysis of coordination development in PRED (population, resources, environment and development) system [16–18], and the assessment and determination of sustainability [19–21]. By contrast, practical studies mostly focus on aspects such as regional coordination development in PRED system [22–24], and the evolution of regional sustainability over time [25–27]. Although the work of some scholars involves the spatial distribution aspect of sustainability [28,29], most work is simply static description with the lack of time content. Currently, research on China's sustainable development is developing rapidly with a consistent deepening trend. As to regional practice, accompanying with the construction advance of national main functional zones, particularly with increasing environmental constraints that guide China's future development, many counties and districts have put an emphasis on promoting the process of sustainable development, and quantitative characterization of sustainable development level and its development process is urgently needed [30]. In order to promote the process of sustainable development towards a new height and to better draw lessons from experience, research on sustainable development field has displayed a gradual shift, namely, from a qualitative approach towards a quantitative approach, with an emphasis on the combination of such shift when factors of regional geography and environment are explicitly considered [31]. Therefore, such process reflects a shift from traditional conceptual models towards specific regional models [32], and a shift from static pattern study towards dynamic process study [33].

Overall, relevant studies have achieved a lot, although some problems and deficiencies still exist. First, studies that focus on the integration of spatial and temporal aspects of sustainable development are still lacking. Also, among existing research results, studies that have taken spatial and temporal dimensions into account simultaneously in order to conduct an integrated study on the evolution of sustainable development are extremely rare. By contrast, most current studies are testing the evolution patterns of sustainable development at temporal scale, although research on sustainable development should ultimately focus on specific geographic regions, or, in other words, research on the evolution of sustainable development at spatial scale is equally important. Second, many problems have been found regarding the assessment system of sustainable development, including mismatches between theory and practice, huge number of index systems, overemphasis on economic assessment, overlook of resource and environment assessment, and the failure of index selection process to truly implement the concept of sustainable development. To a certain extent, the present paper attempts to solve some of these problems by providing some complementary work.

In the present paper, Qingyang City is used as the case study region. Qingyang is located in Loess Plateau zone of eastern Gansu Province, and represents a typical gully region of China's Loess Plateau. The administrative division in Qingyang City includes eight sub-regions, seven counties and one district (Figure 1). It belongs to semiarid region due to the lack of surface water. Thanks to natural conditions and human economic activities, the ecological environment is fragile in this region [34]. In terms of economic structure, agriculture has played a dominant role in Qingyang. In recent years, as Qingyang becomes an important component of Ordos-based national energy industrial base and the energy and chemical industry base of eastern Gansu Province, the contribution of industry to national economy has exceeded that of agriculture. Meanwhile, traditional agriculture-based industrial structure has been broken, leading to the transition between agriculture and industry, as well as the rapid development of local economic structure [35]. However, with the deepening of the transition between the agriculture and industry, contradictions among aspects such as regional population, economy, society, resources and environment become increasingly acute, and man-land relationship has deteriorated rapidly in Qingyang. Also, accompanied with the emergence of a series of ecological and environmental problems, many issues begin to arise from both natural and anthropogenic disturbance. It is not surprising that both the complexity of the integrated system and regional sustainable development have met severe challenges [36]. Worse still, constraints imposed by resources and environment on economic and social development become more and more severe, and contradictions between local resource development and environmental pressure become more and more acute. Equally, the capacity for sustainable development is worrisome. In particular, the lack of water resources and quality in Qingyang has been worsened by a groundwater pollution problem due to oil extraction, which has caused increasingly grim situation of water resource, and seriously restricted economic, social and sustainable development. In fact, a series of core issues related to sustainable development, such as the transition of industry structure, the deterioration of ecological environment, and the transformation of social structure, has been commonly found in the ecologically fragile area of China's Loess Plateau [37]. Therefore, we must stand on the height of sustainable development, and wisely handle the coordinated and balanced development in terms of regional population, economy, society, resources and environment.

**Figure 1.** The administrative map in Qingyang City.

In this paper, the assessment index system and the integrated assessment models with temporal and spatial scale explicitly considered are established, with existing research results and regional characteristics as the basis. Also, the application of GIS technology allows a successful integration of both spatial and temporal dimensions so that a study on synthetic measurement of sustainable development evolution in ecologically fragile area of China's Loess Plateau could be facilitated, the process of two-dimensional analysis and integration could be carried out, and deficiencies in the existing research could be compensated. Overall, the present paper has strong theoretical and practical significance, and could contribute to the understanding, revealing and solving of a series of problems that exist during the process of sustainable development, just as what Qingyang is experiencing now. In this sense, this paper could provide theoretical support and decision-making basis for achieving sustainable development strategies.

#### **2. Methodology**

#### *2.1. The Establishment of the Index System*

The assessment index system of sustainable development is the basis for an integrated assessment of regional sustainable development. As it concerns the division of regional sustainable development, it is generally considered an improvement based on PRED system. According to the conception of PRED, the system is divided into two subsystems (*i.e.*, the economy and the society) in the present study so that the characteristic of regional sustainable development could be reflected more clearly. In a word, based on geographical features as well as actual situations, draw on existing achievements, and followed the integrated, systematic, representative, concise, accessible and practical principles [38,39], the present paper chooses corresponding indexes from five aspects, namely, the population, economy, society, environment and resources, and then establishes the four-category index system that reflects the sustainability level of ecological fragile area in China's Loess Plateau.

According to the criterion layers of scale, structure and health level of the population, 5 indexes are selected for the population system. In addition to these criterion layers, the quality of the population is always counted as the criterion layer, which could reflect the educational status of local people. Since the educational status of local people is reflected in a social system, it is not selected for the population system. Since all the selected indexes are reasonable, the criterion layers could be embodied by index layers, and the characteristics of the population system could be reflected accurately by these criterion layers and index layers.

According to the criterion layers of scale, structure, efficiency and openness of economy, 9 indexes are selected for the economic system. Since 4 criterion layers and 9 index layers are commonly selected for the economic system in the literature, they are both representative and reasonable. The criterion layers can be embodied by index layers relatively. Overall, the characteristics of the economic system could be reflected accurately by these criterion layers and index layers.

According to the criterion layers of quality of people's lives, social development level, and social activity, 9 indexes are selected for the social system. Since 3 criterion layers are commonly selected for the social system in the literature, they are both representative and reasonable. Also, the 9 selected index layers are reasonable. Due to the lack of sufficient data, several indexes are not selected as the criterion layers of social development level (e.g., the percentage of social insurance coverage). Overall, the characteristics of society system can be reflected clearly by these criterion layers and index layers.

According to the criterion layers of natural resources and social resources, 7 indexes are selected for the resource system. Two criterion layers are selected based on the basis of resource classification, and they are both representative and reasonable. However, the index layers are not ideal for the resource system. Due to the lack of sufficient data, several significant indexes are not selected (e.g., the indexes of water resource and energy consumption). Overall, the characteristics of the resource system could be reflected accurately by these criterion layers and index layers.

According to the criterion layers of environmental pollution and protection, 8 indexes are selected for the environmental system. Since 2 criterion layers are commonly used in the literature, they are representative and reasonable. Also, the selected 8 index layers are reasonable. Due to the lack of sufficient data, several indexes are not selected as the criterion layers of environmental protection (e.g., the harnessing percentage of soil and water loss). Overall, the characteristics of the environmental system could be reflected accurately by these criterion layers and index layers.

Meanwhile, index attributes are categorized. If the larger an index value is, regional sustainable development is more likely, such index is then named a positive index, notated as "+". By contrast, if the larger an index value is, regional sustainable development is more unlikely, such index is then named a negative index, notated as "í" (Table 1).


**Table 1.** The assessment index system of sustainable development.

In order to ensure that the data could reflect actual situations of the studied area as much as possible, and to eliminate possible effects of different index dimensions on assessment results, extreme values are standardized so that assessment indexes are dimensionless [40]. The calculation formula is:

Positive index: the larger the index value is, the better it is.

$$Z\_y = \frac{X\_y - X\_{\text{min}}}{X\_{\text{max}} - X\_{\text{min}}} \tag{1}$$

Negative index: the smaller the index value is, the better it is.

$$Z\_y = \frac{X\_{\text{max}} - X\_y}{X\_{\text{max}} - X\_{\text{min}}} \tag{2}$$

In this formula, is the standardized value of ; is the actual value of the administrative unit *i* and indexical stage *j*; and represent the maximum and minimum value, respectively, of the same index at different administrative unit and stage. *ij Z Xij Xij X* max *X* min

#### *2.2. Integrated Assessment Model*

There are lots of models for the assessment of regional sustainable development, and they could be generally divided into 2 types. The first type puts an emphasis on sustainability measurement, which is commonly used to evaluate the level or the capacity of regional sustainable development. The second type puts an emphasis on coordination measurement, which is used to evaluate the coordination development status among different indexes or different systems. Generally speaking, both types tend to put an emphasis on certain aspects of regional sustainable development, and thus are unilateral or non-comprehensive. In this study, an integrated assessment model based on these 2 types is used for a comprehensive evaluation of regional sustainable development. One obvious advantage of this model is that it can be used for a comprehensive assessment of regional sustainable development based on the aspects of sustainability and coordination measurement, and the results of comprehensive measurement could be quantified. Theoretically and empirically, this model is more reasonable than other alternatives. Therefore, this integrated assessment model with 4 components is selected for the present study.

#### 2.2.1. Index of Sustainable Development

Principal component analysis is applied to eliminate possible redundant information among various indexes, to extract important factors, and to calculate the contribution of each index. Then, by using a weighted linear model, the composite score, namely, the index of sustainable development level in each subsystem is calculated [41].

$$F\_m = \sum\_{j=1}^n W\_j X\_j \tag{3}$$

In this formula, *Fm* represents the index of sustainable development level in each subsystem, *m* represents the number of subsystems; *n* represents the number of principal components in each subsystem; *Wj* represents the contribution level of each principal component in each subsystem; *Xj* represents the score of all principal components in each subsystem.

The integrated index value of sustainable development level could reflect the integrated level of regional sustainable development. Based on obtained index values of sustainable development in various subsystems, weighted linear models are used for measurement purpose [41]:

$$SDI = \sum\_{j=1}^{m} \alpha\_j F\_m \tag{4}$$

In this formula, *SDI* represents the integrated index of regional sustainable development; *Fm* represents sustainable development level of each subsystem *Įj* represents the weighted value of each subsystem. In the present study, the coordination degree needs to be measured. Since the conception of coordination degree is based on the equal importance of each subsystem, here it is assumed that all five subsystems are equally important, and each subsystem is assigned the value 0.2.

#### 2.2.2. Sustainability Degree

Sustainability degree is the criterion used to measure regional sustainable development and the sustainability of each subsystem. Usually, the developmental rate of each subsystem is used as the calculation criterion [42]:

$$ISI = \frac{1}{T} \sum\_{i=0}^{T-1} S(t - i) \tag{5}$$

In this formula, *SI* represents the sustainability of regional sustainable development; *S(t-i)* represents the developmental rate of the system during the time span of *t-i*~; ; *T* represents the reference period. *i T* 0,1, , 1

#### 2.2.3. Coordination Degree

Since the developmental process of a system is dynamic, the coordination degree is used to quantify the coordinated development degree of interactive factors within or among systems. Therefore, the coordinated development degree among regional systems could be quantified and determined through the calculation of coordination degree. In the present paper, principal component analysis, simulate regression model and the membership function of fuzzy mathematics are used to assess the coordination degree among subsystems. Based on the development index of each subsystem, regression curve fitting process could provide the predictive value of sustainable development index of each subsystem, namely, the coordination value, and the coordination degree between subsystems could be calculated then using both the coordination value and the membership function of fuzzy mathematics [43].

$$\ln(i/j) = \exp\left\{-\frac{\left(F\_i - \bar{F}\right)^2}{S^2}\right\} \tag{6}$$

In this formula, is the state coordination degree of system relative to system ; is the development index of system (the actual value); is the predictive value of best development degree of system relative to system (the coordinated value); is the variance of system . *ui j* (/ ) *<sup>i</sup> <sup>j</sup> Fi i* ' *F <sup>j</sup> i* <sup>2</sup> *S i*

The coordination degree between two systems, *U(i*, *j)*, could be calculated via coordination degree, with the specific formula shown as following [43]:

$$U(i,j) = \frac{\min\left\{\mu(i'/j), \mu(j/i)\right\}}{\max\left\{\mu(i'/j), \mu(j/i)\right\}}\tag{7}$$

In this formula, represents the coordination degree between system *i* and system *j*; represents the state coordination degree of system *i* relative to system *j*; represents the state coordination degree of system *j* relative to system *i*. , and the more it is close to 1, it is more coordinated between system *i* and system *j*. Otherwise, it is more non-coordinated. *Ui j* (, ) *ui j* (/ ) *uji* ( /) *Ui j* ( , ) 0,1 > @

The coordination degree among more than three systems could be calculated as shown below [43]:

$$U\_k(1,2,...,k) = \frac{\sum\_{l=1}^k U(i \wedge \overline{i}\_{k-1}) \times U\_{k-1}(\overline{i}\_{k-1})}{\sum\_{l=1}^k U\_{k-1}(\overline{i}\_{k-1})}, k = 3,4,...,m\tag{8}$$

In this formula, is the number of subsystem; is the set of any *kí1* subsystem except for system *i*; is the state coordination degree imposed by subsystem *i* on any *kí1* subsystem; is the coordination degree of any *kí1* subsystem except for system *i*. *m <sup>k</sup>* <sup>1</sup> *i* <sup>1</sup> (/ ) *Ui ik* 1 1 ( ) *U i k k*

According to the obtained formula used for calculating the coordination degree, when the actual value is more close to the coordinated value, the larger the coordination degree is, and the entire system is with higher coordination degree. Otherwise, the system is with lower coordination degree. Also, with reference to relevant literature [40], the coordination degree of a system could be categorized as four levels, namely, extremely uncoordinated, uncoordinated, barely coordinated and coordinated (Table 2).

**Table 2.** Metrics for the coordination degree.


#### 2.2.4. Tendency Degree

Tendency degree refers to the maintaining or sustainable ability of a regional system, and is commonly used to measure and determine whether a regional system is developed towards a healthy direction or not. Therefore, it is an integrated indicator and an inherent union of system development, continuity and coordination [42]:

$$S = \alpha\_1 \Box SDI + \alpha\_2 \Box Sl + \alpha\_3 \Box U \tag{9}$$

In this formula, *S* is the tendency degree of sustainable development; *SDI*, *SI* and *U* represent the integrated index of sustainable development degree, sustainability degree and coordination degree. Since the present paper assumes that all three are equally important, *Į*1, *Į*2 and *Į*3 is with same weighted value, and the weighted average method is used to calculate the weighted average with *Į*1 = *Į*2 = *Į*3 = 1/3.

#### *2.3. A Model Based Approach for Analysis of Spatial Structure*

The regional centroid refers to the equilibrium point of an element's torque on its space plane. In practice, with the assumption that a specific region consists of *n* space units, the centroid coordinate of space unit *i* is (*xi*, *yi*), is the element value of space unit *i*, and *x*, *y* represent the spatial coordinate of such element within the entire region. The calculation method is as shown following: *Mi*

$$\infty = \left(\sum\_{l=1}^{n} M\_{l}\chi\_{l}\right) / \left(\sum\_{l=1}^{n} M\_{l}\right) \tag{10}$$

$$\mathcal{Y} = (\sum\_{i=1}^{n} M\_i \mathcal{Y}\_i) / (\sum\_{i=1}^{n} M\_i) \tag{11}$$

In this study, represents element attributes of county *i*, *xi* and *yi* represents the spatial coordinate of the administrative centroid of county; and *i*, *x* and *y* represents the centroid coordinate with relevant attributes. When the centroid attributes of an element are moving through time, the moving direction will point at high-density area, and the deviation distance will indicate the degree of un-equilibrium [44,45]. Also, it can be seen from centroid formula that if the geological locations of every site remain unchanged (*i.e.*, such consistence is assumed unchanged throughout the research period), then factors that directly affect the centroid would be the speed and degree of the centroid attributes. Since the development and change of such attributes would affect regional change of centroid, the latter could clearly reflect the developmental trajectories and spatial differences of such attributes. *Mi*

#### **3. Results and Discussion**

The chosen assessment system and integrated assessment model is used to thoroughly access and analyze the sustainable development status of each county in Qingyang, which, combines with the evolution of spatial structure analysis, reveals the spatial and temporal evolution characteristics and patterns of its regional sustainable development. This present study selects five years, namely, 2000, 2003, 2006, 2009 and 2012 for further analysis. The origin of data used in this study is some local statistical information, including statistical yearbook, statistic bulletin and government report provided by local government.

#### *3.1. Integrated Assessment and Analysis*

As to the aspect of integrated indexes that reflect regional sustainable development, each county in Qingyang has shown a gradual rising trend in recent decade since 2000 (Table 3). Specifically, the rate of increase is relatively fast during the early stage with a high magnitude. However, the climbing space is constrained during the late stage, and the rate of increase becomes slow. This is because the speed of current economic and social development is relatively fast. Although both resources and environmental systems have displayed some degradation trend, the exponential growth rate of sustainable development indexes regarding economic and social systems temporarily exceeds the decreasing rate of sustainable development indexes regarding resources and environmental systems, thus causes such observed phenomena. However, the increase of index values is just temporary with crises hidden behind, which indicates that the increase of integrated degree of regional sustainable development in Qingyang is an abnormal pattern because such development is at the expense of environmental resources.


**Table 3.** The integrated index of sustainable development level.

In terms of the sustainability degree of regional sustainable development, each county in Qingyang has displayed a first increasing and then decreasing trend (Table 4). Before 2009, the overall sustainability degree of every county in Qingyang displayed a rising pattern with an increase of amplitude. After 2009, however, the sustainability degree of every county in Qingyang showed certain level of decrease with a decrease of amplitude. This indicates that during the early stage every county in Qingyang continued to show a certain level of non-deceleration development that is both stable and sustainable. However, during the late stage, the status of sustainable development has been shifted towards the status of unsustainable development due to increasing pressure imposed by rapid economic and social development. Overall, the development in Qingyang is unsustainable.


**Table 4.** The sustainability degree.

As to the coordination aspect of regional sustainable development, although slight fluctuations have been observed in recent decade since 2000, each county in Qingyang has shown varying levels of a declining trend (Table 5). Before 2009, the overall coordination degree of the entire area showed a significant decreasing trend. Specifically, County Huachi and Heshui was changed from the status of coordinated to barely coordinated, whereas the rest was changed from the status of barely coordinated to uncoordinated. Obviously, rapid economic growth has brought enormous pressure on natural resources and environment, and the developmental direction of economic society was opposite to that of natural resources and environment. As a result, coordination degree has been weakened in general, and the coordinated development between some subsystems was gradually deteriorated in particular. Since 2009, the decreasing trend of coordination degree among counties has been constrained to a certain degree, and some counties even have started to show an increasing trend. However, such situation still belongs to the uncoordinated status. This is because in recent decade, Qingyang City government has begun to implement the circular economy as a development strategy. To some extent, such practice has successfully curbed the degradation trend of resources and environmental systems, and thus contributed to a much more coordinated development between subsystems, although it is far less than enough to completely reverse the worsening trend. Overall, the uncoordinated development trend found among Qingyang City's population, economic, social, resources and environmental systems has imposed severe challenges to its regional sustainable development.


**Table 5.** The coordination degree.

As far as the tendency degree of regional sustainable development is concerned, each county in Qingyang has shown a gradual upward trend in recent decade since 2000 (Table 6). Specifically, the rate of increase is relatively fast during the early stage with a high magnitude. However, the rate of increase becomes slow during the late stage. This is because economic and social development is relatively fast, and thus plays an important supportive role in promoting regional sustainable development. Although both resources and environmental systems have displayed some degradation trend, the rapid development of the economy and society has masked the reality, resulting in a gradual increase in the integrated indexes that reflects the status of regional sustainable development. However, the increase of such index values is just temporary, and the increasing trend could not be maintained persistently over time. Actually, the rate of increase becomes slow during the early stage, which has fully demonstrated this point. On the surface, Qingyang City is developing towards sustainable development status, although the truth is that a lot of crises are hidden behind, or, in other words, the current development is by no means sustainable.


**Table 6.** The tendency degree.

#### *3.2. Spatial Structure Analysis*

On the basis of geographic information system and ArcGIS software, the present paper applies centroid displacement method to analyze the evolution of spatial structure regarding regional sustainable development in Qingyang, with spatial gaps of regional sustainable development visualized and quantified. The detailed procedures include: (1) the establishment of spatial databases. On the basis of ArcGIS software, Qingyang administrative map is digitalized (including administrative boundaries and the administrative center), and projection transformation is converted in order to determine the spatial plane coordinate of administrative center for each district. Meanwhile, the attributes of all elements in each district through years are imported; (2) By using the centroid calculation formula and on the basis of relevant element attributes, the spatial centroid of each element through years in Qingyang is calculated. Based on obtained centroid results, the trajectory of the moving centroid for each element is drawn in order to help analyze temporal and spatial evolution patterns.

In Qingyang region, the centroid distribution width of sustainable development integrated index is relatively large with a high magnitude of fluctuations over years (Figure 2). Between 2000 and 2003, the centroid was shifted towards the northeast direction with a high magnitude, and the displacement distance was around 50 km. This is due to the rapid growth rate of regional sustainable development in Huachi and Heshui counties located in the northeastern part of Qingyang. Also, both counties are known for their relatively high level of sustainable development. Between 2003 and 2006, the centroid was shifted towards the southwest direction with a high magnitude. During this period, the growth rate of regional sustainable development in Xifeng, Qingcheng and Zhengyuan is relatively high, all of which are located in the southwestern part of Qingyang with relatively high level of sustainable development. Between 2006 and 2009, the centroid displacement distance was relatively short (e.g., several tens of meters with a high level of overlap). During this period, the spatial distribution pattern of every county and district regarding regional sustainable development did not change significantly. Between 2009 and 2012, the centroid was shifted slightly towards the southeastern direction. This is due to the rapid growth rate of sustainable development in Ning and Zhengning located in the southeastern area of Qingyang.

The centroid of sustainability degree regarding regional sustainable development in Qingyang is with a certain level of fluctuations but limited fluctuation range (Figure 3). All the centroids of sustainability degree are distributed in the area of Qingcheng County. Between 2003 and 2006, the centroid was shifted towards southwest direction with relatively large level of displacement. This is mainly due to relatively high levels of sustainability degree in Xifeng, Qingcheng and Zhengyuan counties, all of which are located in the southwestern part of Qingyang. Between 2006 and 2009, the centroid was shifted slightly towards the north direction, indicating that the rate of increase regarding the sustainability degree is slightly higher in the northern area than that of the southern area, with the latter displaying smaller magnitudes. Between 2009 and 2012, the centroid was shifted towards the southeast direction to some extent. During this period, the sustainability degree of every country and district in Qingyang showed certain level of decrease. However, since the level of decreasing amplitude in Ning and Zhengning in the southeastern parts was relatively small, the overall sustainability degree remained at a relatively large level, and thus the centroid was shifted towards the southeast direction.

**Figure 2.** The centroid distribution of integrated index of sustainable development level.

The centroid distribution range of sustainable development coordination degree in Qingyang region is relatively narrow with small magnitude of fluctuations over time. All the centroids of coordination degree are distributed in the area of Qingcheng County. Obviously, the spatial difference regarding the change of coordination degree is not very large (Figure 3). Between 2000 and 2003, the centroid was slightly shifted towards the southwest direction. Between 2003 and 2006, or between 2006 and 2009, the centroid was shifted towards the northeast direction with relatively small magnitude of displacement. Between 2009 and 2012, the centroid was shifted again towards the southwest direction. Overall, in the recent decade, the spatial distribution pattern regarding the coordination degree in every county and district of Qingyang did not change significantly, and the subtle differences could only be detected along the southwest-northeast direction.

**Figure 3.** The centroid distribution of sustainability degree, coordination degree, and tendency degree.

The centroid of the tendency degree regarding sustainable development in Qingyang displays certain levels of fluctuations with relatively small magnitudes (Figure 3). All the centroids of tendency degree are distributed in the area of Qingcheng County. Between 2003 and 2006, the centroid was shifted towards the southwest direction with relatively large magnitude of displacement. This is due to the rapid growth rate of tendency degree in Xifeng, Qingcheng and Zhengyuan, all of which are located in the southwest area of Qingyang and with high levels of tendency degree. Between 2006 and 2009, the centroid was slightly shifted towards the north direction, indicating that the rate of increase regarding tendency degree is slightly higher in the northern area than that of the southern area, with the latter displaying smaller magnitudes. Between 2009 and 2012, the centroid was shifted towards the southeast direction. This is due to the rapid growth rate of tendency degree in Ning and Zhengning located in the southeast area of Qingyang, and both shared relatively large magnitude of increase.

#### **4. Conclusions**

Since 2000, the development pattern of Qingyang has fitted the mode of unsustainable development in the presence of certain levels of spatial differences. Under current development scenarios with rapid economic and social development, the natural resources and environment are gradually degrading. Because such development is at the expense of natural resources and environment, the coordination degree among population, economy, society, resources and environment is damaged, and the current state of development is basically uncoordinated. Currently, although the level of regional sustainable development and associated tendency degree remains arising, such a trend is solely temporary with severe hidden crises, and the tendency degree begins to show a declining trend. In the long run, such pattern would lead to the overall decrease of all indexes related to regional sustainable development, and Qingyang would not be able to achieve its healthy, long-term sustainable development. In terms of spatial structure, the spatial differences of integrated indexes in terms of regional sustainable development represent the largest ones, followed by the sustainability and tendency degree, and the spatial differences of the coordination degree represent the smallest ones. Also, the spatial differences along the northeast-southwest direction are predominant.

The challenges that Qingyang faces during the critical period of its development are typically found in ecologically fragile area of China's Loess Plateau, known for its non-optimistic sustainable development reality. In such area, water scarcity is a big problem, the ecological environment is fragile, and the poorest or most backward districts could be easily found. To get rid of poverty, local government often emphasizes the rapid growth rate of economy, and is desperate to promote extensive economic growth at the expense of natural resources and environment, causing uncoordinated status among subsystems of regional sustainable development, which would further restrict the extent of regional sustainable development. Meanwhile, with the accelerated pace of industrialization process, the economy has been spurred in such area, which might counterbalance the restraining effect imposed by environmental degradation on regional sustainable development, and enable the level of regional sustainable development to keep increasing. However, such enhancement is abnormal. With the rapid deterioration of the environment and natural resources, their inhibitory effect on economy and society would be expanded, eventually leading to the slow development rate, or the recession, of the entire system.

For the ecologically fragile area of China's Loess Plateau, the only solution is to change the traditional mode of economic development, to follow the guide of ecological economic conception for regional development, to walk the road of circular economy in order to make the population, economy, society, resources and environment systems perfectly integrated, to realize the optimal allocation of regional sustainable development, and, ultimately, to achieve regional sustainable development strategies. Meanwhile, spatial differences deserve more attention, and the characteristics of different regions should be taken into account in order to achieve optimal spatial structure.

#### **Acknowledgments**

We are grateful for financial supports of Natural Science Foundation of China (NO. 41261112; NO. 41301652; NO. 41471116), the Science and Technology Support Program of Gansu Province (NO. 1304FKCA067) and the Program of Promoting the Scientific Research Capability of Young Teachers in Northwest Normal University (NO. NWNU-LKQN-10-20).

#### **Author Contributions**

Chenyu Lu and Yongjin Li designed the study and wrote the paper. Chengpeng Lu and Hengji Li analyzed the data. Chunjuan Wang and Weili Zhu contributed to data collection and processing. All authors have read and approved the final manuscript.

#### **Conflicts of Interest**

The authors declare no conflict of interest.

#### **References**

