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

    黑土旱地改稻田土壤水穩(wěn)性團聚體有機碳和全氮的變化特征

    2020-08-15 03:10:10馬原遲美靜張玉玲范慶峰虞娜鄒洪濤
    中國農(nóng)業(yè)科學(xué) 2020年8期
    關(guān)鍵詞:旱地粒級土壤有機

    馬原,遲美靜,張玉玲, 范慶峰,虞娜,鄒洪濤

    黑土旱地改稻田土壤水穩(wěn)性團聚體有機碳和全氮的變化特征

    馬原,遲美靜,張玉玲, 范慶峰,虞娜,鄒洪濤

    (沈陽農(nóng)業(yè)大學(xué)土地與環(huán)境學(xué)院/農(nóng)業(yè)農(nóng)村部東北耕地保育重點實驗室,沈陽 110866)

    【目的】分析東北黑土旱地改稻田后土壤團聚體組成及其穩(wěn)定性、各粒級團聚體有機碳、全氮含量及其13C、15N自然豐度值的動態(tài)變化,探討旱地改稻田后土壤團聚體有機碳、全氮的賦存能力及穩(wěn)定性,揭示旱地改稻田后土壤團聚體及其有機碳、全氮的演變規(guī)律?!痉椒ā窟x擇東北典型黑土旱地土壤(種植大豆年限大于60年,作為對照)和改種不同年限的稻田土壤(3、5、10、17、20和25年,改稻田前種植作物均為大豆),利用土壤團聚體濕篩分離技術(shù)和穩(wěn)定同位素分析技術(shù),研究旱地改稻田后土壤團聚體有機碳、全氮的動態(tài)變化特征?!窘Y(jié)果】在0—60 cm土層,與對照土壤相比,改種水稻各年限土壤中2—0.25 mm團聚體組成有所減少,0.25—0.053 mm和<0.053 mm團聚體組成有所增加,>2 mm團聚體組成的變化無明顯規(guī)律,但旱地改稻田不同年限均以2—0.053 mm團聚體為主;團聚體平均重量直徑(MWD)與>2 mm團聚體組成之間呈顯著線性正相關(guān)關(guān)系(<0.01),與0.25—0.053 mm、<0.053 mm團聚體組成之間均呈顯著線性負相關(guān)關(guān)系(<0.01或<0.05);水穩(wěn)性團聚體組成變化受水稻種植年限和土層深度的顯著影響,而MWD的變化則受土層深度的顯著影響。與對照土壤相比,在0—40 cm土層,2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降,在改種水稻3—25年間均隨水稻種植年限延長大體上呈增加趨勢。總體上,2—0.25 mm、0.25—0.053 mm團聚體是賦存有機碳和全氮的主要粒級;在0—60 cm土層,>2 mm團聚體有機碳、全氮含量與其團聚體組成之間呈顯著正相關(guān)關(guān)系(<0.01或<0.05),在0—20 cm土層,2—0.25 mm團聚體有機碳、全氮含量與其團聚體組成之間也呈顯著正相關(guān)關(guān)系(<0.01或<0.05);<2 mm團聚體有機碳和全氮含量的變化受水稻種植年限影響顯著,而>0.25 mm團聚體有機碳和全氮含量的變化則受土層深度影響顯著。與對照土壤相比,各粒級團聚體中δ13C在改種水稻3年時均明顯增加,在改種水稻5年時均明顯下降,在改種水稻5—25年間變化不明顯,各粒級團聚體中δ15N在改種水稻25年間均略有下降??傮w上,在改稻田3—25年間,團聚體中δ13C、δ15N的變化受水稻種植年限和土層深度的顯著影響,其數(shù)值均隨粒級的減少而增加,相同年限各粒級團聚體δ13C隨著土層的加深而增大,δ15N無明顯變化規(guī)律?!窘Y(jié)論】東北典型黑土旱地改稻田25年間,土壤中非水穩(wěn)性大團聚體遭受破壞形成了粒徑較小的團聚體,2—0.053 mm水穩(wěn)性團聚體是有機碳、全氮固存的主要載體,較小粒級團聚體賦存的有機碳較為穩(wěn)定,其穩(wěn)定性隨水稻種植年限延長、土層加深而增強。

    黑土;旱地;稻田;水穩(wěn)性團聚體;有機碳;全氮;13C和15N 自然豐度

    0 引言

    【研究意義】土壤團聚體是土壤結(jié)構(gòu)的基礎(chǔ),是由土壤顆粒和各種有機質(zhì)(微生物、植物、動物殘渣及其分泌物)膠結(jié)而成,不同粒級團聚體的分布影響著土壤養(yǎng)分的儲存、土壤生物學(xué)特性和土壤結(jié)構(gòu)[1-2]。因此,土壤團聚體的穩(wěn)定性和分布與土壤有機碳和全氮的固存密切相關(guān)[3],對于研究土壤結(jié)構(gòu)和有機碳、全氮的固存和轉(zhuǎn)化有重要意義[4]。東北黑土有機質(zhì)含量豐富、結(jié)構(gòu)良好。近年來,隨著農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整,黑土旱地改稻田面積逐年增加,這一土地利用方式的改變及土壤環(huán)境條件的變化,影響了土壤有機碳、全氮的數(shù)量[5],進而影響土壤中各粒級水穩(wěn)性團聚體的數(shù)量、分配及其穩(wěn)定性[6]。不同粒級團聚體及其有機碳、全氮的賦存數(shù)量的變化可以用來表征或反映不同環(huán)境條件下土壤質(zhì)量的變化[7-8],因此,研究東北典型黑土旱地改稻田后土壤水穩(wěn)性團聚體及其有機碳、全氮的動態(tài)變化對于揭示黑土有機碳、全氮的演變特征具有重要意義?!厩叭搜芯窟M展】水稻土是長期人為水耕熟化的一類特殊類型土壤,具有較高的固碳潛力。許多研究表明,長期植稻有利于土壤有機碳的固存[9-10],且與旱地土壤相比有機碳含量較高[5,11-14],水稻種植過程主要通過物理保護和腐殖化過程來富集和固存土壤有機碳[10],通過土壤團聚體對有機碳的保護機制可以解釋土壤有機碳的固存和分解[15]。研究發(fā)現(xiàn),太湖地區(qū)水稻土中主要以2—0.25 mm、0.25—0.02 mm團聚體為主[15];浙江平原區(qū)植稻17年土壤中2—0.25 mm團聚體是有機碳的主要載體[16];紅壤荒地開墾種植水稻20年后,水穩(wěn)性大團聚體是有機碳、全氮的主要載體[17];慈溪和上虞地區(qū)植稻300年間,大團聚體有所增多,微團聚體減少,大團聚體是土壤有機碳的主要載體[8]。但也有研究發(fā)現(xiàn),慈溪地區(qū)植稻30—2 000年間,0.25—0.053 mm團聚體逐漸減少,<0.053 mm團聚體增加,土壤有機碳主要分布在0.25—0.053 mm團聚體中[18]。稻田轉(zhuǎn)變?yōu)槭卟说?0年后>0.25 mm團聚體含量及其固碳能力會下降[19];本課題組前期研究發(fā)現(xiàn),東北黑土植稻大于50年的水稻土中主要以2—0.25 mm團聚體為主[20]。土壤13C、15N自然豐度值可以從內(nèi)在反映出土壤有機碳的來源和周轉(zhuǎn)速率[21-24],植被的轉(zhuǎn)變或者新碳的介入會對土壤有機碳、氮的周轉(zhuǎn)產(chǎn)生影響,從而導(dǎo)致各級團聚體δ13C、δ15N發(fā)生改變[25-27]。由此可見,不同土壤、土地利用方式及種植年限對土壤團聚體的影響不盡相同,多數(shù)研究集中在我國亞熱帶地區(qū)。【本研究切入點】課題組在前期研究中已明確東北黑土旱地改稻田后土壤有機碳、全氮的變化特征[11],但關(guān)于東北黑土旱田改稻田后土壤各粒級團聚體有機碳、全氮的變化特征如何?還尚不清楚?!緮M解決的關(guān)鍵問題】利用土壤團聚體濕篩分離技術(shù)和穩(wěn)定同位素分析技術(shù),研究東北黑土區(qū)旱地改稻田后土壤各粒級團聚體組成、團聚體中有機碳、全氮含量及其δ13C、δ15N的變化,闡明黑土區(qū)旱地改稻田后土壤團聚體有機碳、全氮的賦存能力及其穩(wěn)定性,為揭示黑土區(qū)旱地改稻田后土壤團聚體有機碳、全氮的演變規(guī)律提供重要的理論依據(jù)。

    1 材料與方法

    1.1 供試土壤

    供試土壤采自黑龍江省綏化市慶安縣勤勞鎮(zhèn)勤樸村,屬于典型黑土,為黃土性沉積物發(fā)育的土壤。2015 年10月進行實地田間生產(chǎn)調(diào)研,確定連續(xù)種植大豆60年以上的旱地土壤(作為對照土壤)和旱地改稻田年限分別為3、5、10、17、20 和 25年的稻田土壤。供試土壤區(qū)域的氣候和地形(漫崗丘陵區(qū),坡度小于3°)大致相同,旱地改種水稻前種植作物均為大豆,改種水稻前土壤性狀與對照土壤性狀基本一致。旱地改種水稻后每年施用氮、磷、鉀化肥量大致相同,但年限間化學(xué)肥料的品種、用量存在差異,供試土壤不施用有機肥料。

    2015年10月水稻收獲后進行供試田塊土壤樣本采集。將每個相同年限的田塊作為一個采樣區(qū)域,并劃分為面積大約相同的3個采樣單元作為3次采樣重復(fù),每個采樣單元以“S”法布設(shè)5—7點,用土鉆采集布設(shè)點樣本后均勻混合作為一個重復(fù)樣本,混合后土壤樣本裝在塑料箱中(避免破壞原土狀態(tài))??紤]黑土自身腐質(zhì)層比較深厚,同時旱地改稻田年限相對較短,稻田土壤犁底層形成相對較弱,因此確定采樣深度為0—20、20—40和40—60 cm。土樣運回室內(nèi)后,在室溫下自然風(fēng)干,在風(fēng)干過程中用手輕輕地把土塊沿自然結(jié)構(gòu)面掰成直徑約5 cm的小塊,并除去肉眼可以識別的粗根和石塊,風(fēng)干后混合均勻,其中,一部分土樣磨碎過0.15 mm篩備用,一部分土樣過8 mm篩用于制備各粒級團聚體。供試土壤的地理信息及理化性質(zhì)見表1。

    1.2 各級團聚體的分離

    土壤團聚體的分離依據(jù)SIX等[28]提供的土壤團聚體濕篩法,采用團粒分析儀(DIK-2012,日本)進行篩分。稱取50 g自然風(fēng)干土樣,放在2 mm、0.25 mm和0.053 mm組成的套篩上,調(diào)整水面高度使水面剛好沒過土樣,浸泡5 min后,以速度為30 次/min振蕩30 min,篩分結(jié)束后分別將套篩上的團聚體進行收集并轉(zhuǎn)移至蒸發(fā)皿中,并于60 ℃烘干至恒重(約48 h)并稱重,獲得>2 mm、2—0.25 mm、0.25—0.053 mm和<0.053 mm團聚體。并將制備好各粒級團聚體樣品磨碎過0.15 mm篩備用。

    1.3 有機碳、全氮含量及其δ13C、δ15N自然豐度值測定

    原土和各粒級團聚體樣品中的有機碳、全氮含量及其δ13C、δ15N均采用元素分析儀—穩(wěn)定同位素比例質(zhì)譜儀(EA–IRMS,Elementanalysis-Stable100 Isotope Ratio Mass Spectrometer,德國)聯(lián)用測定。

    式中,R、R分別為樣品13C/12C 原子比值和15N/14N原子比值;RR值為0.0112372(以美國南卡羅來納洲白堊紀皮狄組層中的擬箭石化石(Pee Dee Belemnite, PDB)為標準物質(zhì))和0.0036765(以純凈大氣氮(Air Nitrogen, AN)為標準物質(zhì))。

    1.4 數(shù)據(jù)計算

    (1)團聚體平均重量直徑(MWD)

    (2)某粒級團聚體有機碳、全氮含量

    aggregate(aggregate)=con-aggregate)con-aggregate)×aggregate(4)

    式中,aggregate(aggregate)為某粒級團聚體有機碳(全氮)含量(g·kg-1soil),con-aggregate(con-aggregate)為某粒組團聚體有機碳(全氮)測試濃度(g·kg-1aggregate)。

    1.5 數(shù)據(jù)處理

    采用 Excel 2010 和SPSS 19.0 (SPSS Inc., Chicago, IL, USA)進行數(shù)據(jù)處理,采用 Origin9.0 軟件對數(shù)據(jù)進行繪圖。以年限和土層為因素進行雙因素方差分析(two-ways ANOVA),采用Tukey法進行多重比較,若數(shù)據(jù)誤差方差為非齊性,則使用非參數(shù)檢驗(non-parametric test),Kruskal-Wallis單因素方差分析(one-way ANOVA)。若兩個因素存在顯著的交互作用,則進行單因素方差分析(one-way ANOVA),采用Tukey法或Dunnett T3法進行多重比較,并以小寫字母表示差異顯著性(<0.05)。數(shù)據(jù)均為3次重復(fù)的平均值。

    表1 供試土壤的地理信息及理化性質(zhì)

    2 結(jié)果

    2.1 旱地改稻田土壤水穩(wěn)性團聚體組成及平均重量直徑的變化

    由圖1可以看出,旱地改稻田3—25年間,在0—20 cm土層,>2 mm、2—0.25 mm、0.25—0.053 mm和<0.053 mm團聚體組成分別為8.96%—27.11%、27.68%—38.58%、22.29%—37.12%和14.65%—26.05%;在20—40 cm土層,>2 mm、2—0.25 mm、0.25—0.053 mm和<0.053 mm團聚體組成分別為2.87%— 24.37%、24.45%—34.40%、27.38%—38.12%和18.25% —32.32%;在40—60 cm土層,>2 mm、2—0.25 mm、0.25—0.053 mm和<0.053 mm團聚體組成分別為1.84%—7.54%、25.69%—34.41%、33.35%—38.91%和23.00%—33.42%。總體上,在0—60 cm土層,與種植大豆的旱地土壤相比,改種水稻各年限土壤中,2—0.25 mm團聚體組成有所減少,0.25—0.053 mm和<0.053mm團聚體組成有所增加,>2 mm團聚體組成的變化無明顯規(guī)律。另外,水稻種植年限和土層深度之間的交互作用對各粒級團聚體組成均無顯著影響,但水稻種植年限對2—0.25 mm、0.25— 0.053 mm和<0.053 mm團聚體組成影響顯著,土層深度則對>2 mm、0.25—0.053 mm和<0.053 mm團聚體組成影響顯著。由此可見,旱地改稻田,由于每年在水稻種植期間淹水,使土壤中非水穩(wěn)性大團聚體遭受破壞,形成了粒徑較小的團聚體,但均以2—0.25 mm和0.25—0.053 mm團聚體為主;水稻種植年限和土層深度是影響水穩(wěn)性團聚體組成變化的重要因素。

    數(shù)值為平均值±標準誤。* ,**,***,ns分別表示為0.05,0.01,0.001水平顯著和不顯著。下同

    由圖2所示,土層深度是顯著影響土壤團聚體平均重量直徑(MWD)變化的重要因素;總體上,MWD隨土層深度的加深而減少,0—20 cm土層的MWD顯著(<0.01或<0.05)大于20—40 cm和40—60 cm土層的MWD(<0.05)。相關(guān)分析顯示(表2),在3個土層中,MWD與>2 mm團聚體組成之間呈顯著線性正相關(guān)關(guān)系(<0.01),與0.25—0.053 mm、<0.053 mm團聚體組成之間均呈顯著線性負相關(guān)關(guān)系(<0.01或<0.05)。這表明旱地改稻田后,>2 mm團聚體數(shù)量的增加提高了土壤水穩(wěn)性團聚體的穩(wěn)定性,而較小粒級團聚體(<0.25 mm)的增加則降低了土壤水穩(wěn)性團聚體的穩(wěn)定性。

    ***為0.001水平顯著,ns為不顯著

    2.2 旱地改稻田土壤水穩(wěn)性團聚體有機碳(氮)含量的變化

    由表3可以看出,總體上,2—0.25 mm、0.25— 0.053 mm團聚體是賦存有機碳和全氮的主要粒級。與對照土壤相比,在0—20 cm和20—40 cm土層,2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降;在改種水稻3—25年間均隨水稻種植年限延長大體上呈增加趨勢,且在改種水稻25年時,兩個粒級團聚體有機碳和全氮含量與對照土壤的相接近。在40—60 cm土層,2—0.25 mm、0.25—0.053 mm團聚體有機碳和全氮含量在改種水稻3年時均有所下降;兩個粒級團聚體有機碳含量在改種水稻3—17年間大體呈逐年增加趨勢,而在改種水稻17—25年間則呈逐年降低趨勢;2—0.25 mm團聚體全氮含量在改種水稻3—25年間大體呈逐年增加趨勢,0.25—0.053 mm團聚體全氮含量在改種水稻3—17年間有所增加,在改種水稻17—25年間有所下降。>2 mm和<0.053 mm團聚體有機碳和全氮含量在改種水稻后無顯明規(guī)律性。另外,水稻種植年限和土層深度之間的交互作用對各粒級團聚體有機碳和全氮含量均無顯著影響,但水稻種植年限和土層深度對不同粒級團聚體有機碳和全氮含量有重要影響。總體上,水稻種植年限對2—0.25 mm、0.25—0.053 mm和<0.053 mm團聚體有機碳和全氮含量影響顯著,土層深度對>2 mm和2—0.25 mm團聚體有機碳和全氮含量影響顯著。

    表2 各粒級團聚體組成與土壤團聚體平均重量直徑的相關(guān)系數(shù)

    *,**分別表示為0.05,0.01水平顯著。下同

    With< 0.05(*),< 0.01(**)The same as below

    表3 旱地改稻田后不同年限土壤團聚體有機碳和全氮含量

    相關(guān)分析顯示(表4),在3個土層,>2 mm團聚體有機碳、全氮含量與其團聚體組成之間均呈顯著正相關(guān)關(guān)系(<0.01或<0.05)。0—20 cm土層,2—0.25 mm團聚體有機碳、全氮含量與其團聚體組成之間也呈顯著正相關(guān)關(guān)系(<0.01或<0.05)。這說明旱地改種水稻田后,有機碳、全氮對于土壤水穩(wěn)性大團聚體(>0.25 mm),尤其是表層水穩(wěn)性大團聚體(>2 mm)形成具有重要作用,同時大團聚體的形成與穩(wěn)定也有利于有機碳、全氮的賦存。

    表4 各粒級團聚體有機碳、全氮含量與各粒級土壤團聚體的相關(guān)系數(shù)

    2.3 旱地改稻田土壤水穩(wěn)性團聚體中δ13C、δ15N的變化

    由表5可見,與對照土壤相比,各粒級團聚體中δ13C在改種水稻3年時均明顯增加,在改種水稻3—5年間明顯下降,在改種水稻5—25年間變化不明顯;在改種水稻25年間,各粒級團聚體中δ15N值均略有下降??傮w上,在改種水稻3—25年間,團聚體中δ13C、δ15N隨粒級的減小呈增加趨勢,相同年限各粒級團聚體δ13C隨著土層的加深而增大,δ15N的變化無明顯規(guī)律。水稻種植年限顯著影響4個粒級團聚體中δ13C和δ15N,土層深度顯著影響<2 mm粒級團聚體中δ13C和δ15N,二者的交互作用僅顯著影響0.25—0.053 mm團聚體中δ13C和δ15N。這說明隨著團聚體的粒級變小,土壤有機碳(氮)易分解性越小,相對越穩(wěn)定;其穩(wěn)定性隨水稻種植年限延長、土層加深度而增強。

    3 討論

    3.1 旱地改稻田對土壤團聚體組成及其穩(wěn)定性的影響

    土地利用方式和耕作管理(作物類型、施肥、排水和灌溉)是影響土壤結(jié)構(gòu)和特性的最重要因素,并對土壤團聚體的形成和破裂產(chǎn)生影響[29]。Huang等[14]對東北旱地和稻田土壤的研究發(fā)現(xiàn),土壤團聚體主要以2—0.25 mm、0.25—0.053 mm粒級為主,旱地改種稻田27年時稻田土壤中>0.25 mm團聚體顯著高于旱地土壤。Zou等[9]研究表明,隨著種稻年限的延長大團聚體隨之增加,微團聚體隨之減少。李昌新等[11]研究表明,長期玉米和水稻種植下,稻田土壤中>2 mm、<0.053 mm團聚體均顯著高于旱地,而2—0.25 mm團聚體則相反。Zheng等[8]和Pan等[15]研究表明,耕作會引起2—0.25 mm團聚體發(fā)生改變。本研究中,旱地改稻田3年時,2—0.25 mm團聚體下降,<0.25 mm團聚體組成上升;旱地改稻田3—25年間,2—0.25 mm團聚體有所上升,<0.25 mm團聚體有所下降;旱地改稻田25年時2—0.25 mm團聚體低于旱地土壤,<0.25 mm團聚體高于旱地土壤(圖1),這說明黑土在旱地改稻田年限較短時,干濕交替過程[30]和淹水條件[31]會導(dǎo)致大團聚體破裂形成微團聚體,其中2—0.25 mm團聚體對于耕作環(huán)境的變化比較敏感,但隨著水稻種植年限的延長,由于長期淹水會使微團聚體進一步團聚,進而促進大團聚體的形成[15]。

    表5 旱地改稻田后不同年限土壤團聚體中的δ13C和δ15N

    不同小寫字母表示相同粒級團聚體,年限與土層之間存在顯著交互作用達0.05顯著水平

    Indicate the significant (<0.05) differences between years and soil depth by different lowercase letters

    土壤團聚體平均重量直徑(MWD)可以反映出不同粒級團聚體的分布狀況和穩(wěn)定性[1],MWD值越大,團聚體穩(wěn)定性越好[32]。Wang等[19]研究表明,大團聚體在干濕交替過程中易于破裂為微團聚體,并導(dǎo)致團聚體水穩(wěn)性顯著下降。本研究中,MWD值與>2 mm團聚體組成呈顯著正相關(guān),與<0.25 mm團聚體組成呈顯著負相關(guān)(表2);水稻種植年限對MWD影響不顯著,而土層深度則對MWD影響顯著,0—20 cm土層團聚體穩(wěn)定性顯著高于20—60 cm土層(圖2)。這說明旱地改稻田后,由于干濕交替和淹水過程,影響了0—60 cm土層各粒級團聚體的變化,水稻種植年限為17年時,由于0—40 cm土層中>2 mm和2—0.25 mm大團聚體的增加,0.25—0.053 mm和<0.053 mm團聚體的減少(圖1),使0—40 cm土層團聚體水穩(wěn)性有所上升(圖2)。

    3.2 旱地改稻田對土壤團聚體有機碳、全氮含量及其穩(wěn)定性的影響

    土壤團聚體是土壤有機碳、全氮固定[33]和物質(zhì)能量轉(zhuǎn)化[17]的重要場所,不僅可以保護土壤有機質(zhì)也可以延遲其礦化[34],由于各粒級團聚體對土壤有機碳、全氮的保護機制不同,因此對土壤有機碳、全氮的固持也存在差異[33]。大團聚體主要是通過根系、真菌菌絲等膠結(jié)而成。在大多數(shù)情況下,有機物先與其結(jié)合,碳濃度最高且周轉(zhuǎn)率最快[2],而微團聚體受芳香族化合物束縛,碳周轉(zhuǎn)速率較低[35]。李昌新等[10]對種植20多年的南方稻田和旱地土壤的研究表明,稻田和旱地土壤有機碳含量均主要分配在2—0.25 mm和0.25—0.053 mm團聚體中,與旱地土壤相比,稻田土壤降低了2—0.25 mm團聚體有機碳的比例,但提高了其他粒級團聚體有機碳的比例。HUANG等[14]研究表明,旱地轉(zhuǎn)變?yōu)榈咎?7年后顯著提高了2—0.25 mm、0.25—0.053 mm團聚體有機碳含量,但其他粒級團聚體有機碳含量則無顯著變化。徐文靜等[20]對東北黑土區(qū)種稻大于50年的水稻土研究表明,2—0.25 mm團聚體是賦存有機碳的主要載體,具有明顯的固碳能力。許多研究表明,稻田土壤中有機碳含量的變化與大團聚體緊密相關(guān)[10,15,36],各粒級團聚體有機碳含量隨種稻年限的延長而增加[9],但也有研究表明,稻田土壤有機碳可向微團聚體和粉黏粒中轉(zhuǎn)移以利于有機碳的長期賦存[11,14,18]。本研究中,旱地改種稻田后各粒級團聚體中有機碳、全氮含量的變化大致相同,土壤有機碳、全氮含量與土壤團聚體組成之間緊密相關(guān)(表4),在改種水稻年限較短時(3年間),各粒級團聚體中有機碳、全氮含量下降(表3),這主要是由于土地利用方式的改變使土壤中各粒級團聚體遭受不同程度破壞(圖1),進而加速了有機碳(氮)分解和周轉(zhuǎn)速率[37];在改稻田3—25年間,2—0.25 mm、0.25—0.05 mm團聚體有機碳、全氮含量則有所增加(圖3),這主要是由于隨著改種年限延長,淹水條件可延緩?fù)寥烙袡C質(zhì)的分解速率,降低土壤有機碳、氮的礦化速率[38],同時由于大團聚體破裂為微團聚體,使其微團聚體易與有機質(zhì)粘結(jié)[39]。由此可見,東北黑土旱地改稻田25間,2—0.25 mm、0.25—0.053 mm團聚體是土壤有機碳、全氮賦存的主要粒級,對有機碳(氮)的物理保護發(fā)揮重要作用。

    采用δ13C可以研究土壤有機碳的動態(tài)變化,土壤的δ13C隨時間的變化主要受植被輸入的有機碳控制,植物殘體進入土壤后在微生物作用下可發(fā)生礦化分解過程和腐殖化過程,致使土壤的δ13C會與植物本身的δ13C有所不同[24];其次由生物衰變過程控制[22],土壤有機碳分解過程中微生物對12C的優(yōu)先分解會引起土壤δ13C增高,δ13C越大,其有機碳的降解程度越低[25]。大豆和水稻均屬于C3作物,其δ13C介于-29.1‰—-24.2‰[21],水稻的平均δ13C為-26.34‰[4]。本研究中,長期連續(xù)種植大豆的旱地,在改種水稻年限較短時(3年),由于土地利用方式及植被類型的改變,新鮮有機質(zhì)會不同程度進入到不同粒級團聚體中,但又由于水稻單產(chǎn)顯著高于大豆單產(chǎn),對土壤養(yǎng)分的需求大幅度增加,同時大團聚體遭受破壞,加速土壤中易分解有機碳的分解使碳同位素分餾效應(yīng)增加,因此,致使旱地改稻田3年間各粒級團聚體中δ13C均有所增加(表5),這與課題組前期關(guān)于旱地改稻田土壤δ13C的變化結(jié)果相一致[5];在旱地改稻田5—25年,由于大量水稻根系和地上殘體歸還土壤使有機碳得到補充和更新,又由于稻田土壤環(huán)境逐漸趨于穩(wěn)定可使有機碳分解緩慢,進而使各級團聚體中δ13C均有所降低,但隨種稻年限延長變化較為平緩(表5)。LIU等[27]研究發(fā)現(xiàn),新鮮有機碳的賦存一般從大團聚體(>0.25 mm)開始,經(jīng)過分解過程和微生物消耗后,降解的有機碳則被螯合在微團聚體(<0.25 mm)中,植物殘留物的腐殖化作用會伴隨著13C的輕微富集。本研究中,團聚體的粒徑越小,其δ13C越大(表5),這也說明團聚體粒徑越大,賦存新鮮有機碳越多,這一結(jié)果與LIU等[27]研究結(jié)果相類似。另外,相同年限各粒級團聚體中δ13C隨土層加深而增大(表5),這進一步說明,外源新鮮有機碳不能及時地進入補充到深層土壤,同時深層土壤有機碳相對較“老”[23,27],礦化分解速度較慢。

    土壤中δ15N的變異可反映出土壤氮素的轉(zhuǎn)化和遷移特征,其中礦化和硝化過程是導(dǎo)致氮同位素分餾從而影響δ15N的主要原因[40]。研究表明,土壤中δ15N在很大程度受土壤水分的影響[40-41],非淹水條件下土壤中硝化作用較快,而在淹水條件下有利于反硝化作用,礦化作用相對較弱,微生物對氮的固定作用較快。本研究中,旱地改稻田后,由于稻田土壤多處于淹水或含水量較高的狀態(tài),土壤硝化作用受到抑制,大量15N貧化的水稻植株根系或殘體進入土壤[26,35],致使稻田土壤各粒級團聚體的δ15N低于旱地土壤(表5)。

    4 結(jié)論

    東北黑土旱地改種稻田25年間,稻田土壤中2—0.25 mm團聚體有所減少,0.25—0.053 mm和<0.053 mm團聚體有所增加,>2 mm團聚體則無明顯變化規(guī)律,但均以2—0.053 mm團聚體為主;各粒級團聚體有機碳、全氮含量隨改種年限延長呈大致相同的變化趨勢,總體上,2—0.053 mm團聚體是賦存有機碳和全氮的主要粒級;團聚體中δ13C、δ15N均隨粒級的減小呈增加趨勢,相同年限各粒級團聚體δ13C隨著土層的加深而增大,δ15N無明顯變化規(guī)律。因此,東北典型黑土旱地改稻田后,土壤中非水穩(wěn)性大團聚體遭受破壞形成了粒徑較小的團聚體,2—0.053 mm水穩(wěn)性團聚體是有機碳、全氮固存的主要載體,較小粒級團聚體賦存的有機碳較為穩(wěn)定,其穩(wěn)定性隨水稻種植年限延長、土層加深而增大。

    [1] WANG Y, ZHANG J H, ZHANG Z H. Influences of intensive tillage on water-stable aggregate distribution on a steep hillslope., 2015, 151: 82-92.

    [2] HUANG R, LAN M L, LIU J, GAO M. Soil aggregate and organic carbon distribution at dry land soil and paddy soil: the role of different straws returning., 2017, 24(36): 1-11.

    [3] SIX J, BOSSUYT H, DEGRYZE S, DENEF K. A history of research on the link between (Micro) aggregates, soil biota, and soil organic matter dynamics., 2004, 79(1): 7-31.

    [4] LIU Y, HU C, HU W, WANG L, LI Z G, PAN J F, CHEN F. Stable isotope fractionation provides information on carbon dynamics in soil aggregates subjected to different long-term fertilization practices., 2018, 177: 54-60.

    [5] 賈樹海, 張佳楠, 張玉玲. 東北黑土區(qū)旱田改稻田后土壤有機碳、全氮的變化特征, 中國農(nóng)業(yè)科學(xué), 2017, 50(7): 1252-1262.

    JIA S H, ZHANG J N, ZHANG Y L. Changes of the characteristics of soil organic carbon and total nitrogen after conversation from upland to paddy field in black soil region of Northeast China., 2017, 50(7): 1252-1262. (in Chinese)

    [6] ZHAO J S, CHEN S, HU R G, LI Y Y. Aggregate stability and size distribution of red soils under different land uses integrally regulated by soil organic matter, and iron and aluminum oxides., 2017, 167: 73-79.

    [7] 徐香茹, 汪景寬. 土壤團聚體與有機碳穩(wěn)定機制的研究進展. 土壤通報, 2017, 48(6): 1523-1529.

    XU X R, WANG J K. A review on different stabilized mechanisms of soil aggregates and organic carbon.2017, 48(6): 1523- 1529. (in Chinese)

    [8] ZHENG H B, LIU W R, ZHENG J Y, LUO Y, LI R P, WANG H, QI H. Effect of long-term tillage on soil aggregates and aggregate-associated carbon in black soil of Northeast China,, 2018, 13(6): e199523.

    [9] ZOU P, FU J R, CAO Z H, YE J, YU Q G. Aggregate dynamics and associated soil organic matter in topsoils of two 2,000-year paddy soil chronosequences., 2015, 15(3): 510-522.

    [10] WANG P, LIU Y L, LI L Q, CHENG K, ZHENG J F, ZHANG X H, ZHENG J W, JOSEPH S, PAN G X. Long-term rice cultivation stabilizes soil organic carbon and promotes soil microbial activity in a salt marsh derived soil chronosequence., 2015, 5: 15704.

    [11] 李昌新, 黃山, 彭現(xiàn)憲, 黃欠如, 張衛(wèi)建. 南方紅壤稻田與旱地土壤有機碳及其組分的特征差異. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2009, 28(3): 606-611.

    LI C X, HUANG S, PENG X X, HUANG Q R, ZHANG W J. Differences in soil organic carbon fractions between paddy field and upland field in red soil region of South China.2009, 28(3): 606-611. (in Chinese)

    [12] YAN X, ZHOU H, ZHU Q H, WANG X F, ZHANG Y Z, YU X C, PENG X. Carbon sequestration efficiency in paddy soil and upland soil under long-term fertilization in Southern China., 2013, 130: 42-51.

    [13] SUN Y N, HUANG S, YU X C, ZHANG W J. Differences in fertilization impacts on organic carbon content and stability in a paddy and an upland soil in Subtropical China., 2015, 397(1): 189-200.

    [14] HUANG S, PAN X H, GUO J, QIAN C R, ZHANG W J. Differences in soil organic carbon stocks and fraction distributions between rice paddies and upland cropping systems in China., 2014, 14(1): 89-98.

    [15] PAN G X, WU L S, LI L Q, ZHANG X H, GONG W, WOOD Y. Organic carbon stratification and size distribution of three typical paddy soils from Taihu lake region, China., 2008, 20(4): 456-463.

    [16] 毛霞麗, 陸扣萍, 何麗芝, 宋照亮, 徐祖祥, 楊文葉, 徐進, 王海龍. 長期施肥對浙江稻田土壤團聚體及其有機碳分布的影響. 土壤學(xué)報, 2015(4): 828-838.

    MAO X L, LU K P, HE L Z, SONG Z L, XU Z X, YANG W Y, XU J, WANG H L. Effect of long-term fertilizer application on distribution of aggregate-associated organic carbon in paddy doil., 2015(4): 828-838. (in Chinese)

    [17] 陳曉芬, 李忠佩, 劉明, 江春玉. 不同施肥處理對紅壤水稻土團聚體有機碳、氮分布和微生物生物量的影響. 中國農(nóng)業(yè)科學(xué), 2013, 46(5): 950-960.

    CHEN X F, LI Z P, LIU M, JIANG C Y. Effects of different fertilizations on organic carbon and nitrogen contents in water-stable aggregates and microbial biomass content in paddy soil of Subtropical China., 2013, 46(5): 950-960. (in Chinese)

    [18] 王欣欣, 符建榮, 鄒平, 陳維, 葉靜, 俞巧鋼, 姜麗娜, 王強. 長期植稻年限序列水稻土團聚體有機碳分布特征. 應(yīng)用生態(tài)學(xué)報, 2013, 24(3): 719-724.

    WANG X X, FU J R, ZOU P, CHEN W, YE J, YU Q G, JIANG L N, WANG Q. Distribution characteristics of aggregates organic carbon in a paddy soil chronosequence., 2013, 24(3): 719-724. (in Chinese)

    [19] WANG H, GUAN D S, ZHANG R D, CHEN Y J, HU Y T, XIAO L. Soil aggregates and organic carbon affected by the land use change from rice paddy to vegetable field., 2014, 70: 206-211.

    [20] 徐文靜, 叢耀輝, 張玉玲, 段鵬鵬, 范慶鋒, 張玉龍. 黑土區(qū)水稻土水穩(wěn)性團聚體有機碳及其顆粒有機碳的分布特征. 水土保持學(xué)報, 2016, 30(4): 210-215.

    XU W J, CONG Y H, ZHANG Y L, DUAN P P, FAN Q F, ZHANG Y L. Distribution of organic carbon and partivulate organic carbon in water-stable aggregates of paddy soil in black soil area., 2016, 30(4): 210-215. (in Chinese)

    [21] 竇森, 張晉京. 用δ13C值研究土壤有機質(zhì)周轉(zhuǎn)的方法及其評價. 吉林農(nóng)業(yè)大學(xué)學(xué)報, 2001, 23(2): 64-67.

    DOU S, ZHANG J J. Introduction of a method for studying turnover of soil organic matter., 2001, 23(2): 64-67. (in Chinese)

    [22] BAI E, BOUTTON T W, LIU F, WU X B, HALLMARK C T, ARCHER S R. Spatial variation of soil δ13C and its relation to carbon input and soil texture in a subtropical lowland woodland., 2012, 44(1): 102-112.

    [23] PERI P L, LADD B, PEPPER D A, BONSER S P, LAFFAN S W, AMELUNG W. Carbon (δ13C) and nitrogen (δ15N) stable isotope composition in plant and soil in Southern Patagonia's Native Forests., 2012, 18(1): 311-321.

    [24] KULSAWAT W, PORNTEPKASEMSAN B, NOCHIT P. Paddy soil profile dstribution of δ13C subjected to rice straw amendment and burning., 2019, 886: 3-7.

    [25] GUNINA A, KUZYAKOV Y. Pathways of Litter C by Formation of aggregates and SOM density fractions: implications from13C natural abundance., 2014, 71: 95-104.

    [26] LIM S S, KWAK J H, LEE K S, CHANG S X, YOON K S, KIM H Y, CHOI W J. Soil and plant nitrogen pools in paddy and upland ecosystems have contrasting δ15N., 2015, 51(2): 231-239.

    [27] LIU Y, LIU W Z, WU L H, LIU C, WANG L, CHEN F, LI Z G. Soil aggregate-associated organic carbon dynamics subjected to different types of land use: evidence from13C natural abundance., 2018, 122: 295-302.

    [28] SIX J, ELLIOTT E T, PAUSTIAN K, DORAN J W. Aggregation and soil organic matter accumulation in cultivated and native grassland soils., 1998, 62(5): 1367-1377.

    [29] ZHAO J S, CHEN S, HU R G, LI Y L. Aggregate stability and size distribution of red soils under different land uses integrally regulated by soil organic matter, and iron and aluminum oxides., 2017, 167: 73-79.

    [30] VAN VEEN A J, KUIKMAN P J. Soil structural aspects of decomposition of organic matter by micro-organisms., 1990, 11(3): 213-223.

    [31] JIANG X J, XIE D T. Combining ridge with no-tillage in lowland rice-based cropping system: long-term effect on soil and rice yield., 2009, 19(4): 515-522.

    [32] 郝翔翔, 楊春葆, 苑亞茹. 連續(xù)秸稈還田對黑土團聚體中有機碳含量及土壤肥力的影響. 中國農(nóng)學(xué)通報, 2013, 29(35): 263-269.

    HAO X X, YANG C B, YUAN Y R. Effects of continuous straw returning on organic carbon content in aggregates and fertility of black soil., 2013, 29(35): 263-269. (in Chinese)

    [33] 徐虎, 張敬業(yè), 蔡岸冬, 王小利, 張文菊. 外源有機物料碳氮在紅壤團聚體中的殘留特征. 中國農(nóng)業(yè)科學(xué), 2015, 48(23): 4660-4668.

    Xu H, ZHANG J Y, CAI A D,WANG X L, ZHANG W J. Retention characteristic of carbon and nitrogen from amendments in different size aggregates of red soil., 2015, 48(23): 4660-4668. (in Chinese)

    [34] RABBI S M F, WILSON B R, LOCKWOOD P V, DANIEL H, YOUNG I M. Aggregate hierarchy and carbon mineralization in two oxisols of new south wales, Australia., 2015, 146: 193-203.

    [35] 高崇升, 王建國. 黑土農(nóng)田土壤有機碳演變研究進展. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2011, 19(6): 1468-1474.

    Gao C S, Wang J G. A review of researches on evolution of soil organic carbon in mollisols farmland.2011, 19(6): 1468-1474. (in Chinese)

    [36] CHEN Z D, TI J S, CHEN F. Soil aggregates response to tillage and residue management in a double paddy rice soil of the Southern China., 2017, 109(9): 1-12.

    [37] WANG B S, GAO L L, YU W S, WEI X Q, LI J, LI S P, SONG X J, LIANG G P, CAI D X, WU X P. Distribution of soil aggregates and organic carbon in deep soil under long-term conservation tillage with residual retention in dryland., 2019, 11(2): 241-254.

    [38] K?GEL-KNABNER I, AMELUNG W, CAO Z H, FIEDLER S, FRENZEL P, JAHN R, KALBITZ K, K?IBI A, SCHLOTER M. Biogeochemistry of paddy soils., 2010, 157(1): 1-14.

    [39] QIAN J, LIU J J, WANG P F, WANG C, HU J, LI K, LU B, TIAN X, GUAN W Y. Effects of riparian land use changes on soil aggregates and organic carbon., 2018, 112: 82-88.

    [40] 慈恩, 楊林章, 倪九派, 高明, 謝德體. 不同區(qū)域水稻土的氮素分配及δ15N特征. 水土保持學(xué)報, 2009, 23(2): 103-108.

    CI E, YANG L Z, NI J P, GAO M, XIE D T. Distribution and δ15N character of nitrogen in paddy soils located in different regions., 2009, 23(2): 103-108. (in Chinese)

    [41] CHOI W J, RO H M. Differences in isotopic fractionation of nitrogen in water-saturated and unsaturated soils., 2003, 35(3): 483-486.

    Change Characteristics of Organic Carbon and Total Nitrogen in Water-stable Aggregate After Conversion from Upland to Paddy Field in Black Soil

    MA Yuan, CHI MeiJing, ZHANG YuLing, FAN QingFeng, YU Na, ZOU HongTao

    (College of Land and Environment, Shenyang Agricultural University/Key Laboratory of Northeast Arable Land Conservation, Ministry of Agriculture and Rural Affairs, Shenyang 110866)

    【Objective】 The objectives of this study wereto analyze the composition and stability of soil aggregate, the changes of organic carbon (OC), total nitrogen (TN) content,natural abundance of13C and 15N in different-sized aggregates, to explore the sequestration and stability of soil aggregate organic C, TN, and to reveal the evolution of soil aggregate organic C, TN changes after the conversion from upland to paddy field in black soil region of Northeast China.【Method】 Soil samples were collected from upland (soybean planted for over 60 years) in typical black soil and paddy soil with different years (3, 5, 10, 17, 20 and 25 years, soybean was planted in all the fields before conversion to paddy field). The dynamic characteristics of OC and TN in soil aggregates were studied by using wet-sieving method and stable isotope analysis technology.【Result】In the 0-60 cm soil layers, compared with the control treatment, the composition of 2-0.25 mm aggregates in the soil of different years after rice planting was decreased, which of 0.25-0.053 mm and <0.053 mm aggregates was increased. There was no obvious change in the composition of >2 mm aggregates, but the different years of dry land change to paddy fields were dominated by 2-0.053 mm aggregates; the mean weight diameter (MWD) of aggregates was significantly positive correlated with the proportion of >2 mm aggregates (<0.01), and significantly negative correlated with the proportion of 0.25-0.053 mm and <0.053 mm aggregates (<0.01 or<0.05). The change of aggregate composition was significantly affected by different rice planting years and soil depth, whereas the MWD was significantly affected by soil depth. Compared with the control soil, in the 0-40 cm soil layer, the OC and TN contents in the size of 2-0.25 mm and 0.25-0.053 mm aggregates were declined in the 3 years, however there showed increased trend with the extension of rice cultivation in 3-25 years. Generally, OC and TN were mainly accumulated in the 2-0.25 mm and 0.25-0.053 mm aggregates. There existed significant positive correlation between OC and TN contents and aggregate composition in > 2 mm aggregates (<0.01 or<0.05) in the 0-60 cm soil layers, as well as 2-0.25 mm aggregates in the 0-20cm soil layer (<0.01 or<0.05). The OC and TN contents variation in <2 mm aggregates were significantly affected by rice cultivation time, while soil depth significantly affected >0.25 mm aggregates OC and TN contents. Compared with the control soil, the δ13C in each size of aggregates significantly increased in 3 rice planting years and decreased in 5 rice planting years, respectively, while there was no significant change in the 5-25 rice planting years, and the δ15N in all size of aggregates decreased slightly during the 25 years of rice replanting. In general, the δ13C and δ15N of soil in aggregates were significantly affected by rice cultivation time and soil depth, which increased with the decreasing of aggregate size. The δ13C increased with soil depth in the same year, while δ15N had no significant change.【Conclusion】After the conversion from dry land to paddy field for 25 years, non-water-stable macro-aggregates in the soil were damaged and formed into small sized aggregates. The 2-0.053 mm water-stable aggregates were the main carrier of OC and TN sequestration, while OC in small size aggregates more stable, and its stability was increased by the rice cultivation time and soil depth increased.

    black soil; upland; paddy field; water-stable aggregate; organic carbon; total nitrogen;13C and15N natural abundance

    2019-07-12;

    2019-11-06

    國家自然科學(xué)基金面上項目(41571280)、國家自然科學(xué)基金青年基金項目(41101276)

    馬原,E-mail:mayuan275@163.com。通信作者張玉玲,E-mail:zhangyuling@syau.edu.cn

    (責任編輯 李云霞)

    猜你喜歡
    旱地粒級土壤有機
    國外某大型銅礦選礦廠流程考查與分析①
    礦冶工程(2022年6期)2023-01-12 02:15:10
    旱地麥田夏閑期復(fù)種綠肥 保持土壤的可持續(xù)生產(chǎn)力
    山地暗棕壤不同剖面深度的團聚體分布
    旱地冰球運動開展價值的研究
    冰雪運動(2020年6期)2020-07-21 03:28:16
    旱地冰球運動推廣及發(fā)展策略
    冰雪運動(2019年5期)2019-08-24 08:04:48
    西雙版納橡膠林土壤有機碳分布特征研究
    秸稈還田的土壤有機碳周轉(zhuǎn)特征
    土壤有機碳轉(zhuǎn)化研究及其進展
    不同粒級再生骨料取代的混凝土基本性能試驗研究
    AMDIS在土壤有機污染物鑒別中表征性統(tǒng)計量的探究
    777久久人妻少妇嫩草av网站| 欧美激情久久久久久爽电影| 99精品久久久久人妻精品| 亚洲性夜色夜夜综合| 亚洲av美国av| 日韩欧美一区二区三区在线观看| 欧美在线黄色| 亚洲国产精品999在线| 午夜福利18| 女人高潮潮喷娇喘18禁视频| 国产亚洲av高清不卡| 91老司机精品| 日韩欧美国产在线观看| 国产成人av教育| 久久九九热精品免费| 精品久久久久久久末码| 91成年电影在线观看| 免费在线观看完整版高清| av在线播放免费不卡| 欧美日本亚洲视频在线播放| 2021天堂中文幕一二区在线观 | 欧美激情 高清一区二区三区| 草草在线视频免费看| 亚洲第一av免费看| 久久人人精品亚洲av| 久久人人精品亚洲av| 久久久国产欧美日韩av| 国产aⅴ精品一区二区三区波| 麻豆成人午夜福利视频| 啦啦啦 在线观看视频| 国产精品一区二区精品视频观看| 久久久国产成人精品二区| 欧美性猛交黑人性爽| 久久欧美精品欧美久久欧美| 国产精品乱码一区二三区的特点| 中文在线观看免费www的网站 | cao死你这个sao货| 亚洲av第一区精品v没综合| 性色av乱码一区二区三区2| 久久中文字幕一级| 丝袜美腿诱惑在线| 日韩欧美免费精品| 国产成人精品久久二区二区91| 久久精品91无色码中文字幕| 欧美成人一区二区免费高清观看 | 欧美一级毛片孕妇| 亚洲av日韩精品久久久久久密| 久久欧美精品欧美久久欧美| 又黄又爽又免费观看的视频| 国产1区2区3区精品| 99热这里只有精品一区 | 国产精品一区二区免费欧美| 国产成人欧美在线观看| 可以在线观看的亚洲视频| 99久久精品国产亚洲精品| 少妇的丰满在线观看| 中亚洲国语对白在线视频| 亚洲一区二区三区色噜噜| 91大片在线观看| 国产精品电影一区二区三区| 中文字幕av电影在线播放| 国产成人精品久久二区二区免费| 女人爽到高潮嗷嗷叫在线视频| 白带黄色成豆腐渣| 国产成+人综合+亚洲专区| 好看av亚洲va欧美ⅴa在| 久久久久国产精品人妻aⅴ院| 国产精品一区二区精品视频观看| 亚洲精品中文字幕在线视频| 国产av在哪里看| 久久精品国产亚洲av香蕉五月| 日韩欧美三级三区| 黄色视频,在线免费观看| 男女之事视频高清在线观看| 久久久久九九精品影院| 日韩一卡2卡3卡4卡2021年| 成人亚洲精品av一区二区| 九色国产91popny在线| 啦啦啦观看免费观看视频高清| 免费观看精品视频网站| 国产精品亚洲美女久久久| 在线观看免费视频日本深夜| 99久久久亚洲精品蜜臀av| www.自偷自拍.com| 国产又爽黄色视频| 国产精品二区激情视频| 免费看美女性在线毛片视频| 国产黄色小视频在线观看| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 午夜福利在线观看吧| 久久午夜亚洲精品久久| 精品一区二区三区四区五区乱码| 精品午夜福利视频在线观看一区| 最近最新中文字幕大全电影3 | 97碰自拍视频| 日本 av在线| 亚洲精品粉嫩美女一区| 俺也久久电影网| 日本在线视频免费播放| 免费在线观看视频国产中文字幕亚洲| 人人妻,人人澡人人爽秒播| 在线永久观看黄色视频| or卡值多少钱| 一本久久中文字幕| 欧美另类亚洲清纯唯美| 在线观看舔阴道视频| 国产1区2区3区精品| 欧美日韩瑟瑟在线播放| 韩国精品一区二区三区| 男女下面进入的视频免费午夜 | 国产精品久久久av美女十八| 欧美+亚洲+日韩+国产| 国内揄拍国产精品人妻在线 | 国产成人av教育| 国产单亲对白刺激| 国产aⅴ精品一区二区三区波| 亚洲av片天天在线观看| 色老头精品视频在线观看| 亚洲美女黄片视频| 琪琪午夜伦伦电影理论片6080| 国产精品 国内视频| 国产亚洲欧美98| 国产精品亚洲av一区麻豆| 俄罗斯特黄特色一大片| 99国产精品一区二区蜜桃av| 午夜福利视频1000在线观看| 国产日本99.免费观看| 国产高清videossex| 窝窝影院91人妻| 岛国在线观看网站| 老司机靠b影院| 97人妻精品一区二区三区麻豆 | 亚洲性夜色夜夜综合| 麻豆av在线久日| 亚洲人成电影免费在线| 精品久久蜜臀av无| 午夜免费鲁丝| 久久久久国产精品人妻aⅴ院| 日韩免费av在线播放| 亚洲国产看品久久| АⅤ资源中文在线天堂| 日本a在线网址| 欧美色视频一区免费| 成人三级做爰电影| 欧美一级a爱片免费观看看 | 久久久久免费精品人妻一区二区 | 亚洲成av片中文字幕在线观看| 欧美zozozo另类| 午夜老司机福利片| 欧美日韩亚洲国产一区二区在线观看| 最新美女视频免费是黄的| 久久 成人 亚洲| 老司机在亚洲福利影院| 97碰自拍视频| 久久精品国产清高在天天线| 久久久久久久久久黄片| 欧美激情极品国产一区二区三区| 91av网站免费观看| 最近最新免费中文字幕在线| 制服诱惑二区| 国产视频内射| 亚洲成国产人片在线观看| 国产成人啪精品午夜网站| 淫妇啪啪啪对白视频| 一本大道久久a久久精品| 亚洲五月色婷婷综合| 亚洲精品在线观看二区| 日韩欧美一区二区三区在线观看| 变态另类丝袜制服| 免费在线观看完整版高清| 国产一卡二卡三卡精品| 少妇被粗大的猛进出69影院| 在线观看一区二区三区| 欧美在线黄色| 一区福利在线观看| 99re在线观看精品视频| 精华霜和精华液先用哪个| 天天躁狠狠躁夜夜躁狠狠躁| 国产野战对白在线观看| 欧美又色又爽又黄视频| 久久国产精品男人的天堂亚洲| 一a级毛片在线观看| 亚洲男人天堂网一区| 亚洲成av片中文字幕在线观看| 国产一区二区在线av高清观看| 亚洲av熟女| 男人操女人黄网站| 亚洲男人的天堂狠狠| 精品国产国语对白av| 国产亚洲av高清不卡| 在线视频色国产色| 美女高潮喷水抽搐中文字幕| 国产精品自产拍在线观看55亚洲| 亚洲一码二码三码区别大吗| 国产精品av久久久久免费| 亚洲精品美女久久久久99蜜臀| 国产一区二区激情短视频| 日韩大码丰满熟妇| 午夜福利成人在线免费观看| 香蕉久久夜色| 亚洲专区中文字幕在线| 无遮挡黄片免费观看| 国产一卡二卡三卡精品| 日本精品一区二区三区蜜桃| av有码第一页| 国产色视频综合| 18美女黄网站色大片免费观看| av片东京热男人的天堂| 欧美一级a爱片免费观看看 | 亚洲成国产人片在线观看| 色综合欧美亚洲国产小说| √禁漫天堂资源中文www| 正在播放国产对白刺激| 久久热在线av| 免费观看精品视频网站| 热re99久久国产66热| 91在线观看av| 日韩精品中文字幕看吧| 日本成人三级电影网站| 中文亚洲av片在线观看爽| 老司机靠b影院| 99国产极品粉嫩在线观看| 一级a爱片免费观看的视频| 日韩欧美一区视频在线观看| 国产一卡二卡三卡精品| 国产亚洲欧美在线一区二区| 听说在线观看完整版免费高清| 日日干狠狠操夜夜爽| 1024视频免费在线观看| 最近在线观看免费完整版| 免费在线观看亚洲国产| 成年女人毛片免费观看观看9| 欧美日韩亚洲国产一区二区在线观看| 国产片内射在线| 99久久99久久久精品蜜桃| 777久久人妻少妇嫩草av网站| 亚洲国产欧美日韩在线播放| 国产精品久久久人人做人人爽| 99国产综合亚洲精品| 欧美日韩福利视频一区二区| 最近最新中文字幕大全电影3 | 人人妻,人人澡人人爽秒播| 亚洲av成人av| 色精品久久人妻99蜜桃| 非洲黑人性xxxx精品又粗又长| 免费看日本二区| 丰满人妻熟妇乱又伦精品不卡| 精品卡一卡二卡四卡免费| 悠悠久久av| 国产精品,欧美在线| 搡老妇女老女人老熟妇| 桃红色精品国产亚洲av| 91麻豆av在线| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 亚洲精品一区av在线观看| 欧美激情极品国产一区二区三区| 国产色视频综合| 最近最新中文字幕大全电影3 | 丝袜在线中文字幕| 国产又色又爽无遮挡免费看| 免费高清在线观看日韩| 在线观看免费日韩欧美大片| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 在线播放国产精品三级| 日日夜夜操网爽| 欧美国产精品va在线观看不卡| 麻豆成人av在线观看| 日韩精品青青久久久久久| 国产精品香港三级国产av潘金莲| 国产欧美日韩精品亚洲av| 免费看美女性在线毛片视频| 18美女黄网站色大片免费观看| 一进一出抽搐动态| 国产精品久久电影中文字幕| 日日夜夜操网爽| 欧美不卡视频在线免费观看 | 又黄又爽又免费观看的视频| 国产99白浆流出| 日韩成人在线观看一区二区三区| 免费av毛片视频| a级毛片a级免费在线| 午夜福利在线观看吧| 美女 人体艺术 gogo| 亚洲av熟女| 老鸭窝网址在线观看| 免费看日本二区| 人妻丰满熟妇av一区二区三区| 欧美日韩瑟瑟在线播放| 搡老岳熟女国产| 91av网站免费观看| 国产区一区二久久| 国产成+人综合+亚洲专区| 91成人精品电影| 叶爱在线成人免费视频播放| 国产视频内射| 美国免费a级毛片| 久久久久亚洲av毛片大全| 1024香蕉在线观看| 国产精品久久久久久人妻精品电影| 国产亚洲欧美98| 亚洲一区二区三区不卡视频| 在线视频色国产色| 亚洲av熟女| 欧美国产日韩亚洲一区| 搞女人的毛片| 亚洲专区中文字幕在线| 国产精品爽爽va在线观看网站 | 亚洲精品国产一区二区精华液| 一级毛片高清免费大全| 亚洲中文字幕一区二区三区有码在线看 | 久久精品91无色码中文字幕| 一进一出抽搐动态| 级片在线观看| 欧美 亚洲 国产 日韩一| 欧美一区二区精品小视频在线| 久久精品国产亚洲av香蕉五月| 听说在线观看完整版免费高清| 成人手机av| 国产激情久久老熟女| 欧美日本视频| 亚洲成人久久性| 少妇被粗大的猛进出69影院| 精品熟女少妇八av免费久了| 日本在线视频免费播放| 精品少妇一区二区三区视频日本电影| 无遮挡黄片免费观看| 国产一区二区激情短视频| 丝袜人妻中文字幕| 成熟少妇高潮喷水视频| 精品国产国语对白av| 色哟哟哟哟哟哟| 两个人免费观看高清视频| 香蕉丝袜av| 午夜激情福利司机影院| 91大片在线观看| 激情在线观看视频在线高清| 婷婷精品国产亚洲av在线| 亚洲精品美女久久av网站| 久久精品成人免费网站| 亚洲第一欧美日韩一区二区三区| 黄色毛片三级朝国网站| 大型黄色视频在线免费观看| 成人亚洲精品一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区av网在线观看| 黄色 视频免费看| 国产一区二区激情短视频| 中文字幕精品亚洲无线码一区 | 亚洲中文字幕一区二区三区有码在线看 | 亚洲色图 男人天堂 中文字幕| 99精品在免费线老司机午夜| 日本免费一区二区三区高清不卡| 日本精品一区二区三区蜜桃| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3 | 欧美精品亚洲一区二区| 人人妻,人人澡人人爽秒播| 国产不卡一卡二| 久久精品夜夜夜夜夜久久蜜豆 | 女警被强在线播放| 午夜激情福利司机影院| 亚洲avbb在线观看| 日日夜夜操网爽| 亚洲一区中文字幕在线| 悠悠久久av| 午夜免费观看网址| 国产黄色小视频在线观看| 嫁个100分男人电影在线观看| 中文字幕人妻熟女乱码| 日韩有码中文字幕| 99久久国产精品久久久| 亚洲欧美激情综合另类| 国内精品久久久久久久电影| 99久久精品国产亚洲精品| 国产区一区二久久| 可以免费在线观看a视频的电影网站| 久久久久国产一级毛片高清牌| 国产v大片淫在线免费观看| 午夜福利成人在线免费观看| 欧美中文综合在线视频| 欧美国产日韩亚洲一区| 观看免费一级毛片| 午夜激情av网站| av欧美777| 国产成人精品无人区| 人妻久久中文字幕网| 99国产精品一区二区三区| 午夜福利18| 日本免费a在线| 中出人妻视频一区二区| 99国产精品一区二区三区| 黄色女人牲交| 一级片免费观看大全| 狂野欧美激情性xxxx| 久久香蕉精品热| 久久久国产精品麻豆| 国产精品永久免费网站| 成人av一区二区三区在线看| 欧美日韩乱码在线| 九色国产91popny在线| av天堂在线播放| 久久精品91无色码中文字幕| 日日摸夜夜添夜夜添小说| 国产又爽黄色视频| 青草久久国产| 亚洲一区二区三区不卡视频| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 天堂影院成人在线观看| 一区二区三区高清视频在线| 精品国产乱码久久久久久男人| 法律面前人人平等表现在哪些方面| 国产午夜福利久久久久久| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 不卡av一区二区三区| 久久婷婷人人爽人人干人人爱| 久9热在线精品视频| 热re99久久国产66热| 人人妻人人澡欧美一区二区| 日本免费a在线| 欧美久久黑人一区二区| av天堂在线播放| 亚洲精品久久国产高清桃花| 亚洲中文日韩欧美视频| 草草在线视频免费看| 国产色视频综合| 欧美最黄视频在线播放免费| 国产成人啪精品午夜网站| 大型av网站在线播放| 午夜老司机福利片| 精品久久久久久久久久免费视频| 脱女人内裤的视频| 黄色成人免费大全| 国产在线观看jvid| 久久久久久久午夜电影| a在线观看视频网站| 午夜福利成人在线免费观看| 精品卡一卡二卡四卡免费| 不卡一级毛片| 在线永久观看黄色视频| 99久久无色码亚洲精品果冻| 首页视频小说图片口味搜索| 韩国av一区二区三区四区| 男人舔奶头视频| 久久香蕉国产精品| 久久亚洲真实| 日本在线视频免费播放| 色哟哟哟哟哟哟| 18禁美女被吸乳视频| 高潮久久久久久久久久久不卡| 97人妻精品一区二区三区麻豆 | 一进一出抽搐动态| 成人国产综合亚洲| 在线免费观看的www视频| 99久久国产精品久久久| 亚洲午夜精品一区,二区,三区| 19禁男女啪啪无遮挡网站| 国产亚洲欧美精品永久| 黄网站色视频无遮挡免费观看| 视频区欧美日本亚洲| 美女扒开内裤让男人捅视频| 久久久久久九九精品二区国产 | 亚洲狠狠婷婷综合久久图片| 久久伊人香网站| 日本免费一区二区三区高清不卡| 久久中文字幕一级| 国产激情久久老熟女| 欧美三级亚洲精品| 视频区欧美日本亚洲| 亚洲第一欧美日韩一区二区三区| 婷婷精品国产亚洲av| 日本精品一区二区三区蜜桃| 午夜激情福利司机影院| 日韩精品中文字幕看吧| 18禁国产床啪视频网站| bbb黄色大片| 亚洲三区欧美一区| 国产高清有码在线观看视频 | av视频在线观看入口| 在线播放国产精品三级| 一进一出抽搐gif免费好疼| 香蕉丝袜av| 天堂影院成人在线观看| 大型黄色视频在线免费观看| 亚洲av第一区精品v没综合| 观看免费一级毛片| 日韩精品青青久久久久久| 国产成人精品无人区| 夜夜躁狠狠躁天天躁| 国产亚洲精品第一综合不卡| 精品久久久久久久末码| 丁香六月欧美| 日韩 欧美 亚洲 中文字幕| 在线观看免费午夜福利视频| 国语自产精品视频在线第100页| 国产一区二区在线av高清观看| 亚洲avbb在线观看| 国产亚洲精品第一综合不卡| 久久香蕉激情| 一边摸一边做爽爽视频免费| 男女之事视频高清在线观看| 亚洲av中文字字幕乱码综合 | 99久久综合精品五月天人人| 18禁观看日本| 一本久久中文字幕| 亚洲国产精品合色在线| or卡值多少钱| 成人一区二区视频在线观看| 午夜免费观看网址| 欧美成人免费av一区二区三区| 日韩欧美国产一区二区入口| 啪啪无遮挡十八禁网站| 国产亚洲精品综合一区在线观看 | 国产精品香港三级国产av潘金莲| 久久久精品国产亚洲av高清涩受| 欧美三级亚洲精品| 啪啪无遮挡十八禁网站| 一本综合久久免费| 久久香蕉国产精品| 伦理电影免费视频| 久久人人精品亚洲av| 人人妻人人澡人人看| 精品国产超薄肉色丝袜足j| 99riav亚洲国产免费| 黄色成人免费大全| 99国产综合亚洲精品| 日本一本二区三区精品| 一区福利在线观看| 成年女人毛片免费观看观看9| 97碰自拍视频| 一级作爱视频免费观看| 日韩中文字幕欧美一区二区| 一级片免费观看大全| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 欧美最黄视频在线播放免费| 久久精品夜夜夜夜夜久久蜜豆 | 国产不卡一卡二| 午夜老司机福利片| 黄网站色视频无遮挡免费观看| 18禁美女被吸乳视频| 麻豆一二三区av精品| 香蕉av资源在线| 久久精品国产99精品国产亚洲性色| 丝袜在线中文字幕| 亚洲精华国产精华精| 国产伦在线观看视频一区| www.999成人在线观看| 日韩欧美在线二视频| 国产野战对白在线观看| 好男人在线观看高清免费视频 | 亚洲一区高清亚洲精品| 亚洲三区欧美一区| av视频在线观看入口| 精品国产国语对白av| 国产av不卡久久| 在线播放国产精品三级| 亚洲国产精品999在线| 中文字幕精品亚洲无线码一区 | 精华霜和精华液先用哪个| 午夜老司机福利片| 亚洲中文字幕一区二区三区有码在线看 | 色播亚洲综合网| 欧美绝顶高潮抽搐喷水| 亚洲人成77777在线视频| 人人澡人人妻人| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 日韩 欧美 亚洲 中文字幕| 可以免费在线观看a视频的电影网站| 女人高潮潮喷娇喘18禁视频| 很黄的视频免费| 好男人电影高清在线观看| 亚洲人成电影免费在线| 午夜视频精品福利| 久久久国产成人免费| 老司机靠b影院| 香蕉国产在线看| 国产一级毛片七仙女欲春2 | 女性被躁到高潮视频| 国产黄色小视频在线观看| 色播在线永久视频| 亚洲中文av在线| 99国产精品99久久久久| 女性生殖器流出的白浆| 人人妻人人看人人澡| 亚洲性夜色夜夜综合| 国产伦在线观看视频一区| 丝袜美腿诱惑在线| 国产伦在线观看视频一区| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 人人妻人人看人人澡| 淫秽高清视频在线观看| 精品国产国语对白av| 黄色视频,在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 久久久国产精品麻豆| 亚洲精品av麻豆狂野| av有码第一页| 在线观看舔阴道视频| 亚洲aⅴ乱码一区二区在线播放 | 成人一区二区视频在线观看| www.999成人在线观看| 久久青草综合色| av有码第一页| 国产亚洲精品综合一区在线观看 | 宅男免费午夜| 国产高清videossex| 丰满的人妻完整版|