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

    甌江口鳳鱭魚卵、仔稚魚的時(shí)空分布及其與環(huán)境因子關(guān)系

    2022-08-04 15:15:38孫浩奇蔣日進(jìn)陳峰芮銀李霞芳李振華印瑞周永東
    水生生物學(xué)報(bào) 2022年7期
    關(guān)鍵詞:稚魚甌江魚卵

    孫浩奇 蔣日進(jìn) 陳峰 芮銀 李霞芳 李振華 印瑞 周永東

    (1.浙江海洋大學(xué)海洋與漁業(yè)研究所,舟山 316021;2.浙江省海洋水產(chǎn)研究所農(nóng)業(yè)農(nóng)村部重點(diǎn)漁場(chǎng)漁業(yè)資源科學(xué)觀測(cè)實(shí)驗(yàn)站,浙江省海洋漁業(yè)資源可持續(xù)利用技術(shù)研究重點(diǎn)實(shí)驗(yàn)室,舟山 316021)

    鳳鱭(Coilia mystus)隸屬鯡形目(Clupeiformes),鳀科(Engraulidae),鱭屬(Coilia)[1],在我國(guó)東南沿海均有分布,根據(jù)形態(tài)特征差異可大致分為長(zhǎng)江型、閩江型和珠江型三種生態(tài)類群[2]。鳳鱭作為一種短距離溯河洄游性魚類[3],一般在近海分散生活,繁殖期時(shí)成熟的鳳鱭個(gè)體集群洄游至河口進(jìn)行產(chǎn)卵,繁殖活動(dòng)結(jié)束后親體重新回到近海生活[4]。甌江入海河口(以下簡(jiǎn)稱“甌江口”)是鳳鱭主要棲息場(chǎng)所之一,位于浙江南部近岸,是甌江徑流和臺(tái)灣暖流等多水團(tuán)共同影響下的咸淡水交匯區(qū)[5]。甌江口水域多紊流,潮流較為復(fù)雜,其徑流帶來豐富的營(yíng)養(yǎng)鹽和餌料生物,給魚類提供了多元化的棲息地類型[6]。魚卵仔稚魚等對(duì)于棲息環(huán)境的變化較為敏感,外部環(huán)境的微小變化就可能對(duì)早期群體的生存和延續(xù)產(chǎn)生較大的影響。鳳鱭作為甌江口海域主要優(yōu)勢(shì)種,其數(shù)量變化在一定程度上反映出洄游性魚類群落結(jié)構(gòu)對(duì)環(huán)境適應(yīng)的波動(dòng)性[7]。

    基于統(tǒng)計(jì)回歸分析的廣義加性模型(Generalized Additive Model,GAM)是一種分析漁業(yè)資源數(shù)量分布與環(huán)境因子間關(guān)系的模型方法[8,9]。在實(shí)際的資源調(diào)查中,早期群體資源豐度等數(shù)據(jù)可能存在大量0值,資源分布方式不均勻,很難服從某一特定的數(shù)據(jù)分布方式,若采用常規(guī)GAM模型擬合早期群體豐度變化與環(huán)境因子間關(guān)系,可能導(dǎo)致模型擬合結(jié)果與實(shí)際情況誤差較大[10]。兩階段廣義加性模型(Two-stage GAM)通過分階段處理0值數(shù)據(jù),可以較好地看出豐度變化與環(huán)境之間的關(guān)系[11],其在生物群體與環(huán)境關(guān)系的研究中已得到相關(guān)應(yīng)用[12],但在魚卵仔稚魚等早期群體中的研究較少。本研究根據(jù)2015年、2018年和2019年的5—8月在甌江口及其鄰近海域進(jìn)行的鳳鱭資源和環(huán)境調(diào)查數(shù)據(jù),分析鳳鱭早期群體的資源分布,應(yīng)用Two-stage GAM模型解析鳳鱭魚卵、仔稚魚與環(huán)境因子之間的關(guān)系,探討非生物因子對(duì)鳳鱭早期發(fā)育階段生境選擇的影響[13],為探明甌江口鳳鱭產(chǎn)卵場(chǎng)和育幼場(chǎng)提供基礎(chǔ)資料[14],同時(shí)基于2020年5—7月的環(huán)境數(shù)據(jù)進(jìn)行模型預(yù)測(cè),以期為合理預(yù)測(cè)甌江口鳳鱭資源的時(shí)空變動(dòng)提供科學(xué)參考依據(jù)。

    1 材料與方法

    1.1 調(diào)查區(qū)域與采樣方法

    鳳鱭魚卵和仔稚魚資源和環(huán)境數(shù)據(jù)來源于2015年、2018年、2019年5至8月在甌江口內(nèi)進(jìn)行的產(chǎn)卵場(chǎng)選劃項(xiàng)目。共設(shè)置11個(gè)站位,每月的大潮汛和小潮汛之間,在靈昆島至橋下鎮(zhèn)之間水域(27°53′N—28°51′N、120°35′E—120°58′E)調(diào)查采集鳳鱭的魚卵和仔稚魚(圖1),使用淺水I型浮游生物網(wǎng)(網(wǎng)口直徑0.5 m,網(wǎng)目尺寸0.505 mm)進(jìn)行水平拖網(wǎng)采樣。調(diào)查船租用群眾雙拖漁船,單機(jī)功率為56 kW,調(diào)查拖速為2 kn/h,每站拖網(wǎng)時(shí)間10min。每一個(gè)調(diào)查站位的水深、水溫和水鹽等環(huán)境數(shù)據(jù)通過多功能水質(zhì)儀同步測(cè)定。依據(jù)《海洋調(diào)查規(guī)范第六部分-海洋生物調(diào)查》采集、保存樣品,在實(shí)驗(yàn)室根據(jù)形態(tài)特征和分子鑒定技術(shù)對(duì)鳳鱭魚卵和仔稚魚進(jìn)行種類鑒定[15]。根據(jù)各個(gè)站位鳳鱭魚卵和鳳鱭仔稚魚的數(shù)量及濾水量,計(jì)算每個(gè)站位鳳鱭魚卵和仔稚魚的豐度(ind./m3)。

    圖1 甌江口鳳鱭魚卵仔稚魚采樣位點(diǎn)Fig.1 Sampling stations of Coilia mystus eggs and larvae in the Oujiang Estuary

    1.2 初始影響因子的篩選

    產(chǎn)卵場(chǎng)環(huán)境對(duì)于魚類早期群體分布具有顯著影響[16],根據(jù)鳳鱭產(chǎn)卵孵化等相關(guān)研究,在構(gòu)建GAM模型的初始解釋變量中,選取水深(Depth,DEP)、表層水溫(Surface water temperature,ST)、表層水鹽(Surface water salinity,SS)、底層水溫(Bottom water temperature,BT)和底層水鹽(Bottom water salinity,BS)5個(gè)環(huán)境變量,經(jīng)度(Lat)和緯度(Lon)2個(gè)空間變量以及月份(Month)作為時(shí)間變量加入模型進(jìn)行建模。預(yù)測(cè)變量的多重共線性會(huì)對(duì)模型參數(shù)回歸的準(zhǔn)確性產(chǎn)生干擾,本研究通過Pearson相關(guān)分析和方差膨脹因子對(duì)環(huán)境變量及空間變量進(jìn)行共線性診斷。Pearson相關(guān)性閾值設(shè)定為0.7,方差膨脹系數(shù)(VIF)閾值設(shè)定為4,認(rèn)為相關(guān)性和方差膨脹因子高于閾值的變量存在嚴(yán)重的多重共線性[17];通過pearson相關(guān)性分析找出存在多重共線性的變量,去除其中VIF值最大的變量,重復(fù)此步驟至所有預(yù)測(cè)變量的VIF值和Pearson相關(guān)性值均小于閾值。

    1.3 數(shù)據(jù)處理與分析

    棲息地環(huán)境對(duì)生物的分布有十分顯著的影響。為進(jìn)一步了解甌江口鳳鱭產(chǎn)卵場(chǎng)和育幼場(chǎng)的分布,本研究利用兩階段廣義加性模型(Two-stage GAM)分別對(duì)鳳鱭魚卵、仔稚魚的分布與環(huán)境因子的關(guān)系進(jìn)行分析。模型構(gòu)建包含兩部分:第一部分,基于二項(xiàng)分布構(gòu)建GAM1模型,使用logistic鏈接函數(shù)模擬鳳鱭魚卵、仔稚魚出現(xiàn)的概率P與環(huán)境因子的關(guān)系;第二部分,基于高斯分布構(gòu)建GAM2模型,使用identical鏈接函數(shù)模擬鳳鱭魚卵、仔稚魚出現(xiàn)情況下其豐度與環(huán)境因子的關(guān)系,使用對(duì)數(shù)函數(shù)對(duì)豐度進(jìn)行標(biāo)準(zhǔn)化處理[11]。

    式中,P為鳳鱭魚卵、仔稚魚出現(xiàn)概率,Y為鳳鱭魚卵、仔稚魚資源豐度,α為函數(shù)的截距,xi為解釋變量,Si(xi)為各解釋變量的單變量樣條平滑函數(shù),表示殘差。

    根據(jù)赤池信息準(zhǔn)則(Akaike Information Criterion,AIC)和解釋率篩選最優(yōu)模型組合。在模型擬合過程中,采用逐步回歸法依次將解釋變量加入至模型中,至模型的AIC值最小且不再變化,選擇解釋率最高的模型作為最終模型[18]。一般認(rèn)為AIC值越低,偏差解釋率越高,模型擬合效果越好。

    式中,k是參數(shù)個(gè)數(shù),L是似然函數(shù)。

    1.4 模型驗(yàn)證與預(yù)測(cè)

    采用五折交叉驗(yàn)證分別對(duì)GAM1和GAM2的模型預(yù)測(cè)性能和準(zhǔn)確性進(jìn)行評(píng)價(jià)[19]。將總數(shù)據(jù)集隨機(jī)等分為5個(gè)子數(shù)據(jù)集,每次隨機(jī)選擇4個(gè)子數(shù)據(jù)集作為訓(xùn)練集,剩余的1個(gè)子數(shù)據(jù)集作為驗(yàn)證集用于模型預(yù)測(cè)準(zhǔn)確性評(píng)估,重復(fù)計(jì)算100次,對(duì)各模型的準(zhǔn)確性評(píng)估取平均效應(yīng)。

    GAM1模型的預(yù)測(cè)性能評(píng)價(jià)采用受試者工作特征曲下面積(The area under the receiver operating characteristic curve,AUC)。AUC作為靈敏度和特異度在不同閾值下的綜合指標(biāo),不受二項(xiàng)分布類模型預(yù)測(cè)所需閾值的影響,其閾值獨(dú)立性不需要設(shè)定閾值就能對(duì)模型作出判斷,可以較好地反映GAM1模型的預(yù)測(cè)準(zhǔn)確性。AUC變化范圍為(0,1),一般認(rèn)為,AUC值為 0.5—0.7時(shí)模型價(jià)值較低;0.7—0.9時(shí)模型預(yù)測(cè)力較好;>0.9時(shí)預(yù)測(cè)力極好;<0.5模型擬合差于隨機(jī)模型。根據(jù)模型預(yù)測(cè)力,AUC 均值大于0.7的模型可以用于后續(xù)的GAM2模型的構(gòu)建和分析[20]。

    體育(physical education)除去外顯的生物性外,作為一項(xiàng)社會(huì)活動(dòng),其發(fā)展歷程與人類社會(huì)的變遷相伴相生。個(gè)體的實(shí)踐要通過組織化的活動(dòng)才能構(gòu)成具有更豐富社會(huì)意義的行為,民族國(guó)家的興起對(duì)現(xiàn)代體育的發(fā)生與發(fā)展產(chǎn)生了深刻的影響。在競(jìng)技體育之外,民眾自發(fā)地圍繞運(yùn)動(dòng)、健身等展開的社會(huì)活動(dòng)也成為體育社會(huì)學(xué)研究的主要內(nèi)容之一。

    通過GAM1計(jì)算各站位點(diǎn)的出現(xiàn)概率,在0—0.9設(shè)定出現(xiàn)概率的閾值(從0開始以0.1間隔遞增),當(dāng)出現(xiàn)概率低于閾值時(shí)其相對(duì)資源密度為0,出現(xiàn)概率等于或高于閾值時(shí)采用GAM2計(jì)算其相對(duì)資源豐度,對(duì)于出現(xiàn)概率閾值的選擇及GAM2模型的預(yù)測(cè)性能評(píng)價(jià)采用回歸分析法,即通過構(gòu)建驗(yàn)試集的模型預(yù)測(cè)值(y)與實(shí)際觀測(cè)值(y0)間的線性關(guān)系來描述二者的近似程度:

    a和b的均值反映預(yù)測(cè)性能,式中,斜率a=1時(shí)模型預(yù)測(cè)無系統(tǒng)偏差,截距b=0時(shí)預(yù)測(cè)值與實(shí)際觀測(cè)值具有相似的空間特征。同時(shí)計(jì)算決定系數(shù)(R2)作為回歸方程的擬合優(yōu)度[21]。

    根據(jù)2020年5—7月各站點(diǎn)環(huán)境調(diào)查數(shù)據(jù),以0.05°×0.05°網(wǎng)格大小對(duì)研究水域進(jìn)行劃分,獲得每個(gè)網(wǎng)格中心點(diǎn)的坐標(biāo),采用普通克里金插值法計(jì)算網(wǎng)格中心點(diǎn)的環(huán)境數(shù)據(jù)?;赥wo-stage GAM對(duì)2020年5—7月的相對(duì)資源密度分布進(jìn)行預(yù)測(cè),分析不同月份鳳鱭魚卵、仔稚魚的預(yù)測(cè)結(jié)果與實(shí)際調(diào)查站位的時(shí)空分布。

    利用Excel2016進(jìn)行數(shù)據(jù)整理計(jì)算,R 4.0.3軟件和ArcGIS10.2中進(jìn)行模型構(gòu)建和繪圖,其中GAM模型利用R4.0.3中g(shù)am軟件包完成。

    2 結(jié)果

    2.1 GAM模型構(gòu)建

    Pearson相關(guān)性分析和VIF檢驗(yàn)結(jié)果顯示,底溫、底鹽、緯度與其他解釋變量間存在多重共線性問題。逐步舍去超過閾值的解釋變量后,最終選擇水深、表溫、表鹽、經(jīng)度作為解釋變量,用于魚卵和仔稚魚模型的建立和分析(表1)。在構(gòu)建GAM模型過程中,以時(shí)間變量月份作為分類因子,篩選后的水深等解釋變量作為影響因子進(jìn)行建模。依據(jù)AIC最小原則,構(gòu)建鳳鱭魚卵、仔稚魚的Two-stage GAM最優(yōu)模型,結(jié)果如下:

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

    式中,P為鳳鱭魚卵或仔稚魚出現(xiàn)的概率,Y為鳳鱭魚卵或仔稚魚的相對(duì)資源量,Month為月份,ST為表層水溫,SS為表層水鹽,ε為相應(yīng)的GAM模型截距。

    l.GAM1模型最佳預(yù)測(cè)變量組合為Month+ST+SS,模型累積偏差解釋率為23.21%,表溫對(duì)仔稚魚的出現(xiàn)概率具有極顯著影響(P<0.01),貢獻(xiàn)率為9.31%,月份具有顯著影響(P<0.05),貢獻(xiàn)率為5.43%,表鹽貢獻(xiàn)率為8.74%。l.GAM2模型結(jié)果顯示:影響仔稚魚相對(duì)資源密度的變量為表溫和表鹽,模型的累積偏差解釋率為33.3%,其中表溫對(duì)于模型的貢獻(xiàn)率最大,為28.12%,表鹽貢獻(xiàn)率僅為5.20%。

    2.2 影響因子效應(yīng)分析

    溫州甌江口鳳鱭魚卵、仔稚魚在不同月份其分布有一定差異。根據(jù)5—8月各站位平均相對(duì)資源密度結(jié)果顯示(圖2),鳳鱭魚卵基本分布在江心嶼至七都島附近等低鹽淡水區(qū)水域,其中高豐度區(qū)主要集中于七都島上游。研究結(jié)果顯示鳳鱭主要在甌江口上游產(chǎn)卵,產(chǎn)卵時(shí)間主要集中在5—7月,7月甌江水域出現(xiàn)鳳鱭魚卵的概率有明顯下降(圖3a),8月僅部分水域出現(xiàn)魚卵,考慮到實(shí)際調(diào)查中僅2015年8月采集到兩次鳳鱭,當(dāng)年的鳳鱭資源量較其他年份明顯偏高,同時(shí)2018年、2019年8月均未采集到鳳鱭魚卵,推測(cè)甌江口及附近水域的鳳鱭產(chǎn)卵期至8月結(jié)束。

    鳳鱭仔稚魚5—8月在甌江流域均有分布(圖2),7月仔稚魚出現(xiàn)頻率最高(圖3b)。仔稚魚主要集中分布于江心嶼與靈昆島之間水域,資源重心在七都島北部低鹽區(qū)水域;7月與5月、6月相比,低鹽區(qū)仔稚魚相對(duì)資源量較小,七都島及靈昆島附近水域資源量高于甌江其他水域;8月鳳鱭仔稚魚則圍繞靈昆島高鹽區(qū)海域作密集分布,即甌江入??诟浇?在甌江中上游低鹽區(qū)及七都島高低鹽交匯區(qū)相對(duì)資源量較低,靈昆島北部水域仔稚魚相對(duì)資源量高于南部。

    圖2 鳳鱭魚卵、仔稚魚5—8月平均資源密度分布(a.魚卵;b.仔稚魚)Fig.2 Distribution of average resource density of eggs and larvae of Coilia mystus from May to August (a.eggs;b.larvae)

    圖3 月份對(duì)鳳鱭魚卵仔稚魚出現(xiàn)概率的影響(a.魚卵;b.仔稚魚)Fig.3 Effects of spatiotemporal factors on the eggs and larvae distribution of Coilia mystus (a.eggs;b.larvae)

    鳳鱭作為中上層短距離洄游性魚類,表溫鹽度等水環(huán)境因子是影響鳳鱭產(chǎn)卵繁殖及仔稚魚的生長(zhǎng)、分布最重要的環(huán)境因素(圖4)。

    鳳鱭魚卵主要出現(xiàn)在表溫20—32℃,26℃時(shí)鳳鱭魚卵出現(xiàn)概率最高(圖4a),而魚卵的豐度在20—25℃間呈現(xiàn)先降后升的趨勢(shì),在23℃以后魚卵隨表溫升高而增加(圖4b)。仔稚魚的出現(xiàn)概率隨表溫呈先升后降的趨勢(shì),表溫22℃左右出現(xiàn)率最大(圖4c);豐度隨表溫波動(dòng)變化,高峰區(qū)為24—27℃。綜合來看,魚卵適宜表溫為24—26℃,仔稚魚適宜表溫為22—27℃(圖4d)。

    甌江鳳鱭魚卵和仔稚魚的資源變動(dòng)隨鹽度變化趨勢(shì)明顯,魚卵主要集中分布于低鹽淡水區(qū),在低鹽區(qū)的出現(xiàn)概率最高(圖4e),其豐度隨鹽度升高而顯著降低(圖4f),出現(xiàn)概率和豐度的高峰區(qū)集中在表鹽0—1左右,魚卵的鹽度適宜性在1以下。仔稚魚自低鹽淡水區(qū)至入??谙痰粎R水域均有分布,表鹽為7時(shí)仔稚魚出現(xiàn)概率最高(圖4g),資源豐度峰值在表層鹽度10左右,鳳鱭仔稚魚的出現(xiàn)概率和豐度隨表鹽升高而升高(圖4h),仔稚魚的鹽度適宜性為5—13。

    圖4 環(huán)境因子對(duì)鳳鱭魚卵仔稚魚分布的影響(a、b、e、f .魚卵;c、d、g、h.仔稚魚)Fig.4 Effects of environmental factors on the eggs and larvae distribution of Coilia mystus (a,b,e,f.eggs;c,d,g,h.larvae)

    2.3 模型預(yù)測(cè)性能

    交叉驗(yàn)證的結(jié)果顯示(表2),用于分析魚卵、仔稚魚出現(xiàn)概率的GAM1模型的AUC均值大于0.7,模型預(yù)測(cè)性能較好,其中分析魚卵出現(xiàn)概率的e.GAM1模型(AUC均值0.83)優(yōu)于仔稚魚的l.GAM1模型(AUC均值0.71),兩者均可用于后續(xù)的GAM2模型分析。根據(jù)對(duì)GAM2模型的交叉驗(yàn)證線性回歸結(jié)果,e.GAM2模型的線性回歸斜率均值為1.09,截距均值為0.38,決定系數(shù)(R2)均值為0.24;l.GAM2模型的線性回歸斜率a均值為0.54,截距b均值為0.72,決定系數(shù)(R2)均值為0.15。從整體情況來看,Two-stage GAM模型對(duì)于魚卵的擬合效果和預(yù)測(cè)性能方面均優(yōu)于仔稚魚。

    表2 回歸分析統(tǒng)計(jì)結(jié)果Tab.2 Results of regression analysis

    基于GAM1得到的鳳鱭魚卵、仔稚魚的出現(xiàn)概率,設(shè)定不同的出現(xiàn)概率閾值,對(duì)于等于或高于閾值的站位點(diǎn)的相對(duì)資源密度預(yù)測(cè)值采用GAM2計(jì)算,低于閾值的站位點(diǎn)的相對(duì)資源密度定義為0作為其預(yù)測(cè)值,各站位相對(duì)資源密度的預(yù)測(cè)值與實(shí)際值進(jìn)行線性回歸擬合,以選擇適宜的出現(xiàn)概率閾值進(jìn)行相對(duì)資源密度預(yù)測(cè)。不同出現(xiàn)概率的閾值下鳳鱭魚卵、仔稚魚線性回歸的決定系數(shù)結(jié)果顯示(圖5):鳳鱭魚卵出現(xiàn)概率的閾值設(shè)定為0.4時(shí),模型的擬合性能優(yōu)于其他閾值,其線性回歸決定系數(shù)R2為0.541,線性回歸的斜率a為0.754,截距b為0.231;仔稚魚出現(xiàn)概率閾值設(shè)定為0.6時(shí)預(yù)測(cè)值與實(shí)際值的線性回歸決定系數(shù)R2最大,為0.214,其斜率a為0.368,截距b為0.95。因此本文中的鳳鱭魚卵出現(xiàn)概率的閾值為0.4,仔稚魚出現(xiàn)概率的閾值為0.6。

    圖5 出現(xiàn)概率的閾值選擇結(jié)果Fig.5 Threshold selection result of probability of occurrence

    2.4 空間分布預(yù)測(cè)

    基于Two-stage GAM對(duì)2020年5—7月甌江口鳳鱭魚卵和仔稚魚的相對(duì)資源密度的預(yù)測(cè)結(jié)果顯示(圖6):2020年6月鳳鱭魚卵和仔稚魚的適宜棲息水域明顯寬于5月和7月,鳳鱭魚卵預(yù)測(cè)的相對(duì)資源密度高值區(qū)集中在靈昆島東部附近水域以及江心嶼上游,6月份所采集到的鳳鱭魚卵基本位于江心嶼附近水域。實(shí)際調(diào)查中魚卵出現(xiàn)站位較多的月份(5月),多位于模型所預(yù)測(cè)的相對(duì)資源密度高值區(qū)的上游水域。7月的資源調(diào)查未采集到魚卵,模型所預(yù)測(cè)的7月高值區(qū)也基本位于江心嶼上游水域。模型預(yù)測(cè)的仔稚魚相對(duì)資源密度高值區(qū)則集中在七度島附近水域,6月仔稚魚的相對(duì)資源密度高值區(qū)最廣,在設(shè)立的11個(gè)調(diào)查站位中均采集到鳳鱭仔稚魚;除未采集到鳳鱭仔稚魚的5月外,6月和7月的仔稚魚相對(duì)資源密度較高站位基本位于預(yù)測(cè)的高值區(qū)內(nèi)。

    圖6 2020年5—7月鳳鱭魚卵、仔稚魚相對(duì)資源密度的預(yù)測(cè)值與實(shí)際值疊加圖Fig.6 Overlapping maps of prediction and observations of relative resource density of eggs and larvae of Coilia mystus from May to July 2020

    3 討論

    3.1 魚卵仔稚魚的時(shí)空分布

    在長(zhǎng)江口鳳鱭早期群體的相關(guān)研究中[13,22],認(rèn)為長(zhǎng)江口鳳鱭繁殖期從4月開始至10月結(jié)束。甌江口不同月份鳳鱭魚卵和仔稚魚豐度和時(shí)空分布差別明顯(圖2),仔稚魚相較魚卵的峰值出現(xiàn)時(shí)間具有一定的滯后性。與長(zhǎng)江口鳳鱭相比,甌江口鳳鱭繁殖期較短,甌江口鳳鱭的繁殖盛期基本集中在5—7月。鳳鱭繁殖群體從4月開始進(jìn)行洄游性產(chǎn)卵繁殖,在5—7月進(jìn)入繁殖盛期,至8月大部分鳳鱭基本結(jié)束繁殖并進(jìn)入近海索餌。鳳鱭魚卵主要分布在甌江上游淡水區(qū)至七都島附近咸淡水交匯區(qū),6月魚卵相對(duì)資源量最高,集中分布在江心嶼附近。鳳鱭仔稚魚則隨著月份變化從淡水區(qū)向靈昆島附近海域聚集。本研究認(rèn)為造成這種分布的原因主要有三點(diǎn):第一,鳳鱭洄游至淡水區(qū)進(jìn)行繁殖產(chǎn)卵,魚卵適應(yīng)的鹽度生態(tài)位寬度較窄,一般魚卵孵化成仔稚魚所需時(shí)間較短,而仔稚魚相較魚卵對(duì)不同水環(huán)境的適應(yīng)力更強(qiáng),使得仔稚魚較魚卵分布水域更廣泛[23];第二,鳳鱭卵為浮性卵,魚卵、仔稚魚基本隨水流作用進(jìn)行移動(dòng),其分布主要受河口徑流量影響,同時(shí)產(chǎn)卵季節(jié)徑流量的變化對(duì)河口洄游性魚類的繁殖產(chǎn)生一定的影響。入海河口徑流量的改變會(huì)使包括水溫、水鹽、溶解氧、濁度在內(nèi)的環(huán)境因子發(fā)生改變,從而影響產(chǎn)卵群體溯河洄游的數(shù)量,使甌江水域的鳳鱭魚卵、仔稚魚的數(shù)量及分布出現(xiàn)一定程度的波動(dòng)[24],這也是造成鳳鱭魚卵和仔稚魚年度空間分布差異的主要原因之一;第三,甌江口具有豐富的營(yíng)養(yǎng)鹽和浮游生物,加上該海域獨(dú)特的地理位置形成了相對(duì)平靜的海洋環(huán)境,為魚卵、仔稚魚提供了適宜孵化、生長(zhǎng)和索餌的生境[25],使鳳鱭仔稚魚7—8月數(shù)量最高且主要聚集在甌江口靈昆島附近海域。

    3.2 產(chǎn)卵場(chǎng)適宜性的關(guān)鍵影響因子

    魚類不同發(fā)育階段對(duì)環(huán)境的適應(yīng)性不同,棲息環(huán)境的改變和空間生物動(dòng)態(tài)變化是影響其時(shí)空分布的主要因素[26]。研究結(jié)果顯示,表溫、表鹽對(duì)于鳳鱭魚卵、仔稚魚的出現(xiàn)概率或豐度變化具有顯著影響(P<0.05),不同月份下水溫、水鹽等水環(huán)境的變化使鳳鱭早期群體出現(xiàn)明顯的時(shí)空變動(dòng)。

    水溫、鹽度對(duì)魚類的存活、新陳代謝及洄游分布等有不同程度的影響[27]。研究結(jié)果顯示,表溫對(duì)于鳳鱭魚卵、仔稚魚的出現(xiàn)率和豐度均有影響。魚類的生長(zhǎng)發(fā)育具有適溫范圍,當(dāng)水溫達(dá)到一定程度時(shí)刺激魚類親體性腺發(fā)育成熟并產(chǎn)卵[28],在適溫范圍內(nèi),水溫升高可以縮減魚卵孵化所需時(shí)間[29],而水溫過高不利于卵孵化和仔稚魚的生長(zhǎng)發(fā)育,造成魚類早期群體死亡率上升,致使資源量降低。鳳鱭仔稚魚相較魚卵對(duì)表層水溫的適宜性更寬[13],魚卵和仔稚魚的最適表溫在25℃左右,表溫在23—27℃時(shí)最適宜鳳鱭進(jìn)行產(chǎn)卵繁殖,也有利于魚卵孵化和仔稚魚生長(zhǎng)發(fā)育,這與胡麗娟等[13]對(duì)長(zhǎng)江口鳳鱭仔稚魚的適宜水溫研究結(jié)果一致。2015年8水溫(平均26.49℃)處于鳳鱭魚卵的適宜性水溫范圍內(nèi),較其他年份(2018年的31.36℃、2019年的30.04℃)更適宜卵的生存和孵化,這可能是當(dāng)年8月仍采集到魚卵的一個(gè)主要原因。

    河口地區(qū)的鹽度梯度變化十分復(fù)雜,魚類早期發(fā)育階段對(duì)于鹽度的適應(yīng)范圍較窄[30],鹽度影響魚類的生命活動(dòng)主要體現(xiàn)在調(diào)節(jié)體內(nèi)滲透壓平衡、改變腸道酶活性等方式,鹽度過高或過低會(huì)使魚卵孵化率降低,仔稚魚畸形率升高。本研究中鹽度對(duì)于魚卵的分布具有顯著影響,鳳鱭卵作為浮性卵,可以避免在低鹽度水環(huán)境下卵沉底堆積成塊,出現(xiàn)胚胎氧氣供應(yīng)不足等問題。甌江口鳳鱭魚卵主要分布在甌江上游淡水區(qū),說明甌江口鳳鱭群體主要洄游至低鹽水域產(chǎn)卵并完成魚卵孵化過程。與長(zhǎng)江口鳳鱭的產(chǎn)卵場(chǎng)鹽度生態(tài)位寬度相比[13],甌江口海域鳳鱭產(chǎn)卵和孵化的適宜鹽度范圍更窄,造成這種差異的原因可能是由于甌江口相比長(zhǎng)江口的徑流量小,甌江口鳳鱭洄游至最適繁殖的低鹽水域距離和時(shí)間較短。與魚卵相比,鳳鱭仔稚魚對(duì)于鹽度的適應(yīng)范圍廣,魚類早期發(fā)育階段對(duì)鹽度變化適應(yīng)能力較弱,仔稚魚的體液滲透壓在低鹽環(huán)境容易得到平衡[31],利于鳳鱭早期補(bǔ)充群體的攝食生長(zhǎng)。低鹽環(huán)境適宜鳳鱭產(chǎn)卵,對(duì)于魚卵孵化和仔稚魚的生長(zhǎng)有一定的積極作用。

    對(duì)比分析甌江口與長(zhǎng)江口鳳鱭仔稚魚的資源變化可以看出,兩個(gè)不同水域鳳鱭洄游繁殖時(shí)間主要集中在5—7月,對(duì)于水溫、水鹽的適宜范圍大致相同。早期研究將我國(guó)近海洄游性鳳鱭分為珠江型、閩江型和長(zhǎng)江型3個(gè)生態(tài)型,溫州甌江口鳳鱭應(yīng)隸屬于長(zhǎng)江型鳳鱭[6]。而鳳鱭作為短距離河口洄游性魚類,基本在既定的區(qū)域進(jìn)行繁殖發(fā)育,由于地理位置的不同及所處水域環(huán)境變化,不同水域的群體在外部環(huán)境的長(zhǎng)期影響下生態(tài)位寬度發(fā)生改變[32],使得位于甌江口及近海的鳳鱭群體與其他河口群體之間形態(tài)差異逐漸增大,甌江口鳳鱭可能已發(fā)生明顯的地理分化[33]。

    在本研究中魚卵非零值豐度數(shù)據(jù)較少,與模型的標(biāo)準(zhǔn)分布類型契合度不高,采用Two-stage GAM模型對(duì)于分析魚卵時(shí)空分布與環(huán)境變量之間關(guān)系效果優(yōu)于仔稚魚模型,在解決零值數(shù)據(jù)較多引起的模型擬合能力不足等方面有較好的效果。一般認(rèn)為魚類資源出現(xiàn)概率與密度的關(guān)系是相輔相成的,高豐度區(qū)代表更高的發(fā)現(xiàn)率,魚卵和仔稚魚的高頻率出現(xiàn)和高豐度區(qū)在一定程度上可以代表產(chǎn)卵群體及鳳鱭早期群體所需的適宜環(huán)境條件?;诙?xiàng)分布的GAM1模型可以表達(dá)魚卵、仔稚魚的出現(xiàn)概率,從已知數(shù)據(jù)預(yù)測(cè)整體分布特點(diǎn),推斷鳳鱭產(chǎn)卵場(chǎng)和育幼場(chǎng)的理論分布區(qū),基于高斯分布的GAM2模型用于解釋環(huán)境變量對(duì)于豐度變化的影響,可以更好地表達(dá)魚卵、仔稚魚的豐度變化與不同環(huán)境因子間的關(guān)系,在一定程度上反映出魚卵、仔稚魚生長(zhǎng)發(fā)育的生境適宜性需求。模型的整體解釋率不高,一方面可能是選取的自變量覆蓋范圍不足,另一方面缺乏對(duì)生物分布機(jī)制性的解釋,在以后的研究中應(yīng)結(jié)合生理學(xué)和地理水文動(dòng)力等因素做進(jìn)一步分析,了解魚卵、仔稚魚的存活條件,以期準(zhǔn)確解析鳳鱭資源的補(bǔ)充過程和種群的動(dòng)態(tài)變化。

    猜你喜歡
    稚魚甌江魚卵
    高速水流碰撞沖擊對(duì)漂流性魚卵孵化的影響研究
    A new species of the gudgeon genus Microphysogobio Mori, 1934 (Cypriniformes: Cyprinidae) from Zhejiang Province, China
    基于廣義加性模型的天津近海魚卵、仔稚魚群落結(jié)構(gòu)研究
    云南盤仔稚魚口吸盤發(fā)育*
    小金魚孵化記
    溫州甌江北口大橋綠色公路創(chuàng)建實(shí)踐
    上海公路(2019年3期)2019-11-25 07:39:26
    小金魚孵化記
    《甌江叢韻》在浙首發(fā)
    長(zhǎng)江口及其鄰近水域仔稚魚種類組成及分布特征
    云紋石斑魚仔稚魚的攝食習(xí)性與生長(zhǎng)特性
    精品一区二区三区av网在线观看| 久久精品人妻少妇| 91麻豆精品激情在线观看国产| 一级毛片久久久久久久久女| aaaaa片日本免费| 欧美另类亚洲清纯唯美| 女人被狂操c到高潮| 成年免费大片在线观看| 真人一进一出gif抽搐免费| x7x7x7水蜜桃| 午夜免费激情av| 欧美成人性av电影在线观看| 99国产精品一区二区蜜桃av| 国产精品三级大全| 床上黄色一级片| 他把我摸到了高潮在线观看| 亚洲美女视频黄频| 草草在线视频免费看| 亚洲真实伦在线观看| 久久热精品热| 日本黄大片高清| 亚洲成人免费电影在线观看| 悠悠久久av| 亚洲欧美日韩无卡精品| 婷婷六月久久综合丁香| АⅤ资源中文在线天堂| 一夜夜www| 国产精品久久久久久av不卡| 最近中文字幕高清免费大全6 | 黄色女人牲交| 91久久精品国产一区二区三区| 欧美黑人欧美精品刺激| 日日啪夜夜撸| 国产精品久久久久久亚洲av鲁大| 永久网站在线| av在线蜜桃| 99热精品在线国产| 天堂动漫精品| 一区二区三区四区激情视频 | 国产乱人伦免费视频| 亚洲av二区三区四区| 波野结衣二区三区在线| 欧美日韩黄片免| 国产探花极品一区二区| 欧美3d第一页| 真实男女啪啪啪动态图| 少妇的逼水好多| 欧美人与善性xxx| 色综合色国产| 中国美白少妇内射xxxbb| 赤兔流量卡办理| 小蜜桃在线观看免费完整版高清| 午夜视频国产福利| 亚洲三级黄色毛片| 麻豆一二三区av精品| 亚洲自拍偷在线| 人人妻,人人澡人人爽秒播| 午夜福利在线在线| 大型黄色视频在线免费观看| 黄色日韩在线| 国产精品乱码一区二三区的特点| 国产免费男女视频| 国产精品一区二区三区四区免费观看 | 观看美女的网站| 黄色女人牲交| 日韩欧美精品v在线| 久久人人精品亚洲av| 国产午夜福利久久久久久| 在线a可以看的网站| 久久精品91蜜桃| 在线天堂最新版资源| 国产免费一级a男人的天堂| 欧美+亚洲+日韩+国产| 自拍偷自拍亚洲精品老妇| 日韩欧美 国产精品| 欧美三级亚洲精品| 国产一区二区三区视频了| 国产精品综合久久久久久久免费| 欧美日本亚洲视频在线播放| av天堂在线播放| 男人舔女人下体高潮全视频| 亚洲av成人av| 亚洲七黄色美女视频| 一本一本综合久久| 免费看日本二区| 国产av麻豆久久久久久久| 色综合色国产| 国产乱人视频| 国产蜜桃级精品一区二区三区| 日韩强制内射视频| 男插女下体视频免费在线播放| 人妻久久中文字幕网| 观看免费一级毛片| 亚洲中文日韩欧美视频| 国产在线男女| 免费av观看视频| 成人av一区二区三区在线看| 精品久久久久久久久av| 麻豆成人午夜福利视频| 赤兔流量卡办理| 男女啪啪激烈高潮av片| 国产精品三级大全| aaaaa片日本免费| 国产亚洲精品av在线| 嫩草影视91久久| 亚洲熟妇中文字幕五十中出| 国产中年淑女户外野战色| 午夜免费男女啪啪视频观看 | 国产精品av视频在线免费观看| 欧美成人性av电影在线观看| 亚洲天堂国产精品一区在线| 精品一区二区三区人妻视频| 国产精品爽爽va在线观看网站| 在线免费十八禁| 欧美潮喷喷水| 午夜影院日韩av| 日本a在线网址| 亚洲欧美精品综合久久99| 91狼人影院| 久久国产乱子免费精品| 春色校园在线视频观看| 精品久久久久久久人妻蜜臀av| 日韩精品青青久久久久久| 淫秽高清视频在线观看| 成人二区视频| 久久香蕉精品热| 色吧在线观看| 桃红色精品国产亚洲av| 国产亚洲精品综合一区在线观看| 少妇被粗大猛烈的视频| 高清毛片免费观看视频网站| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 丝袜美腿在线中文| 国产综合懂色| 人人妻,人人澡人人爽秒播| 制服丝袜大香蕉在线| 亚洲成人久久爱视频| 小蜜桃在线观看免费完整版高清| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 欧美一区二区亚洲| 久久久久久久久久久丰满 | 精品免费久久久久久久清纯| 亚洲国产精品久久男人天堂| 日韩一本色道免费dvd| 最近视频中文字幕2019在线8| 成熟少妇高潮喷水视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品粉嫩美女一区| 久久精品国产亚洲网站| 亚洲经典国产精华液单| 久久久久久国产a免费观看| 国产亚洲av嫩草精品影院| 淫妇啪啪啪对白视频| 国产乱人视频| 欧美日韩乱码在线| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 色综合站精品国产| 色av中文字幕| 日本-黄色视频高清免费观看| 天堂动漫精品| 婷婷丁香在线五月| 中文字幕av成人在线电影| 极品教师在线免费播放| 国产av一区在线观看免费| 午夜福利高清视频| 久久久久久久精品吃奶| 国产激情偷乱视频一区二区| 他把我摸到了高潮在线观看| 性欧美人与动物交配| 女同久久另类99精品国产91| 午夜亚洲福利在线播放| 99精品在免费线老司机午夜| 一a级毛片在线观看| 毛片女人毛片| 久久6这里有精品| 亚洲久久久久久中文字幕| 免费在线观看影片大全网站| 午夜亚洲福利在线播放| 国内精品久久久久精免费| 欧美日本亚洲视频在线播放| 1024手机看黄色片| 白带黄色成豆腐渣| 国产精品爽爽va在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 日本五十路高清| 免费人成视频x8x8入口观看| 少妇的逼水好多| 一区福利在线观看| 99热6这里只有精品| 亚洲av五月六月丁香网| 91久久精品国产一区二区三区| 九九爱精品视频在线观看| 亚洲国产高清在线一区二区三| 中文字幕熟女人妻在线| 成人综合一区亚洲| 99热网站在线观看| 色视频www国产| 日韩中文字幕欧美一区二区| 天美传媒精品一区二区| 日日夜夜操网爽| 在线观看一区二区三区| 热99re8久久精品国产| 成人国产综合亚洲| 欧美性感艳星| 可以在线观看毛片的网站| 日韩中字成人| 1024手机看黄色片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲无线在线观看| 淫妇啪啪啪对白视频| 久久国产乱子免费精品| 99热6这里只有精品| 变态另类成人亚洲欧美熟女| 美女xxoo啪啪120秒动态图| 极品教师在线视频| 91久久精品电影网| 最后的刺客免费高清国语| 18+在线观看网站| 男女做爰动态图高潮gif福利片| 国产麻豆成人av免费视频| 欧美不卡视频在线免费观看| 亚洲性夜色夜夜综合| 乱码一卡2卡4卡精品| 男人舔奶头视频| 国产精品三级大全| 麻豆国产av国片精品| 精品日产1卡2卡| 欧美+日韩+精品| 亚洲av.av天堂| 亚洲成a人片在线一区二区| 看十八女毛片水多多多| 国产精品女同一区二区软件 | 综合色av麻豆| 高清在线国产一区| 91av网一区二区| 1024手机看黄色片| 中文亚洲av片在线观看爽| 乱人视频在线观看| 最近最新中文字幕大全电影3| 大型黄色视频在线免费观看| 日本免费a在线| 性欧美人与动物交配| av福利片在线观看| 亚洲18禁久久av| 色视频www国产| 两个人的视频大全免费| 禁无遮挡网站| 日本成人三级电影网站| 免费搜索国产男女视频| 一个人看视频在线观看www免费| 少妇的逼好多水| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av在线| 国产亚洲精品久久久com| 国产伦一二天堂av在线观看| 久久久久久大精品| 免费人成在线观看视频色| 最新中文字幕久久久久| 亚洲一区高清亚洲精品| 日韩一本色道免费dvd| 亚洲va在线va天堂va国产| 国产精品亚洲美女久久久| 哪里可以看免费的av片| 亚洲图色成人| 精品久久久久久久久亚洲 | 3wmmmm亚洲av在线观看| 日本a在线网址| 免费搜索国产男女视频| 男人舔女人下体高潮全视频| 久久久久久国产a免费观看| 国产美女午夜福利| 亚洲性久久影院| 欧美日韩精品成人综合77777| 欧美成人一区二区免费高清观看| 亚洲在线观看片| 男女啪啪激烈高潮av片| 国产精品自产拍在线观看55亚洲| 免费看日本二区| 欧美日韩综合久久久久久 | 欧美xxxx黑人xx丫x性爽| 国产精品久久电影中文字幕| 色综合婷婷激情| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 亚洲第一电影网av| 日韩中文字幕欧美一区二区| 久久精品综合一区二区三区| 成人二区视频| 小说图片视频综合网站| 精品久久国产蜜桃| 日本黄色片子视频| .国产精品久久| 久久精品综合一区二区三区| 成人鲁丝片一二三区免费| 国产白丝娇喘喷水9色精品| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美精品v在线| 一区二区三区免费毛片| 国产精品国产高清国产av| 欧美一区二区亚洲| 国产午夜精品论理片| 亚洲内射少妇av| 一本精品99久久精品77| 色吧在线观看| 国产高清三级在线| 亚洲自偷自拍三级| 久久人妻av系列| 日韩精品有码人妻一区| 欧美成人免费av一区二区三区| 日韩欧美国产一区二区入口| 国产精品嫩草影院av在线观看 | 一夜夜www| 中文资源天堂在线| 嫩草影院精品99| 99精品在免费线老司机午夜| 国产黄a三级三级三级人| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 成年女人毛片免费观看观看9| 成人国产一区最新在线观看| 在线免费观看不下载黄p国产 | 99久国产av精品| av天堂在线播放| 热99re8久久精品国产| 99久久久亚洲精品蜜臀av| 国产精华一区二区三区| 国产一级毛片七仙女欲春2| 一本久久中文字幕| 麻豆成人av在线观看| 国产精品久久久久久精品电影| 精品99又大又爽又粗少妇毛片 | 午夜久久久久精精品| 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区免费观看| 国产精品永久免费网站| 高清日韩中文字幕在线| 国产精品一区二区免费欧美| 国产精品久久久久久精品电影| 亚洲美女搞黄在线观看 | 91久久精品国产一区二区成人| 两性午夜刺激爽爽歪歪视频在线观看| 俄罗斯特黄特色一大片| 国产亚洲精品综合一区在线观看| 乱人视频在线观看| 搞女人的毛片| 欧美激情在线99| 国内精品一区二区在线观看| 欧美zozozo另类| 看免费成人av毛片| 亚洲精品影视一区二区三区av| 免费无遮挡裸体视频| 亚洲av免费在线观看| 久久精品国产99精品国产亚洲性色| 亚洲精品亚洲一区二区| 日韩亚洲欧美综合| 身体一侧抽搐| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女| 男人的好看免费观看在线视频| 亚洲成人免费电影在线观看| 国产亚洲av嫩草精品影院| 亚洲国产高清在线一区二区三| 男插女下体视频免费在线播放| 婷婷亚洲欧美| 中国美女看黄片| 女同久久另类99精品国产91| 嫁个100分男人电影在线观看| 又紧又爽又黄一区二区| 亚洲最大成人手机在线| 久久热精品热| 色视频www国产| 精品人妻视频免费看| 老司机午夜福利在线观看视频| 一级a爱片免费观看的视频| 国产单亲对白刺激| 色噜噜av男人的天堂激情| 婷婷丁香在线五月| 精品一区二区免费观看| 久久久精品大字幕| 欧美一区二区国产精品久久精品| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 老司机深夜福利视频在线观看| 日本a在线网址| 少妇高潮的动态图| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 亚洲国产精品久久男人天堂| 国产色爽女视频免费观看| 三级国产精品欧美在线观看| 亚洲狠狠婷婷综合久久图片| 久久国产精品人妻蜜桃| www日本黄色视频网| 51国产日韩欧美| 一级黄色大片毛片| 国内毛片毛片毛片毛片毛片| 又黄又爽又刺激的免费视频.| 欧美bdsm另类| 在线免费观看不下载黄p国产 | 我要看日韩黄色一级片| 淫妇啪啪啪对白视频| 久久这里只有精品中国| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 嫩草影院新地址| 久久久久久久久中文| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 五月伊人婷婷丁香| 亚洲成人久久爱视频| 亚洲国产日韩欧美精品在线观看| 日日撸夜夜添| 久久久精品大字幕| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av不卡在线观看| 国产高清激情床上av| 色尼玛亚洲综合影院| 99热只有精品国产| 可以在线观看毛片的网站| 91久久精品国产一区二区三区| 国产伦一二天堂av在线观看| 欧美成人a在线观看| 女同久久另类99精品国产91| 人妻少妇偷人精品九色| 国产真实伦视频高清在线观看 | 搡女人真爽免费视频火全软件 | 午夜视频国产福利| 一区福利在线观看| 日韩 亚洲 欧美在线| 一级黄色大片毛片| 欧美激情在线99| 国产成人一区二区在线| 91久久精品国产一区二区三区| 联通29元200g的流量卡| 精品午夜福利在线看| 一本精品99久久精品77| 在线观看一区二区三区| 又爽又黄无遮挡网站| 亚洲精品在线观看二区| 国产成人aa在线观看| 亚洲第一区二区三区不卡| 性插视频无遮挡在线免费观看| 国产亚洲精品av在线| 亚洲国产欧美人成| 国内精品久久久久精免费| 国产91精品成人一区二区三区| 99热这里只有精品一区| 日本色播在线视频| aaaaa片日本免费| 国产乱人视频| 1024手机看黄色片| 国产午夜福利久久久久久| 久久天躁狠狠躁夜夜2o2o| 一个人看视频在线观看www免费| 久久久久久久精品吃奶| 国产一区二区三区视频了| 国产中年淑女户外野战色| 精华霜和精华液先用哪个| 国产精品综合久久久久久久免费| x7x7x7水蜜桃| 国产高清激情床上av| 亚洲欧美日韩高清在线视频| 国内毛片毛片毛片毛片毛片| 小说图片视频综合网站| 狂野欧美白嫩少妇大欣赏| 2021天堂中文幕一二区在线观| 精品久久久久久成人av| 国产av不卡久久| 亚洲精华国产精华液的使用体验 | 亚洲无线在线观看| 成人国产综合亚洲| 久9热在线精品视频| av中文乱码字幕在线| 内地一区二区视频在线| 国产男人的电影天堂91| 欧美激情国产日韩精品一区| 男人的好看免费观看在线视频| 日韩一本色道免费dvd| 女人被狂操c到高潮| 99热网站在线观看| 日韩欧美在线乱码| 欧美在线一区亚洲| 亚洲一区二区三区色噜噜| 久久99热6这里只有精品| 国产精品av视频在线免费观看| 免费在线观看日本一区| 69人妻影院| 国产视频一区二区在线看| 18禁黄网站禁片午夜丰满| 成人三级黄色视频| 淫妇啪啪啪对白视频| 国产精品女同一区二区软件 | 国产成人av教育| 一a级毛片在线观看| 国内精品久久久久精免费| 国产黄a三级三级三级人| 九色国产91popny在线| 久久精品影院6| 国产免费男女视频| 国内精品久久久久精免费| 亚洲精品国产成人久久av| 亚洲色图av天堂| 国产人妻一区二区三区在| eeuss影院久久| 人妻丰满熟妇av一区二区三区| 听说在线观看完整版免费高清| .国产精品久久| 亚洲人与动物交配视频| 亚洲美女视频黄频| 少妇的逼水好多| 欧美黑人巨大hd| 亚洲成人久久性| 婷婷丁香在线五月| 国产av不卡久久| 看免费成人av毛片| 成人av在线播放网站| 搡老岳熟女国产| 国内精品久久久久精免费| 国产黄色小视频在线观看| 国产午夜福利久久久久久| 国产69精品久久久久777片| 国产精品,欧美在线| 国产精品人妻久久久久久| 少妇猛男粗大的猛烈进出视频 | 成人综合一区亚洲| 少妇丰满av| 欧美日韩亚洲国产一区二区在线观看| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 可以在线观看的亚洲视频| 久久精品人妻少妇| 欧美3d第一页| 国产蜜桃级精品一区二区三区| 欧美国产日韩亚洲一区| 久久精品国产自在天天线| 欧美极品一区二区三区四区| 美女黄网站色视频| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 亚洲av二区三区四区| 国产欧美日韩一区二区精品| 中文字幕av成人在线电影| x7x7x7水蜜桃| 成人三级黄色视频| 国内久久婷婷六月综合欲色啪| 免费不卡的大黄色大毛片视频在线观看 | 国产主播在线观看一区二区| 别揉我奶头 嗯啊视频| 午夜福利欧美成人| 黄色丝袜av网址大全| 美女黄网站色视频| 精品久久久久久久久久久久久| www.色视频.com| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 亚洲在线自拍视频| 国产精品三级大全| 久久久久久久久久成人| 99国产精品一区二区蜜桃av| 网址你懂的国产日韩在线| 久久6这里有精品| 免费av不卡在线播放| 真人一进一出gif抽搐免费| 亚洲成人久久性| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 一进一出抽搐动态| 神马国产精品三级电影在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 中文字幕免费在线视频6| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 男女做爰动态图高潮gif福利片| 亚洲人成网站高清观看| 国产一区二区亚洲精品在线观看| 亚洲午夜理论影院| 日韩欧美免费精品| 少妇熟女aⅴ在线视频| 亚洲精华国产精华精| 亚洲成人免费电影在线观看| 好男人在线观看高清免费视频| 1000部很黄的大片| 18禁黄网站禁片午夜丰满| 人人妻人人澡欧美一区二区| 干丝袜人妻中文字幕| avwww免费| av在线老鸭窝| 亚洲熟妇中文字幕五十中出| 午夜精品一区二区三区免费看| 久久精品综合一区二区三区| 中国美白少妇内射xxxbb| 深爱激情五月婷婷| 亚洲欧美日韩东京热| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久久久毛片| av女优亚洲男人天堂| 国产精品无大码| 国产在视频线在精品| 国产黄a三级三级三级人| 欧美人与善性xxx| a在线观看视频网站| 97碰自拍视频| 99热这里只有精品一区| 免费高清视频大片| 欧美人与善性xxx|