呂潔華，孫嘉宇，蔡秀亭

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

一、引言與文獻綜述

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

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

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

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

綜合來看，國內(nèi)外學者對林業(yè)全要素生產(chǎn)率的測算及影響因素研究比較充足，但對于林業(yè)綠色全要素生產(chǎn)率的研究仍處于起步階段，缺少對其時空演變的定量分析以及狀態(tài)轉(zhuǎn)移分析?；诖耍疚牟捎贸蔛BM 模型重新測算我國林業(yè)綠色全要素生產(chǎn)率，結合區(qū)域差異性分析及馬爾可夫轉(zhuǎn)移客觀探究其演變歷程并進行長期預測，把握基本演變規(guī)律，為推動林業(yè)“經(jīng)濟生態(tài)化，生態(tài)經(jīng)濟化”發(fā)展提供科學依據(jù)。

二、理論分析與研究假說

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

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

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

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

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

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

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

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

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

（二）變量選取

本文研究的核心變量為林業(yè)綠色全要素生產(chǎn)率，參考劉濤等［21］已有研究，確定林業(yè)綠色全要素生產(chǎn)率測算指標體系（表1）。投入要素包括土地、資本、勞動力和能源，土地、資本和勞動力分別采用林業(yè)用地面積、林業(yè)固定資產(chǎn)投資存量［22］和年末林業(yè)從業(yè)人數(shù)表示，能源投入采用“地區(qū)能源消費總量×地區(qū)林業(yè)總產(chǎn)值/地區(qū)生產(chǎn)總值”方法測算。期望產(chǎn)出包括經(jīng)濟產(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)總值”方法測算林產(chǎn)工業(yè)SO2排放量、林產(chǎn)工業(yè)廢水中COD 排放量和林產(chǎn)工業(yè)固體廢棄物產(chǎn)生量表示。林業(yè)全要素生產(chǎn)率的測算不考慮土地、能源投入以及非期望產(chǎn)出。

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

（二）模型選擇

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

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

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

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

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

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

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

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

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

四、結果與分析

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

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

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

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

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

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

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

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

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

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

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

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

2.區(qū)域差異核密度估計分析 采用高斯核密度函數(shù)明確地區(qū)差異的演變過程，同時采用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 年核密度曲線峰值變動不大，但拖尾現(xiàn)象減弱。2015 年核密度曲線峰值明顯下降且向右移動，寬度增加，表明2015年效率值分布較2010年更加均勻，整體效率值有所提高。2019 年拖尾現(xiàn)象明顯轉(zhuǎn)化為核密度曲線呈現(xiàn)雙峰形態(tài)，各地區(qū)林業(yè)綠色全要素生產(chǎn)率呈現(xiàn)輕微兩極化分布。雙峰峰值較小，寬度較大，效率值集中度降低。

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

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

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

為進一步測算2016—2019年林業(yè)綠色全要素生產(chǎn)率局部莫蘭指數(shù)，本文繪制局部莫蘭指數(shù)散點圖進行區(qū)域空間集聚特征分析。2016—2019 年，多數(shù)地區(qū)處于第一象限（H-H 集聚區(qū)域）和第三象限（LL集聚區(qū)域），空間正自相關性明顯，與前文全局莫蘭指數(shù)分析結果一致，少數(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個水平等級，分別將2004—2019年30個?。ㄊ?、自治區(qū)）的效率值依據(jù)不同等級的范圍進行狀態(tài)歸類，得出處于各等級的個體頻數(shù)，計算狀態(tài)轉(zhuǎn)移概率，最終確定傳統(tǒng)馬爾可夫轉(zhuǎn)移概率矩陣和空間馬爾可夫轉(zhuǎn)移概率矩陣（表4）。

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

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

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

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

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

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

五、主要結論與政策建議

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

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