周 婧，何政偉，張 志，陳 林

（1.成都理工大學(xué) 地球科學(xué)學(xué)院，四川 成都 610059；2.成都理工大學(xué) 旅游與城鄉(xiāng)規(guī)劃學(xué)院，四川 成都 610059）

植被凈初級(jí)生產(chǎn)力（NPP）是指綠色植物在單位時(shí)間和單位面積上所產(chǎn)生的有機(jī)物質(zhì)總量扣除自養(yǎng)呼吸后剩余的部分，是表征植物活動(dòng)的關(guān)鍵因子，可作為評(píng)價(jià)生態(tài)系統(tǒng)結(jié)構(gòu)功能和生物圈承載力的重要指標(biāo)，用以分析區(qū)域生態(tài)環(huán)境本底狀況[1-3]。

1 研究區(qū)概況

川西北江河源區(qū)地處四川省西北部，青藏高原東部，由若爾蓋、石渠、紅原、色達(dá)、阿壩、甘孜、壤塘7個(gè)縣組成，總面積為8.02萬(wàn)km2，總?cè)丝诩s為40萬(wàn)，主要地形為谷地和淺丘。研究區(qū)位于黃河、長(zhǎng)江兩大水系的上游，是重要的水源涵養(yǎng)功能區(qū)和生態(tài)保護(hù)區(qū)，主要?dú)夂蝾?lèi)型為高原亞寒帶氣候以及濕潤(rùn)、半濕潤(rùn)氣候，日照較長(zhǎng)，晝夜溫差大，年均氣溫低[4]。植被以高山草甸為主，森林面積極少，是我國(guó)重要的牧區(qū)之一。由于過(guò)度放牧、氣候條件嚴(yán)峻、地勢(shì)較高等原因?qū)е缕渖鷳B(tài)環(huán)境脆弱性加劇，其水源涵養(yǎng)功能和生物多樣性較易受到破壞，因此實(shí)時(shí)了解其植被生長(zhǎng)狀況對(duì)認(rèn)清該區(qū)域生態(tài)現(xiàn)狀和保護(hù)生態(tài)環(huán)境具有重要意義。

2 NPP反演

2.1 數(shù)據(jù)預(yù)處理

本文選擇MODIS數(shù)據(jù)產(chǎn)品，通過(guò)MRT工具進(jìn)行投影變換、重采樣和掩膜。氣象數(shù)據(jù)選用月均溫度、降水量、相對(duì)濕度和日照百分率等。

2.2 反演過(guò)程

本文利用CASA模型計(jì)算NPP，該模型屬于光能利用率模型，具體計(jì)算過(guò)程為[5-12]：

1）計(jì)算區(qū)域潛在蒸散量Ep，計(jì)算公式為：

2）根據(jù)Ep和月降水量P計(jì)算區(qū)域蒸散量E和月太陽(yáng)凈輻射Rn，計(jì)算公式為：

3）根據(jù)Ep計(jì)算表征植物可利用的有效水分含量的影響程度，計(jì)算公式為：

4）根據(jù)T計(jì)算溫度影響系數(shù)，最大光能利用率εmax使用朱文泉的模擬值。

5）根據(jù)相對(duì)濕度R、溫度t計(jì)算地面水汽壓e，es0=0.107 8 Pa，a、b為常數(shù)。

6）通過(guò)e、日照百分率S計(jì)算最大晴天總輻射量Q0，進(jìn)一步求出 PAR。a=0.229，b=0.679，C1、C2、C3為常數(shù)，φ為緯度，H為海拔。

7）根據(jù)NDVI計(jì)算FPAR，計(jì)算公式為：

8）根據(jù)FPAR計(jì)算NPP，公式為：

3 NPP時(shí)空變化特征

3.1 NPP時(shí)間特征

2006～2009年NPP逐月均值分布如圖1所示，可以看出，研究區(qū)NPP呈較強(qiáng)的季節(jié)性（表1）。以2009年為例，3～5月（春季）水熱條件較差，NPP維持低水平，增幅極低，NPP積累量占全年的16.2%；隨著水熱條件的提高，6～8月（夏季）NPP月積累量急速上升，在7、8月積累量達(dá)到最高值，NPP積累量占全年的63.2%；秋季水熱條件開(kāi)始下降，NPP月積累量開(kāi)始快速下降，NPP積累量占全年的19.3%；12月～次年2月（冬季）是一年中NPP積累量最少的時(shí)段，NPP積累量占全年的1.2%?？傮w來(lái)看，2006～2009年各季NPP所占百分比無(wú)明顯變化，但夏季NPP占比略有下降，春、秋季占比略有增加，可能是由于氣候變暖，導(dǎo)致植被物候期延長(zhǎng)所致。

分析研究區(qū)NPP逐年變化情況可知，2006～2009年NPP年均值在297.56～314.56 gCm-2a-1范圍內(nèi)變化，波動(dòng)不大，其中2008年最低，可能是由于該年水熱條件降低所致。

圖1 研究區(qū)月均NPP變化趨勢(shì)圖

表1 各季NPP占比和年NPP總量

3.2 NPP空間特征

分析研究區(qū)植被NPP空間分布特征可知，NPP空間分布存在明顯差異（圖2），由東向西逐步遞減，可能是因?yàn)闁|部海拔較低，地勢(shì)平坦，而西部地區(qū)海拔較高，且多高山，不利于植被生長(zhǎng)而造成的。同時(shí)，比較各樣點(diǎn)逐月NPP發(fā)現(xiàn)，自東向西NPP峰值出現(xiàn)時(shí)段前移，由8月變?yōu)?月，東部地區(qū)NPP峰值主要集中于8月，而西部地區(qū)NPP峰值基本集中于7月（圖3），這可能與研究區(qū)東西部地形、高程有關(guān)。

圖2 2006～2009年NPP空間分布

圖3 研究區(qū)各樣點(diǎn)逐月NPP變化趨勢(shì)

本文從海拔、坡度、坡向3個(gè)角度來(lái)分析研究區(qū)NPP變化情況，如圖4所示，可以看出，研究區(qū)平均海拔較高，大部分位于3 000～5 000 m范圍內(nèi)，高程自西向東呈遞減趨勢(shì)。在2 503～4 100 m，NPP隨海拔的升高略有下降，但降幅較小；在4 100 m處出現(xiàn)一個(gè)拐點(diǎn)，下降速度明顯加快；在5 440 m處NPP降到最低值，表明NPP與高程可能存在拮抗關(guān)系。分析研究區(qū)NPP與坡度的關(guān)系可知，東、西部地區(qū)的坡度未表現(xiàn)出明顯的差異，地形起伏度不大，超過(guò)90%的地區(qū)位于0～15°范圍內(nèi)；在0～25°范圍內(nèi)，NPP隨坡度的升高略有增長(zhǎng)，但變化不明顯；在25～28.6°范圍內(nèi)，NPP隨坡度的升高快速增長(zhǎng)，達(dá)到404 gCm-2a-1，表明研究區(qū)內(nèi)植被NPP隨坡度升高而增長(zhǎng)。分析研究區(qū)不同坡向的NPP變化趨勢(shì)可知，NPP在東、東南方向達(dá)到最高值，再沿順時(shí)針?lè)较蜻f減，在西、西北方向達(dá)到最低值，然后又沿順時(shí)針?lè)较蜻f增，這可能是由于研究區(qū)處于北半球，所以東、東南坡日照條件較好，植被NPP較高。

研究區(qū)東、西部地區(qū)不同坡度、坡向分布均勻，未見(jiàn)明顯地域差異，但高程則呈明顯的西高東低區(qū)域特征，因此高程差異可能是導(dǎo)致研究區(qū)NPP東高西低以及自東向西NPP峰值時(shí)段前移的重要因子之一。

圖4 海拔、坡度、坡向與NPP的關(guān)系

3.3 NPP區(qū)域增減分析

本文基于像元一元線性回歸方法計(jì)算研究區(qū)NPP絕對(duì)年際變化率，計(jì)算公式為：

式中，θslope為像元在一段時(shí)間內(nèi)一元線性回歸方程的斜率；n為年數(shù)；NPPi為像元點(diǎn)在第i年的NPP。θslope＞0表明 NPP呈增加趨勢(shì)，θslope＜0表明 NPP呈減少趨勢(shì)，|θslope|≈0表明NPP較穩(wěn)定。

標(biāo)準(zhǔn)差分級(jí)法可較好地表現(xiàn)平均值與要素屬性之間的差異，因此本文采用該方法來(lái)對(duì)NPP年際變化率進(jìn)行分級(jí)，共分為7級(jí)（圖5），其中重度NPP減少占研究區(qū)總面積的0.71%、中度NPP減少占5.53%、輕度NPP減少占22.57%、相對(duì)穩(wěn)定區(qū)域占42.43%、輕度NPP增加占20.80%、中度NPP增加占7.19%、重度NPP增加占0.76%。研究區(qū)NPP增長(zhǎng)與減少面積相差不大，但NPP增長(zhǎng)和減少區(qū)域呈明顯的空間分布差異，NPP增加區(qū)域主要集中于地形起伏較大、海拔較高的西部高原山區(qū)，減少區(qū)域則主要出現(xiàn)在較平坦的、易受人類(lèi)活動(dòng)干擾的中部地區(qū)，而東部NPP總體較穩(wěn)定。

圖5 川西北江河源區(qū)2006～2009年NPP增減格局圖

4 NPP對(duì)氣候變化的響應(yīng)

4.1 氣象因子影響的滯后性

氣象因子對(duì)NPP具有重要影響，但其影響具有一定滯后性，因此存在上月氣象因子對(duì)當(dāng)月NPP的貢獻(xiàn)大于當(dāng)月氣象因子對(duì)NPP的貢獻(xiàn)的可能。以2009年8月為例，在研究區(qū)內(nèi)均勻提取樣點(diǎn)，樣點(diǎn)NPP和氣象數(shù)據(jù)如表2所示，將樣點(diǎn)8月NPP分別與7、8月的氣象因子進(jìn)行線性及多項(xiàng)式擬合，分析過(guò)程中去除兩個(gè)明顯異常樣點(diǎn)，以得到更準(zhǔn)確的結(jié)果，結(jié)果見(jiàn)表3。線性和多項(xiàng)式擬合結(jié)果均反映8月NPP與8月溫度建立的回歸方程的R2均低于0.2，但與7月溫度建立的回歸方程的R2均高于0.5，說(shuō)明氣象因子對(duì)NPP的影響具有一定滯后性。7、8月降水量和溫度與8月NPP線性和多項(xiàng)擬合結(jié)果表明，8月NPP主要受7月降水量的影響。

4.2 相關(guān)性分析

溫度和降水量是植被NPP變化的主要驅(qū)動(dòng)因子，如圖6所示，雖然研究區(qū)的月降水量出現(xiàn)雙峰現(xiàn)象，但降水量和氣溫總體呈正態(tài)分布，且二者的時(shí)間變化趨勢(shì)相似。2～8月溫度和降水量均逐漸上升，一般在7、8月達(dá)到峰值，9月開(kāi)始下降，在12月、次年1月降到最低，呈明顯的季節(jié)性特征，夏、秋季降水量較多同時(shí)伴有高溫，春、冬季則為低溫少雨，即研究區(qū)NPP變化具有較強(qiáng)的季節(jié)性差異。分析NPP與溫度、降水量的變化趨勢(shì)可知，從時(shí)間上分析，研究區(qū)月均NPP與月均溫度和降水量的變化趨勢(shì)幾乎一致；從空間上分析，NPP與溫度、降水量的空間分布特征也較一致，由東向西減弱，說(shuō)明溫度和降水量對(duì)NPP的增長(zhǎng)有促進(jìn)作用，但由于氣象因子對(duì)植被NPP影響的滯后性，NPP峰值出現(xiàn)的時(shí)間比溫度、降水量峰值出現(xiàn)時(shí)間約晚一個(gè)月。

表2 NPP和氣象因子樣點(diǎn)數(shù)據(jù)

表3 NPP與氣象因子的線性、多項(xiàng)式回歸結(jié)果

為了更清楚地了解各氣象因子與植被NPP的相關(guān)性，計(jì)算月均溫度、月降水量與月均NPP之間的相關(guān)系數(shù)，如表4所示，可以看出，NPP與溫度、降水量的相關(guān)系數(shù)均較高，其中NPP與降水量的相關(guān)系數(shù)大于0.8，屬高度相關(guān)；NPP與溫度的相關(guān)系數(shù)為0.79，相關(guān)性較強(qiáng)，這一結(jié)果與NPP和氣象因子的擬合結(jié)果一致，因此降水量是研究區(qū)NPP變化的主要驅(qū)動(dòng)因子。

圖6 2006～2009年NPP與氣象因子逐月均值

表4 NPP與氣象因子的相關(guān)系數(shù)

5 結(jié) 語(yǔ)

本文利用MODIS衛(wèi)星數(shù)據(jù)和氣象資料進(jìn)行川西北江河源區(qū)的植被NPP反演，并分析了研究區(qū)NPP的空間分布特征、時(shí)間變化情況以對(duì)溫度、降水量?jī)蓚€(gè)氣象因子的響應(yīng)情況，得到的結(jié)果如下：

1）從時(shí)間上分析，研究區(qū)NPP呈正態(tài)分布，具有明顯的季節(jié)性特征，夏、秋季積累量較多，超過(guò)全年積累量的80%，春、冬季積累量較少，不到全年積累量的20%。2006～2009年，研究區(qū)夏、冬季積累量略有下降，春、秋季積累量略有上升，這一現(xiàn)象可能與研究區(qū)水熱條件的變化有關(guān)。

2）2006～2009年研究區(qū)NPP均值在297.56～314.56 gCm-2a-1之間變動(dòng)，空間分布差異明顯，整體上呈由東向西遞減的趨勢(shì)；且由樣點(diǎn)NPP逐月變化情況可知，研究區(qū)自東向西存在NPP峰值前移的現(xiàn)象。分析研究區(qū)NPP隨高程、坡度、坡向的變化情況可知，高程可能是研究區(qū)NPP分布空間差異和NPP峰值前移的主要影響因子。

3）2006～2009年研究區(qū)NPP增加和減少面積相差不大，但NPP增加區(qū)域主要位于受人類(lèi)活動(dòng)影響較少的西部高原、高山地區(qū)，而NPP減少區(qū)域主要位于易受人類(lèi)活動(dòng)影響的、地形起伏度較低的中部地區(qū)。

4）研究區(qū)月均溫度和降水量時(shí)空變化動(dòng)態(tài)基本一致，從時(shí)間角度來(lái)看，NPP與氣溫、降水量均呈正態(tài)分布，且表現(xiàn)出明顯的季節(jié)性差異；溫度、降水量增加時(shí)NPP也增加，說(shuō)明這兩個(gè)氣象因子對(duì)NPP具有協(xié)同效應(yīng)，但由于氣象因子影響的滯后性，NPP峰值出現(xiàn)時(shí)間比氣象因子峰值出現(xiàn)時(shí)間推遲一個(gè)月。從空間角度來(lái)看，NPP與溫度、降水量總體上均呈東高西低的變化趨勢(shì)。

5）將2009年8月NPP分別與2009年7、8月的氣象因子進(jìn)行線性、多項(xiàng)式擬合發(fā)現(xiàn)，8月NPP與7月降水量建立的回歸方程的R2高于與本月降水量建立的回歸方程的R2，說(shuō)明當(dāng)月NPP主要受上月氣象因子的影響。

6）通過(guò)計(jì)算月均NPP與月均降水量、溫度的相關(guān)系數(shù)發(fā)現(xiàn)，NPP與溫度、降水量的相關(guān)系數(shù)均較高，說(shuō)明這兩個(gè)氣象因子對(duì)NPP變化均有較大影響，且降水量的相關(guān)系數(shù)略高于溫度，表明降水量是研究區(qū)NPP變化最主要的驅(qū)動(dòng)因子。