黃元生，李兆蓓

（華北電力大學(xué) 經(jīng)濟(jì)管理系, 河北 保定 071000）

隨著中國(guó)“碳達(dá)峰”和“碳中和”目標(biāo)的提出，碳減排成為各產(chǎn)業(yè)所面臨的重要問(wèn)題。根據(jù)國(guó)家統(tǒng)計(jì)年鑒的產(chǎn)業(yè)劃分，本文所研究的農(nóng)業(yè)主要指種植業(yè)，不包含林業(yè)、牧業(yè)、漁業(yè)等其他產(chǎn)業(yè)[1]。相比于其他產(chǎn)業(yè)，農(nóng)業(yè)在碳減排方面具有特殊性：一方面，在農(nóng)業(yè)生產(chǎn)過(guò)程中，化肥、農(nóng)藥以及農(nóng)機(jī)的使用都會(huì)產(chǎn)生碳排放；另一方面，農(nóng)業(yè)作物又可以通過(guò)植物細(xì)胞自身代謝過(guò)程吸收二氧化碳，實(shí)現(xiàn)固碳的作用[2]。根據(jù)河北省統(tǒng)計(jì)年鑒的數(shù)據(jù)，2020 年農(nóng)業(yè)在河北省第一產(chǎn)業(yè)GDP 中占比為50.63%，是第一產(chǎn)業(yè)中最主要的產(chǎn)業(yè)。本文對(duì)河北省過(guò)去20 年的農(nóng)業(yè)凈碳匯進(jìn)行測(cè)算，并從內(nèi)部驅(qū)動(dòng)和外部環(huán)境兩個(gè)方面，對(duì)影響河北省凈碳匯的因素進(jìn)行了分析，為河北省發(fā)展農(nóng)業(yè)碳匯提出了針對(duì)性的建議。

一、相關(guān)文獻(xiàn)綜述

近年來(lái)人們對(duì)于碳匯的作用越來(lái)越重視，針對(duì)農(nóng)業(yè)碳匯的研究也越來(lái)越多。部分研究者通過(guò)不同的方法對(duì)農(nóng)業(yè)碳匯進(jìn)行了測(cè)算。

胡歡等在研究中，通過(guò)建立碳源碳匯科學(xué)核算體系，對(duì)河北省碳匯進(jìn)行了測(cè)算，并利用灰色預(yù)測(cè)模型對(duì)2014—2030 年的碳排放強(qiáng)度進(jìn)行預(yù)測(cè)，據(jù)此對(duì)碳減排壓力進(jìn)行了分析[3]。劉愛(ài)玉等在研究中，根據(jù)北京市的土地資源稟賦，對(duì)北京市森林、耕地、草地、濕地的碳匯量進(jìn)行了測(cè)算，加總得到北京市2012 年的碳匯總量，并對(duì)北京市發(fā)展碳匯，控制碳減排提出了政策建議[4]。馬濤對(duì)上海市的碳源與碳匯情況進(jìn)行測(cè)算，認(rèn)為上海市具有很大的碳匯潛力，并提出將農(nóng)業(yè)碳匯納入碳交易市場(chǎng)的建議，以此來(lái)改善環(huán)境[5]。

農(nóng)業(yè)不僅僅是碳匯，在其生產(chǎn)活動(dòng)中同樣會(huì)排放二氧化碳，因此部分學(xué)者在測(cè)算農(nóng)業(yè)碳匯量的基礎(chǔ)上，剔除掉農(nóng)業(yè)生產(chǎn)中的碳排放量，得到凈碳匯值。田云等對(duì)中國(guó)農(nóng)業(yè)凈碳匯狀況進(jìn)行了計(jì)算與分析，并對(duì)其效應(yīng)進(jìn)行了研究，得出結(jié)論：我國(guó)農(nóng)業(yè)凈碳匯量保持上升但區(qū)域凈碳效應(yīng)差異明顯[6]。SHE 等通過(guò)收集中國(guó)主要作物的相關(guān)數(shù)據(jù)，測(cè)算了碳匯量及碳排放量，并計(jì)算了水稻、小麥和玉米三種主要作物的年凈碳匯總量，并根據(jù)測(cè)算結(jié)果對(duì)中國(guó)低碳農(nóng)業(yè)發(fā)展提出了針對(duì)性的政策建議[7]。

農(nóng)業(yè)凈碳匯可以反映出農(nóng)業(yè)對(duì)于碳減排的貢獻(xiàn)度，針對(duì)凈碳匯測(cè)算方法的研究對(duì)于農(nóng)業(yè)碳減排具有重要意義。農(nóng)業(yè)碳排放和碳匯是一個(gè)相對(duì)復(fù)雜系統(tǒng)，受到較多因素的驅(qū)動(dòng)，為了探究農(nóng)業(yè)凈碳匯的驅(qū)動(dòng)因素及其作用效果，部分學(xué)者利用一定方法對(duì)農(nóng)業(yè)凈碳匯的影響因素進(jìn)行分析研究。楊果等先是對(duì)中國(guó)凈碳匯量進(jìn)行了測(cè)算，后基于耦合性模型分析了凈碳匯與農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的關(guān)系[8]。CUI 等利用1997—2017 年中國(guó)31 個(gè)省的面板數(shù)據(jù)，結(jié)合碳排放核算公式、泰爾指數(shù)法和收斂模型，對(duì)碳排放的影響因素進(jìn)行了分析，結(jié)果表明：碳排放總量呈下降趨勢(shì)，其變化主要受到農(nóng)村人均GDP、農(nóng)業(yè)財(cái)政支持水平、農(nóng)業(yè)機(jī)械化程度、農(nóng)業(yè)結(jié)構(gòu)、農(nóng)村人口規(guī)模和城鎮(zhèn)化水平的影響，但是影響效應(yīng)會(huì)根據(jù)區(qū)域不同存在一定的差異[9]。這些研究多數(shù)是針對(duì)少數(shù)指標(biāo)對(duì)于凈碳匯的影響關(guān)系進(jìn)行分析，在指標(biāo)的系統(tǒng)性和全面性上還有深入研究的空間。

結(jié)合現(xiàn)階段國(guó)內(nèi)外農(nóng)業(yè)凈碳匯的研究動(dòng)態(tài)，本文對(duì)河北省農(nóng)業(yè)凈碳匯進(jìn)行了測(cè)算及影響因素分析。本文主要的研究工作及創(chuàng)新點(diǎn)如下：（1）利用碳排放系數(shù)法對(duì)河北省農(nóng)業(yè)碳排放進(jìn)行測(cè)算，采取作物經(jīng)濟(jì)系數(shù)和碳吸收率對(duì)農(nóng)業(yè)碳匯進(jìn)行了測(cè)算，基于此對(duì)河北省2000—2019 年的凈碳匯進(jìn)行測(cè)算及分析；（2）從內(nèi)部驅(qū)動(dòng)和外部環(huán)境兩個(gè)角度選取了若干影響因素，利用VECM 模型和脈沖響應(yīng)對(duì)凈碳匯的影響因素及其作用效果進(jìn)行了分析，并據(jù)此為河北省發(fā)展農(nóng)業(yè)凈碳匯提出參考性的政策建議。

二、研究方法與數(shù)據(jù)來(lái)源

（一） 河北省農(nóng)業(yè)碳排放測(cè)算

本文在參考諸多相關(guān)學(xué)者的碳排放測(cè)算方法后，選擇使用碳排放系數(shù)法[10]來(lái)對(duì)河北省的農(nóng)業(yè)碳排放進(jìn)行測(cè)算，具體公式如下所示：

其中，C表示農(nóng)業(yè)碳排放總量； Ci表示第i 類(lèi)碳源的排放量；Ti表示第i 類(lèi)碳源的規(guī)模；Qi代表第i 類(lèi)碳源的碳排放系數(shù)。

通常農(nóng)業(yè)生產(chǎn)過(guò)程中的碳排放主要來(lái)自于以下幾個(gè)渠道：首先是在農(nóng)業(yè)種植中，生產(chǎn)和使用化肥、農(nóng)藥和農(nóng)用塑料薄膜的過(guò)程中產(chǎn)生的碳排放；其次，在農(nóng)業(yè)的生產(chǎn)過(guò)程中，機(jī)械化作業(yè)燃燒農(nóng)用柴油也會(huì)產(chǎn)生碳排放；此外，土地翻耕會(huì)導(dǎo)致土壤中有機(jī)質(zhì)流失，同樣會(huì)釋放一部分碳源；最后是為灌溉供給能源的燃料燃燒會(huì)產(chǎn)生二氧化碳。根據(jù)IPCC 發(fā)布的數(shù)據(jù)，以上六種碳源的碳排放系數(shù)如表1 所示[11]。

表1 不同碳源碳排放系數(shù)

（二） 河北省農(nóng)業(yè)碳匯測(cè)算

本文對(duì)于農(nóng)業(yè)碳匯的測(cè)算工作主要針對(duì)農(nóng)作物在生長(zhǎng)過(guò)程中所吸收的碳排放，不考慮耕地土壤對(duì)碳的吸收[12]。其中，農(nóng)作物的碳吸收量指植物的生物產(chǎn)量。通過(guò)參考李克讓[13]的研究工作，本文利用農(nóng)作物的經(jīng)濟(jì)產(chǎn)量和碳吸收率測(cè)算不同作物的碳匯量。具體計(jì)算公式如下：

其中，C 為農(nóng)業(yè)碳匯總量；Ci為單一種類(lèi)農(nóng)產(chǎn)品的碳匯量；k 為農(nóng)作物種類(lèi)數(shù)；ci為農(nóng)作物碳吸收率；Yi為農(nóng)作物的經(jīng)濟(jì)產(chǎn)量（本文使用農(nóng)作物產(chǎn)量進(jìn)行表示）；r 為農(nóng)作物經(jīng)濟(jì)產(chǎn)量部分的含水量；HIi為農(nóng)作物的經(jīng)濟(jì)系數(shù)。其中，本文中各類(lèi)農(nóng)作物的碳吸收率及經(jīng)濟(jì)系數(shù)參照王修蘭[14]的研究工作，具體數(shù)值如表2 所示。

表2 經(jīng)濟(jì)作物經(jīng)濟(jì)系數(shù)與碳吸收率匯總表

（三） 數(shù)據(jù)說(shuō)明及來(lái)源

本文所研究的農(nóng)用化肥、農(nóng)藥、農(nóng)用塑料薄膜、農(nóng)用柴油、農(nóng)業(yè)翻耕、灌溉這六種碳源分別使用：農(nóng)用化肥施用折純量、農(nóng)藥使用量、農(nóng)用塑料薄膜使用量、農(nóng)用柴油使用量、農(nóng)作物總播種面積、耕地灌溉面積進(jìn)行表征。

在參考大量文獻(xiàn)與資料以及河北省農(nóng)作物實(shí)際種植情況后，本文選擇水稻、小麥、玉米、棉花、油料、蔬菜、豆類(lèi)、薯類(lèi)、甜菜這九類(lèi)農(nóng)作物的產(chǎn)量來(lái)對(duì)河北省的碳匯量進(jìn)行測(cè)算。

本文根據(jù)河北省2000—2019 年樣本數(shù)據(jù)進(jìn)行分析研究，使用的所有原始數(shù)據(jù)均來(lái)自于國(guó)家統(tǒng)計(jì)局以及河北省統(tǒng)計(jì)年鑒，期間個(gè)別缺失數(shù)據(jù)采用趨勢(shì)外推法進(jìn)行補(bǔ)充。

（四） 實(shí)證方法

1. ADF 檢驗(yàn)

為避免非平穩(wěn)時(shí)間序列的變量在協(xié)整以及因果檢驗(yàn)過(guò)程中出現(xiàn)偽回歸或是因果關(guān)系錯(cuò)誤的情況，在進(jìn)行VAR 模型或VECM 模型分析前，需對(duì)序列進(jìn)行單位根檢驗(yàn)來(lái)判斷序列是否平穩(wěn)。本文使用ADF 檢驗(yàn)判斷單位根存在與否。

2. VECM 動(dòng)態(tài)模型

VECM（向量誤差修正）模型是建立了協(xié)整約束的VAR 模型，作為時(shí)間序列模型既可以表達(dá)時(shí)間序列變量之間的相互作用，同時(shí)又能匯報(bào)出長(zhǎng)期均衡結(jié)果?？紤]到本文選取的外部影響因素和內(nèi)部影響因素與河北省農(nóng)業(yè)凈碳匯之間具有時(shí)間上的自相關(guān)性，且各個(gè)變量之間相互影響、相互解釋?zhuān)疚臉?gòu)建VECM 模型對(duì)農(nóng)業(yè)凈碳匯做進(jìn)一步研究。構(gòu)造VECM 模型前要對(duì)協(xié)整的存在性和模型滯后階數(shù)進(jìn)行檢驗(yàn)。協(xié)整方程為：

VECM 模型為：

上式中，α1和α2表示調(diào)整速度。

3. 脈沖響應(yīng)函數(shù)

脈沖響應(yīng)是分析輸入變量和輸出變量響應(yīng)關(guān)系的函數(shù)。在連續(xù)時(shí)間的系統(tǒng)中，如果輸入的脈沖函數(shù)為δ（t），則系統(tǒng)輸出的函數(shù)h（t）被稱(chēng)為脈沖響應(yīng)函數(shù)。在實(shí)際應(yīng)用中，輸入信號(hào)經(jīng)常使用偽隨機(jī)信號(hào)，輸入輸出的互相關(guān)函數(shù)Ruy（t）可通過(guò)計(jì)算機(jī)獲得，則h（t）=（1/k ）Ruy（t），這種計(jì)算方法最為常見(jiàn)[15]。

三、河北省農(nóng)業(yè)凈碳匯測(cè)算及影響因素分析

（一） 河北省碳排放測(cè)算結(jié)果

根據(jù)河北省農(nóng)業(yè)碳排放的測(cè)算結(jié)果（表3）可以看出：2000—2019 年間整體趨勢(shì)表現(xiàn)為先增后減。在2000—2007 年期間有小幅增長(zhǎng)，于2007 年達(dá)到最大值，這是化肥、農(nóng)藥和農(nóng)膜使用的碳排放的持續(xù)增加導(dǎo)致的，同時(shí)農(nóng)用柴油使用和農(nóng)作物翻耕引起的碳排放均在2007 年達(dá)到峰值。在2007 年后碳排放量越來(lái)越小，呈平穩(wěn)下降態(tài)勢(shì)。這主要是因?yàn)檗r(nóng)業(yè)的碳排放污染逐漸引起了重視，河北省在2007 年開(kāi)始主動(dòng)采取措施，控制農(nóng)業(yè)的碳排放。

表3 河北省農(nóng)業(yè)碳排放總量（萬(wàn)噸）

分別來(lái)看，農(nóng)用化肥、農(nóng)業(yè)翻耕及農(nóng)業(yè)柴油產(chǎn)生的碳排放占比較大，碳排放總量基本與這三種占比較大的碳源增減變化一致，且整體來(lái)看，河北省的碳排放總量受農(nóng)用柴油的影響最大，因其他兩個(gè)占比較大的碳源在20 年間數(shù)值變化整體較小，而農(nóng)用柴油使用所產(chǎn)生的碳排放在20 年間變化較大，影響著碳排放總量。從降低碳排放總量的角度看，可以控制化肥和農(nóng)用柴油的使用量，調(diào)整農(nóng)業(yè)翻耕面積，以進(jìn)一步減少碳排放。

（二） 河北省農(nóng)業(yè)碳匯測(cè)算結(jié)果

由表4 結(jié)果可知，河北省的農(nóng)業(yè)碳匯量在2000—2019 年間呈現(xiàn)穩(wěn)步增長(zhǎng)趨勢(shì)，20 年間增長(zhǎng)了約38%，這主要得益于現(xiàn)代農(nóng)業(yè)生產(chǎn)技術(shù)與水平的整體提高、全國(guó)各地環(huán)境保護(hù)意識(shí)的加強(qiáng)以及河北省多年來(lái)一直重點(diǎn)關(guān)注著地區(qū)內(nèi)的生態(tài)環(huán)境保護(hù)以及發(fā)展綠色農(nóng)業(yè)。

表4 河北省農(nóng)作物碳匯量（萬(wàn)噸）

研究河北省農(nóng)業(yè)碳匯量構(gòu)成可知，小麥和玉米是碳匯貢獻(xiàn)度最高的作物，主要原因在于其在河北省的播種面積較大，經(jīng)濟(jì)產(chǎn)量較大，長(zhǎng)期來(lái)看這兩種農(nóng)作物在以后仍是河北省碳匯量的主要貢獻(xiàn)來(lái)源。此外，河北省的蔬菜產(chǎn)量雖是農(nóng)作物中最大的，但是蔬菜的碳匯量卻不高，原因?yàn)槭卟俗魑锏暮枯^大，在計(jì)算碳匯過(guò)程中含水量越大實(shí)際產(chǎn)生的碳匯量就越小。

通過(guò)觀察20 年間各農(nóng)作物的經(jīng)濟(jì)產(chǎn)量以及碳匯量可以看出，除了小麥、玉米、薯類(lèi)、甜菜的產(chǎn)量整體呈增長(zhǎng)趨勢(shì)外，其他農(nóng)作物產(chǎn)量均有所減少。出現(xiàn)這種情況可能是市場(chǎng)經(jīng)濟(jì)作用的結(jié)果，農(nóng)業(yè)的生產(chǎn)結(jié)構(gòu)開(kāi)始改變，某些作物的產(chǎn)量在某一年后開(kāi)始降低，從而導(dǎo)致碳匯量的降低。但反觀河北省的碳匯總量，在這期間卻一直穩(wěn)步增長(zhǎng)。綜合來(lái)看，部分作物減產(chǎn)并未導(dǎo)致河北省總碳匯的降低，農(nóng)業(yè)總體產(chǎn)量的增長(zhǎng)和種植水平的提高，使得農(nóng)產(chǎn)品在單位面積上可以實(shí)現(xiàn)更高的產(chǎn)量。以上可以說(shuō)明：對(duì)于河北省而言，更新農(nóng)業(yè)技術(shù)、提高種植水平是發(fā)展碳匯的有效途徑。

（三） 河北省農(nóng)業(yè)凈碳匯測(cè)算結(jié)果

由表5 與圖1 可知，近20 年來(lái)河北省的農(nóng)業(yè)凈碳匯整體呈穩(wěn)步上升態(tài)勢(shì)，有輕微浮動(dòng)。結(jié)合上文我們對(duì)近20 年河北省農(nóng)業(yè)碳排放和農(nóng)業(yè)碳匯的分析來(lái)看，農(nóng)業(yè)凈碳匯的波動(dòng)上升趨勢(shì)的主要原因是：農(nóng)業(yè)碳匯穩(wěn)步增長(zhǎng)；相較于碳匯而言，農(nóng)業(yè)碳排放的數(shù)值浮動(dòng)較小。這表明，農(nóng)業(yè)碳匯的增長(zhǎng)是農(nóng)業(yè)凈碳匯增長(zhǎng)的最主要原因。

圖1 河北省農(nóng)業(yè)碳排放、碳匯量、凈碳匯量對(duì)比圖（萬(wàn)噸）

表5 河北省農(nóng)業(yè)凈碳匯量（萬(wàn)噸）

在過(guò)去的20 年間，由于種植科技的進(jìn)步以及耕地單位畝產(chǎn)量的增加，農(nóng)業(yè)碳匯穩(wěn)步提高；同時(shí)，隨著政府和農(nóng)業(yè)工作者減排意識(shí)的提高，耕種形式越來(lái)越綠色環(huán)保，農(nóng)業(yè)碳排放逐漸在最近幾年逐漸降低，因此河北省的農(nóng)業(yè)凈碳匯呈現(xiàn)出波動(dòng)上升的趨勢(shì)。

（四） VECM 動(dòng)態(tài)分析

本文根據(jù)河北省的實(shí)際情況，選取2000—2019 年的相關(guān)數(shù)據(jù)，構(gòu)建了基于協(xié)整約束的向量誤差修正模型，即VECM 模型，對(duì)河北省農(nóng)業(yè)凈碳匯量的影響因素進(jìn)行研究。從全面性的角度出發(fā)，本文選取了外部因素和內(nèi)部因素來(lái)進(jìn)行分析，外部因素有：河北省農(nóng)業(yè)生產(chǎn)總值（NYGDP）、河北省農(nóng)產(chǎn)品總出口（ZCK）；內(nèi)部因素有：河北省農(nóng)業(yè)機(jī)械總動(dòng)力（JX）、河北省除澇面積（CL）、河北省水土流失治理面積（STLS）、河北省受災(zāi)面積（SZMJ）、河北省水資源總量（SZY）、河北省農(nóng)業(yè)用水量（YSL）以及河北省農(nóng)業(yè)凈碳匯（JTH）。

1. 平穩(wěn)性檢驗(yàn)

本文利用ADF 檢驗(yàn)方法對(duì)影響因素時(shí)間序列的平穩(wěn)性進(jìn)行檢驗(yàn)。本文利用EViews 軟件進(jìn)行計(jì)算，計(jì)算結(jié)果為：農(nóng)業(yè)生產(chǎn)總值、農(nóng)產(chǎn)品總出口、農(nóng)業(yè)機(jī)械總動(dòng)力、除澇面積、水土流失治理面積及農(nóng)業(yè)凈碳匯均為一階差分平穩(wěn)序列；而受災(zāi)面積、水資源總量、農(nóng)業(yè)用水量等因素的原序列及一階差分序列均平穩(wěn)。因此，各序列的一階差分序列均平穩(wěn)，符合協(xié)整檢驗(yàn)的前提要求。

2. 協(xié)整檢驗(yàn)

本文使用Johansen 協(xié)整檢驗(yàn)分析變量間的協(xié)整關(guān)系。要先確定最優(yōu)滯后階數(shù)，利用EViews 軟件進(jìn)行計(jì)算，由AIC 和SC 最小準(zhǔn)則和LR 檢驗(yàn)，確定本文的最優(yōu)滯后階數(shù)為1 階。

由于影響因素的時(shí)間序列長(zhǎng)度較短，難以對(duì)多變量同時(shí)檢驗(yàn)，檢驗(yàn)結(jié)果可能會(huì)出現(xiàn)不準(zhǔn)確現(xiàn)象，因此本文使用Johansen 協(xié)整檢驗(yàn)分別對(duì)河北省農(nóng)業(yè)凈碳匯的影響因素進(jìn)行逐個(gè)檢驗(yàn)，檢驗(yàn)的具體結(jié)果如表6 所示。

表6 協(xié)整檢驗(yàn)結(jié)果

本文以農(nóng)業(yè)生產(chǎn)總值（NYGDP）和農(nóng)業(yè)凈碳匯之間的協(xié)整檢驗(yàn)結(jié)果為例：跡統(tǒng)計(jì)量檢驗(yàn)的結(jié)果顯示，第一列的“None”表示檢驗(yàn)原假設(shè)是：“存在零個(gè)協(xié)整關(guān)系”。該假設(shè)下的跡統(tǒng)計(jì)量大于臨界值，表明至少存在一個(gè)協(xié)整關(guān)系。按照該邏輯，繼續(xù)考察，“Atmost1”表示“至多存在1 個(gè)協(xié)整關(guān)系”的原假設(shè)，該假設(shè)下的跡統(tǒng)計(jì)量大于5%的臨界值，表明在5%的置信水平下，農(nóng)業(yè)生產(chǎn)總值與凈碳匯間存在兩個(gè)協(xié)整關(guān)系。以此類(lèi)推，農(nóng)產(chǎn)品總出口、農(nóng)業(yè)機(jī)械總動(dòng)力、除澇面積、水土流失治理面積、受災(zāi)面積、水資源總量、農(nóng)業(yè)用水量與農(nóng)業(yè)凈碳匯之間均存在著2 個(gè)協(xié)整關(guān)系，均存在長(zhǎng)期均衡關(guān)系，可進(jìn)一步建立VECM 模型。

3. 格蘭杰因果檢驗(yàn)

格蘭杰因果關(guān)系檢驗(yàn)（Granger causality test）是一種假設(shè)檢定的統(tǒng)計(jì)方法，檢驗(yàn)一組時(shí)間序列x 是否為另一組時(shí)間序列y 的原因。本研究采用Eviews 10 對(duì)各個(gè)影響因素與凈碳匯量之間的格蘭杰關(guān)系進(jìn)行檢驗(yàn)，表7 為本次格蘭杰因果檢驗(yàn)的結(jié)果。

表7 格蘭杰因果檢驗(yàn)結(jié)果

通過(guò)表7，在5%的置信水平下可以得出：（1）河北省的農(nóng)業(yè)總產(chǎn)值與農(nóng)業(yè)凈碳匯之間存在著雙向的格蘭杰因果關(guān)系，較為顯著；（2）河北省的農(nóng)業(yè)總出口、除澇面積、水土流失治理面積、水資源總量、農(nóng)業(yè)用水量與凈碳匯間均存在著單向的格蘭杰因果關(guān)系，因此這些因素均會(huì)對(duì)凈碳匯造成一定程度的影響。（3）河北省的農(nóng)業(yè)機(jī)械總動(dòng)力、水土流失治理面積與凈碳匯間互相不存在較為顯著的格蘭杰因果關(guān)系。

4. 脈沖響應(yīng)

根據(jù)格蘭杰因果檢驗(yàn)結(jié)果，本節(jié)對(duì)具備雙向或單項(xiàng)格蘭杰因果關(guān)系的指標(biāo)進(jìn)行脈沖響應(yīng)分析，分析這些變量對(duì)于河北省農(nóng)業(yè)凈碳匯的沖擊。我們將總的滯后期數(shù)確定為10 期，外部因素脈沖響應(yīng)結(jié)果如圖2 所示，內(nèi)部因素脈沖響應(yīng)結(jié)果如圖3 所示。

圖2 外部影響因素對(duì)于凈碳匯變動(dòng)的脈沖響應(yīng)

圖3 內(nèi)部影響因素對(duì)于農(nóng)業(yè)凈碳匯變動(dòng)的脈沖響應(yīng)

在脈沖響應(yīng)分析圖中，橫軸代表追蹤期數(shù)，縱軸代表響應(yīng)水平。下面結(jié)合圖2，對(duì)外部因素的脈沖響應(yīng)進(jìn)行分析。

農(nóng)業(yè)生產(chǎn)總值：當(dāng)農(nóng)業(yè)生產(chǎn)總值對(duì)農(nóng)業(yè)凈碳匯作出正向沖擊時(shí)，凈碳匯表現(xiàn)為正向響應(yīng)，效果明顯。從第3 期后開(kāi)始平穩(wěn)下降，在第10 期時(shí)逐漸趨向于零，這說(shuō)明河北省農(nóng)業(yè)生產(chǎn)總值的變動(dòng)在短期內(nèi)會(huì)使得河北省農(nóng)業(yè)凈碳匯呈正向響應(yīng)且增長(zhǎng)較為明顯，長(zhǎng)期來(lái)看，其影響能力會(huì)逐漸減弱。

農(nóng)業(yè)總出口：當(dāng)給定1 單位沖擊，農(nóng)業(yè)凈碳匯先是表現(xiàn)出正向響應(yīng)，于第2 期開(kāi)始呈現(xiàn)負(fù)向響應(yīng)并到達(dá)最小值，后開(kāi)始平穩(wěn)小幅度增長(zhǎng)，在第3 期達(dá)到正值，后延續(xù)到第6 期開(kāi)始繼續(xù)呈現(xiàn)負(fù)向響應(yīng)，響應(yīng)效果基本維持穩(wěn)定。這說(shuō)明河北省農(nóng)產(chǎn)品總出口的變動(dòng)在短期內(nèi)會(huì)使得河北省農(nóng)業(yè)凈碳匯呈增長(zhǎng)態(tài)勢(shì)，長(zhǎng)期來(lái)看，農(nóng)產(chǎn)品總出口的增加可能會(huì)使凈碳匯減少。

結(jié)合圖3，對(duì)內(nèi)部因素的脈沖響應(yīng)進(jìn)行分析：

農(nóng)業(yè)機(jī)械總動(dòng)力：給定1 單位沖擊，凈碳匯先是負(fù)向波動(dòng)，從第2 期開(kāi)始變呈現(xiàn)正向波動(dòng)，在第3 期達(dá)到峰值，此后一直呈現(xiàn)長(zhǎng)期穩(wěn)定的正向影響，表明機(jī)械總動(dòng)力的投入會(huì)使農(nóng)業(yè)凈碳匯有明顯增加，其對(duì)凈碳匯的發(fā)展起著較為明顯的推動(dòng)作用。

水土流失治理面積：此因素對(duì)于凈碳匯的影響雖在第1 期有短期負(fù)相關(guān)，從第2 期轉(zhuǎn)變?yōu)檎蝽憫?yīng)，且長(zhǎng)期處于在穩(wěn)步上升，這表明水土流失治理面積對(duì)于農(nóng)業(yè)凈碳匯的增加具有顯著的積極影響，響應(yīng)程度相對(duì)較大。

受災(zāi)面積：給定1 單位的沖擊，帶給凈碳匯的影響是正向波動(dòng)的且逐漸增強(qiáng)，整體上呈上升態(tài)勢(shì)，由于河北省的受災(zāi)面積是隨著年份逐漸減小的，故說(shuō)明隨著受災(zāi)面積的減少最終會(huì)使得農(nóng)業(yè)凈碳匯量逐漸增加。

除澇面積：給定1 單位的沖擊，凈碳匯先是表現(xiàn)出正向響應(yīng)，在第2 期時(shí)負(fù)向響應(yīng)下降至了最小值，從第3 期開(kāi)始正向波動(dòng)，后期有小幅度波動(dòng)，整體呈現(xiàn)正向響應(yīng)，長(zhǎng)期來(lái)看最終趨向于零。這表明，短期內(nèi)除澇面積的增加可能會(huì)造成農(nóng)業(yè)凈碳匯短暫的降低，但長(zhǎng)期來(lái)看，這一影響因素有利于農(nóng)業(yè)凈碳匯的發(fā)展且作用明顯。

水資源總量：此因素的沖擊使凈碳匯先是正向波動(dòng)，在第2 期達(dá)到峰值，后影響能力開(kāi)始減緩，在第3 期開(kāi)始呈現(xiàn)負(fù)向響應(yīng)且達(dá)到最小值，后開(kāi)始較為平穩(wěn)的小幅正向波動(dòng)，長(zhǎng)期來(lái)看最終逐漸趨向于零。河北省水資源總量的增加短期內(nèi)會(huì)帶來(lái)農(nóng)業(yè)凈碳匯的增加，長(zhǎng)期呈負(fù)向影響。

農(nóng)業(yè)用水量：給定1 單位沖擊，凈碳匯先負(fù)向響應(yīng)，在第2 期開(kāi)始呈現(xiàn)正向響應(yīng)并達(dá)到最大值，但從第4 期開(kāi)始繼續(xù)呈現(xiàn)反向波動(dòng)，整體來(lái)看，農(nóng)業(yè)用水量對(duì)于農(nóng)業(yè)凈碳匯的沖擊呈現(xiàn)負(fù)向響應(yīng)。農(nóng)業(yè)用水量的增加短期內(nèi)會(huì)使得農(nóng)業(yè)凈碳匯增加，但長(zhǎng)期可能會(huì)造成農(nóng)業(yè)凈碳匯量的下降。

四、結(jié)論與建議

（一） 河北省要進(jìn)一步促進(jìn)農(nóng)業(yè)經(jīng)濟(jì)的發(fā)展

通過(guò)本文研究結(jié)果可知，河北省的農(nóng)業(yè)生產(chǎn)總值（NYGDP）對(duì)于農(nóng)業(yè)凈碳匯有著較為顯著的積極影響，其發(fā)展能夠促進(jìn)凈碳匯量的增加。故河北省應(yīng)大力發(fā)展農(nóng)業(yè)經(jīng)濟(jì)，進(jìn)一步提高農(nóng)業(yè)生產(chǎn)總值。政府應(yīng)起到正向的引導(dǎo)作用，首先要加強(qiáng)農(nóng)業(yè)發(fā)展的政策支撐，對(duì)于省內(nèi)的農(nóng)業(yè)如何發(fā)展、發(fā)展速度等做進(jìn)一步的詳細(xì)指示；其次要發(fā)展科技農(nóng)業(yè)、綠色農(nóng)業(yè)、品牌農(nóng)業(yè)、質(zhì)量農(nóng)業(yè)，推動(dòng)全省的農(nóng)業(yè)現(xiàn)代化；最后，持續(xù)優(yōu)化農(nóng)業(yè)生產(chǎn)體系和經(jīng)營(yíng)體系，進(jìn)一步推動(dòng)河北省農(nóng)業(yè)發(fā)展質(zhì)量提升，進(jìn)而促進(jìn)河北省農(nóng)業(yè)生產(chǎn)總值的提高。

（二） 適當(dāng)加強(qiáng)農(nóng)產(chǎn)品的對(duì)外出口

目前，河北省農(nóng)產(chǎn)品在出口貿(mào)易方面還相對(duì)薄弱，農(nóng)產(chǎn)品出口總量相對(duì)較小，且增長(zhǎng)速度較為緩慢。根據(jù)上文的研究可知，農(nóng)產(chǎn)品總出口量（ZCK）會(huì)影響農(nóng)業(yè)凈碳匯的發(fā)展，短期內(nèi)，它對(duì)凈碳匯有著較強(qiáng)的推動(dòng)作用，故合理優(yōu)化農(nóng)產(chǎn)品的產(chǎn)量與質(zhì)量，進(jìn)一步加強(qiáng)對(duì)外出口，這對(duì)于河北省的農(nóng)業(yè)增匯有著顯著的積極影響。建議河北省政府進(jìn)一步出臺(tái)相關(guān)政策，鼓勵(lì)外向型企業(yè)的發(fā)展，從政策與資金角度扶持一些外貿(mào)的龍頭企業(yè)，使其帶動(dòng)省內(nèi)的出口貿(mào)易。河北省還應(yīng)注重農(nóng)產(chǎn)品品牌的創(chuàng)立與注冊(cè)，打造農(nóng)產(chǎn)品知名國(guó)際品牌，如保定醬菜、藁城強(qiáng)筋小麥、新樂(lè)花生、無(wú)極黃瓜等等。此外，還要加強(qiáng)對(duì)于農(nóng)產(chǎn)品的健康安全監(jiān)測(cè)，保證出口產(chǎn)品的健康安全。最后，河北省應(yīng)拓寬農(nóng)產(chǎn)品出口市場(chǎng)，在穩(wěn)定日本、韓國(guó)、印尼、美國(guó)、德國(guó)這五大出口市場(chǎng)的前提下，不斷發(fā)掘新市場(chǎng)。

（三） 不斷提高農(nóng)業(yè)種植科技水平

通過(guò)研究可知，農(nóng)業(yè)機(jī)械總動(dòng)力（JX）對(duì)于農(nóng)業(yè)凈碳匯的發(fā)展具有較為明顯的推動(dòng)作用，投入機(jī)械動(dòng)力越大，凈碳匯量會(huì)呈增長(zhǎng)態(tài)勢(shì)，故農(nóng)業(yè)機(jī)械化水平和整個(gè)種植業(yè)科技水平的提高，對(duì)于發(fā)展農(nóng)業(yè)凈碳匯有著積極意義。河北省應(yīng)樹(shù)立現(xiàn)代農(nóng)業(yè)發(fā)展的新理念，不斷提高農(nóng)業(yè)的科技水平，加大科技投入，促進(jìn)農(nóng)業(yè)的健康穩(wěn)定發(fā)展，在此基礎(chǔ)上繼續(xù)減碳增匯。

（四） 增強(qiáng)水土流失治理、減少受災(zāi)面積

研究表明，水土流失治理面積（STLS）的增加和耕地受災(zāi)面積（SZMJ）的減少能夠有助于農(nóng)業(yè)凈碳匯量的穩(wěn)步增長(zhǎng)。本文認(rèn)為在水土流失治理和減少受災(zāi)面積方面，一是要加強(qiáng)宣傳教育，河北省可以通過(guò)對(duì)于水土保持等法律法規(guī)的宣傳，增強(qiáng)社會(huì)各界的防患意識(shí)；二是河北省政府應(yīng)針對(duì)省內(nèi)實(shí)際情況，制定合理的相關(guān)法律法規(guī)，以政策調(diào)動(dòng)人民參與治理與保護(hù)的積極性；三是在減少省內(nèi)受災(zāi)面積的預(yù)防措施上，要繼續(xù)堅(jiān)持防風(fēng)防洪工程，并進(jìn)一步健全各種災(zāi)害的監(jiān)測(cè)預(yù)報(bào)和信息管理系統(tǒng)。

（五） 合理使用水資源

通過(guò)前文研究可知，河北省的水資源總量（SZY）和農(nóng)業(yè)用水量（YSL）在短期內(nèi)都會(huì)對(duì)農(nóng)業(yè)凈碳匯量產(chǎn)生積極的推動(dòng)作用，但長(zhǎng)期來(lái)看，過(guò)多地投入或不合理使用可能會(huì)使凈碳匯呈現(xiàn)出負(fù)增長(zhǎng)的趨勢(shì)；且河北省自身的水資源較為匱乏，多年來(lái)一直通過(guò)南水北調(diào)工程來(lái)彌補(bǔ)自身水資源不足的情況，截至2021 年12 月，南水北調(diào)工程已為河北省累計(jì)分水超過(guò)130 億立方米。故基于此前提，本文建議河北省要注重合理使用水資源，合理投入農(nóng)業(yè)用水量，切不可一味地過(guò)多投入使用，而是要在節(jié)約用水和減少水資源污染的前提下穩(wěn)步提高農(nóng)業(yè)的有效灌溉水平。此外，河北省還要繼續(xù)合理地、科學(xué)地利用北調(diào)水，建議加強(qiáng)監(jiān)測(cè)與管理，減少水資源浪費(fèi)，發(fā)揮出水資源的最大作用，進(jìn)一步獲得更大的經(jīng)濟(jì)、社會(huì)以及生態(tài)效益。