貴州草海濕地區(qū)域土地利用生態(tài)風(fēng)險評估

任金銅，李仰征，王志紅，歐陽力劍

（貴州工程應(yīng)用技術(shù)學(xué)院生態(tài)工程學(xué)院，貴州 畢節(jié) 551700）

生態(tài)風(fēng)險具有復(fù)雜性、客觀性、不確定性等特征，生態(tài)風(fēng)險評價是在一種或多種外界因素的影響下，對發(fā)生或正在發(fā)生的不利于生態(tài)系統(tǒng)穩(wěn)定和良性循環(huán)的生態(tài)后果進行評估[1]，生態(tài)風(fēng)險評價已成為宏觀生態(tài)管理的一項重要途徑[2]。土地利用作為全球環(huán)境變化的重要組成部分，已成為影響人類社會可持續(xù)發(fā)展的重大問題，土地利用生態(tài)風(fēng)險被用于描述和評價人為活動、自然災(zāi)害或者環(huán)境污染對生態(tài)系統(tǒng)結(jié)構(gòu)和功能所造成的影響[3]，探索土地利用動態(tài)變化基礎(chǔ)上的生態(tài)風(fēng)險，對區(qū)域土地資源合理利用和實現(xiàn)土地生態(tài)環(huán)境動態(tài)平衡具有重要現(xiàn)實意義[4]。

自20世紀80年代以來，國內(nèi)外學(xué)者針對生態(tài)風(fēng)險評價開展了相關(guān)研究工作，有學(xué)者以土地利用為基礎(chǔ)，通過構(gòu)建生態(tài)風(fēng)險指數(shù)模型，對區(qū)域生態(tài)風(fēng)險進行了綜合評價，如任志遠等對榆林市土地生態(tài)風(fēng)險進行了綜合評價[5]，徐羽等分析了鄱陽湖流域土地利用生態(tài)風(fēng)險格局[6]，林媚珍等對中山市土地利用變化及其生態(tài)風(fēng)險響應(yīng)進行了研究[7]；也有學(xué)者基于多源數(shù)據(jù)和相關(guān)模型開展了生態(tài)風(fēng)險評價，如王永芳基于多源數(shù)據(jù)與DPSIR模型對科爾沁沙地進行了生態(tài)風(fēng)險評估[8]，劉曉、蘇維詞等基于RRM模型對三峽庫區(qū)重慶開縣消落區(qū)土地利用生態(tài)風(fēng)險評價[9]，張曉媛、周啟剛等基于綜合模糊評價對三峽庫區(qū)屏障帶重慶段土地利用生態(tài)風(fēng)險評價[10]，岳東霞、曾建軍等基于PSR和熵權(quán)物元可拓模型對甘南高原生態(tài)環(huán)境風(fēng)險評價研究[11]；近年來，景觀生態(tài)風(fēng)險評價理論與方法被高度關(guān)注，基于景觀的生態(tài)風(fēng)險評價主要應(yīng)用于流域[12-13]、行政區(qū)劃[14-15]，以及城市[16-19]等區(qū)域。

貴州草海位于國家級貧困縣威寧縣城附近，隨著經(jīng)濟社會發(fā)展和城鎮(zhèn)化進程的加快，草海及其周邊受人為活動影響較為嚴重，區(qū)域土地利用發(fā)生了明顯變化[20]?；诖耍狙芯恳?S技術(shù)為基礎(chǔ)，以遙感影像解譯獲取近20年來草海濕地區(qū)域土地利用數(shù)據(jù)，通過構(gòu)建生態(tài)風(fēng)險指數(shù)模型，結(jié)合研究區(qū)數(shù)字高程模型（Digital Elevation Model，DEM）對研究區(qū)土地利用生態(tài)風(fēng)險時空變化特征進行了綜合分析，為草海濕地及其周邊生態(tài)環(huán)境保護及土地利用決策提供參考。

1 研究區(qū)概況及數(shù)據(jù)

1.1 研究區(qū)概況

貴州草海位于云貴高原中部頂端的烏蒙山麓腹地，是典型的喀斯特高原湖泊濕地。作為我國為數(shù)不多的亞熱帶高原濕地生態(tài)系統(tǒng)，草海平均海拔2170m，平均水深約2m，具有日照豐富、冬暖夏涼等獨特的高原氣候特征，為我國特有的黑頸鶴等珍稀鳥類及其他候鳥提供了重要的越冬地和遷徙中轉(zhuǎn)站，在“中國生物多樣性保護行動計劃”中被列為Ⅰ級重要濕地，與滇池、青海湖并成為“中國三大高原湖泊”[20-21]。由于地處喀斯特生態(tài)脆弱區(qū)，草海濕地區(qū)域經(jīng)濟、社會發(fā)展在受到自然因素制約的同時，也嚴重影響著區(qū)域生態(tài)環(huán)境問題。2017年8月中央第七環(huán)保督察組反饋草海環(huán)境問題指出“威寧縣城的發(fā)展與草海保護嚴重沖突，‘城進湖退’等問題突出”[22]。

由于草海保護區(qū)緊鄰?fù)幙h城，隨著經(jīng)濟社會的不斷發(fā)展、周邊城鎮(zhèn)化水平的不斷提高，草海濕地區(qū)域人地矛盾、珍稀物種保護與人類社會發(fā)展問題不可回避。為揭示草海濕地區(qū)域土地利用變化導(dǎo)致的生態(tài)風(fēng)險狀況，本研究以草海國家級自然保護區(qū)為中心，研究包括威寧縣小海鎮(zhèn)黑巖洞至草海保護區(qū)一帶的集水區(qū)，涉及草海鎮(zhèn)、小海鎮(zhèn)和雙龍鄉(xiāng)部分區(qū)域，共計53個行政村[23]。

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

本研究所需數(shù)據(jù)主要包括Landsat系列遙感影像、DEM數(shù)據(jù)以及地形圖等。其中，遙感影像數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺①和美國地質(zhì)勘探局網(wǎng)站②提供的軌道號為129/41的Landsat系列影像。DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺③提供的GDEM V2 30m分辨率數(shù)字高程數(shù)據(jù)，地形圖數(shù)據(jù)由威寧縣國土資源局提供，另外輔以中國科學(xué)資源環(huán)境科學(xué)數(shù)據(jù)中心④提供的1:10萬比例尺土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)[24]。

參考《土地利用現(xiàn)狀分類標準》（GBT 21010-2007），結(jié)合研究區(qū)實際和數(shù)據(jù)源特點，將土地利用類型劃分為：耕地、林地、草地、水域和建設(shè)用地等五種類型，通過監(jiān)督分類和人機交互式判讀的方法對不同時期遙感影像數(shù)據(jù)進行解譯，利用混淆矩陣方法對分類結(jié)果進行評價，將分類結(jié)果總體分類精度控制在85%以上。草海濕地區(qū)域土地利用類型解譯結(jié)果如圖1所示。

圖1 草海濕地區(qū)域土地利用解譯結(jié)果圖

2 研究方法

本研究以草海濕地區(qū)域遙感影像解譯出的土地利用類型數(shù)據(jù)為基礎(chǔ)，借助3S技術(shù)方法，利用生態(tài)風(fēng)險指數(shù)，結(jié)合研究區(qū)地形數(shù)據(jù)進行生態(tài)風(fēng)險時空變化特征評價。

2.1 高程和坡度劃分

利用30m分辨率GDEM V2 DEM數(shù)據(jù)，根據(jù)研究區(qū)高程分布情況（2064-2890m），考慮到草海湖區(qū)正常水位2171.7m，本研究利用ArcGIS重分類工具將研究區(qū)劃分為五個高程段，其中2064-2172m作為一個高程段，其他高程段以200m遞增，分別為：2172-2372 m、2372-2572 m、2572-2772 m、2772-2890 m。

由于研究區(qū)地處喀斯特山區(qū)，為進一步研究不同坡度條件下生態(tài)風(fēng)險變化情況，本研究利用Arc?GIS 3D分析工具對DEM數(shù)據(jù)進行處理，獲取研究區(qū)坡度數(shù)據(jù)。結(jié)合研究區(qū)坡度分布情況，將坡度重分類為0°～3°、3°～8°、8°～15°、15°～25°和＞25°五個坡度段。

2.2 生態(tài)風(fēng)險評價模型

根據(jù)研究區(qū)范圍，利用各地類面積比例構(gòu)建生態(tài)風(fēng)險指數(shù)（ERI）評價模型，運用網(wǎng)格采樣法將研究區(qū)劃分為1km×1km的正方形網(wǎng)格作為評價單元，共劃分得到576個采樣區(qū)，計算每個采樣網(wǎng)格的生態(tài)風(fēng)險，作為樣地中心點的生態(tài)風(fēng)險水平[25]。

其中，ERI為生態(tài)風(fēng)險指數(shù)；i為研究區(qū)各土地利用類型，Ai為采樣區(qū)內(nèi)第i種土地利用類型的面積，Wi為采樣區(qū)內(nèi)第i種土地利用類型所反映的生態(tài)風(fēng)險強度參數(shù)。參照相關(guān)研究成果及研究區(qū)實際情況，并經(jīng)咨詢專家打分法，最終確定研究區(qū)各土地利用類型的生態(tài)風(fēng)險強度參數(shù)Wi值依次為：林地0.0427、耕地0.1916、建設(shè)用地0.3934、水域0.1425、草地0.0726、未利用地0.1572[26]。

3 草海濕地區(qū)域生態(tài)風(fēng)險變化特征分析

3.1 生態(tài)風(fēng)險隨時間變化特征

利用ArcGIS統(tǒng)計分析工具結(jié)合公式（1），對1995-2015年間草海濕地區(qū)域576個樣區(qū)的土地利用生態(tài)風(fēng)險指數(shù)進行計算，得出研究區(qū)土地利用生態(tài)風(fēng)險指數(shù)統(tǒng)計表（表1），進一步分析可知：（1）1995年草海濕地區(qū)域生態(tài)風(fēng)險值在0.0423～0.369，均值為0.103，標準差為0.044；2005年生態(tài)風(fēng)險值在0.043～0.373，均值為0.107，標準差為0.045；2015年的生態(tài)風(fēng)險值0.043～0.373，均值為0.110，標準差為0.049。（2）近20年，草海濕地區(qū)域的土地利用生態(tài)風(fēng)險峰值由0.369增加到0.373，均值由0.103增加到0.110，標準差由0.044持續(xù)增加到0.049。由此可知，草海濕地區(qū)域土地利用生態(tài)風(fēng)險指數(shù)有向高風(fēng)險變化趨勢。

表1 研究區(qū)土地利用生態(tài)風(fēng)險指數(shù)統(tǒng)計表

按照生態(tài)風(fēng)險值在空間出現(xiàn)的概率，參考相關(guān)研究、結(jié)合研究區(qū)實際情況，利用GIS標準分類法中的幾何間隔，依次將研究區(qū)生態(tài)風(fēng)險劃分為低生態(tài)風(fēng)險、較低生態(tài)風(fēng)險、中等生態(tài)風(fēng)險、較高生態(tài)風(fēng)險和高生態(tài)風(fēng)險五個等級[25]（圖2）。

圖2 草海濕地區(qū)域土地利用生態(tài)風(fēng)險分布

通過對研究區(qū)生態(tài)風(fēng)險等級結(jié)果進行統(tǒng)計（表2），結(jié)合土地利用類型圖（圖1）分析可知：（1）近20年來，草海濕地區(qū)域以中等和較低生態(tài)風(fēng)險為主，二者總面積所占比例均超過70%，其中中等生態(tài)風(fēng)險占比均超過40%，并呈現(xiàn)持續(xù)下降趨勢；中等和較低生態(tài)風(fēng)險區(qū)域以耕地和草地為主，由于耕地面積的增加導(dǎo)致中等生態(tài)風(fēng)險區(qū)域面積減少。（2）高生態(tài)風(fēng)險區(qū)面積占比較小，均低于1%，但總體呈現(xiàn)上升趨勢；隨著經(jīng)濟社會發(fā)展及城鎮(zhèn)化進程的加快，研究區(qū)內(nèi)建設(shè)用地持續(xù)增加，由于建設(shè)用地生態(tài)風(fēng)險強度參數(shù)最大，導(dǎo)致高生態(tài)風(fēng)險面積持續(xù)增加。（3）較高生態(tài)風(fēng)險區(qū)先減后增，總體呈現(xiàn)下降趨勢；低生態(tài)風(fēng)險區(qū)先增后減，總體呈現(xiàn)增加趨勢。

表2 研究區(qū)生態(tài)風(fēng)險面積統(tǒng)計

3.2 生態(tài)風(fēng)險的空間變化分析

從空間變化的角度，結(jié)合研究區(qū)土地利用生態(tài)風(fēng)險分布（圖2），對草海濕地區(qū)域土地利用生態(tài)風(fēng)險變化情況進行分析可知：（1）近20年來，草海濕地區(qū)域生態(tài)風(fēng)險的空間分布特征較為明顯，其中，威寧縣城附近區(qū)域一直處于高生態(tài)風(fēng)險，草海湖區(qū)和楊家灣水庫區(qū)域處于較高生態(tài)風(fēng)險，中部區(qū)域以中等和較低生態(tài)風(fēng)險區(qū)為主，低生態(tài)風(fēng)險區(qū)主要位于草海濕地區(qū)域周邊以林地為主的山區(qū)。（2）草海濕地區(qū)域高生態(tài)風(fēng)險區(qū)域主要集中在威寧縣城附近，且隨著建設(shè)用地面積持續(xù)增加，高生態(tài)風(fēng)險區(qū)面積持續(xù)增加；較高生態(tài)風(fēng)險區(qū)由連片區(qū)域不斷收縮，導(dǎo)致其空間分布逐漸分散；中等生態(tài)風(fēng)險區(qū)位于較高生態(tài)風(fēng)險區(qū)外圍，以耕地類型為主；低和較低生態(tài)風(fēng)險區(qū)以林地和草地為主，面積先增后減，但隨著退耕還林以及草海周邊水源地保護等政策的實施，使得低和較低生態(tài)風(fēng)險區(qū)面積總體呈增加趨勢。

利用ArcGIS分別將草海濕地區(qū)域1995-2005年和2005-2015年生態(tài)風(fēng)險進行疊加分析，借助EX?CEL數(shù)據(jù)透視表功能，獲取研究區(qū)生態(tài)風(fēng)險轉(zhuǎn)移情況（表3、表4）。經(jīng)分析可知，近20年來研究區(qū)生態(tài)風(fēng)險在空間上發(fā)生不同程度的轉(zhuǎn)換，其中：（1）1995-2005年間，低生態(tài)風(fēng)險主要向較低生態(tài)風(fēng)險轉(zhuǎn)換，轉(zhuǎn)換面積為1615.568hm2；較低生態(tài)風(fēng)險主要向低和中等生態(tài)風(fēng)險轉(zhuǎn)換為主，轉(zhuǎn)換面積分別為1615.568 hm2、1050.162 hm2；中等生態(tài)風(fēng)險主要向較低生態(tài)風(fēng)險轉(zhuǎn)換，轉(zhuǎn)換面積為5014.542 hm2，其中750.243hm2向較高生態(tài)風(fēng)險轉(zhuǎn)換；較高生態(tài)風(fēng)險主要向中等生態(tài)風(fēng)險轉(zhuǎn)換，轉(zhuǎn)換面積為2860.010 hm2，另有86.886 hm2轉(zhuǎn)為了較高生態(tài)風(fēng)險；高生態(tài)風(fēng)險未發(fā)生明顯轉(zhuǎn)移情況。（2）2005-2015年間，低生態(tài)風(fēng)險主要向較低和中等生態(tài)風(fēng)險轉(zhuǎn)移為主，轉(zhuǎn)移面積分別為838.984 hm2、133.799 hm2；較低生態(tài)風(fēng)險主要向中等生態(tài)風(fēng)險轉(zhuǎn)移為主，轉(zhuǎn)移面積有2144.855 hm2，另有190.061 hm2轉(zhuǎn)為低生態(tài)風(fēng)險，7.590 hm2轉(zhuǎn)為較高生態(tài)風(fēng)險；中等生態(tài)風(fēng)險主要向較高生態(tài)風(fēng)險轉(zhuǎn)移為主，2423.245 hm2轉(zhuǎn)為較高生態(tài)風(fēng)險，555.216 hm2轉(zhuǎn)為較低生態(tài)風(fēng)險；較高生態(tài)風(fēng)險中260.447 hm2轉(zhuǎn)為中等生態(tài)風(fēng)險，24.817 hm2轉(zhuǎn)為高生態(tài)風(fēng)險；高生態(tài)風(fēng)險轉(zhuǎn)移依然不明顯，只有0.635 hm2轉(zhuǎn)為較高生態(tài)風(fēng)險。

表3 1995-2005年間研究區(qū)生態(tài)風(fēng)險面積轉(zhuǎn)移矩陣（單位：hm2）

表4 2005-2015年間研究區(qū)生態(tài)風(fēng)險面積轉(zhuǎn)移矩陣（單位：hm2）

3.3 地形特征對生態(tài)風(fēng)險變化影響分析

3.3.1 不同高程帶生態(tài)風(fēng)險變化特征

將研究區(qū)重分類后的高程數(shù)據(jù)與三個時期生態(tài)風(fēng)險分布進行疊加分析，統(tǒng)計不同時期各高程帶上不同生態(tài)風(fēng)險分布情況（表5），結(jié)合公式（1）計算得出研究區(qū)不同高程帶的生態(tài)風(fēng)險指數(shù)（表6）。經(jīng)統(tǒng)計分析可知，近20年來，草海濕地區(qū)域生態(tài)風(fēng)險分布與研究區(qū)不同高程帶具有明顯的相關(guān)性，其中：（1）低生態(tài)風(fēng)險和較低生態(tài)風(fēng)險區(qū)主要位于2172-2372m和2372-2572m高程帶，中等、較高和高生態(tài)風(fēng)險主要位于2172-2372m高程帶，大于2372m高程帶不存在高生態(tài)風(fēng)險，大于2572m高程不存在較高生態(tài)風(fēng)險；結(jié)合圖1進一步分析可知，研究區(qū)位于典型喀斯特山區(qū)，其土地利用類型與地形高程密切相關(guān)，其中，2172-2372m高程帶土地利用類型以耕地和草地為主，且建設(shè)用地主要位于該區(qū)域，而大于2572m高程帶的土地利用以林地為主。

表5 不同高程帶生態(tài)風(fēng)險面積分布（面積單位：hm2）

（2）從不同高程帶的生態(tài)風(fēng)險指數(shù)變化情況來看，20年間草海濕地區(qū)域不同高程帶的生態(tài)風(fēng)險指數(shù)均出現(xiàn)不同程度的增加，其中，高程2372m以下的生態(tài)風(fēng)險指數(shù)增加了0.008，高程2572-2772m的生態(tài)風(fēng)險指數(shù)未發(fā)生明顯變化，大于2772m的生態(tài)風(fēng)險指數(shù)減少了0.001。

表6 不同高程帶生態(tài)風(fēng)險指數(shù)

3.3.2 不同坡度段生態(tài)風(fēng)險變化情況

通過對不同時期草海濕地區(qū)域生態(tài)風(fēng)險與研究區(qū)坡度分類結(jié)果進行疊加分析，得出不同坡度帶的各級生態(tài)風(fēng)險面積分布（表7），結(jié)合公式（1）進一步計算得出不同坡度帶的生態(tài)風(fēng)險指數(shù)（表8）。研究發(fā)現(xiàn)：（1）草海濕地區(qū)域0～3°坡度帶主要分布有草海湖區(qū)和楊家灣水庫，以水域為主，生態(tài)風(fēng)險以中等和較高生態(tài)風(fēng)險為主，該區(qū)域易受人類活動影響，20年間中等生態(tài)風(fēng)險區(qū)域面積先增后減，總體呈增加趨勢，較高生態(tài)風(fēng)險區(qū)面積先減后增，總體呈下降趨勢。（2）高和較高生態(tài)風(fēng)險區(qū)主要位于3～8°和8～15°兩個坡度段，隨著坡度的增加，較高和高生態(tài)風(fēng)險面積持續(xù)下降，該區(qū)域林地分布較低。（3）8～15°坡度帶是各級別生態(tài)風(fēng)險區(qū)域主要分布帶，在該坡度帶，高生態(tài)風(fēng)險區(qū)域面積持續(xù)增加，低生態(tài)風(fēng)險區(qū)則持續(xù)減少，較高生態(tài)風(fēng)險區(qū)面積先減少后增加，總體呈現(xiàn)減少趨勢；較低和中等生態(tài)風(fēng)險區(qū)面積先增加后減少，但總體呈現(xiàn)增加趨勢。（4）低、較低和中等生態(tài)風(fēng)險區(qū)在15～25°坡度帶的分布僅次于8～15°坡度帶，在15～25°坡度帶則是以中等生態(tài)風(fēng)險為主，且高生態(tài)風(fēng)險在該區(qū)域內(nèi)面積持續(xù)增加。（5）大于25°的坡度以較低和中等生態(tài)風(fēng)險為主，低生態(tài)風(fēng)險在該坡度帶面積持續(xù)減少，高生態(tài)風(fēng)險在該坡度帶占比最小，但面積呈增加趨勢。

表7 不同坡度帶生態(tài)風(fēng)險面積分布（面積單位：hm2）

從生態(tài)風(fēng)險指數(shù)變化情況（表8）來看，近20年來草海濕地區(qū)域生態(tài)風(fēng)險指數(shù)均發(fā)生了不同程度的增加，其中，0～3°坡度帶生態(tài)風(fēng)險指數(shù)增加0.008，3～8°和8～15°兩個坡度段增加0.007，15～25°坡度帶則增加0.006，大于25°的坡度帶增加0.005。

表8 不同坡度帶的生態(tài)風(fēng)險指數(shù)變化

4 結(jié)論與討論

本研究以3S技術(shù)為基礎(chǔ)，以遙感影像解譯獲取近20年來研究區(qū)土地利用數(shù)據(jù)，參考前人研究構(gòu)建了生態(tài)風(fēng)險指數(shù)模型，針對研究區(qū)地形特點，結(jié)合研究區(qū)數(shù)字高程模型對研究區(qū)土地利用生態(tài)風(fēng)險時空變化特征進行了綜合分析，得出以下結(jié)論：

（1）近20年來，草海濕地區(qū)域土地利用生態(tài)風(fēng)險指數(shù)有向高風(fēng)險變化趨勢。研究區(qū)內(nèi)以中等和較低生態(tài)風(fēng)險為主，并呈現(xiàn)持續(xù)下降趨勢；高生態(tài)風(fēng)險區(qū)面積占比較小，但總體呈現(xiàn)上升趨勢；較高生態(tài)風(fēng)險區(qū)總體呈現(xiàn)下降趨勢，低生態(tài)風(fēng)險區(qū)總體呈現(xiàn)增加趨勢。

（2）研究區(qū)土地利用生態(tài)風(fēng)險的空間分布特征較為明顯，縣城附近區(qū)域一直處于高生態(tài)風(fēng)險，草海湖區(qū)和楊家灣水庫區(qū)域處于較高生態(tài)風(fēng)險，中部區(qū)域以中等和較低生態(tài)風(fēng)險區(qū)為主，低生態(tài)風(fēng)險區(qū)主要位于草海濕地區(qū)域周邊以林地為主的山區(qū)。

（3）隨著經(jīng)濟社會發(fā)展和城鎮(zhèn)化進程的加快，研究區(qū)土地利用生態(tài)風(fēng)險在空間上發(fā)生不同程度的轉(zhuǎn)換。低生態(tài)風(fēng)險主要向較低生態(tài)風(fēng)險轉(zhuǎn)換，較低生態(tài)風(fēng)險主要向低和中等生態(tài)風(fēng)險轉(zhuǎn)換為主，中等生態(tài)風(fēng)險主要向較低生態(tài)風(fēng)險轉(zhuǎn)換，較高生態(tài)風(fēng)險主要向中等生態(tài)風(fēng)險轉(zhuǎn)換，高生態(tài)風(fēng)險未發(fā)生明顯轉(zhuǎn)移情況。

（4）近20年來，草海濕地區(qū)域土地利用生態(tài)風(fēng)險分布與研究區(qū)不同高程帶和坡度帶具有明顯的相關(guān)性，高程越大土地利用生態(tài)風(fēng)險等級越低，坡度越小土地利用生態(tài)風(fēng)險等級越高。

土地利用對生物多樣性、生態(tài)系統(tǒng)服務(wù)功能及生態(tài)系統(tǒng)資源環(huán)境的承載力有著直接的影響，是近年來全球變化研究的熱點問題之一。基于土地利用的生態(tài)風(fēng)險評估，對于區(qū)域生態(tài)環(huán)境保護及土地利用決策具有重要現(xiàn)實意義。通過對草海濕地區(qū)域近20年來的土地利用生態(tài)風(fēng)險的時空變化特征進行研究可知，人類活動對區(qū)域土地利用生態(tài)風(fēng)險變化起到主導(dǎo)作用，隨著經(jīng)濟社會發(fā)展及城鎮(zhèn)化進程的加快，建設(shè)用地面積急劇增加是導(dǎo)致生態(tài)風(fēng)險向高生態(tài)風(fēng)險變化的主要原因；隨著退耕還林以及草海周邊水源地保護等政策的實施，使得研究區(qū)內(nèi)低和較低生態(tài)風(fēng)險區(qū)面積總體呈增加趨勢。因此，制定科學(xué)合理的土地利用政策對區(qū)域生態(tài)環(huán)境保護具有重要意義。

注釋:

①遙感影像數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺，http://www.gscloud.cn/。

②遙感影像數(shù)據(jù)來源于美國地質(zhì)勘探局網(wǎng)站，http://glovis.usgs.gov/。

③DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺，http://www.gscloud.cn/.地形圖數(shù)據(jù)于1990年12月航攝，1992年4月調(diào)繪。

④土地利用現(xiàn)狀數(shù)據(jù)來源于中國科學(xué)資源環(huán)境科學(xué)數(shù)據(jù)中心，http://www.resdc.cn/。