黃杉杉 馬超

摘 要：基于2000—2018年MODIS-NDVI數(shù)據(jù)（NDVI為歸一化植被指數(shù)）、SRTM DEM數(shù)據(jù)和氣象數(shù)據(jù)，利用趨勢(shì)分析、偏差分析和相關(guān)性分析方法，在時(shí)間和空間尺度上，定量分析了紫金山地區(qū)NDVI的時(shí)空變化規(guī)律以及與氣候因子的響應(yīng)關(guān)系，結(jié)果表明：2000—2018年，紫金山地區(qū)NDVI增長(zhǎng)速率為0.015/a（顯著性水平p<0.01）;紫金山地區(qū)NDVI得到顯著改善，其時(shí)間變化規(guī)律大體是，前9 a完成了由低覆蓋度向中等覆蓋度的轉(zhuǎn)化，后10 a完成了由中等覆蓋度向高覆蓋度轉(zhuǎn)化，19 a區(qū)域NDVI顯著改善面積達(dá)64.85%;研究區(qū)氣候呈暖濕化傾向，19 a氣溫增長(zhǎng)率為0.04 ℃/a，降水量增長(zhǎng)率為14.63 mm/a;年NDVI均值與年降水量呈非常顯著正相關(guān)關(guān)系（相關(guān)系數(shù)r=0.734，顯著性水平p<0.01），76.61%的區(qū)域NDVI與年降水量呈極高顯著正相關(guān)關(guān)系，與年均氣溫相關(guān)性不顯著;自2000年以來，紫金山地區(qū)NDVI整體向好的趨勢(shì)發(fā)展，中低海拔山區(qū)增長(zhǎng)顯著，降水是影響該地區(qū)NDVI變化的主要因素。

關(guān)鍵詞：紫金山地區(qū);歸一化植被指數(shù);時(shí)空變化;氣候變化;相關(guān)性

中圖分類號(hào)：X321;P237 ? 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.12.007

引用格式：黃杉杉，馬超.呂梁紫金山地區(qū)NDVI時(shí)空變化與氣候響應(yīng)[J].人民黃河，2021，43（12）：35-40.

Abstract： Based on MODIS-NDVI data， SRTM DEM data and China meteorological data from 2000 to 2018 and using the methods of trend analysis， deviation analysis and correlation analysis， it quantitatively analyzed the spatial-temporal changes and the response relationship with climate factors of the NDVI in the Zijin Mountain area on spatial and temporal scales. the results show that a） from 2000 to 2018， the annual average growth rate is 0.015/a （significance test p<0.01）; b） the degree of the land desertification in the Zijin Mountain area has been improved significantly. The rules of land cover change are generally that the conversion of desertification land to arid grassland is completed in the first 9 years and the conversion of arid grassland to shrub land is completed in the later 10 years. The area of NDVI has been significantly improved by 64.85% in 19 years; c） the climate shows a tendency to warm and humidify， with the temperature increase rate of 0.04 ℃/a and the precipitation growth rate of 14.63 mm/a in 19 years; d） the annual average NDVI is significantly positively correlation with the annual precipitation （correlation coefficient r=0.734， passed the significance test p<0.01）， and 76.61% of the area has a significant positive correlation with the precipitation， which is not significantly correlated with the annual average temperature. Since 2000， the NDVI of vegetation in the Zijin Mountain area has generally been developed towards a good trend and the growth trend of NDVI in low-mid altitude mountain areas is particularly significant. Climate is the main factor affecting the change of NDVI in this area.

Key words： Zijin Mountain area; normalized difference vegetation index; spatial-temporal change; climate change; correlation

2015年聯(lián)合國(guó)發(fā)展峰會(huì)通過了《改變我們的世界：2030年可持續(xù)發(fā)展議程》（簡(jiǎn)稱《2030議程》），將應(yīng)對(duì)氣候變化與保護(hù)陸地生態(tài)系統(tǒng)列入全球17項(xiàng)可持續(xù)發(fā)展目標(biāo)[1]。通過遙感技術(shù)獲得的歸一化植被指數(shù)（NDVI）作為一種監(jiān)測(cè)區(qū)域植被生長(zhǎng)狀況的有效量化指標(biāo)[2]，是表征陸地生態(tài)系統(tǒng)結(jié)構(gòu)和功能穩(wěn)定性的重要指示器[3]。黃土高原位于半濕潤(rùn)與干旱半干旱過渡帶，土地荒漠化和生態(tài)環(huán)境脆弱是該地區(qū)的典型特征[4]。黃土丘陵溝壑區(qū)是黃土高原乃至全球水土流失最嚴(yán)重的地區(qū)之一，地形破碎，溝壑縱橫[5]。降雨、地形、植被等自然因素是該地區(qū)水土流失的主要因素[6]。研究表明，黃土高原地區(qū)在2000年以后植被覆蓋狀況趨于好轉(zhuǎn)[7-8]，超過50%的區(qū)域得到顯著改善[9-11]。諸多學(xué)者對(duì)黃土高原不同區(qū)域的植被覆蓋狀況進(jìn)行了研究，李雙雙等[12]基于2000—2009年MODIS-NDVI數(shù)據(jù)，采用線性趨勢(shì)分析、Hurst指數(shù)和偏相關(guān)分析法，對(duì)陜甘寧地區(qū)植被覆蓋時(shí)空變化特征進(jìn)行了分析;張雅麗等[13]基于2000—2010年逐月MODIS-NDVI數(shù)據(jù)，采用最大值合成法和多元線性回歸方法，分析了黃土丘陵溝壑區(qū)志丹縣的植被覆蓋變化及其驅(qū)動(dòng)力;趙安周等[14]基于2000—2014年MODIS-NDVI數(shù)據(jù)，利用趨勢(shì)分析、Mann-Kendall檢驗(yàn)、Hurst指數(shù)等方法，對(duì)黃土高原植被覆蓋時(shí)空變化及其原因進(jìn)行了分析。在NDVI對(duì)氣候的響應(yīng)方面，李登科等[15]基于1982—2007年GIMMS和SPOT VGT-NDVI數(shù)據(jù)，采用最大值合成法、趨勢(shì)分析和相關(guān)性分析方法對(duì)陜北黃土丘陵溝壑區(qū)植被覆蓋變化及其對(duì)氣候的響應(yīng)進(jìn)行了分析;張含玉等[16]基于1999—2013年SPOT VGT-NDVI數(shù)據(jù)，采用最大值合成法、一元線性回歸和偏相關(guān)分析法，分析了黃土高原地區(qū)NDVI時(shí)空變化趨勢(shì)及其對(duì)氣候的響應(yīng);孟晗等[17]基于2001—2015年SPOT VGT-NDVI數(shù)據(jù)集，利用最大值合成法、趨勢(shì)分析和相關(guān)性分析方法，對(duì)黃土高原NDVI與氣候因子的空間分布特征及其響應(yīng)關(guān)系進(jìn)行了研究。盡管該領(lǐng)域研究已持續(xù)多年，但已有的研究仍存在不足，主要表現(xiàn)在：多數(shù)研究以黃土高原整體為研究區(qū)，區(qū)域尺度大，忽視了典型小流域的植被覆蓋變化特征，難以指導(dǎo)當(dāng)?shù)厥?、縣以及山地的環(huán)境保護(hù)工作;黃土丘陵溝壑區(qū)具有典型的地域差異性，不同地區(qū)的NDVI、不同植被類型的空間差異性顯著，多數(shù)研究集中于人為因素和氣候條件的共同作用，對(duì)于更深入地剖析自然地理要素對(duì)NDVI的影響的研究相對(duì)薄弱;多數(shù)研究選用NDVI或植被覆蓋度（FVC）作為分析區(qū)域植被狀況的指標(biāo)，指標(biāo)略顯單一。另外，選用時(shí)間序列短、間隔大的遙感數(shù)據(jù)難以揭示植被生態(tài)結(jié)構(gòu)的長(zhǎng)期演化規(guī)律。

呂梁山脈南段紫金山地區(qū)屬于黃土丘陵溝壑區(qū)，是一處獨(dú)立的地理單元，形成天然的內(nèi)陸島嶼型封閉生境，是研究全球生態(tài)環(huán)境變化和內(nèi)陸島嶼生態(tài)演化的天然實(shí)驗(yàn)室。研究采用2000—2018年共19 a的MODIS-NDVI數(shù)據(jù)、氣象數(shù)據(jù)，分析區(qū)域內(nèi)植被類型分布規(guī)律、植被指數(shù)與氣候的相關(guān)性，旨在為相關(guān)研究提供參考。

1 研究區(qū)概括

紫金山地區(qū)南北長(zhǎng)81.5 km，東西寬47.5 km，面積為2 441 km2。研究區(qū)總體呈孤立的穹隆狀隆起，四周為黃土覆蓋，地勢(shì)總體北高南低，主峰海拔1 823 m，最低點(diǎn)為西南角（磧口鎮(zhèn)）的黃河河床，海拔700 m，區(qū)內(nèi)相對(duì)高差為1 123 m。山坡平緩、水系呈放射狀，黃河流經(jīng)臨縣呈弧形轉(zhuǎn)彎。

研究區(qū)屬于溫帶大陸性氣候區(qū)，年平均氣溫為7～9 ℃，年降水量為300～600 mm，年平均風(fēng)速2.0～2.5 m/s，無霜期135～180 d。地貌類型屬于典型的黃土丘陵溝壑區(qū)[18]，溝壑縱橫，溝壑以紫金山為中心向四周呈放射狀分布，水土流失嚴(yán)重。

研究區(qū)森林、灌叢、干旱草地、荒漠、裸地各種覆蓋類型集中分布，構(gòu)成了當(dāng)?shù)氐湫偷纳降厣鷳B(tài)景觀。土壤主要為棕壤、褐土、黃綿土、栗褐土和草甸土[19]。植物按海拔垂直分布的層次較為明顯，既有淺山區(qū)的灌木林帶，也有中山區(qū)的楊樺林帶、落葉松林帶，還有高海拔分布區(qū)的云杉、青杄、紅樺等。森林以油松和落葉松針葉天然次生林為主。灌木樹種主要有胡枝子、虎榛子、繡線菊、沙棘、山桃、黃刺玫、小葉鼠李等，草本植物有苔草、堿草、蒿類[19]。

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

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

研究選用美國(guó)LAADS DAAC發(fā)布的2000—2018年MODIS歸一化植被指數(shù)，其空間分辨率為1 000 m。取植被生長(zhǎng)最旺盛的8月份數(shù)據(jù)，通過MRT軟件數(shù)據(jù)投影方式轉(zhuǎn)化為通用橫軸墨卡托投影（UTM投影），坐標(biāo)系統(tǒng)定義為WGS 1984 UTM Zone 49°N，完成NDVI數(shù)據(jù)的格式轉(zhuǎn)換。

采用2018年美國(guó)陸地資源衛(wèi)星（Landsat 8）的OLI 30 m遙感影像數(shù)據(jù)，通過支持向量機(jī)（SVM）的監(jiān)督分類方法對(duì)遙感影像進(jìn)行土地覆被分類。

數(shù)字高程數(shù)據(jù)源于美國(guó)太空總署（NASA）、美國(guó)國(guó)家影像制圖局（NIMA）以及德國(guó)與意大利航天機(jī)構(gòu)聯(lián)合測(cè)量的航天飛機(jī)雷達(dá)地形測(cè)繪數(shù)字高程模型（SRTM3 DEM），水平分辨率為90 m，垂直分辨率為1 m，主要用于匯水線提取，進(jìn)而劃定研究區(qū)。

氣象數(shù)據(jù)為來自中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心2000—2015年的全國(guó)氣象資料，是基于全國(guó)2 400多個(gè)氣象站點(diǎn)的日觀測(cè)數(shù)據(jù)，通過整理、計(jì)算和空間插值生成的。本研究對(duì)氣象數(shù)據(jù)進(jìn)行投影坐標(biāo)系的轉(zhuǎn)換，投影坐標(biāo)系統(tǒng)定義為WGS 1984 UTM Zone 49° N，并提取研究區(qū)年平均氣溫、降水量空間插值數(shù)據(jù)集（2016—2018年氣象數(shù)據(jù)資料缺少，因此2016—2018年氣溫、降水?dāng)?shù)據(jù)取自研究區(qū)附近氣象站點(diǎn)離石縣的逐年、逐月氣象記錄），用于分析紫金山地區(qū)2000—2018年植被指數(shù)與氣候因素之間的關(guān)系。

2.2 研究方法

3 結(jié)果與分析

3.1 歸一化植被指數(shù)空間分布規(guī)律

取2018年NDVI分布圖（見圖1（a）），對(duì)照2018年紫金山地區(qū)土地覆被分布圖（見圖1（b）），定性表達(dá)紫金山地區(qū)主要土地覆被的空間分布規(guī)律。土地覆被類型與NDVI閾值區(qū)間的相關(guān)分析表明，NDVI高值（0.8≤NDVI≤1.0）主要分布在紫金山，土地覆被類型為林地;NDVI中值（0.4≤NDVI<0.8）主要分布于中海拔地區(qū)，土地覆被類型為草地;低值（0.0≤NDVI<0.4）主要分布于研究區(qū)北部以及黃河沿岸低海拔地區(qū)，土地覆被類型為其他土地。

19 a平均NDVI分布可以定性宏觀表達(dá)紫金山地區(qū)植被景觀的空間變化規(guī)律，由圖1（c）可以看出，NDVI由紫金山向外逐漸遞減，平均NDVI最低值為黃河沿岸的其他土地，最高值為紫金山的林地。19 a平均NDVI整體呈現(xiàn)西部和西北部低，東部高的分布特征，具有明顯的空間差異性，與黃土高原2000—2016年的植被分布特點(diǎn)趨于一致[24]。

為了進(jìn)一步揭示紫金山地區(qū)植被類型的變化趨勢(shì)，采用一元線性回歸方法對(duì)2000—2018年NDVI進(jìn)行逐像元趨勢(shì)分析。圖1（d）顯示，2000—2018年NDVI整體呈明顯改善趨勢(shì)，19 a間NDVI無顯著變化的區(qū)域僅9 km2，占區(qū)域總面積的0.37%;NDVI輕度改善區(qū)域面積約為131 km2，占區(qū)域總面積的5.37%，主要為紫金山和黃河沿岸、湫水河一帶的居民區(qū);NDVI中度改善區(qū)面積約為702 km2，占區(qū)域總面積的28.76%，主要分布于黃河和湫水河沿岸的農(nóng)業(yè)區(qū);64.85%的區(qū)域NDVI呈顯著改善趨勢(shì)，面積約為1 583 km2，主要分布于研究區(qū)中低海拔地區(qū);剩余0.65%的區(qū)域呈高顯著改善趨勢(shì)?？梢?，在這19 a里區(qū)域NDVI顯著改善。

3.2 NDVI時(shí)間變化分析

將2000—2018年紫金山地區(qū)植被NDVI分為5級(jí)，并對(duì)結(jié)果進(jìn)行統(tǒng)計(jì)分析，見表1。表1顯示，2000—2018年紫金山地區(qū)NDVI均值為0.32～0.64，NDVI增長(zhǎng)率為0.015/a（顯著性水平p<0.01），比同為丘陵溝壑區(qū)的皇甫川流域2000—2015年的NDVI增長(zhǎng)率0.004 2/a高[25]。2001年NDVI總值與均值最小，2018年NDVI總值和均值都高于其他年份的。均值偏差分析顯示，以2009年為轉(zhuǎn)折點(diǎn)，2000—2008年的NDVI均值低于2000—2018年NDVI均值（0.47），2009—2018年（除了2010年和2015年）NDVI均值高于2000—2018年NDVI均值，總體上生態(tài)虧缺的年份多于生態(tài)盈余的，兩階段NDVI總量分別占19 a NDVI總量的39.2%和60.8%。

為進(jìn)一步分析紫金山地區(qū)植被覆蓋結(jié)構(gòu)的變化規(guī)律，將NDVI各區(qū)間面積占比繪制成百分比柱狀圖，見圖2。由圖2可知，NDVI為[0.0，0.2）的面積占比極低;NDVI為[0.2，0.4）的面積占比呈波動(dòng)減少趨勢(shì)，其中2005年達(dá)到最大值2 122 km2，占區(qū)域總面積的86.93%，2018年面積僅為28 km2，減少幅度達(dá)98.68%;NDVI為[0.4，0.6）的面積波動(dòng)幅度較大，呈先增大后減小的趨勢(shì)，2009年面積達(dá)到最大，為2 073 km2，2005年面積最小，為255 km2，19 a間最大變動(dòng)幅度為87.70%;NDVI為[0.6，0.8）的面積呈顯著增大趨勢(shì)，2018年面積達(dá)到1 929 km2，占區(qū)域面積的79.02%，19 a間最大變幅為98.39%;NDVI為[0.8，1.0]的面積呈增大趨勢(shì)，但總面積不大。研究區(qū)NDVI時(shí)間變化規(guī)律大體是，前9 a（2000—2008年）完成了低覆蓋度向中等覆蓋度的轉(zhuǎn)化，后10 a（2009—2018年）完成了中等覆蓋度向高覆蓋度的轉(zhuǎn)化，19 a植被覆蓋結(jié)構(gòu)逐漸向好的趨勢(shì)發(fā)展。

3.3 NDVI空間變化分析

為了回溯紫金山地區(qū)NDVI的垂直地帶性變遷，研究設(shè)置了遙感生態(tài)剖面，即在紫金山的最高處，分別對(duì)2000—2018年的遙感影像沿山體部分的東西、南北兩個(gè)方向做剖面（W-E水平剖面、S-N垂直剖面），得到紫金山地區(qū)NDVI時(shí)間序列剖面渲染圖，見圖3，其中藍(lán)色直線為紫金山最高點(diǎn)（海拔1 823 m）。由圖3可以看出，海拔越高植被覆蓋越好，在紫金山山頂區(qū)域NDVI達(dá)到最大值，主要植被類型為林地（0.8≤NDVI≤1.0），2000—2009年紫金山水平地帶性和垂直地帶性NDVI分布無顯著差異，覆被類型主要為其他土地（0.0≤NDVI<0.4）。2010年以后，NDVI發(fā)生了改變，并隨時(shí)間逐漸向高值轉(zhuǎn)化。在W-E剖面變化圖中，NDVI在山體水平地帶分布上東部?jī)?yōu)于西部。2010—2018年，紫金山西部地區(qū)主要植被類型為干旱草地（0.4≤NDVI<0.8），東部地區(qū)植被類型為林地。在N-S剖面變化圖中，NDVI在山體垂直地帶分布上南部?jī)?yōu)于北部。2010—2018年紫金山北部地區(qū)植被類型主要為草地，南部地區(qū)主要為草地和林地。

3.4 NDVI對(duì)氣候的響應(yīng)分析

（1）氣候變化分析。2000—2018年紫金山地區(qū)呈暖濕化趨勢(shì)，見圖4。年均氣溫呈波動(dòng)上升趨勢(shì)，變化平穩(wěn)，其中年平均最高氣溫出現(xiàn)在2017年，為11.3 ℃，年平均最低氣溫為2012年的9.3 ℃。19 a平均氣溫為10.4 ℃，氣溫變化率為0.04 ℃/a。年降水量呈顯著上升趨勢(shì)，變化幅度較大，達(dá)496 mm，其中年降水量最少的年份出現(xiàn)在2005年，為325.2 mm，2016年降水量最多，為821.2 mm。19 a平均降水量為520.9 mm，降水量變化率為14.63 mm/a。

（2）NDVI與氣候變化關(guān)系分析。紫金山地區(qū)逐年平均NDVI與年平均氣溫的皮爾遜（Pearson）相關(guān)系數(shù)r=0.352（顯著性水平p<0.01），總體為不顯著相關(guān);與年降水量的皮爾遜相關(guān)系數(shù)r=0.734（p<0.01），總體為非常顯著正相關(guān)。逐像元分析NDVI與年均氣溫和年降水量的空間相關(guān)性，見圖5。圖5 NDVI與氣候相關(guān)系數(shù)空間分布

根據(jù)顯著性程度以及相關(guān)性，將相關(guān)系數(shù)分為6個(gè)等級(jí)：不顯著負(fù)相關(guān)r為[-0.10，0.00）、不顯著正相關(guān)r為[0.00，0.10）、低顯著正相關(guān)r為[0.10，0.25）、顯著正相關(guān)r為[0.25，0.40）、高顯著正相關(guān)r為[0.40，0.65）、極高顯著正相關(guān)r為[0.65，0.85]。

由圖5（a）可知，NDVI與年平均氣溫大多為顯著正相關(guān)，面積約為1 553 km2，占研究區(qū)總面積的63.62%;紫金山西部22.49%的區(qū)域NDVI與年平均氣溫為高顯著正相關(guān)，面積約為549 km2;東部12.74%的區(qū)域NDVI與年平均氣溫為低顯著正相關(guān)，面積約為311 km2;相關(guān)系數(shù)在不顯著負(fù)相關(guān)、不顯著正相關(guān)和極高顯著正相關(guān)的面積占比非常小，分別占總面積的0.33%、0.74%、0.08%。由圖5（b）可知，NDVI與年降水量大多為極高顯著正相關(guān)關(guān)系，面積為1 870 km2，占研究區(qū)總面積的76.61%;高顯著正相關(guān)主要集中分布在研究區(qū)東北部，面積約為561 km2，占總面積的22.98%;不顯著負(fù)相關(guān)、不顯著正相關(guān)和顯著正相關(guān)面積占比不及1%，面積分別為1、3、6 km2?？梢姡杲邓繉?duì)NDVI的影響十分顯著。

4 結(jié) 語

基于2000—2018年MODIS-NDVI數(shù)據(jù)和紫金山地區(qū)年平均氣溫、年降水量空間插值數(shù)據(jù)，從NDVI空間分布規(guī)律、時(shí)間變化規(guī)律和對(duì)氣候響應(yīng)規(guī)律3方面進(jìn)行分析，主要結(jié)論如下。

（1）2000—2018年紫金山地區(qū)NDVI由東部向西部和西北部逐漸遞減，最高值分布在紫金山，最低值分布在黃河沿岸低海拔地區(qū)，具有顯著的空間差異性。NDVI顯著改善，其中中低海拔山區(qū)NDVI改善最為顯著。

（2）2000—2018年紫金山地區(qū)NDVI發(fā)生了很大變化，呈顯著增大趨勢(shì)，NDVI增長(zhǎng)速率為0.015/a。19 a紫金山地區(qū)植被覆蓋結(jié)構(gòu)由低覆蓋向高覆蓋轉(zhuǎn)變。

（3）19 a紫金山地區(qū)呈暖濕化，氣溫增長(zhǎng)率為0.04 ℃/a，降水量增長(zhǎng)率為14.63 mm/a。其中，氣候因素中的降水因子是影響該地區(qū)NDVI變化的主要因素，年NDVI均值與年降水量呈非常顯著正相關(guān)關(guān)系，與年均氣溫相關(guān)性不顯著。

參考文獻(xiàn)：

[1] United Nations General Assembly. Transforming our World： the 2030 Agenda for Sustainable Development[EB/OL]. [2019-09-25].https：//sustainabledevelopment.un.org/post2015/transforminqourword.

[2] 羅敏，古麗·加帕爾，郭浩，等.2000—2013年塔里木河流域生長(zhǎng)季NDVI時(shí)空變化特征及其影響因素分析[J].自然資源學(xué)報(bào)，2017，32（1）：50-63.

[3] 高江波，焦珂?zhèn)ィ瑓墙B洪.1982—2013年中國(guó)植被NDVI空間異質(zhì)性的氣候影響分析[J].地理學(xué)報(bào)，2019，74（3）：534-543.

[4] 張文強(qiáng)，孫從建，李新功.晉西南黃土高原區(qū)植被覆蓋度變化及其生態(tài)效應(yīng)評(píng)估[J].自然資源學(xué)報(bào)，2019，34（8）：1748-1758.

[5] 楊寒月，張光輝，張寶軍.黃土丘陵溝壑區(qū)溝坡典型植物群落生長(zhǎng)特征[J].水土保持研究，2019，26（2）：62-67.

[6] 郭宇，王樹森，馬迎梅，等.降雨因子對(duì)黃土丘陵溝壑區(qū)不同植被毯措施徑流小區(qū)產(chǎn)流產(chǎn)沙量的影響[J].內(nèi)蒙古林業(yè)科技，2019，45（2）：30-34.

[7] LI Jingjing，PENG Shouzhang，LI Zhi.Detecting and Attributing Vegetation Changes on China’s Loess Plateau[J].Agricultural and Forest Meteorology，2017，247（13）：260-270.

[8] 李銳.黃土高原水土保持工作70年回顧與啟示[J].水土保持通報(bào)，2019，39（6）：298-301.

[9] 高健翎，高燕，馬紅斌，等.黃土高原近70 a水土流失治理特征研究[J].人民黃河，2019，41（11）：72-77.

[10] 高健健，穆興民，孫文義.2000—2012年黃土高原植被覆蓋的時(shí)空變化[J].人民黃河，2015，37（11）：85-91.

[11] 郭永強(qiáng)，王乃江，褚曉升，等.基于Google Earth Engine分析黃土高原植被覆蓋變化及原因[J].中國(guó)環(huán)境科學(xué)，2019，39（11）：4804-4811.

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

[13] 張雅麗，王建軍.基于NDVI的陜北黃土高原丘陵溝壑區(qū)植被覆蓋變化及其驅(qū)動(dòng)力分析[J].林業(yè)經(jīng)濟(jì)，2016，38（12）：21-25，40.

[14] 趙安周，劉憲鋒，朱秀芳，等.2000—2014年黃土高原植被覆蓋時(shí)空變化特征及其歸因[J].中國(guó)環(huán)境科學(xué)，2016，36（5）：1568-1578.

[15] 李登科.陜北黃土高原丘陵溝壑區(qū)植被覆蓋變化及其對(duì)氣候的響應(yīng)[J].西北植物學(xué)報(bào)，2009，29（5）：1007-1015.

[16] 張含玉，方怒放，史志華.黃土高原植被覆蓋時(shí)空變化及其對(duì)氣候因子的響應(yīng)[J].生態(tài)學(xué)報(bào)，2016，36（13）：3960-3968.

[17] 孟晗，黃遠(yuǎn)程，史曉亮.黃土高原地區(qū)2001—2015年植被覆蓋變化及氣候影響因子[J].西北林學(xué)院學(xué)報(bào)，2019，34（1）：211-217.

[18] 楊艷芬，王兵，王國(guó)梁，等.黃土高原生態(tài)分區(qū)及概況[J].生態(tài)學(xué)報(bào)，2019，39（20）：7389-7397.

[19] 郝曉鵬.黑茶山自然保護(hù)區(qū)不同生境下植物群落生物多樣性研究[D].太原：山西大學(xué)，2009：1-65.

[20] 趙維清，李經(jīng)緯，褚琳，等.近10年湖北省植被指數(shù)時(shí)空變化特征及其驅(qū)動(dòng)力[J].生態(tài)學(xué)報(bào)，2019，39（20）：7722-7736.

[21] ZHANG Wei，WANG Lunche，XIANG Feifei，et al.Vegetation Dynamics and the Relations with Climate Change at Multiple Time Scales in the Yangtze River and Yellow River Basin， China[J].Ecological Indicators，2020，110（3）：105892.

[22] 袁金國(guó)，張莎，胡春梓.河北省植被多時(shí)相SPOT-VGT NDVI變化及與氣候因子的關(guān)系[J].科技通報(bào)，2015，31（5）：47-52.

[23] 王麗霞，余東洋，劉招，等.渭河流域NDVI與氣候因子時(shí)空變化及相關(guān)性研究[J].水土保持研究，2019，26（2）：249-254.

[24] 孫銳，陳少輝，蘇紅波.2000—2016年黃土高原不同土地覆蓋類型植被NDVI時(shí)空變化[J].地理科學(xué)進(jìn)展，2019，38（8）：1248-1258.

[25] 黎銘，張會(huì)蘭，孟鋮鋮，等.皇甫川流域2000—2015年植被NDVI時(shí)空變化特征[J].林業(yè)科學(xué)，2019，55（8）：36-44.

【責(zé)任編輯 呂艷梅】