Next Article in Journal
Land Subsidence in a Coal Mining Area Reduced Soil Fertility and Led to Soil Degradation in Arid and Semi-Arid Regions
Next Article in Special Issue
Analysis of the Agricultural Economy and Agricultural Pollution Using the Decoupling Index in Chengdu, China
Previous Article in Journal
The Characteristics of Carbon, Nitrogen and Sulfur Transformation During Cattle Manure Composting—Based on Different Aeration Strategies
Previous Article in Special Issue
Post-Adversities Recovery and Profitability: The Case of Italian Farmers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial-Temporal Characteristics of Agriculture Green Total Factor Productivity in China, 1998–2016: Based on More Sophisticated Calculations of Carbon Emissions

1
Research Center for Economy of Upper Reaches of the Yangtse River/School of Economics, Chongqing Technology and Business University, Chongqing 400067, China
2
Background Operation Center of Credit Card, China Minsheng Bank (Chengdu), Chengdu 610100, China
3
Business School, Hohai University, Nanjing 211100, China
4
School of Mathematics and Statistics, Chongqing Technology and Business University, Chongqing 400067, China
5
School of Public Administration, Sichuan University, Chengdu 610065, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(20), 3932; https://doi.org/10.3390/ijerph16203932
Submission received: 17 September 2019 / Revised: 7 October 2019 / Accepted: 12 October 2019 / Published: 16 October 2019

Abstract

:
Environmental costs should be taken into account when measuring the achievements of China’s agricultural development, since the long-term extensive development of agriculture has caused huge environmental pollution. This study took agricultural carbon emissions as an undesired output to estimate the agricultural development efficiency in 31 provinces of China from 1998 to 2016, based on the green total factor productivity, as assessed by the slacks-based measure directional distance function and constructing the global Malmquist–Luenberger index. We measured agricultural carbon emissions in terms of five aspects: agricultural materials, rice planting, soil, livestock and poultry farming, and straw burning, and then compared the green total factor productivity index and the total factor productivity index. The study came to the following conclusions: (1) the green technology efficiency change was smaller than the technology efficiency change at first, but the gap between them is narrowing with time, such that the former is now larger than the latter; (2) the green technology efficiency was in a declining state and the green technology progress was increasing, promoting the green total factor productivity growth, from 1998 to 2016; and (3) China’s agricultural green total factor productivity increased by 4.2% annually in the east, 3.4% annually in the central region, and 2.5% annually in the west.

1. Introduction

Agriculture is a key sector in China, providing the most basic goods to meet people’s basic needs in order to keep the economy running, and the Chinese government also attaches great importance to the development of agriculture. China’s agriculture has developed rapidly in recent years and achieved great progress [1]. In 2017, China’s total agriculture output value (including crop, livestock, forestry, and fishery) reached 6.8 billion yuan, with 134.9 million hectares of arable land and 361.8 million agricultural practitioners [2]. However, the development of agriculture in China has caused serious pressure on the environment [3], because it produces a large amount of carbon emissions each year, accounting for a high proportion of all carbon emissions—much higher than the United States and the world average [4,5,6,7]. Greenhouse gas emissions (including CO2, CH4, N2O, perfluorocarbons, hydrofluorocarbons, etc.) are often measured by carbon emissions and will threaten human survival and development. Therefore, we should take environmental pollution into account to comprehensively measure the development of China’s agriculture [8].
In fact, many international organizations have begun to attach importance to the relationship between agricultural development and environmental pollution and put forward some relevant measurement frameworks for the development level of green agriculture and suggestions for promoting the development of green agriculture. The meaningful studies in this area include Green Growth Indicators for Agriculture: A Preliminary Assessment by the Organization for Economic Co-operation and Development (2014) [9].
Academia generally measures development in terms of productivity. Productivity refers to the utilization efficiency of human, material and financial resources, which reflects the influence of resource allocation, technological level, and labor force on production activities. From the initial single factor, productivity has evolved into the commonly used total factor productivity (TFP) [10,11]. In the original TFP studies, undesirable outputs such as pollutant emissions were not taken into account. With the deepening of the research, green total factor productivity (GTFP), considering non-expected output, has emerged.
In recent years, there have been many research projects on agricultural TFP, which can be divided into three categories according to the measurement method. The first category is to use the growth accounting method for measurement; this was often used in relatively early studies. Many scholars used the algebraic index method to analyze China’s agricultural productivity [12,13], and many used the Solow residual method [14,15]. These kinds of methods are not outdated, because they are relatively simple and have advantages in some current studies. For example, some scholars make good use of them to calculate TFP with micro data [16,17].
In the second category, data envelopment analysis (DEA) is used to analyze the changes of China’s agricultural productivity. The advantage of this method is that it can deal with multiple input and output variables, and it facilitates the decomposition of agricultural productivity to understand the internal growth momentum of agricultural productivity. DEA has been widely used in empirical analysis [18,19,20]. The third category of research uses stochastic frontier analysis (SFA). According to whether it is necessary to set the specific form of production function in advance, this type of analysis can be further divided into the parametric method and non-parametric method. According to whether the production frontier is affected by random factors, it can also be divided into the stochastic frontier method and deterministic frontier method [21,22]. Similarly, the research on agricultural GTFP can also be divided into four categories according to the measurement method. The first three categories are the same as for agricultural TFPs, and the only difference is that the three methods take environmental pollution as an input variable. The fourth category includes those studies considering pollution as an undesirable output to calculate the GTFP [23,24].
This paper took agricultural carbon emissions (ACEs) as an undesired output to calculate China’s agricultural green total factor productivity. It is difficult to achieve a breakthrough in the research of ACEs in terms of the calculation method. Therefore, more progress tends to be made in the selection of carbon emission sources. West et al. chose agricultural inputs, including fertilizers, agricultural limes, pesticides, agricultural irrigation and seed cultivation, as agricultural carbon emission sources (ACESs) [25]. Johnson et al. argued that ACEs are mainly derived from intestinal fermentation of livestock, manure management, rice growth, and arbitrary disposal of agricultural waste [26]. Other studies, such as those by Tan Q. and Wu Y. et al., determined the sources as having four aspects: rice planting, rice fields, agricultural material input, livestock and poultry farming, and soil management. The index is adequate, but the ACE coefficients were not listed in the paper in detail [27,28].
Like many studies focused on measuring GTFP, this paper takes ACEs as an undesired output. However, a new breakthrough was achieved in making the GTFP more closely reflect reality in the calculation of ACEs. This study calculated ACEs in China by choosing five categories of ACEs: agricultural materials, rice cultivation, N2O emissions due to the damage of the soil surface when planting crops, livestock and poultry farming, and crop straw burning. Contrary to other research, the species of rice, crops on the land and livestock and poultry, and the location of rice, were taken into account when calculating ACEs from rice planting, soil N2O, and livestock and poultry farming, because carbon emissions of those ACESs could be influenced by these factors. In addition, the study added agriculture straw combustion as one of the ACESs, which improved the estimation accuracy for ACEs.
This paper consists of four parts. The first part and the second are the introduction and method, respectively. Section 3 presents the comparison of green total factor productivity index (GTFPI)and total factor productivity index (TFPI), the temporal and spatial characteristics of GTFPI, and the analysis of the spatial correlation of GTFPI in China. The discussion and conclusion are presented in Section 4.

2. Materials and Methods

2.1. Measurement Way of Green Total Factor Productivity

2.1.1. Global Production Possibility Set

China’s 31 provinces (excluding Hong Kong, Macao and Taiwan) are all taken as decision-making units. Suppose each unit uses n kinds of input (x, xRn+) in the production process, and obtains both m kinds of desired output (y, yRm+) and undesired output—ACEs (b). The set of possibilities between inputs, desirable output, and undesired output is defined as environmental production technology [29], which combined with its bounded closure, the strong disposition of input and undesired output, weak disposition of desirable output and undesired output, and zero binding, and the production possible set (P) is expressed as Equation (1):
P ( x ) = { ( x , y , b ) : x   c a n   p r o d u c e   ( y , b ) }
Further, according to the global production possibility set proposed by OH [30], the union set of all production technology sets in the current period can be expressed as Equation (2):
P G ( x ) = P 1 ( x 1 ) P 2 ( x 2 ) P T ( x T )
where G is the global benchmark and T is the period. Under the condition that the remuneration of production technology scale remains unchanged, assume the time period t = 1, 2, ..., T and decision-making unit k = 1, 2, ..., K, the input and output vector is ( x k t + y k t + b k t ) , then the global production possibility set is expressed as:
P G ( x ) = { ( x t + y t + b t ) : t = 1 T k = 1 K Z k y k m t y m t ; t = 1 T k = 1 K Z k b k t = b t ; t = 1 T k = 1 K Z k x k n t x n t }
Z k is the density variable, which represents the weight of each of k decision-making units in the construction if the environmental technical structure.

2.1.2. Slacks-Based Measure (SBM) Directional Distance Function

The related efficiency measurement idea and models belong to radial and linear piecewise measurement theory. The main features of this measurement: strong disposable, namely its linear piecewise frontier sometimes parallel to the horizontal axis or the vertical axis, ensures the efficiency boundary or indifference curve convexity (not bending), but also causes input “congestion” or “slacks”. At the same time, the linear piecewise front obtained by linear programming also violates free disposal—the basic assumption of Neoclassical economics.
Figure 1 illustrates this problem nicely. The frontier of production technology is the linear segment frontier SS’, composed of two efficient production units at point C and point D. The points on the broken line of SS′are the efficient production points of the production technology. Therefore, the production units A and B are the inefficient production points. According to the measurement method of radial technical efficiency, the degree of technical efficiency of A and B can be expressed as OA′/OA and OB′/OB. A′ and B′ on the production front are the efficient reference points of A and B. Like C and D, A′ and B′ are points of technical efficiency 1. Comparing point A′ with point C, we find that the production at point A′ can produce the same output y as the production at point C by continuing to reduce the input of X2, and this condition is called input slacks. If the model is extended to multiple inputs and multiple outputs, there will also be output slacks. Complete efficiency requires neither inefficiency nor input slacks, which is a problem that traditional DEA models have been unable to solve.
The SBM directional distance function can consider the influence of input and output slack variables on efficiency, and it can avoid the bias of traditional radial DEA on efficiency evaluation. According to the non-radial and non-angle SBM efficiency model proposed by TONE [31], both input reduction and output increase are considered. On the basis of Equation (3), the non-radial, non-angle SBM directional distance function containing an undesired output in time t of a decision-making unit k ( x k t + y k t + b k t ) is constructed:
D 0 G ( x k t + y k t + b k t ) = m i n 1 [ 1 N n = 1 N s n x / k n k ] 1 + [ 1 M + 1 ( m = 1 M s m y / y m k + i = 1 I s i b / b i k ) ]
x k n t = t = 1 T k = 1 K Z k t x k n t + s n x , n = 1 , 2 , , N
y k m t = t = 1 T k = 1 K Z k t y k m t s m y , m = 1 , 2 , , M
b k i t = t = 1 T k = 1 K Z k t b k i t + s i b , i = 1 , 2 , , I
In Equations (4)–(7), D 0 G is the average distance between input and output and the production front—the inefficiency degree of input and output. s n x s m y and s i b are the relaxation variables of input, desired output, and undesired output, respectively, all of which are greater than or equal to 0. When s n x = s m y = s i b = 0, the decision-making unit is completely efficient and there is no surplus of input or insufficient output. Otherwise, it means that there is an efficiency loss.

2.1.3. Global Malmquist-Luenberger (GML) Index

The GML index is based on the common global frontier structure of each period, which is multiplicative and with transitivity, which can effectively overcome the defects of Malmquist-Luenberger index and avoid the phenomenon of “technical regression”.
In this paper, the GML index was used to construct the green total factor productivity index (GTFPI), which represented changes in TFP. GTFPI could be decomposed into the green technology efficiency change (GTEC) and green technical progress change (GTPC). The three indexes were established as Equations (8)–(11). If the index was greater than 1, it meant that the variables at stage t+1 were higher than that at stage t, and vice versa.
G T F P I = 1 + D 0 G ( X t , y t , b t ; y t , b t ) 1 + D 0 G ( X t + 1 , y t + 1 , b t + 1 ; y t + 1 , b t + 1 )
G T F P I = G T E C G T P C
G T E C = 1 + D 0 G ( X t , y t , b t ; y t , b t ) 1 + D 0 t + 1 ( X t + 1 , y t + 1 , b t + 1 ; y t + 1 , b t + 1 )
G T P C = 1 + D 0 t + 1 ( X t + 1 , y t + 1 , b t + 1 ; y t + 1 , b t + 1 ) 1 + D 0 G ( X t + 1 , y t , b t + 1 ; y t + 1 , b t + 1 )

2.2. The Selection of Input and Output Indicators

The input and output index system of GTFPI, measured by each province, is shown in Table 1. In this paper, input indexes were selected from eight aspects. In terms of output, the desired output was the actual output value of agriculture, forestry, husbandry and fisheries, calculated for 1997, and ACEs were selected as undesired output indexes. This paper compiled a huge ACE measurement system, which was one of its notable innovations.

2.3. Estimation of Agricultural Carbon Emissions

A measurement system was compiled to calculate ACEs. The carbon emission coefficients represented how much carbon was produced per unit by ACES, and they came from the related natural science research. This study calculated ACEs in China by choosing five categories of ACESs: agricultural materials, rice cultivation, N2O emissions due to the damage to the soil surface when planting crops, livestock and poultry farming, and crop straw burning, and each of the five categories included various different types of ACES.
The ACESs could be converted into ACEs by the Equation (12):
E = E i = Q i × a i
where E represented the total ACEs, E i represented the ACE of the i-th ACES, the quantity of the i-th ACES was represented by Q i and a i represented the carbon emission coefficient of the i-th ACES.
ACESs in this paper emitted CH4 and N2O gases, and they need be converted to standard carbon emissions before calculating ACEs. According to the Intergovernmental Panel on Climate Change (2007) report [32], in terms of the greenhouse effect, 1 ton of CH4 is equivalent to 6.8182 t of C, and one ton of N2O is equivalent to 81.2727 t of C. It could be converted according to Equation (13):
C = C C H 4 + C N 2 O = C H 4 6.8182 + N 2 O 81.2727

2.4. Agricultural Carbon Emission Sources

2.4.1. Carbon Emissions from Agricultural Materials

All research on agricultural carbon emissions considers agricultural materials as one of the ACESs. In fact, in earlier studies, agricultural materials were considered as the only one. According to existing relevant studies [32,33], fertilizer, pesticides, plastic sheeting, diesel oil, and irrigation were considered carbon sources from agricultural materials, and their carbon emission coefficients were 0.8956, 4.9341, 5.18, 0.5927 and 266.48, respectively. The units are all kg C/kg.

2.4.2. Carbon Emissions from Rice Cultivation

Rice cultivation produces CH4. China, as the country with the highest yield of rice in the world, can account for a large proportion of its ACEs coming from CH4. Other dryland crops can also emit CH4 but can partially reabsorb it themselves, so production from other sources was so small that it is ignored in this study. Taking the location and the species of rice into account, the CH4 emission coefficients are shown in Table 2 [34].

2.4.3. N2O Emissions due to Damage to the Soil Surface When Planting Crops

When crops are planted, the surface of the soil is turned over and large amounts of greenhouse gases are released into the air, most notably N2O. In addition, compared with other greenhouse gases, N2O gas has a great potential to contribute to temperature increases and a long retention time, with obvious negative effects. Although CO2 is also released, the amount is not large and some of the CO2 can be absorbed, so it is ignored in this study. Similarly, the variety of crop affects N2O emissions, so we considered six types of crops: paddy rice, winter wheat, spring wheat, soybean, corn, and vegetables, and their N2O emission coefficients were 0.24, 2.05, 0.4, 0.77, 2.532 and 4.21, respectively [35,36]. All the units are g/m2.

2.4.4. Carbon Emissions from Livestock and Poultry Farming

The manure management systems and enteric fermentation during livestock and poultry farming activities produce CH4 and N2O. Before measuring carbon emissions from livestock and poultry farming, the amount of feeding should be properly adjusted because of the different breeding cycle of each kind of livestock or poultry. This study selected the top 12 species of animals in Chinese agriculture: cows, buffalo, cattle, mules, camels, donkeys, horses, live pigs, sheep, goats, rabbits, and poultry.
Pigs, rabbits and poultry have an average life cycle of 200 days, 105 days and 55 days, respectively, and their slaughter rates are all greater than 1, so their numbers are adjusted according to Equation (14):
N i = D L i M i 365
N i represents the average annual feeding amount of livestock or poultry, D L i is the average life cycle, and M i is the annual output.
The slaughter rates of the other selected livestock and poultry are less than 1, so they are adjusted according to Equation (15):
N i = C i , t + C i , t 1 2
N i is the average annual feeding amount, and C i , t and C i , t 1 represent the year-end inventory of years   t and t 1 , respectively.
To calculate carbon emissions from livestock and poultry farming, this paper includes 12 kinds of animals that are the most cultivated in Chinese agriculture, and their CH4 and N2O emission coefficients are listed in Table 3 [33].

2.4.5. Carbon Emissions from Straw Burning

Straw burning, which is very polluting to the environment, has gradually attracted the attention of scholars and is one category of carbon emission sources in this paper. The six main crop straws, rice, wheat, corn, rapeseed, soybean and cotton, were selected as carbon sources of straw burning. According to the existing research [37], the carbon emission coefficients of rice, wheat, corn, rape, soybean and cotton straw are 0.18, 0.16, 0.17, 0.22, 0.15 and 0.13, respectively, and all the units are kg C/kg.

2.5. Sample Selection and Data Sources

The sample includes data from 31 provinces in China (excluding Hong Kong, Macao and Taiwan) from 1997 to 2016, and it was obtained from the China Rural Statistical Yearbook, China Agricultural Yearbook, China Agricultural Statistics, China Animal Husbandry Yearbook, and Statistical Yearbook of China.

2.6. Process

This paper initially calculated ACEs of 31 provinces in China from 1997 to 2016 through the agricultural carbon emission measurement system. Taking the calculated ACEs as the undesired output, GTFPI was calculated by combining the SBM model with the global Malmquist–Luenberger index, according to the input and output index system. In addition, this paper calculated the traditional agricultural total factor productivity index (TFPI); that is, only the actual agricultural, forestry, husbandry and fishery output values were used as the desired output index, without undesired output.
This paper compared GTFPI and its decomposition terms with TFPI and its decomposition terms to analyze the impact of carbon emissions when calculating agricultural TFP, and then analyzed the spatial-temporal characteristics and spatial effects of GTFPI.

3. Results

3.1. Comparison of GTFPI and TFPI

During the sample period, GTFPI was substantially higher than the TFPI. The gap between GTFPI and TFPI was volatile, and it is worth noting that it increased significantly in 2015 and 2016. During the entire sample period, the average annual growth rate of GTFP was 4.07%, and the TFP was only 3.03%.
In order to better analyze the differences in the two indexes, this paper divided the sample period into two stages: 1998–2006 and 2007–2016, then compared them based on the three major regions—east, center, and west, as shown in Figure 2. The differences in GTFPI and TFPI in the three regions were nearly positive, and the average difference between them in 2007–2016 was mostly greater than that in 1998–2006. This indicated that GTFP had a higher growth rate, and this gap continued to widen in the second half of the period, which was the same conclusion as found in many previous studies. In the latter half of the period, the GTFPI in Hebei in the east, Hunan in the central area, and Tibet in the west were much larger than the TFPI. In contrast, the GTFPI in Shanghai and Ningxia were much smaller than the TFPI. We still compared the decomposition terms of the GTFPI and TFPI. It can be seen that the GTEC was smaller than the technology efficiency change (TEC) from 1997 to 2007 and larger than the TEC from 2008 to 2016, and the GTPC was greater than the technical progress change (TPC) at all times, but the gap widened further in the later time period.
Figure 3 shows the difference values between the GTFPI and TFPI, GTEC and TEC, and GTPC and TPC in 1998, 2007, and 2016 in 31 provinces nationwide. First, the dark areas in Figure 2a gradually increased, representing the difference between the GTFPI and TFPI increasing with time, which was the same as found in the analysis above. Many provinces were in the −0.092 to 0.01 stage in 1998, among which the GTFPI of Beijing, Shanghai, Sichuan and Guizhou was lower than the TFPI, and in 2007, the gaps of the four provinces became positive. In 2016, only Qinghai’s GTFPI was lower than its TFPI. Among all the provinces, the gaps of Hebei, Henan and Tibet saw a big increase. For example, Hebei increased from 0.004 in 1998 to 0.649 in 2016, Henan from 0.003 in 1998 to 0.287 in 2016, and Tibet from 0 in 1998 to 0.356 in 2016. In Figure 2b, similarly, the dark areas gradually expanded over time. However, the gaps between the GTEC and TEC in most areas were still negative. In some provinces, the difference gradually narrowed to 0 and eventually turned positive, including Hebei, Heilongjiang, Jiangxi, Hubei, Hunan, Chongqing and Xinjiang. In Figure 2c, except Liaoning, Guangdong, Shanxi, Heilongjiang and Xinjiang, the overall difference between the GTPC and TPC increased gradually. It is worth noting that in 1998, areas with a smaller GTPC than TPC, such as Yunnan, Guizhou, Sichuan and Qinghai in southwest China, were transformed into areas with a greater GTPC than TPC in 2016. In 2016, only Chongqing had a negative balance, and the gap did not exist in Beijing, Tianjin, Liaoning, Shanghai, Jiangsu and Hainan; their GTPC was equal to their TPC.

3.2. Analysis of Temporal and Spatial Characteristics of China’s Agricultural GTFPI

In the sample cycle, China’s agricultural GTFP showed a trend of continuous growth, with an average annual growth of 3.5%. The GTEC was mostly below 1, indicating green efficiency having an overall downward trend, while the GTPC was only less than 1 in 2008, and more than 1 in other years. According to the decomposition results, the growth of China’s agricultural GTFP mainly depended on agricultural green technology progress (GTP), and its contribution rate reached 4.6%. The average GTEC was less than 1, offsetting some of the growth from green technological progress.
In this paper, the sample period was divided into three periods: 1998–2003, 2004–2009, and 2010–2016, shown in Table 4. In any time period, the GTFP in China and the eastern central and western regions of China showed an increasing trend. The GTEC was less than 1 in the whole sample period and nationwide, showing that green technology had a decreasing trend. It should be noted that green technical efficiency (GTE) increased on average between 2004 and 2009 in the east and decreased on average between 1998 and 2003 and between 2010 and 2016. However, GTE in the east gradually increased on average during the whole sample period. GTP increased in three time periods and all regions, which showed that the improvement of China’s agricultural GTFP was mainly driven by GTP.
Taking 1998 as the base period, GTFP, GTE and GTP were all set as 1. We calculated their cumulative value over the sample period. During the period 1998–2003, TFP increased only slightly. As GTE declined, GTP did not fully offset the negative effects, because although China’s agriculture had a high yield in that period, due to market problems and agricultural structural contradictions, the grain could not be sold, which seriously affected the enthusiasm of Chinese farmers in production. From 2004 to 2009 the GTFP grew slowly, with an average annual growth rate of 2.6% which was mainly because the Chinese government has introduced a number of reform policies for rural areas, including abolishing agricultural taxes and subsidizing farmers. From 2010 to 2016, the GTFP showed an annual growth rate of 7.2%. That was because during this period, the Chinese government continued to increase its expenditure on agriculture and vigorously promoted agricultural technology, increasing the role of agricultural technology in promoting GTFP. Throughout the sample period, GTE continued to decline, while GTP continued to increase. This phenomenon showed that the developed technology has not been well used in practice, has not formed a demonstration effect, and has not been widely converted into productivity.
The average agricultural GTFPI in each province was consistent with the national level throughout the sample cycle. Agricultural GTFP was increasing, GTE was decreasing, and GTP was increasing. The agricultural GTFP increase was mainly driven by agricultural GTP. However, there were some differences in the speed of agricultural GTFPI in eastern, central and western China, in that growth was fastest in the east, followed by the central region, and was slowest in the west. The GTFP increased by 4.2% annually in the east region, 3.4% annually in the center, and 2.5% annually in the west.
The study divided the sample period into three periods: 1998–2003, 2004–2009 and 2010–2016, as shown in Figure 4. From the three periods, the fluctuation trend of each sub-period was highly consistent, and the fluctuation range was different. In most provinces, the overall productivity growth rate was above 1 in each sub-period, realizing the growth effect. It is worth noting that Gansu, Qinghai and Ningxia in the western region showed negative growth during the period from 2010 to 2016.
In this paper, 1998, 2004, 2010 and 2016 were selected as observation points for the distribution maps of the GTFPI of various provinces in China, as shown in Figure 5. It can be seen from the figure that some areas have fluctuated greatly over time. Except for some provinces, the overall dark areas gradually expanded, and the GTFPI as a whole was in an incremental state. From the perspective of the spatial dimension, the spatial regional difference was relatively small in 1998, and expanded in 2016. It can be seen that in 2010, the trend of agricultural GTFPI accumulation was relatively obvious, but in the other three years, the high level agricultural GTFPI provinces did not exert an agglomeration effect.
By analyzing the temporal and spatial characteristics of the GTFPI and its decomposition terms, some main conclusions are obtained in this paper: GTE was in decline and GTP was increasing, promoting GTFP growth from 1998 to 2016; and China’s agricultural GTFP is the fastest in the east, followed by the center, and the slowest in the west. These findings were mostly in line with initial expectations. The emergence of the phenomenon of GTE being in decline and GTP increasing was because of a scientific research supply and actual demand disconnection, scientific and technological achievements being idle, and a low conversion rate. Investment in agricultural research and development and the extension of agricultural technology in China have been continuously promoted. However, the transformation of agricultural scientific research results is not timely, and the extension system of agricultural technology and the socialized service system of agriculture responsible for any transformation do not perform effectively.

3.3. The Spatial Correlation Analysis of the GTFPI

Agricultural TFP is determined by a variety of factors. A similar resource endowment and agricultural production environment would inevitably lead to similarity in agricultural GTFP. On the premise of such similarity, the cost of factor flow between neighboring provinces is low, which is beneficial to factor flow and technology diffusion and improved the spatial dependence between provinces. Through the above figure, it was found that GTFPI showed an agglomeration effect among provinces in 2010. Here, it was predicted that GTFPI would have a spatial interaction among provinces.
To test this hypothesis, this paper used the Moran index to measure the spatial correlation of the GTFPI. To ensure the reliability of the conclusion, the first-order adjacency weight matrix and the reciprocal square of distance were used to calculate the Moran index, as shown in Table 5. More than half of the years in which the indexes were calculated using the two weight matrices were not significant, and those years in which the index values were significant were also very small. Therefore, the hypothesis that GTFPI has spatial effects among provinces was not valid.
The results differed from those of many other studies. The first reason for this was that this paper adopted a more reasonable way to calculate ACEs, which was closer to reality. Taking ACEs as an undesired output, the measured GTFPI was naturally different from other studies. Second, compared with the DEA-Malmquist model and SFA model adopted in many studies, the SBM-GML model adopted in this paper made the calculation results of GTFPI more reasonable. Therefore, we believe that such results can better reflect China’s actual national conditions.
First, agriculture is a sector that uses natural resources as production objects, and its production depends largely on the natural endowments of the region. Admittedly, the development of neighboring areas can indirectly affect the agricultural production of the region. However, China has a vast territory and diverse topography. Climate conditions in various regions are also very different, and as this paper analyzed at the provincial level, there are big differences between provinces in resource endowment. Therefore, the neighboring areas have a limited influence, particularly in areas such as Sichuan and Tibet, Gansu, Qinghai, Yunnan, and other provinces bordering them. Due to climatic and geographical reasons, Sichuan is mainly engaged in traditional agriculture, while Tibet is more focused on animal husbandry. The crops cultivated in Gansu, Qinghai and other provinces are also very different from those grown in Sichuan, and the agricultural development situation among provinces is diverse. Therefore, their agricultural production is relatively independent. Second, Chinese provincial governments have the right to formulate a variety of agricultural development policies for their province. Accordingly, this has caused the development of agriculture in each province to be independent to a certain extent.

4. Discussion

Although there has been a lot of research on China’s agricultural green total factor productivity, there are still two points to be improved.
The first relates to undesired output. Most studies take agricultural carbon emissions as an undesired output, but when calculating agricultural carbon emissions, the selection of carbon emission sources is not fully considered. In the early days, a large proportion of research focused on agricultural material as the sole source of carbon emissions. With further development of this area of research, some scholars have added carbon emissions generated by rice to carbon emission sources, but still did not consider that the type and location of rice would also affect carbon emissions. Most scholars did not include soil N2O, as a source of carbon emissions, and those who did failed to take into account that crop varieties grown on the soil also affect N2O emissions. The studies that have looked into carbon emissions from livestock and poultry farming are not comprehensive enough in the selection of the species. Research measuring the agricultural GTFP has tended to focus on the measurement part, and there has not been enough consideration of the choice of ACESs. In recent years, there have been some relatively comprehensive studies, but these still have also omitted some important sources.
The second is the GTFP calculation method. In some studies, ACEs were taken as an input variable and GTFP was calculated with the general DEA model or SFA model. This was reasonable to some extent, but not rigorous enough. By contrast, the SBM directional distance function is the most effective method. However, the distance function fails to effectively deal with the inconsistency of the production frontier of production units in each period, which affects the comparability of inter-period results.
Therefore, this study intended to make breakthroughs in the above two aspects. Firstly, the paper would take the following carbon emission issues into account. The variety (early rice, late rice and in-season rice) and location of the rice affects its carbon emissions. In addition, the types of crops that grow in the soil would also influence the N2O emissions, and different types of livestock and poultry had different coefficients. Furthermore, straw combustion was included into the ACESs to ensure more accurate ACEs. Secondly, the paper adopted the GML index based on the SBM directional distance function to calculate GTFP, because it can not only deal with radial and angle problems causing the GML index to bias effectively, but also realize the global comparability of the production frontier.
Our research leads to the following main conclusions. Firstly, GTFP and TFP were both increasing annually, but the former was growing faster than the latter, and the difference widens later. The GTEC was smaller than the TEC at first, but the gap between them was narrowing with time, until the former became larger than the latter. The GTPC was bigger than the TPC, and the gap had widened over time. Secondly, GTE was in decline and GTP was increasing, promoting GTFP growth from 1997 to 2016. China’s agricultural GTFP increased by 4.2% annually in the east, 3.4% annually in the central region, and 2.5% annually in the west. Thirdly, there was no spatial correlation between China’s provinces in the GTFP index.

5. Conclusions

This paper measured agricultural carbon emissions in terms of five aspects—agricultural materials, rice planting, soil N2O, livestock and poultry farming, and straw burning—and then compared the green total factor productivity index and the total factor productivity index. There are still some deficiencies in this research, however. First, although the selection of ACESs in this paper is more comprehensive and in-depth, it can still only represent the real ACEs to an extent. Due to the progress of related natural sciences, there are still many ACEs that we cannot measure. Secondly, straw burning is taken as one of ACESs, which indeed promotes the research of ACEs measurement. Straw burning has been banned by some local governments in recent years, however, due to policy lag and lax supervision, the actual effect is not obvious, which is the reason for this paper including it in the ACESs. However, the actual quantity of straw burning in each province was not investigated on the spot in this study, so the quantity calculated will be different to the actual value.

Author Contributions

Conceptualization, X.X.; methodology, X.G.; software, X.H.; formal analysis, X.X.; resources, L.C.; data curation, J.H.; writing—original draft preparation, X.X.; writing—review and editing, X.H.; visualization, J.H.; supervision, X.H.; project administration, X.X.; funding acquisition, X.X.

Funding

This research was funded by the MOE Project of Key Research Institute of Humanities and Social Sciences of Research Center for Economy of Upper Reaches of the Yangtse River “Research on agricultural modernization and industrial innovation and development”, grant number CJSYTD201710; National social science foundation of China, grant number 13BGL097; Open subject of Collaborative Innovation Center for Urban Industries Development in Chengdu-Chongqing Economic Zone ”Study on spatiotemporal differences and influencing factors of low carbon agricultural productivity in China”, grant number 1901598; Chongqing Graduate Innovative Research Project, grant number CYS19330, and Chongqing Technology and Business University Graduate Innovative Research Project, grant number yjscxx2019-101-41.

Acknowledgments

We thank Timothy Kyng at Macquarie University for his thoughtful guidance. We also thank Liu Chengjie and Chen Jiping at Chongqing Technology and Business University for their kind help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, J.; Yang, G. Understanding recent challenges and new food policy in China. Glob. Food Secur. 2017, 12, 119–126. [Google Scholar] [CrossRef]
  2. Sheng, Y.; Tian, X.; Qiao, W.; Peng, C. Measuring agricultural total factor productivity in China: Pattern and drivers over the period of 1978–2016. Aust. J. Agric. Resour. Econ. 2019. [Google Scholar] [CrossRef]
  3. Huang, X.; Xu, X.; Wang, Q.; Zhang, L.; Gao, X.; Chen, L. Assessment of Agricultural Carbon Emissions and Their Spatiotemporal Changes in China, 1997–2016. Int. J. Environ. Res. Public Health 2019, 16, 3105. [Google Scholar] [CrossRef] [PubMed]
  4. Angelakoglou, K.; Gaidajis, G. A review of methods contributing to the assessment of the environmental sustainability of industrial systems. J. Clean. Prod. 2015, 108, 725–747. [Google Scholar] [CrossRef]
  5. Xu, X.; Xu, Z.; Chen, L.; Li, C. How Does Industrial Waste Gas Emission Affect Health Care Expenditure in Different Regions of China: An Application of Bayesian Quantile Regression. Int. J. Environ. Res. Public Health 2019, 16, 2748. [Google Scholar] [CrossRef]
  6. Tian, Y.; Li, B.; Zhang, J.B. Research on spatial temporal characteristics and factor decomposition of agricultural carbon emission based on input angle–taking Hubei province for example. Res. Agric. Mod. 2011, 32, 752–755. (In Chinese) [Google Scholar]
  7. Bai, Y.; Deng, X.; Jiang, S.; Zhao, Z.; Miao, Y. Relationship between climate change and low-carbon agricultural production: A case study in Hebei Province, China. Ecol. Indic. 2019, 105, 438–447. [Google Scholar] [CrossRef]
  8. Zhou, Y.; Xu, Y.; Liu, C.; Fang, Z.; Fu, X.; He, M. The Threshold Effect of China’s Financial Development on Green Total Factor Productivity. Sustainability 2019, 11, 3776. [Google Scholar] [CrossRef]
  9. Organization for Economic Co-Operation and Development (OECD). Green Growth Indicators for Agriculture a Preliminary Assessment; Green Growth Studies OECD: Paris, France, 2014. [Google Scholar]
  10. Wu, Y. Research on Agricultural Productivity in China under Low Carbon Emission Constraints: Perspective Based on Spatial Econometrics. Ph.D Thesis, China Agricultural University, Beijing, China, 2018. [Google Scholar]
  11. Li, W. Total Factor Productivity of China’s Agriculture Based on DEA-Malmquist Index Analysis of Time and Space Difference and Influencing Factors. J. Chongqing Technol. Business University Soc. Sci. Ed. 2018, 20, 96–103. (In Chinese) [Google Scholar]
  12. Wong, L.F. Agricultural Productivity in the Socialist Countries; Westviews Press: Boulder, CO, USA, 1986. [Google Scholar] [CrossRef]
  13. Wu, S.; Walker, D.; Devadoss, S. Productivity Growth and its Component in Chinese Agriculture after Reforms. Rev. Dev. Econ. 2001, 5, 375–391. [Google Scholar] [CrossRef]
  14. McmiIlan, J.; Whalley, J.; Zhu, L. The Impact of China’s Economic Reforms on Agricultural Productivity Growth. J. Political Econ. 1989, 97, 781–807. [Google Scholar] [CrossRef]
  15. Lambert, D.K.; Parker, E. Productivity in Chinese Provincial Agriculture. J. Agric. Econ. 2010, 49, 378–392. [Google Scholar] [CrossRef]
  16. Baldoni, E.; Coderoni, S.; Esposti, R. The Complex Farm-Level Relationship between Environmental Performance and Productivity. The Case of Carbon Footprint of Lombardy farms. Environ. Sci. Policy 2018, 89, 73–82. [Google Scholar] [CrossRef]
  17. Baldoni, E.; Coderoni, S.; Esposti, R. The productivity and environment nexus through farm-level data. The Case of Carbon Footprint applied to Italian FADN farms. Bio-Based Appl. Econ. 2017, 6, 119–137. [Google Scholar] [CrossRef]
  18. Coelli, T.J.; Rao, D.S. Total factor productivity growth in agriculture: A Malmquist index analysis of 93 countries, 1980–2000. Agric. Econ. 2005, 32, 115–134. [Google Scholar] [CrossRef]
  19. Coelli, T.J.; Rao, D.S. Implicit Value Shares in Malmquist TFP Index Numbers; CEPA Working Papers No. 4/2001; School of Economics, University of New England: Armidale, Australia, 2001; p. 27. [Google Scholar]
  20. Huang, Z.; Fu, Y.; Liang, Q.; Song, Y.; Xu, X. The Efficiency of Agricultural Marketing Cooperatives in China’s Zhejiang Province. Manag. Decis. Econ. 2013, 34, 272–282. [Google Scholar] [CrossRef]
  21. Bayarsaihan, T.; Coelli, T. Productivity growth in pre- 1990 Mongolian agriculture: Spiraling disaster or emerging success. Agric. Econ. 2003, 28, 121–137. [Google Scholar] [CrossRef]
  22. Fang, W.; Yang, Z.Y.; Liang, J.F. Technical efficiency evaluation and improvement strategy of Guangdong agricultural modern demonstration area based on SFA model. J. South. Agric. 2018, 49, 1249–1255. (In Chinese) [Google Scholar] [CrossRef]
  23. Le, T.L.; Lee, P.; Peng, K.C.; Chung, R.H. Evaluation of total factor productivity and environmental efficiency of agriculture in nine East Asian countries. Agric. Econ. Czech 2019, 65, 249–258. [Google Scholar] [CrossRef]
  24. Dong, G.; Wang, Z.; Mao, X. Production efficiency and GHG emissions reduction potential evaluation in the crop production system based on emergy synthesis and nonseparable undesirable output DEA: A case study in Zhejiang Province, China. PLoS ONE 2018, 13, e0206680. [Google Scholar] [CrossRef]
  25. West, T.O.; Marland, G. A synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: Comparing tillage practices in the United States. Agric. Ecosyst. Environ. 2002, 91, 217–232. [Google Scholar] [CrossRef]
  26. Johnson, J.M.F. Agricultural opportunities to mitigate greenhouse gas emissions. Environ. Pollut. 2007, 150, 107–124. [Google Scholar] [CrossRef] [PubMed]
  27. Tan, Q.-C. Greenhouse Gas Emission in China’s Agriculture: Situation and Challenge. China Popul. Resour. Environ. 2011, 29, 69–75. (In Chinese) [Google Scholar]
  28. Wu, Y.; Feng, K. Spatial-temporal differentiation features and correlation effects of provincial agricultural carbon emissions in China. Environ. Sci. Technol. 2019, 42, 180–190. [Google Scholar]
  29. Fare, R.; Grosskopf, S.; Pasurkajr, C. Environmental Production Functions and Environmental Directional Distance Functions. Energy 2007, 32, 1055–1066. [Google Scholar] [CrossRef]
  30. Oh, D.H. A Global Malmquist-Luenberger Productivity Index. J. Product. Anal. 2010, 34, 183–197. [Google Scholar] [CrossRef]
  31. Tone, K. A Slacks-Based Measure of Efficiency in Data Envelopment Analysis. Eur. J. Oper. Res. 2001, 130, 498–509. [Google Scholar] [CrossRef]
  32. IPCC. Climate Change 2007: Mitigation of Climate Change; Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2008; Volume 45. [Google Scholar] [CrossRef]
  33. Tian, Y.; Zhang, J.B.; He, Y.Y. Research on Spatial-Temporal Characteristics and Driving Factor of Agricultural Carbon Emissions in China. J. Integr. Agric. 2014, 13, 1393–1403. [Google Scholar] [CrossRef] [Green Version]
  34. Min, J.S.; Hu, H. Calculation of greenhouse gases emission from agricultural production in China. China Popul. Resour. Environ. 2012, 22, 21–27. (In Chinese) [Google Scholar]
  35. Xiong, Z.Q.; Xing, G.; Xi, H.T.; Shi, S.L.; Shen, G.Y.; Tu, L.J.; Qian, W. The effects of summer legume crop cultivation on nitrous oxide emissions from upland farmland. Sci. Agric. Sin. 2002, 35, 1104–1108. (In Chinese) [Google Scholar]
  36. Wang, S.B.; Su, W.H. Estimation of nitrous oxide emission and its future change in China. Environ. Sci. 1993, 14, 42–46. (In Chinese) [Google Scholar] [CrossRef]
  37. Liu, L.H.; Jiang, J.Y.; Zong, L.G. Emission inventory of greenhouse gases from agricultural residues combustion: A case study of Jiangsu Province. Environ. Sci. 2011, 32, 1242–1248. [Google Scholar]
Figure 1. Efficiency measure and input slacks.
Figure 1. Efficiency measure and input slacks.
Ijerph 16 03932 g001
Figure 2. The mean gap between green total factor productivity index (GTFPI)and total factor productivity index (TFPI) in different periods in different provinces.
Figure 2. The mean gap between green total factor productivity index (GTFPI)and total factor productivity index (TFPI) in different periods in different provinces.
Ijerph 16 03932 g002
Figure 3. The difference values between the GTFPI and TFPI, GTEC and TEC, and GTPC and TPC in 1998, 2007 and 2016 were distributed in provinces across the country. (a) The gap between the GTFPI and TFPI; (b) the gap between the GTEC and TEC; (c) the gap between the GTPC and TPC.
Figure 3. The difference values between the GTFPI and TFPI, GTEC and TEC, and GTPC and TPC in 1998, 2007 and 2016 were distributed in provinces across the country. (a) The gap between the GTFPI and TFPI; (b) the gap between the GTEC and TEC; (c) the gap between the GTPC and TPC.
Ijerph 16 03932 g003
Figure 4. GTFPI in three periods in each province.
Figure 4. GTFPI in three periods in each province.
Ijerph 16 03932 g004
Figure 5. The spatial distribution of China’s agricultural GTFPI in 1997, 2004, 2010 and 2016.
Figure 5. The spatial distribution of China’s agricultural GTFPI in 1997, 2004, 2010 and 2016.
Ijerph 16 03932 g005
Table 1. Construction of input and output indicators.
Table 1. Construction of input and output indicators.
OutputDesired OutputActual Output Value of Agriculture, Forestry, Husbandry and Fishery (Based on 1997)
Undesired OutputAgricultural Carbon Emissions
InputLaborAgricultural practitioners
LandTotal area sown to crops
MachineryTotal power of agricultural machinery
FertilizerApplication quantity of chemical fertilizer (refractive index)
PesticideConsumption of chemical pesticides
Agricultural filmApplication amount of agricultural film
IrrigationEffective irrigation area
Farm animalsThe number of large cattle at the end of the year
Table 2. CH4 emission coefficients of different rice varieties in China’s provinces (unit: g/m2).
Table 2. CH4 emission coefficients of different rice varieties in China’s provinces (unit: g/m2).
AreaERLRIRAreaERLRIRAreaERLRIR
Beijing0013.23Anhui16.7527.651.24Sichuan6.5518.525.73
Tianjin0011.34Fujian7.7452.643.47Guizhou5.12122.05
Hebei0015.33Jiangxi15.4745.865.42Yunnan2.387.67.25
Shaanxi006.62Shandong0021Tibet006.83
Inner Mongolia008.93Henan0017.85Shaanxi0012.51
Liaoning009.24Hubei17.513958.17Gansu006.83
Jilin005.57Hunan14.7134.156.28Qinghai000
Heilongjiang008.31Guangdong15.0551.657.02Ningxia007.35
Shanghai12.427.553.87Guangxi12.4149.147.78Xinjiang0010.5
Jiangsu16.127.653.55Hainan13.4349.452.29
Zhejiang14.434.557.96Chongqing6.5518.525.73
Note: ER = Early Rice, LR = Late Rice, IR = In-Season Rice.
Table 3. The carbon emission coefficients of major livestock (unit: kg/head/year).
Table 3. The carbon emission coefficients of major livestock (unit: kg/head/year).
ResourcesCH4 Emission CoefficientN2O Emission Coefficient
Enteric Fermentation Manure Management
Cow68161
Buffalo5521.34
Cattle47.811.39
Mule100.91.39
Camel461.921.39
Donkey100.91.39
Horse181.641.39
Live pig13.50.53
Sheep50.150.33
Goat50.170.03
Rabbit0.2540.080.02
Poultry-0.020.02
Table 4. China’s agricultural GTFPI, GTEC, and GTPC in different periods and regions.
Table 4. China’s agricultural GTFPI, GTEC, and GTPC in different periods and regions.
Time1998–20161998–20032004–20092010–2016
WholeGEC0.98980.97890.99720.9923
GTC1.04571.04401.02561.0635
GTFPI1.03481.02191.02281.0553
EastGEC1.00060.99141.01250.9984
GTC1.04161.04131.01451.0656
GTFPI1.04221.03231.02711.0639
CenterGEC0.98190.95810.99540.9912
GTC1.05321.06001.03751.0610
GTFPI1.03411.01561.03271.0516
WestGEC0.98340.98290.98080.9860
GTC1.04251.03291.02851.0632
GTFPI1.02521.01531.00871.0483
Table 5. Global Moran index of the GTFPI from 1998 to 2016.
Table 5. Global Moran index of the GTFPI from 1998 to 2016.
TFP
WqWd
YearMIPVMIPV
1998−0.0970.261−0.0540.243
1999−0.0050.398−0.0670.146
20000.1480.0440.0320.018
2001−0.2570.028−0.0820.077
2002−0.0180.440−0.0130.249
20030.0170.320−0.0210.348
2004−0.1120.203−0.0880.026
20050.0510.183−0.0610.155
2006−0.0560.387−0.0490.254
20070.2300.0120.0850.000
20080.1940.0120.0380.008
20090.0020.364−0.0120.235
20100.0860.1120.0060.085
20110.0180.3320.0290.036
2012−0.0160.440−0.0420.390
20130.0900.1430.0040.132
20140.2160.0110.0150.067
2015−0.1090.162−0.0630.094
2016−0.0400.474−0.0210.353
Note: MI: Moran Index, PV = p-Value, Wq: the first-order adjacent weight matrix, Wd: the geographic distance matrix.

Share and Cite

MDPI and ACS Style

Xu, X.; Huang, X.; Huang, J.; Gao, X.; Chen, L. Spatial-Temporal Characteristics of Agriculture Green Total Factor Productivity in China, 1998–2016: Based on More Sophisticated Calculations of Carbon Emissions. Int. J. Environ. Res. Public Health 2019, 16, 3932. https://doi.org/10.3390/ijerph16203932

AMA Style

Xu X, Huang X, Huang J, Gao X, Chen L. Spatial-Temporal Characteristics of Agriculture Green Total Factor Productivity in China, 1998–2016: Based on More Sophisticated Calculations of Carbon Emissions. International Journal of Environmental Research and Public Health. 2019; 16(20):3932. https://doi.org/10.3390/ijerph16203932

Chicago/Turabian Style

Xu, Xiaocang, Xiuquan Huang, Jun Huang, Xin Gao, and Linhong Chen. 2019. "Spatial-Temporal Characteristics of Agriculture Green Total Factor Productivity in China, 1998–2016: Based on More Sophisticated Calculations of Carbon Emissions" International Journal of Environmental Research and Public Health 16, no. 20: 3932. https://doi.org/10.3390/ijerph16203932

APA Style

Xu, X., Huang, X., Huang, J., Gao, X., & Chen, L. (2019). Spatial-Temporal Characteristics of Agriculture Green Total Factor Productivity in China, 1998–2016: Based on More Sophisticated Calculations of Carbon Emissions. International Journal of Environmental Research and Public Health, 16(20), 3932. https://doi.org/10.3390/ijerph16203932

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