郭力宇 郭昭 王濤 張?jiān)? 李雨思 柳瑩

摘要：為了研究退耕還林還草工程實(shí)施后期，陜北黃土高原地區(qū)地表植被變化特征及氣候因子對(duì)植被變化的影響，采用MODIS NDVI數(shù)據(jù)，通過(guò)像元二分法計(jì)算植被覆蓋度，結(jié)合相關(guān)氣象站點(diǎn)監(jiān)測(cè)數(shù)據(jù)，利用線性回歸分析、相關(guān)系數(shù)及顯著性檢驗(yàn)法，分析了陜北黃土高原植被覆蓋度時(shí)空變化趨勢(shì)及其對(duì)氣溫、降水因子的響應(yīng)。結(jié)果表明：研究區(qū)植被覆蓋度空間分布差異顯著，由北向南植被覆蓋度由低向高過(guò)渡明顯，2009—2016年植被覆蓋等級(jí)的空間轉(zhuǎn)移變化主要發(fā)生在北部和中部地區(qū)，表現(xiàn)出由中低植被覆蓋向中等植被覆蓋轉(zhuǎn)移、中等植被覆蓋向中高植被覆蓋轉(zhuǎn)移的態(tài)勢(shì)。2009—2016年研究區(qū)植被覆蓋度年際變化速率較小，整體呈平緩增加趨勢(shì)，中部和南部部分地區(qū)有減少趨勢(shì)但變化趨勢(shì)不顯著。研究區(qū)內(nèi)年植被覆蓋度與年降水量、年平均氣溫的相關(guān)系數(shù)、偏相關(guān)系數(shù)均較小，相關(guān)性不顯著。不同區(qū)域的植被覆蓋度對(duì)氣候因子的響應(yīng)結(jié)果不同，在空間上正負(fù)相關(guān)性共存，研究區(qū)北部受年降水量影響較大，而中部受年平均氣溫影響較大。研究區(qū)植被覆蓋度對(duì)氣候因子的響應(yīng)存在滯后性，生長(zhǎng)季的植被生長(zhǎng)受同月降水量和前一月平均氣溫的影響較大，且絕大部分地區(qū)均表現(xiàn)為正向相關(guān)。

關(guān)鍵詞：退耕還林還草;植被覆蓋度;動(dòng)態(tài)變化;氣候響應(yīng);陜北黃土高原

中圖分類號(hào)：TP 79文獻(xiàn)標(biāo)志碼：A

DOI：10.13800/j.cnki.xakjdxxb.2019.0219文章編號(hào)：1672-9315（2019）02-0317-10

0引言

植被是反映陸地生態(tài)環(huán)境變化的重要因子之一，在土壤、水分、大氣物質(zhì)和能量循環(huán)過(guò)程中扮演重要角色[1-2]。植被覆蓋度（Fractional Vegetation Coverage，F(xiàn)VC）是衡量地表植被覆蓋信息的重要指標(biāo)，是指植被冠層在地面上垂直投影面積與土地總面積的比值[3]。隨著遙感技術(shù)發(fā)展，歸一化差異植被指數(shù)（Normalized Differential Vegetation Index，NDVI）常用來(lái)監(jiān)測(cè)區(qū)域植被覆蓋變化[4-5]，估算植被覆蓋度[6]。特別是MODIS提供的NDVI數(shù)據(jù)產(chǎn)品具有較高的時(shí)間分辨率[7]，廣泛應(yīng)用于地表植被變化研究中。戚德康等采用MODIS NDVI數(shù)據(jù)研究了2000—2016年朝鮮植被覆蓋變化特征[8];何寶忠等基于MODIS NDVI數(shù)據(jù)，利用像元二分模型估算新疆2005—2015年的植被覆蓋度，研究了新疆不同生態(tài)區(qū)植被覆蓋變化趨勢(shì)[9]。在具體的研究中，線性回歸法[10]、相關(guān)系數(shù)法[11]也已經(jīng)成為研究植被覆蓋變化特征、分析氣候因子與植被覆蓋相關(guān)性的常用研究方法。

陜北黃土高原地區(qū)生態(tài)環(huán)境脆弱，水土流失嚴(yán)重，是我國(guó)生態(tài)環(huán)境建設(shè)與保護(hù)的重點(diǎn)區(qū)域。1999年陜西省開(kāi)始試點(diǎn)推行退耕還林還草工程，隨后國(guó)家正式發(fā)布《退耕還林工程規(guī)劃（2001—2010年）》。自1999—2008年陜北地區(qū)榆林和延安兩地累計(jì)退耕造林面積達(dá)4 721 km2，植被覆蓋度顯著提高，退耕還林還草活動(dòng)成為影響陜北黃土高原植被覆蓋變化的主要驅(qū)動(dòng)因素。

當(dāng)前與陜北黃土高原植被覆蓋變化相關(guān)的研究已有一些成果[12-15]，多集中在退耕還林工程實(shí)施過(guò)程的10年前后，多以年際尺度人類活動(dòng)或自然環(huán)境因子對(duì)植被覆蓋變化影響為主，如范建忠等研究表明受退耕還林還草工程實(shí)施影響，陜北地區(qū)2000—2011年退耕還林區(qū)植被覆蓋度顯著增加[16];韓可欣等研究了2001—2010年陜北黃土高原不同坡度區(qū)植被覆蓋變化特征[17];肖強(qiáng)等研究了2001—2010年黃土高原植被覆蓋變化與氣候因子的相關(guān)性[18]。這些研究為陜北黃土高原生態(tài)環(huán)境建設(shè)發(fā)揮了重要作用，但研究中較少以月際尺度進(jìn)行研究，且很少關(guān)注退耕還林后期（即2001—2010年工程規(guī)劃結(jié)束之后），大范圍人為活動(dòng)減弱背景下，陜北黃土高原植被覆蓋年際和月際尺度上的變化特征。因此，文中選取2009—2016年為研究時(shí)段，利用MODIS NDVI數(shù)據(jù)結(jié)合相關(guān)氣象站監(jiān)測(cè)資料，計(jì)算植被覆蓋度，并分別從年際和月際尺度上分析退耕還林還草后期陜北黃土高原植被變化特征及其對(duì)氣候因子的響應(yīng)，為新一輪退耕還林還草工程實(shí)施與政策制定提供參考。

1研究區(qū)概況

陜北黃土高原地區(qū)是中國(guó)黃土高原的核心區(qū)域，位于107°15′13″E～111°14′40″E，35°04′10″N～39°15′13″N之間，依據(jù)陜西省測(cè)繪地理信息局發(fā)布的《陜西地理省情白皮書(shū)2011》中對(duì)陜北、關(guān)中、陜南的劃分成果，陜北地區(qū)面積約8.3×104 km2，包括榆林和延安兩市的3區(qū)22縣以及銅川市和渭南市北部部分鄉(xiāng)鎮(zhèn)。該區(qū)域?qū)侔敫稍鐪貛Т箨懠撅L(fēng)氣候，處于半干旱和干旱過(guò)渡帶，多年平均溫度8～14 ℃，多年平均降雨量375～613 mm，水熱條件時(shí)空分配不均，呈現(xiàn)出由南而北逐漸減少的變化趨勢(shì)，且降水集中于夏秋2季。植被處于暖溫帶落葉闊葉林和溫帶草原之間，過(guò)渡性明顯，由北向南依次為草原區(qū)、森林草原區(qū)以及落葉闊葉林區(qū)[19]。地勢(shì)西高東低，北部為毛烏素沙漠南緣，地形平坦，植被稀疏，中部為黃土梁峁區(qū)，生態(tài)環(huán)境脆弱，地表植被人類活動(dòng)影響較大，南部山區(qū)植被覆蓋較好。

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

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

研究使用數(shù)據(jù)包括：陜北黃土高原2009—2016年16 d合成的250 m分辨率MODIS NDVI 13Q1Level 3數(shù)據(jù)，下載自http：∥ladsweb.nasa.gov.陜北黃土高原及其附近22個(gè)氣象站點(diǎn)2009—2016年逐月氣溫、降水?dāng)?shù)據(jù)，下載自中國(guó)氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn），站點(diǎn)分布如圖1所示。

采用最大值合成法（Maximum Value Composition，MVC）獲取2009—2016年MODIS NDVI月值和年值影像。采用陜北黃土高原矢量邊界裁剪獲得2009—2016年月值和年值影像。

水熱條件是影響植被生長(zhǎng)的主要因子。文中整理了研究區(qū)及其附近22個(gè)氣象站點(diǎn)逐月平均氣溫、降水量數(shù)據(jù)，獲得月平均和年平均氣溫、月和年降水量數(shù)據(jù)。采用反距離權(quán)重法進(jìn)行空間插值獲取2009—2016年逐月和逐年平均氣溫、降水量空間分布數(shù)據(jù)。

3結(jié)果與分析

3.1植被覆蓋等級(jí)變化特征

基于研究區(qū)特征，將陜北黃土高原植被覆蓋度分為低、中低、中等、中高以及高植被覆蓋度5類等級(jí)，閾值分別為20%，40%，60%，80%，獲得2009—2016年逐年陜北黃土高原植被覆蓋度等級(jí)的空間分布變化（圖2），并統(tǒng)計(jì)分析各年份不同植被覆蓋度等級(jí)的面積變化（圖3）。

從圖2和圖3可以看出，低植被覆蓋度區(qū)與中低植被覆蓋度區(qū)面積變化曲線相似，2009—2013年面積減少趨勢(shì)明顯，2014—2016年先增后減，2015年的峰值分別為11 597，23 892 km2.中等植被覆蓋區(qū)面積最大，2009—2016年面積變化幅度較小，多數(shù)年份在25 000 km2以上，2015年出現(xiàn)低值18 838 km2.中高植被覆蓋區(qū)2009—2013年面積增加明顯，2014—2016年先減少后增加，2015年的低值為13 032 km2.高植被覆蓋區(qū)的面積變化不明顯，面積均值為16 454 km2，占比約19.82%.

陜北黃土高原植被覆蓋由北向南表現(xiàn)出極顯著的過(guò)渡性：研究區(qū)北部以低植被覆蓋與和中低植被覆蓋為主，研究區(qū)中部自北向南，呈現(xiàn)出由中等植被覆蓋向中高植被覆蓋過(guò)渡，南部則為高植被覆蓋集中區(qū)域。2009—2016年研究區(qū)植被覆蓋等級(jí)的空間轉(zhuǎn)移變化主要發(fā)生在北部和中部地區(qū)，除2015年發(fā)生短暫的植被覆蓋減少之外，總體表現(xiàn)為由中低植被覆蓋向中等植被覆蓋轉(zhuǎn)移、中等植被覆蓋向中高植被覆蓋轉(zhuǎn)移的態(tài)勢(shì)。這反映出第1輪退耕還林還草工程實(shí)施結(jié)束后，陜北黃土高原植被覆蓋度仍呈持續(xù)增加過(guò)程。

3.2植被覆蓋動(dòng)態(tài)變化趨勢(shì)

3.2.1月際變化趨勢(shì)

2009—2016年陜北黃土高原月際植被覆蓋度呈周期性變化過(guò)程（圖4），每年2～3月為全年植被覆蓋度最低的月份，小于5%;7～9月為全年植被覆蓋度最高的月份，均在40%以上，其中2012年和2013年7～8月的植被覆蓋度均在55%左右。全年1月、2月、3月、11月以及12月的植被覆蓋度均低于10%，4～9月份為生長(zhǎng)季[26]，植被生長(zhǎng)明顯，植被覆蓋度顯著增加。2009—2016年，植被覆蓋度在45%以上的月份，從2009年的2個(gè)增加到2016年的3個(gè)，其中2014年植被長(zhǎng)勢(shì)最好，有4個(gè)月的植被覆蓋度均在45%以上。月際尺度上，研究區(qū)植被覆蓋度具有平穩(wěn)增長(zhǎng)的趨勢(shì)。

3.2.2年際變化趨勢(shì)

采用一元線性模型在像元尺度計(jì)算2009—2016年研究區(qū)植被覆蓋度的年際變化趨勢(shì)，并統(tǒng)計(jì)像元值的分布頻數(shù)（圖5）?？芍?，研究區(qū)植被覆蓋度的年際變化率處于-0.151～0.159 a-1之間，年際變化率低且年際變化率在-0.04～0.04 a-1之間的區(qū)域面積達(dá)到82 513 km2，占研究區(qū)總面積的99.4%.因而，可將陜北黃土高原植被覆蓋年際變化率劃分為6個(gè)等級(jí)（圖6（a））：b≤-0.04 a-1，-0.04 a-1

統(tǒng)計(jì)不同等級(jí)下的面積，可知：年際變化率小于-0.04 a-1和大于0.04 a-1的區(qū)域面積分別為152 km2和365 km2，占比極小，均不足0.5%;年際變化率在-0.04～-0.02 a-1之間的區(qū)域面積與0.02～0.04 a-1之間的區(qū)域面積的較少，分別為2 578，5 806 km2，占比約3.1%和7%;年際變化率在-0.02～0 a-1之間的區(qū)域面積為31 183 km2，占比約37.56%;而年際變化率在0～0.02 a-1之間的區(qū)域面積最大，達(dá)到42 946 km2，占比約為5172%.

可知，研究區(qū)的植被覆蓋度以增長(zhǎng)為主，中部和南部部分地區(qū)植被覆蓋度呈減少趨勢(shì)，并且8928%的區(qū)域處于-0.02～0.02 a-1之間，總體上植被覆蓋度的年際變化速率較小。

對(duì)年際植被覆蓋度變化趨勢(shì)采用F檢驗(yàn)進(jìn)行顯著性檢驗(yàn)（圖6（b）），根據(jù)檢驗(yàn)結(jié)果分別在置信度α=0.05和α=0.01下將變化趨勢(shì)劃分為5個(gè)等級(jí)：極顯著減少（b<0，P<0.01）、顯著減少（b<0，0.010.05）、顯著增加（b>0，0.010，P<0.01）。

統(tǒng)計(jì)不同等級(jí)區(qū)域面積，可知：極顯著和顯著減少區(qū)域面積為1 325.76 km2，占比不足2%;顯著和極顯著增加區(qū)域面積為7 099.18 km2，占比為8.55%;變化不顯著的區(qū)域面積達(dá)74 604.56 km2，占比為89.85%.綜上可知，2009—2016年陜北黃土高原植被覆蓋度以平緩增加趨勢(shì)為主，部分區(qū)域呈平緩減少趨勢(shì)，少數(shù)地區(qū)呈顯著減少或增加趨勢(shì)，其中植被覆蓋度增加顯著的區(qū)域主要分布在北部和中部，顯著減少區(qū)域呈散點(diǎn)狀分布，集中在城市附近區(qū)域。

3.3植被覆蓋度對(duì)氣候因子的響應(yīng)

采用相關(guān)系數(shù)和偏相關(guān)系數(shù)、t檢驗(yàn)計(jì)算2009—2016年植被覆蓋度與年平均氣溫、年降水量之間的相關(guān)系數(shù)、偏相關(guān)系數(shù)和偏相關(guān)系數(shù)t檢驗(yàn)結(jié)果。統(tǒng)計(jì)相關(guān)系數(shù)和偏相關(guān)系數(shù)的平均值可得表1.t檢驗(yàn)結(jié)果查表可知，置信度α=0.05下，臨界值為2.571，根據(jù)臨界值將相關(guān)系數(shù)劃分為4個(gè)等級(jí)：即顯著負(fù)相關(guān)（t<-2.571）、負(fù)相關(guān)（-2.571≤t<0）、正相關(guān)（0≤t≤2.571）以及顯著正相關(guān)（t>2.571）（圖7）。

coefficient between FVC and mean temperature and precipitation相關(guān)系數(shù)偏相關(guān)系數(shù)相關(guān)系數(shù)偏相關(guān)系數(shù)平均氣溫0.3140.375降水量0.4390.477表1和圖7表明：2009—2016年研究區(qū)植被覆蓋度與年平均氣溫和降水量的相關(guān)系數(shù)平均值分別為0.314和0.439，偏相關(guān)系數(shù)平均值為0375和0.477，表明研究區(qū)年際植被覆蓋度與年降水量的相關(guān)性略高于年平均氣溫，且不同區(qū)域的植被覆蓋度對(duì)氣候因子的響應(yīng)存在較大差異。

植被覆蓋度受年降水量影響呈正相關(guān)的區(qū)域面積為62 558 km2，占總面積的75.34%，其中植被覆蓋度與年降水量呈顯著正相關(guān)的區(qū)域面積為13 392 km2，約占16.13%，主要集中在中部和北部。植被覆蓋度與年平均氣溫呈正相關(guān)的區(qū)域面積小于呈負(fù)相關(guān)的區(qū)域面積，兩者分別為36 373，40 822 km2.其中呈顯著負(fù)相關(guān)的區(qū)域面積為3 041 km2，約占3.66%，主要集中在中部;呈顯著正相關(guān)的區(qū)域分布分散，面積為2 792 km2，約占3.36%.

研究區(qū)內(nèi)的年植被覆蓋度與年降水量以正相關(guān)為主，與年平均氣溫以負(fù)相關(guān)為主，但兩者的相關(guān)系數(shù)、偏相關(guān)系數(shù)均較小，相關(guān)性不顯著。但不同區(qū)域植被覆蓋度對(duì)氣候因子的響應(yīng)不同，其中北部植被覆蓋度受年降水量影響較大，而中部受年平均氣溫影響較大。

3.4植被覆蓋度對(duì)氣候因子響應(yīng)的時(shí)滯性

考慮氣候因子對(duì)植被生長(zhǎng)的影響存在著時(shí)滯性，以2009—2016年植被覆蓋度和月氣象數(shù)據(jù)為基礎(chǔ)，采用生長(zhǎng)季（4～9月份）月植被覆蓋度與同月（4～9月）平均氣溫和降水量及前一月（3～8月）平均氣溫和降水量進(jìn)行相關(guān)系數(shù)和偏相關(guān)系數(shù)計(jì)算，并采用t檢驗(yàn)進(jìn)行顯著性檢驗(yàn)。α=0.05時(shí)，查表可知臨界值為3.182，將相關(guān)系數(shù)和偏相關(guān)系數(shù)的t檢驗(yàn)結(jié)果分為4個(gè)等級(jí)：顯著負(fù)相關(guān)（t<-3.182）、負(fù)相關(guān)（-3.182≤t<0）、正相關(guān)（0≤t≤3.182）以及顯著正相關(guān)（t>3.182），得到圖8，統(tǒng)計(jì)相關(guān)系數(shù)和偏相關(guān)系數(shù)的平均值可得表2.

偏相關(guān)系數(shù)平均氣溫0.5390.855降水量0.6430.473由表2可知，生長(zhǎng)季（4～9月）植被覆蓋度與同月平均氣溫和降水量的相關(guān)系數(shù)平均值分別為0.770和0.802，偏相關(guān)系數(shù)平均值分別為0.539和0.643;與前一月平均氣溫和降水量的相關(guān)系數(shù)平均值分別為0.914和0.731，偏相關(guān)系數(shù)平均值分別為0.855和0.473.表明研究區(qū)生長(zhǎng)季植被覆蓋度受同月降水量的影響大于前一月的降水量，與前一月平均氣溫的偏相關(guān)系數(shù)遠(yuǎn)大于同月的平均氣溫，研究區(qū)的植被覆蓋度對(duì)氣溫的響應(yīng)存在明顯時(shí)滯性，生長(zhǎng)季的植被生長(zhǎng)受當(dāng)月降水量和前一月平均氣溫影響較大。

從圖8可以看出，研究區(qū)生長(zhǎng)季的植被覆蓋度與平均氣溫和降水量的相關(guān)性通過(guò)顯著性檢驗(yàn)的區(qū)域分布差異較大，大部分區(qū)域的植被覆蓋度與同月的平均氣溫和降水量均表現(xiàn)出正相關(guān)，其中：與同月平均氣溫顯著正相關(guān)的區(qū)域面積為9 501 km2，約占11.44%，集中在研究區(qū)南部地區(qū);而呈現(xiàn)負(fù)相關(guān)的區(qū)域面積為11 315 km2，約占1363%，分布在研究區(qū)中部和北部地區(qū);呈顯著負(fù)相關(guān)區(qū)域面積極少，僅為0.02%.與同月降水量顯著正相關(guān)區(qū)域面積為13 657 km2，約占16.45%，集中在研究區(qū)北部和中部的西北地區(qū);呈現(xiàn)負(fù)相關(guān)的區(qū)域面積7 170 km2，約占8.64%，分布于南部地區(qū);呈顯著負(fù)相關(guān)區(qū)域面積極少，不足1 km2.

而生長(zhǎng)季的植被覆蓋度對(duì)前一月的平均氣溫與降水量的響應(yīng)結(jié)果則呈現(xiàn)出較大差異，其中：與前一月平均氣溫的相關(guān)性非常顯著，顯著正相關(guān)區(qū)域主要分布于中部和北部以及南部部分地區(qū)，面積達(dá)47 281 km2，約占56.95%;呈現(xiàn)正相關(guān)區(qū)域面積35 301 km2，約占42.52%，兩者所占比例達(dá)99.47%.與前一月降水量的相關(guān)性不顯著，呈正相關(guān)區(qū)域面積與負(fù)相關(guān)區(qū)域面積相近，空間分布上，呈正相關(guān)區(qū)域集中在研究區(qū)的中部和北部，而南部地區(qū)則呈負(fù)相關(guān)。

3.5討論

退耕還林還草工程實(shí)施后期，人類活動(dòng)減弱背景下，陜北黃土高原植被覆蓋度表現(xiàn)出平緩增長(zhǎng)趨勢(shì)，且與同月平均氣溫、降水和前一月平均氣溫的相關(guān)性明顯，表明退耕還林還草工程已取得顯著成效，陜北黃土高原地區(qū)的生態(tài)環(huán)境呈一定的自我恢復(fù)趨勢(shì)。

此外，陜北黃土高原的南部地區(qū)植被覆蓋度最高，但植被覆蓋度與年降水量、月降水量均呈現(xiàn)出一定的負(fù)相關(guān)，植被生長(zhǎng)對(duì)降水量存在負(fù)向反饋;在研究區(qū)中部和北部地區(qū)的植被覆蓋度較低，卻與降水量呈正相關(guān)，但部分地區(qū)的植被生長(zhǎng)對(duì)年平均氣溫、同月氣溫產(chǎn)生負(fù)向反饋。這種現(xiàn)象的出現(xiàn)，可能與研究區(qū)的植被分布類型有關(guān)，考慮到陜北黃土高原由北向南依次為草原區(qū)、森林草原區(qū)以及落葉闊葉林區(qū)，南部降水量較北部豐富，落葉闊葉林受到較多降水之后，水分不能及時(shí)蒸散反而影響植被生長(zhǎng)，而北部和中部地區(qū)降水量較少，植被類型以草地為主，氣溫較高時(shí)水分過(guò)快蒸散也會(huì)影響植被的生長(zhǎng)。

4結(jié)論

1）陜北黃土高原植被覆蓋度空間分布差異明顯，由北向南植被覆蓋度由低到高過(guò)渡明顯，2009—2016年植被覆蓋等級(jí)的空間轉(zhuǎn)移變化主要發(fā)生在北部和中部地區(qū)，表現(xiàn)出由中低植被覆蓋向中等植被覆蓋轉(zhuǎn)移、中等植被覆蓋向中高植被覆蓋轉(zhuǎn)移的態(tài)勢(shì);

2）2009—2016年陜北黃土高原植被覆蓋整體呈平緩增加趨勢(shì)，中部和南部部分地區(qū)植被覆蓋度呈減少趨勢(shì)，但變化趨勢(shì)皆不顯著，且年際變化率低，89.28%的區(qū)域處于-0.02～0.02 a-1之間;月際尺度上，研究區(qū)植被覆蓋度周期性變化明顯且呈逐漸增加的趨勢(shì);

3）研究區(qū)內(nèi)的年植被覆蓋度與年降水量、年平均氣溫的相關(guān)系數(shù)、偏相關(guān)系數(shù)均較小，相關(guān)性不顯著。不同區(qū)域的植被覆蓋度對(duì)氣候因子的響應(yīng)結(jié)果不同，陜北黃土高原植被覆蓋度與年平均氣溫和年降水量的相關(guān)性在空間上正負(fù)共存，其中年植被覆蓋度與年降水量以正相關(guān)為主，與年平均氣溫以負(fù)相關(guān)為主，研究區(qū)北部受年降水量影響較大，而中部受年平均氣溫影響較大;

4）植被覆蓋度對(duì)氣候因子的響應(yīng)存在滯后性。研究區(qū)植被覆蓋度受同月降水量的影響大于前一月的降水量，并且前一月平均氣溫對(duì)植被覆蓋度的影響遠(yuǎn)大于同月平均氣溫，生長(zhǎng)季的植被生長(zhǎng)受同月降水量和前一月平均氣溫的影響較大，且絕大部分地區(qū)均表現(xiàn)為正向相關(guān)。

參考文獻(xiàn)（References）：

[1]Hilker T，Lyapustin A I，Hall F G，et al.On the measurability of change in Amazon vegetation from MODIS[J].Remote Sensing of Environment，2015，166：233-242.

[2]Ochoa P A，F(xiàn)ries A，Mejía D，et al.Effects of climate，land cover and topography on soil erosion risk in a semiarid basin of the Andes[J].Catena，2016，140：31-42.

[3]甘春英，王兮之，李保生，等.連江流域近18年來(lái)植被覆蓋度變化分析[J].地理科學(xué)，2011，31（8）：1019-1024.

GAN Chunying，WANG Xizhi，LI Baosheng，et al.Changes of vegetation coverage during recent 18 years in Lianjiang river watershed[J].Scientia Geographica Sinica，2011，31（8）：1019-1024.

[4]Liu X，Zhu X，Pan Y，et al.Vegetation dynamics in QinlingDaba mountains in relation to climate factors between 2000 and 2014[J].Journal of Geographical Sciences，2016，26（1）：45-58.

[5]Yang Q，Qin Z，Li W，et al.Temporal and spatial variations of vegetation cover in Hulun buir grassland of inner Mongolia，China[J].Arid Soil Research & Rehabilitation，2012，26（4）：328-343.

[6]張東海，任志遠(yuǎn)，王曉峰，等.基于MODIS的陜西黃土高原植被覆蓋度變化特征及其驅(qū)動(dòng)分析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2013，29（1）：29-35.

ZHANG Donghai，REN Zhiyuan，WANG Xiaofeng，et al.MODIS based analysis of dynamic variation of fractional vegetation coverage of the Loess Plateau of Shaanxi and its driving forces[J].Journal of Ecology & Rural Environment，2013，29（1）：29-35.

[7]宋敏敏，張青峰，吳發(fā)啟，等.基于NDVI的延河流域時(shí)空演變分析[J].水土保持研究，2017，24（4）：6-11.

SONG Minmin，ZHANG Qingfeng，WU Faqi，et al.Spatiotemporal variation of NDVI in Yanhe basin[J].Research of Soil & Water Conservation，2017，24（4）：6-11.

[8]戚德康，南穎，孫博楊，等.2000年以來(lái)朝鮮植被覆蓋變化及其與氣候因子的相關(guān)性研究[J].干旱區(qū)地理，2017，40（4）：839-849.

QI Dekang，NAN Ying，SUN Boyang，et al.Changes of vegetation coverage and its relationship with climate factors in North Korea since 2000[J].Arid Land Geography，2017，40（4）：839-849.

[9]何寶忠，丁建麗，張喆，等.新疆植被覆蓋度趨勢(shì)演變實(shí)驗(yàn)性分析[J].地理學(xué)報(bào)，2016，71（11）：1948-1966.

HE Baozhong，DING Jianli，ZHANG Zhe，et al.Experimental analysis of spatial and temporal dynamics of fractional vegetation cover in Xinjiang[J].Acta Geographica Sinica，2016，71（11）：1948-1966.

[10]Zhou W，Gang C，Chen Y，et al.Grassland coverage interannual variation and its coupling relation with hydrothermal factors in China during 1982—2010[J].Journal of Geographical Sciences，2014，24（4）：593-611.

[11]阿多，趙文吉，宮兆寧，等.1981—2013華北平原氣候時(shí)空變化及其對(duì)植被覆蓋度的影響[J].生態(tài)學(xué)報(bào)，2017，37（2）：576-592.

A Duo，ZHAO Wenji，GONG Zhaoning，et al.Temporal analysis of climate change and its relationship with vegetation cover on the north China plain from 1981 to 2013[J].Acta Ecologica Sinica，2017，37（2）：576-592.

[12]Sun W，Song X，Mu X，et al.Spatiotemporal vegetation cover variations associated with climate change and ecological restoration in the Loess Plateau[J].Agricultural & Forest Meteorology，2015，209-210（1）：87-99.

[13]易浪，任志遠(yuǎn)，張翀，等.黃土高原植被覆蓋變化與氣候和人類活動(dòng)的關(guān)系[J].資源科學(xué)，2014，36（1）：166-174.

YI Lang，REN Zhiyuan，ZHANG Chong，et al.Vegetation cover，climate and human activities on the Loess Plateau[J].Resources Science，2014，36（1）：166-174.

[14]Xie B，Jia X，Qin Z，et al.Vegetation dynamics and climate change on the Loess Plateau，China： 1982—2011[J].Regional Environmental Change，2016，16（6）：1-12.

[15]劉憲鋒，楊勇，任志遠(yuǎn)，等.2000—2009年黃土高原地區(qū)植被覆蓋度時(shí)空變化[J].中國(guó)沙漠，2013，33（4）：1244-1249.

LIU Xianfeng，YANG Yong，REN Zhiyuan，et al.Changes of vegetation coverage in the Loess Plateau in 2000—2009[J].Journal of Desert Research，2013，33（4）：1244-1249.

[16]范建忠，李登科，周輝.陜西省退耕還林工程區(qū)植被覆蓋度的變化分析[J].干旱地區(qū)農(nóng)業(yè)研究，2013，31（4）：207-213.

FAN Jianzhong，LI DengKe，ZHOU Hui.Variation of fractional vegetation coverage in returning cropland to woodland project zones in Shaanxi Province[J].Agricultural Research in the Arid Areas，2013，31（4）：207-213.

[17]韓可欣，袁航.陜北黃土高原植被覆蓋時(shí)空動(dòng)態(tài)研究[J].貴陽(yáng)學(xué)院學(xué)報(bào)，2016，11（4）：24-30.

HAN Kexin，YUAN Hang.Spatial and temporal dynamic of vegetation cover in the Loess Plateau of northern Shaanxi[J].Journal of Guiyang University，2016，11（4）：24-30.

[18]肖強(qiáng)，陶建平，肖洋.黃土高原近10年植被覆蓋的動(dòng)態(tài)變化及驅(qū)動(dòng)力[J].生態(tài)學(xué)報(bào)，2016，36（23）：7594-7602.

XIAO Qiang，Tao Jianping，XIAO Yang.Dynamic vegetation cover change over the past 10 years on the Loess Plateau，China[J].Acta Ecologica Sinica，2016，36（23）：7594-7602.

[19]楊延征，趙鵬祥，郝紅科，等.基于SPOTVGT NDVI的陜北植被覆蓋時(shí)空變化[J].應(yīng)用生態(tài)學(xué)報(bào)，2012，23（7）：1897-1903.

YANG Yanzheng，ZHAO Pengxiang，HAO Hongke，et al.Spatiotemporal variation of vegetation in northern Shaanxi of northwest China based on SPOTVGT NDVI[J].Chinese Journal of Applied Ecology，2012，23（7）：1897-1903.

[20]李苗苗，吳炳方，顏長(zhǎng)珍，等.密云水庫(kù)上游植被覆蓋度的遙感估算[J].資源科學(xué)，2004，26（4）：153-159.

LI Miaomiao，WU Bingfang，YAN Changzhen，et al.Estimation of vegetation fraction in the upper basin of Miyun reservoir by remote sensing[J].Resources Science，2004，26（4）：153-159.

[21]穆少杰，李建龍，陳奕兆，等.2001—2010年內(nèi)蒙古植被覆蓋度時(shí)空變化特征[J].地理學(xué)報(bào)，2012，67（9）：1255-1268.

MU Shaojie，LI Jianlong，CHEN Yizhao，et al.Spatial differences of variations of vegetation coverage in Inner Mongolia during 2001—2010[J].Acta Geographica Sinica，2012，67（9）：1255-1268.

[22]張戈麗，徐興良，周才平，等.近30年來(lái)呼倫貝爾地區(qū)草地植被變化對(duì)氣候變化的響應(yīng)[J].地理學(xué)報(bào)，2011，66（1）：47-58.

ZHANG Geli，XU Xingliang，ZHOU Caiping，et al.Responses of vegetation changes to climatic variations in Hulun Buir grassland in past 30 years[J].Acta Geographica Sinica，2011，66（1）：47-58.

[23]王濤，楊梅煥.榆林地區(qū)植被指數(shù)動(dòng)態(tài)變化及其對(duì)氣候和人類活動(dòng)的響應(yīng)[J].干旱區(qū)研究，2017，34（5）：1133-1140.

WANG Tao，YANG Meihuan.Dynamic change of NDVI and its response to climate change and human activities in Yulin，Shaanxi Province，China[J].Arid Zone Research，2017，34（5）：1133-1140.

[24]陳效逑，王恒.1982—2003年內(nèi)蒙古植被帶和植被覆蓋度的時(shí)空變化[J].地理學(xué)報(bào)，2009，64（1）：84-94.

CHEN Xiaoqiu，WANG Heng.Spatial and temporal variations of vegetation belts and vegetation cover degrees in Inner Mongolia from 1982 to 2003[J].Acta Geographica Sinica，2009，64（1）：84-94.

[25]郭敏杰，張亭亭，張建軍，等.1982—2006年黃土高原地區(qū)植被覆蓋度對(duì)氣候變化的響應(yīng)[J].水土保持研究，2014，21（5）： 35-48.

GUO Minjie，ZHANG Tingting，ZHANG Jianjun，et al.Response of vegetation coverage to climate change in the loess plateau in 1982—2006[J].Research of Soil & Water Conservation，2014，21（5）： 35-48.

[26]宋富強(qiáng)，康慕誼，陳雅如，等.陜北黃土高原植被凈初級(jí)生產(chǎn)力的估算[J].生態(tài)學(xué)雜志，2009，28（11）：2311-2318.

SONG Fuqiang，KANG Muyi，CHEN Yaru，et al.Estimation of vegetation net primary productivity on North Shaanxi loess plateau[J].Chinese Journal of Ecology，2009，28（11）：2311-2318.