劉欣禹，齊衍萍，沙婧婧，徐子鈞，菅瀟揚，杜小媛，溫國義，周瑞佳**

(1.國家海洋局北海環(huán)境監(jiān)測中心，山東青島 266000；2.自然資源部渤海生態(tài)預警與保護修復重點實驗室，山東青島 266000；3.山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室，山東青島 266000)

大型底棲動物作為海洋生態(tài)結(jié)構(gòu)的重要組成部分，通過攝食、掘穴等活動參與碳、氮、磷、硫等元素的循環(huán)，有助于加快底層生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量傳遞[1,2]。由于活動能力弱、棲息地相對穩(wěn)定和對海洋環(huán)境反應(yīng)敏感等特點，所以大型底棲動物群落的物種組成通常會隨著非生物條件和人類干擾的變化而變化，這使得大型底棲動物經(jīng)常被用作評價底棲生境健康狀況的生物指標[2,3]，對海洋生態(tài)環(huán)境研究具有積極意義。

將底棲生物作為海洋生態(tài)環(huán)境監(jiān)測和生態(tài)系統(tǒng)健康評價的生物指標已得到國內(nèi)外科學家的廣泛認可[4]。多樣性指數(shù)、AZTI海洋生物指數(shù)(AMBI)、多變量AMBI (M-AMBI) 均能在一定程度上反映底棲動物受擾動狀態(tài)，可用于評價底棲環(huán)境狀態(tài)[5]。其中，AMBI和M-AMBI是Borja等[6]在海洋生物指數(shù)的基礎(chǔ)上創(chuàng)建的。在我國及歐洲沿海，這3種指數(shù)已被用于港灣、近岸海域以及河口區(qū)的底棲生態(tài)環(huán)境質(zhì)量評價[7-10]。

膠州灣位于山東半島南部，屬于相對獨立的半封閉海灣，女姑口海域位于青島市城陽區(qū)，在膠州灣的東北部，具有典型的泥質(zhì)底質(zhì)。20世紀80年代開始，膠州灣東部淺海海域的養(yǎng)殖業(yè)迅速崛起，90年代周邊城區(qū)工業(yè)和旅游業(yè)等繁榮發(fā)展，近岸海域資源被快速開發(fā)利用[11]。膠州灣環(huán)海高速公路從女姑口上通過，附近灘涂上主要開展底棲型貝類如菲律賓蛤仔(Ruditapesphilippinarum)、縊蟶(Sinonovaculaconstricta)的養(yǎng)殖[12,13]。麻驁等[14]研究膠州灣養(yǎng)殖水域和自然水域的底棲狀況發(fā)現(xiàn)，菲律賓蛤仔養(yǎng)殖活動壓力已造成養(yǎng)殖區(qū)域及其鄰近自然海域的大型底棲群落結(jié)構(gòu)發(fā)生改變。崔正國等[15]對膠州灣底播菲律賓蛤仔的質(zhì)量安全風險進行了評價。其他研究主要基于膠州灣大型底棲動物進行海洋健康狀況評估[16]，針對膠州灣女姑口附近海域調(diào)查相對較少。

本研究基于2018年11月和2019年4月兩個航次膠州灣女姑口鄰近海域的大型底棲動物調(diào)查和環(huán)境要素調(diào)查數(shù)據(jù)，對膠州灣女姑口鄰近海域生物種類組成和生物群落結(jié)構(gòu)進行研究，探討群落結(jié)構(gòu)與環(huán)境因子的相互關(guān)系，并對該海域進行生態(tài)系統(tǒng)健康評價，以了解膠州灣女姑口近海底棲生物資源現(xiàn)狀及生態(tài)健康狀況，為膠州灣整體生態(tài)保護與管理提供數(shù)據(jù)支持。

1 材料與方法

1.1 樣品采集與處理

2018年11月、2019年4月在膠州灣女姑口近海設(shè)置12個站位進行采樣監(jiān)測，站位布設(shè)見圖1。

大型底棲動物用0.1 m2抓斗式采泥器采集，每站采樣4次，樣品用網(wǎng)目為0.5 mm的過篩器分選，用5%福爾馬林溶液固定后，帶回實驗室進行定量定性分析，物種鑒定及名錄登記均通過WoRMS網(wǎng)站進行校對和復核。站位布設(shè)、生物密度及生物量分布圖由ArcGIS 10.6軟件繪制。

測定的環(huán)境參數(shù)包括海水水深(D)、水溫(WT)、鹽度(S)、pH值、溶解氧(DO)、化學需氧量(COD)、磷酸鹽(DIP)、無機氮(DIN)、懸浮物(SS)、表層沉積物底質(zhì)類型及中值粒徑(MD)。

1.2 數(shù)據(jù)處理

1.2.1 群落結(jié)構(gòu)及優(yōu)勢種計算

使用Shannon-Wiener多樣性指數(shù)(H′)[17]、Pielous均勻度指數(shù)(J)[17]和Margalef物種豐富度指數(shù)(d)[17]分析生物群落結(jié)構(gòu)。使用物種優(yōu)勢度指數(shù)(Dominance index,Y)判定優(yōu)勢種[18]。計算公式如下:

(1)

(2)

d=(S-1)/log2N，

(3)

Y=Pi×F，

(4)

式中，Pi=ni/N(ni是第i個物種的個體數(shù)，N是全部物種的個體數(shù))；S為種類數(shù)；F為該種的出現(xiàn)頻率，當Y>0.02時[19]，即為該海域底棲動物優(yōu)勢種。H′值等于0，說明無底棲生物，即嚴重污染；0-1為重污染；1-2為中度污染；2-3為輕度污染；大于3為清潔[17]。

1.2.2 環(huán)境因子相關(guān)分析

使用SPSS 18.0軟件，進行單因素方差分析(One-way ANOVA)，分析各季節(jié)環(huán)境要素的差異性；采用Canoco 5.0軟件對大型底棲動物優(yōu)勢種生物密度和環(huán)境因子進行典范對應(yīng)分析(CCA)，物種矩陣經(jīng)lg(x+1)轉(zhuǎn)換[20]，探究環(huán)境因子對大型底棲動物優(yōu)勢種生物密度變化的影響。

1.2.3 底棲生境質(zhì)量評價

AMBI和M-AMBI采用AMBI V5.0軟件計算，先除去非大型底棲無脊椎動物的種類，其余大型底棲動物的生態(tài)分組根據(jù)文獻[10,21]規(guī)定執(zhí)行。AMBI值為0.0-7.0，數(shù)值越低表示生態(tài)健康狀況越好；M-AMBI值為0.00-1.00，數(shù)值越高表示生態(tài)健康狀況越好。AMBI和M-AMBI的評價標準[10,22]詳見表1。

表1 AMBI和M-AMBI擾動等級、生態(tài)環(huán)境質(zhì)量狀況標準

2 結(jié)果與分析

2.1 群落結(jié)構(gòu)

2.1.1 種類組成與優(yōu)勢種

調(diào)查共發(fā)現(xiàn)82種底棲生物，包括扁形動物、刺胞動物、紐形動物、環(huán)節(jié)動物、軟體動物、節(jié)肢動物、棘皮動物、半索動物和脊索動物9個門類，其中環(huán)節(jié)動物多毛類29種、節(jié)肢動物甲殼類29種、軟體動物18種，以上3個門類種數(shù)占總種數(shù)的92.68%，是女姑口鄰近海域主要的大型底棲動物類群，其他各門類均1種。春季(4月)采集到的大型底棲動物種類較多，為67種，其中環(huán)節(jié)動物27種，軟體動物15種，節(jié)肢動物21種，扁形動物、刺胞動物、紐形動物、棘皮動物各1種。秋季(11月)采集到大型底棲動物46種，其中環(huán)節(jié)動物21種，軟體動物7種，節(jié)肢動物15種，紐形動物、尾索動物、脊索動物各1種。

優(yōu)勢種主要為環(huán)節(jié)動物多毛類和軟體動物，兩個季節(jié)均表現(xiàn)出明顯優(yōu)勢的種為寡鰓齒吻沙蠶(Nephtysoligobranchia)、寡節(jié)甘吻沙蠶(Glycindegurjanovae)和菲律賓蛤仔。其他優(yōu)勢種均為環(huán)節(jié)動物多毛類及小型貝類，膠州灣女姑口鄰近海域大型底棲動物有小型化趨勢。節(jié)肢動物甲殼類在本調(diào)查中共發(fā)現(xiàn)29種，其中秋季15種，其密度占秋季總密度的12.03%，密度占比最高的種為長尾亮鉤蝦(Photislongicaudata)和極地蚤鉤蝦(Pontocratesaltamarinus)，兩種密度占比均為2.95%，站位出現(xiàn)率均為41.67%；春季21種，其密度占春季總密度的13.78%，密度占比最高的種類為東方長眼蝦(Ogyridesorientalis)(密度占比2.23%，站位出現(xiàn)率33.33%)。雖然甲殼類是種類組成的主要類群，但是密度占比未超過14%，站位出現(xiàn)率未超過45%，優(yōu)勢度較低，并沒有出現(xiàn)甲殼類優(yōu)勢種。調(diào)查海域大型底棲動物優(yōu)勢種及優(yōu)勢度指數(shù)見表2。

表2 秋季和春季調(diào)查海域大型底棲動物優(yōu)勢種及優(yōu)勢度指數(shù)

2.1.2 生物密度與生物量分布

調(diào)查海域大型底棲動物平均生物密度為665.83 ind./m2，平均生物量為358.65 g/m2。調(diào)查海域2018年秋季平均生物密度為485.00 ind./m2，平均生物量為235.80 g/m2；2019年春季平均生物密度為846.67 ind./m2，平均生物量為481.60 g/m2。由圖2可知，秋季生物密度高值區(qū)主要位于紅島灘涂附近海域，生物量高值區(qū)位于紅島外側(cè)海域；春季生物密度及生物量的高值區(qū)主要集中在墨水河口附近海域，墨水河口2號站春季調(diào)查發(fā)現(xiàn)大量菲律賓蛤仔，因此生物量生物密度高。在季節(jié)變化方面，春季生物密度、生物量較高，而且種類組成分析發(fā)現(xiàn)，春季軟體動物和節(jié)肢動物種類數(shù)量明顯增多。

圖2 秋季和春季調(diào)查海域大型底棲動物生物密度、生物量分布

2.1.3 群落多樣性特征

由表3可知，調(diào)查海域大型底棲動物春季Shannon-Wiener多樣性指數(shù)(H′)平均值高于秋季，兩季節(jié)的平均多樣性指數(shù)為2.77，說明調(diào)查海域大型底棲動物受到輕度污染；Pielous均勻度指數(shù)(J)略低于秋季，Margalef豐富度指數(shù)(d)高于秋季。

表3 秋季和春季調(diào)查海域大型底棲動物群落多樣性指數(shù)

2.2 大型底棲動物與環(huán)境因子的相關(guān)性分析

調(diào)查海域水深為2.0-7.5 m，平均水深為4.17 m，呈現(xiàn)梯度遞增趨勢；12個調(diào)查站位的底質(zhì)為黏土質(zhì)粉砂；水溫(WT)、鹽度(S)、溶解氧(DO)、化學需氧量(COD)、磷酸鹽(DIP)、懸浮物(SS)季節(jié)間呈極顯著差異(P<0.01)；無機氮(DIN)季節(jié)間呈顯著差異(P<0.05)；季節(jié)間pH值與中值粒徑(MD)無顯著差異。

采用《海水水質(zhì)標準》[23]對調(diào)查海域水質(zhì)進行評價。調(diào)查海域秋季和春季海水中pH值、溶解氧均滿足一類海水水質(zhì)標準要求，化學需氧量、磷酸鹽、無機氮均超一類海水水質(zhì)標準要求。其中，秋季和春季海水化學需氧量全部符合二類水質(zhì)標準；秋季調(diào)查發(fā)現(xiàn)磷酸鹽和無機氮都超四類水質(zhì)標準，超標率均為75.0%；春季調(diào)查發(fā)現(xiàn)磷酸鹽和無機氮都超四類水質(zhì)標準，超標率分別為20.0%和50.0%。

表4 秋季與春季調(diào)查海域各環(huán)境因子平均值

由圖3可知，兩個排序軸的特征值分別為0.58和0.17，種類與環(huán)境因子排序軸的相關(guān)系數(shù)為0.952和0.754，排序可以較好地反映底棲動物與環(huán)境因子間的關(guān)系。鹽度(S)與第一排序軸呈最大正相關(guān)(0.717)，水深(CD)與第一排序軸呈第二正相關(guān)(0.534)，水溫(WT)與第一排序軸呈最大負相關(guān)(-0.754)。說明影響膠州灣女姑口大型底棲動物分布的主要驅(qū)動因子有鹽度、水溫。

N.ol：Nephtys oligobranchia；G.gu：Glycinde gurjanovae；R.ph：Ruditapes philippinarum；H.fi：Heteromastus filiformis；S.sc：Sternaspis sculata；T.mu：Tharyx multifilis；L.ne：Lagis neapolitana；T.la：Theora lata

2.3 AMBI及M-AMBI結(jié)果

2018年秋季 12個調(diào)查站位AMBI為1.41-3.16，均為輕度擾動；M-AMBI為0.46-0.92，站位13生態(tài)環(huán)境質(zhì)量處于中等級，其他站位為優(yōu)或者良(表5)。2019年春季12個調(diào)查站位AMBI為1.58-3.70，10個調(diào)查站位為輕度擾動，2個站位(站位14、站位3)為中度擾動(表5)。站位14第二機會種(EG Ⅳ)占比為77.8%，站位3第二機會種(EG Ⅳ)占比為62.3%。M-AMBI為0.50-0.91，站位14底棲生態(tài)環(huán)境質(zhì)量處于中等級，其他站位屬于良或者優(yōu)(表5)。

表5 秋季與春季調(diào)查海域AMBI及M-AMBI

3 討論

3.1 大型底棲動物與環(huán)境要素的關(guān)系

本研究結(jié)果表明，對膠州灣女姑口附近海域大型底棲動物分布特征造成影響的主要環(huán)境因子為鹽度和水溫。這與陳曉娟等[24]和畢洪生[25]對膠州灣大型底棲動物的研究結(jié)論相似，即水溫與鹽度影響膠州灣大型底棲動物的生物密度分布。崔雯瑤等[26]對膠州灣2018年春季數(shù)據(jù)進行分析，其調(diào)查時間與本研究調(diào)查時間相近，同樣發(fā)現(xiàn)水溫與鹽度會影響膠州灣大型底棲動物生物密度分布，而其他環(huán)境因子無顯著相關(guān)關(guān)系。膠州灣屬于暖溫帶水域，水深較淺，水溫變化劇烈，水溫影響底棲動物的生長和發(fā)育[27]，從而影響生物的種類數(shù)量、生物量和分布范圍[28]。以往調(diào)查中發(fā)現(xiàn)，鹽度對底棲生物的影響主要集中在河口附近區(qū)域，該處因河流淡水的注入，尤其是在雨季，大量注入的淡水引起鹽度發(fā)生較大的變化，從而影響該處的群落組成和生物密度[29,30]。

3.2 膠州灣女姑口附近海域底棲生態(tài)環(huán)境健康狀況

膠州灣女姑口附近海域大型底棲動物秋季Shannon-Wiener多樣性指數(shù)為2.53，春季Shannon-Wiener多樣性指數(shù)為3.01。李新正等[31]2001年調(diào)查發(fā)現(xiàn)，膠州灣底棲動物的物種多樣性指數(shù)春季、秋季平均值分別為4.59和4.76；王金寶等[32]2005-2009年膠州灣調(diào)查，春季Shannon-Wiener多樣性指數(shù)為0.46-3.66，秋季Shannon-Wiener多樣性指數(shù)0.86-3.61；隋吉星等[33]2010年對膠州灣中部海域調(diào)查顯示，膠州灣大型底棲動物春季和秋季Shannon-Wiener多樣性指數(shù)較低，平均值分別為2.59和2.50；全秋梅等[11]2016-2017年調(diào)查膠州灣大型底棲動物春季和秋季Shannon-Wiener多樣性指數(shù)分別為2.61和2.34；丁敬坤等[16]2020年調(diào)查膠州灣發(fā)現(xiàn)，女姑口附近站位Shannon-Wiener多樣性指數(shù)為1.67；近20年調(diào)查顯示，膠州灣大型底棲動物的Shannon-Wiener多樣性指數(shù)逐漸降低，從側(cè)面說明該區(qū)域大型底棲動物受到輕度污染。

本次調(diào)查結(jié)果顯示，膠州灣女姑口附近海域12個站位AMBI秋季為1.41-3.16，春季為1.58-3.70，說明調(diào)查海域皆受到不同程度的干擾，所有站次均處于輕度擾動以上；M-AMBI秋季為0.46-0.92，春季為0.50-0.91，說明調(diào)查海域部分站位的生態(tài)環(huán)境健康狀況處于中等狀態(tài)；Shannon-Wiener多樣性指數(shù)(H′)、AMBI及M-AMBI評價結(jié)果基本一致。本調(diào)查結(jié)果與丁敬坤等[16]評價結(jié)果相似，即膠州灣海域底棲生態(tài)系統(tǒng)受到中等程度的干擾。

菲律賓蛤仔是膠州灣養(yǎng)殖的主要種類，在菲律賓蛤仔較多的站位，群落Shannon-Wiener多樣性指數(shù)和Pielous均勻度指數(shù)降低，小個體的多毛類，如那不勒斯膜帽蟲(L.neapolitana)、多絲獨毛蟲(T.multifilis)等，生物密度升高，此類生物多為機會種，因此導致AMBI變大、M-AMBI變小，王振鐘[34]調(diào)查也得出相似結(jié)論。養(yǎng)殖過程中的營養(yǎng)鹽排放、養(yǎng)殖生物的逃逸、疾病的傳播等給海洋生態(tài)系統(tǒng)帶來較大影響[35]。菲律賓蛤仔自身覓食及挖掘活動對養(yǎng)殖海域的環(huán)境造成一定的擾動，同時其捕獲過程對表層的沉積環(huán)境造成劇烈的物理擾動，對大型底棲動物的生活環(huán)境造成一定的破壞[16]。本次調(diào)查受中度擾動的站位位于菲律賓蛤仔養(yǎng)殖用海附近，該海域的大型底棲動物的健康狀況可能受到養(yǎng)殖活動干擾。

當群落受到中等程度的污染擾動時，個體較大的種的優(yōu)勢度被削弱，生物量占優(yōu)勢的大個體消失，因此在數(shù)量上占優(yōu)勢的種是個體較小的種[35]，膠州灣歷史調(diào)查中已發(fā)現(xiàn)存在大型底棲動物小型化現(xiàn)象[36-38]，畢洪生等[36]曾在1991-1995年對膠州灣的底棲動物進行連續(xù)監(jiān)測，發(fā)現(xiàn)20世紀90年代與80年代相比，小型底泥食性種類豐度上升，主要種類為絲異須蟲(H.filiformis)和方格獨鰓蟲(Tharyxtesselata)；王洪法等[37]2009年在膠州灣調(diào)查發(fā)現(xiàn)優(yōu)勢種為多毛類動物的不倒翁蟲(S.scutata)和軟體動物菲律賓蛤仔；全秋梅等[11]在膠州灣調(diào)查中發(fā)現(xiàn)優(yōu)勢種以多毛類為主，如絲異須蟲、索沙蠶科(Lumbrineriiae)等。在渤海大型底棲動物調(diào)查中也發(fā)現(xiàn)小型化現(xiàn)象，萊州灣自20世紀90年代以來大型底棲動物群落結(jié)構(gòu)已發(fā)生變化，總體呈現(xiàn)小型化趨勢，主要是小型多毛類動物和甲殼類動物取代大個體的棘皮動物和軟體動物[39]。本次膠州灣女姑口調(diào)查中，環(huán)節(jié)動物多毛類種類共29種，在秋季總密度中占比49.65%，在春季占比50.29%，優(yōu)勢種為菲律賓蛤仔、寡鰓齒吻沙蠶和寡節(jié)甘吻沙蠶，其他優(yōu)勢種均為環(huán)節(jié)動物多毛類及小型貝類，與歷史調(diào)查結(jié)果相近，膠州灣女姑口附近種類組成與優(yōu)勢種呈小型化現(xiàn)象。

造成大型底棲動物小型化的原因主要有水體富營養(yǎng)化、過度捕撈及大型底拖網(wǎng)使用等[40,41]。由于近岸水域通常受到較嚴重的環(huán)境脅迫，近海污染、徑流攜帶大量的有機質(zhì)，以及水流和波浪對海底沉積物和底棲生物群落的影響較大，導致近岸的群落主要由一些小型多毛類耐受種組成，這些種多為機會種[27]。本次調(diào)查區(qū)域位于膠州灣東北部，常年平均水深4 m，水動力交換較弱，墨水河、白沙河及李村河入?？谠诒菊{(diào)查區(qū)域附近，調(diào)查發(fā)現(xiàn)重金屬高值區(qū)主要位于東部李村河口[42,43]。青島市生態(tài)環(huán)境狀況公報顯示，2018年青島市膠州灣優(yōu)良水質(zhì)面積比例為73.7%，李村河、墨水河和大沽河入?？诟浇Ｓ蛩|(zhì)較差[44]。2019年膠州灣東北部海域水質(zhì)較差，主要污染物為無機氮[45]。本次調(diào)查秋季磷酸鹽和無機氮都超四類水質(zhì)標準，超標率均為75.0%；春季磷酸鹽和無機氮都超四類水質(zhì)標準，超標率分別為20.0%和50.0%。磷酸鹽和無機氮超標易導致水體富營養(yǎng)化，水體富營養(yǎng)化對大型底棲動物多樣性產(chǎn)生顯著的負面影響[11]，富營養(yǎng)化加速大型底棲動物的種類更新頻率，導致群落組成變化以及群落功能多樣性降低[46]。本次調(diào)查優(yōu)勢種多絲獨毛蟲、絲異須蟲、那不勒斯膜帽蟲在AMBI生態(tài)組中為第二級機會種(EG Ⅳ)，機會種的高生物密度和高出現(xiàn)率導致部分站位生態(tài)環(huán)境質(zhì)量評價較差。大型底棲動物群落結(jié)構(gòu)發(fā)生變動，出現(xiàn)小型化，在一定程度上反映了近年調(diào)查海域受到污染和擾動[39]。

綜上所述，養(yǎng)殖和污染可能導致膠州灣女姑口附近海域大型底棲動物小型化，該海域的底棲生態(tài)環(huán)境健康狀況受到一定擾動。對此，建議加強對膠州灣東北部近岸海域環(huán)境管理，探索有效的綜合治理措施，使膠州灣海洋生態(tài)系統(tǒng)得以持續(xù)健康地發(fā)展。

4 結(jié)論

本次膠州灣女姑口鄰近海域調(diào)查共獲得大型底棲動物82種，以環(huán)節(jié)動物、軟體動物、甲殼動物為主，棘皮動物等其他動物占比很小。兩個調(diào)查季節(jié)中均表現(xiàn)出明顯優(yōu)勢的種為寡鰓齒吻沙蠶、寡節(jié)甘吻沙蠶和菲律賓蛤仔，但是其他優(yōu)勢種均為環(huán)節(jié)動物多毛類及小型貝類，說明膠州灣女姑口鄰近海域大型底棲動物有小型化趨勢，高經(jīng)濟價值種在生物密度和生物量上也有一定程度的下降。典范對應(yīng)分析的結(jié)果顯示，鹽度、水溫是影響大型底棲動物分布特征的主要因子。綜合大型底棲動物群落的Shannon-Wiener多樣性指數(shù)、AMBI及M-AMBI分析結(jié)果，膠州灣女姑口鄰近海域底棲生態(tài)系統(tǒng)受到輕度干擾，紅島灘涂附近海洋生態(tài)環(huán)境質(zhì)量狀況處于一般水平，其他地區(qū)處于優(yōu)良或高等水平，養(yǎng)殖和污染可能導致膠州灣女姑口附近海域大型底棲動物群落結(jié)構(gòu)呈小型化趨勢，并導致該海域的底棲生態(tài)環(huán)境健康狀況受到擾動。膠州灣東北部女姑口海域是人類活動及工業(yè)生活密集區(qū)域，加強該海域的生態(tài)環(huán)境保護具有重要意義。