摘要：通過(guò)研究賀蘭山植被指數(shù)的時(shí)空變化特征，定量評(píng)估賀蘭山植被的影響因素，對(duì)科學(xué)評(píng)價(jià)賀蘭山生態(tài)治理工程效益和區(qū)域生態(tài)系統(tǒng)恢復(fù)及其可持續(xù)發(fā)展提供重要的理論和指導(dǎo)?；谘芯繀^(qū)NDVI數(shù)據(jù)、年均氣溫、降水量數(shù)據(jù)及其他數(shù)據(jù)，采用趨勢(shì)分析和皮爾遜（Pearson）相關(guān)分析、多元回歸殘差分析方法，從時(shí)間和空間尺度研究了賀蘭山2000 — 2020年NDVI的變化特征，分析了影響NDVI變化的主要驅(qū)動(dòng)力。結(jié)果表明，21 a間，賀蘭山NDVI整體呈現(xiàn)波動(dòng)上升趨勢(shì)，NDVI改善趨勢(shì)明顯；空間分布整體呈西南高東北低的變化特征，賀蘭山中部地區(qū)退化趨勢(shì)明顯；降水增加、生態(tài)環(huán)境保護(hù)政策的支持是賀蘭山NDVI變化的主要驅(qū)動(dòng)因素，21 a來(lái)NDVI變化與年均氣溫、累計(jì)降水量變化均呈單峰型變化趨勢(shì)。整體來(lái)看，NDVI與氣溫、累計(jì)降水均呈波動(dòng)上升的變化趨勢(shì)。

關(guān)鍵詞：賀蘭山；NDVI；時(shí)空變化特征

中圖分類號(hào)：TU984.113" " " " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2097-2172（2025）02-0171-06

doi：10.3969/j.issn.2097-2172.2025.02.013

Spatio-temporal Changes and Driving Forces of NDVI in

Helan Mountains from 2000 to 2020

MA Hai 1， ZHANG Qiang 2， LI Jinyu 3， LU Jiaxin 1， LI Yujia 1， CHEN Xiaolan 1， CHEN Zhirong 1

（1. College of Resources， Environment and Life Sciences， Ningxia Normal University， Guyuan Ningxia 756000， China;

2. Secondary School Affiliated to Ningxia Normal University， Guyuan Ningxia 756000， China;

3. Jingyuan County Weather Bureau， Jingyuan Ningxia 756400， China）

Abstract： This study investigated the spatio-temporal variation characteristics of vegetation index （NDVI） in the Helan Mountains， providing a quantitative evaluation of the factors affecting vegetation and offering significant theoretical and practical guidance for evaluating the effectiveness of ecological restoration projects and sustainable development in the region. Using NDVI data， annual average temperature， precipitation data， and other data from the study area， trend analysis， Pearson correlation analysis， and multiple regression residual analysis were employed to study the changes in NDVI from 2000 to 2020 in the Helan Mountains， examining the main driving forces behind the NDVI variations on both temporal and spatial scales. Results showed that over the 21 years， the NDVI in the Helan Mountains exhibited a fluctuating upward trend， with a significant improvement in NDVI. The spatial distribution showed a general pattern of higher values in the southwest and lower values in the northeast， with noticeable degradation in the central part of the Helan Mountains. Increased precipitation and the support of ecological environment protection policies were the main driving factors for NDVI changes in the Helan Mountains. Over the 21 years， NDVI variations showed a unimodal trend with annual average temperature and cumulative precipitation changes. Overall， both NDVI and the annual average temperature， as well as cumulative precipitation， exhibited a fluctuating upward trend.

Key words： Helan Mountain; NDVI; Characteristic of spatio-temporal change

植被覆蓋度是描述陸地生態(tài)系統(tǒng)中植被生長(zhǎng)和變化的重要參數(shù)，其中歸一化植被指數(shù)（Normalized Difference Vegetation Index，NDVI）被廣泛用于植被覆蓋度評(píng)價(jià)中，能夠定量評(píng)價(jià)植被覆蓋度狀況，進(jìn)而反映生態(tài)安全的穩(wěn)定性，對(duì)于評(píng)估生態(tài)安全的穩(wěn)定性等具有重要的指導(dǎo)價(jià)值。習(xí)近平總書記視察寧夏時(shí)指出，寧夏要努力建設(shè)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展先行區(qū)，加強(qiáng)賀蘭山、六盤山、羅山自然保護(hù)區(qū)建設(shè)，統(tǒng)籌生態(tài)保護(hù)修復(fù)和環(huán)境綜合治理，賀蘭山是我國(guó)重要自然地理分界線和西北重要生態(tài)安全屏障。要加強(qiáng)頂層設(shè)計(jì)，狠抓責(zé)任落實(shí)，強(qiáng)化監(jiān)督檢查，堅(jiān)決保護(hù)好賀蘭山生態(tài)［1 ］。寧夏回族自治區(qū)也提出，要實(shí)行最嚴(yán)格的生態(tài)環(huán)境保護(hù)制度，協(xié)同推進(jìn)降碳、節(jié)水、減污、擴(kuò)綠、增長(zhǎng)，全面提升資源生態(tài)系統(tǒng)穩(wěn)定性和生態(tài)服務(wù)功能，努力建設(shè)天藍(lán)、地綠、水美的美麗寧夏［2 ］。

賀蘭山脈位于半干旱區(qū)向干旱區(qū)的過(guò)渡地帶，周邊地區(qū)皆為荒漠與半荒漠，既是內(nèi)蒙古和寧夏兩地的生物資源寶庫(kù)之一，也是生物多樣性保護(hù)的重點(diǎn)區(qū)域［3 ］，具有典型的溫帶干旱區(qū)植被組合，共有常綠針葉林、針闊葉混交林、夏綠闊葉林、常綠針葉灌叢、落葉闊葉灌叢、疏林草原、典型草原、荒漠草原、荒漠、草甸和農(nóng)田等11個(gè)植被型［4 ］。近年來(lái)，賀蘭山的生態(tài)環(huán)境惡化趨勢(shì)基本得到逆轉(zhuǎn)，但生態(tài)環(huán)境脆弱性依然未能從根本上改變。當(dāng)前生態(tài)脆弱區(qū)面臨著既要保護(hù)綠水青山又要?jiǎng)?chuàng)造金山銀山的雙重任務(wù)，生態(tài)保護(hù)與經(jīng)濟(jì)協(xié)同發(fā)展的問(wèn)題日益突出。因此，研究賀蘭山植被指數(shù)的時(shí)空變化特征，定量評(píng)估賀蘭山植被的影響因素，對(duì)科學(xué)評(píng)價(jià)賀蘭山生態(tài)治理工程效益和區(qū)域生態(tài)系統(tǒng)恢復(fù)及其可持續(xù)發(fā)展具有重要的理論和現(xiàn)實(shí)指導(dǎo)意義。

1" "數(shù)據(jù)來(lái)源與方法

1.1" "研究區(qū)概況

賀蘭山處于青藏高原、蒙古高原和黃土高原的交界處，山脈呈南北走向綿延250 km有余，東西寬20～40 km，是我國(guó)重要自然地理分界線和西北重要生態(tài)安全屏障［5 ］。賀蘭山以其獨(dú)特的生態(tài)系統(tǒng)孕育有豐富的自然資源，動(dòng)植物和地下礦藏資源十分豐富，有野生維管植物788種，特有與近特植物60多種；陸生脊椎動(dòng)物352種，其中國(guó)家重點(diǎn)保護(hù)脊椎動(dòng)物有58種。寧夏賀蘭山被稱為當(dāng)?shù)氐摹案赣H山”［4 ］，對(duì)當(dāng)?shù)鼐哂猩鷳B(tài)以及經(jīng)濟(jì)等重要意義。

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

本文應(yīng)用的遙感數(shù)據(jù)采自地理空間數(shù)據(jù)云（https：//ladsweb.nascom.nasa.gov/search/）下載的標(biāo)準(zhǔn)化250 m MOD13A1植被指數(shù)。2000 — 2020年研究區(qū)MOD13A1數(shù)據(jù)產(chǎn)品使用ENVI軟件進(jìn)行輻射校正［6 - 7 ］，將NDVI影響的投影坐標(biāo)與研究區(qū)的投影坐標(biāo)系統(tǒng)相匹配。將預(yù)處理后的遙感影像圖經(jīng)過(guò)2%的拉伸，增強(qiáng)圖像的顯示效果。生成16 d最大化合成NDVI數(shù)據(jù)，投影為Albers（等面積割圓錐投影），WGS84作為大地基準(zhǔn)。應(yīng)用賀蘭山矢量邊界進(jìn)行裁剪，得到研究區(qū)NDVI數(shù)據(jù)［8 ］。氣象數(shù)據(jù)來(lái)自Zenodo網(wǎng)站，該數(shù)據(jù)是新的2000 — 2020年中國(guó)高分辨率（1 km）月度網(wǎng)格化平均氣溫?cái)?shù)據(jù)集和降水量數(shù)據(jù)集，是基于氣象站數(shù)據(jù)的高斯過(guò)程回歸（GPR）方法制作而成的月度柵格數(shù)據(jù)，并通過(guò)氣象站觀測(cè)數(shù)據(jù)進(jìn)行精確度評(píng)估［9 ］，能夠很好地反映地區(qū)的氣溫和降水量時(shí)空分布情況，數(shù)據(jù)坐標(biāo)系為WGS84，文件格式為TIF，柵格單位為℃；降水量單位為mm［10 ］。

1.3" "研究方法

1.3.1" " 線性趨勢(shì)分析" " 賀蘭山2000 — 2020年的植被覆蓋度旨在分析自變量植被覆蓋度隨時(shí)間變量變化的關(guān)系，經(jīng)檢驗(yàn)二者的關(guān)系可近似用一條直線表示，滿足一元線性回歸模型的條件，同時(shí)一元線性回歸可以在空間上模擬每個(gè)柵格的變化趨勢(shì)［11 ］，以單個(gè)像元植被覆蓋度隨時(shí)間變化特征反映整個(gè)空間的變化規(guī)律。其優(yōu)點(diǎn)在于利用不同時(shí)段數(shù)據(jù)值的擬合，消除異常因素對(duì)植被覆蓋的影響，真實(shí)地反映植被覆蓋的演化趨勢(shì)。Slope 為單個(gè)像元 NDVI 回歸方程的斜率反映了在21 a的時(shí)間序列中，賀蘭山地區(qū) NDVI的變化趨勢(shì)。其計(jì)算公式如下［12 ］。

Slope=

式中，NDVIi為第i年的 NDVI值；n為研究時(shí)段跨度（2000 — 2020年共21 a）。若Slope為正值時(shí)，表明隨時(shí)間變化 NDVI升高，區(qū)域植被狀況趨于正向發(fā)展趨勢(shì)且值越大正向趨勢(shì)越明顯；反之則為退化趨勢(shì)。

1.3.2" " 皮爾遜（Pearson）相關(guān)分析" " 主要運(yùn)用Pearson相關(guān)分析來(lái)研究影響因素（氣溫、降水）在賀蘭山與NDVI的相關(guān)性，Pearson相關(guān)系數(shù)是反映等級(jí)相關(guān)程度的統(tǒng)計(jì)分析指標(biāo)［13 ］，是最常用的一種相關(guān)系數(shù)［14 ］。它的值域?yàn)椋?1～1］，其中，相關(guān)程度的強(qiáng)弱可以通過(guò)絕對(duì)值的大小來(lái)描述。將賀蘭山NDVI與驅(qū)動(dòng)因子的相關(guān)性劃分為3個(gè)類型［15 ］：不顯著正相關(guān)、顯著正相關(guān)和極顯著正相關(guān)（表1）。

1.3.3" " 多元回歸殘差分析" " 采用多元回歸殘差分析方法研究人類活動(dòng)（主要為生態(tài)修復(fù)）和氣候變化對(duì)植被NDVI變化的影響及相對(duì)貢獻(xiàn)［16 - 17 ］。該方法主要有以下3個(gè)步驟：基于多年NDVI以及插值后的氣溫和降水量時(shí)間序列數(shù)據(jù)，以NDVI為因變量、以氣溫和降水量為自變量，建立二元線性回歸模型，計(jì)算模型中的各項(xiàng)參數(shù)；基于氣溫、降水量數(shù)據(jù)以及回歸模型的參數(shù)，計(jì)算得到NDVI的預(yù)測(cè)值（NDVICC），用來(lái)表示氣候因素對(duì)植被NDVI的影響；計(jì)算NDVI觀測(cè)值與NDVICC之間的差值，即NDVI殘差（NDVIHA），用來(lái)表示人類活動(dòng)對(duì)植被NDVI的影響［18 ］。計(jì)算公式如下。

NDVICC = a×T+b×P+c

NDVIHA = NDVIobs-NDVICC

式中，NDVICC和NDVIobs分別為基于回歸模型的NDVI預(yù)測(cè)值和基于遙感影像的NDVI觀測(cè)值（無(wú)量綱）；a、b和c為模型參數(shù)；T和P分別為生長(zhǎng)季平均氣溫和累計(jì)降水量；NDVIHA）為殘差。

2" "結(jié)果與分析

2.1" "賀蘭山NDVI時(shí)間變化特征

由圖1可知，從時(shí)間變化分析得出，研究區(qū)內(nèi)NDVI 21 a均值呈現(xiàn)出波浪式增長(zhǎng)，范圍為0.248 2～0.378 4，漲幅約為52.5%，其中2000年NDVI為最小值，2018年為最大值。NDVI高值一直集中在西南，隨著時(shí)間的變化逐漸向東擴(kuò)散（圖2）。通過(guò)水平為0.262 5的顯著性檢驗(yàn)，說(shuō)明賀蘭山NDVI的總體趨勢(shì)在時(shí)間尺度上變化顯著。

2.2" "賀蘭山NDVI空間變化特征

將賀蘭山2000 — 2020年NDVI數(shù)據(jù)在ArcMap進(jìn)行Slope趨勢(shì)分析并進(jìn)行重采樣得到賀蘭山NDVI變化趨勢(shì)圖［19 ］，再將變化趨勢(shì)進(jìn)行分類，分為三個(gè)等級(jí)，分別為植被退化區(qū)（Slope lt; 0）、植被不明顯改善區(qū)（0 lt; Slope≤0.01）、植被明顯改善區(qū)（Slope ＞ 0.01）。由圖3可知，賀蘭山植被覆蓋區(qū)2000 — 2020年NDVI呈現(xiàn)改善趨勢(shì)，植被覆蓋面積占比較大的為植被不明顯改善區(qū)域；植被較明顯改善區(qū)域面積占比較小；其中植被退化區(qū)域大面積分布在賀蘭山中部地區(qū)；對(duì)各個(gè)NDVI變化等級(jí)逐像元進(jìn)行統(tǒng)計(jì)分析，3個(gè)等級(jí)像元數(shù)量比例為5 041∶29 604∶435。

由圖4可知，賀蘭山西南地區(qū)NDVI偏高，東北地區(qū)NDVI明顯偏低，空間分布差異明顯。將研究區(qū)的NDVI劃分為5個(gè)等級(jí)，分別為NDVI lt; 0.2、0.2 ≤NDVI ≤0.4、0.4 lt; NDVI ≤0.6、0.6 lt; NDVI ≤0.8、NDVI" gt; 0.8，對(duì)NDVI 等級(jí)逐像元進(jìn)行統(tǒng)計(jì)分析，5個(gè)等級(jí)像元數(shù)量比例為8 025∶19 425∶" " " " "5 430∶2 159∶40，NDVI" lt; 0.2區(qū)域主要分布于賀蘭山南北兩端，向西南方向NDVI 逐漸增高。

2.3" "賀蘭山NDVI變化特征分析

2.3.1" " NDVI與氣候因素相關(guān)性分析" " 氣溫和降水是影響植被生長(zhǎng)最重要的兩個(gè)氣候因素，植被又通過(guò)生理作用影響著區(qū)域氣候的變化［20 - 21 ］。氣溫通過(guò)影響植被的光合作用、呼吸作用、蒸騰作用、礦物質(zhì)吸收來(lái)改變植被覆蓋度，而降水則通過(guò)影響植被的組成成分、蒸散發(fā)影響植被生長(zhǎng)情況。由圖5可知，賀蘭山2000 — 2020年研究區(qū)年平均氣溫在6.368 2～7.691 3 ℃波動(dòng)增長(zhǎng)，整體變化不大，2017年年平均氣溫最高，2012年年平均氣溫最低。賀蘭山21 a NDVI與平均氣溫變化相似，均呈現(xiàn)波動(dòng)增長(zhǎng)的趨勢(shì)，說(shuō)明具有一定相關(guān)性。對(duì)21 a NDVI與氣溫相關(guān)性進(jìn)行分析可知，平均氣溫和NDVI均值之間的相關(guān)數(shù)值為0.129 ，說(shuō)明賀蘭山平均氣溫和NDVI均值之間有著不顯著的正相關(guān)關(guān)系，所以氣溫升高可在一定程度上促進(jìn)光合作用和提高水分利用效率［22 ］，從而有利于植被生長(zhǎng)，但是這種變化并不顯著。

由圖6可知，賀蘭山21 a NDVI與累計(jì)降水變化較為相似，均呈現(xiàn)出波動(dòng)增長(zhǎng)的趨勢(shì)，說(shuō)明具有一定相關(guān)性。賀蘭山2000 — 2020年研究區(qū)年累計(jì)降水在154.47～288.42 mm波動(dòng)增長(zhǎng)。其中2018年年累計(jì)降水最高，2000年年累計(jì)降水量最低。對(duì)賀蘭山21 a NDVI與降水相關(guān)性進(jìn)行分析可知，累計(jì)降水量和NDVI均值之間的相關(guān)系數(shù)值為0.525，說(shuō)明累計(jì)降水量和NDVI均值之間有著顯著的正相關(guān)關(guān)系，所以降水對(duì)于植被狀況起重要作用。

2.3.2" " NDVI與人類活動(dòng)相關(guān)性分析" " 由圖7、表2可知，人類活動(dòng)（主要為生態(tài)修復(fù)）對(duì)賀蘭山植被NDVI變化的貢獻(xiàn)率為正向作用，其中，人類活動(dòng)的貢獻(xiàn)率在0%～20%、20%～40%、40%～60%范圍的區(qū)域面積較小，像元占比分別為3.4%、7.7%、19.2%；60%～80%、80%～100%范圍的區(qū)域面積較大，像元占比分別為35.1%、34.6%；賀蘭山全域以人類活動(dòng)貢獻(xiàn)率為主，其中貢獻(xiàn)率超過(guò)80%的區(qū)域主要集中在賀蘭山中部地區(qū)，說(shuō)明在賀蘭山植被NDVI變化過(guò)程中人類活動(dòng)起決定作用。

3" "討論與結(jié)論

寧夏賀蘭山21 a來(lái)，NDVI變化整體呈現(xiàn)波浪式增長(zhǎng)趨勢(shì)，從2000、2010、2020年等每10 a最大化合成的NDVI值，高值一直集中在西南，隨著時(shí)間的變化逐漸向東擴(kuò)散，空間分布特征的原因是：賀蘭山東部常年受中亞季風(fēng)的影響，西段受盛行西風(fēng)的影響，水汽的來(lái)源有所不同，降水量由東北向西南逐漸遞減，從年際變化方面分析，近年來(lái)賀蘭山及周邊地區(qū)氣溫升高，降水量增加，從而使氣候由暖干向暖濕轉(zhuǎn)變［23 ］，NDVI變化與氣溫、降水整體均呈波動(dòng)上升的趨勢(shì)，這與李婷婷等［24 ］、劉錟等［25 ］、朱源等［26 ］得出的結(jié)論基本一致。賀蘭山因其地理位置和氣候條件的獨(dú)特性，每年1 — 4月和10 — 12月NDVI基本無(wú)變化，因此在本文分析中不作考慮。

本文以寧夏賀蘭山為研究區(qū)，基于遙感影像數(shù)據(jù)以及氣象數(shù)據(jù)，綜合運(yùn)用MRT中的處理工具和ArcGIS實(shí)現(xiàn)對(duì)賀蘭山NDVI差異分析。采用趨勢(shì)線分析方法，分析了賀蘭山21 a的NDVI的時(shí)空變化特征以及相關(guān)系數(shù)分析了植被指數(shù)與驅(qū)動(dòng)力的相關(guān)性，結(jié)果表明，在時(shí)間尺度上，研究區(qū)21 a來(lái)NDVI呈波浪式增長(zhǎng)，NDVI變化與時(shí)間存在顯著相關(guān)性，從2000年的0.248 2增長(zhǎng)到2018年的0.378 4，植被覆蓋狀況逐年改善；在空間尺度上，賀蘭山植被指數(shù)存在明顯的區(qū)域性差異，呈現(xiàn)出西南高、東北低的特點(diǎn)，賀蘭山中部地區(qū)退化趨勢(shì)明顯；氣候變化和人類活動(dòng)的共同作用是21 a來(lái)賀蘭山植被NDVI變化的主要驅(qū)動(dòng)因素。賀蘭山21 a NDVI的變化與平均氣溫、累計(jì)降水變化比較相似，均呈單峰型變化趨勢(shì)，均呈波動(dòng)上升的變化趨勢(shì)；與氣溫的相關(guān)關(guān)系類型為不顯著正相關(guān)，與降水量的相關(guān)關(guān)系類型為顯著正相關(guān)。人類活動(dòng)對(duì)賀蘭山植被NDVI影響呈正向作用。由此得知，人類活動(dòng)（主要為生態(tài)修復(fù)）對(duì)賀蘭山植被NDVI的作用不可忽視。

參考文獻(xiàn)：

［1］ 郇慶治." 論習(xí)近平生態(tài)文明思想的形成與發(fā)展［J］." 鄱陽(yáng)湖學(xué)刊，2022（4）：5-28；124.

［2］ 李" "錦." 山水寧夏，繁星映雪燦金山［N］." 寧夏日?qǐng)?bào)，2023-01-07（003）.

［3］ 張立漢." 中國(guó)山河全書（上）［M］." 青島：青島出版社，2005.

［4］ 衛(wèi)璽璽，晁鑫艷，鄭景明，等." 賀蘭山東、西側(cè)典型植物群落物種多樣性差異及其影響因子［J］." 生態(tài)環(huán)境學(xué)報(bào)，2024，33（4）：520-530.

［5］ 陳衛(wèi)平." 賀蘭山—銀川盆地景觀格局分析與景觀規(guī)劃［D］." 北京：北京林業(yè)大學(xué)， 2008.

［6］ 張曉濤." 區(qū)域蒸發(fā)蒸騰量的遙感估算［D］." 楊凌：西北農(nóng)林科技大學(xué)，2006.

［7］ 徐根生，何應(yīng)法，付" "煒，等." 基于MOD13Q1的永州近15年植被變化研究［J］." 甘肅農(nóng)業(yè)科技，2015（3）：27-29.

［8］ 魯" "暉." 2000 — 2017年河西地區(qū)植被覆蓋的時(shí)空過(guò)程分析［D］." 蘭州：蘭州大學(xué)，2019.

［9］ 李春艷，張" "慧，于永奎，等." 2022年平羅縣氣象因子對(duì)主要農(nóng)作物生長(zhǎng)發(fā)育的影響［J］." 寒旱農(nóng)業(yè)科學(xué)，2024，3（3）：245-250.

［10］ CARRI？魷N V J， PEREZ-JARAMILLO J， CORDOVEZ V， et al. Pathogen-induced activation of disease-suppressive functions in the endophytic root microbiome［J］. Science. 2019， 366（6465）： 606-612.

［11］ 雷聲劍，張福平，燕玉超，等." 黑河上游植被覆蓋時(shí)空變化特征及其未來(lái)趨勢(shì)［J］." 水土保持通報(bào)，2016，36（3）：159-164；370.

［12］ 馬明國(guó)， 董立新，王雪梅." 過(guò)去21 a中國(guó)西北植被覆蓋動(dòng)態(tài)監(jiān)測(cè)與模擬［J］." 冰川凍土，2003，25（2）：232-236.

［13］ 王曉媛，孫" "嬌." 有機(jī)肥替代氮肥對(duì)銀北灌區(qū)鹽化灌淤土理化性質(zhì)及玉米生長(zhǎng)的影響［J］." 寒旱農(nóng)業(yè)科學(xué)，2024，3（1）：51-56.

［14］ 劉昕怡." 1998 — 2018年桂林市植被變化及對(duì)相關(guān)因子分析［J］." 南寧師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，39（3）：38-46.

［15］ 周" "亮." 基于InVEST模型評(píng)估賀蘭山土地系統(tǒng)碳儲(chǔ)量及生境質(zhì)量［D］." 銀川：寧夏大學(xué)，2022.

［16］ EVANS J， GEERKEN R. Discrimination between climate and human-induced dryland degradation［J］." Journal of Arid Environments， 2004， 57（4）： 535-554.

［17］ WESSELS K J， PRINCE S D， MALHERBE J， et al. Can human-induced land degradation be distinguished from the effects ofrainfall variability？ A case study in South Africa［J］." Journal of Arid Environments， 2007， 68（2）： 271-297.

［18］ 金" "凱，王" "飛，韓劍橋，等." 1982 — 2015年中國(guó)氣候變化和人類活動(dòng)對(duì)植被NDVI變化的影響［J］." 地理學(xué)報(bào)，2020，75（5）：961-974.

［19］ 余" "琳，姚" "昆，劉漢湖，等." 不同高程下的NDVI變化及其與氣候因子的關(guān)系［J］." 物探化探計(jì)算技術(shù)， 2017，39（2）：296-300.

［20］ 仝莉棉." 我國(guó)中緯度夏季風(fēng)邊緣帶植被時(shí)空變化及其影響因素［D］." 蘭州：蘭州大學(xué)，2017.

［21］ 祁玉娟，王曉華." 蘭州市紅古區(qū)草原植被監(jiān)測(cè)研究［J］." 甘肅農(nóng)業(yè)科技，2022，53（8）：75-80.

［22］ 龐國(guó)偉，山琳昕，楊勤科，等." 陜西省不同地貌類型區(qū)植被覆蓋度時(shí)空變化特征及其影響因素［J］." 長(zhǎng)江科學(xué)院院報(bào)，2021，38（3）：51-58；76.

［23］ 施雅風(fēng)，沈永平，李棟梁，等." 中國(guó)西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢(shì)探討［J］." 第四紀(jì)研究，2003（2）：152-164.

［24］ 李婷婷，馬" "超，郭增長(zhǎng)." 基于RSEI模型的賀蘭山長(zhǎng)時(shí)序生態(tài)質(zhì)量評(píng)價(jià)及影響因素分析［J］." 生態(tài)學(xué)雜志，2021，40（4）：1154-1165.

［25］ 劉" "錟，康慕誼." 基于RS和GIS的賀蘭山高山林線初步研究［J］." 自然資源學(xué)報(bào)，2016，31（6）：973-981.

［26］ 朱" "源，王" "志，彭光雄." 基于SPOT-NDVI的賀蘭山植被覆蓋的時(shí)空變化研究［J］." 遙感技術(shù)與應(yīng)用，2012，27（1）：142-148.