晁 敏,王云龍,沈新強
(中國水產科學研究院 東海水產研究所 農業部海洋和河口漁業重點開放實驗室,上海 200090)
近年建在河口或沿海地區的濱海電廠不斷增加,電廠運營在保障海岸帶地區經濟發展的同時,其冷卻水系統的取排水對海區生態環境和資源生物也帶來了較大的負面影響[1],比較直觀的影響為急性熱沖擊和卷載效應[2-4],而經冷卻系統加溫的海水排入海域后同樣具有環境效應[5-6],據研究,電機容量為1800~2000MW 機組的溫排水可使出水口水溫升高8~12℃,而為防止污損生物對冷卻系統的阻塞,氯制劑廣泛應用于電廠冷卻水系統,導致排水口氯的殘留濃度可達 0.1~0.2 mg/L[7],有時可達 0.42~1.0 mg/L[8],因此電廠溫排水同時具有溫度升高和余氯殘留兩個特點。
從20世紀50年代以來,國外研究者就對海水中殘留余氯的生態效應開展了大量研究[9-12],中國自80年代以來也開展了相關影響研究[13-15]。國際上已有余氯標準頒布,如美國 EPA水質標準規定水中總余氯的標準限值為 0.01 mg/L,而歐洲標準更低,為0.005 mg/L,國內尹伊偉[14]早在1992年就曾撰文推薦0.01 mg/L的標準限值,然而至今中國還未制訂溫排水中余氯排放限值標準。
魚類在胚胎發育及仔魚期容易受到水溫[16-17]及污染物[18-20]的影響,而中國目前關于水溫和余氯雙因素對海水魚類胚胎發育和仔魚期影響的研究還較少[21],本研究以黑棘鯛(Acanthopagrus schlegelii)(舊稱黑鯛(Sparus macrocephalus))的前期仔魚為對象,分析溫度、余氯及兩種因子交互作用對仔魚存活率的影響,黑棘鯛具有適溫適鹽性廣、食譜雜、抗病能力較強、生長迅速等生理和生態特點,在中國東南沿海屬經濟魚種,因此是一種比較理想的試驗動物[22],本研究可望為中國將來制訂沿海電廠溫排水溫度和殘留余氯濃度限值提供基礎數據支撐。
試驗在江蘇省海洋水產研究所啟東海水增養殖基地進行。試驗用海水取自基地,經0.45 μm醋酸纖維素濾膜過濾,煮沸、曝氣24 h后備用。取安替福民即次氯酸鈉(分析純)配制余氯母液,由于該母液中氯屬游離態氯,在光作用下易于降解,需避光保存。
試驗設4個水溫條件:18、22,26、30℃,其中18℃為室溫條件下水溫,后 3個水溫條件由水浴控溫; 每個水溫條件下設 0.025、0.05、0.10、0.20、0.40、0.80 mg/L總計6個余氯濃度水平,并以過濾海水為試驗對照,每個濃度設3個平行試驗組。試驗開始前各組預先調節好水溫,加入余氯母液。
以剛孵化出1 d的黑棘鯛前期仔魚為試驗生物,用吸管挑取前期仔魚20~30只置于已調節好水溫、余氯濃度的400 mL燒杯內,并計仔魚放入的時刻為0,于24、48、72、96 h用解剖鏡觀察仔魚存活狀態,當仔魚心臟停止跳動時判定為死亡,試驗期間投喂輪蟲為仔魚的開口餌料。
由于中國電廠冷卻系統加氯方法通常采取周期性連續加氯或每天加氯1 h的方式[23],電廠溫排水中余氯濃度并不是恒定不變的,而且余氯受光照、溫度等因素影響在短短幾小時內可衰減至極低水平,因此為模擬電廠溫排水中余氯分布特征,試驗期間各組試驗液每12 h全部更換,每次海水經水溫預調、加氯,再將仔魚轉移至新的燒杯中。
分別計算24、48、72、96 h仔魚存活率,對仔魚存活率數據進行標準化處理,處理方法如下:

其中Ps為黑棘鯛仔魚標準化死亡率(以100計),P為各組黑棘鯛仔魚存活率(以100計)。
余氯對黑棘鯛仔魚半致死余氯濃度(LC50)計算方法參見文獻[24],采用下式在 Statistica6.0軟件中對各組數據進行曲線擬合:

其中x為余氯濃度,T、W兩個參數由迭代得到,其中T為LC50或EC50,W代表曲線的曲率因為:

當Ps為50時,x=T。
采用 SPSS統計軟件比較各組死亡率是否具有顯著性差異分析,當方差齊時,采用多重比較方差分析方法(LSD),當P< 0.05時統計具有顯著性差異意義,當方差不齊時,對數據進行排秩處理,采用SPSS中廣義線性模型(Generalized linear model,GLM)的多元方差分析中的 SNK (Student-newman-kuels)方法[25],當P< 0.05時統計具有顯著性差異意義。
采用廣義相加模型[26](Generalized additive model,GAM)研究仔魚死亡率與溫度、余氯、曝露時間的關系,模型表達形式如下:

其中:g為仔魚死亡率的連接函數,S為光滑函數,一般采用 3次樣條函數[27-28],其中S1為水溫的光滑函數,S2為余氯的光滑函數,S3為曝露時間的光滑函數,error為隨機誤差。
仔魚死亡數據為二元屬性數據,因此死亡率可采用邏輯連接函數表示,其數學表達式為:

因此廣義相加模型可寫為:

其中,error為隨機誤差。
然而,為更有效地分析各因素間的交互作用,在SPSS中采用雙因素方差分析方法進行研究,當P< 0.05時表明對死亡率具有顯著性影響。
用于試驗的仔魚剛孵化出1 d,仍處于前期仔魚階段,卵黃囊尚未消失,仔魚為內源營養期。6 h時,各組沒有仔魚死亡,其中,18℃組、22℃組仔魚可靈活游動,26℃組的仔魚則沉于燒杯底部,但對接觸反應靈敏,而 30℃組仔魚的背鰭基部出現隆起,發育速度明顯高于其他組,但已有部分仔魚出現尾部彎曲,且對接觸沒有反應,表明運動受到抑制。24 h時,18℃組、22℃組、26℃組沒有仔魚出現死亡,其中,18℃組、22℃組的尾鰭、消化道已形成,26℃組除尾鰭、消化道形成外,部分仔魚已開口,而30℃組全部仔魚運動受到抑制,有 23.4%的仔魚死亡。48 h時,只有18℃組未出現仔魚死亡,22℃組、26℃組、30℃組出現了不同程度的死亡情況,死亡率分別為2.9%、1.3%、98%,30℃組仔魚死亡率增加明顯。72 h時,各組仔魚均有死亡,死亡率分別為0.8%、2.9%、7.7%、98%。96 h時,各組死亡率分別為1.6%、7.0%、7.7%、98%,22℃組、26℃組仔魚死亡率沒有變化,表明該實驗組仔魚的生長和發育已非常適應試驗水溫條件(表 1)。
對各實驗組數據進行多重比較分析,18℃組、22℃組、26℃組 3組仔魚死亡率沒有顯著差別(P>0.05),且死亡率均低于10%,可認為這3個水溫條件對仔魚死亡不存在顯著性影響,而 30℃水溫則對仔魚死亡具有顯著性影響(P< 0.05)。
每種生物對水溫均有一定的耐受范圍,仔魚的耐熱性往往低于幼魚或成魚,曾江寧[29]的研究表明,初孵 1周的黑棘鯛仔魚 24 h的半數致死水溫為30.3℃,與本實驗結果相近,24 h時30℃水溫造成的初孵仔魚死亡率為23.4%。

表1 對照組黑棘鯛仔魚不同水溫條件下的死亡率Tab.1 The lethal rate of Acanthopagrus schlegelii larvae in control group at different water temperature
高濃度余氯可在短時間內造成仔魚死亡,仔魚體形一般為 S形扭曲,或尾部斷裂、甚至蜷縮成球狀。由表2可知,0.8 mg/L的余氯在24 h時對仔魚的死亡率與對照組具有顯著性差異(P<0.05),隨著曝露時間的延長,低濃度余氯對仔魚的死亡率開始與對照組顯現出差異,水溫越高,余氯致毒濃度愈低,如26℃組在48 h時0.05 mg/L的余氯對仔魚的死亡率與對照組具有顯著性差異。
采用公式 1對各組數據進行擬合,得到余氯對仔魚的半致死效應濃度LC50值,列于表3。大多數生物在溫度升高時表現出對更低劑量的余氯敏感,如大型溞(Daphnia magna)、白鰱(Hypophthalmichthys molitrix)[14],虹鱒(Oncorhynchus mykiss)、美洲原銀漢魚(Menidia beryllina)[30],銀鮭(Oncorhynchus kisutch)、胡瓜魚(Osmerus mordax)、斑尾美洲(Notropis hudsonius)[31]等,黑棘鯛仔魚也不例外,從表 3總數據可明顯看出,余氯對黑棘鯛仔魚的 LC50值顯著受到了水溫的協同影響,當水溫達到 30℃時,水溫則成為仔魚主要死亡因素。

表2 各余氯濃度組仔魚死亡率與對照組仔魚死亡率差異性檢驗Tab.2 Significance test of mortality rate of Acanthopagrus schlegelii larvae between rc (residual chlorine) and control group

表3 各水溫條件下余氯對仔魚的LC50值(mg/L)Tab.3 Lethal concentration of residual chlorine to the larvae of Acanthopagrus schlegelii at different temperature
采用 GAM 模型研究仔魚死亡率與水溫、余氯、暴露時間的關系,建立 GAM 模型連接函數所需數據見表4,對response數據,以二元數據賦值,對仔魚未死亡賦值 0,其對應余氯的值選擇余氯對仔魚無觀測影響濃度(None observed effect concentration,簡稱NOEC),即與對照無顯著性差異的余氯濃度,見表 2;對50%死亡的仔魚賦值0.5,其對應余氯的值為LC50。

表4 GAM模型連接函數所需數據Tab.4 Data for link function of GAM

圖1 仔魚死亡率與水溫、余氯、曝露時間的光滑函數擬合結果Fig.1 GAM plots for temperature,residual chlorine and exposure time
GAM分析結果(表5)顯示,仔魚死亡率與水溫、余氯含量、曝露時間并不存在 3次或更復雜分段樣曲線關系(Pr > 0.05),而是近似的分段線性關系,見圖1-1,各因子的光滑函數隨對應因子的變化曲線也具有此特點,從圖1-2可看出,水溫在18~26℃間與S(水溫,4)間具有線性關系,從圖1-3可看出,余氯在0.025~0.4 mg/L 與S(余氯,4)間具有近似線性關系,從圖1-4可看出,96 h曝露時間內時間與仔魚死亡率具有近似線性關系。各因子的光滑函數隨對應因子的變化曲線可進一步解釋仔魚死亡率隨因子變化的關系,當光滑函數值高于0時,所對應的因子值即對仔魚產生明顯致死效應的數值,從圖1-2看出,水溫為 26℃時,S(水溫,4)開始大于 0,水溫由 26℃升至30℃時,S(水溫,4)急劇上升,表明水溫高于 26℃對黑棘鯛初孵仔魚的致死效應明顯,從圖 1-3可看出,余氯濃度高于 0.2 mg/L時,開始對仔魚具有明顯的致死效應,從圖 1-4可看出,在有余氯存在情況下,曝露時間超過 48 h,會對黑棘鯛初孵仔魚產生明顯致死效應。
經多元方差分析結果可知(表 6),水溫、余氯、曝露時間對仔魚死亡率具有顯著性影響(P< 0.05),水溫-余氯、水溫-曝露時間之間具有顯著性交互作用(P<0.05),從圖 2-1和圖 2-2存在線段相交的情況也可說明交互作用的存在,其中,水溫和余氯的交互作用主要發生在低濃度余氯組,余氯和曝露時間的交互作用主要發生在24 h和48 h; 而水溫與曝露時間之間沒有顯著的交互作用(P=0.989>0.05),從圖2-3可看出沒有線段相交的情況。
與多種化學因子如重金屬Cu[32]、Cd[33]、農藥[34]與水溫具有交互作用一樣,余氯與水溫之間也具有交互作用,且一般為協同作用[32],而從本實驗結果看,當水溫高于 30℃或余氯濃度高于 0.8 mg/L時,水溫和余氯間均可成為仔魚死亡主導因素,二者間的交互作用很難被觀測到。
曝露時間與污染物的毒性間通常具有顯著的時間-劑量效應關系,本實驗中,仔魚死亡率與曝露時間之間呈線性關系,隨著曝露時間的延長,余氯的LC50值逐漸減少,與Hall等[35]研究結論相似。

表5 廣義相加模型擬合結果Tab.5 Fitted results of GAM

表6 雙因素方差分析受試者間效應結果Tab.6 Results of tests between subjects effects of univariate analysis

圖2 各因子間交互作用Fig.2 The plot of interaction between factors
18、22、26℃ 3種水溫對黑棘鯛初孵仔魚致死率均低于 10%,與對照組相比不存在顯著性差異,30℃水溫對仔魚死亡率與對照相比具有顯著差異(P< 0.05),24 h時初孵仔魚死亡率 23.4%,48 h時為98%。
24 h時18、22、26、30℃ 4種水溫條件下余氯對黑棘鯛仔魚的的 LC50值分別為 0.816、0.460、0.433、0.319 mg/L,可見水溫升高會增強余氯對仔魚的致死效應,而隨著時間延長,余氯對仔魚的 LC50效應濃度值進一步降低,至96 h,其值分別為0.242、0.211、0.140、<0.025 mg/L。
GAM分析結果顯示,水溫在18~26℃間與仔魚死亡率間具有線性關系,水溫由 26℃升至 30℃時,仔魚死亡率有急劇上升,表明水溫高于 26℃對黑棘鯛初孵仔魚致死效應明顯; 余氯在 0.025~0.4 mg/L與仔魚死亡率間具有近似線性關系,當余氯濃度高于0.2 mg/L時,對仔魚致死效應明顯; 96 h內曝露時間與仔魚死亡率具有近似線性關系,在有余氯存在情況下,曝露時間超過 48 h,會對黑棘鯛初孵仔魚產生明顯致死效應。
經多元方差分析結果可知,水溫、余氯、曝露時間對仔魚死亡率具有顯著性影響(P<0.05),水溫-余氯、水溫-曝露時間之間具有顯著性交互作用(P<0.05)。
致謝:本次實驗得到了江蘇省海洋水產研究所啟東海水增養殖基地的大力支持,同時也得到了東海水產研究所蔣玫副研究員、黃世林、馬繼臻同學的大力幫助,謹此表示感謝。
[1]曾江寧,陳全震,鄭平,等.余氯對水生生物的影響[J].生態學報,2005,25(10):2717-2724.
[2]Bamber R N.The influence of rising background temperature on the effects of marine thermal effluents[J].Journal of Thermal Biology,1995,20(1-2):105-110.
[3]苗綠田.核電站冷卻水系卷載效應的初步研究[J].海洋湖沼通報,1989,1:13-19.
[4]盛連喜,侯文禮.電廠冷卻系統對梭幼魚和對蝦仔蝦卷載效應的初步探討[J].環境科學學報,1994,14(1):47-55.
[5]黃洪輝,王肇鼎,張穗.電廠溫排水中的余氯對鄰近水域生態環境的影響[J].南海研究與開發,1998,2:46-58.
[6]張穗,黃洪輝,陳浩如,等.大亞灣核電站余氯排放對鄰近海域環境的影響[J].海洋環境科學,2000,19(2):14-18.
[7]Bamber R N,Seaby R M H.The effects of power station entrainment passage on three species of marine planktonic crustacean,Acartia tonsa(Copepoda),Crangon crangon(Decapoda) andHomarus gammarus(Decapoda)[J].Marine Environmental Research,2004,57(4):281-294.
[8]張燕,孫英蘭.電廠溫排水中余氯濃度預測[J].海洋科學,2007,31(2):5-8.
[9]Buckley J A.Acute toxicity of residual chlorine in wastewater to coho salmon (Oncorhynchus kisutch) and some resultant hematologic changes[J].Journal of the Fisheries Research Board of Canada,1976,33(12):2854-2856.
[10]Zeitoun I H.Assessment of intermittently chlorinated heated effluents on survival of adult rainbow trout(Salmo gairdneri) at power generating facilities[J].Environmental Science & Technology,1978,12(10):1173-1179.
[11]Giattlna J D,Cherry D S,Cairnsjr J,et a1.Comparison of laboratory and field avoidance behavior of fish in heated chlorinated water[J].Transactions of the American Fisheries Society,1981,110(4):526-535.
[12]Hall L W,Burton D T,Margrey S L,et a1.The effect of acclimation temperature on the interactions of chlorine,△T and exposure duration to eggs,prolarvae and larvae of striped bass,Morone saxatilis[J].Water Research,1983,17(3):309-317.
[13]林可椒,高宏偉.余氯對鱖魚苗的急性毒性試驗[J].淡水漁業,1991,3:24-25.
[14]尹伊偉,溫曉燕,呂頌輝,等.不同溫度下氯對水生生物的毒性及制訂漁業水質標準的探討[J].生態科學,1992,2:41-49.
[15]江志兵,曾江寧,陳全震,等.濱海電廠冷卻水余熱和余氯對中華哲水蚤的影響[J].應用生態學報,2008,19(6):1401-1406.
[16]張海發,劉曉春,王云新,等.溫度、鹽度及pH對斜帶石斑魚受精卵孵化和仔魚活力的影響[J].熱帶海洋學報,2006,25(2):31-36.
[17]關健,柳學周,蘭春燕,等.溫度、鹽度對褐牙鲆(♀)×犬齒牙鲆(♂) 雜交子一代胚胎發育和仔魚存活的影響[J].海洋水產研究,2007,28(3):31-37.
[18]陳民山,范貴旗.勝利原油對海洋魚類胚胎及仔魚的毒性效應[J].海洋環境科學,1991,10(2):1-5.
[19]蔣玫,沈新強,晁敏,等.AOX 漂白廢水對黑鯛魚卵及仔魚的毒性效應[J].環境科學研究,2007,19(2):27-30.
[20]蔣玫,林欽,沈新強.航道疏浚泥懸浮液對黑鯛胚胎的發育及仔魚攝食生長的影響[J].生態環境學報,2009,18(5):1674-1678.
[21]江志兵,廖一波,高愛根.余氯對魚類毒性影響的研究進展[J].海洋學研究,2009,27(4):86-94.
[22]穆景利,王新紅,王淑紅,等.苯并(a)芘影響黑鯛肝臟EROD活性變化的動力學研究[J].海洋科學,2009,33(1):68-71.
[23]劉蘭芬,郝紅,呂光四.電廠溫排水中余氯衰減規律及其影響因素的影響研究[J].水利學報,2004,5:94-98.
[24]Hampel M,Moreno-Garrido I,Sobrino C,et al.Acute toxicity of LAS homologues in marine microalgae:esterase activity and inhibition growth as endpoints of toxicity[J].Ecotoxicology and Environmental Safety,2001,48:287-292.
[25]劉萬里,薛茜,曹明芹,等.用 SPSS實現完全隨機多組比較秩和檢驗的多重比較[J].地方病通報,2007,22(2):27-29.
[26]Taylor J M G,Zhang F,Withers H R.Non-linear effects in dose-time analyses:application of modern statistical techniques[J].Radiotherapy and Oncology,1997,45:133-140.
[27]Cawsey E M,Austin M P,Baker B L.Regional vegetation mapping in Australia:a case study in the practical use of statistical modeling[J].Biodiversity and Conservation,2002,11:2239-2274.
[28]赫曉慧,溫仲明,王金鑫.基于 GAM 模型的延河流域主要草地物種空間分布及其與環境的關系[J].生態學雜志,2008,27(10):1718-1724.
[29]曾江寧.濱海電廠溫排水對亞熱帶海域生態影響的研究[D].杭州:浙江大學,2008.1-157.
[30]Fisher D J,Burton D T,Yonkos L T.The relative acute toxicity of continuous and intermittent exposure of chlorine and bromine to aquatic organisms in the presence and absence of ammonia[J].Water Research,1999,33(3):760-768.
[31]Seegert G L,Brooks A S.The effects of intermittent chlorination on coho salmon,alewife,spottail shiner,and rainbow smelt[J].Transactions of the American Fisheries Society,1978,107(2):346-353.
[32]Gama-Flores J L,Sarma S S S,Nandini S.Interaction among copper toxicity,temperature and salinity on the population dynamicsof Brachionus rotundiformis(Rotifera)[J].Hydrobiologia,2005,181:559-568.
[33]Abdel-lateif H M,Donker M H,Van Straalen N M.Interaction between temperature and cadmium toxicity in the isopodPorcellio scaber[J].Functional Ecology,1998,12(4):521-527.
[34]Fisher S W.Changes in the toxicity of three pesticides as a function of environmental pH and temperature[J].Bulletin of Environmental Contamination and Toxicology,1991,46(2):197-202.
[35]Hall L W,Burton D T,Margrey S.The influence of acclimation temperature on the interactions of chlorine,elevated temperature,and exposure duration for grass shrimp,Palaemonetes pugio[J].Transactions of the American Fisheries Society,1979,108:626-631.