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

    淮河上游南灣湖麥穗魚年齡與生長的研究

    2017-08-16 06:06:56李紅敬王雅平冷秋麗李曉杰李曉鳳于同雷
    水生生物學(xué)報 2017年4期
    關(guān)鍵詞:雌雄體長雌性

    李紅敬 王雅平 冷秋麗 李曉杰 李曉鳳 于同雷 黃 斌

    (信陽師范學(xué)院生命科學(xué)學(xué)院, 信陽 464000)

    淮河上游南灣湖麥穗魚年齡與生長的研究

    李紅敬 王雅平 冷秋麗 李曉杰 李曉鳳 于同雷 黃 斌

    (信陽師范學(xué)院生命科學(xué)學(xué)院, 信陽 464000)

    2014年的3月、5月、7月和12月在淮河上游南灣湖采集麥穗魚(Pseudorasbora parva)樣本532尾, 對麥穗魚的年齡組成與生長進(jìn)行分析。結(jié)果表明樣本的體長分布范圍為35.82—88.28 mm, 平均體長為(61.61±11.8) mm,體重的分布范圍為3.07—59.17 g, 平均體重為(19.23±10.73) g。雄性個體比雌性個體大, 雌雄性比為0.64鯰1。群體的年齡組成為1—3齡, 其中3齡樣本數(shù)量占優(yōu)勢為57.38%。體長與體重的關(guān)系是雌性W=9.602E–5L2.928(R2=0.883); 雄性W=4.487E–5L3.116(R2=0.889), 雌雄樣本間存在顯著性差異(F=5.241, P<0.05)。麥穗魚的鱗徑與體長之間呈線性關(guān)系, 并且雌雄樣本的鱗徑與體長之間的關(guān)系差異性顯著(F=78.405, P<0.05)。生長參數(shù)分別是雌性: L∞=107.005, K=0.246, t0= –0.76; 雄性: L∞=145.254, K=0.181, t0= –0.66。生長拐點是雌性3.607齡對應(yīng)的體長和體重分別為70.46 mm和24.72 g, 雄性5.619齡對應(yīng)的體長和體重分別為98.64 mm和73.53 g。研究結(jié)果表明雌性為勻速生長, 雄性為異速生長; 雄性麥穗魚比雌性麥穗魚的生長速度快。

    南灣湖; 麥穗魚; 年齡; 生長

    麥穗魚(Pseudorasbora parva)隸屬于鯉形目(Cypriniformes), 鯉科(Cyprinidae), 逗亞科(Gobioninae), 麥穗魚屬(Pseudorasbora)。棲息于水淺、流速緩慢水體的中、底層, 尤其是水生植物叢生、富含有機(jī)質(zhì)的池塘、湖泊沿岸帶等水域[1,2], 對不同棲息環(huán)境有著較強(qiáng)的適應(yīng)性; 食雜性[3], 以浮游動物、底棲動物為主要食物[4]。世代周期短、性成熟早、繁殖力強(qiáng)[5]。屬于東亞地區(qū)的土著物種, 廣泛分布于中國大陸[6]、日本[7]、朝鮮[8]、中國臺灣[9]等國家和地區(qū)。20世紀(jì)50年代從黑龍江流域入侵至蘇聯(lián)[10], 此后快速的入侵至歐洲的其他國家, 目前其分布已遍及全球[11]。

    在東亞地區(qū)分布的麥穗魚在生態(tài)系統(tǒng)中占有重要地位, 不僅能為肉食性魚類提供食物資源, 而且能夠維持浮游動植物種群的穩(wěn)定[4,12], 而在其他入侵地, 麥穗魚的入侵給當(dāng)?shù)佤~類生態(tài)帶來了巨大影響[13], 甚至威脅到一些土著種的生存, 如云南的滇池、洱海和撫仙湖就有多種土著魚類因麥穗魚等入侵種而日趨減少或瀕臨滅絕[14]。在英國湖區(qū)麥穗魚在入侵后4年繁殖成為當(dāng)?shù)財?shù)量最多的魚類,而威脅到其他本土魚類的生存[15]。對入侵麥穗魚的管理和控制是目前亟待解決的問題。

    魚類年齡與生長的研究, 是分析和評價魚類種群數(shù)量變動的基本依據(jù)之一。目前在亞洲地區(qū)對于麥穗魚生物學(xué)方面的研究主要集中于體長與體重關(guān)系[16]以及生長參數(shù)的比較[17]等方面, 對于麥穗魚生長規(guī)律等方面的研究還較為缺乏。而在麥穗魚的入侵地對其生物學(xué)方面的研究, 主要是通過比較種群分布以及個體生長的變化情況來確定環(huán)境對其生長的影響[18—21]。目前對于入侵麥穗魚的治理主要是采取生物防治[22]。因此, 我們以分布于東亞大陸淮河上游區(qū)域的南灣湖及其周邊溪流中的麥穗魚為研究對象, 對麥穗魚種群的生長特性開展研究。旨在了解亞熱帶向暖溫帶過渡區(qū)麥穗魚的自然生長規(guī)律, 為入侵麥穗魚的防治和管理以及麥穗魚的應(yīng)用研究提供參考資料, 對淮河流域魚類資源的資料進(jìn)行補(bǔ)充。

    1 材料與方法

    1.1 樣本采集

    南灣湖(32°08′N, 113°58′E)位于淮河上游, 受東亞季風(fēng)氣候的影響, 具有北亞熱帶向暖溫帶過渡的季風(fēng)氣候和山地氣候特征, 與同緯度的其他地區(qū)相比, 冬夏氣溫均偏低, 年均溫15.2℃, 年均降雨量約1140 mm, 年平均蒸發(fā)量約935 m3。

    2014年3月、5月、7月、10月和12月, 在信陽市南灣湖的三個較大入水河口附近(譚家河、董家河以及浉河港), 使用成套刺網(wǎng)在采樣區(qū)內(nèi)對魚類進(jìn)行隨機(jī)采樣。每套刺網(wǎng)由8種網(wǎng)目規(guī)格各異的浮性刺網(wǎng)(網(wǎng)目大小依次為20、25、33、40、50、60、70和80 mm, 均高1.5 m、長20 m) 串聯(lián)而成, 用以捕獲較大體長范圍內(nèi)的魚類個體。日落之前, 在采樣區(qū)下兩套共16條刺網(wǎng), 刺網(wǎng)上端浮子浮力大于底端鉛錘質(zhì)量; 第2天早晨天亮后收網(wǎng), 上岸摘魚,在整個采樣期間, 這兩套刺網(wǎng)按同樣采樣程序用于各個采樣點, 以保障漁獲物數(shù)據(jù)之間的可比性。樣本采集后在新鮮狀態(tài)下運送至實驗室, 用紗布和吸水紙將樣本表面多余的水分擦拭干凈, 使每尾樣本在測量時都保持一樣的狀態(tài), 分別進(jìn)行編號, 用游標(biāo)卡尺測量體長和全長(精確至0.01 mm)。使用FA2104A型電子天平稱量體重(精確至0.01 g)。取每尾樣本的背鰭起點下方至側(cè)線上方的鱗片5—10枚, 用于年齡鑒定。

    1.2 數(shù)據(jù)處理

    體長與體重的關(guān)系采用Keys公式W=aLb對雌雄樣本分別進(jìn)行擬合。其中W表示體重(g)、L表示體長(mm)、a是常數(shù)、b是冪指數(shù)。以體長的自然對數(shù)值作為協(xié)變量進(jìn)行方差分析, 檢驗雌、雄個體間體長與體重的關(guān)系是否存在差異。體長體重方程中指數(shù)b值能反映魚類在不同發(fā)育階段和不同生境中的生長發(fā)育情況, 當(dāng)魚類體長體重方程的b值等于3或接近3時, 體長與體重的生長為等速生長。

    體長與鱗徑的關(guān)系采用不同的鱗徑與體長的關(guān)系方程, 分別為: Einar Lea、Rosa Lee和蒙納斯蒂爾斯基方程Ln=Lrn/R、L=a+bR、以及L=aRb對麥穗魚體長與鱗徑的關(guān)系進(jìn)行擬合, 選擇判定系數(shù)(R2)最大者為其最佳回歸方程; 采用Rosa Lee方程對雌雄樣本的體長進(jìn)行退算。協(xié)方差分析雌雄樣本的鱗徑與體長的關(guān)系。

    數(shù)據(jù)的統(tǒng)計整理使用Microsoft excel2007, 實驗結(jié)果的分析和圖表的制作使用Statistica 8.0。

    2 結(jié)果

    2.1 樣本組成與性比

    使用鱗片能夠成功鑒定出427尾麥穗魚樣本的年齡, 其余105尾樣本由于鱗片破損或者其他原因未能成功鑒定出年齡。對427尾麥穗魚樣本的年齡進(jìn)行統(tǒng)計。其中1齡樣本較少, 占總樣本的6.32%, 2齡樣本占36.3%, 3齡樣本較多占57.38%。

    通過對427尾麥穗魚體長、體重的測量和統(tǒng)計(圖 1、圖 2), 其中雌性167尾, 雄性260尾。性比為0.64鯰1, 用χ2分析比較性比與1鯰1之間存在顯著性差異(P<0.05)。 427尾麥穗魚總樣本的體長范圍為35.82—88.28 mm, 平均體長為(61.61±11.8) mm。雌性樣本的體長為35.82—87.76 mm, 平均體長為(55.99±10.52) mm。雄性樣本的體長為36.83—88.28 mm, 平均體長為(65.22±11.17) mm。雌性樣本和雄性樣本的體長分布存在顯著性差異(Kolmogorov-Smirnov檢驗D=4.635, P<0.05)。

    總樣本的體重范圍為3.07—59.17 g, 平均體重為(19.23±10.73) g。雌性樣本的體重范圍為3.07—49.54 g, 平均體重為(14.08±8.06) g。雄性樣本的體重范圍為3.08—59.17 g, 平均體重為(22.54±10.94) g。雌性樣本和雄性樣本的體重分布存在顯著性差異(Kolmogorov-smirnov檢驗D=4.597, P<0.05)。

    2.2 體長與體重的關(guān)系

    圖 1 麥穗魚總樣本群體及雌、雄樣本體長組成Fig.1 The distribution patter of body length of the total samples, male and female Pseudorasbora parva

    對麥穗魚的體長和體重的關(guān)系進(jìn)行擬合(圖 3),其相關(guān)式為: 雌性W=9.602E–5L2.928R2=0.883; 雄性W=4.487E–5L3.116R2=0.889; 將雌雄樣本的體長、體重進(jìn)行對數(shù)處理, 使其成線性關(guān)系, 用協(xié)方差分析(ANCOVA)檢驗雌雄之間體長、體重關(guān)系存在顯著性差異(F=5.241, P<0.05)。分別將雌雄樣本以及總樣本的長重關(guān)系指數(shù)與勻速生長魚類的b=3進(jìn)行比較, 雌性(b=2.928)與3之間無顯著性差異(t=–1.787, P>0.05), 雄性(b=3.116)與3之間存在顯著性差異(t=9.278, P<0.05), 即雌性為勻速生長, 雄性為異速生長??梢杂肰on Bertalanffy生長方程(VBGF)來描述其生長特性。

    圖 2 麥穗魚總樣本群體及雌、雄樣本體重組成Fig.2 The distribution pattern of body weight of the total samples, male and female Pseudorasbora parva

    圖 3 麥穗魚體長與體重的關(guān)系Fig.3 The relationship between body length and body weight of Pseudorasbora parva

    2.3 體長與鱗徑的關(guān)系

    采用Einar Lea、Rosa Lee和蒙納斯蒂爾斯基這3個方程對鱗徑與體長的關(guān)系進(jìn)行擬合, 并根據(jù)擬合出的方程退算出雌雄樣本的體長(表 1)。通過對3種方法的退算體長與實測體長的比較, 結(jié)果顯示Rosa Lee公式退算出的體長值與實測體長的差異性較小, 因此確定麥穗魚樣本的鱗徑與體長的關(guān)系式為, 雌性: L=16.399+0.053R; 雄性: L=23.501+ 0.0414R。

    協(xié)方差分析顯示雌雄樣本的鱗徑與體長的關(guān)系存在顯著性差異(ANCOVA分析F=78.405, P<0.05)。

    2.4 體長的退算

    根據(jù)鱗徑與體長的關(guān)系, 對麥穗魚雌雄樣本各齡的體長進(jìn)行退算(表 2)。對雌雄實測體長與退算體長分別進(jìn)行χ2檢驗結(jié)果顯示, 雌性P>0.05 (χ2=0.079, df=2, P=0.961), 雄性P>0.05 (χ2=0.108, df=2, P=0.948), 即雌雄樣本的實測體長與退算體長之間不存在顯著性差異。表明選用Rosa Lee公式作為鱗徑與體長的關(guān)系式是可靠的。

    2.5 生長方程

    根據(jù)平均退算體長以及體長與體重的關(guān)系, 計算出各齡的退算體重。根據(jù)ford方程求得雌雄群體的生長參數(shù), 根據(jù)Beverton法從ln(L∞–Lt)=lnL∞+ Kt0–Kt求得理論體長或體重等于0時的年齡t0。各生長參數(shù)如下, 雌性: L∞=107.005, W∞=84.034, K=0.246, t0=–0.76; 雄性: L∞=145.254, W∞=244.989, K=0.181, t0=–0.66。將計算出的生長參數(shù)代入Von Bertalanffy生長方程, 得出麥穗魚雌雄群體體長、體重的生長方程分別為, 雌性: Lt=107.005[1–e–0.246(t+0.76)], Wt=84.034[1–e–0.246(t+0.76)]2.928; 雄性: Lt=145.254 [1–e–0.181(t+0.66)] , Wt=244.989[1–e–0.181(t+0.66)]3.116

    表 1 不同方法退算麥穗魚體長的比較Tab.1 The comparison of backed-calculated length of different methods to Pseudorasbora parva

    F檢驗結(jié)果顯示雌雄的體長生長方程回歸極顯著(P<0.01)。將理論體長與退算體長進(jìn)行χ2檢驗,結(jié)果顯示雌性P>0.05, 雄性P>0.05, 差異性不顯著,說明體長的生長曲線擬合程度可靠。從生長曲線(圖 4)可以看出, 雌雄群體體長的生長都不具拐點,開始生長很快, 但是隨著年齡的增長, 其生長逐漸接近于漸進(jìn)值, 而體重的生長則近似于“S”形生長。

    表 2 麥穗魚各齡退算體長和平均實測體長Tab.2 Backed-calculated Length (BL) and mean observed BL of Pseudorasbora parva

    分別對體長和體重生長方程求一階導(dǎo)數(shù), 即得出其生長速度的方程:

    雌性: dL/dt=26.323e–0.246(t+0.76); dW/dt=60.529e–0.246(t+0.76)[1–e–0.246(t+0.76)]1.928

    雄性: dL/dt=26.291e–0.181(t+0.66); dW/dt=138.173e–0.181(t+0.66)[1–e–0.181(t+0.66)]2.116

    對生長方程求二階導(dǎo)數(shù)可得其生長加速度的方程:

    雌性: d2L/dt2= –6.476e–0.246(t+0.76)

    d2W/dt2=14.89e–0.246(t+0.76)[1–e–0.246(t+0.76)]0.928[2.928e–0.246(t+0.76)–1]

    雄性: d2L/dt2= –4.759e–0.181(t+0.66)

    d2W/dt2=25.009e–0.181(t+0.66)[1–e–0.181(t+0.66)]1.116[3.116e–0.181(t+0.66)–1]

    體長和體重的生長速度和加速度曲線如圖 5—8所示。計算出麥穗魚生長的拐點年齡分別為雌性3.607齡對應(yīng)的體長和體重分別為70.46 mm和24.72 g,雄性5.619齡對應(yīng)的體長和體重分別為98.64 mm和73.53 g。

    圖 4 麥穗魚體長、體重生長曲線Fig.4 The growth curve of Pseudorasbora parva

    3 討論

    圖 5 雄性體長生長速度與加速度曲線Fig.5 Growth in body length as dL/dt and d2L/dt2of male Pseudorasbora parva

    圖 6 雄性體重生長速度與加速度方程Fig.6 Growth in body length as dW/dt and d2W/dt2cruves of male Pseudorasbora parva

    目前, 對麥穗魚的年齡鑒定, 通常使用長度頻度法[17]和鱗片鑒定法[16,19]。我們選取鱗片作為麥穗魚的年齡鑒定材料, 主要是因為麥穗魚的鱗片易處理輪紋清晰。并且多數(shù)小型鯉科魚類都用鱗片鑒定其種群的年齡[23,24]。

    圖 7 雌性體長生長速度與加速度Fig.7 Growth in body weight as dW/dt and d2W/dt2cruves of female Pseudorasbora parva

    圖 8 雌性體重生長速度與加速度Fig.8 Growth in body weight as dW/dt and d2W/dt2of female Pseudorasbora parva

    鑒定結(jié)果顯示南灣湖麥穗魚的最大年齡為3齡,最大體長是88.28 mm, 小于保安湖麥穗魚的最大體長100.5 mm[16]。大于日本九州島麥穗魚的最大體長60 mm[17]。這種差異出現(xiàn)的原因主要是環(huán)境不同導(dǎo)致的。九州島麥穗魚樣本采集自有明海下游的一條人工灌溉溝渠, 溝渠中水生生物數(shù)量較少,不能為魚類的生長提供足夠的食物, 從而影響其生長。而保安湖緯度較低, 水體溫度較高更加適合魚類的生長。而在入侵地布拉迪斯拉發(fā)的麥穗魚最大年齡為5齡[19], 普雷斯帕湖為4齡[18]。在英國2個不同池塘中, 麥穗魚的最大體長分別是76和81 mm[20]。在入侵地麥穗魚的壽命有所延長, 體長相對較小,說明麥穗魚對環(huán)境的適應(yīng)能力強(qiáng), 入侵后并未對種群的生長造成影響。

    研究結(jié)果顯示麥穗魚的雌雄群體的生長方式存在差異, 雌性群體為勻速生長(b=2.928);雄性群體為異速生長(b=3.116)。張?zhí)昧值萚16]研究保安湖麥穗魚的體長與體重的關(guān)系, 雄性麥穗魚樣本在4月、5月和10月屬于勻速生長方式, 其他月份則表現(xiàn)為異速生長的特征; 雌性麥穗魚樣本在4月、5月表現(xiàn)為勻速生長, 其他月份表現(xiàn)為異速生長的特征。而其他地區(qū)的麥穗魚種群的生長表現(xiàn)為異速生長方式[18]。這主要是因為魚類體長與體重的關(guān)系會受到性別、性腺發(fā)育以及季節(jié)的影響[25]。

    我們確定麥穗魚體長與鱗徑的關(guān)系是通過比較三種不同的鱗徑與體長的關(guān)系而篩選出來的最適合的方法。Einar Lea公式認(rèn)為鱗片與體長之間呈正比, Rosa Lee公式認(rèn)為鱗片與體長呈線性關(guān)系,并且多數(shù)魚類的鱗徑與體長都符合這一關(guān)系[26,27]。也有魚類的鱗徑與體長呈曲線關(guān)系[28]我們通過比較不同方法對體長的退算, 從而確定麥穗魚的鱗徑與體長之間呈線性關(guān)系, 并根據(jù)這一關(guān)系對雌雄群體的體長進(jìn)行了退算, 退算結(jié)果略小于實測體長,出現(xiàn)輕微的李氏現(xiàn)象。但是退算值與實測值之間并不存在顯著性差異(P>0.05)說明用線性關(guān)系擬合麥穗魚鱗徑與體長的關(guān)系是可取的, 而其他地區(qū)麥穗魚的研究結(jié)果也證實了其鱗徑與體長之間是線性關(guān)系[16,19]。

    魚類的生長通常會有一個適合其生長參數(shù)的數(shù)學(xué)模型來分析其生長趨勢[29]。von Bertalanffy生長方程是一個主要用來描述年齡與體長之間非線性關(guān)系的魚類生長模型[30]。我們通過von Bertalanffy生長方程計算出麥穗魚的漸進(jìn)體長雄性145.254 mm大于雌性107.005 mm, 但是都要大于所采集麥穗魚樣本的最大實測體長(雄性88.28 mm,雌性87.76 mm), 說明本研究對于麥穗魚的生長參數(shù)的估算是可靠的。K值表示魚類個體的生長接近漸進(jìn)值的速度, 在本研究中雌性麥穗魚的K值(0.246)大于雄性(0.181), 即雌性個體生長到漸進(jìn)體長的速度比雄性快, 也就符合了雌性的漸進(jìn)體長小于雄性的現(xiàn)象。將本研究的麥穗魚的各項生長指標(biāo)與其他地區(qū)的麥穗魚的生長指標(biāo)進(jìn)行比較, 俄羅斯阿穆爾河麥穗魚的K值0.41相對較大, 漸進(jìn)體長是109 mm較小[31], 說明其生長到接近于漸進(jìn)體長的速度更快, 這可能與俄羅斯的氣候有關(guān), 俄羅斯地處北溫帶, 氣候較寒冷, 因此魚類的生長也受到限制。而希臘普雷斯帕湖[18]的麥穗魚屬于入侵種,但是當(dāng)?shù)佧溗媵~的K值0.25和漸進(jìn)體長118 mm與本研究的結(jié)果并沒有太大的差異性, 這也說明了麥穗魚具有較強(qiáng)的適應(yīng)能力, 同時這也可能與希臘的氣候有關(guān), 希臘屬于地中海型氣候, 全年氣溫變化不大, 更加適合魚類的生長。

    由于目前漁獲物的主要組成為2齡和3齡個體,缺少小個體、低年齡的個體, 并且漁獲物中高齡魚個體極少, 這種現(xiàn)象的出現(xiàn), 捕撈的選擇性可能是導(dǎo)致漁獲物1齡個體樣本偏少的主要原因。在采集小型魚類時, 依靠刺網(wǎng)采集樣本, 容易造成漁獲物組成的選擇性, 應(yīng)配合其他網(wǎng)具采集。本研究樣本采集主要在入湖的3個較大河口處采集, 在南灣湖的中心及其他區(qū)域未進(jìn)行采集, 也會造成樣本的選擇性。

    麥穗魚屬于小型鯉科魚類, 其分布范圍廣, 對環(huán)境適應(yīng)能力強(qiáng), 雖然麥穗魚不屬于經(jīng)濟(jì)魚類, 但在水體食物網(wǎng)中起重要作用, 對水體變化敏感, 可以作為水體生態(tài)系統(tǒng)變化的指示生物, 對生態(tài)系統(tǒng)的影響較大, 麥穗魚對維持水生生態(tài)系統(tǒng)的穩(wěn)定具有重要作用, 不僅能夠減緩水藻帶來的水華的危害,而且對浮游動物的生存動態(tài)具有調(diào)控作用。

    我們計算出麥穗魚的生長拐點雌性3.607齡, 雄性5.619齡。在拐點年齡之前, 麥穗魚個體快速生長, 拐點年齡之后生長開始減慢。而我們所采集的麥穗魚樣本年齡都處于拐點年齡之前, 出現(xiàn)這種狀況可能是過度捕撈所致。這表明南灣湖麥穗魚資源利用是不合理的, 過度捕撈使種群結(jié)構(gòu)遭到破壞,資源量日趨減少, 種群出現(xiàn)了明顯的小型化趨勢,嚴(yán)重影響了麥穗魚野生資源的可持續(xù)性。目前, 應(yīng)采取限制電、炸、毒捕撈和捕撈漁具規(guī)格, 以及在繁殖季節(jié)加強(qiáng)保護(hù)禁止捕撈等措施, 以利于麥穗魚資源的恢復(fù)和增殖。

    [1]Asaeda T, Manatunge J.Foraging of a small planktivore (Pseudorasbora parva: Cyprinidae) and its behavioral flexibility in an artificial stream [J].Hydrobiologia, 2005, 549(1): 155—166

    [2]Asaeda T, Manatunge J.Physiological responses of topmouth gudgeon, Pseudorasbora parva, to predator cues and variation of current velocity [J].Aquatic Ecology, 2007, 41(1): 111—118

    [3]Nagata T, Ha J Y, Hanazato T.The predation impact of larval Pseudorasbora parva (Cyprinidae) on zooplankton: a mesocosm experiment [J].Journal of Freshwater Ecology, 2005, 20(4): 757—763

    [4]Zhang T, Cui Y, Li Z, et al.The ecology of topmouth gudgeon, Pseudorasbora parva, in Bao’an Lake–III.Food habits [J].Acta Hydrobiologica Sinica, 1998, 22(Suppl): 157—167 [張?zhí)昧? 崔奕波, 方榕樂, 等.保安湖麥穗魚種群生物學(xué)Ⅲ食性.水生生物學(xué)報, 1998, 22(增刊): 157—167]

    [5]Britton J R, Davies G D, Brazier M, et al.A case study on the population ecology of a topmouth gudgeon (Pseudorasbora parva) population in the UK and the implications for native fish communities [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2007, 17(7): 749—759

    [6]Hu M, Wu Z Q, Liu Y.The fish fauna of mountain streams in the Guanshan National Nature Reserve, Jiangxi, China [J].Environmental Biology of Fishes, 2009, 86(1): 23—27

    [7]Onikura N, Nakajima J.Age, growth and habitat use of the topmouth gudgeon, Pseudorasbora parva in irrigation ditches on northwestern Kyushu Island, Japan [J].Journal of Applied Ichthyology, 2013, 29(1): 186—192

    [8]Jang M H, Yoon J D, Shin J H, et al.Status of freshwater fish around the Korean demilitarized zone and its implications for conservation [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2008, 18(6): 819—828

    [9]Chen I S, Weng C J, Chen Y R, et al.The checklist of inland-water and mangrove fish fauna of Kinmen Island, Fujian Province, Taiwan with comments on ecological conservation of native fishes [J].Journal of Marine Science and Technology, 2013, 21: 316—319

    [10]Britton J R, Davies G D, Brazier M, et al.A case study on the population ecology of a topmouth gudgeon (Pseudorasbora parva) population in the UK and the implications for native fish communities [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2007, 17(7): 749—759

    [11]Gozlan R E, Andreou D, Asaeda T, et al.Pan-continental invasion of Pseudorasbora parva: towards a better understanding of freshwater fish invasions [J].Fish and Fisheries, 2010, 11(4): 315—340

    [12]Chen Y, Li Z, Lin L S.Effects of a small planktivore (Pseudorasbora parva: Cyprinidae) on eutrophication of a shallow eutrophic lake in central China [J].Aquatic Ecosystem Health & Management, 2010, 13(3): 328—334

    [13]Fan L Q, Liu H P, Lin J, et al.Non-native fishes: distribution and assemblage structure in the Ihasa River basin, Tibet, China [J].Acta Hydrobiologica Sinica, 2016, 40(5): 958—967 [范麗卿, 劉海平, 林進(jìn), 等.拉薩河流域外來魚類的分布、群落結(jié)構(gòu)及其與環(huán)境的關(guān)系.水生生物學(xué)報, 2016, 40(5): 958—967]

    [14]Yan Y, Chen Y.Variations in reproductive strategies between one invasive population and two native populations of Pseudorasbora parva [J].Current Zoology, 2009,55(1): 56—60

    [15]Britton J R, Davies G D, Brazier M, et al.A case study on the population ecology of a topmouth gudgeon (Pseudorasbora parva) population in the UK and the implications for native fish communities [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2007, 17(7): 749—759

    [16]Zhang T L, Cui Y B, Fang, R L, et al.Population biology of topmouth gudgeon (Pseudorasbora parva) in Bao’an Lake-I Age and growth [J].Acta Hydrobiologica Sinica, 1998, 22(Suppl): 139—145 [張?zhí)昧? 崔奕波, 方榕樂,等.保安湖麥穗魚種群生物學(xué)Ⅰ年齡與生長.水生生物學(xué)報, 1998, 22(增刊): 139—145]

    [17]Onikura N, Nakajima J.Age, growth and habitat use of the topmouth gudgeon, Pseudorasbora parva in irrigation ditches on northwestern Kyushu Island, Japan [J].Journal of Applied Ichthyology, 2013, 29(1): 186—192

    [18]Rosecchi E, Crivelli A J, Catsadorakis G.The establishment and impact of Pseudorasbora parva, an exotic fish species introduced into Lake Mikri Prespa (north-western Greece) [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 1993, 3(3): 223—231

    [19]Záhorská E, Ková? V, Katina S.Age and growth in a newly-established invasive population of topmouth gudgeon [J].Central European Journal of Biology, 2010, 5(2): 256—261

    [20]Davies G D, Britton J R.Influences of population density, temperature and latitude on the growth of invasive topmouth gudgeon Pseudorasbora parva [J].Ecology of Freshwater Fish, 2015, 24(1): 91—98

    [21]Lemmens P, Mergeay J, Vanhove T, et al.Suppression of invasive topmouth gudgeon Pseudorasbora parva by native pike Esox lucius in ponds [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2015, 25(1): 41—48

    [22]Maceina M J, Boxrucker J, Buckmeier D L, et al.Current status and review of freshwater fish aging procedures used by state and provincial fisheries agencies with recommendations for future directions [J].Fisheries, 2007, 32(7): 329—340

    [23]Quist M C, Jackson Z J, Bower M R, et al.Precision of hard structures used to estimate age of riverine catostomids and cypri-nids in the upper Colorado River basin [J].North American Journal of Fisheries Management, 2007, 27(2): 643—649

    [24]Yan Y Z, Xu Y S, Chu L, et al.Age, growth and reproduction of Sarcocheilichthys nigripinnis from the Qingyi Stream in the Huangshan Mountains [J].Zoological Research, 2012, 33(E1-2): 25—31

    [25]Pauly D.A preliminary compilation of fish length growth parameters [J].Berichte des Instituts für Meereskunde an der Christian-Albrechts Universit?t, Kiel, 1978, 55: 1—200

    [26]Tarkan A S, Copp G H, Zi?ba G, et al.Growth and reproduction of threatened native crucian carp Carassius carassius in small ponds of Epping Forest, south-east England [J].Aquatic Conservation: Marine and Freshwater Ecosystems, 2009, 19(7): 797—805

    [27]Kumar R S, Sarkar U K, Gusain O, et al.Age, growth, population structure and reproductive potential of a vulnerable freshwater mullet, Rhinomugil corsula (Hamilton,1822) from a tropical river betwa in central India [J].Proceedings of the National Academy of Sciences, India Section B: Biological Sciences, 2014, 84(2): 275—286

    [28]Zhou C, Zhu Q, Liu H Z.Analysis of the growth function of the largemouth bronze gudgeon Coreius guichenoti in the Upper Reaches of the Yangtze River [J].Sichuan Journal of Zoology, 2010, 29(4): 510—516 [周燦, 祝茜,劉煥章.長江上游圓口銅魚生長方程的分析.四川動物, 2010, 29(4): 510—516]

    [29]Quince C, Abrams P A, Shuter B J, et al.Biphasic growth in fish I: theoretical foundations [J].Journal of Theoretical Biology, 2008, 254(2): 197—206

    [30]Cope J M, Punt A E.Admitting ageing error when fitting growth curves: an example using the von Bertalanffy growth function with random effects [J].Canadian Journal of Fisheries and Aquatic Sciences, 2007, 64(2): 205—218

    [31]Ivancheva E Y, Ivanchev VP, Sarychev V S.Distribution of stone moroko (Pseudorasbora parva) in the Upper Don basin [J].Russian Journal of Biological Invasions, 2014, 5(3): 151—155

    STUDY ON THE AGE AND GROWTH OF PSEUDORASBORA PARVA FROM NANWAN LAKE UPSTREAM THE HUAIHE RIVER

    LI Hong-Jing, WANG Ya-Ping, LENG Qiu-Li, LI Xiao-Jie, LI Xiao-Feng, YU Tong-Lei and HUANG Bin
    (College of Life Sciences, Xinyang Normal University, Xinyang 464000, China)

    Pseudorasbora parva is a native small cyprinid fish widely distributes in Asia with a considerable population.To investigate the growth, age composition and growth characteristics of Pseudorasbora parva, 532 individuals were collected from three different water inlets of Nanwan Lake in March, May, July and December.The results showed that the body length were ranged from 35.82 to 88.28 mm with the mean of (61.61±11.8) mm, and that the body weight were ranged from 3.07 to 59.17 g with the mean of (19.23±10.73) g.The age were ranged from 1 to 3.3 years and the dominant population with the mean of 57.38% were individuals with 3 years.The relationship between body length and body weight were W=9.602E–5L2.928(R2=0.883) for female and W=4.487E–5L3.116(R2=0.889) for male, which showed a significant difference between male and female (F=5.241, P<0.05).A linear relationship was found between the scale radius and body length, which was significantly affected by gender (F=78.405, P<0.05).Growth parameters were L∞=107.005, K=0.246 and t0= –0.76 for female and L∞=145.254, K=0.181 and t0= –0.66 for male, respectively.The growth inflexion of female (3.607 years) were 70.46 mm for body length and 24.72 g for body weight, and the growth inflexion of male (5.619 years) were 98.64 mm for body length and 73.53 g for body weight.These results suggest that the female was isauxesis while the male was allometric, and that the male had a faster growth rate than female.

    Lake Nanwan; Pseudorasbora parva; Age; Growth

    Q178.1

    A

    1000-3207(2017)04-0835-08

    10.7541/2017.104

    2016-04-05;

    2016-12-14

    國家自然科學(xué)基金(1204306)資助 [Supported by the National Natural Science Foundation of China (1204306)]

    李紅敬(1966—), 男, 河南新野人; 博士; 研究方向為動物生態(tài)學(xué)。E-mail: lhj@xynu.edu.cn

    猜你喜歡
    雌雄體長雌性
    太行雞雌雄鑒別智能技術(shù)的研究與應(yīng)用
    脖子占體長一半的巨龍——新疆巨龍
    史上最長恐龍
    外泌體長鏈非編碼RNA在膀胱癌中的研究進(jìn)展
    連續(xù)超促排卵致腎精不足伴生育力低下雌性小鼠模型制備和比較研究
    河南一種雌性蚜蠅首次記述
    慢性焦慮刺激對成年雌性大鼠性激素水平的影響
    雌雄時代
    Coco薇(2015年12期)2015-12-10 02:40:50
    海州灣及鄰近海域大瀧六線魚的生長、死亡和最適開捕體長研究
    雙酚A對雌性生殖器官的影響及作用機(jī)制
    亚洲不卡免费看| 国产高清三级在线| 美女被艹到高潮喷水动态| 免费无遮挡裸体视频| 国产精品久久久久久亚洲av鲁大| 激情 狠狠 欧美| 毛片女人毛片| av在线老鸭窝| 高清午夜精品一区二区三区 | 久久精品综合一区二区三区| 成人漫画全彩无遮挡| 久久久午夜欧美精品| 色尼玛亚洲综合影院| 免费人成在线观看视频色| 久久久久免费精品人妻一区二区| 全区人妻精品视频| 日本与韩国留学比较| 久久久久久久久大av| 偷拍熟女少妇极品色| 99国产精品一区二区蜜桃av| 亚洲欧美清纯卡通| 亚洲四区av| 热99re8久久精品国产| 干丝袜人妻中文字幕| 高清毛片免费看| 99热这里只有是精品50| 尾随美女入室| 亚洲精品一卡2卡三卡4卡5卡| 国内精品美女久久久久久| 日本五十路高清| 草草在线视频免费看| 久久人人爽人人爽人人片va| 国产精品精品国产色婷婷| 97热精品久久久久久| 久久人人爽人人片av| 人妻久久中文字幕网| 精品乱码久久久久久99久播| 国产精品福利在线免费观看| 99久久精品热视频| 综合色丁香网| 亚洲精品国产av成人精品 | 国内精品久久久久精免费| 97超碰精品成人国产| 久久亚洲国产成人精品v| 国产亚洲精品综合一区在线观看| 亚洲av五月六月丁香网| 国产精品99久久久久久久久| av卡一久久| 久久久久久久久久黄片| 久久婷婷人人爽人人干人人爱| 国产精品av视频在线免费观看| 免费高清视频大片| 免费大片18禁| 亚洲精品影视一区二区三区av| 99在线视频只有这里精品首页| 日本黄色片子视频| 岛国在线免费视频观看| 久久久久九九精品影院| 热99re8久久精品国产| 国产91av在线免费观看| 日本撒尿小便嘘嘘汇集6| 午夜精品国产一区二区电影 | 成人欧美大片| 国产男靠女视频免费网站| 亚洲无线在线观看| 夜夜夜夜夜久久久久| 色5月婷婷丁香| 男人舔奶头视频| 久久久精品94久久精品| 国产精品一二三区在线看| 日日摸夜夜添夜夜添av毛片| avwww免费| 天天躁日日操中文字幕| 亚洲在线自拍视频| 国产精品久久久久久亚洲av鲁大| 丝袜美腿在线中文| 真实男女啪啪啪动态图| 又粗又爽又猛毛片免费看| 久久99热这里只有精品18| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品乱码久久久v下载方式| 精品久久久久久成人av| 99精品在免费线老司机午夜| 国产黄片美女视频| 91午夜精品亚洲一区二区三区| 亚洲精品成人久久久久久| 日本免费a在线| 一级av片app| 大香蕉久久网| 亚洲自偷自拍三级| 免费在线观看成人毛片| 少妇熟女aⅴ在线视频| 日韩中字成人| 午夜免费激情av| 天堂网av新在线| 久久久久久九九精品二区国产| 久久精品综合一区二区三区| 久久草成人影院| 久久鲁丝午夜福利片| eeuss影院久久| 亚洲高清免费不卡视频| 人妻久久中文字幕网| 天堂影院成人在线观看| 99九九线精品视频在线观看视频| av在线天堂中文字幕| 卡戴珊不雅视频在线播放| 国产黄色视频一区二区在线观看 | 国产精品免费一区二区三区在线| 欧美激情在线99| 国产成人aa在线观看| 又爽又黄a免费视频| 啦啦啦韩国在线观看视频| 国产探花极品一区二区| 国产老妇女一区| 美女免费视频网站| 成年女人永久免费观看视频| 中国美女看黄片| 色尼玛亚洲综合影院| 精品熟女少妇av免费看| 舔av片在线| 国产单亲对白刺激| 免费观看在线日韩| 欧美日韩乱码在线| 男女那种视频在线观看| 国产精品综合久久久久久久免费| 日本爱情动作片www.在线观看 | 精品久久久久久久久av| 老司机午夜福利在线观看视频| 久久久久国产精品人妻aⅴ院| 草草在线视频免费看| 亚洲四区av| 国产在线男女| 少妇丰满av| 高清毛片免费观看视频网站| 国产精品电影一区二区三区| 久久精品国产自在天天线| 97热精品久久久久久| 啦啦啦观看免费观看视频高清| 亚洲激情五月婷婷啪啪| 色播亚洲综合网| 国产aⅴ精品一区二区三区波| 1000部很黄的大片| 日本免费a在线| 亚洲国产精品国产精品| 神马国产精品三级电影在线观看| 一边摸一边抽搐一进一小说| 亚洲国产日韩欧美精品在线观看| 高清日韩中文字幕在线| 国产麻豆成人av免费视频| av在线亚洲专区| 国产高清激情床上av| 久久午夜亚洲精品久久| 国产精品乱码一区二三区的特点| АⅤ资源中文在线天堂| 九九爱精品视频在线观看| 国产成人精品久久久久久| 嫩草影院精品99| 女的被弄到高潮叫床怎么办| 悠悠久久av| 欧美日韩乱码在线| 日韩大尺度精品在线看网址| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 国产精品一区二区三区四区免费观看 | 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 日本精品一区二区三区蜜桃| 国产精品久久久久久久久免| 精品一区二区三区人妻视频| 日日摸夜夜添夜夜爱| 国产高清三级在线| 亚洲精品成人久久久久久| 美女cb高潮喷水在线观看| 久久久成人免费电影| 在线观看午夜福利视频| 特级一级黄色大片| 长腿黑丝高跟| 内地一区二区视频在线| 亚洲精品日韩在线中文字幕 | 尤物成人国产欧美一区二区三区| 人人妻人人澡欧美一区二区| 一级毛片我不卡| 赤兔流量卡办理| 男人舔女人下体高潮全视频| 少妇人妻精品综合一区二区 | 不卡一级毛片| 亚洲七黄色美女视频| av在线老鸭窝| 最新中文字幕久久久久| a级毛色黄片| 亚洲成人av在线免费| 五月玫瑰六月丁香| 亚洲精品乱码久久久v下载方式| 国产高清视频在线播放一区| 夜夜爽天天搞| 欧美日韩国产亚洲二区| 国产国拍精品亚洲av在线观看| 啦啦啦啦在线视频资源| 麻豆久久精品国产亚洲av| 国产精品女同一区二区软件| 成人综合一区亚洲| 色综合亚洲欧美另类图片| 能在线免费观看的黄片| 欧美区成人在线视频| 亚洲最大成人av| 午夜福利成人在线免费观看| 国产黄片美女视频| 99热网站在线观看| 狠狠狠狠99中文字幕| 欧美zozozo另类| 麻豆精品久久久久久蜜桃| 免费一级毛片在线播放高清视频| 22中文网久久字幕| 国产黄a三级三级三级人| 简卡轻食公司| 成人一区二区视频在线观看| 亚洲精品日韩在线中文字幕 | a级毛片a级免费在线| 桃色一区二区三区在线观看| 熟妇人妻久久中文字幕3abv| 成年免费大片在线观看| 国产精品乱码一区二三区的特点| 国产爱豆传媒在线观看| 99九九线精品视频在线观看视频| 久久99热6这里只有精品| 熟妇人妻久久中文字幕3abv| av中文乱码字幕在线| 日本免费a在线| 亚洲久久久久久中文字幕| 在线观看免费视频日本深夜| 精品午夜福利在线看| 国产视频内射| 97超碰精品成人国产| www日本黄色视频网| 日韩av在线大香蕉| 淫妇啪啪啪对白视频| 又粗又爽又猛毛片免费看| 女人十人毛片免费观看3o分钟| 国产精品永久免费网站| 在线天堂最新版资源| 免费高清视频大片| 又爽又黄无遮挡网站| h日本视频在线播放| 国产高清视频在线观看网站| 亚洲精品影视一区二区三区av| 女的被弄到高潮叫床怎么办| 两性午夜刺激爽爽歪歪视频在线观看| 美女免费视频网站| 国产精品久久久久久av不卡| 老师上课跳d突然被开到最大视频| 免费观看在线日韩| 综合色av麻豆| 国产单亲对白刺激| 国产高清不卡午夜福利| 久久久久久伊人网av| 午夜福利在线观看吧| 中文字幕久久专区| 国产三级在线视频| 精品久久久久久久末码| 天堂av国产一区二区熟女人妻| 亚洲精品日韩在线中文字幕 | av在线天堂中文字幕| 免费看日本二区| 你懂的网址亚洲精品在线观看 | 毛片女人毛片| 99热精品在线国产| 久久久久久久亚洲中文字幕| 亚洲色图av天堂| 国产高清视频在线观看网站| 桃色一区二区三区在线观看| 欧美不卡视频在线免费观看| 日本精品一区二区三区蜜桃| 国产亚洲欧美98| 亚洲性久久影院| 亚洲美女视频黄频| 熟女人妻精品中文字幕| 熟妇人妻久久中文字幕3abv| 国产蜜桃级精品一区二区三区| 亚洲av免费高清在线观看| 黄色欧美视频在线观看| 国产高清视频在线播放一区| 国产av不卡久久| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av| 我要搜黄色片| 真人做人爱边吃奶动态| 麻豆国产97在线/欧美| 干丝袜人妻中文字幕| 成年版毛片免费区| 国产男靠女视频免费网站| 精品午夜福利在线看| 久久久久九九精品影院| 老女人水多毛片| 少妇被粗大猛烈的视频| 特大巨黑吊av在线直播| 日本色播在线视频| 亚洲精品456在线播放app| 日韩av不卡免费在线播放| 国产亚洲精品久久久com| 大香蕉久久网| 久久6这里有精品| 亚洲成av人片在线播放无| 日本精品一区二区三区蜜桃| 国产成人a区在线观看| 男女做爰动态图高潮gif福利片| 久久99热6这里只有精品| 在线观看av片永久免费下载| 99久国产av精品国产电影| 变态另类成人亚洲欧美熟女| 中国美女看黄片| 色综合色国产| 蜜臀久久99精品久久宅男| 性色avwww在线观看| 看十八女毛片水多多多| 免费人成视频x8x8入口观看| 老熟妇乱子伦视频在线观看| 99国产精品一区二区蜜桃av| 国产午夜精品论理片| 99热网站在线观看| 级片在线观看| 日韩精品青青久久久久久| 99久国产av精品国产电影| 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 日本熟妇午夜| 成人午夜高清在线视频| 人妻丰满熟妇av一区二区三区| 精品人妻偷拍中文字幕| 十八禁网站免费在线| 精品一区二区三区视频在线| 老司机福利观看| .国产精品久久| 亚州av有码| 国产精品美女特级片免费视频播放器| 中国美白少妇内射xxxbb| 免费av毛片视频| 热99re8久久精品国产| 亚洲精品一卡2卡三卡4卡5卡| 国国产精品蜜臀av免费| 18禁在线播放成人免费| 国产精品嫩草影院av在线观看| 久久久久免费精品人妻一区二区| 久久99热这里只有精品18| 最近的中文字幕免费完整| 亚洲欧美日韩东京热| 91麻豆精品激情在线观看国产| 春色校园在线视频观看| 久久亚洲国产成人精品v| 亚洲精品成人久久久久久| 毛片一级片免费看久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产熟女欧美一区二区| 日韩制服骚丝袜av| 久久久久九九精品影院| 欧美日本视频| 日本黄色视频三级网站网址| 菩萨蛮人人尽说江南好唐韦庄 | 久久亚洲国产成人精品v| 99久久中文字幕三级久久日本| 欧美激情在线99| 国产成人a区在线观看| av国产免费在线观看| 99在线视频只有这里精品首页| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看 | 少妇熟女欧美另类| 日本黄色视频三级网站网址| 在线观看美女被高潮喷水网站| 老女人水多毛片| 少妇人妻一区二区三区视频| 好男人在线观看高清免费视频| 国产精品国产高清国产av| 国产精品日韩av在线免费观看| 在线播放国产精品三级| 亚洲美女黄片视频| 简卡轻食公司| 九九在线视频观看精品| 观看美女的网站| 伊人久久精品亚洲午夜| 免费搜索国产男女视频| 欧美日韩综合久久久久久| 高清毛片免费观看视频网站| 国产视频一区二区在线看| av在线播放精品| 国产精品人妻久久久久久| 久久中文看片网| aaaaa片日本免费| 三级男女做爰猛烈吃奶摸视频| 色哟哟·www| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 插阴视频在线观看视频| 中国美白少妇内射xxxbb| 国产精品国产三级国产av玫瑰| 国产精品久久电影中文字幕| 日韩一本色道免费dvd| 乱系列少妇在线播放| 日韩av不卡免费在线播放| 在现免费观看毛片| 国产成人freesex在线 | 日本免费一区二区三区高清不卡| 日韩成人av中文字幕在线观看 | 日日撸夜夜添| 韩国av在线不卡| 99热这里只有是精品50| 神马国产精品三级电影在线观看| 亚洲精品国产av成人精品 | 91在线精品国自产拍蜜月| a级一级毛片免费在线观看| 久久久久久久久久黄片| 欧美高清性xxxxhd video| 久久亚洲精品不卡| 99热这里只有是精品50| 亚洲无线在线观看| 日韩,欧美,国产一区二区三区 | 黄色视频,在线免费观看| 国产免费男女视频| 男女视频在线观看网站免费| 精品一区二区三区视频在线观看免费| 国产麻豆成人av免费视频| 亚洲av五月六月丁香网| 精品久久久久久久久av| 精品一区二区三区视频在线| 韩国av在线不卡| 不卡视频在线观看欧美| 少妇的逼好多水| 又爽又黄a免费视频| 日韩av不卡免费在线播放| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av涩爱 | 免费不卡的大黄色大毛片视频在线观看 | 午夜老司机福利剧场| 伊人久久精品亚洲午夜| av在线天堂中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美国产在线观看| 国产精品久久久久久亚洲av鲁大| 午夜福利视频1000在线观看| 99热6这里只有精品| 亚洲专区国产一区二区| 最近视频中文字幕2019在线8| 最近最新中文字幕大全电影3| 久久久成人免费电影| 精品久久久久久成人av| 97超视频在线观看视频| 午夜视频国产福利| 欧美高清性xxxxhd video| 色综合亚洲欧美另类图片| 天堂网av新在线| 露出奶头的视频| 国产成人freesex在线 | 中文字幕av在线有码专区| 国产黄色小视频在线观看| 免费av不卡在线播放| 卡戴珊不雅视频在线播放| 国产精品久久久久久久电影| 午夜激情欧美在线| 久久久久久九九精品二区国产| 一本一本综合久久| 国产乱人偷精品视频| 日韩,欧美,国产一区二区三区 | 大又大粗又爽又黄少妇毛片口| 久久精品国产自在天天线| 国产亚洲欧美98| 久久欧美精品欧美久久欧美| 性色avwww在线观看| 国产在线男女| 亚洲熟妇熟女久久| 成人美女网站在线观看视频| avwww免费| 亚洲av电影不卡..在线观看| 国产一区亚洲一区在线观看| 欧美日本视频| 国模一区二区三区四区视频| 大型黄色视频在线免费观看| 91精品国产九色| 菩萨蛮人人尽说江南好唐韦庄 | 三级毛片av免费| 亚洲精品在线观看二区| 特大巨黑吊av在线直播| 国产一区亚洲一区在线观看| 亚洲性久久影院| 中国美女看黄片| 天天躁夜夜躁狠狠久久av| 国产精品一二三区在线看| 成人三级黄色视频| 91精品国产九色| 国产精品,欧美在线| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放| 国产高清不卡午夜福利| 欧美xxxx性猛交bbbb| 色噜噜av男人的天堂激情| 亚洲人成网站高清观看| 99九九线精品视频在线观看视频| 亚洲专区国产一区二区| 狂野欧美白嫩少妇大欣赏| 18禁黄网站禁片免费观看直播| 亚洲图色成人| 亚洲无线在线观看| 在线播放无遮挡| 欧美性猛交╳xxx乱大交人| 免费看av在线观看网站| 看黄色毛片网站| 国产成年人精品一区二区| 亚洲不卡免费看| 欧美日韩在线观看h| 99热全是精品| 男人舔奶头视频| 欧美xxxx性猛交bbbb| 国产女主播在线喷水免费视频网站 | 国产蜜桃级精品一区二区三区| 免费一级毛片在线播放高清视频| 日韩av在线大香蕉| 欧美一区二区精品小视频在线| 亚洲欧美精品综合久久99| 亚洲av美国av| 欧美激情国产日韩精品一区| 国产亚洲av嫩草精品影院| 午夜福利在线在线| 天堂影院成人在线观看| 露出奶头的视频| 成人亚洲欧美一区二区av| 国产黄色视频一区二区在线观看 | 亚洲自拍偷在线| 久久久a久久爽久久v久久| 欧美性感艳星| 成年女人看的毛片在线观看| 色视频www国产| 日韩欧美三级三区| 久久精品国产亚洲av涩爱 | 成人美女网站在线观看视频| 一级毛片久久久久久久久女| 永久网站在线| 特大巨黑吊av在线直播| 精品久久久久久久久av| 网址你懂的国产日韩在线| 狂野欧美白嫩少妇大欣赏| 搡老岳熟女国产| 丝袜喷水一区| 日韩强制内射视频| 国产亚洲精品久久久久久毛片| 国产精品一二三区在线看| 丰满人妻一区二区三区视频av| 中国美白少妇内射xxxbb| 在线观看午夜福利视频| 欧美日本亚洲视频在线播放| 久久久精品大字幕| 99国产精品一区二区蜜桃av| 大型黄色视频在线免费观看| 国产乱人偷精品视频| a级毛片a级免费在线| 国产精品亚洲一级av第二区| 欧美3d第一页| 波多野结衣高清作品| 不卡视频在线观看欧美| 成人av在线播放网站| 亚洲美女视频黄频| 少妇的逼水好多| 国产精品久久视频播放| 日本与韩国留学比较| av黄色大香蕉| 综合色丁香网| av.在线天堂| 精品一区二区三区视频在线| 日本色播在线视频| 国产av麻豆久久久久久久| 69av精品久久久久久| 97碰自拍视频| 精品人妻一区二区三区麻豆 | 亚洲图色成人| 国产精品av视频在线免费观看| 成人精品一区二区免费| 淫秽高清视频在线观看| av.在线天堂| 免费看日本二区| 久久人人爽人人爽人人片va| 大香蕉久久网| 亚洲人与动物交配视频| 美女大奶头视频| 淫妇啪啪啪对白视频| av视频在线观看入口| 老熟妇仑乱视频hdxx| 免费黄网站久久成人精品| 最近视频中文字幕2019在线8| 精品一区二区三区视频在线观看免费| 人妻久久中文字幕网| 淫秽高清视频在线观看| 成年女人看的毛片在线观看| 国产91av在线免费观看| 乱人视频在线观看| 国产精品一区二区性色av| 校园春色视频在线观看| 中文字幕熟女人妻在线| 成年女人看的毛片在线观看| 精品99又大又爽又粗少妇毛片| 成人特级av手机在线观看| av天堂中文字幕网| 蜜桃亚洲精品一区二区三区| 99热只有精品国产| 午夜精品国产一区二区电影 | 国产精品一区二区三区四区久久| 国内精品一区二区在线观看| 99热6这里只有精品| 国产精品一区二区三区四区久久| 91在线观看av| 久久午夜亚洲精品久久| 99久久精品热视频| 一个人看视频在线观看www免费| 国产麻豆成人av免费视频| 黑人高潮一二区|