錢田田, 王宏衛(wèi), 余芳瑞, 魏一鳴, 尹玉聰, 薛孟琦

(1. 新疆大學 地理與遙感科學學院, 新疆 烏魯木齊 830017;2. 新疆大學 新疆綠洲生態(tài)自治區(qū)重點實驗室, 新疆 烏魯木齊 830017)

新疆是典型的干旱半干旱區(qū),自然條件較為惡劣.城市作為人口、產(chǎn)業(yè)主要集聚區(qū),是新疆生態(tài)保護的先行區(qū)和主戰(zhàn)場.在生態(tài)環(huán)境脆弱和外來擾動的雙重壓力下,城市生態(tài)保護和治理的核心任務是增強韌性轉(zhuǎn)變,因此,提高城市應對生態(tài)風險的能力成為城市高質(zhì)量發(fā)展的一個重要內(nèi)容,也是韌性城市未來的發(fā)展方向.

韌性作為一種主動探索適應性的調(diào)整變化方式[1],與區(qū)域高質(zhì)量發(fā)展理念存在共通之處[2],在區(qū)域不斷面對外界干擾和沖擊的影響下[3],給區(qū)域發(fā)展提供了新的思路.韌性從工程韌性、生態(tài)韌性發(fā)展到演化韌性(社會生態(tài)韌性).作為城市韌性的重要維度,城市生態(tài)韌性受到國內(nèi)外學者的廣泛關(guān)注.國外生態(tài)韌性的相關(guān)研究主要圍繞韌性理論[4]、形成機制[5]、優(yōu)化路徑[6]等內(nèi)容展開.Holling[7]將生態(tài)韌性定義為生態(tài)系統(tǒng)在應對自然或人為因素變化時所產(chǎn)生的持續(xù)應對性;André等[8]嘗試將生態(tài)韌性與城市治理績效相結(jié)合,構(gòu)建復雜系統(tǒng)下的城市生態(tài)韌性提升體系.國內(nèi)研究則主要涉及韌性演化[9]、指標構(gòu)建[10]和影響因素[11]等方面.夏楚瑜等[12]從抵抗力、適應力和恢復力3個方面構(gòu)建杭州市生態(tài)韌性評估模型,對杭州市生態(tài)韌性時空變化格局進行分析,并提出分區(qū)管理模式;王婷等[13]從我國高質(zhì)量發(fā)展的需求出發(fā),構(gòu)建基于高質(zhì)量發(fā)展的城市生態(tài)韌性測度框架;王松茂等[14]基于演化韌性理論,在評估和分析山東半島城市群城市生態(tài)韌性的時空分異和演進規(guī)律的基礎(chǔ)上,采用障礙度模型識別其主要障礙因子.現(xiàn)有的城市生態(tài)韌性研究限于概念的不統(tǒng)一,側(cè)重分析城市生態(tài)系統(tǒng)在受到外來干擾和沖擊后的抵御和恢復能力,較少有在演進韌性視角下對城市生態(tài)韌性進行深入探討,難以綜合評估城市生態(tài)韌性的內(nèi)涵機制和高質(zhì)量發(fā)展能力.

鑒于此,本文運用驅(qū)動力-壓力-狀態(tài)-影響-響應(DPSIR)模型、熵權(quán)法、綜合指數(shù)法等評估2000-2020年新疆城市的生態(tài)韌性,通過地理探測器探測其影響因子,并提出城市生態(tài)韌性的提升策略.

1 研究數(shù)據(jù)及方法

1.1 研究區(qū)概況

研究區(qū)概況圖,如圖1所示.新疆地廣人稀,城市作為人類活動的主要載體,在具有豐富資源儲備的同時,存在生態(tài)環(huán)境脆弱,生態(tài)承載力低的先天不足問題,社會經(jīng)濟行為受到客觀自然環(huán)境制約,工業(yè)生產(chǎn)和農(nóng)業(yè)活動的環(huán)境依附性需要保持在較低水平.2000-2020年間,新疆在西部大開發(fā)、對口援疆等政策的引導下,經(jīng)濟水平不斷提升,國民生產(chǎn)總值由1 363.56億元上升到13 797.58億元;基礎(chǔ)設(shè)施建設(shè)不斷完善,到2020年,新疆高速公路里程已達到5 412.9 km;生態(tài)環(huán)境遭到破壞,近20年來新疆生態(tài)環(huán)境的主要情況是局部(內(nèi)部綠洲)改良,然而整體逐漸惡化,尤其是綠洲周圍的沙漠過渡帶不斷惡化,對綠洲生態(tài)環(huán)境造成明顯的負面影響.

圖1 研究區(qū)概況圖Fig.1 Overview of study area

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

遵循指標數(shù)據(jù)獲取的可得性、科學性和系統(tǒng)性,運用DPSIR模型構(gòu)建新疆城市生態(tài)韌性評估體系.研究數(shù)據(jù)分別來源于美國航空航天局(NASA)、大氣成分分析組、中國科學院資源環(huán)境科學數(shù)據(jù)中心網(wǎng)、NOAA/NCEI美國國家環(huán)境信息中心、北京大學地理數(shù)據(jù)平臺、Openstreetmap、新疆統(tǒng)計年鑒(2001-2021)、中國縣域統(tǒng)計年鑒(2001-2021)、根據(jù)2001-2021年各地市統(tǒng)計年鑒、縣域統(tǒng)計年鑒和部分縣級以下統(tǒng)計公報進行補充.由于數(shù)據(jù)獲取來源較多,對數(shù)據(jù)進行統(tǒng)一柵格與分辨率處理,重采樣后的數(shù)據(jù)分辨率為1 000 m.

1.3 研究方法

1.3.1 DPSIR模型 DPSIR模型[15]強調(diào)社會經(jīng)濟發(fā)展與環(huán)境之間的關(guān)系.在DPSIR模型中,經(jīng)濟和社會因子作為驅(qū)動力(D),增加或減輕生態(tài)環(huán)境壓力(P),造成城市社會-經(jīng)濟-生態(tài)系統(tǒng)狀態(tài)(S)發(fā)生改變,這些改變對城市資源質(zhì)量、產(chǎn)業(yè)結(jié)構(gòu)、區(qū)域發(fā)展等方面產(chǎn)生影響(I),并使系統(tǒng)以預防、適應或改善的方式做出響應(R).響應又會反過來促進社會經(jīng)濟發(fā)展,調(diào)控環(huán)境壓力,并改變城市當前的狀態(tài).DPSIR模型,如圖2所示.

圖2 DPSIR模型Fig.2 DPSIR model

1) 指標體系構(gòu)建.為準確、客觀地評價新疆城市生態(tài)韌性,參考文獻[16-17]的研究,綜合考慮新疆社會經(jīng)濟發(fā)展、生態(tài)環(huán)境變化、基礎(chǔ)設(shè)施建設(shè)等特征,構(gòu)建DPSIR模型下新疆城市生態(tài)韌性評價指標體系,共計24項指標,如表1所示.

表1 DPSIR模型下新疆城市生態(tài)韌性評價指標體系Tab.1 Evaluation index system of urban ecological resilience in Xinjiang under DPSIR model

2) 數(shù)據(jù)標準化處理.為了消除原始數(shù)據(jù)量綱和數(shù)值間的差異,需要對數(shù)據(jù)進行標準化處理,即

(1)

式(1)中：x′為原始數(shù)據(jù)進行標準化處理后的數(shù)值;xi,j為指標的實際值;max(xi,j)為該指標序列的最大值;min(xi,j)為該指標序列的最小值.

3) 熵權(quán)法.在進行新疆城市生態(tài)韌性評價時,采用客觀賦權(quán)方法中的熵權(quán)法計算指標權(quán)重[19],即

(2)

(3)

dj=1-ej,

(4)

(5)

式(2)～(5)中：pi,j為第j項指標下第i個生態(tài)韌性指標值所占的比重;m為新疆城市總數(shù);ej為第j項指標的熵值;dj為第j項指標的冗余度;wj為第j項指標的權(quán)重;n為指標個數(shù).

4) 新疆城市生態(tài)韌性評價.使用綜合指數(shù)法計算新疆城市生態(tài)韌性,表達式為

(6)

式(6)中：ERI為新疆城市生態(tài)韌性指數(shù);Xj為指標j的標準化值.

1.3.2 全局Moran′sI指數(shù) 全局空間自相關(guān)主要用于分析地理要素在整個區(qū)域內(nèi)的關(guān)聯(lián)特征,通常用全局Moran′sI指數(shù)(IM)來衡量,其表達式為

(7)

1.3.3 地理探測器 地理探測器是由王勁峰等[20]提出的用于空間數(shù)據(jù)探索的方法,因子探測器主要用來檢測地理因素和某種指標的變化在空間上是否具有顯著的一致性,其表達式為

(8)

2 研究結(jié)果與分析

2.1 新疆城市生態(tài)韌性水平演變特征

根據(jù)所建立的新疆城市生態(tài)韌性評價指標體系,對其生態(tài)韌性水平進行評價,分別得到2000年、2010年和2020年新疆城市生態(tài)韌性評價指數(shù).運用 ArcGIS中的自然斷點法對新疆城市生態(tài)韌性進行分等定級,劃分為5類：Ⅰ類為低韌性(ERI<1.691);Ⅱ類為較低韌性(1.691≤ERI≤1.941);Ⅲ類為中等韌性(1.941≤ERI≤2.201);Ⅳ類為較高韌性(2.201≤ERI≤2.556);Ⅴ類為高韌性(ERI>2.556).

2.1.1 時間演變規(guī)律 計算得到的2000-2020年的新疆城市生態(tài)韌性等級劃分,如圖3所示.由圖3可知：新疆城市生態(tài)韌性水平整體得到提升,城市生態(tài)韌性綜合指數(shù)平均值分別為1.885(2000年),2.162(2010年),2.253(2020年).新疆85個城市中生態(tài)韌性水平提升速度最快的是奎屯市,由1.751增加為2.688,而且末縣、塔城市和烏魯木齊縣生態(tài)韌性水平呈現(xiàn)負增長.

(a) 2000年 (b) 2010年 (c) 2020年圖3 2000-2020年新疆城市生態(tài)韌性等級劃分Fig.3 Classification of urban ecological resilience levels in Xinjiang from 2000 to 2020

新疆城市生態(tài)韌性等級變化,如圖4所示.由圖4可知：2000年城市的Ⅰ類生態(tài)韌性水平在2020年大部分優(yōu)化到Ⅱ類和Ⅲ類生態(tài)韌性水平,很難躍遷到Ⅳ類及Ⅴ類生態(tài)韌性水平;2000-2010年,Ⅰ類和Ⅱ類城市數(shù)由53個減少至18個,Ⅲ類、Ⅳ類和Ⅴ類城市數(shù)由32個增加到67個;2010-2020年,Ⅰ類和Ⅱ類城市數(shù)由18個變?yōu)?0個,Ⅲ類城市數(shù)由32個變?yōu)?1個,Ⅳ類城市數(shù)由31個增加到48個,Ⅴ類城市數(shù)由4個增加到6個.產(chǎn)生這一變化的原因主要是20年間新疆城市經(jīng)濟迅速發(fā)展,社會事業(yè)不斷完善,生態(tài)環(huán)境逐步變好.經(jīng)濟上,2000-2020年,新疆85個城市人均GDP總值由7 372元增加到23 845元,作為中國重要的農(nóng)牧業(yè)產(chǎn)地,新疆的農(nóng)業(yè)系統(tǒng)產(chǎn)出結(jié)構(gòu)熵和農(nóng)業(yè)產(chǎn)出系統(tǒng)結(jié)構(gòu)轉(zhuǎn)換率隨著科學技術(shù)水平不斷完善,提升明顯;社會事業(yè)上,地方財政支出由124.96萬元提高至3 480.46萬元,醫(yī)療衛(wèi)生機構(gòu)、福利機構(gòu)、學校和交通道路等公共基礎(chǔ)設(shè)施建設(shè)水平逐步完善;生態(tài)環(huán)境上,單位GDP工業(yè)SO2排放量和單位GDP工業(yè)PM2.5濃度逐步降低,說明在工業(yè)生產(chǎn)上更加注意污染物的排放,但空氣質(zhì)量整體有所下降,土壤侵蝕、植被退化、林草覆蓋率降低等問題不容忽視.

圖4 新疆城市生態(tài)韌性等級變化Fig.4 Changes in urban ecological resilience levels in Xinjiang

從指標體系的5個維度來看,這20年間,驅(qū)動力小幅度增強,說明近20年新疆城市社會經(jīng)濟發(fā)展速度一直保持增長;壓力在顯著增強后又明顯下降,說明2010年后,新疆各城市開始注重排污減污和環(huán)境保護;狀態(tài)不斷變好,說明系統(tǒng)對壓力的轉(zhuǎn)化能力不斷優(yōu)化;影響改變不大,說明在快速發(fā)展的同時,相應的配套設(shè)施也在同步發(fā)展;響應不斷提升,說明基礎(chǔ)設(shè)施和服務能力不斷完善.

2.1.2 空間演變規(guī)律 1) 新疆城市生態(tài)韌性空間相關(guān)性.利用ArcGIS軟件中的空間分析工具,得出2000-2020年新疆城市生態(tài)韌性水平的全局Moran′sI指數(shù),如表2所示.表2中：Z為標準差的倍數(shù);P為顯著性.由表2可知：2000-2020年,新疆城市生態(tài)韌性全局Moran′sI指數(shù)均通過了Z統(tǒng)計量檢驗,Z≥2.58,且P≤0.01,新疆城市生態(tài)韌性之間存在顯著的正相關(guān)關(guān)系,即城市生態(tài)韌性指數(shù)不論高低,在空間上都呈現(xiàn)集聚的特征.

表2 2000-2020年新疆城市生態(tài)韌性全局Moran′s I指數(shù)Tab.2 Global Moran′s I index of urban ecological resilience in Xinjiang from 2000 to 2020

2) 新疆城市生態(tài)韌性空間異質(zhì)性.受制于新疆特殊區(qū)位和自然環(huán)境特征,城市生態(tài)韌性空間異質(zhì)性明顯.在韌性等級劃分中,新疆城市生態(tài)韌性水平呈現(xiàn)出北高南低的特征.生態(tài)韌性水平為Ⅳ類和Ⅴ類的城市主要集中分布于天山北坡經(jīng)濟帶,生態(tài)韌性水平為Ⅰ,Ⅱ,Ⅲ類的城市雖有所提升,但與Ⅳ類和Ⅴ類城市存在一定差異.產(chǎn)生這種現(xiàn)象的主要原因是生態(tài)基礎(chǔ)、經(jīng)濟加持下的政策傾斜.天山北坡經(jīng)濟帶的城市在優(yōu)越的區(qū)位條件與豐富的資源優(yōu)勢的加持下,系統(tǒng)有更強的抗干擾和回穩(wěn)能力.

2020年新疆城市生態(tài)韌性不同維度的空間分異情況,如圖5所示.由圖5可知：從驅(qū)動力維度看,新疆城市驅(qū)動力整體較弱,多集中于等級Ⅱ,Ⅲ,達到等級Ⅴ的只有4個城市,分別為克拉瑪依市、烏魯木齊市、庫爾勒市和伊吾縣,說明這4個城市社會經(jīng)濟發(fā)展的驅(qū)動力較強;從壓力維度看,大部分城市處于高壓力狀態(tài),多處于等級Ⅳ,Ⅴ,僅有若羌縣處于等級Ⅰ,說明該縣城市工業(yè)化水平低,產(chǎn)生的污染較小;從狀態(tài)維度看,大部分城市處于等級Ⅳ,表明在不同驅(qū)動力和壓力影響下,城市生態(tài)系統(tǒng)狀態(tài)表現(xiàn)較好;從影響維度看,等級較為分散,由于自然本底差異,城市生態(tài)系統(tǒng)表現(xiàn)出不同程度的影響水平;從響應維度看,全疆城市整體響應偏低,處于等級Ⅴ的僅有烏魯木齊市,說明新疆城市整體響應能力弱.

(a) ERI (b) 驅(qū)動力 (c) 壓力

2.2 新疆城市生態(tài)韌性影響因素分析

2.2.1 單因子探測分析 采用ArcGIS自然斷點法對新疆城市生態(tài)韌性指標體系中各項指標數(shù)值進行離散化,將離散化后的類型量導入地理探測器中進行因子探測,前8位影響因子的探測結(jié)果,如表3所示.表3中：q為影響因子對城市生態(tài)韌性的決定力.由表3可知：各影響因子對城市生態(tài)韌性的決定力存在較大差異.

表3 前8位影響因子的探測結(jié)果Tab.3 Detection results of top 8 influencing factors

在驅(qū)動力維度中,人均消費水平和城鎮(zhèn)化率能表征城市生態(tài)韌性的動力來源,其中,城鎮(zhèn)化率的q值高于人均消費水平,反映出人口向城市的集中是城市生態(tài)韌性的主要動力來源.城鎮(zhèn)化率和消費水平的提升是城市發(fā)展的主要原因.

在壓力維度中,化肥施用度能夠反映生態(tài)系統(tǒng)所面臨的生態(tài)威脅情況,其q值為0.186,說明新疆城市間受到化肥施用度的影響較大,化肥施用度為城市生態(tài)韌性的主要壓力來源,過度施用化肥對土壤和地下水產(chǎn)生不同程度的污染,直接影響當?shù)刈魑锏漠a(chǎn)量和人們的日常生活.

在狀態(tài)維度中,植被退化和土壤侵蝕能夠反映城市系統(tǒng)的基本情況,其因子決定作用較高,說明植被退化和土壤侵蝕能較大程度地反映新疆城市間的資源環(huán)境情況.

在影響維度中,非農(nóng)產(chǎn)業(yè)區(qū)位商、林草覆蓋率和空氣質(zhì)量指數(shù)能夠反映城市產(chǎn)業(yè)結(jié)構(gòu)、資源質(zhì)量、生態(tài)環(huán)境等方面對生態(tài)韌性的影響,其中,林草覆蓋率因子q值最高,說明林草資源質(zhì)量是城市生態(tài)韌性產(chǎn)生顯著正效應的關(guān)鍵.

2.2.2 交互探測分析 運用地理探測器中的交互探測器探測雙因子對城市生態(tài)韌性的共同作用.新疆城市生態(tài)韌性因子交互探測結(jié)果,如表4所示.圖4中：箭頭表示影響因子兩兩交互的作用值高于單一因子的作用值.

表4 新疆城市生態(tài)韌性因子交互探測結(jié)果Tab.4 Interactive detection results of urban ecological resilience factors in Xinjiang

由表4可知：8個主要影響因子兩兩交互的作用值均高于單一因子的作用值,表明雙因子的共同作用會增加其對生態(tài)韌性水平的解釋力,其中,城鎮(zhèn)化率和林草覆蓋率的交互探測結(jié)果最高,為0.615,表明新疆城市生態(tài)韌性受城鎮(zhèn)化率和林草覆蓋率的共同影響最大.因此,可以從雙因子的交互作用著手,提升城市生態(tài)韌性.

2.3 新疆城市生態(tài)韌性優(yōu)化建議

1) 提升驅(qū)動力.把握天山北坡城市群高質(zhì)量發(fā)展的戰(zhàn)略定位,大力發(fā)展經(jīng)濟,吸引并留住人才,以點帶面,輻射全疆,協(xié)調(diào)推進全疆城市生態(tài)建設(shè),通過技術(shù)溢出、制度溢出和人才溢出,擴大高生態(tài)韌性城市的生態(tài)建設(shè)輻射帶動作用和示范效應,打破區(qū)域城市生態(tài)韌性兩極化困境.

2) 化解壓力.在增強壓力轉(zhuǎn)化的同時,劃定并嚴守生態(tài)保護紅線,加強對社會安全事故及城市節(jié)能減排的管控力度,鼓勵綠色發(fā)展型企業(yè)的轉(zhuǎn)型升級,轉(zhuǎn)變城市高耗能、高污染的發(fā)展模式,堅決不走“先污染,后治理”的老路.

3) 改善狀態(tài).通過制定環(huán)境保護規(guī)章制度和更新產(chǎn)業(yè)結(jié)構(gòu)入手,加快構(gòu)建多元化現(xiàn)代生態(tài)產(chǎn)業(yè)體系,完善劣勢地區(qū)生態(tài)產(chǎn)業(yè)鏈、供應鏈建設(shè),以提高經(jīng)濟穩(wěn)定性,并合理布局應急產(chǎn)業(yè),有效應對各類風險壓力.

4) 優(yōu)化影響.通過對城市用地空間的合理布局,降低社會經(jīng)濟發(fā)展對生態(tài)環(huán)境的壓力,從而有效提升城市生態(tài)系統(tǒng)對壓力的適應能力.

5) 提升響應.加強基礎(chǔ)設(shè)施建設(shè),完善交通道路網(wǎng)絡(luò)體系,提升綜合服務能力,以復合功能為導向布局城市基礎(chǔ)設(shè)施,增強城市生態(tài)承載能力,加大“人才興疆”戰(zhàn)略實施力度,推動科教深度融合創(chuàng)新,加快科技成果資本化、產(chǎn)業(yè)化,有效增強城市生態(tài)系統(tǒng)的學習和轉(zhuǎn)化能力.

3 結(jié)論

2000-2020年新疆發(fā)展迅速,結(jié)合統(tǒng)計數(shù)據(jù)和遙感數(shù)據(jù),建立DPSIR框架下的城市生態(tài)韌性評價指標體系,借助地理探測器對城市生態(tài)韌性影響因子進行探測并分析,得出以下4點結(jié)論.

1) 時間尺度上,隨著生態(tài)理念日益完善和落地,新疆在20年間的城市生態(tài)韌性水平不斷提升,Ⅰ類、Ⅱ類生態(tài)韌性水平的城市明顯減少,Ⅳ類、Ⅴ類生態(tài)韌性水平的城市數(shù)量逐步增多,但跨越式的韌性躍遷尚不明顯.分階段來看,2000-2010年,新疆城市社會經(jīng)濟的發(fā)展給生態(tài)環(huán)境帶來了較大的壓力;2010-2020年,新疆在保持經(jīng)濟快速發(fā)展的同時,開始注重響應措施的提升.

2) 空間分布上,新疆城市生態(tài)韌性之間呈現(xiàn)顯著的空間正相關(guān)性,新疆城市生態(tài)韌性具有明顯的區(qū)域輻射性;新疆城市生態(tài)韌性水平空間異質(zhì)性明顯,天山北坡經(jīng)濟帶中的城市生態(tài)韌性水平相對較高,各時間結(jié)點均表現(xiàn)為北疆城市生態(tài)韌性高,南疆城市生態(tài)韌性低的特征.

3) 影響新疆城市生態(tài)韌性的因子主要為城鎮(zhèn)化率、化肥施用度、植被退化、土壤侵蝕、林草覆蓋率、非農(nóng)產(chǎn)業(yè)區(qū)位商、人均消費水平、空氣質(zhì)量指數(shù);影響因子兩兩交互的作用值均高于單一因子的作用值,其中,城鎮(zhèn)化率和林草覆蓋率的交互探測結(jié)果最高,為0.615.

4) 新疆地域廣闊、自然環(huán)境復雜,在大力發(fā)展經(jīng)濟的同時,應在綠色生態(tài)等思想的引導下,通過合理規(guī)劃城市定位提升驅(qū)動力,借助嚴格劃定生態(tài)保護紅線和發(fā)展綠色企業(yè)等措施化解壓力,結(jié)合制定環(huán)境保護規(guī)章制度和更新產(chǎn)業(yè)結(jié)構(gòu)等策略改善狀態(tài),通過調(diào)整用地空間布局來優(yōu)化影響,并從完善基礎(chǔ)設(shè)施和社會保障等方面提升響應,最終實現(xiàn)新疆城市韌性水平的提升.