基于地理探測器的山東省土地城鎮(zhèn)化時空演變及驅(qū)動因素

劉慶果，張煊赟，莫鈞雄，孫丕苓，3

（1.曲阜師范大學(xué) 地理與旅游學(xué)院，山東 日照 276826；2.重慶師范大學(xué) 地理與旅游學(xué)院，重慶 401331；3.日照市國土空間規(guī)劃與生態(tài)建設(shè)重點(diǎn)實(shí)驗(yàn)室，山東 日照 276826）

0 引言

城鎮(zhèn)化是伴隨工業(yè)化發(fā)展，非農(nóng)產(chǎn)業(yè)在城鎮(zhèn)集聚和農(nóng)村人口向城鎮(zhèn)集中的自然歷史過程，是人類社會發(fā)展的客觀趨勢，是國家現(xiàn)代化的重要標(biāo)志。其中，土地城鎮(zhèn)化是推動城市邁向現(xiàn)代化發(fā)展的空間載體［1］。改革開放以來，區(qū)域土地城鎮(zhèn)化持續(xù)超前發(fā)展和無序蔓延，國土空間格局發(fā)生了劇烈變化，主要表現(xiàn)為建設(shè)空間持續(xù)快速擴(kuò)張，自然生態(tài)空間不斷萎縮，農(nóng)業(yè)耕作空間格局不斷發(fā)生重構(gòu)［2］。進(jìn)入21 世紀(jì)，我國工業(yè)化和城鎮(zhèn)化進(jìn)程加快，城鎮(zhèn)建設(shè)用地快速擴(kuò)張，城鎮(zhèn)人口增加了近3.73 億人，人口城鎮(zhèn)化率提高到59.58%。城鎮(zhèn)化作為土地利用變化的重要表現(xiàn)形式，推動農(nóng)村人口向城鎮(zhèn)地區(qū)集聚［3-5］。因此，厘清土地城鎮(zhèn)化時空格局及驅(qū)動因素對促進(jìn)我國新型城鎮(zhèn)化建設(shè)具有重要意義。

目前，學(xué)術(shù)界對土地城鎮(zhèn)化開展了大量研究，并取得了豐富的研究成果。眾多學(xué)者基于不同的學(xué)科背景，從不同視角開展了土地城鎮(zhèn)化研究，相關(guān)研究主要集中于土地城鎮(zhèn)化水平測度［3，6-9］、土地城鎮(zhèn)化時空格局變化及驅(qū)動因素［10-14］，以及土地城鎮(zhèn)化與人口［15-19］、經(jīng)濟(jì)［20，21］、環(huán)境［22-25］等要素的耦 合協(xié)調(diào)關(guān)系等方面，研究尺度涵蓋國家、區(qū)域、省域、市域、縣域 等行政 單 元［1，5，7，18，26］。由 于 縣 域 是 土 地 利 用 和管理最直接的操作層，是我國土地資源利用、管理和規(guī)劃的最佳尺度，縣域單元成為土地城鎮(zhèn)化研究重點(diǎn)關(guān)注的研究尺度［27］。已有的土地城鎮(zhèn)化研究主要采用綜合指數(shù)法、面積比例法、探索性空間數(shù)據(jù)分析、相關(guān)分析、主成分分析、偏最小二乘法、空間回歸、地理加權(quán)回歸、耦合協(xié)調(diào)模型等方法模型進(jìn)行土地城 鎮(zhèn) 化 水 平 測 度 及 影 響 因 素 分 析［1，14，21，28，29］。其中，綜合指數(shù)法、面積比例法、相關(guān)分析、主成分分析等方法操作較為簡單，應(yīng)用較為廣泛。研究熱點(diǎn)區(qū)涵蓋長三角地區(qū)、東北地區(qū)、中部地區(qū)和環(huán)渤海地區(qū)等經(jīng)濟(jì)發(fā)展較好的地區(qū)。如：王富喜［16］綜合運(yùn)用加權(quán)求和法、彈性系數(shù)法和耦合協(xié)調(diào)度模型對山東半島城市群人口—土地城鎮(zhèn)化質(zhì)量進(jìn)行了測度與協(xié)調(diào)發(fā)展研究；張凱煌等［14］運(yùn)用熵值法和多層線性模型剖析了中國城市土地城鎮(zhèn)化及多層級影響因素；楊喜等［1］基于探索性空間數(shù)據(jù)分析法分析了我國中部地區(qū)土地城鎮(zhèn)化時空特征及影響因素；俞振寧等［28］運(yùn)用ESDA—GWR方法探究了浙江省土地城鎮(zhèn)化空間特征及影響因素問題。綜上，土地城鎮(zhèn)化時空變化的影響因素研究大多以線性分析和相關(guān)分析為主，而對土地城鎮(zhèn)化時空變化影響因素的空間信息關(guān)注相對較少。

本文以山東省為例，運(yùn)用全局自相關(guān)分析和地理探測器模型探討了2000—2020 年山東省土地城鎮(zhèn)化時空格局演變及驅(qū)動因素，以期為該省新型城鎮(zhèn)化建設(shè)提供參考。

1 研究區(qū)概況

山東位于我國黃河下游的華北平原地區(qū)，地處34°22′—38°24′N、114°47′—122°42′E之間，東臨渤海、黃海，與朝鮮半島、日本列島隔海相望，自北向南與河北、河南、安徽、江蘇等省份毗鄰，下轄17 個地市，包括137 個縣（市、區(qū)），土地總面積156367.73km2（圖1）。境內(nèi)地勢中部高、四周低，西南和西北部為低洼平坦，東部和東南部為緩丘起伏。全區(qū)地處溫帶季風(fēng)氣候區(qū)，年平均降水量為550—950mm，年均氣溫為11—14℃，雨熱同期，是我國重要的農(nóng)業(yè)生產(chǎn)基地。2020 年底，山東省總?cè)丝跒?0152.70 萬人，實(shí)現(xiàn)GDP 73129 億元。

圖1 山東省研究區(qū)位Figure 1 Location of Shandong Province

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

2.1 數(shù)據(jù)來源與處理

本文采用的數(shù)據(jù)主要包括城鎮(zhèn)建設(shè)用地和社會經(jīng)濟(jì)數(shù)據(jù)。其中：城鎮(zhèn)建設(shè)用地（含城市用地、建制鎮(zhèn)用地）通過地理空間數(shù)據(jù)云（http：//www.gscloud.cn/）平臺的2000 年、2010 年、2018 年Landsat TM/OLI遙感影像獲取，Kappa系數(shù)超過0.80，能夠滿足研究的需要；社會經(jīng)濟(jì)數(shù)據(jù)主要來源于對應(yīng)年份的《中國縣域統(tǒng)計(jì)年鑒》《山東省統(tǒng)計(jì)年鑒》和各地市統(tǒng)計(jì)年鑒、社會經(jīng)濟(jì)統(tǒng)計(jì)公報(bào)等，部分缺失數(shù)據(jù)采用鄰近年份數(shù)據(jù)代替。為保持統(tǒng)計(jì)口徑的一致性，對GDP、社會固定資產(chǎn)投資，社會消費(fèi)品零售總額，第二、三產(chǎn)業(yè)產(chǎn)值，財(cái)政收入，財(cái)政支出，實(shí)際利用外商投資，出口額等進(jìn)行平滑處理，以消除價(jià)格因素的影響，確保數(shù)據(jù)的可比性［20］?？紤]到山東省行政區(qū)劃的變遷，需要進(jìn)行行政區(qū)域調(diào)整處理。本文以2018 年山東省縣域行政區(qū)劃為基礎(chǔ)，對研究期內(nèi)山東省行政區(qū)劃進(jìn)行調(diào)整。

2.2 研究方法

土地城鎮(zhèn)化水平測度：運(yùn)用縣域城鎮(zhèn)建設(shè)用地面積與縣域總面積的比值表征土地城鎮(zhèn)化水平。計(jì)算公式為［3］：

式中：LU 表示土地城鎮(zhèn)化水平；Su表示城鎮(zhèn)建設(shè)用地面積；Sa表示總面積。

空間自相關(guān)分析：空間自相關(guān)分析是檢驗(yàn)?zāi)骋豢臻g要素或某一空間要素屬性值的空間分布狀況與其臨近空間區(qū)域之間關(guān)聯(lián)程度的方法。本文運(yùn)用全局空間自相關(guān)分析方法（Global Moran′s I 指數(shù)）檢驗(yàn)山東省縣域土地城鎮(zhèn)化水平在空間上是否聚集。計(jì)算公式為［28，29］：

式中：n 為山東省縣區(qū)數(shù)量；Xi與Xj分別表示縣區(qū)i、j 的土地城鎮(zhèn)化水平；Wij為空間權(quán)重矩陣；S2為土地城鎮(zhèn)化水平的方差。Moran′s I 的取值范圍為［-1，1］。若值為正，表明山東省縣區(qū)土地城鎮(zhèn)化呈現(xiàn)空間聚集；若值為負(fù)，表明土地城鎮(zhèn)化呈離散分布；若值為0，則為隨機(jī)分布。

采用局部空間自相關(guān)識別山東省土地城鎮(zhèn)化的空間依賴性和異質(zhì)性，計(jì)算公式為［30］：

影響因子指標(biāo)選?。和恋爻擎?zhèn)化水平的空間格局及演變受到地理區(qū)位、社會經(jīng)濟(jì)因素和政策因素的綜合影響。鑒于區(qū)域政策體制數(shù)據(jù)難以量化，參考已有研究［1，7，18，28，29］，結(jié)合土地城鎮(zhèn)化發(fā)展實(shí)際，從地理區(qū)位、社會因素、經(jīng)濟(jì)因素等方面選擇山東省縣域土地城鎮(zhèn)化時空差異的解釋因子（表1）。本文以土地城鎮(zhèn)化率為因變量，選取距最近地市中心距離、距最近交通道路距離表征地理區(qū)位條件；社會因素方面，選擇城鎮(zhèn)人口密度表征人口規(guī)模；經(jīng)濟(jì)因素方面，以人均GDP、地均固定資產(chǎn)投資額、人均社會消費(fèi)品零售額表征經(jīng)濟(jì)發(fā)展水平，以第二、三產(chǎn)業(yè)比重表征區(qū)域產(chǎn)業(yè)結(jié)構(gòu)，以地均地方財(cái)政收入、地均地方財(cái)政支出表征政府調(diào)控能力，以地均工業(yè)產(chǎn)值表征工業(yè)化水平，以地均外商投資、地均出口額表征開放水平。

表1 土地城鎮(zhèn)化的驅(qū)動因素及其指標(biāo)體系Table 1 Driving factor and its index system of land urbanization

（續(xù)表1）

地理探測器模型：運(yùn)用地理探測器的因子探測和交互作用探測模塊探測山東省縣域土地城鎮(zhèn)化的空間分異特性，并探測各驅(qū)動因素對土地城鎮(zhèn)化空間分異的影響程度及各因子交互作用的影響程度。計(jì)算公式為［30，31］：

式中：q 為土地城鎮(zhèn)化時空格局變化的決定力指標(biāo)，取值范圍為［0，1］，q 值越大，表明驅(qū)動因子對土地城鎮(zhèn)化時空格局變化的影響程度越大；h 為驅(qū)動因子的分類數(shù)；N 為縣域單元數(shù)；σ2為各縣域單元土地城鎮(zhèn)化率的方差。交互探測類型可分為雙因子增強(qiáng)、線性增強(qiáng)、非線性減弱、單因子非線性減弱和獨(dú)立5 類。

3 結(jié)果與分析

3.1 時空變化特征

根據(jù)山東省縣域土城鎮(zhèn)化水平的綜合測度結(jié)果，運(yùn)用ArcGIS自然斷裂點(diǎn)法將其劃分高、較高、中等、較低、低5 個等級（圖2）。2000—2020 年，山東省土地城鎮(zhèn)化水平的等級分布呈現(xiàn)兩極占比低而中間占比高的格局，不同等級城鎮(zhèn)化水平差異顯著。2000 年，山東省土地城鎮(zhèn)化水平總體處于較低水平，土地城鎮(zhèn)化處于高水平等級的縣區(qū)占全省縣區(qū)總數(shù)的5.11%，空間分布較為分散，主要分布于青島市南區(qū)、市北區(qū)、李滄區(qū)和煙臺芝罘區(qū)、濟(jì)南歷下區(qū)、濰坊奎文區(qū)、淄博張店區(qū)等地；處于較高水平等級的縣區(qū)占13.14%，主要分布于濟(jì)南、淄博、濰坊、臨沂、濟(jì)寧、聊城等地市的市區(qū)和周邊縣域；處于低水平等級的縣區(qū)趨于山區(qū)、省域與市域邊緣分布，主要分布泰沂山區(qū)、沂蒙山區(qū)、半島丘陵地區(qū)和黃河三角洲地區(qū)。2010 年山東省土地城鎮(zhèn)化水平得到一定程度的提高，土地城鎮(zhèn)化高水平區(qū)向南部市區(qū)擴(kuò)展，臨沂羅莊區(qū)躍居高水平等級；較高水平區(qū)向西北部擴(kuò)展，城陽區(qū)、萊山區(qū)、環(huán)翠區(qū)、濱城區(qū)、高唐縣等縣區(qū)分別由中等水平提高到較高水平；較低水平及以下等級區(qū)范圍大幅縮減，呈現(xiàn)趨于沂蒙山區(qū)、半島丘陵、黃河三角洲地區(qū)集中分布的態(tài)勢。2020 年，山東省土地城鎮(zhèn)化水平得到大幅提升，處于高水平等級的縣區(qū)數(shù)量增加到19 個，高水平區(qū)向西北部、南部擴(kuò)展，德州德城區(qū)、棗莊市中區(qū)和臨沂河?xùn)|區(qū)、蘭山區(qū)等11 個地市轄區(qū)躍升高水平等級。較高水平區(qū)在聊城境內(nèi)、濟(jì)寧—菏澤和濟(jì)南—淄博—東營交界地區(qū)集中連片分布；較低及以下水平等級范圍大幅縮減；低水平區(qū)則趨于省域邊緣和沂蒙山區(qū)集中分布。

圖2 2000—2020 年山東省土地城鎮(zhèn)化水平空間分布Figure 2 Spatial distribution of land urbanization in Shandng Province，2000—2020

3.2 時空格局演變特征

本文基于GeoDa095i軟件計(jì)算出山東土地城鎮(zhèn)化的Moran′s I 指數(shù)及檢驗(yàn)結(jié)果。山東省土地城鎮(zhèn)化全局空間自相關(guān)檢驗(yàn)結(jié)果通過了顯著性檢驗(yàn)（Z >1.96，P ＜0.01），表明該省土地城鎮(zhèn)化水平存在明顯的空間關(guān)聯(lián)性。2000—2020 年Moran′s I 指數(shù)由0.455 減小到0.320，呈下降趨勢，表明山東省土地城鎮(zhèn)化空間集聚態(tài)勢逐漸減弱。

本文運(yùn)用局部空間自相關(guān)模塊計(jì)算2000—2020年山東省土地城鎮(zhèn)化局域G*i指數(shù)，并繪制研究區(qū)土地城鎮(zhèn)化水平空間集聚圖（圖3）。從圖3 可見，2000—2020 年山東省土地城鎮(zhèn)化空間極化現(xiàn)象顯著，高—高集聚區(qū)不斷擴(kuò)展，并由青島市區(qū)轉(zhuǎn)至濟(jì)南市區(qū)；低—低集聚區(qū)持續(xù)收縮，趨于沂蒙山區(qū)和半島丘陵地區(qū)集中分布。2000 年，土地城鎮(zhèn)化高—高、高—低集聚區(qū)趨于市區(qū)分布，其中：高—高集聚區(qū)分布于青島城陽區(qū)，高—低集聚區(qū)分布于環(huán)翠區(qū)和泰山區(qū)，低—低集聚區(qū)趨于沂蒙山區(qū)、半島丘陵地區(qū)和黃河三角洲集中分布。2020 年，土地城鎮(zhèn)化高值集聚區(qū)分布于濟(jì)南市中區(qū)、天橋區(qū)和歷城區(qū)，低值區(qū)分別于沂蒙山區(qū)、半島丘陵地區(qū)形成兩大低低集聚區(qū)，涵蓋平邑縣、蒙陰縣、沂水縣、沂源縣、萊西市、海陽市、乳山市、棲霞市等17 縣市，其他低—低集聚區(qū)零散分布于市域邊緣的臺兒莊區(qū)、臨邑縣和河口區(qū)。

圖3 2000—2020 年山東省土地城鎮(zhèn)化空間集聚Figure 3 Spatial pattern of land urbanization in Shandong Province，2000-2020

3.3 驅(qū)動因素分析

因子探測結(jié)果：本文運(yùn)用地理探測器模型識別 山東省土地城鎮(zhèn)化時空分異的關(guān)鍵影響因子，結(jié)果如表2 所示。從表2 可見，地理區(qū)位條件、社會因素、經(jīng)濟(jì)因素等方面的因子共同影響了山東省土地城鎮(zhèn)化的時空分異特征。2000 年，各解釋變量中地均固定資產(chǎn)投資的影響力居首位，表明固定資產(chǎn)投資是促進(jìn)土地城鎮(zhèn)化時空分異的主要驅(qū)動力。城鎮(zhèn)人口密度、地均地方財(cái)政收入和地均地方財(cái)政支出分列第二、三、四位，其q 值分別為0.377、0.338、0.313，超過0.3 且通過1%顯著性水平檢驗(yàn)，表明人口規(guī)模與政府調(diào)控能力對土地城鎮(zhèn)化時空分異具有重要驅(qū)動作用。此外，距最近地市距離，第二、三產(chǎn)業(yè)比重，地均外商投資，地均出口額，人均GDP，地均工業(yè)產(chǎn)值的q 值分別為0.295、0.272、0.250、0.189、0.164、0.116，通過1%顯著性水平檢驗(yàn)，表明經(jīng)濟(jì)區(qū)位、產(chǎn)業(yè)結(jié)構(gòu)、開放水平、居民收入水平和工業(yè)發(fā)展水平對土地城鎮(zhèn)化時空分異具有一定影響。地均社會消費(fèi)品零售額和距最近交通道路距離解釋能力較弱，未通過顯著性檢驗(yàn)。2010 年，各解釋變量中城鎮(zhèn)人口密度的影響力躍居首位，其 q 值提升至0.512。地均地方財(cái)政支出、地均固定資產(chǎn)投資額、地均地方財(cái)政收入分居第二、三、四位，其q 值分別為0.427、0.390、0.389，通過1%顯著性水平檢驗(yàn)，表明隨著政府調(diào)控能力的增強(qiáng)和投資水平的提高，促進(jìn)了城鎮(zhèn)建設(shè)用地的快速擴(kuò)張。此外，距最近地市距離，第二、三產(chǎn)業(yè)比重，地均外商投資額，地均出口額，地均工業(yè)產(chǎn)值，人均社會消費(fèi)品零售額等解釋變量通過1%顯著性水平檢驗(yàn)，而距最近交通道路距離解釋能力較弱，未通過顯著性檢驗(yàn)。2020 年各解釋變量中城鎮(zhèn)人口密度的影響強(qiáng)度大于其他解釋變量，其q 值提升至0.583。隨著城鎮(zhèn)化進(jìn)程中教育、醫(yī)療等公共服務(wù)設(shè)施的完善，吸引了大量鄉(xiāng)村人口流向城市，以尋求更好的就業(yè)、教育和醫(yī)療服務(wù)，人口規(guī)模的增長帶動了城鎮(zhèn)擴(kuò)張。第二、三產(chǎn)業(yè)比重，地均地方財(cái)政收入，人均社會消費(fèi)品零售額的影響力度躍居第二至四位，其q 值分別為0.532、0.526、0.468，通過1%顯著性水平檢驗(yàn)，表明產(chǎn)業(yè)結(jié)構(gòu)升級、政府調(diào)控能力提升和居民消費(fèi)水平提高是拉動城鎮(zhèn)化進(jìn)程的重要動力。距最近地市距離、地均地方財(cái)政支出、地均出口額、地均固定資產(chǎn)投資額、地均外商投資、地均工業(yè)產(chǎn)值等驅(qū)動因子通過了1%顯著性水平檢驗(yàn)，距最近地市距離、人均社會消費(fèi)品零售額、地均出口額的影響強(qiáng)度增長較快，表明近10年山東省城市輻射帶動、居民消費(fèi)與出口拉動對促進(jìn)土地城鎮(zhèn)化水平提升作用顯著。

表2 山東省土地城鎮(zhèn)化時空分異的驅(qū)動因素探測結(jié)果Table 2 Detection results of driving factors of land urbanization spatial- temporal differentiation in Shandong Province

總體來看：2000—2020 年人口規(guī)模與政府調(diào)控能力是影響土地城鎮(zhèn)化的主要驅(qū)動因子，城鎮(zhèn)人口數(shù)量增長對城鎮(zhèn)居住、商服、基礎(chǔ)設(shè)施等用地需求的增加帶動了城鎮(zhèn)用地持續(xù)擴(kuò)張；經(jīng)濟(jì)發(fā)展水平與產(chǎn)業(yè)結(jié)構(gòu)對土地城鎮(zhèn)化的影響程度不斷增強(qiáng)，消費(fèi)是拉動經(jīng)濟(jì)增長的主要動力；地均固定資產(chǎn)投資額對土地城鎮(zhèn)化的影響強(qiáng)度逐漸下降，表明山東省土地城鎮(zhèn)化對投資的依賴程度有所降低。

交互探測結(jié)果：山東省土地城鎮(zhèn)化時空分異的驅(qū)動因子交互探測結(jié)果如表3 顯示，驅(qū)動因子交互作用呈現(xiàn)雙因子增強(qiáng)和非線性增強(qiáng)效應(yīng)。2000 年，驅(qū)動因子交互表現(xiàn)為雙因子增強(qiáng)關(guān)系的占48.48%。城鎮(zhèn)人口密度∩人均社會消費(fèi)品零售額交互影響最強(qiáng)，q 值為0.865；人均GDP∩地均固定資產(chǎn)投資、人均社會消費(fèi)品零售額∩地均地方財(cái)政收入交互影響較強(qiáng)，q值分別為0.854、0.850。此外，人均GDP 與產(chǎn)業(yè)結(jié)構(gòu)、政府調(diào)控能力、工業(yè)發(fā)展水平方面因子的交互影響程度均達(dá)到0.8 以上，表明雙因子交互作用對土地城鎮(zhèn)化時空分異的影響程度增強(qiáng)。2010年，驅(qū)動因子交互作用呈雙因子增強(qiáng)關(guān)系的占80.30%。城鎮(zhèn)人口密度∩地均地方財(cái)政支出交互影響最強(qiáng)，q 值為0.710；城鎮(zhèn)人口密度∩地均外商投資、距最近交通道路距離∩城鎮(zhèn)人口密度交互影響較強(qiáng)，q 值分別為0.672、0.638，表明城鎮(zhèn)人口規(guī)模與外商投資水平、交通區(qū)位條件的協(xié)同作用有助于促進(jìn)土地城鎮(zhèn)化水平的提升。2020 年，驅(qū)動因子交互作用呈雙因子增強(qiáng)關(guān)系的占74.24%。城鎮(zhèn)人口密度∩地均出口額、城鎮(zhèn)人口密度∩人均社會消費(fèi)品零售額的交互影響最強(qiáng)，q 值分別為0.824、0.814，表明這一時期城鎮(zhèn)人口規(guī)模與外向型經(jīng)濟(jì)發(fā)展、居民消費(fèi)水平的協(xié)同作用是區(qū)域經(jīng)濟(jì)發(fā)展的重要驅(qū)動力，推動城鎮(zhèn)建設(shè)用地?cái)U(kuò)張。綜上，2000—2020年驅(qū)動因子兩兩交互對土地城鎮(zhèn)化時空分異呈現(xiàn)雙因子增強(qiáng)和非線性增強(qiáng)效應(yīng)。

表3 驅(qū)動因子交互作用探測結(jié)果Table 3 Interaction detection results of natural factors

4 結(jié)論與討論

本文基于遙感監(jiān)測、社會統(tǒng)計(jì)等多源數(shù)據(jù)，運(yùn)用空間自相關(guān)分析和地理探測器模型分析了2000—2020年山東省土地城鎮(zhèn)化時空格局演變及驅(qū)動因素，主要結(jié)論如下：①2000—2020 年，山東省土地城鎮(zhèn)化總體處于中等和較低水平，高水平和低水平等級的縣區(qū)占比較低，較高水平和中等水平、較低水平等級的縣區(qū)占比較高。期間：山東省縣域土地城鎮(zhèn)化水平差異擴(kuò)大，極化顯象顯著。其中：高水平縣區(qū)數(shù)量比例增大，呈現(xiàn)向西北部、南部市區(qū)擴(kuò)展的態(tài)勢，集中分布于青島、濰坊、淄博、濟(jì)南、臨沂等地市的市區(qū)；較高水平由西南向東北擴(kuò)展，趨于市區(qū)周邊和交通沿線地區(qū)分布，集中分布于聊城、菏澤、濟(jì)寧境內(nèi)和濟(jì)南、淄博、東營交界地區(qū)；較低及以下水平等級范圍大幅縮減，相關(guān)縣區(qū)趨于沂蒙山區(qū)、半島丘陵地區(qū)和黃河三角洲沿海地區(qū)集中分布。②2000—2020年，山東省縣域土地城鎮(zhèn)化水平空間集聚特征顯著。局部空間格局特征以低—低集聚為主，以高—低集聚和低—高集聚為輔，其演變特征表現(xiàn)為高—高集聚區(qū)不斷擴(kuò)展，并由青島市區(qū)轉(zhuǎn)至濟(jì)南市區(qū)；低—低集聚區(qū)持續(xù)收縮，趨于沂蒙山區(qū)和半島丘陵地區(qū)集中分布，涵蓋平邑縣、蒙陰縣、沂水縣、沂源縣、萊西市、海陽市、乳山市、棲霞市等縣市。③山東省土地城鎮(zhèn)化水平受地理位置、人口規(guī)模、經(jīng)濟(jì)發(fā)展水平、政府調(diào)控能力、工業(yè)發(fā)展水平、開放水平等因素的共同影響，城鎮(zhèn)人口密度，地均地方財(cái)政收入，地均固定資產(chǎn)投資額，地均地方財(cái)政支出，第二、三產(chǎn)業(yè)比重，距最近地市距離，地均出口額，地均外商投資，人均社會消費(fèi)品零售額，地均工業(yè)產(chǎn)值和人均GDP是影響土地城鎮(zhèn)化水平的重要驅(qū)動因子。地理區(qū)位、人口規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)、政府調(diào)控能力、工業(yè)發(fā)展水平、開放水平等因素對土地城鎮(zhèn)化空間分異的影響增強(qiáng)，人均GDP 和地均固定資產(chǎn)投資的影響減弱。各驅(qū)動因子的交互作用對土地城鎮(zhèn)化時空分異呈現(xiàn)雙因子增強(qiáng)和非線性增強(qiáng)效應(yīng)。

本文基于縣域尺度，運(yùn)用GIS 空間分析和地理探測器模型探討了2000—2020 年山東省土地城鎮(zhèn)化時空格局演變特征及影響因素，有助于厘清快速城鎮(zhèn)化背景下區(qū)域城鎮(zhèn)土地利用的變化趨勢及動力機(jī)制。運(yùn)用地理探測器模型探究土地城鎮(zhèn)化時空分異的驅(qū)動因素，可以有效探究土地城鎮(zhèn)化驅(qū)動因素之間的相互作用機(jī)理。從縣域城鎮(zhèn)面積和人口規(guī)模的影響來看，山東省土地城鎮(zhèn)化水平在未來一定時間內(nèi)仍將持續(xù)提高。隨著生態(tài)文明建設(shè)與糧食安全保障戰(zhàn)略的實(shí)施，資源環(huán)境約束日益增大，山東省土地城鎮(zhèn)化建設(shè)應(yīng)合理引導(dǎo)人口集聚，嚴(yán)格控制土地增量，盤活土地存量，實(shí)現(xiàn)集約化和緊湊化發(fā)展道路。經(jīng)濟(jì)基礎(chǔ)較好的縣區(qū)應(yīng)加強(qiáng)存量土地探潛，推進(jìn)城鎮(zhèn)化建設(shè)由增量擴(kuò)張向存量發(fā)展轉(zhuǎn)變；經(jīng)濟(jì)基礎(chǔ)較薄弱的縣區(qū)應(yīng)充分利用地方資源稟賦，培育特色優(yōu)勢產(chǎn)業(yè)，推進(jìn)產(chǎn)業(yè)結(jié)構(gòu)升級。本研究尚存著不足，有待進(jìn)一步完善。鑒于數(shù)據(jù)的可獲得性和操作性，進(jìn)行土地城鎮(zhèn)化驅(qū)動因素分析時僅選取地均地方財(cái)政收入與支出表征財(cái)政發(fā)展變化對土地城鎮(zhèn)化水平的影響也是值得深入探討的問題。