淡水魚塘葉綠素a與環境因子間影響關系的通徑分析

劉金金,張玉平,張 丹,孫振中,李曉蓓,晏 軍,趙 冉,張 芬,周 軒

(上海市水產研究所,上海 200433)

浮游植物是池塘生態系統中重要的初級生產力，是養殖對象重要的食物來源[1-2]。藻類的生長狀況影響養殖水體的水質演化，其過度繁殖可能產生水華，惡化水質，甚至其分泌的藻毒素會威脅到養殖對象的健康及品質[3-4]。葉綠素a(Chl.a)是浮游植物現存量的重要評價指標[5]，其含量高低與水環境質量密切相關[6]，已有眾多學者在此方面做了大量研究，研究主體主要為湖泊、水庫、河流等水體[7-10]，池塘養殖水體因其半封閉性、人為干涉頻繁等特點，水環境因子復雜多變，與Chl.a之間的關系更加復雜，相關方面的研究相對較少[3-11]。

目前，池塘養殖水體中Chl.a與水環境因子間相關關系的分析方法主要為多元線性回歸分析、灰色關聯分析、主成分分析以及模糊綜合評價[12-14]，其中模糊綜合評價法和灰色聚類分析等方法本身存在對指標評價的主觀性判斷和信息丟失的問題[15]，而通過主成分分析簡化數據信息，剔除部分不適合因子后再逐步建立多元回歸方程的方法所需數據量大，計算繁瑣，且得到的結果較難一致性地描述因子間的關系和規律[16]，因此尋求一種簡便、高效而又科學客觀的分析方法顯得尤為重要。通徑分析是一種可處理較為復雜數據間變量關系的線性回歸分析方法，通徑系數可有效分析相關變量之間影響的直接或間接效應，決定系數可反映出自變量對因變量的作用大小。通徑分析比模糊綜合分析和灰色關聯分析更為全面，比主成分分析法在數據需求量上更少，在計算操作上更加便捷，是一種相對較方便且有效的分析方法，但該方法多用于分類學，在養殖水環境中的應用還極少見。以上海市淡水魚池塘養殖水體水質監測數據為基礎，利用通徑分析方法研究魚塘水體中Chl.a與水環境因子間的關系，探討魚塘水體中影響Chl.a含量變化的主要因素，并建立具有針對性的多元回歸方程，以期能夠更加全面地研究池塘水體中環境因子間的內在影響規律，從而為上海地區淡水魚塘的生態養殖提供一定科學論據。

1 材料與方法

1.1 樣品采集

表1 上海市淡水魚養殖場采樣池塘養殖信息匯總表

1.2 水質測定

1.3 分析方法

利用SPSS statistics 19.0 軟件對水質監測結果進行描述性統計和相關性分析。根據相關性系數(rij)矩陣建立正則方程組，利用高斯消元法得到各環境因子對Chl.a影響的直接通徑系數(Pi)。根據相關矩陣的組成原理，獲得各環境因子通過其他環境因子對Chl.a產生影響的間接通徑系數(∑rijPj)，依次對各水環境因子對Chl.a的影響展開通徑分析，并計算各因子對Chl.a產生影響的決定系數(di)以及通過其他因子影響Chl.a的共同決定系數(dij)。最后通過檢驗偏回歸系數剔除不顯著因子建立多元回歸方程，并對多元回歸模型進行了假設性檢驗。

(1)

相關系數的組成原理：rxiy=Pi+∑rijPj

(2)

(3)

共同決定系數(dij)的計算：

dij=2rij×Pi×Pj

(4)

2 結果與分析

2.1 魚塘水環境概況

魚塘水質因子結果的描述性統計見表2。水溫主要隨氣候變化，同一監測時間內池塘水體間差異不大，水體鹽度和pH也無較大差異，透明度和溶解氧的變化范圍稍大，但變異系數均小于50%，其中pH的數據離散程度最小，表明在淡水魚的養殖周期內水體酸堿性均維持在一個較為穩定的水平；氨氮、亞硝酸鹽氮、硝酸鹽氮、總磷、活性磷等因子的變異系數大于85%，數據離散程度大，其中活性磷的數據離散程度最大。總氮和高錳酸鹽指數的變異系數在50%左右，數值變化范圍相對不大。

表2 淡水魚養殖池塘水環境因子描述性統計結果匯總

2.2 水環境因子間相關系數

表3 淡水魚池塘水環境因子間相關系數矩陣

附：r0.05(1,38)=0.320,r0.01(1,38)=0.413,上標“*”顯著性相關,“**”極顯著性相關.

2.3 通徑分析

2.3.1 關鍵因子的篩選

表4 淡水魚池塘環境因子變化對Chl.a影響的通徑系數匯總

附：** 為極顯著水平，P<0.01;* 為顯著水平，P<0.05.

2.3.2 關鍵因子對Chl.a的通徑分析

表5 淡水魚池塘重要環境因子變化對Chl.a影響的通徑系數

表6 淡水魚池塘重要環境因子變化對Chl.a影響的決定系數匯總

2.4 多元回歸方程的建立

(5)

對多元回歸方程做方差分析(表7)，回歸系數達到極顯著水平(P<0.01)，擬合優度(R=0.784)較高，方程理論上有效。利用所得回歸方程及四個自變量計算Chl.a，將計算值與實測值進行對比，Chl.a的計算值與實測值之間滿足回歸方程：y=0.9889x+2.426，散點示意圖見圖1，表明利用上述多元回歸方程得到的Chl.a的計算值是可接受的。

表7 多元回歸方程的方差分析

圖1 Chl.a理論計算值與實測值間相關性散點分布圖

2.5 Chl.a與環境因子間的關系模型

經過SPSS作相關分析得到池塘水體中各理化因子間相關性，簡化關系圖見2(a)，具有顯著相關性的兩因子間由雙向箭頭連接，具極顯著相關性的兩因子間由加粗雙向箭頭連接。此圖可較直觀地顯示因子間相關性簡況：各因子間相互作用、相互影響、相互制約；水體氮、磷營養鹽指標間的相關性較為密集，反映水體基本狀況的指標(如，t、pH、S等)間相關關系相對簡單。同時此關系圖也存在明顯缺陷：環境因子間的相關性僅是各因子間互相作用的綜合表觀指示，并非兩因子直接作用的結果，且未能簡明扼要地指出影響葉綠素a的關鍵因子和影響程度。

圖2 Chl.a與環境因子間的關系模型及優化

3 討論

3.1 Chl.a與常規理化因子的關系

眾多研究認為水溫是Chl.a含量變化的關鍵因子[18-19]，水溫可通過控制藻類細胞光合作用以及呼吸代謝速率影響Chl.a含量，兩者呈正相關關系，但這些研究的主體多集中在受人為影響較小的自然水域[20-21]，而在人為措施影響較多的養殖池塘中水溫對藻類生長的影響程度相對被弱化[11]。當自然水體中浮游植物大量繁殖時，水體中二氧化碳和氧氣被大量消耗和產生，水體pH和DO含量升高，Chl.a與pH和DO應呈正相關關系[8]，但大量研究發現很多自然水體以及養殖水體中Chl.a含量的變化與pH和DO的關系并不是如此簡單，畢京博等[18]對太湖水體的研究以及陳金玲等[22]對精養蝦池水體的研究結果均顯示Chl.a含量的變化與pH和DO呈負相關關系，與本文結果一致。分析其原因主要為養殖水體中較高含量的懸浮有機物分解以及生物呼吸導致O2的大量消耗，同時產生大量CO2， 掩蓋了藻類光合作用對pH和DO的影響。

3.2 營養關鍵因子

[1]鮑建平,許兆明.太湖浮游植物初級生產力及魚產力估算[J].淡水漁業,1984,14(5):1-5.

[2]鄧金釵,熊邦喜,張小敏,等.池塘主養草魚三種養殖模式浮游植物群落結構的比較[J].淡水漁業,2012,42(5):63-67.

[3]賴子尼,余煜棉,龐世勛,等.鱖魚養殖池塘水體葉綠素a與16項水生態因子的關系[J].中國水產科學,2004,11(5):426-431.

[4]林 超,韓翠敏,游文華,等.不同水體營養條件對刺苦草和密刺苦草生長的影響[J].生態學雜志,2016,35(8):2117-2121.

[5]Pepe M,Giardino C，Borsani G,et al.Relationship between apparent optical properties and photosynthetic in the sub-alpine Lake Iseo [J].Sci Tot Environ,2001,268(1):31-45.

[6]許 濤,王雨春,劉德富,等.三峽水庫香溪河庫灣夏季水華調查[J].生態學雜志,2014,33(3):646-652.

[7]Romo S,Donk E V,Gylstra R,et al.A multivariate analysis of phytoplankton and food web changes in a shallow biomanipulated lake [J].Freshwat Biol,1996,36(3):683-696.

[8]阮曉紅,石曉丹,趙振華,等.蘇州平原河網區淺水湖泊葉綠素a與環境因子的相關關系[J].湖泊科學,2008,20(5):556-562.

[9]江 敏,盧 柳,邢 斌,等.滴水湖水體葉綠素a與水質因子的多元分析[J].上海海洋大學學報,2011,20(3):451-456.

[10]周貝貝,王國祥,徐 瑤,等.南京秦淮河葉綠素a空間分布及其與環境因子的關系[J].湖泊科學,2012,24(2):267-272.

[11]江 敏,余根鼎,戴習林,等.凡納濱對蝦養殖池塘葉綠素a與水質因子的多元回歸分析[J].水產學報,2010,34(11):1712-1718.

[12]梅衛平,江 敏,阮慧慧,等.滴水湖葉綠素a時空分布及其與水質因子的關系[J],生態學雜志,2013,32(5):1249-1254.

[13]李世凱,江 敏,潘若雷,等.羅氏沼蝦養殖塘葉綠素a與水質因子的灰色關聯分析[J].生態與農村環境學報,2014,30(2):257-261.

[14]韋 璐,江 敏,余根鼎,等.凡納濱對蝦養殖塘葉綠素a與水質因子主成分多元線性回歸分析[J].中國水產科學,2012,19(4):620-625.

[15]王 嵐,王亞平,許春雪,等.多元校正分析在水環境地球化學領域中的應用[J].地質通報,2009,28(2/3):257-262.

[16]Noori R,Sabahi M S,Karbassi A R,et al.Multivariate statistical analysis of surface water quality based on correlations and variations in the data set [J].Desalination,2010,260(1-3):129-136.

[17]章宗涉,黃祥飛.淡水浮游生物研究方法[M].北京:科學出版社.1991.

[18]畢京博,鄭 俊,沈玉鳳,等.南太湖入湖口葉綠素a時空變化及其與環境因子的關系[J].水生態學雜志,2012，33(6):7-13.

[19]李瑞嬌,呂元蛟,張 念,等.主養草魚池塘的浮游植物群落結構及其環境影響因素研究[J].淡水漁業,2014,44(2):77-82.

[20]劉輝宇,宋立榮,萬 能,等.滇池水體葉綠素a與相關環境因子的多元分析[J].江漢大學學報(自然科學版),2007,35(4):87-90.

[21]朱佳志,劉明典,黃福江,等.西雙版納南蠟河浮游植物群落結構及其與水環境因子的關系[J].淡水漁業,2015,45(6):39-44,69.

[22]陳金玲,賴秋明,蘇樹葉，等.精養蝦池主要生態因子變化特點與相關性分析[J].南方水產科學,2012,8(4):49-56.

[23]世界經濟合作與發展組織.水體富營養化監測評價與防治[M].北京:中國環境科學出版社.1989.

[24]Redfield A C.The biological control of chemical factors in the environment [J].Sci Prog,1958,46(3):205-221.

[25]Vollenweider R A.Elemental and biochemical composition of plankton biomass;some comments and explorations [J].Arch Hydrobiol,1985,105(1):11-29.

[26]豐茂武,吳云海,馮仕訓,等.不同氮磷比對藻類生長的影響[J].生態環境,2008，17(5):1759-1763.

[27]Michael F P,Julianne D,Pia H M,et al.Effects of modified nutrient concentrations and ratios on the structure and function of the native phytoplankton community in the Neuse River Estuary,North Carolina,USA.[J].Aquatic Ecology,2002,36(3):371-385.

[28]黃 偉,朱旭宇,曾江寧,等.氮磷比對東海浮游植物群落生長影響的微宇宙實驗[J].環境科學,2012,33(6):18323-1838.

[29]Tilman D,Kilham SS,Kilham P.Phytoplankton community ecology:The role of limiting nutrients [J].Annu Rev Ecol Syst,1982,13(1):349-372.

[30]許 海,朱廣偉,秦伯強,等.氮磷比對水華藍藻優勢形成的影響[J].中國環境科學,2011,31(10):1676-1683.