拉市海省級自然保護區生態質量10年變化研究

段衛虎 王一凱 黃 誠 王艷霞 周汝良

（西南林業大學，云南 昆明 650224）

拉市海省級自然保護區位于中國西南部，云南的西北部，是由天然高原湖泊拉市海、文海、文筆水庫和吉子水庫四片組成，其地貌上屬于橫斷山系高山峽谷區。拉市海省級自然保護區類型為內陸高原濕地生態系統，雖然面積不大，但由于其特殊的地理位置和漫長的形成育化過程，孕育了豐富的的動植物資源。拉市海自然保護區中植物共有566種，分為濕地植物和陸生植物。濕地植物有150種，隸屬于39科115屬；陸生植物416種，隸屬于93科253屬。拉市海保護區中共有保護植物11種，其中國家一級保護植物2種，二級保護植物7種。根據中國動物地理區劃，拉市海保護區的動物區系屬東洋界印亞界，西南區西南山地亞區。有鳥類232種，淡水魚類25種，兩棲類14種，爬行類17種。這些動物中，國家一級保護對象8種、國家二級保護對象30種[1-3]。面對全球變暖，人類現代化進程加快，拉市海自然保護區的保護由于特殊地質地貌和脆弱水文環境面臨十分嚴峻的問題。因此，本文通過對拉市海省級自然保護區的十年生態質量變化進行定量化研究，為該保護區自然資源的合理開發利用和生物多樣性保護提供科學的依據[4]。

1 研究區概況

拉市海自然保護區于1998年云南省政府批準成立了麗江拉市海高原濕地省級自然保護區，2002年成立云南麗江拉市海高原濕地自然保護區管理局，其范圍在行政上屬麗江市玉龍納西族自治縣的拉市、白沙、黃山、太安四個鄉（鎮）所轄，包括四個片區。其中，拉市海片區面積為面積5015.0hm2，文海片區面積840.0hm2，吉子水庫片區面積476.0hm2，文筆水庫片區面積192.0hm2。

2 實驗數據及方法

2.1 實驗數據

拉市海自然保護區生態系統質量評估主要利用遙感解譯獲取的2000-2010年逐旬生態系統地表參量柵格數據，包括地上葉面積指數、植被覆蓋度、凈初級生產力，如下表（表1），矢量數據包括拉市海自然保護區生態類型分布圖、保護區范圍邊界圖。

表1 拉市海自然保護區柵格源數據

2.2 算法介紹

本文針對拉市海自然保護區中森林生態系統、草地生態系統、濕地生態系統，分別采用葉面積指數（Leaf Area Index，LAI）、植被覆蓋度（Fraction Vegetation Cover，FVC）、初級凈生產力（Net Primary Production,NPP）及其對應的年均變異指數作為生態質量指標來進行動態評估。

(1)生態指標年均指數計算

森林生態系統年均葉面指數、草地生態系統年均植被覆蓋度、濕地生態系統年均初級凈生產力的計算均按按公式1進行計算：

式中，i為年數，j為旬數，DecLij為第 i年第 j旬影像 LAI、FVC、NPP值。

(2)生態指標年均變異系數

式中，i為年數，j為旬數，DecLij分別對應森林、草地、濕地生態系統的第i年第j旬影像LAI、FVC、NPP值，n分別對應森林、草地、濕地生態系統內影像像元數量。

3 結果與分析

3.1 森林生態系統質量評估

葉面積指數（LAI）：單位土地面積上植物葉片總面積占土地面積的倍數。葉片是植物進行光合作用和與外界進行水汽交換的的主要器官，葉片面積的大小直接影響植物能量合成和物質轉換。LAI是森林生態系統的最重要的參數之一，它不僅被用來估測林分尺度、作為生態過程的主要參數和植物冠層功能過程，而且隨著遙感技術的發展，它被用作進行大面積森林生態系統生產力的估測指標[5－6]。

在本文中，拉市海自然保護區的森林生態系統質量評估采用葉面積指數，數據源為2000～2010年的十年,每年36旬的250m的葉面積指數柵格數據。按公式1，計算拉市海自然保護區的2000～2010年的各年年均葉面積指數。然后將葉面積指數按0-2（低）、2-4（較低）、4-6（中）、6-8（較高）、8-∞（高）五個等級進行分類重采樣。最后統計個等級葉面積指數在拉市海省級自然保護區中所占的面積（如表2）。

從表2分析可以得到，在十年的拉市海省級自然保護區的森林生態系統中，葉面積指數指數處于高等級的占到森林生態系統面積的90%左右，其次為較高等級、低等級，低等級和中等級所占的面積較小。由表2可大略得到，拉市海省級自然保護區森林生態系統的十年的葉面積變化不明顯，為了準確評價其變化程度，引入森林生態系統葉面積指數年變異系數、年均變異系數。

表2 森林生態系統質量年均葉面積指數各等級面積(km2)

表3 森林生態系統葉面積指數年變異系數各等級面積(km2)

圖1 2000-2010年森林生態系統葉面積指數年均變異系數

表1的數據，是按公式2計算森林生態系統葉面積指數年變異系數， 然后按 0-0.2、0.2-0.4、0.4-0.8、0.8-1、1- ∞ 五個等級進行劃分，然后統計各等級變異系數在森林生態系統中所占的面積。結合表1、表2可知，拉市海自然保護區十年的森林生態系統90%面積處于葉面積指數等級高狀態，森林生態系統99%面積的葉面積指數處于變異系數等級的低和較低狀態。由此可知，拉市海自然保護區十年的森林生態系統的質量較好，十年的變化較小。

3.2 草地生態系統質量評估

植被覆蓋度（Fractional Vegetation Cover,簡稱FVC）：植物群落總體或各個體的地上部份的垂直投影面積與樣方面積之比的百分數。反映植被的茂密程度和植物進行光合作用面積的大小。同時，FVC在生態系統的氣候模型中為重要的氣候參數，在水文模型中的空間分布和時間變化動態是能量和水循環的影響因子[7-8]，并且也可以作為土壤侵蝕和水土流失預報的重要因子。FVC適用于草原、草地等生態敏感性較強的生態系統。

本研究采用年均植被覆蓋度及其年均變異系數進行拉市海省級自然保護區草地生態系統質量時空動態評估。數據源為拉市海自然保護區的草地生態系統的2000～2010年，每年36旬的植被覆蓋度柵格數據。先按公式1計算拉市海自然保護區十年的草地生態系統年均植被覆蓋度，再按照公式2、3計算十年的草地生態系統年均變異系數。

表4中的數據是將拉市海自然保護區的濕地生態系統的年均植被覆蓋度按0-20%、20%-40%,40%-60%,60%-80%,80%-100%分為五個等級，然后統計每個等級在濕地生態系統所占的面積。圖2為濕地生態系統植被覆蓋度的年均變異系數，該系數反映了拉市海草地生態系統的質量波動程度。

表4 草地生態系統年均植被覆蓋度各等級面積（km2）

圖2 2000-2010年草地生態系統植被覆蓋度年均變異系數

由表4、圖2分析可知，拉市海自然保護區草地生態系統90%面積的植被覆蓋度質量處于中級、較高級狀態，其次分布在較低級、低級狀態，最后是處于高級狀態，但是高級面積非常小，3%左右。拉市海自然保護區濕地生態系統植被覆蓋度中級、較高級所占的面積，都有大幅度的減小，較高級面積由2000年的83.4375km2減小為2010年的41.0000km2，中級所占的面積有2000年的33.3750km2減小為2010年的8.3750km2，其他等級沒有明顯的變化。

3.3 濕地生態系統質量評估

凈初級生產力則是由光合作用所產生的有機質總量中扣除自養呼吸后的剩余部分。是生產者能用于生長、發育和繁殖的能量值，反映了植物固定和轉化光合產物的效率，也是生態系統中其他生物成員生存和繁衍的物質基礎。反映了植物每年通過光合作用所固定的碳總量，目前NPP研究是全球變化研究的重要內容之一[9-10]。

本文利用凈初級生產力進行拉市海自然保護區濕地生態系統質量評估。本文進行評估的數據源為2000～2010年，每年36旬的凈初級生產力柵格數據，每旬數據為一旬的累積值，單位為0.1g/m2。處理過程為：首先，計算濕地生態系統的年均凈初級生產力，計算公式為公式1。其次，計算濕地生態系統的凈初級生產力年總量，計算公式為公式4。最后，評估濕地生態系統的初級生產力年變化情況，進行年均變異系數計算，計算公式為公式2，公式3。

式中，NT_DecLijk為第i年第j旬影像中第k個像元NPP值，Sk為第k個像元面積。

表5 濕地生態系統年均凈初級生產力各等級面積（km2）

表5中的數據是將計算的得到的SD_AUNi按專家劃分的等級標準分為低、較低、中、較高、高五級，每級對應取值為0-6,6-12,12-18,18-24,24-∞。圖3、圖4為計算得到的濕地生態系統十年的各年凈初級生產力年總量和十年各年的初級生產力年均變異系數。

由表5分析可知，拉市海濕地生態系統面積在4.25～8.6875km2之間，面積較小，并且面積從2000到2010年過程中，濕地面積減小2.9375 km2。但是十年的濕地生態質量均處于高等級良好的狀態。由圖3分析可知，2000到2010年，拉市海自然保護區的濕地生態系統凈初級生產力年總量有了較明顯的上升。由圖4分析可知，拉市海濕地生態系統的質量的波動程度較小，但是結合森林、草地的波動情況，濕地生態系統的質量波動最大。

圖3 2000-2010年濕地生態系統凈初級生產力年總量（t）

圖4 2000-2010年濕地生態系統初級生產力年均變異系數

4 結論與討論

拉市海自然保護區處于青藏高原東南邊陲，古老的地質運動形成特殊的地形地貌、水文環境，這樣的特殊的環境下孕育著世界和我國許多珍惜的瀕危的野生動植物。拉市海自然保護區又作為一個內陸濕地保護區，具有換水周期長，生態系統較脆弱的特點，但是同時面對全球變暖、人類活動干擾強度的急升，我國許多河流、湖泊出現明顯的水量萎縮、季節斷流、甚至干枯的嚴峻形勢。因此，對該保護區生態質量的十年遙感動態評估，為深刻而全面了解拉市海自然保護區近十年的生態系統時空變化奠定基礎。

本文分別針對森林生態系統、草地生態系統、濕地生態系統，應用葉面積指數、植被覆蓋度、凈初級生產力及其變異系數指標，對拉市海省級自然保護區中的主要生態林類型進行針對性、量化的動態評估，為拉市海自然保護區的管理、保護、開發與利用提供更加科學的依據。

在評估的過程中發現，拉市海自然保護區的十年整體生物量有所增加；森林生態系統質量處于高等級，十年沒有明顯變化；草地生態系統處于中、較高等級，但是相對原來面積，明顯的減小。濕地生態系統所占的面積較小，但是面積有所減小。因此，加強對拉市海自然保護區的保護，尤其是草地和濕地生態系統，為野生動植物的生存提供穩定而良好的棲息環境，同時也有利于保護區及周邊居民區的小氣候的良性的可持續發展。

［1］云南植被編寫組.云南植被[M].北京：科學出版社，1987.

［2］王應祥.中國哺乳動物種和亞種分類名錄與分布大全[M].北京：中國林業出版社，2003.

［3］鄭光美.中國鳥類分類與分布名錄[M].北京：科學出版社，2005.

［4］楊國斌.蘭坪云嶺省級自然保護區質量評價研究[J].西部林業科學，2011,40(4)：48-53.

［5］巨文珍，農勝奇.森林生物量研究進展[J].西南林業大學學報，2011,31（2）：78-83.

［6］王希群，馬履一，賈忠奎，徐程揚.葉面積指數的研究和應用進展[J].生態學雜志，2005,24（5）：537-541.

［7］Qi J,Marsett R S,et al.Spatial and temporal dynamics of vegetation in the San Pedro river basin area[J].Agricultural and Forest Mateorology ，2000，105(1)：55-68.

［8］姚云軍，程潔，趙少華，等.基于熱紅外遙感的農田蒸散估算方法研究綜述[J].地球科學進展，2012,27(2)：1308-1318.

［9］王新闖，王世東，張合并.基于MOD17A3的河南省NPP時空格局[J].生態學雜志，2013,32（10）：2797-2805.

［10］張杰，潘曉玲.天山北麓山地-綠洲-荒漠生態系統凈初級生產力空間分布格局及其季節變化[J].干旱區地理，33(1)：78-86.