Integrating Extended Fourier Amplitude Sensitivity Test and Set Pair Analysis for Sustainable Development Evaluation from the View of Uncertainty Analysis

It is of importance but great difficulty to objectively and quantitatively evaluate the sustainable development level, especially in the weight determination process and uncertainty evaluation. The traditional weight determination methods hardly reflect the coupling effect (interaction) among the indices. More importantly, conventional evaluation methods seldom consider the uncertainties of the indices in the index system. Thus, it is indispensable to apply a more comprehensive approach to solve these defects. This paper presents a new method to evaluate the sustainable development level. The approach integrates the advantages of the Extended Fourier Amplitude Sensitivity Test (EFAST) and Set Pair Analysis (SPA) (called EFAST-SPA). The EFAST algorithm is used to determine the indices’ weight, and the SPA is employed to handle the uncertain relations in the evaluation system and to calculate the sustainable development level. A quantitative evaluation on the agricultural sustainable development in the middle reaches of Heihe river has been conducted using the EFAST-SPA method. The results have been compared with the traditional entropy method and it was concluded that EFAST-SPA and entropy are highly in line with the actual development status. In most cases, the EFAST-SPA method can describe the development levels more accurately, which reflects a higher reliability and application value of this proposed approach. Moreover, the presented method deepens the understanding of sustainable development evaluation from the view of uncertainty analysis inside the evaluation system.


Introduction
Sustainability is an integral science which focuses on the dynamic relationship between humans and the environment [1].It especially emphasizes the vulnerability of the coupled human-society-environment system [2].With the rapid expansion of human society, the relationship between human society and the human living environment is in a unstable change [3].This inconstant change exacerbates the uncertainty in the evaluation system and causes several problems for the sustainable development level evaluation.One of change is how to analyze the uncertainties among the indices to identify the real sustainable development status using a reasonable and quantitative method.The evaluation model and indices weighing method selection are two main aspects of the sustainable development evaluation process.Both may bring a great influence on the accurate identification of the dynamic sustainable development status.
Numerous evaluation models are widely used in sustainable development research [4][5][6].While the quantitative evaluation [7] models are becoming mainstream in recent times.The quantitative models can be separated into four types [8].Including the Social-economic model, Ecology model, Systematic model, and some other emerging models [9].The Social-economic model [10] includes the Human Development Index (HDI) [11] and Input-output model [12].These kinds of methods integrate the social-economic indices into a composite index or employ the relationship between consumption and output to evaluate the sustainable development levels.The Social-economic model is easy to handle, but not all the data it needs can be easily obtained in every research area.Thus, it is not considered a universal model.There are some widely used methods in the Ecology model, such as the Ecological Footprint [13] and Emergy Analysis [14].This kind of model often converts the natural resources which are necessary for human survival and ecological restoration into a given unit, and then calculates the natural resource supplies with the same unit.Lastly, sustainable development levels can be determined by quantitatively analyzing the relationship between consumption and supplies.However, this kind of method pays too much attention on the aspect of ecology and less on the social-economic aspect, which may lead to a partial evaluation result.The Systematic model [15] is the most frequently used model in sustainable development evaluation research.This is because it is very comprehensive and easy to conduct.The Systematic model incorporates the social, economic, ecological, and environmental indices into the evaluation system.The sustainable development levels could be calculated via a simple linear weighing method.But a reasonable weight of each index is a high crucial factor for the Systematic model to work out an authentic evaluation result.Nowadays, emerging models which are derived from the interactive of interdisciplinary academic are getting more attention.Including the Neural Network model [16], Genetic Algorithm [17], Projection Pursuit model [18], and Support Vector Machine [19] (SVM) method.These innovative approaches are developing rapidly from the highly interdisciplinary of other research fields.Even there is a high value can be explored in those new methods, they are still not widely used due to the shallowness in current research.
The weighing method is essential for an accurate evaluation of the sustainable development level.A lot of evaluation researches have been carried out using different approaches in the past [20][21][22][23].Among them, several classic theories have been widely applied.Including Analytic Hierarchy Process (AHP), Entropy [24].AHP was proposed by Satty [25] when he researched on how to allocate the electric quantity that depended on the contribution of different industries.Through calculating the decision matrix and the indices' score estimated by experts, the weight of each index can be obtained, but this kind of method is overly relying on the experience of decision-makers.The Entropy method is the most common objective technique used by researchers [26,27].This method usually takes the normalized index entropy as the indices' weight.Even this method can work out precise weights in most cases.But the no-supervision weighing process requires pinpoint accuracy of data.Generally speaking, the traditional weighing methods seldom take the coupling effect among indices into account, which means that the common techniques rarely consider the uncertainties and the dynamic status stated in the beginning.In recent times, some uncertainty analysis approaches have been employed to evaluate the sustainable development level.Fuzzy Logic [28] is a representative solution to handle the uncertainties in evaluation problems which can be conducted without indices' weight.The membership function it contains can be applied to determine the sustainable development level directly.But the calculation process is a bit complicated and it is only applicable when the standard index system exists.It is hard to define a standard index system for dynamic sustainable development.So, the application of this kind method on sustainable development evaluation is not very wide.
The overall objectives of this research are to quantitatively estimate the sustainable development level and to work out a more reasonable classification of evaluation levels through the uncertainty analysis technique.However, lacking comprehensive uncertainty analysis may exist as a common weakness in the previous models and techniques.To overcome the deficiency, this paper proposes a more comprehensive approach integrating the Extended Fourier Amplitude Sensitivity Test (EFAST) [29] and Set Pair Analysis (SPA) [30], namely the EFAST-SPA method, to conduct sustainable development evaluation.The feasibility of applying EFAST to evaluate sustainable development levels has been validated in previous research [31].But only the sustainable development trends which presented evaluation results were analyzed in that research.To provide a more detailed explanation of the evaluation, this paper will convert the evaluation results from a trend description to a quantitative level classification.The classification can make the evaluation results in different areas more comparable.To be brief, two main specific issues are explored in this paper: (1) How to conduct the uncertainty analysis and obtain a quantitative sustainable development level from the EFAST-SPA method?(2) How to verify the accuracy of the evaluation result calculated by this proposed approach?To solve these questions, we completed the following work: Firstly, the EFAST is applied to fix indices' weight because of its advantage in quantitative global sensitivity analysis [31][32][33][34].Secondly, the SPA algorithm, which can analyze the unity, opposition, and uncertainty among different objects [35][36][37], was employed to analyze the uncertain relations in the evaluation system and to serve as evaluation model.Lastly, the evaluation results were compared with the results of the conventional entropy method and the positive analysis was also conducted synchronously.
The remainder of this paper comprises five sections.The materials and evaluation method are described in Section 2. The results and comparison analysis are described in Section 3. The discussion is provided in Section 4, and conclusions are presented in Section 5.

Study Area
Heihe River Basin is the second largest inland river basin in China with an area of 142,300 square kilometers [38].The study area (Figure 1) of this paper is located at the middle reaches of the Heihe River Basin.It is a main agricultural region.There are various landscapes inside the study area [39], of which the oases proportion is 28.9%.However, the Gobi Desert and the mountain region occupies 61.4% and 9.5%, respectively.The lack of water and farmland constrains the agriculture development naturally [40].Additionally, the increasing population and the expansion of human activities in recent years have brought many serious problems [41].Including soil desertification [42], soil salinization, and soil erosion [43].All of those problems have had a negative effect on the ecological environment and on the agriculture sustainable development in this study basin [44,45].

Brief Introduction of the Whole Evaluation Process
Figure 2 illustrates that how uncertainty analysis approaches were applied to evaluate the agricultural sustainable development levels.The evaluation system we constructed is based on the actual condition of the study area.The EFAST and the entropy were used to fix the weights of each index in the evaluation system.Weight comparison between the EFAST and entropy was also conducted.Based on the multi-index comprehensive analysis model and the weights determined by entropy, the development index can be calculated.Meanwhile, with the weights assigned by EFAST and the connection degree which was calculated by SPA, the development index of EFAST-SPA can be obtained.After that, the evaluation results of the subsystem and the entire system were acquired.To avoid repetition, only the subsystem results calculated by entropy and EFAST-SPA were selected to analyze.Positive analysis in this process was conducted to demonstrate the accuracy and the advantages of this new approach when it was compared with the entropy.Afterwards, the evaluation results of the entire system were classified into five levels.Depending on the analysis of evaluation results, the major development problems were revealed in the conclusion section.

Construction of Evaluation System
The evaluation process included many steps.The construction of the indicator system was the first work to be prepared for that process.Numerous international sustainable development evaluation systems were established during the last 30 years and many of them are widely used [46].The latest Sustainable Development Goals [47] system was released by the UN Sustainable

Construction of Evaluation System
The evaluation process included many steps.The construction of the indicator system was the first work to be prepared for that process.Numerous international sustainable development evaluation systems were established during the last 30 years and many of them are widely used [46].The latest Sustainable Development Goals [47] system was released by the UN Sustainable

Construction of Evaluation System
The evaluation process included many steps.The construction of the indicator system was the first work to be prepared for that process.Numerous international sustainable development evaluation systems were established during the last 30 years and many of them are widely used [46].The latest Sustainable Development Goals [47] system was released by the UN Sustainable Development Summit in 2015.However, the international systems are too complex to apply when the research is not on a global scale.While evaluating the sustainable development level of the local region, the local conditions and the problems can vary from one area to another.
The agricultural sustainable development problems in the middle reaches of the Heihe River Basin consisted of four main aspects.Including a long-term fragile ecological environment [48], water shortages [36], undeveloped economy and social development [49], and lower efficiency in agricultural production [50].Depending on these aspects a specific evaluation system (Table 1), which mainly focused on the local problems of the research area, was constructed.The system contained positive and negative indices.For the positive indices, the higher index value means a better development status.For the negative indices, a lower index value means better development status.
Eight negative indices have been set in our evaluation system.Therein, natural population growth rates were designed as negative indices because the population in the research area was large enough.The agricultural population proportion index was negative because a larger agricultural population led to labor redundancy with the progress in agriculture.Obviously, a chemical fertilizer application rate assuming 100% utilization pesticide use intensity, plastic sheeting use, and agricultural disaster rate, it was easy to understand their negative attributes.Besides, the smaller the water table value, the easier access to groundwater.Thus, the water table was set as at a negative index.The more groundwater irrigation in the Heihe middle reaches, the worse impact will appear.So, the groundwater irrigation was regarded as a negative index.
Among the positive indices, the only-child family rate can reflect the proportion of the families which raised only one child.In this research area, a one-child family usually means that this family can get access to higher quality living conditions.Thus, it is logical to define the only-child family rate as a positive index.Combining with the explanation of the negative indices, it is easy to understand the attribute of the rest positive indices.

Index Explanations
The indices selection was on the basis of a great number of related studies [51][52][53][54] and the researcher's experiences.All the index data can be obtained from the statistical yearbook or some specific websites.Each index in the evaluation was directly or closely related to the agricultural development in the study area.Therein, the agricultural economic development level subsystem was designed to focus on economy.Among them, the primary industry output per capita and crop input-output ratio were used to identify the agricultural economic level directly.The second and third industry output per capita were also adopted because these two indices would impede or facilitate the agricultural development to some extent.The livestock input-output ratio index was employed because the stockbreeding was always in conflict with the crop farming.The agricultural productive factors development level subsystem was mainly composed of the agricultural output per capita and the production efficiency.Meat and grain share per capita reflected the basic living substance level for human survival.Net income per capita of farmers was used to measure the farming economic efficiency.The mechanical effective utilization factor and effective irrigation coefficient were chosen as production efficiency indices to identify agricultural mechanization level.Agricultural social development level was designed on the basis of population, basic livelihood protection, and sociocultural development.Natural population growth rate, agricultural population proportion, only-child family rate, and the proportion of agriculture technicians among rural laborers were used to evaluate the labor resources development trend.The hospital beds share per 10 4 population and tap water supply rate were employed to calculate the basic livelihood protection level.Cultural and arts workers per 10 4 population represented the local cultural participation level.Agricultural resources and environment development level subsystem was aimed at agricultural resources consuming and environmental concerns.Therein, arable land per capita and the water table can reflect the resources whose status were germane to agricultural development.The chemical fertilizer application rate assuming 100% utilization, closed forest area, pesticide use intensity, plastic sheeting use, surface water and groundwater irrigation were used to evaluate the environmental impact caused by human activities.The agricultural disaster rate was employed to estimate the environmental influence from nature.

Data Collection and Pre-Processing
The data adopted in this paper was mainly derived from the Zhangye Yearbook (2006-2015) which comprised all the social and economic data.The unbalanced sustainable development states of different counties would make their evaluation comparison more difficult.To overcome this difficulty, the ratio data, which can be obtained after preprocessing, was adopted in this research.For instance, the index namely "Arable land per capita" was the ratio of arable land and total population in the certain area.Besides, some remote sensing data were also used.The surface and groundwater irrigation data were calculated from the "Monthly irrigation datasets (for both surface water and groundwater, 1981-2013) with 30 s spatial resolution over the Heihe River Basin" [55,56].Additionally, from the "Monthly groundwater table depth, soil moisture, and evapotranspiration datasets (1981-2013) with high spatial resolution over the Heihe River Basin" [55,56], the water table data can be extracted.All those remote sensing data were downloaded from the Cold and Arid Regions Science Data Center at Lanzhou.The data available is under strict checking process, which can ensure the data's accuracy and reliability.In this paper we extracted those three indices by administrative boundary.Then the averaged pixel value can be obtained with the tool named "Band Collection Statistics" in Arcgis 10.2 [57].After the extraction we found that even the smallest county Linze still contained 52 pixels, and that the number of pixels can satisfy the accuracy calculation needs.
This research input the index data into the SPA model to evaluate the sustainable development level.The raw data were various in units and cannot be compared directly, so a standardization of those data was the first process before putting them into the model.This paper employed the extreme standardization method [58] as the pre-process approach.

Model Selection
Model selection is an important step to conduct the sustainable development evaluation because it is closely associated with the weighing method.The Systematic model requires a very high accurate weight of index.Thus, this model can be an optimal choice to test the feasibility of EFAST weighing method.
The SPA algorithm can separate the indices into different levels perfectly by quantitatively analyzing the interaction among the indices.Besides, the modified SPA algorithm, which contains the indices weight, can offer a composite index to describe the sustainable development level as the Systematic models do.Thus, the SPA was selected as our evaluation model.

Entropy Evaluation Method
Entropy is the most frequently used method in previous research.Depending on the degree of discreteness of each index, this subjective method can calculate the entropy of the indices.Afterwards, the weight of each index can be obtained through the normalization process of their entropy.
In this research, the evaluation result, which was calculated by the entropy method, can be described using the following formula.
where, T j is the evaluation result of a certain county in its jth year; w i represents the weight of the ith index.

EFAST-SPA Evaluation Method
Fixing a reasonable weight of each index is the key process before evaluating the sustainable development level in the chosen evaluation model.It is hard to quantitatively analyze the complex uncertainties among the indices.That issue raises difficulties in the weight determination process.Besides, accurately classifying the development levels into different categories is another difficulty for the evaluation because the sustainable development status is dynamic and there exists no standard classifications to identify the evaluation results.To solve the weighing problem, the EFAST method was adopted to fix the indices' weight.Based on the analysis of EFAST, both the first-order sensitivity (the main effect, also called MSI) indices, which represent a single index influence on the model output, and the total sensitivity indices (TSI, the total effect), which reflects the main effect and the coupling effect among indices on the model output, can be calculated.The TSI normalization result would offer the weight of each index.The feasibility and advantage of this weighing method has been verified in previous research [31].To provide a reasonable evaluation classification, this research has improved the EFAST with the SPA algorithm to develop a more improved and quantitative evaluation method.The detail of weighing process can be found from a previous work [31] and the combined method can be described subsequently.
All the weighing processes were compiled into Python code which can execute the program efficiently.The weights obtained from this program were adopted into the SPA model to calculate the sustainable development level.If there were two groups whose elements were related from one another, the SPA method would set the elements of those two groups into one-to-one correspondence pairs.The certain-uncertain relations can be analyzed by Equation ( 2): In this expression, the µ represents the connection degree, a, b and c are the degree of identity, uncertainty and opposites, respectively.Additionally, the sum of a, b and c is 1, and the range of them can be defined as 0 ≤ a, b, c ≤ 1. Commonly, the value of opposite coefficient j = −1.The range of the uncertainty coefficient i is (−1, 1) [59].
There always exists various degrees of uncertainty relations between the comparison groups.The expression of the connection degree can be extended into a multi one [60].
In this paper, a standard table, namely the classification of index level (Table 2), was constructed.The boundary of each level was set by equal-interval.For the positive index, the boundary of the first level is the highest in the system, and the boundary of the fourth level is the lowest one.The negative indices' boundary was then set by the opposite conduct.By comparing the standard level data with the original index data, the connection degree can be calculated with Equation (3).The connection degree of positive and negative indices should be calculated by different expressions.The negative index was processed by Equation (4): where Sustainability 2018, 10, 2435 9 of 23 The following expression was used to calculate the connection degree of positive indices. where To simplify the calculation process, the uncertainty coefficient was divided by their range with equal-interval and the value of coefficients were defined as i 1 = 0.5, i 2 = 0, i 3 = −0.5.
After the definitions above, the connection degree of the standard group and the original group can be calculated.Combined with the weights determined by EFAST method, the sustainable development index µ(A, B) can be expressed as the following equation: where, K represents the number of standardization level, µ(A n , B K ) means the connection degree of H(A n , B K ), and ω n is the weight of index n.
The last process was to quantitatively classify the sustainable development index into different levels.The method of classification was introduced as the Table 3: Notes: M and Std are the mean and standard deviation of sustainable development index, respectively.L i represents the ith level.

Comparison of MSI and TSI
According to the weighing process of EFAST, both the MSI and TSI of each index can be calculated.There are a total of six counties in the research area and the MSI value was slight less than the TSI value in all these counties.To be brief, only Ganzhou was selected as a sample to conduct the comparison between MSI and TSI.
From Figure 3, it can easily be seen that the difference between the MSI and TSI in Ganzhou is very small.For instance, the biggest difference between them was found in 2005 and it was only 0.0011, while the smallest difference was only 0.0006.For the entire period, the biggest difference was 0.0012 and the smallest was 0.0004.The difference reflected the coupling effect among the indices to the model output.But those differences were not as significant as expected.These results may be caused by the simple lineal model; because the EFAST algorithm can be more efficient when it is applied in the analysis of nonlinear and high dimension models.Even though, the weak coupling effect still cannot be neglected when an accurate evaluation is needed.Besides, the differences of SIs among each index from 2005-2007 were more significant than the 2008-2013 period.Only a few indices' SIs can be amounted to 80% in the entire SIs space from 2005-2007 because during that period only a small number of indices were paid attention by the farmers or governments.During the 2008-2014 period, more indices were attached importance.That may be because of more comprehensive thoughts of decision-makers.However, the differences of SIs in 2014 were significant again, which may be caused by the policy orientation of the government.
only a small number of indices were paid attention by the farmers or governments.During the 2008-2014 period, more indices were attached importance.That may be because of more comprehensive thoughts of decision-makers.However, the differences of SIs in 2014 were significant again, which may be caused by the policy orientation of the government.Legend: the sequence number is in accordance with Table 1.
Figure 3. SIs of each index in Ganzhou.Legend: the sequence number is in accordance with Table 1.

Weights Determined by Entropy and EFSAT
The agricultural sustainable development levels were variable because different areas would be confronted with different problems.So, in this paper, different weights were set for the same index in different areas.The weights assigned by the entropy and EFAST are shown in Figure 4.

Weights Determined by Entropy and EFSAT
The agricultural sustainable development levels were variable because different areas would be confronted with different problems.So, in this paper, different weights were set for the same index in different areas.The weights assigned by the entropy and EFAST are shown in Figure 4.   From Figure 4, it can be clearly seen that high conformity exists between the weights assigned by entropy and EFSAT methods.Taking Ganzhou as an example, closed forest area, the proportion of agriculture technician among rural laborers, meat share per capita, livestock input-output ratio, and pesticide use intensity were regarded as the top five indices which had a higher weight than other indices in entropy method.While according to the EFAST, the top five indices included meat share per capita, livestock input-output ratio, the proportion of agriculture technician among rural laborers, closed forest area, and cultural and arts workers per 10 population.Obviously, we can find that four indices were the same in the top five of both methods.When the top five indices of the six counties were concluded by entropy and EFAST, there were also 18 of the same indices amongst the 30 indices.Entropy is a widely used approach in weight determination research, the 60% overlap rate of higher weight indices shows the feasibility of EFAST.Further validation would need to be conducted via the evaluation results comparison of these two methods.
Weight determination is a process to screen out the important indices to the target from the index system and fix the important indices with high weights.The more discrete the weights, the better the evaluation results will be.The dispersed degree of entropy and EFAST is exhibited in Table 4.In this table it was found that all the dispersed degrees in EFAST were greater than that in entropy.The smallest difference existed in Sunan, however, the ratio between EFAST and entropy was still high as 1.41.From Figure 4, it can be clearly seen that high conformity exists between the weights assigned by entropy and EFSAT methods.Taking Ganzhou as an example, closed forest area, the proportion of agriculture technician among rural laborers, meat share per capita, livestock input-output ratio, and pesticide use intensity were regarded as the top five indices which had a higher weight than other indices in entropy method.While according to the EFAST, the top five indices included meat share per capita, livestock input-output ratio, the proportion of agriculture technician among rural laborers, closed forest area, and cultural and arts workers per 10 4 population.Obviously, we can find that four indices were the same in the top five of both methods.When the top five indices of the six counties were concluded by entropy and EFAST, there were also 18 of the same indices amongst the 30 indices.Entropy is a widely used approach in weight determination research, the 60% overlap rate of higher weight indices shows the feasibility of EFAST.Further validation would need to be conducted via the evaluation results comparison of these two methods.
Weight determination is a process to screen out the important indices to the target from the index system and fix the important indices with high weights.The more discrete the weights, the better the evaluation results will be.The dispersed degree of entropy and EFAST is exhibited in Table 4.In this table it was found that all the dispersed degrees in EFAST were greater than that in entropy.The smallest difference existed in Sunan, however, the ratio between EFAST and entropy was still high as 1.41.

Analysis of Subsystem Evaluation Results
The evaluation results of the subsystem reflect the balance degree of the four subsystems during their development phase.The Subsystem Development Index (SDI) was calculated with no classification of levels.The reason was that we only want to analyze the development trend on each subsystem, and the SDI with no level classification can be more intuitional for analysis via a straight-line graph.To avoid repetition, only Ganzhou and Minle were selected to analyze the difference of SDI calculated by entropy and EFAST-SPA.
From Figures 5 and 6, it can be seen that there was a high consistence between the SDI marked with EFAST-SPA and entropy.But it was noted that subsystems in Ganzhou experienced an unbalanced development from 2005 to 2014.The agricultural economic development level showed a slow development speed during that period.Meanwhile, the agricultural social subsystem developed rapidly at the first half research period, but its trend declined in the following years.Especially in 2013, the agricultural social SDI nearly decreased to the level of 2005.The agricultural productive factors, resources, and environment subsystem experienced a circuitous rise during the entire research period.Thus it can be concluded that the development pace among those four subsystems in Ganzhou was not in line and the pace did not satisfy the balance demand of sustainable development.
Similarly, there existed an unbalanced development trend among these subsystems in Linze.The agricultural productive factors and the agricultural social subsystem showed a rapid development in the first several years and then stayed steady or subtlety declined in the following years.Besides, the agricultural economic development level expressed a very slow development speed during the entire time.However, agricultural resources and environment development level even experienced a sharp downward trend from 2009 to 2012, this worsening ecological environment phenomenon violated the concept of sustainable development.
Overall, both the entropy and the EFAST-SPA have evaluated the subsystem sustainable development level objectively.Entropy is a widely accepted evaluation method.Thus, this paper selected its results as our reference standard to verify the EFAST-SPA method.Thus, similar evaluation results between entropy and EFAST-SPA further verified the accuracy of our proposed method.

Analysis of Subsystem Evaluation Results
The evaluation results of the subsystem reflect the balance degree of the four subsystems during their development phase.The Subsystem Development Index (SDI) was calculated with no classification of levels.The reason was that we only want to analyze the development trend on each subsystem, and the SDI with no level classification can be more intuitional for analysis via a straightline graph.To avoid repetition, only Ganzhou and Minle were selected to analyze the difference of SDI calculated by entropy and EFAST-SPA.
From Figure 5 and 6, it can be seen that there was a high consistence between the SDI marked with EFAST-SPA and entropy.But it was noted that subsystems in Ganzhou experienced an unbalanced development from 2005 to 2014.The agricultural economic development level showed a slow development speed during that period.Meanwhile, the agricultural social subsystem developed rapidly at the first half research period, but its trend declined in the following years.Especially in 2013, the agricultural social SDI nearly decreased to the level of 2005.The agricultural productive factors, resources, and environment subsystem experienced a circuitous rise during the entire research period.Thus it can be concluded that the development pace among those four subsystems in Ganzhou was not in line and the pace did not satisfy the balance demand of sustainable development.
Similarly, there existed an unbalanced development trend among these subsystems in Linze.The agricultural productive factors and the agricultural social subsystem showed a rapid development in the first several years and then stayed steady or subtlety declined in the following years.Besides, the agricultural economic development level expressed a very slow development speed during the entire time.However, agricultural resources and environment development level even experienced a sharp downward trend from 2009 to 2012, this worsening ecological environment phenomenon violated the concept of sustainable development.
Overall, both the entropy and the EFAST-SPA have evaluated the subsystem sustainable development level objectively.Entropy is a widely accepted evaluation method.Thus, this paper selected its results as our reference standard to verify the EFAST-SPA method.Thus, similar evaluation results between entropy and EFAST-SPA further verified the accuracy of our proposed method.

Comprehensive Evaluation Results and Positive Analysis
This section would exhibit the comprehensive development index (CDI) calculated by entropy and EFSAT-SPA methods (see Figure 7).According to the sustainable development level classification method introduced in Section 2.4.3, the evaluation results of each area were classified into different levels (see Figure 8).

Positive Analysis of Entropy and EFAST-SPA
Figure 7 illustrated the CDI of those six counties.The trend it showed also reflected the consistence pattern between the entropy and the EFAST-SPA approach on the whole study area.Only

Comprehensive Evaluation Results and Positive Analysis
This section would exhibit the comprehensive development index (CDI) calculated by entropy and EFSAT-SPA methods (see Figure 7).According to the sustainable development level classification method introduced in Section 2.4.3, the evaluation results of each area were classified into different levels (see Figure 8).

Comprehensive Evaluation Results and Positive Analysis
This section would exhibit the comprehensive development index (CDI) calculated by entropy and EFSAT-SPA methods (see Figure 7).According to the sustainable development level classification method introduced in Section 2.4.3, the evaluation results of each area were classified into different levels (see Figure 8).

Positive Analysis of Entropy and EFAST-SPA
Figure 7 illustrated the CDI of those six counties.The trend it showed also reflected the consistence pattern between the entropy and the EFAST-SPA approach on the whole study area.Only The trend it showed also reflected the consistence pattern between the entropy and the EFAST-SPA approach on the whole study area.Only a small part of these dots showed obvious differences.The inconsistent year should be emphasized in the positive analysis.For instance, the CDI in 2013 calculated by entropy in Ganzhou showed a significant downward trend, but the EFAST-SPA showed a rising trend.Compared with 2012 (Table 5), many positive indices decreased significantly, but few negative indices improved.So, the evaluation results of entropy were more reasonable.However, in 2006, the negative indices in Ganzhou, such as pesticide use intensity and plastic sheeting use, experienced a weak rise, and that weak rise lead to a significant improvement in the positive indices.This improvement can probably bring a rise in CDI.Thus, the rising trend of the EFAST-SPA was more accurate than the entropy for 2006.Notes: the number of index is consistent with Table 2.
During 2006-2008, CDI in Minle also showed a difference between entropy and EFAST-SPA.From Table 6, it can be noted that most positive indices in 2006 improved compared to 2005.For the negative indices, only the Natural population growth rate increased.But the benefits that came from the positive indices compensated for the loss that the negative indices brought.Therefore, the rising trend in entropy was correct through the positive analysis.However, the significant rising trend of 2007 and 2008 in the EFSAT-SPA method can be more accurate than the entropy, because, in comparison with 2007, most positive indices in 2008 experienced an improvement, and meanwhile, most negative indices decreased to a variable extent.Notes: the number of index is consistent with Table 2.
There also existed some differences between entropy and EFAST-SPA in Shandan during 2012 and 2013.As shown in Table 7, 15 positive indices remained the same or decreased in 2013, and 6 negative indices remained unchanged or increased, which could probably cause a sharp downward trend in the CDI.This downward trend was not evidently expressed in the entropy method.Notes: the number of index is consistent with Table 2.
Only six obvious inconsistent years existed among the sixty points in Figure 7. Two of the inconsistent years showed that the entropy was more accurate.The remaining four stated that the EFAST-SPA performed better.The entropy often provided a narrow range CDI.So, it cannot perform as well as EFAST-SPA when radical change occurs.That was the main reason why EFAST-SPA performed better in the other four inconsistent years.In this paper, the 90% consistent results can verify the feasibility of our proposed method.Overall, the EFAST-SPA method may cause some deviations in this research.However, this small deviation cannot affect its application.

Sustainable Development Level Classification of Entropy and EFAST-SPA
The evaluation results were classified into five levels and shown in Figure 8.The difference in levels between entropy and EFAST-SPA was not significant.In this section, only the classifications of EFSTA-SPA were selected to rank the sustainable development level of the six counties from 2005-2014.The ranking was calculated by summing up the levels of each county.The smaller the sum, the better the sustainable status.Thus, the ranking was listed as: Sunan (20) > Shandan (36) > Gaotai (40) > Linze (41) > Minle (44) > Ganzhou (45).Only six obvious inconsistent years existed among the sixty points in Figure 7. Two of the inconsistent years showed that the entropy was more accurate.The remaining four stated that the EFAST-SPA performed better.The entropy often provided a narrow range CDI.So, it cannot perform as well as EFAST-SPA when radical change occurs.That was the main reason why EFAST-SPA performed better in the other four inconsistent years.In this paper, the 90% consistent results can verify the feasibility of our proposed method.Overall, the EFAST-SPA method may cause some deviations in this research.However, this small deviation cannot affect its application.

Sustainable Development Level Classification of Entropy and EFAST-SPA
The evaluation results were classified into five levels and shown in Figure 8.The difference in levels between entropy and EFAST-SPA was not significant.In this section, only the classifications of EFSTA-SPA were selected to rank the sustainable development level of the six counties from 2005-2014.The ranking was calculated by summing up the levels of each county.The smaller the sum, the better the sustainable status.Thus, the ranking was listed as: Sunan (20) > Shandan (36) > Gaotai (40) > Linze (41) > Minle (44) > Ganzhou (45).

Sustainable Development Level Classification on Spatial Map
To demonstrate the level classification more directly and analyze the sustainable development relationship spatially, a spatial map containing the classification could be a useful tool for the analysis.To keep the paper reasonably concise, only the classifications in 2005, 2010, and 2014 (see  were exhibited as representatives.
The color in the maps showed that all the six counties were becoming greener from 2005 to 2014.Which reflected a better development trend during the research period.Besides, the sustainable development status was better when the area was further towards Ganzhou.This phenomenon existed in all the maps.

Comparison of the Weights Assigned by Entropy and EFAST
A reasonable weight for the index is a primary and key process for evaluation.Few previous studies were based on the analysis of the coupling effect among the indices.Applying the global sensitivity analysis method aims at offering a deeper understanding of weight determination and providing more reasonable weights for the evaluation indices.This paper found that the EFAST algorithm, which considered the coupling effect among the indices, was feasible.Furthermore, the EFAST algorithm can be able to fix a more comprehensive and reasonable weight than the traditional entropy method.That is because a more discrete weight range can efficiently separate the indices into different important levels.Additionally, this subjective method has been widely used in the other models and its prudent calculation has been verified.

Analysis of the Unbalance Development in Subsystems
The purpose of setting subsystems is to analyze the balance degree of a sustainable development level from a micro-scale.This paper analyzed the four subsystems in each county and found various problems in different counties.Therein, the unbalanced development pace of the subsystems was a common defect among all of these six counties.Due to the results calculated by entropy, EFAST-SPA showed the similar trend.Only the EFAST-SPA results would be analyzed in this section.The following analyses would reveal the major problems in each county and these problems can guide the decision-maker to make more reasonable decisions [61].
The most serious problem in Ganzhou was the agricultural social subsystem, because this subsystem developed weakly, while the other three subsystems developed well during the research period.Especially the agricultural sources and environment section, the SDI of this subsystem in 2014 was approximately 2.4 times what it was in 2010.The agricultural social subsystem also developed a bit more slowly than other subsystems in Minle.However, when it was compared with other counties, the other three subsystems improved at a slow speed, too.The unbalance trend occurred in Gaotai in the first half of the research period, but the development pace was becoming more balanced between 2010 and 2014.The agricultural and social section was the only subsystem which improved significantly during 2005-2011.However, its SDI showed a bit of a decrease during 2012-2014.Moreover, the SDI of the rest of the subsystems nearly remained the same over the entire ten year period.The agricultural resources and environment in Sunan decreased distinctly when compared with the rest subsystems.It obvious that high-speed development was obtained at the cost of excess resource exploitation.
The agricultural productive factors subsystem was the only low-speed development section in Linze.Thus, the mechanization of agriculture could be an efficient solution to improve the sustainable development level in this area.Overall, the information deduced from the evaluation can offer a clear understanding of the development problems the counties faced.This deeper understanding could help the decision-makers to regulate a more targeted ordinance to solve the specific problems.

Analysis of EFAST-SPA Sustainable Development Levels
Positive analysis demonstrated that the EFAST-SPA results were in line with the authentic development status in most cases.Besides, the span of CDI calculated by EFAST-SPA was wider than that calculated by entropy, which indicated that the EFAST-SPA can evaluate the agricultural sustainable development level more fully.This advantage could offer more detained information which the decision-maker needs.The following analysis was based on the EFAST-SPA evaluation results.
The sustainable development levels in Sunan retained a good development status during the entire research period.However, the status of Shandan was changing between weak sustainable and basic sustainable.Besides, Linze and Gaotai kept the same development pace and both of them have been improved to basic sustainable from unsustainable.Additionally, these two counties experienced the greatest improvement among all of the six counties.The counties Ganzhou and Minle still remained in the weak sustainable status at 2014; they were the areas who needed a better improvement urgently.
From the perspective of spatial analysis, we can find that the status could be better if the county was located further from Ganzhou.Two reasons could be used to explain that result.Firstly, Ganzhou is the central area of economy, politics, and culture.The pursuit of economy improvement may have been a bit excessive during the research period, and that probably lead a bad impact on the local ecology and environment.Sustainable development is a comprehensive concept which considers not only the economy but also includes the ecology and environment.Thus, the worsening condition of ecology and environment will cause a low sustainable development level.The development style of Ganzhou will also influence the nearby counties the most.Secondly, natural conditions are another factor which can affect agricultural development.The altitude in Sunan ranges from 1700-5565 m, which is the highest in all of the six counties.This high elevation has reduced human activity in this area.Less population will bring less impact on the environment and ecology.Besides, the small population also contributes to improve per capita resources and economic growth.In conclusion, reducing the intensity of human activities will benefit the agricultural sustainable development.

Challenges of EFST-SPA Evaluation Method for Future Research
Although the EFAST-SPA method can identify the sustainable development status successfully.The evaluation process can be further improved from the following perspectives.Firstly, the evaluation system we constructed was on the basis of our research experiences and some indices were removed because of data unavailability.So, an improvement in this system would make it more representative and drive the evaluation result more authentic.Secondly, for convenient calculation purposes, the uncertainty coefficients in the SPA model were simplified by dividing them at equal intervals.Even though this is a common conduction, we still assume that improved uncertainty coefficients can offer more accurate results.Thirdly, there is no standard level boundary of the indices to follow because the standard boundaries are various in different study areas.Thus, the standard boundary we proposed is only appropriate for our research.Defining a universal standard boundary of the evaluation indices is a challenging topic for the further research.

Conclusions
In this research, two uncertainty analysis methods (EFAST and SPA) were integrated to evaluate the agricultural sustainable development level in the middle reaches of the Heihe River Basin in consideration of the uncertain relationships among the indices and the dynamic development process.The final result of sustainable development level during 2005-1014 was Sunan > Shandan > Gaotai > Linze > Minle > Ganzhou.This result was in nearly complete accordance with that of the traditional entropy method.The similar result can verify the feasibility of the EFAST-SPA method.More importantly, the method we proposed can describe the development patterns more completely than the entropy method, which indicated that our method can grasp the microchanges of the dynamic development.Besides, the positive analysis which was used to verify the accuracy of our results also proved that our evaluation was in line with the local actual development level in most cases.Thus, our result can be regarded as an accurate evaluation under the index system in this research.Meanwhile the accuracy of this result also reversely verified the reasonability of the weights assigned by our proposed method.In conclusion, the EFAST-SPA technique can provide objective and reasonable weight for the indices via a more comprehensive technique, it also can offer a quantitative and accurate evaluation for the sustainable development evaluation in the view of uncertainty analysis.Thus, we believe that the EFAST-SPA method can be applied as a feasible and reliable method for sustainable development evaluation research.

Sustainability 2018 , 22 Figure 1 .Figure 2 .
Figure 1.Map of the middle reaches in Heihe river showing its six municipalities and corresponding population size, land cover, and major rivers.

Figure 1 .
Figure 1.Map of the middle reaches in Heihe river showing its six municipalities and corresponding population size, land cover, and major rivers.

Sustainability 2018 , 22 Figure 1 .
Figure 1.Map of the middle reaches in Heihe river showing its six municipalities and corresponding population size, land cover, and major rivers.

Figure 2 .
Figure 2. Strategy for agriculture sustainable development evaluation using the EFAST-SPA method.

Figure 2 .
Figure 2. Strategy for agriculture sustainable development evaluation using the EFAST-SPA method.

Figure 3 .
Figure 3. SIs of each index in Ganzhou.Legend: the sequence number is in accordance with Table1. 0

Figure 8 .
Figure 8. Sustainable development level classification of entropy and EFSAT-SPA.

Figure 8 .
Figure 8. Sustainable development level classification of entropy and EFSAT-SPA.

Sustainability 2018 , 22 Figure 9 .
Figure 9.The agricultural sustainable development level of the six counties in 2005.

Figure 10 .
Figure 10.The agricultural sustainable development level of the six counties in 2010.

Figure 9 . 22 Figure 9 .
Figure 9.The agricultural sustainable development level of the six counties in 2005.

Figure 10 .
Figure 10.The agricultural sustainable development level of the six counties in 2010.Figure 10.The agricultural sustainable development level of the six counties in 2010.

Figure 10 .
Figure 10.The agricultural sustainable development level of the six counties in 2010.Figure 10.The agricultural sustainable development level of the six counties in 2010.

Figure 10 .
Figure 10.The agricultural sustainable development level of the six counties in 2010.

Figure 11 .
Figure 11.The agricultural sustainable development level of the six counties in 2014.Figure 11.The agricultural sustainable development level of the six counties in 2014.

Figure 11 .
Figure 11.The agricultural sustainable development level of the six counties in 2014.Figure 11.The agricultural sustainable development level of the six counties in 2014.

Table 1 .
Evaluation system of agricultural sustainable development in the middle reaches of the Heihe River Basin.

Table 2 .
Standard level of the evaluation system.
Notes: The negative indexes are marked with "*".

Table 3 .
The classification of the agricultural sustainable development level in Heihe middle reaches.

Table 4 .
The dispersed degree of weights in the six counties.

Table 4 .
The dispersed degree of weights in the six counties.

Table 7 .
The 2012 and 2013 index value in Shandan.