摘要 為明確延邊朝鮮族自治州（以下簡稱延邊州）植被指數(shù)變化規(guī)律及其對氣候因子的響應(yīng)機制，利用2000—2022年的MOD13Q1數(shù)據(jù)，結(jié)合降水量、氣溫等資料，通過一元線性回歸法、Theil-Sen Median 趨勢分析、Mann-Kendall檢驗法和相關(guān)系數(shù)分析法等，從時間和空間2個尺度分析延邊州NDVI變化特征及其與氣溫和降水量的響應(yīng)程度。結(jié)果表明：時間變化上來看，2000—2022年延邊州植被NDVI多年平均值為0.539，呈波動上升趨勢，上升斜率為0.003 6/年。各縣（市）植被NDVI存在差異，但總體均呈現(xiàn)相似的波動上升趨勢。其中，和龍市年均NDVI最高，為0.560，延吉市最低，為0.515。圖們市NDVI年均值上升趨勢最為明顯，為0.004 3/年，和龍市NDVI年均值上升趨勢最慢，為0.003 1/年。空間分布上來看，延邊州植被整體覆蓋狀況較好，植被NDVI年均值大于0.6的區(qū)域面積占總面積的22.93%，NDVI年均值在0.4～0.6的占67.16%；植被NDVI年均值小于0.4的面積僅占9.91%，主要分布在各縣（市）中心及周邊區(qū)域。植被NDVI呈上升趨勢的面積占研究區(qū)總面積的94.23%，其中，呈極顯著上升趨勢的面積占71.80%，呈下降趨勢的面積僅占研究區(qū)面積的2.48%，零星分布在部分耕地區(qū)域。植被NDVI對氣候因子響應(yīng)方面，延邊州大部分地區(qū)NDVI與氣溫、降水量呈正相關(guān)，約有67.23%的地區(qū)與降水量呈顯著正相關(guān)，大于與氣溫呈顯著正相關(guān)區(qū)域（59.49%），表明全州范圍內(nèi)NDVI與降水量和氣溫均以正向促進為主，且與降水量的相關(guān)性強于氣溫。

關(guān)鍵詞 NDVI；時空變化特征；氣候因子；響應(yīng)程度；延邊州

中圖分類號 Q948.112 文獻標識碼 A 文章編號 0517-6611（2024）21-0187-05

doi：10.3969/j.issn.0517-6611.2024.21.039

開放科學(xué)（資源服務(wù)）標識碼（OSID）：

The Spatiotemporal Variation Characteristics of NDVI in Yanbian Korean Autonomous Prefecture and Its Relationship with Climate Factors

WANG Jing-jie，ZHAO Jing，LIU Yan-bin et al

（Meteorological Bureau of Yanbian Korean Autonomous Prefecture，Yanji，Jilin 133001）

Abstract In order to clarify the variation law of vegetation index in Yanbian Korean Autonomous Prefecture and its response mechanism to climate factors，MOD13Q1 data from 2000 to 2022 were used，combined with precipitation，temperature and other data.Through univariate linear regression，Theil Sen Median trend analysis，Mann-Kendall test and correlation coefficient analysis，the NDVI variation characteristics of Yanbian Prefecture and their response degree to temperature and precipitation were analyzed from two scales of time and space.The results showed that in terms of temporal variation，the average NDVI of vegetation in Yanbian Prefecture from 2000 to 2022 was 0.539，showing a fluctuating upward trend with an upward slope of 0.003 6 per year.There were differences in vegetation NDVI among various counties and cities，but overall they showed a similar fluctuating upward trend.Among them，Helong City had the highest annual NDVI value at 0.560，while Yanji City had the lowest at 0.515.The annual average increase trend of NDVI in Tumen City was the most obvious，at 0.004 3/year，while the annual average increase trend of NDVI in Helong City was the slowest，at 0.003 1/year.In terms of spatial distribution，the overall vegetation coverage in Yanbian area was good，with 22.93% of the total area having an annual average NDVI greater than 0.6，ranging from 0.4 to 0.6，accounting for 67.16%;the area with vegetation NDVI less than 0.4 accounted for only 9.91%，mainly distributed in the centers and surrounding areas of various counties and cities.The area with vegetation NDVI showing an upward trend accounted for 94.23% of the total area of the study area，of which 71.8% showed a highly significant upward trend，and only 2.48% showed a downward trend，sporadically distributed in some cultivated land areas.In terms of the response of vegetation NDVI to climate factors，in most areas of Yanbian Prefecture，NDVI was positively correlated with temperature and precipitation，with about 67.23% of areas showing a significant positive correlation with precipitation，which was greater than the significant positive correlation with temperature （59.49%）.This indicated that NDVI was mainly positively promoted with precipitation and temperature throughout the state，and the correlation with precipitation was stronger than that with temperature.

Key words NDVI;Spatiotemporal variation characteristics;Climate factors;Response level;Yanbian Prefecture

作者簡介 王警潔（1983—），女，內(nèi)蒙古寧城人，工程師，從事氣象服務(wù)工作與研究。

收稿日期 2023-12-26

植被作為陸地生態(tài)系統(tǒng)的主體，在調(diào)節(jié)碳平衡、物質(zhì)能量循環(huán)及維持生態(tài)系統(tǒng)多樣性等方面發(fā)揮著重要作用，是全球氣候變化的“指示器”［1-2］。分析長時間序列的植被指數(shù)變化是監(jiān)測植被覆蓋變化的主要手段。在眾多植被指數(shù)中，歸一化植被指數(shù)（normalized difference vegetation index，NDVI）被證實能夠真實評價植被覆蓋狀況，已被廣泛應(yīng)用于全球和區(qū)域尺度的植被覆蓋變化監(jiān)測。近年來，國內(nèi)外眾多學(xué)者基于長時間序列NDVI數(shù)據(jù)源對不同區(qū)域的植被覆蓋動態(tài)變化及其影響因素進行了研究［3-15］。如毛德華等［3］利用1982—2006年GIMMS NDVI數(shù)據(jù)和MODIS NDVI數(shù)據(jù)分析了東北多年凍土區(qū)植被覆蓋的時空特征，并結(jié)合同期氣象觀測數(shù)據(jù)和土地利用/覆被數(shù)據(jù)，探討凍土區(qū)植被覆蓋對氣候變化和土地覆被變化的響應(yīng)，發(fā)現(xiàn)植被NDVI與氣溫和降水均具有較顯著的相關(guān)性，不同的土地利用類型間的轉(zhuǎn)變對植被NDVI的大小和空間分布產(chǎn)生了重要影響；劉家福等［4］基于MODIS NDVI和GIMMS NDVI數(shù)據(jù)集，研究了東北黑土區(qū)1982—2016年植被生長動態(tài)，指出平原地區(qū)植被NDVI與氣溫主要呈顯著正相關(guān)，而平原地區(qū)邊緣和山地地區(qū)的植被NDVI與降水量呈顯著正相關(guān)；石淞等［5］分析了2000—2020年東北地區(qū)植被時空演變特征及其對氣候變化的響應(yīng)，指出近年來東北地區(qū)NDVI均呈上升趨勢，與氣溫、降水呈正相關(guān)，其中降水對NDVI的影響作用最強。上述學(xué)者的研究多集中在省域以上的大尺度范圍，但對延邊地區(qū)小尺度近20年植被指數(shù)評估的研究鮮見報道。基于此，筆者利用2000—2022年長時間序列MOD13Q1數(shù)據(jù)，以整個延邊朝鮮族自治州（以下簡稱延邊州）作為研究對象，分析植被指數(shù)NDVI的時空變化和氣溫、降水量對其的響應(yīng)程度，為該區(qū)域生態(tài)與環(huán)境保護、國土資源的開發(fā)建設(shè)提供科學(xué)依據(jù)。

1 資料與方法

1.1 數(shù)據(jù)源及預(yù)處理

1.1.1 遙感數(shù)據(jù)。2000—2022年NDVI數(shù)據(jù)來源于美國國家航空航天局的MOD13Q1產(chǎn)品，時間和空間分辨率分別為16 d和250 m。使用MRT軟件對MOD13Q1時間序列數(shù)據(jù)進行鑲嵌、重采樣、投影轉(zhuǎn)換和格式轉(zhuǎn)換等批處理；應(yīng)用最大值合成法得到月NDVI 時間序列，利用均值法獲取2000—2022年年平均NDVI時間序列。

1.1.2 矢量文件。來源于全國地理信息資源目錄服務(wù)系統(tǒng)官網(wǎng)1∶100萬矢量數(shù)據(jù)，數(shù)據(jù)來源可靠，精度高。

1.1.3 氣象數(shù)據(jù)。來源于氣象大數(shù)據(jù)云平臺“天擎”，選取2000—2022年延邊州9個氣象站點的月平均氣溫和月降水量數(shù)據(jù)，通過克里金插值方法生成柵格數(shù)據(jù)集。

1.2 研究方法 該研究采用一元線性回歸法計算2000—2022年延邊州植被NDVI時間變化特征。采用Theil-Sen Median趨勢分析和Mann-Kendall檢驗法判斷延邊州植被NDVI時間序列變化趨勢和顯著程度。采用皮爾遜相關(guān)分析法估算延邊州NDVI與降水量和氣溫的相關(guān)性，從而分析氣候因素對全州植被NDVI變化的響應(yīng)。

2 結(jié)果與分析

2.1 植被NDVI時間變化特征 采用一元線性回歸方法對延邊州年均植被NDVI時間序列進行趨勢分析，結(jié)果如圖1和表1所示。從圖1可以看出，延邊州及各縣（市）NDVI年均值存在差異，但整體波動趨勢相似。2000—2022年延邊州植被NDVI多年平均值為0.539，最大值出現(xiàn)在2019年，為0.595，最小值出現(xiàn)在2010年，為0.497。和龍市NDVI年均值最高，為0.560；汪清縣次之，為0.559；延吉市最低，為0.515。

從表1可以看出，延邊州及各縣（市）NDVI年均值均呈現(xiàn)波動上升趨勢，延邊州上升斜率為0.003 6/年。圖們市植被NDVI上升趨勢最為明顯，上升斜率為0.004 3/年，汪清縣植被NDVI上升趨勢次之，上升斜率為0.004 0/年，和龍市植被NDVI上升趨勢最慢，上升斜率為0.003 1/年。

2.2 植被NDVI空間分布特征 利用ArcGIS軟件計算2000—2022年延邊州NDVI多年均值，得到延邊州植被NDVI多年均值空間分布，如圖2所示。延邊州植被覆蓋狀況整體較好，但植被NDVI在空間上呈現(xiàn)明顯地域差異，大部分地區(qū)的植被NDVI年均值較高，局部地區(qū)植被NDVI年均值較低。從圖2可以看出，植被NDVI年均值大于0.6的區(qū)域面積占總面積的22.93%，主要分布在汪清中部、敦化北部、安圖南部等區(qū)域；延邊州大部分地區(qū)植被NDVI年均值在0.4～0.6，其面積占總面積的67.16%；植被NDVI年均值小于0.4的面積僅占9.91%，主要分布在各縣（市）中心及周邊區(qū)域。由此可知，植被NDVI低值區(qū)主要分布在城市化水平較高的市中心及其周邊地區(qū)，而植被NDVI高值區(qū)主要位于林區(qū)，原因可能與以上地區(qū)的土地植被覆蓋類型、人類活動強度以及城市擴張速度有關(guān)。

2.3 植被NDVI空間變化趨勢 采用Theil-Sen Median趨勢分析方法，基于像元尺度計算植被NDVI的變化趨勢，揭示延邊州植被NDVI變化的空間差異。結(jié)果如圖3所示，2000—2022年延邊州植被NDVI變化斜率為-0.027～0.020/年。整體上，植被NDVI呈上升趨勢的面積遠遠大于呈下降趨勢的面積，東部呈上升趨勢的面積大于西部地區(qū)，南部呈上升趨勢的面積大于北部地區(qū)。

應(yīng)用Mann-Kendall檢驗法，基于像元尺度，計算植被NDVI變化趨勢的顯著性，并對趨勢類別進行統(tǒng)計，結(jié)果如圖4和表2所示。植被NDVI呈上升趨勢的面積占研究區(qū)總面積94.23%，其中，呈極顯著上升趨勢的面積占71.80%，植被NDVI無變化的面積占研究區(qū)總面積的3.29%，呈下降趨勢的面積僅占研究區(qū)面積的2.48%，零星分布在部分耕地區(qū)域。以上結(jié)果表明，23年間延邊州整體植被覆蓋明顯向好，這與延邊州近年來實施的天然林保護工程、退耕還林還（草）工程、國家重點生態(tài)功能區(qū)的建設(shè)等有很大的相關(guān)性。

2.4 植被NDVI和氣候因子相關(guān)性分析 采用Pearson相關(guān)分析估算植被NDVI與氣候因子（年平均氣溫、年降水量）的相關(guān)性，并進行顯著性檢驗，進而分析氣候因素對延邊州植被NDVI的影響。

從NDVI年均值與年平均氣溫的相關(guān)系數(shù)及顯著性檢驗空間分布（圖5）可以看出，延邊州植被NDVI與氣溫在空間上呈現(xiàn)正相關(guān)趨勢，其中呈顯著正相關(guān)的面積占總面積的

59.49%，主要集中在敦化北部和南部、安圖北部、汪清中部等地區(qū)，說明氣溫是影響上述地區(qū)NDVI的主要因素之一。NDVI與氣溫在空間上呈不顯著正相關(guān)的面積占總面積的38.75%，主要分布在琿春大部、敦化中部、延龍圖交界地區(qū)、安圖西南部等地區(qū)。NDVI與氣溫在空間上呈負相關(guān)的面積僅占總面積的1.76%，且分布離散。

從NDVI年均值與年降水量的相關(guān)系數(shù)及顯著性檢驗空間分布（圖6）可以看出，研究區(qū)大部分地區(qū)NDVI與降水量呈正相關(guān)趨勢，其中呈顯著正相關(guān)的區(qū)域主要分布在琿春密江、英安、春化、汪清羅子溝、天橋嶺、春陽、大興溝等地區(qū)，其面積占總面積的67.23%，說明降水量是影響上述地區(qū)NDVI的主要因素之一。不顯著正相關(guān)的面積占總面積的30.18%，主要集中在敦化額穆、雁鳴湖、大石頭、大蒲柴、和龍西城、安圖二道白河等地區(qū)。NDVI與降水量在空間上呈負相關(guān)的面積僅占總面積的2.59%，零星分布在敦化額穆北部、江南、翰章、安圖二道白河中部、琿春三家子、哈達門鄉(xiāng)等地區(qū)。

綜上所述，延邊州大部分地區(qū)NDVI與氣溫、降水量呈正相關(guān)，約有67.23%的地區(qū)與降水量呈顯著正相關(guān)，大于與氣溫呈顯著正相關(guān)區(qū)域（59.49%），表明全州范圍內(nèi)NDVI與降水量和氣溫均以正向促進為主，且與降水量的相關(guān)性強于氣溫。

3 結(jié)論

該研究利用一元回歸分析法、Theil-Sen Median趨勢分析、Mann-Kendall顯著性檢驗和相關(guān)系數(shù)法，基于2000—2022年的長時間序列MOD13Q1數(shù)據(jù)，揭示了延邊州植被NDVI時空演變特征，并探討其與氣溫、降水量等氣候因子之間的關(guān)系。主要研究結(jié)論如下：

（1）時間變化上來看，2000—2022年延邊州植被NDVI多年平均值為0.539，呈波動上升趨勢，上升斜率為0.003 6/年。各縣（市）植被NDVI年均值存在差異，但總體均呈現(xiàn)相似的波動上升趨勢。其中，和龍市NDVI年均值最高，為0.560，延吉市最低，為0.515。圖們市NDVI年均值上升趨勢最為明顯，上升斜率為0.004 3/年，和龍市NDVI年均值上升趨勢最慢，上升斜率為0.003 1/年。

（2）空間分布上來看，延邊州植被整體覆蓋狀況較好，植被NDVI年均值大于0.6的區(qū)域面積占總面積的22.93%，在0.4～0.6的占67.16%；植被NDVI年均值小于0.4的面積僅占9.91%，主要分布在各縣（市）中心及周邊區(qū)域。植被NDVI呈上升趨勢的面積占研究區(qū)總面積94.23%，其中，呈極顯著上升趨勢的面積占71.80%，呈下降趨勢的面積僅占研究區(qū)面積的2.48%，零星分布在部分耕地區(qū)域。

（3）植被NDVI對氣候因子響應(yīng)方面，延邊州大部分地區(qū)NDVI與氣溫、降水量呈正相關(guān)，約有67.23%的地區(qū)與降水量呈顯著正相關(guān)，大于與氣溫呈顯著正相關(guān)區(qū)域（59.49%），表明全州范圍內(nèi)NDVI與降水量和氣溫均以正向促進為主，且與降水量的相關(guān)性強于氣溫。

參考文獻

［1］

龐鑫，劉珺.氣候變化對亞洲地區(qū)植被NDVI變化的影響［J］.自然資源遙感，2023，35（2）：295-305.

［2］ 崔林麗，史軍，肖風(fēng)勁，等.中國東部NDVI的變化趨勢及其與氣候因子的相關(guān)分析［J］.資源科學(xué)，2010，32（1）：124-131.

［3］ 毛德華，王宗明，宋開山，等.東北多年凍土區(qū)植被NDVI變化及其對氣候變化和土地覆被變化的響應(yīng)［J］.中國環(huán)境科學(xué)，2011，31（2）：283-292.

［4］ 劉家福，馬帥，李帥，等.1982—2016年東北黑土區(qū)植被NDVI動態(tài)及其對氣候變化的響應(yīng)［J］.生態(tài)學(xué)報，2018，38（21）：7647-7652.

［5］ 石淞，李文，丁一書，等.東北地區(qū)植被時空演變及影響因素分析［J］.中國環(huán)境科學(xué)，2023，43（1）：276-289.

［6］ 王景發(fā).吉林省NDVI時空變化：以2010—2019年為例［J］.農(nóng)業(yè)與技術(shù)，2020，40（13）：103-106.

［7］ 桑婧.1981—2017年黃河流域內(nèi)蒙古段氣候及植被生態(tài)變化狀況［J］.內(nèi)蒙古氣象，2021（1）：33-38.

［8］ 程兀杰，孟妮娜，蔡昕楠，等.陜西省NDVI時空變化及其對氣候和人類活動的響應(yīng)［J］.人民黃河，2023，45（4）：28-34.

［9］ 林楠，姜琦剛，陳永良.近20年來延邊州NDVI變化與氣候關(guān)聯(lián)分析研究［C］//中國地球物理學(xué)會.中國地球物理·2009.合肥：中國科學(xué)技術(shù)大學(xué)出版社，2009：709-710.

［10］ 董葉輝，南穎，劉志鋒，等.近十年延邊朝鮮族自治州植被變化特征分析［J］.延邊大學(xué)學(xué)報（自然科學(xué)版），2010，36（3）：269-274.

［11］ 劉琪，景海濤，劉盼盼，等.太行山區(qū)NDVI時空變化及其與氣候因子的關(guān)系［J］.西北林學(xué)院學(xué)報，2023，38（1）：1-9，33.

［12］ 金麗娜，李雄飛，曲靜，等.秦嶺山區(qū)NDVI與氣候因子時空變化特征及NDVI對氣候的響應(yīng)［J］.甘肅科學(xué)學(xué)報，2022，34（3）：65-72.

［13］ 蒲泓君，邵懷勇，賴云波，等.四川省2000～2015年NDVI動態(tài)變化及其對氣候波動的響應(yīng)［J］.四川環(huán)境，2022，41（6）：107-115.

［14］ 張曉明.江西省NDVI時空變化特征及其與氣候因子相關(guān)性分析［J］.長江信息通信，2022，35（6）：7-12.

［15］ 黃啟芬，張娜，陳澤，等.2000—2018年貴州植被NDVI對氣候變化和人類活動的響應(yīng)［J］.貴州科學(xué)，2023，41（2）：48-53.