青土湖綠洲植被指數(shù)適用性及優(yōu)選研究

楊麗娟 黃峰 潘若云

摘要 [目的]為精準(zhǔn)評(píng)價(jià)青土湖綠洲植被恢復(fù)情況，開展植被指數(shù)適用性及優(yōu)選研究。[方法]基于2013—2019年植被生長(zhǎng)旺季的Landsat8-OLI遙感數(shù)據(jù)，分析歸一化植被指數(shù)、比值植被指數(shù)、大氣阻抗植被指數(shù)等8種常用植被指數(shù);采用趨勢(shì)分析法研究各植被指數(shù)的變化特征，通過(guò)評(píng)價(jià)植被指數(shù)時(shí)空變化特征與綠洲實(shí)際恢復(fù)特征的一致性研究各植被指數(shù)的適用性;基于各植被指數(shù)的變化率及趨勢(shì)顯著性，評(píng)價(jià)不同植被指數(shù)檢測(cè)青土湖綠洲恢復(fù)的敏感性。[結(jié)果]增強(qiáng)植被指數(shù)、歸一化植被指數(shù)、土壤調(diào)節(jié)指數(shù)和修改型土壤調(diào)節(jié)指數(shù)對(duì)青土湖綠洲植被覆蓋變化相對(duì)敏感，是適用于青土湖綠洲、能有效檢測(cè)綠洲植被覆蓋變化的植被指數(shù)。[結(jié)論]研究成果可為青土湖生態(tài)輸水調(diào)控及綠洲恢復(fù)管理提供一定的技術(shù)支撐和科學(xué)參考。

關(guān)鍵詞 植被指數(shù);適用性;敏感性;生態(tài)輸水;青土湖

中圖分類號(hào) Q948 ?文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2020）19-0065-04

doi：10.3969/j.issn.0517-6611.2020.19.018

Abstract [Objective] In order to evaluate the vegetation restoration in Qingtu Oasis accurately，the study on the applicability and optimization of vegetation index in Qingtu Oasis was carried out.[Method]Based on the Landsat8OLI remote sensing data of the flourishing vegetation season from 2013 to 2019，the normalized difference vegetation index，ratio vegetation index，atmospherically resistant vegetation index and the other 5 vegetation indexes were analyzed.Then the trend analysis was used to study the variation characteristics of each vegetation index.The applicability of each vegetation index was studied by evaluating the consistency between the temporal and spatial variation characteristics of vegetation index and the actual recovery characteristics of oasis.Based on the change rate and trend significance of each vegetation index，the sensitivity of different vegetation indexes to detect the recovery of Qingtu Lake Oasis was evaluated.[Result]The enhanced vegetation index，normalized difference vegetation index，soiladjusted vegetation index and modified soiladjusted vegetation index were relatively sensitive to the change of vegetation cover in Qingtu Oasis，which were suitable for effectively detecting the change of vegetation cover in Qingtu Oasis.[Conclusion]The research results can provide certain technical support and scientific reference for ecological water conveyance regulation and oasis restoration management of the Qingtu Lake.

Key words Vegetation index;Applicability;Sensitivity;Ecological water conveyance;Qingtu Lake

基金項(xiàng)目 國(guó)家重點(diǎn)研發(fā)計(jì)劃（2017YFC0406101）。

作者簡(jiǎn)介 楊麗娟（1999—），女，江蘇南通人，從事水務(wù)工程研究。*通信作者，副教授，博士，碩士生導(dǎo)師，從事文水資源、生態(tài)水文研究。

收稿日期 2020-03-31

植被指數(shù)（vegetation index，VI）是對(duì)地表植被狀況的簡(jiǎn)單有效的度量，也是評(píng)價(jià)干旱區(qū)生態(tài)輸水效應(yīng)的常用指標(biāo)[1]。當(dāng)前已有歸一化植被指數(shù)、土壤調(diào)節(jié)植被指數(shù)、比值植被指數(shù)等上百種植被指數(shù)產(chǎn)品應(yīng)用于科研生產(chǎn)實(shí)踐，不同的植被指數(shù)由于地理環(huán)境、植被覆蓋度、植被生長(zhǎng)狀況等因素的影響，對(duì)于不同的地區(qū)有不同的適用性和敏感性[2] 。郭玉川等[3] 研究了歸一化植被指數(shù)、增強(qiáng)植被指數(shù)等4種指數(shù)反演植被覆蓋度的精度，建議選擇歸一化植被指數(shù)以評(píng)估塔里木河生態(tài)輸水后的生態(tài)恢復(fù)程度;林茜等[4] 在塔里木河流域通過(guò)評(píng)估歸一化植被指數(shù)、比值植被指數(shù)等和植被覆蓋度的相關(guān)性，建議選擇歸一化植被指數(shù)以獲取植被信息。

我國(guó)河西走廊石羊河流域尾閭青土湖自生態(tài)輸水以來(lái)，水域及綠洲面積在一定程度上得到恢復(fù)，整體植被覆蓋度有所增加。

鑒于不同植被指數(shù)在不同區(qū)域識(shí)別植被特征的適用性差異以及檢測(cè)植被覆被變化的敏感性差異，為科學(xué)準(zhǔn)確地評(píng)估石羊河流域生態(tài)輸水對(duì)其尾閭青土湖綠洲的恢復(fù)效應(yīng)，開展植被指數(shù)適用性及優(yōu)選研究。筆者選取在干旱、半干旱地區(qū)較為常用的8種植被指數(shù)開展比較研究，遴選適用于檢測(cè)青土湖地區(qū)植被覆蓋變化的指數(shù)，為生態(tài)輸水調(diào)控提供參考，也為其他流域合理選取植被指數(shù)、科學(xué)監(jiān)測(cè)植被覆蓋變化提供借鑒。

1 資料與方法

1.1 研究區(qū)概況

青土湖綠洲位于石羊河流域尾閭（圖1），地處甘肅省民勤縣東北部，在巴丹吉林沙漠東南部、騰格里沙漠西部，地理位置39°04′～39°09′N、103°36′～103°39′E[5] 。青土湖地區(qū)年均氣溫7.8 ℃，降水稀少且主要集中在7—9月，多年平均降水量?jī)H110 mm[6] ，全年盛行西北風(fēng)，屬于溫帶大陸性干旱氣候[7] 。

青土湖曾是石羊河流域最大的湖泊，水域面積超過(guò)400 km2。后來(lái)因氣候變化和流域上、中游不合理用水等因素影響，湖泊面積開始萎縮，加上紅崖山水庫(kù)的修建，青土湖地表水補(bǔ)給源逐漸消失，于1959年徹底干涸[8] 。青土湖干涸導(dǎo)致民勤地區(qū)失去了天然的生態(tài)屏障，植被減少、荒漠化加劇，生態(tài)系統(tǒng)嚴(yán)重退化[6] 。為促進(jìn)當(dāng)?shù)厣鷳B(tài)恢復(fù)，2010年9月石羊河流域綜合治理工程有計(jì)劃地向下游生態(tài)輸水，青土湖開始形成人工季節(jié)性水面，區(qū)域環(huán)境質(zhì)量明顯好轉(zhuǎn)[9] 。

1.2 數(shù)據(jù)選取與方法

1.2.1 遙感影像。

選取青土湖生態(tài)輸水后2013—2019年逐年的Landsat8-OLI遙感數(shù)據(jù)，下載于美國(guó)地質(zhì)調(diào)查局提供的網(wǎng)站https：//earthexplorer.usgs.gov/。衛(wèi)星訪問(wèn)地球的周期為16 d，TM1～TM9的9個(gè)波段空間分辨率均為30 m。根據(jù)青土湖綠洲的經(jīng)緯度范圍，選擇的OLI數(shù)據(jù)行號(hào)為131，列號(hào)為33。綜合考慮云量信息和植被信息，收集整理7、8月份的遙感影像資料。利用ENVI 5.3軟件預(yù)處理遙感影像，主要包括輻射定標(biāo)、大氣校正、圖像裁剪、去除異常值等。

1.2.2 植被指數(shù)。

選取在干旱、半干旱地區(qū)較為常用的歸一化植被指數(shù)、比值植被指數(shù)、大氣阻抗植被指數(shù)等8種植被指數(shù)，其計(jì)算公式如表1所示，其中，NIR、RED、GREEN和BLUE分別是近紅外波段（0.845～0.885 μm）、紅光波段（0.630～0.680 μm）、綠光波段（0.525～0.600 μm）和藍(lán)光波段（0.450～0.453 μm）的反射率;L是土壤亮度校正因子，隨綠色植被數(shù)量的變化而變化，采用Huete[13]建議的L最佳取值0.5;C1、C2是大氣修正參數(shù);l為土壤調(diào)節(jié)參數(shù)，取值分別為6.0、7.5和1.0[17] ;γ為光學(xué)路徑的效應(yīng)因子，取決于氣溶膠的類型，采用Kaufman等[12]推薦的取值1。NDVI、ARVI、SAVI、MSAVI2、NDGI和EVI取值為-1～1，非正值代表裸地或者水體，正值表示有綠色植被，且越接近1表示植被覆蓋度越大[12-15] 。RVI取值在0～2和8～30時(shí)無(wú)綠色植被覆蓋，在2～8時(shí)表示有植被覆蓋，且數(shù)值越接近8表示植被覆蓋度越大[11] 。DVI取值隨植被覆蓋的增加而增加[17]。

1.2.3 適用性及敏感性評(píng)價(jià)。

針對(duì)各個(gè)植被指數(shù)，采用一元線性回歸法逐像元擬合其變化趨勢(shì)，以綜合反映研究區(qū)域植被的時(shí)空演變情況。根據(jù)一元線性擬合線的斜率s與回歸分析顯著性P值將變化趨勢(shì)分為5個(gè)等級(jí)：顯著下降（s<0，P≤0.01）、下降（s<0，0.010.05）、上升（s>0，0.010，P≤0.01）[18] 。采用Kappa系數(shù)衡量不同植被指數(shù)分析成果的一致性，識(shí)別性能相近的植被指數(shù);通過(guò)對(duì)比分析各植被指數(shù)所表征的植被覆蓋變化趨勢(shì)與植被實(shí)際動(dòng)態(tài)變化特征的一致性，評(píng)估植被指數(shù)的適用性。由Kappa系數(shù)值表示一致性的程度：一致性程度極差（Kappa<0）、一致性程度微弱（0≤Kappa≤0.20）、一致性程度弱（0.20

2 結(jié)果與分析

2.1 植被指數(shù)適用性分析

生態(tài)輸水以來(lái)，青土湖綠洲植被呈較為顯著的恢復(fù)趨勢(shì)。植被覆蓋區(qū)域由湖區(qū)向外緣荒漠區(qū)擴(kuò)展，總體呈擴(kuò)大趨勢(shì);距離湖區(qū)越近，植被覆蓋密度增加越顯著;距離湖區(qū)較遠(yuǎn)的區(qū)域植被覆蓋密度增加不明顯[20-21] 。NDVI、RVI、SAVI、EVI和MSAVI2的分析成果與青土湖綠洲實(shí)際恢復(fù)情況較為一致（圖2），各指標(biāo)分析成果之間的Kappa系數(shù)都在0.6～1.0，一致性程度達(dá)到高、極高水平（表2）。ARVI、NDGI、DVI的分析成果與青土湖綠洲實(shí)際恢復(fù)情況差異較大（圖2），且ARVI、NDGI、DVI、NDVI、RVI、SAVI、EVI、MSAVI2分析成果的Kappa系數(shù)都低于0.6，一致性較低，表明ARVI、NDGI、DVI在青土湖綠洲的適用性相對(duì)較差。其中，ARVI反映植被覆蓋在外圍荒漠區(qū)呈顯著下降趨勢(shì)，NDGI反映83%的區(qū)域植被覆蓋呈顯著下降趨勢(shì)，DVI在植被覆蓋較高的水域附近無(wú)明顯變化。

2.2 植被指數(shù)敏感性分析

各植被指數(shù)空間均值序列趨勢(shì)分析成果（圖3）表明，各指數(shù)均呈顯著上升趨勢(shì)。其中，EVI、NDVI、SAVI、MSAVI2的趨勢(shì)變化率基本一致，RVI、ARVI、DVI、NDGI的趨勢(shì)變化率相對(duì)較低，表明EVI、NDVI、SAVI和MSAVI2檢測(cè)青土湖綠洲恢復(fù)變化特征的性能基本一致，RVI、ARVI、DVI和NDGI檢測(cè)青土湖綠洲恢復(fù)變化特征的敏感性相對(duì)較低。

由于DVI對(duì)土壤背景的變化敏感，比較適合在低植被覆蓋區(qū)監(jiān)測(cè)生態(tài)環(huán)境的變化。當(dāng)植被覆蓋度較低時(shí)，DVI隨生物量的增加而增加，較高時(shí)DVI對(duì)植被的靈敏度下降。青土湖湖區(qū)附近自生態(tài)輸水以來(lái)植被覆蓋高且增加趨勢(shì)顯著[20] ，DVI在此情況下易飽和，因此DVI對(duì)湖區(qū)附近的敏感性較低。ARVI使用藍(lán)光波段和紅光波段的組合矯正大氣散射的影響，在該研究中外圍荒漠低植被覆蓋區(qū)的應(yīng)用效果相對(duì)較差，可能在一定程度上受裸地背景的影響。NDGI選用紅光波段和綠光波段的亮度值作為波段組合的信息源，兩者亮度差異值較小，應(yīng)用于檢測(cè)青土湖綠洲植被恢復(fù)的效果相對(duì)較差。因此DVI、ARVI和NDGI三者在青土湖綠洲的適用性和敏感性不足。

3 結(jié)論

生態(tài)輸水以來(lái)，青土湖綠洲的植被覆蓋區(qū)域由湖區(qū)向周圍荒漠?dāng)U大，總體覆蓋密度呈增加趨勢(shì)。NDVI、RVI、SAVI、EVI和MSAVI2能客觀反映這一變化趨勢(shì)，應(yīng)用于評(píng)價(jià)青土湖綠洲生態(tài)恢復(fù)的適用性較好;ARVI、NDGI、DVI未能客觀反映青土湖綠洲恢復(fù)特征，其適用性相對(duì)較差。綜合植被指數(shù)空間均值變化率及趨勢(shì)顯著性評(píng)價(jià)，EVI、NDVI、SAVI、MSAVI2檢測(cè)青土湖綠洲植被覆蓋變化的敏感性相對(duì)較高，RVI的敏感性相對(duì)較弱。建議在青土湖生態(tài)恢復(fù)評(píng)價(jià)工作中選用EVI、NDVI等植被指數(shù)。

參考文獻(xiàn)

[1] 郭鈮.植被指數(shù)及其研究進(jìn)展[J].干旱氣象，2003，21（4）：71-75.

[2] 賀金蘭，張明月，田尉霞.基于遙感影像的不同植被指數(shù)比較研究[J].科技創(chuàng)新與應(yīng)用，2017（6）：36-37.

[3] 郭玉川，何英，李霞.基于MODIS的干旱區(qū)植被覆蓋度反演及植被指數(shù)優(yōu)選[J].國(guó)土資源遙感，2011（2）：115-118.

[4] 林茜，郭飛，黃昌春，等.基于MODIS的植被指數(shù)對(duì)估測(cè)塔里木河中下游植被蓋度的適用性[J].草業(yè)科學(xué)，2016，33（12）：2434-2441.

[5] 張華，張玉紅，張改改.民勤綠洲青土湖植被優(yōu)勢(shì)種地上生物量估算[J].干旱區(qū)地理，2020，43（1）：201-210.

[6] 陳政融，劉世增，劉淑娟，等.蘆葦和白刺空間格局對(duì)青土湖生態(tài)輸水的響應(yīng)[J].草業(yè)科學(xué)，2015，32（12）：1960-1968.

[7] 趙強(qiáng)，王乃昂，李秀梅，等.青土湖地區(qū)9500a BP以來(lái)的環(huán)境變化研究[J].冰川凍土，2005，27（3）：352-359.

[8] YANG X M，LIU S Z，YANG T B，et al.Spatialtemporal dynamics of desert vegetation and its responses to climatic variations over the last three decades：A case study of Hexi region in Northwest China[J].Journal of arid land，2016，8（4）：556-568.

[9] 劉淑娟，袁宏波，劉世增，等.石羊河尾閭水面形成區(qū)土壤顆粒的分形特征[J].水土保持通報(bào)，2013，33（6）：285-289.

[10] ROUSE J W，JR，HAAS R H，SCHELL J A，et al.Monitoring vegetation systems in the Great Plains with ERTS[C]//Proceedings of third earth resources technology satellite1 symposium.Washington：NASA，1974：310-317.

[11] JORDAN C F.Derivation of leafarea index from quality of light on the foresr floor[J].Ecology，1969，50（4）：663-666.

[12] KAUFMAN Y J，TANRE D.Atmospherically resistant vegetation index （ARVI） for EOSMODIS[J].IEEE Transactions on Geoscience and Remote Sensing，1992，30（2）：261-270.

[13] HUETE A R.A Soil Adjusted Vegetation Index（SAVI）[J].Remote Sens Environ，1988，25（3）：295-309.

[14] JIANG Z Y，HUETE A R，LI J，et al.Interpretation of the modified soiladjusted vegetation index isolines in redNIR reflectance space[J].Journal of applied remote sensing，2007，1：1-12.

[15] RHEW I C，VANDER STOEP A，KEARNEY A，et al.Validation of the normalized difference vegetation index as a measure of neighborhood greenness[J].Annals of epidemiology，2011，21（12）：946-952.

[16] BECKER F，CHOUDHURY B J.Relative sensitivity of normalized difference vegetation Index （NDVI） and microwave polarization difference Index（MPDI） for vegetation and desertification monitoring[J].Remote sensing of environment，1988，24（2）：297-311.

[17] LIU H Q，HUETE A R.A feedback based modification of the NDVI to minimize canopy background and atmosphere noise[J].IEEE Trans Geosci Remote Sensing，1995，33：457-465.

[18] WOLTERS M，GARBUTT A，BAKKER J P.Saltmarsh restoration：Evaluating the success of deembankments in northwest Europe[J].Biological conservation，2005，123（2）：249-268.

[19] LANDIS J R，KOCH G G.The measurement of observer agreement for categorical data[J].Biometrics，1977，33（1）：159-174.

[20] 趙軍，楊建霞，朱國(guó)鋒.生態(tài)輸水對(duì)青土湖周邊區(qū)域植被覆蓋度的影響[J].干旱區(qū)研究，2018，35（6）：1251-1261.

[21] CHUNYU X Z，HUANG F，XIA Z Q，et al.Assessing the ecological effects of water transport to a lake in arid regions：A case study of Qingtu Lake in Shiyang River basin，Northwest China[J].International journal of environmental research and public health，2019，16（1）：1-14.