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

    施肥措施對復(fù)墾土壤團(tuán)聚體碳氮含量和作物產(chǎn)量的影響

    2020-12-03 02:18:34曹寒冰謝鈞宇強(qiáng)久次仁洪堅平荊耀棟孟會生
    農(nóng)業(yè)工程學(xué)報 2020年18期

    曹寒冰,謝鈞宇,2,強(qiáng)久次仁,郭 璐,洪堅平,荊耀棟,孟會生

    施肥措施對復(fù)墾土壤團(tuán)聚體碳氮含量和作物產(chǎn)量的影響

    曹寒冰1,謝鈞宇1,2,強(qiáng)久次仁3,郭 璐4,洪堅平1,荊耀棟1,孟會生1※

    (1. 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院,太谷 030801;2. 山西省土壤環(huán)境與養(yǎng)分資源重點實驗室,太原 030031;3.西藏山南地區(qū)錯那縣農(nóng)牧局,山南 856700;4. 山東省濱州畜牧獸醫(yī)研究院,濱州 256600)

    研究復(fù)墾后不同施肥措施下有機(jī)碳(OC)和全氮(TN)在水穩(wěn)性團(tuán)聚體及粉黏粒組分中的分布特征,以期深入理解不同施肥措施下土壤有機(jī)碳的固持機(jī)制。以生土和連續(xù)6年不同施肥措施的復(fù)墾土壤為研究對象,采集0~20 cm耕層土壤樣品,利用濕篩法進(jìn)行土壤粒徑分組,分析大粒徑大團(tuán)聚體(> 2 mm)、小粒徑大團(tuán)聚體(>0.25~2 mm)、微團(tuán)聚體(0.053~0.25 mm)和粉黏粒組分(< 0.053 mm)中OC和TN含量,判斷各粒徑團(tuán)聚體及粉黏粒組分中有機(jī)碳儲量的驅(qū)動因素,探究團(tuán)聚體及粉黏粒組分中有機(jī)碳含量與作物產(chǎn)量之間的關(guān)系。試驗設(shè)不施肥(CK)、施氮磷鉀化肥(NPK)、單施有機(jī)肥(M)和有機(jī)無機(jī)肥配施(MNPK)4個處理。結(jié)果表明:1)整個試驗周期(2008—2013年),同CK相比, NPK、M以及MNPK處理均顯著提高了玉米籽粒產(chǎn)量,且以MNPK處理的效果最顯著,分別提高了79.49%、116.07%和113.85%。2)大團(tuán)聚體和微團(tuán)聚體中OC和TN含量相近,總體高于粉黏粒組分。同生土相比,CK、NPK和M處理均顯著提高了>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分中OC含量。MNPK處理顯著提高了各粒徑團(tuán)聚體和< 0.053 mm組分中OC含量。此外,CK、NPK、M和MNPK處理均顯著提高了各粒徑團(tuán)聚體中TN含量;大團(tuán)聚體和微團(tuán)聚體中有機(jī)碳和全氮富集系數(shù)相近,總體高于粉黏粒組分。M以及MNPK處理均顯著提高了各粒徑團(tuán)聚體及粉黏粒組分(除> 2 mm團(tuán)聚體外)的C/N比。3)土壤有機(jī)碳儲量受> 2 mm團(tuán)聚體的百分比數(shù)量驅(qū)動,其余各粒徑團(tuán)聚體及粉黏粒組分中有機(jī)碳含量對土壤有機(jī)碳儲量的貢獻(xiàn)較大。4)作物產(chǎn)量隨著各粒徑團(tuán)聚體及粉黏粒組分中有機(jī)碳含量的增加呈線性增加,說明該區(qū)域土壤有機(jī)碳庫仍未達(dá)到飽和,且以0.053~0.25 mm團(tuán)聚體中有機(jī)碳含量與作物產(chǎn)量關(guān)系最緊密。因此,在采煤塌陷區(qū)的低肥力復(fù)墾土壤上采取有機(jī)無機(jī)肥配施,對土壤質(zhì)量及農(nóng)業(yè)生態(tài)持續(xù)性具有積極的作用。

    復(fù)墾;土壤;水穩(wěn)性團(tuán)聚體;有機(jī)碳;全氮;作物產(chǎn)量

    0 引 言

    團(tuán)聚體作為土壤的重要組成部分,直接參與土壤的物理、化學(xué)和生物過程,其數(shù)量和質(zhì)量對于維持土壤肥力、改良土壤質(zhì)量和調(diào)節(jié)生態(tài)功能等方面發(fā)揮著積極作用[1-2]。有機(jī)碳作為土壤質(zhì)量的核心,是評價土壤肥力的一個關(guān)鍵指標(biāo),與團(tuán)聚體相互作用緊密,有機(jī)碳與礦物顆粒膠結(jié)形成穩(wěn)定的團(tuán)聚體[3-4],同時,團(tuán)聚體的物理保護(hù)機(jī)制將“禁錮”在其中的有機(jī)碳與微生物之間的空間“隔離”,從而減少有機(jī)碳的損失[5]。施肥、耕作、土地利用方式的改變以及種植制度等因素均會影響土壤團(tuán)聚體中有機(jī)碳、全氮的分布,且以施肥對它們的影響最顯著[6-7]。眾所周知,土壤是地球陸地生態(tài)系統(tǒng)中最大的有機(jī)碳庫,其庫容微小的變化都會對大氣中的碳濃度產(chǎn)生巨大的影響[8]。因此,通過合理的田間管理措施維持并增強(qiáng)土壤有機(jī)碳固存對于緩解溫室氣體排放、提升土壤質(zhì)量、促進(jìn)作物高產(chǎn)穩(wěn)產(chǎn)有重要意義。

    農(nóng)田土壤是受人類活動影響較頻繁的土壤類型之一[9],近年來國內(nèi)外就施肥對土壤團(tuán)聚體中有機(jī)碳和全氮含量的影響已有不少報道,但是研究結(jié)果不盡相同。蘇慧清等[10]利用棕壤36 a定位試驗發(fā)現(xiàn),有機(jī)碳和全氮主要分布在粉黏粒組分(< 0.053 mm)中,長期施化肥顯著提高了< 0.053 mm組分中有機(jī)碳和全氮含量,卻顯著降低了大團(tuán)聚體(0.5~1 mm)中有機(jī)碳和全氮含量,而無論是單施有機(jī)肥還是有機(jī)無機(jī)肥配施均顯著提高了0.5~1 mm團(tuán)聚體中碳氮含量。但是Xie等[11]通過中國西北地區(qū)塿土21 a定位肥料試驗研究發(fā)現(xiàn),有機(jī)碳和全氮主要固存在微團(tuán)聚體(< 0.25 mm)中,化肥及有機(jī)無機(jī)肥配施處理均顯著提高了大團(tuán)聚體(> 2、>0.25~2 mm)和微團(tuán)聚體(< 0.25 mm)中有機(jī)碳和全氮含量。此外,張藝等[12]報道,連續(xù)30 a在紅壤性水稻土上施有機(jī)肥,有機(jī)碳主要貯存在大粒徑大團(tuán)聚體(> 2 mm)中,而長期施化肥肥后顯著提高了小粒徑大團(tuán)聚體(>0.25~2 mm)中有機(jī)碳儲量。由此可見,施肥對不同土壤類型團(tuán)聚體中有機(jī)碳、全氮含量的影響并不一致,不同粒徑的團(tuán)聚體對協(xié)調(diào)、供應(yīng)土壤養(yǎng)分,增強(qiáng)土壤固碳能力等方面發(fā)揮的作用不盡相同[9]。造成研究結(jié)果之間的差異,可能與施肥歷史、土壤性質(zhì)、種植體系、氣候條件等因素有關(guān)。因此,有必要針對性地開展系統(tǒng)研究。

    山西省煤炭資源尤其豐富,是中國開采量比例最大的省份之一,據(jù)統(tǒng)計,每年原煤產(chǎn)量能達(dá)到全國的1/4以上。隨著時代的發(fā)展,煤炭資源開采在支撐社會經(jīng)濟(jì)發(fā)展的同時,也對土地資源和生態(tài)環(huán)境造成了持久而嚴(yán)重的負(fù)面影響[13]。據(jù)統(tǒng)計,山西因煤炭開采而損毀的耕地已超11萬hm2。采煤塌陷區(qū)復(fù)墾土壤是被井工開采破壞耕地后而進(jìn)行重構(gòu)的一種土壤類型,其層次和結(jié)構(gòu)均已發(fā)生了明顯的變化,由于重構(gòu)土壤就近取自于采煤塌陷區(qū)周邊黃土質(zhì)生土,其物理結(jié)構(gòu)、化學(xué)性質(zhì)以及微生物類群均處于初期構(gòu)建階段,與一般農(nóng)田土壤有顯著差異。因此,將破壞的土地復(fù)墾和重新利用,是緩解該區(qū)域人地矛盾的關(guān)鍵途徑,甚至可以應(yīng)對中國糧食安全面臨的嚴(yán)峻挑戰(zhàn)[14]。前人針對復(fù)墾土壤有機(jī)碳固存機(jī)制的研究主要在團(tuán)聚體中有機(jī)碳和全氮分布對施肥的響應(yīng)[15-16],而就驅(qū)動因素及其與作物產(chǎn)量的關(guān)系仍不清楚。因此,本研究以山西省東南部長治市襄垣縣復(fù)墾6 a的土壤為研究對象,分析不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳和全氮含量的變化特征,評估有機(jī)碳和全氮在其中的富集狀況,判斷土壤有機(jī)碳儲量變化的驅(qū)動因子,并探究作物產(chǎn)量與各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量的關(guān)系,以期為深入認(rèn)知復(fù)墾生態(tài)工程和培肥改良措施對土壤結(jié)構(gòu)與質(zhì)量的影響過程和機(jī)制提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗地概況

    試驗地位于山西省長治市襄垣縣王橋鎮(zhèn)(36°28′11.95″N,113°00′52.57″E,平均海拔980 m),屬潞安集團(tuán)五陽煤礦井田范圍,煤礦開采導(dǎo)致地面沉陷后地形呈馬鞍狀,最大落差為4~5 m,馬鞍狀峰距為150~180 m,水澆地變成了旱薄地,土地生產(chǎn)力嚴(yán)重下降。該區(qū)屬于低山黃土丘陵區(qū),屬于暖溫帶半濕潤大陸性季風(fēng)氣候,年平均氣溫9.5 ℃,年均降水量532.8 mm。供試土壤類型為石灰性褐土,黃土母質(zhì)。試驗開始前0~20 cm耕層土壤基本理化性質(zhì)如表1所示。圖1為研究區(qū)位置示意圖。

    表1 2008年耕層土壤(0~20 cm)基本理化性質(zhì)

    圖1 試驗區(qū)位圖

    1.2 試驗設(shè)計

    試驗開始于2008年,選擇塌陷年限相同且塌陷后地形基本一致的塌陷地,采用混推復(fù)墾的方式,將表土剝離50~100 cm,利用挖深墊淺、就地平整法對其進(jìn)行復(fù)墾,種植作物為春玉米(大豐30號),種植密度為60 000株/hm2,并于每年5月左右播種,10月左右收獲,且在收獲時將玉米秸稈全部粉碎、翻壓還田。從2013年開始進(jìn)行第6年復(fù)墾,不同處理下0~20 cm耕層土壤的基本理化性狀見表2。

    表2 復(fù)墾6a時不同處理下0~20 cm土層的基本性質(zhì)

    本研究包括4個處理,分別是不施肥(CK)、平衡施肥(NPK)、單施有機(jī)肥(M)、有機(jī)無機(jī)肥配施(MNPK)。供試肥料為尿素(含N 46%)、過磷酸鈣(含P2O512%)和氯化鉀(含K2O 60%),氮、磷、鉀肥和有機(jī)肥均作為基肥在玉米播種前一周一次性撒施,隨即翻耕入土。此外,另取未復(fù)墾生土(即0年)作為培肥效果對照。每個處理按照完全隨機(jī)區(qū)組設(shè)計,重復(fù)3次,小區(qū)面積為10 m × 5 m = 50 m2。除了CK處理外,各施肥處理依據(jù)等質(zhì)量養(yǎng)分施肥原則(表3)。試驗供試的有機(jī)肥為完全腐熟的雞糞,其中含有機(jī)碳:14.97%、氮:1.68%、磷:2.46%和鉀:1.35%。

    1.3 樣品采集與分析

    春玉米收獲前一天(2018年9月29日),用定制的環(huán)刀(高度為10 cm、直徑為10 cm)采集0~20 cm土層的原狀土樣,每個小區(qū)采集5個樣點,然后混合成1個樣品。帶回實驗室后,沿土壤結(jié)構(gòu)小心地掰成< 1 cm的土塊,剔除動植物殘體及石塊,過8 mm篩后,于通風(fēng)干燥處風(fēng)干,裝入塑封袋中備用。

    表3 不同農(nóng)田管理措施下肥料用量

    同時,春玉米收獲后,每個小區(qū)用土鉆(高度為20 cm、內(nèi)徑約為2.5 cm)采集0~20 cm土層土壤樣品6個并混合,裝入塑封袋中帶回實驗室,在陰涼通風(fēng)處自然風(fēng)干后剔除石塊、根系等雜物,研磨后過0.15 mm篩,保存在塑封袋中測定土壤基本理化性狀。

    采用Elliott等[17]提出的濕篩方法分離得到水穩(wěn)性團(tuán)聚體及粉黏粒組分,具體操作方法如下:首先準(zhǔn)備好2、0.25和0.053 mm孔徑的篩子,并按照大粒徑在上、小粒徑在下的順序擺放好,其次稱取50 g過8 mm篩的風(fēng)干土樣置于此套篩上,然后緩慢地將整套篩子放在裝滿2/3蒸餾水的桶中心,浸泡5 min,上下震蕩5 min,振幅為3 cm,最后將各個篩子中的土樣用蒸餾水小心地洗入已知重量的鋁盒中,同時向桶中加入0.2 mol/L CaCl2絮凝,再離心收集這部分土樣,上述所有收集的土樣于50 ℃烘至恒定的質(zhì)量,冷卻后稱質(zhì)量,即可得到> 2、>0.25~2、0.053~0.25 mm粒徑的團(tuán)聚體和< 0.053 mm的粉黏粒組分,磨細(xì)過0.15 mm篩裝入塑封袋中用于后續(xù)測定OC和TN含量。

    1.4 測定項目及方法

    土壤及各粒徑土壤團(tuán)聚體中有機(jī)碳含量采用重鉻酸鉀-容量法測定[18],全氮含量采用K2SO4-CuSO4-Se催化劑和濃硫酸消煮,半微量開氏法測定[18],速效磷采用碳酸氫鈉浸提經(jīng)連續(xù)流動分析儀測定[18],速效鉀采用醋酸銨浸提——火焰光度法測定[18]。

    1.5 計算方法

    土壤團(tuán)聚體及粉黏粒組分的有機(jī)碳富集系數(shù)和全氮富集系數(shù)計算公式如下:

    式中EE分別代表土壤團(tuán)聚體及粉黏粒組分的有機(jī)碳富集系數(shù)和全氮富集系數(shù);SOCa為土壤團(tuán)聚體及粉黏粒組分中有機(jī)碳含量,g/kg;TNa為土壤團(tuán)聚體及粉黏粒組分中全氮含量,g/kg。SOC為土壤有機(jī)碳含量,g/kg,TN為土壤全氮含量,g/kg。

    為了明確團(tuán)聚體對土壤總有機(jī)碳儲量的貢獻(xiàn)率,本文采用如下方法進(jìn)行計算:

    式中1和2分別表示引起土壤有機(jī)碳儲量變化的2個驅(qū)動因子,其中1是由團(tuán)聚體中有機(jī)碳含量變化所驅(qū)動的土壤有機(jī)碳儲量變化,2是由團(tuán)聚體的數(shù)量變化所驅(qū)動的土壤有機(jī)碳儲量變化。如果1>2,說明土壤有機(jī)碳儲量的變化是由該粒徑團(tuán)聚體或粉黏粒組分中有機(jī)碳含量的變化所引起的;反之,則說明是由該粒徑團(tuán)聚體或粉黏粒組分的數(shù)量所引起的。為生土的團(tuán)聚體分布比例,%;△為各處理下復(fù)墾后引起的團(tuán)聚體分布比例的變化值,%;為各處理下復(fù)墾后團(tuán)聚體中有機(jī)碳含量,g/kg;,△為各處理下復(fù)墾后引起的團(tuán)聚體中有機(jī)碳含量的變化差值,g/kg。

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

    試驗結(jié)果采用Microsoft Excel 2010軟件進(jìn)行數(shù)據(jù)整理,不同處理下作物產(chǎn)量、水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳和全氮含量、EE值以及C/N比采用LSD法進(jìn)行差異性檢驗。利用Origin 8.1 軟件作圖。

    2 結(jié)果與分析

    2.1 不同施肥措施下玉米籽粒產(chǎn)量

    整個試驗周期(2008—2013年),各處理中玉米籽粒平均產(chǎn)量的變幅范圍分別從725.70 kg/hm2到3 745.70 kg/hm2(CK),1 440.75 kg/hm2到6 218.80 kg/hm2(NPK),2 577.30 kg/hm2到7 944.70 kg/hm2(M),2 259.60 kg/hm2到8 156.65 kg/hm2(MNPK),平均產(chǎn)量的變化范圍從3 086.47 kg/hm2到6 668.87 kg/hm2。同對照(CK)相比,施化肥(NPK)、單施有機(jī)肥(M)以及有機(jī)無機(jī)肥配施(MNPK)均顯著提高了玉米籽粒產(chǎn)量,且以MNPK處理的效果最顯著(表4),分別比對照提高了79.49%、116.07%和113.85%。由此可見,連續(xù)多年有機(jī)無機(jī)肥配施有利于作物增產(chǎn)。

    表4 連續(xù)6 a復(fù)墾不同施肥措施下玉米籽粒產(chǎn)量

    注:不同小寫字母表示同一年份不同處理間達(dá)5%顯著水平。

    Note: Different lowercase letters for the same year mean significant differences between treatments at the 5% probability level.

    2.2 不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分的分布比例

    施肥措施及施肥年限下土壤團(tuán)聚體分布比例的兩因素方差分析結(jié)果顯示(表5),施肥措施并沒有影響各粒徑團(tuán)聚體及粉黏粒組分的分布比例,而復(fù)墾年限對各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的分布比例產(chǎn)生了顯著影響。隨著復(fù)墾時間的延長,各處理中0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分的分布比例均呈降低的趨勢。復(fù)墾6 a CK處理中> 2 mm團(tuán)聚體的分布比例是0年的3倍。

    表5 不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分的分布比例

    注:不同小寫字母表示同一粒徑年際間達(dá)5%顯著水平。為復(fù)墾年限,為施肥處理,×為復(fù)墾年限與施肥處理的交互效應(yīng)。**為< 0.01,ns為不顯著。

    Note: Different lowercase letters for the same aggregate size fraction mean significant differences between years at the 5% probability level.means reclamation years,means fertilization treatments,×means interactive effects of reclamation years and fertilization years. ** means< 0.01, ns means no significant differences.

    2.3 土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量

    不同施肥措施下各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量總體高于生土,且施有機(jī)肥處理(M和MNPK)的水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量增幅更大(圖 2a)。而就不同粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分而言,各粒徑團(tuán)聚體中有機(jī)碳含量相近,總體高于粉黏粒組分。

    同生土相比,CK和NPK處理均顯著提高了>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分中有機(jī)碳含量,增幅分別為87.07%~97.96%、66.19%~92.09%和92.80%~97.60%。M處理亦顯著提高了>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分中有機(jī)碳含量,增幅分別為95.24%、148.20%和112.00%。MNPK處理顯著提高了> 2、>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分中有機(jī)碳含量,分別為生土的2.52倍、2.61倍、2.35倍和1.93倍。

    綜上,MNPK處理對于提高大團(tuán)聚體(> 2和>0.25~2 mm)中有機(jī)碳含量效果最明顯,M處理對于提高微團(tuán)聚體(0.053~0.25 mm)和粉黏粒組分(< 0.053 mm)中有機(jī)碳效果最明顯。

    2.4 土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中全氮含量

    不同施肥措施下各粒徑水穩(wěn)性團(tuán)聚體中全氮含量均顯著高于生土,變化趨勢與有機(jī)碳一致(圖2 b)。但是所有處理較生土均對< 0.053 mm組分中全氮含量無顯著影響。同生土相比,CK和NPK處理均顯著提高了> 2、>0.25~2和0.053~0.25 mm團(tuán)聚體中全氮含量,增幅分別為41.83%~43.66%、46.78%~52.35%和29.95%~30.43%。M和MNPK處理亦均顯著提高了> 2、>0.25~2和0.053~0.25 mm團(tuán)聚體中全氮含量,增幅分別為52.81%~63.80%、50.43%~62.96%和39.61%~44.93%。

    綜上,MNPK處理對于提高大團(tuán)聚體(> 2 mm和>0.25~2 mm)中全氮含量效果最明顯,M處理對于提高微團(tuán)聚體(0.053~0.25 mm)中全氮含量效果最明顯。

    2.5 土壤團(tuán)聚體及粉黏粒組分有機(jī)碳和全氮富集系數(shù)

    長期不同施肥措施在一定程度上影響了各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的有機(jī)碳富集系數(shù)(E值)(表 6)??傮w看來,大團(tuán)聚體和微團(tuán)聚體的E值相近(除CK處理外),而粉黏粒組分的E值小于各粒徑團(tuán)聚體。同生土相比,CK處理顯著降低了> 2 mm團(tuán)聚體的E值,達(dá)46.94%。MNPK處理顯著降低了< 0.053 mm組分的E值,達(dá)23.08%。而其余施肥處理對各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的E值無顯著影響。

    注:不同小寫字母表示同一粒徑不同處理間達(dá)5%顯著水平;不同大寫字母表示同一處理不同粒徑間達(dá)5%顯著水平。

    長期不同施肥措施對> 2 mm和>0.25~2 mm團(tuán)聚體的全氮富集系數(shù)(E值)產(chǎn)生了顯著影響,而對0.053~0.25 mm團(tuán)聚體及< 0.053 mm組分的E值無顯著影響(表 6)。同生土相比,CK和NPK處理顯著提高了> 2 mm和>0.25~2 mm團(tuán)聚體的E值,增幅分別為17.00%~19.21%和19.70%~28.14%。M和MNPK處理亦顯著提高了>2 mm和>0.25~2 mm團(tuán)聚體的E值,增幅分別為13.33%~34.58%和11.68%~33.37%,且以MNPK處理的增幅最大。

    2.6 土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分的碳氮比

    不同施肥措施顯著影響了各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的C/N比(除> 2 mm團(tuán)聚體)(表7)。同生土相比,CK處理顯著提高了>0.25~2 mm團(tuán)聚體和< 0.053 mm組分的C/N比,分別提高了34.10%和79.66%。NPK處理顯著提高了0.053~0.25 mm團(tuán)聚體的C/N比,達(dá)46.11%。M和MNPK處理均顯著提高了>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分的C/N比,增幅分別為29.81%~60.81%、67.70%~71.32%和81.32%~97.45%。

    表6 不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分的有機(jī)碳富集系數(shù)和全氮富集系數(shù)

    注:不同小寫字母表示同一粒徑不同處理間達(dá)5%顯著水平;不同大寫字母表示同一處理不同粒徑間達(dá)5%顯著水平。

    Note: Different lowercase letters for the same aggregate size fraction mean significant differences among treatments at the 5% probability level. Different uppercase letters for the same treatment mean significant differences between aggregate size fractions at the 5% probability level.

    總之,M和MNPK處理均顯著提高了各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分(除> 2 mm團(tuán)聚體外)的C/N比。

    表7 不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分的碳氮比

    2.7 不同施肥措施下復(fù)墾土壤團(tuán)聚體及粉黏粒組分有機(jī)碳儲量的驅(qū)動因子

    由圖3可知,經(jīng)過不同的施肥措施管理農(nóng)田后,>2 mm團(tuán)聚體的2值顯著高于1值,而>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分的2值均顯著低于1值。其中,在> 2 mm團(tuán)聚體中,各處理下的2值均顯著高于生土,1值亦均顯著高于生土(除CK處理外)。>0.25~2 mm團(tuán)聚體中各處理下的1值和2值亦顯著高于生土。但是在0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分中,各處理下的1值顯著高于生土,而2值顯著低于生土。

    注:F1表示團(tuán)聚體中有機(jī)碳含量的變化;F2表示團(tuán)聚體數(shù)量的變化。

    2.8 不同施肥措施下土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量與作物產(chǎn)量的關(guān)系

    由圖4 a可知,隨著土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量的增加,玉米產(chǎn)量顯著增加。當(dāng)土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量每增加1 g/kg時,玉米產(chǎn)量增加2 517 kg/hm2。進(jìn)一步分析后發(fā)現(xiàn),0.053~0.25 mm團(tuán)聚體中有機(jī)碳含量與作物籽粒產(chǎn)量之間的關(guān)系最密切(圖4 d),當(dāng)該粒徑團(tuán)聚體中有機(jī)碳含量增加1 g/kg時,作物產(chǎn)量可以增加3 075 kg/hm2;其次是> 2 mm團(tuán)聚體與作物產(chǎn)量的關(guān)系較為密切(圖4 b),當(dāng)該粒徑團(tuán)聚體中有機(jī)碳含量增加1 g/kg時,作物產(chǎn)量可以增加2 245 kg/hm2。然而>0.25~2 mm團(tuán)聚體(圖4 c)和<0.053 mm組分(圖4 e)中有機(jī)碳含量對作物產(chǎn)量無顯著影響。

    圖4 土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量與復(fù)墾第6年作物產(chǎn)量的關(guān)系

    3 討 論

    3.1 不同施肥措施對復(fù)墾土壤水穩(wěn)性團(tuán)聚體中有機(jī)碳固存的影響

    土壤團(tuán)聚體作為土壤有機(jī)質(zhì)儲存的倉庫:表層土壤有機(jī)碳90%以上貯存在團(tuán)聚體中,團(tuán)聚體的團(tuán)聚作用對實現(xiàn)農(nóng)田土壤有機(jī)碳固存具有重要意義[5]。本研究結(jié)果顯示,不同施肥措施下復(fù)墾土壤各粒徑水穩(wěn)性團(tuán)聚體中有機(jī)碳和全氮含量總體均顯著高于生土,且隨著外源碳投入量的增加呈上升趨勢,特別是有機(jī)無機(jī)肥配施處理(MNPK)對大團(tuán)聚體(> 2 mm和>0.25~2 mm)中有機(jī)碳和全氮含量有明顯的提升效果(圖2)。而粉黏粒組分(< 0.053 mm)中有機(jī)碳和全氮含量維持在一個相對穩(wěn)定的水平,無明顯變化。此外,大團(tuán)聚體和微團(tuán)聚體中有機(jī)碳和全氮含量相近,但是總體高于< 0.053 mm組分的。這與多數(shù)研究結(jié)果類似[10,19]。他們也發(fā)現(xiàn),長期有機(jī)無機(jī)肥配施(MNPK)較不施肥(CK)均顯著提高了紅壤和棕壤上各粒徑水穩(wěn)性團(tuán)聚體中有機(jī)碳或全氮含量,尤其是> 0.25 mm團(tuán)聚體。不難理解,外源有機(jī)質(zhì)的輸入增加了土壤有機(jī)質(zhì)含量,有機(jī)質(zhì)中的膠結(jié)物質(zhì)與土壤顆粒形成團(tuán)聚體,繼而增加了團(tuán)聚體中有機(jī)碳含量[20]。雖然生土、CK和NPK處理下有機(jī)碳投入量均來源于作物殘茬,但是CK和NPK處理下作物產(chǎn)量顯著高于生土,且NPK處理對作物增產(chǎn)的效果優(yōu)于CK處理(表4),所以歸還到農(nóng)田中的碳投入量也多,而M以及MNPK處理下土壤有機(jī)碳投入量除了來源于作物的,還有相當(dāng)一部分是來源于本身—有機(jī)肥。故施有機(jī)肥的處理對團(tuán)聚體中有機(jī)碳及全氮含量的提升效果最優(yōu)。曾有研究報道大團(tuán)聚體中有機(jī)碳含量與土壤肥力呈正相關(guān)[21],因此MNPK處理是當(dāng)前復(fù)墾土壤土質(zhì)提升的最佳措施。

    不同施肥措施之間< 0.053 mm組分中有機(jī)碳和全氮含量無明顯變化,說明該組分已接近或達(dá)到了碳氮飽和水平。苑亞茹等[22]在黑土上也報道,同施化肥相比,連續(xù)11 a不同劑量的有機(jī)肥與化肥配施對<0.053 mm組分中有機(jī)碳含量無顯著影響。Kool等[23]在團(tuán)聚體層級發(fā)育模型的基礎(chǔ)上提出了土壤有機(jī)碳的等級飽和模型,即團(tuán)聚體中有機(jī)碳含量并沒有隨著外源碳投入量的增加呈線性趨勢,最小粒徑的團(tuán)聚體首先飽和,最大粒徑的團(tuán)聚體最后飽和,最終土壤有機(jī)碳庫達(dá)到飽和。各粒徑水穩(wěn)性團(tuán)聚體中有機(jī)碳和全氮含量之間無顯著差異,這與Li等[24]的研究結(jié)果一致,可能是因為團(tuán)聚體中有機(jī)碳和全氮的“濃縮效應(yīng)”[24]。

    前人通過計算元素富集系數(shù)(Element enrich factors,)來比較不同施肥措施下各粒徑團(tuán)聚體及粉黏粒組分中碳氮的分布狀況。當(dāng)> 1時,說明土壤組分富集碳和氮,反之,當(dāng)< 1時,說明土壤組分中這2種元素已損失[25-26]。本研究結(jié)果表明,生土和復(fù)墾土壤中各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的E值均小于1(表6),而E值接近或大于1。同生土相比,MNPK處理總體提高了大團(tuán)聚體的EE值,但降低了< 0.053 mm組分的EE值(表6)。說明該區(qū)域土壤中活性有機(jī)碳首先發(fā)生了礦化,且礦化程度較高,氮素在土壤中的累積程度高于有機(jī)碳。這一方面是因為玉米生長于高溫多雨的夏季,集中降水沖刷及強(qiáng)光暴曬均會加快土壤有機(jī)碳的分解;另一方面是因為生土或者復(fù)墾土壤上,自身土壤肥力的限制,再加上外源補(bǔ)充碳不足導(dǎo)致的。值得注意的是,氮素的大量累積一部分被作物吸收利用,一部分因揮發(fā)或淋溶而損失,還有一部分便殘留在土壤中,至地下水和土體的環(huán)境風(fēng)險增加[27],有關(guān)這方面的研究有待于后續(xù)開展。

    有機(jī)肥的輸入為土壤中微生物的活動提供了能量,導(dǎo)致了很強(qiáng)的正激發(fā)效應(yīng)[28],使得< 0.053 mm組分中的碳氮已損失(表6),此時粉黏粒組分中的有機(jī)碳和全氮已經(jīng)向團(tuán)聚體中轉(zhuǎn)移。可見,本研究所在區(qū)域的土壤上,團(tuán)聚體中有機(jī)碳和全氮受化學(xué)保護(hù)作用是有限的,而有機(jī)碳和全氮主要受團(tuán)聚體的物理保護(hù)作用。這與多數(shù)人的研究結(jié)果一致[29-31]。邵興芳[29]在黑土旱地的團(tuán)聚體氮礦化研究結(jié)果表明,< 0.053 mm組分中氮礦化量明顯高于各大小粒徑團(tuán)聚體的。楊洪波等[30]也發(fā)現(xiàn)了同樣的現(xiàn)象:< 0.053 mm組分比大團(tuán)聚體中含有更多的易礦化有機(jī)氮。Craswell等[31]將不同粒徑土壤團(tuán)聚體培養(yǎng)后發(fā)現(xiàn),隨著團(tuán)聚體粒徑的減小,其中的氮礦化率逐漸增大。所以,MNPK處理下< 0.053 mm組分中的碳氮易先發(fā)生礦化現(xiàn)象。

    3.2 不同施肥措施對復(fù)墾土壤水穩(wěn)性團(tuán)聚體中C/N比的影響

    土壤C/N比值是評價有機(jī)物分解程度的一個指標(biāo),比值越低,說明有機(jī)物的分解程度越高,可供微生物利用的碳源較少,微生物的活性降低,從而影響?zhàn)B分的有效性;反之則表示有機(jī)碳累積程度高[27]。本研究顯示,M和MNPK處理較生土均顯著提高了各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的C/N比(除> 2 mm團(tuán)聚體外)(表7),這與蘇慧清等[10]的研究結(jié)果一致。但是不同施肥措施下復(fù)墾土壤各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分的C/N比均低于一般農(nóng)田,這可能是因為礦區(qū)復(fù)墾土壤由于其物理結(jié)構(gòu)、化學(xué)性質(zhì)以及微生物類群均處于初期構(gòu)建階段,培肥時間還較短,與一般農(nóng)田土壤有所差異。因此,可通過增加土壤中的有機(jī)碳含量來提高土壤的C/N比,前人研究報道,減氮運籌下增施有機(jī)肥,能夠使有機(jī)碳增加9%~54%[32-33]。

    3.3 不同施肥措施下復(fù)墾土壤團(tuán)聚體中有機(jī)碳固存的驅(qū)動因素

    土壤各粒級團(tuán)聚體中有機(jī)碳含量是土壤有機(jī)物質(zhì)平衡和轉(zhuǎn)化速率的微觀表征,對于土壤肥力和土壤碳固存具有雙重意義[21]。本研究結(jié)果顯示,長期復(fù)墾后土壤有機(jī)碳儲量增加,是因為表層土壤> 2 mm團(tuán)聚體的數(shù)量增加,其次才依賴于該粒徑團(tuán)聚體中有機(jī)碳含量的增加(圖 2 a)。正如本研究中表5所示,隨著復(fù)墾時間的延長,各施肥處理下大團(tuán)聚體的數(shù)量呈增加趨勢。對于>0.25~2 、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分來講,它們中的有機(jī)碳含量對土壤有機(jī)碳儲量提升的作用大于它們的數(shù)量。這也潛在說明復(fù)墾過程中團(tuán)聚體中有機(jī)碳發(fā)生了再分配過程。王興等[21]的研究也表明,植被恢復(fù)可以加速有機(jī)碳從細(xì)顆粒組分向較大粒徑的組分中轉(zhuǎn)移,即有機(jī)碳的在分配過程,這種轉(zhuǎn)化主要是由于在植被恢復(fù)過程中加速了細(xì)顆粒組分易于和土壤大團(tuán)聚體結(jié)合,從而使有機(jī)碳受到保護(hù),不易被微生物所分解。

    3.4 不同施肥措施對復(fù)墾土壤上作物產(chǎn)量的影響

    作物產(chǎn)量是農(nóng)田土壤生產(chǎn)力的最直接體現(xiàn),是農(nóng)田土壤肥力變化的最有力指示物。本研究結(jié)果表明,同生土相比,不同施肥措施均顯著提高了作物產(chǎn)量,且以MNPK處理效果最顯著(表4)。與Xie等[34]在塿土上的研究結(jié)果一致,連續(xù)35 a有機(jī)無機(jī)肥配施顯著提高了小麥和玉米籽粒產(chǎn)量。有機(jī)無機(jī)肥配施能夠優(yōu)化土壤結(jié)構(gòu),補(bǔ)充和平衡土壤養(yǎng)分,提高土壤肥力,最終提升作物產(chǎn)量。

    本研究結(jié)果還顯示,作物產(chǎn)量隨著土壤水穩(wěn)性團(tuán)聚體及粉粘黏粒組分中有機(jī)碳含量的增加呈線性增加,且以0.053~0.25 mm組分中有機(jī)碳含量與作物產(chǎn)量之間的關(guān)系最密切,其次是> 2 mm團(tuán)聚體,如前所述,本研究所在區(qū)域的土壤上,團(tuán)聚體對有機(jī)碳的物理保護(hù)作用發(fā)揮著主要作用,提升了土壤肥力,進(jìn)而促進(jìn)了作物增產(chǎn)。研究報道土壤有機(jī)碳含量與作物產(chǎn)量之間呈線性相關(guān)關(guān)系[35],說明該區(qū)域土壤有機(jī)碳庫還未達(dá)到飽和水平,這主要是因為該土壤初始有機(jī)碳水平較低[36],因此,仍然可以繼續(xù)向該土壤上補(bǔ)充碳源。關(guān)于復(fù)墾褐土有機(jī)碳的飽和點,尚有待于進(jìn)一步深入研究。

    4 結(jié) 論

    1)整個試驗周期(2008—2013年),同CK相比, NPK、M以及MNPK處理均顯著提高了玉米籽粒產(chǎn)量,且以MNPK處理的效果最顯著,分別提高了79.49%、116.07%和113.85%。

    2)同生土相比,不同施肥措施均顯著提高了復(fù)墾土壤各粒徑水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳和全氮含量(除< 0.053 mm組分中全氮含量外),且以有機(jī)無機(jī)肥配施處理(MNPK)對提高大團(tuán)聚體(> 2 mm和>0.25~2 mm)中有機(jī)碳和全氮含量效果最明顯,單施有機(jī)肥處理(M)對提高微團(tuán)聚體(0.053~0.25 mm)中有機(jī)碳和全氮效果最明顯。各粒徑團(tuán)聚體中有機(jī)碳和全氮含量相近,總體高于粉黏粒組分(< 0.053 mm)。

    3)> 2 mm粒徑團(tuán)聚體對復(fù)墾土壤有機(jī)碳儲量提升的貢獻(xiàn)在于其數(shù)量的增加,而>0.25~2、0.053~0.25 mm團(tuán)聚體和< 0.053 mm組分對土壤有機(jī)碳儲量提升的貢獻(xiàn)在于其有機(jī)碳含量的增加。

    4)作物產(chǎn)量隨著復(fù)墾土壤水穩(wěn)性團(tuán)聚體及粉黏粒組分中有機(jī)碳含量的增加呈線性增加,說明土壤有機(jī)碳庫仍未飽和。當(dāng)土壤有機(jī)碳含量增加1 g/kg時,作物產(chǎn)量提高2 517 kg/hm2。且以0.053~0.25 mm團(tuán)聚體中有機(jī)碳含量與作物產(chǎn)量的關(guān)系最密切。

    [1]Haydu-houdeshell C, Graham R, Hendrix P, et al. Soil aggregate stability under chaparral species in southern California[J]. Geoderma, 2018, 310: 201-208.

    [2]張迪,姜佰文,梁世鵬,等. 草甸黑土團(tuán)聚體穩(wěn)定性對耕作與炭基肥施用的響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(14):125-132. Zhang Di, Jiang Baiwen, Liang Shipeng, et al. Responsive of aggregate stability of meadow black soil to different tillage practices and carbon-based fertilizers[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(14): 125-132. (in Chinese with English abstract)

    [3]李文軍,楊基峰,彭保發(fā),等. 施肥對洞庭湖平原水稻土團(tuán)聚體特征及其有機(jī)碳分布的影響[J]. 中國農(nóng)業(yè)科學(xué),2014,47(20):4007-4015. Li Wenjun, Yang Jifeng, Peng Baofa, et al. Effects of fertilization on aggregate characteristics and organic carbon distribution in a paddy soil in Dongting Lake Plain of China[J]. Scientia Agricultura Sinica, 2014, 47(20): 4007-4015. (in Chinese with English abstract)

    [4]劉振東,李貴春,周穎,等. 無機(jī)肥配施糞肥對華北褐土團(tuán)聚體分布及有機(jī)碳含量的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2013,32(11):2239-2245. Liu Zhendong, Li Guichun, Zhou Ying, et al. The effect of fertilizer management practices on distribution of aggregates and SOC[J]. Journal of Agro-Environment Science, 2013, 32(11): 2239-2245. (in Chinese with English abstract)

    [5]李昊昱,孟兆良,龐黨偉,等. 周年秸稈還田對農(nóng)田土壤固碳及冬小麥-夏玉米產(chǎn)量的影響[J]. 作物學(xué)報,2019,45(6):893-903. Li Haoyu, Meng Zhaoliang, Pang Dangwei, et al. Effect of annual straw return model on soil carbon sequestration and crop yields in winter wheat-summer maize rotation farmland[J]. Acta Agronomica Sinica, 2019, 45(6): 893-903. (in Chinese with English abstract)

    [6]Graf-rosenfellner M, Cierjacks A, Kleinschmit B, et al. Soil formation and its implications for stabilization of soil organic matter in the riparian zone[J]. Catena, 2016, 139: 9-18.

    [7]邵慧蕓,李紫玥,劉丹,等. 有機(jī)肥施用量對土壤有機(jī)碳組分和團(tuán)聚體穩(wěn)定性的影響[J]. 環(huán)境科學(xué),2019,40(10):4691-4699. Shao Huiyun, Li Ziyue, Liu Dan, et al. Effects of manure application rates on the soil carbon fractions and aggregate stability[J]. Environmental Science, 2019, 40(10): 4691-4699. (in Chinese with English abstract)

    [8]Stockmann U, Adams M A, Crawford J W, et al. The knowns, known unknowns and unknowns of sequestration of soil organic carbon[J]. Agriculture, Ecosystems & Environment, 2013, 164: 80-99.

    [9]劉杰,馬艷婷,王憲玲,等. 渭北旱塬土地利用方式對土壤團(tuán)聚體穩(wěn)定性及其有機(jī)碳的影響[J]. 環(huán)境科學(xué),2019,40(7):3361-3368. Liu Jie, Ma Yanting, Wang Xianling, et al. Impact of land use type on the stability and organic carbon content of soil aggregates in the Weibei dryland[J]. Environmental Science, 2019, 40(7): 3361-3368. (in Chinese with English abstract)

    [10]蘇慧清,韓曉日,楊勁峰,等. 長期施肥棕壤團(tuán)聚體分布及其碳氮含量變化[J]. 植物營養(yǎng)與肥料學(xué)報,2017,23(4):924-932. Su Huiqing, Han Xiaori, Yang Jinfeng, et al. Effect of long-term fertilization on distribution of aggregates and organic carbon and total nitrogen contents in a brown soil[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(4): 924-932. (in Chinese with English abstract)

    [11]Xie J, Hou M, Zhou Y, et al. Carbon sequestration and mineralization of aggregate-associated carbon in an intensively cultivated anthrosol in north China as affected by long term fertilization[J]. Geoderma, 2017, 296: 1-9.

    [12]張藝,戴齊,尹力初,等. 后續(xù)施肥措施改變對水稻土團(tuán)聚體有機(jī)碳分布及其周轉(zhuǎn)的影響[J]. 土壤,2017,49(5):969-976. Zhang Yi, Dai Qi, Yi Lichu, et al. Effects of following-up fertilization reforming on distribution and turnover of aggregate-associated organic carbon in paddy soils[J]. Soils, 2017, 49(5): 969-976. (in Chinese with English abstract)

    [13]李晉川,白中科,柴書杰,等. 平朔露天煤礦土地復(fù)墾與生態(tài)重建技術(shù)研究[J]. 科技導(dǎo)報,2009,27(17):30-34. Li Jinchuan, Bai Zhongke, Chai Shujie, et al. Study on technology of land reclamation and ecological rehabilitation of waste land in Pingshuo surface mine[J]. Science & Technology Review, 2009, 27(17): 30-34. (in Chinese with English abstract)

    [14]姚致遠(yuǎn). 旱地豆科綠肥提升土壤碳氮儲量及降低環(huán)境代價的潛力與機(jī)制[D]. 楊凌:西北農(nóng)林科技大學(xué),2019. Yao Zhiyuan. Potential and Mechanisms of Dryland Leguminous Green Manure for Improving Soil Carbon/Nitrogen Stocks and Reducing Environmental Costs[D]. Yangling: Northwest A&F University, 2019. (in Chinese with English abstract)

    [15]何冰,李廷亮,栗麗,等. 采煤塌陷區(qū)復(fù)墾土壤團(tuán)聚體碳氮分布對施肥的響應(yīng)[J]. 水土保持學(xué)報,2018,32(4):184-196. He Bing, Li Tingliang, Li Li, et al. Response of carbon and nitrogen distribution of reclaimed soil aggregates to fertilizers in coal mining subsidence area[J]. Journal of Soil and Water Conservation, 2018, 32(4): 184-196. (in Chinese with English abstract)

    [16]梁利寶,馮鵬艷,許劍敏. 施肥對采煤塌陷復(fù)墾土壤團(tuán)聚體組成及其碳、氮分布的影響[J]. 灌溉排水學(xué)報,2019,38(7):45-51. Liang Libao, Feng Pengyan, Xu Jianmin. Efficacy of fertilization in improving soil aggregation, carbon and nitrogen in soil reclaimed from subsided areas caused by coal mining[J]. Journal of Irrigation and Drainage, 2019, 38(7): 45-51. (in Chinese with English abstract)

    [17]Elliott E T. Aggregate structure and carbon, nitrogen, and phosphorus in native and cultivated soils[J]. Soil Science Society of America Journal, 1986, 50(3): 627-633.

    [18]鮑士旦. 土壤農(nóng)化分析[M]. 北京:中國農(nóng)業(yè)出版社, 2000.

    [19]邸佳穎,劉小粉,杜章留,等. 長期施肥對紅壤性水稻土團(tuán)聚體穩(wěn)定性及固碳特征的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報,2014,22(10):1129-1138. Di Jiaying, Liu Xiaofen, Du Zhangliu, et al. Influences of long-term organic and chemical fertilization on soil aggregation and associated organic carbon fractions in a red paddy soil[J]. Journal of Chinese Eco-Agriculture, 2014, 22(10): 1129-1138. (in Chinese with English abstract)

    [20]Sodhi G P S, Beri V, Benbi D K. Soil aggregation and distribution of carbon and nitrogen in different fractions under long-term application of compost in rice-wheat system[J]. Soil and Tillage Research, 2009, 103(2): 412-418.

    [21]王興,鐘澤坤,張欣怡,等. 長期撂荒恢復(fù)土壤團(tuán)聚體組成與有機(jī)碳分布關(guān)系[J]. 環(huán)境科學(xué),2020,41(5):2416-2424. Wang Xing, Zhong Zekun, Zhang Xinyi, et al. Relational between the composition of soil aggregates and the distribution of organic carbon under long-term abandoned restoration[J]. Environmental Science, 2019, 2020, 41(5): 2416-2424. (in Chinese with English abstract)

    [22]苑亞茹,鄒文秀,郝翔翔,等. 黑土團(tuán)聚體結(jié)合碳對不同有機(jī)肥施用量的響應(yīng)[J]. 生態(tài)學(xué)報,2019,39(9):3235-3242. Wan Yaru, Zou Wenxiu, Hao Xiangxiang, et al. Response of aggregate-associated carbon in mollisols to application rate of organic manure[J]. Acta Ecologica Sinica, 2019, 39(9): 3235-3242. (in Chinese with English abstract)

    [23]Kool D M, Chung H, Tata K R, et al. Hierarchical saturation of soil carbon pools near a natural CO2spring[J]. Global Change Biology, 2007, 13(6): 1282-1293.

    [24]Li H, Han X, Wang F, et al. Impact of soil management on organic carbon content and aggregate stability[J]. Communications in Soil Science and Plant Analysis, 2007, 38(13-14): 1673-1690.

    [25]Chen Y, Zhang X, He H, et al. Carbon and nitrogen pools in different aggregates of a Chinese Mollisol as influenced by long-term fertilization[J]. Journal of Soils and Sediments, 2010, 10(6): 1018-1026.

    [26]Amelung W, Zech W, Zhang X, et al. Carbon, nitrogen and sulfur pools in particle size fractions as influenced by climate[J]. Soil Science Society of America Journal, 1998, 62(1): 172-181.

    [27]林誠,鄭祥洲,郭寶玲,等. 亞熱帶地區(qū)不同種植年限果園土壤團(tuán)聚體結(jié)構(gòu)及有機(jī)碳、氮分布特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2019,38(7):1560-1566. Lin Cheng, Zheng Xiangzhou, Guo Baoling, et al. Characteristics of soil aggregate structure and organic carbon nitrogen distribution in orchards of different planting age in subtropical areas[J]. Journal of Agro-Enviroment Science, 2019, 38(7): 1560-1566. (in Chinese with English abstract)

    [28]Kuzyakov Y. Priming effects: Interactions between living and dead organic matter[J]. Soil Biology and Biochemistry, 2010, 42(9): 1363-1371.

    [29]邵興芳. 長期有機(jī)培肥模式下黑土團(tuán)聚體碳氮積累及礦化特征[D]. 武漢:武漢理工大學(xué),2014. Shao Xingfang. Carbon and Nitrogen Accumulation and Mineralization in Aggregates under Long-term Manure Fertilization Practices on Black Soil[D]. Wuhan: Wuhan University of Technology, 2014. (in Chinese with English abstract)

    [30]楊洪波,申艷,徐明崗,等. 長期施肥下黑土不同團(tuán)聚體氮組分的植物有效性差異[J]. 植物營養(yǎng)與肥料學(xué)報,2018,24(6):1581-1587. Yang Hongbo, Shen Yan, Xu Minggang, et al. Difference in plant availability of nitrogen in aggregates of black soil under long-term fertilization[J]. Plant Nutrition and Fertilizer Science, 2018, 24(6): 1581-1587. (in Chinese with English abstract)

    [31]Craswell E T, Saffigna P G, Waring S A. The mineralization of organic nitrogen in dry soil aggregates of different sizes[J]. Plant and Soil, 1970, 33(1–3): 383–392.

    [32]Pan G, Zhou P, Li Z, et al. Combined inorganic/organic fertilization enhances N efficiency and increases rice productivity through organic carbon accumulation in a rice paddy from the Tai Lake region, China[J]. Agriculture, Ecosystems & Environment, 2009, 131(3-4): 274-280.

    [33]Li Z, Liu M, Wu X, et al. Effects of long-term chemical fertilization and organic amendments on dynamics of soil organic C and total N in paddy soil derived from barren land in subtropical China[J]. Soil and Tillage Research, 2010, 106(2): 268-274.

    [34]Xie J, Peng B, Wang R, et al. Responses of crop productivity and physical protection of organic carbon by macroaggregates to long-term fertilization of an anthrosol[J]. European Journal of Soil Science, 2018, 69(3): 555-567.

    [35]Lal R. Beyond copenhagen: Mitigating climate change and achieving food security through soil carbon sequestration[J]. Food Security, 2010, 2(2): 169–177.

    [36]Wang R, Zhou J, Xie J, et al. Carbon sequestration in irrigated and rain-fed cropping systems under long-term fertilization regimes[J]. Journal of Soil Science and Plant Nutrition, DOI. org /10. 1007/s42729-020-00181-6.

    Effects of fertilization regimes on carbon and nitrogen contents of aggregates and maize yield in reclaimed soils

    Cao Hanbing1, Xie Junyu1,2, QiangJiu Ciren3, Guo Lu4, Hong Jianping1, Jing Yaodong1, Meng Huisheng1※

    (1.,,030801,; 2.,030031,; 3.,856700,; 4.,256600,)

    This study aims to investigate the distribution characteristics of organic carbon (OC) and total nitrogen (TN) in water-stable aggregates, silt and clay fractions, in order to gain a deep understanding the sequestration mechanism of soil organic carbon (SOC) under various fertilization regimes. Topsoil samples (0-20 cm) were collected from the raw soil and reclaimed soil under different fertilization regimes of six consecutive ages. A wet sieving method was used to fractionate the particle size in soil, including the large macroaggregates (> 2 mm), small macroaggregates (0.25-2 mm), microaggregates (0.053-0.25 mm), as well as the silt- and clay-sized fractions (< 0.053 mm). The contents of OC and TN were analyzed for the water-stable aggregates, the silt and clay fractions. A systematic experiment was also carried out to reveal the driving factors for the changes of SOC storage, while, the relationship between crop yields and the OC content in the water-stable aggregates, the silt and clay fractions. Four treatments were selected, including no fertilizer (control, CK), nitrogen, phosphorus, and potassium fertilizer (NPK), single application of organic manure (M), and the combined application of organic manure and synthetic fertilizers (MNPK). The results showed that the contents of OC and TN were similar in macroaggregates and microaggregates, particularly higher than that in the silt and clay fractions. Compared with the raw soil, the CK, NPK and M treatments significantly increased the OC contents in the 0.25-2 mm, 0.053-0.25 mm aggregates, and the < 0.053 mm fractions, with the increase of 97.96%, 66.19% and 97.60% (CK), 87.07%, 92.09%, 92.80% (NPK), and 95.24%, 148.20%, 112.00% (M), respectively. The MNPK treatment significantly increased the OC contents in > 2, 0.25-2, 0.053-0.35 mm aggregates, and the < 0.053 mm fractions, with the increase of 151.88%, 160.54%, 135.25%, and 92.80%, respectively. In addition, the CK, NPK, M and MNPK treatments significantly increased the TN contents in > 2 mm, 0.25-2 mm and 0.053-0.25 mm aggregates, with the increase of 43.66%, 46.78%, and 29.95% (CK), 41.83%, 52.35%, 30.43% (NPK), 52.81%, 50.43%, 44.93% (M), and 63.80%, 62.96%, 39.61% (MNPK), respectively. The enrichment factors of organic carbon (E) and total nitrogen (E) were similar in the macroaggregates and microaggregates, generally higher than those in silt and clay fractions. Both the M and MNPK treatments significantly increased the C/N ratios of water-stable aggregates, and the silt and clay fractions (except for > 2 mm aggregates). The OC sequestrated in the > 2 mm aggregates was driven by the proportions of > 2 mm, whereas, the OC sequestrated in the other water-stable aggregates, and the silt and clay fractions was driven by the OC contents. The crop yields increased linearly as the increase of OC contents in the water-stable aggregates, and the silt and clay fractions, indicating that the SOC pool in this area has not saturated, especially for the < 0.053 mm fraction. Therefore, the application of organic manure combined with NPK (MNPK) to the reclaimed soil with low fertility can be an effective way to improve soil performance in coal mining subsidence areas.

    reclamation; soils; water-stable aggregates; organic carbon; total nitrogen; crop yield

    曹寒冰,謝鈞宇,強(qiáng)久次仁,等. 施肥措施對復(fù)墾土壤團(tuán)聚體碳氮含量和作物產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(18):135-143.doi:10.11975/j.issn.1002-6819.2020.18.017 http://www.tcsae.org

    Cao Hanbing, Xie Junyu, QiangJiu Ciren, et al. Effects of fertilization regimes on carbon and nitrogen contents of aggregates and maize yield in reclaimed soils[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(18): 135-143. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.18.017 http://www.tcsae.org

    2020-06-18

    2020-08-18

    國家自然科學(xué)基金(41807102;U1710255-3)、山西省高等學(xué)校科技創(chuàng)新項目(2019L0363)、山西農(nóng)業(yè)大學(xué)科技創(chuàng)新基金(青年科技創(chuàng)新)項目(2019004)和山西省優(yōu)秀博士來晉工作獎勵(SXYBKY201805)資助

    曹寒冰,博士,講師,研究方向為高效施肥與土壤質(zhì)量提升。Email:caohanbing119@163.com

    孟會生,博士,副教授,研究方向為礦區(qū)土壤治理與生物修復(fù)。Email:huishengmeng@126.com

    10.11975/j.issn.1002-6819.2020.18.017

    S157.1

    A

    1002-6819(2020)-18-0135-09

    12—13女人毛片做爰片一| 国产欧美日韩综合在线一区二区| av又黄又爽大尺度在线免费看| 自拍欧美九色日韩亚洲蝌蚪91| 50天的宝宝边吃奶边哭怎么回事| 97人妻天天添夜夜摸| 国产精品影院久久| 久久精品人人爽人人爽视色| 中文字幕人妻熟女乱码| 欧美久久黑人一区二区| 女同久久另类99精品国产91| 亚洲第一av免费看| 满18在线观看网站| 国产亚洲欧美精品永久| 又黄又粗又硬又大视频| 国产一区有黄有色的免费视频| 大香蕉久久网| 视频区图区小说| 50天的宝宝边吃奶边哭怎么回事| 不卡一级毛片| a级片在线免费高清观看视频| 国产精品av久久久久免费| 午夜精品久久久久久毛片777| 高清在线国产一区| 欧美日韩av久久| 国产精品麻豆人妻色哟哟久久| 色94色欧美一区二区| 一区在线观看完整版| 国产人伦9x9x在线观看| 建设人人有责人人尽责人人享有的| 国产成人系列免费观看| 亚洲七黄色美女视频| 99热网站在线观看| 亚洲自偷自拍图片 自拍| 久久精品国产a三级三级三级| 精品人妻熟女毛片av久久网站| 午夜福利在线观看吧| 午夜激情久久久久久久| 97人妻天天添夜夜摸| 欧美在线黄色| 母亲3免费完整高清在线观看| 一本色道久久久久久精品综合| 色播在线永久视频| 精品亚洲成国产av| 久久人妻熟女aⅴ| 菩萨蛮人人尽说江南好唐韦庄| 91精品国产国语对白视频| 另类精品久久| 国产一区二区激情短视频| 欧美乱妇无乱码| 欧美大码av| 国产精品免费视频内射| 欧美激情 高清一区二区三区| 色精品久久人妻99蜜桃| 欧美国产精品一级二级三级| 亚洲成人免费电影在线观看| 久久毛片免费看一区二区三区| 人成视频在线观看免费观看| 男女高潮啪啪啪动态图| 777久久人妻少妇嫩草av网站| 搡老岳熟女国产| 久久人妻福利社区极品人妻图片| 成人三级做爰电影| 久热爱精品视频在线9| 在线观看免费视频网站a站| 亚洲七黄色美女视频| 中文字幕色久视频| 精品欧美一区二区三区在线| 亚洲国产欧美在线一区| 国产精品电影一区二区三区 | 天天躁日日躁夜夜躁夜夜| 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女 | 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 欧美日韩精品网址| 一区福利在线观看| 欧美亚洲 丝袜 人妻 在线| 久久久久精品国产欧美久久久| 国产精品一区二区精品视频观看| 亚洲精品久久成人aⅴ小说| 国产成人精品久久二区二区91| 成年女人毛片免费观看观看9 | 精品高清国产在线一区| 涩涩av久久男人的天堂| 久久国产精品大桥未久av| 亚洲伊人久久精品综合| 欧美性长视频在线观看| 欧美黑人精品巨大| 中亚洲国语对白在线视频| 欧美成人免费av一区二区三区 | 淫妇啪啪啪对白视频| 黄色视频在线播放观看不卡| 天天影视国产精品| 午夜福利一区二区在线看| 宅男免费午夜| 国产老妇伦熟女老妇高清| 正在播放国产对白刺激| 热re99久久国产66热| 亚洲精品av麻豆狂野| 亚洲成人手机| 在线观看舔阴道视频| 五月天丁香电影| 老司机靠b影院| 午夜福利欧美成人| 亚洲av电影在线进入| 夜夜爽天天搞| 亚洲国产av影院在线观看| 91大片在线观看| 日韩免费高清中文字幕av| 久热爱精品视频在线9| 午夜免费成人在线视频| av视频免费观看在线观看| 亚洲精品粉嫩美女一区| 日韩人妻精品一区2区三区| 久久九九热精品免费| 黄色片一级片一级黄色片| 国产精品欧美亚洲77777| 一区二区三区乱码不卡18| 水蜜桃什么品种好| 日韩免费av在线播放| 伊人久久大香线蕉亚洲五| 亚洲国产毛片av蜜桃av| 考比视频在线观看| 国产在线视频一区二区| 飞空精品影院首页| 欧美日韩亚洲高清精品| 国产野战对白在线观看| 久久精品国产亚洲av香蕉五月 | 女人爽到高潮嗷嗷叫在线视频| 99国产精品一区二区三区| 激情视频va一区二区三区| 中文字幕制服av| 亚洲欧美色中文字幕在线| 免费高清在线观看日韩| 国产精品偷伦视频观看了| 国产精品一区二区在线观看99| 美女午夜性视频免费| 美女高潮到喷水免费观看| 国产成人影院久久av| 精品国产乱码久久久久久小说| 欧美大码av| 女性被躁到高潮视频| 免费在线观看影片大全网站| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| 亚洲av国产av综合av卡| 亚洲第一欧美日韩一区二区三区 | 日日夜夜操网爽| 动漫黄色视频在线观看| 日韩精品免费视频一区二区三区| 精品人妻熟女毛片av久久网站| 老汉色av国产亚洲站长工具| 超色免费av| 国产精品欧美亚洲77777| 欧美日韩国产mv在线观看视频| 亚洲中文字幕日韩| 久久精品aⅴ一区二区三区四区| 亚洲人成电影观看| 一区二区三区乱码不卡18| 一区福利在线观看| 69精品国产乱码久久久| 亚洲欧美一区二区三区久久| a级片在线免费高清观看视频| 中文字幕人妻熟女乱码| 手机成人av网站| 宅男免费午夜| 国产精品亚洲av一区麻豆| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 狠狠婷婷综合久久久久久88av| 国产xxxxx性猛交| 国产一区二区三区综合在线观看| 母亲3免费完整高清在线观看| 成人三级做爰电影| 精品少妇内射三级| 黄网站色视频无遮挡免费观看| 久久毛片免费看一区二区三区| 这个男人来自地球电影免费观看| 成人国语在线视频| 啦啦啦视频在线资源免费观看| 中文亚洲av片在线观看爽 | 激情视频va一区二区三区| 成年动漫av网址| 伊人久久大香线蕉亚洲五| 免费少妇av软件| 国产av精品麻豆| 亚洲人成电影免费在线| 国产aⅴ精品一区二区三区波| 亚洲欧美日韩高清在线视频 | 欧美精品人与动牲交sv欧美| 国产单亲对白刺激| 天天躁日日躁夜夜躁夜夜| 国产深夜福利视频在线观看| 欧美日本中文国产一区发布| 久久性视频一级片| 久久午夜综合久久蜜桃| 国产亚洲精品一区二区www | 美女高潮到喷水免费观看| 亚洲国产欧美在线一区| 自线自在国产av| 久久人妻熟女aⅴ| 欧美久久黑人一区二区| 国产不卡一卡二| 欧美乱妇无乱码| 一本久久精品| 亚洲精品国产精品久久久不卡| 老司机午夜福利在线观看视频 | 天天操日日干夜夜撸| 我的亚洲天堂| 色在线成人网| 亚洲天堂av无毛| 老司机深夜福利视频在线观看| 亚洲av美国av| 亚洲精品国产一区二区精华液| 丝袜在线中文字幕| 一级,二级,三级黄色视频| 美女福利国产在线| 高清视频免费观看一区二区| 老司机午夜福利在线观看视频 | 久久精品国产亚洲av高清一级| 欧美另类亚洲清纯唯美| 亚洲色图av天堂| 80岁老熟妇乱子伦牲交| a在线观看视频网站| 精品国内亚洲2022精品成人 | 高清欧美精品videossex| 极品少妇高潮喷水抽搐| 国产一区二区三区在线臀色熟女 | 国产免费av片在线观看野外av| 母亲3免费完整高清在线观看| 中文字幕最新亚洲高清| 亚洲中文字幕日韩| 一二三四社区在线视频社区8| 午夜精品久久久久久毛片777| 啦啦啦视频在线资源免费观看| 91字幕亚洲| 国产一区二区三区在线臀色熟女 | 久久av网站| 欧美乱妇无乱码| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| 亚洲av日韩精品久久久久久密| 国产成人一区二区三区免费视频网站| 亚洲欧美一区二区三区久久| 天天影视国产精品| 母亲3免费完整高清在线观看| 欧美大码av| 丝袜在线中文字幕| 久久国产精品大桥未久av| av片东京热男人的天堂| 12—13女人毛片做爰片一| 交换朋友夫妻互换小说| 国产精品久久久久成人av| 日本精品一区二区三区蜜桃| 黑人猛操日本美女一级片| 国产在线一区二区三区精| 超色免费av| 不卡av一区二区三区| 黄色视频不卡| 在线播放国产精品三级| 国产淫语在线视频| 久久久久视频综合| 日本五十路高清| 中文欧美无线码| 在线十欧美十亚洲十日本专区| 捣出白浆h1v1| 中文字幕最新亚洲高清| 久久精品aⅴ一区二区三区四区| aaaaa片日本免费| 国产一区二区三区视频了| 国产精品久久久久久人妻精品电影 | 精品久久久精品久久久| 深夜精品福利| 91大片在线观看| 日本五十路高清| 国产精品 国内视频| 国产亚洲欧美精品永久| 国产精品熟女久久久久浪| 国产不卡av网站在线观看| 免费人妻精品一区二区三区视频| 亚洲情色 制服丝袜| 欧美成人午夜精品| 欧美午夜高清在线| 国产精品成人在线| 成人精品一区二区免费| 国产男靠女视频免费网站| 一夜夜www| a在线观看视频网站| 脱女人内裤的视频| 日本五十路高清| 狠狠狠狠99中文字幕| 久久久精品94久久精品| 欧美 日韩 精品 国产| videosex国产| 一本色道久久久久久精品综合| 一二三四在线观看免费中文在| 精品国产一区二区三区久久久樱花| 91老司机精品| 人妻一区二区av| 免费久久久久久久精品成人欧美视频| 久久国产精品男人的天堂亚洲| bbb黄色大片| 亚洲精品在线美女| 欧美日韩成人在线一区二区| 欧美日韩av久久| 肉色欧美久久久久久久蜜桃| 欧美乱码精品一区二区三区| 黑丝袜美女国产一区| 亚洲 欧美一区二区三区| 国产不卡av网站在线观看| 色综合欧美亚洲国产小说| 久久精品国产a三级三级三级| 国产欧美日韩一区二区精品| 欧美亚洲日本最大视频资源| tocl精华| 999久久久精品免费观看国产| 一级,二级,三级黄色视频| 老司机深夜福利视频在线观看| 国产午夜精品久久久久久| 国产精品一区二区免费欧美| 人妻一区二区av| 国产成人啪精品午夜网站| 亚洲中文字幕日韩| 亚洲欧美一区二区三区黑人| 两性午夜刺激爽爽歪歪视频在线观看 | 激情视频va一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲成a人片在线一区二区| 天天躁夜夜躁狠狠躁躁| 一进一出抽搐动态| 好男人电影高清在线观看| 男人操女人黄网站| 国产高清视频在线播放一区| 久9热在线精品视频| 国产亚洲一区二区精品| 亚洲人成电影观看| 精品国产超薄肉色丝袜足j| 色精品久久人妻99蜜桃| 女同久久另类99精品国产91| 高清毛片免费观看视频网站 | 亚洲成国产人片在线观看| 亚洲伊人色综图| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 麻豆国产av国片精品| 999久久久精品免费观看国产| 高清毛片免费观看视频网站 | 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 午夜激情久久久久久久| 18禁观看日本| 国产欧美日韩一区二区精品| 亚洲国产成人一精品久久久| 国产成人一区二区三区免费视频网站| 国产在线一区二区三区精| 亚洲色图 男人天堂 中文字幕| 中文字幕精品免费在线观看视频| 国产亚洲精品第一综合不卡| 国产亚洲一区二区精品| 精品国产一区二区三区久久久樱花| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 黄色毛片三级朝国网站| 男女无遮挡免费网站观看| 国产一区二区在线观看av| 12—13女人毛片做爰片一| 18在线观看网站| 亚洲午夜精品一区,二区,三区| 国产三级黄色录像| 正在播放国产对白刺激| 99久久人妻综合| 国产黄频视频在线观看| 亚洲自偷自拍图片 自拍| 久久久水蜜桃国产精品网| 欧美大码av| 久久国产精品男人的天堂亚洲| 三级毛片av免费| 999久久久精品免费观看国产| 51午夜福利影视在线观看| 国产成人av激情在线播放| 国产在线观看jvid| 久久精品91无色码中文字幕| 欧美 亚洲 国产 日韩一| 亚洲成av片中文字幕在线观看| 国产单亲对白刺激| 美国免费a级毛片| 十八禁网站网址无遮挡| 欧美中文综合在线视频| 中文字幕制服av| 久久久久久久精品吃奶| 热re99久久国产66热| 亚洲av美国av| www.精华液| 波多野结衣av一区二区av| 男人操女人黄网站| 黑人巨大精品欧美一区二区mp4| 99精国产麻豆久久婷婷| 天堂8中文在线网| 80岁老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆| 欧美 日韩 精品 国产| 中亚洲国语对白在线视频| 99九九在线精品视频| 国产91精品成人一区二区三区 | 在线十欧美十亚洲十日本专区| 亚洲,欧美精品.| 高清av免费在线| 亚洲av成人不卡在线观看播放网| 欧美激情高清一区二区三区| 国产在线精品亚洲第一网站| 新久久久久国产一级毛片| 大香蕉久久成人网| 99国产精品一区二区三区| 女人久久www免费人成看片| 中文字幕最新亚洲高清| 啦啦啦免费观看视频1| 国产麻豆69| 亚洲精品中文字幕在线视频| 国产精品 国内视频| 我要看黄色一级片免费的| 大香蕉久久网| 久久热在线av| 国产精品香港三级国产av潘金莲| bbb黄色大片| 激情视频va一区二区三区| av网站免费在线观看视频| 一区二区三区乱码不卡18| 久久久水蜜桃国产精品网| av天堂久久9| 亚洲国产看品久久| 亚洲伊人久久精品综合| 丝袜喷水一区| 在线播放国产精品三级| 日本黄色日本黄色录像| 久久精品国产综合久久久| 欧美激情高清一区二区三区| 久久av网站| 成人18禁高潮啪啪吃奶动态图| 高清视频免费观看一区二区| 丝袜美足系列| 久久99一区二区三区| 黄色成人免费大全| 亚洲熟女精品中文字幕| 午夜免费成人在线视频| 免费久久久久久久精品成人欧美视频| 亚洲男人天堂网一区| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| 水蜜桃什么品种好| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区蜜桃| 18禁裸乳无遮挡动漫免费视频| 男女无遮挡免费网站观看| 国产精品一区二区在线不卡| 大香蕉久久网| 热99国产精品久久久久久7| 男女边摸边吃奶| 天堂俺去俺来也www色官网| 69精品国产乱码久久久| 在线观看人妻少妇| 桃花免费在线播放| 免费看a级黄色片| 99久久国产精品久久久| 热re99久久精品国产66热6| 久久精品91无色码中文字幕| 国产精品98久久久久久宅男小说| 男女下面插进去视频免费观看| av电影中文网址| 91成年电影在线观看| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 久久人妻av系列| 久久天堂一区二区三区四区| 久久久久精品国产欧美久久久| 成年动漫av网址| 老司机福利观看| 色在线成人网| 亚洲精品乱久久久久久| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 国产亚洲精品第一综合不卡| 欧美一级毛片孕妇| 欧美精品高潮呻吟av久久| 操美女的视频在线观看| 亚洲美女黄片视频| 老司机午夜福利在线观看视频 | 91av网站免费观看| 免费在线观看视频国产中文字幕亚洲| 香蕉丝袜av| 久久精品亚洲精品国产色婷小说| 18禁美女被吸乳视频| 国产精品亚洲一级av第二区| 99精品久久久久人妻精品| 一个人免费看片子| 国产日韩一区二区三区精品不卡| 黑丝袜美女国产一区| 国产精品99久久99久久久不卡| 欧美精品一区二区免费开放| 久久久国产欧美日韩av| 精品第一国产精品| 午夜福利免费观看在线| 老熟妇乱子伦视频在线观看| 欧美日韩一级在线毛片| 久久狼人影院| 一级毛片精品| cao死你这个sao货| 12—13女人毛片做爰片一| 日韩中文字幕视频在线看片| 中文字幕人妻丝袜一区二区| 亚洲专区字幕在线| 两性夫妻黄色片| 午夜福利视频在线观看免费| 国产精品一区二区在线不卡| 国产av精品麻豆| 在线观看免费视频日本深夜| 99热网站在线观看| 久久中文看片网| 在线观看舔阴道视频| 又大又爽又粗| videos熟女内射| 成人18禁高潮啪啪吃奶动态图| 国产精品自产拍在线观看55亚洲 | 亚洲视频免费观看视频| netflix在线观看网站| 国产精品 国内视频| 亚洲 国产 在线| 国产精品欧美亚洲77777| 亚洲中文av在线| 国产精品一区二区精品视频观看| 国产成人一区二区三区免费视频网站| av视频免费观看在线观看| 热re99久久国产66热| 亚洲精品自拍成人| 乱人伦中国视频| 久久ye,这里只有精品| 黄色视频,在线免费观看| 侵犯人妻中文字幕一二三四区| 视频区图区小说| 欧美变态另类bdsm刘玥| 王馨瑶露胸无遮挡在线观看| 欧美日韩成人在线一区二区| 一个人免费看片子| 怎么达到女性高潮| 人成视频在线观看免费观看| 大型黄色视频在线免费观看| 一级毛片电影观看| 黑丝袜美女国产一区| 黄网站色视频无遮挡免费观看| 欧美黄色淫秽网站| 一二三四社区在线视频社区8| 女性被躁到高潮视频| 日韩有码中文字幕| 国产免费av片在线观看野外av| 9191精品国产免费久久| 大香蕉久久成人网| 日本欧美视频一区| 国产1区2区3区精品| 老鸭窝网址在线观看| 99久久人妻综合| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 12—13女人毛片做爰片一| 国产精品免费大片| 老司机影院毛片| 久久久精品区二区三区| 国产成人欧美在线观看 | av福利片在线| 日韩欧美一区视频在线观看| 亚洲中文日韩欧美视频| 亚洲av国产av综合av卡| 一个人免费看片子| 日日夜夜操网爽| 色精品久久人妻99蜜桃| 91成人精品电影| 黄片播放在线免费| 人成视频在线观看免费观看| 国产精品影院久久| 91老司机精品| 天天操日日干夜夜撸| 国产在线一区二区三区精| 国产在线精品亚洲第一网站| 国产成人精品久久二区二区91| 伊人久久大香线蕉亚洲五| 91大片在线观看| 亚洲精品中文字幕一二三四区 | 亚洲avbb在线观看| 免费高清在线观看日韩| 热re99久久精品国产66热6| 欧美日韩亚洲国产一区二区在线观看 | 久久狼人影院| 王馨瑶露胸无遮挡在线观看| 午夜福利在线观看吧| 欧美 日韩 精品 国产| 日韩大片免费观看网站| 久久久欧美国产精品| 黄色毛片三级朝国网站| 肉色欧美久久久久久久蜜桃| 国产精品1区2区在线观看. | 丝袜美腿诱惑在线| 老熟妇乱子伦视频在线观看| 香蕉国产在线看| 99国产极品粉嫩在线观看| 高清av免费在线| 99re6热这里在线精品视频| 99九九在线精品视频| 久久久久久久久久久久大奶| 婷婷成人精品国产| 操美女的视频在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| 亚洲七黄色美女视频| 热re99久久精品国产66热6| 操美女的视频在线观看|