• <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型中子星與黑洞表面溫度分布的研究
    亚洲精品一卡2卡三卡4卡5卡| 精品一区二区三区视频在线观看免费 | 久久国产亚洲av麻豆专区| 国产乱人伦免费视频| 免费看a级黄色片| 国产精品二区激情视频| a级毛片在线看网站| 久久久国产成人精品二区 | 久久久久久久精品吃奶| aaaaa片日本免费| 精品国产乱子伦一区二区三区| 好看av亚洲va欧美ⅴa在| 嫩草影视91久久| 深夜精品福利| 久久精品成人免费网站| www国产在线视频色| 欧美大码av| 97人妻天天添夜夜摸| 欧美日韩亚洲高清精品| 老熟妇仑乱视频hdxx| 激情视频va一区二区三区| 女性生殖器流出的白浆| 纯流量卡能插随身wifi吗| 搡老乐熟女国产| 亚洲男人的天堂狠狠| 欧美激情高清一区二区三区| 99久久99久久久精品蜜桃| 国产精品一区二区免费欧美| 在线观看午夜福利视频| 69精品国产乱码久久久| 欧美色视频一区免费| 99在线视频只有这里精品首页| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国语在线视频| 黄色怎么调成土黄色| 国产精品亚洲一级av第二区| av欧美777| 国产成人精品久久二区二区91| 久久久国产成人精品二区 | 女生性感内裤真人,穿戴方法视频| 在线观看免费日韩欧美大片| 国产精品一区二区在线不卡| 亚洲欧美精品综合一区二区三区| 男女下面进入的视频免费午夜 | 国产一区在线观看成人免费| 丝袜美腿诱惑在线| 国产麻豆69| 看片在线看免费视频| 99久久久亚洲精品蜜臀av| 亚洲精品美女久久久久99蜜臀| 在线观看免费日韩欧美大片| 亚洲情色 制服丝袜| 久久香蕉精品热| 亚洲成人国产一区在线观看| 麻豆av在线久日| 天堂影院成人在线观看| 国产亚洲av高清不卡| 国产欧美日韩一区二区三| 电影成人av| 夫妻午夜视频| 水蜜桃什么品种好| 欧美黑人精品巨大| 一级片免费观看大全| 欧美中文综合在线视频| 97人妻天天添夜夜摸| 久久精品国产亚洲av香蕉五月| 人成视频在线观看免费观看| 成人三级做爰电影| 色婷婷久久久亚洲欧美| 国产又爽黄色视频| 亚洲aⅴ乱码一区二区在线播放 | 久久久水蜜桃国产精品网| 真人一进一出gif抽搐免费| 黄片播放在线免费| 热re99久久精品国产66热6| 丰满的人妻完整版| 视频在线观看一区二区三区| 色尼玛亚洲综合影院| 亚洲专区国产一区二区| 精品国产美女av久久久久小说| 高清毛片免费观看视频网站 | 精品一品国产午夜福利视频| 夜夜爽天天搞| 成人18禁在线播放| 麻豆成人av在线观看| 精品国产美女av久久久久小说| 十八禁人妻一区二区| 久久人妻av系列| 天堂俺去俺来也www色官网| 日韩欧美三级三区| 男人操女人黄网站| 精品国内亚洲2022精品成人| 国产精品综合久久久久久久免费 | 欧美日韩乱码在线| 久久久久国产一级毛片高清牌| 色婷婷久久久亚洲欧美| 丁香欧美五月| 91av网站免费观看| 国产97色在线日韩免费| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人精品中文字幕电影 | 久久精品国产综合久久久| 亚洲欧美日韩无卡精品| 国产伦人伦偷精品视频| 99久久综合精品五月天人人| 老司机在亚洲福利影院| 91麻豆精品激情在线观看国产 | 亚洲自拍偷在线| 久久99一区二区三区| 免费av中文字幕在线| 九色亚洲精品在线播放| 91麻豆精品激情在线观看国产 | 亚洲国产精品999在线| 少妇裸体淫交视频免费看高清 | 国产亚洲精品第一综合不卡| 亚洲欧美一区二区三区黑人| 日本三级黄在线观看| 又紧又爽又黄一区二区| 91精品三级在线观看| 国产国语露脸激情在线看| av电影中文网址| 国产真人三级小视频在线观看| 午夜福利在线观看吧| 久久热在线av| ponron亚洲| 俄罗斯特黄特色一大片| 中出人妻视频一区二区| 黄色成人免费大全| 亚洲av美国av| 午夜精品久久久久久毛片777| 99精品在免费线老司机午夜| 新久久久久国产一级毛片| 在线看a的网站| 国产精品 欧美亚洲| 一级毛片高清免费大全| 欧美日韩瑟瑟在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久99一区二区三区| tocl精华| 首页视频小说图片口味搜索| 亚洲欧美日韩高清在线视频| 亚洲国产精品合色在线| 老司机靠b影院| 一区二区日韩欧美中文字幕| av在线播放免费不卡| 女性生殖器流出的白浆| 久久性视频一级片| 精品一区二区三区四区五区乱码| svipshipincom国产片| 韩国av一区二区三区四区| 日韩欧美国产一区二区入口| 啦啦啦 在线观看视频| 亚洲av熟女| 高清在线国产一区| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 人妻久久中文字幕网| 丰满饥渴人妻一区二区三| 亚洲av成人av| 免费在线观看亚洲国产| a级毛片黄视频| 男女床上黄色一级片免费看| 免费在线观看亚洲国产| 国产精品99久久99久久久不卡| 88av欧美| 日本vs欧美在线观看视频| 日韩精品中文字幕看吧| 日韩欧美三级三区| 亚洲欧美激情综合另类| av免费在线观看网站| 一二三四社区在线视频社区8| 欧美日韩一级在线毛片| 91麻豆av在线| 老汉色av国产亚洲站长工具| 亚洲av美国av| 欧美成人性av电影在线观看| 一级片免费观看大全| 久久久久久亚洲精品国产蜜桃av| svipshipincom国产片| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| www.自偷自拍.com| 亚洲精品国产精品久久久不卡| 国产精品免费一区二区三区在线| 国产一区二区三区在线臀色熟女 | 亚洲色图 男人天堂 中文字幕| 法律面前人人平等表现在哪些方面| 男男h啪啪无遮挡| 精品一区二区三区视频在线观看免费 | 99国产综合亚洲精品| 多毛熟女@视频| 999精品在线视频| 国产精品av久久久久免费| 大码成人一级视频| 久久久久久免费高清国产稀缺| 琪琪午夜伦伦电影理论片6080| 久久久久国产一级毛片高清牌| 久久久久久免费高清国产稀缺| 欧美日韩亚洲综合一区二区三区_| 级片在线观看| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看 | 国产成人一区二区三区免费视频网站| 99国产精品一区二区三区| 深夜精品福利| 成年女人毛片免费观看观看9| 妹子高潮喷水视频| 中文字幕另类日韩欧美亚洲嫩草| svipshipincom国产片| 自线自在国产av| 免费人成视频x8x8入口观看| 黄色怎么调成土黄色| 一级a爱片免费观看的视频| tocl精华| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 精品久久蜜臀av无| 久久婷婷成人综合色麻豆| 在线看a的网站| 窝窝影院91人妻| 国产不卡一卡二| 别揉我奶头~嗯~啊~动态视频| 免费在线观看视频国产中文字幕亚洲| 亚洲久久久国产精品| 久久久久国产一级毛片高清牌| 久久婷婷成人综合色麻豆| 最近最新中文字幕大全免费视频| 免费观看精品视频网站| 一a级毛片在线观看| 一进一出抽搐动态| 亚洲专区国产一区二区| 极品教师在线免费播放| 最新美女视频免费是黄的| 精品人妻1区二区| 高清毛片免费观看视频网站 | 日本黄色日本黄色录像| 亚洲欧美精品综合久久99| 老熟妇仑乱视频hdxx| 亚洲黑人精品在线| 91精品国产国语对白视频| 在线观看免费视频网站a站| 一进一出抽搐动态| 女同久久另类99精品国产91| 黄频高清免费视频| 黄色怎么调成土黄色| 99re在线观看精品视频| 国产成人一区二区三区免费视频网站| 亚洲专区国产一区二区| 99精国产麻豆久久婷婷| 黄片播放在线免费| 男人舔女人下体高潮全视频| 嫁个100分男人电影在线观看| xxxhd国产人妻xxx| 午夜两性在线视频| 国产91精品成人一区二区三区| 免费观看人在逋| 亚洲欧美精品综合久久99| 黄网站色视频无遮挡免费观看| 成年人免费黄色播放视频| av超薄肉色丝袜交足视频| 中出人妻视频一区二区| 级片在线观看| 男男h啪啪无遮挡| 久久婷婷成人综合色麻豆| 搡老岳熟女国产| 国产av一区二区精品久久| 精品久久久久久久毛片微露脸| 日本a在线网址| 午夜精品久久久久久毛片777| 在线播放国产精品三级| 怎么达到女性高潮| 亚洲精品成人av观看孕妇| 亚洲国产精品sss在线观看 | 成年人黄色毛片网站| 欧美久久黑人一区二区| 精品国产亚洲在线| av视频免费观看在线观看| 男女之事视频高清在线观看| 亚洲欧美精品综合久久99| 757午夜福利合集在线观看| 日本欧美视频一区| cao死你这个sao货| 亚洲九九香蕉| 中出人妻视频一区二区| 成人国产一区最新在线观看| 超色免费av| 久久精品亚洲熟妇少妇任你| 国产精品日韩av在线免费观看 | 999久久久精品免费观看国产| 国产成人精品无人区| 亚洲午夜理论影院| 99国产精品免费福利视频| 美女扒开内裤让男人捅视频| 韩国精品一区二区三区| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| av有码第一页| 丰满迷人的少妇在线观看| 国产又色又爽无遮挡免费看| 精品久久久久久电影网| 亚洲精品中文字幕在线视频| 精品福利永久在线观看| 香蕉久久夜色| 亚洲中文av在线| 最近最新中文字幕大全电影3 | 精品国产国语对白av| 日韩欧美三级三区| 欧美中文综合在线视频| 天天影视国产精品| 国产麻豆69| 欧美日韩黄片免| 亚洲欧美一区二区三区黑人| av国产精品久久久久影院| 午夜a级毛片| 成人国语在线视频| 伊人久久大香线蕉亚洲五| 中亚洲国语对白在线视频| 日本三级黄在线观看| 老汉色∧v一级毛片| 男女之事视频高清在线观看| 丰满人妻熟妇乱又伦精品不卡| 少妇粗大呻吟视频| 1024视频免费在线观看| 国产乱人伦免费视频| 亚洲精品中文字幕一二三四区| 国产成人精品久久二区二区免费| 天堂动漫精品| 黑人欧美特级aaaaaa片| 国产91精品成人一区二区三区| 国产亚洲欧美精品永久| 亚洲成人精品中文字幕电影 | 欧美亚洲日本最大视频资源| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| av在线天堂中文字幕 | 精品乱码久久久久久99久播| 亚洲男人的天堂狠狠| 国产深夜福利视频在线观看| 欧美色视频一区免费| 妹子高潮喷水视频| 中文欧美无线码| 在线观看66精品国产| 精品福利永久在线观看| 色老头精品视频在线观看| 51午夜福利影视在线观看| 啦啦啦在线免费观看视频4| 亚洲国产看品久久| 巨乳人妻的诱惑在线观看| av在线播放免费不卡| 视频在线观看一区二区三区| 亚洲中文日韩欧美视频| av视频免费观看在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 美国免费a级毛片| 久久精品影院6| √禁漫天堂资源中文www| 超碰97精品在线观看| 久久精品国产清高在天天线| 超色免费av| 精品卡一卡二卡四卡免费| av视频免费观看在线观看| 久久精品成人免费网站| 欧美 亚洲 国产 日韩一| 国产又色又爽无遮挡免费看| 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看 | 十分钟在线观看高清视频www| 韩国av一区二区三区四区| 99精品欧美一区二区三区四区| 91麻豆av在线| 亚洲专区中文字幕在线| 1024视频免费在线观看| 热re99久久国产66热| 精品人妻在线不人妻| 女人被狂操c到高潮| 激情视频va一区二区三区| 性色av乱码一区二区三区2| 法律面前人人平等表现在哪些方面| 国产av在哪里看| 日韩av在线大香蕉| 人妻丰满熟妇av一区二区三区| 在线观看www视频免费| 中亚洲国语对白在线视频| 久久国产亚洲av麻豆专区| 婷婷六月久久综合丁香| 中文字幕人妻熟女乱码| 国产精品99久久99久久久不卡| 欧美日韩乱码在线| 久久精品国产综合久久久| 天堂√8在线中文| 精品一区二区三区av网在线观看| 欧美在线黄色| a级毛片在线看网站| 超碰成人久久| 国产精品久久久人人做人人爽| 久9热在线精品视频| 热99国产精品久久久久久7| 婷婷六月久久综合丁香| 99久久国产精品久久久| 亚洲精品一卡2卡三卡4卡5卡| av有码第一页| 男女高潮啪啪啪动态图| 精品国产一区二区三区四区第35| 午夜福利免费观看在线| 精品日产1卡2卡| 精品人妻1区二区| 久久久久久久久久久久大奶| 国产精品一区二区在线不卡| 欧美+亚洲+日韩+国产| 免费人成视频x8x8入口观看| 99精品在免费线老司机午夜| 天天影视国产精品| 成人手机av| 热99国产精品久久久久久7| 亚洲久久久国产精品| 老司机福利观看| 亚洲精品在线美女| 国产伦人伦偷精品视频| 国产成人av激情在线播放| 亚洲精品一二三| 亚洲男人的天堂狠狠| 久久狼人影院| 欧美丝袜亚洲另类 | 级片在线观看| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av| 免费久久久久久久精品成人欧美视频| 久久久久精品国产欧美久久久| 亚洲国产毛片av蜜桃av| 真人一进一出gif抽搐免费| 欧美激情久久久久久爽电影 | 国产精品久久久人人做人人爽| 如日韩欧美国产精品一区二区三区| 国产精品99久久99久久久不卡| 丝袜美腿诱惑在线| 啦啦啦在线免费观看视频4| 18禁国产床啪视频网站| 欧洲精品卡2卡3卡4卡5卡区| 在线观看一区二区三区激情| 99香蕉大伊视频| 妹子高潮喷水视频| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 国产伦人伦偷精品视频| 91在线观看av| 亚洲男人天堂网一区| 淫妇啪啪啪对白视频| 成人免费观看视频高清| 老司机亚洲免费影院| 这个男人来自地球电影免费观看| 久久婷婷成人综合色麻豆| 久久这里只有精品19| 国产野战对白在线观看| 免费av中文字幕在线| 欧美黑人欧美精品刺激| 看黄色毛片网站| 在线看a的网站| 国产精品亚洲一级av第二区| 超色免费av| 少妇 在线观看| 亚洲精华国产精华精| 久久香蕉精品热| www.www免费av| 水蜜桃什么品种好| 这个男人来自地球电影免费观看| 999久久久国产精品视频| 成人国语在线视频| 国产欧美日韩一区二区三| 亚洲精品av麻豆狂野| 香蕉丝袜av| 一二三四在线观看免费中文在| 真人一进一出gif抽搐免费| 午夜免费成人在线视频| 国产激情久久老熟女| 99热只有精品国产| 精品卡一卡二卡四卡免费| 亚洲一码二码三码区别大吗| 欧美日韩黄片免| 色婷婷av一区二区三区视频| 久热这里只有精品99| 淫妇啪啪啪对白视频| 亚洲狠狠婷婷综合久久图片| 欧美日韩精品网址| 国产熟女午夜一区二区三区| 91成人精品电影| 多毛熟女@视频| 欧美日本中文国产一区发布| 亚洲一区二区三区不卡视频| 国产黄色免费在线视频| 久久人人97超碰香蕉20202| 久久精品国产99精品国产亚洲性色 | 亚洲欧美激情在线| 精品一区二区三区四区五区乱码| 午夜福利欧美成人| 日韩大尺度精品在线看网址 | 中文字幕色久视频| 91麻豆精品激情在线观看国产 | 男女下面插进去视频免费观看| 一边摸一边做爽爽视频免费| 日本wwww免费看| av在线播放免费不卡| 两个人看的免费小视频| av在线天堂中文字幕 | 亚洲一码二码三码区别大吗| 激情视频va一区二区三区| 精品一区二区三区四区五区乱码| 精品一品国产午夜福利视频| 高潮久久久久久久久久久不卡| 久久中文字幕一级| 变态另类成人亚洲欧美熟女 | 亚洲中文日韩欧美视频| 大型av网站在线播放| 亚洲精品国产区一区二| 两人在一起打扑克的视频| 亚洲片人在线观看| 日韩大尺度精品在线看网址 | 这个男人来自地球电影免费观看| 亚洲av五月六月丁香网| 免费在线观看黄色视频的| 亚洲国产毛片av蜜桃av| 久久人妻av系列| 久久中文字幕一级| 久久精品91蜜桃| 18禁黄网站禁片午夜丰满| 男女下面进入的视频免费午夜 | 欧美日韩福利视频一区二区| 国产精品影院久久| 一级片'在线观看视频| 亚洲精品美女久久久久99蜜臀| 黄色怎么调成土黄色| 一夜夜www| 久久香蕉精品热| 精品无人区乱码1区二区| 男女做爰动态图高潮gif福利片 | 欧美性长视频在线观看| 欧美色视频一区免费| 免费少妇av软件| av国产精品久久久久影院| 高清黄色对白视频在线免费看| 久久婷婷成人综合色麻豆| 色综合欧美亚洲国产小说| 国产精品二区激情视频| 免费日韩欧美在线观看| 在线观看免费日韩欧美大片| 涩涩av久久男人的天堂| 国产在线精品亚洲第一网站| 一个人免费在线观看的高清视频| 国产精品98久久久久久宅男小说| 精品无人区乱码1区二区| 男女做爰动态图高潮gif福利片 | 久久久久国产精品人妻aⅴ院| 欧美成狂野欧美在线观看| 亚洲国产毛片av蜜桃av| 欧美最黄视频在线播放免费 | 日日爽夜夜爽网站| 香蕉丝袜av| 日本黄色视频三级网站网址| 久久精品国产99精品国产亚洲性色 | 久久婷婷成人综合色麻豆| 欧美一级毛片孕妇| 黄色视频不卡| 不卡av一区二区三区| 免费人成视频x8x8入口观看| 成年版毛片免费区| 村上凉子中文字幕在线| 亚洲av电影在线进入| 亚洲专区中文字幕在线| 久久久精品国产亚洲av高清涩受| 99国产极品粉嫩在线观看| 在线观看免费午夜福利视频| 91字幕亚洲| 免费在线观看黄色视频的| ponron亚洲| 免费少妇av软件| 免费在线观看黄色视频的| 老熟妇仑乱视频hdxx| 成人黄色视频免费在线看| 老汉色∧v一级毛片| 免费在线观看黄色视频的| 亚洲精品在线观看二区| 亚洲成a人片在线一区二区| 99国产极品粉嫩在线观看| 十分钟在线观看高清视频www| 久久国产精品男人的天堂亚洲| 婷婷丁香在线五月| 国产精品免费一区二区三区在线| 黄色怎么调成土黄色| 亚洲精品国产精品久久久不卡| 在线国产一区二区在线| 久久久久久久午夜电影 | 亚洲国产精品999在线| 精品日产1卡2卡| 亚洲成人国产一区在线观看| 黄色毛片三级朝国网站| 国产精品av久久久久免费| 日韩国内少妇激情av| 一级片免费观看大全| 天堂俺去俺来也www色官网| 精品熟女少妇八av免费久了| 村上凉子中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美三级三区| 久久国产亚洲av麻豆专区| 欧美激情高清一区二区三区| 国产精品一区二区精品视频观看| 国产精华一区二区三区| 国产欧美日韩一区二区三区在线| 国产97色在线日韩免费| av欧美777| 久久国产精品人妻蜜桃| 淫秽高清视频在线观看|