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

    地下水位波動帶三氮遷移轉(zhuǎn)化過程研究進展

    2021-03-19 05:29:14王金生何柱錕
    水文地質(zhì)工程地質(zhì) 2021年2期
    關(guān)鍵詞:硝化波動水位

    劉 鑫,左 銳,王金生,何柱錕,李 橋

    (1.北京師范大學(xué)水科學(xué)研究院,北京 100875;2.地下水污染控制與修復(fù)教育部工程研究中心,北京 100875)

    地下水是我國的主要供水水源,其分布廣泛、供水量穩(wěn)定,在城鄉(xiāng)居民生活用水、農(nóng)業(yè)用水等方面發(fā)揮著重要作用。某一地區(qū)地下水資源量的變化一般與開發(fā)利用程度、豐枯水期、跨區(qū)域調(diào)水、壓采等因素有關(guān),直接表現(xiàn)為地下水位的波動。受天然背景條件和人為污染影響,地下水水質(zhì)在不同區(qū)域表現(xiàn)出明顯的差異性。水位和水質(zhì)是評價地下水資源的重要指標。以往將二者結(jié)合的研究中,主要關(guān)注地下水位的恢復(fù)對生態(tài)系統(tǒng)的改善[1],和地下水水質(zhì)安全的地下水位閾值,側(cè)重于環(huán)境地質(zhì)災(zāi)害的水位問題[2-3]。但是,早期的研究缺乏對于水位變化水質(zhì)響應(yīng)過程的總結(jié)。

    近些年來,地下水位波動過程水質(zhì)的響應(yīng)逐漸得到了一些學(xué)者的關(guān)注。“地下水位波動”、“水位波動帶”等關(guān)鍵詞逐漸出現(xiàn)在諸多的研究中[4-6]。有研究表明,當?shù)叵滤幌陆禃r,包氣帶范圍擴大,氧化作用增強,有機質(zhì)分解,而且包氣帶介質(zhì)特性和水動力條件的變化導(dǎo)致發(fā)生諸如吸附、溶解沉淀等物理化學(xué)作用;地下水位升高時,先前存在于包氣帶的污染物可能溶入地下水中,溶濾作用增強[7]。除此之外,也有學(xué)者探究了水位波動對土壤水分特征規(guī)律影響[5,7],以及關(guān)注地下水中特征組分的遷移轉(zhuǎn)化過程,例如在地下水位波動過程中,研究三氮[6]、鐵錳[4]、砷[8]等特征組分演化規(guī)律及遷移轉(zhuǎn)化過程,從而探究水位波動條件下污染物的歸趨。前人研究沒有針對水位波動帶環(huán)境指標的變化探索規(guī)律。

    作為地下水典型污染物,三氮(硝酸鹽、亞硝酸鹽、銨氮)在很多地區(qū)的土壤、地下水中廣泛存在。三氮參與自然界的氮循環(huán),主要作用包括硝化作用、反硝化作用和異化還原作用等,這些作用需要在微生物的介導(dǎo)下完成[9-11]。這方面的研究已經(jīng)有較好的總結(jié)和綜述,近幾年來,出現(xiàn)了針對地下水位波動帶開展三氮遷移轉(zhuǎn)化的研究[5-6]。地下水位波動引起波動帶土-水環(huán)境要素改變,對三氮的轉(zhuǎn)化過程影響顯著。水位波動也可能使包氣帶中的三氮浸溶進入地下水中,并在淺層地下水中遷移,進而導(dǎo)致地下水中的三氮污染。目前,對于水位波動帶三氮的遷移轉(zhuǎn)化過程和機理尚不明確,缺乏深入研究。因此,探究地下水位波動帶三氮的遷移轉(zhuǎn)化過程及機理具有重要意義。

    本文在綜述國內(nèi)外有關(guān)文獻的基礎(chǔ)上,就地下水位波動過程及其影響、水位波動帶三氮的遷移轉(zhuǎn)化過程和生物地球化學(xué)響應(yīng)分別進行總結(jié),以期為更深入的探討水位波動帶三氮遷移轉(zhuǎn)化的機理研究提供參考,也為地下水水位閾值管理與三氮污染控制提供科學(xué)依據(jù)。

    1 地下水位波動過程

    1.1 地下水位波動過程及影響因素

    地下水位隨氣候和人類活動的影響而上下波動,在最高水位和最低水位之間的區(qū)域稱地下水位波動帶[12]。地下水位波動包括水位上升、下降兩個過程,是地下水資源量增加或減少最直接的表現(xiàn)形式。它反映了地下水補給、徑流、排泄的動態(tài)過程。造成地下水位波動的原因有很多,包括自然條件下,豐水期的降雨補給、潮汐過程、客水、地表水補給地下水過程等,以及人類活動影響下,人工開采、跨區(qū)域調(diào)水、回灌補給、農(nóng)業(yè)灌溉、水資源保護工程和水利工程建設(shè)等[11,13-15]。例如,2014年以來,受潮白河道人工回灌技術(shù)實施、南水北調(diào)中線工程回補等影響,北京地區(qū)地下水位由整體下降轉(zhuǎn)變?yōu)椴糠值貐^(qū)上升,以懷柔應(yīng)急水源地為例,其地下水位上升最大幅度超過11 m[14]。在集中開采的地下水水源地,受到間歇性開采影響,開采期開采井附近地下水位隨之下降,停采期地下水位逐步恢復(fù),地下水位反復(fù)下降與回升,形成地下水位大幅波動[4]。洱海受豐水期降水補給和側(cè)向徑流補給,近岸區(qū)地下水位呈現(xiàn)明顯的季節(jié)性波動,最大變幅達到20 m[15]。王玉剛等[16]通過研究三工河流域尾閭綠洲水資源動態(tài)特征,發(fā)現(xiàn)區(qū)域輸水灌溉是地下水水位抬升的直接動因。

    1.2 地下水位波動帶關(guān)鍵環(huán)境要素變化

    目前,地下水位波動研究主要關(guān)注第四系松散孔隙地層,包氣帶巖性以砂土、粉土、黏土為主,多位于平原區(qū)、河漫灘、階地以及山前沖洪積扇等。埋深較淺的潛水含水層與地表河流水力聯(lián)系密切,同時受大氣降水等補給,潛水位呈現(xiàn)規(guī)律性的波動,形成地下水位波動帶。

    地下水位波動帶水動力過程不斷變化。在地下水位不斷上升和下降的過程中,潛水面逐漸上升和下移,水位標高不斷增加和下降,可能導(dǎo)致區(qū)域含水層水力梯度改變,地下水交換速率增加[17]。相比于水位靜態(tài)平衡,水位波動可能攜帶更多的包氣帶物質(zhì)溶解到地下水中,從而改變地下水的化學(xué)組成。長期地下水位波動過程,會導(dǎo)致波動帶內(nèi)氧化還原條件顯著變化。在水位靜態(tài)平衡時,包氣帶至含水層垂向剖面通常從氧化到逐漸還原,并隨地下水位波動而變化。在地下水位上升時,波動帶由氧化環(huán)境變?yōu)檫€原環(huán)境;在地下水位下降時,波動帶由還原環(huán)境變?yōu)檠趸h(huán)境[18]。地下水位上升將飽和帶的孔隙水向上攜帶,從而將非飽和帶內(nèi)的空氣截留,在靜水壓力作用下,溶解到地下水中,使地下水中溶解氧增加[17]。反之,地下水位下降會導(dǎo)致之前已經(jīng)飽和的土壤介質(zhì)暴露在孔隙空氣中,使得土壤空隙中溶解氧增加而趨于氧化環(huán)境[19]。Kohfahl 等[17]研究水位波動帶氧氣通量時發(fā)現(xiàn),一定幅度的水位波動是孔隙潛水含水層充氧的重要途徑,且增大地下水位波動幅度,波動帶氧化層厚度也會隨之擴大。Bohlke 等[20]探究了美國安諾卡沙平原地下水位回升后,包氣帶水巖作用過程,包氣帶厚度1~4 m,巖性以中砂為主。水位回升導(dǎo)致包氣帶水巖作用范圍增加,沿著水流路徑由淺至深,氧化還原電位(Eh)降低。隨著Eh 進一步降低,厭氧菌厭氧分解有機物并釋放出有機酸引起pH 降低。探究水巖相互作用對 NO-3和 SO24-分布的影響,認為FeS2和其他含亞鐵相的離子是反硝化的主要控制因素。

    當包氣帶存在細顆粒介質(zhì)(如黏性土)時,水位波動過程影響著波動帶中的吸附-解吸以及微生物功能[21-22]。在低水位期,地下水中污染物的稀釋、吸附占主導(dǎo)作用,同時,地下水位下降引起水位波動帶通氣量增加,促進微生物降解作用。在高水位期,水位上升增加地下水與包氣帶中污染物的接觸范圍,以介質(zhì)中污染物的解吸作用為主,污染物再次進入地下水,同時形成厭氧環(huán)境,好氧微生物作用受到抑制[21]。當前研究表明,水位波動導(dǎo)致氧化還原條件改變會影響波動帶的生物地球化學(xué)和微生物動力學(xué)特征,波動帶內(nèi)微生物群落具有更大的功能多樣性和靈活性,可以適應(yīng)不斷變化的含水飽和度和氧化還原電位[18]。水位的波動有利于有機質(zhì)的遷移,以及土壤中營養(yǎng)物質(zhì)和氣體的擴散,波動帶有機質(zhì)得到補給,所以環(huán)境中微生物的營養(yǎng)供給更充足[23]。同時,由于土壤微生物活動與生物氣體(CO2、CH4、N2O 等)的動力學(xué)密切相關(guān),地下水位的波動也可能會調(diào)節(jié)與大氣的氣體交換[24]。水位波動過程也影響著波動帶土壤含水率、基質(zhì)勢,使其呈現(xiàn)周期性的改變[18]。最后,地下水位的改變和氧化還原條件的相關(guān)變化可能會使波動帶具有獨特的地球化學(xué)和礦物學(xué)特征,如綠銹族混合價鐵礦物的存在[18]。

    因此,地下水位波動過程中水動力條件和生物地球化學(xué)特征的不斷響應(yīng),波動帶介質(zhì)中土壤含水率、氧化還原電位、氧氣和有機質(zhì)含量等表現(xiàn)出一定的分帶規(guī)律,沿著水流路徑會形成特定的物理、化學(xué)梯度,進而形成特定的生物群落分帶格局,并發(fā)生著不同的生物地球化學(xué)過程。同時,水位波動影響著地下水水化學(xué)作用,如吸附、解吸作用的發(fā)生,從而控制污染物的遷移轉(zhuǎn)化。

    2 地下水位波動帶三氮的遷移轉(zhuǎn)化過程

    2.1 地下水中的“三氮”污染特征

    作為地下水中的典型污染物,三氮的來源、遷移及轉(zhuǎn)化規(guī)律已有大量研究[6,11,25-27]。地下水中氮的來源主要包括土壤有機氮礦化、地表富營養(yǎng)水體補給、農(nóng)業(yè)氮施肥、生活污水排放等。 NH+4由于帶正電荷,容易被帶負電荷的土壤顆粒吸附,控制其在地下水中的遷移過程。 NO-3由于水溶性強,吸附能力弱,很容易通過大氣降水或灌溉補給,直接發(fā)生垂直遷移進入地下水中,造成淺層地下水污染。一些實驗和現(xiàn)場研究表明,NH+4進入土壤后,在包氣帶土層的淺表層被迅速吸附而減少,在不同土壤剖面上 NH+4的起始濃度隨土壤剖面深度的增加而下降,NO-3在地下水中含量迅速增加[26-27]。地下水中的 NO-2極不穩(wěn)定,容易被氧化成為 NO-3。地下水中 NO-2含量總是小于 NO-3、NH+4的含量,但由于其化學(xué)性質(zhì)和環(huán)境毒性大,常被視為重要的氮污染標志[26]。自然界中氮循環(huán)過程見圖1,三氮參與的轉(zhuǎn)化過程包括硝化作用、反硝化作用和異化(同化)還原作用等。其他的過程也影響著氮的轉(zhuǎn)化,例如有機氮礦化作用、厭氧氨氧化作用等。三氮的轉(zhuǎn)化受地下水化學(xué)組分、微生物群落結(jié)構(gòu)、氮素種類及其濃度、溶解氧、有機碳含量等因素的影響,同時還受到土壤介質(zhì)物理性質(zhì),地下水流速與滯留時間等因素間接影響[11]。

    2.2 地下水位波動帶三氮的遷移和轉(zhuǎn)化過程

    隨著地下水位波動,受波動帶內(nèi)氮素原始累積量、波動帶介質(zhì)顆粒大小、水動力條件影響,三氮呈現(xiàn)不同的遷移特征[28]。水位上升對包氣帶污染物存在強烈的浸溶過程,有利于在地下水中的側(cè)向遷移。Rasiah 等[28]連續(xù)2 個豐枯水期對加拿大Mulgrave 河流域39 個監(jiān)測井中硝酸鹽、氨氮監(jiān)測發(fā)現(xiàn),地下水中硝態(tài)氮和氨氮的濃度隨著雨季地下水位增加而增加,在雨季高峰期出現(xiàn)波動,在雨停后隨水位下降而下降,高水位期濃度可增加至低水位期的近50 倍,且氮的浸出量隨氮肥施加量的增加而增加。豐水期地下水對地表河流側(cè)向補給時,豐水期地表水中52%的氮污染增長來自于地下水[28]。這與zhang[29]室內(nèi)試驗的結(jié)果一致,她認為水位波動促進介質(zhì)中可溶性氮素的解吸和溶濾,使其更充分的溶解進入地下水中,并隨著水位降低過程地下水不斷向外遷移。Huang 等[25]研究得到了同樣的結(jié)論,并認為過量施用氮肥,波動帶中的氮肥積累已經(jīng)造成了西南地區(qū)地下水中氮素含量超標,并存在氮素進一步淋溶污染的潛在風險。由此可見,當包氣帶存在三氮的累積性污染,降雨、灌溉、人工補給等導(dǎo)致地下水位波動時,三氮存在向地下水中遷移的可能,且隨著地下水的側(cè)向徑流,遷移進入其他水體。水位波動帶中氮的積累量和補給水源的入滲速率及入滲量是直接影響三氮遷移過程的關(guān)鍵因素,垂直入滲補給和蒸發(fā)蒸騰速率控制的地下水位波動影響著地下水的氮負荷[30]。水位波動帶土壤孔隙度、黏粒含量等也是三氮遷移的影響因素[30-31]。

    波動帶內(nèi)環(huán)境指標呈現(xiàn)規(guī)律性變化,同時刺激微生物活動,從而影響三氮在水位波動帶內(nèi)的轉(zhuǎn)化行為。目前水位波動對污染物轉(zhuǎn)化過程影響的研究多聚焦于對氧化還原環(huán)境十分敏感的組分,如鐵、錳、砷等的轉(zhuǎn)化規(guī)律[4,8],針對三氮研究并不多見。對于水位波動帶三氮轉(zhuǎn)化過程的研究集中在水位波動過程三氮含量的變化規(guī)律和與環(huán)境指標相關(guān)性分析,從微生物群落結(jié)構(gòu)和功能基因豐度層面探究其機理。如李翔[6]通過開展土柱對比實驗,探討不同水位波動幅度下氮素運移規(guī)律。結(jié)果表明,水位上升階段,硝化作用減弱,反硝化作用增強,土柱內(nèi)硝態(tài)氮濃度減小,銨態(tài)氮濃度增大;水位下降階段,土柱內(nèi)硝化作用增強,NO-3濃度隨之增大。這與楊洋[5]室內(nèi)試驗的結(jié)果一致。水位變化幅度較大,NO-3的變化幅度也較大。

    氧化還原電位是影響反硝化作用的主要因素。潛水位抬升將導(dǎo)致厭氧環(huán)境范圍擴大,有利于反硝化作用進行,最終導(dǎo)致氨氮濃度增加[31]。杜濤[32]認為地下水位回升后,地下水中總氮含量將會一定幅度的降低,有利于降低地下水硝酸鹽污染程度。Li 等[27]利用實測資料研究了地下水位波動對地下水三氮含量的影響,研究表明,氧化還原電位和硝態(tài)氮含量與地下水埋深有較強的相關(guān)性。硝態(tài)氮與銨態(tài)氮濃度呈顯著的負相關(guān),pH 值與硝態(tài)氮濃度也呈顯著負相關(guān)。地下水位變化影響著土壤水分的理化性質(zhì),并進一步對三氮在土壤中的遷移產(chǎn)生顯著影響。Zhang 等[29]通過室內(nèi)砂箱試驗探究水位波動帶三氮變化規(guī)律及微生物的響應(yīng)過程,結(jié)果表明,隨著地下水位下降和土壤剖面持續(xù)干燥,波動帶內(nèi)反硝化作用減弱,硝化作用增強,波動帶內(nèi)三氮含量變化與微生物群落結(jié)構(gòu)及功能基因顯著相關(guān)。

    綜上,地下水位波動促進波動帶內(nèi)三氮浸溶進入地下水,且有利于其在地下水中的遷移,存在造成地下水水質(zhì)污染的可能。地下水位波動影響三氮轉(zhuǎn)化過程,總體趨勢:水位上升,硝化作用減弱,反硝化作用增強;水位下降,硝化作用增強,反硝化作用減弱。

    3 地下水位波動帶三氮的生物地球化學(xué)過程

    3.1 三氮轉(zhuǎn)化的微生物驅(qū)動

    微生物是三氮生物地球化學(xué)過程的重要媒介,對水位波動帶三氮遷移轉(zhuǎn)化起著關(guān)鍵性的作用。近些年來,研究自然環(huán)境(土壤、地表水、地下水)中參與氮循環(huán)的功能微生物種類和多樣性是環(huán)境微生物領(lǐng)域的研究熱點[33-34]。近十年來,分子生物學(xué)檢測技術(shù)得到長足發(fā)展。除了PCR 擴增技術(shù)、變性梯度凝膠電泳技術(shù)(DGGE)、末端限制性片段長度多態(tài)性分析技術(shù)等,基因芯片、高通量測序技術(shù)也廣泛用于三氮轉(zhuǎn)化微生物的分子檢測研究中。隨著生物地球化學(xué)循環(huán)研究方法與分子生物學(xué)檢測技術(shù)的廣泛結(jié)合,氮循環(huán)的機制和過程認識更加深入。三氮參與氮循環(huán)的多個關(guān)鍵過程,包括硝化作用、反硝化作用、異化還原作用及同化還原作用等,參與作用的微生物功能基因見圖1,具體功能基因信息見表1。KEGG(Kyoto Encyclopedia of Genes and Genomes,http://www.genome.jp/kegg/)基因數(shù)據(jù)庫是系統(tǒng)分析基因功能,聯(lián)系基因組信息和功能信息的大型知識庫,KEGG 數(shù)據(jù)庫注釋下的氮循環(huán)通路為ko00910。

    3.1.1 硝化作用

    硝化過程是指 NH+4氧化為 NO-3的過程(NH+4→NH2OH→ NO-2→ NO-3)。多數(shù)研究表明,硝化作用需要分步完成,包括將 NH+4氧化為 NO-2的亞硝化過程和NO-2氧 化為 NO-3的硝化過程[13],參與這兩個過程的氨氧化菌和亞硝酸鹽氧化菌統(tǒng)稱為硝化作用菌,常見的硝化作用菌包括亞硝酸菌屬(Nitrosomonas)、紅球菌屬(Rhodococcus)、假單胞菌屬(Pseudomonas)等。氨氧化是硝化作用的關(guān)鍵步驟,主要由氨氧化古菌(AOA)和氨氧化細菌(AOB)參與完成,氨單加氧酶基因(amoA)是氨氧化過程的標志基因。因此,amoA基因是硝化作用中最受關(guān)注的基因,常用于分析氨氧化古菌(AOA)和氨氧化細菌(AOB)的豐度和多樣性,進而標記硝化作用進程。目前研究表明,AOA、AOB 廣泛存在于土壤環(huán)境中,amoA基因和AOA/AOB 比值隨土壤深度增加而增大[35]。氨氮的濃度、氧氣含量、以及pH 值等環(huán)境因素會影響AOA 和AOB 在氨氧化過程的貢獻,從而導(dǎo)致兩種氨氧化菌的生態(tài)位分離[36-37]。Gurbry-Rangin等[37]利用高通量測序?qū)τ?7 個土樣的氨氧化古菌(AOA)的多樣性進行分析,發(fā)現(xiàn)氨氧化古菌(AOA)的多樣性隨pH 增加顯著增加,amoA豐度隨pH 變化而變化;氨氧化古菌(AOA)更適宜生存在低氮、強酸和高溫條件,氨氧化細菌(AOB)更易存在于高氮、中性和堿性條件。由此可見環(huán)境因素能夠主導(dǎo)AOA、AOB 的功能,進而影響硝化作用過程。此外,土壤濕度、總有機碳、總氮含量對 AOA 和 AOB 均有顯著影響[38]。除amoA基因外,羥胺氧化還原酶基因hao、亞硝酸鹽氧化還原酶基因nxrA也是參與硝化作用的重要功能基因。

    3.1.2 反硝化作用

    硝化-反硝化作用是耦合發(fā)生的過程,也是控制三氮遷移轉(zhuǎn)化的關(guān)鍵作用。反硝化作用是將 NO-3還原成N2的過程(NO-3→ NO-2→NO→N2O→N2)。反硝化作用是一個比較復(fù)雜的兼性厭氧過程,普遍認為細菌是反硝化優(yōu)勢種群。常見的反硝化功能菌多為兼性厭氧化能異氧菌,主要包括芽孢桿菌屬(Bacillus)、假單胞菌屬(Pseudomonas)、產(chǎn)堿桿菌屬(Alcalogenes)、微球菌(Micrococcus)、副球菌屬(Paracoccus)等。

    硝酸鹽還原酶控制著 NO-3還原為 NO-2的過程,標記基因narG是硝酸鹽還原酶的標志性基因。 NO-2還原為NO 的過程是反硝化作用的標志性反應(yīng),相應(yīng)的亞硝酸鹽還原酶是執(zhí)行該步驟的限速酶,其常用的標志基因有nirK和nirS,其豐度常被用來表征及預(yù)測N2O 的釋放[39]。氧化亞氮還原酶控制著N2O 還原為N2的過程,是反硝化作用的最后關(guān)鍵環(huán)節(jié),此酶的標志基因nosZ是研究反硝化微生物的功能標志基因[33]。傳統(tǒng)觀點認為反硝化作用只發(fā)生在無氧或微氧條件下。但20世紀80年代研究發(fā)現(xiàn),好氧反硝化菌具有周質(zhì)型還原酶,該酶活性對于氧氣不敏感,可以驅(qū)動反硝化作用在好氧環(huán)境中進行,至今,好氧反硝化菌僅在細菌中發(fā)現(xiàn)[40]。Mosier 等[41]分析了美國舊金山灣沉積物中含nirK和nirS的反硝化微生物空間變化特征,結(jié)果表明沉積物中反硝化能力與nirS基因的豐度高度相關(guān)。Huang 等[42]分析珠江水體反硝化相關(guān)的功能基因,發(fā)現(xiàn)nirS與nirK主要聚集在水體深度0~10 cm 處,nosZ分布于深度5~15 cm 處,narG集中分布在15 cm以下。在高硝酸鹽濃度影響下,narG和nirS的豐度較高;在貧營養(yǎng)環(huán)境中nosZ具有更高的活性。土壤沉積物中無機氮的含量(包括 NO-3、NO-2和NH+4)、有機質(zhì)的含量、溶解氧和氧化還原電位都將影響反硝化相關(guān)基因的分布。

    表1 參與三氮轉(zhuǎn)化的主要功能基因信息Table 1 Main functional gene information of the three-nitrogen transformation

    3.1.3 異化還原作用與同化還原作用

    異化還原(DRNA)和同化還原作用都指將 NO-3還原成 NH+4的過程(NO-3→ NO-2→ NH+4)。有研究表明與發(fā)酵有關(guān)的梭狀芽孢桿菌(Clostridia)、脫硫弧菌(Desulfovibrio)、假單胞菌(Pseudomonas)在硝酸鹽缺乏條件下都可以進行DRNA 過程。同化還原與異化還原作用的不同之處在于同化還原作用需要NAD(P)H 作為電子供體,且存在有機氮中間體[43]。由于傳統(tǒng)分析方法很難辨別異化還原過程的產(chǎn)物從而識別這一過程,且在厭氧條件下,反硝化作用更容易發(fā)生,所以對異化還原過程關(guān)注度較低。近些年來,通過識別功能基因nrfA(亞硝酸鹽異化還原酶標記基因)可以較好地標記異化還原過程。與nirB相比,nrfA基因在微生物中更加廣泛的存在。有研究表明,nrfA基因相對于nirB基因?qū)Φ蜐舛?NO-3更敏感[44]。由于DNRA過程需要分解有機質(zhì),在碳受限情況下,DNRA 過程會被抑制。隨土層從上而下,由氧化環(huán)境變?yōu)檫€原環(huán)境,nrfA的分布受 NH+4濃度影響,nrfA基因豐度減小,但相對豐度增加[42]。有研究表明,DRNA 反應(yīng)速率隨溫度升高而降低,且其反應(yīng)過程與pH 具有一定關(guān)系[45]。硝酸鹽同化還原酶的標記基因為narB基因。

    3.2 水位波動帶三氮轉(zhuǎn)化的微生物響應(yīng)

    目前針對于水位波動帶三氮轉(zhuǎn)化相關(guān)微生物的研究不是很多,且主要聚焦于硝化和反硝化作用的相關(guān)微生物功能菌和功能基因,對于其他過程研究較少。趙磊等[12]研究了農(nóng)田水位波動帶中氮循環(huán)功能菌群分布特征,發(fā)現(xiàn)在波動帶中上層,由于施肥影響氨氮含量較高,氧氣含量豐富,硝化菌群占絕對優(yōu)勢,利于發(fā)生硝化作用;硝態(tài)氮不易被土壤吸附,容易向下遷移聚集在波動帶下層,硝態(tài)氮增加以及波動帶下層偏還原的條件致使反硝化細菌豐富,反硝化作用加劇。波動帶與異化還原作用相關(guān)的NrfA基因也在厭氧環(huán)境中比在好氧環(huán)境中豐度更大[11]。Liu 等[36]研究表明,氧化還原電位與AOA、AOB 豐度呈顯著的正相關(guān)關(guān)系,水位周期性上升,土壤中氧化還原電位降低,AOA 和AOB 的amoA基因拷貝數(shù)均有所下降,這可能會減弱波動帶內(nèi)的氨氧化過程,從而削弱硝化作用。且與AOB 相比,AOA 在水位波動帶具有更強的適應(yīng)能力,尤其是在低氧條件下具有更強的耐受性。這與Leininger 等[35]研究結(jié)果一致。在深度0~30 cm的農(nóng)業(yè)土壤中,AOA 與AOB 所占比值在1 000 以下;在30~40 cm 缺氧時,比值上升到3 000 以上[35]。崔榮陽等[46-47]模擬地下水位升降過程,分析不同水位條件下土壤剖面AOA-amoA、AOB-amoA、nirK、nirS、nosZ基因豐度的變化特征,研究結(jié)果表明,水位波動對以上功能基因豐度均有顯著影響,AOA-amoA、AOBamoA基因豐度隨著地下水位下降而增加,nirK、nirS、nosZ基因豐度隨著地下水位下降而降低。且AOAamoA和nirS基因?qū)λ簧蹈舾?,分別在水位波動過程硝化與反硝化作用中占主導(dǎo)地位。水位波動過程造成的土壤水分變化也將影響部分功能基因豐度,劉若萱等[48]通過室內(nèi)培養(yǎng)實驗探究不同土壤水分條件對硝化和反硝化的影響。結(jié)果表明,nirS、nirK及AOB-amoA的豐度均隨著水分增加而增加,當nirS、nirK及AOB-amoA的豐度最高時,硝化和反硝化活性也最高,表明這些基因在土壤硝化和反硝化過程中起了重要作用。Zhang 等[29]通過室內(nèi)試驗探究了淺層地下水位波動過程三氮轉(zhuǎn)化功能基因的響應(yīng),結(jié)果顯示樣本中變形菌門和酸桿菌門是波動帶的優(yōu)勢菌門。當?shù)叵滤幱诘臀唬▌訋С掷m(xù)干燥,反硝化功能基因豐度逐漸降低,硝化功能基因逐漸增加。隨著土壤深度增加,微生物功能基因豐度降低,土壤中水分含量、總氮、總磷、氨氮、硝酸鹽以及pH 等化學(xué)指標與波動帶中微生物群落結(jié)構(gòu)顯著相關(guān),nosZ/narG和nosZ/(nirK+nirS)是控制水位波動過程三氮轉(zhuǎn)化的關(guān)鍵基因組。

    綜上,水位波動過程直接影響波動帶內(nèi)氧化還原條件以及化學(xué)組分的再分布,微生物群落結(jié)構(gòu)和功能基因隨之響應(yīng),進而影響三氮的轉(zhuǎn)化。普遍規(guī)律有:波動帶上層,硝化菌群為優(yōu)勢菌群,下層偏還原的條件致使反硝化菌群富集。水位上升過程,與硝化作用相關(guān)的AOA 和AOB 的amoA基因豐度降低,與反硝化作用相關(guān)的nirK、nirS、nosZ基因豐度增加;水位下降過程,AOA 和AOB 的amoA基因豐度增加,nirK、nirS、nosZ基因豐度降低。

    4 結(jié)論與展望

    針對地下水位波動帶三氮遷移轉(zhuǎn)化過程進行系統(tǒng)綜述,對地下水位波動帶關(guān)鍵環(huán)境要素變化特征、三氮的遷移轉(zhuǎn)化規(guī)律以及微生物功能菌群和功能基因響應(yīng)等方面的研究進行了歸納和總結(jié)?,F(xiàn)有研究表明:(1)水位波動帶中環(huán)境指標如土壤含水率、氧化還原電位、溶解氧和有機質(zhì)含量等均表現(xiàn)出一定的分帶性規(guī)律。(2)隨著地下水位波動,包氣帶中的三氮易浸溶進入地下水并發(fā)生遷移。且地下水位上升,硝化作用減弱,反硝化作用增強,地下水位下降,硝化作用增強,反硝化作用減弱。(3)波動帶三氮轉(zhuǎn)化相關(guān)的微生物菌群結(jié)構(gòu)和功能基因更多樣化,呈現(xiàn)一定的分布特征。如波動帶上層,硝化菌群為優(yōu)勢菌群,下層偏還原的條件致使反硝化菌群富集等。

    目前,雖然國內(nèi)外對水位波動帶三氮遷移轉(zhuǎn)化過程已經(jīng)取得了一定的研究成果,但主要還是聚焦于某些化學(xué)指標(如氧化還原電位、土壤含水率、pH 值等)含量、微生物功能基因變化以及多因素間相關(guān)性展開討論,對于水位波動條件下三氮遷移轉(zhuǎn)化的機理探討仍較為薄弱。且目前的研究主要針對硝化作用和反硝化作用,對于其他的三氮轉(zhuǎn)化作用聚焦較少,參與作用的微生物標志性功能基因的研究也較少。為完善水位波動帶三氮遷移轉(zhuǎn)化過程研究,進一步研究中需要關(guān)注:(1)水位波動帶三氮轉(zhuǎn)化及與微生物作用機理。將水化學(xué)演化分析與分子生物學(xué)高通量測序技術(shù)相結(jié)合,在了解波動帶環(huán)境指標變化的前提下,深入研究波動帶三氮轉(zhuǎn)化機理及微生物功能基因響應(yīng),探究影響三氮轉(zhuǎn)化的關(guān)鍵功能基因。(2)除關(guān)注硝化、反硝化作用外,增加異化還原、同化還原和厭氧氨氧化等作用過程的研究,擴展補充波動帶三氮轉(zhuǎn)化的更多過程,同時識別以上過程的微生物標志性功能基因。(3)探究不同情境(降雨、客水補給、人工開采)、不同影響因素(波動幅度、波動周期)的水位波動過程。細化分析不同情境水位波動影響的三氮遷移轉(zhuǎn)化,對比不同波動幅度、波動周期對三氮轉(zhuǎn)化的影響,識別水位波動帶三氮轉(zhuǎn)化過程的關(guān)鍵影響要素。筆者認為,深化認識水位波動帶這一特殊的環(huán)境條件下三氮的遷移轉(zhuǎn)化過程,可以為地下水位閾值管理提供科學(xué)依據(jù),同時可以有效緩解由于水資源不合理開采或者補給引起的地下水中的三氮污染問題。

    猜你喜歡
    硝化波動水位
    羊肉價回穩(wěn) 后期不會大幅波動
    微風里優(yōu)美地波動
    中國化肥信息(2019年3期)2019-04-25 01:56:16
    MBBR中進水有機負荷對短程硝化反硝化的影響
    干濕法SO2排放波動對比及分析
    厭氧氨氧化與反硝化耦合脫氮除碳研究Ⅰ:
    基于MFAC-PID的核電站蒸汽發(fā)生器水位控制
    海水反硝化和厭氧氨氧化速率同步測定的15N示蹤法及其應(yīng)用
    基于PLC的水位控制系統(tǒng)的設(shè)計與研究
    河南科技(2014年4期)2014-02-27 14:07:11
    短程硝化反硝化脫氮技術(shù)的研究進展
    河南科技(2010年1期)2010-10-26 07:53:50
    久久人人爽人人片av| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 亚洲国产日韩一区二区| av在线app专区| 97精品久久久久久久久久精品| 一二三四中文在线观看免费高清| 久久国产亚洲av麻豆专区| 精品久久久久久电影网| 美女主播在线视频| 男女啪啪激烈高潮av片| 中文字幕人妻丝袜制服| 亚洲在久久综合| 国产淫语在线视频| 久久久久久人妻| 国产成人aa在线观看| 嫩草影院入口| 精品少妇久久久久久888优播| 又大又黄又爽视频免费| 日产精品乱码卡一卡2卡三| 叶爱在线成人免费视频播放| 深夜精品福利| 亚洲三级黄色毛片| 极品少妇高潮喷水抽搐| 热re99久久精品国产66热6| 热99国产精品久久久久久7| 国产视频首页在线观看| av又黄又爽大尺度在线免费看| 国产 一区精品| 欧美精品一区二区大全| 老司机影院毛片| 欧美变态另类bdsm刘玥| 涩涩av久久男人的天堂| 在线观看www视频免费| 春色校园在线视频观看| 亚洲欧洲精品一区二区精品久久久 | 中文字幕人妻熟女乱码| 欧美精品一区二区大全| 国产av精品麻豆| 亚洲精品美女久久久久99蜜臀 | av片东京热男人的天堂| 人人妻人人澡人人看| 日产精品乱码卡一卡2卡三| 丝袜美腿诱惑在线| 亚洲成av片中文字幕在线观看 | 免费大片黄手机在线观看| 日日啪夜夜爽| 热99国产精品久久久久久7| 亚洲,一卡二卡三卡| 日日爽夜夜爽网站| 少妇猛男粗大的猛烈进出视频| 精品一区在线观看国产| 国产亚洲午夜精品一区二区久久| 日韩一卡2卡3卡4卡2021年| 视频在线观看一区二区三区| 岛国毛片在线播放| 亚洲成av片中文字幕在线观看 | 永久网站在线| 日韩熟女老妇一区二区性免费视频| 精品视频人人做人人爽| www.自偷自拍.com| 十八禁高潮呻吟视频| 免费看av在线观看网站| 少妇熟女欧美另类| 高清欧美精品videossex| 久久午夜福利片| 校园人妻丝袜中文字幕| 男女边吃奶边做爰视频| 色吧在线观看| 午夜av观看不卡| 一级片免费观看大全| 一级毛片黄色毛片免费观看视频| 黄色一级大片看看| 在线天堂中文资源库| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费看不卡的av| 久久国产精品男人的天堂亚洲| 美女视频免费永久观看网站| 97人妻天天添夜夜摸| 一级a爱视频在线免费观看| 国产精品熟女久久久久浪| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 91精品三级在线观看| 国产av一区二区精品久久| 巨乳人妻的诱惑在线观看| 免费黄频网站在线观看国产| 中文字幕制服av| 18禁观看日本| 久久久精品区二区三区| 最近2019中文字幕mv第一页| 女性被躁到高潮视频| 在现免费观看毛片| 亚洲欧美一区二区三区久久| 午夜福利视频精品| 中文字幕av电影在线播放| 国产成人免费无遮挡视频| 久久99热这里只频精品6学生| 欧美另类一区| av国产精品久久久久影院| 亚洲精品国产av成人精品| 中文字幕精品免费在线观看视频| 婷婷色综合www| 99久久精品国产国产毛片| 热99国产精品久久久久久7| 午夜福利影视在线免费观看| 丰满少妇做爰视频| 亚洲国产最新在线播放| 久久国产亚洲av麻豆专区| av在线播放精品| 天堂8中文在线网| 咕卡用的链子| 一本久久精品| 日韩欧美精品免费久久| 亚洲国产最新在线播放| 五月天丁香电影| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区国产| 97精品久久久久久久久久精品| 97精品久久久久久久久久精品| 青春草国产在线视频| 在线免费观看不下载黄p国产| 欧美日韩精品网址| 久久久久精品性色| www日本在线高清视频| 亚洲四区av| 一级毛片我不卡| 亚洲精品国产av成人精品| 久久国内精品自在自线图片| 韩国av在线不卡| 亚洲五月色婷婷综合| 黄网站色视频无遮挡免费观看| 国产深夜福利视频在线观看| 少妇精品久久久久久久| 欧美精品高潮呻吟av久久| 国产xxxxx性猛交| 国产xxxxx性猛交| a级片在线免费高清观看视频| 另类亚洲欧美激情| 国产精品二区激情视频| 国产白丝娇喘喷水9色精品| 免费观看在线日韩| 你懂的网址亚洲精品在线观看| 制服丝袜香蕉在线| 母亲3免费完整高清在线观看 | 又大又黄又爽视频免费| 日韩一卡2卡3卡4卡2021年| 18禁裸乳无遮挡动漫免费视频| 久久久久国产一级毛片高清牌| 国产综合精华液| 亚洲国产欧美网| 成年人午夜在线观看视频| 久久久久精品性色| 丰满饥渴人妻一区二区三| 黄网站色视频无遮挡免费观看| 欧美成人午夜精品| 丝袜在线中文字幕| 久久久久久久国产电影| 十八禁高潮呻吟视频| 一级爰片在线观看| 国产精品香港三级国产av潘金莲 | 大片电影免费在线观看免费| 黑人欧美特级aaaaaa片| 啦啦啦在线免费观看视频4| av电影中文网址| 最近的中文字幕免费完整| 久久久精品区二区三区| 国产xxxxx性猛交| 18+在线观看网站| 伦精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 最近中文字幕2019免费版| 精品人妻偷拍中文字幕| 亚洲一区二区三区欧美精品| 黄片无遮挡物在线观看| 日日撸夜夜添| 欧美另类一区| 久久人人爽av亚洲精品天堂| 亚洲成人手机| 水蜜桃什么品种好| 一区二区三区激情视频| 乱人伦中国视频| 亚洲国产成人一精品久久久| tube8黄色片| 国产不卡av网站在线观看| 五月开心婷婷网| 91国产中文字幕| 国产精品人妻久久久影院| 激情五月婷婷亚洲| 国产av国产精品国产| 日韩免费高清中文字幕av| 免费日韩欧美在线观看| 欧美亚洲日本最大视频资源| 亚洲av在线观看美女高潮| 久久这里有精品视频免费| 人人妻人人爽人人添夜夜欢视频| av有码第一页| 亚洲精品aⅴ在线观看| 亚洲精品视频女| 久久这里只有精品19| 午夜激情久久久久久久| 国产综合精华液| 国产成人av激情在线播放| 一二三四中文在线观看免费高清| 一本久久精品| 日日爽夜夜爽网站| 亚洲精品国产一区二区精华液| 91aial.com中文字幕在线观看| 久久久久国产网址| 亚洲国产毛片av蜜桃av| 国产精品久久久久久精品电影小说| 亚洲情色 制服丝袜| 最黄视频免费看| 欧美xxⅹ黑人| 亚洲成av片中文字幕在线观看 | 免费在线观看视频国产中文字幕亚洲 | 人妻系列 视频| 男人爽女人下面视频在线观看| 一区二区三区乱码不卡18| 日韩成人av中文字幕在线观看| 日韩大片免费观看网站| www.精华液| 亚洲欧美清纯卡通| 成人二区视频| 涩涩av久久男人的天堂| 国产精品二区激情视频| 高清黄色对白视频在线免费看| 日本wwww免费看| 色哟哟·www| 一区二区三区四区激情视频| 亚洲欧美一区二区三区国产| 女人高潮潮喷娇喘18禁视频| 久久鲁丝午夜福利片| 亚洲经典国产精华液单| 18在线观看网站| 国产97色在线日韩免费| 热re99久久国产66热| 又黄又粗又硬又大视频| 最近的中文字幕免费完整| 国产淫语在线视频| 久久久国产精品麻豆| 男女免费视频国产| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 青春草国产在线视频| 亚洲精品久久成人aⅴ小说| 国产成人精品久久久久久| 视频区图区小说| 成人漫画全彩无遮挡| 一本大道久久a久久精品| 另类精品久久| 欧美精品高潮呻吟av久久| 亚洲人成77777在线视频| 搡女人真爽免费视频火全软件| 一区二区三区四区激情视频| 一边摸一边做爽爽视频免费| 亚洲精品一区蜜桃| 亚洲成av片中文字幕在线观看 | 日韩av在线免费看完整版不卡| av免费在线看不卡| av天堂久久9| 国产精品久久久久成人av| 18禁动态无遮挡网站| 久久毛片免费看一区二区三区| 午夜福利在线免费观看网站| 国产极品粉嫩免费观看在线| 免费观看a级毛片全部| 久久免费观看电影| 天堂中文最新版在线下载| 欧美最新免费一区二区三区| 久久女婷五月综合色啪小说| 99精国产麻豆久久婷婷| 欧美日韩一区二区视频在线观看视频在线| 男女下面插进去视频免费观看| 亚洲av电影在线进入| 欧美人与善性xxx| 亚洲在久久综合| 国产一区二区三区综合在线观看| 另类亚洲欧美激情| 久久毛片免费看一区二区三区| 成人午夜精彩视频在线观看| 欧美av亚洲av综合av国产av | 王馨瑶露胸无遮挡在线观看| 在线观看www视频免费| 少妇被粗大的猛进出69影院| 制服丝袜香蕉在线| 免费播放大片免费观看视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 精品视频人人做人人爽| 波野结衣二区三区在线| 久久99蜜桃精品久久| 久久午夜综合久久蜜桃| 国产老妇伦熟女老妇高清| 女人高潮潮喷娇喘18禁视频| 欧美激情极品国产一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲精品乱久久久久久| 97人妻天天添夜夜摸| 亚洲欧美色中文字幕在线| 亚洲成色77777| 少妇被粗大猛烈的视频| 久久影院123| 亚洲三区欧美一区| 亚洲国产精品999| 满18在线观看网站| 大码成人一级视频| 最近手机中文字幕大全| 国产在线一区二区三区精| 久久久欧美国产精品| 国产成人精品久久久久久| 亚洲av男天堂| 久久97久久精品| 超碰97精品在线观看| 日韩电影二区| 老司机影院毛片| 国产精品秋霞免费鲁丝片| 一区二区三区四区激情视频| 国产亚洲一区二区精品| 精品人妻在线不人妻| 在线看a的网站| 中国三级夫妇交换| 啦啦啦在线免费观看视频4| 如何舔出高潮| 亚洲欧美成人精品一区二区| 麻豆av在线久日| 日韩av免费高清视频| 精品少妇内射三级| 高清在线视频一区二区三区| 岛国毛片在线播放| 汤姆久久久久久久影院中文字幕| 男女高潮啪啪啪动态图| 在现免费观看毛片| 久久精品国产a三级三级三级| 90打野战视频偷拍视频| 69精品国产乱码久久久| 国产精品一区二区在线不卡| 国产亚洲一区二区精品| 老熟女久久久| 亚洲第一区二区三区不卡| 亚洲av成人精品一二三区| 校园人妻丝袜中文字幕| 一区二区三区激情视频| 久久久欧美国产精品| 天天躁夜夜躁狠狠久久av| av一本久久久久| 欧美日韩av久久| 国产又色又爽无遮挡免| 成人毛片a级毛片在线播放| 精品久久蜜臀av无| 热99国产精品久久久久久7| 免费少妇av软件| 天堂中文最新版在线下载| 69精品国产乱码久久久| 国产精品女同一区二区软件| 麻豆乱淫一区二区| 久久毛片免费看一区二区三区| 亚洲精品视频女| av天堂久久9| 久久人人爽人人片av| 久久狼人影院| 美女高潮到喷水免费观看| 熟女电影av网| 日韩av在线免费看完整版不卡| 日本午夜av视频| xxx大片免费视频| 好男人视频免费观看在线| 我要看黄色一级片免费的| 久久午夜福利片| 曰老女人黄片| 亚洲成人一二三区av| videos熟女内射| h视频一区二区三区| 中文字幕人妻丝袜一区二区 | 欧美精品亚洲一区二区| 久久久精品区二区三区| 韩国精品一区二区三区| 久久婷婷青草| 亚洲经典国产精华液单| 亚洲精品久久久久久婷婷小说| 亚洲国产av新网站| 亚洲国产最新在线播放| 赤兔流量卡办理| 日本午夜av视频| 婷婷色综合www| 考比视频在线观看| 一区二区av电影网| 亚洲av免费高清在线观看| 欧美97在线视频| 18在线观看网站| av在线播放精品| 国产日韩欧美亚洲二区| 久久久亚洲精品成人影院| 精品少妇久久久久久888优播| 韩国精品一区二区三区| 高清黄色对白视频在线免费看| 伦理电影免费视频| 国产欧美日韩一区二区三区在线| 亚洲一码二码三码区别大吗| 伦理电影免费视频| 婷婷色av中文字幕| 美女高潮到喷水免费观看| 一级毛片 在线播放| 亚洲国产精品成人久久小说| 国产精品麻豆人妻色哟哟久久| 亚洲男人天堂网一区| 热re99久久国产66热| 在线观看www视频免费| av有码第一页| 成人国产av品久久久| 日本欧美视频一区| 男女边摸边吃奶| 久久精品国产综合久久久| 亚洲精品久久成人aⅴ小说| 亚洲精品久久久久久婷婷小说| 黄色怎么调成土黄色| 十八禁高潮呻吟视频| 伊人亚洲综合成人网| 国产日韩一区二区三区精品不卡| 高清欧美精品videossex| 欧美97在线视频| 91aial.com中文字幕在线观看| 春色校园在线视频观看| 另类亚洲欧美激情| 亚洲天堂av无毛| 午夜91福利影院| 日本欧美视频一区| 欧美日韩成人在线一区二区| 9191精品国产免费久久| h视频一区二区三区| 街头女战士在线观看网站| 大片免费播放器 马上看| 熟妇人妻不卡中文字幕| 天天操日日干夜夜撸| 激情五月婷婷亚洲| 久久久久久人妻| 99国产综合亚洲精品| 激情五月婷婷亚洲| 国产野战对白在线观看| 国产极品天堂在线| 卡戴珊不雅视频在线播放| 久久狼人影院| 国产黄频视频在线观看| 人妻人人澡人人爽人人| 国产午夜精品一二区理论片| 日韩伦理黄色片| 曰老女人黄片| 国产欧美亚洲国产| 一级片'在线观看视频| 另类精品久久| 中文字幕另类日韩欧美亚洲嫩草| 国产男女内射视频| 美女视频免费永久观看网站| 亚洲国产精品一区二区三区在线| 97在线视频观看| 超碰97精品在线观看| 黑丝袜美女国产一区| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 欧美精品人与动牲交sv欧美| 色视频在线一区二区三区| 男人添女人高潮全过程视频| 亚洲精品,欧美精品| 人妻人人澡人人爽人人| a级片在线免费高清观看视频| a级片在线免费高清观看视频| 免费黄频网站在线观看国产| 女性生殖器流出的白浆| 国产黄频视频在线观看| 性色avwww在线观看| 97在线人人人人妻| 欧美精品一区二区大全| 久久ye,这里只有精品| 久久久久久久大尺度免费视频| 欧美日韩一区二区视频在线观看视频在线| 伦理电影免费视频| 99久久人妻综合| 午夜福利网站1000一区二区三区| 亚洲国产精品一区三区| 久久精品国产a三级三级三级| 人成视频在线观看免费观看| 啦啦啦中文免费视频观看日本| 极品少妇高潮喷水抽搐| 亚洲一码二码三码区别大吗| 男女边吃奶边做爰视频| 免费看av在线观看网站| 寂寞人妻少妇视频99o| 午夜91福利影院| 久久精品国产亚洲av涩爱| 99久久中文字幕三级久久日本| 在线观看免费高清a一片| 亚洲国产看品久久| 又大又黄又爽视频免费| 国精品久久久久久国模美| 丁香六月天网| 日韩精品有码人妻一区| 免费黄色在线免费观看| 亚洲激情五月婷婷啪啪| 亚洲精品日本国产第一区| 国产精品.久久久| 国产精品一区二区在线不卡| 一个人免费看片子| 国产 一区精品| 水蜜桃什么品种好| 国产av精品麻豆| 啦啦啦在线免费观看视频4| 一级爰片在线观看| 亚洲综合色惰| 美女视频免费永久观看网站| 人成视频在线观看免费观看| 午夜日韩欧美国产| 亚洲,欧美,日韩| 国产精品三级大全| 波多野结衣av一区二区av| 人人妻人人澡人人看| 亚洲av综合色区一区| 天天躁夜夜躁狠狠久久av| 天天躁日日躁夜夜躁夜夜| 捣出白浆h1v1| 国产片内射在线| 国产日韩欧美在线精品| 日本wwww免费看| 男人舔女人的私密视频| 国产精品偷伦视频观看了| 97人妻天天添夜夜摸| 黄频高清免费视频| 午夜日韩欧美国产| 久久国产精品大桥未久av| 啦啦啦中文免费视频观看日本| 亚洲视频免费观看视频| 国产男人的电影天堂91| 男女午夜视频在线观看| 精品国产乱码久久久久久小说| 国产成人a∨麻豆精品| 尾随美女入室| 国产视频首页在线观看| 男人添女人高潮全过程视频| av免费在线看不卡| 男人舔女人的私密视频| 纵有疾风起免费观看全集完整版| 成人国产av品久久久| 午夜福利在线观看免费完整高清在| 久久 成人 亚洲| 亚洲久久久国产精品| 精品少妇久久久久久888优播| 赤兔流量卡办理| 国产男人的电影天堂91| 涩涩av久久男人的天堂| 一区二区av电影网| 精品99又大又爽又粗少妇毛片| 2022亚洲国产成人精品| 国产 精品1| 日韩av不卡免费在线播放| 国产精品99久久99久久久不卡 | 午夜福利在线免费观看网站| 免费观看无遮挡的男女| 国产在线视频一区二区| 久久综合国产亚洲精品| 两个人免费观看高清视频| 啦啦啦在线观看免费高清www| 91成人精品电影| 各种免费的搞黄视频| 新久久久久国产一级毛片| 99热网站在线观看| 丝瓜视频免费看黄片| 国产欧美日韩综合在线一区二区| 婷婷成人精品国产| 中文欧美无线码| 一边摸一边做爽爽视频免费| 日韩一本色道免费dvd| 少妇 在线观看| 国产极品天堂在线| 侵犯人妻中文字幕一二三四区| 91午夜精品亚洲一区二区三区| 热99国产精品久久久久久7| 中文字幕人妻丝袜制服| 一级毛片我不卡| 99国产综合亚洲精品| 99九九在线精品视频| 色网站视频免费| 欧美日韩国产mv在线观看视频| 在线天堂最新版资源| 最近手机中文字幕大全| 国产有黄有色有爽视频| 日本欧美国产在线视频| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| a级毛片在线看网站| 九色亚洲精品在线播放| 麻豆av在线久日| 精品少妇一区二区三区视频日本电影 | 久久久久人妻精品一区果冻| 久久久欧美国产精品| 久久ye,这里只有精品| 免费看不卡的av| 黑人欧美特级aaaaaa片| 精品少妇黑人巨大在线播放| 日本爱情动作片www.在线观看| 国产一区二区三区综合在线观看| 欧美中文综合在线视频| 赤兔流量卡办理| 制服人妻中文乱码| 亚洲欧美成人精品一区二区| 久久影院123| 亚洲精品美女久久av网站| av在线播放精品| 久久久久久久久久久免费av| 王馨瑶露胸无遮挡在线观看| 熟女电影av网| 超碰成人久久| 国产福利在线免费观看视频| 亚洲欧美一区二区三区国产| 国产精品成人在线| 亚洲国产欧美日韩在线播放|