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

    不同稻草還田技術(shù)對(duì)煙-稻輪作系統(tǒng)土壤養(yǎng)分、有機(jī)碳及微生物多樣性的影響

    2023-07-18 08:49:32黃錦文李日坤陳志誠(chéng)張汴泓雷涵潘睿欣楊銘榆潘美清唐莉娜
    中國(guó)水稻科學(xué) 2023年4期
    關(guān)鍵詞:輪作根際稻草

    黃錦文 李日坤 陳志誠(chéng) 張汴泓 雷涵 潘睿欣 楊銘榆 潘美清 唐莉娜

    不同稻草還田技術(shù)對(duì)煙-稻輪作系統(tǒng)土壤養(yǎng)分、有機(jī)碳及微生物多樣性的影響

    黃錦文1, 2李日坤1, 2陳志誠(chéng)3張汴泓1, 2雷涵1, 2潘睿欣1, 2楊銘榆1, 2潘美清1, 2唐莉娜3, *

    (1福建農(nóng)林大學(xué)農(nóng)學(xué)院/作物遺傳育種與綜合利用教育部重點(diǎn)實(shí)驗(yàn)室, 福州 350002;2作物生態(tài)與分子生理學(xué)福建省高校重點(diǎn)實(shí)驗(yàn)室, 福州 350002;3福建省煙草專賣局煙草科學(xué)研究所, 福州 350003;*通信聯(lián)系人, email: 704142780@qq.com)

    【目的】探討不同稻草還田技術(shù)措施對(duì)土壤特性的影響及根際微生物修復(fù)和改良土壤的內(nèi)在機(jī)制,探索適合福建煙-稻輪作區(qū)的稻草高效還田技術(shù),以期為煙-稻輪作區(qū)作物優(yōu)質(zhì)高效栽培提供依據(jù)。【方法】以烤煙“云煙87”和晚稻“甬優(yōu)1540”為材料,通過(guò)不同的稻草還田配施技術(shù)處理,進(jìn)行了2年的田間定位試驗(yàn),分析比較了煙-稻輪作系統(tǒng)耕層土壤養(yǎng)分、根際微生物和作物產(chǎn)量的變化規(guī)律以及根際微生物與土壤環(huán)境因子間的內(nèi)在關(guān)系?!窘Y(jié)果】產(chǎn)量分析表明,連續(xù)2年稻草還田后,各處理產(chǎn)量均比對(duì)照顯著增加,其中配施石灰和尿素的T2處理最顯著,其煙草和水稻產(chǎn)量分別比對(duì)照提高25.67%、11.49%;其次是配施腐稈靈和尿素的T3處理,純稻草T1還田處理效果最差;T1處理第1年稻草還田后煙草產(chǎn)量與對(duì)照差異不顯著,第2年煙草和水稻產(chǎn)量分別比對(duì)照提高3.77%和5.90%。耕層土壤養(yǎng)分特性分析表明,經(jīng)連續(xù)2年的稻草還田處理后,耕層土壤養(yǎng)分與有機(jī)碳含量均顯著增加,其中T2處理效果最好,其次是T3處理。T2處理堿解氮、有效磷、速效鉀、有機(jī)碳和可溶性有機(jī)碳含量分別比對(duì)照提高了7.09%、5.97%、3.01%、11.32%和5.47%。根際微生物群落分析表明,與對(duì)照相比,T2、T3處理下碳源代謝相關(guān)的、自毒物質(zhì)降解有關(guān)的節(jié)桿菌屬()和纖維素分解有關(guān)的假散囊菌屬()豐度均顯著上升,而與發(fā)病相關(guān)的腐質(zhì)霉屬()和青霉菌屬()相對(duì)豐度則顯著下降。相關(guān)分析表明,根際有益菌豐度與土壤養(yǎng)分因子多為正相關(guān),而致病菌的豐度則與土壤養(yǎng)分因子負(fù)相關(guān)?!窘Y(jié)論】稻草還田配施石灰或腐稈靈和適量速效氮肥,有利于提高耕層養(yǎng)分含量,促進(jìn)固氮等有益微生物生長(zhǎng),抑制病原菌,為煙-稻輪作提供良好的土壤環(huán)境,進(jìn)而提高作物產(chǎn)量。

    稻草還田;煙-稻輪作;土壤養(yǎng)分;根際微生物;群落多樣性

    福建是中國(guó)優(yōu)質(zhì)煙區(qū)之一,有適宜煙草種植的優(yōu)越生態(tài)條件和悠久種煙史。長(zhǎng)期以煙-稻復(fù)種連作為主。前人研究表明,合理的煙-稻輪作通過(guò)干濕交替引起的土壤氧化和還原作用的交替,有利于調(diào)節(jié)土壤有機(jī)質(zhì)的積累和分解,改善土壤結(jié)構(gòu),從而保持和提高土壤肥力[1]。然而,20世紀(jì)80年代以來(lái),隨著作物高產(chǎn)品種的推廣,化肥用量日益增長(zhǎng),煙區(qū)長(zhǎng)期單一的煙-稻復(fù)種連作耕作制度,導(dǎo)致了土壤物理結(jié)構(gòu)逐漸變差,土壤中微生物區(qū)系平衡被破壞,養(yǎng)分失衡,大量氮肥的投入直接加快了土壤酸化進(jìn)程[2]。土壤質(zhì)量問(wèn)題已經(jīng)成為限制福建煙-稻輪作區(qū)作物生產(chǎn)持續(xù)健康發(fā)展的瓶頸。

    稻草是可被直接利用的可再生纖維素資源,含有大量的有機(jī)碳、纖維素、氨基酸和豐富的礦質(zhì)元素。稻草還田是土壤改良的主要措施之一。前人研究表明,稻草還田對(duì)于提高土壤微生物功能多樣性和群落代謝能力[3-5],改善土壤的理化性狀[6-8],增加土壤肥力[9-11],提高作物產(chǎn)量[12-14]等方面均產(chǎn)生積極影響。但在福建煙-稻輪作區(qū),由于烤煙大田移栽早,導(dǎo)致前茬水稻的稻草還田后溶田時(shí)間縮短,傳統(tǒng)的簡(jiǎn)易稻草還田技術(shù)往往使稻草難以完全腐爛分解,影響到煙株的早生快發(fā),若遇到爛冬天氣,還嚴(yán)重影響耕作,因此,近幾年福建煙-稻輪作區(qū)稻草還田面積在不斷縮減。隨著分子生物學(xué)技術(shù)在土壤環(huán)境領(lǐng)域的應(yīng)用,根際微生物的研究越來(lái)越受到重視。研究發(fā)現(xiàn)根際微生物中一些重要的功能微生物作為主要的共生功能體參與植物根系養(yǎng)分轉(zhuǎn)化,對(duì)調(diào)節(jié)有機(jī)質(zhì)、生產(chǎn)力、土壤碳動(dòng)態(tài)和養(yǎng)分循環(huán)起著關(guān)鍵作用[15]。據(jù)此,本研究通過(guò)連續(xù)2年的不同稻草還田配施處理,結(jié)合根際土壤微生物群落特征分析,考查土壤養(yǎng)分及有機(jī)碳的變化規(guī)律,探討不同稻草還田技術(shù)對(duì)土壤特性的影響及根際微生物修復(fù)和改良土壤的內(nèi)在機(jī)制,探索適合福建煙-稻輪作區(qū)的稻草高效還田技術(shù),以期為煙-稻輪作區(qū)作物優(yōu)質(zhì)高效栽培提供依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)區(qū)域概況

    選用福建主栽烤煙品種云煙87及晚季水稻甬優(yōu)1540為供試材料,于2019年10月―2021年10月在福建省煙草科學(xué)研究所宦溪基地(119°36′86″E,26°17′33″N)進(jìn)行2年連續(xù)跟蹤試驗(yàn)。該區(qū)域?qū)儆趤啛釒ШQ笮约撅L(fēng)氣候,參試土壤為黏壤質(zhì)水稻土,土壤基本理化性狀態(tài)如下:pH值為4.96,有機(jī)質(zhì)含量27.65 g/kg,全氮含量1.68 /kg,全磷含量0.62 g/kg,全鉀含量21.58 g/kg,堿解氮含量96.37 mg/kg,速效磷含量51.10 mg/kg,速效鉀含量190.06 mg/kg。

    1.2 試驗(yàn)設(shè)計(jì)

    田間定位試驗(yàn)共設(shè)置4個(gè)處理:T0,無(wú)稻草還田(對(duì)照);T1,純稻草還田(即傳統(tǒng)稻草還田方式);T2,稻草還田+石灰(1500 kg/hm2)+尿素(75 kg/hm2);T3,稻草還田+腐稈靈(30 kg/hm2)+尿素(75 kg/hm2)。每個(gè)處理設(shè)置3個(gè)重復(fù),共12個(gè)小區(qū),隨機(jī)區(qū)組排列。每個(gè)小區(qū)的面積為120 m2,烤煙種植密度為1.2 m×0.5 m;晚季水稻種植密度0.2 m×0.2 m。小區(qū)之間筑田埂并用塑料薄膜包埋隔開。腐稈靈由廣東高明金葵子科技有限公司生產(chǎn)(有效活菌數(shù)≥5×107/g)。

    稻草還田方法:水稻收獲后,將稻草切成3~5段,均勻撒在田間,T2處理在稻草面上撒施石灰和尿素,T3處理撒施腐稈靈和尿素。然后灌水浸泡3~5 d后,用旋耕機(jī)打田,將稻草打爛,田間保持7~10 cm水層,淹漚35 d左右。整地前15~20 d,開好排水溝讓水分自然落干,耕翻起壟用于種植煙草。2020年2月25日煙苗移栽,7月3日終采,采收結(jié)束后種水稻。晚稻于2020年7月16日插秧,10月23日成熟收割,割后按照試驗(yàn)要求進(jìn)行第2輪的稻草還田處理,2021年重復(fù)上一年的試驗(yàn)??緹熂八揪串?dāng)?shù)貎?yōu)質(zhì)、高產(chǎn)栽培要求進(jìn)行水肥管理,除試驗(yàn)處理外,其他施肥量和管理方式一致。

    1.3 測(cè)定項(xiàng)目與方法

    1.3.1 作物產(chǎn)量測(cè)定

    煙葉成熟分期采收、烘烤,烤后煙葉按小區(qū)分別測(cè)產(chǎn),計(jì)算各小區(qū)產(chǎn)量及上等煙比例。

    水稻成熟期,每個(gè)小區(qū)根據(jù)平均有效穗數(shù)取有代表性水稻5株,考查有效穗數(shù)、每穗總粒數(shù)及實(shí)粒數(shù)、千粒重,計(jì)算結(jié)實(shí)率,各小區(qū)田塊進(jìn)行人工收割測(cè)產(chǎn)。

    1.3.2 土壤養(yǎng)分及有機(jī)碳測(cè)定

    耕層土壤樣本采集,于烤煙旺長(zhǎng)期及水稻成熟收獲后每個(gè)小區(qū)按照S型均勻取10個(gè)點(diǎn),將各點(diǎn)土壤充分混合后,做好標(biāo)記,帶回實(shí)驗(yàn)室風(fēng)干,用于土壤養(yǎng)分、有機(jī)碳及可溶性有機(jī)碳的測(cè)定。

    土壤堿解氮、有效磷、速效鉀含量測(cè)定[16]:堿解氮的測(cè)定采用堿解擴(kuò)散法;速效磷采用碳酸氫鈉浸提-鉬銻抗比色法測(cè)定;速效鉀含量的測(cè)定采用醋酸銨浸提的火焰光度計(jì)法。

    土壤有機(jī)碳(soil organic carbon,SOC)、可溶性有機(jī)碳(dissolved organic carbon,DOC)含量測(cè)定:土壤有機(jī)碳的測(cè)定采用重鉻酸鉀-濃硫酸外加熱法[16];可溶性有機(jī)碳參照占新華方法進(jìn)行測(cè)定[17]。

    DOC有效率[18-19]=(DOC/SOC)×100%。

    1.3.3 煙草根際土壤微生物群落調(diào)查

    樣品采集:于第2年烤煙旺長(zhǎng)期,每個(gè)小區(qū)隨機(jī)取3棵煙株,連根拔起,先將根系周圍土樣輕微抖動(dòng)剝離,留下根系周圍1 mm的土壤,然后在磷酸鹽緩沖液(PBS)中浸泡5 min,用鑷子輕微攪動(dòng)便于根土分離,然后將根系取出,余下土壤溶液置于4℃、10 000離心機(jī)中離心30 s,沉淀土壤作為根際樣品。

    土壤總DNA的提取和高通量測(cè)序分析:土壤總DNA采用土壤總DNA提取試劑盒(PowerSoil DNA Isolation Kit)進(jìn)行提取,采用1%凝膠電泳檢測(cè)后,經(jīng)Nanodrop 2000C測(cè)定總DNA濃度。將檢測(cè)合格的土壤總DNA樣品送至北京奧維森基因科技有限公司進(jìn)行高通量測(cè)序分析,每組樣品包含3個(gè)重復(fù)。細(xì)菌測(cè)序引物338F/806R擴(kuò)增可變區(qū)V3-V4區(qū),真菌使用引物ITS1F-ITS2F擴(kuò)增內(nèi)部轉(zhuǎn)錄間隔區(qū)(ITS1)。采用Illumina Miseq 2500對(duì)純化后的PCR產(chǎn)物進(jìn)行測(cè)序,并進(jìn)行相應(yīng)生物信息學(xué)分析。

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

    測(cè)序數(shù)據(jù)分析:原始數(shù)據(jù)先去除Barcode、Primer序列后使用FLASH軟件進(jìn)行拼接,得到原始Tags數(shù)據(jù)(raw_tags),對(duì)Tags數(shù)據(jù)進(jìn)行過(guò)濾,得到高質(zhì)量Tags數(shù)據(jù)(clean_tags)。使用usearch軟件按照97%相似性序列進(jìn)行操作分類單元(OUT)聚類(不含單序列),并選取OTU的代表性序列,采用RDP Classifier算法和Silva數(shù)據(jù)庫(kù)Release 119對(duì)OTU代表序列進(jìn)行比對(duì)分析,并在各個(gè)水平(Phylum,Class,Order,F(xiàn)amily,Genus,Species)注釋其群落的物種信息。對(duì)各樣本中的OTU豐度進(jìn)行均一化處理,使用Qiime 1.8計(jì)算observed- species,Chao1,Shannon,goods-coverage,PD_whole_tree指數(shù),使用R語(yǔ)言2.15.3繪制稀釋曲線并進(jìn)行Alpha多樣性指數(shù)組間差異分析。使用R語(yǔ)言進(jìn)行PCoA統(tǒng)計(jì)分析和作圖。

    采用Excel 2016和SPSS 20.0軟件進(jìn)行數(shù)據(jù)分析,最小顯著差異法(LSD)比較不同數(shù)據(jù)組間的差異,Origin 2021b SR2、R語(yǔ)言作圖。

    2 結(jié)果與分析

    2.1 煙-稻輪作系統(tǒng)作物產(chǎn)量與產(chǎn)值變化

    2.1.1 烤煙經(jīng)濟(jì)性狀表現(xiàn)

    由表1可知,烤煙2年產(chǎn)量均表現(xiàn)為T2>T3>T1>T0,年度間差異顯著。第1年傳統(tǒng)純稻草還田T1處理產(chǎn)量與對(duì)照T0差異不顯著,其他處理顯著高于T0;第2年各處理產(chǎn)量差異均達(dá)顯著水平。連續(xù)2年的稻草還田后,配施石灰的T2處理和配施腐稈靈的T3處理分別比對(duì)照T0增產(chǎn)25.67%和9.48%,比T1分別增產(chǎn)21.10%和5.50%,T1處理比對(duì)照增產(chǎn)3.77%。各處理產(chǎn)值表現(xiàn)與產(chǎn)量一致,也是T2處理最高,其次是T3處理。均價(jià)和上中等煙比例方面,連續(xù)2年稻草還田后,T2與T3處理顯著高于T1和T0處理,T1處理顯著高于對(duì)照T0。

    表1 煙-稻輪作系統(tǒng)中烤煙經(jīng)濟(jì)性狀的變化

    表中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差;同列數(shù)字后不同字母表示同一年份的處理間在0.05水平差異顯著。*< 0.05;**< 0.01;NS,>0.05。下同。

    Values are means ± standard deviation. Different lowercase letters in each column indicate significant difference between different treatments in the same year at 0.05 level. *< 0.05;**< 0.01;NS,>0.05. The same below.

    2.1.2 水稻產(chǎn)量及其構(gòu)成因素變化

    水稻產(chǎn)量分析表明(表2),各處理產(chǎn)量存在年度間差異,第2年各稻草還田處理產(chǎn)量比第1年均有所增加。2年試驗(yàn)均表現(xiàn)為配施石灰的T2處理及配施腐稈靈的T3處理實(shí)際產(chǎn)量都顯著高于對(duì)照T0和傳統(tǒng)的純稻草還田T1處理,T1處理也顯著高于對(duì)照。連續(xù)2年稻草還田后,T1、T2與T3處理實(shí)際產(chǎn)量分別比對(duì)照T0增產(chǎn)5.90%、11.49%和10.40%。從產(chǎn)量構(gòu)成因素分析,稻草還田主要是增加了單位面積的有效穗數(shù),其次是穗粒數(shù)的增加。

    表2 煙-稻輪作系統(tǒng)中水稻產(chǎn)量及其構(gòu)成因素的變化

    2.2 煙-稻輪作系統(tǒng)耕層土壤養(yǎng)分、有機(jī)碳、可溶性有機(jī)碳含量及其有效率變化

    2.2.1 耕層土壤速效養(yǎng)分含量

    表3表明,連續(xù)2年稻草還田后,各季作物收獲后土壤堿解氮均表現(xiàn)為T2>T3>T1>T0,各處理均顯著高于對(duì)照T0,但T2與T3處理間差異不顯著;有效磷方面,各處理均顯著高于對(duì)照T0,T2處理有效磷增加最多,其次是T1和T3,T1、T2和T3處理分別比對(duì)照T0提高3.68%、5.97%和3.82%。速效鉀方面,經(jīng)2年4茬作物種植后各處理均顯著高于對(duì)照T0,但各稻草還田處理間差異不顯著。

    表3 耕層土壤養(yǎng)分、有機(jī)碳、可溶性有機(jī)碳含量及其有效率變化(2021年)

    DOC有效率=可溶性有機(jī)碳/土壤有機(jī)碳×100%。下同。

    DOC-SOC ratio, Dissolved organic carbon/Soil organic carbon×100%. The same below.

    2.2.2 耕層土壤有機(jī)碳、可溶性有機(jī)碳含量及其有效率變化

    秸稈中含有豐富的碳源,還田后有利于促進(jìn)土壤礦化作用。表3表明,連續(xù)2年稻草還田及4茬作物的種植后,各處理SOC含量均比對(duì)照有顯著提高,T1、T2、T3處理分別比對(duì)照T0提高了4.64%、11.32%、10.48%。

    各稻草還田處理均有助于增加土壤DOC含量(表3),以T2處理效果最顯著,T2處理分別比對(duì)照T0和T1處理提高5.47%、2.08%,其次是T1和T3處理,分別比對(duì)照T0提高3.32%和3.14%。

    DOC有效率即土壤中可溶性有機(jī)碳占土壤總有機(jī)碳的百分比,它反映了土壤有機(jī)碳質(zhì)量[18-19]。表3表明,與各處理DOC含量變化不同的是,各處理DOC有效率均顯著低于對(duì)照T0,T1、T2、T3處理DOC有效率分別比對(duì)照T0下降2.44%、7.31%和7.31%。

    2.3 根際土壤微生物群落變化及其與環(huán)境因子的關(guān)系

    2.3.1 根際土壤微生物群落Alpha多樣性分析

    煙草旺長(zhǎng)期根際微生物群落稀釋曲線(圖1)是采用97%相似度的操作分類單元(OUT)對(duì)所測(cè)定土壤樣本進(jìn)行稀釋,當(dāng)曲線趨向平坦時(shí),說(shuō)明測(cè)序數(shù)據(jù)量足夠大,可以反映樣品中絕大多數(shù)的微生物信息[20]。圖1-A(細(xì)菌)和圖1-B(真菌)均表明,隨著測(cè)序深度的增加,各處理的OTU數(shù)逐漸增長(zhǎng),隨后趨于平緩,說(shuō)明測(cè)序數(shù)據(jù)合理。

    測(cè)序深度指數(shù)(Coverage)為樣本文庫(kù)覆蓋度。表4表明,細(xì)菌和真菌各樣本測(cè)序覆蓋度均大于95%,樣本測(cè)序結(jié)果已經(jīng)包含樣本中絕大部分的物種,可進(jìn)行后續(xù)分析。從細(xì)菌指數(shù)多樣性看,各處理辛普森指數(shù)差異不顯著,豐富度指數(shù)Chao1、多樣性指數(shù)Observed species和香農(nóng)指數(shù)(Shannon index)均是T2處理最高,其次是T1處理;真菌指數(shù)多樣性分析表明(表4),各處理間豐富度指數(shù)Chao1和多樣性指數(shù)Observed species差異不顯著;香農(nóng)指數(shù)均是T1處理最高,其次是T2處理。

    韋恩圖(Venn Graph)反映了處理間特有及共有操作分類單元(OUT)數(shù)量。圖2-A表明,處理T0、T1、T2和T3分別檢測(cè)到2413、2503、2599、2331個(gè)細(xì)菌OTU,其中獨(dú)有OTU分別是76、95、146、55個(gè),共有數(shù)量為2332。圖2-B表明,T0、T1、T2、T3各處理獨(dú)有真菌OTU分別有90、84、94、68個(gè),共有數(shù)量為365。

    A―細(xì)菌稀釋曲線;B―真菌稀釋曲線。

    表4 烤煙旺長(zhǎng)期根際土壤細(xì)菌和真菌的α-多樣性指數(shù)

    A―細(xì)菌維恩圖; B―真菌維恩圖。

    2.3.2 根際土壤微生物群落PCoA分析

    采用主坐標(biāo)分析法(Principal Co-ordinates Analysis,PCoA)分析根際土壤細(xì)菌群落結(jié)構(gòu)(圖3-A),第一主坐標(biāo)(PC1)的貢獻(xiàn)率為32.20%,第二主坐標(biāo)(PC2)的貢獻(xiàn)率為18.95%,T1和T3在負(fù)半軸聚集,T0和T2在正半軸聚集,不同處理間區(qū)分明顯,說(shuō)明各處理下根際土壤細(xì)菌群落存在明顯差異。真菌群落結(jié)構(gòu)PCoA分析可知(圖3-B),第一主坐標(biāo)(PC1)的貢獻(xiàn)率為58.37%,第二主坐標(biāo)(PC2)的貢獻(xiàn)率為21.31%,T0、T3在正半軸聚集,T1、T2則在負(fù)半軸聚集,不同處理間區(qū)分明顯,說(shuō)明不同處理對(duì)真菌群落結(jié)構(gòu)有顯著影響。

    2.3.3 根際微生物群落結(jié)構(gòu)變化

    從圖4-A細(xì)菌門水平上變化可知,一共鑒定到47個(gè)門類,其中相對(duì)豐度占比排前10位的依次是變形菌門(Proteobacteria)、綠彎菌門(Chloroflexi)、酸桿菌門(Acidobacteriota)、放線菌門(Actinobacteriota)、髕骨菌門(Patescibacteria)、擬桿菌門(Bacteroidota)、芽單胞菌門(Gemmatimonadota)、疣微菌門(Verrucomicrobiota)、Myxococcota、WPS-2。T1和T3處理較T0處理變形菌門(Proteobacteria)相對(duì)豐度顯著提高,T2和T3處理較T0處理放線菌門(Actinobacteriota)相對(duì)豐度顯著提高。圖4-B顯示,一共鑒定到520個(gè)細(xì)菌屬類,其中相對(duì)豐度超過(guò)1%的有unidentified、uncultured、uncultured_bacterium羅河桿菌屬()、、伯克霍爾德氏菌屬()、uncultured_、節(jié)桿菌屬()、不動(dòng)桿菌屬()、慢生根瘤菌屬()、鞘脂單胞屬()和。其中,T2、T3處理、節(jié)桿菌屬()的相對(duì)豐度較T0和T1處理顯著提高,T2、T3處理的鞘脂單胞屬(Sphingomonas)、Candidatus_Solibacter屬較T0處理有明顯提高。根際土壤真菌結(jié)構(gòu)分析顯示(圖5-A),所有真菌共包含13個(gè)門,主要是子囊菌門(Ascomycota)、擔(dān)子菌門(Basidiomycota)、霉菌門(Mortierellomycota)、unidentified、壺菌門(Chytridiomycota)、羅茲菌門(Rozellomycota)、球囊菌門(Glomeromycota)、油壺菌門(Olpidiomycota)、梳霉門(Kickxellomycota)等。子囊菌門(Ascomycota)相對(duì)豐度占總量的80%以上,T3與T0相比,子囊菌門(Ascomycota)豐度顯著提高,而T1和T2差異不顯著。與T0相比,擔(dān)子菌門在T1、T2、T3中相對(duì)豐度均顯著上升。真菌群落在屬水平上的差異上(圖5-B),共鑒定到310個(gè)真菌屬類,其中相對(duì)豐度超過(guò)1%的有腐質(zhì)霉屬()、、、葡萄孢屬()、被孢霉屬()、青霉菌屬()、鐮刀菌屬()、和。與對(duì)照T0相比,T1、T2、T3處理腐質(zhì)霉屬()、青霉菌屬()相對(duì)豐度均顯著下降,、鐮刀菌屬()豐度均有上升,其中豐度差異達(dá)到顯著水平。

    A―細(xì)菌; B―真菌。

    A―門水平;B―屬水平。

    2.3.4 根際土壤細(xì)菌、真菌主要菌屬功能分析

    挑選相對(duì)豐度大于1%的已鑒定的根際微生物優(yōu)勢(shì)菌屬進(jìn)行功能分析(表5),細(xì)菌菌屬8個(gè),真菌菌屬6個(gè)。在細(xì)菌方面,相比于T0、T1處理,T2、T3處理顯著提高了具有碳源代謝功能的[21]及具有降解自毒物質(zhì)功能的節(jié)桿菌屬()[22]豐度;相比于對(duì)照T0,T1、T2與T3處理均顯著提高了具有分解有機(jī)質(zhì)功能的[23]豐度;T1與T0處理相比,二者菌屬豐度差異不顯著,T1處理的節(jié)桿菌屬()豐度顯著低于對(duì)照;真菌方面,與T0相比,各稻草還田處理均顯著降低了腐質(zhì)霉屬(,與秸稈腐爛及植物病害發(fā)生相關(guān))[24]的豐度、青霉菌屬(,與植物病害發(fā)生有關(guān))[25]的豐度,顯著提升了具有纖維素分解功能的假散囊菌()[26]豐度。

    A―門水平;B―屬水平。

    2.3.5 根際土壤優(yōu)勢(shì)微生物(屬水平)與土壤養(yǎng)分、有機(jī)碳相關(guān)性分析

    根際土壤優(yōu)勢(shì)微生物(屬水平)與環(huán)境因子的相關(guān)性分析見表6。細(xì)菌方面,與碳源代謝相關(guān)的菌屬與有效磷(AP)極顯著正相關(guān)(0.790**);具有降解自毒物質(zhì)功能的節(jié)桿菌屬()與土壤堿解氮(AN)和速效鉀(AK)呈極顯著正相關(guān),相關(guān)系數(shù)分別為0.856**和0.718**;分解有機(jī)質(zhì)有關(guān)的菌屬與堿解氮(AN)顯著負(fù)相關(guān)(?0.604*)。在真菌方面,除了和菌外,其他優(yōu)勢(shì)菌屬均與環(huán)境因子間呈顯著性相關(guān),其中具有分解纖維素功能的菌屬與堿解氮(AN)、速效磷(AP)及有機(jī)碳(SOC)呈極顯著或顯著正相關(guān);(葡萄孢屬)與堿解氮(AN)呈極顯著負(fù)相關(guān)(?0.702*);與病害發(fā)生有關(guān)的(青霉菌屬)與有效磷(AP)和有機(jī)碳(SOC)呈顯著負(fù)相關(guān);有助于纖維素分解的(鐮刀菌屬)則與堿解氮(AN)、有效磷(AP)、速效鉀(AK)、和有機(jī)碳(SOC)呈極顯著或顯著正相關(guān)。可見,土壤優(yōu)勢(shì)微生物不管是有益菌還是致病菌均與土壤速效養(yǎng)分及有機(jī)碳含量密切相關(guān)。

    表5 根際土壤優(yōu)勢(shì)細(xì)菌、真菌菌屬功能分析

    表6 根際優(yōu)勢(shì)微生物(屬水平)與耕層土壤養(yǎng)分、有機(jī)碳的相關(guān)性

    3 討論

    稻草還田能夠把水稻吸收的大部分元素歸還土壤中,從而有利作物對(duì)養(yǎng)分的吸收,提高作物產(chǎn)量。在本研究中,第1年稻草還田后,T2和T3處理均可以顯著增加煙草和水稻產(chǎn)量,而傳統(tǒng)的純稻草還田T1處理下煙草產(chǎn)量與對(duì)照差異并不顯著;但第2年稻草還田后,各處理烤煙和水稻產(chǎn)量均顯著高于T0,上下兩茬作物均為T2處理產(chǎn)量最高,其次是T3處理。王雪仁等[13]研究發(fā)現(xiàn),稻草直接翻埋還田后反而降低了烤煙產(chǎn)量但可減輕田間青枯病發(fā)病率,而稻草堆漚發(fā)酵后還田可顯著增加烤煙產(chǎn)量。黃平娜等[14]研究表明,稻草翻壓還田效果優(yōu)于稻草覆蓋還田,對(duì)烤煙產(chǎn)量的影響以第2年、第3的效果優(yōu)于第1年??梢?,作物生長(zhǎng)對(duì)稻草還田的響應(yīng)因還田處理技術(shù)及還田時(shí)間長(zhǎng)短而異。在本研究中,T2與T3處理效果優(yōu)于傳統(tǒng)純稻草還田T1處理,這是由于T2和T3處理分別加入了石灰和腐稈靈有利于加速稻草的腐爛分解,且這兩個(gè)處理在稻草還田時(shí)均配施了氮素以補(bǔ)充微生物對(duì)氮素的需求,因而有利微生物生長(zhǎng),增加土壤養(yǎng)分,進(jìn)而促進(jìn)作物生長(zhǎng)[5]。而純稻草還田的T1處理由于稻草腐解緩慢,再加上煙草生長(zhǎng)過(guò)程中大量分解稻草的微生物需要消耗土壤氮源,容易出現(xiàn)與作物爭(zhēng)氮的現(xiàn)象,因而短期內(nèi)反而不利煙草生長(zhǎng)。

    稻草還田是改良土壤質(zhì)地、提高綜合生產(chǎn)能力的重要措施。經(jīng)連續(xù)2年的稻草還田及4茬作物種植后,耕層土壤養(yǎng)分測(cè)定結(jié)果表明,速效養(yǎng)分供應(yīng)方面均顯著提升,這與前人研究一致[6];T2、T3處理耕層土壤各速效養(yǎng)分含量均顯著高于對(duì)照,其中T2處理效果最佳,其堿解氮、有效磷和速效鉀分別比對(duì)照T0提高了7.09%、5.97%和3.01%;其次是T3處理。這可能是由于新鮮秸稈還田后,施用石灰可活化土壤中礦化態(tài)氮素、促進(jìn)秸稈有機(jī)物分解[12],而施用腐稈靈,加速秸稈腐熟后形成的腐殖質(zhì)具有溶磷作用,則有助于堿解氮和有效磷等養(yǎng)分積累[8]。同時(shí),T2與T3處理在稻草還田時(shí)添加適量速效氮肥,有利稻草還田前期土壤分解菌大量快速繁殖,形成優(yōu)勢(shì)菌群,進(jìn)而加快稻草的分解。

    土壤有機(jī)碳(SOC)是碳庫(kù)中最活躍的組成部分,是農(nóng)作物生長(zhǎng)發(fā)育的根基,也是土壤有機(jī)質(zhì)和各種養(yǎng)分良性循環(huán)的基礎(chǔ)[30]。水稻秸稈有機(jī)碳含量豐富,秸稈還田是影響有機(jī)碳轉(zhuǎn)化和積累的重要農(nóng)田管理措施之一[31]。本研究發(fā)現(xiàn),經(jīng)過(guò)2年稻草還田及4茬作物的種植后T1、T2、T3處理SOC分別比對(duì)照T0提高了4.64%、11.32%、10.48%??扇苄杂袡C(jī)碳(DOC)作為碳循環(huán)中的重要組成部分,在 SOC形成中起著重要作用。本研究中,與對(duì)照相比,T1、T2與T3處理均顯著提高土壤DOC含量,表現(xiàn)為T2>T3>T1>T0。這是因?yàn)榈静葸€田后大量的營(yíng)養(yǎng)被土壤中的微生物吸收利用,尤其是加石灰或腐稈靈的處理其微生物生長(zhǎng)更快,進(jìn)而加快了秸稈的腐解,提升了土壤DOC的數(shù)量[32]。DOC有效率即土壤中DOC占SOC的百分比[18],主要反映土壤有機(jī)碳的穩(wěn)定性及損失情況,分配比越高土壤有機(jī)碳的活性越大,穩(wěn)定性也越差[33]。該比例雖然很低,但其作為微生物生長(zhǎng)的主要能源,在碳素循環(huán)方面具有重要作用[34]。在本研究中,各處理DOC有效率變化與其DOC含量變化呈現(xiàn)相反趨勢(shì),即各處理該比值均顯著低于對(duì)照T0,表現(xiàn)為T0>T1>T3>T2趨勢(shì)。這可能是因?yàn)榈静葸€田能刺激微生物生長(zhǎng),可溶性有機(jī)物因具有較高的生物降解性,很容易被微生物作為首選基質(zhì)底物而消耗掉[35],而加石灰或腐稈靈的稻草還田,由于土壤中微生物數(shù)量過(guò)于龐大,因此導(dǎo)致大量的可溶性有機(jī)物被消耗[36],進(jìn)而導(dǎo)致其DOC有效率降低。

    植物與環(huán)境之間的交互作用關(guān)鍵在于根際,根際微生物聚集在根系周圍,將有機(jī)物轉(zhuǎn)變?yōu)闊o(wú)機(jī)物,為植物提供養(yǎng)分,因此,根際微生物對(duì)土壤環(huán)境的調(diào)控,已成為當(dāng)前農(nóng)田土壤研究的熱點(diǎn)問(wèn)題[37, 38]。本研究表明,加石灰的T2處理不管是細(xì)菌還是真菌群落其多樣性均最突出。從群落構(gòu)成分析,在細(xì)菌門水平上一共鑒定到47個(gè)門類,其中變形菌門(Proteobacteria)相對(duì)豐度最高,為優(yōu)勢(shì)菌門,是細(xì)菌中最大的一門,包括很多病原菌,也包括很多可以進(jìn)行固氮的細(xì)菌;T2和T3處理較T0處理放線菌門(Actinobacteriota)相對(duì)豐度顯著提高,放線菌門主要促進(jìn)土壤動(dòng)植物殘?bào)w腐爛,同時(shí)對(duì)自然界的氮素循環(huán)也起著一定的作用[39];而在真菌門水平上,共包含13個(gè)門,其中子囊菌門(Ascomycota)相對(duì)豐度占總量的80%以上,與T0相比,擔(dān)子菌門(Basidiomycota)在T1、T2、T3中相對(duì)豐度均顯著上升。Carniel等[40]研究發(fā)現(xiàn),擔(dān)子菌門(Basidiomycota)和子囊菌門(Ascomycota)通過(guò)分泌漆酶和木質(zhì)素過(guò)氧化物酶降解秸稈中的木質(zhì)素。在細(xì)菌屬水平上,一共鑒定到520個(gè)屬類,其中,豐度占比前10位的細(xì)菌屬中,T1和T3處理的羅河桿菌屬()的相對(duì)豐度較T0顯著提高,T2和T3處理的節(jié)桿菌屬()較T0顯著提高。有研究表明節(jié)桿菌屬()具有降解連作障礙產(chǎn)生的自毒物質(zhì)作用[23]。在真菌屬水平上共鑒定到310個(gè)真菌屬類,在占比前10位的真菌屬中,與T0處理相比,T1、T2、T3處理腐質(zhì)霉屬()、青霉菌屬()相對(duì)豐度均顯著下降,、鐮刀菌屬()相對(duì)豐度顯著上升。主要菌屬功能分析發(fā)現(xiàn)(表5),和鐮刀菌屬()具有較強(qiáng)的纖維素分解能力[26];而腐質(zhì)霉屬()是腐生菌,同時(shí)也與多種疾病的發(fā)生有關(guān)[24];青霉菌屬()與多種病害的產(chǎn)生有關(guān)[25]。由此可見,T2與T3處理一方面促進(jìn)了根際土壤分解秸稈和固氮有益菌群的增加,另一方面還抑制了病害相關(guān)的致病菌生長(zhǎng)。進(jìn)一步的相關(guān)分析發(fā)現(xiàn),碳源代謝有關(guān)的[21]具有降解降解自毒物質(zhì)功能的(節(jié)桿菌屬)[22]分解有機(jī)質(zhì)功能的[23]分解纖維素功能的和鐮刀菌屬[26]與土壤AN、AP、AK和SOC等養(yǎng)分因子呈顯著或極顯著正相關(guān)(表6),而與病害發(fā)生有關(guān)的腐質(zhì)霉屬[24]和青霉菌屬[25]與土壤AK和SOC等因子呈顯著或極顯著負(fù)相關(guān)(表6)??梢姡H有益菌的增加與大部分的土壤養(yǎng)分因子呈正相關(guān)關(guān)系,而致病性有害菌屬則與土壤養(yǎng)分呈負(fù)相關(guān)關(guān)系。

    [1] 熊云明, 黃國(guó)勤, 王淑彬, 劉隆旺. 稻田輪作對(duì)土壤理化性狀和作物產(chǎn)量的影響[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2004(4): 42-45.

    Xiong Y M, Huang G Q, Wang S B, Liu L W. Effect of crop rotation in paddy field on soil physical and chemical characteristics and crop yield[J]., 2004(4): 42-45. (in Chinese with English abstract)

    [2] 向鵬華, 單雪華, 黃銀章, 郭維, 龍世平. 煙-稻復(fù)種連作年限對(duì)土壤理化性狀及煙葉產(chǎn)量與品質(zhì)的影響[J]. 中國(guó)土壤與肥料, 2016(5): 105-109.

    Xiang P H, Shan X H, Huang Y Z, Guo W, Long S P. Effects of tobacco-rice continuous cropping years on soil physicochemical properties and tobacco yield and quality[J]., 2016 (5): 105-109. (in Chinese with English abstract)

    [3] 周文新, 陳冬林, 卜毓堅(jiān), 屠乃美. 稻草還田對(duì)土壤微生物群落功能多樣性的影響[J]. 環(huán)境科學(xué)學(xué)報(bào), 2008(2): 326-330.

    Zhou W X, Chen D L, Bu Y J, Tu N M. Effects of rice-straw returning to the field on the metabolic diversity of soil microbial communities[J]., 2008(2): 326-330. (in Chinese with English abstract)

    [4] Yan S Y, Song J M, Fan J S, Yan C, Dong S K, Ma C M, Gong Z P. Changes in soil organic carbon fractions and microbial community under rice straw return in Northeast China[J]., 2020, 22, e00962: 1-12.

    [5] 倪國(guó)榮, 涂國(guó)全, 魏賽金, 吳建富, 潘曉華. 稻草還田配施催腐菌劑對(duì)晚稻根際土壤微生物與酶活性及產(chǎn)量的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2012, 31(1): 149-154.

    Ni G R, Tu G Q, Wei S J, Wu J F, Shi Q H. Effects of straw-returning using agent on microbe and enzyme activity in rhizosphere soils and yield of late rice[J]., 2012, 31(1): 149-154. (in Chinese with English abstract)

    [6] 黃平娜, 秦道珠, 龍懷玉, 張認(rèn)連, 雷秋良, 楊全柳. 稻草還田對(duì)煙田速效養(yǎng)分變化及煙葉產(chǎn)量品質(zhì)的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2008, 24(12): 294-297.

    Huang P N, Qin D Z, Long H Y, Zhang R L, Lei Q L, Yang Q L. Straw also available nutrient-changes in the tobacco fields and tobacco yield and quality[J]., 2008, 24(12): 294-297. (in Chinese with English abstract)

    [7] 陸文龍, 趙標(biāo), 五毛毛. 施用不同方式處理的秸稈對(duì)土壤磷形態(tài)分布的影響[J]. 江蘇農(nóng)業(yè)科學(xué), 2018, 46(2): 232-234.

    Lu W L, Zhao B, Wu M M.Effects of different treatments of straw on soil phosphorus speciation distribution[J]., 2018, 46(2): 232-234. (in Chinese with English abstract)

    [8] 易鎮(zhèn)邪, 符呈祥, 褚百鳳, 付紅, 陳冬林, 周文新, 屠乃美. 快腐劑處理還田早稻秸稈對(duì)晚季土壤化學(xué)與生物學(xué)特性的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2012, 28(27): 94-98.

    Yi Z X, Fu C X, Chu B F, Fu H, Chen D L, Zhou W X, Tu N M. Effect of treating returned early rice straw by decomposition accelerant on soil chemical and biological characteristics in late season[J]., 2012, 28(27): 94-98. (in Chinese with English abstract)

    [9] 吳海勇, 李明德, 劉瓊峰, 吳小丹. 稻草不同途徑還田對(duì)土壤結(jié)構(gòu)及有機(jī)質(zhì)特征的影響[J]. 土壤通報(bào), 2012,43(4):836-841.

    Wu H Y, Li M D, Liu Q F, Wu X D. Effects of different modes of straw returning on soil structure and character of soil organic matter[J]., 2012, 43(4): 836-841. (in Chinese with English abstract)

    [10] 李明德, 肖漢乾, 湯海濤, 張一揚(yáng), 湯睿, 靳志麗, 張淑貞, 張四偉. 稻草還田對(duì)煙田土壤性狀和煙草產(chǎn)量及品質(zhì)的影響[J]. 中國(guó)土壤與肥料, 2006(6): 41-44.

    Li M D, Xiao H Q, Tang H T, Zhang Y Y, Tang R, Jin Z L, Zhang S Z, Zhang S W. Effect of incorporation straw into field on soil properties and tobacco yield and quality[J]., 2006(6): 41-44. (in Chinese with English abstract)

    [11] Zhao S C, He P, Qiu S J, Jia L L, Liu M C, Jin J Y, Johnston A M. Long-term effects of potassium fertilization and straw return on soil potassium levels and crop yields in north-central China[J]., 2014, 169: 116-122.

    [12] 廖萍, 劉磊, 何宇軒, 唐剛, 張俊, 曾勇軍, 吳自明, 黃山. 施石灰和秸稈還田對(duì)雙季稻產(chǎn)量和氮素吸收的互作效應(yīng)[J]. 作物學(xué)報(bào), 2020, 46(1): 84-92.

    Liao P, Liu L, He Y X, Tang G, Zhang J, Zeng Y J, Wu Z M, Huang S. Interactive effects of liming and straw incorporation on yield and nitrogen uptake in a double rice cropping system[J]., 2020, 46(1): 84-92. (in Chinese with English abstract)

    [13] 王雪仁, 張瀛, 張珊珊, 林建麒. 不同方式稻草還田對(duì)烤煙產(chǎn)質(zhì)量的影響[J]. 安徽農(nóng)業(yè)科學(xué), 2017, 45(27): 41-44.

    Wang X R, Zhang Y, Zhang S S, Lin J L. Effects of different ways of straw returning to soil on yield and quality of flue-cured tobacco[J]., 2017, 45(27): 41-44. (in Chinese with English abstract)

    [14] 黃平娜, 秦道珠, 龍懷玉, 劉淑軍,黃晶. 連續(xù)3年稻草還田對(duì)烤煙產(chǎn)量品質(zhì)及后茬晚稻產(chǎn)量影響[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2010, 26(11): 194-199.

    Huang P N, Qin D Z, Long H Y, Liu S J, Huang J. Effects of rice straw for three consecutive years to the yields and qualities of flue-cured tobacco and stalk return field late rice[J]., 2010, 26(11): 194-199. (in Chinese with English abstract)

    [15] 劉王鎖, 李海泉, 何毅, 黃業(yè)蕓, 邱開陽(yáng), 謝應(yīng)忠. 根際微生物對(duì)植物與土壤交互調(diào)控的研究進(jìn)展[J]. 中國(guó)土壤與肥料, 2021(5): 318-327.

    Liu W S, Li H Q, He Y, Huang Y Y, Qiu K Y, Xie Y Z. Research progress on the interaction regulation of plant and soil by rhizosphere microbiome[J]., 2021(5): 318-327. (in Chinese with English abstract)

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

    Bao S D. Soil and Agricultural Chemistry Analysis[M]. Beijing: China Agriculture Press, 2010. (in Chinese)

    [17] 占新華, 周立祥. 土壤溶液和水體中水溶性有機(jī)碳的比色測(cè)定[J]. 中國(guó)環(huán)境科學(xué), 2002(5):433-437.

    Zhan X H, Zhou L X. Colorimetric determination of dissolved organic carbon in soil solution and water environment[J]. China Environmental Science, 2002(5):433-437. (in Chinese with English abstract)

    [18] 周吉祥, 張賀, 楊靜, 李桂花, 張建峰. 連續(xù)施用土壤改良劑對(duì)沙質(zhì)潮土肥力及活性有機(jī)碳組分的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(16): 3307-3318.

    Zhou J X, Zhang H, Yang J, Li G H, Zhang J F. Effects of continuous application of soil amendments on fluvo-aquic soil fertility and active organic carbon components[J]., 2020, 53(16): 3307-3318. (in Chinese with English abstract)

    [19] Li H, Yang S, Semenov M V, Yao F, Ye J, Bu R C, Ma R A, Lin J J, Kurganova I, Wang X G, Deng Y, Kravchenko I, Jiang Y, Kuzyakov Y. Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community[J]., 2021, 27(12): 1-17. doi: 10.1111/gcb.15593.

    [20] Wang Y, Sheng H F, He Y, Wu J Y, Jiang Y X, Tam N F, Zhou H W. Comparison of the levels of bacterial diversity in freshwater, intertidal wetland, and marine sediments by using millions of illumina tags[J]., 2012, 78(23): 8264-8271.

    [21] Yang S D, Xiao J, Liang T, He W Z, Tan H W. Response of soil biological properties and bacterial diversity to different levels of nitrogen application in sugarcane fields[J]., 2021, 11(1): 172

    [22] 李敏, 王桂蓮, 馬璐, 張琇. 節(jié)桿菌降解阿魏酸的效能[J]. 微生物學(xué)通報(bào), 2021, 48(5): 1550-1559.

    Li M, Wang G L, Ma L, Zhang X. Degradation of ferulic acid bysp.[J]., 2021, 48(5): 1550-1559. (in Chinese with English abstract)

    [23] 黎妍妍, 馮吉, 王林, 付裕、李錫宏. 萬(wàn)壽菊-煙草輪作調(diào)理植煙土壤細(xì)菌群落結(jié)構(gòu)的作用[J]. 中國(guó)煙草科學(xué), 2021, 42(1): 14-19.

    Li Y Y, Feng J, Wang L, Fu Y, Li X H. Effects of marigold-tobacco rotation on bacterial community in tobacco rhizosphere soil[J]., 2021, 42(1): 14-19. (in Chinese with English abstract)

    [24] 楊金燕, 姜于蘭, 楊亞曦, 曾琛, 王肸芃. 腐質(zhì)霉屬真菌分類的研究進(jìn)展[J]. 貴州農(nóng)業(yè)科學(xué), 2015, 43(8): 126-130.

    Yang J Y, Jiang Y L, Yang Y X, Zeng C, Wang X P. Advances in taxonomy of Humicola genera[J]., 2015, 43(8): 126-130. (in Chinese with English abstract)

    [25] 戴蓬博, 張榮, 孫廣宇. 中國(guó)蘋果病害病原菌物名錄[J]. 菌物學(xué)報(bào), 2021,40(4): 936-964.

    Dai P B, Zhang R, Sun G Y. A checklist of pathogenic fungi on apple in China[J]., 40(4): 936-964. (in Chinese with English abstract)

    [26] Doveri F, Sarrocco S, Vannacci G. Studies on three rare coprophilous plectomycetes from Italy[J]., 2013, 124: 279-300.

    [27] Heuvel R N V D, Biezen E V D, Jetten M S M, Hefting M M, Kartal B. Denitrification at pH 4 by a soil-derived Rhodanobacter-dominated community[J]., 2010, 12(12): 3264-3271.

    [28] 張珂飛,鐘永嘉, 孫麗莉, 廖紅. 植物有益伯克霍爾德氏菌的研究進(jìn)展及其在農(nóng)業(yè)中的應(yīng)用[J]. 微生物學(xué)報(bào), 2021, 61(8): 2205-2218.

    Zhang K F, Zhong Y J, Sun L L, Liao H. Plant-associated beneficial[J]., 2021, 61(8): 2205-2218. (in Chinese with English abstract)

    [29] 張科, 李臻, 鄭瑤, 麻紅星, 劉夢(mèng)含, 丁慧杰, 王瑜, 劉麗, 夏西超. 河南葉縣巖鹽可培養(yǎng)中度嗜鹽菌的多樣性[J]. 微生物學(xué)通報(bào), 2020, 47(12): 3987-3997.

    Zhang K, Li Z, Zheng Y, Ma H X, Liu M H, Ding H J, Wang Y, Liu L, Xia X C. Biodiversity of culturable moderate halophilic bacteria of rock salt in Yexian county, Henan Province[J]., 2020, 47(12): 3987-3997. (in Chinese with English abstract)

    [30] Schlesinger W H. Evidence from chrono sequence studies for a low carbon-storage potential of soils[J]., 1990, 348(6298): 232-234.

    [31] 趙惠麗, 董金琎, 師江瀾, 徐苗, 田霄鴻. 秸稈還田模式對(duì)小麥-玉米輪作體系土壤有機(jī)碳固存的影響[J]. 土壤學(xué)報(bào), 2021, 58(1): 213-224.

    Zhao H L, Dong J J, Shi J L, Shi J L, Xu M, Tian X H. Effect of straw returning mode on soil organic carbon sequestration[J]., 2021, 58(1): 213-224. (in Chinese with English abstract)

    [32] 王保君, 王偉, 胡乃娟, 顧澤海, 陳兵林, 張政文, 徐蔣來(lái), 朱利群. 麥秸還田下水氮管理對(duì)稻田土壤養(yǎng)分、酶活性及碳庫(kù)的短期影響[J]. 核農(nóng)學(xué)報(bào), 2016, 30(5): 957-964.

    Wang B J, Wang W, Hu N J, Gu Z H, Chen B L, Zhang Z W, Xu J L, Zhu L Q. Short-term effect of different water and nitrogen managements on paddy soil nutrient, enzyme activity and carbon pool under wheat straw-returning fields[J]., 2016, 30(5): 957-964. (in Chinese with English abstract)

    [33] Jandl R, Sollins P. Water-extractable soil carbon in relation to the belowground carbon cycle[J]., 1997, 25: 196-201.

    [34] Chantigny M H. Dissolved and water-extractable organic matter in soils: A review on the influence of land use and management practices[J]., 2003, 113: 357-380.

    [35] Chen H L, Zhou J M, Xiao B H. Characterization of dissolved organic matter derived from rice straw at different stages of decay[J]., 2010, 10(5): 915-922.

    [36] Kalbitz K, Solinger S, Park J H, Michalzik B, Matzner E. Controls on the dynamics of dissolved organic matter in soils: A review[J]., 2000, 165(4): 277-304.

    [37] Garcia J, Kao-Kniffin J. Microbial group dynamics in plant rhizospheres and their implications on nutrient cycling[J]., 2023, 14: 1186322. doi: 10.3389/ fmicb.2018.01516.

    [38] Li X Z, Rui J P, Xiong J B, Li J B, He Z L, Zhou J Z, Yannarell A C, Mackie R I. Functional potential of soil microbial communities in the maize rhizosphere[J]., 2014, 9(11): e112609.

    [39] 吳晉元, 職曉陽(yáng), 李巖, 關(guān)統(tǒng)偉, 唐蜀昆, 徐麗華, 李文均. 云南江城和黑井鹽礦沉積物未培養(yǎng)放線菌多樣性比較[J]. 微生物學(xué)報(bào), 2008(10): 1550-1555.

    Wu J Y, Zhi X Y, Li Y, Guan T W, Tang S K, Xu L H, Li W J. Comparison of actinobacterial diversity in Jiangcheng and Heijing saline mines in Yunnan by using culture-independent approach[J]., 2008(10): 1550-1555. (in Chinese with English abstract)

    [40] Carniel F C, Fortuna L, Zanelli D, Tretiach M. Graphene environmental biodegradation: Wood degrading and saprotrophic fungi oxidize few-layer graphene[J]., 2021: 125553.

    Effects of Straw Returning Techniques on Soil Nutrients, Organic Carbon and Microbial Diversity in Tobacco-rice Rotation System

    HUANG Jinwen1, 2, LI Rikun1, 2, CHEN Zhicheng3, ZHANG Bianhong1, 2, LEI Han1, 2, PAN Ruixin1, 2, YANG Mingyu1, 2, PAN Meiqing1, 2, TANG Lina3, *

    (Key Laboratory for Genetic Breeding and Multiple Utilization of Crops, Ministry of Education/College of Agriculture, Fujian Agriculture and Forestry University, Fuzhou 350002, China; Key Laboratory of Crop Ecology and Molecular Physiology, Fuzhou 350002, China; Fujian Provincial Tobacco Agriculture Research Institute, Fuzhou 350003, China; Corresponding author, email: 704142780@qq.com)

    【Objective】It is of great significance to elucidate the effects of different rice-straw returning techniques on soil properties and the underlying mechanisms of rhizosphere microbial remediation and soil improvement, and to find an efficient straw returning technology suitable for the tobacco-rice rotation area in Fujian, so as to lay a basis for high-quality and efficient crop cultivation.【Method】Using flue-cured tobaccoYunyan 87and late riceYongyou 1540as materials, a 2-year field experimentwas carried out to reveal changes in soil nutrient properties and microbial communities in tobacco-rice rotation system and their relationships under various straw returning modes.【Results】With rice straw returning for two consecutive years, grain yield of each treatment was significantly increased, among which, T2treatment with lime and urea was the best. The yield of tobacco and rice under T2treatment increased by 25.67% and 11.49% respectively compared with T0, followed by T3treatment with Fuganlin and urea, and T1treatment with straw returning.In the first year, the yield of tobacco in T1treatment was not significantly different from that of the control, but in the second year, the yields of tobacco and rice were 3.77% and 5.90% higher than those of the control, respectively. After two consecutive years of straw returning, the nutrient and organic carbon contents of topsoil in each treatment were significantly increased, T2treatment ranked first, followed by T3treatment.The contents of alkali-hydrolyzed nitrogen, available phosphorus, available potassium, soil organic carbon and dissolved organic carbon in T2treatment increased by 7.09%, 5.97%, 3.01%, 11.32% and 5.47%, respectively compared with the control.The rhizosphere microbial community analysis showed that the relative abundance ofrelated to carbon source metabolism,related to the degradation of toxic substances andrelated to cellulolysis increased significantly in T2and T3treatments compared with the control, while the abundance ofandassociated with pathogenesis decreased significantly.Further correlation analysis found that the increase of beneficial bacteria in the rhizosphere were positively correlated with soil nutrients, while the abundance of pathogenic bacteria were negatively correlated with soil nutrients.【Conclusion】Rice-straw returning with lime or Fuganlin and appropriate amount of urea is conducive to increasing the nutrient contents in the topsoil, promoting the growth of beneficial microorganisms such as nitrogen-fixing bacteria, inhibiting pathogenic bacteria, and creating a good soil environment for tobacco-rice rotation, and thus improving crop yield.

    straw returning; tobacco-rice rotation; soil nutrients; rhizosphere microbes; community diversity

    10.16819/j.1001-7216.2023.220901

    2022-09-01;

    2022-10-09。

    中國(guó)煙草總公司福建省公司科技項(xiàng)目(20193500002400143)。

    猜你喜歡
    輪作根際稻草
    稻草熊
    根際微生物對(duì)植物與土壤交互調(diào)控的研究進(jìn)展
    農(nóng)民瓜菜“接茬輪作”增效益
    黃花蒿葉水提物對(duì)三七根際尖孢鐮刀菌生長(zhǎng)的抑制作用
    輪作休耕,讓土地歇一會(huì)兒
    耕地輪作效益調(diào)查報(bào)告
    促植物生長(zhǎng)根際細(xì)菌HG28-5對(duì)黃瓜苗期生長(zhǎng)及根際土壤微生態(tài)的影響
    兩根稻草
    稻草可喂牛 喂前需處理
    最后一根稻草
    日韩制服丝袜自拍偷拍| 一个人免费在线观看的高清视频| cao死你这个sao货| 国产单亲对白刺激| 黑人操中国人逼视频| 亚洲五月色婷婷综合| 精品亚洲成a人片在线观看| 国产在线免费精品| 久久这里只有精品19| 99久久精品国产亚洲精品| www.熟女人妻精品国产| 国产亚洲av高清不卡| 免费少妇av软件| 国产在线视频一区二区| 女同久久另类99精品国产91| 两个人免费观看高清视频| 国产精品av久久久久免费| 一区二区三区精品91| 午夜成年电影在线免费观看| 九色亚洲精品在线播放| 国产淫语在线视频| 亚洲精华国产精华精| 成人国语在线视频| 一本大道久久a久久精品| 久久久久久久大尺度免费视频| 久热这里只有精品99| 男女免费视频国产| 成人18禁高潮啪啪吃奶动态图| 国产精品国产高清国产av | 欧美激情 高清一区二区三区| 久久久久久久精品吃奶| av网站在线播放免费| 咕卡用的链子| 久久久国产欧美日韩av| 99国产综合亚洲精品| 黄色片一级片一级黄色片| 亚洲欧美日韩高清在线视频 | av福利片在线| 久久 成人 亚洲| 美女视频免费永久观看网站| 日本黄色日本黄色录像| 99精品久久久久人妻精品| 久久久久久久久免费视频了| 亚洲精品国产区一区二| 黄片小视频在线播放| 亚洲av日韩在线播放| 十分钟在线观看高清视频www| 国产成人免费无遮挡视频| 另类精品久久| 一本一本久久a久久精品综合妖精| 一边摸一边抽搐一进一小说 | 色婷婷久久久亚洲欧美| 人妻 亚洲 视频| 大型av网站在线播放| 亚洲综合色网址| 搡老熟女国产l中国老女人| 狠狠精品人妻久久久久久综合| 久久久久精品人妻al黑| bbb黄色大片| xxxhd国产人妻xxx| 乱人伦中国视频| 久久精品人人爽人人爽视色| 菩萨蛮人人尽说江南好唐韦庄| 亚洲午夜理论影院| 精品一区二区三卡| 91国产中文字幕| 久久久久视频综合| 777久久人妻少妇嫩草av网站| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 欧美国产精品一级二级三级| 精品国产亚洲在线| 捣出白浆h1v1| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 精品亚洲成国产av| 亚洲黑人精品在线| 成人18禁高潮啪啪吃奶动态图| 日韩一区二区三区影片| 十八禁网站免费在线| 国产av又大| 午夜免费鲁丝| av天堂在线播放| 国产精品自产拍在线观看55亚洲 | 一级毛片女人18水好多| 国产欧美日韩一区二区三区在线| 成人亚洲精品一区在线观看| 成年女人毛片免费观看观看9 | 成人特级黄色片久久久久久久 | 中文字幕色久视频| 水蜜桃什么品种好| 久久中文看片网| 久久九九热精品免费| 亚洲欧美日韩高清在线视频 | 亚洲免费av在线视频| 一本综合久久免费| 亚洲国产看品久久| 国产精品偷伦视频观看了| 黄片播放在线免费| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡| 国产又爽黄色视频| 一边摸一边抽搐一进一小说 | 老熟妇乱子伦视频在线观看| 久久毛片免费看一区二区三区| 麻豆av在线久日| 午夜成年电影在线免费观看| 国产精品久久久久久精品古装| 91精品三级在线观看| 亚洲色图 男人天堂 中文字幕| 国产99久久九九免费精品| 少妇粗大呻吟视频| 国产av精品麻豆| 亚洲免费av在线视频| 91九色精品人成在线观看| 日韩成人在线观看一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 久久久欧美国产精品| 99riav亚洲国产免费| 欧美久久黑人一区二区| 激情视频va一区二区三区| 欧美黑人欧美精品刺激| 最黄视频免费看| 亚洲av片天天在线观看| 悠悠久久av| 国产三级黄色录像| 日韩大片免费观看网站| 亚洲欧美激情在线| 亚洲中文av在线| 成人国语在线视频| 亚洲精品一卡2卡三卡4卡5卡| 深夜精品福利| 亚洲成国产人片在线观看| 九色亚洲精品在线播放| 日本黄色日本黄色录像| 91麻豆av在线| 少妇粗大呻吟视频| 99国产精品一区二区蜜桃av | 激情视频va一区二区三区| 大香蕉久久网| 国产精品免费大片| 在线亚洲精品国产二区图片欧美| 亚洲国产精品一区二区三区在线| 久久青草综合色| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 热re99久久国产66热| 亚洲第一av免费看| 国产黄频视频在线观看| 一本大道久久a久久精品| 久久精品国产亚洲av香蕉五月 | 亚洲 欧美一区二区三区| 美女高潮喷水抽搐中文字幕| 国产麻豆69| 国产精品国产高清国产av | 欧美成狂野欧美在线观看| 大片免费播放器 马上看| 天天影视国产精品| 久久香蕉激情| 激情在线观看视频在线高清 | 视频在线观看一区二区三区| 一级毛片女人18水好多| 欧美日韩亚洲高清精品| 老司机影院毛片| 国产伦人伦偷精品视频| av一本久久久久| 纯流量卡能插随身wifi吗| 99在线人妻在线中文字幕 | 欧美变态另类bdsm刘玥| 别揉我奶头~嗯~啊~动态视频| 国产一区二区 视频在线| 侵犯人妻中文字幕一二三四区| 欧美黄色片欧美黄色片| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 免费观看av网站的网址| 国精品久久久久久国模美| 日本欧美视频一区| 国产精品成人在线| 男女高潮啪啪啪动态图| 亚洲成人免费av在线播放| 俄罗斯特黄特色一大片| 亚洲精品国产色婷婷电影| 欧美日韩国产mv在线观看视频| 国产亚洲精品第一综合不卡| videosex国产| 久久久国产一区二区| 欧美亚洲日本最大视频资源| 露出奶头的视频| 欧美精品一区二区免费开放| 岛国在线观看网站| 精品人妻在线不人妻| 亚洲天堂av无毛| 国产精品一区二区在线不卡| 又紧又爽又黄一区二区| 久久人妻熟女aⅴ| 日本wwww免费看| 亚洲专区国产一区二区| 久久狼人影院| 1024视频免费在线观看| 一本久久精品| 80岁老熟妇乱子伦牲交| 老司机午夜十八禁免费视频| 亚洲情色 制服丝袜| 999精品在线视频| 亚洲欧美一区二区三区黑人| 丁香六月欧美| 黄片小视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自偷自拍图片 自拍| 99re6热这里在线精品视频| 国产熟女午夜一区二区三区| 大片免费播放器 马上看| 午夜福利欧美成人| 久久天躁狠狠躁夜夜2o2o| 男女高潮啪啪啪动态图| 久久久久久久大尺度免费视频| 久久久久精品人妻al黑| 可以免费在线观看a视频的电影网站| 老司机深夜福利视频在线观看| 国产真人三级小视频在线观看| 女人久久www免费人成看片| 欧美国产精品一级二级三级| 天堂俺去俺来也www色官网| 人妻一区二区av| 久9热在线精品视频| 怎么达到女性高潮| 午夜福利乱码中文字幕| 亚洲熟女精品中文字幕| 99久久99久久久精品蜜桃| kizo精华| 久久香蕉激情| 女人爽到高潮嗷嗷叫在线视频| 国产免费av片在线观看野外av| 日本av免费视频播放| 国产1区2区3区精品| 99国产极品粉嫩在线观看| 免费观看av网站的网址| 一级毛片女人18水好多| 亚洲精品在线观看二区| 交换朋友夫妻互换小说| 精品视频人人做人人爽| 久久中文字幕一级| 久久久国产一区二区| 国产精品久久久久久精品电影小说| 99re在线观看精品视频| 久久精品成人免费网站| 一本大道久久a久久精品| 99re在线观看精品视频| 国产精品 国内视频| 精品一区二区三区av网在线观看 | 久久99一区二区三区| 亚洲欧洲日产国产| 黄色 视频免费看| 免费高清在线观看日韩| 可以免费在线观看a视频的电影网站| 12—13女人毛片做爰片一| 人人妻人人澡人人爽人人夜夜| 老司机午夜十八禁免费视频| 一本综合久久免费| 一级黄色大片毛片| 日本黄色视频三级网站网址 | 欧美精品一区二区大全| 一区二区三区精品91| 日韩中文字幕视频在线看片| 久久久久久久久免费视频了| 女性被躁到高潮视频| 99久久人妻综合| 久久天堂一区二区三区四区| 免费女性裸体啪啪无遮挡网站| 男人操女人黄网站| 黄色丝袜av网址大全| av不卡在线播放| 91成人精品电影| 久久久久精品人妻al黑| 精品福利永久在线观看| 欧美日韩福利视频一区二区| 中文字幕另类日韩欧美亚洲嫩草| 18禁黄网站禁片午夜丰满| 亚洲欧洲精品一区二区精品久久久| av不卡在线播放| 国产精品影院久久| 亚洲男人天堂网一区| 国产精品 国内视频| 午夜福利免费观看在线| 久久久国产成人免费| 法律面前人人平等表现在哪些方面| 国产精品久久久久成人av| 新久久久久国产一级毛片| 一级,二级,三级黄色视频| 男女高潮啪啪啪动态图| 久久亚洲精品不卡| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频| 满18在线观看网站| 高清欧美精品videossex| 美女高潮喷水抽搐中文字幕| 可以免费在线观看a视频的电影网站| 伊人久久大香线蕉亚洲五| 肉色欧美久久久久久久蜜桃| 人人妻人人添人人爽欧美一区卜| 亚洲美女黄片视频| 午夜老司机福利片| 久久国产精品影院| 极品教师在线免费播放| 在线观看免费高清a一片| 黄色成人免费大全| 日韩视频一区二区在线观看| 搡老乐熟女国产| 国产日韩欧美视频二区| 电影成人av| 日日摸夜夜添夜夜添小说| 欧美亚洲日本最大视频资源| 午夜福利乱码中文字幕| 真人做人爱边吃奶动态| 侵犯人妻中文字幕一二三四区| 精品一区二区三区视频在线观看免费 | 免费少妇av软件| 91精品三级在线观看| 欧美在线一区亚洲| 久久这里只有精品19| 搡老乐熟女国产| 丰满少妇做爰视频| 老熟妇仑乱视频hdxx| 精品视频人人做人人爽| avwww免费| 国产欧美亚洲国产| 亚洲人成电影观看| 久久久久视频综合| 久久青草综合色| 美女国产高潮福利片在线看| 麻豆乱淫一区二区| 欧美黑人精品巨大| 日韩熟女老妇一区二区性免费视频| 香蕉国产在线看| 日韩欧美一区视频在线观看| 亚洲性夜色夜夜综合| av天堂久久9| 嫩草影视91久久| 久久久久精品人妻al黑| 日本av手机在线免费观看| 国产xxxxx性猛交| 国产av又大| 99九九在线精品视频| 18禁观看日本| 这个男人来自地球电影免费观看| 成年人免费黄色播放视频| 两人在一起打扑克的视频| 亚洲成人手机| 1024香蕉在线观看| 久久久久国内视频| 久久久久网色| 黄色视频不卡| 欧美日韩成人在线一区二区| 最新的欧美精品一区二区| 国产精品欧美亚洲77777| 99久久99久久久精品蜜桃| h视频一区二区三区| 亚洲精品自拍成人| 天天躁日日躁夜夜躁夜夜| 丁香六月天网| 18禁裸乳无遮挡动漫免费视频| 黄网站色视频无遮挡免费观看| 99riav亚洲国产免费| 王馨瑶露胸无遮挡在线观看| 一个人免费在线观看的高清视频| 午夜免费鲁丝| 国产一区二区三区综合在线观看| 久久久久精品国产欧美久久久| 少妇被粗大的猛进出69影院| 成年女人毛片免费观看观看9 | 九色亚洲精品在线播放| 日韩视频一区二区在线观看| 国产激情久久老熟女| 欧美黑人精品巨大| 欧美黄色淫秽网站| 搡老岳熟女国产| 国产极品粉嫩免费观看在线| 婷婷成人精品国产| 露出奶头的视频| 国产av又大| 日韩欧美一区二区三区在线观看 | 搡老乐熟女国产| 久热爱精品视频在线9| 在线观看免费视频网站a站| 国产精品av久久久久免费| 欧美人与性动交α欧美精品济南到| 搡老熟女国产l中国老女人| 少妇粗大呻吟视频| 成年人午夜在线观看视频| 激情视频va一区二区三区| 黄色视频,在线免费观看| 亚洲第一欧美日韩一区二区三区 | 黄网站色视频无遮挡免费观看| 免费少妇av软件| 91国产中文字幕| 日日夜夜操网爽| 日韩欧美国产一区二区入口| videosex国产| 亚洲成a人片在线一区二区| 99国产精品一区二区三区| 国产精品久久久人人做人人爽| 99热网站在线观看| 曰老女人黄片| 亚洲久久久国产精品| 三上悠亚av全集在线观看| 性高湖久久久久久久久免费观看| 久久久久国内视频| 国产激情久久老熟女| 一本久久精品| 日韩制服丝袜自拍偷拍| 国产黄色免费在线视频| www.熟女人妻精品国产| 欧美日韩黄片免| 久久亚洲真实| 91老司机精品| 伊人久久大香线蕉亚洲五| 久久中文字幕一级| 亚洲欧美日韩高清在线视频 | 亚洲av日韩精品久久久久久密| 夫妻午夜视频| 热re99久久精品国产66热6| 精品久久久久久久毛片微露脸| 亚洲色图综合在线观看| 亚洲久久久国产精品| 两人在一起打扑克的视频| 亚洲第一av免费看| 色婷婷av一区二区三区视频| 国产淫语在线视频| 国产97色在线日韩免费| 国产主播在线观看一区二区| 国产精品亚洲av一区麻豆| 精品人妻1区二区| 久久天躁狠狠躁夜夜2o2o| 久久青草综合色| 亚洲国产av新网站| 女性被躁到高潮视频| 91成人精品电影| 成人手机av| 美女福利国产在线| 国产又色又爽无遮挡免费看| 欧美成狂野欧美在线观看| 一级毛片女人18水好多| 精品亚洲成国产av| 99riav亚洲国产免费| av不卡在线播放| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 欧美 亚洲 国产 日韩一| 午夜福利视频在线观看免费| 国产av精品麻豆| 国产精品熟女久久久久浪| 免费观看人在逋| 黄色成人免费大全| 国产精品久久久久久人妻精品电影 | 亚洲欧美日韩另类电影网站| av网站免费在线观看视频| 一区二区三区乱码不卡18| 啪啪无遮挡十八禁网站| 亚洲一区中文字幕在线| 欧美乱妇无乱码| 日日夜夜操网爽| 欧美午夜高清在线| 亚洲精品国产色婷婷电影| 亚洲 国产 在线| 日韩大片免费观看网站| 午夜福利视频在线观看免费| 国产高清国产精品国产三级| 色在线成人网| 亚洲精品国产色婷婷电影| 少妇粗大呻吟视频| 少妇的丰满在线观看| 乱人伦中国视频| 大片电影免费在线观看免费| 欧美在线黄色| av一本久久久久| 99久久人妻综合| 久久久久国内视频| 亚洲专区字幕在线| 免费人妻精品一区二区三区视频| 久久精品亚洲熟妇少妇任你| 久久性视频一级片| 亚洲五月色婷婷综合| 一进一出抽搐动态| 999久久久精品免费观看国产| 精品久久久久久电影网| 日韩一区二区三区影片| 免费看十八禁软件| 丰满人妻熟妇乱又伦精品不卡| 99re在线观看精品视频| 日韩一卡2卡3卡4卡2021年| 亚洲第一av免费看| 青青草视频在线视频观看| 两个人免费观看高清视频| 国产在线精品亚洲第一网站| 国产91精品成人一区二区三区 | 成人国产av品久久久| 精品一区二区三卡| 久久国产精品人妻蜜桃| 国产在线观看jvid| 午夜免费鲁丝| 蜜桃在线观看..| 真人做人爱边吃奶动态| 欧美激情极品国产一区二区三区| 成年人午夜在线观看视频| 满18在线观看网站| 少妇粗大呻吟视频| 视频区图区小说| 成年动漫av网址| 国产亚洲午夜精品一区二区久久| 美女高潮到喷水免费观看| 真人做人爱边吃奶动态| 日韩人妻精品一区2区三区| 亚洲精品av麻豆狂野| 国产福利在线免费观看视频| 不卡av一区二区三区| 香蕉国产在线看| 乱人伦中国视频| 桃花免费在线播放| 久久狼人影院| 高清在线国产一区| 久久精品人人爽人人爽视色| 日韩免费高清中文字幕av| 亚洲av国产av综合av卡| 国产一区二区三区视频了| 中文字幕色久视频| 一区福利在线观看| 亚洲中文日韩欧美视频| 免费高清在线观看日韩| 午夜两性在线视频| 国产精品秋霞免费鲁丝片| 午夜福利视频在线观看免费| 午夜激情久久久久久久| 成年人免费黄色播放视频| 在线观看免费高清a一片| 男女免费视频国产| 亚洲五月色婷婷综合| 久久精品亚洲av国产电影网| 亚洲国产毛片av蜜桃av| 免费久久久久久久精品成人欧美视频| 少妇被粗大的猛进出69影院| 久久精品国产综合久久久| 91av网站免费观看| 黄网站色视频无遮挡免费观看| 如日韩欧美国产精品一区二区三区| 国产一区有黄有色的免费视频| 久久久久视频综合| 啦啦啦 在线观看视频| 下体分泌物呈黄色| 汤姆久久久久久久影院中文字幕| 久久久国产欧美日韩av| 欧美成人免费av一区二区三区 | 精品视频人人做人人爽| √禁漫天堂资源中文www| 国产精品一区二区在线观看99| 午夜福利影视在线免费观看| 极品教师在线免费播放| 成年版毛片免费区| 日韩视频一区二区在线观看| 日本黄色视频三级网站网址 | www.熟女人妻精品国产| 亚洲成人国产一区在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲av欧美aⅴ国产| 韩国精品一区二区三区| 老司机午夜十八禁免费视频| 久久av网站| 老司机深夜福利视频在线观看| 色精品久久人妻99蜜桃| 99国产极品粉嫩在线观看| 大码成人一级视频| av有码第一页| 亚洲专区字幕在线| 99香蕉大伊视频| 国产深夜福利视频在线观看| 最近最新中文字幕大全电影3 | 国产在线视频一区二区| 在线 av 中文字幕| 午夜老司机福利片| 日韩欧美三级三区| 99国产精品99久久久久| 欧美日韩成人在线一区二区| 纯流量卡能插随身wifi吗| 高清在线国产一区| 欧美黄色片欧美黄色片| 国产老妇伦熟女老妇高清| 最新的欧美精品一区二区| 99国产极品粉嫩在线观看| 天堂8中文在线网| 国产三级黄色录像| 日韩一区二区三区影片| 操出白浆在线播放| 免费黄频网站在线观看国产| √禁漫天堂资源中文www| 精品久久久精品久久久| 精品乱码久久久久久99久播| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 超碰97精品在线观看| 欧美中文综合在线视频| 18禁国产床啪视频网站| 女性被躁到高潮视频| 嫁个100分男人电影在线观看| 欧美在线黄色| 国产亚洲午夜精品一区二区久久| 国产精品av久久久久免费| 777久久人妻少妇嫩草av网站| 久久久久久亚洲精品国产蜜桃av| 国产不卡一卡二| 国产成人精品久久二区二区免费| 久久九九热精品免费|