廣西大學資源環(huán)境與材料學院 曾恒，孫翔

一、引言

牡蠣是世界范圍的養(yǎng)殖品種，也是中國最主要的養(yǎng)殖貝類之一。截至2019年，近海養(yǎng)殖牡蠣產(chǎn)量已經(jīng)達到了5225595噸，占海水貝類養(yǎng)殖產(chǎn)量的38.47%。近江牡蠣的生物沉積作用使營養(yǎng)物質滯留在近江牡蠣養(yǎng)殖區(qū)及周邊海域的底質中[1]，導致近江牡蠣養(yǎng)殖區(qū)域底質的有機物富集效應。近江牡蠣同化作用后產(chǎn)生的糞便及過剩攝食產(chǎn)生的假糞形成生物沉積物，其中一部分經(jīng)過微生物的礦化分解再生為營養(yǎng)鹽[2]進入水體，一部分生物沉積物沉積在養(yǎng)殖區(qū)域及周邊海域的底質中被埋藏。另外一部分有機氮進入沉積物缺氧層，在發(fā)硝化細菌的作用下被還原為N2[2]。欽州灣作為全國最大的近江牡蠣天然苗種繁殖區(qū)[3]，關于近海養(yǎng)殖活動對海域水質變化的研究少有報道[4]。本研究通過2019年夏季和秋季對欽州灣養(yǎng)殖區(qū)的葉綠素a及無機氮進行監(jiān)測，分析欽州灣養(yǎng)殖區(qū)域與非養(yǎng)殖區(qū)域的葉綠素a濃度、無機氮濃度及其變化趨勢，為探究近江牡蠣的生物沉積作用對欽州灣水質的影響提供了科學依據(jù)。

二、數(shù)據(jù)與方法

本研究基于2019年4月欽州灣的水質監(jiān)測數(shù)據(jù)，分析欽州灣水質的空間分布特征。將普通克里金插值法與反距離加權法進行精度比較后，選用精度更高、制圖效果更好的普通克里金插值法對監(jiān)測站點數(shù)據(jù)進行插值分析，并利用皮爾遜相關系數(shù)進行葉綠素a和無機氮相關性分析。

圖1 欽州灣營養(yǎng)鹽流動圖

三、結果與討論

（一）葉綠素a分布特征

欽州灣2019年4月份葉綠素a濃度分布特征見圖2。葉綠素a濃度范圍為5.01-8.64μg/L。區(qū)域分布上，養(yǎng)殖區(qū)域葉綠素a濃度明顯低于非養(yǎng)殖區(qū)，浮筏養(yǎng)殖區(qū)的葉綠素a濃度明顯低于沉排養(yǎng)殖區(qū)，茅尾海葉綠素a濃度高于灣外濃度。近江牡蠣為濾食性生物，它能大量濾食水體中的浮游植物和顆粒物質，從而對浮游植物產(chǎn)生下行控制，葉綠素a指標間接反映了海域的浮游植物生物量，因而養(yǎng)殖區(qū)葉綠素濃度明顯低于非養(yǎng)殖區(qū)。欽州灣灣區(qū)呈葫蘆狀，灣口有欽江、茅嶺江等江河的淡水注入，使得灣區(qū)海水的鹽度呈現(xiàn)從灣頂?shù)綖晨谥饾u遞減的梯度分布。每年5-7月，近江牡蠣養(yǎng)殖的采苗階段將會在灣頂鹽度較低的沙井沉排養(yǎng)殖區(qū)進行。由于本次選取的研究時間為4月，外海育肥期的近江牡蠣被拉回內灣七十二涇進行產(chǎn)卵，茅尾海內部的采苗活動沒有開展，導致浮筏養(yǎng)殖區(qū)的葉綠素a濃度明顯低于沉排養(yǎng)殖區(qū)。欽州灣屬于半封閉海灣，水動力條件較弱，營養(yǎng)鹽不斷累積于茅尾海。同時受茅嶺江、欽江和大風江等河流輸入的影響，大量營養(yǎng)物質輸入促進了浮游植物的大規(guī)模繁殖。同時，在茅尾海內部，只有茅嶺江保苗區(qū)進行近江牡蠣的養(yǎng)殖活動，近江牡蠣的下行控制得到減弱[5]，進而導致尾海葉綠素a濃度高于外灣濃度。

圖2 欽州灣2019年4月葉綠素a濃度分布特征圖

（二）無機氮分布特征

欽州灣2019年4月份無機氮(DIN)濃度分布特征見圖3。無機氮(DIN)濃度范圍為27.66-36.98μmol/L。區(qū)域分布上，茅嶺江和欽江的入海區(qū)域無機氮的濃度高于茅尾海其余區(qū)域。針對灣頸七十二徑區(qū)域，靠近灣外的地區(qū)無機氮的濃度高于靠近灣外的區(qū)域。這表明隨著欽州灣降雨量和徑流量的增加，沿岸流域農耕施肥中大量的化肥和居民生活排污水大量流入欽州灣，且保持較高濃度值，表明廣西北部灣近岸海域是所有營養(yǎng)鹽的匯入點[6]，該時期營養(yǎng)鹽的陸源輸入是欽州灣的主要污染源。同時七十二徑作為近江牡蠣的養(yǎng)殖區(qū)域，氨氮是近江牡蠣的主要排泄物，表明了近江牡蠣的呼吸排泄活動對養(yǎng)殖區(qū)域產(chǎn)生了不可忽視的影響。

圖3 欽州灣2019年4月無機氮濃度分布特征圖

（三）營養(yǎng)鹽結構及其環(huán)境因子的關系

海洋環(huán)境是由化學、物理等要素組成的復合體，各種環(huán)境因素隨著物理條件或者環(huán)境條件的變化而變化。本文通過分析4月欽州灣葉綠素a和無機氮的相關性系數(shù)，顯示相關性為0.3868。結果表明，欽州灣葉綠素a的濃度和無機氮濃度的相關性不顯著。

（四）結論

本文探討了2019年欽州灣海區(qū)葉綠素a、無機氮的時空分布特征，主要得出以下結論:（1）區(qū)域分布上，養(yǎng)殖區(qū)域葉綠素a濃度明顯低于非養(yǎng)殖區(qū)，浮筏養(yǎng)殖區(qū)的葉綠素a濃度明顯低于沉排養(yǎng)殖區(qū)，茅尾海葉綠素a濃度高于灣外濃度。（2）茅嶺江和欽江的入海區(qū)域無機氮的濃度高于茅尾海區(qū)域。針對灣頸七十二徑區(qū)域區(qū)域，靠近灣外的地區(qū)無機氮的濃度高于靠近灣外的區(qū)域。（3）通過皮爾遜相關系數(shù)，葉綠素a和無機氮不存在相關性。

盡管本研究對欽州灣2019年4月葉綠素a和無機氮濃度進行了分析，但是影響海洋環(huán)境的影響因子很多，在今后條件允許的情況下，本研究會繼續(xù)開展欽州灣理化參數(shù)(DO、pH和COD等參數(shù))和營養(yǎng)鹽指標(包括NH4+-N、NO2—-N、NO3—-N-等參數(shù))的監(jiān)測分析，為探究欽州灣近江牡蠣的生物沉積作用對欽州灣水質的影響提供科學依據(jù)。