任曉明，劉陽(yáng)，徐賓鐸，張崇良，任一平,2，程遠(yuǎn)，薛瑩*

( 1. 中國(guó)海洋大學(xué) 水產(chǎn)學(xué)院，山東 青島 266003；2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室 海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實(shí)驗(yàn)室，山東 青島 266237；3. 青島職業(yè)技術(shù)學(xué)院 信息學(xué)院，山東 青島 266555；4. 近海（大連）生態(tài)發(fā)展有限公司，遼寧 大連 116023)

1 引言

海州灣是位于黃海南部沿岸的開放型海灣，地處南北氣候過渡帶，具有明顯的季風(fēng)氣候特點(diǎn)，年降水量近1 000 mm。海州灣及鄰近海域水質(zhì)肥沃，是多種魚類和經(jīng)濟(jì)無脊椎動(dòng)物的產(chǎn)卵場(chǎng)和育幼場(chǎng)。然而，由于近年來捕撈強(qiáng)度增大，海州灣及鄰近海域生態(tài)環(huán)境惡化，生物多樣性和生態(tài)環(huán)境質(zhì)量面臨前所未有的威脅，漁業(yè)資源衰退日益嚴(yán)重，許多傳統(tǒng)經(jīng)濟(jì)魚類目前已難形成漁汛[1]。研究表明，傳統(tǒng)的單魚種的漁業(yè)管理策略已經(jīng)不能適應(yīng)現(xiàn)代漁業(yè)建設(shè)的需要，要采取基于生態(tài)系統(tǒng)的漁業(yè)管理方式，遏制漁業(yè)資源衰退和生態(tài)惡化的趨勢(shì)[2]。為了更加有效地開展?jié)O業(yè)資源開發(fā)及管理工作，促進(jìn)該海域生態(tài)系統(tǒng)功能的恢復(fù)和漁業(yè)資源的可持續(xù)發(fā)展，對(duì)海州灣及鄰近海域生態(tài)系統(tǒng)結(jié)構(gòu)和功能的研究就顯得十分必要。

Ecopath模型[3-4]根據(jù)營(yíng)養(yǎng)動(dòng)力學(xué)原理，構(gòu)建以物質(zhì)平衡為基礎(chǔ)的生態(tài)系統(tǒng)模型，描述平衡生態(tài)系統(tǒng)中各功能群的生物生產(chǎn)和能量流動(dòng)，量化各種生態(tài)系統(tǒng)指標(biāo)，反映生態(tài)系統(tǒng)結(jié)構(gòu)和功能的特征及其變化趨勢(shì)。最初由Polovina[5]使用這種方法來評(píng)估水域生態(tài)系統(tǒng)在穩(wěn)定狀態(tài)下的系統(tǒng)組成生物量和食物消耗，之后結(jié)合Ulanowicz[6]的能量分析生態(tài)學(xué)理論，逐漸發(fā)展成為一種生態(tài)系統(tǒng)能量流動(dòng)分析方法[7]。目前，Ecopath模型已在國(guó)內(nèi)外諸多水域得到廣泛的應(yīng)用和研究[8-10]，為實(shí)現(xiàn)基于生態(tài)系統(tǒng)的漁業(yè)管理提供了必要的數(shù)據(jù)支撐和理論依據(jù)。

本研究基于2018年春季和秋季在海州灣及鄰近海域開展的漁業(yè)資源底拖網(wǎng)調(diào)查，通過構(gòu)建海州灣及鄰近水域生態(tài)系統(tǒng)的Ecopath模型，分析該生態(tài)系統(tǒng)的營(yíng)養(yǎng)相互關(guān)系和特征，旨在為海州灣及鄰近海域漁業(yè)資源的可持續(xù)利用和科學(xué)管理提供理論依據(jù)。

2 材料方法

2.1 研究海域

樣品采自2018年春季（4月）和秋季（9月），在海州灣及鄰近海域（34°20′～35°40′N、119°20′～121°10′ E）進(jìn)行漁業(yè)資源底拖網(wǎng)調(diào)查（圖1）。采用分層隨機(jī)取樣的方法設(shè)計(jì)調(diào)查站位，結(jié)合水深和底質(zhì)等因素的差異將調(diào)查海域劃分為A、B、C、D、E 共5個(gè)區(qū)域，將每個(gè)區(qū)域按照經(jīng)緯度10′×10′劃分成小區(qū)，各區(qū)域隨機(jī)選取的調(diào)查站位數(shù)量分別為2、4、2、7、3，共18個(gè)調(diào)查站位。調(diào)查用船為220 kW單拖漁船，拖速為2.0 kn左右，每站平均拖網(wǎng)時(shí)間約為1 h。調(diào)查網(wǎng)口寬度約為12 m，囊網(wǎng)網(wǎng)目為17 mm。

2.2 Ecopath模型建模

2.2.1 模型原理

Ecopath模型廣泛應(yīng)用于水生生態(tài)系統(tǒng)研究，為基于生態(tài)系統(tǒng)的漁業(yè)管理提供基礎(chǔ)信息。該模型定義生態(tài)系統(tǒng)由一系列生態(tài)關(guān)聯(lián)的功能群構(gòu)成，這些功能群包含1個(gè)或多個(gè)具有相似的生態(tài)功能或種群動(dòng)態(tài)的物種，所有功能群成分基本覆蓋生態(tài)系統(tǒng)能流的全過程。該模型包含1組線性方程，描述了1段特定時(shí)間內(nèi)的系統(tǒng)能量輸入和輸出平衡，公式如下：

式中，Bi和 Bj分別表示功能群 i和 j的生物量，（P/B）i代表功能群i的生產(chǎn)量與生物量的比值，（Q/B）j表示功能群j的消耗量與生物量的比值，EEi（Ecotrophic Efficiency）是指生態(tài)營(yíng)養(yǎng)效率，DCij為被捕食組i占捕食組j的總捕食量的比例，Yi為漁獲量，BAi為生物量累積，Ei為凈遷出量。

對(duì)于每個(gè)功能群，食物組成矩陣DCij及B、P/B、Q/B和EE 4個(gè)基本參數(shù)中的任意3個(gè)都必須輸入，以構(gòu)建模型，其他參數(shù)可以由模型估計(jì)得出。

2.2.2 功能群的劃分及參數(shù)確定

根據(jù)生物種類攝食習(xí)性、生態(tài)學(xué)特征等方面的相似性劃分功能群，同時(shí)將一些具有重要經(jīng)濟(jì)價(jià)值或生態(tài)功能的單一物種設(shè)定為功能群，將海州灣及鄰近海域生態(tài)系統(tǒng)定義為26個(gè)功能群，基本覆蓋該生態(tài)系統(tǒng)能量流動(dòng)的全過程（表1）。

B主要依據(jù)本次漁業(yè)資源底拖網(wǎng)調(diào)查數(shù)據(jù)和參考文獻(xiàn)[11-13]；P/B和Q/B參考發(fā)表的文獻(xiàn)或相似生態(tài)系統(tǒng)的Ecopath模型[10,12-17]；功能群的食物組成主要來源于海州灣及鄰近海域采樣魚類的胃含物分析和相關(guān)文獻(xiàn)數(shù)據(jù)[18-28]。

2.2.3 模型調(diào)試

根據(jù)營(yíng)養(yǎng)動(dòng)力學(xué)原理，Ecopath模型的調(diào)試是為了平衡生態(tài)系統(tǒng)的輸入和輸出，模型平衡需要滿足EE值介于0～1之間。模型參數(shù)化估計(jì)首次運(yùn)行后，會(huì)有一些功能群的EE大于1（不平衡功能群），需要反復(fù)調(diào)整不平衡功能群的各項(xiàng)參數(shù)，以達(dá)到模型平衡狀態(tài)。另外，還需要考慮總效率GE值，表示為生產(chǎn)量和消耗量的比值（P/Q）。在許多情況下，GE值應(yīng)介于 0.1～0.3 之間[3]。

圖1 海州灣及鄰近海域調(diào)查區(qū)域Fig.1 The bottom trawl survey areas in Haizhou Bay and adjacent waters

表1 海州灣及鄰近海域生態(tài)系統(tǒng)Ecopath模型的功能群及主要種類組成Table1 Functional groups and main species based on Ecopath model in the Haizhou Bay and adjacent waters

通過Pedigree指數(shù)[4]分析數(shù)據(jù)來源的可靠性和模型質(zhì)量，量化模型輸入?yún)?shù)的不確定性，對(duì)模型中的輸入?yún)?shù)的來源及質(zhì)量進(jìn)行分析。Pedigree指數(shù)范圍為0～1.0，1.0代表數(shù)據(jù)質(zhì)量較高，通過精確采樣獲得；0代表數(shù)據(jù)來源模糊，數(shù)據(jù)參考其他模型或文獻(xiàn)等。

2.3 混合營(yíng)養(yǎng)效應(yīng)分析

通過混合營(yíng)養(yǎng)效應(yīng)分析 (Mixed Trophic Impact Analysis, MTI)[29]研究目標(biāo)海域生態(tài)系統(tǒng)各功能群之間的營(yíng)養(yǎng)相互關(guān)系。MTI分析通過構(gòu)建n×n的矩陣M來計(jì)算，其元素mij表示建模食物網(wǎng)中的每對(duì)功能群，該分析描述了功能群i生物量的極小變動(dòng)對(duì)生態(tài)系統(tǒng)其他功能群j的營(yíng)養(yǎng)影響—包括直接的（捕食）和間接的（競(jìng)爭(zhēng)）影響[4,30-31]。mij值的正或負(fù)代表由于影響組i生物量的極小變動(dòng)造成被影響組j生物量的增加或減少。矩陣M的負(fù)元素表示負(fù)影響，即捕食者對(duì)被捕食者功能群產(chǎn)生的直接或間接的影響，同理，正元素表示被捕食者對(duì)捕食者功能群產(chǎn)生正影響。因此，矩陣M的正負(fù)元素與上行效應(yīng)和下行效應(yīng)相關(guān)聯(lián)[31]。

3 結(jié)果

3.1 海州灣及鄰近海域生態(tài)系統(tǒng)模型

從表2可以看出，海州灣及鄰近海域生態(tài)系統(tǒng)的營(yíng)養(yǎng)級(jí)范圍從1.00到4.19。營(yíng)養(yǎng)級(jí)2.00～2.99包括所有的底棲無脊椎動(dòng)物和浮游動(dòng)物；魚類，頭足類和口蝦蛄的營(yíng)養(yǎng)級(jí)均大于3.00，方氏云鳚在魚類中營(yíng)養(yǎng)級(jí)最低，為3.22，魚食性魚類的營(yíng)養(yǎng)級(jí)最高，為4.19。從圖2中可以看出，海州灣及鄰近海域生態(tài)系統(tǒng)營(yíng)養(yǎng)流動(dòng)以碎屑和牧食為基礎(chǔ)。該海域生態(tài)系統(tǒng)的能量主要來源于碎屑和浮游植物，其中47%來源于碎屑，53%來源于浮游植物。

評(píng)價(jià)本模型整體質(zhì)量的Pedigree 指數(shù)為0.660，與Morissette[32]等評(píng)價(jià)的其他50個(gè)不同的生態(tài)系統(tǒng)的Pedigree 指數(shù)（0.164～0.676）相比，處于較合理的范圍，表明本模型輸入?yún)?shù)的可靠性好，模型可信度較高。

表2 海州灣及鄰近海域生態(tài)系統(tǒng)Ecopath模型的基本參數(shù)Table2 Basic input data and estimated parameters for the Haizhou Bay and adjacent waters Ecopath model

3.2 營(yíng)養(yǎng)相互關(guān)系

圖2 海州灣及鄰近海域生態(tài)系統(tǒng)食物網(wǎng)Fig.2 Food web of the Haizhou Bay and adjacent waters

圖3 海州灣及鄰近海域生態(tài)系統(tǒng)各功能群之間的營(yíng)養(yǎng)關(guān)系Fig.3 Trophic relations of functional groups in the Haizhou Bay and adjacent waters ecosystem

海州灣及鄰近海域各功能群間的混合營(yíng)養(yǎng)效應(yīng)分析如圖3所示，藍(lán)色表示正影響，即該功能群生物量的增加有利于被影響功能群生物量的增加；紅色表示負(fù)影響，即該功能群生物量的增加造成被影響功能群生物量的減少；顏色的深淺表示影響程度的強(qiáng)弱。結(jié)果表明，各功能群生物量的增加加大了對(duì)餌料資源的種內(nèi)競(jìng)爭(zhēng)，均會(huì)對(duì)自身產(chǎn)生負(fù)影響；浮游植物和有機(jī)碎屑作為被捕食者，對(duì)大多數(shù)功能群產(chǎn)生正影響；浮游動(dòng)物和其他軟體動(dòng)物受到初級(jí)生產(chǎn)者和捕食者的雙重影響，在能量傳遞中發(fā)揮著關(guān)鍵作用，對(duì)生態(tài)系統(tǒng)的影響較強(qiáng)烈。其他蝦類生物量的增加對(duì)其他底棲動(dòng)物、戴氏赤蝦、其他蟹類和其他蝦/魚食性魚類產(chǎn)生顯著負(fù)影響，影響值分別為-0.449、-0.332、-0.317和-0.28；漁業(yè)捕撈對(duì)大部分漁業(yè)生物功能群有明顯負(fù)影響（圖3）。表3列出了該海域生態(tài)系統(tǒng)內(nèi)前10位的上行效應(yīng)與下行效應(yīng)，其中，大型頭足類、其他蝦/魚食性魚類、魚食性魚類對(duì)食物網(wǎng)產(chǎn)生的下行效應(yīng)分別為0.949、0.919、0.905；浮游植物、其他底棲動(dòng)物對(duì)食物網(wǎng)產(chǎn)生的上行效應(yīng)分別為0.967、0.825。

3.3 生態(tài)系統(tǒng)的總體特征

在Ecopath模型中，通過網(wǎng)絡(luò)分析功能計(jì)算的系統(tǒng)穩(wěn)定性和成熟度等多個(gè)生態(tài)系統(tǒng)參數(shù)，是評(píng)價(jià)生態(tài)系統(tǒng)結(jié)構(gòu)和功能的重要指標(biāo)。表4為現(xiàn)階段海州灣及鄰近海域生態(tài)系統(tǒng)特征參數(shù)表，其中系統(tǒng)總流量是衡量生態(tài)系統(tǒng)規(guī)模的重要指標(biāo)，是生態(tài)系統(tǒng)各功能群總輸出量、總呼吸量、總消耗量和流入碎屑總量的總和，海州灣及鄰近海域生態(tài)系統(tǒng)的總流量為4 790.691 t/（km2·a），總 消 耗 量 為 542.975 t/（km2·a）， 流 入 碎 屑 總 量 為2 045.675 t/（km2·a），分別占系統(tǒng)總流量的 11.34% 和42.71%，說明生態(tài)系統(tǒng)有較多能量并為被充分利用，而是轉(zhuǎn)入了系統(tǒng)再循環(huán)。系統(tǒng)的總能量轉(zhuǎn)換效率為12.63%，其中來自初級(jí)生產(chǎn)者能量流的平均轉(zhuǎn)換效率是12.33%，而來自碎屑能量流的平均轉(zhuǎn)換效率是13.20%。系統(tǒng)凈生產(chǎn)量為1891.732 t/（km2·a），總初級(jí)生產(chǎn)量為2202.041 t/（km2·a）?？偝跫?jí)生產(chǎn)量與總呼吸量的比值(Total Primary Production/Total Respiration,TPP/TR)是描述生態(tài)系統(tǒng)成熟度的關(guān)鍵指標(biāo)，海州灣及鄰近海域生態(tài)系統(tǒng)的TPP/TR為7.069。系統(tǒng)的連接指數(shù)（Connectance Index, CI）和系統(tǒng)雜食指數(shù)（System Omnivory Index, SOI）分別為 0.429 和 0.204，系統(tǒng)循環(huán)指數(shù)（Finn Cycling Index, FCI）為 1.392%（表 4）。

表3 海州灣及鄰近海域生態(tài)系統(tǒng)內(nèi)主要營(yíng)養(yǎng)效應(yīng)（前10位）Table3 The top ten trophic effects of the Haizhou Bay and adjacent waters ecosystems

4 討論

本研究基于Ecopath模型的構(gòu)建，從物質(zhì)能量平衡的角度分析了海州灣及鄰近海域生態(tài)系統(tǒng)的營(yíng)養(yǎng)結(jié)構(gòu)和生態(tài)系統(tǒng)特征，研究發(fā)現(xiàn)該生態(tài)系統(tǒng)功能群的營(yíng)養(yǎng)級(jí)范圍從1.00（浮游植物和碎屑）到4.19（魚食性魚類）。魚食性魚類具有最高的營(yíng)養(yǎng)級(jí)，攝食大量魚類，是海州灣食物網(wǎng)中的高營(yíng)養(yǎng)級(jí)捕食者[33]。而方氏云鳚則是海州灣食物網(wǎng)中營(yíng)養(yǎng)級(jí)最低的魚類，雙殼類在其食物中所占比例最高[33]。

表4 海州灣及鄰近海域生態(tài)系統(tǒng)的總體特征參數(shù)Table4 General characteristic parameters for the Haizhou Bay and adjacent waters ecosystem

根據(jù)生態(tài)系統(tǒng)成熟度及穩(wěn)定性理論[34]，在生態(tài)系統(tǒng)發(fā)育初期，大部分能量被用于物種的生長(zhǎng)發(fā)育，較多功能群的生產(chǎn)量大于呼吸量（TPP/TR＞1）。隨著生態(tài)系統(tǒng)逐步演替，不斷聚集生物量，用于維持呼吸作用的能量逐漸增多，凈生產(chǎn)量降低，在生態(tài)系統(tǒng)將要達(dá)到成熟狀態(tài)時(shí)，總生物量接近最大值，因此對(duì)于一個(gè)成熟的生態(tài)系統(tǒng)來說，TPP/TR趨近1，而TPP/B較小。海州灣及鄰近海域生態(tài)系統(tǒng)TPP/TR為7.096，明顯高于1，總初級(jí)生產(chǎn)量超過了總呼吸量，意味著該生態(tài)系統(tǒng)中有許多營(yíng)養(yǎng)物質(zhì)未被利用，TPP/B為56.866，表明該水域生態(tài)系統(tǒng)尚未發(fā)展成熟，與王騰等[14]、張碩等[15]對(duì)海州灣的研究結(jié)論一致。循環(huán)指數(shù)（FCI）代表生態(tài)系統(tǒng)被循環(huán)利用的程度，海州灣及鄰近海域生態(tài)系統(tǒng)的再循環(huán)率為1.392%。在生態(tài)系統(tǒng)演替發(fā)展過程中，物種逐漸增多，各物種間的捕食關(guān)系也由線性逐漸趨向網(wǎng)狀，食物網(wǎng)結(jié)構(gòu)更加復(fù)雜。CI和SOI是表示生態(tài)系統(tǒng)內(nèi)部連接復(fù)雜程度的指標(biāo)，值越大代表系統(tǒng)內(nèi)部越復(fù)雜，功能群之間的聯(lián)系就越強(qiáng)、越穩(wěn)定，系統(tǒng)抵抗外界環(huán)境變動(dòng)的能力越強(qiáng)?，F(xiàn)階段海州灣及鄰近海域生態(tài)系統(tǒng)的SOI和CI分別為0.204和0.429，王騰等[14]對(duì)該海域2013年生態(tài)系統(tǒng)的SOI和CI研究結(jié)果分別為0.174和0.415，兩者基本一致。另外，該結(jié)果與其他未成熟生態(tài)系統(tǒng)也相似，如膠州灣[10]、黃河口[35]。說明該系統(tǒng)的食物網(wǎng)結(jié)構(gòu)相對(duì)簡(jiǎn)單，系統(tǒng)處于不穩(wěn)定狀態(tài)，同時(shí)具有的凈生產(chǎn)量較高，為1891.732 t/（km·a）?？偟膩碚f，海州灣及鄰近海域生態(tài)系統(tǒng)食物網(wǎng)連接相對(duì)簡(jiǎn)單，容易受外界擾動(dòng)影響。分析其原因，可能是近年來海州灣及鄰近海域受捕撈和海洋開發(fā)活動(dòng)等因素的影響，生態(tài)環(huán)境變化較大，魚類多樣性降低，資源衰退[36]，導(dǎo)致食物網(wǎng)結(jié)構(gòu)簡(jiǎn)單。

目前，食物矩陣的分析結(jié)果大多由胃含物分析得來，而胃含物分析方法存在一些不足[37-39]，不能反映物種長(zhǎng)時(shí)期內(nèi)攝食的情況，因而會(huì)影響該模型的結(jié)果。輸入準(zhǔn)確的食物矩陣數(shù)據(jù)是提高模型準(zhǔn)確性的關(guān)鍵。隨著技術(shù)發(fā)展，穩(wěn)定同位素方法已經(jīng)成為一種分析食物網(wǎng)營(yíng)養(yǎng)結(jié)構(gòu)的重要手段[40]，這種方法能提供一段時(shí)間內(nèi)物種積累的營(yíng)養(yǎng)數(shù)據(jù)、食物來源以及物種在食物網(wǎng)中的位置[41]。因此，在接下來的研究中，需要結(jié)合穩(wěn)定同位素方法，調(diào)整食物矩陣，提高模型質(zhì)量。

通過對(duì)海州灣及鄰近海域生態(tài)系統(tǒng)模型進(jìn)行研究，能夠了解該海域的營(yíng)養(yǎng)結(jié)構(gòu)和系統(tǒng)發(fā)育狀況，有助于該海域生態(tài)系統(tǒng)的管理規(guī)劃，為漁業(yè)資源的可持續(xù)管理提供科學(xué)的理論指導(dǎo)。為了防止海州灣及鄰近海域漁業(yè)資源的進(jìn)一步衰退，實(shí)現(xiàn)漁業(yè)的可持續(xù)發(fā)展，提出以下幾點(diǎn)建議：（1）減少漁業(yè)燃料補(bǔ)貼，并增加個(gè)人轉(zhuǎn)行從事其他行業(yè)的補(bǔ)貼，降低捕撈強(qiáng)度；（2）加大在休漁期間非法捕魚的懲罰力度，與配額捕撈制度相結(jié)合，避免休漁期過后出現(xiàn)過度競(jìng)爭(zhēng)；（3）發(fā)展近海養(yǎng)殖業(yè)，增加海產(chǎn)品供應(yīng)，減輕沿海漁業(yè)資源的捕撈壓力。加大對(duì)該海域生態(tài)環(huán)境的保護(hù)，維護(hù)生態(tài)系統(tǒng)的健康和穩(wěn)定，最終實(shí)現(xiàn)海洋漁業(yè)的可持續(xù)發(fā)展。