王 宇，段永紅，白 杰

（山西農(nóng)業(yè)大學資源環(huán)境學院，山西太谷030801）

五臺山地區(qū)的動植物資源豐富，區(qū)域內(nèi)植物種類達600余種，植被的垂直地帶性明顯，植被種類以草地為主，有華北地區(qū)較為完整的高山草甸草原生態(tài)系統(tǒng)[1]。五臺山自然保護區(qū)的建立對于研究這一地區(qū)的生物多樣性，保障區(qū)域內(nèi)高山草甸生態(tài)系統(tǒng)的穩(wěn)定性，維持這一地區(qū)的生態(tài)平衡都有重要的意義。近年來，人類對保護區(qū)及周邊地區(qū)的開發(fā)力度逐年增強，風景區(qū)內(nèi)的旅游活動及景點的開發(fā)也越來越頻繁，這些活動都對保護區(qū)內(nèi)的植被造成了極大的威脅。對保護區(qū)內(nèi)破壞狀況的調(diào)查對于制定五臺山地區(qū)的可持續(xù)發(fā)展戰(zhàn)略，治理周邊區(qū)域的開發(fā)亂象，發(fā)展生態(tài)旅游事業(yè)都有積極的作用。

遙感技術(shù)因其覆蓋面積大、實時性強、獲取速度快等優(yōu)點被廣泛應用于農(nóng)業(yè)、林業(yè)、地質(zhì)、水文、軍事、環(huán)保等領(lǐng)域。我國應用遙感方法對植被的調(diào)查始于20世紀20年代[2]，遙感大面積同步觀測的特點及其與GPS，GIS的綜合運用使其在這方面的研究越來越多，研究方向也越來越廣泛。20世紀90年代末，遙感技術(shù)被應用于保護區(qū)的調(diào)查與動態(tài)監(jiān)測[3]，這一領(lǐng)域的研究主要分為：對保護區(qū)土地利用動態(tài)變化的研究[4]；對保護區(qū)植被變化的研究[5]；對保護區(qū)人為活動的研究[6]。隨著定量遙感的迅速發(fā)展，使用遙感技術(shù)可以有效地對林業(yè)資源和植被變化實施快速監(jiān)測，并可進一步分析評估其誘因[7]。

研究使用 1990，1995，2000，2005，2010，2015 年6期TM影像分析了近25 a來五臺山地區(qū)的植被變化情況。本研究有利于對五臺山的植被變化狀況進行動態(tài)監(jiān)測，對于研究造成五臺山植被破壞的原因并據(jù)此提出對策有重要的意義。

1 研究區(qū)與數(shù)據(jù)來源

1.1 研究區(qū)概況

五臺山位于山西省東北部，屬于太行山系的一支，地域范圍包括五臺縣全境、繁峙縣南山區(qū)、原平市東山區(qū)、定襄縣東北山區(qū)、盂縣北山區(qū)和河北省阜平縣西山區(qū)，面積約有6 530 km2，北臺頂為華北地區(qū)最高峰，海拔3 061 m[8]。五臺山是首批國家級風景名勝區(qū)、首批國家森林公園、首批5A級景區(qū)[9]，還是我國最早建立的2個草地類的保護區(qū)之一。本研究以五臺山5座臺峰及其山體延伸范圍內(nèi)的山區(qū)為研究對象，研究區(qū)東西最長距離約為100 km，南北最長距離約為80 km，總面積約為6萬hm2。區(qū)域內(nèi)山體多為土石山區(qū)，氣候?qū)儆谂瘻貛駶櫚霛駶櫄夂?，臺頂2 600 m以上屬于高寒氣候[10]。全年平均氣溫為-5℃，極端最低氣溫可達-39.5℃，極端最高氣溫只有20℃。平均年降水量828.5 mm，主要集中在夏季，區(qū)域降雨量受地形影響較大，南坡多于北坡，全年無霜期僅有72 d左右，冬天長達270 d。五臺山的森林多分布于海拔2 000 m左右的山坡上，主要喬木有落葉松、油松、青楊、白樺、云杉等[11]，林帶上下分布著亞高山草甸、林緣草地、灌叢草地和山地草原[12]。

1.2 數(shù)據(jù)來源

本研究的主要數(shù)據(jù)為6期TM影像，影像獲取時間分別為1990年9月16日，1995年6月26日，2000年5月6日，2005年6月21日，2010年5月2日，2015年7月3日，研究中還用到了空間分辨率為30 m的DEM數(shù)據(jù)。研究中用ENVI 5.1進行影像的預處理及植被覆蓋度的計算，用ArcGIS 10.0進行地形分析以及制圖工作。

2 研究方法

2.1 研究數(shù)據(jù)預處理

本研究數(shù)據(jù)的預處理包括對遙感影像的預處理和對多源數(shù)據(jù)投影坐標的變換2個部分。首先需要對6期TM影像進行大氣校正，利用地形圖對其進行幾何校正；然后以TM影像為參照，對DEM數(shù)據(jù)進行坐標投影的變換；最后用研究區(qū)的矢量邊界對TM影像和DEM數(shù)據(jù)進行裁剪。

2.2 研究區(qū)植被覆蓋度估算

目前運用遙感技術(shù)提取區(qū)域內(nèi)植被覆蓋信息的方法主要分為經(jīng)驗模型法、植被指數(shù)法以及混合像元分解法[13-15]。本研究運用植被指數(shù)法對1990—2015年的6期TM影像進行植被覆蓋信息提取。選用歸一化植被指數(shù)（NDVI）對植被覆蓋度進行估算，歸一化植被指數(shù)不僅能部分消除地形、大氣、陰影的影響，而且能有效區(qū)分植被與巖石、裸土，是植被覆蓋度的最佳指示因子[16]。在計算歸一化植被指數(shù)時，選用對綠色植物吸收較強的紅波段及對綠色植物高反射、高透射的近紅外波段，即TM影像中的3波段與4波段。歸一化植被指數(shù)的計算公式如下。

式中，DNNIR，DNR分別為TM影像中紅波段與近紅外波段的反射率，即3波段與4波段的反射率。

目前，利用歸一化植被指數(shù)估算植被覆蓋度最常用的方法為像元二分模型法，也是利用遙感影像估算大面積植被覆蓋狀況最有效的途徑之一[17]。該模型假設(shè)影像每個像元的信息都是由植被信息和土壤信息2個部分組成的，則植被信息占像元的面積比例就是該像元的植被覆蓋度f[18]。設(shè)裸土的像元信息為Ssoil，純植被覆蓋區(qū)域的像元信息為Sveg，則觀測到的影像信息S表示如下。

整個研究區(qū)的植被覆蓋度F計算公式如下。

使用歸一化植被指數(shù)計算植被覆蓋度時取累計頻率為5%的NDVI值為NDVIsoil，取累計頻率為95%的NDVI值為NDVIveg。由此得出植被覆蓋度的估算公式如下。

2.3 研究區(qū)植被覆蓋度變化趨勢分析

利用線性傾向趨勢估計對1990—2015年間研究區(qū)的植被覆蓋度進行趨勢分析，即選取年平均植被覆蓋度值在像元尺度進行一元線性回歸，通過對植被覆蓋度均值的回歸運算分析研究區(qū)植被覆蓋度的變化趨勢[19]。

式中，Sslope為回歸斜率，n為監(jiān)測年數(shù)；Fi為第i年的植被覆蓋度值。

回歸斜率可以反映出1990—2015年植被覆蓋度的變化趨勢與幅度，當斜率大于0時，表明植被覆蓋度呈上升趨勢，當斜率小于0時，表明植被覆蓋度呈下降趨勢，且斜率的絕對值越大，變化幅度越大?；谠撔甭士梢杂嬎阊芯繀^(qū)內(nèi)的植被覆蓋度變化率，計算公式如下。

植被覆蓋度變化率＝Sslope/植被覆蓋度均值×10×100% （6）

3 結(jié)果與分析

3.1 研究區(qū)植被覆蓋度分級結(jié)果

以研究區(qū)1990—2015年生長季的TM-NDVI數(shù)據(jù)為基礎(chǔ)，利用像元二分模型估算1990—2015年五臺山核心區(qū)的植被覆蓋度。根據(jù)水利部2008年頒布的《土壤侵蝕分類分級標準（SL 190—2007）》，將不同的水土流失等級對應不同的植被覆蓋度[20]。結(jié)合五臺山地區(qū)的實際情況，將植被覆蓋度大小分為5級：＜30%（低植被覆蓋度）、30%～45%（中低植被覆蓋度）、45%～60%（中等植被覆蓋度）、60%～75%（中高植被覆蓋度）和＞75%（高植被覆蓋度）。分級結(jié)果如圖1所示。

從圖1可以看出，五臺山地區(qū)的植被覆蓋度分級分布有較強的空間規(guī)律性，總體呈現(xiàn)中間高、四周低的趨勢。這與五臺山地區(qū)的地形狀況有關(guān)，中間為山區(qū)，森林、草甸面積大，受人類活動擾動較小；周邊為居民聚集區(qū)，森林、草甸面積小，受人類活動擾動大。從各年份的植被覆蓋度情況來看，1990，1995年研究區(qū)中的山區(qū)多為高植被覆蓋度區(qū)域，2000，2005，2010年研究區(qū)以中高分辨率區(qū)域為主，但中低、低分辨率區(qū)域逐漸增多，到2015年時研究區(qū)內(nèi)的中低、低分辨率區(qū)域達到最多，高分辨率區(qū)域也是6期影像中最少的。

3.2 研究區(qū)植被覆蓋度變化

表1 1990—2015年植被覆蓋度平均值

表1為6期影像植被覆蓋度平均值的統(tǒng)計結(jié)果，可以看出，2015年與1990年相比，植被覆蓋度平均值下降約21百分點，2000—2010年植被覆蓋度變化最明顯，10 a約降低了15.06%。從圖2可以看出，研究區(qū)的植被覆蓋度呈現(xiàn)不斷下降的趨勢，且2000年以后的下降速度明顯加快。基于表1中的數(shù)據(jù)對其進行變化趨勢分析，得出回歸斜率約為-0.009，印證了研究區(qū)域植被覆蓋度在這25 a中有明顯的降低趨勢，其植被覆蓋度變化率約為15.3%。

遙感變化檢測是通過比值或差值運算得到不同影像的變化區(qū)域及變化狀況的方法，圖3為1990—2015年間植被覆蓋度變化檢測結(jié)果。從1990，2015年的變化檢測結(jié)果可以看出，研究區(qū)植被覆蓋度在25 a中變化明顯，大部分區(qū)域都有不同程度的降低。西臺頂西側(cè)、南臺頂南側(cè)降低的最為明顯，北臺頂北側(cè)的變化相對較小，東臺頂東側(cè)的變化最小。從各時間階段的變化來看，2000—2005年研究區(qū)植被覆蓋度的降低幅度最大，變化區(qū)域也最為廣泛；1990—1995年降低幅度最小，而且沒有變化量明顯較高的區(qū)域。從各區(qū)域變化時間來看，南臺頂南側(cè)的植被破壞主要發(fā)生在2000—2005年；西臺頂西側(cè)與北臺頂北側(cè)的變化主要發(fā)生在2010—2015年，其他區(qū)域的植被破壞在各時間段比較平均，沒有植被覆蓋度快速降低的明顯時間節(jié)點。另外，從圖3可以看出，這25 a間以臺懷鎮(zhèn)為中心的五臺山旅游區(qū)植被覆蓋度降低幅度也較大，大部分區(qū)域的降低幅度在30%以上，幾座臺峰的變化尤為明顯。

3.3 人為活動對區(qū)域植被的影響分析

經(jīng)過對研究區(qū)的實地走訪、調(diào)查，發(fā)現(xiàn)造成該區(qū)域植被破壞的原因主要有3點：（1）五臺山旅游區(qū)內(nèi)旅游開發(fā)對高山植被的破壞；（2）研究區(qū)內(nèi)的礦山開采活動對林地和草地的破壞；（3）旅游開發(fā)和礦山開發(fā)伴隨的道路建設(shè)對林地和草地的破壞。

旅游開發(fā)的影響區(qū)域為臺懷鎮(zhèn)周邊，主要為5座臺峰，造成植被破壞的原因包括寺廟建設(shè)、停車場等旅游基礎(chǔ)設(shè)施的建設(shè)及頻繁的旅游活動。雖然這些建筑的單個面積和總面積都較小，但多分布于亞高山草甸中，由于亞高山草甸的生態(tài)比較脆弱，破壞以后的恢復周期較長，因此，旅游開發(fā)可能會導致5座臺峰臺頂及周邊海拔較高區(qū)域的草甸造成難以修復的破壞。近年來，五臺山的旅游開發(fā)力度越來越大，5個臺頂?shù)乃聫R都在擴建，通往臺頂?shù)牡缆芬苍诓粩鄶U寬，這些開發(fā)活動都對臺頂及周邊的亞高山草甸造成了較大的破壞。

礦山開采區(qū)域分布于五臺山風景區(qū)的外圍，但由于開采區(qū)的面積大，開采活動對植被、水體、土壤等自然因素都有嚴重影響，因此，礦山開采對五臺山的生態(tài)環(huán)境影響最深遠。五臺山周邊的礦山開采區(qū)距離臺懷鎮(zhèn)直線距離在10～50 km，南臺頂南側(cè)的五臺縣境內(nèi)、南臺頂西側(cè)的代縣境內(nèi)以及西臺頂西側(cè)的繁峙縣境內(nèi)是分布最聚集，面積最大的區(qū)域。由于礦山開采區(qū)的分布區(qū)域分散、面積大、海拔范圍大，因此，對區(qū)域內(nèi)的林地和草甸影響都較大。

道路建設(shè)的目的是為旅游開發(fā)、礦山開采及周邊居民服務，因此，多分布于旅游區(qū)、居民區(qū)和開采區(qū)的周邊。旅游區(qū)中的道路寬度約5～15 m，海拔越高，道路的寬度越窄。這一部分道路主要影響5座臺峰的植被，由于臺頂部分區(qū)域的土質(zhì)較為稀松，因此，在這些區(qū)域道路的破壞區(qū)域遠超道路本身范圍。采礦區(qū)周邊的道路，由于要滿足大型貨車的運行，且建設(shè)過程不規(guī)范，因此，寬度在5～30 m范圍內(nèi)不等。這部分道路破壞的植被主要為五臺山周邊山區(qū)的林地和部分海拔較高區(qū)域的草地。在調(diào)查過程中還發(fā)現(xiàn)，五臺山周邊的風力發(fā)電站建設(shè)過程中的道路建設(shè)也造成了較大面積的植被破壞。風力發(fā)電站主要分布于西臺、中臺、北臺西側(cè)部分海拔較高的山體頂部，這些區(qū)域的植被多為草甸，這些區(qū)域的道路影響的植被也多為草甸。道路本身的面積雖然不大，但由于它本身的分布范圍廣、海拔跨度大，且道路周邊的人類活動頻繁，因此，對植被影響較大。

4 結(jié)論與討論

1990—2015年，研究區(qū)的植被覆蓋度降低約21%，降低趨勢明顯，其中在2000—2010年植被覆蓋度降幅最大，年平均降幅達到1.5%。植被覆蓋度變化最大的區(qū)域為南臺西南側(cè)與西臺西側(cè)，大部分區(qū)域的降幅都達到了30%以上。野外實地調(diào)查發(fā)現(xiàn)，南臺西南側(cè)和西臺西側(cè)都分布有大面積的礦山開采破壞區(qū)，這些區(qū)域的植被覆蓋度大幅降低可能與礦山開采活動有關(guān)。五臺山的5座臺峰位于五臺山草地自然保護區(qū)的范圍內(nèi)，自1990年以來該區(qū)域的植被覆蓋度也明顯降低，且近年來降幅不斷增大，5個臺頂?shù)穆糜伍_發(fā)活動可能是造成這一變化的重要原因。此外，旅游開發(fā)、礦山開采過程中伴隨的道路建設(shè)也對五臺山地區(qū)的植被有較大的影響。綜上所述，五臺山地區(qū)的人為開發(fā)、開采活動是造成該區(qū)域近年來植被覆蓋度持續(xù)降低的重要原因。

針對上述的五臺山植被變化現(xiàn)狀及其原因，特提出以下建議。在五臺山的5座臺峰高海拔區(qū)域發(fā)展生態(tài)旅游，減少人為旅游開發(fā)活動對草甸的影響；加強對該區(qū)域，尤其是亞高山草甸區(qū)域道路建設(shè)的規(guī)劃，盡量避免不必要的道路建設(shè)，減少道路建設(shè)對植被的影響；加強對礦山開采區(qū)的規(guī)劃和管理，避免不規(guī)范甚至違法的開采活動；加強對開采區(qū)的生態(tài)恢復工作，避免對區(qū)域生態(tài)造成不可逆的破壞；加大對五臺山區(qū)域植被及生態(tài)的監(jiān)測力度，保障對區(qū)域植被變化信息的及時獲取。

參考文獻：

［1］樊文華，郭先龍，池寶亮，等.五臺山草地自然保護區(qū)草地資源的開發(fā)利用[J].中國草地，1999（2）：13-16.

［2］李偉濤.高分辨率遙感森林植被分類提取研究[D].北京：北京林業(yè)大學，2016.

［3］王坎，唐瑤，王玲玲.基于遙感技術(shù)的十八里長峽自然保護區(qū)保護效果監(jiān)測與評估[J].華中師范大學學報，2015，49（6）：929-935.

［4］ BURNETT C，F(xiàn)ALL A，TOMPPO E.Monitoring current status of and trends in boreal forest landuse in Russian KARELIA[J].Conservation Ecology，2003，7（2）：1850-1851.

［5］JUSOFF K，SETIAWANI.Quantifying deforestation in a permanent forest reserve using vectorised Landset TM[J].Journal of Tropical Forest Science ，2003，15（4）：570-582.

［6］劉文軍.基于遙感技術(shù)的土地利用/土地覆蓋現(xiàn)狀信息提取研究：以艾比湖濕地保護區(qū)為例[J].新疆環(huán)境保護，2007（3）：1-4.

［7］白振平，劉洪利.霧靈山植被變化遙感監(jiān)測[J].首都師范大學學報，2003，24（4）：59-62.

［8］張利利，陳建強.山西五臺山風景區(qū)生態(tài)旅游構(gòu)想[J].安徽農(nóng)業(yè)科學，2011，39（5）：2903-2905，2913.

［9］程占紅，趙蒙.旅游干擾對五臺山植被景觀區(qū)植物多樣性的影響[J].干旱區(qū)地理，2012，35（4）：578-586.

［10］張建彪，閆美芳，上官鐵梁.五臺山亞高山草甸的β多樣性研究[J].西北植物學報，2006，26（2）：389-392.

［11］劉宏芳，袁建英.五臺山植被景觀格局分析[J].山西農(nóng)業(yè)大學學報，2012，32（5）：437-441.

［12］袁建英.五臺山地區(qū)植被景觀生態(tài)研究[D].太原：山西大學，2001.

［13］牛寶茹，劉俊蓉，王政偉.干旱區(qū)植被覆蓋度提取模型的建立[J].地球信息科學，2005，7（1）：84-97.

［14］王宇琛，李松鳴，王霄，等.汾河流域2001—2013年植被覆蓋時空變化及其對氣候因子的響應 [J].山西農(nóng)業(yè)科學，2016，44（5）：640-645.

［15］周湘山，孫保平，李錦榮，等.內(nèi)蒙古卓資縣土地利用與植被覆蓋動態(tài)變化[J].山西農(nóng)業(yè)科學，2011，39（9）：981-986.

［16］王國芳.基于遙感技術(shù)的植被覆蓋度的動態(tài)監(jiān)測[J].山西農(nóng)業(yè)科學，2015，43（5）：592-595.

［17］劉宇.基于遙感的若爾蓋濕地保護區(qū)沙化草地時空變化特征分析[D].雅安：四川農(nóng)業(yè)大學，2014.

［18］任唯敏.植被指數(shù)計算區(qū)域植被覆蓋度的適用性研究 [D].西安：西北大學，2012.

［19］俞世榮，郭藹平.農(nóng)業(yè)研究中的多元回歸分析和矩陣淺介[J].天津農(nóng)業(yè)科學，1983（1）：34-44.

［20］詹起林.基于森林植被指數(shù)的快速生態(tài)評價：以九寨溝景區(qū)為例[D].上海：上海師范大學，2012.