Next Article in Journal
Using Knowledge Granularity Entropy to Measure Eco-Environmental Impacts of Land Cover Changes in ASEAN from 2001 to 2020
Next Article in Special Issue
Study of the Spatial Spillover Effects of the Efficiency of Agricultural Product Circulation in Provinces along the Belt and Road under the Green Total Factor Productivity Framework
Previous Article in Journal
Assessment of Variability in Hydrological Droughts Using the Improved Innovative Trend Analysis Method
Previous Article in Special Issue
Research on Agricultural Total Production Factors in the Bohai Bay Area Based on Agricultural Data from 2011 to 2020
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Dynamics of Agricultural Sustainability Assessment: A Study across 30 Chinese Provinces

1
Management College, Ocean University of China, Qingdao 266100, China
2
School of Business and Tourism Management, Yunnan University, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(11), 9066; https://doi.org/10.3390/su15119066
Submission received: 4 May 2023 / Revised: 30 May 2023 / Accepted: 31 May 2023 / Published: 3 June 2023

Abstract

:
Agricultural sustainability is crucial for ensuring food security, promoting economic development, maintaining ecological balance, and achieving sustainable development goals. In this study, based on relevant theories of agricultural sustainability, an analytical framework is constructed for agricultural sustainability encompassing economic, resource, environmental, and social dimensions. The Analytic Network Process (ANP) method is employed to determine indicator weights and assess the spatiotemporal changes in agricultural sustainability levels across Chinese provinces. The findings reveal that environmental quality is the primary dimension for assessing agricultural sustainability, and the significance of the rural social development dimension is continuously increasing. Although the sustainability levels have significantly improved in various regions of China, there remain issues of development imbalance and instability. In conclusion, this paper offers a comprehensive understanding of the spatiotemporal changes in agricultural sustainability across Chinese provinces, providing valuable insights for policymakers and researchers.

1. Introduction

Agriculture is a fundamental pillar of human civilization, shaping the way societies have evolved and interacted throughout history. In the contemporary world, agricultural sustainability development has emerged as a crucial aspect of global progress, significantly impacting food security, economic growth, and ecological balance [1]. As the global population continues to rise and the ramifications of climate change become increasingly evident, the intensified demand for agriculture leads to mounting pressures on natural resources and ecosystems [2]. Consequently, achieving sustainability in global agriculture is essential to meet current human development needs and ensure long-term sustainable growth [3].
Sustainable agriculture refers to a mode of production that aims to meet current and future human needs while increasing resource efficiency, protecting natural resources and ecosystems, and promoting socio-economic well-being [4]. On a macro level, sustainable agriculture contributes to global food security and poverty eradication [5]. Moreover, sustainable agricultural practices can help mitigate climate change by reducing greenhouse gas emissions, enhancing carbon sequestration, and promoting biodiversity and ecosystem protection [6,7]. As a result, improving agricultural sustainability can directly advance several targets of the United Nations’ 2030 Sustainable Development Goals (SDGs), primarily involving SDG1 “No Poverty”, SDG2 “Zero Hunger”, SDG4 “Quality Education”, SDG5 “Gender Equality”, SDG8 “Decent Work and Economic Growth”, SDG12 “Responsible Consumption and Production”, SDG13 “Climate Action”, SDG14 “Life Below Water”, and SDG15 “Life on Land” [8]. Given the multifaceted nature of agricultural sustainability and its intersection with ecological, environmental, social, and resource aspects, a comprehensive evaluation of agricultural sustainability is paramount for informing policy decisions and promoting sustainable development worldwide.
Over the past few decades, various efforts have been made to assess agricultural sustainability, ranging from simplistic indicator-based approaches to intricate integrated assessment models [9,10,11,12]. On the indicator level, initial evaluations of agricultural sustainability focus primarily on the ecological impacts of agricultural production [13], gradually expanding to encompass economic, environmental, and social aspects [14,15], with some studies emphasizing the role of resource utilization [16]. At the model level, numerous assessment models have been developed for gauging agricultural sustainability, such as the Sustainability Assessment of Farming and the Environment (SAFE) [17], Sustainability Assessment (MOTIFS) [18], Sustainability Assessment for Food and Agriculture (SAFA) system [19], Agricultural Environmental Indicators (AEI) [20], the Multi-Criteria Decision Analysis (MCDA)-based method [21,22] and some complex assessment models based on the generalized combination rule [22] and deep learning [23].
In the literature, indicators or models for agricultural sustainability assessment may differ due to variations in time, space, and theoretical considerations, with most involving the standardization of indicator data and weight calculation [6,24,25,26,27]. Previous studies have mostly employed evaluation methods such as the Entropy weight method [28], Principal Component Analysis (PCA) [29], Dematel [30], and Analytic Hierarchy Process (AHP) [31] for calculating indicator weights in agricultural sustainability assessments. While these studies have significantly contributed to understanding agricultural sustainability, many of these approaches fail to capture the full complexity and interconnectedness, ignoring the interactions among ecological, resource, environmental, and social factors [32,33]. Addressing these limitations, this study constructs a comprehensive agricultural sustainability evaluation system using the Analytic Network Process (ANP) methodology. By acknowledging the interconnected nature of agricultural systems, this approach provides a more integrated and detailed understanding of agricultural sustainability, offering a more robust framework for its assessment.
The Analytic Network Process (ANP) is a multi-criteria decision making method developed by Saaty based on the foundation of the Analytic Hierarchy Process (AHP), allowing for the comprehensive evaluation of complex systems with multiple interdependent dimensions [34]. The primary feature of ANP lies in providing a framework that considers the interrelationships among all evaluation indicators (criteria) within and across clusters [32]. ANP has been employed by numerous researchers to address various decision-making problems, proving its suitability for constructing models with evaluation criteria and dimensions encompassing intricate interactions [35]. Consequently, it has been widely applied in risk assessment [36], environmental management [37], and various evaluation problems [38]. This study uses the ANP method to thoroughly explore the interplay among economic, resource, environmental, and social dimensions while evaluating agricultural sustainability levels.
By integrating the ANP methodology, this study seeks to develop a holistic agricultural sustainability assessment framework that captures the complexity and interconnections among various dimensions of agricultural sustainability, offering a more robust and comprehensive evaluation tool for policymakers and researchers. Using China as a case study, we investigate the spatiotemporal dynamics of agricultural sustainability, providing valuable insights for those seeking to advance sustainable agriculture worldwide. Our findings contribute to a broader understanding of agricultural sustainability assessment and management, emphasizing the importance of a holistic approach that considers the interactions among economic, resource, environmental, and social dimensions.
The remainder of this study is organized as follows. The next section describes the methods, starting with the research area and data sources, then the construction of the evaluation index system, determination of index weights using the ANP method, and calculation of agricultural sustainability levels. In Section 3, we analyze the results, examining index weights and conducting temporal and spatial comparative analyses of agricultural sustainability. Section 4 focuses on the discussion, comparing our findings with previous studies, exploring the implications of our results for agricultural policy and planning, and identifying limitations and potential improvements for future research. The conclusion summarizes our findings and offers policy recommendations for enhancing agricultural sustainability within sustainable development goals and global ecological, environmental, and resource management.

2. Materials and Methods

2.1. Research Area and Data Source

Situated in East Asia, China spans from 73°40  to 135°2  east longitude and 3°52  to 53°33  north latitude, encompassing a diverse array of climates and abundant natural resources (Figure 1). China’s agriculture provides a supply of staple foods for approximately 22% of the global population, establishing itself as one of the world’s largest food producers and consumers [39]. Investigating agricultural sustainability in China offers valuable insights that contribute to advancing research and practice in global sustainable agriculture.
The primary data used in this paper are the agricultural and rural-related data of 30 provinces, autonomous regions, and municipalities in mainland China from 2005 to 2020. Due to data availability and completeness, this study does not include Shanghai, Taiwan, Macau, and Hong Kong. All data are from the China Statistical Yearbook and the China Rural Statistical Yearbook released by the National and Local Bureau of Statistics of China. All the data are carried out at the provincial level because all indicators at the provincial level are the complete data, which can ensure this paper’s scientific and comprehensive results.

2.2. Construction of Evaluation Index System

In order to construct a comprehensive agricultural sustainability evaluation system, we integrated the multiple connotations of agricultural sustainability into different discourses. Globally, politics prioritize sustainable agriculture techniques and economic growth [40], while academia examines resource use and environmental preservation [41]. Notably, social welfare is increasingly recognized in both domains [42]. Based on the above-mentioned theoretical analysis, this study synthesizes the focal points of both domains, developing a comprehensive agricultural sustainability assessment framework encompassing four dimensions: agricultural economy, resource utilization, environmental quality, and rural society. The details are as follows.
In the agricultural economy dimension, specifically, five key indicators are considered: per labor gross agricultural value (e1), per labor grain output (e2), agricultural mechanization level (e3), agricultural electrification level (e4), and land productivity (e5) [43,44]. These indicators reflect the agricultural sector’s overall economic vitality and directly influence the potential for sustainable growth. Higher values in these indicators suggest a more robust agricultural economy conducive to long-term sustainability. Enhanced agricultural electrification (e4), for example, fosters the adoption of advanced machinery and irrigation systems, catalyzing production efficiency, economic growth, and, consequently, agricultural sustainability [45].
Transitioning to the resource utilization dimension, we focus on four indicators: per capita arable land area (e6), per capita water resource availability (e7), rural power supply level (e8), and effective irrigation rate (e9) [43,46]. Efficient resource utilization efficiency is the basis for achieving sustainable agricultural development, reducing resource depletion and achieving environmental protection. Higher scores in these indicators denote more sustainable resource use and bolster agricultural activities without compromising the needs of future generations.
For the environmental quality dimension, we examine five indicators: fertilizer use intensity (e10), pesticide use intensity (e11), forest coverage rate (e12), air quality index (e13), and soil and water conservation area (e14) [6,46]. These indicators illustrate the ecological impact of agricultural practices and underscore the necessity of environmentally friendly approaches. Reduced values in e10 and e11 and increased values in e12, e13, and e14 imply enhanced environmental quality—an indispensable component for the long-term viability of agriculture and rural ecosystems.
Finally, the rural society dimension incorporates four indicators: per capita disposable income of rural residents (e15), Engel’s coefficient of rural residents (e16), rural healthcare level (e17), and education level of rural residents (e18) [44,47]. These epitomize the overall welfare of the rural populace, integral to sustainable agriculture. High disposable income and education levels in rural areas foster investment in and comprehension of sustainable agricultural practice, while robust rural healthcare sustains a capable workforce [26]. A lower Engel coefficient indicates improved rural living standards and a transition toward sustainable lifestyles. Collectively, these indicators buttress the long-term sustainability of agriculture.
Based on the above analysis and guided by systematicity, scientificity, and data availability principles, this paper constructs a comprehensive agricultural sustainability evaluation system with 18 secondary indicators covering the four dimensions of the economic, resource, environmental, and social aspects of sustainable agricultural development, as shown in Table 1.

2.3. Index Weights Determination Using ANP

In order to take into account the complex interactions between the dimensions of the agricultural sustainability evaluation indicators, this study uses the Analytic Network Process (ANP) method to evaluate indicator weights. The calculation process and basic steps are as follows.

2.3.1. Network Model Construction

The ANP network structure is divided into control and network layers. The control layer determines the system’s general objective and the corresponding criterion, while the construction of the relationship between the network layer elements is the core of the ANP method. In this study, the control layer has only one element, agricultural sustainability, and thus was omitted. The network layer is constructed using the previous agricultural sustainability evaluation index system. The clusters in the network layer and the nodes under each cluster correspond to the four dimensions and indicators of the agricultural economy, resource use, environmental quality, and rural society in the evaluation system (Table 1), respectively.
In the network layer, the nodes are interdependent and dominated by each other, and the nodes and clusters are not internally independent, forming a network structure of mutual influence [48]. The dependency relationships may not also occur between nodes from different clusters but between those included in the same cluster. Within the agricultural sustainability evaluation system, for example, the node “e1 (Gross agricultural product of rural population)” from the “Agricultural Economy” cluster is both affected by “e4 (Agricultural electrification level)” from the same cluster and “e9 (Effective irrigation rate)” from the “Resource utilization” cluster, and so on.
Based on the above analysis, a network model of agricultural sustainability evaluation was constructed in this study, as shown in Figure 2. In this network model, the influence relationships between different nodes are indicated by arrows: arrows pointing to other clusters indicate that there are nodes under the original cluster that influence the nodes of other clusters, and arrows pointing to their self-clusters indicate that there are nodes under the same cluster that influence each other.

2.3.2. Calculation of Weights

Due to the complex relationship between the elements and the excessive calculation of ANP method, this paper uses Superdecisions software to calculate the weights in the agricultural sustainability evaluation index system. The specific process and principles are as follows.
  • step 1: conduct the pairwise comparisons.
    Expert opinions were sought to perform pairwise comparisons of the indicators within each cluster and between clusters, reflecting their relative importance in determining agricultural sustainability. The pairwise comparisons were quantified using a nine-point scale, where 1 represented equal importance and 9 indicated extreme importance. Inverses of these values were assigned for reciprocal comparisons [32,49]. These values were subsequently input into the SuperDecisions software for further computations.
  • step 2: calculate the unweighted supermatrix.
    Based on the scoring results in step 1, firstly, the maximum eigenvalues are calculated, and then the eigenvectors are normalized. 2. Please carefully check variable formatting (italic, bold, subscript, uppercase, etc.) throughout the manuscript to ensure the formatting is consistent and revise if needed.
    R j l R i 1 R i 2 R i n Normalized Eigenvector R i 1 e 11 e 12 e 1 n Ω i 1 j l R i 2 e 21 e 22 e 2 n Ω i 2 j l R i n e n 1 e n 2 e n n Ω i n j l
    Then, using the eigenvalue method, calculate the ranking vector ( Ω i 1 j l , Ω i 2 j l , , Ω i n j l ) and complete the consistency check.  W i j  is denoted as:
    W i j = Ω i 1 j 1 Ω i 1 j 2 Ω i 1 j n Ω i 2 j 1 Ω i 2 j 2 Ω i 2 j n Ω i n j 1 Ω i n j 2 Ω i n j n
    The column vector  W i j  represents the influence ranking of element  R i n i  in system  R i  on element  R j j  in system  R j . If the elements in  R i  have no influence on the elements in  R j , then  Ω i j = 0  ( i = 1 , 2 , , N j = 1 , 2 , , N ).
    Finally, the unweighted supermatrix of elemental interactions under the control level of the total system is obtained.
    W = r 11 r 12 r 1 n 1 r 21 r 22 r 2 n 2 r j 1 r j 2 r j n j r n 1 r n 2 r N N r 11 r 12 r 1 n r 21 r 22 r 2 n r i 1 r i 2 r i n r N 1 r N 2 r N N [    Ω 11       Ω 12         Ω 1 j     Ω 1 N       Ω 21       Ω 22         Ω 2 j     Ω 2 N                               Ω i 1       Ω i 2         Ω i j     Ω i N                               Ω N 1       Ω N 2         Ω N j     Ω N N    ]
  • step 3: calculate the weighted supermatrix.
    In the above supermatrix, all  W i j  matrices are based on  R j n  as the sub-criterion, and the ranking vectors are obtained by pairwise comparisons of elements in  R i . Although each column in  W i j  is column-normalized, W is not normalized. Therefore, it is necessary to compare the importance of each element group  R i  ( i = 1 , 2 , , N ) under the control layer for the sub-criterion  R j  ( j = 1 , 2 , , N ) and pass the consistency check.
    R i R 1 R 2 R N Normalized Eigenvector R 1 r 11 r 12 r 1 N a 1 j R 2 r 21 r 22 r 2 N a 2 j R N r N 1 r N 2 r N N a N j
    According to the relative weight value  A  of  R i , calculate the eigenvector to obtain the following weighted matrix  A :
    A = α 11 α 12 α 1 j α 1 n α 21 α 22 α 2 j α 2 n α i 1 α i 2 α i j α i n α N 1 α N 2 α N N α n n
    Weight the elements in the supermatrix to obtain  W ¯ = [ Ω i j ] n × n , where  w ¯ i j = a i j Ω i j , ( i = 1 , 2 , , N j = 1 , 2 , , N ). W is the weighted supermatrix, and its characteristic is that the sum of the columns is 1.
    W ¯ = a 11 Ω 11 a 12 Ω 12 a 1 j Ω 1 j a 1 n Ω 1 n a 21 Ω 21 a 22 Ω 22 a 2 j Ω 2 j a 2 n Ω 2 n a i 1 Ω i 1 a i 2 Ω i 2 a i j Ω i j a i n Ω i n a n 1 Ω n 1 a n 2 Ω n 2 a n j Ω n j a n n Ω n n
  • step 4: Calculate the weights of the elements
    By using the weighted supermatrix and the idea of normalization, the limit supermatrix  W  is obtained through consecutive multiplication, and its vector is the weight vector  W  of the element  R i n i :
    W = ( ω 11 , ω 12 , , ω n n ) Γ
    Here,  Γ  represents the number of matrix multiplications performed until the convergence criterion is satisfied. Thus, the weight vectors of each element  R i j  are as follows.
    W i = ( ω i 1 , ω i 2 , , ω i n ) Γ

2.4. Calculation of the Evaluation Value

2.4.1. Normalization of Indicator Data

Due to our agricultural sustainability evaluation system involving many aspects of economic resources, environment, and society, the indicators are scattered, some of which are negative. In order to eliminate the influence of direction and dimensionality, this paper adopts a normalization method to process the data. Specifically, we postulate the application of our agricultural sustainability indicator framework to assess the sustainability levels of T distinct regions. For the t-th region, the manifestation of agricultural sustainability pertaining to the n-th indicator, denoted as  e n , is represented by  A t n , where  t = 1 , 2 , , T , and  n = 1 , 2 , , 18 . The specific equations are as follows:
X t n = A t n min A 1 n , , A T n max A 1 n , , A T n min A 1 n , , A T n
X t n = max A 1 n , , A T n A t n max A 1 n , , A T n min A 1 n , , A T n
Here, Equation (9) is employed for the normalization of positive indicators, while Equation (10) is utilized for the normalization of negative indicators. It is important to note that the aforementioned normalization equations are exclusively applicable for spatial analysis, which is employed to compare the agricultural sustainability development levels of distinct regions at the same point in time. When comparing temporal differences, we normalize the data for the same region at different points in time. In this case, the new equations can be expressed in Equations (11) and (12), where Y represents the total duration of the data, y denotes the ordinal year with an interval of one year, and other variables remain unchanged. Unless otherwise specified, the following analysis will evaluate agricultural sustainability development levels for different regions from both spatial and temporal perspectives.
X y n = A y n min A 1 n , , A Y n max A 1 n , , A Y n min A 1 n , , A Y n
X y n = max A 1 n , , A Y n A y n max A 1 n , , A Y n min A 1 n , , A Y n

2.4.2. Calculation and Processing of Evaluation Value

After standardizing all the collected data and combining the indicator weights and normalized data for each region, the development levels of different dimensions of regional agricultural sustainability ( S t i ) can be obtained, as shown in Equation (14).
S t i = n = 1 M i w i · x t i , t = 1 , 2 , , T
where i represents the four dimensions within the agricultural sustainable development evaluation system, while M denotes the number of component indicators for each dimension. The level of regional integrated agricultural sustainability( S t ) is:
S t = i = 1 4 s t i ,
In order to improve spatiotemporal analysis, the results of agricultural sustainability evaluation are further processed in this paper, and the processing is as follows:
  • For the temporal analysis, to comprehensively capture the overall development and internal disparities in agricultural sustainability, we conducted calculations of the mean ( S y ¯ ) and coefficient of variation ( C V ) for the results (as shown in Equations (15) and (16)). The  S y ¯  represents the central tendency of the data, while the  C V  quantifies the extent of variability relative to the mean, thereby shedding light on the degree of heterogeneity among the studied objects [50,51].
    S y ¯ = 1 n i = 1 n S i , y
    C V = σ S y ¯ × 100
    where y denotes a specific year,  S i , y  represents the agricultural sustainability level of region i in a specific year, n signifies the number of observed regions, and  σ  denotes the standard deviation of the dataset.
    Additionally, this study decomposes the agricultural sustainability level into dimensions to better reflect the changes in each dimension and their contributions to the overall level of agricultural sustainability. By calculating the contribution margin( C R ), we can clearly distinguish the strengths and weaknesses of the dimensions that contribute to agricultural sustainability and make targeted improvements [52,53]; the formula is as follows:
    C R t i = S t i / S t , i = 1 , 2 , , 4 , t = 1 , 2 , , T
  • In the spatial analysis, this study employs the Natural Breaks Classification(NBC) method to further process the agricultural sustainability levels of various regions. This method is based on the distribution characteristics of the data, attempting to group similar data values into the same category while maximizing the differences between different categories [54]. In this study, we classified the level of agricultural sustainability using NBC method according to 5 classification standards in the ecological field: low level, relatively low level, medium level, relatively high level, and high level [55]. The Jenks optimization algorithm is employed to determine the optimal breaking points for the classification [56]. The detailed calculation principle is described below, with Equations (18) and (19) representing the optimization algorithm:
    S S D c , k = min 1 j n k + 1 S S D c 1 , j 1 + 1 n j + 1 i = j n k + 1 ( S i S ¯ j ) 2
    S ¯ j = 1 n j + 1 i = j n S i
    In the above equations,  S S D c , k  represents the sum of squared deviations between classes for c classes, wherein  c = 5 , and k observations; n denotes the number of observed regions;  S i  stands for each individual agricultural sustainability level; and  S ¯ j  is the mean value from j to n.

3. Result

3.1. Index Weights Results

In this paper, the index weights of the agricultural sustainability evaluation system were determined by the ANP method. Following the procedures and formulas delineated in Section 2.2, this study employs the SuperDecisions software to determine the weights of the various elements (Figure 3).
At the indicator level, the weight of an indicator reflects its relative importance compared to all the indicators within the evaluation system. For a given evaluation indicator, the higher the weight, the greater its impact on the final agricultural sustainability assessment value. As illustrated in Figure 3, it is evident that Forest coverage rate (e12, weight: 0.1537), Air quality index (e13, weight: 0.1132), and Per capita disposable income (e15, weight: 0.1305) are the three most crucial indicators. In contrast, the Agricultural electrification level (e4) and Rural healthcare level (e17) are the indicators with the most negligible impact on the agricultural sustainability assessment, with weights of only 0.0104 and 0.0117, respectively. Aside from indicators such as Agricultural mechanization level (e3) and Fertilizer use intensity (e10), which also have relatively low impact weights, most indicator weights lie within a range of 0.3 to 0.7.
At the dimension level, according to Figure 3, the weight of Environmental Quality (C3) is significantly higher than the other three dimensions—at 0.3999—accounting for 40% of the total. The Agricultural Economy dimension follows, with a weight of 0.2138. Resource Utilization and Rural Society exhibit similar weights, at 0.1923 and 0.1939, respectively. This suggests that Environmental quality is the primary evaluation dimension in the agricultural sustainability assessment process, with the other three dimensions having relatively less influence.

3.2. Temporal Analysis of Agricultural Sustainability

Upon obtaining the index weights, the data are processed using the time normalization method presented in Section 2.3 (Equations (11) and (12)), followed by calculating the agricultural sustainability levels of Chinese provinces from 2005 to 2020 using Equations (13) and (14). Further processing is conducted based on these results (Equations (15)–(17)) and shown in Section 3.3.1 and Section 3.3.2. In accordance with the paper’s formatting and presentation requirements, the specific results of the agricultural sustainability levels for Chinese provinces are not shown here. For detailed information, refer to Table A1 in the Appendix A.

3.2.1. Overall Temporal Variation Analysis

Figure 4 illustrates the overall temporal variation of agricultural sustainability by province in China. From 2005 to 2020, there is a general trend of increasing agricultural sustainability levels in all provinces.
As evident from Figure 4a, the agricultural sustainability levels of various Chinese provinces exhibit similar trends—displaying significant increases except for a substantial decline in 2013. It is important to note that the rate of increase in agricultural sustainability levels varies across provinces. Some provinces experienced relatively steady growth, while others exhibited more fluctuating patterns. For instance, Beijing exhibited a more consistent increase throughout the period, while provinces such as Liaoning and Jilin showed more variation in their trends. Provinces such as Anhui, Hubei, Hunan, and Sichuan demonstrate the most significant increases in agricultural sustainability, with growth rates exceeding 93%.
As indicated by Figure 4b, the average agricultural sustainability levels of Chinese provinces exhibit a rapid growth trend, rising from 0.2464 in 2005 to 0.8407 in 2020. This upward trajectory was interrupted by fluctuations in 2012 and 2013 when the average value declined to 0.416 before rebounding to 0.518. On the other hand, the coefficient of variation displays a more complex pattern. Initially, the values showed a pronounced decrease; however, this trend was disrupted by a notable increase in 2013, with the coefficient reaching its peak at 0.160. Following this apex, the coefficient gradually declined, stabilizing in 2019 and 2020. The trends in mean values and CVs indicate that while the agricultural sustainability capacities of Chinese provinces have been rapidly improving, there is a convergence effect of low-efficiency regions catching up with high-efficiency regions.

3.2.2. Dimensional Temporal Variation Analysis

Figure 5 illustrates the temporal changes in the mean values and contribution rates of different dimensions of agricultural sustainability across regions. Between 2005 and 2020, although the scores of each dimension have been continuously increasing, the contribution rates of these dimensions have experienced various changes.
As seen in Figure 5a, the sustainability levels of the four dimensions—agricultural economy, resource utilization, environmental quality, and rural society—all exhibit a continuous upward trend. Among these, the rural society dimension displays the most significant change, increasing more than 20-fold from 0.0077 in 2005 to 0.1856 in 2020. Environmental quality consistently remains the most critical component among the four indicators, with a score of 0.3571 in 2020, significantly higher than the other dimensions. Agricultural economy and resource utilization dimensions demonstrate similar, stable upward trends, rising from 0.0304 and 0.0541 in 2005 to 0.1603 and 0.1377 in 2020, respectively.
According to Figure 5b, the contribution rates of the four dimensions—Agriculture economy, Resource utilization, Environmental quality, and Rural society—exhibit different changes during the process but ultimately stabilize. Between 2005 and 2020, the proportion of the Agriculture economy dimension of the overall agricultural sustainability demonstrates a generally increasing trend, rising from 12% in 2005 to its peak of 25% in 2013, then stabilizing between 17% and 20%. Conversely, the Resource utilization dimension shows a declining trend, decreasing from its peak of 25% in 2008 to 16% in 2020. Throughout the entire period, the Environmental quality dimension consistently maintains the highest proportion among the four categories, despite experiencing a decline of more than 20%, ultimately remaining between 40% and 42%. Lastly, the contribution of the Rural society dimension to agricultural sustainability significantly increases, rising from 3% in 2005 to eventually stabilizing around 22%.

3.3. Spatial Analysis of Agricultural Sustainability

This section employs the same processing approach as Section 3.3, utilizing the spatial normalization method (Equations (9) and (10)) to process the data. Combined with the index weights from Section 3.2, each province’s relative agricultural sustainability scores are calculated using Equations (13) and (14) for each year. Finally, the NBC method is applied for classification and comparison. Similarly, detailed information can be found in Appendix A in Table A2.

3.3.1. Overall Spatial Variation Analysis

In order to compare more intuitively the spatial differences in the level of agricultural sustainability among different provinces, this study uses ArcGIS software to create maps of agricultural sustainability levels for Chinese provinces in 2005, 2010, 2015, and 2020, with a 5-year interval between each map. The resulting maps are presented in Figure 6.
As shown in Figure 6, the differences in agricultural sustainability among many regions remain relatively stable. Between 2005 and 2020, Zhejiang and Fujian’s relative agricultural sustainability levels have consistently been higher than those of other regions. Beijing, Jiangxi, and Guangdong maintain stable positions in the second tier. Sichuan and Tibet are among the regions with medium agricultural sustainability levels. Henan and Anhui exhibit relatively low agricultural sustainability, while Gansu and Qinghai rank at the bottom.
Apart from these regions, the agricultural sustainability levels of other areas have exhibited varying changes over the 16 years. For example, the relative agricultural sustainability ratings of provinces such as Hebei and Inner Mongolia have displayed different declining trends, with Jilin experiencing the most significant drop (from a high level in 2005 to a medium level in 2020). In contrast, the agricultural sustainability ratings of Heilongjiang, Hunan, and four other provinces have risen. With the exception of Jiangsu, which exhibited relatively large fluctuations in its agricultural sustainability (low level in 2005, medium level in 2010, high level in 2015, and medium level in 2020), nine areas, including Tianjin, Shanxi, and Shandong, have fluctuated between adjacent ratings.

3.3.2. Dimentional Spatial Variation Analysis

Figure 7 displays the spatial differences among various dimensions of agricultural sustainability in 2005, 2010, 2015, and 2020. This table provides an intuitive representation of the comparative advantages and disparities among different dimensions across regions at different time points.
As shown in Figure 7, during the period of 2005–2020, some regions demonstrated remarkable advantages in single dimensions. For example, Beijing’s sustainability level in the rural society dimension was significantly higher than other regions, Tibet led substantially in the resources utilization dimension, and Yunnan consistently held an advantage in the environmental quality dimension. In addition, some regions gradually developed their own advantageous dimensions over the 16-year period. For instance, Heilongjiang surpassed Jiangsu and Jilin in the agricultural economy dimension after 2010 and maintained a leading position. Fujian’s sustainability level in the resources utilization dimension gradually exceeded Xinjiang and ranked only behind Tibet in 2020.
Furthermore, Figure 7 also reveals the unbalanced development of agricultural sustainability across various regions in China. Generally, development in the environmental quality dimension is the most balanced among all provinces, with universally high sustainability levels and relatively small gaps. In the resources utilization dimension, by contrast, most regions exhibit the lowest sustainability scores among the four dimensions, with only a few exceptions such as Tibet, Xinjiang, and Fujian. A similar phenomenon is observed in the rural society dimension. Lastly, the agricultural economy dimension is the most distinctive, showing a divergent trend. It evolves from initially high scores across most regions to a gradual concentration in advantageous areas such as Heilongjiang, Jiangsu, Shandong, and Henan, with relative declines in other regions such as Zhejiang and Fujian.

4. Discussion

In this study, we employ the ANP method to determine index weights for an agricultural sustainability evaluation system and assess the spatiotemporal changes in agricultural sustainability levels across Chinese provinces from 2005 to 2020. Our findings are generally consistent with previous research but also show some differences. Based on these findings, we offer new insights into the constituent factors, development levels, and future policy recommendations for agricultural sustainability in China.
Considering the index weight results, our analysis emphasizes the importance of environmental quality in driving agricultural sustainability outcomes, which diverges from existing studies focusing more on the agricultural economic dimension [15,28]. The potential reasons for this difference may be attributed to our consideration of the interrelationships among different dimensions and indicators when calculating weights and variations in experts’ subjective judgments. Moreover, although the per capita disposable income in rural areas is significantly higher than other dimension indicators at the index level, the overall weight of the social dimension does not show a significant difference compared to the agricultural economy and resource utilization dimensions. This insight suggests that, while prioritizing environmental protection in promoting sustainable agricultural development, it is essential to pay attention to the coordinated development of the agricultural economy, resource utilization, and rural society [57].
In terms of temporal trends, our analysis first reconfirms the findings of Liu et al. (2019) [58] and Hu et al. (2022) [59] that the level of agricultural sustainability has significantly improved in all regions of China, indicating that China’s agricultural policies and practices are effective in promoting sustainable agriculture. However, the considerable decline in agricultural sustainability observed in various provinces in 2013 may suggest the impact of external factors, such as the peak severity of adverse weather conditions such as smog that year [60]. This implies that while maintaining agricultural sustainability growth, regions should focus on enhancing their response to external shocks and improving the stability of agricultural sustainability growth. Secondly, our study confirms the convergence effect of low-efficiency regions catching up with high-efficiency regions in terms of agricultural sustainability across China [61], highlighting the importance of targeted policies and interventions in addressing regional disparities. Lastly, the analysis of the temporal changes in each dimension reveals that the increasing importance of the rural social dimension reflects the significant role of rural development in enhancing agricultural sustainability [62]. Simultaneously, the persistent high contribution of the environmental quality dimension reiterates the necessity of prioritizing environmental protection in agricultural policies and practices [63].
Regarding spatial disparity, overall, some regions, such as Zhejiang and Fujian, consistently exhibit higher levels of agricultural sustainability, while others, such as Gansu and Qinghai, perform poorly, and some regions show considerable fluctuations in rankings. This indicates that agricultural sustainability levels in China still present specific imbalances and instabilities in development [64,65]. It is worth noting that regions such as Shandong and Henan, which are considered areas with high agricultural sustainability in some studies [65], are only ranked as moderate or even relatively low in our findings. This could be due to our evaluation process placing greater emphasis on the environmental quality dimension and using many per capita indicators in constructing the index system, which weakens these regions’ overall leading advantages in economic, resource, and social dimensions. At the dimension level, some regions demonstrate distinct advantages in specific dimensions, such as Beijing in the rural social dimension, Tibet in resource utilization, and Yunnan in environmental quality. The dimensional advantages of regions such as Heilongjiang and Fujian have gradually been established. These findings emphasize the importance of leveraging regional strengths and addressing weaknesses to promote balanced and sustainable agricultural development [66]. In policy formulation and agricultural practices, it is crucial to recognize these differences and strive for balanced development across all aspects. By capitalizing on the strengths of advantageous dimensions and addressing the shortcomings of weaker ones, we can work towards achieving sustainable, resilient, and inclusive agricultural systems that contribute to the well-being of both humanity and our planet.

5. Conclusions

In conclusion, our study contributes to a broad understanding of agricultural sustainability assessment and management by providing a comprehensive and detailed evaluation system based on the ANP methodology, which examines the spatial and temporal variation in agricultural sustainability across regions in China over the last 16 years. Our findings provide valuable insights and policy recommendations for improving agricultural sustainability in China, and they can serve as a valuable reference for other countries seeking to advance sustainable agricultural development.
The main contributions of this study can be summarized as follows. First, by considering the various dimensions and indicators of agricultural sustainability as an interacting whole, our analysis provides a more comprehensive and accurate assessment of China’s agricultural sustainability. Second, our examination of the different dimensions reveals the necessity for adopting a sustainable development model centered on environmental protection, with coordinated growth in agriculture, resource utilization, and social aspects after reaching a specific scale of economic development. Third, our analysis of the spatiotemporal differences across China’s regions offers a theoretical basis for targeted policymaking and interventions in specific areas.
However, due to data limitations, this study omits the evaluation of a few regions and cannot incorporate all relevant evaluation indicators. Furthermore, the assessment method employed in this research is based on subjective and relative perspectives and does not delve deeply into the influencing factors. Therefore, objectively measuring agricultural sustainability and investigating the external impact factors, internal driving forces, and the effectiveness of policies and measures affecting agricultural sustainability are directions for our future research.

Author Contributions

Conceptualization, Y.-W.D. and Y.-P.F.; writing—original draft, Y.-P.F.; writing—review and editing, Y.-W.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major Program of National Social Science Foundation of China grant number No. 18ZDA055.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: http://www.stats.gov.cn/sj/ (accessed on 30 March 2023).

Acknowledgments

We would like to express our gratitude to Wang Y.-C. and Cao Z.-X. for their invaluable assistance in the presentation of this article. Their expertise and dedication have been crucial to the success of this research. We deeply appreciate their generosity and support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Tables of Results

Table A1. Temporal evolution of agricultural sustainability levels (2005–2020).
Table A1. Temporal evolution of agricultural sustainability levels (2005–2020).
Region2005200620072008200920102011201220132014201520162017201820192020
Beijing0.35600.36140.37970.51610.53310.53610.57100.59370.41080.47560.52690.55060.56950.60070.56420.6677
Tianjin0.23480.26780.28950.28380.34700.35280.41190.44300.32220.47410.53940.58360.63060.66040.70780.7903
Hebei0.22840.25180.31520.34020.38840.41040.49000.53760.33250.46880.55530.59190.60150.69200.73270.8286
Shanxi0.19130.23120.26040.28580.37330.41180.49750.50590.47170.57340.58390.66350.65580.70130.72740.8363
Neimeng0.25480.25820.28670.34390.41640.45190.50820.52670.44310.56350.60890.61430.59390.70140.77240.7950
Liaoning0.36220.34310.33040.35660.43870.51090.59640.67920.46330.45010.53720.62300.64700.68930.77930.8088
Jilin0.34430.30480.30220.32030.40730.49090.50250.53120.36730.43550.46890.58500.60140.74950.71870.7847
Heilongjiang0.19780.20470.18050.25290.37680.39750.43710.47740.40360.50820.51820.62170.63320.82790.84190.8821
Jiangsu0.19670.22380.27330.29770.48270.49150.55550.60660.41630.45610.53230.59850.62840.73640.69570.7845
Xinjiang0.23300.22170.26840.25140.34860.41080.41780.45720.43670.59250.66980.64300.63980.67850.73680.7263
Zhejiang0.21810.21050.22520.25940.46060.50150.49260.56790.47710.56190.62470.65220.66120.78100.76460.7906
Anhui0.23910.23880.26490.30080.36190.47380.45690.52080.42370.48320.63790.71790.70890.76740.80620.9397
Fujian0.25380.21530.25100.26060.39770.44020.44930.50680.52130.51560.59900.66870.65030.74680.78220.8068
Jiangxi0.23290.22940.23260.26330.35240.43660.41750.48110.39190.56660.61970.65350.65410.78980.80020.8255
Shandong0.22430.23140.28310.30890.30760.38010.44950.47310.29650.50210.55760.61330.66570.72220.75240.8867
Henan0.25140.24410.29640.31480.40390.43400.47540.51510.33570.50430.55660.57920.63340.67540.71630.8645
Hubei0.21420.19550.23300.28320.43270.46590.44260.57910.51010.53110.62060.69410.73810.75830.80310.9395
Hunan0.16280.18510.20980.27240.37910.37430.41720.47230.42840.57660.65170.69990.72780.75010.82490.9548
Guangdong0.22260.22360.19310.31450.39900.44040.44560.49300.36340.48030.52760.57900.56430.68110.69570.7623
Guangxi0.22470.25150.23030.31260.37600.39690.40050.46450.39670.54100.61370.66050.68770.81330.80240.8600
Hainan0.25980.31890.28000.44120.52340.51190.54800.57320.49760.51050.52950.64050.61420.74050.68280.7770
Chongqing0.18360.17990.27310.31350.35860.40220.43160.47020.34630.54590.58400.65040.67770.80580.77160.8597
Sichuan0.28630.19910.26520.28380.34090.39430.44090.49630.43240.55000.60550.66210.71680.79550.84540.9329
Guizhou0.15050.13890.18540.23770.32850.39080.34340.42990.31560.55670.64730.68150.71080.82630.81630.8760
Yunnan0.24500.23650.28600.28060.30610.35650.41840.38530.35030.57210.60230.69440.73920.79040.79750.9180
Tibet0.40380.33020.36960.42910.45320.48990.50550.50150.53620.51820.49970.55380.69210.79360.85940.8062
Shaanxi0.32390.31820.30700.30440.48410.51950.58270.59660.40520.44540.55660.53620.57450.66740.77170.8296
Gansu0.23190.11060.20810.24280.39660.43580.50710.58310.53070.56760.60370.60210.64740.70590.83120.9244
Qinghai0.27250.23420.18680.33070.51510.52630.52820.50360.39160.51550.53620.55110.64270.71380.82400.8998
Ningxia0.19250.22900.22880.25690.40030.44540.47790.57190.44760.56920.57640.58750.59870.71190.80800.8643
Xinjiang0.23300.22170.26840.25140.34860.41080.41780.45720.43670.59250.66980.64300.63980.67850.73680.7263
Table A2. Spatial differences of agricultural sustainability levels (2005–2020).
Table A2. Spatial differences of agricultural sustainability levels (2005–2020).
Region2005200620072008200920102011201220132014201520162017201820192020
Beijing0.43470.46600.43090.49250.49900.49730.48260.44750.45160.44190.44290.42920.44090.42190.41710.4422
Tianjin0.40930.42470.40550.37930.38560.39160.40730.37450.37210.38290.39400.38710.38300.32980.36760.3755
Hebei0.38670.41020.39080.41090.41750.42190.41960.40640.33760.34590.36850.34460.33630.28270.33020.3303
Shanxi0.25890.29790.27570.29320.29780.31550.31340.31990.31150.31960.31920.30430.27320.24410.27120.2658
Neimeng0.44340.45950.46440.47870.47880.48690.49110.49100.46090.46010.46670.45500.43430.42170.45640.4352
Liaoning0.46470.47290.46260.46100.46250.47060.48670.47450.44270.40860.41300.41860.41890.38200.41910.3989
Jilin0.52190.51830.48940.50270.48160.49060.50810.48360.46940.45950.44720.46580.45470.44130.44660.4310
Heilongjiang0.47120.49400.48460.50310.50450.52240.53430.53290.53380.53660.51820.53280.52510.58170.53860.5169
Jiangsu0.35830.38750.37900.38690.42910.42100.43520.42250.42810.42410.43070.42290.43550.36820.41750.4351
Zhejiang0.53650.54870.52960.53480.54270.53570.54920.54190.55030.54490.54740.54290.54930.53400.54750.5741
Anhui0.32400.34350.31640.34340.34340.36640.35430.36940.36160.32360.36300.36000.34830.28910.34500.3680
Fujian0.55240.54650.54480.53410.53880.53860.54340.55080.57580.56000.56850.57240.56190.52950.56540.5703
Jiangxi0.46970.47970.48490.49230.49000.48420.49230.46240.45940.48980.49040.47700.46850.41360.45950.4644
Shandong0.35240.41210.39730.38620.39100.40250.41020.39750.33960.34300.33460.33480.34480.27200.32120.3376
Henan0.34380.37230.37150.38090.39950.38690.38150.36830.34250.34280.33120.31980.32590.24280.31110.3299
Hubei0.36140.37120.35640.38160.41390.40750.39770.43550.43190.40630.40910.41580.42670.36950.40550.4266
Hunan0.36140.39780.40390.43690.45800.45180.46050.43850.43600.45280.46020.45170.45320.38440.44170.4607
Guangdong0.44740.46730.43240.45700.46630.47290.47840.47420.44130.45630.45930.44930.43790.40320.43120.4528
Guangxi0.40730.40630.38280.39080.40880.39270.39840.39500.39630.40520.41590.41610.41200.40750.41290.4194
Hainan0.41430.43260.41340.44200.44250.43830.44810.45120.45240.43940.43850.44070.43420.42090.43360.4289
Chongqing0.30400.28280.27350.28500.29580.31080.32040.32980.30170.32440.34250.33740.33460.31920.34910.3653
Sichuan0.40350.39780.41200.40310.42010.42070.42980.38440.40970.39920.39210.38300.39660.36010.41420.4124
Guizhou0.34040.33020.33430.33770.31930.36980.35320.37070.32940.36780.39070.39110.39050.40110.39500.3973
Yunnan0.41990.41820.42610.41390.41480.41710.44400.43210.44660.47200.46960.47210.47430.45930.47600.4827
Tibet0.39930.39150.38620.39680.40150.40420.41420.41050.41780.38240.38100.37490.40700.40010.43880.4135
Shaanxi0.36680.37650.36270.36980.38140.39290.39430.38090.36530.35160.37990.33670.32860.31180.35860.3586
Gansu0.22980.21690.23950.21130.21530.21550.22050.22950.28820.29100.28980.26800.26780.24930.29920.2991
Qinghai0.23760.23220.20630.21200.22270.25820.25190.23560.24360.25110.26090.23570.26390.25670.29590.2835
Ningxia0.31080.31520.30160.30610.31440.32520.32140.31330.31980.31660.30990.29210.28200.27200.32390.3057
Xinjiang0.29090.30530.28180.27570.29510.31640.31040.31160.35570.34810.34830.33240.32350.31230.33470.3108

References

  1. Piñeiro, V.; Arias, J.; Dürr, J.; Elverdin, P.; Ibáñez, A.M.; Kinengyere, A.; Opazo, C.M.; Owoo, N.; Page, J.R.; Prager, S.D.; et al. A scoping review on incentives for adoption of sustainable agricultural practices and their outcomes. Nat. Sustain. 2020, 3, 809–820. [Google Scholar] [CrossRef]
  2. MacLaren, C.; Mead, A.; van Balen, D.; Claessens, L.; Etana, A.; de Haan, J.; Haagsma, W.; Jäck, O.; Keller, T.; Labuschagne, J.; et al. Long-term evidence for ecological intensification as a pathway to sustainable agriculture. Nat. Sustain. 2022, 5, 770–779. [Google Scholar] [CrossRef]
  3. Viana, C.M.; Freire, D.; Abrantes, P.; Rocha, J.; Pereira, P. Agricultural land systems importance for supporting food security and sustainable development goals: A systematic review. Sci. Total. Environ. 2022, 806, 150718. [Google Scholar] [CrossRef]
  4. Legg, W. Sustainable Agricultural Development for Food Security and Nutrition: What Roles for Livestock? Technical Report No. 1916-2017-1398; Food and Agriculture Organization of the United Nations and the High Level Panel of Experts on Food Security and Nutrition: Washington, DC, USA, 2017. [Google Scholar]
  5. FAO. The State of Food and Agriculture 2020: Overcoming Water Challenges in Agriculture; Number 2020 in The State of Food and Agriculture (SOFA); FAO: Rome, Italy, 2020. [Google Scholar] [CrossRef]
  6. Streimikis, J.; Baležentis, T. Agricultural sustainability assessment framework integrating sustainable development goals and interlinked priorities of environmental, climate and agriculture policies. Sustain. Dev. 2020, 28, 1702–1712. [Google Scholar] [CrossRef]
  7. Jing, J.; Cong, W.F.; Bezemer, T.M. Legacies at work: Plant-soil-microbiome interactions underpinning agricultural sustainability. Trends Plant Sci. 2022, 27, 781–792. [Google Scholar] [CrossRef]
  8. Sachs, J.; Kroll, C.; Lafortune, G.; Fuller, G.; Woelm, F. Sustainable Development Report 2022; Cambridge University Press: Cambridge, UK, 2022. [Google Scholar]
  9. Binder, C.R.; Feola, G.; Steinberger, J.K. Considering the normative, systemic and procedural dimensions in indicator-based sustainability assessments in agriculture. Environ. Impact Assess. Rev. 2010, 30, 71–81. [Google Scholar] [CrossRef]
  10. Marchand, F.; Debruyne, L.; Triste, L.; Gerrard, C.; Padel, S.; Lauwers, L. Key characteristics for tool choice in indicator-based sustainability assessment at farm level. Ecol. Soc. 2014, 19, 46. [Google Scholar] [CrossRef] [Green Version]
  11. de Olde, E.M.; Oudshoorn, F.W.; Sorensen, C.A.G.; Bokkers, E.A.M.; de Boer, I.J.M. Assessing sustainability at farm-level: Lessons learned from a comparison of tools in practice. Ecol. Indic. 2016, 66, 391–404. [Google Scholar] [CrossRef]
  12. Riera, A.; Duluins, O.; Schuster, M.; Baret, P.V. Accounting for diversity while assessing sustainability: Insights from the Walloon bovine sectors. Agron. Sustain. Dev. 2023, 43, 30. [Google Scholar] [CrossRef]
  13. von Wirén-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. Lv, Y.; Zhang, C.; Ma, J.; Yun, W.; Gao, L.; Li, P. Sustainability Assessment of Smallholder Farmland Systems: Healthy Farmland System Assessment Framework. Sustainability 2019, 11, 4525. [Google Scholar] [CrossRef] [Green Version]
  15. Abdar, Z.K.; Amirtaimoori, S.; Mehrjerdi, M.R.Z.; Boshrabadi, H.M. A composite index for assessment of agricultural sustainability: The case of Iran. Environ. Sci. Pollut. Res. 2022, 29, 47337–47349. [Google Scholar] [CrossRef]
  16. Peng, W.; Zheng, H.; Robinson, B.E.; Li, C.; Li, R. Comparing the importance of farming resource endowments and agricultural livelihood diversification for agricultural sustainability from the perspective of the food–energy–water nexus. J. Clean. Prod. 2022, 380, 135193. [Google Scholar] [CrossRef]
  17. Van Cauwenbergh, N.; Biala, K.; Bielders, C.; Brouckaert, V.; Franchois, L.; Garcia Cidad, V.; 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]
  18. Meul, M.; Van Passel, S.; Nevens, F.; Dessein, J.; Rogge, E.; Mulier, A.; Van Hauwermeiren, A. MOTIFS: A monitoring tool for integrated farm sustainability. Agron. Sustain. Dev. 2008, 28, 321–332. [Google Scholar] [CrossRef] [Green Version]
  19. Scialabba, N. SAFA Guidelines: Sustainability Assessment of Food and Agriculture Systems. Version 3.0; Food and Agriculture Organization of the United Nations: Wanshington, DC, USA, 2013. [Google Scholar]
  20. Marta-Costa, A.A.; da Silva, E.L.S. Methods and Procedures for Building Sustainable Farming Systems: Application in the European Context; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  21. Talukder, B.; Blay-Palmer, A.; Hipel, K.W.; vanLoon, G.W. Elimination Method of Multi-Criteria Decision Analysis (MCDA): A Simple Methodological Approach for Assessing Agricultural Sustainability. Sustainability 2017, 9, 287. [Google Scholar] [CrossRef] [Green Version]
  22. Du, Y.W.; Zhong, J.J. Dynamic multicriteria group decision-making method with automatic reliability and weight calculation. Inf. Sci. 2023, 634, 400–422. [Google Scholar] [CrossRef]
  23. Delnevo, G.; Girau, R.; Ceccarini, C.; Prandi, C. A Deep Learning and Social IoT Approach for Plants Disease Prediction Toward a Sustainable Agriculture. IEEE Internet Things J. 2022, 9, 7243–7250. [Google Scholar] [CrossRef]
  24. Talukder, B.; Hipel, K.W.; vanLoon, G.W. Using multi-criteria decision analysis for assessing sustainability of agricultural systems. Sustain. Dev. 2018, 26, 781–799. [Google Scholar] [CrossRef]
  25. Scown, M.W.; Winkler, K.J.; Nicholas, K.A. Aligning research with policy and practice for sustainable agricultural land systems in Europe. Proc. Natl. Acad. Sci. USA 2019, 116, 4911–4916. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, S.; Yang, Y.; Zheng, H.; Mi, C.; Ma, T.; Shi, R. A framework for assessing sustainable agriculture and rural development: A case study of the Beijing-Tianjin-Hebei region, China. Environ. Impact Assess. Rev. 2022, 97, 106861. [Google Scholar] [CrossRef]
  27. Du, Y.W.; Shen, X.L. Group hierarchical DEMATEL method for reaching consensus. Comput. Ind. Eng. 2023, 175, 108842. [Google Scholar] [CrossRef]
  28. Qin, Y.; He, J.; Wei, M.; Du, X. Challenges Threatening Agricultural Sustainability in Central Asia: Status and Prospect. Int. J. Environ. Res. Public Health 2022, 19, 6200. [Google Scholar] [CrossRef] [PubMed]
  29. Dong, F.; Mitchell, P.D.; Colquhoun, J. Measuring farm sustainability using data envelope analysis with principal components: The case of Wisconsin cranberry. J. Environ. Manag. 2015, 147, 175–183. [Google Scholar] [CrossRef]
  30. Du, Y.W.; Li, X.X. Hierarchical DEMATEL method for complex systems. Expert Syst. Appl. 2021, 167, 113871. [Google Scholar] [CrossRef]
  31. Tzouramani, I.; Mantziaris, S.; Karanikolas, P. Assessing Sustainability Performance at the Farm Level: Examples from Greek Agricultural Systems. Sustainability 2020, 12, 2929. [Google Scholar] [CrossRef] [Green Version]
  32. Saaty, T.L. Decision Making with Dependence and Feedback: The Analytic Network Process: The Organization and Prioritization of Complexity; RWS Publications: Pittsburgh, PA, USA, 1996. [Google Scholar]
  33. Jharkharia, S.; Shankar, R. Selection of logistics service provider: An analytic network process (ANP) approach. Omega 2007, 35, 274–289. [Google Scholar] [CrossRef]
  34. Saaty, T.L. Theory and Applications of the Analytic Network Process: Decision Making with Benefits, Opportunities, Costs, and Risks; RWS Publications: Pittsburgh, PA, USA, 2005. [Google Scholar]
  35. Saaty, T.L. The Modern Science of Multicriteria Decision Making and Its Practical Applications: The AHP/ANP Approach. Oper. Res. 2013, 61, 1101–1118. [Google Scholar] [CrossRef]
  36. Erol, H.; Dikmen, I.; Atasoy, G.; Birgonul, M.T. An analytic network process model for risk quantification of mega construction projects. Expert Syst. Appl. 2022, 191, 116215. [Google Scholar] [CrossRef]
  37. Lin, W.C. Maritime Environment Assessment and Management Using through Balanced Scorecard by Using DEMATEL and ANP Technique. Int. J. Environ. Res. Public Health 2022, 19, 2873. [Google Scholar] [CrossRef]
  38. Kheybari, S.; Rezaie, F.M.; Farazmand, H. Analytic network process: An overview of applications. Appl. Math. Comput. 2020, 367, 124780. [Google Scholar] [CrossRef]
  39. Ghose, B. Food security and food self-sufficiency in China: From past to 2050. Food Energy Secur. 2014, 3, 86–95. [Google Scholar] [CrossRef]
  40. Garnett, T.; Appleby, M.C.; Balmford, A.; Bateman, I.J.; Benton, T.G.; Bloomer, P.; Burlingame, B.; Dawkins, M.; Dolan, L.; Fraser, D.; et al. Sustainable Intensification in Agriculture: Premises and Policies. Science 2013, 341, 33–34. [Google Scholar] [CrossRef] [PubMed]
  41. Altieri, M.A. Agroecology: The Science of Sustainable Agriculture, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar] [CrossRef]
  42. Janker, J.; Mann, S.; Rist, S. What is Sustainable Agriculture? Critical Analysis of the International Political Discourse. Sustainability 2018, 10, 4707. [Google Scholar] [CrossRef] [Green Version]
  43. Guo, H.; Gu, F.; Peng, Y.; Deng, X.; Guo, L. Does Digital Inclusive Finance Effectively Promote Agricultural Green Development?—A Case Study of China. Int. J. Environ. Res. Public Health 2022, 19, 6982. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, Z.; Huang, L.; Yin, L.; Wang, Z.; Zheng, D. Evaluation of Sustainable and Analysis of Influencing Factors for Agriculture Sector: Evidence from Jiangsu Province, China. Front. Environ. Sci. 2022, 10, 836002. [Google Scholar] [CrossRef]
  45. Vysochyna, A.; Stoyanets, N.; Mentel, G.; Olejarz, T. Environmental Determinants of a Country’s Food Security in Short-Term and Long-Term Perspectives. Sustainability 2020, 12, 4090. [Google Scholar] [CrossRef]
  46. Nowak, A.; Krukowski, A.; Różańska-Boczula, M. Assessment of Sustainability in Agriculture of the European Union Countries. Agronomy 2019, 9, 890. [Google Scholar] [CrossRef] [Green Version]
  47. Lynch, J.; Donnellan, T.; Finn, J.A.; Dillon, E.; Ryan, M. Potential development of Irish agricultural sustainability indicators for current and future policy evaluation needs. J. Environ. Manag. 2019, 230, 434–445. [Google Scholar] [CrossRef]
  48. Saaty, T.L. Fundamentals of the analytic network process—Dependence and feedback in decision-making with a single network. J. Syst. Sci. Syst. Eng. 2004, 13, 129–157. [Google Scholar] [CrossRef]
  49. Li, Q.; Fan, C. Evaluation of Hydraulic-Tunnel-Lining Durability Based on ANP and Cloud-Model-Improved Matter-Element Theory. Sustainability 2022, 14, 11801. [Google Scholar] [CrossRef]
  50. Asfaw, A.; Simane, B.; Hassen, A.; Bantider, A. Variability and time series trend analysis of rainfall and temperature in northcentral Ethiopia: A case study in Woleka sub-basin. Weather. Clim. Extrem. 2018, 19, 29–41. [Google Scholar] [CrossRef]
  51. Jalilibal, Z.; Amiri, A.; Castagliola, P.; Khoo, M.B.C. Monitoring the coefficient of variation: A literature review. Comput. Ind. Eng. 2021, 161, 107600. [Google Scholar] [CrossRef]
  52. Du, Y.W.; Song, Q. Ecological Security Evaluation for Marine Ranching Based on the PLTS/ANP Method: A Case Study of Rongcheng. Complexity 2022, 2022, 2134296. [Google Scholar] [CrossRef]
  53. Tong, W.; Qin, G.; Xu, Y.; Fu, C. Sustainable growth, input factors, and technological progress in agriculture: Evidence from 1990 to 2020 in China. Front. Environ. Sci. 2023, 10, 1040356. [Google Scholar] [CrossRef]
  54. Cai, J.; Huang, B.; Song, Y. Using multi-source geospatial big data to identify the structure of polycentric cities. Remote. Sens. Environ. 2017, 202, 210–221. [Google Scholar] [CrossRef]
  55. Hou, K.; Li, X.; Wang, J.j.; Zhang, J. An analysis of the impact on land use and ecological vulnerability of the policy of returning farmland to forest in Yan’an, China. Environ. Sci. Pollut. Res. 2016, 23, 4670–4680. [Google Scholar] [CrossRef]
  56. Aditian, A.; Kubota, T.; Shinohara, Y. Comparison of GIS-based landslide susceptibility models using frequency ratio, logistic regression, and artificial neural network in a tertiary region of Ambon, Indonesia. Geomorphology 2018, 318, 101–111. [Google Scholar] [CrossRef]
  57. Liu, F.; Wang, C.; Luo, M.; Zhou, S.; Liu, C. An investigation of the coupling coordination of a regional agricultural economics-ecology-society composite based on a data-driven approach. Ecol. Indic. 2022, 143, 109363. [Google Scholar] [CrossRef]
  58. 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. [Google Scholar] [CrossRef] [Green Version]
  59. Hu, J.; Lyu, J.; Zhang, X. Evaluating Agricultural Sustainability and Green GDP in China: An Emergy Analysis. Int. J. Environ. Res. Public Health 2022, 19, 16735. [Google Scholar] [CrossRef] [PubMed]
  60. Chen, S.; Guo, J.; Song, L.; Cohen, J.B.; Wang, Y. Intra-seasonal differences in the atmospheric systems contributing to interannual variations of autumn haze pollution in the North China Plain. Theor. Appl. Climatol. 2020, 141, 389–403. [Google Scholar] [CrossRef]
  61. Zhang, Y.; Long, H.; Chen, S.; Ma, L.; Gan, M. The development of multifunctional agriculture in farming regions of China: Convergence or divergence? Land Use Policy 2023, 127, 106576. [Google Scholar] [CrossRef]
  62. Zhu, C.; Lin, Y.; Zhang, J.; Gan, M.; Xu, H.; Li, W.; Yuan, S.; Wang, K. Exploring the relationship between rural transition and agricultural eco-environment using a coupling analysis: A case study of Zhejiang Province, China. Ecol. Indic. 2021, 127, 107733. [Google Scholar] [CrossRef]
  63. Fusco, G.; Campobasso, F.; Laureti, L.; Frittelli, M.; Valente, D.; Petrosillo, I. The environmental impact of agriculture: An instrument to support public policy. Ecol. Indic. 2023, 147, 109961. [Google Scholar] [CrossRef]
  64. Yang, T.; Sun, Y.; Li, X.; Li, Q. An ecosystem elasticity perspective of paddy ecosystem sustainability evaluation: The case of China. J. Clean. Prod. 2021, 295, 126292. [Google Scholar] [CrossRef]
  65. Liang, X.; Gong, Q.; Li, S.; Huang, S.; Guo, G. Regional agricultural sustainability assessment in China based on a developed model. Environ. Dev. Sustain. 2022. [Google Scholar] [CrossRef]
  66. Bao, B.; Jiang, A.; Jin, S. What Drives the Fluctuations of “Green” Development in China’s Agricultural Sector? An Entropy Method Approach. Pol. J. Environ. Stud. 2022, 31, 3491–3507. [Google Scholar] [CrossRef]
Figure 1. The spatial distribution features of China’s three regions.
Figure 1. The spatial distribution features of China’s three regions.
Sustainability 15 09066 g001
Figure 2. Network model of agricultural sustainability.
Figure 2. Network model of agricultural sustainability.
Sustainability 15 09066 g002
Figure 3. The index weights of agricultural sustainability evaluation system.
Figure 3. The index weights of agricultural sustainability evaluation system.
Sustainability 15 09066 g003
Figure 4. Temporal variations in regional agricultural sustainability.
Figure 4. Temporal variations in regional agricultural sustainability.
Sustainability 15 09066 g004
Figure 5. Dimensional temporal variation in agricultural sustainability.
Figure 5. Dimensional temporal variation in agricultural sustainability.
Sustainability 15 09066 g005
Figure 6. Spatial variations in regional agricultural sustainability.
Figure 6. Spatial variations in regional agricultural sustainability.
Sustainability 15 09066 g006
Figure 7. Dimensional spatial variation in agricultural sustainability.
Figure 7. Dimensional spatial variation in agricultural sustainability.
Sustainability 15 09066 g007
Table 1. Agricultural sustainability evaluation index system.
Table 1. Agricultural sustainability evaluation index system.
DimensionsIndexesDescriptionUnitsDirectionReferences
Agricultural
Economy
e1: Per labor agricultural gross valueAgricultural, forestry, animal husbandry, and fishery output value/rural populationYuan/rural labor+[15,44]
Agricultural
Economy
e2: Per labor grain outputTotal grain output/rural laborkg/person+[15]
e3: Agricultural mechanization levelTotal agricultural machinery power/cultivated land areakW/hm 2 +[15]
e4: Agricultural electrification levelRural electricity consumption/rural populationkW·h/person+[43,45]
e5: Land productivityTotal grain output/cultivated land areakg/hm 2 +[18]
Resource
Utilization
e6: Per capita arable land areaTotal arable land area at year-end/total population at year-endhm 2 /person+[46]
e7: Per capita water resource availabilityTotal water resources availability/total populationm 3 /person+[44]
e8: Rural power supply levelTotal rural power supply generation/rural population%+[26]
e9: Effective irrigation rateEffective irrigation area/sown area%+[12,26]
Environmental
Quality
e10: Fertilizer use intensityFertilizer consumption/sown areakg/hm 2 -[15,43]
e11: Pesticide use intensityPesticide consumption/sown areakg/hm 2 -[12,31]
e12: Forest coverage rateForest area/land area%+[12,15]
e13: Air quality indexDays of good air quality/365/+[18]
e14: Soil and water conservation areaSoil and water conservation areakhm 2 +[6]
Rural Societye15: Per capita disposable incomePer capita disposable income of rural residentsYuan/person+[18,26]
e16: Engel’s coefficient of rural residentsFood expenditure cost/total rural resident consumption expenditure/-[44]
e17: Rural healthcare levelAverage health personnel per thousand rural population/+[15]
e18: Education level of rural residentsAverage years of educationyears+[24,26]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Du, Y.-W.; Fan, Y.-P. Spatiotemporal Dynamics of Agricultural Sustainability Assessment: A Study across 30 Chinese Provinces. Sustainability 2023, 15, 9066. https://doi.org/10.3390/su15119066

AMA Style

Du Y-W, Fan Y-P. Spatiotemporal Dynamics of Agricultural Sustainability Assessment: A Study across 30 Chinese Provinces. Sustainability. 2023; 15(11):9066. https://doi.org/10.3390/su15119066

Chicago/Turabian Style

Du, Yuan-Wei, and Yi-Pin Fan. 2023. "Spatiotemporal Dynamics of Agricultural Sustainability Assessment: A Study across 30 Chinese Provinces" Sustainability 15, no. 11: 9066. https://doi.org/10.3390/su15119066

APA Style

Du, Y. -W., & Fan, Y. -P. (2023). Spatiotemporal Dynamics of Agricultural Sustainability Assessment: A Study across 30 Chinese Provinces. Sustainability, 15(11), 9066. https://doi.org/10.3390/su15119066

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