Effects of agricultural technical efficiency on agricultural carbon emission: Based on spatial spillover effect and threshold effect analysis
-
摘要: 农业技术是促进农业产业发展的根本力量, 探究其对农业碳排放的影响机制, 对实现我国“双碳”目标具有重要意义。本文基于2001—2020年我国31个省、直辖市和自治区(港澳台地区以外)的面板数据, 使用随机前沿模型对农业技术效率进行测算, 并对各地区农业碳排放总量与强度进行核算, 构建空间杜宾模型和以农业技术效率为门槛变量的门槛模型, 探究农业技术效率和农业碳排放的空间效应与非线性关系。结果表明: 全国农业碳排放总量与强度近年来呈下降趋势。中部地区农业碳排放总量高于东西部地区, 东部地区农业技术效率高于中西部地区, 而农业碳排放强度则低于中西部地区。农业碳排放强度与农业技术效率具有空间自相关性, 并表现为集聚特征, 集聚类型以高高聚聚和低低聚集为主。农业碳排放强度具有正向空间溢出效应, 而农业技术效率对农业碳排放强度则表现为负向空间溢出, 此外城镇化、人力资本水平和人均耕地面积对农业碳排放强度具有负向影响, 农业经济发展水平和农业受灾程度为正向影响。农业技术效率与农业碳排放强度存在双门槛效应, 当农业技术效率达到“拐点”后, 其对农业碳排放强度的影响转变为负向, 当进一步提升农业技术效率水平后, 其影响力会因边际效应递减而减弱。本研究为探索实现“双碳”目标的路径提供理论基础与政策依据。Abstract: Global warming, caused by the greenhouse effect, has triggered numerous unprecedented extreme weather events globally. Agricultural technology is the fundamental force that promotes the development of the agricultural industry. Studying the impact mechanism of agricultural technology on agricultural carbon emissions will help transform traditional agriculture into ecological, green, and low-carbon modern agriculture, and it is of great significance to the realization of carbon neutrality and carbon peaks. This study used panel data from 31 provinces and cities in China from 2001 to 2020. First, the stochastic frontier model was used to extend existing research from a broad and narrow sense of agricultural technical progress to agricultural technical efficiency. The total agricultural carbon emissions and intensity of agricultural carbon emissions were then calculated and compared. Finally, we constructed the spatial Dubin model and the threshold model with agricultural technical efficiency as the threshold variable, which revealed the spatial effect and non-linear relationship between agricultural technical efficiency and agricultural carbon emissions. The results showed that the total and intensity of agricultural carbon emissions had decreased in recent years. Central China had more agricultural carbon emissions than eastern and western China, and eastern China had a higher technical efficiency of agriculture and a lower carbon emission intensity of agriculture than central and western China. Agricultural carbon emission intensity and technical efficiency had spatial autocorrelation and agglomeration characteristics, and high-high clustering and low-low clustering are the main factors among the provinces. Agricultural carbon emission intensity had a positive spatial spillover effect on itself, but agricultural technical efficiency had a negative spatial spillover effect, which was conducive to the overall reduction of agricultural carbon emissions. Additionally, urbanization, human capital level, and per capita cultivated land area also had negative effects on agricultural carbon emission intensity, but the level of agricultural economic development and the degree of agricultural disaster had positive effects. There was a double threshold effect between agricultural technical efficiency and agricultural carbon emission intensity, which meant that when agricultural technical efficiency reached the “inflection point”, its impact on agricultural carbon emission intensity became negative, and after the level of agricultural technical efficiency was further improved, its influence weakened due to the diminishing marginal effect. Most existing research began with a broad or narrow definition of technological progress, but this study used technical efficiency as the research object after the decomposing technological progress in a broad sense, which further validated the indisputable and decisive role of technological progress in agricultural energy conservation and emission reduction. This study provides a theoretical and policy basis for exploring the path to achieving the “double carbon” goal.
-
随着温室效应的不断加剧, 由温室效应带来的气候变暖在全球范围内引发了前所未有的极端天气事件, 因此世界各国开始谋求合作共同应对随之而来的危机, 中国也提出了要在2030年前二氧化碳排放达到高峰, 力争在2060年实现碳中和的碳达峰碳中和(“双碳”)目标[1]。中国是一个农业大国, 种植业的二氧化碳排放达到全国二氧化碳排放总量的16%~17%[2], 已经成为了不可忽视的重要碳排放源。因此寻找低碳农业的实现路径, 减少农业碳排放对实现我国的“双碳”目标具有重要意义。而在众多影响农业碳排放的因素中, 农业技术已经成了促进节能减排不可忽视的重要力量[3]。2022年中央一号文件指出要研发应用减碳增汇型农业技术, 从而在发展规划层面肯定了农业技术在减少农业碳排放上起到的重要作用。从近代以来人类社会的发展史来看, 技术进步会刺激生产增加能源消耗, 但同时农业技术进步能够提高能源利用效率, 降低农业碳排放强度[4], “双碳”目标约束下的减排政策也会倒逼农业技术的进步。当前中国正处于农业技术快速进步的历史时期, 各地区之间的农业技术水平差距较大, 西部地区的农业技术水平仍然有较大的进步空间。在此背景下, 为促进农业低碳绿色发展, 实现农业发展方式的转变, 就必须发挥农业技术在其中的关键性作用, 同时也是推动农业高质量发展的必然选择。
农业碳排放是指在农业生产过程中, 化肥农药的使用、能源的消耗和土地翻耕灌溉过程中直接或间接造成的温室气体(主要指二氧化碳)的排放[5]。众多学者在农业碳排放相关领域进行了大量研究, 如农业碳排放的测算[6-7]和影响农业碳排放的因素[8-9]。当前有关农业技术与农业碳排放关系的研究, 主要可以分为两个方面: 一是农业技术进步会显著增加农业碳排放总量, 这是由于技术进步往往会促进经济增长, 通过投入农药化肥和能源消耗等途径增加了农业碳排放[10-11]。二是农业技术进步有助于降低农业碳排放强度[4,10,12-13], 这是因为技术进步将改变不同要素之间的边际替代率, 同时提高了能源要素的利用效率, 进而影响农业碳排放强度。
农业全要素生产率(TFP)可被分解为农业技术进步(TC)、农业技术效率(TEC)、农业规模效应(SEC)和要素重置效应(FAEC)[14]。其中, 农业全要素生产率属于广义的农业技术进步; 被分解后的农业技术进步属于狭义的农业技术进步, 是指代表最先进生产技术的生产前沿面整体向外移动, 即采用新技术使得以同样的投入组合生产出更多的产出; 农业技术效率指在现有农业技术水平下, 协调各种农业资源要素, 促使现有的技术水平潜能更大程度地释放; 农业规模效应指农业生产规模的扩大带来的经济效益的提高; 要素重置效应指生产要素在市场供给方之间发生变动和重新配置[15]。现有的研究大多从农业技术进步和规模效应的角度探究其与农业碳排放的关系, 其中农业技术进步的相关研究有: 李成龙等[14]研究认为农业技术进步并不总是有利于减少农业碳排放, 其中机械型技术进步促进了农业碳排放强度的提高, 而生物型技术进步降低了农业碳排放强度。徐清华等[16]研究认为农业机械化降低了本地农业碳排放的同时, 又具有空间负向溢出效应。规模效应的研究有: 李波等[5]通过Kaya恒等式变形对农业碳排放影响因素进行分解研究, 结果表明劳动力规模对碳排放量有一定的抑制作用, 而农业经济发展则会推动农业碳排放。要素重置效应涉及各生产部门、产业组织和生产消费端的要素流动, 其对农业碳排放的影响较为复杂, 目前少有相关研究成果, 而对农业技术效率的研究也仅停留在农业技术效率与农业碳排放存在负脱钩关系[17]的研究层次上, 缺少对其影响机制进行进一步的实证分析。
综上所述, 当前对农业技术进步与农业碳排放的研究, 一方面大多数学者往往局限于广义的农业技术进步与狭义的农业技术进步对农业碳排放影响研究, 甚至在研究中不能有效地区分广义与狭义的技术进步。另一方面缺少农业技术效率和农业碳排放关系的深入研究, 使得这一研究领域长期以来处于空白状态。因此本文将以空间视角研究分析农业技术效率对农业碳排放的影响, 同时构建门槛模型以探究二者的非线性关系, 填补相关研究的空白, 为政府达到农业碳排放减排目标提供理论依据和针对性的政策建议。
1. 理论分析与研究假设
本研究选取狭义农业(种植业)作为研究对象。参考李波等[5]、陈军娟等[17]和何炫蕾等[18]学者的研究, 将农业碳排放定义为在农业生产过程中产生的碳排放, 农业碳排放强度定义为每万元农业总产值引起的农业碳排放。当前对于农业碳排放影响要素研究基本认为农业碳排放具有空间溢出效应[19]。基于此, 农业碳排放强度也应当具有空间溢出效应, 且其效应为正向影响。综上所述提出如下假设:
假设1: 农业碳排放强度具有空间溢出效应, 并表现出正向效应。
农业技术效率指农业现实产出与理论最大可能产量的比率, 其数值大小反映的是生产者利用现有技术的有效程度, 取值范围介于0至1之间, 数值越大表明农业技术效率越高[20]。当前学界对于技术效率的计算大多采用随机前沿模型, 例如杨增旭等[21]利用随机前沿生产函数计算了我国小麦(Triticum aestivum)和玉米(Zea mays)化肥施用的技术效率。汪小勤等[22]运用随机前沿模型验证了农业公共投资能够促进农业技术效率提升。因此本文对比其他生产函数, 选择使用随机前沿模型计算农业技术效率。
现有的研究已经证明了农业技术效率与农业碳排放有高度相关性[17], 而农业技术的传播往往限制较低, 容易形成技术集聚。高鸣等[23]通过实证研究后发现农业技术效率不但不会随着距离的扩大而降低, 甚至会通过规模效应使技术效率获得进一步提升, 证明了农业技术效率存在空间自相关。张思麒等[24]研究认为“前沿型”技术进步具有正向溢出效应, “追随型”技术进步对经济发展水平相对落后的地区产生的贡献更大。基于此, 提出如下假设:
假设2: 农业技术效率具有显著的空间相关性, 对邻近省份的农业碳排放具有负向空间溢出效应。
技术进步对农业碳排放的影响并非一成不变, 一些学者的研究指出不同路径下技术进步对农业碳排放的影响各不相同[13]。同时技术效率的提升意味着各投入要素间的协调得到了优化, 而这种协调的优化并不等同于要素投入的减少, 甚至可能促进要素的投入。此外, 杨义武等[25]认为农业技术进步与技术效率同时具有边际效应, 边际效应的存在也使得农业技术效率对碳排放产生非线性影响。Wang[26]认为门槛模型能够有效解决变量间的非线性关系问题。因此本文采用Hansen[27]提出的面板数据门槛回归模型, 探究农业技术效率与农业碳排放强度存在的非线性相关关系, 相应的假设如下:
假设3: 农业技术效率与农业碳排放强度具有非线性相关关系, 并存在门槛效应。
2. 研究方法
2.1 变量选取与度量
2.1.1 被解释变量: 农业碳排放
本文将农业碳排放设定为农业碳排放强度, 含义为每万元农业GDP所产生的碳排放量, 为避免价格因素带来的影响, 本文以2001年为基准的不变价格种植业总值表示农业GDP。参考学者李波等[5]的研究, 农业碳排放主要来源为农业生产活动中化肥、农药和农膜的生产与使用、农业机械的燃料消耗、农地翻耕和农业灌溉导致的土壤有机碳流失。据此可构建农业碳排放计算公式如下所示:
$$ C = \sum {{C_i} = \sum {{T_i} \times {\delta _i}} } $$ (1) 式中: C为农业碳排放总量; Ci为对应碳源的农业碳排放量; Ti为对应碳排放源的量; δi为对应碳源的碳排放系数, 其具体数值如表1所示。
表 1 各类农业碳排放源的碳排放系数及其来源Table 1. Carbon emission coefficient of each carbon source of agriculture and its reference source碳源
Carbon source计算方法
Computational method碳排放系数
Carbon emission coefficient数据参考来源
Reference source化肥 Chemical fertilizer 实际化肥投入量 Actual fertilizer input 0.859 kg∙kg−1 West, et al[28], ORNL[29] 农药 Pesticide 实际农药投入量 Actual pesticide inputs 4.934 kg∙kg−1 ORNL[29] 柴油 Diesel 农业机械消耗柴油量 Diesel fuel consumption by agricultural machinery 0.593 kg∙kg−1 IREEA 农膜 Agricultural film 农用塑料薄膜使用量 Amount of plastic film used for agriculture 5.18 kg∙kg−1 IPCC[30] 农地翻耕 Agricultural land tilling 农作物实际播种总面积 Actual total area sown of crops 312.6 kg∙hm−2 CAB[31] 农业灌溉 Agricultural irrigation 农作物实际灌溉面积 Actual irrigated area of crops 20.476 kg∙hm−2 RDRCH[5] ORNL: 美国橡树岭国家实验室; IREEA: 南京农业大学农业资源与生态环境研究所; IPCC: 政府间气候变化专门委员会; CAB: 中国农业大学农学与生物技术学院; RDRCH: 湖北农村发展研究中心。ORNL: Oak Ridge National Laboratory, USA; IREEA: Institute of Resource, Ecosystem and Environment of Agriculture, Nanjing Agricultural University; IPCC: Intergovernmental Panel on Climate Change; CAB: College of Agronomy and Biotechnology, China Agricultural University; RDRCH: Rural Development Research Center of Hubei. 因各地区农业规模存在较大差异性, 仅以农业碳排放总量衡量地区农业碳排放水平缺少说服力, 而农业碳排放与农业经济发展水平密不可分, 因此本文选择将农业碳排放强度作为被解释变量。其计算公式为:
$$ {\rm{CI}}_{it}=C_{it}/{\rm{GCP}}_{it} $$ (2) 式中: i为省份, t为时间, C
it为农业碳排放总量, CIit为农业碳排放强度, GCPit为不变价格的种植业总产值。 2.1.2 解释变量: 农业技术效率
与要求技术中性、规模报酬不变的传统柯布-道格拉斯(C-D)生产函数和不变替代弹性生产函数(CES)不同, 超越对数生产函数允许各投入要素之间的替代弹性是可变的, 并能反映出投入要素的交叉效应对产出变量的影响[32]。因此本文通过构建超越对数生产函数的随机前沿模型计算各地区历年农业技术效率。借鉴白俊红等[33]和李成龙等[14]的模型构建思路, 构建具体的模型如下:
$$ \begin{split} &\qquad \ln {Y_{it}} = {\alpha _0} + \sum\limits_n {{\alpha _n}\ln {x_{nit}}} + \frac{1}{2}\sum\limits_n {\sum\limits_m {{\alpha _{nm}}\ln {x_{nit}}\ln {x_{mit}} + } } \\& \sum\limits_n {{\alpha _{nt}}T\ln {x_{nit}} + } {\alpha _t}T + \frac{1}{2}{\alpha _{tt}}{T^2} + {\nu _{it}} - {\mu _{it}} \\[-18pt] \end{split}$$ (3) 式中: i为省份; t为时间; Yit为农业总产出, 本文以不变价格的种植业总产值表示(亿元); xnit、xmit分别为第n个投入量和第m个投入量, 本文选取的投入要素为第一产业从业人员数(万人)、农作物播种总面积(hm2)、农业机械总动力(万kW)、化肥施用量(万t); T表示技术效率变化的时间趋势; α表示各参数系数; νit表示随机误差项; µit代表技术无效因素, 二者互相独立并服从正态分布。基于此, 农业技术效率的计算公式为:
$$ \mathrm{ATE}_{i t}=\exp \left(-\mu_{i t}\right) $$ (4) 式中: ATEit为农业技术效率, 取值范围为0<ATEit≤1, 数值越大表明农业技术效率越高。
2.1.3 控制变量
参照现有文献[5,8-9,11,14,16,19], 总结归纳影响农业碳排放的影响因素, 将农业经济发展水平、人力资本水平、农业受灾程度和人均耕地面积等4个变量设置为控制变量。为了更贴近农村劳动力现实情况, 本文选取各省份农村地区初中及其以上学历人数与农村地区6岁以上总人数的比率作为衡量人力资本水平的指标。
2.2 数据说明及描述性统计
考虑到数据的可获得性与实际情况, 本文以港澳台地区以外, 中国31个省、直辖市和自治区2001—2020年共20年的面板数据作为分析样本, 相关统计数据来源于《中国统计年鉴》《中国农村统计年鉴》《中国人口统计年鉴》《中国人口与就业统计年鉴》和各省(直辖市、自治区)统计年鉴, 在不影响结果的前提下对极少部分缺失数据使用插值法补齐, 总结归纳各变量设定及描述性统计情况如表2所示。
表 2 农业碳排放研究相关变量说明及描述性统计Table 2. Description and descriptive statistics of related variables of agricultural carbon emission research变量类型
Type of
variable指标名称
Name of
index变量符号
Variable
symbol测算方法
Calculating
method单位
Unit均值
Mean
value标准差
Standard
deviation被解释变量
Explained variable农业碳排放强度
Carbon emission intensity
of agricultureCarbon I 农业碳排放/不变价格的种植业总产值
Agricultural carbon emissions /
gross plantation value at constant pricet·(104 ¥)−1 0.51 0.28 解释变量
Explaining variable农业技术效率
Agricultural technical
efficiencyATE 基于超越对数生产模型计算
Based on translog production model— 0.19 0.10 控制变量
Control variable城镇化率
Urbanization rateUrban 城镇常住人口/常住总人口
Permanent urban population /
total permanent population% 0.50 0.17 农业经济发展水平
Level of agricultural
economic developmentEy 农业总产值/农业从业人员
Total agricultural output value /
agricultural employees×104 ¥·capita−1 1.79 1.30 人力资本水平
Level of human capitalEdu 农村地区初中及其以上学历人数/
农村地区6岁以上人数
Number of people with junior high school
education or above in rural areas / number
of people over 6 years old in rural areas% 0.05 0.01 农业受灾程度
Degree of agricultural
damageDisaster 受灾面积/农作物播种面积
Affected area / crop sown area% 0.22 0.16 人均耕地面积
Arable land per capitaArea 农作物播种面积/第一产业从业人数
Crop sown area / number of
workers in primary industryhm2·capita−1 0.63 0.34 样本时间
Sample timeYear 年份(2001—2020年)
Year (2001—2020)— — — 研究对象
Object of studyProvince 省份(31个省、直辖市、自治区)
Provinces (31 provinces, cities, autonomous regions)
— — — 2.3 模型建立
2.3.1 空间自相关检验
对变量的空间相关性分析, 首先必须明确变量是否存在空间效应, 而检验空间效应最常用的方法是计算Moran’s I指数。因此构建全局Moran’s I指数模型, 对农业碳排放强度和农业技术效率的空间自相关性进行检验。公式如下:
$$ \begin{split} I = \frac{{n\displaystyle \sum\limits_{i = 1}^n {\displaystyle \sum\limits_{j = 1}^n {{w_{ij}}\left( {{x_i} - \overline x } \right)} } \;\left( {{x_j} - \overline x } \right)}}{{\displaystyle \sum\limits_{i = 1}^n {\displaystyle \sum\limits_{j = 1}^n {{w_{ij}}{{\displaystyle \sum\limits_{j = 1}^n {\left( {{x_i} - \overline x } \right)}^2 }}} } }} \end{split} $$ (5) 式中: I为全局Moran’s I指数; n为选取省份的个数, 本文n=31; xi、xj为i省份和j省份的农业碳排放强度和农业技术效率;
$\overline x $ 为样本平均值; wij为邻接空间权重矩阵, 当i与j不相邻时设置为0, 反之则为1。全局Moran’s I指数服从正态分布, 其取值范围为[−1, 1], 当Moran’s I指数为0时, 说明变量间不存在空间相关性; 当Moran’s I指数大于0时, 说明变量间存在空间正相关性; 当Moran’s I指数小于0时, 说明变量间存在空间负相关性; Moran’s I指数的绝对值越大, 说明其相关性越强。2.3.2 空间面板回归模型
$$ \left\{ \begin{gathered} {y_{it}} = \tau {\kern 1pt} {y_{i,t - 1}} + \rho {\kern 1pt} {{w'}_i}{y_t} + {{x'}_{it}}\beta + {{w'}_i}{X_t}\delta + {u_i} + {\gamma _t} + {\varepsilon _{it}} \\ {\varepsilon _{it}} = \lambda m_i^t{\varepsilon _t} + {\nu _{it}} \\ \end{gathered} \right. $$ (6) 式中: τ、ρ、β均为系数; y
i,t−1为被解释变量的一阶滞后; ${{w'}_i}{X_t}\delta $ 为解释变量的空间滞后;${{w'}_i}$ 为空间权重矩阵W的第i行; γt为时间效应;$m_i^t $ 为扰动项空间权重矩阵M的第i行; ui为区域i的个体效应; εit为残差; εt为空间滞后误差项; vit为随机误差项。为考虑不同情况下农业碳排放空间溢出效应的差异性, 本文构建3种空间权重矩阵, 其具体含义如下所示:
$$ 邻接空间矩阵: {w}_{1}= \left\{ \begin{array}{l} 1\text{,} \;\;i 与 {j} 相邻 \\ 0\text{,}\;\;i与 {j} 不相邻 \end{array} \right. $$ (7) $$ \begin{split} \\ 地理距离矩阵: {w}_{2}= \left\{ \begin{array}{l} \dfrac{1}{{d}^{2}}\text{,} \qquad i\ne j\\ 0\text{,} \qquad \;\;\; i=j \end{array} \right. \end{split} $$ (8) 式中: d为地区i与j的距离, 本文将各省省会之间的距离作为衡量距离的指标。
$$ 经济距离矩阵: {w}_{3}= \left\{ \begin{array}{l} \dfrac{1}{\left|\overline{{y}_{i}}-\overline{{y}_{j}}\right|}\text{ ,}\;\; i\ne j\\ 0\text{ ,}\;\; i=j \end{array} \right. $$ (9) 式中:
$\overline{y_i} $ 表示地区i的人均GDP, 考虑到本文的研究数据为2001到2020年的面板数据, 因此以20年人均GDP的均值作为衡量指标。空间面板回归模型的具体形式一般根据LM检验、LR检验、Wald检验和Hausman检验的结果来确定。其中, Hausman检验用于确定模型使用固定效应或随机效应, 本文Hausman检验的原假设为“支持模型使用随机效应”; LM检验用于确定模型是否存在空间自相关性, 其原假设为“模型不存在空间自相关”; LR检验与Wald检验用于确定模型是否使用空间误差模型(SEM)或空间滞后模型(SLM), 其原假设为“空间杜宾模型(SDM)不能简化为空间误差模型或空间滞后模型”。
根据以上检验结果, 对公式(6)进行改进补充。如果模型为空间杜宾模型, 则λ=0; 若模型为空间滞后模型, 则λ=0且δ=0; 若模型为空间误差模型, 则τ=ρ=0且δ=0。
2.3.3 门槛回归模型
效率的提升并不代表投入要素的减少, 甚至可能促进要素的投入, 此外技术效率具有边际效应, 因而农业技术效率对农业碳排放的影响可能是非线性的。综上所述, 本文参考Hansen[27]的面板回归模型, 依照本文理论分析中的研究假设, 将农业碳排放强度(Carbon I)作为被解释变量, 农业技术效率(ATE)同时作为门槛变量和核心解释变量, 构建门槛回归模型, 公式如(10)所示:
$$ \begin{split} &\qquad {\rm{Carbon}}\,\,{\rm{I}} = {\omega _i} + {\alpha _{\text{1}}}{\rm{Contro}}{{\rm{l}}_{it}} + {\beta _1}{\rm{ATE}} \cdot D\left( {{\rm{ATE}} < {\theta _1}} \right) + \\ & \; {\beta _2}{\rm{ATE}} \cdot D( {{\theta _1} \leqslant{\rm{ ATE}} < {\theta _2}} ) + \cdots + {\beta _n}{\rm{ATE}} \cdot D( {\theta _{n - 1}} \leqslant {\rm{ATE}} < \\ & {\theta _n} ) + {\beta _{n + 1}}{\rm{ATE}} \cdot D\left( {{\rm{ATE}} > {\theta _n}} \right) + {\mu _{it}} \\[-10pt] \end{split}$$ (10) 式中: θ1, θ2,
$\cdots ,$ θn表示n个不同水平的门槛值; Control表示控制变量; D(·)为指标函数; μit为随机扰动项。3. 结果与分析
3.1 农业碳排放变化规律
依据公式(1)、(2), 按照国家统计局关于东中西部的划分, 将北京、天津、河北、辽宁、上海、江苏、浙江、福建、山东、广东和海南划分为东部地区, 山西、吉林、黑龙江、安徽、江西、河南、湖北和湖南为中部地区, 内蒙古、广西、重庆、四川、贵州、云南、西藏、陕西、甘肃、青海、宁夏和新疆为西部地区, 计算2001—2020年中国不同地区农业碳排放总量和农业碳排放强度, 所得农业碳排放总量和农业碳排放强度分区域趋势变化图如图1所示。
![]()
图 1 2001—2020年中国不同地区农业碳排放总量(a)和农业碳排放强度(b)Figure 1. Total amount (a) and intensity (b) of agricultural carbon emissions in different regions of China from 2001 to 2020从图1可以看出, 2001—2020年各地区农业碳排放总量表现为上升趋势, 但自2015年以后有所降低, 而农业碳排放强度则一直保持下降趋势。分区域来看, 中部地区的农业碳排放总量高于东部与西部, 这是因为中部地区包括河南、黑龙江、吉林在内均为农业大省, 拉高了农业碳排放总量。农业碳排放强度则是东部地区低于中西部地区, 这主要是因为东部地区农业技术水平、城镇化率等显著降低农业碳排放强度的指标优于中西部地区。以2020年为例, 农业碳排放总量(万t)最高的10个省(直辖市、自治区)为河南(1237.51)、山东(961.28)、黑龙江(828.12)、安徽(683.90)、河北(681.71)、江苏(648.55)、新疆(624.42)、四川(613.10)、湖北(608.25)和湖南(589.34)。农业碳排放强度(t∙万元−1)最高的10个省(直辖市、自治区)为内蒙古(0.48)、吉林(0.44)、安徽(0.34)、新疆(0.31)、宁夏(0.31)、山西(0.30)、甘肃(0.27)、上海(0.25)、青海(0.25)和黑龙江(0.25)。
结合以上数据可知, 以农作物播种总面积衡量的农业大省所产生的农业碳排放总量相对较高, 中部地区的碳排放总量较高, 西部地区的农业碳排放强度较高。中国地域广大, 幅员辽阔, 各地区的农业生产要素禀赋各不相同, 是造成这些差异的主要原因。因此各地区应当依据自己的农业发展情况, 实事求是地制定合理的碳减排目标。
3.2 农业技术效率变化规律
构建超越对数生产函数, 依据公式(3)和(4)计算历年各区域农业技术效率, 结果如图2所示。由图2结果可知, 2001—2020年各地区农业技术效率水平不断提高, 东部地区农业技术效率显著高于中西部地区, 而中部地区仅略微高于西部地区, 且二者均低于全国平均水平。分省市来看, 20年平均农业技术效率高于全国平均水平(0.186)的省市从大到小依次为: 北京(0.381)、上海(0.286)、浙江(0.280)、天津(0.238)、福建(0.231)、广东(0.224)、辽宁(0.223)、海南(0.223)、山东(0.211)、河北(0.206)和江苏(0.196); 低于全国平均水平从小到大依次为: 内蒙古(0.124)、安徽(0.129)、宁夏(0.137)、青海(0.137)、贵州(0.139)、山西(0.141)、云南(0.141)、江西(0.152)、广西(0.155)、湖北(0.155)、黑龙江(0.156)、陕西(0.157)、重庆(0.157)、甘肃(0.162)、河南(0.164)、四川(0.168)、吉林(0.169)、湖南(0.175)、西藏(0.78)和新疆(0.180)。
![]()
图 2 2001—2020年中国不同地区农业技术效率Figure 2. Agricultural technical efficiencies in different regions of China from 2001 to 2020对比分析各地区农业技术效率差异可知, 东部沿海省市的农业技术水平远高于内陆省份。与此同时, 随着农村劳动力转移进城, 中西部地区农业生产水平却依旧维持在较低水平, 而未来中国的农业发展方式也要从过去的粗放式低效益转变为产业化、生态化、国际化和地区化。在这样的背景下, 中西部地区要积极学习借鉴东部沿海省市的发展经验, 引进新型农业机械和信息化管理技术, 才能从根本上解决当前农业发展中存在的问题。考虑到北京、上海、天津第一产业产值分别位列全国各省(直辖市、自治区)倒数第2、第1和第4, 浙江省在提高农业技术效率上的经验更具有参考价值, 值得其他各省市学习。
3.3 农业碳排放与农业技术效率的空间自相关性检验结果分析
基于空间面板数据, 利用Stata 16软件参照公式(5)计算2001—2020年农业碳排放强度与农业技术效率的全局Moran’s I指数和2020年农业碳排放强度与农业技术效率的局部Moran’s I指数, 结果分别如表3和图3所示。
表 3 2001—2020年农业碳排放强度和农业技术效率的全局Moran’s I指数及检验Table 3. Moran’s I indexes and test of agricultural carbon emission intensity and agricultural technical efficiency from 2001 to 2020年份
Year农业碳排放强度
Agricultural carbon
emission intensity年份
Year农业技术效率
Agricultural technical
efficiencyMoran’s I Z P Moran’s I Z P 2001 0.278 2.636 0.004 2001 0.338 3.562 0.000 2002 0.224 2.196 0.014 2002 0.345 3.592 0.000 2003 0.237 2.287 0.011 2003 0.351 3.619 0.000 2004 0.220 2.151 0.016 2004 0.356 3.643 0.000 2005 0.185 1.840 0.033 2005 0.362 3.665 0.000 2006 0.185 1.866 0.031 2006 0.366 3.684 0.000 2007 0.221 2.131 0.017 2007 0.370 3.701 0.000 2008 0.118 1.267 0.103 2008 0.374 3.717 0.000 2009 0.198 1.971 0.024 2009 0.378 3.73 0.000 2010 0.208 2.056 0.020 2010 0.381 3.742 0.000 2011 0.195 1.931 0.027 2011 0.384 3.752 0.000 2012 0.211 2.059 0.020 2012 0.386 3.761 0.000 2013 0.203 1.999 0.023 2013 0.389 3.769 0.000 2014 0.177 1.826 0.034 2014 0.391 3.776 0.000 2015 0.229 2.258 0.012 2015 0.393 3.782 0.000 2016 0.249 2.485 0.006 2016 0.394 3.787 0.000 2017 0.277 2.644 0.004 2017 0.396 3.792 0.000 2018 0.331 3.189 0.001 2018 0.397 3.795 0.000 2019 0.320 3.103 0.001 2019 0.399 3.798 0.000 2020 0.368 3.519 0.000 2020 0.400 3.801 0.000 ![]()
图 3 2020年农业碳排放强度(a)与农业技术效率(b) Moran’s I散点图1: 北京; 2: 天津; 3: 河北; 4: 山西; 5: 内蒙古; 6: 辽宁; 7: 吉林; 8: 黑龙江; 9: 上海; 10: 江苏; 11: 浙江; 12: 安徽; 13: 福建; 14: 江西; 15: 山东; 16: 河南; 17: 湖北; 18: 湖南; 19: 广东; 20: 广西; 21: 海南; 22: 重庆; 23: 四川; 24: 贵州; 25: 云南; 26: 西藏; 27: 陕西; 28: 甘肃; 29: 青海; 30: 宁夏; 31: 新疆。Figure 3. Moran’s I scatter diagrams of agricultural carbon emission intensity and agricultural technical efficiency in 20201: Beijing; 2: Tianjin; 3: Hebei; 4: Shanxi; 5: Inner Mongolia; 6: Liaoning; 7: Jilin; 8: Heilongjiang; 9: Shanghai; 10: Jiangsu; 11: Zhejiang; 12: Anhui; 13: Fujian; 14: Jiangxi; 15: Shandong; 16: Henan; 17: Hubei; 18: Hunan; 19: Guangdong; 20: Guangxi; 21: Hainan; 22: Chongqing; 23: Sichuan; 24: Guizhou; 25: Yunnan; 26: Tibet; 27: Shaanxi; 28: Gansu; 29: Qinghai; 30: Ningxia; 31: Xinjiang.由表3可知, 2001—2020年农业碳排放强度与农业技术效率的全局Moran’s I指数均为正, 且通过P<0.05显著性检验, 说明农业碳排放强度与农业技术效率存在显著空间相关性并表现出集聚的特征。从时间变化来看, 农业碳排放强度的Moran’s I值随着时间推移表现为上下波动后上升, 说明农业碳排放强度的聚集程度呈增强趋势; 农业技术效率的Moran’s I值随着时间推移表现持续升高趋势, 说明农业技术效率的聚集程度在逐渐增强。
Moran’s I散点图可形象表现地区间的空间相关性, 其含义为第1象限内的省份农业碳排放强度或农业技术效率较高, 形成了高高集聚; 第2象限内的省份自身水平较低而周边水平较高, 形成了低高聚集; 第3象限内的省份自身水平较高而周边水平较低, 形成了高低聚集; 第4象限内的省份自身水平较低的同时周边水平同样较低, 形成了低低聚集。从结果来看, 农业碳排放强度形成高高聚集的省份共有10个, 低高聚集的有5个, 高低聚集的有4个, 低低聚集的有12个, 主要以高高聚集和低低聚集为主; 农业技术效率形成高高聚集的有9个, 低高聚集的有2个, 高低聚集的有2个, 低低聚集的有18个, 同样以高高聚集和低低聚集为主。
上述结果表明, 农业碳排放强度与农业技术效率显著存在空间相关性, 聚集类型均以高高聚集和低低聚集为主。根据上文空间面板数据模型构建的前提条件, 本文将进一步对农业技术效率对农业碳排放的空间溢出效应进行检验分析。
3.4 农业技术效率对农业碳排放强度的空间面板回归模型选择
根据空间自相关检验的结果, 农业碳排放强度与农业技术效率的全局Moran’s I指数显著且不等于0, 说明变量存在空间相关性, 可以对其进行进一步分析, 确定其所使用的空间面板回归模型。分别对3种空间权重矩阵的空间面板回归模型进行LM检验、LR检验、Wald检验和Hausman检验, 所得结果如表4所示。
表 4 农业技术效率对农业碳排放强度空间面板回归模型的LM检验、LR检验、Wald检验和Hausman检验结果Table 4. Results of LM test, LR test, Wald test and Hausman test of the spatial panel regression model of agricultural technical efficiency to agricultural carbon emission intensity检验类别
Test category检验项目
Inspection itemW1 W2 W3 LM检验 LM test LM(误差)检验 LM (error) test 431.232*** 426.437*** 64.628*** 稳健LM(误差)检验 Robust LM (error) test 300.825*** 370.055*** 50.390*** LM(滞后)检验 LM (lag) test 134.108*** 77.351*** 15.260*** 稳健LM(滞后)检验 Robust LM (lag) test 3.701* 20.969*** 1.023 LR检验 LR test SDM与SLM的卡方检验 SDM&SLM chi2 88.50*** 68.00*** 46.97*** SDM与SEM的卡方检验 SDM&SEM chi2 64.73*** 44.35*** 49.50*** Wald检验 Wald test SDM与SLM的卡方检验 SDM&SLM chi2 27.73*** 36.70*** 31.27*** SDM与SEM的卡方检验 SDM&SEM chi2 24.76*** 21.59*** 18.58*** Hausman检验 Hausman test 显著性卡方检验 Prob>=chi2 54.41*** 73.37*** 124.75*** 拟合度 R2 时间效应 Time 0.0750 0.4271 0.0271 个体效应 Ind 0.8809 0.9098 0.8895 混合效应 Both 0.7442 0.7408 0.7247 模型选择 Model select 个体固定效应
的SDM模型
SDM model of individual
fixed effects个体固定效应
的SDM模型
SDM model of individual
fixed effects个体固定效应
的SDM模型
SDM model of individual
fixed effects***和*分别表示在1%和10%水平显著; W1、W2和W3分别代表邻接空间矩阵、地理距离矩阵和经济距离矩阵; LM检验、LR检验和Wald检验分别代表拉格朗日乘子检验、似然比检验和沃尔德检验; SLM、SEM和SDM分别代表空间滞后模型、空间误差模型和空间杜宾模型; Chi2为卡方检验。*** and * indicate significance at 1% and 10% levels, respectively. W1, W2 and W3 represent adjacency spatial matrix, geographical distance matrix and economic distance matrix, respectively. LM test, LR test and Wald test represent Lagrange multiplier test, likelihood ratio test and Wald test, respectively. SLM, SEM and SDM represent spatial lag model, spatial error model and spatial Durbin model, respectively. Chi2 is the chi-square test. 由表4可知, 邻接空间矩阵、地理距离矩阵和经济距离矩阵的空间面板模型均通过LM检验拒绝原假设, 表明模型存在空间交互效应, 且空间误差模型显著水平优于空间滞后模型; LR检验与Wald检验结果均表明拒绝原假设, 应将模型设置为空间杜宾模型; Hausman检验与拟合度分析结果表明应选择个体固定效应模型。
综上所述, 根据LM检验、LR检验和Hausman检验结果, 确定3种空间权重矩阵的空间面板模型均为个体固定效应的空间杜宾模型, 构建模型公式如下所示:
$$ \begin{split} \\ & \qquad {\rm{Carbon}}\;{{\rm{I}}_{it}} = \rho {W_j}{\rm{Carbon}}\;{{\rm{I}}_{it}} \; + \, {\beta _1}{\rm{AT}}{{\rm{E}}_{it}} \,+ \, {\beta _2}{\rm{Urba}}{n_{it}} + \\ & { \; \beta _3}{\rm{E}}{{\rm{y}}_{it}} \, + \, {\beta _4}{\rm{Ed}}{{\rm{u}}_{it}} \, + \,{\beta _5}{\rm{Disaster}}{_{it}} \, + {\beta _6}{\rm{Are}}{{\rm{a}}_{it}} + \;{\theta _1}{W_j}{\rm{AT}}{{\rm{E}}_{it}} + \\ & {\theta _2}{W_j}{\rm{Urba}}{{\rm{n}}_{it}} \; + \;{\theta _3}{W_j}{\rm{E}}{{\rm{y}}_{it}} \;+ \, {\theta _4}{W_j}{\rm{Ed}}{{\rm{u}}_{it}} \;+ \; {\theta _5}{W_j}{\rm{Disaster}}{_{it}} + \\ & {\theta _6}{W_j}{\rm{Are}}{{\rm{a}}_{it}} + {\mu _i} + {\varepsilon _{it}} \\[-10pt] \end{split} $$ (11) 式中: i为省份; t为时间; ρ为空间回归相关系数, 表示相邻地区的被解释变量对本地区被解释变量的影响; β1~β6表示本地区解释变量对本地区被解释变量的影响程度; θ1~θ6表示相邻地区解释变量对本地区被解释变量的影响程度; Wj为空间权重矩阵, j=1,2,3, 分别表示邻接空间矩阵、地理距离矩阵和经济距离矩阵; Acarbon为被解释变量; ATE、Urban、Ey、Edu、Disaster、Area为解释变量, 具体意义见表2; μi为i地区的个体效应; εit为随机误差项。
3.5 农业技术效率对农业碳排放强度的空间杜宾模型回归结果分析
将3种矩阵的模型均设置为个体固定效应的空间杜宾模型, 使用Stata 16软件对其分别进行回归分析, 所得结果分别如表5、表6和表7所示。
表 5 农业技术效率对农业碳排放强度的邻接空间矩阵空间杜宾模型回归结果及效应分解Table 5. Regression results and effect decomposition of the adjacency space matrix spatial Durbin model of agricultural technical efficiency to agricultural carbon emission intensity变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 6.203***(0.452) −7.081***(0.473) 5.529***(0.431) −7.747***(0521) −2.217***(0.393) Urban −0.361*** (0.107) 0.095(0.128) −0.390***(0.095) −0.282(0.188) −0.672***(0.176) Ey −0.036*** (0.008) 0.090***(0.012) −0.021**(0.008) 0.160***(0.025) 0.139***(0.029) Edu 0.473(1.138) −4.457***(1.529) −0.416(1.077) −9.736***(2.763) −10.151***(3.004) Disaster 0.079***(0.030) 0.049(0.050) 0.099***(0.030) 0.222**(1.103) 0.321***(0.113) Area 0.007***(0.002) −0.016***(0.004) 0.004*(0.002) −0.029***(0.008) −0.024***(0.009) rho 0.602***(0.034) sigma2_e 0.005***(0.0003) ***、**和*分别表示在1%、5%和10%水平显著; Main意为本地区解释变量对本地区被解释变量的影响程度系数β, Wx意为相邻地区解释变量对本地区被解释变量的影响程度系数θ, Spatial意为相邻地区的被解释变量对本地区被解释变量的影响系数ρ, Variance意为差异系数; 直接效应、间接效应与总效应分别表示利用偏微分法的无偏估计结果, 即本地区解释变量对本地区被解释变量的影响程度, 相邻地区解释变量对本地区被解释变量的影响程度和直接效应与间接效应的总和; ATE、Urban、Ey、Edu、Disaster、Area、rho和sigma2_e分别代表农业技术效率、城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积、空间自相关系数ρ和个体效应的特异误差; 括号内为标准差。***, ** and * indicate significance at 1%, 5% and 10% levels, respectively. Main means the coefficient β of the influence degree of the explanatory variables in the local region on the explained variables in the local region; Wx means the coefficient θ of the influence degree of the explanatory variables in neighboring regions on the explained variables in the local region; Spatial means the influence coefficient ρ of the explained variables in neighboring regions on the explained variables in the local region; and Variance means the difference coefficient. Direct effect, indirect effect and total effect respectively represent the unbiased estimation results using partial differential method, that is, the degree of influence of local explanatory variables on local explained variables, the degree of influence of neighboring explanatory variables on local explained variables and the sum of direct and indirect effects. ATE, Urban, Ey, Edu, Disaster, Area, rho and sigma2_e represent agricultural technical efficiency, urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capital cultivated land area, spatial autocorrelation coefficient ρ and individual effect specific error, respectively. Standard deviations are in parentheses. 表 6 农业技术效率对农业碳排放强度的地理距离矩阵空间杜宾模型回归结果及效应分解Table 6. Geographical distance matrix spatial Durbin model regression results and effect decomposition of agricultural technical efficiency to agricultural carbon emission intensity变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 7.110***(0.390) −7.914***(0.418) 6.672***(0.386) −9.018***(0.530) −2.346***(0.468) Urban −0.363***(0.096) 0.185(0.120) −0.374***(0.088) −0.143(0.199) −0.517***(0.185) Ey −0.040***(0.007) 0.075***(0.015) −0.032***(0.008) 0.139***(0.042) 0.107**(0.045) Edu 0.076(1.040) −3.762**(1.486) −0.468(0.995) −10.322***(3.365) −10.790***(3.579) Disaster 0.072***(0.027) 0.001(0.060) 0.079***(0.026) 0.133(0.154) 0.213(0.159) Area 0.009***(0.002) −0.005(0.056) 0.010***(0.023) 0.002(0.015) 0.011(0.017) rho 0.654***(0.037) sigma2_e 0.004***(0.0002) 同表5注释。Note the same as in Table 5. 表 7 农业技术效率对农业碳排放强度的经济距离矩阵空间杜宾模型回归结果及效应分解Table 7. Economic distance matrix spatial Durbin model regression results and effect decomposition of agricultural technical efficiency on agricultural carbon emission intensity变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 6.393***(0.379) −7.741***(0.389) 5.837***(0.370) −8.755***(0.446) −2.918***(0.384) Urban −0.444***(0.093) 0.252**(0.126) −0.445***(0.085) 0.025(0.180) −0.420**(0.166) Ey −0.033***(0.008) 0.126***(0.016) −0.018**(0.008) 0.222***(0.031) 0.204***(0.033) Edu 3.181***(1.098) −10.768***(1.732) 1.943*(1.002) −18.431***(2.876) −16.488***(2.945) Disaster 0.102***(0.027) 0.028(0.054) 0.114***(0.028) 0.169(0.112) 0.282**(0.128) Area 0.010***(0.002) −0.018***(0.004) 0.008***(0.002) −0.027***(0.009) −0.019*(0.010) rho 0.537***(0.035) sigma2_e 0.005***(0.0003) 同表5注释。Note the same in Table 5. 表5、表6和表7的结果表明, 3种空间权重矩阵在个体固定效应空间杜宾模型下, 农业碳排放强度空间自回归项rho的系数为正, 且均通过了1%显著性检验, 说明农业碳排放强度具有正向的空间溢出效应, 本省农业碳排放强度的提高会促进邻省的农业碳排放强度的上升, 据此假设1得到验证。此外, 根据LeSage等[34]的研究, 若被解释变量的空间自相关rho系数显著且不为0, 则不能直接使用空间矩阵估计系数Wx的回归系数衡量解释变量对被解释变量的空间溢出效应, 应当使用偏微分法得到解释变量的无偏估计结果, 即直接效应、间接效应和总效应。
1)农业技术效率。农业技术效率在3个空间权重矩阵模型中, 直接效应均表现为显著正向影响, 间接效应和总效应则表现为显著负向影响。农业技术效率的提高会刺激农业要素的投入, 此时每一单位的要素投入将会带来更多的农业产出, 表现为边际效应递增。而农业碳排放也将随着农业产量的提升而增大, 因此农业技术效率的提升将促进农业碳排放。而间接效应为负, 则证明了农业技术具有“示范效应”, 邻近省份学习先进的农业技术, 可以有效降低该省的农业碳排放强度, 并且从总效应结果上看, 农业技术效率的提高归根结底有利于抑制农业碳排放, 据此假设2得到验证。
2)城镇化。城镇化在3个空间权重矩阵模型中, 直接效应和总效应表现为显著负向影响, 间接效应未表现出显著性。就业城镇化的含义是农村地区的人口转移至城镇就业从而实现的城镇化, 就业人口的城镇化促进了农业要素投入的集约化, 从而降低了农业碳排放强度。因此, 提高就业城镇化率, 能够降低农业碳排放强度。
3)农业经济发展水平。农业经济发展水平在3个空间权重矩阵模型中, 直接效应表现为显著负向影响, 而间接效应和总效应则表现为显著正向影响。农业经济发展水平的提升意味着要素投入减少的同时反而有更多的产出, 这意味着同等产出水平下所投入的要素量减少了, 这使得农业碳排放也降低了。然而本地区农业经济发展水平的提高并不能降低相邻地区农业碳排放强度, 这可能是因为存在“挤出”效应, 农业经济发展水平较高的地区占据了较多的“低碳资源”, 使得邻近地区在降低农业碳排放强度上不具有相对优势。总体上看, 农业经济发展水平对农业碳排放强度有促进作用。
4)人力资本水平。人力资本水平在3个空间权重矩阵模型中, 间接效应和总效应表现为显著负向影响, 直接效应则未表现出一致影响方向。总的来看, 人力资本水平的提升意味着劳动力素质的提升, 农业劳动者将更善于利用先进的农业生产技术与管理手段, 促使农业生产在更低的要素投入下拥有更高的产量, 节约了能源消耗, 促进了农业碳排放强度降低。此外, 邻近地区的农业从业人员也通过学习先进技术, 使人力资本水平对邻省农业碳排放强度具有负向的空间溢出效应。因此, 人力资本水平越高, 农业碳排放强度越低。
5)农业受灾程度。农业受灾程度在3个空间权重矩阵模型中, 直接效应、间接效应和总效应均表现为正向影响。农业耕作受到破坏后, 将在短时间内投入大量要素对其进行重建, 从而大幅度提高碳排放, 因此其对农业碳排放强度的影响为正, 而农业灾害具有规模范围性, 往往连片出现, 因此表现出正向空间溢出效应。总体而言, 农业受灾程度越深, 农业碳排放强度越高。
6)人均耕地面积。人均耕地面积在3个空间权重矩阵模型中, 直接效应均表现为显著正向, 地理距离矩阵的间接效应与总效应表现为正向但未显著, 邻接空间矩阵与经济距离矩阵的间接效应与总效应表现为显著负向。人均耕地面积的提升意味着农业耕地面积、能源消耗和农资消耗的增加, 从而促进了本地区农业碳排放强度的增加, 而本地区人均耕地面积的增加使得相邻地区的农业要素被大量转移, 使得相邻地区的农业碳排放有所降低。从整体来看, 人均耕地面积的增加能够降低农业碳排放强度。
3.6 农业碳排放研究中变量的多重共线性检验
为了避免变量间存在影响研究结果的严重多重共线性, 本文使用Pearson相关系数法和VIF检验对其进行多重共线性检验, 结果如表8所示。
表 8 农业碳排放研究中变量的Pearson相关系数及VIF检验结果Table 8. Pearson correlation coefficient and VIF test results of variables in the study of agricultural carbon emissions变量
VariableUrban Ey Edu Disaster Area Urban 1.000 Ey 0.551*** 1.000 Edu 0.590*** 0.537*** 1.000 Disaster −0.363*** −0.446*** −0.239*** 1.000 Area 0.285*** 0.399*** 0.296** −0.030 1.000 VIF 2.66 2.35 2.05 1.35 1.26 ***、**分别表示在1%、5%水平显著; Urban、Ey、Edu、Disaster、Area和VIF分别表示城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积和方差膨胀因子。***, ** indicate significance at 1%, 5% levels, respectively. Urban, Ey, Edu, Disaster, Area and VIF represent the urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capital cultivated land area and variance inflation factor. 模型中可能存在严重多重共线性的Pearson相关系数上限值为0.7[35], 而本文回归模型变量的Pearson相关系数绝对值均小于0.6。此外, 方差膨胀因子(VIF)判断是否存在多重共线性的标准为: 1)VIF的最大值小于10; 2)VIF的均值大于1。本文VIF检验结果显示, 变量中VIF最大值为2.66<10, VIF平均值为1.81>1, 满足不存在严重多重共线性的两个条件。综上所述, 本文各变量间不存在严重多重共线性。
3.7 农业碳排放研究中变量的单位根检验与协整检验
针对模型回归中可能出现的伪回归, 参考郭庆旺等[36]的方法, 对各变量进行单位根检验和协整检验。面板数据单位根检验根据面板数据的特征共有5种检验方法, 鉴于本文所使用的面板数据为短面板数据, 应当使用HT检验和IPS检验, 本文额外使用LLC检验意在提高检验结果的可靠性, 所得单位根检验结果如表9所示。
表 9 农业碳排放研究中变量的单位根检验结果Table 9. Results of unit root tests for variables in agricultural carbon emissions research变量 Variable HT检验 HT test IPS检验 IPS test LLC检验 LLC test 平稳性 Stationarity Statistic P Statistic P Statistic P Carbon I 0.5601 0.0045 −4.9399 0.0000 −6.6489 0.000 平稳 Steady ATE 0.9396* 0.0000 −29.9033 0.0000 −3.1050 0.0010 平稳 Steady Urban 0.4386 0.0000 −6.3996 0.0000 −20.9000 0.000 平稳 Steady Ey 0.3116* 0.0000 −7.5412** 0.0000 −6.4405* 0.0000 平稳 Steady Edu 0.1314 0.0000 −7.5447 0.0000 −5.7000 0.0000 平稳 Steady Disaster −0.1156 0.0000 −12.6652 0.0000 −18.5207 0.0000 平稳 Steady Area −0.0036* 0.0000 −5.8085* 0.0000 −13.7496* 0.0000 平稳 Steady *表示进行了一阶差分, **表示进行了二阶差分; HT检验、IPS检验、LLC检验分别表示Harris-Tzavalis检验、Im, Pesaran and Shin检验和Levin-Lin-Chu 检验; Carbon I、ATE、Urban、Ey、Edu、Disaster和Area分别代表农业碳排放强度、农业技术效率、城镇化率、农业经济发展水平、人力资本水平、农业受灾程度和人均耕地面积。* means the first difference, ** means the second difference. HT test, IPS test and LLC test represent Harris-Tzavalis test, Im, Pesaran and Shin test, and Levin-Lin-Chu test. Carbon I, ATE, Urban, Ey, Edu, Disaster and Area represent agricultural carbon emission intensity, agricultural technical efficiency, urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, and per capital cultivated land area, respectively. 根据表9结果, 各变量均通过了单位根检验, 表现出变量的平稳性, 但因对部分变量进行了一阶和二阶差分处理, 而在原始模型中本文仍旧使用原序列进行回归, 因此需对面板数据模型进行协整检验, 所得结果如表10所示。协整检验结果显示Pedroni检验、Kao检验和Westerlund检验均高度显著通过协整检验。综合单位根检验和协整检验结果, 说明变量存在长期均衡关系, 不存在伪回归情况。
表 10 农业碳排放研究中变量的协整检验结果Table 10. Cointegration test results of the variables in the study of agricultural carbon emissions检验类别 Kind of inspection Statistic P 佩德罗尼检验 Pedroni test MPP test 8.0859 0.0000 PP test −4.5221 0.0000 ADF test −3.1695 0.0008 卡奥检验 Kao test MDF test −2.9334 0.0017 DF test −4.4456 0.0000 ADF test −3.8192 0.0001 维斯特隆德检验 Westerlund test VR 3.3257 0.0004 MPP test、PP test、ADF test、MDF test、DF test和VR分别表示修正的菲利普斯-佩荣检验、菲利普斯-佩荣检验、增广迪基-富勒检验、修正的迪基-富勒检验、迪基-富勒检验和方差比。MPP test, PP test, ADF test, MDF test, DF test and VR represent Modified Phillips-Perron test, Phillips-Perron test, Augmented Dickey-Fuller test, Modified Dickey-Fuller test, Dickey-Fuller test and Variance ratio, respectively. 3.8 农业技术效率对农业碳排放强度的门槛效应检验
在门槛效应回归前应对其进行Hausman检验, 在固定效应和随机效应之间做出选择, 结果表明固定效应更适合本模型。因此使用固定效应模型进行门槛效应检验, 确定门槛个数和门槛效应是否存在, 在stata16中输入门槛指令, 并使用Bootstrap反复抽样300次, 所得检验结果如表11所示。
表 11 农业技术效率对农业碳排放强度的门槛效应检验结果Table 11. Test results of threshold effect of agricultural technical efficiency on agricultural carbon emission intensity门槛检验 Threshold test F 临界值 Critical value 10% 5% 1% 单一门槛 Single threshold 197.94*** 46.5580 51.5043 65.6828 双重门槛 Double threshold 118.84*** 31.6226 35.4662 45.8873 三重门槛 Triple threshold 76.39 81.7964 94.8090 111.1802 ***表示在1%水平显著。*** indicates significance at 1% level. 门槛效应检验结果显示模型通过单一门槛和双重门槛检验, 未通过三重门槛检验, 表明农业技术效率与农业碳排放强度存在双重门槛效应。
3.9 农业技术效率对农业碳排放强度的门槛模型估计结果分析
根据门槛效应检验结果可知, 回归模型为双重门槛模型, 据此将公式(10)设置为双门槛, 并对面板数据进行门槛回归, 所得结果如表12所示。
表 12 农业技术效率对农业碳排放强度的门槛模型回归估计结果Table 12. Threshold model regression estimation results of agricultural technical efficiency on agricultural carbon emission intensity变量
Variable系数估计值
Coefficient
estimation标准误
Standard
errorStatistic Urban −0.480*** 0.074 −6.46 Ey −0.011 0.009 −1.15 Edu 0.369 1.008 0.37 Disaster 0.1078*** 0.033 3.30 Area −0.003 0.003 −1.28 ATE (ATE≤0.0746) 0.496 0.309 1.61 ATE (0.0746<ATE≤0.2590) −2.122*** 1.150 −14.16 ATE (ATE>0.2590) −1.538*** 0.130 −11.84 Constant 1.082*** 0.056 19.30 Observations 620 F-value 542.31 R-squared 0.8819 ***表示在1%水平显著; Urban、Ey、Edu、Disaster、Area和ATE分别代表城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积和农业技术效率。*** indicates significance at 1% level. Urban, Ey, Edu, Disaster, Area and ATE represent urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capita cultivated land area and agricultural technical efficiency, respectively. 门槛模型的回归结果表明, 当农业技术效率(ATE)低于第一门槛值0.0746时, 农业技术效率的系数估计值为0.496, 这说明农业技术效率的提升能够提高农业碳排放强度, 但此时结果并不显著; 当农业技术效率水平介于0.0746和0.2590之间时, 此时的系数估计值为−2.122, 这说明当农业技术效率越过拐点后, 其对农业碳排放强度的影响将转变为负向; 当农业技术效率提高到0.2590以上时, 系数估计值会降低到−1.538, 说明农业技术效率达到转变为边际效益递减的临界值后, 农业技术效率对农业碳排放仍然有减排作用, 但其影响力已经不如上一个门槛阶段, 据此假设3得到验证。
以门槛值作为分类标准, 将2015年31个省、直辖市和自治区(港澳台地区以外)农业技术效率进行分类, 所得结果如表13所示。分类结果显示, 技术效率较高(ATE>0.2590)的地区均为东部沿海地区, 中国各省均跨过了第一阶段门槛区间, 进入到农业技术效率有效促进农业碳减排的发展区间。
表 13 农业技术效率门槛值及省份分布(2015年)Table 13. Agricultural technical efficiency threshold and provincial distribution (2015)门槛值及区间
Threshold value and interval省份
Province低 Low (ATE≤0.0746) 中等 Medium
(0.0746<ATE≤0.2590)新疆、西藏、湖南、吉林、四川、河南、甘肃、重庆、陕西、黑龙江、湖北、广西、江西、
云南、山西、贵州、青海、宁夏、安徽、内蒙古
Xinjiang, Tibet, Hunan, Jilin, Sichuan, Henan, Gansu, Chongqing, Shaanxi, Heilongjiang, Hubei, Guangxi, Jiangxi,
Yunnan, Shanxi, Guizhou, Qinghai, Ningxia, Anhui, Inner Mongolia高 High
(ATE>0.2590)北京、上海 、浙江、天津 、福建、广东 、辽宁、海南、山东、河北 、江苏
Beijing, Shanghai, Zhejiang, Tianjin, Fujian, Guangdong, Liaoning, Hainan, Shandong, Hebei, Jiangsu综上所述, 农业技术效率的提升会促进农业碳排放强度的降低, 但其影响是非线性的。在不同区间内, 农业技术效率对农业碳排放的影响程度存在差异。其中, 当农业技术效率低于0.0746时, 技术效率对农业碳排放强度的影响为负向; 当高于0.0746后, 将转变为显著负向效应; 进一步提高农业技术效率并高于0.2590后, 负向效应会减弱。结合表13的分类结果, 截至2015年中国所有省份已经跨过了“拐点”, 未到达第二门槛值的省份应当继续提升农业技术效率, 发挥农业技术在低碳农业发展中的关键性作用。
4. 讨论与结论
4.1 讨论
本文基于港澳台地区以外, 中国31个省、直辖市和自治区2001—2020年的面板数据, 对农业碳排放总量、农业碳排放强度和农业技术效率进行了测算, 剖析了农业碳排放强度和农业技术效率的空间相关性, 利用空间杜宾模型和门槛模型, 分别探究了农业技术效率对农业碳排放强度的空间溢出效应和农业技术效率与农业碳排放强度的非线性关系。在YANG等[4]、杨钧[11]、李成龙等[14]和徐清华等[16]学者的研究基础上, 本文从技术效率的角度出发, 进一步探究了农业技术效率与农业碳排放强度的空间效应和非线性关系。现有的研究大多将研究范围局限于广义的技术进步或狭义的技术进步, 缺少对技术效率进行深入探索, 本文将广义的技术进步分解后的技术效率作为核心, 进一步验证了技术进步对农业碳排放的空间影响机制与边际效应, 为探索实现“双碳”目标的路径提供理论基础与政策依据。
4.2 结论
本文主要结论: 1)总体上各地区农业碳排放总量和强度均表现出下降趋势, 空间分布上中部地区农业碳排放总量高于东西部地区, 东部地区农业碳排放强度低于中西部地区。农业技术效率水平则随时间逐年提高, 东部地区农业技术效率显著高于中西部地区, 而中部地区仅略微高于西部地区。2)农业技术效率与农业碳排放强度具有高度空间相关性并表现出集聚的特征, 聚集类型均以高高聚集和低低聚集为主。3)空间杜宾模型结果表明, 本省农业碳排放强度的提高会促进邻省的农业碳排放强度的上升, 具有正向空间溢出效应, 而农业技术效率对农业碳排放强度则表现为负向空间溢出, 并在总体上有利于抑制农业碳排放。此外城镇化、人力资本水平和人均耕地面积对农业碳排放强度具有负向影响, 农业经济发展水平和农业受灾程度则为正向影响。4)农业技术效率与农业碳排放强度存在双门槛效应, 当农业技术效率越过“拐点”(0.0746)后, 其对农业碳排放强度的影响转变为负向, 当进一步提升农业技术效率水平后, 其影响力会因边际效应递减而减弱。
对策建议: 1)各省市应依照自身农业发展实际情况, 制定符合本省省情的碳减排目标。不拘泥于农业碳排放总量的变化, 而是将农业碳排放强度作为衡量减排成效的标准。2)积极推动建立农业数据库和农作物检测系统, 应用地理信息遥感技术和计算机网络虚拟技术等数字化农业技术, 发挥农业技术的关键性作用, 提高农业技术效率, 推进新型农业技术革命, 推动传统农业迈向现代化农业。3)重视农业生产技术的传播, 将东部地区先进的农业技术手段推广到中西部地区, 促进中西部地区农业生产方式从原本的粗放型转变为产业化、生态化、集约化和低碳化。
-
![]()
图 1 2001—2020年中国不同地区农业碳排放总量(a)和农业碳排放强度(b)
Figure 1. Total amount (a) and intensity (b) of agricultural carbon emissions in different regions of China from 2001 to 2020
![]()
图 2 2001—2020年中国不同地区农业技术效率
Figure 2. Agricultural technical efficiencies in different regions of China from 2001 to 2020
![]()
图 3 2020年农业碳排放强度(a)与农业技术效率(b) Moran’s I散点图
1: 北京; 2: 天津; 3: 河北; 4: 山西; 5: 内蒙古; 6: 辽宁; 7: 吉林; 8: 黑龙江; 9: 上海; 10: 江苏; 11: 浙江; 12: 安徽; 13: 福建; 14: 江西; 15: 山东; 16: 河南; 17: 湖北; 18: 湖南; 19: 广东; 20: 广西; 21: 海南; 22: 重庆; 23: 四川; 24: 贵州; 25: 云南; 26: 西藏; 27: 陕西; 28: 甘肃; 29: 青海; 30: 宁夏; 31: 新疆。
Figure 3. Moran’s I scatter diagrams of agricultural carbon emission intensity and agricultural technical efficiency in 2020
1: Beijing; 2: Tianjin; 3: Hebei; 4: Shanxi; 5: Inner Mongolia; 6: Liaoning; 7: Jilin; 8: Heilongjiang; 9: Shanghai; 10: Jiangsu; 11: Zhejiang; 12: Anhui; 13: Fujian; 14: Jiangxi; 15: Shandong; 16: Henan; 17: Hubei; 18: Hunan; 19: Guangdong; 20: Guangxi; 21: Hainan; 22: Chongqing; 23: Sichuan; 24: Guizhou; 25: Yunnan; 26: Tibet; 27: Shaanxi; 28: Gansu; 29: Qinghai; 30: Ningxia; 31: Xinjiang.
表 1 各类农业碳排放源的碳排放系数及其来源
Table 1 Carbon emission coefficient of each carbon source of agriculture and its reference source
碳源
Carbon source计算方法
Computational method碳排放系数
Carbon emission coefficient数据参考来源
Reference source化肥 Chemical fertilizer 实际化肥投入量 Actual fertilizer input 0.859 kg∙kg−1 West, et al[28], ORNL[29] 农药 Pesticide 实际农药投入量 Actual pesticide inputs 4.934 kg∙kg−1 ORNL[29] 柴油 Diesel 农业机械消耗柴油量 Diesel fuel consumption by agricultural machinery 0.593 kg∙kg−1 IREEA 农膜 Agricultural film 农用塑料薄膜使用量 Amount of plastic film used for agriculture 5.18 kg∙kg−1 IPCC[30] 农地翻耕 Agricultural land tilling 农作物实际播种总面积 Actual total area sown of crops 312.6 kg∙hm−2 CAB[31] 农业灌溉 Agricultural irrigation 农作物实际灌溉面积 Actual irrigated area of crops 20.476 kg∙hm−2 RDRCH[5] ORNL: 美国橡树岭国家实验室; IREEA: 南京农业大学农业资源与生态环境研究所; IPCC: 政府间气候变化专门委员会; CAB: 中国农业大学农学与生物技术学院; RDRCH: 湖北农村发展研究中心。ORNL: Oak Ridge National Laboratory, USA; IREEA: Institute of Resource, Ecosystem and Environment of Agriculture, Nanjing Agricultural University; IPCC: Intergovernmental Panel on Climate Change; CAB: College of Agronomy and Biotechnology, China Agricultural University; RDRCH: Rural Development Research Center of Hubei. 
下载: 导出CSV
表 2 农业碳排放研究相关变量说明及描述性统计
Table 2 Description and descriptive statistics of related variables of agricultural carbon emission research
变量类型
Type of
variable指标名称
Name of
index变量符号
Variable
symbol测算方法
Calculating
method单位
Unit均值
Mean
value标准差
Standard
deviation被解释变量
Explained variable农业碳排放强度
Carbon emission intensity
of agricultureCarbon I 农业碳排放/不变价格的种植业总产值
Agricultural carbon emissions /
gross plantation value at constant pricet·(104 ¥)−1 0.51 0.28 解释变量
Explaining variable农业技术效率
Agricultural technical
efficiencyATE 基于超越对数生产模型计算
Based on translog production model— 0.19 0.10 控制变量
Control variable城镇化率
Urbanization rateUrban 城镇常住人口/常住总人口
Permanent urban population /
total permanent population% 0.50 0.17 农业经济发展水平
Level of agricultural
economic developmentEy 农业总产值/农业从业人员
Total agricultural output value /
agricultural employees×104 ¥·capita−1 1.79 1.30 人力资本水平
Level of human capitalEdu 农村地区初中及其以上学历人数/
农村地区6岁以上人数
Number of people with junior high school
education or above in rural areas / number
of people over 6 years old in rural areas% 0.05 0.01 农业受灾程度
Degree of agricultural
damageDisaster 受灾面积/农作物播种面积
Affected area / crop sown area% 0.22 0.16 人均耕地面积
Arable land per capitaArea 农作物播种面积/第一产业从业人数
Crop sown area / number of
workers in primary industryhm2·capita−1 0.63 0.34 样本时间
Sample timeYear 年份(2001—2020年)
Year (2001—2020)— — — 研究对象
Object of studyProvince 省份(31个省、直辖市、自治区)
Provinces (31 provinces, cities, autonomous regions)
— — —
下载: 导出CSV
表 3 2001—2020年农业碳排放强度和农业技术效率的全局Moran’s I指数及检验
Table 3 Moran’s I indexes and test of agricultural carbon emission intensity and agricultural technical efficiency from 2001 to 2020
年份
Year农业碳排放强度
Agricultural carbon
emission intensity年份
Year农业技术效率
Agricultural technical
efficiencyMoran’s I Z P Moran’s I Z P 2001 0.278 2.636 0.004 2001 0.338 3.562 0.000 2002 0.224 2.196 0.014 2002 0.345 3.592 0.000 2003 0.237 2.287 0.011 2003 0.351 3.619 0.000 2004 0.220 2.151 0.016 2004 0.356 3.643 0.000 2005 0.185 1.840 0.033 2005 0.362 3.665 0.000 2006 0.185 1.866 0.031 2006 0.366 3.684 0.000 2007 0.221 2.131 0.017 2007 0.370 3.701 0.000 2008 0.118 1.267 0.103 2008 0.374 3.717 0.000 2009 0.198 1.971 0.024 2009 0.378 3.73 0.000 2010 0.208 2.056 0.020 2010 0.381 3.742 0.000 2011 0.195 1.931 0.027 2011 0.384 3.752 0.000 2012 0.211 2.059 0.020 2012 0.386 3.761 0.000 2013 0.203 1.999 0.023 2013 0.389 3.769 0.000 2014 0.177 1.826 0.034 2014 0.391 3.776 0.000 2015 0.229 2.258 0.012 2015 0.393 3.782 0.000 2016 0.249 2.485 0.006 2016 0.394 3.787 0.000 2017 0.277 2.644 0.004 2017 0.396 3.792 0.000 2018 0.331 3.189 0.001 2018 0.397 3.795 0.000 2019 0.320 3.103 0.001 2019 0.399 3.798 0.000 2020 0.368 3.519 0.000 2020 0.400 3.801 0.000
下载: 导出CSV
表 4 农业技术效率对农业碳排放强度空间面板回归模型的LM检验、LR检验、Wald检验和Hausman检验结果
Table 4 Results of LM test, LR test, Wald test and Hausman test of the spatial panel regression model of agricultural technical efficiency to agricultural carbon emission intensity
检验类别
Test category检验项目
Inspection itemW1 W2 W3 LM检验 LM test LM(误差)检验 LM (error) test 431.232*** 426.437*** 64.628*** 稳健LM(误差)检验 Robust LM (error) test 300.825*** 370.055*** 50.390*** LM(滞后)检验 LM (lag) test 134.108*** 77.351*** 15.260*** 稳健LM(滞后)检验 Robust LM (lag) test 3.701* 20.969*** 1.023 LR检验 LR test SDM与SLM的卡方检验 SDM&SLM chi2 88.50*** 68.00*** 46.97*** SDM与SEM的卡方检验 SDM&SEM chi2 64.73*** 44.35*** 49.50*** Wald检验 Wald test SDM与SLM的卡方检验 SDM&SLM chi2 27.73*** 36.70*** 31.27*** SDM与SEM的卡方检验 SDM&SEM chi2 24.76*** 21.59*** 18.58*** Hausman检验 Hausman test 显著性卡方检验 Prob>=chi2 54.41*** 73.37*** 124.75*** 拟合度 R2 时间效应 Time 0.0750 0.4271 0.0271 个体效应 Ind 0.8809 0.9098 0.8895 混合效应 Both 0.7442 0.7408 0.7247 模型选择 Model select 个体固定效应
的SDM模型
SDM model of individual
fixed effects个体固定效应
的SDM模型
SDM model of individual
fixed effects个体固定效应
的SDM模型
SDM model of individual
fixed effects***和*分别表示在1%和10%水平显著; W1、W2和W3分别代表邻接空间矩阵、地理距离矩阵和经济距离矩阵; LM检验、LR检验和Wald检验分别代表拉格朗日乘子检验、似然比检验和沃尔德检验; SLM、SEM和SDM分别代表空间滞后模型、空间误差模型和空间杜宾模型; Chi2为卡方检验。*** and * indicate significance at 1% and 10% levels, respectively. W1, W2 and W3 represent adjacency spatial matrix, geographical distance matrix and economic distance matrix, respectively. LM test, LR test and Wald test represent Lagrange multiplier test, likelihood ratio test and Wald test, respectively. SLM, SEM and SDM represent spatial lag model, spatial error model and spatial Durbin model, respectively. Chi2 is the chi-square test.
下载: 导出CSV
表 5 农业技术效率对农业碳排放强度的邻接空间矩阵空间杜宾模型回归结果及效应分解
Table 5 Regression results and effect decomposition of the adjacency space matrix spatial Durbin model of agricultural technical efficiency to agricultural carbon emission intensity
变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 6.203***(0.452) −7.081***(0.473) 5.529***(0.431) −7.747***(0521) −2.217***(0.393) Urban −0.361*** (0.107) 0.095(0.128) −0.390***(0.095) −0.282(0.188) −0.672***(0.176) Ey −0.036*** (0.008) 0.090***(0.012) −0.021**(0.008) 0.160***(0.025) 0.139***(0.029) Edu 0.473(1.138) −4.457***(1.529) −0.416(1.077) −9.736***(2.763) −10.151***(3.004) Disaster 0.079***(0.030) 0.049(0.050) 0.099***(0.030) 0.222**(1.103) 0.321***(0.113) Area 0.007***(0.002) −0.016***(0.004) 0.004*(0.002) −0.029***(0.008) −0.024***(0.009) rho 0.602***(0.034) sigma2_e 0.005***(0.0003) ***、**和*分别表示在1%、5%和10%水平显著; Main意为本地区解释变量对本地区被解释变量的影响程度系数β, Wx意为相邻地区解释变量对本地区被解释变量的影响程度系数θ, Spatial意为相邻地区的被解释变量对本地区被解释变量的影响系数ρ, Variance意为差异系数; 直接效应、间接效应与总效应分别表示利用偏微分法的无偏估计结果, 即本地区解释变量对本地区被解释变量的影响程度, 相邻地区解释变量对本地区被解释变量的影响程度和直接效应与间接效应的总和; ATE、Urban、Ey、Edu、Disaster、Area、rho和sigma2_e分别代表农业技术效率、城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积、空间自相关系数ρ和个体效应的特异误差; 括号内为标准差。***, ** and * indicate significance at 1%, 5% and 10% levels, respectively. Main means the coefficient β of the influence degree of the explanatory variables in the local region on the explained variables in the local region; Wx means the coefficient θ of the influence degree of the explanatory variables in neighboring regions on the explained variables in the local region; Spatial means the influence coefficient ρ of the explained variables in neighboring regions on the explained variables in the local region; and Variance means the difference coefficient. Direct effect, indirect effect and total effect respectively represent the unbiased estimation results using partial differential method, that is, the degree of influence of local explanatory variables on local explained variables, the degree of influence of neighboring explanatory variables on local explained variables and the sum of direct and indirect effects. ATE, Urban, Ey, Edu, Disaster, Area, rho and sigma2_e represent agricultural technical efficiency, urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capital cultivated land area, spatial autocorrelation coefficient ρ and individual effect specific error, respectively. Standard deviations are in parentheses.
下载: 导出CSV
表 6 农业技术效率对农业碳排放强度的地理距离矩阵空间杜宾模型回归结果及效应分解
Table 6 Geographical distance matrix spatial Durbin model regression results and effect decomposition of agricultural technical efficiency to agricultural carbon emission intensity
变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 7.110***(0.390) −7.914***(0.418) 6.672***(0.386) −9.018***(0.530) −2.346***(0.468) Urban −0.363***(0.096) 0.185(0.120) −0.374***(0.088) −0.143(0.199) −0.517***(0.185) Ey −0.040***(0.007) 0.075***(0.015) −0.032***(0.008) 0.139***(0.042) 0.107**(0.045) Edu 0.076(1.040) −3.762**(1.486) −0.468(0.995) −10.322***(3.365) −10.790***(3.579) Disaster 0.072***(0.027) 0.001(0.060) 0.079***(0.026) 0.133(0.154) 0.213(0.159) Area 0.009***(0.002) −0.005(0.056) 0.010***(0.023) 0.002(0.015) 0.011(0.017) rho 0.654***(0.037) sigma2_e 0.004***(0.0002) 同表5注释。Note the same as in Table 5.
下载: 导出CSV
表 7 农业技术效率对农业碳排放强度的经济距离矩阵空间杜宾模型回归结果及效应分解
Table 7 Economic distance matrix spatial Durbin model regression results and effect decomposition of agricultural technical efficiency on agricultural carbon emission intensity
变量
Variable模型估计系数
Main空间矩阵估计系数
Wx空间自相关估计系数
Spatial差异系数
Variance直接效应
Direct effect间接效应
Indirect effect总效应
Total effectATE 6.393***(0.379) −7.741***(0.389) 5.837***(0.370) −8.755***(0.446) −2.918***(0.384) Urban −0.444***(0.093) 0.252**(0.126) −0.445***(0.085) 0.025(0.180) −0.420**(0.166) Ey −0.033***(0.008) 0.126***(0.016) −0.018**(0.008) 0.222***(0.031) 0.204***(0.033) Edu 3.181***(1.098) −10.768***(1.732) 1.943*(1.002) −18.431***(2.876) −16.488***(2.945) Disaster 0.102***(0.027) 0.028(0.054) 0.114***(0.028) 0.169(0.112) 0.282**(0.128) Area 0.010***(0.002) −0.018***(0.004) 0.008***(0.002) −0.027***(0.009) −0.019*(0.010) rho 0.537***(0.035) sigma2_e 0.005***(0.0003) 同表5注释。Note the same in Table 5.
下载: 导出CSV
表 8 农业碳排放研究中变量的Pearson相关系数及VIF检验结果
Table 8 Pearson correlation coefficient and VIF test results of variables in the study of agricultural carbon emissions
变量
VariableUrban Ey Edu Disaster Area Urban 1.000 Ey 0.551*** 1.000 Edu 0.590*** 0.537*** 1.000 Disaster −0.363*** −0.446*** −0.239*** 1.000 Area 0.285*** 0.399*** 0.296** −0.030 1.000 VIF 2.66 2.35 2.05 1.35 1.26 ***、**分别表示在1%、5%水平显著; Urban、Ey、Edu、Disaster、Area和VIF分别表示城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积和方差膨胀因子。***, ** indicate significance at 1%, 5% levels, respectively. Urban, Ey, Edu, Disaster, Area and VIF represent the urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capital cultivated land area and variance inflation factor.
下载: 导出CSV
表 9 农业碳排放研究中变量的单位根检验结果
Table 9 Results of unit root tests for variables in agricultural carbon emissions research
变量 Variable HT检验 HT test IPS检验 IPS test LLC检验 LLC test 平稳性 Stationarity Statistic P Statistic P Statistic P Carbon I 0.5601 0.0045 −4.9399 0.0000 −6.6489 0.000 平稳 Steady ATE 0.9396* 0.0000 −29.9033 0.0000 −3.1050 0.0010 平稳 Steady Urban 0.4386 0.0000 −6.3996 0.0000 −20.9000 0.000 平稳 Steady Ey 0.3116* 0.0000 −7.5412** 0.0000 −6.4405* 0.0000 平稳 Steady Edu 0.1314 0.0000 −7.5447 0.0000 −5.7000 0.0000 平稳 Steady Disaster −0.1156 0.0000 −12.6652 0.0000 −18.5207 0.0000 平稳 Steady Area −0.0036* 0.0000 −5.8085* 0.0000 −13.7496* 0.0000 平稳 Steady *表示进行了一阶差分, **表示进行了二阶差分; HT检验、IPS检验、LLC检验分别表示Harris-Tzavalis检验、Im, Pesaran and Shin检验和Levin-Lin-Chu 检验; Carbon I、ATE、Urban、Ey、Edu、Disaster和Area分别代表农业碳排放强度、农业技术效率、城镇化率、农业经济发展水平、人力资本水平、农业受灾程度和人均耕地面积。* means the first difference, ** means the second difference. HT test, IPS test and LLC test represent Harris-Tzavalis test, Im, Pesaran and Shin test, and Levin-Lin-Chu test. Carbon I, ATE, Urban, Ey, Edu, Disaster and Area represent agricultural carbon emission intensity, agricultural technical efficiency, urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, and per capital cultivated land area, respectively.
下载: 导出CSV
表 10 农业碳排放研究中变量的协整检验结果
Table 10 Cointegration test results of the variables in the study of agricultural carbon emissions
检验类别 Kind of inspection Statistic P 佩德罗尼检验 Pedroni test MPP test 8.0859 0.0000 PP test −4.5221 0.0000 ADF test −3.1695 0.0008 卡奥检验 Kao test MDF test −2.9334 0.0017 DF test −4.4456 0.0000 ADF test −3.8192 0.0001 维斯特隆德检验 Westerlund test VR 3.3257 0.0004 MPP test、PP test、ADF test、MDF test、DF test和VR分别表示修正的菲利普斯-佩荣检验、菲利普斯-佩荣检验、增广迪基-富勒检验、修正的迪基-富勒检验、迪基-富勒检验和方差比。MPP test, PP test, ADF test, MDF test, DF test and VR represent Modified Phillips-Perron test, Phillips-Perron test, Augmented Dickey-Fuller test, Modified Dickey-Fuller test, Dickey-Fuller test and Variance ratio, respectively.
下载: 导出CSV
表 11 农业技术效率对农业碳排放强度的门槛效应检验结果
Table 11 Test results of threshold effect of agricultural technical efficiency on agricultural carbon emission intensity
门槛检验 Threshold test F 临界值 Critical value 10% 5% 1% 单一门槛 Single threshold 197.94*** 46.5580 51.5043 65.6828 双重门槛 Double threshold 118.84*** 31.6226 35.4662 45.8873 三重门槛 Triple threshold 76.39 81.7964 94.8090 111.1802 ***表示在1%水平显著。*** indicates significance at 1% level.
下载: 导出CSV
表 12 农业技术效率对农业碳排放强度的门槛模型回归估计结果
Table 12 Threshold model regression estimation results of agricultural technical efficiency on agricultural carbon emission intensity
变量
Variable系数估计值
Coefficient
estimation标准误
Standard
errorStatistic Urban −0.480*** 0.074 −6.46 Ey −0.011 0.009 −1.15 Edu 0.369 1.008 0.37 Disaster 0.1078*** 0.033 3.30 Area −0.003 0.003 −1.28 ATE (ATE≤0.0746) 0.496 0.309 1.61 ATE (0.0746<ATE≤0.2590) −2.122*** 1.150 −14.16 ATE (ATE>0.2590) −1.538*** 0.130 −11.84 Constant 1.082*** 0.056 19.30 Observations 620 F-value 542.31 R-squared 0.8819 ***表示在1%水平显著; Urban、Ey、Edu、Disaster、Area和ATE分别代表城镇化率、农业经济发展水平、人力资本水平、农业受灾程度、人均耕地面积和农业技术效率。*** indicates significance at 1% level. Urban, Ey, Edu, Disaster, Area and ATE represent urbanization rate, agricultural economic development level, human capital level, agricultural disaster degree, per capita cultivated land area and agricultural technical efficiency, respectively.
下载: 导出CSV
表 13 农业技术效率门槛值及省份分布(2015年)
Table 13 Agricultural technical efficiency threshold and provincial distribution (2015)
门槛值及区间
Threshold value and interval省份
Province低 Low (ATE≤0.0746) 中等 Medium
(0.0746<ATE≤0.2590)新疆、西藏、湖南、吉林、四川、河南、甘肃、重庆、陕西、黑龙江、湖北、广西、江西、
云南、山西、贵州、青海、宁夏、安徽、内蒙古
Xinjiang, Tibet, Hunan, Jilin, Sichuan, Henan, Gansu, Chongqing, Shaanxi, Heilongjiang, Hubei, Guangxi, Jiangxi,
Yunnan, Shanxi, Guizhou, Qinghai, Ningxia, Anhui, Inner Mongolia高 High
(ATE>0.2590)北京、上海 、浙江、天津 、福建、广东 、辽宁、海南、山东、河北 、江苏
Beijing, Shanghai, Zhejiang, Tianjin, Fujian, Guangdong, Liaoning, Hainan, Shandong, Hebei, Jiangsu
下载: 导出CSV
-
[1] 王一鸣. 百年大变局、高质量发展与构建新发展格局[J]. 管理世界, 2020, 36(12): 1−13 WANG Y M. Changes unseen in a century, high-quality development, and the construction of a new development pattern[J]. Management World, 2020, 36(12): 1−13
[2] 黄杰, 孙自敏. 中国种植业碳生产率的区域差异及分布动态演进[J]. 农业技术经济, 2022(7): 109−127 HUANG J, SUN Z M. Regional differences and distribution dynamics of carbon productivity in planting industry in China[J]. Journal of Agrotechnical Economics, 2022(7): 109−127
[3] 李劼, 徐晋涛. 我国农业低碳技术的减排潜力分析[J]. 农业经济问题, 2022, 43(3): 117−135 LI J, XU J T. Analyses of carbon reduction potential of low carbon technologies in China[J]. Issues in Agricultural Economy, 2022, 43(3): 117−135
[4] YANG L S, LI Z. Technology advance and the carbon dioxide emission in China - Empirical research based on the rebound effect[J]. Energy Policy, 2017, 101(2): 150−161
[5] 李波, 张俊飚, 李海鹏. 中国农业碳排放时空特征及影响因素分解[J]. 中国人口·资源与环境, 2011, 21(8): 80−86 LI B, ZHANG J B, LI H P. Research on spatial-temporal characteristics and affecting factors decomposition of agricultural carbon emission in China[J]. China Population, Resources and Environment, 2011, 21(8): 80−86
[6] 刘杨, 刘鸿斌. 山东省农业碳排放特征、影响因素及达峰分析[J]. 中国生态农业学报(中英文), 2022, 30(4): 558−569 LIU Y, LIU H B. Characteristics, influencing factors and peak analysis of agricultural carbon emissions in Shandong Province[J]. Chinese Journal of Eco-Agriculture, 2022, 30(4): 558−569
[7] ZHOU P, ANG B W, HAN J Y. Total factor carbon emission performance: a Malmquist index analysis[J]. Energy Economics, 2010, 32(1): 194−201
[8] 张永强, 田媛, 王珧, 等. 农村人力资本、农业技术进步与农业碳排放[J]. 科技管理研究, 2019, 39(14): 266−274 ZHANG Y Q, TIAN Y, WANG Y, et al. Rural human capital, agricultural technology progress and agricultural carbon emissions[J]. Science and Technology Management Research, 2019, 39(14): 266−274
[9] 庞丽. 我国农业碳排放的区域差异与影响因素分析[J]. 干旱区资源与环境, 2014, 28(12): 1−7 PANG L. Empirical study of regional carbon emissions of agriculture in China[J]. Journal of Arid Land Resources and Environment, 2014, 28(12): 1−7
[10] 邓祥征, 钟海玥, 白雪梅, 等. 中国西部城镇化可持续发展路径的探讨[J]. 中国人口·资源与环境, 2013, 23(10): 24−30 DENG X Z, ZHONG H Y, BAI X M, et al. Path of sustainable urbanization in Western China[J]. China Population, Resources and Environment, 2013, 23(10): 24−30
[11] 杨钧. 农业技术进步对农业碳排放的影响−中国省级数据的检验[J]. 软科学, 2013, 27(10): 116−120 YANG J. The effects of technological advances on agricultural carbon emission — Evidence from Chinese provincial data[J]. Soft Science, 2013, 27(10): 116−120
[12] VALIN H, HAVLÍK P, MOSNIER A, et al. Agricultural productivity and greenhouse gas emissions: trade-offs or synergies between mitigation and food security?[J]. Environmental Research Letters, 2013, 8(3): 035019 doi: 10.1088/1748-9326/8/3/035019
[13] 周勇, 林源源. 技术进步对能源消费回报效应的估算[J]. 经济学家, 2007(2): 45−52 doi: 10.3969/j.issn.1003-5656.2007.02.007 ZHOU Y, LIN Y Y. The estimation of technological progress on the energy consumption returns effects[J]. Economist, 2007(2): 45−52 doi: 10.3969/j.issn.1003-5656.2007.02.007
[14] 李成龙, 周宏. 农业技术进步与碳排放强度关系−不同影响路径下的实证分析[J]. 中国农业大学学报, 2020, 25(11): 162−171 doi: 10.11841/j.issn.1007-4333.2020.11.17 LI C L, ZHOU H. Relationship between agricultural technology progress and carbon emission intensity: an empirical analysis under different influence paths[J]. Journal of China Agricultural University, 2020, 25(11): 162−171 doi: 10.11841/j.issn.1007-4333.2020.11.17
[15] 丁超. 三角贸易及其对中国劳动生产率的影响研究[D]. 厦门: 厦门大学, 2014 DING C. Study on present situation of triangle trade and its influence on China’s labor productivity[D]. Xiamen: Xiamen University, 2014
[16] 徐清华, 张广胜. 农业机械化对农业碳排放强度影响的空间溢出效应−基于282个城市面板数据的实证[J]. 中国人口·资源与环境, 2022, 32(4): 23−33 XU Q H, ZHANG G S. Spatial spillover effect of agricultural mechanization on agricultural carbon emission intensity: an empirical analysis of panel data from 282 cities[J]. China Population, Resources and Environment, 2022, 32(4): 23−33
[17] 陈军娟, 燕振刚, 李薇, 等. 基于系统动力学的民勤绿洲农业系统碳排放仿真模拟研究[J]. 西南农业学报, 2022, 35(6): 1432−1440 CHEN J J, YAN Z G, LI W, et al. Simulation study on carbon emission of Minqin oasis agricultural system based on system dynamics[J]. Southwest China Journal of Agricultural Sciences, 2022, 35(6): 1432−1440
[18] 何炫蕾, 陈兴鹏, 庞家幸. 基于LMDI的兰州市农业碳排放现状及影响因素分析[J]. 中国农业大学学报, 2018, 23(7): 150−158 HE X L, CHEN X P, PANG J X. Current situation and influencing factors of agricultural carbon emissions in Lanzhou based on LMDI[J]. Journal of China Agricultural University, 2018, 23(7): 150−158
[19] 贺青, 张虎, 张俊飚. 农业产业聚集对农业碳排放的非线性影响[J]. 统计与决策, 2021, 37(9): 75−78 HE Q, ZHANG H, ZHANG J B. Nonlinear effects of agricultural industry aggregation on agricultural carbon emissions[J]. Statistics & Decision, 2021, 37(9): 75−78
[20] 张抗私, 史策. 认知能力、技术进步与就业极化[J]. 现代财经(天津财经大学学报), 2022, 42(5): 95−113 ZHANG K S, SHI C. Cognitive ability, technological progress and employment polarization[J]. Modern Finance and Economics - Journal of Tianjin University of Finance and Economics, 2022, 42(5): 95−113
[21] 杨增旭, 韩洪云. 化肥施用技术效率及影响因素−基于小麦和玉米的实证分析[J]. 中国农业大学学报, 2011, 16(1): 140−147 YANG Z X, HAN H Y. Technical efficiency of fertilizer and its influencing factors: based on wheat and corn empirical study[J]. Journal of China Agricultural University, 2011, 16(1): 140−147
[22] 汪小勤, 姜涛. 基于农业公共投资视角的中国农业技术效率分析[J]. 中国农村经济, 2009(5): 79−86 WANG X Q, JIANG T. Analysis of China’s agricultural technical efficiency from the perspective of agricultural public investment[J]. Chinese Rural Economy, 2009(5): 79−86
[23] 高鸣, 宋洪远. 粮食生产技术效率的空间收敛及功能区差异−兼论技术扩散的空间涟漪效应[J]. 管理世界, 2014(7): 83−92 GAO M, SONG H Y. Spatial convergence and functional zone differences of technical efficiency of grain production: a discussion on the spatial ripple effect of technology diffusion[J]. Management World, 2014(7): 83−92
[24] 张思麒, 刘导波. 技术进步视角下中国产业结构高级化格局及影响因素[J]. 经济地理, 2022, 42(5): 104−113 ZHANG S Q, LIU D B. Spatial pattern and influencing factors of industrial structure advancement in China based on technological progress[J]. Economic Geography, 2022, 42(5): 104−113
[25] 杨义武, 林万龙, 张莉琴. 农业技术进步、技术效率与粮食生产−来自中国省级面板数据的经验分析[J]. 农业技术经济, 2017(5): 46−56 YANG Y W, LIN W L, ZHANG L Q. Agricultural technical progress, technical efficiency and grain production: an empirical analysis from provincial panel data in China[J]. Journal of Agrotechnical Economics, 2017(5): 46−56
[26] WANG Q Y. Fixed-effect panel threshold model using stata[J]. The Stata Journal, 2015, 15(1): 121−134 doi: 10.1177/1536867X1501500108
[27] HANSEN B E. Threshold effects in non-dynamic panels: Estimation, testing, and inference[J]. Journal of Econometrics, 1999, 93(2): 345−368 doi: 10.1016/S0304-4076(99)00025-1
[28] 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[J]. Agriculture Ecosystems & Environment, 2002, 91(1/3): 217−232
[29] 智静, 高吉喜. 中国城乡居民食品消费碳排放对比分析[J]. 地理科学进展, 2009, 28(3): 429−434 ZHI J, GAO J X. Comparative analysis of carbon emissions from food consumption of urban and rural residents in China[J]. Progress in Geography, 2009, 28(3): 429−434
[30] IPCC. Climate Change 2007: The Physical Science Basis: Working Group I Contribution to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change[M]. New York: Cambridge University Press, 2007
[31] 伍芬琳, 李琳, 张海林, 等. 保护性耕作对农田生态系统净碳释放量的影响[J]. 生态学杂志, 2007(12): 2035−2039 WU F L, LI L, ZHANG H L, et al. Effects of conservation tillage on net carbon release from farmland ecosystem[J]. Chinese Journal of Ecology, 2007(12): 2035−2039
[32] 傅晓霞, 吴利学. 前沿分析方法在中国经济增长核算中的适用性[J]. 世界经济, 2007, 30(7): 56−66 FU X X, WU L X. Applicability of frontier analysis method in China’s economic growth accounting[J]. The Journal of World Economy, 2007, 30(7): 56−66
[33] 白俊红, 江可申, 李婧. 应用随机前沿模型评测中国区域研发创新效率[J]. 管理世界, 2009(10): 51−61 BAI J H, JIANG K S, LI J. Stochastic frontier model is applied to evaluate the innovation efficiency of regional R&D in China[J]. Management World, 2009(10): 51−61
[34] LESAGE J, PACE R K. Introduction to Spatial Econometrics[M]. Boca Raton: CRC Press, 2009: 155–165
[35] 赵昕东, 刘成坤. 人口老龄化对制造业结构升级的作用机制研究−基于中介效应模型的检验[J]. 中国软科学, 2019(3): 153−163 ZHAO X D, LIU C K. Research on the mechanism of population aging to the upgrading of manufacturing structure — Based on the test of mediation effect model[J]. China Soft Science, 2019(3): 153−163
[36] 郭庆旺, 贾俊雪. 中国全要素生产率的估算: 1979—2004[J]. 经济研究, 2005, 40(6): 51−60 GUO Q W, JIA J X. Estimating total factor productivity in China: 1979–2004[J]. Economic Research Journal, 2005, 40(6): 51−60
-
期刊类型引用(27)
1. 张凯丽，张风丽. “双碳”目标下兵团农业碳排放核算及驱动因素研究. 新疆农垦经济. 2025(02): 73-83 . 
百度学术
2. 李国志，江梦滢，颜诗旋. 数字普惠金融对农业碳排放的影响研究. 甘肃金融. 2024(03): 8-14 . 百度学术
3. 周佩萱. 动员、合作与整合：村庄产业发展能力提升机制研究——以陕西省裕村苍术产业发展为例. 山东行政学院学报. 2024(02): 73-80 . 百度学术
4. 刘小妍. 长江经济带城乡基本公共服务均等化对城乡融合的影响研究. 浙江工贸职业技术学院学报. 2024(01): 85-92 . 百度学术
5. 周晨阳，张凌宇，韦宏耀. 农业保险的增产减排效应研究：能否兼顾低碳发展. 金融理论与实践. 2024(03): 106-118 . 百度学术
6. 李宽，史磊. 农村产业融合对农业碳排放的影响：机制路径及空间溢出效应分析. 中国农业资源与区划. 2024(04): 1-14 . 百度学术
7. 夏文浩，霍瑜，逯渊，王超毅. 新疆农业碳排放的时空差异与空间溢出效应分析. 干旱区地理. 2024(06): 1084-1096 . 百度学术
8. 赵金彩，杜永猛，任世鑫，段锴丰，任芮莹. 乡村振兴对黄河流域农业碳排放的空间溢出及门槛效应研究. 中国生态农业学报(中英文). 2024(10): 1766-1779 . 本站查看
9. 刘建华，闫静，王慧扬，施天乐. 黄河流域农业科技创新对农业碳排放强度影响研究. 生态经济. 2024(12): 118-127 . 百度学术
10. 李宽，史磊，张弘. 我国新型农业经营主体发展对农业碳排放强度的影响：“减碳效应”或“增碳效应”. 农业技术经济. 2024(11): 51-73 . 百度学术
11. 卢泓钢，王益如，吴伟光. 中国农业减污降碳协同增效的水平测度及其时空收敛性研究——基于“减污、降碳、扩绿、增长”的四维分析. 生态经济. 2024(12): 186-195 . 百度学术
12. 刘浩，刘芮萌，欧阳博强. “双碳”背景下韶关市农业碳中和效应的测度与动态评价. 安徽农业科学. 2024(23): 205-210 . 百度学术
13. 喻芳玲，鲍宏礼，吴学兵. 新型城镇化对农业碳排放的空间效应研究. 安徽农业科学. 2024(24): 203-207 . 百度学术
14. 张灵蕤，刘辉，邓岚，李群. “双碳”目标下我国农林业碳排放效率的时空演变及影响因素分析. 林业经济. 2024(08): 59-83 . 百度学术
15. 白鹏飞，黄文茜. 新质生产力对农业碳排放强度的影响研究. 江西农业学报. 2024(12): 80-92 . 百度学术
16. 尹婧文，申澳，刘晴兮，邓桂雨，付帅. 基于LMDI模型的江西省农业水资源利用——碳排放耦合分析及碳减排优化策略研究. 科技经济市场. 2023(03): 100-102 . 百度学术
17. 李菲菲，徐绘薇，王舒泓，崔金栋. 吉林省农业碳排放测度分析及影响因素研究. 综合智慧能源. 2023(08): 36-43 . 百度学术
18. 张千禧，曹智，王介勇. 双碳视角下乡村地域系统碳效应研究进展. 中国生态农业学报(中英文). 2023(08): 1256-1274 . 本站查看
19. 胡池群，马晓钰，刘家民. 城乡融合对农业碳排放强度的影响研究. 农业现代化研究. 2023(04): 668-679 . 百度学术
20. 蔡祖梅，黄森慰，李学渊，洪月苇，刘彦. 农业机械化水平对农业碳排放强度的空间溢出效应和门槛效应分析. 福建农林大学学报(哲学社会科学版). 2023(05): 21-29 . 百度学术
21. 陈卫洪，钱海慧，颜光耀. 畜牧业产业集聚对畜牧业碳排放的影响——基于库兹涅茨曲线与空间间接效应分析. 黑龙江畜牧兽医. 2023(20): 10-21+28 . 百度学术
22. 朱前涛，秦昊. 财政支农强度、农业碳排放及空间溢出效应. 天津商业大学学报. 2023(06): 10-18 . 百度学术
23. 朱嘉晴，秦会艳，张梦春. 基于时变参数C-D生产函数的江西省农业碳排放驱动因素及脱钩效应. 应用生态学报. 2023(11): 3085-3094 . 百度学术
24. 刘凡，周书灵. 碳排放视角下的农业生态效率研究. 宿州学院学报. 2023(10): 31-37 . 百度学术
25. 熊飞雪，赵星磊，郭子毅，朱述斌. 土地整治对农业碳排放的影响研究——基于高标准农田建设政策的准自然实验. 中国生态农业学报(中英文). 2023(12): 2022-2032 . 本站查看
26. 朱前涛，秦昊. 土地生产率与农业碳排放：非线性关系及空间溢出. 农林经济管理学报. 2023(06): 770-779 . 百度学术
27. 程娜，桑一铭，李博文. “双碳”目标下中国“碳解锁”发展研究. 改革. 2023(12): 151-162 . 百度学术
其他类型引用(21)
计量
- 文章访问数: 1267
- HTML全文浏览量: 316
- PDF下载量: 122
- 被引次数: 48
