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

黃河三門峽庫區浮游植物群落結構空間分布特征及多樣性

2023-10-07 18:36:19單凱馬紅亮周子俊婁廣艷黃玉芳
人民黃河 2023年10期
關鍵詞:群落結構

單凱 馬紅亮 周子俊 婁廣艷 黃玉芳

摘 要:三門峽水庫是重要飲用水源地,為了解黃河三門峽水庫浮游植物情況,揭示水庫修建與運用對浮游植物群落特征的影響,進而判定其富營養化風險,于2019 年7—8 月、2020 年6—7 月對該庫區不同河段(自然河道段、過渡段、庫區段)的浮游植物群落結構進行研究,計算了浮游植物組成、密度、生物量、Margalef 豐富度指數、Shannon-Wiener 多樣性指數、Pielou 均勻度指數等。結果表明:三門峽庫區浮游植物共6 門39 種,群落組成以硅藻占絕對優勢,其次為綠藻、藍藻;浮游植物密度、生物量、Margalef 豐富度指數在不同河段均呈現顯著空間異質性,三門峽水庫修建與運行顯著改變了浮游植物群落特征;三門峽庫區水體總體為中營養,其中自然河道段為貧營養,過渡段和庫區段均為中營養。

關鍵詞:浮游植物;群落結構;生物多樣性;富營養化;三門峽庫區

中圖分類號:X522 文獻標志碼:A doi:10.3969/ j.issn.1000-1379.2023.10.014

引用格式:單凱,馬紅亮,周子俊,等.黃河三門峽庫區浮游植物群落結構空間分布特征及多樣性[ J].人民黃河,2023,45(10):78-83,89.

浮游植物是水生態系統的重要組成部分,是魚類等水生動物的天然餌料,在能量流動、物質循環和信息傳遞中具有重要作用[1] ,其種類組成及多樣性水平直接影響水生態系統的結構和功能。浮游植物具有生命力強、繁殖速度快、對環境變化反應靈敏等特點,其群落結構變化與水環境因子密切相關[2] ,因此常用于河流、湖庫等水體水質評價、預警等[3-5] 。此外,Margalef豐富度指數、Shannon-Wiener 多樣性指數、Pielou 均勻度指數、藻類密度等浮游植物群落參數常用于水生態狀況評價和富營養化評估[6-10] 。同時浮游植物群落組成可以指示湖庫的營養狀態,如硅藻在貧營養水體中占據優勢,藍藻、綠藻和隱藻門等在中營養湖庫中占有一定比例, 而綠藻和藍藻在富營養湖庫中占據優勢[11] 。

浮游植物群落結構、演替規律受諸多因素影響。不同地區、水域、季節引起浮游植物群落結構變化的環境因子差異顯著。水溫、溶解氧和高錳酸鹽指數是影響撫仙湖浮游植物群落結構的關鍵因子[12] ,浮游植物密度和生物量與水溫顯著正相關,與總磷、總氮、磷酸鹽、氨氮、硅酸鹽顯著負相關[5] 。三峽庫區干流浮游植物群落結構在魚類生活史不同時期對環境因子的響應不同,浮游植物群落在繁殖期主要受硝酸氮影響,在育肥期主要受水溫、總磷和懸浮物影響,在越冬期主要受氨氮、pH 值和溶解氧影響[13] 。除上述水環境因子外,水動力、水文、地形地貌、氣象條件等也是引起浮游植物群落結構變化的重要因素。笪文怡等[14] 認為對于新安江水庫這種大型貧-中營養水庫,浮游植物群落結構不僅受來水營養鹽負荷的影響,而且受水文和氣象條件影響。

三門峽水利樞紐是黃河干流上興建的第一座綜合性大型水利樞紐,對黃河水沙調節具有重要意義。三門峽水庫的建設運用,導致河段流速驟然減緩,河流的強烈淤積、侵蝕作用改變了原來的河谷形態[15] ,水生生物生境水動力和水文條件發生了較大變化,且庫區河段位于渭河、汾河等水質污染嚴重支流的入黃口下游,河段水質受到一定影響。目前,對三門峽庫區河段的研究多集中在濕地、魚類及泥沙等領域[16-17] ,對浮游生物尤其是浮游植物群落結構的研究鮮有報道,而水庫建設導致的水動力和水文條件變化對浮游植物群落結構變化的影響仍不甚明晰。因此,本文以三門峽庫區浮游植物為研究對象,將庫區河段分為自然河道段、過渡段和庫區段,分析浮游植物群落結構、生物多樣性及其分布特征,闡明水庫建設對浮游植物群落的影響,揭示水庫富營養化風險,以期為三門峽水庫富營養化控制和水源地保護提供參考。

1 采樣斷面設置及研究方法

三門峽水庫位于陜、晉、豫三省交界處,以防洪為主,兼顧防凌、灌溉、發電、供水,是黃河流域控制性骨干水庫之一。該水庫工程于1957 年4 月動工興建,1960 年9 月開始蓄水運用,庫區范圍涉及黃河龍門、渭河華縣、汾河河津和北洛河頭4 個水文站至大壩區間的干支流。1974 年開始采用“蓄清排渾” 運用,2003 年起非汛期運用水位不超過318 m,經常運用的庫區基本為潼關—三門峽大壩[18] ,本文以此為研究范圍。三門峽水庫是三門峽市主要水源地,供水主要滿足三門峽市工農業生產、生活用水需要,水庫水質變化直接影響三門峽市供水安全。

1.1 采樣斷面及采樣方法

2019 年7—8 月、2020 年6—7 月對三門峽庫區河段進行了浮游植物采樣調查工作。根據三門峽水庫非汛期運用水位不超過318 m 的控制條件,將庫區河段分為自然河道段(潼關—大禹渡)、過渡段(大禹渡—北村)和庫區段(北村—三門峽大壩)3 個部分,其中:自然河道段長45.64 km,主槽寬940 m、灘地寬2 390m、庫岸間距3 330 m,汛期沖刷,非汛期淤積,汛期小水時也發生淤積;過渡段長25.58 km,主槽寬690 m、灘地寬2 760 m、庫岸間距3 450 m,該河段在非汛期形成淤積三角洲;庫區段長42.28 km,主槽寬550 m、灘地寬1 500 m、庫岸間距2 050 m,該庫段已經形成高灘深槽的單一河勢,主河槽比較穩定,汛期中小洪水一般不出槽,非汛期受蓄水影響,風浪對高岸淘刷。選擇有代表性水域設置采樣斷面,共設置10 個采樣斷面,見表1、圖1。

浮游植物采集、鑒定按浮游植物常規調查方法進行。用25 號浮游植物網采集浮游植物,裝入100 mL標本瓶后加入1.5 mL 魯哥氏液固定,用于定性分析。此外,每個采樣斷面用采水器取水樣1 L,加入15 mL魯哥氏液固定,靜置沉淀48 h 后去除上清液,濃縮為30 mL 放入定量瓶保存,用于定量分析[19-20] 。采用0.1mL 吸管吸取0.1 mL 濃縮后的樣品于計數框內,蓋上蓋玻片,在光學顯微鏡(10 倍目鏡、40 倍物鏡)下采用目鏡視野法進行浮游植物計數與鑒定。每個樣品計數2 片,取平均值。生物量計算采用體積換算法[20] ,將每種浮游植物的個體數換算成生物量(浮游植物的體積質量接近1)。

2.2 浮游植物密度及生物量

三門峽庫區各采樣斷面浮游植物密度為9.65 萬~179.25 萬個/ L,平均密度為65.10 萬個/ L。整體上,三門峽庫區硅藻門密度占絕對優勢,其次為綠藻門、藍藻門。三門峽庫區浮游植物密度順河向逐步增大,且在空間分布上存在顯著差異,見圖3、表4。三門峽庫區浮游植物生物量為0.164 2~1.846 9 mg/ L,各采樣斷面平均生物量為0.758 1 mg/ L。整體上,三門峽庫區硅藻門、綠藻門生物量占優勢地位,其次為甲藻門。與密度的分布規律相似,浮游植物生物量順河向呈增大趨勢,其中庫區段生物量顯著高于自然河道段,見圖3、表4。從浮游植物密度對應的營養狀態來看,三門峽庫區總體上為中營養,其中自然河道段為貧營養,過渡段及庫區段均為中營養。

各采樣斷面硅藻門相對密度占比最大,平均為77.67%;綠藻門次之,平均為14.70%;藍藻門相對密度占比相對較小,為5.10%。從圖4 可以看出,自然河道段S1、S2、S3 斷面藍藻門相對密度占比高于綠藻門的,藍藻門與綠藻門的相對密度占比具有此消彼長的關系。

各采樣斷面硅藻門生物量占比最高,平均占比為77.87%;綠藻門次之,平均占比為13.49%;甲藻門相對較低,平均占比為5.84%。從圖5 可知,甲藻門在自然河道段沒有分布,在過渡段及庫區段出現并且生物量占比較大, S10 采樣斷面甲藻門生物量占比達25.05%。

2.3 浮游植物多樣性

三門峽庫區浮游植物多樣性指標見圖6、表5。三門峽庫區Margalef 豐富度指數Dm 為1.480~2.615,平均值為1.975。順河向Dm呈增大趨勢,過渡段及庫區段Dm顯著高于自然河道段Dm。從Margalef 豐富度指數看, 三門峽庫區各河段營養狀態均為富營養。

Shannon-Wiener 多樣性指數H′為1.987~2.531,平均值為2.216。從Shannon-Wiener 多樣性指數H′空間分布特征來看,與Dm變化趨勢略有不同,過渡段最大,自然河道段最小,但各采樣斷面Shannon-Wiener 多樣性指數H′并沒有顯著差異。從Shannon-Wiener 多樣性指數H′看,三門峽庫區各河段污染程度均為中污染,各河段營養狀態均為中營養。Pielou 均勻度指數J′為0.588~0.764,平均值為0.677,與Dm 的空間變化趨勢相反,順河向J′逐漸減小,但各采樣斷面之間J′并沒有顯著差異。從Pielou 均勻度指數J′看,三門峽庫區污染程度為輕污染或無污染,各河段營養狀態均為貧營養。

3 討論

三門峽庫區浮游植物群落種類組成以硅藻、綠藻為主,這與黃河流域龍羊峽水庫、李家峽水庫、積石峽水庫、小浪底水庫等河道型水庫,以及沙湖、烏梁素海、東平湖等自然湖泊浮游植物種類組成較為一致[22-26] 。

三門峽庫區綠藻和硅藻種類占比分別為46.15%、28.21%,種類組成表現為綠藻高于硅藻,與黃河自然河道相反(硅藻門占比為46. 78%,綠藻門占比為33.29%[27] )。與此類似,丁一桐等[27] 研究發現,黃河干流水庫影響河段表現出不同于自然河道的浮游植物群落特征。這可能歸因于水庫建設引起庫區水文、水動力以及水環境變化,進而影響藻類群落特征[28-30] 。

針對黃河多泥沙河流特征,相關研究表明水流急、泥沙量大的水體硅藻占比相對較大[22] 。三門峽水庫建設改變了河段原有的河流形態,受水庫回水頂托,三門峽庫區水體流速變緩[31] ,水體中泥沙沉積導致庫區水體含沙量較自然河道顯著降低,透明度高的靜水適宜綠藻生存,從而使三門峽庫區綠藻種類占比高于硅藻。

從密度及生物量來看,三門峽庫區浮游植物平均密度為65.10 萬個/ L,平均生物量為0.758 1 mg/ L,遠低于沙湖[24] 、烏梁素海[25] 、東平湖[32] 的。浮游植物密度及生物量的增大是水體營養狀態的表現形式之一。一般而言,水體富營養化程度越高,浮游植物密度及生物量越大[32] 。與黃河干流高海拔地區的龍羊峽水庫、李家峽水庫、積石峽水庫[22] 相比,三門峽庫區浮游植物密度和生物量均較高,主要原因是黃河龍羊峽至積石峽段位于黃河上游高海拔地區,氣候寒冷,水溫低,抑制了浮游植物的生長繁殖[22,33] ,而三門峽庫區位于黃河中下游河段,海拔低,水溫相對較高,氮、磷等營養鹽濃度較高[34-35] ,有利于浮游植物的生長繁殖。

三門峽水庫的修建顯著改變了原有河道水動力、水文情勢以及水環境等條件,非汛期水位不超過318 m 的運用模式將庫區河段分為自然河道段、過渡段和庫區段,其浮游植物群落在上述河段呈現出顯著的空間異質性。具體而言,三門峽庫區浮游植物密度、生物量及物種豐富度指數Dm沿河流方向均具有增大趨勢,其中自然河道段和庫區段分布具有顯著的空間差異性。此外,在水體光照條件較差時,即使其他水體條件均有利于浮游植物的生存,浮游植物的生物量仍然可能下降[36] 。三門峽庫區這種空間分異性的產生可能與庫區水文、水動力條件以及泥沙條件有關。三門峽庫區自然河道段水文狀況類似于自然河流水體,水流急、泥沙含量相對較大,其生態系統接近自然河流生態系統,具有浮游植物生物量低的特點[37] 。庫區段水體類似于靜水湖泊,流速極為緩慢,水體因泥沙沉積而透明度大,浮游植物能夠利用的有效光源多[38] ,穩定的環境條件適宜藻類生長[22] 。總體而言,庫區蓄水對浮游植物生長具有促進作用,三門峽水庫的修建及運用顯著改變了浮游植物的群落結構。

根據浮游植物密度對應的營養狀態,三門峽庫區總體為中營養,其中自然河道段為貧營養,過渡段和庫區段均為中營養,主要原因可能是,三門峽庫區3 個分段的水體流速、水溫、含沙量等不同,順河向流速明顯減緩,污染物質滯留時間延長[13] 導致一定富營養化風險。根據生物多樣性指數判斷,三門峽庫區為富營養(根據Dm值),雖然多樣性指數因其所反映群落信息變化的不確定性而不能準確描述富營養化過程[39] ,但表明三門峽庫區有出現富營養化的風險。三門峽庫區Shannon-Wiener 物種多樣性指數H′表明庫區水質污染程度為中污染,Pielou 均勻度指數J′表明庫區水質污染程度為輕污染或無污染,即三門峽庫區水質尚未達到重污染的程度。但根據歷史監測資料,三門峽庫區污染物中氨氮和總磷超標情況較為多見[34] ,部分季節甚至出現Ⅴ類和劣Ⅴ類水,說明采用生物多樣性指數評價三門峽庫區水質存在局限性。

4 結論

三門峽庫區浮游植物種類以綠藻和硅藻為主,在組成上與黃河干流自然河道水體不同,主要表現為三門峽庫區硅藻種類占比小于綠藻,而自然河道水體硅藻種類多于綠藻。三門峽水庫修建與運行顯著改變了原有的水動力條件、水文情勢以及水環境,浮游植物密度及生物量等群落結構特征發生顯著變化,表現為庫區蓄水對浮游植物生長具有促進作用,為魚類提供天然餌料資源。根據浮游植物密度對應的營養狀態,三門峽庫區總體為中營養,其中自然河道段為貧營養,過渡段和庫區段均為中營養,這可能與自然河道段水體流速較快、含沙量高有關。其他多樣性指標表明在三門峽水庫的管理工作中需進一步關注庫區富營養化風險。

參考文獻:

[1] 鄭娜,楊丹丹,蒲松,等.沙壩水庫浮游生物資源調查研究[J].現代農業科技,2021(7):161-166,169.

[2] 周緒申,李娜,孫博聞,等.白洋淀浮游生物群落結構的季節變化及其與環境因子的關系[J].水利水電技術(中英文),2021,52(8):110-119.

[3] KAMENIR Y,DUBINSKY Z,ZOHARY T.Phytoplankton SizeStructure Stability in a Meso?Eutrophic Subtropical Lake[J].Hydrobiologia,2004,520(1-3):89-104.

[4] WU N C,SCHMALZ B,FOHRER N.Development and Testingof a Phytoplankton Index of Biotic Integrity (P-IBI) for aGerman Lowland River [J]. Ecological Indicators, 2011,13(1):158-167.

[5] 余業鑫,李艷,向羅京,等.漢江下游干支流浮游植物群落特征及其對水質的指示評價[J].中國環境監測,2022,38(1):124-135.

[6] 況琪軍,馬沛明,胡征宇,等.湖泊富營養化的藻類生物學評價與治理研究進展[J].安全與環境學報,2005,5(2):87-91.

[7] SABATER S,ARTIGAS J,DURDN C,et al.Longitudinal De?velopment of Chlorophyll and Phytoplankton Assemblages ina Regulated Large River (the Ebro River)[J].Science of theTotal Environment,2008,404(1):196-206.

[8] BORICS G,V?RB?R? G,GRIGORSZKY I,et al.A New E?valuation Technique of Potamo?Plankton for the Assessmentof the Ecological Status of Rivers[J].Large Rivers,2007,17(3-4):465-486.

[9] CELLAMARE M,MORIN S,COSTE M,et al.Ecological As? sessment of French Atlantic Lakes Based on Phytoplankton,Phytob?Enthos and Macrophytes[J].Environmental Monitoringand Assessment,2012,184(8):4685-4708.

[10] 錢奎梅,陳宇煒,宋曉蘭.太湖浮游植物優勢種長期演化與富營養化進程的關系[J].生態科學,2008,27(2):65-70.

[11] 屈寧,鄧建明,張禎,等.2015—2020 年洪澤湖浮游植物群落結構及其環境影響因子[J].環境科學,2022,43(6):3097-3105.

[12] 吉正元,劉紹俊,施艷峰,等.撫仙湖浮游植物群落結構周年變化及水質生物評價[J].海洋湖沼通報,2019(3):78-86.

[13] 魏念,余麗梅,杜開開,等.三峽庫區干流浮游植物群落結構及環境影響因子相關性分析[J].長江流域資源與環境,2022,31(3):615-623.

[14] 笪文怡,朱廣偉,吳志旭,等.2002—2017 年千島湖浮游植物群落結構變化及其影響因素[J].湖泊科學,2019,31(5):1320-1333.

[15] 談皓,潘明強,李繼偉,等.黃河潼關至三門峽大壩河段治理措施探討[J].水利與建筑工程學報,2007,5(4):83-86.

[16] 張軍燕,張建軍,沈紅保,等.小浪底水庫調水調沙對壺口至三門峽段魚類群落結構的影響[J].生態學雜志,2012,31(10):2613-2618.

[17] 王富強,張紅璐,趙衡,等.三門峽庫區濕地水資源利用效用評價[J].人民黃河,2021,43(3):69-73,101.

[18] 吳保生,鄭珊,沈逸.三門峽水庫沖淤與“318 運用”的影響[J].水利水電技術,2020,51(11):1-12.

[19] 章宗涉,黃祥飛.淡水浮游生物研究方法[M].北京:科學出版社,1991:340-344.

[20] 中華人民共和國農業部.淡水浮游生物調查技術規范:SC/ T 9402—2010[S].北京:中國農業出版社,2011:4-6.

[21] 王朝暉,韓博平,胡韌,等.廣東省典型水庫浮游植物群落特征與富營養化研究[J].生態學雜志,2005,24(4):402-405,409.

[22] 陳燕琴,申志新,劉玉婷.黃河上游龍羊峽至積石峽段春秋季浮游植物群落結構及多樣性評價[J].水生態學雜志,2014,35(5):52-59.

[23] 翟江輝.小浪底水庫水質和微生物多樣性特征研究[D].鄭州:河南農業大學,2012:37-41.

[24] 邱小琮,趙紅雪.寧夏沙湖浮游植物群落結構及多樣性研究[J].水生態學雜志,2011,32(1):20-26.

[25] 金業.烏梁素海浮游植物功能群時空變異性及其影響因子研究[D].呼和浩特:內蒙古大學,2019:17-22.

[26] 董貫倉,冷春梅,叢旭日,等.南水北調東線工程運行3年后東平湖浮游植物群落特征及環境驅動因子[J].湖泊科學,2022,34(1):61-73.

[27] 丁一桐,潘保柱,趙耿楠,等.黃河干流全河段浮游植物群落特征與水質生物評價[J].中國環境科學,2021,41(2):891-901.

[28] 王孟,鄔紅娟,馬經安,等.長江流域大型水庫富營養化特征及成因分析[J].長江流域資源與環境,2004,13(5):477-481.

[29] LANGE T R,RADA R G.Community Dynamics of Phyto?plankton in Atypical Navigation Pool in the Upper Missis?sippi River [ J]. The Journal of the Iowa Academy ofScience,1993,100(1):21-27.

[30] HUFF D R.Phytoplankton Communities in Navigation PoolNo.7 of the Upper Mississippi River [J]. Hydrobiologia,1986,136(1):47-56.

[31] 金雙彥,張堯旺,王劍非.三門峽庫區淤積物粒徑分析[J].中國水土保持,2006(8):16-18.

[32] 陳璽,曹秀云,周易勇.小浪底水庫與東平湖浮游植物分布的異質性[J].環境科學與技術,2013,36(增刊1):50-59.

[33] 王恒山,劉信利,王玫.青海省黃河上游羊曲段水生生物調查研究[J].河北漁業,2016(9):17-21.

[34] 劉茹,屠新武.黃河三門峽庫區水質評價與演變勢分析[J].西北水力發電,2005,21(1):55-59.

[35] 馮健,周懷東,彭文啟,等.三門峽水庫水質模型與水環境容量分析[J].人民黃河,2014,36(3):35-36,39.

[36] ADMIRAAL W,BREEBAART L,TUBBING G M J,et al.Plankton in the River Rhine: Structural and FunctionalChanges Observed During Downstream Transport [ J ].Journal of Plankton Research,1992,14(10):1351-1368.

[37] REYNOLDS C S,DESCY J P,PADIS?K J.Are Phytoplank?ton Dynamics in Rivers so Different from Those in ShallowLake? [J].Hydrobiologia,1994,289(1-3):1-7.

[38] 何書晗,歐陽添,趙璐,等.三峽庫區支流浮游植物群落穩定性及其驅動因子分析[J].環境科學,2021,42(7):3242-3252.

[39] 劉書宇,馬放,張建祺.景觀水體富營養化模擬過程中藻類演替及多樣性指數研究[J].環境科學學報,2007,27(2):337-341.

【責任編輯 呂艷梅】

猜你喜歡
群落結構
海南省紅樹林區浮游動物多樣性的初步研究
海南省紅樹林區大型底棲動物多樣性的初步研究
長豐鰱苗種培育對池塘浮游動物群落結構的影響
河北漁業(2015年12期)2015-12-21 14:37:31
濟南流域夏季大型底棲動物群落結構調查及水質評價
河北漁業(2015年11期)2015-11-18 16:39:54
渭河陜西段浮游植物群落結構及水質評價
河北漁業(2015年11期)2015-11-18 16:37:49
濟南流域春季浮游植物多樣性的研究
河北漁業(2015年10期)2015-10-15 11:14:57
貴州草海濕地大型底棲動物群落結構及水質生物學評價
濟南流域輪蟲群落結構及多樣性的研究
河北漁業(2015年7期)2015-07-21 11:06:52
黃河蘭州市區段浮游動物群落結構調查及水質初步評價
河北漁業(2015年6期)2015-06-24 11:04:37
安徽瓦埠湖湖濱帶植物群落結構研究
主站蜘蛛池模板: 亚洲最大在线观看| 久久国产高清视频| AV不卡在线永久免费观看| 国产一级二级三级毛片| 人禽伦免费交视频网页播放| 国产精品私拍在线爆乳| 久久6免费视频| 日韩天堂视频| 久爱午夜精品免费视频| 尤物亚洲最大AV无码网站| 毛片免费在线| 国产AV毛片| 国产精品高清国产三级囯产AV| 国产成人综合日韩精品无码首页 | 91区国产福利在线观看午夜| 精品人妻无码中字系列| 欧美啪啪精品| 国产香蕉97碰碰视频VA碰碰看| 国产精品久久久久鬼色| 国产又色又刺激高潮免费看| 亚洲高清在线天堂精品| 精品久久久无码专区中文字幕| 久久免费视频6| 欧美一区国产| 精品色综合| 国产亚洲精| 久久久久久久蜜桃| 六月婷婷综合| 中文无码毛片又爽又刺激| 狠狠色噜噜狠狠狠狠色综合久| 国内精品九九久久久精品| 亚洲国产天堂久久综合226114| 67194亚洲无码| 亚洲欧美日韩视频一区| 免费欧美一级| 国产成人精品一区二区秒拍1o| 久久精品一卡日本电影| 黄色三级网站免费| AV老司机AV天堂| 呦系列视频一区二区三区| 亚洲一级色| 欧美日韩国产成人高清视频| 在线中文字幕日韩| 国产精品污污在线观看网站| 玖玖精品在线| 亚洲黄色成人| 人人91人人澡人人妻人人爽 | 亚洲精品无码专区在线观看| 久久久久久久久久国产精品| 无码丝袜人妻| 亚洲欧美日韩高清综合678| 人妻无码一区二区视频| 欧美日韩第二页| 98精品全国免费观看视频| 欧美精品1区| 亚洲无码高清视频在线观看| 99中文字幕亚洲一区二区| 国产浮力第一页永久地址| 亚洲天堂.com| 97人人做人人爽香蕉精品| 亚洲一区二区三区麻豆| 波多野衣结在线精品二区| 欧美性爱精品一区二区三区 | 中文字幕在线看| 国产jizzjizz视频| 亚洲人在线| 国产va免费精品观看| 国产va在线观看免费| 亚洲一级毛片在线观| 波多野结衣一二三| 欧美亚洲国产一区| 婷婷亚洲综合五月天在线| 无码电影在线观看| 国产精品区网红主播在线观看| 色婷婷成人| 永久免费AⅤ无码网站在线观看| 亚洲人成成无码网WWW| 久久久久人妻一区精品色奶水 | 97国产一区二区精品久久呦| 国产欧美在线| 亚洲最大看欧美片网站地址| 国产精品亚洲精品爽爽|