焦崗湖沉水植被對水體反射光譜的影響

王杰，崔玉環，李健，王旦，王晟

(1.安徽大學 資源與環境工程學院，合肥 230601；2.安徽大學 濕地生態保護與修復安徽省重點實驗室，合肥 230601；3.安徽農業大學 理學院，合肥 230036)

0 引言

近年，隨著內陸湖泊水環境問題的突出，利用衛星遙感技術對湖泊水質狀況開展實時監測的需求越來越迫切。水質遙感的核心問題就是建立水體表觀反射率與水質參數間的關系，而水體表觀反射率往往會受到水體生物、大氣散射及水氣界面折射等方面的干擾，進而降低了水體內部物質濃度的反演精度[1]。在內陸湖泊中存在著大量的光學淺水區，湖泊底質對水中太陽光線存在著一定的吸收或反射作用，對湖泊離水輻射信號產生了重要影響。

常見的湖泊底質類型有草型、淤泥型、巖石型以及混合型等，存在著不同的光學特性，導致光學淺水區水底反射光對湖泊離水輻射信號產生的影響存在明顯差異[2-4]。與其他底質類型不同，草型底質屬于波段選擇性吸收與散射型，反射特性與波長的關系較為密切。已有研究表明，淺水湖泊沉水植被冠層對水體反射光譜的影響程度因冠層水下深度(簡稱冠層水深)、植被蓋度、水體透明度及遙感波段而異[5-6]。

目前許多研究者利用遙感技術對Laurentian湖、太湖等淺水湖泊的沉水植被生長特征進行識別與監測[7-10]，并探究了水環境因素(如懸浮物、葉綠素a)對沉水植被冠層光譜特征的影響[11-12]。在太湖地區的研究發現，受到沉水植被生長高度、水質條件等因素的影響，草型湖泊底質反射對水面遙感反射率的影響程度存在較大的空間差異[13]。在海南新村港海草底質光學淺水區的光譜測試分析表明，水體反射光譜在550 nm處形成反射峰，隨海水深度增加發生藍移，隨葉面積指數增加發生紅移[14]。然而，針對內陸湖泊光學淺水區中沉水植被生長參數(冠層蓋度、冠層水下深度)對水體表觀反射光譜的影響機理還需要開展進一步的研究。

鑒于此，本文選擇淮河流域焦崗湖為研究區，通過湖泊水體反射光譜、沉水植被生長參數的現場測定，結合水色要素的室內測定結果，探究湖底沉水植被生長狀態對焦崗湖光學淺水區水體反射光譜的影響。開展該研究對提高內陸湖泊水質參數遙感反演精度、改進現有的沉水植被遙感信息提取算法具有重要理論和現實意義。

1 材料與方法

1.1 研究區概況

焦崗湖地跨安徽省淮南市鳳臺縣與阜陽市潁上縣，位于淮河干流北岸，平均水深為0.8 m，屬于淮河流域的淺水型淡水湖泊，為國家級濕地公園、國家4A級旅游景區。該湖泊近似橢圓形，在平水期湖面面積為31.3 km2，東西長度為15 km，南北最寬距離為5 km。湖泊匯水總面積480 km2，納潁上境的古沙河、老墩溝，鳳臺境的濁溝、花水澗來水，經湖區調蓄后，由便民溝經焦崗閘注入淮河。

湖中西部為養殖區，面積近13.0 km2，水體較為渾濁，懸浮物含量在80 mg/L左右，沉水植被生長較少；湖泊東部為非養殖區，水體較清澈，湖泊底質為淤泥，湖底沉水植被生長狀況良好，優勢物種為菹草，伴生物種為狐尾藻、金魚藻、荇菜等，其蓋度在冬季可達到80%～97%[15]。

1.2 水樣采集與測試

本研究采樣工作在2017年12月23—24日開展，天氣晴朗無云，水面平靜。由于焦崗湖沉水植被生長在湖泊非養殖區，光譜采樣點主要分布在湖泊東南部菹草生長區，按照相近蓋度不同冠層水深、相近冠層水深不同蓋度的采樣方案，共布設了10個沉水植被生長區水體光譜樣點(L2～L11)。為與其他類型水體作比較，本文在淤泥底質型光學淺水區、光學深水區共布設了4個樣點(L1、L12、L13、L14)，具體測定參數見表1。

表1 焦崗湖水體采樣點實測參數

2 結果分析

2.1 無植被生長區水體光譜特征

與沉水植被生長區水體反射光譜作比較，本文在湖泊無植被生長區選擇4個樣點，實地測定了無植被生長區水體反射光譜。其中，L1、L12、L14為光學深水區樣點，L13為光學淺水區樣點。

由圖1可看出，焦崗湖光學深水區水體呈現出典型的水體光譜特征，在可見光部分的綠波段0.58 um附近形成了明顯的反射峰，在0.67 um處形成吸收谷，而在近紅外波段反射率明顯下降，在0.7 um、0.8 um附近形成了2個小的反射峰。由于焦崗湖水體中葉綠素含量介于14～26 ug/L，空間差異不明顯(表1)，懸浮物含量成為影響光學深水區水體反射光譜的主要因素。即隨著懸浮物含量的增加，在0.4～0.58 um水體光譜反射率呈減小的趨勢，這與葉綠素對其吸收作用的影響有關，而在0.58～0.9 um光譜反射率呈增大的趨勢，反映了近紅外波段對懸浮物含量變化更為敏感。

圖1 焦崗湖光學深水區水體反射光譜

本文測定了焦崗湖光學淺水區(L13)的水體表面及對應淤泥底質的反射光譜(圖2)。受淤泥底質對入射太陽光線吸收作用的影響，無植被生長的光學淺水區水體反射光譜表現出不同于光學深水區的光譜特征，在可見光部分的紅波段(0.58～0.71 um)有寬的反射峰，在近紅外波段反射率下降，且在0.8 um處有小的反射峰。其中，焦崗湖光學淺水區在0.4～0.61 um光譜反射率要低于光學深水，而在0.61～0.9 um要高于光學深水，反射率最大差值達到0.04，出現在0.8 um附近。

圖2 無植被生長光學淺水區水體反射光譜

由圖中看出，無植被生長光學淺水區光譜反射率介于光學深水與淤泥底質反射光譜之間，這反映了湖泊淤泥底質對水體反射光譜的影響在不同的太陽光譜波段存在一定差異，在藍綠波段以吸收作用為主，在紅光、近紅外波段以反射作用為主。

2.2 不同冠層蓋度條件下的水體光譜特征

除了淤泥型底質，草型底質對湖泊光學淺水區水體反射光譜也有著一定的影響。本文根據焦崗湖菹草生長區的水體反射光譜(樣點L2、L5、L6、L9、L10)，分析在同一冠層水深(20 cm)不同冠層蓋度條件下的水體反射光譜特征。

圖3顯示了在同一冠層水深下不同蓋度沉水植被生長的焦崗湖水體反射光譜曲線。從圖中可以看出，在焦崗湖地區沉水植被生長的水體反射光譜特征與光學深水區(L1)光譜特征相似，在綠光波段0.58 um附近形成了明顯的反射峰，在近紅外波段反射率明顯下降，在0.7 um、0.8 um附近形成了2個小的反射峰。然而，由于沉水植被冠層的選擇性吸收作用[2]，可見光波段的反射率低于光學深水區光譜反射率，且隨著沉水植被蓋度的增大，水體光譜反射率整體呈增大的趨勢，其中在0.45～0.7 um、0.8～0.9 um水體反射光譜對沉水植被蓋度變化最為敏感。

圖3 同一冠層水深不同植被蓋度下水體反射光譜

2.3 不同冠層水深條件下的水體光譜特征

已有研究表明，沉水植被冠層水深也是影響湖泊光學淺水區水體反射光譜的重要因素[12，14，18]。本文測定了相同植被蓋度(80%)不同冠層水深的水體反射光譜(樣點L3、L4、L7、L8、L11)，分析冠層水深變化對光學淺水區水體反射光譜的影響。

由圖4可以看出，在蓋度相同的情況下，隨著植被冠層水深的減少，沉水植被冠層逐漸接近水面，水體光譜反射率在可見光波段(0.4～0.7 um)呈顯著增大的趨勢，而在近紅外波段(0.8～0.9 um)變化趨勢不明顯。當沉水植被冠層生長到水面時，水體反射光譜呈現出近似于陸生植被的反射光譜特征，反映了沉水植被冠層水深對焦崗湖水體反射光譜具有一定的貢獻度，且其貢獻度隨著植被冠層水深的減少而增大。這與馬榮華等[11]在太湖的東部湖灣區利用Hydrolight水體輻射傳輸模型模擬的研究結論相一致。

圖4 相同蓋度不同冠層水深下水體反射光譜

3 討論

3.1 懸浮物濃度及透明度對光學淺水區水體反射光譜的影響

由于研究區水體葉綠素含量較低(表1)，懸浮物含量成為影響湖泊水體反射光譜的重要因素。本文通過分析光學淺水區水體懸浮物含量以及透明度與實測光譜反射率的關系，探究懸浮物含量、透明度對光學淺水區水體反射光譜的影響(圖5)。

圖5 水體懸浮物含量、透明度與水體光譜反射率的相關關系

由圖5看出，在沉水植被生長的湖泊光學淺水區，懸浮物含量與水體光譜反射率總體呈正相關，但其相關性并不顯著，相關系數最高僅達到0.4。這可能是由于沉水植被生長對水體懸浮物顆粒有著強烈吸附作用，導致懸浮物含量較低，相對于沉水植被冠層，懸浮物顆粒的散射作用對光學淺水區水體反射光譜特征的影響并不明顯。

水體透明度與光譜反射率呈負相關關系，且在0.4～0.68 um二者間相關性顯著，相關系數最大達到0.78左右。這是由于透明度能綜合反映水體中懸浮物、葉綠素、CDOM等參數對水體光學性質的影響[5]，可以作為水體反射光譜的指示性指標。

3.2 底質類型對光學淺水區水體反射光譜的影響

在內陸湖泊光學淺水區，湖泊底質是影響水體反射光譜的重要因素，其反射特性對于水底反射光對水面遙感信號的貢獻度有著顯著的影響[19]。本文所研究的焦崗湖地區，湖泊底質包括淤泥型與草型這2種類型。為探究草型底質對水體反射反射光譜的影響，本研究根據湖泊水體樣點L12(光學深水)、L13(淤泥底質光學淺水)、L4(草型底質光學淺水)光譜測定結果，比較并分析不同湖泊底質及光學深水區水體反射光譜特征(圖6)。

圖6 不同湖泊底質類型水體反射光譜

圖6顯示了焦崗湖不同底質類型的光學淺水區以及光學深水區水體反射光譜特征。從圖中可看出，淤泥型底質、草型底質對水體反射光譜的影響存在著明顯差異。在可見光波段0.4～0.6 um，淤泥型底質、草型底質水體反射率均小于湖泊光學深水區，這說明2種類型底質在可見光波段對水體入射光線具有吸收作用，但草型底質的吸收作用明顯強于淤泥型底質；在0.6～0.9 um波譜范圍內，湖泊光學深水區水體反射率高于草型底質水體，低于淤泥型底質水體，最大反射率差值位于0.7 um附近，這反映了草型底質對水體入射光線的選擇性吸收作用，而淤泥型底質則具有較強的反射作用。以上結論在Mobley、周冠華等的研究[2-3]中也得以證實。

3.3 沉水植被生長對光學淺水反射光譜的影響

為進一步分析沉水植被生長對光學淺水區水體反射光譜的影響，本研究分析沉水植被蓋度、冠層水深與水體光譜反射率的關系，探究沉水植被蓋度、冠層水深對水體反射光譜的影響(圖7)。

圖7 沉水植被蓋度、冠層水深與水體光譜反射率的相關關系

從圖7看出，在草型底質湖泊光學淺水區，沉水植被蓋度、冠層水深均對水體光譜反射率存在一定的影響。其中，在相同冠層水深條件下，沉水植被蓋度與水體光譜反射率呈正相關，在0.67 um附近相關系數達到0.81；在相同植被蓋度情況下，隨著冠層水深的減小，水體光譜反射率呈增加的趨勢，在0.7～0.9 um處相關系數較高，最大值達到0.78左右。這反映了焦崗湖沉水植被冠層對入射太陽光具有一定的透射與折射作用[20]，即隨著植被蓋度增大或者冠層水深的減少，沉水植被冠層對太陽光的透射與折射作用加強，其對水體反射光譜的貢獻度增加，進而導致水體光譜反射率的增大。

3.4 構建光學淺水區水體光譜綜合模型

在沉水植被生長的湖泊光學淺水區，水體反射光譜受水體懸浮物含量、葉綠素含量、沉水植被生長參數、底質類型等因素綜合作用的影響。本文基于電磁波與水體相互作用的輻射傳輸理論，在Lee等[21]提出的生物光學模型基礎上，構建了適用于光學淺水沉水植被生長區的水體反射光譜綜合模型。該模型將水體離水輻射信號分解為來自水體自身散射、底質反射兩部分(式(1))。

(1)

(2)

(3)

ρ(λ)=F(T，FC，R，λ)

(4)

4 結束語

本文根據2017年冬季淮河流域焦崗湖光學淺水區水體反射光譜與沉水植被生長參數的實測數據，結合采樣點水質參數實驗室測定，探究沉水植被生長參數對淺水型湖泊光學淺水區水體反射光譜特征的影響，初步得出以下結論：

①在無沉水植被生長區，焦崗湖光學淺水區水體反射光譜介于淤泥底質與光學深水之間，在可見光波段光譜反射率要低于光學深水，而在近紅外波段要高于光學深水，反射率最大差值出現在0.8 um附近，反映了湖泊淤泥底質對水體光譜反射率的影響在不同光譜通道存在明顯差異。

②不同湖泊底質類型對水體反射光譜的影響存在明顯差異。在可見光波段，草型與淤泥型底質對水體入射光線都具有吸收作用，草型底質的吸收作用明顯強于淤泥型底質；在近紅外波段，草型底質、淤泥型底質對水體入射光線具有一定的反射作用，淤泥型底質的反射作用明顯強于草型底質。

③在沉水植被生長區，焦崗湖光學淺水區水體反射光譜特征與光學深水區相似，但整體光譜反射率低于光學深水區，沉水植被蓋度、冠層水深均對水體反射光譜存在影響。植被蓋度與水體光譜反射率呈正相關，在0.67 um附近相關性達到最高；冠層水深與水體光譜反射率呈負相關，在0.88 um處相關性最高。這反映了沉水植被冠層對入射太陽光具有一定的反射作用，隨著植被蓋度的增大或者冠層水深的減少，沉水植被冠層在近紅外波段反射作用對水體反射光譜的貢獻度增大，導致水體反射光譜的增大。

本文根據淮河流域焦崗湖實測光譜數據探究了沉水植物菹草生長對光學淺水區水體反射光譜特征的影響規律。然而，湖泊沉水植被對水體反射光譜的影響規律因沉水植物優勢群落類型以及水體自身環境(如懸浮物濃度、葉綠素含量)的不同可能會存在一定的差異[2，11，18]，因此這些因素對湖泊光學淺水區反射光譜所產生的影響還有待于進一步的研究。