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

    濕地蘆葦根結(jié)合好氣細(xì)菌群落時(shí)空分布及其與水質(zhì)因子的關(guān)系

    2013-09-15 06:32:26薇1郭逍宇1霏1
    生態(tài)學(xué)報(bào) 2013年5期
    關(guān)鍵詞:凈化系統(tǒng)類群蘆葦

    熊 薇1,郭逍宇1,*,趙 霏1

    (1.首都師范大學(xué)資源環(huán)境與旅游學(xué)院,北京100048;2.北京市城市環(huán)境過程與數(shù)字模擬重點(diǎn)實(shí)驗(yàn)室-省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地,北京100048;3.三維信息獲取與應(yīng)用教育部重點(diǎn)實(shí)驗(yàn)室,北京100048;4.資源環(huán)境與地理信息系統(tǒng)北京市重點(diǎn)實(shí)驗(yàn)室,北京100048)

    再生水作為一種穩(wěn)定的再生水源,在緩解北京市濕地水資源供需矛盾、減輕濕地水環(huán)境污染方面發(fā)揮了重要作用[1-3],據(jù)統(tǒng)計(jì)到2010年北京市再生水利用6.8億m3,其中河湖環(huán)境利用2.1億m3,占再生水利用總量的30.9%[4],但再生水中豐富的氮、磷含量決定了再生水利用必然會(huì)產(chǎn)生許多負(fù)面環(huán)境效應(yīng),再加上城市濕地多為封閉性水體,水流速度較緩,水體更新周期較長(zhǎng),水體生態(tài)系統(tǒng)結(jié)構(gòu)較為單一,使得水體富營(yíng)養(yǎng)化成為濕地再生水利用的最大障礙[5]。

    城市人工濕地不僅具有重要的景觀作用,而且具有良好的環(huán)境污染修復(fù)能力,能有效去除水體氮、磷、各種有機(jī)物質(zhì)、重金屬氧化物,降低水體生化需氧量(BOD)、總懸浮固體(TSS)含量[6-7],植物在人工濕地氮污染修復(fù)過程中扮演重要角色[8],Oostrom等認(rèn)為濕地中13%的氮是通過基質(zhì)的過濾和植物吸收作用被去除的,濕地植物的凈化作用主要是通過植物群落根際分泌物及凋落物的分解影響微生物群落的功能,進(jìn)而影響系統(tǒng)的凈化效能[9]。Li等采用基因文庫技術(shù)對(duì)北京翠湖人工濕地蘆葦和香蒲根內(nèi)生細(xì)菌群落結(jié)構(gòu)分別進(jìn)行研究,結(jié)果表明濕地植物根內(nèi)生菌群落對(duì)濕地氮、磷、硫以及有機(jī)復(fù)合物具重要的降解作用,對(duì)濕地修復(fù)具重要意義[10-11]。Atif等運(yùn)用變性梯度凝膠電泳(DGGE)技術(shù)研究了兩個(gè)典型垂直流人工濕地植物凋落物及沉積物中氮去除功能細(xì)菌的群落結(jié)構(gòu),結(jié)果表明功能細(xì)菌群落因基質(zhì)不同及采樣點(diǎn)位置不同而具有顯著差異,功能細(xì)菌群落多樣性與濕地凈化系統(tǒng)功能作用相關(guān)[12]。大量研究表明濕地植物根圍細(xì)菌多樣性與水體物理化學(xué)性質(zhì)、濕地凈化效率、人工濕地系統(tǒng)穩(wěn)定性具顯著相關(guān)性[5,13-14],且具顯著的季節(jié)動(dòng)態(tài)特征[15-16],但從環(huán)境因素角度揭示濕地微生物群落結(jié)構(gòu)差異機(jī)制的研究尚鮮見報(bào)道。

    借助于生態(tài)學(xué)多樣性測(cè)度方法研究復(fù)雜多樣的微生物群落結(jié)構(gòu)已經(jīng)較為成熟[17-19],而濕地微生物多樣性信息主要依靠培養(yǎng)法和非培養(yǎng)法如基于聚合酶鏈?zhǔn)椒磻?yīng)(PCR)的限制性片段長(zhǎng)度多態(tài)性性分析(ARDRA)、變性梯度凝膠電泳(DGGE)、熒光原位雜交(FISH)技術(shù)獲得。盡管傳統(tǒng)的微生物培養(yǎng)法對(duì)于研究微生物群落結(jié)構(gòu)有不足之處,但培養(yǎng)法不僅對(duì)微生物資源的開發(fā)與利用具重要作用,而且可通過培養(yǎng)基的選擇作用有針對(duì)性的反映人工濕地中具有不同分解作用的細(xì)菌類群,如在普通牛肉膏蛋白胨培養(yǎng)基上生長(zhǎng)的細(xì)菌可能代表了對(duì)含氮有機(jī)污染物具有較強(qiáng)分解能力的微生物類群,而有機(jī)氮素如蛋白質(zhì)、多肽、氨基酸和尿素等是富營(yíng)養(yǎng)水體的主要污染元素,這一類群的細(xì)菌可能在氮污染水體的治理上具有重要的應(yīng)用潛力[20],對(duì)研究再生水補(bǔ)水濕地微生物凈化機(jī)理具有重要作用?;诖吮狙芯恳员本W林匹克公園再生水人工濕地凈化系統(tǒng)為例,借助細(xì)菌培養(yǎng)與ARDRA相結(jié)合的技術(shù)分析研究夏秋兩季的凈化系統(tǒng)不同功能區(qū)濕地蘆葦根結(jié)合好氣細(xì)菌的數(shù)量、群落結(jié)構(gòu)及其多樣性,并借助典范對(duì)應(yīng)分析(CCA)研究水質(zhì)變量對(duì)群落結(jié)構(gòu)的影響,揭示影響細(xì)菌群落結(jié)構(gòu)特征變異的主要水質(zhì)參數(shù)[21],為再生水城市濕地利用提供合理的科學(xué)依據(jù)。

    1 研究區(qū)概況

    奧林匹克公園人工濕地凈化系統(tǒng)地處北京市北部,位于奧林匹克公園園內(nèi),于2007年開始運(yùn)行。該區(qū)屬暖溫帶半干旱半濕潤(rùn)大陸性季風(fēng)氣候。多年平均氣溫11—12℃,極端最高氣溫41.6℃,極端最低氣溫-21.2℃。多年平均水面蒸發(fā)量在1200 mm左右,多年平均降雨量約600 mm。降雨年內(nèi)分配不均,多集中在汛期6—9月,約占全年的85%。該凈化系統(tǒng)以兩個(gè)復(fù)合垂直流人工濕地為核心,以氧化塘為補(bǔ)充,形成多層次的凈化系統(tǒng),將來自清河再生水廠的再生水凈化處理后作為奧林匹克公園主湖補(bǔ)充水源。其中,復(fù)合垂直流人工濕地由垂直向下流和垂直向上流兩部分濕地構(gòu)成。濕地系統(tǒng)內(nèi)水生植物主要包括蘆葦、香蒲、水蔥、千屈菜等,植被覆蓋度較高,可達(dá)70%—90%。凈化系統(tǒng)各處理區(qū)分布及水流方向如圖1[22]。再生水作為補(bǔ)充水源首先進(jìn)入復(fù)合垂直流人工濕地再生水處理區(qū)(RW),同時(shí)主湖水循環(huán)至復(fù)合垂直流人工濕地循環(huán)水處理區(qū)(CW),再生水處理出水經(jīng)氧化塘區(qū)(OP)深化處理后進(jìn)入主湖,主湖水停留后會(huì)進(jìn)行下一次循環(huán)。

    圖1 北京奧林匹克公園再生水人工濕地凈化系統(tǒng)示意圖Fig.1 Location and distribution of the purification system in Beijing Olympic Park

    2 材料與方法

    2.1 蘆葦根系樣品采集及水質(zhì)測(cè)定

    選取研究區(qū)不同處理區(qū)典型濕地植物蘆葦為研究對(duì)象,于2010年11月9日、2011年6月11日在人工濕地凈化系統(tǒng)中RW、CW、OP功能區(qū)分別獲得蘆葦根系樣品,新鮮表面附著基質(zhì)的蘆葦根系樣品置于無菌三角瓶中密封,立即帶回實(shí)驗(yàn)室處理分析。并對(duì)樣品進(jìn)行編號(hào),夏季編號(hào)為SRW、SCW、SOP,秋季各功能區(qū)樣品編號(hào)為ARW、ACW、AOP。于夏、秋兩季分別在3個(gè)功能區(qū)均勻布3個(gè)樣點(diǎn),采集水樣,分析與氮污染相關(guān)的水體富營(yíng)養(yǎng)化指標(biāo)。水質(zhì)樣品低溫儲(chǔ)藏帶回實(shí)驗(yàn)室處理。其中,水體化學(xué)需氧量(COD)、總磷(TP)采用國(guó)標(biāo)法測(cè)定,NH+4、NO-3濃度由離子色譜儀(ICS-1500型、ICS-90型,美國(guó)Dionex)測(cè)定,水體總氮(TN)采用耶拿Multi N/CUV分析儀測(cè)定。其余水質(zhì)參數(shù),包括葉綠素a(CHL)、溶解氧(DO)、氯離子(Cl-)、藍(lán)綠藻(PCY)、鹽度(Sal)、總?cè)芙庑怨腆w(TDS)、電導(dǎo)系數(shù)(SPC)、氧化還原電位(ORP)由水質(zhì)儀(Hydrolab Datasonde5 5 X,the United States)現(xiàn)場(chǎng)測(cè)定,測(cè)定結(jié)果見表1。

    表1 夏秋兩季水質(zhì)理化性質(zhì)Table 1 Water quality of constructed wetland in autumn and summer

    2.2 根結(jié)合菌細(xì)菌培養(yǎng)、計(jì)數(shù)、分離純化

    將新鮮表面附著基質(zhì)的蘆葦根系置于震蕩儀震蕩去除表面附著物,取表面清洗干凈的蘆葦根系1 g,再用無菌水清洗3次,無菌濾紙吸取多余水分,置于滅菌研缽內(nèi),加入少量無菌石英砂及少量無菌水研磨后定容100 mL。此液體依次稀釋涂板,采用LB培養(yǎng)基培養(yǎng),每梯度3個(gè)平行,在30℃條件下倒置培養(yǎng)48 h,記錄菌落數(shù)目,計(jì)算可培養(yǎng)細(xì)菌的平均數(shù)目。選擇菌落數(shù)在30—300范圍的稀釋度平板,將所有單菌落分別轉(zhuǎn)接劃線LB培養(yǎng)基平板上,分離純化。

    2.3 擴(kuò)增16SrDNA限制性分析

    2.3.1 PCR模板制備

    用無菌牙簽挑取少量已純化的菌體重懸于20μL無菌ddH2O中,置于沸水浴中煮沸10 min,立即置于冰上5 min,4℃,10000 rpm離心5 min,置于冰上,使用時(shí)取2.0μL上清液做模板[23-24]。

    2.3.2 16SrDNA擴(kuò)增

    采用細(xì)菌通用引物765f:5'-GGAGCAAACAGGATTTGATA-3',和1492r:5'-GGTTACCTTGTTACGACTT-3'[25]對(duì)細(xì)菌DNA進(jìn)行擴(kuò)增。25μL擴(kuò)增體系包含:12.5μL 2×Taq PCR Master Mix(天根生化科技有限公司,KT201-12),1μL 10μmol/L引物765f和1μL 10μmol/L 1492r(生工生物工程(上海)有限公司),2μL DNA模板,8.5μL無菌水。PCR反應(yīng)程序94℃,5 min;30個(gè)循環(huán)包括94℃,30 s,55℃,1 min,72℃,1 min;72℃ 延伸10 min;4℃保存。1.0%的瓊脂糖凝膠檢測(cè)PCR產(chǎn)物。

    2.3.3 擴(kuò)增片段酶切分型

    用限制性內(nèi)切酶Hae III對(duì)PCR產(chǎn)物酶切,10μL酶切體系:5μL PCR產(chǎn)物,0.3μL Hae III(10 U/μL),1 μL 10×buffer,3.7μL雙蒸水。37℃反應(yīng)4 h。酶切產(chǎn)物用2.0%瓊脂糖凝膠電泳,1×TAE電泳緩沖液,100 V電泳60 min。比較分析酶切電泳圖譜。

    2.4 統(tǒng)計(jì)數(shù)據(jù)分析

    據(jù)酶切圖譜信息,將酶切圖譜完全一致的菌落劃分為1個(gè)操作分類單元(operational taxonomic unit,OTU)[26],據(jù)此估計(jì)樣品中存在的最低限度的細(xì)菌數(shù)目,并運(yùn)用SPSS18.0將菌落聚類。依此分類分析細(xì)菌群落結(jié)構(gòu)及多樣性分析。利用Shannon-Weaver指數(shù)(H')評(píng)價(jià)細(xì)菌群落綜合多樣性,以EI計(jì)算均勻度指數(shù),以Margale指數(shù)(D')評(píng)價(jià)群落豐富度。CCA排序圖由Canoco for Windows 4.5和TWINSPAN軟件包共同完成,其余圖表由Microsoft Excel 2007完成。多樣性評(píng)價(jià)指數(shù)計(jì)算公式如下:

    式中,Pi=ni/N,S是圖譜的分類數(shù),ni是第i類包含的圖譜,N為總圖譜數(shù)。

    3 結(jié)果與分析

    3.1 細(xì)菌群落密度分析

    人工濕地微生物群落密度與污水凈化效率具有密切關(guān)系,微生物數(shù)量與污染物去除率呈正相關(guān)[14,20,27]。群落密度分析結(jié)果見表2,從季節(jié)變異看,夏季蘆葦根結(jié)合微生物群落密度較秋季增加76.37%—96.70%,具有明顯的季節(jié)效應(yīng)。從濕地類型差異來看,再生水和循環(huán)水經(jīng)復(fù)合垂直流人工濕地處理進(jìn)入氧化塘后,微生物群落密度減少。從補(bǔ)水水質(zhì)差異看,再生水處理區(qū)微生物群落密度高于循環(huán)水處理區(qū),夏秋兩季在不同功能區(qū)表現(xiàn)出相同的分布格局。為進(jìn)一步分析環(huán)境因子對(duì)蘆葦根結(jié)合細(xì)菌群落密度的影響,采用spearman秩相關(guān)系數(shù)分析(表4)水環(huán)境因子對(duì)群落密度的影響,結(jié)果顯示凈化系統(tǒng)水體濃度(0.934**,P﹤0.001)和ORP(0.886*,P﹤0.005)與群落密度顯著相關(guān),是影響細(xì)菌密度時(shí)空差異的主要因素。

    表2 蘆葦根結(jié)合細(xì)菌群落密度Table 2 Community density of the different water treatment units

    3.2 蘆葦根結(jié)合細(xì)菌群落結(jié)構(gòu)的時(shí)空變異分析

    通過聚類分析,將獲得138株分離物劃分為29個(gè)OTUs,每個(gè)OTUs至少代表一種細(xì)菌,表明研究區(qū)至少分布29種不同細(xì)菌(圖2)。其中,16個(gè)OTUs屬于夏季,14個(gè)OTUs屬于秋季,OTU28是夏秋兩季唯一的共有OTU,該類群從夏季優(yōu)勢(shì)類群過渡為秋季偶見類群。夏季16個(gè)OTUs中,11個(gè)OTUs屬于再生水處理區(qū),8個(gè)OTUs屬于循環(huán)水處理區(qū),4個(gè)OTUs屬于氧化塘,其中OTU1,OTU28是夏季3個(gè)功能區(qū)共有類群,OTU1從再生水處理區(qū)和循環(huán)水處理區(qū)的優(yōu)勢(shì)類群過渡為氧化塘的非優(yōu)勢(shì)類群,OTU28從再生水區(qū)和氧化塘的優(yōu)勢(shì)類群過渡為循環(huán)水區(qū)的偶見類群。11個(gè)OTUs屬于各功能區(qū)特有OTU,特有OTU均為非優(yōu)勢(shì)類群或偶見類群,其中8個(gè)OTUs屬于再生水處理區(qū),3個(gè)OTUs屬于循環(huán)水處理區(qū)。秋季14個(gè)OTUs中,6個(gè)OTUs屬于再生水處理區(qū),5個(gè)OTUs屬于循環(huán)水處理區(qū),10個(gè)OTUs屬于氧化塘,其中,OTU12,OTU18是秋季3個(gè)功能區(qū)的共有類群,OTU12從再生水處理區(qū)和循環(huán)水處理區(qū)的偶見類群過渡為氧化塘的非優(yōu)勢(shì)類群,OTU18是秋季3個(gè)功能區(qū)的共有優(yōu)勢(shì)類群。9個(gè)OTUs屬于各功能區(qū)特有OTU,特有OTU均為非優(yōu)勢(shì)類群或偶見類群,其中2個(gè)OTUs屬于再生水處理區(qū),2個(gè)OTUs屬于循環(huán)水處理區(qū),5個(gè)OTUs屬于氧化塘??梢娂竟?jié)動(dòng)態(tài)對(duì)微生物種群類型分布的響應(yīng)體現(xiàn)在夏秋兩季微生物豐富度差別不大,但夏季群落均勻性高于秋季,及共優(yōu)種比例增高,優(yōu)勢(shì)種群類型改變。

    圖2 不同OTUs分離菌株在總分離菌株中的分布Fig.2 Bacterial number distribution of OTUs in all bacterial number

    根據(jù)單個(gè)OTUs包含菌株與該樣品中總菌株數(shù)的比例,對(duì)菌株情況進(jìn)行了統(tǒng)計(jì)(表3)。季節(jié)動(dòng)態(tài)對(duì)蘆葦根結(jié)合好氣細(xì)菌群落結(jié)構(gòu)影響主要體現(xiàn)在:其一,夏季較秋季優(yōu)勢(shì)類群、非優(yōu)勢(shì)類群及偶見類群3大類群的OTU類型多樣化;其二,夏季較秋季偶見類群數(shù)量百分比呈下降趨勢(shì)。蘆葦根結(jié)合好氣細(xì)菌群落結(jié)構(gòu)空間差異夏季體現(xiàn)為:其一,復(fù)合垂直流人工濕地(再生水處理區(qū)、循環(huán)水處理區(qū))較氧化塘非優(yōu)勢(shì)類群和偶見類群類型多樣化;其二,復(fù)合垂直流人工濕地(再生水處理區(qū)、循環(huán)水處理區(qū))較氧化塘非優(yōu)勢(shì)類群及偶見類群數(shù)量百分比增加,其三,不同功能區(qū)之間優(yōu)勢(shì)類群差異增大,再生水處理區(qū)為OTU1和OTU28,循環(huán)水處理區(qū)為OTU1,氧化塘為OTU7和OTU28。三者共同作用導(dǎo)致復(fù)合垂直流人工濕地豐富度及均勻度均增高。蘆葦根結(jié)合好氣細(xì)菌群落結(jié)構(gòu)空間差異秋季體現(xiàn)為:其一,氧化塘較復(fù)合垂直流人工濕地豐富度增高,均勻度變化規(guī)律不明顯;其二,不同功能區(qū)之間優(yōu)勢(shì)類群趨于一致,均為OTU18。綜合不同OTUs分離菌株在總分離菌株中的分布及細(xì)菌群落結(jié)構(gòu)統(tǒng)計(jì)分析結(jié)果,多樣化的類型及其較低頻次的偶見類群共同作用使得夏季較秋季微生物群落結(jié)構(gòu)的豐富度及均勻度均呈現(xiàn)升高趨勢(shì)。共有的優(yōu)勢(shì)類群多度降低,特有的偶見類群增加,同時(shí)偶見類群的多度減少。偶見種類群增加主要體現(xiàn)在特有種類群增多。

    表3 凈化系統(tǒng)可培養(yǎng)蘆葦根結(jié)合細(xì)菌群落結(jié)構(gòu)的時(shí)空差異Table 3 Spatial variability of bacterial community structure in reclaimed water management system

    為解釋水質(zhì)變量對(duì)群落結(jié)構(gòu)分布格局的影響,采用Spearman秩相關(guān)分析群落組成類群的數(shù)量百分比與主要水質(zhì)指標(biāo)之間的關(guān)系。結(jié)果如表4,優(yōu)勢(shì)類群所占百分比與PCY濃度相關(guān)(r=0.771*,P<0.05),偶見類群所占百分比與ORP(r=0.771*,P<0.05)、TP(r=0.771*,P<0.01)、TN(r=0.771*,P<0.05)、NH+4(r=0.829*,P<0.05)相關(guān)。相關(guān)分析結(jié)果表明細(xì)菌優(yōu)勢(shì)類群是水質(zhì)凈化的主要貢獻(xiàn)者,優(yōu)勢(shì)類群對(duì)消除水體富營(yíng)養(yǎng)化起重要作用,與水體富營(yíng)養(yǎng)化緊密相關(guān)。在一定條件下,細(xì)菌群落優(yōu)勢(shì)類群百分比隨水體PCY濃度增高而增加。細(xì)菌偶見類群一般因環(huán)境擾動(dòng)導(dǎo)致新物種的產(chǎn)生或者已存物種的消亡而出現(xiàn),濕地水體營(yíng)養(yǎng)水平及氧化還原環(huán)境的改變是偶見類群出現(xiàn)的擾動(dòng)因素,其改變導(dǎo)致新的適應(yīng)性類群出現(xiàn)及非偶見類群數(shù)量減少,與偶見類群緊密相關(guān)。

    表4 水質(zhì)指標(biāo)與蘆葦根結(jié)合細(xì)菌群落特征的相關(guān)性分析Table 4 Correlation analysis of Water quality and bacterial community characteristics

    3.3 蘆葦根結(jié)合細(xì)菌群落多樣性分析

    采用多樣性指數(shù)進(jìn)一步分析蘆葦根結(jié)合細(xì)菌群落的時(shí)空變異,分析結(jié)果見表5。其中,Shannon-Weaver指數(shù)受群落物種豐富度影響較大,Margale指數(shù)反映群落物種豐富度,均勻度指數(shù)反應(yīng)群落物種均勻性大小。從表中可以看出,夏季較秋季復(fù)合垂直流人工濕地好氣蘆葦根結(jié)合細(xì)菌群落Shannon-Weaver指數(shù)、豐富度指數(shù)顯著增高,均勻度指數(shù)變化不大,氧化塘區(qū)好氣細(xì)菌群落均勻度指數(shù)增加,Shannon-Weaver指數(shù)及Margale指數(shù)均減少,復(fù)合垂直流人工濕地與氧化塘細(xì)菌多樣性具不同的季節(jié)差異。夏季再生水凈化系統(tǒng)中,再生水處理區(qū)具最高Shannon-Weaver指數(shù)、Margale指數(shù),氧化塘區(qū)Shannon-Weaver指數(shù)、Margale指數(shù)處于最低水平,但是在秋季再生水凈化系統(tǒng)中氧化塘區(qū)具最高Shannon-Weaver指數(shù)、Margale指數(shù),循環(huán)水區(qū)Shannon-Weaver指數(shù)、Margale指數(shù)具最低水平??偟膩砜?,氧化塘細(xì)菌群落多樣性受季節(jié)動(dòng)態(tài)影響最小,而復(fù)合垂直流人工濕地細(xì)菌群落多樣性季節(jié)變異較大。再生水處理區(qū)較循環(huán)水處理區(qū)夏秋兩季細(xì)菌群落多樣,豐富度及均勻性均增高。

    表5 秋夏兩季可培養(yǎng)蘆葦根結(jié)合細(xì)菌多樣性分析Table 5 Diversity of bacterial community of three water treatment units and two seasons

    3.4 CCA排序分析

    為進(jìn)一步分析再生水濕地凈化系統(tǒng)細(xì)菌群落結(jié)構(gòu)與水質(zhì)參數(shù)響應(yīng)關(guān)系,結(jié)合運(yùn)用典范對(duì)應(yīng)分析(CCA)排序法和雙向指示種分類法(TWINSPAN),分析夏秋兩季細(xì)菌群落結(jié)構(gòu)特征及其成因。以6樣方內(nèi)29個(gè)OTUs在群落中權(quán)重組成的物種矩陣和12種環(huán)境參數(shù)組成的環(huán)境變量矩陣為基礎(chǔ)數(shù)據(jù),應(yīng)用CCA排序,得到樣方、物種、環(huán)境因子相互關(guān)系的二維排序圖,并應(yīng)用TWINSPAN將29個(gè)OTUs分為7類(圖3)。排序結(jié)果表明,CCA排序圖第一軸AX1和第二軸AX2的特征值累計(jì)占總特征值的65.06%,排序圖包含了大部分的信息,排序效果良好。在CCA排序圖中,物種用三角形符號(hào)表示,物種之間夾角余弦值越大表示它們?cè)较嚓P(guān)。樣方點(diǎn)用圓圈表示,樣方之間的距離越近則表示樣方間越相似。各環(huán)境因子用箭頭表示,箭頭的長(zhǎng)短表示某個(gè)環(huán)境因子與種群分布相關(guān)性的大小,線越長(zhǎng),說明相關(guān)性越大,反之則小。從圖5箭頭的連線長(zhǎng)度可以看出,選取的12個(gè)與氮污染相關(guān)的水環(huán)境變量對(duì)蘆葦根結(jié)合好氣細(xì)菌群落結(jié)構(gòu)空間分布狀況都有一定程度的影響,其中,TDS、pH值、DO、Sal、TP、ORP對(duì)好氣細(xì)菌群落結(jié)構(gòu)時(shí)空分布影響較大,TN較NH+4、NO-3對(duì)群落結(jié)構(gòu)時(shí)空變異貢獻(xiàn)率高。箭頭方向表示環(huán)境變量之間或者環(huán)境變量與排序軸的之間的關(guān)系,夾角越小,相關(guān)性越大[16],其中NH+4、TN、ORP、TP與第1軸關(guān)系密切,NO-3、TDS、pH值、Sal與第2軸密切相關(guān)。

    圖3 29 OTUs的CCA雙軸排序圖,其中a為夏秋兩季物種與環(huán)境因子的CCA排序圖,b是對(duì)a圖中第4類的放大,是秋季特有物種與水質(zhì)的CCA排序圖Fig.3 Bi-plot(species and environmental variables)of CCA for summer and autumn survey on bacterial diversity in reclaimed water purification system(a),Bi-plot(species and environmental variables)of CCA for autumn survey on diversity of unique bacteria,which is the additional information of(b)

    依據(jù)TWINSPAN對(duì)物種分類結(jié)果,將圖3中29 OTUs圈為4類,其中夏季分為三類,分別是:如圖第1類,包含Ⅳ(SRW)樣方的偶見類群及部分非優(yōu)勢(shì)類群;第2類,包含Ⅳ、Ⅴ(SRW、SCW)樣方共有種群,且為優(yōu)勢(shì)種或非優(yōu)勢(shì)種;第3類,包含Ⅴ、Ⅵ(SCW、SOP)樣方的共有種及Ⅴ樣方的獨(dú)有種,秋季分為一類,即第4類,第4類包括秋季特有類群。圖3中第2類臨近排序中心,在凈化系統(tǒng)中具相對(duì)穩(wěn)定豐度,這些類群對(duì)環(huán)境適應(yīng)性較強(qiáng),生態(tài)幅比較廣泛,在濕地物質(zhì)和能量代謝中發(fā)揮重要作用,如,OTU28(B)在夏季3個(gè)功能區(qū)均出現(xiàn),為夏季再生水處理區(qū)優(yōu)勢(shì)種,是水質(zhì)凈化的重要類群,第1類遠(yuǎn)離排序中心,均為夏季再生水處理區(qū)獨(dú)有種,易受環(huán)境因子影響,較少頻率出現(xiàn)于再生水凈化系統(tǒng)中,對(duì)水體凈化不起主要作用。夏季種群沿第2軸分散分布于第1軸的負(fù)半軸,而秋季種群集中分布于第1軸的正半軸,細(xì)菌群落結(jié)構(gòu)具顯著季節(jié)差異。從圖3還可以看出,、TN、ORP、TP是導(dǎo)致凈化系統(tǒng)可培養(yǎng)細(xì)菌群落季節(jié)演替的主要水質(zhì)變量,、TDS、pH、Sal是導(dǎo)致夏季可培養(yǎng)細(xì)菌群落結(jié)構(gòu)空間演替的主要因素。第2類,包括優(yōu)勢(shì)類群OTU1(a)OTU28(B)與PCY和CHL密切相關(guān),對(duì)水體富營(yíng)養(yǎng)化去除其主要作用類群。第3類,包括樣方SOP優(yōu)勢(shì)類群OTU7(g)及樣方SCW偶見類群OTU3(c)、OTU10(j)、OTU11(k)、OTU16(p)與DO具有密切的關(guān)系,DO是改善濕地氧化還原環(huán)境的重要環(huán)境因子,氧化塘與復(fù)合垂直流濕地基質(zhì)中具不同氧氣分布特征,導(dǎo)致各區(qū)細(xì)菌群落結(jié)構(gòu)發(fā)生顯著改變。第1類與測(cè)定的水體富營(yíng)養(yǎng)化環(huán)境因子不具顯著關(guān)系,說明再生水特有的性質(zhì)影響了該區(qū)根結(jié)合細(xì)菌群落的偶見類群的分布格局,需要進(jìn)一步結(jié)合細(xì)菌測(cè)序結(jié)果分析。從圖3可以看出,秋季各樣點(diǎn)及蘆葦根結(jié)合可培養(yǎng)細(xì)菌群落結(jié)構(gòu)梯度排序是各水質(zhì)變量共同作用的結(jié)果,這與秋季各處理區(qū)凈化效率降低,各功能區(qū)優(yōu)勢(shì)類群趨于一致有緊密聯(lián)系。結(jié)合分類結(jié)果可以看出,第2類和第4類類群具相對(duì)較高豐度,受環(huán)境因子影響較小,如,O(r)是秋季凈化系統(tǒng)3個(gè)功能區(qū)共有優(yōu)勢(shì)種,OTU12(l)在秋季凈化系統(tǒng)3個(gè)功能區(qū)中均有出現(xiàn),但其與水體富營(yíng)養(yǎng)化環(huán)境因子不具密切關(guān)系,這說明秋季優(yōu)勢(shì)類群的水質(zhì)凈化作用較小。第1類較偏好高TN、TP濃度,第3類較偏好高DO、NH+4、Sal濃度,這兩類與水體富營(yíng)養(yǎng)化環(huán)境因子密切相關(guān),但這些類群不占優(yōu)勢(shì)。

    4 討論與結(jié)論

    大量研究表明微生物在人工濕地中起核心作用,研究濕地微生物多樣性的技術(shù)方法很多,而培養(yǎng)法只記錄適合在培養(yǎng)基上生長(zhǎng)的微生物群落,該法只能反映自然環(huán)境少部分的微生物,不能全面反映濕地中細(xì)菌群落的真實(shí)情況。但是,培養(yǎng)法的優(yōu)勢(shì)在于利用不同的培養(yǎng)基,可以選擇性的反映濕地中對(duì)不同資源具較強(qiáng)利用能力的微生物群落[20]。如本研究中利用LB培養(yǎng)基分離的細(xì)菌可以反映對(duì)再生水中氮素如蛋白質(zhì)、多肽、氨基酸、尿素等及有機(jī)碳具較強(qiáng)分解能力的細(xì)菌類群[28-29],是修復(fù)水體富營(yíng)養(yǎng)化的優(yōu)勢(shì)類群,分析其時(shí)空變異規(guī)律及其與水富營(yíng)養(yǎng)化因子的關(guān)系,對(duì)于再生水補(bǔ)水濕地水生態(tài)修復(fù)機(jī)制具有重要的意義。

    4.1 凈化系統(tǒng)可培養(yǎng)蘆葦根結(jié)合細(xì)菌群落特征及多樣性時(shí)空差異

    微生物與環(huán)境因子之間有密切的響應(yīng)關(guān)系,許多研究表明微生物群落結(jié)構(gòu)特征是微生物群落與環(huán)境因子長(zhǎng)期響應(yīng)的結(jié)果,受植物生長(zhǎng)狀況和代謝活動(dòng)、氣候季節(jié)變化、污水處理工藝、污水水體水質(zhì)狀況等多種因素影響,基于多因子共同作用微生物群落結(jié)構(gòu)及多樣性呈現(xiàn)差異[28,30]。

    (1)凈化系統(tǒng)蘆葦根結(jié)合好氣細(xì)菌群落密度、群落組成特征、群落多樣性均表現(xiàn)出顯著的季節(jié)差異。夏季較秋季凈化系統(tǒng)蘆葦根結(jié)合好氣細(xì)菌群落密度增加,這與夏季較秋季凈化系統(tǒng)細(xì)菌群落優(yōu)勢(shì)類群百分比增大,非優(yōu)勢(shì)類群及偶見類群百分比減少相關(guān)。夏季較秋季復(fù)合垂直流人工濕地細(xì)菌群落多樣性、豐富度、均勻度一致呈現(xiàn)增大趨勢(shì),這與夏季較秋季該類型人工濕地細(xì)菌群落優(yōu)勢(shì)類群、非優(yōu)勢(shì)類群和偶見類群種類數(shù)均增加密切相關(guān)。Florent Chazarenc等采用Biolog Ecoplate的方法對(duì)人工濕地微生物群落時(shí)空差異進(jìn)行分析,結(jié)果表明季節(jié)動(dòng)態(tài)是影響微生物群落多樣性最主要因素之一[15]。高程等對(duì)北京城區(qū)不同水體可培養(yǎng)細(xì)菌的研究表明夏季細(xì)菌數(shù)量高于秋季[20]。王曉丹等同時(shí)采用用LB固體培養(yǎng)基培養(yǎng)法和PCR-DGGE技術(shù)對(duì)密云水庫細(xì)菌群落結(jié)構(gòu)進(jìn)行研究,結(jié)果表明季節(jié)差異導(dǎo)致群落結(jié)構(gòu)發(fā)生改變[28]。陸開宏等通過對(duì)富營(yíng)養(yǎng)水體植物根際微生物群落特征分析發(fā)現(xiàn)因季節(jié)差異,植物的生長(zhǎng)狀況及代謝改變,植物根際微生物群落結(jié)構(gòu)出現(xiàn)差異,并影響水體污染物的去除效率[30]。由前人研究可知,在濕地生態(tài)系統(tǒng)或水體生態(tài)系統(tǒng)中,采用培養(yǎng)法或非培養(yǎng)法,細(xì)菌群落結(jié)構(gòu)研究結(jié)果均表明細(xì)菌群落受季節(jié)影響,具明顯的季節(jié)效應(yīng)。在再生水濕地凈化系統(tǒng)中,夏季向秋季的季節(jié)轉(zhuǎn)變不僅降低了凈化系統(tǒng)細(xì)菌生長(zhǎng)溫度,影響了細(xì)菌生長(zhǎng)繁殖及活性,而且導(dǎo)致濕地蘆葦枯萎凋零,植物功能減弱,根系分泌物質(zhì)能力、泌氧功能降低,植物-微生物的互作去污作用降低,依托于植物、微生物、基質(zhì)共同作用凈化水質(zhì)的人工濕地生態(tài)功能受到影響,水質(zhì)凈化效率降低,進(jìn)一步促使水體有機(jī)負(fù)荷及污染物的積累,水質(zhì)惡化,細(xì)菌密度及多樣性降低,群落結(jié)構(gòu)產(chǎn)生差異[31-32]。

    (2)凈化系統(tǒng)蘆葦根結(jié)合可培養(yǎng)細(xì)菌群落密度、細(xì)菌群落構(gòu)成特征、群落多樣性均表現(xiàn)出顯著的空間差異。夏季各處理區(qū)群落結(jié)構(gòu)差異明顯高于秋季,夏季凈化系統(tǒng)空間差異具體表現(xiàn)為復(fù)合垂直流濕地較氧化塘好氣細(xì)菌群落密度、多樣性、豐富度增加,均勻度減少,這與復(fù)合垂直流濕地非優(yōu)勢(shì)類群和偶見類群類型、數(shù)量百分比增加相關(guān);再生水處理區(qū)較循環(huán)水處理區(qū)可培養(yǎng)細(xì)菌群落密度、多樣性、豐富度、均勻度增加,這與再生水凈化區(qū)優(yōu)勢(shì)類群、非優(yōu)勢(shì)類群類型增加,非優(yōu)勢(shì)類群數(shù)量百分比增加,偶見類群量百分比減少相關(guān)。秋季凈化系統(tǒng)空間分布特征具體表現(xiàn)為:3個(gè)不同功能區(qū)具有相同優(yōu)勢(shì)類群,但3個(gè)功能區(qū)優(yōu)勢(shì)類群百分比呈現(xiàn)顯著差異,而偶見類群百分比在3各功能區(qū)具有趨同效應(yīng);復(fù)合垂直流濕地較氧化塘細(xì)菌群落多樣性及豐富度均減少,均勻度變化規(guī)律不明顯,這與復(fù)合垂直流濕地非優(yōu)勢(shì)類群、偶見類群類型減少相關(guān);再生水處理區(qū)較循環(huán)水處理區(qū)可培養(yǎng)細(xì)菌群落密度、多樣性、豐富度、均勻度增加,這與再生水處理區(qū)非優(yōu)勢(shì)類群類型數(shù)、百分比數(shù)均增加相關(guān),夏秋兩季分析結(jié)果一致表明再生水處理區(qū)具最高的群落密度和細(xì)菌群落多樣性。前人研究表明濕地微生物數(shù)量與其凈化功能之間存在顯著關(guān)系,微生物數(shù)量越多污染物的去除效率越高[27]。魏成等研究表明提高濕地系統(tǒng)的多樣性對(duì)提高濕地系統(tǒng)的凈化效率和效率的穩(wěn)定性具重要意義[14]。濕地工藝流程與其凈化功能之間同樣存在顯著關(guān)系,復(fù)合垂直流人工濕地特殊結(jié)構(gòu)及進(jìn)水類型有利于水體營(yíng)養(yǎng)及溶解氧分布、好養(yǎng)-厭氧-好養(yǎng)功能層形成,有利于微生物生長(zhǎng)繁殖及活性,促進(jìn)復(fù)雜有機(jī)物、氮、磷去除,提高污水凈化效率[33]。其次,進(jìn)水水質(zhì)的差異是影響微生物及其凈化功能的主要因素,再生水中豐富的氮、磷等營(yíng)養(yǎng)元素利于微生物數(shù)量增加,并改變微生物群落結(jié)構(gòu)及其多樣性[23]。本研究中,復(fù)合垂直流濕地先進(jìn)工藝特征及再生水特殊水質(zhì)是影響再生水補(bǔ)水濕地系統(tǒng)中不同處理區(qū)微生物群落結(jié)構(gòu)差異的主要原因。

    4.2 凈化系統(tǒng)蘆葦根結(jié)合可培養(yǎng)細(xì)菌群落結(jié)構(gòu)與水環(huán)境的相互關(guān)系

    據(jù)文獻(xiàn),細(xì)菌群落結(jié)構(gòu)與水質(zhì)之間的關(guān)系還沒有明確定論,研究方法也多采用相關(guān)性分析法。已有相關(guān)研究表明水質(zhì)狀況對(duì)微生物群落產(chǎn)生重要影響的同時(shí)微生物適應(yīng)環(huán)境并改變水體環(huán)境,微生物群落多樣性呈現(xiàn)出與水環(huán)境質(zhì)量協(xié)同變化的特征[34]。梁威等的研究發(fā)現(xiàn)濕地微生物活動(dòng)是污水中TN及COD的主要去除途徑[27],其隨后研究表明污水中BOD的去除率與復(fù)合垂直流濕地細(xì)菌總數(shù)顯著相關(guān)[35]。CCA和TWINSPAN是兩種具有很強(qiáng)解釋能力的多元統(tǒng)計(jì)方法,它們生態(tài)學(xué)中成功運(yùn)用,已被認(rèn)為是最有效的群落定量分析方法[36]。近年來有學(xué)者嘗試將CCA排序和多元回歸結(jié)合,這為分析微生物群落物種對(duì)水質(zhì)變量的聯(lián)合反應(yīng)提供了通用框架。YAN等采用CCA分析了污水處理廠微生物群落結(jié)構(gòu)的季節(jié)演替,結(jié)果表明COD、、KN、C/N/P值是重要的微生物生長(zhǎng)限制因子,DO是微生物群落演替的關(guān)鍵因子[16]。Shiah&Ducklow通過對(duì)細(xì)菌群落構(gòu)成的控制研究揭示出溫度是冬季、秋季、春季的細(xì)菌群落結(jié)構(gòu)的控制因子,而無機(jī)營(yíng)養(yǎng)限制及基質(zhì)是夏季細(xì)菌群落結(jié)構(gòu)的控制因子[37]。另一些研究,包括一個(gè)中型實(shí)驗(yàn),也支持同樣的觀點(diǎn),認(rèn)為有機(jī)碳、無機(jī)營(yíng)養(yǎng)是夏季細(xì)菌群落的重要限制因子[38-39]。微生物的呼吸作用和發(fā)酵過程是有機(jī)污染物的重要去除機(jī)制,其去除過程依賴于濕地氧化還原條件,高ORP為微生物提供氧化環(huán)境,有利于好養(yǎng)過程,如硝化作用、有機(jī)碳的降解,相反,低ORP提供還原環(huán)境,有利于厭氧過程,如硫酸鹽降解、甲烷化。

    本研究分析結(jié)果表明凈化系統(tǒng)蘆葦根結(jié)合好氣細(xì)菌群落密度與水體濃度及ORP相關(guān),群落構(gòu)成與水體PCY、TP、TN和濃度及ORP相關(guān),ORP是DO、有機(jī)物和含氮污染物三者氧化還原性的綜合體現(xiàn),寬泛的ORP范圍有利于多種污染物的去除[8],是好養(yǎng)微生物通過氨化作用對(duì)水體有機(jī)氮降解的產(chǎn)物。夏季凈化系統(tǒng)較秋季因具較高的細(xì)菌群落密度,群落多樣性,凈化效率提高,對(duì)有機(jī)物的去除效率提高,雖然夏季較秋季水體具較高的DO濃度,但較秋季ORP依然處于較低水平,這也是ORP與可培養(yǎng)細(xì)菌群落密度負(fù)相關(guān)的原因,較高的DO濃度有利于細(xì)菌的氨化作用,這是夏季較秋季水體具低水平COD、TDS及TN,同時(shí)與好氣細(xì)菌群落密度正相關(guān)的原因。PCY濃度是水體富營(yíng)養(yǎng)化的重要水質(zhì)指標(biāo),其生長(zhǎng)繁殖需要充足的營(yíng)養(yǎng)和優(yōu)越的物理?xiàng)l件[40],PCY濃度與優(yōu)勢(shì)種百分比相關(guān),說明優(yōu)勢(shì)種對(duì)消除水體富營(yíng)養(yǎng)化,對(duì)去除氮污染起主要作用,ORP及氮磷含量及N/P比值改變導(dǎo)致偶見種類型改變,改變好氣細(xì)菌群落種群豐富度。因此ORP、濃度既影響好氣微生物群落數(shù)量,同時(shí)改變微生物群落結(jié)構(gòu)。由CCA分析表明、COD、ORP、TP是蘆葦可培養(yǎng)根結(jié)合細(xì)菌群落結(jié)構(gòu)季節(jié)差異的驅(qū)動(dòng)因子。在驅(qū)動(dòng)因子中,、ORP、TP與好氣細(xì)菌偶見種百分比相關(guān),說明該濕地凈化系統(tǒng)細(xì)菌群落季節(jié)差異不僅體現(xiàn)在優(yōu)勢(shì)種群的改變上,同時(shí)體現(xiàn)在偶見種群的多少上,偶見種群的百分比與濕地凈化系統(tǒng)的穩(wěn)定性及凈化效率密切相關(guān)。而這種微生物群落群落結(jié)構(gòu)季節(jié)差異很大程度上可以由夏秋植物生長(zhǎng)狀況及代謝活動(dòng)不同,微生物生長(zhǎng)繁殖及活性不同來解釋。在夏季,蘆葦生長(zhǎng)狀況良好,根系的泌氧功能有利于根際氧化還原電位提高,碳(COD)、氮()、磷(TP)及高氧化還原電位促進(jìn)微生物生長(zhǎng)繁殖及活性,細(xì)菌群落度增加,促進(jìn)水體凈化,濕地系統(tǒng)凈化效率提高。在夏季濕地較高凈化效率的基礎(chǔ)之上,功能區(qū)間水質(zhì)差異增大,微生物群落結(jié)構(gòu)空間差異增大,在排序圖中夏季種群分散分布;相反,在秋季,微生物活動(dòng)減弱,植物枯萎,水體碳氮磷濃度增加,氧化還原電位降低,高有機(jī)負(fù)荷抑制微生物生長(zhǎng)及活性[31-32],細(xì)菌群落密度減少,濕地凈化系統(tǒng)凈化效率降低,好氣細(xì)菌群落結(jié)構(gòu)空間差異不明顯,在排序圖中秋季種群集中分布。本研究結(jié)果表明水質(zhì)變量TDS、pH值、Sal是導(dǎo)致夏季好氣細(xì)菌群落結(jié)構(gòu)空間差異的主要水質(zhì)變量。SRW較SCW和SOP具較低TDS、pH值、Sal水平,在SCW和SOP細(xì)菌群落更偏好于較高TDS、Sal濃度及較高pH值。而在SRW,細(xì)菌群落與富營(yíng)養(yǎng)化程度相關(guān)。這可能與濕地功能分布相關(guān)。

    據(jù)酶切分類結(jié)果,分別于29 OTUs中隨機(jī)選取分離物進(jìn)行16s rDNA序列測(cè)定,同源性比對(duì)表明:夏季凈化系統(tǒng)中優(yōu)勢(shì)好氣細(xì)菌為葡萄球?qū)?Staphylococcus)、腸桿菌屬(Enterobacteriaceae),微小桿菌屬(Exiguobacterium)和芽孢桿菌屬(Bacillus.),氣單胞菌屬(Aeromonas),其中氣單胞菌屬分布于夏季凈化系統(tǒng)各功能區(qū)中。秋季人工濕地凈化系統(tǒng)主要好氣細(xì)菌為假單胞菌屬(Pseudomonas),黃單胞菌屬(Xanthomonadaceae),芽孢桿菌屬(Bacillus),其中假單胞菌屬為3個(gè)處理區(qū)共有的優(yōu)勢(shì)種。據(jù)已有報(bào)道稱,氣單胞菌在濕地系統(tǒng)氮循環(huán)過程中起重要作用,其生長(zhǎng)受銨離子、氮、磷濃度影響明顯[41],具有將轉(zhuǎn)化為的能力[42],腸桿菌耐污能力強(qiáng),具固氮能力,并抗植物病原體,同時(shí)在磷循環(huán)中起重要作用[10]。葡萄球?qū)偌?xì)菌多為機(jī)會(huì)致病病原菌,對(duì)固氮菌有協(xié)同促進(jìn)作用[43]。秋季主要細(xì)菌種類:假單胞菌屬,芽孢桿菌屬,及黃單胞菌屬,這幾種菌屬的微生物生長(zhǎng)快速,其中假單胞菌被頻繁的從濕地植物根系周圍分離到,不僅以碳?xì)浠衔餅闋I(yíng)養(yǎng)來源,在去除碳源污染方面起主要作用,還是氨化細(xì)菌的優(yōu)勢(shì)菌屬,具很好地去除有機(jī)氮的效果,芽孢桿菌對(duì)大分子結(jié)構(gòu)的有機(jī)物具破壞作用,有效促進(jìn)有機(jī)物的分解[44]??傊J葦根結(jié)合菌促進(jìn)了水體凈化過程中的植物修復(fù)作用,尤其在氮磷循環(huán)及碳化合物降解去除方面。但是,限于培養(yǎng)法不能直接的反映蘆葦根結(jié)合菌的全部信息,結(jié)合非培養(yǎng)方法分析奧林匹克公園再生水凈化濕地蘆葦根結(jié)合菌對(duì)水體凈化的機(jī)制原理是下一步的工作重點(diǎn)。

    [1]Zhou L B,Han L B,Su D R.Advance in research on application of reclaimed water for turf-grass irrigation.Water Saving Irrigation.Water Saving Irrigation,2004,5(1):17-20.

    [2]Guo X Y,Gong H L,Zhao W J.Effects of recycled water irrigation on bacterial community compositions in lawn soil.Acta Microbiologica Sinica,2010,50(11):1519-1524.

    [3]Yang H,Abbaspour K C.Analysis of wastewater reuse potential in Beijing.Desalination,2007,212(1/3):238-250.

    [4]Wen J.Beijing water authority issued the 12thFive-Year Plan about reclaimed water.China News Service,2011-01-28.http://www.chinanews.com/ny/2011/01-28/2817700.shtml.

    [5]Zhu J N,Wang X D,Zhai Z H,Ma W L,Li R Q,Wang X L,Li Y H.Changes of bacterial community structure on reusing domestic sewage of daoxianghujing hotel to landscape water.Environmental Science,2010,31(5):1299-1306.

    [6]Zhao SM,Hu N,Chen Z J,Zhao B,Liang Y X.Bioremediation of reclaimed wastewater used as landscape water by using the denitrifying bacterium Bacillus cereus.Bulletin of Environmental Contamination and Toxicology,2009,83(3):337-340.

    [7]Thurston JA,Gerba C P,F(xiàn)oster K E,Karpiscak M M,Gerba C P.Fate of indicator microorganisms,giardia and cryptosporidium in subsurface flow constructed wetlands.Water Research,2001,35(6):1547-1551.

    [8]Faulwetter J L,Gagono V,Sundberg C,Chazarenc F,Burr M D,Brisson J,Camper A K,Stein OR.Microbial processes influencing performance of treatment wetlands:a review.Ecological Engineering,2009,35(6):987-1004.

    [9]Kennedy A C.Rhizosphere//Sylvia D M,Hartel P G,F(xiàn)uhrman J J,and Zuberer Z A,eds.Principals and Applications of Soil Microbiology.Upper Saddle River:Pearson Prentice Hall Publishing,2005:309-407.

    [10]Li Y H,Zhu JN,Zhai Z H.Endophytic bacterial diversity in roots of Phragmites australisin constructed Beijing Cuihu Wetland(China).FEMS Microbiology Letters,2010,309(1):84-93.

    [11]Li Y H,Liu Y F,Liu Y,Zhu J N,Zhang Q.Endophytic bacterial diversity in roots of Typha angustifolia L.in the constructed Beijing Cuihu Wetland(China).Research in Microbiology,2011,162(2):124-131.

    [12]Mustafa A,Scholz M.Characterization of microbial communities transforming and removing nitrogen in wetlands.Wetlands,2011,31(3):583-592.

    [13]Ravit B,Ehrenfeld J G,Haeggblom M M,Bartels M.The effects of drainage and nitrogen enrichment on Phragmites australis,Spartina alterniflora,and their root-associated microbial communities.Wetlands,2007,27(4):915-927.

    [14]Wei C,Liu P.Relationship between wastewater purification and diversity of rhizosphere microorganism in the constructed wetland.Journal of Agro-Environment Science,2008,27(6):2401-2406.

    [15]Chazarenc F,Brisson J,Merlin G.Seasonal and spatial changes of microorganism communities in constructed wetlands:a community level physiological profiling analysis.International Journal of Chemical Engineering,2010,doi:10.1155/2010/490240.

    [16]Yan Q M,Zhang X X,Zhang T,F(xiàn)ang H H P.Seasonal microbial community shift in a saline sewage treatment plant.Frontiers of Environmental Science and Engineering in China,2011,5(1):40-47.

    [17]Kirk J L,Beaudette L A,Hart M,Moutoglis P,Klironomos J N,Lee H,Trevors J T.Methods of studying soil microbial diversity.Journal of Microbiological Methods,2004,58(2):169-188.

    [18]Bohannan B JM,Hughes J.New approaches to analyzing microbial biodiversity data.Current Opinion in Microbiology,2003,6(3):282-287.

    [19]Hughes J B,Hellmann J J,Ricketts T H,Bohannnan B J.Counting the uncountable:statistical approaches to estimating microbial diversity.Applied and Environmental Microbiology,2001,67(10):4399-4406.

    [20]Gao C,Huang M R,Tao S,Sun X,Li J,Guo L D.Seasonal dynamics of culturable bacterium numbers in freshwater bodies of different water quality in Beijing.Acta Ecologica Sinica,2011,31(4):1157-1163.

    [21]Hejcmanovā-NezerkováP,Hejcman M.A canonical correspondence analysis(CCA)of the vegetation-environment relationships in Sudanese savannah,Senegal.South African Journal of Botany,2006,72(2):256-262.

    [22]Deng Z Z.On Water System and Rainwater Harvesting System in Beijing Olympic Green.Beijing:China Water Power Press,2009.

    [23]Guo X Y,Dong Z,Gong H L.Influence of reclaimed water irrigation on microbial community on lawn soil.China Environmental Science,2006,26(4):482-485

    [24]Sakai M,Matsunaga N,Matsuka A,Kanazawa S.Application of T-RFLPanalisis to the study of bacterial community structure in the rhizosphere.Soil Science and Plant Nutrition,2001,47(4):773-778.

    [25]Dong X Z.Common Bacteria Identification Manual.Beijing:Chinese TCM Publishing House,2000.

    [26]Sun L.Endophytic Bacteria and Root-Associated Bacteria with Rice by Culture-Independent and Culture-Dependent Approaches[D].Beijing:Capital Normal University,2006.

    [27]Liang W,Wu Z B,Zhou Q H,Cheng SP,F(xiàn)u G P.Analysis of substrate microorganisms in the constructed wetland and their correlation with wastewater purification effects.China Environmental Science,20r02,22(3):282-285.

    [28]Wang X D,Zhai Z H,Zhao S.The bacterial diversity in different seasons of Miyun Resevoir.Acta Ecologica Sinica,2009,29(7):3919-3927.

    [29]Li F D,Yu Z N,He SH.Agricultural Microbiology Experiment Technology.Beijing:China Agriculture Press,1996.

    [30]Lu K H,Hu Z Y,Liang JJ,Zhu JY.Characteristics of rhizosphere microbial community structure of two aquatic plants in eutrophic waters.China Environmental Science,2010,30(11):1508-1515.

    [31]Huang D F,Li T,Lu B.Pollutants removal and analysis of structure changes of microbial community in integrated vertical-flow constructed wetland.Research of Environmental Science,2007,20(6):137-141.

    [32]Wang X D,Zhai Z H,Zhao S,Li R Q,Ma W L,Li Y H.Effect of free surface flow wetland and subsurface flow wetland on bacterial diversity in Beijing Cuihu Wetland Park.Environmental Science,2009,30(1):280-288.

    [33]Zhou Q H,He F,Zhang L P,Wang Y F,Wu Z B.Characteristics of the microbial communities in the integrated vertical-flow constructed wetlands.Journal of Environmental Sciences,2009,21(9):1261-1267.

    [34]Wang Z Q,Wang W D,Zhu G B,Yin C Q.A comparative study on the diversity of rhizospheric bacteria community structure in constructed wetland and natural wetland with reed domination.Acta Ecologica Sinica,2011,31(16):4489-4498.

    [35]Liang W,Wu Z B,Zhan F C,Deng J Q.Seasonal variation of macrophytes root-zone microorganisms and purification effect in the constructed wetland system.Journal of Lake Sciences,2004,16(4):312-317.

    [36]ter Braak C J F,Verdonschot P F M.Canonical correspondence analysis and related multivariate methods in aquatic ecology.Aquatic Sciences,1995,57(3):254-289.

    [37]Shiah F K,Ducklow H W.Temperature regulation of heterotrophic bacterioplankton abundance,production,and specific growth rate in Chesapeake Bay.Limnology Oceanography,1994,39(6):1243-1258.

    [38]Rivkin R B,Anderson M R.Inorganic nutrient limitation of oceanic bacterioplankton.Limnol Oceanogr,1997,42(4):730-740.

    [39] ?vre?s L,Bourne D,Sandaa R A,Casamayor E O,Benlloch S,Goddard V,Smerdon G,Heldal M,Thingstad T F.Response of bacterial and viral communities to nutrient manipulations in seawater mesocosms.Aquatic Microbial Ecology,1997,31(2):109-121..

    [40]Kong F X,Gao G.Hypothesis on cyanobacteria bloom-forming mechanism in large shallow eutrophic lakes,Acta Ecologica Sinica,2005,25(3):589-595.

    [41]Minana-Galbis D,F(xiàn)arfan M,F(xiàn)uste M C,Lore J G.Aeromonas bivalvium sp.nov.,isolated from bivalve mollusks.Int.J.Syst.International Journal of Systematic and Evolutionary,2007,57(3):582-587.

    [42]Pianetti A,Bruscolini F,Rocchi M B,Sabatini L,Citterio B.Influence of different concentrations of nitrogen and phosphorous on Aeromonas spp.growth.Igiene E SanitàPubblica,2006,62(6):609-622.

    [43]Holguin G,Guzman M A,Bashan Y.Two new nitrogen-fixing bacteria from the rhizosphere of mangrove trees:their isolation,identification and in vitro interaction with rhizosphere Staphylococcus sp..FEMSMicrobiology Ecology Letters,1992,101(3):207-216.

    [44]Chen B Q,Yin C Q.The ecological distributions of bacteria and algae in the stimulation wetland system.Acta Ecologica Sinica,1998,18(6):634-639.

    參考文獻(xiàn):

    [1] 周陸波,韓烈保.再生水灌溉草坪綠地的研究進(jìn)展.節(jié)水灌溉,2004,5:17-20.

    [2] 郭逍宇,宮輝力,趙文吉.再生水灌溉對(duì)草坪根際可培養(yǎng)細(xì)菌組成影響.微生物學(xué)報(bào),2010,50(11):1519-1524.

    [4] 文靜.北京水務(wù)局發(fā)布再生水十二五規(guī)劃.中國(guó)新聞網(wǎng),2011-01-28.http://www.chinanews.com/ny/2011/01-28/2817700.shtml.

    [5] 朱競(jìng)男,王曉丹,翟振華,馬文林,李榮旗,王雪蓮,李艷紅.稻香湖景酒店景觀再生水生產(chǎn)中的細(xì)菌群落結(jié)構(gòu)變化.環(huán)境科學(xué),2010,31(5):1299-1306.

    [14] 魏成,劉平.人工濕地污水凈化效率與根際微生物群落多樣性的相關(guān)性研究.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2008,27(6):2401-2406.

    [20] 高程,黃滿榮,陶爽,孫翔,黎景,郭良棟.北京城區(qū)不同水質(zhì)水體可培養(yǎng)細(xì)菌數(shù)量的季節(jié)動(dòng)態(tài)變化.生態(tài)學(xué)報(bào),2011,31(4):1157-1163.

    [22] 鄧卓智.北京奧林匹克公園水系及雨洪利用系統(tǒng)研究、設(shè)計(jì)及示范.北京:中國(guó)水利水電出版社,2009.

    [23] 郭逍宇,董志,宮輝力.再生水灌溉對(duì)草坪土壤微生物群落的影響.中國(guó)環(huán)境科學(xué),2006,26(4):482-485.

    [25] 東秀珠.常見細(xì)菌系統(tǒng)鑒定手冊(cè).北京:中國(guó)中醫(yī)藥出版社,2000.

    [26] 孫磊,非培養(yǎng)方法和培養(yǎng)方法對(duì)水稻內(nèi)生細(xì)菌和根結(jié)合細(xì)菌的研究[D].北京:首都師范大學(xué),2006.

    [27] 梁威,吳振斌,周巧紅,成水平,付貴萍.構(gòu)建濕地基質(zhì)微生物與凈化效果及相關(guān)分析.中國(guó)環(huán)境科學(xué),2002,22(3):282-285.

    [28] 王曉丹,翟振華,趙爽.密云水庫不同季節(jié)細(xì)菌群落多樣性.生態(tài)學(xué)報(bào),2009,29(7):3919-3927.

    [29] 李阜棣,喻子牛,何紹紅.農(nóng)業(yè)微生物學(xué)實(shí)驗(yàn)技術(shù).北京:中國(guó)農(nóng)業(yè)出版社,1996.

    [30] 陸開宏,胡智勇,梁晶晶,朱津永.富營(yíng)養(yǎng)水體中2種水生植物的根際微生物群落特征.中國(guó)環(huán)境科學(xué),2010,30(11):1508-1515.

    [31] 黃德鋒,李田,陸斌.復(fù)合垂直流人工濕地污染物去除及微生物群落結(jié)構(gòu)的PCR-DGGE分析.環(huán)境科學(xué)研究,2007,20(6):137-141.

    [32] 王曉丹,翟振華,趙爽,李榮旗,馬文林,李艷紅.北京翠湖表流和潛流濕地對(duì)細(xì)菌多樣性的影響.環(huán)境科學(xué),2009,30(1):280-288.

    [34] 汪仲瓊,王為東,祝貴兵,尹澄清.人工和天然濕地蘆葦根際土壤細(xì)菌群落結(jié)構(gòu)多樣性的比較.生態(tài)學(xué)報(bào),2011,31(16):4489-4498.

    [35] 梁威,吳振斌,詹發(fā)萃,鄧家齊.人工濕地植物根區(qū)微生物與凈化效果的季節(jié)變化.湖泊科學(xué),2004,16(4):312-317.

    猜你喜歡
    凈化系統(tǒng)類群蘆葦
    石磨豆腐
    蘆葦
    黃河之聲(2021年19期)2021-02-24 03:25:24
    薏苡種質(zhì)資源ISSR分子標(biāo)記篩選及親緣關(guān)系分析
    汽車空氣質(zhì)量智能檢測(cè)凈化系統(tǒng)
    基于單片機(jī)的室內(nèi)空氣凈化系統(tǒng)
    電子制作(2018年9期)2018-08-04 03:30:52
    蘆葦
    歲月(2018年2期)2018-02-28 20:40:58
    黑石頂自然保護(hù)區(qū)土壤動(dòng)物功能類群對(duì)季節(jié)動(dòng)態(tài)的響應(yīng)
    蘆葦筏
    幼兒100(2016年28期)2016-02-28 21:26:18
    高爐煤氣凈化系統(tǒng)發(fā)展前沿技術(shù)簡(jiǎn)況
    2011年春夏季黃海和東海微型浮游動(dòng)物類群組成及其攝食的研究
    亚洲av五月六月丁香网| 久久久久久伊人网av| 床上黄色一级片| 欧美成人一区二区免费高清观看| 亚洲男人的天堂狠狠| 亚洲午夜理论影院| 自拍偷自拍亚洲精品老妇| 国产蜜桃级精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 嫁个100分男人电影在线观看| 极品教师在线视频| 国内精品一区二区在线观看| 午夜a级毛片| 欧美日韩中文字幕国产精品一区二区三区| 久久久久九九精品影院| 男女那种视频在线观看| 全区人妻精品视频| 欧美黑人欧美精品刺激| 精品欧美国产一区二区三| 色av中文字幕| 婷婷精品国产亚洲av| 亚洲成人中文字幕在线播放| 精品免费久久久久久久清纯| 亚洲av不卡在线观看| 亚洲成av人片在线播放无| 欧美一区二区精品小视频在线| 久久久久久久久久黄片| 麻豆成人午夜福利视频| 99精品在免费线老司机午夜| 日韩一本色道免费dvd| 亚洲熟妇中文字幕五十中出| 麻豆国产av国片精品| 久久人妻av系列| 黄色配什么色好看| 精华霜和精华液先用哪个| 99久久精品热视频| 美女免费视频网站| 日韩精品青青久久久久久| 最新中文字幕久久久久| 国产蜜桃级精品一区二区三区| 国产高清激情床上av| 自拍偷自拍亚洲精品老妇| 一区二区三区高清视频在线| 日本五十路高清| 亚洲一区二区三区色噜噜| 国产精品亚洲一级av第二区| 亚洲国产精品合色在线| 亚洲第一区二区三区不卡| 亚洲中文字幕日韩| 男女啪啪激烈高潮av片| 婷婷色综合大香蕉| 精品乱码久久久久久99久播| 国产精品女同一区二区软件 | 亚洲精品影视一区二区三区av| 久久九九热精品免费| 日本撒尿小便嘘嘘汇集6| 老司机福利观看| 99热精品在线国产| 在线播放无遮挡| 一级黄色大片毛片| 成人毛片a级毛片在线播放| 日本成人三级电影网站| 免费无遮挡裸体视频| 国产精品综合久久久久久久免费| 搡老岳熟女国产| 男女边吃奶边做爰视频| 日本爱情动作片www.在线观看 | 亚洲成人免费电影在线观看| 真人做人爱边吃奶动态| 国产精品久久久久久久久免| 一进一出抽搐动态| 精品人妻熟女av久视频| 91麻豆av在线| 真实男女啪啪啪动态图| 日韩在线高清观看一区二区三区 | 草草在线视频免费看| 国产精品人妻久久久影院| 欧美性猛交╳xxx乱大交人| 日韩欧美精品免费久久| 成年免费大片在线观看| 在线播放国产精品三级| 国产一级毛片七仙女欲春2| 在线免费观看的www视频| 岛国在线免费视频观看| av女优亚洲男人天堂| 狂野欧美激情性xxxx在线观看| 熟妇人妻久久中文字幕3abv| 久久久久国产精品人妻aⅴ院| 麻豆久久精品国产亚洲av| 精品人妻视频免费看| 99久久中文字幕三级久久日本| 精品午夜福利在线看| 欧美成人性av电影在线观看| 国产麻豆成人av免费视频| 久久精品国产亚洲av天美| 色综合婷婷激情| 不卡一级毛片| 国产高清激情床上av| 极品教师在线免费播放| 久久午夜亚洲精品久久| 中文字幕免费在线视频6| 波多野结衣巨乳人妻| 欧美三级亚洲精品| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 国产真实乱freesex| 色尼玛亚洲综合影院| 搡老熟女国产l中国老女人| 在线天堂最新版资源| 国产一区二区在线观看日韩| 国产精品爽爽va在线观看网站| 天天躁日日操中文字幕| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 久久久久久久久久成人| 永久网站在线| 亚洲美女视频黄频| 日韩欧美免费精品| 国产人妻一区二区三区在| 亚洲在线自拍视频| 午夜日韩欧美国产| 人妻夜夜爽99麻豆av| 草草在线视频免费看| 观看美女的网站| 欧美xxxx性猛交bbbb| 国产aⅴ精品一区二区三区波| 亚洲va日本ⅴa欧美va伊人久久| 亚洲 国产 在线| 国产精品人妻久久久影院| 欧美成人a在线观看| 婷婷丁香在线五月| 成人综合一区亚洲| 国产亚洲精品久久久久久毛片| 国产三级中文精品| 最近视频中文字幕2019在线8| 亚洲五月天丁香| 午夜激情欧美在线| 一级av片app| 一卡2卡三卡四卡精品乱码亚洲| 乱系列少妇在线播放| 免费黄网站久久成人精品| av在线亚洲专区| 国产免费一级a男人的天堂| av天堂在线播放| 久久中文看片网| 国产高清有码在线观看视频| 国产中年淑女户外野战色| 日本一本二区三区精品| 成人av一区二区三区在线看| 天美传媒精品一区二区| 中文字幕人妻熟人妻熟丝袜美| 国产高清视频在线播放一区| 我要搜黄色片| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 天天躁日日操中文字幕| 免费在线观看日本一区| 国产一区二区在线观看日韩| 亚洲av免费高清在线观看| 成人一区二区视频在线观看| 国产一区二区三区av在线 | 热99在线观看视频| 日韩亚洲欧美综合| 国产精品99久久久久久久久| 色吧在线观看| xxxwww97欧美| 小说图片视频综合网站| 日韩欧美在线二视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一及| 亚洲精品456在线播放app | 中文字幕精品亚洲无线码一区| 成人永久免费在线观看视频| 亚洲一级一片aⅴ在线观看| 成年人黄色毛片网站| 在线观看舔阴道视频| www日本黄色视频网| 18禁裸乳无遮挡免费网站照片| 联通29元200g的流量卡| 村上凉子中文字幕在线| 俺也久久电影网| 国产亚洲精品综合一区在线观看| 国产一级毛片七仙女欲春2| 久久久精品欧美日韩精品| 国产精品久久久久久久久免| 亚洲成人免费电影在线观看| 最近中文字幕高清免费大全6 | 婷婷精品国产亚洲av在线| 狠狠狠狠99中文字幕| 乱码一卡2卡4卡精品| 国内少妇人妻偷人精品xxx网站| 极品教师在线免费播放| 中文字幕久久专区| 国产蜜桃级精品一区二区三区| 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 免费av毛片视频| 51国产日韩欧美| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 中出人妻视频一区二区| 亚洲专区中文字幕在线| 日本一二三区视频观看| 亚洲色图av天堂| 香蕉av资源在线| 乱系列少妇在线播放| 成人国产综合亚洲| 国产国拍精品亚洲av在线观看| 日本黄色片子视频| ponron亚洲| 91在线精品国自产拍蜜月| 精品久久久噜噜| 色在线成人网| 最近在线观看免费完整版| 在线观看舔阴道视频| 国产主播在线观看一区二区| 最近最新免费中文字幕在线| 国产在线精品亚洲第一网站| 国产午夜福利久久久久久| 欧美bdsm另类| 亚洲无线在线观看| 在线免费十八禁| 国产色爽女视频免费观看| 小说图片视频综合网站| 国产精品,欧美在线| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 国产精品98久久久久久宅男小说| 久久精品夜夜夜夜夜久久蜜豆| 免费人成视频x8x8入口观看| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 一区福利在线观看| 色播亚洲综合网| 亚洲美女搞黄在线观看 | 久久香蕉精品热| 国产精品国产三级国产av玫瑰| 亚洲四区av| 国产在线男女| 深夜精品福利| 国产久久久一区二区三区| 亚洲成人免费电影在线观看| 九色成人免费人妻av| 亚洲电影在线观看av| 欧美日本亚洲视频在线播放| 69人妻影院| 乱系列少妇在线播放| 国产久久久一区二区三区| 欧美日韩瑟瑟在线播放| 国产三级中文精品| 在线免费观看的www视频| 观看美女的网站| 久久精品国产亚洲网站| 日韩欧美在线乱码| av中文乱码字幕在线| 在线观看美女被高潮喷水网站| 极品教师在线视频| 在线观看66精品国产| 亚洲av成人av| 国产色爽女视频免费观看| 婷婷色综合大香蕉| 国产69精品久久久久777片| 在线观看免费视频日本深夜| 久久久久久久久久黄片| 亚洲va在线va天堂va国产| 欧美日韩乱码在线| 亚洲av日韩精品久久久久久密| 国产v大片淫在线免费观看| 欧美最新免费一区二区三区| 色哟哟·www| 久久精品国产99精品国产亚洲性色| 特级一级黄色大片| 99久久精品热视频| 国产成人aa在线观看| 亚洲无线观看免费| 他把我摸到了高潮在线观看| 能在线免费观看的黄片| 国产伦一二天堂av在线观看| 99久久中文字幕三级久久日本| 精品无人区乱码1区二区| 伦理电影大哥的女人| 亚洲内射少妇av| 欧美成人免费av一区二区三区| 12—13女人毛片做爰片一| 天美传媒精品一区二区| 亚洲精品乱码久久久v下载方式| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品粉嫩美女一区| 很黄的视频免费| 国产成人一区二区在线| 午夜福利欧美成人| 日韩欧美免费精品| 久久精品国产鲁丝片午夜精品 | a级毛片a级免费在线| 啦啦啦啦在线视频资源| 少妇高潮的动态图| 日本爱情动作片www.在线观看 | 在线观看av片永久免费下载| 国产亚洲av嫩草精品影院| 欧美精品啪啪一区二区三区| 欧美性感艳星| 人人妻人人看人人澡| 长腿黑丝高跟| 中文在线观看免费www的网站| 久久国产乱子免费精品| 欧美丝袜亚洲另类 | 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 欧美黑人巨大hd| 伦精品一区二区三区| 欧美区成人在线视频| 免费在线观看日本一区| 三级男女做爰猛烈吃奶摸视频| 色吧在线观看| 性色avwww在线观看| 又黄又爽又刺激的免费视频.| 人妻夜夜爽99麻豆av| 国产精品国产三级国产av玫瑰| 性插视频无遮挡在线免费观看| 成熟少妇高潮喷水视频| 久久精品91蜜桃| 国产免费av片在线观看野外av| 天天躁日日操中文字幕| 男人狂女人下面高潮的视频| 熟妇人妻久久中文字幕3abv| 五月伊人婷婷丁香| 国产精品自产拍在线观看55亚洲| 国产69精品久久久久777片| 日韩精品中文字幕看吧| 久久欧美精品欧美久久欧美| 亚洲国产欧洲综合997久久,| 国产精品久久久久久av不卡| 亚洲第一电影网av| 成人国产一区最新在线观看| 中国美女看黄片| 亚洲成av人片在线播放无| 九九久久精品国产亚洲av麻豆| 人妻少妇偷人精品九色| 在线观看舔阴道视频| 亚洲精品成人久久久久久| 男插女下体视频免费在线播放| 国产av麻豆久久久久久久| 一a级毛片在线观看| 午夜精品久久久久久毛片777| 欧美人与善性xxx| 久久精品久久久久久噜噜老黄 | 日韩欧美精品免费久久| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 十八禁网站免费在线| 干丝袜人妻中文字幕| 999久久久精品免费观看国产| 最近视频中文字幕2019在线8| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 亚洲av第一区精品v没综合| 窝窝影院91人妻| 草草在线视频免费看| 最后的刺客免费高清国语| 91在线精品国自产拍蜜月| 99久久成人亚洲精品观看| 亚洲精品一区av在线观看| 岛国在线免费视频观看| 国产 一区 欧美 日韩| 国产一区二区三区视频了| 亚洲av免费在线观看| 欧美色欧美亚洲另类二区| 久久久久国产精品人妻aⅴ院| 欧美精品国产亚洲| 午夜福利成人在线免费观看| ponron亚洲| 亚洲精品影视一区二区三区av| 久久久久久久久中文| 国产私拍福利视频在线观看| 国产麻豆成人av免费视频| 在线观看av片永久免费下载| 亚洲欧美激情综合另类| 2021天堂中文幕一二区在线观| 三级毛片av免费| 欧美+亚洲+日韩+国产| 国产中年淑女户外野战色| 午夜福利在线在线| 最新在线观看一区二区三区| 国产一区二区在线观看日韩| 亚洲第一电影网av| 亚洲自偷自拍三级| 99热精品在线国产| 舔av片在线| 欧美高清成人免费视频www| 欧美性猛交黑人性爽| 免费一级毛片在线播放高清视频| 久久久久久大精品| 国产免费男女视频| 久久久久久久精品吃奶| 三级男女做爰猛烈吃奶摸视频| 亚洲精品久久国产高清桃花| 免费观看在线日韩| 麻豆国产av国片精品| 精品人妻偷拍中文字幕| 麻豆av噜噜一区二区三区| 人人妻,人人澡人人爽秒播| 国产成人av教育| 午夜福利在线观看免费完整高清在 | 在线观看午夜福利视频| 波野结衣二区三区在线| 国产精品久久久久久av不卡| .国产精品久久| 国产精品久久久久久av不卡| 变态另类成人亚洲欧美熟女| 最后的刺客免费高清国语| 免费av不卡在线播放| 91久久精品国产一区二区成人| 看免费成人av毛片| 亚洲aⅴ乱码一区二区在线播放| 中文字幕高清在线视频| 亚洲av日韩精品久久久久久密| av专区在线播放| 十八禁网站免费在线| 精品久久久噜噜| 国产精品人妻久久久影院| 久久久久久久亚洲中文字幕| 亚洲av成人av| 看片在线看免费视频| 欧美精品国产亚洲| 窝窝影院91人妻| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 长腿黑丝高跟| 看片在线看免费视频| 欧美一区二区国产精品久久精品| 日韩欧美在线乱码| 免费大片18禁| 国内精品一区二区在线观看| 欧美成人a在线观看| 黄色日韩在线| 亚洲欧美清纯卡通| 午夜福利高清视频| 亚洲成人久久爱视频| 婷婷精品国产亚洲av在线| 最新在线观看一区二区三区| 久久国产乱子免费精品| 精品99又大又爽又粗少妇毛片 | 国产久久久一区二区三区| 国产伦精品一区二区三区视频9| 久久这里只有精品中国| 久久久久久久久久久丰满 | 在线观看午夜福利视频| 日韩中文字幕欧美一区二区| 久久九九热精品免费| 成人毛片a级毛片在线播放| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 一级黄片播放器| 亚洲国产精品sss在线观看| 日本在线视频免费播放| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| eeuss影院久久| 国产aⅴ精品一区二区三区波| 精品久久国产蜜桃| 伦精品一区二区三区| 日本 欧美在线| av在线蜜桃| 午夜激情福利司机影院| 久久精品久久久久久噜噜老黄 | 欧美一区二区亚洲| 变态另类成人亚洲欧美熟女| 久久九九热精品免费| 此物有八面人人有两片| 免费黄网站久久成人精品| 国产高清视频在线观看网站| 91狼人影院| 色精品久久人妻99蜜桃| 国内少妇人妻偷人精品xxx网站| 干丝袜人妻中文字幕| 精品一区二区免费观看| 看片在线看免费视频| 91在线精品国自产拍蜜月| 熟女电影av网| 久久久久性生活片| 国产蜜桃级精品一区二区三区| 免费在线观看日本一区| 99九九线精品视频在线观看视频| 91狼人影院| 搡老妇女老女人老熟妇| 特级一级黄色大片| 99九九线精品视频在线观看视频| 中文资源天堂在线| 久久久久久久午夜电影| 国产亚洲精品综合一区在线观看| 听说在线观看完整版免费高清| 成年版毛片免费区| 人妻少妇偷人精品九色| 午夜福利欧美成人| 啦啦啦韩国在线观看视频| 我要看日韩黄色一级片| 在线a可以看的网站| 日韩 亚洲 欧美在线| 国产一区二区在线观看日韩| 亚洲三级黄色毛片| 国内毛片毛片毛片毛片毛片| 亚洲成人免费电影在线观看| 麻豆精品久久久久久蜜桃| 午夜a级毛片| 精品国内亚洲2022精品成人| 欧美成人性av电影在线观看| 色哟哟哟哟哟哟| 又爽又黄无遮挡网站| 国产精品免费一区二区三区在线| 亚洲国产日韩欧美精品在线观看| 久久久久国内视频| 婷婷精品国产亚洲av在线| 色尼玛亚洲综合影院| 久久久久久久午夜电影| 亚洲精品一卡2卡三卡4卡5卡| 熟女人妻精品中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 欧美最黄视频在线播放免费| 一夜夜www| 999久久久精品免费观看国产| 女人被狂操c到高潮| 日韩欧美三级三区| 亚洲国产日韩欧美精品在线观看| 99riav亚洲国产免费| 99国产极品粉嫩在线观看| 高清在线国产一区| 国产一区二区三区在线臀色熟女| 国产乱人伦免费视频| 一进一出抽搐动态| 国产蜜桃级精品一区二区三区| 波多野结衣巨乳人妻| 日本 欧美在线| 精品一区二区三区人妻视频| 亚洲成人久久爱视频| 亚洲av中文av极速乱 | 国产熟女欧美一区二区| 免费人成视频x8x8入口观看| 日韩,欧美,国产一区二区三区 | 成人av在线播放网站| 国产高潮美女av| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 久久久久性生活片| 国产一区二区三区av在线 | 少妇裸体淫交视频免费看高清| 成人特级黄色片久久久久久久| 亚洲人成伊人成综合网2020| 我的女老师完整版在线观看| 亚洲,欧美,日韩| 给我免费播放毛片高清在线观看| 国产精品福利在线免费观看| 老师上课跳d突然被开到最大视频| av在线天堂中文字幕| 欧美日韩综合久久久久久 | 亚洲成人久久性| 一卡2卡三卡四卡精品乱码亚洲| 国产一区二区激情短视频| 婷婷精品国产亚洲av| 亚洲欧美激情综合另类| 日韩欧美三级三区| 国产精品乱码一区二三区的特点| 久久久色成人| 成人国产麻豆网| 国产激情偷乱视频一区二区| videossex国产| 亚洲 国产 在线| 天堂影院成人在线观看| 国产精品精品国产色婷婷| 欧美日韩综合久久久久久 | 成人国产一区最新在线观看| 天堂影院成人在线观看| 午夜亚洲福利在线播放| 国产高潮美女av| 色尼玛亚洲综合影院| 99九九线精品视频在线观看视频| 欧美日韩黄片免| 极品教师在线视频| 琪琪午夜伦伦电影理论片6080| 亚洲经典国产精华液单| 久久99热6这里只有精品| 国产精品人妻久久久影院| 又黄又爽又免费观看的视频| 国产成人av教育| 亚洲一级一片aⅴ在线观看| 天天一区二区日本电影三级| 麻豆成人午夜福利视频| 99久久精品热视频| 国产伦在线观看视频一区| 久久九九热精品免费| 少妇丰满av| 91在线观看av| 欧美色视频一区免费| 97人妻精品一区二区三区麻豆| 国产激情偷乱视频一区二区| 性色avwww在线观看| 亚洲成人中文字幕在线播放| 窝窝影院91人妻| 日韩欧美精品免费久久| 搡老岳熟女国产| 赤兔流量卡办理| 99riav亚洲国产免费| 亚洲精品日韩av片在线观看| 真人做人爱边吃奶动态| 春色校园在线视频观看| 亚洲精品成人久久久久久| 中文字幕高清在线视频| 给我免费播放毛片高清在线观看| avwww免费| 亚洲无线观看免费| 成人av在线播放网站| 我的女老师完整版在线观看| 国内精品一区二区在线观看| 韩国av在线不卡| 天堂√8在线中文| 国产乱人伦免费视频|