張成龍　魯成樹　曹海濤

摘要：文章以皖江城市帶為研究對(duì)象，基于1990年、1995年、2000年、2005年、2010年、2015年六期遙感解譯數(shù)據(jù)，借助Arcgis和GS+軟件對(duì)相應(yīng)數(shù)據(jù)進(jìn)行分析處理，并利用空間自相關(guān)進(jìn)一步確定生態(tài)風(fēng)險(xiǎn)指數(shù)在空間上的依賴性。研究表明：1990～2015年土地利用類型減少最多的為耕地，且主要轉(zhuǎn)換為建設(shè)用地;生態(tài)風(fēng)險(xiǎn)等級(jí)主要以中等和較低等級(jí)為主;1990～2015年morans I指數(shù)不斷增加，且逐漸趨向于1，表明生態(tài)風(fēng)險(xiǎn)指數(shù)空間相關(guān)程度明顯，且有逐步增強(qiáng)的趨勢(shì)。

關(guān)鍵詞：皖江城市帶;土地利用變化;生態(tài)風(fēng)險(xiǎn);空間自相關(guān)

隨著社會(huì)經(jīng)濟(jì)高速發(fā)展的需求以及人口數(shù)量不斷增長(zhǎng)的現(xiàn)實(shí)，對(duì)自然資源過(guò)度的消耗產(chǎn)生了一系列生態(tài)環(huán)境問(wèn)題，已經(jīng)嚴(yán)重威脅了自身的安全與發(fā)展。因此，如何改變當(dāng)下困境成為人們關(guān)注的焦點(diǎn)，生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)作為改變這一困境的重要評(píng)價(jià)方法，也越來(lái)越受到人們的關(guān)注。生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)是生態(tài)系統(tǒng)及其組成部分在受到外界壓力的情況下，對(duì)其可能發(fā)生或正在發(fā)生的不利生態(tài)影響的過(guò)程進(jìn)行評(píng)估的一種方法。近年來(lái)，土地利用變化與生態(tài)風(fēng)險(xiǎn)之間相關(guān)性的研究越來(lái)越多，土地利用變化已經(jīng)作為一種人為干擾源出現(xiàn)在生態(tài)系統(tǒng)中，而生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)作為當(dāng)今風(fēng)險(xiǎn)管理的一種重要科學(xué)手段，也為土地利用生態(tài)效應(yīng)的定量化評(píng)估帶來(lái)了新思路。

當(dāng)前，針對(duì)土地利用生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，在研究方法上主要有半變異函數(shù)分析以及克里格差值、空間自相關(guān)、轉(zhuǎn)移矩陣方法等。在研究模式上主要有兩種，一是基于傳統(tǒng)的風(fēng)險(xiǎn)源匯;二是直接從景觀生態(tài)學(xué)的角度評(píng)價(jià)生態(tài)風(fēng)險(xiǎn)。在研究尺度上主要集中在流域與城鎮(zhèn)上，如高永年以及曾勇分別對(duì)太湖流域以及呼和浩特市區(qū)生態(tài)風(fēng)險(xiǎn)進(jìn)行了評(píng)價(jià);而對(duì)于更大尺度上的經(jīng)濟(jì)帶研究甚少。皖江城市帶是作為泛長(zhǎng)三角地區(qū)的重要組成部分，在中西部產(chǎn)業(yè)轉(zhuǎn)移承接上具有重大戰(zhàn)略意義，對(duì)皖江城市帶土地利用變化的深入研究有利于揭示高速城鎮(zhèn)化背景下土地利用格局的時(shí)空變化特征以及對(duì)生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的影響。

一、研究區(qū)域概況

皖江城市帶總面積約76000km2，包含合肥市、蕪湖市、安慶市、滁州市、馬鞍山市、銅陵市、池州市、宣城市8個(gè)地級(jí)市以及六安市的金安區(qū)以及舒城縣，人口約為3058萬(wàn)人。皖江城市帶是安徽省僅有的以產(chǎn)業(yè)轉(zhuǎn)移為主線的區(qū)域戰(zhàn)略規(guī)劃，它于2010年8月被國(guó)務(wù)院正式納入國(guó)家發(fā)展戰(zhàn)略中，在中西部崛起戰(zhàn)略地位舉足輕重。

二、研究數(shù)據(jù)與研究方法

（一）研究數(shù)據(jù)

該研究所涉及的數(shù)據(jù)有兩類：土地利用數(shù)據(jù)與社會(huì)經(jīng)濟(jì)數(shù)據(jù)。土地利用數(shù)據(jù)包括1990年、1995年、2000年、2005年、2010年及2015年，共計(jì)6期，來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)中心對(duì)研究區(qū)域1990～2015年遙感影像的解譯，柵格大小為1km×1km。

（二）研究方法

1. 土地用途類型劃分

本研究參照胡金龍的研究成果，并結(jié)合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究的現(xiàn)實(shí)需要，將皖江城市帶土地利用類型劃分為：耕地、林地、草地、水域、建設(shè)用地、未利用地共計(jì)6類。

2. 區(qū)域生態(tài)風(fēng)險(xiǎn)分析

首先根據(jù)皖江城市帶的面積，形狀、地理位置等影響因素，采取面狀評(píng)價(jià)與點(diǎn)狀評(píng)價(jià)相結(jié)合的方法，參照相關(guān)研究劃定生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)單元，通過(guò)Arcgis中“漁網(wǎng)工具”設(shè)定不同尺度得評(píng)價(jià)單元，最終選取選取10km×10km的網(wǎng)格評(píng)價(jià)樣地，共計(jì)889個(gè)。通過(guò)生態(tài)風(fēng)險(xiǎn)公式計(jì)算出每個(gè)評(píng)價(jià)樣地的生態(tài)風(fēng)險(xiǎn)指數(shù)，并劃分風(fēng)險(xiǎn)等級(jí)，并將其作為評(píng)價(jià)單元中心點(diǎn)的風(fēng)險(xiǎn)水平。對(duì)于評(píng)價(jià)地塊邊緣不足10km×10km的網(wǎng)格評(píng)價(jià)樣地，則按照實(shí)際情況計(jì)算各地類所占的百分比。生態(tài)風(fēng)險(xiǎn)如式（1）所示：

ERI=（1）

其中，ERI表示評(píng)價(jià)單元的生態(tài)風(fēng)險(xiǎn)指數(shù);Ai表示評(píng)價(jià)單元內(nèi)第i類土地利用非類型所占的面積;Wi表示第i類土地利用類型所反映的生態(tài)風(fēng)險(xiǎn)的參數(shù);A表示第i個(gè)評(píng)價(jià)單元的總面積。本研究權(quán)重的確定參考胡金龍的研究，耕地、林地、草地、水域、建設(shè)用地、未利用地的權(quán)重分別為0.1916、0.0427、0.0726、0.1425、0.3934、0.1572。

3. 地統(tǒng)計(jì)法

基于上述計(jì)算出的各個(gè)評(píng)價(jià)單元的生態(tài)風(fēng)險(xiǎn)指數(shù)以及生態(tài)風(fēng)險(xiǎn)指數(shù)的空間特性，采用半方差的分析方法對(duì)評(píng)價(jià)單元的生態(tài)風(fēng)險(xiǎn)指數(shù)進(jìn)行插值，最終得到整個(gè)研究區(qū)域生態(tài)風(fēng)險(xiǎn)分布圖，如式（2）所示。

y（h）=[Z（xi+h）-Z（xi）]2（2）

其中，r（h）表示辦方差：Z（xi+h）與Z（xi）分別表述系統(tǒng)在xi+h和xi處的風(fēng)險(xiǎn)值;n（h）表示評(píng)價(jià)單元間距為h樣本對(duì)數(shù)總和：h為評(píng)價(jià)單元間距。運(yùn)用Arcgis中空間模型對(duì)半方差函數(shù)進(jìn)行擬合，再此基礎(chǔ)上對(duì)皖江城市帶風(fēng)險(xiǎn)值進(jìn)行插值。

4. 空間自相關(guān)

空間自相關(guān)分析主要是為了確定變量在空間上的依賴性，其主要思想源自地理學(xué)第一定律。本研究利用morans I指數(shù)來(lái)測(cè)度所有生態(tài)單元的生態(tài)風(fēng)險(xiǎn)值在整體上的自相關(guān)，此外，采用空間關(guān)聯(lián)局域指標(biāo)（LISA）檢測(cè)觀測(cè)值的高值或低值的局部空間集聚。

三、結(jié)果分析

（一）土地利用變化分析

土地轉(zhuǎn)移矩陣可以描述一定階段各土地利用類型相互轉(zhuǎn)換的方向和數(shù)量。1990～2015年，減少最多的地類為耕地，共減少1359.5km2;增加最多的地類為建設(shè)用地，共增加1430.68km2。在土地利用轉(zhuǎn)移方面，耕地主要轉(zhuǎn)化為建設(shè)用地，轉(zhuǎn)移面積為1335.05km2，轉(zhuǎn)移率為3.32%;其次是水域和林地，分別轉(zhuǎn)移了196.32km2和59.5km2，主要原因是耕地被開發(fā)進(jìn)行養(yǎng)殖業(yè)發(fā)展和“退耕還林”政策的施行。建設(shè)用地主要轉(zhuǎn)化為耕地，轉(zhuǎn)移面積為100.56km2，轉(zhuǎn)移率為2.54%，這與我國(guó)實(shí)行的“增減掛鉤”政策密不可分。其他各地類相互之間轉(zhuǎn)移不明顯，具體情況如表1所示。總體上看，主要表現(xiàn)為建設(shè)用地、水域的增加，耕地、林地、草地、未利用地的減少，這與我國(guó)這些年大力追求經(jīng)濟(jì)發(fā)展的現(xiàn)實(shí)保持一致。

皖江城市帶土地利用類型以耕地為主，歷年占比在52.79%左右。1990～2015年耕地、林地、草地、未利用地面積持續(xù)減少，其中耕地面積減少最多，為1359.50km2，占減少面積的87.2%;其次為林地，減少面積為187.07km2，占減小面積的12%;草地和未利用地減少面積相對(duì)較小，共12.45km2，占減少面積的0.8%。建設(shè)用地與水域面積持續(xù)上升，建設(shè)用地增加1430.68km2，占增加面積的91.77%;水域面積增加128.34km2，占增加面積的8.23%。

從土地利用變化時(shí)空轉(zhuǎn)移來(lái)看，1990～2015年轉(zhuǎn)出面積最多為耕地，主要轉(zhuǎn)換為建設(shè)用地、水域和林地;耕地是建設(shè)用地增加的主要土地類型來(lái)源;未利用地隨著人們利用程度的加深，總體面積在逐漸減少，且主要轉(zhuǎn)換為耕地。

從生態(tài)風(fēng)險(xiǎn)指數(shù)和等級(jí)分布上看，1990年、1995年、2000年、2005年、2010年、2015年歷年生態(tài)風(fēng)險(xiǎn)指數(shù)區(qū)間分別[0.04，0.30]、[0.04，0.32]、[0.04，0.33]、[0.04，0.33]、[0.04，0.34]、[0.04，0.35]，生態(tài)風(fēng)險(xiǎn)指數(shù)最大值由0.30增加到0.35，增幅為16.67%;生態(tài)風(fēng)險(xiǎn)等級(jí)主要呈現(xiàn)出在合肥市和蕪湖市兩個(gè)“主副”中心迅速擴(kuò)散，在低等級(jí)風(fēng)險(xiǎn)區(qū)域基本不變的形態(tài)，這也與歷年生態(tài)風(fēng)險(xiǎn)指數(shù)最小值為0.04不變的現(xiàn)實(shí)保持一致。

從生態(tài)風(fēng)險(xiǎn)等級(jí)變化來(lái)看，1990～2015年，生態(tài)單元風(fēng)險(xiǎn)等級(jí)有增有減，且等級(jí)增加的面積多于等級(jí)減少的面積;等級(jí)增加最大幅度為三個(gè)等級(jí)，等級(jí)減少最大幅度為兩個(gè)等級(jí)，研究區(qū)域整體生態(tài)風(fēng)險(xiǎn)水平在增加。從生態(tài)風(fēng)險(xiǎn)等級(jí)面積轉(zhuǎn)入與轉(zhuǎn)出來(lái)看，較低生態(tài)風(fēng)險(xiǎn)區(qū)與中生態(tài)風(fēng)險(xiǎn)區(qū)是五個(gè)風(fēng)險(xiǎn)區(qū)內(nèi)最活躍的兩個(gè)風(fēng)險(xiǎn)區(qū);高生態(tài)風(fēng)險(xiǎn)區(qū)的轉(zhuǎn)入與轉(zhuǎn)出面積都相對(duì)較小，是五個(gè)風(fēng)險(xiǎn)區(qū)內(nèi)變化最小的風(fēng)險(xiǎn)區(qū);高生態(tài)風(fēng)險(xiǎn)區(qū)面積的轉(zhuǎn)入只來(lái)源中生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)，且其轉(zhuǎn)出面積全部為較高生態(tài)風(fēng)險(xiǎn)區(qū)。

從生態(tài)風(fēng)險(xiǎn)空間分析來(lái)看，歷年morans I值在（0.8，1）之間，且隨著時(shí)間的推移morans I值在不斷增長(zhǎng)，表明生態(tài)風(fēng)險(xiǎn)指數(shù)空間相關(guān)程度明顯，且有逐步增強(qiáng)的趨勢(shì)。從局部自相關(guān)分析來(lái)看，1990～2015年皖江城市帶生態(tài)風(fēng)險(xiǎn)指數(shù)主要呈現(xiàn)高—高和低-低聚集分布。高—高聚集主要分布在合肥市大部分區(qū)域以及蕪湖市、滁州市、馬鞍山市小部分區(qū)域，且有逐漸增多的趨勢(shì);低-低聚集主要分布在安慶市、宣城市以及池州市的山區(qū)地帶，整體上保持穩(wěn)定態(tài)勢(shì)。

參考文獻(xiàn)：

[1]傅伯杰，陳利頂，馬克明.黃土丘陵區(qū)小流域土地利用變化對(duì)生態(tài)環(huán)境的影響——以延安市羊圈溝流域?yàn)槔齕J].地理學(xué)報(bào)，1999，54（03）：241-246.

[2]李春燕，南靈.陜西省土地生態(tài)安全動(dòng)態(tài)評(píng)價(jià)及障礙因子診斷[J].中國(guó)土地科學(xué)，2015，29（04）：72-81.

[3]周利軍，張雪萍，陳設(shè).扎龍自然保護(hù)區(qū)土地利用變化與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].自然災(zāi)害學(xué)報(bào)，2009，18（02）：186-190.

[4]周汝佳，張永戰(zhàn)，何華春.基于土地利用變化的鹽城海岸帶生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 地理研究，2016，35（06）：1017-1028.

[5]朱會(huì)義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論[J].地理學(xué)報(bào)， 2003，58（05）：643-650.

[6]劉曉，蘇維詞，王錚，等.基于 RRM 模型的三峽庫(kù)區(qū)重慶開縣消落區(qū)土地利用生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].環(huán)境科學(xué)學(xué)報(bào)，2012，32（01）：248-256.

[7]謝花林.基于景觀結(jié)構(gòu)和空間統(tǒng)計(jì)學(xué)的區(qū)域生態(tài)風(fēng)險(xiǎn)分析[J].生態(tài)學(xué)報(bào)， 2008，28（10）：5021-5026.

[8]王玲.基于土地利用變化的上海市生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].農(nóng)業(yè)科學(xué)研究，2018，39（01）：10-16.

[9]高永年，高俊鋒，許妍.太湖流域水生態(tài)功能區(qū)土地利用變化的景觀生態(tài)風(fēng)險(xiǎn)效應(yīng)[J].自然資源學(xué)報(bào)，2010，25（07）：1088-1096.

[10]曾勇.區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)以呼和浩特市區(qū)為例[J].生態(tài)學(xué)報(bào)，2010，30（03）： 668-673.

[11]胡金龍，周志翔，騰明君，等.基于土地利用變化的典型喀斯特流域生態(tài)風(fēng)險(xiǎn)評(píng)估[J].應(yīng)用生態(tài)學(xué)報(bào)，2017，28（06）：2003-2012.

[12]趙潔巖，李陽(yáng)兵，邵景安.基于土地利用變化的三峽庫(kù)區(qū)小流域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].自然資源學(xué)報(bào)，2013，28（06）：944-956.

[13]李哈濱，王政權(quán)，王慶成.空間異質(zhì)性定量研究理論與方法[J].應(yīng)用生態(tài)學(xué)報(bào)， 1998，9（06）：651-657.

（作者單位：安徽師范大學(xué)地理與旅游學(xué)院）