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

    不同根系分泌物對土壤N2O排放及同位素特征值的影響

    2020-06-03 08:02:36莊姍林偉丁軍軍鄭欠寇馨月李巧珍李玉中
    中國農(nóng)業(yè)科學(xué) 2020年9期

    莊姍,林偉,丁軍軍,鄭欠,寇馨月,李巧珍,李玉中,2

    不同根系分泌物對土壤N2O排放及同位素特征值的影響

    莊姍1,林偉1,丁軍軍1,鄭欠1,寇馨月1,李巧珍1,李玉中1,2

    (1中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所/農(nóng)業(yè)農(nóng)村部旱地節(jié)水農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室,北京 100081;2中國農(nóng)業(yè)科學(xué)院環(huán)境穩(wěn)定同位素實(shí)驗(yàn)室,北京 100081)

    【】探究植物根系分泌的主要組分(有機(jī)酸、氨基酸、糖類)對土壤N2O排放及其微生物過程的影響,為選擇適宜的植物進(jìn)而控制土壤N2O排放提供支撐。通過室內(nèi)試驗(yàn)分別添加草酸、絲氨酸、葡萄糖于土壤中模擬根系的3種主要分泌物,每種分泌物設(shè)置兩個(gè)濃度水平:低濃度(150 μg C·d-1)和高濃度(300 μg C·d-1),另設(shè)置添加蒸餾水的對照組,共7個(gè)處理。將土壤置于120 mL玻璃瓶中進(jìn)行培養(yǎng),24 h內(nèi)采集氣體樣品7次,每次培養(yǎng)2 h,獲取N2O排放速率、日累積排放量和同位素特征值(15Nbulk18O和SP(site preference,SP=15Nα-15Nβ))。添加3種根系分泌物組分后,土壤N2O排放速率均逐漸升高,且均高于對照。高濃度處理組N2O累積排放量為:葡萄糖((3.2±1.3)mg·kg-1·d-1)處理>絲氨酸((2.6±0.5)mg·kg-1·d-1)處理>草酸((1.4±0.2)mg·kg-1·d-1)處理,低濃度處理組為:草酸((2.7±1.3)mg·kg-1·d-1)處理>絲氨酸((1.8±0.4)mg·kg-1·d-1)處理>葡萄糖((1.6±0.8)mg·kg-1·d-1)處理;添加根系分泌物的不同處理間土壤N2O的18O值無明顯差異,并穩(wěn)定在24.1‰—25.6‰,且均顯著高于對照((20.1±1.5)‰);土壤N2O的15Nbulk值與添加根系分泌物的種類有關(guān),其中草酸處理組為(-20.06±2.22)‰、絲氨酸處理組為(-22.33±1.10)‰、葡萄糖處理組為(-13.86±1.11)‰、對照組為(-23.14±3.72)‰。各處理土壤N2O的SP值的變化范圍為13.13‰—15.03‰,根系分泌物濃度越高,SP值越低。綜合分析不同處理4個(gè)指標(biāo)(N2O排放速率、N2O的15Nbulk18O和SP值)的不同時(shí)刻的檢測值與日均值的校正系數(shù),添加根系分泌物后第16小時(shí)各處理4個(gè)指標(biāo)的校正系數(shù)最接近于1。在NH+ 4-300 mg N·kg-1的土壤環(huán)境下根系分泌物促進(jìn)N2O的排放,且在培養(yǎng)期間(24 h)土壤N2O排放速率逐漸升高。高濃度處理組葡萄糖對土壤N2O排放速率促進(jìn)效果最強(qiáng),低濃度處理組草酸對土壤N2O排放速率促進(jìn)效果最強(qiáng)。與對照組相比,根系分泌物的添加使N2O的18O值顯著升高;與對照組相比,葡萄糖的添加使15Nbulk值顯著升高。根系分泌物濃度越高,反硝化作用對N2O的貢獻(xiàn)越大。

    根系分泌物;硝化作用;反硝化作用;N2O;同位素特征值

    0 引言

    【研究意義】N2O是一種重要的溫室氣體,能夠破壞平流層中的臭氧層且參與大氣中光化學(xué)反應(yīng),從而進(jìn)一步加劇溫室效應(yīng)[1]。土壤是N2O的主要排放源[1-2],約占全球N2O排放總量的65%[3]。自20世紀(jì)初合成氨技術(shù)發(fā)明以來,氮肥被大量的投入到農(nóng)業(yè)生產(chǎn)中,在增加作物產(chǎn)量的同時(shí)也使得N2O排放量增加了19%[4-6]。而菜地具有高水、高肥、肥料利用率低等特點(diǎn)[7],造成了N2O的大量排放。因此研究菜地N2O排放對減緩全球變暖具有重要的理論和現(xiàn)實(shí)意義?!厩叭搜芯窟M(jìn)展】硝化作用和反硝化作用是N2O產(chǎn)生的主要過程[8]。傳統(tǒng)的研究方法主要有乙炔抑制法和同位素標(biāo)記法等,這些方法在實(shí)際操作過程中都有一些局限性,如抑制劑(乙炔)和標(biāo)記物在系統(tǒng)中擴(kuò)散不均,易受系統(tǒng)干擾等[9-10]。而自然豐度條件下的N2O同位素位嗜值(SP值)不受N2O前體同位素的影響,進(jìn)而克服土壤與N作用的環(huán)境阻礙而直接評價(jià)N2O產(chǎn)生的微生物過程,是估測N2O來源及貢獻(xiàn)的一種有效方法[11]。盡管N2O的15Nbulk和18O值易受底物同位素影響,但也能夠在一定程度上反應(yīng)微生物過程,其與SP值相結(jié)合進(jìn)行分析有助于更加全面了解N2O的排放行為。N2O排放的影響因素很多,其中對土壤孔隙含水量(WFPS)[12-13]、pH[14]、O2濃度[15]、溫度[16-18]等研究較多。然而,植物通過根系分泌物對N2O的產(chǎn)生也具有重要的影響。植物凈光合固定的碳5%—25%通過根系分泌物流失到土壤中[19-21],并以有機(jī)酸,氨基酸,糖類和其他類型分子的形式釋放,是土壤重要的碳源,影響土壤硝化和反硝化作用[22-23];根系分泌物還可以作為電子供體還原NO- 3和NO- 2[24],改變N2O產(chǎn)生的電子傳遞途徑;還可通過改變根際周圍微生物的群落結(jié)構(gòu)影響N2O的產(chǎn)生[25]。草酸、絲氨酸、葡萄糖分別是根系分泌物中普遍存在并且含量較多的一種低分子有機(jī)酸、氨基酸、糖類[26],它們在N2O產(chǎn)生過程中的作用機(jī)理存在一定差異。研究發(fā)現(xiàn)草酸顯著促進(jìn)土壤中有機(jī)態(tài)氮和有機(jī)無機(jī)復(fù)合態(tài)氮的礦化速度,促進(jìn)氮的釋放[27]。也有研究表明草酸促進(jìn)N2O還原基因的復(fù)制,對N2O的還原產(chǎn)生影響[28]。葡萄糖也能夠促進(jìn)氮的分解,但這種促進(jìn)作用明顯小于草酸,此外,優(yōu)先利用葡萄糖的微生物可能導(dǎo)致氮分解率下降[28]。絲氨酸不僅是細(xì)胞蛋白質(zhì)的重要組成部分,也是微生物的能量來源,在土壤中可以礦化為NH+ 4,改變土壤的C/N和pH[29]?!颈狙芯壳腥朦c(diǎn)】不同植物根系分泌物的組成和含量不同,大量研究表明N2O的排放與根系分泌物的組成和含量有關(guān)[24-25,29-30],但鮮有系統(tǒng)研究不同組分單獨(dú)作用對N2O產(chǎn)生過程的影響。又由于根系分泌物組成和土壤環(huán)境的復(fù)雜,在自然條件下很難做到定量提取研究[31],所以通過模擬根系分泌物組分和含量制備根系分泌物溶液添加到土壤中是研究根系分泌物作用下土壤N2O排放特征的重要方法。盡管實(shí)驗(yàn)室環(huán)境能夠控制土壤溫度、濕度等環(huán)境因子以減少環(huán)境差異對N2O排放的影響,但添加根系分泌物后微生物對各根系分泌物的相互作用和響應(yīng)時(shí)間是影響N2O排放量的關(guān)鍵,N2O排放的時(shí)空變異性對定量分析不同根系分泌物組分對N2O產(chǎn)生作用的影響仍無法排除。因此研究不同根系分泌物組分對土壤N2O排放的影響,同時(shí)確定合理的氣體取樣時(shí)間,對于提高試驗(yàn)效率,深入研究根系分泌物對N2O產(chǎn)生在時(shí)間和空間上的作用至關(guān)重要?!緮M解決的關(guān)鍵問題】通過對根系分泌物中3種主要組分(有機(jī)酸、絲氨酸、糖類)的主要成分(草酸、絲氨酸、葡萄糖)作用下土壤N2O排放和同位素特征的綜合分析,揭示不同根系分泌物對土壤N2O產(chǎn)生過程的影響及其作用機(jī)制,為選擇適宜的植物進(jìn)而減緩?fù)寥繬2O排放提供支撐。

    1 材料與方法

    1.1 土壤采集

    土壤采集于中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所順義實(shí)驗(yàn)基地,東經(jīng)40°05′2″,北緯116°55′19″。基地為典型的暖溫帶半濕潤大陸性季風(fēng)氣候,年平均日照2 684 h,年平均氣溫12.5℃,年平均降水量623.5 mm,無霜期195 d左右。土壤類型為潮褐土,土壤密度1.40 g·cm-3,有機(jī)質(zhì)16.48 g·kg-1,全氮1.10 g·kg-1, 全磷0.68 g·kg-1,全鉀20.31 g·kg-1,pH 8.7。

    土壤采集于2018年8月(雨后2 d),取0—20 cm表層土壤于陰涼處自然風(fēng)干。除去砂石、植物根系等雜質(zhì),過2 mm篩,室內(nèi)保存?zhèn)溆谩?/p>

    1.2 試驗(yàn)方法

    選擇3種代表性有機(jī)物,即草酸代表有機(jī)酸,絲氨酸代表氨基酸,葡萄糖代表糖類,以此為研究對象。依據(jù)白菜根系分泌速率[26],設(shè)置2個(gè)添加水平:150 μg C·d-1(低濃度)和300 μg C·d-1(高濃度),另添加蒸餾水為空白對照,共7個(gè)處理,每次添加溶液體積為0.85 mL,每種處理設(shè)3次重復(fù)。用120 mL密封性良好的血清瓶作培養(yǎng)瓶,稱取40 g風(fēng)干土放入120 mL血清瓶中,共21個(gè)。用蒸餾水調(diào)節(jié)土壤孔隙含水量(WFPS:water-filled-pore-space)至75%左右(稱重法),放入25℃恒溫培養(yǎng)箱黑暗預(yù)培養(yǎng)。根據(jù)培養(yǎng)前重量每4天補(bǔ)充一次水分(日蒸發(fā)量約為0.85 mL),共補(bǔ)充3次。第12天補(bǔ)充水分的同時(shí)施入(NH4)2SO4溶液(300 mgN·kg-1干土,足量氮使微生物對根系分泌物敏感,且使土壤產(chǎn)生的N2O同位素豐度在儀器測結(jié)果最準(zhǔn)確的范圍內(nèi))。在第14天用2.5 mL一次性注射器依次吸取草酸、絲氨酸、葡萄糖溶液和蒸餾水,插入土壤緩慢推動(dòng)注射器使溶液緩慢送出,作為正式培養(yǎng)第1天,之后每隔24 h添加一次。盡管絲氨酸能夠提供N素,但是培養(yǎng)期間添加絲氨酸N素的總量不足土壤總氮的3%,對N2O排放量的影響較小。

    1.3 主要藥品及規(guī)格

    草酸(H2C2O4·2H2O):分析純,北京化工廠;D-無水葡萄糖(C6H12O6):>99.9%,MYM Biological Technology;絲氨酸(C3H7NO3):99%,北京百靈威科技有限公司;硫酸銨((NH 4)2SO4):分析純,天津市福晨化學(xué)試劑廠。

    1.4 氣體樣品采集

    在正式培養(yǎng)的第3天,添加根系分泌物溶液后的第2、4、8、12、16、20和24小時(shí)取樣,共采集7次。采樣前將培養(yǎng)瓶蓋上橡膠塞和鋁蓋,并用密封器壓緊密封,培養(yǎng)2 h。其他時(shí)間保持培養(yǎng)瓶敞口。用50 mL針筒抽取30 mL氦氣注入培養(yǎng)瓶中,來回抽動(dòng)針筒重復(fù)3次,將氣體混勻,再將針筒抽到30 mL刻度處,待針筒不回彈后拔出并迅速插入已抽成真空的20 mL血清瓶,做好標(biāo)記。取樣的同時(shí)抽取實(shí)驗(yàn)室內(nèi)空氣作為N2O的環(huán)境背景值(3次重復(fù))。

    1.5 樣品測定與方法:

    1.5.1 N2O同位素特征值測定與校正 穩(wěn)定同位素質(zhì)譜儀(Delta V Plus‐Precon, Thermo Fisher Scientific, Bremen, Germany)用來測定N2O同位素特征值:15Nα15Nbulk18O計(jì)算公式如下[32]:

    18O(‰)(2)

    式中,=15N/14N,=18O/16O。sa表示樣品,st表示標(biāo)準(zhǔn)物。15N和18O的標(biāo)準(zhǔn)分別是大氣中的N2和維也納標(biāo)準(zhǔn)平均海洋水(Vienna Standard Mean Ocean Water),校準(zhǔn)用標(biāo)準(zhǔn)參考?xì)怏w(美國液化空氣特種氣體有限責(zé)任公司(Air Liquide America, Specialty Gases LLC))。

    由于測定的N2O為空氣中和土壤產(chǎn)生的N2O的混合氣,受空氣同位素值的影響,各檢測結(jié)果(15Nα15Nbulk18O)需按照如下公式校正[33]:

    式中,soil-emitted代表土壤排放N2O的同位素值;mix、mix分別表示土壤與空氣混合氣體的同位素值和N2O濃度;air、air分別表示空氣背景的同位素值和N2O濃度(0.63 mg·m-3。mix-air為土壤排放的N2O濃度。

    由于土壤N2O排放具有時(shí)間變異性,N2O的日平均同位素特征值用通量的加權(quán)平均值來表示[34]:

    式中,wa表示同位素的日通量加權(quán)平均值。e表示第-1取樣到次取樣時(shí)間內(nèi)土壤排放N2O的量(時(shí)刻的排放速率×?xí)r間間隔(2 h 或4 h))。δ表示次取樣檢測的同位素特征值。

    1.5.2 N2O來源確定的理論基礎(chǔ) N2O為直線型分子:Nβ-Nα=O。根據(jù)兩個(gè)氮原子與氧原子的相對位置不同,則中間的氮原子命名為α氮原子,末端的氮原子為β氮原子[32]。自然豐度條件下N2O分子中15N原子出現(xiàn)在α位和β位的差異,稱為同位素位嗜值(site preference),即SP值,計(jì)算公式如下[35-36]:

    =15Nα-15Nβ(6)

    在不同微生物群落或酶參與下的不同反應(yīng)過程產(chǎn)生的N2O的值不同,這是值用于土壤N2O溯源研究的重要理論基礎(chǔ)。純細(xì)菌培養(yǎng)下硝化作用和反硝化作用的值是已知的,且其值不受底物氮氧同位素的影響。但自然狀況下這兩種過程往往同時(shí)存在[32],為測定氣體樣品中N2O的同位素值,通常利用同位素二源混合模型計(jì)算硝化作用和反硝化作用對N2O的貢獻(xiàn)[37]。公式如下:

    N=1-D(8)

    式中,N和D分別表示N2O來自土壤中硝化和反硝化作用的貢獻(xiàn)比例N表示純細(xì)菌培養(yǎng)條件下硝化作用的值,為33‰[38-39];D表示純細(xì)菌培養(yǎng)條件下反硝化作用的值,為0‰[40-41];sa表示測得樣品的SP值,樣品值大于33‰或小于0‰均按33‰或0‰計(jì)算。

    1.5.3 N2O排放速率的計(jì)算 采集的氣體用穩(wěn)定同位素質(zhì)譜儀(Delta V Plus-Preco,Thermo Fisher Scientific,Bremen,Germany)結(jié)合痕量氣體濃縮系統(tǒng)測定N2O氣體峰面積,根據(jù)N2O在樣品中峰面積與環(huán)境背景中峰面積的比值,計(jì)算樣品N2O濃度,根據(jù)培養(yǎng)前后瓶內(nèi)N2O的濃度變化計(jì)算單位時(shí)間內(nèi)土壤N2O的排放速率(培養(yǎng)前以及環(huán)境背景氣體N2O的體積分?jǐn)?shù)為320 nmol·mol-1)。公式如下[42]:

    式中,表示N2O排放速率(μg·kg·h-1),22.4為標(biāo)準(zhǔn)大氣壓下的摩爾體積(L·mol-1)。為標(biāo)準(zhǔn)狀態(tài)下N2O氣體的摩爾質(zhì)量(44 g·mol-1)。1000為單位換算。表示培養(yǎng)瓶的體積,0.09 L。sa、air分別表示培養(yǎng)后瓶內(nèi)氣體和空氣中N2O濃度(nmol·mol-1)。Δ為培養(yǎng)時(shí)間,2 h。為培養(yǎng)溫度,25℃。

    1.5.4 N2O通量和同位素值的校正系數(shù) 各采樣時(shí)刻N(yùn)2O排放速率和日平均排放速率的校正,各采樣時(shí)刻測得的同位素值和同位素的日均值之間的校正公式如下[43-44]:

    式中,F(=2、4、8、12、16、20、24)表示各時(shí)間的平均通量。C為校正系數(shù),ave為N2O日平均排放速率、15Nbulk、18O和SP的日平均值,其中15Nbulk、18O和SP的日平均值為各時(shí)刻排放速率的加權(quán)平均值(根據(jù)公式(5)計(jì)算)。F為實(shí)測的第次實(shí)測通量,或第次取樣檢測的同位素值。校正系數(shù)越接近1,則說明該時(shí)刻測量的N2O排放速率和同位素值越接近日平均值,測量結(jié)果越具有代表性,即可采用該時(shí)刻作為最佳的取樣時(shí)間[45]。

    1.6 數(shù)據(jù)處理與分析

    采用Microsoft Excel·2016對試驗(yàn)數(shù)據(jù)進(jìn)行處理和做圖。用SPSS·25.0(SPSS Inc., Chicago, IL, USA)進(jìn)行統(tǒng)計(jì)分析。各處理間N2O累計(jì)排放量,18O,15Nbulk和SP用單因素方差分析(One-Way ANOVA),LSD法進(jìn)行多重比較。雙因素方差分析(Two-Way ANOVA)比較根系分泌物組分和濃度及其交互作用對N2O排放量和同位素特征值的影響。用一元線性回歸分析最佳采樣時(shí)間的檢測值和日均值間的關(guān)系。

    2 結(jié)果

    2.1 N2O排放速率的日變化

    如圖1所示,在取樣期間所有處理的N2O排放速率均呈上升趨勢,添加根系分泌物處理組均高于對照。在培養(yǎng)時(shí)間內(nèi)添加低濃度、高濃度草酸處理土壤N2O排放速率迅速升高,分別從13.21 μg·kg-1·h-1增加到226.28 μg·kg-1·h-1和從17.84 μg·kg-1·h-1增加到97.21 μg·kg-1·h-1,添加草酸后土壤N2O排放速率,低濃度處理組高于高濃度處理組;葡萄糖處理中低濃度組低于高濃度組,低濃度絲氨酸、高濃度絲氨酸、低濃度葡萄糖和高濃度葡萄糖處理的土壤N2O排放速率分別從70.72 μg·kg-1·h-1增加到125.79 μg·kg-1·h-1、50.82 μg·kg-1·h-1增加到171.71 μg·kg-1·h-1、20.34 μg·kg-1·h-1增加126.69 μg·kg-1·h-1和20.46 μg·kg-1·h-1增加到185.21 μg·kg-1·h-1。

    雙因素方差分析顯示N2O累計(jì)排放量受根系分泌物和濃度共同作用的影響(<0.05)(表1)。N2O累計(jì)排放量低濃度處理組:草酸處理>絲氨酸處理>葡萄糖處理,而在高濃度處理組:葡萄糖處理>絲氨酸處理>草酸處理。單因素方差分析顯示,低濃度草酸、高濃度絲氨酸和高濃度葡萄糖處理土壤N2O累計(jì)排放速率顯著高于對照(<0.05)。增加根系分泌物濃度時(shí),僅高濃度葡萄糖處理N2O累計(jì)排放量顯著高于低濃度葡萄糖處理(<0.05)。

    誤差線表示標(biāo)準(zhǔn)誤差(n=3)。下同

    表1 N2O累計(jì)排放量和同位素日加權(quán)平均值

    CL:150 μg C·d-1草酸處理;CH:300 μg C·d-1草酸處理;AL:150 μg C·d-1絲氨酸處理;AH:300 μg·C·d-1絲氨酸處理;PL:150 μg C·d-1葡萄糖處理;PH:300 μg·C·d-1葡萄糖處理;Control:蒸餾水處理(對照)。下同

    CL: 150 μg C·d-1Oxalic acid treatment; CH: 300 μg·C·d-1Oxalic acid treatment; AL:150 μg C·d-1serine treatment; AH:300 μg C·d-1serine Treatment; PL:150 μg C·d-1glucose treatment; PH: 300 μg C·d-1glucose treatment; Control: Distilled water treatment. The same as below

    18O,15Nbulk和SP的值為各處理不同采樣時(shí)間同位素值的加權(quán)平均值,由公式(5)計(jì)算。下同

    The weighted average of18O,15Nbulkand SP were calculated by equation (5). The same as below

    同列中小寫字母表示單因素方差分析(One-Way ANOVA)顯著性<0.05。*表示雙因素方差分析(Two-Way ANOVA)顯著性<0.05

    The values followed by different little letters in the same column indicate the significant differences between treatments at<0.05 level using One-Way ANOVA. * indicate the significant differences between treatments at<0.05 level using Two-Way ANOVA

    2.2 N2O同位素值的日變化

    不同處理土壤N2O同位素特征值不同(圖2)。添加根系分必物處理的18O變化基本一致,穩(wěn)定在23.43‰—25.90‰;而對照則在20.69‰—22.13‰。3種不同根系分泌物處理N2O的15Nbulk值不同,葡萄糖處理組15Nbulk值最高,變化范圍分別是:-15.22‰— 13.55 ‰(低濃度)和-14.82‰—-10.70‰(高濃度)。低濃度絲氨酸處理15Nbulk值的范圍是-26.56‰— 22.07‰,低于高濃度絲氨酸處理(-22.21‰— -16.66‰)。兩濃度草酸處理15Nbulk值差異較小,低、高濃度的變化范圍分別是-21.89‰—-15.27‰和-23.97‰—-16.75‰。3種根系分泌物濃度越高,SP值越低,低濃度草酸處理SP值從16.27‰增加至17.65‰,高濃度草酸處理SP值從16.69‰降至12.67‰;兩濃度的絲氨酸處理SP值逐漸升高,低濃度絲氨酸和高濃度絲氨酸處理變化范圍是:11.64‰—18.50‰和11.45‰—14.83‰;兩濃度葡萄糖處理SP值先升高后降低,極值分別出現(xiàn)在添加葡萄糖后的第8小時(shí)(18.79‰)和第12小時(shí)(17.26‰)。

    雙因素方差分析顯示,N2O的15Nbulk受添加根系分泌物成分的影響(<0.05),其中葡萄糖處理15Nbulk值最高((-13.86±1.11)‰)。SP值受添加根系分泌物濃度的影響(<0.05),濃度越高,SP值越低。單因素方差分析顯示添加根系分泌物處理N2O的18O值顯著高于對照(<0.05),但不同組分處理間無差異;添加葡萄糖和絲氨酸的不同濃度間N2O的SP值均無明顯差異;添加草酸后土壤N2O的SP值,低濃度處理組顯著高于高濃度處理組(<0.05)。

    圖2 添加不同根系分泌物后N2O同位素特征值(δ18O, δ15Nbulk, SP)隨時(shí)間的變化

    2.3 不同處理反硝化作用對N2O排放的貢獻(xiàn)

    通過二元混合模型(公式(6),(7))計(jì)算各處理反硝化作用和硝化作用對土壤N2O排放的相對貢獻(xiàn)率如圖3所示??瞻讓φ仗幚?、低濃度草酸、高濃度草酸、低濃度絲氨酸、高濃度絲氨酸、低濃度葡萄糖、高濃度葡萄糖處理反硝化作用的貢獻(xiàn)分別是:60%、48%、60%、52%、57%、51%、54%。添加草酸、絲氨酸、葡萄糖的濃度越高,反硝化作用越強(qiáng)。

    2.4 各時(shí)刻N(yùn)2O排放速率和同位素特征值的校正系數(shù)

    N2O排放速率和同位素特征值隨時(shí)間有一定的變化,添加根系分泌物后各時(shí)刻的檢測值和日平均值的校正系數(shù)如圖4所示。不同指標(biāo)的校正系數(shù)的波動(dòng)范圍不同,N2O排放速率的校正系數(shù)波動(dòng)較大(0.49—8.46),同位素值的校正系數(shù)波動(dòng)較?。?.79—1.87),其中各指標(biāo)前8 h校正系數(shù)的波動(dòng)較大,后逐漸減小。N2O排放速率的校正系數(shù)為1時(shí)的點(diǎn)在添加根系分泌物后的第16 h;N2O的18O值24 h內(nèi)比較穩(wěn)定,各處理校正系數(shù)穩(wěn)定在1附近;N2O的15Nbulk值的校正系數(shù)呈先增高后降低的趨勢,校正系數(shù)為1的點(diǎn)出現(xiàn)在添加根系分泌物后的第2—8小時(shí)和第16—20小時(shí);SP值的校正系數(shù)為1的點(diǎn)出現(xiàn)在添加根系分泌物后的第16—20小時(shí)。

    圖3 不同處理反硝化作用-硝化作用對N2O排放的貢獻(xiàn)

    圖4 添加根系分泌物后各時(shí)刻N(yùn)2O排放速率,穩(wěn)定同位素值(δ18O, δ15Nbulk,SP)的校正系數(shù)

    3 討論

    3.1 根系分泌物組分和濃度對N2O排放量及其同位素特征值的影響

    本研究發(fā)現(xiàn)添加根系分泌物處理土壤N2O排放量高于對照,說明根系分泌物的添加會(huì)促進(jìn)土壤N2O的排放,這和以往研究結(jié)果一致[46-47]。植物根系分泌物使植物通過光合作用固定的碳外流[48],并以有機(jī)碳的形式釋放到土壤中。微生物從根系分泌物中獲取能量,使土壤酶數(shù)量增大,活性增強(qiáng),導(dǎo)致N2O排放量的增加[49]。植物根系分泌物也可作為碳源改變微生物細(xì)胞內(nèi)N2O產(chǎn)生過程的電子傳遞途徑[50],并促進(jìn)N2O的產(chǎn)生。微生物對不同根系分泌物的響應(yīng)機(jī)制不同,從而導(dǎo)致N2O產(chǎn)生路徑也會(huì)有差異,進(jìn)而造成N2O排放量的差異。

    本研究發(fā)現(xiàn)添加3種根系分泌物處理間土壤N2O的18O值無明顯差異,且均顯著高于對照(<0.05)。還原為N2O的反硝化過程(→→NO→ N2O)可看作O逐漸去除的過程,16O比18O更容易斷裂使得產(chǎn)物N2O中富集18O從而產(chǎn)生同位素分餾[38];同時(shí),N2O還原(N2O→N2)過程也會(huì)導(dǎo)致剩余N2O的18O值升高[51-52]。而在還原過程中產(chǎn)生的中間產(chǎn)物可能會(huì)與環(huán)境中的H2O發(fā)生O交換,從而減弱在還原過程中發(fā)生的分餾效應(yīng)[38]。以上反硝化過程中的每一個(gè)還原步驟均是在不同酶的催化下發(fā)生,且導(dǎo)致不同的O同位素分餾[53]。因此根系分泌物處理組土壤N2O的18O值高于對照組可能是由于根系分泌物的添加促進(jìn)了反硝化作用的進(jìn)行,使反硝化還原酶活性改變所致;且3種根系分泌物可能共同作用了某一種對O同位素分餾起關(guān)鍵作用的酶而導(dǎo)致處理間N2O的18O值無差異。另外,添加3種根系分泌物可能造成了不同程度的N2O向N2的還原。

    各根系分泌物處理N2O的15Nbulk值不同,15Nbulk值受根系分泌物組分的影響。在肥料充足的土壤條件下,N2O的15Nbulk值的升高和降低是由多種原因造成的,同位素的分餾可能是重要的原因。硝化作用的產(chǎn)物可作為反硝化作用的底物,根據(jù)同位素分餾原理,微生物優(yōu)先利用輕同位素而導(dǎo)致產(chǎn)物發(fā)生重同位素貧化[54]使由反硝化過程產(chǎn)生的N2O的15Nbulk值降低;其次,N2O還原可能對15Nbulk值產(chǎn)生重要影響,N2O還原過程中14N-O鍵比15N-O鍵更容易斷裂,導(dǎo)致剩余N2O的15Nbulk值升高;所以添加根系分泌物通過促進(jìn)反硝化過程和N2O還原可能會(huì)導(dǎo)致N2O的15Nbulk值升高[30]。又因?yàn)橥寥乐械牟煌⑸镞^程導(dǎo)致的N同位素分餾不同[55],這使得不同根系分泌物處理的土壤N2O的15Nbulk值有差異。由于葡萄糖會(huì)促進(jìn)N2O還原過程,該理論已在本課題組其他研究中得到驗(yàn)證[56]。此外,Speir等[57]通過15N標(biāo)記肥料也發(fā)現(xiàn),葡萄糖為易有效碳,而易有效碳的添加會(huì)促進(jìn)反硝化過程進(jìn)行完全,促進(jìn)N2O向N2的還原過程,這可能是葡萄糖處理組N2O的15Nbulk值顯著高于對照的原因(<0.05)。

    SP值表征的是直線型分子N2O(N-N-O)的分子內(nèi)15N的位置嗜好性差異,該值和微生物過程有關(guān),是一個(gè)被國際認(rèn)可的N2O溯源研究工具[35]。本研究發(fā)現(xiàn)同一根系分泌物組分作用下,SP值隨根系分泌物濃度的升高而降低。由于純細(xì)菌培養(yǎng)下硝化作用下產(chǎn)生的N2O的SP值為33‰[38-39],反硝化作用下產(chǎn)生的N2O的SP值為0‰[40-41]。本研究所有處理SP值介于兩者之間則說明產(chǎn)生的N2O來源于硝化作用和反硝化作用的共同作用。根據(jù)硝化和反硝化的二元混合模型可知SP值的升高和降低分別代表土壤微生物過程趨于硝化作用方向和反硝化作用進(jìn)行,SP值越低,則反硝化作用越強(qiáng)。根系分泌物的添加為土壤提供碳源,從而促進(jìn)了異養(yǎng)微生物作用下的反硝化過程。在適宜范圍內(nèi),反硝化作用隨著土壤碳源濃度的升高而增強(qiáng),則產(chǎn)生的N2O的SP值也會(huì)因此而降低。

    3.2 基于N2O排放速率及同位素特征值確定最佳取樣時(shí)間

    在采取手動(dòng)N2O氣體觀測中,通常采用1次的檢測結(jié)果而作為1日甚至多日的平均情況,所以采樣時(shí)間的選擇尤為重要。根據(jù)研究需要,基于不同指標(biāo)的取樣時(shí)間不同,N2O排放速率、N2O的15Nbulk及SP值的最佳取樣時(shí)間分別為添加根系分泌物后的第12—16 h、4—8 h或16—20 h和16 h。由于N2O的18O值隨著培養(yǎng)時(shí)間無明顯變化,因此,取樣時(shí)間對N2O的18O值的影響不大。綜合N2O的排放速率和各穩(wěn)定同位素值的取樣時(shí)間,選擇添加根系分泌物后的第16 h為最佳取樣時(shí)間,從而可以對各個(gè)指標(biāo)進(jìn)行綜合分析,更有利于量化不同微生物過程對土壤N2O排放的貢獻(xiàn)率。對各處理最佳取樣時(shí)間下得到的各個(gè)指標(biāo)的結(jié)果分別與日平均值做相關(guān)分析顯示(圖5)。N2O排放速率:2=0.93(<0.001),N2O的18O值:2=0.71(<0.001),N2O的15Nbulk值:2=0.91(<0.05),SP值:2=0.91(<0.001)。因此,確定不同處理組的最佳取樣時(shí)間是添加根系分泌物后的第16小時(shí)。此外,該時(shí)刻同位素值相對穩(wěn)定,說明土壤內(nèi)微生物過程穩(wěn)定,利于分析N2O產(chǎn)生的微生物過程,因此結(jié)果最具代表性。

    圖5 最佳取樣時(shí)間N2O排放速率、穩(wěn)定同位素值(δ18O, δ15Nbulk, SP)和日平均值的比較

    3.3 研究展望

    鑒于穩(wěn)定同位素技術(shù)在區(qū)分N2O產(chǎn)生的微生物過程中存在優(yōu)勢,未來將繼續(xù)利用同位素技術(shù)圍繞不同根系分泌物組分及其組合作用下對土壤N2O的貢獻(xiàn)規(guī)律及其產(chǎn)生的多種微生物途徑進(jìn)行研究。由于N2O還原對同位素值的影響較大,且促進(jìn)N2O還原為N2是減少農(nóng)業(yè)土壤N2O排放的努力方向,所以未來研究將N2O還原納入考慮,同時(shí)引入國際穩(wěn)定同位素先進(jìn)技術(shù)和方法,如15Nbulk-18O模型和15Nbulk-SP模型,通過對反應(yīng)底物濃度、底物同位素和前體同位素圖譜的綜合分析,進(jìn)一步闡明植物根系分泌物對土壤N2O產(chǎn)生的作用機(jī)制。

    4 結(jié)論

    通過外源添加根系分泌物3種主要組分,研究其對N2O排放高峰期的排放速率和穩(wěn)定同位素特征值的影響。在施加NH+ 4的土壤中添加草酸、絲氨酸、葡萄糖促進(jìn)土壤N2O的排放,高濃度組處理葡萄糖的促進(jìn)效果最強(qiáng),而低濃度組處理草酸的促進(jìn)效果最強(qiáng);根系分泌物作用下土壤排放N2O的18O的值升高;根系分泌物不同組分作用下N2O的15Nbulk值不同,葡萄糖的添加使15Nbulk值顯著升高;各處理N2O的SP值的范圍為13.13‰—15.03‰,根系分泌物濃度越高SP值越低;反硝化作用對N2O的貢獻(xiàn)越強(qiáng)。通過對所測指標(biāo)(N2O排放速率、15Nbulk、18O和SP值)的校正系數(shù)和回歸分析顯示添加根系分泌物后的最佳取樣時(shí)間為添加根系分泌物后的第16小時(shí)。

    [1] Chapuis Lardy L, Wrage N, Metay A, Chotte J L, Bernoux M. Soils, a sink for N2O? A review.2010, 13(1): 1-17.

    [2] Conrad R. Soil microorganisms as controllers of atmospheric trace gases (H2, CO, CH4, OCS, N2O, and NO).1996, 60(4): 609-640.

    [3] SEITZINGER S P, KROEZE C, STYLES R V. Global distribution of N2O emissions from aquatic systems: natural emissions and anthropogenic effects.2000, 2(3/4): 267-279.

    [4] NEVISON C, HOLLAND E. A reexamination of the impact of anthropogenically fixed nitrogen on atmospheric N2O and the stratospheric O3layer.. 1997, 102(D21): 25519-25536.

    [5] 韓麗梅, 鞠會(huì)艷, 楊振明. 兩種基因型大豆根分泌物對大豆根腐病菌的化感作用. 應(yīng)用生態(tài)學(xué)報(bào), 2005(1): 137-141.

    HAN L M, JU H Y, YANG Z M. Allelopathy of root exudates from two genotypes soybeans on root rot pathogenic fungi.2005(1): 137-141. (in Chinese)

    [6] STAVI I, LAL R. Agriculture and greenhouse gases, a common tragedy. A review.2013, 33(2): 275-289.

    [7] 胡春勝, 張玉銘, 秦樹平, 王玉英, 李欣曉, 董文旭. 華北平原農(nóng)田生態(tài)系統(tǒng)氮素過程及其環(huán)境效應(yīng)研究. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(10): 1501-1514.

    HU C S, ZHANG Y M, QIN S P, WANG Y Y, LI X X, DONG W X. Nitrogen processes and related environmental effects on agroecosystem

    in the North China Plain.2018, 26(10): 1501-1514. (in Chinese)

    [8] MAN-YOUNG J, REINHARD W, DEULLAE M, ANETTE G, SOO JE P, JONG GEOL, SO JEONG K, SUNG KEUN RHEE. Isotopic signatures of N2O produced by ammonia-oxidizing archaea from soils.2013, 8(5): 1115-1125.

    [9] Bernot M J, Dodds W K, Gardner W S, mccarthy M J, Sobolev D, Tank L. Comparing denitrification estimates for a texas estuary by using acetylene inhibition and membrane inlet mass spectrometry., 2003, 69(10): 5950-5956.

    [10] Sutka R L, Ostrom N E, Ostrom P H, Breznak J A, Gandhi H, Pitt A J, Li F. Distinguishing nitrous oxide production from nitrification and denitrification on the basis of isotopomer abundances.2006 72(1): 638-644.

    [11] BAGGS E, MAIR L, MAHMOOD S. A review of stable isotope techniques for N2O source partitioning in soils: recent progress, remaining challenges and future considerations.2007, 22(11): 1664-1672.

    [12] 丁軍軍, 張薇, 李玉中, 林偉, 徐春英, 李巧珍.不同灌溉量對華北平原菜地N2O排放及其來源的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(7): 2269-2276.

    DING J J, ZHANG W, LI Y Z, LIN W, XU C Y, LI Q Z. Effects of soil water condition on N2O emission and its sources in vegetable farmland of North China Plain.2017, 28(7): 2269-2276. (in Chinese)

    [13] 鄭欠, 丁軍軍, 李玉中, 林偉, 徐春英, 李巧珍, 毛麗麗. 土壤含水量對硝化和反硝化過程N(yùn)2O排放及同位素特征值的影響. 中國農(nóng)業(yè)科學(xué), 2017, 50(24): 4747-4758.

    ZHENG Q, DING J J, LI Y Z, LIN W, XU C Y, LI Q Z, MAO L L. The effects of soil water content on N2O emissions and isotopic signature of nitrification and denitrification.2017, 50(24): 4747-4758. (in Chinese)

    [14] ?UHEL J, ?IMEK M, LAΜGHLIN R J, BRU D, CHENEBY D. WATSON CJ, PHILIPPOT L. Insights into the effect of soil pH on N2O and N2emissions and denitrifier community size and activity.2010, 76(6): 1870-1878.

    [15] HANG-WEI H, DELI C, JI-ZHENG H. Microbial regulation of terrestrial nitrous oxide formation: understanding the biological pathways for prediction of emission rates.2015, 39(5): 729-749.

    [16] SHUKRA RAJ P, OHKYUNG C, SAMIR KUMAR K, KARTIK C, SUNGPYO K, JAW WOO L. Effects of temperature on nitrous oxide (N2O) emission from intensive aquaculture system.2015, 518-519: 16-23.

    [17] LAI T V, FARQUHARSON R, DENTON M D. High soil temperatures alter the rates of nitrification, denitrification and associated N2O emissions.2019,19(5): 2176-2189.

    [18] BOUWMAN A F, BEUSEN A H W, GRIFFIOEN J, VAN GROENIGEN J W, HEFTING M M, OENEMA O, VAN PUIJENBROEK P J T M, SEITZINGE R S, SLOMP C P, STEHFEST E. Global trends and uncertainties in terrestrial denitrification and N2O emissions.2013, 368(1621): 91-97.

    [19] JONES D L, NGUYEN C, FINLAY R D. Carbon flow in the rhizosphere: carbon trading at the soil–root interface.2009, 321(1/2): 5-33.

    [20] BAIS H P, WEIR T L, PERRY L G, GILROY S, VIVANCO J M. The role of root exudates in rhizosphere interactions with plants and other organisms.,2006, 57(1): 233-266.

    [21] HAICHAR F E Z, SANTAELLA C, HEULIN T, ACHOUAK W. Root exudates mediated interactions belowground.2014, 77(7): 69-80.

    [22] BATEMAN E J, BAGGS E M. Contributions of nitrification and denitrification to N2O emissions from soils at different water-filled pore space.2005, 41(6): 379-388.

    [23] WU H, WANG X, HE X, ZHANG S, LIANG R, SHEN J. Effects of root exudates on denitrifier gene abundance, community structure and activity in a micro-polluted constructed wetland.2017, 598: 697-703.

    [24] 馬舒坦, 顏曉元. 甲酸鹽和葡萄糖對兩種土壤N2O排放的刺激作用. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2019, 38(1): 241-248.

    MA S T, YAN X Y. Effect of formate and glucose organic carbon on N2O emission from two soils.2019, 38(1): 241-248. (in Chinese)

    [25] LAPALANGARICA FUENTES A, MANRUBIA M, GILES M E, MITCHELL S, DANIELL T J. Effect of model root exudate on denitrifier community dynamics and activity at different water-filled pore space levels in a fertilised soil.2018, 120: 70-79.

    [26] 徐衛(wèi)紅, 劉吉振, 黃河, 熊治廷. 高鋅脅迫下不同大白菜品種生長、Zn吸收及根系分泌物的研究. 中國農(nóng)學(xué)通報(bào), 2006, 22(8): 458.

    XU W H, LIU J Z, HUANG H, XIONG Z Y. Study of Zn stress on plant growth, Zn uptake and root exudates in different cultivars of Chinese Cabbage.2006, 22(8): 458. (in Chinese)

    [27] YUAN Y S, ZHAO W Q, ZHANG Z L, XIAO J, LI D D, LIU Q, YIN H J. Impacts of oxalic acid and glucose additions on N transformation in microcosms via artificial roots.2018, 121: 16-23.

    [28] BOL R, R?CKMANN T, BLACKWELL M, YAMULKI S. Influence of flooding on15N,18O,115N and215N signatures of N2O released from estuarine soils-a laboratory experiment using tidal flooding chambers.2004, 18(14): 1561-1568.

    [29] HU G Q, HE H B, ZHANG W, ZHAO J S, CUI J H, LI B, ZHANG X D. The transformation and renewal of soil amino acids induced by the availability of extraneous C and N.2016, 96: 86-96.

    [30] HENRY S, TEXIER S, HALLET S, BRU D, DAMBREVILLE C, CHèNEBY D, BIZOUARD F, GERMON J C, PHILIPPOT L. Disentangling the rhizosphere effect on nitrate reducers and denitrifiers: insight into the role of root exudates.2010, 10(11): 3082-3092.

    [31] OBURGER E, JONES D L. Sampling root exudates – mission impossible?2018, 6:116-133.

    [32] TOYODA S, YOSHIDA N. Determination of nitrogen isotopomers of nitrous oxide on a modified isotope ratio mass spectrometer., 1999, 71(20): 4711-4718.

    [33] OPDYKE M R, OSTROM N E, OSTROM P H. Evidence for the predominance of denitrification as a source of N2O in temperate agricultural soils based on isotopologue measurements.2009, 23(4). DOI:10.1029/2009GB003523.

    [34] 房福力, 李玉中. 基于同位素特征的華北平原菜地N2O排放監(jiān)測中取樣時(shí)間的探討. 植物營養(yǎng)與肥料學(xué)報(bào), 2016, 22(4): 978-987.

    FANG F L, LI Y Z. Preliminary research on N2O sampling time based on isotopomer signature measurement of vegetable fields in the North China Plain.2016, 22(4): 978-987. (in Chinese)

    [35] Toyoda S, Mutobe H, Yamagishi H, Yoshida N, Tanji Y. Fractionation of N2O isotopomers during production by denitrifier.2005, 37(8): 1535-1545.

    [36] DECOCK C, SIX J. How reliable is the intramolecular distribution of15N in N2O to source partition N2O emitted from soil?2013, 65: 114-127.

    [37] ZHANG W, LI Y Z, XU C Y, LI Q Z, LIN W. Isotope signatures of N2O emitted from vegetable soil: Ammonia oxidation drives N2O production in NH+ 4-fertilized soil of North China., 2016, 6: 29257.

    [38] Toyoda S, Mutobe H, Yamagishi H, Yoshida N, Tanji Y. Fractionation of N2O isotopomers during production by denitrifier.2005, 37(8): 1535-1545.

    [39] Jung M Y, Well R, Min D, Giesemann A, Rhee S K. Isotopic signatures of N2O produced by ammonia-oxidizing archaea from soils.2013, 8(5): 1115-1125.

    [40] Yamazaki T, Hozuki T, Arai K, Toyoda S, Koba K, Fujiwara T. Isotopomeric characterization of nitrous oxide produced by reaction of enzymes extracted from nitrifying and denitrifying bacteria.2014, 11(10): 2679-2689.

    [41] BRAKER G, LEWICKASZCZEBAK D. Dual isotope and isotopomer signatures of nitrous oxide from fungal denitrification a pure culture study.2014, 28(17): 1893-1903.

    [42] 張仲新, 李玉娥, 華珞, 萬運(yùn)帆, 姜寧寧. 不同施肥量對設(shè)施菜地N2O排放通量的影響. 農(nóng)業(yè)工程學(xué)報(bào), 2010, 26(5): 269-275.

    ZHANG Z X, LI Y E, HUA L, WAN Y F, JIANG N N. Effects of different fertilizer levels on N2O flux from protected vegetable land.2010, 26(5): 269-275. (in Chinese)

    [43] 田慎重, 寧堂原, 遲淑筠, 王瑜, 王丙文, 韓惠芳, 李成慶, 李增嘉. 不同耕作措施的溫室氣體排放日變化及最佳觀測時(shí)間. 生態(tài)學(xué)報(bào), 2012, 32(3): 879-888.

    Tian s Z, ning t Y, chi s Y, wang y, wang b W, han h F, li c Q, li z J. Diurnal variations of the greenhouse gases emission and their optimal observation duration under different tillage systems.2012, 32(3): 879-888. (in Chinese)

    [44] 徐鈺, 劉兆輝, 石璟, 魏建林, 李國生, 王梅, 江麗華. 北方設(shè)施菜地土壤N2O排放通量日變化及最佳觀測時(shí)間確定. 中國農(nóng)業(yè)氣象, 2016, 37(5): 505-512.

    XU Y, LIU Z H, SHI J, WEI J L, LI G S, WANG M, JIANG L H. Diurnal variation characteristic of nitrous oxide from greenhouse vegetable soil during emission peak and its optimal observation duration.2016, 37(5): 505-512. (in Chinese)

    [45] 鄭循華, 王明星, 王躍思. 太湖地區(qū)農(nóng)田NO排放不連續(xù)測量最佳時(shí)間. 環(huán)境科學(xué), 2000, 21(1): 1-6.

    Zheng x H, wang m X, wang y S. Diurnal variation of NO emission from an ecosystem and the optimum observation time for intermittent flux measurement.2000, 21(1): 1-6. (in Chinese)

    [46] 王琳, 周曉麗, 馬銀麗, 巨曉棠, 吉艷芝, 張麗娟. 銨態(tài)氮源和碳源對土壤N2O、CO2釋放的影響. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào), 2016, 33(1): 23-28.

    WANG L, ZHOU X L, MA Y L, JU X T, JI Y Z , ZHANG L J. Effect of ammonium nitrogen source and carbon source on the CO2and N2O emissions of soil.2016, 33(1): 23-28. (in Chinese)

    [47] 李貴桐, 李保國, 黃元仿. 碳源與底物對不同層次土壤產(chǎn)生N2O能力的影響. 土壤與環(huán)境, 2002, 11(3): 227-231.

    LI G T, LI B G, HUANG Y F. Effects of carbon and substrate on the N2O productivity of different soil layers.2002, 11(3): 227-231. (in Chinese)

    [48] NGUYEN C. Rhizodeposition of organic C by plants: mechanisms and controls.2003, 23: 375-396.

    [49] HENRY S, TEXIER S, HALLET S, BRU D, DAMBREVILLE C, CHèNEBY D, BIZOUARD F, GERMON J C, PHILIPPOT L. Disentangling the rhizosphere effect on nitrate reducers and denitrifiers: Insight into the role of root exudates., 2010, 10(11): 3082-3092.

    [50] Mothapo N, Chen H, Cubeta M A, Grossman J M, Fuller F, Shi W. Phylogenetic, taxonomic and functional diversity of fungal denitrifiers and associated N2O production efficacy.2015, 83: 160-175.

    [51] Kool D M, Wrage N, Oenema O, Dolfing J, Van Groenigen J W. Oxygen exchange between (de)nitrification intermediates and H2O and its implications for source determination of NO- 3, and N2O: a review.2010, 21(22): 3569-3578.

    [52] CASCIOTTI K L, SIGMAN D M, HASTINGS M G, BOHLKE J K, HILKERT A. Measurement of the oxygen isotopic composition of nitrate in seawater and freshwater using the denitrifier method.2002, 74(19): 4905-4912.

    [53] ROHE L, WELL R, LEWICKA SZCZEBAK D. Use of oxygen isotopes to differentiate between nitrous oxide produced by fungi or bacteria during denitrification.2017, 31(16): 1297-1312.

    [54] PéREZ T, TRUMBORE S. E, TYLER S C, MATSON P A, ORTIZ MONASTERIO I, RAHN T, GRIFFITH D W T. Identifying the agricultural imprint on the global N2O budget using stable isotopes.2001, 106(9): 9869-9878.

    [55] Hyodo A, Malghani S, Zhou Y, Mushinski R M, West J B. Biochar amendment suppresses N2O emissions but has no impact on15N site preference in an anaerobic soil.2019, 33(2):165-175.

    [56] 鄭欠. 土壤水分和pH對N2O排放及同位素特征值影響的機(jī)理研究[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2018.

    ZHENG Q. The study on mechanism of N2O emissions and isotopic signature under different soil water content and pH[D]. Beijing: Chinese Academy of Agricultural Sciences, 2018. (in Chinese)

    [57] SPEIR T W, KETTLES H A, MORE R D. Aerobic emissions of N2O and N2from soil cores: Factors influencing production from13N-labelled NO3-and NH4+.1995, 27(10): 1299-1306.

    [58] Henderson S L, Dandie C E, Goyer C, Patten C, Zebarth B J, Burton D L, Trevors J T. Changes in denitrifier abundance, denitrification gene mRNA levels, nitrous oxide emissions, and denitrification in anoxic soil microcosms amended with glucose and plant residues.2010, 76: 2155-2164.

    [59] 賈俊仙, 李忠佩, 車玉萍. 添加葡萄糖對不同肥力紅壤性水稻土氮素轉(zhuǎn)化的影響. 中國農(nóng)業(yè)科學(xué), 2010, 43(8): 1617-1624.

    JIA J X, LI Z P, CHE Y P. Effects of glucose addition on N transformations in paddy soils with a gradient of organic C content in subtropical China., 2010, 43(8): 1617-1624. (in Chinese)

    [60] 崔敏, 冉煒, 沈其榮. 水溶性有機(jī)質(zhì)對土壤硝化作用過程的影響. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2006, 22(3): 45-50.

    CUI M, RAN W, SHEN Q R. Effects of dissolved organic matter on nitrification in soil., 2006, 22(3): 45-50. (in Chinese)

    [61] 王風(fēng), 陳思, 楊厚花, 沈仕洲, 張克強(qiáng), 王曉光. 葡萄糖添加對室溫和凍結(jié)過程土壤N2O排放特征影響. 生態(tài)科學(xué), 2017, 36(3): 31-35.

    WANG F, CHEN S, YANG H H, SHEN S Z, ZHANG K Q, WANG X G. Effect of glucose addition on N2O emission from three types of cultivated soils under ambient and freezing temperature.2017, 36(3): 31-35. (in Chinese)

    [62] 姜寧寧, 李玉娥, 華珞, 萬運(yùn)帆, 石生偉. 不同氮源及秸稈添加對菜地土壤N2O排放影響. 土壤通報(bào), 2012, 43(1): 219-223.

    JANG N N, LI Y E, HUA L, WAN Y F, SHI S W. Effect of different nitrogen sources and straw adding on N2O emission from vegetable soil.2012, 43(1): 219-223. (in Chinese)

    [63] 林偉, 張薇, 李玉中, 徐春英, 徐春英, 李巧珍, 鄭欠. 有機(jī)肥與無機(jī)肥配施對菜地土壤N2O排放及其來源的影響. 農(nóng)業(yè)工程學(xué)報(bào), 2016, 32(19): 148-153.

    LIN W, ZHANG W, LI Y Z, XU C Y, LI Q Z, ZHENG Q. Effects of combined application of manure and inorganic fertilizer on N2O emissions and sources in vegetable soils.2016, 32(19): 148-153. ( in Chinese)

    [64] 林偉, 丁軍軍, 李玉中, 徐春英, 李巧珍, 鄭欠, 莊姍. 有機(jī)肥和無機(jī)肥對菜地土壤N2O排放及其來源的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(5): 100-108. (in Chinese)

    LIN W, DING J J, LI Y Z, XU C Y, LI Q Z, ZHENG Q, ZHUANG S. Effects of organic and inorganic fertilizers on emission and sources of N2O in vegetable soils.2018, 29(5): 100-108. (in Chinese)

    [65] MCLAIN J E T, MARTENS D A. Nitrous oxide flux from soil amino acid mineralization., 2005, 37(2): 289-299.

    [66] CLARHOLM M, SKYLLBERG U, ROSLING A. Organic acid induced release of nutrients from metal-stabilized soil organic matter- The unbutton model.2015, 84: 168-176.

    [67] Keiluweit M, Boμgoure J J, Nico P S, Pett Ridge J, Kleber M. Mineral protection of soil carbon counteracted by root exudates.2015, 5(6): 588-595.

    [68] MALY S, KRáLOVEC J, HAMPEL D. Effects of long-term mineral fertilization on microbial biomass, microbial activity, and the presence of R and K -strategists in soil.2009, 45(7): 753-760.

    [69] LIU B, BAKKEN L R, FROSTEGARD A. Denitrification gene pools, transcription and kinetics of NO, N2O and N2production as affected by soil pH.2010, 72(3): 407-417.

    Effects of Different Root Exudates on Soil N2O Emissions and Isotopic Signature

    ZHUANG Shan1, LIN Wei1, DING JunJun1, ZHENG Qian1, KOU XinYue1, LI QiaoZhen1, LI YuZhong1,2

    (1Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences/Key Laboratory of Dryland Farming Agriculture, Ministry of Agriculture and Rural Affairs, Beijing 100081;2Environmental Stable Isotope Laboratory, Chinese Academy of Agricultural Sciences, Beijing 100081)

    【】The aim of this study was to investigate the effects of the main plant root exudates (organic acids, amino acids, sugars) on N2O emission and its microbial processes, so as to provide a support for selecting suitable plants to control soil N2O emissions. 【】 Three main components of root exudates, including oxalic acid, serine and glucose, were added in the soil, and two levels of concentration were set for each component: low concentration (150 μg C·d-1) and high concentration (300 μg C·d-1). There were totally 7 treatments with a control treatment treated with distilled water, and all treatments were incubated in 120 mL glass bottles in incubator. The gas samples were sampled 7 times within 24 hours after 2-hour incubation each time. After sampling, the N2O emission rate, daily cumulative emission and isotope signatures (15Nbulk,18O and SP (intramolecular15N site preference,15Nα-15Nβ)) were measured, then the optimal sampling time was determined according to their daily variation rules. 【】 The N2O emission rate of soils increased gradually after adding three components of root exudates, which were higher than the control treatment. The cumulative emission of N2O in the high concentration treatment was: Glucose treatments ((3.2±1.3) mg·kg-1·d-1)>Serine treatments ((2.6±0.5) mg·kg-1·d-1)>Oxalic acid treatments((1.4±0.2) mg·kg-1·d-1), low concentration treatment: Oxalic acid ((2.7±1.3) mg·kg-1·d-1)>Serine ((1.8±0.4) mg·kg-1·d-1)>Glucose ((1.6±0.8) mg·kg-1·d-1); the values of18O of N2O were not different among different root exudate treatments and were stable at 24.1‰-25.6‰, but significantly higher than the control treatment ((20.1±1.5) ‰);the15Nbulkvalue of N2O was related to the component of root exudates added, which was (-20.06±2.22) ‰ of oxalic acid treatment, (-22.33±1.10) ‰ of serine treatment, (-13.86±1.11) ‰ treated of glucose treatment, and (-23.14±3.72) ‰ of the control treatment. The SP value of N2O of each treatment ranged from 13.13‰ to 15.03‰, and the higher the root exudate concentration, the lower the SP value; after a comprehensive analysis of the correction coefficients of four indexes (N2O emission rate, the value of15Nbulk,18O and SP) at each sampling time and their daily mean values of 7 treatments, the correction coefficients of all treatments were closest to 1 at the 16th hour after the addition of root exudates.【】 In the soil environment with NH+ 4- 300 mg N·kg-1, the root exudates promoted N2O emission and the N2O emission rate increased gradually during the culture time (24 hours). The promotion effect of glucose in high concentration group was the strongest, while that of oxalic acid in low concentration group was the strongest. Compared with thecontrol treatment, the addition of root exudates significantly increased the18O value of N2O; the addition of glucose significantly increased the15Nbulkvalue of N2O. The higher the concentration of root exudates, the stronger the contribution of denitrification to N2O was detected.

    root exudates; nitrification; denitrification; N2O; isotopic signature

    10.3864/j.issn.0578-1752.2020.09.013

    2019-08-02;

    2019-11-06

    國家自然科學(xué)基金(41473004,41701308)、國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0201702)

    莊姍,Tel:18519190629;E-mail:489300802@qq.com。通信作者李玉中,Tel:18810871629;E-mail:liyuzhong@caas.cn

    (責(zé)任編輯 李云霞)

    欧美绝顶高潮抽搐喷水| 在线观看免费午夜福利视频| 国产精品亚洲av一区麻豆| 日本撒尿小便嘘嘘汇集6| 国内久久婷婷六月综合欲色啪| √禁漫天堂资源中文www| 女人爽到高潮嗷嗷叫在线视频| 日本免费一区二区三区高清不卡 | 自拍欧美九色日韩亚洲蝌蚪91| 国产成人精品久久二区二区免费| 亚洲全国av大片| 亚洲男人天堂网一区| 电影成人av| 一卡2卡三卡四卡精品乱码亚洲| 一本综合久久免费| 日本欧美视频一区| 精品国产乱码久久久久久男人| 一进一出抽搐动态| 99国产精品99久久久久| 亚洲五月色婷婷综合| av在线播放免费不卡| 日本免费a在线| 美女 人体艺术 gogo| 脱女人内裤的视频| 不卡一级毛片| av有码第一页| 操美女的视频在线观看| 又黄又粗又硬又大视频| 欧美在线黄色| 亚洲av熟女| 最新在线观看一区二区三区| 久久久久国产精品人妻aⅴ院| 婷婷精品国产亚洲av在线| 亚洲av美国av| 黄片大片在线免费观看| 精品国产乱子伦一区二区三区| 国产精品亚洲一级av第二区| 岛国视频午夜一区免费看| 亚洲精品久久国产高清桃花| 久久影院123| 又紧又爽又黄一区二区| 岛国视频午夜一区免费看| 午夜免费激情av| 亚洲第一电影网av| 999久久久精品免费观看国产| 国产私拍福利视频在线观看| 午夜免费激情av| 久久国产精品人妻蜜桃| 夜夜看夜夜爽夜夜摸| 国产精品久久久人人做人人爽| av超薄肉色丝袜交足视频| 欧美另类亚洲清纯唯美| 中文亚洲av片在线观看爽| 婷婷六月久久综合丁香| 桃红色精品国产亚洲av| 久久国产精品人妻蜜桃| 天天躁夜夜躁狠狠躁躁| 久久中文字幕一级| 老司机午夜福利在线观看视频| 久久精品亚洲精品国产色婷小说| 日韩欧美一区视频在线观看| 亚洲五月婷婷丁香| 成人欧美大片| 最好的美女福利视频网| 国产成人精品久久二区二区91| 99久久综合精品五月天人人| 亚洲国产精品合色在线| 欧美午夜高清在线| 欧美色视频一区免费| 国产成人av激情在线播放| 精品国产亚洲在线| 国产精品爽爽va在线观看网站 | 亚洲第一欧美日韩一区二区三区| 精品福利观看| 亚洲成av片中文字幕在线观看| 后天国语完整版免费观看| 50天的宝宝边吃奶边哭怎么回事| 乱人伦中国视频| 免费高清在线观看日韩| 国产激情久久老熟女| 九色亚洲精品在线播放| 国产一卡二卡三卡精品| 无人区码免费观看不卡| 久久久久国产精品人妻aⅴ院| av在线播放免费不卡| 久久精品91无色码中文字幕| 亚洲第一电影网av| 午夜福利高清视频| 叶爱在线成人免费视频播放| 午夜成年电影在线免费观看| 老汉色∧v一级毛片| 黄色视频不卡| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 精品人妻在线不人妻| 久久精品91无色码中文字幕| 国产1区2区3区精品| 自线自在国产av| 国产精品爽爽va在线观看网站 | 久久午夜综合久久蜜桃| 97碰自拍视频| 国产精品 国内视频| 亚洲天堂国产精品一区在线| 欧美日韩黄片免| 亚洲精品粉嫩美女一区| 极品人妻少妇av视频| 欧美乱色亚洲激情| 日本 av在线| 18美女黄网站色大片免费观看| 久久国产精品人妻蜜桃| 国产男靠女视频免费网站| 免费观看人在逋| 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| av在线天堂中文字幕| 久久精品国产99精品国产亚洲性色 | 精品一区二区三区视频在线观看免费| 嫩草影院精品99| 久久草成人影院| 亚洲av五月六月丁香网| 久久久水蜜桃国产精品网| 非洲黑人性xxxx精品又粗又长| 欧美乱码精品一区二区三区| 国产精品久久久人人做人人爽| 91国产中文字幕| 日日摸夜夜添夜夜添小说| 国产av在哪里看| 欧美日韩乱码在线| 91大片在线观看| 午夜视频精品福利| 欧美激情高清一区二区三区| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三区在线| 日日干狠狠操夜夜爽| 国产精品永久免费网站| 此物有八面人人有两片| 久久精品91无色码中文字幕| 亚洲av成人一区二区三| 午夜福利一区二区在线看| 国产av又大| 看黄色毛片网站| 亚洲成人精品中文字幕电影| 高清毛片免费观看视频网站| 国产熟女xx| 亚洲片人在线观看| 欧美日韩福利视频一区二区| 成年版毛片免费区| 国产精品九九99| 操美女的视频在线观看| 日本黄色视频三级网站网址| 午夜亚洲福利在线播放| 最近最新中文字幕大全电影3 | 日本欧美视频一区| av福利片在线| 久热爱精品视频在线9| 免费高清视频大片| 国产成人欧美在线观看| 老司机靠b影院| 免费高清在线观看日韩| 成人18禁高潮啪啪吃奶动态图| 国产成人欧美| 母亲3免费完整高清在线观看| 黄片播放在线免费| 国产在线观看jvid| 欧美日韩一级在线毛片| 老司机在亚洲福利影院| 国产一区二区三区综合在线观看| 精品日产1卡2卡| 亚洲人成网站在线播放欧美日韩| 日韩欧美免费精品| 欧美最黄视频在线播放免费| 久久久久久久午夜电影| 亚洲国产中文字幕在线视频| 亚洲 国产 在线| 亚洲av第一区精品v没综合| 久久影院123| www.自偷自拍.com| 日韩大码丰满熟妇| 亚洲自偷自拍图片 自拍| 宅男免费午夜| 91麻豆av在线| 天天躁夜夜躁狠狠躁躁| 99riav亚洲国产免费| 波多野结衣巨乳人妻| 丁香欧美五月| av福利片在线| 一边摸一边抽搐一进一出视频| 亚洲成av片中文字幕在线观看| 欧美黄色淫秽网站| 亚洲精品中文字幕在线视频| av超薄肉色丝袜交足视频| 日日爽夜夜爽网站| 巨乳人妻的诱惑在线观看| 国产av又大| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美一区二区三区在线观看| 欧美黄色淫秽网站| 纯流量卡能插随身wifi吗| 日本a在线网址| 久久人人爽av亚洲精品天堂| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| 亚洲专区国产一区二区| 大陆偷拍与自拍| 天堂影院成人在线观看| 在线视频色国产色| 国产色视频综合| 欧美日本视频| 88av欧美| 亚洲国产精品久久男人天堂| 在线视频色国产色| 午夜老司机福利片| 999久久久精品免费观看国产| 欧美激情高清一区二区三区| 精品国产乱子伦一区二区三区| 成人三级黄色视频| 国产精品美女特级片免费视频播放器 | videosex国产| 久久久久久久久免费视频了| 国产精品香港三级国产av潘金莲| 国产区一区二久久| 女人被躁到高潮嗷嗷叫费观| 一本久久中文字幕| 欧美不卡视频在线免费观看 | 久久中文看片网| 成人国语在线视频| 999久久久国产精品视频| 色精品久久人妻99蜜桃| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产亚洲在线| 电影成人av| 欧美日本视频| 亚洲中文av在线| www.熟女人妻精品国产| 国产三级在线视频| 国产欧美日韩综合在线一区二区| 国产精品免费视频内射| 大型av网站在线播放| 日日爽夜夜爽网站| 天堂动漫精品| 国产精品一区二区在线不卡| 在线观看免费视频网站a站| 在线观看免费午夜福利视频| 一边摸一边做爽爽视频免费| 亚洲aⅴ乱码一区二区在线播放 | 波多野结衣巨乳人妻| 亚洲视频免费观看视频| 中出人妻视频一区二区| 国产精品亚洲一级av第二区| 亚洲,欧美精品.| 美女午夜性视频免费| avwww免费| 日韩高清综合在线| 精品第一国产精品| 亚洲中文字幕一区二区三区有码在线看 | 国产免费男女视频| 日韩精品免费视频一区二区三区| 色精品久久人妻99蜜桃| 国产精品一区二区免费欧美| 老司机靠b影院| 亚洲人成77777在线视频| 老司机午夜福利在线观看视频| 九色亚洲精品在线播放| 欧美激情 高清一区二区三区| 一区在线观看完整版| 看片在线看免费视频| 免费看a级黄色片| 国产成人精品无人区| 午夜久久久久精精品| 99久久久亚洲精品蜜臀av| 嫩草影院精品99| 精品久久久久久成人av| 午夜两性在线视频| 俄罗斯特黄特色一大片| 色综合婷婷激情| 宅男免费午夜| 国产人伦9x9x在线观看| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 久久久久亚洲av毛片大全| 欧美成狂野欧美在线观看| 美女国产高潮福利片在线看| 亚洲av第一区精品v没综合| 色综合婷婷激情| 亚洲精品美女久久久久99蜜臀| 欧美另类亚洲清纯唯美| 日本欧美视频一区| 久久久久国产一级毛片高清牌| 韩国精品一区二区三区| 黄色片一级片一级黄色片| 精品一区二区三区视频在线观看免费| 亚洲人成电影观看| 欧美精品亚洲一区二区| 18禁美女被吸乳视频| 麻豆一二三区av精品| 欧美黄色片欧美黄色片| 亚洲情色 制服丝袜| 亚洲片人在线观看| 国产成人啪精品午夜网站| 成人三级黄色视频| 桃色一区二区三区在线观看| 黑人巨大精品欧美一区二区蜜桃| 女警被强在线播放| 欧美亚洲日本最大视频资源| 成年人黄色毛片网站| 国产乱人伦免费视频| 亚洲av第一区精品v没综合| 欧美日本视频| 国产精品精品国产色婷婷| 久久青草综合色| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 国产精品香港三级国产av潘金莲| 熟女少妇亚洲综合色aaa.| 在线视频色国产色| 97超级碰碰碰精品色视频在线观看| 午夜精品国产一区二区电影| 真人一进一出gif抽搐免费| 侵犯人妻中文字幕一二三四区| 久久人妻熟女aⅴ| 欧美日本视频| 欧美激情高清一区二区三区| 免费在线观看视频国产中文字幕亚洲| 国内毛片毛片毛片毛片毛片| 丝袜人妻中文字幕| 在线观看免费视频日本深夜| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久久人妻蜜臀av | 欧洲精品卡2卡3卡4卡5卡区| 黄色女人牲交| 制服人妻中文乱码| 一区二区日韩欧美中文字幕| 色尼玛亚洲综合影院| 国产又色又爽无遮挡免费看| 国产精品久久电影中文字幕| 少妇被粗大的猛进出69影院| 国产精品 欧美亚洲| av欧美777| 99久久精品国产亚洲精品| 亚洲国产欧美日韩在线播放| 亚洲人成电影免费在线| 老熟妇仑乱视频hdxx| 制服诱惑二区| 一进一出抽搐gif免费好疼| 制服诱惑二区| 午夜福利高清视频| 久久这里只有精品19| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 此物有八面人人有两片| 一进一出抽搐动态| 亚洲精品在线观看二区| 首页视频小说图片口味搜索| 成人三级黄色视频| 免费久久久久久久精品成人欧美视频| 免费少妇av软件| 国产精品亚洲av一区麻豆| 午夜福利18| 欧美色欧美亚洲另类二区 | 校园春色视频在线观看| av在线播放免费不卡| 国产亚洲欧美精品永久| 少妇被粗大的猛进出69影院| 久久国产精品人妻蜜桃| 国产一区二区三区综合在线观看| 99国产精品99久久久久| 久久中文字幕人妻熟女| 色婷婷久久久亚洲欧美| 老鸭窝网址在线观看| av免费在线观看网站| 日本 欧美在线| 国语自产精品视频在线第100页| 亚洲国产精品成人综合色| 在线十欧美十亚洲十日本专区| 亚洲精品中文字幕一二三四区| 757午夜福利合集在线观看| 国产成人精品无人区| 亚洲av五月六月丁香网| 日韩大码丰满熟妇| 色哟哟哟哟哟哟| 一区在线观看完整版| 久久久久久人人人人人| 国内精品久久久久久久电影| 成人三级黄色视频| 日韩欧美在线二视频| 国产免费av片在线观看野外av| 国产片内射在线| 午夜免费观看网址| 脱女人内裤的视频| 久久久精品欧美日韩精品| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 曰老女人黄片| a在线观看视频网站| 一区二区三区激情视频| 欧美色视频一区免费| ponron亚洲| 国产麻豆成人av免费视频| 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 美女大奶头视频| 婷婷丁香在线五月| 免费av毛片视频| 最新美女视频免费是黄的| 亚洲成av片中文字幕在线观看| 国产激情欧美一区二区| 国产一区二区三区在线臀色熟女| 久久精品亚洲熟妇少妇任你| 视频在线观看一区二区三区| 999精品在线视频| 好男人在线观看高清免费视频 | 亚洲狠狠婷婷综合久久图片| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 夜夜看夜夜爽夜夜摸| 69精品国产乱码久久久| 午夜免费成人在线视频| 亚洲第一欧美日韩一区二区三区| 色综合亚洲欧美另类图片| 久久人人爽av亚洲精品天堂| 午夜福利免费观看在线| 色婷婷久久久亚洲欧美| 亚洲人成电影观看| 日本免费一区二区三区高清不卡 | 国产蜜桃级精品一区二区三区| 九色国产91popny在线| 天堂√8在线中文| 在线国产一区二区在线| 一夜夜www| 欧美不卡视频在线免费观看 | 免费人成视频x8x8入口观看| 精品不卡国产一区二区三区| 欧美色视频一区免费| 乱人伦中国视频| 亚洲第一欧美日韩一区二区三区| 午夜免费观看网址| 久久国产亚洲av麻豆专区| 久久中文字幕一级| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| 婷婷精品国产亚洲av在线| 婷婷丁香在线五月| 日韩欧美三级三区| 亚洲av片天天在线观看| 国产高清有码在线观看视频 | 国产成人av激情在线播放| 国产精品日韩av在线免费观看 | 日韩免费av在线播放| 欧美成狂野欧美在线观看| 成人国产一区最新在线观看| 欧美成人免费av一区二区三区| 色综合亚洲欧美另类图片| 脱女人内裤的视频| 少妇裸体淫交视频免费看高清 | 最新美女视频免费是黄的| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品一区二区www| 免费人成视频x8x8入口观看| 亚洲精品中文字幕在线视频| 女人精品久久久久毛片| 他把我摸到了高潮在线观看| 久久 成人 亚洲| 日本a在线网址| www日本在线高清视频| 国产精品av久久久久免费| 亚洲av成人一区二区三| 啪啪无遮挡十八禁网站| 亚洲欧美激情在线| 嫁个100分男人电影在线观看| 国语自产精品视频在线第100页| 如日韩欧美国产精品一区二区三区| 黄色 视频免费看| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| 久久青草综合色| 久久中文字幕一级| 久久香蕉精品热| 无人区码免费观看不卡| av欧美777| 中文字幕久久专区| 51午夜福利影视在线观看| 又紧又爽又黄一区二区| 热99re8久久精品国产| 国产不卡一卡二| 午夜免费鲁丝| 免费女性裸体啪啪无遮挡网站| 一级黄色大片毛片| 国产成人一区二区三区免费视频网站| av欧美777| 啦啦啦韩国在线观看视频| 免费女性裸体啪啪无遮挡网站| 欧美成人午夜精品| x7x7x7水蜜桃| 免费av毛片视频| 男男h啪啪无遮挡| 欧美在线一区亚洲| 一进一出抽搐gif免费好疼| 久久久久国产一级毛片高清牌| 精品乱码久久久久久99久播| 国产亚洲精品av在线| 久久天躁狠狠躁夜夜2o2o| 亚洲专区中文字幕在线| 亚洲av电影不卡..在线观看| www.999成人在线观看| 熟妇人妻久久中文字幕3abv| 在线观看免费日韩欧美大片| 国产三级在线视频| 久久影院123| 国产精品精品国产色婷婷| 又紧又爽又黄一区二区| 老司机午夜十八禁免费视频| 欧美午夜高清在线| 韩国av一区二区三区四区| 这个男人来自地球电影免费观看| 嫩草影视91久久| 多毛熟女@视频| 琪琪午夜伦伦电影理论片6080| 曰老女人黄片| 国产国语露脸激情在线看| 日本欧美视频一区| 国产精品香港三级国产av潘金莲| 中文字幕久久专区| 满18在线观看网站| 午夜视频精品福利| 国产激情久久老熟女| 国产三级黄色录像| 九色亚洲精品在线播放| 久久久精品国产亚洲av高清涩受| 日韩欧美一区二区三区在线观看| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 99国产极品粉嫩在线观看| 黄色丝袜av网址大全| 嫩草影院精品99| 亚洲熟女毛片儿| 日本vs欧美在线观看视频| 久久天堂一区二区三区四区| 搡老熟女国产l中国老女人| 午夜免费激情av| 村上凉子中文字幕在线| 夜夜躁狠狠躁天天躁| 亚洲精品中文字幕一二三四区| 黄色 视频免费看| 国产亚洲欧美精品永久| 国产精品日韩av在线免费观看 | 国产成人av教育| 日韩欧美一区视频在线观看| 亚洲av成人一区二区三| 男人的好看免费观看在线视频 | 亚洲 国产 在线| 老司机午夜福利在线观看视频| 一二三四社区在线视频社区8| 亚洲自偷自拍图片 自拍| 好男人在线观看高清免费视频 | 久久天躁狠狠躁夜夜2o2o| 99国产极品粉嫩在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩欧美免费精品| 黄色 视频免费看| 免费久久久久久久精品成人欧美视频| 成人免费观看视频高清| 欧美成人午夜精品| 可以免费在线观看a视频的电影网站| 一二三四社区在线视频社区8| 色老头精品视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一区福利在线观看| 中文字幕高清在线视频| 久久久久九九精品影院| 巨乳人妻的诱惑在线观看| 国产1区2区3区精品| 视频在线观看一区二区三区| 亚洲成人精品中文字幕电影| 亚洲色图综合在线观看| 亚洲国产欧美一区二区综合| 亚洲精品一区av在线观看| 50天的宝宝边吃奶边哭怎么回事| cao死你这个sao货| 国产精品影院久久| 黄片播放在线免费| 亚洲国产精品成人综合色| 99国产极品粉嫩在线观看| 亚洲在线自拍视频| 啪啪无遮挡十八禁网站| 女人高潮潮喷娇喘18禁视频| 国产精品一区二区免费欧美| 99国产精品一区二区蜜桃av| 伊人久久大香线蕉亚洲五| 国产日韩一区二区三区精品不卡| 狠狠狠狠99中文字幕| 黑人巨大精品欧美一区二区蜜桃| 两个人看的免费小视频| 亚洲国产精品成人综合色| 婷婷精品国产亚洲av在线| 亚洲欧美日韩高清在线视频| av网站免费在线观看视频| 亚洲伊人色综图| 50天的宝宝边吃奶边哭怎么回事| av网站免费在线观看视频| 国产野战对白在线观看| 亚洲色图综合在线观看| 日韩大码丰满熟妇| 美女午夜性视频免费| 亚洲国产欧美网| 亚洲色图 男人天堂 中文字幕| 国产主播在线观看一区二区| 国产伦一二天堂av在线观看| 在线观看免费午夜福利视频|