藺星娜，牛健植，賈京偉，駱紫藤，杜 捷，2

植被作為聯(lián)系土壤、大氣和水分的自然紐帶［1］，對(duì)土壤侵蝕的防治、水資源的調(diào)節(jié)、水質(zhì)的改善、區(qū)域小氣候的改善及水土保持功能的發(fā)揮具有十分重要的作用;因此，對(duì)植被進(jìn)行動(dòng)態(tài)監(jiān)測(cè)，并分析其覆蓋度變化及影響因素，已經(jīng)成為全球變化研究的一個(gè)重要領(lǐng)域［2-5］。歸一化植被指數(shù)(Normal Difference Vegetation Index，NDVI)能夠廣泛應(yīng)用于植被覆蓋監(jiān)測(cè)、土地覆蓋分類(lèi)和作物估產(chǎn)等眾多研究領(lǐng)域［6］，與植被分布密度成線(xiàn)性相關(guān)，是評(píng)估地表植被覆蓋和生長(zhǎng)狀況的一個(gè)簡(jiǎn)單、有效的度量參數(shù)［7］。目前常用的 NDVI系列數(shù)據(jù)包括 NOAA/AVHRR數(shù)據(jù)(分辨率為8 km×8 km)、EOS/MODIS數(shù)據(jù)(分辨率為1 km×1 km)，以及SPOT Vegetation NDVI數(shù)據(jù)(分辨率為1 km ×1 km)［8］。與其他數(shù)據(jù)相比，SPOT Vegetation數(shù)據(jù)具有空間分辨率較高、剔除了強(qiáng)水汽吸收帶等優(yōu)勢(shì)，廣泛應(yīng)用于植被變化監(jiān)測(cè)。

植被覆蓋變化是氣候和人類(lèi)活動(dòng)共同作用的結(jié)果［9］。其中:降水和溫度通過(guò)影響植物光合作用、呼吸作用及土壤有機(jī)碳分解等影響植物生長(zhǎng)和分布，導(dǎo)致植被覆蓋度變化;人類(lèi)活動(dòng)通過(guò)改變土地利用狀況，直接改變生態(tài)系統(tǒng)的類(lèi)型、結(jié)構(gòu)和功能［10］以及植樹(shù)造林等活動(dòng)，從而改變植被覆蓋度。目前，國(guó)內(nèi)外采用遙感手段，動(dòng)態(tài)監(jiān)測(cè)植被變化，并大量研究植被變化與氣候因子及人類(lèi)活動(dòng)的關(guān)系，認(rèn)為植被覆蓋度與氣溫、降水有很強(qiáng)的相關(guān)關(guān)系，氣溫越高，降水量越大，植被覆蓋度越大，植被覆蓋度對(duì)氣溫與降水量的響應(yīng)有很強(qiáng)的季節(jié)性;同時(shí)發(fā)現(xiàn)，月植被覆蓋度對(duì)月平均降水量的響應(yīng)普遍滯后1個(gè)月［11-14］;土地利用類(lèi)型之間的轉(zhuǎn)化對(duì)不同等級(jí)植被覆蓋度的變化有很大影響，且植被覆蓋度的變化在不同研究地對(duì)土地利用的響應(yīng)程度不同［15-18］。還有學(xué)者認(rèn)為，土地利用的變化較氣候變化更容易影響植被覆蓋度［10］;還有研究者采用殘差法［19］、修正NDVI法［20］、主成分分析［21］、因子分析及降水利用率(RUE)［22］及 NDVI與 RUE 相結(jié)合［23］等方法，剔除氣候?qū)χ脖蛔兓挠绊?，進(jìn)而研究人類(lèi)活動(dòng)對(duì)植被變化的影響。其中，NDVI與RUE相結(jié)合的方法，可以剔除降雨量變化對(duì)植被覆蓋變化的影響，判斷人類(lèi)活動(dòng)對(duì)地表的干擾程度排［23］。

清水河流域位于河北省張家口市，主要分布于崇禮縣，作為2022年冬奧會(huì)的主辦地區(qū)之一，區(qū)域水土資源的有效利用、植物資源景觀生態(tài)功能的高效發(fā)揮，將會(huì)保證冬奧會(huì)成功舉辦;但該區(qū)域面臨土少石多，水資源匱乏，西溝區(qū)域土壤侵蝕嚴(yán)重，植被亟待恢復(fù)等重要問(wèn)題。因此，研究植被覆蓋度的時(shí)空變化特征及其影響因子，對(duì)當(dāng)?shù)刂脖唤ㄔO(shè)及水土保持工作的繼續(xù)開(kāi)展有著十分重要的意義。筆者通過(guò)遙感技術(shù)，研究清水河流域多年植被覆蓋度的時(shí)空動(dòng)態(tài)變化，結(jié)合NDVI與RUE的變化趨勢(shì)，分析其驅(qū)動(dòng)因子，為該地區(qū)植被恢復(fù)、生態(tài)系統(tǒng)穩(wěn)定、景觀功能完善提供科學(xué)依據(jù)和決策服務(wù)。

1 研究區(qū)概況

清水河上游流域位于 E 114°47'～115°31'，N 40°49'～41°16'，地處河北省張家口市大境門(mén)外的紅旗橋斷面以上，包括橋西區(qū)的部分區(qū)域小西溝和崇禮縣大部分區(qū)域，面積 2 095 km2［24］;屬于壩上壩下過(guò)渡型山區(qū)，海拔813～2 174 m，東西差異顯著;屬于東亞大陸性季風(fēng)氣候，中溫帶亞干旱區(qū)，干旱指數(shù)為2.5～3.0，年平均降水量488 mm。土壤類(lèi)型主要是棕壤土與褐土;其中，棕壤土主要分布在流域北半部，土質(zhì)較好，褐土主要分布在流域南部，沿河兩岸土壤主要為草甸土［25］。研究區(qū)位置如圖1。

圖1 研究區(qū)示意圖Fig．1 Location of Study Area

2 研究方法

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

1999—2015年，研究區(qū)的NDVI數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:∥www．resdc．cn/)“中國(guó)長(zhǎng)時(shí)間序列1 km植被指數(shù)(NDVI)數(shù)據(jù)集”，該數(shù)據(jù)集是基于SPOT-Vegetation 1 km每10天合成的最大化NDVI，利用最大值合成法 MVC(Maximun Value Composites)，進(jìn)一步合成的中國(guó)1 km、每年12個(gè)月份的植被指數(shù)數(shù)據(jù)集，筆者利用1999—2015年共17年12個(gè)月份的植被指數(shù)數(shù)據(jù)，應(yīng)用ArcGIS軟件平臺(tái)，計(jì)算得出該研究區(qū)1999—2015年每年的年最大植被指數(shù)，該指數(shù)可反映當(dāng)年地表植被生長(zhǎng)最好的狀況［26］。氣象數(shù)據(jù)是基于全國(guó)2 400多個(gè)氣象站點(diǎn)日觀測(cè)數(shù)據(jù)，通過(guò)整理、計(jì)算和空間插值處理生成，并裁剪所得，空間分辨率為1 km×1 km。土地覆蓋數(shù)據(jù)來(lái)源于中國(guó)1∶10萬(wàn)比例尺土地利用現(xiàn)狀遙感監(jiān)測(cè)數(shù)據(jù)庫(kù)，該數(shù)據(jù)集的數(shù)據(jù)是以各期Landsat TM/ETM遙感影像為數(shù)據(jù)源，通過(guò)人工目視解譯生成。

RUE是指地面凈初級(jí)生產(chǎn)力(ANPP)與降水量(P)之間的比值。大量研究表明，NDVI與ANPP之間存在很好的線(xiàn)性關(guān)系［23］，筆者研究利用年累計(jì)NDVI(∑ NDVI)計(jì)算 RUE，即 RUE= ∑ NDVI/P，從而獲得1999—2015年的RUE柵格圖。

2.2 數(shù)據(jù)分析方法

2.2.1 趨勢(shì)分析 一元線(xiàn)性回歸分析可以在每個(gè)像元的基礎(chǔ)上，模擬1999—2015年間，研究區(qū)各年份NDVI、RUE的變化趨勢(shì)，進(jìn)而可以反映研究區(qū)的植被隨時(shí)間的變化規(guī)律，其計(jì)算公式［6，21］如下:式中:S為像元NDVI、RUE的回歸斜率值，可表示NDVI、RUE在研究期間的總體變化趨勢(shì);n為NDVI、RUE的年跨度，此處n=17;Ni代表第i年最大NDVI或RUE值。S＞0，表示 NDVI、RUE隨時(shí)間變化呈上升趨勢(shì)，且S越大，表示上升趨勢(shì)越明顯［6］;反之，表示 NDVI、RUE 隨時(shí)間變化呈下降趨勢(shì)。

趨勢(shì)相關(guān)系數(shù)可以定量給出某種要素隨時(shí)間序列的升降程度，當(dāng)趨勢(shì)相關(guān)系數(shù)＞0時(shí)，表示該要素在時(shí)間序列上，有線(xiàn)性增加的趨勢(shì);反之，表示有線(xiàn)性下降趨勢(shì)。趨勢(shì)相關(guān)系數(shù)的計(jì)算公式如下:

式中:r為NDVI、RUE隨時(shí)間的趨勢(shì)相關(guān)系數(shù);Yi為第i年的年份;Y和N分別為Y和N的多年平均值。

趨勢(shì)顯著性檢驗(yàn)，用于檢驗(yàn)?zāi)撤N趨勢(shì)是否有意義，其顯著性?xún)H代表趨勢(shì)性變化置信程度的高低，與變化快慢無(wú)關(guān)［21］，其計(jì)算公式［27］為

2.2.2 偏相關(guān)分析 偏相關(guān)分析是指當(dāng)2個(gè)變量同時(shí)與第3個(gè)變量相關(guān)時(shí)，可將第3個(gè)變量的影響剔除，只分析另外2個(gè)變量之間的相關(guān)性［6，21］。計(jì)算偏相關(guān)系數(shù)前，需要計(jì)算相關(guān)系數(shù)，計(jì)算公式如下:

式中:rxy為變量x和y的相關(guān)系數(shù);x和y分別為x和y的均值。

偏相關(guān)系數(shù)計(jì)算公式如下:

式中:rxy，z為變量z固定后，變量x和y的偏相關(guān)系數(shù)，即為剔除 z的影響后，x、y 的相關(guān)性;rxy、rxz、ryz分別表示變量x與變量y、變量x與變量z、變量y與變量z的相關(guān)系數(shù)。

偏相關(guān)系數(shù)的顯著性檢驗(yàn)，采用t檢驗(yàn)法，其統(tǒng)計(jì)量公式為

式中m為自由度個(gè)數(shù)，此處m=2。

3 結(jié)果與分析

3.1 張家口清水河上游流域NDVI的時(shí)空變化

3.1.1 植被年際變化 圖2顯示，1999—2015年清水河流域的NDVI值在0.3～0.9之間，且NDVI在研究年限內(nèi)不斷波動(dòng)，并呈上升趨勢(shì)，其變化趨勢(shì)與海河流域NDVI的變化趨勢(shì)相一致［6］，說(shuō)明研究區(qū)在近17年植被有所改善。不同土地利用類(lèi)型之間的NDVI隨年份的波動(dòng)趨勢(shì)基本相同，其中灘地、疏林地、城鎮(zhèn)及農(nóng)村的波動(dòng)幅度較大。主要原因在于這3類(lèi)土地利用類(lèi)型主要為人類(lèi)的活動(dòng)區(qū)域，且具有植被類(lèi)型較單一［28］、分布較零散［29］的特點(diǎn)。不同土地利用之間NDVI值從大到小依次為疏林地、有林地、灌木林地、草地、旱地、城鎮(zhèn)及農(nóng)村、灘地。

圖2 1999—2015年研究區(qū)植被指數(shù)NDVI的年際變化Fig．2 Change of NDVI of the study area in 1999—2015

3.1.2 植被空間分布 研究區(qū)1999—2015年平均植被指數(shù)NDVI的空間分布見(jiàn)圖3?？梢缘贸銮逅由嫌瘟饔蛑脖桓采w有很強(qiáng)的分布特征，植被覆蓋自東向西逐漸降低，東溝流域的植被明顯好于西溝植被，且溝道處的植被覆蓋較低。造成該分布的主要原因?yàn)?研究區(qū)降水量空間變化大(圖4)，降水量為366～582 mm，降水量自東向西逐漸減少，溝道的降水量明顯少于周?chē)貐^(qū);同時(shí)，研究區(qū)西溝的主要土地利用類(lèi)型為灌木林地和旱地，旱地及城鎮(zhèn)、農(nóng)村主要分布在溝道處(圖5)，故研究區(qū)NDVI空間分布差異明顯。

3.1.3 植被時(shí)空變化 通過(guò)對(duì)研究區(qū)域的像元進(jìn)行趨勢(shì)分析，得到清水河上游流域近17年NDVI的空間變化趨勢(shì)(圖6)。根據(jù)趨勢(shì)顯著性檢驗(yàn)及等間距法，對(duì)研究區(qū)各像元的回歸斜率進(jìn)行分類(lèi)，可將研究區(qū)分為:明顯退化、無(wú)顯著變化、輕度改善、中度改善和明顯改善5類(lèi)，具體分類(lèi)標(biāo)準(zhǔn)如表1。結(jié)果表明，清水河上游流域85%為植被改善區(qū)。其中:8.5%為明顯改善區(qū)，主要分布于崇禮縣高家營(yíng)鎮(zhèn)西溝部分、石嘴子鄉(xiāng)以及紅旗營(yíng)鄉(xiāng)的草地、旱地及灌木林地;71.6%為中度改善區(qū)，主要分布于西溝流域、正溝流域以及東溝上游區(qū)域;研究區(qū)的明顯退化區(qū)位于崇禮縣城周邊。

3.2 植被變化的氣候驅(qū)動(dòng)力分析

圖3 研究區(qū)年均NDVI的空間分布Fig．3 Spatial distribution of annual NDVI in the study area

圖4 研究區(qū)年均降水量的空間分布Fig．4 Spatial distribution of annual precipitation in the study area

圖5 研究區(qū)土地利用現(xiàn)狀圖Fig．5 Land use map of the study area

在中國(guó)北方地區(qū)，植被覆蓋度在很大程度上受到氣候條件的影響。其中，降水和氣溫是影響植被生長(zhǎng)的2個(gè)重要?dú)夂蛞蜃?，通過(guò)對(duì)比1999—2015年研究區(qū)整體NDVI與年降水量、年平均氣溫的變化趨勢(shì)，可以得出研究區(qū)整體的NDVI與年降水量的波動(dòng)趨勢(shì)完全相同，年降水量降低，NDVI降低(圖7和圖8)。

通過(guò)對(duì)1999—2015年研究區(qū)各像元的NDVI與年降水量、年平均氣溫進(jìn)行偏相關(guān)性分析發(fā)現(xiàn)，1999—2015年研究區(qū)NDVI與年降水量的偏相關(guān)系數(shù)為0.335±0.073，其中，顯著相關(guān)(P＞0.05)只占1.36%，NDVI與年平均氣溫的偏相關(guān)系數(shù)為0.055±0.101，均不顯著相關(guān)(P＜0.05)。研究表明，研究區(qū)1999—2015年植被變化的主要驅(qū)動(dòng)因子，并不是氣候因子。

圖6 1999—2015年研究區(qū)NDVI變化格局Fig．6 Spatial distribution of NDVI change in thestudy area from 1999 to 2015

表1 1999—2015年研究區(qū)植被變化趨勢(shì)分類(lèi)標(biāo)準(zhǔn)Tab．1 Classification standard for NDVI changes in the study area from 1999 to 2015

圖7 1999—2015年研究區(qū)NDVI與年降水量的變化Fig．7 Changes of annual precipitation and NDVI of the study area in 1999—2015

3.3 植被變化的人為驅(qū)動(dòng)力分析

筆者結(jié)合RUE與NDVI的變化趨勢(shì)，研究人類(lèi)活動(dòng)對(duì)研究區(qū)植被覆蓋的干擾程度。根據(jù)李輝霞等［23］分類(lèi)標(biāo)準(zhǔn)，可將清水河上游流域人類(lèi)干擾強(qiáng)度分為6類(lèi):強(qiáng)度負(fù)向干擾、輕度負(fù)向干擾、無(wú)明顯人為干擾、輕度正向干擾、中度正向干擾及強(qiáng)度正向干擾。圖9顯示，研究區(qū)69%屬于中度正向干擾區(qū)，16%屬于強(qiáng)度正向干擾區(qū)域，14%屬于無(wú)明顯人為干擾區(qū)。

圖8 1999—2015年研究區(qū)NDVI與年平均氣溫之間的變化Fig．8 Changes of annual average temperature and vegetation value of study area in 1999—2015

圖9 1999—2015年研究區(qū)人類(lèi)活動(dòng)強(qiáng)度分區(qū)圖Fig．9 Results of human activity intensity assessment of the study area in 1999—2015

通過(guò)對(duì)比研究區(qū)NDVI變化格局(圖6)與人類(lèi)活動(dòng)強(qiáng)度分區(qū)(圖9)，得出研究區(qū)NDVI明顯退化區(qū)與人類(lèi)強(qiáng)度負(fù)向干擾區(qū)一致，NDVI無(wú)顯著變化區(qū)由人類(lèi)活動(dòng)輕度負(fù)向干擾區(qū)及無(wú)明顯人為干擾區(qū)組成，NDVI明顯改善區(qū)88%為人類(lèi)活動(dòng)強(qiáng)度正向干擾區(qū)，人類(lèi)活動(dòng)中度正向干擾區(qū)的NDVI大部分表現(xiàn)為輕度改善及中度改善。

研究表明，1999—2015年研究區(qū)NDVI變化主要受人類(lèi)活動(dòng)的干擾，且人類(lèi)活動(dòng)對(duì)研究區(qū)植被改善起到很大的促進(jìn)作用。自1999年，退耕還林工程配套荒山荒地造林工程、張家口塞北林場(chǎng)建設(shè)工程、城鄉(xiāng)綠化工程、新民居“百村綠化”工程和“雙百里”綠化工程在研究區(qū)相繼開(kāi)展，2010年起，京津冀風(fēng)沙治理工程、清水河上游綜合治理工程、頭道營(yíng)至多樂(lè)美地公路沿線(xiàn)綠化工程、退耕還林補(bǔ)植補(bǔ)造工程及之后的奧運(yùn)綠化工程等的實(shí)施，對(duì)清水河上游流域進(jìn)行全覆蓋式的植被改善。西溝高家營(yíng)鎮(zhèn)、石嘴子鄉(xiāng)和驛馬圖鄉(xiāng)的植被條件較差，水土流失嚴(yán)重［30，31］，故 該區(qū)域?yàn)楫?dāng) 地水 土保 持的 重 點(diǎn) 治 理區(qū)［30］，經(jīng)封山育林、人工造林后，植被明顯改善;同時(shí)，崇禮縣城的開(kāi)發(fā)建設(shè)，使得縣城周?chē)暮档刂脖幻黠@退化。在研究區(qū)人為中度正向干擾區(qū)，東溝下游部分主要為有林地、灌木林地和草地，且植被生長(zhǎng)較好［29］:因此，在人類(lèi)的中度干擾下，該區(qū)植被表現(xiàn)為輕度改善，而其他區(qū)域在人類(lèi)的中度干擾下，表現(xiàn)為中度改善。

4 結(jié)論與討論

清水河上游流域植被從東向西逐漸減少，且溝道處植被明顯比周?chē)脖徊?。主要原因在于研究區(qū)的降水量從東南向西北逐漸減少，且西溝海拔明顯低于東溝，而旱地主要分布在地勢(shì)較低，且平緩的溝道處。

1999—2015年，清水河上游流域的植被不斷波動(dòng)，并呈現(xiàn)上升趨勢(shì)，通過(guò)分析研究區(qū)NDVI與年降水量、年平均氣溫的變化，得出該區(qū)植被的波動(dòng)主要受年降水量波動(dòng)的影響;通過(guò)分析研究區(qū)各像元1999—2015年NDVI與年降水量、年平均氣溫的偏相關(guān)性，得出該區(qū)NDVI與年降水量、年平均氣溫均無(wú)顯著偏相關(guān)性，主要原因在于該區(qū)NDVI在研究年限內(nèi)的變化，主要由人類(lèi)活動(dòng)干擾造成。

1999—2015年，清水河上游流域85%的區(qū)域?yàn)槿藶檎蚋蓴_區(qū)，說(shuō)明在國(guó)家實(shí)施退耕還林、京津冀風(fēng)沙源治理之后，研究區(qū)進(jìn)行大面積的植樹(shù)造林及封山保護(hù)，從而使得研究區(qū)85%的區(qū)域植被顯著改善。在人類(lèi)正向干擾作用下，植被改善后的研究區(qū)依然呈現(xiàn)東溝植被好于西溝，說(shuō)明氣候因子依然是限制研究區(qū)植被空間分布的條件。

［1］ SUN Hongyu，WANG Changyao，NIU Zheng，et al．A-nalysis of the vegetation cover change and the relationship between NDVI and environmental factors by using NOAA time series data［J］．Journal of Remote Sensing，1998，266(2):153．

［2］ 崔林麗，史軍，肖風(fēng)勁，等．中國(guó)東部NDVI的變化趨勢(shì)及其與氣候因子的相關(guān)分析［J］．資源科學(xué)，2010，32(1):124．CUI Linli，SHI Jun，XIAO Fengjing，et al．Variation trends in vegetation NDVI and its correlation with climatic factors in eastern China［J］．Resources Science，2010，32(1):124．

［3］ 李春暉，楊志峰．黃河流域NDVI時(shí)空變化及其與降水/徑流關(guān)系［J］．地理研究，2004，23(6):753．LI Chunhui，YANG Zhifeng．Spatio-temporal changes of NDVI and their relations with precipitation and runoff in the Yellow River Basin［J］．Geographical Research，2004，23(6):753．

［4］ 賈艷紅，趙傳燕，南忠仁．西北干旱區(qū)黑河下游植被覆蓋變化研究綜述［J］．地理科學(xué)進(jìn)展，2007，26(4):64．JIA Yanhong，ZHAO Chuanyan，NAN Zhongren．Review of study on vegetation cover change in the lower reaches of Heihe River in Northwest China arid area［J］．Progress in Geography，2007，26(4):64．

［5］ ZHAO Maosheng，F(xiàn)U Congbin，YAN Xiaodong，et al．Study on the relationship between different ecosystems and climate in China using NOAA/AVHRR data［J］．Acta Geographica Sinica，2001，56(3):287．

［6］ 王永財(cái)，孫艷玲，王中良．1998—2011年海河流域植被覆蓋變化及氣候因子驅(qū)動(dòng)分析［J］．資源科學(xué)，2014，36(3):594．WANGYongcai，SUN Yanling，WANGZhognliang．Spatial-temporal change in vegetation cover and climate factor drivers of variation in the Haihe River Basin 1998－2011［J］．Resources Science，2014，36(3):594．

［7］ 楊紹鍔，廖雪萍，譚裕模，等．廣西近十年植被NDVI變化及其對(duì)降水的響應(yīng)特征分析［J］．西南農(nóng)業(yè)學(xué)報(bào)，2013，26(2):766．YANG Shaoe，LIAOXueping，TAN Yumo，et al．Analysis of Guangxi vegetation NDVI variations and response to precipitation in recent decade［J］．Southwest China Journal of Agricultural Sciences，2013，26(2):766．

［8］ 朱源，王志，彭光雄．基于SPOT-NDVI的賀蘭山植被覆蓋的時(shí)空變化研究［J］．遙感技術(shù)與應(yīng)用，2012，27(1):142．ZHU Yuan，WANGZhi，PENG Guangxiong．The temporal and spatial changes of vegetation cover on Helan ranges based on SPOT-NDVI［J］．Remote Sensing Technology and Application，2012，27(1):142．

［9］ YANG Shengtian，ZHOU Xu，LIU Xiaoyan，et al．Impacts of land use on vegetation coverage in the rich and coarse sediment area of Yellow River Basin［J］．Acta Geographica Sinica，2014，69(1):31．

［10］劉軍會(huì)，高吉喜．氣候和土地利用變化對(duì)中國(guó)北方農(nóng)牧交錯(cuò)帶植被覆蓋變化的影響［J］．應(yīng)用生態(tài)學(xué)報(bào)，2008，19(9):2016．LIU Junhui，GAO Jixi．Effects of climate and land use change on the changes of vegetation coverage in farmingpastoral ecotone of Northern China［J］．Chinese Journal of Applied Ecology，2008，19(9):2016．

［11］ POCCARD Y R I．RICHARD Y，POCCARD I．A statistical study of NDVI sensitivity to seasonal and interannual rainfall variations in Southern Africa［J］．International Journal of Remote Sensing，1998，19(15):2907．

［12］ SCANLON T M，ALBERTSON J D，CAYLOR K K，et al．Determining land surface fractional cover from NDVI and rainfall time series for a savanna ecosystem［J］．Remote Sensing of Environment，2002，82(2/3):376．

［13］陳曉光，李劍萍，李志軍，等．青海湖地區(qū)植被覆蓋及其與氣溫降水變化的關(guān)系［J］．中國(guó)沙漠，2007，27(5):797．CHEN Xiaoguang，LI Jianping LI Zhijun，et al．Vegetation Coverage and Its Relation with Climate Change in Qinghai Lake Area［J］．Journal of Desert Research，2007，27(5):797．

［14］武永利，李智才，王云峰，等．山西典型生態(tài)區(qū)植被指數(shù)(NDVI)對(duì)氣候變化的響應(yīng)［J］．生態(tài)學(xué)雜志，2009，28(5):925．WU Yongli，LI Zhicai，WANG Yunfeng，et al．Responses of vegetation index(NDVI)in typical ecological areas of Shanxi Province to climate change［J］．Chinese Journal of Ecology，2009，28(5):925．

［15］李錦榮，孫保平，凌俠，等．退耕還林工程前后土地利用與植被蓋度動(dòng)態(tài)變化研究:以安塞縣為例［J］．內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2010，31(3):130．LI Jinrong，SUN Baoping，LING Xia，et al．Study on land use and dynamic changes of vegetation cover in and after the project of returning farmland to forests:taking the Ansai County as example［J］．Journal of Inner Mongolia Agricultural University(Natural Science Edition)，2010，31(3):130．

［16］毛德華，王宗明，宋開(kāi)山，等．東北多年凍土區(qū)植被NDVI變化及其對(duì)氣候變化和土地覆被變化的響應(yīng)［J］．中國(guó)環(huán)境科學(xué)，2011，31(2):283．MAO Dehua，WANG Zongming，SONG Kaishan，et al．The vegetation NDVI variation and its responses to climate change and LUCC from 1982 to 2006 year in northeast permafrost region［J］．China Environmental Science，2011，31(2):283．

［17］ IMMERZEEL W W，QUIROZ R A．Understanding precipitation patterns and land use interaction in Tibet using harmonic analysis of SPOT VGT-S10 NDVI time series［J］．International Journal of Remote Sensing，2005，26(11):2281．

［18］ WANG Fei，WANG Xuan，ZHAO Ying，et al．Correlation analysis of NDVI dynamics and hydro-meteorological variables in growth period for four land use types of a water scarce area［J］．Earth Science Informatics，2014，7(3):187．

［19］劉斌，孫艷玲，王中良，等．華北地區(qū)植被覆蓋變化及其影響因子的相對(duì)作用分析［J］．自然資源學(xué)報(bào)，2015(1):12．LIU Bin，SUN Yanling，WANG Zhongliang，et al．Analysis of the vegetation cover change and the relative role of Its influencing factors in North China［J］．Journal of Nat-ural Resources，2015，30(1):12．

［20］李輝霞，劉國(guó)華，傅伯杰．基于NDVI的三江源地區(qū)植被生長(zhǎng)對(duì)氣候變化和人類(lèi)活動(dòng)的響應(yīng)研究［J］．生態(tài)學(xué)報(bào)，2011，31(19):5495．LI Huixia，LIU Guohua，F(xiàn)U Bojie．Response of vegetation to climate change and human activity based on NDVI in the Three-River Headwaters region［C］．Acta Ecologica Sinica，2011，31(19):5495．

［21］穆少杰，李建龍，陳奕兆，等．2001—2010年內(nèi)蒙古植被覆蓋度時(shí)空變化特征［J］．地理學(xué)報(bào)，2012，67(9):1255．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．

［22］彭飛，王濤，薛嫻．基于RUE的人類(lèi)活動(dòng)對(duì)沙漠化地區(qū)植被影響研究:以科爾沁地區(qū)為例［J］．中國(guó)沙漠，2010，30(4):896．PENG Fei，WANG Tao，XUE Xian．Study of human impact on vegetation in desertified regions based on rainfall use efficiency:a case of Horqin region，Inner Mongolia［J］．Journal of Desert Research，2010，30(4):896．

［23］李輝霞，周紅藝，魏興琥．基于RUE和NDVI的人類(lèi)活動(dòng)對(duì)植被干擾強(qiáng)度分析:以桂西北為例［J］．中國(guó)沙漠，2014，34(3):927．LI Huixia，ZHOU Hongyi，WEI Xinghu．Analysis of the impact of human disturbance on vegetation based on RUE and NDVI:a case study in Northwest Guangxi，China［J］．Journal of Desert Research，2014，34(3):927．

［24］侯美英．抓好清水河上游綜合治理建設(shè)生態(tài)宜居城市——張家口市清水河上游水土保持綜合治理措施淺析［J］．中國(guó)化工貿(mào)易，2013(12):333．HOU Meiying．Grasp the comprehensive management of the upper reaches of Qingshui River and build an ecological livable city．Analysis of comprehensive measures for soil and water conservation in the upper reaches of Zhangjiakou，Qingshui River．［J］．chemical industry of China，2013(12):333．

［25］劉三龍．清水河張家口以上流域徑流變化分析［C］．中國(guó)水利學(xué)會(huì)2010學(xué)術(shù)年會(huì)，貴陽(yáng)，2010:215．LIU Sanlong．Runoff variation analysis over Zhangjiakou basin in Qingshui River［C］．Proceedings of the 2010 Annual Conference of China Society of water conservancy，Guiyang，2010:215．

［26］李曉兵，王瑛，李克讓?zhuān)甆DVI對(duì)降水季節(jié)性和年際變化的敏感性［J］．地理學(xué)報(bào)，2000，55(b11):82．LI XiaoBing，WANGYing，LI Kerang．NDVI Sensitivity to Seasonal and Interannual Rainfall Variations in Northern China［J］．Acta Geographica Sinica，2000，55(b11):82．

［27］劉振元，張杰，陳立．青藏高原植被指數(shù)最新變化特征及其與氣候因子的關(guān)系［J］．氣候與環(huán)境研究，2017，22(3):289．LIU Zhenyuan， ZHANG Jie， CHEN Li． The latest change in the Qinghai-Tibetan Plateau vegetation index and its relationship with climate factors［J］．Climatic ＆Environmental Research，2017，22(3):289．

［28］賴(lài)程艷．生態(tài)城鎮(zhèn)建設(shè)中植被的配置［J］．經(jīng)濟(jì)技術(shù)協(xié)作信息，2011(9):111．LAI Chengyan．Configuration of vegetation in the ecological construction［J］．Economic and Technological Cooperation Information，2011(9):111．

［29］蔣軒，周堅(jiān)華．針對(duì)分層分類(lèi)和多描述符空間的城鎮(zhèn)植被群分類(lèi)［J］．遙感信息，2015(2):43．JIANG Xun，ZHOU Jianhua，Classification of urban vegetation population through hierarchical classification in multi-descriptor space［J］．Remote Sensing Information，2015(2):43．

［30］班富孝．清水河流域水文特性［J］．南水北調(diào)與水利科技，2000(1):11．BAN Fuxiao．Hydrological characteristics of Qingshui River Basin［J］．South to North Water Transfer and water conservancy science and technology，2000(1):11．

［31］班富孝．清水河流域上游水土流失及其治理措施探討［J］．地下水，2010，32(3):171．BAN Fuxiao．Study on soil erosion and its control measures in the upper reaches of Qingshui River Basin［J］．Groundwater，2010，32(3):171．