梁爽 張吳平 畢如田 高培霞 王國(guó)芳

摘 要： 為科學(xué)有效地管理和調(diào)控植被資源，解決水土流失、植被減少等生態(tài)問題，該研究依據(jù)山西省2005—2015年MOD17A3H數(shù)據(jù)，利用ARCGIS、ENVI等軟件，運(yùn)用統(tǒng)計(jì)學(xué)分析方法，揭示了山西省植被NPP時(shí)空分布變化規(guī)律及對(duì)氣候、人為等影響因素的響應(yīng)特征。結(jié)果表明：山西省植被NPP平均值為326.5 g（C）·m-2·a-1，其中草地、耕地、灌叢和林地的NPP多年平均值依次為300.3、353.6、366.5和390.1 g（C）·m-2·a-1；植被NPP總體波動(dòng)增大、變化為顯著、極顯著、顯著增大區(qū)域面積比例達(dá)56.33%，集中在山西省西部；植被NPP極顯著、顯著減小區(qū)域面積集中在山西省東南角，占比為2.22%；草地NPP變化速率最大，耕地大于灌叢，林地最??；植被NPP平均值和降水平均值之間表現(xiàn)為呈顯著正相關(guān)。基于柵格單元值計(jì)算，全省17.01%的區(qū)域植被NPP與降水之間表現(xiàn)為顯著或極顯著正相關(guān)，集中在山西省北部；全省3.66%的區(qū)域植被NPP與氣溫之間表現(xiàn)為顯著或極顯著負(fù)相關(guān)，集中在山西省中部。這表明2005—2015年，山西省植被NPP總體呈好轉(zhuǎn)趨勢(shì)；不同植被對(duì)人類活動(dòng)及環(huán)境變化的響應(yīng)有所差異，草地、耕地生態(tài)結(jié)構(gòu)穩(wěn)定性較弱，NPP變化明顯，灌叢和林地穩(wěn)定性較強(qiáng)，NPP數(shù)值穩(wěn)定；植被NPP與降水之間呈顯著正相關(guān)，與氣溫之間呈負(fù)相關(guān)，氣候因子整體上促進(jìn)植被NPP增大，人為因素整體上抑制植被NPP增大。

關(guān)鍵詞： 山西省， 植被凈初級(jí)生產(chǎn)力， 土地利用類型， MODIS， 驅(qū)動(dòng)力

中圖分類號(hào)： Q948 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1000-3142（2018）08-1005-10

Abstract： In order to solve the ecological problems like soil erosion， vegetation reduction by scientifically and effectively management of vegetation resources. In this paper， the data of MOD17A3H in Shanxi Province were used to analyze the temporal and spatial distribution of vegetation net primary productivity （NPP） in Shanxi Province during 2005—2015 by using the software such as ARCGIS and ENVI. The climate， land use types and elevation data were also used to analyze the response characteristics of vegetation NPP to climate and other influencing factors by using statistical method. The results were as follows： The average vegetation NPP of Shanxi Province during 2005—2015 was 326.5 g（C）·m-2·a-1， and the average NPP of grassland， farmland， shrubs and forestland during 2005—2015 were 300.3， 353.6， 366.5 and 390.1 g（C）·m-2·a-1； The overall fluctuation of NPP in vegetation in Shanxi Province increased significantly with a significant or extremely significant increase in the proportion of regional area reached 56.33%， mainly concentrated in the western part of Shanxi Province； The significant or extremely significant reduction of vegetation NPP area was mainly concentrated in the southeastern corner of Shanxi Province， accounting for 2.22% of the total area； The rate of change of NPP in grassland was the largest， the farmland was greater than the shrub， and the forest land was the smallest； There was a significant positive correlation between the average of vegetation NPP and the average of precipitation in Shanxi Province. Based on calculated pixel value， the area of vegetation NPP and precipitation had significant or extremely significant positive relation in the total study area， accounting for 17.01%， mainly concentrated in northern Shanxi Province； The area of vegetation NPP and temperature had significant and extremely significant negative correlation relation in the total study area， accounting for 3.66%， mainly concentrated in central Shanxi Province. The results indicated that the vegetation NPP of Shanxi Province during 2005—2015 was improving； The response of different vegetation to human activities and environmental factors was different， the stability of grassland and farmland was weak， and the NPP changes obviously， the stability of shrub and forestland was strong， and NPP was stable； There was a significant positive correlation between vegetation NPP and precipitation， and a negative correlation between vegetation NPP and temperature， as a whole， the climate factor promotes the rise of NPP， while the human factor is mainly manifested as the inhibitory effect on NPP changes.

Key words： Shanxi Province， vegetation net primary productivity （NPP）， land use type， MODIS， driving factor

植被凈初級(jí)生產(chǎn)力（net primary productivity，NPP）是指單位面積的植被在單位時(shí)間內(nèi)憑借自身光合作用產(chǎn)生的有機(jī)物扣去自養(yǎng)呼吸消耗后余下的部分（Lieth & Whittaker，1975），體現(xiàn)了植被將無機(jī)碳固定、轉(zhuǎn)化為有機(jī)碳的效率。NPP是植被自身和其他生物在生態(tài)系統(tǒng)中生存、發(fā)育及繁衍的物質(zhì)基礎(chǔ)，它是植被自身和外界環(huán)境因素間相互作用、相互影響所產(chǎn)生的結(jié)果，是評(píng)價(jià)生態(tài)環(huán)境優(yōu)劣的重要指標(biāo)之一。因此，分析各種影響因素與植被NPP之間的相關(guān)性，進(jìn)而改善生態(tài)、推動(dòng)可持續(xù)發(fā)展的研究在生態(tài)問題日益嚴(yán)峻的今天顯得尤為重要。自Ebermayer于1876年開始研究植被初級(jí)凈生產(chǎn)力以來，不同國(guó)家的專家學(xué)者通過大量的實(shí)驗(yàn)研究，提出了各種計(jì)算和測(cè)量植被NPP的方法，在模型模擬植被NPP方面取得了一定成果。自20世紀(jì)80年代起，專家們?yōu)榱藦浹a(bǔ)植被NPP經(jīng)驗(yàn)?zāi)Ｐ偷娜毕荩芯砍隽嘶谥脖簧锾卣鞯闹脖籒PP過程模型。隨著3S（地理信息科學(xué)、全球定位系統(tǒng)、遙感）技術(shù)的高速發(fā)展，植被NPP過程模型得到了充分的普及與應(yīng)用（孫力煒，2013）；20世紀(jì)90年代后期，以孫睿和朱啟疆（2001）、樸世龍等（2001）、陳利軍等（2002）為代表的國(guó)內(nèi)學(xué)者們?cè)贚andsat、MODIS等數(shù)據(jù)大規(guī)模進(jìn)入國(guó)內(nèi)的基礎(chǔ)上大力促進(jìn)了國(guó)內(nèi)關(guān)于植被NPP的分析研究。

隨著技術(shù)日趨成熟、數(shù)據(jù)日益完善，植被NPP的研究也日漸成熟，在改善生態(tài)環(huán)境、推動(dòng)可持續(xù)發(fā)展方面發(fā)揮著重要作用。山西省地處中國(guó)北部，屬溫帶大陸性季風(fēng)氣候，是典型的黃土高原，境內(nèi)有呂梁、太行兩座山脈，植被類型豐富。山西省礦產(chǎn)資源豐富，近年隨著經(jīng)濟(jì)發(fā)展，人類活動(dòng)加強(qiáng)，資源開采過度使得土壤污染、植被減少、水土流失等生態(tài)問題日益嚴(yán)峻，潛在的資源得不到科學(xué)合理地開發(fā)利用，經(jīng)濟(jì)可持續(xù)發(fā)展受到嚴(yán)重影響，植被NPP作為評(píng)價(jià)生態(tài)系統(tǒng)的重要指標(biāo)，而對(duì)于山西省植被NPP狀況的研究相對(duì)匱乏。因此，本研究利用山西省2005—2015年MO17A3H數(shù)據(jù)，綜合氣溫、降水等數(shù)據(jù)，運(yùn)用統(tǒng)計(jì)學(xué)分析方法，揭示山西省植被NPP時(shí)空分布格局及其變化規(guī)律，分析植被NPP在不同環(huán)境條件（區(qū)域、植被類型、高程）下對(duì)氣候變化及人類活動(dòng)的響應(yīng)特征，為山西省植被保護(hù)和合理開發(fā)利用提供科學(xué)依據(jù)。

1 研究區(qū)概況

山西省位于中國(guó)北部，總面積15.67 萬km2，是典型的黃土高原地形，山區(qū)面積比例達(dá)到80%以上，地勢(shì)上東北高西南低，境內(nèi)大部分地區(qū)高程超過1 500 m。山西省屬于溫帶大陸性季風(fēng)氣候，年平均氣溫最低為3 ℃，最高為14 ℃，晝夜溫度差異明顯，南北溫度差異明顯。全省冬季日平均氣溫均在0 ℃以下，夏季全省氣溫較高，7月份日平均氣溫介于21～26 ℃之間。全省年總降水量為400～650 mm，降雨量不平均，6月份至8月份雨多且集中，降水量達(dá)到全年總降水量的60%，全省地形地貌復(fù)雜多變，山區(qū)眾多，對(duì)降雨的影響很大，不同地區(qū)降水差異性顯著。

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

2.1 數(shù)據(jù)來源及處理

NPP數(shù)據(jù)為美國(guó)NASA EOS/MODIS 2005-2015年的MOD17A3H數(shù)據(jù)，其利用MODIS衛(wèi)星數(shù)據(jù)結(jié)合模型運(yùn)算出NPP數(shù)據(jù)，精度高、易操作，空間分辨率為500 m，時(shí)間分辨率為1 a。

高程數(shù)據(jù)來源于中國(guó)科學(xué)院地理空間數(shù)據(jù)云（http：//www.gscloud.cn），空間分辨率為30 m。

土地利用類型數(shù)據(jù)為2010年美國(guó)MODIS的MCD12Q1數(shù)據(jù)，空間分辨率達(dá)到500 m。

氣象數(shù)據(jù)來自中國(guó)氣象科學(xué)數(shù)據(jù)網(wǎng)（http：//cdc.cma.gov.cn），采用反距離權(quán)重法（IDW）對(duì)2005—2015年山西省及其周邊?。ㄊ校?0個(gè)氣象站點(diǎn)逐年的溫度、降水?dāng)?shù)據(jù)進(jìn)行插值運(yùn)算。

為避免因數(shù)據(jù)之間坐標(biāo)、分辨率等屬性不同而導(dǎo)致誤差，利用ARCGIS軟件將所有數(shù)據(jù)統(tǒng)一轉(zhuǎn)換為WGS84/UTM坐標(biāo)系、500 m分辨率后，由ARCGIS和ENVI軟件進(jìn)行后續(xù)的研究分析。

2.2 研究方法

2.2.1 植被NPP變化趨勢(shì)研究 采用線性回歸方法研究2005—2015年山西省每個(gè)柵格單元內(nèi)植被NPP隨時(shí)間變化趨勢(shì)（穆少杰等，2013；朱鋒等，2010；吳珊珊等，2016）。

將檢驗(yàn)結(jié)果分為極顯著負(fù)相關(guān)、顯著負(fù)相關(guān)、相關(guān)性不顯著、顯著正相關(guān)和極顯著正相關(guān)5種類型（李曉榮等，2017）。

2.2.3 植被NPP分區(qū)統(tǒng)計(jì)研究 基于高程數(shù)據(jù)將研究區(qū)域劃分成不同梯度區(qū)域，結(jié)合土地利用類型和植被NPP年平均值的數(shù)據(jù)，對(duì)該區(qū)域植被NPP分類研究，分析高程梯度上植被NPP分布變化特征及人為、氣候等因素與植被NPP間的相關(guān)性。

3 結(jié)果與分析

3.1 NPP空間分布特征

利用山西省MOD17A3H數(shù)據(jù)，根據(jù)土地利用類型數(shù)據(jù)進(jìn)行分類研究，發(fā)現(xiàn)山西省植被NPP 2005—2015年間平均值為326.5 g（C）·m-2·a-1，其中93.53%的區(qū)域植被NPP值在200～500 g（C）·m-2·a-1之間，草地NPP平均值為300.3 g（C）·m-2·a-1，耕地NPP平均值為353.6 g（C）·m-2·a-1，灌叢NPP平均值為366.5 g（C）·m-2·a-1，林地NPP平均值為390.1 g（C）·m-2·a-1。

從圖1可以看出，山西省植被NPP空間分布整體上呈現(xiàn)北低南高，西低東高的特征。結(jié)合自然條件、地形地貌及人為活動(dòng)等因素分析可知，山西省水熱等自然條件自北向南、自西向東逐漸變好，境內(nèi)呂梁、太行兩座山脈主要分布在東部和南部，山區(qū)植被覆蓋度高，受人為因素影響較小，植被NPP數(shù)值高；山西省煤炭資源豐富，境內(nèi)礦區(qū)眾多，而西部和北部礦區(qū)多為露天礦區(qū)，其開采方式會(huì)對(duì)植被造成一定程度的破壞，植被覆蓋度較低，所以植被NPP值較低。因此山西省植被NPP總體上呈現(xiàn)北低南高，西低東高的分布特征。

3.2 植被NPP年際變化分析

從圖3可以看出，2005—2015年，山西省植被NPP總體波動(dòng)增大，不同類型植被NPP的年際變化趨勢(shì)和山西省植被NPP總體變化趨勢(shì)大體相同，表明是某種大范圍影響因素影響了山西省植被NPP總體變化趨勢(shì)。

從圖4可以看出，2005—2015年，山西省86.07%的區(qū)域植被NPP年際變化率介于0～15 g（C）·m-2·a-1之間，87.73%的區(qū)域植被NPP年際變化率為正值，表明山西植被NPP總體呈增大趨勢(shì)。山西省中西部地區(qū)植被NPP增大較快，植被NPP減小區(qū)域則集中分布在東南部。從圖5可以看出，植被NPP整體顯著改變，植被NPP顯著增大及極顯著增大區(qū)域面積占總面積的56.33%，主要集中在山西省西部；植被NPP顯著減小以及極顯著減小區(qū)域主要集中在山西省東南角。草地NPP顯著或者極顯著增大區(qū)域面積占草地區(qū)域總面積比例最大，為69.20%，遠(yuǎn)高于其他類型植被；林地NPP顯著或者極顯著減小區(qū)域面積占林地區(qū)域總面積比例最大，為7.56%?？傮w上林地和灌叢生態(tài)系統(tǒng)結(jié)構(gòu)穩(wěn)定，對(duì)環(huán)境變化響應(yīng)不顯著，而以草本科植物為主構(gòu)成的草地、耕地生態(tài)系統(tǒng)結(jié)構(gòu)相對(duì)簡(jiǎn)單，受環(huán)境變化影響顯著，當(dāng)環(huán)境變好之后NPP會(huì)顯著增大。

從圖4和圖5可以看出，植被NPP快速增大的區(qū)域，其植被NPP多表現(xiàn)為顯著或極顯著增大，表明研究期間該區(qū)域植被狀況明顯好轉(zhuǎn)；植被NPP減小的區(qū)域，其植被NPP多表現(xiàn)為變化不顯著，表明研究期間植被NPP減小區(qū)域的植被無明顯惡化。植被NPP顯著或極顯著增大區(qū)域以草地與耕地居多，而顯著或極顯著減小區(qū)域則以林地居多，NPP減小區(qū)域大部分分布在太行山內(nèi)，該區(qū)域內(nèi)植被NPP較低且植被NPP 變化速率也較小，表明該區(qū)域長(zhǎng)期受到負(fù)面因素影響。植被NPP顯著減小區(qū)域面積小、分布零散，而氣候因素影響廣泛，表明氣候因素不是造成此部分植被變差的主要原因，結(jié)合太行山近些年開發(fā)歷史分析，推測(cè)人為因素是部分植被NPP顯著減小的重要原因之一。

3.3 植被NPP和氣候因素之間相關(guān)性分析

基于研究區(qū)域平均值進(jìn)行計(jì)算可知，2005—2015年山西省降水與植被NPP之間表現(xiàn)為顯著正相關(guān)關(guān)系（P<0.05），相關(guān)性系數(shù)為0.366；氣溫與植被NPP相關(guān)性不顯著，相關(guān)系數(shù)為-0.223。結(jié)合圖6和圖7可以看出，降水與植被NPP之間多為體現(xiàn)為正相關(guān)性，其中17.01%的區(qū)域表現(xiàn)為顯著或者極顯著正相關(guān)，集中在山西省北部；氣溫與植被NPP之間多體現(xiàn)為負(fù)相關(guān)性，其中3.66%的區(qū)域表現(xiàn)為顯著或者極顯著負(fù)相關(guān)，集中在山西省中部；基于不同類型植被NPP與氣候因素之間偏相關(guān)性系數(shù)顯著性檢驗(yàn)結(jié)果圖（圖8），結(jié)合植被生態(tài)特征分析可知；不同類型植被NPP與降水、氣溫之間偏相關(guān)性的顯著性檢驗(yàn)結(jié)果和山西省整體檢驗(yàn)結(jié)果大致相同，其中草地生態(tài)系統(tǒng)結(jié)構(gòu)簡(jiǎn)單多變，對(duì)降水、溫度，尤其對(duì)降水的變化響應(yīng)最為顯著，林地生態(tài)系統(tǒng)結(jié)構(gòu)穩(wěn)定，對(duì)降水、溫度的變化響應(yīng)最不明顯。結(jié)合植被NPP變化研究結(jié)果能夠看出，部分植被NPP變化顯著區(qū)域，其植被NPP與降水、氣溫之間多表現(xiàn)顯著相關(guān)，表明氣候因素是造成植被NPP顯著變化的因素之一。

上述研究發(fā)現(xiàn)，降水多數(shù)促進(jìn)植被NPP增大，而氣溫多數(shù)抑制植被NPP增大，這正是干旱脅迫現(xiàn)象的體現(xiàn)，表明干旱是抑制山西省植被NPP增加的重要因素之一（史海萍，2006）。但氣溫和降水僅在整體水平上顯著影響植被NPP變化，通過逐象元研究分析發(fā)現(xiàn)，其對(duì)植被NPP顯著影響的區(qū)域面積比例很低，表明氣溫和降水僅是影響山西省植被NPP變化的多種因素中的兩種因素。

3.4 高程梯度上植被NPP變化特征

山西省境內(nèi)山脈連綿起伏、溝壑縱橫交錯(cuò)，全區(qū)高程梯度大，不同高程地區(qū)的環(huán)境及植被類型差異顯著?；诟叱虜?shù)據(jù)對(duì)山西省2005—2015年植被NPP數(shù)據(jù)進(jìn)行分區(qū)統(tǒng)計(jì)。從圖9可以看出，隨著高程增加，年均植被NPP呈現(xiàn)波動(dòng)增大，不同類型植被NPP隨高程變化趨勢(shì)和研究區(qū)域總體趨勢(shì)大致相同，在高程介于700～800 m之間區(qū)域，研究區(qū)域整體植被NPP以及不同類型植被NPP都有一個(gè)明顯減小的趨勢(shì)，結(jié)土地利用類型數(shù)據(jù)分析，山西省中部城鎮(zhèn)地區(qū)集中在高程700～800 m處，人類生產(chǎn)活動(dòng)強(qiáng)烈，抑制了植被改善。

4 討論與結(jié)論

本研究結(jié)果表明，山西省植被NPP的空間分布變化主要取決于空間位置、地形地貌、開發(fā)程度和人為管理等因素，這些因素通過控制水熱、土壤等因子影響著植被NPP的空間分布以及變化特征， 其中干旱是抑制植被NPP增加的重要因素之一。受溫帶大陸季風(fēng)性氣候影響，山西省雨季和熱季同期，但降水量遠(yuǎn)小于蒸發(fā)量，容易形成干旱脅迫現(xiàn)象，抑制植被NPP的增加；草地植被以及耕地植被受人類影響較強(qiáng)，NPP變化明顯，林地多位于山區(qū)，受人類影響較小，NPP數(shù)值較為穩(wěn)定。

本研究中，山西省植被整體上呈好轉(zhuǎn)趨勢(shì)，與國(guó)內(nèi)許多專家學(xué)者的研究結(jié)果相同，黨躍軍等（2015）的研究發(fā)現(xiàn)山西省在2003—2012年10 a間植被呈好轉(zhuǎn)趨勢(shì)，草地、灌叢、林地均有改善。武永利等（2008）對(duì)山西省植被指數(shù)時(shí)空變化特征的研究發(fā)現(xiàn)山西省近25 a間植被指數(shù)波動(dòng)上升，植被總體呈好轉(zhuǎn)趨勢(shì)。杜加強(qiáng)等（2015）關(guān)于新疆地區(qū)植被NDVI的研究得出植被從1982—2012年逐漸改善；李輝霞等（2011）對(duì)三江源區(qū)域NDVI的研究得出植被從2000—2014年逐漸變好；李曉榮等（2017）對(duì)太行山植被NPP的研究中，發(fā)現(xiàn)植被NPP在2000—2014年間呈好轉(zhuǎn)趨勢(shì)，但有些區(qū)域在研究期間植被并未得到改善。張莎（2015）基于長(zhǎng)時(shí)間多系列遙感數(shù)據(jù)對(duì)河北植被NPP的研究得出河北省大部分區(qū)域植被NPP在2000—2014年間逐漸減小。為此阿多等（2017）提出，不同地區(qū)植被的變化規(guī)律之間存在明顯的差異性，所以，不同區(qū)域的植被在面臨相同的環(huán)境改變條件下，體現(xiàn)出的響應(yīng)特征也會(huì)表現(xiàn)出較明顯的差異。本研究目前僅對(duì)MODIS NPP數(shù)據(jù)進(jìn)行了分析應(yīng)用，缺少實(shí)地的考察認(rèn)識(shí)，結(jié)合社會(huì)、人為、經(jīng)濟(jì)方面的研究還不夠。對(duì)人類活動(dòng)造成植被NPP變化的研究缺乏定量研究，實(shí)際的應(yīng)用意義不大，這是在以后研究中需要改進(jìn)完善的。MODIS數(shù)據(jù)的精度還存在一定的問題，需要綜合考慮多方面的影響因素，并且結(jié)合研究區(qū)域的具體情況，這樣才能在本土化研究中取得更準(zhǔn)確的研究結(jié)果。目前，不同植被NPP模型之間在精度和準(zhǔn)度方面存在較大差異，其研究結(jié)果缺乏可比性；有學(xué)者認(rèn)為綜合多種算法的集合預(yù)估方法能夠有效地提高模型模擬精度，其模擬結(jié)果明顯好于單一算法，這將會(huì)是成為今后關(guān)于植被NPP研究的一個(gè)重要前進(jìn)方向。

綜上所述，得出結(jié)論如下：

（1）山西省植被NPP 2005—2015年間平均值為326.5 g（C）·m-2·a-1，其中93.53%的區(qū)域植被NPP值在200～500 g（C）·m-2·a-1之間，草地、耕地、灌叢和林地的NPP平均值分別為300.3、353.6、366.5、390.1 g（C）·m-2·a-1。（2）2005—2015年，山西省植被NPP整體波動(dòng)增大、變化顯著，顯著增大及極顯著增大區(qū)域面積比例達(dá)56.33%，主要集中在山西省西部；植被NPP顯著減小以及極顯著減小區(qū)域主要分布于山西省東南角，占比僅為2.22%。不同類型植被NPP的年際變化趨勢(shì)與研究區(qū)域植被NPP總體變化趨勢(shì)大致相同。（3）草地NPP變化速率最大，耕地大于灌叢，而林地NPP變化速率最小。林地生態(tài)結(jié)構(gòu)穩(wěn)定，NPP變化穩(wěn)定；草地生態(tài)系統(tǒng)結(jié)構(gòu)簡(jiǎn)單，適應(yīng)能力較強(qiáng)，對(duì)環(huán)境變化響應(yīng)更迅速，NPP變化明顯。（4）山西省植被NPP與氣溫因素在整體水平上相關(guān)性不顯著，在象元水平上表現(xiàn)為負(fù)相關(guān)性，其中3.66%的區(qū)域表現(xiàn)為顯著或者極顯著負(fù)相關(guān)，主要集中在山西省中部；山西省植被NPP與降水因素在整體水平上表現(xiàn)為顯著正相關(guān)，在象元水平上表現(xiàn)為正相關(guān)性，其中17.01%的區(qū)域表現(xiàn)為顯著或者極顯著正相關(guān)，主要集中在山西省北部。（5）研究期間，氣候因素整體上促進(jìn)了山西省植被NPP上升，而人為因素主要抑制植被NPP上升。

參考文獻(xiàn)：

A D， ZHAO WJ， GONG ZN， et al， 2017. Temporal analysis of climate change and its relationship with vegetation cover on the North China plain from 1981 to 2013 [J]. Acta Ecol Sin， 37（2）： 576-592. [阿多， 趙文吉， 宮兆寧， 等， 2017. 1981-2013華北平原氣候時(shí)空變化及其對(duì)植被覆蓋度的影響 [J]. 生態(tài)學(xué)報(bào)， 37（2）：576-592.]

ANAYETI A， SHI QD， LIU M， et al， 2014. The vegetation cover changes of regress analysis on using time-serial images of remote sensing [J]. J Xinjiang Univ （Nat Sci Ed）， 34（3）： 341-344. [阿斯姆古麗·阿納耶提， 師慶東， 劉曼， 等， 2014. 長(zhǎng)時(shí)間序列遙感數(shù)據(jù)回歸分析支持下的植被變化分析 [J]. 新疆大學(xué)學(xué)報(bào)（自然科學(xué)版）， 31（3）：341-344.]

CHEN LJ， LIU GH， LI HG， 2002. Estimating net primary productivity of terrestrial vegetation in China using remote sensing [J]. J Remote Sens， 6（2）： 129-135. [陳利軍， 劉高煥， 勵(lì)惠國(guó)， 2002. 中國(guó)植被凈第一性生產(chǎn)力遙感動(dòng)態(tài)監(jiān)測(cè) [J]. 遙感學(xué)報(bào)， 6（2）：129-135.]

DU JQ， AHATI J， ZHAO CX， et al， 2015. Dynamic changes in vegetation NDVI from 1982 to 2012 and its response to climate change and human activities in Xinjiang， China [J]. Chin J Appl Ecol， 26（12）： 3567-3578. [杜加強(qiáng)， 賈爾恒·阿哈提， 趙晨曦， 等， 2015. 1982-2012年新疆植被NDVI的動(dòng)態(tài)變化及其對(duì)氣候變化和人類活動(dòng)的響應(yīng) [J]. 應(yīng)用生態(tài)學(xué)報(bào)， 26（12）：3567-3578.]

DANG YJ， WANG LX， YAN JX， 2015. Spatio-temporal dynamic of vegetation and its response to climatic factors in Shanxi Province from 2003 to 2012 [J]. Res Soil Water Conserv， 22（2）： 235-240. [黨躍軍， 王禮霄， 嚴(yán)俊霞， 2015. 山西省2003-2012年植被時(shí)空變化格局及對(duì)氣候因子的響應(yīng) [J]. 水土保持研究， 22（2）：235-240.]

LIETH H， WHITTAKER RH， 1975. Primary Productivity of the Biosphere [M]. Berlin Heidelberg： Springer-Verlag.

LI HX， LIU GH， FU BJ， 2011. Response of vegetation to climate change and human activity based on NDVI in the Three-River Headwaters region [J]. Acta Ecol Sin， 31（19）： 5495-5504. [李輝霞， 劉國(guó)華， 傅伯杰， 2011. 基于NDVI的三江源地區(qū)植被生長(zhǎng)對(duì)氣候變化和人類活動(dòng)的響應(yīng)研究 [J]. 生態(tài)學(xué)報(bào)， 31（19）：5495-5504.]

LI XR， GAO H， HAN LP， et al， 2017. Spatio-temporal variations NPP and the driving factors in Taihang Mountain Area [J]. Chin J Ecol-Agric， 25（4）： 498-508. [李曉榮， 高會(huì)， 韓立樸， 等， 2017. 太行山植被NPP時(shí)空變化特征以及驅(qū)動(dòng)力分析 [J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)， 25（4）：498-508.]

MU SJ， LI JL， ZHOU W， et al， 2013. Spatial-temporal distribution of net primary productivity and its relationship with climate factors inInner Mongolia from 2001 to 2010 [J]. Acta Ecol Sin， 33（12）：3752-3764. [穆少杰， 李建龍， 周偉， 等， 2013. 2001-2010年內(nèi)蒙古植被凈初級(jí)生產(chǎn)力的時(shí)空格局及其與氣候的關(guān)系 [J]. 生態(tài)學(xué)報(bào)， 33（12）：3752-3764.]

PIAO SL， FANG JY， GUO QY， 2001. Application of CASA model to the estimation on of Chinese terrestrial net primary productivity [J]. Chin J Plant Ecol， 25（5）： 603-608. [樸世龍， 方精云， 郭慶華， 2001. 利用CASA模型估算我國(guó)植被凈第一性生產(chǎn)力 [J]. 植物生態(tài)學(xué)報(bào)， 25（5）：603-608.]

SHI HP， 2006. Research of the climate of Shanxi Province and the impact on the vegetation ecosystems [D]. Lanzhou： Lanzhou University： 34-56. [史海萍， 2006. 山西氣候變化及植被生態(tài)系統(tǒng)的響應(yīng)研究 [D]. 蘭州： 蘭州大學(xué)：34-56.]

SUN LW， 2013. Temporal and spatial variation of vegetable net primary productivity and influence of climate change and human activity in Qilian mountains [D].Lanzhou： Northwest Normal University. [孫力煒， 2013. 祁連山區(qū)植被凈第一性生產(chǎn)力的時(shí)空分布特征及氣候變化和人類活動(dòng)的影響 [D]. 蘭州： 西北師范大學(xué).]

SUN R， ZHU QJ， 2001. Effect of climate of terrestrial net primary productivity in China [J]. J Remote Sens， 5（1）： 58-61. [孫睿， 朱啟疆， 2001. 氣候變化對(duì)中國(guó)陸地植被凈第一性生產(chǎn)力影響的初步研究 [J]. 遙感學(xué)報(bào)， 5（1）：58-61.]

WU SS， YAO ZJ， JIANG LG， et al， 2016. The spatial-temporal variations and hydrological effects of vegetation NPP based on MODIS in the source of the Yangtze River [J]. J Nat Resourc， 31（1）： 39-51. [吳珊珊， 姚治君， 姜麗光， 等， 2016. 基于MODIS的長(zhǎng)江源植被NPP時(shí)空變化特征及其水文效應(yīng) [J]. 自然資源學(xué)報(bào)， 31（1）：39-51.]

WU YL， LUAN Q， ZHAO YQ， et al， 2008. Analyses on spatial-temporal change characteristics of NDVI during the past 25 years in Shanxi [J]. Ecol Environ， 17（6）： 2330-2335. [武永利， 欒青， 趙永強(qiáng)， 等， 2008. 近25年山西植被指數(shù)時(shí)空變化特征分析 [J]. 生態(tài)環(huán)境， 17（6）：2330-2335.]

YUAN W， CAI WW， LIU D， et al， 2014. Satellite-based vegetation production models of terrestrial ecosystem： An overview [J]. Advan Earth Sci， 29（5）： 541-550. [袁文平， 蔡文文， 劉丹， 等， 2014. 陸地生態(tài)系統(tǒng)植被生產(chǎn)力遙感模型研究進(jìn)展 [J]. 地球科學(xué)進(jìn)展， 29（5）：541-550.]

ZHANG S， 2015. Research on spatio-temporal distribution of vegetation net primary productivity in Hebei Province based on longtime-series multi-source remote sensing data [D]. Shijiazhuang： Hebei Normal University. [張莎， 2015. 基于長(zhǎng)時(shí)間序列多源遙感數(shù)據(jù)的河北省植被 NPP時(shí)空分布研究 [D]. 石家莊：河北師范大學(xué).]

ZHU F， LIU ZM， WANG ZM， et al， 2010. Temporal-spatial characteristic and factors influencing crop NPP across northeastern China [J]. Resourc Sci， 32（11）： 2079-2084. [朱鋒， 劉志明， 王宗明， 等， 2010. 東北地區(qū)農(nóng)田凈初級(jí)生產(chǎn)力時(shí)空特征及其影響因素分析 [J]. 資源科學(xué)， 32（11）：2079-2084.]

ZHU WQ， PAN YZ， ZHANG JS， 2007. Estimation of net primary productivity of Chinese terrestrial vegetation based on remote sensing [J]. Chin J Plant Ecol， 31（3）： 413-424. [朱文泉， 潘耀忠， 張錦水， 2007. 中國(guó)陸地植被凈初級(jí)生產(chǎn)力遙感估算 [J]. 植物生態(tài)學(xué)報(bào)， 31（3）：413-424.]