王 靜，萬(wàn)紅蓮*，張 翀

(1.寶雞文理學(xué)院 陜西省災(zāi)害監(jiān)測(cè)與機(jī)理模擬重點(diǎn)實(shí)驗(yàn)室，陜西 寶雞 721013；2.寶雞文理學(xué)院 地理與環(huán)境學(xué)院，陜西 寶雞 721013)

植被和氣候變化是影響自然生態(tài)系統(tǒng)能量收支平衡的直接因素，也是人類生存環(huán)境的重要組成部分[1]。IPCC第四次工作報(bào)告指出，近百年來(lái)全球氣溫變化總趨勢(shì)發(fā)生了明顯變化[2]，我國(guó)氣候變化趨勢(shì)與全球趨勢(shì)基本一致。在這樣的背景下，掌握陸地植被覆蓋年際變化以及對(duì)氣候要素的響應(yīng)規(guī)律，對(duì)評(píng)價(jià)陸地生態(tài)系統(tǒng)的環(huán)境質(zhì)量、調(diào)節(jié)植被覆蓋與人類活動(dòng)的協(xié)調(diào)發(fā)展具有重要的理論和實(shí)際意義[3]。以MODIS(中分辨率成像光譜儀)為光學(xué)遙感儀器獲得的NDVI數(shù)據(jù)(即歸一化植被指數(shù)，Normalized Difference Vegetation Index)具有廉價(jià)性和高時(shí)間分辨率，已經(jīng)成為目前應(yīng)用最為廣泛、前景最為廣闊的遙感數(shù)據(jù)之一[4-5]。而NPP，即植物的凈初級(jí)生產(chǎn)力(net primary productivity)，指單位面積單位時(shí)間內(nèi)綠色植物呼吸后剩下的所累積的能量或所生產(chǎn)的有機(jī)物質(zhì)，是生態(tài)系統(tǒng)碳平衡評(píng)估的基礎(chǔ)，也是表征植被活動(dòng)和評(píng)價(jià)陸地生態(tài)系統(tǒng)穩(wěn)定性的關(guān)鍵參數(shù)[6-7]。我國(guó)人多地少且資源相對(duì)匱乏，對(duì)NDVI和NPP的研究就顯得更為迫切和重要。國(guó)內(nèi)外學(xué)者基于兩種數(shù)據(jù)對(duì)氣候因子的響應(yīng)規(guī)律做了大量研究，塞內(nèi)加爾草原凈初級(jí)生產(chǎn)力(NPP)或生物量和作物產(chǎn)量的長(zhǎng)期變化趨勢(shì)不僅受降水的影響，而且與土地覆蓋和土地利用變化具有顯著的相關(guān)性[8]；5～9月是中國(guó)東北樣帶NPP發(fā)展的集中時(shí)段，占年NPP的89.8%。夏季(6～8月)NPP的變化是NPP年際變化的主要影響因素，占年NPP的65.9%[9]；降水不僅是內(nèi)蒙古草原NPP的重要?dú)夂蛴绊懸蜃?，而且還具有明顯的累積效應(yīng)[10]。相對(duì)于氣溫和降水量因子，人類活動(dòng)是黃土高原植被NPP變化的主要影響因素，且對(duì)NPP有雙重影響[11]。

寶雞地區(qū)地勢(shì)復(fù)雜，植被類型多樣，分縣(市)區(qū)(以下簡(jiǎn)稱縣區(qū))詳細(xì)研究該地區(qū)植被覆蓋變化與氣候響應(yīng)對(duì)渭河沿岸的農(nóng)田建設(shè)、秦嶺以及黃土高原地區(qū)植被恢復(fù)具有重要的現(xiàn)實(shí)指導(dǎo)意義。對(duì)于該區(qū)域的相關(guān)研究，眾多學(xué)者只是將其作為秦嶺或西北地區(qū)以及黃土高原的一部分來(lái)進(jìn)行植被NDVI的相關(guān)分析[12-14]，且秦超[15]、翟雅倩[16]整體研究了近年來(lái)寶雞市的植被覆蓋度動(dòng)態(tài)變化特征，發(fā)現(xiàn)寶雞市整體植被覆蓋度高，以中覆蓋度、中高覆蓋度和高覆蓋度為主，各行政區(qū)域植被覆蓋度差異性顯著。但這兩位學(xué)者對(duì)局部自然環(huán)境特征的表征沒(méi)有細(xì)化。而就寶雞地區(qū)植被覆蓋內(nèi)部差異及基于NPP數(shù)據(jù)的相關(guān)性分析的研究鮮見(jiàn)。

本研究基于MODIS NDVI和NPP數(shù)據(jù)，結(jié)合該地區(qū)11個(gè)氣象站點(diǎn)(其中，寶雞縣2003年撤縣設(shè)立陳倉(cāng)區(qū)，本文仍沿用舊制寶雞縣)氣溫和降水量的實(shí)測(cè)數(shù)據(jù)，分析了年際尺度上區(qū)內(nèi)植被覆蓋與氣候因子的相關(guān)性及其對(duì)氣候因子的響應(yīng)，以期對(duì)區(qū)域農(nóng)業(yè)發(fā)展提供理論參考。

1 研究區(qū)概況

寶雞位于陜西省(關(guān)中平原)西部，地處東經(jīng)106°18′～108°03′，北緯33°35′～35°06′。東與咸陽(yáng)市相連，南與漢中市相接，西、西北分別與甘肅省天水市和平?jīng)鍪信?，寶雞地區(qū)包括11個(gè)氣象站點(diǎn)(圖1)，總面積1.8萬(wàn)km2。寶雞屬暖溫帶半濕潤(rùn)氣候，受東亞季風(fēng)(包括高原季風(fēng))控制，造就了冬季寒冷干燥，夏季溫?zé)岫嘤昱c炎熱干燥交替，春季多變少雨且升溫迅速，秋季多陰雨連綿且降溫迅速的全年氣候變化特征。氣溫、降水量季相較為分明。年平均氣溫13 ℃左右，年降水量710～1000 mm，7、8、9這3個(gè)月降水集中，占全年降水的50%左右；寶雞地區(qū)是關(guān)中降水量最多的地區(qū)。寶雞地區(qū)生境條件和植被類型多樣，主要包括落葉闊葉林、針葉林、草原、灌叢和栽培植被等。森林覆被率達(dá)36%～42%，但分布不均，以秦嶺和關(guān)山為集中分布區(qū)域，關(guān)中和臺(tái)塬地區(qū)主要為栽培植被。

圖1 寶雞地區(qū)概況及氣象站點(diǎn)分布

2 資料與方法

2.1 數(shù)據(jù)來(lái)源

表1為本研究主要數(shù)據(jù)源信息，包括氣象數(shù)據(jù)(11個(gè)氣象站點(diǎn))、MODIS NDVI、NPP和高程DEM數(shù)據(jù)。研究中站點(diǎn)選取分別是：寶雞市、寶雞縣、鳳翔縣、岐山縣、眉縣、扶風(fēng)縣、隴縣、千陽(yáng)縣、麟游縣、鳳縣和太白縣，本研究也采用此11個(gè)站點(diǎn)作為行政區(qū)域進(jìn)行下文研究。

2.2 研究方法

2.2.1 趨勢(shì)分析 通過(guò)模擬每個(gè)柵格的變化趨勢(shì)和單個(gè)像元時(shí)間變化特征來(lái)反映整個(gè)空間變化規(guī)律，進(jìn)而表征區(qū)域時(shí)空格局演變[17]。Slope為NDVI在各個(gè)像元點(diǎn)的變化趨勢(shì)線斜率。若Slope值為正，則區(qū)域植被NDVI隨時(shí)間變化呈上升趨勢(shì)，且數(shù)值大小與植被覆蓋度增加趨勢(shì)強(qiáng)弱呈正比例關(guān)系；若Slope值為負(fù)，則表示植被指數(shù)隨時(shí)間變化呈下降趨勢(shì)。其計(jì)算公式如下：

(1)

公式(1)中：n為年份，i為年序列，NDVIi為第i年的最大化NDVI值。

表1 研究所用數(shù)據(jù)源

2.2.2 相關(guān)分析法 運(yùn)用相關(guān)分析法分析植被覆蓋變化對(duì)氣候因子的響應(yīng)，結(jié)合ArcGIS10.1和Matlab軟件，計(jì)算NDVI與氣溫和降水量間的相關(guān)系數(shù)，采用NDVI影像與經(jīng)反距離權(quán)重插值法獲取年均溫和年降水量的柵格數(shù)據(jù)，繪制NDVI與年均溫和年降水的相關(guān)系數(shù)空間分布圖，NDVI與氣候因子的空間關(guān)系通過(guò)顯著性檢驗(yàn)來(lái)進(jìn)行分析。計(jì)算公式如下[18]：

(2)

3 結(jié)果與分析

3.1 寶雞地區(qū)植被覆蓋變化分析

經(jīng)ArcGIS區(qū)域統(tǒng)計(jì)功能獲取年均NDVI，表示寶雞地區(qū)當(dāng)年植被覆蓋水平(圖2)。2001～2013年NDVI變化情況為：2006～2007、2008～2010、2011～2012年緩慢增長(zhǎng)，2001～2002、2003～2004年飛速增長(zhǎng)，2004～2006、2007～2008、2010～2011年緩慢下降，2002～2003年迅速下降，2012～2013年變化較不明顯?？傮w來(lái)看，NDVI以0.032/10 a的速度呈上升趨勢(shì)。三北防護(hù)林工程區(qū)1982～2006年近25 a植被覆蓋的平均增速為0.007/10 a，與此相比，寶雞地區(qū)植被覆蓋增速較快[19]，且NDVI值在2010年達(dá)到最高值0.572。2001年寶雞地區(qū)平均NDVI值為0.5274，2013年增加至0.5668，增速為7.47%。說(shuō)明寶雞地區(qū)近13 a來(lái)的植被覆蓋有增無(wú)減。

3.2 植被覆蓋的空間分布特征

寶雞地區(qū)2001～2013年多年平均NDVI空間分布如圖3所示，整體上呈從東北向西南遞增的特征。由于太白山的山地地貌類型和其為秦嶺山脈主峰的原因，以及太白河和水河的河水補(bǔ)給，寶雞地區(qū)植被覆蓋最高值集中分布在太白縣南部。千河沿岸以及渭河流域?qū)氹u段的寶雞市、寶雞縣、鳳翔、岐山南北沿岸的“帶狀”地區(qū)，植被覆蓋度最低，主要原因是這些地區(qū)人口集聚分布，經(jīng)濟(jì)發(fā)展最為迅速，隨著城市化的推動(dòng)，城市擴(kuò)建對(duì)于土地需求逐日擴(kuò)大，以0.49%/a的速度增加，建筑面積的迅速增加極大地削減了植被面積。此外，大量的農(nóng)田建設(shè)和農(nóng)作物種植呈現(xiàn)出植被覆蓋明顯的季節(jié)分布特征，從1999年陜西地區(qū)退耕還林試點(diǎn)實(shí)施以來(lái)，寶雞地區(qū)的農(nóng)作物種植、綠化(包括城市綠化和公園綠化)及林地面積分別以0.11%/a的速度減少和以1.07%/a、0.21%/a的速度增加，而研究區(qū)植被覆蓋面積主要依靠于農(nóng)作物種植，農(nóng)作物播種面積占總植被覆蓋面積的34.44%，是林業(yè)面積的一半，因而該區(qū)域植被覆蓋表現(xiàn)最低。

圖2 2001～2013年寶雞地區(qū)年均NDVI變化趨勢(shì)

3.3 植被覆蓋空間格局演變態(tài)勢(shì)

寶雞地區(qū)2001～2013年植被覆蓋變化具有明顯的區(qū)域差異(圖4)。由于近年來(lái)寶雞地區(qū)正在大力實(shí)施退耕還林還草和保護(hù)森林的政策[20]，植被恢復(fù)較好的區(qū)域廣泛分布于千河沿岸的隴縣、千陽(yáng)、寶雞縣地區(qū)，漆水河以北的麟游地區(qū)，嘉陵江沿岸的鳳縣西南地區(qū)，狀況逐年好轉(zhuǎn)。而近年來(lái)城市化水平逐年提升，城市土地?cái)U(kuò)張需求速度有增無(wú)減，所以大量的城區(qū)建筑用地占用了耕地及草地資源，使得千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸的植被覆蓋度降低。

為了更好地表征寶雞地區(qū)植被恢復(fù)狀況，參照已有研究[17,21-23]，將Slope分為如圖4、表2所示的7個(gè)等級(jí)。結(jié)果顯示，嚴(yán)重退化面積1 km2，中度退化面積117 km2，分別占總面積的0.006%和0.646%，兩個(gè)等級(jí)的植被覆蓋分布狀況集中在千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸“帶狀”地區(qū)；輕微退化面積772 km2，占總面積的4.262%，主要分布在千河以西的隴縣地區(qū)、渭河沿岸南部及太白河和水河之間的太白縣地區(qū)；基本不變面積4426 km2，輕度改善面積6126 km2，分別占總面積的24.434%和33.819%，零星分布于除千河以西的隴縣地區(qū)、渭河沿岸南部以外的寶雞其他地區(qū)；中度改善面積4488 km2，占總面積的24.776%，明顯改善面積2184 km2，占總面積的12.057%，兩者集中于千河沿岸的隴縣、千陽(yáng)、寶雞縣地區(qū)、漆水河以北的麟游地區(qū)、嘉陵江沿岸的鳳縣西南部地區(qū)。綜上所述，近13 a寶雞地區(qū)植被覆蓋狀況呈穩(wěn)中上升的趨勢(shì)，這主要是因?yàn)閷氹u地區(qū)實(shí)施“退耕還林還草”、森林保護(hù)等政策，使得該地綠化城市面積和造林面積分別以0.91%/a和0.73%/a的增速改善，以及當(dāng)?shù)厝藗儽Ｗo(hù)環(huán)境意識(shí)的不斷增強(qiáng)，植被覆蓋狀況逐年見(jiàn)好，由此生態(tài)環(huán)境也得到逐年改善[24]。

圖3 年均NDVI變化斜率趨勢(shì)圖

3.4 植被覆蓋與氣溫和降水的關(guān)系

NDVI數(shù)據(jù)的提取是建立在氣候要素對(duì)植被前期生長(zhǎng)作用的基礎(chǔ)之上的，因而研究NDVI和氣溫、降水量在植被生長(zhǎng)期內(nèi)的關(guān)系具有重要的實(shí)際意義[13-14,25-27]。由圖6和表3可知，近13 a來(lái)寶雞地區(qū)植被覆蓋與氣溫呈正相關(guān)的區(qū)域像元占總像元的80.99%，降水量呈正相關(guān)的區(qū)域像元占79.01%，表明年均NDVI與年均溫和年降水的相關(guān)性整體上表現(xiàn)均較好。

植被覆蓋與氣溫呈正相關(guān)的區(qū)域像元數(shù)占總像元的80.99%，除鳳縣大部分地區(qū)和千陽(yáng)中西部地區(qū)外，廣泛分布于寶雞地區(qū)各個(gè)縣區(qū)。由于以高山草甸為主體的關(guān)山草原位于隴縣西南部、太白縣中部至東部，秦嶺主峰穿其境而過(guò)且有太白河和水河水分的補(bǔ)給，因此其相應(yīng)區(qū)域相關(guān)系數(shù)最大；千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸“帶狀”地區(qū)，地勢(shì)平坦、氣溫四季分明，土壤有機(jī)質(zhì)的分解在夏季氣溫升高的條件下加速，對(duì)常溫地區(qū)植被凈初級(jí)生產(chǎn)力的提高起到了有效的促進(jìn)作用[28-30]，又加之大量植被綠化面積的建設(shè)，所以其相關(guān)系數(shù)最大(表3)。植被覆蓋與氣溫呈負(fù)相關(guān)的區(qū)域集中分布在鳳縣大部分地區(qū)和千陽(yáng)中西部地區(qū)，原因可能是鳳縣地處秦嶺山地地區(qū)，且紫柏山坐落其境，另外有嘉陵江和褒河水分的補(bǔ)給。吳岳支脈箭筈嶺坐落于千陽(yáng)南部，兩個(gè)地區(qū)均海拔較高、年均溫較低。當(dāng)植物生長(zhǎng)的最適溫度低于夏季氣溫時(shí)，植株內(nèi)Ru Bis CO(核酮糖-1,5-二磷酸羧化酶/加氧酶)的羧化反應(yīng)小于其加氧反應(yīng)，加速了凈光合作用的下降速度[31]；同時(shí)，土壤水分蒸發(fā)和植物蒸騰作用會(huì)因溫度的升高而加劇，嚴(yán)重影響了植被生長(zhǎng)[32]?？傮w而言，介于0.3～0.9相關(guān)系數(shù)的像元均勻分布在寶雞各個(gè)地區(qū)，說(shuō)明寶雞地區(qū)植被覆蓋與氣溫的相關(guān)性表現(xiàn)較好。

圖4 年均NDVI變化斜率分級(jí)圖

植被覆蓋與降水量呈負(fù)相關(guān)的區(qū)域集中在渭河流域?qū)氹u段的千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸“帶狀”地區(qū)(表3)，這些地區(qū)具有較強(qiáng)的降水量，而隨著城市化進(jìn)程的加快，城市擴(kuò)展帶來(lái)建筑面積擴(kuò)大，植被覆蓋度不斷降低，因此植被覆蓋與降水呈負(fù)相關(guān)；除此之外，寶雞其他各個(gè)縣區(qū)的植被覆蓋與降水量均呈正相關(guān)，其中鳳縣大面積地區(qū)呈非常顯著正相關(guān)。這是因?yàn)樗侵参锷L(zhǎng)過(guò)程中的決定性因素，植物的生長(zhǎng)狀況會(huì)直接受到降水量的影響。降水量是影響我國(guó)干旱變化最主要的因素[33]，降水量對(duì)植物生長(zhǎng)的限制程度尤以年降水量少的干旱、半干旱地區(qū)表現(xiàn)最為明顯[10]。綜上可見(jiàn)，寶雞地區(qū)植被覆蓋與降水量的相關(guān)性表現(xiàn)較好。

表2 2001～2013年寶雞地區(qū)年均NDVI變化趨勢(shì)分級(jí)統(tǒng)計(jì)

表3 2001～2013年寶雞地區(qū)年均NDVI與氣溫和降水量相關(guān)性的像元比例

3.5 植被恢復(fù)與NPP響應(yīng)關(guān)系分析

寶雞地區(qū)年均NPP為668.0831 g·C/(m2·a)，受降水量影響(年平均降水量529.89 mm)，2001、2002、2004年千陽(yáng)的NPP值均呈現(xiàn)最低，鳳翔在2003、2005～2010共7年的NPP均為最低。而11個(gè)行政區(qū)域2001～2010年NPP的最高值均出現(xiàn)在鳳縣。表4數(shù)據(jù)顯示，2001～2010年對(duì)年均NPP影響較大的因子為NDVI。在12個(gè)分析結(jié)果中，以麟游的顯著性最高，達(dá)到0.01顯著水平，岐山、隴縣、千陽(yáng)、鳳縣次之，達(dá)到0.05顯著水平，且均為正相關(guān)。依據(jù)生態(tài)學(xué)相關(guān)理論，NDVI與NPP具有直接或間接的相互關(guān)系，在區(qū)域和生態(tài)系統(tǒng)模型構(gòu)建中，NDVI常被直接或間接地用于計(jì)算植被NPP[34-35]。因此，NDVI與NPP具有同步增減的變化趨勢(shì)，即隨著植被覆蓋率的增加(減少)而NPP相應(yīng)上升(下降)。

氣溫和降水量與NPP的相關(guān)性均較差，均未達(dá)到顯著性水平，氣溫除鳳翔、麟游地區(qū)外其他地區(qū)均呈負(fù)相關(guān)，降水除岐山外其余地區(qū)均呈正相關(guān)。相關(guān)系數(shù)的顯著性均不如前50 a的結(jié)果[36]顯著，這一數(shù)據(jù)分析結(jié)果也在個(gè)別站點(diǎn)變化的差異較大。2000年以來(lái)國(guó)家推行天然林保護(hù)政策，寶雞地區(qū)特別是渭河沿岸及支流周邊區(qū)域陸續(xù)實(shí)行了包括封山育林、加強(qiáng)水源涵養(yǎng)及天然次生林保護(hù)、生態(tài)修復(fù)等退耕還林還草工程，寶雞地區(qū)在陜西省退耕還林試行初期(1999年)森林覆蓋率達(dá)34.30%，之后以1.43%/5 a的速度增長(zhǎng)，森林單位蓄積量也以385 m3/5 a的速度增加。與此同時(shí)，氣溫和降水等氣象因子也在逐年發(fā)生變化，在如此復(fù)雜的人類活動(dòng)和氣候變化的綜合作用下，本研究選取渭河流域中上游的寶雞地區(qū)，研究范圍較小，更具典型代表意義，對(duì)進(jìn)一步理清植被恢復(fù)與NPP變化的響應(yīng)關(guān)系具有良好的科學(xué)和實(shí)際意義。

表4 2001～2010年不同行政區(qū)域NPP與NDVI及氣候因子間相關(guān)系數(shù)統(tǒng)計(jì)

注：*表示在0.05水平上顯著相關(guān)；**表示在0.01水平上顯著相關(guān)。

4 討論

寶雞地區(qū)相比于干旱地區(qū)，水熱條件適中，由于水資源有一定保障，因而對(duì)于植被生長(zhǎng)而言溫度更易成為其主要驅(qū)動(dòng)因素[37-38]。理論上而言，日積溫伴隨著溫度的上升而增加，對(duì)大多數(shù)植被生長(zhǎng)顯然有積極作用，鳳縣大部分地區(qū)(秦嶺、紫柏山、嘉陵江、褒河)和千陽(yáng)中西部地區(qū)(吳岳支脈箭筈嶺)植被覆蓋程度大幅提高，生態(tài)恢復(fù)效果最為明顯，但植被與氣溫卻呈現(xiàn)負(fù)相關(guān)關(guān)系。可見(jiàn)，國(guó)家退耕還林政策的實(shí)施使這些山地地區(qū)林草用地增加，單株植被的長(zhǎng)勢(shì)與整體植被覆蓋程度的關(guān)系不大，因此，溫度的驅(qū)動(dòng)作用很難體現(xiàn)出來(lái)。同時(shí)，寶雞具有以山地、丘陵為主的復(fù)雜地形，在海拔較高地區(qū)淺層地下水儲(chǔ)存不足，天然降水是植被生長(zhǎng)需水依賴的主要來(lái)源，從而形成了寶雞地區(qū)植被與降水量的相關(guān)性中東部低于北、西、南部的格局；而在千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸“帶狀”地區(qū)，由于地勢(shì)平坦，人類活動(dòng)在城鎮(zhèn)化進(jìn)程中強(qiáng)度增大，植被生長(zhǎng)受其間接影響與降水量呈負(fù)相關(guān)發(fā)展趨勢(shì)。因此，植被對(duì)溫度和降水的負(fù)相關(guān)并不代表需求為負(fù)，而應(yīng)考慮退耕還林、城市建設(shè)、土地利用等人類活動(dòng)的客觀因素影響。

5 結(jié)論

“退耕還林還草”、森林保護(hù)等政策實(shí)施以來(lái)，寶雞地區(qū)綠化城市面積增速(0.91%/a)和造林面積增速(0.73%/a)較好，表現(xiàn)在千河沿岸的隴縣、千陽(yáng)縣、寶雞縣沿岸地區(qū)，基本不變面積為4426 km2，占總面積的24.434%；明顯改善2184 km2，占12.057%，生態(tài)環(huán)境逐年改善明顯。

寶雞全區(qū)植被覆蓋與氣溫較顯著正相關(guān)像元(0.6～0.9)占15.016%，與降水量較顯著正相關(guān)像元占16.872%。鳳縣大部分地區(qū)和千陽(yáng)中西部地區(qū)植被覆蓋與年均溫呈現(xiàn)負(fù)相關(guān)，但鳳縣大部分地區(qū)與年降水量呈正相關(guān)，且正相關(guān)系數(shù)顯著；千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風(fēng)南北沿岸“帶狀”地區(qū)，植被覆蓋與年降水量呈負(fù)相關(guān)，但此地帶植被覆蓋與氣溫呈顯著正相關(guān)，以千河與渭河交叉口的寶雞縣東北部與鳳翔縣東南部交接地區(qū)顯著性最為明顯。

退耕后10 a NDVI與NPP呈良好的正相關(guān)性，以麟游縣最為顯著(0.01水平上顯著)，即凈初級(jí)生產(chǎn)力與植被覆蓋具有同步增減的變化趨勢(shì)。

[1] 方精云,宋永昌,劉鴻雁,等.植被氣候關(guān)系與我國(guó)的植被分區(qū)[J].植物學(xué)報(bào),2002,44(9):1105-1122.

[2] 趙宗慈,王紹武,羅勇.IPCC成立以來(lái)對(duì)溫度升高的評(píng)估與預(yù)估[J].氣候變化研究進(jìn)展,2007(3):183-184.

[3] Jacqueline G F, Antoine G, Niklaus E Z. Tree line shifts in the Swiss Alps: climate change or land abandonment?[J]. Journal of Vegetation Science, 2007, 18: 571-582.

[4] 馬志勇,張慧,張海軍.基于MODIS杭州地區(qū)植被變化分析[J].安徽農(nóng)業(yè)科學(xué),2006,34(18):4618-4619.

[5] 韓愛(ài)惠.用MODIS植被指數(shù)的時(shí)間序列分析提取荒漠化敏感區(qū)域的方法[J].林業(yè)資源管理,2004(1):57-60.

[6] 朱文泉,陳云浩,徐丹,等.陸地植被凈初級(jí)生產(chǎn)力計(jì)算模型研究進(jìn)展[J].生態(tài)學(xué)雜志,2005,24(3):296-300.

[7] Ruimy A, Saugier B, Dedieu G. Methodology for the estimation of terrestrial net primary production from remotely sensed data[J]. J Geophys Res, 1994(D): 5263-5283.

[8] Tottrup C, Rasmussen M S. Mapping long-term changes in savannah crop productivity in Senegal through trend analysis of time series of remote sensing data [J]. Agriculture, Ecosystems & Environments, 2004, 103(3): 545-560.

[9] Zhu W Q, Pan Y Z, Liu X, et al. Spatio-temporal distribution of net primary productivity along the northeast China transect and its response to climatic change[J]. Journal of Forestry Research, 2006,17(2): 93-98.

[10] 龍慧靈,李曉兵,黃玲梅,等.內(nèi)蒙古草原生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力及其與氣候的關(guān)系[J].植物生態(tài)學(xué)報(bào),2010,34(7):781-791.

[11] 周夏飛,馬國(guó)霞,曹國(guó)志,等.2001～2013年黃土高原植被凈初級(jí)生產(chǎn)力時(shí)空變化及其歸因[J].安徽農(nóng)業(yè)科學(xué),2017,45(14):48-53.

[12] 楊尚武,張勃.基于SPOT NDVI的甘肅河?xùn)|植被覆蓋變化及其對(duì)氣候因子的響應(yīng)[J].生態(tài)學(xué)雜志,2014,33(2):455-461.

[13] 崔曉臨,白紅英,王濤.秦嶺地區(qū)植被NDVI海拔梯度差異及其氣溫響應(yīng)[J].資源科學(xué),2013,35(3):618-626.

[14] 戴聲佩,張勃,王海軍,等.中國(guó)西北地區(qū)植被時(shí)空演變特征及其對(duì)氣候變化的響應(yīng)[J].遙感技術(shù)與應(yīng)用,2010,25(1):69-76.

[15] 秦超,李君軼,陳宏飛.近20年寶雞市植被覆蓋度動(dòng)態(tài)變化及驅(qū)動(dòng)力分析[J].山東農(nóng)業(yè)科學(xué),2014,46(9):98-105.

[16] 翟雅倩,周旗,張翀,等.近15年寶雞市植被覆蓋度變化研究[J].寶雞文理學(xué)院學(xué)報(bào):自然科學(xué)版,2016,36(4):31-36.

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

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

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

[20] 周德成,趙淑清,朱超.退耕還林工程對(duì)黃土高原土地利用/覆被變化的影響:以陜西省安塞縣為例[J].自然資源學(xué)報(bào),2011,26(11):1866-1878.

[21] 王震,閆文德,劉曙光,等.基于MODIS數(shù)據(jù)的中國(guó)三種主要土地類型變化的空間特征分析[J].生態(tài)學(xué)報(bào),2017,37(10):1-7.

[22] 張君,延軍平.1982～2013年陜西不同植被類型NDVI變化特征分析[J].干旱區(qū)資源與環(huán)境,2017,31(4):86-92.

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

[24] 宋德明.陜西省地理志[M].西安:陜西人民出版社,2000:221- 262.

[25] 金麗娜.西安近10年植被覆蓋變化及其與氣候因子的相關(guān)性研究[J].甘肅科學(xué)學(xué)報(bào),2014,26(3):23-26.

[26] 蒲蕾,任志遠(yuǎn).陜西省不同地區(qū)NDVI變化與氣候因子的關(guān)系及響應(yīng)研究[J].水土保持通報(bào),2013,32(2):265-275.

[27] 方精云,樸世龍,賀金生,等.近20年中國(guó)植被活動(dòng)在增強(qiáng)[J].中國(guó)科學(xué)C輯,2003,33(6):554-565.

[28] Gamer W. Comparing global models of terrestrial net primary productivity (NPP): Overview and key results[J]. Global Change Biology, 1995, 5(1): 1-15.

[29] Myneni R B, Keeling C D, Tucker C J, et al. Increased plant growth in the northern high latitudes from 1981 to 1991[J]. Nature, 1997, 386(6626): 698-702.

[30] Zhang X Y, Hu Y F, Zhuang D F, et al. The spatial pattern and differentiation of NDVI in Mongolia Plateau[J]. Geographical Research, 2009, 28(1): 10-18.

[31] 白文龍,張福平,倪海燕,等.關(guān)中地區(qū)植被覆蓋變化及其對(duì)氣候因子的響應(yīng)研究[J].農(nóng)業(yè)現(xiàn)代化研究,2013,34(1):104-108.

[32] 賈東于,文軍,張?zhí)锰?等.黃土高原降水對(duì)土壤含水量和導(dǎo)熱率的影響[J].高原氣象,2014,33(3):712-720.

[33] 翟盤(pán)茂,鄒旭愷.1951～2003年中國(guó)氣溫和降水變化及其對(duì)干早的影響[J].氣象變化研究進(jìn)展,2005,1(1):16-18.

[34] 周珺.基于遙感數(shù)據(jù)的重慶市凈初級(jí)生產(chǎn)力(NPP)時(shí)空特征研究[D].重慶:西南大學(xué),2013.

[35] 朱文泉,潘耀忠,何浩,等.中國(guó)典型植被最大光利用率模擬[J].科學(xué)通報(bào),2006,51(6):700-710.

[36] 蔣沖,王飛,穆興民.渭河流域相對(duì)濕度與植被覆蓋變化及其響應(yīng)關(guān)系[J].水土保持通報(bào),2012,32(5):122-127.

[37] 崔林麗,史軍,楊引明,等.中國(guó)東部植被NDVI對(duì)氣溫和降水的旬響應(yīng)特征[J].地理學(xué)報(bào),2009,64(7):850-860.

[38] 王靜,萬(wàn)紅蓮,張翀. 2001-2013年寶雞地區(qū)植被NDVI對(duì)氣溫和降水量的季節(jié)響應(yīng)特征[J].水土保持通報(bào),2017,37(5):235-240,349.