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

東平湖酶活性與溫室氣體通量及其影響因素

2022-02-13 11:50:04劉濤曹起孟張菊楊麗偉張智博鄧煥廣
人民黃河 2022年2期
關鍵詞:影響因素

劉濤 曹起孟 張菊 楊麗偉 張智博 鄧煥廣

摘 要:為研究東平湖沉積物酶活性與水體溫室氣體排放通量特征及其影響因素,2018年9月布設23個采樣點,采集表層沉積物和表層水樣品,分析沉積物中脲酶、蔗糖酶和堿性磷酸酶活性,采用頂空法和雙層擴散模型估算水氣界面N2O、CO2排放通量;采用相關分析、多元線性回歸法探討影響酶活性與溫室氣體排放通量的主要因子。結果表明:沉積物脲酶、蔗糖酶和堿性磷酸酶活性平均值分別為1.379 mg/(g·d)、22.633 μg/(g·d)和0.307 mg/(g·d),N2O、CO2排放通量分別為1.18 μmol/(m2·h)、19.91 mmol/(m2·h);酶活性在湖區東部(老湖碼頭)較高、在大汶河入湖口和湖心區較低,而N2O、CO2排放通量在大汶河入湖口較高、出湖口較低;沉積物總有機碳與3種酶活性均極顯著正相關,是影響酶活性的主要因子之一;酶活性與影響因子的回歸模型可分別解釋脲酶、蔗糖酶和堿性磷酸酶活性變異程度的43.6%、73.4%和88.8%;溫室氣體排放通量僅與表層水體部分理化指標顯著相關,表層水體無機氮濃度、pH值和DO是影響N2O排放通量的主要因子,而CO2排放通量主要受水深的影響。

關鍵詞:酶活性;溫室氣體;排放通量;影響因素;東平湖

中圖分類號:X524;TV211 文獻標志碼:A

doi:10.3969/j.issn.1000-1379.2022.02.021

引用格式:劉濤,曹起孟,張菊,等.東平湖酶活性與溫室氣體通量及其影響因素[J].人民黃河,2022,44(2):106-111.

Abstract: In order to study the characteristics and influencing factors of sediment enzyme activity and greenhouse gas emission flux in Dongping Lake, surface sediment and water samples were collected simultaneously at 23 sampling sites around the lake in September 2018. The activities of urease, sucrase and alkaline phosphatase in sediments were analyzed, and the N2O and CO2 emission fluxes at the water-air interface were analyzed and estimated by the headspace method and double-layered diffusion model. The main influencing factors were explored by conducting Pearson correlation analysis and multiple linear regression. The results show that the activities of urease, sucrase and alkaline phosphatase are 1.379 mg/(g·d), 22.633 μg/(g·d) and 0.307 mg/(g·d) respectively. The fluxes of N2O and CO2 are 1.18 μmol/(m2·h) and 19.91 mmol/(m2·h) respectively. The sediment enzyme activity is significantly higher in the east of the lake (near Laohu Wharf), and lower in the inlet of Dawen River and central area of the lake;while the N2O and CO2 emission fluxes are higher in the inlet of Dawen River and lower in the outlet of the lake. TOC concentrations in the sediments are significantly positively correlated with the enzyme activities, indicating that sediment TOC is one of the main influencing factors. The regression models of enzyme activity and the main influencing factors can explain 43.6%, 73.4% and 88.8% of the variation degree of urease, sucrase and alkaline phosphatase activity respectively. N2O and CO2 emission fluxes are only significantly correlated with part of physicochemical indexes of the surface water;and specifically, inorganic nitrogen concentrations, water pH and DO are the main factors that affecting N2O emission fluxes, while CO2 emission fluxes are mainly affected by water depth.

Key words: enzyme activity;greenhouse gas;emission flux;influencing factor;Dongping Lake

湖泊濕地分布廣泛、水文條件特殊,在陸地生態系統碳氮循環中具有重要作用,是溫室氣體的“源”與“匯”[1]。沉積物酶主要來自微生物細胞,在調節沉積物碳、氮、磷等元素生物地球化學循環中發揮著關鍵作用,酶活性可直接表征沉積物微生物活性,也可以作為人類活動影響下沉積物營養物質含量的敏感指示因子[2-3]。目前,有關沉積物酶活性時空分布特征及其與營養鹽和水體富營養化關系的研究較多[2,4-5],但對沉積物酶活性和溫室氣體排放及其影響因素的研究較少[6]。

東平湖位于山東省西南部,常年水面面積為124 km2,是山東省第二大淡水湖泊[7],也是南水北調東線工程的主要調蓄湖和山東省西水東送工程的水源地。目前東平湖整體呈中營養至富營養狀態,隨著對東平湖外源污染的控制和調水工程的實施,內源污染物釋放成為影響東平湖水質的主要因素[8-9]。已有研究多關注東平湖水質、沉積物分布等[8,10-14],對于影響水體碳、氮、磷等元素循環的因子如微生物活性,以及反映其環境效應的指標如溫室氣體排放通量等的研究相對較少[15-16]。研究通過對東平湖沉積物酶活性和水氣界面溫室氣體排放通量的測定,分析其空間分布特征,并根據二者與水體理化性質的統計學分析,探討影響東平湖沉積物酶活性和溫室氣體排放通量的主要環境因子,以期為進一步研究微生物活性在湖泊濕地生物地球化學循環中的作用和東平湖水生態環境保護提供理論依據。

1 材料與方法

1.1 樣品采集和預處理

依據《湖泊沉積物調查規范》[17],按照均勻性和可達性原則在東平湖設置23個采樣點,見圖1。2018年9月在各采樣點使用卡蓋式采水器采集表層水樣,一部分置于500 mL聚乙烯塑料瓶內放入便攜式冷藏箱中保存,另一部分置于有機玻璃管密封后用于溫室氣體質量濃度的測定[18]。使用溶解氧儀現場測定水溫和溶解氧(DO)含量,采用便攜式氧化還原電位儀測定水體氧化還原電位(Eh),采用塞氏盤測定水體透明度,采用自制水深測定裝置測定水深。使用一次性注射器采集空氣樣品,儲存于氣體采樣袋;采用抓斗式采泥器采集表層(0~10 cm)沉積物置于自封袋中放入冷藏箱保存。采集樣品的同時記錄各采樣點周圍環境信息,并用GPS定位。

實驗室采用頂空法預處理有機玻璃管水樣后測定CO2和N2O的質量濃度[18]。取部分水樣過0.45 μm濾膜置于聚乙烯塑料瓶內,用于氨氮(NH3-N)、硝態氮(NO-3-N)、亞硝態氮(NO-2-N)和溶解性有機碳的測定,濾膜用于葉綠素a的測定。沉積物于陰涼通風處自然風干后研磨過60目篩備用。

1.2 樣品的分析和測定

取適量過篩后沉積物樣品測定酶活性[19],蔗糖酶(sucrase)、脲酶(urease)、堿性磷酸酶(alkaline phosphatase)活性分別采用3,5-二硝基水楊酸比色法、苯酚鈉-次氯酸鈉比色法、磷酸苯二鈉比色法測定。氣體中CO2、N2O質量濃度采用氣相色譜儀測定。水體中NH3-N、NO-3-N、NO-2-N、葉綠素a質量濃度分別采用納氏試劑分光光度法、紫外分光光度法、N-(1-萘基)-乙二胺分光光度法和丙酮分光光度法測定[20],水中總有機碳采用總有機碳分析儀測定。沉積物中總磷(TP)和總有機碳分別采用酸溶-鉬銻抗比色法[21]和重鉻酸鉀氧化外加熱法測定。

1.3 數據處理與統計

東平湖水面寬闊,流速緩慢,滿足雙層擴散模型的假設條件,因此采用雙層模型計算水氣界面CO2和N2O的排放通量[22]。數據的統計分析采用Excel 2019和SPSS 25.0,并利用ArcGIS10.6進行克里金插值分析和繪圖。

2 結 果

2.1 東平湖水體基本理化性質

東平湖表層水體DO、NH3-N質量濃度均能達到《地表水環境質量標準》(GB 3838—2002)Ⅲ類水標準(見表1),但總無機氮平均質量濃度為1.174 mg/L,有73.9%的點位總無機氮質量濃度超過Ⅲ類水標準。變異系數(CV)是反映變異程度的一個統計量,根據Wilding[23]對變異程度的分類,表層水體pH值、氧化還原電位、溶解性有機碳和沉積物pH值、TP為小變異(CV<15%),表層水體DO、透明度、NH3-N、沉積物氧化還原電位和水深為中等變異(15%≤CV≤36%),而表層水體葉綠素a、NO-3-N、NO-2-N和沉積物總有機碳為高度變異(CV>36%)。水體無機氮質量濃度與田翠翠等[24]2016年10月的調查結果基本一致,沉積物總磷質量濃度低于田忠志等[10]的研究結果(平均值為0.795 mg/g),一定程度上反映了近10 a來東平湖水體磷負荷降低,而沉積物總有機碳的含量和范圍則與張智博等[11-12]的研究結果較為一致。

2.2 東平湖沉積物酶活性與溫室氣體排放通量特征

東平湖沉積物脲酶、蔗糖酶和堿性磷酸酶活性的平均值分別為1.379 mg/(g·d)、22.633 μg/(g·d)和0.307 mg/(g·d),其空間差異均較為顯著(CV>36%),且以蔗糖酶的變異系數最大,見表2。東平湖水氣界面N2O和CO2排放通量的平均值分別為1.18 μmol/(m2·h)和19.91 mmol/(m2·h),N2O排放通量空間變異大(CV>36%),而CO2排放通量空間變異小(CV<15%)。

進一步對東平湖沉積物酶活性和溫室氣體排放通量的空間分布特征進行分析,結果見圖2。3種酶活性的空間分布相似,均表現為湖區東部(老湖碼頭附近水域)酶活性顯著高于其他水域,此外脲酶和堿性磷酸酶活性在出湖口較高。然而各酶活性低值的分布有明顯差異,脲酶活性在大汶河入湖口附近最低,而蔗糖酶活性在湖區西部和湖心區處于較低水平,湖心區堿性磷酸酶活性明顯較低,此外,3種酶活性的最低值均出現在大汶河入湖口附近。N2O和CO2排放通量的空間分布亦表現出一定的相似性和規律性,其最大值均位于大汶河入湖口,而最小值均位于出湖口附近。N2O排放通量在大汶河入湖口和湖區東部較大,自東向西顯著減小,且出湖口水域N2O排放通量亦明顯減小。CO2排放通量表現為自大汶河入湖口向湖心至出湖口逐漸降低。

2.3 環境因子對沉積物酶活性的影響

堿性磷酸酶與脲酶和蔗糖酶活性均極顯著正相關,但脲酶和蔗糖酶活性無顯著相關性,見表3。3種酶活性均與沉積物總有機碳極顯著正相關,與水深、表層水體pH值和DO負相關,與沉積物TP正相關。堿性磷酸酶與水深顯著負相關;蔗糖酶與表層水體pH值、DO顯著負相關;脲酶與沉積物TP顯著正相關,與沉積物氧化還原電位顯著負相關。

采用多元線性回歸模型進一步分析影響沉積物酶活性的環境因子,結果見表4。脲酶活性主要受沉積物總有機碳TOC影響;蔗糖酶活性受沉積物總有機碳TOC、水體pH值pHw和水深H影響,其影響程度表現為pHw>H>TOC;堿性磷酸酶活性受沉積物總有機碳TOC、水體透明度SD、水體pH值pHw和NH3-N(公式中用ρ(NH3-N)表示)的影響,且其影響程度表現為SD>TOC>pHw>ρ(NH3-N)。各回歸方程可分別解釋脲酶、蔗糖酶和堿性磷酸酶活性變異程度的43.6%、73.4%和88.8%。

2.4 酶活性和環境因子對溫室氣體排放通量的影響

由表5可知:東平湖N2O排放通量與表層水體pH值極顯著負相關,與DO顯著負相關,與表層水體NO-2-N質量濃度顯著正相關;CO2排放通量與水深顯著正相關。N2O、CO2排放通量與沉積物理化性質和酶活性均無顯著相關性。

3 討 論

3.1 東平湖沉積物酶活性影響因素

湖泊沉積物酶活性的空間異質性與湖泊富營養化水平、沉積物的污染程度和沉積環境有關,富營養化水平高、沉積物污染嚴重、水動力條件差、水體交換能力弱的點位酶活性高[25]。由此分析,東平湖沉積物酶活性在湖區東部較高與該水域受到湖岸(老湖碼頭)人類活動的影響有關,外源性污染物被沉積物吸附[26],酶活性增強;脲酶和堿性磷酸酶活性在出湖口較高與出湖口水域狹窄且受水閘控制水體流動緩慢、水體自凈能力弱有關;大汶河入湖口水動力作用強,湖心區遠離湖岸,受湖岸帶人類活動影響小,水質較好,因此沉積物酶活性較低。

脲酶、蔗糖酶和堿性磷酸酶均屬于水解酶類,它們通過水解機制將大分子有機物轉化為寡聚或單分子物質,寡聚和單分子物質再被微生物同化,因此沉積物酶活性常與有機物含量具有顯著正相關性[27]。本研究中,3種酶活性均與沉積物總有機碳具有極顯著正相關關系,可表征沉積物中有機物的累積狀況。堿性磷酸酶能將有機磷化合物催化分解為無機磷化合物,并為微生物提供維持其生物活性的磷,從而影響脲酶和蔗糖酶活性[25,27],本研究中沉積物堿性磷酸酶活性與脲酶、蔗糖酶活性均極顯著正相關,與黎楊等[27]的研究結果一致。蔗糖酶可以表征碳循環,其活性反映碳代謝強度與規律,而脲酶是氮素轉化的關鍵酶,其活性反映沉積物有機氮向有效氮的轉化能力。本研究中,脲酶活性除與沉積物總有機碳極顯著正相關外,還與沉積物TP顯著正相關,與沉積物氧化還原電位顯著負相關,表明脲酶可在一定程度上表征沉積物質量。同時,3種酶活性還與表層水體pH值、透明度、NH3-N相關,說明酶活性與湖泊水質有一定關系,可作為分析湖泊水質的指標。

3.2 東平湖水氣界面溫室氣體排放影響因素

與沉積物酶活性的空間分布特征不同,東平湖N2O和CO2排放通量在大汶河入湖口較大,而在出湖口較小。由于大汶河流域內社會經濟發展和人口增長速度快,因此大量污染物質(特別是營養物質)經大汶河匯入東平湖,對東平湖水質尤其是入湖口水質造成較大影響[28],大量營養物質的輸入促進了水體溫室氣體的產生和排放[29]。

水體N2O主要來源于硝化和反硝化微生物的脫氮反應過程,其產生取決于水體N的轉化速率及環境條件[30],而NO-2-N是硝化和反硝化過程的中間產物,所以東平湖表層水體N2O的排放通量與水體NO-2-N質量濃度具有顯著正相關關系。同時,東平湖表層水體DO含量較高,使NO-3-N和NO-2-N還原受阻,水體反硝化作用受到抑制[31]。水體硝化作用反應底物NH3-N質量濃度較低,N2O排放通量與NH3-N質量濃度正相關,但相關性不顯著,也說明水體中N2O除了來源于水體硝化作用外,還來源于沉積物的反硝化作用[32]。東平湖N2O排放通量與表層沉積物脲酶活性正相關,在一定程度上反映了表層沉積物中耦合硝化—反硝化作用的存在。

表層水體pH值直接影響水體CO2的質量濃度,并通過影響微生物活性間接影響N2O的產生和排放[1,23]。東平湖湖水偏堿性,游離的CO2會轉變為碳酸鹽,由此導致CO2排放通量降低;較高的pH值可抑制微生物的活性,導致有機質分解速率減小,N2O產生量因細菌所需底物供應不足而減少,因此CO2和N2O的排放通量均與表層水體pH值負相關。N2O和CO2均與水深正相關,而水深與表層水體pH值負相關,與透明度、葉綠素a、溶解性有機碳均顯著負相關,與NO-3-N極顯著正相關,因此水深可間接反映表層水體的理化性質對溫室氣體排放的影響。

東平湖水氣界面溫室氣體排放通量與沉積物酶活性無顯著相關性,可能原因是,沉積物酶活性對沉積物中溫室氣體產生和排放的影響大于其對表層水體中溫室氣體的影響,同時研究測定的酶種類較少,缺少對影響沉積物中氮素轉化的關鍵酶如硝酸還原酶、亞硝酸還原酶的分析,需在今后的研究中補充完善。另外,不同酶的活性具有較為顯著的時空變異性[2],本研究僅分析了沉積物酶活性的空間分布特征,今后還需要進一步對酶活性的季節性變化進行分析,并結合微生物的組成和溫度等因素深入研究酶活性變化的關鍵影響因子及其對溫室氣體產生和排放的影響。

4 結 語

東平湖沉積物酶活性空間變異大,均表現為湖區東部(老湖碼頭附近)酶活性較高,而大汶河入湖口和湖心區酶活性較低,反映了湖岸區人類活動和水動力條件對沉積物酶活性的影響;N2O和CO2排放通量表現為大汶河入湖口較大,而出湖口較小,反映了大汶河營養物質輸入對東平湖溫室氣體排放的影響。沉積物總有機碳是影響東平湖沉積物酶活性的主要因子之一,此外蔗糖酶活性還受表層水體pH值和水深的影響,堿性磷酸酶還受到透明度、水體pH值和NH3-N的影響;酶活性與影響因子的回歸方程可分別解釋脲酶、蔗糖酶和堿性磷酸酶活性變異程度的43.6%、73.4%和88.8%。東平湖N2O和CO2排放通量主要受表層水體理化性質的影響,N2O排放通量的主要影響因子為表層水體無機氮濃度、pH值和DO,而CO2排放通量主要受水深的影響。

參考文獻:

[1] 楊平,仝川.淡水水生生態系統溫室氣體排放的主要途徑及影響因素研究進展[J].生態學報,2015,35(20):6868-6880.

[2] 曲靜怡,楊統一,陳芳艷,等.古運河(鎮江段)沉積物酶活性及其與氮磷陰離子的相關性[J].環境工程學報,2014,8(3):1057-1061.

[3] HILL B H,ELONEN C M,JICHA T M,et al.Sediment Microbial Enzyme Activity as an Indicator of Nutrient Limitation in the Great Rivers of the Upper Mississippi River Basin[J].Biogeochemistry,2010,97(2-3):195-209.

[4] 甘茂林,張景平,江志堅,等.大亞灣沉積物中可溶性有機氮與蛋白酶和脲酶活性特征及其影響因素[J].海洋環境科學,2016,35(6):806-813.

[5] 趙振華,曹晶晶,劉月利.秦淮河沉積物中營養鹽及相關酶活性時空變異特征及其響應關系研究[J].南京大學學報(自然科學),2016,52(1):36-45.

[6] 李小飛,林賢彪,程呂,等.崇明東灘沉積物胞外酶活性與溫室氣體排放潛力及其影響因素[J].水土保持學報,2015,29(3):202-207.

[7] 東平湖志編纂委員會.東平湖志[M].濟南:山東大學出版社,1993:8-13.

[8] 胡尊芳,宋印勝,孫建峰,等.東平湖枯豐水期水質健康風險評價[J].水電能源科學,2016,34(9):31-34.

[9] 劉學利,姚昕,董杰,等.東平湖水環境研究進展與展望[J].中國農學通報,2016,32(11):82-87.

[10] 田忠志,邢友華,姜瑞雪,等.東平湖表層沉積物中磷的形態分布特征研究[J].長江流域資源與環境,2010,19(6):719-723.

[11] 張智博,劉濤,伍青山,等.東平湖濕地土壤有機碳分布特征及其影響因素[J].人民黃河,2019,41(7):92-96,147.

[12] 張菊,何振芳,董杰,等.東平湖表層沉積物重金屬的空間分布及污染評價[J].生態環境學報,2016,25(10):1699-1706.

[13] 陳麗,劉衍君,曹建榮,等.東平湖表層沉積物中多氯聯苯來源初探[J].人民黃河,2020,42(5):81-86.

[14] 張菊,鄧煥廣,陳詩越.東平湖湖心區沉積物組成及重金屬污染評價[J].人民黃河,2010,32(7):47-48.

[15] 鄧煥廣,劉濤,魯長娟,等.山東省東平湖菹草(Potamogeton crispus)腐爛過程中水體溫室氣體溶存濃度及擴散通量變化[J].湖泊科學,2020,32(5):1484-1495.

[16] DENG HG,ZHANG J,WU J J,et al.Biological Denitrification in a Macrophytic Lake:Implications for Macrophytes-Dominated Lake Management in the North of China[J].Environmental Science and Pollution Research,2020,27(34):1-12.

[17] 范成新.湖泊沉積物調查規范[M].北京:科學出版社,2018:16-20.

[18] 鄧煥廣,張智博,張菊,等.東魚河春季沉積物反硝化脫氮作用與N2O排放研究[J].水土保持學報,2019,33(1):283-287.

[19] 關松蔭.土壤酶及其研究法[M].北京:農業出版社,1986:274-312.

[20] 國家環境保護總局《水和廢水監測分析方法》編委會.水和廢水監測分析方法[M].4版.北京:中國環境科學出版社,2002:266-280.

[21] 魯如坤.土壤農業化學分析方法[M].北京:中國農業科技出版社,2000:56-57.

[22] 鄧煥廣,張智博,劉濤,等.城市湖泊不同水生植被區水體溫室氣體溶存濃度及其影響因素[J].湖泊科學,2019,31(4):1055-1063.

[23] WILDING L P.Spatial Variability:Its Documentation,Accommodation and Implication to Soil Surveys[C]//NIELSEN D R,BOUMA J.Soil Spatial Variability.Wageningen:Proceedings of a Workshop of the ISSS and the SSSA,1985:166-194.

[24] 田翠翠,吳幸強,馮閃閃,等.東平湖沉水植物分布格局及其與環境因子的關系[J].環境科學與技術,2018,41(11):15-20.

[25] 袁軼君,畢永紅,朱孔賢,等.三峽水庫沉積物中堿性磷酸酶的活性[J].環境科學與技術,2014,37(1):60-64.

[26] 鄧煥廣,張菊,張智博,等.河岸滲濾系統對城市降雨徑流中氮磷去除效率的中試研究[J].聊城大學學報(自然科學版),2019,32(3):89-98.

[27] 黎楊,夏品華,葛皓,等.高原湖泊濕地草海沉積物胞外酶與有機碳的時空動態[J].中國環境科學,2017,37(7):2723-2730.

[28] 師吉華,李秀啟,董貫倉,等.東平湖水環境質量綜合評價[J].大連海洋大學學報,2015,30(4):391-397.

[29] 鄧煥廣,張智博,劉濤,等.城市河岸帶綠地N2O和CH4排放對硝態氮輸入的短期響應研究[J].聊城大學學報(自然科學版),2020,33(2):97-104.

[30] 邢濤,李俊雄,李彬彬,等.水生植物對草型富營養化湖泊氣態氮排放及沉積物氮去除的影響[J].生態學雜志,2018,37(3):771-778.

[31] DENG H G,WANG D Q,CHEN Z L,et al.Vertical Dissolved Inorganic Nitrogen Fluxes in Marsh and Mudflat Areas of the Yangtze Estuary[J].Journal of Environmental Quality,2014,43(2):745-752.

[32] 徐會顯,姜星宇,姚曉龍,等.鄱陽湖水體氧化亞氮排放特征及影響因素[J].湖泊科學,2016,28(5):972-981.

【責任編輯 呂艷梅】

猜你喜歡
影響因素
房地產經濟波動的影響因素及對策
零售銀行如何贏得客戶忠誠度
醫保政策對醫療服務價格影響因素的探討
東林煤礦保護層開采瓦斯抽采影響因素分析
影響農村婦女政治參與的因素分析
高新技術企業創新績效影響因素的探索與研究
水驅油效率影響因素研究進展
突發事件下應急物資保障能力影響因素研究
中國市場(2016年36期)2016-10-19 03:54:01
環衛工人生存狀況的調查分析
中國市場(2016年35期)2016-10-19 02:30:10
農業生產性服務業需求影響因素分析
商(2016年27期)2016-10-17 07:09:07
主站蜘蛛池模板: 国产在线98福利播放视频免费| 久草性视频| 亚洲无码高清一区| 国产91在线|日本| 人与鲁专区| 波多野结衣久久精品| 91麻豆国产精品91久久久| 亚洲日韩第九十九页| 草草线在成年免费视频2| 97成人在线视频| www精品久久| 成人国产免费| 丝袜美女被出水视频一区| 亚洲无限乱码| 欧美精品H在线播放| 91久久国产成人免费观看| 五月综合色婷婷| 99中文字幕亚洲一区二区| 欧洲av毛片| 国产精品亚洲专区一区| 免费xxxxx在线观看网站| 国产午夜精品一区二区三区软件| 国产欧美日韩一区二区视频在线| 中国精品久久| 无码福利视频| 欧美 亚洲 日韩 国产| 亚洲精品天堂在线观看| 狠狠五月天中文字幕| 欧美成a人片在线观看| 成人午夜视频网站| 亚洲视频影院| 亚洲伊人久久精品影院| 日韩精品毛片人妻AV不卡| 国产日本一线在线观看免费| 成AV人片一区二区三区久久| 久久精品国产电影| A级毛片无码久久精品免费| 亚洲欧州色色免费AV| 国产精品一区二区在线播放| 国产精品亚欧美一区二区三区| 日本午夜精品一本在线观看| 亚洲精品麻豆| 欧美日韩导航| 26uuu国产精品视频| 伊人久综合| 国产乱人免费视频| 一区二区三区在线不卡免费| 伊人久久久久久久久久| 91亚洲国产视频| 欧美啪啪精品| 欧美精品1区2区| 永久免费无码日韩视频| 毛片网站在线看| 91小视频在线观看| 欧美久久网| 毛片在线播放a| 日本成人在线不卡视频| 国产亚洲精品自在线| 99视频精品在线观看| 欧美午夜精品| 亚洲成人精品| 国产福利免费视频| 欧美成人午夜视频| 欧美成人二区| 久久国产精品影院| 国产资源站| 成人国产精品一级毛片天堂| 国产精品永久不卡免费视频| 国产成年无码AⅤ片在线| 久久9966精品国产免费| 扒开粉嫩的小缝隙喷白浆视频| 亚洲成人精品久久| 亚洲色婷婷一区二区| 国产va在线| 亚洲三级视频在线观看| 亚洲精品日产AⅤ| 777午夜精品电影免费看| 亚洲中文字幕久久精品无码一区| 国产成人福利在线视老湿机| 久久视精品| 日韩午夜福利在线观看| 成人韩免费网站|