宋偉宏，程慧波

(甘肅省環(huán)境監(jiān)測(cè)中心站，甘肅蘭州 730000)

土地利用/覆蓋變化是全球環(huán)境變化研究的前沿和熱點(diǎn)問(wèn)題[1-3]。隨著土地利用/覆被變化研究的深入，政府、基層組織和居民對(duì)生態(tài)環(huán)境變化的響應(yīng)成為研究熱點(diǎn)[4-6]。由于對(duì)生態(tài)環(huán)境變化的響應(yīng)關(guān)系著區(qū)域的可持續(xù)能力，在一些脆弱區(qū)域，當(dāng)基層組織和居民對(duì)生態(tài)環(huán)境有著消極響應(yīng)的時(shí)候，政府的積極響應(yīng)對(duì)于區(qū)域的可持續(xù)發(fā)展尤為重要[7-10]。甘肅祁連山國(guó)家級(jí)自然保護(hù)區(qū)地處青藏、蒙新、黃土三大高原交匯地帶的祁連山北麓，屬森林生態(tài)系統(tǒng)類(lèi)型的自然保護(hù)區(qū)。祁連山是青海、甘肅兩省的天然分界線，分布有豐富的冰川、雪山、森林、草地和濕地資源，是我國(guó)高原生態(tài)安全屏障的重要組成部分，保障著全國(guó)，特別是西北地區(qū)的生態(tài)安全，在維護(hù)整個(gè)河西陸地生態(tài)系統(tǒng)平衡以及保障河西地區(qū)經(jīng)濟(jì)社會(huì)和諧發(fā)展中起著決定性的作用[11-13]。近年來(lái)，受自然因素和人為因素的影響，祁連山地區(qū)冰川、雪線不斷退縮，天然林草植被不斷退化，水源涵養(yǎng)功能呈現(xiàn)萎縮，荒漠化日趨擴(kuò)大，生物多樣性遭到威脅。開(kāi)展祁連山自然保護(hù)區(qū)生態(tài)環(huán)境監(jiān)測(cè)對(duì)保護(hù)祁連山生態(tài)環(huán)境具有重要的意義。

1 資料與方法

1.1研究區(qū)概況甘肅祁連山國(guó)家級(jí)自然保護(hù)區(qū)位于祁連山北坡中、東段，地理位置位于97°25′～103°46′E、36°43′～39°36′N(xiāo)，南沿祁連山主脈與青海省接壤，西至肅南縣界與肅北蒙古族自治縣相鄰，東至天祝縣界與永登縣相連，北至祁連山森林分布下線與河西走廊相鄰。祁連山國(guó)家級(jí)自然保護(hù)區(qū)屬大陸性高寒半濕潤(rùn)山地氣候，表現(xiàn)為冬季長(zhǎng)而寒冷干燥，夏季短而溫涼濕潤(rùn)，全年降水量主要集中在5—9月，該區(qū)由淺山地帶向深山地帶氣溫遞減、雨量遞增，高山寒冷而陰濕，淺山地帶熱而干燥。隨著山區(qū)海拔的升高，各氣候要素發(fā)生有規(guī)律自下而上的變化，呈明顯的山地垂直氣候帶，自下而上為淺山荒漠草原氣候帶、淺山干草原氣候帶、中山森林草原氣候帶、亞高山灌叢草甸氣候帶、高山冰雪植被氣候帶。

1.2數(shù)據(jù)與方法

1.2.1數(shù)據(jù)來(lái)源與處理。為建立2000—2016年甘肅省祁連山自然保護(hù)區(qū)時(shí)空數(shù)據(jù)庫(kù)，從USGS網(wǎng)站下載了102景Landsat TM/ETM+/OLI遙感影像，作為甘肅祁連山自然保護(hù)區(qū)的基礎(chǔ)研究數(shù)據(jù)。EVI數(shù)據(jù)來(lái)源美國(guó)航天局NASA數(shù)據(jù)共享網(wǎng)(https://ladsweb.modaps.eosdis.nasa.gov/search/)。通過(guò)植被指數(shù)算法獲得，該算法基于質(zhì)量、云和觀測(cè)幾何濾波器進(jìn)行數(shù)據(jù)濾波，只有高質(zhì)量的、無(wú)云的濾波數(shù)據(jù)被用來(lái)進(jìn)行合成，合成的目的是選擇最佳的觀測(cè)數(shù)據(jù)，這種方法能保證植被指數(shù)值的時(shí)空一致性。

植被覆蓋度是基于植被指數(shù)利用像元二分法模型計(jì)算的，計(jì)算公式為：

VFC=(EVI-EVIsoil)/(EVIveg-EVIsoil)

(1)

其中， EVIsoil為完全是裸土或無(wú)植被覆蓋區(qū)域的EVI值，EVIveg則代表完全被植被所覆蓋的像元的EVI值，即純植被像元的EVI值。

1.2.2研究方法。采用決策樹(shù)分析方法，通過(guò)采用人工與自動(dòng)相結(jié)合的方式，對(duì)于光譜劃分機(jī)理清楚的類(lèi)型采用人工建樹(shù)方法，對(duì)于類(lèi)型的光譜變化比較大、規(guī)律不清楚的類(lèi)型采用自動(dòng)方法(最鄰近方法)，并且通過(guò)地面核查提高土地覆蓋分類(lèi)精度。考慮到土地覆被人工目視解譯自身存在的誤差影像，為科學(xué)評(píng)判土地覆被變化狀況，采用以下公式提取土地覆被面積變化狀態(tài)：

ΔS=Sj-Si

式中，ΔS為兩期地物類(lèi)型面積的插值，i和j均為年份；ΔE為兩期地物類(lèi)型矢量化數(shù)據(jù)誤差；Pi和Pj分別為i和j時(shí)刻地物類(lèi)型多邊形的周長(zhǎng)；Hc為遙感影像像元的一半，取值為15 m。

2 結(jié)果與分析

2.1祁連山自然保護(hù)區(qū)土地利用/覆蓋以2016年度覆蓋全甘肅省祁連山國(guó)家級(jí)自然保護(hù)區(qū)的Landsat8 OLI影像為主要數(shù)據(jù)源，通過(guò)室內(nèi)解譯并結(jié)合野外核查，得出2016年度甘肅省祁連山國(guó)家級(jí)自然保護(hù)區(qū)土地利用/土地覆蓋狀況(圖1)。祁連山獨(dú)特的自然條件和人為活動(dòng)的影響，決定了各生態(tài)系統(tǒng)的構(gòu)成及面積相差較大，如表1所示。2016年整個(gè)保護(hù)區(qū)生態(tài)系統(tǒng)組成以草地為主，占總面積的52.9%；其次是荒漠，占22.4%；森林與荒漠接近，占21.9%；農(nóng)田占1.2%，冰川積雪占1.1%，濕地占0.3%；人工生態(tài)系統(tǒng)最小，僅占0.2%。

圖1 2016年祁連山國(guó)家級(jí)自然保護(hù)區(qū)土地覆被現(xiàn)狀Fig.1 Status of land cover in Qilian Mountain National Nature Reserve in 2016

從圖2可看出，森林主要以落葉闊葉灌叢為主，面積473 933.6 hm2，占森林面積的80.6%。落葉闊葉灌叢主要分布在實(shí)驗(yàn)區(qū)，面積占該類(lèi)林地總面積的44.4%；其次分布在核心區(qū)，占26.2%；分布在外圍保護(hù)地帶的面積最少，僅占8.5%。常綠針葉林面積104 733.4 hm2，占森林面積的17.8%。常綠針葉林主要分布在實(shí)驗(yàn)區(qū)，占該類(lèi)林地面積的64.7%；緩沖區(qū)占20.6%，外圍保護(hù)地帶最小，僅占1.2%。其他林地面積很小，主要零星分布在自然保護(hù)區(qū)的實(shí)驗(yàn)區(qū)、外圍保護(hù)地帶。

從圖3可看出，草地構(gòu)成中沒(méi)有占主導(dǎo)地位的類(lèi)型。溫性草原最多，面積386 154.0 hm2，也僅占草地總面積的20.2%；溫性草原主要分布在實(shí)驗(yàn)區(qū)，占該類(lèi)草地面積的55.0%，分布在外圍保護(hù)地帶的占36.0%，分布在核心區(qū)的最少，僅占4.1%。高寒草原面積351 879.9 hm2，主要分布在外圍保護(hù)地帶，占該類(lèi)草地面積的37.7%，分布在實(shí)驗(yàn)區(qū)的占22.1%，分布在核心區(qū)的占20.5%，分布在緩沖區(qū)的占19.7%。溫性草甸面積286 814.3 hm2，主要分布在實(shí)驗(yàn)區(qū)，占該類(lèi)草地面積的51.6%；分布在外圍保護(hù)地帶的占22.5%，分布在緩沖區(qū)的占14.3%，分布在核心區(qū)的占11.6%。高寒草甸面積73 735.5 hm2，主要分布在實(shí)驗(yàn)區(qū)，占該類(lèi)草地面積的48.3%；分布在核心區(qū)的占24.5%，分布在緩沖區(qū)的占18.0%；分布在外圍保護(hù)地帶的最少，僅占9.2%。稀疏草地面積318 725.2 hm2，主要分布在外圍保護(hù)地帶，占該類(lèi)草地面積的53.9%，分布在實(shí)驗(yàn)區(qū)的占28.2%，分布在核心區(qū)的占10.0%；分布在緩沖區(qū)的最少，占7.9%。

圖2 2016年祁連山國(guó)家級(jí)自然保護(hù)區(qū)森林覆被現(xiàn)狀Fig.2 Status of forest cover in Qilian Mountain National Nature Reserve in 2016

圖3 2016年祁連山國(guó)家級(jí)自然保護(hù)區(qū)草地覆被現(xiàn)狀Fig.3 Status of grassland cover in Qilian Mountain National Nature Reserve in 2016

人工生態(tài)系統(tǒng)主要是指由人為集中活動(dòng)區(qū)域而強(qiáng)烈改變了地表自然狀態(tài)的有人為活動(dòng)的區(qū)域。祁連山自然保護(hù)區(qū)人工生態(tài)系統(tǒng)主要是建設(shè)用地，面積1 943.4 hm2，分布在外圍保護(hù)地帶的建設(shè)用地占該類(lèi)總面積的54.8%，實(shí)驗(yàn)區(qū)的占43.7%，緩沖區(qū)的占1.3%，核心區(qū)的占0.2%。交通用地面積377.4 hm2，其中分布在實(shí)驗(yàn)區(qū)的占53.5%，在外圍保護(hù)地帶的占40.4%，緩沖區(qū)的占6.1%，核心區(qū)沒(méi)有明顯的交通用地。采礦場(chǎng)面積583.0 hm2，分布在外圍保護(hù)地帶的占73.0%，在實(shí)驗(yàn)區(qū)的占27.0%。祁連山自然保護(hù)區(qū)土地覆被面積如表1所示。

2.2祁連山保護(hù)區(qū)土地覆被變化近16年來(lái)祁連山國(guó)家級(jí)自然保護(hù)區(qū)主要生態(tài)系統(tǒng)整體處于穩(wěn)定狀態(tài)，草地生態(tài)系統(tǒng)是外圍保護(hù)地帶由于祁連山國(guó)家級(jí)自然保護(hù)區(qū)地處山區(qū)、人為大規(guī)?；顒?dòng)干擾少，其中2000—2005年變化幅度相對(duì)于大于后面2個(gè)時(shí)期(表2)。2000—2005年自然保護(hù)區(qū)整體變化趨勢(shì)是森林和草地生態(tài)系統(tǒng)及冰川積雪面積減少，農(nóng)田、濕地、荒漠和人工生態(tài)系統(tǒng)面積增加(表3)。

2000—2005年森林生態(tài)系統(tǒng)面積凈減少475.9 hm2，減少的森林主要是高山灌叢變?yōu)椴莸?，其次是變?yōu)槁愕?，其他減少變化量很小。草地凈減少839.1 hm2，主要是被開(kāi)墾為耕地，其次是由于植被退化變?yōu)槁愕?，其他減少量很小。農(nóng)田生態(tài)系統(tǒng)凈增加1 635.7 hm2，增加的耕地主要來(lái)源于草地開(kāi)墾，其次是裸土地開(kāi)發(fā)為耕地。濕地凈增加680.2 hm2，增加的濕地主要來(lái)源于荒漠和草地生態(tài)系統(tǒng)。人工生態(tài)系統(tǒng)增加163.1 hm2，主要來(lái)源于荒漠和草地生態(tài)系統(tǒng)?；哪鷳B(tài)系統(tǒng)增加1 257.4 hm2，其增減變化主要是高山裸巖與冰川積雪之間的相互轉(zhuǎn)換。冰川積雪凈減少2 421.4 hm2，其增減變化與荒漠生態(tài)系統(tǒng)相似。

相對(duì)于前5 年，祁連山國(guó)家級(jí)自然保護(hù)區(qū)在2005—2010年變化相對(duì)較小，其中森林、農(nóng)田、濕地、荒漠生態(tài)系統(tǒng)和冰川積雪面積凈減少，草地和人工生態(tài)系統(tǒng)面積凈增加。各生態(tài)系統(tǒng)在此期間變化量小，沒(méi)有顯著變化(表4)。

2010—2016年草地、農(nóng)田和荒漠生態(tài)系統(tǒng)面積凈減少，森林、濕地、人工生態(tài)系統(tǒng)和冰川積雪面積凈增加。在此期間除冰川積雪與高山裸巖之間相互轉(zhuǎn)換的面積較大外，其他類(lèi)型生態(tài)系統(tǒng)變化的面積都較小(表5)。

2.3祁連山保護(hù)區(qū)EVI變化甘肅省祁連山國(guó)家級(jí)自然保護(hù)區(qū)植被覆蓋指數(shù)(EVI)空間分布如圖4所示。紅色-藍(lán)色-紫色代表EVI植被指數(shù)由低往高逐步上升，同色系中顏色越淺代表其值越高。祁連山國(guó)家級(jí)自然保護(hù)區(qū)年EVI總體上呈現(xiàn)片狀分布特征，且存在一定的空間分布差異。EVI高值即植被覆蓋指數(shù)高的區(qū)域主要分布在自然保護(hù)區(qū)的東南區(qū)域，西北區(qū)域EVI值普遍較低。每年EVI植被指數(shù)分布大體一致，說(shuō)明甘肅省祁連山國(guó)家級(jí)自然保護(hù)區(qū)EVI空間變化差異不大。

根據(jù)EVI，計(jì)算得出祁連山自然保護(hù)區(qū)植被覆蓋度。2000—2016年植被覆蓋度變化中，2004年植被覆蓋度最低，為25.98%。整體上植被覆蓋度有上升的態(tài)勢(shì)特征。年植被覆蓋度變化與生長(zhǎng)季植被覆蓋度變化基本一致。祁連山外圍保護(hù)地帶植被覆蓋度在祁連山自然保護(hù)區(qū)的核心區(qū)最低，如圖5所示。

表12016年祁連山國(guó)家級(jí)自然保護(hù)區(qū)土地覆被面積及比例

Table1AreaandproportionoflandcoverofQilianMountainNationalNatureReservein2016

生態(tài)系統(tǒng)Ecosystem土地覆被類(lèi)型Land cover type面積Area∥hm2比例Proportion∥%森林Forest落葉闊葉林139.90常綠針葉林1 047 334.017.8針闊混交林4 045.80.7落葉闊葉灌叢473 933.680.6稀疏林25.10稀疏灌叢4 860.50.8灌木園地162.90合計(jì)587 901.2100草地Grassland溫性草原386 154.027.2高寒草原351 879.924.8溫性草甸286 814.320.2高寒草甸73 735.55.2稀疏草地318 725.222.5合計(jì)1 417 308.9100農(nóng)田Farmland旱地32 700.4100濕地Wetland草本濕地6 578.071.3湖泊579.76.3水庫(kù)/坑塘1 257.013.6河流808.68.8合計(jì)9 223.3100人工Artificial建設(shè)用地1 943.466.9交通用地377.413.0采礦場(chǎng)583.020.1合計(jì)2 903.8100荒漠Desert裸巖562 803.993.6戈壁2 442.00.4裸土35 932.26.0鹽堿地7.00合計(jì)601 185.1100冰川積雪Glacier and snow29 437.4100

表2 祁連山國(guó)家級(jí)自然保護(hù)區(qū)土地覆被面積

表3 2000—2005年生態(tài)系統(tǒng)變化轉(zhuǎn)換矩陣

表4 2005—2010年生態(tài)系統(tǒng)變化轉(zhuǎn)換矩陣

表5 2010—2016年生態(tài)系統(tǒng)變化轉(zhuǎn)換矩陣

圖4 2000—2016年祁連山自然保護(hù)區(qū)EVI空間分布Fig.4 EVI spatial distribution of Qilian Mountain Nature Reserve from 2000 to 2016

圖5 2000—2016年祁連山自然保護(hù)區(qū)植被覆蓋度空間分布Fig.5 Spatial distribution of vegetation coverage in Qilian Mountain Nature Reserve from 2000 to 2016

3 結(jié)論

(1)祁連山獨(dú)特的自然條件和人為活動(dòng)的影響，決定了各土地覆被的構(gòu)成及面積相差較大。2016年整個(gè)保護(hù)區(qū)生態(tài)系統(tǒng)組成以草地為主，占總面積的52.9%；其次是荒漠，占22.4%，森林占21.9%；農(nóng)田占1.2%，冰川積雪占1.1%，濕地占0.3%；人工生態(tài)系統(tǒng)最小，僅占0.2%。

(2)森林主要以落葉闊葉灌叢為主，面積473 933.6 hm2，占森林面積的80.6%；落葉闊葉灌叢主要分布在實(shí)驗(yàn)區(qū)，面積占該類(lèi)林地總面積的44.4%；其次分布在核心區(qū)，占26.2%；分布在外圍保護(hù)地帶的面積最少，僅占8.5%。溫性草原最多，面積386 154.0 hm2，也僅占草地總面積的20.2%。

(3)近16年來(lái)祁連山國(guó)家級(jí)自然保護(hù)區(qū)主要生態(tài)系統(tǒng)整體處于穩(wěn)定狀態(tài)，2000—2005年自然保護(hù)區(qū)整體變化趨勢(shì)是森林和草地生態(tài)系統(tǒng)及冰川積雪面積減少，農(nóng)田、濕地、荒漠和人工生態(tài)系統(tǒng)面積增加。2005—2010年，森林、農(nóng)田、濕地、荒漠生態(tài)系統(tǒng)和冰川積雪面積凈減少，草地和人工生態(tài)系統(tǒng)面積凈增加。2010—2016年，草地、農(nóng)田和荒漠生態(tài)系統(tǒng)面積凈減少，森林、濕地、人工生態(tài)系統(tǒng)和冰川積雪面積是凈增加。

(4)祁連山自然保護(hù)區(qū)年EVI總體上呈現(xiàn)片狀分布特征，且存在一定的空間分布差異。EVI高值即植被覆蓋指數(shù)高的區(qū)域主要分布在自然保護(hù)區(qū)的東南區(qū)域，西北區(qū)域EVI值普遍較低。每年EVI植被指數(shù)分布大體一致。2000—2016年植被覆蓋度變化中，2004年植被覆蓋度最低，為25.98%。整體上植被覆蓋度有上升的態(tài)勢(shì)特征。年植被覆蓋度變化與生長(zhǎng)季植被覆蓋度變化基本一致。