近30年來青藏高原高寒草地NDVI動態(tài)變化對自然及人為因子的響應(yīng)

張 江,袁旻舒,張 婧,李函微,王潔儀,張 賢,鞠佩君,蔣海波,陳 槐,朱求安,3,4,*

1 西北農(nóng)林科技大學(xué)林學(xué)院,生態(tài)預(yù)測與全球變化研究中心,楊凌 712100 2 中國科學(xué)院成都生物研究所山地生態(tài)恢復(fù)與生物資源利用重點實驗室,成都 610041 3 河海大學(xué)水文水資源學(xué)院,南京 210098 4 國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心,國家科技基礎(chǔ)條件平臺,北京 100101

氣候變化和人類活動對陸地生態(tài)系統(tǒng)的影響一直是全球氣候變化研究中的熱點[1-2]。草地生態(tài)系統(tǒng)約占陸地總面積的1/5,是陸地生態(tài)系統(tǒng)重要組成部分,在調(diào)節(jié)氣候、水土保持、防風固沙等方面發(fā)揮著重要作用[3-4]。有研究表明,高海拔、高緯度地區(qū)草地植被動態(tài)對氣候變化及人類活動的反應(yīng)較為敏感[5-6],主要表現(xiàn)為植被生長狀態(tài)及分布范圍的變化[7]。植被的分布由于受到氣候模式的制約,在氣候暖化過程中植被可能會向兩極高緯度地區(qū)及高海拔地區(qū)移動,進而改變物種的分布范圍[8]。與此同時,人口遷移、畜牧業(yè)發(fā)展及土地利用變化等人類活動的增強也會影響植被的生長狀態(tài)及分布范圍[9]。

青藏高原是全球海拔最高的獨特地域單元[10],平均海拔超過4000 m,素有“世界第三極”之稱,屬于氣候變化敏感區(qū)和生態(tài)脆弱帶[11]。草地是青藏高原主要的植被類型,總面積約為1.59×106km2。因此,青藏高原被認為是研究草地植被動態(tài)對氣候變化及人類活動響應(yīng)機制的理想場所[12]。一些研究表明,在過去的三十多年,隨著氣候變化和人類活動的加劇,青藏高原草地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能發(fā)生了巨大變化,草地呈現(xiàn)出不同程度的退化且日趨惡化,已嚴重阻礙了當?shù)厣鐣?jīng)濟的可持續(xù)發(fā)展[13-14]。也有一些研究表明,青藏高原草地退化趨勢已有所緩解,呈現(xiàn)出“整體改善、局部退化”的趨勢[7, 15-18]。在青藏高原草地植被變化的研究中,氣溫和降水被認為是導(dǎo)致草地植被變化的主要氣候因子,而放牧及人口數(shù)量變化引起的其他人類活動(如土地利用變化、植樹造林等)則被認為是導(dǎo)致草地植被變化的主要人為因子[19-20]。

歸一化植被指數(shù)(Normalized difference vegetation index, NDVI)因能有效地反映植被覆蓋程度、生長狀況、生物量及凈初級生產(chǎn)力等植被參數(shù),被廣泛應(yīng)用于大尺度植被動態(tài)研究中[21-22]。例如卓嘎等[23]利用NDVI數(shù)據(jù)分析了西藏地區(qū)植被的分布狀況及變化趨勢；董璐等[24]利用NDVI數(shù)據(jù)研究了新疆地區(qū)植被生長對氣溫的響應(yīng)；李鵬[25]利用NDVI數(shù)據(jù)研究了青藏高原植被物候?qū)O端氣候的響應(yīng)。目前在區(qū)域尺度上,關(guān)于氣候變化對青藏高原植被動態(tài)影響的研究較多,例如楊元合[12]、Gao等[26]、劉軍會等[16]利用NDVI數(shù)據(jù)研究分析了青藏高原草地植被動態(tài)與氣候因子氣溫、降水之間的關(guān)系。而關(guān)于人類活動對青藏高原草地植被動態(tài)影響的研究多為樣點實驗,例如周華坤等[27]、宋磊等[28]在青藏高原高寒草原上進行放牧試驗來研究放牧強度對植被群落的高度、蓋度、生物量及多樣性的影響,在區(qū)域尺度上的研究相對較少。因此在區(qū)域尺度上綜合分析氣候變化、人類活動對青藏高原高寒草地植被變化影響的研究亟待開展。

本研究基于32年(1982—2013年)青藏高原143個縣區(qū)生長季草地的NDVI、平均氣溫、降水、夏季牧場放牧強度及人口數(shù)量等數(shù)據(jù),采用趨勢分析法和面板數(shù)據(jù)模型(Panel data model)在區(qū)域尺度上研究了青藏高原高寒草地NDVI動態(tài)變化及其對氣候變化和人類活動的響應(yīng),以期為青藏高原地區(qū)生態(tài)環(huán)境保護和區(qū)域經(jīng)濟發(fā)展提供理論參考依據(jù)。

1 研究方法

1.1 研究區(qū)域

青藏高原位于中國西南部邊緣(75—105°E, 27.5—37.5°N),總面積約為2.57×106km2,境內(nèi)有河流、冰川、湖泊等多種地貌,地形復(fù)雜多變,亦有“世界屋脊”之稱,區(qū)域氣候易受東亞季風、印度季風、西南季風及西風環(huán)流系統(tǒng)的影響,水熱分布差異明顯[10-11]。本研究基于青藏高原境內(nèi)的143個縣級行政單元,主要包括青海省、西藏自治區(qū)及四川省部分地區(qū)(圖1)。

圖1 研究區(qū)域

1.2 數(shù)據(jù)來源與預(yù)處理

1.2.1青藏高原草地空間分布數(shù)據(jù)

植被空間分布數(shù)據(jù)來源于歐洲航空局(European Space Agency, ESA(http://www.esa.int/ESA)),時間跨度為1992—2013年,空間分辨率為300 m。將草地植被類型從數(shù)據(jù)集中分離得到1992—2013年青藏高原草地分布數(shù)據(jù)。本研究的時間范圍為1982—2013年,由于缺少1982—1991年植被空間分布數(shù)據(jù),因此1982—1991年的草地分布數(shù)據(jù)采用1992年青藏高原的草地分布數(shù)據(jù)補充完整。將獲得的草地分布數(shù)據(jù)與各縣區(qū)疊加,得到各縣區(qū)草地分布,進而統(tǒng)計得到1982—2013年各縣區(qū)的草地面積變化序列。

1.2.2NDVI數(shù)據(jù)

本研究選用美國航空航天局(National Aeronautics and Space Administration, NASA)全球監(jiān)測與模型研究組(Global Inventory Monitoring and Modelling Studies, GIMMS)提供的1982年1月至2013年12月,每十五天合成的空間分辨率為8 km的NDVI3g(GIMMS NDVI Version3(https://ecocast.arc.nasa.gov/data/pub/gimms/))數(shù)據(jù)集來描述青藏高原高寒草地的生長狀況及覆蓋程度。該數(shù)據(jù)集消除了火山爆發(fā)、太陽高度角、云覆蓋等因素的影響,已被廣泛應(yīng)用于植被動態(tài)的研究中[29-30]。每月的NDVI采用了最大值合成法(MVC),該方法可以進一步消除云霧、大氣等因素的干擾,提高數(shù)據(jù)的精度[31]。青藏高原冬季時期草地大部分被積雪覆蓋,因此本研究使用生長季平均NDVI來反映草地生長狀態(tài),以消除積雪覆蓋對NDVI的影響,生長季為每年5月至9月[12]。將獲得的青藏高原生長季平均NDVI數(shù)據(jù)與青藏高原各縣區(qū)草地分布圖進行疊加,得到各縣區(qū)草地NDVI的變化序列。

1.2.3氣候數(shù)據(jù)

為分析氣候變化對青藏高原高寒草地NDVI變化的影響,研究選取了對植被生長影響最為密切的生長季平均溫度和生長季降水作為影響草地NDVI變化主要的氣候因子。基于全國標準氣象站的數(shù)據(jù)(中國氣象局(http://www.cma.gov.cn/)),利用薄板樣條插值法[32]對氣溫和降水進行插值,生成了全國氣溫和降水數(shù)據(jù),分辨率為10 km。該數(shù)據(jù)集在其他研究中得到了驗證與應(yīng)用[33-34]。從上述全國氣象數(shù)據(jù)中提取青藏高原生長季平均氣溫和生長季降水數(shù)據(jù)，并將其分別與青藏高原各縣區(qū)草地分布數(shù)據(jù)進行疊加,統(tǒng)計得到1982—2013年各縣區(qū)草地分布區(qū)的生長季平均氣溫與降水量的變化序列。

1.2.4人類活動數(shù)據(jù)

為了分析人類活動對青藏高原高寒草地NDVI變化的影響,研究選取了影響高寒草地植被變化人類活動因子中最重要的兩個指標,即反映畜牧業(yè)發(fā)展規(guī)模、直接影響草地生長狀態(tài)的放牧強度以及代表城鎮(zhèn)規(guī)模、直接影響放牧強度變化的人口數(shù)量。1982—2013年各縣區(qū)畜牧種類及數(shù)量和人口數(shù)據(jù)來源于各省、市、縣(志)統(tǒng)計年鑒(國家數(shù)字電子圖書館(http://www.nlc.cn/))。

在計算放牧強度時,為消除不同畜牧種類及數(shù)量之間的差異，均以綿羊為基準單位進行羊單位轉(zhuǎn)換,轉(zhuǎn)換比值參考《草業(yè)科學(xué)概論》[35]。為了與生長季草地NDVI、氣溫、降水數(shù)據(jù)相對應(yīng),研究計算了各縣區(qū)夏季牧場的放牧強度,其計算公式如下：

(1)

式中,GIi、Ui和Ai分別為第i年該縣區(qū)的放牧強度(羊單位/hm2)、羊單位(只/個)和草地面積(hm2)；R和S分別為可利用草地面積占比和夏季牧場占比。其中各縣區(qū)可利用草地面積、夏季牧場占比數(shù)據(jù)均來源于各省、市、縣(志)統(tǒng)計年鑒。

1.3 分析方法

1.3.1放牧強度分區(qū)

為研究不同放牧強度對草地NDVI的影響,本研究將各縣區(qū)的放牧強度擬劃分為輕、中、重度三個放牧等級區(qū)。由于青藏高原各縣區(qū)草地載畜能力不同,放牧強度的閾值劃分(即單位面積上牲畜的數(shù)量)也存在差異。通過文獻收集,對各研究點放牧強度(輕、中、重度)的定義閾值進行了統(tǒng)計(表1),研究點的分布如圖2。將計算的各縣放牧強度的均值與表1中放牧強度的閾值進行對比,對各縣區(qū)進行等級劃分。若同一個縣區(qū)存在多個放牧實驗,放牧強度的閾值取其均值,區(qū)域內(nèi)缺少放牧實驗時,則參考鄰近縣區(qū)。

表1 放牧實驗中放牧強度閾值和數(shù)據(jù)來源

圖2 青藏高原輕、中、重度放牧分區(qū)示意圖

1.3.2趨勢分析

草地NDVI、氣候因子及人為因子在青藏高原及各縣區(qū)的變化趨勢采用趨勢分析法。計算公式如下：

(2)

式中,k為斜率；Yi為第i年各變量值；i為年變量；n為年數(shù)。

在空間尺度中,根據(jù)斜率k及P值大小,將各變量變化趨勢劃分為-3、-2、-1、+1、+2和+3六個等級,分別代表極顯著下降(k<0,P≤0.01)、顯著下降(k<0, 0.010.05)、上升不顯著(k>0,P>0.05)、顯著上升(k>0, 0.010,P≤0.01)。

1.3.3面板數(shù)據(jù)模型

面板數(shù)據(jù)Panel Data即空間和時間同時擴展的混合數(shù)據(jù)集合,是指由變量Y關(guān)于N個不同對象在T個觀測時期所得到二維結(jié)構(gòu)數(shù)據(jù)。其優(yōu)勢在于,面板數(shù)據(jù)能同時反映某一時期各個個體數(shù)據(jù)的規(guī)律和每個個體隨時間變化的規(guī)律。與截面數(shù)據(jù)模型相比,面板數(shù)據(jù)模型可以控制不可觀測變量所導(dǎo)致的最小二乘法(OLS)估計的偏差,使得模型設(shè)定更加合理、模型參數(shù)的樣本估計量更加準確；與時間序列模型相比,該模型擴大了樣本信息、降低了變量間的共線性,提高了估計量的有效性[72]。面板分析選取了直接影響草地NDVI變化的氣候因子氣溫、降水及人類活動因子放牧強度。人口數(shù)量作為間接影響草地NDVI的因子也被引入模型(主要通過影響放牧數(shù)量、退耕還草等方式間接影響草地NDVI)。面板數(shù)據(jù)模型設(shè)定如下：

NDVIit=αi+βi1Tempit+βi2Precit+βi3GIit+βi4Popuit+μit

(3)

式中,NDVIit、Tempit、Precit、GIit和Popuit分別表示第i個縣區(qū)在第t年生長季的草地NDVI、生長季平均氣溫(℃)、生長季降水(mm)、夏季牧場的放牧強度(羊單位/hm2)及人口(百萬)；αi、βi1、βi2、βi3和βi4分別表示第i個縣區(qū)在第t年的截距項及4個自變量的系數(shù)；i為縣區(qū)個數(shù)；t為不同年份；μit為隨機擾動項。

面板數(shù)據(jù)模型分析中,采用z-score標準化法對各變量進行標準化處理,求得各變量的相對貢獻;采用未標準化的數(shù)據(jù)分析各變量單位變化量對NDVI的影響。為避免面板模型出現(xiàn)偽回歸現(xiàn)象,采用ADF單位根檢驗法[73]對面板數(shù)據(jù)進行單位根檢驗。面板數(shù)據(jù)有混合效應(yīng)、固定效應(yīng)和隨機效應(yīng)3種模型[74],模型的確立采用F檢驗和Hausman檢驗[75-76]。數(shù)據(jù)分析表明，ADF單位根檢驗結(jié)果拒絕其零假設(shè)(P<0.01)表明數(shù)據(jù)在時間序列上整體呈現(xiàn)平穩(wěn)性不存在單位根，F(xiàn)檢驗和Hausman檢驗結(jié)果均拒絕零假設(shè)(P<0.01),因此本研究面板分析在時間和截面上設(shè)定為固定效應(yīng)模型。

2 結(jié)果

2.1 青藏高原生長季草地NDVI、氣候因子及人為因子的時空變化

趨勢分析結(jié)果顯示,1982—2013年青藏高原生長季草地NDVI的變化范圍在0.431—0.471之間,草地NDVI整體上呈增長趨勢,但未達顯著水平(P>0.05)。青藏高原大部分縣區(qū)生長季草地NDVI呈增長趨勢,少部分縣區(qū)呈下降趨勢,各縣區(qū)的變化趨勢無明顯空間異質(zhì)性(圖3)。

32年來,青藏高原生長季平均氣溫呈顯著上升趨勢,增長率約為每10年0.420℃。各縣區(qū)生長季平均氣溫增長趨勢顯著,無明顯的空間異質(zhì)性。青藏高原生長季降水量整體上呈增長趨勢,最大、最小降水量分別為481.179 mm和381.635 mm。各縣區(qū)生長季降水量的變化趨勢呈明顯的空間異質(zhì)性,青藏高原西部和北部降水量顯著增加,而東部和南部降水量增長趨勢不明顯且部分縣區(qū)呈下降趨勢(圖3)。

1982—2013年青藏高原放牧數(shù)量逐年上升,夏季牧場放牧強度增長趨勢顯著,每年每公頃增長約0.014個羊單位。各縣區(qū)夏季牧場放牧強度的變化趨勢存在一定的空間異質(zhì)性,其中青藏高原中部和北部地區(qū)放牧強度呈下降趨勢,而東南部地區(qū)放牧強度呈顯著增長趨勢。32年來青藏高原人口增長趨勢顯著,總?cè)丝谟?982年的817萬增長至2013年的1200萬,每年大約增長13萬人。除部分縣區(qū)人口呈現(xiàn)下降趨勢外,其他各縣區(qū)人口增長趨勢顯著,且無明顯空間異質(zhì)性(圖3)。

圖3 1982—2013年青藏高原生長季草地NDVI、氣溫、降水、放牧強度及人口的年際變化與趨勢變化圖

2.2 各變量因子對草地NDVI的影響

分別對不同的放牧強度分區(qū)及青藏高原進行了面板模型分析,4個分析情形下模型結(jié)果均具有較高的擬合優(yōu)度,青藏高原全區(qū)以及輕度、中度和重度放牧區(qū)的模型決定系數(shù)分別為0.56、0.77、0.48和0.41(表2)。在輕度放牧區(qū),生長季草地NDVI與氣溫、降水及放牧強度之間均呈極顯著正相關(guān),與人口之間相關(guān)性不顯著。氣溫、降水及放牧強度每增加1單位,草地NDVI相應(yīng)增加0.0047、0.0007和0.0421,相對貢獻系數(shù)分別為0.0576、0.6751和0.2927。在中度放牧區(qū),生長季草地NDVI與氣溫、降水及放牧強度之間呈極顯著正相關(guān),與人口之間相關(guān)性仍不顯著。放牧強度對草地NDVI變化的貢獻系數(shù)減小至0.0073,而氣溫的貢獻系數(shù)上升至0.0234。在重度放牧區(qū),生長季草地NDVI與氣溫、降水呈極顯著正相關(guān),與放牧強度之間則呈極顯著負相關(guān),相對貢獻系數(shù)分別為0.2635、0.6889和-0.1525,同時放牧強度每增加1單位,草地NDVI相應(yīng)減少0.0016。全區(qū)分析結(jié)果表明,生長季草地NDVI與氣溫、降水之間呈極顯著正相關(guān),與放牧強度呈負相關(guān),但未達到顯著水平(P>0.05)。放牧強度每增加1個單位,草地NDVI減少0.0001。降水的相對貢獻系數(shù)最大,溫度次之,分別為0.6963和0.2810。

表2 面板數(shù)據(jù)模型估計結(jié)果

1982—2013年,生長季平均氣溫在不同放牧區(qū)中均呈顯著增長趨勢(P<0.01),生長季降水量在輕度放牧區(qū)呈顯著增長趨勢(P<0.01),但在中度、重度及全區(qū)中,生長季降水增長趨勢未達到顯著水平(P>0.05)(表3)。32年來,輕度放牧區(qū)生長季氣溫與降水量分別增加了約1.632℃和43.488 mm,兩者的增加使草地NDVI分別增加了0.0077和0.0304;放牧強度減少了0.096,導(dǎo)致草地NDVI變化了0.0054。在中度放牧區(qū)中,生長季平均氣溫、降水及放牧強度的增加分別導(dǎo)致草地NDVI增加了0.0300、0.0271和0.0014。在重度放牧區(qū)中,溫度與降水的增加使草地NDVI分別增加了0.0222和0.0238,但放牧強度的增加導(dǎo)致草地NDVI減少了約0.0015。對于整個青藏高原地區(qū),1982—2013年青藏高原草地生長季平均氣溫、降水量及放牧強度的增長率分別為每年0.042℃、1.056 mm和0.014個羊單位/hm2(圖3和表3),溫度和降水的增加使草地NDVI分別增加了約0.0317和0.0270,放牧強度的增加則導(dǎo)致草地NDVI減少了約0.00005。

表3 1982—2013年氣溫、降水、放牧強度及草地NDVI的變化量

3 分析與討論

3.1 青藏高原氣候及人類活動的變化

青藏高原作為氣候變化的敏感區(qū)對氣候變暖已表現(xiàn)出諸多響應(yīng)。近幾十年來,青藏高原氣溫顯著上升,降水量略微增加,總體呈現(xiàn)出“暖濕化”趨勢[77-79]。本研究結(jié)果與其類似,青藏高原生長季平均氣溫、生長季降水量整體上均呈增長趨勢。東部及南部地區(qū)氣溫增長趨勢顯著,降水增長趨勢不明顯,區(qū)域氣候呈“暖干化”趨勢；青藏高原西部和北部地區(qū)氣溫、降水增長趨勢顯著,區(qū)域氣候呈“暖濕化”趨勢(圖3)。大尺度環(huán)流系統(tǒng)會影響水汽的分布[80],在過去的幾十年青藏高原地區(qū)西風、東亞季風、印度季風的增強可能改變了青藏高原地區(qū)的水汽分布[81-82],導(dǎo)致了青藏高原東南部和西北部區(qū)域氣候差異的出現(xiàn)。

青藏高原是我國重要的草地畜牧業(yè)基地[83],放牧是其重要的生產(chǎn)方式之一[7]。近幾十年來,青藏高原人口顯著增加[84]。大規(guī)模人口對物質(zhì)資源(奶、肉、毛皮制品)的需求量不斷加大使得放牧數(shù)量也隨之增加,進而導(dǎo)致放牧強度逐年上升[85]。但青藏高原地域廣闊各縣區(qū)人口分布、經(jīng)濟發(fā)展水平、區(qū)域政策及放牧管理方式的差異可能導(dǎo)致了放牧強度變化趨勢存在一定的空間異質(zhì)性。青藏高原東部、南部縣區(qū)人口、牲畜數(shù)量較多且草地分布較少,多屬于重度放牧區(qū),而西北部人口、牲畜數(shù)量較少且草地分布較廣,多屬于輕度放牧區(qū)。

3.2 草地NDVI變化對氣候變化及人類活動的響應(yīng)

1982—2013年,青藏高原生長季草地NDVI整體呈增長趨勢,草地NDVI除在部分縣區(qū)呈下降趨勢外,均呈增長趨勢(圖3)。劉軍同等[16]利用NDVI研究青藏高原植被覆蓋變化時發(fā)現(xiàn)青藏高原植被覆蓋度呈“整體升高、局部退化”的趨勢,且改善區(qū)域面積遠大于退化區(qū)域面積。本研究結(jié)果與前人研究結(jié)果基本一致,青藏高原高寒草地生態(tài)系統(tǒng)呈現(xiàn)出“整體改善,局部退化”狀態(tài)[16, 18, 86-87]。

植被的生長狀態(tài)及分布會受到溫度、降水等氣候因子的影響[88]。氣候條件的改變會使植被的種類、分布發(fā)生相應(yīng)的改變以適應(yīng)新的氣候特點[7]。一些研究表明,氣溫的升高會引起植物生長期提前、生長季的加速以及生長期延長,進而會導(dǎo)致NDVI的增加[89-90]。此外,溫度和降水可以通過生物地球化學(xué)反饋影響土壤養(yǎng)分及水分的有效性,進而影響植被的生長狀況[91-92]。在本研究中,面板分析結(jié)果均表明生長季平均氣溫、降水與草地NDVI呈極顯著正相關(guān)(表2),且趨勢分析顯示1982—2013年青藏高原生長季平均氣溫、降水量均呈增長趨勢(圖3)。生長季溫度和降水量的增加可能提高了土壤養(yǎng)分和水分的有效性,加速了草地的生長進而導(dǎo)致了草地NDVI的增加[12]。魏彥強等[7]利用NDVI數(shù)據(jù)在研究青藏高原植被帶變化對氣溫、降水變化的響應(yīng)中也發(fā)現(xiàn)“暖溫化”對植被狀況的改善具有明顯的促進作用。

植被覆蓋變化除受氣候變化影響外,還與人類活動息息相關(guān)[93]。草地植被變化極易受放牧、土地利用變化等人類活動的影響。一些研究表明,放牧強度不斷增加會導(dǎo)致植被的多樣性和生物量顯著降低[69]。許岳飛等[40]以青藏高原高山嵩草草甸為對象得出隨著放牧牦牛強度的增加,草甸群落的蓋度、群落物種數(shù)、物種多樣性均顯著降低；孫大帥[53]對青藏高原東部高寒草甸植被的調(diào)查研究表明隨著放牧強度的增加,植被生物量和蓋度顯著降低；周興民等[41]在海北高寒草甸生態(tài)系統(tǒng)站長期放牧試驗表明,重度放牧組中牲畜對植被的過度啃食和踐踏作用導(dǎo)致了優(yōu)良牧草占比大幅度下降,植被地上地下生物量也隨放牧強度的增大而線性下降,與此同時,放牧踐踏的加強導(dǎo)致了土壤板結(jié)、土壤有機質(zhì)和全氮含量的下降。但也有研究表明,草原植物的凈初級生產(chǎn)力、生物多樣性及生物量并不是在不利用的情況下最高,而是在適牧條件下最高,即所謂的“優(yōu)化放牧理論”和“中度干擾理論”[94]。植物生物量及多樣性隨著放牧強度的增加會呈現(xiàn)先上升后下降的趨勢[95]。例如澤讓東科等[96]關(guān)于牦牛放牧對青藏高原草地影響的研究發(fā)現(xiàn),隨著放牧強度的增加,群落生物量及草地植被生物量會呈現(xiàn)先上升后下降的趨勢,適度放牧可以刺激草地植物的生長和提高草地生物多樣性。1982—2013年青藏高原各縣區(qū)夏季牧場放牧強度增長趨勢存在一定空間異質(zhì)性,這可能導(dǎo)致了放牧強度與生長季草地NDVI在不同的區(qū)域中呈現(xiàn)出不同的相關(guān)關(guān)系。在輕、中度放牧區(qū)中草地NDVI與放牧強度呈極顯著正相關(guān),放牧強度的適當增加可能刺激了草地植物的生長、增加了植物的生物量,進而促進了草地NDVI的增加。在重度放牧區(qū)中草地NDVI與放牧強度呈極顯著負相關(guān),長期高強度的放牧可能使得草地的生物量下降、多樣性降低,加上放牧踩踏作用的增強可能導(dǎo)致了重度放牧區(qū)草地NDVI的減少。結(jié)果與“中度干擾理論”及“放牧優(yōu)化理論”相符,即中度或適度放牧能刺激植物產(chǎn)量、改善植被生長狀態(tài)、增加植物的地上生物量。在中度或適度放牧干擾下,草地生態(tài)系統(tǒng)中植被物種的種間關(guān)系也發(fā)生改變,競爭優(yōu)勢種在生長一定程度上因牲畜啃食與踐踏作用而受到抑制,種間競爭作用減弱,競爭優(yōu)勢種和放牧忍受種共存將可能導(dǎo)致群落內(nèi)的植物多樣性與生物量的提高。但隨著放牧強度的持續(xù)增加,干擾作用的增強使得耐受性較差的物種無法生存,最終導(dǎo)致了物種多樣性的降低、生物量的減少[53]。

青藏高原植被的變化受到水資源、地形地貌、氣候及人類活動等因素的影響[10, 18]。一些研究認為氣候變化在青藏高原植被變化中起主導(dǎo)作用[97-98]。也有研究認為過度放牧、土地利用變化、濫挖濫采等人類活動才是導(dǎo)致青藏高原植被變化的主導(dǎo)因子[27, 84]。在輕度放牧區(qū)中,氣溫、降水與NDVI呈顯著正相關(guān),雖然氣溫與降水每增加1個單位對草地NDVI的貢獻小于放牧強度,但32年來氣溫與降雨增長趨勢顯著,其增加量所導(dǎo)致草地NDVI的增加量均大于放牧強度導(dǎo)致草地NDVI的增加量。在長時間序列上,降水對草地NDVI的促進作用更強,主導(dǎo)了青藏高原生長季草地NDVI的變化。在中度放牧區(qū)中,氣溫增長趨勢顯著,可能導(dǎo)致了溫度對草地NDVI促進作用的增強,而降雨增長趨勢不顯著,可能導(dǎo)致了降水對草地NDVI促進作用的減弱。在重度放牧區(qū)及整個青藏高原地區(qū),盡管放牧強度的增加會導(dǎo)致草地NDVI的減少,但32年來降雨、溫度增加所導(dǎo)致草地NDVI的增長量均大于放牧強度導(dǎo)致NDVI的減少量,因此青藏高原生長季草地NDVI呈現(xiàn)出上升趨勢。氣候因子氣溫與降水的增加均會促進草地植被增長,人類活動放牧強度的大小對植被生長則具有雙重影響,兩者長期的共同作用導(dǎo)致了青藏高原高寒草地植被的變化。

4 結(jié)論

本研究以青藏高原143個縣區(qū)32年的植被、氣候、放牧等數(shù)據(jù)為基礎(chǔ),利用趨勢分析、面板數(shù)據(jù)模型對青藏高原生長季草地NDVI、氣候及人類活動因子綜合分析得到以下結(jié)論：

(1)1982—2013年,青藏高原生長季平均氣溫與降水量整體增加,區(qū)域氣候呈現(xiàn)“暖濕化”趨勢。

(2)青藏高原生長季草地NDVI整體呈緩慢增長趨勢,草地生長狀態(tài)呈現(xiàn)出“整體改善、局部退化”趨勢。

(3)氣候因子主導(dǎo)了青藏高原生長季草地NDVI的變化。氣候變暖、降水量的增加改善了青藏高原草地植被的生長狀態(tài),促進了草地NDVI的增加。放牧強度的不斷增加會導(dǎo)致草地NDVI的減少,但其導(dǎo)致草地NDVI減少量小于氣候變化導(dǎo)致草地NDVI的增加量,因此青藏高原草地NDVI整體上呈現(xiàn)出增長趨勢。