• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于棲息地指數(shù)的東太平洋黃鰭金槍魚(yú)漁場(chǎng)預(yù)報(bào)

    2016-04-13 00:51:06趙海龍陳新軍方學(xué)燕
    生態(tài)學(xué)報(bào) 2016年3期

    趙海龍, 陳新軍, 方學(xué)燕

    1上海海洋大學(xué)海洋科學(xué)學(xué)院, 上海 201306

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

    3 上海海洋大學(xué)大洋漁業(yè)可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室, 上?!?01306

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

    ?

    基于棲息地指數(shù)的東太平洋黃鰭金槍魚(yú)漁場(chǎng)預(yù)報(bào)

    趙海龍1,4, 陳新軍1,2,3,4,*, 方學(xué)燕1,4

    1上海海洋大學(xué)海洋科學(xué)學(xué)院, 上海201306

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

    3 上海海洋大學(xué)大洋漁業(yè)可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室, 上海201306

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

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

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

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

    1材料與方法

    1.1數(shù)據(jù)來(lái)源

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

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

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

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

    1.3棲息地模型的建立

    通常在作業(yè)漁船下鉤之前,船長(zhǎng)會(huì)根據(jù)探漁儀映像、海洋環(huán)境狀況、周?chē)鷿O船作業(yè)情況進(jìn)行綜合判斷,使得作業(yè)漁船往往會(huì)集中在某一區(qū)域生產(chǎn),同時(shí)作業(yè)漁船之間會(huì)產(chǎn)生相互干擾,從而影響CPUE值[13-14]。因此,本研究以產(chǎn)量作為表征棲息地指數(shù)的指標(biāo)。分1—3月、4—6月、7—9月和10—12月各季節(jié)分別建立基于SST、SSH的適應(yīng)性指數(shù)(SI)模型。本文假定最高產(chǎn)量(PROmax)為黃鰭金槍魚(yú)資源分布最多的海域,認(rèn)定其棲息地指數(shù)為1;作業(yè)產(chǎn)量為0時(shí),則認(rèn)定是長(zhǎng)鰭金槍魚(yú)資源分布較少的海域,認(rèn)定其 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]計(jì)算棲息地綜合指數(shù)(HSI),HSI 在 0(不適宜)到 1(最適宜)之間變化。算術(shù)平均值法在金槍魚(yú)棲息地指數(shù)模型中得到較好應(yīng)用[12,17-18]。計(jì)算公式為:

    HSI= (SISST+SISSH)/2

    (2)

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

    1.4中心漁場(chǎng)預(yù)報(bào)驗(yàn)證

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

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

    2結(jié)果

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

    分析認(rèn)為,高產(chǎn)(月產(chǎn)量超過(guò) 200 t)分布在2—3月和6月,產(chǎn)量最高為 2月,達(dá)到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 黃鰭金槍魚(yú)延繩釣產(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 黃鰭金槍魚(yú)延繩釣月產(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)量最高時(shí)的SSH 為 0.3 m,1 月和 3月CPUE 最高時(shí)的SSH為 0.3 m,2月CPUE 最高時(shí) 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)量最高時(shí)的SSH為0.6m,7月和8月CPUE最高的SSH為0.5 m;9月和10月產(chǎn)量最高時(shí)的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 黃鰭金槍魚(yú)月產(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 曲線模型見(jiàn)表 1,在統(tǒng)計(jì)上均為顯著(P<0.05)。

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

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

    2.5模型驗(yàn)證分析

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

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

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

    2.6漁情預(yù)報(bào)驗(yàn)證結(jié)果

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

    表3 漁場(chǎng)預(yù)報(bào)結(jié)果統(tǒng)計(jì)

    3討論與分析

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

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

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

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

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

    參考文獻(xià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ù)對(duì)印度洋大眼金槍魚(yú)分布模式的研究. 水產(chǎn)學(xué)報(bào), 2007, 31(6):805- 812.

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

    [14]陳新軍. 漁業(yè)資源經(jīng)濟(jì)學(xué). 北京:中國(guó)農(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ú)中心漁場(chǎng)預(yù)報(bào)的實(shí)現(xiàn)及驗(yàn)證. 上海海洋大學(xué)學(xué)報(bào), 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]崔雪森, 樊偉, 張晶. 太平洋黃鰭金槍魚(yú)延繩釣漁獲分布及漁場(chǎng)水溫淺析. 海洋通報(bào), 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

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

    趙海龍, 陳新軍, 方學(xué)燕.基于棲息地指數(shù)的東太平洋黃鰭金槍魚(yú)漁場(chǎng)預(yù)報(bào).生態(tài)學(xué)報(bào),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.

    日本一本二区三区精品| 麻豆一二三区av精品| 亚洲人成网站高清观看| 国语自产精品视频在线第100页| 又粗又爽又猛毛片免费看| 国产私拍福利视频在线观看| 校园人妻丝袜中文字幕| 久久草成人影院| 国产激情偷乱视频一区二区| 亚洲av中文字字幕乱码综合| 亚洲av中文字字幕乱码综合| 99在线视频只有这里精品首页| 18禁在线播放成人免费| 性色avwww在线观看| 午夜爱爱视频在线播放| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久久久久| 国产在视频线在精品| 亚洲人成网站在线播| 久久久久久久久久久丰满| 丝袜喷水一区| 天堂网av新在线| 免费人成在线观看视频色| 日韩欧美国产在线观看| 久久久久久久亚洲中文字幕| 久久久欧美国产精品| 免费观看在线日韩| 久久精品国产自在天天线| 国产一区二区亚洲精品在线观看| 六月丁香七月| 男插女下体视频免费在线播放| 男人的好看免费观看在线视频| 青春草亚洲视频在线观看| 乱人视频在线观看| 淫秽高清视频在线观看| or卡值多少钱| 日本av手机在线免费观看| 搡老妇女老女人老熟妇| 综合色丁香网| 亚洲av第一区精品v没综合| 观看美女的网站| 综合色丁香网| 亚洲av熟女| 欧美zozozo另类| 欧美激情在线99| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色| 午夜免费男女啪啪视频观看| 99久久久亚洲精品蜜臀av| 26uuu在线亚洲综合色| 欧美性猛交╳xxx乱大交人| 99在线视频只有这里精品首页| 尤物成人国产欧美一区二区三区| 欧美日韩综合久久久久久| 神马国产精品三级电影在线观看| 色哟哟·www| 插阴视频在线观看视频| 午夜精品一区二区三区免费看| 成人亚洲欧美一区二区av| 亚洲乱码一区二区免费版| 少妇的逼好多水| 亚洲av第一区精品v没综合| 久久99热6这里只有精品| a级一级毛片免费在线观看| 亚洲电影在线观看av| www.色视频.com| 一夜夜www| 成年免费大片在线观看| 日韩国内少妇激情av| 亚洲真实伦在线观看| 欧美一区二区国产精品久久精品| 秋霞在线观看毛片| 麻豆成人av视频| 99国产极品粉嫩在线观看| 国产成人福利小说| 成人av在线播放网站| 午夜精品一区二区三区免费看| 日本黄色视频三级网站网址| 久久人人爽人人爽人人片va| 国产黄色视频一区二区在线观看 | 人人妻人人看人人澡| 91久久精品电影网| 床上黄色一级片| 日韩亚洲欧美综合| 久久这里有精品视频免费| 我的老师免费观看完整版| 少妇的逼好多水| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验 | 特大巨黑吊av在线直播| 特级一级黄色大片| av国产免费在线观看| 欧美bdsm另类| 午夜福利在线观看吧| 日韩亚洲欧美综合| 中文字幕久久专区| 日韩精品青青久久久久久| 色综合亚洲欧美另类图片| 亚洲无线在线观看| 亚洲国产精品成人综合色| 一边摸一边抽搐一进一小说| 在现免费观看毛片| 一级av片app| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 啦啦啦啦在线视频资源| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 美女内射精品一级片tv| 精品久久久久久久久久免费视频| 国产三级中文精品| 欧美zozozo另类| 成年av动漫网址| .国产精品久久| 99热网站在线观看| 1024手机看黄色片| 尾随美女入室| 最新中文字幕久久久久| 国产成人影院久久av| 亚洲天堂国产精品一区在线| 国产精品一区二区性色av| 国产av不卡久久| 1024手机看黄色片| 搞女人的毛片| 99久久精品热视频| 狂野欧美白嫩少妇大欣赏| 波多野结衣巨乳人妻| 国产亚洲5aaaaa淫片| 免费大片18禁| 最近最新中文字幕大全电影3| 女同久久另类99精品国产91| 色综合站精品国产| 欧美极品一区二区三区四区| 天堂中文最新版在线下载 | 成人亚洲精品av一区二区| 人人妻人人澡人人爽人人夜夜 | 成人亚洲欧美一区二区av| 国产一区亚洲一区在线观看| 啦啦啦啦在线视频资源| 在线天堂最新版资源| 国产亚洲91精品色在线| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 麻豆乱淫一区二区| .国产精品久久| 在线免费观看的www视频| 久久久久网色| 亚洲欧美日韩东京热| 国产不卡一卡二| 欧美xxxx黑人xx丫x性爽| 国产成人a区在线观看| 国产乱人视频| 免费观看人在逋| 午夜精品国产一区二区电影 | 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久成人av| 中文字幕av成人在线电影| 搡老妇女老女人老熟妇| 亚洲在线自拍视频| 深夜a级毛片| videossex国产| 亚洲av成人精品一区久久| 久久精品人妻少妇| 丝袜美腿在线中文| 啦啦啦韩国在线观看视频| 国产极品天堂在线| 日日干狠狠操夜夜爽| 国产综合懂色| 人人妻人人澡欧美一区二区| 春色校园在线视频观看| 91精品一卡2卡3卡4卡| 欧美bdsm另类| 插阴视频在线观看视频| 欧美激情国产日韩精品一区| 国产精品人妻久久久影院| 国产成人aa在线观看| 熟妇人妻久久中文字幕3abv| 亚洲欧美成人综合另类久久久 | 日韩视频在线欧美| 日韩精品有码人妻一区| 男人的好看免费观看在线视频| 赤兔流量卡办理| 精品人妻视频免费看| 99久久无色码亚洲精品果冻| 亚洲av中文字字幕乱码综合| 亚洲真实伦在线观看| 色综合站精品国产| 1024手机看黄色片| or卡值多少钱| 精品国产三级普通话版| 日韩欧美精品v在线| 美女黄网站色视频| 亚洲四区av| 国内久久婷婷六月综合欲色啪| 在线播放无遮挡| 人人妻人人澡人人爽人人夜夜 | 天堂√8在线中文| 两个人视频免费观看高清| 最后的刺客免费高清国语| 91精品国产九色| 亚洲在线自拍视频| 久久久久久九九精品二区国产| 国产av不卡久久| 51国产日韩欧美| 国产成人一区二区在线| 大型黄色视频在线免费观看| 国产精品人妻久久久影院| 国产视频首页在线观看| 中文字幕免费在线视频6| 日韩大尺度精品在线看网址| 一进一出抽搐gif免费好疼| 国产成年人精品一区二区| 亚洲欧美日韩高清专用| 国产成人精品一,二区 | 三级毛片av免费| 亚洲精品乱码久久久v下载方式| 99久久九九国产精品国产免费| 亚洲第一区二区三区不卡| 国产男人的电影天堂91| 日本av手机在线免费观看| 夜夜夜夜夜久久久久| 日韩在线高清观看一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 1024手机看黄色片| 天堂网av新在线| 蜜臀久久99精品久久宅男| 丰满乱子伦码专区| 久久午夜福利片| 久久久久久久久久久丰满| 日韩国内少妇激情av| 99久久久亚洲精品蜜臀av| 亚洲精品色激情综合| 亚洲va在线va天堂va国产| 九九爱精品视频在线观看| 看十八女毛片水多多多| 五月伊人婷婷丁香| 午夜亚洲福利在线播放| 久久精品久久久久久久性| 色5月婷婷丁香| 2021天堂中文幕一二区在线观| 国国产精品蜜臀av免费| 国产美女午夜福利| 啦啦啦韩国在线观看视频| 久久久久久九九精品二区国产| 变态另类丝袜制服| 最近2019中文字幕mv第一页| 亚洲成av人片在线播放无| 国产精品综合久久久久久久免费| 久久精品久久久久久噜噜老黄 | 久久久久久久亚洲中文字幕| 国产精品av视频在线免费观看| 能在线免费观看的黄片| 久久久久国产网址| 一级二级三级毛片免费看| 亚洲人成网站高清观看| 婷婷六月久久综合丁香| 91精品国产九色| 国产高潮美女av| 岛国在线免费视频观看| 91精品一卡2卡3卡4卡| 亚州av有码| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 一级毛片我不卡| 久久久久久久久中文| 男人舔女人下体高潮全视频| 精品一区二区三区视频在线| 人体艺术视频欧美日本| 69人妻影院| 看片在线看免费视频| 国产精品久久久久久精品电影| 麻豆乱淫一区二区| 欧美精品一区二区大全| 日日摸夜夜添夜夜爱| 欧美变态另类bdsm刘玥| 久久精品影院6| 蜜桃久久精品国产亚洲av| 日本五十路高清| 午夜福利视频1000在线观看| .国产精品久久| 99在线人妻在线中文字幕| 国产v大片淫在线免费观看| 国产免费一级a男人的天堂| 女人十人毛片免费观看3o分钟| 啦啦啦观看免费观看视频高清| 黄片wwwwww| 亚洲成人久久爱视频| av福利片在线观看| 卡戴珊不雅视频在线播放| 在线播放无遮挡| 长腿黑丝高跟| 亚洲人成网站高清观看| 99久久中文字幕三级久久日本| 99久国产av精品| 久久久午夜欧美精品| 99久国产av精品国产电影| 夫妻性生交免费视频一级片| 亚洲七黄色美女视频| 国产色婷婷99| 亚洲av免费在线观看| 免费黄网站久久成人精品| 日本五十路高清| 给我免费播放毛片高清在线观看| 黄色日韩在线| 日本av手机在线免费观看| 赤兔流量卡办理| 国产精品一区二区三区四区久久| 国产又黄又爽又无遮挡在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲激情五月婷婷啪啪| 国产成人影院久久av| 久久鲁丝午夜福利片| 狂野欧美激情性xxxx在线观看| 欧美一级a爱片免费观看看| 久久99热这里只有精品18| 国产单亲对白刺激| 只有这里有精品99| 一进一出抽搐动态| 悠悠久久av| 一本久久精品| av福利片在线观看| 性插视频无遮挡在线免费观看| 内射极品少妇av片p| 欧美性猛交黑人性爽| 十八禁国产超污无遮挡网站| 国产蜜桃级精品一区二区三区| 国产v大片淫在线免费观看| 蜜臀久久99精品久久宅男| 国产精品美女特级片免费视频播放器| 国产高潮美女av| 日韩亚洲欧美综合| 国产高潮美女av| 老司机影院成人| 国产午夜福利久久久久久| 性色avwww在线观看| 免费av不卡在线播放| 69av精品久久久久久| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 国产白丝娇喘喷水9色精品| 一区二区三区四区激情视频 | 97热精品久久久久久| 白带黄色成豆腐渣| 18禁在线无遮挡免费观看视频| 国产精品一区二区三区四区免费观看| 国产成人a∨麻豆精品| 老熟妇乱子伦视频在线观看| 欧美成人a在线观看| 国内精品一区二区在线观看| 91在线精品国自产拍蜜月| 国产探花极品一区二区| 丝袜美腿在线中文| 一级黄色大片毛片| 男插女下体视频免费在线播放| 国产高清不卡午夜福利| 久久久久九九精品影院| 亚洲国产高清在线一区二区三| 国产成人影院久久av| 美女脱内裤让男人舔精品视频 | 人妻制服诱惑在线中文字幕| 免费搜索国产男女视频| 在线观看66精品国产| 亚洲五月天丁香| .国产精品久久| 美女内射精品一级片tv| 夜夜看夜夜爽夜夜摸| 高清毛片免费观看视频网站| 色哟哟哟哟哟哟| 好男人在线观看高清免费视频| 亚洲中文字幕一区二区三区有码在线看| 99久久精品国产国产毛片| 国产成人精品久久久久久| 男插女下体视频免费在线播放| 国产老妇女一区| 99国产精品一区二区蜜桃av| 97超碰精品成人国产| 麻豆国产97在线/欧美| 亚洲国产精品久久男人天堂| 国产精品久久久久久久电影| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 久久久色成人| 国产精品一区二区在线观看99 | 亚洲av免费高清在线观看| 亚洲av成人精品一区久久| 成年女人看的毛片在线观看| 我的女老师完整版在线观看| 一区福利在线观看| 毛片女人毛片| 天堂网av新在线| 国产蜜桃级精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 麻豆国产97在线/欧美| 91精品国产九色| 此物有八面人人有两片| 亚洲av熟女| 日韩大尺度精品在线看网址| 亚洲欧美精品专区久久| 免费观看人在逋| 亚洲人成网站在线播| 国内精品美女久久久久久| 日韩强制内射视频| 床上黄色一级片| 久99久视频精品免费| 午夜爱爱视频在线播放| 精品日产1卡2卡| 色综合色国产| 欧美日韩精品成人综合77777| 99久久无色码亚洲精品果冻| 免费av观看视频| 欧美日本亚洲视频在线播放| 久久精品综合一区二区三区| 性插视频无遮挡在线免费观看| 女人被狂操c到高潮| 老熟妇乱子伦视频在线观看| 亚洲av男天堂| 日日撸夜夜添| 免费观看的影片在线观看| 欧美变态另类bdsm刘玥| eeuss影院久久| 久久久色成人| 日韩大尺度精品在线看网址| 男女下面进入的视频免费午夜| 麻豆成人午夜福利视频| 亚洲第一电影网av| 日韩一本色道免费dvd| 久久久久久国产a免费观看| 99热全是精品| 国产亚洲精品久久久久久毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产探花极品一区二区| 婷婷色综合大香蕉| 久久午夜福利片| 国产免费男女视频| 国产精品国产三级国产av玫瑰| 久久久色成人| 亚洲性久久影院| 亚洲精品色激情综合| 国产视频内射| 免费观看精品视频网站| 国产高潮美女av| 久久久久九九精品影院| 美女 人体艺术 gogo| 99精品在免费线老司机午夜| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 久久久国产成人免费| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 亚洲av男天堂| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 欧美人与善性xxx| 亚洲国产欧美在线一区| 色视频www国产| 亚洲激情五月婷婷啪啪| 欧美激情在线99| 啦啦啦观看免费观看视频高清| eeuss影院久久| 久久国内精品自在自线图片| 卡戴珊不雅视频在线播放| 成人一区二区视频在线观看| 国产女主播在线喷水免费视频网站 | 国产精品1区2区在线观看.| 国产蜜桃级精品一区二区三区| 久久精品夜色国产| 国产综合懂色| 日韩国内少妇激情av| 1000部很黄的大片| 欧美一区二区国产精品久久精品| 国产麻豆成人av免费视频| 久久久久久久久久久免费av| 亚洲av成人av| 免费在线观看成人毛片| 成年女人看的毛片在线观看| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 国产 一区 欧美 日韩| 久久精品人妻少妇| 亚洲欧美日韩卡通动漫| 亚洲三级黄色毛片| 亚洲成人中文字幕在线播放| 久久久久久大精品| 亚洲人与动物交配视频| 91狼人影院| 爱豆传媒免费全集在线观看| 男人狂女人下面高潮的视频| 欧美三级亚洲精品| 99在线人妻在线中文字幕| videossex国产| 少妇熟女aⅴ在线视频| 国产精品三级大全| 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 国产免费一级a男人的天堂| 午夜精品在线福利| 国产av一区在线观看免费| 日韩欧美 国产精品| 亚洲婷婷狠狠爱综合网| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 欧美最黄视频在线播放免费| 亚洲乱码一区二区免费版| 亚洲精品久久国产高清桃花| 又爽又黄a免费视频| 天天躁夜夜躁狠狠久久av| 校园人妻丝袜中文字幕| 亚洲成a人片在线一区二区| av福利片在线观看| 日韩三级伦理在线观看| 韩国av在线不卡| 亚洲成人av在线免费| 久久久久久国产a免费观看| 久久久精品大字幕| 此物有八面人人有两片| 在线观看午夜福利视频| 久久久久免费精品人妻一区二区| 午夜福利在线观看免费完整高清在 | 国产精品,欧美在线| 久久久久久久久久久丰满| 日韩av不卡免费在线播放| 六月丁香七月| 黄色一级大片看看| 国产伦一二天堂av在线观看| 成人国产麻豆网| 久久久国产成人免费| 在线a可以看的网站| 神马国产精品三级电影在线观看| 亚洲欧美中文字幕日韩二区| 两个人的视频大全免费| 18禁在线播放成人免费| 成人一区二区视频在线观看| 国产激情偷乱视频一区二区| 国产 一区精品| 国产人妻一区二区三区在| 白带黄色成豆腐渣| 亚洲国产欧美人成| 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 六月丁香七月| 成人欧美大片| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久久久成人av| 久久久午夜欧美精品| 天天躁日日操中文字幕| 国产午夜精品一二区理论片| 精品久久久久久久久av| 人妻久久中文字幕网| 1000部很黄的大片| 给我免费播放毛片高清在线观看| 不卡视频在线观看欧美| 成人午夜高清在线视频| 可以在线观看毛片的网站| 精品一区二区免费观看| 内地一区二区视频在线| av天堂在线播放| 美女高潮的动态| 国产极品天堂在线| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 久久午夜福利片| 波野结衣二区三区在线| 我的女老师完整版在线观看| 久久九九热精品免费| h日本视频在线播放| 欧美性感艳星| 亚洲无线观看免费| 精品久久久久久久久av| av专区在线播放| 亚洲一级一片aⅴ在线观看| 国产中年淑女户外野战色| 欧美区成人在线视频| www.色视频.com| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 欧美成人免费av一区二区三区| 最好的美女福利视频网| 国产69精品久久久久777片| 久久精品人妻少妇| 亚洲最大成人手机在线| 国产精品一区二区性色av| 性色avwww在线观看| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 人妻制服诱惑在线中文字幕| 一级二级三级毛片免费看| 大又大粗又爽又黄少妇毛片口| 国产成人91sexporn| 欧美成人精品欧美一级黄| 1000部很黄的大片| 一边摸一边抽搐一进一小说| 看十八女毛片水多多多| 久久久久九九精品影院| 久久久久久久久久久免费av| 国模一区二区三区四区视频| 伦理电影大哥的女人| 好男人视频免费观看在线| 老熟妇乱子伦视频在线观看| 亚洲高清免费不卡视频| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 狂野欧美激情性xxxx在线观看| 国产色婷婷99| 精品午夜福利在线看| 真实男女啪啪啪动态图| 亚洲自偷自拍三级| 最近最新中文字幕大全电影3| 爱豆传媒免费全集在线观看| 国产精品永久免费网站| 中文精品一卡2卡3卡4更新|