基于GIS和RS的砒砂巖區生態環境質量綜合評價

姚俊娜, 秦 奮

(河南大學 黃河中下游數字地理技術實驗室, 河南 開封 475000)

基于GIS和RS的砒砂巖區生態環境質量綜合評價

姚俊娜, 秦 奮

(河南大學 黃河中下游數字地理技術實驗室, 河南 開封 475000)

以Landsat TM影像數據、數字高程模型(DEM)、土壤類型數據、基礎地圖和統計資料為數據源，以GIS和RS為技術手段，采用主成分分析法評價模型，對砒砂巖區2000年和2010年的生態環境質量狀況進行了綜合評價。結果表明：砒砂巖區2000年和2010年的生態環境質量總體情況一般，但與2000年的生態環境質量情況相比，2010年有明顯的改善。GIS和RS相結合的手段是進行大范圍、條件復雜和數據難以獲取地區生態環境質量綜合評價的有效手段。

地理信息系統; 遙感; 生態環境質量; 綜合評價; 砒砂巖

砒砂巖分布區是黃河粗泥沙的主要來源區，被稱為“世界上水土流失之最”和“環境癌癥”。近年來，為了改善該區域的生態環境質量，國家和地方投入了大量的人力和物力，并且高標準實施了“國家水土保持重點建設工程”、“黃河上中游水土保持重點防治工程”、“黃土高原地區水土保持淤地壩試點工程”、“晉陜蒙砒砂巖去沙棘生態工程”等一大批生態建設重點工程，該區域的生態環境質量得到了明顯的改善。但是砒砂巖地區生態環境依舊非常脆弱，生態環境問題依舊異常突出，特別是水資源短缺、水土流失、土地荒漠化等生態環境問題逐步加劇，生態環境治理工作面臨嚴重的挑戰。

采用GIS和RS技術進行生態環境質量評價，國內外許多學者進行了研究和嘗試。早在1994年加拿大就建立了生態監測與評估網絡，并利用該網絡獲得的長期監測數據對生態環境進行了評估[1];Basso等基于GIS和RS對意大利南部Agri流域的環境脆弱性進行了評價[2]；國內學者中鄭新奇[3]、王宏偉[4]、劉建[5]等都曾基于GIS和RS技術分別對區域、流域和縣級行政區的生態環境質量進行了綜合評價。但是他們的生態環境質量評價工作，以行政區和流域邊界為研究邊界的居多，涉及到像砒砂巖地區這種跨流域、跨行政區的研究相對較少。本文采用GIS和RS技術進行砒砂巖區生態環境質量綜合評價，探索生態環境治理工作取得的成果和不足，為以后的工作提供意見和建議。

1 研究區概況

砒砂巖分布區位于東經108°45′—111°31′、北緯38°10′—40°10′的晉、陜、蒙三省交界處，集中分布在內蒙古自治區鄂爾多斯市的東勝區、準格爾旗、伊金霍洛旗、達拉特旗、杭錦旗以及陜西省的神木、府谷兩縣，還有山西省的河曲、保德兩縣，在內蒙古的清水河縣也有零星的分布，分布區總面積約1.67萬km2[6]。砒砂巖分布區人口密度相對較低，截至2010年，研究區內人口約69萬。砒砂巖分布區地處內陸，屬典型的干旱、半干旱大陸性氣候，年平均氣溫6～9℃，降水保證率低，并且年際變化較大，汛期主要集中在6—9月，多年平均降水量為306～453 mm。研究區地勢西高東低，海拔高度73～1 625 m，境內大部分地區溝網縱橫密布，地表被切割呈支離破碎狀，地貌以風沙地貌和丘陵溝壑地貌為主。地表土壤主要是黃土、礫土等三種，植被種類較多，以耐旱耐瘠、蓄水保土的植被為主，如：沙棘、檸條、油松、油蒿等。

2 數據來源和評價指標體系的構建

2.1 數據來源

砒砂巖分布幅員遼闊、地質地貌情況復雜、溝道縱橫、河網密布、植被種類多樣，生態環境質量評價工作復雜而艱難。本文考慮到研究區的具體情況，在進行生態環境質量綜合評價時以多元遙感數據、地形圖及其它專題數據作為數據源，采用遙感和地理信息技術進行數據處理，構建研究區環境質量評價數據庫。表1是本研究所用到的數據及其來源。

表1 本研究所用數據及其數據來源

2.2 指標體系

生態環境系統是一個龐大而復雜的綜合系統，包含了自然、社會等很多因素。生態環境質量綜合評價指標體系是生態環境質量評價的基本尺度和衡量標準，是生態環境綜合評價的根本條件和理論基礎。生態環境質量評價指標的選擇應遵循科學性、代表性、全面性和可操作性的原則[7]。本文在參考周小成[8]和朱遠輝等[9]研究的基礎上，建立包括水熱氣象、地形地貌、土地覆蓋和土壤侵蝕在內的4個一級指標，其中又容納了歸一化水汽指數(NDMI)、歸一化差異水體指數NDWI、積溫指數、降雨量指數、高程指數、坡度指數、坡向指數、植被指數NDVI、土壤亮度指數、土地利用和土壤侵蝕強度在內的11個二級指標，見表2。

表2 砒砂巖區生態環境質量評價指標以及計算方法

2.2.1 水熱氣象因子 影響生態環境狀況的水熱氣象因子是多方面的。本研究充分利用遙感數據快速而準確提取水熱氣象因子的優勢，從遙感數據中提取歸一化水汽指數和歸一化差異水體指數，方法見表2。積溫和降雨量指數提取時采用的數據源為研究區周圍57個氣象站點的觀測數據。在積溫指數的計算過程中，首先以57個氣象站2000年和2010年的積溫數據為分析樣本，以積溫為因變量，以氣象站的經度、緯度和海拔高度為自變量進行多元線性回歸，分別得到2000年和2010年積溫與海拔高度和經緯度的回歸方程，然后再進行插值和殘差的空間內插等工作，最終獲得研究區的積溫數據。降雨量指數的計算參考孟慶香等[10]的研究成果，直接以氣象站的觀測數據為基礎，采用克里金插值方法進行空間內插而得到。

2.2.2 地形地貌因子 地形地貌影響著土壤和植被的形成和發育過程，也是造成小氣候差異的重要原因。本文選取地形中最具有代表性的高程、坡度和坡向作為地形地貌因子指標。高程指數、坡度指數和坡向指數因子的提取直接以研究區經過處理的DEM數據為基礎，經過標準化處理而獲得，方法見表2。

2.2.3 土地覆蓋因子 土地覆蓋情況的差異是造成環境狀況差異的重要因素，本文選取歸一化植被指數、土壤亮度指數和土地利用情況三個能夠反映研究區植被覆蓋情況、土地裸化程度和土地利用情況的指標作為土地覆蓋因子指標。提取到的歸一化植被指數、土壤亮度指數和土地利用數據年際間的差異與采用的遙感影像的月份有著直接的聯系，為了消除影像時間差帶來的影響，本研究所采用的影像均是7月份。歸一化植被指數和土壤亮度指數的提取方法見表2。考慮到研究區的實際情況，本研究將研究區土地利用類型劃分為6大類，分別為：耕地、林地、草地、建設用地、水域和裸地。在土地利用分類體系的基礎上，采用目視解譯和手動糾正的方法獲得研究區的土地利用數據。

2.2.4 土壤侵蝕因子 土壤侵蝕是全球性的主要問題之一，不但導致土地退化、土地生產力降低，影響農業生產和食物安全，而且隨徑流泥沙運移的污染物質對異地生態、環境、人類生存和社會經濟發展帶來嚴重影響。砒砂巖分布區存在著嚴重的水土流失問題。目前土壤侵蝕模型方面的研究成果很多，本文結合研究區實際情況及其可操作性，參考齊清等[11]的研究，選用RUSLE計算土壤侵蝕量，見式(1)，得到研究區2000年和2010年的土壤侵蝕量。

A=RKLSCP

(1)

式中：A——土壤侵蝕量；R——降雨侵蝕因子；K——土壤可侵蝕性因子；L——坡長因子；S——坡度因子；C——作物覆蓋與管理因子；P——水土保持措施因子。

3 砒砂巖區生態環境質量綜合評價

3.1 評價模型的選擇

目前國內外應用的生態環境質量評價方法主要有層次分析法、人工神經網絡評價法、灰色綜合評估法、物元分析法、模糊綜合評價法、主成分分析法等[12]，本研究選取主成分分析法建立評價模型。主成分分析法是一種基于線性變換，找到一組新的替換指標,從而實現降維和評估的方法[13]。這種分析方法在保證原始數據信息損失最小的前提下，經過線性變換和舍棄部分信息，以少數的綜合變量取代原有的多維變量，這樣既抓到了重點，又簡化了工作，并且在整個過程中不再需要專家打分，使評價的結果更加客觀和真實。

3.2 數據的標準化

進行生態環境質量綜合評價所采用的指標數據類型復雜、來源不一致并且量綱也不統一，這就使得生態環境質量評價的各項指標之間不具有直接的可比性，不能直接進行生態環境質量評價。因此首先應對評價指標進行標準化處理，從而消除量綱以及量綱單位的不同帶來的不可公度性[14]。本文所用到的標準化公式見式(2)。

(2)

式中：Ai——第i個單項指標的標準化值；Xi——第i個單項指標的原始值；Xmin——研究區域內該指標的最小值；Xmax——研究區域內該指標的最大值。

3.3基于主成分分析的砒砂巖區生態環境質量綜合評價

利用標準化后的各評價指標數據在ArcGIS中執行空間主成分分析，提取出累計貢獻率達到85%的主成分因子，具體過程如下：

(1) 將標準化后的歸一化水汽指數NDMI、降雨量指數、歸一化差異水體指數NDWI、積溫指數進行主成分分析，生成水熱氣象因子；將標準化后的高程指數、坡度指數、坡向指數進行主成分分析，生成地形地貌因子；將標準化后的植被指數NDVI、土壤亮度指數NDSI和土地利用數據進行主成分分析，生成土地覆蓋因子。

(2) 在上述分析的基礎上，將水熱氣象因子、地形地貌因子、土地覆蓋因子、土壤侵蝕因子進行第二次主成分分析，并且按照公式(3)計算出砒砂巖區生態環境質量綜合指數。

(3)

式中：E——生態環境質量綜合指數；E4CI——由三項地形指標進行主成分分析提取的主因子；Wci——對應主因子的貢獻率；r——提取的主因子的數量。

由主成分分析的結果可知(表3)，在對水熱氣象因子、地形地貌因子以及土地覆蓋因子的相關指標進行第一次主成分分析時，前兩個主成分的累積貢獻率均達到92%以上，因此計算水熱氣象因子、土地覆蓋因子和地形地貌因子時均選取前兩個主成分，就可以保證保存了大部分的數據信息，滿足分析的需要，由于篇幅有限，第一次主成分分析的特征向量在此不再列出。在對水熱氣象因子、地形地貌因子、土地覆蓋因子和土壤侵蝕因子進行第二次主成分分析后可知，前三個主成分的累計貢獻率已經達到92%以上，因此在進行生態環境質量指數計算時，選取前三個主成分，就可以達到較高的可信度，滿足分析的需要。有第二次主成分分析特征向量值(表4)的大小可知，土地覆蓋因子、土壤侵蝕因子和地形地貌因子的特征向量值比較大，說明砒砂巖分布區生態環境質量受土地覆蓋、地形地貌和土壤侵蝕的影響比較大，受水熱氣象的影響相對較小。

表3 主成分分析的特征值、貢獻率和累計貢獻率

表4 各主成分特征向量

4 生態環境質量評價結果分級和原因分析

4.1生態環境質量指數分級

生態環境質量綜合指數代表環境質量的狀況，為了便于比較分析，將環境質量指數進行分級處理，把生態環境綜合評價結果劃分為良好(>70)、較好(60～70)、一般(50～60)、較差(40～50)、惡劣(≤40)共5級，見表5。不同等級綜合指數的分級空間分布特征，體現了生態環境狀況的區域性差異，經過以上分析處理得到的砒砂巖區2000年和2010年生態環境質量分級圖分別見附圖10，分級結果面積統計見表6。

表5 砒砂巖區生態環境質量評價結果分級標準

表6 砒砂巖區生態環境質量評價分級結果統計對比表

4.2 生態環境質量評價結果分析

4.2.1 生態環境質量變化信息的提取 只有合理地提取生態環境的變化信息才能更好地研究本研究區生態環境在時間上和空間上的變化情況。本文參考土地利用轉移矩陣的計算方法[15]，提取出生態環境質量等級變化圖，并將其分為生態環境質量等級不變、提高和降低三類，等級提高的面積2 330.89 m2，等級降低的面積1 636.21 m2，等級未變化的面積13 024.80 m2，等級提高和等級降低的面積只差694.69 m2。

4.2.2 結果分析 將兩個時期生態環境質量提高和下降的面積比較，如提高的面積大于下降的面積，則認為后一時期的生態環境質量總體優于前一時期的生態環境質量。整體上來看，砒砂巖分布區的生態環境質量總體一般，2000年和2010年生態環境質量良好所占的比例不足5%，而生態環境質量一般以下的區域所占據的比例超過40%。

從時間變化上來看，2010年研究區中生態環境惡劣的面積占研究區總面積的0.03%，相比2000年的0.3%，有了相應的提高；生態環境較差的面積2010年所占的比例為2.28%，相對2000年的5.12%有了明顯的降低；2010年生態環境質量一般的面積占研究區總面積的39.39%，與2000年相比基本保持不變；研究區域2010年生態環境質量較好的區域面積占總面積的56.05%，與2000年的50.72%相比有了較大幅度的提高。

從空間分布來看，一方面，砒砂巖分布區生態環境質量的共同特點是中部好于東、西部，北部好于南部；另一方面，砒砂巖分布區的生態環境質量等級分布不平衡，具有明顯的區域分布特點，生態環境質量較好的地區分布在海拔相對居中的鄂爾多斯市的準格爾旗和伊金霍洛旗，以及達拉特旗的東南部地區，生態環境質量狀況較差和惡劣的區域2000年主要分布在南部和東北部，而在2010年生態環境狀況較差和惡劣的區域主要分布在西北部和東部。2010年和2000年相比研究區的東南部和西北部的生態環境狀況有著明顯的好轉，這可能與當地的環保措施有著直接的關系，也反映了砒砂巖區生態環境建設工作取得了一定的成效。

5 結 論

(1) 砒砂巖分布區2000—2010年生態環境質量狀況整體較差，生態環境質量等級以一般和較差為主，從空間分布上來看，砒砂巖分布區生態環境質量中部好于東、西兩端，北部好于南部。

(2) 采用主成分分析法作為砒砂巖區生態環境質量綜合評價模型，避免了人為賦權重值帶來的誤差，并且保留了絕大部分的信息，分析過程更加簡便、快捷和準確。

(3) GIS與RS技術相結合的技術手段，能夠使評價數據易于獲取，并且信息量豐富，非常適合大范圍、條件復雜的區域生態環境質量評價工作，并且評價的結果能夠客觀準確。

[1] 陳彩虹,沈翠新.現代信息技術與生態環境質量評價[J].經濟林研究,2003,21(4):131-133.

[2] Basso F, Bove E, Dumontet S, et al. Evaluating environmental sensitivity at the basin scale through the use of geographic information systems and remotely sensed data: an example covering the Agri basin(Southern Italy)[J]. Catena,2000,40(1):19-35.

[3] 鄭新奇,王愛萍.基于RS與GIS的區域生態環境質量綜合評價研究:以山東省為例[J].環境科學學報,2000,20(4):89-493.

[4] 王宏偉,張小雷,喬木,等.基于GIS的伊犁河流域生態環境質量評價與動態分析[J].干旱區地理,2008,31(2):215-221.

[5] 劉建,何宗,張治清.基于GIS的蒲江縣生態環境綜合評價研究[J].水土保持研究,2012,19(5):130-132.

[6] 王愿昌,吳永紅,寇權,等.砒砂巖分布范圍界定與類型區劃分[J].中國水土保持科學,2007,5(1):14-18.

[7] 林積泉,王伯鐸,馬俊杰,等.小流域治理環境質量綜合評價指標體系研究[J].水土保持研究,2005,12(1):69-71.

[8] 周小成,汪小欽,江洪,等.九龍江流域生態環境質量遙感評價與分析[J].地球信息科學學報,2009,11(2):231-236.

[9] 朱遠輝,劉凱,艾彬,等.新豐江流域生態環境質量遙感評價[J].熱帶地理,2013,33(2):133-140.

[10] 孟慶香,劉國彬,楊勤科.基于GIS的黃土高原氣象要素空間插值方法[J].水土保持研究,2010,17(1):10-14.

[11] 齊清,王天明,寇曉軍,等.基于GIS的黃土高原小流域土壤侵蝕定量評價[J].水土保持研究,2009,16(3):1-5,37.

[12] 鄧春光,張曉麗,崔文超.基于“3S”技術的生態環境質量評價研究進展[J].林業調查規劃,2007,32(3):14-17.

[13] 馮利華.環境質量的主成分分析[J].數學的實踐與認識,2003,33(8):32-35.

[14] 麻素挺,湯潔,林年豐.基于GIS與RS多源空間信息的吉林西部生態環境綜合評價[J].資源科學,2004,26(4):140-145.

[15] 劉瑞,朱道林.基于轉移矩陣的土地利用變化信息挖掘方法探討[J].資源科學,2010,32(8):1544-1550.

ComprehensiveAssessmentonEco-environmentalQualityoftheAreaofSandstonBasedonRSandGIS

YAO Jun-na, QING Fen

(LaboratoryofGeospatialTechnologyfortheMiddleandLowerReachesofYellowRiver,He′nanUniversity,Kaifeng,He′nan475000,China)

In this paper, Landsat TM image data, digital elevation model DEM, soil type data, base maps and statistical information were used as data source. By using GIS and RS techniques, this paper comprehensively evaluated the ecological environment quality of the sandstone area in 2000 and 2010 by using principal component analysis evaluation model. The results showed that eco-environmental quality of the sandstone area in 2000 and 2010 was ordinary, but it had been improved in 2010 compared with 2000. Assessing eco-environmental quality by using GIS and RS is an effective means for the area with large range, complicated condition, or difficulty to gain the necessary data.

GIS; RS; eco-environmental quality; comprehensive assessment; sandstone

2013-12-02

：2013-12-25

十二五國家科技支撐計劃項目(黃河中游砒砂巖區抗蝕促生技術集成與示范2013BAC05B01)

姚俊娜(1987—),女,河南洛陽人,碩士,主要研究方向為地理性息系統應用。E-mail:yaojunna_07@163.com

秦奮(1966—),男,河南開封人,教授,博士,GIS軟件研究與應用、虛擬地理環境、數據集成與共享。E-mail:qinfun@126.com

X171.1;TP79

：A

：1005-3409(2014)06-0193-05