梁彥慶，官洋溢，黃志英，焦新穎 (.河北師范大學資源與環(huán)境科學學院，河北 石家莊 00024；2.河北省環(huán)境變化遙感識別技術創(chuàng)新中心，河北 石家莊 00024；3.中國科學院新疆生態(tài)與地理研究所荒漠與綠洲生態(tài)國家重點實驗室，新疆 烏魯木齊 8300；4.新疆大學資源與環(huán)境科學學院，新疆 烏魯木齊 830046；.河北地質大學土地科學與空間規(guī)劃學院，河北 石家莊 0003；6.中國科學院大學，北京 00049)

中國經濟的高速發(fā)展給土地生態(tài)環(huán)境造成了巨大壓力，土地資源綠色利用與經濟增長的合理匹配成為實現(xiàn)國家可持續(xù)發(fā)展的重要前提。隨著城市化進程的不斷加快，由土地利用引發(fā)的資源環(huán)境保護與經濟發(fā)展之間的矛盾沖突愈加嚴峻，這與我國當前綠色、協(xié)調的發(fā)展理念不符。土地綠色利用與生產已成為生態(tài)文明建設新時期土地利用的新模式、新需求，將綠色發(fā)展理念貫徹到土地利用過程中，在土地總量有限、生態(tài)環(huán)境脆弱的前提下，通過土地綠色利用轉型減少資源投入與污染排放，實現(xiàn)土地利用中經濟、社會和環(huán)境綜合效益的最優(yōu)協(xié)調發(fā)展，成為繼土地集約利用之后的一種土地利用新理念[1-3]。

學者們圍繞土地利用與城鎮(zhèn)化、區(qū)域一體化的耦合關系進行了相關研究，大致分為兩類：一是全國、區(qū)域、城市等不同尺度下的研究。如張紅鳳等[4]從時點狀態(tài)和時段變化的角度對我國各省份的城鎮(zhèn)化發(fā)展和土地集約利用水平的耦合特征進行了分析；陳萬旭等[5]采用空間計量經濟學模型探測長江中游土地利用轉型與城鎮(zhèn)化之間的空間關系；邢穎等[6]對都勻市土地利用變化與經濟增長的協(xié)調性進行了考察。二是時空效應、協(xié)調發(fā)展程度評價等不同視角下的研究。如李曉慶等[7]基于耦合協(xié)調發(fā)展度及空間自相關模型探究長江兩大城市群新型城鎮(zhèn)化與土地集約利用的協(xié)調時空分異特征；苗苗等[8]探討長江中游城市群土地利用與社會-經濟-自然系統(tǒng)的內在耦合協(xié)調性；張志等[9]從空間尺度分析了湖北省各城市土地集約利用程度與經濟發(fā)展的耦合協(xié)調關系。上述研究大多注重解決經濟、社會或環(huán)境某一方面的問題，著眼于全面考慮資源節(jié)約、環(huán)境保護、生態(tài)可持續(xù)的土地綠色利用研究較少，且土地綠色利用與經濟增長的相關研究也局限于某個因素的單項研究，缺乏對兩者的耦合分析，特別是對耦合協(xié)調時空異質性及其驅動因素的深入探討較少。

京津冀城市群作為我國北方人居環(huán)境的重要保障地帶，土地生態(tài)保護與經濟發(fā)展等矛盾突出，基于2006、2011和2016年3期截面數(shù)據，探索土地綠色利用與經濟增長耦合協(xié)調發(fā)展的時空分異規(guī)律及其影響因素，可為實現(xiàn)京津冀區(qū)域土地資源高效綠色利用與經濟社會可持續(xù)發(fā)展提供借鑒和參考。

1 研究區(qū)域概況與數(shù)據來源

1.1 研究區(qū)域概況

京津冀地區(qū)是我國北方最重要的經濟中心，地處36°01′～42°37′ N，113°04′～119°53′ E，總面積約為21.80萬km2，包括北京、天津和河北3個省級行政區(qū)，轄13個城市，擁有良好的區(qū)位條件和資源稟賦，現(xiàn)已成為繼長三角、珠三角之后的中國第三大經濟增長極地區(qū)。2016年，總人口超過1億人，地區(qū)生產總值達到74 612.6億元，占全國的10%。然而，由于京津冀不同發(fā)展定位城市的土地利用模式、結構不同，造成區(qū)域間經濟發(fā)展不平衡、土地利用效益低下、生態(tài)環(huán)境污染等問題一直比較突出。為保證京津冀協(xié)同可持續(xù)發(fā)展，亟需重點研究和解決土地資源高效綠色利用與經濟發(fā)展的協(xié)調問題。

1.2 數(shù)據來源

數(shù)據來源于2007、2012和2017年的《北京統(tǒng)計年鑒》《天津統(tǒng)計年鑒》《河北統(tǒng)計年鑒》《中國城市統(tǒng)計年鑒》和各市統(tǒng)計公報，無法直接獲取的各類地均資源投入和經濟密度等指標值根據相關指標計算得到。

2 研究方法與指標體系

2.1 研究方法

2.1.1超效率SBM模型

土地綠色利用水平可采用考慮資源與環(huán)境代價的綜合土地利用效率來衡量?；诜瞧谕a出的SBM模型于2001年由TONE[10]提出，能夠解決傳統(tǒng)數(shù)據包絡模型造成的投入產出變量松弛性問題，但難以進一步區(qū)分有效決策單元間的效率差異，而超效率DEA模型可實現(xiàn)對所有效率值大于1的決策單元(DMU)的區(qū)分，因此運用結合2種模型優(yōu)勢的超效率SBM模型對綜合土地利用效率進行測算[10]，不僅可以恰當處理非期望產出，同時能夠解決非有效單元和有效單元的排序問題。模型公式為

(1)

(2)

式(1)～(2)中，ρ為測度的土地綠色利用水平值；假設有n個DMU，每個DMU由投入m、期望產出r1和非期望產出r2構成，其投入、期望產出和非期望產出向量分別為x∈Rm、yd∈Rr1、yu∈Rr2；λj為約束條件。

2.1.2標準差橢圓法

標準差橢圓法最早由Lefever提出，通過重心、長半軸、短半軸、方位角等構成要素能夠揭示土地綠色利用與經濟增長的空間分布特征。其中，重心表示區(qū)域要素空間分布的相對位置，長短軸之比用來描述要素空間分布形態(tài)，方位角表示空間分布的主趨勢方向[11]。平均中心坐標、x軸標準差、y軸標準差計算公式分別為

(3)

(4)

(5)

(6)

2.1.3耦合協(xié)調度模型

耦合協(xié)調度模型用來描述系統(tǒng)之間相互作用、相互影響的關聯(lián)性[12-13]，可將土地綠色利用和經濟增長2個系統(tǒng)相互作用、協(xié)調發(fā)展的程度定義為土地綠色利用-經濟增長系統(tǒng)耦合協(xié)調度。表達式為

(7)

(8)

(9)

T=αU1+βU2，

(10)

(11)

式(7)～(11)中，U1和U2分別為土地綠色利用系統(tǒng)和經濟增長系統(tǒng)的綜合得分值；pi和qi為對應指標權重；xi和yi分別為各指標對應的標準化值；m和n分別為2個系統(tǒng)指標總數(shù)；C為系統(tǒng)之間的耦合度；T為系統(tǒng)之間的綜合調和指數(shù)；α、β為待定系數(shù)，反映2個子系統(tǒng)對總系統(tǒng)貢獻的重要程度，參照文獻[14]，取α=β=1/2。

采用均勻分布函數(shù)法[15-16]，將京津冀土地綠色利用-經濟增長耦合協(xié)調度(D)劃分為如下等級：0～0.19、>0.19～0.29、>0.29～0.39、>0.39～0.49分別表示嚴重、中度、輕度與瀕臨失調；>0.49～0.59、>0.59～0.69、>0.69～0.79、>0.79～0.89、>0.89～1.00分別表示勉強、初級、中級、良好與優(yōu)質協(xié)調。

2.1.4地理探測器

地理探測器是識別地理事物空間分異性及其驅動力的統(tǒng)計學方法[17]，用于檢驗影響因素與地理事物的變化在空間上是否具有顯著一致性。其中因子探測器用來探測京津冀土地綠色利用與經濟增長耦合協(xié)調的空間分異性以及影響因素對其的解釋力，大小用q值來衡量，表達式為

(12)

式(12)中，h為變量Y或者X的分層；Nh和N分別為h層和全區(qū)的單元數(shù)；σh2和σ2分別為h層和全區(qū)Y值的方差。q的取值在[0,1]之間，數(shù)值越大，說明影響因子對土地綠色利用與經濟增長耦合協(xié)調的影響越大。

交互探測器(P)用來識別影響因子兩兩共同作用時是否會增加或減弱對土地綠色利用與經濟增長耦合協(xié)調的解釋力。P(x1∩x2)max[P(x1),P(x2)]，代表因子x1與x2交互后雙線性加強；P(x1∩x2)>P(x1)+P(x2)，代表因子x1與x2交互后非線性加強；P(x1∩x2)=P(x1)+P(x2)，代表因子x1與x2相互獨立。

2.2 指標體系構建

土地綠色利用是在一定約束條件下實現(xiàn)多目標并同時達到最優(yōu)水平的土地利用方式組合，不僅需要增加土地資源的產出強度，還要考慮在節(jié)能減排約束條件下，降低污染物排放水平，從而實現(xiàn)土地資源的綠色可持續(xù)利用，因此在指標選擇上立足于土地利用過程中生態(tài)水平的提升能力，重點關注減少資源消耗和污染排放2個方面[18]。經濟增長是實現(xiàn)經濟可持續(xù)發(fā)展的關鍵，大部分學者認為從經濟增長的有效性、穩(wěn)定性、持續(xù)性、創(chuàng)新性、分享性等方面能較好地反映經濟增長內涵，因此在指標選擇上主要關注經濟增長結構性、穩(wěn)定性、福利性等方面[19]?；谕恋鼐G色利用和經濟增長的內涵，對黃鑫等[20]、梁流濤等[2]學者研究中使用頻率較高的指標進行統(tǒng)計，并結合京津冀發(fā)展實際以及數(shù)據可獲取性、連續(xù)性等因素，構建京津冀土地綠色利用與經濟增長評價指標體系(表1)。其中土地綠色利用由資源投入、污染排放(非期望產出)、經濟產出(期望產出)3個一級指標，8個二級指標構成；經濟增長由經濟發(fā)展、國內外貿易和經濟增長效益3個一級指標，8個二級指標構成。為減少人為因素對權重的影響，土地綠色利用和經濟增長各指標權重采用熵值法確定，權重越小表明評價指標變異系數(shù)越小，提供的有效信息越少(表1)。

表1 京津冀土地綠色利用與經濟增長評價指標體系

3 結果與分析

3.1 土地綠色利用與經濟增長的時空分布特征

運用DEA-Solver Pro5軟件，選擇規(guī)模報酬不變情況下考慮非期望產出的超效率SBM模型，對京津冀13市2006—2016年的土地綠色利用水平值進行測算(圖1)。結果顯示，2006—2016年，京津冀土地綠色利用水平均值從高到低依次為北京、承德、滄州、保定、張家口、天津、唐山、石家莊、邯鄲、衡水、廊坊、邢臺和秦皇島，大致呈西高東低、北高南低的分布格局。不同城市間的土地綠色利用水平差別較大，北京、承德、滄州、保定的土地綠色利用水平均值分別為2.67、1.44、1.48、1.09，高于城市群0.91的平均水平，表明這些城市土地綠色利用水平較高效，其中北京作為京津冀城市群的核心城市，一直強調低消耗、低污染、高效益的集約型經濟增長模式，注重對土地的健康有序利用；承德、滄州、保定近年加強工業(yè)轉型和產業(yè)結構調整，發(fā)展環(huán)境友好型工業(yè)，對提高土地綠色利用水平均起到了積極作用；其余城市由于產業(yè)耗能相對較高，工業(yè)污染物排放量相對較大，土地綠色利用水平均值普遍偏低，僅張家口略高于城市群平均水平，達到0.95，秦皇島最低，僅為0.28。

圖1 2006—2016年土地綠色利用水平平均值變化

運用多因素綜合評價法測算得到京津冀13市2006—2016年的經濟增長水平值，根據其值對經濟增長動態(tài)變化進行分析。由圖2可知，2006—2016年京津冀經濟增長水平均值存在顯著的區(qū)域非均衡發(fā)展特征，在空間上呈現(xiàn)東高西低的分布格局。其中北京、天津的經濟增長水平均處于區(qū)域領先地位，經濟增長水平均值分別為0.90和0.64，高于區(qū)域平均水平；唐山、石家莊、秦皇島、廊坊等地次之；而衡水為生態(tài)休閑旅游型城市，均值最低，僅為0.03?？傮w來看，除天津外，其余各市經濟增長綜合水平總體呈下降趨勢。分階段來看，2006—2011年呈下降趨勢的有北京、石家莊、唐山、廊坊、邢臺和衡水6市，其余7市均呈上升趨勢；2011—2016年呈下降趨勢的有承德、張家口、邢臺、唐山、保定、邯鄲，其余7市呈上升趨勢。

圖2 2006—2016年經濟增長水平變化

土地綠色利用與經濟增長標準差橢圓的參數(shù)變化見表2，空間分布見圖3。

表2 土地綠色利用與經濟增長的標準差橢圓參數(shù)變化

2006—2016年京津冀地區(qū)土地綠色利用標準差橢圓整體呈北(偏東)—南(偏西)的分布格局，但在2個時間段內分別呈現(xiàn)向東北、西北方向的移動路徑。2006—2011年長短軸之比由1.92增加到2.00，表明北(偏東)—南(偏西)方向呈擴張趨勢，2011—2016年長短軸之比由2.00減小到1.86，表明北(偏西)—南(偏東)方向呈擴張趨勢，且擴張速度快于2006—2011年，區(qū)域間差距逐漸拉大，方位角轉動方向由順時針轉為逆時針，轉動較明顯，表明城市群西北部地區(qū)對生態(tài)環(huán)境的重視程度加強，土地綠色利用水平得到逐步提高。

2006—2016年京津冀13市經濟增長的空間分布整體呈北(偏東)—南(偏西)方向，移動路徑范圍較小，2006—2011年經濟增長重心位置由39.48° N，116.61° E向東北方向移至39.52° N，116.71° E，2011—2016年向北移動至39.54° N，116.71° E。近10 a長短軸之比由1.75減小到1.54，表明東北—西南方向呈持續(xù)收縮趨勢，方位角由36.55°逐漸變?yōu)?9.33°，呈順時針方向旋轉，表明城市群東北部地區(qū)的經濟發(fā)展水平較高，發(fā)展速度較快，東北部與西南部地區(qū)的經濟增長水平差距明顯拉大。

綜上所述，2006—2016年，京津冀各城市土地綠色利用與經濟增長時序變化差異顯著，土地綠色利用水平較高區(qū)域集中在北京、承德等中、北部，隨時間變化其主方向由北(偏東)—南(偏西)方向變?yōu)楸?偏西)—南(偏東)方向，而城市群西北部地區(qū)呈現(xiàn)“追趕”趨勢，土地利用效率逐步提高；經濟增長水平較高區(qū)域集中分布在北京、唐山等中、東部區(qū)域，其主方向始終為北(偏東)—南(偏西)方向，其余地區(qū)發(fā)展相對較弱。在空間演變上，兩系統(tǒng)標準差橢圓總體呈現(xiàn)東北—西南格局分布，京津冀東北部地區(qū)的土地綠色利用與經濟增長態(tài)勢好于西南地區(qū)。

3.2 土地綠色利用與經濟增長耦合協(xié)調的空間異質性

為避免數(shù)據量綱的影響，先對土地綠色利用和經濟增長的原始數(shù)據進行標準化處理，再進行耦合協(xié)調度測算，根據測算結果對京津冀近年來各市的土地綠色利用與經濟增長耦合協(xié)調水平進行等級劃分，共分為中度失調、輕度失調、瀕臨失調、勉強協(xié)調和初級協(xié)調5個等級(圖4)。

2006—2016年間，京津冀各城市土地綠色利用與經濟增長耦合協(xié)調度呈現(xiàn)緩慢上升趨勢，耦合協(xié)調度大致在0.20～0.65之間，只有少數(shù)城市在個別年份高于0.60。京津冀整體協(xié)調關系有一定改善，從0.34上升到0.41，由輕度失調轉變?yōu)闉l臨失調。

圖3 土地綠色利用與經濟增長標準差橢圓分布Fig.3 Elliptical distribution of standard deviation of green land use and economic growth

圖4 土地綠色利用與經濟增長耦合協(xié)調度的時空演變格局Fig.4 Spatial-temporal evolution patterns of coupling coordination degree between green land use and economic growth

對各個階段展開具體分析：(1)中度失調階段，2006年有張家口、承德、衡水、邢臺和邯鄲；2011年變?yōu)?個，邯鄲、張家口由于系統(tǒng)關聯(lián)度增強，變?yōu)檩p度失調，新增保定；2016年僅有3個，分別為承德、衡水和邢臺。(2)輕度失調，2006年處于該階段的有6個，分別為秦皇島、唐山、廊坊、石家莊、保定和滄州；2011年新增張家口和邯鄲，表明張家口和邯鄲的協(xié)調度較2006年有所提升，保定降為中度失調階段；2016年減少為5個，分別為張家口、保定、滄州、邯鄲和秦皇島。(3)瀕臨失調階段，2006年僅有天津；2011年無城市處于此階段；2016年有3個，分別為唐山、廊坊和石家莊，表明其土地綠色利用與經濟增長之間的關聯(lián)度增強，這些城市在發(fā)展經濟的同時也注重土地利用綠色化和對資源、環(huán)境的保護。(4)勉強協(xié)調階段，2006年未有城市處于該階段，2011和2016年處于此階段的均僅有天津。(5)初級協(xié)調階段。2006—2016年處于該階段的均僅有北京。

從圖5可知，2006—2016年間，京津冀13個城市土地綠色利用與經濟增長的耦合協(xié)調度均得到不同程度提升，但呈現(xiàn)區(qū)域發(fā)展不平衡的明顯特征，耦合協(xié)調度由高到低依次為北京、天津、石家莊、唐山、廊坊、秦皇島、滄州、保定、張家口、邯鄲、衡水、邢臺和承德。北京、天津土地綠色利用與經濟增長的耦合協(xié)調度遙遙領先，與河北省11個城市拉開差距；而從河北省內部看，石家莊位居首位，唐山、廊坊等環(huán)京津、沿海區(qū)域次之，冀西北和冀南區(qū)域城市的耦合協(xié)調度較差。通過比較土地綠色利用與經濟增長系統(tǒng)的綜合得分值可以看出，北京、天津的土地綠色利用落后于經濟增長，表現(xiàn)為土地綠色利用滯后型狀態(tài)，而河北省各市相反，表現(xiàn)為經濟增長滯后型狀態(tài)，表明京津冀城市群尚未達到土地綠色利用與經濟增長高協(xié)調、高質量的發(fā)展階段，其中北京、天津的土地利用仍較為粗放，需要前瞻性的國土空間規(guī)劃進行指導，而河北省各市的經濟發(fā)展水平較差，還需加大力度優(yōu)化產業(yè)結構，實現(xiàn)產業(yè)綠色化高端化發(fā)展。

圖5 土地綠色利用與經濟增長耦合協(xié)調度演進態(tài)勢

總體來看，近10 a來京津冀13個城市耦合協(xié)調發(fā)展水平不斷提升，形成以北京、天津為中心的輻射擴散趨勢，表明京津冀在經濟發(fā)展的同時，注重土地集約利用，提高資源利用效率和減少環(huán)境污染，逐步走上健康的可持續(xù)發(fā)展道路，但需要注意的是冀西北和冀南等區(qū)域耦合協(xié)調度處于輕度失調階段的城市占比依然很大。對于不同類型城市，未來應有所側重地提高土地綠色利用水平并優(yōu)化調整經濟結構，有效促進不同地區(qū)經濟的健康增長，從而促進土地綠色利用與經濟增長的全面協(xié)調發(fā)展。

3.3 土地綠色利用與經濟增長耦合協(xié)調的影響因素

選取與當前實際發(fā)展狀況最為相近的2016年為研究時點，運用地理探測器軟件中的因子探測器與交互探測器來探測各驅動因子對京津冀土地綠色利用與經濟增長的影響程度。借鑒文獻[21-22]并結合數(shù)據可得性和典型性，選取單位GDP能耗(X1)、人口密度(X2)、第三產業(yè)占GDP比重(X3)、城鎮(zhèn)化率(X4)、建成區(qū)綠化率(X5)、專利申請授權量(X6)和工業(yè)固體廢物綜合利用率(X7)共7個指標作為土地綠色利用與經濟增長協(xié)調發(fā)展的驅動因子，計算得到各指標的解釋力大小q分別為0.486 6、0.399 8、0.552 2、0.799 4、0.340 3、0.614 8、0.229 6?？梢?，京津冀土地綠色利用與經濟增長協(xié)調的主導驅動因子及其解釋力差異較為明顯。

整體來看，因子探測分析結果中所有驅動因子的解釋力q值均處于顯著范圍內，但其大小具有明顯差距。其中，城鎮(zhèn)化率q值最大，明顯區(qū)別于其他因子，表明該因子對土地綠色利用-經濟增長協(xié)調發(fā)展的解釋力具有主導作用，因此屬于主導驅動因子；其次為專利申請授權量、第三產業(yè)占GDP比重，其q值均高達0.5以上，這些因子涉及科技創(chuàng)新、產業(yè)結構等方面，能較為全面地反映城市的可持續(xù)和健康發(fā)展狀況，對于土地綠色利用與經濟增長的協(xié)調發(fā)展具有重要促進作用，屬于高水平影響因子；而單位GDP能耗、人口密度、建成區(qū)綠化率和工業(yè)固體廢物綜合利用率的影響程度屬于較低水平，q值均處于0.5以下，由于不同城市在經濟發(fā)展中各有側重，對綠化率、工業(yè)固體廢物綜合利用率的要求差異較大，導致對系統(tǒng)協(xié)調發(fā)展的綜合影響力較弱；能耗、人口的增加會導致土地、水電等各種資源需求量和污染物排放增加等問題，這些因子對土地綠色利用與經濟增長協(xié)調發(fā)展的影響力均較弱，屬于一般影響因子。

從表3可知，各驅動因子間存在顯著交互作用，且任意兩因子交互作用后的作用力顯著高于單因子，具有雙因子增強或非線性增強特征，表明產業(yè)結構、生態(tài)投入、科技創(chuàng)新、城市化水平等方面的綜合作用要優(yōu)于單項因子的獨立作用。整體來看，產業(yè)結構與科技創(chuàng)新、環(huán)境規(guī)制對土地綠色利用與經濟增長協(xié)調發(fā)展的交互影響力顯著區(qū)別于其他影響因子的交互影響力。其中第三產業(yè)占GDP比重與工業(yè)固體廢物綜合利用率的交互作用值最大，為0.955 6；第三產業(yè)占GDP比重與專利申請授權量的交互作用值次之，為0.911 1。

表3 土地綠色利用-經濟增長耦合協(xié)調影響因子交互探測結果

綜上所述，各驅動因子對京津冀土地綠色利用與經濟增長協(xié)調發(fā)展均具有一定影響，且影響力差異較大，未來要在快速發(fā)展經濟的同時，發(fā)揮高水平驅動因子的主導作用，同時注重與低水平驅動因子之間的融合，增強因子間的綜合影響水平，逐步提高土地綠色利用與經濟增長的協(xié)調發(fā)展水平。

4 結論

運用超效率SBM模型、標準差橢圓法、耦合協(xié)調度模型、地理探測器等方法，探討了京津冀13個城市2006—2016年土地綠色利用與經濟增長的時空變化規(guī)律、耦合協(xié)調的時空異質性及其影響因素，結論如下：

(1)近10 a京津冀各城市土地綠色利用與經濟增長水平隨時間變化呈現(xiàn)顯著的區(qū)域差異。其中北京、承德等中、北部區(qū)域的土地綠色利用水平較高，均值均在1.00以上；北京、天津等中、東部地區(qū)的經濟增長水平較高，均值均在0.20以上。土地綠色利用發(fā)展主方向由北(偏東)—南(偏西)變?yōu)楸?偏西)—南(偏東)，經濟增長主方向始終為北(偏東)—南(偏西)，其移動趨勢分別為東北—西北向、東北—北向，但在空間上2個系統(tǒng)總體呈東北—西南格局分布格局。

(2)由耦合協(xié)調度測算結果可知，京津冀城市群的土地綠色利用與經濟增長耦合協(xié)調度在0.34～0.41之間，隨時間變化整體協(xié)調度不斷提升，由輕度失調轉變?yōu)闉l臨失調，但依舊處于較低水平。耦合協(xié)調水平在空間分布上具有集中性與穩(wěn)定性，但存在明顯的區(qū)域差異，表現(xiàn)出顯著的空間異質性，即北京、天津遙遙領先，環(huán)京津、沿海區(qū)域次之，冀西北和冀南區(qū)域較差，其中北京、天津表現(xiàn)為土地綠色利用滯后型狀態(tài)，河北省各市表現(xiàn)為經濟增長滯后型狀態(tài)。

(3)土地綠色利用與經濟增長協(xié)調發(fā)展受城鎮(zhèn)化率、專利申請授權量、第三產業(yè)占GDP比重等多因子共同驅動，其中城鎮(zhèn)化率的解釋力最大，且任意兩因子的交互作用均比單因子作用強。

由于京津冀土地綠色利用-經濟增長耦合協(xié)調性較弱，未來應繼續(xù)運用政府政策調控機制，協(xié)調好土地資源利用與經濟增長之間的關系，推動以綠色化為重點的技術改造，減少能源損耗與環(huán)境污染問題。結合城市功能定位，需因地制宜地對不同區(qū)域采取差別化的管理措施，如北京、天津需轉變經濟增長方式，優(yōu)化產業(yè)結構，提高土地綠色利用效率，同時將科技創(chuàng)新作為推動經濟發(fā)展的根本動力，減少對資源的依賴，降低污染物排放；河北省應推動城市土地資源利用方式向節(jié)約、集約的綠色化方向發(fā)展，進一步推動經濟增長向低投入、低排放和高產出的綠色發(fā)展方式轉變，破解保障發(fā)展與保護資源的兩難局面。

對2006—2016年京津冀土地綠色利用與經濟增長耦合協(xié)調的時空異質性及其影響因素進行分析，可為土地生態(tài)利用與經濟協(xié)調發(fā)展提供新的思路，同時對其他城市群土地綠色利用研究具有一定參考價值。未來可通過甄選優(yōu)化評價指標體系，加強遙感數(shù)據、土地利用空間數(shù)據等多元數(shù)據的運用與探索，實現(xiàn)對研究結果空間異質性的深入分析。目前將土地綠色利用與經濟增長結合起來的研究較少，在對耦合協(xié)調發(fā)展的時空變化特征進行分析評價的基礎上，仍需對2個系統(tǒng)間的交互作用進行系統(tǒng)、全面的研究。