張智韜，杜瑞麒，楊 帥，楊 寧，魏廣飛，姚志華，邱元霖

·農(nóng)業(yè)信息與電氣技術(shù)·

水鹽交互作用對(duì)河套灌區(qū)土壤光譜特征的影響

張智韜1,2，杜瑞麒1,2，楊 帥1，楊 寧1，魏廣飛1，姚志華1，邱元霖1

（1. 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，楊凌 712100；2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室，楊凌 712100）

為探究土壤水鹽交互作用對(duì)Sentinel-2衛(wèi)星光譜特征的影響，該研究以內(nèi)蒙古河套灌區(qū)沙壕渠灌域?yàn)檠芯繀^(qū)域，分別在2018年和2019年的4—5月共采集280個(gè)裸土期表層土壤樣本，測(cè)定其土壤含水率和含鹽量，并獲取同步的Sentinel-2衛(wèi)星遙感數(shù)據(jù)，構(gòu)建基于土壤水鹽-反射率原理的土壤光譜特征理論模型，在此基礎(chǔ)上結(jié)合土壤水鹽交互作用構(gòu)建水鹽交互模型，并比較2種模型對(duì)土壤光譜的模擬效果，分析土壤水鹽交互作用對(duì)土壤光譜估算的影響。結(jié)果表明：1）土壤水鹽交互作用對(duì)光譜的影響因波段類型和水鹽含量的不同而有所不同。在可見光范圍上影響相對(duì)較弱，其作用范圍為-0.11～0.29；在近紅外和短波紅外范圍上影響相對(duì)較強(qiáng)，其作用范圍為-0.35～0.61；當(dāng)水分或鹽分中某個(gè)含量較高時(shí)對(duì)光譜影響較弱，在水鹽含量程度相似時(shí)影響較強(qiáng)。2）與土壤光譜特征理論模型相比，水鹽交互模型能明顯地改善土壤光譜的模擬效果，能將模擬相關(guān)系數(shù)由0.14～0.44提升到0.29～0.59，均方根誤差由0.032～0.082降低到0.029～0.068。該研究結(jié)果揭示了鹽分和水分對(duì)光譜特征的干擾過程，為土壤鹽分的估算提供策略與方法，對(duì)實(shí)現(xiàn)區(qū)域尺度上土壤鹽分的精準(zhǔn)監(jiān)測(cè)具有一定的意義。

土壤；水分；鹽分；交互；Sentinel-2衛(wèi)星；河套灌區(qū)；光譜特征

0 引 言

土壤鹽漬化是全世界普遍發(fā)生的現(xiàn)象，在降雨較少且蒸發(fā)量較大的地區(qū)尤為明顯，會(huì)嚴(yán)重遏制當(dāng)?shù)氐霓r(nóng)業(yè)生產(chǎn)發(fā)展以及影響其生態(tài)系統(tǒng)[1-4]。衛(wèi)星遙感能快速而又準(zhǔn)確地監(jiān)測(cè)到土壤鹽漬動(dòng)態(tài)變化，為土地治理和生態(tài)恢復(fù)提供有效的信息和管理方法[5]。

利用衛(wèi)星遙感技術(shù)對(duì)土壤鹽分進(jìn)行反演的研究已經(jīng)非常深入[6-10]。Wei等[11]利用多光譜傳感器獲取的遙感影像數(shù)據(jù)，建立基于不同變量組合的土壤鹽分反演模型，并根據(jù)反演結(jié)果篩選出最佳鹽分反演模型；Abderrazak等[12]通過構(gòu)建新的光譜指數(shù)對(duì)摩洛哥中部塔德拉平原土壤鹽漬化進(jìn)行時(shí)空監(jiān)測(cè)；Jiang等[13]利用土壤鹽分在不同深度之間的內(nèi)在聯(lián)系，來監(jiān)測(cè)新疆干旱區(qū)綠洲不同深度土壤鹽分。雖然上述研究已經(jīng)成功地將光學(xué)遙感應(yīng)用于土壤鹽分的反演，但仍存在精度不足的問題。由于水分吸收帶的范圍和鹽分在光譜上有較大的重合，使得當(dāng)研究區(qū)域內(nèi)的土壤水分在時(shí)空上具有很強(qiáng)的動(dòng)態(tài)變化時(shí)，土壤水分和鹽分對(duì)土壤反射光譜有相似的影響，從而導(dǎo)致土壤水分給土壤鹽分反演帶來一定干擾[14]。因此，有必要研究土壤水鹽對(duì)土壤光譜特征的影響過程與作用方式。

針對(duì)此問題，Lobell等[15]根據(jù)土壤在不同水分條件下的土壤光譜特征來構(gòu)建土壤水分與土壤光譜反射率之間的數(shù)學(xué)模型，提出兩者之間存在指數(shù)關(guān)系；Wang等[16]在Lobell等[15]的基礎(chǔ)上通過建立土壤水鹽與光譜特征之間的數(shù)學(xué)模型分析不同水分條件下的土壤鹽分對(duì)土壤光譜特征的作用機(jī)制；Yang等[17]根據(jù)上述研究結(jié)果，利用土壤水鹽-反射率原理構(gòu)建了模擬土壤反射率的數(shù)學(xué)模型，并分析土壤在不同水鹽條件下的光譜特征，從而得出了土壤水分和鹽分是以相互獨(dú)立的方式對(duì)土壤光譜產(chǎn)生影響的結(jié)論，但模擬精度不高，這可能是研究中未考慮水鹽之間的交互作用，使得水鹽與光譜特征之間的影響機(jī)制存在一定的局限性。因此，研究水鹽交互作用對(duì)土壤光譜特征的影響具有重要意義。與此同時(shí)，上述研究涉及的試驗(yàn)均是在室內(nèi)進(jìn)行，通過嚴(yán)格控制條件達(dá)到理想狀態(tài)，并且研究所用到的光譜測(cè)定儀器均為手持設(shè)備，無法像光學(xué)衛(wèi)星一樣來提供大區(qū)域且長(zhǎng)時(shí)間的光譜信息，這些問題使得研究中的發(fā)現(xiàn)與結(jié)論可能無法反映出室外的真實(shí)情況，在野外實(shí)際應(yīng)用上存在一定困難。相比于其他衛(wèi)星[18-22]，于2015年發(fā)射的Sentinel-2衛(wèi)星攜帶著一枚波段范圍在443～2 190 nm的多光譜儀器，可提供周期為5 d、分辨率上至10 m的光譜圖像。優(yōu)秀的光譜能力，較短的回訪周期以及較高的空間分辨率使得Sentinel-2衛(wèi)星成為了目前全球環(huán)境監(jiān)測(cè)的重要組成部分[23]，為研究區(qū)域尺度上的水鹽交互提供了可能與機(jī)遇。目前而言，Sentinel-2衛(wèi)星在水鹽交互上的研究較少。

由于河套灌區(qū)存在降水少蒸發(fā)大的特點(diǎn)，使得區(qū)域內(nèi)土鹽漬化土壤面積已達(dá)到39萬hm2，土壤鹽漬化的現(xiàn)象非常嚴(yán)重。因此，本文以Sentinel-2衛(wèi)星采集的研究區(qū)遙感圖像為數(shù)據(jù)源，河套灌區(qū)不同鹽漬化程度的表層土壤為研究對(duì)象，根據(jù)土壤水鹽-反射率的原理構(gòu)建理論模型，在此基礎(chǔ)上考慮水鹽交互作用構(gòu)建水鹽交互模型，使用2種模型分別模擬土壤在不同水鹽條件下的光譜反射率，通過比較2個(gè)模型的模擬效果來研究水鹽交互作用對(duì)土壤光譜特征的影響，以期為提高反演土壤鹽分精度提供參考。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)域位于河套灌區(qū)解放閘灌域下的沙壕渠（107°06′～107°11′E，40°54′～41°60′N），該灌區(qū)位于內(nèi)蒙古西部，東接包頭市，西到烏蘭布和沙漠，南臨黃河，北至陰山，其灌溉面積有0.688億hm2，是中國(guó)的三大灌區(qū)之一（圖1a）。河套灌區(qū)屬于典型的大陸性氣候，年降雨量為158 mm，年蒸發(fā)量超過2 000 mm，土質(zhì)以粉質(zhì)黏壤土為主。

圖1 研究區(qū)域及采樣點(diǎn)分布

1.2 數(shù)據(jù)來源

1.2.1 Sentinel-2衛(wèi)星數(shù)據(jù)

Sentinel-2號(hào)衛(wèi)星攜帶1枚多光譜成像儀，高度為786 km，可覆蓋13個(gè)光譜波段，幅寬達(dá)290 km。地面分辨率分別為10、20和60 m、每顆衛(wèi)星的重訪周期為10 d，2顆互補(bǔ)，重訪周期為5 d。在USGS網(wǎng)站下載指定日期的Sentinal-2衛(wèi)星影像（2018年4月30日、2018年5月13日、2019年4月16日和2019年5月18日），利用Sen2cor軟件對(duì)衛(wèi)星影像進(jìn)行大氣校正，并在SNAP里針對(duì)特定波段導(dǎo)出成tif格式，最終在ENVI軟件中讀取數(shù)據(jù)。

1.2.2 土壤水分與鹽分?jǐn)?shù)據(jù)

本研究于2018—2019年在河套灌區(qū)沙壕渠進(jìn)行共計(jì)4次土壤采樣工作，采樣時(shí)間分別為2018年4月30日、2018年5月13日、2019年4月16日和2019年5月18日，每次采樣均布置了70個(gè)采樣點(diǎn)，總計(jì)280個(gè)采樣樣本。采樣點(diǎn)涵蓋了不同鹽分梯度水平，采樣點(diǎn)的深度在0～10 cm，考慮到衛(wèi)星分辨率的問題，采樣是由五點(diǎn)法完成，采樣完畢時(shí)記錄相應(yīng)的GPS位置，結(jié)果如圖1b所示。

在野外采集土樣時(shí)，每個(gè)土樣在剔除雜質(zhì)后放入鋁盒內(nèi)，并記錄鋁盒的編號(hào)和質(zhì)量。在土樣帶回室內(nèi)時(shí)，將沒加蓋的鋁盒放入干燥箱內(nèi)，在105 ℃、24 h恒溫條件下用干燥法測(cè)得土壤樣本含水率；烘干之后進(jìn)行研磨，經(jīng)過質(zhì)量比1∶5進(jìn)行配置溶液靜置沉淀后，采用電導(dǎo)率儀（DDS-307A，上海佑科儀器公司）進(jìn)行溶液電導(dǎo)率的測(cè)定（dS/m）。

根據(jù)測(cè)定結(jié)果，可得出采樣點(diǎn)土壤含水率和含鹽量數(shù)據(jù)的統(tǒng)計(jì)分析情況，結(jié)果如表1所示。

表1 土壤鹽分和水分統(tǒng)計(jì)

1.3 研究方法

1.3.1 土壤光譜特征理論模型

為了能夠表示土壤水分和鹽分對(duì)土壤光譜特征的綜合影響，Lobell等[15-17]構(gòu)建出如下所示的量化土壤水分、鹽分和土壤反射率之間定量關(guān)系的模型

式中為受水分和鹽分影響下的土壤表面反射率；R為不受水分和鹽分影響的反射率；SS為土壤電導(dǎo)率（表征鹽分），dS/m；SM為土壤含水率，g/g；表示因土壤水分含量改變而引起土壤表面反射率改變的速率，無量綱；表示因土壤含鹽量改變而引起土壤表面反射率改變的速率，無量綱。R、和會(huì)因?yàn)椴ㄩL(zhǎng)的不同而有所不同。

對(duì)于土壤含水率為SM1、含鹽量為SS1的某土壤樣本來說，其波長(zhǎng)為的理論土壤光譜反射率1可以表示為

同樣地，對(duì)于土壤含水率為SM2、含鹽量為SS2的土壤樣本來說，其波長(zhǎng)為的理論土壤光譜反射率2可表示為

將式（2）和式（3）兩邊同時(shí)進(jìn)行對(duì)數(shù)計(jì)算并相減，可以得到式（4），即

通過式（5）即可求出式（2）中的參數(shù)，即土壤水分含量改變時(shí)所引起土壤反射率改變的速率。

通過式（6）即可求出式（2）中的參數(shù)，即土壤鹽分含量改變時(shí)所引起土壤反射率改變的速率。

綜上，通過所獲取的土壤樣本水分和鹽分?jǐn)?shù)據(jù)，結(jié)合式（5）和式（6）即可計(jì)算出參數(shù)和參數(shù)，使得可通過式（1）構(gòu)建出土壤在不同水鹽條件下的理論光譜特征模型。

1.3.2 土壤水鹽交互作用的定義

Lobell等[15-17]的研究表明，盡管通過式（1）能夠反映土壤水分、鹽分與土壤光譜之間的關(guān)系，但模擬結(jié)果仍存在偏差，這說明該模型在表達(dá)土壤水分和鹽分對(duì)土壤光譜影響上仍具有局限性。從式（1）可看出，水分和鹽分是以乘積的形式來表征對(duì)土壤光譜的作用，并未考慮水鹽之間是否會(huì)以交互作用形式來對(duì)光譜特征產(chǎn)生影響，這可能是造成估算偏差的原因。為此，本文在式（1）的基礎(chǔ)上，結(jié)合水鹽交互作用來量化土壤水分、鹽分與土壤光譜特征之間的關(guān)系，其具體表達(dá)為

通過式（8）得出不同水鹽條件下的水鹽交互作用，即為

1.3.3 用特征量來表達(dá)水鹽交互作用

根據(jù)交互作用在統(tǒng)計(jì)學(xué)上的定義，土壤的水鹽交互作用可解釋為，土壤含水率相同時(shí)，土壤光譜會(huì)因鹽分水平不同而發(fā)生變化，同樣地，土壤鹽分相同時(shí)，土壤光譜會(huì)因含水率水平不同而發(fā)生變化。為進(jìn)一步描述土壤水鹽的交互作用，采用特征量來進(jìn)行量化，其具體如下步驟：

1）繪制土壤水鹽分布的散點(diǎn)圖

為消除土壤水分和鹽分單位差異以及能更好地對(duì)兩者進(jìn)行比較，本文對(duì)土壤水鹽數(shù)據(jù)進(jìn)行均一化計(jì)算，其計(jì)算公式如下

式中SMmax和SSmax分別表示所有土壤樣本中水分含量和含鹽量的最大值，SMmin和SSmin分別表示所有土壤樣本中水分含量和含鹽量的最小值，SMi和SSi分別為第個(gè)土壤樣本的水分含量和鹽分含量。通過式（10）和式（11）可看出，經(jīng)均一化后的土壤水分和鹽分的數(shù)值范圍都在0～1之間，其值越大則代表其水平越高。

2）計(jì)算鹽水比線

3）特征量的計(jì)算

土壤的水鹽交互作用實(shí)際上是水中的氫氧根離子和鹽分中的離子（如碳酸氫根離子，硫酸根離子，氯離子）之間的相互作用過程，但最為直觀的體現(xiàn)則是水分和鹽分含量的相對(duì)比例，故本文選取特征量來量化這一特征，其具體的計(jì)算方式如圖2所示。

圖2 土壤水鹽分布散點(diǎn)圖

4）水鹽交互作用與特征量

1.3.4 水鹽交互模型的構(gòu)建

由式（9）和式（13）可得，在理論模型基礎(chǔ)上，通過結(jié)合土壤水鹽交互作用可構(gòu)建出水鹽交互模型，其具體表達(dá)為

總體上，對(duì)于指定的來說，當(dāng)參數(shù)a、b、Ro以及函數(shù)關(guān)系f已知時(shí)，即可求出土壤在任何水鹽條件下的土壤反射率，為了使文章的思路和計(jì)算過程更加清晰明了，本文繪制了相關(guān)的技術(shù)路線圖，其具體流程如圖3所示。

1.3.5 評(píng)估方法

本文通過相關(guān)系數(shù)（），均方根誤差（Root Mean Square Error，RMSE）以及平均絕對(duì)誤差（Mean Absolute Error，MAE）來評(píng)估模型的模擬效果[11]，其中越大，RMSE和MAE越小則表明模型的模擬效果越好。

2 結(jié)果與分析

2.1 不同水鹽含量下的土壤光譜特征

根據(jù)土壤樣本的含水率和電導(dǎo)率數(shù)據(jù)，篩選出電導(dǎo)率相同但含水率不同的土壤樣本，并將這些樣本對(duì)應(yīng)的Sentinel-2衛(wèi)星的光譜圖像繪制出相應(yīng)的光譜反射率的折線圖，以此分析土壤在同一鹽分含量時(shí)，不同的水分含量與土壤光譜特征之間的關(guān)系，結(jié)果如圖4所示。

圖4 不同電導(dǎo)率下土壤光譜特征隨含水率變化

一般而言，波段反射率與水分含量之間呈負(fù)相關(guān)關(guān)系，即隨著水分含量的增加，波段的反射率會(huì)隨著下降。由圖4可知，當(dāng)電導(dǎo)率EC為30.3 dS/m時(shí)，波段的反射率與水分含量之間的關(guān)系符合這一規(guī)律，但在EC為9.3、21.5和132 dS/m時(shí)，波段的反射率與水分含量之間的關(guān)系不完全符合這一規(guī)律。當(dāng)EC為9.3和21.5 dS/m時(shí)，波段的反射率在含水率0.151~0.222 g/g范圍上會(huì)隨著水分含量的增加而增加；當(dāng)EC為132 dS/m時(shí)，波段的反射率在含水率0.127~0.231 g/g范圍上會(huì)隨著水分含量的增加而降低，符合波段反射率與水分含量之間呈負(fù)相關(guān)關(guān)系的規(guī)律，但通過比較波段反射率在含水率0.118和0.127g/g時(shí)的相對(duì)大小可看出，波段反射率并不會(huì)因?yàn)樗趾康慕档投兴黾?。上述分析表明，波段的反射率與水分含量之間的關(guān)系會(huì)因?yàn)辂}分水平的不同而有所不同，并非一直符合波段反射率與水分含量之間呈負(fù)相關(guān)關(guān)系的規(guī)律。

由圖4還可知，波段反射率與水分含量之間的關(guān)系會(huì)隨著波段的不同而有所不同，波長(zhǎng)在0.496～0.835m范圍上的波段反射率會(huì)隨著波長(zhǎng)的增加而增加；在0.835～2.202m范圍上，波段反射率會(huì)隨著波長(zhǎng)的增加而降低，與此同時(shí)，其增加或降低的趨勢(shì)會(huì)隨著鹽分水平的增加而逐漸放緩，這說明，波段反射率與水分含量之間的關(guān)系在不同波段上的差異性也會(huì)因鹽分水平的不同而有所不同，因此，在不考慮鹽分條件下，通過利用不同波段對(duì)水分敏感性的差異來構(gòu)建能夠反映土壤水分含量的波段差異化指數(shù)，其結(jié)果會(huì)存在偏差。

同樣地，根據(jù)土壤樣本的含水率和電導(dǎo)率的大小進(jìn)行分類，劃分出在滿足含水率相同條件下，具有不同電導(dǎo)率的土壤樣本，并根據(jù)這些樣本對(duì)應(yīng)的Sentinel-2衛(wèi)星的光譜圖像繪制出相應(yīng)的光譜反射率的折線圖，以此分析土壤在同一含水率時(shí)，不同鹽分含量與土壤光譜特征之間的關(guān)系，其結(jié)果如圖5所示。

研究表明[15-17]，波段反射率與鹽分含量之間呈正相關(guān)關(guān)系，即隨著鹽分含量的增加，波段的反射率會(huì)隨著下降。但由圖5可知，土壤樣本所呈現(xiàn)出波段反射率與鹽分含量之間的關(guān)系只在部分鹽分含量范圍內(nèi)符合這一規(guī)律。土壤含水率為0.146 g/g、電導(dǎo)率為29～97 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而增加，而電導(dǎo)率為108～244 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而降低；當(dāng)含水率為0.163 g/g時(shí)，電導(dǎo)率為43～154 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而增加，而電導(dǎo)率為154～211 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而降低；含水率為0.204 g/g、電導(dǎo)率為86～120 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而增加，而電導(dǎo)率為120～276 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而降低；當(dāng)含水率為0.234 g/g時(shí)，電導(dǎo)率為26～96 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而增加，而電導(dǎo)率為96～130 dS/m時(shí)，波段反射率會(huì)因鹽分含量的增加而降低。

通過比較不同水分含量條件下的波段反射率與鹽分含量之間的關(guān)系，即圖5a～圖5d不同的折線圖，可以發(fā)現(xiàn)在水分含量保持不變的條件下，波段的反射率并不會(huì)隨著鹽分含量的改變而出現(xiàn)較大的幅度變化，只有在個(gè)別鹽分含量下，波段的反射率要遠(yuǎn)遠(yuǎn)大于其他鹽分含量下波段的反射率，例如當(dāng)含水率為0.204 g/g時(shí)，0.496～0.945m的波段在電導(dǎo)率為120 dS/m時(shí)的反射率要遠(yuǎn)大于電導(dǎo)率為26、86和276 dS/m下的反射率，而在電導(dǎo)率為26、86和276 dS/m之間的反射率相差并不大，這表明水分會(huì)弱化波段反射率對(duì)鹽分的敏感性，使得波段反射率在不同鹽分條件下不會(huì)相差太大，但對(duì)于波段反射率在個(gè)別鹽分含量上會(huì)發(fā)生較大變化的現(xiàn)象，這或許可以解釋為不同水分條件都會(huì)對(duì)應(yīng)一個(gè)鹽分范圍，在此范圍內(nèi)的波段反射率對(duì)鹽分敏感性要大于水分。

2.2 基于土壤水鹽-反射率原理的理論模型

2.2.1 模型的構(gòu)建

根據(jù)所取土壤樣本的鹽分含量和水分含量數(shù)據(jù)，篩選出鹽分間隔相同下的土壤樣本數(shù)據(jù)集W1-W5，并結(jié)合式（4）進(jìn)行擬合，以此計(jì)算出土壤水鹽-反射率理論模型中的參數(shù)，其擬合的結(jié)果如表1和表2所示。

由表1和表2可知，在多數(shù)波段上，不同鹽分間隔所對(duì)應(yīng)的相關(guān)系數(shù)較為接近，但常數(shù)項(xiàng)有較大差別。從式（4）和式（5）可看出，常數(shù)項(xiàng)能反映出擬合誤差的相對(duì)大小，其值越大，說明擬合效果越不穩(wěn)定，誤差會(huì)越大，因此，在結(jié)合波段在不同鹽分間隔下的相關(guān)系數(shù)和常數(shù)項(xiàng)后，選出不同波段在計(jì)算參數(shù)上所對(duì)應(yīng)的最佳鹽分間隔，并計(jì)算出相應(yīng)的參數(shù)。同理，根據(jù)所取土壤樣本的鹽分含量和水分含量數(shù)據(jù)，篩選出水分間隔相同下的土壤樣本數(shù)據(jù)集S1和S2，并結(jié)合式（5）進(jìn)行擬合，根據(jù)擬合的相關(guān)系數(shù)以及常數(shù)項(xiàng)2，選出不同波段在計(jì)算參數(shù)上的最佳水分間隔，并計(jì)算出相應(yīng)的參數(shù)。根據(jù)上述過程，最終可得到理論模型中的參數(shù)和參數(shù)，其結(jié)果如表3所示。

表1 土壤光譜特征理論模型中參數(shù)a的擬合效果

表2 土壤光譜特征理論模型中擬合的常數(shù)項(xiàng)e1

表3 土壤光譜特征理論模型中參數(shù)a和參數(shù)b

2.2.2 基于理論模型的反射率估計(jì)值

根據(jù)表3的結(jié)果可計(jì)算出波段在不同水鹽水平下的理論反射率，為評(píng)估基于土壤水鹽-反射率原理的理論模型應(yīng)用效果，將各個(gè)波段的理論反射率分別與對(duì)應(yīng)在Sentinal-2衛(wèi)星上的波段反射率建立相關(guān)性分析，并計(jì)算出RMSE和MAE，其結(jié)果如圖6所示。由圖6可知，二者相關(guān)系數(shù)為0.14～0.44，均方根誤差為0.032～0.082。隨著波長(zhǎng)的增加，波段的理論估計(jì)值與Sentinal-2衛(wèi)星上的波段實(shí)測(cè)值之間的相關(guān)系數(shù)會(huì)逐漸減小，RMSE和MAE都有不同程度地上升，這表明隨著波長(zhǎng)的增加，使得通過理論公式來估計(jì)不同水鹽條件下的土壤反射率的效果會(huì)降低。從不同波段所對(duì)應(yīng)的相關(guān)系數(shù)來看，波段2，波段3和波段4的理論估計(jì)值與Sentinal-2衛(wèi)星上的波段實(shí)測(cè)值之間的相關(guān)性較好，說明在Sentinal-2衛(wèi)星上，通過式（1）來估計(jì)不同水鹽水平下的土壤反射率在可見光波段范圍上效果較好，其中在藍(lán)色波段上相關(guān)系數(shù)最高，達(dá)到了0.44。

圖6 基于土壤水鹽-反射率理論估算的反射率與Sentinal-2衛(wèi)星反射率的比較分析

2.3 水鹽交互作用

2.3.1 水鹽交互作用在不同水鹽水平下的特征

在計(jì)算出波段在不同水鹽水平下的理論反射率之后，可根據(jù)式（9）計(jì)算出水鹽交互作用，以此分析每個(gè)波段在不同水鹽條件下水鹽交互作用的特點(diǎn)，其結(jié)果如圖7所示。

圖7 水鹽交互作用在不同水鹽條件下的特征

由圖7可知，水鹽交互作用隨著波段的不同，其對(duì)應(yīng)的區(qū)間范圍和分布特點(diǎn)會(huì)有所不同。從水鹽交互作用區(qū)間上限來看，水鹽作用的區(qū)間上限總體上是隨著波長(zhǎng)的增加而增加，水鹽交互作用的最大值在波段2～波段4上較低（0.15～0.29），在波段8，波段11和波段12上較高（0.43～0.44），在波段5～波段7和波段9上最大（0.50～0.61）；從水鹽交互作用區(qū)間下限來看，除波段9以外，其余各個(gè)波段對(duì)應(yīng)的水鹽交互作用最小值相差不大，主要集中在-0.1左右。從水鹽交互作用的分布特點(diǎn)來看，在波段2～波段5和波段8上，水鹽交互作用呈現(xiàn)出在兩邊大，中間小的分布特點(diǎn)，即水鹽交互作用在水分較低鹽分較高和水分較高鹽分較低的區(qū)域上較大，鹽分水分中間水平的區(qū)域上較小，與此同時(shí)，在波段2～波段5上，隨著波長(zhǎng)的增加，水鹽交互作用較大的范圍會(huì)逐漸增加；在波段6和波段7上，水鹽交互作用分布較為均勻，主要在水分鹽分較高或較低的區(qū)域出現(xiàn)較大的情況；在波段9、波段11和波段12上，水鹽交互作用普遍較大，只有在水分鹽分都較高的區(qū)域上表現(xiàn)較低。由水鹽交互作用的定義可知，當(dāng)波段的水鹽交互作用越大，說明土壤水鹽-反射率理論模型應(yīng)用效果越差，水鹽對(duì)該波段反射率的共同作用程度越高。通過比較不同波段所對(duì)應(yīng)的水鹽交互作用區(qū)間范圍和分布特點(diǎn)，可看出波段類型和水鹽條件是影響水鹽交互作用的2個(gè)重要因素，水鹽交互作用在可見光范圍上整體表現(xiàn)較弱，在近紅外和短波紅外范圍上整體表現(xiàn)較強(qiáng)，但無論在哪個(gè)波段上，水鹽交互作用會(huì)因水鹽含量的不同而有所差別。因此，通過考慮水鹽交互作用來模擬光譜特征會(huì)改善土壤水鹽-反射率理論模型的模擬效果。

2.3.2 水鹽交互作用與特征量的相關(guān)性分析

根據(jù)土壤樣本相應(yīng)的水分均一化，鹽分均一化以及鹽水比，劃分出鹽水比大于1的土壤樣本數(shù)據(jù)，并計(jì)算出相應(yīng)的特征量和不同波段下的水鹽交互作用，將不同波段的水鹽交互作用與特征量之間進(jìn)行相關(guān)性分析，并計(jì)算出擬合誤差RMSE，其結(jié)果如圖8所示。

由圖8可知，從水鹽交互作用與特征量之間的相關(guān)系數(shù)來看，波段2～波段8所顯示出的水鹽交互作用與特征量之間的相關(guān)性要好于或至少等于波段9，波段11和波段12，其相關(guān)性在0.26～0.49，這表明當(dāng)鹽水比大于1時(shí)，即鹽分起主導(dǎo)作用時(shí)，在波段2～波段8上能夠通過鹽水比特征的特征量來較好地表示相應(yīng)波段下的水鹽交互作用。同樣地，針對(duì)鹽水比小于1時(shí)的土壤樣本數(shù)據(jù)計(jì)算出相應(yīng)的特征量和不同波段下的水鹽交互作用，對(duì)每個(gè)波段的水鹽交互作用與特征量之間進(jìn)行相關(guān)性分析，并計(jì)算出擬合誤差RMSE，其結(jié)果如圖9所示。

圖8 鹽水比大于1下的水鹽交互作用（(,)）與特征量（）之間的相關(guān)分析

Fig.8 Correlation analysis between water-salt interaction ((,)) and characteristic quantity () in condition of brine ratio greater than 1

圖9 鹽水比小于1下的水鹽交互作用與特征量之間的相關(guān)分析

由圖9可知，從水鹽交互作用與特征量之間的相關(guān)系數(shù)來看，波段6～波段9和波段11所顯示出的水鹽交互作用與特征量之間的相關(guān)性要好于波段2～波段5和波段12，其相關(guān)系數(shù)在0.41～0.54之間，這表明當(dāng)鹽水比小于1時(shí)，即水分起主導(dǎo)作用時(shí)，在波段6～波段9和波段11上能夠通過具有鹽水比特征的特征量來較好地表示相應(yīng)波段下的水鹽交互作用。

通過比較圖8和圖9可以看出，當(dāng)鹽水比大于1和小于1時(shí)，水鹽交互作用與特征量之間的相關(guān)性較好時(shí)所對(duì)應(yīng)的波段范圍是不同的，且鹽水比小于1時(shí)不同波段的水鹽交互作用與特征量之間的相關(guān)性都要好于鹽水比大于1時(shí)的情況，這表明具有鹽水比特征的特征量與波段的水鹽交互作用之間的相關(guān)性大小會(huì)隨著水分或者鹽分主導(dǎo)地位的改變而有所改變，即水分和鹽分分別起主導(dǎo)作用時(shí)，Sentinal-2衛(wèi)星對(duì)水鹽比例的敏感波段范圍會(huì)有所不同。

2.4 基于土壤水鹽-反射率原理的水鹽交互模型

結(jié)合土壤水鹽-反射率原理和水鹽交互作用構(gòu)建水鹽交互模型，并計(jì)算每個(gè)波段的估計(jì)值，將其與對(duì)應(yīng)在Sentinal-2衛(wèi)星上的波段反射率建立相關(guān)性分析，計(jì)算出RMSE和MAE，以此分析水鹽交互作用對(duì)土壤光譜估算的影響，結(jié)果如圖10所示。

由圖10可知，從波段2～波段12，水鹽交互模型的估計(jì)值與Sentinal-2衛(wèi)星上的波段實(shí)測(cè)值之間的相關(guān)系數(shù)會(huì)逐漸減小，但在部分波段上有所回升，其范圍在0.29～0.59；RMSE逐漸上升，從0.029上升到0.068；MAE在波段2～波段8上波動(dòng)較小，波動(dòng)范圍在11%～14%左右，在波段9，波段11和波段12上波動(dòng)較大，從16.18%上升到21.44%。

圖10 基于水鹽交互模型估算的反射率與Sentinal-2衛(wèi)星反射率的比較

由圖10可知，從水鹽交互模型的估計(jì)值與Sentinal-2衛(wèi)星波段實(shí)測(cè)值之間的相關(guān)系數(shù)來看，波段2和波段3的相關(guān)性最為顯著，相關(guān)系數(shù)分別為0.59和0.54，其次是波段4、波段7和波段8，其相關(guān)系數(shù)在0.43～0.48之間，而波段5、波段6和波段9～波段12所對(duì)應(yīng)的水鹽交互下的估計(jì)值與Sentinal-2衛(wèi)星波段實(shí)測(cè)值之間的相關(guān)性較低，其相關(guān)系數(shù)在0.29～0.38之間；從水鹽交互下的估計(jì)值與Sentinal-2衛(wèi)星波段實(shí)測(cè)值之間的RMSE來看，波段2和波段3的RMSE最小，分別為0.029和0.031，在波段4到波段9上，各自波段的RMSE差別不大，其范圍在0.043～0.053，而波段11和波段12的RMSE較大，分別為0.063和0.068。從水鹽交互模型估計(jì)值與Sentinal-2衛(wèi)星波段實(shí)測(cè)值之間的MAE來看，波段8、波段3和波段4的MAE較小，分別為11.05%、12.26%和12.11%，其次是波段2、波段5、波段6和波段7，MAE范圍為13.03%～14.22%，最后是波段9、波段11和波段12，MAE分別為16.18%、20.28%和21.44%。

通過上述的分析可以得出，將土壤水鹽-反射率理論模型和水鹽交互作用結(jié)合起來計(jì)算出的估計(jì)值和在Sentinal-2衛(wèi)星波段實(shí)測(cè)值之間的關(guān)系在波段2~波段8上要比在波段9、波段11和波段12上更為密切，且擬合效果較好，因此，在Sentinal-2衛(wèi)星波段2~波段8的波長(zhǎng)范圍上，可以通過構(gòu)建考慮水鹽交互作用下的水鹽-反射率模型來估計(jì)土壤在不同水鹽條件及相應(yīng)水鹽交互作用下所產(chǎn)生的土壤反射率。

3 討 論

3.1 水鹽相互作用對(duì)土壤光譜特征的影響

相比于理論模型中的波段估計(jì)值，在土壤水鹽-反射率原理基礎(chǔ)上考慮水鹽交互作用下計(jì)算得出波段估計(jì)值之后，可針對(duì)不同水鹽條件下水鹽的比例特點(diǎn)來進(jìn)行光譜估計(jì)，為了能夠直觀地看到水鹽交互作用對(duì)光譜特征的影響，將理論模型和水鹽交互模型的模擬結(jié)果進(jìn)行比較，其結(jié)果如表4。

表4 水鹽交互作用對(duì)土壤光譜特征的作用

由表4可知，無論是相關(guān)系數(shù)還是均方根誤差，考慮水鹽交互模型的模擬結(jié)果均要優(yōu)于理論模型的模擬結(jié)果，這表明考慮水鹽交互作用能夠提高光譜估計(jì)的精度，但從提高精度的程度來看，均方根誤差RMSE在不同波段上的表現(xiàn)較為穩(wěn)定，而相關(guān)系數(shù)在不同波段上的表現(xiàn)有所不同。在波段2～波段5上，考慮水鹽交互作用的土壤水鹽-反射率模型的相關(guān)系數(shù)較高，相比于土壤水鹽-反射率理論模型，有一定的提升；在波段6～波段8和波段11～波段12上，考慮水鹽交互作用的土壤水鹽-反射率模型的相關(guān)系數(shù)不如在波段2～波段4上高，但相比于土壤水鹽-反射率理論模型的表現(xiàn)，有很明顯的提升，較好地改善土壤水鹽-反射率理論模型在波段6～波段8和波段11～波段12上模擬效果差的問題；在波段9上，水鹽交互模型雖能提高光譜估計(jì)的精度，但模擬效果較低，不如其余波段的模擬效果好。上述分析表明，考慮水鹽交互作用能夠提高波段反射率與土壤水分和鹽分之間關(guān)系的顯著性，尤其是在可見光和近紅外波段上（波段2～波段8），該波長(zhǎng)范圍的反射率能較好地反映出相應(yīng)的土壤水分和鹽分含量，然而，在波段9和波段11~波段12上，該范圍的波段反射率與土壤水分和鹽分關(guān)系較差，對(duì)土壤水分和鹽分含量并不敏感。

由2.3節(jié)可知，Sentinal-2衛(wèi)星對(duì)水鹽比例的敏感波段范圍會(huì)隨著水分或者鹽分主導(dǎo)地位的改變而有所改變，因此將考慮水鹽交互作用和理論的土壤水鹽-反射率模型的模擬結(jié)果在鹽水比大于1和小于1這2種情況下分別進(jìn)行比較，以此分析水鹽交互作用在不同鹽水比例條件下對(duì)光譜特征的影響，其結(jié)果如表5所示。

由表5可知，從相關(guān)系數(shù)和均方根誤差來看，考慮水鹽交互作用的土壤水鹽-反射率模型的模擬結(jié)果在鹽水比大于1和小于1這2種情況下均好于理論的土壤水鹽-反射率模型的模擬結(jié)果，但2個(gè)模型的模擬結(jié)果在波段上的表現(xiàn)會(huì)因?yàn)辂}水比情況的不同而有所差別。當(dāng)鹽水比大于1時(shí)，相比于理論的土壤水鹽-反射率模型，考慮水鹽交互作用的土壤水鹽-反射率模型的光譜估計(jì)精度在波段2～波段5和波段8上提升較大，在其余波段上提升并不明顯；當(dāng)鹽水比小于1時(shí)，理論的土壤水鹽-反射率模型的模擬效果在所有波段上普遍較差，相比于理論的土壤水鹽-反射率模型，考慮水鹽交互作用的土壤水鹽-反射率模型的光譜估計(jì)精度普遍提升，但其模擬效果在波段6～波段9上較好。上述分析表明，無論是鹽分還是水分含量占主導(dǎo)地位，在可見光和近紅外波段上，水鹽交互作用能夠明顯提高波段反射率的估計(jì)精度，但隨著水鹽比例開始增加時(shí)，模型會(huì)在波長(zhǎng)大的波段上取得較好的估計(jì)精度，這表明水分含量比例增加時(shí)，在可見光波段上，波長(zhǎng)較大的波段對(duì)土壤水分和鹽分會(huì)更加敏感。需要注意的是，當(dāng)鹽水比小于1時(shí)，土壤水鹽-反射率理論模型的估計(jì)精度要遠(yuǎn)遠(yuǎn)低于鹽水比大于1時(shí)的估計(jì)精度，說明當(dāng)水鹽比例增加時(shí)，土壤水鹽-反射率理論模型具有一定的局限性，而考慮水鹽交互作用的土壤水鹽-反射率模型能夠克服這一問題，取得較好的估計(jì)精度。

表5 水鹽交互作用在不同水鹽比例下對(duì)土壤光譜特征的影響

3.2 水鹽相互作用的實(shí)際意義

根據(jù)土壤水鹽-反射率的理論估計(jì)值、考慮水鹽交互的估計(jì)值和Sentinal-2衛(wèi)星實(shí)測(cè)值，可以計(jì)算出波段2的理論估計(jì)和考慮水鹽交互的估計(jì)的相對(duì)誤差分布，結(jié)果如圖11所示。

圖11 波段2的理論估計(jì)與考慮水鹽交互估計(jì)的反射率相對(duì)誤差分布

由圖11可知，無論是理論模型還是水鹽交互模型，其模擬結(jié)果都能表明土壤光譜特征會(huì)因水鹽含量的不同而會(huì)有所變化。通過參考光譜估計(jì)的相對(duì)誤差分布圖，能為Sentinal-2衛(wèi)星反演土壤水分或鹽分含量方面上提供參考，并且考慮水鹽交互作用能夠進(jìn)一步提高反演精度，從而解決反演精度較差的問題，增強(qiáng)Sentinal-2衛(wèi)星在監(jiān)測(cè)土壤水分或鹽分變化上的能力，同時(shí)也為其他光譜傳感器在監(jiān)測(cè)土壤水分或鹽分上提供了思路。

本研究根據(jù)土壤水鹽-反射率原理構(gòu)建了模擬土壤光譜特征的理論模型，并在此基礎(chǔ)上定義水鹽交互作用來構(gòu)建水鹽交互模型。相比于理論模型，水鹽交互模型能明顯地改善模擬光譜特征的效果，能更加準(zhǔn)確地描述土壤水鹽對(duì)光譜特征的影響機(jī)制，為從光學(xué)遙感技術(shù)中提高反演土壤鹽分的精度提供了策略與方法，但其研究所用到的土壤樣本均為裸土期，并未涉及到植被覆蓋的情況。當(dāng)作物處于生長(zhǎng)時(shí)期，土壤和植被都會(huì)對(duì)土壤光譜特征產(chǎn)生影響，使得很難將土壤和植被剝離開來單獨(dú)研究土壤水鹽對(duì)光譜特征的影響。不過很多研究結(jié)果表明，通過計(jì)算反映植被狀態(tài)的光譜指數(shù)并構(gòu)建其反演模型是反演土壤鹽分的有效方法，這說明在一定情況下，植被狀態(tài)是影響光譜特征的主要因素。與正常情況相比，當(dāng)作物受到水分脅迫和鹽分脅迫時(shí)，其狀態(tài)的改變會(huì)在光譜特征上有所體現(xiàn)。因此，當(dāng)土壤表面有植被覆蓋時(shí)，可從植被生長(zhǎng)發(fā)育與水鹽的響應(yīng)機(jī)理上來研究水鹽交互作用對(duì)光譜特征影響及模擬，使得在土壤水鹽對(duì)光譜特征的影響機(jī)制上的研究更加深入，通過完善土壤水鹽交互在不同時(shí)期下的表現(xiàn)來構(gòu)建一個(gè)全面的理論體系。

4 結(jié) 論

1）土壤水分和鹽分之間存在交互作用并能以此方式來影響光譜特征。與理論模型相比，通過考慮水鹽交互作用構(gòu)建的土壤水鹽交互模型能明顯地改善模擬效果，相關(guān)系數(shù)由0.14～0.44提升到0.29～0.59，均方根誤差由0.032～0.082降低到0.029～0.068。與此同時(shí)，與理論模型相比，土壤水鹽交互模型的相對(duì)誤差的區(qū)間范圍更小，低誤差所占據(jù)的面積更大，表明考慮水鹽交互作用能夠進(jìn)一步提高反演精度，從而解決反演精度較差的問題，增強(qiáng)Sentinal-2衛(wèi)星在監(jiān)測(cè)土壤水分或鹽分變化上的能力。因此，土壤水鹽交互模型在土壤水分或鹽分的監(jiān)測(cè)上具有更大應(yīng)用潛力。

2）土壤水鹽交互作用影響的強(qiáng)弱會(huì)因波段類型和水鹽含量的不同而有所不同。由圖7可知，在可見光范圍上影響相對(duì)較弱，其作用范圍為-0.11～0.29；在近紅外和短波紅外范圍上影響相對(duì)較強(qiáng)，其作用范圍為-0.35～0.61。

3）基于水鹽含量比例構(gòu)建的特征量能夠反映出土壤的水鹽交互作用，但兩者之間的相關(guān)性會(huì)隨著鹽水比的改變而有所改變。當(dāng)鹽水比大于1時(shí)，水鹽交互作用與特征量之間的相關(guān)性主要集中在0.26～0.49；當(dāng)鹽水比小于1時(shí)，水鹽交互作用與特征量之間的相關(guān)性主要集中在0.35～0.54。

[1]Du Lei, Zheng Zicheng, Li Tingxuan, et al.Effects of irrigation frequency on transportation and accumulation regularity of greenhouse soil salt during different growth stages of pepper[J]. Scientia Horticulturae, 2019, 256: 108568.

[2]Taha Gorji, Aylin Yildirim, Nikou Hamzehpour, et al. Soil salinity analysis of Urmia Lake Basin using Landsat-8 OLI and Sentinel-2A based spectral indices and electrical conductivity measurements[J]. Ecological Indicators, 2020, 112: 106173.

[3]Zheng Wei, Xue Dongmei, Li Xiangzhen, et al. The responses and adaptations of microbial communities to salinity in farmland soils: A molecular ecological network analysis[J]. Applied Soil Ecology, 2017, 120: 239-246.

[4]Oliveira D L G D, Campos C F D, Tales T, et al. Salt-affected soils of the coastal plains in Rio Grande do Sul, Brazil[J]. Geoderma Regional, 2018, 14: e00186-.

[5]黃權(quán)中，徐旭，呂玲嬌，et al. 基于遙感反演河套灌區(qū)土壤鹽分分布及對(duì)作物生長(zhǎng)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(1)：102-109. Huang Quanzhong, Xu Xu, Lǚ Lingjiao, et al. Soil salinity distribution based on remote sensing and its effect on crop growth in Hetao Irrigation District[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(1): 102-109. (in Chinese with English abstract)

[6]Jin Pingbin, Li Pingheng, Wang Quan, et al. Developing and applying novel spectral feature parameters for classifying soil salt types in arid land[J]. Ecological Indicators, 2015, 54: 116-123.

[7]Said Nawar, Henning Buddenbaum, Joachim Hill. Estimation of soil salinity using three quantitative methods based on visible and near-infrared reflectance spectroscopy: A case study from Egypt[J]. Arabian Journal of Geosciences, 2014. doi:10.1007/s12517-014-1580-y.

[8]Peng J, Ji W, Ma Z, et al. Predicting total dissolved salts and soluble ion concentrations in agricultural soils using portable visible near-infrared and mid-infrared spectrometers[J]. Biosystems Engineering, 2015, 152: 94-103.

[9]Rossel R A V, Webster R. Predicting soil properties from the Australian soil visible–near infrared spectroscopic database[J]. European Journal of Soil Science, 2012, 63(6): 848-860.

[10]Yue J, Tian J, Tian Q, et al. Development of soil moisture indices from differences in water absorption between shortwave-infrared bands[J]. Isprs Journal of Photogrammetry & Remote Sensing, 2019, 154: 216-230.

[11]Wei G, Li Y, Zhang Z, et al. Estimation of soil salt content by combining UAV-borne multispectral sensor and machine learning algorithms[J]. PeerJ, 2020, 8(2): e9087.

[12]Abderrazak, El Harti, Lhissou, et al. Spatiotemporal monitoring of soil salinization in irrigated Tadla Plain (Morocco) using satellite spectral indices[J]. International Journal of Applied Earth Observation & Geoinformation, 2016, 50:64-73.

[13]Jiang H, Shu H. Optical remote-sensing data based research on detecting soil salinity at different depth in an arid-area oasis, Xinjiang, China[J/OL]. Earth Science Informatics, 2018. doi:10.1007/s12145-018-0358-2.

[14]Periasamy S, Shanmugam R S. Multispectral and microwave remote sensing models to survey soil moisture and salinity[J]. Land Degradation & Development, 2017, 28(4): 1412-1425.

[15]Lobell D B, Asner G P. Moisture Effects on soil reflectance[J]. Society of America Journal, 2002, 66(3): 722-727.

[16]Wang Q, Li P, Chen X. Modeling salinity effects on soil reflectance under various moisture conditions and its inverse application: A laboratory experiment[J]. Geoderma, 2012, 170: 103-111.

[17]Yang X, Yu Y. Estimating soil salinity under various moisture conditions: an experimental study[J]. IEEE Transactions on Geoscience & Remote Sensing, 2017, 55(5): 2525-2533.

[18]Filion, Rébecca, Bernier M, et al. Remote sensing for mapping soil moisture and drainage potential in semi-arid regions: Applications to the Campidano plain of Sardinia, Italy[J]. Science of the Total Environment, 2016, 543: 862-876..

[19]Wu Weicheng, Ahmad S Mhaimeed, Waleed M Al-Shafie, et al. Mapping soil salinity changes using remote sensing in Central Iraq[J]. Geoderma Regional, 2014, 2/3: 21-31.

[20]Bian J, Zhang Z, Chen J, et al. Simplified evaluation of cotton water stress using high resolution unmanned aerial vehicle thermal imagery[J]. Remote Sensing, 2019, 11(3): 267.

[21]Tanr?kut E, Karaer A, Celik O, et al. Role of endometrial concentrations of heavy metals (cadmium, lead, mercury and arsenic) in the aetiology of unexplained infertility[J]. European Journal of Obstetrics & Gynecology & Reproductive Biology, 2014, 179: 187-190.

[22]Klaus Schneider-Zapp, Manuel Cubero-Castan, Dai Shi, et al. A new method to determine multi-angular reflectance factor from lightweight multispectral cameras with sky sensor in a target-less workflow applicable to UAV[J]. Remote Sensing of Environment, 2019, 229: 60-68.

[23]Reza A, Jalal A, Claudia N, et al. Synergetic use of Sentinel-1 and Sentinel-2 data for soil moisture mapping at plot scale[J]. Remote Sensing, 2018, 10(8):1285.

Effects of water-salt interaction on soil spectral characteristics in Hetao Irrigation Areas of Inner Mongolia, China

Zhang Zhitao1,2, Du Ruiqi1,2, Yang Shuai1, Yang Ning1, Wei Guangfei1, Yao Zhihua1, Qiu Yuanlin1

(1.,,712100,; 2.,,712100,)

Hetao irrigation area is characterized by less precipitation and more evaporation, which makes the sever soil salinization. Satellite based on remote sensing technology can quickly and accurately monitor the dynamic changes of soil salinization, and provide effective information and methods for land management and ecological restoration. Soil moisture and salinity have similar effects on soil spectrum. In order to explore the influence of soil water-salt interaction on the spectral characteristics of Sentinel-2 satellite, a total of 280 surface soil samples were collected from April to May in 2018 and 2019 in Shahao canal irrigation area of Hetao Irrigation Area in Inner Mongolia, China. Soil moisture content and electrical conductivity were measured, and Sentinel-2 satellite remote sensing data were obtained to construct the soil water-salt-reflectance model based on theory of soil spectral characteristics. Combined with soil water-salt interaction, the water-salt interaction model was constructed, and the simulation effects of the two models on soil spectrum were compared, and the influence of soil water-salt interaction on soil spectral estimation was analyzed. The results showed that: 1) the water-salt interaction was affected by soil moisture and salinity. In the visible light range, the effect range was from -0.11 to 0.29. In the near infrared and short wave infrared ranges, the effect was relatively strong with the range from -0.35 to 0.61. When the soil moisture or electrical conductivity were high, the influence was weak. When the soil moisture or electrical conductivity was similar, the effect was strong. 2) Compared with the theoretical model of soil spectral characteristics, the water-salt interaction model significantly improved the simulation accuracy of soil spectrum, increased the simulation correlation coefficient from 0.14-0.44 to 0.29-0.59, and reduced the root mean square error from 0.032-0.082 to 0.029-0.068. At the same time, compared with the theoretical model, the range of relative error of soil water-salt interaction model was smaller, and the area occupied by low error was larger, indicating that the relative error distribution estimated by reference spectrum can provide valuable information for the retrieval of soil moisture or salt content by Sentinal-2 satellite, and the inversion can be further improved by considering water- salt interaction. 3) The characteristic quantities based on the brine ratio could reflect the water-salt interaction of soil, but the correlation between them changed with the change of brine ratio. When the brine ratio was greater than 1, the correlation between water-salt interaction and characteristic quantity was mainly 0.26-0.49; when the brine ratio was less than 1, the correlation between water-salt interaction and characteristic quantity was mainly 0.35-0.54. This study reveals the interference process of salt and water on spectral characteristics, provides strategies and methods for soil salt estimation, and has certain significance for realizing accurate monitoring of soil salinity at regional scale. However, vegetation coverage is not involved in the study. Therefore, it is necessary to study the spectral characteristics of water-salt interaction from vegetation growth and development and water-salt response mechanism in the future.

soils; water; salinity; interaction; sentinel-2 satellite; Hetao Irrigation Area; spectral characteristics

張智韜，杜瑞麒，楊帥，等. 水鹽交互作用對(duì)河套灌區(qū)土壤光譜特征的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(18)：153-164.doi：10.11975/j.issn.1002-6819.2020.18.019 http://www.tcsae.org

Zhang Zhitao, Du Ruiqi, Yang Shuai, et al. Effects of water-salt interaction on soil spectral characteristics in Hetao Irrigation Areas of Inner Mongolia, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(18): 153-164. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.18.019 http://www.tcsae.org

2020-07-03

2019-09-10

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2017YFC0403302、2016YFD0200700）；楊凌示范區(qū)科技計(jì)劃項(xiàng)目（2018GY-03）

張智韜，博士，副教授，主要從事遙感技術(shù)在節(jié)水灌溉及水資源中的應(yīng)用研究。Email：zhitaozhang@126.com

10.11975/j.issn.1002-6819.2020.18.019

X833；S271；S155.1

A

1002-6819(2020)-18-0153-12