劉旭東，于建釗，張曉紅，崔文連，汪進(jìn)生，孫立娥，王艷玲，賈勝華

1.青島市環(huán)境監(jiān)測中心站，山東 青島 266003 2.中國環(huán)境監(jiān)測總站，國家環(huán)境保護(hù)環(huán)境監(jiān)測質(zhì)量控制重點(diǎn)實(shí)驗(yàn)室，北京 100012 3.浙江省河海測繪院，浙江 杭州 310008

生產(chǎn)力屬于生態(tài)系統(tǒng)生態(tài)學(xué)的部分，其基本定義為“單位時間、單位面積或體積生物有機(jī)體有機(jī)物質(zhì) (或能量) 的增加量”[1-4]，一般分為初級生產(chǎn)力和次級生產(chǎn)力。次級生產(chǎn)力指的是動物和異養(yǎng)微生物通過生長和繁殖而增加的生物量或儲存的能量。從能量學(xué)的角度來看，次級生產(chǎn)力與攝食量、同化量、排泄量以及呼吸量存在著密切關(guān)系[5]。大型底棲動物是海洋環(huán)境中的一個重要的生態(tài)類群，它在海洋生態(tài)系的能流和物流中占有十分重要的地位。以下2個特點(diǎn)使得底棲動物生產(chǎn)力的研究受到關(guān)注：①與浮游生物相比，大型底棲動物由于移動能力弱，活動范圍有限，不易受短時或者隨機(jī)性的短暫環(huán)境變化影響，反而能夠更準(zhǔn)確地反映出所處環(huán)境的長期、宏觀變化[6]；②大型底棲動物生活在泥水界面交界處，取食浮游生物、底棲藻類和有機(jī)碎屑等，本身又被經(jīng)濟(jì)水產(chǎn)動物 (如魚類) 所食，此處產(chǎn)生大量的物質(zhì)和能量的流動，在整個底棲生態(tài)系統(tǒng)中的物質(zhì)循環(huán)和能量流動過程中起著承上啟下的作用[7]。因此研究海洋大型底棲動物的次級生產(chǎn)力，不僅可幫助人們了解大型底棲動物的功能變化，而且對研究生態(tài)系統(tǒng)的物流、能流過程和海洋生物資源的持續(xù)利用以及評價海洋環(huán)境壓力也具有重要意義[8-9]。

海洋大型底棲動物群落次級生產(chǎn)力的估算發(fā)展于20 世紀(jì) 70 年代，采用的多是一些經(jīng)典算法，如同生群法(包括減員累計法、增長累計法、瞬時增長法、Allen 曲線法)和非同生群法(體長頻率法)，但有些數(shù)據(jù)不易獲得，估算過程復(fù)雜費(fèi)力。BREY[9]根據(jù)WARWICK[10]的建議，采用對數(shù)轉(zhuǎn)換的多元線性回歸法，給出了Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ停瑢⑷菀撰@得的參數(shù)(如生物量、密度)與底棲動物的次級生產(chǎn)力聯(lián)系起來。TUMBIOLO等[11]曾運(yùn)用此公式計算了包括多毛類、軟體動物、甲殼類和棘皮動物在內(nèi)的34個不同海域大型底棲動物的次級生產(chǎn)力，驗(yàn)證了Brey(1990)模型能夠較好地擬合這些數(shù)據(jù)，但也同時指出其估計偏差與水溫和水深等環(huán)境變量密切相關(guān)。中國于2001年引入此模型[12]，目前國內(nèi)大型底棲動物次級生產(chǎn)力估算基本采用Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ蚚13-31]。此后國外學(xué)者又提出部分大型底棲動物次級生產(chǎn)力估算模型[11,32-36]，并在模型中添加了環(huán)境因子(如水溫、水深等參數(shù))。但這些模型或是建模數(shù)據(jù)范圍較小，不一定適用于低緯度海域[33]，或是采用生物壽命等不易確定的參數(shù)，難以應(yīng)用到實(shí)踐[32]。DOLBETH等[35]曾在2005年評估了在不同環(huán)境下幾種模型的估計效果，發(fā)現(xiàn)Brey(2001)[34]模型是最好的選擇。此外，BREY于2012年在 Brey(2001)估算模型的基礎(chǔ)上增加了大型底棲動物食性和移動性等易觀察的數(shù)據(jù)(共計20個參數(shù))，相較以前應(yīng)用的線性模型提出了改進(jìn)的人工神經(jīng)網(wǎng)絡(luò)模型Brey(2012)[36]，使估算結(jié)果更加真實(shí)可靠。國內(nèi)目前對其他模型的應(yīng)用不多，Brey(2001)模型[34]僅有幾篇相關(guān)研究成果[8,30,37]，Tumbiolo和 Downing 模型[11]報道更少[30]，關(guān)于Brey(2012)模型，筆者只查閱到一篇文獻(xiàn)[8]。

膠州灣是一個典型的半封閉海灣，在生態(tài)類型上更是高度人為干預(yù)下的復(fù)合生態(tài)系統(tǒng)的典型實(shí)例，膠州灣是中國海洋生態(tài)學(xué)研究開展最早的海域之一，對于膠州灣大型底棲動物次級生產(chǎn)力的研究報道始于2005年，迄今筆者只查閱到4篇相關(guān)文獻(xiàn)[13-14,24,38]，文獻(xiàn)中均采用了Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ偷挠嬎惴椒?，通過梳理發(fā)現(xiàn)存在3個問題：①在其計算過程中，不論是運(yùn)用逐種累計法還是位點(diǎn)平均計算法，采用的均是Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ偷钠骄鶇?shù)，并未考慮底棲動物類群分類，應(yīng)該根據(jù)底棲動物不同的類群選擇相應(yīng)的參數(shù)，這樣更符合Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ偷挠嬎惴绞?，估計結(jié)果也更加準(zhǔn)確，目前國內(nèi)已有部分研究采用這種更準(zhǔn)確的計算方式[8,27,30-31];②次級生產(chǎn)力可以用濕重、干重、無灰分干重、碳含量、氮含量或能量含量等來表示,但從理論上來說，按能量測定是最真實(shí)、最合理、最具有可比性的，特別是在進(jìn)行生態(tài)系統(tǒng)的能流分析時更應(yīng)如此[5]，對于膠州灣次級生產(chǎn)力的報道，目前均是采用的去灰干重的含量來表示，未能換算到能量含量水平;③Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ褪窃缙诘难芯拷Y(jié)果，其公式簡單，僅包含生物量和生物密度2個參數(shù)，并沒有包含環(huán)境及底棲動物類群的功能性特點(diǎn)。近年來無論是研究的深度還是廣度均與國外研究存在較大差距[39]。基于以上3個原因，筆者根據(jù)2015—2017年連續(xù)3年對膠州灣大型底棲動物的監(jiān)測，選取Brey(1990) 經(jīng)驗(yàn)?zāi)Ｐ蚚9]，Brey(2001)模型[34，40]以及Brey(2012)模型[36,40]分別估算了膠州灣海域大型底棲動物群落的次級生產(chǎn)力，并與膠州灣的歷史數(shù)據(jù)進(jìn)行對比，以期為該海域大型底棲動物的功能研究及資源利用提供參考。

1 調(diào)查與分析方法

1.1 研究區(qū)域與采樣點(diǎn)布設(shè)

研究區(qū)域范圍為膠州灣海域，避開了灣口處于航道位置及灣口北部中沙礁處附近不易采集樣品的區(qū)域，共布設(shè)了12個采樣點(diǎn)(圖1)，于2015—2017年每年調(diào)查春、夏、秋3個航次，調(diào)查時間分別為每年的4月下旬至5月、7月下旬至8月、10月。其中，2015年夏季航次由于潮水問題，QD07及QD08站未能成功采集樣品，研究總計使用3年9個航次共106個采樣點(diǎn)的大型底棲動物采樣數(shù)據(jù)。采樣使用取樣面積為0.05 m2的蚌式采泥器，每個采樣點(diǎn)取樣4次，用0.5 mm孔徑的網(wǎng)篩分選樣品。樣品的處理、保存、計數(shù)等均按《海洋調(diào)查規(guī)范》[12]操作，為避免混淆，研究所描述的生物量如無特殊說明均為濕重生物量(WM)。

圖1 采樣點(diǎn)設(shè)置示意圖Fig.1 Sketch map of sampling stations

1.2 數(shù)據(jù)處理及分析

1.2.1 Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ蚚9]

lgP=a+b1×lgB-b2×lgW

式中：P為年均大型底棲動物次級生產(chǎn)力，g/(m2·a)；B為年均去灰干重生物量(濕重的9/50作為去灰干重，即濕重的1/5作為干重，干重的9/10作為去灰干重，該公式所描述的生物量按去灰干重計)，g/m2；W為年均個體重量，g/ind.。不同類群的系數(shù)見表1。為了最終和另外2個模型結(jié)果對比，將該計算單位換算為kJ/(m2·a)，換算比例為1 g去灰干重=0.518 g碳，1 g碳=45.7 kJ能量[41]。同時，為了和膠州灣歷史調(diào)查數(shù)據(jù)對比，采用了平均(表1中其他類)參數(shù)另行計算一次。

表 1 Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ椭胁煌惾旱?3 個系數(shù)Table 1 Three coefficients for different phyla in Brey′s empirical formula

1.2.2 Brey(2001)模型[34,40]

該模型目前為4-04版本，于2004年發(fā)布，相較2001模型微調(diào)了部分參數(shù)值，具體公式為

Lg(P/B)=7.947-2.294 lg(M)-2 409.856/(T+273)+0.168/D+0.194 SubT+0.180 InEpi+0.277 MoEpi+0.174 Taxon1-0.188 Taxon2+0.33Taxon3-0.062 Habitat1+582.851 lg(M)/(T+273)

式中：P為大型底棲動物年均次級生產(chǎn)力，kJ/(m2·a)；M為年平均個體體重能值，kJ；T為平均底層水溫，℃；D為平均水深，m；B為年平均生物量對應(yīng)的能值，kJ/m2；SubT:潮下帶值為1，潮間帶為0；InEpi:底內(nèi)生物值為1，底表生物為0；MoEpi:移動生物值為1，固著生物為0；Taxon1:環(huán)節(jié)動物或甲殼動物值為1，其他類群為0；Taxon2：棘皮動物值為1，其他類群為0；Taxon3:昆蟲值為1，其他類群為0；Habitat1:湖泊生境值為1，其他生境值為0。

1.2.3 Brey(2012)模型[36,40]

Lg(P/B)=a0+a1×H1+a2×H2

H1=TanH(b0+b1×lg(M)+b2×1/T+b3×lg(D)+b4×Mollusca…+b20×Exploited)

H2=TanH(c0+c1×lg(M)+c2×1/T+c3×lg(D)+c4×Mollusca…+c20×Exploited)

式中：P為大型底棲動物年均次級生產(chǎn)力，kJ/(m2·a)；B為年平均生物量對應(yīng)的能值,kJ/m2；M為年平均個體體重能值,J；T為平均底層水溫,℃；D為平均水深,m；b4～b8(c4～c8)為生物類群參數(shù)，分別對應(yīng)軟體動物、環(huán)節(jié)動物、甲殼類、棘皮動物及昆蟲類；b9～b12(c9～c12)為棲息類型及運(yùn)動方式參數(shù)，分別對應(yīng)是否底內(nèi)生物、是否固著生物、是否只爬行、是否兼性游泳生物；b13～b15(c13～c15)為食性參數(shù)，分別對應(yīng)植食、雜食及肉食動物；b16～b18(c16～c18)為棲息環(huán)境參數(shù)，分別對應(yīng)湖泊、河流、海洋；b19(c19)為是否為潮間帶，b20(c20)為魚類種群，通常默認(rèn)設(shè)置為0。

1.2.4 大型底棲動物能值的轉(zhuǎn)換關(guān)系

根據(jù)模型對數(shù)據(jù)的要求，需將生物量濕重數(shù)據(jù)轉(zhuǎn)換為能值，此項(xiàng)數(shù)據(jù)來自于一個關(guān)于水生生物身體成分的全球數(shù)據(jù)庫[40,42]。它涵蓋了體重(濕、干、無灰干重)，身體成分(蛋白質(zhì)、脂質(zhì)、碳水化合物)，元素(C、N、P)和能量含量(J)之間的比率，目前為Conversion04版本，涵蓋了共計3 444種生物的5 678個數(shù)據(jù)信息。

依據(jù)膠州灣3年調(diào)查的生物種類情況， 生物類群、運(yùn)動方式、食性等參數(shù)輸入數(shù)據(jù)庫中計算得到該次分析需要用到的轉(zhuǎn)換關(guān)系(表2)。其中未能檢索到螠蟲動物門的多皺無吻螠(Arhynchiterugosum)及節(jié)肢動物門的日本長尾蟲(Aspeudesnipponicus)的信息，由于此2種生物在調(diào)查中出現(xiàn)率低，且并非膠州灣內(nèi)密度及生物量占據(jù)優(yōu)勢地位的類群，故采用與其分類及演化地位相近的星蟲動物門和等足目的基準(zhǔn)能量值來替代。

表2 大型底棲動物濕重生物量和能值的轉(zhuǎn)換關(guān)系Table 2 Conversion factors of the wet weight and energy of macrozoobenthos

1.2.5 數(shù)據(jù)分析方法

事先將生產(chǎn)力結(jié)果進(jìn)行l(wèi)g(X+1)轉(zhuǎn)換，并通過了正態(tài)分布和方差齊性檢驗(yàn)。不同模型計算的生產(chǎn)力結(jié)果比較采用配對t檢驗(yàn)，不同季節(jié)生產(chǎn)力的比較采用方差分析(ANOVA)。分析結(jié)果由SPSS 19.0完成，采樣點(diǎn)及平面分布圖由ArcGis10.2軟件繪制。

2 結(jié)果

2.1 種類組成情況

歷次調(diào)查共獲11門19綱41目111科共191種底棲動物(附錄A)，其中，環(huán)節(jié)動物門多毛類出現(xiàn)的種類數(shù)最多，共46種，占底棲動物總種類組成的36.8%；其次為節(jié)肢動物甲殼類，共38種，占30.4%；軟體動物共出現(xiàn)24種，占19.2%；棘皮動物出現(xiàn)8種，占6.4%；其余門類種數(shù)較少，合計占7.2%。3年調(diào)查中，生物種類出現(xiàn)數(shù)每年基本一致，各季節(jié)種類數(shù)有所不同，春夏季生物種類較多，秋季較少(圖2)。

圖2 調(diào)查海域的大型底棲動物組成Fig.2 Composition of each phylum of macrozoobenthos in Jiaozhou Bay

2.2 不同類群的密度及生物量

從調(diào)查數(shù)據(jù)來看，膠州灣不同調(diào)查航次的大型底棲動物類群的密度、生物量分布相差較大(表3，圖3和圖4)，春、夏、秋季大型底棲動物的平均密度分別為595.3、505.3、277.5個/m2，平均密度為459.4個/m2，多毛類和軟體動物占據(jù)密度組成主要部分，平均密度分別為200.1、166.8個/m2；春、夏、秋季平均生物量分別為1 078.5、1 011.9、394.8 g/m2，平均生物量為828.4 g/m2，軟體動物是占據(jù)生物量組成的主要部分，平均生物量為786.8 g/m2。膠州灣大型底棲動物的密度和生物量均呈現(xiàn)春、夏季較高，秋季低的趨勢，在不同季節(jié)中，密度變化與多毛類(r=0.705，P=0.034)和軟體動物(r=0.794，P=0.011)的密度相關(guān)，生物量的變化則與軟體動物(主要為菲律賓蛤仔)的生物量變化有關(guān)(r=0.996，P<0.001)。

表3 膠州灣大型底棲動物主要類群密度和生物量Table 3 The density and biomass of main groups of macrozoobenthos at different investigations

圖3 大型底棲動物主要類群生物密度的變化分布Fig.3 The density variation of main groups of macrozoobenthos at different investigations

圖4 大型底棲動物主要類群生物量的變化分布Fig.4 The biomass variation of main groups of macrozoobenthos at different investigations

2.3 3種模型的次級生產(chǎn)力估計

去灰干重計，Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ头诸惾簠?shù)估算膠州灣海域年均次級生產(chǎn)力為43.20 g/(m2·a)，換算后為1 022.66 kJ/(m2·a)；Brey(2001)模型估算膠州灣海域年均次級生產(chǎn)力為837.81 kJ/(m2·a)；Brey(2012)模型估算膠州灣海域年均次級生產(chǎn)力為986.84 kJ/(m2·a)。各采樣點(diǎn)3種估算方法的分布趨勢基本一致(圖5)。配對t檢驗(yàn)分析結(jié)果顯示,Brey(2012)與Brey(2001)模型估算結(jié)果并無顯著性差異(t=0.397，P=0.692)；Brey(2012)與Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ凸浪憬Y(jié)果差異顯著(t=5.561，P<0.001)；Brey(2001)與Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ凸浪憬Y(jié)果差異顯著(t=6.067，P<0.001)。

此外，采用Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ推骄鶇?shù)估算膠州灣海域年均次級生產(chǎn)力，結(jié)果為54.71 g/(m2·a)，高于按類群估算為43.20 g/(m2·a)的結(jié)果(t=9.690，P<0.001)。

圖5 不同估算方法的膠州灣大型底棲動物次級生產(chǎn)力的變化趨勢Fig.5 The change in the secondary production of macrozoobenthos among different models

3 討論

3.1 不同海域次級生產(chǎn)力的比較

在進(jìn)行底棲動物群落歷史對比分析時,由于以往調(diào)查者的目的和側(cè)重點(diǎn)有所差異, 調(diào)查采樣點(diǎn)和監(jiān)測水平不一致，調(diào)查結(jié)果與采樣航次、采樣點(diǎn)是否均勻覆蓋、季節(jié)性有很大關(guān)系。該研究結(jié)果與對膠州灣次級生產(chǎn)力的歷史研究[13-14,24,38]結(jié)果整體趨勢一致，即膠州灣大型底棲動物次級生產(chǎn)力由灣口到灣北部的灣頂呈現(xiàn)升高趨勢，菲律賓蛤仔是次級生產(chǎn)力的主要貢獻(xiàn)者，菲律賓蛤仔資源的貢獻(xiàn)使得膠州灣大型底棲動物的次級生產(chǎn)力遠(yuǎn)高于其他海域的次級生產(chǎn)力(表4)。

表4 膠州灣大型底棲動物平均次級生產(chǎn)力與其他海域的比較Table 4 The secondary production of macrozoobenthos in Jiaozhou Bay and the comparison with other seas

研究所調(diào)查的膠州灣海域大型底棲動物的次級生產(chǎn)力較歷史資料高有以下幾個原因：①調(diào)查的采樣點(diǎn)均分布于膠州灣內(nèi)，并避開了水深較深的航道區(qū)域及不易采樣的中砂礁位置，與其他調(diào)查航次采樣點(diǎn)相比，處于其高值區(qū)[13,38]；②次級生產(chǎn)力與水溫密切相關(guān)，升高溫度可加快底棲動物的生長發(fā)育速度，進(jìn)而提高生產(chǎn)力[5]。研究未能調(diào)查冬季航次，因此在估算年均生產(chǎn)力時可能會導(dǎo)致結(jié)果偏高，推薦根據(jù)實(shí)際情況采用按季節(jié)估算的結(jié)果；③菲律賓蛤仔是次級生產(chǎn)力的主要貢獻(xiàn)者,該次調(diào)查中有很多采樣點(diǎn)位于菲律賓蛤仔養(yǎng)殖區(qū)，可采到大量的菲律賓蛤仔(如在J6采樣點(diǎn)，2016年春季航次中菲律賓蛤仔的密度達(dá)到1 060個/m2)，因此模型的估計結(jié)果受菲律賓蛤仔資源的波動影響。

3.2 合理模型的選擇

經(jīng)驗(yàn)?zāi)Ｐ椭皇菍θ郝浯渭壣a(chǎn)力的估算，并非真實(shí)值，研究所采用的3個模型逐步復(fù)雜精細(xì)化，這與近年來對于次級生產(chǎn)力研究的發(fā)展相關(guān)。其差異主要表現(xiàn)在2個方面：①模型考慮的參數(shù)由2個增加到20個(如增加了底棲動物的食性、運(yùn)動方式情況)，同時增加了環(huán)境因子;②能值轉(zhuǎn)換以前所有底棲動物均采用一個轉(zhuǎn)化基準(zhǔn)，隨著研究的進(jìn)一步深入，現(xiàn)在精細(xì)到更細(xì)的類群專屬轉(zhuǎn)換比例。研究顯示所用的3種經(jīng)驗(yàn)?zāi)Ｐ偷玫降拇渭壣a(chǎn)力估算存在差異，Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ凸浪憬Y(jié)果偏高。而關(guān)于Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ团cBrey(2001)模型估算結(jié)果差異性顯著的情況，均在國內(nèi)的3篇報道[8,30,37]中出現(xiàn)，一般認(rèn)為Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ凸街形瓷婕碍h(huán)境參數(shù)，實(shí)質(zhì)上僅給出特定環(huán)境條件下群落次級生產(chǎn)力的平均水平，這可能使 Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ凸綄Σ煌竟?jié)和不同環(huán)境的次級生產(chǎn)力估算效果不理想[8,30,37]。此外在膠州灣海域中采用Brey(1990)經(jīng)驗(yàn)?zāi)Ｐ筒话搭惾簞澐钟嬎愕那闆r，其結(jié)果高于分類的計算結(jié)果，因?yàn)檐涹w動物菲律賓蛤仔是膠州灣大型底棲動物次級生產(chǎn)力的主要貢獻(xiàn)者，其計算系數(shù)2種方式有較大差別，在該研究中不按類群劃分會導(dǎo)致計算結(jié)果高估26.7%。

從建模方法看，Brey(1990) 經(jīng)驗(yàn)?zāi)Ｐ蚚9]與Brey(2001)模型[34，40]采用的為多元線性回歸法(MLM)，而Brey(2012)模型[36,40]采用人工神經(jīng)網(wǎng)絡(luò)法(ANN)。MLM對獨(dú)立變量間的交互性非常敏感，而ANN可以對自變量和因變量間復(fù)雜的、非線性的、間斷性的關(guān)系建模。大型底棲動物群落次級生產(chǎn)力估計上，MLM法優(yōu)于ANN法[36]。因此，雖然3種模型估算的膠州灣大型底棲動物次級生產(chǎn)力的時空變化趨勢基本一致，Brey(2001)模型與Brey(2012)模型的估算結(jié)果無顯著性差異，該研究仍推薦采用Brey(2012)模型進(jìn)行膠州灣大型底棲動物次級生產(chǎn)力的估計，并采用Brey(2012)模型的計算結(jié)果來討論膠州灣大型底棲動物次級生產(chǎn)力的時空分布。

3.3 膠州灣大型底棲動物次級生產(chǎn)力的時空分布

由于存在影響生態(tài)系統(tǒng)的多重壓力因素(如富營養(yǎng)化、污染、棲息地喪失、氣候變化、入侵物種引進(jìn)等)，可以通過對種群的靜態(tài)和動態(tài)分析檢測到這些變化，次級生產(chǎn)力在實(shí)際應(yīng)用中提供的信息遠(yuǎn)比單獨(dú)的密度或生物量等其他靜態(tài)測量方法所提供的信息多，它結(jié)合了種群生態(tài)性能的靜態(tài)和動態(tài)組成部分，并不僅僅服務(wù)于能量流動和利用效率的研究[43]，如圖6所示，2017年春季具有較高的密度和生物量水平，但生產(chǎn)力卻處于平均水平；2015年夏季生物量處于一般水平，但生產(chǎn)力卻處于最高水平；2015年春季密度處于較高水平，但生產(chǎn)力卻處于低水平。

注：得分為1～0，分別代表每個參數(shù)在3年研究中的最高和最低值。圖6 不同調(diào)查時期的密度、生物量及次級生產(chǎn)力的得分Fig.6 Scores for the secondary production, density and biomass for the macrozoobenthos of Jiaozhou Bay from 2015 to 2017

膠州灣大型底棲動物具有較高的生產(chǎn)力水平(一般也意味著高經(jīng)濟(jì)水平)，但是高水平的生產(chǎn)力并不一定代表一個更健康的生態(tài)系統(tǒng)，一些影響可能導(dǎo)致整體底棲動物產(chǎn)量的增加，因?yàn)榄h(huán)境的變化有利于高產(chǎn)機(jī)會主義物種，但同時也降低了種群的復(fù)雜性(降低均勻度和豐富度)[43]，正確認(rèn)識膠州灣大型底棲動物的次級生產(chǎn)力來源對研究膠州灣海域的大型底棲動物群落生態(tài)環(huán)境狀況也具有重要意義。通過該研究發(fā)現(xiàn)，膠州灣大型底棲動物各類群所占比例相比于歷史數(shù)據(jù)并沒有發(fā)生太大變化，該研究的結(jié)果是環(huán)節(jié)動物多毛類(36.8%)>節(jié)肢動物甲殼類(30.4%)>軟體動物(19.2%)；徐兆東[24]2014年研究結(jié)果為環(huán)節(jié)動物多毛類(41.8%)>節(jié)肢動物甲殼類(29.3%)>軟體動物(20.5%)；于海燕等[44]1998—2001年的研究結(jié)果為多毛類(41.3%)>節(jié)肢動物甲殼類(28.6%)>軟體動物(18.3%)，孫道遠(yuǎn)等[45]1991—1994年的研究結(jié)果為多毛類(39.0%)>節(jié)肢動物甲殼類(28.5%)>軟體動物(20.0%)。對比發(fā)現(xiàn)膠州灣歷次調(diào)查大型底棲動物類群順序相同，所占比例極其接近，說明膠州灣近30年的大型底棲動物在類群組成上并沒有發(fā)生明顯變化，處于穩(wěn)定狀態(tài)。膠州灣大型底棲動物主要的次級生產(chǎn)力貢獻(xiàn)來源于菲律賓蛤仔。

膠州灣海域大型底棲動物次級生產(chǎn)力呈現(xiàn)出夏季>春季>冬季的變化趨勢(表5)，季節(jié)間差異顯著(F2,36=4.351，P=0.021)，春、夏、秋季的估算結(jié)果分別為880.17、1 537.64、458.33 kJ/(m2·a)。主要是夏季水溫較高，菲律賓蛤仔較活躍，從而具備較高的生物量和密度。在空間分布上(圖7)，高次級生產(chǎn)力的區(qū)域?yàn)樘幱跒稠斨胁考t島以南到向西偏移至大沽河口和紅石崖附近，這也與以往研究一致[13-14,24,38]?？傮w來看，膠州灣大型底棲動物次級生產(chǎn)力具備較高的水平，且具有顯著的時空變化，菲律賓蛤仔資源的變化是造成波動的最主要因素。

表5 不同季節(jié)膠州灣大型底棲動物的次級生產(chǎn)力Table 5 The secondary production of macrozoobenthos in Jiaozhou Bay at different seasons

圖7 膠州灣大型底棲動物的次級生產(chǎn)力分布Fig.7 The secondary production distribution of macrozoobenthos in Jiaozhou Bay

4 結(jié)論

膠州灣近30年的大型底棲動物在類群組成上并沒有發(fā)生明顯變化，處于穩(wěn)定狀態(tài)。膠州灣大型底棲動物次級生產(chǎn)力由灣口到灣北部的灣頂呈現(xiàn)升高趨勢，菲律賓蛤仔是次級生產(chǎn)力的主要貢獻(xiàn)者，也導(dǎo)致膠州灣大型底棲動物的次級生產(chǎn)力遠(yuǎn)高于其他海域。

相較其他估算模型，推薦使用Brey(2012)模型估算膠州灣海域大型底棲動物的次級生產(chǎn)力，膠州灣海域春、夏、秋季的估算結(jié)果分別為880.17、1 537.64、458.33 kJ/(m2·a)。