王龍燕,閔文武,王金樂,陳飛雄,3,周其椿,吳俁學
(1 貴州省生物研究所,貴州 貴陽 550009;2 貴州省農業科學院水產研究所,貴州 貴陽 550025;3 貴州省特種水產工程技術中心,貴州 貴陽 550025)
集裝箱循環水養殖系統(CRAS)是以集裝箱為養殖載體,應用高新技術有效控制養殖環境和養殖過程,實現分區養殖、易位處理和循環利用,具有節地、節水和高產等優點,是促進漁業生產高質量發展的有效技術模式[1]。應用該系統養殖羅非魚(Oreochromismossambicus)、草魚(Ctenopharyngodonidella)、大鱗鲃(Barbuscapito)、黃顙魚(Pelteobagrusfulvidraco)、鯉魚(Cyprinuscarpio)等品種取得了良好的經濟和社會效益[2-6]。但集約化、高密度的養殖必然帶來污染物的大量排放,識別系統排放尾水中的主要污染物并提高污染物的處理效率成為當前研究的熱點[7-9]。“集裝箱+生態池塘”是集養殖生產和尾水處理的綜合養殖系統,具有雙高效的優勢,是農業農村部2018—2020年主推的養殖模式之一[10]。
浮游植物是水生態系統中主要生物群落之一,是構成水生態系統中多元化營養鏈的第一步,是浮游動物和底棲動物等水生動物的主要攝食對象[11-12],是水環境健康狀況的重要指示生物[13-14]。在池塘循環流水養殖和人工濕地-池塘復合養殖系統等綜合養殖模式中關于浮游植物的研究已有報道[15-18],“集裝箱+生態池塘”循環水養殖系統中浮游植物的研究尚未見報道。
通過對系統內水質指標和浮游植物進行檢測,分析系統尾水污染狀況,探索浮游植物變動規律,為科學管理和調節系統水質質量,提高系統養殖效益提供理論指導,為推廣應用集裝箱循環水養殖新模式提供科學依據和數據支撐。
試驗基地集裝箱(6.1 m×2.8 m×2.4 m)共20個,養殖水體容積為25 m3/個,養殖魚類品種為鯉魚、草魚、大口黑鱸(Micropterussalmoides)、黃顙魚、斑點叉尾鮰(IetalurusPunetaus)和長吻鮠(Leiocassislongirostris)。系統水流方向為集裝箱養殖區→物理過濾區→一級沉淀池→二級沉淀池→生態凈化池塘→集裝箱養殖區的循環(圖1)。系統每天早上6時開始每隔3 h循環處理一次,每次1 h,每次通過2臺轉鼓微濾機輪流進行過濾,每天循環6次把養殖水體全部更換,每次循環回用水量的70%為生態凈化池塘水體, 30%為水庫新水。

圖1 “集裝箱+生態池塘”循環水養殖系統水流方向示意圖
物理過濾區配置2臺120目轉鼓微濾機,單臺水處理量為80 m3/h。生態池塘由一級沉淀池、二級沉淀池和生態凈化塘組成,一級沉淀池(18 m2)由3個逐級降低的階梯式水池(長×寬×深:6 m×1 m×1 m)連接組成,一級沉淀池尾水逐級漫水進入二級沉淀池,二級沉淀池(長×寬×深:6 m×4 m×1.5 m,24 m2)由鏤空磚墻圍砌而成,二級沉淀池與生態凈化塘頂高程相同,水通過磚墻縫隙進入生態凈化塘,生態凈化塘面積4 000 m2,平均水深2.0 m。一級沉淀池種植有粉綠狐尾藻(Myriophyllumaquaticum),布滿沉淀池;二級沉淀池采用浮筏種植有粉綠狐尾藻和豆瓣菜(Nasturtiumofficinale),分別占二級沉淀池面積的80%和5%,一級和二級沉淀池中水生植物在試驗期內未進行采收。生態凈化塘內放養鰱(Hypophthalmichthysmolitrix)和鳙(Hypophthalmichthysnobilis),放養密度為鰱(300 g/尾)100尾/畝(1畝=0.067 hm2),鳙(200 g/尾)100尾/畝,試驗期內未見死亡,也未進行捕撈。

運用浮游植物物種數、密度、生物量、Shannon-wiener 多樣性指數和Margalef豐富度指數進行浮游植物多樣性及群落結構特征分析,運用污染物去除率(E)進行尾水治理效果分析[24],采用單因子污染物指數法(Pij)和總污染指數(S)識別尾水中主要污染物和污染程度[25],浮游植物優勢物種根據優勢度指數(Y)判定[26],Y>0.02即為優勢種,Y、E、Pij和S計算公式如下:
Y=Fi×Pi
(1)
(2)
(3)
S=∑pij
(4)
式中:Y為優勢度指數,Fi為第i個物種出現頻率(出現樣點數/總采樣點數);Pi為第i個物種相對豐度(i物種豐度/總物種豐度);E為污染物去除率,%;CA和CB分別為系統A和B點污染物質量濃度(mg/L)。Pij為污染物i的標準指數,Cij為污染物i在監測點j的實測質量濃度(mg/L);Csi為水質參數i的地表水水質標準,根據河段水功能區劃分,本研究依據國家標準[27]Ⅲ類水質標準進行計算。S為總污染物指數。
運用PRIMER 5計算浮游植物Shannon-wiener 多樣性指數和Margalef豐富度指數,運用Excel office 2019和SPSS 19.0進行數據統計處理,A、B兩點水質及浮游植物群落結構差異采用方差分析和t檢驗。運用CANOCO 4.5 分析浮游植物與環境因子之間的關系,先用浮游植物密度進行去趨勢對應分析(Detrended correspondence analysis,DCA),結果中Lengths of gradient第一軸值大于4時,選擇典范對應分析(Canonical Correspondence Analysis,CCA),介于3~4之間選用CCA和線性模型的冗余分析(Redundancy analysis,RDA)均可,小于3時選擇RDA。選用Y>0.000 1的物種進行DCA、CCA和RDA分析,分析前對生物數據和環境數據(除pH外)進行lg(x+1)轉換,運用前選法(Forward selection)和999次蒙特卡洛置換檢驗(Monte Carlo permutation test)識別各環境因子的貢獻率。
系統T、pH和DO變化情況如表1所示, A和B的T整體表現為從1月至8月不斷升高并達到最高值,9月至12月不斷降低,夏季和冬季水溫差異顯著(P<0.05),A點T比B點高0.6℃,整體增加4.23%,但差異不顯著。pH和DO月變化差異不顯著,尾水通過凈化后pH和DO均有不同程度升高,分別升高5.17%和23.49%。

表1 系統T、pH和DO月變化情況

圖2 A和B污染物質量濃度變化及去除率


表2 系統和TP污染指數