喀斯特山區生態環境脆弱性綜合評價
——以貴州省黔中地區為例

安宏鋒，徐浩，安寧，王之明，左思藝

(1.貴州省環境監測中心站，貴州貴陽　550001； 2.貴州師范大學山地資源與環境遙感重點實驗室，貴州貴陽　550001)

?

喀斯特山區生態環境脆弱性綜合評價
——以貴州省黔中地區為例

安宏鋒1，徐浩1，安寧2，王之明1，左思藝1

(1.貴州省環境監測中心站，貴州貴陽550001； 2.貴州師范大學山地資源與環境遙感重點實驗室，貴州貴陽550001)

摘要：從自然狀態、人文壓力方面探討建立喀斯特山區生態環境脆弱性指標評價體系，通過變異系數法確定了指標權重，收集了2010年貴州省黔中地區的自然、人文、遙感數據，并以此為基礎對該地區的生態環境脆弱性進行了綜合評價。結果表明：生態環境重度及以上脆弱類型主要分布在裸地、耕地，中度及以下脆弱類型主要分布在草地、灌木林地、有林地。黔中地區生態極強度脆弱區主要分布在北部區域；而中度和強度脆弱區則沿極強度脆弱區分布；輕度脆弱區集中在中部、南部；生態環境脆弱區主要發生在坡度25°以下的區域。

關鍵詞：喀斯特山區；生態環境；脆弱性

自美國學者Clement提出生態過渡帶概念以來，生態脆弱性開始受到廣泛關注和研究，理論認知得到了豐富和發展。王小丹等[1-2]所總結的生態脆弱性的概念在不同學科領域和地域范圍內有不同理解，仍未達成共識。本文趨向于脆弱性即脆弱程度，作為一種過渡性的表示方法，其運用于生態環境評價則表示某一地區生態系統或環境受到干擾時，易從一種狀態轉變為另一種狀態能力的強弱。喀斯特作為生態脆弱區的典型性代表，峰叢-洼地、峰林-谷地以及水系二重系統發育，致使石漠化、水土流失嚴重且人地關系矛盾突出[3-5]。本文從自然狀態、人文壓力兩個主要方面對生態脆弱進行綜合性研究，探討貴州省黔中地區的生態脆弱性強度分布情況及其影響要素，并提出對生態環境保護和治理的相關對策和建議。

1研究區概況

黔中地區位于貴州省中部，東經105°43′25.641″～108°11′47.184″，北緯25°39′3.534″～27°30′58.187″，包括貴陽市六區一市三縣、安順市及黔南州部分區縣，共計20個縣(市、區)、250個鄉鎮。區域國土面積24 317 km2，約占全省的14%，大部分屬于云貴高原的喀斯特丘陵地貌發育區，海拔最低447 m，最高2 177 m，屬亞熱帶季風區。區域內主要水系有烏江水系、北盤江水系和紅水河水系，水資源十分豐富。土地利用類型以林地、耕地、草地為主，林地約占50%，耕地約占30%，草地約占15%。2010年黔中地區人口共計764.78萬人，約占貴州省總人口的22%。城鎮基本分布在山間平地。特殊的喀斯特地形、地貌發育使該地區可利用土地資源較少，人地關系面臨嚴峻挑戰，生態趨于惡化，生態環境脆弱性加劇。隨著經濟社會發展，部分地區人口聚集、城市化加快，治理和改善區域脆弱的生態環境已經被提上社會經濟可持續發展的議事日程。

2研究方法

2.1評價指標因子選擇

生態脆弱性是自然環境與人類活動共同作用的產物。自然要素分為地形、氣候、水文、土壤和植被；人類活動主要體現為次生環境效應。為擺脫以往以大區域為評價單元的概述性分析方法，實現以像元為單元對生態環境脆弱程度的更精確評價，依據周永娟[6]等的理念，從自然狀態和人文壓力兩個方面構建綜合的評價指標體系，最終以遙感圖像像元為基本評價單元對黔中地區生態環境脆弱性進行評價，因此本文選取自然狀態因子和人文壓力與生態環境脆弱性緊密聯系的9個因子作為評價指標，如表1所示。

表1　生態脆弱評價體系Table 1　The evaluation system of ecological vulnerability

2.2數據來源及處理

2.2.1數據來源

以黔中地區2010年30 m分辨率的TM遙感影像數據和DEM數據、2010年第六次人口普查數據、2010年中國區域經濟統計年鑒、《貴州土種志》和貴州省土壤圖、地貌圖和水文地質圖以及2010年土地利用數據等為基礎數據源。

2.2.2數據處理

(1)植被覆蓋度指數

歸一化植被指數和植被覆蓋度指數是反映綠色植被覆蓋情況、生態環境質量的重要指標之一，由于計算簡單、操作方便，在各種環境監測項目中被廣泛運用于反映自然生態環境的好壞。歸一化植被指數的計算公式如下：

NDVI=(NIR-R)/(NIR+R)

(1)

式中，NIR為影像的近紅外波段；R為影像的紅光波段；NDVI為歸一化植被指數，其取值范圍為-1～1。根據歸一化植被指數計算方法，通過進一步計算可得植被覆蓋度指數，其公式如下：

Fv=(NDVI-NDVImin)/(NDVImax-NDVImin)

(2)

式中，Fv為植被覆蓋度指數，其值域范圍為0～1；NDVI為式(1)所得；NDVImax和NDVImin分別為歸一化植被指數的最大值和最小值。

(2)基巖裸露率

石漠化是巖溶地區土地退化的極端形式。作為表征石漠化強度的主要因子，基巖裸露率因子采用Huang[7]等模仿歸一化植被指數提出的歸一化巖石指數和張曉倫[8]等基于NDRI像元二分模型提取基巖裸露率的TM波段運算方法，自動提取基巖裸露信息，其計算公式為：

NDRI=(SWIR-R)/(SWIR+R)

(3)

Fr=(NDRI-NDRImin)/(NDRImax-NDRImin)

(4)

式中，SWIR表示短波紅外波段(TM影像第五波段)；R是影像的紅光波段；NDRI為歸一化巖石指數；Fr為基巖裸露率；NDRImax是歸一化巖石指數的最大值；NDRImin為歸一化巖石指數最小值。

(3)土壤侵蝕

采用通用土壤流失方程對黔中地區的土壤侵蝕量進行評價，其公式為：

USLE=R·K·LS·C·P

(5)

式中，R表示降雨侵蝕力因子，采用地球系統科學數據共享網數據；C為植被覆蓋與管理因子，其值利用蔡崇發[9]等建立的C值與植被覆蓋度指數的回歸方程并結合實際情況進行判斷，如水體與建設用地以及水田沒有土壤侵蝕，則C值為0，裸地的土壤侵蝕最強C值為1，其他土地類型的C值為方程算出的平均值代替；P因子、K因子參考楊廣斌[10-11]結合貴州實際情況對水土保持措施因子P、土壤可蝕性因子K的賦值，其對應的類型取值情況如表2、表3所示。

表2　黔中地區各地類P、C取值Table 2　The P and C value table of land types in central region of Guizhou province

表3　黔中地區各土壤類型K取值Table 3　The K value table of soil types in central region of Guizhou province

對于LS坡長-坡度因子，可由30 m空間分辨率的DEM數據提取坡度數據，然后根據公式求得。坡度因子計算公式為：

(6)

坡長因子則引用Wischmeier和Smith提出的坡長因子計算方法，其公式為：

L=(λ/22.13)m

(7)

(8)

(9)

式中，θ為坡度，°；λ為坡長，m。

將上述各因子數據轉化為30 m×30 m的柵格數據，代入公式(5)計算得出黔中地區的土壤侵蝕圖，經標準化處理后如圖1所示。

圖1　黔中地區標準化土壤侵蝕圖Fig.1　The map of standardized soil erosion in central region of Guizhou province

2.2.3變異系數法權重確定

數學分析方法和權重分析方法在地理分析中得

到了廣泛的應用研究。本文采用變異系數法來客觀確定各評價因子指標的權重值。變異系數法是直接利用各項指標所包含的信息，通過計算得到指標的權重，是一種客觀賦權的方法。

為了更好地實現生態脆弱性的評價，將各項評價指標進行離差標準化處理，然后將處理后的數據運用變異系數法公式求取各指標權重值。其標準化處理公式為：

(10)

根據各指標標準化處理后的柵格數據統計其均值和標準差，并通過變異系數法公式算各指標的變異系數和權重值，其公式如下：

(11)

(12)

將表4所得到的各指標的權重值代入公式(13)得到黔中地區生態脆弱性情況，并在地理信息系統軟件中用自然段點法把數據分為輕度脆弱、中度脆弱、強度脆弱和極強度脆弱區，如圖2所示。其中，水域和水田兩種土地利用類型作為掩膜區，默認為非生態脆弱。各脆弱性強度類型面積如表5所示。

(13)

式中，Wi為i項指標的權重值；Xi為i項指標的標準化數值。

表4　黔中地區變異系數方法求取權重系數表Table 4　The index weight was determined by the coefficient of variation method in central region of Guizhou province

圖2　黔中地區生態脆弱強度Fig.2　The ecological vulnerability intensity in central region of Guizhou province表5　黔中地區脆弱性強度類型面積統計Table 5　The area of vulnerability intensity in central region of Guizhou province

脆弱性強度類型面積/km2面積百分比/%非生態脆弱2244.289.23輕度脆弱5392.5322.18中度脆弱7738.2531.82強度脆弱6234.7925.64極強度脆弱2707.1611.13總計24317100

3結果分析

3.1黔中地區生態環境脆弱性現狀

黔中地區生態環境脆弱性按分布面積從大到小依次是中度脆弱、強度脆弱、輕度脆弱、極強度脆弱、非生態脆弱類型，以輕度脆弱、中度脆弱、強度脆弱為主，三者之和所占比例達到80%左右，中度和強度脆弱面積占57%左右。結果顯示，黔中地區生態脆弱性現狀不容樂觀。

3.2黔中地區生態環境脆弱區空間分布

從圖2來看，生態極強度脆弱區主要分布在黔中北部開陽縣、息烽縣大部以及平壩縣靠花溪區部分、清鎮市以西、修文縣北部和福泉市、白云區的小部分區域；而中度和強度脆弱區域則沿極強度脆弱區域分布，中度脆弱區域所占面積比例較小；輕度脆弱區集中分布在同屬于“花生殼”狀的喀斯特地貌黔中地區中部、南部，這些區域開墾率低，對自然狀態造成的人為擾動較小，脆弱性相對較輕。

通過疊加地貌數據，統計分析得出黔中地區脆弱性等級地貌類型面積比例，如表6所示。在山地、丘陵地貌區則相對較多，而在臺地和盆地地區，生態脆弱區的種類分布相對很少；輕度和中度脆弱在山地上分布較多，而強度和極強度脆弱則在丘陵和山地上發育相對較多，是因為在純度較低的可溶性碳酸鹽巖上發育的山地、丘陵相對高度較小，坡度不是很大的區域有農耕基礎，人為活動較多(坡耕地、建筑等)，使得強度及以上脆弱區域面積比例較大。

表6　黔中地區脆弱性等級地貌類型面積比例Table 6　The area ratio of land form types by the vulnerability level in central region of Guizhou province

3.3黔中地區生態環境脆弱性的主導因素

由表1和表4所示黔中地區生態脆弱性各評價指標的權重值可知，影響生態脆弱性的因子主要為土壤侵蝕因子、單位糧食產量和人均耕地面積，三者權重貢獻率為0.56；其次是降水、土層厚度以及開墾率，權重貢獻率為0.27。自然環境狀態因子占據了0.63的貢獻率，人文壓力因子則占據了0.37。自然狀態因子是生態脆弱的基礎，而不合理的人類活動則會加劇生態脆弱的過程，因而生態脆弱是人與自然環境共同作用的結果。

3.4生態脆弱性與土地利用類型的關系

統計不同土地利用類型的生態脆弱性數據，如表7所示，為黔中地區各土地利用類型中的脆弱性面積比例。由于本文不探討水體生態脆弱性，而水田長期存在人為管理和耕作，不存在生態脆弱，因而對于水域和水田兩種類型不加以統計。

表7　黔中地區各土地利用類型中的脆弱性面積比例Table 7　The vulnerability area ratio of land types in central region of Guizhou province

由表7可知，裸地無植被覆蓋，土壤和礫石裸露在外面，水土流失最為嚴重，因而其發生中度及以上生態脆弱性的面積比例最大；旱地作為農業生產的重要場所，由于人們不合理的土地利用方式容易造成水土流失、生態脆弱加劇；其他林地(包括園地、果林、苗圃等)作為耕地的一種轉型利用方式，既實現了對生態環境的保護，又實現了坡耕地土地利用的轉型升級和經濟生產增值，是坡耕地優化利用的一種良好方法，但其生態結構單一，強度生態脆弱發生的可能性也不容忽視，應加強防護和管理；城鄉用地以及工礦、交通用地經過很好的工程處理，可改變性很小，因而強度及以上脆弱發生率相對較小；草地、灌木林地、有林地因人類影響和破壞相對較小，且屬于一種生態蓄積結果，如果遭到嚴重破壞，特別是在喀斯特上發育的綠色植被受破壞后將很難恢復，強度生態脆弱發生率最高。

3.5生態脆弱性與坡度的關系

人類活動致使生態脆弱最突出的表現就是農業耕種。貴州是一個處于山區的農業大省，擁有大面積的坡耕地，因而坡度是造成該地區生態系統脆弱的自然相關因子。采用水土保持中的坡度劃分原則(微坡為<5°，較緩坡為5°～8°，緩坡為8°～15°，較陡坡為15°～25°，陡坡為25°～35°，極陡坡為>35°)劃分坡度，并統計得到生態脆弱性在坡度類型的面積情況，如表8所示。

由表8可知，生態脆弱性主要發生在較陡坡及其以上區域，即25°以上的區域，其中以緩坡和較陡坡為主，達到各脆弱性等級面積的0.55以上，而在15°以下區域也分布有強度和極強度脆弱性，且面積比達到的0.33，說明除自然條件外，人類大量不合理的活動(耕作、開采、建筑等)造成了嚴重的生態脆弱性，因而應強化對人類活動區域的管理。

表8　黔中地區各生態脆弱性中坡度類型面積比例Table 8　The area ratio of slope types by the vulnerability level in central region of Guizhou province

4結論及建議

(1)黔中地區生態環境脆弱性以輕度、中度、強度脆弱為主，中度和強度脆弱面積占57%左右，地區生態環境脆弱，必須對黔中地區的生態環境保護給予足夠重視。

(2)在土層厚度的數據選擇上，本文使用了土壤的平均厚度數據，若能收集到更加精確的土層厚度數據，評價結果會更加準確。由于黔中地區地形起伏大，采用30 m空間分辨率的DEM數據計算坡度和坡長存在較大誤差，利用此值計算的土壤侵蝕的誤差也較大，若能收集到更加精確的DEM數據也會使評價結果更加準確。缺乏鄉鎮單位糧食產量、人均GDP數據，這些數據是以縣域為單位進行統計，對整個評價結果的精度會帶來一定誤差。

(3)自然狀態因子是生態脆弱的基礎，而不合理的人類活動則會加劇生態脆弱的過程，因而生態脆弱是人與自然環境共同作用的結果。人類不合理的開發利用，如在喀斯特山丘、坡地上不合理開墾耕地，致使喀斯特山區原本稀薄的土層，或無土壤覆蓋，或巖石裸露，造成生態更加脆弱。

(4)裸地治理、坡耕地的水土保持生態建設和土地利用類型優化可以很好地實現生態脆弱性治理，與經濟發展相協調。山地、丘陵地區可以在適宜利用的坡耕地上引導種植果林、經濟林等，或加強工程治理，如在坡耕地區域采取林耕相間的耕作模式，加強對裸地、坡耕地的水土保持生態建設和合理利用，保護土壤，才能實現生態和生產的良性循環發展，改善生態環境脆弱性。

(5)生態脆弱性在區域空間分布上與坡度、人為活動、可溶性碳酸性鹽巖的分布等有密切聯系。其中，坡度較小、人類活動較多的山地和丘陵地區生態脆弱度較大；而在純碳酸巖鹽的喀斯特地貌上，極強度生態脆弱發育較多，中度和強度脆弱則廣泛分布在亞喀斯特地貌上。因此，加強對不同純度碳酸巖鹽發育的喀斯特地貌環境的保護和治理研究，對改善地區生態脆弱性有著重要意義。

參考文獻(References)：

[1]王小丹, 鐘祥浩. 生態環境脆弱性概念的若干問題探討[J]. 山地學報, 2003, 21(增刊1): 21- 25.

[2]李鶴, 張平宇, 程葉青. 脆弱性的概念及其評價方法[J]. 地理科學進展, 2008, 27 (2): 18- 25.

[3]謝丹, 劉小麗, 劉毅, 等. 云貴可持續發展定位挑戰與對策[J]. 環境影響評價, 2014(2):29- 31.

[4]劉小麗, 劉毅, 任景明, 等. 云貴地區生態環境現狀及演變態勢風險分析[J]. 環境影響評價, 2015, 37(1): 27- 30.

[5]李天威, 任景明, 曾思育, 等. 云貴地區水環境現狀特點及未來難點解析[J]. 環境影響評價, 2015, 37(3): 40- 44.

[6]周永娟, 王效科, 歐陽志云. 生態系統脆弱性研究[J]. 生態經濟, 2009(11): 165- 189.

[7]Huang Qiuhao, Cao Yunlong. Mapping Karst Rock in Southwest China[J]. Mountain Research and Development, 2009, 29(1): 14- 20.

[8]張曉倫, 甘淑. 基于NDRI像元二分模型的石漠化信息提取研究[J]. 新技術新工藝, 2014(1): 72- 75.

[9]蔡崇發, 丁樹文, 史志華, 等. 應用USLE模型與地理信息系統IDRISI預測小流域土壤侵蝕的研究[J]. 水土保持學報, 2000, 14(2): 19- 24.

[10]楊廣斌. 基于網格數據的喀斯特山區土壤侵蝕研究[D]. 貴陽: 貴州師范大學, 2004: 36- 44.

[11]楊廣斌, 王濟, 蔡雄飛, 等. 喀斯特山區土壤侵蝕評價及數據模擬[M]. 北京: 氣象出版社, 2013: 75- 80.

收稿日期：2015-07-02

基金項目：2015年貴州省地下水基礎環境狀況調查評估(2015A080)

作者簡介：安宏鋒(1982—)，男，河南洛陽人，工程師，碩士，主要從事生態環境監測，E-mail：hongfeng_1225@163.com 通訊作者：徐浩(1975—)，男，貴州天柱人，高級工程師，本科，研究方向為環境監測與評價，E-mail：1026022117@qq.com

DOI:10.14068/j.ceia.2016.04.014

中圖分類號：X826

文獻標識碼：A

文章編號：2095-6444(2016)04-0051-06

Comprehensive Evaluation of Eco-environmental Vulnerability in Karst Mountain Areas—Central Region of Guizhou Province as An Example

AN Hong-feng1, XU Hao1, AN Ning2, WANG Zhi-ming1, ZUO Si-yi1

(1.Guizhou Environmental Monitoring Center Station, Guiyang 550001, China; 2. Key Laboratory of Mountain Resources and Environmental Remote Sensing, Guizhou Normal University, Guiyang 550001, China)

Abstract:From the aspects of natural environment and effect of human activities, the evaluation system of eco-environmental vulnerability index in karst mountain areas was discussed in this paper. The index weight was determined by the coefficient of variation method. Based on the natural environment, humanity and remote sensing (RS) data, the eco-environmental vulnerability in central region of Guizhou Province was evaluated. Results showed that the areas with severe and above eco-environmental vulnerability were mainly distributed in bare land and cropland, while those of moderate and below were distributed in grassland, shrub land and forest. In central region of Guizhou Province, the extremely severe vulnerable area was mainly distributed in the northern region; the moderate and severe vulnerable areas were distributed along the extremely severe vulnerable area; the mildly vulnerable area were distributed in the central and southern regions; the eco-environmental vulnerable areas were mainly distributed in regions with the slope less than 25 degrees.

Key words:karst mountain area; eco-environment; vulnerability