劉櫟杉， 延軍平， 李雙雙

(陜西師范大學(xué) 旅游與環(huán)境學(xué)院， 陜西 西安 710062)

近年來，隨著“三北防護(hù)林”、“退耕還林還草”等防沙治沙工程實(shí)施的不斷深入，中國(guó)北方地區(qū)沙漠化較之20世紀(jì)末期有所減輕[1]，荒漠化整體擴(kuò)展趨勢(shì)已得到初步抑制[2-4]。黃森旺等[5]通過近25 a三北防護(hù)林工程區(qū)土地退化及驅(qū)動(dòng)力分析，發(fā)現(xiàn)三北防護(hù)林建設(shè)區(qū)土地退化整體呈減弱趨勢(shì)，大部分省區(qū)退化程度都有所減弱，其中尤以內(nèi)蒙古中東部、青海東北部、新疆北部最為明顯。

青藏高原作為亞洲乃至北半球氣候變化的“感應(yīng)器”和“敏感區(qū)”，對(duì)中國(guó)乃至亞洲生態(tài)安全具有重要的屏障作用[6-10]。青海作為青藏高原的一部分，地理環(huán)境復(fù)雜，生態(tài)環(huán)境脆弱，是氣候變化的敏感區(qū)，也是東部季風(fēng)區(qū)和青藏高原區(qū)的過渡帶，國(guó)家“退耕還林還草”重點(diǎn)建設(shè)區(qū)，長(zhǎng)江、黃河、瀾滄江3大河流的發(fā)源地。地表植被覆蓋動(dòng)態(tài)變化研究對(duì)該地區(qū)環(huán)境保護(hù)和生態(tài)建設(shè)具有重要的意義。

無論是全年還是生長(zhǎng)季，青藏高原植被變化與氣溫變化相關(guān)性明顯高于降水，高原植被對(duì)溫度變化具有較強(qiáng)的敏感性[11]。Lean等[12]提出用統(tǒng)計(jì)方法模擬全球溫度序列，共考慮4個(gè)因素，即人類活動(dòng)(包括溫室氣體及氣溶膠)、太陽輻射、火山氣溶膠及ENSO。王紹武等[13]在此基礎(chǔ)上進(jìn)一步分析，認(rèn)為1999—2008年有La Nina 事件，太陽輻射則由強(qiáng)到弱，但人類活動(dòng)仍然呈線性增長(zhǎng)。ENSO 及太陽輻射帶來的降冷，在相當(dāng)大程度上抵消了人類活動(dòng)造成的變暖。但是，這些自然因素并未改變氣候變化的趨勢(shì)。2001—2010年仍是全球有觀測(cè)記錄以來最暖的10 a，這一點(diǎn)在青藏高原體現(xiàn)的尤為明顯[14]。氣溫持續(xù)增高導(dǎo)致活動(dòng)積溫增加，一方面有利于青海東部濕潤(rùn)地區(qū)植被生長(zhǎng)，另一方面卻加劇了青海西北地區(qū)干旱化趨勢(shì)，影響植被覆蓋狀況改善。在全球變暖背景下，關(guān)注青海省生態(tài)恢復(fù)時(shí)空格局響應(yīng)特征，對(duì)區(qū)域生態(tài)建設(shè)具有重要的實(shí)踐意義。

基于此，本文選取青海省作為研究區(qū)域，利用2000—2009年遙感數(shù)據(jù)，以像元為基本計(jì)算單元，分析氣候變暖背景下“退耕還林還草”實(shí)施10 a來高原植被覆蓋的時(shí)空變化特征，以期為區(qū)域可持續(xù)發(fā)展、生態(tài)修復(fù)工程規(guī)劃以及生態(tài)環(huán)境保護(hù)提供一些理論依據(jù)。

1 研究區(qū)概況

青海省地域遼闊，主體位于青藏高原東北部，西北內(nèi)陸腹地，處于3大自然區(qū)的交匯地帶。地理坐標(biāo)為北緯31°39′—39°19′，東經(jīng)89°35′—103°04′，總面積7.21×105km2。由于地質(zhì)構(gòu)造控制，地貌基本格局呈現(xiàn)北西西—南東東走向，大地貌單元基本上沿緯向呈帶狀分布，自北向南依次是祁連山—阿爾金山、柴達(dá)木盆地—共和盆地—河湟谷地—黃南低地、東昆侖山脈、青南高原、唐古拉山脈。這些近乎東西走向的高大山脈，成為重要的自然地理分界線和行政區(qū)劃的界山。東部河湟谷地是黃土高原與青藏高原的過渡地帶，受季風(fēng)氣候影響，自然條件相對(duì)較好，也是我國(guó)第一階梯向第二階梯過渡地帶，蘊(yùn)藏著豐富的水力資源。西北部的柴達(dá)木盆地，干旱少雨，具有典型的荒漠景觀。南部為青南高原，是我國(guó)三江源生態(tài)重點(diǎn)保護(hù)區(qū)，具有高寒特征的景觀特色[15]。在全球變化和人類活動(dòng)的綜合影響下，青藏高原呈現(xiàn)出生態(tài)系統(tǒng)穩(wěn)定性降低、資源環(huán)境壓力增大等問題，突出表現(xiàn)為冰川退縮顯著，土地退化形勢(shì)嚴(yán)峻，水土流失加劇，生物多樣性威脅加大與珍稀生物資源減少，自然災(zāi)害增多等[10]。

2 資料來源與研究方法

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

本研究中遙感數(shù)據(jù)來源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)，其中NDVI數(shù)據(jù)為2000—2009年16 d最大值合成Terra—MODIS L3級(jí)產(chǎn)品數(shù)據(jù)，分辨率為1 km。選取8月作為研究時(shí)段，此時(shí)段區(qū)域植被覆蓋達(dá)到一年中最高、云量較少且具有一定持續(xù)性，可以較好反映區(qū)域植被時(shí)空變化特征。由于青海省區(qū)域面積較大，Terra—MODIS L3數(shù)據(jù)需要兩幅影像進(jìn)行區(qū)域全覆蓋(h25v05，h26v05)，使用GIS 9.3軟件將兩幅影像拼接和重投影，并利用青海省各州行政邊界，提取青海省區(qū)域及各州年植被覆蓋數(shù)據(jù)。NDVI數(shù)字高程數(shù)據(jù)(DEM)為分辨率90 m SRTM。

2.2 研究方法

2.2.1 趨勢(shì)分析 一元線性回歸分析可以模擬每個(gè)柵格的變化趨勢(shì)，以單個(gè)像元時(shí)間變化特征反映整個(gè)空間變化規(guī)律，綜合反映區(qū)域時(shí)空格局演變[16-17]。其中Slope為像元NDVI回歸方程的斜率，若Slope值為正，表示隨時(shí)間變化植被指數(shù)升高，區(qū)域植被覆蓋度存在增加趨勢(shì)，且數(shù)值愈大植被覆蓋度增加趨勢(shì)愈明顯；反之，Slope值為負(fù)，表示隨時(shí)間變化植被指數(shù)呈下降趨勢(shì)。其計(jì)算公式如下：

(1)

式中：S——像元NDVI回歸方程的斜率；n——年跨度；i——年序號(hào)； NDVIi——第i年NDVI值。下同。

為了更好地評(píng)價(jià)青海省植被恢復(fù)狀況，參照已有研究[17]，將Slope分為7個(gè)等級(jí)：嚴(yán)重退化、輕微退化、中度退化、基本不變、輕微改善、中度改善、明顯改善(表1)。

表1 植被覆蓋變化趨勢(shì)分類標(biāo)準(zhǔn)

2.2.2 標(biāo)準(zhǔn)差分析 標(biāo)準(zhǔn)差表示個(gè)數(shù)據(jù)變量偏離常態(tài)的距離的平均數(shù)，是離差平方和和平均后的方根，能反映一個(gè)數(shù)據(jù)集的離散程度，其值越大，說明該地區(qū)在研究時(shí)段內(nèi)各像元NDVI距離平均值越遠(yuǎn)，即該段時(shí)間植被覆蓋的年際變化較大，植被穩(wěn)定性差[18]。結(jié)合區(qū)域自然地理特征，如柴達(dá)木盆地植被稀少，其大幅波動(dòng)變化可能性相對(duì)較小，可以進(jìn)一步驗(yàn)證穩(wěn)定性和趨勢(shì)分析結(jié)論的可靠性：

(2)

3 結(jié)果與分析

3.1 青海省植被覆蓋年際變化特征

2000—2009年青海省植被覆蓋呈現(xiàn)明顯的增加趨勢(shì)，增速為0.018/10 a，快于三北防護(hù)林工程區(qū)1982—2006年近25 a植被覆蓋平均增速0.007/10 a[19]，但慢于陜甘寧黃土高原區(qū)2000—2009年植被覆蓋平均增速0.032/10 a[20]。2000年植被NDVI為0.317，2009年上升到0.340，增長(zhǎng)率為7.2%。其中，海東地區(qū)植被覆蓋改善最為明顯，10 a間增加了18.1%，其次是海南州16.5%，格爾木11.4%，海西州10.5%。海北州和果洛州植被條件相對(duì)較好，其增長(zhǎng)率與其他地區(qū)相比顯得比較緩慢，10 a分別增長(zhǎng)了0.9%和1.3%(表2)。青海植被恢復(fù)大致可以分為兩個(gè)階段，第一階段為2000—2005年植被快速增長(zhǎng)，即“退耕還林草”工程實(shí)施前6 a，青海省植被NDVI由2000年0.317上升到2005年0.349，達(dá)到10 a最高；第二階段為2005—2009年，青海省植被呈現(xiàn)波動(dòng)下降，其中2005—2006年和2007—2008年NDVI下降速率分別為4.92%，3.18%(圖1)。

表2 青海省植被覆蓋變化統(tǒng)計(jì)特征

圖1 2000-2009年青海省NDVI變化曲線

3.2 青海省植被覆蓋空間格局演變

3.2.1 青海植被格局基本特征 青海省植被格局基本特征表現(xiàn)在以下兩個(gè)方面。

(1)在空間上，受水熱、地形等自然地理?xiàng)l件限制，青海省“南高北低，東高西低”的植被格局基本穩(wěn)定(圖2)。柴達(dá)木盆地、青南高原西北形成兩大植被相對(duì)低區(qū)。祁連山中東段—河湟谷地，黃南低地—果洛州黃河源區(qū)，格爾木南部—玉樹州形成3大植被相對(duì)高區(qū)。

(2)青海省植被覆蓋最高的區(qū)域?yàn)辄S南州，2009年植被NDVI為0.702，其次是西寧地區(qū)為0.635，植被覆蓋相對(duì)較低的區(qū)域?yàn)楹Ｎ髦?、格爾木，植被NDVI分別為0.151，0.262，其中海西州植被覆蓋僅為黃南州的1/4。

圖2 2000-2009年青海省植被覆蓋空間格局變化

3.2.2 植被覆蓋空間格局演變態(tài)勢(shì) 實(shí)施“退耕還林還草”10 a間，青海省植被覆蓋整體呈增加趨勢(shì)(表3)，生態(tài)退化初步得到遏制，但其空間分布存在差異(附圖6)。空間格局演變特征主要表現(xiàn)為：(1)呈退化趨勢(shì)區(qū)域比重為18.40%，集中分布于祁連山中東段和青南高原三江源地區(qū)；(2)基本不變區(qū)域比重為29.15%，集中分布于柴達(dá)木盆地中部和西北部風(fēng)蝕殘丘、流沙區(qū)，東昆侖山南部可可西里高原亦有零星分布；(3)輕微改善區(qū)比重最大(32.66%)，除柴達(dá)木盆地中部和西北部均有分布，相對(duì)比較集中的區(qū)域有柴達(dá)木盆地邊緣和東昆侖山南側(cè)；(4)中度和明顯改善區(qū)比重為19.80%，連片分布于柴達(dá)木盆地東南邊緣，青海湖盆地、茶卡—共和盆地、河湟谷地、黃南低地，在南部青南高原、巴顏喀拉山和阿尼瑪卿山西側(cè)亦有帶狀分布；(5)明顯改善區(qū)域相對(duì)較少，比重為6.48%，空間分布主要集中在兩個(gè)地區(qū)：柴達(dá)木盆地東南邊緣形成一條南北向“W”型綠色屏障(98°經(jīng)線周邊)；青海東部低山、丘陵、盆地區(qū)。

表3 2000-2009年青海省植被覆蓋演化趨勢(shì)

3.2.3 植被覆蓋空間格局穩(wěn)定性 青海省2000—2009年植被覆蓋標(biāo)準(zhǔn)差介于0.0～0.05，穩(wěn)定性存在明顯地域差異(圖3)?？臻g格局主要表現(xiàn)為“東南波動(dòng)，西北穩(wěn)定，高原溫帶波動(dòng)，高原寒帶、亞寒帶穩(wěn)定”。其中，柴達(dá)木盆地、可可西里高原植被變化相對(duì)穩(wěn)定，冷龍嶺、大通河谷地、阿尼瑪卿山西側(cè)、青南高原東南部變化幅度居中，青海東部中低山地、丘陵、盆地地區(qū)變化幅度最為明顯。

青海省以昆侖山為界，北部大部分屬于高原溫帶，降水量由東向西遞減，干旱程度逐漸增強(qiáng)，植被分布依次為半干旱草原、干旱草原、半荒漠草原、荒漠草原；南部屬于高原亞寒帶，降水量自東南向西北遞減，植被由東南向西北依次為寒溫針葉林、高山灌叢草甸、高寒草甸、高寒草原、高寒荒漠化草原。柴達(dá)木盆地、青南高原西北植被相對(duì)較低，加之氣候寒冷，植被以半荒漠草原、荒漠草原為主，人類活動(dòng)影響較少，所以近10 a來植被覆蓋呈現(xiàn)出低位穩(wěn)定變化的特征；青海東部中低山地、丘陵、盆地地區(qū)，人為因素驅(qū)動(dòng)明顯，這里為我國(guó)“退耕還林草”工程實(shí)施的重點(diǎn)區(qū)域，氣候相對(duì)適宜，降水量相對(duì)較多，植被易恢復(fù)，因此，區(qū)域植被覆蓋表現(xiàn)出波動(dòng)上升的特征。

圖3 2000-2009年青海省植被覆蓋變化穩(wěn)定程度

4 結(jié) 論

(1)2000—2009年青海省植被覆蓋呈現(xiàn)明顯的增加趨勢(shì)，增速為0.018/10 a，快于三北防護(hù)林工程區(qū)1982—2006年近25 a植被覆蓋平均增速0.007/10 a，但慢于陜甘寧黃土高原區(qū)2000—2009年植被覆蓋平均增速0.032/10 a。

2000—2009年青海省植被恢復(fù)具有階段性，大致可以分為兩個(gè)階段：第一次階段為2000—2005年植被快速增長(zhǎng)；第二階段為2005—2009年植被呈現(xiàn)波動(dòng)下降。

2000—2009年青海省海東地區(qū)植被覆蓋改善最為明顯，10 a間增加了18.1%，其次是海南州16.5%，格爾木11.4%，海西州10.5%。祁連山東段的海北州和黃河源的果洛州植被條件相對(duì)較好，其增長(zhǎng)率相對(duì)比較緩慢，10 a間分別增長(zhǎng)了0.9%和1.3%。

(2)受水熱、地形等自然地理?xiàng)l件限制，青海省“南高北低、東高西低”的植被格局基本穩(wěn)定。柴達(dá)木盆地、青南高原西北形成兩大植被相對(duì)低區(qū)，“祁連山中東段—河湟谷地”，“黃南低地—果洛州黃河源區(qū)”，“格爾木南部—玉樹州”形成3大植被相對(duì)高區(qū)。

(3)2000—2009年青海省植被覆蓋呈退化趨勢(shì)區(qū)域比重為18.40%，集中分布于祁連山中東段和青南高原三江源地區(qū)；明顯改善區(qū)相對(duì)較少(6.48%)，空間上集中分布在兩個(gè)地區(qū)：一是柴達(dá)木盆地東南邊緣形成一條南北向“W”型綠色屏障(98°經(jīng)線)，另一是青海東部低山、丘陵、盆地區(qū)。

(4)青海省2000—2009年植被覆蓋標(biāo)準(zhǔn)差介于0.0～0.05，穩(wěn)定性存在明顯地域差異?？臻g格局主要表現(xiàn)為“東南波動(dòng)，西北穩(wěn)定，高原溫帶波動(dòng)，高原寒帶、亞寒帶穩(wěn)定”。

氣候變暖背景下，青海省生態(tài)恢復(fù)的可持續(xù)性依然是未來研究的熱點(diǎn)話題，還有許多工作需要進(jìn)一步探索。如進(jìn)一步評(píng)估重點(diǎn)和典型區(qū)域(青南高原三江源地區(qū)、祁連山東段山區(qū)、湟水谷地、黃南低地等)生態(tài)系統(tǒng)恢復(fù)適應(yīng)性和穩(wěn)定性；進(jìn)一步量化氣候變化和人類活動(dòng)對(duì)青海省植被恢復(fù)影響，系統(tǒng)評(píng)價(jià)兩者在生態(tài)恢復(fù)過程中的相對(duì)作用；進(jìn)一步探討青海省生態(tài)系統(tǒng)良性恢復(fù)后系統(tǒng)服務(wù)功能，關(guān)注生態(tài)建設(shè)向生態(tài)發(fā)展的轉(zhuǎn)型，促進(jìn)區(qū)域經(jīng)濟(jì)和生態(tài)良性互動(dòng)發(fā)展。

[參考文獻(xiàn)]

[1]國(guó)家林業(yè)局.第五次中國(guó)荒漠化和沙漠化公報(bào)[Z].2011.

[2]王濤,吳薇,薛嫻,等.近50年來中國(guó)北方沙漠化土地的時(shí)空變化[J].地理學(xué)報(bào),2004,59(2):203-212.

[3]王濤,陳廣庭,趙哈林,等.中國(guó)北方沙漠化過程及其防治研究的新進(jìn)展[J].中國(guó)沙漠,2006,26(4):507-516.

[4]許端陽,李春蕾,莊大方,等.氣候變化和人類活動(dòng)在沙漠化過程中相對(duì)作用評(píng)價(jià)綜述[J].地理學(xué)報(bào),2011,66(1):68-76.

[5]黃森旺,李曉松,吳炳方,等.近25年三北防護(hù)林工程區(qū)土地退化及驅(qū)動(dòng)力分析[J].地理學(xué)報(bào),2012,67(5):589-598.

[6]李輝霞,劉國(guó)華,傅伯杰.基于NDVI的三江源地區(qū)植被生長(zhǎng)對(duì)氣候變化和人類活動(dòng)的響應(yīng)研究[J].生態(tài)學(xué)報(bào),2011,31(19):5495-5504.

[7]劉紀(jì)遠(yuǎn),徐新良,邵全琴.近30年來青海三江源地區(qū)草地退化的時(shí)空特征[J].地理學(xué)報(bào),2008,63(4):364-376.

[8]張鐿鋰,劉林山,擺萬奇,等.黃河源地區(qū)草地退化空間特征[J].地理學(xué)報(bào),2006,61(1):4-12.

[9]董立新,王文科,孔金玲,等.黃河上游瑪多縣生態(tài)環(huán)境變化遙感監(jiān)測(cè)及成因分析[J].水土保持通報(bào),2005,25(4):68-72.

[10]孫鴻烈,鄭度,姚檀棟,等.青藏高原國(guó)家生態(tài)安全屏障保護(hù)與建設(shè)[J].地理學(xué)報(bào),2012,67(1):3-12.

[11]張戈麗,歐陽華,張憲洲,等.基于生態(tài)地理分區(qū)的青藏高原植被覆被變化及其對(duì)氣候變化的響應(yīng)[J].地理研究,2010,29(11):2004-2016.

[12]Lean J L, Rind D H. How natural and anthropogenic influences alter global and regional surface temperatures: 1889 to 2006 [J]. Geophys. Res. Lett., 2008, 35, L18701, doi: 10.1029/2008 GL 034864.

[13]王紹武,羅勇,唐國(guó)利,等.近10年全球變暖停滯了嗎?[J].氣候變化研究進(jìn)展,2010,6(2):95-99.

[14]唐國(guó)利,羅勇,黃建斌,等.氣候變暖在繼續(xù)[J].氣候變化研究進(jìn)展,2012,8(4):235-242.

[15]張忠孝.青海地理[M].西寧：青海人民出版社,2004.

[16]馬明國(guó),王建,王雪梅.基于遙感的植被年際變化及其與氣候關(guān)系研究進(jìn)展[J].遙感學(xué)報(bào),2006,10(3):421-431.

[17]宋怡,馬明國(guó).基于SPOT VEGETATION數(shù)據(jù)的中國(guó)西北植被覆蓋變化分析[J].中國(guó)沙漠,2007,27(1):89-93.

[18]徐建華.現(xiàn)代地理學(xué)中的數(shù)學(xué)方法[M].北京:高等教育出版社,2002.

[19]王強(qiáng),張勃,戴聲佩,等.基于GIMMS AVHRR NDVI數(shù)據(jù)的三北防護(hù)林工程區(qū)植被覆蓋動(dòng)態(tài)變化[J].資源科學(xué), 2011,33(8):1613-1620.

[20]李雙雙,延軍平,萬佳.近10年陜甘寧黃土高原區(qū)植被覆蓋時(shí)空變化特征[J].地理學(xué)報(bào),2012,67(7):960-970.