高藝寧 王宏亮 趙萌莉

(1.內(nèi)蒙古師范大學(xué) 科技處，呼和浩特 010022；2.北京大學(xué) 城市規(guī)劃與設(shè)計(jì)學(xué)院，深圳 518055；3.內(nèi)蒙古農(nóng)業(yè)大學(xué) 草原與資源環(huán)境學(xué)院，呼和浩特 010010)

草原植被作為陸地表層復(fù)雜生態(tài)系統(tǒng)的重要組成部分，為人類提供多樣化的生態(tài)服務(wù)功能[1]。其中，植被初級(jí)凈生產(chǎn)能力(net primary productivity，NPP)具有較高的生物生產(chǎn)和碳循環(huán)功能[2-3]。草原植被NPP能夠反映草地植物群落在自然環(huán)境下的生產(chǎn)能力[4]，是表征草地植被群落生長(zhǎng)動(dòng)態(tài)和陸地植被質(zhì)量的重要參數(shù)[5]，也是衡量陸地生態(tài)系統(tǒng)應(yīng)對(duì)氣候變化的重要指標(biāo)[6]。在草地生態(tài)系統(tǒng)動(dòng)態(tài)過程中，植被NPP與氣溫、降水等氣候因子關(guān)系密切，其動(dòng)態(tài)變化反饋了植物代謝在溫度和水分條件下的變化過程[7]，較好地反映氣候特征對(duì)地表植被覆蓋狀況的影響[8]，對(duì)于區(qū)域環(huán)境具有重要的指示作用[9]。由于植被NPP在全球氣候變化與土地利用/覆被變化中的重要作用，國(guó)際地圈/生物圈計(jì)劃(IGBP)、全球氣候變化與陸地生態(tài)系統(tǒng)(GCTE)以及京都協(xié)定書均把植被NPP動(dòng)態(tài)變化列入核心研究計(jì)劃[10-11]，我國(guó)“美麗中國(guó)生態(tài)文明建設(shè)科技工程”和“人地系統(tǒng)耦合”等生態(tài)系統(tǒng)綜合性研究也將草原陸地生態(tài)系統(tǒng)NPP列為熱點(diǎn)問題[12]。因此，研究草原植被NPP時(shí)空動(dòng)態(tài)分布及其與氣候因子的關(guān)系，有助于深入分析草原陸地生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)。

尺度依賴是生態(tài)學(xué)研究的關(guān)鍵問題之一[13]，傳統(tǒng)野外調(diào)查方法僅適用于斑塊等小尺度的局地研究，難以滿足景觀區(qū)域等大尺度的綜合測(cè)量。因此，采用模型模擬植被NPP成為一種重要而又被廣泛應(yīng)用的科學(xué)方法。目前，用于測(cè)算植被NPP的模型方法可劃分為氣候統(tǒng)計(jì)模型、生態(tài)系統(tǒng)過程模型和光能利用率模型等。應(yīng)用上述方法，學(xué)者們針對(duì)氣候變化采用遙感數(shù)據(jù)對(duì)植被NPP時(shí)空變化特征進(jìn)行了廣泛研究。如侯雨等[14]基于光能利用率模型分析了1982—2000年我國(guó)不同植被類型NPP的季節(jié)性和年際性變化規(guī)律，認(rèn)為植被NPP年內(nèi)季節(jié)性變化規(guī)律明顯；劉剛等[15]采用植被凈初級(jí)生產(chǎn)力模型對(duì)2001—2014年中國(guó)植被NPP進(jìn)行估算，得到了我國(guó)植被NPP空間分布取決于植被分布和氣候條件的結(jié)論；苑全治等[16]借助植被動(dòng)態(tài)模型IBIS，對(duì)過去50年氣候變化下中國(guó)潛在植被NPP的動(dòng)態(tài)變化進(jìn)行模擬，判斷了氣候變化下植被NPP的脆弱性。圍繞我國(guó)不同區(qū)劃特點(diǎn)的植被NPP研究也較為豐富。如張雪蕾等[17]對(duì)石羊河流域植被NPP進(jìn)行估算；李登秋等[18]分析了江西省林地生態(tài)系統(tǒng)NPP時(shí)空變化特征；劉夏等[19]對(duì)三江平原濕地NPP變化進(jìn)行評(píng)價(jià)。研究發(fā)現(xiàn)，植被NPP變化研究多集中在全國(guó)、流域或省份等大尺度的空間格局分析，針對(duì)縣域等中小尺度研究相對(duì)較少[20]。研究對(duì)象主要以林地、濕地等生態(tài)系統(tǒng)為主，側(cè)重草原生態(tài)系統(tǒng)的研究相對(duì)不足。基于此，拓展荒漠草原植被NPP時(shí)空特征分析，模擬縣域尺度植被NPP對(duì)氣候變化的響應(yīng)成為本研究的主旨。

四子王旗地處內(nèi)蒙古農(nóng)牧交錯(cuò)帶，是溫帶草原向干旱荒漠過渡的典型荒漠草原區(qū)。該區(qū)植被類型多樣，草本植物豐富，為我國(guó)畜牧業(yè)生產(chǎn)提供了重要的資料來源。隨著近年來經(jīng)濟(jì)發(fā)展和人口增長(zhǎng)的不斷施壓，連續(xù)性地超載過牧、礦產(chǎn)開發(fā)以及不合理地土地利用，使得該區(qū)草原生態(tài)問題較為突出，表現(xiàn)為植被覆蓋度降低、草原荒漠化加劇以及牧草產(chǎn)量下降等。為進(jìn)一步評(píng)估草原陸地生態(tài)系統(tǒng)植被覆蓋狀況，開展典型荒漠草原植被NPP研究具有重要的應(yīng)用價(jià)值，不僅有利于草原生態(tài)修復(fù)與生產(chǎn)力提升，還為開展縣域植被NPP模擬提供了案例支持。據(jù)此，以四子王旗為研究區(qū)域，基于遙感、氣象和景觀分類數(shù)據(jù)，利用CASA模型模擬1987—2016年植被NPP的時(shí)空變化，旨在分析植被NPP與年均氣溫和年降水量的相關(guān)關(guān)系。

1 材料與方法

1.1 研究區(qū)概況

四子王旗隸屬于內(nèi)蒙古自治區(qū)烏蘭察布市境內(nèi)，地理坐標(biāo)為110°19′53″～112°59′37″ E和41°11′32″～43°22′31″ N，東與錫林郭勒盟蘇尼特右旗毗鄰，南與呼和浩特市武川縣交界，西與包頭市達(dá)爾罕茂明安聯(lián)合旗相連，北與蒙古國(guó)接壤，轄區(qū)總面積為24 036.13 km2。四子王旗地處中溫帶大陸性季風(fēng)氣候區(qū)，氣候干冷、日照充足，降水量少。年均日照3 094 h，平均氣溫4.4 ℃，年平均無霜期138 d，年平均降水量368.8 mm，主要集中在夏季，具有四季分明、雨熱同季、無霜期短等氣候特點(diǎn)。

圖1 四子王旗地理位置及景觀分類圖Fig.1 Geographical location and landscape classification map of Siziwang Banner

1.2 遙感數(shù)據(jù)與處理

1.2.1植被NDVI數(shù)據(jù)

植被NDVI來自美國(guó)國(guó)家航天航空局(NASA)全球監(jiān)測(cè)與模型研究組(GIMMS)提供的GIMMS NDVI數(shù)據(jù)和MODIS NDVI數(shù)據(jù)，數(shù)據(jù)源為http：∥ecocast.arc.nasa.gov和http：∥modis.gsfc.nasa.gov。其中，1987—2006年時(shí)間序列數(shù)據(jù)采用GIMMS NDVI半月最大合成數(shù)據(jù)，應(yīng)用ArcGIS10.5中Resample工具進(jìn)行數(shù)據(jù)重采樣；2007—2016年時(shí)間序列數(shù)據(jù)采用MODIS NDVI數(shù)據(jù)中的MOD13A2級(jí)數(shù)據(jù)，利用MODIS網(wǎng)站提供的MRT處理軟件對(duì)數(shù)據(jù)進(jìn)行投影變換、拼接裁切和格式轉(zhuǎn)換。具體是將植被NDVI數(shù)據(jù)統(tǒng)一地理坐標(biāo)系為WGS_1984，投影為Albers Equal_Area Conic，并轉(zhuǎn)換HDF格式為TIFF格式文件，得到逐月的NDVI數(shù)據(jù)，空間分辨率為 1 km×1 km。由于2種NDVI源自不同傳感器類型，可采用GIMMS NDVI和MODIS NDVI重疊期數(shù)據(jù)，基于均方根誤差RMSE和相關(guān)系數(shù)R2進(jìn)行數(shù)據(jù)一致性檢驗(yàn)[6]。本研究對(duì)2000年2月—2006年12月的2組數(shù)據(jù)進(jìn)行線性回歸，得到擬合方程YMODIS=1.211XGIMMS-0.034，其中R2=0.818，RMSE為0.228，說明數(shù)據(jù)具有顯著的一致性(P<0.01)。由于MODIS NDVI值>GIMMS NDVI值，進(jìn)而借助回歸方程對(duì)GIMMS NDVI進(jìn)行數(shù)據(jù)校正。

1.2.2氣象數(shù)據(jù)

氣象數(shù)據(jù)資料采用中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http：∥cdc.cma.gov.cn)提供的1987—2016年四子王旗及周邊市縣氣象站的月平均溫度和月降水?dāng)?shù)據(jù)，包括二連浩特、滿都拉、朱日和、達(dá)茂、化德、包頭、呼和浩特和集寧等氣象站。通過ArcGIS10.5軟件中Geostatistical Analyst模塊分別對(duì)氣象數(shù)據(jù)進(jìn)行空間插值，本研究采用Kriging插值，獲取與NDVI柵格大小一致、投影相同的數(shù)據(jù)像元，經(jīng)過掩膜裁切后獲得四子王旗氣溫和降水柵格數(shù)據(jù)。

1.2.3景觀分類

景觀分類來源于土地利用覆被類型，是依據(jù)內(nèi)蒙古自治區(qū)土地調(diào)查規(guī)劃院提供的2016年度四子王旗土地利用現(xiàn)狀數(shù)據(jù)，按照《土地利用分類標(biāo)準(zhǔn)》劃分農(nóng)田、林地、草地、水域、人工和荒漠等6類景觀。在景觀分類的基礎(chǔ)上，結(jié)合四子王旗自然特征和遙感影像，并輔助野外考察記錄進(jìn)行矢量化制圖。

1.3 研究方法

1.3.1CASA模型

CASA(Carnegie-Ames-Stanford Approach)模型是利用光能利用率原理的過程模型，通過植物吸收光合有效輻射(APAR)和實(shí)際光能利用率(ε)來計(jì)算植被NPP的模型[21]。植被吸收的光合有效輻射取決于太陽(yáng)總輻射和植被對(duì)光合有效輻射的吸收比例；實(shí)際光能利用率是植被把所吸收的光合有效輻射轉(zhuǎn)化為有機(jī)碳的效率，溫度和水分對(duì)其影響較為明顯。具體計(jì)算公式如下：

NPP(x,t)=APAR(x,t)×ε(x,t)

(1)

APAR(x,t)=SOL(x,t)×FPAR(x,t)×0.5

(2)

ε(x,t)=Tε1(x,t)×Tε2(x,t)×Wε(x,t)×εmax

(3)

式中：APAR(x,t)表示像元x在t月吸收的光合有效輻射，MJ/m2；ε(x,t)表示像元在t月的實(shí)際光能利用率，g/MJ(C)；SOL(x,t)表示像元x在t月太陽(yáng)總輻射量，MJ/m2；常數(shù)0.5表示植被所能利用的太陽(yáng)有效輻射占太陽(yáng)總輻射的比例；FPAR(x,t)為植被層對(duì)入射光合有效輻射的吸收比。SOL(x,t)采用聯(lián)合國(guó)糧食與農(nóng)業(yè)組織FAO的《The ETo calculator reference manual》[22]獲取，F(xiàn)PAR(x,t)與NDVI、植被指數(shù)SR之間存在較好的線性關(guān)系，可以通過MOD13A2產(chǎn)品提取歸一化植被指數(shù)對(duì)FPAR進(jìn)行估算[23]。εmax表示植被在理想條件下對(duì)光合有效輻射的利用效率，根據(jù)朱文泉等[24]模擬的最大光能利用率確定草地εmax值為0.542 g/MJ(C)。Tε1(x,t)表示低溫和高溫對(duì)光能利用率的影響；Tε2(x,t)表示環(huán)境氣溫從最適溫度向低溫或高溫變化時(shí)植被光能利用率逐步變小的趨勢(shì)。

(4)

(5)

式中：Topt(x,t)代表植被生長(zhǎng)最適氣溫，即選擇一年中NDVI最大月份的當(dāng)月平均氣溫；若某月均溫T(x,t)≤-10 ℃時(shí)，則Tε1(x,t)取值為0；若某月均溫T(x,t)比最適氣溫Topt(x,t)高10 ℃或低13 ℃時(shí)，則Tε2(x,t)等于最適氣溫月Topt(x,t)值的1/2。

(6)

式中：Wε(x,t)表示水分對(duì)植被光能利用率的影響，隨著環(huán)境水分的增加，Wε(x,t)將逐步增大，取值0.5～1.0。EET(x,t)是通過降水量和植被凈輻射計(jì)算的實(shí)際蒸發(fā)量；PET(x,t)是潛在蒸發(fā)量。

1.3.2植被NPP年際變化率

植被NPP采用一元線性回歸分析1987—2016年四子王旗植被NPP的時(shí)空動(dòng)態(tài)，單個(gè)像元多年回歸方程中趨勢(shì)線斜率即為年際變化率[25]。計(jì)算公式如下：

(7)

式中：θslope為趨勢(shì)斜率；n為監(jiān)測(cè)時(shí)間段的年數(shù)；NPPi為第i年的植被NPP。利用NPP序列與時(shí)間序列相關(guān)關(guān)系判斷NPP年際變化的顯著性，斜率為正表示上升，反之則表示下降。

1.3.3相關(guān)性計(jì)算

植被NPP與氣候因子相關(guān)性計(jì)算是基于植被NPP像元與年均氣溫、年降水量柵格像元的空間分析，采用相關(guān)系數(shù)法進(jìn)行公式計(jì)算。公式如下：

(8)

2 結(jié)果與分析

2.1 植被NPP空間分布特征

1987—2016年四子王旗植被NPP值為144.52 g/(m2·年)(C)，年均總量值為3.47 Tg(C)，約占內(nèi)蒙古自治區(qū)植被NPP年均值和年均總量的42.51%和1.06%。由于荒漠草原面積約為2.40×104km2，占內(nèi)蒙古土地總面積的2.03%，因此荒漠草原植被NPP年均值和年均總量明顯低于內(nèi)蒙古植被NPP的整體水平[25]。在分布上，研究區(qū)植被NPP的空間差異明顯，高值分布在南部區(qū)，低值位于北部區(qū)，呈現(xiàn)出南高北低的分布特征。由于南部屬于山地丘陵區(qū)，大量的林木、灌叢密集生長(zhǎng)，使其植被NPP明顯高于北部區(qū)荒漠沙地的植被覆蓋。另外，部分南部區(qū)域因旗政府所在地烏蘭花鎮(zhèn)及鄰近鄉(xiāng)鎮(zhèn)在城鎮(zhèn)建設(shè)和人類活動(dòng)方面的頻繁擾動(dòng)，形成了南部局地植被NPP低值分布的空間格局。在分類上(圖2)，植被NPP的一般生態(tài)特征為林地景觀>農(nóng)田景觀>草地景觀>荒漠景觀，而基于景觀分類的植被NPP年際統(tǒng)計(jì)結(jié)果為：草地景觀(2.38 Tg(C))>林地景觀(0.52 Tg(C))>農(nóng)田景觀(0.51 Tg(C))>荒漠景觀(0.01 Tg(C))，說明荒漠化草地是研究區(qū)的主導(dǎo)景觀類型，植被NPP的空間分布與景觀分類的分布面積基本一致。由于短花針茅是四子王旗草地景觀的建群種，分布面積較大，因此以短花針茅為優(yōu)勢(shì)種群的荒漠化草地在植被NPP上占據(jù)明顯優(yōu)勢(shì)。

圖2 四子王旗不同植被NPP占比情況圖Fig.2 Proportion of vegetation NPP for different types in Siziwang Banner

2.2 植被NPP年際變化特征

從總體變化上看(圖3)，1987—2016年四子王旗植被NPP年際變化呈現(xiàn)出波動(dòng)上升的趨勢(shì)。其中，最低值出現(xiàn)在1995年的3.09 Tg(C)，低于多年平均值9.65%；2016年植被NPP總量最高，約為3.69 Tg(C)，高于多年平均值7.90%；其余各年份植被NPP年際總量均處于多年平均值±6%的浮動(dòng)范圍。從年際增長(zhǎng)上看，大部分地區(qū)的植被NPP均有所增長(zhǎng)，平均增長(zhǎng)值為3.77 g/(m2·年)(C)，具有增長(zhǎng)趨勢(shì)的區(qū)域約占研究區(qū)總面積的94.23%。其中，增長(zhǎng)值>20 g/(m2·年)(C)的區(qū)域約占研究區(qū)面積比例為0.03%；增長(zhǎng)值為10～20 g/(m2·年)(C)的區(qū)域約占研究區(qū)面積比例為4.80%；增長(zhǎng)值為5～10 g/(m2·年)(C)的區(qū)域約占研究區(qū)面積比例為1.85%；增長(zhǎng)值<5 g/(m2·年)(C)的區(qū)域占研究區(qū)面積比例為87.55%。從景觀分類上看，不同景觀類型的植被NPP年際波動(dòng)差異明顯，草地景觀NPP年際波動(dòng)較大，林地景觀、農(nóng)田景觀和荒漠景觀的波動(dòng)較小。結(jié)合增減面積比來看，研究末期除林地景觀植被NPP增減面積比為2.770外，其余農(nóng)田、草地和荒漠景觀的植被NPP增減面積比分別是0.908、0.657和0.949，表明林地覆蓋面積增加明顯，其他景觀覆蓋面積略有縮減。分析發(fā)現(xiàn)，植被NPP增加的原因與當(dāng)?shù)厣鷳B(tài)保護(hù)政策密不可分，近30年的退耕還林工程、京津風(fēng)沙源治理工程以及草原荒漠化防治工程對(duì)該區(qū)生態(tài)環(huán)境的改善起到了修復(fù)作用，成為林地景觀植被NPP增加的主要原因；但人口的增長(zhǎng)和經(jīng)濟(jì)的發(fā)展，使得大量草地和部分耕地轉(zhuǎn)變成建設(shè)用地，導(dǎo)致研究區(qū)對(duì)應(yīng)的景觀植被NPP呈動(dòng)態(tài)減少的變化。

圖3 四子王旗植被NPP總量變化趨勢(shì)圖Fig.3 Trend of total vegetation NPP in Siziwang Banner from 1987 to 2016

2.3 植被NPP與氣候因子的關(guān)系

1987—2016年四子王旗年均氣溫為3.12 ～6.97 ℃，平均值為4.48 ℃，并具有北高南低的溫度分布特征；研究期內(nèi)年降水量為195.2～566.7 mm，平均值為320.5 mm，呈現(xiàn)出由南向北遞減的降雨變化趨勢(shì)?；谙嚓P(guān)性分析發(fā)現(xiàn)，研究區(qū)植被NPP與年均氣溫、年降水量均呈顯著正相關(guān)關(guān)系(P≤0.05)，相關(guān)系數(shù)分別為0.061和0.296。不同植被NPP與年均氣溫的相關(guān)性差異明顯，農(nóng)田景觀和林地景觀的相關(guān)系數(shù)為正，草地景觀和荒漠景觀相關(guān)系數(shù)為負(fù)。通過相關(guān)性空間顯示(圖4)，二者呈正相關(guān)的區(qū)域占比約為90.53%，廣泛分布于北部牧區(qū)荒漠草原；呈極顯著正相關(guān)的區(qū)域占比約為9.97%，主要分布于山地丘陵區(qū)?；诓煌脖籒PP與年降水量均呈正相關(guān)關(guān)系，在通過顯著性水平檢驗(yàn)的前提下，農(nóng)田景觀植被NPP相關(guān)系數(shù)最小(0.003)，林地景觀植被NPP相關(guān)系數(shù)最大(0.234)。通過相關(guān)性空間顯示可知，二者呈正相關(guān)的區(qū)域占研究區(qū)面積的94.48%，廣泛分布于研究區(qū)各鄉(xiāng)鎮(zhèn)、蘇木區(qū)域。其中，南部丘陵區(qū)域植被NPP與年降水量呈極顯著正相關(guān)，占研究區(qū)面積的12.10%，呈極顯著負(fù)相關(guān)的區(qū)域僅分布于紅格爾蘇木中心，占比約為1.53%。研究發(fā)現(xiàn)，年降水量是影響研究區(qū)植被NPP的主要?dú)夂蛞蜃?。一方面由于氣候因子所引發(fā)的雨熱分異直接影響植被NPP的空間分布，另一方面不同植被的生長(zhǎng)特性對(duì)生境條件的反饋差異，引起不同景觀植被NPP對(duì)氣溫和降水量的差異化響應(yīng)。

圖4 四子王旗1987—2016年植被NPP與年均氣溫、年降水量的關(guān)系Fig.4 Correlation of annual NPP and annual rainfall and annual temperature from 1987 to 2016 in Siziwang Banner

3 討論與結(jié)論

3.1 討論

采用CASA模型模擬了1987—2016年荒漠草原植被NPP的時(shí)空變化，并與氣候因子進(jìn)行了相關(guān)性分析，研究發(fā)現(xiàn)荒漠草原植被NPP低于內(nèi)蒙古植被NPP的整體水平，且空間分布具有南高北低的變化特征，這與韓芳[26]的研究結(jié)果一致，但植被NPP的明顯差異與不同研究單元的自然條件和覆被狀況存有較大關(guān)系。李剛等[27]發(fā)現(xiàn)植被NPP與水熱條件緊密相關(guān)，尤其與降水量關(guān)系更為密切，這與本研究結(jié)果相一致。原因是研究區(qū)處于典型干旱區(qū)，以高溫少雨為特征的干暖化氣候使得降水量成為制約植被NPP增長(zhǎng)的主導(dǎo)因子。另外，本研究受不確定性因素和模型變量的影響較為明顯。就不確定性因素而言，由于數(shù)據(jù)獲取源的限制，NDVI數(shù)據(jù)空間精度為1 km×1 km，對(duì)于縣域尺度的空間模擬略顯不足；植被NPP測(cè)算與研究者的認(rèn)知水平關(guān)系較大，其模型選擇、參數(shù)設(shè)定和模擬過程始終存在人為判定影響。在模型變量方面，本研究?jī)H采用年均氣溫和年降水量來表征氣候因子，而氣候條件還包括太陽(yáng)輻射、風(fēng)蝕、日照以及氣壓等因素，存有氣候因子變量選擇的局限性；由于植被NPP時(shí)空變化還與土壤條件和地形地貌等地理環(huán)境因素相關(guān)，未來針對(duì)研究區(qū)植被NPP的相關(guān)性分析有待進(jìn)一步深入。

3.2 結(jié)論

基于遙感數(shù)據(jù)、氣象數(shù)據(jù)和景觀分類數(shù)據(jù)，通過CASA模型模擬了1987—2016年荒漠草原植被NPP的時(shí)空變化及其與年均氣溫和年降水量的相關(guān)關(guān)系，得到如下結(jié)論：

1)1987—2016年四子王旗年均植被NPP為144.52 g/(m2·年)(C)，植被NPP年均總量為3.47 Tg(C)。從空間分布上看，研究區(qū)植被NPP的空間差異明顯，表現(xiàn)為南高北低的分布特征；從分類統(tǒng)計(jì)上看，草地景觀>林地景觀>農(nóng)田景觀>荒漠景觀，說明荒漠化草地是研究區(qū)的主導(dǎo)景觀類型。

2)研究區(qū)植被NPP年際變化呈波動(dòng)上升的趨勢(shì)，波動(dòng)為3.09～3.69 Tg(C)。其中，約有94.23%地區(qū)植被NPP呈增長(zhǎng)趨勢(shì)，平均增長(zhǎng)值為3.77 g/(m2·年)(C)；另有，農(nóng)田、林地、草地和荒漠植被NPP增減面積比分別為0.908、2.770、0.657和0.949，表明植被NPP既受生態(tài)保護(hù)政策的影響，又與當(dāng)?shù)亟?jīng)濟(jì)社會(huì)發(fā)展相關(guān)。

3)本研究植被NPP與年均氣溫和年降水量均呈顯著正相關(guān)關(guān)系，相關(guān)系數(shù)分別為0.061和0.296；基于相關(guān)性空間顯示結(jié)果，呈正相關(guān)的分布面積所在比例分別為90.53%和94.48%，表明降水量是影響植被NPP的主要?dú)夂蛞蜃印?/p>