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

    小流域地下水與地表水微生物組成結(jié)構(gòu)差異及其成因分析

    2012-12-25 02:09:02黃小蘭陳建耀周世寧謝麗純楊雪韻蔣華波董林垚江西師范大學地理與環(huán)境學院鄱陽湖濕地與流域研究教育部重點實驗室江西南昌00中山大學地理科學與規(guī)劃學院廣東廣州51075中山大學生命科學學院廣東廣州51075
    中國環(huán)境科學 2012年9期
    關(guān)鍵詞:井水基巖花崗巖

    黃小蘭,陳建耀,周世寧,謝麗純,楊雪韻,蔣華波,董林垚 (1.江西師范大學地理與環(huán)境學院,鄱陽湖濕地與流域研究教育部重點實驗室,江西 南昌 00;.中山大學地理科學與規(guī)劃學院,廣東 廣州 51075;.中山大學生命科學學院,廣東 廣州 51075)

    小流域地下水與地表水微生物組成結(jié)構(gòu)差異及其成因分析

    黃小蘭1,2,陳建耀2*,周世寧3,謝麗純2,楊雪韻2,蔣華波2,董林垚2(1.江西師范大學地理與環(huán)境學院,鄱陽湖濕地與流域研究教育部重點實驗室,江西 南昌 330022;2.中山大學地理科學與規(guī)劃學院,廣東 廣州 510275;3.中山大學生命科學學院,廣東 廣州 510275)

    對中山大學濱海水循環(huán)綜合試驗基地的河道水和地下水進行采樣,運用16S rRNA基因分析方法,研究微生物的群落組成、結(jié)構(gòu),同時,通過Mothur軟件分析不同采樣點的系統(tǒng)發(fā)育多樣性和文庫稀釋曲線.結(jié)果表明,上下游河道水和地下水均以變形菌(Proteobacteria),厚壁菌(Firmicutes),古菌(Archaea)為優(yōu)勢種群,變形菌在上游堰、下游堰、R1、R2、R3分別占各自總數(shù)的29.7%、30.3%、28.1%、20.9%、14.2%;厚壁菌在上游堰、下游堰、R1、R2、R3分別占各自總數(shù)的9.7%、20.3%、19.8%、8.6%、39.6%;古菌(不可培養(yǎng)的泉古菌和廣古菌)在上游堰、下游堰、R1、R2、R3分別占各自總數(shù)的14.2%、12.3%、11.7%、12.1%、3.5%.表明地下水、地表水微生物組成結(jié)構(gòu)有一定的趨同性.另一方面,酸桿菌(Acidobacteria)、硝化螺旋菌(Nitrospira)、擬桿菌(Bacteroidetes) 和Candidate division OPx等在不同采樣點的豐度各異,不同埋藏條件下的地下水和地表水多樣性指數(shù)與稀釋曲線有所不同,揭示小流域的微生物組成、結(jié)構(gòu)具有空間差異.研究顯示,受周圍水文地質(zhì)條件、補給差異等影響,微生物種群隨之演化,微生物特征與水文地質(zhì)、化學離子密切相關(guān).

    小流域;地下水;地表水;微生物;差異;成因

    地表水-地下水的相互作用在全球水循環(huán)中發(fā)揮著重要的作用.國內(nèi)外學者在地表水與地下水的相互作用[1-2]的研究取得了顯著成果,隨著多學科方法的綜合應用及交叉研究的加強,微生物-水-巖相互作用成為區(qū)域地下水地球化學研究的熱點,其中物理、化學、生物的耦合作用尤為重要[3].當前地表-地下水的相互作用幾乎涵蓋了水循環(huán)過程中降雨[4]、截留[5]、蒸發(fā)[6]、下滲[7]等數(shù)學模擬和物理過程,以及滲濾作用[8]、污染物遷移[9]、氯氟烴(CFCs)測年[10]等化學過程,而生物過程研究的內(nèi)容主要出現(xiàn)在綜述性文章或者針對某種污染的生物修復[11]方面,原位地表水與地下水相互作用的生物過程研究則較少且偏宏觀[12].基于此,本研究從微生物的角度出發(fā),通過小流域地下水與地表水微生物組成結(jié)構(gòu)差異及其成因分析,研究地表-地下水微生物分布特征與水文地質(zhì)、化學離子之間的協(xié)同聯(lián)系,探討華南濱?;◢弾r小流域地下水與河道地表水的交互、演化和環(huán)境作用,以期為流域水資源規(guī)劃、生物多樣性和流域生態(tài)系統(tǒng)管理等研究提供參考.

    1 材料與方法

    1.1 研究區(qū)概況

    珠海市位于北緯21°48′~22°27′、東經(jīng)113°03′~114°19′.年降水量 1770~2300mm,屬亞熱帶季風氣候區(qū),海洋對本地氣候的調(diào)節(jié)作用十分明顯.選擇位于珠海市區(qū)的中山大學濱海水循環(huán)綜合試驗基地為研究對象,基地以雨水、地表水、土壤水、地下水、海水的水分相互轉(zhuǎn)化與物質(zhì)遷移為理論框架,設立了相關(guān)的試驗觀測設施,整個試驗設施分空中、地表、地下3大模塊,本研究對象的地表模塊為試驗小流域的上游堰和下游堰(表1,圖 1),其基巖類型為強風化花崗巖,土壤主要是花崗巖風化土壤,土壤厚度 0~3m,該流域的地質(zhì)特點在珠江三角洲濱海地區(qū)十分典型[13];地下模塊為R1~R3觀測井(表1,圖1),根據(jù)陜西工程勘察研究院的實地鉆探結(jié)果,R井組自上而下垂向分布有第四系全新統(tǒng)海陸交互相沉積層(Q4mc)、燕山期花崗巖(γ52),地下水為上層潛水和承壓基巖裂隙水.校區(qū)內(nèi)地下水主要受大氣降水、地下徑流補給與附近地表水體滲透補給,以蒸發(fā)及地下水徑流的方式排泄,研究區(qū)上游堰和下游堰平均流速3~5m/s,地下水流速介于110~132m/a之間.采樣時間在2010年6月大雨過后,屬雨季.

    圖1 河水-地下水水力聯(lián)系概念Fig.1 The water power correlation conception diagraph between river water and groundwater

    1.2 研究方法

    野外調(diào)研與實驗室分析相結(jié)合.野外采樣點現(xiàn)場采用GPS定位,用Horiba EC/pH計即時測量水電導率(EC)和pH值等(表1);實驗室分析先在現(xiàn)場用無菌聚乙烯瓶采集水樣 2L/井,隨即經(jīng)0.8μm和0.22μm濾膜2級過濾,濾后的水樣保存好送實驗室進行化學離子分析,將 0.22μm濾膜放入消毒后的冷凍存管中,于冰盒中保存,當天運回實驗室-20℃以下冰存.

    實驗室分析方面,采用酸堿滴定法輔以儀器設備分析主要化學組成成分(Ca2+、Mg2+、K+、Na+、Cl-、CO3-、SO42-、HC、Sidiss),然后用GW-Chart軟件作出Piper圖進行徑流不同采樣點的化學組成比較;同時,運用16S rRNA基因分析方法[14-15]檢測地下水與河道地表水的微生物種群,即首先采取改進的玻璃微珠法[16]提取水樣微生物總DNA,利用美國OMEGA公司膠回收試劑盒純化,置于-20℃保存?zhèn)溆?隨后按以下步驟進行實驗:(1)PCR(聚合酶鏈式反應)擴增(擴增采用的引物見表 2):原核類群 PCR反應過程為94℃變性4min,然后94℃,45s,52~55℃ 45s,72℃1min循環(huán)30次,72℃延伸10min;古菌PCR反應過程為 94℃變性 4min,94℃ 1min,50~60℃1min,72℃ 1min循環(huán) 26次,72℃延伸 10min. (2)PCR產(chǎn)物克隆及RFLP(限制性片斷長度多態(tài)性)分析:以2μL陽性克隆的菌落PCR產(chǎn)物為模板,分別選用RsaI(GT’ AC)和MspI(C’ CGG) 2種內(nèi)切酶對 DNA沉淀進行完全酶切,酶切產(chǎn)物用3%瓊脂糖凝膠電泳進行檢測,分析不同陽性克隆所產(chǎn)生的PCR產(chǎn)物經(jīng)酶切后的譜型,相同的譜型作為一個OTU,將具有不同譜型的克隆子送上海生工生物技術(shù)公司完成序列測定.(3)通過BLAST 程序?qū)y序結(jié)果輸入 GenBank (http://www.ncbi.nlm.nih.gov/)中進行序列相似性比對分析,通過Mothur和Excel軟件生成物種多樣性指數(shù)和稀釋曲線.

    表1 樣本基本信息Table 1 Basic information of the samples

    表2 PCR擴增采用的引物[17-19]Table 2 Primers used for PCR

    2 結(jié)果與討論

    水樣微生物 16S rDNA的部分序列在Genbank核苷酸數(shù)據(jù)庫中接受號為 HQ596988-HQ597006(網(wǎng)址為 http://www.ncbi.nlm.nih.gov/ nuccore/HQ596988- HQ597006).

    2.1 微生物種群與多樣性

    2.1.1 微生物優(yōu)勢種群 數(shù)理統(tǒng)計分析中,按占測序總數(shù) 5%以上[20]的同類微生物確定為優(yōu)勢種群(表3),試驗流域河道水和R井組地下水的微生物均以變形菌門(Proteobacteria),厚壁菌門(Firmicutes),古菌(Archaea)為優(yōu)勢種群,但同為優(yōu)勢群的變形菌門(Proteobacteria),在上游堰分布以α、 β、δ-proteobacteria為優(yōu)勢種群;在下游堰分布以γ、α、β-proteobacteria為優(yōu)勢種群;R1以 α、β、γ-proteobacteria 為 優(yōu) 勢 種 群 ;R2 以 β、δ-proteobacteria為優(yōu)勢種群;R3以 γ、βproteobacteria為優(yōu)勢種群;Archaea除R3外,其余4個水樣均屬優(yōu)勢微生物種群.Candidate division OPx(Candidate division OP3、OP8、OP11等總稱為OPx,x為數(shù)字)在上游堰、下游堰、R2井水樣中大量存在,而R1和R3井中卻沒有發(fā)現(xiàn).硝化螺旋菌(Nitrospirae)在R2和R3井有,而在R1井沒有檢測出,相應地,反硝化細菌(Denitrifying bacterium)在R2和R3井中也有少量存在.藍細菌在R2、下游堰、上游堰均有分布,而R1和R3井沒有發(fā)現(xiàn).

    表3 樣本優(yōu)勢微生物種群Table 3 Predominant microbial population of the samples

    圖2 各樣點的微生物聚類分析Fig.2 Cluster analysis of Ur, Ds, R1, R2, R3

    用MOTHUR軟件的UPGMA算法對微生物種群進行聚類分析(圖2),從UPGMA樹看出,上游堰河道水與R2地下水之間關(guān)系較密切,下游堰與R1關(guān)系較近,而R3單獨一枝,說明13.2m深的地下水微生物與其他4個樣點之間的差異更為明顯.

    系統(tǒng)發(fā)育多樣性分析運用MOTHUR軟件分析5個采樣點的系統(tǒng)發(fā)育多樣性(圖3),并生成16S rDNA克隆文庫的稀釋曲線(圖4).

    圖3 河道水和地下水的系統(tǒng)發(fā)育多樣性分析Fig.3 The richness/diversity estimators from groundwater and surface water samples

    結(jié)果表明,上游堰、下游堰和R2井水樣的香農(nóng)指數(shù)(Shannon)大于R1、R3井,Simpson、Chao1和 Ace指數(shù)也驗證了前者的微生物多樣性指數(shù)高于后者.

    圖4 16S rDNA克隆文庫的稀釋曲線Fig.4 The rarefaction curves from groundwater and surface water samples

    當克隆文庫的稀釋曲線趨于平緩或者趨近平臺期時,就可以認為文庫的庫容已經(jīng)足夠[21]. R3井和R1井的16S rDNA克隆文庫的稀釋曲線首先趨于平緩,說明R3井和R1井的多樣性指數(shù)雖然較低,但文庫已基本能夠代表樣品中微生物的多樣性.理論上,R1井樣品來自淺層非承壓含水層,與地表非飽和帶聯(lián)系密切,受地表水(河道)與雨水補給,可供微生物生長繁殖的營養(yǎng)物豐足,因而 R1井水的微生物多樣性指數(shù)應該也較高,其主要微生物種群與河道水基本一致,只是各種群的豐度稍有差異,這在本次采樣分析的前1次采樣結(jié)果中已有驗證(圖 5).本次實驗結(jié)果出現(xiàn)R1井水微生物多樣性指數(shù)低的原因可能是現(xiàn)場過濾水樣時耗費時間稍長,導致優(yōu)勢微生物生長過快而抑制了其他微生物的檢出.

    2.2 成因分析

    2.2.1 水文地質(zhì)成因 一方面,地下水、地表水微生物組成結(jié)構(gòu)有一定的趨同性;另一方面,受水文地質(zhì)條件,補給差異等影響,各層微生物種群出現(xiàn)差異并體現(xiàn)地表-地下水之關(guān)聯(lián).R井組所監(jiān)測的地下水既包括潛水也包括承壓的基巖裂隙水,R1井地質(zhì)構(gòu)造以細砂為主,屬第四系全新統(tǒng)海陸交互相沉積層(Q4mc);R2井以中風化花崗巖為主,R3井以微風化花崗巖為主,為燕山期花崗巖(γ52).一般地,沉積層具有較豐富的有機質(zhì),其中所含的活性有機碳較高,可供微生物利用的物質(zhì)也相對較高.林立虹等[22]在距地表2.8km深處巖石裂隙的水脈中,發(fā)現(xiàn)僅僅依靠地質(zhì)作用產(chǎn)生的硫酸鹽、氫氣等能量而生存的微生物種群.R3觀測井打入基巖,由于細菌硫酸鹽還原作用(2CH2O+SO42?→H2S+2HCO3?)而富硫化氫,這種環(huán)境顯然不適合好氧微生物發(fā)展,在這缺氧的富硫化氫環(huán)境中,R3井水發(fā)現(xiàn)大量的Bacillus芽孢桿菌(屬于厚壁菌門 Firmicutes),芽胞是細菌在極端環(huán)境下產(chǎn)生的一種抗逆性強的構(gòu)造,莢膜也是細菌在特殊環(huán)境下產(chǎn)生的一種可使細菌免受惡劣環(huán)境損傷的細菌附屬物[23],R3井為基巖裂隙水,其水量受補給源控制也不夠豐足,處于還原環(huán)境,在現(xiàn)場采水樣時有明顯的臭雞蛋味,因此在R3井水的分析結(jié)果中,各種芽孢桿菌占絕對優(yōu)勢,而其他微生物種群明顯少于另外 4個水樣.R2井表現(xiàn)一定的承壓與缺氧環(huán)境,相比而言,R2井水的微生物多樣性較豐富,可能與其含水層包涵更多的花崗巖風化破碎層因而與河道水聯(lián)通有關(guān).

    上游堰和下游堰的基巖類型為強風化花崗巖,土壤主要是花崗巖風化土壤,其河道水的補給源除地表徑流、壤中流之外,基流也是重要組成部分.曾松青等[24-25]利用同位素研究中山大學濱海水循環(huán)綜合試驗基地的花崗巖小流域,結(jié)果表明,基流占小流域全年徑流量的比例為 68%~74%.這與 Katsuyama等[26]發(fā)現(xiàn)一個花崗巖小流域的基巖地下水對河岸帶有大量的補給一致.因此,盡管采樣前下大雨,上游堰和下游堰河道水均有部分來自地下水.在微生物組成結(jié)構(gòu)上,表現(xiàn)為上游堰河道水與 R2井水較為接近,均是古菌占最大優(yōu)勢,微生物香農(nóng)指數(shù)也相近;下游堰河道水與 R1井水較為接近,均是厚壁菌門(Firmicutes)占最大優(yōu)勢,這從它們的水力聯(lián)系和地質(zhì)構(gòu)造上(圖 5)可得以解釋:上游堰與基巖風化殼(燕山期花崗巖(γ52))連通,R2井埋藏條件也主要為燕山期花崗巖風化殼(γ52);下游堰與第四系全新統(tǒng)海陸交互相沉積層(Q4mc)聯(lián)系密切,R1井埋藏條件也主要為第四系全新統(tǒng)海陸交互相沉積層(Q4mc),與微生物聚類結(jié)果一致.

    對比分析試驗小組于 6月份實測上下游堰水位與R1~R3井水位變化資料(圖6),上游堰的水位變化幅度大,依次是下游堰、R1、R2、R3井.R1井埋藏條件為海陸交互相沉積層,與下游堰河床類似,水位變化幅度自然較大;R2井與R3井的主要差別在于前者包涵更多的花崗巖風化破碎層,而后者破碎層較少,但有一定裂隙發(fā)育.因而,河道水和R2井地下水比R3井地下水有更強的水力聯(lián)系;R3井地下水則比較平穩(wěn),受地表水影響較小.從時間來看,井水對河道水漲落無明顯的延遲,即地下水與河道水微生物的異同不受延遲影響.微生物分布方面,R2井Archaea為最優(yōu)勢種群,與上游堰一致,其余的優(yōu)勢種群如β-proteobacteria,Candidate division OPx, Firmicutes, δ-proteobacteria等,也與上游堰基本相同,但硝化細菌(Nitrospira)在R2井為優(yōu)勢微生物種群而在上游堰卻未發(fā)現(xiàn),是因為硝化細菌大多數(shù)為專性化能自養(yǎng)型,只有少數(shù)為兼性自養(yǎng)型,大多數(shù)好氧,怕光.R1井Firmicutes為最優(yōu)勢種群,與下游堰一致,其他如 α, β, γ-proteobacteria、Archaea等也與下游堰基本類同,但R1井有優(yōu)勢種群酸桿菌(Acidobacteria)而下游堰只有少量酸桿菌存在,這與 R1井處于果園,與地表有較多的枯葉腐爛相關(guān).R3井地下水處于還原環(huán)境,因而屬于厚壁菌門(Firmicutes)的各種芽孢桿菌占微生物總數(shù)的 39.6%,為絕對優(yōu)勢,其余優(yōu)勢種群γ-proteobacteria,β-proteobacteria,Nitrospira,Bacte roidetes所占比例明顯較少,處于5%~10%之間.

    圖5 R1井地下水與河道水主要微生物種群Fig.5 Microbial communities of analyzed clones from R1 and surface water samples

    圖6 地表水與地下水采樣點水位變化比較Fig.6 The comparing of water-level fluctuations between river water and groundwater

    2.2.2 化學離子成因 地表地下水的微生物組 成、結(jié)構(gòu)差異也表現(xiàn)與周邊環(huán)境的化學離子組成相適應,例如,硝化螺旋菌(Nitrospirae)在R2和R3井地下水中有,而在 R1井水中沒有檢測出,硝化螺旋菌對氮的硝化和亞硝酸鹽氧化起作用[27].原先R井所在位置有工人居住,受雞、鴨、羊等排泄物污染,但自 2009年始,中山大學珠海校區(qū)美化清潔校園環(huán)境,將果園的居民遷出,R井附近基本保持自然生態(tài)環(huán)境,最淺層井 R1由雨水和周邊地表水入滲補給后原污染物已慢慢稀釋和降解,R2和R3井為基巖裂隙水,原先雞、鴨等排泄污染物尚存留于基巖風化殼層中,因此,隨著周邊環(huán)境的日益變化,微生物組成、結(jié)構(gòu)不斷地改變.

    圖7 研究區(qū)地表水和地下水主要化學離子Piper圖Fig.7 Piper diagram for river water and groundwater in the Zhuhai Campus catchment

    試驗流域基巖由單一的黑云母花崗巖組成,徑流的化學組成屬于Sidiss-Na+-HCO3?-Cl?型,以溶解Si、Na+、HCO3?和Cl-為主要成分,Ca2+、SO42?則變異系數(shù)較大[28],微生物種群中有少量硫酸鹽還原菌(sulfate-reducing bacterium)存在.在補給區(qū)R井組,地下水徑流條件較好,地下水在含水層中停留時間短,水交替迅速,含水層中易溶組分如Cl-、SO42-、Na+、K+等不斷被淋濾并由地下徑流帶走,離子以 Ca2+、HCO3-占優(yōu)勢,形成低礦化度的Ca2+- Na+-HCO3--Cl?型地下水,R井組地下水檢 測 出 屬 于 γ-proteobacteria 的 海 桿 菌(Marinobacter)和屬于β-proteobacteria的脫氯單胞菌屬(Dechloromonas),說明 R井組地下水有Na+、Cl?等離子循環(huán).高全洲等[29]研究表明,珠江三角洲濱海地區(qū)自然化學風化、海鹽沉降和人類活動對徑流化學組成有近于相等的貢獻水平,此外,自然化學風化過程所產(chǎn)生的 HCO3-離子響應洪水的稀釋效應,海鹽沉降通量也響應徑流的稀釋效應.從主要化學離子Piper圖(圖7)也可以看出,雨季尤其是暴雨時上下游河道水的離子成分更接近雨水,稀釋效應顯著.盡管如此,河道水與井水中的優(yōu)勢微生物種群一致性較高,不受這種稀釋效應影響,此特征表明可從微生物角度辨識地下水/地表水的相互作用.本次水樣中自然化學風化和海鹽沉降的貢獻份額都相對減小,酸沉降的貢獻稍有增加,如現(xiàn)場測定的 pH值偏弱酸性(表1),這在微生物組成、結(jié)構(gòu)中也有所反映,表現(xiàn)在本次測序的微生物種群中酸桿菌(Acidobacteria,大多由嗜酸菌組成)均有存在,尤其是在R1井中已達到了優(yōu)勢種群水平.

    3 結(jié)論

    3.1 微生物組成結(jié)構(gòu)方面,淺層潛水與地表非飽和帶(土壤)聯(lián)系密切,受地表水(河道)與雨水補給,一般為氧化環(huán)境,其中的微生物組成結(jié)構(gòu)與下游河道水相似,表現(xiàn)為 Firmicutes為最優(yōu)勢種群,其他如,α, β, γ-proteobacteria、Archaea等也與下游河道水基本類同;處于花崗巖破碎帶的基巖裂隙水,其略帶承壓的含水層為主要聯(lián)系通道與上游河道水交互,相應的微生物組成結(jié)構(gòu)與河道水有較大的相似性,表現(xiàn)為 Archaea為最優(yōu)勢種群,其余的優(yōu)勢種群如β-proteobacteria,Candidate division OPx,Firmicutes,δ-proteobacteria等,也與上游堰基本相同,且微生物多樣性指數(shù)和稀釋曲線也較為接近;處于裂隙發(fā)育但破碎層較少的基巖裂隙水,其承壓含水層受補給源控制水量不夠豐足,地下水位比較平穩(wěn),處于還原環(huán)境,各種芽孢桿菌占絕對優(yōu)勢,而其他微生物種群較少.

    3.2 受水文地質(zhì)和化學離子等綜合影響,酸桿菌在淺層潛水中表現(xiàn)為優(yōu)勢菌,藍細菌僅在花崗巖破碎帶的基巖裂隙水和河道水中系統(tǒng)發(fā)育同源性相近,硝化與反硝化菌主要在基巖裂隙水中存在,這些差異的存在表明在每個獨特地質(zhì)構(gòu)造的含水層中,地下水與地表水的交互程度不同,微生物的組成結(jié)構(gòu)隨不同地理時空的物理化學差異而變化并且適應.

    [1] VanderVelde Y, deRooij G H, Torfs P J J F. Catchment-scale non-linear groundwater-surface water interactions in densely drained lowland catchments [J]. Hydrology and Earth System Sicence, 2009,13(10):1867-1885.

    [2] Bailly-Comte V, Jourde H, Pistre S. Conceptualization and classification of groundwater -surface water hydrodynamic interactions in karst watersheds: Case of the karst watershed of the Coulazou River (Southern France) [J]. Journal of Hydrology, 2009,376:456-462.

    [3] 滕彥國,左 銳,王金生,等.區(qū)域地下水演化的地球化學研究進展 [J]. 水科學進展, 2010,21(1):127-135.

    [4] Vidal J P, Wade S D. Multimodel projections of catchment-scale precipitation regime [J]. Journal of Hydrology, 2008,353: 143-158.

    [5] Schellekens J. The interception and runoff generating processes in the bisley catchment, luquillo experimental forest, Puerto Rico [J]. Physics and Chemistry of the Earth, 2000,25(13):659-664.

    [6] Zhang Y Q, Chiew F H S, Zhang L, et al. Estimating catchment evaporation and runoff using MODIS leaf area index and the Penman–Monteith Equation [J]. Water Resources Research, DOI: 10.1029/2007WR006563.

    [7] Yimer F, Messing I, Ledin S, et al. Effects of different land use types on infiltration capacity in a catchment in the highlands of Ethiopia [J]. Soil Use and Management, 2008,24(4):344-349.

    [8] Fukada T, Hiscock K M, Dennis P F, et al. A dual isotope approach to identify denitrification in groundwater at a river-bank infiltration site[J]. Water Research, 2003,37:3070-3078.

    [9] Willhelm S R, Schiff S L, Robertson W D. Biochemical evolution of domestic waste water in septic systems: 2.Application of conceptual model in sandy aquifers [J]. Ground Water, 1996,34: 853-864.

    [10] 秦大軍.地下水CFC定年方法及應用 [J]. 地下水, 2005,27(6): 435-437.

    [11] 李昊翔.螺旋藻對重金屬的耐受性和吸附研究 [D]. 杭州:浙江大學, 2005.

    [12] 滕彥國,左 銳,王金生.地表水-地下水的交錯帶及其生態(tài)功能[J]. 地球與環(huán)境, 2007,35(1):1-8.

    [13] 孫 濤.新編華南花崗巖分布圖及其說明 [J]. 地質(zhì)通報, 2006,25(3):332-335.

    [14] 黃小蘭,陳建耀,周世寧,等.珠海市海-陸交錯帶水環(huán)境原核生物多樣性 [J]. 應用生態(tài)學報, 2010,21(2):452-457.

    [15] 黃小蘭,陳建耀,謝麗純,等.小流域地下水流微生物空間差異分析[J].環(huán)境科學,2010,31(10): 2299-2304.

    [16] 肖 凱,曹理想,陸勇軍,等.廣東省金山溫泉沉積物微生物多樣性初步分析 [J]. 微生物學報, 2008,48(6):717-724.

    [17] David M W, Michael J F, Stephen C N, et al. natural view of microbial biodiversity within hot spring cyanbacterial mat communities [J]. Microbiol. Mol. Biol. Rev., 1998,62:1353-1370.

    [18] Godon J J, Zumstein E, Dahert P, et al. Molecular microbial diversity of an anaerobic digestor as determined by small-subunit rDNA sequence analysis [J]. Appl. Environ. Microbiol., 1997, 63:2802-2813.

    [19] Nakamura Y, Kaneko T, Sato S, et a1. Complete genome structure of the thermophilic cyanobacteria Therraosynechococcuselongatus BP-1(Supplement) [J]. DNA Res., 2002,93(4):135-148.

    [20] Huang X L, Chen J Y, Zhou S N, et al. Vertical microbiological variation of a coastal aquifer in southern China [J]. Water Science and Technology, 2012,65(4):703-712.

    [21] 關(guān)統(tǒng)偉,吳晉元,職曉陽,等.硝爾庫勒湖沉積物中非培養(yǎng)放線菌多樣性 [J]. 微生物學報, 2008,48(7):851-856.

    [22] Lin L H, Wang P L, Douglas R, et al. Long-term sustain ability of a high-energy, low-diversity crustal biome [J]. Science, 2006,314(5798):479-482.

    [23] 李 鵬,李志會,金琳琳,等.山東肥城寒武紀竹葉狀石灰?guī)r菌群分析 [J]. 中國微生態(tài)學雜志, 2009,21(8):721-726.

    [24] 曾松青,陳建耀,付叢生.同位素在小流域基流計算中的應用研究 [J]. 水文, 2010,30(2):20-24.

    [25] 付叢生,陳建耀,曾松青,等.濱海地區(qū)潮汐對地下水位變化影響的統(tǒng)計學分析 [J]. 水利學報, 2008,39(12):1365-1376.

    [26] Katsuyama M, Ohte N, Kabeya N. Effects of bedrock permeability on hillslope and riparian groundwater dynamics in a weathered granite catchment [J]. Water Resources Research, 2005, DOI: 10.1029/2004WR003275.

    [27] 黃小蘭,陳建耀.微生物應用于污水污泥處理的研究 [J]. 亞熱帶資源與環(huán)境學報, 2010,5(1):48-55.

    [28] 趙新鋒.城市化地下水環(huán)境效應研究 [D]. 廣州:中山大學, 2008.

    [29] 高全洲,陶貞.華南濱?;◢弾r丘陵的化學風化與化學徑流 [J].中國科學D輯:地球科學, 2010,40(6):758-767.

    Cause on the differences of microorganism between small watershed surface water and underground water.

    HUANG Xiao-lan1,2, CHEN Jian-yao2*, ZHOU Shi-ning3, XIE Li-chun2, YANG Xue-yun2, JIANG Hua-bo2, DONG Lin-yao2(1.Key Laboratory of Poyang Lake Wetland and Watershed Research, Ministry of Education, School of Geography and Environment, Jiangxi Normal University, Nanchang 330022, China;2.School of Geography and Planning, Sun Yat-sen University, Guangzhou 510275, China;3.School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, China). China Environmental Science, 2012,32(9):1647~1654

    The surface and underground water of Sun Yat-sen University water cycle experiment base was sampled and analyzed by 16S rRNA gene analytical methods to acquire the data of microorganism communities’ distribution, then, the rarefaction curves and richness/diversity estimators were calculated and compared by Mothur software. Microbiological communities in groundwater were similar to those in surface water , with the predominant group of Proteobacterium, Firmicutes, Uncultured Archaeon(uncultured Crenarchaeote and Euryarchaeote). Proteobacteria made up 29.7%、30.3%、28.1%、20.9%、14.2% of the microbial communities from Shangyouyan, Xiayouyan, R1, R2, R3 , respectively. The same samples from Shangyouyan, Xiayouyan, R1, R2, R3 also consisted of 9.7%、20.3%、19.8%、8.6%、39.6% Firmicutes and 14.2%、12.3%、11.7%、12.1%、3.5% Archaea, respectively. On the other hand, the data of richness and diversity estimators such as Chao1, Ace, Shannon, and Simpson indicated that microorganism diversities were different in groundwater and surface water, which showed the difference among those microorganism’s characteristics. The study concluded that microorganism evolves in response to the change of the hydrogeological and ionic conditions that are influenced by the water environment and the recharging of water resources.

    small watershed;groundwater;surface water;microorganism;difference;cause

    2012-02-09

    國家“973”項目(2012CB417003);廣東省自然科學基金項目(9251027501000021);江西省教育廳科技項目(GJJ12186);廣東水利科技創(chuàng)新與推廣項目(2009-2011);江西省博士后科研項目擇優(yōu)資助(2012)

    * 責任作者, 教授, chenjyao@mail.sysu.edu.cn

    X172

    A

    1000-6923(2012)09-1647-08

    黃小蘭(1968-),女,江西贛州人,副研究員,博士, 主要研究方向為城市化水環(huán)境效應.發(fā)表論文20余篇.

    猜你喜歡
    井水基巖花崗巖
    井水的奧秘
    花崗巖
    抗剝落劑TR-500S改善花崗巖混合料路用性能研究
    石油瀝青(2020年1期)2020-05-25 06:53:56
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    基于改進物元的大壩基巖安全評價
    河北省基巖熱儲開發(fā)利用前景
    盤一井水氡異常落實與分析
    黃驊埕古1井水位異常研究
    花崗巖儲集層隨鉆評價方法及應用
    裝飾塊料
    建筑與預算(2014年2期)2014-10-14 13:28:42
    在线永久观看黄色视频| 国产精品亚洲av一区麻豆| 亚洲 欧美 日韩 在线 免费| 99在线视频只有这里精品首页| 亚洲中文av在线| 亚洲电影在线观看av| 久久久久久免费高清国产稀缺| 日韩大尺度精品在线看网址| 99热这里只有是精品50| 国产亚洲欧美98| 在线看三级毛片| 最新在线观看一区二区三区| 欧美色欧美亚洲另类二区| 又爽又黄无遮挡网站| 欧美黄色淫秽网站| 精品午夜福利视频在线观看一区| 97人妻精品一区二区三区麻豆| 好男人电影高清在线观看| 日韩欧美国产一区二区入口| 日本免费一区二区三区高清不卡| 一级黄色大片毛片| 亚洲熟女毛片儿| 国产爱豆传媒在线观看 | 男人舔女人下体高潮全视频| 成年版毛片免费区| 777久久人妻少妇嫩草av网站| 少妇熟女aⅴ在线视频| 久久久久免费精品人妻一区二区| 国产1区2区3区精品| 看免费av毛片| 精品久久久久久久久久免费视频| 91大片在线观看| 精品国产乱子伦一区二区三区| 久久国产精品人妻蜜桃| 欧美一级a爱片免费观看看 | 精华霜和精华液先用哪个| 亚洲自拍偷在线| 欧美在线黄色| 国产av在哪里看| 国产亚洲精品第一综合不卡| 老熟妇乱子伦视频在线观看| 欧美日韩精品网址| 久久国产精品影院| 国产精品99久久99久久久不卡| 在线a可以看的网站| 小说图片视频综合网站| 一边摸一边做爽爽视频免费| 日本黄色视频三级网站网址| 精品一区二区三区视频在线观看免费| 亚洲人成网站在线播放欧美日韩| 亚洲成人免费电影在线观看| 欧美日韩黄片免| 国产乱人伦免费视频| 成年版毛片免费区| 黑人巨大精品欧美一区二区mp4| 午夜激情av网站| 国产精品亚洲一级av第二区| 免费在线观看影片大全网站| 欧美日韩国产亚洲二区| 亚洲欧美精品综合久久99| 成年女人毛片免费观看观看9| 69av精品久久久久久| 在线观看免费日韩欧美大片| 欧美色视频一区免费| 成人手机av| 99久久国产精品久久久| 夜夜夜夜夜久久久久| 久久久国产成人精品二区| 成熟少妇高潮喷水视频| 伦理电影免费视频| 中文字幕久久专区| 精品国内亚洲2022精品成人| 中文字幕av在线有码专区| 亚洲国产欧美网| 又黄又粗又硬又大视频| 手机成人av网站| 麻豆成人av在线观看| 国语自产精品视频在线第100页| 国产视频一区二区在线看| 日韩中文字幕欧美一区二区| 久久精品91蜜桃| 亚洲专区中文字幕在线| 午夜影院日韩av| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩无卡精品| 啦啦啦免费观看视频1| 欧美在线黄色| 草草在线视频免费看| 亚洲,欧美精品.| 亚洲精品国产精品久久久不卡| 亚洲国产欧美人成| 国产精品一区二区三区四区久久| 国产精品爽爽va在线观看网站| 亚洲av第一区精品v没综合| 人成视频在线观看免费观看| 黄片大片在线免费观看| 一个人观看的视频www高清免费观看 | 在线十欧美十亚洲十日本专区| 欧美高清成人免费视频www| 久久久久免费精品人妻一区二区| 久久久久免费精品人妻一区二区| 欧美乱色亚洲激情| 国产精品亚洲美女久久久| 又紧又爽又黄一区二区| 日本黄大片高清| 97碰自拍视频| 国产在线观看jvid| 国产精品免费一区二区三区在线| 五月玫瑰六月丁香| 嫩草影院精品99| 欧美zozozo另类| 欧美成人一区二区免费高清观看 | 精品欧美国产一区二区三| 国产一区二区激情短视频| or卡值多少钱| 18禁美女被吸乳视频| 黄片大片在线免费观看| 香蕉国产在线看| 麻豆国产av国片精品| 麻豆成人午夜福利视频| 成年免费大片在线观看| 级片在线观看| 一级毛片高清免费大全| 搞女人的毛片| 色哟哟哟哟哟哟| 神马国产精品三级电影在线观看 | 91国产中文字幕| 午夜老司机福利片| 欧美黄色淫秽网站| 国产高清videossex| 波多野结衣巨乳人妻| 欧美性长视频在线观看| 国产亚洲av高清不卡| 波多野结衣高清无吗| 少妇人妻一区二区三区视频| 久久亚洲精品不卡| 88av欧美| 制服丝袜大香蕉在线| 一边摸一边抽搐一进一小说| 亚洲人成电影免费在线| 国产视频内射| √禁漫天堂资源中文www| 日本一本二区三区精品| 麻豆一二三区av精品| 亚洲国产欧美网| 在线a可以看的网站| 每晚都被弄得嗷嗷叫到高潮| 国产又色又爽无遮挡免费看| 岛国视频午夜一区免费看| 露出奶头的视频| 午夜福利在线观看吧| 久久这里只有精品19| 日韩欧美在线二视频| 50天的宝宝边吃奶边哭怎么回事| 欧美最黄视频在线播放免费| 香蕉久久夜色| 麻豆成人av在线观看| av天堂在线播放| 国产精品,欧美在线| АⅤ资源中文在线天堂| 久久久国产精品麻豆| 欧美黑人巨大hd| 日本撒尿小便嘘嘘汇集6| 国产人伦9x9x在线观看| 麻豆一二三区av精品| 精品福利观看| 伊人久久大香线蕉亚洲五| 亚洲成人中文字幕在线播放| 精品乱码久久久久久99久播| 在线观看美女被高潮喷水网站 | 成人三级黄色视频| 午夜a级毛片| 国产精品一及| 久久久久久大精品| 国产成人精品久久二区二区91| 国产亚洲av嫩草精品影院| 国产亚洲欧美98| 99热这里只有是精品50| 熟妇人妻久久中文字幕3abv| 在线a可以看的网站| 日本免费一区二区三区高清不卡| 国产主播在线观看一区二区| 中文字幕人妻丝袜一区二区| 在线观看美女被高潮喷水网站 | 久久久水蜜桃国产精品网| 欧美一区二区国产精品久久精品 | 首页视频小说图片口味搜索| 欧美另类亚洲清纯唯美| 看免费av毛片| 国产伦在线观看视频一区| 欧美三级亚洲精品| 搡老岳熟女国产| videosex国产| 日韩av在线大香蕉| 久久这里只有精品19| 高潮久久久久久久久久久不卡| 欧美 亚洲 国产 日韩一| 精品国产超薄肉色丝袜足j| 一区福利在线观看| 国产精品九九99| 大型av网站在线播放| 欧美极品一区二区三区四区| 免费观看精品视频网站| 国产成+人综合+亚洲专区| 悠悠久久av| 一区福利在线观看| 99国产综合亚洲精品| 精品欧美国产一区二区三| 亚洲精品国产一区二区精华液| 久久精品国产综合久久久| 999精品在线视频| 在线视频色国产色| 亚洲熟妇中文字幕五十中出| 啪啪无遮挡十八禁网站| 久久久精品国产亚洲av高清涩受| 18禁黄网站禁片午夜丰满| 久久香蕉国产精品| 色综合婷婷激情| 亚洲性夜色夜夜综合| 在线观看免费视频日本深夜| 亚洲欧洲精品一区二区精品久久久| www.熟女人妻精品国产| 亚洲午夜理论影院| 精品国产亚洲在线| 国产精品久久久久久精品电影| 国产精品九九99| 毛片女人毛片| 高清在线国产一区| 好看av亚洲va欧美ⅴa在| 三级国产精品欧美在线观看 | 男人舔女人的私密视频| 亚洲精品色激情综合| 久久天躁狠狠躁夜夜2o2o| 成年免费大片在线观看| 波多野结衣巨乳人妻| 悠悠久久av| 波多野结衣高清作品| 国产黄色小视频在线观看| 亚洲成a人片在线一区二区| 男男h啪啪无遮挡| 又紧又爽又黄一区二区| 日韩精品免费视频一区二区三区| 亚洲av熟女| 亚洲人成伊人成综合网2020| 观看免费一级毛片| 黄色成人免费大全| 国产精品av视频在线免费观看| 日本a在线网址| 老司机午夜福利在线观看视频| 国产成人精品久久二区二区免费| 精品国产超薄肉色丝袜足j| avwww免费| 免费观看人在逋| aaaaa片日本免费| 成人欧美大片| 少妇粗大呻吟视频| 久久久水蜜桃国产精品网| 亚洲成av人片免费观看| 亚洲精品国产一区二区精华液| 日韩国内少妇激情av| 亚洲成人中文字幕在线播放| 亚洲精品中文字幕一二三四区| 在线看三级毛片| 黑人操中国人逼视频| 免费在线观看影片大全网站| 99热只有精品国产| 久久久久九九精品影院| 亚洲天堂国产精品一区在线| 欧美成人午夜精品| av欧美777| 免费人成视频x8x8入口观看| 可以免费在线观看a视频的电影网站| 免费在线观看视频国产中文字幕亚洲| 午夜福利在线在线| or卡值多少钱| e午夜精品久久久久久久| 一级毛片精品| 亚洲中文字幕日韩| 香蕉av资源在线| 18美女黄网站色大片免费观看| 日本撒尿小便嘘嘘汇集6| 窝窝影院91人妻| 日本黄大片高清| 亚洲18禁久久av| 日韩高清综合在线| 亚洲国产中文字幕在线视频| 欧美不卡视频在线免费观看 | 麻豆成人午夜福利视频| 亚洲人成网站高清观看| 小说图片视频综合网站| 亚洲精品久久成人aⅴ小说| 午夜亚洲福利在线播放| 国产69精品久久久久777片 | 看免费av毛片| 又粗又爽又猛毛片免费看| 亚洲精品美女久久久久99蜜臀| 国产熟女午夜一区二区三区| 高清在线国产一区| 亚洲av五月六月丁香网| 国产精品久久久av美女十八| 午夜福利免费观看在线| x7x7x7水蜜桃| 看黄色毛片网站| 日韩精品青青久久久久久| 国产黄a三级三级三级人| 国产av又大| 怎么达到女性高潮| 欧美+亚洲+日韩+国产| 亚洲狠狠婷婷综合久久图片| 99久久国产精品久久久| 99国产精品99久久久久| 99久久久亚洲精品蜜臀av| 国产精品一区二区三区四区久久| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| 在线观看免费视频日本深夜| 久久精品国产清高在天天线| 人成视频在线观看免费观看| 欧美中文日本在线观看视频| 午夜日韩欧美国产| 999精品在线视频| 久久精品aⅴ一区二区三区四区| 黄色视频,在线免费观看| 最近在线观看免费完整版| 欧美一级毛片孕妇| 日韩精品中文字幕看吧| 99re在线观看精品视频| av视频在线观看入口| 国产激情久久老熟女| 一级黄色大片毛片| 午夜免费观看网址| 一本一本综合久久| 中文字幕久久专区| 午夜福利在线在线| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 国产精品一区二区精品视频观看| 色综合亚洲欧美另类图片| 久久性视频一级片| 久久国产精品人妻蜜桃| 亚洲欧美日韩高清在线视频| 午夜视频精品福利| 亚洲熟妇中文字幕五十中出| 精品熟女少妇八av免费久了| 9191精品国产免费久久| 日本五十路高清| 欧美在线一区亚洲| 性欧美人与动物交配| 亚洲国产精品合色在线| 亚洲美女视频黄频| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 男人舔女人的私密视频| 午夜亚洲福利在线播放| 国产精品永久免费网站| 精品一区二区三区视频在线观看免费| 国产精品日韩av在线免费观看| 国产成人av教育| 搞女人的毛片| 久久这里只有精品中国| 国产精品久久视频播放| 亚洲av成人精品一区久久| 国产亚洲精品久久久久久毛片| 波多野结衣高清作品| 国产真实乱freesex| 欧美乱色亚洲激情| 国产高清激情床上av| 午夜免费激情av| 亚洲中文日韩欧美视频| 天天躁狠狠躁夜夜躁狠狠躁| 69av精品久久久久久| 欧美大码av| 国产高清视频在线播放一区| 国产三级黄色录像| bbb黄色大片| 好看av亚洲va欧美ⅴa在| 午夜a级毛片| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 亚洲全国av大片| 两个人视频免费观看高清| 久久国产精品人妻蜜桃| 午夜视频精品福利| 亚洲性夜色夜夜综合| 精品欧美一区二区三区在线| 国产亚洲av高清不卡| 一二三四社区在线视频社区8| 很黄的视频免费| 欧美不卡视频在线免费观看 | 亚洲一区高清亚洲精品| 亚洲国产精品合色在线| 国产一级毛片七仙女欲春2| www.999成人在线观看| 欧美在线一区亚洲| 成熟少妇高潮喷水视频| 三级国产精品欧美在线观看 | 亚洲精品美女久久久久99蜜臀| 午夜老司机福利片| 亚洲最大成人中文| 99热只有精品国产| 午夜老司机福利片| 超碰成人久久| 免费av毛片视频| 亚洲人成电影免费在线| 亚洲成av人片在线播放无| 18禁黄网站禁片免费观看直播| www.www免费av| 天堂av国产一区二区熟女人妻 | 1024香蕉在线观看| 日本撒尿小便嘘嘘汇集6| 欧美乱码精品一区二区三区| 97人妻精品一区二区三区麻豆| 91国产中文字幕| 韩国av一区二区三区四区| 国产精品av视频在线免费观看| 两性夫妻黄色片| 国产熟女xx| 国产爱豆传媒在线观看 | 久久国产精品人妻蜜桃| 亚洲 国产 在线| 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 日韩欧美国产一区二区入口| 久久天躁狠狠躁夜夜2o2o| 精品高清国产在线一区| 久久婷婷人人爽人人干人人爱| 国产欧美日韩一区二区精品| 久久精品亚洲精品国产色婷小说| 久久亚洲精品不卡| av福利片在线| 又黄又爽又免费观看的视频| 亚洲精品美女久久久久99蜜臀| www国产在线视频色| 国产精品亚洲一级av第二区| 亚洲av日韩精品久久久久久密| 色综合欧美亚洲国产小说| 999精品在线视频| 日本a在线网址| 啦啦啦观看免费观看视频高清| 2021天堂中文幕一二区在线观| a级毛片在线看网站| 中文字幕高清在线视频| 他把我摸到了高潮在线观看| 精品一区二区三区四区五区乱码| 小说图片视频综合网站| 嫩草影院精品99| 正在播放国产对白刺激| 午夜免费成人在线视频| 青草久久国产| 岛国在线观看网站| 国产三级黄色录像| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 国产成人影院久久av| svipshipincom国产片| 久久久精品国产亚洲av高清涩受| 中文亚洲av片在线观看爽| 日本a在线网址| 亚洲国产日韩欧美精品在线观看 | 国产激情偷乱视频一区二区| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 欧美三级亚洲精品| av视频在线观看入口| 人人妻人人看人人澡| 岛国视频午夜一区免费看| 老司机午夜福利在线观看视频| 天天躁夜夜躁狠狠躁躁| 男插女下体视频免费在线播放| 天天一区二区日本电影三级| 国产三级黄色录像| 国产高清videossex| 色噜噜av男人的天堂激情| 国内精品久久久久久久电影| 国产99白浆流出| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 久久久水蜜桃国产精品网| 成人18禁在线播放| 精品久久久久久,| 国产精华一区二区三区| 激情在线观看视频在线高清| 久久中文字幕人妻熟女| 中文字幕久久专区| 久久天堂一区二区三区四区| 特大巨黑吊av在线直播| 97碰自拍视频| 成年人黄色毛片网站| 夜夜看夜夜爽夜夜摸| 国产精品九九99| 最新美女视频免费是黄的| 亚洲片人在线观看| 女人被狂操c到高潮| 久久久久亚洲av毛片大全| 久9热在线精品视频| 午夜两性在线视频| 夜夜爽天天搞| 欧美性长视频在线观看| 亚洲九九香蕉| 50天的宝宝边吃奶边哭怎么回事| 久久精品国产清高在天天线| 久久伊人香网站| 日韩精品免费视频一区二区三区| 成人国语在线视频| 国产三级中文精品| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美成人午夜精品| 午夜福利在线观看吧| 一二三四在线观看免费中文在| 岛国在线免费视频观看| 亚洲午夜理论影院| 九色成人免费人妻av| 久久久久久久久免费视频了| 国产精品一及| 欧美国产日韩亚洲一区| 99久久精品热视频| 欧美成人免费av一区二区三区| 午夜影院日韩av| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| 1024手机看黄色片| 身体一侧抽搐| 男插女下体视频免费在线播放| 日本成人三级电影网站| 国产精品久久久久久人妻精品电影| 欧美激情久久久久久爽电影| 51午夜福利影视在线观看| 国产精品亚洲美女久久久| or卡值多少钱| 亚洲午夜理论影院| 90打野战视频偷拍视频| 国产激情偷乱视频一区二区| 日韩三级视频一区二区三区| 中文字幕久久专区| 亚洲熟妇熟女久久| 一级毛片精品| 国产精品日韩av在线免费观看| 精品国内亚洲2022精品成人| 国产成年人精品一区二区| 一本精品99久久精品77| 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 日韩欧美一区二区三区在线观看| 黄色丝袜av网址大全| 在线永久观看黄色视频| 亚洲精品国产一区二区精华液| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 国产高清激情床上av| 国产高清videossex| 宅男免费午夜| 精品一区二区三区av网在线观看| 亚洲男人天堂网一区| 国产99久久九九免费精品| 88av欧美| 精品国产乱子伦一区二区三区| 国产成年人精品一区二区| xxx96com| 天堂动漫精品| 国产精品野战在线观看| 琪琪午夜伦伦电影理论片6080| 精品高清国产在线一区| 午夜福利免费观看在线| 一本综合久久免费| 亚洲精品中文字幕一二三四区| 欧美绝顶高潮抽搐喷水| 最近最新中文字幕大全免费视频| 成年版毛片免费区| 色精品久久人妻99蜜桃| 久久精品91蜜桃| 后天国语完整版免费观看| 亚洲一区二区三区色噜噜| 国产欧美日韩一区二区精品| 国产精品综合久久久久久久免费| 久久久国产成人免费| 大型黄色视频在线免费观看| 亚洲成a人片在线一区二区| 在线视频色国产色| 久久久久国产一级毛片高清牌| 亚洲国产看品久久| 国产久久久一区二区三区| 在线免费观看的www视频| 日日夜夜操网爽| 国产一区二区三区在线臀色熟女| 女人爽到高潮嗷嗷叫在线视频| 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 欧美黑人巨大hd| 亚洲精品在线美女| 麻豆成人午夜福利视频| 久久精品人妻少妇| 久久久久性生活片| 在线视频色国产色| 亚洲国产中文字幕在线视频| 欧美午夜高清在线| 高潮久久久久久久久久久不卡| av有码第一页| 成人永久免费在线观看视频| 色播亚洲综合网| 亚洲人成网站在线播放欧美日韩| 制服诱惑二区| 嫩草影视91久久| 两人在一起打扑克的视频| 日韩精品免费视频一区二区三区| 久久久久精品国产欧美久久久| 波多野结衣高清作品| 国产成人精品久久二区二区免费| 我要搜黄色片| 婷婷精品国产亚洲av| 日本一二三区视频观看| 久久久久久人人人人人| x7x7x7水蜜桃| 日韩精品免费视频一区二区三区| 久久热在线av| 琪琪午夜伦伦电影理论片6080|