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

    水沙環(huán)境變化對季節(jié)性多沙河流沉積物菌群特征的影響

    2017-11-23 02:50:49周孝德
    中國環(huán)境科學(xué) 2017年11期
    關(guān)鍵詞:豐水期渭河樣點

    劉 睿,周孝德

    ?

    水沙環(huán)境變化對季節(jié)性多沙河流沉積物菌群特征的影響

    劉 睿,周孝德*

    (西安理工大學(xué)西北旱區(qū)生態(tài)水利工程國家重點實驗室,陜西西安 710048)

    為辨析沉積物粒度特征、氮磷含量及上覆水水質(zhì)3類環(huán)境因子,對季節(jié)性多泥沙河流沉積物細菌群落多樣性和結(jié)構(gòu)特征變化的影響,采用T-RFLP技術(shù)分析了不同水文時期渭河(陜西段)流域沉積物細菌群落結(jié)構(gòu)變化特征,并通過冗余分析和蒙特卡羅檢驗識別了不同時期的關(guān)鍵環(huán)境驅(qū)動因子.結(jié)果表明,不同水文時期渭河水沙環(huán)境差異顯著.豐水期沉積物TN含量升高了1~2個數(shù)量級,且空間差異明顯減小(RSD=17%).沉積物中的氮磷比受上覆水氮磷比影響較大,兩者之間在豐水期呈現(xiàn)出一定的正相關(guān)關(guān)系(0.687,P=0.6).渭河干流沉積物細菌群落Shannon多樣性指數(shù)在2.02~3.18之間,不同時期的沿程變化特征各異.相對平水期和枯水期,豐水期渭河沉積物細菌多樣性的空間差異最大(RSD=37.4%),不同點位優(yōu)勢群落結(jié)構(gòu)組成的Bray-Curtis相似度最低(17.4%).PAT比對表明優(yōu)勢細菌大多歸屬于變形菌門(Proteobacteria)和擬桿菌門(Bacteroidetes).沉積物氮磷含量對細菌群落變化的影響,在平水期和枯水期較為直接.上覆水水質(zhì)對不同水文時期沉積物菌群變化解釋度的貢獻占比超過60%,對豐水期沉積物微生物環(huán)境的影響尤其顯著.

    渭河;沉積物;細菌群落;T-RFLP

    沉積物是河流水生態(tài)環(huán)境的重要組成部分.其間以細菌為主的微生物,是河流底部生物地球化學(xué)過程的重要參與者[1-3].在不同的水力條件和理化狀態(tài)下,河流底部沉積物與上覆水進行著不同程度的化學(xué)物質(zhì)[4-5]與生物信息[6]交換,從而持續(xù)影響和塑造著河床沉積物的微生物生態(tài)環(huán)境.

    對于季節(jié)性多泥沙河流,不同水文時期的河道水力特征、泥沙行為與顆粒特性迥異,影響沉積物微生物群落的環(huán)境因子更為復(fù)雜.基于分子微生物技術(shù),國內(nèi)外學(xué)者對全球不同區(qū)域的地表水體沉積物中,細菌群落結(jié)構(gòu)的垂向和水平分布及季節(jié)變化特征進行了大量研究[7-9],并分析了沉積物化學(xué)特性,如有機物[10]、氮[11]、磷[12-13]、重金屬[14-15]、pH值[8]、土質(zhì)類型[16]等對細菌群落的影響.但相關(guān)研究對水文及泥沙含量時空差異顯著的季節(jié)性多沙河流關(guān)注較少,此類河流的沉積物化學(xué)性質(zhì)、粒度特征及上覆水水質(zhì)同沉積物細菌群落特征的相關(guān)關(guān)系,以及上述各類環(huán)境因子對細菌群落變化影響程度的相對大小尚未明晰.

    本文以我國西北地區(qū)典型季節(jié)性多泥沙河流——渭河為例,采用末端限制性片段長度多態(tài)性技術(shù)(T-RFLP),分析了渭河干流及南北兩典型支流的沉積物細菌群落結(jié)構(gòu)及多樣性的時空變化特征,結(jié)合限制性排序分析,探究了不同水文時期影響群落的關(guān)鍵環(huán)境因子的差異,比較了沉積物氮磷、沉積物粒度特征、上覆水水質(zhì)三類環(huán)境因子對細菌群落變化的貢獻量,旨在為同類河流水生態(tài)環(huán)境相關(guān)研究和管理保護提供參考和依據(jù).

    1 材料與方法

    1.1 研究區(qū)域與采樣方法

    陜西渭河流域地處半干旱半濕潤地區(qū),冬季干燥寒冷,夏季炎熱多雨,不同季節(jié)水沙變化懸殊,年徑流量的56%~70%和輸沙量的60%~ 90%集中在汛期7~10月[17-18].渭河南北支流的水沙特征迥異,其南岸秦嶺北麓是黃河中下游重要的洪水來源區(qū),北岸黃土高原則為黃河流域的主要產(chǎn)沙區(qū)[19].考慮到渭河陜西段沿程環(huán)境變化及兩岸支流水沙特征差異,本研究選取的渭河(陜西段)流域,涵蓋了渭河中游(寶雞峽至咸陽鐵橋)和下游(咸陽鐵橋至潼關(guān)入黃口)干流,以及黑河(南)與涇河(北)兩條重要支流.渭河干流采樣點W1~W3屬中游河段,分別處于寶雞、咸陽武功、咸陽興平行政區(qū);W4~W6屬下游河段,其中W4、W5處于西安市、W6處于潼關(guān)行政區(qū).W2與W4分別位于黑河與涇河匯入之前的干流河段,W3與W5分別位于黑河與涇河入?yún)R之后的干流河段.具體采樣點位置如圖1所示.

    圖1 采樣點地理位置

    河流表層沉積物樣本(約0~10cm)使用彼得遜抓泥斗于河流斷面中心10m范圍內(nèi)采集,每個斷面采集3~5次,分裝于500mL無菌廣口聚乙烯瓶.采集沉積物之前,用有機玻璃采水器在沉積物采集點上方采集水樣3L用于水質(zhì)檢測.沉積物和水樣經(jīng)冰上運輸至實驗室4℃保存,并盡快進行后續(xù)分析.

    1.2 環(huán)境理化指標分析方法

    充分冷凍干燥的沉積物樣本,經(jīng)混勻、研磨、過篩后,采用Foss8400全自動凱氏定氮儀測定總氮(TN)含量,采用瑞典FI515流動注射分析儀測定氨氮(NH4+-N)和硝態(tài)氮(NO3--N)含量,采用HClO4-H2SO4消煮、鉬銻抗比色法測定全磷(TP)含量,并通過Mastersizer 2000激光粒度儀測定凍干沉積物的顆粒粒度分布.水質(zhì)相關(guān)指標及分析方法與之前對渭河浮游菌群的研究相同[20].

    1.3 T-RFLP實驗方法

    沉積物樣本中的總DNA使用PowerSoil ? DNA Isolation Kit(Mo Bio)提取.選擇5'端帶FAM熒光標記的細菌通用引物27F(5'-GAGTTT- GATCMTGGCTCAG-3')/1492R(5'-GGTTACC- TTGTTACGACTT-3')進行PCR擴增.擴增體系和反應(yīng)條件采用之前土壤和地表水相似研究中的方法[21-22].擴增產(chǎn)物經(jīng)純化(TIANquick Mini Purification Kit)后,采用限制性內(nèi)切酶HhaⅠ(TaKaRa)酶切2h,70℃滅活15min送檢.

    1.4 數(shù)據(jù)分析方法

    將T-RFLP峰譜圖掃圖結(jié)果中50~500bp的片段作為有效片段[22].使用PRIMER 5.0統(tǒng)計各樣本中細菌群落的物種豐度(),均勻度()和Shannon多樣性(),并計算樣本群落間Bray- Curtis相似度,完成群落聚類分析.通過在線比對工具Phylogenetic Assignment Tool(PAT) (http: //mica.ibest.uidaho.edu/pat.php),識別不同片段代表的細菌類群[23].采用Canoco 4.5對細菌群落限制性片段(T-RF)相對豐度>0.05的部分進行去趨勢分析(DCA)[24],結(jié)果表明各季節(jié)優(yōu)勢細菌群落數(shù)據(jù)均適合冗余分析(RDA).最后通過手動選擇,在候選指標中依次選出不存在共線性的環(huán)境因子[25],并對各個入選指標和最終排序結(jié)果進行蒙特卡羅檢驗.

    2. 結(jié)果與討論

    2.1 沉積物粒度特征、氮磷含量及上覆水水質(zhì)

    渭河各樣點沉積物顆粒粒度特征指標如表1所示,其中粒度多樣性(PSD)通過Shannon多樣性指數(shù)公式計算得出[26].西安樣點(W4)沉積物同其它樣點相比,中位粒徑明顯偏大.這同宋進喜等[27]2011年在相近河段的采樣分析結(jié)果相一致.除W4樣點外,非汛期渭河流域沉積物中值粒徑(D50)均值為33.1μm,沉積物構(gòu)成以中粉砂和細粉砂為主[28].由于降水較少,枯水期采集到的沉積物粒度構(gòu)成通??梢苑€(wěn)定較長時間(6~7個月)[29].對比豐枯兩季各樣點的粒度分布差異,可以看出豐水期雨洪對渭河陜西段沉積物的影響存在明顯的分段特征.結(jié)合豐水期W1樣點上覆水中極高的顆粒物濃度(TSS=16.33g/L)可以看出,上游雨洪攜帶著大量泥沙.這些上游來沙使寶雞峽下游W1樣點沉積物中的細顆粒和粒度多樣性顯著增加,體積平均粒徑顯著減小.進入中游(寶雞峽-咸陽鐵橋)寬淺河段后,河流中的泥沙逐步沉積(同期上覆水TSS經(jīng)此河段降低83%),使相應(yīng)樣點(W2和W3)沉積物的體積平均粒徑和中位粒徑明顯增大.進入渭河下游蜿蜒河段后,顆粒進一步淤積,相應(yīng)樣點(W5和W6)沉積物的比表面積減小,顆粒平均粒徑變化程度降低.

    豐水期渭河流域沉積物TN含量較平、枯兩季提升了1~2個數(shù)量級,氮磷比(N/P)也顯著增加(表2).同時,流域內(nèi)各樣點沉積物TN與N/P的空間相對標準偏差,在平、枯兩季均高于88%,而在豐水期顯著減小至17%和48%.結(jié)合對沉積物粒度特征沿程變化的分析可以看出,豐水期上游泥沙在渭河陜西段的沿程沉積過程,顯著改變了沉積物的化學(xué)性質(zhì),使流域沉積物趨于均質(zhì)化.渭河干流與涇河樣點沉積物TP含量的最低值出現(xiàn)在枯水期,且此時的空間異質(zhì)性(RSD=65%)較其它時期更高.上覆水的實測結(jié)果顯示,枯水期各樣點水體TP含量與平水期相近,而溶解態(tài)總磷(DTP)較平水期平均增加了一倍.冬季水生動植物殘體礦化可能引發(fā)沉積物氧化還原條件改變[30],從而導(dǎo)致沉積物發(fā)生磷解吸,使含磷量降低.

    表1 豐枯兩季各樣點沉積物顆粒粒度特征

    Pearson相關(guān)性分析結(jié)果表明,豐水期沉積物粒度特征與TP含量存在顯著相關(guān)關(guān)系.顆粒比表面積()與TP含量顯著負相關(guān)(-0.844,= 0.008),中值粒徑(D50)與TP含量顯著正相關(guān)(0.867,=0.005).沉積物氮磷比與上覆水氮磷比之間,僅在豐水期存在中度正相關(guān)關(guān)系(相關(guān)系數(shù)為0.687),但顯著性水平不高(=0.6).

    表2 沉積物氮磷及上覆水水質(zhì)特征

    注:Temp-溫度,SpCond-電導(dǎo)率,Sal-鹽度,TrubSC-濁度,PCY-藻細胞,RSD-相對標準偏差.

    2.2 沉積物細菌群落多樣性的時空分布

    三季24個沉積物樣點共檢出限制性片段(T-RF)97條.渭河中下游流域沉積物細菌物種豐度的空間差異大于群落均勻度和多樣性,相對標準偏差(RSD)呈現(xiàn)豐水期>枯水期>平水期的關(guān)系(表3).枯水期渭河干流物種豐度(21)和多樣性(2.73)最高,豐水期最低.這與東平湖沉積物細菌多樣性特征的季節(jié)變化規(guī)律一致[31].平水期支流與干流的物種數(shù)、均勻度和多樣性最為接近;枯水期和豐水期的物種數(shù)和多樣性均表現(xiàn)為黑河高于干流,涇河低于干流的關(guān)系.整體來看,黑河樣點沉積物細菌多樣性較高,穩(wěn)定在2.8以上,而涇河樣點沉積物細菌多樣性季節(jié)波動劇烈,不同時期的多樣性均值為1.75.作為渭河南、北兩岸支流的代表,黑河流域土地以林地和農(nóng)田為主,河流含沙量小、有機污染少、水質(zhì)狀況較好[32-33];而涇河流經(jīng)黃土高原地區(qū),含沙量大,并接受大量的城市污水和工業(yè)廢水,水質(zhì)較差[34-35].黑河與涇河沉積物微生物多樣性的差異,與水文特征與流域污染狀況相對應(yīng),體現(xiàn)了渭河南北兩岸支流迥異的水生態(tài)環(huán)境.

    不同時期渭河干流和的沿程變化趨勢相似,物種數(shù)較高的點位,多樣性也較高.平水期渭河干流沉積物細菌Shannon多樣性指數(shù)變化范圍在2.37~2.82之間(圖2).中游河段(W1~W3)多樣性穩(wěn)定而略有增加,經(jīng)西安臨潼河段(W4、W5)多樣性下降到最低值后,至下游潼關(guān)W6樣點恢復(fù)至中游水平.枯水期渭河干流沉積物細菌多樣性自上游至下游逐步減小(2.12~3.18),豐水期則沿程表現(xiàn)出較大的上下波動(2.02~2.89).渭河干流沉積物氮、磷含量在豐水期趨于空間均質(zhì)化,但細菌群落多樣性卻表現(xiàn)出高于其它時期的空間差異(RSD=37.4%).這些均表明了水沙特征變化劇烈的季節(jié)性河流中,沉積物細菌群落變化驅(qū)動因子的復(fù)雜性.

    圖2 不同水文時期沉積物細菌群落Shannon多樣性指數(shù)沿程變化趨勢

    Fig.2 Shannon index variation of sediment bacterial communities along Weihe River

    表3 河流域沉積物細菌群落多樣性統(tǒng)計

    2.3 優(yōu)勢細菌群落構(gòu)成及聚類分析

    將各樣點中相對豐度高于5%的所有T-RFs集合視為優(yōu)勢細菌群落,以不同顏色表示優(yōu)勢T-RF在樣點中的相對豐度,并計算各樣點Bray-Curtis相似度,對細菌群落進行聚類分析,得到的結(jié)果如圖3所示.平水期、枯水期和豐水期流域內(nèi)各樣點檢測出的優(yōu)勢T-RF種類分別為26、21和31種.豐水期沉積物中的優(yōu)勢細菌種類最豐富,細菌群落組成的空間相似度最低(17.4%);枯水期沉積物優(yōu)勢菌種最少,菌群結(jié)構(gòu)的空間相似度最高(35.8%).

    優(yōu)勢片段202bp和204bp出現(xiàn)在枯水期除涇河(J)外的所有樣點中,相對豐度分別達到6.0%~28.1%和4.8%~14.0%,是枯水期渭河中下游沉積物普遍存在的優(yōu)勢細菌類群.經(jīng)PAT比對,202bp和204bp T-RF通常來自β-變形菌綱中的菌群,前者多歸屬于纖毛菌屬(sp.),后者可能涉及叢毛單胞菌屬(sp.)、食酸菌屬(sp.)、代爾夫特菌屬(sp.)等.87bp和99bp出現(xiàn)在枯水期超過50%的樣本中,其中99bp在下游樣點(W5和W6)中的豐度達到26%以上,是枯水期渭河下游河段的關(guān)鍵優(yōu)勢菌種.經(jīng)比對,87bp T-RF指示的菌群極可能為擬桿菌門黃桿菌屬中的噬冷黃桿菌().99bp T-RF則可能來自ε-變形菌綱(-)中的彎曲菌屬(sp.)、沃林氏菌屬(sp.)、螺桿菌屬(sp.),γ-變形菌綱(Proteobacteria)中的硫微螺菌屬(sp.)、假單胞菌屬(sp.)、曼氏桿菌屬(sp.),以及擬桿菌門中的擬桿菌屬(sp.)、列文虎克菌屬(sp.)和黃桿菌屬(sp.)等.

    57、58、63、77、93bp是不同季節(jié)都有檢出的優(yōu)勢片段.其中57bp和63bp經(jīng)常出現(xiàn)在渭河中游(W1、W2)與黑河(H)樣點,93bp則檢出于不同季節(jié)的兩支流和咸陽(W3)、西安(W4)樣點.在比對結(jié)果中,58bp T-RF指示的菌群常見于土壤、沉積物及白蟻和蚯蚓消化道,從屬于δ-變形菌綱(δ)和α-變形菌綱(α).77bp T-RF指示的菌群多鑒定為放線菌門下的節(jié)細菌屬(sp.)和微顫菌屬(sp.).93bp T-RF廣泛分布在深海、湖泊、污水處理系統(tǒng)等不同類型的環(huán)境中,菌群主要涉及ε-變形菌綱中的螺桿菌屬(sp.),δ-變形菌綱中的暗桿菌屬(sp.)、地桿菌屬(sp.),α-變形菌綱的鞘脂單胞菌屬(sp.)和sp.菌屬,以及擬桿菌門中的黃桿菌屬、地桿菌屬和噬纖維菌屬(sp.)等.

    圖3 沉積物優(yōu)勢細菌群落T-RF構(gòu)成及聚類關(guān)系(N-平水期,D-枯水期,W-豐水期)

    51bp和52bp幾乎僅在涇河樣點中出現(xiàn),且在涇河沉積物細菌群落中的相對豐度均超過70%,是該樣點的絕對優(yōu)勢菌群.它們曾經(jīng)被檢出于厭氧膜生物反應(yīng)器[36]和濕地香蒲根際環(huán)境[37]中,可能指示了對污染物降解具有積極意義的菌群.

    統(tǒng)計不同樣點和不同時期的群落相似度均值(表4),結(jié)果表明W1、W3和H樣點在不同季節(jié)的Bray-Curtis相似度最高(均>35%),沉積物優(yōu)勢細菌群落結(jié)構(gòu)最為穩(wěn)定.涇河樣點相似度均值最低(僅5.8%),沉積物優(yōu)勢細菌群落結(jié)構(gòu)的季節(jié)差異最大.不同時期流域內(nèi)各樣點優(yōu)勢菌群相似度均值的大小關(guān)系為:枯水期>平水期>豐水期,表明豐水期水沙過程增大了沉積物優(yōu)勢群落結(jié)構(gòu)在流域內(nèi)的空間差異,對河流沉積物微生物生態(tài)過程具有重要影響.

    表4 Bray-Curtis相似度均值(%)

    2.4 不同水文時期關(guān)鍵環(huán)境因子的識別

    在沉積物氮磷指標(TNs、TPs、N/Ps)的基礎(chǔ)上,依次加入沉積物顆粒粒度特征指標(,[4,3],D[3,2],d10,d50,d90,PSD)和上覆水水質(zhì)指標(TSS、IMn、CODCr、TN、TP、TOC、DIMn、DCODCr、NO2-N、NO3--N、NH4+-N、DTN、DTP、C/N、C/P、N/P、溶解性碳氮比DC/N、溶解性碳磷比DC/P、溶解性氮磷比DN/P、Temp、pH值、SpCond、Sal、DO、CHL、PCY)作為解釋變量,參與CCA或RDA分析,得出各類指標對不同水文時期渭河流域沉積物優(yōu)勢細菌群落構(gòu)成(優(yōu)勢T-RFs)及群落多樣性指數(shù)(、、)變化的解釋程度(表5).從結(jié)果中可以看出,僅通過沉積物氮磷指標,對枯水期優(yōu)勢菌群結(jié)構(gòu)變化的解釋度可以達到74.3%;但對平水期和豐水期優(yōu)勢菌群變化的解釋度在50%以下,且分析結(jié)果未能通過蒙特卡羅檢驗(>0.05).由于枯水期渭河流量與流速較低,河流上覆水與沉積物的物質(zhì)能量交互減弱,沉積物微生物群落變化主要受沉積物自身性質(zhì)影響.而在平水期和豐水期,渭河水量和含沙量增大,河道內(nèi)物質(zhì)遷移轉(zhuǎn)化行為多樣而活躍,影響沉積物細菌群落變化的環(huán)境因素更加復(fù)雜.考慮沉積物顆粒特性的影響后,環(huán)境因子對沉積物優(yōu)勢細菌群落變化的解釋度均超過50%;進一步加入上覆水指標后,解釋度可以達到86%以上,同時,排序結(jié)果的顯著性水平也明顯提高.這一結(jié)果表明,對于季節(jié)性多泥沙河流,除過沉積物氮磷等化學(xué)特性之外,顆粒粒度與上覆水水質(zhì)也在很大程度上影響著沉積物的細菌群落變化.

    將標準化后的沉積物細菌群落多樣性指標與優(yōu)勢T-RF相對豐度數(shù)據(jù)放入同一個排序空間,從全部環(huán)境指標中手動選擇共同解釋度最高,且不存在共線性關(guān)系的最優(yōu)解釋因子組合,結(jié)果如圖4所示.平水期沉積物總磷(TPs)、氮磷比(N/P s)與上覆水COD、溶解態(tài)硝酸鹽氮(NO3--N)、溫度(Temp)和pH值對沉積物群落特征變化的共同解釋度為96.2% (=0.018).pH值和N/P s對物種變化的單獨解釋度最高,分別為27%(=0.032)和22.5%(=0.030).且pH在平水期與沉積物細菌群落均勻度顯著正相關(guān)(0.775,=0.024).氮磷作為構(gòu)成生物體的重要生源物質(zhì),在不同類型地表水的沉積物中,往往是與細菌群落分布和多樣性高低相關(guān)性最強的環(huán)境因子[11-13].而對夏季太湖不同植物類型湖區(qū)的研究發(fā)現(xiàn)[8],pH值是除有機物含量外,對沉積物細菌群落變化影響最顯著的指標.酸堿度一方面可能反映離子和金屬等其它環(huán)境因子的變化狀態(tài),另一方面也可以直接對細菌的新陳代謝過程產(chǎn)生影響[38-39].此外,從各樣點在RDA排序圖上的分布和距離可以看出,平水期渭河中游(W1、W2、W3)沉積物優(yōu)勢細菌群落結(jié)構(gòu)和多樣性狀態(tài)較為相似,而同下游(W4、W5、W6)各樣點存在明顯差異.

    枯水期沉積物總氮(TNs)、總磷(TPs)、體積平均粒徑([4,3]),及上覆水溶解性硝酸鹽氮(NO3--N)、溶解性碳氮比(DC/N)和葉綠素a(Chl a)構(gòu)成的排序空間對沉積物優(yōu)勢細菌群落變化解釋度為97.0%(=0.01).其中,DC/N和Chl a的單獨解釋度分別為42.2%(=0.062)和29.0%(=0.008).DC/N與沉積物細菌群落物種豐度(-0.755,=0.03)均勻度(-0.853,=0.007)和多樣性指數(shù)(-0.827,=0.011)均呈顯著負相關(guān)關(guān)系.上覆水NO3--N與99bp的相對豐度顯著負相關(guān)(-0.960,=0.04),Chl a與58bp的相對豐度顯著負相關(guān)(-0.998,=0.039).渭河流域枯水期實測水體Chl a平均濃度僅為平水期的五分之一.冬季藻類的消亡可以使氮素重新沉降,刺激表層沉積物中的細菌生長,增加氮循環(huán)相關(guān)細菌群落的多樣性,但對具體菌種相對豐度的影響不一[40-41].沉積物體積平均粒徑對細菌群落變化的解釋程度,超過TNs和TPs排在環(huán)境因子中的第三位,表明枯水期沉積物粒度特征對渭河流域沉積物生態(tài)環(huán)境的重要性.在對海洋環(huán)境的研究中,也注意到了粒度特征帶來的影響,認為沙質(zhì)和泥質(zhì)沉積物上的物理化學(xué)過程不盡相同[14,42],而顆粒的大小決定著細菌附著程度和攜氧能力[16],從而導(dǎo)致了微生物群落的差異.排序空間中,渭河干流及黑河樣點主要沿第二排序軸表現(xiàn)出梯度變化,且?guī)缀醢囱爻添樞蜃韵露吓挪?表現(xiàn)出枯水期河流沉積物理化及微生物生態(tài)變化的空間連續(xù)性.

    表5 不同環(huán)境變量對優(yōu)勢群落構(gòu)成和群落多樣性指數(shù)的解釋度

    注:P1為第一排序軸蒙特卡羅檢驗結(jié)果,P2為整個排序空間的檢驗結(jié)果.通過手動選擇獲得共同解釋度最高,且不存在共線性關(guān)系的最優(yōu)解釋因子組合,入選指標依解釋度由高到低排序.

    豐水期入選的環(huán)境因子全部為上覆水水質(zhì)指標,未過濾水樣碳磷比(C/P)、COD、藻細胞數(shù)(PCY)、溶解性氨氮(NH4+-N)和溶解性碳氮比(DC/N)共同解釋了94.6% (=0.064)的沉積物細菌群落變化.其中,DC/N和NH4+-N是單獨解釋度最高的指標,分別解釋了38.0%(=0.132)和18.1%(=0.180)的沉積物物種變化,但未能達到統(tǒng)計學(xué)意義上的顯著性.上覆水溶解性碳氮比與沉積物細菌群落均勻度和多樣性指數(shù)顯著負相關(guān),相關(guān)系數(shù)分別為-0.881(=0.009)和-0.854 (=0.015).對渭河不同季節(jié)的實測研究發(fā)現(xiàn)[43],夏季沉積物和間隙水中的氮磷含量,同河流補給地下水的交換通量正相關(guān),即河流上覆水中的氮磷在這一時期向沉積物發(fā)生了遷移.豐水期渭河河道內(nèi)垂向的物質(zhì)遷移過程(河流補給地下水和泥沙沉積),增強了水質(zhì)對沉積物細菌群落的影響.與枯水期相似,豐水期除涇河樣點外的其它各樣點依然沿第二軸分布,但分布順序與空間位置并不存在對應(yīng)關(guān)系.渭河中游流域樣點(W1、W2、W3、H)的環(huán)境及物種狀況相對接近,而西安段樣點(W4、W5、J)之間的環(huán)境和微生物狀態(tài)差異顯著.

    圖4 不同水文時期沉積物細菌群落空間變化與環(huán)境因子的RDA分析

    圖5 不同時期各類環(huán)境因子對解釋沉積物細菌群落變化的貢獻程度

    不同水文時期的冗余分析結(jié)果表明,渭河沉積物細菌群落的關(guān)鍵環(huán)境驅(qū)動因子,存在季節(jié)性變化又具有一定的共性(圖5).沉積物氮磷指標TPs、TNs或N/Ps,是平水期和枯水期沉積物菌群變化的關(guān)鍵解釋因子,在整體解釋變量中的貢獻度分別為37%和14%.沉積物粒度特征對細菌群落的影響僅在枯水期有直接的體現(xiàn).上覆水水質(zhì)對沉積物細菌群落變化的解釋貢獻度,在不同季節(jié)均超過60%.其中,溶解態(tài)無機氮(NO3--N和NH4+-N)對沉積物細菌群落的影響,在3個時期中均有不同形式的體現(xiàn).浮游植物相關(guān)的Chl a和PCY指標對優(yōu)勢菌群變化,在枯水期和豐水期表現(xiàn)出不可忽視的作用.而平水期和豐水期共有的解釋變量為上覆水COD.豐水期RDA空間中的環(huán)境因子全部為上覆水水質(zhì)指標,表明了這一時期渭河水體垂向交互作用的劇烈程度,以及上覆水環(huán)境對季節(jié)性多泥沙河流沉積物細菌群落變化的重要影響.總體來說,沉積物氮磷含量的影響比顆粒粒度特征更為顯著,上覆水溶解性無機氮含量與溶解性碳氮比對渭河流域沉積物優(yōu)勢細菌群落結(jié)構(gòu)變化具有重要意義.

    上述結(jié)果表明,對于季節(jié)性多泥沙河流,尤其在水沙量大、河道垂向物質(zhì)能量交互劇烈的水文時期,除過沉積物氮磷等化學(xué)特性之外,沉積物顆粒粒度與上覆水水質(zhì)對沉積物的細菌群落變化的影響作用也不可忽視.

    3 結(jié)論

    3.1 平枯兩季渭河流域沉積物以中粉砂和細粉砂為主,TP平均含量高于TN.豐水期渭河上游來水攜沙在中下游河段發(fā)生逐步沉積,使渭河(陜西段)沉積物粒度特征與氮磷含量顯著改變.沉積物TN含量在豐水期升高了1~2個數(shù)量級,且空間差異明顯減小(RSD=17%).較其它時期,豐水期沉積物氮、磷含量與上覆水氮、磷濃度之間的相關(guān)性系數(shù)最高(0.687,=0.6).

    3.2 渭河流域沉積物樣本中,基于T-RFLP技術(shù)共檢出16s rDNA限制性片段(T-RF)97條.干流沉積物細菌群落多樣性在在2.37~2.82之間,不同時期的沿程變化特征各異.渭河豐水期水沙過程的改變,使沉積物氮磷性質(zhì)趨于空間均質(zhì)化,卻增加了細菌群落特征的空間差異,Shannon多樣性空間相對標準偏差達到37.4%,優(yōu)勢菌群Bray- Curtis相似度為17.4%.表明了季節(jié)性多沙河流沉積物細菌群落變化的復(fù)雜性,及汛期水沙變化對河流沉積物微生物生態(tài)環(huán)境的重大影響.流域優(yōu)勢菌群以變形菌門和擬桿菌門為主.不同時期不同河段的沉積物優(yōu)勢菌群在綱水平上表現(xiàn)出較大的差異.

    3.3 在沉積物自身氮磷指標的基礎(chǔ)上考慮沉積物顆粒特性和上覆水水質(zhì)的影響后,對渭河流域沉積物細菌群落結(jié)構(gòu)和多樣性變化的解釋度均有增加,且排序結(jié)果的顯著性水平也明顯提高.超過94%的沉積物細菌群落結(jié)構(gòu)及多樣性變化,與沉積物及上覆水理化指標相關(guān),但不同水文時期菌群變化的關(guān)鍵驅(qū)動因子不同.沉積物氮磷及粒度特征對細菌群落變化的影響,在平水期和枯水期表現(xiàn)得更為直接;而上覆水水質(zhì)對不同時期沉積物菌群變化的解釋度貢獻量均超過60%,對豐水期沉積物細菌群落的影響尤其顯著.

    [1] Brune A, Frenzel P, Cypionka H. Life at the oxic-anoxic interface: microbial activities and adaptations [J]. Fems Microbiology Reviews, 2000,24:691–700.

    [2] Droppo I G, Krishnappan B G, Lawrence J R. Microbial interactions with naturally occurring hydrophobic sediments: Influence on sediment and associated contaminant mobility [J]. Water Research, 2016,(92):121-130.

    [3] 鮑林林,王曉燕,陳永娟,等.北運河沉積物中主要脫氮功能微生物的群落特征 [J]. 中國環(huán)境科學(xué), 2016,36(5):1520-1529.

    [4] Liu T Z, Wang H, Zhang Z, et al. Application of synthetic iron- oxide coated zeolite for the pollution control of river sediments [J]. Chemosphere, 2017,180:160-168.

    [5] 汪 淼,王圣瑞,焦立新,等.滇池沉積物內(nèi)源氮釋放風(fēng)險及控制分區(qū) [J]. 中國環(huán)境科學(xué), 2016,36(3):798-807.

    [6] Russo S A, Hunn J, Characklis GW. Considering Bacteria- Sediment Associations in Microbial Fate and Transport Modeling [J]. Journal of Environmental Engineering, 2011,137(8):697-706.

    [7] Xu P, Laura G. Leff. Longitudinal changes in the benthic bacterial community of the Mahoning River (Ohio, U.S.A.)[J]. Hydrobiologia, 2004,522:329–335.

    [8] Shao K Q, Gao G, Wang Y P, et al. Vertical diversity of sediment bacterial communities in two different trophic states of the eutrophic Lake Taihu, China [J]. Journal of Environmental Sciences, 2013,25(6):1186–1194.

    [9] 鄭艷玲,侯立軍,陸 敏,等.崇明東灘夏冬季表層沉積物細菌多樣性研究 [J]. 中國環(huán)境科學(xué), 2012,32(2):300-310.

    [10] 王 凱,鄒 立,高冬梅,等.黃河口潮灘春季細菌群落的分布特征及其影響因素研究 [J]. 中國海洋大學(xué)學(xué)報, 2016,46(1): 108-115.

    [11] 王永霞.云南高原程海湖沉積物中的細菌多樣性研究[D]. 昆明:云南大學(xué), 2014.

    [12] 宋洪寧.東平湖沉積物細菌多樣性分析[D]. 泰安:山東農(nóng)業(yè)大學(xué), 2010.

    [13] 魯克帥.南四湖沉積物細菌多樣性研究[D]. 曲阜:曲阜師范大學(xué), 2013.

    [14] Wang L P, Liu L S, Zheng B H, et al. Analysis of the bacterial community in the two typical intertidal sediments of Bohai Bay, China by pyrosequencing [J]. Marine Pollution Bulletin, 2013,72:181–187.

    [15] 劉愛菊,王洪海,潘嘉芬,等.孝婦河表層沉積物中重金屬賦存形態(tài)與微生物群落組成 [J]. 中國環(huán)境科學(xué), 2010,30(8):1103- 1109.

    [16] Verma P, Raghavan R V, Jeon C O, et al. Complex bacterial communities in the deep-sea sediments of the Bay of Bengal and volcanic Barren Island in the Andaman Sea [J]. Marine Genomics, 2017,31:33–41.

    [17] 孫 悅,李棟梁,朱擁軍.渭河徑流變化及其對氣候變化與人類活動的響應(yīng)研究進展 [J]. 干旱氣象, 2013,31(2):396-405.

    [18] 查小春,延軍平.全球變化下秦嶺南北河流徑流泥沙比較分析 [J]. 地理科學(xué), 2002,22(4):403-407.

    [19] 劉曉瓊,劉彥隨,李同昇.基于小波多尺度變換的渭河水沙演變規(guī)律研究 [J]. 地理科學(xué), 2015,35(2):211-217.

    [20] 劉 睿,吳 巍,周孝德,等.渭河浮游細菌群落結(jié)構(gòu)特征及其關(guān)鍵驅(qū)動因子 [J]. 環(huán)境科學(xué)學(xué)報, 2017,37(3):934-944.

    [21] Huang Y, Zou L, Zhang S Y, et al. Comparison of bacterioplankton communities in threeheavily polluted streams in China [J]. Biomedical and Environmental Sciences, 2011,24(2): 140-145.

    [22] Xie S G, Sun W M, Luo C L, et al. Stable Isotope Probing Identifies Novel m-Xylene Degraders in Soil Microcosms from Contaminated and Uncontaminated Sites [J]. Water Air and Soil Pollution, 2010,212(1-4):113-122.

    [23] Kent A D, Smith D J, Benson B J, et al. Web-Based Phylogenetic Assignment Tool for Analysis of Terminal Restriction Fragment Length Polymorphism Profiles of Microbial Communities [J]. Applied and Environmental Microbiology, 2003,69(11):6768-6776.

    [24] 王 萃,鮑林林,王曉燕.密云水庫底泥和庫濱區(qū)土壤中氨氧化細菌的多樣性和豐度 [J]. 環(huán)境科學(xué)學(xué)報, 2013,33(12):3334- 3340.

    [25] O'Brien R M. .Caution regarding rules of thumb for variance inflation factors [J]. Quality & Quantity, 2007,41(5):673-690.

    [26] Leduc D, Rowden A A, Probert P K, et al. Further evidence for the effect of particle-size diversity on deep-sea benthic biodiversity [J]. Deep Sea Research Part I: Oceanographic Research Papers, 2012,63:164-169.

    [27] 宋進喜,于 芳,王 珍.渭河陜西段河床沉積物的粒度參數(shù)分析 [J]. 南水北調(diào)與水利科技, 2013,11(4):75-78.

    [28] 李 勇,李海燕,趙應(yīng)權(quán).沉積物粒度特征及其對環(huán)境的指示意義——以濠河為例 [J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2015, 45(3):918-925.

    [29] 楊小剛,宋進喜,陳 佳,等.渭河陜西段潛流帶沉積物重金屬變化初步分析 [J]. 環(huán)境科學(xué)學(xué)報, 2014,34(8):2051-2061.

    [30] 張臺凡,宋進喜,楊小剛,等.渭河陜西段沉積物中總磷、總氮時空分布特征及其影響因素研究 [J]. 環(huán)境科學(xué)學(xué)報, 2015,35(5): 1393-1399.

    [31] 宋洪寧,杜秉海,張明巖,等.環(huán)境因素對東平湖沉積物細菌群落結(jié)構(gòu)的影響[J]. 微生物學(xué)報, 2010,50(8):1065-1071.

    [32] 劉鐵龍.陜西省渭河主要支流水量與調(diào)控能力狀況分析 [J]. 陜西水利, 2015,(4):15-19.

    [33] 牛珊珊.2005—2010年渭河主要支流水質(zhì)變化趨勢 [J]. 人民黃河, 2014,36(10):89-91.

    [34] 李家科,李懷恩,董 雯,等.渭河關(guān)中段典型支流非點源污染監(jiān)測與負荷估算 [J]. 環(huán)境科學(xué)學(xué)報, 2011,31(7):1470-1478.

    [35] 于 芳,宋進喜,殷旭旺,等.涇河陜西段水體污染特征分析 [J]. 北京師范大學(xué)學(xué)報(自然科學(xué)版), 2013,49(4):421-424.

    [36] 隋立新.不同溫度下微生物群落結(jié)構(gòu)與厭氧膜生物反應(yīng)器膜污染之間關(guān)系的研究 [D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2014.

    [37] 趙 霏.再生水濕地香蒲根際細菌多樣性研究—以北京市永定河城市景觀補水段為例 [D]. 北京:首都師范大學(xué), 2014.

    [38] Koski-Vahala J, Hartikainen H, Tallberg P. Phosphorus mobilization from various sediment pools in response to increased pH and silicate concentration [J]. Journal of Environmental Quality, 2001,30(2):546–552.

    [39] Yannarell A C, Triplett E W. Geographic and environmental sources of variation in lake bacterial community composition [J]. Applied and Environmental Microbiology, 2005,71(1):227–239.

    [40] 左新宇,梁運祥.微囊藻與硝化細菌在不同種群密度條件下的相互作用 [J]. 環(huán)境科學(xué)與技術(shù), 2013,36(12):65-70.

    [41] 黃 睿,沈 烽,羅 娟,等.藍藻水華消亡對湖泊表層沉積物中氨氧化細菌豐度和群落結(jié)構(gòu)的影響 [J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2015,31(3):334-339.

    [42] de Beer D, Wenzh?fer F, Ferdelman T, et al. Transport and mineralization rates in North Sea sandy intertidal sediments, Sylt-Romo Basin, Wadden Sea [J]. Limnology & Oceanography, 2005,50:113–127.

    [43] 王 珍.渭河(陜西段)潛流帶水交換對TN、TP遷移轉(zhuǎn)化的影響研究 [D]. 西安:西北大學(xué), 2014.

    致謝:感謝北京大學(xué)環(huán)境學(xué)院謝曙光老師對本研究中分子微生物分析方法的指導(dǎo).

    Influences of sediment characteristics and overlying water quality on sediment bacterial communities in a seasonal sandy river.

    LIU Rui, ZHOU Xiao-de*

    (State Key Laboratory of Ecological Water Conservancy in the Northwest Arid Area, Xi’an University of Technology, Xi’an 710048, China)., 2017,37(11):4342~4352

    The seasonal change of the diversity and structure of sediment bacterial communities in sandy Weihe River and the influences of sediment particle characteristic, sediment nutrient content and the quality of overlaying water were investigated using terminal-restriction fragment length polymorphism (T-RFLP) analysis. The results of redundancy analysis (RDA) revealed significant seasonal variations of sediment characteristics and water quality. Total nitrogen content in sediment increased one order of magnitude with a decrease of spatial heterogeneity (RSD=17%) in wet season, while N: P in sediment was moderately correlated with N:P in overlying water. Based on T-RFLP, the Shannon index of sediment bacterial communities in mainstream ranged between 2.02 and 3.18. Spatial change pattern of bacterial diversity was also considerably different between hydrological seasons. Sediment bacterial communities in wet season had the strongest heterogeneity (RSD=37.4%) and the lowest Bray-Curtis similarity (17.4%). Most of the dominant fragments belong to species within the phyla of Proteobacteria and Bacteroidetes. Properties of sediment strongly influenced bacterial communities in normal and dry season, while the quality of overlying water had a close relationship with sediment bacteria community in wet season. All these results indicated the complexity of sediment bacterial community variation in seasonally sandy river, and significant impacts of sediment characteristics and overlying water quality.

    Weihe River;sediment;bacterial community;T-RFLP

    X172

    A

    1000-6923(2017)11-4342-11

    劉 睿(1985-),女,河南泌陽人,西安理工大學(xué)博士研究生,研究方向為環(huán)境微生物.發(fā)表論文5篇.

    2017-04-16

    教育部高等學(xué)校博士學(xué)科點專項科研基金(20136118120022);陜西省科技統(tǒng)籌創(chuàng)新工程計劃項目(2013SZS02- P01)

    *責(zé)任作者, 教授, zhouxd@xaut.edu.cn

    猜你喜歡
    豐水期渭河樣點
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    牛欄江流域浮游動物群落結(jié)構(gòu)特征及驅(qū)動因子
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調(diào)查分析
    小麥條銹病田間為害損失的初步分析
    湖北植保(2022年4期)2022-08-23 10:51:52
    基于空間模擬退火算法的最優(yōu)土壤采樣尺度選擇研究①
    土壤(2021年1期)2021-03-23 07:29:06
    退耕還林工程對渭河洪峰流量的影響
    基于分融策略的土壤采樣設(shè)計方法*
    渭河故道恢復(fù)濕地生態(tài)的引水方案創(chuàng)新
    滇池流域入湖河流水文形貌特征對豐水期大型底棲動物群落結(jié)構(gòu)的影響
    1996—2009年渭河干流氮素污染特征
    在线观看66精品国产| 特级一级黄色大片| 亚洲色图av天堂| 国语自产精品视频在线第100页| 特级一级黄色大片| 国内精品美女久久久久久| 神马国产精品三级电影在线观看| 黄片wwwwww| 亚洲av福利一区| 国产高清国产精品国产三级 | 色吧在线观看| 中文精品一卡2卡3卡4更新| av在线观看视频网站免费| 亚洲在线自拍视频| 国产久久久一区二区三区| 久久人妻av系列| 青春草视频在线免费观看| 国产午夜精品一二区理论片| 免费无遮挡裸体视频| 日本黄色视频三级网站网址| 日产精品乱码卡一卡2卡三| 精华霜和精华液先用哪个| 丰满人妻一区二区三区视频av| 丰满乱子伦码专区| 日韩高清综合在线| 久久精品国产99精品国产亚洲性色| 18禁动态无遮挡网站| 观看免费一级毛片| 免费人成在线观看视频色| 国产又色又爽无遮挡免| 国产在线男女| av天堂中文字幕网| 中文天堂在线官网| 搡女人真爽免费视频火全软件| 国产精品一区二区三区四区久久| 一二三四中文在线观看免费高清| 午夜爱爱视频在线播放| 亚洲精品一区蜜桃| 亚洲性久久影院| www.av在线官网国产| 亚洲中文字幕日韩| 亚洲成人精品中文字幕电影| 久久婷婷人人爽人人干人人爱| 国产成人精品婷婷| 国产精品女同一区二区软件| 国产精品精品国产色婷婷| 性插视频无遮挡在线免费观看| 亚洲欧美日韩高清专用| 超碰97精品在线观看| 尤物成人国产欧美一区二区三区| 少妇熟女aⅴ在线视频| 日本黄色视频三级网站网址| 精华霜和精华液先用哪个| 中文精品一卡2卡3卡4更新| 一级毛片久久久久久久久女| 最后的刺客免费高清国语| 中文乱码字字幕精品一区二区三区 | 亚洲av.av天堂| av在线蜜桃| 在线观看av片永久免费下载| av免费在线看不卡| 综合色丁香网| 国产精品久久久久久精品电影小说 | 狂野欧美激情性xxxx在线观看| 六月丁香七月| 国产亚洲一区二区精品| 国产又黄又爽又无遮挡在线| 国产精品一区二区三区四区免费观看| 成人毛片a级毛片在线播放| 天堂√8在线中文| 在线观看美女被高潮喷水网站| 亚洲欧美日韩高清专用| av女优亚洲男人天堂| 大又大粗又爽又黄少妇毛片口| 免费看光身美女| 九九在线视频观看精品| 国产一区二区在线av高清观看| 美女内射精品一级片tv| 久久草成人影院| 久久精品国产99精品国产亚洲性色| 欧美97在线视频| 伦理电影大哥的女人| 精品熟女少妇av免费看| 久久久久久国产a免费观看| 嫩草影院入口| 特大巨黑吊av在线直播| 亚洲丝袜综合中文字幕| 在线a可以看的网站| 又粗又爽又猛毛片免费看| 成人三级黄色视频| 精品久久久久久久久av| 久久午夜福利片| 秋霞在线观看毛片| 久久精品熟女亚洲av麻豆精品 | 国产精品野战在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产自在天天线| or卡值多少钱| 国产亚洲av嫩草精品影院| 成人鲁丝片一二三区免费| 国产高清三级在线| 九九在线视频观看精品| 色吧在线观看| 国产精品久久久久久久久免| 国产毛片a区久久久久| 亚洲精品久久久久久婷婷小说 | 亚洲美女视频黄频| 少妇人妻一区二区三区视频| 在线观看美女被高潮喷水网站| 国产人妻一区二区三区在| 久久久久久久久久久丰满| 一个人看视频在线观看www免费| av.在线天堂| 看非洲黑人一级黄片| 久久精品久久久久久噜噜老黄 | 十八禁国产超污无遮挡网站| 日韩,欧美,国产一区二区三区 | 亚洲va在线va天堂va国产| 国产亚洲一区二区精品| 91狼人影院| 精品不卡国产一区二区三区| 美女国产视频在线观看| 最近中文字幕高清免费大全6| 99热6这里只有精品| 国产高清三级在线| 日本黄大片高清| 亚洲精品成人久久久久久| 长腿黑丝高跟| 亚洲精品国产av成人精品| 午夜福利成人在线免费观看| 国产精品乱码一区二三区的特点| 日韩欧美精品v在线| 国产黄a三级三级三级人| 国产精品久久久久久av不卡| 中文字幕亚洲精品专区| 日本熟妇午夜| 中文欧美无线码| 少妇的逼好多水| 日韩高清综合在线| 国产在线一区二区三区精 | 国产精品.久久久| 一区二区三区高清视频在线| av又黄又爽大尺度在线免费看 | 干丝袜人妻中文字幕| 尤物成人国产欧美一区二区三区| 国产淫语在线视频| 欧美区成人在线视频| 久久99精品国语久久久| 久久久久久久午夜电影| 性色avwww在线观看| 极品教师在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲自拍偷在线| 亚洲欧洲国产日韩| 欧美激情在线99| 在现免费观看毛片| 日本一本二区三区精品| 久久精品91蜜桃| 黑人高潮一二区| 免费搜索国产男女视频| 久久精品久久久久久久性| 免费一级毛片在线播放高清视频| 日韩欧美 国产精品| 国产精品一及| 久久韩国三级中文字幕| 三级国产精品片| 联通29元200g的流量卡| 久久这里有精品视频免费| 老师上课跳d突然被开到最大视频| 欧美另类亚洲清纯唯美| 亚洲久久久久久中文字幕| 99热这里只有是精品50| 午夜免费男女啪啪视频观看| 十八禁国产超污无遮挡网站| 国产成人午夜福利电影在线观看| 国产午夜精品一二区理论片| 两性午夜刺激爽爽歪歪视频在线观看| 免费av不卡在线播放| 一区二区三区高清视频在线| 2022亚洲国产成人精品| 国产国拍精品亚洲av在线观看| 91精品伊人久久大香线蕉| 女的被弄到高潮叫床怎么办| 午夜福利网站1000一区二区三区| 国产一区二区亚洲精品在线观看| 亚洲五月天丁香| 免费一级毛片在线播放高清视频| 寂寞人妻少妇视频99o| 七月丁香在线播放| 在现免费观看毛片| 中文在线观看免费www的网站| 两个人的视频大全免费| 一边亲一边摸免费视频| 国产久久久一区二区三区| 人人妻人人澡人人爽人人夜夜 | 亚洲精品aⅴ在线观看| 国产精品国产三级国产av玫瑰| 午夜福利在线观看吧| 亚洲丝袜综合中文字幕| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 国产中年淑女户外野战色| 看片在线看免费视频| 国产午夜福利久久久久久| 亚洲一级一片aⅴ在线观看| 国产亚洲精品久久久com| 一级爰片在线观看| 男女啪啪激烈高潮av片| 淫秽高清视频在线观看| 女人被狂操c到高潮| 国产黄片美女视频| 色综合亚洲欧美另类图片| 亚洲五月天丁香| av黄色大香蕉| 色网站视频免费| 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 色尼玛亚洲综合影院| 精品国产一区二区三区久久久樱花 | 波多野结衣巨乳人妻| 免费观看a级毛片全部| 六月丁香七月| 丝袜喷水一区| 日本-黄色视频高清免费观看| 亚洲一级一片aⅴ在线观看| 亚洲四区av| 精品久久久噜噜| 国产精品一区二区性色av| 长腿黑丝高跟| 青春草国产在线视频| 国产午夜精品论理片| 国产午夜精品久久久久久一区二区三区| 嫩草影院新地址| 18禁动态无遮挡网站| 看十八女毛片水多多多| 国产成人91sexporn| 久久国内精品自在自线图片| 天堂网av新在线| 寂寞人妻少妇视频99o| 国产成人精品婷婷| 好男人在线观看高清免费视频| 欧美一级a爱片免费观看看| 1000部很黄的大片| 色哟哟·www| 国产真实乱freesex| 欧美97在线视频| 男的添女的下面高潮视频| 激情 狠狠 欧美| 国产69精品久久久久777片| 午夜视频国产福利| 少妇丰满av| 卡戴珊不雅视频在线播放| 别揉我奶头 嗯啊视频| 欧美最新免费一区二区三区| 婷婷六月久久综合丁香| 午夜福利在线观看吧| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 国产av不卡久久| 欧美xxxx性猛交bbbb| 中文字幕av在线有码专区| 国产三级在线视频| av视频在线观看入口| 又爽又黄a免费视频| 插逼视频在线观看| 亚洲熟妇中文字幕五十中出| 最近的中文字幕免费完整| 免费av观看视频| 精品人妻视频免费看| 美女内射精品一级片tv| 夜夜爽夜夜爽视频| 国产亚洲5aaaaa淫片| 视频中文字幕在线观看| 蜜桃久久精品国产亚洲av| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 秋霞在线观看毛片| 国产精品乱码一区二三区的特点| 国产成人91sexporn| 免费观看的影片在线观看| 黄色欧美视频在线观看| 在线观看av片永久免费下载| 91在线精品国自产拍蜜月| 又爽又黄a免费视频| 日韩欧美 国产精品| or卡值多少钱| 日本三级黄在线观看| 三级国产精品欧美在线观看| 亚洲av福利一区| 国产精品久久久久久av不卡| 成人综合一区亚洲| 亚洲激情五月婷婷啪啪| 搡老妇女老女人老熟妇| 尤物成人国产欧美一区二区三区| 男的添女的下面高潮视频| 在线观看美女被高潮喷水网站| 亚洲图色成人| 我的老师免费观看完整版| 亚州av有码| 久久久久久久午夜电影| 国产精品av视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 村上凉子中文字幕在线| 丰满人妻一区二区三区视频av| 综合色丁香网| 亚洲最大成人手机在线| 成人二区视频| 99久久成人亚洲精品观看| 99久久精品一区二区三区| 特级一级黄色大片| 啦啦啦观看免费观看视频高清| 欧美日本视频| av专区在线播放| 国内精品一区二区在线观看| 亚洲内射少妇av| 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费| 午夜福利在线观看吧| av专区在线播放| 亚洲精品成人久久久久久| 国产成人aa在线观看| 伦精品一区二区三区| 久久99热这里只有精品18| 国产老妇伦熟女老妇高清| 高清视频免费观看一区二区 | 91aial.com中文字幕在线观看| 最后的刺客免费高清国语| 看黄色毛片网站| av国产免费在线观看| 午夜亚洲福利在线播放| 老司机福利观看| 国产精品国产三级专区第一集| 男女国产视频网站| 日本一二三区视频观看| 干丝袜人妻中文字幕| 麻豆一二三区av精品| 久久精品人妻少妇| 成年女人永久免费观看视频| 亚洲精品国产av成人精品| 看片在线看免费视频| 蜜臀久久99精品久久宅男| 国产精品日韩av在线免费观看| 一个人免费在线观看电影| 日本欧美国产在线视频| 男人狂女人下面高潮的视频| 久久精品熟女亚洲av麻豆精品 | 国产成人aa在线观看| 岛国毛片在线播放| 亚洲精品国产成人久久av| 国产精品久久电影中文字幕| 久久精品人妻少妇| 最近手机中文字幕大全| 国产黄片视频在线免费观看| 亚洲在线观看片| 国产一区二区三区av在线| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 少妇的逼好多水| 国产淫片久久久久久久久| 91狼人影院| 我的女老师完整版在线观看| 建设人人有责人人尽责人人享有的 | 国产精品一区二区三区四区免费观看| 99在线人妻在线中文字幕| 蜜桃亚洲精品一区二区三区| 中文亚洲av片在线观看爽| 性色avwww在线观看| av在线亚洲专区| 色噜噜av男人的天堂激情| 狂野欧美白嫩少妇大欣赏| av天堂中文字幕网| av在线播放精品| www.av在线官网国产| 成年女人永久免费观看视频| 国产综合懂色| 身体一侧抽搐| 国产人妻一区二区三区在| 在线a可以看的网站| av在线天堂中文字幕| 国产在视频线在精品| 精品一区二区免费观看| 桃色一区二区三区在线观看| 日韩欧美 国产精品| 岛国在线免费视频观看| 日本一本二区三区精品| 日韩精品有码人妻一区| 日本黄色片子视频| 国内精品美女久久久久久| av卡一久久| 久久久欧美国产精品| 草草在线视频免费看| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 色噜噜av男人的天堂激情| 人妻少妇偷人精品九色| 高清午夜精品一区二区三区| 亚洲av一区综合| 夜夜爽夜夜爽视频| 日韩制服骚丝袜av| 干丝袜人妻中文字幕| 成人毛片60女人毛片免费| 国产精品久久久久久av不卡| 国产伦在线观看视频一区| 亚州av有码| 亚洲欧美精品综合久久99| 身体一侧抽搐| 性插视频无遮挡在线免费观看| 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| 韩国av在线不卡| 高清午夜精品一区二区三区| 成人午夜高清在线视频| 午夜视频国产福利| 丰满人妻一区二区三区视频av| 久久精品久久久久久久性| 寂寞人妻少妇视频99o| 久久6这里有精品| 欧美变态另类bdsm刘玥| 免费av毛片视频| 精品无人区乱码1区二区| 天堂中文最新版在线下载 | 丰满人妻一区二区三区视频av| 国产视频内射| 亚洲av日韩在线播放| 卡戴珊不雅视频在线播放| 尾随美女入室| 日韩欧美三级三区| 欧美成人免费av一区二区三区| 男人舔奶头视频| 能在线免费看毛片的网站| kizo精华| 99久久中文字幕三级久久日本| 建设人人有责人人尽责人人享有的 | 日本免费a在线| 青春草视频在线免费观看| 日韩中字成人| 熟女人妻精品中文字幕| 中文亚洲av片在线观看爽| 69av精品久久久久久| av专区在线播放| 国产精品av视频在线免费观看| 嫩草影院新地址| 亚洲精华国产精华液的使用体验| 3wmmmm亚洲av在线观看| 国产人妻一区二区三区在| av在线观看视频网站免费| 18禁动态无遮挡网站| 国产成人91sexporn| 欧美97在线视频| 嫩草影院精品99| 一个人看视频在线观看www免费| 六月丁香七月| 国产精品乱码一区二三区的特点| 国产亚洲最大av| 级片在线观看| 久久久精品94久久精品| 免费看光身美女| 波野结衣二区三区在线| 国产精品一区二区三区四区久久| 免费黄网站久久成人精品| 国产av在哪里看| 国产精品嫩草影院av在线观看| 黄色一级大片看看| 久久久精品欧美日韩精品| 美女被艹到高潮喷水动态| 免费不卡的大黄色大毛片视频在线观看 | 高清日韩中文字幕在线| 久久久久久久国产电影| 亚洲自拍偷在线| 久久精品国产亚洲网站| 午夜日本视频在线| 午夜精品在线福利| 亚洲经典国产精华液单| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久人妻蜜臀av| 日韩欧美三级三区| 亚洲图色成人| 欧美bdsm另类| 精品99又大又爽又粗少妇毛片| 精品人妻偷拍中文字幕| 亚洲综合精品二区| 久久鲁丝午夜福利片| 国产高清国产精品国产三级 | 精华霜和精华液先用哪个| 欧美高清性xxxxhd video| 97人妻精品一区二区三区麻豆| 国产成人a∨麻豆精品| 亚洲,欧美,日韩| 亚洲av电影在线观看一区二区三区 | 日日摸夜夜添夜夜添av毛片| 国产精品三级大全| 久久久a久久爽久久v久久| 热99在线观看视频| 成年av动漫网址| 中文欧美无线码| 在线免费观看不下载黄p国产| 天堂中文最新版在线下载 | 免费在线观看成人毛片| 午夜福利网站1000一区二区三区| 久久久久九九精品影院| 黄色日韩在线| 欧美激情在线99| 内射极品少妇av片p| 亚洲精品国产成人久久av| a级一级毛片免费在线观看| 国产精品人妻久久久影院| 国产精品久久久久久久久免| av福利片在线观看| 久久6这里有精品| 久久精品国产亚洲av涩爱| 久久综合国产亚洲精品| 小蜜桃在线观看免费完整版高清| 欧美一区二区国产精品久久精品| 国产色爽女视频免费观看| 免费电影在线观看免费观看| 精品久久久久久成人av| 村上凉子中文字幕在线| 超碰97精品在线观看| 欧美三级亚洲精品| 亚洲四区av| 国产精品一区二区三区四区免费观看| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 国产av不卡久久| 99热这里只有精品一区| 久久久亚洲精品成人影院| 好男人在线观看高清免费视频| 精品熟女少妇av免费看| 能在线免费看毛片的网站| 三级男女做爰猛烈吃奶摸视频| 嫩草影院入口| 国产黄片美女视频| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 校园人妻丝袜中文字幕| 亚洲国产日韩欧美精品在线观看| 国产亚洲最大av| 免费一级毛片在线播放高清视频| 99在线人妻在线中文字幕| 成人综合一区亚洲| 精品久久久久久久久av| 欧美成人午夜免费资源| 国产 一区精品| 婷婷色av中文字幕| 长腿黑丝高跟| 成年版毛片免费区| 69av精品久久久久久| 特级一级黄色大片| 超碰av人人做人人爽久久| 看黄色毛片网站| 男女国产视频网站| 夫妻性生交免费视频一级片| 最后的刺客免费高清国语| 成人亚洲欧美一区二区av| 99热这里只有是精品在线观看| 亚洲真实伦在线观看| 51国产日韩欧美| 久久久久久久久大av| 中文字幕亚洲精品专区| 国产精品日韩av在线免费观看| 免费av观看视频| 青春草亚洲视频在线观看| 久久精品91蜜桃| 久久久国产成人免费| 久久精品91蜜桃| 国产单亲对白刺激| 久久精品久久久久久噜噜老黄 | 一区二区三区高清视频在线| 欧美极品一区二区三区四区| 国产精品一及| 啦啦啦观看免费观看视频高清| 免费观看性生交大片5| 国产色婷婷99| 亚洲国产精品专区欧美| 国产一区有黄有色的免费视频 | 国产精品一区二区在线观看99 | 一级毛片电影观看 | 欧美又色又爽又黄视频| 久99久视频精品免费| 男的添女的下面高潮视频| 亚洲精品一区蜜桃| 国产一区亚洲一区在线观看| 99热这里只有是精品50| 嫩草影院新地址| 三级男女做爰猛烈吃奶摸视频| 2021少妇久久久久久久久久久| 人体艺术视频欧美日本| 一本一本综合久久| 男的添女的下面高潮视频| 国产淫片久久久久久久久| 男人舔女人下体高潮全视频| 真实男女啪啪啪动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | h日本视频在线播放| 一边摸一边抽搐一进一小说| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 99久久精品热视频| 国产69精品久久久久777片| 日本wwww免费看| 一级爰片在线观看| 视频中文字幕在线观看| 成人亚洲欧美一区二区av| 欧美性猛交╳xxx乱大交人| 成年女人看的毛片在线观看| 久久久国产成人免费| 日产精品乱码卡一卡2卡三| 国产精品麻豆人妻色哟哟久久 | 免费播放大片免费观看视频在线观看 | 国产精品美女特级片免费视频播放器| 亚洲精品,欧美精品| 女人被狂操c到高潮| 国国产精品蜜臀av免费| 久久久久免费精品人妻一区二区| 亚洲国产精品成人久久小说|