趙慧 韓穎娟 張承明 張學藝 李夢華 王艷萍 馬力文 趙金龍

摘要：歸一化植被指數(shù)（NDVI）是表征植被生長狀況的重要因子。基于2000年以來的MODIS NDVI數(shù)據(jù)，運用Mann-Kendall檢驗、趨勢分析、Hurst指數(shù)等方法，研究寧夏年最大NDVI時空變化特征及其與氣候、人類活動的關系。結(jié)果表明：（1）寧夏多年平均最大NDVI呈現(xiàn)南北高、中間低的空間分布特征。（2）19年最大NDVI極顯著上升，增長速率為0.008/年，70.3%的地區(qū)顯著增大（P<0.05），植被生長狀況明顯改善，未來27.8%的區(qū)域最大NDVI將顯著改善。（3）最大NDVI與春季、夏季和年降水量相關性較顯著，與平均氣溫相關性較弱。最大NDVI與地方公共財政收入和地方公共財政支出呈極顯著正相關關系（P<0.01）;大部分地區(qū)NDVI與人口及農(nóng)作物總播種面積正相關，固原市與二者呈極顯著負相關關系（P<0.01）。生態(tài)建設成效尚不穩(wěn)固，需持續(xù)實施生態(tài)工程，改善植被生長狀況。

關鍵詞：NDVI;時空特征;植被變化;氣候變化;人類活動

中圖分類號： Q948.11;TP79 ?文獻標志碼： A ?文章編號：1002-1302（2021）14-0220-06

植被是地球生態(tài)系統(tǒng)的重要組成部分，在地表一系列物質(zhì)與能量交換過程中扮演著關鍵角色，對人類生產(chǎn)生活影響重大。植被生長對氣候變化及人類活動的響應一直是陸地生態(tài)系統(tǒng)的研究重點之一，通過分析植被時空變化特征及其與氣候因子、人類活動的關系，能夠了解影響植被生長的關鍵因素，對應對氣候變化具有重要意義，同時能為生態(tài)建設、植被恢復提供科學指導。近年來，植被指數(shù)成為描述植被生長狀態(tài)的重要指標，廣泛應用于土地利用[1]、植被分類[2]、植被凈初級生產(chǎn)力[3]、森林病蟲害[4]、農(nóng)作物產(chǎn)量預估[5]等研究。其中，歸一化植被指數(shù)（Normalized Differential Vegetation Index，NDVI）是一種能夠反映植被光合作用強度的植被指數(shù)，是近紅外波段與紅光波段的反射率之差與兩者之和的比值，可有效表征地面植被覆蓋情況[6]，是揭示區(qū)域植被生長時空格局的重要指標。

寧夏回族自治區(qū)地處我國西北地區(qū)東部，位于我國北方農(nóng)牧交錯帶，水土流失嚴重，生態(tài)環(huán)境敏感脆弱。2000年以來，寧夏依托三北防護林、天然林保護和退耕還林等國家重大林業(yè)工程，組織實施了封山禁牧、防沙治沙等一系列生態(tài)工程，實施生態(tài)補償制度等，取得了較好的生態(tài)效益。目前，針對全國、西北及寧夏地區(qū)研究分析了AVHRR[7]、MODIS[7-8]、GIMMS[9]、SPOT-VGT[10-12]、NDVI數(shù)據(jù)的時空變化特征及其與氣候因子的關系，表明寧夏植被生長受降水因素影響較大。但已有成果研究時段主要集中在2015年以前，且大多基于較低分辨率的NDVI數(shù)據(jù)，而近年來的研究多集中在更大區(qū)域，專門針對寧夏NDVI變化情況及其與氣象因子、人類活動的關系的探索較少。因此，本研究基于2000—2018年MODIS NDVI 250 m分辨率數(shù)據(jù)，分析了寧夏各區(qū)域NDVI的變化特征及其與氣候因子、人類活動的相關關系，探究影響寧夏植被變化的驅(qū)動因素，以期為寧夏乃至西北地區(qū)的植被生態(tài)恢復治理提供科學決策依據(jù)。

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

1.1 研究區(qū)概況

寧夏（35°14′～39°23′N，104°17′～107°39′E）地處黃河中上游，三面環(huán)沙，大部分區(qū)域?qū)贉貛Т箨懶詺夂?，年平均降水量?00 mm左右，年蒸發(fā)量近 2 000 mm。按照氣候、地形、植被等因素自北向南分為灌區(qū)（引黃灌溉區(qū)）、中部干旱帶及南部山區(qū)。灌區(qū)農(nóng)業(yè)資源豐富，以灌溉綠洲為主;中部干旱帶氣候干燥，植被類型以荒漠化草原為主;南部山區(qū)屬于黃土高原，降水資源較為豐富，植被類型以森林、草原為主。

1.2 數(shù)據(jù)來源及處理

遙感數(shù)據(jù)來源于MOD13Q1 NDVI數(shù)據(jù)，分辨率為250 m。應用MRT工具進行重投影，將全年23期NDVI數(shù)據(jù)進行最大值合成，每個像元選取年內(nèi)最大值作為該像元的值，得到2000—2018年每年的最大NDVI（NDVImax），批量裁切得到研究區(qū)數(shù)據(jù)。地面氣象數(shù)據(jù)來自寧夏氣象綜合數(shù)據(jù)庫，包括24個地面氣象站19年逐日降水量、平均氣溫的數(shù)據(jù)統(tǒng)計得到春夏秋季以及全年值（春季3—5月，夏季6—8月，秋季9—11月），并用ArcGIS進行克里金插值得到柵格數(shù)據(jù)。社會經(jīng)濟數(shù)據(jù)來自《寧夏統(tǒng)計年鑒》，包括地方公共財政收入、地方公共財政支出、人口和農(nóng)作物總播種面積。

1.3 研究方法

1.3.1 Mann-Kendall檢驗 Mann-Kendall檢驗能夠判斷時間序列的變化趨勢。通過計算得到按照時間序列（n個數(shù)據(jù)）正序得到的統(tǒng)計量序列，UFk（k=1，2，…，n），再按照時間序列逆序計算，使得統(tǒng)計量UBk=-UFk，繪制UFk、UBk 2條曲線，判識曲線交點，比較交點、統(tǒng)計量序列及給定顯著性水平臨界值線的關系，判斷時間序列的突變點及變化趨勢[13]。

1.3.2 趨勢分析法 用一元線性回歸得到NDVImax與時間的關系，應用最小二乘法估計NDVImax的變化趨勢，即斜率：

b=n×∑ni=1（i×NDVIi）-∑ni=1i∑ni=1NDVIin×∑ni=1i2-（∑ni=1i）2。（1）

式中：b為趨勢線的斜率;n為19;i為1～19;NDVIi為第i年最大歸一化植被指數(shù)。b表征NDVI的變化速率，當b>0時，NDVImax呈增大趨勢，植被改善;當b<0時，NDVImax呈減小趨勢，植被退化。對斜率b進行顯著性檢驗，變化趨勢分為極顯著（P<0.01）、顯著（P<0.05）、不顯著（P≥0.05）。

1.3.3 Hurst指數(shù) Hurst指數(shù)能夠描述時間序列變化的持續(xù)性。計算數(shù)據(jù)的差分序列，進而計算均值、標準差、累積離差、極差等，得到Hurst指數(shù)[14]。Hurst=0.50，表明時間序列前后的變化無關;Hurst>0.50，時間序列前后變化趨勢相同，系統(tǒng)發(fā)展有持續(xù)性;Hurst<0.50，表明時間序列前后變化趨勢相反，系統(tǒng)發(fā)展具有反持續(xù)性。Hurst越接近0，反持續(xù)性越強，越接近1，持續(xù)性越強。

1.3.4 相關性分析 應用Matlab、SPSS分析近19年寧夏NDVImax與降水量、平均氣溫柵格數(shù)據(jù)以及社會經(jīng)濟指標的相關性，對相關性進行顯著性檢驗。P<0.05表示NDVImax與因子顯著相關。

2 結(jié)果與分析

2.1 NDVImax時空特征

2.1.1 NDVImax空間分布特征 2000—2018年NDVImax平均值高值區(qū)主要位于灌區(qū)北部和南部山區(qū)，中部干旱帶NDVImax較小。六盤山林區(qū)NDVImax最大，黃河流域、清水河流域、南華山林區(qū)、賀蘭山林區(qū)次之，騰格里沙漠東南緣及環(huán)香山地區(qū)較小。灌區(qū)、南部山區(qū)NDVImax>0.50的面積均超過區(qū)域面積25%，灌區(qū)黃河流經(jīng)，水資源豐富，南部山區(qū)降水較多，林草長勢較好。中部干旱帶NDVImax>0.50的區(qū)域僅占6.5%。南部山區(qū)植被總體好于灌區(qū)和中部干旱帶（表1）。

2.1.2 NDVImax時間變化特征 2000—2018年寧夏平均NDVImax呈極顯著上升，增長速率為0.008/年，平均值從2000年的0.29增長至2018年的0.50。2000年以來，NDVImax≤0.25的面積波動減小，NDVImax在0.50～1.00的面積占比以13.9%/年的速率極顯著增加（圖1）。

研究區(qū)2000—2018年NDVImax Mann-Kendall檢驗顯示，寧夏、灌區(qū)、中部干旱帶NDVImax 2009年起發(fā)生突變，南部山區(qū)2008年起發(fā)生突變。灌區(qū)、中部干旱帶、南部山區(qū)分別從2010、2011、2011年起顯著增長。銀川市NDVImax突變年份最早，為2006年，中衛(wèi)市最晚，為2010年。從突變年起，5市NDVImax均呈增長趨勢。石嘴山和銀川地區(qū)NDVImax呈顯著增長的起始年份最早，分別為2007、2008年;中衛(wèi)市最晚，為2012年。中衛(wèi)市建于2004年，為5市最晚。NDVImax突變年份的早晚可能與城市發(fā)展水平及隨之而來的綠化建設有關。

2000—2018年寧夏NDVImax斜率在-0.050～0.056之間。95.1%地區(qū)b>0，NDVImax呈增大趨勢，4.9%地區(qū)呈減小趨勢。斜率通過0.05水平顯著性檢驗區(qū)均達各區(qū)域的50%，高值區(qū)主要位于灌區(qū)北部、大羅山、沿黃河沿清水河流域、吳忠中衛(wèi)市南部及固原市，顯著增加占全區(qū)70.3%，極顯著增加占53.8%，斜率低值區(qū)占全區(qū)的1.5%，主要分布在銀川市、利通區(qū)、沙坡頭區(qū)、中寧縣城區(qū)，NDVImax顯著減小與城市發(fā)展有關。無b=0的區(qū)域。寧夏大部分地區(qū)NDVImax增大，植被狀況總體改善明顯（表2、圖2-a）。

像元尺度2000—2018年寧夏37.7%的區(qū)域Hurst>0.50，NDVImax變化持續(xù)性較強區(qū)主要分布在鹽池市中部、中衛(wèi)市南部以及固原市。其余大部分地區(qū)NDVImax變化呈反持續(xù)性，40.5%區(qū)域Hurst值為0.4～0.5。區(qū)域尺度NDVImax持續(xù)性分析顯示，寧夏Hurst指數(shù)為0.47，中部干旱帶、南部山區(qū)、固原市、中衛(wèi)市Hurst>0.50，NDVI有一定的持續(xù)性;灌區(qū)、銀川市、石嘴山市、吳忠市Hurst<0.50，NDVImax有反持續(xù)性。除石嘴山市外，各區(qū)域 |Hurst-0.5|≤0.10，表明NDVImax持續(xù)性和反持續(xù)性均較弱（表2）。

將NDVImax空間變化趨勢與Hurst指數(shù)分布圖疊加，得到持續(xù)顯著改善（Hurst>0.5，b>0），持續(xù)顯著退化（Hurst>0.5，b<0）、反持續(xù)顯著改善（Hurst<0.5，b>0）、反持續(xù)顯著退化（Hurst<0.5，b<0）的分布。未來NDVImax將顯著改善的區(qū)域占27.8%（持續(xù)顯著改善27.2%，反持續(xù)顯著退化0.6%）（圖2-b）。

2.2 植被與氣候變化的相關性分析

寧夏2000以來氣候發(fā)生了明顯變化。19年降水量增加速率為3.87 mm/年（R=0.474，P<0.05），增濕顯著，3區(qū)5市增濕不顯著。NDVImax與春夏秋季以及年降水量以正相關為主（圖3-a～圖3-d，表3）。春季固原市大部分地區(qū)降水量與NDVImax顯著正相關，相關系數(shù)R較大，賀蘭山、銀川中部、吳忠東部、中衛(wèi)西部和南部部分地區(qū)NDVImax與春季降水顯著正相關。NDVImax與夏季降水顯著正相關的相關系數(shù)高值區(qū)主要位于賀蘭山東麓、沙坡頭區(qū)西部和香山東南部部分地區(qū);固原市夏季降水充足，南部NDVImax與夏季降水相關性不顯著，此時降水不是植被生長的主要限制因子。大部分地區(qū)NDVImax與秋季降水量相關系數(shù)未通過顯著性檢驗。63.5%區(qū)域的NDVImax與年總降水量顯著相關，相關系數(shù)為0.50～1.00且通過顯著性檢驗的區(qū)域主要位于賀蘭山東麓、灌區(qū)南部、中部干旱帶南部以及南部山區(qū)北部。

19年寧夏平均氣溫呈升高趨勢，增溫速率為0.028 ℃/年（R=0.401，P<0.1），增溫不顯著。銀川、中衛(wèi)市增溫明顯，增溫速率分別為0.031 ℃/年（R=0.458，P<0.05）、0.049 ℃/年（R=0.588，P<0.01）。寧夏NDVImax與平均氣溫顯著相關的面積較?。▓D3-e～圖3-h，表3）。春季與NDVImax呈顯著正相關的區(qū)域零星分布，相關系數(shù)大部分為0.25～0.75，春季氣溫高，則草地植被返青早，春播作物苗情好，NDVImax隨溫度升高而增大。NDVImax與夏季平均氣溫顯著負相關主要分布在石嘴山市，顯著正相關主要在銀川市中部和中衛(wèi)市北部。NDVImax與秋季平均氣溫顯著負相關占全區(qū)的2.4%。NDVImax與年平均氣溫顯著相關主要位于北部4市，其中負相關區(qū)主要在銀川市北部。寧夏大部分地區(qū)NDVImax與春夏秋季、年平均氣溫無顯著相關關系。

2.3 植被與人類活動的相關性分析

經(jīng)濟社會發(fā)展情況一定程度上可以表征人類活動的強度。地方公共財政收入、地方公共財政支出、人口、農(nóng)作物總播種面積等社會經(jīng)濟發(fā)展指標與NDVImax的相關系數(shù)，可以表示植被與人類活動的相關性（表4）。地方公共財政收入和地方公共財政支出反映了區(qū)域經(jīng)濟發(fā)展水平，能夠在一定程度上表征區(qū)域?qū)Σ焕麣夂蛞蛩赜绊懙倪m應補償能力。區(qū)域經(jīng)濟向好使得其在氣候變化過程中具有較好的適應能力，且可能在城市綠化、農(nóng)田水利設施等生態(tài)建設方面有更高的投入，因此，NDVImax與地方公共財政收入和支出呈極顯著正相關關系，且相關系數(shù)均在0.70以上。

1982年以來，寧夏先后組織實施4次政策性移民搬遷?！笆濉逼陂g重點實施的中南部生態(tài)移民攻堅工程，將中南部山區(qū)貧困程度深、生存條件差、發(fā)展難度大、不適宜人類生存地區(qū)生活最困難的群眾，進行生態(tài)搬遷安置。固原市是移民重點遷出區(qū)。2000年以來，固原市人口極顯著減少，人口密度降低，加上全面實施退耕還林、圍欄禁牧政策，固原農(nóng)作物播種總面積極顯著減小，草地林地面積增加，NDVImax呈極顯著增大趨勢，NDVImax與人口、農(nóng)作物播種總面積呈極顯著負相關關系。其他市人口極顯著增多，NDVImax與人口成極顯著正相關，與吳忠市農(nóng)作物總播種面積無顯著相關關系，與其他市總播種面積顯著正相關。

3 討論與結(jié)論

3.1 討論

本研究旨在明確寧夏歸一化植被指數(shù)變化特征，了解植被與氣候及人類活動的相關性。寧夏NDVImax總體上增長，且空間上呈現(xiàn)南北高、中間低的分布特征，該結(jié)果與中國[9，15]、西北地區(qū)[12]、黃土高原[16-17]、寧夏[18]、賀蘭山[19]等區(qū)域的相關研究結(jié)論基本一致。本研究重點針對寧夏3區(qū)5市NDVImax變化特征進行對比分析，提高了精細化程度。研究顯示寧夏大部分地區(qū)NDVImax具有弱持續(xù)性，與趙安周等針對黃土高原的研究[20]基本一致。相比于平均氣溫，降水量與寧夏NDVImax的相關性更顯著，即水分因子是寧夏植被生長的主要限制因子，該研究結(jié)果與黃悅悅等的研究結(jié)論[18，21]基本一致，所用衛(wèi)星數(shù)據(jù)分辨率更高、氣象站點數(shù)更多，分析精度得到提升。近19年植被有所退化的區(qū)域主要集中在部分市（縣、區(qū)）城區(qū)，可能位于城市擴張過程中土地利用類型轉(zhuǎn)變的地區(qū)。有研究表明，人口城鎮(zhèn)化發(fā)展、人口素質(zhì)提高等能夠減緩人口增加對植被覆蓋的負面影響[22]。寧夏大部分地區(qū)NDVImax與人口以及農(nóng)作物總播種面積呈正相關關系，可能得益于人口城鎮(zhèn)化發(fā)展、綠化技術進步、環(huán)保意識提高、農(nóng)業(yè)技術及灌溉設施改善等因素，但研究時段較短，具體原因有待進一步探討。

當前寧夏NDVImax變化普遍具有弱持續(xù)性，存在生態(tài)退化風險，需要繼續(xù)依托生態(tài)工程的大力實施，提高NDVI變化穩(wěn)定性，保護生態(tài)建設成果。同時，需要在城市化建設中，合理規(guī)劃，在土地利用類型轉(zhuǎn)換的同時，提高城市綠化覆蓋率，與人工增雨、優(yōu)化水資源利用、生態(tài)工程、生態(tài)補償相結(jié)合，開展綜合生態(tài)治理。

3.2 結(jié)論

本研究基于MOD13Q1 NDVI數(shù)據(jù)、氣象數(shù)據(jù)、社會經(jīng)濟數(shù)據(jù)等，分析了寧夏NDVImax時空變化特征，研究了NDVImax與氣候、人類活動的相關關系。研究結(jié)果表明：（1）寧夏多年平均NDVImax呈現(xiàn)出南北高、中間低的分布特征。南部山區(qū)植被好于灌區(qū)、中部干旱帶。2000—2018年，寧夏平均NDVImax極顯著增大，增長速率為0.08/年（P<0.01）;均值為0.384。（2）19年間，寧夏70.3%的地區(qū)NDVImax顯著增大。植被有所退化的區(qū)域主要集中在銀川市、大武口區(qū)、利通區(qū)、沙坡頭區(qū)和中寧縣城區(qū)?；谙裨叨菻urst指數(shù)和NDVImax線性傾向率分析，寧夏未來27.8%的區(qū)域植被將顯著改善。NDVImax持續(xù)性較弱。（3）全區(qū)大部分地區(qū)NDVImax隨夏季、年降水量增加而增大，與秋季降水相關性較弱;NDVImax與春夏秋季、年平均氣溫顯著正相關面積較小。降水與寧夏NDVImax相關性更強，主要為正相關。NDVImax與地方公共財政收入和地方公共財政支出呈極顯著正相關關系;大部分地區(qū)NDVImax與人口以及農(nóng)作物總播種面積呈正相關關系，固原市NDVImax與二者呈極顯著負相關關系。

生態(tài)工程是人為影響植被生長的重要驅(qū)動因素[23-25]，如退耕還林還草、封山禁牧等生態(tài)工程通過改善土壤質(zhì)量[26]、水土流失[27]、草原植物群落結(jié)構(gòu)[28]及提高產(chǎn)草量[29]、植被總初級生產(chǎn)力[30]等方式使植被得到恢復，NDVI增大。由于數(shù)據(jù)限制，文中未作深入探討，需結(jié)合造林面積、綠化面積等數(shù)據(jù)進一步探討生態(tài)工程對植被的影響。

參考文獻：

[1]朱永森，曾永年，張 猛. 基于HJ衛(wèi)星數(shù)據(jù)與面向?qū)ο蠓诸惖耐恋乩?覆蓋信息提取[J]. 農(nóng)業(yè)工程學報，2017，33（14）：258-265.

[2]Wardlow B D，Egbert S L，Kastens J H. Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S.Central Great Plains[J]. Remote Sensing of Environment，2007，108（3）：290-310.

[3]朱文泉，潘耀忠，張錦水. 中國陸地植被凈初級生產(chǎn)力遙感估算[J]. 植物生態(tài)學報，2007，31（3）：413-424.

[4]郭仲偉，吳朝陽，汪簫悅. 基于衛(wèi)星遙感數(shù)據(jù)的森林病蟲害監(jiān)測與評價[J]. 地理研究，2019，38（4）：831-843.

[5]Mkhabela M S，Bullock P，Raj S，et al. Crop yield forecasting on the Canadian Prairies using MODIS NDVI data[J]. Agricultural and Forest Meteorology，2011，151（3）：385-393.

[6]張學霞，葛全勝，鄭景云. 近50年北京植被對全球變暖的響應及其時效——基于遙感數(shù)據(jù)和物候資料的分析[J]. 生態(tài)學雜志，2005，24（2）：123-130.

[7]陳曉光，李劍萍，韓穎娟，等. 寧夏近20年來植被覆蓋度及其與氣溫降水的關系[J]. 生態(tài)學雜志，2007，26（9）：1375-1383.

[8]趙卓文，張連蓬，李 行，等. 基于MOD13Q1數(shù)據(jù)的寧夏生長季植被動態(tài)監(jiān)測[J]. 地理科學進展，2017，36（6）：741-752.

[9]劉 可，杜靈通，侯 靜，等. 近30年中國陸地生態(tài)系統(tǒng)NDVI時空變化特征[J]. 生態(tài)學報，2018，38（6）：1885-1896.

[10]杜靈通，田慶久. 寧夏回族自治區(qū)NDVI的時空變化特征研究[J]. 水土保持通報，2011，31（5）：208-214.

[11]杜靈通，田慶久. 寧夏植被覆蓋動態(tài)變化及與氣候因子的關系[J]. 中國沙漠，2012，32（5）：1479-1485.

[12]周夢甜，李 軍，朱康文. 西北地區(qū)NDVI變化與氣候因子的響應關系研究[J]. 水土保持研究，2015，22（3）：182-187.

[13]魏鳳英. 現(xiàn)代氣候統(tǒng)計診斷與預測技術[M]. 北京：氣象出版社，1999.

[14]陳彥光. 基于Matlab的地理數(shù)據(jù)分析[M]. 北京：高等教育出版社，2012.

[15]劉憲鋒，朱秀芳，潘耀忠，等. 1982—2012年中國植被覆蓋時空變化特征[J]. 生態(tài)學報，2015，35（16）：5331-5342.

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

[17]Ren Y J，Lü Y H，F(xiàn)u B J，et al. Driving factors of land change in Chinas Loess Plateau：quantification using geographically weighted regression and management implications[J]. Remote Sensing，2020，12（3）：453.

[18]黃悅悅，楊 東，馮 磊. 2000—2016年寧夏植被覆蓋度的時空變化及其驅(qū)動力[J]. 生態(tài)學雜志，2019，38（8）：2515-2523.

[19]朱 源，王 志，彭光雄.基于SPOT-NDVI的賀蘭山植被覆蓋的時空變化研究[J]. 遙感技術與應用，2012，27（1）：142-148.

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

[21]李 欣，王連喜，李 琪，等. 寧夏近25年植被指數(shù)變化及其與氣候的關系[J]. 干旱區(qū)資源與環(huán)境，2011，25（9）：161-166.

[22]李 超. 人口總量與結(jié)構(gòu)變化影響植被覆蓋的實證研究[J]. 廣東農(nóng)業(yè)科學，2015，42（12）：186-192.

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

[24]史娜娜，肖能文，王 琦，等. 錫林郭勒植被NDVI時空變化及其驅(qū)動力定量分析[J]. 植物生態(tài)學報，2019，43（4）：331-341.

[25]錢 琛，陳海濱，侯現(xiàn)慧.生態(tài)建設背景下陜西省植被覆蓋時空變化及其影響因素研究[J]. 干旱區(qū)地理，2020，43（2）：425-433.

[26]彭文英，張科利，陳 瑤，等. 黃土坡耕地退耕還林后土壤性質(zhì)變化研究[J]. 自然資源學報，2005，20（2）：272-278.

[27]王怡菲，姚順波，鄧元杰. 渭河流域水土流失治理效率的時空格局演化與影響因素[J]. 地理科學，2019，39（5）：836-846.

[28]高 露，張圣微，朱仲元，等. 放牧對干旱半干旱草原植物群落結(jié)構(gòu)和生態(tài)功能的影響[J]. 水土保持研究，2019，26（6）：205-211.

[29]趙 鋼，曹子龍，李青豐. 春季禁牧對內(nèi)蒙古草原植被的影響[J]. 草地學報，2003，11（2）：183-188.

[30]尤南山，董金瑋，肖 桐，等. 退耕還林還草工程對黃土高原植被總初級生產(chǎn)力的影響[J]. 地理科學，2020，40（2）：315-323.