栗向陽　王磊　張宇欣

摘要：在我國(guó)經(jīng)濟(jì)發(fā)展逐漸轉(zhuǎn)向高質(zhì)量增長(zhǎng)階段的背景下，研究市域單元城鎮(zhèn)化的協(xié)調(diào)發(fā)展程度對(duì)轉(zhuǎn)變城市經(jīng)濟(jì)增長(zhǎng)方式和推動(dòng)城市高質(zhì)量發(fā)展具有重要的現(xiàn)實(shí)意義.基于2003—2019年中國(guó)285個(gè)城市的面板數(shù)據(jù)，利用耦合協(xié)調(diào)度模型和探索性時(shí)空數(shù)據(jù)分析方法從城鎮(zhèn)化水平和城鎮(zhèn)化質(zhì)量角度對(duì)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的時(shí)空格局特征進(jìn)行識(shí)別和刻畫，并借助空間面板杜賓模型分析我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的驅(qū)動(dòng)因素及其空間效應(yīng).研究表明：研究期間內(nèi)，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平較低但呈持續(xù)上升特征，城際差異呈持續(xù)縮小特征，空間格局呈“東高西低”和“南高北低”的地帶性特征.Morans I指數(shù)表明中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展呈現(xiàn)出空間集聚分布，且集聚程度隨時(shí)間不斷增強(qiáng)；時(shí)間路徑分析表明城鎮(zhèn)化協(xié)調(diào)發(fā)展局部空間結(jié)構(gòu)和空間依賴方向上穩(wěn)定性較強(qiáng)，其空間格局演變呈動(dòng)態(tài)整合的特征；時(shí)空躍遷分析表明城鎮(zhèn)化協(xié)調(diào)發(fā)展具有明顯的路徑依賴和空間鎖定特征.經(jīng)濟(jì)發(fā)展、人口規(guī)模、人力資本、科技創(chuàng)新、產(chǎn)業(yè)結(jié)構(gòu)、政府能力對(duì)本地城市城鎮(zhèn)化協(xié)調(diào)發(fā)展均具有正向直接效應(yīng)；人力資本、科技創(chuàng)新和產(chǎn)業(yè)結(jié)構(gòu)對(duì)鄰域城鎮(zhèn)化協(xié)調(diào)發(fā)展具有負(fù)向溢出效應(yīng)，人口規(guī)模對(duì)鄰域城鎮(zhèn)化協(xié)調(diào)發(fā)展具有正向溢出效應(yīng).在未來城市發(fā)展過程中，要綜合考慮城鎮(zhèn)化協(xié)調(diào)發(fā)展的區(qū)域差異和空間溢出效應(yīng)，因地制宜地制定推動(dòng)城鎮(zhèn)化高質(zhì)量發(fā)展的指導(dǎo)意見和政策方針.

關(guān)鍵詞：城鎮(zhèn)化協(xié)調(diào)發(fā)展；城鎮(zhèn)化水平和質(zhì)量；耦合協(xié)調(diào)度模型；空間計(jì)量模型

中圖分類號(hào)：F 126文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1001－988Ⅹ（2024）02－0110－13

Spatial－temporal evolution of coupled and coordinatedurbanization and its influencing factors

—From the perspective of urbanization level and quality

LI Xiang－yang WANG Lei ，ZHANG Yu－xin3

Abstract：Under the background that Chinas economic development is gradually turning to the stage of high－quality growth，it is of great practical significance to study the coordinated development degree of urban unit urbanization for transforming the mode of urban economic growth and promoting the high－quality development of the city.Based on the panel data of 285 cities in China from 2003 to 2019，the coupling coordination degree model and exploratory spatial－temporal data analysis method

is used to identify and characterize the spatial－temporal pattern characteristics of coupled and coordinated urbanization development level from the perspective of urbanization level and urbanization quality，and analyzes the driving factors and spatial effects of coupled and coordinated urbanization development in China with the help of spatial panel Durbin model.The results show that：During the study period，the coupled coordinated development of Urbanization in China is at a low level，but it continues to rise，and the inter－city differences continue to narrow.The spatial pattern shows a zonal characteristic of“high in the east and low in the west”and“high in the south and low in the north”.Morans I index shows that the coupled coordinated development of Urbanization in China presents a spatial agglomeration distribution，and the agglomeration degree increases with time.The time path analysis shows that the local spatial structure and spatial dependence of coordinated urbanization development are stable，and its spatial pattern evolution shows the dynamic characteristics of spatial integration.The spatial－temporal transition analysis shows that the coordinated development of urbanization has obvious characteristics of path dependence and space lock.Economic development，population size，human capital，scientific and technological innovation，industrial structure and government capacity have positive and direct effects on the coordinated development of local urbanization.Human capital，scientific and technological innovation and industrial structure have negative spillover effects on the coordinated development of neighborhood urbanization，while population size has positive spillover effects on the coordinated development of neighborhood urbanization.In the process of future urban development，regional differences and spatial spillover effects of coordinated urbanization development should be considered comprehensively，and guidelines and policies for promoting high－quality urbanization development should be formulated according to local conditions.

Key words：urban coordinated development;level and quality of urbanization;coupling coordination degree model;spatial econometric model

城鎮(zhèn)化是一個(gè)伴隨著人口結(jié)構(gòu)轉(zhuǎn)變、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、土地利用規(guī)模擴(kuò)張、經(jīng)濟(jì)增長(zhǎng)方式轉(zhuǎn)變和城市環(huán)境質(zhì)量提升的復(fù)雜動(dòng)態(tài)過程［1－2］.其水平和質(zhì)量的提升，不僅有利于拉動(dòng)消費(fèi)內(nèi)需增長(zhǎng)，實(shí)現(xiàn)我國(guó)循環(huán)經(jīng)濟(jì)的持續(xù)健康發(fā)展；也有助于推動(dòng)工業(yè)化、信息化和農(nóng)業(yè)現(xiàn)代化的協(xié)同發(fā)展，加快推進(jìn)社會(huì)主義現(xiàn)代化進(jìn)程.自改革開放以來，伴隨著工業(yè)化的推進(jìn)，我國(guó)的城鎮(zhèn)化發(fā)展也取得了巨大成就.若以城鎮(zhèn)常住人口與社會(huì)總?cè)丝诘谋戎刈鳛楹饬恐笜?biāo)，我國(guó)的城鎮(zhèn)化率從1978年的17.92%提升到2020年的63.89%，年均增速1.095個(gè)百分點(diǎn)；城鎮(zhèn)常住人口也由1.7億增加到9.02億［2－4］.從全球范圍看，我國(guó)的城鎮(zhèn)化水平也遠(yuǎn)超過世界銀行公布的55.3%的世界平均標(biāo)準(zhǔn).中國(guó)的城鎮(zhèn)化發(fā)展，不僅決定著其自身經(jīng)濟(jì)發(fā)展動(dòng)態(tài)和走向，也決定著世界城鎮(zhèn)化發(fā)展的未來和前景，影響著世界經(jīng)濟(jì)的發(fā)展格局［5］.經(jīng)濟(jì)學(xué)家Stiglitz［6－7］就曾高度評(píng)價(jià)過中國(guó)的城鎮(zhèn)化發(fā)展，認(rèn)為其與以美國(guó)為首的新技術(shù)革命是影響21世紀(jì)世界經(jīng)濟(jì)發(fā)展最為深刻的事件.另一方面，由于國(guó)情和發(fā)展的需要，我國(guó)在追求城鎮(zhèn)化規(guī)?？焖贁U(kuò)張的同時(shí)忽略了城鎮(zhèn)化的發(fā)展質(zhì)量，使得我國(guó)的城鎮(zhèn)化發(fā)展出現(xiàn)了水平與質(zhì)量不相協(xié)調(diào)的現(xiàn)象，表現(xiàn)為城鎮(zhèn)化各子系統(tǒng)失衡發(fā)展（例如土地城鎮(zhèn)化快于人口城鎮(zhèn)化）、城市空間結(jié)構(gòu)分布不合理、城市管理和治理水平效率低下，進(jìn)而導(dǎo)致城市無序擴(kuò)張、生態(tài)環(huán)境惡化、交通嚴(yán)重?fù)矶碌纫幌盗小俺鞘胁　眴栴}.

當(dāng)前，我國(guó)正處于社會(huì)經(jīng)濟(jì)轉(zhuǎn)型升級(jí)的重要時(shí)期，城鎮(zhèn)化的發(fā)展也逐漸由加速推進(jìn)階段過渡到高質(zhì)量發(fā)展階段.為了克服城鎮(zhèn)化快速發(fā)展過程中存在的發(fā)展不平衡、發(fā)展不協(xié)調(diào)的問題，積極發(fā)揮城鎮(zhèn)化對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展的引擎作用，黨中央和國(guó)務(wù)院也提出了若干政策和建議.2014年3月，國(guó)務(wù)院印發(fā)《國(guó)家新型城鎮(zhèn)化規(guī)劃（2014—2020年）》，指出要走以人為本、布局高效、生態(tài)文明的新型城鎮(zhèn)化道路，不斷提升城市的可持續(xù)發(fā)展水平，并積極推進(jìn)城鎮(zhèn)化水平和質(zhì)量的穩(wěn)步提升和協(xié)同發(fā)展.黨的“十八大”、“十九大”和“二十大”報(bào)告連續(xù)提出要通過綠色低碳、集約高效的發(fā)展方式來持續(xù)提升我國(guó)的城鎮(zhèn)化發(fā)展質(zhì)量.《中華人民共和國(guó)國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展第十四個(gè)五年規(guī)劃和2035年遠(yuǎn)景目標(biāo)綱要》更是明確要求“要加快轉(zhuǎn)變城市發(fā)展方式，推進(jìn)新型城市建設(shè)，統(tǒng)籌城市規(guī)劃管理，提高城市治理水平，通過實(shí)施城市更新行動(dòng)，推動(dòng)城市空間結(jié)構(gòu)優(yōu)化和品質(zhì)提升”.“十四五”規(guī)劃期間將是我國(guó)推進(jìn)新型城鎮(zhèn)化高質(zhì)量發(fā)展的關(guān)鍵時(shí)期，在此背景下，如何協(xié)調(diào)城鎮(zhèn)化的發(fā)展已成為學(xué)者們關(guān)注的焦點(diǎn).

縱觀現(xiàn)有文獻(xiàn)，關(guān)于城鎮(zhèn)化協(xié)調(diào)發(fā)展的研究主要集中在以下幾個(gè)方面.一是城鎮(zhèn)化內(nèi)部各子系統(tǒng)的耦合協(xié)調(diào)發(fā)展.學(xué)者們從人口城鎮(zhèn)化、土地城鎮(zhèn)化和經(jīng)濟(jì)城鎮(zhèn)化等方面出發(fā)，分析城鎮(zhèn)化各系統(tǒng)的協(xié)調(diào)匹配程度、演變模式和時(shí)空差異［2，8－12］，并著重探究了經(jīng)濟(jì)發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)、政府行政、戶籍管制、財(cái)政體制、土地制度等因素對(duì)人口和土地城鎮(zhèn)化失衡發(fā)展的影響機(jī)制［13－18］，進(jìn)而提出城鎮(zhèn)化協(xié)調(diào)發(fā)展的優(yōu)化路徑和策略［19－21］.二是城鎮(zhèn)化與農(nóng)業(yè)化、工業(yè)化、生態(tài)環(huán)境等相關(guān)系統(tǒng)的耦合協(xié)調(diào)發(fā)展.鄧宗兵等［22］利用耦合協(xié)調(diào)度模型研究了長(zhǎng)江經(jīng)濟(jì)帶生態(tài)文明建設(shè)和新型城鎮(zhèn)化協(xié)調(diào)發(fā)展的時(shí)空分布差異，并從對(duì)外開放程度、經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)等方面探討了影響機(jī)制；王玉娟等［23］以全國(guó)284個(gè)地級(jí)及以上城市為例，通過研究發(fā)現(xiàn)新型城鎮(zhèn)化和低碳發(fā)展水平具有彼此強(qiáng)化的內(nèi)生互動(dòng)性，并且這種交互性受內(nèi)生和外生驅(qū)動(dòng)因素的影響非常顯著；李建豹等［24］從江蘇省的城鎮(zhèn)化和CO2排放水平的耦合協(xié)調(diào)水平進(jìn)行分析，發(fā)現(xiàn)其協(xié)調(diào)度隨著時(shí)間演進(jìn)呈上升趨勢(shì)，并且“北高南低，東高西低”的空間格局較為穩(wěn)定.目前，關(guān)于城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的研究已取得了豐碩成果，但在以下幾個(gè)方面仍存在不足，總結(jié)如下：① 現(xiàn)有學(xué)者主要從人口、經(jīng)濟(jì)、土地、生態(tài)等方面對(duì)中國(guó)城鎮(zhèn)化發(fā)展進(jìn)行研究，并重點(diǎn)關(guān)注人口城鎮(zhèn)化、土地城鎮(zhèn)化等城鎮(zhèn)化子系統(tǒng)的耦合協(xié)調(diào)性，較少?gòu)某擎?zhèn)化水平和質(zhì)量維度探究城鎮(zhèn)化子系統(tǒng)之間的相互關(guān)系.② 作為推動(dòng)區(qū)域協(xié)調(diào)發(fā)展和城鄉(xiāng)融合發(fā)展的重要引擎，城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展是一個(gè)時(shí)空維度并行的過程，而針對(duì)其路徑演化的研究往往將時(shí)間和空間維度割裂開，需要進(jìn)一步關(guān)注城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的時(shí)空關(guān)聯(lián)和時(shí)空交互特征.③ 已有研究更多地關(guān)注城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的測(cè)度和評(píng)價(jià)，而對(duì)影響城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的因素和機(jī)制探討不足.有鑒于此，本文從城鎮(zhèn)化的規(guī)模和質(zhì)量?jī)蓚€(gè)方面出發(fā)，將城鎮(zhèn)化系統(tǒng)分解為城鎮(zhèn)化水平和城鎮(zhèn)化質(zhì)量?jī)蓚€(gè)子系統(tǒng).在此基礎(chǔ)上，以中國(guó)285個(gè)地級(jí)及以上城市為研究對(duì)象，通過引入耦合協(xié)調(diào)度模型、探索性時(shí)空數(shù)據(jù)分析框架和空間計(jì)量模型等方法，分析中國(guó)城鎮(zhèn)化發(fā)展的耦合關(guān)系及其時(shí)空演變格局，并探究城鎮(zhèn)化耦合協(xié)調(diào)水平的影響因素和空間效應(yīng)，旨在為推動(dòng)我國(guó)新型城鎮(zhèn)化的高質(zhì)量發(fā)展提供有益參考和借鑒.

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

1.1 研究方法

1.1.1 城鎮(zhèn)化水平和質(zhì)量的測(cè)度

上述兩個(gè)方面對(duì)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展水平進(jìn)行分析，并將城鎮(zhèn)化發(fā)展劃分為城鎮(zhèn)化水平和質(zhì)量?jī)蓚€(gè)子系統(tǒng).其中，城鎮(zhèn)化水平從數(shù)量上反映了城鎮(zhèn)化發(fā)展情況；參考相關(guān)研究［3，25］，本文采用城鎮(zhèn)人口占市轄區(qū)總?cè)丝诘谋戎睾统擎?zhèn)建設(shè)用地占市轄區(qū)總土地的比重來刻畫城鎮(zhèn)化水平這一子系統(tǒng).作為一個(gè)內(nèi)涵豐富的綜合性概念，城鎮(zhèn)化質(zhì)量的具體涵義包括經(jīng)濟(jì)發(fā)展質(zhì)量、居民生活質(zhì)量、社會(huì)發(fā)展質(zhì)量、公共服務(wù)質(zhì)量和生態(tài)環(huán)境質(zhì)量等五個(gè)方面.其中經(jīng)濟(jì)發(fā)展質(zhì)量是城鎮(zhèn)化發(fā)展的物質(zhì)基礎(chǔ)，居民生活質(zhì)量是城鎮(zhèn)化發(fā)展的目的，社會(huì)發(fā)展質(zhì)量、公共服務(wù)質(zhì)量和生態(tài)環(huán)境質(zhì)量是城鎮(zhèn)化發(fā)展的保障.因此在參考相關(guān)研究［1，3，14－25］的基礎(chǔ)上本文從上述五個(gè)方面選取18項(xiàng)指標(biāo)衡量城鎮(zhèn)化質(zhì)量發(fā)展子系統(tǒng)（具體的指標(biāo)選取參見表1）.本文采用熵值法測(cè)算各城市的城鎮(zhèn)化水平和質(zhì)量.首先采用極差法對(duì)各指標(biāo)的原始數(shù)據(jù)進(jìn)行歸一化處理，其次利用熵值法計(jì)算各指標(biāo)的信息熵并確定各指標(biāo)的權(quán)重，最后采用多目標(biāo)加權(quán)求和法分別計(jì)算各城市的城鎮(zhèn)化水平和質(zhì)量.

其中，uit表示城市i在t年的城鎮(zhèn)化水平和質(zhì)量，其值越大，表示城市的城鎮(zhèn)化水平和質(zhì)量越高，反之則表示城市的城鎮(zhèn)化水平和質(zhì)量越低;pijt為城市i在t年第j項(xiàng)指標(biāo)標(biāo)準(zhǔn)化處理后的值;wj為第j指標(biāo)的權(quán)重.

1.1.2 耦合協(xié)調(diào)度模型

耦合協(xié)調(diào)度是測(cè)度兩個(gè)及以上系統(tǒng)之間相互作用和聯(lián)系程度的指標(biāo)，本文利用城鎮(zhèn)化水平和質(zhì)量耦合分析城鎮(zhèn)化的協(xié)調(diào)發(fā)展水平［26－27］.構(gòu)建耦合協(xié)調(diào)度模型如下：

其中，A為耦合度，u1表示城鎮(zhèn)化水平，u2表示城鎮(zhèn)化質(zhì)量；T為綜合協(xié)調(diào)指數(shù)，a和b分別為城鎮(zhèn)化水平和質(zhì)量對(duì)其貢獻(xiàn)份額，一般認(rèn)為在城鎮(zhèn)化進(jìn)程中，兩者的重要程度相當(dāng)，故取a=b=0.5；C表示城鎮(zhèn)化水平和質(zhì)量的耦合協(xié)調(diào)度，0≤C≤1，值越大，城鎮(zhèn)化的協(xié)調(diào)發(fā)展水平越高，借鑒相關(guān)研究［25，28－29］，本研究將城鎮(zhèn)化水平和質(zhì)量的協(xié)調(diào)等級(jí)分為10類.

1.1.3 探索性時(shí)空數(shù)據(jù)分析（ESTDA）

基于地理空間分異規(guī)律的研究在時(shí)間和空間維度上常處于割裂狀態(tài)，為了彌補(bǔ)這一缺陷，Rey等學(xué)者綜合考慮了時(shí)空因素的耦合性和連續(xù)性特征并提出了探索性時(shí)空數(shù)據(jù)分析（ESTDA）框架［30］.本文引入ESTDA分析框架并利用全局（Morans I）空間自相關(guān)、局部空間自相關(guān)（LISA）、LISA時(shí)間路徑、時(shí)空躍遷等方法對(duì)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的時(shí)空交互特征進(jìn)行分析.利用Morans I指數(shù)衡量城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的空間集聚狀況.通過LISA時(shí)間路徑揭示城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的局部時(shí)空動(dòng)態(tài)演變特征和時(shí)空交互變化特征，其特征量包括相對(duì)長(zhǎng)度（Γi）、彎曲度（Δi）和移動(dòng)方向（θi）.計(jì)算公式如下［31］：

其中，n為研究單元數(shù)量，本文取值為285；d（Li，t，Li，t+1）為研究單元i在莫蘭散點(diǎn)圖中從第t年到第t+1年的移動(dòng)距離，d（Li，1，Li，T）為研究單元i從研究初期到末期的移動(dòng)距離；Γi為時(shí)間路徑相對(duì)長(zhǎng)度，其值越大表明城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的局部空間結(jié)構(gòu)動(dòng)態(tài)性越強(qiáng)；Δi為時(shí)間路徑彎曲度，其值越大城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展在變化過程中的波動(dòng)性越強(qiáng)，相鄰區(qū)域?qū)Τ鞘衖的空間溢出或極化作用越明顯；θi為時(shí)間路徑移動(dòng)方向，劃分為4種類型，其中，0～90°方向（90～270°方向）表示城市及相鄰地區(qū)的耦合協(xié)調(diào)發(fā)展水平呈正（負(fù)）向協(xié)同增長(zhǎng)， 90～180°方向（270～360°方向）表示本地城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展水平增長(zhǎng)趨勢(shì)較低（高）而其相鄰城市的增長(zhǎng)趨勢(shì)較高（低）.

借助時(shí)空躍遷方法探究我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展局部鄰域空間關(guān)系的時(shí)空轉(zhuǎn)移特征，并將其分為Ⅰ、Ⅱ、Ⅲ和Ⅵ四種類型，分類標(biāo)準(zhǔn)和依據(jù)參照相關(guān)研究［32－36］.則城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的空間凝聚度可表示為：

其中，F(xiàn)0，t表示在t研究時(shí)段內(nèi)，本地區(qū)與鄰域均未發(fā)生躍遷的城市數(shù)量；m為可能發(fā)生躍遷的數(shù)量，在本研究中取值為16×285=4560；St取值為［0，1］，其值越大，表明城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的空間穩(wěn)定性越強(qiáng)，路徑依賴性特征越明顯.

1.1.4 空間計(jì)量模型

城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展水平具有較強(qiáng)的空間效應(yīng)，而傳統(tǒng)的計(jì)量分析方法未將這種空間交互效應(yīng)納入到模型中，會(huì)導(dǎo)致估計(jì)結(jié)果出現(xiàn)偏差［37－38］.本文以城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展水平為被解釋變量，構(gòu)建如公式9所示的空間計(jì)量模型：

其中，W為空間權(quán)重矩陣，本文基于Rook法則構(gòu)建鄰接矩陣，并以反距離權(quán)重矩陣作為替代權(quán)重進(jìn)行穩(wěn)健性分析；λ為空間誤差系數(shù)，ρ為空間自回歸系數(shù)；Wuit為空間誤差項(xiàng)，WCODit為被解釋變量的空間滯后項(xiàng)；CODit為城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展程度；GDPit為地區(qū)經(jīng)濟(jì)發(fā)展水平，采用城市人均GDP增速表示；POPit為城市人口規(guī)模，以各城市人口密度的自然對(duì)數(shù)表征；HUMit為人力資本，以萬人普通高等學(xué)校在校學(xué)生數(shù)和專任教師數(shù)通過熵值法［39－40］計(jì)算求得；INOVit為城市創(chuàng)新水平，利用萬人專利授權(quán)數(shù)的自然對(duì)數(shù)測(cè)算；STRUit為產(chǎn)業(yè)結(jié)構(gòu)水平，采用第二、三產(chǎn)業(yè)產(chǎn)值占地區(qū)生產(chǎn)總值的比重表示；GOVit為政府能力，以人均地方財(cái)政預(yù)算收入和支出利用熵值法求取；εit為隨機(jī)擾動(dòng)項(xiàng).

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

考慮到數(shù)據(jù)的可獲得性并消除新冠疫情對(duì)本研究可能造成的干擾，本文選取全國(guó)285個(gè)地級(jí)及以上城市2003—2019年的面板數(shù)據(jù)作為研究樣本，所涉及的數(shù)據(jù)包括社會(huì)經(jīng)濟(jì)數(shù)據(jù)和地理空間數(shù)據(jù).社會(huì)經(jīng)濟(jì)數(shù)據(jù)主要來源于各省、直轄市的統(tǒng)計(jì)年鑒、中國(guó)城市統(tǒng)計(jì)年鑒以及EPS數(shù)據(jù)平臺(tái)，缺失數(shù)據(jù)采用加權(quán)平均法和趨勢(shì)外推法補(bǔ)齊.地理空間數(shù)據(jù)包括行政邊界數(shù)據(jù)和夜間燈光遙感數(shù)據(jù).其中，各級(jí)行政邊界數(shù)據(jù)來源于國(guó)家基礎(chǔ)地理信息中心.夜間燈光數(shù)據(jù)則來源于哈佛大學(xué)數(shù)據(jù)庫(kù)中的“類NPP－VIIRS夜間燈光數(shù)據(jù)集”［41］，融合了DMSP－OLS和NPP－VIIRS兩種夜間燈光遙感數(shù)據(jù)，并經(jīng)過幾何校正、過飽和校正、時(shí)序校正和數(shù)據(jù)合成處理，能夠長(zhǎng)時(shí)間監(jiān)測(cè)城鎮(zhèn)化動(dòng)態(tài)變化過程.參考相關(guān)文獻(xiàn)［42－43］，計(jì)算城市燈光平均亮度值和燈光面積占轄區(qū)總面積的比重，并將其分別作為城鎮(zhèn)化質(zhì)量和水平的替代指標(biāo)以進(jìn)行穩(wěn)健性檢驗(yàn).相關(guān)變量的描述性統(tǒng)計(jì)結(jié)果如表2所示.

2 城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的時(shí)空特征

2.1 城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的時(shí)空演化格局

利用耦合協(xié)調(diào)度模型測(cè)算2003—2019年中國(guó)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展水平，并利用ArcGIS進(jìn)行可視化分析得到其空間分布圖（圖1）.2003—2019年，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展程度處于較低水平，但整體呈現(xiàn)波動(dòng)上升態(tài)勢(shì).從各年份協(xié)調(diào)程度的均值來看，2003年中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)程度的平均水平為0.1424，2019年其值上升為0.2844，年均提升6.23%；從各等級(jí)城市數(shù)量分布來講，2003年中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展程度在中等及以上水平的城市有5個(gè)（以0.5為劃分依據(jù)），占全部城市的比重為1.75%，2019年中等及以上水平的城市增加至24個(gè)，占比達(dá)到8.42%.長(zhǎng)期以來，我國(guó)城鎮(zhèn)化的發(fā)展以規(guī)模擴(kuò)張為主，使得以高質(zhì)高效為特征的城鎮(zhèn)化質(zhì)量的提升遠(yuǎn)遠(yuǎn)滯后于城鎮(zhèn)化水平的增強(qiáng)，進(jìn)而導(dǎo)致城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展一直處于較低水平.2003—2019年，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的空間分布變化趨勢(shì)明顯且整體呈現(xiàn)“東高西低”和“南高北低”的地帶性分異特征.研究期間內(nèi)，城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的高值區(qū)表現(xiàn)出“零星散點(diǎn)狀分布—沿海帶狀分布”的變化特點(diǎn)，高值點(diǎn)集中分布在珠三角、長(zhǎng)三角和京津冀等地區(qū)，主要包括深圳、東莞、上海、廈門、蘇州、中山等城市；耦合協(xié)調(diào)發(fā)展水平的低值區(qū)呈片狀分布，低值點(diǎn)主要分布在我國(guó)的甘肅、寧夏、四川、云南、黑龍江等省份和地區(qū)，其中隴南、固原、中衛(wèi)、張掖、安康等城市的城鎮(zhèn)化不協(xié)調(diào)程度尤為凸顯.從區(qū)域差異來看，研究期間內(nèi)，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的Theil指數(shù)由研究初期的0.4609下降至期末的0.2907，年均下降2.31%，表明我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的區(qū)域差異整體呈現(xiàn)波動(dòng)縮小態(tài)勢(shì).究其原因，伴隨著我國(guó)城鎮(zhèn)化質(zhì)量的不斷提升，城鎮(zhèn)化耦合協(xié)調(diào)程度整體也得到提升，使得城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的總體差異有所下降；東、中、西部三個(gè)地區(qū)及地區(qū)間的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展差異對(duì)我國(guó)整體區(qū)域差異的貢獻(xiàn)度由大到小依次為東部地區(qū)>西部地區(qū)>中部地區(qū)>地區(qū)間，說明地區(qū)內(nèi)差異是造成城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展整體差異的主要原因，同時(shí)也表明我國(guó)東部地區(qū)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的區(qū)內(nèi)差異最大，西部地區(qū)次之，中部地區(qū)最小.

2.2 城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的時(shí)空交互特征

2.2.1 空間自相關(guān)分析

基于Rook鄰近規(guī)則構(gòu)建空間權(quán)重矩陣，并利用ArcGIS軟件計(jì)算2003—2019年中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的Morans I指數(shù)（表3）.研究期間內(nèi)，城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的全局Morans I指數(shù)均為正值，且通過了1%的顯著性檢驗(yàn)，表明中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平在全域范圍內(nèi)具有較強(qiáng)的正相關(guān)性，即在空間上趨向于集聚分布.從時(shí)序上來看，2003—2019年城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的Morans I指數(shù)不斷增大，2003年Morans I值為0.50，2019年增大為0.60，說明城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的“馬太效應(yīng)”愈加明顯，其空間依賴程度和空間集聚水平逐步增強(qiáng).

2.2.2 LISA時(shí)間路徑分析

利用公式分別計(jì)算中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平時(shí)空路徑的相對(duì)長(zhǎng)度、彎曲度及移動(dòng)方向，并采用自然斷點(diǎn)法對(duì)其進(jìn)行分級(jí)，結(jié)果如圖2所示.由相對(duì)長(zhǎng)度來看，2003—2019年，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平相對(duì)長(zhǎng)度大于1的城市數(shù)量為107個(gè)，占城市總數(shù)的37.54%，說明城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的整體空間結(jié)構(gòu)較為穩(wěn)定；相對(duì)長(zhǎng)度整體呈現(xiàn)以城市群（圈）為核心向外圍遞減的變化趨勢(shì)，表明城市群尤其是中心城市城鎮(zhèn)化耦合協(xié)調(diào)程度的局部空間結(jié)構(gòu)具有較強(qiáng)的動(dòng)態(tài)性，究其原因，城市群的地理位置、交通區(qū)位、資源稟賦、政策條件等較為優(yōu)越，其城鎮(zhèn)化質(zhì)量的提升相較于其它地區(qū)也尤為明顯，因此這些地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的提升也最快.從彎曲度來看，研究期間內(nèi)，全國(guó)285個(gè)城市的彎曲度均大于1，表明中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)程度在發(fā)展過程中具有較強(qiáng)的空間依賴性；其中，呈現(xiàn)較高彎曲度的城市（前10位）主要分布在我國(guó)中部地區(qū)及河北、山東、陜西、廣西等省份，如延安、青島、安康、陽江、石家莊等，表明這些地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)水平的變化在空間依賴方向上具有很強(qiáng)的波動(dòng)性，易受鄰域空間極化或溢出效應(yīng)的影響；彎曲度較低的城市有東莞、齊齊哈爾、宿州和菏澤等，這些地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)水平具有較強(qiáng)的路徑依賴特征，受臨近區(qū)域空間作用的影響較小，其空間依賴變遷過程較為穩(wěn)定.從移動(dòng)方向上分析，2003—2019年，時(shí)間路徑移動(dòng)方向在0～90°方向和90～270°方向的城市有230個(gè)，占城市總數(shù)的80.70%，表明中國(guó)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展以協(xié)同增長(zhǎng)為主，整體呈現(xiàn)空間整合的動(dòng)態(tài)性特征；其中，正向協(xié)同增長(zhǎng)的城市有113個(gè)，占協(xié)同增長(zhǎng)城市數(shù)的49.13%，主要分布在東部沿海及安徽省、河南省、湖北省、貴州省、重慶市和四川省的部分地區(qū)，這些城市的城鎮(zhèn)化耦合協(xié)調(diào)水平呈現(xiàn)協(xié)同高增長(zhǎng)特征；負(fù)向協(xié)同增長(zhǎng)的城市有117個(gè)，占協(xié)同增長(zhǎng)城市數(shù)的50.87%，集中分布于西部?jī)?nèi)陸和東北地區(qū)，該地區(qū)呈現(xiàn)協(xié)同低增長(zhǎng)特征，同時(shí)也說明我國(guó)的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平仍處于低速增長(zhǎng)階段.

2.2.3 LISA時(shí)空躍遷分析

本研究利用時(shí)空躍遷和轉(zhuǎn)移概率矩陣分析中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展局部空間關(guān)聯(lián)類型的轉(zhuǎn)移特征和演變過程.由表4可知，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展局部空間類型未發(fā)生躍遷（Ⅰ型）的數(shù)量為4355，占全部可能發(fā)生躍遷數(shù)量的95.5%，即城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展局部Morans I的空間凝聚度0.955，局部空間類型發(fā)生時(shí)空躍遷的數(shù)量為205，占全部可能發(fā)生躍遷數(shù)量的4.5%，局部Morans I的空間離散度為0.045，表明我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的局部區(qū)域結(jié)構(gòu)極為穩(wěn)定，其分布具有較強(qiáng)的空間鎖定特征；其中，在發(fā)生躍遷的城市中，各類型時(shí)空躍遷的數(shù)量依次為Ⅱ型（102）、Ⅲ型（97）、Ⅵ型（6），發(fā)生的概率分別為0.022，0.021和0.001，即城市自身發(fā)生躍遷而鄰域未發(fā)生的情形在我國(guó)最顯著，其次是城市自身未發(fā)生躍遷而鄰域發(fā)生躍遷的情景，而城市自身和鄰域均發(fā)生躍遷的情形比較稀少.具體來看，在Ⅰ型中時(shí)空躍遷概率最大的是HH→HH（0.972）和LL→LL（0.974），說明大部分城市穩(wěn)定地保持在“協(xié)調(diào)型”或“失衡型”的空間關(guān)聯(lián)狀態(tài)下；在Ⅱ型中時(shí)空躍遷概率最大的是LH→HH（0.089），反映出部分城市由“洼地型”向“協(xié)調(diào)型”轉(zhuǎn)變，城市的耦合協(xié)調(diào)水平在發(fā)展過程中受到鄰域的影響逐漸趨同于周邊城市；在Ⅲ型中時(shí)空躍遷概率最大的是LH→LL（0.074），反映出部分城市的耦合協(xié)調(diào)發(fā)展水平由“洼地型”向“失衡型”轉(zhuǎn)變；在Ⅵ型中時(shí)空躍遷概率最大的是LL→HH（0.002），隨著城鎮(zhèn)化的推進(jìn)，部分地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)水平與鄰域協(xié)同提升，局部空間關(guān)聯(lián)狀態(tài)由“失衡型”轉(zhuǎn)變?yōu)椤皡f(xié)調(diào)型”.

3 城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的影響因素及空間效應(yīng)分析

3.1 模型選擇

本文采用空間計(jì)量模型進(jìn)行分析.進(jìn)一步地，利用LM和穩(wěn)健性LM檢驗(yàn)識(shí)別空間關(guān)聯(lián)形式，LM、robust LM檢驗(yàn)的統(tǒng)計(jì)量均通過了1%水平的顯著性檢驗(yàn)，表明空間誤差項(xiàng)和空間滯后項(xiàng)可能同時(shí)存在，傾向于選擇空間杜賓模型；同時(shí)，LR檢驗(yàn)和Wald檢驗(yàn)的統(tǒng)計(jì)量均在1%水平的顯著，表明SDM模型不能退化為SEM模型或者SAR模型，因此選擇空間杜賓模型.此外，Hausman檢驗(yàn)和聯(lián)合顯著性檢驗(yàn)的檢驗(yàn)結(jié)果也表明時(shí)間和個(gè)體固定效應(yīng)同時(shí)存在，所以本研究最終選用時(shí)空雙向固定效應(yīng)空間杜賓模型對(duì)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平的影響因素及其空間效應(yīng)進(jìn)行分析.

3.2 影響因素分析

空間面板杜賓模型（SPDM）的回歸分析結(jié)果如表6，被解釋變量空間滯后項(xiàng)的回歸系數(shù)ρ為正（0.102），并通過了1%的顯著性檢驗(yàn)，表明城鎮(zhèn)化的耦合協(xié)調(diào)水平具有顯著的正向空間溢出效應(yīng)，即城市自身城鎮(zhèn)化耦合協(xié)調(diào)程度的提升會(huì)促進(jìn)其鄰域城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展.同時(shí)這一分析結(jié)果也驗(yàn)證了將空間交互效應(yīng)納入到計(jì)量模型的合理性和必要性，采用空間計(jì)量模型分析具有良好的穩(wěn)健性.本研究對(duì)影響城鎮(zhèn)化耦合協(xié)調(diào)水平的因素進(jìn)行分析，并將其分解為直接效應(yīng)和間接效應(yīng)，結(jié)果如表7所示.

經(jīng)濟(jì)發(fā)展水平的直接效應(yīng)系數(shù)為0.006，在1%水平上顯著，間接效應(yīng)系數(shù)為-0.001，未通過顯著性水平檢驗(yàn)，表明城市本身經(jīng)濟(jì)發(fā)展水平的提升會(huì)促進(jìn)本地城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展，這與金丹等［9］、周艷等［11］、張鵬巖［28］等的研究結(jié)論一致，但尚未對(duì)鄰域城鎮(zhèn)化耦合協(xié)調(diào)程度形成顯著的空間溢出作用.究其原因，經(jīng)濟(jì)增長(zhǎng)是城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展空間格局形成的內(nèi)生驅(qū)動(dòng)力，一方面，經(jīng)濟(jì)的發(fā)展必然伴隨著外來資金的注入和大規(guī)模的人口集聚，使得居住用地、工業(yè)用地等城市建設(shè)用地規(guī)模不斷擴(kuò)張，進(jìn)而導(dǎo)致城鎮(zhèn)化水平持續(xù)擴(kuò)大；另一方面，經(jīng)濟(jì)增長(zhǎng)會(huì)增加政府社會(huì)固定資產(chǎn)投資和對(duì)公共設(shè)施的投資，提高教育醫(yī)療、文體娛樂、交通通信等公共服務(wù)設(shè)施水平，持續(xù)推動(dòng)城市環(huán)境和質(zhì)量提升，吸引人才等優(yōu)質(zhì)資源向城市集聚，進(jìn)一步使得城鎮(zhèn)化質(zhì)量不斷提升；整體來說，經(jīng)濟(jì)發(fā)展會(huì)通過改變城鎮(zhèn)化水平和質(zhì)量的相互關(guān)系來影響城市自身的城鎮(zhèn)化耦合協(xié)調(diào)程度.

人口規(guī)模的直接效應(yīng)系數(shù)為0.073，間接效應(yīng)系數(shù)為0.034，兩者均至少通過了10%水平的顯著性檢驗(yàn)，表明人口規(guī)模的擴(kuò)張不僅會(huì)推動(dòng)本地城市的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展，這與范擎宇等［10］的研究結(jié)論較為一致，并且對(duì)鄰域城鎮(zhèn)化耦合協(xié)調(diào)程度具有正向溢出效應(yīng).主要原因是人口的規(guī)模集聚會(huì)倒逼城鎮(zhèn)化發(fā)展走集約化道路，在一定程度上改變長(zhǎng)期以來我國(guó)存在的“人口城鎮(zhèn)化滯后于土地城鎮(zhèn)化”的不合理現(xiàn)象以及“重體量不重質(zhì)量”的城鎮(zhèn)化發(fā)展模式，改善城市土地利用方式、有效提升土地利用效率，進(jìn)而通過提升城鎮(zhèn)化質(zhì)量來促進(jìn)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展.對(duì)于鄰域地區(qū)，由于地域范圍內(nèi)存在人口的自由流動(dòng)和再分布，城市自身人口規(guī)模的擴(kuò)大會(huì)間接提高鄰近城市的人口集聚程度，從而提升鄰域地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平.但整體來講，城市人口集聚效應(yīng)對(duì)其本身的城鎮(zhèn)化耦合協(xié)調(diào)水平的帶動(dòng)作用更為明顯.

人力資本的直接效應(yīng)系數(shù)為0.168，間接效應(yīng)系數(shù)為-0.075，并且均至少通過了5%水平的顯著性檢驗(yàn)，表明本地城市人力資本水平的提高對(duì)其自身城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展具有推動(dòng)作用，這與吳一凡等［8］的研究結(jié)果較一致，但對(duì)其鄰近地區(qū)具有負(fù)向溢出效應(yīng)，即抑制了鄰域城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展.對(duì)于城市來說，人力資本會(huì)通過人力資本規(guī)模擴(kuò)大和人力資本外部性對(duì)本地城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展產(chǎn)生正向促進(jìn)作用，人力資本規(guī)模的擴(kuò)大會(huì)增加高素質(zhì)人才比例，引導(dǎo)人才供給與產(chǎn)業(yè)發(fā)展相匹配，促進(jìn)當(dāng)?shù)禺a(chǎn)業(yè)轉(zhuǎn)型升級(jí)，實(shí)現(xiàn)城市人口結(jié)構(gòu)和產(chǎn)業(yè)結(jié)構(gòu)的匹配協(xié)調(diào)；人力資本外部性會(huì)促進(jìn)信息、知識(shí)的外溢，加快人才、勞動(dòng)力的集聚，從而實(shí)現(xiàn)城市綜合發(fā)展水平和質(zhì)量的提升.同時(shí)，人力資本這一因素對(duì)周邊區(qū)域城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展具有較強(qiáng)的空間負(fù)外部性，原因是人才這一生產(chǎn)性要素具有較強(qiáng)的流動(dòng)性，人力資本水平高的城市會(huì)對(duì)鄰近地區(qū)的高技能、高素質(zhì)人才產(chǎn)生虹吸效應(yīng)，間接抑制鄰域地區(qū)的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展.

創(chuàng)新水平的直接效應(yīng)系數(shù)為0.006，在1%水平上顯著，間接效應(yīng)系數(shù)為-0.004，通過了10%水平的顯著性檢驗(yàn)，表明科技創(chuàng)新水平的提升會(huì)促進(jìn)城市自身城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展，與范擎宇等［10］的研究結(jié)論一致，但會(huì)通過負(fù)向溢出效應(yīng)抑制鄰域城鎮(zhèn)化耦合協(xié)調(diào)程度.科技創(chuàng)新作為人力資本集聚的產(chǎn)出和體現(xiàn)，其對(duì)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的影響機(jī)理與人力資本相似，從城市自身來分析，科技創(chuàng)新不僅可以通過推動(dòng)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)、優(yōu)化城市空間布局來提高資源和能源的利用效率，也可以通過培育環(huán)保型的現(xiàn)代產(chǎn)業(yè)企業(yè)來降低對(duì)資源和能源的消耗.總的來說，科技創(chuàng)新主要是通過優(yōu)化城市產(chǎn)業(yè)結(jié)構(gòu)和空間布局，間接提升城市生態(tài)環(huán)境質(zhì)量，進(jìn)而促進(jìn)城鎮(zhèn)化的協(xié)調(diào)發(fā)展.從鄰域角度分析，科技創(chuàng)新水平高的地區(qū)也會(huì)通過虹吸效應(yīng)對(duì)鄰域地區(qū)產(chǎn)生負(fù)外部性，使得周邊城市的城鎮(zhèn)化耦合協(xié)調(diào)程度降低.

產(chǎn)業(yè)結(jié)構(gòu)的直接效應(yīng)系數(shù)為0.105，間接效應(yīng)系數(shù)為-0.130，兩者均通過了1%水平的顯著性檢驗(yàn)，表明產(chǎn)業(yè)結(jié)構(gòu)的升級(jí)有助于城市自身城鎮(zhèn)化耦合協(xié)調(diào)程度的提升，這與周艷等［11］、張鵬巖等［28］、朱高立等［44］的研究結(jié)果較一致，但對(duì)鄰域城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展具有空間負(fù)外部性.從城市本身分析，產(chǎn)業(yè)結(jié)構(gòu)升級(jí)意味著城市自身的產(chǎn)業(yè)體系逐漸轉(zhuǎn)變?yōu)橐愿咝录夹g(shù)產(chǎn)業(yè)、先進(jìn)制造業(yè)、現(xiàn)代服務(wù)業(yè)等資本密集型和技術(shù)密集型產(chǎn)業(yè)為主，城市的經(jīng)濟(jì)發(fā)展水平和活力不斷增強(qiáng)，城鎮(zhèn)化質(zhì)量得到持續(xù)提升.同時(shí)，產(chǎn)業(yè)結(jié)構(gòu)水平的提高也會(huì)增加城市就業(yè)空間，吸引高素質(zhì)人才和高技能勞動(dòng)力不斷向城市集聚，使得城市自身的城鎮(zhèn)化水平和質(zhì)量均得到了提升并在更深層次上相協(xié)調(diào).而對(duì)于鄰近區(qū)域，隨著城市自身產(chǎn)業(yè)結(jié)構(gòu)的升級(jí)和產(chǎn)業(yè)體系的調(diào)整，城市傾向于將其資源密集型產(chǎn)業(yè)和能源依賴型企業(yè)就近轉(zhuǎn)移給鄰近城市，造成鄰域地區(qū)能源資源浪費(fèi)、生態(tài)環(huán)境質(zhì)量下降，從而抑制了周邊區(qū)域的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展.

政府能力的直接效應(yīng)系數(shù)為0.171，通過了1%水平的顯著性檢驗(yàn)，間接效應(yīng)系數(shù)為0.011，未通過顯著性檢驗(yàn)，表明政府能力的提高對(duì)本地城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展具有推動(dòng)作用，這與郭付友等［27］的研究結(jié)論一致，但其對(duì)鄰近地域的城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展并無顯著影響.對(duì)于城市自身來說，政府能力尤其是財(cái)政能力的提升，會(huì)通過增加科技、教育、醫(yī)療等的支出提高城市公共服務(wù)供給水平，進(jìn)而提升城市人居環(huán)境質(zhì)量和人民生活質(zhì)量，吸引人口等資源要素向城市流動(dòng)，同步拉動(dòng)城鎮(zhèn)化水平和質(zhì)量提升；同時(shí)，政府財(cái)政的增加也有助于增加社會(huì)固定資產(chǎn)投資，并從短期的增加生產(chǎn)要素投入、提升投資需求和長(zhǎng)期的彌補(bǔ)社會(huì)存量、擴(kuò)大社會(huì)再生產(chǎn)拉動(dòng)宏觀經(jīng)濟(jì)增長(zhǎng)［45］，實(shí)現(xiàn)城市人口集聚、建成區(qū)面積擴(kuò)大以及土地利用強(qiáng)度和效率提升，從對(duì)城鎮(zhèn)化水平和質(zhì)量的協(xié)調(diào)發(fā)展產(chǎn)生積極影響.

3.3 穩(wěn)健性分析

隨著數(shù)據(jù)科學(xué)的發(fā)展，以遙感數(shù)據(jù)為代表的地理空間大數(shù)據(jù)越來越多的被學(xué)者們所采用，成為反映和感知人類活動(dòng)的重要數(shù)據(jù)來源.其中，以美國(guó)國(guó)家海洋和大氣管理局（NOAA）發(fā)布的夜間燈光數(shù)據(jù)集，不僅能從微觀尺度刻畫人類的社會(huì)經(jīng)濟(jì)活動(dòng)，也為長(zhǎng)時(shí)間監(jiān)測(cè)區(qū)域城鎮(zhèn)化發(fā)展提供了可能.因此在參考Li和Wang［46］的研究基礎(chǔ)上，本文利用城市燈光平均亮度值和燈光面積占轄區(qū)總面積的比重作為城鎮(zhèn)化質(zhì)量和城鎮(zhèn)化水平的替代指標(biāo).利用耦合協(xié)調(diào)度模型和空間杜賓模型對(duì)本文的研究結(jié)論進(jìn)行穩(wěn)健性分析.另一方面，考慮到不同的空間權(quán)重矩陣所代表的空間相互關(guān)系有所不同，從而導(dǎo)致研究的分析結(jié)果存在差異.因此，本文進(jìn)一步將鄰接權(quán)重矩陣替換為反距離空間權(quán)重矩陣并對(duì)空間計(jì)量模型進(jìn)行估計(jì)，以驗(yàn)證不同空間權(quán)重矩陣下實(shí)證結(jié)果的一致性和保證研究結(jié)果的可靠性.表8和表9的穩(wěn)健性檢驗(yàn)計(jì)算結(jié)果顯示無論是更換被解釋變量還是替換空間權(quán)重矩陣，被解釋變量空間滯后項(xiàng)、解釋變量空間滯后項(xiàng)和空間誤差項(xiàng)系數(shù)的正負(fù)號(hào)和顯著性均與4.2的基礎(chǔ)回歸分析結(jié)果保持一致，說明本文的實(shí)證分析結(jié)果具有較強(qiáng)的穩(wěn)健性.

4 結(jié)論與討論

4.1 討論

本研究分析了2003—2019年中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的時(shí)空演化格局、時(shí)空交互特征和影響因素機(jī)制.與以往城鎮(zhèn)化協(xié)調(diào)發(fā)展的研究相比［9－10，29］，本文從城鎮(zhèn)化水平和質(zhì)量?jī)蓚€(gè)維度進(jìn)行耦合協(xié)調(diào)度分析，能較全面地揭示城鎮(zhèn)化發(fā)展的內(nèi)涵.從時(shí)空演變特征來看，城鎮(zhèn)化耦合協(xié)調(diào)度整體不斷提升，高值集聚分布在珠三角、長(zhǎng)三角和京津冀等地區(qū)，低值集聚分布在中、西部和東北地區(qū)，空間集聚特征和趨同特征明顯，與金丹等［9］研究結(jié)果相同.從城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的影響因素分析，經(jīng)濟(jì)發(fā)展、人口規(guī)模、政府能力、科技創(chuàng)新等因素對(duì)城鎮(zhèn)化協(xié)調(diào)程度的直接效應(yīng)均較顯著，與范擎宇等［10］的研究結(jié)論相似，但上述因素的間接效應(yīng)為負(fù)或不顯著，與前人的研究［10］不相一致，究其原因，現(xiàn)有研究?jī)A向于以長(zhǎng)三角等沿海發(fā)達(dá)地區(qū)為研究對(duì)象，這些區(qū)域的城鎮(zhèn)化發(fā)展已具規(guī)模，影響城鎮(zhèn)化協(xié)調(diào)發(fā)展的各項(xiàng)因素具有明顯的空間溢出效應(yīng)，而本文以全國(guó)285座城市為樣本，大部分城市尚未形成規(guī)模效應(yīng)，對(duì)周邊地區(qū)的虹吸效應(yīng)強(qiáng)于溢出效應(yīng)，這也間接說明本文研究結(jié)論的普適性和可信度較高.

本文從城鎮(zhèn)化水平和質(zhì)量?jī)蓚€(gè)方面考察城鎮(zhèn)化的協(xié)調(diào)發(fā)展水平，并將空間溢出效應(yīng)納入到研究之中，對(duì)于推動(dòng)我國(guó)城鎮(zhèn)化的協(xié)調(diào)發(fā)展具有重要的參考價(jià)值.城鎮(zhèn)化發(fā)展具有多維性，涵蓋領(lǐng)域和內(nèi)容諸多.本文雖然選取了諸多指標(biāo)對(duì)城鎮(zhèn)化水平和質(zhì)量的耦合協(xié)調(diào)發(fā)展水平進(jìn)行測(cè)度，并利用夜間燈光數(shù)據(jù)從燈光面積占比和燈光平均亮度代替城鎮(zhèn)化水平和質(zhì)量進(jìn)行穩(wěn)健性分析，但在未來研究過程中還需進(jìn)一步融合諸如POI等數(shù)據(jù)集對(duì)城鎮(zhèn)化的協(xié)調(diào)發(fā)展進(jìn)行測(cè)度，使研究指標(biāo)的選取更加豐富和多樣.此外，由于省域單元的研究較多及縣域單元數(shù)據(jù)獲取難度大，本文僅從市域單元對(duì)城鎮(zhèn)化協(xié)調(diào)發(fā)展的時(shí)空演化和空間效應(yīng)進(jìn)行分析，這使得研究結(jié)論的實(shí)踐指導(dǎo)建議在省域和縣域尺度的適用性有所降低.在不同尺度下城鎮(zhèn)化協(xié)調(diào)發(fā)展的內(nèi)在驅(qū)動(dòng)機(jī)制可能存在差異，未來的研究中將進(jìn)一步分析省域和縣域尺度的城鎮(zhèn)化協(xié)調(diào)發(fā)展，并深入探討不同尺度下城鎮(zhèn)化協(xié)調(diào)發(fā)展時(shí)空演化和形成機(jī)制的內(nèi)在差異，以便獲取更具針對(duì)性的研究結(jié)論和政策啟示.

4.2 結(jié)論

本文以2003—2019年中國(guó)285個(gè)地級(jí)及以上城市的面板數(shù)據(jù)為研究樣本，利用耦合協(xié)調(diào)度模型和ESTDA方法從城鎮(zhèn)化水平和質(zhì)量視角探討城鎮(zhèn)化協(xié)調(diào)發(fā)展的時(shí)空演化格局和交互特征，并借助空間面板杜賓模型探究城鎮(zhèn)化協(xié)調(diào)發(fā)展的影響因素及空間效應(yīng)，得出以下主要結(jié)論：

1）從時(shí)空演化格局來看，研究期間內(nèi)，中國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展水平較低，但總體呈現(xiàn)波動(dòng)上升趨勢(shì)，協(xié)調(diào)度的提升空間較大；城鎮(zhèn)化的耦合協(xié)調(diào)空間格局具有明顯的“東高西低”和“南高北低”的地帶性特征，高值區(qū)主要分布在珠三角、長(zhǎng)三角和京津冀等城市群；市域間城鎮(zhèn)化協(xié)調(diào)發(fā)展的整體差異不斷縮小，三大經(jīng)濟(jì)區(qū)由大到小依次為東部、西部和中部地區(qū).

2）從時(shí)空交互特征來看，中國(guó)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展具有正相關(guān)性，呈現(xiàn)空間集聚分布，并且集聚程度在時(shí)間趨勢(shì)上遞增；時(shí)間路徑分析表明我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展具有較強(qiáng)的空間依賴性和穩(wěn)定的空間結(jié)構(gòu)，其空間格局演變呈現(xiàn)空間整合的動(dòng)態(tài)性特征；時(shí)空躍遷分析表明我國(guó)城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的分布格局具有明顯的路徑依賴和空間鎖定特征.

3）空間面板杜賓模型的分析結(jié)果表明，經(jīng)濟(jì)發(fā)展、人口規(guī)模、人力資本、科技創(chuàng)新、產(chǎn)業(yè)結(jié)構(gòu)、政府能力對(duì)本地城市的城鎮(zhèn)化協(xié)調(diào)發(fā)展均具有促進(jìn)作用，其中人力資本、產(chǎn)業(yè)結(jié)構(gòu)及政府能力是影響城鎮(zhèn)化協(xié)調(diào)發(fā)展程度的關(guān)鍵因素；人口規(guī)模對(duì)鄰域城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間溢出效應(yīng)為正，人力資本、科技創(chuàng)新和產(chǎn)業(yè)結(jié)構(gòu)的空間溢出效應(yīng)為負(fù)，而經(jīng)濟(jì)發(fā)展水平和政府能力的空間效應(yīng)不顯著.

根據(jù)以上的研究結(jié)論，結(jié)合新型城鎮(zhèn)化發(fā)展戰(zhàn)略規(guī)劃實(shí)施以來城市建設(shè)過程中存在的問題，本文提出如下政策建議.

針對(duì)當(dāng)前我國(guó)城鎮(zhèn)化發(fā)展“重形式輕內(nèi)涵”而導(dǎo)致的城鎮(zhèn)化水平和質(zhì)量不相協(xié)調(diào)的現(xiàn)狀，各地應(yīng)嚴(yán)格按照國(guó)家新型城鎮(zhèn)化發(fā)展戰(zhàn)略規(guī)劃目標(biāo)，積極推進(jìn)以人為核心的城鎮(zhèn)化發(fā)展，提高城市治理能力和管理服務(wù)水平，著力解決人口城鎮(zhèn)化滯后于土地城鎮(zhèn)化的問題，從而推動(dòng)我國(guó)城鎮(zhèn)化發(fā)展質(zhì)量與快速擴(kuò)張的城鎮(zhèn)化水平相適應(yīng).

考慮到中國(guó)的城鎮(zhèn)化建設(shè)具有明顯的區(qū)域異質(zhì)性特征，在未來的城市發(fā)展過程中，應(yīng)因地制宜地制定推動(dòng)城鎮(zhèn)化協(xié)調(diào)發(fā)展的指導(dǎo)意見和政策方針.對(duì)于經(jīng)濟(jì)發(fā)展水平較高的東部先發(fā)地區(qū)，要根據(jù)城市主體功能定位，大力發(fā)展戰(zhàn)略性新興產(chǎn)業(yè)和現(xiàn)代服務(wù)業(yè)，并通過優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)、增強(qiáng)創(chuàng)新能力、提高經(jīng)濟(jì)效益、降低能源消耗等方式推動(dòng)城鎮(zhèn)化發(fā)展階段由以規(guī)模擴(kuò)展為主轉(zhuǎn)變?yōu)橐再|(zhì)量提升為主.

實(shí)證分析結(jié)果表明經(jīng)濟(jì)發(fā)展等因素能顯著促進(jìn)城鎮(zhèn)化的耦合協(xié)調(diào)發(fā)展，因此在城鎮(zhèn)化建設(shè)過程中，一是要加快經(jīng)濟(jì)發(fā)展速度、提高經(jīng)濟(jì)發(fā)展質(zhì)量，通過增加基礎(chǔ)設(shè)施投資和公共服務(wù)投資，吸引人才等優(yōu)質(zhì)資源向城鎮(zhèn)集聚；二是要合理控制城市人口規(guī)模，既避免城市規(guī)模較小造成的城市規(guī)模效應(yīng)不足，又要避免城市規(guī)模過大導(dǎo)致的城市無序蔓延發(fā)展，從而提高城市土地利用效率和城市發(fā)展質(zhì)量；三是要依托城市人才、科技等資源優(yōu)勢(shì)，增強(qiáng)城市創(chuàng)新能力，提高城市資源利用效率，從而發(fā)揮創(chuàng)新對(duì)城鎮(zhèn)化發(fā)展的引領(lǐng)作用；四是要調(diào)整城市產(chǎn)業(yè)結(jié)構(gòu)和布局，依據(jù)城市資源稟賦和比較優(yōu)勢(shì)，培育和形成以高新技術(shù)產(chǎn)業(yè)、先進(jìn)制造業(yè)等為主的現(xiàn)代化產(chǎn)業(yè)體系；五是要提高政府治理體系和治理能力現(xiàn)代化水平，完善城市治理結(jié)構(gòu)，創(chuàng)新城市治理模式.

要從交通、市場(chǎng)等方面加快區(qū)域一體化建設(shè)，強(qiáng)化城市間勞動(dòng)力、資本、產(chǎn)業(yè)等各項(xiàng)資源要素的聯(lián)系程度，從而加強(qiáng)區(qū)域城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展的空間聯(lián)動(dòng)效應(yīng).要通過優(yōu)化城市空間結(jié)構(gòu)體系來增強(qiáng)直轄市等中心城市的輻射帶動(dòng)作用，從而發(fā)揮中心城市對(duì)我國(guó)城鎮(zhèn)化發(fā)展的支撐和引領(lǐng)作用.最后，從人才、科創(chuàng)、產(chǎn)業(yè)等方面構(gòu)建區(qū)域協(xié)調(diào)發(fā)展和利益補(bǔ)償機(jī)制，降低中心城市對(duì)中小城市的“虹吸效應(yīng)”，推動(dòng)區(qū)域城鎮(zhèn)化的協(xié)同發(fā)展.

參考文獻(xiàn)：

［1］王怡睿，黃煌，石培基.中國(guó)城鎮(zhèn)化質(zhì)量時(shí)空演變研究［J］.經(jīng)濟(jì)地理，2017，37（1）：90.

［2］范擎宇，楊山，胡信.耦合視角下長(zhǎng)三角地區(qū)城鎮(zhèn)化協(xié)調(diào)度的時(shí)空特征及交互機(jī)制［J］.地理研究，2020，39（2）：289.

［3］賴行健.我國(guó)城鎮(zhèn)化質(zhì)量與水平協(xié)調(diào)發(fā)展關(guān)系研究［J］.開發(fā)研究，2017（5）：34.

［4］王建康，谷國(guó)鋒，姚麗，等.中國(guó)新型城鎮(zhèn)化的空間格局演變及影響因素分析——基于285個(gè)地級(jí)市的面板數(shù)據(jù)［J］.地理科學(xué)，2016，36（1）：63.

［5］方創(chuàng)琳.中國(guó)新型城鎮(zhèn)化高質(zhì)量發(fā)展的規(guī)律性與重點(diǎn)方向［J］.地理研究，2019，38（1）：13.

［6］張榮天，焦華富.中國(guó)新型城鎮(zhèn)化研究綜述與展望［J］.世界地理研究，2016，25（1）：59.

［7］于猛，宋家寧.我國(guó)城鎮(zhèn)化模式研究綜述——基于規(guī)模選擇和典型模式推廣的分析討論［J］.中國(guó)土地，2013（3）：54.

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

［9］金丹，戴林琳.中國(guó)人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)發(fā)展的時(shí)空特征與驅(qū)動(dòng)因素［J］.中國(guó)土地科學(xué)，2021，35（6）：74.

［10］范擎宇，楊山.長(zhǎng)三角地區(qū)城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間特征及形成機(jī)理［J］.地理科學(xué)進(jìn)展，2021，40（1）：124.

［11］周艷，黃賢金，徐國(guó)良，等.長(zhǎng)三角城市土地?cái)U(kuò)張與人口增長(zhǎng)耦合態(tài)勢(shì)及其驅(qū)動(dòng)機(jī)制［J］.地理研究，2016，35（2）：313.

［12］徐靜嫻.西部地區(qū)人口城市化與土地城市化耦合協(xié)調(diào)發(fā)展研究［J］.湖北農(nóng)業(yè)科學(xué)，2021，60（17）：139.

［13］劉娟，鄭欽玉，郭銳利，等.重慶市人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)發(fā)展評(píng)價(jià)［J］.西南師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，37（11）：66.

［14］張春梅，張小林，吳啟焰，等.城鎮(zhèn)化質(zhì)量與城鎮(zhèn)化規(guī)模的協(xié)調(diào)性研究——以江蘇省為例［J］.地理科學(xué)，2013，33（1）：16.

［15］陶然，曹廣忠.“空間城鎮(zhèn)化”、“人口城鎮(zhèn)化”的不匹配與政策組合應(yīng)對(duì)［J］.改革，2008（10）：83.

［16］范進(jìn)，趙定濤.土地城鎮(zhèn)化與人口城鎮(zhèn)化協(xié)調(diào)性測(cè)定及其影響因素［J］.經(jīng)濟(jì)學(xué)家，2012（5）：61.

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

［18］唐鵬，周來友，石曉平.地方政府對(duì)土地財(cái)政依賴的影響因素研究——基于中國(guó)1998—2010年的省際面板數(shù)據(jù)分析［J］.資源科學(xué)，2014，36（7）：1374.

［19］LIU Y，ZHOU G，LIU D，et al.The interaction of population，industry and land in process of urbanization in China：a case study in Jilin Province［J］.Chinese Geographical Science，2018，28（3）：529.

［20］SUN Y，CUI Y，GIL－LAFENTE A M.Analyzing the coupling coordination among economic，social，and environmental benefits of urban infrastructure：case study of four chinese autonomous municipalities［J］.Mathematical Problems in Engineering，2018，2018：828038.

［21］李曉慶，姜博，初楠臣.城市群土地集約利用與新型城鎮(zhèn)化耦合協(xié)調(diào)分析——以我國(guó)三大城市群為例［J］.長(zhǎng)江流域資源與環(huán)境，2017，26（7）：983.

［22］趙建吉，劉巖，朱亞坤，等.黃河流域新型城鎮(zhèn)化與生態(tài)環(huán)境耦合的時(shí)空格局及影響因素［J］.資源科學(xué)，2020，42（1）：159.

［23］王玉娟，江成濤，蔣長(zhǎng)流.新型城鎮(zhèn)化與低碳發(fā)展能夠協(xié)調(diào)推進(jìn)嗎？——基于284個(gè)地級(jí)及以上城市的實(shí)證研究［J］.財(cái)貿(mào)研究，2021，32（9）：32.

［24］李建豹，黃賢金，揣小偉，等.江蘇省人口城鎮(zhèn)化與能源消費(fèi)CO2排放耦合協(xié)調(diào)度時(shí)空格局及影響因素［J］.經(jīng)濟(jì)地理，2021，41（5）：57.

［25］李江蘇，王曉蕊，苗長(zhǎng)虹，等.城鎮(zhèn)化水平與城鎮(zhèn)化質(zhì)量協(xié)調(diào)度分析——以河南省為例［J］.經(jīng)濟(jì)地理，2014，34（10）：70.

［26］呂添貴，吳次芳，游和遠(yuǎn).鄱陽湖生態(tài)經(jīng)濟(jì)區(qū)水土資源與經(jīng)濟(jì)發(fā)展耦合分析及優(yōu)化路徑［J］.中國(guó)土地科學(xué)，2013，27（9）：3.

［27］郭付友，李誠(chéng)固，陳才，等.2003年以來東北地區(qū)人口城鎮(zhèn)化與土地城鎮(zhèn)化時(shí)空耦合特征［J］.經(jīng)濟(jì)地理，2015，35（9）：49.

［28］張鵬巖，楊丹，李二玲，等.人口城鎮(zhèn)化與土地城鎮(zhèn)化的耦合協(xié)調(diào)關(guān)系——以中原經(jīng)濟(jì)區(qū)為例［J］.經(jīng)濟(jì)地理，2017，37（8）：145.

［29］呂添貴，吳次芳，李洪義，等.人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)性測(cè)度及優(yōu)化——以南昌市為例［J］.地理科學(xué)，2016，36（2）：239.

［30］ZHANG Y，PAN J，ZHANG Y，et al.Spatial－temporal characteristics and decoupling effects of Chinas carbon footprint based on multi－source data［J］.Journal of Geographical Sciences，2021，31（3）：327.

［31］張欣，潘竟虎.中國(guó)城市蔓延時(shí)空動(dòng)態(tài)識(shí)別及驅(qū)動(dòng)因素探測(cè)［J］.人文地理，2021，36（4）：114.

［32］許靜，楊永春.“一帶一路”背景下中國(guó)西部城市經(jīng)濟(jì)全球化水平的空間演化［J］.經(jīng)濟(jì)地理，2021，41（7）：20.

［33］張榮天，張小林，陸建飛，等.我國(guó)鄉(xiāng)村轉(zhuǎn)型發(fā)展時(shí)空分異格局與影響機(jī)制分析［J］.人文地理，2021，36（3）：138.

［34］REY S J，MURRAY A T，ANSELIN L.Visualizing regional income distribution dynamics［J］.Letters in Spatial and Resource Sciences，2011，4（1）：81.

［35］蘇聰文，鄧宗兵，李莉萍，等.中國(guó)水生態(tài)文明發(fā)展水平的空間格局及收斂性［J］.自然資源學(xué)報(bào)，2021，36（5）：1282.

［36］馬慧強(qiáng)，廉倩文，韓增林，等.基本公共服務(wù)—城鎮(zhèn)化—區(qū)域經(jīng)濟(jì)耦合協(xié)調(diào)發(fā)展時(shí)空演化［J］.經(jīng)濟(jì)地理，2020，40（5）：19.

［37］尹鵬，王富喜，段佩利.中國(guó)基本公共服務(wù)效率與城鎮(zhèn)化質(zhì)量的時(shí)空耦合關(guān)系研究［J］.地理科學(xué)，2021，41（4）：571.

［38］CHEN Z，YU B，YANG C，et al.An extended time－series（2000—2018） of global NPP－VIIRS－like nighttime light data［J］.Earth System Science Data，2021，13（3）：889.

［39］SMALL C，ELVIDGE C D.Night on earth：Mapping decadal changes of anthropogenic night light in Asia［J］.International Journal of Applied Earth Observation & Geoinformation，2013，22：40.

［40］邵帥，李欣，曹建華，等.中國(guó)霧霾污染治理的經(jīng)濟(jì)政策選擇——基于空間溢出效應(yīng)的視角［J］.經(jīng)濟(jì)研究，2016，51（9）：73.

［41］朱高立，鄒偉，王雪琪.經(jīng)濟(jì)結(jié)構(gòu)調(diào)整對(duì)人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)性的影響差異［J］.中國(guó)人口·資源與環(huán)境，2018，28（5）：93.

［42］宋麗智.我國(guó)固定資產(chǎn)投資與經(jīng)濟(jì)增長(zhǎng)關(guān)系再檢驗(yàn)：1980—2010年［J］.宏觀經(jīng)濟(jì)研究，2011（11）：17.

［43］LI X，WANG L.Does administrative division Adjustment promote low－carbon city development？Empirical evidence from the“revoke county to urban district”in China［J］.Environmental Science and Pollution Research，2023，30（5）：11542.

（責(zé)任編輯 武維寧）

收稿日期：2023－11－21;修改稿收到日期：2024－03－06

基金項(xiàng)目：國(guó)家社科基金重大項(xiàng)目（19YJA630079）；教育部人文社科研究規(guī)劃基金項(xiàng)目（19YJA630079）

作者簡(jiǎn)介：栗向陽（1995—），河南駐馬店人，助理研究員，博士.主要研究方向?yàn)槌鞘泻蛥^(qū)域發(fā)展與治理.E－mail：lixiangyang@whu.edu.cn

*通信聯(lián)系人，講師，博士.主要研究方向?yàn)槿颦h(huán)境變遷.E－mail：zhangyx@xynu.edu.cn