**1. Introduction**

Agriculture plays a central role in sustainable development, which determines the relationship between the global economy, society, and the natural world [1]. Since China implemented the reform and opening policy in 1978, China's agriculture has developed rapidly, producing considerable economic benefits and meeting the growing food demand of people [2,3]. China has done remarkable work in feeding 22% of the world's population with only 9% of the world's arable land, making an important contribution to international food security [4]. At the same time, China has also paid a huge environmental price. Environmental problems such as water eutrophication, soil acidification, air pollution, the reduction of biodiversity, etc., have become increasingly prominent, which greatly hinders the green and sustainable development of agriculture [5–7]. Therefore, China needs to coordinate the relationship between agricultural development and environmental protection urgently and promote the green transformation of agriculture.

In recent years, China has realized the importance of agricultural green development. Different ways have been adopted in different provinces in the process of agricultural green development, which have achieved varying degrees of effectiveness. However, how to realize the sustainable utilization of resources and sustainable development of agriculture while protecting the environment is still an urgent problem to be solved.

Agriculture Green Development is the future development direction of agriculture for China and the world [8–10]. In order to achieve green development in agriculture, it is necessary to make a quantitative evaluation of the current agricultural development level. In previous studies, in order to evaluate the sustainable agricultural development, many scholars have constructed different evaluation index systems for sustainable agricultural development around three dimensions, including society, economy, and ecology, evaluating the sustainable development of agriculture from the perspectives of farms, regions, and countries [11–16]. In recent years, under the background of green transformation of agriculture, scholars have studied from different perspectives on the construction of evaluation index system and evaluation methods of agricultural green development. For instance, Wei et al. selected 14 indicators from four dimensions of resource conservation, ecological conservation, environmental friendliness, and quality efficiency to construct China's agricultural green development index, then evaluating the level of agricultural green development in China [17]. Gong et al. selected 10 indicators from three dimensions of low-carbon economy, economic growth, and safe supply to construct a development index to evaluate the level of agricultural green development [18]. More scholars have used principal component analysis [19], entropy weight method [20–22], or analytic hierarchy process [23] to evaluate the level of agricultural green development in different provinces of China.

The overall objective of Agriculture Green Development is to coordinate "green" with "development" to realize the transformation of current agriculture with high resource consumption and high environmental costs into a green agriculture and countryside with high productivity, high resource use efficiency, and low environmental impact [10]. However, the resource and environment cost of agricultural production in China is too high. Irrational fertilization leads to serious acidification of croplands, which poses a threat to agricultural production and hinders China' s green development of agriculture seriously [5,24,25]. In areas with scarce water resources, the quality and safety of agricultural products cannot be guaranteed, which will also hinder the green development of agriculture [26]. Tu et al. pointed out that the green development of agriculture is not only reflected in the low agricultural non-point source pollution, but also in the high efficiency of agricultural production [27]. Inconsistent agricultural development conditions, different ecological environments, and different resource endowments in various regions have an impact on agricultural production and the level of agricultural development [28,29]. In addition, changes in institutions and policies, the level of urbanization, and the level of science and technology also have an impact on agricultural production efficiency [30]. Among many influencing factors, technological progress has been regarded as the main driving factor of agricultural production [31]. Meanwhile, the green development of agriculture puts forward higher requirements for technological progress, which requires the integration of multiple disciplines and technological innovation to realize the green development of agriculture [32,33].

Through reviewing literatures, since the concept of agricultural green development is relatively new, most studies mainly focus on two aspects: The first is the construction and level measurement of the evaluation index system of agricultural green development, and the second is the driving or restrictive factors that affect agricultural green development. The research on the spatial characteristics of agricultural green development and its influencing factors is still blank. Therefore, this article first evaluates the status of agricultural green development level in China through improved entropy weight method to measure the level of agricultural green development. Then, in order to study the spatial correlation and heterogeneity of agricultural green development in China, the article conducted an empirical analysis on the influencing factors of agricultural green development with panel data from 31 provinces in China. This study provides data support and empirical basis for the proposal and formulation of agricultural green development policy in different regions of China.

### **2. Materials and Methods**

### *2.1. Materials*

#### 2.1.1. Construction of Evaluation Index System for Agricultural Green Development Level

According to the definition of the level of agricultural green development and the research of other scholars, this article selects 16 indicators from four perspectives, including social development, economic benefits, resource input, and ecological environment to construct an evaluation index system for the level of agricultural green development in China, as shown in Table 1. All data are derived from the *China Statistical Yearbook*, *China Rural Statistical Yearbook*, *China Environment Statistical Yearbook*, and *Provincial Statistical Yearbook* from 2008 to 2019.



### 2.1.2. Influencing Factors of Agricultural Green Development

This paper used the level of agricultural green development in each province as the explained variable, and *urb*, *mech*, *rnd*, and *land* as the explanatory variables to analyze the impact of each variable on the green development of agriculture. The specific meaning is as follows:

*urb*: A variable that represents the level of urbanization, which is reflected by the proportion of urban population in the total population. On the one hand, urbanization can gather social resources in a certain space and accelerate the transformation and upgrading of industrial structure. On the other hand, the improvement of urbanization level can increase the emission of environmental pollution and inhibit the green development of agriculture [42,43].

*mech*: A variable that represents the level of mechanization, which is reflected by the total power of agricultural machinery. The improvement of the level of mechanization is conducive to improving the efficiency of agricultural production and has a positive impact on the green development of agriculture [27].

*rnd*: A variable that represents the investment in scientific and technological research and development, which is reflected by the scientific and technological expenditure. The investment in scientific research and development can promote technological progress, improve agricultural production efficiency, and promote the green development of agriculture [30,44].

*land*: A variable that represents the area of arable land, which is reflected by the arable land area. Arable land is an indispensable element in agricultural production. For one thing, arable land area represents the resource endowment of a certain area, which can increase crop yields and promote the green development of agriculture. For another, the area of arable land decreases year by year while the utilization and output rate of arable land are increased, which has a positive impact on the green development of agriculture.

#### *2.2. Entropy Weight Method*

In the application of comprehensive evaluation methods, there are multiple evaluation methods, which can be divided into subjective weighting method and objective weighting method when determining the weight of each index. Subjective weighting methods tend to rely on subjective judgments and lacks certain objectivity, which makes the evaluation result subjectively random. The objective method uses the value of indicators to calculate the weight without considering human factors. Entropy weight method determines each index's weight on the basis of amount of information. Among the objective methods, entropy weight method has been widely used in many fields, which will make the evaluation more accurate when dealing with large amounts of data.

Therefore, we selected the entropy weight method of objective weighting methods to determine the weight of various indicators through information entropy, and then made a comprehensive evaluation of the level of agricultural green development in China. Data selected were panel data containing both time series and cross-section. In order to guarantee the comparability among different years, time variable was introduced to improve the entropy weight method. The formula of improved entropy weight method is expressed as Equations (1)–(8) [45].


Standardization of positive indicators:

$$Z\_{\lambda ij} = \frac{\mathbf{x}\_{\lambda ij} - \mathbf{x}\_{\text{min}}}{\mathbf{x}\_{\text{max}} - \mathbf{x}\_{\text{min}}} \tag{1}$$

Standardization of negative indicators:

$$Z\_{\lambda \text{ij}} = \frac{\mathbf{x\_{\text{max}}} - \mathbf{x\_{\lambda \text{ij}}}}{\mathbf{x\_{\text{max}}} - \mathbf{x\_{\text{min}}}} \tag{2}$$

Step 3 Normalization of indicators, as is shown in Equation (3):

$$P\_{\lambda ij} = \frac{Z\_{\lambda ij}}{\sum\_{\lambda=1}^{h} \sum\_{i=1}^{n} Z\_{\lambda ij}} \tag{3}$$

Step 4 Calculate the entropy of each indicator, as is shown in Equations (4) and (5):

$$E\_j = -k \sum\_{\lambda=1}^{h} \sum\_{i=1}^{n} P\_{\lambda ij} \ln P\_{\lambda ij} \tag{4}$$

$$k = \frac{1}{\ln(h \times n)}\tag{5}$$

Step 5 Calculate the redundancy of the entropy values of each indicator, as is shown in Equation (6):

$$D\_{\bar{j}} = 1 - E\_{\bar{j}} \tag{6}$$

Step 6 Calculate the weight of each indicator, as is shown in Equation (7):

$$\mathcal{W}\_{\dot{j}} = \frac{D\_{\dot{j}}}{\sum\_{j=1}^{m} D\_{\dot{j}}} \tag{7}$$

Step 7 Calculate the level of agricultural green development in each province in each year, as is shown in Equation (8).

$$\mathbb{C}\_{\lambda i} = Z\_{\lambda i j} \rtimes \mathcal{W}\_{\dot{j}} \tag{8}$$

#### *2.3. The Spatial Autocorrelation Analysis Method*

#### 2.3.1. Global Spatial Autocorrelation

Before studying the spillover effects of agricultural green development level in China, it is necessary to test whether there exits statistical spatial correlation among the provincial agricultural green development level. It is common in the literature to measure the global spatial autocorrelation with Moran s *I* index. The formula of Moran s *I* is calculated in Equation (9) [46].

$$\text{Morran's } I = \frac{n}{\sum\_{i=1}^{n} (\mathbf{x}\_i - \overline{\mathbf{x}})^2} \frac{\sum\_{i=1}^{n} \sum\_{j=1}^{n} w\_{ij} (\mathbf{x}\_i - \overline{\mathbf{x}}) (\mathbf{x}\_j - \overline{\mathbf{x}})}{\sum\_{i=1}^{n} \sum\_{j=1}^{n} w\_{ij}} \tag{9}$$

where *xi* denotes the observed value of the agricultural green development level in province *i*; *xj* denotes the observed value of the agricultural green development level in province *j*; *x* denotes the mean of observed values of the agricultural green development level; *n* denotes the total provincial number; *wij* denotes the spatial weight matrix. We used 0-1 adjacent space weight matrix as the spatial weight matrix in this article. The definition is shown in Equation (10).

$$w\_{ij} = \begin{cases} 1, \text{ prove since } i \text{ adjacent to proveine } j \\ 0, \text{ otherwise} \end{cases}, i = 1, 2, \dots, n; j = 1, 2, \dots, n \tag{10}$$

Moran s *I* index ranges from −1 to 1. The Moran s *I* index shows a positive spatial autocorrelation if it is between 0 and 1, which represents that the "high value" is adjacent to the "high value", and the "low value" is adjacent to the "low value". The Moran s *I* index shows a negative spatial autocorrelation if it is between −1 and 0, which represents that the "high value" is adjacent to the "low value". There is no spatial autocorrelation if the Moran s *I* index is 0 [47].

#### 2.3.2. Local Spatial Autocorrelation

The meaning of the local Moran index is similar to that of the global Moran index, which is expressed as Equation (11). If Moran s *Ii* > 0, it means that the "high(low) value" in region *i* is surrounded by the "high(low) value". Otherwise, if Moran s *Ii* < 0, it means that the "high(low) value" in region *i* is surrounded by the "low(high) value" [47].

$$\text{Moran's } I\_i = \frac{n^2}{\sum\_{i=1}^n (\mathbf{x}\_i - \overline{\mathbf{x}})^2} \frac{(\mathbf{x}\_i - \overline{\mathbf{x}}) \sum\_{i=1}^n \sum\_{j=1}^n w\_{ij} (\mathbf{x}\_j - \overline{\mathbf{x}})}{\sum\_{i=1}^n \sum\_{j=1}^n w\_{ij}} \tag{11}$$

## *2.4. Spatial Panel Data Model*

Spatial econometric models mainly include Spatial Durbin Model (SDM), Spatial Autoregressive Model (SAR), and Spatial Error Model (SEM). The SDM model is more general, which is a general model including SAR model and SEM model. The model is expressed as Equation (12) [44].

$$Y = \rho WY + X\beta + \theta WX + \varepsilon \tag{12}$$

where ρ represents the spatial autocorrelation coefficient, *W* represents the spatial weight matrix, *X* represents the explanatory variables, *WX* represents the spatial lag of explanatory variables, β represents the regression coefficient of the explanatory variables, and ε represents the random error term.

When θ = 0, the SDM model will be simplified to the SAR model. The model is expressed as Equation (13).

$$Y = \rho \mathcal{W}Y + X\beta + \varepsilon \tag{13}$$

When θ + ρβ= 0, the SDM model will be simplified to the SEM model. The model is expressed as Equations (14) and (15).

$$Y = X\beta + u\tag{14}$$

$$
\mu = \lambda \mathcal{W} \mu + \varepsilon \tag{15}
$$

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

#### *3.1. Calculation of Agricultural Green Development Level*

According to the basic principle of entropy weight method, the article calculated the index weight of each index step by step, and then calculated the level of agricultural green development of 31 provinces (cities and municipalities) from 2007 to 2018 through weighted calculation. The results are shown in Table 2.

It can be seen from Table 2 that from 2007 to 2018, the overall level of agriculture green development in China shows a slow upward trend. The average national level of agricultural green development fluctuated from 0.429 to 0.468, with an average annual growth rate of 0.76%. Provincial level of the agricultural green development shows an upward trend, but there exist differences among provinces. Among the top 10 provinces in terms of average level are 8 provinces in the eastern region, while among the bottom 10 provinces in terms of average level are 7 provinces in the western region. It can be seen that the provinces with a higher level of agricultural green development in China are concentrated in the eastern area, and the provinces with a lower level of development are concentrated in the western region. The three major regions have significant regional differences in the level of agricultural green development, showing the spatial characteristics of Eastern > Central > Western on the whole. In addition, the range of provincial agricultural green development level decreased from 0.341 to 0.274, indicating that the gap of agricultural green development level among provinces in China is decreasing.



#### *3.2. Spatial Correlation Analysis*

#### 3.2.1. Global Autocorrelation Test

The results of the Moran s *I* index are shown in Table 3. It shows that the Moran s *I* index of provincial agricultural green development level in China is positive at the significance of 5% from 2007 to 2018. It indicates that provincial agricultural green development has an obvious positive autocorrelation in the spatial distribution. The level of provincial agricultural green development presents a spatial agglomeration effect. A positive spatial correlation means that the characteristics of adjacent provinces are similar, that is, provinces with higher level tend to be close to those with higher level, or provinces with lower level tend to be relatively close to those with lower level. Therefore, there is a spatial correlation between the provincial level of agricultural green development in China on the whole, which means that the level of agricultural green development shows significant spatial agglomeration characteristics in space.

**Table 3.** Global Moran's *I* of AGD level during 2007–2018.


#### 3.2.2. Local Correlation Test

In order to further test the spatial relationship of agricultural green development in China and investigate the distribution characteristics of agricultural green development in local areas, the article draws the local Moran index scatter plot of some years (2007, 2011, 2015, 2018) based on the 0–1 adjacent spatial weight matrix, as is shown in Figure 1a–d. The presentation of Moran scatter plot is to divide the space into four quadrant units, and different quadrant units represent the relationship between the local spatial unit and its adjacent spatial unit. Specifically, the first quadrant is a high-high (H-H) agglomeration area, which means that the region with a high level is surrounded by regions with high level. The second quadrant is a low-high (L-H) agglomeration area, which means that the region with a low level is surrounded by regions with high level. The third quadrant is a low-low (L-L) agglomeration area, which means that the regions with a low level is surrounded by regions with low level. The fourth quadrant is a high-low (H-L) agglomeration area, which means that the region with a high level is surrounded by regions with low level.

**Figure 1.** Local Moran index scatter plot of AGD level.

In order to understand the agglomeration or dispersion of provinces more intuitively, the data of each quadrant were summarized and presented in the form shown in Table 4. From the perspective of time dynamics, the level of green agricultural development in various provinces has changed, but in the selected four years, most of the provinces are located in the H-H and L-L agglomeration areas. From a long-term perspective, the provinces located in each quadrant have relatively little change. Provinces located in the H-H agglomeration area for a long time are mostly provinces in the eastern coastal area, while the provinces located in the L-L agglomeration area are mostly provinces in the western region. It indicates that agricultural green development in China exists as a local spatial correlation, and the level of agricultural green development in different provinces exists spatially heterogeneously.




**Table 4.** *Cont.*

### *3.3. Spatial Econometric Analysis*

#### 3.3.1. Model Test and Selection

Before the spatial econometric regression, it is necessary to analyze the panel data with the ordinary least squares (OLS) model. First, we used F test to compare the fixed effect model with the mixed effect model without considering spatial correlation, and the result showed that the fixed effect model was better than the mixed effect model. Then, we used the Breusch and Pagan (B&P) test to compare the random effect model and the mixed effect model, and the result showed that the random effect model was better than the mixed effect model. Finally, we used Hausman test to compare the random effect model and the fixed effect model, and the result showed that the fixed effect model was better than the random effect model. The test results are shown in Table 5. Through the above three tests, we found that the fixed effect model is more suitable. Therefore, the fixed effect model was used for further spatial econometric model selection to analyze the influencing factors of agricultural green development level in China.

**Table 5.** Panel data related effect diagnosis results.


Considering the spatial correlation, we used LM test and Robust LM test to compare the Spatial Autoregressive Model (SAR) with the Spatial Error Model (SEM), as is shown in Table 6. It can be seen from the LM results that the LM statistics of both the SAR model and the SEM model are at the significance level of 5%, rejecting the null hypothesis that there is no spatial lag term and spatial error term. The result indicates that the two models can be used to analyze the influencing factors of agricultural green development. From the perspective of robust LM test results, the SAR model failed to pass the significance test, indicating that the SEM model is better than the SAR model.


**Table 6.** LM test results.

In order to make the estimation results more accurate, based on the fixed effect model estimation results of Spatial Durbin Model (SDM), we used Wald test and likelihood ratio (LR) test to determine whether the SDM model with fixed effect needs to be simplified to SAR model or SEM model. The test results are shown in Table 7, where the null hypothesis of test 1 is θ= 0, and the null hypothesis of test 2 is θ= −ρβ. From the test results, we can see that the statistics of the two tests are both at the significance of 1% and reject the null hypothesis. The results mean that there may be bias when only using the SAR model or the SEM model to analyze the spatial spillover effect of influencing factors of agricultural green development, and the SDM model is the suitable spatial econometric model.

**Table 7.** Model comparison test results.


In this article, we listed the regression results of the fixed effect model under the SDM model, the SAR model, the SEM model, and the OLS model for comparative analysis, and the results are shown in Table 8. According to the results in Table 8, compared with the SAR model, the SEM model, and the OLS model, the SDM model considering spatial lag effect terms and spatial error effect terms has the best fitting effect, indicating that spatial geographical factors cannot be ignored when analyzing the influencing factors of agricultural green development level. After comparing the significance of each variable, the significance of spatial effect, and the Log-likelihood, it is found that the spatial effect of the three spatial econometric models is at the significance of 1%. However, the independent variable coefficient of the SDM model is more significant, and the Log-likelihood value is the largest, which further confirms that the SDM model is the optimal estimation model.

In order to clarify the type of the selected model, we made a further test on the fixed effect model under the Spatial Durbin model, exploring the choice of spatial fixed effect model, time fixed effect model, or spatial and time double fixed effect model, as is shown in Table 9. According to the results in Table 9, although the fitting effect under the time fixed effect is lower than that under the spatial fixed effect, the coefficient of each variable of the estimated result under the time fixed effect is significant. Therefore, this article selected the time fixed effect model under the Spatial Durbin Model for analysis.


**Table 8.** Model comparison estimation results.

Note: Standard errors in parentheses \* *p* < 0.10, \*\* *p* < 0.05, \*\*\* *p* < 0.01.



Note: Standard errors in parentheses \* *p* < 0.10, \*\* *p* < 0.05, \*\*\* *p* < 0.01.

### 3.3.2. Results and Discussion

Due to the vast territory of China, there are obvious differences in agricultural resource endowments and supporting policies in different provinces. Therefore, there are obvious regional differences in the factors that affecting the level of agricultural green development. In order to analyze the spatial characteristics and spatial spillover effect of the factors affecting the level of agricultural green development in different regions, China can be divided into three areas based on the geographical and administrative location, i.e., the eastern, central, and western areas, as is shown in Table 10 [48]. In this study, 31 provinces (municipalities, autonomous regions) were considered exclusive of Hong Kong, Macao, and Taiwan for the data absence.



The article constructs the Spatial Durbin Model with time fixed effect on the influencing factors of agricultural green development in the nation, the eastern area, the central area, and the western area for analysis, and the estimated results are shown in Table 11.


**Table 11.** Estimation of the Spatial Durbin Model (SDM) model of influencing factors of AGD level.

Note: Standard errors in parentheses \*\* *p* < 0.05, \*\*\* *p* < 0.01.

From the national perspective, the four influencing factors of urbanization level, agricultural mechanization level, scientific and technological R&D investment, and arable land area have significant impacts on the level of agricultural green development in China at the significance level of 1%. Among the four factors, urbanization level, agricultural mechanization level, and scientific and technological R&D investment have a positive impact on the level of agricultural green development, while arable land area has a negative impact on the level of agricultural green development. Meanwhile, the factors in the adjacent regions have a significant impact on the provincial agricultural green development level. Among the four factors, the urbanization level and agricultural mechanization level in the adjacent provinces have a significant negative impact on the provincial agricultural green development level, while the scientific and technological R&D investment and arable land area in the adjacent provinces have a significant positive impact on the provincial agricultural green development level.

From the regional level, the impact of various factors on the level of agricultural green development in different regions is different, and the impact direction is also different. Therefore, we adopted the horizontal comparison method to analyze each variable in the eastern, central, and western areas to illustrate the spatial heterogeneity.

Considering the impact of the variable "*lnurb*" on the level of agricultural green development, the urbanization level has a significant impact on the green development of agriculture in the eastern, central, and western areas at a significance level of 1%, but there are differences in the impact direction. The coefficient of "*lnurb*" is significantly negative in the eastern and central areas and is significantly positive in the western area. The eastern and western areas have superior geographical location and are economically developed. The improvement of urbanization level releases the surplus labor force in rural areas to engage in the secondary and tertiary industries, which increases the emission of environmental pollution, thus restraining the green development of agriculture. The urbanization level in the western area is in the initial stage, and the improvement of urbanization level gathers some social resources, which is conducive to accelerating the promotion of the green development of agriculture.

Considering the impact of the variable "*lnmech*" on the level of agricultural green development, the impact of agricultural mechanization level on the green development of agriculture in the eastern and western areas is positive at the significance level of 5% and 1%, respectively. In the process of agricultural production, agricultural machinery can improve agricultural production efficiency, improve land productivity, save labor, and promote the green development of agriculture.

Regarding the impact of the variable "*lnrnd*" on the level of agricultural green development, the impact of scientific and technological R&D investment on the green development of agriculture in the eastern and western areas is positive at the significance level of 1%. Investment in scientific and technological research and development can effectively promote technological progress, and technological progress provides guarantee for the entire process of agricultural production, as well as improving the agricultural production efficiency in the meantime, which can promote the green development of agriculture.

With regards to the impact of the variable "*lnland*" on the level of agricultural green development, the impact of arable land area on the green development of agriculture in the eastern, central, and western areas is negative at the significance level of 1%. In recent years, with the acceleration of industrialization and urbanization, the arable land area has a trend of continuous decrease, but the increase of the utilization rate and output rate of arable land has played a significant positive role in the green development of agriculture.

Considering the impact of the variable "*W*×*lnurb*" on the level of agricultural green development, the urbanization level in adjacent provinces is significantly negative in the western area at a significance level of 1%, indicating that the regional spillover effect of urbanization on agricultural green development is negative. The improvement of urbanization level in adjacent provinces has a restraining effect on the local green development of agriculture for the western area. The possible reason is that the development of urbanization in the western area can attract agricultural surplus labor in adjacent provinces, which has led to a decline in the level of agricultural development in the local area.

Concerning the impact of the variable "*W*×*lnmech*" on the level of agricultural green development, the agricultural mechanization level in adjacent provinces is significantly negative in the eastern and central areas at a significance level of 5% and 1%, respectively, indicating that the regional spillover effect

of agricultural mechanization level on agricultural green development is negative. The improvement of agricultural mechanization level in adjacent provinces has a restraining effect on the local green development of agriculture for the eastern and central area.

Regarding the impact of the variable "*W*×*lnrnd*" on the level of agricultural green development, the scientific and technological R&D investment in adjacent provinces is significantly negative in the central area at a significance level of 5%, indicating that the regional spillover effect of scientific and technological R&D investment on agricultural green development is negative. The increase of scientific and technological R&D investment in adjacent provinces has a restraining effect on the local green development of agriculture for the central area.

Finally, with regards to the impact of the variable "*W*×*lnland*" on the level of agricultural green development, the arable land area in adjacent provinces is significantly negative in the central area at a significance level of 5%, indicating that the regional spillover effect of arable land area on agricultural green development is negative. The increase of arable land area in adjacent provinces has a restraining effect on the local green development of agriculture for the central area.

#### **4. Conclusions and Policy Implications**

#### *4.1. Conclusions*

Firstly, we evaluated the level of agricultural green development in China from 2007 to 2018, and the results showed that the level of agricultural green development in China presented a fluctuating upward trend, and there are differences among provinces. On the whole, the level of agricultural green development in the three major regions presented the characteristics of Eastern > Central > Western.

Secondly, we explored the spatial relationship of the level of agricultural green development among provinces in China. The results showed that China's provincial agricultural green development (AGD) level has an obvious positive autocorrelation in spatial distribution, showing the significant characteristics of spatial agglomeration in space

Finally, the Spatial Durbin Model was used to analyze the influencing factors of agricultural green development in China. The results showed that, from the national perspective, the urbanization level, agricultural mechanization level, scientific and technological R&D investment, and the arable land area had significant impact on the level of agricultural green development in the local and adjacent areas. From the regional perspective, in the eastern area, the four factors had a significant impact on the local green development agriculture, and the agricultural mechanization level in adjacent areas had a negative significant impact on the local green development of agriculture. In the central area, the urbanization level and arable land area had a significant impact on the local green development of agriculture, and the mechanization level, scientific and technological R&D investment, and arable land area in adjacent areas had a negative significant impact on the local green development of agriculture. In the western area, the four factors had a significant impact on the local green development of agriculture, and the urbanization level in adjacent areas had a significant negative impact on the local green development of agriculture.

### *4.2. Policy Implications*

On the basis of the above results, several relevant policy implications can be derived, as follows:

First, agricultural green development is a systematic and complex issue. Therefore, China should introduce an authoritative and scientific agricultural green development level evaluation system and strengthen the monitoring of various data such as agricultural resources and the environment. In this way, China can make dynamic evaluation of agricultural green development in various provinces, which is helpful for the timely issue policies to promote the green transformation of agriculture. Second, China should accelerate the integration of urban and rural areas, promote the integrated development of the three major industries, and optimize the structure; in addition, paying more

attention to the preference for agriculture, so as to realize green efficiency in agricultural production and stable growth of farmers' income.

Third, China should increase the scientific and technological R&D investment in the whole process of agricultural production, including agricultural machinery, agricultural intelligent equipment, and water-saving and energy-saving equipment, etc.

Fourth, the stability of the amount of arable land resources is the most basic element to ensure national food security. It is necessary to establish a strict protection system to protect the red line of 1.8 billion mu of arable land. While ensuring the quantity of arable land resources, the quality of soil should also be improved to ensure the stability and quality of arable land.

### **5. Limitations and Future Work**

There are some limitations in this article. First, this study focuses on the estimation of the agricultural green development and its influencing factors in China based on the static Spatial Durbin Model, but the long-term dynamic association may be different or more complex. Therefore, the dynamic Spatial Durbin Model with time series is also worth discussing. Second, the spatial heterogeneity of China's agricultural green development is only divided into three regions of the eastern, the central, and the west, therefore there may be a spatial spillover effect bias problem in this article. If the geographical regions of China are divided into more regions, the influence of various factors on different regions may be changed.

In the future research, additional aspects should be considered. First, we will incorporate the time series factor into the model and discuss the time and spatial heterogeneity of agricultural green development in China with the dynamic panel data under the consideration of time and space. Second, we will subdivide the geographical regions of China from three major regions into six geographical regions to make the research conclusions more accurate.

**Author Contributions:** Conceptualization, H.G. and C.P.; methodology, H.G. and S.X.; software, H.G. and S.X.; validation, H.G. and S.X.; formal analysis, H.G.; resources, H.G.; data curation, C.P.; writing—original draft preparation, H.G.; writing—review and editing, S.X. and C.P.; visualization, H.G.; supervision, C.P.; project administration, H.G.; funding acquisition, H.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Social Science Fund Project of Jilin Province, China, grant No.2017B032; MOE (Ministry of Education in China) Project of Humanities and Social Sciences, grant/award No.18YJC630128; Social Science Fund Project of Jilin Province, China, grant No.2018BS33; and Social Science Fund Project of "the 13th Five-Year" of Education Department of Jilin Province, China, grant No.JJKH20190736SK.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
