其中,沉水植物干重生物量用沉水植物濕重生物量乘0.08表示[17],藻類干重生物量通過Chl.a濃度乘70再乘湖泊平均水深得到[18].
根據(jù)Aizaki修正的營養(yǎng)狀態(tài)指數(shù)(trophic state index,TSIM)來評價研究水域的富營養(yǎng)化狀態(tài),并用加權(quán)平均計算處理,其中Chl.a、TP和透明度的比重分別為54%、29.7%和16.3%,評價標準為:TSIM<37為貧營養(yǎng),37≤TSIM<53為中營養(yǎng),53≤TSIM<65為富營養(yǎng),TSIM≥65為重富營養(yǎng)[19-21].
1.6 統(tǒng)計分析
采用One-way ANOVA 對不同區(qū)域光學特性和衰減因子的差異性進行分析. 線性回歸運用SPSS 20.0進行分析,采用向后多元回歸方法篩選不同區(qū)域內(nèi)的主導衰減因子. 采用平均相對誤差(MAPE)和均方根誤差(RMSE)兩個指標對多元回歸模型精度驗證結(jié)果進行評價,計算公式為:
(5)
(6)
式中,n為樣點數(shù)目,Esi和Oi分別為第i個樣點的模型估算值和實測值.
2 結(jié)果與分析
2.1 光合有效輻射垂向特征
由于PAR測量數(shù)據(jù)較多,因此在3種區(qū)域選擇了6個采樣點(1#、3#、4#、5#、6#、8#)的PAR垂向分布特征及指數(shù)回歸結(jié)果進行展示(圖2). 在現(xiàn)場測量過程中,明顯的異常值被舍棄,重新進行測量. 各采樣點的指數(shù)擬合決定系數(shù)(R2)都達到了0.9以上,其中,無植物區(qū)域擬合效果較好,沉水植物分布區(qū)域內(nèi)擬合效果較差,測量時水下植物的遮擋是影響擬合效果的主要原因,天氣狀況、波浪也可能對結(jié)果產(chǎn)生一定影響.
圖2 6個代表點位PAR強度垂直分布特征Fig.2 Vertical distribution characteristics of PAR intensity of six typical sampling sites
2.2 光學特性與影響因子的空間差異及線性擬合
盡管東太湖面積很小,但不同區(qū)域kd(PAR)的差異比較大(表1). 8月份各采樣點的kd(PAR)均值在1.04~6.35 m-1之間變化,極差為5.31 m-1,其中6#、7#、8#點位kd(PAR)均值較小,范圍為1.04~1.16 m-1,該區(qū)域為沉水植物茂盛區(qū);1#、2#、5#、10#點位kd(PAR)均值較大,范圍為1.39~3.47 m-1,該區(qū)域為沉水植物稀疏區(qū);3#、4#、9#點位kd(PAR)均值最大,范圍為5.97~6.35 m-1,該區(qū)域為無植物區(qū).
表1 東太湖kd(PAR)的計算結(jié)果
如圖3所示,東太湖kd(PAR)的空間變化規(guī)律為:植物茂盛區(qū)(1.09 m-1)<植物稀疏區(qū)(2.50 m-1)<無植物區(qū)(6.19 m-1).kd(PAR)在植物茂盛區(qū)最小,變化范圍為0.73~1.89 m-1,平均值僅為1.09 m-1,顯著低于植物稀疏區(qū)和無植物區(qū)(P<0.05);kd(PAR)在無植物區(qū)最大,變化范圍為4.41~11.80 m-1,平均值為6.19 m-1,顯著高于植物茂盛區(qū)和植物稀疏區(qū)(P<0.05). 真光層深度與kd(PAR)相反,空間變化規(guī)律為無植物區(qū)(0.78 m)<植物稀疏區(qū)(2.25 m)<植物茂盛區(qū)(4.38 m),變化范圍為0.39~6.31 m,3種區(qū)域間存在顯著性差異(P<0.05).
圖3 東太湖水體kd(PAR)、真光層深度、ISS、OSS、葉綠素a和CDOM吸收系數(shù)箱線圖(實線上、下端代表最大值和最小值;矩形的上、下端代表上四分位數(shù)和下四分位數(shù);中間橫線為中位線;小正方形為均值;圖中字母表示Tamhane’s T2多重比較檢驗結(jié)果,字母不同表示差異顯著(P<0.05))Fig.3 Box plot of kd(PAR), euphotic depth, suspended particulate inorganic matter, suspended particulate organic matter, chlorophyll-a and CDOM absorption coefficient in eastern Lake Taihu
ISS和OSS濃度的變化規(guī)律一致:植物茂盛區(qū)(ISS:1.95 mg/L,OSS:1.63 mg/L)<植物稀疏區(qū)(ISS:11.91 mg/L,OSS:4.00 mg/L)<無植物區(qū)(ISS:60.85 mg/L,OSS:9.29 mg/L),二者濃度范圍分別為0.33~119.27 mg/L、0.51~17.31 mg/L,3種區(qū)域間存在顯著性差異(P<0.05);Chl.a濃度變化規(guī)律為植物茂盛區(qū)(9.61 μg/L)<植物稀疏區(qū)(22.08 μg/L)<無植物區(qū)(23.74 μg/L),范圍在4.79~58.21 μg/L之間,植物稀疏區(qū)與植物茂盛區(qū)差異顯著(P<0.05),而與無植物區(qū)沒有顯著差異(P>0.05);αCDOM(440)變化規(guī)律為植物稀疏區(qū)(0.63 m-1)<植物茂盛區(qū)(0.69 m-1)<無植物區(qū)(0.76 m-1),變化范圍為0.10~2.59 m-1,且3種區(qū)域無顯著性差異(P>0.05).
與無植物區(qū)相比,沉水植物茂盛區(qū)kd(PAR)減少了83.4%,光衰減因子中,ISS、OSS、Chl.a濃度和αCDOM(440)分別減少了96.8%、82.5%、59.5%和9.2%.
分別將整個東太湖和3種區(qū)域的kd(PAR)與ISS、OSS、Chl.a濃度和CDOM進行線性擬合,結(jié)果發(fā)現(xiàn)東太湖kd(PAR)與ISS濃度有很好的線性關(guān)系,R2為0.947(圖4a),與OSS濃度的線性關(guān)系也比較好,R2為0.853(圖4b),但與Chl.a濃度的線性關(guān)系較差,R2為0.151(圖4c),與CDOM吸收系數(shù)無顯著線性關(guān)系,R2為0.022(圖4d). 茂盛區(qū)kd(PAR)與各影響因子均沒有顯著線性關(guān)系(圖5);稀疏區(qū)kd(PAR)與ISS、OSS濃度呈顯著正線性關(guān)系,與Chl.a濃度、CDOM吸收系數(shù)無顯著線性關(guān)系(圖6);無植物區(qū)kd(PAR)與ISS、OSS濃度呈顯著正線性關(guān)系,與Chl.a濃度呈顯著負線性關(guān)系,與CDOM吸收系數(shù)無顯著線性關(guān)系(圖7).
圖4 東太湖kd(PAR)與無機懸浮物(a)、有機懸浮物(b)、葉綠素a(c)和CDOM吸收系數(shù)(d)的線性擬合Fig.4 Linear fitting between kd(PAR) and suspended particulate inorganic matter (a), suspended particulate organic matter (b), chlorophyll-a (c) and CDOM absorption coefficient (d) in eastern Lake Taihu
2.3 不同區(qū)域的主導影響因子
上節(jié)分析表明,東太湖不同區(qū)域的kd(PAR)差異顯著,衰減因子也有較大的區(qū)別. 通過向后法進行多元線性回歸,對每種區(qū)域光衰減的主導因子進行分析. 沉水植物茂盛區(qū)回歸模型先后剔除了αCDOM(440)、Chl.a、ISS和OSS濃度,表明影響該區(qū)域水下光衰減的因素主要是水體本身的吸收和散射;沉水植物稀疏區(qū)的kd(PAR)與衰減因子的表達式見表2模型1,模型先后剔除了OSS濃度和αCDOM(440),保留了ISS和Chl.a濃度,說明該區(qū)域水下光照衰減由ISS和Chl.a濃度共同主導,模型精度較高,擬合結(jié)果R2為0.883,驗證結(jié)果MAPE和RMSE分別為15.8%和0.412;表2中模型2為無植物區(qū)域的kd(PAR)回歸方程,模型先后剔除了OSS濃度、αCDOM(440)和Chl.a濃度,僅保留了ISS濃度,因此,該區(qū)域的ISS濃度是水下光照衰減的主導因子,模型精度也比較高,擬合結(jié)果R2為0.840,驗證結(jié)果MAPE和RMSE分別為6.3%和0.636.
圖5 植物茂盛區(qū)kd(PAR)與無機懸浮物(a)、有機懸浮物(b)、葉綠素a(c)和CDOM吸收系數(shù)(d)的線性擬合Fig.5 Linear fitting between kd(PAR)and suspended particulate inorganic matter (a), suspended particulate organic matter (b), chlorophyll-a (c) and CDOM absorption coefficient (d) in submerged plants sparse region
圖6 植物稀疏區(qū)kd(PAR)與無機懸浮物(a)、有機懸浮物(b)、葉綠素a(c)和CDOM吸收系數(shù)(d)的線性擬合Fig.6 Linear fitting between kd(PAR) and suspended particulate inorganic matter (a), suspended particulate organic matter (b), chlorophyll-a (c) and CDOM absorption coefficient (d) in submerged plants lush region
圖7 無植物區(qū)kd(PAR)與無機懸浮物(a)、有機懸浮物(b)、葉綠素a(c)和CDOM吸收系數(shù)(d)的線性擬合Fig.7 Linear fitting between kd(PAR) and suspended particulate inorganic matter (a), suspended particulate organic matter (b), chlorophyll-a (c) and CDOM absorption coefficient (d) in no-plant region
表2 多元線性回歸模型驗證結(jié)果的誤差統(tǒng)計
2.4 東太湖草藻干重比和區(qū)域富營養(yǎng)化評價
表3記錄了采樣期間各采樣點沉水植物生物量和優(yōu)勢種,根據(jù)沉水植物生物量和葉綠素a濃度數(shù)據(jù),計算出各采樣點草藻干重比,見表4. 采用TSIM方法評價區(qū)域富營養(yǎng)化狀態(tài),見表5. 總體上講,東太湖目前處于由草型湖泊向藻型湖泊過渡的階段,并且湖泊處于富營養(yǎng)水平,相比以前的草型穩(wěn)態(tài)[22],東太湖面臨著向藻型穩(wěn)態(tài)轉(zhuǎn)變的趨勢. 其中,沉水植物茂盛區(qū)草藻干重比大于100,屬于草型穩(wěn)態(tài),區(qū)域處于中營養(yǎng)水平;稀疏區(qū)草藻干重比在1~100之間,屬于草藻共生狀態(tài),區(qū)域處于富營養(yǎng)水平;無植物區(qū)位于行洪通道區(qū)域,沒有沉水植物分布,屬于藻型穩(wěn)態(tài),處于重富營養(yǎng)水平.
表3 沉水植物生物量和優(yōu)勢種
表4 東太湖各采樣點草藻干重比
表5 3種區(qū)域富營養(yǎng)化指標值及評價結(jié)果
3 討論
東太湖面積雖小,但本研究發(fā)現(xiàn),不同區(qū)域的kd(PAR)存在顯著性差異,各采樣點均值在1.04~6.35 m-1之間變化. 其中沉水植物分布茂盛的區(qū)域kd(PAR)最小,主要是因為沉水植物一方面可以減小波浪的沖刷作用、防止底泥懸浮、降低水體中的懸浮物濃度,另一方面沉水植物能夠分泌化感物質(zhì)抑制藻類的生長,降低Chl.a濃度[23-24]. 采樣期間臺風對該區(qū)域基本沒有影響,kd(PAR)變化不明顯. 行洪通道區(qū)域kd(PAR)最大,區(qū)域內(nèi)沒有沉水植物,流速和波浪比較大,水體渾濁. 該區(qū)域承接大太湖的流水,水體懸浮物濃度比較高,再加上風浪的擾動,懸浮物沉降較慢,風浪較大時底泥還會懸浮,因此懸浮物濃度最高. 沉水植物相對稀疏的區(qū)域基本位于圍網(wǎng)拆除區(qū)域內(nèi)部,其中1#點在東山沿岸附近,該區(qū)域的kd(PAR)較大. 圍網(wǎng)拆除區(qū)域的點位會受行洪通道區(qū)域的影響,中心區(qū)域波浪較大時會向兩邊傳播,雖然沉水植物可以固著底泥,但波浪引起的切應力大于底泥起動應力時就會造成底泥的再懸浮,臺風過后可以觀察到水體懸浮物濃度明顯增加.
雖然kd(PAR)屬于表觀光學特性,但是主要由固有光學特性決定,盡管它受到入射光場變化的影響,但并不十分敏感[25]. 張運林等[26]通過在不同表面光強下得出的光學衰減系數(shù)分析發(fā)現(xiàn),其基本上與表面光強無關(guān),不隨光場條件的改變而改變,主要受制于水中的組成物質(zhì). 不同湖泊光衰減的主要因素也有區(qū)別. van Duin等估算了荷蘭20個淺水湖泊中各成分對衰減系數(shù)的貢獻,發(fā)現(xiàn)在大多數(shù)湖泊,光衰減主要是由藻類生物量引起的[27]. 美國佛羅里達洲的一個大型富營養(yǎng)化淺水湖泊Apopka湖中,Chl.a對光衰減的貢獻率甚至超過50%. Erikson等研究發(fā)現(xiàn)Xolotlán湖的水下光照深度主要受Chl.a濃度的調(diào)節(jié)[28]. 還有一些研究認為懸浮物濃度是影響光衰減的重要因素[29-30].
然而我們的研究發(fā)現(xiàn),東太湖不同區(qū)域的主導衰減因子有所不同. 在沉水植物分布茂盛的區(qū)域,光合有效輻射衰減系數(shù)模型剔除了所有變量,區(qū)域內(nèi)光衰減系數(shù)較小,主要是水體本身對光的吸收和散射. 在沉水植物分布稀疏的區(qū)域,ISS和Chl.a濃度較高,二者共同主導了光照的衰減,該區(qū)域?qū)儆诓菰骞采鸂顟B(tài). 雖然也有沉水植物分布,但是生物量比較小,對底泥和藻類的抑制作用相對減弱,在外力的擾動下容易發(fā)生底泥懸浮,釋放營養(yǎng)鹽,為藻類生長提供營養(yǎng)物質(zhì). 以往對太湖的研究表明,水下光照衰減主要由懸浮物引起[11,30],但我們對東太湖沉水植物稀疏區(qū)域的研究發(fā)現(xiàn),Chl.a也是影響水下光照衰減的重要因素之一,在生態(tài)修復實踐中不容忽視. 近幾年沉水植物的逐漸減少表明東太湖目前可能接近穩(wěn)定轉(zhuǎn)變的關(guān)鍵階段,從大型水生植物主導的清澈狀態(tài)向浮游植物主導的渾濁狀態(tài)轉(zhuǎn)變將導致許多淺水湖泊水質(zhì)持續(xù)不良,然而,一旦轉(zhuǎn)變后再恢復大型水生植物主導的清澈狀態(tài)往往更加困難[9]. 因此需要重點關(guān)注此區(qū)域,減小波浪對底泥擾動造成的底泥懸浮和營養(yǎng)鹽釋放,降低湖水中懸浮物和Chl.a濃度. 在無植物區(qū)域,ISS濃度要比Chl.a濃度大得多,是影響光衰減的主要因素. 該區(qū)域底泥容易再懸浮,在水流的輸移作用下泥沙會向植物區(qū)輸送,進而影響到兩側(cè)沉水植物的生長. 3個區(qū)域的OSS和CDOM影響都比較小,一方面是因為東太湖OSS和CDOM濃度比較低,懸浮物主要由ISS構(gòu)成,OSS僅占到總懸浮物的17%,導致引起水下光照衰減的影響無法顯現(xiàn),另一方面CDOM對光的吸收在400 nm以下的紫外短波部分貢獻比較高,500 nm以上的光衰減主要來源于懸浮物和Chl.a.
真光層深度是湖泊物理、化學、生物綜合作用的關(guān)鍵指標,在湖泊生態(tài)修復和水質(zhì)改善方面發(fā)揮著重要作用[31]. 真光層深度與水深的比值可以作為劃分是否可恢復沉水植物的標準,比值大于1時,表明水下光照環(huán)境可滿足沉水植物的生長需求. 根據(jù)真光層深度與透明度的關(guān)系,在實際測量中可以通過測量湖水透明度來判斷區(qū)域內(nèi)沉水植物能否生存. 本文根據(jù)計算的真光層深度和測量的透明度數(shù)據(jù),建立了回歸方程,結(jié)果表明在置信區(qū)間為95%的水平下,真光層深度和水體透明度呈顯著的正線性關(guān)系:Deu=2.70SD(R2=0.969,P≤0.0001,n=160). 真光層深度大約是透明度的2.7倍,在實際應用中,可通過測量透明度反算出真光層深度,然后比較真光層深度與水深的大小,當真光層深度≥水深時,表明可以恢復沉水植物. 采樣期間東太湖平均水深1.85 m,即將透明度提高到0.7 m以上,就有可能恢復沉水植物,因此當?shù)厣鷳B(tài)修復工程中可借鑒此數(shù)據(jù). 總的來說,東太湖光照衰減的主要因素是ISS和Chl.a,水下光照的分布決定了沉水植物的分布情況. 目前,東太湖正處于從草型湖區(qū)到藻型湖區(qū)過渡的階段,尤其在沉水植物稀疏區(qū)域,圍網(wǎng)拆除后增加的水動力對沉水植物生長造成了不利影響. 另外,該區(qū)域長有許多浮葉植物,阻礙了光照在水下的傳播,進一步影響了沉水植物的光照需求. 因此,建議可采取圍隔等消浪措施減小水動力,防止底泥懸浮,提高水體透明度,同時加強水生植物管理,適度收割浮葉植物,改善水下光照環(huán)境,保證沉水植物生長良好,促進湖水水體的良性循環(huán). 鑒于大太湖藻類輸入的影響,有關(guān)部門需要及時關(guān)注和監(jiān)測東太湖入口處——東茭嘴附近的水質(zhì)情況,采取相關(guān)措施阻止外源性藻類入侵東太湖,削弱藻類與沉水植物的競爭力.
4 結(jié)論
1)東太湖kd(PAR)范圍為0.73~11.8 m-1,空間變化規(guī)律為:植物茂盛區(qū)(1.09 m-1)<植物稀疏區(qū)(2.50 m-1)<無植物區(qū)(6.19 m-1),對應的真光層深度范圍為0.39~6.31 m. 東太湖光學特性空間差異顯著. 衰減因子中,ISS和OSS濃度空間差異顯著(P<0.05),CDOM吸收系數(shù)無顯著性差異(P>0.05),沉水植物稀疏區(qū)葉綠素a濃度顯著高于植物茂盛區(qū)(P<0.05),但與無植物區(qū)沒有顯著差異(P>0.05). 沉水植物的分布影響著區(qū)域內(nèi)的水質(zhì)因子變化.
2)kd(PAR)與ISS、OSS濃度的線性擬合較好,與Chl.a、CDOM擬合較差;沉水植物茂盛區(qū)衰減主要是水體本身對光的吸收和散射作用,稀疏區(qū)屬于草藻共生狀態(tài),由Chl.a和ISS共同主導,無植物區(qū)域主導衰減因子為ISS,整個區(qū)域的OSS和CDOM對光照衰減基本無影響. ISS和Chl.a濃度是影響該區(qū)域光照衰減的主要因子.
3)東太湖正處于由草型穩(wěn)態(tài)向藻型穩(wěn)態(tài)過渡的階段,整個湖泊屬于富營養(yǎng)水平,其中沉水植物茂盛區(qū)屬于草型穩(wěn)態(tài),處于中營養(yǎng)水平;稀疏區(qū)屬于草藻共生狀態(tài),處于富營養(yǎng)水平;無植物區(qū)屬于藻型穩(wěn)態(tài),處于重富營養(yǎng)水平.
4)東太湖真光層深度大約是透明度的2.7倍.