榆林地區水土流失敏感性時空特征

李團勝，史小惠，康歡歡，石玉瓊，閆　穎

(1.長安大學 地球科學與資源學院，陜西 西安　710054;2.山西農業大學 林學院， 山西 太谷　030801)

?

·環境科學·

榆林地區水土流失敏感性時空特征

李團勝1，史小惠1，康歡歡1，石玉瓊2，閆穎1

(1.長安大學 地球科學與資源學院，陜西 西安710054;2.山西農業大學 林學院， 山西 太谷030801)

用10km×10km的網格,把榆林地區劃分為527個單元樣地,采用加權和法,計算每個單元的水土流失敏感性指數,并在Arcgis中進行Kriging空間插值,分析水土流失敏感性時空特征。結果表明：全區無不敏感級別。2000年中度敏感區面積占全榆林面積的31.33%。高度敏感區面積占51.808%。極度敏感區占16.862%。2005年中度敏感區面積占35.876%。高度敏感區面積占48.153%。極度敏感區占15.971%。2010年輕度敏感區占0.229%。中度敏感區面積占28.047%。高度敏感區面積占44.176%。極度敏感區面積占27.548%?？傮w來說,東部水土流失敏感性較高，西部敏感性較低。從2000年到2010年,中度敏感性和高度敏感性面積減少了,而極度敏感性面積增加了。

水土流失;敏感性;土地覆蓋;地形起伏;土壤質地;榆林

生態敏感性是在現有自然環境背景下, 人類活動干擾和自然環境變化導致發生區域生態或環境問題的難易程度及其可能性大小[1],早在1980年,Suffling就提出了從生態敏感到生態擾動的指標[2],之后,Horne研究了澳大利亞雨林對選擇性伐木的生態敏感性[3],Rodriguez研究了大陸架生態敏感性[4],近年來,Renetzeder研究了奧地利農業景觀生境對氣候變化的敏感性[5],Thomas以熱帶地區為例研究了景觀對快速的環境變化的敏感性[6]。國內的研究也逐步增多,取得了一系列成果。目前對生態敏感性的研究主要集中在對單一環境要素的研究,如酸沉降、鹽漬化、沙漠化、土壤侵蝕等,并呈現出從單一生態敏感性因子研究向綜合敏感性研究發展的趨勢[7]。多數單要素研究是研究水土流失或土壤侵蝕敏感性的,而且研究集中在黃土高原[8-10]、東北地區[11-13]以及西南石質山區[14-18],另外還有針對整個國家[19]、 西藏[20]、鄂爾多斯[21]、桐柏山——大別山區[22]、海南[23]的研究。水土流失是受多種因素影響的,有的因素和人類活動密切相關,如土地利用結構、水土保持措施等,還有一些因素,如氣候、地形等屬于自然因素,雖然也受人為活動影響,但主要受制于自然因素。從目前的研究來看,大多數研究是按照國家環境保護部的《生態功能區劃技術暫行規程》中的生態系統敏感性評價方法進行的,也有按照水土流失通用方程進行評價的,研究方法很不統一。榆林地區位于黃土高原和毛烏素沙漠的交接地帶,生態環境脆弱,水土流失嚴重,大體以長城為界,西北部是屬于毛烏素沙漠的風沙灘地區,東南部是屬于黃土高原的黃土丘陵溝壑區,本文基于GIS,利用加權和法,分析了榆林地區水土流失敏感性的時空特征。

1　研究區域概況

榆林地區位于陜西省最北部,東經107°15′到111°15′, 北緯36°57′到39°35′之間, 總面積43 578 km2。地貌上可分為兩大單元,即位于毛烏素沙漠南緣的風沙灘地區和黃土丘陵溝壑區,二者大體以長城為界,毛烏素沙漠南緣的風沙灘地區屬鄂爾多斯高原的一部分,黃土丘陵溝壑區屬黃土高原的一部分。榆林氣候屬溫帶半干旱大陸性季風氣候,四季分明,年平均氣溫10℃,年平均降水量400 mm左右,年內降水量多集中于6-9月份。榆林地區是陜西省唯一包括外流河和內陸河流域的地市級區域。外流河全部屬于黃河水系,內陸河分布在神木、定邊縣北部的沙漠閉流區。榆林主要土類有風沙土、黃綿土、新積土、灰鈣土、栗鈣土、黑壚土、鹽土等。地帶性植被為干草原和森林草原。主要植被有森林、干草原、灌叢、草甸、沙生植被、鹽生植被等。榆林礦產資源極其豐富,有煤、天然氣、石油、鹽等,目前為陜北能源化工基地的核心區域,榆林轄1區11縣。

2　研究方法

在2000年、2005年和2010年遙感影像解譯得到的3期土地覆蓋圖的基礎上,把研究區域用10km×10km的網格進行劃分,共分為527塊單元樣地,計算每塊樣地的水土流失敏感性指數,然后在ArcGIS中用Kriging插值法進行空間插值,得到榆林地區水土流失敏感性指數圖,在進行水土流失敏感性評價時,考慮土壤質地、地形起伏以及土地覆蓋3個因子,由于缺乏降水量資料,同時土地覆蓋本身間接地反映了降水量的大小,因而沒有考慮降水因素。首先把土壤質地、地形起伏以及土地覆蓋3個因子進行敏感性分級,各分為不敏感(1)、輕度敏感(2)、中度敏感(3)、高度敏感(4)和極度敏感(5)共5級(表1),然后計算每一樣地內每一因子的敏感性指數,其計算公式如下：

(1)

式中,sij為第j塊樣地水土流失某一因子i的敏感性指數,sij值大表明敏感性高;cijk為j塊樣地內第k級所占面積的比例,wik為因子i第k級的權重,不敏感的權重為0.05,輕度敏感的權重為0.1,中度敏感的權重為0.2,高度敏感的權重為0.25,極敏感的權重為0.4。按照上述公式分別計算出土壤質地、植被、地形起伏度等3個單因子敏感性指數。其中地形起伏敏感性指數的計算是根據榆林地區DEM數據,以5 km×5 km移動窗口為標準,獲取得到地形起伏度,并參考文獻[12]中的劃分指標劃分為5個等級,按照等級權重指數計算得出527塊樣地的地形起伏度敏感性指數。

表1　榆林地區土壤質地、地形起伏和土地覆蓋敏感性分級標準

再計算每一樣地內水土流失敏感性指數,其公式如下:

(2)

式中:xj為第j塊樣地水土流失敏感性指數,xj值大表明敏感性高;sij為j塊樣地i因子敏感性指數,wi為i因子的權重,其中地形起伏度因子的權重為0.4,土壤質地因子的權重為0.2,土地覆蓋因子的權重為0.4。

在作水土流失敏感性指數插值圖時先采用自然斷點法,把敏感性指數分為10級作圖,以顯示敏感性指數的空間分布,然后根據計算結果,結合該研究區實際,把敏感性分為5級，以反映敏感性相對大小。敏感性指數小于0.1為不敏感,[0.1～0.2)為輕度敏感,[0.2～0.24)為中度敏感,[0.24～0.28)為高度敏感,敏感性指數等于大于0.28為極度敏感。

3　結果與分析

2000年,榆林地區水土流失敏感性指數在0.186 3～0.323之間(圖1a),東南部及東北部地區敏感性指數較高,西部及中部地區敏感性指數較低。子洲縣、清澗縣、綏德縣、米脂縣和吳堡縣,其敏感性指數在0.220 7～0.312 8之間;神木縣的高家堡鎮、欄桿堡鎮、大柳塔鎮等,其敏感性指數在0.238 3～0.319 6之間;府谷縣的大昌汗、老高川等,其敏感性指數在0.262～0.295之間。其余地區敏感性指數較低。

2005年,榆林地區水土流失敏感性指數在0.171 5～0.319 4之間(圖1b),東南部、中部榆陽區及東北部敏感性指數較高;西部及中東部地區指數較低。子洲縣、清澗縣、綏德縣、米脂縣和吳堡縣其敏感性指數在0.215～0.319 4之間;神木縣的高家堡鎮、欄桿堡鎮、錦界鎮、爾林兔鎮、孫家岔鄉、大柳塔鎮等,其敏感性指數在0.239 5～0.309 7之間;府谷縣的大昌汗、麻鎮、黃甫、墻頭等,其敏感性指數在0.255 2～0.305 8之間;榆陽區的榆陽鄉、青云鄉、劉千河鄉、牛家梁鎮等,其敏感性指數在0.237 2～0.298 8之間,這些地區水土流失敏感性指數較高,其余地區敏感性指數較低。

2010年,榆林地區水土流失敏感性指數在0.147 2～0.335 2之間(圖1c),東南部及東北部地區敏感性指數較高,西部及中部地區敏感性指數較低。子洲縣、清澗縣、綏德縣、米脂縣和吳堡縣敏感性指數在0.228 3～0.335 2之間;神木縣的欄桿堡鎮、王家墩、武家莊、大柳塔鎮、孫家岔鄉等,其敏感性指數在0.266 9～0.305 4之間;府谷縣的大昌汗、老高川、廟溝門等,其敏感性指數在0.267 1～0.303 7之間;榆陽區的榆陽鄉、青云鄉等,其敏感性指數在0.250 7～0.302 7之間。

從敏感性分級圖(圖2)可看出,榆林地區主要是中度敏感、高度敏感和極度敏感3個級別,無不敏感級別。2000年,中度敏感區的面積占全區域面積的31.33%,主要分布在西部的定邊和靖邊,以及東部的佳縣大部分地區。高度敏感區的面積占51.808%,主要分布在東北部的神木和府谷、中部的榆陽和橫山、東南部的米脂、吳堡。極度敏感區面積占16.862%,主要分布在東南部的清澗、綏德、子洲,東北部神木的大柳塔、欄桿堡,府谷的大昌汗、麻鎮、墻頭等地。2005年中度敏感區的面積占全區域面積的35.876%, 主要分布在西部的定邊和靖邊及東部的佳縣,中部的橫山。高度敏感區的面積占48.153%,分布同2000年,只是中部的橫山縣2000年主要是高度敏感區,而2005年有相當一部分轉化為中度敏感區。極度敏感區占15.971%,主要分布在東南部的清澗、綏德、子洲,東北部神木的大柳塔、爾林兔、府谷的大昌汗、麻鎮、墻頭等地。2010年,輕度敏感區占全區域面積的0.229%,主要分布在定邊縣的西北部。中度敏感區的面積占全區域面積的28.047%,主要分布在定邊和靖邊。高度敏感區面積占44.176%,主要分布在佳縣、橫山、榆陽以及神木。極度敏感區面積占27.548%,主要分布在東南部的清澗、綏德、子洲、米脂、吳堡,東北部神木的大柳塔、孫家岔、永興、欄桿堡、瓦羅等地,府谷的大部分地區。

圖1　榆林地區2000年(a)、2005年(b)和2010年(c)水土流失敏感性指數的空間分布Fig.1　Spatial distribution of soil and water loss sensitivity index of Yulin in 2000(a), 2005(b) and 2010(c)

可見,從2000年到2005年,中度敏感區面積略有上升,而高度敏感區面積和極度敏感區面積略有下降。從2005年到2010年,中度敏感區面積和高度敏感區面積有所下降,而極度敏感區增加較多。從2000年到2010年總體情況來看,中度敏感區和高度敏感區面積有所減少,而極度敏感區面積有所增加,尤其是府谷和神木的極度敏感區面積增加較多,東南部的極度敏感區面積也有所增加。且極度敏感區主要在黃土丘陵溝壑區,總體來看,黃土丘陵溝壑區的敏感性增強了,變化最明顯的就是佳縣、橫山縣和府谷縣,佳縣大部分區域在2000年和2005年時是中度敏感范圍,而到2010年,全縣范圍處于高度敏感區范圍。橫山縣2000年中度敏感區范圍比較小,2005年中度敏感區范圍擴大,高度敏感區范圍縮小,而到2010年,全縣絕大部分處于高度敏感區范圍,中度敏感區范圍縮小。府谷縣在2000及2005年大部分地區處于高度敏感區,而到了2005年,大部分地區處于極度敏感區。

圖2　榆林地區2000年(a)、2005年(b)和2010年(c)水土流失敏感性分布Fig.2　Spatial distribution of soil and water loss sensitivity of Yulin in 2000(a), 2005(b) and 2010(c)

榆林地區氣候干旱,降水稀少,生態環境脆弱,除了自然原因外,就是人為因素,作為人為因素,對榆林地區影響最大的主要是兩個方面,一是退耕還林還草政策的實施,另一個是陜北能源化工基地的建設。原國家計委于2003年3月正式批準陜北能源化工基地在榆林啟動建設以來,榆林就大規模招商引資,目前已形成榆林經濟開發區、神府經濟開發區,榆神煤化工業園區、榆橫煤化工業園區、府谷煤電化工業園區、綏米佳鹽化工業園區、靖邊煤油氣綜合利用產業園及吳堡煤焦化工業園區的“兩區六園”格局。2010年以前的10年內,累計建成各類重點項目70個,完成投資1 260億元,投資在100億元以上的在建、前期大型產業項目達12個,僅2009年一年就組織實施了90個重點建設項目,50個重大前期項目,有6個項目建成投產。根據榆林工業發展空間布局示意圖,榆林工業園區多達26個,其中2個位于定邊和靖邊縣境內,其余24個都集中在其他縣區?？梢娮?003年以來,榆林的能源化工基地建設如火如荼。退耕還林還草,有利于生態環境的改善,而大規模建設與礦產資源的開發,特別是近10年來,隨著神木錦界、店塔、大柳塔、府谷新民、廟溝門、榆陽魚河、鎮川、牛家梁和曹家灘等26個大型工業園區、開發區建設突飛猛進,290個能源重化工企業的重大項目及基礎設施建設紛紛上馬,幾乎遍布整個區域,由于盲目資源開發及園區建設,使本來就非常脆弱的生態環境雪上加霜,據研究,該地區僅煤田開發一項就造成1.73×104hm2的植被被毀,新增沙漠化土地近2.0×104hm2。另據典型調查,榆陽區小壕兔鄉因礦區大量抽取地下水,居民原有2～3m深的飲用水井如今挖至10m多深才能見到水,一些梁峁區的防護林因地下水位下降出現枯死的現象嚴重。盡管陜北能源化工基地建設是本著把資源開發好,把生態保護好的原則,但礦產資源開發尤其是地下礦產資源的開發不可避免的會破壞地下水體系,也因而必然殃及到地表水系統,同時大量的開發建設者蜂擁而至,無疑增加了該區域的負擔,使原本承載力比較小的自然環境更是不堪重負,勢必增加水土流失的風險,增強水土流失的敏感性。因此,如何協調脆弱的生態與開發建設之間的矛盾,如何走可持續發展的道路,是各級部門應該深思的問題。

4　結　論

1)2000年,榆林地區水土流失敏感性指數在0.186 3～0.323之間,2005年,水土流失敏感性指數在0.171 5～0.319 4之間,2010年,水土流失敏感性指數在0.147 2～0.335 2之間。

2)從2000年到2010年,榆林地區無不敏感區,主要屬于中度、高度和極度敏感區,中度敏感性和高度敏感性面積減少了,極度敏感區面積增加了?？傮w說來,東北部和東部敏感性高于西部。礦產資源的開發以及能源化工基地的建設會進一步增強水土流失的敏感性。

3)退耕還林還草政策的實施增大了林草覆蓋率,有利于水土保持,使水土流失敏感性降低,而礦產資源的開發和能源化工基地建設等人為活動的增強又增大了水土流失的敏感性。

[1]歐陽志云,王效科,苗鴻. 中國生態環境敏感性及其區域差異規律研究[J]. 生態學報,2000,20(1):9-12.

[2]SUFFLING R. An index of ecological sensitivity to disturbance, based on ecosystem age, and related to landscape diversity[J]. Journal of Environmental Management, 1980, 10(3): 253-263.

[3]HORNE R, HICKEY J. Ecological sensitivity of Australian rain-forests to selective logging[J]. Australian Journal of Ecology, 1991, 16(1): 119-129.

[4]RODRIGUEZE V L. Ecological sensitivity atlas of the argentine continental shelf [J]. International Hydrographic Review, 1992, 69(2): 47-53.

[5]RENETZEDER C, KNOFLACHER M, LOIBL W, et al. Are habitats of Austrian agricultural landscapes sensitive to climate change?[J]. Landscape and Urban Planning, 2010, 98:150-159.

[6]MICHAEL F. Thomas. Landscape sensitivity to rapid environmental change-a Quaternary perspective with examples from tropical areas[J]. Catena, 2004,55: 107-124.

[7]顏磊,許學工,謝正磊,等.北京市域生態敏感性綜合評價[J].生態學報,2009,29(6): 3117-3125.

[8]李國強,陳利頂,高啟晨,等.黃土高原地區西氣東輸工程沿線水土流失敏感性評價[J]. 水土保持學報,2003,17 (6):55-58.

[9]LIU L, LIU X H. Sensitivity analysis of soil erosion in the Northern Loess Plateau[J]. Procedia Environmental Sciences，2010,2:134-148.

[10] 劉康,康艷,曹明明,等. 基于GIS的陜西省水土流失敏感性評價[J]. 水土保持學報,2004,18 (5):168-170.

[11] 陳建軍,張樹文,李洪星，等. 吉林省土壤侵蝕敏感性評價[J]. 水土保持通報,2005,25 (3):49-53.

[12] 萬忠成,王治江,董麗新，等.遼寧省生態系統敏感性評價[J]. 生態學雜志,2006,25 (6):677-681.

[13] 吳志軍,劉淼,王子楠. 基于GIS的鞍山市水土流失敏感性評價[J]. 遼寧工程技術大學學報(自然科學版),2008,28(Sup):131-133.

[14] 盧遠; 華璀; 周興. 基于GIS的廣西土壤侵蝕敏感性評價[J]. 水土保持研究,2007,14 (1):98-100.

[15] 莫斌,朱波,王玉寬,等. 重慶市土壤侵蝕敏感性評價[J].水土保持通報,2004, 24(5): 45-59.

[16] 李陽兵, 邵景安, 王世杰, 等.基于巖溶生態系統特性的水土流失敏感性評價[J].山地學報,2007,25(6):671- 677.

[17] 羅俊,王克林,陳洪松,等. 桂西北喀斯特地區水土流失敏感性評價[J].長江流域資源與環境,2009,18(6):579-584.

[18] 凡非得,王克林,熊鷹,等.西南喀斯特區域水土流失敏感性評價及其空間分異特征[J]. 生態學報,2011,31(21):6353-6362.

[19] 王效科,歐陽志云,肖寒,等.中國水土流失敏感性分布規律及其區劃研究[J]. 生態學報,2001,21(1):14-19.

[20] 王小丹,鐘祥浩,范建容.西藏水土流失敏感性評價及其空間分布規律[J]. 地理學報,2004,59(2):183-188.

[21] 王再嵐,李政海,馬中,等. 鄂爾多斯東勝地區不同生態功能區的土壤侵蝕敏感性[J].生態學報,2009,29(1):484-491.

[22] ZHANG Rong-hua, LIU Xia, GARY C. Heathman,Xiaoyou Yao,Xuli Hu,Guangcan Zhang. Assessment of soil erosion sensitivity and analysis of sensitivity factors in the Tongbai-Dabie mountainous area of China[J]. Catena，2013，101：92-98.

[23] QIU Peng-hua, XU Song-jun, XIE Gen-zong, et al. Analysis of the ecological vulnerability of the western Hainan Island based on its landscape pattern and ecosystem sensitivity[J]. ACTA ECOLOGICA SINICA, 2007,27(4):1257-1264.

[24] WANG Xiao-ming, BIAN Zheng-fu. The Implications of Ecological Sensitivity on Exploitation of Unutilized Land: A Case Study in Ji′Nan City, China[J]. Procedia Environmental Sciences,2011, 10: 275-281.

(編輯徐象平)

Spatial-temporal characteristics of sensitivity to soil and water loss in Yulin, China

LI Tuan-sheng1, SHI Xiao-hui1, KANG Huan-huan1, SHI Yu-qiong2,YAN Ying1

(1.College of Earth Science and Resources, Chang′an University, Xi′an 710054, China;2.College of Forestry, Shanxi Agricultural University, Taigu 030801, China)

Yulin was divided into 527 units by 10km×10km grids. The index of sensitivity to soil and water loss for each unit was calculated by weighted sum method, and then interpolated by Kriging in Arcgis in order to analyze the spatial-temporal characteristics of sensitivity to soil and water loss in Yulin. The results showed that there is no insensitive level in Yulin. In 2000, 31.33% of the total area was moderate level, 51.808% of the total area being high level, 16.862% of the total area being extreme level. In 2005, 35 876% of the total area was moderate level, and 48.153% of the total area being high level, 15.971% of the total area being extreme level. In 2010, 0.229% of the total area was mild level, 28.047% of the total area being moderate level, 44.176% of the total area being high level, 27.548% of the total area, extreme level. All in all, sensitivity to soil and water loss is high in the east and low in the west part. The area of moderate and high sensitivity to soil and water loss decreased, and that of extreme sensitivity increased from 2000 to 2010.

soil and water loss; sensitivity; land cover; relief; soil texture; Yulin

2015-03-09

國家自然科學基金資助項目(31140042)；陜西省自然科學基礎研究計劃基金資助項目(2009JM5004)；陜西省軟科學計劃基金資助項目(2007RK21); 長安大學中央高?；鹳Y助項目(2013G1271102)

李團勝，男，陜西武功人，博士，教授，從事生態環境與土地利用研究。

P941.73；Q149

A

10.16152/j.cnki.xdxbzr.2016-02-020