金則澎，毛峰，程乾，李軍，張軒豪

(1.浙江工商大學(xué) 旅游與城鄉(xiāng)規(guī)劃學(xué)院，杭州 310018；2.杭州電子科技大學(xué) 管理學(xué)院，杭州 310018；3.重慶師范大學(xué) 地理與旅游學(xué)院，重慶 400047)

0 引言

近幾十年來，由于我國經(jīng)濟(jì)的飛速發(fā)展導(dǎo)致湖泊資源開發(fā)利用不斷增加，在人為因素和自然因素的影響下，湖泊水環(huán)境逐年惡化[1]。淡水湖泊作為人類重要的飲用水源，也面臨著農(nóng)業(yè)、氣候變化和人類活動帶來的污染問題[2]。有色可溶解有機(jī)物(colored dissolved organic matter，CDOM)是溶解有機(jī)物(dissolved organic matter，DOM)的生物活性成分，CDOM在400～500 nm處吸收高達(dá)90%的水下太陽輻射，可以保護(hù)水下生態(tài)系統(tǒng)免受有害紫外線輻射[3-4]。但是，高濃度的CDOM對太陽光的吸收能力較強(qiáng)，會抑制水下植物的光合作用，這對淡水生態(tài)系統(tǒng)的初級生產(chǎn)產(chǎn)生了重大影響[5-6]。而且CDOM在水中具有吸附特性，會影響水污染物的分布與遷移[7]。此外，CDOM對飲用水處理過程也會帶來許多負(fù)面影響，例如CDOM會降低氧化劑和消毒劑的效力，并在氧化過程中產(chǎn)生不良的副產(chǎn)物[8]。

了解淡水中CDOM的來源、濃度和循環(huán)對于管理水生資源和預(yù)測環(huán)境變化非常重要。監(jiān)測湖泊中的CDOM，研究CDOM的時(shí)空分布能有效評價(jià)水質(zhì)，并制定出保護(hù)水資源的計(jì)劃。水質(zhì)監(jiān)測分為傳統(tǒng)理化監(jiān)測、生物監(jiān)測與遙感監(jiān)測三類，常規(guī)水質(zhì)監(jiān)測的效率不高，而且不能顯示湖泊整體的水質(zhì)情況[9-10]。目前人們已經(jīng)開發(fā)了不同類型的檢索算法來反演水質(zhì)CDOM，如經(jīng)驗(yàn)法、半解析法、矩陣反演法和機(jī)器學(xué)習(xí)法。然而，這些模型并不能很好地適用于復(fù)雜水域，特別是低濃度CDOM的內(nèi)陸水域。最近，Chen等[11]用休倫湖薩吉諾灣的數(shù)據(jù)對其中15種CDOM檢索方法進(jìn)行了比較，發(fā)現(xiàn)大多數(shù)算法在CDOM濃度較高時(shí)往往會低估其濃度，在CDOM濃度較低時(shí)往往會高估其濃度。因?yàn)樗{(lán)色光譜范圍內(nèi)的色胺吸收與由葉綠素a(chlorophyll a，Chl-a)和總懸浮物質(zhì)(total suspended matter，TSM)引起的吸收是重疊的[12]。因此，迫切需要開發(fā)一種適用于計(jì)算低濃度水域CDOM值的算法。自上個(gè)世紀(jì)以來，人工神經(jīng)網(wǎng)絡(luò)(artificial neural network，ANN)就被用于CDOM反演，并顯示出處理復(fù)雜內(nèi)陸水域的潛力。同時(shí)，基于機(jī)器學(xué)習(xí)的方法也被用于CDOM反演，并獲得了良好的效果[13]。目前已有很多基于機(jī)器學(xué)習(xí)的方法，如隨機(jī)森林回歸(random forest regression，RFR)、梯度提升回歸樹(gradient boosting regression tree，GBRT)、支持向量機(jī)回歸(support vector regression，SVR)和高斯過程回歸(gaussian process regression，GPR)，但哪一種更適合淡水湖泊中的低濃度CDOM檢測尚不得而知。GBRT的優(yōu)點(diǎn)是它可以擬合復(fù)雜的非線性關(guān)系，但不需要事先進(jìn)行數(shù)據(jù)轉(zhuǎn)換或消除孤立點(diǎn)掃描[14]。因此，在本研究中，對GBRT模型進(jìn)行了測試，并將其與其他幾種機(jī)器學(xué)習(xí)方法進(jìn)行了比較，結(jié)果表明GBRT在測試的模型中確實(shí)是最好的。

一些面向海洋的海洋水色傳感器，具有粗糙的空間分辨率，不適合在大多數(shù)河流、湖泊進(jìn)行遙感應(yīng)用。然而面向陸地的衛(wèi)星傳感器，如Landsat系列，具有很好的空間分辨率(Landsat-8為30 m)，一些基于Landsat-8的波段比模型已被提出用于低濃度的CDOM(CDOM在440 nm的吸收系數(shù)在0.066～1.242 m-1范圍內(nèi))監(jiān)測。但是Landsat系列衛(wèi)星的光譜分辨率有限，不能很精確地確定最優(yōu)波段比。因此，研究人員更加關(guān)注GF-5傳感器，GF-5衛(wèi)星影像涵蓋了可見光至短波紅外(400～2 500 nm)波段，空間分辨率為 30 m，在可見光波段的光譜分辨率可達(dá)5 nm，能提供更加精確的波段，有利于確定更加精確的波段比模型。

本文的目標(biāo)是：探索GBRT模型和GF-5衛(wèi)星影像反演千島湖低濃度CDOM的能力；描繪2019年每個(gè)季度千島湖水質(zhì)CDOM濃度的時(shí)空變化；討論影響千島湖水質(zhì)CDOM時(shí)空格局的環(huán)境因子。為實(shí)現(xiàn)這些目標(biāo)，將2018年4月至2019年9月從千島湖收集的實(shí)測數(shù)據(jù)用于GBRT算法。在2019年，利用四個(gè)GF-5衛(wèi)星影像和千島湖水文和氣象數(shù)據(jù)研究CDOM與降雨量、氣壓、風(fēng)速等環(huán)境因素的耦合關(guān)系。

1 數(shù)據(jù)和方法

1.1 研究區(qū)域

千島湖(圖1)位于浙江省西部和安徽省南部，水域面積為580 km2，有1 078個(gè)島嶼，面積大于2 500 m2，平均深度34 m，最大容量19 267 m3，如果儲水量處于108 m的正常值，則流域面積為10 480 km2[15]。記錄的年平均氣溫為17 ℃，年平均降水量為1 636.5 mm(1961—2014年)。千島湖水域內(nèi)的用地類型以林業(yè)用地為主，其次是農(nóng)業(yè)和城鎮(zhèn)建設(shè)用地(圖2)。千島湖水域水質(zhì)優(yōu)良，水體透明度高，是中國長江三角洲地區(qū)的主要飲用水源，為周圍至少一千萬人口提供了飲用水。然而，人們對千島湖水域的不合理開發(fā)引發(fā)了大量的水土流失，隨著降雨流入湖泊的工農(nóng)業(yè)和生活廢水等污染物的增多，千島湖水域的水質(zhì)情況出現(xiàn)下降的趨勢，湖泊水環(huán)境問題越來越突出[16-17]。

根據(jù)以往研究的結(jié)果，將千島湖劃分為五個(gè)次區(qū)域，以揭示千島湖不同水生環(huán)境中CDOM的時(shí)空變化：東北、西北、西南、湖泊的中部和東南部地區(qū)。以千島湖水域周邊5 km為緩沖區(qū)，進(jìn)行了土地類型分類(圖2)。

1.2 實(shí)地測量

在2018年4月19日、6月18日、2019年9月8日的三次實(shí)地取樣活動中，共從千島湖采集了40個(gè)水樣(圖1)。使用標(biāo)準(zhǔn)取水瓶在表面(深度約10 cm以下)收集水樣，在環(huán)境水溫下保存在琥珀瓶(聚丙烯250 ml)中，24小時(shí)內(nèi)運(yùn)往實(shí)驗(yàn)室分析，測量濃度等水質(zhì)參數(shù)。在每個(gè)采樣位置，使用ASD(analytical spectral devices)測量地表水光譜。FieldSpec地物光譜儀(波長范圍為325～1 075 nm，間隔為1 nm)在觀測點(diǎn)測量了水面輻亮度(Lt)、天空光輻亮度(Li)和水面下行輻照度(Ed)。指向水面的Lt傳感器的觀測角是太陽高度角40°、太陽方位角90°方向。指向天空的Li傳感器采用與Lt傳感器相同的太陽方位角，觀測的太陽天頂角是40°。為了減小不確定度，每個(gè)取樣點(diǎn)同時(shí)測量了十條光譜，并選擇光譜平均值來計(jì)算遙感反射率(Rrs)，表達(dá)如式(1)所示。

(1)

式中：水表面反射因子ρ的值受風(fēng)速影響，在風(fēng)速可忽略的平靜天氣下取值為0.022，風(fēng)速5 m/s時(shí)取值0.025；風(fēng)速10 m/s時(shí)，取值0.026～0.028。本研究考慮到千島湖的地理?xiàng)l件和水文特征，設(shè)置ρ=0.024 5。

在實(shí)驗(yàn)室低壓條件下(<5ATM)，通過GF/F玻璃微纖維膜(0.45 μm)對地表水樣品進(jìn)行過濾，以測定表面水樣的色度吸收。通過Cray-60分光輻射度計(jì)測量了200～800 nm的CDOM吸光度，測量結(jié)果為1 cm長，基線校正為Milli-Q baseline。最終，吸收系數(shù)的計(jì)算如式(2)所示。

(2)

式中：A(λ)是Cary-60分光輻射度計(jì)測量的CDOM吸光度；Path是用過的試管的路徑長度(1 cm)。在這項(xiàng)研究中，CDOM濃度是通過其在440 nm處的吸收系數(shù)來參數(shù)化的。CDOM濃度如表1所示。

表1 千島湖2018—2019年CDOM(440)野外實(shí)測數(shù)據(jù)統(tǒng)計(jì)

為了使GF-5影像導(dǎo)出aCDOM(440)，使用其光譜響應(yīng)函數(shù)RSR(relative spectral respone)將研究區(qū)實(shí)測光譜模擬出各個(gè)波段光譜，計(jì)算如式(3)所示。

(3)

式中：Rrs(Bk)是GF-5影像的第k(k=2，3，4，…，n)個(gè)波段的遙感反射率，第k個(gè)波段的Rrs(Bk)是由λm到λn波段范圍的Rrsfield積分得到;Rrsfield是現(xiàn)場測量的光譜。

1.3 水文和氣象數(shù)據(jù)

分別從屯溪和漁梁兩個(gè)水文站以及中國氣象數(shù)據(jù)網(wǎng)(http：//data.cma.cn/)收集了2019年千島湖降雨量、氣壓、風(fēng)速等氣象數(shù)據(jù)。

1.4 GF-5影像數(shù)據(jù)獲取和處理

GF-5衛(wèi)星以高空間分辨率(30 m)，高光譜分辨率(5～10 nm)的方式，對地球表面進(jìn)行觀測。2019年，共采集了六張無云GF-5衛(wèi)星影像，開發(fā)了CDOM反演算法，并對千島湖的CDOM時(shí)空變化進(jìn)行了分析。使用了一張與野外實(shí)測日期相同的2019年9月8日的GF-5衛(wèi)星影像驗(yàn)證算法。利用ENVI5.3中的大氣校正模塊(FLAASH)進(jìn)行大氣校正，得到表面反射率(Rt)。在此基礎(chǔ)上，采用5*5均值濾波，降低了圖像的不確定性。ENVI中執(zhí)行FLAASH大氣校正如式(4)所示。

(4)

式中：L是傳感器接收的像素輻亮度值；ρe是像素及其周邊像素的平均表面反射率；S是天體反照率；La是被大氣散射進(jìn)入傳感器的輻亮度；系數(shù)A和C分別取決于大氣和地理位置條件，與表面無關(guān)；大氣校正后輸出的是表面反射率影像Rt。最后計(jì)算遙感反射率(Rrs)如式(5)所示。

(5)

式中：Lr(θ，φ)是天頂角θ和方位角φ觀測角下的水表面輻亮度；Ed是向下輻照度。這兩個(gè)參數(shù)使用Hydrolight估算。計(jì)算這兩個(gè)物理量的過程中，需要在Hydrolight軟件中輸入一些參數(shù)，如太陽天頂角、影像拍攝時(shí)間、地理位置、風(fēng)速和云覆蓋百分比等。這些參數(shù)可從GF-5影像數(shù)據(jù)的頭文件和國家氣象中心提供的數(shù)據(jù)中獲取。其他需要輸入的參數(shù)(例如水質(zhì)參數(shù)濃度和水深值)對Lr(θ，φ)和Ed的大小無影響而被設(shè)置成缺省值或是任意值。

2 算法實(shí)現(xiàn)

2.1 算法背景

集成學(xué)習(xí)算法是一種強(qiáng)大的機(jī)器學(xué)習(xí)技術(shù)，已成功地用于各種分類和回歸問題，以提供更準(zhǔn)確的響應(yīng)變量估計(jì)[18]。在本研究中選擇GBRT模型，因?yàn)樗軌蜻m應(yīng)復(fù)雜的非線性關(guān)系，而且不需要事先的數(shù)據(jù)轉(zhuǎn)換或消除孤立點(diǎn)掃描。GBRT模型使用CART樹作為弱分類器，需要多次迭代。新生成的回歸樹在每次迭代中都會擬合上一棵樹的誤差，采用梯度下降法在每次迭代中向損失函數(shù)的負(fù)梯度移動，使損失函數(shù)下降。一般來說，GBRT模型的每一次迭代都會產(chǎn)生一個(gè)弱分類器，其精度不高，但集成弱分類器可以達(dá)到更高的精度。每個(gè)回歸樹的預(yù)測結(jié)果的加權(quán)和就是預(yù)測值(圖3)。

圖3 GBRT模型構(gòu)建過程

GBRT模型表達(dá)如式(6)所示。

(6)

式中：n是弱學(xué)習(xí)器的數(shù)量；θi是系數(shù)(減少過擬合)；fi是弱學(xué)習(xí)器；Fm是最終的一般模型。GBRT模型有幾個(gè)重要的參數(shù)：每個(gè)弱學(xué)習(xí)器的最大深度(一般不超過五個(gè))、弱學(xué)習(xí)器的最大數(shù)量，以及學(xué)習(xí)率(更高的學(xué)習(xí)率意味著更強(qiáng)的修正，使模型更加復(fù)雜)。所以需要選擇適合的參數(shù)來提高模型的性能。

2.2 算法架構(gòu)

涉及藍(lán)、綠和紅域遙感的波段比算法已被廣泛用于淡水湖的CDOM遙感[19]。此外，與使用單波段反射率相比，波段比模型可以減少大氣校正的更多不確定性[20]。本研究以B25(Rrs493)、B41(Rrs561)和B56(Rrs625)組成的GF-5影像波段比作為GBRT模型的輸入，通過決定系數(shù)(R2)、均方根誤差(RMSE)、偏差和平均絕對百分比誤差(MAPE)四個(gè)指標(biāo)確定最佳組合。

現(xiàn)場測量的40個(gè)樣本被隨機(jī)分成兩個(gè)獨(dú)立的數(shù)據(jù)集，它們的樣本比例約為3∶1，這是機(jī)器學(xué)習(xí)中常用的比例。因此，GBRT的訓(xùn)練和驗(yàn)證數(shù)據(jù)集分別包含30個(gè)和10個(gè)現(xiàn)場測量樣本?；A(chǔ)學(xué)習(xí)者的數(shù)量為40人，學(xué)習(xí)率為0.1，每個(gè)學(xué)習(xí)者的最大深度為5，其他參數(shù)默認(rèn)配置。

2.3 與其他CDOM檢索算法的比較

為了驗(yàn)證GBRT在估算淡水環(huán)境中低濃度CDOM方面的優(yōu)勢，還對以往的兩種傳統(tǒng)模型進(jìn)行了比較。對數(shù)模型，參見式(7)，其中aCDOM(440)的范圍為0.51～25.1 m-1[21]。另一個(gè)模型是一元二次多項(xiàng)式模型，見式(8)，其中aCDOM(440)的范圍為0.644～1.413 m-1[22]。第二個(gè)模型在水體類型和CDOM范圍方面與構(gòu)建模型更具可比性。

(7)

(8)

式中：B25、B41和B56是GF-5影像數(shù)據(jù)的波段。對于不同的研究地點(diǎn)，參數(shù)D、E和F需要重新校準(zhǔn)。

3 結(jié)果分析

3.1 野外測量的CDOM和地表光譜

在以前的許多研究地點(diǎn)，內(nèi)陸水域中的CDOM濃度范圍很大。例如，在明尼蘇達(dá)州的15個(gè)湖泊中，aCDOM(440)的范圍為0.6～19.4 m-1[23]。與這些CDOM濃度變化較大的湖泊相比，2019年千島湖的CDOM濃度變化相對較低。野外實(shí)測的aCDOM(440)范圍為0.103～1.366 m-1，平均值為0.53 m-1。野外實(shí)測數(shù)據(jù)顯示出明顯的季節(jié)變化，4月份的aCDOM(440)在0.103～0.507 m-1(平均0.253 m-1)范圍內(nèi)，而9月份的變化范圍為0.196～1.366 m-1(平均0.637 m-1)，說明秋季的CDOM濃度高于早春。

千島湖實(shí)測的遙感反射率如圖4所示，顯示了復(fù)雜內(nèi)陸淡水的典型光譜特征。浮游植物對葉綠素和胡蘿卜素的低吸收和低散射產(chǎn)生了570 nm的峰值[24]。670 nm處的小反射谷可能是由葉綠素a對紅光波段強(qiáng)烈吸收引起的[25]。在千島湖，由藻類色素和水體自身吸收較小引起的700 nm附近的葉綠素?zé)晒獠⒉伙@著。

圖4 千島湖實(shí)測地表光譜(Rrs)

3.2 模型評估和比較

為了確定GBRT模型的最佳輸入，對所有可能的波段比組合進(jìn)行了計(jì)算?？紤]四個(gè)指標(biāo)(R2、MAPE、RMSE、偏差)，最終選擇了B25/B41和B56/B25來反演CDOM濃度。最佳的GBRT模型的R2為0.90和0.93，MAPE為0.34%和15.53%，RMSE為0.35 m-1和0.16 m-1，培訓(xùn)和驗(yàn)證數(shù)據(jù)集的偏差為-0.004和-0.032，測量到的訓(xùn)練和驗(yàn)證數(shù)據(jù)集最佳模型的CDOM濃度散點(diǎn)圖如圖5(a)和圖5(b)所示。結(jié)果表明，該算法的性能可以用于反演千島湖的CDOM濃度。

圖5 模型驗(yàn)證結(jié)果

本研究利用在千島湖的實(shí)測資料，對前人建立的兩個(gè)模型進(jìn)行了檢驗(yàn)，研究發(fā)現(xiàn)對數(shù)模型(R2=0.357 4，RMSE=0.394 5 m-1)(圖6(a))和多項(xiàng)式模型(R2=0.048 4，RMSE=0.414 7 m-1)(圖6(b))。在千島湖不能很好地?cái)M合。將GBRT與RFR、SVR、GPR進(jìn)行了比較，所有模型均選擇波段比B25/B41和B56/B25作為輸入。表2顯示了統(tǒng)計(jì)分析的結(jié)果。結(jié)果表明，支持向量回歸對訓(xùn)練數(shù)據(jù)集(R2=0.44，MAPE=64.5%)和驗(yàn)證數(shù)據(jù)集(R2=0.27，MAPE=68.3%)都表現(xiàn)不佳，而高斯過程回歸對驗(yàn)證數(shù)據(jù)集產(chǎn)生了過擬合(R2=0.02，MAPE=86.3%)。與隨機(jī)森林回歸相比，對于驗(yàn)證數(shù)據(jù)集，梯度提升回歸樹的R2增加了20%，而MAPE降低了26.17%。

圖6 檢驗(yàn)千島湖低濃度CDOM反演的兩個(gè)經(jīng)驗(yàn)?zāi)Ｐ?/p>

表2 四種典型機(jī)器學(xué)習(xí)算法的比較

3.3 CDOM的時(shí)空變化

將訓(xùn)練好的GBRT模型應(yīng)用于2019年的四幅GF-5影像，以估計(jì)千島湖的CDOM濃度，如圖7所示。為了進(jìn)一步分析CDOM的時(shí)空變化，根據(jù)每個(gè)區(qū)域的有效像素估計(jì)了各個(gè)子區(qū)域每個(gè)季度的平均aCDOM(440) 濃度(圖8)。結(jié)果表明，2019年千島湖CDOM的季平均濃度較低，為春季0.61 m-1、夏季0.65 m-1、秋季0.75 m-1和冬季0.54 m-1，CDOM濃度季節(jié)性差異較為顯著，其中秋季明顯高于其他三季。四幅反演圖像中，CDOM濃度高值區(qū)與千島湖流域用地分類圖(圖2)中的建筑/道路、農(nóng)田用地相契合，表明了人類活動對于水體CDOM濃度影響較大。千島湖西北區(qū)和西南區(qū)的CDOM含量明顯高于其他地區(qū)。西北區(qū)的新安江占千島湖總徑流量的60%，西南區(qū)的武強(qiáng)溪占千島湖總徑流量的25%。眾所周知，徑流會將陸地腐殖質(zhì)帶入湖泊，使得大量的CDOM通過武強(qiáng)溪和新安江匯入千島湖，西北和西南區(qū)河口的水流流速比中部地區(qū)的水流流速大，河口區(qū)域的水體充分混合，使得水體分層不明顯，因?yàn)樽贤夤鉄o法輕易穿透表層水體，因此光漂白速率會比那些分層明顯的中部區(qū)域低，從而導(dǎo)致西北和西南區(qū)的CDOM濃度相對較高。

圖8 千島湖季平均aCDOM(440)濃度

3.4 CDOM與水文氣象因子的關(guān)系

以前的一些研究表明，降水量對水中CDOM濃度有重大影響[26]。CDOM濃度與降水量的R2為0.620 8，表明降水量對千島湖CDOM濃度的變化影響顯著。千島湖地區(qū)年間降水量不均勻，水位有明顯落差，其平均水深為37 m，最大水深為90 m。豐水期的降水量比枯水期的降水量高得多，高度動態(tài)變化的水位會使得湖泊底質(zhì)變化，而水體底部有機(jī)質(zhì)的釋放是水中CDOM的來源之一。豐水期水流速度大，河流中的CDOM會流入主湖區(qū)；而枯水期水域的水位低，水流速度不大，從河流流入的CDOM就會減少。夏季和秋季正是千島湖的豐水期，因此其CDOM濃度高于春季和冬季。然而，在本研究中，千島湖的CDOM濃度峰值出現(xiàn)在秋季，而不是在降雨量最大的夏季(圖9)。

圖9 千島湖季平均降雨量aCDOM(440)之間的相關(guān)關(guān)系

眾所周知，光化學(xué)降解和光漂白是去除CDOM的重要機(jī)制。例如，在太湖進(jìn)行的一項(xiàng)為期12 d的實(shí)驗(yàn)發(fā)現(xiàn)，當(dāng)CDOM暴露在太陽輻射下時(shí)，CDOM的吸收量aCDOM(355)和aCDOM(280)分別下降了29.8%和20.8%[27]。2019年，千島湖在夏季(5—7月)經(jīng)歷了高強(qiáng)度的UV-B輻射，這顯著增加地表水的光漂白強(qiáng)度，從而消除部分降雨和徑流帶來的CDOM。因此，2019年夏季千島湖反演的CDOM濃度并不太高。相比之下，秋季的CDOM濃度較高可能是由于以下三個(gè)原因。①經(jīng)過連續(xù)多個(gè)月的降雨后，地表徑流和地下水逐漸將豐富的土壤有機(jī)質(zhì)釋放到水中。②秋季太陽輻射減弱，減緩了CDOM的光漂白。③秋季越來越多的落葉腐爛，直接增加了土壤和水中的有機(jī)質(zhì)含量；冬季是一年中降雨量最少的季節(jié)，加上溫度下降導(dǎo)致的水面結(jié)冰，徑流量大大減少，CDOM濃度處于一年中最低；春季由于降雨量的回升，以及冰雪融化加大了地表徑流，導(dǎo)致外源有機(jī)物的輸入加大，CDOM濃度相比較冬季有所上升。

分析衛(wèi)星獲取的CDOM濃度與水文和氣象因子(氣壓和風(fēng)速)之間的相關(guān)性發(fā)現(xiàn)，CDOM濃度與氣壓和風(fēng)速的R2分別為0.082 5和0.096 4(圖10)，表明這些因素對千島湖CDOM濃度的變化沒有顯著影響。

圖10 千島湖季平均氣壓、風(fēng)速以及與aCDOM(440)之間的相關(guān)關(guān)系

4 結(jié)束語

利用梯度提升回歸樹算法，結(jié)合GF-5高光譜衛(wèi)星影像建立了千島湖低濃度CDOM遙感模型。經(jīng)比較驗(yàn)證，模型性能良好，RMSE=0.1 m-1，MAPE=19.0%。然后，應(yīng)用GBRT模型預(yù)測了2019年千島湖CDOM濃度的時(shí)空變化。千島湖整體aCDOM(440)較低，最低為0.005 m-1，最高為1.472 m-1，CDOM濃度季節(jié)差異較為顯著，秋季>夏季>春季>冬季。各個(gè)季節(jié)高CDOM濃度部分一般位于湖的邊緣區(qū)，主要是入水口以及湖周圍與人類活動密切相關(guān)的地方，如西北和西南區(qū)的入水口附近的CDOM濃度較高。最后結(jié)合水文數(shù)據(jù)，發(fā)現(xiàn)CDOM濃度與氣壓和風(fēng)速的相關(guān)性很低(R2分別為0.082 5和0.096 4)，表明這些氣象因素對千島湖CDOM濃度的變化影響不大。

由于現(xiàn)場測量和觀測僅在一年內(nèi)進(jìn)行，更多的環(huán)境因素或事件可能會改變CDOM和水質(zhì)，因此GBRT模型仍需進(jìn)一步改進(jìn)。本研究表明，GBRT模型和GF-5影像能準(zhǔn)確監(jiān)測低濃度湖泊中CDOM的時(shí)空變化，可以更好地了解河流、湖泊系統(tǒng)中可溶解有機(jī)質(zhì)與其環(huán)境因子之間的耦合關(guān)系。