趙海龍, 陳新軍, 方學燕

1上海海洋大學海洋科學學院, 上海　201306

2 國家遠洋漁業(yè)工程技術(shù)研究中心, 上?！?01306

3 上海海洋大學大洋漁業(yè)可持續(xù)開發(fā)省部共建教育部重點實驗室, 上?！?01306

4 遠洋漁業(yè)協(xié)同創(chuàng)新中心, 上?！?01306

?

基于棲息地指數(shù)的東太平洋黃鰭金槍魚漁場預(yù)報

趙海龍1,4, 陳新軍1,2,3,4，*, 方學燕1,4

1上海海洋大學海洋科學學院, 上海201306

2 國家遠洋漁業(yè)工程技術(shù)研究中心, 上海201306

3 上海海洋大學大洋漁業(yè)可持續(xù)開發(fā)省部共建教育部重點實驗室, 上海201306

4 遠洋漁業(yè)協(xié)同創(chuàng)新中心, 上海201306

摘要:黃鰭金槍魚是東太平洋海域重要的金槍魚種類之一，也是我國金槍魚延繩釣的主要捕撈對象之一。根據(jù)2011年東太平洋海域(20°N—35°S、85°W—155°W)延繩釣生產(chǎn)統(tǒng)計數(shù)據(jù)，結(jié)合表溫(SST)和海面高度(SSH)的遙感數(shù)據(jù)，采用頻次分析法獲得黃鰭金槍魚分布的SST和SSH適宜范圍；運用一元非線性回歸方法，以漁獲量為適應(yīng)性指數(shù)，按季度分別建立了基于SST和SSH的長鰭金槍魚棲息地適應(yīng)性指數(shù)，采用算術(shù)平均法獲得基于SST和SSH環(huán)境因子的棲息地指數(shù)綜合模型，并用2012年各月實際作業(yè)漁場進行驗證。結(jié)果顯示，黃鰭金槍魚漁場多分布在SST為24—29℃、SSH為0.3—0.7 m的海域。采用一元非線性回歸建立的各因子適應(yīng)性指數(shù)模型在統(tǒng)計上均為顯著(P<0.05)。經(jīng)與2012年實際生產(chǎn)情況比較，作業(yè)漁場預(yù)報準確性達66%以上。研究獲得棲息地指數(shù)模型可為金槍魚延繩釣漁船尋找中心漁場提供參考。

關(guān)鍵詞：黃鰭金槍魚; 漁情預(yù)報; 東太平洋; 棲息地指數(shù); 海洋環(huán)境因子

黃鰭金槍魚(Thunnusalbacares)是一種高度洄游的大洋性魚類，在東太平洋是延繩釣的重要捕撈對象[1]。影響黃鰭金槍魚漁場分布的環(huán)境因子很多，涉及水溫、水溫垂直結(jié)構(gòu)、葉綠素和溶解氧等[2- 10]。Romena[8]指出黃鰭金槍魚的中心漁場分布受20 ℃等溫線影響。Song等[9]分析得出黃鰭金槍魚的垂直分布和溫躍層有關(guān)。但Maury等[10]指出，在延繩釣漁業(yè)中，表溫對黃鰭金槍魚的單位捕撈努力量漁獲量(CPUE) 影響較小，因為其主要分布在深層海域。綜合前人研究成果，認為影響黃鰭金槍魚漁場分布的因素是多方面的，但是目前的研究往往采用單一因子進行分析[5- 10]，各因素沒有進行有機整合，也沒有建立相應(yīng)的漁場預(yù)報模型。棲息地指數(shù)(HSI)常被用來描述野生動物的棲息地質(zhì)量，目前已被用于魚類漁場的預(yù)報中[11- 12]，它能夠把各種環(huán)境因子有機的整合在一起。為此，本研究利用我國東太平洋黃鰭金槍魚延繩釣生產(chǎn)統(tǒng)計數(shù)據(jù)，對黃鰭金槍魚漁場的空間分布特征進行分析，探討海洋環(huán)境因子與黃鰭金槍魚產(chǎn)量和CPUE的關(guān)系，并建立黃鰭金槍魚棲息地適應(yīng)性模型，為東太平洋黃鰭金槍魚中心漁場預(yù)報提供參考。

1材料與方法

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

生產(chǎn)統(tǒng)計資料來自上海金優(yōu)遠洋漁業(yè)有限公司2011—2012年在東太平洋公海金槍魚延繩釣生產(chǎn)數(shù)據(jù)。其中2011年的數(shù)據(jù)用于黃鰭金槍魚棲息地研究，2012年的數(shù)據(jù)用于棲息地指數(shù)模型驗證。共6艘生產(chǎn)船，每艘船噸位均為157 t，主機功率均為407kW，冷海水保鮮。作業(yè)海域為東太平洋海域(20°N—35°S、85°W—155°W)。數(shù)據(jù)包括作業(yè)時間、作業(yè)位置、黃鰭金槍魚產(chǎn)量(kg)和鉤數(shù)。

海表溫(SST)數(shù)據(jù)來自美國國家航空航天局的遙感資料(http://poet.jpl.nasa.gov/)，其數(shù)據(jù)空間分辨率為1°×1°。海面高度(SSH，m)的數(shù)據(jù)來自哥倫比亞大學衛(wèi)星遙感網(wǎng)站(http://iridl.ldeo.columbia.edu/docfind/databrief/cat-ocean.html)，空間分辨率為 1°×1°。時間和區(qū)域與生產(chǎn)統(tǒng)計數(shù)據(jù)相同。

1.2漁場分布與海表溫、海面高度的關(guān)系分析

首先以經(jīng)緯度1°×1°為空間統(tǒng)計單位，按月對其作業(yè)位置、產(chǎn)量和放鉤數(shù)進行初步統(tǒng)計，并計算平均每千鉤產(chǎn)量(CPUE，kg/千鉤)。利用頻度分析法按SST 1 ℃為組距、SSH 0.1 m為組距，分別分析各月產(chǎn)量、CPUE與SST和SSH的關(guān)系，獲得各月作業(yè)漁場SST和SSH的最適范圍。

1.3棲息地模型的建立

通常在作業(yè)漁船下鉤之前，船長會根據(jù)探漁儀映像、海洋環(huán)境狀況、周圍漁船作業(yè)情況進行綜合判斷，使得作業(yè)漁船往往會集中在某一區(qū)域生產(chǎn)，同時作業(yè)漁船之間會產(chǎn)生相互干擾，從而影響CPUE值[13-14]。因此，本研究以產(chǎn)量作為表征棲息地指數(shù)的指標。分1—3月、4—6月、7—9月和10—12月各季節(jié)分別建立基于SST、SSH的適應(yīng)性指數(shù)(SI)模型。本文假定最高產(chǎn)量(PROmax)為黃鰭金槍魚資源分布最多的海域，認定其棲息地指數(shù)為1；作業(yè)產(chǎn)量為0時，則認定是長鰭金槍魚資源分布較少的海域，認定其 HSI為0[15]。

(1)

式中，SIi為i季度的適應(yīng)性指數(shù)；PROi,max為i季度的最大產(chǎn)量。

利用一元非線性回歸分別建立SST與SI、SSH與SI 之間的關(guān)系模型，利用 DPS7.5 軟件求解。利用算術(shù)平均法(arithmetic mean model，AMM)[11,16]計算棲息地綜合指數(shù)(HSI)，HSI 在 0(不適宜)到 1(最適宜)之間變化。算術(shù)平均值法在金槍魚棲息地指數(shù)模型中得到較好應(yīng)用[12,17-18]。計算公式為：

HSI= (SISST+SISSH)/2

(2)

式中，SISST、SISSH分別為基于SST、SSH的適應(yīng)性指數(shù)。

1.4中心漁場預(yù)報驗證

(1)將 2012 年生產(chǎn)統(tǒng)計數(shù)據(jù)(PRO)和棲息地指數(shù)(HSI)均分為 5 個級別[19]。將產(chǎn)量統(tǒng)計數(shù)據(jù)(PRO)采用自然邊界法(Natural Breaks)[20]進行劃分，0 t≤PRO<500 t，記為等級1；500 t≤PRO<1000 t，記為等級2；1000 t≤PRO<5000 t，記為等級3；5000 t≤PRO<10000 t，記為等級4；PRO>10000 t，記為等級5。同樣，HSI也劃分為 5 個等級，即：0.0≤HSI <0.1，記為等級1；0.1≤HSI <0.3，記為等級2；0.3≤HSI <0.5，記為等級3；0.5≤HSI <0.7，記為等級4；0.7≤HSI≤1.0，記為等級5。

(2)對于同一個作業(yè)漁區(qū)(1°×1°)，如果其產(chǎn)量數(shù)據(jù)級別與棲息地指數(shù)級別相同或相差之絕對值小于2，則認為模型預(yù)報正確；如果級別相差之絕對值大于等于2，則認為模型預(yù)報不正確[19]。

2結(jié)果

2.1產(chǎn)量和 CPUE 的逐月分布

分析認為，高產(chǎn)(月產(chǎn)量超過 200 t)分布在2—3月和6月，產(chǎn)量最高為 2月，達到250 t 以上(圖1)，占全年總產(chǎn)量的 12.6%，其 CPUE 為 79.5 kg/千鉤。產(chǎn)量最低的為 4 月，僅為 77.8 t，占全年總產(chǎn)量的 3.9%，其 CPUE 為 31.2 kg/千鉤(圖1)。

圖1　黃鰭金槍魚延繩釣產(chǎn)量和單位捕撈努力量漁獲量(CPUE)月變化Fig.1　Monthly catch and catch per fishing unit (CPUE) of yellowfin tuna

2.2產(chǎn)量、CPUE 與SST 的關(guān)系

由圖2可知，1—3 月產(chǎn)量和 CPUE 較高的SST 為 24—27 ℃；4—6月產(chǎn)量和CPUE 較高的SST為 25—29 ℃；7—9月份產(chǎn)量與CPUE較高的SST基本為26—28 ℃，此外發(fā)現(xiàn)在SST為24 ℃以下海域，其產(chǎn)量和CPUE分布有些不一致。10—12月產(chǎn)量和CPUE較高的SST為25—28 ℃。

圖2　黃鰭金槍魚延繩釣月產(chǎn)量和 CPUE 與表溫關(guān)系Fig.2　Monthly catch and CPUE of Yellowfin tuna with the changes of sea surface temperature柱狀為產(chǎn)量；線段為CPUE；SST:表溫Sea surface temperature

2.3產(chǎn)量、CPUE 與SSH的關(guān)系

由圖3可知，1—3 月產(chǎn)量最高時的SSH 為 0.3 m，1 月和 3月CPUE 最高時的SSH為 0.3 m，2月CPUE 最高時 SSH為 0.2 m； 4—5月產(chǎn)量和 CPUE 較高的 SSH 為 0.3—0.8 m；6月產(chǎn)量和 CPUE最高的 SSH 均為 0.6 m；7月和8月產(chǎn)量最高時的SSH為0.6m，7月和8月CPUE最高的SSH為0.5 m；9月和10月產(chǎn)量最高時的SSH為0.4 m，但其CPUE最高的SSH分別為0.3 m和0.5 m；11月和12月產(chǎn)量和CPUE均最高的SSH分別為0.6 m和0.4 m。

圖3　黃鰭金槍魚月產(chǎn)量和 CPUE 與 海面高度的關(guān)系Fig.3　Relationship between monthly catch and CPUE of yellowfin tuna and sea surface height柱狀為產(chǎn)量；線段為CPUE；SSH:海面高度 sea surface height

2.4HSI 模型的建立

利用一元非線性回歸模型擬合以 SST和SSH 為基礎(chǔ)的 SI 曲線。擬合的 SI 曲線模型見表 1，在統(tǒng)計上均為顯著(P<0.05)。

表1　基于單因子的各月黃鰭金槍魚棲息地指數(shù)模型

SSH:海面高度sea surface height；SST:表溫sea surface temperature

2.5模型驗證分析

表2可反映出各月份HSI為0.6以上時的作業(yè)點個數(shù)，其各月占全月作業(yè)點數(shù)的比重在60%—75%。HSI為0.6以上時，全年作業(yè)點個數(shù)為265個，占總作業(yè)點數(shù)的66.42%，因此模型能夠較好的反映東太平洋黃鰭金槍魚漁場的分布情況。

表2　各月份不同棲息地指數(shù)(HSI)值下的實際作業(yè)點個數(shù)

HSI:棲息地指數(shù)habitat suitability index

2.6漁情預(yù)報驗證結(jié)果

由表 3可知，1—3 月份作業(yè)漁區(qū)數(shù)為35個，利用HSI模型預(yù)報正確的為21個，預(yù)報準確率為60%；4—6月份作業(yè)漁區(qū)數(shù)為24個，預(yù)報準確率為62%；7—9月份作業(yè)漁區(qū)數(shù)為28個，預(yù)報準確率為68%；10—12月份作業(yè)漁區(qū)41個，預(yù)報準確率為71%。1—12 月份平均預(yù)報準確率為66%。

表3　漁場預(yù)報結(jié)果統(tǒng)計

3討論與分析

(1)溫度是影響海洋魚類活動最重要的環(huán)境因子之一，直接或間接地影響到魚類資源量的分布、洄游和空間集群等[21]，本文著重對東太平洋黃鰭金槍魚漁場分布及其與SST關(guān)系進行了分析，并建立了基于SST的SI模型；通過產(chǎn)量、CPUE分析了其作業(yè)漁場的季節(jié)變化，得出了黃鰭金槍魚漁場空間分布的一些初步規(guī)律，作業(yè)漁場多分布在SST為24—29 ℃的海域，約占總漁獲量的96%以上，上述SST范圍可作為全年各月中心漁場分布的指標之一。盡管所分析的數(shù)據(jù)來源和空間尺度不一，但本研究結(jié)果與崔雪森等[22]的研究結(jié)論基本一致。研究結(jié)果可為漁業(yè)生產(chǎn)提供參考。

(2)SSH雖然沒有像水溫那樣顯著地影響著魚類活動，但是SSH也是影響魚類洄游、集群和分布的關(guān)鍵性因素之一[13]。本文分析了黃鰭金槍魚漁場分布及其與SSH關(guān)系，并建立基于SSH的SI模型，作業(yè)漁場多分布在SSH為0.3—0.7 m的海域，約占總漁獲量的95%以上，上述SSH范圍可作為全年各月中心漁場分布的指標之一。

(3)本研究根據(jù)已建立的棲息地指數(shù)模型，利用 SST和SSH兩個海洋環(huán)境因子對東太平洋黃鰭金槍魚漁場進行了預(yù)測。根據(jù) 2012年各月實際產(chǎn)量分布與理論計算獲得 HSI 分析，其平均漁場預(yù)報精度達到了66%。在所有月份中，實際作業(yè)漁場的范圍基本上落在漁情預(yù)報的理論范圍內(nèi)(表3)。因此，本研究所建立的漁情預(yù)報模型用來預(yù)測東太平洋黃鰭金槍魚漁場是可行的。當然，漁情預(yù)報的精度和檢驗方法還有進一步改進的地方，比如在模型構(gòu)建應(yīng)考慮更多的環(huán)境因子，并且要不斷更新生產(chǎn)統(tǒng)計數(shù)據(jù)來完善漁情預(yù)報模型。研究結(jié)果對指導(dǎo)漁業(yè)生產(chǎn)節(jié)省生產(chǎn)成本有實際意義。

(4)通過分析CPUE和產(chǎn)量與SST、SSH的關(guān)系，發(fā)現(xiàn)最適的SST和SSH的范圍并不完全一致。因為在商業(yè)化作業(yè)海域，船長一方面通常是依靠探魚儀來尋找中心漁場，同時往往根據(jù)船隊的作業(yè)情況來尋找漁場，因此往往會集中在一個海域進行集中作業(yè)。因此，在作業(yè)集中的海域產(chǎn)量往往會比較高，同時也影響到CPUE值；相反，一些作業(yè)漁船少的區(qū)域，其CPUE有可能較高。因此，在今后的研究中需要針對不同作業(yè)類型和不同漁業(yè)，研究采用CPUE和產(chǎn)量哪一個更為合適。

盡管上述建立的東太平洋黃鰭金槍魚棲息地指數(shù)模型有較高的精確度，但是黃鰭金槍魚漁場還與500 m以內(nèi)水層的水溫結(jié)構(gòu)、葉綠素濃度、溫躍層、餌料生物分布、鋒面以及 厄爾尼諾-南方濤動(ENSO)等有關(guān)[2- 4,8- 9]，這些因素對黃鰭金槍魚漁場均有一定的影響。因此，在以后的研究中需綜合考慮上述因子，以進一步完善黃鰭金槍魚棲息地模型，為海洋漁業(yè)生產(chǎn)提供科學參考。

參考文獻(References):

[1]Avise J C, Lansman R A. Polymorphism of mitochondrial DNA in populations of higher animals // Nei M, Koehn R K. Evolution of Genes and Proteins. Sunderland, Mass:Sinauer Associates Inc., 1983:147- 164.

[2]Holland K N, Brill R W, Chang R K C. Horizontal and vertical movements of yellowfin and bigeye tuna associated with fish aggregating devices. Fisheries Bulletin, 1990, 88(3):493- 507.

[3]Brill R W, Block B A, Boggs C H, Bigelow K A, Freund E V, Marcinek D J. Horizontal movements and depth distribution of large adult yellowfin tuna (Thunnusalbacares) near the Hawaiian Islands, recorded using ultrasonic telemetry:implications for the physiological ecology of pelagic fishes. Marine Biology, 1999, 133(3):395- 408.

[4]Dagorn L, Holland K N, Hallier J P, Taquet M, Moreno G, Sancho G, Itano D G, Aumeeruddy R, Girard C, Million J, Fonteneau. Deep diving behavior observed in yellowfin tuna (Thunnusalbacares). Aquatic Living Resources, 2006, 19(1):85- 88.

[5]Lee P F, Chen I C, Tseng W N. Distribution patterns of three dominant tuna species in the indian ocean [EB/OL]. [2011-01-28]. http://proceedings.esri.com/library/userconf/ proc99/proceed/papers/pap564/p564.htm.

[6]Mohri M, Nishida T. Consideration on distribution of adult yellowfin tuna (Thunnusalbacares) in the Indian Ocean based on Japanese tuna longline fisheries and survey information. IOTC Proceedings, 2000, 3:276- 282.

[7]Marsac F. Changes in depth of yellowfin tuna habitat in the Indian Ocean:an historical perspective 1955-2001. IOTC Proceedings, 2002, 5:450- 458.

[8]Romena N A. Factors affecting distribution of adult yellowfin tuna (Thunnusalbacares) and its reproductive ecology in the Indian Ocean based on Japanese tuna longline fisheries and survey information. IOTC Proceedings, 2001, 4:336- 389.

[9]Song L M, Zhang Y, Xu L X, Jiang W X, Wang J Q. Environmental preferences of longlining for yellowfin tuna (Thunnusalbacares) in the tropical high seas of the Indian Ocean. Fisheries Oceanography, 2008, 17(4):239- 253.

[10]Maury O, Gascuel D, Marsac F, Fonteneau A, Rosa A L D. Hierarchical interpretation of nonlinear relationships linking yellowfin tuna (Thunnusalbacares) distribution to the environment in the Atlantic Ocean. Canadian Journal of Fisheries and Aquatic Sciences, 2001, 58(3):458- 469.

[11]Brown S K, Buja K R, Jury S H, Monaco M E. Habitat suitability index models for eight fish and invertebrate species in Casco and Sheepscot Bays, Maine. North American Journal of Fisheries Management, 2000, 20(2):408- 435.

[12]馮波, 陳新軍, 許柳雄. 應(yīng)用棲息地指數(shù)對印度洋大眼金槍魚分布模式的研究. 水產(chǎn)學報, 2007, 31(6):805- 812.

[13]陳新軍. 漁業(yè)資源與漁場學. 北京:海洋出版社, 2004.

[14]陳新軍. 漁業(yè)資源經(jīng)濟學. 北京:中國農(nóng)業(yè)出版社, 2004.

[15]Mohri M, Nishida T. Seasonal changes in bigeye tuna fishing areas in relation to the oceanographic parameters in the Indian Ocean. Journal of National Fisheries University, 1999, 47(2):43- 54.

[16]Hess G R, Bay J M. A regional assessment of windbreak habitat suitability. Environmental Monitoring and Assessment, 2000, 61(2):239- 256.

[17]Wang C H, Wang S B. Assessment of South Pacific albacore stock (Thunnusalalunga) by improved Schaefer model. Journal of Ocean University of China, 2006, 5(2):106- 114.

[18]Zainuddin M, Saitoh K, Saitoh S. Albacore (Thunnusalalunga) fishing ground in relation to oceanographic conditions in the western North Pacific Ocean using remotely sensed satellite data. Fisheries Oceanography, 2008, 17(2):61- 73.

[19]高峰, 陳新軍, 范江濤, 雷林, 官文江. 西南大西洋阿根廷滑柔魚中心漁場預(yù)報的實現(xiàn)及驗證. 上海海洋大學學報, 2011, 20(5):754- 758.

[20]Smith R M. Comparing traditional methods for selecting class intervals on choropleth maps. The Professional Geographer, 1986, 38(1):62- 67.

[22]崔雪森, 樊偉, 張晶. 太平洋黃鰭金槍魚延繩釣漁獲分布及漁場水溫淺析. 海洋通報, 2005, 24(5):54- 59.

Forecasting fishing ground ofYellowfintunaintheEasternPacificOceanbasedonthehabitatsuitabilityindex

ZHAOHailong1,4,CHENXinjun1,2,3,4,*,FANGXueyan1,4

1 College of Marine Sciences of Shanghai Ocean University, Shanghai 201306, China 2 National Engineering Research Center for Oceanic Fisheries, Shanghai Ocean University, Shanghai 201306, China 3 The Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China

4CollaborativeInnovationCenterforDistant-waterFisheries,Shanghai201306,China

Abstract:Yellowfin tuna, Thunnus albacares, is a species of tuna found in the eastern Pacific Ocean, and has been one of the main targets of purse seine and longline fisheries since the 1970s. Chinese longline tuna fleets started to capture this species in 1999. According to the catch statistics from the Food and Agriculture Organization of the United Nations (FAO), during the period between 2008 and 2012, the annual catch of T. albacares ranged from 185 000 t to 260 000 t with an average of 213 000 t. The formation of T. albacares fishing grounds is considered to be complicated, and may be affected by many environmental factors. Many methods have been used to predict the location of fishing grounds. These include habitat suitability index (HSI), which is generally used to describe the quality of fish habitat, but recently has been applied to predict the location of fishing grounds. HSI models can be also used to inform fishery management and fish conservation. In this study, a HSI model was developed to predict the locations of T. albacares fishing grounds in the eastern Pacific using sea surface temperature (SST) and sea surface height (SSH) as explanatory environmental variables. The suitable ranges of SST and SSH were estimated using the frequency analysis method. Catch data for T. albacares were used as a suitability index, and quarterly suitability curves based on SST and SSH were derived using non-linear regression. Catch data were obtained from the Chinese longline fishery operating in the eastern Pacific Ocean (20°N—30°S and 85°W—155°W) in 2011. The spatial resolution of catch and environmental data is 1 degree latitude by 1 degree longitude, and data were recorded monthly. The HSI model was set up using an arithmetic mean model (AMM). The model was validated using separate a data set (2012 catch data from the same area). The T. albacares fishing grounds were mainly located in the waters with SST 24—29°C and SSH 0.3—0.7 m. The suitability index model for each factor (SST and SSH) was significant (P < 0.05). The accuracy with which fishing grounds were predicted for each quarter in 2012 varied from 60% to 71%, with an average of 66%. The actual T. albacares fishing grounds were almost all located in the forecast areas in all months. The HIS forecasting model developed in this study could provide valuable information for finding T. albacares fishing grounds in the eastern Pacific Ocean. However, in order to better forecast T. albacares fishing grounds using HIS models, more environmental factors should be included in the model, such as water temperature structure, sea front, Chlorophyll-a, and El Nio-Southern Oscillation (ENSO) index.

Key Words:Thunnus albacares; fishing ground forecasting; eastern pacific; habitat suitability index; environmental data

DOI:10.5846/stxb201405130975

*通訊作者

Corresponding author.E-mail:xjchen@shou.edu.cn

收稿日期:2014- 05- 13; 網(wǎng)絡(luò)出版日期:2015- 06- 11

基金項目:國家863計劃(2012AA092303); 國家發(fā)改委產(chǎn)業(yè)化專項(2159999); 上海市科技創(chuàng)新行動計劃(12231203900)

趙海龍, 陳新軍, 方學燕.基于棲息地指數(shù)的東太平洋黃鰭金槍魚漁場預(yù)報.生態(tài)學報,2016,36(3):778- 785.

Zhao H L, Chen X J, Fang X Y.Forecasting fishing ground of Yellowfin tuna in the Eastern Pacific Ocean based on the habitat suitability index.Acta Ecologica Sinica,2016,36(3):778- 785.