地形因子影響典型氣候帶植被指數(shù)的定量歸因

摘 要：對(duì)不同氣候帶地形因子影響歸一化植被指數(shù)（NDVI）的變化進(jìn)行分析研究，對(duì)于探索水土災(zāi)害、生態(tài)研究以及區(qū)域環(huán)境質(zhì)量等具有重要價(jià)值。以我國(guó)7個(gè)典型氣候帶下的縣域范圍為例，基于Landsat 8遙感數(shù)據(jù)與DEM數(shù)據(jù)，從地理探測(cè)器與地理加權(quán)回歸視角，探索中國(guó)典型氣候帶的小區(qū)域范圍尺度下地形因子與NDVI變化的關(guān)系與規(guī)律。研究表明：（1）NDVI排名為：漠河市（北溫帶）gt;白沙黎族自治縣（中/北熱帶）gt;浦北縣（南亞熱帶）gt;魯?shù)榭h（南亞熱帶）gt;方山縣（南溫帶）gt;貴德縣（高原氣候帶）gt;巴林左旗（中溫帶）。（2）地形因子對(duì)植被指數(shù)的解釋力排名為：坡向因子lt;坡度因子lt;高程因子，地形因子之間大部分不存在顯著性差異，不同研究區(qū)的地形因子對(duì)NDVI的影響并非簡(jiǎn)單獨(dú)立，而是通過(guò)非線(xiàn)性增強(qiáng)的過(guò)程來(lái)起作用。（3）地形因子的局部擬合空間格局與植被覆蓋空間分布格局相似，地形對(duì)植被的空間分布有著重要影響。（4）根據(jù)調(diào)整擬合優(yōu)度，各研究區(qū)地形因子與植被的擬合度排名為：貴德縣gt;方山縣gt;白沙黎族自治縣gt;浦北縣gt;漠河市g(shù)t;魯?shù)榭hgt;巴林左旗，地形因子對(duì)植被指數(shù)擬合程度隨地域變化存在顯著差異。

關(guān)鍵詞：地理探測(cè)器；地理加權(quán)回歸；地形因子；歸一化植被指數(shù)；氣候帶

中圖分類(lèi)號(hào)：P931；Q948.2"" 文獻(xiàn)標(biāo)志碼：A""" 文章編號(hào)：1673-5072（2024）05-0457-09

近年來(lái)，在全球氣候變暖的大背景和人類(lèi)活動(dòng)的影響之下，干旱、暴雨、泥石流、高溫等災(zāi)害頻發(fā)，對(duì)中國(guó)的生態(tài)環(huán)境與經(jīng)濟(jì)發(fā)展帶來(lái)巨大傷害。綠色植物是人類(lèi)賴(lài)以生存的物質(zhì)，是陸地生態(tài)環(huán)境系統(tǒng)的主要組成結(jié)構(gòu)，對(duì)保護(hù)和維持全世界的生態(tài)環(huán)境與綠色可持續(xù)發(fā)展起著不可磨滅的作用［1-2］。植被覆蓋度能夠反映區(qū)域植被的生長(zhǎng)狀況，大區(qū)域范圍植被覆蓋的變化體現(xiàn)了自然氣候演變和人類(lèi)活動(dòng)對(duì)生態(tài)環(huán)境的作用［3］。國(guó)內(nèi)外對(duì)區(qū)域植被指數(shù)的研究主要有比值植被指數(shù)（RVI）、歸一化植被指數(shù)（NDVI）、優(yōu)化的比值植被指數(shù)（MSR）等［4-7］。NDVI作為反映植被生長(zhǎng)狀況的因子［8］，能有效反映年度植被生長(zhǎng)所達(dá)到的最佳狀況，且與地上生物量有較強(qiáng)的相關(guān)關(guān)系［9］，在植被覆蓋研究中被廣泛應(yīng)用。區(qū)域內(nèi)的植被生長(zhǎng)狀況受氣候、地形、土壤成分、人類(lèi)活動(dòng)等多重因素的共同作用［10-14］，這些影響因素中，地形及其空間分布格局直觀(guān)而深刻地揭示了地表構(gòu)造變形及各種地表過(guò)程的相互作用［15］，對(duì)光照、降水、土壤、溫度的空間分布具有重要影響［16］，是植被空間分布的重要影響作用之一［17］。

近年來(lái)，地形與植被的相關(guān)研究正逐漸受到國(guó)內(nèi)外專(zhuān)家的關(guān)注。地形對(duì)植被的影響存在著顯著的空間異質(zhì)性。不同區(qū)域植被變化的地形效應(yīng)存在著多樣性：李娟和龔純偉［18］研究發(fā)現(xiàn)，高原氣候帶的祁連山植被覆蓋變化在高程、坡度等地形條件下差異明顯，坡向的地形效應(yīng)不明顯；劉梁美子等［19］得出亞熱帶氣候帶的黔桂喀斯特地貌的植被覆蓋變化與坡度、高程、坡向都呈現(xiàn)出地形分異。不同地域的植被在陰、陽(yáng)坡的變化存在差異性：陳佳敏等［20］發(fā)現(xiàn)延安市陰坡的植被覆蓋、高度、物種多樣性都優(yōu)于陽(yáng)坡；龍映豪等［21］發(fā)現(xiàn)貴州開(kāi)陽(yáng)縣陰坡與陽(yáng)坡的植被覆蓋狀況隨不同地形因子的變化而變化。不同疆域的地形因子對(duì)植被的解釋力存在異質(zhì)性：銀朵朵和王艷慧［22］探究發(fā)現(xiàn)，溫帶大陸性半干旱季風(fēng)氣候帶的地形因子對(duì)植被覆蓋率影響顯著，植被覆蓋度隨海拔增加而逐漸增高，而植被覆蓋度在坡度的分異明顯；李夢(mèng)華等［23］研究發(fā)現(xiàn)，寧夏植被覆蓋度空間分布受海拔的影響程度高于坡度；邵嘉豪等［24］研究認(rèn)為，山西省的地形中坡度因子的解釋力高于高程因子。由此可見(jiàn)，研究區(qū)域不同、空間尺度不同，地形因子對(duì)植被空間分異存在一定的規(guī)律性與差異性［25］。

目前，地形因子對(duì)植被的歸因研究較多且趨于成熟，但現(xiàn)有的研究多注重于單一區(qū)域、大尺度，而對(duì)小尺度、多區(qū)域的研究較少［12-14］。小尺度研究區(qū)的自然影響變量較弱，其氣候、降水等諸多因子在區(qū)域內(nèi)變化較小。小區(qū)域植物群落物種多樣性的變化與分布最常見(jiàn)的形式就是植物群落隨著高程和地形的變化發(fā)生相應(yīng)的變化［26］。因此，在地理探測(cè)與地理加權(quán)聯(lián)合視角下，本研究圍繞中國(guó)典型氣候帶的7個(gè)代表地區(qū)，通過(guò)提取當(dāng)?shù)氐腘DVI與同研究視角下出現(xiàn)最多的3種地形因子（高程、坡度、坡向），探索不同地域的地形因子對(duì)植被分布的驅(qū)動(dòng)方式與特點(diǎn)，為研究區(qū)的災(zāi)害防治、植被保護(hù)和生態(tài)治理提供一定參考。

1 研究區(qū)概況

綜合考慮相對(duì)海拔超過(guò)800 m、地形多樣、人口密度較少、具有不同氣候特點(diǎn)等因素，最終選取了海南省白沙黎族自治縣、廣西壯族自治區(qū)蒲北縣、云南省魯?shù)榭h、青海省貴德縣、山西省方山縣、內(nèi)蒙古自治區(qū)巴林左旗、黑龍江省漠河市共7個(gè)縣域（100°58′—124°12′E，18°56′—53°33′N(xiāo)），從最南端的中熱帶到最北端的北溫帶，從沿海的南亞熱帶到內(nèi)陸的高原氣候帶，研究涉及的氣候帶豐富多樣，各研究區(qū)內(nèi)地形地貌相對(duì)復(fù)雜，詳情如圖1、表1所示。

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

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

多光譜數(shù)據(jù)主要來(lái)自于地理空間數(shù)據(jù)云（http：//www.gscloud.cn/）的美國(guó)Landsat 8國(guó)家遙感衛(wèi)星的OLI（陸地成像儀）2020年或2021年的成像數(shù)據(jù)，包括9個(gè)波段，主要選取7—9月，云量均小于6%的遙感衛(wèi)星數(shù)據(jù)，采用其近紅外波段（5波段）和紅光波段（4波段），空間分辨率為30 m。數(shù)字高程模型（DEM）來(lái)源于USUG網(wǎng)站的12.5 m空間分辨率的數(shù)字高程數(shù)據(jù)。中國(guó)行政邊界等矢量數(shù)據(jù)（審圖號(hào)：GS（2020）4619號(hào)）與中國(guó)氣候帶數(shù)據(jù)均來(lái)自于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（http：//www.resdc.cn）。

2.2 數(shù)據(jù)預(yù)處理

通過(guò)ENVI遙感數(shù)據(jù)處理軟件對(duì)符合季節(jié)與云量的Landsat 8數(shù)據(jù)進(jìn)行初步的輻射定標(biāo)、大氣校正、圖像增強(qiáng)、NDVI計(jì)算后，對(duì)獲得的NDVI數(shù)據(jù)進(jìn)行裁剪，然后利用ArcGIS軟件進(jìn)行NDVI最大值合成（MVC）。利用數(shù)字高程模型提取出地形因子數(shù)據(jù)（高程、坡度、坡向）后，按照自然間斷法將各研究區(qū)域的地形因子分為9類(lèi)。利用漁網(wǎng)建立各研究區(qū)100×100行的漁網(wǎng)點(diǎn)數(shù)據(jù)，使用行政邊界數(shù)據(jù)對(duì)其進(jìn)行篩選，以獲取各研究區(qū)的采樣點(diǎn)，數(shù)量為4000～4500個(gè)。然后，對(duì)采樣點(diǎn)的植被指數(shù)和重分類(lèi)后的地形因子數(shù)據(jù)進(jìn)行信息采集統(tǒng)計(jì)。

2.3 研究方法

2.3.1 歸一化植被指數(shù)計(jì)算

NDVI反映一定地域內(nèi)的植被生長(zhǎng)情況，對(duì)遙感影像的紅光（4波段）、近紅光（5波段）進(jìn)行波段運(yùn)算，利用綠色植物強(qiáng)吸收可見(jiàn)光紅光波段（0.6～0.7 μm）以及對(duì)植被具有高反射近紅外波段（0.7～1.1 μm）的特性，通過(guò)變換，增強(qiáng)植被信號(hào)，以削弱噪聲組合而成［27］。

2.3.2 地理探測(cè)器

地理探測(cè)器是基于驅(qū)動(dòng)因子、風(fēng)險(xiǎn)探測(cè)、生態(tài)與交互的一種探索空間異質(zhì)性與自然驅(qū)動(dòng)因子交互的空間統(tǒng)計(jì)模型［28］。

空間分異與地形因子探測(cè)：通過(guò)樣方點(diǎn)對(duì)重分類(lèi)的地形因子和植被指數(shù)進(jìn)行信息提取，解釋不同地形因子分級(jí)區(qū)域?qū)χ脖恢笖?shù)的影響程度。自然因子解釋力公式為

PD=1-∑Lh=1Nhσ2hNσ2=1-SSWSST，（2）

式中：PD為地形因子對(duì)植被指數(shù)（NDVI）的解釋力，值域范圍在［0，1］，值越趨近于1，其地形因子對(duì)植被生長(zhǎng)和空間分布影響程度越大；h=1，2，3，…L，是植被指數(shù)與地形因子因變量的分級(jí)數(shù)；L為各研究區(qū)內(nèi)地形因子的總分級(jí)數(shù)；Nh和N分別為層h和區(qū)域的單元數(shù)；σ2h和σ2分別是層h和區(qū)域內(nèi)樣方提取點(diǎn)的植被指數(shù)的方差；SSW和SST分別為層內(nèi)方差之和和區(qū)域總方差。

生態(tài)探測(cè)：比較2個(gè)地形因子對(duì)NDVI的空間分布狀況是否存在顯著的差異與區(qū)別，以F統(tǒng)計(jì)量來(lái)表示。

F=Nx1×（Nx2-1）×SSWx1Nx2×（Nx1-1）×SSWx2，（3）

式中：Nx1及Nx2分別表示2個(gè)地形因子的樣本數(shù)量；SSWx1和SSWx2分別表示由2個(gè)地形因子形成分層的層內(nèi)方差之和［29］。

地形因子交互探測(cè)：對(duì)地形因子PD（x1）、PD（x2）與PD（x1∩x2）進(jìn)行對(duì)比，評(píng)價(jià)2個(gè)地形因子對(duì)NDVI的交互影響作用（非線(xiàn)性減弱、單因子非線(xiàn)性減弱、互相增強(qiáng)、非線(xiàn)性增強(qiáng)、獨(dú)立增強(qiáng)）［12］。

2.3.3 地理加權(quán)回歸

地理加權(quán)回歸模型是一種非常有效的揭示空間非平穩(wěn)性和空間依賴(lài)的方法［30］。利用地理加權(quán)回歸模型探究地形因子同NDVI的變化趨勢(shì)與空間分布上存在的規(guī)律。公式為

yi=ρ0（u0+v0）+∑kj=1ρj（ui+vi）xij+εi，（4）

式中：yi為樣本i的NDVI變化擬合值，（ui+vi）為第i個(gè)樣點(diǎn)的空間地理坐標(biāo)，ρ0（u0+v0）是樣點(diǎn)的常數(shù)項(xiàng)估計(jì)數(shù)值，ρj（ui+vi）是樣點(diǎn)i的第j個(gè)回歸參數(shù)值，xij是第j個(gè)自變量在樣點(diǎn)i里面的數(shù)值，εi為服從均值為0的正態(tài)分布誤差［31］。

經(jīng)過(guò)驗(yàn)證，各區(qū)域的地形因子之間的VIF值均小于5。在參數(shù)選擇方面，采用核類(lèi)型為FIXED（“高斯”核）與交叉驗(yàn)證準(zhǔn)則的CV法，以確定核的范圍。

地理加權(quán)回歸分析后，以相等間隔法對(duì)局部擬合度（Local R2）進(jìn)行分級(jí)可視化，探索各研究區(qū)的地形因子（高程、坡度、坡向）與植被指數(shù)之間的空間分布規(guī)律和特點(diǎn)，即局部擬合優(yōu)度值越高，地形與植被的擬合程度越高，地形對(duì)植被的影響越大。

3 結(jié) 果

3.1 歸一化植被指數(shù)的空間分布特征

將NDVI分為低（－1，0.2］、中低（0.2，0.4］、中（0.4，0.6］、高（0.6，1］4個(gè)等級(jí)，各研究區(qū)植被等級(jí)狀況如圖2所示：漠河市、白沙黎族自治縣、浦北縣和魯?shù)榭h的植被狀況較好，其高植被覆蓋比例均超過(guò)69%；而方山縣、貴德縣和巴林左旗的植被狀況相對(duì)較差，其中巴林左旗的植被覆蓋最差，低和中低植被覆蓋率占比超過(guò)80%。總體而言，各研究區(qū)的植被指數(shù)呈現(xiàn)出由南向北、由沿海向內(nèi)陸遞減的趨勢(shì)。值得注意的是，雖然貴德縣的緯度較方山縣低，但植被覆蓋狀況卻較差。這是因?yàn)榍嗖馗咴瓪夂驇У奶厥獾乩砗铜h(huán)境條件所造成的，貴德縣海拔較高、地處內(nèi)陸且土壤貧瘠，限制了植被的生長(zhǎng)和發(fā)展。盡管漠河市位于我國(guó)最北端的內(nèi)陸地區(qū)，但與中/北熱帶的白沙黎族自治縣與南亞熱帶的浦北縣等南方沿海地區(qū)相比，其植被覆蓋狀況卻更好，這得益于當(dāng)?shù)厝丝谙∩?、土壤肥沃、水源豐富等多種原因共同促進(jìn)了當(dāng)?shù)刂脖坏纳L(zhǎng)。

各研究區(qū)的NDVI分布如圖3所示：白沙黎族自治縣、浦北縣、魯?shù)榭h、巴林左旗、漠河市植被覆蓋度隨著距離低植被覆蓋中心的距離增加而逐漸增加，形成了明顯的空間梯度。貴德縣和方山縣的植被覆蓋度呈現(xiàn)出線(xiàn)性分布的特征，以河流與低谷向四周逐漸擴(kuò)散并呈現(xiàn)逐漸變化的分層現(xiàn)象。

3.2 基于地理探測(cè)器的地形因子歸因分析

3.2.1 空間分異與地形因子探測(cè)分析

如圖4所示：各研究區(qū)的高程、坡度和坡向的PD均值分別為0.133、0.085和0.061，即地形因子對(duì)植被指數(shù)的解釋力排名為坡向因子lt;坡度因子lt;高程因子。以高程因子為第一解釋因子的地區(qū)有白沙黎族自治

縣、貴德縣、巴林左旗與漠河市，以坡度因子為第一解釋因子的地區(qū)僅有方山縣，而以坡向因子為第一解釋因子的地區(qū)有浦北縣與魯?shù)榭h。浦北縣的坡向因子和方山縣的坡度因子的解釋力都明顯高于高程因子對(duì)植被的解釋力。各地形探測(cè)因子之間的變化都較為顯著，位于高原氣候帶的貴德縣各地形因子對(duì)NDVI的影響最為明顯，地形因子PD值的和為0.475，遠(yuǎn)超第2位的方山縣（0.379）。中亞熱帶的魯?shù)榭h的各地形因子之間PD值變化相差不大，且對(duì)植被指數(shù)的影響都較低。

3.2.2 生態(tài)探測(cè)分析

結(jié)果顯示（表2）：研究區(qū)的地形探測(cè)因子之間大部分沒(méi)有顯著差異，僅位于高原氣候帶的貴德縣和南溫帶的方山縣的地形因子在生態(tài)探測(cè)中顯示出顯著差異。具體而言，貴德縣的高程與坡度、坡向?qū)χ脖恢笖?shù)NDVI的影響存在顯著差異，而方山縣的坡度與坡向?qū)χ脖恢笖?shù)NDVI的影響也存在顯著差異。

3.2.3 地形因子交互探測(cè)分析

如表3所示：地形因子對(duì)NDVI具有交互作用，呈現(xiàn)出非線(xiàn)性增強(qiáng)的互動(dòng)效應(yīng)。這表明不同研究區(qū)的地形因子對(duì)NDVI的影響并非簡(jiǎn)單獨(dú)立，而是通過(guò)協(xié)同促進(jìn)的過(guò)程來(lái)起作用。在各地形因子的交互作用中，貴德縣的因子交互作用最顯著，特別是高程∩坡向的交互作用（0.567）最為突出。另外，不同研究區(qū)的因子交互作用也存在顯著差異。貴德縣、巴林左旗、白沙黎族自治縣和漠河市的因子交互作用以高程∩坡向?yàn)樽罡?，而浦北縣和方山縣以坡度∩坡向的交互作用最高，魯?shù)榭h則以高程∩坡度的交互作用最高。

3.3 基于地理加權(quán)回歸的地形因子歸因分析

各研究區(qū)的地形因子與植被的局部擬合度排序?yàn)楦咴瓪夂驇зF德縣（0.387）gt;南溫帶方山縣（0.365）gt;中/北熱帶白沙黎族自治縣（0.307）gt;南亞熱帶浦北縣（0.270）gt;北溫帶漠河市（0.224）gt;中亞熱帶魯?shù)榭h（0.118）gt;中溫帶巴林左旗（0.100）。

如圖5所示：貴德縣與方山縣的局部擬合分布均沿河流向外擴(kuò)散，白沙黎族自治縣和浦北縣呈現(xiàn)散點(diǎn)狀聚集，局部擬合分布格局與NDVI的空間分布格局相似。結(jié)合圖1中DEM數(shù)據(jù)可以看出，各研究區(qū)的地形變化與局部擬合度空間分布存在一定的關(guān)聯(lián)，如貴德縣與方山縣的沿河流分布的河谷與河流階地，白沙黎族自治縣和浦北縣散點(diǎn)狀分布的平地與谷底與圖5的空間分布格局對(duì)應(yīng)。

4 討 論

從7個(gè)研究區(qū)的植被指數(shù)狀況可以看出，不同氣候帶的植被覆蓋狀況整體由南向北，由沿海向內(nèi)陸呈下降趨勢(shì)的空間格局，與羅爽等［32］得到的中國(guó)植被空間分布格局結(jié)果相似。這與中國(guó)特定的地理環(huán)境與空間位置有一定聯(lián)系。中國(guó)擁有遼闊的緯度范圍（從南緯3°到北緯54°）與獨(dú)特的海陸位置（從綿延的海岸線(xiàn)再到西北部的崇山峻嶺），這導(dǎo)致了中國(guó)廣泛的氣候變化。南部地區(qū)位于亞熱帶和熱帶，氣候溫暖潮濕；中部地區(qū)屬于溫帶，四季分明；北部地區(qū)則屬于寒帶，冬季漫長(zhǎng)嚴(yán)寒。同時(shí)，高山脈與高原的地貌特征對(duì)氣候形成

起到重要的作用，高山脈和高原會(huì)影響氣候系統(tǒng)的穩(wěn)定性，產(chǎn)生垂直氣候帶和地形風(fēng)，導(dǎo)致氣溫、降水等氣候要素在不同高度和區(qū)域之間發(fā)生明顯變化，從而形成了中國(guó)特有植被空間變化的趨勢(shì)。

從各地形因子的均值看，高程因子對(duì)植被指數(shù)的解釋力最大。多數(shù)氣候帶高程因子（貴德縣、漠河市、白沙黎族自治縣）對(duì)植被指數(shù)變化有顯著作用，與彭文甫等［12］、王睿卿等［33］的研究較為吻合。山地的氣候條件及土壤狀況等自然因素的變化由高程決定，隨著海拔的升高，溫度與降水等出現(xiàn)相應(yīng)變化，從而形成了明顯的植被垂直帶［18］，這是多數(shù)研究中高程因子對(duì)植被空間的變化影響作用更顯著的可能原因。山西省方山縣則是以坡向因子發(fā)揮更大的作用，與鄧椿等［34］、邵嘉豪等［24］對(duì)地形因子探測(cè)結(jié)果一致。山西省位于南溫帶的黃土丘陵區(qū)，溝壑縱橫，氣溫偏低，降水偏少，自然植被包括草地與針葉林。結(jié)合遙感影像、地形數(shù)據(jù)與局部擬合空間分布數(shù)據(jù)，發(fā)現(xiàn)地勢(shì)平緩的地區(qū)種植有大量的作物。夏天，方山縣的農(nóng)作物生長(zhǎng)狀況常優(yōu)于該地區(qū)的自然植被，干擾了垂直地帶性的規(guī)律，使得該地區(qū)的坡度因子對(duì)植被的影響作用更顯著。巴林左旗與魯?shù)榭h的地形因子對(duì)植被因子的解釋力顯著低于其他研究區(qū)，揭示了巴林左旗與魯?shù)榭h的地形因子在NDVI的空間分布與形成過(guò)程中的弱影響機(jī)制。

根據(jù)交互探測(cè)發(fā)現(xiàn)，不同氣候帶的地形因子之間都不是簡(jiǎn)單獨(dú)立的對(duì)植被指數(shù)NDVI起作用，而是非線(xiàn)性增強(qiáng)促進(jìn)的過(guò)程，與許多專(zhuān)家學(xué)者得出的地形因子與植被間的交互結(jié)果相似［34-36］。

目前地形對(duì)植被的解釋力與空間分異的研究正趨近成熟，但植被的空間格局與生長(zhǎng)狀況不僅受地形因子的驅(qū)動(dòng)，同樣受溫度、土壤、水、人類(lèi)活動(dòng)、地質(zhì)活動(dòng)［37］等因素的影響。本研究通過(guò)約束空間范圍和增加典型樣區(qū)以增強(qiáng)地形因子的影響力，能較好地反應(yīng)地形因子對(duì)植被的影響。但仍然忽略了非自然等因素的作用，后續(xù)將繼續(xù)加強(qiáng)對(duì)小尺度、多區(qū)域的非自然因素與地形因子的綜合研究，為植被的區(qū)域生態(tài)價(jià)值與災(zāi)害防護(hù)提供更科學(xué)的指導(dǎo)與參考。

參考文獻(xiàn)：

［1］ FORKEL M，CARVALHAIS N，RDENBECK C，et al.Enhanced seasonal CO2 exchange caused by amplified plant productivity in northern ecosystems［J］.Science，2016，351（6274）：696-699.

［2］ PENG J，LIU Z，LIU Y，et al.Trend analysis of vegetation dynamics in Qinghai-Tibet Plateau using Hurst Exponent［J］.Ecological Indicators，2012，14（1）：28-39.

［3］ PICKUP G，BASTIN G N，CHEWINGS V H.Identifying trends in land degradation in non-equilibrium rangelands［J］.Journal of Applied Ecology，1998，35（3）：365-377.

［4］ 駱社周，程峰，王方建，等.基于TM遙感數(shù)據(jù)的西藏林芝地區(qū)葉面積指數(shù)反演［J］.遙感技術(shù)與應(yīng)用，2012，27（5）：740-745.

［5］ 朱緒超，袁國(guó)富，易小波，等.基于Landsat 8 OLI影像的塔里木河下游河岸林葉面積指數(shù)反演［J］.干旱區(qū)地理，2014，37（6）：1248-1256.

［6］ ZHENG Y D，XU Y C，YAN H J，et al.Analysis of Influencing factors in wheat/maize canopy RVI and NDVI acquisition using ground-based remote sensing system［J］.Spectroscopy and Spectral Analysis，2021，41（8）：2578-2585.

［7］ 安德帥，徐丹丹，劉月，等.高光譜與擬合多光譜植被指數(shù)反演武夷山亞高山草甸L(zhǎng)AI的對(duì)比研究［J］.生態(tài)科學(xué)，2022，41（5）：187-196.

［8］ ZHAO L，DAI A，DONG B.Changes in global vegetation activity and its driving factors during 1982-2013［J］.Agricultural and Forest Meteorology，2018，249：198-209.

［9］ FU G，SUN W，LI S W，et al.Modeling aboveground biomass using MODIS images and climatic data in grasslands on the Tibetan Plateau［J］.Journal of Resources and Ecology，2017，8（1）：42-49.

［10］LUAN J K，LIU D F，HUANG Q，et al.Analysis of the spatial-temporal change and impact factors of the vegetation index in Yulin，Shaanxi Province，in the last 17 years［J］.Acta Ecologica Sinica，2018，38（8）：2780-2790.

［11］WANG T，BAI H Y.Variation of vegetation NDVI in response to climate changes and human activities in Qinling Mountains［J］.Mountain Research，2017，35（6）：778-789.

［12］彭文甫，張冬梅，羅艷玫，等.自然因子對(duì)四川植被NDVI變化的地理探測(cè)［J］.地理學(xué)報(bào)，2019，74（9）：1758-1776.

［13］呂泳潔，丁文廣，鄧喆，等.基于地理探測(cè)器的甘肅省植被覆蓋時(shí)空變化及解釋力分析［J］.地球環(huán)境學(xué)報(bào)，2022，13（2）：185-195.

［14］廖成浩，曾艷，姚昆，等.江西植被EVI時(shí)空動(dòng)態(tài)監(jiān)測(cè)及驅(qū)動(dòng)因素分析［J］.云南大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，44（5）：981-989.

［15］羅明良，湯國(guó)安，劉學(xué)軍.地形形態(tài)要素研究綜述［C/OL］//中國(guó)地理學(xué)會(huì)2007年學(xué)術(shù)年會(huì)論文集.https：//xueshu.baidu.com/usercenter/paper/show？paperid=8d62cf840d8da22bf3822d6e47dbb953amp;site=xueshu_seamp;hitarticle=1.

［16］傅抱璞.地形和海拔高度對(duì)降水的影響［J］.地理學(xué)報(bào)，1992，47（4）：302-314.

［17］FLORINSKY I V ，KURYAKOVA G A .Influence of topography on some vegetation cover properties［J］.Catena，1996，27（2）：123-141.

［18］李娟，龔純偉.祁連山國(guó)家公園植被覆蓋變化地形分異效應(yīng)［J］.水土保持通報(bào)，2021，41（3）：228-237.

［19］劉梁美子，占車(chē)生，胡實(shí)，等.黔桂喀斯特山區(qū)年NDVI變化的影響因素研究［J］.地理科學(xué)進(jìn)展，2019，38（11）：1783-1792.

［20］陳佳敏，王寧，趙富王，等.黃土丘陵區(qū)溝沿線(xiàn)不同微地形植物群落特征［J］.中國(guó)水土保持科學(xué)，2022，20（5）：31-38.

［21］龍映豪，李旭東，程?hào)|亞.山區(qū)植被覆蓋度變化的地形分異特征：以貴州開(kāi)陽(yáng)縣為例［J］.水土保持研究，2022，29（5）：268-275，282.

［22］銀朵朵，王艷慧.溫帶大陸性半干旱季風(fēng)氣候帶植被覆蓋度時(shí)空變化及其地形分異研究［J］.生態(tài)學(xué)報(bào)，2021，41（3）：1158-1167.

［23］李夢(mèng)華，韓穎娟，趙慧，等.基于地理探測(cè)器的寧夏植被覆蓋度時(shí)空變化特征及其驅(qū)動(dòng)因子分析［J］.生態(tài)環(huán)境學(xué)報(bào)，2022，31（7）：1317-1325.

［24］邵嘉豪，李晶，閆星光，等.基于地理探測(cè)器的山西省2000～2020年NPP時(shí)空變化特征及解釋力分析［J］.環(huán)境科學(xué)，2023，44（1）：312-322.

［25］李曉靜，張斌，劉守江，等.涼山州喜德縣植被蓋度與地形因子的空間耦合［J］.西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，43（2）：129-136.

［26］馬梓策，于紅博，曹聰明，等.中國(guó)植被覆蓋度時(shí)空特征及其影響因素分析［J］.長(zhǎng)江流域資源與環(huán)境，2020，29（6）：1310-1321.

［27］張薇，萬(wàn)虹麟，陳晨，等.基于歸一化植被指數(shù)的植被覆蓋度時(shí)空變化研究：以大運(yùn)河文化帶滄州段為例［J］.河北水利電力學(xué)院學(xué)報(bào)，2021，31（4）：72-78.

［28］王勁峰，徐成東.地理探測(cè)器：原理與展望［J］.地理學(xué)報(bào)，2017，72（1）：116-134.

［29］HEIN L，DE R N，HIERNAUX P，et al.Desertification in the Sahel：Towards better accounting for ecosystem dynamics in the interpretation of remote sensing images［J］.Journal of Arid Environments，2011，75（11）：1164-1172.

［30］覃文忠，王建梅，劉妙龍.地理加權(quán)回歸分析空間數(shù)據(jù)的空間非平穩(wěn)性［J］.遼寧師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，28（4）：476-479.

［31］李晶晶，閆慶武，胡苗苗.基于地理加權(quán)回歸模型的能源“金三角”地區(qū)植被時(shí)空演變及主導(dǎo)因素分析［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2018，34（8）：700-708.

［32］羅爽，劉會(huì)玉，龔海波.1982—2018年中國(guó)植被覆蓋變化非線(xiàn)性趨勢(shì)及其格局分析［J］.生態(tài)學(xué)報(bào)，2022，42（20）：8331-8342.

［33］王睿卿，蔣曉輝，聶桐.基于地理探測(cè)器的珠江流域NDVI時(shí)空變化及解釋力分析［J］.人民珠江，2022，43（7）：61-73.

［34］鄧椿，蔣曉輝，聶桐，等.山西省植被覆蓋度多因子探測(cè)特征與解釋力分析［J］.環(huán)境科學(xué)與技術(shù)，2022，45（2）：182-191.

［35］張思源，聶瑩，張海燕，等.基于地理探測(cè)器的內(nèi)蒙古植被NDVI時(shí)空變化與解釋力分析［J］.草地學(xué)報(bào)，2020，28（5）：1460-1472.

［36］姜萍，胡列群，肖靜，等.新疆植被NDVI時(shí)空變化及定量歸因［J］.水土保持研究，2022，29（2）：212-220.

［37］劉守江，朱兆棋，陳曦，等.龍門(mén)山銀廠(chǎng)溝區(qū)域5·12地震前后災(zāi)害與生態(tài)研究進(jìn)展與展望［J］.西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，42（3）：239-245.

Quantitative Attribution of Topographic Factors InfluencingVegetation Indices in Typical Climatic Zones：

Based on Geographic Detector and GWR

YIN Jia-yuan，LUO Ming-liang，BAI Lei-chao，WU Qiu-sheng

（a.School of Geographical Sciences，b.Jialing River Basin Research Institute，China West Normal University，

Nanchong Sichuan 637009，China）

Abstract：The research on the changes of normalised different vegetation index （NDVI） influenced by topographic factors in different climatic zones is of great value for the exploration of soil and water hazards，ecological research，and regional environmental quality.Taking the county areas in seven typical climate zones as examples and the Landsat 8 remote sensing data and DEM data as the basis，this study explores the relationship between topographic factors and changes of NDVI as well as its variation patterns at the small area scale of the typical climate zones in China through geographic detector and geographic weighted regression.The results are as follows：（1）The NDVI is ranked as：Mohe city （the northern temperate zone）gt;Baisha Lizu autonomous county （the central/northern tropical zone）gt;Pubei county （the southern subtropical zone）gt;Ludian county （the southern subtropical zone）gt;Fangshan county （the southern temperate zone）gt;Guide county （the plateau climatic zone）gt;Bairin left banner （the middle temperate zone）.（2）In terms of their explanatory power for the vegetation index，the topographic factors are ranked as：slope direction factorlt;slope factorlt;elevation factor，and most of the topographic factors show no significant difference；in different study areas，topographic factors exert its influence on NDVI through nonlinear enhancement rather than in a simple and independent way.（3）The local fitted spatial pattern of topographic factors is similar to the spatial distribution pattern of vegetation cover，so topography has an important influence on the spatial distribution of vegetation.（4） According to the adjustment of the goodness of fit，the ranking of the fitting of topographic factors on vegetation in each study area is：Guide countygt;Fangshan countygt;Baisha Lizu autonomous countygt;Pubei countygt;Mohe citygt;Ludian countygt;Bairin left banner，and the fitting degree of topographic factors on vegetation index varies significantly with the geographic area.

Keywords：geographic detector；geographic weighted regression；topographic factors；NDVI；climatic zones