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

    基于GAM模型研究水溫垂直結(jié)構(gòu)對熱帶中西太平洋黃鰭金槍魚漁獲率的影響

    2021-07-06 08:10:20楊勝龍范秀梅吳祖立伍玉梅戴陽
    海洋學(xué)報(bào) 2021年4期
    關(guān)鍵詞:黃鰭溫躍層環(huán)境變量

    楊勝龍,范秀梅,吳祖立,伍玉梅,戴陽*

    (1.中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所農(nóng)業(yè)農(nóng)村部遠(yuǎn)洋與極地漁業(yè)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室,上海200090;2.中國水產(chǎn)科學(xué)研究院漁業(yè)資源與遙感信息技術(shù)重點(diǎn)開放實(shí)驗(yàn)室,上海200090)

    1 引言

    黃鰭金槍魚(Thunnusobesus)是遠(yuǎn)洋漁業(yè)重要經(jīng)濟(jì)魚種,其索餌水層分布受海洋環(huán)境和餌料分布的共同影響[1],進(jìn)而影響漁船捕撈效率[2–3],因此有大量關(guān)于黃鰭金槍魚垂直游動水層的研究[4–7]。海洋環(huán)境還會影響金槍魚及其餌料生物的適宜垂直棲息空間選擇,故大量研究從海洋環(huán)境對金槍魚適宜生境的影響角度進(jìn)行分析[7–9]。Deary等[7]研究了黃鰭金槍魚在最適宜水溫(17.5 ℃)、溶解氧(107.4 mol/kg)和極限水溫(12℃)、溶解氧(43.6 mol/kg)條件下,其垂直活動空間分布特征,結(jié)果表明,太平洋黃鰭金槍魚垂直活動水層存在明顯的空間分布特征。而在溶解氧不受限制的海域,黃鰭金槍魚垂直游動主要受水溫垂直結(jié)構(gòu)的影響[10],其大部分時(shí)間分布在低于海表溫度8℃以內(nèi)的海域(距海表水溫8℃,定義為△8℃)[11]。目前,相關(guān)文獻(xiàn)多側(cè)重于黃鰭金槍魚水層分布及環(huán)境影響的研究,且有關(guān)環(huán)境對延繩釣漁獲率影響分析多采用表層遙感數(shù)據(jù)[12],次表層環(huán)境數(shù)據(jù)對延繩釣漁獲率影響研究報(bào)道較少[13]。

    在相同的漁業(yè)資源量條件下,金槍魚棲息的生境不同會導(dǎo)致金槍魚分布密度不同。在金槍魚適宜棲息的垂直空間被壓縮區(qū)域,魚群分布相對密集,有利于表層漁船捕撈,漁獲率高;反之漁獲率低[8–9]。因此,金槍魚棲息生境對單位捕撈努力量漁獲量(Catch Per Unite Effort,CPUE)與資源量的回歸曲線有重要影響[8–9]。依據(jù)魚群活動水層改變作業(yè)深度有可能提高漁獲率,金槍魚游動水層對商業(yè)CPUE值亦有影響[1]。因此,僅了解表層環(huán)境變量對CPUE的水平空間影響是不夠的,還應(yīng)該考慮金槍魚適宜棲息水層[9]。研究中西太平洋黃鰭金槍魚適宜生境空間對延繩釣漁獲率影響的文獻(xiàn)較少,因此本文在前人研究基礎(chǔ)上,以次表層環(huán)境變量作為黃鰭金槍魚適宜生境范圍指標(biāo),用模型分析垂直水溫結(jié)構(gòu)對熱帶中西太平洋黃鰭金槍魚延繩釣漁獲率的影響,為延繩釣金槍魚捕撈作業(yè)和資源養(yǎng)護(hù)提供理論支持。

    2 方法與材料

    2.1 數(shù)據(jù)和處理方法

    本 文 的 研 究 區(qū) 域 為25°S~25°N,130°E~130°W(圖1),漁業(yè)數(shù)據(jù)采用中西太平洋漁業(yè)委員會(Western and Central Pacific Fishery Comm ission,WCPFC)2007?2017年的延繩釣生產(chǎn)統(tǒng)計(jì)數(shù)據(jù)(分辨率:5°×5°),漁獲率采用CPUE(單位:尾/千鉤)來表示,按下式計(jì)算:

    圖1 2007?2017年中西太平洋延繩釣黃鰭金槍魚CPUE分布Fig.1 CPUE distribution of longline yellow fin tuna in western and central Pacific from 2007 to 2017

    式中,CPUE(i,j)、Nfish(i,j)、Nhook(i,j)分別為第i個(gè)經(jīng)度、第j個(gè)緯度處方格的月平均值、月釣獲總尾數(shù)和月總投鉤數(shù)。

    Deary等[7]綜合大量文獻(xiàn),認(rèn)為10.5 ℃為太平洋黃鰭金槍魚生理容忍水溫下限值,17.5 ℃為黃鰭金槍魚適宜水溫值。標(biāo)志放流結(jié)果表明,東太平洋黃鰭金槍魚白天頻繁進(jìn)入12℃等溫線覓食?!?℃被認(rèn)為是限制黃鰭金槍魚垂直分布的關(guān)鍵因子[10]。溫躍層是影響金槍魚適宜垂直游動和索餌的關(guān)鍵環(huán)境因子[1]。因此,本文選擇次表層水溫(12~18℃)、垂直溫差(△8℃)和溫躍層作為影響黃鰭金槍魚漁獲率的備選環(huán)境變量。同時(shí),考慮到水溫和溫躍層對延繩釣漁獲率的影響有可能是交叉的,因此,本文環(huán)境變量選擇了△8℃等溫線深度與溫躍層下界深度的差值(sdc)以及17℃和18℃等溫線深度和溫躍層下界深度的差值(sdc1、sdc2),綜合分析了中西太平洋黃鰭金槍魚適宜生境空間對漁獲率的影響。

    采用Argo浮標(biāo)數(shù)據(jù)(來源于http://www.argo.org.cn/)計(jì)算次表層水溫、垂直溫差和溫躍層月平均產(chǎn)品數(shù)據(jù)(分辨率:1°×1°),計(jì)算方法參考文獻(xiàn)[13]。對所有的環(huán)境產(chǎn)品數(shù)據(jù),按月平均,即計(jì)算2007?2017多年月平均數(shù)據(jù)。同時(shí)為與黃鰭金槍魚延繩釣捕撈數(shù)據(jù)匹配,將所有分辨率為1°×1°的網(wǎng)格環(huán)境參數(shù)值轉(zhuǎn)換成5°×5°的空間分辨率。

    2.2 模型方法

    GAM模型可以模擬分析輸入因子對CPUE的非線性影響作用[12]。CPUE的變化可能是變量的交叉影響結(jié)果,包括隨時(shí)空變化的環(huán)境變量和資源量。本文構(gòu)建時(shí)空GAM模型和環(huán)境GAM模型,分別討論CPUE時(shí)空分布特征和環(huán)境變量對CPUE的影響。時(shí)間變量有年和月,空間變量有經(jīng)度和緯度。次表層環(huán)境變量有次表層水溫(12~18℃)、垂直溫差(△8℃)和溫躍層以及△8℃和17℃與溫躍層下界深度差,共計(jì)17個(gè)輸入變量。所有變量自由度為4,采用R軟件包的mgcv函數(shù)庫構(gòu)建GAM模型進(jìn)行計(jì)算。采用逐步添加變量的方法構(gòu)建模型,以0.05 顯著性水平選擇影響變量,最優(yōu)模型選擇采用AIC(Akaike Information Criterion)值。公式(1)計(jì)算的CPUE值有部分為0,一般給CPUE加上1個(gè)常數(shù)。本文對所有CPUE加上常數(shù)1后取對數(shù)變換,構(gòu)建的GAM模型方程為

    式中,e為模型設(shè)差;a為常數(shù)。

    采用特征根判定法對輸入變量多重共線性進(jìn)行檢驗(yàn)。依據(jù)首次構(gòu)建的環(huán)境GAM模型解釋方差,剔除解釋方差較小的變量,直到kappa系數(shù)小于1 000,最后構(gòu)建最終的環(huán)境GAM模型。

    3 結(jié)果

    3.1 GAM模型結(jié)果

    時(shí)空變量GAM模型結(jié)果見表1 和表2 ,模型AIC值為19 965.13 ;偏差解釋率為45.3 %(表2)。從方差和F值看,經(jīng)度(lon)和緯度(lat)對中西太平洋黃鰭金槍魚延繩釣漁獲率影響很大,年(year)和月(month)對延繩釣漁獲率影響相比很?。ū?)。

    表1 中西太平洋黃鰭金槍魚時(shí)空變量GAM模型F-檢驗(yàn)值Table 1 F-test value of the spatio-tem poral GAM m odel of yellow fin tuna in western and central Pacific

    表2 中西太平洋黃鰭金槍魚時(shí)空變量GAM模型統(tǒng)計(jì)參數(shù)Table 2 Statistical characteristics of the spatio-tem poral GAM m odel of yellow fin tuna in western and central Pacific

    環(huán)境變量剔除了13~17℃等溫線深度、△8℃等溫線深度、溫躍層下界與17℃深度差(sdc1)和溫躍層下界與18℃深度差(sdc2)等變量,保留了溫躍層上界深度(upsd)、溫躍層上界溫度(upwd)、溫躍層下界深度(downsd)、溫躍層下界溫度(downwd)、溫躍層厚度(hd)、溫躍層強(qiáng)度(qd)等6個(gè)溫躍層參數(shù)以及12℃深度(D12)、18℃深度(D18)和△8℃等溫線深度與溫躍層下界深度差(sdc)。此時(shí)kappa系數(shù)是236,不存在嚴(yán)重的多重共線性。

    基于保留的變量構(gòu)建的環(huán)境GAM模型結(jié)果見表3 和表4 。在p=0.05 的顯著檢驗(yàn)水平下,所有環(huán)境變量通過顯著性檢驗(yàn),此時(shí)模型AIC值最小,為23 784.19。次表層環(huán)境變量對中西太平洋黃鰭金槍魚延繩釣漁獲率存在非線性影響。所有次表層環(huán)境變量對中西太平洋黃鰭金槍魚延繩釣CPUE偏差解釋率為28.1 %,模型擬合的決定系數(shù)為0.2 7 9 。從方差和F值看,溫躍層上界溫度和深度、溫躍層下界深度、18℃等溫線深度、△8℃等溫線深度及其和溫躍層下界深度差對延繩釣漁獲率影響較大。溫躍層厚度、強(qiáng)度、溫躍層下界溫度和12℃等溫線深度對中西太平洋黃鰭金槍魚延繩釣漁獲率影響較低。

    表3 次表層環(huán)境變量GAM模型統(tǒng)計(jì)參數(shù)Table3 Statistical characteristics of the spatio-tem poral GAM model for subsurface environm ental variables

    表4 次表層環(huán)境變量GAM模型檢驗(yàn)值Table 4 GAM m odel test value of subsurface environm ental variables

    3.2 CPUE時(shí)空分布特征

    圖2 表明,熱帶中西太平洋黃鰭金槍魚延繩釣CPUE影響在2012年之前維持在較低的水平,2012年之后CPUE開始上升。不同月份熱帶中西太平洋黃鰭金槍魚延繩釣CPUE呈現(xiàn)穹頂狀分布。CPUE從1月一直增長至6月,達(dá)到全年最高,之后逐漸下降。熱帶中西太平洋黃鰭金槍魚延繩釣CPUE在赤道以南的值要高于赤道以北的值,高值CPUE集中出現(xiàn)在10°S。隨著緯度向北,黃鰭金槍魚延繩釣CPUE值逐漸遞減。圖2 d顯示,熱帶中西太平洋西部區(qū)域黃鰭金槍魚延繩釣CPUE要大于東部區(qū)域。CPUE最大出現(xiàn)在140°E區(qū)域,隨著經(jīng)度由西向東,CPUE值逐漸變小。對CPUE影響密切的空間區(qū)域在10°S,140°E附近。

    圖2 時(shí)空變量對中西太平洋黃鰭金槍魚延繩釣CPUE的影響Fig.2 Theeffectsof spatial-temporalpredictorson CPUE of western and central Pacific

    3.3 次表層環(huán)境因子對CPUE影響

    在p等于0.05 的顯著檢驗(yàn)水平下,溫躍層對熱帶中西太平洋黃鰭金槍魚延繩釣CPUE存在顯著影響(圖3)。溫躍層上界深度和溫度對中西太平洋黃鰭金槍魚延繩釣CPUE影響呈現(xiàn)遞增趨勢,隨著溫躍層上界溫度和深度值變大,延繩釣CPUE逐漸遞增,對延繩釣CPUE影響密切的溫度和深度范圍分別為27~28℃和70~90m(圖3 a,圖3 b)。溫躍層下界深度在250m以淺時(shí),對延繩釣CPUE影響是微弱震蕩,在250m時(shí)達(dá)到最大,此時(shí)下界深度對延繩釣CPUE影響也是最密切的;之后隨著下界深度的變大,延繩釣CPUE值快速降低(圖3 c)。溫躍層下界溫度對延繩釣CPUE影響呈現(xiàn)非線性遞減趨勢。隨著溫躍層下界溫度變大,延繩釣CPUE值逐漸變低。溫躍層下界溫度對延繩釣CPUE影響密切的溫度范圍為13~16℃。溫躍層厚度對延繩釣CPUE影響呈現(xiàn)先震蕩遞減再增加的趨勢。溫躍層厚度對延繩釣CPUE影響密切的深度范圍為40~70m(圖3 e)。溫躍層強(qiáng)度對延繩釣CPUE影響呈現(xiàn)震蕩遞減趨勢,溫躍層強(qiáng)度在0.1 ℃/m時(shí)對延繩釣CPUE影響密切,之后隨著溫躍層強(qiáng)度變大,延繩釣CPUE逐漸變?。▓D3 f)。

    12℃和18℃等溫線深度對熱帶中西太平洋黃鰭金槍魚延繩釣CPUE影響相似,對延繩釣CPUE影響呈現(xiàn)先震蕩后遞增趨勢(圖3 g,圖3 h)。12℃等溫線深度對延繩釣CPUE影響密切的深度在340m深度左右;18℃等溫線深度影響密切的深度在230m深度左右?!?℃等溫線深度與溫躍層下界深度的差值對熱帶中西太平洋黃鰭金槍魚延繩釣CPUE影響呈現(xiàn)先快速遞減后緩慢增加,在深度差為70m左右影響最密切但影響最低(圖3 i)。

    4 討論

    4.1 因子對漁獲率影響顯著性分析

    海洋生態(tài)系統(tǒng)的改變會影響金槍魚的水平?垂直分布、洄游、幼魚存活率和補(bǔ)充強(qiáng)度,從而在不同的空間尺度上影響漁業(yè)捕撈和漁獲率[14]。金槍魚的垂直分布主要受海洋環(huán)境和餌料生物影響,海洋環(huán)境主要包括水溫、溫躍層和溶解氧[1]。海洋環(huán)境和餌料生物呈現(xiàn)時(shí)空變化特征,因此金槍魚漁船捕撈作業(yè)和漁獲率也都呈現(xiàn)時(shí)空變化的規(guī)律,年、月和經(jīng)緯度是CPUE時(shí)空變化的反映。金槍魚的垂直分布水層是否和延繩釣投鉤深度水層吻合,直接影響延繩釣金槍魚捕撈效率[2–3]。因此,本文從次表層環(huán)境影響黃鰭金槍魚垂直游動,進(jìn)而影響延繩釣漁獲率的角度,構(gòu)建了次表層環(huán)境變量對CPUE的非線性響應(yīng)模型,討論了水溫垂直結(jié)構(gòu)對延繩釣黃鰭金槍魚CPUE的影響。

    CPUE變量可能是許多相互交叉的因子的函數(shù),包括環(huán)境變量的時(shí)空變化、捕撈能力和資源量。漁業(yè)資源空間分布實(shí)際受環(huán)境變量影響,環(huán)境變量的空間分布驅(qū)動漁業(yè)資源時(shí)空變動。構(gòu)建2個(gè)模型,分別探討了中西太平洋延繩釣黃鰭金槍魚CPUE時(shí)空分布特征和環(huán)境變量對CPUE的影響。第一個(gè)GAM模型只包含時(shí)空變量,不包括環(huán)境變量,模型結(jié)果揭示了CPUE時(shí)空分布特征。第二個(gè)模型只包含環(huán)境變量,揭示了環(huán)境變量本身分布對CPUE的影響。Briand等[14]采用相似方法分析了南太平洋長鰭金槍魚延繩釣漁獲率變化。時(shí)空模型方差和F值顯示(表1 ,表2 ),CPUE呈現(xiàn)明顯的空間分布特征,CPUE高值主要出現(xiàn)在西南區(qū)域。CPUE隨時(shí)間變化并不明顯,表明各年和各月CPUE差異不大。

    環(huán)境模型結(jié)果表明,溫躍層上界溫度和深度、溫躍層下界深度、18℃等溫線深度、△8℃等溫線深度及其和溫躍層下界深度差對延繩釣黃鰭金槍魚漁獲率有影響。溫躍層是水溫垂直結(jié)構(gòu)中一個(gè)關(guān)鍵的環(huán)境因素,它不僅影響著金槍魚的垂直水層分布,也影響著金槍魚的餌料生物垂直水層分布[1]。晚上黃鰭金槍魚棲息在溫躍層以上暖水層,白天黃鰭金槍魚下游到溫躍層內(nèi)部靠近下界區(qū)域覓食深水散射層(Deep Scattering layer,DSL)上層海洋生物[15]。溫躍層上界深度和溫度以及溫躍層下界深度對漁獲率影響較大。黃鰭金槍魚有能力突破溫躍層進(jìn)入深水層,但這種行為很少[4–5]。這可能是本文模型溫躍層強(qiáng)度、厚度和溫躍層下界溫度相比影響較小原因。次表層溫度對熱帶中西太平洋黃鰭金槍魚延繩釣CPUE也有影響。Deary等[7]以17.5 ℃為黃鰭金槍魚適宜棲息水溫,分析了其適宜分布的垂直深度,同時(shí)以10.5 ℃為太平洋黃鰭金槍魚容忍水溫的下限值,繪制了其最大垂直下游深度空間分布。18℃接近太平洋黃鰭金槍魚適宜棲息水溫是其影響相對較大原因?!?℃被認(rèn)為是限制黃鰭金槍魚垂直下游的溫差[10],因此△8℃等溫線深度與溫躍層下界深度的差對漁獲率影響顯著。

    4.2 關(guān)鍵次表層環(huán)境對黃鰭金槍魚延繩釣漁獲率的影響

    延繩釣黃鰭金槍魚捕獲深水覓食的散射層生物,其捕撈效率不僅與漁業(yè)資源空間分布有關(guān),還與投鉤深度是否與黃鰭金槍魚索餌活動水層吻合及餌料生物分布密度有關(guān)[2]。延繩釣漁船投鉤深度一般保持不變,因此深水覓食的黃鰭金槍魚索餌水層分布對表層漁船捕撈效率有重要影響。

    影響黃鰭金槍魚索餌水層分布的2個(gè)關(guān)鍵因子是餌料生物深度[15]和環(huán)境因子[16–17]。黃鰭金槍魚水平空間和垂直水層分布均受餌料生物影響。延繩釣黃鰭金槍魚漁獲率與餌料生物密度在水平空間呈正相關(guān)關(guān)系[1]。在餌料生物分布區(qū)域,其水層分布對黃鰭金槍魚索餌水層和表層漁船作業(yè)有重要影響。研究表明,白天黃鰭金槍魚會頻繁下游到深層冷水區(qū)域覓食DSL上層的海洋生物[15],索餌水層多在溫躍層內(nèi)部靠近下界區(qū)域[18],因此DSL的水層分布直接影響黃鰭金槍魚索餌水層分布。在太平洋,季節(jié)性DSL晝夜垂直分布和溫躍層有關(guān),DSL上層的海洋生物白天分布在溫躍層以下區(qū)域,晚上分布在溫躍層以上水域[19]。太平洋DSL分布在200~500m水層,其水層深度與溫躍層深度呈正相關(guān)。當(dāng)溫躍層變淺時(shí),DSL垂直分布也變淺,黃鰭金槍魚索餌水層變淺;反之,溫躍層變深時(shí),DSL分布水層也變深[19],黃鰭金槍魚索餌水層變深。溫躍層通過影響DSL水層對黃鰭金槍魚索餌水層有直接和間接影響。因此,不同的海域因溫躍層分布不同,黃鰭金槍魚白天索餌水層會有空間差異,對延繩釣漁獲率有重要影響。此外,局部過于密集的生物分布并不利于延繩釣漁船捕撈。金槍魚類都依靠視覺進(jìn)行捕食,在餌料生物過于密集的水層,延繩釣投鉤的誘餌并不占優(yōu)勢,黃鰭金槍魚會更優(yōu)先捕食餌料生物,導(dǎo)致漁獲率低[1]。

    黃鰭金槍魚覓食時(shí),下游的水層深度不但受限于自身生理功能,還受外界環(huán)境因子影響。黃鰭金槍魚下游最大深度受其生理承受的最低水溫和溶解氧限制。實(shí)驗(yàn)表明,在周邊水溫急劇下降超過10℃時(shí),黃鰭金槍魚不能有效提升自身肌體溫度,因此△8℃被認(rèn)為是限制黃鰭金槍魚垂直分布的關(guān)鍵因子[10]。本研究結(jié)果表明,△8℃等溫線深度與溫躍層下界深度的差對漁獲率影響更顯著,表明△8℃等溫線深度對CPUE的作用和溫躍層下界深度的作用是交叉的。

    本文模型結(jié)果表明,溫躍層上界溫度和深度、溫躍層下界深度、18℃等溫線深度、△8℃等溫線深度及其和溫躍層下界深度差對延繩釣漁獲率影響較大,是影響熱帶中西太平洋黃鰭金槍魚延繩釣漁獲率的關(guān)鍵環(huán)境因子。溫躍層下界深度通過影響?zhàn)D料生物水層分布,進(jìn)而影響黃鰭金槍魚索餌水層。圖3 表明,在溫躍層下界深度逐漸變大到250m時(shí),CPUE都比較大,這與中西太平洋延繩釣投鉤深度主要覆蓋70~300m[20]有關(guān)。在溫躍層下界深度為250m上下時(shí),黃鰭金槍魚在該水層附近覓食餌料生物,索餌水層和延繩釣多數(shù)釣鉤深度吻合。因此在250m深度區(qū)域,延繩釣黃鰭金槍魚漁獲率最大。下界深度超過300m時(shí),黃鰭金槍魚垂直活動空間大、分布不集中和索餌水層深都不利于表層漁船捕撈,漁獲率低。18℃等溫線多分布在溫躍層內(nèi)部,揭示黃鰭金槍魚適宜棲息的水溫深度在投鉤的主要深度區(qū)間,棲息深度和投鉤深度吻合時(shí)有較好的漁獲率。溫躍層上界溫度不影響黃鰭金槍魚垂直水層分布,但黃鰭金槍魚是暖水性魚類,棲息和產(chǎn)卵需要在一定的水溫之上才可進(jìn)行。厄爾尼諾和拉尼娜事件影響著太平洋的溫躍層變化,太平洋金槍魚空間分布發(fā)生改變[7]。中西太平洋黃鰭金槍魚圍網(wǎng)產(chǎn)量在29℃最高[21]。Lan等[12]采用海洋表層環(huán)境變量分析了中西太平洋黃鰭金槍魚CPUE變化并進(jìn)行預(yù)測。延繩釣CPUE隨著海表溫度升高而變大,影響密切區(qū)域的溫度為25~30℃。本文模型結(jié)果揭示,溫躍層上界溫度對CPUE影響密切的溫度為27~28℃,與黃鰭金槍魚喜好棲息在表層28℃有關(guān)。溫躍層上界深度為70~90m。在中西太平洋,黃鰭金槍魚幼魚夜間分布都淺于100m,但比東太平洋黃鰭金槍魚垂直分布深,全年多數(shù)月份都在溫躍層之上水域,少數(shù)時(shí)間游到溫躍層以下[22]。

    圖3 次表層環(huán)境變量對中西太平洋黃鰭金槍魚延繩釣CPUE的影響Fig.3 The effects of subsurface environmental variables on CPUE of yellow tuna in the western and central Pacific

    4.3 次表層環(huán)境對黃鰭金槍魚延繩釣CPUE標(biāo)準(zhǔn)化的啟示

    CPUE經(jīng)常作為資源豐度指標(biāo)用于漁業(yè)資源分析,CPUE的商業(yè)生產(chǎn)數(shù)據(jù)受各種因素影響(如時(shí)間、空間、漁船參數(shù)和環(huán)境變量等)[23]。因此,名義CPUE值并不能很好地反映真實(shí)的資源情況[3],一般會對CPUE進(jìn)行標(biāo)準(zhǔn)化[24]。在不同的水溫垂直結(jié)構(gòu)水域,CPUE對資源量的回歸直線斜率是不同的[8?9]。當(dāng)標(biāo)準(zhǔn)化CPUE沒有考慮金槍魚生活習(xí)性時(shí),也會產(chǎn)生估計(jì)偏差[25],從而無法準(zhǔn)確對漁場和資源進(jìn)行評估。如在法屬波利尼亞馬克薩斯群島海域,調(diào)查發(fā)現(xiàn),金槍魚資源量雖然低于塞舌爾群島海域,但調(diào)查的CPUE值要高于后者。因?yàn)樵隈R克薩斯群島海域大眼金槍魚垂直活動水層在250m以淺,和調(diào)查投鉤深度吻合,因此大量的商業(yè)延繩釣漁船在此作業(yè)[1],符合通常定義的中心漁場標(biāo)準(zhǔn)。因此,分析次表層環(huán)境對延繩釣漁獲率影響是非常有意義的。

    Hinton和Nakano[26]考慮到漁船的作業(yè)效率與金槍魚在水團(tuán)中分布的交叉作用對捕撈作業(yè)有影響,將海洋環(huán)境對黃鰭金槍魚的行為限制加入到CPUE標(biāo)準(zhǔn)化中,提出了基于生境空間的CPUE標(biāo)準(zhǔn)化方法,其中1個(gè)關(guān)鍵的變量就是黃鰭金槍魚在每個(gè)溫度差內(nèi)分布的時(shí)間。研究證實(shí),結(jié)合金槍魚垂直習(xí)性選擇和金槍魚上鉤深度,能有效地降低相對豐度估計(jì)時(shí)的不確定性[26]。

    GAM模型常被用于CPUE標(biāo)準(zhǔn)化[24]。Lan等[12]采用海洋表層環(huán)境變量構(gòu)建了GAM模型,模型偏差解釋率為33.8%。本文采用GAM模型構(gòu)建了次表層環(huán)境變量對CPUE的非線性關(guān)系,模型偏差解釋率為28.1%。模型結(jié)果表明,次表層環(huán)境變量對CPUE有顯著影響,也證實(shí)在漁業(yè)CPUE標(biāo)準(zhǔn)化中,應(yīng)該考慮和納入影響金槍魚垂直分布的次表層環(huán)境變量,使標(biāo)準(zhǔn)化CPUE更加合理。

    5 不足和展望

    GAM模型的結(jié)果揭示了各個(gè)變量對漁獲率的影響。水溫是影響黃鰭金槍魚水平和垂直分布的關(guān)鍵環(huán)境因子,但其他環(huán)境因子也有重要影響,如溶解氧也會影響黃鰭金槍魚垂直分布,葉綠素會影響?zhàn)D料生物分布,從而影響黃鰭金槍魚的空間分布。由于沒有足夠的溶解氧數(shù)據(jù),本文沒有對其做相關(guān)分析。餌料生物分布對黃鰭金槍魚垂直分布的影響也有重要作用,這需要未來做更多的調(diào)查,獲取更多的數(shù)據(jù),多角度結(jié)合研究。

    猜你喜歡
    黃鰭溫躍層環(huán)境變量
    基于最大熵模型的云南思茅松潛在分布區(qū)
    黃鰭鯛選育新品系“中珠1 號”養(yǎng)殖試驗(yàn)
    黃鰭金槍魚深水網(wǎng)箱養(yǎng)殖技術(shù)初探
    不同規(guī)格黃鰭鯛的形態(tài)性狀與體質(zhì)量的相關(guān)性研究
    多AUV溫躍層觀測方法研究
    全球變暖背景下赤道太平洋溫躍層的快慢變化特征與機(jī)制*
    黃鰭鯛苗種淡化技術(shù)研究
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    徹底弄懂Windows 10環(huán)境變量
    熱帶太平洋溫躍層深度的年代際變化特征及原因*
    久久精品国产综合久久久| 99香蕉大伊视频| av网站在线播放免费| 久久午夜综合久久蜜桃| 9热在线视频观看99| 国产欧美日韩一区二区三区在线| 亚洲成av片中文字幕在线观看| 视频区欧美日本亚洲| 母亲3免费完整高清在线观看| 日本撒尿小便嘘嘘汇集6| 正在播放国产对白刺激| 亚洲精华国产精华精| 操美女的视频在线观看| 91精品国产国语对白视频| 欧美日韩av久久| 国产av在哪里看| 18禁国产床啪视频网站| 欧美在线一区亚洲| 精品日产1卡2卡| 在线国产一区二区在线| 大型黄色视频在线免费观看| 国产免费现黄频在线看| 国产成人一区二区三区免费视频网站| 丁香六月欧美| 成人特级黄色片久久久久久久| 国产野战对白在线观看| 久久久精品欧美日韩精品| 久久久久九九精品影院| 午夜精品久久久久久毛片777| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 天堂俺去俺来也www色官网| 亚洲精品国产区一区二| 欧美亚洲日本最大视频资源| 自线自在国产av| 少妇粗大呻吟视频| 国产精品亚洲一级av第二区| 久久香蕉精品热| 老司机靠b影院| 身体一侧抽搐| 91字幕亚洲| www日本在线高清视频| 欧美精品亚洲一区二区| 亚洲精品中文字幕在线视频| 国产一区在线观看成人免费| 亚洲第一青青草原| 国产精品久久视频播放| 久久人妻av系列| 九色亚洲精品在线播放| 国产精品国产av在线观看| 啪啪无遮挡十八禁网站| 黑人猛操日本美女一级片| 国产精品乱码一区二三区的特点 | 亚洲欧美日韩无卡精品| 精品免费久久久久久久清纯| 成年女人毛片免费观看观看9| 久久久久久久久中文| 三上悠亚av全集在线观看| 天堂动漫精品| 操出白浆在线播放| 日日摸夜夜添夜夜添小说| 国产99白浆流出| 1024香蕉在线观看| 欧美精品一区二区免费开放| 日韩大尺度精品在线看网址 | 久热这里只有精品99| 一夜夜www| 亚洲免费av在线视频| 看片在线看免费视频| 搡老熟女国产l中国老女人| 大香蕉久久成人网| 亚洲精品在线观看二区| 亚洲伊人色综图| 男人舔女人的私密视频| a级毛片黄视频| 亚洲精品在线美女| 婷婷精品国产亚洲av在线| 午夜a级毛片| 精品久久久久久久毛片微露脸| 国产精品久久久久久人妻精品电影| 午夜精品在线福利| 精品久久久久久电影网| 一区在线观看完整版| 午夜精品在线福利| 夜夜夜夜夜久久久久| 国产亚洲欧美精品永久| 国产成人av激情在线播放| 热re99久久国产66热| 亚洲色图av天堂| 国产成人欧美| 国产野战对白在线观看| 国产一卡二卡三卡精品| 国产精品九九99| 欧美成人性av电影在线观看| 黄片大片在线免费观看| 操美女的视频在线观看| 国产高清国产精品国产三级| 国产日韩一区二区三区精品不卡| 午夜精品国产一区二区电影| 欧美日本亚洲视频在线播放| 99热只有精品国产| 亚洲av第一区精品v没综合| 成人黄色视频免费在线看| 欧美另类亚洲清纯唯美| 欧美国产精品va在线观看不卡| av天堂在线播放| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美一区二区综合| 欧美黄色淫秽网站| 久久影院123| 久久伊人香网站| 亚洲 欧美一区二区三区| 超碰成人久久| 身体一侧抽搐| 99re在线观看精品视频| 成年人免费黄色播放视频| 日本a在线网址| 亚洲色图 男人天堂 中文字幕| 亚洲成a人片在线一区二区| 伦理电影免费视频| netflix在线观看网站| 亚洲欧美激情在线| 成人国产一区最新在线观看| 色噜噜av男人的天堂激情| 欧美黄色片欧美黄色片| 99热6这里只有精品| 看黄色毛片网站| 国产男靠女视频免费网站| 18禁在线播放成人免费| 色尼玛亚洲综合影院| 夜夜看夜夜爽夜夜摸| 欧美黑人欧美精品刺激| 久久久国产成人免费| 99久久精品一区二区三区| 嫩草影视91久久| 偷拍熟女少妇极品色| 欧美成人免费av一区二区三区| 欧美激情久久久久久爽电影| 99热只有精品国产| 免费在线观看成人毛片| 三级男女做爰猛烈吃奶摸视频| 国产91精品成人一区二区三区| 男人和女人高潮做爰伦理| 丁香六月欧美| 少妇人妻一区二区三区视频| 一本久久中文字幕| 人人妻人人看人人澡| 午夜免费成人在线视频| www.www免费av| 99热精品在线国产| 婷婷色综合大香蕉| 精品欧美国产一区二区三| 免费一级毛片在线播放高清视频| 一级a爱片免费观看的视频| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添小说| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 国产 在线| 国产精品女同一区二区软件 | 亚洲在线观看片| 日韩精品青青久久久久久| 欧美性猛交╳xxx乱大交人| 免费看日本二区| 亚洲av五月六月丁香网| 丰满人妻一区二区三区视频av| 欧美色视频一区免费| 真人做人爱边吃奶动态| 国产免费av片在线观看野外av| 久久久色成人| 久久久久国产精品人妻aⅴ院| 黄色女人牲交| 热99re8久久精品国产| 日韩国内少妇激情av| 久久这里只有精品中国| 99热6这里只有精品| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| 97超级碰碰碰精品色视频在线观看| netflix在线观看网站| 在线观看一区二区三区| 三级毛片av免费| 好看av亚洲va欧美ⅴa在| 国产高清激情床上av| 男人和女人高潮做爰伦理| netflix在线观看网站| 日本 av在线| 国产精华一区二区三区| 日韩国内少妇激情av| 我的女老师完整版在线观看| 色视频www国产| 看十八女毛片水多多多| 九九热线精品视视频播放| 久久久国产成人精品二区| av天堂中文字幕网| 哪里可以看免费的av片| 亚洲美女搞黄在线观看 | 搞女人的毛片| 欧美精品国产亚洲| 99在线视频只有这里精品首页| 国产久久久一区二区三区| 深爱激情五月婷婷| 国产一区二区三区在线臀色熟女| 欧美日韩乱码在线| 禁无遮挡网站| 欧美极品一区二区三区四区| 亚洲,欧美,日韩| 自拍偷自拍亚洲精品老妇| 一二三四社区在线视频社区8| 99国产综合亚洲精品| 一边摸一边抽搐一进一小说| 一个人观看的视频www高清免费观看| 小蜜桃在线观看免费完整版高清| 欧美乱色亚洲激情| 在线天堂最新版资源| 18禁黄网站禁片免费观看直播| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久com| 99久久精品国产亚洲精品| 看片在线看免费视频| 精品人妻偷拍中文字幕| 亚洲18禁久久av| 亚洲欧美清纯卡通| av在线老鸭窝| 久久精品91蜜桃| 国产私拍福利视频在线观看| 欧美成人一区二区免费高清观看| 国产精品久久久久久久久免 | 日本成人三级电影网站| 一a级毛片在线观看| 岛国在线免费视频观看| 99久久精品国产亚洲精品| 国产精品99久久久久久久久| 午夜福利18| 听说在线观看完整版免费高清| 亚洲avbb在线观看| av在线观看视频网站免费| 亚洲精品影视一区二区三区av| 免费观看精品视频网站| 免费av毛片视频| 欧美色欧美亚洲另类二区| 中文亚洲av片在线观看爽| 舔av片在线| av黄色大香蕉| 在线观看一区二区三区| 亚洲一区高清亚洲精品| 老熟妇仑乱视频hdxx| 日日干狠狠操夜夜爽| 亚洲在线自拍视频| 亚洲综合色惰| 亚洲人成网站在线播| 一进一出抽搐gif免费好疼| 亚洲av五月六月丁香网| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| www.熟女人妻精品国产| 亚洲成人中文字幕在线播放| 国产美女午夜福利| 身体一侧抽搐| 桃红色精品国产亚洲av| 欧美最新免费一区二区三区 | 成人美女网站在线观看视频| 亚洲精品日韩av片在线观看| 少妇人妻精品综合一区二区 | 搡老妇女老女人老熟妇| 亚洲午夜理论影院| 国产亚洲欧美在线一区二区| 亚洲人成网站高清观看| 婷婷六月久久综合丁香| 美女高潮的动态| 精品一区二区三区视频在线观看免费| 男女视频在线观看网站免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇中文字幕五十中出| 国产久久久一区二区三区| 亚洲真实伦在线观看| 欧美成人一区二区免费高清观看| 成人国产综合亚洲| 天美传媒精品一区二区| 99精品在免费线老司机午夜| 国产私拍福利视频在线观看| 久久精品国产亚洲av天美| 日本与韩国留学比较| 久9热在线精品视频| 久久精品国产亚洲av涩爱 | 国产精品99久久久久久久久| 亚洲av电影在线进入| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久久末码| 特级一级黄色大片| 亚洲最大成人av| 天堂网av新在线| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 婷婷精品国产亚洲av在线| 午夜福利视频1000在线观看| 欧美最新免费一区二区三区 | 免费一级毛片在线播放高清视频| 亚洲熟妇中文字幕五十中出| 在线播放国产精品三级| 婷婷六月久久综合丁香| 99热这里只有精品一区| eeuss影院久久| 欧美色欧美亚洲另类二区| 每晚都被弄得嗷嗷叫到高潮| 久久久色成人| 亚洲第一区二区三区不卡| 成人午夜高清在线视频| 国产亚洲精品综合一区在线观看| 亚洲av日韩精品久久久久久密| 最近最新中文字幕大全电影3| 日韩成人在线观看一区二区三区| 麻豆av噜噜一区二区三区| 身体一侧抽搐| 精品欧美国产一区二区三| 麻豆成人av在线观看| 久久6这里有精品| 小蜜桃在线观看免费完整版高清| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 久久婷婷人人爽人人干人人爱| 两人在一起打扑克的视频| 久久热精品热| 18禁黄网站禁片午夜丰满| 熟女电影av网| 亚洲五月婷婷丁香| 久久精品91蜜桃| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 欧美日韩中文字幕国产精品一区二区三区| 人妻丰满熟妇av一区二区三区| 国产私拍福利视频在线观看| 老女人水多毛片| 久久精品国产亚洲av香蕉五月| 国产v大片淫在线免费观看| 亚洲精品一区av在线观看| 成人性生交大片免费视频hd| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 国产精品1区2区在线观看.| 日本免费a在线| 国产在线男女| 国产三级黄色录像| 窝窝影院91人妻| www.999成人在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品人妻熟女av久视频| 在线观看av片永久免费下载| 三级毛片av免费| 免费在线观看影片大全网站| 欧美日韩福利视频一区二区| 久久久久久九九精品二区国产| 欧美成人性av电影在线观看| 国产aⅴ精品一区二区三区波| 1000部很黄的大片| 日韩欧美三级三区| 精品熟女少妇八av免费久了| 亚洲va日本ⅴa欧美va伊人久久| 丰满人妻一区二区三区视频av| 很黄的视频免费| 国产黄片美女视频| 精品久久国产蜜桃| 美女高潮的动态| 国产亚洲精品久久久久久毛片| 久久久久久久午夜电影| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| 国产一区二区亚洲精品在线观看| 午夜激情福利司机影院| 少妇熟女aⅴ在线视频| 久久中文看片网| 亚洲欧美清纯卡通| 亚洲美女搞黄在线观看 | 2021天堂中文幕一二区在线观| 国产黄a三级三级三级人| 国产亚洲av嫩草精品影院| 内地一区二区视频在线| 免费看日本二区| a在线观看视频网站| 国产成人aa在线观看| 欧洲精品卡2卡3卡4卡5卡区| 成人高潮视频无遮挡免费网站| 最好的美女福利视频网| 国产亚洲欧美在线一区二区| 国产亚洲精品久久久com| 欧美潮喷喷水| 国产三级中文精品| 精品乱码久久久久久99久播| 成人无遮挡网站| 一级av片app| 欧美+日韩+精品| 精品一区二区免费观看| 热99re8久久精品国产| 可以在线观看的亚洲视频| 老司机福利观看| 国产精品1区2区在线观看.| 日韩欧美精品v在线| 最近最新免费中文字幕在线| 国产探花在线观看一区二区| 欧美在线黄色| 亚洲精品粉嫩美女一区| 国产亚洲欧美98| 熟女电影av网| av在线天堂中文字幕| 一本精品99久久精品77| 国产三级在线视频| 国模一区二区三区四区视频| 亚洲在线观看片| 精品人妻视频免费看| 精品欧美国产一区二区三| 国产探花极品一区二区| 男女床上黄色一级片免费看| 又爽又黄无遮挡网站| 999久久久精品免费观看国产| 亚洲av第一区精品v没综合| 免费在线观看日本一区| 精品久久久久久久久亚洲 | 一个人免费在线观看电影| 亚洲第一电影网av| 欧美+亚洲+日韩+国产| 午夜福利在线观看吧| 亚洲最大成人中文| 99在线人妻在线中文字幕| 亚洲av不卡在线观看| 全区人妻精品视频| 午夜a级毛片| 男人舔女人下体高潮全视频| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久成人| 变态另类丝袜制服| 丰满人妻熟妇乱又伦精品不卡| 97碰自拍视频| 一个人观看的视频www高清免费观看| 国产精品电影一区二区三区| 五月伊人婷婷丁香| 麻豆国产av国片精品| 国产精品乱码一区二三区的特点| 国产国拍精品亚洲av在线观看| 日本在线视频免费播放| 精品人妻1区二区| 久久久久九九精品影院| 国内精品一区二区在线观看| 亚洲熟妇熟女久久| 丰满人妻熟妇乱又伦精品不卡| 国产野战对白在线观看| 久久久国产成人免费| 精品99又大又爽又粗少妇毛片 | 欧美bdsm另类| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 最新在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 欧美极品一区二区三区四区| 色吧在线观看| 国产精品野战在线观看| 免费观看人在逋| 少妇高潮的动态图| 在线观看午夜福利视频| 不卡一级毛片| 国产中年淑女户外野战色| 午夜免费成人在线视频| 特大巨黑吊av在线直播| 九九久久精品国产亚洲av麻豆| 国产单亲对白刺激| 欧美一区二区精品小视频在线| 久久婷婷人人爽人人干人人爱| 国产精品精品国产色婷婷| 欧洲精品卡2卡3卡4卡5卡区| 99久久久亚洲精品蜜臀av| 免费搜索国产男女视频| 欧美最黄视频在线播放免费| 国产高清激情床上av| 日韩有码中文字幕| 全区人妻精品视频| 日韩av在线大香蕉| 国产精品人妻久久久久久| 亚洲激情在线av| 老司机福利观看| 精品人妻1区二区| 别揉我奶头~嗯~啊~动态视频| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 国产视频一区二区在线看| 成人午夜高清在线视频| 欧美日本视频| 久久精品影院6| 中文字幕人妻熟人妻熟丝袜美| 五月伊人婷婷丁香| 国产视频内射| 国内精品美女久久久久久| 欧美成人一区二区免费高清观看| 免费看光身美女| 18禁在线播放成人免费| a级一级毛片免费在线观看| 香蕉av资源在线| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 99精品久久久久人妻精品| 国产午夜福利久久久久久| 性色av乱码一区二区三区2| 极品教师在线免费播放| 亚洲在线观看片| 一区福利在线观看| 在线看三级毛片| 小说图片视频综合网站| 欧美在线一区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 88av欧美| 国产aⅴ精品一区二区三区波| 婷婷六月久久综合丁香| 国产v大片淫在线免费观看| 亚洲真实伦在线观看| 国产成人av教育| 精品一区二区三区人妻视频| 国产精品久久久久久人妻精品电影| 伊人久久精品亚洲午夜| 中文字幕熟女人妻在线| 中文字幕高清在线视频| 国产精品日韩av在线免费观看| 国产av麻豆久久久久久久| 99在线人妻在线中文字幕| 窝窝影院91人妻| 国产精品三级大全| 在线播放无遮挡| 欧美成人一区二区免费高清观看| 欧美+亚洲+日韩+国产| 欧美高清性xxxxhd video| 国产欧美日韩精品亚洲av| 国产伦在线观看视频一区| 97热精品久久久久久| 亚洲专区中文字幕在线| 日本三级黄在线观看| 精品久久久久久久久av| 亚洲最大成人中文| 午夜精品久久久久久毛片777| 五月伊人婷婷丁香| 99精品在免费线老司机午夜| 亚洲aⅴ乱码一区二区在线播放| 99久久久亚洲精品蜜臀av| 色综合欧美亚洲国产小说| 国产成年人精品一区二区| 搞女人的毛片| 1024手机看黄色片| 少妇被粗大猛烈的视频| 国产av麻豆久久久久久久| 免费搜索国产男女视频| 国产三级中文精品| 天堂影院成人在线观看| 午夜福利免费观看在线| 欧美区成人在线视频| 99热这里只有是精品在线观看 | 人妻久久中文字幕网| 亚洲国产欧洲综合997久久,| 亚洲国产精品合色在线| 熟女人妻精品中文字幕| 亚洲美女视频黄频| 小蜜桃在线观看免费完整版高清| 欧美成人一区二区免费高清观看| 最好的美女福利视频网| 淫妇啪啪啪对白视频| 国产午夜精品论理片| 久久九九热精品免费| 久久人人精品亚洲av| 天堂网av新在线| 亚洲欧美日韩卡通动漫| 麻豆成人av在线观看| 99热这里只有是精品在线观看 | 91九色精品人成在线观看| www.www免费av| 国产高潮美女av| 色综合欧美亚洲国产小说| 女生性感内裤真人,穿戴方法视频| 国产精品98久久久久久宅男小说| 狂野欧美白嫩少妇大欣赏| 男女做爰动态图高潮gif福利片| 免费av观看视频| 最后的刺客免费高清国语| 亚洲自拍偷在线| 亚洲在线观看片| 亚洲黑人精品在线| 午夜免费成人在线视频| 淫妇啪啪啪对白视频| 亚洲国产精品成人综合色| 在线观看av片永久免费下载| 久久久久精品国产欧美久久久| 免费看美女性在线毛片视频| 免费在线观看影片大全网站| 日本三级黄在线观看| 俄罗斯特黄特色一大片| 国产精品一区二区性色av| 欧美高清性xxxxhd video| 90打野战视频偷拍视频| 九九在线视频观看精品| 成年女人毛片免费观看观看9| 免费高清视频大片| 黄色日韩在线| 在线看三级毛片| 十八禁网站免费在线| 日韩精品中文字幕看吧| x7x7x7水蜜桃| 亚洲,欧美精品.| 成人国产综合亚洲| 国产日本99.免费观看| 五月玫瑰六月丁香| 网址你懂的国产日韩在线| 成年免费大片在线观看| 床上黄色一级片| av福利片在线观看| 免费看日本二区| a级毛片免费高清观看在线播放| 亚洲电影在线观看av| 精品久久久久久成人av|