膠州灣女姑口海域大型底棲動(dòng)物群落特征與生態(tài)健康評(píng)價(jià)*

劉欣禹，齊衍萍，沙婧婧，徐子鈞，菅瀟揚(yáng)，杜小媛，溫國(guó)義，周瑞佳**

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

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

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

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

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

1 材料與方法

1.1 樣品采集與處理

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

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

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

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

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

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

(1)

(2)

d=(S-1)/log2N，

(3)

Y=Pi×F，

(4)

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

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

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

1.2.3 底棲生境質(zhì)量評(píng)價(jià)

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

表1 AMBI和M-AMBI擾動(dòng)等級(jí)、生態(tài)環(huán)境質(zhì)量狀況標(biāo)準(zhǔn)

2 結(jié)果與分析

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

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

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

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

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

2.1.2 生物密度與生物量分布

調(diào)查海域大型底棲動(dòng)物平均生物密度為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號(hào)站春季調(diào)查發(fā)現(xiàn)大量菲律賓蛤仔，因此生物量生物密度高。在季節(jié)變化方面，春季生物密度、生物量較高，而且種類組成分析發(fā)現(xiàn)，春季軟體動(dòng)物和節(jié)肢動(dòng)物種類數(shù)量明顯增多。

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

2.1.3 群落多樣性特征

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

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

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

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

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

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

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

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

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

3 討論

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

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

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

膠州灣女姑口附近海域大型底棲動(dòng)物秋季Shannon-Wiener多樣性指數(shù)為2.53，春季Shannon-Wiener多樣性指數(shù)為3.01。李新正等[31]2001年調(diào)查發(fā)現(xiàn)，膠州灣底棲動(dòng)物的物種多樣性指數(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年對(duì)膠州灣中部海域調(diào)查顯示，膠州灣大型底棲動(dòng)物春季和秋季Shannon-Wiener多樣性指數(shù)較低，平均值分別為2.59和2.50；全秋梅等[11]2016-2017年調(diào)查膠州灣大型底棲動(dòng)物春季和秋季Shannon-Wiener多樣性指數(shù)分別為2.61和2.34；丁敬坤等[16]2020年調(diào)查膠州灣發(fā)現(xiàn)，女姑口附近站位Shannon-Wiener多樣性指數(shù)為1.67；近20年調(diào)查顯示，膠州灣大型底棲動(dòng)物的Shannon-Wiener多樣性指數(shù)逐漸降低，從側(cè)面說(shuō)明該區(qū)域大型底棲動(dòng)物受到輕度污染。

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

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

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

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

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

4 結(jié)論

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