東北地區(qū)國土開發(fā)強度的時空格局及其驅(qū)動機制

趙 雲(yún) 泰，王 亞 男，呂 曉*

(1.中國國土勘測規(guī)劃院，北京 100035；2.東北大學文法學院，遼寧省自然資源廳土地保護與利用重點實驗室，遼寧 沈陽 110169)

0 引言

國土空間開發(fā)是基于區(qū)域資源數(shù)量及有效性，強調(diào)以人為本，重視區(qū)域社會經(jīng)濟發(fā)展和人居環(huán)境建設，對建設用地進行高強度、大規(guī)模開發(fā)的復雜多維立體的人類活動[1,2]。隨著城鎮(zhèn)化、工業(yè)化快速發(fā)展，中國所面臨的資源枯竭、環(huán)境惡化、區(qū)域發(fā)展失調(diào)及國土空間開發(fā)失衡等問題日益嚴重[3,4]，“不平衡不充分”問題在我國國土開發(fā)過程中尤為突出，影響人與自然協(xié)調(diào)發(fā)展[5]。因此，探究區(qū)域國土開發(fā)強度的時空演變特征及驅(qū)動機制，對優(yōu)化國土空間利用結(jié)構(gòu)、構(gòu)建基于資源環(huán)境承載力的國土空間開發(fā)格局、促進生態(tài)文明建設至關(guān)重要。

目前，國內(nèi)外關(guān)于國土開發(fā)強度的研究十分豐富，不少學者運用耦合協(xié)調(diào)度模型研究區(qū)域內(nèi)國土開發(fā)強度與資源環(huán)境承載力、土地利用效益和資源環(huán)境水平的互動演化規(guī)律[6-9]；趙亞莉等不僅描述區(qū)域國土開發(fā)強度的演變格局，還運用地理探測器、地理加權(quán)回歸模型等計量方法研究經(jīng)濟、產(chǎn)業(yè)、投資等因素對國土開發(fā)強度的影響機制[10-15]；楊浩等在國土空間開發(fā)的基礎(chǔ)上，進一步探討了土地開發(fā)均衡度及國土空間開發(fā)利用效率[2,16-18]?？傮w看，目前多數(shù)研究重點關(guān)注快速城鎮(zhèn)化區(qū)域和經(jīng)濟相對發(fā)達地區(qū)，對東北地區(qū)的區(qū)域性問題關(guān)注不夠，對國土開發(fā)強度演化特征的描述也有待進一步豐富?；诖?，本文基于東北地區(qū)40個城市(地區(qū)、盟)的土地利用數(shù)據(jù)和社會經(jīng)濟數(shù)據(jù)，運用泰爾指數(shù)、莫蘭指數(shù)、熱點分析揭示1990-2018年東北地區(qū)國土開發(fā)強度的時空演變格局，并從社會經(jīng)濟、自然環(huán)境、財政投資等方面建立指標體系，利用地理探測器開展東北地區(qū)國土開發(fā)強度的驅(qū)動機制分析，以期為揭示東北地區(qū)國土空間開發(fā)格局的內(nèi)在機制、優(yōu)化其國土空間開發(fā)格局提供科學參考。

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

1.1 研究區(qū)概況

本文研究區(qū)包括黑龍江、吉林、遼寧三省陸域及內(nèi)蒙古東四盟地區(qū)(115°E～135°E，38°N～55°N)，區(qū)域國土面積126萬km2(圖1)。東北地區(qū)屬于溫帶大陸性季風氣候，四季分明，7-10月降水較多[19]；地形以平原、高原、山地為主，內(nèi)有三江平原、遼河平原、松嫩平原，外有長白山、大小興安嶺、內(nèi)蒙古高原作為天然屏障，地勢整體中間低、兩側(cè)高；境內(nèi)河網(wǎng)密布，但水資源分布集中在東部地區(qū)。東北地區(qū)擁有豐富的礦產(chǎn)、森林資源和肥沃的黑土，不僅是我國工業(yè)發(fā)展的搖籃，更是我國重要的商品糧基地和生態(tài)安全屏障，陸域與俄羅斯、蒙古、朝鮮相鄰，與韓國、日本隔海相望，地緣優(yōu)勢突出，可依托“一帶一路”戰(zhàn)略擴大對外開放優(yōu)勢。

圖1 研究區(qū)行政區(qū)劃及高程示意

改革開放后，依靠計劃經(jīng)濟時期國家重點投資建設而興起的東北老工業(yè)基地一度陷入市場困境，產(chǎn)生“東北現(xiàn)象”。對此，2003年中央頒布并實施“振興東北等老工業(yè)基地戰(zhàn)略”，在該戰(zhàn)略實施的10年內(nèi)，東北地區(qū)經(jīng)濟發(fā)展迅速，領(lǐng)先于全國經(jīng)濟平均水平。但近年來受體制和產(chǎn)業(yè)結(jié)構(gòu)影響，東北地區(qū)經(jīng)濟發(fā)展相對滯后，區(qū)域發(fā)展差異逐漸增大，同時出現(xiàn)建設用地無序開發(fā)現(xiàn)象加劇、耕地蠶食嚴重、土地利用效益低下等問題，被稱為“新東北現(xiàn)象”[20-22]。為恢復東北老工業(yè)基地的發(fā)展活力、推動產(chǎn)業(yè)轉(zhuǎn)型升級，中央于2016年實施新一輪東北振興戰(zhàn)略。如何協(xié)同推進經(jīng)濟穩(wěn)定復興、糧食長久安全、自然生態(tài)保護，成為東北地區(qū)亟待完成的重要任務。因此，優(yōu)化東北地區(qū)國土空間開發(fā)結(jié)構(gòu)，優(yōu)化配置各類空間資源，對促進東北地區(qū)按照主體功能區(qū)規(guī)劃發(fā)展、實現(xiàn)經(jīng)濟—社會—生態(tài)效益的聯(lián)動與統(tǒng)一至關(guān)重要[20]。

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

研究所需土地利用數(shù)據(jù)來源于中國科學院資源環(huán)境科學與數(shù)據(jù)中心基于Landsat遙感影像通過人工目視解譯生成的1 km土地利用柵格數(shù)據(jù)，利用ArcGIS 10.2對其進行矢量化處理，而后利用“面積制表”功能提取1990年、1995年、2000年、2005年、2010年、2015年、2018年7個研究時點東北地區(qū)各地市各地類面積。社會經(jīng)濟發(fā)展數(shù)據(jù)主要來源于中國縣域統(tǒng)計年鑒及相關(guān)省市統(tǒng)計年鑒。

1.3 研究方法

1.3.1 國土開發(fā)強度指數(shù)和泰爾指數(shù) 國土開發(fā)強度指數(shù)(LDI)(式(1))表征區(qū)域內(nèi)建設用地面積CLA占區(qū)域總面積TA的比例，反映區(qū)域內(nèi)用地結(jié)構(gòu)和對人口、經(jīng)濟活動的承載度[6]，狹義上指土地開發(fā)程度，側(cè)重對建設用地(本文指城鎮(zhèn)用地、農(nóng)村居民點等)的描述，注重建設用地的布局發(fā)展、強度約束及分區(qū)規(guī)劃；泰爾指數(shù)(式(2)-式(5))能將區(qū)域總差異分為組內(nèi)差異和組間差異，反映區(qū)域內(nèi)和區(qū)域間差異的演變趨勢和規(guī)律[23]。

LDI=CLA/TA

(1)

(2)

(3)

(4)

Za=Ra/R,Zb=Rb/R,Zj=(Cj/C)·(Rai/R)

(5)

式中：R為東北地區(qū)總差異；Rai為區(qū)域內(nèi)基本單元國土開發(fā)強度差異；Ra為組內(nèi)差異(省內(nèi)差異加權(quán)求和);Rb為組間差異；Za、Zb、Zj分別為區(qū)域內(nèi)差異、區(qū)域間差異和子區(qū)域差異對總差異的貢獻率；Ci、Cj分別為各市和各省建設用地面積；Cji、Tji分別為j省i市建設用地面積和土地總面積；Ti、Tj分別為各市和各省土地總面積；C、T分別為研究區(qū)域建設用地總面積和土地總面積；t為j省地市i的個數(shù)。

1.3.2 空間自相關(guān)分析 運用GeoDa軟件計算全局莫蘭指數(shù)(Moran′sI)(式(6))，以表征東北地區(qū)國土開發(fā)強度的空間相關(guān)性[24]，從整體上揭示國土開發(fā)布局的集聚情況[25]；引入局域莫蘭指數(shù)Ii(式(7)、式(8))對國土開發(fā)強度進行局域自相關(guān)分析[26]。

(6)

(7)

(8)

1.3.3 熱點分析 為顯示東北地區(qū)1990-2018年國土開發(fā)強度在空間上的高值集聚區(qū)和低值集聚區(qū)，本文運用Getis-OrdGi*方法模擬具有統(tǒng)計量顯著性的冷熱點空間集聚[27]。

1.3.4 地理探測器 地理探測器用于探測地理現(xiàn)象及其驅(qū)動因素空間分異性[28]。東北地區(qū)國土開發(fā)強度時空演變驅(qū)動因子涉及多方面，且作用規(guī)律復雜多變，因此，本文運用地理探測器分析不同驅(qū)動因素對國土開發(fā)強度的影響程度，計算公式如下：

(9)

1.3.5 國土開發(fā)強度驅(qū)動機制指標體系構(gòu)建 為探討東北地區(qū)國土開發(fā)強度時空演變驅(qū)動機制，從人口、經(jīng)濟、產(chǎn)業(yè)、投資、政策、生態(tài)六方面構(gòu)建指標體系(表1)。人口的遷移與集聚體現(xiàn)了社會經(jīng)濟發(fā)展對建設用地的需求，采用人口城鎮(zhèn)化率、人口密度表示人口對國土開發(fā)強度的驅(qū)動作用[29]；國土開發(fā)的主要目標是實現(xiàn)經(jīng)濟增長，采用地均國內(nèi)生產(chǎn)總值、地均財政收入表示經(jīng)濟對國土開發(fā)強度的影響[24]；產(chǎn)業(yè)結(jié)構(gòu)升級以及二、三產(chǎn)業(yè)產(chǎn)值增量驅(qū)動建設用地擴張，用于表示產(chǎn)業(yè)對國土開發(fā)強度的影響[11]；投資強度決定土地開發(fā)的產(chǎn)出效益，采用地均固定資產(chǎn)投資、地均財政支出表征投資對國土開發(fā)強度的影響[12]；土地利用政策的核心之一是保護耕地，保障永久基本農(nóng)田的面積，采用耕地比重表示政策對國土開發(fā)強度的影響[21]；區(qū)域自然資源稟賦與生態(tài)系統(tǒng)功能影響國土開發(fā)強度，選取綠地覆蓋率、水資源總量表征生態(tài)對國土開發(fā)強度的影響[30]。

表1 國土開發(fā)強度時空格局演變影響因素

2 結(jié)果分析

2.1 國土開發(fā)強度的時空演化及其區(qū)域差異

2.1.1 國土開發(fā)強度的時空演變趨勢 由圖2可知，1990-2018年東北地區(qū)國土開發(fā)強度整體呈穩(wěn)定增長態(tài)勢，2010-2018年增速變快，但整體國土開發(fā)強度仍低于全國同期水平。從空間演變看，40個地市的國土開發(fā)強度空間分異明顯，總體呈東南高、西北低趨勢。國土開發(fā)強度較高的區(qū)域逐漸由東南沿海的遼中南城市群向北部擴展至哈長城市群沿線，分布格局與國家戰(zhàn)略發(fā)展的核心區(qū)域吻合[19]，且東部沿海地區(qū)國土開發(fā)強度日益增加，而內(nèi)蒙古東四盟及黑龍江北部地區(qū)國土開發(fā)強度仍較低。

圖2 東北地區(qū)國土開發(fā)強度時空演變格局

2.1.2 國土開發(fā)強度差異測算及分解 利用泰爾指數(shù)計算1990年、1995年、2000年、2005年、2010年、2015年、2018年7個研究時點東北地區(qū)的總體差異以及組間、組內(nèi)、區(qū)域內(nèi)差異(圖3)。

(1)東北地區(qū)國土開發(fā)強度區(qū)域差異明顯。如圖3a所示，7個研究時點的總體泰爾指數(shù)(R)平均值為0.41，說明東北地區(qū)國土開發(fā)強度差異明顯。從時空層面上看，1990-2018年國土開發(fā)強度的總體泰爾指數(shù)變化相對平穩(wěn)，說明伴隨著新型城鎮(zhèn)化建設以及東北振興戰(zhàn)略的實施，東北地區(qū)國土開發(fā)強度差異未發(fā)生明顯變化。東北地區(qū)組內(nèi)差異(Ra)和組間差異(Rb)均較大，且在1990-2015年未發(fā)生波動，二者分別為0.22、0.19，對總差異的貢獻率分別為54%、46%，說明期間東北地區(qū)省間及省內(nèi)國土開發(fā)強度差異均較大；但2015-2018年組內(nèi)差異下降而組間差異上升，且二者變化幅度大致相同，導致組間差異大于組內(nèi)差異，說明期間黑、吉、遼及內(nèi)蒙古東四盟四區(qū)域間的國土開發(fā)強度差異是總差異的主要來源，四區(qū)域間非均衡發(fā)展趨勢加重。

(2)區(qū)域內(nèi)國土開發(fā)強度極不平衡。如圖3b所示，黑、吉、遼及內(nèi)蒙古東四盟四區(qū)域內(nèi)部的泰爾指數(shù)在1990-2015年變化相對平穩(wěn)，在2018年略有下降，說明四區(qū)域內(nèi)部國土開發(fā)強度差異略有收斂，這與組內(nèi)差異在2018年變小的趨勢相吻合。四區(qū)域相比較可知，遼寧的泰爾指數(shù)最小，對總體泰爾指數(shù)的貢獻率較小，省內(nèi)各地市國土開發(fā)強度趨于均衡；吉林的泰爾指數(shù)最大，與黑龍江和內(nèi)蒙古東四盟相近，說明三區(qū)域各地市國土開發(fā)強度差異較大，內(nèi)部發(fā)展極不平衡。

圖3 東北地區(qū)國土開發(fā)強度泰爾指數(shù)演變趨勢

2.2 國土開發(fā)強度的空間自相關(guān)分析

2.2.1 全域空間自相關(guān)分析 7個研究時點的全局莫蘭指數(shù)值(0.55～0.58)通過了1%的顯著性檢驗(表2)，說明東北地區(qū)國土開發(fā)強度表現(xiàn)出極顯著的空間集聚特征。2015-2018年全局莫蘭指數(shù)呈上升趨勢，說明隨著經(jīng)濟發(fā)展，東北地區(qū)國土開發(fā)強度的空間相關(guān)性增強，高值區(qū)和低值區(qū)更加集聚，各地市之間國土開發(fā)強度的相互影響及聯(lián)動效應增強，這也與圖3a組內(nèi)差異縮小、組間差異增大的結(jié)果相吻合。綜合看，東北地區(qū)區(qū)域間國土開發(fā)強度差異一直存在并可能進一步擴大，空間發(fā)展失衡現(xiàn)象將長期存在。統(tǒng)計分析7個研究時點東北地區(qū)國土開發(fā)強度Moran′sI散點圖每個象限的地市數(shù)量可知，L-L地市數(shù)量(2018年21個，其他年份均為20個)明顯多于其他象限，占東北地區(qū)總地市數(shù)量的1/2；H-H地市數(shù)量一直為14，而L-H(2015年為6個，其他年份均為5個)和H-L(1990-2010年為1個，2010年后減少為0)地市數(shù)量較少，表明東北地區(qū)國土開發(fā)強度較低的地市較多，高值區(qū)和低值區(qū)集聚現(xiàn)象明顯，即總體差異明顯，這與泰爾指數(shù)所得結(jié)論吻合。兩輪東北振興戰(zhàn)略實施后，東北地區(qū)地市間經(jīng)濟、技術(shù)交流增多，但并未實現(xiàn)區(qū)域全面振興，社會發(fā)展格局依然固化，經(jīng)濟基礎(chǔ)薄弱地市依然相對落后，這說明東北地區(qū)經(jīng)濟發(fā)展失衡的大格局仍未改變，實現(xiàn)協(xié)調(diào)發(fā)展任重道遠。

表2 東北地區(qū)國土開發(fā)強度的全局莫蘭指數(shù)

2.2.2 局域空間自相關(guān)分析 為更好地比較東北振興戰(zhàn)略實施后國土開發(fā)強度空間相關(guān)性變化格局，利用GeoDa軟件繪制7個研究時點國土開發(fā)強度LISA圖(圖4)，據(jù)此分析東北地區(qū)國土開發(fā)強度的局域相關(guān)性。由圖4可知，東北地區(qū)國土開發(fā)強度表現(xiàn)出較明顯的南北空間分異格局，中部城市空間相關(guān)性不顯著。低值區(qū)主要集中在黑龍江省北部大興安嶺—黑河—伊春—鶴崗—佳木斯—七臺河一帶以及內(nèi)蒙古呼倫貝爾等城市，高值區(qū)主要集中在遼寧省南部錦州—盤錦—鞍山—營口—大連一帶。第一輪東北振興戰(zhàn)略實施后，低值集聚區(qū)向東移動，所含城市變多。2013年“新東北現(xiàn)象”出現(xiàn)后，2015年低值集聚區(qū)擴張至東北地區(qū)北部，國土開發(fā)強度空間差異進一步擴大，區(qū)域發(fā)展極度失衡。針對此現(xiàn)象，2016年實施新一輪東北振興戰(zhàn)略，2018年低值集聚區(qū)有所減少，但區(qū)域發(fā)展不平衡問題依然存在。

圖4 東北地區(qū)國土開發(fā)強度LISA圖

2.3 熱點分析

利用ArcGIS 10.2軟件熱點分析功能，從90%、95%、99% 3個置信區(qū)間分析7個研究時點東北地區(qū)國土開發(fā)強度的冷熱點空間分布格局及其演變趨勢，觀察高值區(qū)與低值區(qū)的集聚位置，制作冷熱點空間分布圖(圖5)?？梢钥闯?，東北地區(qū)國土開發(fā)強度的熱點區(qū)域主要集中在沈陽、錦州、遼陽、盤錦、鞍山、營口、大連等遼中南地區(qū)，截至2018年，該地區(qū)土地總面積為5.84萬km2，占東北地區(qū)土地總面積的5%，建設用地面積為7 476 km2，占東北地區(qū)建設用地總面積的20%，經(jīng)濟總量占東北地區(qū)的30%，說明省會和沿海城市及其周邊城市經(jīng)濟發(fā)展狀況相對較好，建設用地開發(fā)程度較高，是東北地區(qū)經(jīng)濟發(fā)展的核心區(qū)域。東北地區(qū)國土開發(fā)強度的冷點區(qū)域有擴散趨勢，主要由黑龍江北部的大興安嶺、黑河、伊春、佳木斯及內(nèi)蒙古呼倫貝爾等城市擴散至七臺河、齊齊哈爾及延邊州，這些地區(qū)土地利用類型以耕地、林地、草地為主，2018年土地總面積為54.07萬km2，占東北地區(qū)土地總面積的44%，是熱點地區(qū)的9.3倍，但建設用地面積僅為6 347 km2，占東北地區(qū)建設用地總面積的17%，經(jīng)濟總量僅占東北地區(qū)的8.7%，遠低于熱點地區(qū)，這說明經(jīng)濟發(fā)展落后地區(qū)建設用地開發(fā)程度也較低，東北地區(qū)國土開發(fā)強度空間分異格局與其區(qū)域經(jīng)濟空間分異格局一致。

圖5 東北地區(qū)國土開發(fā)強度熱點分析

2.4 驅(qū)動機制分析

運用地理探測器中的“因子探測器”對2000-2018年東北地區(qū)40個地市指標變量進行探測，所得結(jié)果均在1%水平上顯著，并根據(jù)每項指標的權(quán)重計算出人口、經(jīng)濟、產(chǎn)業(yè)、投資、政策、生態(tài)6項要素對國土開發(fā)強度驅(qū)動效應的變化(表3)。整體看，2000年人口、經(jīng)濟、投資及政策是影響東北地區(qū)國土空間開發(fā)格局的主導因素，而產(chǎn)業(yè)和生態(tài)的影響相對較弱；2018年除政策要素外，其余要素對東北地區(qū)國土空間開發(fā)格局的驅(qū)動力均在波動中保持平穩(wěn)或增加趨勢，尤以經(jīng)濟、產(chǎn)業(yè)兩要素的變化顯著。

表3 東北地區(qū)不同指標因子對國土開發(fā)強度的影響探測結(jié)果

2.4.1 人口、生態(tài)對國土開發(fā)強度的影響 研究期內(nèi)人口和生態(tài)要素對東北地區(qū)國土開發(fā)強度的驅(qū)動效應均呈現(xiàn)先增后減狀態(tài)，變化幅度較小。東北地區(qū)40個地市平均城鎮(zhèn)化率在2000年達到0.61，2005年降至0.48，此后出現(xiàn)人口外流現(xiàn)象，建設用地需求相對減弱，所以人口對建設用地的驅(qū)動作用呈回落趨勢。生態(tài)資源條件較好地區(qū)建設用地擴張是以犧牲林地、草地等為代價，該區(qū)域建設用地擴張受到一定程度的限制，所以生態(tài)因素對東北地區(qū)國土開發(fā)強度的驅(qū)動作用相對穩(wěn)定。

2.4.2 經(jīng)濟發(fā)展水平對國土開發(fā)強度的影響 經(jīng)濟發(fā)展水平對東北地區(qū)國土開發(fā)強度的影響不斷增大，成為主導因素。經(jīng)濟發(fā)展水平較高的地區(qū)，生產(chǎn)活動及經(jīng)濟開發(fā)區(qū)建設推動其他類型土地轉(zhuǎn)向建設用地，此外，政府通過土地出讓等方式增加財政收入，間接擴大了建設用地規(guī)模。

2.4.3 投資對國土開發(fā)強度的影響 投資對東北地區(qū)國土開發(fā)強度的影響較大，在研究期內(nèi)呈先升后降的波動變化趨勢。單位面積固定資產(chǎn)投資越高，建設用地開發(fā)利用程度就越高，而政府為進一步擴大經(jīng)濟發(fā)展規(guī)模，會加大對區(qū)域財政投資力度，使區(qū)域?qū)θ丝?、企業(yè)、資本等要素的吸引和集聚作用加強，增大了建設用地需求，導致建設用地開發(fā)強度進一步增加[12]。但東北地區(qū)體制較為僵化，對長期單一的財力刺激敏感度大幅下降，導致投資對國土開發(fā)強度的驅(qū)動效應減弱，后經(jīng)改革雖有所好轉(zhuǎn)，但并未根治。

2.4.4 產(chǎn)業(yè)對國土開發(fā)強度的影響 產(chǎn)業(yè)對東北地區(qū)國土開發(fā)強度的影響較小，但整體呈上升趨勢。二、三產(chǎn)業(yè)發(fā)展帶來更多就業(yè)機會，促進人口集聚，從而增加居住、商服、交通用地需求，促進建設用地開發(fā)；同時，產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化也會提高城市土地開發(fā)強度[13]。

2.4.5 政策對國土開發(fā)強度的影響 政策對東北地區(qū)國土開發(fā)強度的影響維持在較高水平，但在研究期內(nèi)呈下降趨勢，特別是土地利用政策對東北地區(qū)國土開發(fā)強度的影響逐漸減弱。東北地區(qū)作為我國糧食主產(chǎn)地，承擔著維護國家糧食安全的重任，農(nóng)業(yè)現(xiàn)代化發(fā)展需要大量的勞動力和資金投入，限制了農(nóng)村人口外流，從而增加了農(nóng)村建設用地需求[21]；同時，耕地的規(guī)?；?、集約化發(fā)展對基礎(chǔ)設施配套建設需求也不斷提高。因此，土地利用政策對東北地區(qū)國土開發(fā)強度的驅(qū)動力較高。但在2008年金融危機發(fā)生后，東北地區(qū)為促進經(jīng)濟發(fā)展，耕地保護壓力相對變小，土地利用政策對國土開發(fā)強度的影響隨之下降。此外，為緩解“東北現(xiàn)象”帶來的負面影響，政府先后圈定遼寧沿海、長吉圖及哈長城市群作為重點培育區(qū)域，也促進了東北地區(qū)建設用地擴張。

總之，國土開發(fā)強度的提升主要表現(xiàn)為建設用地擴張，而建設用地擴張及其空間格局演變是多種因素綜合作用的結(jié)果，驅(qū)動因子具有多元性和復雜性。隨著人口、經(jīng)濟、產(chǎn)業(yè)的發(fā)展，對建設用地的需求不斷加大，刺激建設用地擴張。對土地開發(fā)的投資成為建設用地擴張的資本動力，會吸引大量人口和企業(yè)集聚；政府根據(jù)區(qū)域特征差異所制定的主體功能區(qū)規(guī)劃以及頒布的相關(guān)政策，是建設用地合理擴張的依據(jù)；良好的自然資源本底條件及生態(tài)環(huán)境降低了建設用地的擴張成本，促進了建設用地開發(fā)[31]。因此，適宜的生態(tài)環(huán)境、雄厚的資金及政策的合理引導為滿足建設用地需求提供了保障(圖6)。

圖6 東北地區(qū)國土開發(fā)驅(qū)動機制

3 結(jié)論與討論

本研究基于土地利用數(shù)據(jù)，利用泰爾指數(shù)、莫蘭指數(shù)、熱點分析及地理探測器方法，分析了1990-2018年東北地區(qū)國土開發(fā)強度的時空格局及其驅(qū)動機制，結(jié)果表明：1)東北地區(qū)國土開發(fā)強度總體差異較大，發(fā)展不平衡不充分問題顯著，組間差異和組內(nèi)差異均較大，且四區(qū)域的組間差異有上升趨勢。從區(qū)域內(nèi)部看，遼寧各地市之間國土開發(fā)強度差異相對較小，而黑龍江、吉林、內(nèi)蒙古東四盟內(nèi)部各地市之間差異較大。2)1990-2018年東北地區(qū)國土開發(fā)強度總體呈增長態(tài)勢，但增長速度相對緩慢。整體呈東南高、西北低的態(tài)勢，低值區(qū)占比較多，主要集中在黑龍江、內(nèi)蒙古東四盟的北部，而高值區(qū)主要集中在遼中南城市群和哈長城市群，這與國家發(fā)展戰(zhàn)略布局相吻合。高值區(qū)與低值區(qū)的集聚性呈增強趨勢，國土開發(fā)強度空間分布的非均衡趨勢更顯著，這與東北地區(qū)經(jīng)濟發(fā)展的區(qū)域差異格局相一致。3)人口、生態(tài)對東北地區(qū)國土開發(fā)強度的驅(qū)動力在波動中保持平穩(wěn)，政策對國土開發(fā)強度的驅(qū)動力處于較高水平，但在持續(xù)下降，經(jīng)濟、產(chǎn)業(yè)、投資對國土開發(fā)強度的驅(qū)動力不斷上升，尤以經(jīng)濟、產(chǎn)業(yè)的驅(qū)動力上升趨勢顯著。總體看，東北地區(qū)區(qū)域間經(jīng)濟發(fā)展不平衡不協(xié)調(diào)是國土開發(fā)強度不均衡及開發(fā)利用效益不高的主要原因，研究結(jié)果可在理論和實踐上為優(yōu)化東北地區(qū)國土空間開發(fā)格局、促進資源要素有效流動、激活要素生產(chǎn)力提供參考。

改革開放以來，在各種經(jīng)濟發(fā)展機遇下，東北地區(qū)表現(xiàn)出的遲鈍、固化、發(fā)展活力不足致使其地域空間結(jié)構(gòu)僵化，經(jīng)濟格局和國土開發(fā)強度空間格局基本保持穩(wěn)定[19]。如何進一步活化要素資源配置，打破“區(qū)域—要素”流通限制，激發(fā)內(nèi)生活力，成為下一步的研究重點。東北地區(qū)作為中國重要的老工業(yè)基地、商品糧主要產(chǎn)地以及國家重要的生態(tài)安全屏障，要以東北振興為機遇，加快工業(yè)化、城鎮(zhèn)化及農(nóng)業(yè)現(xiàn)代化建設步閥[32]，進行合理的主體功能區(qū)定位，緩解城鎮(zhèn)開發(fā)與耕地保護之間的矛盾[33]，打造國土空間開發(fā)保護新格局。對農(nóng)產(chǎn)品主產(chǎn)區(qū)要加強財政轉(zhuǎn)移支付力度，勇?lián)鷩壹Z食安全“壓艙石”的重任；要把握區(qū)域相對優(yōu)勢，優(yōu)化經(jīng)濟發(fā)展路徑，引進科技發(fā)展人才，增強科技創(chuàng)新能力，借“一帶一路”規(guī)劃做好向北開放的窗口，實現(xiàn)開放型經(jīng)濟轉(zhuǎn)型；要改革僵化的體制機制，擴大就業(yè)，提高基礎(chǔ)設施建設水平與公共服務供給能力，消除資源枯竭型城市的貧困，在生態(tài)功能區(qū)加強生態(tài)環(huán)境保護和修復力度，努力實現(xiàn)東北地區(qū)經(jīng)濟—社會—生態(tài)效益的聯(lián)動發(fā)展。