楊喜 盧新海

摘要 從空間效應(yīng)視角下分析城市土地城鎮(zhèn)化驅(qū)動(dòng)因素可以為優(yōu)化城市土地資源可持續(xù)開(kāi)發(fā)利用提供新參考?；?003—2017年中國(guó)283座城市面板數(shù)據(jù)，構(gòu)建經(jīng)濟(jì)地理空間權(quán)重矩陣，借助空間面板杜賓模型及其偏分分解方法分析中國(guó)城市土地城鎮(zhèn)化驅(qū)動(dòng)因素。研究表明：①研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化水平呈持續(xù)上升特征，城際差異呈持續(xù)擴(kuò)大特征，并伴隨著極化現(xiàn)象。②中國(guó)城市土地城鎮(zhèn)化在地理空間上并非隨機(jī)分布，而是存在顯著正向全局空間自相關(guān)性，局域空間格局特征表現(xiàn)為高值集聚、低值塌陷、低值集聚和高值凸起四種類型，數(shù)量上以HH和LL為主，以LH和HL為輔，并且存在“低者恒低、高者恒高”的空間格局特征。③中國(guó)城市土地城鎮(zhèn)化存在顯著正向空間溢出效應(yīng)，本地城市土地城鎮(zhèn)化水平提升會(huì)驅(qū)動(dòng)鄰近城市土地城鎮(zhèn)化水平提升，這打破了傳統(tǒng)計(jì)量研究中樣本相互獨(dú)立的基本假設(shè)，只有將空間效應(yīng)納入城市土地城鎮(zhèn)化的計(jì)量回歸模型之中，才能合理解釋城市土地城鎮(zhèn)化驅(qū)動(dòng)因素。④經(jīng)濟(jì)發(fā)展、第二產(chǎn)業(yè)和土地財(cái)政對(duì)本地城市土地城鎮(zhèn)化具有正向直接效應(yīng);經(jīng)濟(jì)發(fā)展、人口城鎮(zhèn)化和對(duì)外開(kāi)放對(duì)鄰近城市土地城鎮(zhèn)化具有負(fù)向溢出效應(yīng)，土地財(cái)政對(duì)鄰近城市土地城鎮(zhèn)化具有正向溢出效應(yīng)。據(jù)此，管控城市土地城鎮(zhèn)化不可“孤城而為”，需要綜合考慮鄰近城市之間空間溢出效應(yīng)的存在，發(fā)揮城際之間在推進(jìn)土地城鎮(zhèn)化理性發(fā)展中的合作聯(lián)動(dòng)治理作用。

關(guān)鍵詞 空間效應(yīng);土地城鎮(zhèn)化;經(jīng)濟(jì)地理空間權(quán)重;空間面板杜賓模型

土地城鎮(zhèn)化是城市生產(chǎn)、生活與生態(tài)關(guān)系在城市空間上的綜合映射，是推動(dòng)城市邁向現(xiàn)代化發(fā)展的空間載體。然而，城市土地城鎮(zhèn)化超前發(fā)展以及無(wú)序蔓延備受詬病，卻又在各地屢禁不止，由此導(dǎo)致城市生產(chǎn)空間、生活空間和生態(tài)空間之間的矛盾日益加劇，對(duì)城市可持續(xù)發(fā)展構(gòu)成了威脅。根據(jù)《中國(guó)城市建設(shè)統(tǒng)計(jì)年鑒》統(tǒng)計(jì)發(fā)現(xiàn)，2003—2017年中國(guó)城市建設(shè)用地面積由從28 308 km2增長(zhǎng)至56 225 km2，增長(zhǎng)了98.62%。2019年11月中辦、國(guó)辦印發(fā)《關(guān)于在國(guó)土空間規(guī)劃中統(tǒng)籌劃定落實(shí)三條控制線的指導(dǎo)意見(jiàn)》，明確提出將城鎮(zhèn)開(kāi)發(fā)邊界、永久基本農(nóng)田和生態(tài)保護(hù)紅線三條控制線作為調(diào)整經(jīng)濟(jì)結(jié)構(gòu)、規(guī)劃產(chǎn)業(yè)發(fā)展、推進(jìn)城鎮(zhèn)化不可逾越的紅線[1]。在三大國(guó)土空間紅線強(qiáng)力管控下，傳統(tǒng)依靠外延空間無(wú)序擴(kuò)張的城市發(fā)展路徑將會(huì)被遏制。但各城市發(fā)展并非孤立地存在，城市土地城鎮(zhèn)化快速發(fā)展也并非一城之因，人流、物流、信息流和資金流在空間流動(dòng)過(guò)程中將不同的城市串聯(lián)在一起，既往研究基于樣本相互獨(dú)立的假設(shè)與現(xiàn)實(shí)不符。因此，從空間效應(yīng)視角下加強(qiáng)對(duì)中國(guó)城市土地城鎮(zhèn)化驅(qū)動(dòng)因素的研究，對(duì)有效治理城市土地城鎮(zhèn)化超前發(fā)展和促進(jìn)城市土地資源可持續(xù)開(kāi)發(fā)利用具有十分重要的意義。

1 文獻(xiàn)綜述

城鎮(zhèn)化是推動(dòng)人類文明發(fā)展進(jìn)步的重要驅(qū)動(dòng)力，城鎮(zhèn)化的核心是人口城鎮(zhèn)化，城鎮(zhèn)化的載體是土地城鎮(zhèn)化，作為城鎮(zhèn)化一體兩面之一的土地城鎮(zhèn)化受到了研究者廣泛關(guān)注。當(dāng)前學(xué)者圍繞土地城鎮(zhèn)化的研究主要經(jīng)歷了從空間特征分析到與人口城鎮(zhèn)化關(guān)系分析再到驅(qū)動(dòng)因素分析的研究脈絡(luò)。一是土地城鎮(zhèn)化空間特征，如盧斌瑩等[2]借助GIS分析了西安市南郊土地城市化空間特征;楊洋等[3]利用夜間燈光數(shù)據(jù)對(duì)環(huán)渤海地區(qū)土地城鎮(zhèn)空間特征進(jìn)行了分析;崔許鋒[4]分析了云南省土地城鎮(zhèn)化空間異質(zhì)性特征。二是土地城鎮(zhèn)化與人口城鎮(zhèn)關(guān)系，如劉瓊等[5]、于立等[6]分別從不同的時(shí)間階段視角來(lái)研究土地城鎮(zhèn)化與人口城鎮(zhèn)化協(xié)調(diào)關(guān)系;李寶禮等[7]、許芬[8]、李子聯(lián)[9]探究了土地城鎮(zhèn)化快于人口城鎮(zhèn)化的原因;王成新等[10]從結(jié)構(gòu)視角探討了土地城市化與人口城市化異速增長(zhǎng)的原因。三是土地城鎮(zhèn)化驅(qū)動(dòng)因素，如李永樂(lè)等[11]、吳一凡等[12]采用最小二乘回歸計(jì)量模型（ordinary least squares，OLS）分析了房?jī)r(jià)、經(jīng)濟(jì)發(fā)展水平和人口規(guī)模等因素對(duì)土地城鎮(zhèn)化的影響;張立新等[13]采用偏最小二乘回歸計(jì)量模型（partial least squares， PLS）分析了經(jīng)濟(jì)發(fā)展水平、人口規(guī)模和產(chǎn)業(yè)結(jié)構(gòu)升級(jí)等因素對(duì)長(zhǎng)江經(jīng)濟(jì)帶土地城鎮(zhèn)化的影響;劉耀林等[14]采用Logistic計(jì)量模型分析了城鎮(zhèn)化水平、產(chǎn)業(yè)結(jié)構(gòu)和地形條件等因素對(duì)湖北省城鄉(xiāng)建設(shè)用地城鎮(zhèn)化的影響;王曉云等[15]采用了Tobit計(jì)量模型分析了經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)和外資利用規(guī)模等因素對(duì)甘肅省土地城鎮(zhèn)化的影響;盧新海等[16]采用了空間誤差模型（spatial error model，SEM）分析了經(jīng)濟(jì)發(fā)展、人口規(guī)模和產(chǎn)業(yè)結(jié)構(gòu)等因素對(duì)中部地區(qū)土地城鎮(zhèn)化的影響，雖然該研究考慮到了空間效應(yīng)，但是該研究尚未將被解釋變量和解釋變量的空間效應(yīng)共同納入模型分析之中，并且該研究只分析了截面數(shù)據(jù)上的誤差項(xiàng)空間效應(yīng)。通過(guò)文獻(xiàn)梳理可知，當(dāng)前學(xué)界對(duì)土地城鎮(zhèn)化的研究已經(jīng)取得了較為豐富的成果，但是還存在進(jìn)一步推進(jìn)研究的方面：一是，在研究尺度上，現(xiàn)有研究較多關(guān)注單個(gè)城市尺度、省級(jí)尺度以及區(qū)域級(jí)尺度，而較少關(guān)注全國(guó)地級(jí)市尺度土地城鎮(zhèn)化的研究，造成對(duì)全國(guó)地級(jí)市層面土地城鎮(zhèn)化空間格局特征和驅(qū)動(dòng)因素把握不足;二是，在驅(qū)動(dòng)因素方法選擇上，現(xiàn)有文獻(xiàn)多采用未考慮空間效應(yīng)的基于樣本相互獨(dú)立假設(shè)的研究方法，較少使用空間面板計(jì)量模型來(lái)揭示土地城鎮(zhèn)化驅(qū)動(dòng)因素，忽略了鄰近城市在土地資源管理、產(chǎn)業(yè)布局和經(jīng)濟(jì)發(fā)展政策與策略所存在的空間互動(dòng)效應(yīng)，導(dǎo)致對(duì)土地城鎮(zhèn)化驅(qū)動(dòng)因素模型估計(jì)上存在偏誤;三是，在研究重點(diǎn)上，現(xiàn)有文獻(xiàn)鮮有關(guān)注土地城鎮(zhèn)化空間溢出效應(yīng)的研究，也鮮有文獻(xiàn)采用空間面板計(jì)量模型中的偏微分分解方法來(lái)分析土地城鎮(zhèn)化各驅(qū)動(dòng)因素的直接效應(yīng)與溢出效應(yīng)。

基于此，本文在借鑒現(xiàn)有文獻(xiàn)研究基礎(chǔ)上，采集中國(guó)283座城市2003—2017年面板數(shù)據(jù)，利用空間自相關(guān)來(lái)識(shí)別中國(guó)城市土地城鎮(zhèn)化全局和局域空間格局特征，借助空間面板杜賓模型對(duì)中國(guó)城市土地城鎮(zhèn)化空間溢出效應(yīng)進(jìn)行參數(shù)化估計(jì)，并利用偏微分分解方法對(duì)城市土地城鎮(zhèn)化各驅(qū)動(dòng)因素的空間效應(yīng)進(jìn)行分解分析，以期為城市以及城市之間空間治理和土地城鎮(zhèn)化理性發(fā)展提供參考借鑒。

2 研究方法

2.1 空間自相關(guān)

（1）全局空間自相關(guān)。全局空間自相關(guān)是從全域視角對(duì)研究對(duì)象空間相關(guān)性進(jìn)行度量的指標(biāo)，本文采用全局Morans I指數(shù)來(lái)度量城市土地城鎮(zhèn)化全局空間格局特征，其取值范圍為[-1，1]，若大于0表示空間正相關(guān)，若小于0表示空間負(fù)相關(guān)，若等于0表示屬性數(shù)據(jù)在空間上隨機(jī)分布，不具有空間自相關(guān)。對(duì)于空間自相關(guān)顯著性水平可通過(guò)標(biāo)準(zhǔn)化統(tǒng)計(jì)量閾值Z來(lái)進(jìn)行檢驗(yàn)，當(dāng)Z>1.96或P<0.05，表明通過(guò)顯著性水平檢驗(yàn)，意味著地理屬性觀測(cè)值存在顯著的空間自相關(guān)，函數(shù)公式參考文獻(xiàn)[17]。

式中，n為樣本城市數(shù)量;xi和xj分別表示城市i和城市j的觀測(cè)值;為觀測(cè)值平均數(shù);S2為樣本方差;wij為城市i和城市j在空間上的鄰近關(guān)系。

（2）局部空間自相關(guān)。局部空間自相關(guān)是從局域視角對(duì)研究對(duì)象空間相關(guān)性進(jìn)行度量的指標(biāo)，測(cè)度的是本地區(qū)域觀測(cè)值與周圍鄰近區(qū)域觀測(cè)值之間的相似程度，可以用來(lái)識(shí)別局域空間格局集聚和離散特征。本文采用局域Morans Ii指數(shù)來(lái)識(shí)別城市土地城鎮(zhèn)化局域空間格局特征，具體可分為四種集聚類型：HH（High-High）、LH（Low-High）、LL（Low-Low）和HL（High-Low），函數(shù)公式參考文獻(xiàn)[17]，式中變量含義與式（1）相同。

2.2 空間計(jì)量模型

（1）空間面板計(jì)量模型設(shè)定?？臻g計(jì)量模型將長(zhǎng)期被傳統(tǒng)計(jì)量模型所忽略的空間效應(yīng)納入模型中，可用來(lái)識(shí)別不同地理空間單元之間的空間交互關(guān)系?，F(xiàn)階段空間面板計(jì)量模型使用最多和比較成熟的有三種設(shè)定形式。

①空間面板滯后模型（spatial pcnel lag model，SPLM）描述的是空間實(shí)質(zhì)相關(guān)，空間效應(yīng)主要以被解釋變量滯后形式存在于計(jì)量模型設(shè)定中。SPLM模型設(shè)定為：

②空間面板誤差模型（spatial panel error model，SPEM）描述的是空間擾動(dòng)相關(guān)，空間效應(yīng)主要以誤差滯后形式存在于計(jì)量模型設(shè)定中。SPEM模型設(shè)定為：

③空間面板杜賓模型（spatial panel Dubin model，SPDM）同時(shí)將被解釋變量和解釋變量的空間效應(yīng)納入計(jì)量模型設(shè)定中。SPDM模型設(shè)定為：

式中，Yit為被解釋變量的觀測(cè)值，在本文為城市土地城鎮(zhèn)化;Xit為解釋變量的觀測(cè)值，在本文為影響被解釋變量變動(dòng)的各影響因素;ρ為被解釋變量的空間回歸系數(shù);λ為誤差項(xiàng)的空間回歸系數(shù);θ為解釋變量的空間回歸系數(shù);β為解釋變量的回歸系數(shù);μi、υt分別表示空間固定效應(yīng)和時(shí)間固定效應(yīng);εit為服從獨(dú)立同分布的隨機(jī)誤差項(xiàng)，表示未被納入計(jì)量模型中的其他因素;Wij為空間權(quán)重矩陣。當(dāng)θ=0、ρ≠0時(shí)，SPDM模型就會(huì)退化為SPLM模型;當(dāng)θ+ρβ=0時(shí)，SPDM模型就會(huì)退化為SPEM模型。

（2）空間權(quán)重矩陣構(gòu)建。本文借鑒譚術(shù)魁等[18]空間權(quán)重構(gòu)建思路，將城市之間的地理特征和經(jīng)濟(jì)特征同時(shí)考慮在內(nèi)來(lái)構(gòu)建綜合的經(jīng)濟(jì)地理空間權(quán)重矩陣。構(gòu)建此空間權(quán)重矩陣的主要依據(jù)是，一方面，根據(jù)Tobler[19]地理學(xué)第一定律，地區(qū)之間相互影響的關(guān)系會(huì)隨地理距離增大而減小，由此構(gòu)建基于地理距離衰減的地理距離矩陣，該權(quán)重可以很好地體現(xiàn)即使區(qū)域之間空間上不鄰但也會(huì)具有要素流動(dòng)的客觀現(xiàn)實(shí)情況，同時(shí)該權(quán)重也可以把鄰近關(guān)系隨距離衰減的最本質(zhì)特征反映出來(lái)。另一方面，城市之間互動(dòng)影響的能力并不是完全對(duì)等的，經(jīng)濟(jì)實(shí)力強(qiáng)的城市往往對(duì)其他城市影響更大、影響范圍更廣，由此根據(jù)城市之間經(jīng)濟(jì)實(shí)力上的差異，賦予各城市不同的經(jīng)濟(jì)空間權(quán)重。綜合以上考慮構(gòu)建經(jīng)濟(jì)地理空間權(quán)重矩陣：

式中，Wij為經(jīng)濟(jì)地理空間權(quán)重矩陣;wij為地理距離矩陣;diag（·）為對(duì)角矩陣;d 為城市之間的地理距離;i為第i個(gè)城市在研究期內(nèi)從時(shí)間t0到t1時(shí)段地區(qū)經(jīng)濟(jì)產(chǎn)值的平均值;為研究期內(nèi)從時(shí)間t0到t1時(shí)段所有城市地區(qū)經(jīng)濟(jì)產(chǎn)值的平均值;n為城市總數(shù)量。

（3）空間效應(yīng)分解。由于變量之間存在空間上的相關(guān)性，故此參考學(xué)術(shù)界常規(guī)做法采取極大似然估計(jì)（maximum likelihood estimation，MLE）方法對(duì)空間計(jì)量模型進(jìn)行估計(jì)。同時(shí)根據(jù)Lesage等[20]的偏微分方法將解釋變量對(duì)被解釋變量的空間效應(yīng)分解為直接效應(yīng)和間接效應(yīng)，將SPDM模型移項(xiàng)整理成一般形式：

式中，分解矩陣中主對(duì)角線上的元素用來(lái)捕捉直接效應(yīng)，分解矩陣中非對(duì)角線上的元素用來(lái)捕捉間接效應(yīng)。由偏微分分解公式可以發(fā)現(xiàn)，某一特定單元解釋變量的變動(dòng)，不僅會(huì)對(duì)該單元自身被解釋變量的變動(dòng)產(chǎn)生影響，而且同時(shí)還會(huì)對(duì)其他單元被解釋變量的變動(dòng)產(chǎn)生影響，前者稱為直接效應(yīng)，后者稱為間接效應(yīng)（空間溢出效應(yīng)）。如果ρ=0并且θk=0，則不存在空間溢出效應(yīng)。

3 變量選取與數(shù)據(jù)來(lái)源

3.1 變量選取

土地城鎮(zhèn)化最為直接的表現(xiàn)就是城市用地對(duì)非城市用地的占用和替代，是城市建設(shè)用地?cái)U(kuò)張和蔓延的過(guò)程，本文借鑒吳一凡等[12]、盧新海等[16] 和孫焱林等[21]研究方法，采用城市建設(shè)用地面積與城市總面積之比來(lái)表征城市土地城鎮(zhèn)化。

土地城鎮(zhèn)化驅(qū)動(dòng)因素的選擇，綜合考慮地級(jí)市城市特性以及數(shù)據(jù)的可獲得性與一致性，并參考借鑒前人研究的基礎(chǔ)，從經(jīng)濟(jì)發(fā)展、人口城鎮(zhèn)化、產(chǎn)業(yè)結(jié)構(gòu)、土地財(cái)政和對(duì)外開(kāi)放來(lái)選擇解釋變量。①城市經(jīng)濟(jì)發(fā)展（Pgdp）。城市經(jīng)濟(jì)發(fā)展是城市實(shí)力的綜合體現(xiàn)，是推動(dòng)城市土地城鎮(zhèn)化發(fā)展的基礎(chǔ)動(dòng)力，同時(shí)也是城市之間相互影響能力大小的主要因素，借鑒張立新等[13]方法采用人均實(shí)際GDP表示。②城市人口城鎮(zhèn)化（Pur）。城市人口的增長(zhǎng)加大城市空間承載的壓力，會(huì)迫使城市向外拓展發(fā)展空間，借鑒楊孟禹等[22]方法采用市轄區(qū)人口與城市總?cè)丝跀?shù)之比表示。③城市產(chǎn)業(yè)結(jié)構(gòu)（Sei）。第二產(chǎn)業(yè)具有土地資源消耗大的特征，城市第二產(chǎn)業(yè)占比上升也會(huì)推動(dòng)土地城鎮(zhèn)化的擴(kuò)張，借鑒李成宇等[23]方法采用第二產(chǎn)業(yè)產(chǎn)值與GDP之比表示。④城市土地財(cái)政 （Lfd）。土地財(cái)政可以反映出地方政府“以地生財(cái)”的能力，是推動(dòng)土地城鎮(zhèn)化的重要資金來(lái)源，土地財(cái)政也會(huì)存在示范效應(yīng)，城市之間會(huì)競(jìng)相模仿，借鑒邵朝對(duì)等[24]方法采用人均實(shí)際土地出讓收入表示。⑤城市對(duì)外開(kāi)放（Fdi）。學(xué)界通常用外商投資來(lái)衡量城市對(duì)外開(kāi)放水平，外商投資的引進(jìn)可能存在兩種效應(yīng)，一是其帶來(lái)的擠出效應(yīng)，外商投資憑借其技術(shù)和管理優(yōu)勢(shì)對(duì)本地土地資源消耗大的產(chǎn)業(yè)產(chǎn)生擠出效應(yīng)，緩解土地城鎮(zhèn)化快速蔓延;二是城市之間競(jìng)爭(zhēng)效應(yīng)，城市在吸引外資上存在以“以地引資”，加大土地資源消耗推動(dòng)了土地城鎮(zhèn)化。借鑒周琳等[25]方法采用 FDI與GDP之比表示。為了增強(qiáng)數(shù)據(jù)的可比性，利用GDP指數(shù)對(duì)GDP數(shù)據(jù)做平減處理，其他價(jià)格數(shù)據(jù)利用CPI指數(shù)進(jìn)行銷賬處理，利用匯率將外商投資數(shù)據(jù)調(diào)整為人民幣計(jì)價(jià)。

3.2 數(shù)據(jù)來(lái)源

由于部分城市數(shù)據(jù)缺失較多，為了保持?jǐn)?shù)據(jù)前后的一致性和連貫性，最終以中國(guó)大陸283座地級(jí)及以上城市為研究對(duì)象，變量數(shù)據(jù)主要來(lái)源于《中國(guó)城市統(tǒng)計(jì)年鑒（2003—2017）》和《中國(guó)建設(shè)統(tǒng)計(jì)年鑒（2003—2017）》，部分缺失數(shù)據(jù)通過(guò)查詢各省市統(tǒng)計(jì)年鑒和統(tǒng)計(jì)公報(bào)補(bǔ)充。由于城市GDP指數(shù)數(shù)據(jù)和CPI指數(shù)數(shù)據(jù)缺失較多，采用城市所在省級(jí)數(shù)據(jù)來(lái)替代。匯率數(shù)據(jù)和歷年各省GDP指數(shù)和各省CPI指數(shù)數(shù)據(jù)來(lái)源于國(guó)家統(tǒng)計(jì)局。各城市地理空間位置坐標(biāo)信息借助ArcGIS軟件提取，最終整理獲得2003—2017年中國(guó)283座城市面板數(shù)據(jù)。

4 實(shí)證結(jié)果與分析

4.1 時(shí)序演進(jìn)特征分析

圖1為借助Eviews10軟件繪制的中國(guó)城市土地城鎮(zhèn)化核密度曲線圖，借此來(lái)刻畫中國(guó)城市土地城鎮(zhèn)化時(shí)序演進(jìn)特征。①?gòu)暮嗣芏惹€位置來(lái)看，2003—2017年曲線位置整體呈現(xiàn)出向右逐漸遷移態(tài)勢(shì)，表明在研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化水平呈現(xiàn)出持續(xù)上升演進(jìn)特征。②從曲線主峰波峰高度來(lái)看，2003—2017年呈現(xiàn)出逐漸下降態(tài)勢(shì)，由尖峰形態(tài)逐漸向?qū)挿逍螒B(tài)演變，表明在研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化城際差異呈現(xiàn)持續(xù)擴(kuò)大化演進(jìn)特征。③從曲線波峰數(shù)量來(lái)看，2003—2017年呈現(xiàn)出以一個(gè)主峰為主和右側(cè)隆起次峰為輔的形態(tài)，表明在研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化呈現(xiàn)出了極化特征。④從曲線拖尾來(lái)看，2003—2017年右側(cè)拖尾出現(xiàn)了延長(zhǎng)、抬厚趨勢(shì)，表明在研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化處于高值區(qū)城市數(shù)量比例呈現(xiàn)出擴(kuò)大演進(jìn)特征。

4.2 空間格局特征分析

（1）全局空間格局特征。由表1可知，2003—2017年中國(guó)城市土地城鎮(zhèn)化全局Morans I指數(shù)處于0.406～0.427 之間，并且Z-Value >1.96，P-Value<0.01，各年都通過(guò)1%顯著性水平檢驗(yàn)，表明中國(guó)城市土地城鎮(zhèn)化在地理空間上并非隨機(jī)分布，而是存在顯著正向全局空間自相關(guān)性。從總體上看，全局Morans I指數(shù)呈現(xiàn)上升趨勢(shì)，具體而言，2003—2011年呈上升趨勢(shì)，2011—2015年呈下降趨勢(shì)，2015—2017年呈上升趨勢(shì)，表明研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化全局空間關(guān)聯(lián)呈現(xiàn)由弱相關(guān)向強(qiáng)相關(guān)演變特征，但年際間存在一定的波動(dòng)性，城市土地城鎮(zhèn)化全局空間格局關(guān)系特征尚未形成穩(wěn)態(tài)。通過(guò)全局空間自相關(guān)，可以發(fā)現(xiàn)傳統(tǒng)研究中基于樣本相互獨(dú)立的基本假設(shè)與現(xiàn)實(shí)不符，因此需要將空間效應(yīng)納入城市土地城鎮(zhèn)化的計(jì)量回歸模型之中。

（2）局域空間格局特征。由圖2可知，2003—2017年中國(guó)城市土地城鎮(zhèn)化局域空間格局特征具體可歸為四類：第一類高值集聚型 （HH），該類型表示本地城市土地城鎮(zhèn)化水平高，其周圍鄰近城市土地城鎮(zhèn)化水平高，呈現(xiàn)為“中心高，四周高”的高水平均值空間關(guān)聯(lián)集聚狀態(tài)。第二類低值塌陷型 （LH），該類型表示本地城市土地城鎮(zhèn)化水平低，但其周圍鄰近城市土地城鎮(zhèn)化水平高，呈現(xiàn)為“中心低，四周高”的空間非均衡關(guān)聯(lián)集聚狀態(tài)。第三類低值集聚型（LL），該類型表示本地城市土地城鎮(zhèn)化水平低，其周圍鄰近城市土地城鎮(zhèn)化水平也低，呈現(xiàn)為“中心低，四周低”的低水平均值空間關(guān)聯(lián)集聚狀態(tài)。第四類高值凸起型（HL），該類型表示本地城市土地城鎮(zhèn)化水平高，但其周圍鄰近城市土地城鎮(zhèn)化水平低呈現(xiàn)為“中心高，四周低”的空間非均衡關(guān)聯(lián)集聚狀態(tài)。從數(shù)量上來(lái)看，中國(guó)城市土地城鎮(zhèn)化水平在研究期內(nèi)所呈現(xiàn)的是以高值集聚和低值集聚為主，以低值塌陷和高值凸起為輔。城市土地城鎮(zhèn)化水平四種集聚類型走勢(shì)上具有較強(qiáng)的路徑依賴性，但各類型之間大致呈小幅波動(dòng)的平行走勢(shì)，類型之間很難發(fā)生突破性的大躍遷，即當(dāng)鄰近城市土地城鎮(zhèn)化處于高（低）水平值時(shí)，本地也更容易成為高（低）水平值城市，或者本地也更難向低（高）水平值城市轉(zhuǎn)移，由此城市土地城鎮(zhèn)化容易形成“低者恒低、高者恒高”的局域空間格局特征。主要原因是城市土地城鎮(zhèn)化地理空間上存在鄰近同伴效應(yīng)和空間溢出效應(yīng)，鄰近城市在經(jīng)濟(jì)發(fā)展中土地資源管理和產(chǎn)業(yè)布局上所采取的政策措施和策略行為上具有互動(dòng)模仿、學(xué)習(xí)和競(jìng)爭(zhēng)效應(yīng)，經(jīng)濟(jì)結(jié)構(gòu)上存在趨于或互補(bǔ)現(xiàn)象。

4.3 空間計(jì)量模型檢驗(yàn)與選擇

空間計(jì)量包括空間滯后（spatial lag model，SLM）空間誤差（spatial error model，SEM）以及空間杜賓（spatial Dubin model，SDM）有三大基礎(chǔ)模型，而在具體模型的選擇上還需要進(jìn)行一系列的檢驗(yàn)。由表2可知，在LM檢驗(yàn)中 LM-lag 和LM-error檢驗(yàn)都通過(guò)了1%顯著性水平檢驗(yàn)，在Robust LM檢驗(yàn)中LM-lag未通過(guò)10%顯著性水平檢驗(yàn)，LM-error檢驗(yàn)通過(guò)了1%顯著性水平檢驗(yàn)，表明空間滯后模型和空間誤差模型可能同時(shí)存在，那么就傾向于選擇空間杜賓模型[20]。此時(shí)還需要對(duì)空間杜賓模型進(jìn)行檢驗(yàn)，檢驗(yàn)其是否可以退化為空間滯后模型或空間誤差模型，其中Wald檢驗(yàn)和LR檢驗(yàn)都在1%顯著性水平上拒絕了空間杜賓模型可以退化為空間滯后模型和空間誤差模型的原假設(shè)，因此應(yīng)當(dāng)選擇空間杜賓模型[26]。同時(shí)Hausman 檢驗(yàn)在1%顯著性水平上拒絕隨機(jī)效應(yīng)的原假設(shè)，所以最終采用固定效應(yīng)空間面板杜賓模型。

根據(jù)表3空間面板杜賓模型估計(jì)結(jié)果可以發(fā)現(xiàn)，被解釋變量空間回歸系數(shù)ρ在三種固定效應(yīng)估計(jì)下系數(shù)都為正，并且均通過(guò)了1%顯著性水平檢驗(yàn)，結(jié)果較為穩(wěn)健，由此可以發(fā)現(xiàn)中國(guó)城市土地城鎮(zhèn)化存在顯著的正向空間溢出效應(yīng)，本地城市土地城鎮(zhèn)化對(duì)其他鄰近城市土地城鎮(zhèn)化具有顯著的正向影響，這也進(jìn)一步驗(yàn)證了前面的統(tǒng)計(jì)分析所得結(jié)果和將空間效應(yīng)納入計(jì)量模型中的合理性和必要性。在三種固定效應(yīng)估計(jì)中，時(shí)空固定效應(yīng)模型擬合優(yōu)度（R2）和對(duì)數(shù)似然值（log-likelihood）最大，并且sigma2也是最小。綜合上述檢驗(yàn)結(jié)果，本文最終選擇更為一般的時(shí)空雙固定效應(yīng)空間面板杜賓模型作為最終的分析模型。但是空間杜賓模型所得出的結(jié)果，解釋變量系數(shù)并不直接表示邊際效應(yīng)，也不表示直接效應(yīng)和空間溢出效應(yīng)[27]，所以并不能采取傳統(tǒng)的非空間模型解釋系數(shù)的方法對(duì)其加以解釋[20]，因而其并不是本文研究的主要分析對(duì)象。

4.4 空間效應(yīng)分解分析

表4為在時(shí)空雙固定效應(yīng)下空間面板杜賓基礎(chǔ)上，采取偏微分方法對(duì)城市土地城鎮(zhèn)化各驅(qū)動(dòng)因素進(jìn)行的空間效應(yīng)分解結(jié)果。

（1）經(jīng)濟(jì)發(fā)展水平直接效應(yīng)系數(shù)為0.435 8，溢出效應(yīng)系數(shù)為-0.225 7，并且二者都通過(guò)了1%顯著性水平檢驗(yàn)。表明本地城市經(jīng)濟(jì)發(fā)展水平提升會(huì)推動(dòng)本地城市土地城鎮(zhèn)化發(fā)展，這與高金龍等[28]研究結(jié)論一致，但是會(huì)通過(guò)負(fù)向溢出效應(yīng)抑制鄰近城市土地城鎮(zhèn)化發(fā)展。主要原因是經(jīng)濟(jì)發(fā)展水平高的城市，經(jīng)濟(jì)增長(zhǎng)中具有自我強(qiáng)化功能，投資回報(bào)往往也較高，容易成為資本、勞動(dòng)要素流入的目的地，需要有更多土地資源作為城市發(fā)展空間承載體，提高城市在發(fā)展過(guò)程中對(duì)土地資源的需求，推動(dòng)了城市空間生長(zhǎng)。同時(shí)本地城市經(jīng)濟(jì)發(fā)展水平提升，對(duì)鄰近經(jīng)濟(jì)發(fā)展較弱城市資本、人力等流動(dòng)性生產(chǎn)要素產(chǎn)生虹吸效應(yīng)，而鄰近城市生產(chǎn)要素的流出會(huì)降低對(duì)土地資源的需求，抑制了鄰近城市空間的進(jìn)一步擴(kuò)展。

（2）人口城鎮(zhèn)化水平直接效應(yīng)系數(shù)為0.037 9，未通過(guò)10%顯著性水平檢驗(yàn)，溢出效應(yīng)系數(shù)為-0.267 2，通過(guò)了1%顯著性水平檢驗(yàn)。表明本地城市人口城鎮(zhèn)化水平提升對(duì)本地城市土地城鎮(zhèn)化尚未形成顯著的促進(jìn)作用，但會(huì)通過(guò)負(fù)向溢出效應(yīng)抑制鄰近城市土地城鎮(zhèn)化的發(fā)展。主要原因是當(dāng)前城市土地城鎮(zhèn)化并非由人口城鎮(zhèn)化帶來(lái)的自然引致需求，存在人口城鎮(zhèn)化滯后于土地城鎮(zhèn)化[9]和“化地不化人”[29]的城鎮(zhèn)化現(xiàn)象，以及傳統(tǒng)“重體量不重質(zhì)量”的發(fā)展模式造成了城鎮(zhèn)化演化的無(wú)序[30]，由此造成人口城鎮(zhèn)化對(duì)本地土地城鎮(zhèn)化促進(jìn)作用并未通過(guò)統(tǒng)計(jì)學(xué)意義上的檢驗(yàn)。同時(shí)在人口自由流動(dòng)和人口增長(zhǎng)率下降情況下，本地人口城鎮(zhèn)化水平的提升，鄰近城市人口就會(huì)相應(yīng)的流出，抑制了鄰近城市由人口效應(yīng)帶來(lái)驅(qū)動(dòng)土地城鎮(zhèn)化的基礎(chǔ)動(dòng)力。

（3）產(chǎn)業(yè)結(jié)構(gòu)水平直接效應(yīng)系數(shù)為0.129 9，通過(guò)了1%顯著性水平檢驗(yàn)，溢出效應(yīng)系數(shù)為0.036 6，未通過(guò)10%顯著性水平檢驗(yàn)。表明本地城市第二產(chǎn)業(yè)結(jié)構(gòu)水平提升會(huì)推動(dòng)本地城市土地城鎮(zhèn)化發(fā)展，這與呂立剛等[31]研究結(jié)論一致，對(duì)鄰近城市土地城鎮(zhèn)化尚未形成顯著的推動(dòng)作用。主要原因是第二產(chǎn)業(yè)比重增大，意味著第三產(chǎn)業(yè)比重縮小，而第二產(chǎn)業(yè)往往是土地資源需求和土地資源消耗大的產(chǎn)業(yè)，其比重提升會(huì)進(jìn)一步推動(dòng)城市土地城鎮(zhèn)化發(fā)展。同時(shí)第二產(chǎn)業(yè)由于投資的地域性和固定性特征表現(xiàn)為較強(qiáng)的本地效應(yīng)，造成其對(duì)鄰近城市土地城鎮(zhèn)化推動(dòng)作用并未通過(guò)統(tǒng)計(jì)學(xué)意義上的檢驗(yàn)。

（4）土地財(cái)政水平直接效應(yīng)系數(shù)為0.016 0，溢出效應(yīng)系數(shù)為0.026 7，并且二者都通過(guò)了5%顯著性水平檢驗(yàn)。表明本地城市土地財(cái)政水平提升不僅會(huì)推動(dòng)本地城市土地城鎮(zhèn)化發(fā)展，這與劉顏[32]研究結(jié)論一致，而且還會(huì)通過(guò)正向溢出效應(yīng)推動(dòng)鄰近城市土地城鎮(zhèn)化發(fā)展。主要原因是土地財(cái)政是地方政府收入重要來(lái)源和緩解財(cái)政壓力重要手段[33]，政府通過(guò)低價(jià)征收土地高價(jià)賣出土地“以地生財(cái)”模式樂(lè)此不疲，土地財(cái)政的存在推動(dòng)了城市空間的擴(kuò)張，是驅(qū)動(dòng)土地城鎮(zhèn)化快速發(fā)展的根源。同時(shí)本地城市低價(jià)征地高價(jià)賣地的“低征高賣”的生財(cái)模式容易通過(guò)示范、模仿和學(xué)習(xí)效應(yīng)對(duì)鄰近城市進(jìn)行空間傳導(dǎo)，推動(dòng)了鄰近城市土地城鎮(zhèn)化發(fā)展。

（5）對(duì)外開(kāi)放水平直接效應(yīng)系數(shù)為-0.004 4，未通過(guò)10%顯著性水平檢驗(yàn)，溢出效應(yīng)系數(shù)為-0.056 2，通過(guò)了1%顯著性水平檢驗(yàn)。表明本地城市外商投資水平的提升對(duì)本地城市土地城鎮(zhèn)化并無(wú)顯著的影響，但會(huì)通過(guò)負(fù)向溢出效應(yīng)抑制鄰近城市土地城鎮(zhèn)化的發(fā)展，這與盧新海等[34]研究結(jié)論較為一致。主要原因是城市之間在外商引資上存在競(jìng)爭(zhēng)效應(yīng)，存在以低價(jià)出讓土地使用權(quán)為代價(jià)作為吸引外資的籌碼，增加了外資對(duì)土地資源占用，推動(dòng)了土地城鎮(zhèn)化的擴(kuò)張，但外資的進(jìn)駐也可能會(huì)憑借其先進(jìn)技術(shù)和管理經(jīng)驗(yàn)提高土地資源利用效率，并對(duì)本地土地資源消耗大的產(chǎn)業(yè)產(chǎn)生擠出效應(yīng)[35]，抑制土地城鎮(zhèn)化無(wú)序蔓延，而這其中錯(cuò)綜復(fù)雜的作用機(jī)制，造成對(duì)外開(kāi)放水平對(duì)本地城市土地城鎮(zhèn)化未通過(guò)統(tǒng)計(jì)學(xué)意義上的檢驗(yàn)。同時(shí)吸引外資多的往往是經(jīng)濟(jì)實(shí)力較強(qiáng)和地理區(qū)位條件較為優(yōu)越的地區(qū)[36-38]，本地城市吸引外資多就會(huì)阻斷經(jīng)濟(jì)實(shí)力較弱和地理區(qū)位條件較差的城市以犧牲土地資源為代價(jià)來(lái)達(dá)到引進(jìn)外資的目的，從而抑制了其他城市土地城鎮(zhèn)化的擴(kuò)張。

5 結(jié)論與啟示

本文采集中國(guó)283座城市2003—2017年面板數(shù)據(jù)，利用空間自相關(guān)和借助空間面板杜賓模型及其偏微分分解方法，從空間效應(yīng)視角對(duì)中國(guó)城市土地城鎮(zhèn)化驅(qū)動(dòng)因素進(jìn)行分析，得出以下主要結(jié)論與啟示。

（1）從時(shí)序特征來(lái)看，研究期內(nèi)中國(guó)城市土地城鎮(zhèn)化水平呈現(xiàn)出持續(xù)上升演進(jìn)特征，土地城鎮(zhèn)化城市之間差異呈持續(xù)擴(kuò)大特征，并伴隨著極化現(xiàn)象，土地城鎮(zhèn)化處于高值區(qū)城市數(shù)量比例呈現(xiàn)出擴(kuò)大特征。

（2）從空間特征來(lái)看，中國(guó)城市土地城鎮(zhèn)化在空間分布上并非隨機(jī)狀態(tài)，而是存在顯著正向全局空間自相關(guān)性，但這種全局空間關(guān)系尚未形成穩(wěn)態(tài)，年際間存在一定的波動(dòng)性。局域空間格局特征表現(xiàn)為高值集聚、低值集聚、低值塌陷和高值凸起四種集聚類型，數(shù)量上以HH和LL為主，以LH和HL為輔，并且存在“低者恒低、高者恒高”局域空間格局特征。

（3）從SPDM估計(jì)結(jié)果來(lái)看，中國(guó)城市土地城鎮(zhèn)化空間回歸系數(shù)在三種空間杜賓模型固定效應(yīng)估計(jì)下系數(shù)都顯著為正，本地城市土地城鎮(zhèn)化水平提升會(huì)促使其他鄰近城市土地城鎮(zhèn)化水平提升，城市土地城鎮(zhèn)化空間溢出效應(yīng)的存在，打破了傳統(tǒng)研究中基于樣本相互獨(dú)立的基本假設(shè)，將空間效應(yīng)納入城市土地城鎮(zhèn)化的計(jì)量回歸模型之中，才能合理解釋城市土地城鎮(zhèn)化驅(qū)動(dòng)因素。

（4）從SPDM空間效應(yīng)分解來(lái)看，經(jīng)濟(jì)發(fā)展水平對(duì)本地城市土地城鎮(zhèn)化具有推動(dòng)作用，對(duì)鄰近城市土地城鎮(zhèn)化具有負(fù)向溢出效應(yīng)。人口城鎮(zhèn)化水平對(duì)本地城市土地城鎮(zhèn)化并無(wú)統(tǒng)計(jì)學(xué)意義上的推動(dòng)作用，對(duì)鄰近城市土地城鎮(zhèn)化具有負(fù)向溢出效應(yīng)。第二產(chǎn)業(yè)結(jié)構(gòu)水平對(duì)本地城市土地城鎮(zhèn)化具有推動(dòng)作用，對(duì)鄰近城市土地城化并無(wú)統(tǒng)計(jì)學(xué)意義上溢出效應(yīng)。土地財(cái)政水平對(duì)本地城市土地城鎮(zhèn)化具有推動(dòng)作用，對(duì)鄰近城市土地城鎮(zhèn)化具有正向溢出效應(yīng)。對(duì)外開(kāi)放水平對(duì)本地城市土地城鎮(zhèn)化并無(wú)統(tǒng)計(jì)學(xué)意義上的推動(dòng)作用，對(duì)鄰近城市土地城鎮(zhèn)化具有負(fù)向溢出效應(yīng)。

（5）政策啟示如下：第一，城市土地城鎮(zhèn)化空間溢出效應(yīng)的存在決定了城市在土地資源開(kāi)發(fā)利用中不可“孤城而為”，在區(qū)域一體化發(fā)展背景下，綜合考慮本地城市與鄰近城市之間的空間互動(dòng)關(guān)系，加強(qiáng)城市之間在土地城鎮(zhèn)化上的聯(lián)動(dòng)管控作用[39-40]。第二，推進(jìn)國(guó)土空間治理能力現(xiàn)代化，推進(jìn)土地供給側(cè)結(jié)構(gòu)性改革，城市發(fā)展中需要摒棄以城市空間擴(kuò)展為主的外延式土地城鎮(zhèn)化，推進(jìn)以人為本的內(nèi)涵式人口城鎮(zhèn)化發(fā)展，加強(qiáng)對(duì)城市現(xiàn)有存量土地的挖潛。第三，推進(jìn)城市新舊動(dòng)能轉(zhuǎn)換和產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級(jí)，將城市經(jīng)濟(jì)增長(zhǎng)從“樓宇經(jīng)濟(jì)”中釋放出來(lái)，逐步擺脫對(duì)土地財(cái)政的依賴，加強(qiáng)對(duì)外商投資管理，避免城市之間以土地價(jià)格戰(zhàn)手段來(lái)引進(jìn)外資，利用其先進(jìn)管理技術(shù)和經(jīng)驗(yàn)，提高城市土地利用效率。第四，根據(jù)城市實(shí)際需求實(shí)施差異化的建設(shè)用地配給政策，在理性推進(jìn)土地城鎮(zhèn)化過(guò)程中既需要兼顧短期目標(biāo)與長(zhǎng)遠(yuǎn)目標(biāo)相統(tǒng)一，也需要綜合考慮經(jīng)濟(jì)效益、社會(huì)效益和生態(tài)效益與城市空間承載能力協(xié)調(diào)統(tǒng)一，從而為城市健康可持續(xù)發(fā)展提供基礎(chǔ)保障。

參考文獻(xiàn)

[1]新華社. 中共中央辦公廳 國(guó)務(wù)院辦公廳印發(fā)《關(guān)于在國(guó)土空間規(guī)劃中統(tǒng)籌劃定落實(shí)三條控制線的指導(dǎo)意見(jiàn)》[EB/OL].2019-11-01[2019-12-02]. http：//www.gov.cn/xinwen/2019-11/01/content_5447654.htm.

[2]盧斌瑩，陳正江，侯春紅，等.基于GIS的西安市南郊土地城市化空間發(fā)展研究[J].人文地理，2010，25（1）：56-59.

[3]楊洋，黃慶旭，章立玲.基于DMSP/OLS夜間燈光數(shù)據(jù)的土地城鎮(zhèn)化水平時(shí)空測(cè)度研究：以環(huán)渤海地區(qū)為例[J].經(jīng)濟(jì)地理，2015，35（2）：141-148，168.

[4]崔許鋒.民族地區(qū)的人口城鎮(zhèn)化與土地城鎮(zhèn)化：非均衡性與空間異質(zhì)性[J].中國(guó)人口·資源與環(huán)境，2014，24（8）：63-72.

[5]劉瓊，杜曉航，盛業(yè)旭.基于階段對(duì)比的中國(guó)人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)關(guān)系[J].中國(guó)人口·資源與環(huán)境，2018，28（1）：26-34.

[6]陳春，于立，張銳杰，等.中國(guó)城鎮(zhèn)化加速階段中期土地城鎮(zhèn)化與人口城鎮(zhèn)化的協(xié)調(diào)程度[J].長(zhǎng)江流域資源與環(huán)境，2016，25（11）：1654-1662.

[7]李寶禮，胡雪萍.金融發(fā)展會(huì)造成人口城鎮(zhèn)化滯后于土地城鎮(zhèn)化嗎：基于安徽省16個(gè)地級(jí)市面板數(shù)據(jù)的研究[J].華東經(jīng)濟(jì)管理，2014，28（12）：18-23.

[8]許芬.中國(guó)人口城鎮(zhèn)化滯后土地城鎮(zhèn)化的根源：基于城鄉(xiāng)間要素不平等交換視角的分析[J].城市問(wèn)題，2016（7）：12-17.

[9]李子聯(lián).人口城鎮(zhèn)化滯后于土地城鎮(zhèn)化之謎：來(lái)自中國(guó)省際面板數(shù)據(jù)的解釋[J].中國(guó)人口·資源與環(huán)境，2013，23（11）：94-101.

[10]王成新，王波濤，王翔宇.基于結(jié)構(gòu)視角的中國(guó)人口城市化與土地城市化異速增長(zhǎng)研究[J].中國(guó)人口·資源與環(huán)境，2016，26（8）：135-141.

[11]李永樂(lè)，舒幫榮，吳群.房地產(chǎn)價(jià)格與土地城鎮(zhèn)化：傳導(dǎo)機(jī)制與實(shí)證研究[J].中國(guó)土地科學(xué)，2013，27（11）：52-58.

[12]吳一凡，劉彥隨，李裕瑞.中國(guó)人口與土地城鎮(zhèn)化時(shí)空耦合特征及驅(qū)動(dòng)機(jī)制[J].地理學(xué)報(bào)，2018，73（10）：1865-1879.

[13]張立新，朱道林，杜挺，等.長(zhǎng)江經(jīng)濟(jì)帶土地城鎮(zhèn)化時(shí)空格局及其驅(qū)動(dòng)力研究[J].長(zhǎng)江流域資源與環(huán)境，2017，26（9）：1295-1303.

[14]劉耀林，李紀(jì)偉，侯賀平，等.湖北省城鄉(xiāng)建設(shè)用地城鎮(zhèn)化率及其影響因素[J].地理研究，2014，33（1）：132-142.

[15]王曉云，楊秀平，張雪梅.基于DEA-Tobit兩步法的城鎮(zhèn)化效率評(píng)價(jià)及其影響因素：從人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)發(fā)展的視角[J].生態(tài)經(jīng)濟(jì)，2017，33（5）：29-34.

[16]盧新海，柯楠，匡兵，等.中部地區(qū)土地城鎮(zhèn)化水平差異的時(shí)空特征及影響因素[J].經(jīng)濟(jì)地理，2019，39（4）：192-198.

[17]李慶.空間相關(guān)性對(duì)各省市生態(tài)文明建設(shè)的影響分析[J].中國(guó)人口·資源與環(huán)境，2019，29（9）：91-98.

[18]譚術(shù)魁，饒映雪，朱祥波.土地投入對(duì)中國(guó)經(jīng)濟(jì)增長(zhǎng)的影響[J].中國(guó)人口·資源與環(huán)境，2012，22（9）：61-67.

[19]TOBLER W R. A computer movie simulating urban growth in the Detroit region[J]. Economic geography， 1970， 46（1）：234-240.

[20]LESAGE J P， PACE R K. Introduction to spatial econometrics[M]. Boca Raton： CRC Press， 2009：45-75.

[21]孫焱林，覃飛，陳亞會(huì)，等.基于空間計(jì)量的地級(jí)市政府間城鎮(zhèn)化競(jìng)爭(zhēng)分析[J].中國(guó)人口·資源與環(huán)境，2018，28（5）：115-122.

[22]楊孟禹，張可云.中國(guó)城市擴(kuò)張的空間競(jìng)爭(zhēng)實(shí)證分析[J].經(jīng)濟(jì)理論與經(jīng)濟(jì)管理，2016（9）：100-112.

[23]李成宇，張士強(qiáng)，張偉.中國(guó)省際工業(yè)生態(tài)效率空間分布及影響因素研究[J].地理科學(xué)，2018，38（12）：1970-1978.

[24]邵朝對(duì)，蘇丹妮，鄧宏圖.房?jī)r(jià)、土地財(cái)政與城市集聚特征：中國(guó)式城市發(fā)展之路[J].管理世界，2016（2）：19-31，187.

[25]周琳，范建雙，虞曉芬.政府間競(jìng)爭(zhēng)影響城市土地市場(chǎng)化水平的雙邊效應(yīng)研究：基于財(cái)政競(jìng)爭(zhēng)和引資競(jìng)爭(zhēng)的不同作用[J].中國(guó)土地科學(xué)，2019，33（5）：60-68.

[26]ELHORST J P. Spatial econometrics： from cross-sectional data to spatial panels[M]. Heidelberg： Springer，2014： 63-67.

[27]伍駿騫，方師樂(lè)，李谷成，等.中國(guó)農(nóng)業(yè)機(jī)械化發(fā)展水平對(duì)糧食產(chǎn)量的空間溢出效應(yīng)分析 [J].中國(guó)農(nóng)村經(jīng)濟(jì)，2017（6）：44-57.

[28]高金龍，包菁薇，劉彥隨，等.中國(guó)縣域土地城鎮(zhèn)化的區(qū)域差異及其影響因素[J].地理學(xué)報(bào)，2018，73（12）：2329-2344.

[29]文貫中，熊金武.化地不化人的城市化符合中國(guó)國(guó)情嗎：人口密集型的“老浦西”和土地資本密集型的“新浦東”的歷史比較[J].城市規(guī)劃，2012，36（4）：18-24.

[30]朱莊瑞，藏波.房地產(chǎn)市場(chǎng)促進(jìn)城鎮(zhèn)化建設(shè)的作用機(jī)理與協(xié)調(diào)發(fā)展研究[J].中國(guó)人口·資源與環(huán)境，2016，26（9）：116-122.

[31]呂立剛，隋雪艷，汪翔，等.江蘇省土地城鎮(zhèn)化的空間分異及其主導(dǎo)因素探測(cè)[J].人文地理，2018，33（4）：88-94，112.

[32]劉顏.土地財(cái)政對(duì)中國(guó)城市空間擴(kuò)張的影響：基于動(dòng)態(tài)空間計(jì)量模型的實(shí)證檢驗(yàn)[J].財(cái)經(jīng)理論與實(shí)踐，2019，40（3）：9-14.

[33]安勇，原玉廷.土地財(cái)政、扭曲效應(yīng)與區(qū)域創(chuàng)新效率[J].中國(guó)土地科學(xué)，2019，33（8）：36-42，52.

[34]盧新海，唐一峰，匡兵.長(zhǎng)江中游城市群城市土地利用效率空間溢出效應(yīng)研究[J].長(zhǎng)江流域資源與環(huán)境，2018，27（2）：252-261.

[35]劉騫文，閆笑.地方政府“土地引資”背景下的FDI擠入擠出效應(yīng)研究[J].財(cái)經(jīng)研究，2016，42（1）：17-29.

[36]袁曉玲，呂文凱.從“資源引致”向“效率引致”：基于政府效率、引資優(yōu)惠及溢出效應(yīng)對(duì)FDI的影響分析[J].現(xiàn)代經(jīng)濟(jì)探討，2019（7）：10-18.

[37]張林，張軍龍，劉凱，等.環(huán)境約束下山東省城鎮(zhèn)化效率及影響因素[J].山東師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2019，34（3）：295-304.

[38]徐超.代合治.山東省新型城鎮(zhèn)化的時(shí)間演變與空間差異研究[J].山東師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2018，33（1）：105-111.

[39]張紅鳳，王鶴鳴，何旭.基于改進(jìn)引力模型的山東省城市空間聯(lián)系與格局劃分[J].山東財(cái)經(jīng)大學(xué)學(xué)報(bào)，2019，31（3）：110-120.

[40]孔艷芳.城鎮(zhèn)化是否縮小了中國(guó)城鄉(xiāng)收入差距：基于直接影響與空間溢出效應(yīng)的經(jīng)驗(yàn)論證[J].山東財(cái)經(jīng)大學(xué)學(xué)報(bào)，2019，31（4）：87-98.