• <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)計量的探究
    少妇被粗大猛烈的视频| 夫妻午夜视频| 美女内射精品一级片tv| 久久久久久久国产电影| 久久av网站| 两个人的视频大全免费| 久久久久久伊人网av| 国精品久久久久久国模美| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| 日本免费在线观看一区| 亚洲av成人精品一二三区| 视频中文字幕在线观看| 亚洲av综合色区一区| 午夜免费男女啪啪视频观看| 午夜日本视频在线| 男女国产视频网站| 少妇熟女欧美另类| 联通29元200g的流量卡| 国产女主播在线喷水免费视频网站| 一级毛片 在线播放| 91精品伊人久久大香线蕉| 免费看日本二区| 中国美白少妇内射xxxbb| 久久久久久九九精品二区国产| 777米奇影视久久| 日韩欧美精品免费久久| 亚洲人成网站高清观看| 日本黄大片高清| 成人国产av品久久久| 日本黄大片高清| 日日撸夜夜添| 热99国产精品久久久久久7| 国产乱人偷精品视频| 少妇高潮的动态图| 精品一区二区免费观看| 久久久久久人妻| 99re6热这里在线精品视频| 日日啪夜夜撸| 少妇裸体淫交视频免费看高清| a级毛色黄片| 久久久欧美国产精品| 国产永久视频网站| 免费av中文字幕在线| 久久久久网色| 联通29元200g的流量卡| 久久精品人妻少妇| 久久综合国产亚洲精品| 国产精品伦人一区二区| 欧美激情极品国产一区二区三区 | 亚洲国产精品国产精品| 久久国产精品大桥未久av | 亚洲成人手机| 又爽又黄a免费视频| 亚洲成色77777| 国产精品av视频在线免费观看| 观看美女的网站| 热re99久久精品国产66热6| 国产成人午夜福利电影在线观看| 亚洲精品视频女| 伦理电影免费视频| 美女视频免费永久观看网站| 伦理电影大哥的女人| 日韩伦理黄色片| 国产av国产精品国产| 美女视频免费永久观看网站| 久久国产精品大桥未久av | 国产精品久久久久成人av| 国产大屁股一区二区在线视频| av福利片在线观看| 十分钟在线观看高清视频www | 少妇人妻久久综合中文| 搡老乐熟女国产| 少妇人妻精品综合一区二区| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 男女边吃奶边做爰视频| 超碰97精品在线观看| 国产伦理片在线播放av一区| 国产探花极品一区二区| 午夜福利影视在线免费观看| 国产无遮挡羞羞视频在线观看| 狠狠精品人妻久久久久久综合| 免费人成在线观看视频色| 1000部很黄的大片| 成人漫画全彩无遮挡| 三级国产精品片| 国产毛片在线视频| 少妇人妻 视频| 美女福利国产在线 | 久久久久久九九精品二区国产| 精品一区在线观看国产| 精品99又大又爽又粗少妇毛片| 性高湖久久久久久久久免费观看| 黑人高潮一二区| 亚洲欧美中文字幕日韩二区| 亚洲精品一二三| 中文欧美无线码| 日本爱情动作片www.在线观看| 晚上一个人看的免费电影| 大香蕉久久网| 大码成人一级视频| 99热全是精品| 国产在线免费精品| 久久久久视频综合| 久久精品夜色国产| 欧美激情极品国产一区二区三区 | 国产午夜精品一二区理论片| 亚洲精品日韩av片在线观看| 久久韩国三级中文字幕| 国产成人精品一,二区| 两个人的视频大全免费| 国产综合精华液| 女的被弄到高潮叫床怎么办| 欧美最新免费一区二区三区| 国产国拍精品亚洲av在线观看| 国产美女午夜福利| 不卡视频在线观看欧美| 美女国产视频在线观看| 一本一本综合久久| 97超碰精品成人国产| 国产真实伦视频高清在线观看| av在线播放精品| 亚洲图色成人| 欧美精品一区二区大全| 日本av手机在线免费观看| 久久久久久久国产电影| 亚洲欧美精品自产自拍| 国产精品欧美亚洲77777| 日本黄色片子视频| 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 免费在线观看成人毛片| .国产精品久久| 久久久久久伊人网av| 久久韩国三级中文字幕| 免费看日本二区| 伦精品一区二区三区| 国内揄拍国产精品人妻在线| 色婷婷久久久亚洲欧美| 伦精品一区二区三区| 亚洲av成人精品一二三区| 毛片女人毛片| 亚洲,欧美,日韩| 成人国产麻豆网| 亚洲精品国产av成人精品| 国产淫语在线视频| 日产精品乱码卡一卡2卡三| 青青草视频在线视频观看| 免费av不卡在线播放| 男女下面进入的视频免费午夜| 国精品久久久久久国模美| 熟妇人妻不卡中文字幕| 成人免费观看视频高清| 最新中文字幕久久久久| 亚洲精品国产色婷婷电影| 亚洲精品国产av成人精品| 极品教师在线视频| 国产又色又爽无遮挡免| 三级国产精品欧美在线观看| 久久久久久伊人网av| 国产免费视频播放在线视频| 国产精品一二三区在线看| 成年女人在线观看亚洲视频| 黄色配什么色好看| 日日啪夜夜撸| 精品久久国产蜜桃| 国产精品久久久久成人av| av在线app专区| 交换朋友夫妻互换小说| av国产精品久久久久影院| 最近最新中文字幕大全电影3| 亚洲欧美一区二区三区黑人 | 久久99热6这里只有精品| 网址你懂的国产日韩在线| 国产一区亚洲一区在线观看| 欧美极品一区二区三区四区| av网站免费在线观看视频| 深爱激情五月婷婷| 国产精品麻豆人妻色哟哟久久| 男女边摸边吃奶| 欧美zozozo另类| 亚洲国产毛片av蜜桃av| 亚洲伊人久久精品综合| 国产成人免费无遮挡视频| 国产一区有黄有色的免费视频| 男人舔奶头视频| 国产高清有码在线观看视频| 中文字幕久久专区| 亚洲熟女精品中文字幕| 99久久精品国产国产毛片| 18+在线观看网站| 最新中文字幕久久久久| 人妻少妇偷人精品九色| 黄色日韩在线| 久久影院123| 免费黄网站久久成人精品| 午夜福利网站1000一区二区三区| 久久久久视频综合| 中文字幕精品免费在线观看视频 | 边亲边吃奶的免费视频| 777米奇影视久久| 国产中年淑女户外野战色| 国产精品99久久久久久久久| 天天躁夜夜躁狠狠久久av| 亚洲精品色激情综合| 王馨瑶露胸无遮挡在线观看| 嘟嘟电影网在线观看| 国产午夜精品一二区理论片| 中文天堂在线官网| 亚洲怡红院男人天堂| 亚洲国产最新在线播放| 老师上课跳d突然被开到最大视频| 日韩一本色道免费dvd| 在线观看免费高清a一片| 一级毛片aaaaaa免费看小| 少妇丰满av| 亚洲国产精品成人久久小说| 岛国毛片在线播放| 久久久久人妻精品一区果冻| 一级爰片在线观看| 久久久久精品久久久久真实原创| 熟女人妻精品中文字幕| 大陆偷拍与自拍| 噜噜噜噜噜久久久久久91| av黄色大香蕉| 亚洲成人中文字幕在线播放| 97精品久久久久久久久久精品| 看十八女毛片水多多多| 我要看日韩黄色一级片| 国产日韩欧美亚洲二区| 日产精品乱码卡一卡2卡三| 美女内射精品一级片tv| 中文资源天堂在线| 一本色道久久久久久精品综合| 亚洲人与动物交配视频| 欧美区成人在线视频| 亚洲av.av天堂| 亚洲aⅴ乱码一区二区在线播放| 99久久中文字幕三级久久日本| 春色校园在线视频观看| 极品少妇高潮喷水抽搐| 亚洲成人av在线免费| 99久久中文字幕三级久久日本| 国产免费一区二区三区四区乱码| 亚洲精品一区蜜桃| 天堂中文最新版在线下载| 九色成人免费人妻av| 亚洲国产色片| 午夜免费男女啪啪视频观看| 成人毛片a级毛片在线播放| 七月丁香在线播放| 美女福利国产在线 | 美女高潮的动态| 欧美日本视频| 婷婷色麻豆天堂久久| 国产国拍精品亚洲av在线观看| 免费观看的影片在线观看| 国产乱人偷精品视频| 看十八女毛片水多多多| 亚洲怡红院男人天堂| 97热精品久久久久久| 国产大屁股一区二区在线视频| 人体艺术视频欧美日本| 一级毛片久久久久久久久女| 国产av一区二区精品久久 | 久热久热在线精品观看| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区国产| 国产高潮美女av| 日韩欧美精品免费久久| .国产精品久久| 色视频在线一区二区三区| 久久婷婷青草| 国产亚洲91精品色在线| 最近最新中文字幕大全电影3| 日韩中字成人| 国产v大片淫在线免费观看| 91午夜精品亚洲一区二区三区| 日本一二三区视频观看| 九九在线视频观看精品| 国产成人a∨麻豆精品| 亚洲精品一二三| 国产成人精品婷婷| 日本欧美视频一区| 丰满人妻一区二区三区视频av| 美女高潮的动态| 国产永久视频网站| 丝袜脚勾引网站| 色网站视频免费| 日本av手机在线免费观看| 少妇的逼水好多| 国产一区亚洲一区在线观看| 高清在线视频一区二区三区| 亚洲av成人精品一二三区| 免费看av在线观看网站| 蜜臀久久99精品久久宅男| 老司机影院成人| 亚洲成人av在线免费| 女人十人毛片免费观看3o分钟| 中文资源天堂在线| 青春草亚洲视频在线观看| 日韩电影二区| 免费久久久久久久精品成人欧美视频 | 精品亚洲成a人片在线观看 | 国产精品国产av在线观看| 久久久久视频综合| 国产一区有黄有色的免费视频| 搡女人真爽免费视频火全软件| 国产精品伦人一区二区| 亚洲天堂av无毛| 国产视频内射| av播播在线观看一区| 水蜜桃什么品种好| 纵有疾风起免费观看全集完整版| 精品一区二区三区视频在线| 日韩欧美 国产精品| 日韩免费高清中文字幕av| 亚洲国产精品999| 国产伦精品一区二区三区视频9| 在线观看国产h片| 欧美日韩视频高清一区二区三区二| 国产无遮挡羞羞视频在线观看| 色网站视频免费| 久久久久人妻精品一区果冻| 乱码一卡2卡4卡精品| 国产精品国产三级专区第一集| 欧美成人a在线观看| 精品久久久噜噜| 少妇丰满av| 建设人人有责人人尽责人人享有的 | 不卡视频在线观看欧美| 又粗又硬又长又爽又黄的视频| 亚洲精品456在线播放app| 日韩欧美 国产精品| 免费观看无遮挡的男女| 最近最新中文字幕大全电影3| av播播在线观看一区| 亚洲国产精品成人久久小说| av福利片在线观看| 最近中文字幕高清免费大全6| 日韩一区二区视频免费看| 少妇人妻 视频| 国产69精品久久久久777片| 精品熟女少妇av免费看| 丰满迷人的少妇在线观看| 精品亚洲乱码少妇综合久久| 国产永久视频网站| 久久6这里有精品| 亚洲人成网站在线观看播放| 我的女老师完整版在线观看| 国产片特级美女逼逼视频| 亚洲久久久国产精品| 亚州av有码| 少妇熟女欧美另类| 免费av不卡在线播放| 国产亚洲午夜精品一区二区久久| 免费黄色在线免费观看| 伦理电影大哥的女人| 777米奇影视久久| 色视频www国产| 欧美三级亚洲精品| 青青草视频在线视频观看| 最近最新中文字幕大全电影3| 蜜臀久久99精品久久宅男| 亚洲高清免费不卡视频| 亚洲精品国产av蜜桃| 身体一侧抽搐| 国产黄片视频在线免费观看| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 简卡轻食公司| 午夜日本视频在线| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 亚洲美女搞黄在线观看| 亚洲av男天堂| 国产免费又黄又爽又色| 夫妻午夜视频| 中文资源天堂在线| 啦啦啦中文免费视频观看日本| 内射极品少妇av片p| 嫩草影院新地址| 国产探花极品一区二区| 国产精品久久久久久av不卡| 高清午夜精品一区二区三区| 亚洲四区av| 国产精品一区二区性色av| 精品一区二区三区视频在线| 亚洲av国产av综合av卡| 色视频www国产| 国产 一区精品| 久久久久久久大尺度免费视频| 99热国产这里只有精品6| 内地一区二区视频在线| 亚洲伊人久久精品综合| 日本欧美视频一区| 男女国产视频网站| 久久精品国产亚洲av涩爱| 日韩成人av中文字幕在线观看| 日韩,欧美,国产一区二区三区| 国产高清有码在线观看视频| 成人亚洲欧美一区二区av| 能在线免费看毛片的网站| 亚洲电影在线观看av| 精品一区二区三区视频在线| 汤姆久久久久久久影院中文字幕| 91精品国产九色| 久久精品久久久久久噜噜老黄| 国产爽快片一区二区三区| 日本与韩国留学比较| 色婷婷久久久亚洲欧美| av在线老鸭窝| 国产成人freesex在线| 欧美一区二区亚洲| 18+在线观看网站| 99久久精品热视频| 国产精品久久久久久久久免| 亚洲精品,欧美精品| 在线免费观看不下载黄p国产| 国产精品久久久久久久电影| 久久97久久精品| 国产毛片在线视频| 99视频精品全部免费 在线| 一级毛片电影观看| 国产一区有黄有色的免费视频| 日本wwww免费看| 99热6这里只有精品| 大话2 男鬼变身卡| 嫩草影院入口| 久久久欧美国产精品| 日韩中字成人| 天堂俺去俺来也www色官网| 成人毛片a级毛片在线播放| 亚洲va在线va天堂va国产| 岛国毛片在线播放| 亚洲国产精品一区三区| 日韩一区二区视频免费看| 爱豆传媒免费全集在线观看| 国产久久久一区二区三区| 日本与韩国留学比较| 国产成人aa在线观看| 亚洲欧美精品自产自拍| 日韩欧美精品免费久久| 三级经典国产精品| 日本爱情动作片www.在线观看| 欧美高清成人免费视频www| 高清av免费在线| 又黄又爽又刺激的免费视频.| 一区二区av电影网| 纯流量卡能插随身wifi吗| 少妇人妻一区二区三区视频| 亚洲精品成人av观看孕妇| 午夜福利在线观看免费完整高清在| 天天躁日日操中文字幕| 国产伦在线观看视频一区| 91精品国产国语对白视频| 九九在线视频观看精品| 97精品久久久久久久久久精品| 午夜福利在线在线| 尾随美女入室| 777米奇影视久久| 99久久人妻综合| 精品久久久久久电影网| 久久国产精品男人的天堂亚洲 | 国产乱人偷精品视频| 大码成人一级视频| 亚洲国产精品999| 国产精品久久久久久av不卡| 亚洲精品aⅴ在线观看| 99热全是精品| 日韩不卡一区二区三区视频在线| 国产成人a区在线观看| 国产视频首页在线观看| 成人免费观看视频高清| 精品一区二区三区视频在线| 国产 一区精品| 中文字幕人妻熟人妻熟丝袜美| 色哟哟·www| 久久久久久久久久成人| 亚洲综合精品二区| 亚洲精品乱码久久久久久按摩| 欧美日韩亚洲高清精品| 97超视频在线观看视频| 国产精品人妻久久久久久| 免费观看无遮挡的男女| 身体一侧抽搐| 国产一级毛片在线| 99久久中文字幕三级久久日本| 赤兔流量卡办理| 91在线精品国自产拍蜜月| 日本免费在线观看一区| 久久99热6这里只有精品| 99久久精品热视频| 午夜视频国产福利| 美女福利国产在线 | 一级毛片我不卡| 91午夜精品亚洲一区二区三区| 亚洲国产成人一精品久久久| 久久人人爽人人爽人人片va| 麻豆成人av视频| 最后的刺客免费高清国语| 亚洲激情五月婷婷啪啪| 内地一区二区视频在线| 多毛熟女@视频| 国产成人精品福利久久| 国产片特级美女逼逼视频| 久久国产乱子免费精品| 中文字幕人妻熟人妻熟丝袜美| 伦理电影免费视频| 网址你懂的国产日韩在线| 国产人妻一区二区三区在| 观看av在线不卡| 国产免费又黄又爽又色| 日本爱情动作片www.在线观看| 欧美三级亚洲精品| 美女高潮的动态| 中国三级夫妇交换| 人人妻人人爽人人添夜夜欢视频 | 麻豆乱淫一区二区| 久久久久久伊人网av| 18+在线观看网站| 国产免费一区二区三区四区乱码| 亚洲一区二区三区欧美精品| 在线免费十八禁| 国产黄频视频在线观看| 极品少妇高潮喷水抽搐| 久久影院123| 大陆偷拍与自拍| 免费在线观看成人毛片| 国产成人免费观看mmmm| 亚洲国产精品999| 日韩免费高清中文字幕av| 只有这里有精品99| 久久久a久久爽久久v久久| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 国产精品一二三区在线看| 日韩大片免费观看网站| 国产精品99久久久久久久久| 亚洲av免费高清在线观看| 日本与韩国留学比较| tube8黄色片| 欧美性感艳星| 夜夜看夜夜爽夜夜摸| 免费看av在线观看网站| 亚洲国产欧美在线一区| 免费大片黄手机在线观看| 五月天丁香电影| 看非洲黑人一级黄片| h视频一区二区三区| 午夜福利视频精品| 午夜日本视频在线| 乱码一卡2卡4卡精品| 久久综合国产亚洲精品| 大话2 男鬼变身卡| 国产欧美日韩一区二区三区在线 | 国产精品一及| 欧美人与善性xxx| 毛片女人毛片| 国产一区二区在线观看日韩| 久久久精品免费免费高清| 美女内射精品一级片tv| 免费人妻精品一区二区三区视频| 在线观看人妻少妇| 一边亲一边摸免费视频| 国产成人午夜福利电影在线观看| 日韩大片免费观看网站| 亚洲,一卡二卡三卡| 不卡视频在线观看欧美| 91久久精品国产一区二区成人| 大陆偷拍与自拍| 婷婷色综合www| 日韩av在线免费看完整版不卡| 韩国高清视频一区二区三区| 日韩伦理黄色片| 国产高清三级在线| 少妇人妻久久综合中文| 美女中出高潮动态图| 久久久久人妻精品一区果冻| 看免费成人av毛片| 亚洲欧美清纯卡通| 美女主播在线视频| 国产在线男女| 一级毛片我不卡| 哪个播放器可以免费观看大片| 97在线人人人人妻| 日韩av免费高清视频| 中文资源天堂在线| 日本爱情动作片www.在线观看| 国产白丝娇喘喷水9色精品| 美女内射精品一级片tv| 欧美三级亚洲精品| 亚洲精品第二区| 国产一区二区在线观看日韩| 在线观看人妻少妇| 欧美另类一区| 亚洲国产日韩一区二区| 18禁裸乳无遮挡免费网站照片| 午夜福利网站1000一区二区三区| 韩国av在线不卡| 日韩中字成人| 在线观看一区二区三区激情| 多毛熟女@视频| 久久鲁丝午夜福利片| 成人无遮挡网站| 久久久亚洲精品成人影院| 如何舔出高潮| 久久久久视频综合| 亚洲精品一二三| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人 | 国产成人精品婷婷| 国产伦理片在线播放av一区| 成人漫画全彩无遮挡|