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

    長(zhǎng)期化肥減量配施秸稈對(duì)水稻土壤微生物多樣性和功能的影響

    2024-11-07 00:00:00曾憲楠王曉慧孫羽王麒宋秋來(lái)梁全喜
    江蘇農(nóng)業(yè)科學(xué) 2024年18期

    摘要:為探究長(zhǎng)期化肥減量配施秸稈對(duì)水稻土壤微生物結(jié)構(gòu)和功能的影響,以東北黑土為試驗(yàn)對(duì)象,基于黑龍江省農(nóng)業(yè)科學(xué)院長(zhǎng)期水稻秸稈還田和化肥減施平臺(tái),利用高通量測(cè)序技術(shù)對(duì)化肥不同施加量[N1(173 kg/hm2)、N2(133 kg/hm2)、N3(93 kg/hm2)、N4(53 kg/hm2)、N5(0 kg/hm2)]配施水稻秸稈(所有處理秸稈還田量均為7 500 kg/hm2)土壤細(xì)菌和真菌群落結(jié)構(gòu)和功能進(jìn)行研究,分析土壤微生物群落結(jié)構(gòu)和功能與土壤理化性質(zhì)之間的關(guān)系。結(jié)果表明,不同化肥施加量配施水稻秸稈顯著改變了土壤全氮、全磷和速效鉀的含量。秸稈還田下化肥施加量處理對(duì)土壤細(xì)菌香農(nóng)指數(shù)和Chao1指數(shù)影響顯著,但對(duì)土壤真菌α多樣性指數(shù)并未產(chǎn)生顯著影響。主坐標(biāo)結(jié)果表明,秸稈還田下化肥施加量處理可對(duì)土壤細(xì)菌和真菌的群落結(jié)構(gòu)產(chǎn)生顯著影響。酸桿菌門(mén)、放線菌門(mén)和擬桿菌門(mén)是相對(duì)豐度最高的3個(gè)細(xì)菌門(mén)。擔(dān)子菌門(mén)、子囊菌門(mén)和壺菌門(mén)是相對(duì)豐度最高的3個(gè)真菌門(mén)?;贔APROTAX分析土壤細(xì)菌的功能,所有處理中土壤細(xì)菌的主要功能是chemoheterotrophy(20.71%)>aerobic_chemoheterotrophy(11.37%)>nitrification(7.54%)?;贔unGuild分析土壤真菌功能,所有處理中土壤真菌的主要功能是Saprotroph>Pathotroph>Symbiotroph。結(jié)構(gòu)方程表明,在秸稈還田基礎(chǔ)上氮肥的減少會(huì)使土壤中的N和P累積,但會(huì)使土壤有機(jī)質(zhì)含量增多,并通過(guò)增加土壤有機(jī)質(zhì)的含量影響土壤細(xì)菌多樣性?;侍砑恿繒?huì)改變土壤pH值而影響土壤真菌多樣性。土壤細(xì)菌多樣性與土壤全氮、硝態(tài)氮、總鉀和速效鉀含量呈正相關(guān),而土壤真菌多樣性與速效磷、有機(jī)質(zhì)含572t9bJMg5f41x9u2hkQdsdMOWXWD8EHmC6hXfzKuN8=量呈正相關(guān)。本研究結(jié)果表明減少氮肥施加是提高土壤微生物多樣性,改善土壤質(zhì)量的有效手段。

    關(guān)鍵詞:秸稈;水稻土;高通量測(cè)序;結(jié)構(gòu)方程模型;網(wǎng)絡(luò)分析

    中圖分類(lèi)號(hào):S511.06 文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1002-1302(2024)18-0251-10

    收稿日期:2023-10-31

    基金項(xiàng)目:黑龍江省省屬科研院所科研業(yè)務(wù)費(fèi)(編號(hào):CZKYF2021-2-C027)。

    作者簡(jiǎn)介:曾憲楠(1985—),女,黑龍江綏化人,碩士,副研究員,主要從事保護(hù)性耕作研究。E-mail:zengxiannanzxn@163.com。

    水稻是主要的糧食作物,水稻全球的種植面積超過(guò)1.62億hm2,是全球半數(shù)以上人口的主食1-2,也是我國(guó)主要的糧食作物[3。隨著經(jīng)濟(jì)的快速發(fā)展和為了滿(mǎn)足人口對(duì)糧食的最基本需求,需要提高水稻產(chǎn)量。追求產(chǎn)量增加導(dǎo)致人們過(guò)度施用化肥,我國(guó)已成為世界上最大的化肥消費(fèi)國(guó)4。長(zhǎng)期大量、單純地使用化肥導(dǎo)致溫室氣體排放增加5-6,引起氮磷水污染7,土壤退化8,降低化肥使用效率,影響作物產(chǎn)量9,農(nóng)業(yè)的可持續(xù)發(fā)展受到嚴(yán)重威脅。保障國(guó)家糧食安全,推進(jìn)農(nóng)業(yè)綠色可持續(xù)發(fā)展,保持土壤肥力的前提下,化肥的合理施用已成為農(nóng)業(yè)生產(chǎn)的研究熱點(diǎn)。

    我國(guó)作為農(nóng)業(yè)大國(guó),具有豐富的農(nóng)作物秸稈資源,每年農(nóng)作物秸稈產(chǎn)出量將近7億t[10-11。農(nóng)作物秸稈富含農(nóng)作物生長(zhǎng)所需的主要營(yíng)養(yǎng)元素,秸稈還田可改良土壤結(jié)構(gòu),提高土壤肥力12-13。已有學(xué)者對(duì)化肥減量配施秸稈進(jìn)行研究。吳立鵬等研究發(fā)現(xiàn),與常規(guī)施化肥比較,秸稈還田與減量氮肥配施使稻田土壤有機(jī)碳、可溶性有機(jī)碳含量增加[14。對(duì)于小麥產(chǎn)量,與常規(guī)施化肥比較,秸稈還田與減量氮肥配施對(duì)小麥產(chǎn)量無(wú)負(fù)面影響15。秸稈還田與減量氮肥配施是減少化肥用量的一種有效措施。

    土壤中微生物種類(lèi)豐富且數(shù)量眾多,在土壤生態(tài)功能,尤其是元素循環(huán)過(guò)程中扮演著重要角色,對(duì)植物的健康起到重要作用[16-17。土壤微生物對(duì)環(huán)境變化較為敏感,微生物群落的改變從一定程度上反映了土壤生態(tài)功能變化。因此,研究長(zhǎng)期化肥減施配秸稈還田對(duì)土壤微生物群落的影響,對(duì)評(píng)價(jià)秸稈還田和配施化肥對(duì)土壤生態(tài)的影響以及尋找環(huán)境友好型的農(nóng)業(yè)措施具有重要的指導(dǎo)意義。已有的研究發(fā)現(xiàn)使用化肥減施和秸稈還田會(huì)顯著改變土壤微生物學(xué)特性。Huang等的研究表明,秸稈與化肥減量配施可以提高紫色土土壤氨氧化細(xì)菌和氨氧化古菌的amoA基因多樣性18。楊濱娟等研究秸稈還田與不同比例化肥配施對(duì)稻田土壤根際微生物的影響,發(fā)現(xiàn)秸稈還田處理能夠增加固氮菌的數(shù)量[19??梢?jiàn)不同農(nóng)作物的秸稈、不同的土壤類(lèi)型對(duì)土壤微生物群落結(jié)構(gòu)和豐度的影響不一致。

    黑土具有土質(zhì)肥沃、保肥、適宜農(nóng)作物生長(zhǎng)等良好的農(nóng)業(yè)生產(chǎn)條件,黑龍江省位于我國(guó)黑土帶上,已成為我國(guó)重要的商品糧生產(chǎn)基地,也是我國(guó)糧食安全的“壓艙石”[20。對(duì)水稻秸稈還田的研究較多,以往的研究主要集中探討了秸稈還田對(duì)作物產(chǎn)量21-22或土壤理化性質(zhì)的影響23-24,但是關(guān)于水稻秸稈還田與化肥減量配施后,東北黑土土壤微生物群落結(jié)構(gòu)和多樣性變化的研究較少。為此,以東北黑土為對(duì)象,基于黑龍江省農(nóng)業(yè)科學(xué)院長(zhǎng)期水稻秸稈還田和化肥減施平臺(tái),應(yīng)用高通量測(cè)序技術(shù),研究水稻秸稈與化肥減量配施對(duì)土壤細(xì)菌和真菌群落結(jié)構(gòu)的影響及與土壤理化性質(zhì)的關(guān)系,以期為寒地稻田合理施肥提供科學(xué)依據(jù)。

    1 材料與方法

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

    試驗(yàn)于2017年在黑龍江省農(nóng)業(yè)科學(xué)院現(xiàn)代農(nóng)業(yè)示范區(qū)(45°49′N(xiāo),126°48′E,海拔117 m)進(jìn)行。該地區(qū)屬中溫帶大陸性季風(fēng)氣候,冬夏季節(jié)短,年平均氣溫為4.5 ℃,年平均氣溫≥10 ℃有效積溫2 700 ℃以上,作物年生長(zhǎng)周期約為155 d,年平均降水量為569.1 mm,降水主要集中在6—9月。供試土壤類(lèi)型為黑土。0~20 cm土層土壤養(yǎng)分含量如下:土壤有機(jī)質(zhì)含量為26.50 g/kg,全氮含量為2.01 g/kg,有效氮含量為79.56 mg/kg,有效磷含量為55.84 mg/kg,速效鉀含量為168.42 mg/kg。

    1.2 試驗(yàn)材料

    供試水稻品種:龍稻21(培育單位為黑龍江省農(nóng)業(yè)科學(xué)院耕作栽培研究所);供試秸稈為單季水稻殘?bào)w,生長(zhǎng)周期完成后,冬季不留茬。秋季收獲水稻,將試驗(yàn)區(qū)所有秸稈搬出區(qū)外,人工粉碎成約5 cm,秸稈按照不同處理用量,均勻拋撒,翻埋殘茬秸稈,越冬;翌年春季泡田后使用攪漿平地機(jī)平整土地,插秧。田間管理措施同當(dāng)?shù)爻R?guī)管理一致。

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

    采用單因素試驗(yàn)設(shè)計(jì),共有5個(gè)氮肥施入量,分別為N1(173 kg/hm2)、N2(133 kg/hm2)、N3(93 kg/hm2)、N4(53 kg/hm2)、N5(0 kg/hm2)(表1)。所有處理秸稈還田量均為7 500 kg/hm2。所有處理的磷、鉀施入量分別為純磷46 kg/hm2,純鉀75 kg/hm2(依據(jù)當(dāng)?shù)厥┓史椒坝昧浚?/p>

    1.4 樣本采集

    于2023年水稻收獲期進(jìn)行土壤樣品采集。用土鉆在每個(gè)小區(qū)按照S型隨機(jī)采集5點(diǎn),采集0~20 cm土層土壤,5點(diǎn)混合為1個(gè)復(fù)合樣品。去除土壤中的石塊和植物殘茬,用2 mm的篩子篩土。一部分土壤樣品經(jīng)風(fēng)干后用于測(cè)定土壤的理化性質(zhì)。另一部分樣品在-80 ℃冰箱中冷凍保存,用于微生物DNA的提取和群落分析。

    1.5 測(cè)定方法

    1.5.1 土壤化學(xué)性質(zhì)測(cè)定

    全氮含量采用凱氏定氮法測(cè)定,堿解氮含量采用堿解擴(kuò)散法測(cè)定,銨態(tài)氮含量采用KCl浸提蒸餾-凱氏定氮法測(cè)定,硝態(tài)氮含量采用酚二磺酸法測(cè)定,土壤總有機(jī)碳含量采用重鉻酸鉀法測(cè)定,總磷含量采用鉬銻抗分光光度法測(cè)定,速效磷含量采用Olsen法測(cè)定,速效鉀含量采用火焰光度法測(cè)定,pH值采用電位法測(cè)定[25

    1.5.2 微生物測(cè)序

    使用Mobio12888 power soil extraction kit提取土壤總DNA。本研究通過(guò)Illumina MiSeq測(cè)序平臺(tái)檢測(cè)土壤細(xì)菌16S rRNA基因片段序列,引物是16S rRNA基因V4區(qū)的通用引物515F(正向引物:5′-GTGCCAGCMGCCGCGGTAA-3′)和806R(反向引物,5′-GACTACHVGGGTWTCTAAT-3′)。真菌使用ITS rDNA上的ITS1-ITS2可變區(qū)ITS1(5′-CTTGGTCATTTAGAGGAAGTAA-3′) 和ITS2 (5′-GCTGCGTTCATCGATGC-3′)。每個(gè)樣品含有特有的6 bp條形碼序列被添加到引物中,用于區(qū)分多個(gè)樣品。PCR反應(yīng)體系共25 μL:含有 2.5 μL TransStart Buffer、2 μL dNTP、1 μL 正向/反向引物 (10 ng/μL)和1 μL DNA(30 ng/μL)模版,17.5 μL去離子水。 PCR條件如下:94 ℃預(yù)變性5 min;94 ℃變性30 s,55 ℃ 30 s,72 ℃45 s,30個(gè)循環(huán);72 ℃最終延伸10 min。PCR產(chǎn)物經(jīng)2%瓊脂糖電泳檢測(cè),并使用AxyPrep DNA純化試劑盒進(jìn)行純化。每個(gè)樣品分別進(jìn)行3個(gè)PCR重復(fù),然后等量合并3個(gè)PCR樣品,并在Illumina Miseq v3平臺(tái) (2×300 bp) 上進(jìn)行PE300雙端測(cè)序。

    1.5.3 高通量測(cè)序數(shù)據(jù)分析

    使用北京百邁克生物信息云平臺(tái)上的QIIME1(版本1.81)軟件對(duì)原始序列進(jìn)行分析。

    原始序列按照以下標(biāo)準(zhǔn)進(jìn)行質(zhì)量過(guò)濾:去除平均得分<20或長(zhǎng)度<200 bp的原始序列。 使用 PEAR軟件(版本0.9.8)對(duì)雙端測(cè)序序列進(jìn)行拼合。使用Usearch軟件(版本7.1)去除嵌合體。使用UPARSE軟件以 97% 的相似度生成操作分類(lèi)單位 (OTU)并對(duì)OTU的代表性序列進(jìn)行分類(lèi)學(xué)注釋?zhuān)眯哦乳撝禐?.7,細(xì)菌比對(duì)數(shù)據(jù)庫(kù)為SILVA(v.3.2),真菌比對(duì)數(shù)據(jù)庫(kù)為UNITE(v.8.2)。所有序列按照單個(gè)樣本的最低序列數(shù)對(duì)序列進(jìn)行歸一化后進(jìn)行后續(xù)的統(tǒng)計(jì)分析。

    1.6 數(shù)據(jù)分析

    基于OTU表計(jì)算所有樣品的α多樣性指數(shù)和β多樣性指數(shù)(主坐標(biāo)分析),使用R軟件(v.3.2.5)的Vegan軟件包計(jì)算。采用FAPROTAX功能預(yù)測(cè)分析環(huán)境樣品中細(xì)菌分為80多個(gè)功能分組(如硝酸鹽呼吸、產(chǎn)甲烷、發(fā)酵、植物病原等)。采用FunGuild功能預(yù)測(cè)將環(huán)境樣品中真菌分為病理營(yíng)養(yǎng)型(pathotroph)、共生營(yíng)養(yǎng)型(symbiotroph)、腐生營(yíng)養(yǎng)型(saprotroph) 3個(gè)營(yíng)養(yǎng)型。利用Spearman相關(guān)方法分析了土壤微生物與養(yǎng)分含量之間的關(guān)系,并構(gòu)建相關(guān)矩陣。利用通徑系數(shù)(也稱(chēng)為結(jié)構(gòu)方程模型)分析了N添加、微生物多樣性、土壤理化性質(zhì)之間的關(guān)系。使用AMOS 22.0(IBM Corporation,美國(guó))軟件對(duì)SEM進(jìn)行分析。使用SPSS 22.0進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析。分別使用Kolmogorov-Smirnov和Levene檢驗(yàn)分析了數(shù)據(jù)分布的正態(tài)性和異方差特征。使用Origin 8.0軟件繪制圖表。

    2 結(jié)果與分析

    2.1 長(zhǎng)期化肥減量配施秸稈對(duì)水稻土壤化學(xué)性質(zhì)的影響

    在秸稈還田下,不同化肥施加量?jī)H對(duì)土壤全氮、全磷和速效鉀含量有顯著影響(圖1)。其中,土壤全氮含量隨著化肥施加量減少呈先下降后上升再下降的趨勢(shì)。土壤全磷含量隨著化肥施加量的減少呈先上升后下降趨勢(shì),在N2處理達(dá)到最高值。土壤速效鉀含量呈現(xiàn)的規(guī)律為N4>N2>N1>N3>N5,N4與其他4個(gè)處理差異顯著。

    2.2 長(zhǎng)期化肥減量配施秸稈對(duì)水稻土壤微生物α多樣性的影響

    在秸稈還田下不同化肥施入量對(duì)土壤細(xì)菌香農(nóng)指數(shù)和Chao1指數(shù)的影響顯著(表2)。Chao1指數(shù)呈現(xiàn)的規(guī)律為N4>N5>N1>N2>N3,香農(nóng)指數(shù)呈現(xiàn)的規(guī)律是N5>N4>N1>N2>N3。但是土壤真菌香農(nóng)指數(shù)和Chao1指數(shù)在不同化肥施加量處理下差異并不顯著(表2)。

    2.3 長(zhǎng)期化肥減量配施秸稈對(duì)水稻土壤微生物β多樣性的影響

    基于Bray-Curtis距離的非度量多維排列NMDS,結(jié)果表明不同處理下土壤細(xì)菌(R2=0.19,圖2-A)和真菌群落組成(R2=0.31,圖2-B)均存在顯著差異。

    2.4 長(zhǎng)期化肥減量配施秸稈對(duì)水稻土細(xì)菌和真菌群落結(jié)構(gòu)組成的影響

    細(xì)菌共有OTU數(shù)量為1 969個(gè),N1處理特有OTU數(shù)目為123個(gè),而N5特有種OTU數(shù)目為101個(gè)(圖3-A)。真菌共有OTU數(shù)量為360個(gè),N1處理特有OTU數(shù)目為256個(gè),而N5特有種OTU數(shù)目為189個(gè)(圖3-B)。通過(guò)GraPhlAn可視化結(jié)果所示,細(xì)菌門(mén)水平上節(jié)點(diǎn)大小表現(xiàn)前3位為酸桿菌門(mén)(Acidobacteria)>放線菌門(mén)(Actinobacteria)>擬桿菌門(mén)(Bacteroidetes)(圖3-E)。真菌門(mén)水平上節(jié)點(diǎn)大小表現(xiàn)為擔(dān)子菌門(mén)(Basidiomycota)>子囊菌門(mén)(Ascomycota)>壺菌門(mén)(Chytridiomycota)(圖3-F)。細(xì)菌相對(duì)豐度屬水平前3位分別是uncultured_bacterium_f_Nitrosomonadaceae、uncultured_bacterium_f_Anaerolineaceae、Sphingomonas(圖3-C)。其中,N1處理的uncultured_bacterium_f_Nitrosomonadaceae和uncultured_bacterium_f_Nitrosomonadaceae高于其他處理組(圖3-C)。而Sphingomonas在N5處理中最高。真菌相對(duì)豐度屬水平前3位分別是Cladosporium、Mrakia、Schizothecium(圖3-D)。其中,在N3處理的Cladosporium達(dá)到最低相對(duì)豐度。在N4處理Mrakia為最低值,而Schizothecium為最高值(圖3-D)。

    2.5 相關(guān)性網(wǎng)絡(luò)與關(guān)鍵種

    微生物隨機(jī)分子生態(tài)網(wǎng)絡(luò)(MENA)結(jié)果顯示,土壤細(xì)菌中共有45個(gè)nodes,100個(gè)edges,modularity是0.460,graph_diameter是16.62,graph_density是0.101,clustering_coefficient是0.364,正相關(guān)68%,負(fù)相關(guān)32%(圖2-C)。其中,與周?chē)锓N聯(lián)系緊密并且主要呈正相關(guān)的主要是Proteobacteria門(mén),主要與周?chē)锓N呈負(fù)相關(guān)的是Bacteroidetes(圖2-C)。土壤真菌中共有50個(gè)nodes,100個(gè)edges,modularity是0.473,graph_diameter是19.59,graph_density是0.082,clustering_coefficient是0.309,正相關(guān)85%,負(fù)相關(guān)15%(圖2-D)。其中,與周?chē)锓N聯(lián)系緊密并且主要呈正相關(guān)的主要是Ascomycota,主要與周?chē)锓N呈負(fù)相關(guān)的是Basidiomycota(圖2-D)。Zi-pi圖結(jié)果表示土壤細(xì)菌中有5個(gè)關(guān)鍵屬(Connectors),分別為uncultured_Acidobacteria_bacterium、uncultured_bacterium_c_OPB35_soil_group、Sphingomonas、uncultured_bacterium_o_Gaiellales、Pseudarthrobacter。其中,Pseudarthrobacter、Sphingomonas在N5處理的相對(duì)豐度最高,分別為0.85%和2.47%。土壤真菌中有5個(gè)關(guān)鍵屬(Connectors),分別為Pyrenochaetopsis、Vishniacozyma、Stenella、Ramichloridium、Hortaea。其中,Pyrenochaetopsis和Hortaea在N4處理中相對(duì)豐度最高。

    2.6 微生物之功能預(yù)測(cè)

    由圖4-A可知,基于FAPROTAX分析土壤細(xì)菌的功能,所有土壤中相對(duì)豐度前3位的是chemoheterotrophy、aerobic_chemoheterotrophy、nitrification。chemoheterotrophy在N5處理達(dá)到最高值20.71%。nitrification、aerobic_ammonia_oxidation等功能相對(duì)豐度在N4處理達(dá)到最高值,分別為11.37%、7.54%?;贔unGuild分析土壤真菌功能,總體上3種營(yíng)養(yǎng)類(lèi)型豐度從高到低分別是嗜酸>病理營(yíng)養(yǎng)>共生營(yíng)養(yǎng)(圖4-B)。進(jìn)一步對(duì)真菌功能進(jìn)行分析,結(jié)果表明相對(duì)豐度前3位的是腐生真菌、蘭科腐生、植物病原體,且隨著化肥施加量的降低,腐生功能相對(duì)豐度先降低再升高再降低,而病原體表現(xiàn)出先升高再降低再升高的趨勢(shì),關(guān)鍵拐點(diǎn)均為N3和N4處理(圖4-C)。

    2.7 相關(guān)性分析

    結(jié)構(gòu)方程結(jié)果表明,在秸稈還田基礎(chǔ)上氮肥的減少會(huì)使土壤中的N和P累積、土壤細(xì)菌和真菌的多樣性、pH值減少,但土壤有機(jī)質(zhì)含量增多,并且通過(guò)增加土壤有機(jī)質(zhì)的含量影響土壤細(xì)菌多樣性(圖5-A)。 化肥添加量會(huì)改變土壤pH值而影響土壤真菌多樣性(圖5-A)。土壤細(xì)菌多樣性與土壤全氮、硝態(tài)氮、總鉀和速效鉀含量呈正相關(guān),而土壤真菌多樣性與速效磷、有機(jī)質(zhì)含量呈正相關(guān)(圖5-B)。

    3 討論與結(jié)論

    3.1 少量氮肥施加結(jié)合秸稈還田提高土壤微生物多樣性

    本研究發(fā)現(xiàn),化肥施加會(huì)顯著影響細(xì)菌α多樣性但是對(duì)真菌alpha多樣性影響不顯著(表2),并且少量添加氮肥(53 kg/hm2)和不添加氮肥(0 kg)下土壤中細(xì)菌明顯高于其他化肥添加量組(表2)。這是由于一方面適當(dāng)?shù)牡适┘邮谴龠M(jìn)微生物活動(dòng)的關(guān)鍵。少量的氮肥不僅可以為土壤中的微生物提供額外的氮源,促進(jìn)微生物的代謝活動(dòng)和生長(zhǎng),還可以改變土壤中的生態(tài)環(huán)境,刺激一些特定類(lèi)型微生物的生長(zhǎng)和多樣性,豐富土壤微生物群落,這與之前的研究結(jié)果[26-27一致。有研究表明單獨(dú)添加化肥可顯著改變土壤真菌多樣性,少量添加會(huì)增加土壤真菌多樣性,而大量添加會(huì)顯著抑制土壤真菌的多樣性28。本研究由于秸稈還田和化肥同時(shí)添加,秸稈的施加導(dǎo)致土壤中有機(jī)碳含量和種類(lèi)都發(fā)生了改變,反而抵消了化肥添加所導(dǎo)致的土壤養(yǎng)分之間的化學(xué)計(jì)量的改變,因此土壤真菌并未發(fā)生改變。馬玉穎等的研究結(jié)果表明,有機(jī)肥和秸稈還田會(huì)顯著提高砂姜黑土真菌多樣性[29。這表明土壤真菌對(duì)秸稈還田和化肥添加的響應(yīng)是不同的,這也與秸稈的質(zhì)量、土壤類(lèi)型等緊密相關(guān)29。秸稈還田能夠增加土壤的有機(jī)質(zhì)含量,豐富土壤中的碳源,為微生物提供多樣的有機(jī)質(zhì)和營(yíng)養(yǎng)物質(zhì)。秸稈還田還能夠改善土壤結(jié)構(gòu),增加團(tuán)聚體數(shù)量,增加土壤的通透性和空氣含量,為微生物提供適宜生長(zhǎng)的空間??傊?,在秸稈還田的基礎(chǔ)上只需要少量或者不施加氮肥即可維持土壤中細(xì)菌和真菌的多樣和豐富度。

    3.2 少量氮肥施加結(jié)合秸稈還田對(duì)土壤微生物組成的影響

    本研究發(fā)現(xiàn),優(yōu)勢(shì)細(xì)菌在未添加氮肥的處理中含量最高,Sphingomonas在未添加氮肥(0 kg)的處理組相對(duì)豐度最高(圖3-C)。據(jù)報(bào)道Sphingomonas能夠利用多種有機(jī)污染物作為碳源和能源,如多環(huán)芳烴、農(nóng)藥、染料等,從而減少這些物質(zhì)對(duì)植物和環(huán)境的危害、促進(jìn)植物生長(zhǎng),這有助于糧食安全生產(chǎn)和糧食高產(chǎn)。此外,uncultured_bacterium_f_Nitrosomonadaceae隨著氮肥減施呈先下降后上升的趨勢(shì),這表明本身土壤的硝化細(xì)菌也能夠滿(mǎn)足農(nóng)田土壤的氨氧化過(guò)程累積氮的需求。特別是,Sphingomonas既是優(yōu)勢(shì)屬也是關(guān)鍵屬。Pseudarthrobacter、Sphingomonas等關(guān)鍵屬與周?chē)锓N聯(lián)系密切,并且他們也在無(wú)氮肥添加處理中相對(duì)豐度最高。

    然而,優(yōu)勢(shì)真菌在少量氮肥的處理中表現(xiàn)優(yōu)異。Cladosporium是一種致病菌,而其相對(duì)豐度在N3和N4處理中較低。而Schizothecium相對(duì)豐度在N4處理中最高,Schizothecium能夠提高植物活力、植物鮮重、根鮮重和根干重,這使得植物能夠向土壤中分泌更多促進(jìn)微生物多樣性的物質(zhì)。此外,Hortaea、Pyrenochaetopsis屬等真菌與周?chē)锓N聯(lián)系緊密,并且也在低氮肥添加處理中相對(duì)豐度最高。據(jù)報(bào)道,Hortaea能有效防治植物和真菌性土傳病害,同時(shí)可使植物葉部的細(xì)菌和真菌病害明顯減少。Pyrenochaetopsis能夠直接參與調(diào)節(jié)土壤碳循環(huán)中的土壤呼吸作用和土壤酶活性,并且對(duì)土壤難溶性碳(有機(jī)質(zhì)等)具有較強(qiáng)的溶解能力[30??傊诮斩掃€田基礎(chǔ)上,減少氮肥施加量能夠調(diào)動(dòng)土壤中優(yōu)勢(shì)和關(guān)鍵物種來(lái)促進(jìn)碳氮循環(huán)、生物降解、病害防治、污染修復(fù)等。

    3.3 少量氮肥施加結(jié)合秸稈還田促進(jìn)土壤養(yǎng)分的累積

    在秸稈還田的基礎(chǔ)上,土壤中總氮、堿解氮含量隨著氮肥的減施呈波動(dòng)上升,而硝態(tài)氮和銨態(tài)氮含量波動(dòng)下降。這可能是由于少量的氮肥會(huì)促進(jìn)土壤的固氮作用,大量的氮肥反而會(huì)起抑制作用。一部分氮會(huì)被土壤團(tuán)聚體或有機(jī)質(zhì)吸附,降低氮的流失速率,從而增加土壤的氮含量;少量的氮肥可以促進(jìn)土壤中氮循環(huán)過(guò)程,增強(qiáng)土壤中氮的微生物分解作用,本研究中硝化細(xì)菌的豐度在少量施加氮肥處理中也明顯增多(圖3),這些微生物能夠分解秸稈,氮有效性通過(guò)優(yōu)先改變秸稈碳同化微生物群落介導(dǎo)土壤有機(jī)質(zhì)的啟動(dòng)效應(yīng)。此外,少量的氮肥可能導(dǎo)致土壤pH值的變化,使土壤中的氮更容易形成穩(wěn)定的化合物,不容易流失,從而增加土壤氮含量[31。

    在秸稈還田基礎(chǔ)上,在少量施加氮肥的處理中土壤有效鉀含量和土壤pH值達(dá)到最高值。這是由于秸稈還田可改善土壤結(jié)構(gòu),增加土壤的有機(jī)質(zhì)含量,改善土壤的質(zhì)地,這有助于提高土壤的保水保肥能力,有利于磷、鉀等營(yíng)養(yǎng)物質(zhì)的積累。同時(shí),有機(jī)質(zhì)含量與土壤中鉀含量顯著相關(guān)(圖5),一方面秸稈還田通過(guò)增加大團(tuán)聚體中的土壤腐殖酸來(lái)提高對(duì)鉀的吸附;另一方面,秸稈的摻入能夠在一定程度上提高土壤中磷鉀的有效性,使其更易被植物吸收[32

    施用少量氮肥可以為作物生長(zhǎng)提供必要的營(yíng)養(yǎng)物質(zhì),促進(jìn)植物對(duì)磷和鉀的吸收,提高植物對(duì)這些養(yǎng)分的利用效率,有助于土壤中磷鉀的累積。此外,施加少量的氮肥可促進(jìn)微生物活動(dòng),土壤中的微生物會(huì)參與分解釋放出養(yǎng)分,包括磷和鉀,將這些養(yǎng)分轉(zhuǎn)化為植物可吸收的形式。本試驗(yàn)結(jié)果表明,細(xì)菌主要通過(guò)調(diào)節(jié)土壤中鉀含量,而真菌主要通過(guò)調(diào)節(jié)磷含量來(lái)改善土壤質(zhì)量(圖5)??傊诮斩掃€田的基礎(chǔ)上應(yīng)當(dāng)注意適量施用氮肥,避免過(guò)量氮肥造成負(fù)面影響,以實(shí)現(xiàn)土壤養(yǎng)分、pH值的合理調(diào)控。

    綜上所述,在秸稈還田基礎(chǔ)上,少量化肥或者不使用化肥能夠提高土壤細(xì)菌和真菌的多樣性和微生物之間的聯(lián)系,促進(jìn)優(yōu)勢(shì)屬和關(guān)鍵屬發(fā)揮降解污染、促進(jìn)植物生長(zhǎng)、氮素循環(huán)等功能,提高土壤磷鉀的含量。因此,減少氮肥施加是提高土壤微生物多樣性,改善土壤質(zhì)量的有效手段。本試驗(yàn)結(jié)果可為秸稈還田替代化肥技術(shù)的應(yīng)用提供技術(shù)支持,為化肥減施的農(nóng)業(yè)可持續(xù)發(fā)展提供理論依據(jù)。

    本研究發(fā)現(xiàn),東北黑土區(qū)水稻秸稈條件下,不同化肥添加量對(duì)土壤微生物多樣性的影響并不一致。不同化肥添加量對(duì)土壤細(xì)菌α多樣性的影響要顯著高于真菌α多樣性,但均顯著改變了土壤細(xì)菌和真菌β多樣性。不同化肥施加量顯著改變了土壤全氮、全磷、有效鉀的含量,但是對(duì)其他理化性質(zhì)影響不顯著。不同化肥添加量并不會(huì)導(dǎo)致土壤細(xì)菌和真菌的優(yōu)勢(shì)菌發(fā)生改變,但會(huì)顯著改變其相對(duì)豐度,細(xì)菌門(mén)中酸桿菌門(mén)、放線菌門(mén)和擬桿菌門(mén)相對(duì)豐度最高;真菌門(mén)中擔(dān)子菌門(mén)、子囊菌門(mén)和壺菌門(mén)相對(duì)豐度最高。土壤細(xì)菌群落組成和多樣性主要受到土壤全氮、硝態(tài)氮、總鉀和速效鉀含量的影響,土壤真菌群落結(jié)構(gòu)和多樣性主要受到土壤pH值的影響。綜上,化肥減施下秸稈還田會(huì)改變土壤微生物的結(jié)構(gòu)和多樣性,而常規(guī)水稻秸稈還田可以提高土壤微生物的多樣性,對(duì)保護(hù)土壤的健康和穩(wěn)定具有積極意義。

    參考文獻(xiàn):

    [1]擺小蓉,閔煒?lè)?,石亞飛,等. 干旱對(duì)不同水稻萌發(fā)種子中逆境生理及相應(yīng)基因表達(dá)的影響[J]. 華北農(nóng)學(xué)報(bào),2023,38(4):101-111.

    [2]盧鈺霞,張?jiān)倬?,田志? 關(guān)聯(lián)分析在水稻產(chǎn)量性狀遺傳研究中的應(yīng)用[J]. 浙江農(nóng)業(yè)科學(xué),2021,62(7):1293-1298.

    [3]吳玉德,張 旭,關(guān)法春,等. 寒地水稻秸稈還田對(duì)土壤生態(tài)環(huán)境影響的研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué),2023,51(15):1-8.

    [4]王俊芹,苑甜甜. 中國(guó)農(nóng)業(yè)綠色發(fā)展政策演進(jìn)及政策工具分析[J]. 河北學(xué)刊,2023,43(2):130-139.

    [5]Linquist B A,Adviento-Borbe M A,Pittelkow C M,et al. Fertilizer management practices and greenhouse gas emissions from rice systems:a quantitative review and analysis[J]. Field Crops Research,2012,135:10-21.

    [6]Yao Z S,Zheng X H,Dong H B,et al. A 3-year record of N2O and CH4 emissions from a sandy loam paddy during rice seasons as affected by different nitrogen application rates[J]. Agriculture,Ecosystems & Environment,2012,152:1-9.

    [7]文高輝,王夏玙,胡賢輝,等. 耕地規(guī)模、細(xì)碎化對(duì)化肥面源污染的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文),2023,31(7):1081-1093.

    [8]Islam M R,Singh Chauhan P,Kim Y,et al. Community level functional diversity and enzyme activities in paddy soils under different long-term fertilizer management practices[J]. Biology and Fertility of Soils,2011,47(5):599-604.

    [9]Kumar A,Yadav D S. Long-term effects of fertilizers on the soil fertility and productivity of a rice-wheat system[J]. Journal of Agronomy and Crop Science,2001,186(1):47-54.

    [10]姚 遠(yuǎn),馬 泉,鄭國(guó)利,等. 麥/稻秸稈還田氮在后茬稻/麥植株中的分配去向探析[J]. 麥類(lèi)作物學(xué)報(bào),2023,43(8):1047-1055.

    [11]譚 娟,陳 楠,董 偉,等. 玉米秸稈還田量對(duì)小麥生理生態(tài)特征和產(chǎn)量的影響[J]. 安徽農(nóng)業(yè)科學(xué),2019,47(18):41-42,45.

    [12]劉安凱,史登林,王小利,等. 秸稈和生物炭還田對(duì)稻田土壤有機(jī)碳礦化和水稻產(chǎn)質(zhì)量的影響[J]. 山地農(nóng)業(yè)生物學(xué)報(bào),2021,40(4):38-45.

    [13]章力干,石心怡,王玉寶,等. 秸稈還田對(duì)中國(guó)主要糧食作物病害影響的Meta分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(21):93-100.

    [14]吳立鵬,張士榮,婁金華,等. 秸稈還田與優(yōu)化施氮對(duì)稻田土壤碳氮含量及產(chǎn)量的影響[J]. 華北農(nóng)學(xué)報(bào),2019,34(4):158-166.

    [15]吳玉紅,郝興順,田霄鴻,等. 秸稈還田與化肥減量配施對(duì)稻茬麥土壤養(yǎng)分、酶活性及產(chǎn)量影響[J]. 西南農(nóng)業(yè)學(xué)報(bào),2018,31(5):998-1005.

    [16]Hu X J,Gu H D,Liu J J,et al. High variation of fungal communities and associated potential plant pathogens induced by long-term addition of N fertilizers rather than P,K fertilizers:a case study in a Mollisol field[J]. Soil Ecology Letters,2022,4(4):348-361.

    [17]YaoQ,Liu J J,Yu Z H,et al. Changes of bacterial community compositions after three years of biochar application in a black soil of northeast China[J]. Applied Soil Ecology,2017,113:11-21.

    [18]Huang R,Gao M,Wang Y Y,et al. Response of the soil N2O emission and ammonia-oxidizing microorganism community to the maize straw return with reducing fertilizer in purple soil[J]. Environmental Science,2019,40(1):401-411.

    [19]楊濱娟,黃國(guó)勤,錢(qián)海燕. 秸稈還田配施化肥對(duì)土壤溫度、根際微生物及酶活性的影響[J]. 土TM2lMVAkHphEEJtaMKz/v7tOH9xbZsF1DZCu7/uOrYo=壤學(xué)報(bào),2014,51(1):150-157.

    [20]韓曉增,鄒文秀. 東北黑土地保護(hù)利用研究足跡與科技研發(fā)展望[J]. 土壤學(xué)報(bào),2021,58(6):1341-1358.

    [21]孟祥宇,冉 成,李金平,等. 秸稈還田對(duì)東北黑土稻區(qū)土壤養(yǎng)分及水稻產(chǎn)量的影響[J]. 東北農(nóng)業(yè)科學(xué),2023,48(2):64-67,89.

    [22]張宇杰,王志強(qiáng),馬 鵬,等. 麥稈還田下水氮耦合對(duì)水稻氮素吸收利用及產(chǎn)量的影響[J]. 中國(guó)水稻科學(xué),2022,36(4):388-398.

    [23]孫園園,張 橋,孫永健,等. 秸稈還田與水氮管理對(duì)水稻氮素利用及土壤理化性質(zhì)的影響[J]. 湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,48(1):65-74.

    [24]赫 兵,李 超,嚴(yán)永峰,等. 水稻秸稈秋季水耙漿還田對(duì)土壤及水稻性狀的影響[J]. 中國(guó)稻米,2023,29(4):66-71.

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

    [26]陳凱鵬,韋菊嫻,陳 丹,等. 化肥減量配施秸稈對(duì)雙季稻田固氮微生物群落的影響[J]. 環(huán)境科學(xué),2022,43(10):4725-4734.

    [27]黃 容,高 明,王鎣燕,等. 紫色土N2O排放及氨氧化微生物群落結(jié)構(gòu)對(duì)玉米秸稈與化肥減量配施的響應(yīng)[J]. 環(huán)境科學(xué),2019,40(1):401-411.

    [28]Zhu C,Ling N,Guo J J,et al. Impacts of fertilization regimes on arbuscular mycorrhizal fungal (AMF) community composition were correlated with organic matter composition in maize rhizosphere soil[J]. Frontiers in Microbiology,2016,7:1840.

    [29]馬玉穎,張煥朝,項(xiàng)興佳,等. 長(zhǎng)期施肥對(duì)砂姜黑土叢枝菌根真菌群落的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào),2018,29(10):3398-3406.

    [30]伍龍梅,張 彬,黃 慶,等. 華南稻區(qū)不同水稻品種耐淹出苗能力及秧苗素質(zhì)比較[J]. 廣東農(nóng)業(yè)科學(xué),2021,48(10):141-147.

    [31]趙子婧,孫建平,戴相林,等. 秸稈還田結(jié)合減量施肥對(duì)水稻產(chǎn)量和土壤養(yǎng)分的影響[J]. 江蘇農(nóng)業(yè)科學(xué),2022,50(10):66-71.

    [32]Yuan L,Zhang Z C,Cao X C,et al. Responses of rice production,milled rice quality and soil properties to various nitrogen inputs and rice straw incorporation under continuous plastic film mulching cultivation[J]. Field Crops Research,2014,155:164-171.

    国产一区亚洲一区在线观看| 精品一品国产午夜福利视频| 国产精品久久久av美女十八| av在线app专区| 色94色欧美一区二区| 国产伦理片在线播放av一区| 国产精品麻豆人妻色哟哟久久| 嫩草影视91久久| 国产日韩欧美亚洲二区| 国产xxxxx性猛交| 日韩电影二区| av视频免费观看在线观看| 精品一品国产午夜福利视频| 1024香蕉在线观看| 日本黄色日本黄色录像| 国产在线视频一区二区| 人人妻人人澡人人爽人人夜夜| 成年人免费黄色播放视频| 国产免费现黄频在线看| 亚洲第一区二区三区不卡| 久久免费观看电影| 久久久久久久久久久久大奶| 亚洲精品久久久久久婷婷小说| 99精品久久久久人妻精品| av在线观看视频网站免费| 亚洲美女视频黄频| 成人黄色视频免费在线看| 桃花免费在线播放| 日本猛色少妇xxxxx猛交久久| 欧美日韩视频高清一区二区三区二| netflix在线观看网站| 水蜜桃什么品种好| 99热全是精品| 韩国精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美黄色片欧美黄色片| 最近中文字幕2019免费版| 日韩av在线免费看完整版不卡| 操出白浆在线播放| 亚洲av福利一区| 成人漫画全彩无遮挡| 亚洲综合精品二区| 免费黄色在线免费观看| 老司机靠b影院| 一级爰片在线观看| 国产成人一区二区在线| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲av国产电影网| 亚洲国产av影院在线观看| 亚洲精品国产色婷婷电影| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 中文字幕制服av| 欧美av亚洲av综合av国产av | 亚洲精品视频女| 在线 av 中文字幕| 国产精品一区二区精品视频观看| 美女视频免费永久观看网站| 精品亚洲成国产av| 欧美老熟妇乱子伦牲交| 在线天堂最新版资源| 亚洲熟女精品中文字幕| 狂野欧美激情性xxxx| 国产精品一区二区在线不卡| 精品国产国语对白av| 国产亚洲av高清不卡| 伊人久久国产一区二区| 国产不卡av网站在线观看| 十八禁网站网址无遮挡| 亚洲精品aⅴ在线观看| 观看av在线不卡| av网站在线播放免费| 亚洲成人手机| 日韩免费高清中文字幕av| 国产乱来视频区| tube8黄色片| 高清视频免费观看一区二区| 中文天堂在线官网| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 高清不卡的av网站| 久久久久国产精品人妻一区二区| www.自偷自拍.com| 久久久精品免费免费高清| 90打野战视频偷拍视频| 久久久久久久久久久久大奶| av网站在线播放免费| 亚洲四区av| 少妇精品久久久久久久| videosex国产| 免费高清在线观看视频在线观看| 日韩av免费高清视频| 欧美日韩视频高清一区二区三区二| 国产成人精品久久二区二区91 | 一级片免费观看大全| 亚洲色图 男人天堂 中文字幕| 各种免费的搞黄视频| 丁香六月欧美| 亚洲三区欧美一区| 国产日韩欧美在线精品| 不卡视频在线观看欧美| 男女之事视频高清在线观看 | 人体艺术视频欧美日本| 一区二区三区精品91| 十分钟在线观看高清视频www| 久久99热这里只频精品6学生| 国产精品久久久久久精品电影小说| 99久久人妻综合| 国产探花极品一区二区| 天天躁夜夜躁狠狠久久av| 九色亚洲精品在线播放| 少妇人妻精品综合一区二区| 一个人免费看片子| av有码第一页| 十八禁人妻一区二区| 色综合欧美亚洲国产小说| 丰满乱子伦码专区| 99久久人妻综合| 五月开心婷婷网| 波多野结衣一区麻豆| 亚洲精品国产一区二区精华液| 亚洲成色77777| 亚洲av在线观看美女高潮| 国产视频首页在线观看| 男人舔女人的私密视频| 国产在线免费精品| 精品国产一区二区三区久久久樱花| 男女床上黄色一级片免费看| 两性夫妻黄色片| 欧美在线一区亚洲| 亚洲欧洲国产日韩| 在线观看国产h片| 欧美国产精品va在线观看不卡| www.自偷自拍.com| 老汉色av国产亚洲站长工具| 日韩一本色道免费dvd| 国产精品亚洲av一区麻豆 | 岛国毛片在线播放| 亚洲精品自拍成人| av又黄又爽大尺度在线免费看| netflix在线观看网站| 国产精品蜜桃在线观看| 蜜桃在线观看..| 高清黄色对白视频在线免费看| 亚洲免费av在线视频| 91精品伊人久久大香线蕉| 丝袜脚勾引网站| 王馨瑶露胸无遮挡在线观看| 日本欧美国产在线视频| 一级a爱视频在线免费观看| 在线观看免费视频网站a站| 久久鲁丝午夜福利片| 亚洲图色成人| 国产精品二区激情视频| 交换朋友夫妻互换小说| 极品人妻少妇av视频| 日本av手机在线免费观看| 少妇人妻精品综合一区二区| 韩国高清视频一区二区三区| 国产精品一二三区在线看| 久久精品aⅴ一区二区三区四区| 男女无遮挡免费网站观看| 国产日韩欧美在线精品| 99精国产麻豆久久婷婷| 人体艺术视频欧美日本| 欧美成人精品欧美一级黄| 亚洲视频免费观看视频| 在线观看免费视频网站a站| 制服丝袜香蕉在线| 一级毛片黄色毛片免费观看视频| 国产极品粉嫩免费观看在线| 亚洲色图 男人天堂 中文字幕| 成年女人毛片免费观看观看9 | a 毛片基地| 日韩 亚洲 欧美在线| 日韩 欧美 亚洲 中文字幕| 久久精品熟女亚洲av麻豆精品| a 毛片基地| 国产一区二区三区av在线| 国产97色在线日韩免费| 在线观看免费视频网站a站| 校园人妻丝袜中文字幕| 一区二区av电影网| 1024香蕉在线观看| 亚洲国产欧美网| 丰满乱子伦码专区| 免费高清在线观看日韩| 波多野结衣一区麻豆| 精品少妇久久久久久888优播| 亚洲欧美激情在线| 另类亚洲欧美激情| 亚洲激情五月婷婷啪啪| 美女扒开内裤让男人捅视频| 不卡av一区二区三区| tube8黄色片| 天堂中文最新版在线下载| e午夜精品久久久久久久| 亚洲成人手机| 乱人伦中国视频| 最近中文字幕2019免费版| 亚洲免费av在线视频| 一本大道久久a久久精品| 成年人午夜在线观看视频| 十八禁高潮呻吟视频| 啦啦啦 在线观看视频| 高清黄色对白视频在线免费看| 精品午夜福利在线看| xxx大片免费视频| 国产一区亚洲一区在线观看| 2021少妇久久久久久久久久久| av在线老鸭窝| xxxhd国产人妻xxx| 王馨瑶露胸无遮挡在线观看| 久久久久久久国产电影| 天天躁狠狠躁夜夜躁狠狠躁| 波野结衣二区三区在线| 免费在线观看黄色视频的| 午夜老司机福利片| 久久鲁丝午夜福利片| 在线天堂最新版资源| 五月开心婷婷网| 热re99久久精品国产66热6| 久久97久久精品| 毛片一级片免费看久久久久| 亚洲国产精品一区三区| 最新在线观看一区二区三区 | 国产精品一二三区在线看| 可以免费在线观看a视频的电影网站 | 欧美精品亚洲一区二区| 99九九在线精品视频| 欧美少妇被猛烈插入视频| 亚洲精品乱久久久久久| 大话2 男鬼变身卡| 亚洲精品美女久久久久99蜜臀 | 亚洲av中文av极速乱| 国产成人一区二区在线| 天天躁夜夜躁狠狠躁躁| 亚洲精品成人av观看孕妇| e午夜精品久久久久久久| 亚洲国产欧美一区二区综合| 不卡视频在线观看欧美| 国产有黄有色有爽视频| 成人手机av| 两性夫妻黄色片| 国产精品麻豆人妻色哟哟久久| 国产亚洲一区二区精品| 18禁裸乳无遮挡动漫免费视频| 婷婷色麻豆天堂久久| 国产精品一二三区在线看| 亚洲成人一二三区av| 亚洲欧美日韩另类电影网站| 国产片特级美女逼逼视频| 亚洲国产精品成人久久小说| 亚洲免费av在线视频| 汤姆久久久久久久影院中文字幕| 国产亚洲欧美精品永久| 欧美人与性动交α欧美精品济南到| 少妇的丰满在线观看| 青春草国产在线视频| av福利片在线| 肉色欧美久久久久久久蜜桃| 亚洲国产看品久久| 亚洲自偷自拍图片 自拍| 久久久精品94久久精品| 国产精品久久久人人做人人爽| 少妇 在线观看| 精品一区二区免费观看| 欧美精品一区二区免费开放| 成人毛片60女人毛片免费| 丁香六月天网| 国产成人精品无人区| 久久久国产精品麻豆| 哪个播放器可以免费观看大片| 国产av码专区亚洲av| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说| 亚洲欧美激情在线| 亚洲,欧美,日韩| 1024香蕉在线观看| 搡老岳熟女国产| 国产 一区精品| 亚洲精品在线美女| 日日啪夜夜爽| 国产精品一二三区在线看| 一边摸一边抽搐一进一出视频| 久久精品亚洲熟妇少妇任你| 午夜日韩欧美国产| 观看美女的网站| 国产精品二区激情视频| 九色亚洲精品在线播放| 国产精品久久久久成人av| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产av蜜桃| 亚洲国产最新在线播放| videos熟女内射| 亚洲男人天堂网一区| av.在线天堂| 国产精品一二三区在线看| 久久精品久久久久久噜噜老黄| 欧美日韩视频精品一区| 欧美人与善性xxx| 日本午夜av视频| 久久久久久久久免费视频了| 中文字幕精品免费在线观看视频| 男人添女人高潮全过程视频| 国产男女内射视频| 欧美久久黑人一区二区| 欧美另类一区| 欧美成人午夜精品| 国产精品人妻久久久影院| 精品少妇内射三级| 精品亚洲乱码少妇综合久久| 久久女婷五月综合色啪小说| 香蕉国产在线看| 综合色丁香网| 999久久久国产精品视频| 欧美在线一区亚洲| 黄色一级大片看看| 这个男人来自地球电影免费观看 | 你懂的网址亚洲精品在线观看| 亚洲国产av影院在线观看| 亚洲成人免费av在线播放| 国产亚洲午夜精品一区二区久久| 热re99久久国产66热| 亚洲国产欧美日韩在线播放| 国产精品一区二区精品视频观看| 久久99热这里只频精品6学生| 国产成人a∨麻豆精品| 丰满乱子伦码专区| 久久毛片免费看一区二区三区| 欧美xxⅹ黑人| 国产成人精品福利久久| 国产精品国产av在线观看| 欧美激情高清一区二区三区 | 国产免费福利视频在线观看| 免费黄网站久久成人精品| 丝瓜视频免费看黄片| 亚洲四区av| 久久天堂一区二区三区四区| 久久韩国三级中文字幕| 不卡av一区二区三区| 青春草国产在线视频| 又粗又硬又长又爽又黄的视频| 国产精品无大码| 久久久久久人妻| 天天躁狠狠躁夜夜躁狠狠躁| 国产av一区二区精品久久| 国产精品香港三级国产av潘金莲 | 国产福利在线免费观看视频| 人体艺术视频欧美日本| 肉色欧美久久久久久久蜜桃| 久久热在线av| 啦啦啦视频在线资源免费观看| 1024视频免费在线观看| 国产成人欧美在线观看 | 天天躁夜夜躁狠狠躁躁| 高清视频免费观看一区二区| 亚洲五月色婷婷综合| 久久av网站| 99re6热这里在线精品视频| 午夜激情久久久久久久| 2018国产大陆天天弄谢| 不卡av一区二区三区| 亚洲国产av影院在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品一二三| 免费看不卡的av| 国产精品蜜桃在线观看| 国产精品一区二区在线不卡| 一区二区三区激情视频| 欧美人与性动交α欧美软件| 亚洲精品日本国产第一区| 七月丁香在线播放| 久久av网站| a级毛片在线看网站| 交换朋友夫妻互换小说| 亚洲色图 男人天堂 中文字幕| 一区二区日韩欧美中文字幕| 国产精品国产av在线观看| 亚洲成人免费av在线播放| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 七月丁香在线播放| 日韩av不卡免费在线播放| 波多野结衣av一区二区av| 国产一卡二卡三卡精品 | 多毛熟女@视频| 亚洲在久久综合| 亚洲av成人精品一二三区| 亚洲av电影在线观看一区二区三区| 国产精品免费大片| 黄频高清免费视频| 极品少妇高潮喷水抽搐| 久久韩国三级中文字幕| av.在线天堂| 国产成人精品久久久久久| 国产一区亚洲一区在线观看| 男女高潮啪啪啪动态图| 婷婷色综合www| 免费黄色在线免费观看| 久久久国产欧美日韩av| 麻豆乱淫一区二区| 亚洲 欧美一区二区三区| 久久久久久久精品精品| 又粗又硬又长又爽又黄的视频| 18在线观看网站| 久久久久久久久久久免费av| 婷婷色av中文字幕| netflix在线观看网站| 高清在线视频一区二区三区| 一个人免费看片子| 嫩草影院入口| 中文乱码字字幕精品一区二区三区| 色播在线永久视频| 好男人视频免费观看在线| 国产野战对白在线观看| 欧美日韩亚洲高清精品| 精品免费久久久久久久清纯 | 久久久亚洲精品成人影院| 中国三级夫妇交换| 国产精品久久久久久人妻精品电影 | 亚洲欧美日韩另类电影网站| 99久久99久久久精品蜜桃| 欧美黑人欧美精品刺激| 日韩av免费高清视频| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| av福利片在线| 久久影院123| 一边摸一边做爽爽视频免费| 亚洲成人手机| 狂野欧美激情性bbbbbb| 久久久久久久精品精品| 精品久久久久久电影网| 黄色视频在线播放观看不卡| 国产精品 国内视频| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 日本一区二区免费在线视频| 夜夜骑夜夜射夜夜干| 97人妻天天添夜夜摸| 精品一区二区三卡| 天天躁夜夜躁狠狠久久av| 一级毛片电影观看| 少妇人妻久久综合中文| 各种免费的搞黄视频| 不卡视频在线观看欧美| 老熟女久久久| 国产精品免费视频内射| 国产深夜福利视频在线观看| av天堂久久9| 卡戴珊不雅视频在线播放| 久久久欧美国产精品| 国产欧美亚洲国产| 国产成人精品久久二区二区91 | 男女下面插进去视频免费观看| 男女床上黄色一级片免费看| 性色av一级| 国产精品嫩草影院av在线观看| 最近的中文字幕免费完整| 侵犯人妻中文字幕一二三四区| 两个人免费观看高清视频| www.av在线官网国产| 美女主播在线视频| 狠狠婷婷综合久久久久久88av| 最近中文字幕2019免费版| 青春草国产在线视频| 亚洲欧美一区二区三区久久| 中文字幕高清在线视频| 亚洲精品日本国产第一区| 美女中出高潮动态图| 国产精品久久久久久人妻精品电影 | 成人亚洲欧美一区二区av| 男女下面插进去视频免费观看| 国产男人的电影天堂91| 一区二区三区四区激情视频| 无遮挡黄片免费观看| 大陆偷拍与自拍| 国产视频首页在线观看| 老熟女久久久| 777米奇影视久久| 亚洲av国产av综合av卡| 国产xxxxx性猛交| 久久久久精品性色| av有码第一页| 高清av免费在线| 一本大道久久a久久精品| 九色亚洲精品在线播放| 女人久久www免费人成看片| 久久99一区二区三区| 国产熟女午夜一区二区三区| 久久ye,这里只有精品| 桃花免费在线播放| 欧美另类一区| 免费日韩欧美在线观看| 人人妻人人澡人人看| 超碰成人久久| 婷婷色av中文字幕| 99re6热这里在线精品视频| 51午夜福利影视在线观看| 亚洲熟女毛片儿| 欧美97在线视频| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 欧美亚洲日本最大视频资源| 一级黄片播放器| 美女中出高潮动态图| 欧美 亚洲 国产 日韩一| 80岁老熟妇乱子伦牲交| 国产精品免费大片| 91精品国产国语对白视频| 亚洲成人手机| 亚洲 欧美一区二区三区| 热re99久久国产66热| 麻豆乱淫一区二区| 悠悠久久av| 黑人巨大精品欧美一区二区蜜桃| 另类亚洲欧美激情| 香蕉丝袜av| 亚洲精品久久久久久婷婷小说| 岛国毛片在线播放| 久久综合国产亚洲精品| 啦啦啦在线观看免费高清www| 麻豆乱淫一区二区| 久久亚洲国产成人精品v| 精品国产乱码久久久久久男人| 免费在线观看视频国产中文字幕亚洲 | 另类亚洲欧美激情| 欧美日韩一级在线毛片| 亚洲五月色婷婷综合| 黑人猛操日本美女一级片| 免费在线观看完整版高清| 免费观看a级毛片全部| 一本色道久久久久久精品综合| 亚洲国产最新在线播放| 国产无遮挡羞羞视频在线观看| 亚洲精品一区蜜桃| 中文字幕人妻熟女乱码| 久久久久久久大尺度免费视频| 国产成人精品久久二区二区91 | 桃花免费在线播放| 亚洲婷婷狠狠爱综合网| 欧美黄色片欧美黄色片| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲日产国产| av卡一久久| 狂野欧美激情性bbbbbb| 精品一区二区免费观看| 亚洲色图综合在线观看| 我要看黄色一级片免费的| 一本一本久久a久久精品综合妖精| 纵有疾风起免费观看全集完整版| 伦理电影大哥的女人| 国产麻豆69| 18禁动态无遮挡网站| 欧美黑人欧美精品刺激| 亚洲精品aⅴ在线观看| 观看av在线不卡| 啦啦啦视频在线资源免费观看| 久久精品aⅴ一区二区三区四区| 国产一区亚洲一区在线观看| 最新在线观看一区二区三区 | 精品第一国产精品| av在线app专区| 一区二区三区精品91| 母亲3免费完整高清在线观看| 777久久人妻少妇嫩草av网站| 亚洲精华国产精华液的使用体验| 亚洲第一青青草原| 亚洲一码二码三码区别大吗| 少妇的丰满在线观看| 中文字幕制服av| 丝袜美腿诱惑在线| 制服人妻中文乱码| 成年av动漫网址| 男女下面插进去视频免费观看| 如日韩欧美国产精品一区二区三区| 一本一本久久a久久精品综合妖精| 看十八女毛片水多多多| 满18在线观看网站| 日本wwww免费看| 免费av中文字幕在线| 亚洲欧美中文字幕日韩二区| 妹子高潮喷水视频| 成年动漫av网址| 免费高清在线观看视频在线观看| 亚洲一区中文字幕在线| 日韩 欧美 亚洲 中文字幕| 日韩人妻精品一区2区三区| 看免费av毛片| 美国免费a级毛片| 亚洲精品中文字幕在线视频| 亚洲av日韩精品久久久久久密 | 视频区图区小说| 国产精品 国内视频| 亚洲婷婷狠狠爱综合网| 国产一区有黄有色的免费视频| 大香蕉久久成人网| 久久久久久人人人人人| 久久精品亚洲av国产电影网| 亚洲精品第二区| 国产高清不卡午夜福利| 亚洲 欧美一区二区三区| 老熟女久久久| 亚洲精品美女久久久久99蜜臀 | 国产av一区二区精品久久| 免费看不卡的av| 国产精品 欧美亚洲| 三上悠亚av全集在线观看| 丝袜在线中文字幕| 国产精品一区二区在线观看99| 天堂8中文在线网| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区av在线| 久久精品亚洲熟妇少妇任你|