天然林保護工程對中國林業(yè)全要素生產(chǎn)率的影響
——基于省級面板數(shù)據(jù)實證檢驗

伍海泉，印 慧，曹玉昆

（東北林業(yè)大學經(jīng)濟管理學院，哈爾濱 150040）

據(jù)第九次全國森林資源清查（2014～2018年）結果，我國森林面積為2.2 億hm2，森林覆蓋率為22.96%，森林蓄積量為175.6 億m3，實現(xiàn)30年連續(xù)保持森林面積、森林蓄積量“雙增長”[1]。大氣二氧化碳（CO2）濃度升高是影響陸地生態(tài)系統(tǒng)可持續(xù)性的主要全球變化因子之一[2]，而森林則是陸地生態(tài)系統(tǒng)中重要的碳氮儲存庫。應繼續(xù)加大力度推進生態(tài)文明建設，堅持節(jié)約自然資源與保護生態(tài)環(huán)境。20 世紀80年代以來，我國提出并實施多項生態(tài)計劃，包括退耕還林工程、長江和珠江流域及農(nóng)田防護林體系建設工程、天然林保護工程等。1998年長江流域和松花江、嫩江流域特大洪災后，黨中央、國務院決定在云南等12 個省（區(qū)、市）重點國有林區(qū)開展天然林保護工程試點。2000年《長江上游、黃河上中游地區(qū)天然林資源保護工程實施方案》經(jīng)國務院批準實施，《東北、內(nèi)蒙古等重點國有林區(qū)天然林資源保護工程實施方案》隨即出臺。至此，天然林保護工程在我國17 個?。▍^(qū)、市）重點林區(qū)實施，具體包括長江上游地區(qū)云南、四川、貴州、重慶、湖北、西藏6 ?。▍^(qū)、市），黃河上中游地區(qū)陜西、甘肅、青海、寧夏、內(nèi)蒙古、山西、河南7 ?。▍^(qū)），以及黑龍江、吉林、海南、新疆4?。▍^(qū)）。天然林保護工程實施是針對我國長期以來天然林資源過度消耗引起生態(tài)環(huán)境惡化，國家從社會經(jīng)濟可持續(xù)發(fā)展戰(zhàn)略高度做出的重大決策[3]。天然林保護工程作為全國最大林業(yè)生態(tài)保護工程之一，實施20年來，中央財政總投入已超過4 000 億元[4]，林區(qū)自然資源與社會面貌發(fā)生巨大變化。

目前，天然林保護一期工程已結束，二期工程進入關鍵時段。如何客觀評價天然林保護工程實施效果，哪些因素對其影響顯著，是該項生態(tài)保護工程尚待解決的問題。以往研究主要集中在以下三個方面：一是天然林保護工程實踐綜述，研究內(nèi)容包括天然林保護工程實施手段和取得成績。陳麗榮等通過繪制可視化知識圖譜，展示天然林保護工程研究領域時間分布和研究熱點演進，揭示我國天然林保護工程研究整體脈絡[5]；曹玉昆等通過對我國天然林保護工程效益評價研究文獻梳理，系統(tǒng)分析天然林保護工程效益評價主要特點和方法[6]。二是天然林保護工程績效和價值研究，例如秦利等通過構建天然林保護工程政策實施績效綜合評價指標體系，采用模糊綜合評價法對伊春林區(qū)天保工程二期政策實施績效進行綜合評價[7]；崔方茹等利用靖邊縣天保工程相關資料和既有成果，構建森林生態(tài)系統(tǒng)服務功能評估指標體系，評估靖邊縣天保工程森林生態(tài)系統(tǒng)服務功能價值[8]。三是天保工程相關影響研究。如臧良震運用西部11 個省級面板數(shù)據(jù)，采用面板向量自回歸和脈沖響應分析研究工程建設對農(nóng)民收入產(chǎn)生影響[9]；翟緒軍等采用黑龍江省國有林區(qū)林業(yè)產(chǎn)業(yè)數(shù)據(jù)，分析實施天然林保護工程對國有林區(qū)林業(yè)產(chǎn)業(yè)總體結構變動產(chǎn)生的影響[10]。已有報道多是通過文獻歸納整理對天然林保護工程定性評價分析為主，缺乏以計量核算為主的定量評價，且研究對象或面板數(shù)據(jù)多集中于西部、東北部等區(qū)域天然林保護工程，針對全國天然林保護工程政策實施效果評價尤其與林業(yè)全要素生產(chǎn)率相結合研究較少。因此，本文使用多期雙重差分法（Difference in difference，DID），運用我國1997～2018年省級面板數(shù)據(jù)，基于林業(yè)全要素生產(chǎn)率視角，天然林保護工程項目實施效果作出評價，為我國森林資源經(jīng)營管理和生態(tài)資源保護提供參考。

1 模型設定

1.1 方法選取與模型設定

為研究天然林保護工程對林業(yè)全要素生產(chǎn)率政策效應，本文選取雙重差分法（DID）對政策實施地區(qū)作驗證。雙重差分模型適用基礎是先將觀測樣本準確分為試驗組和對照組，本文選擇實施天然林保護工程省份為“試驗組”，其他省份作為相應“對照組”，設置DID回歸方程如下：

上述公式被解釋變量中，設置政策影響變量為yit；觀測組別為d i，分別代表試驗組和對照組；觀測年份為dt；μit是殘差項。在借鑒相關研究前提下，本文在原基礎模型中加入部分控制變量。拓展后DID方程如下:

方程中，下標p代表省份，t代表觀測年份。在Treatpt賦值中，甘肅、貴州、海南、河南、黑龍江、湖北、吉林、內(nèi)蒙古、寧夏、青海、山西、陜西、四川、西藏、新疆、云南、重慶地區(qū)為試驗組1，其他地區(qū)為對照組0。Yearpt為天保工程實施年份，1998年及以后賦值為1，其他為0；Xpt為控制變量；γp為地區(qū)固定效應；ηt為年份固定效應；μpt為隨機擾動項。β1為重點關注系數(shù)，Treatpt×Yearpt則是反映政策凈效應交乘變量。其次，參數(shù)β2表示無天然林保護工程政策試驗組和對照組之間目標變量時間變化情況，參數(shù)β3表示試驗組與對照組之間不隨時間變動差異，式（2）中包含γp地區(qū)虛擬變量和ηt年份虛擬變量。

1.2 指標選取與數(shù)據(jù)來源

本文被解釋變量為通過全要素生產(chǎn)率動態(tài)變化的Malmquist指數(shù)DEA模型測度我國31個省級林業(yè)全要素生產(chǎn)率，Malmquist指數(shù)以t期為基期，公式如下：

上述公式中，MALM為從t到t+1期林業(yè)全要素生產(chǎn)率、PECH為從t到t+1 期純技術效率變動指數(shù)、SECH為從t到t+1 期規(guī)模效率變動指數(shù)、TECHCH為從t到t+1 期技術進步變動指數(shù)、EF?FCH為從t到t+1 期技術效率變動指數(shù)。xt和xt+1分別是t到t+1期單位投入向量，yt和yt+1分別是t到t+1 期單位產(chǎn)出向量，Dt(xt,yt)、Dt+1(xt+1,yt+1)分別是以t和t+1 時期技術為參照距離函數(shù)。上述公式中下標v是可變規(guī)模報酬假設，c 是不變規(guī)模報酬假設，計算結果如表1所示。

表1 1997～2018年我國31個省級林業(yè)全要素生產(chǎn)率Malmquist指數(shù)及分解Table1 Malmquist index and decomposition of average total factor productivity of forestry in 31 provinces in China from 1997 to 2018

本文參考相關研究成果后選取系列控制變量是在兼顧數(shù)據(jù)指標可得性前提下，基于天然林保護工程實施制度、管理和環(huán)境3個維度，確定經(jīng)濟發(fā)展水平、城鎮(zhèn)化水平、對外開放水平、科技發(fā)展水平、產(chǎn)業(yè)現(xiàn)代化程度、政府干預程度、自然災害情況、林業(yè)人力資本水平等變量。確定變量值為經(jīng)計算處理比率值，分別采用地區(qū)人均GDP對數(shù)代表經(jīng)濟發(fā)展水平，常住城鎮(zhèn)人口占總人口比重的城鎮(zhèn)化率衡量城鎮(zhèn)化水平，進出口貿(mào)易額比重代表對外開放程度，研究與試驗發(fā)展經(jīng)費投入代表科技發(fā)展水平，林業(yè)第二、三產(chǎn)業(yè)占林業(yè)產(chǎn)業(yè)總產(chǎn)值比重代表產(chǎn)業(yè)結構，中央政府投資占林業(yè)總投資比重代表政府干預程度，森林病蟲害發(fā)生面積代表自然災害情況，在崗職工年平均工資代表林業(yè)人力資本投入水平。數(shù)據(jù)主要來源于1998～2018年林業(yè)統(tǒng)計年鑒和國家統(tǒng)計局數(shù)據(jù)庫。本文選擇1998～2018年我國31 個省級地區(qū)，合計5 859 個樣本開展經(jīng)驗驗證，具體變量描述性統(tǒng)計如表2所示。

表2 1997～2018年樣本描述性統(tǒng)計Table 2 Descriptive statistics of samples from 1997 to 2018

續(xù)表

2 實證結果

2.1 實證結果與差分分析

本文在模型測算中逐步加入控制變量以研究天然林保護工程政策凈效應，估計結果如表3 所示。在天然林保護工程實施制度、管理和環(huán)境等控制變量逐步加入后，解釋變量系數(shù)符號和顯著性并未發(fā)生本質變化，說明估計結果穩(wěn)健性良好。由于實證結果為0.0866～0.1057，在5%置信水平顯著，說明天然林保護工程實施顯著提高山西、內(nèi)蒙古、吉林等17 個政策實施省級地區(qū)林業(yè)全要素生產(chǎn)效率。在影響天然林保護工程提高生產(chǎn)率控制變量中，產(chǎn)業(yè)結構、人力資本、科技發(fā)展水平和經(jīng)濟發(fā)展水平影響顯著。其中，人力資本投入程度在5%置信水平顯著，產(chǎn)業(yè)結構、科技和經(jīng)濟發(fā)展水平在10%置信水平顯著。分析結果顯示控制變量中人力資本投入程度最為顯著，主要是因隨著科技發(fā)展水平提高和天然林保護工程資金投入力度加大，人均生產(chǎn)效率得到顯著提高，直接推動林業(yè)全生產(chǎn)要素提高；另一方面，林區(qū)經(jīng)濟發(fā)展水平提高帶來工資收入水平上升，人均勞動力成本加大，也促使相關企業(yè)和經(jīng)營主體提高管理水平，實現(xiàn)減員增效。同時，產(chǎn)業(yè)結構、科技和經(jīng)濟發(fā)展水平提升一定程度上對沖禁止商業(yè)性砍伐對林業(yè)企業(yè)收入影響，開辟新財源，為提升林業(yè)全要素生產(chǎn)率提供資金支持、物質基礎和科技保障。與人力資本、產(chǎn)業(yè)結構、科技發(fā)展水平和經(jīng)濟發(fā)展水平相比，城鎮(zhèn)化水平、對外開放程度、政府干預和自然災害情況在林業(yè)全要素生產(chǎn)效率中影響不顯著。從區(qū)域經(jīng)濟結構特點看，天然林保護工程實施加快第三產(chǎn)業(yè)占比提高，但由于林區(qū)氣候、土地等資源限制以及《森林法》等相關法律法規(guī)規(guī)制，歷史上第一產(chǎn)業(yè)所占比例較低，農(nóng)業(yè)人口比重相對較低。天然林保護工程實施引起農(nóng)業(yè)人口向城鎮(zhèn)人口轉化空間有限，并不足構成顯著影響。同時，改革開發(fā)以來，我國對木材等林業(yè)產(chǎn)品進口數(shù)量逐年提升，但林業(yè)相關資金、勞動力和科技等生產(chǎn)要素跨國界流動并不顯著。對外開放程度變化主要受林區(qū)內(nèi)生經(jīng)濟增長和外部貿(mào)易環(huán)境影響，與林業(yè)全要素生產(chǎn)率相關程度低。天然林保護工程實施區(qū)域大部分為國有林區(qū)，政府干預程度較高，政策實施并未顯著提高政府干預水平。天然林保護工程實施增加政府對林區(qū)發(fā)展資金投入，但資金相當比例被用于消耗性和保障性支出，受林業(yè)管理水平提升影響不顯著。

表3 DID 回歸結果分析Table 3 Regression results

2.2 分組驗證結果與差分分析

為進一步研究天然林保護工程實施過程中資源配置效率和影響因素作用機制，本文將各變量再次分組。首先將我國31 個省域按照1997～2018年年均森林覆蓋率進行排名，然后將森林覆蓋率排名前15且實施天然林保護工程省級地區(qū)，即海南、云南、黑龍江、吉林、陜西、湖北、四川、重慶8個地區(qū)分為試驗組，其余年均森林覆蓋率未排全國前15 且實施天然林保護工程省級地區(qū)，即貴州、內(nèi)蒙古、河南、山西、西藏、寧夏、青海、新疆、甘肅9個地區(qū)為對照組。首先對森林覆蓋率排名前15 名8 個省級地區(qū)即試驗組進行異質性分析（見表4），回歸結果說明，森林資源覆蓋率高與天然林保護工程對該地區(qū)林業(yè)全要素生產(chǎn)率影響有正向作用，政策實施效果影響更為顯著。

表4 異質性分析和森林資源覆蓋率分組回歸結果Table 4 Heterogeneity analysis and grouping regression results of forest resource coverage

式中，minori為虛擬變量，若i省是森林資源排名前15，則minori=1，否則為0。

森林覆蓋率是反映森林資源豐富程度重要參考指標，天然林保護工程對林業(yè)全要素生產(chǎn)率影響與森林覆蓋率顯著正相關，說明政策實施后所在地區(qū)森林蓄積、面積、覆蓋率顯著提高，森林結構趨于合理、森林質量明顯改善[11]。同時天然林保護工程投入巨額資金，顯著加強政策實施省域森林生態(tài)系統(tǒng)、濕地生態(tài)系統(tǒng)、荒漠生態(tài)系統(tǒng)建設和生物多樣性保護，促進該地區(qū)森林資源和林業(yè)產(chǎn)業(yè)協(xié)調(diào)發(fā)展。

3 結論與政策啟示

3.1 結論

本文以我國31個省域林業(yè)全要素生產(chǎn)率（TFP）Malmquist 指數(shù)為被解釋變量，基于天然林保護工程實施政策、管理和環(huán)境等維度，將天然林保護工程這一重點生態(tài)工程實施視為一項政策沖擊，通過DID模型對政策實施省域和未實施省域分組測算分析，引入經(jīng)濟發(fā)展水平、城鎮(zhèn)化水平、對外開放水平、科技發(fā)展水平、產(chǎn)業(yè)現(xiàn)代化程度、政府干預程度、自然災害情況、林業(yè)人力資本水平等其他變量因素。結果表明，天然林保護工程實施整體上對林業(yè)全要素生產(chǎn)率影響顯著。人力投入減少、產(chǎn)業(yè)結構合理化、科技投入和經(jīng)濟發(fā)展水平提高，可有效提升要素資源在林業(yè)生產(chǎn)之間配置效率，城鎮(zhèn)化水平和對外開放程度提高、政府干預和自然災害情況對林業(yè)全要素生產(chǎn)效率影響并不顯著。同時，森林覆蓋率高地區(qū)，天然林保護工程對其林業(yè)全要素生產(chǎn)率提升影響更加顯著。充分說明天然林保護工程政策舉措和資金投入，有效促進政策實施省域森林結構更合理、森林資源更豐富。

3.2 政策啟示

①全面總結天然林保護工程實施經(jīng)驗，適度擴大工程實施范圍。優(yōu)化天保工程資金支出結構，同時積極填補森林撫育資金缺口，適度上調(diào)森林撫育補貼標準，通過擴大森林撫育補貼規(guī)模和范圍，逐步提高森林資源質量加快低質森林資源改造速度。鼓勵開發(fā)林業(yè)生態(tài)產(chǎn)品，走綠色發(fā)展之路，加快經(jīng)濟發(fā)展[12]，保護林業(yè)生態(tài)多樣性，開展“山水林田湖草”整體開發(fā)，提升林區(qū)生態(tài)產(chǎn)品開發(fā)利用效率。

②積極培養(yǎng)壯大林業(yè)科研隊伍，加速林業(yè)科技成果轉化。建設培養(yǎng)專業(yè)化的林業(yè)技術應用和推廣團隊，應用現(xiàn)代化信息技術和媒體渠道，幫助基層技術人員掌握多種技術和裝備。科學編制人才培養(yǎng)計劃，創(chuàng)新人才使用機制，培養(yǎng)高素質林業(yè)工程管理人員、高層次科研人員、高水平技術人員、高技能產(chǎn)業(yè)工人。完善營林管護技術手段，支持利用衛(wèi)星遙感探測、無人機巡航等技術手段，擴大人均管護面積，提高人均工作效率，全面保護天然林。

③加強森林資源生態(tài)保護，提升地區(qū)森林覆蓋率。要加大管護和執(zhí)法力度，采取有效手段保護天然林；加強人工林封育管理，加強退化林分改造，提高有林地郁閉度和生物多樣性。同時加大國有林區(qū)基礎設施建設支持力度，將電網(wǎng)與飲水安全以及管護站點用房等基礎設施建設納入上級政府建設規(guī)劃統(tǒng)籌安排，同時將國有林區(qū)道路按屬性納入相關公路網(wǎng)規(guī)劃。依法依規(guī)嚴管森林資源，將非林業(yè)部門管理林地納入森林資源統(tǒng)計管理，加強林地監(jiān)管，確保林地面積只增不減。