王迎賓
(1.浙江海洋大學(xué) 水產(chǎn)學(xué)院,浙江 舟山 316022)
在不斷增加的捕撈強(qiáng)度的影響下,漁業(yè)資源遭受著過(guò)度捕撈的壓力,資源的可持續(xù)利用面臨挑戰(zhàn)[1]。隨著漁業(yè)生物育苗和養(yǎng)殖技術(shù)的快速發(fā)展,人類已經(jīng)可以通過(guò)采取增殖放流等人為介導(dǎo)手段,大規(guī)模地對(duì)自然生態(tài)開展人工干預(yù)。從最初單純以增加商業(yè)捕撈對(duì)象資源量為目的增殖放流,到如今以構(gòu)建健康的半人工生態(tài)系統(tǒng)為目的海洋牧場(chǎng)建設(shè),使?jié)O業(yè)資源得以延續(xù)。國(guó)際《海洋科學(xué)百科全書》對(duì)“海洋牧場(chǎng)”有一個(gè)簡(jiǎn)單而明確的定義[2?3],即海洋牧場(chǎng)通常是指資源增殖。由此可見,海洋牧場(chǎng)的建設(shè)與管理離不開資源增殖。但與最初的僅以增加產(chǎn)量為目的的增殖放流有所不同,如今人類對(duì)人為介導(dǎo)所產(chǎn)生影響的關(guān)注也逐漸延伸到資源的可持續(xù)性,例如,海洋牧場(chǎng)海域內(nèi)增殖目標(biāo)種的最可持續(xù)產(chǎn)量(Maximum Sustainable Yield,MSY)水平,就是眾多學(xué)者和牧場(chǎng)經(jīng)營(yíng)者十分關(guān)注的問(wèn)題之一。
MSY 被定義為在一定的海域范圍內(nèi)和現(xiàn)有環(huán)境條件下,可以連續(xù)獲得的最大平均產(chǎn)量,長(zhǎng)期以來(lái)一直被全球許多漁業(yè)視為重要的管理目標(biāo)。20 世紀(jì)50 年代,隨著剩余產(chǎn)量模型的發(fā)展,MSY 被廣泛使用,其概念在1982 年《聯(lián)合國(guó)海洋法公約》中被采納,而且作為法律條款寫入諸多漁業(yè)發(fā)達(dá)國(guó)家的漁業(yè)法規(guī)當(dāng)中,許多國(guó)際組織也將MSY 作為主要的管理目標(biāo)之一。剩余產(chǎn)量模型自提出之時(shí),就受到了廣泛歡迎,很多學(xué)者也對(duì)傳統(tǒng)模型進(jìn)行了完善和改進(jìn)。Thiaw 等[4]將上升流等環(huán)境因子加入剩余產(chǎn)量模型,對(duì)塞內(nèi)加爾白蝦(Penaeus notialis)資源進(jìn)行了評(píng)估。Reid 等[5]使用基于狀態(tài)空間的剩余產(chǎn)量模型,估算了尼皮貢湖白鮭(Coregonus clupeaformis)的主要生物學(xué)參考點(diǎn)。Shibata 等[6]提出了對(duì)傳統(tǒng)剩余產(chǎn)量模型進(jìn)行了調(diào)整,并基于此分析了季節(jié)間在兩個(gè)區(qū)域來(lái)回游動(dòng)的漁業(yè)資源的產(chǎn)量變化。然而從生物群體角度看,傳統(tǒng)的以及上述改進(jìn)的剩余產(chǎn)量模型均只考慮了生物群體自身的繁衍,而無(wú)法解決像海洋牧場(chǎng)中受增殖放流影響較大的生物群體MSY 的估算問(wèn)題。針對(duì)這一問(wèn)題,國(guó)內(nèi)外至今尚未見到有相關(guān)研究報(bào)道。
海洋牧場(chǎng)海域界限比較清晰,并且其中的定棲性種類移動(dòng)性較弱,使得資源的環(huán)境容納量相對(duì)穩(wěn)定。相反,開放海域生產(chǎn)作業(yè)的邊界不夠明確,而且在實(shí)施增殖放流后,放流種類的棲息海域可能有所擴(kuò)展,特別是游動(dòng)性較強(qiáng)的種類,從而使研究海域的環(huán)境容納量難以確定。因此,本研究針對(duì)海洋牧場(chǎng)海域定棲性種類,建立適用于增殖放流影響下的剩余產(chǎn)量模型,并通過(guò)模擬,估算不同增殖及捕撈策略下MSY 的評(píng)估結(jié)果,以期為增殖放流背景下的資源可持續(xù)利用提供理論依據(jù)。
剩余產(chǎn)量模型是現(xiàn)代漁業(yè)資源評(píng)估和管理的主要理論模型之一,其將每年種群的生長(zhǎng)、補(bǔ)充和自然死亡結(jié)合起來(lái)作為種群的年剩余產(chǎn)量,表示除了捕撈以外的種群年綜合生產(chǎn)力,因此捕撈和剩余產(chǎn)量是影響種群生物量的兩個(gè)因素[7]。Schafer 提出的剩余產(chǎn)量模型,是最早被廣泛采用的產(chǎn)量模型,該模型假定生物量增長(zhǎng)率與生物量的關(guān)系服從邏輯斯蒂(Logistic)分布[8?9]。它的具體公式為
式中,Bt表示第t年的資源生物量;r為內(nèi)稟增長(zhǎng)率;K為環(huán)境容納量;Ct表示第t年的產(chǎn)量。
在非平衡狀態(tài)下,離散型剩余產(chǎn)量模型的一般形式為
當(dāng)不考慮捕撈時(shí),上述模型通過(guò)r來(lái)實(shí)現(xiàn)群體自身的邏輯斯蒂增長(zhǎng),即生物量的年際增長(zhǎng)服從S 型曲線(圖1)。當(dāng)考慮捕撈時(shí),人們希望找到能保證生物量可持續(xù)時(shí)的最大產(chǎn)量,即MSY。上述模型中,r表征的是資源棲息海域生物量的自身增長(zhǎng),也就是生物通過(guò)自身繁殖、生長(zhǎng)等作用,呈現(xiàn)出的生物量增長(zhǎng)率。因此,當(dāng)有外部因素介入,例如增殖放流,上述模型就有局限性了。因?yàn)榉帕鞯膫€(gè)體雖然也增加了群體的生物量,但其增長(zhǎng)特征無(wú)法用表示群體增長(zhǎng)率的r進(jìn)行描述。因此,須在公式(2)中引入其他參數(shù),表示增殖放流對(duì)群體生物量變化帶來(lái)的影響。
基于此考慮,在公式(2)的基礎(chǔ)上加入描述增殖群體增長(zhǎng)特征的參數(shù)—有效增殖率(e),用以表示增殖放流的群體生長(zhǎng)至可捕撈規(guī)格的生物量比例。于是,得到增殖剩余產(chǎn)量模型公式
式中,Bt表示第t年原存群體(海域原本存在的群體)生物量;St表示第t年增殖放流苗種全部長(zhǎng)到可捕成體規(guī)格時(shí)的生物量。該模型將海域中原存群體與增殖群體對(duì)生物量的影響分開考慮。因?yàn)樵鲋橙后w并未改變?cè)S驅(qū)ι锕逃械娜菁{能力,所以模型中的K值大小并未改變。因此,該模型的曲線將同樣呈現(xiàn)S 型,并趨于同樣的極值。
圖1 海洋牧場(chǎng)刺參(Stichopus japonicus)增殖模擬生物量和產(chǎn)量變化Fig.1 Simulation changes of biomass and yield of Stichopus japonicus enhancement in marine ranching
為了分析增殖剩余產(chǎn)量模型在不同增殖放流策略下的使用效果,本研究基于增殖定棲性種類刺參(Stichopus japonicus)的資料進(jìn)行模擬研究,主要考慮定棲性種類活動(dòng)范圍小,生活邊界比較清晰。明確的生活范圍邊界,是剩余產(chǎn)量模型應(yīng)用的前提,估算得到的MSY 與該范圍也是相對(duì)應(yīng)的。因此,以定棲性種類為例,所分析的問(wèn)題比較簡(jiǎn)單,誤差也相對(duì)較小。
增殖刺參群體20 年的生物量和產(chǎn)量的變動(dòng)過(guò)程如圖1 所示?;谠鲋呈S喈a(chǎn)量模型,模擬比較無(wú)增殖放流和每年增殖放流量穩(wěn)定的情況下,產(chǎn)量與生物量,以及產(chǎn)量與捕撈努力量之間的關(guān)系。對(duì)于每種增殖策略,分析其MSY、捕撈努力量(EMSY)和生物量(BMSY)的變動(dòng)特征。同時(shí),對(duì)于不同水平的年增殖放流量(0.36 t、0.48 t 和0.60 t),以及年增殖放流量持續(xù)變動(dòng)的情況,也進(jìn)行了模擬分析,用于對(duì)比。模擬使用相關(guān)參數(shù)如表1所示。
為了探討不同捕撈強(qiáng)度對(duì)增殖剩余產(chǎn)量模型評(píng)估結(jié)果產(chǎn)生的影響,在不同增殖放流數(shù)量基礎(chǔ)上(年增殖放流量穩(wěn)定、持續(xù)增加、持續(xù)減少和隨機(jī)4 種增殖策略,增殖放流均從第8 年開始,相關(guān)參數(shù)如表1所示),基于現(xiàn)實(shí)捕撈努力量和產(chǎn)量數(shù)據(jù),模擬形成了產(chǎn)量的3 種變化趨勢(shì),即持續(xù)增加、持續(xù)下降和先增后降(圖2)??紤]到產(chǎn)量和捕撈努力量的觀測(cè)值可能存在誤差,在模擬過(guò)程中對(duì)產(chǎn)量和捕撈努力量均加入隨機(jī)誤差,ε~N(0,σ2),其中σ等于產(chǎn)量或捕撈努力量均值的10%。模擬中參數(shù)的估算使用觀測(cè)誤差法,每種增殖策略和捕撈策略情況均模擬重復(fù)1 000 次,然后繪制箱型圖來(lái)呈現(xiàn)模擬結(jié)果的分布情況。
與傳統(tǒng)的Schaefer 剩余產(chǎn)量模型相比,增殖剩余產(chǎn)量模型增加了1 個(gè)參數(shù)和1 個(gè)變量,即有效增殖率e和年增殖放流苗種全部長(zhǎng)到可捕成體規(guī)格時(shí)的生物量S。增殖剩余產(chǎn)量模型是在Schaefer 模型基礎(chǔ)上提出的,因此其適用條件和假設(shè)與傳統(tǒng)模型既有相似之處,也有所不同。首先,在有限生態(tài)系統(tǒng)中,資源群體的資源生物量可增長(zhǎng)到接近該生態(tài)系統(tǒng)的最大容納量,并且該容納量不因增殖放流而發(fā)生變化;第二,環(huán)境容納量接近于初始資源群體;第三,生物量的增長(zhǎng)服從邏輯斯蒂的S 型曲線模式;第四,增殖放流群體與海域中原存群體之間會(huì)產(chǎn)生抑制作用,每年增殖放流量越大,原存群體的最大容納量越小,同時(shí)原存群體的MSY 以及相應(yīng)的EMSY越小。根據(jù)增殖剩余產(chǎn)量模型計(jì)算得到的生物量年際變化曲線雖然呈現(xiàn)出S 型,但與公式(2)的曲線相比,其曲線逼近極值的速度加快了,生物量增長(zhǎng)率(dBt/dt)的拐點(diǎn)也提前出現(xiàn)(圖3)。由于增殖剩余產(chǎn)量模型的結(jié)構(gòu)與傳統(tǒng)的生物量邏輯斯蒂生長(zhǎng)模型相比有所變化,其相應(yīng)的一些關(guān)鍵參考點(diǎn)計(jì)算公式也有所改變。MSY 表示在保證資源可持續(xù)利用的前提下,每年可獲得的最大產(chǎn)量。因此,增殖剩余產(chǎn)量模型需滿足
表1 以刺參(Stichopus japonicus)增殖為基礎(chǔ)的海洋牧場(chǎng)定棲性種類模擬參數(shù)Table 1 Simulation parameters of Stichopus japonicus enhancement in marine ranching
在每年增殖放流量穩(wěn)定的前提下,將產(chǎn)量對(duì)生物量求導(dǎo),并令其等于0,即
圖2 模擬4 種捕撈努力量和產(chǎn)量的變化趨勢(shì)Fig.2 Trends of fishing effort and yield of the four simulated scenarios
圖3 有無(wú)增殖放流情況下漁業(yè)資源群體邏輯斯蒂增長(zhǎng)示意圖Fig.3 Schematic diagram of logistic growth of fishery stock with or without enhancement
則
此時(shí),對(duì)應(yīng)的MSY 的計(jì)算公式為
可見,與傳統(tǒng)的Schafer 剩余產(chǎn)量模型相比,使用增殖剩余產(chǎn)量模型求算的MSY,除了與參數(shù)r和K有關(guān),還受到增殖放流的影響。換言之,可將增殖放流背景下的MSY 視為來(lái)自兩個(gè)方面,即海域原存群體和增殖放流群體。
對(duì)公式(4)進(jìn)行變換得到
式中,It表示第t年的CPUE;q表示可捕系數(shù)。進(jìn)一步變換得到
在此設(shè)定公式(10)等于0,即意味著dBt/dt=0。當(dāng)每年增殖數(shù)量穩(wěn)定,并且給定每年的捕撈努力量大小時(shí),公式(10)便是1 個(gè)一元二次方程,可以計(jì)算得到dBt/dt=0時(shí)的產(chǎn)量,即給定年增殖放流數(shù)量和捕撈努力量時(shí)的MSY 值。這樣,可以在固定每年增殖放流數(shù)量情況下,通過(guò)設(shè)定不同的捕撈努力量,得到年增殖放流量一定時(shí)的MSY 及對(duì)應(yīng)的EMSY。相反,同樣可以固定每年的捕撈努力量,而設(shè)定不同的年放流量,此時(shí)的持續(xù)產(chǎn)量和對(duì)應(yīng)的捕撈努力量將隨著年增殖放流量的增加而持續(xù)升高。
由此可以知,增殖剩余產(chǎn)量模型對(duì)應(yīng)的MSY 曲線將不再是公式(2)得到的拋物線,而是1 條偏態(tài)曲線。其中,產(chǎn)量與生物量的關(guān)系曲線為負(fù)偏態(tài)(圖4),而產(chǎn)量與捕撈努力量的關(guān)系曲線為正偏態(tài)(圖5)。表示在增殖放流背景下,獲得MSY 時(shí)的生物量低于無(wú)增殖放流時(shí)的情況,而獲得MSY 時(shí)的捕撈努力量則高于無(wú)增殖放流時(shí)(圖4,圖5)。圖4 顯示,增殖剩余產(chǎn)量模型的曲線與縱軸的交點(diǎn)大于0,即理論上即使海域中該漁業(yè)資源的生物量為0,也可以依靠增殖放流獲得一定的產(chǎn)量。由圖5 可知,當(dāng)產(chǎn)量達(dá)到MSY后,繼續(xù)增加捕撈努力量,由于增殖放流對(duì)生物量的補(bǔ)充作用,持續(xù)產(chǎn)量曲線下降速度較為緩慢。圖4 曲線中對(duì)應(yīng)的BMSY和MSY 可使用公式(6)和(7)計(jì)算得到,而圖5 中曲線中的EMSY可通過(guò)求解公式(10)中Et的解而得到。當(dāng)無(wú)增殖放流時(shí),取St=0 即可。
公式(6)和(7)中相應(yīng)參數(shù)的求解方法與傳統(tǒng)模型相似,可以使用觀測(cè)誤差法或者過(guò)程誤差法。其中觀測(cè)誤差法假設(shè)生物量動(dòng)態(tài)變化模型是準(zhǔn)確的,誤差出現(xiàn)在單位捕撈努力量漁獲量的觀測(cè)值上,其參數(shù)求算過(guò)程如下:
圖4 有無(wú)增殖放流情況下產(chǎn)量?生物量關(guān)系示意圖Fig.4 Schematic diagram of the relationship of yield-biomass with or without enhancement
圖5 有無(wú)增殖放流情況下產(chǎn)量?捕撈努力量關(guān)系示意圖Fig.5 Schematic diagram of the relationship of yield-fishing effort with or without enhancement
過(guò)程誤差法假設(shè)各觀測(cè)數(shù)據(jù)是準(zhǔn)確的,誤差源于生物量動(dòng)態(tài)變化模型,參數(shù)求算過(guò)程如下:
公式(13)和(16)分別為觀測(cè)誤差法和過(guò)程誤差法的目標(biāo)函數(shù),通過(guò)最大似然法可對(duì)上述r、k、e和q等參數(shù)進(jìn)行求解,進(jìn)而使用公式(6)和(7)求算相應(yīng)的BMSY和MSY。
通過(guò)計(jì)算,當(dāng)沒(méi)有增殖放流時(shí),MSY 約為12 t,對(duì)應(yīng)的BMSY和EMSY分別約為和34 h。而當(dāng)增殖放流存在時(shí),MSY 的值也有所增加,即增殖放流群體提供了更多的產(chǎn)量(以圖5 為例,年增殖放流苗種0.48 t 時(shí),在BMSY和EMSY分別約等于71.1 t 和45 h 時(shí),得到的MSY 約為12.7 t,)。因此,增殖放流對(duì)承受更大的捕撈強(qiáng)度以及獲得更多的產(chǎn)量起到了積極的作用。
當(dāng)假定海域中初始資源生物量為80 t,且每年產(chǎn)量恒定在9 t 時(shí),無(wú)增殖放流情況下的漁業(yè)資源生物量緩慢增加,最大值無(wú)法達(dá)到K的水平,而是趨于104 t左右。每年增殖放流苗種穩(wěn)定在0.48 t 時(shí)(按照苗種平均體重2 g、成體平均體重100 g、有效增殖率0.245計(jì)算,將有5.88 t 的成體可供捕撈),生物量快速增加,并逐漸穩(wěn)定在124 t 左右(圖6a)。當(dāng)增殖放流量從每年0.2 t 開始逐年增加時(shí),雖然開始生物量的增長(zhǎng)率不是最高,但生物量在持續(xù)地緩慢增加,并且隨著放流數(shù)量不斷加大,生物量逐漸逼近K(圖6a)。若初始生物量不變,年產(chǎn)量增加到12 t 時(shí),無(wú)增殖放流情況下到第20 年時(shí),資源生物量將降至0。年增殖放流苗種穩(wěn)定在0.48 t 的情況,生物量緩慢增加至94 t左右。而當(dāng)增殖放流數(shù)量逐漸增加,在開始的幾年里,由于放流量較小,無(wú)法彌補(bǔ)產(chǎn)量,而出現(xiàn)生物量下降的趨勢(shì),但隨著增殖放流量的逐年增加,生物量開始攀升,并逐漸逼近K(圖6b)。由此可見,不同增殖放流量和產(chǎn)量水平的組合,將直接影響海域中漁業(yè)資源生物量的水平,以及MSY 的大小。
圖6 捕撈作用下不同增殖放流數(shù)量對(duì)漁業(yè)資源生物量的影響Fig.6 Effects of different stock enhancement numbers on fishery resource biomass under fishing
上述呈現(xiàn)的是將海域中原存群體與增殖群體結(jié)合在一起考慮的結(jié)果??梢韵胂?,實(shí)施增殖放流后,增殖群體將對(duì)原存群體產(chǎn)生抑制作用。若將兩者分開考慮,單獨(dú)分析海域中原存群體,可以更加清晰地看到這種抑制作用。增殖剩余產(chǎn)量模型(公式(3))中,Bt表示第t年原存群體生物量(第t?1 年增殖放流群體中存活下來(lái)的個(gè)體與當(dāng)年原存群體合并后,成為第t年的原存群體),而St表示第t年增殖放流苗種長(zhǎng)到可捕成體規(guī)格時(shí)的生物量。圖7 顯示了年增殖放流量的變化對(duì)海域中原存群體所能提供產(chǎn)量的影響情況。當(dāng)無(wú)增殖放流時(shí),MSY 出現(xiàn)在K/2時(shí),對(duì)應(yīng)的EMSY約為34 h(圖7)。而隨著增殖放流量的增加,原存群體能提供的MSY 越來(lái)越小,并且對(duì)應(yīng)的EMSY也逐漸降低(年增殖放流苗種為0.36 t、0.48 t 和0.60 t時(shí),原存群體對(duì)應(yīng)的MSY 和EMSY分別為10.1 t 和24 h、9.5 t 和21 h 以及9.2 t 和20 h)(圖7)。原存群體出現(xiàn)MSY時(shí)對(duì)應(yīng)的生物量仍舊是其對(duì)應(yīng)的Ki/2(Ki代表不同年增殖放流數(shù)量情況下,原存群體所對(duì)應(yīng)的環(huán)境容納量,并且Ki<K)。由此可知,增殖放流導(dǎo)致原存群體的K變小,顯示了增殖放流對(duì)原存群體的抑制作用。
圖7 不同增殖放流數(shù)量對(duì)海域中原存漁業(yè)資源產(chǎn)量?捕撈努力量的影響Fig.7 The effects of different stock enhancement numbers on the relationship between yield-fishing effort for the fishery resources that original existed in the sea area
本研究基于觀測(cè)誤差方法,使用增殖剩余產(chǎn)量模型對(duì)相關(guān)參數(shù)進(jìn)行求解(詳見公式(11))。不同增殖放流及捕撈策略下BMSY和MSY 的模擬估算結(jié)果如圖8 所示(刪除了部分異常估算值)??傮w而言,隨機(jī)增殖放流策略得到的BMSY和MSY 模擬估算結(jié)果變異性最小,且結(jié)果準(zhǔn)確性最高,4 種捕撈過(guò)程得到的結(jié)果基本落在真實(shí)值區(qū)間內(nèi)(圖8d)。其他3 種放流策略情況下,BMSY和MSY 的評(píng)估結(jié)果分別有被低估和高估的趨勢(shì),但程度不一致(圖8a 至圖8c),其中產(chǎn)量持續(xù)增加的捕撈過(guò)程得到的模擬評(píng)估結(jié)果中值的誤差最大,其余3 種捕撈過(guò)程模擬評(píng)估結(jié)果的中值與真實(shí)值較為接近。穩(wěn)定和持續(xù)增加的增殖放流兩種策略中,產(chǎn)量持續(xù)增加的捕撈策略對(duì)BMSY的評(píng)估結(jié)果均低于其他捕撈策略(圖8a1和8b1),而MSY 的估算結(jié)果則均明顯高于其他捕撈策略(圖8a2 和8b2)。對(duì)于持續(xù)減少的增殖放流策略,當(dāng)產(chǎn)量持續(xù)增加時(shí),MSY 的估算結(jié)果雖然較其他捕撈策略更為離散,但中值差異不大(圖8c2)。對(duì)比4 種增殖策略,持續(xù)增加的增殖放流得到了最小的BMSY中值估計(jì)值(圖8b1)和最大MSY 中值估計(jì)值(圖8b2)。
圖9 展示了不同的生物量與放流量及捕撈努力量與放流量組合水平下的MSY 估計(jì)值。理論上,即使生物量等于0,當(dāng)有增殖放流的持續(xù)補(bǔ)充時(shí),也能維持低水平的MSY(圖9a)。隨著生物量增加,MSY水平也隨之增加。當(dāng)生物量在環(huán)境容納量的30%~60%(圖9a 中B在50~100 t)時(shí),MSY 保持在較高水平。當(dāng)生物量繼續(xù)增加,MSY 開始下降。當(dāng)生物量接近環(huán)境容納量時(shí),MSY 趨于0。隨著捕撈努力量的增加,不同增殖放流量水平的MSY 均呈現(xiàn)先增高后下降的趨勢(shì)(圖9b)。在同一捕撈努力量水平下,增殖放流量增加時(shí),MSY 先降低再升高,并且隨著增殖放流量的持續(xù)增加而不斷升高(圖9b)。當(dāng)增殖放流量逐漸增加,得到MSY 所對(duì)應(yīng)的生物量和捕撈努力量變化不明顯。
圖8 不同增殖和捕撈策略下BMSY(左)和MSY(右)的模擬估算結(jié)果Fig.8 The simulation estimate results of BMSY(left) and MSY (right) under different stock enhancement and fishing strategies
圖9 不同生物量和捕撈努力量與不同年增殖放流量水平組合下MSY 的模擬估算值Fig.9 The simulation estimates of MSY under combinations of different biomass and fishing effort levels with different stock enhancement levels
本研究在傳統(tǒng)Schaefer 剩余產(chǎn)量模型基礎(chǔ)上,將漁業(yè)資源生物量邏輯斯蒂增長(zhǎng)特征以及增殖放流對(duì)生物量的影響相結(jié)合,提出了適用于定棲性種類的增殖剩余產(chǎn)量模型,并模擬不同的捕撈和增殖策略情況,對(duì)其MSY 等特征值進(jìn)行了估算。該模型與傳統(tǒng)Schaefer模型最大的不同就在于,增加了表示增殖效果的參數(shù)有效增殖率e。增殖剩余產(chǎn)量模型中,r與e兩個(gè)參數(shù)均用于描述漁業(yè)資源生物量的增長(zhǎng)率情況。雖然r和e所代表的原存群體和增殖群體最終的增長(zhǎng)結(jié)果都通過(guò)Bt來(lái)體現(xiàn),但當(dāng)?shù)趖+1 年新的增殖發(fā)生之前,模型是將原存群體和增殖群體分開考慮。此外,值得注意的是,本研究在不同增殖和捕撈策略下模擬得到的參數(shù)r與e的相關(guān)系數(shù)平均值等于?0.68,相關(guān)性達(dá)到了統(tǒng)計(jì)學(xué)顯著性水平(p<0.05)。由此可知,參數(shù)r和e兩者之間存在顯著的負(fù)相關(guān),即原存群體和增殖群體兩者的增長(zhǎng)存在相互抑制作用,兩個(gè)群體中一個(gè)增長(zhǎng)加快必然導(dǎo)致另一個(gè)增長(zhǎng)減慢。因此,在使用增殖剩余產(chǎn)量模型時(shí),若能事先獲得關(guān)于參數(shù)r和e的相關(guān)信息,或者固定其中1 個(gè)參數(shù)值,將有助于獲得準(zhǔn)確的參數(shù)估算結(jié)果。
增殖剩余產(chǎn)量模型對(duì)參數(shù)r和e的表達(dá)以及原存群體和增殖群體兩者的相互抑制性,關(guān)系到如何對(duì)兩個(gè)群體開展分析的問(wèn)題,即是將兩個(gè)群體綜合在一起考慮,還是將兩者分別考慮。增殖剩余產(chǎn)量模型的最終結(jié)果是將兩個(gè)群體綜合在一起呈現(xiàn)的。圖6 表明,增殖群體確實(shí)能夠使產(chǎn)量持續(xù)并維持在較高的水平,特別是當(dāng)產(chǎn)量過(guò)高,生物量開始下降時(shí),此時(shí)增加放流量能夠?qū)σ虿稉贫档偷纳锪科鸬窖a(bǔ)償作用。然而,如果將兩個(gè)群體分開考慮,則能夠清晰看到兩者之間相互抑制的效應(yīng)。如圖7 所示,隨著每年增殖量的增加,原存群體所維持的MSY 水平逐漸降低,并且EMSY也隨之減小,說(shuō)明增殖放流實(shí)施后,原存群體提供產(chǎn)量的作用被逐漸弱化。以無(wú)增殖放流和每年增殖放流0.48 t 苗種兩種情況為例,當(dāng)沒(méi)有增殖時(shí),原存群體的MSY=11.91 t,相應(yīng)的BMSY=84.13 t。而當(dāng)每年增殖放流0.48 t 苗種以后,總的MSY=12.72 t,此時(shí)原存群體的BMSY=71.14 t。假如在無(wú)增殖放流時(shí),把生物量維持在71.14 t(即年放流0.48 t 時(shí)原存群體的BMSY)的水平上,此時(shí)MSY=9.51 t,說(shuō)明增殖群體彌補(bǔ)了3.21 t 的產(chǎn)量。而這3.21 t 的產(chǎn)量來(lái)自增殖群體和原存群體兩方面,即增殖苗種生長(zhǎng)至可供捕撈尺寸時(shí)的群體以及從84.13 t 降到71.14 t 后減少的12.99 t的原存群體,這兩部分加在一起提供了比無(wú)增殖時(shí)的MSY 多出來(lái)的3.21 t。
在增殖放流作用影響下,通常高產(chǎn)量情況下評(píng)估得到的MSY 值要高于其他情況,特別是增殖放流量持續(xù)增加的時(shí)候(圖8)。這主要受到模型求解參數(shù)結(jié)果的影響。增殖剩余產(chǎn)量模型對(duì)參數(shù)r和e的求算均有不同程度偏高的趨勢(shì),其中,參數(shù)e的偏高程度略大于r。當(dāng)增殖放流量持續(xù)增加時(shí),生物量不斷增長(zhǎng),因此參數(shù)r和e高估趨勢(shì)增加,導(dǎo)致MSY 值變大(見公式(7))。又由于e的高估程度大于r,因此BMSY值相對(duì)變?。ㄒ姽剑?))。對(duì)于產(chǎn)量逐年增加的情況,當(dāng)增殖放流量也持續(xù)增加時(shí),此時(shí)體現(xiàn)生物量的變化率較大,因此參數(shù)r和e的估計(jì)值較大,MSY也較大;而當(dāng)增殖放流量持續(xù)減少時(shí),生物量的變化率減小,則參數(shù)r和e估計(jì)偏高程度降低,MSY 也相應(yīng)減?。▓D8b,圖8c)。隨機(jī)增殖放流和產(chǎn)量先增后降的策略得到的BMSY和MSY 模擬估算結(jié)果變異性最?。▓D8),說(shuō)明同類數(shù)據(jù)的對(duì)比性有利于模型評(píng)估得到穩(wěn)定的結(jié)果。對(duì)比圖8b 和8c 可發(fā)現(xiàn),當(dāng)增殖放流與產(chǎn)量變化趨勢(shì)一致時(shí),得到MSY 的評(píng)估結(jié)果變異性較小,而相反時(shí)則較大。說(shuō)明非同類數(shù)據(jù)變化相一致時(shí),有助于提高模型參數(shù)估算的穩(wěn)定性。
從資源整體來(lái)看,實(shí)施增殖放流量的大小對(duì)MSY水平的影響并不明顯(圖9a),表明增殖剩余產(chǎn)量模型中eS所表示的增殖放流群體可補(bǔ)充到資源群體并被捕撈利用的數(shù)量基本是穩(wěn)定的,即S增大則e將減小。少量的增殖放流,效果反而不如無(wú)增殖放流(圖9b),顯示了此時(shí)增殖放流對(duì)原存群體的抑制作用大于促進(jìn)作用。當(dāng)增殖放流量逐漸增加時(shí),促進(jìn)作用開始顯現(xiàn),此時(shí)生物量的變化受增殖放流作用的影響越來(lái)越顯著。
本研究中對(duì)參數(shù)的估算使用了觀測(cè)誤差法,實(shí)踐證明觀測(cè)誤差法具有一定的優(yōu)越性[9?10]。然而,誤差的來(lái)源形式是多種多樣的,不僅包括過(guò)程誤差和觀測(cè)誤差,還存在許多其他形式的誤差類型[11]。傳統(tǒng)剩余產(chǎn)量模型中的參數(shù)(r、K和q)之間存在著高度的相關(guān)性[12],本研究提出的增殖剩余產(chǎn)量模型,參數(shù)之間同樣存在相關(guān)性。事實(shí)上,在許多漁業(yè)資源評(píng)估的模型中都存在參數(shù)相關(guān)的情況,這難免對(duì)參數(shù)估算結(jié)果產(chǎn)生影響,但影響程度有多大,則需進(jìn)一步的分析與探討。此外,增殖剩余產(chǎn)量模型中相應(yīng)變量,如B、S和C之間也會(huì)產(chǎn)生相互影響,特別是每年的增殖放流量發(fā)生較大變化的時(shí)候,通過(guò)公式(5)得到的BMSY求算結(jié)果將會(huì)存在誤差。因此,公式(6)的結(jié)果應(yīng)以每年的增殖放流量穩(wěn)定為前提,或者使用數(shù)值模擬的方式估算MSY,以避免多變量之間復(fù)雜的關(guān)系。
為了進(jìn)一步豐富剩余產(chǎn)量模型的應(yīng)用,許多學(xué)者在傳統(tǒng)模型基礎(chǔ)上演變出了系列模型[13-16]。而在漁業(yè)資源評(píng)估工作不斷發(fā)展過(guò)程中,模型的發(fā)展也是內(nèi)在的需求。本研究提出了一個(gè)適用于增殖放流背景下的定棲性種類剩余產(chǎn)量模型。選擇定棲性種類,主要是考慮其移動(dòng)性較弱,當(dāng)增殖放流實(shí)施后,不會(huì)出現(xiàn)棲息海域范圍發(fā)生變化,從而引起環(huán)境容納量改變的情況,這樣就使模型參數(shù)不至于過(guò)于復(fù)雜。受此條件的限制,該模型比較適用于海洋牧場(chǎng)內(nèi)增殖的定棲性種類,或者那些在比較清晰的作業(yè)邊界內(nèi),且不會(huì)超出邊界范圍的增殖目標(biāo)種。在今后的研究中,我們將針對(duì)這一點(diǎn)對(duì)模型作進(jìn)一步的調(diào)整,使其適用范圍更加廣泛。與傳統(tǒng)剩余產(chǎn)量模型相似,增殖剩余產(chǎn)量模型對(duì)參數(shù)的評(píng)估結(jié)果容易受到參數(shù)所設(shè)置初始值的影響,在后續(xù)的研究中,我們將考慮開展敏感性評(píng)估和根據(jù)實(shí)際觀測(cè)進(jìn)行驗(yàn)證。此外,模型中新增加的參數(shù)e會(huì)受到其他因素的影響,比如環(huán)境因子、人工魚礁集魚效果等等,這些在本研究中未進(jìn)行細(xì)化,因此也同樣將是我們未來(lái)研究的重點(diǎn)內(nèi)容。