吳 寧，戴永務(wù)，林偉明，王強強

(1. 福建農(nóng)林大學(xué) 經(jīng)濟管理學(xué)院，福州 350002；2.福建農(nóng)林大學(xué) 金山學(xué)院，福州 3500023.福建高校特色新型智庫集體林業(yè)改革發(fā)展研究中心，福州 350002)

林業(yè)產(chǎn)業(yè)是中國國民經(jīng)濟基礎(chǔ)產(chǎn)業(yè)，發(fā)展林業(yè)是把綠水青山變成金山銀山的重要途徑。截止2020年，中國林業(yè)產(chǎn)業(yè)總產(chǎn)值達7.55萬億元，林產(chǎn)品對外貿(mào)易總額達 1 600 億美元，創(chuàng)造 3 400 多萬個就業(yè)崗位[1]。近年來中國林業(yè)產(chǎn)業(yè)發(fā)展取得了矚目的成績，但各地區(qū)林業(yè)產(chǎn)業(yè)發(fā)展?fàn)顩r差異較大。已有學(xué)者測度省域林業(yè)產(chǎn)業(yè)集聚水平[2-4]，也有學(xué)者按不同區(qū)域分別測度了西南省份、東北國有林區(qū)、長江經(jīng)濟帶各省林業(yè)產(chǎn)業(yè)集聚水平[5-7]，結(jié)果表明，林業(yè)產(chǎn)業(yè)區(qū)域發(fā)展異質(zhì)性顯著，形成了五大各具特色林業(yè)產(chǎn)業(yè)集聚[8]。探究林業(yè)產(chǎn)業(yè)區(qū)域發(fā)展異質(zhì)性的根源，新結(jié)構(gòu)經(jīng)濟學(xué)提出一個國家或經(jīng)濟地區(qū)在一定的時間段的要素稟賦是給定的，認為經(jīng)濟結(jié)構(gòu)內(nèi)生決定于要素稟賦結(jié)構(gòu)[9]。當(dāng)?shù)氐囊赝度虢Y(jié)構(gòu)與產(chǎn)業(yè)空間分布集聚態(tài)勢密切相關(guān)[10]，土地、勞動、資本三種要素在空間上的分布離散或者集中直接影響了產(chǎn)業(yè)的區(qū)位選擇[11]，可見要素稟賦結(jié)構(gòu)對中國林業(yè)產(chǎn)業(yè)發(fā)展格局具有重大影響。隨著時間的推移要素稟賦結(jié)構(gòu)也會發(fā)生變化，在林業(yè)產(chǎn)業(yè)發(fā)展過程中，隨著要素投入的種類、數(shù)量、質(zhì)量等發(fā)生變動，導(dǎo)致要素稟賦結(jié)構(gòu)發(fā)生改變[12-14]。另外，林業(yè)要素稟賦結(jié)構(gòu)在林業(yè)技術(shù)水平進步和林業(yè)要素相對價格變動等因素的共同影響作用下也會產(chǎn)生流動[15]?，F(xiàn)有研究中，學(xué)者們對林業(yè)要素稟賦結(jié)構(gòu)時空演變的關(guān)注不夠，缺乏從空間的角度對中國林業(yè)要素稟賦結(jié)構(gòu)變化的空間特征進行深入分析。那么，中國各地區(qū)林業(yè)要素稟賦結(jié)構(gòu)的不斷變化是否存在空間上的相互依賴性？區(qū)域間存在何種差異？林業(yè)要素稟賦結(jié)構(gòu)發(fā)生變化的根本原因主要有哪些？為了正確回答上述研究問題，本文以中國省級(直轄市、自治區(qū))面板數(shù)據(jù)為實驗樣本，分析中國林業(yè)要素稟賦結(jié)構(gòu)變化內(nèi)在形成原因，并通過構(gòu)建林業(yè)要素稟賦結(jié)構(gòu)指數(shù)，采用具有探索性的空間數(shù)據(jù)分析方法對中國林業(yè)要素稟賦結(jié)構(gòu)的整體時空異質(zhì)性變化進行深入研究分析，期望對中國林業(yè)要素稟賦結(jié)構(gòu)的時空異質(zhì)分布變化形成一個更深入的科學(xué)認識，為政府制定科學(xué)合理的林業(yè)要素投入方案提供信息參考。

1 研究方法與數(shù)據(jù)來源

1.1 林業(yè)要素稟賦結(jié)構(gòu)變化的計算方法

國際上慣用資源稟賦系數(shù)來衡量一個區(qū)域某種資源豐裕度[16]。其中，林業(yè)各要素稟賦系數(shù)的計算公式如下：

Sit=(sit/st)/(yit/yt)

(1)

Lit=(lit/lt)/(yit/yt)

(2)

Kit=(kit/kt)/(yit/yt)

(3)

其中，sit、lit、kit分別指的是i地區(qū)t時期的林業(yè)土地、林業(yè)勞動力、林業(yè)資本；st、lt、kt分別代表全國t時期林業(yè)土地、林業(yè)勞動力、林業(yè)資本；yit代表i地區(qū)t時期的林業(yè)產(chǎn)業(yè)總產(chǎn)值、yt代表t時期全國林業(yè)產(chǎn)業(yè)總產(chǎn)值；Sit、Lit、Kit分別指的是i地區(qū)的林業(yè)土地要素稟賦系數(shù)、林業(yè)勞動力要素稟賦系數(shù)、林業(yè)資本要素稟賦系數(shù)。利用林業(yè)各要素稟賦系數(shù)可以構(gòu)建反映各地區(qū)要素投入結(jié)構(gòu)的要素結(jié)構(gòu)指數(shù)[16]。林業(yè)要素稟賦結(jié)構(gòu)指數(shù)(ESI)如下：

ESISit=Sit/(Sit+Kit+Lit)

(4)

ESILit=Lit/(Sit+Kit+Lit)

(5)

ESIKit=Kit/(Sit+Kit+Lit)

(6)

其中，ESISit、ESILit、ESIKit分別指的是i地區(qū)t時期的林業(yè)土地要素結(jié)構(gòu)指數(shù)、林業(yè)勞動力要素結(jié)構(gòu)指數(shù)、林業(yè)資本要素結(jié)構(gòu)指數(shù)。

1.2 林業(yè)要素稟賦結(jié)構(gòu)變化的空間特性檢驗方法

為了研究中國林業(yè)要素稟賦結(jié)構(gòu)的時空相關(guān)性及異質(zhì)性，本文主要采用具有探索性的空間數(shù)據(jù)分析方法來對其進行計量檢驗。一般情況下，Morans’I分為全局性和局域性兩類，其作用各異。其中，從總體上反映了林業(yè)要素稟賦結(jié)構(gòu)的空間特性的全局性的Morans’I函數(shù)表達式定義為：

(7)

局域性Morans’I及相應(yīng)的莫蘭散點圖則可對林業(yè)要素稟賦結(jié)構(gòu)的局域空間相互關(guān)系及異質(zhì)性特征進行分析，局域性Morans’I的表達式如下：

(8)

1.3 數(shù)據(jù)來源

本文的數(shù)據(jù)主要包含林業(yè)產(chǎn)出及林業(yè)要素兩大類。本文應(yīng)用各地區(qū)林業(yè)產(chǎn)業(yè)總產(chǎn)值表示林業(yè)產(chǎn)出變量。林業(yè)要素包括土地、勞動力和資本三類。森林面積更能反應(yīng)林業(yè)土地面積的變化，因此以森林面積表征林業(yè)土地要素狀況；其中勞動力數(shù)據(jù)指的是林業(yè)系統(tǒng)各地區(qū)從業(yè)人員年末總?cè)藬?shù)；對于林業(yè)資本要素投入，本文確定基期2003年的資本存量為當(dāng)年的固定資產(chǎn)投入額除以10%作為當(dāng)年的資本存量，以后各期的資本存量按照永續(xù)盤存法推算，對應(yīng)每年固定資產(chǎn)的折舊率為9.6%[17]。上述指標(biāo)數(shù)據(jù)來源于2003—2018年的《中國林業(yè)統(tǒng)計年鑒》。

2 結(jié)果與分析

2.1 林業(yè)要素稟賦結(jié)構(gòu)演變特征

為了從宏觀上把握中國林業(yè)要素稟賦結(jié)構(gòu)的總體變化，本文依據(jù)林業(yè)要素稟賦系數(shù)測算出2003—2018年林業(yè)要素稟賦結(jié)構(gòu)指數(shù)(詳見表1)。林業(yè)土地要素結(jié)構(gòu)指數(shù)均值由2003年的0.317增至2018年的0.379，增長0.062；林業(yè)勞動力要素結(jié)構(gòu)指數(shù)均值由2003年0.277增至0.355，增長0.078；林業(yè)資本要素結(jié)構(gòu)指數(shù)均值由2003年的0.406減至2018年的0.266，降低0.140?？梢娏謽I(yè)土地、勞動力要素結(jié)構(gòu)指數(shù)呈現(xiàn)出穩(wěn)健增長趨勢，而林業(yè)資本要素結(jié)果指數(shù)呈現(xiàn)明顯降低趨勢。從區(qū)域林業(yè)要素稟賦結(jié)構(gòu)指數(shù)變化狀況分析看，2003—2018中國華北地區(qū)資本要素稟賦結(jié)構(gòu)最高，勞動力要素稟賦結(jié)構(gòu)指數(shù)與林業(yè)土地結(jié)構(gòu)指數(shù)相當(dāng)。東北地區(qū)林業(yè)勞動力要素結(jié)構(gòu)指數(shù)較高，林業(yè)土地結(jié)構(gòu)指數(shù)與林業(yè)資本結(jié)構(gòu)指數(shù)呈現(xiàn)下降趨勢。西南地區(qū)的林業(yè)土地要素結(jié)構(gòu)指數(shù)較高，林業(yè)土地要素結(jié)構(gòu)指數(shù)在0.492上下浮動，勞動力要素結(jié)構(gòu)指數(shù)呈現(xiàn)上升趨勢，資本結(jié)構(gòu)指數(shù)呈現(xiàn)下降趨勢?？傮w來看，中國林業(yè)土地結(jié)構(gòu)指數(shù)較高的地區(qū)逐漸由東北地區(qū)向西南地區(qū)轉(zhuǎn)移，東北地區(qū)聚集了較多的林業(yè)勞動力，華東地區(qū)、華中地區(qū)、華南地區(qū)、西北地區(qū)林業(yè)各要素稟賦結(jié)構(gòu)指數(shù)相當(dāng)，但各要素稟賦結(jié)構(gòu)指數(shù)隨著時間的推移存在不同程度波動，可以預(yù)測各地區(qū)之間存在要素流動。

表1 2003—2018林業(yè)要素稟賦結(jié)構(gòu)指數(shù)

續(xù)表1

2.2 林業(yè)要素稟賦結(jié)構(gòu)的全局自相關(guān)性

圖1 2003—2018年中國林業(yè)要素稟賦結(jié)構(gòu)指數(shù)全局性莫蘭指數(shù)

本文利用GEODA空間分析工具，依據(jù)式(7)計算出了2003—2018年中國林業(yè)要素結(jié)構(gòu)指數(shù)的全局性Morans’I(圖1)。從圖1可以看出，2003—2018年中國林業(yè)要素稟賦結(jié)構(gòu)存在顯著的空間正相關(guān)性。其中，林業(yè)土地要素結(jié)構(gòu)指數(shù)的全局性Morans’I均值為0.202，林業(yè)勞動力要素結(jié)構(gòu)指數(shù)的全局性Morans’I均值為0.326，林業(yè)資本要素結(jié)構(gòu)指數(shù)的全局性Morans’I均值為0.194?？梢?，林業(yè)勞動力要素稟賦的空間正相關(guān)性最為突出，即勞動力要素的空間集聚性最強。根據(jù)變化趨勢進行分析，林業(yè)土地要素結(jié)構(gòu)指數(shù)和林業(yè)資本要素結(jié)構(gòu)指數(shù)的全局自相關(guān)性變化趨勢高度一致，均呈現(xiàn)下降趨勢，林業(yè)勞動力的全局自相關(guān)性波動較大且總體上呈現(xiàn)下降趨勢。該結(jié)果可以說明，土地要素本身具有較強的空間位置特征，很難去改變其具體地理空間位置，且林業(yè)資本投入其穩(wěn)定性較好，因此林業(yè)土地、資本要素結(jié)構(gòu)全局性Morans’I呈現(xiàn)平滑的發(fā)展趨勢，未出現(xiàn)大波動。

2.3 林業(yè)要素稟賦結(jié)構(gòu)的時空異質(zhì)性

全局性Morans’I只能從總體上說明中國林業(yè)要素稟賦結(jié)構(gòu)存在空間自相關(guān)性，不能反映各地區(qū)之間的林業(yè)要素稟賦結(jié)構(gòu)的異質(zhì)性特征。本文選取具有代表性年份的莫蘭散點圖，來分析局部林業(yè)要素稟賦的空間異質(zhì)性特征。

2.3.1 林業(yè)土地要素結(jié)構(gòu)的時空異質(zhì)性

由于全局自相關(guān)性數(shù)據(jù)分析結(jié)果發(fā)現(xiàn)，中國林業(yè)土地要素結(jié)構(gòu)指數(shù)在2007年前后發(fā)生了較大幅度波動，所以本文選擇了2003年、2007年、2018年的三個數(shù)據(jù)來分析其區(qū)域林業(yè)土地要素空間集聚變動情況。2003—2018年期間部分地區(qū)的林業(yè)土地要素結(jié)構(gòu)指數(shù)特征發(fā)生了較大變化，對比三個時間點發(fā)現(xiàn)，西藏的林業(yè)土地要素結(jié)構(gòu)指數(shù)的空間相關(guān)類型由2003年的低低集聚類型逐漸轉(zhuǎn)變?yōu)?007年的高高集聚類型，與上文西南地區(qū)林業(yè)要素稟賦結(jié)構(gòu)指數(shù)描述的特征一致。2003—2018年期間中國各地區(qū)林業(yè)土地要素結(jié)構(gòu)指數(shù)的變化在時空表現(xiàn)出正相關(guān)的變化特點，但也仍然存在一定的時空異質(zhì)性，林業(yè)土地要素結(jié)構(gòu)的空間差異呈現(xiàn)出擴大趨勢。

圖2 2003年土地結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖3 2007年土地結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖4 2018年土地結(jié)構(gòu)指數(shù)的莫蘭散點圖

2.3.2 林業(yè)勞動力要素結(jié)構(gòu)的時空異質(zhì)性

由于中國林業(yè)勞動力要素稟賦結(jié)構(gòu)指數(shù)在2007年達到最大值，據(jù)此本文分別選取了2003年、2007年及2018年全國林業(yè)勞動力要素結(jié)構(gòu)指數(shù)的數(shù)據(jù)對其空間集聚特性進行綜合分析(詳見圖5、圖6、圖7)。總體上，中國林業(yè)地區(qū)林業(yè)勞動力要素結(jié)構(gòu)指數(shù)也呈現(xiàn)出正相關(guān)空間依賴性。林業(yè)要素稟賦結(jié)構(gòu)指數(shù)在2003較均勻分布在四個象限，至2018年集中分布在一三象限，林業(yè)勞動力要素空間聚類顯著，且林業(yè)勞動力結(jié)構(gòu)指數(shù)相近的地區(qū)在空間上更容易集聚。從三個時間節(jié)點來分析，部分地區(qū)的林業(yè)勞動力要素集聚類型發(fā)生了巨大的變化。滬、鄂、渝、浙、桂五個地區(qū)的集聚類型由低高集聚轉(zhuǎn)向低低集聚。由此充分說明了無論是區(qū)域內(nèi)還是區(qū)域間，林業(yè)勞動力要素流動狀態(tài)活躍，林業(yè)勞動力空間集聚類型發(fā)生變遷，林業(yè)勞動力空間異質(zhì)性顯著。

圖5 2003年勞動力結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖6 2007年勞動力結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖7 2018年勞動力結(jié)構(gòu)指數(shù)的莫蘭散點圖

2.3.3 林業(yè)資本要素結(jié)構(gòu)的時空異質(zhì)性

根據(jù)全局莫蘭指數(shù)可知林業(yè)資本要素結(jié)構(gòu)全局莫蘭指數(shù)在2007年時達到最大均值，據(jù)此本文分別選取了2003年、2007年及2018年林業(yè)資本要素結(jié)構(gòu)指數(shù)對林業(yè)資本空間集聚類型進行分析。結(jié)合統(tǒng)計圖8、圖9、圖10可以發(fā)現(xiàn)，2003—2018年間的中國林業(yè)資本要素結(jié)構(gòu)指數(shù)分布逐漸形成向一二三象限分布的趨勢，全國資本要素呈現(xiàn)出低低集聚類型包圍的趨勢。該結(jié)果說明，2003—2018年期間林業(yè)資本要素結(jié)構(gòu)指數(shù)的分布表現(xiàn)出區(qū)域的時空異質(zhì)性，且林業(yè)資本要素的低低集聚類型在逐漸擴大，與林業(yè)資本要素稟賦結(jié)構(gòu)指數(shù)減少的趨勢相符。

圖8 2003年林業(yè)資本結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖9 2007年林業(yè)資本結(jié)構(gòu)指數(shù)的莫蘭散點圖

圖10 2018年林業(yè)資本結(jié)構(gòu)指數(shù)的莫蘭散點圖

3 研究結(jié)論及政策啟示

本文通過構(gòu)建林業(yè)要素稟賦結(jié)構(gòu)指數(shù)分析了2003—2018年以來中國林業(yè)要素稟賦結(jié)構(gòu)的發(fā)展演變過程，應(yīng)用具有探索性的空間數(shù)據(jù)分析方法對中國林業(yè)要素結(jié)構(gòu)稟賦指數(shù)結(jié)構(gòu)的變化空間自相關(guān)性及其空間異質(zhì)性問題進行了深入研究分析，可以得到如下結(jié)論：①林業(yè)土地、勞動力要素稟賦結(jié)構(gòu)指數(shù)呈現(xiàn)出穩(wěn)健增長趨勢，而林業(yè)資本要素結(jié)構(gòu)指數(shù)下降明顯；東北地區(qū)勞動力結(jié)構(gòu)指數(shù)較大，說明東北地區(qū)從事林業(yè)產(chǎn)業(yè)勞動者較多；中國林業(yè)土地結(jié)構(gòu)指數(shù)較高的地區(qū)逐漸由東北地區(qū)向西南地區(qū)轉(zhuǎn)移；②林業(yè)要素稟賦結(jié)構(gòu)指數(shù)存在顯著的空間正相關(guān)性，林業(yè)勞動力要素稟賦結(jié)構(gòu)指數(shù)的空間正相關(guān)性最為突出。根據(jù)變化趨勢進行分析，林業(yè)土地要素結(jié)構(gòu)指數(shù)和林業(yè)資本要素結(jié)構(gòu)指數(shù)的全局自相關(guān)性變化趨勢高度一致，且均呈現(xiàn)下降趨勢；③林業(yè)要素稟賦結(jié)構(gòu)指數(shù)呈現(xiàn)出局部的空間異質(zhì)性，林業(yè)土地要素結(jié)構(gòu)指數(shù)的空間差異呈現(xiàn)出擴大趨勢，部分地區(qū)空間集聚類型發(fā)生了明顯變化；林業(yè)勞動力要素在時空上呈現(xiàn)出向高高集聚類型或低低集聚類型兩極分化的發(fā)展趨勢；林業(yè)資本要素結(jié)構(gòu)指數(shù)分布逐漸趨向一二三象限分布，且林業(yè)資本要素的低低集聚類型在逐漸擴大。

林業(yè)要素稟賦結(jié)構(gòu)決定了中國林業(yè)產(chǎn)業(yè)貿(mào)易的比較優(yōu)勢，林業(yè)要素稟賦結(jié)構(gòu)發(fā)生變化正是改變了當(dāng)前中國林業(yè)相關(guān)產(chǎn)業(yè)的貿(mào)易比較優(yōu)勢之處所在。中國林產(chǎn)品對外貿(mào)易市場結(jié)構(gòu)也應(yīng)及時做出一些相應(yīng)的優(yōu)化調(diào)整，才能持續(xù)有效地不斷提升當(dāng)前中國林產(chǎn)品的國際市場競爭力。各省和地區(qū)應(yīng)根據(jù)自身要素稟賦結(jié)構(gòu)特征采取具有針對性的政策安排，優(yōu)化林業(yè)要素稟賦市場，促進中國林業(yè)要素稟賦結(jié)構(gòu)的不斷優(yōu)化轉(zhuǎn)型升級。