徐麗華，謝德體

(西南大學(xué)資源環(huán)境學(xué)院，重慶 400716)

【研究意義】土壤有機(jī)質(zhì) (SOM, Soil Organic Matter) 含量是衡量土壤質(zhì)量和作物長(zhǎng)勢(shì)的關(guān)鍵參數(shù)[1-2]。高光譜技術(shù)由于能快速、無(wú)損地監(jiān)測(cè)SOM的變化，而在SOM的預(yù)測(cè)中得到了廣泛的應(yīng)用[3-4]。【前人研究進(jìn)展】受氣候、母質(zhì)、地形、生物等成土因素、成土年齡以及人類活動(dòng)的影響，不同土壤的理化性質(zhì)有明顯的差別[5]。而土壤的理化性質(zhì)，如SOM、土壤水分含量、氧化鐵含量、機(jī)械組成等因素也是影響土壤光譜特性的重要因素[6]。對(duì)荒漠SOM預(yù)測(cè)的研究顯示，無(wú)人干擾區(qū)和人為干擾區(qū)建立的SOM含量反演模型的精度并不相同[4]。盛浩等人的研究發(fā)現(xiàn)，底土的溶解性有機(jī)質(zhì)的數(shù)量和化學(xué)結(jié)構(gòu)、土壤光譜的曲線形狀及特征值對(duì)土地利用方式的變化比較敏感[7]；但也有研究顯示不同的土壤養(yǎng)分對(duì)土地利用方式變化的敏感度并不相同[8]。【本研究切入點(diǎn)】本研究擬利用來(lái)自園地、旱地、水田3種不同土地利用方式的土壤樣本，嘗試建立不同土地利用方式下，SOM含量高光譜反演模型?！緮M解決的關(guān)鍵問(wèn)題】探討SOM含量反演精度對(duì)不同土地利用類型的響應(yīng)，為分析人為活動(dòng)對(duì)SOM含量的影響提供一定的參考。

圖1 采樣點(diǎn)分布

1 材料與方法

1.1 土壤樣本的野外采集及實(shí)驗(yàn)室化學(xué)分析

在三峽庫(kù)區(qū)王家溝小流域內(nèi)，采集0～20 cm耕層土壤樣本共111個(gè)(圖1)。其中，林地、園地、旱地、水田4種利用方式下的土壤樣本數(shù)分別為6、43、29和33。在實(shí)驗(yàn)室內(nèi)對(duì)土壤樣本進(jìn)行風(fēng)干、研磨、過(guò)篩處理。處理好的樣本被分成2份，分別用于 SOM含量和光譜反射率的測(cè)定。用重鉻酸鉀容量法測(cè)定樣本的SOM含量[9]。

圖2 土壤樣本的反射光譜曲線

1.2 光譜測(cè)量及光譜預(yù)處理

采用FieldSpec Pro FR 野外光譜儀(光譜范圍為350～2500 nm)，對(duì)同一土壤樣本進(jìn)行10次重復(fù)測(cè)量。進(jìn)行樣本有效性的檢查和剔除異常值后取多次測(cè)量的平均值作為該土壤樣本的反射光譜。為了剔除光譜曲線的起始處和末端噪聲的影響，采用400～2450 nm作為光譜分析波段(圖2)。在此基礎(chǔ)上，對(duì)圖2的光譜進(jìn)行了斷點(diǎn)擬合和去包絡(luò)預(yù)處理，以便消除反射光譜的階躍性變化和突出光譜特征(圖3)。

圖3 預(yù)處理后的反射光譜曲線

表1 SOM含量的描述性統(tǒng)計(jì)(g/kg)

1.3 建模方法及模型精度評(píng)價(jià)指標(biāo)

文中采用偏最小二乘回歸(Partial Least Squares Regression, PLSR)[10]法對(duì)SOM含量進(jìn)行建模，并用相關(guān)系數(shù)R對(duì)模型精度進(jìn)行評(píng)價(jià)。相關(guān)系數(shù)一般的判斷標(biāo)準(zhǔn)是[11]： |R|<0.3，則2個(gè)變量算作無(wú)相關(guān)； 0.3≤|R|<0.5，則表示2個(gè)變量之間呈現(xiàn)低度相關(guān)性； 0.5≤|R|<0.8，表示2個(gè)變量具有顯著的相關(guān)性； 0.8≤|R|<1，表示2個(gè)變量具有高度相關(guān)性；R=1，表示2個(gè)變量之間完全線性相關(guān)。

(1)

2 結(jié)果與分析

2.1 不同土地利用方式下的SOM統(tǒng)計(jì)特征

因?yàn)閬?lái)自林地的樣本數(shù)量過(guò)少，所以文中只對(duì)來(lái)自于園地、旱地、水田的樣本進(jìn)行了建模分析。對(duì)來(lái)自園地、旱地、水田的土壤樣本有機(jī)質(zhì)含量進(jìn)行統(tǒng)計(jì)分析，用3倍標(biāo)準(zhǔn)差的檢驗(yàn)原則對(duì)數(shù)據(jù)進(jìn)行離群值檢驗(yàn)。在園地和水田樣本中未發(fā)現(xiàn)離群值，旱地樣本中發(fā)現(xiàn)1個(gè)離群值，在分析時(shí)剔除該離群值。對(duì)不同利用方式下的SOM含量進(jìn)行描述性統(tǒng)計(jì)(表1)。水田和旱地都屬于耕地，受人為干擾的程度和時(shí)間均高于園地，其SOM含量的均值也高于園地；研究區(qū)內(nèi)的水田與旱地均采用輪作的方式，受人為干擾程度基本相同，但水田SOM含量高于旱地，可能是由于水田長(zhǎng)期處于水淹情況下，有機(jī)質(zhì)分解慢的原因。3種土地利用方式下土壤的SOM含量的變異系數(shù)比較接近，但園地土壤的變異系數(shù)低于水田和旱地。

圖4 不同利用類型土壤有機(jī)質(zhì)含量與去包絡(luò)后的光譜反射率的相關(guān)關(guān)系

圖中和分別表示標(biāo)定樣本和驗(yàn)證樣本

2.2 不同土地利用方式的特征波段的選擇

對(duì)經(jīng)過(guò)預(yù)處理的光譜反射率與園地、旱地、水田樣本的SOM含量進(jìn)行相關(guān)分析(圖4)。園地、旱地、水田 SOM含量與光譜反射率的相關(guān)系數(shù)的變化范圍分別是：-0.416～0.465,-0.480～0.490,-0.619～0.575。水田SOM含量與土壤光譜反射率之間相關(guān)系數(shù)最大絕對(duì)值為0.619，二者之間的相關(guān)性較高；與水田相比，旱地SOM含量與光譜反射率的相關(guān)性在整個(gè)波段范圍內(nèi)普遍較低，但旱地SOM含量與光譜反射率的相關(guān)性普遍高于園地。

2.3 不同土地利用方式下的SOM含量預(yù)測(cè)精度

將園地、旱地、水田土壤樣本的2/3的樣本用于標(biāo)定，1/3用于驗(yàn)證，樣本隨機(jī)選取。其中園地、旱地、水田土壤樣本的標(biāo)定集樣本數(shù)分別為29、18和22；園地、旱地、水田土壤樣本的驗(yàn)證集樣本數(shù)分別為14、10和11。采用留一法對(duì)模型進(jìn)行交叉驗(yàn)證，剔除異常值。將剔除異常值后的樣本數(shù)據(jù)用于模型的定標(biāo)和驗(yàn)證。

將所選特征波段的反射率作為輸入變量，分別對(duì)園地、旱地、水田3種利用方式下的SOM含量進(jìn)行PLS建模，并利用獲得的模型預(yù)測(cè)值與SOM含量實(shí)測(cè)值繪制散點(diǎn)圖(圖5)。

從圖5可以看出，對(duì)于園地、旱地和水田土壤有機(jī)質(zhì)含量的反演，無(wú)論是標(biāo)定集還是驗(yàn)證集，其預(yù)測(cè)值和實(shí)測(cè)值的相關(guān)系數(shù)都位于0.5～0.8 (其值分別為0.730和0.626、0.771和0.769、0.865和0.776)，這說(shuō)明不同利用方式下SOM含量反演模型獲得的預(yù)測(cè)值和實(shí)測(cè)值之間均具有顯著的相關(guān)性,但三種不同土地利用下的SOM含量反演精度并不相同。其中，無(wú)論是標(biāo)定集還是驗(yàn)證集，水田SOM含量反演模型獲得的預(yù)測(cè)值和實(shí)測(cè)值的相關(guān)系數(shù)都是最高的，其全體水田樣本預(yù)測(cè)值和實(shí)測(cè)值之間的總相關(guān)系數(shù)為0.838；園地土壤預(yù)測(cè)值和實(shí)測(cè)值的相關(guān)系數(shù)均為最低(總體相關(guān)系數(shù)為0.700)。旱地SOM含量預(yù)測(cè)值和實(shí)測(cè)值的總體相關(guān)系數(shù)為0.781，其反演精度居中。這說(shuō)明，本研究中土地利用方式對(duì)SOM含量反演影響較大。水田和旱地土壤耕作時(shí)間和耕作強(qiáng)度高于園地土壤，其SOM含量反演精度也高于園地SOM含量的反演精度。主要原因可能是由于土地利用方式的不同，使得人類對(duì)不同土壤的干擾時(shí)間和干擾強(qiáng)度有所差異，進(jìn)而造成了不同利用下的土壤光譜和SOM含量的差異，從而影響SOM含量的反演精度。從2.1中可以看出，園地、旱地、水田3種利用方式下SOM含量的變異系數(shù)從大到小分別是水田、旱地、園地。這個(gè)變化趨勢(shì)與SOM含量預(yù)測(cè)精度的排列順序一致，這說(shuō)明SOM含量的反演精度與土壤有機(jī)含量的變異系數(shù)是有一定關(guān)系的。在一定程度上，如果SOM含量分布過(guò)于集中，往往可能難以獲得較高的反演精度。

3 討 論

本研究中，園地和旱地土壤的類型均為紫色土，其反演模型的精度均低于水田SOM反演模型的精度(水田土壤類型為水稻土)。說(shuō)明雖然處于相同的氣候條件下，不同土壤類型的SOM含量預(yù)測(cè)精度是有所差異的。許多學(xué)者在研究不同土壤類型的SOM含量高光譜反演中也得到了相同的結(jié)論。如楊揚(yáng)在用高山草甸土和高山草原土的原始反射光譜分別建立偏最小二乘回歸模型時(shí)發(fā)現(xiàn)，高山草原土的反演模型精度高于高山草甸土[12]。盧艷麗在對(duì)東北平原8種土壤類型的有機(jī)質(zhì)含量進(jìn)行高光譜反演中也發(fā)現(xiàn)，不同土壤類型的反演模型的精度并不相同[13]。因此，土壤類型是影響SOM含量反演精度的一個(gè)重要因素。

但即使是來(lái)自于同一土壤類型的樣本，由于土地利用方式不同，其SOM含量模型的反演精度也不一定相同。研究中，園地和旱地具有相同的土壤類型，由于利用方式不同，其土壤樣本SOM含量與光譜反射率的相關(guān)性不僅存在差異，其SOM含量的反演精度也并不相同。為了驗(yàn)證SOM含量反演模型對(duì)土地利用方式的敏感性，將園地與旱地的土壤樣本合并，選擇4個(gè)特征波段(λ=683,1146,1151,1152 nm)建立紫色土SOM含量反演模型，其標(biāo)定集和驗(yàn)證集獲得的實(shí)測(cè)值和預(yù)測(cè)值的相關(guān)系數(shù)分別為0.501,0.574，總相關(guān)系數(shù)為0.528。其精度均低于按園地和旱地利用類型建立的SOM含量反演模型精度。因此，較單一的土地利用方式而言，混合不同利用方式的樣本建立的SOM含量反演模型的精度較低；與土壤類型相比，SOM含量的反演精度對(duì)土地利用方式更為敏感。其主要原因可能是由于對(duì)不同利用方式的土地有不同的施肥量和耕作措施，會(huì)在一定程度上改變土壤的理化性質(zhì)[2]，使得SOM含量發(fā)生改變，所以土壤光譜反射特征對(duì)土地利用方式的變化會(huì)產(chǎn)生敏感的響應(yīng)[7]。盛浩等人的研究表明，當(dāng)天然林地轉(zhuǎn)變?yōu)槠渌恋乩梅绞胶?，?huì)帶來(lái)底土溶解性有機(jī)質(zhì)的損失和光譜曲線的變化[7]。這些研究也進(jìn)一步說(shuō)明土地利用方式是影響SOM反演模型精度的一個(gè)重要因素。薛莉紅對(duì)稻田、桃園、菜地3種不同土地利用方式下SOM含量的高光譜反演研究中卻得到了不同的結(jié)論[8]。對(duì)于不同供試樣本得到不同的結(jié)論，其原因可能與不同土地利用類型下的施肥種類(有機(jī)肥還是無(wú)機(jī)肥)、施肥量、受人為活動(dòng)擾動(dòng)的程度有關(guān)。

4 結(jié) 論

園地、旱地、水田3種土地利用方式下的SOM含量反演模型具有不同的精度，其中水田SOM含量反演模型的精度最高，園地SOM含量反演精度最低，旱地SOM含量反演精度居中；即使是同一土壤類型，但是由于土地利用方式不同，其SOM含量預(yù)測(cè)模型精度也并不相同， 旱地和園地的土壤類型均為紫色土，但旱地SOM反演模型精度高于園地；較單一的土地利用方式而言，混合不同土地利用方式的樣本建立的SOM含量預(yù)測(cè)模型的精度較低；與土壤類型相比，SOM含量的預(yù)測(cè)精度對(duì)土地利用方式更為敏感。今后的研究可進(jìn)一步從施肥類型及種植年限等方面考慮人為活動(dòng)對(duì)SOM含量反演精度的影響。