薛建春，張安錄

引言

當前是“城鄉(xiāng)分割”走向“城鄉(xiāng)融合”發(fā)展的關鍵時期，“粗放的”土地城鎮(zhèn)化和“半集約式”的人口城鎮(zhèn)化成為中國城鎮(zhèn)化發(fā)展的主要特征。城鎮(zhèn)化作為一個復雜系統(tǒng)過程，土地資源的集約節(jié)約利用是其核心要義，土地城鎮(zhèn)化關注遷移數量，而人口城鎮(zhèn)化則相對關注遷移質量，無論人口或者土地都是生產建設的核心要素，城鎮(zhèn)化過程中雖然土地利用轉變了其基本屬性，但是用地效率不能得到提升，導致微觀上土地資源要素價格扭曲，宏觀上土地城鎮(zhèn)化、人口城鎮(zhèn)化與土地利用相脫節(jié)，人口要素變?yōu)槿丝趬毫?，城?zhèn)化發(fā)展必將陷入困局。新經濟常態(tài)下城鄉(xiāng)二元結構問題不僅體現在投資與收入差異，還涉及社會保障與土地資源同權同價等多方面，城鎮(zhèn)化發(fā)展不僅影響了生產要素的投入，而且還包含了要素投入后的“余值”，因此，從城鎮(zhèn)化發(fā)展的多個角度分析其對提升城市土地利用全要素生產率（Total Factor Productivity，TFP）的影響機制，可以破解城鎮(zhèn)化發(fā)展的內部牽制以及人口、土地要素投入后對城市土地利用全要素生產率的“余值”效應，尋找人口城鎮(zhèn)化與土地城鎮(zhèn)化協同發(fā)展的政策實施方向。

綜述現有文獻，就人口城鎮(zhèn)化與土地城鎮(zhèn)化而言，學者們多從指標測度、[1](p67-75)發(fā)展階段、[2](p12-18)驅動因素、[3](p138-143,p149)協調分析、[4](p67-75)[5](p1345-1354)競爭分析[6](p115-122)角度展開，然而城鎮(zhèn)化作為一個復雜開放系統(tǒng)，其發(fā)展過程直接影響著城市生產要素的流動，土地利用不局限于資源配置與利用效率本身，技術進步與規(guī)模效率更能促進土地利用全要素生產率的提升，當人口數量的激增不能匹配生產力的發(fā)展，城市系統(tǒng)生產率中的有效產出不能歸因于生產要素部分未被適當考慮，學者們開始注意到人口城鎮(zhèn)化和城市土地擴張對城市土地利用全要素生產率的影響分析。長江經濟帶農村勞動力轉移促進了林業(yè)全要素生產率的提升，[7](p1-15)新型城鎮(zhèn)化通過推動產業(yè)結構轉型、能源消費結構優(yōu)化顯著促進綠色全要素生產率的提升，[8](p116-119)城鎮(zhèn)化與物流綠色全要素生產率之間存在非線性關系，并在交通基礎設施提升過程中由負向影響轉為正向影響。[9](p123-132)還有學者從土地利用與管理制度出發(fā)，分析土地利用效率與城鎮(zhèn)化之間存在的協同與拮抗效應，[10](p27-32)單獨分析建設用地利用效率的空間非均衡性，[11](p154-159)探索城市蔓延對全要素生產率的影響機制與現實效應，[12](p90-100)這些研究為本文提供了理論基礎，但研究區(qū)域主要集中在人口稠密、土地富饒且稀少的東南地區(qū)，對于人少地廣且土地相對貧瘠的黃河流域研究較少，本研究以黃河流域69 個城市（州）為研究對象，通過厘清人口城鎮(zhèn)化與土地城鎮(zhèn)化之間影響城市土地利用全要素生產率的內部掣肘機制，探尋人口城鎮(zhèn)化與土地城鎮(zhèn)化對土地利用全要素生產率的影響效應，挖掘生態(tài)脆弱且人稀地廣的黃河流域城市經濟增長源泉，為政府制定高質量可持續(xù)發(fā)展政策提供合理建議。

一、理論分析與模型構建

城鎮(zhèn)化是各類生產要素在城鄉(xiāng)內部分配或者集聚的復雜過程，[13](p146-156)國內學者經常將城鎮(zhèn)化分為人口城鎮(zhèn)化和土地城鎮(zhèn)化進行研究，人口是城鎮(zhèn)化的核心，土地是城鎮(zhèn)化的載體，[14]工業(yè)發(fā)展與人口城鎮(zhèn)化、土地城鎮(zhèn)化之間的關聯效應直接影響區(qū)域環(huán)境狀況。[15](p145-149)土地利用全要素生產率是土地利用過程中借助投入產出分析得到的土地集約節(jié)約利用指標，是產出增長率超出要素投入增長率的“余值”部分，本研究中的土地利用全要素生產率不僅包含人口、土地投入要素，而且包含技術進步、能力實現、規(guī)模效應導致的產出增加。眾多學者[16](p54-61,p81)[17](p95-101)[18](p50-57)依據古典經濟學思想，選用土地、資本和勞動力作為投入要素計算土地利用全要素生產率，這也體現了城鎮(zhèn)化與土地利用之間的緊密關聯性。城鎮(zhèn)化對土地利用的影響機制簡述如下：

（一）人口城鎮(zhèn)化對城市土地利用全要素生產率的影響機制

人口城鎮(zhèn)化是農村人口遷移進入城市并失去耕作土地參與城鎮(zhèn)產業(yè)建設的過程，體現了城鎮(zhèn)化過程中人口的流動。第一，城市二三產業(yè)主要集中在城鎮(zhèn)區(qū)，隨著二三產業(yè)結構比例增高以及產業(yè)結構的升級，產業(yè)發(fā)展提供了大量的就業(yè)崗位，吸引農業(yè)人口向城鎮(zhèn)流動，導致二三產業(yè)人力投入提高，一方面增加了土地利用全要素生產率的投入要素量，另一方面增加了產業(yè)產值，帶動城市土地利用全要素生產率提升；第二，伴隨經濟新常態(tài)以及科技進步對農牧業(yè)的影響，農牧業(yè)耕種解放了更多的勞動力，也讓更多的傳統(tǒng)手工業(yè)者尋找到新的供給側需求，農村剩余勞動力開始流向有需求市場的城鎮(zhèn)，促進新產業(yè)新業(yè)態(tài)的積極發(fā)展，滿足供給側結構市場需要，優(yōu)化產業(yè)結構，引發(fā)經濟結構升級，提升城市土地利用全要素生產率；第三，中國實施多年的城鄉(xiāng)二元結構，導致城鎮(zhèn)人均收入持續(xù)高于農牧民收入，同時戶籍壁壘引發(fā)的社會保障差異，資源配置制度引發(fā)的城鄉(xiāng)資本差異，這些城鄉(xiāng)邊際產出差異對農村居民轉向城鎮(zhèn)形成“拉力”，[19](p69-70)城市基礎設施等在日漸增大的人口壓力下，資金投入逐步加大，從資本角度拉升了土地利用全要素生產率，人口城鎮(zhèn)化對土地利用全要素生產率產生引致效應；第四，日益便捷的交通設施以及交通成本也為人口流動提供了便捷，受城市虹吸效應影響，城市周邊居住的農民選擇城鎮(zhèn)工作，往返于城鄉(xiāng)之間，帶動土地利用全要素生產率。

（二）土地城鎮(zhèn)化對土地利用全要素生產率的影響機制

土地城鎮(zhèn)化是農業(yè)用地轉變?yōu)槌擎?zhèn)建設用地的轉型過程，衡量了土地利用轉移的效果。改革開放初期，我國土地城鎮(zhèn)化速度較快，隨著經濟的快速發(fā)展與習近平生態(tài)文明思想的提出，黨的十八大以后土地生態(tài)效益逐漸被大家認識，土地城鎮(zhèn)化速度趨緩。但不斷增加的城市建設用地需求導致城市建成區(qū)面積持續(xù)擴大，各類土地逐步向城市建設用地轉變，首先，土地作為基本的投入要素，城市建設用地擴張帶動產業(yè)發(fā)展和產業(yè)結構優(yōu)化，直接影響土地利用全要素生產率提升，同時，科學技術進步也保證了投入土地的產出增加；其次，在中央和地方政府形成的政治集權、經濟分權的委托—代理關系下，[20](p287-295)地方政府借助“土地發(fā)動機”實現城市生產總值提升和政府財政收益，使得土地征收功能異化，這種商業(yè)化的土地資源調控手段對土地利用全要素生產率產生拮抗效應。[21](p69-78,p107)城鎮(zhèn)化對土地利用全要素生產率的影響機制如圖1所示：

圖1 人口城鎮(zhèn)化和土地城鎮(zhèn)化對土地利用TFP影響機制

（三）模型構建

1.城市土地利用全要素生產率測度模型

借鑒學者們對土地利用全要素生產率測度的研究方法，[16](p54-61，p81)本研究選用包含非期望產出的Super-SBM 模型與Malmquist 指數相結合，選用窗口寬度為10、窗口偏移量為0 的參比方式計算黃河流域城市土地利用全要素生產率指數，這種設置方式可以保證每個決策單元計算時都參考同一前沿面，實現決策單元之間的有效比較，而且數據包絡方法（Data Envelopment Analysis，DEA）不需要給出具體的生產函數形式，同時還考慮了抑制土地利用全要素生產率的非期望因素，對于DEA有效的決策單元還可以繼續(xù)區(qū)分其生產率值的大小。其規(guī)劃模型表示如下：

公式中，ρj表示j城市的效率值，每個城市單元有m種投入和q種產出，xij和yrj分別為城市單元j自身的輸入和輸出，Xij和Yrj分別為城市單元整體的輸入和輸出，且表示過度投入和表示產出不足。

考慮到規(guī)模報酬可變條件下部分Super-SBM模型可能存在無可行解的問題，所以選用規(guī)模報酬不變條件下基于投入導向的Super-SBM 模型與Malmquist 指數模型結合計算黃河流域城市土地利用全要素生產率指數。根據RF^are提出的幾何平均值方法，全要素生產率指數可以進一步分解為規(guī)模效率變化（EC）和技術進步變化（TC），其公式表示為：

公式中，MI 表示土地利用全要素生產率指數，由于全要素生產率指數測算的是兩個生產率的比值，因此采用累乘的方法計算各決策單元各年的土地利用全要素生產率、規(guī)模效率和技術進步值進行后期分析。

2.城鎮(zhèn)化對城市土地利用全要素生產率的門檻效應

本研究從定量角度分析人口城鎮(zhèn)化和土地城鎮(zhèn)化對土地利用全要素生產率的影響效應，由于城鎮(zhèn)化與土地利用之間存在動態(tài)效應，且目前的人地掛鉤制度、城鄉(xiāng)建設用地增減掛鉤制度等都使得城鎮(zhèn)化與土地利用之間的關系愈加復雜，因此將城鎮(zhèn)化的二次項加入模型中，考察兩者之間是否存在非線性關系。設定模型如下：

公式中，tfpit表示i城市t時期的城市土地利用全要素生產率，luit表示i城市t時期土地城鎮(zhèn)化，puit表示i城市t時期人口城鎮(zhèn)化表示土地城鎮(zhèn)化與人口城鎮(zhèn)化的二次項，lnconit表示模型中的控制變量，εit表示模型的擾動項。

其次，為了識別人口城鎮(zhèn)化和土地城鎮(zhèn)化對土地利用全要素生產率影響的拐點，選用門檻模型作進一步分析，雙門檻模型設定如下：

公式中，γl、γh分別代表土地城鎮(zhèn)化的兩個門檻值，σl、σh分別代表人口城鎮(zhèn)化的兩個門檻值，ηit和δit表示門檻模型的擾動項，其他變量與公式（4）中相同。

二、變量選取與數據來源

（一）研究區(qū)域概況

黃河流域自西向東途經青海、四川、甘肅、寧夏、內蒙古、山西、陜西、河南、山東9個省區(qū)69個城市（州），流域內69 個城市（州）的常住人口21749.9萬人，占全國人口總量的15.54%，土地面積218.3萬平方公里，占國土面積22.74%，2020 年流域內城市人口城鎮(zhèn)化均值41.11%，低于國家水平43.37%，GDP 總量133754.4 億元，占全國經濟總量的14%，人均GDP 水平低于全國水平。2019 年習近平總書記號召推動黃河流域高質量發(fā)展，堅持生態(tài)優(yōu)先和綠色發(fā)展，這對流域內的城鎮(zhèn)化建設與土地利用提出了全新挑戰(zhàn)，一方面土地是城鎮(zhèn)建設的載體，但黃河流域內未利用土地21.72 萬平方公里，約占流域面積的11%，雖地廣但可利用性差；另一方面人口城鎮(zhèn)化要求城鎮(zhèn)土地具有更高經濟承載力與生態(tài)承載力，但現實是流域內人才流失較嚴重，甚至直接影響了經濟與生態(tài)建設。因此，研究黃河流域69個城市（州）人口城鎮(zhèn)化、土地城鎮(zhèn)化對土地利用全要素生產率的影響效應，可以更深入理解人口城鎮(zhèn)化與土地城鎮(zhèn)化的協同發(fā)展方向，以及對黃河流域高質量發(fā)展的推進。

（二）變量選取

1.城市土地利用全要素生產率測度

根據古典經濟學思想，選擇土地投入、資本投入和勞動力投入作為土地利用全要素生產率測算的投入要素，其中土地投入指標選用城市建設用地面積（包括城鎮(zhèn)村及工礦用地和交通運輸用地）表示，資本投入指標選用地均固定資產投資總額表示，勞動力投入指標選用第二、三產業(yè)從業(yè)人數表示，產出要素包含期望產出與非期望產出，二、三產業(yè)產值表示期望產出的經濟效益，[22]在崗職工年平均工資表示期望產出的社會效益，土地利用碳排放量表示非期望產出的生態(tài)效益。

為了簡化計算，土地利用碳排放量利用系數代入法計算，以城市能源消費量作為建設用地碳排放量基礎數據，采用能源消耗總量與每噸標準煤綜合碳排放系數的乘積計入，而每噸標準煤碳排放系數則主要考慮10 種能源類型，煤炭、焦炭、石油、原油、汽油、煤油、柴油、燃料油、液化石油氣和天然氣，通過計算研究區(qū)域各年的能源消耗碳排放系數①此處省略了研究區(qū)間內根據69 個城市（州）10 種能源消耗量占比作為權重得出的綜合碳排放系數過程。，最終確定系數為0.743t/t標準煤；其他地類碳排放量主要計入耕地、園地、林地、草地、水域和其他未利用地，各類能源、土地碳排放系數如表1所示。而且，利用GDP指數將各類價值指標平減到以2011年為基期的數值，減少價格對指標的影響。

表1 不同類型能源、土地碳排放系數

2.城鎮(zhèn)化變量與控制變量

本研究采用2011—2020 年黃河流域9 省69 個城市（州）的面板數據作為研究樣本，主要解釋變量說明如下：

（1）城鎮(zhèn)化水平

盡管復雜的復合序參考體系可以提高估計參量的測算精度和穩(wěn)定性，但為了避免割裂參量與研究內涵之間的直接聯系，本研究選用非農人口占城市常住人口比重表示人口城鎮(zhèn)化水平（pu），選用城市建設用地面積占行政區(qū)劃土地面積比重表示土地城鎮(zhèn)化水平（lu）。

（2）控制變量

①消費水平（consum），用城市人均年消費支出表示。②產業(yè)結構（ratio），考慮我國目前正處于工業(yè)化發(fā)展階段，因此選用二產業(yè)產值與三產業(yè)產值比重表示城市產業(yè)結構的優(yōu)化情況，值越小產業(yè)結構越高級。③經濟增長（bgdp），采用城市每年總GDP與行政區(qū)劃土地面積比值表示。④政府規(guī)模（brevenue），用城市的財政總收入與行政區(qū)劃土地面積比值表示。對所有經濟變量按2011 年不變價進行平減，各變量的描述性統(tǒng)計見表2。

（三）數據來源

考慮到MaxDEA 軟件對數據要求的完整性，城市土地利用全要素生產率測度所需的投入、產出變量數據均來源于黃河流域各地區(qū)的《城市統(tǒng)計年鑒》《能源統(tǒng)計年鑒》《中國國土資源年鑒》和土地利用變化數據（2011—2020），統(tǒng)計分析中采用的各變量數據主要來源于Wind 數據庫和EPS 數據庫，個別缺失數據使用插值法獲得。

三、結果分析

（一）黃河流域城市土地利用全要素生產率結果分析

2011—2020 年黃河流域內城市土地利用全要素生產率呈現增長趨勢，10 年間69 個城市（州）平均值1.504，且2011—2012年間增速較緩慢，2013—2020 年增速較快，技術進步也呈現增長態(tài)勢，但規(guī)模效率則顯示微弱下降，說明土地利用全要素生產率增長主要依賴技術進步而非規(guī)模效率。從各城市10 年間土地利用全要素生產率的平均值分析，僅有2 個城市低于1.0，處于非有效狀態(tài)，其他67 個城市（州）均處于有效狀態(tài)，而且黃河下游鄭州市的土地利用全要素生產率最大（2.8918）；從規(guī)模效率的平均值分析，35個城市（州）處于規(guī)模效率有效狀態(tài)，且主要分布在山西、陜西、山東和上游的青海省內，但下游菏澤市的規(guī)模效率值最大（1.4066）；從技術進步的平均值分析，所有城市的技術進步值均大于1.0，而且技術進步的增長明顯高于全要素生產率的增長，鄭州市的技術進步值最高，更進一步說明技術進步有效促進了城市土地利用全要素生產率。（見圖2）

圖2 黃河流域內城市土地利用全要素生產率、規(guī)模效率、技術進步各年平均值

（二）黃河流域人口城鎮(zhèn)化與土地城鎮(zhèn)化演變分析

用ArcGIS分別繪制2011年和2020年的人口城鎮(zhèn)化與土地城鎮(zhèn)化時空差異圖（圖3），可以發(fā)現黃河流域上游的人口城鎮(zhèn)化水平與土地城鎮(zhèn)化水平低于中下游，特別是黃河“幾字彎”城市的人口城鎮(zhèn)化發(fā)展顯示減弱狀態(tài)，土地城鎮(zhèn)化發(fā)展最快的是黃河下游的山東省。10 年中69 個城市（州）的人口城鎮(zhèn)化平均值為34.27%，土地城鎮(zhèn)化平均值8.42%。其中，烏海市的人口城鎮(zhèn)化水平最高（86.53%），定西市的人口城鎮(zhèn)化水平最低（12.54%），這是因為內蒙古烏海市是一個僅設區(qū)不設縣的地級城市，因此人口城鎮(zhèn)化水平相對較高。鄭州市的土地城鎮(zhèn)化水平最高（27.73%），玉樹藏族自治州的土地城鎮(zhèn)化水平最低（0.1%），說明黃河下游城市對土地的開發(fā)強度遠高于上游地區(qū)。

圖3 黃河流域2011、2020年人口城鎮(zhèn)化、土地城鎮(zhèn)化時空變化

從時間趨勢看，10 年間有58 個城市（州）的人口城鎮(zhèn)化水平處于上升態(tài)勢，11個城市（州）則處于下降態(tài)勢，這主要是因為城市人口的凈流出導致，比較典型的如內蒙古的包頭、鄂爾多斯市，青海省的海東、海北藏族自治州、黃南藏族自治州、果洛藏族自治州，陜西的渭南、延安等，甚至寧夏的首府銀川，也出現人口流失現象。這些城市大多沒有良好的大學資源和強有力的人才引進政策，為了尋找更好的工作機會與發(fā)展機會，外出就讀的青年放棄回鄉(xiāng)工作，更直接加重出生率降低，人口數量銳減。從流域人口城鎮(zhèn)化整體水平看，2011—2020年黃河流域城市人口城鎮(zhèn)化均值從35.75% 提升至41.11%，增長5.51%低于同期全國人口城鎮(zhèn)化增長率12.6%。反觀流域內土地城鎮(zhèn)化水平，所有城市均處于上升態(tài)勢，10年間黃河流域城市建設用地面積總共增長了7300km2，其中，鄭州市2011—2020年建設用地增長了約410km2，占鄭州市行政區(qū)劃面積的5.49%；因此，土地城鎮(zhèn)化是黃河流域城市城鎮(zhèn)化的主要結果，造成這種局面，一方面因為土地是重要的生產要素，土地擴張可以帶動GDP 增長，另一方面由于政府間財政競爭，各城市爭相效仿出讓土地滿足財政需求，沒有考慮與之相適應的人口城鎮(zhèn)化匹配，因此“鬼城”應運而生。[23](p160)

為了更進一步說明人口城鎮(zhèn)化與土地城鎮(zhèn)化之間的發(fā)展關系，借鑒城市用地增長彈性系數概念，利用城鎮(zhèn)用地增長率與城鎮(zhèn)人口增長率的比值表示，并采用中國城市規(guī)劃設計研究院在《2000 年城鎮(zhèn)用地預測綜合報告》中提出的1.12控制值進行衡量。[24]10 年間黃河流域內城市用地增長彈性系數僅有3 個地區(qū)（海南藏族自治州、蘭州和東營市）大于1.12，顯示土地城鎮(zhèn)化發(fā)展過快，其他66 個城市的人口、土地城鎮(zhèn)化發(fā)展均處于較合理范圍。

（三）城鎮(zhèn)化對土地利用全要素生產率的影響效應分析

首先利用面板數據選擇混合效應模型、固定效應模型和隨機效應模型進行分析，Hausman檢驗顯示卡方統(tǒng)計量在1%水平顯著拒絕隨機效應的原假設，城鎮(zhèn)化對土地利用全要素生產率、技術進步回歸采用固定效應模型，對規(guī)模效率回歸采用隨機效應模型（表3）。

表3 2011—2020年城鎮(zhèn)化對土地利用全要素生產率、技術效率、技術進步的回歸結果

1.流域內總體效應分析

無論混合效應模型、隨機效應模型還是固定效應模型，人口城鎮(zhèn)化與土地利用全要素生產率均呈現倒“U”型關系，土地城鎮(zhèn)化與土地利用全要素生產率呈現“U”型關系，且都在1%水平顯著；根據Hausman檢驗結果，選用固定效應模型系數分析變量之間的影響關系，比較而言，土地城鎮(zhèn)化對土地利用全要素生產率的影響效果要大于人口城鎮(zhèn)化，而且土地城鎮(zhèn)化和人口城鎮(zhèn)化對土地利用全要素生產率呈現相反影響效果，1%顯著水平下土地城鎮(zhèn)化每增加1%，土地利用全要素生產率減少0.2716，人口城鎮(zhèn)化每增加1%，土地利用全要素生產率增加0.0603。從各控制變量分析，人均消費水平與經濟增長促進土地利用全要素生產率增長，產業(yè)結構與政府規(guī)模則抑制土地利用全要素生產率增長，模型中所有控制變量對土地利用全要素生產率的影響效應均在1%水平顯著，說明土地利用全要素生產率存在消費效應和經濟效應，而且，城市中二產產值與三產產值比重越大，越形成抑制效應。

2.城鎮(zhèn)化對土地利用規(guī)模效率與技術進步的影響分析

繼續(xù)分析人口城鎮(zhèn)化和土地城鎮(zhèn)化對土地利用規(guī)模效率和技術進步的影響，Hausman 檢驗顯示1%水平顯著拒絕對技術進步的隨機效應假設，但無法顯著拒絕對規(guī)模效率隨機效應模型的原假設，因此對規(guī)模效率的回歸分析采用隨機效應模型，對技術進步的回歸分析采用固定效應模型。表3 結果顯示，人口城鎮(zhèn)化和土地城鎮(zhèn)化對規(guī)模效率無顯著影響，但是從系數值分析，人口城鎮(zhèn)化與土地利用規(guī)模效率、技術進步均呈現倒“U”型關系，僅土地城鎮(zhèn)化與技術進步呈現“U”型關系，且土地城鎮(zhèn)化相較于人口城鎮(zhèn)化對土地利用技術進步的影響效果更強烈，土地城鎮(zhèn)化每增加1%技術進步下降0.3315。同時，人均消費水平與經濟增長在1%水平顯著促進技術進步，產業(yè)結構與政府規(guī)模也在1%水平顯著抑制技術進步，這與對土地利用全要素生產率的影響效果相同，更進一步說明黃河流域土地利用全要素生產率提升主要是因為技術進步引起。

3.城鎮(zhèn)化對土地利用全要素生產率的區(qū)域差異分析

為了更深入分析黃河流域城鎮(zhèn)化對不同區(qū)域土地利用全要素生產率的影響效應，將黃河流域劃分為上中下游三個區(qū)域，上游、中游、下游分別有25、26 和18 個城市（州），分析其區(qū)域差異性。從表4 可以發(fā)現，黃河流域上游的人口城鎮(zhèn)化對土地利用全要素生產率沒有顯著影響，土地城鎮(zhèn)化與土地利用全要素生產率在1%水平顯著呈現倒“U”型關系；中游城市的人口城鎮(zhèn)化、土地城鎮(zhèn)化均與土地利用全要素生產率在1%水平顯著且呈現倒“U”型關系；下游城市的人口城鎮(zhèn)化與土地利用全要素生產率在1%水平顯著且呈現倒“U”型關系，土地城鎮(zhèn)化與土地利用全要素生產率在1%水平顯著且呈現“U”型關系，可見流域內城鎮(zhèn)化對土地利用全要素生產率影響存在區(qū)域差異性。下游城市的土地城鎮(zhèn)化水平（17.70%）明顯高于上中游（3.45%、6.77%），雖然建設用地投入增多，但由于土地利用全要素生產率是考慮了包含投入要素及以外的其他因素，因此建設用地擴張并沒有直接促進土地利用全要素生產率的增長，同時也說明土地城鎮(zhèn)化可能是一個門檻變量，當其高于某個門檻值的時候會導致人口城鎮(zhèn)化對土地利用全要素生產率呈現相反影響效果。再從各控制變量分析，上中游城市的消費效應顯著促進土地利用全要素生產率，但下游影響效果不顯著，經濟增長在中下游顯著促進土地利用全要素生產率，但在上游城市中卻顯示抑制作用，上中下游產業(yè)結構（二產比重增加）均顯著抑制土地利用全要素生產率，政府規(guī)模則僅在下游顯著抑制土地利用全要素生產率?？傮w而言，流域內各區(qū)域存在土地利用全要素生產率的消費效應正向影響與產業(yè)結構工業(yè)化的逆向影響。

表4 黃河流域上中下游城鎮(zhèn)化對土地利用全要素生產率影響效應差異分析

（四）人口城鎮(zhèn)化與土地城鎮(zhèn)化影響土地利用全要素生產率的門檻效應

1.人口城鎮(zhèn)化影響土地利用全要素生產率存在土地城鎮(zhèn)化的雙門檻效應

城鎮(zhèn)化發(fā)展始于城鎮(zhèn)建設用地擴張，2003年以來，隨著西部大開發(fā)、中部崛起等戰(zhàn)略實施，中國土地供給開始傾向東北和中西部地區(qū)，土地城鎮(zhèn)化增長速度快于人口城鎮(zhèn)化，二者的協同發(fā)展也成為學者們研究的重心，[25](p1654-1662)[26](p129-135)但是二者并行發(fā)展過程中對土地利用全要素生產率是否存在門檻效應，很少有學者探究，基于前述城鎮(zhèn)化對土地利用全要素生產率影響的區(qū)域差異性，本節(jié)利用模型（5）探索人口城鎮(zhèn)化影響土地利用全要素生產率是否存在土地城鎮(zhèn)化的門檻效應。門檻檢驗首先探究不同門檻區(qū)間內的參數空間是否顯著不同，其次借助Bootstrap 方法對門檻的估計值和實際值進行一致性檢驗，[27](p15-28)分別通過F統(tǒng)計量的P值與似然比LR值進行檢驗，檢驗結果見表5。

表5 門檻效應檢驗與門檻值、置信區(qū)間

檢驗結果顯示，F統(tǒng)計量與P值均拒絕不存在門檻效應的原假設，因此，人口城鎮(zhèn)化影響土地利用全要素生產率過程中，土地城鎮(zhèn)化存在兩個門檻值，分別為0.48%與20.37%。為了更好理解門檻值的估計與置信區(qū)間的構造過程，選用LR=7.3523 表示95%的置信水平，繪制門檻模型的似然比函數圖，圖4更直觀顯示存在兩個門檻值。

圖4 門檻估計值與置信區(qū)間

表6 給出了門檻回歸的結果，人口城鎮(zhèn)化與土地利用全要素生產率之間的關系被土地城鎮(zhèn)化這一門檻變量劃分為三個區(qū)間，且不同區(qū)間影響差異顯著，當土地城鎮(zhèn)化處于較低水平時，人口城鎮(zhèn)化對土地利用全要素生產率在1%水平顯著負向影響，說明土地城鎮(zhèn)化水平較低時，人口城鎮(zhèn)化抑制土地利用全要素生產率增長，這是因為人口壓力在有限的土地資源上造成承載力不足，導致土地利用全要素生產率不能提升。當土地城鎮(zhèn)化逐步提升邁過第一門檻值后，人口城鎮(zhèn)化對土地利用全要素生產率的影響轉為正向促進作用，但在統(tǒng)計上不顯著。這是因為當建設用地擴張后，流入城市的人口不會成為承載壓力，更多功能是作為勞動力要素的投入，因此促進了城市土地利用全要素生產率的提升。但是當土地城鎮(zhèn)化越過第二個門檻值后，人口城鎮(zhèn)化對土地利用全要素生產率的影響系數在10%水平顯著為負，抑制作用代替了第二階段的促進作用，主要原因是隨著土地城鎮(zhèn)化的快速發(fā)展，一方面要求人口數量的快速遷移，另一方面對人口城鎮(zhèn)化發(fā)展的質量提出更高要求，如果流入城市的人口僅是數量上的增加沒有質量上的提高，人口的壓力功能大于勞動力要素功能抑制了土地利用全要素生產率提升。為了保證研究結論的可靠性，本研究采用兩種方式進行穩(wěn)健性檢驗，首先考慮替換解釋變量的方法進行穩(wěn)健性檢驗，將土地城鎮(zhèn)化指標用市轄區(qū)建設用地面積與行政區(qū)劃面積比重表示，得到表6中的穩(wěn)健性檢驗（I）。其次，考慮增加模型中的控制變量，看回歸結果是否穩(wěn)健，檢驗結果見表6 中的穩(wěn)健性檢驗（II）。可以發(fā)現，兩種穩(wěn)健性檢驗結果都顯示人口城鎮(zhèn)化影響土地利用全要素生產率時存在土地城鎮(zhèn)化雙門檻效應。分析各控制變量的影響效果，消費水平與經濟增長依舊顯著促進城市土地利用全要素生產率，產業(yè)結構（二產/三產比重）則抑制土地利用全要素生產率。

表6 黃河流域人口城鎮(zhèn)化影響土地利用全要素生產率的門檻效應及穩(wěn)健性檢驗

2.土地城鎮(zhèn)化影響土地利用全要素生產率存在人口城鎮(zhèn)化的單門檻效應

同理，分析土地城鎮(zhèn)化影響土地利用全要素生產率是否存在人口城鎮(zhèn)化的門檻效應。采用模型（6）進行回歸分析，結果顯示存在單門檻效應（表7）。當人口城鎮(zhèn)化低于68.31%時，土地城鎮(zhèn)化對土地利用全要素生產率的正向促進作用為0.085，當人口城鎮(zhèn)化高于68.31%時，土地城鎮(zhèn)化對土地利用全要素生產率的促進作用增強。同樣采用替換解釋變量的方法和增加控制變量方法進行穩(wěn)健性檢驗，穩(wěn)健性檢驗（I）用非農人口與總人口的比值替換人口城鎮(zhèn)化，穩(wěn)健性檢驗（II）增加固定資產投資控制變量，可以發(fā)現人口城鎮(zhèn)化的單門檻效應依舊存在。

表7 黃河流域土地城鎮(zhèn)化影響土地利用全要素生產率的門檻效應及穩(wěn)健性檢驗

四、結論與建議

通過分析人口城鎮(zhèn)化、土地城鎮(zhèn)化與土地利用全要素生產率的影響機制，初步判定城鎮(zhèn)化與土地利用全要素生產率之間存在非線性關系，隨后利用固定效應模型和門檻模型定量分析，得出如下結論：

（一）2011—2020 年黃河流域城市土地利用全要素生產率呈現增長趨勢，且2013年后增長速度更快，97.1%的城市土地利用全要素生產率大于1.0，且土地利用全要素生產率主要由于技術進步引起。今后土地利用過程中，不僅要考慮土地資源的合理配置，更要注重技術進步對土地利用的提升影響。

（二）2011—2020 年間黃河流域全部城市的土地城鎮(zhèn)化水平處于上升態(tài)勢，10年里建設用地面積總共增長了7300km2；11 個城市（州）的人口城鎮(zhèn)化水平處于下降趨勢，說明流域內城市存在人口流失現象，但從建設用地彈性系數分析，絕大部分城市處于合理水平。

（三）人口城鎮(zhèn)化與土地利用全要素生產率之間存在倒“U”型關系，土地城鎮(zhèn)化與土地利用全要素生產率之間存在“U”型關系，證實了城鎮(zhèn)化與土地利用全要素生產率之間存在非線性關系，而且人口城鎮(zhèn)化、土地城鎮(zhèn)化與土地利用全要素生產率之間的非線性關系存在區(qū)域差異性，黃河流域上、中游內人口城鎮(zhèn)化、土地城鎮(zhèn)化與土地利用全要素生產率呈現倒“U”型關系，但下游土地城鎮(zhèn)化與土地利用全要素生產率則呈現“U”型關系。

（四）以土地城鎮(zhèn)化作為門檻變量，分析人口城鎮(zhèn)化對土地利用全要素生產率的影響關系，發(fā)現人口城鎮(zhèn)化對土地利用全要素生產率的非線性影響存在雙門檻效應，土地城鎮(zhèn)化初級階段，人口城鎮(zhèn)化抑制土地利用全要素生產率提升，土地城鎮(zhèn)化進入中級階段以后，人口城鎮(zhèn)化促進土地利用全要素生產率提升，當土地城鎮(zhèn)化越過第二個門檻值后，人口城鎮(zhèn)化再次呈現對土地利用全要素生產率的負向影響。但是分析土地城鎮(zhèn)化對土地利用全要素生產率的影響過程中，人口城鎮(zhèn)化僅存在單門檻效應，且影響方向相同。

基于以上分析結論得出如下政策啟示：

（一）積極推進黃河流域人口城鎮(zhèn)化高質量發(fā)展，控制土地城鎮(zhèn)化高速發(fā)展。首先，人口和土地作為土地利用的投入要素，如果僅追求數量的提升，不能滿足土地利用過程中對技術進步的需求，人口要素轉化為人口壓力；其次，土地城鎮(zhèn)化發(fā)展速度過快（超越第二門檻值），人口城鎮(zhèn)化不再促進土地利用全要素生產率，因此，流域內可以借助發(fā)展高端產業(yè)，努力提升教育水平，創(chuàng)造更多就業(yè)機會，營造良好生態(tài)環(huán)境，吸引大學生創(chuàng)業(yè)，集聚高層次人才，引導新增城鎮(zhèn)人口高質量發(fā)展。同時，積極盤活現有城鎮(zhèn)低效利用土地和存量土地，避免土地城鎮(zhèn)化過速發(fā)展。

（二）建立以常住人口為標準的“人地掛鉤”思想，城鎮(zhèn)化過程中匹配適宜的建設用地指標，預留生態(tài)用地，避免“人地脫鉤”。當城鎮(zhèn)建設用地增加規(guī)模與新增常住人口規(guī)模實現良性互動時，[28](p143-149)建設用地彈性系數趨于合理，土地城鎮(zhèn)化水平處于兩個門檻值中間，保證人口城鎮(zhèn)化對土地利用全要素生產率的促進效用。

（三）加強黃河流域產業(yè)結構高級化建設。黃河流域三產產值相比二產產值比重偏低，通過改善營商環(huán)境，發(fā)展現代服務業(yè)、生態(tài)旅游業(yè)等，加強城鄉(xiāng)溝通性基礎設施建設，帶動城鄉(xiāng)之間的互動，實現產業(yè)結構的升級，提升居民的消費能力，最終縮小城鄉(xiāng)收入差距，全面促進經濟增長，提高城市土地利用全要素生產率。