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

    中西太平洋金槍魚圍網(wǎng)漁場(chǎng)聚類及其原因分析

    2016-12-22 01:53:54胡啟偉陳新軍徐良琦余景
    海洋學(xué)報(bào) 2016年12期
    關(guān)鍵詞:漁場(chǎng)金槍魚時(shí)間尺度

    胡啟偉,陳新軍,2,3,4*,徐良琦,余景

    (1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海 201306;2.大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室,上海201306;3.國(guó)家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心,上海 201306;4.遠(yuǎn)洋漁業(yè)協(xié)同創(chuàng)新中心,上海 201306;5.上海水產(chǎn)集團(tuán)總公司,上海 200090;6.中國(guó)水產(chǎn)科學(xué)院 南海水產(chǎn)研究所,廣東 廣州 510300)

    ?

    中西太平洋金槍魚圍網(wǎng)漁場(chǎng)聚類及其原因分析

    胡啟偉1,陳新軍1,2,3,4*,徐良琦5,余景6

    (1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海 201306;2.大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室,上海201306;3.國(guó)家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心,上海 201306;4.遠(yuǎn)洋漁業(yè)協(xié)同創(chuàng)新中心,上海 201306;5.上海水產(chǎn)集團(tuán)總公司,上海 200090;6.中國(guó)水產(chǎn)科學(xué)院 南海水產(chǎn)研究所,廣東 廣州 510300)

    本文根據(jù)1995—2010年我國(guó)中西太平洋金槍魚圍網(wǎng)生產(chǎn)統(tǒng)計(jì)數(shù)據(jù),按年和月不同時(shí)間空間分辨率對(duì)金槍魚圍網(wǎng)漁場(chǎng)進(jìn)行聚類分析,劃定不同漁場(chǎng)類型;同時(shí)結(jié)合海表溫度(SST)及Nio3.4區(qū)指數(shù),探討其漁場(chǎng)類型形成的原因。研究表明,在月為時(shí)間尺度下氣候異常事件(El Nio 和La Nia事件)發(fā)生頻次與漁場(chǎng)聚類結(jié)果類別相關(guān)聯(lián),在1—12月發(fā)生氣候異常事件頻數(shù)分布可劃分為以下幾個(gè)階段:1—3月、4—6、7—9月、10—12月,與月時(shí)間尺度下金槍魚圍網(wǎng)中心漁場(chǎng)4種聚類結(jié)果的時(shí)間范圍具有一致性。通過(guò)漁場(chǎng)重心聚類結(jié)果和El Nio 和La Nia事件分類統(tǒng)計(jì)對(duì)比發(fā)現(xiàn),在年時(shí)間尺度下,氣候異常事件的類型與聚類結(jié)果相關(guān)聯(lián),聚類結(jié)果同一類別包含的年份發(fā)生的氣候異常事件具有一致性,即在同一類別下1995年、1997年為強(qiáng)El Nio 年;1998年、2007、年、2009年為正常年份;2010年為強(qiáng)La Nia年;1999年、2000年、2001年為L(zhǎng)a Nia年;1996年、2008年為L(zhǎng)a Nia年;2002年、2004年為El Nio 年。研究認(rèn)為,中心漁場(chǎng)的年際聚類變化與El Nio、La Nia事件的發(fā)生分布具有很強(qiáng)的相關(guān)性,因此可以利用El Nio、La Nia指標(biāo)來(lái)預(yù)測(cè)漁場(chǎng)的年間和月份間的變化。

    中西太平洋;金槍魚圍網(wǎng)漁場(chǎng);聚類分析;海表面溫度;厄爾尼諾;拉尼娜

    1 引言

    金槍魚漁業(yè)是眾多太平洋島國(guó)和地區(qū)重要的生活和經(jīng)濟(jì)來(lái)源[1],我國(guó)的金槍魚圍網(wǎng)主要分布在中西太平洋島國(guó)附近海域[2]。海洋環(huán)境對(duì)中西太平洋金槍魚圍網(wǎng)漁場(chǎng)有著重要影響[3],研究漁場(chǎng)類型對(duì)指導(dǎo)科學(xué)生產(chǎn)具有重要意義。全球氣候異常也會(huì)使其資源和漁場(chǎng)發(fā)生波動(dòng),特別是厄爾尼諾(El Nio)和拉尼娜(La Nia)事件[4]。有研究表明,厄爾尼諾-南方濤動(dòng)現(xiàn)象對(duì)中西太平洋鰹魚圍網(wǎng)漁場(chǎng)的空間分布有顯著影響[5]。但氣候變化對(duì)中西太平洋金槍魚漁業(yè)資源豐度卻有不同的意見(jiàn)[1,6],Lehodey等[6]研究認(rèn)為,氣候變化會(huì)減少該區(qū)域金槍魚漁業(yè)產(chǎn)量,而Matear等[1]認(rèn)為沒(méi)有太大影響,有關(guān)El Nio 和La Nia 事件對(duì)金槍魚漁場(chǎng)空間位置和資源豐度分布影響機(jī)理有待進(jìn)一步研究。因此,本研究利用長(zhǎng)時(shí)間序列的中西太平洋金槍魚圍網(wǎng)漁場(chǎng)數(shù)據(jù),結(jié)合海洋表面溫度(SST)和El Nio/La Nia事件,運(yùn)用聚類和空間疊加分析方法分析中西太平洋金槍魚漁場(chǎng)分布及其環(huán)境狀況,掌握其資源空間分布的變化規(guī)律,為科學(xué)指導(dǎo)我國(guó)金槍魚圍網(wǎng)生產(chǎn)提供科學(xué)依據(jù)。

    2 材料與方法

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

    采用1995—2010年1-12月中西太平洋鰹魚生產(chǎn)統(tǒng)計(jì)數(shù)據(jù)以及對(duì)應(yīng)的海表溫度(SST)環(huán)境數(shù)據(jù),時(shí)間分辨率為月,空間分辨率為經(jīng)緯度5°×5°。經(jīng)緯范圍為0°~15°S,125°~180°E,其中生產(chǎn)數(shù)據(jù)包括月份、作業(yè)次數(shù)以及鰹魚自由群產(chǎn)量(單位:t)和鰹魚流木群(單位:t),CPUE為每月的捕撈產(chǎn)量(單位:t/月)。SST數(shù)據(jù)來(lái)自美國(guó)NOAA的OceanWatch數(shù)據(jù)庫(kù)(http://oceanwatch.pifsc.noaa.gov/las/servlets/dataset)。

    2.2 研究方法

    (1)CPUE計(jì)算。CPUE為每月的單位作業(yè)網(wǎng)次捕撈產(chǎn)量。

    (2)不同時(shí)間尺度的漁場(chǎng)聚類。從年和月兩種尺度對(duì)以CPUE為指標(biāo)的金槍魚圍網(wǎng)漁場(chǎng)進(jìn)行聚類分析。第一種情況,以 12個(gè)月份CPUE為指標(biāo)對(duì)不同年份CPUE類別進(jìn)行系統(tǒng)聚類。第二種情況,按照16年每年12個(gè)月份CPUE為指標(biāo)對(duì)年份進(jìn)行系統(tǒng)聚類。

    (3)聚類方法。系統(tǒng)聚類將n個(gè)CPUE樣品各自作為一類,規(guī)定樣品之間的距離和類與類之間的空間距離,將距離近的兩類并為一個(gè)新類,計(jì)算新類與其他類之間的空間距離,重復(fù)進(jìn)行兩個(gè)空間距離最近的類合并,每次減少一類,直至所有的樣品合并成為一類。最后形成一個(gè)親疏關(guān)系圖譜(聚類樹(shù)圖形),通常從圖上能清晰的看出應(yīng)分為幾類及每類包含的變量。由于在系統(tǒng)聚類中類平均很好地利用了所有樣品之間的信息,在很多情況下它被認(rèn)為是一種比較好的系統(tǒng)聚類方法。因此,本研究用系統(tǒng)聚類中的類平均法聚類(average linkage method)進(jìn)行兩種尺度下的聚類分析[7-8]。

    定義類與類之間的聚類為兩類最近樣品間的距離[8],即:

    (1)

    若某一類GK與GL聚成新的一類,記為GM,類GM與任意已有類GJ之間的距離為[8]:

    (2)

    類與類之間的平均距離定義為CPUE樣品之間平方距離的平均值。GK與GL之間的平方距離為[8]:

    (3)

    類間平方距離的遞推公式為[14]:

    (4)

    (4)不同類別漁場(chǎng)分布與海洋環(huán)境關(guān)系分析。上述聚類劃分后,重點(diǎn)探討不同類型(本研究中可以取一個(gè)代表性年份或者月份)CPUE與海洋環(huán)境的關(guān)系。利用R軟件繪出每年各月CPUE空間分布熱圖(HeatMap),通過(guò)熱圖顏色的變化分布可以直觀的觀察到CPUE每年每月的變化趨勢(shì)。此外,色差也能反應(yīng)出“冷熱”月份和年份。

    (5)z分?jǐn)?shù)(z-score),也稱標(biāo)準(zhǔn)分?jǐn)?shù)(standard score)是一個(gè)數(shù)與平均數(shù)的差再除以標(biāo)準(zhǔn)差的過(guò)程。本研究用z分?jǐn)?shù)表示:“一個(gè)給的時(shí)間尺度下CPUE分?jǐn)?shù)距離平均數(shù)多少個(gè)標(biāo)準(zhǔn)差”在平均數(shù)之上的分?jǐn)?shù)會(huì)得到一個(gè)正的標(biāo)準(zhǔn)分?jǐn)?shù),在平均數(shù)之下的分?jǐn)?shù)會(huì)得到一個(gè)負(fù)的標(biāo)準(zhǔn)分?jǐn)?shù)。z分?jǐn)?shù)是一種可以看出某分?jǐn)?shù)在分布中相對(duì)位置的方法。z分?jǐn)?shù)能夠真實(shí)的反應(yīng)一個(gè)分?jǐn)?shù)距離平均數(shù)的相對(duì)標(biāo)準(zhǔn)距離[9]。

    3 結(jié)果

    3.1 CPUE的時(shí)間變化

    由圖1可知,在月的時(shí)間尺度下,1995-2010年1-12月CPUE平均值變化不大,穩(wěn)定在15 t/月左右,但每月CPUE離散變化表現(xiàn)出明顯不同,總體來(lái)看2-7月的離散度要小于其他月份。由圖2可知,在年的時(shí)間尺度下,1995—2010年均CPUE值變化較大。1997年為16年間的最小值,平均CPUE為10.3 t/網(wǎng)次;2008年為16年間最大值,平均為19.8 t/網(wǎng)次??傮w來(lái)看(圖1,圖2),16年間年均CPUE的變化幅度遠(yuǎn)比月平均CPUE大。

    圖1 1995—2010年各月單位捕撈努力量漁獲量分布Fig.1 Monthly CPUE distribution during 1995 to 2010

    圖2 1995—2010年每年12個(gè)月份單位捕撈努力量漁獲量分布Fig.2 Monthly CPUE distribution in each year from 1995 to 2010

    (1)基于月份的金槍魚圍網(wǎng)漁場(chǎng)聚類

    根據(jù)R軟件繪制月熱圖和聚類樹(shù)(圖3a),漁場(chǎng)聚類可分為4類:1—3月、4—7月、8—9月及10—12月結(jié)果(表1)。以月為時(shí)間尺度1—12月發(fā)生氣候異常事件(El Nio 和La Nia事件)頻數(shù)的分布(圖4和圖5),可劃分為以下幾個(gè)階段:1—3月、4—6月、7—9月、10—12月,與金槍魚圍網(wǎng)中心漁場(chǎng)聚類結(jié)果基本具有一致性。此外,由熱圖可知在月時(shí)間尺度下(圖3a),z分?jǐn)?shù)范圍在(-3,3)之間漁場(chǎng)中心變化較大,但分布比較集中在均值附近,只有少部分月份漁場(chǎng)中心變異較大。

    (2)基于年份的金槍魚圍網(wǎng)漁場(chǎng)聚類

    在年時(shí)間尺度下(圖3b),漁場(chǎng)聚類可分為6類,第一類:1995年、1997年,第二類:1998年、2007年、2009年,第三類:2005年、2010年,第四類:1999年、2000年、2001年,第五類:1996年、2006年、2008年,第六類:2002年、2003年、2004年結(jié)果(表1)。通過(guò)漁場(chǎng)中心聚類結(jié)果和El Nio 和La Nia事件(圖3)分類統(tǒng)計(jì)對(duì)比發(fā)現(xiàn)第一類:1995年、1997年為強(qiáng)厄爾尼諾年,第二類:1998年、2007年、2009年為正常年份,第三類:2010年為強(qiáng)拉尼娜年,第四類:1999年、2000年、2001年為拉尼娜年,第五類:1996年、2008年為拉尼娜年,第六類:2002年、2004年為厄爾尼諾年。因此,中心漁場(chǎng)的年際聚類變化與El Nio、La Nia事件的發(fā)生分布基本具有一致性。另外,由熱圖分析(圖3b)可知z分?jǐn)?shù)范圍在(-2,2)之間漁場(chǎng)中心仍存在變化,變化幅度小于月時(shí)間尺度,但分布較發(fā)散,大多數(shù)年份的漁場(chǎng)中心都偏離平均值。這可能是因?yàn)槟陼r(shí)間尺度大,漁場(chǎng)中心受到不同尺度環(huán)境要素的影響。

    另外,在聚類的結(jié)果中,通過(guò)其(圖3a,圖3b)下方的聚類色標(biāo)可發(fā)現(xiàn)在年時(shí)間尺度下1995、1997年;月時(shí)間尺度下8、9月漁場(chǎng)中心差異較大。

    圖3 不同時(shí)間尺度下單位捕撈努力量漁獲量聚類和熱圖分布Fig.3 The clustering and heatmap distribution of CPUE at different time scales

    時(shí)間尺度類別分類結(jié)果年11995年、1997年21998年、2007年、2009年32005年、2010年41999年、2000年、2001年51996年、2006年、2008年62002年、2003年、2004年月18月、9月210月、11月、12月31月、2月、3月44月、5月、6月、7月

    圖4 1995年1月-2010年12月Nio3.4 區(qū)SSTA 時(shí)間序列分布Fig.4 The time series distribution of SSTA in the area Nio3.4 from January of 1995 to December of 2010

    圖5 1995—2010年12個(gè)月份El Nio 和La Nia事件發(fā)生次數(shù)Fig.5 The frequency of El Nio and La Nia event of 12 months from 1995-2010

    3.3 漁場(chǎng)類型與SST的關(guān)系

    在月時(shí)間尺度下,由圖1和6a可知,CPUE變化基本和SST變化一致。其中SST在6、7、8、9月份較其他月份較低,波動(dòng)幅度偏高,同期CPUE較其他月份出現(xiàn)下降。在年時(shí)間尺度下,由圖2、圖6b可知,在1995-2002年SST年際變化幅度較其余年份大,同期CPUE波動(dòng)范圍較大;在2003-2010年SST變化幅度較小,同期CPUE相對(duì)穩(wěn)定。另外氣候異常變化對(duì)CPUE影響較低,例如1997年、2008年是都是氣候異常年份(圖4),同期年SST和年CPUE都是出現(xiàn)較大幅度的波動(dòng)。因此,在不同的時(shí)間尺度下都發(fā)現(xiàn)SST的變化與CPUE的變化密切相關(guān),氣候異常對(duì)CPUE的影響顯著。

    4 討論與分析

    水溫是影響魚類活動(dòng)最重要的環(huán)境因子之一,魚類的分布、洄游遷移和集群等會(huì)直接或間接地受到環(huán)境溫度的限制[10—11]。本文通過(guò)不同時(shí)間尺度下CPUE和SST時(shí)間序列分析(圖1,圖2,圖6)發(fā)現(xiàn),當(dāng)SST值在27.0~30.0℃時(shí),CPUE變化趨勢(shì)和SST變化趨勢(shì)基本一致。

    月份時(shí)間序列中5—9月SST值明顯低于其他月份分布在27.5~28.5℃之間(圖6b),對(duì)應(yīng)時(shí)間階段CPUE值也出現(xiàn)下降,其值分布在12.5~17.5 t 之間(圖1)。這反映出中西太平洋的金槍魚生長(zhǎng)有一定的適宜溫度,27.5~28.5℃并不適合它的生長(zhǎng),唐峰華等[3]研究表明中西太平洋的金槍魚生長(zhǎng)最適宜的海表溫度為29~31℃[3]。

    圖6 1995—2010年不同時(shí)間尺度下海表面溫度分布Fig.6 The distribution of SST at different time scales during January of 1995 to December of 2010

    圖7 2009年10月(a)、11月(b)、12月(c)中西太平洋金槍魚CPUE空間分布與SST關(guān)系Fig.7 CPUE distribution of tuna and relationship with SST in the west-central Pacificin October(a),November(b) and Decmber(c) of 2009

    圖8 2010年10月(a)、11月(b)、12月(c)中西太平洋金槍魚CPUE空間分布與SST關(guān)系Fig.8 CPUE distribution of tuna and relationship with SST in the west-central Pacificin October(a),November(b) and December(c) of 2010

    在年時(shí)間尺度下,整體來(lái)看1995—2002年SST值(在28~29℃之間波動(dòng))明顯低于2003—2010年SST值(在28.5~29.5℃之間波動(dòng)),年CPUE值的變化基本和SST值變化一致。這一方面反映出了20世紀(jì)以來(lái)受全球氣候變暖影響中西太平洋海域海表溫度增加,另一方面也表明了氣候變化影響著中西太平洋金槍魚漁業(yè)[12]。

    根據(jù)郭愛(ài)等[14]、胡奎偉等[15]、唐峰華等[3]研究結(jié)果表明,中西太平洋金槍魚圍網(wǎng)漁場(chǎng)主要分布在SST值為28~32℃的海域,最適SST值為29~31℃,而El Nio和La Nia的發(fā)生引起SST的變化是造成金槍魚圍網(wǎng)漁場(chǎng)中心變動(dòng)的主要原因之一,另外海水溫度的變化還會(huì)帶來(lái)海水葉綠素、營(yíng)養(yǎng)鹽等物質(zhì)的變動(dòng)[1],也會(huì)對(duì)中西太平洋金槍魚圍網(wǎng)漁場(chǎng)中心的變動(dòng)產(chǎn)生影響。

    [1] Matear R J,Chamberlain M A,Sun C,et al. Climate change projection for the western tropical Pacific Ocean using a high-resolution ocean model:implications for tuna fisheries[J]. Deep-Sea Research Part Ⅱ:Topical Studies in Oceanography,2015,113:22-46.

    [2] 朱江峰,戴小杰. 中西太平洋金槍魚圍網(wǎng)漁業(yè)現(xiàn)狀及我國(guó)發(fā)展對(duì)策[J]. 中國(guó)漁業(yè)經(jīng)濟(jì),2009,27(1):115-119.

    Zhu Jiangfeng,Dai Xiaojie. China development stagey on tuna fishery in western and central Pacific[J]. Chinese Fisheries Economics,2009,27(1):115-119.

    [3] 唐峰華,崔雪森,楊勝龍,等. 海洋環(huán)境對(duì)中西太平洋金槍魚圍網(wǎng)漁場(chǎng)影響的GIS時(shí)空分析[J]. 南方水產(chǎn)科學(xué),2014,10(2):18-26.

    Tang Fenghua,Cui Xuesen,Yang Shenglong,et al. GIS analysis on effect of temporal and spatial patterns of marine environment on purse seine fishery in the western and central Pacific[J]. South China Fisheries Science,2014,10(2):18-26.

    [4] 徐冰,陳新軍,田思泉,等. 厄爾尼諾和拉尼娜事件對(duì)秘魯外海莖柔魚漁場(chǎng)分布的影響[J]. 水產(chǎn)學(xué)報(bào),2012,36(5):696-707.

    Xu Bing,Chen Xinjun,Tian Siquan,et al. Effects of El Nio/La Nia on distribution of fishing ground ofDosidicusgigasoff Peru waters[J]. Journal of Fisheries of China,2012,36(5):696-707.

    [5] 周甦芳. 厄爾尼諾-南方濤動(dòng)現(xiàn)象對(duì)中西太平洋鰹魚圍網(wǎng)漁場(chǎng)的影響[J]. 中國(guó)水產(chǎn)科學(xué),2005,12(6):739-744.

    Zhou Sufang. Impacts of the El Nio-Southern Oscillation on skipjack tuna purse-seine fishing grounds in the Western and Central Pacific Ocean[J]. Journal of Fishery Sciences of China,2005,12(6):739-744.

    [6] Lehodey P,Senina I,Nicol S,et al. Modelling the impact of climate change on South Pacific albacore tuna[J]. Deep-Sea Research Part Ⅱ:Topical Studies in Oceanography,2015,113:246-259.

    [7] 杜云艷,周成虎,邵全琴,等. 東海區(qū)海洋漁業(yè)資源環(huán)境的空間聚類分析[J]. 高技術(shù)通訊,2002,12(1):91-95.

    Du Yunyan,Zhou Chenghu,Shao Quanqin,et al. Extracting spatio-temporal patterns from ocean fishery dataset in East China Sea using spatial cluster[J]. Chinese High Technology Letters,2002,12(1):91-95.

    [8] 謝中華. MATLAB統(tǒng)計(jì)分析與應(yīng)用:40個(gè)案例分析[M]. 北京:北京航空航天大學(xué)出版社,2010:292-295.

    Xie Zhonghua. MATLAB Statistical Analysis and Application of Case Studies 40[M]. Beijing:Beijing University of Aeronautics and Astronautics Press,2010:292-295.

    [9] Runyon R P,Haber A,Coleman K A. Behavioral statistics:the Core[M]. Core,1994.

    [10] 郭愛(ài),陳新軍. ENSO與中西太平洋金槍魚圍網(wǎng)資源豐度及其漁場(chǎng)變動(dòng)的關(guān)系[J]. 海洋漁業(yè),2005,27(4):338-342.

    Guo Ai,Chen Xinjun. The relationship between ENSO and tuna purse-seine resource abundance and fishing grounds distribution in the Western and Central Pacific Ocean[J]. Marine Fisheries,2005,27(4):338-342.

    [11] 胡奎偉,許柳雄,陳新軍,等. 海洋遙感在漁場(chǎng)分析中的研究進(jìn)展[J]. 中國(guó)水產(chǎn)科學(xué),2012,19(6):1079-1088.

    Hu Kuiwei,Xu Liuxiong,Chen Xinjun,et al. Application of ocean remote sensing on fishing ground analysis:a review[J]. Journal of Fishery Sciences of China,2012,19(6):1079-1088.

    [12] 王少琴,許柳雄,朱國(guó)平,等. 中西太平洋金槍魚圍網(wǎng)的黃鰭金槍魚CPUE時(shí)空分布及其與環(huán)境因子的關(guān)系[J]. 大連海洋大學(xué)學(xué)報(bào),2014,29(3):303-308.

    Wang Shaoqin,Xu Liuxiong,Zhu Guoping,et al. Spatial-temporal profiles of CPUE and relations to environmental factors for yellowfin tunaThunnusalbacoresfrom purse-seine fishery in Western and Central Pacific Ocean[J]. Journal of Dalian Ocean University,2014,29(3):303-308.

    [13] 靳少非,樊偉. 鰹魚資源開(kāi)發(fā)利用研究現(xiàn)狀及未來(lái)氣候變化背景下研究展望[J]. 漁業(yè)信息與戰(zhàn)略,2014,29(4):272-279.

    Jin Shaofei,F(xiàn)an Wei. Review and perspectives on skipjack tuna fishery under global climate change[J]. Fishery Information & Strategy,2014,29(4):272-279.

    [14] 郭愛(ài),陳新軍,范江濤. 中西太平洋鰹魚時(shí)空分布及其與ENSO關(guān)系探討[J]. 水產(chǎn)科學(xué),2010,29(10):591-596.

    Guo Ai,Chen Xinjun,F(xiàn)an Jiangtao. Spatial and temporal distribution of skipjack and relationship with ENSO in the West-central Pacific Ocean[J]. Fisheries Science,2010,29(10):591-596.

    [15] 胡奎偉,朱國(guó)平,王學(xué)昉,等. 中西太平洋鰹魚豐度的時(shí)空分布及其與表溫的關(guān)系[J]. 海洋漁業(yè),2011,33(4):417-422.

    Hu Kuiwei,Zhu Guoping,Wang Xuefang,et al. Spatio-temporal distribution of skipjack tuna (Katsuwonuspelamis) abundance and its relationship with sea surface temperature in Western and Central Pacific Ocean[J]. Marine Fisheries,2011,33(4):417-422.

    Cluster analysis of tuna purse seine fishery in the Western and Central Pacific

    Hu Qiwei1,Chen Xinjun1,2,3,4,Xu Liangqi5,Yu Jing6

    (1.CollegeofMarineSciences,ShanghaiOceanUniversity,Shanghai201306,China;2.TheKeyLaboratoryofSustainableExploitationofOceanicFisheriesResources,ShanghaiOceanUniversity,MinistryofEducation,Shanghai201306,China; 3.NationalDistant-waterFisheriesEngineeringResearchCenter,Shanghai201306,China; 4.CollaborativeInnovationCenterforDistant-waterFisheries,Shanghai201306,China;5.ShanghaiFisheriesCompany,Shanghai200090,China;6.SouthChinaSeaFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Guangzhou510300,China)

    According to production statistics from tuna purse seine fishery during 1995 to 2010,cluster analysis of purse seine fishery by year and month different time scale is analyzed. Combined with sea surface temperature (SST) and Nio3.4 zone index,the causes of different fishing ground type are discussed. The studies shown that the frequency distribution of abnormal weather events (El Nio and La Nia events)are closely related to the class of clustering results of fishing ground in a monthly time scale,which can be divided into the following phases,i.e. January to March,April to June,July to September and October to December,this results is basically consistent with the time range of clustering results of fishing ground. Compared with center of clustering fishing ground and classification statistical of El Nio and La Nia events,it is found that the type of abnormal climate event associated with the results of clustering. Under the same category,the first category is in 1995 and 1997 with strong El Nio years,the second category is 1998,2007 and 2009 with normal year,the third category is 2010 with a strong La Nia year,the fourth category is 1999,2000 and 2001 with La Nia years,the fifth category is 1996 and 2008 with La Nia Year,and category VI is 2002 and 2004 with El Nio years. It is concluded that the changing annual fishing grounds has a strong correlation with the occurrence of El Nio and La Nia events,the index of the El Nio and La Nia can be used to predict the changes of fishing grounds at year and month scales.

    the Western and Central Pacific; tuna purse seine fishery class; cluster analysis; sea surface temperature; El Nio; La Nia

    2016-06-06;

    2016-07-22。

    海洋局公益性行業(yè)專項(xiàng)(20155014);上海市科技創(chuàng)新計(jì)劃(15DZ1202200)。

    胡啟偉(1992—),男,河南省固始市人,主要研究方向?yàn)闈O業(yè)遙感。E-mail:1140629881@qq.com

    *通信作者:陳新軍(1967—),男,教授。E-mail:xjchen@shou.edu.cn

    10.3969/j.issn.0253-4193.2016.12.007

    S931

    A

    0253-4193(2016)12-0066-10

    胡啟偉,陳新軍,徐良琦,等. 中西太平洋金槍魚圍網(wǎng)漁場(chǎng)聚類及其原因分析[J].海洋學(xué)報(bào),2016,38(12):66—75,

    Hu Qiwei,Chen Xinju,Xu Liangqi,et al. Cluster analysis of tuna purse seine fishery in the Western and Central Pacific [J]. Haiyang Xuebao,2016,38(12):66—75,doi:10.3969/j.issn.0253-4193.2016.12.007

    猜你喜歡
    漁場(chǎng)金槍魚時(shí)間尺度
    時(shí)間尺度上非完整系統(tǒng)的Noether準(zhǔn)對(duì)稱性與守恒量
    金槍魚淹死事件
    可以加速的金槍魚機(jī)器人
    軍事文摘(2021年22期)2022-01-18 06:22:00
    《關(guān)于加強(qiáng)遠(yuǎn)洋魷釣漁船作業(yè)管理的通知》公布
    時(shí)間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    金槍魚與海豚
    交直流混合微電網(wǎng)多時(shí)間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    埃及超大規(guī)模養(yǎng)魚場(chǎng)
    連續(xù)三年保持50%增長(zhǎng),入駐通威、新希望養(yǎng)殖基地,這家漁機(jī)企業(yè)將打造未來(lái)無(wú)人漁場(chǎng)
    國(guó)內(nèi)新型遠(yuǎn)洋金槍魚圍網(wǎng)漁船首航
    商周刊(2018年19期)2018-12-06 09:49:41
    高清日韩中文字幕在线| 国产淫片久久久久久久久 | 国产97色在线日韩免费| av专区在线播放| 又紧又爽又黄一区二区| 丁香欧美五月| 精品免费久久久久久久清纯| 操出白浆在线播放| 日日摸夜夜添夜夜添小说| 久久国产精品影院| 欧美bdsm另类| 色播亚洲综合网| 久久久久性生活片| 精品国内亚洲2022精品成人| 精品国产超薄肉色丝袜足j| 琪琪午夜伦伦电影理论片6080| 在线观看一区二区三区| 91字幕亚洲| 成人特级黄色片久久久久久久| 又紧又爽又黄一区二区| 老司机在亚洲福利影院| 国产精品永久免费网站| x7x7x7水蜜桃| 色精品久久人妻99蜜桃| 夜夜爽天天搞| 夜夜爽天天搞| 成年人黄色毛片网站| 国产成人aa在线观看| 一边摸一边抽搐一进一小说| 亚洲av第一区精品v没综合| 69av精品久久久久久| 天天躁日日操中文字幕| 久久亚洲真实| 日本一二三区视频观看| www.www免费av| 亚洲av免费高清在线观看| 国产欧美日韩一区二区精品| 精品熟女少妇八av免费久了| 在线看三级毛片| 成人鲁丝片一二三区免费| svipshipincom国产片| 欧美日本视频| 丁香欧美五月| 色精品久久人妻99蜜桃| 亚洲在线观看片| 国产视频一区二区在线看| 97人妻精品一区二区三区麻豆| 国产精品自产拍在线观看55亚洲| 成年免费大片在线观看| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 久久6这里有精品| 老司机午夜福利在线观看视频| 搡老熟女国产l中国老女人| 亚洲精品在线观看二区| 久久久精品欧美日韩精品| 国产成人aa在线观看| 99久久精品一区二区三区| 国产高潮美女av| 亚洲18禁久久av| 青草久久国产| 国产色爽女视频免费观看| 日本成人三级电影网站| 婷婷六月久久综合丁香| 亚洲精华国产精华精| 精品日产1卡2卡| 一个人看视频在线观看www免费 | 久久伊人香网站| 波野结衣二区三区在线 | 97超级碰碰碰精品色视频在线观看| 中文字幕av在线有码专区| 日韩人妻高清精品专区| 亚洲av一区综合| 欧美中文日本在线观看视频| 久久国产乱子伦精品免费另类| 午夜a级毛片| 亚洲av日韩精品久久久久久密| 亚洲av免费在线观看| 黄色女人牲交| 色噜噜av男人的天堂激情| 搞女人的毛片| 欧美激情久久久久久爽电影| 九九在线视频观看精品| 不卡一级毛片| 午夜福利18| 久久精品国产亚洲av香蕉五月| 国产精品一区二区三区四区久久| 99国产综合亚洲精品| 精品久久久久久久末码| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 国产一区在线观看成人免费| 国产精华一区二区三区| 国产激情偷乱视频一区二区| 午夜影院日韩av| 香蕉丝袜av| 亚洲精品456在线播放app | 亚洲在线观看片| 国产爱豆传媒在线观看| 我要搜黄色片| 欧美国产日韩亚洲一区| 老汉色∧v一级毛片| 亚洲精品粉嫩美女一区| 一级黄片播放器| 国产亚洲精品av在线| 午夜影院日韩av| 久久6这里有精品| 国产精品 欧美亚洲| 三级男女做爰猛烈吃奶摸视频| 国产精品av视频在线免费观看| 日本 欧美在线| 国产色爽女视频免费观看| 人妻夜夜爽99麻豆av| АⅤ资源中文在线天堂| 桃色一区二区三区在线观看| 婷婷亚洲欧美| www.999成人在线观看| 51国产日韩欧美| 天堂av国产一区二区熟女人妻| 男人舔奶头视频| av在线天堂中文字幕| 欧美3d第一页| 亚洲av五月六月丁香网| 深夜精品福利| 国产麻豆成人av免费视频| 在线观看舔阴道视频| 色在线成人网| 在线观看午夜福利视频| 色综合婷婷激情| 国产午夜精品久久久久久一区二区三区 | 岛国在线免费视频观看| 日韩欧美国产在线观看| 亚洲成人久久性| 综合色av麻豆| 无人区码免费观看不卡| 欧美一级a爱片免费观看看| 亚洲 欧美 日韩 在线 免费| eeuss影院久久| 最近最新中文字幕大全免费视频| 亚洲av成人精品一区久久| 久久久久久九九精品二区国产| 露出奶头的视频| 女人被狂操c到高潮| 国语自产精品视频在线第100页| 91麻豆精品激情在线观看国产| 亚洲av免费高清在线观看| 麻豆国产97在线/欧美| 亚洲激情在线av| 免费搜索国产男女视频| 窝窝影院91人妻| 最好的美女福利视频网| 亚洲精华国产精华精| 欧美一级a爱片免费观看看| 波多野结衣高清无吗| 欧美成人性av电影在线观看| 久久国产精品人妻蜜桃| 精品久久久久久,| 母亲3免费完整高清在线观看| 99热这里只有是精品50| 91av网一区二区| 色精品久久人妻99蜜桃| 在线播放无遮挡| 男女下面进入的视频免费午夜| 欧美成狂野欧美在线观看| 一级毛片高清免费大全| 久久婷婷人人爽人人干人人爱| 一本久久中文字幕| 欧美乱色亚洲激情| www.999成人在线观看| 桃色一区二区三区在线观看| 精品一区二区三区视频在线观看免费| 欧美一区二区亚洲| 男女那种视频在线观看| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 99热精品在线国产| 国产主播在线观看一区二区| 99在线视频只有这里精品首页| 在线播放无遮挡| 婷婷六月久久综合丁香| 中文字幕精品亚洲无线码一区| 亚洲va日本ⅴa欧美va伊人久久| 精品国产美女av久久久久小说| 一级作爱视频免费观看| h日本视频在线播放| 波多野结衣高清作品| 桃色一区二区三区在线观看| av在线蜜桃| 日韩欧美免费精品| 精品国内亚洲2022精品成人| 久久精品国产亚洲av香蕉五月| 别揉我奶头~嗯~啊~动态视频| 无限看片的www在线观看| 国产精品一及| av片东京热男人的天堂| 欧美日韩精品网址| 偷拍熟女少妇极品色| 欧美国产日韩亚洲一区| 日韩欧美精品v在线| 在线观看午夜福利视频| 免费看光身美女| 99精品欧美一区二区三区四区| 在线看三级毛片| 国产老妇女一区| 舔av片在线| www.熟女人妻精品国产| 嫩草影视91久久| 69av精品久久久久久| 免费人成在线观看视频色| 亚洲精品日韩av片在线观看 | 99久久成人亚洲精品观看| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美精品免费久久 | 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| 一区二区三区国产精品乱码| 久久中文看片网| 欧美色欧美亚洲另类二区| 黄色成人免费大全| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 亚洲乱码一区二区免费版| 成人无遮挡网站| 久久精品国产综合久久久| 99久久精品热视频| 国内精品久久久久久久电影| 国产午夜福利久久久久久| 深爱激情五月婷婷| 精品久久久久久久人妻蜜臀av| 精品国产三级普通话版| 欧美黄色淫秽网站| 国内精品美女久久久久久| av视频在线观看入口| 国产欧美日韩精品一区二区| 美女黄网站色视频| 亚洲国产精品久久男人天堂| 国产真实伦视频高清在线观看 | 久久6这里有精品| 成人欧美大片| 久久精品夜夜夜夜夜久久蜜豆| bbb黄色大片| av天堂中文字幕网| 狠狠狠狠99中文字幕| 精品国内亚洲2022精品成人| 精品乱码久久久久久99久播| svipshipincom国产片| 国产av在哪里看| 亚洲五月婷婷丁香| 亚洲专区国产一区二区| 狠狠狠狠99中文字幕| 亚洲av中文字字幕乱码综合| 在线免费观看的www视频| 老司机午夜福利在线观看视频| 91久久精品电影网| 一边摸一边抽搐一进一小说| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 国产精品一及| 国产精品乱码一区二三区的特点| 亚洲国产精品sss在线观看| 国产一区二区亚洲精品在线观看| 两人在一起打扑克的视频| 91在线观看av| 少妇人妻一区二区三区视频| 听说在线观看完整版免费高清| 日韩有码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三区av网在线观看| 手机成人av网站| 久久婷婷人人爽人人干人人爱| 嫩草影院入口| 国产精品久久久久久久久免 | 好男人在线观看高清免费视频| 成人精品一区二区免费| 成年女人看的毛片在线观看| 色哟哟哟哟哟哟| 床上黄色一级片| 欧美黄色淫秽网站| 精品福利观看| 他把我摸到了高潮在线观看| tocl精华| 一级毛片女人18水好多| 久久精品国产综合久久久| 久久久久久大精品| 操出白浆在线播放| 搡老熟女国产l中国老女人| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩卡通动漫| 99久久综合精品五月天人人| 精品不卡国产一区二区三区| 国产精品三级大全| 在线免费观看不下载黄p国产 | 亚洲18禁久久av| 免费在线观看亚洲国产| 一个人看的www免费观看视频| 国产精品久久久人人做人人爽| 69av精品久久久久久| 99久久成人亚洲精品观看| 久久香蕉精品热| www.999成人在线观看| 亚洲成av人片免费观看| 最后的刺客免费高清国语| 国产伦人伦偷精品视频| av片东京热男人的天堂| 欧美激情久久久久久爽电影| 欧美av亚洲av综合av国产av| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 黑人欧美特级aaaaaa片| 国产99白浆流出| 五月玫瑰六月丁香| 精品熟女少妇八av免费久了| 美女被艹到高潮喷水动态| 人人妻人人澡欧美一区二区| 欧美性猛交╳xxx乱大交人| 给我免费播放毛片高清在线观看| 亚洲国产中文字幕在线视频| 悠悠久久av| 亚洲一区高清亚洲精品| 日韩大尺度精品在线看网址| 亚洲专区国产一区二区| 欧美午夜高清在线| 一区二区三区免费毛片| 天堂影院成人在线观看| 女同久久另类99精品国产91| 午夜福利18| 最近视频中文字幕2019在线8| 久久九九热精品免费| 亚洲av电影不卡..在线观看| 国产午夜精品论理片| 好男人电影高清在线观看| 免费高清视频大片| 日韩高清综合在线| 亚洲人成网站在线播| 无限看片的www在线观看| 精品人妻1区二区| 亚洲熟妇熟女久久| 久久性视频一级片| 欧洲精品卡2卡3卡4卡5卡区| 国产精品99久久久久久久久| 男人舔奶头视频| 成人特级av手机在线观看| 成年人黄色毛片网站| 天堂影院成人在线观看| 亚洲色图av天堂| 丁香欧美五月| 超碰av人人做人人爽久久 | 一级黄片播放器| 国产av麻豆久久久久久久| 不卡一级毛片| 丰满人妻一区二区三区视频av | 久久伊人香网站| 国产高潮美女av| 欧美一区二区国产精品久久精品| www国产在线视频色| 白带黄色成豆腐渣| 国产精品久久久久久精品电影| 十八禁人妻一区二区| 麻豆一二三区av精品| 亚洲成人中文字幕在线播放| 欧美成人a在线观看| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 天堂动漫精品| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 国产免费一级a男人的天堂| 国产成年人精品一区二区| 十八禁人妻一区二区| 久久国产精品人妻蜜桃| 亚洲电影在线观看av| 国产精品久久久久久精品电影| 久久久久九九精品影院| 一本精品99久久精品77| 国内精品美女久久久久久| 国产亚洲欧美在线一区二区| 美女cb高潮喷水在线观看| 国内精品久久久久精免费| 国产亚洲欧美98| 欧美精品啪啪一区二区三区| 啦啦啦免费观看视频1| 欧美zozozo另类| 国产成人av激情在线播放| 久久人人精品亚洲av| 禁无遮挡网站| 级片在线观看| 99久久久亚洲精品蜜臀av| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 国产黄片美女视频| 哪里可以看免费的av片| 国产视频内射| 亚洲色图av天堂| 午夜a级毛片| 欧美黑人巨大hd| av黄色大香蕉| 成人特级黄色片久久久久久久| 88av欧美| 国产成人aa在线观看| 在线观看舔阴道视频| 亚洲精品一区av在线观看| 亚洲黑人精品在线| 国产在线精品亚洲第一网站| 日韩欧美 国产精品| 欧美日韩乱码在线| 99久久成人亚洲精品观看| 成年女人看的毛片在线观看| 麻豆国产97在线/欧美| 两人在一起打扑克的视频| 欧美精品啪啪一区二区三区| 男女下面进入的视频免费午夜| 久久天躁狠狠躁夜夜2o2o| 观看免费一级毛片| 午夜福利高清视频| 偷拍熟女少妇极品色| 亚洲精品乱码久久久v下载方式 | 丰满乱子伦码专区| 国产精品久久久久久久久免 | aaaaa片日本免费| 国产伦一二天堂av在线观看| 日韩欧美一区二区三区在线观看| 禁无遮挡网站| 国内精品一区二区在线观看| 99精品欧美一区二区三区四区| 91麻豆av在线| 成人av一区二区三区在线看| 超碰av人人做人人爽久久 | 色综合婷婷激情| 午夜免费观看网址| 在线观看免费午夜福利视频| 久久久久免费精品人妻一区二区| 人人妻人人澡欧美一区二区| 高清毛片免费观看视频网站| 国产伦一二天堂av在线观看| 精品国产美女av久久久久小说| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 色综合欧美亚洲国产小说| 91久久精品电影网| 亚洲黑人精品在线| 99久久99久久久精品蜜桃| 国产精品一及| 国产淫片久久久久久久久 | 日韩亚洲欧美综合| 中文资源天堂在线| 午夜福利18| 99riav亚洲国产免费| 欧美av亚洲av综合av国产av| 88av欧美| 亚洲欧美激情综合另类| 成年版毛片免费区| 亚洲内射少妇av| 精品日产1卡2卡| 两个人的视频大全免费| 日本免费a在线| xxx96com| 欧美av亚洲av综合av国产av| 毛片女人毛片| 欧美最新免费一区二区三区 | 欧美日韩黄片免| 少妇裸体淫交视频免费看高清| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 午夜免费激情av| 免费av观看视频| 97超视频在线观看视频| 成熟少妇高潮喷水视频| 男女午夜视频在线观看| 国产精品嫩草影院av在线观看 | 日本与韩国留学比较| 真人一进一出gif抽搐免费| 欧美不卡视频在线免费观看| 午夜福利免费观看在线| 真实男女啪啪啪动态图| 国产主播在线观看一区二区| 很黄的视频免费| 网址你懂的国产日韩在线| 19禁男女啪啪无遮挡网站| 老司机福利观看| 一区二区三区高清视频在线| 69人妻影院| 啦啦啦韩国在线观看视频| 亚洲精品成人久久久久久| 午夜福利视频1000在线观看| 天堂av国产一区二区熟女人妻| 国内揄拍国产精品人妻在线| 成人性生交大片免费视频hd| 国产精品一区二区三区四区免费观看 | 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 国产一区二区激情短视频| 亚洲不卡免费看| 大型黄色视频在线免费观看| 精品久久久久久久末码| 熟妇人妻久久中文字幕3abv| 欧美黑人欧美精品刺激| 精品日产1卡2卡| 最好的美女福利视频网| 久久久国产成人免费| 欧美日韩黄片免| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av涩爱 | 夜夜躁狠狠躁天天躁| 91av网一区二区| 久99久视频精品免费| 变态另类成人亚洲欧美熟女| 国产精品自产拍在线观看55亚洲| 淫妇啪啪啪对白视频| 中文字幕人成人乱码亚洲影| 有码 亚洲区| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 国模一区二区三区四区视频| 国产伦一二天堂av在线观看| 哪里可以看免费的av片| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 午夜福利成人在线免费观看| 亚洲自拍偷在线| 国产视频内射| 欧美+亚洲+日韩+国产| 亚洲成人久久爱视频| 人妻夜夜爽99麻豆av| 俺也久久电影网| 黄色日韩在线| 国产伦人伦偷精品视频| 成人精品一区二区免费| 亚洲国产欧洲综合997久久,| 日韩欧美国产在线观看| avwww免费| 欧美中文日本在线观看视频| 成人av一区二区三区在线看| a在线观看视频网站| 亚洲人成网站在线播放欧美日韩| 成人午夜高清在线视频| 国产精品乱码一区二三区的特点| 国产黄色小视频在线观看| 中文资源天堂在线| or卡值多少钱| 日韩欧美在线乱码| 国产一区二区亚洲精品在线观看| 日本一本二区三区精品| 日本 av在线| 91麻豆av在线| 最好的美女福利视频网| 一进一出抽搐动态| 色哟哟哟哟哟哟| 男插女下体视频免费在线播放| 亚洲av五月六月丁香网| 欧美日韩乱码在线| 亚洲国产欧美网| 国产一区二区激情短视频| 国产成人a区在线观看| 欧美午夜高清在线| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 桃色一区二区三区在线观看| 伊人久久大香线蕉亚洲五| 成人性生交大片免费视频hd| 日本黄色视频三级网站网址| 亚洲国产精品成人综合色| 国产精品乱码一区二三区的特点| 欧美成人性av电影在线观看| 欧美日韩瑟瑟在线播放| 免费搜索国产男女视频| 国产精品98久久久久久宅男小说| 亚洲精品影视一区二区三区av| 99视频精品全部免费 在线| 国产精品久久视频播放| 岛国在线观看网站| 中亚洲国语对白在线视频| 国产精品久久久久久久电影 | 免费一级毛片在线播放高清视频| 国产国拍精品亚洲av在线观看 | 欧美乱色亚洲激情| 综合色av麻豆| 9191精品国产免费久久| 90打野战视频偷拍视频| 国产av麻豆久久久久久久| 天堂√8在线中文| 美女被艹到高潮喷水动态| 亚洲国产欧美人成| 小蜜桃在线观看免费完整版高清| 国产一区二区激情短视频| 国产av麻豆久久久久久久| 久久久久久国产a免费观看| 国产精品野战在线观看| 欧美日韩中文字幕国产精品一区二区三区| av女优亚洲男人天堂| 欧美最新免费一区二区三区 | 美女大奶头视频| 欧美中文日本在线观看视频| 搡老岳熟女国产| 婷婷丁香在线五月| 最新中文字幕久久久久| 国产av一区在线观看免费| 丰满人妻一区二区三区视频av | 欧美区成人在线视频| 99riav亚洲国产免费| 熟女电影av网| 一本久久中文字幕| 成人av一区二区三区在线看| 三级毛片av免费| 日韩欧美精品v在线| 午夜福利成人在线免费观看| 男女做爰动态图高潮gif福利片| 国产在视频线在精品| 丰满人妻一区二区三区视频av | 精品国产亚洲在线| 亚洲最大成人手机在线| 精品一区二区三区人妻视频| 午夜精品久久久久久毛片777| 日本免费一区二区三区高清不卡| 国产成人av教育|