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

    秘魯外海莖柔魚索餌場棲息地研究

    2023-05-18 05:00:10劉必林陳新軍田思泉
    水產(chǎn)科學(xué) 2023年3期
    關(guān)鍵詞:漁獲量加性表面溫度

    劉必林,陳新軍,田思泉

    ( 1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海 201306; 2.國家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心,上海 201306; 3.大洋漁業(yè)資源可持續(xù)教育部重點(diǎn)實(shí)驗(yàn)室,上海 201306; 4.農(nóng)業(yè)農(nóng)村部大洋漁業(yè)開發(fā)重點(diǎn)實(shí)驗(yàn)室,上海 201306; 5.農(nóng)業(yè)農(nóng)村部大洋漁業(yè)資源環(huán)境科學(xué)觀測實(shí)驗(yàn)站,上海 201306 )

    莖柔魚(Dosidicusgigas)為大洋性淺海種,廣泛分布于東太平洋北美洲北部N 40°至智利南部S 47°,是世界上最主要的頭足類資源之一,在我國遠(yuǎn)洋魷釣漁業(yè)中占據(jù)著極為重要的地位,其年產(chǎn)量約占我國遠(yuǎn)洋魷釣產(chǎn)量的1/3以上。我國于2001年6—9月在秘魯外海對莖柔魚資源進(jìn)行了首次生產(chǎn)性調(diào)查,并取得了成功,最高年產(chǎn)量突破20萬t。廣義加性模型最先由Hastie等[1]提出,之后被廣泛應(yīng)用到漁業(yè)資源的研究中[2-9]。兩步廣義加性模型作為廣義加性模型的擴(kuò)展[10-11],能夠在建模過程中保留漁業(yè)生產(chǎn)中產(chǎn)量為零的數(shù)據(jù),因此近年來被廣泛用于預(yù)測漁業(yè)資源的棲息地分布。

    目前有關(guān)莖柔魚的棲息地研究主要采用棲息地適宜性指數(shù)模型[12-16],而該模型無法解決產(chǎn)量為零的數(shù)據(jù)。筆者根據(jù)2006—2010年秘魯外海莖柔魚魷釣生產(chǎn)數(shù)據(jù)以及海洋環(huán)境數(shù)據(jù),利用兩步廣義加性模型分析莖柔魚適宜索餌場棲息環(huán)境,預(yù)測其索餌場棲息地并進(jìn)行驗(yàn)證。

    1 材料與方法

    1.1 漁業(yè)生產(chǎn)數(shù)據(jù)獲取

    2006—2010年東南太平洋W 74°~86°,S 8°~30°海域莖柔魚魷釣生產(chǎn)數(shù)據(jù)由上海海洋大學(xué)魷釣技術(shù)組提供,包括生產(chǎn)日期(每年1—12月)、經(jīng)度、緯度、產(chǎn)量、作業(yè)船數(shù)(年均200余艘)和平均日產(chǎn)量。其時(shí)間分辨率為d,空間分辨率為0.5°×0.5°。單位捕撈努力量漁獲量(CPUE)計(jì)算公式如下:

    (1)

    式中,∑m為0.5°×0.5°內(nèi)所有漁船的產(chǎn)量(t),∑t為0.5°×0.5°范圍內(nèi)所有漁船的總作業(yè)時(shí)間(d)。

    1.2 環(huán)境數(shù)據(jù)獲取

    莖柔魚洄游特性明顯,其資源量會發(fā)生季節(jié)性變化,從而呈現(xiàn)季節(jié)性分布,因此筆者按春、夏、秋、冬4個(gè)季節(jié)研究其分布與主要環(huán)境因子的關(guān)系。海表面溫度和海表面高度是影響莖柔魚空間分布的重要指標(biāo)[12-16]。此外,葉綠素a代表初級生產(chǎn)力,其含量高低直接影響莖柔魚幼魚的豐度,進(jìn)而影響成魚的補(bǔ)充量。環(huán)境數(shù)據(jù)取自網(wǎng)站http://oceanwatch.pifsc.noaa.gov/las/servlets/dataset,其中:海表面溫度時(shí)間分辨率為周,空間分辨率為0.1°×0.1°;海表面高度時(shí)間分辨率為周,空間分辨率為0.3°×0.3°;葉綠素a時(shí)間分辨率為周,空間分辨率為0.05°×0.05°。為了與生產(chǎn)數(shù)據(jù)相匹配,環(huán)境數(shù)據(jù)按年、月、日的時(shí)間分辨率及0.5°×0.5°的空間分辨率進(jìn)行整合。

    1.3 索餌場推測

    1.3.1 模型建立

    采用兩步廣義加性模型[10-11,17]預(yù)測和分析春、夏、秋、冬4個(gè)季節(jié)的莖柔魚的分布。廣義加性模型1利用二項(xiàng)式分布的邏輯連接方程估算莖柔魚存在的概率(P),廣義加性模型2利用高斯分布估算取對數(shù)后的莖柔魚單位捕撈努力量漁獲量(y),最后取對數(shù)的綜合單位捕撈努力量漁獲量(豐度)(D)等于廣義加性模型1和廣義加性模型2兩者所得概率相乘。廣義加性模型1(GAM1)和廣義加性模型2(GAM2)計(jì)算方程如下:

    GAM1: logit (P)=s(Lon)+s(Lat)+s(TSS)+s(hSS)+s(ρChla)+s(Inter)+ε

    (2)

    GAM2: ln (y)=s(Lon)+s(Lat)+s(TSS)+s(hSS)+s(ρChla)+s(Inter)+ε

    (3)

    ln (D)=p×ln (y)

    (4)

    式中,Lon為經(jīng)度(°),Lat為緯度(°),TSS為海表面溫度(℃),hSS為海表面高度(cm),ρChla為葉綠素a質(zhì)量濃度(mg/m3),Inter為交互項(xiàng),s為平滑函數(shù)。

    對樣本經(jīng)度、緯度、單位捕撈努力量漁獲量、海表面溫度、海表面高度、葉綠素a等6個(gè)因子進(jìn)行相關(guān)性分析,結(jié)果顯示,經(jīng)度與緯度相關(guān)性達(dá)到0.82,因此兩者當(dāng)中只選擇緯度放入模型中計(jì)算。

    廣義加性模型估算采用R軟件mgcv程序包。運(yùn)用卡方檢驗(yàn)和F統(tǒng)計(jì)檢驗(yàn)對廣義加性模型1和廣義加性模型2的主效應(yīng)因子及交互效應(yīng)因子進(jìn)行檢驗(yàn)。

    1.3.2 模型驗(yàn)證

    兩步廣義加性模型估算的結(jié)果用來預(yù)測各捕撈地點(diǎn)的莖柔魚的豐度。預(yù)測的莖柔魚豐度(D′)與觀察豐度(D)回歸方程如下:

    ln (D)=a+b×ln (D′)

    (5)

    式中,系數(shù)a代表系統(tǒng)誤差,系數(shù)b越接近1表示預(yù)測值與觀察越接近,預(yù)測越準(zhǔn)確。系數(shù)檢驗(yàn)采用雙尾檢驗(yàn),H0:a=0,H0:b=1。

    1.3.3 模型預(yù)測

    為了預(yù)測適合的索餌場,將生產(chǎn)海區(qū)的環(huán)境數(shù)據(jù)代入所建立的兩步廣義加性模型,以此生成對應(yīng)的豐度,然后利用ArcGIS 9.0軟件繪圖。

    2 結(jié) 果

    2.1 影響因子的選擇

    春、夏、秋、冬廣義加性模型1解釋莖柔魚存在概率的離差為24.1%~42.4%,廣義加性模型2解釋莖柔魚單位捕撈努力量漁獲量的離差為23.3%~41.7%。春季,用于廣義加性模型1分析的所有空間和環(huán)境因子均不顯著,用于廣義加性模型2分析的緯度因子顯著(表1);夏季,用于廣義加性模型1分析的海表面溫度和單位捕撈努力量漁獲量因子顯著,用于廣義加性模型2分析的緯度、海表面溫度和葉綠素a顯著;秋季,用于廣義加性模型1分析的緯度和海表面溫度因子顯著,用于廣義加性模型2分析的海表面溫度、葉綠素a和海表面高度因子顯著;冬季,用于廣義加性模型1分析的緯度因子顯著,用于廣義加性模型2分析的緯度和海表面高度因子顯著。

    表1 莖柔魚4個(gè)季節(jié)存在概率(廣義加性模型1)和單位捕撈努力量漁獲量(廣義加性模型2)分析結(jié)果Tab.1 Summary results of the ‘GAM1=probability of presence model’ and ‘GAM2=CPUE model’ for Jumbo flying squid D. gigas in the four modeling seasons

    2.2 模型擬合診斷

    模型擬合診斷結(jié)果顯示,殘差的數(shù)據(jù)點(diǎn)分布在正態(tài)Q-Q圖上基本呈線性重合(圖1),散點(diǎn)分布圖比較均勻,表明單位捕撈努力量漁獲量數(shù)據(jù)適合用廣義加性模型分析。

    圖1 春季(a)、夏季(b)、秋季(c)、冬季(d)莖柔魚單位捕撈努力量漁獲量擬合優(yōu)度的Q-Q診斷圖Fig.1 Q-Q plots for the goodness-of-fit GAM model to the CPUE data of jumbo flying squid D. gigas for spring(a), summer(b),autumn(c) and winter(d)

    2.3 索餌場棲息環(huán)境

    春季,空間分布對單位捕撈努力量漁獲量影響表明,主要作業(yè)漁場位于S 11°~15°,隨著緯度降低,對單位捕撈努力量漁獲量的影響明顯下降(P<0.001,表1、圖2a)。從單位捕撈努力量漁獲量與海表面溫度關(guān)系可以看出,作業(yè)海表面溫度為16.5~20.5 ℃,海表面溫度對單位捕撈努力量漁獲量的影響略呈上升趨勢(P>0.05,表1、圖2b)。單位捕撈努力量漁獲量與葉綠素a質(zhì)量濃度(0.2~0.4 mg/m3)關(guān)系顯示,葉綠素a質(zhì)量濃度對單位捕撈努力量漁獲量影響不顯著(P>0.05,表1、圖2c)。海表面高度對單位捕撈努力量漁獲量的影響顯示,主要海表面高度為20~30 cm,海表面高度對單位捕撈努力量漁獲量的影響不顯著(P>0.05,表1、圖2d)。

    圖2 春季索餌場各空間和環(huán)境要素對單位捕撈努力量漁獲量的相對效應(yīng)Fig.2 Relative effects of spatial and environmental factors on CPUE in feeding ground in spring虛線為95%置信區(qū)間;橫軸內(nèi)側(cè)刻度表示數(shù)據(jù)點(diǎn)相對密度;下同.Dashed lines indicate 95% confidence intervals; the relative density of data points are shown by the rug on the abscissa; et sequentia.

    夏季,空間分布對單位捕撈努力量漁獲量影響表明,主要作業(yè)漁場位于S 12°~17°,隨著緯度的降低,對單位捕撈努力量漁獲量的影響逐漸減小(P<0.001,表1、圖3a)。從單位捕撈努力量漁獲量與海表面溫度關(guān)系可以看出,主要作業(yè)海表面溫度為20.5~25.0 ℃,海表面溫度在20.5~21.5 ℃時(shí)對單位捕撈努力量漁獲量的影響顯著(P<0.001,表1、圖3b)。單位捕撈努力量漁獲量與葉綠素a關(guān)系顯示,葉綠素a質(zhì)量濃度在0.2~0.3 mg/m3對單位捕撈努力量漁獲量的影響呈升高趨勢(P<0.05,表1、圖3c)。海表面高度對單位捕撈努力量漁獲量的影響顯示,主要海表面高度為24~31 cm,海表面高度對單位捕撈努力量漁獲量的影響不顯著(P>0.05,表1、圖3d)。

    圖3 夏季索餌場各空間和環(huán)境要素對單位捕撈努力量漁獲量的相對效應(yīng)Fig.3 Relative effects of spatial and environmental factors on CPUE in feeding ground in summer

    秋季,空間分布對單位捕撈努力量漁獲量影響表明,作業(yè)漁場范圍擴(kuò)大到S 12°~30°,雖然隨著緯度的降低,對單位捕撈努力量漁獲量的影響逐漸下降,但是由于95%置信區(qū)間過大而存在較大的不確定性(P>0.05,表1、圖4a)。從單位捕撈努力量漁獲量與海表面溫度關(guān)系可以看出,主要作業(yè)海表面溫度為20.5~26.0 ℃,海表面溫度對單位捕撈努力量漁獲量的影響明顯逐漸下降(P<0.001,表1、圖4b)。

    圖4 秋季索餌場各空間和環(huán)境要素對單位捕撈努力量漁獲量的相對效應(yīng)Fig.4 Relative effects of spatial and environmental factors on CPUE in feeding ground in autumn

    單位捕撈努力量漁獲量與葉綠素a關(guān)系顯示,葉綠素a質(zhì)量濃度對單位捕撈努力量漁獲量的影響逐漸升高(P<0.05,表1、圖4c)。海表面高度對單位捕撈努力量漁獲量的影響顯示,主要海表面高度為20~32 cm,海表面高度在28~32 cm時(shí)對單位捕撈努力量漁獲量的影響最大(P<0.001,表1、圖4d)。

    冬季,空間分布對單位捕撈努力量漁獲量影響表明,雖然主要作業(yè)漁場位于S 12°~18°,隨著緯度的降低,對單位捕撈努力量漁獲量的影響呈明顯升高的趨勢(P<0.001,表1、圖5a)。由單位捕撈努力量漁獲量與海表面溫度關(guān)系可以看出,主要作業(yè)海表面溫度為17~20 ℃,海表面溫度對單位捕撈努力量漁獲量的影響不顯著(P>0.05,表1、圖5b)。單位捕撈努力量漁獲量與葉綠素a質(zhì)量濃度關(guān)系顯示,葉綠素a質(zhì)量濃度對單位捕撈努力量漁獲量的影響不顯著(P>0.05,表1、圖5c)。海表面高度對與單位捕撈努力量漁獲量關(guān)系顯示,10~26 cm對單位捕撈努力量漁獲量的影響呈下降趨勢,主要海表面高度為20~32 cm,對單位捕撈努力量漁獲量的影響呈拋物線狀,即海表面高度在20~28 cm段對單位捕撈努力量漁獲量的影響不顯著,在28~32 cm對單位捕撈努力量漁獲量的影響顯著(P<0.001,表1、圖5d)。

    圖5 冬季索餌場各空間和環(huán)境要素對單位捕撈努力量漁獲量的相對效應(yīng)Fig.5 Relative effects of spatial and environmental factors on CPUE in feeding ground in winter

    2.4 索餌場位置預(yù)測

    春季,大部分海域單位捕撈努力量漁獲量都很低,S 15°,W 77°附近海域單位捕撈努力量漁獲量較高,可能是莖柔魚適合的索餌場(圖6a);夏季,S 22°以北海域單位捕撈努力量漁獲量低,S 22°以南單位捕撈努力量漁獲量較高,S 24~26°,W 75~78°海域單位捕撈努力量漁獲量最高,是莖柔魚適合的索餌場(圖6b);秋季秘魯沿岸S 12~16°和智利中北部S 28°附近海域單位捕撈努力量漁獲量較高,是莖柔魚的適合索餌場(圖6c);冬季秘魯沿岸高單位捕撈努力量漁獲量海區(qū)基本不變,而智利沿岸則消失(圖6d)。

    圖6 預(yù)測2006—2010年春季(a)、夏季(b)、秋季(c)和冬季(d)秘魯莖柔魚單位捕撈努力量漁獲量分布Fig.6 Predicted CPUE distributions of Jumbo flying squid D. gigas in spring (a), summer (b), autumn (c) and winter (d) from 2006 to 2010

    2.5 模型驗(yàn)證

    單位捕撈努力量漁獲量觀察值與預(yù)測值回歸分析顯示,截距接近0(P<0.05),斜率接近1(P<0.05),說明模型預(yù)測效果較好(表2)。

    表2 莖柔魚4個(gè)季節(jié)觀察單位捕撈努力量漁獲量與預(yù)測單位捕撈努力量漁獲量回歸分析結(jié)果Tab.2 Summary results of the regression analysis of the observed CPUE and the predicted CPUE for the four modeling seasons

    3 討 論

    3.1 廣義加性模型優(yōu)越性

    近年來,統(tǒng)計(jì)回歸模型被廣泛應(yīng)用于魚類的棲息地預(yù)測及與環(huán)境關(guān)系的分析研究中[11,17-20]。廣義加性模型最早由Hastie等[1]提出,是最常見的回歸模型[3]。國外學(xué)者較早地將廣義加性模型應(yīng)用到漁業(yè)資源的研究中[3-6]。近年來,國內(nèi)學(xué)者也開始將廣義加性模型應(yīng)用到漁業(yè)資源與環(huán)境關(guān)系的研究中[7-9]。研究認(rèn)為,廣義加性模型因其靈活性,要比其他線性和非線性模型更適合用來研究魚類分布與海洋環(huán)境的關(guān)系[21-22]。兩步廣義加性模型作為廣義加性模型的擴(kuò)展,它分兩步對漁業(yè)資源進(jìn)行分析,首先分析哪里有或者哪里沒有漁業(yè)資源的分布,其次再分析有資源分布的地方的分布密度是多少[10-11]。該方法的優(yōu)點(diǎn)在于能夠在建模過程中保留漁業(yè)生產(chǎn)中產(chǎn)量為零的數(shù)據(jù)。

    3.2 單位捕撈努力量漁獲量與環(huán)境因子關(guān)系

    莖柔魚為大洋性廣溫種,適宜生存溫度15~28 ℃[23],所能忍受的極限溫度下限為4 ℃,上限為32 ℃[24]。已有的研究表明,海表面溫度是影響莖柔魚資源分布的關(guān)鍵因子[12-15],北半球溫度24~29 ℃、南半球溫度17~23 ℃時(shí)莖柔魚資源密度最高[23,25-30]。夏季和秋季作業(yè)海表面溫度在20~26 ℃,廣義加性模型分析顯示,海表面溫度變化對莖柔魚資源密度影響明顯,21 ℃附近單位捕撈努力量漁獲量最高;冬季和春季作業(yè)海表面溫度主要在16~20 ℃,廣義加性模型分析顯示,海表面溫度變化對莖柔魚資源密度無明顯影響。葉綠素a含量的高低反應(yīng)初級生產(chǎn)力狀況,因此將其看作是指示莖柔魚產(chǎn)卵場的重要環(huán)境因子[31],然而對成魚的分析發(fā)現(xiàn),葉綠素a含量對單位捕撈努力量漁獲量的影響不是特別明顯,這是因?yàn)樽鳂I(yè)漁場通常位于索餌場,這一結(jié)論與胡振明等[12]在秘魯外海的研究結(jié)論一致。

    3.3 索餌場預(yù)測

    通過單位捕撈努力量漁獲量與海洋環(huán)境關(guān)系的分析,可以推測春季莖柔魚的索餌場位于S 15°,W 77°附近海域,而生產(chǎn)數(shù)據(jù)表明這一區(qū)域正是單位捕撈努力量漁獲量最高的海區(qū)。已有的研究認(rèn)為,莖柔魚在加利福尼亞至智利北部370~460 km海域內(nèi)的資源密度最為豐富[24],筆者推測莖柔魚秋季索餌場主要位于秘魯沿岸S 12°~16°和智利中北部沿岸S 28°附近,冬季位于秘魯沿岸S 13°~18°,而夏季索餌場位于智利北部的S 24°~26°公海海域。

    4 結(jié) 論

    春季莖柔魚索餌場位于S 15°,W 77°附近海域,夏季位于S 24°~26°,W 75°~78°海域,秋季位于秘魯沿岸S 12°~16°和智利中部S 28°附近海域,冬季位于秘魯沿岸S 13°~18°海域。莖柔魚屬高度洄游的大洋性頭足類,其資源、漁場的變動與海洋環(huán)境及氣候變化的關(guān)系復(fù)雜,尤其在東太平洋海域受厄爾尼諾/拉尼娜直接影響,而本研究中用于索餌場棲息地研究的環(huán)境數(shù)據(jù)只包括海表面溫度、海表面高度和葉綠素a,缺少鹽度、恩索指數(shù)等環(huán)境數(shù)據(jù),加上用于分析的漁業(yè)數(shù)據(jù)時(shí)間序列太短,因此筆者只是得到了初步的研究結(jié)果。今后將在獲得多年的漁業(yè)生產(chǎn)數(shù)據(jù)和足夠的環(huán)境數(shù)據(jù)的基礎(chǔ)上,進(jìn)一步系統(tǒng)深入研究。

    猜你喜歡
    漁獲量加性表面溫度
    東南太平洋公海魷魚釣漁業(yè)狀況分析
    ?2?4[u]-加性循環(huán)碼
    世界竹筴魚資源開發(fā)現(xiàn)狀及其建議
    結(jié)合注意力機(jī)制的區(qū)域型海表面溫度預(yù)報(bào)算法
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關(guān)系:超可加性、次可加性或不可加性
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關(guān)系:超可加性、次可加性或不可加性
    鄱陽湖湖口刀鱭幼魚出現(xiàn)的時(shí)間及漁獲量
    熱電池新型隔熱結(jié)構(gòu)設(shè)計(jì)及表面溫度研究
    基于加性指標(biāo)的網(wǎng)絡(luò)斷層掃描的研究
    Kerr型中子星與黑洞表面溫度分布的研究
    日韩欧美在线乱码| 国内精品久久久久久久电影| 88av欧美| 久久中文字幕人妻熟女| videosex国产| 午夜福利视频1000在线观看| 国产日本99.免费观看| 成人三级做爰电影| 成人18禁高潮啪啪吃奶动态图| 99精品在免费线老司机午夜| 成人国产综合亚洲| 免费无遮挡裸体视频| 欧美久久黑人一区二区| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| 国产精品一区二区三区四区久久| 搡老岳熟女国产| 熟妇人妻久久中文字幕3abv| 国产精品亚洲av一区麻豆| 亚洲九九香蕉| 国产精品,欧美在线| 亚洲五月天丁香| 精品一区二区三区视频在线观看免费| 久久精品aⅴ一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 国产片内射在线| 69av精品久久久久久| 中文字幕熟女人妻在线| 美女高潮喷水抽搐中文字幕| 啦啦啦免费观看视频1| 91麻豆精品激情在线观看国产| 久久久久久国产a免费观看| 久久精品国产亚洲av高清一级| 亚洲成av人片免费观看| 成熟少妇高潮喷水视频| 亚洲av成人一区二区三| 亚洲精品在线美女| 久久亚洲精品不卡| videosex国产| 亚洲欧美日韩高清专用| 欧美另类亚洲清纯唯美| 日韩欧美在线乱码| 国产高清视频在线播放一区| 久久精品夜夜夜夜夜久久蜜豆 | 欧美高清成人免费视频www| 亚洲av五月六月丁香网| 国产精品久久久久久亚洲av鲁大| 欧美黄色片欧美黄色片| 搡老妇女老女人老熟妇| 最近最新中文字幕大全免费视频| 97超级碰碰碰精品色视频在线观看| 久久伊人香网站| 精品久久久久久久久久免费视频| 黄色成人免费大全| 久久久久亚洲av毛片大全| 国产精品日韩av在线免费观看| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 久久久久亚洲av毛片大全| 91成年电影在线观看| 狂野欧美白嫩少妇大欣赏| 一级毛片精品| 2021天堂中文幕一二区在线观| 久久久精品欧美日韩精品| 亚洲欧美激情综合另类| 天天一区二区日本电影三级| 国产欧美日韩一区二区三| 这个男人来自地球电影免费观看| 亚洲欧美日韩高清在线视频| 99国产极品粉嫩在线观看| 午夜福利免费观看在线| 国产97色在线日韩免费| 激情在线观看视频在线高清| 国产精品一及| 一本一本综合久久| 制服人妻中文乱码| 免费无遮挡裸体视频| 婷婷精品国产亚洲av在线| 1024香蕉在线观看| а√天堂www在线а√下载| 亚洲一卡2卡3卡4卡5卡精品中文| 正在播放国产对白刺激| 精品一区二区三区av网在线观看| 欧美不卡视频在线免费观看 | 国产欧美日韩一区二区精品| 亚洲av电影不卡..在线观看| 91九色精品人成在线观看| 欧美三级亚洲精品| 婷婷六月久久综合丁香| 午夜亚洲福利在线播放| 欧美色视频一区免费| 91麻豆av在线| 午夜成年电影在线免费观看| 91在线观看av| 少妇粗大呻吟视频| 国产精品爽爽va在线观看网站| 禁无遮挡网站| 亚洲欧洲精品一区二区精品久久久| 中文资源天堂在线| 又粗又爽又猛毛片免费看| 亚洲,欧美精品.| 亚洲精品中文字幕一二三四区| 最近视频中文字幕2019在线8| 亚洲五月婷婷丁香| 日韩欧美国产在线观看| 99国产精品一区二区三区| 国产爱豆传媒在线观看 | 母亲3免费完整高清在线观看| 日日爽夜夜爽网站| 热99re8久久精品国产| 男人的好看免费观看在线视频 | 欧美黑人欧美精品刺激| 丁香欧美五月| 亚洲自拍偷在线| 日本黄大片高清| 久久久久亚洲av毛片大全| 变态另类丝袜制服| 亚洲色图 男人天堂 中文字幕| 在线播放国产精品三级| 午夜激情av网站| 国产精品久久久久久亚洲av鲁大| a级毛片a级免费在线| 精品欧美国产一区二区三| svipshipincom国产片| 国产伦在线观看视频一区| 天天躁夜夜躁狠狠躁躁| 欧美日韩中文字幕国产精品一区二区三区| 久久香蕉激情| xxx96com| 三级男女做爰猛烈吃奶摸视频| 久久国产乱子伦精品免费另类| 19禁男女啪啪无遮挡网站| 一区二区三区高清视频在线| 欧美在线黄色| 久久亚洲精品不卡| 国产成人精品无人区| 国产亚洲精品久久久久5区| 18禁黄网站禁片午夜丰满| xxx96com| ponron亚洲| 欧美又色又爽又黄视频| 757午夜福利合集在线观看| 久久午夜综合久久蜜桃| 欧美中文综合在线视频| x7x7x7水蜜桃| 18禁黄网站禁片午夜丰满| 国产视频内射| 日本撒尿小便嘘嘘汇集6| 精品久久久久久成人av| 免费无遮挡裸体视频| 亚洲成人中文字幕在线播放| 日韩欧美免费精品| 国产麻豆成人av免费视频| 免费在线观看日本一区| 一区二区三区高清视频在线| 国产av不卡久久| 十八禁人妻一区二区| 老司机福利观看| 国产伦一二天堂av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美一区二区三区黑人| 18禁黄网站禁片午夜丰满| 啦啦啦韩国在线观看视频| 亚洲中文字幕日韩| 亚洲精品美女久久久久99蜜臀| 亚洲精品久久国产高清桃花| 国产精品永久免费网站| 亚洲午夜精品一区,二区,三区| 国产精品久久久av美女十八| 美女午夜性视频免费| 亚洲成人久久性| 久热爱精品视频在线9| 婷婷精品国产亚洲av| 一区二区三区国产精品乱码| 亚洲av片天天在线观看| 中文在线观看免费www的网站 | 国产精品永久免费网站| 热99re8久久精品国产| 麻豆成人午夜福利视频| 最近视频中文字幕2019在线8| 在线免费观看的www视频| ponron亚洲| 欧美一区二区精品小视频在线| 在线a可以看的网站| 久久婷婷成人综合色麻豆| 此物有八面人人有两片| 成人三级黄色视频| 看片在线看免费视频| 精品第一国产精品| 在线观看免费日韩欧美大片| 日韩高清综合在线| 精华霜和精华液先用哪个| 精品国产亚洲在线| 亚洲一区中文字幕在线| 国产亚洲av高清不卡| 12—13女人毛片做爰片一| 老司机靠b影院| 2021天堂中文幕一二区在线观| 久久久水蜜桃国产精品网| 国产av又大| 免费在线观看日本一区| 亚洲 国产 在线| 久久久久久大精品| 亚洲人与动物交配视频| 在线观看舔阴道视频| 黑人巨大精品欧美一区二区mp4| 久久久久久人人人人人| 久99久视频精品免费| 日韩av在线大香蕉| 一级作爱视频免费观看| 亚洲国产精品久久男人天堂| 亚洲人成网站高清观看| 91老司机精品| 久久伊人香网站| 看片在线看免费视频| 国语自产精品视频在线第100页| 国产一区在线观看成人免费| 日本一二三区视频观看| 亚洲国产精品成人综合色| 欧美人与性动交α欧美精品济南到| 国产69精品久久久久777片 | 老熟妇仑乱视频hdxx| 久久精品国产亚洲av香蕉五月| 久久久久久大精品| 国产精品 国内视频| 国产真人三级小视频在线观看| www.999成人在线观看| 一边摸一边抽搐一进一小说| 很黄的视频免费| 欧美性猛交黑人性爽| 日韩高清综合在线| 成人三级做爰电影| 欧美乱色亚洲激情| 亚洲电影在线观看av| 在线观看舔阴道视频| 亚洲成人免费电影在线观看| 国产精品综合久久久久久久免费| 老熟妇仑乱视频hdxx| 久久久水蜜桃国产精品网| 制服诱惑二区| 高潮久久久久久久久久久不卡| 国产精品久久久久久久电影 | 亚洲av成人精品一区久久| av免费在线观看网站| 人人妻,人人澡人人爽秒播| 国产av麻豆久久久久久久| 在线十欧美十亚洲十日本专区| 99久久久亚洲精品蜜臀av| 久久精品夜夜夜夜夜久久蜜豆 | 一级片免费观看大全| 成人三级做爰电影| 国产精品乱码一区二三区的特点| 久久人人精品亚洲av| 99精品欧美一区二区三区四区| 久久人妻福利社区极品人妻图片| 97超级碰碰碰精品色视频在线观看| 亚洲av中文字字幕乱码综合| 精品电影一区二区在线| 久久久精品欧美日韩精品| 99热6这里只有精品| 精品少妇一区二区三区视频日本电影| 18禁美女被吸乳视频| 亚洲 欧美一区二区三区| 亚洲片人在线观看| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 成年版毛片免费区| 亚洲自拍偷在线| 国产精品久久久久久人妻精品电影| 精品日产1卡2卡| 亚洲欧美日韩东京热| 国产成人系列免费观看| 久久久久国内视频| 午夜精品一区二区三区免费看| 午夜福利高清视频| 黄色视频,在线免费观看| 国产亚洲欧美98| 啦啦啦韩国在线观看视频| 精品一区二区三区视频在线观看免费| 五月玫瑰六月丁香| 1024视频免费在线观看| 亚洲av电影不卡..在线观看| 欧美日本视频| 精品电影一区二区在线| 90打野战视频偷拍视频| 成人亚洲精品av一区二区| 超碰成人久久| 两个人免费观看高清视频| 欧美乱妇无乱码| 怎么达到女性高潮| 亚洲专区中文字幕在线| 亚洲人成77777在线视频| 美女 人体艺术 gogo| 超碰成人久久| 他把我摸到了高潮在线观看| 日日干狠狠操夜夜爽| 一级作爱视频免费观看| 在线播放国产精品三级| 777久久人妻少妇嫩草av网站| 无遮挡黄片免费观看| 在线观看美女被高潮喷水网站 | 给我免费播放毛片高清在线观看| 精品国内亚洲2022精品成人| 久久精品人妻少妇| 香蕉丝袜av| 看片在线看免费视频| 88av欧美| 91老司机精品| 亚洲熟妇中文字幕五十中出| 国产精品久久电影中文字幕| 麻豆久久精品国产亚洲av| 五月玫瑰六月丁香| 成人手机av| 中文在线观看免费www的网站 | 亚洲av电影在线进入| 老司机午夜十八禁免费视频| 亚洲国产日韩欧美精品在线观看 | 精品高清国产在线一区| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 一区二区三区国产精品乱码| 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 精品免费久久久久久久清纯| 午夜福利视频1000在线观看| 啦啦啦韩国在线观看视频| 亚洲国产精品sss在线观看| 十八禁人妻一区二区| 看黄色毛片网站| 久久精品夜夜夜夜夜久久蜜豆 | 人人妻人人看人人澡| 久久久久久人人人人人| 国产91精品成人一区二区三区| 大型av网站在线播放| 亚洲欧美精品综合久久99| 亚洲午夜精品一区,二区,三区| 亚洲成人中文字幕在线播放| 好男人在线观看高清免费视频| 757午夜福利合集在线观看| 亚洲熟妇中文字幕五十中出| 成人三级做爰电影| 久久久久久九九精品二区国产 | 欧美日韩精品网址| 国产精品免费一区二区三区在线| 久久性视频一级片| 久久精品亚洲精品国产色婷小说| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 三级国产精品欧美在线观看 | 桃红色精品国产亚洲av| 伦理电影免费视频| 神马国产精品三级电影在线观看 | 欧美最黄视频在线播放免费| 此物有八面人人有两片| 久久中文字幕人妻熟女| 午夜福利成人在线免费观看| 两个人的视频大全免费| 一个人免费在线观看电影 | 日韩欧美在线乱码| 高清在线国产一区| 亚洲av五月六月丁香网| 国产成年人精品一区二区| 久久午夜亚洲精品久久| 欧美大码av| 又黄又爽又免费观看的视频| 国产精品一及| 国产单亲对白刺激| 久久精品91蜜桃| 亚洲无线在线观看| 成人18禁在线播放| 亚洲 国产 在线| 日韩高清综合在线| 欧美日韩国产亚洲二区| 男女下面进入的视频免费午夜| 中亚洲国语对白在线视频| 天堂影院成人在线观看| 精品电影一区二区在线| 男女之事视频高清在线观看| 精品国产乱子伦一区二区三区| 美女午夜性视频免费| 在线观看免费日韩欧美大片| 欧美性长视频在线观看| 99精品欧美一区二区三区四区| 午夜久久久久精精品| 亚洲激情在线av| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| 国内久久婷婷六月综合欲色啪| 精品国产美女av久久久久小说| 日韩欧美国产在线观看| 国产欧美日韩精品亚洲av| 国产激情欧美一区二区| 成在线人永久免费视频| 国产熟女xx| 亚洲精品粉嫩美女一区| 看黄色毛片网站| 99国产综合亚洲精品| 亚洲国产日韩欧美精品在线观看 | 亚洲在线自拍视频| 久久久久亚洲av毛片大全| 麻豆av在线久日| 黄色毛片三级朝国网站| 一区福利在线观看| 在线看三级毛片| 小说图片视频综合网站| 亚洲精品国产精品久久久不卡| 成熟少妇高潮喷水视频| 成人国产一区最新在线观看| 久久久久国内视频| 国产免费男女视频| 18禁国产床啪视频网站| 精品欧美国产一区二区三| 少妇熟女aⅴ在线视频| 男女床上黄色一级片免费看| 国产视频内射| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看| 国产精品一区二区三区四区免费观看 | 男女视频在线观看网站免费 | 久久精品国产清高在天天线| av片东京热男人的天堂| 精品国产亚洲在线| 国产成人精品久久二区二区91| 一本精品99久久精品77| 日韩欧美精品v在线| 99在线视频只有这里精品首页| 国产精品一区二区三区四区免费观看 | 国产精品自产拍在线观看55亚洲| 亚洲在线自拍视频| 国产探花在线观看一区二区| 精品人妻1区二区| 十八禁人妻一区二区| 老司机深夜福利视频在线观看| 国产午夜精品久久久久久| 日本三级黄在线观看| 日韩欧美精品v在线| 亚洲人与动物交配视频| 妹子高潮喷水视频| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 国产精品爽爽va在线观看网站| 午夜精品一区二区三区免费看| 少妇被粗大的猛进出69影院| 久久久国产成人精品二区| 久久久久久人人人人人| 禁无遮挡网站| 精品少妇一区二区三区视频日本电影| 国产成人精品久久二区二区91| 国产亚洲精品一区二区www| 久久久久久久午夜电影| 国产在线精品亚洲第一网站| 亚洲美女视频黄频| 久久久久久久精品吃奶| 国内精品久久久久久久电影| 亚洲中文字幕一区二区三区有码在线看 | av免费在线观看网站| 久久久久久大精品| xxx96com| 婷婷六月久久综合丁香| 黄色片一级片一级黄色片| av福利片在线观看| 久久久久久国产a免费观看| 亚洲av熟女| 国产午夜精品论理片| 黄色视频不卡| 欧美大码av| 一进一出好大好爽视频| 欧美极品一区二区三区四区| 日韩欧美精品v在线| e午夜精品久久久久久久| 国产亚洲av嫩草精品影院| 黄色毛片三级朝国网站| 少妇人妻一区二区三区视频| 在线看三级毛片| 悠悠久久av| 成人av在线播放网站| 成人特级黄色片久久久久久久| 精品国产乱码久久久久久男人| 午夜精品一区二区三区免费看| 久久精品国产清高在天天线| tocl精华| 一本大道久久a久久精品| 久久这里只有精品中国| 黄色女人牲交| 亚洲一区高清亚洲精品| 欧美三级亚洲精品| 久久精品91无色码中文字幕| 日韩国内少妇激情av| 日本免费a在线| 久久久久久久久中文| 免费观看人在逋| 国产亚洲av高清不卡| 成人av在线播放网站| 搡老岳熟女国产| 亚洲人成网站高清观看| 香蕉丝袜av| 一本久久中文字幕| 午夜免费成人在线视频| 精品久久久久久久久久久久久| 国产精品日韩av在线免费观看| 国内精品一区二区在线观看| 99在线视频只有这里精品首页| 在线视频色国产色| 亚洲人成网站高清观看| 免费观看人在逋| 一区二区三区激情视频| 欧美丝袜亚洲另类 | tocl精华| 50天的宝宝边吃奶边哭怎么回事| 老熟妇仑乱视频hdxx| 麻豆av在线久日| 一级毛片精品| 午夜免费观看网址| 中文字幕最新亚洲高清| 丰满的人妻完整版| 欧美久久黑人一区二区| 国产精品av久久久久免费| 亚洲av中文字字幕乱码综合| 亚洲成人国产一区在线观看| av超薄肉色丝袜交足视频| 国产亚洲精品综合一区在线观看 | 久久久久性生活片| 亚洲一区二区三区色噜噜| 毛片女人毛片| 日本成人三级电影网站| 日本a在线网址| 亚洲一区二区三区不卡视频| 亚洲精品美女久久av网站| 可以在线观看毛片的网站| 国产av不卡久久| 天堂影院成人在线观看| 一进一出好大好爽视频| 制服人妻中文乱码| 国产精品久久久久久精品电影| 欧美一级a爱片免费观看看 | 三级男女做爰猛烈吃奶摸视频| 国产成人av教育| 日本 av在线| 男男h啪啪无遮挡| 麻豆av在线久日| 婷婷精品国产亚洲av在线| 级片在线观看| 日韩欧美免费精品| 日韩中文字幕欧美一区二区| 两性夫妻黄色片| 99精品在免费线老司机午夜| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| 婷婷六月久久综合丁香| 19禁男女啪啪无遮挡网站| 怎么达到女性高潮| 久久久国产精品麻豆| 午夜精品一区二区三区免费看| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 亚洲 欧美 日韩 在线 免费| 99国产精品一区二区蜜桃av| 欧美日韩瑟瑟在线播放| 久久久水蜜桃国产精品网| 亚洲欧美日韩高清在线视频| 天堂动漫精品| 精品国产乱子伦一区二区三区| 亚洲av美国av| 国产精品 欧美亚洲| 两个人免费观看高清视频| 精品日产1卡2卡| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 男女午夜视频在线观看| 亚洲人成伊人成综合网2020| ponron亚洲| 久久九九热精品免费| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放 | 91麻豆精品激情在线观看国产| 国产熟女xx| 亚洲成av人片免费观看| 男人舔女人下体高潮全视频| 在线国产一区二区在线| 国产亚洲av嫩草精品影院| 中文字幕人成人乱码亚洲影| 国产在线精品亚洲第一网站| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 亚洲男人天堂网一区| 亚洲人成网站高清观看| 久久久久久国产a免费观看| 老司机深夜福利视频在线观看| 又大又爽又粗| 一区二区三区激情视频| 久久久久国内视频| 成人特级黄色片久久久久久久| 亚洲成av人片免费观看| 此物有八面人人有两片| 男女那种视频在线观看| 亚洲国产精品成人综合色| 日本五十路高清| 啦啦啦观看免费观看视频高清| 日本 欧美在线| av中文乱码字幕在线| 一级毛片女人18水好多| 麻豆av在线久日| 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 国产精品精品国产色婷婷| 久久久久国产精品人妻aⅴ院| 黄片大片在线免费观看| 操出白浆在线播放| 午夜福利高清视频| 中出人妻视频一区二区| 99热只有精品国产| 亚洲中文日韩欧美视频| 免费搜索国产男女视频| 午夜视频精品福利|