王 靜,萬紅蓮*,張 翀
(1.寶雞文理學院 陜西省災(zāi)害監(jiān)測與機理模擬重點實驗室,陜西 寶雞 721013;2.寶雞文理學院 地理與環(huán)境學院,陜西 寶雞 721013)
植被和氣候變化是影響自然生態(tài)系統(tǒng)能量收支平衡的直接因素,也是人類生存環(huán)境的重要組成部分[1]。IPCC第四次工作報告指出,近百年來全球氣溫變化總趨勢發(fā)生了明顯變化[2],我國氣候變化趨勢與全球趨勢基本一致。在這樣的背景下,掌握陸地植被覆蓋年際變化以及對氣候要素的響應(yīng)規(guī)律,對評價陸地生態(tài)系統(tǒng)的環(huán)境質(zhì)量、調(diào)節(jié)植被覆蓋與人類活動的協(xié)調(diào)發(fā)展具有重要的理論和實際意義[3]。以MODIS(中分辨率成像光譜儀)為光學遙感儀器獲得的NDVI數(shù)據(jù)(即歸一化植被指數(shù),Normalized Difference Vegetation Index)具有廉價性和高時間分辨率,已經(jīng)成為目前應(yīng)用最為廣泛、前景最為廣闊的遙感數(shù)據(jù)之一[4-5]。而NPP,即植物的凈初級生產(chǎn)力(net primary productivity),指單位面積單位時間內(nèi)綠色植物呼吸后剩下的所累積的能量或所生產(chǎn)的有機物質(zhì),是生態(tài)系統(tǒng)碳平衡評估的基礎(chǔ),也是表征植被活動和評價陸地生態(tài)系統(tǒng)穩(wěn)定性的關(guān)鍵參數(shù)[6-7]。我國人多地少且資源相對匱乏,對NDVI和NPP的研究就顯得更為迫切和重要。國內(nèi)外學者基于兩種數(shù)據(jù)對氣候因子的響應(yīng)規(guī)律做了大量研究,塞內(nèi)加爾草原凈初級生產(chǎn)力(NPP)或生物量和作物產(chǎn)量的長期變化趨勢不僅受降水的影響,而且與土地覆蓋和土地利用變化具有顯著的相關(guān)性[8];5~9月是中國東北樣帶NPP發(fā)展的集中時段,占年NPP的89.8%。夏季(6~8月)NPP的變化是NPP年際變化的主要影響因素,占年NPP的65.9%[9];降水不僅是內(nèi)蒙古草原NPP的重要氣候影響因子,而且還具有明顯的累積效應(yīng)[10]。相對于氣溫和降水量因子,人類活動是黃土高原植被NPP變化的主要影響因素,且對NPP有雙重影響[11]。
寶雞地區(qū)地勢復(fù)雜,植被類型多樣,分縣(市)區(qū)(以下簡稱縣區(qū))詳細研究該地區(qū)植被覆蓋變化與氣候響應(yīng)對渭河沿岸的農(nóng)田建設(shè)、秦嶺以及黃土高原地區(qū)植被恢復(fù)具有重要的現(xiàn)實指導(dǎo)意義。對于該區(qū)域的相關(guān)研究,眾多學者只是將其作為秦嶺或西北地區(qū)以及黃土高原的一部分來進行植被NDVI的相關(guān)分析[12-14],且秦超[15]、翟雅倩[16]整體研究了近年來寶雞市的植被覆蓋度動態(tài)變化特征,發(fā)現(xiàn)寶雞市整體植被覆蓋度高,以中覆蓋度、中高覆蓋度和高覆蓋度為主,各行政區(qū)域植被覆蓋度差異性顯著。但這兩位學者對局部自然環(huán)境特征的表征沒有細化。而就寶雞地區(qū)植被覆蓋內(nèi)部差異及基于NPP數(shù)據(jù)的相關(guān)性分析的研究鮮見。
本研究基于MODIS NDVI和NPP數(shù)據(jù),結(jié)合該地區(qū)11個氣象站點(其中,寶雞縣2003年撤縣設(shè)立陳倉區(qū),本文仍沿用舊制寶雞縣)氣溫和降水量的實測數(shù)據(jù),分析了年際尺度上區(qū)內(nèi)植被覆蓋與氣候因子的相關(guān)性及其對氣候因子的響應(yīng),以期對區(qū)域農(nóng)業(yè)發(fā)展提供理論參考。
寶雞位于陜西省(關(guān)中平原)西部,地處東經(jīng)106°18′~108°03′,北緯33°35′~35°06′。東與咸陽市相連,南與漢中市相接,西、西北分別與甘肅省天水市和平?jīng)鍪信?,寶雞地區(qū)包括11個氣象站點(圖1),總面積1.8萬km2。寶雞屬暖溫帶半濕潤氣候,受東亞季風(包括高原季風)控制,造就了冬季寒冷干燥,夏季溫熱多雨與炎熱干燥交替,春季多變少雨且升溫迅速,秋季多陰雨連綿且降溫迅速的全年氣候變化特征。氣溫、降水量季相較為分明。年平均氣溫13 ℃左右,年降水量710~1000 mm,7、8、9這3個月降水集中,占全年降水的50%左右;寶雞地區(qū)是關(guān)中降水量最多的地區(qū)。寶雞地區(qū)生境條件和植被類型多樣,主要包括落葉闊葉林、針葉林、草原、灌叢和栽培植被等。森林覆被率達36%~42%,但分布不均,以秦嶺和關(guān)山為集中分布區(qū)域,關(guān)中和臺塬地區(qū)主要為栽培植被。
圖1 寶雞地區(qū)概況及氣象站點分布
表1為本研究主要數(shù)據(jù)源信息,包括氣象數(shù)據(jù)(11個氣象站點)、MODIS NDVI、NPP和高程DEM數(shù)據(jù)。研究中站點選取分別是:寶雞市、寶雞縣、鳳翔縣、岐山縣、眉縣、扶風縣、隴縣、千陽縣、麟游縣、鳳縣和太白縣,本研究也采用此11個站點作為行政區(qū)域進行下文研究。
2.2.1 趨勢分析 通過模擬每個柵格的變化趨勢和單個像元時間變化特征來反映整個空間變化規(guī)律,進而表征區(qū)域時空格局演變[17]。Slope為NDVI在各個像元點的變化趨勢線斜率。若Slope值為正,則區(qū)域植被NDVI隨時間變化呈上升趨勢,且數(shù)值大小與植被覆蓋度增加趨勢強弱呈正比例關(guān)系;若Slope值為負,則表示植被指數(shù)隨時間變化呈下降趨勢。其計算公式如下:
(1)
公式(1)中:n為年份,i為年序列,NDVIi為第i年的最大化NDVI值。
表1 研究所用數(shù)據(jù)源
2.2.2 相關(guān)分析法 運用相關(guān)分析法分析植被覆蓋變化對氣候因子的響應(yīng),結(jié)合ArcGIS10.1和Matlab軟件,計算NDVI與氣溫和降水量間的相關(guān)系數(shù),采用NDVI影像與經(jīng)反距離權(quán)重插值法獲取年均溫和年降水量的柵格數(shù)據(jù),繪制NDVI與年均溫和年降水的相關(guān)系數(shù)空間分布圖,NDVI與氣候因子的空間關(guān)系通過顯著性檢驗來進行分析。計算公式如下[18]:
(2)
經(jīng)ArcGIS區(qū)域統(tǒng)計功能獲取年均NDVI,表示寶雞地區(qū)當年植被覆蓋水平(圖2)。2001~2013年NDVI變化情況為:2006~2007、2008~2010、2011~2012年緩慢增長,2001~2002、2003~2004年飛速增長,2004~2006、2007~2008、2010~2011年緩慢下降,2002~2003年迅速下降,2012~2013年變化較不明顯??傮w來看,NDVI以0.032/10 a的速度呈上升趨勢。三北防護林工程區(qū)1982~2006年近25 a植被覆蓋的平均增速為0.007/10 a,與此相比,寶雞地區(qū)植被覆蓋增速較快[19],且NDVI值在2010年達到最高值0.572。2001年寶雞地區(qū)平均NDVI值為0.5274,2013年增加至0.5668,增速為7.47%。說明寶雞地區(qū)近13 a來的植被覆蓋有增無減。
寶雞地區(qū)2001~2013年多年平均NDVI空間分布如圖3所示,整體上呈從東北向西南遞增的特征。由于太白山的山地地貌類型和其為秦嶺山脈主峰的原因,以及太白河和水河的河水補給,寶雞地區(qū)植被覆蓋最高值集中分布在太白縣南部。千河沿岸以及渭河流域?qū)氹u段的寶雞市、寶雞縣、鳳翔、岐山南北沿岸的“帶狀”地區(qū),植被覆蓋度最低,主要原因是這些地區(qū)人口集聚分布,經(jīng)濟發(fā)展最為迅速,隨著城市化的推動,城市擴建對于土地需求逐日擴大,以0.49%/a的速度增加,建筑面積的迅速增加極大地削減了植被面積。此外,大量的農(nóng)田建設(shè)和農(nóng)作物種植呈現(xiàn)出植被覆蓋明顯的季節(jié)分布特征,從1999年陜西地區(qū)退耕還林試點實施以來,寶雞地區(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變化趨勢
寶雞地區(qū)2001~2013年植被覆蓋變化具有明顯的區(qū)域差異(圖4)。由于近年來寶雞地區(qū)正在大力實施退耕還林還草和保護森林的政策[20],植被恢復(fù)較好的區(qū)域廣泛分布于千河沿岸的隴縣、千陽、寶雞縣地區(qū),漆水河以北的麟游地區(qū),嘉陵江沿岸的鳳縣西南地區(qū),狀況逐年好轉(zhuǎn)。而近年來城市化水平逐年提升,城市土地擴張需求速度有增無減,所以大量的城區(qū)建筑用地占用了耕地及草地資源,使得千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸的植被覆蓋度降低。
為了更好地表征寶雞地區(qū)植被恢復(fù)狀況,參照已有研究[17,21-23],將Slope分為如圖4、表2所示的7個等級。結(jié)果顯示,嚴重退化面積1 km2,中度退化面積117 km2,分別占總面積的0.006%和0.646%,兩個等級的植被覆蓋分布狀況集中在千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸“帶狀”地區(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%,兩者集中于千河沿岸的隴縣、千陽、寶雞縣地區(qū)、漆水河以北的麟游地區(qū)、嘉陵江沿岸的鳳縣西南部地區(qū)。綜上所述,近13 a寶雞地區(qū)植被覆蓋狀況呈穩(wěn)中上升的趨勢,這主要是因為寶雞地區(qū)實施“退耕還林還草”、森林保護等政策,使得該地綠化城市面積和造林面積分別以0.91%/a和0.73%/a的增速改善,以及當?shù)厝藗儽Wo環(huán)境意識的不斷增強,植被覆蓋狀況逐年見好,由此生態(tài)環(huán)境也得到逐年改善[24]。
圖3 年均NDVI變化斜率趨勢圖
NDVI數(shù)據(jù)的提取是建立在氣候要素對植被前期生長作用的基礎(chǔ)之上的,因而研究NDVI和氣溫、降水量在植被生長期內(nèi)的關(guān)系具有重要的實際意義[13-14,25-27]。由圖6和表3可知,近13 a來寶雞地區(qū)植被覆蓋與氣溫呈正相關(guān)的區(qū)域像元占總像元的80.99%,降水量呈正相關(guān)的區(qū)域像元占79.01%,表明年均NDVI與年均溫和年降水的相關(guān)性整體上表現(xiàn)均較好。
植被覆蓋與氣溫呈正相關(guān)的區(qū)域像元數(shù)占總像元的80.99%,除鳳縣大部分地區(qū)和千陽中西部地區(qū)外,廣泛分布于寶雞地區(qū)各個縣區(qū)。由于以高山草甸為主體的關(guān)山草原位于隴縣西南部、太白縣中部至東部,秦嶺主峰穿其境而過且有太白河和水河水分的補給,因此其相應(yīng)區(qū)域相關(guān)系數(shù)最大;千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸“帶狀”地區(qū),地勢平坦、氣溫四季分明,土壤有機質(zhì)的分解在夏季氣溫升高的條件下加速,對常溫地區(qū)植被凈初級生產(chǎn)力的提高起到了有效的促進作用[28-30],又加之大量植被綠化面積的建設(shè),所以其相關(guān)系數(shù)最大(表3)。植被覆蓋與氣溫呈負相關(guān)的區(qū)域集中分布在鳳縣大部分地區(qū)和千陽中西部地區(qū),原因可能是鳳縣地處秦嶺山地地區(qū),且紫柏山坐落其境,另外有嘉陵江和褒河水分的補給。吳岳支脈箭筈嶺坐落于千陽南部,兩個地區(qū)均海拔較高、年均溫較低。當植物生長的最適溫度低于夏季氣溫時,植株內(nèi)Ru Bis CO(核酮糖-1,5-二磷酸羧化酶/加氧酶)的羧化反應(yīng)小于其加氧反應(yīng),加速了凈光合作用的下降速度[31];同時,土壤水分蒸發(fā)和植物蒸騰作用會因溫度的升高而加劇,嚴重影響了植被生長[32]??傮w而言,介于0.3~0.9相關(guān)系數(shù)的像元均勻分布在寶雞各個地區(qū),說明寶雞地區(qū)植被覆蓋與氣溫的相關(guān)性表現(xiàn)較好。
圖4 年均NDVI變化斜率分級圖
植被覆蓋與降水量呈負相關(guān)的區(qū)域集中在渭河流域?qū)氹u段的千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸“帶狀”地區(qū)(表3),這些地區(qū)具有較強的降水量,而隨著城市化進程的加快,城市擴展帶來建筑面積擴大,植被覆蓋度不斷降低,因此植被覆蓋與降水呈負相關(guān);除此之外,寶雞其他各個縣區(qū)的植被覆蓋與降水量均呈正相關(guān),其中鳳縣大面積地區(qū)呈非常顯著正相關(guān)。這是因為水是植物生長過程中的決定性因素,植物的生長狀況會直接受到降水量的影響。降水量是影響我國干旱變化最主要的因素[33],降水量對植物生長的限制程度尤以年降水量少的干旱、半干旱地區(qū)表現(xiàn)最為明顯[10]。綜上可見,寶雞地區(qū)植被覆蓋與降水量的相關(guān)性表現(xiàn)較好。
表2 2001~2013年寶雞地區(qū)年均NDVI變化趨勢分級統(tǒng)計
表3 2001~2013年寶雞地區(qū)年均NDVI與氣溫和降水量相關(guān)性的像元比例
寶雞地區(qū)年均NPP為668.0831 g·C/(m2·a),受降水量影響(年平均降水量529.89 mm),2001、2002、2004年千陽的NPP值均呈現(xiàn)最低,鳳翔在2003、2005~2010共7年的NPP均為最低。而11個行政區(qū)域2001~2010年NPP的最高值均出現(xiàn)在鳳縣。表4數(shù)據(jù)顯示,2001~2010年對年均NPP影響較大的因子為NDVI。在12個分析結(jié)果中,以麟游的顯著性最高,達到0.01顯著水平,岐山、隴縣、千陽、鳳縣次之,達到0.05顯著水平,且均為正相關(guān)。依據(jù)生態(tài)學相關(guān)理論,NDVI與NPP具有直接或間接的相互關(guān)系,在區(qū)域和生態(tài)系統(tǒng)模型構(gòu)建中,NDVI常被直接或間接地用于計算植被NPP[34-35]。因此,NDVI與NPP具有同步增減的變化趨勢,即隨著植被覆蓋率的增加(減少)而NPP相應(yīng)上升(下降)。
氣溫和降水量與NPP的相關(guān)性均較差,均未達到顯著性水平,氣溫除鳳翔、麟游地區(qū)外其他地區(qū)均呈負相關(guān),降水除岐山外其余地區(qū)均呈正相關(guān)。相關(guān)系數(shù)的顯著性均不如前50 a的結(jié)果[36]顯著,這一數(shù)據(jù)分析結(jié)果也在個別站點變化的差異較大。2000年以來國家推行天然林保護政策,寶雞地區(qū)特別是渭河沿岸及支流周邊區(qū)域陸續(xù)實行了包括封山育林、加強水源涵養(yǎng)及天然次生林保護、生態(tài)修復(fù)等退耕還林還草工程,寶雞地區(qū)在陜西省退耕還林試行初期(1999年)森林覆蓋率達34.30%,之后以1.43%/5 a的速度增長,森林單位蓄積量也以385 m3/5 a的速度增加。與此同時,氣溫和降水等氣象因子也在逐年發(fā)生變化,在如此復(fù)雜的人類活動和氣候變化的綜合作用下,本研究選取渭河流域中上游的寶雞地區(qū),研究范圍較小,更具典型代表意義,對進一步理清植被恢復(fù)與NPP變化的響應(yīng)關(guān)系具有良好的科學和實際意義。
表4 2001~2010年不同行政區(qū)域NPP與NDVI及氣候因子間相關(guān)系數(shù)統(tǒng)計
注:*表示在0.05水平上顯著相關(guān);**表示在0.01水平上顯著相關(guān)。
寶雞地區(qū)相比于干旱地區(qū),水熱條件適中,由于水資源有一定保障,因而對于植被生長而言溫度更易成為其主要驅(qū)動因素[37-38]。理論上而言,日積溫伴隨著溫度的上升而增加,對大多數(shù)植被生長顯然有積極作用,鳳縣大部分地區(qū)(秦嶺、紫柏山、嘉陵江、褒河)和千陽中西部地區(qū)(吳岳支脈箭筈嶺)植被覆蓋程度大幅提高,生態(tài)恢復(fù)效果最為明顯,但植被與氣溫卻呈現(xiàn)負相關(guān)關(guān)系??梢姡瑖彝烁€林政策的實施使這些山地地區(qū)林草用地增加,單株植被的長勢與整體植被覆蓋程度的關(guān)系不大,因此,溫度的驅(qū)動作用很難體現(xiàn)出來。同時,寶雞具有以山地、丘陵為主的復(fù)雜地形,在海拔較高地區(qū)淺層地下水儲存不足,天然降水是植被生長需水依賴的主要來源,從而形成了寶雞地區(qū)植被與降水量的相關(guān)性中東部低于北、西、南部的格局;而在千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸“帶狀”地區(qū),由于地勢平坦,人類活動在城鎮(zhèn)化進程中強度增大,植被生長受其間接影響與降水量呈負相關(guān)發(fā)展趨勢。因此,植被對溫度和降水的負相關(guān)并不代表需求為負,而應(yīng)考慮退耕還林、城市建設(shè)、土地利用等人類活動的客觀因素影響。
“退耕還林還草”、森林保護等政策實施以來,寶雞地區(qū)綠化城市面積增速(0.91%/a)和造林面積增速(0.73%/a)較好,表現(xiàn)在千河沿岸的隴縣、千陽縣、寶雞縣沿岸地區(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ū)和千陽中西部地區(qū)植被覆蓋與年均溫呈現(xiàn)負相關(guān),但鳳縣大部分地區(qū)與年降水量呈正相關(guān),且正相關(guān)系數(shù)顯著;千河與渭河交叉口沿渭河向東方向的寶雞縣、鳳翔、岐山、眉縣、扶風南北沿岸“帶狀”地區(qū),植被覆蓋與年降水量呈負相關(guān),但此地帶植被覆蓋與氣溫呈顯著正相關(guān),以千河與渭河交叉口的寶雞縣東北部與鳳翔縣東南部交接地區(qū)顯著性最為明顯。
退耕后10 a NDVI與NPP呈良好的正相關(guān)性,以麟游縣最為顯著(0.01水平上顯著),即凈初級生產(chǎn)力與植被覆蓋具有同步增減的變化趨勢。
[1] 方精云,宋永昌,劉鴻雁,等.植被氣候關(guān)系與我國的植被分區(qū)[J].植物學報,2002,44(9):1105-1122.
[2] 趙宗慈,王紹武,羅勇.IPCC成立以來對溫度升高的評估與預(yù)估[J].氣候變化研究進展,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è)科學,2006,34(18):4618-4619.
[5] 韓愛惠.用MODIS植被指數(shù)的時間序列分析提取荒漠化敏感區(qū)域的方法[J].林業(yè)資源管理,2004(1):57-60.
[6] 朱文泉,陳云浩,徐丹,等.陸地植被凈初級生產(chǎn)力計算模型研究進展[J].生態(tài)學雜志,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)凈初級生產(chǎn)力及其與氣候的關(guān)系[J].植物生態(tài)學報,2010,34(7):781-791.
[11] 周夏飛,馬國霞,曹國志,等.2001~2013年黃土高原植被凈初級生產(chǎn)力時空變化及其歸因[J].安徽農(nóng)業(yè)科學,2017,45(14):48-53.
[12] 楊尚武,張勃.基于SPOT NDVI的甘肅河東植被覆蓋變化及其對氣候因子的響應(yīng)[J].生態(tài)學雜志,2014,33(2):455-461.
[13] 崔曉臨,白紅英,王濤.秦嶺地區(qū)植被NDVI海拔梯度差異及其氣溫響應(yīng)[J].資源科學,2013,35(3):618-626.
[14] 戴聲佩,張勃,王海軍,等.中國西北地區(qū)植被時空演變特征及其對氣候變化的響應(yīng)[J].遙感技術(shù)與應(yīng)用,2010,25(1):69-76.
[15] 秦超,李君軼,陳宏飛.近20年寶雞市植被覆蓋度動態(tài)變化及驅(qū)動力分析[J].山東農(nóng)業(yè)科學,2014,46(9):98-105.
[16] 翟雅倩,周旗,張翀,等.近15年寶雞市植被覆蓋度變化研究[J].寶雞文理學院學報:自然科學版,2016,36(4):31-36.
[17] 李雙雙,延軍平,萬佳.近10年陜甘寧黃土高原區(qū)植被覆蓋時空變化特征[J].地理學報,2012,67(7):960-970.
[18] 徐建華.現(xiàn)代地理學中的數(shù)學方法[M].北京:高等教育出版社,2002:37-68.
[19] 王強,張勃,戴聲佩,等.基于GIMMSAVHRR NDVI數(shù)據(jù)的三北防護林工程區(qū)植被覆蓋動態(tài)變化[J].資源科學, 2011,33(8):1613-1620.
[20] 周德成,趙淑清,朱超.退耕還林工程對黃土高原土地利用/覆被變化的影響:以陜西省安塞縣為例[J].自然資源學報,2011,26(11):1866-1878.
[21] 王震,閆文德,劉曙光,等.基于MODIS數(shù)據(jù)的中國三種主要土地類型變化的空間特征分析[J].生態(tài)學報,2017,37(10):1-7.
[22] 張君,延軍平.1982~2013年陜西不同植被類型NDVI變化特征分析[J].干旱區(qū)資源與環(huán)境,2017,31(4):86-92.
[23] 宋怡,馬明國.基于SPOT VEGETATION數(shù)據(jù)的中國西北植被覆蓋變化分析[J].中國沙漠,2007,27(1):89-93.
[24] 宋德明.陜西省地理志[M].西安:陜西人民出版社,2000:221- 262.
[25] 金麗娜.西安近10年植被覆蓋變化及其與氣候因子的相關(guān)性研究[J].甘肅科學學報,2014,26(3):23-26.
[26] 蒲蕾,任志遠.陜西省不同地區(qū)NDVI變化與氣候因子的關(guān)系及響應(yīng)研究[J].水土保持通報,2013,32(2):265-275.
[27] 方精云,樸世龍,賀金生,等.近20年中國植被活動在增強[J].中國科學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ū)植被覆蓋變化及其對氣候因子的響應(yīng)研究[J].農(nóng)業(yè)現(xiàn)代化研究,2013,34(1):104-108.
[32] 賈東于,文軍,張?zhí)锰?等.黃土高原降水對土壤含水量和導(dǎo)熱率的影響[J].高原氣象,2014,33(3):712-720.
[33] 翟盤茂,鄒旭愷.1951~2003年中國氣溫和降水變化及其對干早的影響[J].氣象變化研究進展,2005,1(1):16-18.
[34] 周珺.基于遙感數(shù)據(jù)的重慶市凈初級生產(chǎn)力(NPP)時空特征研究[D].重慶:西南大學,2013.
[35] 朱文泉,潘耀忠,何浩,等.中國典型植被最大光利用率模擬[J].科學通報,2006,51(6):700-710.
[36] 蔣沖,王飛,穆興民.渭河流域相對濕度與植被覆蓋變化及其響應(yīng)關(guān)系[J].水土保持通報,2012,32(5):122-127.
[37] 崔林麗,史軍,楊引明,等.中國東部植被NDVI對氣溫和降水的旬響應(yīng)特征[J].地理學報,2009,64(7):850-860.
[38] 王靜,萬紅蓮,張翀. 2001-2013年寶雞地區(qū)植被NDVI對氣溫和降水量的季節(jié)響應(yīng)特征[J].水土保持通報,2017,37(5):235-240,349.