中國(guó)林業(yè)綠色全要素生產(chǎn)率的時(shí)空演變分析

呂潔華，孫嘉宇，蔡秀亭

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

一、引言與文獻(xiàn)綜述

生態(tài)文明建設(shè)是我國(guó)全面建成社會(huì)主義現(xiàn)代化強(qiáng)國(guó)的重要一步，是關(guān)系經(jīng)濟(jì)社會(huì)、人民福祉和民族未來(lái)的長(zhǎng)遠(yuǎn)大計(jì)。林業(yè)產(chǎn)業(yè)作為綠色產(chǎn)業(yè)，在完善傳統(tǒng)生產(chǎn)結(jié)構(gòu)單一模式后，為經(jīng)濟(jì)建設(shè)和生態(tài)建設(shè)作出重要貢獻(xiàn)。然而，長(zhǎng)期以來(lái)我國(guó)林業(yè)發(fā)展存在的技術(shù)水平落后、生產(chǎn)效率低下、從業(yè)人員環(huán)保意識(shí)薄弱等問(wèn)題嚴(yán)重制約林業(yè)經(jīng)濟(jì)的持續(xù)穩(wěn)定發(fā)展［1］。黨的十八屆五中全會(huì)以后，“綠色”成為新發(fā)展理念之一，經(jīng)濟(jì)發(fā)展與資源環(huán)境保護(hù)相協(xié)調(diào)成為新的發(fā)展潮流?！笆奈濉币?guī)劃提出：林業(yè)發(fā)展要踐行“綠水青山就是金山銀山”理念，堅(jiān)持高質(zhì)量發(fā)展鮮明導(dǎo)向，走“經(jīng)濟(jì)生態(tài)化、生態(tài)經(jīng)濟(jì)化”現(xiàn)代林業(yè)發(fā)展道路。而我國(guó)林業(yè)發(fā)展在依賴大量投入資源實(shí)現(xiàn)經(jīng)濟(jì)騰飛后步入滯緩期，全要素生產(chǎn)率開始成為林業(yè)實(shí)現(xiàn)長(zhǎng)期增長(zhǎng)的根本，主要包括技術(shù)創(chuàng)新、規(guī)模效應(yīng)等要素投入以外影響效率的動(dòng)力來(lái)源。綠色全要素生產(chǎn)率逐漸成為考察經(jīng)濟(jì)效率的重要指標(biāo)，而對(duì)林業(yè)綠色全要素生產(chǎn)率的全面把握將有助于穩(wěn)定發(fā)展資源節(jié)約型、環(huán)境友好型林業(yè)產(chǎn)業(yè)。

全要素生產(chǎn)率最初表示單一要素生產(chǎn)率，逐步演變?yōu)槎喾N投入要素對(duì)經(jīng)濟(jì)增長(zhǎng)的貢獻(xiàn)率，能更加準(zhǔn)確地衡量經(jīng)濟(jì)活動(dòng)效率［2］；而后基于經(jīng)濟(jì)假設(shè)與實(shí)際的符合程度，全要素生產(chǎn)率具體定義為要素投入對(duì)經(jīng)濟(jì)增長(zhǎng)貢獻(xiàn)以外的其他因素的貢獻(xiàn)量［3］。林業(yè)發(fā)展進(jìn)入要素驅(qū)動(dòng)增長(zhǎng)的瓶頸時(shí)期，全要素生產(chǎn)率成為推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)的重要著力點(diǎn)。史常亮等［4］研究發(fā)現(xiàn)，中國(guó)林業(yè)全要素生產(chǎn)率水平較低，主要由于林業(yè)生產(chǎn)技術(shù)效率低下，技術(shù)無(wú)效率現(xiàn)象普遍存在，也表明林業(yè)技術(shù)進(jìn)步是驅(qū)動(dòng)林業(yè)發(fā)展的主要因素［5］。而陳向華等［6］研究黑龍江省林業(yè)生產(chǎn)效率時(shí)得出不同結(jié)論，即全要素生產(chǎn)率低下主要是由于林區(qū)技術(shù)落后所致。地區(qū)間資源稟賦、經(jīng)濟(jì)實(shí)力、技術(shù)創(chuàng)新水平等不同勢(shì)必導(dǎo)致林業(yè)生產(chǎn)效率的差異。徐瑋等［7］研究發(fā)現(xiàn)我國(guó)中部地區(qū)屬于林業(yè)較發(fā)達(dá)地區(qū)，但未達(dá)到效率最優(yōu)，而其余地區(qū)林業(yè)效率處于最佳狀態(tài)；基于多重效益目標(biāo)，黃安勝等［8］研究發(fā)現(xiàn)，地區(qū)間林業(yè)全要素生產(chǎn)率及其分解水平迥乎不同，主要源于潛在的驅(qū)動(dòng)因素差異。同時(shí)，全國(guó)范圍內(nèi)產(chǎn)業(yè)網(wǎng)絡(luò)的完善，在一定程度上打破產(chǎn)業(yè)間市場(chǎng)壁壘，促使林業(yè)全要素生產(chǎn)率的空間差異降低，呈收斂發(fā)展趨勢(shì)［9］。此外，林業(yè)產(chǎn)業(yè)的空間溢出效應(yīng)以及相似的經(jīng)濟(jì)、資源、政策條件會(huì)推動(dòng)林業(yè)在空間上形成集群，全要素生產(chǎn)率水平呈現(xiàn)集聚分布［10］，而林業(yè)產(chǎn)業(yè)空間集聚可進(jìn)一步影響地區(qū)林業(yè)全要素生產(chǎn)率水平［11-12］。

為推動(dòng)林業(yè)全要素生產(chǎn)率，促進(jìn)林業(yè)經(jīng)濟(jì)發(fā)展，學(xué)者們廣泛搜尋影響林業(yè)全要素生產(chǎn)率的因素，深入挖掘其內(nèi)部作用機(jī)制。曹玉昆等［13］研究發(fā)現(xiàn)，國(guó)家對(duì)林業(yè)的財(cái)政支持能有效推動(dòng)林業(yè)全要素生產(chǎn)率增長(zhǎng)；董婭楠等［14］研究認(rèn)為，林業(yè)全要素生產(chǎn)率受當(dāng)期FDI 水平影響不顯著，而對(duì)于滯后一期的FDI 具有正向反饋，表現(xiàn)出外商直接投資對(duì)林業(yè)產(chǎn)業(yè)的影響具有一定時(shí)滯性；揭昌亮等［15］深入分析進(jìn)出口貿(mào)易對(duì)林業(yè)生產(chǎn)效率的影響，發(fā)現(xiàn)產(chǎn)業(yè)在參與國(guó)際貿(mào)易活動(dòng)過(guò)程中可以改變企業(yè)生產(chǎn)技術(shù)和規(guī)模經(jīng)濟(jì)，進(jìn)而影響全要素生產(chǎn)率。另外，劉東等［16］還探究林業(yè)發(fā)展基礎(chǔ)水平、人力資本水平、造林護(hù)林投入等因素對(duì)林業(yè)全要素生產(chǎn)率的影響。

隨著社會(huì)對(duì)生態(tài)環(huán)境關(guān)注度提高，全要素生產(chǎn)率與生態(tài)環(huán)境的結(jié)合更符合高質(zhì)量發(fā)展要義，綠色全要素生產(chǎn)率也廣泛用于效率測(cè)算［17］。學(xué)者們較早展開對(duì)林業(yè)可持續(xù)發(fā)展的相關(guān)研究，如郎曉娟等［18］從綠色經(jīng)濟(jì)視角下分析發(fā)現(xiàn)，產(chǎn)業(yè)的“綠色投資”“綠色創(chuàng)新”“綠色新政”將會(huì)提升林業(yè)全要素生產(chǎn)率，推動(dòng)產(chǎn)業(yè)高質(zhì)量發(fā)展；劉清泉等［19］從可持續(xù)發(fā)展視角探究廣東省林業(yè)全要素生產(chǎn)率和地區(qū)經(jīng)濟(jì)的協(xié)同發(fā)展關(guān)系。然而，直接將生態(tài)環(huán)境資源要素加入到林業(yè)全要素生產(chǎn)率測(cè)算的研究尚處于起步階段。沈偉航等［20］發(fā)現(xiàn)林產(chǎn)工業(yè)的綠色全要素生產(chǎn)率增速減緩，主要原因是純效率的下降；劉濤等［21］創(chuàng)新性測(cè)算中國(guó)林業(yè)綠色全要素生產(chǎn)率，發(fā)現(xiàn)其在效率和技術(shù)進(jìn)步雙驅(qū)動(dòng)下呈增長(zhǎng)趨勢(shì)，但在空間范圍內(nèi)卻呈現(xiàn)非均衡分布，并進(jìn)一步分析政府支持等因素對(duì)地區(qū)林業(yè)綠色全要素生產(chǎn)率的影響。

綜合來(lái)看，國(guó)內(nèi)外學(xué)者對(duì)林業(yè)全要素生產(chǎn)率的測(cè)算及影響因素研究比較充足，但對(duì)于林業(yè)綠色全要素生產(chǎn)率的研究仍處于起步階段，缺少對(duì)其時(shí)空演變的定量分析以及狀態(tài)轉(zhuǎn)移分析。基于此，本文采用超效率SBM 模型重新測(cè)算我國(guó)林業(yè)綠色全要素生產(chǎn)率，結(jié)合區(qū)域差異性分析及馬爾可夫轉(zhuǎn)移客觀探究其演變歷程并進(jìn)行長(zhǎng)期預(yù)測(cè)，把握基本演變規(guī)律，為推動(dòng)林業(yè)“經(jīng)濟(jì)生態(tài)化，生態(tài)經(jīng)濟(jì)化”發(fā)展提供科學(xué)依據(jù)。

二、理論分析與研究假說(shuō)

20世紀(jì)30年代，柯布和道格拉斯提出的C-D 生產(chǎn)函數(shù)揭示經(jīng)濟(jì)活動(dòng)的要素投入和產(chǎn)出關(guān)系，歸納資本、勞動(dòng)力以及除此以外的其他部分因素為實(shí)現(xiàn)經(jīng)濟(jì)活動(dòng)的源泉。在此基礎(chǔ)上，索羅于20世紀(jì)50年代提出具有規(guī)模報(bào)酬不變特性的總量生產(chǎn)函數(shù)，確定資本及勞動(dòng)力以外的因素為全要素生產(chǎn)率，并將其歸結(jié)為技術(shù)進(jìn)步。20 世紀(jì)80 年代中期，內(nèi)生增長(zhǎng)理論進(jìn)一步確定內(nèi)生的技術(shù)進(jìn)步是實(shí)現(xiàn)經(jīng)濟(jì)長(zhǎng)期增長(zhǎng)的決定因素，也表明全要素生產(chǎn)率的增長(zhǎng)是實(shí)現(xiàn)經(jīng)濟(jì)長(zhǎng)期穩(wěn)定增長(zhǎng)的根本動(dòng)力來(lái)源。而林業(yè)綠色全要素生產(chǎn)率綜合考慮林業(yè)經(jīng)濟(jì)發(fā)展動(dòng)能與生態(tài)資源環(huán)境保護(hù)的思想，是林業(yè)產(chǎn)業(yè)實(shí)現(xiàn)綠色可持續(xù)發(fā)展的重要測(cè)度指標(biāo)。

由于林業(yè)產(chǎn)業(yè)具有高度的資源環(huán)境依賴屬性，同時(shí)各區(qū)域林業(yè)產(chǎn)業(yè)在后期的生產(chǎn)活動(dòng)中自發(fā)產(chǎn)生知識(shí)信息交流，使得林業(yè)綠色全要素生產(chǎn)率在長(zhǎng)期內(nèi)呈現(xiàn)出發(fā)散或收斂態(tài)勢(shì)。因此，假定林業(yè)綠色全要素生產(chǎn)率長(zhǎng)期內(nèi)有區(qū)域內(nèi)和區(qū)域間兩條演變路徑。

地理學(xué)第一定律揭示任何事物在空間范圍內(nèi)都存在關(guān)聯(lián)，且關(guān)聯(lián)度與相鄰距離呈反比關(guān)系，即距離越近，關(guān)聯(lián)程度越高。林業(yè)發(fā)展對(duì)基礎(chǔ)林業(yè)產(chǎn)業(yè)的依賴性較高，即林業(yè)產(chǎn)業(yè)廠商傾向于選擇靠近林業(yè)資源豐富的區(qū)位，由此誘導(dǎo)林業(yè)產(chǎn)業(yè)在資源充裕位置形成集群。同時(shí)，林業(yè)企業(yè)在競(jìng)爭(zhēng)發(fā)展過(guò)程中存在知識(shí)和信息交流、技術(shù)學(xué)習(xí)模仿等行為，使區(qū)域內(nèi)企業(yè)關(guān)聯(lián)度提高，帶動(dòng)區(qū)域內(nèi)企業(yè)經(jīng)濟(jì)運(yùn)作模式、生產(chǎn)效率逐漸趨同。最終林業(yè)發(fā)展表現(xiàn)為，地區(qū)林業(yè)全要素生產(chǎn)率對(duì)周邊區(qū)域林業(yè)綠色全要素生產(chǎn)率提高具有溢出效應(yīng)，區(qū)域內(nèi)林業(yè)綠色全要素生產(chǎn)率差異減小，從而呈現(xiàn)出空間范圍內(nèi)的相似發(fā)展，演變成時(shí)序上的收斂趨勢(shì)?；诖?，提出第一個(gè)研究假說(shuō)：

H1：林業(yè)綠色全要素生產(chǎn)率在區(qū)域內(nèi)由于空間溢出效應(yīng)呈現(xiàn)收斂態(tài)勢(shì)。

地理學(xué)第二定律的提出同時(shí)證實(shí)，事物的空間隔離性造成事物在地理范圍內(nèi)的差異。由于各區(qū)域在資源稟賦、經(jīng)濟(jì)發(fā)展水平、技術(shù)創(chuàng)新能力、生態(tài)環(huán)境及政策制度等方面存在明顯差異，勢(shì)必促使林業(yè)產(chǎn)業(yè)向不同方向發(fā)展，導(dǎo)致區(qū)域間產(chǎn)業(yè)經(jīng)濟(jì)運(yùn)作模式、生產(chǎn)效率的差異。在空間范圍內(nèi)，不同林業(yè)全要素生產(chǎn)率水平的企業(yè)最終出現(xiàn)差異性集聚分布，區(qū)域間林業(yè)綠色全要素生產(chǎn)率呈現(xiàn)時(shí)序發(fā)散趨勢(shì)。基于此，提出第二個(gè)研究假說(shuō)：

H2：林業(yè)綠色全要素生產(chǎn)率在區(qū)域間由于空間異質(zhì)性呈現(xiàn)發(fā)散態(tài)勢(shì)。

三、數(shù)據(jù)來(lái)源、變量選取與模型選擇

（一）數(shù)據(jù)來(lái)源

由于《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》于2004年開始對(duì)我國(guó)環(huán)境和林業(yè)土地等指標(biāo)進(jìn)行系統(tǒng)的核算，因此為保證數(shù)據(jù)的測(cè)算口徑一致，本文選取的研究期間為2004—2019 年。由于我國(guó)香港等地區(qū)數(shù)據(jù)缺失，因此只考察其余30個(gè)省（市、自治區(qū)），數(shù)據(jù)主要來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)林業(yè)統(tǒng)計(jì)年鑒》《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》《中國(guó)能源統(tǒng)計(jì)年鑒》等。其中，固定資產(chǎn)投資額采用固定資產(chǎn)投資指數(shù)平減，其他指數(shù)采用居民消費(fèi)價(jià)格指數(shù)平減。

（二）變量選取

本文研究的核心變量為林業(yè)綠色全要素生產(chǎn)率，參考劉濤等［21］已有研究，確定林業(yè)綠色全要素生產(chǎn)率測(cè)算指標(biāo)體系（表1）。投入要素包括土地、資本、勞動(dòng)力和能源，土地、資本和勞動(dòng)力分別采用林業(yè)用地面積、林業(yè)固定資產(chǎn)投資存量［22］和年末林業(yè)從業(yè)人數(shù)表示，能源投入采用“地區(qū)能源消費(fèi)總量×地區(qū)林業(yè)總產(chǎn)值/地區(qū)生產(chǎn)總值”方法測(cè)算。期望產(chǎn)出包括經(jīng)濟(jì)產(chǎn)出和生態(tài)產(chǎn)出，分別用地區(qū)林業(yè)生產(chǎn)總值和地區(qū)造林面積表示。非期望產(chǎn)出主要為環(huán)境污染物排放量（產(chǎn)出量），采用“地區(qū)工業(yè)污染物排放量（產(chǎn)生量）×地區(qū)林業(yè)第二產(chǎn)業(yè)產(chǎn)值/工業(yè)生產(chǎn)總值”方法測(cè)算林產(chǎn)工業(yè)SO2排放量、林產(chǎn)工業(yè)廢水中COD 排放量和林產(chǎn)工業(yè)固體廢棄物產(chǎn)生量表示。林業(yè)全要素生產(chǎn)率的測(cè)算不考慮土地、能源投入以及非期望產(chǎn)出。

表1 林業(yè)綠色全要素生產(chǎn)率測(cè)算指標(biāo)體系

（二）模型選擇

1.超效率SBM 模型 超效率SBM 模型綜合考慮影響決策單位效率水平的各投入產(chǎn)出要素之間的非線性關(guān)系，確定決策單位的效率水平。相較于傳統(tǒng)數(shù)據(jù)包絡(luò)模型，超效率SBM 模型將非期望產(chǎn)出要素考慮到?jīng)Q策單位的效率測(cè)算，同時(shí)可以非徑向調(diào)整投入產(chǎn)出的松弛變量，最大程度改善效率值。綜合比較，超效率SBM 模型獨(dú)有的特性更加符合林業(yè)綠色全要素生產(chǎn)率測(cè)算需求，因此本文采用超效率SBM 模型測(cè)算林業(yè)綠色全要素生產(chǎn)率，將決策單位不同效率水平加以明確區(qū)分。超效率SBM 模型具體表示為：-

2.泰爾指數(shù) 本文采用泰爾指數(shù)度量各地區(qū)林業(yè)綠色全要素生產(chǎn)率在區(qū)域間的差異，具體表達(dá)式為：

式（2）中，T表示區(qū)域間差異的總體泰爾指數(shù)，n表示研究區(qū)域的數(shù)量，δi表示i地區(qū)的林業(yè)綠色全要素生產(chǎn)率，-δ為其平均值。當(dāng)將n個(gè)地區(qū)劃分為K組，則每組的地區(qū)數(shù)量為nk，總差異可劃分為組間差異（Tb）和組內(nèi)差異（Tm）。

3.核密度函數(shù) 核密度函數(shù)是基于核函數(shù)通過(guò)采用平滑估計(jì)以確定隨機(jī)變量概率密度的分布形態(tài)的方法［23］，采用核密度函數(shù)研究林業(yè)綠色全要素生產(chǎn)率區(qū)域差異演變狀況。具體表達(dá)式為：

式（3）中，K(*)為核密度函數(shù)，h表示帶寬，最大帶寬的確定應(yīng)滿足積分局方誤差值最小。

4.莫蘭指數(shù) 為探究林業(yè)全要素生產(chǎn)率的空間相關(guān)性特征引入莫蘭指數(shù)，全局莫蘭指數(shù)表達(dá)式為：

5.馬爾可夫鏈 馬爾可夫轉(zhuǎn)移是一個(gè)隨機(jī)過(guò)程，其假定研究對(duì)象狀態(tài)只與前一期狀態(tài)有關(guān)。而林業(yè)發(fā)展效率狀態(tài)會(huì)在前一期效率水平上進(jìn)行投入產(chǎn)出調(diào)整以達(dá)到期望效率值，而更前期的效率水平對(duì)其影響較不明顯，因此可以假定林業(yè)綠色全要素生產(chǎn)率符合馬爾可夫轉(zhuǎn)移的應(yīng)用前提。馬爾可夫轉(zhuǎn)移根據(jù)數(shù)據(jù)分布將指標(biāo)值分為l個(gè)狀態(tài)類別，再根據(jù)不同狀態(tài)轉(zhuǎn)移頻數(shù)確定轉(zhuǎn)移概率。具體計(jì)算公式為：

式（6）中，pαβ表示t時(shí)期狀態(tài)為α，后一時(shí)期為β的概率，nαβ表示t時(shí)期林業(yè)綠色全要素生產(chǎn)率狀態(tài)為α，后一時(shí)期為β狀態(tài)的頻數(shù)；nα表示t時(shí)期狀態(tài)為α的個(gè)體數(shù)量?？臻g馬爾可夫鏈可考慮相鄰區(qū)域狀態(tài)影響，采用Rook 臨界方法確定空間權(quán)重，計(jì)算空間滯后類型，根據(jù)同樣計(jì)算概率的方法最終確定l個(gè)（l×l）的條件轉(zhuǎn)移概率矩陣?；隈R爾可夫轉(zhuǎn)移矩陣，對(duì)林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移進(jìn)行長(zhǎng)期預(yù)測(cè)，具體公式為：

式（7）中，S表示長(zhǎng)期馬爾可夫轉(zhuǎn)移概率矩陣，M是已知轉(zhuǎn)移概率矩陣。

四、結(jié)果與分析

（一）林業(yè)綠色全要素生產(chǎn)率與林業(yè)全要素生產(chǎn)率描述性分析

采用DEA Solver 軟件中超效率SBM 模型測(cè)算我國(guó)30 個(gè)?。ㄊ?、自治區(qū)）林業(yè)效率值，由于篇幅有限，此處僅列示2004、2010、2015 以及2019 年代表性年份測(cè)算值（表2）。

表2 代表性年份林業(yè)綠色全要素水平

從表2 看出，2004 年我國(guó)大部分地區(qū)林業(yè)綠色全要素生產(chǎn)率處于低水平和中低水平，而處于中高水平和高水平的省份多集中在黃河流域，包括寧夏、陜西和上海。2010 年相對(duì)于2004 年大部分地區(qū)效率值出現(xiàn)輕微波動(dòng)，但仍多處于低水平。其中甘肅、寧夏和陜西效率值明顯下降，青海和重慶效率值明顯升高，效率值在空間上呈現(xiàn)無(wú)規(guī)律分布。2015 年與2010 年相比，中東部地區(qū)以及云南效率值明顯提高，效率值在空間范圍內(nèi)出現(xiàn)較明顯區(qū)域分布差異。2019 年大部分地區(qū)效率值改善，其中青海、重慶和浙江效率值明顯提高，而效率值較低地區(qū)主要分布在東北地區(qū)和西南地區(qū)以及部分西部地區(qū)。

對(duì)比全國(guó)平均水平（圖1）可得，兩者整體呈現(xiàn)增長(zhǎng)趨勢(shì)，表明我國(guó)林業(yè)生產(chǎn)效率逐漸提高。而林業(yè)綠色全要素生產(chǎn)率明顯高于林業(yè)全要素生產(chǎn)率，且波動(dòng)明顯。表明林業(yè)資源能源投入所帶來(lái)的經(jīng)濟(jì)效益高于發(fā)展所導(dǎo)致的環(huán)境低效率，進(jìn)一步說(shuō)明林業(yè)綠色全要素生產(chǎn)率更能真實(shí)反映林業(yè)發(fā)展現(xiàn)狀。

圖1 全國(guó)林業(yè)全要素生產(chǎn)率及林業(yè)全要素生產(chǎn)率平均水平

對(duì)比四年林業(yè)綠色全要素生產(chǎn)率空間變動(dòng)發(fā)現(xiàn)，優(yōu)勢(shì)區(qū)域從黃河流域轉(zhuǎn)移到中東部地區(qū)，經(jīng)濟(jì)崛起帶動(dòng)林業(yè)發(fā)展，空間集聚效應(yīng)逐漸呈現(xiàn)。林業(yè)產(chǎn)業(yè)經(jīng)過(guò)2010—2019 年十年發(fā)展努力，多數(shù)地區(qū)效率值扭轉(zhuǎn)起初不利局勢(shì)。但東北及西南部分地區(qū)效率值仍較低，表明該部分地區(qū)生產(chǎn)效率改善及整體均衡發(fā)展具有較大提升空間。對(duì)比2004 與2019 年林業(yè)綠色全要素生產(chǎn)率，大部分地區(qū)效率值向更高水平轉(zhuǎn)移，表明林業(yè)正在向高質(zhì)量發(fā)展階段推進(jìn)。西部部分地區(qū)及海南有所降低，主要受制于地區(qū)經(jīng)濟(jì)技術(shù)和自然條件水平。東北和西南地區(qū)效率值雖然始終處于低水平，但整體效率值升高，表明該部分地區(qū)具有突破低效率水平的內(nèi)在發(fā)展?jié)摿Α?/p>

（二）林業(yè)綠色全要素生產(chǎn)率區(qū)域差異分析

1.區(qū)域差異收斂性分析 采用Stata15 軟件計(jì)算2004—2019 年林業(yè)綠色全要素生產(chǎn)率泰爾指數(shù)，并依據(jù)我國(guó)四大經(jīng)濟(jì)區(qū)域進(jìn)行分組，測(cè)算組間和組內(nèi)泰爾指數(shù)（圖2）。從總體泰爾指數(shù)看出，2004—2009 年區(qū)域間差異呈現(xiàn)“鋸齒狀”波動(dòng)，2009 年后波動(dòng)幅度明顯減小。2015—2019 年指數(shù)值達(dá)到平穩(wěn)狀態(tài)，但仍存在發(fā)散趨勢(shì)。而2004—2015 年的組內(nèi)泰爾指數(shù)與整體泰爾指數(shù)的變動(dòng)趨勢(shì)保持一致，表明此期間差異主要來(lái)自組內(nèi)差異，經(jīng)濟(jì)區(qū)域間差異較小。2015—2019 年組間泰爾指數(shù)明顯升高，達(dá)到與組內(nèi)泰爾指數(shù)相當(dāng)?shù)乃?，表明要?shí)現(xiàn)長(zhǎng)期內(nèi)收斂需控制經(jīng)濟(jì)區(qū)域間林業(yè)綠色全要素生產(chǎn)率差異。

圖2 林業(yè)綠色全要素生產(chǎn)率區(qū)域差異

從圖3 各經(jīng)濟(jì)區(qū)域?qū)傮w泰爾指數(shù)貢獻(xiàn)率來(lái)看，中部和東北地區(qū)對(duì)區(qū)域總差異的貢獻(xiàn)程度較穩(wěn)定，東北地區(qū)貢獻(xiàn)最小。在2004—2011 年，東部地區(qū)和西部地區(qū)指數(shù)貢獻(xiàn)率處于互補(bǔ)波動(dòng)狀態(tài)，是總體差異主要來(lái)源。2011 年以后西部地區(qū)的差異不斷減小到與中部地區(qū)貢獻(xiàn)率相近的水平，而東部地區(qū)的貢獻(xiàn)率仍處于上升趨勢(shì)，是未來(lái)區(qū)域差異的主要來(lái)源。

圖3 各經(jīng)濟(jì)地區(qū)泰爾指數(shù)貢獻(xiàn)率

2.區(qū)域差異核密度估計(jì)分析 采用高斯核密度函數(shù)明確地區(qū)差異的演變過(guò)程，同時(shí)采用Stata15軟件繪制2004、2010、2015和2019年的林業(yè)綠色全要素生產(chǎn)率核密度圖（圖4）。

圖4 代表性年份林業(yè)綠色全要素生產(chǎn)率核密度圖

2004 年核密度曲線呈現(xiàn)單峰形式，且峰值處于低于0.5 區(qū)間，表明多數(shù)地區(qū)效率值較低。但核密度曲線存在向右拖尾現(xiàn)象，表明少部分地區(qū)效率值較高。2004—2010 年核密度曲線峰值變動(dòng)不大，但拖尾現(xiàn)象減弱。2015 年核密度曲線峰值明顯下降且向右移動(dòng)，寬度增加，表明2015年效率值分布較2010年更加均勻，整體效率值有所提高。2019 年拖尾現(xiàn)象明顯轉(zhuǎn)化為核密度曲線呈現(xiàn)雙峰形態(tài)，各地區(qū)林業(yè)綠色全要素生產(chǎn)率呈現(xiàn)輕微兩極化分布。雙峰峰值較小，寬度較大，效率值集中度降低。

（三）林業(yè)綠色全要素生產(chǎn)率空間自相關(guān)分析

表3 計(jì)算林業(yè)綠色全要素生產(chǎn)率的全局莫蘭指數(shù)。數(shù)據(jù)顯示，在2004—2012 年，只有2005 年存在輕微的空間正自相關(guān)性，其余年份林業(yè)綠色全要素生產(chǎn)率大多呈現(xiàn)隨機(jī)分散分布狀態(tài)。而在2014—2019年，空間正自相關(guān)明顯增強(qiáng)，林業(yè)發(fā)展的空間集聚特征不斷呈現(xiàn)。

表3 林業(yè)綠色全要素生產(chǎn)率全局莫蘭指數(shù)

為進(jìn)一步測(cè)算2016—2019年林業(yè)綠色全要素生產(chǎn)率局部莫蘭指數(shù)，本文繪制局部莫蘭指數(shù)散點(diǎn)圖進(jìn)行區(qū)域空間集聚特征分析。2016—2019 年，多數(shù)地區(qū)處于第一象限（H-H 集聚區(qū)域）和第三象限（LL集聚區(qū)域），空間正自相關(guān)性明顯，與前文全局莫蘭指數(shù)分析結(jié)果一致，少數(shù)地區(qū)如廣西、浙江、海南處于第二象限（L-H 集聚區(qū)域）。處于第一象限的地區(qū)多為中東部地區(qū)，而處于第三象限的地區(qū)多為西部地區(qū)和東北部地區(qū)，表明效率值較高的地區(qū)集聚在中東部地區(qū)，而效率值較低的地區(qū)集聚在西部以及東北地區(qū)，林業(yè)綠色全要素生產(chǎn)率在空間內(nèi)呈現(xiàn)兩極化分布。

（四）林業(yè)綠色全要素生產(chǎn)率狀態(tài)轉(zhuǎn)移分析

1.短期區(qū)域狀態(tài)轉(zhuǎn)移分析 根據(jù)林業(yè)綠色全要素生產(chǎn)率的分布情況將效率值分為低（Ⅰ）、中低（Ⅱ）、中高（Ⅲ）和高（Ⅳ）4個(gè)水平等級(jí)，分別將2004—2019年30個(gè)?。ㄊ?、自治區(qū)）的效率值依據(jù)不同等級(jí)的范圍進(jìn)行狀態(tài)歸類，得出處于各等級(jí)的個(gè)體頻數(shù)，計(jì)算狀態(tài)轉(zhuǎn)移概率，最終確定傳統(tǒng)馬爾可夫轉(zhuǎn)移概率矩陣和空間馬爾可夫轉(zhuǎn)移概率矩陣（表4）。

表4 林業(yè)綠色全要素生產(chǎn)率傳統(tǒng)和空間馬爾可夫轉(zhuǎn)移概率矩陣 %

對(duì)傳統(tǒng)馬爾可夫轉(zhuǎn)移概率矩陣進(jìn)行分析發(fā)現(xiàn)，馬爾可夫轉(zhuǎn)移概率矩陣主對(duì)角線上概率值均高于非主對(duì)角線上的概率值，類型Ⅰ、Ⅱ、Ⅲ、Ⅳ保持其初始狀態(tài)的概率分別為86.19%、64.46%、48.53%和58.82%，均大于轉(zhuǎn)移為其他狀態(tài)類型的概率。其中類型Ⅰ保持初始狀態(tài)的概率最高，表現(xiàn)低效率水平地區(qū)林業(yè)綠色全要素生產(chǎn)率的持續(xù)性。這表明林業(yè)綠色全要素生產(chǎn)率保留其初始狀態(tài)的概率較高，林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移存在“俱樂部效應(yīng)”。另外，轉(zhuǎn)移概率矩陣中效率值向上轉(zhuǎn)移的概率高于向下轉(zhuǎn)移的概率，如矩陣中類型Ⅱ轉(zhuǎn)移為類型Ⅲ和Ⅳ的概率比例為23.15%，明顯高于轉(zhuǎn)移為類型Ⅰ的概率（12.4%），其他狀態(tài)類型類似。表明林業(yè)綠色全要素生產(chǎn)率具有向好發(fā)展趨勢(shì)，我國(guó)林業(yè)傾向于向更高效率水平發(fā)展。再次，矩陣中鄰近主對(duì)角線的概率值高于非鄰近的概率值，如矩陣中類型Ⅰ轉(zhuǎn)移為類型Ⅱ的概率為10.48%，大于轉(zhuǎn)移為類型Ⅲ和類型Ⅳ的概率（3.33%），此規(guī)律同樣適用于向下轉(zhuǎn)移的情況。這表明林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移對(duì)初始狀態(tài)類型具有一定的依賴性，難以實(shí)現(xiàn)跨級(jí)別轉(zhuǎn)移。

將相鄰區(qū)域效率狀態(tài)考慮到林業(yè)綠色全要素生產(chǎn)率狀態(tài)轉(zhuǎn)移概率測(cè)算中，構(gòu)建空間馬爾可夫轉(zhuǎn)移矩陣發(fā)現(xiàn)，在相應(yīng)類型空間滯后狀態(tài)影響下，林業(yè)綠色全要素生產(chǎn)率保持其初始狀態(tài)的概率升高，如空間滯后類型為Ⅰ時(shí)，類型Ⅰ保持初始狀態(tài)的概率提高到87.13%。這表明在相應(yīng)空間滯后狀態(tài)類型影響下，“俱樂部效應(yīng)”增強(qiáng)，相似的生產(chǎn)效率水平將導(dǎo)致區(qū)域間信息、技術(shù)趨同，林業(yè)綠色全要素生產(chǎn)率狀態(tài)持續(xù)性增強(qiáng)。另外以空間滯后等級(jí)高的區(qū)域?yàn)猷?，效率值向上轉(zhuǎn)移的概率升高，以空間滯后等級(jí)低的區(qū)域?yàn)猷?，效率值向下轉(zhuǎn)移的概率升高。如在不考慮空間滯后時(shí)，類型Ⅰ向上轉(zhuǎn)移為類型Ⅱ、Ⅲ、Ⅳ的總概率為13.81%，在空間滯后類型為Ⅱ的情況下，相應(yīng)轉(zhuǎn)移概率升高到16.09%。表明林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移存在空間溢出效應(yīng)，鄰近區(qū)域效率狀態(tài)對(duì)區(qū)域狀態(tài)演變同時(shí)具有正向和負(fù)向溢出影響。在空間滯后影響下，林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移對(duì)初始狀態(tài)仍具有較強(qiáng)依賴性。如空間滯后類型為Ⅰ的情況下，類型Ⅰ向類型Ⅱ轉(zhuǎn)移的概率（10.89%）高于向類型Ⅲ和Ⅳ轉(zhuǎn)移的總概率（1.98%）。這表明林業(yè)發(fā)展的流動(dòng)性較弱，低效率水平地區(qū)短期內(nèi)仍難以向高端狀態(tài)推進(jìn)。

2.長(zhǎng)期區(qū)域狀態(tài)空間演變預(yù)測(cè) 為探究林業(yè)綠色全要素生產(chǎn)率狀態(tài)的長(zhǎng)期變化，采用R 軟件對(duì)馬爾可夫轉(zhuǎn)移概率矩陣進(jìn)行長(zhǎng)期預(yù)測(cè)，相對(duì)于短期狀態(tài)轉(zhuǎn)移規(guī)律存在一定的異同之處。

長(zhǎng)期內(nèi)林業(yè)綠色全要素生產(chǎn)率狀態(tài)轉(zhuǎn)移仍難以實(shí)現(xiàn)趨同，呈現(xiàn)向多個(gè)狀態(tài)類型轉(zhuǎn)移的趨勢(shì)。從表5 可以看出，在不考慮空間滯后情況下，四種類型效率值轉(zhuǎn)移呈現(xiàn)較均勻分布，在空間滯后類型為Ⅰ的情況下，長(zhǎng)期轉(zhuǎn)移概率集中在類型Ⅰ和類型Ⅱ；林業(yè)綠色全要素生產(chǎn)率長(zhǎng)期轉(zhuǎn)移具有空間收斂性，即效率值趨于向空間滯后類型方向轉(zhuǎn)移。表現(xiàn)在矩陣中，狀態(tài)轉(zhuǎn)移為空間滯后類型的概率最高。如在空間滯后類型為Ⅰ的情況下，長(zhǎng)期轉(zhuǎn)移為Ⅰ類型的概率最高（46.11%），表明長(zhǎng)期內(nèi)相鄰區(qū)域生產(chǎn)技術(shù)的溢出影響將使地區(qū)林業(yè)生產(chǎn)效率向空間滯后類型收斂，促使其呈現(xiàn)空間集聚分布。另外類型Ⅲ長(zhǎng)期內(nèi)趨于向類型Ⅰ轉(zhuǎn)移，表明林業(yè)發(fā)展處于中高水平的地區(qū)在空間上過(guò)度集聚可能導(dǎo)致惡性競(jìng)爭(zhēng)、市場(chǎng)失靈，進(jìn)而影響產(chǎn)業(yè)生產(chǎn)效率。這警示處于中高水平的林業(yè)產(chǎn)業(yè)地區(qū)，若不及時(shí)進(jìn)行產(chǎn)業(yè)調(diào)整，長(zhǎng)期內(nèi)仍可能陷入產(chǎn)業(yè)低效率發(fā)展陷阱；長(zhǎng)期內(nèi)，林業(yè)綠色全要素生產(chǎn)率轉(zhuǎn)移具有空間溢出性，表現(xiàn)為當(dāng)空間滯后類型為低水平或中低水平時(shí)，狀態(tài)轉(zhuǎn)移概率將集中在低水平或中低水平方向。而空間滯后類型為高水平時(shí)，狀態(tài)轉(zhuǎn)移概率將集中在高水平方向。該現(xiàn)象表明空間溢出性使得林業(yè)綠色全要素生產(chǎn)率在空間上呈現(xiàn)集聚分布。如空間滯后類型為Ⅰ的情況下，長(zhǎng)期狀態(tài)轉(zhuǎn)移概率集中在向類型Ⅰ和類型Ⅱ轉(zhuǎn)移，概率分別為46.11%和28.67%。相似結(jié)論適用于空間滯后類型為Ⅱ和Ⅳ的情況。

表5 林業(yè)綠色全要素生產(chǎn)率長(zhǎng)期馬爾可夫轉(zhuǎn)移概率預(yù)測(cè)矩陣 %

五、主要結(jié)論與政策建議

采用超效率SBM 模型測(cè)算2004—2019 年我國(guó)30 個(gè)省（市、自治區(qū)）的林業(yè)綠色全要素生產(chǎn)率，在此基礎(chǔ)上進(jìn)行時(shí)空演變及狀態(tài)轉(zhuǎn)移分析。得到以下研究結(jié)論：第一，我國(guó)多數(shù)地區(qū)林業(yè)綠色全要素生產(chǎn)率實(shí)現(xiàn)了向更高效率水平轉(zhuǎn)移，優(yōu)勢(shì)區(qū)域隨經(jīng)濟(jì)重心轉(zhuǎn)向中東部地區(qū)。第二，林業(yè)綠色全要素生產(chǎn)率地區(qū)間差異波動(dòng)下降，且差異主要來(lái)自經(jīng)濟(jì)區(qū)域內(nèi)部，但經(jīng)濟(jì)區(qū)域間差異影響不斷提高。第三，2013 年后林業(yè)綠色全要素生產(chǎn)率空間正相關(guān)性明顯，高效率區(qū)域聚集在中東部地區(qū)，低效率區(qū)域集聚在西部和東北地區(qū)，呈現(xiàn)出明顯的兩極分化局勢(shì)。第四，從狀態(tài)轉(zhuǎn)移分析看出，短期內(nèi)林業(yè)綠色全要素生產(chǎn)率具有向好發(fā)展趨勢(shì)，但長(zhǎng)期仍難以趨同；林業(yè)綠色全要素生產(chǎn)率狀態(tài)轉(zhuǎn)移存在“俱樂部效應(yīng)”，對(duì)初始狀態(tài)的依賴使得狀態(tài)轉(zhuǎn)移流動(dòng)性較弱，而由于長(zhǎng)期收斂，空間集聚分布明顯，但過(guò)度空間集聚將導(dǎo)致中高林業(yè)綠色全要素生產(chǎn)率水平區(qū)域向低水平轉(zhuǎn)移。此外，林業(yè)綠色全要素生產(chǎn)率狀態(tài)轉(zhuǎn)移存在空間溢出性。

基于上述研究結(jié)論，提出如下政策建議：首先，協(xié)同區(qū)域經(jīng)濟(jì)發(fā)展，共建良好宏觀經(jīng)濟(jì)環(huán)境。推進(jìn)東北地區(qū)振興和西部大開發(fā)戰(zhàn)略，加速地區(qū)優(yōu)勢(shì)資源引進(jìn)，完善地區(qū)體制建設(shè)。同時(shí)要加強(qiáng)林業(yè)政策支持和善林業(yè)發(fā)展基礎(chǔ)建設(shè)，完善林業(yè)產(chǎn)權(quán)保護(hù)相關(guān)法律法規(guī)，提高林場(chǎng)主生產(chǎn)積極性。其次，推動(dòng)權(quán)利適度下放，引導(dǎo)區(qū)域創(chuàng)新競(jìng)爭(zhēng)，有效利用產(chǎn)業(yè)集聚優(yōu)勢(shì)。政府要嚴(yán)格監(jiān)督林業(yè)產(chǎn)業(yè)市場(chǎng)價(jià)格，明確賞罰機(jī)制。最后，促進(jìn)區(qū)域交流，強(qiáng)化林業(yè)產(chǎn)業(yè)正向溢出作用，深化“以強(qiáng)扶弱”發(fā)展理念。推動(dòng)落實(shí)中東部地區(qū)人才資源、信息技術(shù)向西部及東北部地區(qū)引進(jìn)的優(yōu)惠政策，完善區(qū)域間產(chǎn)業(yè)鏈，加強(qiáng)林業(yè)產(chǎn)業(yè)市場(chǎng)建設(shè)，協(xié)調(diào)推進(jìn)地區(qū)林業(yè)產(chǎn)業(yè)發(fā)展。