氣候變化對東北大西洋漁獲物組成、多樣性和營養(yǎng)級的影響

陳爽，陳新軍,2,3,4,5

( 1. 上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海 201306；2. 青島海洋科學(xué)與技術(shù)試點國家實驗室 海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室，山東 青島 266237；3. 農(nóng)業(yè)農(nóng)村部 大洋漁業(yè)開發(fā)重點實驗室，上海 201306；4. 國家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心，上海 201306；5. 大洋漁業(yè)資源可持續(xù)開發(fā)教育部重點實驗室，上海 201306)

1 引言

氣候變化對海洋漁業(yè)的影響是一個復(fù)雜的生態(tài)過程，魚類對海洋環(huán)境的改變有最直接的響應(yīng)，其中包括魚類個體的生長、死亡、繁殖、洄游分布等[1–3]。不僅如此，海洋魚類作為海洋生態(tài)系統(tǒng)中最重要的組成部分，其種群資源的波動和分布范圍的變化也會間接影響生態(tài)系統(tǒng)的食物鏈構(gòu)成、物種的多樣性和生物營養(yǎng)級水平[4–5]。在近年來的研究中，更多時空尺度上的氣候變化現(xiàn)象，例如太平洋十年濤動（Pacific Decadal Oscillation, PDO）、厄爾尼諾?南方濤動（El Ni?o-Southern Oscillation, ENSO）等，也給海洋漁業(yè)資源的可持續(xù)開發(fā)與科學(xué)管理帶來了新的挑戰(zhàn)[6–7]。

東北大西洋海域連接北冰洋和赤道地區(qū)，跨越了亞熱帶、溫帶和寒帶，包括大范圍的深海區(qū)，大陸架淺海和大量河口、海灣，地理位置優(yōu)越，適合各種類型的魚類生長，形成了許多重要的漁場，也是世界上漁業(yè)資源最豐富的地區(qū)。根據(jù)政府間氣候變化專門委員會（IPCC）報告，在全球變暖的背景下，處于高緯度地區(qū)的東北大西洋氣候變化速率將遠(yuǎn)遠(yuǎn)高于全球平均水平，其漁業(yè)生產(chǎn)力將面臨更大的挑戰(zhàn)[8]。對東北大西洋漁業(yè)資源，大部分學(xué)者集中于氣候變化對魚類資源量、洄游分布和群落結(jié)構(gòu)等方面的研究，例如Sirabella[9]運用主成分分析和典范對應(yīng)分析發(fā)現(xiàn)在北大西洋濤動（North Atlantic Oscillation, NAO）發(fā)生時，海水溫度與當(dāng)年大西洋鱈（Gadus morhua）的補(bǔ)充群體的生物量關(guān)系最為密切；Drinkwater[10]通過預(yù)測模型發(fā)現(xiàn)水溫與北海、愛爾蘭海、凱爾特海的大西洋鱈魚群體補(bǔ)充量呈負(fù)相關(guān)，與巴倫支海、冰島的大西洋鱈魚群體補(bǔ)充量呈正相關(guān)；在洄游分布方面，Perry等[11]發(fā)現(xiàn)北海漁場多個魚類種群都發(fā)生顯著的向北移動，Rose[12]發(fā)現(xiàn)了水溫與毛鱗魚（Osmerus mordax）移動的距離存在一個用對數(shù)函數(shù)描述的非線性關(guān)系，水溫每上升1℃，毛鱗魚分布范圍在緯度方向上向北移動275 km；在群落變化上，F(xiàn)ossheim[13]發(fā)現(xiàn)近年來隨著氣候變暖，巴倫支海原有的魚類群落發(fā)生了顯著變化，出現(xiàn)了許多來自低緯度海域的魚類，海域的生產(chǎn)力也顯著提高；Stenevik和Sundby[14]發(fā)現(xiàn)北海的魚類群落中小型中上層魚類正在慢慢占據(jù)優(yōu)勢地位，主要包括大西洋鯡（Clupea harengus）和大西洋鯖（Scomber scombrus）等。近年來東北大西洋漁業(yè)資源變化受氣候變化的影響較明顯，故本研究基于聯(lián)合國糧食與農(nóng)業(yè)組織（FAO）提供的漁獲物統(tǒng)計數(shù)據(jù)，從漁獲物組成水平（多樣性和營養(yǎng)級）探討東北大西洋海域漁業(yè)生產(chǎn)的結(jié)構(gòu)演變，尋求它們與氣候、海洋環(huán)境因素之間的相互聯(lián)系，為進(jìn)一步研究東北大西洋漁業(yè)資源與氣候變化的關(guān)系提供科學(xué)依據(jù)。

2 材料與方法

2.1 數(shù)據(jù)來源

2.1.1 漁業(yè)數(shù)據(jù)

東北大西洋漁業(yè)數(shù)據(jù)選取FAO每年統(tǒng)計的漁獲生產(chǎn)數(shù)據(jù)，獲取網(wǎng)址為http://www.fao.org/fishery/statistics/global-capture-production/query/en，生產(chǎn)捕撈區(qū)域為FAO劃分的27漁區(qū)（36°N～90°N，42°W～69°E），包括了東北大西洋海域和北冰洋部分海域，時間跨度為1982?2016年，漁獲產(chǎn)量單位為t。

2.1.2 環(huán)境數(shù)據(jù)和氣候指數(shù)

本研究中環(huán)境數(shù)據(jù)選取海表面溫度（SST）、海表面鹽度（SSS）、海表面高度（SSH）、海冰密集度（Sea Ice Concentration, SIC），這些數(shù)據(jù)來源于亞太數(shù)據(jù)研究中心網(wǎng)站（http://apdrc.soest.hawaii.edu/index.php），時間分辨率為月，空間分辨率為1°×1°，時間跨度為1982?2016年。

東北大西洋地理環(huán)境特殊，存在有多重時間尺度的海洋氣候變化。本研究選取了該海域最典型的氣候變化特征NAO、大西洋年代際濤動（Atlantic Multidecadal Oscillation, AMO）和北極濤動（Actical Oscillation, AO）作為氣候變化指標(biāo)。氣候變化指數(shù)來源于美國國家海洋和大氣管理局網(wǎng)站（https://www.esrl.noaa.gov/psd/data/climateindices/list/），時 間 分 辨 率 為月，選取時間跨度為1982?2016年。

2.2 研究方法

2.2.1 數(shù)據(jù)預(yù)處理

根據(jù)1982?2016年漁獲產(chǎn)量統(tǒng)計數(shù)據(jù)共得到399種物種，其中有海洋魚類308種，甲殼類生物39種，軟體類生物52種。相關(guān)漁獲物營養(yǎng)級參考世界Fishbase數(shù)據(jù)庫（http://www.fishbase.org/search.php）和Sealifebase數(shù)據(jù)庫（http://www.sealifebase.org/search.php）中所提供的營養(yǎng)級信息。

將得到的環(huán)境數(shù)據(jù)和氣候指數(shù)按年進(jìn)行平均，得到東北大西洋各年的環(huán)境數(shù)據(jù)和氣候指數(shù)平均值。

2.2.2 生物多樣性指數(shù)

生物多樣性指數(shù)能夠反映組成結(jié)構(gòu)的均勻度和豐富度，本研究選用Simpson多樣性指數(shù)（D）做為指標(biāo)[15]，其計算公式為

式中，Yi為每年第i個物種的漁獲產(chǎn)量，Y為當(dāng)年的總產(chǎn)量，D的值介于0～1之間。

2.2.3 平均營養(yǎng)級

漁獲物的平均營養(yǎng)級水平（MTL）可以作為漁業(yè)生態(tài)系統(tǒng)結(jié)構(gòu)和功能變化的重要指標(biāo)[16]。本研究中漁獲物平均營養(yǎng)級依據(jù)Pauly等[17]的方法進(jìn)行計算，其公式為

式中，TLi為第i個物種的營養(yǎng)級。

2.2.4 主成分變化分析

本研究中1982?2016年共出現(xiàn)有記錄的商業(yè)性捕撈魚種308種，對年平均產(chǎn)量10000 t以上的魚種進(jìn)行篩選，共有包括甲殼類和軟體類在內(nèi)的漁獲物種類46種（表1）。由于考慮到漁獲物產(chǎn)量的變化會受到人為捕撈的影響，與之相比，漁獲物產(chǎn)量所占比例在捕撈干擾下相對穩(wěn)定，因此選用主要的漁獲物產(chǎn)量所占比例來進(jìn)行漁獲物組成變化分析。

表1 漁獲物組成第1和第2主成分載荷值Table 1 Loadings on the first and second principal components from the analysis of catch composition

主成分分析的方法能夠?qū)Χ嗑S的變量進(jìn)行降維處理，提取能夠反映變量間相互關(guān)系的共同特征，在漁業(yè)生態(tài)領(lǐng)域有較廣泛的應(yīng)用[18]。對漁獲量占比、環(huán)境和氣候因子分別進(jìn)行主成分分析，提取漁獲物組成和東北大西洋氣候環(huán)境變化的主成分變化特征值（Principal Component Score，PCs），并檢驗它們是否存在相關(guān)性，具有較高相關(guān)性的主成分特征值視為環(huán)境和氣候變化下漁獲物組成變化特征。

續(xù)表1

2.2.5 漁獲物組成和氣候環(huán)境因素之間的關(guān)系

將2.2.4中與環(huán)境、氣候因子主成分特征值具有較高相關(guān)性的漁獲物組成變化主成分特征值作為指標(biāo)，各環(huán)境因子和氣候指數(shù)作為解釋變量，并建立廣義可加模型（Generalized Additive Models，GAM），分析漁獲物組成變化與各環(huán)境因子和氣候因子間的關(guān)系。GAM模型的表達(dá)式為

式中，g為連接函數(shù)，α為模型的截距，s為非線性的平滑函數(shù)，X1、X2、···、Xi為解釋變量，ε為誤差項。

采用基于赤池信息準(zhǔn)則（Akaike Information Criterion, AIC）的逐步回歸方法檢驗?zāi)Ｐ偷臄M合程度，一般AIC值越小，擬合程度越好。AIC值計算公式為

式中，k表示模型獨立參數(shù)個數(shù)；L表示模型的極大似然函數(shù)。

3 結(jié)果

3.1 漁獲物組成、多樣性和營養(yǎng)級變化特征

由圖1可知，東北大西洋漁獲物多樣性指數(shù)（Diversity Index, DI）在1990年之前呈持續(xù)上升的趨勢，在1990年達(dá)到最大值0.92，1990?2002年間開始緩慢下降，2002年之后急劇下降，并且在2002?2010年間都處于較低水平，2010年之后趨于平穩(wěn)（圖1）。平均營養(yǎng)級水平（Mean Trophic Level, MTL）則在2002年之前未出現(xiàn)較大波動，一直處于較低水平，平均值為3.59，2002年后開始急劇上升，2002?2010年間出現(xiàn)了與同時期DI相反的變化，2010年后平均營養(yǎng)級持續(xù)升高（圖1）。

圖1 1982?2016年多樣性指數(shù)和平均營養(yǎng)級變化Fig. 1 The variations of diversity index and mean trophic level during 1982 to 2016

對46種漁獲物的產(chǎn)量占比進(jìn)行主成分分析，得到前兩位的主成分方差累積解釋率達(dá)到52.4%。其中第一位、第二位主成分的解釋率分別為35.3%、17.1%（圖2）。通過主成分分析載荷值可知（表1），以黍鯡、大西洋鯡為主的小型中上層魚類和以挪威長臀鱈、牙鱈、青鱈為主的中下層魚類是第一主成分變化的主要貢獻(xiàn)者。而大西洋鱈、大西洋鯖等魚類是第二主成分變化的主要貢獻(xiàn)者。

圖2 1982?2016年漁獲物組成主成分得分年際變化Fig. 2 The variations of the principal component scores for the catch composition during 1982 to 2016

3.2 漁獲物組成與氣候環(huán)境因素之間的關(guān)系

對東北大西洋SST、SSS、SSH、SIC、NAO、AMO和AO等氣候、環(huán)境因子進(jìn)行主成分分析，其前兩位的累計解釋率為78.7%，其中第一主成分PCI解釋率為50.7%，SST、SSH等環(huán)境因子的載荷值較高，顯示了區(qū)域環(huán)境因素的變化特征；第二主成分PCII的解釋率為28.0%，以AO、NAO為主的長時間尺度氣候變化具有較高的載荷值，與氣候變化較為相關(guān)（圖3，表2）。將漁獲物組成變化特征DI、MTL、PC1、PC2與該海域的氣候環(huán)境主成分變化特征PCI、PCII進(jìn)行相關(guān)性檢驗后發(fā)現(xiàn)，漁獲物組成變化PC1與氣候環(huán)境變化PCI、PCII均呈顯著相關(guān)（p＜0.05），DI和MTL都與PCI具有顯著相關(guān)性，而PC2與氣候環(huán)境變化PCI、PCII相關(guān)性并不顯著（表3）。因此DI和MTL與海洋環(huán)境因素有較高的相關(guān)性，與NAO、AO等氣候因素并不相關(guān)。漁獲物組成PC1則可以較好的解釋氣候、環(huán)境變化對其的影響。

圖3 1982?2016年東北大西洋氣候、環(huán)境因素主成分得分年際變化Fig. 3 The variations of the principal component scores for the climatic and environmental factors during 1982 to 2016

表2 氣候、環(huán)境因素第1和第2主成分載荷值Table 2 Loadings on the first and second principal components from the analysis of the climatic and environmental factors

表3 氣候、環(huán)境因素主成分變化與各類漁獲物組成指數(shù)相關(guān)性系數(shù)Table 3 The correlation coefficients of climatic and environmental factors principal component scores and different kinds of catch composition indexes

通過GAM模型對漁獲物組成變化PC1與海域氣候、環(huán)境因子進(jìn)行擬合，結(jié)果如表4所示，選擇AIC值最小的模型為最優(yōu)模型，因此最終模型可表示為：PC1～s(SST)+s(SSS)+s(SSH)+s(SIC)+s(NAO)。

該模型對漁獲物組成PC1的累計解釋率達(dá)到94.3%，其中貢獻(xiàn)率最大的因子是SST，解釋率為74.8%，其他有影響的環(huán)境因子包括SSH、SSS、SIC，氣候因素中NAO也對PC1有影響。

表4 漁獲物組成PC1與各氣候、環(huán)境因素的GAM模型擬合Table 4 GAM models fitted to the first principal component scores of catch composition (PC1) and climatic and environmental factors

由氣候、環(huán)境因子與漁獲物組成變化特征PC1之間的關(guān)系（圖4）可知，SST與PC1呈現(xiàn)負(fù)相關(guān)，海表水溫在5.8～6.2℃之間時，漁獲物組成變化幅度較明顯；SSS與PC1幾乎呈線性負(fù)相關(guān)，SSS為33.55°時置信區(qū)間最??；SSH在?0.73～?0.71 m時對漁獲物組成影響并不顯著，當(dāng)SSH達(dá)到?0.71 m甚至更高時漁獲物結(jié)構(gòu)會受到較大影響；SIC和NAO對PC1的影響趨勢相似，幾乎都呈線性正相關(guān)。

圖4 氣候、環(huán)境因素對漁獲物組成PC1的影響Fig. 4 Effects of climatic and environmental factors on the first principal component scores of catch composition (PC1)

4 討論與分析

4.1 漁獲物多樣性和平均營養(yǎng)級的變化

不同漁獲物組成指標(biāo)的變化能夠反映東北大西洋1982?2016年間漁獲物組成大致經(jīng)歷了3個時期的變化。第1時期為1982?2002年，DI和MTL都在波動中呈現(xiàn)輕微的下降趨勢；第2時期為2002?2010年， DI和MTL在該時期分別發(fā)生了急劇的下降和上升，并在2006年后開始恢復(fù)；第3個時期為2010年之后DI趨于平穩(wěn)，MTL在該時期不斷升高。

MTL的變化往往被認(rèn)為是人為捕撈和環(huán)境變化共同影響之下的結(jié)果[19]，Pauly等[20]發(fā)現(xiàn)全球漁獲物平均營養(yǎng)級從20世紀(jì)70年代以來呈現(xiàn)穩(wěn)步的下降趨勢，這種變化遵循“捕撈對象沿著海洋食物網(wǎng)向下級移動（Fishing Down Marine Food Webs）”的模式。東北大西洋漁獲物MTL在1982?2002年間持續(xù)下降的現(xiàn)象則與這種模式相符，根據(jù)FAO的統(tǒng)計報告，東北大西洋是漁業(yè)資源過度開發(fā)最嚴(yán)重的海區(qū)，自20世紀(jì)50年代以來，大西洋鱈、黑線鱈、大西洋鯡等魚種經(jīng)歷了高強(qiáng)度的捕撈，導(dǎo)致資源衰退甚至多個種群面臨崩潰，同時玉筋魚、毛鱗魚等低營養(yǎng)級的漁業(yè)資源開始開發(fā)，漁獲量不斷增加[21]。研究發(fā)現(xiàn)，2006年之后的MTL出現(xiàn)了顯著的增長，這一變化過程和該海域漁業(yè)資源恢復(fù)與新物種的開發(fā)有著必然的聯(lián)系，嚴(yán)格的捕撈控制和氣候變暖使東北大西洋鱈科魚類逐漸出現(xiàn)資源恢復(fù)的跡象[22–23]，隨著捕撈活動增加，傳統(tǒng)的高營養(yǎng)級魚類產(chǎn)量占比也開始回升。2003?2007年間對藍(lán)鱈的高強(qiáng)度捕撈也是MTL不斷上升的主要因素之一[24]。

東北大西洋漁獲物組成多樣性整體上呈不斷下降的趨勢，這一變化受多種因素影響。高強(qiáng)度的人為捕撈導(dǎo)致漁獲物的產(chǎn)量和種類不斷減少，往往會造成組成結(jié)構(gòu)的多樣性降低[25–26]。尤其是上世紀(jì)對該海域主要經(jīng)濟(jì)魚種鱈魚、鯡魚等資源的過度開發(fā)，導(dǎo)致這些魚類難以再形成漁汛，這些主要捕撈對象在總漁獲物中的占比不斷下降，而失去天敵的中小型上層魚類迅速成為優(yōu)勢魚種并愈發(fā)突出，漁獲物占比失衡嚴(yán)重，造成多樣性不斷降低的局面。除此以外氣候變化也會在一定程度上影響群落的穩(wěn)定性。

4.2 氣候變化對漁獲物組成的影響

主成分分析能夠從多個變量中提取共同的變化特征，將變量間復(fù)雜多變的結(jié)構(gòu)轉(zhuǎn)化為幾組特征值來表征，在群落變化的研究中有較廣泛的應(yīng)用[18,27]。在本研究中，漁獲物結(jié)構(gòu)PC1主要表征大西洋鯡、黍鯡等小型中上層魚類和挪威長臀鱈、牙鱈、青鱈、平鲉等中下層高營養(yǎng)級魚類的變化情況，PC2則主要表征了大西洋鱈、大西洋鯖、歐洲無須鱈這三種經(jīng)濟(jì)魚類的變化。東北大西洋海域氣候、環(huán)境因子的主成分分析結(jié)果顯示，PCI主要表示環(huán)境因子的變化特征，PCII主要表示NAO、AO為主的年際氣候變化現(xiàn)象。結(jié)合各主成分分析結(jié)果的特征值發(fā)現(xiàn)，PC1在1995年發(fā)生了由正值到負(fù)值的轉(zhuǎn)變，而海域環(huán)境因子變化PCI和氣候變化PCII也基本在同一時期內(nèi)發(fā)生了由負(fù)轉(zhuǎn)正和由正轉(zhuǎn)負(fù)的特征變化，相關(guān)性檢驗表明，PC1與二者均具有相關(guān)性，這一結(jié)果在某種程度上表明漁獲物結(jié)構(gòu)的主成分變化與氣候變化有必然聯(lián)系。而PC2的變化與氣候環(huán)境之間并無直接聯(lián)系，猜測PC2在1994年前后由正值轉(zhuǎn)向負(fù)值的變化，主要是長期的過度捕撈導(dǎo)致的大西洋鱈魚資源量下降造成，而2008?2009年前后由負(fù)轉(zhuǎn)正的變化與北歐國家開始大力開發(fā)鯖魚資源有關(guān)[28]。

為探究不同氣候、環(huán)境因子對漁獲物組成結(jié)構(gòu)產(chǎn)生的影響，選取GAM模型分析不同的因子的影響程度和影響機(jī)制，發(fā)現(xiàn)SST是漁獲物組成結(jié)構(gòu)變化最主要的影響因子，而SSH、SSS、SIC、NAO也在一定程度上影響結(jié)構(gòu)變化。組成結(jié)構(gòu)變化與SST最為密切，SST與PC1有顯著的負(fù)效應(yīng)，即當(dāng)水溫升高時大西洋鯡、黍鯡等小型中上層魚類在漁獲物組成中的占比上升，鱈魚、鲆鰈類中下層魚類在漁獲物組成中的占比會顯著下降，反之亦然。由于大多數(shù)的海洋魚類是變溫動物，它們對外界環(huán)境溫度的變化有最直接的反應(yīng)，適宜的升溫可以提高魚類的生長速度和食物轉(zhuǎn)化率，當(dāng)水溫超過其適宜范圍時，魚類自身的生理機(jī)能將限制其耐熱性，從而導(dǎo)致魚類為尋求適宜的棲息地而改變原有的分布范圍或者無法承受超過自身耐熱的溫度而死亡，甚至滅絕[29–30]。溫度的升高對北大西洋絕大多數(shù)冷水性魚類的生活史[31]、繁殖[32]、補(bǔ)充[10]都有著負(fù)面的影響，還會引起暖水性、溫水性魚類向北遷移，冷水性魚類的棲息地面積縮減[11]，最終導(dǎo)致該海域漁獲物組成結(jié)構(gòu)發(fā)生變化。研究表明，SSH與海洋表層水團(tuán)運輸有關(guān)，SSH變動影響浮游生物和仔稚魚運輸[33–34]，而海冰是比目魚、鱈魚產(chǎn)卵和覓食的主要場所，海冰的消融會影響這些冷水性物種的產(chǎn)量[35]。除了海洋環(huán)境因素對漁獲產(chǎn)量有著直接影響，氣候變化也是北大西洋生態(tài)系統(tǒng)中最關(guān)鍵的一環(huán)。氣候的躍變往往伴隨著生態(tài)系統(tǒng)的格局發(fā)生轉(zhuǎn)換，這種變換的機(jī)制是通過食物鏈自上而下的傳遞造成的[36–37]。在北大西洋，NAO引發(fā)寒、暖流強(qiáng)弱交替決定了該海域的浮游生物組成與豐度，對魚類補(bǔ)充量和資源變動影響深刻[38–39]。因此，NAO也可以作為東北大西洋漁獲物組成結(jié)構(gòu)變化的重要指標(biāo)。

4.3 后續(xù)研究應(yīng)考慮的問題

由于數(shù)據(jù)來源的局限性，本研究中選取東北大西洋漁獲產(chǎn)量數(shù)據(jù)作為漁獲物組成的指標(biāo)，雖具有一定的參考價值但也有一些不足，原因如下：不同的魚類對氣候變化有著不同的響應(yīng)模式，即使是在同一氣候背景或捕撈強(qiáng)度下，暖水性物種勢必在氣候變暖的大趨勢下遭受的捕撈強(qiáng)度要大于冷水性物種；其次，東北大西洋海域的漁業(yè)管理策略對許多經(jīng)濟(jì)型魚種的開發(fā)起著決定性作用，其漁獲產(chǎn)量往往與該海域歷來實施的總允許可捕量有一定的關(guān)系。雖然在后期構(gòu)建模型時，選擇使用漁獲產(chǎn)量占比變化來表征結(jié)構(gòu)組成上的差異，能減少其因捕撈強(qiáng)度等其他因素的干擾，但無法完全消除。在本研究中，選取的研究區(qū)域是FAO劃定的27漁區(qū)，這個區(qū)域既包含了北大西洋又包含了北冰洋的部分海域，由于跨越的緯度較大，因此不同區(qū)域的魚類受到的環(huán)境因素的影響程度也有區(qū)別，同樣是海冰，對于較寒冷海域的鱈和鲆鰈類產(chǎn)生的影響更為明顯，而對鯡、毛鱗魚等中上層魚類影響較顯著的主要是海表溫度，這也有可能影響研究結(jié)果。除此以外，選取的遙感數(shù)據(jù)對海洋表層與底層魚類的影響程度存在差異，在本研究中默認(rèn)具有相同程度的影響，在今后的研究中可以選擇不同溫躍層深度或海流的環(huán)境因子來更深入的進(jìn)行研究。