999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

岱衢洋魚類索餌場仔稚魚的生態位及其分化研究

2020-09-10 20:53:46樊紫薇蔣日進徐開達張琳琳印瑞王好學
海洋開發與管理 2020年4期

樊紫薇 蔣日進 徐開達 張琳琳 印瑞 王好學

摘要:為掌握重要魚類索餌場仔稚魚的生態分布,促進海洋漁業的發展,文章基于2011—2012年岱衢洋海域4個航次的定量采樣,運用生態位寬度、生態位重疊和典范對應分析等方法,分析該海域仔稚魚群落與海水環境因子的關系。研究結果表明:岱衢洋海域仔稚魚的生態位寬度可分為廣生態位、中生態位和窄生態位,其中廣生態位種對資源環境的適應能力強且分布范圍較廣,而窄生態位種對資源環境的適應能力較弱且分布不均勻;生態位寬度與相對重要性指數呈極顯著的正相關關系;仔稚魚群落生態位重疊指數大于0.6的種對占總種對數的24.24%,生態位重疊程度較低;仔稚魚分布存在明顯的生態分化現象,仔稚魚的生態位特征和生態習性受溫度、鹽度、懸浮物和葉綠素a等環境因子的影響較大。

關鍵詞:岱衢洋;仔稚魚;生態位寬度;生態位重疊;生態位分化

Abstract:Based on the quantitative sampling of 4 voyages in the Daiquyang from 2011 to 2012,using the niche breadth,niche overlap and canonical correspondence analysis,the relationship between the community of fish larvae and water environment factors was analyzed.The results showed that the niche breadth value of the fish larvae can be divided into niche breadth groups,middle niche breadth groups and narrow niche breadth groups.Among them,the species with wide niche which have strong adaptability to the environment and a wide range of distribution,narrow niche species have poor adaptability to the environment.Pearson correlation analysis showed that the niche breadth value was significantly correlated with the relative importance index.The pairs which the niche overlap above 0.6 accounting for 24.24% of the total number,and the degree of ecological overlap was low.The fish larvae showed obvious ecological differentiation.The niche characteristics changes and ecological habit of fish larvae are greatly affected by environmental factors such as temperature,salinity,suspended and chlorophyll a matter.

Key words:Daiquyang,Fish larvae,Niche breadth,Niche overlap,Niche differentiation

0 引言

生態位是指物種繁殖、發育和生存等資源環境特征的組合或排列[1],是生態學研究的重要思想之一[2]。有關生態位的理論已被廣泛應用于陸生生物群落的研究[3-8],近年來生態位的研究對象已逐步擴展到海洋生態系統。Korsu等[9]運用生態位測度分析芬蘭北部溪流土著種和入侵種魚類的競爭對棲息地資源環境的利用狀況;Galparsoro等[10]通過生態位因子分析法,對西班牙北部比斯開灣內的歐洲龍蝦(Homarus gammarus)的適宜棲息地模型進行預測,并繪制棲息地適宜性分布圖;于振海等[11]的研究表明,生態位寬度值可更好地反映物種生物量在不同尺度的分布狀況;何雄波等[12]認為閩江口主要魚類生態位的變動受季節變化的影響比較明顯;彭松耀等[13]運用生態位寬度和重疊,對南黃海大型底棲動物利用資源環境的狀況進行分析,表明物種生態位的動態特征與物種自身的生態習性和棲息環境有密切聯系。

上述研究表明,生態位理論在分析海洋生物對資源環境的利用狀況、種間競爭關系和種群群落地位等方面具有重要的應用意義和進一步的拓展空間。其中,浮游生物生態位可反映物種和種群在群落中的時空位置及其功能和地位[14-15],是理解浮游生物與生態環境協同適應的重要理論基礎,對于認知海洋生態系統的客觀發展規律具有深遠的意義,然而目前有關浮游生物生態位的研究較少[16-18]。

仔稚魚是海洋生態系統的重要組成部分,是魚類資源的補充群體,是海洋生物食物鏈中的主要餌料生物,在海洋生態系統物質能量傳遞中扮演關鍵角色,是漁業生物資源可持續利用的重要基礎[19-20]。同時,仔稚魚的數量分布和生長等與海洋環境變化密切相關。為更好地了解岱衢洋海域仔稚魚的資源生態狀況,本研究進行4個航次的定量采樣,綜合分析該海域魚類資源的早期補充機制以及仔稚魚的生態位及其分化特征,以期為該海域海洋生態系統的穩定和漁業水域環境的修復提供科學依據。

1 材料和方法

1.1 調查站位和樣品采集

岱衢洋位于舟山群島的長涂島和衢山島之間,優越的海洋環境造就著名的“岱衢族”大黃魚(Larimichthys crocea)的產卵場和索餌場,同時是帶魚(Trichiurus japonicus)等經濟魚類的棲息地[21]。

本研究在岱衢洋海域設置10個調查站位,分別于2011年5月(春)、8月(夏)、11月(秋)和2012年2月(冬)進行4個航次的仔稚魚調查(圖1)。

仔稚魚的采集方法依據《海洋調查規范:海洋生物調查》(GB/T 12763.6—1991)[22]。利用大型浮游生物網(口徑80 cm、長280 cm、網目尺寸0.505 cm)逐站進行表層水平拖網采樣,每站拖10 min,拖速為2 n mile/h,網口均配備水平流量計;采用多功能水質檢測儀(JFE-AAQ171),對海水的溫度(T)、鹽度(S)、懸浮物(SPM)、葉綠素a(Chl a)、活性硅酸鹽(As)和活性磷酸鹽(Ni)等相關環境因子數據進行同步采集;采集的仔稚魚樣品用5%福爾馬林溶液固定后帶回實驗室,根據仔稚魚的形態特征進行種類鑒定[23-24]并計數。

1.2 數據處理

1.2.1 優勢種

岱衢洋海域仔稚魚的優勢種采用相對重要性指數[25],計算公式為:

1.2.2 生態位寬度和生態位重疊

生態位寬度指數(Shannon-Weaver指數[26])的計算公式為:

生態位寬度指數(Bi)的取值范圍為0~R,其值越大,表示該物種的生態位寬度越大;生態位重疊指數(Qik)的取值范圍為0~1。

1.2.3 生態位分化

運用Canoco 5.0軟件分析仔稚魚與環境因子之間的關系并繪制排序圖,以揭示仔稚魚在該海域的生態位分化狀態。

對物種數據進行除趨勢對應分析(DCA),根據每個軸的梯度長度(LGA)選擇適宜的排序方法:當LGA<3時,選擇冗余分析(RDA);當LGA>4時,選擇典范對應分析(CCA);當3

采用Excel 2016軟件計算相對重要性指數,采用R語言計算生態位寬度指數和生態位重疊指數;采用Primer 5.0軟件對生態位寬度進行聚類分析,采用SPSS 17.0軟件進行Pearson相關性分析;在數據分析前對原始數據進行標準化處理。

2 結果

2.1 種類組成

2011年共鑒定出仔稚魚29種(不包括1個未定名物種),隸屬于9目17科,其中有2種鑒定到科。鱸形目種類最多,為7科15種;鯡形目種類次之,為1科4種。

本研究選取IRI>10的仔稚魚種群作為研究對象。其中,優勢種包括鳀(Engraulis japonicus)、中華小公魚(Stolephorus chinensis)和矛尾蝦虎魚(Chaeturichthys stigmatias)3種,常見種包括(Liza haematocheilus)、拉氏狼牙蝦虎魚(Odontamblyopus lacepedii)和竿蝦虎魚(Luciogobius guttatus)3種,一般種包括黃姑魚(Nibea albiflora)、帶魚(Trichiurus lepturus)、六絲鈍尾蝦虎魚(Amblychaeturichthys hexanema)、斑頭六線魚(Agrammus agrammus)、焦氏舌鰨(Cynoglossus joyneri)和前頜間銀魚(Hemisalanx prognathus)6種(表1)。

2.2 生態位寬度

根據生態位寬度可將仔稚魚分為3類,即廣生態位種(Bi≥1.0)、中生態位種(0.4≤Bi<1.0)和窄生態位種(0

由于黃姑魚的生態位寬度為0,無法進行聚類分析。對其余11種仔稚魚的生態位寬度進行聚類分析,可將其劃分為3組:①第一組為窄生態位種,包括竿蝦虎魚1種;②第二組為中生態位種,包括帶魚、斑頭六線魚和焦氏舌鰨3種;③第三組為廣生態位種,包括拉氏狼牙蝦虎魚、矛尾蝦虎魚、、鳀、中華小公魚、六絲鈍尾蝦虎魚和前頜間銀魚7種。各種類的生態位寬度越接近,其相似性也越接近(圖3)。

由圖2和圖3可以看出,仔稚魚種類的聚類分析結果與生態位寬度的分布相一致。

2.3 生態位寬度與相對重要性指數的關系

對岱衢洋海域仔稚魚的相對重要性指數和生態位寬度進行Pearson相關性分析,二者之間存在極顯著的正相關關系(F=0.741,P=0.006<0.01)(圖4)。

2.4 生態位重疊

岱衢洋海域仔稚魚的生態位重疊指數如表2所示。

由表2可以看出,岱衢洋海域仔稚魚物種之間的生態位重疊不均衡。仔稚魚種對共有66對,生態位重疊指數的范圍為0~0.968。其中:Qik>0.6的種對有16對,占總種對數的24.24%;0.3≤Qik≤0.6的種對有22對,占總種對數的33.33%;Qik<0.3的種對有28對,占總種對數的42.42%;重疊指數為0的種對有6對。黃姑魚(A4)和焦氏舌鰨(A11)的重疊指數最大(0.968),中華小公魚(A2)和帶魚(A5)的重疊指數其次(0.906),表明這些種類之間的資源環境利用能力具有相似性。總體而言,岱衢洋海域仔稚魚的生態位重疊程度較低。

2.5 生態位分化

岱衢洋海域仔稚魚的群落結構受懸浮物和活性磷酸鹽的影響較大。根據仔稚魚與環境因子的CCA分析結果,第一排序軸(CCA1,特征值為0.99)和第二排序軸(CCA2,特征值為0.85)表現出很高的物種-環境相關性(CCA1和CCA2與環境因子的相關性分別為0.99和0.93)。根據顯著性檢驗結果,第一排序軸和所有排序軸均呈現極顯著性差異(第一排序軸:F=1.3,P=0.001;所有排序軸:F=7.0,P=0.002)。其中,第一排序軸與懸浮物和活性磷酸鹽呈較大正相關關系,而與活性硅酸鹽和溫度呈較大負相關關系。中華小公魚、黃姑魚、帶魚、六絲鈍尾蝦虎魚、拉氏狼牙蝦虎魚和焦氏舌鰨等仔稚魚受溫度和鹽度的影響較大,而鳀、、矛尾蝦虎魚和竿蝦虎魚等仔稚魚受葉綠素a的影響較大。在排序圖中,受不同環境因子的影響,仔稚魚的群落分布呈明顯的兩極分化狀態(圖5)。

3 討論

3.1 生態位寬度

生態位是物種和種群對資源環境生態適應的結果,為某區域內物種群落的生態動態特征提供簡單和可靠的解釋[29]。生態位寬度是反映物種和種群利用資源環境多樣化程度的重要指標[30]。

鳀、中華小公魚和矛尾蝦虎魚是岱衢洋海域的優勢種,且為廣生態位種。一方面,該海域的魚類資源結構正向小個體、低營養級和低經濟價值的中上層魚類種群演替,群落中的優勢種如大黃魚、小黃魚和帶魚等經濟種類逐漸被鳀、小公魚和蝦虎魚等餌料性魚類所取代[31];另一方面,該海域位于亞熱帶季風性氣候區,同時受臺灣暖流的影響,優勢種為暖水種和暖溫種,其均為廣生態位種,分布均勻且活動范圍較大,這些種類的分布與其利用資源環境的能力和環境適應性有較大關系。

根據Pearson相關性分析,岱衢洋仔稚魚的相對重要性指數與生態位寬度存在極顯著的正相關關系,但二者并非為完全的正相關關系,如常見種和拉氏狼牙蝦虎魚的生態位寬度均大于優勢種鳀和中華小公魚,這可能與計算方法和考慮的生態特征有所不同有關。相對重要性指數與物種的豐度、生物量和出現站位占比等相關,但生態位寬度不僅考慮上述因素,而且綜合考慮不同站位之間物種的分布差異。

根據聚類分析,可將仔稚魚的生態位寬度分為3組,其中第三組的生態位寬度最大,主要包括鳀、、矛尾蝦虎魚和拉氏狼牙蝦虎魚等分布廣和環境適應性強的物種;第二組的生態位寬度較大,但斑頭六線魚、帶魚和焦氏舌鰨在各站位的數量分布不均勻;第一組的生態位寬度最小,這是由于近年來受環境惡化和過度捕撈的影響,產卵場破壞嚴重,黃姑魚的資源量減少[32],而竿蝦虎魚生長緩慢,物種在各站位的數量分布極不均勻。李顯森等[33]認為生態位寬度的聚類分析可反映群落中物種的分布格局,即生態位寬度越大,該物種的數量分布越均勻。這與本研究結果一致。

3.2 生態位重疊

生態位重疊體現物種對同一生境的利用和物種間的配置狀況,可反映物種對資源環境利用的相似程度和競爭程度[34]。根據對生態位寬度的分析,生態位寬度較大的物種對資源環境的利用能力和競爭能力較強,反之則較弱。然而本研究發現,生態位重疊指數最大的種對,其生態位寬度并不一定較大,例如:黃姑魚(A4)和焦氏舌鰨(A11)的生態位重疊指數最大(0.968),而前者的生態位寬度為0,后者為中生態位種。黃姑魚和焦氏舌鰨的生態位重疊指數最大的原因在于二者具有極為相似的棲息環境和相同的生態習性:二者均為暖溫性中下層魚類;黃姑魚棲息于泥沙底質較淺的沿岸海域,最適宜的生長水溫為24℃~28℃,產卵期為5—7月[35];焦氏舌鰨生活于砂泥底質海域,最適宜的生長水溫為24℃~28℃,魚卵和仔稚魚出現于5—9月[36]。此外,黃姑魚(A4)和拉氏狼牙蝦虎魚(A9)的生態位重疊指數為0,這是由于:①拉氏狼牙蝦虎魚棲息于水深為2~8 m的泥沙底質海域,于每年2—4月和7—9月在咸淡水水域產卵,與黃姑魚的產卵時間不同;②黃姑魚為洄游性種類,而拉氏狼牙蝦虎魚為地方性種類,二者生態習性不同。

同時,生態位重疊指數與物種分布環境位點的重合情況密切相關[37]。Wathne等[38]認為當生態位重疊指數大于0.6時,種間生態位重疊較顯著,且物種分布環境位點的重合度較高。岱衢洋海域仔稚魚的生態位重疊指數差異較大,物種之間存在不同程度的重疊:①鳀、中華小公魚、和矛尾蝦虎魚等廣生態位種與其他物種的生態位重疊指數大于0.6的概率較高;②生態位重疊指數為0的種對有6對,表明這些物種分布環境位點的差異較大,這與物種自身的生態特征和棲息環境等密切相關。

3.3 生態位分化

海洋生態系統與資源環境之間的相互作用基于多重的時間和空間尺度。相對于其他海洋生物群落而言,仔稚魚的自主運動能力較弱,具有隨波逐流的特點,因此環境因素對仔稚魚的生態分布發揮主導作用[17]。

沿第一排序軸,仔稚魚在溫度、鹽度、懸浮物和葉綠素a 4個梯度上的分化十分明顯:溫度是影響魚類游泳能力的重要因素,魚類通過調整自身的生理和生化過程以適應溫度的變化,從而影響其游泳能力[39];鹽度不僅影響魚類的代謝活動和滲透壓調節,而且影響魚類的性腺發育、組織結構和胚胎孵化等生理功能[40];葉綠素a的濃度與浮游植物的分布密切相關,而浮游植物是仔稚魚的餌料生物,從而間接影響仔稚魚的數量[41]。中華小公魚是中上層魚類,對溫度的反應非常敏感;黃姑魚和帶魚是洄游性種類,春、夏季沿岸海水升溫,魚類由深水區向近岸海域洄游,而秋、冬季則由近岸海域向深水區洄游;鳀和等在早期發育階段的食物來源與浮游植物密切相關,因此受葉綠素a濃度的影響較大。

典范對應分析(CCA)以環境因子為主體,通過排序來判斷物種在生境中的生態分化情況,可對生態位及其重疊無法合理解釋的結果做出補充說明[42]。例如:①黃姑魚(A4)和焦氏舌鰨(A11)對資源環境利用的相似程度較低,但二者均分布在溫度較高的海域,具有相似的資源環境需求,因此二者的生態位寬度不同而生態位重疊指數較大;②六絲鈍尾蝦虎魚(A6)分布于溫度和鹽度較高的海域,而前頜間銀魚(A12)受懸浮物的影響較大,因此二者雖均屬廣生態位種,但存在空間分布的差異,生態位重疊指數較小;③2011年岱衢洋海域水體呈中度污染和輕度富營養化狀態[43],因此大部分仔稚魚生活在活性硅酸鹽和活性磷酸鹽濃度較低的海域;④前頜間銀魚(A12)為廣生態位種,但與其他種類的生態位重疊指數均較低,CCA排序圖顯示其分布于懸浮物濃度較高的海域,而其他種類分布于懸浮物濃度較低的海域,生態適應性的差異使其存在空間分布的差異,并呈現生態位分化的特征;⑤在排序軸上分布差異較大的物種,其生態位重疊指數均較低。綜上所述,受不同環境因子的影響,仔稚魚的分布呈現明顯的生態分化現象。

參考文獻

[1] 鄭光美.鳥類學[M].北京:北京師范大學出版社,1995.

[2] 張金屯.數量生態學[M].北京:科學出版社,2004.

[3] CODY M.Niche theory and plant growth form[J].Vegetatio,1991,97:39-55.

[4] DECANS T,MARGERIE P,RENAULT J,et al.Niche overlap and species assemblage dynamics in an ageing pasture gradient in north-western France[J].Acta Oecologica,2011,37(3):212-219.

[5] 劉金福,洪偉.格氏栲群落生態學研究:格氏栲林主要種群生態位的研究[J].生態學報,1999,19(3):59-64.

[6] 蘇志堯,吳大榮,陳北光.粵北天然林優勢種群生態位研究[J].應用生態學報,2003,14(1):25-29.

[7] 劉貴河,王國杰,汪詩平,等.內蒙古典型草原主要草食動物食性及其營養生態位研究:以大針茅群落為例[J].草地學報,2013,21(3):439-445.

[8] 張晶,趙成章,任悅,等.張掖國家濕地公園優勢鳥類種群生態位研究[J].生態學報,2018,38(6):2213-2220.

[9] KORSU K,HEINO J,HUUSKO A,et al.Specific niche characteristics facilitate the invasion of an alien fish invader in boreal streams[J].International Journal of Ecology,2012,1:1-10.

[10] GALPARSORO I,BALD J,LIRIA P,et al.Predicting suitable habitat for the European lobster(Homarus gammarus),on the Basque continental shelf(Bay of Biscay),using Ecological-Niche Factor Analysis[J].Ecological Modelling,2009,220(4):556-567.

[11] 于振海,金顯仕,李顯森.黃海中南部主要魚種的生態位分析[J].漁業科學進展,2010,31(6):1-8.

[12] 何雄波,李軍,沈忱,等.閩江口主要漁獲魚類的生態位寬度與重疊[J].應用生態學報,2018,29(9):3085-3092.

[13] 彭松耀,李新正,王洪法,等.南黃海春季大型底棲動物優勢種生態位[J].生態學報,2015,35(6):1917-1928.

[14] ROSENTHAL G.Selecting target species to evaluate the success of wet grassland restoration[J].Agriculture,Ecosystems and Environment,2003,98(1/3):227-246.

[15] JORGE S A,TOWNSEND P.Interpretation of models of fundamental ecological niches and species′distributional areas[J].Biodiversity Informatics,2005,2:1-10.

[16] 徐曉群,曾江寧,陳全震,等.浙江三門灣浮游動物優勢種空間生態位[J].應用生態學報,2013,24(3):818-824.

[17] 梁淼,姜倩,孫麗艷,等.曹妃甸近岸海域大、中型浮游動物優勢種空間生態位研究[J].生態環境學報,2018,27(7):1241-1250.

[18] 李興,李建茹,李暢游.內蒙古烏梁素海浮游植物優勢種的生態位分析[J].水生態學雜志,2017,38(6):40-47.

[19] 趙靜,章守宇,汪振華,等.嵊泗人工魚礁海域魚類和大型無脊椎動物群落結構特征[J].中國水產科學,2010,17(5):1045-1056.

[20] 戴媛媛,侯純強,楊森,等.天津海域人工魚礁區浮游動物群落結構及其與環境因子的相關性研究[J].海洋湖沼通報,2018(5):163-170.

[21] 張洪亮,徐開達,朱增軍,等.岱衢洋魚類資源及其群落多樣性的季節變化[J].水產學報,2012,36(4):601-607.

[22] 國家技術監督局.海洋調查規范:海洋生物調查(GB/T 12763.6-1991)[S].北京:中國標準出版社,1992.

[23] 張仁齋,陸穗芳,趙傳絪,等.中國近海魚卵與仔魚[M].上海:上海科學技術出版社,1985.

[24] 萬瑞景,張仁齋.中國近海及其鄰近海域魚卵與仔稚魚[M].上海:上海科學技術出版社,2016.

[25] PINKAS L,OLIPHANT M S,IVERSON I L K.Food habits of albacore,Bluefin tuna,and bonito in California waters[J].Fish Bulletin,1971,1971(152):1-105.

[26] SHANNON C E,WEAVER W.The mathematical theory of communication[M].Chicago:University of Illinois Press,1963.

[27] PIANKA E R.The structure of lizard communities[J].Annual Review of Ecology and Systematics,1973,4:53-74.

[28] PETR J L.Multivariate analysis of ecological data using CANOCO[M].New York:Cambridge University Press,2003.

[29] HUTCHINSON G E.Cold spring harbor symposium of quantitative biology[J].Concluding Remarks,1957,22:415-427.

[30] 宋永昌.植被生態學[M].第二版.北京:高等教育出版社,2017.

[31] 趙淑江,呂寶強,李汝偉,等.物種滅絕背景下東海漁業資源衰退原因分析[J].中國科學:地球科學,2015,45(11):1628-1640.

[32] 逯云召,宓慧菁,于燕光,等.天津地區黃姑魚人工繁殖技術[J].河北漁業,2018(2):43-45,56.

[33] 李顯森,于振海,孫珊,等.長江口及其毗鄰海域魚類群落優勢種的生態位寬度與重疊[J].應用生態學報,2013,24(8):2353-2359.

[34] SPIES T A,FRANKLIN J F,KLOPSCH M.Canopy gaps in Douglas-fir forests of the Cascade Mountains[J].Canadian Journal of Forest Research,1990,20:649-658.

[35] 陳慶凱.黃姑魚規模化人工育苗及養殖技術的研究[D].廈門:集美大學,2016.

[36] 楊東萊,吳光宗,龐鴻艷.渤海灣的半滑舌鰨及焦氏舌鰨的魚卵和仔、稚魚的季節分布[J].生態學雜志,1984(3):30-33.

[37] 葛寶明,鮑毅新,鄭祥,等.靈昆島潮間帶大型底棲動物群落結構與生態位分析[J].生態學報,2005,25(11):245-251.

[38] WATHNE J A,HAUG T,LYDERSEN C.Prey preference and niche overlap of ringed seals Phoca hispida and harp seals P.groenlandica in the Barents Sea[J].Marine Ecology Progress Series,2000,194:233-239.

[39] 徐革鋒,牟振波,韓英,等.溫度和流速對細鱗鮭(Brachymystax lenok)幼魚游泳能力的影響及其與呼吸代謝的關系[J].海洋與湖沼,2015,46(4):786-792.

[40] 童燕,陳立僑,莊平,等.急性鹽度脅迫對施氏鱘的皮質醇、代謝反應及滲透調節的影響[J].水產學報,2007(S1):38-44.

[41] 高生泉,林以安,金明明,等.山東半島以南水域營養鹽的時空分布及其與魚產卵場的關系[J].海洋學報(中文版),2003(S2):157-166.

[42] 鄭挺,林元燒,曹文清,等.北部灣北部生態系統結構與功能:浮游動物空間生態位及其分化[J].生態學報,2014,34(13):3635-3649.

[43] 李振華,張亞洲,周永東,等.岱衢洋海域2010-2011年網采浮游植物群落特征與水質評價[J].浙江大學學報(理學版),2017,44(3):327-338.

主站蜘蛛池模板: 日本福利视频网站| 玖玖免费视频在线观看| 国产成人一区免费观看| 一本视频精品中文字幕| 日本人妻丰满熟妇区| P尤物久久99国产综合精品| 国产精品13页| 中国毛片网| 最新亚洲人成无码网站欣赏网| 狠狠综合久久久久综| 尤物国产在线| 无码一区中文字幕| 日本精品视频| 少妇高潮惨叫久久久久久| 99久久婷婷国产综合精| 亚洲中文久久精品无玛| 日韩二区三区无| 欧美一区二区自偷自拍视频| 国产欧美日韩免费| 一区二区三区精品视频在线观看| 爆操波多野结衣| 日韩a级片视频| 天天操天天噜| 日韩天堂视频| 福利在线一区| 亚洲欧美精品一中文字幕| 91久久精品国产| www.亚洲一区二区三区| 亚洲国产精品成人久久综合影院| 中文字幕伦视频| 国产精品hd在线播放| 亚洲第一视频区| 国产精品白浆无码流出在线看| 亚洲人成成无码网WWW| 国产精品人成在线播放| 亚洲综合二区| a天堂视频在线| 国产国产人成免费视频77777| 女人毛片a级大学毛片免费| 99视频在线精品免费观看6| 成人av手机在线观看| 久久99国产视频| 国产无码在线调教| 日韩毛片免费| 中文国产成人精品久久| 又猛又黄又爽无遮挡的视频网站| 国产色偷丝袜婷婷无码麻豆制服| 欧美精品亚洲日韩a| 人妻丰满熟妇AV无码区| a毛片在线播放| 国产精品久久久久久久久| 99在线观看国产| 国产后式a一视频| 亚洲AV无码乱码在线观看代蜜桃| 亚洲性影院| 高清码无在线看| 无码专区第一页| 国产中文在线亚洲精品官网| 少妇被粗大的猛烈进出免费视频| 男女男精品视频| 天天干天天色综合网| 日本一本在线视频| 红杏AV在线无码| 99er精品视频| 手机看片1024久久精品你懂的| 欧美精品高清| 欧美69视频在线| 欧美日韩福利| 国产黄在线免费观看| www中文字幕在线观看| 国产日本欧美亚洲精品视| 99热亚洲精品6码| 99尹人香蕉国产免费天天拍| 成人免费网站久久久| 一区二区影院| 无码精品国产dvd在线观看9久| 国产综合网站| 国产在线精品99一区不卡| 国产最新无码专区在线| 99无码中文字幕视频| 亚洲视频a| 在线免费不卡视频|