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

    秸稈還田對黑土亞表層微生物群落結(jié)構(gòu)的影響特征及原因分析

    2020-03-03 13:54:38董建新李玉義逄煥成
    農(nóng)業(yè)工程學報 2020年1期
    關鍵詞:高量表層群落

    叢 萍,王 婧,董建新,李玉義,劉 娜,逄煥成

    秸稈還田對黑土亞表層微生物群落結(jié)構(gòu)的影響特征及原因分析

    叢 萍1,2,王 婧1,董建新2,李玉義1,劉 娜1,逄煥成1※

    (1. 中國農(nóng)業(yè)科學院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081;2. 中國農(nóng)業(yè)科學院煙草研究所,青島 266100)

    為明確切碎秸稈與秸稈顆粒對黑土亞表層土壤微生物群落結(jié)構(gòu)的影響效應,從而評價不同秸稈還田方式對亞表層的培肥效果,該研究于2016—2018年在東北黑土區(qū)進行一次性玉米秸稈深埋還田試驗,設置切碎秸稈低量(QS1)、切碎秸稈高量(QS5)、秸稈顆粒低量(KL1)與秸稈顆粒高量(KL5)4種秸稈還田處理,并與秸稈不還田(CK)進行對比,于每年玉米收獲季對土壤理化指標及微生物菌群結(jié)構(gòu)進行監(jiān)測。結(jié)果表明,1)秸稈還田第1年,切碎秸稈處理顯著提高土壤總磷脂脂肪酸含量及真菌摩爾百分數(shù),其高量處理較CK最高增加71.0%和120.5%,而秸稈顆粒處理對細菌、革蘭氏陽性菌和革蘭氏陰性菌的摩爾百分數(shù)增幅更顯著,其高量處理最高增加41.6%、29.7%和26.3%;還田第2年高量處理顯著提高各菌群磷脂脂肪酸含量,且切碎高量處理的真菌摩爾百分數(shù)含量顯著高于顆粒高量處理21.0%;還田第3年僅高量處理下的菌群結(jié)構(gòu)有顯著分異。2)還田初期切碎秸稈處理顯著提高真菌:細菌比值,而低量還田則對革蘭氏陽性菌和革蘭氏陰性菌比有提高,隨還田時間的增加,高量還田比值提高更顯著,利于長期維持生態(tài)系統(tǒng)穩(wěn)定性。3)秸稈高量還田可顯著改變土壤理化因子水平,這是導致微生物群落結(jié)構(gòu)分異的重要原因,其顯著影響因子隨還田年限而更替:第1年,土壤容重、酸堿度、全氮和碳氮比為顯著(<0.05)影響因素;第2年,土壤含水量、有機碳、碳氮比和土壤容重為極顯著影響因素(<0.01);第3年僅有機碳為顯著因素(<0.05)。切碎秸稈高量還田處理微生物群落結(jié)構(gòu)分布與CK區(qū)分最為明顯,對土壤真菌群落的調(diào)控能力更強,更適宜于東北黑土亞表層肥力的提升。

    土壤;秸稈;微生物群落結(jié)構(gòu);深還田;亞表層

    0 引 言

    土壤微生物是土壤生態(tài)系統(tǒng)中最活躍的具有生命特征的重要組成成分,參與生態(tài)系統(tǒng)中的各個養(yǎng)分生化循環(huán)過程[1-2],土壤微生物多樣性以及結(jié)構(gòu)變化已成為衡量土壤肥力質(zhì)量的重要指標[3-4]。秸稈還田是當前增加農(nóng)地有機物料投入的重要手段之一,可以培肥土壤、提升地力。秸稈施入到土壤中對微生物群落結(jié)構(gòu)及多樣性有重要影響,一方面作為高碳源有機物質(zhì)為土壤微生物生命活動提供能量,另一方面,秸稈還田后改變了土壤的孔隙結(jié)構(gòu)及水分狀況等,影響了微生物的生存環(huán)境[5-6]??梢姡斩掃€田易于影響微生物生境中的敏感性因子,研究土壤中的微生物變化有利于更好地評價秸稈還田的培肥效果。

    亞表層是位于心土層(20~50 cm)范圍內(nèi)的重要耕作層以及作物根系分布層,位于約20~40 cm,對土壤碳固定以及作物產(chǎn)量提高有重要作用[7]。然而當前東北黑土犁底層增厚變硬(增厚至7~11 cm)[8]、亞表層有機碳含量銳減(較墾前降低40%~60%[9]、黑土層變薄(近40%已小于30 cm)[10-11]等問題嚴峻,亟需有效措施進行培肥改良。秸稈深還田是培肥亞表層的重要措施[12-14],它直接將秸稈施入30~40 cm土層,使土壤容重降低9.9%~6.2%[13],有機碳增加10%~20%,并且改善效果在一定程度上隨著秸稈用量的增加而增大[14]。然而,當前有關秸稈還田對亞表層微生物群落結(jié)構(gòu)影響的研究甚少,土壤生物學肥力的改善效果未知。

    研究表明,秸稈用量、秸稈形態(tài)以及還田年限等對土壤微生物均有不同影響。秸稈用量因耕作方式以及氣候類型而異,東北玉米大豆輪作區(qū)免耕條件下60%的秸稈還田量更有利于提高微生物活性[15],而在南方水稻產(chǎn)區(qū)耕翻條件下2/3還田量處理的土壤微生物活度最高,少免耕條件下1/3還田量處理最高[16]。秸稈形態(tài)越細小,粉碎化程度越高,越有利于加快微生物的腐解進程[17],并且能夠影響細菌及真菌的群落組成[18]。在中國北方地區(qū),長期玉米秸稈還田能夠顯著提高真菌生物量,而對細菌無顯著影響[19],但是短期的秸稈還田則能夠增加土壤總磷脂脂肪酸(PLFA)數(shù)量、真菌以及細菌生物量[20]。但上述有關秸稈還田對土壤微生物影響的研究往往只針對于0~20 cm耕層,且還田量大多低于全量還田,致使玉米秸稈的肥料化利用率仍然偏低。收割機切碎后的秸稈直接還田操作相對簡便,是農(nóng)民常用的秸稈還田形式嗎,但秸稈形態(tài)較為粗放,可能無法在短期內(nèi)快速培肥,而新型的秸稈顆粒肥是用常規(guī)秸稈細粉后擠壓制成,具有體積小、施用方便、養(yǎng)分釋放率高的優(yōu)點,且前期試驗驗證其具有快速培肥耕層土壤的效應[21-22]。對于土壤亞表層而言,其氧氣以及水分條件等均異于耕層[23],當不同形態(tài)或者是大量秸稈施入亞表層后,碳源的增加及其他理化性狀的改變必然引起土壤微生物群落的變化,但其影響特征以及改變原因尚不清晰。因此,本研究將收割機切碎的秸稈(QS)以及秸稈顆粒(KL)2種形態(tài)的秸稈實施深還田,在前期研究的基礎上設置高、低兩種倍量,并與秸稈不還田(CK)進行對比。研究目的在于:一是探究不同秸稈形態(tài)以及倍量還田下亞表層土壤細菌及真菌群落結(jié)構(gòu)的變化,二是明確東北黑土區(qū)影響亞表層微生物群落變化的主要環(huán)境因子,以期從土壤微生物學角度為東北黑土亞表層的培肥方式提供指導。

    1 材料和方法

    1.1 試驗地概況

    田間定位試驗于2015-2018年在吉林省公主嶺市朝陽坡鎮(zhèn)(43°36′N,124°43′E)進行,該地位于我國東北黃金玉米生產(chǎn)帶上,地處吉林省中西部,東遼河中游右岸。當?shù)匾杂牮B(yǎng)農(nóng)業(yè)為主,種植作物主要為玉米、大豆。屬溫帶大陸性季風氣候,年平均氣溫5.6 ℃,年均降水量594.8 mm,無霜期144 d。試驗地土壤類型為典型黑土,耕作方式以旋耕、翻耕為主。土壤基礎理化性狀如表1。

    表1 2015年10月試驗地土壤基本理化性狀

    1.2 試驗材料

    本試驗設置玉米切碎秸稈、玉米秸稈顆粒2種不同秸稈形態(tài),以當?shù)亟斩挸A窟€田的用量為基準(玉米收獲后剩余地上部的總量,每年約為15 t/hm2),分別設置玉米秸稈切碎還田(QS1, 15 t/hm2)、玉米秸稈切碎還田(QS5, 75 t/hm2)、玉米秸稈顆粒還田(KL1, 15 t/hm2)、玉米秸稈顆粒還田(KL5, 75 t/hm2)4種秸稈還田處理,并以秸稈不還田(CK)為對照。每個處理3次重復,共設15個小區(qū),隨機區(qū)組排列,每個小區(qū)為長4 m、寬4 m(6壟)共計16 m2見方,小區(qū)與小區(qū)之間留出1 m距離,避免因農(nóng)機具拖曳秸稈而對小區(qū)處理造成影響。試驗布置前需收集當?shù)赜衩资斋@時隨收割機切碎的玉米秸稈,長度約10~15 cm,將收集起的秸稈分為2部分,一部分風干后用作切碎秸稈還田處理(QS),另一部分收集起來的玉米秸稈用HC-2000型粉碎機粉碎后過2 mm篩,按照30%~35%比例添加蒸餾水,攪拌均勻后用FTHBCX350型環(huán)模秸稈顆粒機擠壓制得直徑4 mm、長度4~6 cm的棒狀固體,即秸稈顆粒(KL),不同形態(tài)秸稈的主要理化性狀如表2所示。

    表2 不同形態(tài)秸稈的主要理化性狀

    1.3 試驗設計

    在2015年10月前茬玉米收獲后實施秸稈翻埋,將收集到的切碎秸稈以及秸稈顆粒按照規(guī)定用量均勻鋪于各個小區(qū)內(nèi),由拖拉機牽引翻埋犁(功率>92 kW)統(tǒng)一對整片試驗田進行翻埋耕作,從而避免對相鄰小區(qū)的碾壓。秸稈翻埋犁作業(yè)深度約為30~40 cm,可將絕大部分切碎秸稈以及秸稈顆粒翻入20~40 cm亞表層范圍。此后3年試驗期間不再進行任何秸稈還田,每年玉米播種前需清除地表殘留秸稈。當?shù)胤N植制度為玉米連作一年一熟制,供試玉米品種為奇玉8號,通常為每年5月初進行播種,種植密度為67 500 株/hm2,播種前進行側(cè)面深施基肥,施用量為純N 225 kg/hm2、P2O575 kg/hm2和K2O 225 kg/hm2。其中磷肥、鉀肥一次播種,氮肥30%基施,70%拔節(jié)施,各處理的田間管理與普通高產(chǎn)玉米相似,并于每年10月初收獲。

    1.4 測定指標與方法

    1.4.1 土壤樣品采集

    2016-2018年連續(xù)3年在夏玉米成熟期(10月初),在距離小區(qū)邊緣0.5 m的內(nèi)部依據(jù)五點取樣原則用不銹鋼土鉆取20~40 cm土層的土樣,混勻作為一個重復。將土壤樣品從土鉆中取出后,剔除可見的掉落物和根系,將土塊輕輕捏碎、混勻,在田間過2 mm篩[24],一部分新鮮土樣立即放在保溫盒內(nèi)用冰袋保存運回實驗室,冷凍干燥后放在?80 ℃保存,用于磷脂脂肪酸(PLFA)法測定微生物群落結(jié)構(gòu);另一部分土樣放入自封袋中帶回實驗室風干,用于測定土壤理化因子指標。

    1.4.2 磷脂脂肪酸(PLFA)的測定

    土壤微生物群落采取磷脂脂肪酸(PLFA)方法進行測定。磷脂脂肪酸的提取方法主要參照Frostegard等[25]的方法。采用美國MIDI Sherlock微生物鑒定系統(tǒng)平臺(Version 4.5, MIDI, Inc, Newark, DE)的PLFA模塊,以氫氣為載氣,以安捷倫Agilent 6890N氣象色譜儀及FID檢測器為硬件平臺,色譜柱為Agilent 19091B?102(25 m× 200m× 0.33m),每個樣品進樣量為1L。根據(jù)MIDI平臺規(guī)范,所有測試均采用0.5 mg/mL十九酸甲酯標準品(冷泉港生物科技股份有限公司)進行校正。脂肪酸定量用峰面積和內(nèi)標曲線法表征,單位nmol/g。

    以磷脂脂肪酸16:0、17:0、16:1ω5c、16:1ω7c、16:1ω9c、17:1ω8c、18:1ω5c、18:1ω7c、a15:0、a17:0、cy17:0、cy19:0ω8c、i14:0、i15:0、i16:0、i17:0和i19:0表征細菌[26]。用18:2ω6,9c、18:1ω9c和18:3ω6c(6,9,12)表征真菌[27]。以16:0(10Me)、17:0(10Me)和18:0(10Me)表征放線菌。以14:0、17:0、18:0、i14:0、i15:0、i16:0、i17:0、a15:0、a17:0和cy17:0、cy19:0ω8c、16:1ω5c、16:1ω9c、16:1ω7c、17:1ω8c、18:1ω5c、18:1ω7c分別表征革蘭氏陽性菌(G+)和革蘭氏陰性菌(G-)[28]。

    1.4.3 環(huán)境因子的測定

    秸稈還田后直接改變土壤容重(soil bulk density,SBD,g/cm3)、土壤含水量(soil water content,SWC,g/kg)、pH值、土壤有機碳含量(soil organic carbon,SOC,g/kg)、土壤全氮(total nitrogen,TN,g/kg)以及土壤C/N,因此選擇此6項作為主要的環(huán)境因子。其中,SBD在田間運用環(huán)刀法對20~40 cm土層進行測定;SWC則運用經(jīng)典烘干法對20~40 cm土層進行測定;土壤pH運用FE38-FiveEasyPlus?(Mettler-Toledo, Switzerland)型pH計進行測定;土壤有機質(zhì)(soil organic matter,SOM,g/kg)采用重鉻酸鉀外加熱法測定,再根據(jù)SOM(g/kg)=SOC(g/kg)×1.724來換算成有機碳的含量;TN采用凱氏定氮法進行測定;C/N則是SOC與TN的比值。

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

    采用Microsoft Excel 2016對試驗數(shù)據(jù)進行整理,運用SAS 9.4軟件中的單因素方差分析(one-way ANOVA)比較不同處理間的差異(LSD,=0.05)。采用Canoco 5.0軟件對不同處理間微生物群落分布差異進行主成分分析(principal component analysis, PCA),對微生物群落結(jié)構(gòu)變化與環(huán)境因子的相互關系進行冗余分析(redundancy analysis, RDA)。

    2 結(jié)果與分析

    2.1 秸稈還田對亞表層土壤微生物群落的影響

    2.1.1 亞表層土壤微生物磷脂脂肪酸(PLFA)的變化特征

    在黑土亞表層中微生物群落數(shù)量的大小順序為細菌最大,真菌與放線菌次之。與CK相比,秸稈還田前2年顯著增加亞表層土壤磷脂脂肪酸(PLFA)總量以及細菌、真菌的摩爾百分數(shù),但能夠降低放線菌摩爾百分數(shù),降低幅度隨還田時間的增加略有減?。ū?)。主要表現(xiàn)如下:

    表3 2016-2018年不同處理下亞表層土壤(20~40 cm)PLFA總量以及各菌群PLFA變化

    注:CK為秸稈不還田;QS1為玉米秸稈切碎還田,15 t·hm-2;KL1為玉米秸稈顆粒還田,15 t·hm-2;QS5為玉米秸稈切碎還田,75 t·hm-2;KL5為玉米秸稈顆粒還田,75 t·hm-2。小寫字母代表不同處理間進行差異比較,差異分析在<0.05的顯著性水平,下同。

    Note: CK, no straw return; QS1, chopped straw return with 15 t·hm-2; KL1, pelletized straw return with 15 t·hm-2;QS5, chopped straw return with 75 t·hm-2; KL5, pelletized straw return with 75 t·hm-2. Different lowercase letters represent differences between treatments at the level of 0.05, the same below.

    2016年,PLFA總量變幅在18.96~32.36 nmol/g,與CK相比,QS1與QS5處理分別顯著增加亞表層土壤PLFA含量45.5%(<0.05)與71.0%(<0.05)。秸稈還田能顯著增加細菌摩爾百分數(shù),增幅為14.5%~41.6%(<0.05),以秸稈顆粒最為顯著;切碎秸稈還田能顯著增加真菌摩爾百分數(shù),QS5與QS1分別較CK提高120.5%(<0.05)與54.7%(<0.05),且QS5處理顯著高于QS1;秸稈還田能降低放線菌的摩爾百分數(shù),其中KL5降低幅度最大,約為29.4%(<0.05);KL5處理下G+與G-的摩爾百分數(shù)最高,分別高于CK 29.7%(<0.05)與26.3%(<0.05),而KL1次之,QS5最低。由此可見,在還田第1年秸稈形態(tài)對PLFA總量以及菌群結(jié)構(gòu)影響顯著,秸稈顆粒能夠顯著增加細菌摩爾百分數(shù),而切碎秸稈更易引起真菌摩爾百分數(shù)的顯著提高,且切碎秸稈對真菌的增加程度要高于秸稈顆粒之于細菌。

    2017年,PLFA總量變幅在15.83~20.50 nmol/g,與CK相比,QS5與KL5處理分別顯著增加土壤PLFA總量29.5%(<0.05)與16.2%(<0.05),但二者無顯著差異;秸稈還田處理下細菌摩爾百分數(shù)含量均顯著高于CK,表現(xiàn)為QS5>KL5與QS1>KL1,增幅為6.0%~19.9%(<0.05);從真菌含量變化看,QS5與KL5分別顯著高于CK 74.5%(<0.05)與44.2%(<0.05),且QS5亦顯著高于KL5 21.0%(<0.05);而對于放線菌而言,KL5與QS5仍表現(xiàn)出顯著降低趨勢,低于CK約24.3%;秸稈高量還田較CK顯著增加革蘭氏陽性菌(G+)與革蘭氏陰性菌(G-)的摩爾百分數(shù)約6.3%(<0.05)與13.2%(<0.05),但KL5與QS5之間無顯著差異。由此可見,秸稈還田第2年PLFA總量以及各菌群變化受秸稈用量的影響顯著,秸稈高量還田能夠顯著增加PLFA總量以及細菌、真菌、G+與G-的摩爾百分數(shù),顯著降低放線菌的摩爾百分數(shù)。

    2018年,PLFA總量變幅在15.04~18.56 nmol/g,但秸稈還田對PLFA總量及細菌、真菌與G+的摩爾百分數(shù)無顯著影響。其中,秸稈高量還田能顯著降低放線菌的摩爾百分數(shù)約17.9%(<0.05),KL5亦能顯著降低G-20.1%(<0.05)??梢婋S著還田時間增加,秸稈對土壤微生物群落結(jié)構(gòu)的影響逐漸減小,微生物群落結(jié)構(gòu)有回歸于最初結(jié)構(gòu)組成的趨勢。

    2.1.2 亞表層真菌:細菌比值以及G+:G-比值的變化特征

    土壤中真菌:細菌比值及G+:G-比值常用來表示土壤生態(tài)系統(tǒng)的穩(wěn)定性,比值越大,土壤生態(tài)系統(tǒng)越穩(wěn)定[29]。其中真菌:細菌比值的提高可增加土壤固碳能力,而G+:G-比值高則表明腐殖化程度較高的碳源較多[30]。從真菌:細菌比值來看(圖1a),與CK相比,2016年切碎秸稈還田提高了真菌:細菌比值,而秸稈顆粒還田降低了真菌:細菌比值,其中QS5處理下土壤亞耕層真菌:細菌比值為0.34,顯著高于CK的0.19(<0.05),KL5僅為0.10,顯著低于CK(<0.05);2017年,QS5與KL5真菌:細菌比值均顯著高于CK,其中QS5最高為0.23(<0.05),KL5為0.20(<0.05);2018年各處理真菌:細菌比值均與CK無顯著差異。隨著秸稈還田時間的增加,切碎秸稈還田真菌:細菌比值逐漸降低,而秸稈顆粒還田則逐漸上升,秸稈高量還田更易提高真菌:細菌比值。

    從G+:G-比值來看(圖1b),2016年,QS1的G+:G-比值為1.89,顯著高于CK的1.65(<0.05),而QS5為1.26,顯著低于CK(<0.05);2017年,各處理雖與CK相比無顯著差異,但KL1的G+:G-比值最高為1.89,顯著高于QS5與KL5(<0.05);2018年,雖然秸稈還田處理均較CK提高了G+:G-比值,但無顯著差異。由此可見,G+:G-比值的變化受秸稈形態(tài)與秸稈用量的綜合影響,還田前期高量處理在一定程度上抑制了G+:G-比值的增加,但隨著還田時間的延長,秸稈高量還田的G+:G-比值逐漸提高,這是因為G-對新鮮有機碳的輸入反應更敏感,起初占據(jù)優(yōu)勢地位[31],而到后期有機碳轉(zhuǎn)化為穩(wěn)定碳源,G+開始增長。

    注:不同字母代表相同年份下不同處理間差異顯著(P<0.05)。

    2.1.3 亞表層各菌群PLFA的分布特征

    對亞表層土壤中所檢測出的所有PLFA單體的摩爾數(shù)進行主成分(PCA)分析(圖2、圖3和圖4)。2016年,PC1與PC2分別占據(jù)了總變異的25.23%和24.34%,QS5處理的樣本在第一主成分上明顯與CK區(qū)分開,而QS1處理的樣本在第二主成分上明顯與其他處理的樣本區(qū)分開(圖2a)。其中飽和脂肪酸(14:0、15:0、16:0、17:0、20:0、23:0、i13:0、i14:0、i16:0和i17:0)、環(huán)丙烷脂肪酸(17:0cy ω7c)、甲基支鏈脂肪酸(17:0 (10Me)和16:0(10Me))以及部分單不飽和脂肪酸(16:1 ω7c、16:1 ω5c、17:1 ω8c和18:1 ω9c)和多不飽和脂肪酸(22:5 ω6c、20:5 ω3c、19:3 ω3c和18:2 ω6c)在第一主成分上載荷值較高且為正值(圖2b),而在第一主成分的載荷值的負值區(qū)域僅有部分單不飽和脂肪酸(14:1 ω9c、14:1 ω8c、19:1 ω8c、20:1 ω9c和16:1 ω9c)和多不飽和脂肪酸(21:3 ω3c、19:3 ω6c和16:3 ω6c)分布。從第二主成分上來看,在載荷值為正值部分主要分布飽和脂肪酸及甲基支鏈脂肪酸,而載荷值為負值部分主要分布不飽和脂肪酸和環(huán)丙烷脂肪酸。其中不飽和脂肪酸多為革蘭氏陰性菌和真菌的標記脂肪酸,而飽和脂肪酸、環(huán)丙烷脂肪酸以及甲基支鏈脂肪酸多為革蘭氏陽性菌和放線菌的標記脂肪酸,其中環(huán)丙烷脂肪酸在厭氧環(huán)境下廣泛存在,可作為厭氧環(huán)境指標[32]。

    2017年的PCA分析顯示,PC1與PC2分別占據(jù)了總變異的36.66%和20.66%,在第一主成分上QS5樣本與KL5和QS1樣本與CK和KL1樣本均有較好的區(qū)分,而在第二主成分上,QS1樣本則與其他處理的樣本有較好區(qū)分(圖3a)。從第一主成分來劃分,土壤中多數(shù)PLFA單體的載荷值為正值,僅有少數(shù)多不飽和脂肪酸(21:3 ω3c和18:2 ω6c)以及甲基支鏈脂肪酸(17:1 ω7c(10Me)和16:0(10Me))的載荷值為負值;從第二主成分來劃分,甲基支鏈脂肪酸(18:1 ω7c(10Me)和19:1 ω7c(10Me))與單不飽和脂肪酸(22:1 ω3c和19:1 ω8c)載荷值為正值,而載荷值為負值部分多為飽和脂肪酸(12:0、17:0、a14:0和a16:0)分布(圖3b)。

    2018年的PCA分析顯示,PC1與PC2分別占據(jù)了總變異的28.48%和21.79%,在第一主成分上QS5樣本與CK有明顯區(qū)分,在第二主成分上KL1、KL5和QS1樣本與CK的樣本有較好的區(qū)分(圖4a)。從第一主成分來劃分,多不飽和脂肪酸(24:1 ω3c、20:5 ω3c和20:3 ω6c)以及飽和脂肪酸(12:0、a13:0和i15:1 ω6c)的載荷值為正值,甲基支鏈脂肪酸(18:0(10Me)和16:0(10Me))以及部分飽和脂肪酸和小部分不飽和脂肪酸的載荷值為負值;而飽和脂肪酸與不飽和脂肪酸在第二主成分上均占有相當?shù)恼撦d荷值(圖4b)。

    圖3 2017年不同處理下亞表層土壤微生物群落結(jié)構(gòu)的主成分分析

    圖4 2018年不同處理下亞表層土壤微生物群落結(jié)構(gòu)的主成分分析

    2.2 土壤微生物群落結(jié)構(gòu)與土壤環(huán)境因子的關系

    2.2.1 秸稈還田后土壤環(huán)境因子的變化

    土壤環(huán)境因子的變化是影響土壤微生物群落結(jié)構(gòu)組成的重要因素。表4表明,秸稈還田1 a,即2016年,QS5對土壤容重(SBD)、土壤含水量(SWC)及全氮(TN)的改善效果最顯著,幅度均>20%(<0.05),而KL5對土壤有機碳(SOC)及碳氮比(C/N)的提升最顯著,分別達21.3%(<0.05)與13.2%(<0.05),且秸稈還田對土壤pH有普遍降低作用。秸稈還田2 a,低量還田對各因素的影響程度降低,僅高量還田對各因素仍有顯著影響,對土壤pH開始表現(xiàn)出提升效應,且QS5顯著提高SOC達27.3%(<0.05)。秸稈還田3 a,僅秸稈高量還田對各因子有顯著影響,但影響程度減弱??梢?,秸稈一次性高量還田能夠多年影響土壤環(huán)境因子,對SWC、SOC、TN、C/N的提高具有正效應,對SBD的增加具有負效應,對土壤pH則表現(xiàn)為前期顯著降低后期顯著提高。

    表4 秸稈深埋還田對亞表層土壤主要環(huán)境因子的影響

    2.2.2 土壤微生物群落PLFA與土壤環(huán)境因子的冗余分析

    對20~40 cm亞表層土壤微生物群落PLFA與土壤環(huán)境因子進行冗余分析(圖5),以揭示秸稈還田后引起微生物群落改變的重要原因。RDA分析表明,亞表層土壤環(huán)境因子對土壤微生物群落結(jié)構(gòu)有重要影響,2016年,兩個排序軸解釋了總變異約80%,說明這兩個排序軸在一定程度上能反映土壤環(huán)境因子對微生物群落影響的大部分信息,其中SWC、C/N(<0.05)、TN(<0.05)、pH(<0.05)和SBD(<0.05)是影響微生物群落結(jié)構(gòu)的重要因素,且影響程度依次增大。2017年,兩個排序軸解釋了總變異約70%,其中TN(<0.05)、SBD(<0.01)、C/N(<0.01)、SOC(<0.01)和SWC(<0.01)是影響微生物群落的重要因素,且影響程度依次增大。2018年,2個排序軸對總變異的解釋度不足50%,表明隨著時間的延長,土壤環(huán)境因子對微生物群落結(jié)構(gòu)的影響減弱,只有SOC對微生物群落結(jié)構(gòu)有顯著(<0.05)影響。可見,對微生物群落影響的主要環(huán)境因子隨著還田時間而變化,初始階段微生物受土壤通氣性、酸堿性及C/N影響較大,中期受土壤水分、通氣性及碳源狀況影響較大,后期則受碳源數(shù)量影響大。

    圖5 2016、2017與2018年不同處理下亞表層土壤微生物群落結(jié)構(gòu)與土壤環(huán)境因子的冗余分析

    3 討 論

    3.1 不同秸稈還田方式對亞表層土壤微生物群落結(jié)構(gòu)的影響

    土壤亞表層位于心土層范圍內(nèi),是經(jīng)常被忽視的重要耕作層和根系分布層,研究其微生物菌群結(jié)構(gòu)變化對于指導土壤生物學肥力的提高具有重要意義[7]。秸稈作為一種天然生物質(zhì)資源,施入到土壤中為微生物提供了豐富的碳源和氮源,促進了微生物的生長和繁殖,必然引起微生物群落結(jié)構(gòu)的變化[33],許多研究證實,秸稈翻埋較秸稈覆蓋還田更能顯著提高微生物生物量[34]。本研究中將秸稈深還至亞表層范圍,同時以年為時間尺度研究發(fā)現(xiàn),還田后第一年土壤微生物群落變化幅度最大,土壤磷脂脂肪酸(PLFA)總量顯著提高,相比較而言,對真菌的影響幅度高于細菌,且真菌的變幅僅與秸稈形態(tài)有關,與秸稈用量無關,表現(xiàn)為切碎秸稈處理的土壤真菌摩爾百分數(shù)約為秸稈顆粒的2.5~3.5倍,這表明切碎秸稈在一定時期內(nèi)更利于真菌生長,提高真菌:細菌比值。這主要有兩方面原因,一是切碎秸稈體積較大,施入到土壤中后呈穿插分布狀,形成更大的孔隙結(jié)構(gòu),使氣相占比增加,而土壤大孔隙數(shù)量與真菌的相對豐度顯著相關[35],通氣性好的土壤環(huán)境更利于真菌的生長[36];二是玉米秸稈的表皮組織主要為硅細胞和木栓細胞,還有較多硅化物和蠟質(zhì)[37],主要依賴于真菌進行分解[38],而秸稈顆粒的粉碎化程度高,表皮組織被高度粉碎,供真菌利用的碳源減少[39],導致顆粒處理下真菌占比未有顯著提高。細菌的變幅不及真菌顯著則主要是因為細菌主要出現(xiàn)在中、小孔隙中[40],多為好氧型或兼性厭氧型。秸稈還田2 a后高量還田較低量還田顯著提高了土壤PLFA總量及細菌、真菌、革蘭氏陽性菌(G+)與革蘭氏陰性菌(G-)的摩爾百分數(shù),高量還田優(yōu)勢突顯,這說明成倍增加秸稈用量能夠較長時間維持土壤微生物的活躍狀態(tài),因其延長了腐解期,提高了碳源供應量[14]。然而到秸稈還田第三年,各處理的菌群變化與CK差異不顯著,表明如若沒有外界的反復或長期刺激,微生物群落結(jié)構(gòu)最終回歸于擾動前的穩(wěn)定狀態(tài),許多文獻也表明,土壤微生物群落結(jié)構(gòu)的明顯改變往往發(fā)生在長期施肥或秸稈年年還田的處理中[19,28]。另外,本研究中雖然真菌摩爾百分數(shù)有提高,但仍未在亞表層中形成以真菌為優(yōu)勢種群的微生物群落,細菌仍為優(yōu)勢種群,孫冰潔[35]研究指出,長期免耕和壟作下不能形成真菌優(yōu)勢種群,但長期秋翻下5~10 cm和10~20 cm土層真菌的優(yōu)勢地位明顯,這一方面是與秸稈殘體處于不同的分解階段有關[41],另一方面也受長期試驗影響。秸稈深還田后普遍降低了放線菌PLFA含量,這可能與放線菌適宜較高的pH環(huán)境有關[42]。另外,秸稈還田雖然在一定程度上提高了細菌、真菌的摩爾百分數(shù),但是對微生物群落多樣性卻未見顯著提高,表5根據(jù)微生物各菌群摩爾百分數(shù)計算出香濃指數(shù)(shannon),以表征微生物群落物種豐富度[43],然而秸稈還田的shannon指數(shù)均低于CK,這是因為秸稈還田后,腐解秸稈的菌群相對增多,逐漸占據(jù)優(yōu)勢地位,從而使微生物群落多樣性下降[44]。

    表5 不同秸稈還田處理下亞表層土壤微生物群落物種豐富度(shannon指數(shù))

    3.2 亞表層土壤微生物群落結(jié)構(gòu)變化的原因分析

    秸稈還田在一定程度上改變了微生物的群落分布,在不同秸稈形態(tài)處理中,以切碎秸稈處理與CK差異較大;在不同秸稈用量處理中,以高量還田差異較大。這主要是因為不同秸稈還田處理引起土壤環(huán)境因子的改變,最終促使微生物群落分布特征出現(xiàn)差異[45]。本研究發(fā)現(xiàn),秸稈還田對土壤環(huán)境因子的影響逐年改變,從而使影響微生物群落的主要要素發(fā)生更替。秸稈還田1 a,微生物群落結(jié)構(gòu)的分異要素以秸稈形態(tài)為主,土壤容重(SBD)、pH、全氮(TN)與碳氮比(C/N)是影響微生物群落結(jié)構(gòu)的顯著因素(<0.05),該階段碳源充足,腐解菌的增殖以及活性以通氣狀況和氮源為主要限制因素,秸稈翻埋進入亞表層改變了其孔隙結(jié)構(gòu),使亞表層半?yún)捬醐h(huán)境中的氧氣增多,對于多數(shù)好氧的腐解菌群[46]來說利于其大量繁殖,然而玉米秸稈為高C/N植物殘體,降解時易與土壤爭奪氮素,故除能反映土壤孔隙結(jié)構(gòu)的SBD外,TN也成為影響菌群分布的顯著因素。pH在調(diào)節(jié)微生物群落結(jié)構(gòu)及多樣性上也具有重要作用[47],多數(shù)微生物適宜的pH在6.5~7.5,秸稈還田第1年土壤pH的降低以及可利用碳源的復雜性高均在一定程度上利于真菌生存[38-39]。切碎秸稈的土壤菌群PLFA數(shù)量最高,尤其是真菌數(shù)量的顯著增加,表明切碎秸稈對微生物群落結(jié)構(gòu)的調(diào)控能力更強,這與其對土壤理化因子的顯著影響有關。秸稈還田2 a,微生物群落結(jié)構(gòu)的分異要素以秸稈用量為主,土壤含水量(SWC)、土壤有機碳(SOC)、C/N與SBD是影響微生物群落結(jié)構(gòu)的極顯著因素(<0.01),TN為顯著因素(<0.05),這主要是因為隨著秸稈的腐解,其保水能力逐漸降低[48],水分逐漸成為菌群生長的限制因子;低量秸稈還田在第一年已腐解大部分,微生物可分解利用的新鮮碳源減少,而高量還田為微生物生命活動提供了充足碳源以及較高C/N,促使微生物結(jié)構(gòu)的進一步改變。秸稈還田第三年,土壤微生物群落結(jié)構(gòu)對土壤環(huán)境因子的響應不敏感,這與王曉玥等[43]的研究相似,這是因為3 a來的農(nóng)機具壓實、淋溶作用以及秸稈的腐解共同削弱了秸稈還田對土壤理化因子的影響,各處理的理化因子水平與CK差異縮?。ū?),其中僅SOC仍具有顯著影響(<0.05),表明了一次性秸稈高量還田后SOC的緩慢累積可能是導致后期真菌含量提升的主要原因。需要指出的是,前期試驗證實秸稈顆粒雖然較切碎秸稈有更高的養(yǎng)分釋放速率,但是本研究中發(fā)現(xiàn)其還至亞表層時對SOC、TN及微生物群落結(jié)構(gòu)的調(diào)控能力不及切碎秸稈,這與在我們在黃淮海地區(qū)的研究[49]不一致,究其原因可能有二,一是黑土地區(qū)以雨養(yǎng)農(nóng)業(yè)為主,降水量偏少,秸稈顆粒只有吸水后才能膨散開,才有利于養(yǎng)分釋放;二是當?shù)鼐暥容^高,黑土亞表層冬季出現(xiàn)凍融現(xiàn)象,秸稈顆粒受凍融影響后不易膨散。由此看來,秸稈顆粒在黑土地區(qū)的施用可能更適宜于耕層淺還田,耕層土壤最先吸納利用水分,更易發(fā)揮秸稈顆??纱罅窟€田以及快速釋放養(yǎng)分的優(yōu)勢。

    4 結(jié) 論

    兩種形態(tài)的玉米秸稈深還田可不同程度地提高亞表層土壤磷脂脂肪酸(PLFA)總量、細菌及真菌摩爾百分數(shù),但對放線菌有降低作用。

    切碎秸稈處理較顆粒秸稈處理更利于促進真菌群落的形成。隨還田時間的增加,高量還田的革蘭氏陽性菌(G+):革蘭氏陰性菌(G-)提高更顯著,利于長期維持生態(tài)系統(tǒng)穩(wěn)定性。

    土壤理化因子的改變是引起土壤微生物群落結(jié)構(gòu)分異的重要原因。第1年,土壤容重(SBD)、pH、全氮(TN)與碳氮比(C/N)是引起微生物群落結(jié)構(gòu)分異的顯著性因子(<0.05);第2年,土壤含水量(SWC)、土壤有機碳(SOC)、C/N與SBD為極顯著因子(<0.01);第3年,SOC為顯著性因子(<0.05)。

    總之,QS5處理對黑土亞表層理化因子的改善最為顯著,對真菌群落結(jié)構(gòu)的調(diào)控能力更強,更適宜于東北黑土亞表層土壤生物學肥力的改善。

    [1] Breulmann M, Masyutenko N P, Kogut B M, et al. Short-term bioavailability of carbon in soil organic matter fractions of different particle sizes and densities in grassland ecosystems[J]. Science of Total Environment, 2014, 497: 29-37.

    [2] Cusack D F, Silver W L, Torn M S, et al. Changes in microbial community characteristics and soil organic matter with nitrogen additions in two tropical forests[J]. Ecology, 2011, 92(3): 621-632.

    [3] Garbeva P, Veen J A V, Elsas J D V. Microbial diversity in soil: selection of microbial populations by plant and soil type and implications for soil suppressiveness[J]. Annual Review of Phytopathology, 2004, 42: 243-270.

    [4] Janvier C, Villeneuve F, Alabouvette C, et al. Soil health through soil disease suppression: which strategy from descriptors to indicators[J]. Soil Biology & Biochemistry, 2007, 39(1): 1-23.

    [5] Allison V J, Miller R M, Jastrow J D, et al. Changes in soil microbial community structure in a tallgrass prairie chronosequence[J]. Soil Science Society of America Journal, 2005, 69(5): 1412.

    [6] Buyer J S, Teasdale J R, Roberts D P, et al. Factors affecting soil microbial community structure in tomato cropping systems[J]. Soil Biology & Biochemistry, 2010, 42(5): 831-841.

    [7] Kirkegaard J A, Lilley J M, Howe G N, et al. Impact of subsoil water use on wheat yield[J]. Australia Journal of Agriculture Research, 2007, 58(4): 303-315.

    [8] 劉爽,何文清,嚴昌榮. 不同耕作措施對旱地農(nóng)田土壤物理特性的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2010,28(2):65-69.

    Liu Shuang, He Wenqing, Yan Changrong. Effects of different tillage managements on soil physical properties in dryland[J]. Agricultural Research in the Arid Areas, 2010, 28(2): 65-69. (in Chinese with English abstract)

    [9] Liu Xiaobing, Han Xiaozeng, Song Chunyu, et al. Soil organic carbon dynamics in black soils of China under different agricultural management systems[J]. Communications in Soil Science and Plant Analysis, 2003, 34(7/8): 973-984.

    [10] 龔子同,張之一,張甘霖. 草原土壤:分布、分類與演化[J]. 土壤,2009,41(4):505-511.

    Gong Zitong, Zhang Zhiyi, Zhang Ganlin. Grassland soils: Distribution, classification and evolution[J]. Soils, 2009, 41(4): 505-511. (in Chinese with English abstract)

    [11] 劉武仁,鄭金玉,羅洋,等. 東北黑土區(qū)發(fā)展保護性耕作可行性分析[J]. 吉林農(nóng)業(yè)科學,2008,33(3):3-4,13.

    Liu Wuren, Zheng Jinyu, Luo Yang, et al. Feasible analysis of developing conservation tillage in black soil zone of Northeastern China[J]. Journal of Jilin Agricultural Sciences, 2008, 33(3): 3-4,13. (in Chinese with English abstract)

    [12] 竇森. 玉米秸稈“富集深還”與土壤亞表層培肥[J]. 植物營養(yǎng)與肥料學報,2017,23(6):1670-1675.

    Dou Sen. Improving subsoil fertility through a new technology of continuous in belt and deep incorporation of corn stover[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(6): 1670-1675. (in Chinese with English abstract)

    [13] 韓曉增,鄒文秀,王鳳仙,等. 黑土肥沃耕層構(gòu)建效應[J].應用生態(tài)學報,2009,20(12):2996-3002.

    Han Xiaozeng, Zou Wenxiu, Wang Fengxian, et al. Construction effect of fertile cultivated layer in black soil[J]. Chinese Journal of Applied Ecology, 2009, 20(12): 2996-3002. (in Chinese with English abstract)

    [14] 叢萍,李玉義,王婧,等. 不同用量秸稈深還對提高亞表層土壤肥力質(zhì)量的影響[J]. 植物營養(yǎng)與肥料學報,2019. doi:10.11674/zwyf.19043

    Cong Ping, Li Yuyi, Wang Jing, et al. The feasibility of improving the subsoil fertility quality by incremental-deep straw application[J]. Journal of Plant Nutrition and Fertilizer, 2019. doi:10.11674/zwyf.19043 (in Chinese with English abstract)

    [15] 蔡麗君,張敬濤,蓋志佳,等. 免耕條件下秸稈還田量對土壤酶活性的影響[J]. 土壤通報,2015,46(5):1127-1132.

    Cai Lijun, Zhang Jingtao, Gai Zhijia, et al. Effect of the amount of stalk return to field on soil enzyme activities under no-tillage[J]. Chinese Journal of Soil Science, 2015, 46(5): 1127-1132. (in Chinese with English abstract)

    [16] 陳冬林,易鎮(zhèn)邪,周文新,等. 不同土壤耕作方式下秸稈還田量對晚稻土壤養(yǎng)分與微生物的影響[J]. 環(huán)境科學學報,2010,30(8):1722-1728.

    Chen Donglin, Yi Zenxie, Zhou Wenxin, et al. Effects of straw return on soil nutrients and microorganisms in late rice under different soil tillage systems[J]. Acta Scientiae Circumstantiae, 2010, 30(8): 1722-1728. (in Chinese with English abstract)

    [17] Koullas D P, Christakopoulos P, Kekos D, et al. Correlating the effect of pretreatment on the enzymatic hydrolysis of straw[J]. Biotechnology and Bioengineering, 1992, 39(1): 113-116.

    [18] Wang Jingjing, Zhang Huiwen, Li Xinyu, et al. Effects of tillage and residue incorporation on composition and abundance of microbial communities of a fluvo-aquic soil[J]. European Journal of Soil Biology, 2014, 65: 70-78.

    [19] Zhao Shicheng, Li Kejiang, Zhou Wei, et al. Changes in soil microbial community, enzyme activities and organic matter fractions under long-term straw return in north-central China[J]. Agriculture Ecology Environment, 2016, 216: 82-88.

    [20] Guo Lijin, Zhang Zhisheng, Wang Dandan, et al. Effects of short-term conservation management practices on soil organic carbon fractions and microbial community composition under a rice-wheat rotation system[J]. Biology and Fertility of Soils, 2014, 51(1): 65-75.

    [21] 王婧,張莉,逄煥成,等. 秸稈顆?;€田加速腐解速率提高培肥效果[J]. 農(nóng)業(yè)工程學報,2017,33(6):177-183.

    Wang Jing, Zhang Li, Pang Huancheng, et al. Returning granulated straw for accelerating decomposition rate and improving soil fertility[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(6): 177-183. (in Chinese with English abstract)

    [22] Zhang Li, Wang Jing, Pang Huancheng, et al. Effects of pelletized straw on soil nutrient properties in relation to crop yield[J]. Soil Use and Management, doi: 10.1111/sum.12450.

    [23] 矯麗娜,李志洪,殷程程,等. 高量秸稈不同深度還田對黑土有機質(zhì)組成和酶活性的影響[J]. 土壤學報,2015,52(3):665-672.

    Jiao Lina, Li Zhihong, Yin Chengcheng, et al. Effect of incorporation of crop straw on composition of soil organic matter and enzyme activity in black soil relative to depth and rate of the incorporation[J]. Acta Pedologica Sinica, 2015, 52(3): 665-672. (in Chinese with English abstract)

    [24] Ren Chengjie, Zhao Fazhu, Kang Di, et al. Linkages of C:N:P stoichiometry and bacterial community in soil following afforestation of former farmland[J]. Forest Ecology and Management, 2016, 376: 59-66.

    [25] Frostegard A, Tunlid A, Baath E. Microbial biomass measured as total lipid phosphate in soils of different organic content[J]. Journal of Microbialogical Methods, 1991, 14(3): 151-163.

    [26] Green C T, Scow K M. Analysis of phospholipid fatty acids (PLFA) to characterize microbial communities in aquifers[J]. Hydrogeology Journal, 2000, 8(1): 126-141.

    [27] Ai Chao, Liang Guoqing, Sun Jingwen, et al. The alleviation of acid soil stress in rice by inorganic or organic ameliorants is associated with changes in soil enzyme activity and microbial community composition[J]. Biology and Fertility of Soils, 2015, 51(4): 465-477.

    [28] 張倩. 長期施肥下稻麥輪作體系土壤團聚體碳氮轉(zhuǎn)化特征[D]. 北京:中國農(nóng)業(yè)科學院,2017.

    Zhang Qian. Transformation of Carbon and Nitrogen in Soil Aggregates as Affected by Long-term Fertilization in a Rice-wheat Rotation System[D]. Beijing: Chinese Academy of Agricultural Science, 2017. (in Chinese with English abstract)

    [29] De Vries F T, Hoffland E, Eekeren N V, et al. Fungal/bacterial ratios in grasslands with contrasting nitrogen management[J]. Soil Biology and Biochemistry, 2006, 38(8): 2092-2100.

    [30] Fierer N, Schimel J P, Holden P A, et al. Variations in microbial community composition through two soil depth profiles[J]. Soil Biology and Biochemistry, 2003, 35(1): 167-176.

    [31] 李倩,李曉秀,吳會軍,等. 不同氣候和施肥條件下保護性耕作對農(nóng)田土壤碳氮儲量的影響[J]. 植物營養(yǎng)與肥料學報,2018,24(6):1539-1549.

    Li Qian, Li Xiaoxiu, Wu Huijun, et al. Effects of conservation tillage practices on soil carbon and nitrogen stocks in farmland under different climatic types and fertilization conditions[J]. Journal of Plant Nutrition and Fertilizers, 2018, 24(6): 1539-1549. (in Chinese with English abstract)

    [32] Vestal J R, White D C. Lipid analys is in microbial eology: Quantitative approaches to the study of microbial communities[J]. Biology Seienee, 1989, 39: 535-541.

    [33] Lou Yilai, Liang Wenju, Xu Minggang, et al. Straw coverage alleviates seasonal variability of the topsoil microbial biomass and activity[J]. Catena, 2011, 86(2): 117-120.

    [34] 于寒,梁烜赫,張玉秋,等. 不同秸稈還田方式對玉米根際土壤微生物及酶活性的影響[J]. 農(nóng)業(yè)資源與環(huán)境學報,2015,32(3):305-311.

    Yu Han, Liang Xuanhe, Zhang Yuqiu, et al. Effects of different straw returning modes on the soil microorganism and enzyme activity on corn yield[J]. Journal of Agricultural Resources and Environment, 2015, 32(3): 305-311. (in Chinese with English abstract)

    [35] 孫冰潔. 不同耕作方式下土壤微生物在黑土有機碳固定中的作用研究[D]. 沈陽:中國科學院研究生院(東北地理與農(nóng)業(yè)生態(tài)研究所),2016.

    Sun Bingjie. Effects of Soil Microorganisms on Organic Carbon Sequestration in Black Soil under Different Tillage Practices[D]. Shenyang: Northeast Institute of Geography and Agri-Ecology, Chinese Academy of Sciences, 2016. (in Chinese with English abstract)

    [36] Nunan N, Wu K, Young I M, et al. Spatial distribution of bacterial communities and their relationships with the micro–architecture of soil[J]. FEMS Microbiology Ecology, 2003, 44(2): 203-215.

    [37] 徐廣印,沈勝強,胡建軍,等. 秸稈冷態(tài)壓縮成型微觀結(jié)構(gòu)變化的實驗研究[J]. 太陽能學報,2010,31(3):273-278.

    Xu Guangyin, Shen Shengqiang, Hu Jianjun, et al. Experimental study on the micro structure changes in the process of cold molding with straw[J]. Acta Energiae Solaris Sinica, 2010, 31(3): 273-278. (in Chinese with English abstract)

    [38] Poll C, Marhan S, Ingwersen J, et al. Dynamics of litter carbon turnover and microbial abundance in a rye detritusphere[J]. Soil Biology & Biochemistry, 2008, 40(6): 1306-1321.

    [39] Wei Meng, Hu Guoqing, Wang Hui, et al. 35 years of manure and chemical fertilizer application alters soil microbial community composition in a Fluvo-aquic soil in Northern China[J]. European Journal of Soil Biology, 2017, 82: 27-34.

    [40] Badin A L, Mustafa T, Bertrand C, et al. Microbial communities of urban stormwater sediments: The phylogenetic structure of bacterial communities varies with porosity[J]. FEMS Microbiology Ecology, 2012, 81(2): 324-338.

    [41] Helgason B L, Gregorich E G, Janzen H H, et al. Long–term microbial retention of residue C is site–specific and depends on residue placement[J]. Soil Biology & Biochemisty, 2014, 68: 231-240.

    [42] 張斯童,蘭雪,李哲,等. 微生物降解玉米秸稈的研究進展[J]. 吉林農(nóng)業(yè)大學學報,2016,38(5):517-522.

    Zhang Sitong, Lan Xue, Li Zhe, et al. Research progress of microbial degradation of corn straw[J]. Journal of Jilin Agricultural University, 2016, 38(5): 517-522. (in Chinese with English abstract)

    [43] 王曉玥,蔣瑀霽,隋躍宇,等. 田間條件下小麥和玉米秸稈腐解過程中微生物群落的變化:BIOLOG分析[J]. 土壤學報,2012,49(5):1003-1011.

    Wang Xiaoyue, Jiang Yuji, Sui Yueyu, et al. Changes of microbial communities during decomposition of wheat and maize straw: analysis by biology[J]. Acta Pedologica Sinica, 2012, 49(5): 1003-1011. (in Chinese with English abstract)

    [44] 張紅. 植物秸稈腐解特性及微生物多樣性的響應研究[D].西安:西北農(nóng)林科技大學,2017.

    Zhang Hong. The Response of Plant Straws Decomposition Characteristics and Microbial Diversity[D]. Xi’an: Northeast A&F University, 2017. (in Chinese with English abstract)

    [45] 王芳. 秸稈腐解的微生物群落演替過程研究[D]. 西安:西北農(nóng)林科技大學,2019.

    Wang Fang. The Study on the Microbial Community Succession of Straw Decomposition[D]. Xi’an: Northeast A&F University, 2019. (in Chinese with English abstract)

    [46] 王秀紅,李欣欣,史向遠,等. 好氧堆肥微生物代謝多樣性及其細菌群落結(jié)構(gòu)[J]. 環(huán)境科學研究,2018,31(8):1457-1463.

    Wang Xiuhong, Li Xinxin, Shi Xiangyuan, et al. Microbial metabolism diversity and bacterial flora structure during aerobic composting[J]. Research of Environmental Sciences, 2018, 31(8): 1457-1463. (in Chinese with English abstract)

    [47] Meril? P?ivi, MalmivaaraL?ms? Minna, Spetz Peter, et al. Soil organic matter quality as a link between microbial community structure and vegetation composition along a successional gradient in a boreal forest[J]. Applied Soil Ecology, 2010, 46(2): 259-267.

    [48] Xie Lihua, Liu Mingzhu, Ni Boli, et al. Utilization of wheat straw for the preparation of coated controlled-release fertilizer with the function of water retention[J]. Journal of Agricultural and Food Chemistry, 2012, 60(28): 6921-6928.

    [49] 叢萍,李玉義,高志娟,等. 秸稈顆?;吡窟€田快速提高土壤有機碳含量及小麥玉米產(chǎn)量[J]. 農(nóng)業(yè)工程學報,2019,35(1):148-156.

    Cong Ping, Li Yuyi, Gao Zhijuan, et al. High dosage of pelletized straw returning rapidly improving soil organic carbon content and wheat-maize yield[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 148-156. (in Chinese with English abstract)

    Effects and analysis of straw returning on subsoil microbial community structure in black soil

    Cong Ping1,2, Wang Jing1, Dong Jianxin2, Li Yuyi1, Liu Na1, Pang Huancheng1※

    (1.100081,; 2.266100,)

    The subsoil layer is an important plow layer between topsoil and bottom soil. It is also an important crop root distribution area. However, long-term rotary tillage caused a thick plow pan in the subsoil layer of black soil in Northeast China, which seriously affected the operation of water, fertilizer, gas, and heat in the soil. As is known to all, straw returning is a win-win measure that can not only fertilize the soil but also ensure the utilization of waste resources. Chopped straw is a common way for local farmers to return to the field, while the pelletized straw return is a novel straw return management and is made of straw after highly comminuted and extruded. Our previous research proved that these two forms of deep straw return had a good effect on subsoil physicochemical properties. However, their effects on soil microbial community structure of the subsoil layer were still unknown. To clarified the effect of chopped straw and pelletized straw on soil microbial community structure in subsoil layer of the black soil, then, evaluated the effect of maize straw returning on fertilization of subsoil layer, a one-time deep straw returning experiment was conducted in the black soil of Northeast China from 2016 to 2018. There were five treatments in this experiment, including chopped straw with low dosage (QS1), chopped straw with high dosage (QS5), pelletized straw with low dosage (KL1), pelletized straw with high dosage (KL5), and no straw returning (CK). The contents of soil microbial flora phospholipid fatty acid, soil bulk density, soil water content, pH, soil organic carbon and total nitrogen were measured in the maize harvest season. The results showed that the microbial flora total phospholipid fatty acid content and fungi phospholipid fatty acid content were significantly increased by chopped straw treatments in the first year of straw returning, with a high dosage increasing up to 71.0% and 120.5%. The phospholipid fatty acid content of bacteria, Gram-positive bacteria and Gram-negative bacteria increased more significantly by pelletized treatment, with high dosage treatment up to 41.6%, 29.7% and 26.3%, and the phospholipid fatty acid content of fungi significantly increased in the second year. The phospholipid fatty acid content of each flora, especially the fungal phospholipid fatty acid content of chopped straw with high dosage treatment, was significantly higher than that of pelletized straw with high dosage treatment (21.0%), and only the flora with high dosage straw return had significant changes in the third year. Chopped straw treatment significantly increased the fungi-bacteria ratio in the early stage of straw returning, while low dosage straw returning increased the Gram-positive bacteria: Gram-negative bacteria ratio. With time increasing, the ratio of high dosage straw return increased more significantly, which was conducive to long-term maintenance of ecosystem stability. High dosage of straw returning could significantly change the level of soil physical and chemical factors, which was an important reason for the differentiation of microbial community structure. The significant soil factors changed with the years of straw returning. In the first year, the significant factors were soil bulk density, pH, total nitrogen and C:N ratio (<0.05), in the second year, the significant factors were soil water content, soil organic carbon, C:N ratio and soil bulk density (<0.01), and soil organic carbon was the only significant factor in the third year (<0.05). In conclusion, chopped straw with high dosage treatment had the most obvious difference between microbial community structure and CK, and had a stronger ability to regulate subsoil fungal community. It was more suitable for the improvement of the subsoil fertility of black soil and promoted the resource utilization of straw in Northeast China.

    soils; straw; microbial community structure; deep returning; subsoil

    叢 萍,王 婧,董建新,李玉義,劉 娜,逄煥成. 秸稈還田對黑土亞表層微生物群落結(jié)構(gòu)的影響特征及原因分析[J]. 農(nóng)業(yè)工程學報,2020,36(1):109-118.doi:10.11975/j.issn.1002-6819.2020.01.013 http://www.tcsae.org

    Cong Ping, Wang Jing, Dong Jianxin, Li Yuyi, Liu Na, Pang Huancheng. Effects and analysis of straw returning on subsoil microbial community structure in black soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(1): 109-118. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.01.013 http://www.tcsae.org

    2019-09-02

    2019-12-23

    國家重點研發(fā)計劃2016YFD0300804,旱地保護性耕作及其農(nóng)藝化機械作業(yè)技術(2016~2020)

    叢 萍,助理研究員,博士,從事土壤耕作與培肥技術研究。Email:cpqdjz@126.com

    逄煥成,研究員,博士,從事合理耕層構(gòu)建以及鹽堿地改良利用。Email:panghuancheng@caas.ac.cn

    10.11975/j.issn.1002-6819.2020.01.013

    S154.3

    A

    1002-6819(2020)-01-0109-10

    猜你喜歡
    高量表層群落
    大學生牙齦炎齦上菌斑的微生物群落
    半潛式平臺表層卡套管處理與認識
    海洋石油(2021年3期)2021-11-05 07:43:10
    合成微生物群落在發(fā)酵食品中的應用研究
    長期施用污泥對土壤-蘿卜系統(tǒng)重金屬積累及土壤養(yǎng)分含量的影響
    水體表層沉積物對磷的吸收及釋放研究進展
    山中宰相高量成
    金沙江文藝(2017年2期)2017-11-24 16:57:48
    為相為國理事,退休為民造福
    百家講壇(2017年11期)2017-06-28 20:29:56
    氬弧熔覆原位合成Ti(C,N)-WC增強鎳基表層復合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    春季和夏季巢湖浮游生物群落組成及其動態(tài)分析
    超聲波光整強化40Cr表層顯微硬度研究
    黄色成人免费大全| 丰满的人妻完整版| 成人午夜高清在线视频| 91麻豆av在线| 国产成+人综合+亚洲专区| 悠悠久久av| 日韩欧美一区二区三区在线观看| 日韩av在线大香蕉| 狂野欧美激情性xxxx| 日日夜夜操网爽| 亚洲 欧美 日韩 在线 免费| 亚洲激情在线av| 欧美在线一区亚洲| 日本三级黄在线观看| 搡老熟女国产l中国老女人| 我的老师免费观看完整版| 91麻豆精品激情在线观看国产| 欧美zozozo另类| 桃色一区二区三区在线观看| 国产成+人综合+亚洲专区| 午夜免费成人在线视频| 久久中文字幕一级| 免费看日本二区| 国产亚洲精品av在线| 精品电影一区二区在线| 国产成人av教育| 久久这里只有精品19| 午夜精品一区二区三区免费看| 特级一级黄色大片| 女人高潮潮喷娇喘18禁视频| 搡老妇女老女人老熟妇| 国产成人av教育| 日日爽夜夜爽网站| 好男人电影高清在线观看| 免费在线观看完整版高清| 男女视频在线观看网站免费 | 97人妻精品一区二区三区麻豆| 久久精品国产亚洲av香蕉五月| 听说在线观看完整版免费高清| 真人一进一出gif抽搐免费| 精品午夜福利视频在线观看一区| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 19禁男女啪啪无遮挡网站| 国产三级黄色录像| 亚洲熟妇中文字幕五十中出| 黑人巨大精品欧美一区二区mp4| 黑人巨大精品欧美一区二区mp4| 日本五十路高清| 亚洲国产看品久久| 日韩欧美一区二区三区在线观看| 黄色视频不卡| 制服丝袜大香蕉在线| 欧美一区二区国产精品久久精品 | 成人午夜高清在线视频| 可以在线观看的亚洲视频| 在线看三级毛片| 欧美极品一区二区三区四区| 女人被狂操c到高潮| 成人亚洲精品av一区二区| 久久草成人影院| 亚洲欧美精品综合久久99| 免费看日本二区| 精品国内亚洲2022精品成人| 国产日本99.免费观看| 操出白浆在线播放| 成熟少妇高潮喷水视频| svipshipincom国产片| 九色国产91popny在线| 亚洲午夜理论影院| 国产亚洲精品av在线| 国产精品99久久99久久久不卡| 在线观看免费日韩欧美大片| 亚洲成人久久爱视频| 国产区一区二久久| 麻豆一二三区av精品| 久久久久久免费高清国产稀缺| 日本精品一区二区三区蜜桃| 欧美色视频一区免费| 久久人妻av系列| 丁香六月欧美| 国产精品久久久久久亚洲av鲁大| 亚洲激情在线av| 1024视频免费在线观看| 久久草成人影院| 欧美 亚洲 国产 日韩一| 精品人妻1区二区| 久久精品综合一区二区三区| 久久香蕉激情| 亚洲中文日韩欧美视频| 欧美一区二区精品小视频在线| 男人舔女人下体高潮全视频| 亚洲精品久久成人aⅴ小说| 精品久久蜜臀av无| 国产成人av激情在线播放| 法律面前人人平等表现在哪些方面| 美女扒开内裤让男人捅视频| 亚洲国产欧美网| 麻豆一二三区av精品| 黑人欧美特级aaaaaa片| 两个人看的免费小视频| 免费在线观看日本一区| 免费在线观看日本一区| 一级毛片高清免费大全| 国产探花在线观看一区二区| 窝窝影院91人妻| 亚洲专区国产一区二区| 少妇粗大呻吟视频| 日韩欧美精品v在线| 亚洲乱码一区二区免费版| 丁香欧美五月| 国产精品久久久久久精品电影| 日本一二三区视频观看| 国产成人精品久久二区二区免费| 搡老岳熟女国产| 日日干狠狠操夜夜爽| 色av中文字幕| 精品国产超薄肉色丝袜足j| 天堂动漫精品| xxxwww97欧美| 黑人巨大精品欧美一区二区mp4| 国产成年人精品一区二区| 成年人黄色毛片网站| 亚洲在线自拍视频| 亚洲全国av大片| 久久精品综合一区二区三区| 免费在线观看成人毛片| 在线十欧美十亚洲十日本专区| 黄色女人牲交| 亚洲人成电影免费在线| 亚洲av中文字字幕乱码综合| 99re在线观看精品视频| 麻豆国产97在线/欧美 | 国产99久久九九免费精品| 国产单亲对白刺激| 最近在线观看免费完整版| 青草久久国产| 嫁个100分男人电影在线观看| 久9热在线精品视频| 成人亚洲精品av一区二区| 日韩欧美国产在线观看| 亚洲国产日韩欧美精品在线观看 | 99精品久久久久人妻精品| 久久久久国产一级毛片高清牌| 精品国产乱子伦一区二区三区| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 欧美日韩亚洲综合一区二区三区_| 麻豆成人午夜福利视频| 欧美大码av| 看黄色毛片网站| 国产午夜精品论理片| 精品久久久久久成人av| 一本精品99久久精品77| 人妻丰满熟妇av一区二区三区| 曰老女人黄片| 全区人妻精品视频| 村上凉子中文字幕在线| 女生性感内裤真人,穿戴方法视频| 欧美性猛交黑人性爽| 天堂√8在线中文| 黄色视频,在线免费观看| 亚洲av第一区精品v没综合| 两个人视频免费观看高清| 亚洲精华国产精华精| 在线观看午夜福利视频| 亚洲欧美日韩高清专用| 成人18禁在线播放| 成人精品一区二区免费| 亚洲成人免费电影在线观看| 国产在线观看jvid| 精品乱码久久久久久99久播| 日本a在线网址| 国产区一区二久久| 亚洲在线自拍视频| 美女午夜性视频免费| 久久久久九九精品影院| 窝窝影院91人妻| 伊人久久大香线蕉亚洲五| 国产精品av久久久久免费| 校园春色视频在线观看| 精品电影一区二区在线| a级毛片在线看网站| 两个人看的免费小视频| 在线视频色国产色| 99久久无色码亚洲精品果冻| 脱女人内裤的视频| 18禁观看日本| 这个男人来自地球电影免费观看| 欧美最黄视频在线播放免费| 男人舔女人下体高潮全视频| 日韩有码中文字幕| 亚洲av五月六月丁香网| 久久久久久久午夜电影| 日本三级黄在线观看| 久久香蕉精品热| 天堂√8在线中文| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 一边摸一边抽搐一进一小说| 色av中文字幕| 在线观看舔阴道视频| 精品人妻1区二区| 久久久久久大精品| 国产熟女午夜一区二区三区| 国产又黄又爽又无遮挡在线| 亚洲欧美日韩东京热| 悠悠久久av| 1024香蕉在线观看| 日韩av在线大香蕉| 成人国语在线视频| 欧美大码av| 免费一级毛片在线播放高清视频| 麻豆成人av在线观看| 大型黄色视频在线免费观看| 午夜两性在线视频| 色哟哟哟哟哟哟| 香蕉av资源在线| 国产精品精品国产色婷婷| 午夜精品在线福利| 久9热在线精品视频| 午夜精品一区二区三区免费看| 欧美不卡视频在线免费观看 | 在线免费观看的www视频| 国产成人aa在线观看| 99国产精品一区二区三区| 国产精品,欧美在线| 很黄的视频免费| 国产精品98久久久久久宅男小说| 久久人妻福利社区极品人妻图片| 成人欧美大片| 亚洲成人国产一区在线观看| 亚洲男人的天堂狠狠| 嫩草影视91久久| 日本 欧美在线| 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月| 很黄的视频免费| 亚洲中文字幕日韩| 亚洲人成77777在线视频| 久久伊人香网站| 亚洲精品在线观看二区| 最近在线观看免费完整版| 亚洲色图av天堂| 欧美性猛交黑人性爽| 国产人伦9x9x在线观看| 美女免费视频网站| 夜夜夜夜夜久久久久| 在线a可以看的网站| 久久久精品欧美日韩精品| 精品午夜福利视频在线观看一区| 亚洲av成人不卡在线观看播放网| 又大又爽又粗| www日本在线高清视频| 亚洲,欧美精品.| 两人在一起打扑克的视频| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区激情短视频| 夜夜夜夜夜久久久久| 欧美性长视频在线观看| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 成人三级做爰电影| 亚洲国产欧美网| 欧美最黄视频在线播放免费| 999精品在线视频| 亚洲18禁久久av| 午夜亚洲福利在线播放| 日韩高清综合在线| 黑人欧美特级aaaaaa片| 老熟妇乱子伦视频在线观看| 精华霜和精华液先用哪个| 舔av片在线| 久久久精品国产亚洲av高清涩受| 国产成人av激情在线播放| 中文亚洲av片在线观看爽| 丰满的人妻完整版| 最近最新免费中文字幕在线| 亚洲avbb在线观看| 国产99白浆流出| 国产精品久久久久久人妻精品电影| 看片在线看免费视频| 日韩欧美一区二区三区在线观看| 亚洲avbb在线观看| 国产不卡一卡二| 女警被强在线播放| 老汉色∧v一级毛片| 国产成人精品无人区| 日韩欧美在线二视频| 999久久久精品免费观看国产| 天堂√8在线中文| 真人一进一出gif抽搐免费| 午夜成年电影在线免费观看| 国产免费男女视频| 欧美三级亚洲精品| 日本 av在线| 日韩av在线大香蕉| 色av中文字幕| 亚洲午夜理论影院| 此物有八面人人有两片| 久久香蕉国产精品| 亚洲欧美日韩高清专用| 黄色a级毛片大全视频| 在线观看66精品国产| 国内少妇人妻偷人精品xxx网站 | 久久精品综合一区二区三区| 精品熟女少妇八av免费久了| 欧美乱色亚洲激情| 美女午夜性视频免费| 日本一二三区视频观看| 丝袜人妻中文字幕| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 国产精品一区二区精品视频观看| 亚洲成人中文字幕在线播放| 免费观看人在逋| 日本一二三区视频观看| 久久热在线av| 色哟哟哟哟哟哟| 少妇人妻一区二区三区视频| 亚洲欧美一区二区三区黑人| 国产成人啪精品午夜网站| 国产一级毛片七仙女欲春2| 国产高清激情床上av| 少妇熟女aⅴ在线视频| 日韩精品免费视频一区二区三区| 国产欧美日韩精品亚洲av| 丝袜美腿诱惑在线| 欧美日韩亚洲国产一区二区在线观看| 国产av麻豆久久久久久久| 一级黄色大片毛片| 日韩三级视频一区二区三区| 亚洲成av人片免费观看| 欧美日韩福利视频一区二区| 男女那种视频在线观看| 国产精品影院久久| 久久久国产精品麻豆| 国产成人av教育| 免费一级毛片在线播放高清视频| 日韩免费av在线播放| 久久久精品欧美日韩精品| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 在线视频色国产色| 俺也久久电影网| 国产又黄又爽又无遮挡在线| 日韩大尺度精品在线看网址| 亚洲真实伦在线观看| 欧美黑人欧美精品刺激| 91成年电影在线观看| 久久精品成人免费网站| 91九色精品人成在线观看| 美女免费视频网站| 色综合站精品国产| 国产久久久一区二区三区| 小说图片视频综合网站| av在线天堂中文字幕| 午夜久久久久精精品| 亚洲精品在线观看二区| 婷婷亚洲欧美| 婷婷六月久久综合丁香| 可以在线观看的亚洲视频| 国产午夜精品论理片| 婷婷精品国产亚洲av在线| 国产成人精品无人区| 老司机福利观看| 香蕉久久夜色| 国产成人系列免费观看| av视频在线观看入口| 亚洲一区高清亚洲精品| 免费在线观看视频国产中文字幕亚洲| 欧美午夜高清在线| 91国产中文字幕| 熟女少妇亚洲综合色aaa.| 久久精品91蜜桃| 久久中文看片网| 免费搜索国产男女视频| 制服人妻中文乱码| 国产成人系列免费观看| 成人三级做爰电影| 少妇被粗大的猛进出69影院| 日韩欧美国产在线观看| 99热只有精品国产| 欧美一级a爱片免费观看看 | 国产激情久久老熟女| www日本黄色视频网| 狂野欧美激情性xxxx| 午夜影院日韩av| 狂野欧美白嫩少妇大欣赏| 亚洲精品中文字幕在线视频| 一级作爱视频免费观看| 欧美三级亚洲精品| 成人手机av| 757午夜福利合集在线观看| 51午夜福利影视在线观看| 搞女人的毛片| 精品免费久久久久久久清纯| 欧美日韩精品网址| 在线观看午夜福利视频| 欧美三级亚洲精品| 后天国语完整版免费观看| 999精品在线视频| 韩国av一区二区三区四区| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 亚洲熟妇熟女久久| 麻豆一二三区av精品| av有码第一页| 国产av麻豆久久久久久久| 女生性感内裤真人,穿戴方法视频| 免费观看精品视频网站| 国产亚洲欧美在线一区二区| 一级毛片高清免费大全| 一区福利在线观看| 夜夜夜夜夜久久久久| 一区二区三区高清视频在线| 亚洲熟妇熟女久久| 欧美一级毛片孕妇| 老熟妇仑乱视频hdxx| 99国产综合亚洲精品| 成人手机av| 国产av一区在线观看免费| 999久久久精品免费观看国产| 19禁男女啪啪无遮挡网站| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 成人午夜高清在线视频| 精品国内亚洲2022精品成人| 国产在线观看jvid| 久9热在线精品视频| 色老头精品视频在线观看| 床上黄色一级片| 成人亚洲精品av一区二区| 天天躁夜夜躁狠狠躁躁| 午夜福利在线观看吧| 国产午夜精品论理片| 天堂√8在线中文| 国产精品 欧美亚洲| 国产三级在线视频| 69av精品久久久久久| 成人欧美大片| 欧美午夜高清在线| www国产在线视频色| 亚洲专区中文字幕在线| 俺也久久电影网| 国产精品国产高清国产av| 久久久久久国产a免费观看| 制服人妻中文乱码| 国产主播在线观看一区二区| 国产精品 欧美亚洲| 国产蜜桃级精品一区二区三区| 中文在线观看免费www的网站 | 日日爽夜夜爽网站| 淫秽高清视频在线观看| 亚洲av片天天在线观看| 黑人巨大精品欧美一区二区mp4| 热99re8久久精品国产| 亚洲 欧美 日韩 在线 免费| 一级片免费观看大全| 男女那种视频在线观看| 2021天堂中文幕一二区在线观| 午夜精品在线福利| 在线a可以看的网站| 国产精品九九99| 午夜成年电影在线免费观看| 18禁黄网站禁片免费观看直播| 亚洲免费av在线视频| av天堂在线播放| 两人在一起打扑克的视频| 天天躁夜夜躁狠狠躁躁| 最近在线观看免费完整版| 欧美成人午夜精品| cao死你这个sao货| 一区二区三区国产精品乱码| 亚洲男人天堂网一区| 国产精品av久久久久免费| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 午夜老司机福利片| 在线观看免费视频日本深夜| 国产精品 国内视频| 日本熟妇午夜| 久久人妻福利社区极品人妻图片| 天堂av国产一区二区熟女人妻 | 一级片免费观看大全| 亚洲av五月六月丁香网| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 丝袜美腿诱惑在线| 天天一区二区日本电影三级| 俄罗斯特黄特色一大片| 一本大道久久a久久精品| 精品国产乱码久久久久久男人| 欧美绝顶高潮抽搐喷水| 露出奶头的视频| 亚洲av片天天在线观看| or卡值多少钱| 日日爽夜夜爽网站| 深夜精品福利| 婷婷丁香在线五月| 此物有八面人人有两片| 国产精华一区二区三区| 又黄又爽又免费观看的视频| 嫩草影院精品99| 精品欧美国产一区二区三| 精品久久久久久成人av| 99久久久亚洲精品蜜臀av| 狠狠狠狠99中文字幕| 免费观看精品视频网站| 深夜精品福利| 亚洲人成网站高清观看| 亚洲成人国产一区在线观看| 在线播放国产精品三级| 99国产精品一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 男女视频在线观看网站免费 | 国内精品一区二区在线观看| 精品电影一区二区在线| 国产精品久久久久久人妻精品电影| 99久久99久久久精品蜜桃| 国产av不卡久久| 淫妇啪啪啪对白视频| 五月伊人婷婷丁香| 久9热在线精品视频| 麻豆成人av在线观看| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频| 男人舔女人下体高潮全视频| 国产成人欧美在线观看| 国内少妇人妻偷人精品xxx网站 | 搡老妇女老女人老熟妇| 9191精品国产免费久久| 哪里可以看免费的av片| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 在线看三级毛片| 亚洲精品一区av在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲真实伦在线观看| 桃色一区二区三区在线观看| 精品人妻1区二区| aaaaa片日本免费| 色哟哟哟哟哟哟| 亚洲成av人片在线播放无| 男人舔奶头视频| 两性夫妻黄色片| 手机成人av网站| 51午夜福利影视在线观看| 欧美乱码精品一区二区三区| 精品乱码久久久久久99久播| 精品久久久久久成人av| 亚洲国产看品久久| 最新在线观看一区二区三区| 久久精品人妻少妇| 精品一区二区三区四区五区乱码| 嫩草影视91久久| 国产熟女午夜一区二区三区| 国产一级毛片七仙女欲春2| av福利片在线| 巨乳人妻的诱惑在线观看| 国产成人av激情在线播放| 精品国产乱码久久久久久男人| 又黄又爽又免费观看的视频| 成人av在线播放网站| 中文字幕人成人乱码亚洲影| 高潮久久久久久久久久久不卡| 国产成人精品久久二区二区91| 精品电影一区二区在线| 国产久久久一区二区三区| 成在线人永久免费视频| 麻豆成人av在线观看| 97超级碰碰碰精品色视频在线观看| 久久天堂一区二区三区四区| 又黄又粗又硬又大视频| 午夜精品久久久久久毛片777| 亚洲国产看品久久| 国产成人精品久久二区二区免费| 十八禁人妻一区二区| 国产亚洲精品第一综合不卡| 中文字幕久久专区| 久久亚洲真实| 国产精品影院久久| 国产成人av激情在线播放| 日韩av在线大香蕉| 午夜免费成人在线视频| 国产成人av激情在线播放| 日韩av在线大香蕉| av中文乱码字幕在线| 日本一本二区三区精品| 欧美在线一区亚洲| av天堂在线播放| 亚洲天堂国产精品一区在线| 日本a在线网址| 嫩草影视91久久| 男女那种视频在线观看| 99国产精品一区二区蜜桃av| 神马国产精品三级电影在线观看 | 日韩三级视频一区二区三区| 国产精品98久久久久久宅男小说| 久久久久久亚洲精品国产蜜桃av| 老汉色av国产亚洲站长工具| av福利片在线观看| 久久精品国产99精品国产亚洲性色| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 97超级碰碰碰精品色视频在线观看| 免费高清视频大片| 在线a可以看的网站| 老司机午夜福利在线观看视频| 色综合站精品国产| 动漫黄色视频在线观看| 国产成人啪精品午夜网站| 香蕉丝袜av| 午夜视频精品福利| 亚洲第一电影网av|