Next Article in Journal
The Effect of Airport Self-Service Characteristics on Passengers’ Perceived Value, Satisfaction, and Behavioral Intention: Based on the SOR Model
Next Article in Special Issue
Stakeholder Roles and Perspectives on Sedimentation Management in Small-Scale Irrigation Schemes in Ethiopia
Previous Article in Journal
Feasibility of the Fintech Industry as an Innovation Platform for Sustainable Economic Growth in Korea
Previous Article in Special Issue
Designing Optimum Water-Saving Policy in China Using Quantity and Price Control Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating Agricultural Sustainability Based on the Water–Energy–Food Nexus in the Chenmengquan Irrigation District of China

1
College of Agricultural Engineering, Hohai University, Nanjing 210098, China
2
College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China
3
Shandong Hydrology and Water Resource Bureau of Yellow River Conservancy, Jinan 250000, China
4
Institute for Water Quality and Resource Management, TU Wien, Vienna 1040, Austria
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(19), 5350; https://0-doi-org.brum.beds.ac.uk/10.3390/su11195350
Submission received: 14 August 2019 / Revised: 15 September 2019 / Accepted: 20 September 2019 / Published: 27 September 2019

Abstract

:
Agriculture is one of the largest consumers of water and energy. This paper evaluated the agricultural sustainability of the Chenmengquan irrigation district of China based on the water–energy–food nexus. One objective weighting method and one subjective weighting method were integrated, based on game theory, and a matter–element model was constructed to evaluate agricultural sustainability for the research region. The sensitivity of each index to the evaluation class was also analyzed. The results showed that agricultural sustainability was moderate in 2006–2012 and high in 2012–2015. The indexes, which represent water-use efficiency and yield per unit area of crops, had higher sensitivities in the context of the present case study. The results also indicated that agricultural sustainability had a comparatively positive trend between 2012 and 2015, and that pesticide utilization was the most important issue for agricultural sustainability. The approach of using the combination of a weighting method, based upon game theory, and the use of the matter–element model provides a guide for the evaluation of agricultural sustainability.

1. Introduction

The issue of agricultural sustainability has received considerable critical attention in the world with the changing environment and the growing population [1]. In the case of limited arable land resources, in order to satisfy the expected demands for food, agriculture must be transformed, in order to feed a growing global population and provide the basis for economic growth and poverty reduction [2]. More than 80% of agricultural production growth by 2050 is expected to come from increased productivity on presently cultivated land [3]. However, intensive agricultural practices have resulted in serious environmental degradation, including groundwater level decline, the loss of biodiversity, and soil and water pollution [4]. Improving agricultural sustainability is the inevitable choice for agricultural development in the future. Operationalizing the assessment of agricultural sustainability is an important premise in order to promote it. The assessment of agricultural sustainability aims to provide an efficient tool for decision-makers to determine which actions should be preferentially taken in an attempt to promote agricultural sustainability [5].
Studies over the past two decades have developed a wide variety of methodologies for evaluating agricultural sustainability. Most of methodologies have a similar structure constructed by index selection and data gathering, reference values definition, weighting indexes, and index aggregation [6,7,8,9,10,11,12]. However, there are huge differences between selections of indexes. Some methods take three dimensions, accounting for economic, social, and environmental issues in the frame of a holistic vision of agricultural sustainability [6,7]. These methods select more quantitative indexes to define a methodology, but some indexes of economic and social dimensions (e.g., the income of agricultural producers) cannot sufficiently reflect the agricultural sustainability; those indexes depend more on the economic situation in the specific region and at the specific time. Some methods establish the index system based on general attributes (e.g., the productivity, stability, and compatibility) of agricultural sustainability. Those kinds of index systems have a good thinking structure regarding the agricultural sustainability as a whole [8]. Due to the highly complex nature of systemic indicators, those indexes remain qualitative rather than quantitative parameters [13]. These qualitative indexes must be artificially assigned numerical values for any quantitative evaluation of agricultural sustainability. In this paper, indexes were selected from the perspective of agricultural intensity based on the concept of the water–energy–food nexus (WEFN).
The Bonn conference of 2011 promoted the concept of the water–energy–food nexus (WEFN) in the international discussion of sustainable development [14,15]. Water, energy, and food are closely linked in agricultural practices. Agriculture has accounted for 70% of the world’s total freshwater extraction, making it the largest water user. Food production and supply chains consume about 30% of global energy [16]. The WEFN has emerged as a new perspective to agricultural sustainability, with its conceptual relevance, and pragmatic potential emphasized by many policymakers [17]. The aim of this research is not to analyze the impact of water resources and energy consumption on agricultural production, but to evaluate the status of agricultural sustainability from an integrated perspective. Based on the WEFN, the level of agricultural water use, the intensity of agricultural energy use and food productivity are regarded as evaluation objectives. Taking advantage of “degree” instead of “quantity” to select indexes provides a useful approach for sustainability assessment across different scales.
Many mathematical models were applied to quantitative evaluation of sustainability, such as the matter–element model [18], the artificial neural network model [19], and the gray correlation model [20]. Agricultural sustainability as the evaluation object can be regarded as a composite whole with different aspects. The attribute diversity of agricultural sustainability can be assumed as an incompatibility problem. Therefore, the matter–element model is an appropriate mathematical tool for assessing agricultural sustainability. Considering that different indexes have different impacts on agricultural sustainability, this paper also adopts a comprehensive weighting method which combines the G1 method and the entropy method to determine the weights of indexes. The G1 method is a subjective method. It is improved according to the analytic hierarchy process (AHP). That method overcomes the large amount of calculations and consistency checks [21]. The entropy method was proposed by Shannon and Weaver [22]. The entropy method is a method to determine the weight of the index according to the amount of information provided by the observational values of each index. The comprehensive weighting method can reach consistency or compromise between the subjective weights and the objective weights [23].
The purpose of this study was to develop a methodology to evaluate agricultural sustainability in a given region with intensive agricultural practice. The index system was established based on the WEFN. A comprehensive weighting method was used to weight indexes. The matter–element model was applied to reference values’ definitions and indexes’ aggregations. Through the analysis of the results calculated by the matter–element model, the key information affecting agricultural sustainability can be identified.

2. Case Study

The Chenmengquan irrigation district (117.17°–117.34° E, 36.70°–36.89° N) is located in Jinan City, the capital of Shandong province in China (Figure 1). The study area is important for agricultural production in Jinan City. The total population of the irrigation district is 223,000, of which the agricultural population is 185,000. The total area of the irrigation district is approximately 20,000 ha. The main crops in the irrigation district are wheat, maize, and vegetables. With the continuous adjustment of the proportion of different crop acreages, the vegetable planting area has increased substantially. The study area is located in the mid-latitude zone and belongs to the warm, temperate, semi-humid monsoon climate. Rainfall varies greatly in different seasons. The average annual precipitation is 673 mm (1956–2016) and 70% of the annual precipitation is concentrated between June and September [24].
Irrigation plays an important role for agriculture in the study area. Irrigation water is mainly supplied by the Yellow River and groundwater. With increased use of water for agriculture the exploitation of groundwater gradually exceeded the recharge, resulting in the continuous decline of the groundwater level in 2000–2006. Since 2006, the administrative department has promoted the protection of groundwater. However, there are still other problems like soil and water pollution caused by the excessive application of fertilizer and pesticide.

3. Methods

Based on the above introduction, the steps of the proposed methodology are shown in Figure 2.

3.1. Index System

To assess agricultural sustainability, representative indexes are identified in Table 1 through a comprehensive perspective based on WEFN (ci represents index number, i = 1, 2, …, 10).

3.1.1. Water

Irrigation plays an important role in increasing the output of agricultural products. However, it also faces the practical problems of improving the efficiency of agricultural water use and controlling the total amount of water used. Cao et al. proposed an evaluation index for agricultural water use effect based on water footprint (WF) and irrigation development [25]. The indexes c1 (agricultural blue water proportion), c2 (water-use efficiency), and c3 (irrigation proportion of arable land), which represent agricultural water use effects, were used as indexes in the water sector.
The index c1 (agricultural blue water proportion) is based on the concept of WF proposed by Hoekstra [26]. It indicates the dependence of agriculture on irrigation water. The index expression for c1 is given in Equation (1), as the ratio between agricultural blue water footprint (AWFb) and the agricultural water footprint (AWF).
c 1 = A W F b / A W F
AWFb and AWF are calculated as in Equations (2) and (3)
A W F b [ m 3 ] = c r = 1 c r n E T b , c r × 10 × S c r
A W F [ m 3 ]   = c r = 1 c r n ( E T b , c r + E T g , c r ) × 10 × S c r
ETb,cr and ETg,cr (in mm) are the blue and green evapotranspiration components of a given crop, which are computed by the CropWat® software. Scr (in ha) is the total area cultivated of a given crop.
The index c2 (water-use efficiency) is essential regarding the water consumption. It has a high correlation with the infrastructural performance of the irrigation network, the farm management techniques, and irrigation methods. The index c3 (irrigation proportion of arable land) is the ratio of effective irrigated area to arable land area. The effective irrigated area is the area of farmland which is relatively flat, has a certain water source and irrigation facilities, and can be irrigated normally in an average year. The data of c2 (water-use efficiency) and c3 (irrigation proportion of arable land) were derived from the Statistical Yearbook of Jinan City.

3.1.2. Energy

The energy consumed by agriculture can be divided into direct energy and indirect energy [27]. Direct energy includes fuel and electricity mainly used for irrigation and agricultural machinery’s operation. Indirect energy refers to the energy used to produce fertilizers and pesticides. In this research, energy consumption intensity is expressed by four indicators: c4 (energy utilization for irrigation amount per arable land), c5 (agricultural machinery power per arable land), c6 (fertilizer utilization amount per arable land), and c7 (pesticide utilization amount per arable land).
The index c4 (energy utilization for irrigation amount per arable land) is given in Equation (4)
c 4 = E i r / S
where S represents the area of the arable land and Eir represents energy utilization for irrigation, which refers to the energy requirement from-source-to-field. It not only depends on the specific water source, but also on the distribution of irrigation systems, the characteristics of the conveyance with their respective water losses and operational pressures [28]. For simplicity, the energy requirement for conveyance was not included, hence the water source was assumed to be on-farm [29]. For surface irrigation, the energy for surface water transport was mainly driven by gravity. Consequently, energy used for groundwater abstraction was computed as Eir, shown in Equation (5)
E i r = V o l u m e g × H t / ( 367 × η )
H t = H l i f t + H p r
where Volumeg represents the gross amount of water from a groundwater source. η represents the pump efficiency; and Ht, Hlift, and Hpr represent the total hydraulic head, the elevation head, and the nominal operating pressure of the irrigation system, respectively. The data of c5 (agricultural machinery power per arable land), c6 (fertilizer utilization amount per arable land), and c7 (pesticide utilization amount per arable land) were derived from the Statistical Yearbook of Jinan City.

3.1.3. Food

The indexes of the food sector were formulated according to the actual agriculture production in the given area. Wheat, maize, and vegetables make up more than 90% of the arable land at the study site. Because different crops have different attributes, it is not rigorous to use the total output value or total yield of all crops as the index. In this research, the index c8 (yield per unit area of wheat), c9 (yield per unit area of maize), and c10 (yield per unit area of vegetables) were used to represent the food productivities of the arable land. The data of these indexes were calculated from data on gross yields and planting areas given in the Statistical Yearbook of Jinan City.

3.2. Calculation of the Weights

In a comprehensive evaluation, the weight reflects the significant performance of each basic indicator or dimension criterion [30,31]. Major factors act on sustainability in different degrees. Thus, each major factor must be assigned with a weight. Selecting the appropriate weight determination method is the key step of the comprehensive evaluation work. Weight calculation methods primarily include objective weighting methods and subjective evaluation methods. Objective weighting methods are represented by the entropy method [32] and the coefficient of variation (CV) method [33]. Subjective weighting methods are represented by the analytic hierarchy process (AHP) method [34] and the Delphi method [35].
Each type of those methods has its own advantages and disadvantages. Objective weighting methods depend on objective data completely, whereas subjective weighting methods are influenced by human-subjective experience inevitably. The combination weighting method [36] combines the subjective weighting method with the objective method based on game theory, so that the weight can reflect the subjective and objective information simultaneously. The central idea is to minimize the deviations of different weights calculated by different methods. In this study, subjective weights were determined by the G1 method, whereas objective weights were acquired by the entropy method.

3.2.1. The G1 Method

The G1 method (order relation analysis) is an improved method based on the analytic hierarchy process (AHP) [21]. This method does not require a consistency check, so a large amount of calculations can be avoided. The method consists of three steps:
Firstly, the order of the evaluation indexes has to be determined
c 1 c 2 c n
The number of the evaluation indexes is n. “ ” indicates that the index on the left side of the symbol is more important.
Secondly, the importance ratio assignments of adjacent evaluation indexes are determined.
w i 1 1 w i 1 = r i ( i = 2 , 3 , , n )
The reference assignments for r i are shown in Table 2.
Finally, the weights of indexes are calculated
w i 1 = ( 1 + i = 2 n k = i n r k ) 1
w i 1 1 = w i 1 · r i   ( i = n , n 1 , , 2 )
where w n 1 , w i 1 , and w i 1 represent the weights of the indexes cn, ci−1, and ci using the G1 method, respectively.

3.2.2. The Entropy Method

The entropy method is based on objective data. Therefore, it can avoid personal subjective interference. Supposing the number of evaluation objects is m and the number of evaluation indexes is n, the original data matrix can be expressed as: X = {xij}m×n (j = 1, 2, 3, …, m; I = 1, 2, 3, …, n).
In the process of the calculation, negative and extreme values cannot be used directly because of the concept of logarithm and entropy. Thus, the index data needs to be normalized in advance. The normalized indexes can be divided into positive type and negative type, as shown in Equations (11) and (12):
y i j = x i j min ( x i j ) max ( x i j ) min ( x i j )
y i j = max ( x i j ) x i j max ( x i j ) min ( x i j )
The entropy measure is calculated using the following equation:
E i = ln ( m ) 1 j = 1 m p i j × ln p i j
where p i j = y i j j = 1 m y i j , Ei is the entropy of the ith index. It is supposed that when p i j = 0 , ln p i j = 0 .
The weights of indexes are obtained as:
w i 2 = 1 E i n i = 1 n E i
where w i 2 denotes the weight of the index xi using entropy method.

3.2.3. Combination Weighting Rule

If L types of methods are used to assign weights to an index; a weight vector set can be constructed: wkT = {wk1, wk2, …, wkn}, k = 1,2, …, L. The linear combination of L weight vectors can be represented by Equation (15):
w = k = 1 L α k w k T
where w is a possible weight vector obtained through L types of methods, and α k is the coefficient of the kth method.
Searching for the most ideal weight vector can be aided by optimizing the coefficients α k in the Formula (15). Based on game theory [37,38], the goal of the optimization is to minimize the deviation between w and the various w k T :
min k = 1 L α k w k T w u 2 u = 1 , 2 , , L
The optimal first-order derivative conditions can be converted to the following linear system of equations:
[ w 1 · w 1 T w 1 · w 2 T w 1 · w L T w 2 · w 1 T w 2 · w 2 T w 2 · w L T w L · w 1 T w L · w 2 T w L · w L T ] [ α 1 α 2 α 3 α 4 ] = [ w 1 · w 1 T w 2 · w 2 T w L · w L T ]
α 1 , α 2 , , α L can be computed mathematically and normalized, as in Equation (18):
α k = α k / k = 1 L α k
The combination weights of indexes can be obtained as:
w = k = 1 L α k w k T
where w represents ( w 1 , w 2 , …, w i , …, w n ).

3.3. Matter–Element Model

The matter–element theory has been widely applied to mathematics and system theory [39]. The matter–element model consists of objects, characteristics, and values [40]. Therefore, agricultural sustainability as the evaluation object can be regarded as a composite whole with different characteristics and corresponding values. The attribute diversity of agricultural sustainability can be assumed as an incompatibility problem. The extension method aims to solve this type of problem by combining indexes with different attributes to deduce a comprehensive result [36].

3.3.1. Determining the Matter–Element

The matter–element is the basic element to describe objects (R), which is constituted of a given object (N), a characteristic (C), and the value of the characteristic (V). The matter–element is defined according to by Equation (20):
R = ( N , C , V ) = [ N c 1 v 1 c 2 v 2 c i v i c n v n ]   ( i = 1 , 2 , , n )

3.3.2. Determining the Classical Domain and the Controlled Domain Matter–Element Matrix

The classical domain matter–element matrix can be expressed as Equation (21):
R t = ( N t , C , V t i ) = [ N t c 1 a t 1 , b t 1 c 2 a t 2 , b t 2 c i a t i , b t i c n a t n , b t n ]
where t represents the tth class of agricultural sustainability. Rt is the matrix of the classical domain in the tth class. Nt is the agricultural sustainability in tth class. Vti is the value range of the ith index in the tth class, with the lower bound of ati and the upper bound of bti.
The controlled domain matter–element matrix can be expressed as Equation (22):
R p = ( N , C , V p i ) = [ N p c 1 a p 1 , b p 1 c 2 a p 2 , b p 2 c i a p i , b p i c n a p n , b p n ]
where Rp represents the matter–element of the controlled domain. Vpi is the value range of the ith index in all classes, with the lower bound of api and the upper bound of bpi.

3.3.3. Correlation Degree between Each Index and Each Class

The correlation degree determines whether one index belongs to a sustainability class and how much it belongs there. The value of the correlation degree was calculated according to the following functions.
K t ( x i j ) = { ρ ( x i j , V t i ) | V t i |   ,   x i j V t i ρ ( x i j , V t i ) ρ ( x i , V p i ) ρ ( x i , V t i )   ,   x i j V t i
where
ρ ( x i j , V t i ) = | x i j 1 2 ( a t i + b t i ) | 1 2 ( b t i a t i )
ρ ( x i j , V p i ) = | x i j 1 2 ( a p i + b p i ) | 1 2 ( b p i a p i )
where xij is the value of the ith index of the jth object. Kt(xij) is the value of correlation degree between the ith index of the jth object and the tth class. cij belongs to the tth class if Kt (xij) = max [Kt(xij)] (t = 1, 2, …, o).
Combined with the corresponding weights according to Section 3.2, the integrated correlation degree Kt (Nj) between Nj and the tth class is shown in Equation (26).
K t ( N j ) = i = 1 n w i K t ( x i j )   ( t = 1 , 2 , , o )
where Nj is the jth object and o represents the number of all the classes. In conclusion, Nj belongs to the tth class if Kt (Nj) = max [Kt (Nj)] (t = 1, 2, …, o).

4. Results

4.1. The Weights of Indexes

Subjective weights were calculated by the G1 method, as stated in Section 3.2.1. The order of the evaluation indexes is c4, c2, c7, c6, c1, c10, c8, c9, c5, and finally, c3. The importance ratio assignments of the adjacent evaluations are c4/c2 = 1.2; c2/c7 = 1; c7/c6 = 1; c6/c1 = 1.2; c1/c10 = 1; c10/c8 = 1; c8/c9 = 1.2; c9/c5 = 1; c5/c3 = 1. It should be noted that the G1 method is a subjective method, so that the order and order and ratios of the indexes have to be defined according to study site characteristics and research objectives on a case-by-case basis.
Objective weights were calculated by the entropy method, as stated in Section 3.2.2. According to the combination weighting rule, α 1 = 0.915 and α 2 = 0.085 . Subjective weights, objective weights, and synthetic weights of indexes are shown in Table 3.

4.2. The Classical Domain and the Controlled Domain

Agricultural sustainability classes reflect the degree of agricultural sustainability [41]. The degree of agricultural sustainability was divided into five classes: very low, low, moderate, high, and very high. Thus, as described in Section 3.3.2, the classical domain and the controlled domain are listed in Table 4. The classical domain of each index was determined by data from Jinan over the past 30 years. The controlled domain was determined by minimum or maximum values of the indexes of the research site.

4.3. Determination of the Sustainability Class

For each object to be evaluated, the correlation between each index and each evaluation class can be calculated by Equations (23)–(25). For example, Table 5 shows the correlation degree and evaluation class of each index in 2015. K1(c1) = −0.050, K2(c1) = −0.023, K3(c1) = 0.062, K4(c1) = −0.021, and K5(c1) = −0.038. Thus, c1 belongs to the moderate class because K3(c1) is the maximum of all of the classes.
The correlation degree between each index and evaluation class is shown in Figure 3. Evaluation results of each index from 2006 to 2015 are shown in Table 6. For simplicity, I, II, III, IV, and V represent very low, low, moderate, high, and very high, respectively.
Combined with the corresponding weights in Table 3, the integrated correlation degree between each year and evaluation class was calculated by Equation (26), as shown in Figure 4. Thus, the final evaluation results of every year were obtained as shown in Table 7.

5. Analysis and Discussion

The level of agricultural sustainability was promoted from class III to class IV in 2013 (Table 7), and the level was steady in class III before 2013. A higher degree of correlation to a given class indicates that the state in a given evaluation class is more stable, or that there is a greater tendency to convert to this class. Thus, it is worth noting that the integrated correlation degree of class IV was continuously increasing in 2013–2015 (Figure 2). It indicates that the state of agricultural sustainability in class IV was more and more stable. In the meantime, in regard to the other classes, the integrated correlation degree of class V was also growing, whereas the integrated correlation degree of class III was constant. Cumulatively, reflecting that agricultural sustainability was in a comparatively positive trend.
The above-mentioned analysis path can be extended to analyze the sustainable development of each index. The index c1 (agricultural blue water proportion) was switching back and forth between class II and class III. Besides, there was no obvious upward trend for the correlation degree of class IV. The evaluation result of c3 (irrigation proportion of arable land) was one class higher than that of c1. The trends of those two indexes were similar. This shows that the dependencies of agriculture on irrigation and irrigation intensity have not been effectively alleviated in recent years. The index c2 (water-use efficiency) ranked class IV in all the years and the correlation degree to class V kept increasing. This shows that water-use efficiency was improved. The improvement of water-use efficiency benefited from the renovation of irrigation infrastructure and the optimization of irrigation methods (i.e., sprinkler or drip irrigation) in the Chenmengquan irrigation district. The index c4 (energy utilization for irrigation amount per arable land) experienced a leap from class I to class IV in 10 years, which can be attributed to the increase of the groundwater table caused by the regional groundwater protection by the administrative department [42]. The indexes c5 (agricultural machinery power per arable land) and c6 (fertilizer utilization amount per arable land) were both in class III all the time and their correlation degrees to class III were significantly higher than those to other classes. This suggests that the sustainability classes of c5 and c6 will not change in the short term. The amount of fertilizer utilization did not decrease effectively. Except for 2006, c7 (pesticide utilization amount per arable land) was always in class II and the correlations to the two adjacent classes were still much lower than class II. It performed the worst of all the indexes. Regarding the use of pesticide, adequate attention should be paid by the administrative department. In the period 2010–2015, c8 (yield per unit area of wheat) and c10 (yield per unit area of vegetable) were both in class IV, and both of their correlation degrees with class IV, were getting close to the correlation degrees of class V. Those results showed that the productivity of wheat and vegetables was gradually increasing. However, the productivity of maize was getting worse, which can be demonstrated by the fact that c9 (yield per unit area of maize) dropped down from class IV to class III after 2012.
In addition, the sensitivity of the indexes was analyzed in this research. The simplest form of sensitivity analysis includes varying one index by a given amount at a time in the model and examining the impact of such change on the result of the model [43]. Consider that agricultural sustainability consists of multiple indexes which have different attributes. Each index value was given certain changes (−5%, −10%, +5%, and +10%) in double-directions [44]. The evaluation result of agricultural sustainability in 2015 did not change with the variations of the indexes. However, the evaluation results of c2, c8, c9, and c10 were changed, as shown in Table 8. Compared with other indexes, these four indexes had higher sensitivities. The variations of c8 (yield per unit area of wheat), c9 (yield per unit area of maize), c10 (yield per unit area of vegetable), and c2 (water-use efficiency) had a strong influence on the evaluation results of corresponding indexes.
It is worth noting that the proposed procedure cannot be generalized to other contexts and should be considered site-specific. The indexes need to be selected according to the characteristics of the agriculture in other research regions.

6. Conclusions

This paper proposes an index-based approach to evaluate agricultural sustainability in a given area based on the concept of the WEFN. The index system was established in the dimensions of agricultural water use, the intensity of agricultural energy use, and food productivity. The index system can not only reflect the agricultural sustainability comprehensively, but also avoid the influence of scale variation. In this premise, the evaluation approach integrated the combination weighting method and the matter–element model. The matter–element model was applied to quantify the fuzzy concept of agricultural sustainability. The combination weighting method gave each index a more scientific weight. The methodology was implemented in a case research, on an irrigation district in Jinan City of China, which is characterized by highly intensive agricultural practices. The results show the evaluation classes of agricultural sustainability in the study area. From 2006 to 2012, the level of agricultural sustainability was class III (moderate). From 2013 to 2015, the level of agricultural sustainability was class IV (high). At the same time, the results indicate the trend of the state change and the single index information. The agricultural sustainability was in a comparatively positive trend after 2012. Excessive use of pesticides was a shortcoming of agricultural practices in the study area. The results derived from the sensitivity analysis indicated that c2 (water-use efficiency), c8 (yield per unit area of wheat), c9 (yield per unit area of maize), and c10 (yield per unit area of vegetable) had higher sensitivities. The research results could be applied both for scientific and decision-making purposes. The methodology is suitable to evaluate agricultural sustainability.

Author Contributions

C.L. conceived the original research and drafted the manuscript; Z.Z. guided the execution of the research; S.L. and B.F. assisted in developing the research design; Q.L. contributed to the collection and processing of data; J.T. contributed to editing the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agathokleous, E.; Calabrese, E.J. Hormesis can enhance agricultural sustainability in a changing world. Global Food Sec. 2019, 20, 150–155. [Google Scholar] [CrossRef]
  2. FAO. Climate-Smart Agriculture Sourcebook. Available online: http://www.fao.org/3/i3325e/i3325e.pdf (accessed on 22 September 2019).
  3. FAO. The State of World’s Land and Water Resources for Food and Agriculture. Available online: http://www.fao.org/docrep/017/i1688e/i1688e.pdf (accessed on 22 September 2019).
  4. Xue, J.; Huo, Z.; Huang, Q.; Wang, F.; Boll, J.; Huang, G.; Qu, Z. Assessing sustainability of agricultural water saving in an arid area with shallow groundwater. Irrig. Drain. 2019, 68, 205–217. [Google Scholar] [CrossRef]
  5. Byomkesh, T.; Blay-Palmer, A. Comparison of Methods to Assess Agricultural Sustainability. In Sustainable Agriculture Reviews 13; Springer: Cham, Switzerland, 2017. [Google Scholar]
  6. Gomez-Limon, J.A.; Sanchez-Fernandez, G. Empirical evaluation of agricultural sustainability using composite indicators. Ecol. Econ. 2010, 69, 1062–1075. [Google Scholar] [CrossRef]
  7. Kamali, F.P.; Borges, J.A.R.; Meuwissen, M.P.M.; de Boer, I.J.M.; Lansink, A. Sustainability assessment of agricultural systems: The validity of expert opinion and robustness of a multi-criteria analysis. Agric. Syst. 2017, 157, 118–128. [Google Scholar] [CrossRef]
  8. Talukder, B.; Blay-Palmer, A.; Hipel, K.; Van Loon, G. Elimination Method of Multi-Criteria Decision Analysis (MCDA): A Simple Methodological Approach for Assessing Agricultural Sustainability. Sustainability 2017, 9, 287. [Google Scholar] [CrossRef]
  9. Bertocchi, M.; Demartini, E.; Marescotti, M.E. Ranking farms using quantitative indicators of sustainability: The 4Agro method. In 2nd International Symposium New Metropolitan Perspectives—Strategic Planning, Spatial Planning, Economic Programs and Decision Support Tools, through the Implementation of Horizon/Europe 2020; Calabro, F., Della Spina, L., Eds.; Elsevier Science Bv: Amsterdam, The Netherlands, 2016; pp. 726–732. [Google Scholar]
  10. Talukder, B.; Hipel, K.W. The PROMETHEE Framework for Comparing the Sustainability of Agricultural Systems. Resources 2018, 7, 74. [Google Scholar] [CrossRef]
  11. Van Cauwenbergh, N.; Biala, K.; Bielders, C.; Brouckaert, V.; Franchois, L.; Cidad, V.G.; Hermy, M.; Mathijs, E.; Muys, B.; Reijnders, J.; et al. SAFE—A hierarchical framework for assessing the sustainability of agricultural systems. Agric. Ecosyst. Environ. 2007, 120, 229–242. [Google Scholar] [CrossRef]
  12. Wohlfender-Bühler, D.; Feusthuber, E.; Wäger, R.; Mann, S.; Aubry, S.J. Genetically modified crops in Switzerland: Implications for agrosystem sustainability evidenced by multi-criteria model. Agron. Sustain. Dev. 2016, 36, 33. [Google Scholar] [CrossRef]
  13. Von Wiren-Lehr, S. Sustainability in agriculture—An evaluation of principal goal-oriented concepts to close the gap between theory and practice. Agric. Ecosyst. Environ. 2001, 84, 115–129. [Google Scholar] [CrossRef]
  14. Hoff, H. Understanding the Nexus: Background Paper for the BONN2011 Nexus Conference; SEI: Stockholm, Sweden, 2011; 51p. [Google Scholar]
  15. Li, M.; Fu, Q.; Singh, V.P.; Ji, Y.; Liu, D.; Zhang, C.; Li, T. An optimal modelling approach for managing agricultural water-energy-food nexus under uncertainty. Sci. Total Environ. 2019, 651, 1416–1434. [Google Scholar] [CrossRef]
  16. FAO. Walking the Nexus Talk: Assessing the Water-Energy-Food Nexus. Available online: http://www.fao.org/3/a-i3959e.pdf (accessed on 22 September 2019).
  17. Smajgl, A.; Ward, J.; Pluschke, L. The water-food-energy Nexus—Realising a new paradigm. J. Hydrol. 2016, 533, 533–540. [Google Scholar] [CrossRef]
  18. Han, H.; Li, H.; Zhang, K. Urban Water Ecosystem Health Evaluation Based on the Improved Fuzzy Matter-Element Extension Assessment Model: Case Study from Zhengzhou City, China. Math. Prob. Eng. 2019, 2019, 1–14. [Google Scholar] [CrossRef]
  19. Gebler, D.; Wiegleb, G.; Szoszkiewicz, K. Integrating river hydromorphology and water quality into ecological status modelling by artificial neural networks. Water Res. 2018, 139, 395–405. [Google Scholar] [CrossRef] [PubMed]
  20. Sun, D.; Wu, J.; Zhang, F.; Su, W.; Hui, H. Evaluating Water Resource Security in Karst Areas Using DPSIRM Modeling, Gray Correlation, and Matter-Element Analysis. Sustainability 2018, 10, 3934. [Google Scholar] [CrossRef]
  21. Qian, C.; Zhang, M.; Chen, Y.; Wang, R. A Quantitative Judgement Method for Safety Admittance of Facilities in Chemical Industrial Parks based on G1-Variation Coefficient Method. Procedia Eng. 2014, 84, 223–232. [Google Scholar] [CrossRef] [Green Version]
  22. Shemshadi, A.; Shirazi, H.; Toreihi, M.; Tarokh, M.J. A fuzzy VIKOR method for supplier selection based on entropy measure for objective weighting. Expert Syst. Applic. 2011, 38, 12160–12167. [Google Scholar] [CrossRef]
  23. Dong, L.; Gu, X.; Wu, X.; Liao, H. An improved MULTIMOORA method with combined weights and its application in assessing the innovative ability of universities. Expert Syst. 2019, 36, e12362. [Google Scholar] [CrossRef]
  24. Wen, C.; Dong, W.; Meng, Y.; Li, C.; Zhang, Q. Application of a loose coupling model for assessing the impact of land-cover changes on groundwater recharge in the Jinan spring area, China. Environ. Earth Sci. 2019, 78, 13. [Google Scholar] [CrossRef]
  25. Cao, X.; Ren, J.; Wu, M.; Guo, X.; Wang, Z.; Wang, W. Effective use rate of generalized water resources assessment and to improve agricultural water use efficiency evaluation index system. Ecol. Indic. 2018, 86, 58–66. [Google Scholar] [CrossRef]
  26. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232. [Google Scholar] [CrossRef]
  27. Pachauri, S.; Spreng, D. Direct and indirect energy requirements of households in India. Energy Policy 2002, 30, 511–523. [Google Scholar] [CrossRef]
  28. De Vito, R.; Portoghese, I.; Pagano, A.; Fratino, U.; Vurro, M. An index-based approach for the sustainability assessment of irrigation practice based on the water-energy-food nexus framework. Adv. Water Resour. 2017, 110, 423–436. [Google Scholar] [CrossRef]
  29. Daccache, A.; Ciurana, J.S.; Rodriguez Diaz, J.A.; Knox, J.W. Water and energy footprint of irrigated agriculture in the Mediterranean region. Environ. Res. Lett. 2014, 9, 124014. [Google Scholar] [CrossRef]
  30. Wang, M.; Tang, D.; Bai, Y.; Xia, Z. A compound cloud model for harmoniousness assessment of water allocation. Environ. Earth Sci. 2016, 75, 11. [Google Scholar] [CrossRef]
  31. Talukder, B.; Hipel, K.W.; VanLoon, G.W. Developing Composite Indicators for Agricultural Sustainability Assessment: Effect of Normalization and Aggregation Techniques. Resources 2017, 6, 66. [Google Scholar] [CrossRef]
  32. Islam, A.R.M.T.; Ahmed, N.; Bodrud-Doza, M.; Chu, R. Characterizing groundwater quality ranks for drinking purposes in Sylhet district, Bangladesh, using entropy method, spatial autocorrelation index, and geostatistics. Environ. Sci.Pollut. Res. 2017, 24, 26350–26374. [Google Scholar] [CrossRef] [PubMed]
  33. Busing, F.M.T.A.; Groenen, P.J.K.; Heiser, W.J. Avoiding degeneracy in multidimensional unfolding by penalizing on the coefficient of variation. Psychometrika 2005, 70, 71–98. [Google Scholar] [CrossRef]
  34. Pourghasemi, H.R.; Pradhan, B.; Gokceoglu, C. Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran. Nat. Hazards 2012, 63, 965–996. [Google Scholar] [CrossRef]
  35. Okoli, C.; Pawlowski, S.D. The Delphi method as a research tool: An example, design considerations and applications. Inf. Manag. 2004, 42, 15–29. [Google Scholar] [CrossRef]
  36. Zhang, W.; Chen, J.-P.; Wang, Q.; An, Y.; Qian, X.; Xiang, L.; He, L. Susceptibility analysis of large-scale debris flows based on combination weighting and extension methods. Nat. Hazards 2013, 66, 1073–1100. [Google Scholar] [CrossRef]
  37. Szeszlér, D. Hitting a path: A generalization of weighted connectivity via game theory. J. Combin. Optim. 2019, 38, 72–85. [Google Scholar] [CrossRef]
  38. Madani, K. Game theory and water resources. J. Hydrol. 2010, 381, 225–238. [Google Scholar] [CrossRef]
  39. Liu, D.; Zou, Z. Water quality evaluation based on improved fuzzy matter-element method. J. Environ. Sci. 2012, 24, 1210–1216. [Google Scholar] [CrossRef]
  40. He, Y.-X.; Dai, A.-Y.; Zhu, J.; He, H.-Y.; Li, F. Risk assessment of urban network planning in china based on the matter-element model and extension analysis. Int. J. Electric. Power Energy Syst. 2011, 33, 775–782. [Google Scholar] [CrossRef]
  41. Gong, J.; Liu, Y.; Chen, W. Land suitability evaluation for development using a matter-element model: A case study in Zengcheng, Guangzhou, China. Land Use Policy 2012, 29, 464–472. [Google Scholar] [CrossRef]
  42. Yang, R.; Cui, B.; Zhao, H.; Lei, X. An Integrated Act for Groundwater Protection of Jinan City, China. Procedia Environ. Sci. 2010, 2, 1745–1749. [Google Scholar] [CrossRef] [Green Version]
  43. Nemati, H.; Shokri, M.R.; Ramezanpour, Z.; Ebrahimi Pour, G.H.; Muxika, I.; Borja, Á. Sensitivity of indicators matters when using aggregation methods to assess marine environmental status. Mar. Poll. Bull. 2018, 128, 234–239. [Google Scholar] [CrossRef]
  44. Alvarez, M.C.; Franco, A.; Pérez-Domínguez, R.; Elliott, M. Sensitivity analysis to explore responsiveness and dynamic range of multi-metric fish-based indices for assessing the ecological status of estuaries and lagoons. Hydrobiologia 2013, 704, 347–362. [Google Scholar] [CrossRef]
Figure 1. A location map of the Chenmengquan irrigation district.
Figure 1. A location map of the Chenmengquan irrigation district.
Sustainability 11 05350 g001
Figure 2. The steps of the methodology.
Figure 2. The steps of the methodology.
Sustainability 11 05350 g002
Figure 3. The correlation degree between each index and each evaluation class. (a) The correlation degree between c1 and each class; (b) the correlation degree between c2 and each class; (c) the correlation degree between c3 and each class; (d) the correlation degree between c4 and each class; (e) the correlation degree between c5 and each class; (f) the correlation degree between c6 and each class; (g) the correlation degree between c7 and each class; (h) the correlation degree between c8 and each class; (i) the correlation degree between c9 and each class; (j) the correlation degree between c10 and each class.
Figure 3. The correlation degree between each index and each evaluation class. (a) The correlation degree between c1 and each class; (b) the correlation degree between c2 and each class; (c) the correlation degree between c3 and each class; (d) the correlation degree between c4 and each class; (e) the correlation degree between c5 and each class; (f) the correlation degree between c6 and each class; (g) the correlation degree between c7 and each class; (h) the correlation degree between c8 and each class; (i) the correlation degree between c9 and each class; (j) the correlation degree between c10 and each class.
Sustainability 11 05350 g003
Figure 4. The integrated correlation degree between each year and evaluation class.
Figure 4. The integrated correlation degree between each year and evaluation class.
Sustainability 11 05350 g004
Table 1. The index system based on water–energy–food nexus.
Table 1. The index system based on water–energy–food nexus.
SectorIndex (unit)
Water c1 Agricultural blue water proportion
c2 Water-use efficiency
c3 Irrigation proportion of arable land
Energyc4 Energy utilization for irrigation amount per arable land (kWh/ha)
c5 Agricultural machinery power per arable land (kW/ha)
c6 Fertilizer utilization amount per arable land(kg/ha)
c7 Pesticide utilization amount per arable land(kg/ha)
Foodc8 Yield per unit area of Wheat(t/ha)
c9 Yield per unit area of Maize(t/ha)
c10 Yield per unit area of Vegetable(t/ha)
Table 2. The importance ratio assignments of adjacent evaluation indexes.
Table 2. The importance ratio assignments of adjacent evaluation indexes.
riInstruction
1.0The index ci−1 and index ci are equally important
1.2The index ci−1 is slightly more important than index ci
1.4The index ci−1 is fairly more important than index ci
1.6The index ci−1 is strongly more important than index ci
2.0The index ci−1 is extremely more important than index ci
Table 3. The weights of indexes.
Table 3. The weights of indexes.
Indexc1c2c3c4c5c6c7c8c9c10
wi10.0950.1140.0790.1370.0790.1140.1140.0950.0790.095
wi20.1330.0840.1420.1930.0550.0720.0050.1530.0960.068
wi0.130.0870.1370.1880.0570.0750.0140.1480.0940.07
Table 4. The classical domain and the controlled domain.
Table 4. The classical domain and the controlled domain.
IndexThe Classical DomainThe Controlled Domain
Very LowLowModerateHighVery High
c1[0.8, 1][0.6, 0.8][0.4, 0.6][0.2, 0.4][0, 0.2][0, 1]
c2[0, 0.2][0.2, 0.4][0.4, 0.6][0.6, 0.8][0.8, 1][0, 1]
c3[0, 0.2][0.2, 0.4][0.4, 0.6][0.6, 0.8][0.8, 1][0, 1]
c4[1.6, 2][1.2, 1.6][0.8, 1.2][0.4, 0.8][0, 0.4][0, 2]
c5[25, 30][20, 25][15, 20][10, 15][0, 10][0, 30]
c6[800, 1000][600, 800][400, 600][200, 400][0, 200][0, 1000]
c7[16, 20][12, 16][8, 12][4, 8][0, 4][0, 20]
c8[2, 3][3, 4][4, 5][5, 6][6, 7][2, 7]
c9[3, 4][4, 5][5, 6][6, 7][7, 8][3, 8]
c10[30, 42][42, 54][54, 66][66, 78][78, 90][30, 90]
Table 5. The correlation degree and evaluation class of each index in 2015.
Table 5. The correlation degree and evaluation class of each index in 2015.
IndexVery LowLowModerateHighVery HighCLASSES
c1−0.050−0.0230.062−0.021−0.038moderate
c2−0.038−0.0220.022−0.009−0.022moderate
c3−0.090−0.074−0.0430.050−0.010high
c4−0.111−0.085−0.0340.069−0.024high
c5−0.020−0.0060.015−0.014−0.017moderate
c6−0.028−0.0110.032−0.014−0.024moderate
c7−0.0040.005−0.003−0.007−0.007low
c8−0.100−0.085−0.0530.042−0.008high
c9−0.035−0.0140.039−0.018−0.030moderate
c10−0.044−0.036−0.0190.032−0.006high
Table 6. The evaluation class of each index in 2006–2015.
Table 6. The evaluation class of each index in 2006–2015.
Index2006200720082009201020112012201320142015
c1
c2
c3
c4
c5
c6
c7
c8
c9
c10
Table 7. The integrated correlation degree and evaluation class of 2006–2015.
Table 7. The integrated correlation degree and evaluation class of 2006–2015.
ObjectVery LowLowModerateHighVery HighClasses
2006−0.374−0.270−0.105−0.172−0.293moderate
2007−0.310−0.348−0.112−0.197−0.312moderate
2008−0.447−0.2460.130−0.026−0.250moderate
2009−0.460−0.2360.067−0.014−0.246moderate
2010−0.427−0.2270.189−0.091−0.265moderate
2011−0.420−0.2020.205−0.106−0.274moderate
2012−0.477−0.280.03630.032−0.227moderate
2013−0.480−0.2770.0210.049−0.226high
2014−0.490−0.2950.0130.062−0.217high
2015−0.521−0.3510.0180.110−0.186high
Table 8. Evaluation results with index variations in 2015.
Table 8. Evaluation results with index variations in 2015.
IndexVariation Range
−5%−10%+5%+10%0
c2IIIIIIIIIIVIII
c8IVIVIV
c9IIIIIIIIIIIIII
c10IVIIIIVIV

Share and Cite

MDPI and ACS Style

Liu, C.; Zhang, Z.; Liu, S.; Liu, Q.; Feng, B.; Tanzer, J. Evaluating Agricultural Sustainability Based on the Water–Energy–Food Nexus in the Chenmengquan Irrigation District of China. Sustainability 2019, 11, 5350. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195350

AMA Style

Liu C, Zhang Z, Liu S, Liu Q, Feng B, Tanzer J. Evaluating Agricultural Sustainability Based on the Water–Energy–Food Nexus in the Chenmengquan Irrigation District of China. Sustainability. 2019; 11(19):5350. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195350

Chicago/Turabian Style

Liu, Chang, Zhanyu Zhang, Shuya Liu, Qiaoyuan Liu, Baoping Feng, and Julia Tanzer. 2019. "Evaluating Agricultural Sustainability Based on the Water–Energy–Food Nexus in the Chenmengquan Irrigation District of China" Sustainability 11, no. 19: 5350. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop