太行山區(qū)土地利用沖突測度及其時空特征

薛朝浪,趙宇鸞

(貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院,貴州 貴陽 550001)

1 引言

土地資源是人類社會進(jìn)行物質(zhì)生產(chǎn)的空間載體與生產(chǎn)資料,土地利用的實(shí)質(zhì)是人類對地球表面空間格局的改造行為與作為生產(chǎn)要素的價值索取,短時期內(nèi)變化受區(qū)域經(jīng)濟(jì)社會活動驅(qū)動作用顯著[1]。隨著我國城鄉(xiāng)空間結(jié)構(gòu)的急劇變遷、生產(chǎn)方式的深刻變革,城鎮(zhèn)化與工業(yè)化的發(fā)展推動城鎮(zhèn)人口集聚愈發(fā)顯著,進(jìn)而導(dǎo)致土地開發(fā)外延擴(kuò)張迅速,趨使“生產(chǎn)—生活—生態(tài)”用地時空格局嬗變[2]。因此,在土地資源數(shù)量固定且有限的情況下,城鎮(zhèn)建設(shè)脅迫農(nóng)業(yè)生產(chǎn)與生態(tài)服務(wù)空間萎縮加劇,區(qū)域土地系統(tǒng)利用結(jié)構(gòu)和功能失衡逐漸凸顯,土地利用沖突日益成為制約區(qū)域可持續(xù)發(fā)展和國土空間功能優(yōu)化的瓶頸[3,4]。

目前,土地利用沖突問題受到了國內(nèi)外專家學(xué)者的普遍關(guān)注[3-5]。眾多學(xué)者在土地利用沖突的概念內(nèi)涵[3,4,6]、分類與特征[7-9]、影響機(jī)制[10-12]、診斷模型[13-18]等方面進(jìn)行了大量研究,研究對象或范圍主要關(guān)注城市群[14,16]與大城市邊緣區(qū)[11,15]、礦糧復(fù)合區(qū)[13,17,18]、自然保護(hù)區(qū)[11,19]、風(fēng)景名勝區(qū)[20]的土地利用沖突問題。然而,學(xué)者們對復(fù)雜性和異質(zhì)性較大的山區(qū)土地利用沖突研究還略顯不足[21],尤其是缺乏對整個山區(qū)內(nèi)部不同行政單元土地利用沖突的空間分異特征的探討,因此不利于針對性地提出山區(qū)土地可持續(xù)利用的對策與緩解土地利用沖突具體措施。

太行山區(qū)是我國京津冀地區(qū)國土資源開發(fā)的重要組成部分之一。隨著京津冀協(xié)同發(fā)展戰(zhàn)略的進(jìn)一步推進(jìn),在華北平原經(jīng)濟(jì)發(fā)展加速融合的背景下,太行山區(qū)土地資源利用的沖突問題不容忽視。因此,研究以地理區(qū)位優(yōu)勢突出、土地資源利用較緊張的太行山區(qū)為例,借鑒并嘗試構(gòu)建基于PSR評估框架的土地利用沖突強(qiáng)度診斷模型,通過探索性空間數(shù)據(jù)分析方法(ESDA),對1990—2015年太行山區(qū)的土地利用沖突水平進(jìn)行了定量測度,并對各縣級行政單元土地利用沖突強(qiáng)度演變的時空分異特征進(jìn)行了探討,以期能為我國山區(qū)土地利用沖突強(qiáng)度測度研究提供案例補(bǔ)充,同時為我國太行山區(qū)土地資源的可持續(xù)利用、區(qū)域協(xié)調(diào)發(fā)展、國土空間功能優(yōu)化提供決策依據(jù)。

2 數(shù)據(jù)與方法

2.1 研究區(qū)概況

太行山北起北京西山,向南延伸至河南與山西交界的黃河北岸,西至汾河谷地,東臨華北平原[22],涉及北京、山西、河北、河南4省(市)101個縣級行政區(qū),總面積約13.69×104km2(圖1)。太行山是我國黃土高原和華北平原的地理分界線,受地質(zhì)構(gòu)造影響,西部斷陷抬升形成了山西高原,東部斷陷下降形成了華北平原,斷層?xùn)|西兩側(cè)最大高差達(dá)1500—2000m,平均海拔約861m。太行山區(qū)是太行山地及其與之相鄰的丘陵、平原地區(qū),地勢西高東低,地形起伏大,區(qū)域內(nèi)亞高山、中山、低山、盆地、丘陵等山地類型多樣;氣候?qū)倥瘻貛Т箨懶约撅L(fēng)氣候,多年平均氣溫11.1℃,多年平均降水511.1mm,夏季降水集中,具有雨熱同期、四季分明的特征。太行山區(qū)地處京津冀、太原與中原城市群交界地帶,地理區(qū)位優(yōu)勢突出,資源型城市數(shù)量多、分布廣,境內(nèi)有峰峰、井陘、臨城等礦區(qū),礦產(chǎn)資源豐富;太行山區(qū)作為華北地區(qū)的重要生態(tài)屏障,區(qū)域生態(tài)空間功能效益優(yōu)勢明顯,有小五臺山、歷山等國家自然保護(hù)區(qū),滹沱河、桑干河等河流自西向東匯入渤海,是華北平原的重要水源地。隨著我國經(jīng)濟(jì)社會的快速發(fā)展,太行山區(qū)的城鄉(xiāng)建設(shè)面貌顯著提升,但同時也面臨著經(jīng)濟(jì)建設(shè)與生態(tài)保護(hù)的雙重壓力,土地資源的可持續(xù)利用成為區(qū)域協(xié)調(diào)發(fā)展普遍關(guān)注的問題。

圖1 研究區(qū)示意圖

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

本研究中，土地利用數(shù)據(jù)(1990年、2000年、2010年、2015年)由國家重點(diǎn)基礎(chǔ)研究“973”項目“山區(qū)國土空間功能優(yōu)化與調(diào)控對策”課題組提供。該數(shù)據(jù)基于TM影像解譯而來,采用Googel高分影像及野外調(diào)查驗(yàn)證,總精度達(dá)到85%以上?？h級行政單元界線等矢量數(shù)據(jù)來源于中國科學(xué)院山地表生過程與生態(tài)調(diào)控重點(diǎn)實(shí)驗(yàn)室。1990年、2000年、2010年的經(jīng)濟(jì)社會統(tǒng)計數(shù)據(jù)獲取自中國科學(xué)院資源環(huán)境數(shù)據(jù)中心,2015年的經(jīng)濟(jì)社會統(tǒng)計數(shù)據(jù)由各省(市、縣、區(qū))統(tǒng)計年鑒和統(tǒng)計公報整理獲得。本文利用ArcGIS軟件提取太行山區(qū)各縣級行政單元的不同土地利用類型數(shù)量,并與同時期的經(jīng)濟(jì)社會數(shù)據(jù)構(gòu)建土地利用沖突分析數(shù)據(jù)庫,使用Geoda軟件進(jìn)行空間關(guān)聯(lián)性測度。

2.3 測度方法與指標(biāo)構(gòu)建

土地利用沖突指數(shù):土地利用沖突主要是指經(jīng)濟(jì)社會發(fā)展過程中土地系統(tǒng)利用不盡協(xié)調(diào),社會、經(jīng)濟(jì)、生態(tài)效益互為矛盾的現(xiàn)象[4]。可用土地利用沖突綜合指數(shù)(Land Use Conflict Index,LUCI)來定量表征區(qū)域土地利用沖突的整體發(fā)展程度與水平,具體由壓力指數(shù)(TP)、狀態(tài)指數(shù)(TS)和響應(yīng)指數(shù)(TR)構(gòu)成[21](表1)。

表1 土地利用沖突綜合指數(shù)計算方法

注:r1,r2,…,rn為每類指數(shù)對應(yīng)的權(quán)重;wi,j,k為單個指標(biāo)對應(yīng)的權(quán)重;U(Xi,j,k)為單個指標(biāo)的標(biāo)準(zhǔn)化值;i、j和k分別表示各指標(biāo)的序號。

指標(biāo)體系構(gòu)建:本文基于PSR評估框架,并結(jié)合層次分析法,制定了指標(biāo)構(gòu)建應(yīng)遵循的幾項主要原則:厘清評價的目標(biāo)導(dǎo)向并兼顧單項指標(biāo)系統(tǒng)完整,使各指標(biāo)間既相互獨(dú)立又整體關(guān)聯(lián);參考已有的指標(biāo)體系[13,21],依據(jù)科學(xué)性自上而下編制指標(biāo)分類;梳理各指標(biāo)的典型代表性,刪減冗余指標(biāo),征詢領(lǐng)域內(nèi)專家意見并采用Delphi法[23]確定各指標(biāo)權(quán)重,由此構(gòu)建適宜太行山區(qū)土地利用沖突評價的指標(biāo)體系(表2)。

表2 土地利用沖突強(qiáng)度診斷的PSR評價指標(biāo)體系構(gòu)建

具體來看,各指標(biāo)的選取應(yīng)與太行山區(qū)的經(jīng)濟(jì)社會發(fā)展要素相匹配,通過分析區(qū)域土地利用變化驅(qū)動力,對反映區(qū)域土地系統(tǒng)利用脅迫程度、利用格局、策略行為的影響因子進(jìn)行篩選。經(jīng)濟(jì)發(fā)展與人口增長是土地利用沖突發(fā)生與發(fā)展的根本原因,山區(qū)則因地域的特殊性而更加顯著,因此選取耕地復(fù)種指數(shù)、農(nóng)業(yè)機(jī)械總動力、水果面積占種植業(yè)比重、大牲畜年末存欄數(shù)、鄉(xiāng)村勞動力總量、農(nóng)民人均年收入為土地利用沖突壓力指標(biāo),用以從農(nóng)牧業(yè)發(fā)展?fàn)顩r、城鄉(xiāng)人口結(jié)構(gòu)、經(jīng)濟(jì)收入水平方面量化土地系統(tǒng)承載壓力;土地利用結(jié)構(gòu)是土地資源具有多宜性特征的表現(xiàn)形式,而功能作用是土地作為生產(chǎn)要素的本質(zhì)體現(xiàn),不同社會團(tuán)體與個人根據(jù)自身利益需求,在一定時期和一定程度上對土地利用形態(tài)進(jìn)行改造,以滿足利益需求最大化,因此土地利用狀態(tài)指標(biāo)從用地結(jié)構(gòu)和利用程度方面進(jìn)行考量,選取墾殖系數(shù)、建設(shè)用地比例、勞均耕地面積、未利用地面積比值、土地利用綜合程度指數(shù)來反映區(qū)域土地利用狀況。同時,各利益相關(guān)方為提高土地利用效益、發(fā)揮土地基本功能效用、維持土地生態(tài)效益和可持續(xù)性,進(jìn)而采取工程技術(shù)手段或政策措施等方式對土地資源的利用加以管制。由于政策性因素難以量化,研究選取灌溉面積占耕地比例、化肥投入水平、林草地比例作為土地利用沖突響應(yīng)指標(biāo)。

指標(biāo)中，土地利用綜合程度指數(shù)的計算公式為[24]:

(1)

式中,La為土地利用綜合程度指數(shù);Ci為第i級土地利用程度分級面積百分比;Ai為第i級土地利用類型分級指數(shù),各土地利用類型分級指數(shù)按受人類活動影響程度及生態(tài)效益進(jìn)行賦值:未利用地為1,林草水用地為2,耕地為3,城鎮(zhèn)、農(nóng)村居民點(diǎn)、工礦等建設(shè)用地為4。

由于各指標(biāo)具有不同的衡量單位,為使其評估結(jié)果具有可靠性,需對指標(biāo)原始數(shù)值進(jìn)行標(biāo)準(zhǔn)化處理。依據(jù)各指標(biāo)對評價結(jié)果的正負(fù)向作用,采用max-min標(biāo)準(zhǔn)化方法,公式為[25]:

(2)

式中,xij為i個評價單元的第j個指標(biāo)數(shù)值；i=1,2,…,n;j=1,2,…,m。

2.4 空間關(guān)聯(lián)分析

Global Moran′s I指數(shù)計算公式為[27]:

(3)

(4)

(5)

(6)

式中,E(I)為Moran′s I的期望值;Var(I)為方差。當(dāng)|ZScore|>1.96(α=0.05)時,屬性值在空間分布上存在顯著的空間關(guān)聯(lián)性。Moran′s I的值域?yàn)閇-1,1],可知Moran′s I指數(shù)為一個從-1—1之間連續(xù)變化的指數(shù)。若Moran′s I指數(shù)通過顯著性水平檢驗(yàn),則存在3個空間關(guān)聯(lián)性判別階段:①當(dāng)Moran′s I>0時,存在正的空間關(guān)聯(lián)性,地理要素呈空間集聚態(tài)勢,值越大，空間集聚趨勢越顯著;②當(dāng)Moran′s I<0時,存在負(fù)的空間關(guān)聯(lián)性,地理要素呈空間離散狀態(tài),即空間異質(zhì)性增強(qiáng),值越小，空間異質(zhì)性越顯著;③當(dāng)Moran′s I=0時,則呈隨機(jī)分布,不存在空間關(guān)聯(lián)性。

(7)

(8)

3 結(jié)果及分析

3.1 土地利用沖突強(qiáng)度變化

根據(jù)表1中各指數(shù)公式,計算得到1990—2015年太行山區(qū)各土地利用沖突指數(shù)量化值(圖2)。從壓力指數(shù)來看,1990—2015年土地利用壓力指數(shù)呈倒“U”型變化態(tài)勢。1990—2010年從0.287緩慢上升至0.434,上升幅度為51.22%,年均增加2.56%;而2015年土地利用壓力大幅跌至0.211,較2010年降低51.38%,年均降幅達(dá)10.28%?？芍?0年土地利用增長緩慢,后五年下降快速,土地系統(tǒng)壓力有所緩解。從狀態(tài)指數(shù)來看,土地利用狀態(tài)變化態(tài)勢隨時間發(fā)展呈“N”型上下波動,1990—2000年從0.160緩慢上升到0.208,2010年大幅降到0.129,反映出1990—2010年該時段內(nèi)土地利用結(jié)構(gòu)有較大調(diào)整,區(qū)域建設(shè)用地增長顯著,耕地及草地大量減少(圖3)。2015年土地利用狀態(tài)指數(shù)又緩慢增至0.147,呈平穩(wěn)上升趨勢,這可能與近年來經(jīng)濟(jì)社會的快速發(fā)展趨使區(qū)域用地需求增長有關(guān)。從響應(yīng)指數(shù)來看,1990—2000年變化較穩(wěn)定,土地利用響應(yīng)指數(shù)值僅增加0.003,2000—2015年呈緩慢下降趨勢,降低了0.057,反映太行山區(qū)土地利用保護(hù)和生態(tài)保育水平進(jìn)一步提高,土地生態(tài)環(huán)境進(jìn)一步改善,這可能與國家開展退耕還林、生態(tài)建設(shè)及農(nóng)業(yè)水利基礎(chǔ)設(shè)施建設(shè)等因素有關(guān)。

圖2 1990—2015年太行山區(qū)各土地利用沖突指數(shù)變化狀況

圖3 1990—2015年太行山區(qū)各土地利用類型變化狀況

太行山區(qū)的土地利用沖突綜合指數(shù)總體呈下降趨勢,區(qū)域用地沖突及矛盾有所緩解。具體來看,1990—2010年土地利用沖突指數(shù)呈緩慢上升趨勢,由0.219上升至0.282,增幅28.77%,年均增加1.44%,表明區(qū)域土地系統(tǒng)利用沖突逐漸上升;2010年后土地利用沖突指數(shù)快速降到2015年的0.166,降幅達(dá)41.13%,年均降低8.22%,土地系統(tǒng)利用沖突狀況大幅緩解?？梢?1990—2015年太行山區(qū)土地利用沖突綜合指數(shù)總體呈減緩趨勢,而在時間段上的變化有所區(qū)別,前20年受區(qū)域經(jīng)濟(jì)社會開發(fā)的影響，土地利用沖突積累上升；后五年可能受到國家宏觀政策調(diào)控和土地可持續(xù)性利用認(rèn)識的影響,在土地生態(tài)保護(hù)措施、土地用途管制、鄉(xiāng)城人口遷移等因素作用下,土地系統(tǒng)利用矛盾進(jìn)一步持續(xù)向好的方向轉(zhuǎn)變。

3.2 土地利用強(qiáng)度變化的空間分異特征

總體空間分異特征:利用GeoDa軟件分別計算太行山區(qū)各土地利用沖突指數(shù)的Global Moran′s I值，見表3。

表3 太行山區(qū)不同時段各土地利用沖突指數(shù)的Global Moran′s I值

按顯著性水平95%(α=0.05)的判斷標(biāo)準(zhǔn)對各指數(shù)的Global Moran′s I值進(jìn)行檢驗(yàn),其值在理論上達(dá)到了數(shù)理統(tǒng)計意義的顯著性水平,可知各指數(shù)均具有空間集聚性?？傮w而言,1990—2015年的25年間,土地利用沖突壓力指數(shù)空間集聚特征最為顯著,Global Moran′s I值為0.3059;綜合指數(shù)和狀態(tài)指數(shù)的Global Moran′s I值次之,分別為0.2811、0.1739;響應(yīng)指數(shù)空間集聚性較小,Global Moran′s I值僅為0.0969。具體來看,土地利用沖突壓力指數(shù)的Global Moran′s I值從時段Ⅰ的0.2677上升至?xí)r段Ⅱ的0.3202,增加了0.0525,表明區(qū)域土地系統(tǒng)利用壓力空間集聚趨勢正逐漸增強(qiáng),而時段Ⅲ該指數(shù)的Global Moran′s I值反降到0.1344,減少0.1858,減幅58.03%,空間集聚性大幅減弱;土地利用沖突狀態(tài)指數(shù)的Global Moran′s I值由時段Ⅰ的0.3058減少到時段Ⅱ的0.2388,減少0.067,空間集聚性小幅降低;時段Ⅲ該指數(shù)進(jìn)一步快速減到0.0049,相比時段Ⅱ減少0.2339,降幅達(dá)97.95%,說明狀態(tài)指數(shù)空間集聚性正隨時間的推移趨于向異質(zhì)性轉(zhuǎn)變,反映出區(qū)域土地利用結(jié)構(gòu)空間分異愈發(fā)顯著;土地利用響應(yīng)指數(shù)的Global Moran′s I值由時段Ⅰ的0.2275減少到時段Ⅱ的0.1519,時段Ⅲ再減到0.0023,同樣說明該指數(shù)空間集聚趨勢在減弱。

總體來看,太行山區(qū)土地利用沖突綜合指數(shù)在空間分布上呈集聚趨勢,但該趨勢隨著時間的推移而逐漸減弱。綜合指數(shù)在時段Ⅰ、時段Ⅱ、時段Ⅲ的Global Moran′s I值分別為0.2591、0.2929、0.1767,其動態(tài)變化過程呈先增后減的路徑,時段Ⅰ—Ⅱ緩慢增加,增加了0.0338,到時段Ⅲ減少了0.1162,說明1990—2010年太行山區(qū)經(jīng)濟(jì)社會發(fā)展受資源型產(chǎn)業(yè)影響,“先城后礦”模式導(dǎo)致工業(yè)化與城鎮(zhèn)化建設(shè)過程中的土地利用沖突加劇;而2010—2015年由于資源型城市轉(zhuǎn)型發(fā)展與產(chǎn)業(yè)結(jié)構(gòu)調(diào)整,加之國家土地用途管制措施大力實(shí)施,區(qū)域土地利用沖突強(qiáng)度有所減緩,隨著綠色發(fā)展政策推行,未來太行山區(qū)的土地利用沖突可能將進(jìn)一步緩解。

從圖4可知,熱點(diǎn)區(qū)總體呈“山西—河南—河北”的遷移走勢,說明山西省的土地系統(tǒng)承載壓力有所緩解,河北省的土地系統(tǒng)承載壓力正在逐步增加。具體來看,1990—2000年,太行山區(qū)的土地利用沖突壓力指數(shù)熱點(diǎn)區(qū)主要集聚于山西省長治盆地及其周邊縣域,這與該時期區(qū)域的資源型經(jīng)濟(jì)發(fā)展有關(guān)。長治地區(qū)礦產(chǎn)資源豐富、資源型城市眾多,趨使城鎮(zhèn)建設(shè)及經(jīng)濟(jì)發(fā)展快速,土地系統(tǒng)利用壓力增大;2000—2010年,太行山區(qū)土地利用沖突壓力指數(shù)的熱點(diǎn)區(qū)逐漸遷移到山西與河南兩省交界的黃河經(jīng)濟(jì)帶,在工業(yè)經(jīng)濟(jì)建設(shè)引導(dǎo)下,區(qū)域鋼鐵、化工、冶金等主導(dǎo)優(yōu)勢產(chǎn)業(yè)驅(qū)動城市化、工業(yè)化進(jìn)程發(fā)展快速,土地系統(tǒng)壓力上升顯著;2010—2015年,太行山區(qū)土地利用沖突壓力指數(shù)的熱點(diǎn)區(qū)進(jìn)一步遷移至河北省中南部地區(qū),尤其以石家莊及其周邊縣域最為顯著,這與京津地區(qū)產(chǎn)業(yè)轉(zhuǎn)移及人口疏解有關(guān)。隨著京津冀協(xié)同發(fā)展戰(zhàn)略進(jìn)一步推進(jìn),區(qū)域交通流、人口流、物質(zhì)流等各種流交匯加速,趨使產(chǎn)業(yè)發(fā)展承載用地和居住用地需求旺盛,各業(yè)用地矛盾加劇。冷點(diǎn)區(qū)趨于分散,主要集聚于太行山區(qū)地形崎嶇且經(jīng)濟(jì)發(fā)展落后的縣域,這可能與自身經(jīng)濟(jì)發(fā)展后勁不足、域外城鎮(zhèn)化吸引山區(qū)人口外遷有關(guān)。

圖4 太行山區(qū)土地利用沖突壓力指數(shù)的時空演化特征

太行山區(qū)土地利用沖突狀態(tài)指數(shù)的冷、熱點(diǎn)演化趨勢見圖5。從圖5可見,1990—2000年、2000—2010年太行山區(qū)土地利用狀態(tài)指數(shù)的冷、熱點(diǎn)區(qū)趨于集中,熱點(diǎn)區(qū)主要從河北邯鄲、邢臺、石家莊、保定、張家口等地級市所轄縣域向山西長治、運(yùn)城、陽泉等地級市所轄縣域遷移；冷點(diǎn)區(qū)的遷移趨勢與熱點(diǎn)區(qū)相反,從山西長治、運(yùn)城、陽泉等地級市所轄縣域向河北邯鄲、邢臺、石家莊、保定和河南安陽、鶴壁、新鄉(xiāng)、洛陽等地級市所轄縣域遷移。河北、山西兩省的土地利用系統(tǒng)狀態(tài)變化較劇烈。2010—2015年,該區(qū)域的土地利用沖突狀態(tài)指數(shù)冷、熱點(diǎn)區(qū)趨于分散,熱點(diǎn)區(qū)主要分布于石家莊和邢臺部分縣域,而冷點(diǎn)區(qū)主要分布在長治盆地周邊和晉北地區(qū)的部分縣域。

圖5 太行山區(qū)土地利用沖突狀態(tài)指數(shù)的時空演化特征

太行山區(qū)土地利用沖突響應(yīng)指數(shù)的冷、熱點(diǎn)演化趨勢見圖6。從圖6可見,1990—2000年該區(qū)域熱點(diǎn)區(qū)主要集聚于山西陽泉、忻州、大同、運(yùn)城等所轄縣域,2000—2010年逐漸遷移至河北張家口、保定與河南新鄉(xiāng)、焦作所轄地區(qū),2010—2015年分布范圍有所縮小,以邢臺、晉中、安陽所轄縣域較集中。冷點(diǎn)區(qū)趨于分散,其范圍主要集中在山西境內(nèi),以山地區(qū)縣域?yàn)橹鳌?/p>

圖6 太行山區(qū)土地利用沖突響應(yīng)指數(shù)的時空演化特征

總體來說,太行山區(qū)土地利用沖突綜合指數(shù)的熱點(diǎn)區(qū)、次熱區(qū)趨于向京津冀地區(qū)遷移,而冷點(diǎn)區(qū)、次冷區(qū)則趨于向山西廣大區(qū)域轉(zhuǎn)移(圖7)。具體來看,1990—2000年太行山區(qū)土地系統(tǒng)利用綜合指數(shù)的熱點(diǎn)區(qū)主要集中分布在長治地區(qū)及其周邊區(qū)域,這與區(qū)域資源經(jīng)濟(jì)的建設(shè)有關(guān);隨著產(chǎn)業(yè)改革和政策導(dǎo)向指引,2000—2010年土地利用熱點(diǎn)區(qū)逐漸遷移至山西與河南兩省交界的運(yùn)城、洛陽、安陽、鶴壁、新鄉(xiāng)、焦作等所轄縣域;2010—2015年進(jìn)一步遷移至石家莊及其與京津相鄰的區(qū)域,這與京津冀協(xié)同發(fā)展戰(zhàn)略背景下的工業(yè)化建設(shè)及城鎮(zhèn)化人口集聚驅(qū)動土地需求增加有關(guān)。

圖7 太行山區(qū)土地利用沖突綜合指數(shù)的時空演化特征

4 結(jié)論與討論

4.1 結(jié)論

研究通過構(gòu)建基于PSR評估框架的土地利用沖突強(qiáng)度診斷模型,對太行山區(qū)的土地利用沖突水平進(jìn)行了定量測度,并對各縣級行政單元土地利用沖突強(qiáng)度演變的時空分異特征進(jìn)行了探討,主要結(jié)論為:①1990—2015年太行山區(qū)的土地利用沖突強(qiáng)度有所緩解,綜合指數(shù)從0.219下降到0.166,呈進(jìn)一步減緩趨勢。其中,1990—2010年土地利用綜合指數(shù)呈緩慢上升態(tài)勢,由0.219上升到0.282,土地利用沖突強(qiáng)度增加,受區(qū)域資源型經(jīng)濟(jì)發(fā)展驅(qū)動作用顯著;2010—2015年,綜合指數(shù)由0.282下降到0.166,土地利用沖突強(qiáng)度大幅減緩,降幅達(dá)41.13%,受國家綠色發(fā)展政策及農(nóng)業(yè)水利等基礎(chǔ)設(shè)施建設(shè)、退耕還林、鄉(xiāng)城人口遷移等因素影響,土地利用沖突可能逐步緩解。②1990—2015年太行山區(qū)土地利用沖突的壓力指數(shù)、狀態(tài)指數(shù)、響應(yīng)指數(shù)及綜合指數(shù)的Global Moran′s I值分別為0.3059、0.1739、0.0969、0.2811,均通過了α=0.05的顯著性水平檢驗(yàn),壓力指數(shù)空間集聚特征最為顯著,響應(yīng)指數(shù)空間集聚性較小。土地利用各指數(shù)時段Ⅲ的Global Moran′s I值均小于時段Ⅱ的值,說明太行山區(qū)時空格局上的土地利用沖突空間異質(zhì)性逐漸增強(qiáng)。③1990—2015年太行山區(qū)各縣級行政單元的土地利用沖突綜合指數(shù)集聚特征顯著。土地利用熱點(diǎn)區(qū)、次熱區(qū)呈山西—河南—河北的遷移走勢,冷點(diǎn)區(qū)、次冷區(qū)趨于向山西廣大區(qū)域轉(zhuǎn)移。土地利用熱點(diǎn)區(qū)趨于向京津冀地區(qū)遷移,以石家莊、張家口、保定所轄縣域最為顯著,受京津冀一體化建設(shè)推動土地需求增加影響較大;土地利用冷點(diǎn)區(qū)主要分布在太行山區(qū)地形崎嶇且經(jīng)濟(jì)發(fā)展落后的縣域,這可能與自身經(jīng)濟(jì)發(fā)展后勁不足,域外城鎮(zhèn)化吸引山區(qū)人口外遷有關(guān)。

4.2 討論

在京津冀協(xié)同發(fā)展戰(zhàn)略的背景下,未來太行山區(qū)的土地利用沖突問題仍將進(jìn)一步凸顯。作為京津地區(qū)產(chǎn)業(yè)轉(zhuǎn)移及人口疏解承載地的河北省太行山區(qū)域,地理區(qū)位優(yōu)勢突出、土地資源豐富,在城鎮(zhèn)化發(fā)展、產(chǎn)業(yè)轉(zhuǎn)移承接地建設(shè)等方面將發(fā)揮重要作用。根據(jù)研究結(jié)論,建議進(jìn)一步關(guān)注太行山低山丘陵區(qū)的土地利用沖突和山區(qū)國土空間功能優(yōu)化問題,加強(qiáng)“京津冀+晉”產(chǎn)業(yè)協(xié)同發(fā)展規(guī)劃及生態(tài)聯(lián)同保護(hù)機(jī)制的研究和建設(shè),開展以燕山—太行山區(qū)為主體的京冀晉豫四省生態(tài)經(jīng)濟(jì)走廊建設(shè)的研究與探討。同時,受太行山區(qū)經(jīng)濟(jì)發(fā)展落后的縣域人口外遷和農(nóng)戶生計轉(zhuǎn)型的影響,太行山區(qū)土地利用的沖突有所緩解,但地形、發(fā)展底子薄等多種制約因素疊加的影響,太行山區(qū)的土地利用沖突問題不容忽視,特別是環(huán)京津貧困帶內(nèi)集中連片特困地區(qū)的土地利用沖突問題在精準(zhǔn)脫貧中應(yīng)給予關(guān)注?？傊?隨著城市化與工業(yè)化的快速發(fā)展,太行山區(qū)的土地資源可持續(xù)利用及土地利用沖突緩解仍將是該地區(qū)協(xié)調(diào)發(fā)展需要關(guān)注的焦點(diǎn)。