王傳勝, 李秋秋

(1.中國(guó)科學(xué)院 區(qū)域可持續(xù)發(fā)展分析與模擬重點(diǎn)實(shí)驗(yàn)室, 北京 100101;2.中國(guó)科學(xué)院大學(xué), 北京 100049; 3.中國(guó)科學(xué)院 地理科學(xué)與資源研究所, 北京 100101)

?

2000年以來(lái)西秦嶺地區(qū)植被覆蓋變化特征
——以隴南市為例

王傳勝1,3, 李秋秋1,2,3

(1.中國(guó)科學(xué)院 區(qū)域可持續(xù)發(fā)展分析與模擬重點(diǎn)實(shí)驗(yàn)室, 北京 100101;2.中國(guó)科學(xué)院大學(xué), 北京 100049; 3.中國(guó)科學(xué)院 地理科學(xué)與資源研究所, 北京 100101)

利用NDVI數(shù)據(jù)，以隴南市為案例區(qū)，研究了2000年以來(lái)西秦嶺地區(qū)植被覆蓋度變化的時(shí)空特征。采集2000—2010年6—10月NDVI數(shù)據(jù)，利用像元二分法設(shè)計(jì)了植被覆蓋度(Vc)的計(jì)算方法，結(jié)果顯示，2000—2010年隴南市植被覆蓋度變化呈基本穩(wěn)定態(tài)勢(shì)，與同期溫度和降雨的變化顯著相關(guān)。利用一元線性回歸和標(biāo)準(zhǔn)差方法，分析了Vc的時(shí)空變化特征，結(jié)果顯示，土地面積將近90%的區(qū)域表現(xiàn)為穩(wěn)定趨勢(shì)，Vc的年際波動(dòng)較小，空間分布與林地的分布區(qū)域基本吻合；將近1/10的區(qū)域表現(xiàn)為上升趨勢(shì)，Vc的年際波動(dòng)較大，主要分布在徽成盆地、西禮盆地、白龍江河谷陽(yáng)面山坡等區(qū)域，與耕地、草地分布區(qū)域基本吻合，表明退耕還林工程對(duì)植被覆蓋度增加的貢獻(xiàn)顯著。論文為20世紀(jì)末以來(lái)西秦嶺地區(qū)生態(tài)建設(shè)的成效提供了證據(jù)，也為當(dāng)?shù)厣罨鷳B(tài)環(huán)境建設(shè)和可持續(xù)發(fā)展提供了借鑒。

植被覆蓋度; 時(shí)空變化特征; 西秦嶺

植被覆蓋和植被活動(dòng)及其變化長(zhǎng)期以來(lái)一直是地球陸地生態(tài)系統(tǒng)研究的主體內(nèi)容。近年來(lái)，在全球變化領(lǐng)域，通過(guò)了解大氣組成、土地利用和氣候變化對(duì)植被的影響與反饋，探討植被對(duì)全球變化的響應(yīng)[1-4]；另一方面，在可持續(xù)發(fā)展領(lǐng)域，通過(guò)了解區(qū)域植被覆蓋度和植被活動(dòng)的變化，分析人類活動(dòng)對(duì)植被生態(tài)系統(tǒng)的影響[5]。近20 a來(lái)，利用衛(wèi)星遙感數(shù)據(jù)紅光和紅外波段的不同組合進(jìn)行植被研究，取得了較好的效果[6-7]。在眾多的植被指數(shù)中，NDVI(Normalized Differential Vegetation Index，歸一化植被指數(shù))是應(yīng)用最為廣泛的一種，已有大量的文獻(xiàn)研究了我國(guó)不同區(qū)域植被變化的時(shí)空特征及其與氣候、人類活動(dòng)的關(guān)系[8-11]，研究發(fā)現(xiàn)NDVI與降水、氣溫等自然氣候因子[12]和城鎮(zhèn)發(fā)展、退耕還林等人類活動(dòng)[13-14]有較好相關(guān)性。在現(xiàn)有研究中，崔曉臨等[15]利用2000—2009年NDVI數(shù)據(jù)，研究了秦嶺東部地區(qū)(陜西、河南)植被指數(shù)的變化，結(jié)論認(rèn)為植被覆蓋呈逐年增加態(tài)勢(shì)，同時(shí)又進(jìn)一步研究了海拔梯度差異的氣溫響應(yīng)，認(rèn)為高海拔地區(qū)植被生態(tài)系統(tǒng)更易受到全球氣候變化的影響。還有其他研究涉及到秦嶺部分地區(qū)，結(jié)果認(rèn)為隴南山區(qū)植被呈增加趨勢(shì)[16]。

本文重點(diǎn)關(guān)注秦嶺西部地區(qū)，即嘉陵江干流以西的西秦嶺部分，在行政上屬甘肅省，包括隴南市、天水和定西2市的渭河以南部分縣，以及甘南州的迭部、舟曲和臨潭3縣，因西秦嶺主要部分為隴南市，故本研究以隴南市為案例區(qū)域。隴南市為全國(guó)扶貧重點(diǎn)地區(qū)，同時(shí)也是全國(guó)水土流失防治和生態(tài)建設(shè)的重點(diǎn)地區(qū)。據(jù)隴南市林業(yè)局提供的資料，自1999年以來(lái)，通過(guò)實(shí)施退耕還林、天然林保護(hù)、長(zhǎng)防林、生態(tài)公益林、自然保護(hù)區(qū)等生態(tài)工程，累計(jì)造林45.33萬(wàn)hm2，其中累計(jì)退耕面積8.24萬(wàn)hm2，全市森林覆蓋率由38.9%增加到42%。本研究通過(guò)植被覆蓋度變化的時(shí)空特征揭示，一方面為驗(yàn)證20世紀(jì)末以來(lái)生態(tài)建設(shè)的成效提供依據(jù)，另一方面也為當(dāng)?shù)厣罨鷳B(tài)環(huán)境建設(shè)和可持續(xù)發(fā)展提供借鑒。

1　研究區(qū)概況

隴南市經(jīng)緯度為32°35′45″—34°32′00″N，104°01′19″ —106°35′20″E，在地貌上，北部屬黃土高原丘陵溝壑區(qū)，西部為岷山山脈高山嶺谷區(qū)，東南部為秦巴山地高山峽谷區(qū)；在氣候上，為亞熱帶向暖溫帶過(guò)渡地區(qū)，年平均氣溫5～16℃，年降雨量400～1 000 mm；在流域單元上，是甘肅省唯一的長(zhǎng)江流域地區(qū)；在植被上，為常綠—落葉針闊混交林區(qū)域。

隴南市與陜、川兩省接壤，轄1區(qū)8縣，土地總面積2.79萬(wàn)km2，2012年常住人口256.95萬(wàn)人，將近占全省比重的1/10。境內(nèi)山地廣布，氣候垂直分異顯著，海拔高度介于550～4 187 m之間，坡度25°以上面積占國(guó)土面積的43%[17]，為全國(guó)水土流失和泥石流災(zāi)害的重點(diǎn)防治區(qū)，也是全國(guó)退耕還林的重點(diǎn)實(shí)施區(qū)域。據(jù)隴南市土地局資料，全市將近2/3的土地利用類型為林草地，其中，林地面積占全市面積的49.67%，牧草地面積占14.98%，耕地和園地面積占20.48%，建設(shè)用地面積占1.66%，水域面積占1.05%，其他用地面積占12.15%。耕地主要分布在東部的徽(縣)成(縣)盆地和北部的西(和)禮(縣)盆地，以及武都區(qū)的白龍江干流河谷和文縣的白水江干流河谷區(qū)域。

2　數(shù)據(jù)來(lái)源與研究方法

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

植被覆蓋度變化分析主要使用NDVI數(shù)據(jù)，數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站的MODIS中國(guó)合成產(chǎn)品，提取基于Terra衛(wèi)星的2000—2010年NDVI中國(guó)1 km格網(wǎng)的NDVI月合成產(chǎn)品，該產(chǎn)品是取每月的最大值為月度實(shí)際值，坐標(biāo)系為WGS84?？紤]到夏秋季是植被活動(dòng)的旺盛時(shí)期，因此選擇6—10月份的NDVI值作為背景數(shù)據(jù)，在ARCGIS軟件支持下取得6—10月的NDVI最大值，作為NDVI的年度考察值。

輔助分析數(shù)據(jù)包括：(1) 地形圖和土地利用柵格數(shù)據(jù)，前者來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站提供的ASTER GDEM DEM數(shù)據(jù)，后者來(lái)源于中國(guó)科學(xué)院提供的2000年遙感數(shù)據(jù)，兩項(xiàng)數(shù)據(jù)主要用于研究區(qū)域一般特征分析和植被覆蓋度變化的空間特征分析；(2) 氣象數(shù)據(jù)，來(lái)源于國(guó)家氣象局提供的武都測(cè)站2000—2010年度氣溫和降雨數(shù)據(jù)，主要用于和植被覆蓋度變化進(jìn)行相關(guān)分析；(3) 統(tǒng)計(jì)數(shù)據(jù)，主要來(lái)源于實(shí)地調(diào)研的專業(yè)部門數(shù)據(jù)和相關(guān)年度的統(tǒng)計(jì)年鑒數(shù)據(jù)。

2.2研究方法

2.2.1植被覆蓋度植被覆蓋度的計(jì)算通常采用像元二分法[18-20]，這是一種簡(jiǎn)單實(shí)用的遙感估算模型，其原理是假設(shè)一個(gè)影像像元的地表覆蓋分為有植被和無(wú)植被兩部分，遙感傳感器觀測(cè)到的光譜信息也是這兩部分要素的線性加權(quán)組合，權(quán)重是兩要素的面積比重，植被要素的面積比重就是像元中植被部分的權(quán)重，因此，植被面積比重的計(jì)算公式為[21-22]：Pveg=(S-Ssoil)/(Sveg-Ssoil)，其中，Pveg為觀測(cè)區(qū)域植被面積的比重，S為混合像元的遙感信息值，Sveg是完全植被覆蓋像元的遙感信息值，Ssoil為完全無(wú)植被覆蓋像元的遙感信息值。

本文以2000—2010年NDVI的年度值為遙感信息值，定義遙感估算的植被覆蓋面積比重為植被覆蓋度Vc(Vegetation Coverage)，則根據(jù)上述原理，Vc的計(jì)算公式為：

Vc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

(1)

式中：NDVI——像元的實(shí)際值；NDVIveg——完全植被覆蓋像元的遙感信息值；NDVIsoil——完全無(wú)植被覆蓋像元的遙感信息值。理論上，NDVIveg——純植被像元；NDVIsoil——純土壤像元，通常被分別取值為影像象元中的最大值與最小值；但實(shí)際上，因地表不同地物之間的相互影響，圖像不可避免地存在噪聲，因而實(shí)際研究中并不是直接取影像中的最大值和最小值，而是根據(jù)影像像元值的累積頻率分布情況，設(shè)定一定的置信度區(qū)間，取置信度區(qū)間的最大值和最小值。一般情況下，如沒(méi)有足夠的實(shí)地觀測(cè)的樣點(diǎn)數(shù)據(jù)，置信區(qū)間通常設(shè)定為5%～95%[23]。因此，本文以NDVI像元值累積分布5%的值為影像中NDVI最小值，低于此值的Vc值均取0；NDVI像元值累積分布95%的為像元中最大值，高于此值的Vc值均取1；介于其中的用公式(1)計(jì)算結(jié)果，利用ArcGIS的reclassify模塊和raster calculator模塊取得最終結(jié)果。

2.2.2一元線性回歸一元線性回歸方法主要用來(lái)測(cè)度Vc隨時(shí)間變化的趨勢(shì)[24-26]，通過(guò)單個(gè)像元Vc隨時(shí)間的變化趨勢(shì)模擬刻畫整個(gè)區(qū)域植被覆蓋度變化的空間態(tài)勢(shì)，回歸方程的斜率是Vc隨時(shí)間變化的主要參數(shù)，計(jì)算公式為：

(2)

式中：Vcslp——方程的斜率；i，n——Vc自2000年起的年度序號(hào)，共11 a，i∈(1，11)，n=11。Vcslp>0，植被覆蓋度呈上升趨勢(shì)，值越大，表明植被活動(dòng)強(qiáng)度越大；Vcslp<0，植被覆蓋度呈下降趨勢(shì)，值越小，表明植被活動(dòng)強(qiáng)度下降的趨勢(shì)越顯著。

2.2.3標(biāo)準(zhǔn)差標(biāo)準(zhǔn)差方法主要用來(lái)測(cè)度各年度Vc像素值的離散程度或變化幅度，本文也用來(lái)和植被覆蓋度斜率變化進(jìn)行時(shí)空組合分析，前者表示植被覆蓋度變化的幅度，而后者表征變化方向。本研究中，標(biāo)準(zhǔn)差的計(jì)算使用ArcGIS空間分析工具中的cell statistics模塊，選擇2000—2010年各年Vc的標(biāo)準(zhǔn)差輸出形式得到具體的結(jié)果。

3　結(jié)果與分析

3.1全域植被覆蓋度變化的時(shí)序特征

利用像元二分法計(jì)算2000—2010年隴南市各年的Vc，分析Vc隨時(shí)間變化的特點(diǎn)。圖1顯示，2000年以來(lái)，隴南市Vc趨勢(shì)線較為平穩(wěn)，植被覆蓋度圍繞0.64左右上下波動(dòng)，近年來(lái)趨于穩(wěn)定。

3.2植被覆蓋度變化趨勢(shì)的空間分異

3.2.1變化趨勢(shì)劃分因?yàn)橹脖坏淖兓嬖谧陨碜匀蛔兓拈撝捣秶虼藢?duì)于Vcslp是上升、穩(wěn)定還是下降的判斷，首先需要確定其自然變動(dòng)的閾值區(qū)。目前基本采用兩種方法，一種是選擇植被覆蓋度為0的訓(xùn)練區(qū)，看Vc變化的值域區(qū)間，根據(jù)國(guó)內(nèi)學(xué)者在青藏高原的研究[27]，值域-0.2～0.2為植被自然變動(dòng)的范圍；另一種方法借助統(tǒng)計(jì)學(xué)方法，如采用Jenks′natural breaks斷點(diǎn)分級(jí)，確定樣本的自然斷點(diǎn)[24]。本文采用前一種方法，因隴南市地處河流上游地區(qū)，河流為線狀要素，也無(wú)面積較大的湖泊，無(wú)法滿足訓(xùn)練區(qū)的精度要求，考慮到隴南市城鎮(zhèn)建設(shè)用地區(qū)域原本狹小，近幾年在城鎮(zhèn)擴(kuò)張中綠化覆蓋率提高不大，在土地利用圖中面積最大的可達(dá)2 km2，故選擇城鎮(zhèn)建設(shè)用地作為訓(xùn)練區(qū)。根據(jù)對(duì)2000年土地利用數(shù)據(jù)中建設(shè)用地Vc斜率柵格數(shù)據(jù)分布的特征值分析，Vc介于-0.027 1～0.018 61，和第一種方法中確定的值域相當(dāng)，因此將此范圍作為Vc變化相對(duì)穩(wěn)定的閾值區(qū)間。

表1顯示，近十年來(lái)將近90%的區(qū)域Vc變化趨勢(shì)相對(duì)穩(wěn)定，將近1/10的區(qū)域表現(xiàn)為上升趨勢(shì)，面積2 563.23 km2，1.33%的區(qū)域?yàn)榻档挖厔?shì)，面積372.08 km2。圖2顯示，植被覆蓋度增加趨勢(shì)區(qū)域主要分布在中北部黃土高原地貌類型區(qū)，大體上呈面狀或片狀集中分布于徽成盆地區(qū)、西禮盆地區(qū)、武都區(qū)白龍江干流河谷的陽(yáng)面山坡區(qū)，這些區(qū)域的土地利用類型以耕地和草地為主，且調(diào)研資料顯示，這些區(qū)域也是退耕還林面積較大的區(qū)域，在一定程度上可以證實(shí)退耕還林工程對(duì)植被覆蓋度增加的貢獻(xiàn)顯著。降低趨勢(shì)的區(qū)域一是南部文縣的河谷區(qū)域，尤其是白水江上游；二是西和、宕昌兩縣之間的山區(qū)。

圖1　隴南市Vc時(shí)序變化趨勢(shì)

變化趨勢(shì)閾值范圍像元數(shù)實(shí)際面積/km2面積比重/%穩(wěn)定趨勢(shì)-0.0271～0.018616951124964.6989.48上升趨勢(shì)0.01861～0.084671372563.239.19下降趨勢(shì)-0.0989～-0.02711036372.081.33

3.2.2變化趨勢(shì)的年度波動(dòng)為便于和植被覆蓋度變化的斜率進(jìn)行比較分析，采用Jenks′natural breaks方法將Vc標(biāo)準(zhǔn)差分為三級(jí)，結(jié)果見(jiàn)表2，圖3。表2顯示，標(biāo)準(zhǔn)差最大值為0.371 7，標(biāo)準(zhǔn)差較大的區(qū)域，即介于0.141 4-0.371 7的區(qū)域柵格單元數(shù)有10 646個(gè)，面積占13.70%，標(biāo)準(zhǔn)差較小的0-0.080 2的區(qū)域面積接近全部區(qū)域的一半，說(shuō)明大多數(shù)區(qū)域植被覆蓋度變化范圍不大。圖3顯示，標(biāo)準(zhǔn)差的空間分布態(tài)勢(shì)和斜率的空間分布態(tài)勢(shì)類似，標(biāo)準(zhǔn)差較大的區(qū)域主要分布在斜率變化較大(上升或下降)的區(qū)域，如徽成盆地區(qū)、西禮盆地區(qū)、白龍江和白水江兩岸，且主要為耕地和草地的分布區(qū)域。

圖2　Vc上升和下降趨勢(shì)的空間分布

標(biāo)準(zhǔn)差閾值范圍柵格數(shù)面積/km2面積比重/%0～0.08023782313584.0348.690.0802～0.14142921510492.4937.610.1414～0.3717106463823.4813.70

圖3　Vc變化的標(biāo)準(zhǔn)差分布

3.2.3變化趨勢(shì)的時(shí)空組合分析按照表3的分類結(jié)果，把標(biāo)準(zhǔn)差值域進(jìn)一步分為兩類，利用ArcGIS的柵格計(jì)算模塊對(duì)植被覆蓋度變化進(jìn)行斜率和標(biāo)準(zhǔn)差的組合分類，結(jié)果分為6類區(qū)域。標(biāo)準(zhǔn)差小則表明像元的Vc年際變化較小，植被覆蓋度處于平穩(wěn)的改善、穩(wěn)定或退化趨勢(shì)；標(biāo)準(zhǔn)差大即像元的Vc年際波動(dòng)較大，植被覆蓋度處于急劇或大幅波動(dòng)的改善、穩(wěn)定或退化趨勢(shì)。表3顯示，研究區(qū)大體處于Vc小范圍變化的穩(wěn)定態(tài)勢(shì)中；植被覆蓋度改善且變化范圍大的地區(qū)可視為生態(tài)效益明顯或起伏中轉(zhuǎn)好，此類地區(qū)占研究區(qū)面積5.86%，只有少數(shù)地區(qū)處于植被覆蓋度在大范圍波動(dòng)中退化的狀態(tài)。附圖11顯示，斜率穩(wěn)定趨勢(shì)地區(qū)中，標(biāo)準(zhǔn)差小的區(qū)域多呈面狀連續(xù)分布態(tài)勢(shì)，標(biāo)準(zhǔn)差大的區(qū)域多呈零星的點(diǎn)狀或條帶狀分布；上升趨勢(shì)地區(qū)中則相反，標(biāo)準(zhǔn)差大的區(qū)域多呈面狀、連續(xù)性分布，標(biāo)準(zhǔn)差小的區(qū)域多為零星的點(diǎn)狀分布。

表3　Vc變化的空間分類

3.3植被覆蓋度變化的影響因子分析

2000—2010年，隴南市降雨量和溫度年際波動(dòng)較大。從Vc與氣候要素特征值變化的關(guān)系來(lái)看，隴南市植被覆蓋度與同期平均氣溫呈顯著負(fù)相關(guān)，而和降雨量呈顯著正相關(guān)，氣溫、降水的變化對(duì)Vc變化的貢獻(xiàn)達(dá)到54.7%，43.5%。

有大量研究論證，植被覆蓋既受降水氣溫等自然生長(zhǎng)環(huán)境的影響，同時(shí)也受到人為活動(dòng)的作用，如退耕還林、天然林保護(hù)等生態(tài)環(huán)境建設(shè)[28-31]。在上文植被覆蓋改善區(qū)域中，徽成盆地區(qū)、西禮盆地區(qū)、武都區(qū)白龍江干流河谷的陽(yáng)面山坡區(qū)也是隴南市退耕還林工程實(shí)施較多的區(qū)域，故而隴南市植被蓋度變化的影響因子可以歸納為氣候因子變化與退耕還林的影響。

4　結(jié) 論

(1) 2000年以來(lái)，隴南市植被覆蓋度變化總體相對(duì)穩(wěn)定。時(shí)間序列的分析表明，NDVI總體呈上升態(tài)勢(shì)，Vc表現(xiàn)相對(duì)平穩(wěn)。植被覆蓋度變化和6—8月平均氣溫、降雨量有較高相關(guān)性，表明隴南市氣溫、降雨變化對(duì)植被覆蓋度影響較大。

(2) 植被覆蓋度空間變化態(tài)勢(shì)和時(shí)間序列的分析結(jié)果基本一致，將近90%的區(qū)域表現(xiàn)為穩(wěn)定趨勢(shì)，將近1/10的區(qū)域表現(xiàn)為上升趨勢(shì)。從空間變化態(tài)勢(shì)上看，上升趨勢(shì)區(qū)域多數(shù)呈面狀或片狀的集中分布態(tài)勢(shì)，下降趨勢(shì)區(qū)域多呈局部或零星的分散分布態(tài)勢(shì)；從年度變化的空間態(tài)勢(shì)來(lái)看，植被覆蓋度顯著上升或下降區(qū)域，Vc的年際波動(dòng)較大，空間分布與耕地、草地分布區(qū)域基本吻合；而植被覆蓋度變化相對(duì)穩(wěn)定區(qū)域，Vc的年際波動(dòng)較小，空間分布與林地的分布區(qū)域基本吻合。植被覆蓋度上升區(qū)域主要分布在徽成盆地、西禮盆地、白龍江河谷陽(yáng)面山坡等區(qū)域，也是退耕還林面積較大的區(qū)域，一定程度上說(shuō)明退耕還林工程對(duì)植被覆蓋度增加的貢獻(xiàn)顯著。

(3) NDVI對(duì)揭示研究區(qū)域的植被覆蓋度變化有較好的效果，但因研究區(qū)域植被的垂直分異顯著，植被隨氣候變化的時(shí)空差異較大，目前采用的數(shù)據(jù)，時(shí)間和空間的分辨率依然不夠，特別是大比例尺的植被數(shù)據(jù)和土地利用數(shù)據(jù)不足，增加了研究結(jié)果的不確定性，這在下一步的深化研究中需要補(bǔ)充。另外，對(duì)人類活動(dòng)的刻畫，特別是退耕還林以來(lái)植被覆蓋度增加的原因解釋，需要更詳盡的實(shí)地調(diào)研，進(jìn)一步豐富研究?jī)?nèi)容。

[1]樸世龍,方精云.1982—1999年我國(guó)陸地植被活動(dòng)對(duì)氣候變化響應(yīng)的季節(jié)差異[J].地理學(xué)報(bào),2003,58(1):119-125.

[2]Cramer W P, Leemans R. Assessing impacts of climate change on vegetation using climate classification systems[M].London: Chopman and Holl,1993.

[3]周廣勝,王玉輝.全球變化與氣候—植被分類研究和展望[J].科學(xué)通報(bào),1999,24(44):2587-2593.

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

[5]田慶久,閔祥軍.植被指數(shù)研究進(jìn)展.地球科學(xué)進(jìn)展[J],1998,13(4):327-333.

[6]賈坤,姚云軍,魏香琴,等.植被蓋度遙感估算研究進(jìn)展[J].地球科學(xué)進(jìn)展,2013,28(7):774-782.

[7]韋振鋒,王德光,張翀,等.近12年陜甘寧黃土高原區(qū)植被物候時(shí)空變化特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2014,30(4):423-429.

[8]張瑋,張鐿鋰,王兆鋒,等.珠穆朗瑪峰自然保護(hù)區(qū)植被變化分析[J].地理科學(xué)進(jìn)展,2006,25(3):12-22.

[9]韓秀珍,李三妹,羅敬寧,等.近20 a中國(guó)植被時(shí)空變化研究[J].干旱區(qū)研究,2008,25(6):753-759.

[10]孫艷玲,郭鵬,延曉冬,等.內(nèi)蒙古植被覆蓋變化及其與氣候、人類活動(dòng)的關(guān)系[J].自然資源科學(xué),2010,25(3):407-414.

[11]鄭有飛,劉宏舉,吳榮軍,等.貴州省NDVI變化及其與主要?dú)夂蛞蜃拥南嚓P(guān)性[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2009,25(1):12-17.

[12]徐興奎,陳紅, Levy Jason K.氣候變暖背景下青藏高原植被覆蓋特征的時(shí)空變化及其成因分析[J].科學(xué)通報(bào),2008,53(4):456-462.

[13]王靜,郭鈮,蔡迪花,等.瑪曲縣草地退牧還草工程效果評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2009,29(3):1276-1284.

[14]周洪建,王靜愛(ài),岳耀杰,等.人類活動(dòng)對(duì)植被退化/恢復(fù)影響的空間格局:以陜西省為例[J].生態(tài)學(xué)報(bào),2009,29(9):4847-4856.

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

[16]楊尚武,張勃,趙一飛,等.1998—2011年甘肅河?xùn)|地區(qū)NDVI時(shí)空變化[J].干旱區(qū)研究,2014,31(1):74-79.

[17]陳行,朱珊珊,李秋秋,等.西秦嶺地區(qū)農(nóng)戶生計(jì)改變對(duì)生態(tài)建設(shè)的影響[J].山地學(xué)報(bào),2014,32(6):662-670.

[18]Leprieur C, Verstraete M M, Pinty B. Evaluation of the performance of various vegetation indices to retrieve vegetation cover from AVHRR data[J]. Remote Sensing Reviews, 1994,10(4):265-284.

[19]陳晉,陳云浩,何春陽(yáng),等.基于土地覆蓋分類的植被覆蓋率估算亞像元模型與應(yīng)用[J].遙感學(xué)報(bào),2001,5(6):416-422.

[20]Zribi M, Le Hegarat-Mascle S, Taconet O, et al. Derivation of wild vegetation cover density in semi-arid regions: ERS2/SAR evaluation[J]. International Journal of Remote Sensing, 2003,24(6):1335-1352.

[21]劉廣峰,吳波,范文義,等.基于像元二分模型的沙漠化地區(qū)植被覆蓋度提取:以毛烏素沙地為例[J].水土保持研究,2007,14(2):268-271.

[22]郭芬芬,范建容,嚴(yán)冬,等.基于像元二分模型的昌都縣植被蓋度遙感估算[J].中國(guó)水土保持,2010(5):65-67.

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

[24]于伯華,呂昌河,呂婷婷,等.青藏高原植被覆蓋變化的地域分異特征[J].地理科學(xué)進(jìn)展,2009,28(3),391-397.

[25]梁四海,陳江,金曉梅,等.近21 a青藏高原植被覆蓋變化規(guī)律[J].地球科學(xué)進(jìn)展,2007,22(1):33-40.

[26]馬明國(guó),王建,王雪梅.基于遙感的植被年際變化及其與氣候關(guān)系研究進(jìn)展[J].遙感學(xué)報(bào),2006,10(3):421-431.

[27]Ding Mingjun, Zhang Yili, Shen Zhenxi, et al. Land cover change along the Qinghai Tibet Highway and Railway from 1981 to 2001[J]. Journal of Geographical Sciences, 2006,16(4):387-395.

[28]廖清飛,張?chǎng)?馬全,等.青海省東部農(nóng)業(yè)區(qū)植被覆蓋時(shí)空演變遙感監(jiān)測(cè)與分析[J].生態(tài)學(xué)報(bào),2014,34(20):5936-5943.

[29]孫智輝,曹雪梅,李新亞,等.氣候變化和人類活動(dòng)對(duì)吳起土壤侵蝕的影響[J].水土保持研究,2009,16(6):30-34,39.

[30]田海靜,曹春香,戴晟懋,等.準(zhǔn)格爾旗植被覆蓋度變化的時(shí)間序列遙感監(jiān)測(cè)[J].地球信息科學(xué)學(xué)報(bào),2014,16(1):126-133.

[31]高彥凈,謝余初,錢大文,等.甘肅白龍江流域植被覆蓋度及景觀格局變化[J].水土保持研究,2015,22(1):181-187.

Study on Spatiotemporal Variations of Vegetation Cover in the Western Qinling Mountain Areas During the Period from 2000 to 2010—A Case of Longnan of Gansu, China

WANG Chuansheng1,3, LI Qiuqiu1,2,3

(1.Key Laboratory of Regional Sustainable Development Modeling, CAS, Beijing 100101, China; 2.University of Chinese AcademyofSciences,Beijing100049,China; 3.InstituteofGeographicSciencesandNaturalResourcesResearch,CAS,Beijing100101,China)

As the main topic of terrestrial ecosystem study, vegetation cover change gained extensive attention in the process of global climate change and sustainable development recently. Based on MODIS NDVI data from June to October during the period from 2000 to 2010, taking Longnan City as a case area, we developed the calculation method of vegetation cover (Vc) by using Pixel Dichotomy model and analyzed the spatiotemporal variations of vegetation cover in the western Qinling Mountain areas by using simple linear regression and standard deviation method. The results showed that vegetation cover remained stable and significantly correlated with temperature and precipitation during the decade. The vegetation cover of 90% of study area showed the stability with small annual variation and also was consistent with the spatial distribution of forest land; while the vegetation cover in other 10% of study areas showed the growing tendency with significant variations and also was consistent with the spatial distribution of farmland and grassland, especially in Huicheng Basin, Xili Basin and adret slope of Bailong River Valley, which indicates that the Grain for Green project has the great contribution to vegetation cover increase. This study proved the effect of ecological construction in the western Qinling Mountain areas since late last century. All the findings can also provide

for local ecological environment construction and sustainable development.

vegetation cover; spatiotemporal variation; western Qinling Mountain

2015-02-08

2015-04-10

國(guó)家自然科學(xué)基金資助項(xiàng)目(41171109);中國(guó)科學(xué)院重點(diǎn)部署項(xiàng)目(KZZD-EW-06)

王傳勝(1965—),男,甘肅蘭州人,副研究員,博士,研究方向:區(qū)域規(guī)劃與資源環(huán)境保護(hù)。E-mail:wangcs@igsnrr.ac.cn

TP753; K921/927

A

1005-3409(2016)02-0308-05