朱峰,石志洲,凌信文,夏永健,李陽,翁怡嬋,劉溢馨
(1.國家海洋局閩東海洋環(huán)境監(jiān)測中心站,福建 寧德 352100;2.海洋赤潮災(zāi)害立體監(jiān)測技術(shù)與應(yīng)用國家海洋局重點實驗室,上海 200081)
三沙灣位于福建省寧德市東南部,是一個半封閉形的港灣,海域面積714 km2(黃標等,2002)。三沙灣海區(qū)東北部毗鄰霞浦縣的溪南鎮(zhèn),是全國有名的“大黃魚之鄉(xiāng)”,中部有官井洋野生大黃魚的產(chǎn)卵場和三都島、青山島兩個主要的網(wǎng)箱養(yǎng)殖大黃魚基地,西北部的福安和蕉城區(qū)沿岸也有大量的網(wǎng)箱養(yǎng)殖區(qū),整個海區(qū)內(nèi)有22 萬個網(wǎng)箱,網(wǎng)箱養(yǎng)殖的品種主要為大黃魚、鱸魚、鯛科魚類等,其中大黃魚養(yǎng)殖占60%左右(蘇躍中,2009)。其中在青山島周圍不到10 hm2的海域上分布著約10 萬個網(wǎng)箱,網(wǎng)箱的規(guī)格一般為4 m×4 m×7 m,育苗的養(yǎng)殖密度為每箱1×104尾,成魚的養(yǎng)殖密度約為每箱1 500 尾(廖紅梅等,2011)。近年來,隨著水產(chǎn)養(yǎng)殖業(yè)的迅速發(fā)展,三沙灣海區(qū)內(nèi)的網(wǎng)箱數(shù)量逐年增加(宋蘭芬,2007)。在夏秋季節(jié),富營養(yǎng)化的水體中溶解氧會被大量消耗,各種細菌病和寄生蟲病的病原體(如弧菌類、刺激隱核蟲、本尼登蟲等) 也會在網(wǎng)箱內(nèi)大量滋生(林克冰 等,1999),若不及時采取措施改善網(wǎng)箱內(nèi)的水質(zhì)環(huán)境,極易造成養(yǎng)殖魚類的突發(fā)性疾病并大批死亡,例如2009年5月至8月,福建省寧德市網(wǎng)箱養(yǎng)殖大黃魚相繼發(fā)生刺激隱核蟲病,造成成魚損失近萬噸,魚苗損失1 億多尾,經(jīng)濟損失約1.87 億元。本文根據(jù)2011年5月和8月的監(jiān)測資料,對三沙灣網(wǎng)箱養(yǎng)殖區(qū)和非網(wǎng)箱養(yǎng)殖區(qū)的水質(zhì)現(xiàn)狀進行分析對比,研究了水質(zhì)環(huán)境對網(wǎng)箱養(yǎng)殖病害的影響;對水質(zhì)條件和魚病發(fā)生時的死亡率進行相關(guān)性分析和研究,得到三沙灣網(wǎng)箱養(yǎng)殖病害和水文條件的相關(guān)性;分析了網(wǎng)箱養(yǎng)殖對底棲生物群落結(jié)構(gòu)的影響。通過本研究,力圖為網(wǎng)箱養(yǎng)殖病害的預(yù)防和防治提供基礎(chǔ)數(shù)據(jù),并為網(wǎng)箱養(yǎng)殖的可持續(xù)發(fā)展以及海區(qū)的環(huán)境保護提供科學依據(jù)。
在三沙灣海區(qū)內(nèi)共設(shè)置了9 個監(jiān)測站點,其中4 個位于網(wǎng)箱養(yǎng)殖區(qū)(3、4、6、8),4 個位于非網(wǎng)箱養(yǎng)殖區(qū)(2、5、7、9),1 號點安置了哈希hydrolab DS5X 多功能水質(zhì)自動監(jiān)測探頭,位于青山島和斗帽島之間網(wǎng)箱養(yǎng)殖最密集的地區(qū),其站點分布詳見圖1。三沙灣的北面有閩東第一大江賽江又叫白馬河入海,西北面有“五江三溪”之一的霍童溪注入三沙灣,兩條河流帶來的淡水對2、3、4、6 站點的鹽度有顯著的影響。
圖1 2011年三沙灣采樣站點分布圖
2011年5月16-19日和8月18-20日(即網(wǎng)箱養(yǎng)殖大黃魚正常期和細菌病暴發(fā)期) 分別對各個監(jiān)測站點采樣一次,采樣水深為0.5 m。8 個站點共同的監(jiān)測項目有水溫、鹽度、透明度、pH、葉綠素-a、溶解氧(DO)、化學需氧量(COD)、亞硝酸鹽(NO2-N)、硝酸鹽(NO3-N)、銨鹽(NH4-N)、活性磷酸鹽(PO4-P) 和弧菌總數(shù)。
在細菌病暴發(fā)、病魚大量死亡期間(8月23日-9月10日),利用經(jīng)過校準后的哈希hydrolab DS5X 多功能水質(zhì)自動監(jiān)測探頭,于每日14 點觀測溫度、鹽度、pH、葉綠素、溶解氧等水質(zhì)指標,同時在該站點隨機選擇6 個網(wǎng)箱,每日統(tǒng)計1 齡大黃魚的平均死亡量。8月份還有3 個站位(4、7、8) 用采泥器采集底棲生物一次。樣品采集和分析方法按照海洋監(jiān)測規(guī)范(GB 17378-2007,2008)執(zhí)行。
水質(zhì)評價標準執(zhí)行海水水質(zhì)標準第二類標準值(pH 7.5-8.5,COD ≤3 mg/L,DIN ≤0.30 mg/L,PO4-P≤0.030 mg/L) (GB 30972-1997, 1998)。Redfield 提出海水中的平均N/P 是15,而浮游植物N/P 的組成為16∶1,用近海海水進行生物培養(yǎng)發(fā)現(xiàn),當N/P<8 時浮游植物生長受N 限制,N/P>30時受P 限制。水體營養(yǎng)級的標準按文獻(郭衛(wèi)東等,1998) 采用的標準(表1)。
海水富營養(yǎng)程度的評價采用鄒景忠提出的方法(林榮根,1996)。
其中各要素的單位是mg/L,當E≥1,既表示該海區(qū)富營養(yǎng)化,數(shù)值越大既富營養(yǎng)化程度越高。
使用t 檢驗法分析同一月份養(yǎng)殖區(qū)和非養(yǎng)殖區(qū)的差異顯著性,以及不同季節(jié)的同一用海類型的海區(qū)差異顯著性(各因素顯著水平設(shè)置為P<0.05)(GB 17378.2-2007,2008)。
利用SPSS11.0 軟件統(tǒng)計各水質(zhì)因素與死亡數(shù)量的均值、標準差和變異系數(shù),并用雙因素分析法對所得數(shù)據(jù)進行Pearson 相關(guān)性分析(各因素顯著水平均設(shè)置為P<0.05)。
采用PRIMER 6.0 統(tǒng)計軟件包計算3 個站位的底棲生物Shannon-Wiever 物種多樣性指數(shù)(H′)、Margalef 物種豐度指數(shù)(d)、Pielous 物種均勻度指數(shù)(J) 和豐度/生物量比較曲線,即ABC 曲線(Abundance/Biomass Curves)。根據(jù)生態(tài)演替理論,未受擾動的群落,趨向于由少量的大個體、長生活史的物種主導,生物量曲線完全在豐度曲線之上。在中度干擾的群落將會被短生活史和高種群增長率的物種占據(jù),此時生物量曲線和豐度曲線較為接近。在受嚴重擾動的群落,群落由小個體機會性物種占據(jù),豐度曲線將完全位于生物量曲線之上(Vergnon et al,2006)。
表1 海水營養(yǎng)級的劃分標準
細菌病暴發(fā)前后各站點的監(jiān)測結(jié)果見表2。
2.2.1 水質(zhì)狀況比較
通過對5月和8月8 個監(jiān)測站點兩次監(jiān)測數(shù)據(jù)的分析顯示(見表2):5月份網(wǎng)箱內(nèi)的水溫比對照海區(qū)略低,8月份網(wǎng)箱內(nèi)水溫比對照海區(qū)略高;而鹽度網(wǎng)箱區(qū)比對照組都略低,其中5月份有差異(P<0.05)。這是由于網(wǎng)箱一般都會選擇在沿岸水體較淺的地方,水溫變化幅度會比中心海區(qū)快。沿岸河流和淡水的注入,也會使網(wǎng)箱養(yǎng)殖區(qū)的3、4、6號點海水的鹽度降低。網(wǎng)箱內(nèi)的DO 略低,但差異并不顯著(P>0.05)。8月份弧菌數(shù)量大于對照組(P<0.05),這是因為網(wǎng)箱養(yǎng)殖區(qū)養(yǎng)殖魚類排泄物和殘餌等有機養(yǎng)料相對豐富,各種細菌病的病原體也會滋生;網(wǎng)箱內(nèi)無機氮和無機磷的含量均大于對照海區(qū),表現(xiàn)出一定程度的富營養(yǎng)化。
表2 細菌病暴發(fā)前(5月) 和發(fā)病期間(8月) 各站點水質(zhì)調(diào)查表
2.2.2 N/P 值比較
養(yǎng)殖區(qū)的N/P 值最接近Redfield 值,非養(yǎng)殖區(qū)的N/P 值明顯高于16∶1。在兩個月份N/P 值都是網(wǎng)箱區(qū)與對照區(qū)有差異(P<0.05),這是由于5月養(yǎng)殖區(qū)的磷酸鹽含量顯著高于對照區(qū)(P<0.05) 所致,而8月份主要是由于養(yǎng)殖區(qū)的無機氮和磷酸鹽含量都高于對照區(qū)而造成的??梢娋W(wǎng)箱養(yǎng)殖給海區(qū)帶來的主要污染物是活性磷酸鹽和無機氮。
2.2.3 富營養(yǎng)程度比較
海區(qū)普遍富營養(yǎng)化,網(wǎng)箱區(qū)更為嚴重(表2)。養(yǎng)殖區(qū)絕大部分水質(zhì)都是富營養(yǎng),只有夏季8 號點的為中度營養(yǎng)。對照區(qū)水質(zhì)稍好,靠近白馬河口的非養(yǎng)殖區(qū)2 號點為富營養(yǎng),位于三沙灣中央的7 號點為磷中等限制潛在性富營養(yǎng),位于東吾洋的5 號點為中度營養(yǎng),靠近東沖口的9 號點富營養(yǎng)化指數(shù)最低,為中度營養(yǎng)。網(wǎng)箱區(qū)除了夏季的8 號點,均不符合二類水質(zhì),非養(yǎng)殖區(qū)營養(yǎng)鹽濃度稍低,但是除了9 號點,也達不到二類水質(zhì)標準。網(wǎng)箱內(nèi)的無機氮和磷酸鹽都是富營養(yǎng)級水平,對照區(qū)磷酸鹽偏向中度營養(yǎng)并且5月份平均E 值為還未達到臨界值。
2.3.1 水質(zhì)狀況比較
從5月到8月,養(yǎng)殖區(qū)和非養(yǎng)殖區(qū)在溫度升高的同時,DO 下降、葉綠素含量和弧菌數(shù)量升高(P<0.05),這是因為8月份養(yǎng)殖魚類及海洋生物生長旺盛,雖然浮游植物的光合作用會產(chǎn)生大量的溶解氧,但由于高密度的網(wǎng)箱養(yǎng)殖魚類,以及微生物對有機物的降解作用,溶解氧被大量消耗(Holmer,2008),另外水溫的升高也會降低海水溶解氧的飽和濃度;浮游植物大量繁盛,葉綠素含量上升;排泄物和殘餌等有機養(yǎng)料相對豐富,此時水溫更適合微生物繁殖,同時各種細菌病的病原體也會滋生。養(yǎng)殖海區(qū)有顯著性季節(jié)變化的還有8月份鹽度和亞硝酸鹽顯著升高,pH 值顯著降低(P<0.05);非養(yǎng)殖海區(qū)8月份的COD 和富營養(yǎng)鹽指數(shù)也顯著的升高(P<0.05),可見網(wǎng)箱養(yǎng)殖和陸源污染物在夏季對海區(qū)的影響更為廣泛。
2.3.2 N/P 值以及富營養(yǎng)程度的比較
對照本監(jiān)測結(jié)果(見表2),5月份和8月份海區(qū)的N/P 值絕大部分在8 至30 之間。通過t 檢驗,兩個月份的N/P 值差異不顯著(P<0.05)。從表2可以看出,8月份無論網(wǎng)箱內(nèi)還是對照區(qū)的富營養(yǎng)化程度比5月份更加嚴重。由此可以看出,8月份更為活躍的網(wǎng)箱養(yǎng)殖業(yè)和陸源污染給水體帶來明顯的營養(yǎng)鹽污染。
2011年8月三都澳海區(qū)內(nèi)網(wǎng)箱養(yǎng)殖的大黃魚等普遍暴發(fā)了刺激隱核蟲病,細菌病爆發(fā)期間hydrolab DS5X 水質(zhì)探頭每日14 點的數(shù)據(jù)經(jīng)校準后如表3 所示。
海水水質(zhì)和死亡量各因素均值、標準差和變異系數(shù)見表4。
各因素間相關(guān)系數(shù)列于表5。網(wǎng)箱養(yǎng)殖大黃魚在細菌病暴發(fā)時期內(nèi)的死亡量與海水溫度呈極顯著正相關(guān)(P<0.01),與鹽度呈極顯著負相關(guān)(P<0.01),與pH、葉綠素和溶解氧無顯著相關(guān)。
表4 各因素均值、標準差和變異系數(shù)統(tǒng)計表
表5 各水質(zhì)因素與死亡數(shù)間的相關(guān)系數(shù)表
大潮期間(8月27-9月1日,農(nóng)歷28 至初4),海水交換速率快,受外海高鹽海水的影響鹽度升高,網(wǎng)箱內(nèi)和對照區(qū)的溫度和鹽度變化趨于平衡,死亡數(shù)量下降;小潮期間(9月2-9月10日,農(nóng)歷初5 至13),海水交換速率緩慢,網(wǎng)箱內(nèi)(淺水海區(qū)) 的水溫要比對照區(qū)略高,由于沿岸淡水河流的注入,也使沿岸的網(wǎng)箱養(yǎng)殖海區(qū)鹽度下降,正逢夏季海水富營養(yǎng)化程度和細菌密度較高,魚類死亡數(shù)量增加。
2.5.1 底棲生物的種類組成和優(yōu)勢種
三沙灣海區(qū)有底棲生物3 門17 種,其中環(huán)節(jié)動物種類最多,有9 種,占總種數(shù)的52.9 %;甲殼動物6 種,占總種數(shù)的35.3%;棘皮動物2 種,占總種數(shù)的11.8%。
三沙灣海域底棲生物豐度優(yōu)勢種有4 種(優(yōu)勢度(Y) ≥0.05)。
優(yōu)勢種分別為多毛類的半球索沙蠶(Lumbrineris spharocephala),甲殼動物的日本鼓蝦(Alpheus japonicus)、變態(tài)蟳(Charybdis variegata),棘皮動物的馬氏刺蛇尾(Ophiothrix marenzelleri) (表6)。
2.5.2 各站位底棲生物的群落多樣性
分析生物群落的多樣性一般考慮群落中物種的豐富性和異質(zhì)性。三沙灣各站位的底棲生物的種類數(shù)、棲息密度、生物量、物種多樣性指數(shù)、均勻度和豐度見表7 所示。
表6 底棲生物優(yōu)勢種
由于魚類不能攝食和吸收全部的餌料,大約只有1/5 的餌料被有效利用,殘餌和魚糞沉降到海底成為底質(zhì)污染的主要來源(董雙林等,2000),易造成養(yǎng)殖區(qū)沉積物有機碳、硫化物含量和底層海水DO 濃度的超標(黃洪輝等,2007),TP 及各形態(tài)P 含量的高值都出現(xiàn)在青山島和斗帽島之間、三都島南岸及礁頭、城澳北岸等大黃魚養(yǎng)殖海域,低值出現(xiàn)在非養(yǎng)殖區(qū)(吳芳等,2005)。從表7 中可得本次監(jiān)測結(jié)果三都澳海區(qū)非養(yǎng)殖區(qū)的7 號站位的底棲生物種類數(shù)最高,物種多樣性指數(shù)和豐度也明顯高于養(yǎng)殖區(qū)(P<0.05)。
表7 養(yǎng)殖區(qū)和對照區(qū)海域各站位底棲生物的種類數(shù)、生物量、密度、物種多樣性指數(shù)、均勻度指數(shù)和豐富度指數(shù)
2.5.3 ABC 曲線
通過三沙灣海區(qū)3 個站位底棲生物的豐度/生物量比較曲線(圖2),可以看出3 個站位底棲生物的生物量和豐度曲線相距較近且生物量曲線在豐度曲線之上,說明了三沙灣海區(qū)的底棲生物均受到較小程度的擾動。非養(yǎng)殖區(qū)的7 號點豐度和生物量相差的距離要比養(yǎng)殖海區(qū)的4 號點和8 號點更為接近,這說明非養(yǎng)殖區(qū)對底棲生物的擾動更強烈。非養(yǎng)殖海區(qū)的干擾主要來自航道清淤和輪船行駛帶來的底質(zhì)攪動和翻耕,對大型底棲生物有嚴重的影響。養(yǎng)殖區(qū)雖然會降低底層水質(zhì)溶解氧、增加底質(zhì)有機物(黃洪輝等,2007;吳芳等,2005),并最終影響大型底棲動物的豐富度、多樣性等(Pearson et al,1978) ,但是也杜絕了航道清淤等其他的人為擾動,對保護大型底棲生物有一定的積極意義(韓慶喜等,2011)。
圖2 三沙灣海區(qū)底棲生物豐度/生物量比較曲線
三沙灣大部分海區(qū)已呈現(xiàn)富營養(yǎng)化,主要超標污染物是無機氮和磷酸鹽。養(yǎng)殖區(qū)和對照區(qū)的對比顯示,高密度網(wǎng)箱養(yǎng)殖改變了海區(qū)的水質(zhì)環(huán)境,使海域細菌濃度升高,營養(yǎng)鹽特別是活性磷酸鹽濃度升高,N/P 值降低,富營養(yǎng)化程度加劇。由于季節(jié)變化和高密度網(wǎng)箱養(yǎng)殖的影響,8月份水溫升高的同時,pH 和DO 值顯著降低,弧菌密度升高。網(wǎng)箱內(nèi)的魚類體質(zhì)下降,攝食量降低。由于此時正處于高溫期,網(wǎng)箱內(nèi)的水體富營養(yǎng)化程度比較高,浮游生物大量繁殖,細菌密度明顯增多,由刺激隱核蟲病引起的繼發(fā)性細菌性感染愈發(fā)嚴重。由水質(zhì)環(huán)境和網(wǎng)箱養(yǎng)殖魚類死亡率的相關(guān)性分析可以得出,魚類病情的危害程度與三沙灣水文環(huán)境有一定的關(guān)系。大潮期間死亡率下降,小潮期間死亡率上升。
本次三沙灣底棲生物調(diào)查的優(yōu)勢種類為多毛類和甲殼動物,網(wǎng)箱養(yǎng)殖降低了底棲生物的種類和物種多樣性,但是對保護大型底棲生物有一定的積極意義。
針對三都澳海區(qū)網(wǎng)箱養(yǎng)殖魚類高溫期細菌病頻發(fā)的現(xiàn)狀,為改善養(yǎng)殖區(qū)自然環(huán)境,促進網(wǎng)箱養(yǎng)殖業(yè)的可持續(xù)健康發(fā)展,現(xiàn)提出以下幾點措施:1、規(guī)劃網(wǎng)箱布局,降低養(yǎng)殖密度,合理投餌(廖紅梅 等,2011;Cai et al,2007),加強環(huán)保宣傳工作增強漁民的衛(wèi)生意識,減輕網(wǎng)箱內(nèi)的富營養(yǎng)化程度。2、增加養(yǎng)殖容量,開發(fā)推廣健康養(yǎng)殖技術(shù)。漁排內(nèi)部要適當留出空網(wǎng)箱以加快網(wǎng)箱內(nèi)部區(qū)域的水流循環(huán),網(wǎng)箱內(nèi)可適當混養(yǎng)貝類、甲殼類、大型藻類或濾食性魚類,可有效地降低水體的富營養(yǎng)化程度,促進網(wǎng)箱內(nèi)營養(yǎng)物質(zhì)的可循環(huán)再利用(董雙林等,1998;蔡惠文 等,2009)。3、加強海洋環(huán)境監(jiān)測和監(jiān)督管理(Mohammad et al,2008)。統(tǒng)籌海陸污染物排放,控制污染物總量。國家海洋局東海分局目前在現(xiàn)有三沙灣增養(yǎng)殖區(qū)大面監(jiān)測的基礎(chǔ)上逐步增加水質(zhì)在線監(jiān)控探頭和遠紅外視頻監(jiān)控探頭,可以實時反映海水養(yǎng)殖污染狀況,以建立海洋環(huán)境預(yù)警系統(tǒng),并監(jiān)測主要江河和排污口的污染物入海通量。為進一步研究三沙灣網(wǎng)箱養(yǎng)殖容量,科學劃分海洋功能區(qū)域提供參考資料。
Cai H,Sun Y,2007. Management of marine cage aquaculture:environmental carrying capacity method based on dry feed conversion rate.Environment Science Pollution Research,14(7):463-469.
Holmer M,2008. Aquaculture in the Ecosystem. Scotland:Springer:1-46.
Mohammad A,F(xiàn)iras O,Tariq A,et al,2008. Neighboring cage fish farming affecting water and seabed quality of the Jordanian northern gulf of Aqaba,Red Sea.Ocean Science Journal,43(1):9-16.
Pearson T H,Rosenberg R,1978. Macrobenthic succession in relation to organic enrichment and pollution of the marine environment. Oceanographic Marine Biology Annual Review,16:229-311.
Vergnon R,Blanchard F,2006. Evaluation of trawling disturbance on macrobenthic invertebrate com munities in the Bay of Biscay,F(xiàn)rance:Abundance Biomass Comparison (ABC method). Aquatic Living Resources,19:219-228.
GB 17378-2007,2008.海洋監(jiān)測規(guī)范.北京:中國標準出版社.
GB 30972-1997,1998.海水水質(zhì)標準.北京:中國標準出版社
GB 17378.2-2007,2008.海洋監(jiān)測規(guī)范.北京:中國標準出版社.
蔡惠文,任永華,孫英蘭,等,2009.海水養(yǎng)殖環(huán)境容量研究進展.海洋通報,28(2):109-115.
董雙林,李德尚,潘克厚,1998.論海水養(yǎng)殖的養(yǎng)殖容量.中國海洋大學學報(自然科學版),28(2):253-258.
董雙林,潘克厚,Brockmann U,2000.海水養(yǎng)殖對沿岸生態(tài)環(huán)境影響的研究進展.青島海洋大學學報(自然科學版),30(4):575-582.
郭衛(wèi)東,章小明,楊逸萍,等,1998.中國近岸海域潛在性富營養(yǎng)化程度的評價.臺灣海峽,17(1):64-70.
韓慶喜,李寶泉,韓秋影,等,2011.漁業(yè)捕撈對威海港附近海域底上大型底棲群落結(jié)構(gòu)影響的初步研究.海洋通報,30(2):121-126.
黃標,錢魯閩,劉家富,等,2002.福建三都澳水產(chǎn)養(yǎng)殖區(qū)水體氮磷含量及潛在性富營養(yǎng)化程度分析.臺灣海峽,21(4):411-415.
黃洪輝,林欽,甘居利,等,2007.大鵬澳海水魚類網(wǎng)箱養(yǎng)殖對沉積環(huán)境的影響.農(nóng)業(yè)環(huán)境科學學報,26(1):75-80.
廖紅梅,林培華,高健,2011.寧德市大黃魚成魚養(yǎng)殖現(xiàn)狀及對策.貴州農(nóng)業(yè)科學,39(1):165-168
林克冰,周宸,劉家富,等,1999.海水網(wǎng)箱養(yǎng)殖大黃魚病原菌研究.海洋科學,128(4):58-62.
林榮根,1996.海水富營養(yǎng)化水平評價方法淺析.海洋環(huán)境科學,15(2):28-31.
宋蘭芬,2007.福建海水養(yǎng)殖產(chǎn)量與污染關(guān)系的研究.碩士學位論文.福建:福建農(nóng)林大學
蘇躍中,2009.福建省主要養(yǎng)殖區(qū)海水網(wǎng)箱養(yǎng)殖魚類刺激隱核蟲病的調(diào)查及防控對策.水產(chǎn)科技情報,36(1):4-7.
吳芳,郭衛(wèi)東,鄭佩如,等.2005.三都澳養(yǎng)殖海域沉積物中P 的形態(tài)分布特征.海洋環(huán)境科學,24(4):24-27.