摘 要：【目的】為較全面掌握川西高原植被NDVI動態(tài)變化及驅(qū)動因素?！痉椒ā恳訫ODIS-NDVI為基礎(chǔ)數(shù)據(jù)，借助Sen+MK和地理探測器模型對其2000—2020年變化規(guī)律展開探討?！窘Y(jié)果】①川西高原植被NDVI整體呈現(xiàn)自西南向東北方向傾斜，植被NDVI數(shù)值整體相對更高；②近20年，全域植被NDVI數(shù)值無顯著變化趨勢的地區(qū)全域面積占比約75%；③高程、年平均氣溫和土壤類型主導(dǎo)了川西高原植被NDVI的空間布局，解釋力在35%以上?！窘Y(jié)論】川西高原近20年植被整體生長狀態(tài)良好，自然因素對其植被空間布局的影響顯著。

關(guān)鍵詞：植被NDVI；地理探測器；川西高原

中圖分類號：Q948" " "文獻(xiàn)標(biāo)志碼：A" " 文章編號：1003-5168（2023）13-0094-05

DOI：10.19968/j.cnki.hnkj.1003-5168.2023.13.019

Study on NDVI Changes of Vegetation in Western Sichuan Plateau

CHENG Bo ZHOU Xianyang YAO Kun DIDI Keniu YU Lei FAN Xiao

（School of Resources and Environment， Xichang University， Xichang 615000，China）

Abstract： [Purposes] This paper aims to comprehensively understand the NDVI dynamic changes and driving factors of vegetation in the western Sichuan plateau. [Methods] This study conducted a discussion on its change rules from 2000 to 2020 based on MODIS-NDVI data， using Sen+MK and geographic detector models. [Findings] The results are as follows： ①The NDVI values of vegetation on the western Sichuan plateau as a whole show a slope from southwest to northeast， and the NDVI values of vegetation are relatively higher as a whole; ②In the past 20 years， the area with no significant change trend in NDVI values of vegetation throughout the region accounts for about 75%；③Elevation， annual average temperature， and soil type are the dominant factors affecting the NDVI spatial layout of vegetation in the study area， with an explanatory power of over 35%. [Conclusions] The overall growth of vegetation on the western Sichuan plateau in the past 20 years is good， and natural factors have a significant impact on its spatial layout of vegetation.

Keywords： vegetation NDVI; geographic detectors; western sichuan plateau

0 引言

植被是陸地生態(tài)系統(tǒng)的重要組成部分之一，是實現(xiàn)土壤、大氣和能量進(jìn)行物質(zhì)交換的重要“紐帶”［1-2］。受全球氣候變暖、降雨豐沛季節(jié)性分布不均、不合理資源開采等自然和人類活動的干擾，區(qū)域植被生長狀態(tài)和覆蓋程度均發(fā)生著顯著變化［3］。因此，及時、全面了解區(qū)域植被NDVI（Normalized Difference Vegetation Index，歸一化植被指數(shù)）的變化規(guī)律，并對影響其改變的因素展開探討，有助于區(qū)域植被的保護(hù)與恢復(fù)［4-5］。

近年來，隨著遙感技術(shù)的不斷發(fā)展和完善，利用遙感獲取植被指數(shù)完成其變化的監(jiān)測分析已成為當(dāng)前研究熱點之一［6］。植被NDVI作為描述植被生長狀態(tài)與覆蓋度的綜合指標(biāo)之一，在相關(guān)領(lǐng)域被廣泛應(yīng)用［7］。王麗霞等［8］基于分辨率為500 m的MODIS數(shù)據(jù)完成渭河流域2000—2015年植被NDVI時空變化規(guī)律探究，分析了氣候因子對其變化的影響；鄭勇等［9］基于NDVI數(shù)據(jù)實現(xiàn)川西高原近20年植被覆蓋變化與驅(qū)動分析，發(fā)現(xiàn)18年內(nèi)區(qū)域植被覆蓋整體呈現(xiàn)好轉(zhuǎn)的發(fā)展趨勢；張華等［10］基于地理探測器完成祁連山地區(qū)2000—2019年植被NDVI的變化規(guī)律與驅(qū)動，發(fā)現(xiàn)其近20年內(nèi)植被處于逐漸好轉(zhuǎn)的發(fā)展趨勢，自然因子是影響NDVI變化的主要因素；王葉蘭等［11］對川西高原2001—2021年植被NDVI的時空變化規(guī)律，發(fā)現(xiàn)全域大部分地區(qū)植被NDVI數(shù)值呈升高的變化趨勢，高程對其影響作用最顯著。以上成果均對區(qū)域植被NDVI的變化具有重要價值。

川西高原是四川生態(tài)西部高原的重要建設(shè)區(qū)，加強對該地區(qū)植被覆蓋及生長狀態(tài)的監(jiān)測十分必要。目前針對川西高原開展的研究相對較少，因此本研究基于MODIS-NDVI，結(jié)合Sen氏斜率、MK檢驗及地理探測器，探討該地區(qū)近20年植被NDVI變化規(guī)律及驅(qū)動因素。

1 研究區(qū)概況

川西高原地處青藏高原東南緣，主要包括甘孜和阿壩州兩個地級市，面積約2.3×106" km2。地勢呈現(xiàn)自西北向東南傾斜的狀態(tài)，西北相比東南地區(qū)海拔整體相對較高，為典型的高原氣候，冬季寒冷漫長，雨旱分明。

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

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

研究所用的植被NDVI為MOD13Q1產(chǎn)品成果，重訪周期為16天，分辨率為250 m。對基礎(chǔ)數(shù)據(jù)進(jìn)行投影變換、鑲嵌、裁剪等預(yù)處理，為消除云、霧和氣溶膠等噪聲對數(shù)據(jù)精度的干擾，采用最大合成對數(shù)據(jù)進(jìn)行進(jìn)一步處理，最終采用的各年植被NDVI成果為年最大合成結(jié)果。

2.2 研究方法

2.2.1 Sen氏斜率。為解釋2000—2020年川西北高原植被NDVI隨時間的變化規(guī)律，以Sen氏斜率對其進(jìn)行計算[12]，表達(dá)式如式（1）。

2.2.2 MK檢驗。相比傳統(tǒng)的F或者T數(shù)據(jù)顯著性檢驗方法，MK檢驗不需要數(shù)據(jù)滿足特定的分布形態(tài)，能夠有效降低像元數(shù)值缺失或異常對計算結(jié)果準(zhǔn)確性的干擾［13］。

研究結(jié)合參考前人成果［13-14］，將植被NDVI斜率的顯著性檢驗結(jié)果分為：極顯著退化（Slopelt;0，P≤0.01）、顯著退化（Slopelt;0，0.01≤P≤0.05）、無顯著退化（Slopelt;0，P≥0.05）、無顯著改善（Slope＞0，P≥0.05），顯著改善（Slope＞0，0.01≤P≤0.05）和極顯著改善（Slope＞0，P＜0.01）。

2.2.3 地理探測器。地理探測器作為分析變量空間分異并對其背后驅(qū)動力進(jìn)行解釋的工具，其不需要自變量和因變量之間滿足某種特定假設(shè)關(guān)系，已被廣泛應(yīng)用于地理現(xiàn)象影響因素的探索分析［15-17］。因子探測器用于分析各因素對植被NDVI空間布局的影響作用，模型如式（2）和式（3）。

3 結(jié)果分析

3.1 分布格局

川西高原植被NDVI平均值分布如圖1所示，分析可知：該地區(qū)植被覆蓋在空間上整體呈現(xiàn)出明顯的分布差異性，呈現(xiàn)出自西南向東北傾斜方向，植被整體生長狀態(tài)逐漸變化好，覆蓋狀況整體有相對較好地發(fā)展趨勢。

參考前人成果［12］及NDVI分布規(guī)律，以等間距分類將該地區(qū)植被NDVI分為5個等級。結(jié)合圖1可知，低和中低植被NDVI覆蓋區(qū)全域面積占比僅約3.22%，結(jié)合地形、土地和植被等數(shù)據(jù)，它們主要分布于海拔幾乎接近5 000 m以上的高山地區(qū)，這些地域地勢險峻，主要被裸巖石質(zhì)地、沙地或裸土地覆蓋，除部分區(qū)域分布有極少許高山稀疏植被外，大多地區(qū)幾乎為無植被生長區(qū)；中植被NDVI覆蓋區(qū)全域面積占比約7.70%，它們主要集中分布于石渠縣北部、新龍縣西部和東南部、稻城縣東南與理塘縣界交界處，此外在九龍縣、康定市和道孚縣等部分地區(qū)也有零散分布，結(jié)合相關(guān)資料分析可知，這些區(qū)域海拔高度范圍在4 500～5 000 m之間，主要分布的植被種類包括高寒草甸、灌叢和高山稀疏植被等，土地景觀類型以低覆蓋度草地和灌木林地為主，同時少部分區(qū)域還有中覆蓋度草地；相比下，中高和高植被NDVI覆蓋區(qū)在全域的分布范圍最廣全域面積占比達(dá)89.08%，在全域海拔小于4 500 m的地區(qū)均有分布，受區(qū)域海拔整體呈現(xiàn)出西北向東南逐漸降低的趨勢變化的影響，這些范圍內(nèi)東部植被覆蓋狀況整體較好。

3.2 趨勢發(fā)展

結(jié)合Sen計算結(jié)果，川西高原植被NDVI斜率及顯著性檢驗如圖2所示。分析可知，2000—2020年，川西高原植被NDVI斜率計算數(shù)值為負(fù)數(shù)的像元數(shù)量全域占比約29.65%，它們大部分零散分布于高原西北部、北部、南部和東南部地區(qū)；約70.35%的地區(qū)植被NDVI斜率數(shù)值大于零，客觀程度表明在整個研究期內(nèi)川西高原植生長整體呈現(xiàn)良好的發(fā)展態(tài)勢。

結(jié)合顯著性檢驗對該地區(qū)植被NDVI斜率變化結(jié)果展開進(jìn)一步分析可知：整個研究時段內(nèi)，全域植被NDVI呈現(xiàn)數(shù)值極顯著（退化或改善）改變的地區(qū)全域面積占比約6.62%，而在極顯著改變的地區(qū)中有0.64%的區(qū)域植被NDVI呈現(xiàn)退化的發(fā)展?fàn)顟B(tài)，植被NDVI極顯著退化區(qū)主要分布于石渠和色達(dá)縣中部，而極顯著改善區(qū)則主要分布于高原東北部的若爾蓋、阿壩、紅原、九寨溝及松潘等縣；植被NDVI顯著改善或退化地區(qū)全域面積占比約11.19%，它們主要環(huán)繞極顯著改變區(qū)形成緩沖帶，分布范圍不斷向外進(jìn)行擴(kuò)張；相比之下，21年內(nèi)植被NDVI斜率呈現(xiàn)無顯著退化或改善地區(qū)在全域分布范圍最廣，全域面積占比達(dá)到75.47%，面積分布已在全域占據(jù)絕對優(yōu)勢；整個2000—2020年，川西高原植被NDVI數(shù)值呈現(xiàn)極顯著、顯著改善及無顯著改善區(qū)域地區(qū)已占據(jù)全域總面積的82.72%。

3.3 影響因素分析

為了解各影響因素對川西高原植被NDVI空間布局的影響，利用因子探測器對其進(jìn)行計算，結(jié)果見表1。

分析可知，各因子對川西北高原植被NDVI空間布局的影響作用也呈現(xiàn)顯著差異，解釋力從大到小分別為高程＞平均氣溫＞土壤類型＞植被類型＞土地利用類型＞降雨量＞坡度＞GDP＞坡向＞人口密度。高程對研究區(qū)植被NDVI空間布局的解釋力最強，q的數(shù)值達(dá)到0.472 2位列各影響因素之首，該地區(qū)為青藏高原東南緣，是“世界屋脊”的重要組成部分，海拔數(shù)值的差異對氣溫、日照強度及時長等自然環(huán)境均有顯著影響，因此它在本研究選取的指標(biāo)中對植被NDVI布局的影響作用最明顯。平均氣溫和土壤類型對植被NDVI布局的影響作用基本相似，q值分別為0.356 6和0.341 1，解釋力均在35%左右。適宜的溫度對促進(jìn)植被的生長起到積極促進(jìn)作用，相反溫度過高或過低反而會對植被生長狀態(tài)產(chǎn)生抑制作用，該地區(qū)平均氣溫呈現(xiàn)東部和南部相對較高的特點，川西高原北部和西部地區(qū)平均氣溫相對較低，年均氣溫主要維持在10 ℃以下不利于植被生長，而東部和南部年平均氣溫卻在10 ℃～25 ℃之間是較適宜植被生長的區(qū)間，因此東部和南部地區(qū)植被生長狀態(tài)相比下更好，這也促成了平均氣溫會對植被NDVI的布局產(chǎn)生影響。土壤作為植被生長的基礎(chǔ)，作為物質(zhì)能量交換的重要基質(zhì)場所，整個高原土壤類型多達(dá)24種且呈現(xiàn)顯著的垂直地帶性，隨著海拔的增高土壤的營養(yǎng)成分和類型也發(fā)生急劇改變，從適宜其生長的類型向抑制其生長的種類轉(zhuǎn)變，因此土壤類型對植被NDVI空間布局影響作用較顯著。植被類型和土地利用類型對植被NDVI布局的影響作用也基本相似，q值分別為0.228 8和0.216 5，解釋力均在20%左右。不同類型的植被對生境的要求有著顯著的差異，生長狀態(tài)自然也會存在顯著區(qū)別，而土地類型是人類活動作為人類活動的典型代表，全域超過70%的地區(qū)分布著林草地，林草地分布范圍若發(fā)生改變則植被生長則會變化。除去以上提交的指標(biāo)外，剩余指標(biāo)中降雨量和坡度的解釋力分別為5.89%和1.45%影響作用均非常有限，而GDP、坡向和人口密度的影響作用則幾乎可忽略不計。

4 結(jié)論

川西高原是長江上游重要的生態(tài)屏障，本研究以MODIS數(shù)據(jù)為基礎(chǔ)，對其2000—2020年植被NDVI的空間分布、變化規(guī)律及布局影響因素展開分析，主要結(jié)果如下。

①受先天自然環(huán)境的影響，川西高原植被NDVI空間布局變化規(guī)律顯著，全域整體呈現(xiàn)東部相比西部、北部相比南部植被NDVI數(shù)值較高的特征。

②21年內(nèi)，植被NDVI無顯著改變、顯著改善和極顯著改善區(qū)在全域面積占比已超過85%，川西高原植被NDVI整體呈現(xiàn)良好的發(fā)展態(tài)勢。

③各驅(qū)動因素對川西高原植被NDVI空間布局影響作用差異顯著。高程對植被NDVI空間布局的解釋力均在45%以上，平均氣溫和土壤類型解釋力在35%左右，植被類型和土地利用類型的解釋力在20%～30%之間，剩余其余指標(biāo)的解釋力均在10%以下。

參考文獻(xiàn)：

［1］石淞，李文，翟育涔，等.中國東北虎豹國家公園植被NDVI時空變化及原因探究［J/OL］.南京林業(yè)大學(xué)學(xué)報（自然科學(xué)版）：1-12［2023-03-22］.http：//kns.cnki.net/kcms/detail/32.1161.S.20230320.1610.004.html.

［2］王葉蘭，楊鑫，郝利娜.2001—2021年川西高原植被NDVI時空變化及影響因素分析［J/OL］.自然資源遙感：1-9［2023-03-22］.http：//kns.cnki.net/kcms/detail/10.1759.p.20230207.1452.010.html.

［3］繆利，陸晴，劉根林，等.1999—2019年青藏高原不同植被類型NDVI時空演變特征及其對氣候因子的響應(yīng)［J］.水土保持研究，2023，30（1）：97-105.

［4］胡志強，宋孝玉，覃琳，等.2001—2020年內(nèi)陸干旱草原縣域尺度植被NDVI時空變化特征及其驅(qū)動因素［J］.水土保持通報，2022，42（5）：213-221.

［5］歐陽習(xí)軍，董曉華，魏榕，等.青藏高原植被生長季NDVI時空變化及對氣候因子的響應(yīng)分析［J］.水土保持研究，2023，30（2）：220-229.

［6］徐勇，戴強玉，黃雯婷，等.2000—2020年西南地區(qū)植被NDVI時空變化及驅(qū)動機(jī)制探究［J］.環(huán)境科學(xué)，2023，44（1）：323-335.

［7］徐勇，鄭志威，郭振東，等.2000—2020年長江流域植被NDVI動態(tài)變化及影響因素探測［J］.環(huán)境科學(xué)，2022，43（7）：3730-3740.

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

［9］鄭勇，楊武年，劉沖，等.川西高原近20a植被覆蓋變化遙感動態(tài)監(jiān)測及驅(qū)動力分析［J］.遙感技術(shù)與應(yīng)用，2020，35（6）：1447-1456.

［10］張華，李明，宋金岳，等.基于地理探測器的祁連山國家公園植被NDVI變化驅(qū)動因素分析［J］.生態(tài)學(xué)志，2021，40（8）：2530-2540.

［11］王葉蘭，楊鑫，郝利娜.川西高原植被物候及其對氣候變化的響應(yīng)［J/OL］.長江科學(xué)院院報：1-10［2023-03-22］.http：//kns.cnki.net/kcms/detail/42.1171.tv.20220715.1632.002.html.

［12］陶帥，鄺婷婷，彭文甫，等.2000—2015年長江上游NDVI時空變化及驅(qū)動力：以宜賓市為例［J］.生態(tài)學(xué)報，2020，40（14）：5029-5043.

［13］向珈瑤，彭文甫，陶帥.2000—2020年岷江上游植被NDVI時空變化及其地形響應(yīng)［J］.長江流域資源與環(huán)境，2022，31（7）：1534-1547.

［14］涂又，姜亮亮，劉睿，等.1982—2015年中國植被NDVI時空變化特征及其驅(qū)動分析［J］.農(nóng)業(yè)工程學(xué)報，2021，37（22）：75-84.

［15］丁永康，葉婷，陳康.基于地理探測器的滹沱河流域植被覆蓋時空變化與驅(qū)動力分析［J］.中國生態(tài)農(nóng)業(yè)學(xué)報（中英文），2022，30（11）：1737-1749.

［16］高思琦，董國濤，蔣曉輝，等.基于地理探測器的三江源植被變化及自然驅(qū)動因子分析［J］.水土保持研究，2022，29（4）：336-343.

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