海洋大型底棲動物群落次級生產(chǎn)力估算模型研究進展及應用評價

賈勝華，壽鹿，廖一波，曾江寧，高愛根，湯雁濱

（國家海洋局第二海洋研究所國家海洋局海洋生態(tài)系統(tǒng)與生物地球化學重點實驗室，浙江杭州310012）

海洋大型底棲動物群落次級生產(chǎn)力估算模型研究進展及應用評價

賈勝華，壽鹿，廖一波，曾江寧，高愛根，湯雁濱

（國家海洋局第二海洋研究所國家海洋局海洋生態(tài)系統(tǒng)與生物地球化學重點實驗室，浙江杭州310012）

大型底棲動物是海洋生態(tài)系統(tǒng)的重要組成部分，在生態(tài)系統(tǒng)物流和能流中占有重要地位。定量估算海洋大型底棲動物次級生產(chǎn)力在其功能研究中具有重要意義。本文綜述了海洋大型底棲動物次級生產(chǎn)力估算模型研究的3個階段：起步階段（1979-1990）、發(fā)展與完善（1990-2001）、成熟并廣泛應用（2001-至今）。國內關于海洋大型底棲動物次級生產(chǎn)力估算基本都采用Brey（1990）經(jīng)驗公式，采用經(jīng)驗公式的方式可分為3類：逐種計算、按站位計算和按類群計算，逐種和按類群計算最符合Brey（1990）給出的參考步驟。不同方法估算同一海域次級生產(chǎn)力結果不同，同一估算方法估算不同生境不同群落次級生產(chǎn)力結果偏差也不同。今后的工作應注重種群次級生產(chǎn)力研究，積累大量基礎數(shù)據(jù)，以便建立適宜我國特定海域的大型底棲動物次級生產(chǎn)力估算模型。

海洋大型底棲動物；群落次級生產(chǎn)力；估算模型；杭州灣潮間帶

次級生產(chǎn)力指動物通過生長、繁殖而增加的生物量或儲存的能量，是衡量生物生產(chǎn)能力的量值（Brey et al，1990）。海洋大型底棲動物群落次級生產(chǎn)力指單位時間內（通常為1年）研究海域內所有大型底棲動物通過生長和繁殖而使其有機物質或能量的增加量。一般用干重（dry weight，DW）、去灰干重（ash free dry weight，AFDW）或能值表示，單位分別為g（DW）/（m2·a）、g（AFDW）/（m2·a）和kJ/（m2·a）。影響海洋大型底棲動物群落次級生產(chǎn)力的主要生物因素有大型底棲動物的種類組成、棲息密度和生物量等，環(huán)境因素有水深、底溫、溶解氧、葉綠素a等（Brey et al，1998；李新正等，2005；嚴娟等，2012）。研究海洋大型底棲動物群落的次級生產(chǎn)力，不僅對研究底棲動物落自身變化機制具有重要意義（李新正等，2005），而且可為深入研究海洋生態(tài)系統(tǒng)初級、次級、終極生產(chǎn)力之間的物質循環(huán)與能量流動的動態(tài)變化規(guī)律和海洋生物資源的持續(xù)利用（Tumbiolo et al，1994）以及評價海洋環(huán)境壓力和生態(tài)系統(tǒng)健康（Dolbeth et al，2005；韓慶喜等，2011；劉曉收等，2014）提供幫助。

1 研究進展

海洋大型底棲動物群落次級生產(chǎn)力估算研究可分為三個階段。第一階段：起步階段（1979-1990），以Brey（1990）經(jīng)驗公式發(fā)表為標志；第二階段：發(fā)展與完善階段（1990-2001），以Brey（2001）估算模型建立為標志；第三階段：成熟并廣泛應用階段（2001-至今），以國內學者（于子山等，2001）首次引入經(jīng)驗公式對渤海大型底棲動物群落次級生產(chǎn)力進行初步研究為標志。

1.1 起步階段（1979-1990）

海洋大型底棲動物群落次級生產(chǎn)力估算起始于20世紀70年代，但在90年代以前，采用的多是一些經(jīng)典算法如同生群法（包括減員累計法、增長累計法、瞬時增長法、Allen曲線法）和非同生群法（體長頻率法），結果雖較為準確，但估算過程耗時費力、復雜漫長。人們希望建立次級生產(chǎn)力與生物量、棲息密度、溫度、水深等易測定參數(shù)的關系式（Waters et al，1977），但關系式是否簡單而可靠，并能準確估算群落次級生產(chǎn)力（Robertson et al，1979），Mann（1969）認為對大多數(shù)生物而言P/B值與生活史關系密切。隨后，陸續(xù)出現(xiàn)P/B值與生物生活史關系的報道（Brey et al，1990）。Robertson（1979）在前人對海洋大型底棲動物種群次級生產(chǎn)力研究的基礎上，基于物種生活史首次提出估算海洋大型底棲動物次級生產(chǎn)力的公式：

其中，P為大型底棲動物次級生產(chǎn)力，單位為g（DW）/（m2·a）；B為年平均干重生物量，單位為g（DW）/m2；L為生物壽命。

但物種生活史特征數(shù)據(jù)較難獲取，而且野外生物壽命更多取決于其他生物捕食或人類過度開發(fā)造成的死亡而不是絕對的最大壽命；受環(huán)境影響，某些一年多代的生物，其壽命為變量。因此該公式實踐中很難應用（Brey et al，1990；Tumbiolo et al，1994；Dolbeth et al，2005；Cusson et al，2005）。于是，Schwinghamer（1986）提出P/B值與平均個體體重的關系式：

其中，W為平均個體體重能值，單位為kJ/ind。

相對生物壽命，體重數(shù)據(jù)較穩(wěn)定且易測定（Tumbiolo et al，1994；Cusson et al，2005）。以上是早期的經(jīng)驗模型，P/B值與單一參數(shù)如生物壽命、平均個體體重之間的關系。Brey（1990）和Edgar（1990）建立了兩個參數(shù)的經(jīng)驗模型。Edgar（1990）基于41種大型底棲動物生產(chǎn)力數(shù)據(jù)資料提出P（μg/day）與去灰干重生物量B（μg）、水溫T（℃）之間的關系式：

但該公式并沒有得到廣泛應用。Brey（1990）根據(jù)Warwick（1980）“群落次級生產(chǎn)力也許可以通過經(jīng)驗關系估算，而不需要估算每個種群的生產(chǎn)力”的建議，采用對數(shù)轉換的多元線性回歸法，對來自北緯和南緯12°-70°、水深300 m以內海域大型底棲動物P、P/B與B、W之間的關系進行了分析，并將生產(chǎn)力的原始數(shù)值與通過經(jīng)驗公式的估計值比較，結果顯示：除去生物量較小的群落，采用經(jīng)驗公式估算次級生產(chǎn)力的偏差不超過20%。認為在群落水平上可以借助P和B、W之間的經(jīng)驗公式來估算群落次級生產(chǎn)力。群落次級生產(chǎn)力估算的經(jīng)驗公式：

其中：P為大型底棲動物年次級生產(chǎn)力，單位為g（AFDW）/（m2·a）；B為大型底棲動物年平均去灰干重，單位為g（AFDW）/m2；W為大型底棲動物個體的年平均去灰干重，單位為g（AFDW）/ ind。經(jīng)驗公式對淺海（水深小于300 m）的大型底棲動物群落生產(chǎn)力粗略估算很有用，準確性與物種組成有關，種類數(shù)越多，估算結果會越好。Tumbiolo等（1994）用來自全球34個不同海域、125種海洋無脊椎動物的次級生產(chǎn)力資料驗證了該公式，結果表明Brey（1990）經(jīng)驗公式與大型底棲動物群落實際次級生產(chǎn)力吻合度較高。至此，一個既實用又能較為準確估算海洋大型底棲動物群落次級生產(chǎn)力的公式誕生了。

1.2 發(fā)展與完善階段（1990-2001）

人們對淡水大型底棲動物次級生產(chǎn)力的研究總是走在海洋前面。法國科研人員發(fā)現(xiàn)：河流與湖泊的底棲生物生產(chǎn)力與年最大生物量的比例約為3∶2，水溫、水深、溶解氧等環(huán)境因子會影響大型底棲動物的生產(chǎn)力（龔志軍等，2001）。

Morin和Bourassa（1992）提出溪流中底棲動物次級生產(chǎn)力的估算公式：

其中，P為大型底棲動物次級生產(chǎn)力，單位為g（DW）/（m2·a）；B為年平均干重生物量，單位為g（DW）/m2；M為平均個體質量（g DW），單位為g（DW）/ind；T為底層水溫（℃）。

Benke（1993）在Plant等（1989）對淡水無脊椎生物研究的基礎上提出河流中底棲生物次級生產(chǎn)力的估算公式：

其中T為溫度，W為個體平均重量。以上是考慮3個參數(shù)的估算模型。溫度通過影響動物代謝速率改變次級生產(chǎn)力，水深通過影響食物質量改變次級生產(chǎn)力，并且相同的水深，食物質量因水溫不同而產(chǎn)生差異（Tumbiolo et al，1994），水深對次級生產(chǎn)力的影響在海洋中可能高于淡水（Plante et al，1989）。Tumbiolo等（1994）利用海洋大型底棲動物次級生產(chǎn)力的數(shù)據(jù)檢驗了Plante等（1989）淡水估算模型與Brey（1990）經(jīng)驗公式，認為淡水估算模型不適宜海洋次級生產(chǎn)力的估算，而Brey（1990）經(jīng)驗公式雖然估算效果較好但其偏差與水溫和水深密切相關，在一些溫度極高或極低和一些水深極淺或極深的海域，Brey（1990）經(jīng)驗公式估算值可達到超過10倍的偏差。并參考該淡水估算模型首次將環(huán)境因子考慮到海洋次級生產(chǎn)力的估算中，建立了一個適用范圍更廣的估算模型：

其中，P為大型底棲動物次級生產(chǎn)力，單位為g（DW）/（m2·a）；B為年平均干重生物量，單位為g（DW）/m2；Tb為底層水溫（℃）；Z為平均水深（單位為m，潮間帶水深為0 m）；Wm為群落最大個體重量，單位為mg（DW）。該模型中濕重、干重、去灰干重和能值的轉換關系見表1（Brey et al，2001；張崇良等，2011；紀瑩璐等，2015）。

表1 大型底棲動物重量轉換關系

但Tumbiolo等（1994）估算模型的建模數(shù)據(jù)主要集中于北緯50°-60°的歐洲海域，不一定適合低緯度的海域。有研究表明，該模型估算結果偏低（Dolbeth et al，2005；張崇良等，2011）。Sprung（1993）提出如果樣品中物種數(shù)少時，經(jīng)驗公式的估算誤差可能很大，而且由于大型底棲動物的種類及所處環(huán)境溫度不同，其生物生長和代謝速率不同，影響次級生產(chǎn)力，因此次級生產(chǎn)力估算也要考慮其種類組成和底層水溫。Brey（2001）提出隨著樣本物種數(shù)增多，經(jīng)驗公式估算誤差會迅速下降，回答了Sprung（1993）提出的問題，同時改進了經(jīng)驗公式，考慮了水深、水溫、大型底棲動物種類組成以及與底質的關系提出群落次級生產(chǎn)力的估算模型，見公式（8）和表2。Brey（2001）估算模型不適合估算單個種群的生產(chǎn)力（偏差在-60%~150%），幾個種群的估算值合并偏差會迅速減?。˙rey et al，2001）。該模型被認為是目前估算海洋大型底棲動物次級生產(chǎn)力最好的模型（Dolbeth et al，2005；Cusson et al，2005）。

表2 Brey（2001）估算模型參數(shù)

1.3 成熟并廣泛應用階段（2001-至今）

在Brey（1990）經(jīng)驗公式提出10年后，國內學者（于子山等，2001）首次采用該經(jīng)驗公式對渤海大型底棲動物群落次級生產(chǎn)力進行初步研究，并將大型底棲動物個體的年平均去灰干重（W）轉換為棲息密度（A）和以去灰干重表示的生物量（B）：

國內學者在引入Brey（1990）經(jīng)驗公式估算大型底棲動物群落次級生產(chǎn)力上可以分為三類：第一類是逐種計算（于子山等，2001；王金輝等，2006；金亮等，2007；袁偉等，2007；金亮等，2009；王宗興等，2011；杜永芬等，2012；方平福等，2013；吳辰等，2013；紀瑩璐等，2015），第二類是按站位（或斷面或樣點）計算（李新正等，2005；呂小梅等，2008；周進等，2008；周細平等，2008；劉勇等，2008；方少華等，2009；張海波等，2010；焦海峰等，2011；周福芳等，2012；周進等，2012；嚴娟等，2012），第三類是按類群計算（杜飛雁等，2008；田勝艷等，2010；王志忠等，2012；蔡立哲等，2012；劉坤等，2015）。

逐種計算：設調查海域布設m（j=1，2…m）個站位（潮間帶為布設m個樣點），每個站位（樣點）共有n（i=1，2…n）個物種。則調查海域大型底棲動物的次級生產(chǎn)力P為：

其中，Pj為第j個站位的次級生產(chǎn)力，Pi為第j個站位第i種的次級生產(chǎn)力。Ai和Bi代表第j個站位第i種的棲息密度和去灰干重。如果Ai和Bi分別代表第j個站位第i種的一個季節(jié)的棲息密度和去灰干重，則P為該季節(jié)的次級生產(chǎn)力；如果Ai和Bi分別代表第j個站位第i種的年平均密度和去灰干重，則P為年平均次級生產(chǎn)力。

按站位（潮間帶為斷面或樣點）計算：計算每個站位（斷面或樣點）大型底棲動物的棲息密度和生物量，然后代入Brey（1990）經(jīng)驗公式得到每個站位（斷面或樣點）大型底棲動物的生產(chǎn)力。如果每個站位（斷面或樣點）大型底棲動物的棲息密度和生物量是一個季節(jié)（或一個航次）的數(shù)值，則結果為該季節(jié)（或該航次）一個站位（斷面或樣點）大型底棲動物的次級生產(chǎn)力；如果每個站位（斷面或樣點）大型底棲動物的棲息密度和生物量是一年（或幾個航次）的平均值，則結果為該站位（斷面或樣點）大型底棲動物的年平均次級生產(chǎn)力。

按類群計算：將大型底棲動物分為多毛類、軟體動物、甲殼動物、棘皮動物和其他類動物共5個類群，求出每個站位（斷面或樣點）各類群的棲息密度和生物量，分別計算其次級生產(chǎn)力。各類群的生產(chǎn)力相加得到該站位（斷面或樣點）的生產(chǎn)力。計算公式如下：

不同類群模型系數(shù)a，b1和b2見表3。其中：P為大型底棲動物年次級生產(chǎn)力，單位為g（AFDW）/（m2·a）；B為大型底棲動物年平均去灰干重，單位為g（AFDW）/m2；W為大型底棲動物個體的年平均去灰干重，單位為g（AFDW）/ind。

表3 Brey（1990）經(jīng)驗公式中不同類群的3個系數(shù)值

該公式可轉換為：

見表4。

表4 Brey（1990）轉換公式中不同類群的3個系數(shù)值

Brey（1990）給出采用經(jīng)驗公式估算群落次級生產(chǎn)力的參考步驟為：（1）計算年所有物種的平均生物量B和年平均個體體重W，（2）通過經(jīng)驗公式計算每一個種的估算值，盡可能使用按類群的經(jīng)驗公式，（3）把所有值加起來得到群落的產(chǎn)量?？梢钥闯?，逐種計算和按類群計算符合該步驟，按站位計算是對該步驟的一種簡化。

在國內開始嘗試利用Brey（1990）經(jīng)驗公式的時候，國外誕生了Brey（2001）估算模型，見公式（8）和表2。張崇良（2011）首次引入Brey（2001）估算模型對膠州灣潮間帶大型底棲動物群落次級生產(chǎn)力進行研究，紀瑩露（2015）估算遼東灣西部海域時也采用該模型。具體算法如下：將所有航次的大型底棲動物棲息密度和生物量取平均值作為年平均棲息密度和年平均生物量，然后利用該模型逐種計算年平均次級生產(chǎn)力，每站所有物種平均次級生產(chǎn)力相加得到該站年平均次級生產(chǎn)力。從網(wǎng)站（http：//www.thomas-brey.de/ science/virtualhandbook/）下載操作手冊，按說明進行數(shù)據(jù)處理。此外，Brey（2012）在Brey（2001）估算模型的基礎上增加大型底棲動物食性和移動性參數(shù)提出改進模型，由于模型比較復雜，可登陸上述網(wǎng)站下載計算表單。

在次級生產(chǎn)力估算模型研究發(fā)展的第3階段，國外開始朝向特定生境中大型底棲動物次級生產(chǎn)力估算模型研究發(fā)展。Cusson等（2005）分析了全球海洋大型底棲動物次級生產(chǎn)力的分布模式，提出了特定生境中次級生產(chǎn)力估算模型，將估算模型中水深的適用范圍由0~300 m提高到0~930 m，見表5。Petracco（2012）采用9種估算模型研究了沙灘生境大型底棲動物次級生產(chǎn)力，以評價不同模型對該生境次級生產(chǎn)力估算的適宜性，結果表明Cusson等（2005）沙質估算模型最優(yōu)。但該模型中生物壽命是解釋P/B值很重要的變量，而生物壽命不是經(jīng)驗模型的優(yōu)選參數(shù)（Cusson et al，2005），因此該模型的實用性還有待論證。

表5 Cusson等（2005）估算模型參數(shù)

2 存在的問題

2.1 模型使用限制條件及自身缺陷

經(jīng)驗模型可有效估算大型底棲動物次級生產(chǎn)力（Dolbeth et al，2005），但使用時應注意其應用范圍和限制條件，Brey（1990）經(jīng)驗公式建模數(shù)據(jù)大部分來自北半球溫帶淺海和潮間帶，采用年平均生物量和年平均個體體重為參數(shù)，適合估算大型底棲動物年平均次級生產(chǎn)力（張崇良等，2011）。Tumbiolo等（1994）估算模型建模數(shù)據(jù)主要集中于北緯50°-60°的歐洲海域，不一定適合低緯度海域。高生物量單一物種主導的群落、物種數(shù)小于5的群落及暴露于極端環(huán)境的群落，其次級生產(chǎn)力估算均不適宜采用經(jīng)驗模型，否則偏差較大（Brey et al，1990；Dolbeth et al，2005）。

現(xiàn)存經(jīng)驗模型考慮的環(huán)境因子僅有水深、水溫，但影響次級生產(chǎn)力的環(huán)境因子有很多，如溶解氧、底質、葉綠素a、水體總氮、總磷等（Brey et al，1990；龔志軍等，2001；嚴娟等，2012），而且特定海域特定季節(jié)次級生產(chǎn)力的影響因素也不同。從經(jīng)驗模型的研究進展可以看出，建模者試圖建立適用范圍更廣、包含物種更多的模型，但由于特定海域次級生產(chǎn)力的影響因素不同，建立普適性的經(jīng)驗模型較為困難。因此，可以參照Cusson等（2005）的研究，建立特定海域或特定生境次級生產(chǎn)力估算模型。而我國海洋大型底棲動物次級生產(chǎn)力對單個種群基礎研究較少，建立適合我國特定海域的次級生產(chǎn)力經(jīng)驗模型尚有一定困難。

2.2 不同方法估算同一海域次級生產(chǎn)力結果不同

國內大型底棲動物群落次級生產(chǎn)力估算基本采用Brey（1990）經(jīng)驗公式，但在應用時存在3種類型，而Brey（2001）估算模型只有兩篇報道（張崇良等，2011；紀瑩露等，2015）采用。為比較不同估算方法產(chǎn)生的結果差異，下面以2006年9月和2007年5月杭州灣潮間帶大型底棲動物調查資料為基礎，分別采用經(jīng)驗公式逐種、按站位和按類群以及Brey（2001）估算模型計算該海域大型底棲動物群落次級生產(chǎn)力。2006年9月和2007年5月杭州灣潮間帶大型底棲動物調查資料由浙江省“908專項”調查課題組提供，圖1為采樣斷面分布圖。

圖1 潮間帶大型底棲動物采樣斷面分布圖

將春秋兩航次大型底棲動物豐度和生物量取平均值作為年平均豐度和年平均生物量，然后采用采用經(jīng)驗公式逐種、按站位與按類群分別計算各斷面年平均次級生產(chǎn)力，濕重的9/50為去灰干重（李新正等，2005），所有斷面次級生產(chǎn)力的平均值為該海域大型底棲動物群落的年平均次級生產(chǎn)力。同時采用Brey（2001）估算模型計算大型底棲動物群落次級生產(chǎn)力，并按照表1轉換系數(shù)，將估算模型的結果轉換為去灰干重（紀瑩璐等，2015）。采用SPSS 16.0軟件對大型底棲動物次級生產(chǎn)力進行以斷面和不同估算模型為因子的雙因素方差分析，顯著性水平設置為α=0.05。

不同方法的估算結果見表6。杭州灣潮間帶大型底棲動物年平均次級生產(chǎn)力在1.09~3.35 g（AFDW）/（m2·a）之間，在全國近岸海域處于較低水平。經(jīng)驗公式的估算結果低于Brey（2001）估算模型的結果，經(jīng)驗公式3種計算方法之間相比較，按類群計算結果最小，按站位計算結果最大，逐種計算的結果在二者之間。雙因素方差分析表明（表7），不同斷面對大型底棲動物次級生產(chǎn)力有極顯著影響（P＜0.01），不同估算模型對大型底棲動物次級生產(chǎn)力有顯著影響（P＜0.05）。不同估算方法間大型底棲動物次級生產(chǎn)力均數(shù)多重比較結果表明（表8），Brey（1990）經(jīng)驗公式3種計算方法間差異不顯著，Brey（1990）經(jīng)驗公式和Brey（2001）估算模型間的次級生產(chǎn)力存在顯著差異（P＜0.05）。

模型中自變量解釋因變量變異均達到75%以上（Brey et al，1990；Brey et al，2001），表明以上估算模型本身均是合理的。Brey（1990）經(jīng)驗公式的3種算法中，按類群計算最符合Brey（1990）給出的參考步驟，其次是逐種計算，按站位計算是對經(jīng)驗公式的簡化，3種計算方法間沒有顯著性差異。Dolbeth（2005）和Cusson（2005）的研究表明，Brey（2001）估算模型結果偏差最小。因此，建議采用Brey（2001）估算模型計算河口海域大型底棲動物群落次級生產(chǎn)力。另外，由于估算值與真實值的偏差隨機分布，不同估算值合并偏差則可相互抵消，因此采用不同模型估算同一海域大型底棲生物群落次級生產(chǎn)力，則估算結果確定的范圍和平均值應是比較可靠的。

表6 杭州灣潮間帶各斷面大型底棲動物次級生產(chǎn)力

表7 不同斷面、估算方法對次級生產(chǎn)力的方差分析結果

表8 不同估算方法間次級生產(chǎn)力均數(shù)的兩兩比較

2.3 同一估算方法估算不同生境不同群落次級生產(chǎn)力結果偏差不同

不同生境大型底棲動物群落因種類組成和數(shù)量組成不同造成次級生產(chǎn)力存在明顯差異，如次級生產(chǎn)力海草床海域＞軟泥底質海域＞潮間帶海域（Tumbiolo et al，1994）、潮間帶大葉藻海域＞泥沙底質海域（Dolbeth et al，2005）、硬質底質＞泥沙底質（Cusson et al，2005），而同一估算模型計算不同生境不同群落底棲動物次級生產(chǎn)力產(chǎn)生的結果偏差往往也不同。

一般將經(jīng)典算法結果作為次級生產(chǎn)力精確計算值，以此評價不同估算模型的結果偏差（Brey et al，2001；Dolbeth et al，2005；Cusson et al，2005）。Brey（1990）采用經(jīng)驗公式計算了4個不同群落底棲動物次級生產(chǎn)力，估算值偏差分別為-11%，-5%，+1%，+66%；Dolbeth（2005）采用6種次級生產(chǎn)力估算模型分別計算葡萄牙蒙德古河口潮間帶大葉藻床海域和泥沙海域兩種生境5個河口種的生產(chǎn)力，發(fā)現(xiàn)Brey（2001）估算模型估算偏差最穩(wěn)定且最小（-4%~+10%），Robertson（1979）模型估算偏差在10%~49%，Schwinghamer（1986）估算模型估算大葉藻床海域底棲動物次級生產(chǎn)力偏差較小，而Tumbiolo等（1994）估算模型估算泥沙海域底棲動物次級生產(chǎn)力偏差較小。估算模型的優(yōu)劣與底棲動物種類組成和模型采用的參數(shù)密切相關（Brey et al，1990；Dolbeth et al，2005），因此特定海域大型底棲動物次級生產(chǎn)力估算模型的選擇需考慮其底棲動物組成特點并結合環(huán)境因子分析，而不同估算模型對特定海域次級生產(chǎn)力估算的適用性有待進一步討論。

3 總結與展望

估算模型是對海洋大型底棲動物次級生產(chǎn)力這一生態(tài)過程的模擬和預測，極大地方便了次級生產(chǎn)力的估算。但模型只是提供了一個估算值，沒有采用經(jīng)典算法計算那樣精確。海洋大型底棲動物群落次級生產(chǎn)力估算模型研究可分為3個階段：起步階段（1979-1990）、發(fā)展與完善階段（1990-2001）、成熟并廣泛應用階段（2001-至今）。

經(jīng)驗模型可有效估算大型底棲動物次級生產(chǎn)力，但使用時應注意其應用范圍和限制條件。影響次級生產(chǎn)力的環(huán)境因子有很多，但現(xiàn)存經(jīng)驗模型僅考慮水深、水溫兩個環(huán)境因子。不同方法估算同一海域次級生產(chǎn)力結果不同，同一估算方法估算不同生境不同群落次級生產(chǎn)力結果偏差也不同，應建立特定海域或特定生境的次級生產(chǎn)力估算模型。

我國海洋大型底棲動物群落次級生產(chǎn)力估算模型研究缺少對單個種群進行大量基礎研究的階段，而是直接跳躍到采用經(jīng)驗模型估算的階段。采用經(jīng)典算法如同生群法和體長頻率法估算單個種群次級生產(chǎn)力的研究很少（吳寶鈴等，1992；周一兵等，1995；林巋璇等，2008），沒有足夠基礎數(shù)據(jù)建模。從國外引進的群落次級生產(chǎn)力估算模型是通過單個種群生產(chǎn)力資料建模得到，由于中國海洋環(huán)境特征、開發(fā)與管理方式、大型底棲動物種類組成與國外存在一定差異，不能單純照搬國外模型，如Brey（1990）模型建模數(shù)據(jù)主要集中在北緯30°-70°、水深50 m以內淺海，對渤海、黃海、東?？赡苓m用，但對緯度較低、水深較大的南海不一定適用；Tumbiolo等（1994）模型建模數(shù)據(jù)主要集中于北緯50°-60°的西歐海域，不適于我國海域次級生產(chǎn)力估算。因此，一方面需加強我國海域大型底棲動物單個種群次級生產(chǎn)力的研究，另一方面應借鑒國際主流次級生產(chǎn)力估算模型的相關經(jīng)驗，結合中國海生態(tài)系統(tǒng)的特點，根據(jù)所獲基礎資料建立適于我國特定海域或特定生境的次級生產(chǎn)力估算模型，在時空尺度上，分析不同區(qū)域、不同季節(jié)大型底棲動物的次級生產(chǎn)力狀況。次級生產(chǎn)力與生物棲息的環(huán)境密切相關，對特定海域大型底棲動物群落次級生產(chǎn)力估算較準確的做法是先獲得單個種群的生產(chǎn)力數(shù)據(jù)，然后建立的模型可以對該海域次級生產(chǎn)力這一生態(tài)過程做出較好的預測。

B enee A,1993.Concepts and patterns of invertebrate production in running waters.Proceedings-international Association of Theoretical and Applied Limnology,25:15-25.

Brey T,1990.Estimating productivity of macrobenthic invertebrates from biomass and mean individual weight.Archive of Fishery and Marine Research,32(4):329-343.

Brey T,2001.Population dynamics in benthic invertebrates.http:// thomas-brey.de/science/virtualhandbook.

Brey T,2012.A multi-parameter artificial neural network model to estimate macrobenthic invertebrate productivity and production. Limnol.Oceanogr.Methods 10,581-589.

Brey T,Gerdes D,1998.High Antarctic macrobenthic community production.Journal of Experimental Marine Biology and Ecology, 231(2):191-200.

Cusson M,Bourget E,2005.Global patterns of macroinvertebrate production in marine benthic habitats.Marine Ecology Progress Series,297:1-14.

Dolbeth M,Lilleb? A,Cardoso P,et al,2005.Annual production of estuarine fauna in different environmental conditions:an evaluation of the estimation methods.Journal of Experimental Marine Biology and Ecology,326(2):115-127.

Edgar G J,1990.The use of the size structure of benthic macrofaunal communities to estimate faunal biomass and secondary production. Journal of Experimental Marine Biology and Ecology,137(3):195-214.

Mann K H,1969.The dynamics of aquatic ecosystems.http://www.books .goole.com.

Morin A,Bourassa N,1992.Modèles empiriques de la productionannuelle et du rapport P/B d'invertébrés benthiques d'eau courante. Canadian Journal of Fisheries and Aquatic Sciences,49(3):532-539.

Petracco M,Cardoso R S,Corbisier T N,et al,2012.Secondary production of sandy beach macrofauna:An evaluation of predictive models.Estuarine,Coastal and Shelf Science,115(10):359-365.

Plante C,Downing J A,1989.Production of freshwater invertebrate populations in lakes.Canadian Journal of Fisheries and Aquatic Sciences,46(9):1489-1498.

Robertson A I,1979.The relationship between annual production: biomass ratios and lifespans for marine macrobenthos.Oecologia,38 (2):193-202.

Schwinghamer P,Hargrave B,Peer D,et al,1986.Partitioning of production and respiration among size groups of organisms in an intertidal benthic community.Marine Ecology Progress Series,31 (2):131-142.

Sprung M,1993.Estimating macrobenthic secondary productor from body weight and biomass:a field test in a non-boreal intertidal habitat. Marine Ecology-Progress Series,100:103-103.

Tumbiolo M L,Downing J A,1994.An empirical model for the prediction of secondary production in marine benthic invertebrate populations. Marine Ecology Progress Series,114:165-174.

Warwick R,1980.Population dynamics and secondary production of benthos.Marine benthic dynamics,11.

Waters T F,1977.Secondary production in inland waters.Advances in ecological research,10:91-164.

蔡立哲，許鵬，傅素晶，等，2012.湛江高橋紅樹林和鹽沼濕地的大型底棲動物次級生產(chǎn)力.應用生態(tài)學報，23（4）：965-971.

杜飛雁，王雪輝，李純厚，等，2008.大亞灣大型底棲動物生產(chǎn)力變化特征.應用生態(tài)學報，19（4）：873-880.

杜永芬，高抒，于子山，等，2012.福建羅源灣潮間帶大型底棲動物的次級生產(chǎn)力.應用生態(tài)學報，23（7）：1904-1912.

方平福，章書聲，鮑毅新，等，2013.溫州天河潮間帶大型底棲動物的次級生產(chǎn)力.生態(tài)學雜志，32（1）：106-113.

方少華，呂小梅，張躍平，等，2009.湄洲灣東吳港區(qū)附近潮間帶大型底棲動物的時空分布及次級生產(chǎn)力.臺灣海峽，28（3）：392-398.

龔志軍，謝平，閻云君，2001.底棲動物次級生產(chǎn)力研究的理論與方法.湖泊科學，13（1）：79-88.

韓慶喜，李寶泉，韓秋影，等，2011.漁業(yè)捕撈對威海港附近海域底上大型底棲群落結構影響的初步研究.海洋通報，30（2）：121-126.

紀瑩璐，趙寧，楊傳平，等，2015.遼東灣西部海域大型底棲動物群落次級生產(chǎn)力初探.中國海洋大學學報：自然科學版，45（4）：53-58.|

焦海峰，施慧雄，尤仲杰，等，2011.浙江漁山列島巖礁潮間帶大型底棲動物次級生產(chǎn)力.應用生態(tài)學報，22（8）：2173-2178.

金亮，蔡立哲，周細平，等，2007.深圳灣北岸泥灘大型底棲動物次級生產(chǎn)力研究.臺灣海峽，26（3）：415-421.

金亮，林秀春，蔡立哲，2009.湄洲灣貝類養(yǎng)殖灘涂大型底棲動物次級生產(chǎn)力研究.海洋環(huán)境科學，28（1）：8-11.

李新正，王洪法，張寶琳，2005.膠州灣大型底棲動物次級生產(chǎn)力初探.海洋與湖沼，36（6）：49-55.

林巋璇，韓潔，林旭吟，等，2008.廈門潮間帶小頭蟲（Capitella capitata）的種群動態(tài)及次級生產(chǎn)力研究.北京師范大學學報：自然科學版，44（3）：314-317.

劉坤，林和山，王建軍，等，2015.廈門近岸海域大型底棲動物次級生產(chǎn)力.生態(tài)學雜志，34（12）：3409-3415.

劉曉收，趙瑞，華爾，等，2014.萊州灣夏季大型底棲動物群落結構特征及其與歷史資料的比較.海洋通報，33（3）：283-292.

劉勇，線薇薇，孫世春，等，2008.長江口及其鄰近海域大型底棲動物生物量、豐度和次級生產(chǎn)力的初步研究.中國海洋大學學報：自然科學版，38（5）：749-756.

呂小梅，方少華，張躍平，等，2008.福建海壇海峽潮間帶大型底棲動物群落結構及次級生產(chǎn)力.動物學報，54（3）：428-435.

田勝艷，張文亮，于子山，等，2010.膠州灣大型底棲動物的豐度，生物量和生產(chǎn)量研究.海洋科學，34（6）：81-87.

王金輝，楊春旺，孫亞偉，等，2006.象山港大型底棲動物的生物多樣性和次級生產(chǎn)力研究.天津農(nóng)學院學報，13（2）：24-28.

王志忠，張金路，陳述江，等，2012.黃河入?？诔遍g帶大型底棲動物群落組成及次級生產(chǎn)力.海洋環(huán)境科學，31（5）：657-661.

王宗興，孫丕喜，劉彩霞，等，2011.桑溝灣大型底棲動物的次級生產(chǎn)力.應用與環(huán)境生物學報，17（4）：495-498.

吳寶鈴，丘建文，1992.青島多齒圍沙蠶的生產(chǎn)量.生態(tài)學報，12（1）：61-67.

吳辰，蔡立哲，陳昕韡，等，2013.福建漳江口紅樹林和鹽沼濕地大型底棲動物次級生產(chǎn)力研究.廈門大學學報：自然科學版，52（2）：259-266.

嚴娟，莊平，侯俊利，等，2012.長江口潮間帶大型底棲動物次級生產(chǎn)力及其影響因子.應用與環(huán)境生物學報，18（6）：935-942.

于子山，張志南，韓潔，2001.渤海大型底棲動物次級生產(chǎn)力的初步研究.青島海洋大學學報：自然科學版，31（6）：867-871.

袁偉，張志南，于子山，2007.膠州灣西部海域大型底棲動物次級生產(chǎn)力初步研究.應用生態(tài)學報，18（1）：145-150.

張崇良，徐賓鐸，任一平，等，2011.膠州灣潮間帶大型底棲動物次級生產(chǎn)力的時空變化.生態(tài)學報，31（17）：5071-5080.

張海波，蔡燕紅，王薇，2010.象山港潮間帶大型底棲動物次級生產(chǎn)力初探.寧波大學學報：理工版，23（4）：26-30.

周福芳，史秀華，邱國玉，等，2012.深圳灣不同生境濕地大型底棲動物次級生產(chǎn)力的比較研究.生態(tài)學報，32（20）：6511-6519.

周進，紀煒煒，2012.三都澳大型底棲動物次級生產(chǎn)力.海洋漁業(yè)，34（1）：32-38.

周進，李新正，李寶泉，2008.黃海中華哲水蚤度夏區(qū)大型底棲動物的次級生產(chǎn)力.動物學報，54（3）：436-441.

周細平，蔡立哲，梁俊彥，等，2008.廈門海域大型底棲動物次級生產(chǎn)力的初步研究.廈門大學學報：自然科學版，47（6）：902-906.

周一兵，謝祚渾，1995.蝦池中日本刺沙蠶的次級生產(chǎn)力研究.水產(chǎn)學報，19（2）：140-150.

（本文編輯：袁澤軼）

Research progress and application evaluation of predictive models for secondary production of marine macrobenthos community

JIA Sheng-hua,SHOU Lu,LIAO Yi-bo,ZENG Jiang-ning,GAO Ai-gen,TANG Yan-bin
（Keylaboratory of Marine EcosystemandBiogeochemistry,Second Instituteof Oceanography, State Oceanic Administration,Hangzhou 310012,China）

The benthic fauna is an important component of the marine ecosystem,and plays a vital role in material cycle and energy flow of marine ecosystem.It is of great significance in the marine ecosystem to accurately estimate the secondary production of macrobenthos community.In this study,a broad literature review of predictive models for secondary production of marine macrobenthos community was conducted and the three study stages were outlined including the initial stage（from 1979 to 1990）,development and perfection stage（from 1990 to 2001）,mature and wide-spreading stage（2001 to the present）.In general,Brey's（1990）empirical formula has been widely used in the estimation of secondary production of marine macrobenthos community in China,and the way of using it can be divided into three levels,computing secondary production based on species,stations and groups.Of those,computation based on species or groups is most appropriate for Brey's（1990）reference steps.Secondary production estimated by different empirical models in the same area is different, and the deviations produced by the same empirical model in the different habitats or communities are also different.The future research work of empirical model to estimate the secondary production of marine macrobenthos community should focus on the research of the secondary production of populations,so we can establish a proper empirical model for a certain sea area in China on the basis of accumulating a lot of basic data.

marine macrobenthos;secondary production of community;empirical models;intertidal zone;Hangzhou Bay

Q179.4

A

1001-6932（2017）04-0370-09

10.11840/j.issn.1001-6392.2017.04.002

2016-03-24；

2016-06-13

海洋公益性行業(yè)科研專項（201305009；201305043；201405007；201505004-3）；國家海洋局第二海洋研究所基本科研業(yè)務費專項（JG1311；JG1312)。

賈勝華（1990-），碩士研究生。主要從事海洋生態(tài)學/大型底棲生物研究。電子郵箱：jiashenghuahh@163.com。

曾江寧，博士，研究員。電子郵箱：jiangningz@126.com。