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

    中街山列島海域曼氏無針烏賊棲息地偏向性研究

    2022-12-12 11:41:40周永東劉姝含胡揚(yáng)杰徐開達(dá)
    關(guān)鍵詞:環(huán)境模型

    吳 天,梁 君,周永東,劉姝含,胡揚(yáng)杰,趙 瑞,徐開達(dá)

    (浙江海洋大學(xué)海洋與漁業(yè)研究所,浙江省海洋水產(chǎn)研究所,浙江省海洋漁業(yè)資源可持續(xù)利用技術(shù)研究重點(diǎn)實(shí)驗(yàn)室,浙江舟山 316021)

    曼氏無針烏賊Sepiella maindroni俗稱墨魚,是分布廣泛的淺海種類,隸屬軟體動(dòng)物門Mollusca、頭足綱Cephalopoda、烏賊目Sepiida、烏賊科Sepiidae、無針烏賊屬Sepiella,主要分布在日本海、馬來群島海域和印度洋東海岸,在我國(guó)東海和黃海均有分布,其中心分布區(qū)在浙江近海和閩東海域。歷史上東海區(qū)曼氏無針烏賊最高年產(chǎn)超過7×104t,曾與大黃魚Larimichthys crocea、小黃魚Larimichthys polyactis、帶魚Trichiurus lepturus并稱為我國(guó)四大海產(chǎn)[1-3]。但是20 世紀(jì)80 年代以來,由于過度捕撈和對(duì)資源的不合理利用,加之全球變暖造成棲息地退化,使得曼氏無針烏賊資源出現(xiàn)嚴(yán)重衰退現(xiàn)象,資源變動(dòng)較大[4]。中街山列島海域位于舟山漁場(chǎng)的中心區(qū)域,受沿岸水與外洋水的影響,水體交換快、餌料充足,是我國(guó)主要經(jīng)濟(jì)魚類生長(zhǎng)、索餌、繁殖的重要場(chǎng)所,歷史上素有“魚倉(cāng)”之美譽(yù),曾是曼氏無針烏賊的重要產(chǎn)卵場(chǎng)之一[5]。

    棲息地是指生物個(gè)體和種群居住的場(chǎng)所,是由生物出現(xiàn)在環(huán)境中的空間范圍和環(huán)境條件綜合構(gòu)成,棲息地的選擇是指生物為了達(dá)到生存目的而尋找相對(duì)適宜棲息地的過程[6]。棲息地的研究一直以來都是動(dòng)物生態(tài)學(xué)的研究重點(diǎn),當(dāng)棲息地質(zhì)量下降時(shí),會(huì)影響物種的數(shù)量,甚至導(dǎo)致物種滅絕[7]。在漁業(yè)資源研究中,棲息地的研究?jī)?nèi)容主要是生物棲息環(huán)境的變化對(duì)生物活動(dòng)的影響,生物對(duì)棲息地的選擇是受多種環(huán)境因子共同作用的結(jié)果,單因子分析并不能很好地解釋生物對(duì)棲息地的選擇機(jī)制,而廣義可加模型(GAM)能夠直接處理響應(yīng)變量與解釋變量之間的非線性關(guān)系,通過非參數(shù)的方法檢驗(yàn)數(shù)據(jù)結(jié)構(gòu)[8-9]。本研究基于2018-2021 年中街山列島海域漁業(yè)調(diào)查數(shù)據(jù)和環(huán)境數(shù)據(jù),通過GAM 構(gòu)建非參數(shù)的擬合模型,探究中街山列島海域曼氏無針烏賊資源分布與環(huán)境因子之間的關(guān)系,了解曼氏無針烏賊的棲息地選擇機(jī)制,旨在為曼氏無針烏賊的種群保護(hù)及產(chǎn)卵場(chǎng)生物多樣性保護(hù)提供理論基礎(chǔ),為浙江漁場(chǎng)修復(fù)振興計(jì)劃的實(shí)施提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域和站位布設(shè)

    于2018-2021 年春季(4 月)和秋季(10 月)在中街山列島海域進(jìn)行了8 個(gè)航次的底拖網(wǎng)漁業(yè)資源調(diào)查,調(diào)查范圍如圖1 所示。本研究中站位設(shè)置采取區(qū)域隨機(jī)分布的原則,沿潮流主軸(西北-東南)共設(shè)置了16 個(gè)站位。

    圖1 調(diào)查站位示意圖Fig.1 Sketch of survey station

    1.2 采樣方法和樣品處理

    調(diào)查船為58.8 kW 單拖漁船,拖速為2 kn,拖曳時(shí)間為30 min,調(diào)查網(wǎng)具網(wǎng)口周長(zhǎng)25 m,囊網(wǎng)網(wǎng)目尺寸3.0 cm。在每一次調(diào)查時(shí),使用溫鹽深測(cè)定儀(CTD)同步采集環(huán)境數(shù)據(jù),并在各站位采集水樣。共測(cè)定水深、表層海水水溫、表層海水鹽度、底層海水水溫、底層海水鹽度、葉綠素a、溶解氧、懸浮物、無機(jī)氮和活性磷酸鹽10 類環(huán)境數(shù)據(jù)。樣品采集、處理和分析等均按照《海洋漁業(yè)資源調(diào)查規(guī)范》[10]的規(guī)定執(zhí)行。

    1.3 數(shù)據(jù)處理和方法分析

    1.3.1 資源密度和生物量估算

    曼氏無針烏賊的資源密度估算和生物量估算按照以下公式[11]計(jì)算:

    式中:D 為現(xiàn)存資源密度,ind.·km-2;C 為每小時(shí)取樣面積內(nèi)的漁獲量,kg;q 為網(wǎng)具漁獲率,由于曼氏無針烏賊屬于沿岸暖水底棲種,q 取0.5;A 為網(wǎng)具每小時(shí)掃海面積,km2;B 為生物量,kg·km-2;為個(gè)體平均質(zhì)量,kg·ind-2。

    1.3.2 環(huán)境因子的相關(guān)性分析

    在建模中,需要對(duì)環(huán)境因子進(jìn)行相關(guān)性分析,克服共線性較高的解釋變量對(duì)參數(shù)預(yù)估準(zhǔn)確性的干擾,當(dāng)2 個(gè)環(huán)境因子間存在共線關(guān)系時(shí),保留1 個(gè)變量作為解釋變量[12]。本研究采用Pearson 相關(guān)性分析對(duì)于環(huán)境因子之間可能存在共線性關(guān)系的環(huán)境因子進(jìn)行相關(guān)性分析,并通過方差膨脹因子(variance inflation factors,VIF)對(duì)所有候選變量進(jìn)行檢驗(yàn)。

    1.3.3 廣義可加模型(GAM)分析

    本研究以曼氏無針烏賊資源分布為響應(yīng)變量,其中季節(jié)作為分類變量,各環(huán)境因子為解釋變量。由于曼氏無針烏賊的零值站位出現(xiàn)的次數(shù)較多,假設(shè)曼氏無針烏賊分布模式為二項(xiàng)分布,GAM 模型選擇logistic 為鏈接函數(shù)模擬中街山列島海域曼氏無針烏賊資源分布同環(huán)境因子的關(guān)系[12-13]。GAM 的表達(dá)式為:

    式中:P 為曼氏無針烏賊出現(xiàn)的概率;α 表示截距項(xiàng);Si 表示樣條平滑函數(shù);xi表示第i 個(gè)解釋變量;ε 表示誤差項(xiàng)。

    GAM 模型是在加性模型基礎(chǔ)上提出的非參數(shù)擬合模型,能夠直接處理響應(yīng)變量與多個(gè)解釋變量之間的非線性關(guān)系。在模型的擬合過程中,首先確定曼氏無針烏賊資源分布與單一環(huán)境因子的關(guān)系,通過逐步回歸法將篩選后環(huán)境因子加入模型中,然后利用赤池信息準(zhǔn)則(akaike information criterion,AIC)逐步檢驗(yàn)加入不同環(huán)境因子后的模型擬合度,通過F 檢驗(yàn)來評(píng)估各因子的顯著性,以檢驗(yàn)?zāi)P椭忻總€(gè)環(huán)境因子的重要性,最終遴選出最優(yōu)模型[13-14]。

    數(shù)據(jù)通過Excel 軟件計(jì)算和整理,曼氏無針烏賊資源分布圖采用ArcGIS 10.2 軟件繪制,Pearson 相關(guān)性分析和VIF 檢驗(yàn)通過SPSS 23.0 實(shí)現(xiàn),GAM 模型構(gòu)建與檢驗(yàn)及可視化分析通過R 4.2.0 實(shí)現(xiàn)。

    2 結(jié)果

    2.1 曼氏無針烏賊時(shí)空分布特征

    曼氏無針烏賊資源空間分布季節(jié)變化明顯。春季,曼氏無針烏賊主要分布于島礁外側(cè),東北部居多(圖2)。秋季,曼氏無針烏賊主要分布于島礁內(nèi)的東南和西南方向(圖3)。時(shí)間因子月份對(duì)曼氏無針烏賊出現(xiàn)概率有明顯影響,2018-2021 年曼氏無針烏賊的出現(xiàn)概率秋季明顯高于春季。但春季曼氏無針烏賊的體質(zhì)量高于秋季(圖4)。

    圖2 2018-2021 春季曼氏無針烏賊尾數(shù)密度和生物量空間分布圖Fig.2 2018-2021 spring distribution of S.maindroni tail density and biomass

    圖3 2018-2021 秋季曼氏無針烏賊尾數(shù)密度和生物量空間分布圖Fig.3 2018-2021 autumn distribution of S.maindroni tail density and biomass

    圖4 2018-2021 曼氏無針烏賊尾數(shù)密度和生物量對(duì)照?qǐng)D(a.尾數(shù)密度;b.生物量)Fig.4 2018-2021 S.maindroni comparison of tail density and biomass(a.tail density;b.biomass)

    2.2 關(guān)鍵環(huán)境因子識(shí)別

    對(duì)于海洋生物,考慮的主要環(huán)境因子有溫度、鹽度、海表高度和葉綠素a等[5]。依據(jù)Pearson 相關(guān)性分析以及可視化結(jié)果(圖5)顯示,春季表層水溫和底層水溫、表層水溫和底層鹽度、底層水溫和底層鹽度、底層鹽度和活性磷酸鹽、葉綠素a和活性磷酸鹽存在顯著相關(guān)(P<0.05),秋季表層水溫、表層鹽度、底層水溫、底層鹽度之間均存在顯著相關(guān),無機(jī)氮和活性磷酸鹽之間存在顯著相關(guān),同時(shí)對(duì)各解釋變量進(jìn)行多重線性檢驗(yàn),方差膨脹系數(shù)閾值設(shè)定為4,去掉相關(guān)性和方差膨脹系數(shù)均高于閾值的變量,因此本文選取水深、表層水溫、表層鹽度、懸浮物、葉綠素a、溶解氧和無機(jī)氮作為春季曼氏無針烏賊棲息地模型構(gòu)建的解釋變量;水深、懸浮物、底層水溫、葉綠素a、溶解氧和活性磷酸鹽作為秋季曼氏無針烏賊棲息地模型構(gòu)建的解釋變量。

    圖5 春季和秋季不同環(huán)境因子相關(guān)性分析熱圖(a 為春季;b 為秋季)Fig.5 Correlation analysis of different environmental factors in spring and autumn(a.spring;b.autumn)

    2.3 GAM 模型構(gòu)建

    在GAM 模型中加入篩選后的環(huán)境因子,依據(jù)AIC 最小原則構(gòu)建曼氏無針烏賊棲息地選擇的最優(yōu)模型。結(jié)果如下:

    春季曼氏無針烏賊棲息地選擇的最優(yōu)模型:

    秋季曼氏無針烏賊棲息地選擇的最優(yōu)模型:

    式中:P為曼氏無針烏賊出現(xiàn)的概率;Depth 為水深,m;Chl-a為葉綠素a,μg·L-1;DO 為溶解氧,mg·L-1;ε 為相應(yīng)的GAM 模型截距(表1)。

    春季(GAM1)模型最佳預(yù)測(cè)變量組合為Depth+Chl-a(表1),模型累計(jì)偏差解釋率為70.35%,水深和葉綠素a 對(duì)曼氏無針烏賊的出現(xiàn)概率具有顯著影響(P<0.05),其貢獻(xiàn)率分別為23.94%和46.41%(表2)。秋季(GAM2)模型最佳預(yù)測(cè)變量為DO(表1),模型累計(jì)偏差解釋率為24.93%(表2)。

    表1 GAM 模型變量篩選及影響因子的參數(shù)分析Tab.1 Variable screening process for GAM and parameters analysis

    表2 模型擬合結(jié)果及各因子參數(shù)Tab.2 Result from models fitting and parameters of each factor

    模型結(jié)果顯示春季(GAM1)水深、葉綠素a 是影響曼氏無針烏賊棲息地選擇的主要因素。春季曼氏無針烏賊主要出現(xiàn)的水深范圍在25~35 m 之間,呈先上升后下降的趨勢(shì),30 m 時(shí)曼氏無針烏賊出現(xiàn)的概率最高(圖6a)。葉綠素a 在0.5~1 μg·L-1時(shí)曼氏無針烏賊出現(xiàn)概率為上升趨勢(shì),當(dāng)葉綠素a 超過1 μg·L-1時(shí)曼氏無針烏賊出現(xiàn)的概率開始下降(圖6b)。秋季,影響曼氏無針烏賊出現(xiàn)概率的單一環(huán)境因子為溶解氧,范圍在5.6~6.2 mg·L-1之間,呈先上升后下降的趨勢(shì)(圖6c)。

    圖6 GAM1 和GAM2 環(huán)境因子對(duì)曼氏無針烏賊分布的影響Fig.6 Effects of GAM1 and GAM2 of environmental factors on the S.maindroni

    3 討論

    3.1 曼氏無針烏賊的時(shí)空分布變化

    中街山列島海域海洋性氣候明顯,春季在季風(fēng)的影響下,東海表層水進(jìn)入浙江近岸,與陸地徑流(長(zhǎng)江、錢塘江)交匯,形成鋒區(qū),漁場(chǎng)近岸水溫較高,向東部水溫逐步遞減,近岸水與外側(cè)水鹽度差距較大。秋季季風(fēng)轉(zhuǎn)換,偏北風(fēng)頻率的增加,漁場(chǎng)內(nèi)受冷暖水團(tuán)交匯影響較強(qiáng),由近岸往東水溫逐步升高,鹽度自西向東逐步增大[15]。曼氏無針烏賊屬沿岸暖水底棲種,生命周期僅有1a,春季由外海低溫高鹽的越冬場(chǎng)游向近岸島礁附近高溫低鹽的產(chǎn)卵場(chǎng)繁殖產(chǎn)卵,秋季隨季風(fēng)轉(zhuǎn)換,開始向外海高鹽海域移動(dòng),進(jìn)行索餌、越冬洄游[16-19]。研究結(jié)果顯示,春季(4 月)曼氏無針烏賊主要分布于島礁外側(cè)的東北部,秋季(10 月)曼氏無針烏賊主要分布于島礁內(nèi)部偏南方向,這與倪正雅等[17]和李成林等[20]的研究結(jié)果相一致。故不同季節(jié)曼氏無針烏賊的時(shí)空分布與海水溫度、海水鹽度和洄游路徑有關(guān)。

    本研究中,秋季曼氏無針烏賊的尾數(shù)密度明顯大于春季,但春季的生物量卻遠(yuǎn)高于秋季。春季為曼氏無針烏賊的繁殖期,主要為個(gè)體較大的曼氏無針烏賊產(chǎn)卵親體,而秋季為曼氏無針烏賊的補(bǔ)充群體,故春季烏賊的體質(zhì)量要高于秋季。而造成尾數(shù)密度秋季大于春季的原因可能與增殖放流和海洋保護(hù)區(qū)的資源養(yǎng)護(hù)作用有關(guān)。從2013 年開始,在中街山列島海洋保護(hù)區(qū)內(nèi)實(shí)施了大規(guī)模曼氏無針烏賊增殖放流,累計(jì)放流受精卵達(dá)7 835×104ind.,一定程度上增加了曼氏無針烏賊的補(bǔ)充群體[21]。中街山列島國(guó)家級(jí)海洋特別保護(hù)區(qū)的設(shè)立有利于海洋生物的資源保護(hù)和修復(fù),烏賊增殖礁的投放,也有利于自然海域曼氏無針烏賊受精卵的附著,提高了曼氏無針烏賊的附卵量,增加補(bǔ)充群體[22]。

    3.2 棲息地選擇的關(guān)鍵影響因子

    不同生活史階段對(duì)于生境的需求也有所不同,模型擬合結(jié)果顯示春季曼氏無針烏賊的分布受水深和葉綠素a2 個(gè)關(guān)鍵因子的影響,秋季曼氏無針烏賊的分布受單一因子溶解氧的影響。

    曼氏無針烏賊主要分布于近海淺水海域,其產(chǎn)卵場(chǎng)位于島礁周圍50 m 范圍內(nèi),在繁殖期時(shí)從60~80 m的越冬場(chǎng)游向10~40 m 的產(chǎn)卵場(chǎng)作生殖洄游[3]。李星頡等[23]的研究表明,中街山列島海域曼氏無針烏賊的附卵水深范圍在10~20 m。本研究通過GAM 模型顯示春季曼氏無針烏賊會(huì)選擇水深范圍在25~35 m 的海域作為產(chǎn)卵場(chǎng)和育幼場(chǎng),在水深為30 m 時(shí),曼氏無針烏賊的資源密度達(dá)到最高值,原因可能是本研究基于二項(xiàng)分布的GAM 模型來表達(dá)曼氏無針烏賊的出現(xiàn)概率,從已知數(shù)據(jù)預(yù)測(cè)整體分布特點(diǎn),推斷曼氏無針烏賊出現(xiàn)理論分布區(qū),而前人研究以產(chǎn)卵附著基(柳珊瑚)上受精卵的附著率為調(diào)查內(nèi)容,故本研究與前人研究結(jié)果略有差異。

    葉綠素a的濃度值反映了海域初級(jí)生產(chǎn)者的現(xiàn)存量,能直接反映海區(qū)的肥癟程度,是養(yǎng)育生物資源的直接指標(biāo),也能間接反映海區(qū)初級(jí)生產(chǎn)力的高低[24-26]。葉綠素a作為浮游植物豐富度的重要指標(biāo),能為浮游生物提供有利的生長(zhǎng)條件[27]。葉綠素a濃度與頭足類的漁場(chǎng)分布有一定關(guān)聯(lián)性,如劍尖槍烏賊Uroteuthis edulis、神戶烏賊Sepia kobiensis、金烏賊Sepia esculenta、長(zhǎng)蛸Octopus variabilis等漁場(chǎng)分布均受到葉綠素a濃度值的影響[28]。根據(jù)GAM 模型分析的結(jié)果,葉綠素a對(duì)曼氏無針烏賊的分布影響較大,隨著葉綠素a濃度的增加曼氏無針烏賊的資源密度也呈上升趨勢(shì),浮游生物的增多,可以為烏賊提供充足的餌料[29]。成體曼氏無針烏賊主要以小型魚類、甲殼動(dòng)物和少量軟體動(dòng)物為主,幼烏賊主要以浮游動(dòng)物,如箭蟲、橈足類等,特別是夜光藻可作為烏賊漁汛的直觀指標(biāo)[23,30]。

    海水中充裕的溶解氧是維持海洋異養(yǎng)生物生命活動(dòng)、避免水體有機(jī)物發(fā)生厭氧分解而散發(fā)臭味的必要條件,故海水中溶解氧被作為保障海洋生物生存的水化學(xué)指標(biāo)和水體的感官指標(biāo)[31]。李星頡等[23]對(duì)中街山列島海域曼氏無針烏賊產(chǎn)卵區(qū)的生態(tài)環(huán)境進(jìn)行調(diào)查,結(jié)果顯示4-6 月自然海域曼氏無針烏賊產(chǎn)卵場(chǎng)的溶解氧為7.0~8.8 mg·L-1。蔣霞敏等[32]和劉振勇等[33]的研究表明曼氏無針烏賊室內(nèi)孵化期的溶解氧為4.3~4.8 mg·L-1,在養(yǎng)殖過程中溶解氧應(yīng)維持在5.0 mg·L-1以上。頭足類的所有活動(dòng)均需要消耗氧氣,且曼氏無針烏賊的耗氧率與水溫成正比例關(guān)系,個(gè)體體重成反比例關(guān)系,在低溶解氧的狀態(tài)下,會(huì)造成曼氏無針烏賊生長(zhǎng)緩慢,使其新陳代謝率降低,影響個(gè)體活動(dòng)[34-35]。本研究結(jié)果顯示,春季中街山列島海域平均溶解氧為7.0 mg·L-1,這與李星頡等[23]的研究結(jié)果一致。模型擬合結(jié)果顯示秋季曼氏無針烏賊偏向棲息于溶解氧濃度在5.5~6.5 mg·L-1之間的海域,同春季相比產(chǎn)卵區(qū)溶氧含量有所降低,蔣國(guó)昌等[36]的研究結(jié)果表明,東海區(qū)海水氧含量季節(jié)變化明顯,春季含量最高,秋季含量最低。故季節(jié)變化同產(chǎn)卵區(qū)溶解氧濃度變化有明顯關(guān)系。此外,近年來曼氏無針烏賊小型化明顯,人工養(yǎng)殖的曼氏無針烏賊個(gè)體普遍偏小,生長(zhǎng)速度快,壽命周期為4 個(gè)月左右,且人工養(yǎng)殖環(huán)境中溶氧可調(diào)控,故自然海域和養(yǎng)殖環(huán)境下曼氏無針烏賊所需溶氧有所不同[37]。

    在構(gòu)建曼氏無針烏賊棲息偏向性模型時(shí),除水深、水溫、鹽度、葉綠素a、溶解氧等環(huán)境因子外,餌料生物與捕食者、天然附著基柳珊瑚與海藻、底質(zhì)類型、流速等均可對(duì)曼氏無針烏賊棲息地的選擇產(chǎn)生影響。因此,今后還需考慮生物與非生物因素共同影響下曼氏無針烏賊棲息地的選擇,以期進(jìn)一步提高模型的預(yù)測(cè)能力,全面客觀地分析曼氏無針烏賊棲息地選擇偏好與分布特征。

    猜你喜歡
    環(huán)境模型
    一半模型
    長(zhǎng)期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    孕期遠(yuǎn)離容易致畸的環(huán)境
    不能改變環(huán)境,那就改變心境
    環(huán)境
    孕期遠(yuǎn)離容易致畸的環(huán)境
    3D打印中的模型分割與打包
    国产成人精品在线电影| 日日爽夜夜爽网站| 少妇 在线观看| 精品亚洲成国产av| videosex国产| 99香蕉大伊视频| 欧美少妇被猛烈插入视频| 一级a做视频免费观看| 99九九在线精品视频| 国产精品蜜桃在线观看| 成年人免费黄色播放视频| 一级片免费观看大全| 国产欧美另类精品又又久久亚洲欧美| 欧美国产精品va在线观看不卡| 国产在线一区二区三区精| 久久精品国产亚洲av涩爱| 18在线观看网站| xxx大片免费视频| 一本久久精品| 丁香六月天网| 亚洲一码二码三码区别大吗| 亚洲欧洲日产国产| 中文乱码字字幕精品一区二区三区| av在线app专区| videosex国产| 少妇人妻久久综合中文| 国产精品久久久久久av不卡| 中文字幕免费在线视频6| 国产色婷婷99| 中文天堂在线官网| 日韩制服骚丝袜av| 性色avwww在线观看| 成年av动漫网址| 91精品三级在线观看| 亚洲精品久久久久久婷婷小说| 日韩av免费高清视频| 2021少妇久久久久久久久久久| 97在线人人人人妻| 亚洲av日韩在线播放| 亚洲av男天堂| 2021少妇久久久久久久久久久| videos熟女内射| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 亚洲精品国产av成人精品| 高清在线视频一区二区三区| 亚洲精品美女久久久久99蜜臀 | 精品福利永久在线观看| 男女下面插进去视频免费观看 | 久久久精品区二区三区| 久久人人97超碰香蕉20202| 亚洲三级黄色毛片| 人人妻人人爽人人添夜夜欢视频| 久久精品人人爽人人爽视色| 中文天堂在线官网| 大香蕉久久成人网| 亚洲av.av天堂| 亚洲国产精品一区三区| 中文字幕亚洲精品专区| 亚洲国产成人一精品久久久| 又黄又爽又刺激的免费视频.| 狠狠精品人妻久久久久久综合| 91成人精品电影| 人妻一区二区av| 久久99蜜桃精品久久| 欧美性感艳星| 18禁观看日本| 中文字幕人妻丝袜制服| 国产成人免费观看mmmm| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 美女大奶头黄色视频| 一边亲一边摸免费视频| 建设人人有责人人尽责人人享有的| 亚洲欧美成人精品一区二区| 中文字幕人妻熟女乱码| 99久久综合免费| 纯流量卡能插随身wifi吗| 看免费成人av毛片| 少妇人妻久久综合中文| 欧美精品亚洲一区二区| 欧美+日韩+精品| 五月伊人婷婷丁香| 777米奇影视久久| 人人澡人人妻人| 夜夜骑夜夜射夜夜干| 老熟女久久久| 成人亚洲欧美一区二区av| 美女国产视频在线观看| 久久精品国产综合久久久 | 一级片'在线观看视频| 老司机影院成人| 青春草视频在线免费观看| 国产精品 国内视频| 久久ye,这里只有精品| 99热这里只有是精品在线观看| 精品一品国产午夜福利视频| 日本av免费视频播放| xxxhd国产人妻xxx| 国产爽快片一区二区三区| 看十八女毛片水多多多| 最近最新中文字幕免费大全7| 婷婷色综合www| 精品视频人人做人人爽| 国产精品一国产av| 亚洲国产最新在线播放| 成年女人在线观看亚洲视频| 中文欧美无线码| 国产综合精华液| 如何舔出高潮| 熟女电影av网| 伦理电影免费视频| 超碰97精品在线观看| 高清视频免费观看一区二区| 菩萨蛮人人尽说江南好唐韦庄| 最后的刺客免费高清国语| 又黄又爽又刺激的免费视频.| 国产女主播在线喷水免费视频网站| 久久久精品94久久精品| 久久久久久人妻| 欧美日韩综合久久久久久| 18禁裸乳无遮挡动漫免费视频| 美女国产高潮福利片在线看| 亚洲av在线观看美女高潮| 老司机影院毛片| 亚洲一级一片aⅴ在线观看| 色婷婷久久久亚洲欧美| 十八禁网站网址无遮挡| 中国三级夫妇交换| 日韩一区二区三区影片| 国产av精品麻豆| 国产国拍精品亚洲av在线观看| 最近的中文字幕免费完整| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品无大码| 日日爽夜夜爽网站| 青春草视频在线免费观看| 日产精品乱码卡一卡2卡三| 男女高潮啪啪啪动态图| 国产深夜福利视频在线观看| 一级片免费观看大全| 一级毛片电影观看| 老熟女久久久| 亚洲精品久久久久久婷婷小说| 亚洲,一卡二卡三卡| av在线app专区| 欧美日韩一区二区视频在线观看视频在线| 久久久久人妻精品一区果冻| 新久久久久国产一级毛片| 女人久久www免费人成看片| 精品国产乱码久久久久久小说| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 成人国产av品久久久| 国内精品宾馆在线| 免费黄网站久久成人精品| 女的被弄到高潮叫床怎么办| 青春草视频在线免费观看| 国产免费现黄频在线看| 中文欧美无线码| 国产激情久久老熟女| 一级毛片我不卡| 如何舔出高潮| 亚洲成av片中文字幕在线观看 | 精品熟女少妇av免费看| 国产欧美日韩综合在线一区二区| av女优亚洲男人天堂| 欧美+日韩+精品| 国产成人免费观看mmmm| 蜜桃国产av成人99| 久久av网站| 日韩人妻精品一区2区三区| 激情五月婷婷亚洲| 在线亚洲精品国产二区图片欧美| 男女啪啪激烈高潮av片| 大片免费播放器 马上看| 免费黄频网站在线观看国产| 亚洲美女黄色视频免费看| 咕卡用的链子| 精品少妇黑人巨大在线播放| 国产精品人妻久久久影院| a 毛片基地| 久久久久精品人妻al黑| 夫妻午夜视频| 日本91视频免费播放| 男女啪啪激烈高潮av片| 午夜福利视频精品| 国产精品人妻久久久久久| 日韩精品有码人妻一区| 不卡视频在线观看欧美| 女人精品久久久久毛片| 欧美亚洲日本最大视频资源| 男的添女的下面高潮视频| 久久精品久久精品一区二区三区| 成年女人在线观看亚洲视频| 色婷婷av一区二区三区视频| 99热这里只有是精品在线观看| 久久久久久久久久人人人人人人| 久久热在线av| 97在线人人人人妻| 观看av在线不卡| 在线 av 中文字幕| 99国产精品免费福利视频| 18禁裸乳无遮挡动漫免费视频| a级毛色黄片| 99国产综合亚洲精品| 激情五月婷婷亚洲| 免费少妇av软件| 日本av免费视频播放| 高清不卡的av网站| 宅男免费午夜| 高清av免费在线| 国产成人欧美| 亚洲精品日韩在线中文字幕| 精品国产国语对白av| 欧美变态另类bdsm刘玥| 色吧在线观看| 在线观看www视频免费| 国产精品久久久久久久电影| 亚洲,一卡二卡三卡| 成年动漫av网址| 免费在线观看黄色视频的| 熟女人妻精品中文字幕| 免费观看在线日韩| 91精品伊人久久大香线蕉| 一区二区日韩欧美中文字幕 | 麻豆精品久久久久久蜜桃| 国产高清不卡午夜福利| 亚洲av电影在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 日本色播在线视频| 久久精品夜色国产| 寂寞人妻少妇视频99o| 人人妻人人添人人爽欧美一区卜| 少妇的逼水好多| 亚洲国产av影院在线观看| 久久久久国产精品人妻一区二区| 制服丝袜香蕉在线| 精品少妇内射三级| 欧美激情国产日韩精品一区| 大码成人一级视频| 新久久久久国产一级毛片| 亚洲美女黄色视频免费看| 免费日韩欧美在线观看| 免费av不卡在线播放| 中文字幕最新亚洲高清| 男女午夜视频在线观看 | 亚洲国产精品国产精品| av免费观看日本| 国产精品一区www在线观看| 色94色欧美一区二区| 青春草亚洲视频在线观看| 久久久国产欧美日韩av| 日本wwww免费看| 黄网站色视频无遮挡免费观看| 亚洲欧美一区二区三区黑人 | 高清毛片免费看| 汤姆久久久久久久影院中文字幕| 亚洲国产精品成人久久小说| 中文字幕另类日韩欧美亚洲嫩草| 男男h啪啪无遮挡| 国产精品人妻久久久久久| 一本—道久久a久久精品蜜桃钙片| 黄色 视频免费看| 亚洲av欧美aⅴ国产| 午夜91福利影院| 亚洲成av片中文字幕在线观看 | 啦啦啦中文免费视频观看日本| 在线观看免费视频网站a站| 又大又黄又爽视频免费| 欧美日韩亚洲高清精品| 日本与韩国留学比较| 日韩欧美精品免费久久| 性高湖久久久久久久久免费观看| 人人妻人人添人人爽欧美一区卜| 日本欧美视频一区| 最近最新中文字幕免费大全7| 欧美成人精品欧美一级黄| 日韩av不卡免费在线播放| 丰满少妇做爰视频| 伦理电影大哥的女人| 青春草亚洲视频在线观看| 亚洲国产看品久久| 成人国语在线视频| 男女午夜视频在线观看 | 日韩视频在线欧美| 中国三级夫妇交换| 乱码一卡2卡4卡精品| 日韩伦理黄色片| a级毛色黄片| 交换朋友夫妻互换小说| 亚洲精品日本国产第一区| 久久女婷五月综合色啪小说| 久久 成人 亚洲| 高清视频免费观看一区二区| 美女xxoo啪啪120秒动态图| 日韩电影二区| 国产成人精品婷婷| 欧美精品高潮呻吟av久久| 国产极品粉嫩免费观看在线| 国产欧美日韩一区二区三区在线| 美女主播在线视频| 丝袜在线中文字幕| 国产精品无大码| 亚洲五月色婷婷综合| 制服丝袜香蕉在线| 99久国产av精品国产电影| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91| 熟妇人妻不卡中文字幕| 最黄视频免费看| 人妻系列 视频| 亚洲一码二码三码区别大吗| 国产老妇伦熟女老妇高清| 国产精品久久久久久久电影| 婷婷色综合www| 亚洲av电影在线进入| 啦啦啦中文免费视频观看日本| 18禁国产床啪视频网站| 狂野欧美激情性xxxx在线观看| 在线观看免费视频网站a站| 精品久久蜜臀av无| 巨乳人妻的诱惑在线观看| 亚洲av男天堂| 一个人免费看片子| 我的女老师完整版在线观看| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 国产色爽女视频免费观看| 少妇人妻久久综合中文| 日韩制服丝袜自拍偷拍| 精品国产一区二区久久| 内地一区二区视频在线| 亚洲国产av影院在线观看| 中文字幕制服av| 香蕉国产在线看| 久久久精品区二区三区| 如日韩欧美国产精品一区二区三区| 人妻系列 视频| 日韩av在线免费看完整版不卡| 国产片特级美女逼逼视频| 久久人妻熟女aⅴ| 人妻人人澡人人爽人人| 精品熟女少妇av免费看| 久久久久网色| 一级片免费观看大全| av女优亚洲男人天堂| 久久99一区二区三区| 高清欧美精品videossex| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 人妻 亚洲 视频| 国产激情久久老熟女| 亚洲第一av免费看| 国产男人的电影天堂91| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 青春草视频在线免费观看| 中文精品一卡2卡3卡4更新| 午夜影院在线不卡| 国产激情久久老熟女| 日韩精品免费视频一区二区三区 | 日日啪夜夜爽| 我的女老师完整版在线观看| 在线观看一区二区三区激情| 欧美精品国产亚洲| 两个人免费观看高清视频| 2018国产大陆天天弄谢| 午夜福利乱码中文字幕| 亚洲精品国产色婷婷电影| 午夜免费男女啪啪视频观看| 国产男女内射视频| 久久婷婷青草| 国产精品一区二区在线不卡| 人人妻人人爽人人添夜夜欢视频| 一区二区三区四区激情视频| 波野结衣二区三区在线| 我要看黄色一级片免费的| 在线天堂最新版资源| av电影中文网址| 国产有黄有色有爽视频| 高清欧美精品videossex| 永久网站在线| 9热在线视频观看99| 亚洲精品aⅴ在线观看| 黄色视频在线播放观看不卡| 欧美亚洲 丝袜 人妻 在线| 久久久精品区二区三区| 久久精品国产自在天天线| 中文字幕免费在线视频6| 一区二区av电影网| 国产日韩欧美亚洲二区| 精品午夜福利在线看| 成人综合一区亚洲| 欧美另类一区| 亚洲国产精品999| 校园人妻丝袜中文字幕| 成人亚洲精品一区在线观看| 一本久久精品| 精品福利永久在线观看| 51国产日韩欧美| 丝袜人妻中文字幕| 欧美 日韩 精品 国产| 国产欧美日韩综合在线一区二区| 国产亚洲精品第一综合不卡 | 国产伦理片在线播放av一区| 国产av国产精品国产| 久久久久久久久久久免费av| 欧美人与性动交α欧美软件 | 亚洲精品456在线播放app| 免费观看a级毛片全部| 中文字幕制服av| 国产一区二区在线观看av| 十八禁高潮呻吟视频| 各种免费的搞黄视频| 韩国精品一区二区三区 | 欧美3d第一页| 男女免费视频国产| 美女福利国产在线| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 麻豆乱淫一区二区| 肉色欧美久久久久久久蜜桃| 欧美精品av麻豆av| 亚洲第一av免费看| 人人妻人人澡人人看| 国产 精品1| 中文字幕另类日韩欧美亚洲嫩草| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| 男女啪啪激烈高潮av片| 一区二区三区乱码不卡18| 欧美激情国产日韩精品一区| 久久人人爽人人爽人人片va| 久久精品国产鲁丝片午夜精品| 在线观看免费日韩欧美大片| 国产黄色视频一区二区在线观看| 99精国产麻豆久久婷婷| 欧美国产精品va在线观看不卡| 久久久亚洲精品成人影院| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 国产亚洲精品第一综合不卡 | 天堂8中文在线网| 99re6热这里在线精品视频| 美女福利国产在线| 在线看a的网站| 少妇 在线观看| 男人操女人黄网站| 精品亚洲乱码少妇综合久久| 久久人人爽人人爽人人片va| 久热这里只有精品99| 国产综合精华液| 色94色欧美一区二区| 2021少妇久久久久久久久久久| 高清不卡的av网站| 99视频精品全部免费 在线| 中文乱码字字幕精品一区二区三区| 18在线观看网站| 丝袜在线中文字幕| 中国三级夫妇交换| kizo精华| 亚洲婷婷狠狠爱综合网| 成人毛片a级毛片在线播放| 国产日韩欧美视频二区| 搡女人真爽免费视频火全软件| 在线观看一区二区三区激情| 国产日韩一区二区三区精品不卡| 90打野战视频偷拍视频| 久久久久精品性色| 一区在线观看完整版| 97在线人人人人妻| www.熟女人妻精品国产 | 日本欧美国产在线视频| 只有这里有精品99| 日日啪夜夜爽| 激情五月婷婷亚洲| 国产精品.久久久| 又黄又粗又硬又大视频| 日产精品乱码卡一卡2卡三| 天天躁夜夜躁狠狠躁躁| 大码成人一级视频| 久久久久久久久久成人| 亚洲综合精品二区| 这个男人来自地球电影免费观看 | 丰满饥渴人妻一区二区三| 日韩精品有码人妻一区| 精品视频人人做人人爽| 日本爱情动作片www.在线观看| 精品国产露脸久久av麻豆| 亚洲人与动物交配视频| 少妇被粗大猛烈的视频| 精品国产一区二区三区四区第35| 人人妻人人添人人爽欧美一区卜| 全区人妻精品视频| 久久99热6这里只有精品| 精品一区二区免费观看| 日韩欧美一区视频在线观看| 一区二区三区四区激情视频| 亚洲av日韩在线播放| 妹子高潮喷水视频| 亚洲成国产人片在线观看| 欧美老熟妇乱子伦牲交| 9191精品国产免费久久| 亚洲av国产av综合av卡| 亚洲精品一区蜜桃| 国产精品欧美亚洲77777| 九草在线视频观看| av女优亚洲男人天堂| 国产熟女午夜一区二区三区| 国产xxxxx性猛交| 午夜免费鲁丝| 天堂中文最新版在线下载| 国产精品国产三级国产专区5o| 国产麻豆69| 国产亚洲精品第一综合不卡 | 97在线人人人人妻| 久久久久久伊人网av| 久热这里只有精品99| 国产成人免费无遮挡视频| 纵有疾风起免费观看全集完整版| 飞空精品影院首页| 2018国产大陆天天弄谢| 欧美激情极品国产一区二区三区 | 制服诱惑二区| 久久国产精品男人的天堂亚洲 | 狠狠婷婷综合久久久久久88av| 久久人人爽人人片av| 亚洲综合色惰| 欧美激情 高清一区二区三区| 性色avwww在线观看| 欧美日韩成人在线一区二区| 涩涩av久久男人的天堂| 久久综合国产亚洲精品| 国产精品久久久久久久电影| 热re99久久精品国产66热6| 女人精品久久久久毛片| 人人妻人人爽人人添夜夜欢视频| 久久国内精品自在自线图片| 日本色播在线视频| 日韩电影二区| 国产成人a∨麻豆精品| 男女啪啪激烈高潮av片| 亚洲精品一区蜜桃| 在线观看国产h片| 插逼视频在线观看| 超碰97精品在线观看| 亚洲av福利一区| 黄色配什么色好看| 亚洲人成网站在线观看播放| 亚洲成人手机| 1024视频免费在线观看| 男女午夜视频在线观看 | 人体艺术视频欧美日本| 国产日韩欧美视频二区| 国产精品国产三级国产av玫瑰| 久久久久久久久久成人| 成年人午夜在线观看视频| 免费在线观看黄色视频的| 精品人妻在线不人妻| 91精品国产国语对白视频| 91午夜精品亚洲一区二区三区| 亚洲国产精品国产精品| 国产精品女同一区二区软件| 亚洲国产精品一区二区三区在线| 亚洲国产成人一精品久久久| 97在线视频观看| 99精国产麻豆久久婷婷| 婷婷色av中文字幕| 免费久久久久久久精品成人欧美视频 | 欧美激情极品国产一区二区三区 | 国产精品无大码| 精品一区二区三卡| av在线播放精品| 天天躁夜夜躁狠狠久久av| 丝瓜视频免费看黄片| 亚洲av电影在线观看一区二区三区| 日本-黄色视频高清免费观看| 欧美成人精品欧美一级黄| 国产亚洲午夜精品一区二区久久| 亚洲成人av在线免费| 国产精品久久久久久精品古装| 2018国产大陆天天弄谢| 亚洲人成77777在线视频| 建设人人有责人人尽责人人享有的| 日韩中文字幕视频在线看片| 九草在线视频观看| 一级黄片播放器| 国产精品人妻久久久影院| 女人久久www免费人成看片| 少妇熟女欧美另类| 五月玫瑰六月丁香| 久久久久久久久久久免费av| 青春草国产在线视频| 国产成人精品婷婷| 亚洲人成77777在线视频| 妹子高潮喷水视频| 日本猛色少妇xxxxx猛交久久| 色网站视频免费| 蜜桃国产av成人99| 97在线视频观看| 肉色欧美久久久久久久蜜桃| 国产欧美日韩一区二区三区在线| 搡老乐熟女国产| 寂寞人妻少妇视频99o| 日韩制服丝袜自拍偷拍| 久久人人爽av亚洲精品天堂| 男的添女的下面高潮视频| 国产亚洲精品第一综合不卡 | 满18在线观看网站| 亚洲av在线观看美女高潮| 黄色怎么调成土黄色|