白雪蓮,喬榮榮,季樹(shù)新,閆志堅(jiān),常學(xué)禮,*,趙文智

1 中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院,中國(guó)生態(tài)系統(tǒng)研究網(wǎng)絡(luò)臨澤內(nèi)陸河流域研究站,中國(guó)科學(xué)院內(nèi)陸河流域生態(tài)水文重點(diǎn)實(shí)驗(yàn)室,蘭州 730000 2 魯東大學(xué)資源與環(huán)境工程學(xué)院,煙臺(tái) 264025 3 中國(guó)農(nóng)業(yè)科學(xué)院草原研究所,呼和浩特 010010 4 中國(guó)科學(xué)院大學(xué),北京 100049

非光合植被(Non-photosynthetic Vegetation, NPV)是陸地生態(tài)系統(tǒng)重要組成部分,是衡量地表植被覆蓋狀況的重要指標(biāo)之一[1]。NPV主要指枯草、枯葉、枯枝、枯干、作物留茬和凋落物以及植物活體非光合部分等[2],是生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動(dòng)重要環(huán)節(jié)之一,特別是枯死部分不僅能夠增加土壤有機(jī)質(zhì)改善土壤質(zhì)地,還能夠有效減少地表徑流和土壤風(fēng)蝕[3—5],并在草原火災(zāi)監(jiān)測(cè)方面起著關(guān)鍵作用[6—7],及時(shí)、準(zhǔn)確的掌握NPV動(dòng)態(tài)對(duì)于水土保持、風(fēng)沙天氣的減緩、沙塵源的識(shí)別和草原生態(tài)系統(tǒng)保護(hù)有至關(guān)重要的意義。

隨著遙感技術(shù)提高及易操作、低成本和可周期性重復(fù)特點(diǎn)使其被廣泛應(yīng)用到區(qū)域植被動(dòng)態(tài)監(jiān)測(cè)[8—10],并應(yīng)用一些指數(shù)來(lái)表征農(nóng)田NPV取得了大量結(jié)果,如歸一化指數(shù)(Normalized Difference Index, NDI)、歸一化差異衰敗指數(shù)(Normalized Difference Senescent Vegetation Index, NDSVI)、土壤調(diào)整作物茬指數(shù)(Soil Adjusted Corn Residue Index, SACRI)和修正土壤調(diào)整作物茬指數(shù)(MSACRI)等,但是估算誤差較大[11—13]。與上述研究途徑不同,Elvidge(1988)利用機(jī)載可見(jiàn)紅外成像分光計(jì)2180—2220 nm和2310—2380 nm波段數(shù)據(jù)計(jì)算并提出了纖維素指數(shù)概念,同時(shí)還發(fā)現(xiàn)了灌木分布區(qū)反射光譜在2100 nm處有明顯光譜吸收槽[14]。Daughtry等于1996年在此基礎(chǔ)上提出了纖維吸收指數(shù)(Cellulose Absorption Index,CAI),涉及波段涵蓋了Elvidge的發(fā)現(xiàn),分別為2000—2050 nm,2080—2130 nm和2190—2240 nm[15]。CAI提出后,大量的驗(yàn)證與應(yīng)用研究表明與其他類(lèi)似指數(shù)相比,CAI可以更好反映非光合植被蓋度(fNPV)[16—20]。但是,在MODIS數(shù)據(jù)36個(gè)波段中,缺少與 CAI計(jì)算所對(duì)應(yīng)波段(2000—2050 nm,2080—2130 nm和2190—2240 nm),因此,如果能在MODIS數(shù)據(jù)中發(fā)現(xiàn)與CAI計(jì)算波長(zhǎng)密切相關(guān)波段成為NPV研究關(guān)鍵。針對(duì)這一問(wèn)題,Cao等根據(jù) MODIS波段特征提出DFI能很好反映地表枯落物狀態(tài),其研究結(jié)果在國(guó)際遙感雜志(International Journal of Remote Sensing)發(fā)表[18—19]。此后,相關(guān)研究直接將DFI應(yīng)用到fNPV估算中[11, 21—22],但是DFI是否可以準(zhǔn)確反應(yīng)NPV纖維素含量狀況？DFI計(jì)算所采用波段與CAI反應(yīng)敏感波長(zhǎng)存在偏差的事實(shí)如何補(bǔ)償?shù)冗@些細(xì)節(jié)問(wèn)題一直被忽視。因此,完成CAI和DFI關(guān)系分析是實(shí)現(xiàn)fNPV更加精確估算的關(guān)鍵。在已有研究中CAI和DFI的關(guān)系是通過(guò)構(gòu)建像元三分法進(jìn)行間接比較分析[2, 21, 23—24],缺乏對(duì)二者之間的直接實(shí)驗(yàn)驗(yàn)證。從NPV解譯方法來(lái)看,覆蓋度估算主要有非光合植被-裸地(NPV-BS)二分模型、光合植被-非光合植被-裸地(PV-NPV-BS)三分模型和光合植被-非光合植被-裸地-陰影(PV-NPV-BS-Shadow)四分模型等[2, 9, 21, 23]。這些方法在理論上較好詮釋了地表植被覆蓋構(gòu)成和外界因素對(duì)蓋度估算影響,但已有研究也明確指出了上述方法存在問(wèn)題[25—26],如像元三分和四分模型在植被枯黃期、PV蓋度(fPV)限制(<30%)條件下才能較好估算fNPV等。此外,一些相關(guān)研究中關(guān)注了裸地可能會(huì)對(duì)fNPV提取產(chǎn)生影響,但是在蓋度計(jì)算中對(duì)于是否應(yīng)該剔除掉裸地值缺少直接的對(duì)比分析[9, 12, 21, 23]。

從有關(guān)大尺度NPV研究現(xiàn)狀來(lái)看[18—19, 23—24, 26—27],目前處在將NPV作為一個(gè)均質(zhì)整體分析階段,對(duì)其由NPV蓋度和不同生活型植物凋落組成等導(dǎo)致空間異質(zhì)性缺乏考慮[2, 16, 27—29]。但是,已有研究表明灌木fNPV估算誤差高于草本和農(nóng)作物[30],明顯看出木質(zhì)材料和草本材料存在差異。中國(guó)對(duì)于NPV的研究主要集中于錫林郭勒草原,并把該區(qū)看作典型草原,忽略了該區(qū)草甸草原和沙地草場(chǎng)異質(zhì)性以及木質(zhì)NPV(灌木、半灌木)和草本NPV在光譜反射方面的差異,在一定程度上影響fNPV估算精度。另外,現(xiàn)在對(duì)于植被蓋度較低的荒漠草原fNPV研究較薄弱,有必要弄清此類(lèi)草原fNPV分布狀況。同時(shí),基于不同生活型植物NPV差異導(dǎo)致空間異質(zhì)性研究相對(duì)薄弱[16],尤其是對(duì)不同生活型植物fNPV、CAI和DFI關(guān)系尚不明朗,急需開(kāi)展此方面更加深入的研究。

因此,本文研究目的為：(1)構(gòu)建fNPV、CAI和DFI的關(guān)系,估算不同草原類(lèi)型fNPV,解釋不同草原(草場(chǎng))類(lèi)型fNPV空間異質(zhì)性特征；(2)對(duì)比DFI直接估算fNPV與CAI間接估算fNPV差異,弄清DFI是否可以代替CAI參與fNPV估算,提高fNPV估算精度；(3)對(duì)比剔除裸地前后的估算結(jié)果,解釋土壤背景值對(duì)fNPV估算模型和估算結(jié)果的影響；(4)為利用非生長(zhǎng)季遙感數(shù)據(jù)識(shí)別區(qū)域fNPV并評(píng)估區(qū)域土壤風(fēng)蝕、草地沙化程度識(shí)別等提供基礎(chǔ)支持。

1 研究區(qū)概況與方法

1.1 研究區(qū)概況

地面反射光譜測(cè)定在中國(guó)農(nóng)業(yè)科學(xué)院草原研究所鄂爾多斯沙地試驗(yàn)站完成(圖1),位于內(nèi)蒙古自治區(qū)鄂爾多斯市達(dá)拉特旗,地處庫(kù)布齊沙漠東段；地理位置為110°0′51″E,40°18′53″N；海拔高度為1040 m。氣候特點(diǎn)為冬季寒冷漫長(zhǎng),夏季溫暖短促,降水較少；年均溫度6℃,最低氣溫-32.3℃,最高氣溫38.3℃；年降水量250—350 mm之間,集中在7、8月,屬于溫帶大陸性氣候。主要植物種有灌木北沙柳(Salixpsammophila)和檸條錦雞兒(CaraganakorshinskiiKom.),半灌木油蒿(Artemisiaordosica),草本植物黃蒿(Artemisiascopariawaldst.etKit)、狗尾草(Setariaviridis)、沙蓬(Agriophyllumsquarrosum)和霧冰藜(Bassiadasyphylla)等。根據(jù)內(nèi)蒙古植被類(lèi)型圖和草原實(shí)地調(diào)查(圖1),案例分析區(qū)域分別選擇錫林郭勒盟以多年生草本為主的西烏珠穆沁旗草甸草原和以多年生+半灌木為主的阿巴嘎旗典型草原區(qū)、以草本和灌木混合為主的鄂爾多斯十大孔兌沙地草場(chǎng)以及以草本、半灌木和灌木混合為主的巴彥淖爾市烏拉特中旗荒漠草原。其中,為最大可能規(guī)避異質(zhì)性產(chǎn)生的干擾,在所選縣(旗)域中用2015年LUCC數(shù)據(jù)中的草地將典型區(qū)選出。

圖1 研究區(qū)位置Fig.1 Location of study area

1.2 研究方法

1.2.1實(shí)驗(yàn)材料

依據(jù)沙地草場(chǎng)植物生長(zhǎng)型組成特點(diǎn),分別選擇灌木(北沙柳、檸條錦雞兒)、半灌木(油蒿)和草本(黃蒿、狗尾草、沙蓬)三種類(lèi)型進(jìn)行反射光譜測(cè)定。

1.2.2地面反射光譜數(shù)據(jù)與遙感數(shù)據(jù)

地面反射光譜數(shù)據(jù)測(cè)定時(shí)間為2019年5月9日—13日,此時(shí)可減少PV干擾。選用美國(guó)ASD公司Field Spec 4高分辨率光譜儀,波長(zhǎng)范圍為 350—2500 nm,其中為減小裸地對(duì)反射率影響,每種植被光譜測(cè)定時(shí)采用同一裸地作為背景。測(cè)定時(shí),探頭垂直向下,視場(chǎng)角為 25°,與實(shí)測(cè)目標(biāo)地物距離為0.67 m,對(duì)應(yīng)地面樣本為直徑 30 cm 樣圓。光譜測(cè)定時(shí)間為10—14時(shí)之間,對(duì)同一個(gè)樣圓植物反射光譜測(cè)定5次,并用相匹配數(shù)據(jù)處理軟件ViewSpec Pro導(dǎo)出所測(cè)光譜值,對(duì)應(yīng)MODIS數(shù)據(jù)波段B1(620—670 nm)、B2(841—876 nm)、B6(1628—1652 nm)、B7(2105—2155 nm)和計(jì)算CAI所需R2.0(2000—2050 nm)、R2.1(2080—2130 nm)和R2.2(2190—2240 nm),計(jì)算平均值。為保證光譜測(cè)定準(zhǔn)確性,每隔10—15 min利用標(biāo)準(zhǔn)白板進(jìn)行校光。測(cè)定期間,光照條件較穩(wěn)定、風(fēng)力弱,對(duì)光譜反射值干擾極小。

選用MODIS數(shù)據(jù)進(jìn)行典型研究區(qū)fNPV估算,數(shù)據(jù)來(lái)源于NASA地表過(guò)程數(shù)據(jù)中心(LP DAAC),下載地址為https://ladsweb.modaps.eosdis.nasa.gov/,數(shù)據(jù)類(lèi)型為MCD43A4,該數(shù)據(jù)為基于16天觀測(cè)數(shù)據(jù)生成的日產(chǎn)品數(shù)據(jù),已經(jīng)經(jīng)過(guò)了地形校正、輻射定標(biāo)和大氣校正。在本研究中,選用2019年4月15日數(shù)據(jù),原始數(shù)據(jù)為HDF格式,通過(guò)MRT(Modis Reprojection Tool)進(jìn)行轉(zhuǎn)換,將1、2、6和7波段數(shù)據(jù)由HDF格式轉(zhuǎn)換為T(mén)IFF格式。

1.2.3fNPV提取

由于傳統(tǒng)目視估測(cè)植被蓋度比照相法精度低,所以植被蓋度測(cè)定采用相機(jī)拍照,然后利用 ENVI 5.3通過(guò)監(jiān)督分類(lèi)和目視解譯結(jié)合方法進(jìn)行蓋度提取。拍照時(shí),擋住陽(yáng)光減小陰影對(duì)蓋度解譯誤差,相機(jī)鏡頭垂直向下位于樣圓中心正上方,調(diào)整距離樣圓高度,保證樣圓完全處于相機(jī)成像范圍內(nèi)。

1.2.4指數(shù)計(jì)算與模型構(gòu)建

本文主要選取CAI和DFI反應(yīng)NPV變化特征,CAI計(jì)算公式為[15]：

CAI=[0.5×(R2.0+R2.2)-R2.1]×100

(1)

式中,對(duì)于ASD高光譜數(shù)據(jù),R2.0為2000—2050 nm波段平均反射率,R2.1為 2080—2130 nm波段平均反射率,R2.2為2190—2240 nm波段平均反射率。為了擴(kuò)大CAI值差異性,乘以系數(shù)100。

DFI計(jì)算公式為[18]：

(2)

式中,B1,B2,B6,B7分別代表MODIS第1、2、6和7波段。為了擴(kuò)大DFI值差異性,乘以系數(shù)100。

依據(jù)不同生活型植被構(gòu)建7種模型,分別為草本、半灌木、灌木、草本+半灌木+灌木、草本+半灌木、草本+灌木和半灌木+灌木。

1.2.5光譜特征

在條件控制上,選擇fNPV基本相同時(shí)草本、半灌木、灌木和總體(fNPV分別為61.24%、61.52%、60.84%和61.08%,最大值和最小值相差僅為0.68%)的反射光譜。去除受大氣影響嚴(yán)重波段和水汽吸收波段,保留350—1300、1450—1750和 2000—2350 nm波段范圍[31],呈現(xiàn)光譜特征曲線(圖2)。

圖2 不同材料關(guān)鍵波段NPV反射光譜曲線Fig.2 The spectral reflectance curves of NPV in key band 灰色柱狀為MODIS波段波長(zhǎng)范圍,紅色虛線方框?yàn)镃AI計(jì)算所需波長(zhǎng)范圍

1.2.6fNPV,CAI 和DFI關(guān)系分析

fNPV和CAI、CAI和DFI、DFI和fNPV關(guān)系分析采用線性擬合分析和P值檢驗(yàn)法。其中,主要采用決定系數(shù)(R2)來(lái)評(píng)估DFI和fNPV擬合程度。顯著性檢查閾值采用P= 0.05和P= 0.001,當(dāng)P>0.05時(shí)不顯著,0.001

1.2.7應(yīng)用案例

實(shí)地考察顯示不同草原區(qū)植被組成結(jié)構(gòu)不同,不同草原類(lèi)型區(qū)fNPV估算采用各區(qū)域優(yōu)勢(shì)植被類(lèi)型所建立模型,故在草甸草原區(qū)選用草本模型,典型草原區(qū)選用草本和半灌木混合模型,荒漠草原區(qū)選用草本、半灌木和灌木混合模型,沙地草場(chǎng)選用草本和灌木混合模型進(jìn)行區(qū)域fNPV反演。其中,裸地蓋度值(fBS)采用樣方遮蔭拍照后解譯計(jì)算,使用遙感數(shù)據(jù)計(jì)算得到的fNPV值減去fBS值即為實(shí)際fNPV值。

2 結(jié)果

2.1 不同NPV反射光譜特征

在蓋度相同情況下(圖2),草本、半灌木、灌木和總體光譜反射率具有相似變化趨勢(shì),但是光譜反射率具有差別。350—1300 nm和1450—1750 nm波長(zhǎng)范圍,草本和灌木NPV光譜反射率大于總體,半灌木小于總體。2000—2300 nm波長(zhǎng)范圍,不同生活型NPV光譜反射率差異較小,但草本NPV光譜反射率小于總體,半灌木大于總體。

從不同波段來(lái)看(圖2),B1波段光譜反射率由高到低為草本、灌木、總體和半灌木(20.75%、18.93%、18.49%和15.78%),灌木和總體光譜反射率差值最小(0.44%)；B2波段由高到低為草本、灌木、總體和半灌木(28.36%、28.33%、26.2%和21.93%),總體和半灌木光譜反射率差值最大(4.27%)；B6波段由高到低為草本、灌木、總體和半灌木(39.87%、39.58%、38.9%和37.26%),灌木和總體光譜反射率差值最小(0.58%),總體和半灌木光譜反射率差值最大(1.64%)；B7波段由高到低為半灌木、灌木、總體和草本(36.06%、35.69%、35.57%和34.95%),灌木和總體光譜反射率差值最小(0.12%),總體和草本光譜反射率差值最大(0.62%)。其中,B1和B6波段草本和半灌木光譜反射率都為最高和最低,B2波段半灌木光譜反射率最低,；B7波段草本和半灌木光譜反射率為最低和最高,總體和灌木在2149 nm處產(chǎn)生交叉,波長(zhǎng)為2105—2149 nm時(shí),灌木光譜反射率高于總體,波長(zhǎng)為2149—2155 nm時(shí),灌木光譜反射率低于總體。在B1和B6波段易于分辨草本和半灌木,B2波段最易于區(qū)分半灌木與其他生活型植物,在B7波段易于區(qū)分草本和半灌木以及灌木和總體。

從fNPV,CAI 和DFI分布來(lái)看(圖3),fNPV隨著CAI 和DFI的增加在不斷增加。DFI變化在-9.43—19.48之間,平均值為3.76；CAI變化在-4.5—3.6之間,平均值為-0.23。

圖3 不同生活型植物fNPV、CAI和DFI的3D空間圖 Fig.3 Distribution of fNPV, CAI and DFI in different non-photosynthetic vegetation types

2.2 fNPV,CAI 和DFI關(guān)系

從CAI 和DFI關(guān)系模型來(lái)看(表1),CAI隨著DFI的增加呈增加趨勢(shì),二者呈極顯著正相關(guān)(P<0.001)。其中,半灌木變化率最大,為0.3294；半灌木和灌木混合模型變化率最小,為0.1816。草本、半灌木和灌木混合模型變化率為0.2323,草本和半灌木混合模型變化率為0.2485,草本和灌木混合模型變化率為0.246。

從CAI 和fNPV關(guān)系模型來(lái)看(表1),fNPV隨著CAI的增加呈增加的趨勢(shì),二者呈極顯著正相關(guān)(P<0.001)。其中,半灌木變化率最大,為21.03；草本模型變化率最小,為12.82。草本、半灌木和灌木混合模型變化率為13.64,草本和半灌木混合模型變化率為13.38,草本和灌木混合模型變化率為13.14。

表1 不同生活型植物與組成的fNPV和CAI、CAI和DFI的關(guān)系

從fNPV1和DFI關(guān)系來(lái)看(表2),變化率最大和最小的是半灌木和灌木,分別為6.927和2.9678；草本+半灌木+灌木、草本+半灌木和草本+灌木區(qū)域的變化率分別為3.1688、3.3237和3.2332。其中,fBS值最大和最小的為半灌木+灌木和草本,分別為28.1001和8.0507；草本+半灌木+灌木、草本+半灌木和草本+灌木區(qū)域的fBS值分別為21.0085、17.2148和17.8186。

表2 不同生活型植物與組成的fNPV和DFI的關(guān)系以及各模型中fBS值

從fNPV2和DFI關(guān)系來(lái)看(表2),變化率最大和最小的是半灌木和草本+半灌木,分別為7.3905和3.3245；草本+半灌木+灌木、草本+灌木和草本區(qū)域的變化率分別為3.453、3.9402和3.8019。其中,fBS值最大和最小的是半灌木+灌木和草本,為22.4315和5.9261；草本+半灌木+灌木、草本+半灌木和草本+灌木模型的fBS值分別為18.651、17.1992和14.4532。

2.3 不同草原區(qū)案例分析

從包含裸地的不同草原區(qū)平均fNPV來(lái)看(表3,圖4),fNPV1在4種草原區(qū)平均值為67.33%,fNPV2在4種草原區(qū)平均值為68.47%,fNPV1比f(wàn)NPV2小1.14%。其中,草甸草原區(qū)fNPV最大,fNPV1為84%,fNPV2為86.08%,二者相差2.08%；荒漠草原區(qū)fNPV最小,fNPV1為53.78%,fNPV2為54.36%,二者相差0.58%。典型草原區(qū)fNPV1和fNPV2差值最小,為0.01%,沙地草場(chǎng)fNPV1和fNPV2差值為1.91%。

從剔除裸地后的不同草原區(qū)平均fNPV來(lái)看(表3,圖4),fNPV1在4種草原區(qū)平均值為51.3%,fNPV2在4種草原區(qū)平均值為54.42%,fNPV1比f(wàn)NPV2小3.11%。其中,草甸草原區(qū)fNPV最大,fNPV1為75.95%,fNPV2為80.15%,二者相差4.2%；荒漠草原區(qū)fNPV最小,fNPV1為32.77%,fNPV2為35.71%,二者相差2.94%。典型草原區(qū)fNPV1和fNPV2差值最小,為0.03%,沙地草場(chǎng)fNPV1和fNPV2差值最大,為5.28%。

表3 不同草原區(qū)平均fNPV

3 討論

fPV和fNPV能較好的反應(yīng)區(qū)域植被生物量和草原植被生長(zhǎng)狀況,尤其是fNPV估算對(duì)于草原防火和風(fēng)蝕估算具有重要意義,被廣泛應(yīng)用到風(fēng)蝕計(jì)算中。fNPV受空間異質(zhì)性、季節(jié)、裸地性質(zhì)以及提取方法和數(shù)據(jù)等眾多因素影響[2, 9, 21, 23, 25—26]。從本文研究來(lái)看,植物類(lèi)型差異、裸地干擾和建立模型方法對(duì)fNPV影響較大。

首先,從NPV光譜特征來(lái)看,不同生活型NPV光譜變化趨勢(shì)與已有研究結(jié)果枯落物變化趨勢(shì)一致(圖2),這說(shuō)明光譜變化趨勢(shì)受植物類(lèi)型影響微弱[18, 21, 23, 26]。但是同樣植被蓋度條件下,NPV在不同波長(zhǎng)范圍表現(xiàn)出一定差異說(shuō)明不同生活型NPV光譜反射率不同(圖2),原因在于不同植物在纖維素吸收方面可能存在差異,尤其是在2100 nm附近纖維素與木質(zhì)素吸收容易受水分影響[25, 28]。利用這四個(gè)波段光譜反射率差異可以區(qū)別草本、半灌木、灌木和總體等不同生活型NPV。

其次, DFI和CAI呈極顯著正相關(guān)(P<0.001)關(guān)系(圖3,表2),說(shuō)明DFI在難以獲取CAI所需的波段時(shí)可以代替CAI進(jìn)行植被蓋度計(jì)算；通過(guò)對(duì)比f(wàn)NPV1和fNPV2,二者在典型草原區(qū)差異最小(0.01%和0.03%)(表4),但在其他區(qū)域差異較大,說(shuō)明直接應(yīng)用DFI計(jì)算fNPV可能受植被類(lèi)型和空間異質(zhì)性等因子影響,也說(shuō)明在典型草原中應(yīng)用DFI直接計(jì)算fNPV受限程度較小。DFI和fNPV存在極好線性擬合關(guān)系(圖3),呈極顯著正相關(guān)關(guān)系(P<0.001),這一結(jié)果與已有研究結(jié)果一致[18, 24, 27]。但從草本+半灌木+灌木混合來(lái)看,本研究DFI和fNPV2決定系數(shù)R2(0.65)略小于Cao等和Chai等決定系數(shù)R2(0.95和0.68)[19, 27]。原因是本研究采樣區(qū)鄂爾多斯沙地與錫林郭勒草原的裸地類(lèi)型和裸地含水量等自然條件存在異質(zhì)性,進(jìn)而導(dǎo)致NPV光譜反射率存在差異[28],具體表現(xiàn)在DFI,本研究DFI(平均為3.87%)小于錫林郭勒草原(最小為3.97%)[21]；二是本研究所選物種包括草本、半灌木和灌木,所選NPV比Cao等所選草本和灌木更加豐富[19]。與Ren等對(duì)內(nèi)蒙古荒漠草原CAI和fNPV決定系數(shù)R2相差較小(R2= 0.67,n= 155)[11],在這種樣本數(shù)量基本相同情況下,出現(xiàn)R2相差較小現(xiàn)象的原因可能是CAI和DFI高度相關(guān)。從具體生活型NPV的DFI與fNPV關(guān)系來(lái)看(表3),半灌木斜率(6.927和7.3905)與草本+半灌木+灌木混合模型斜率(3.1688和3.453)差值最大(3.7582和3.9375),這說(shuō)明不同植被類(lèi)型DFI與fNPV關(guān)系存在較大差異,也說(shuō)明半灌木DFI和fNPV關(guān)系特征最明顯,在遙感監(jiān)測(cè)過(guò)程中易于區(qū)分。對(duì)比f(wàn)NPV空間分布來(lái)看(圖4),NPV空間異質(zhì)性和生活型植物分布比例差異說(shuō)明草原區(qū)不同區(qū)域應(yīng)參考不同擬合模型[18],不可將NPV作為均質(zhì)材料處理。

圖4 不同草原類(lèi)型區(qū)fNPV空間分布Fig.4 Distribution of fNPV in different grassland areas A為fNPV1剔除裸地,B為fNPV2剔除裸地,C為fNPV1包含裸地,D為fNPV2包含裸地；1為草甸草原,2為典型草原,3為荒漠草原,4為沙地草場(chǎng)

最后,通過(guò)對(duì)比剔除裸地前后結(jié)果說(shuō)明在一定程度上裸地反射率對(duì)fNPV估算產(chǎn)生干擾。對(duì)比王光鎮(zhèn)等應(yīng)用NDVI-DFI三分模型所得西烏珠穆沁旗4月fNPV高達(dá)77.44%的結(jié)果[26],本研究草甸草原區(qū)(西烏珠穆沁旗)剔除裸地的fNPV1更準(zhǔn)確(剔除裸地的fNPV1和fNPV2為75.95%和80.15%,與已有結(jié)果僅相差1.49%和2.71%；包含裸地fNPV1和fNPV2為84%和86.08%,與已有結(jié)果相差6.56%和8.64%)。本結(jié)果fNPV1偏小1.49%的原因可能是遙感數(shù)據(jù)選用時(shí)間為4月15日,該時(shí)間比4月5日返青植被多。本結(jié)果高于柴國(guó)奇等應(yīng)用Sentinel- 2數(shù)據(jù)對(duì)錫林郭勒4月初fNPV為56%的結(jié)果[32],導(dǎo)致這種結(jié)果的原因是受草原類(lèi)型影響,本研究選取草甸草原與柴國(guó)奇等選取典型草原枯落物數(shù)量差異導(dǎo)致蓋度差異。柴國(guó)奇[32]等結(jié)果與本研究典型草原結(jié)果(剔除裸地fNPV1和fNPV2為50.61%和50.64%,包含裸地的fNPV1和fNPV2為67.83%和67.84%)相差較小(5.39%和5.36%,11.83%和11.84%),這也證實(shí)了剔除裸地后的fNPV結(jié)果更準(zhǔn)確。已有研究通過(guò)端元值進(jìn)行蓋度提取時(shí)[6—7, 17, 23],也充分考慮了裸地的影響,所以在建立模型應(yīng)用時(shí),裸地是必須考慮的重要影響因子。

4 結(jié)論

本文基于實(shí)測(cè)數(shù)據(jù)和MODIS數(shù)據(jù)對(duì)不同生活型植物NPV光譜反射特征進(jìn)行了分析,構(gòu)建了不同生活型植物DFI、CAI和fNPV模型并對(duì)不同草原區(qū)fNPV進(jìn)行了估算,主要得到以下結(jié)論：

(1)草本、半灌木、灌木和總體NPV光譜反射率具有相似變化趨勢(shì),但反射值存在顯著差異,波長(zhǎng)范圍為350—1300 nm和1450—1750 nm時(shí),草本和灌木NPV光譜反射率大于總體,半灌木小于總體；波長(zhǎng)范圍為2000—2300 nm時(shí),草本NPV光譜反射率小于總體,半灌木大于總體。

(2)DFI和CAI、CAI和fNPV、DFI和fNPV均呈顯著正相關(guān)關(guān)系(P<0.001),DFI可以代替CAI進(jìn)行fNPV的估算。

(3)不同草原區(qū)fNPV存在異質(zhì)性,草甸草原、典型草原、沙地草場(chǎng)和荒漠草原分別為75.95%、50.61%、45.88%和32.77%。DFI在典型草原區(qū)估算fNPV效果最好(fNPV1和fNPV2平均差值為3.11%,最小差值為0.03%)。

(4)在進(jìn)行fNPV估測(cè)時(shí),不同草原區(qū)應(yīng)參考不同關(guān)系模型,不可將NPV作為均質(zhì)材料處理,還要剔除裸地干擾。