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

    不同施肥處理對玉米-小麥輪作土壤微生物群落功能多樣性的影響

    2016-12-24 05:31:34邢鵬飛武曉森高圣超李洪杰趙同凱周曉琳沈德龍孫軍德
    微生物學(xué)雜志 2016年1期
    關(guān)鍵詞:均勻度氮磷碳源

    邢鵬飛, 武曉森, 高圣超, 李洪杰, 趙同凱, 周曉琳, 沈德龍*, 孫軍德

    (1.沈陽農(nóng)業(yè)大學(xué) 土地與環(huán)境學(xué)院,遼寧 沈陽 110866;2.中國農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081;3.農(nóng)業(yè)部微生物產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室,北京 100081;4.河北農(nóng)業(yè)大學(xué) 生命科學(xué)學(xué)院,河北 保定 071000;5.德州市農(nóng)業(yè)科學(xué)研究院,山東 德州 253015)

    ?

    不同施肥處理對玉米-小麥輪作土壤微生物群落功能多樣性的影響

    邢鵬飛1,2,3, 武曉森2,3,4, 高圣超2,3, 李洪杰5, 趙同凱5, 周曉琳5, 沈德龍2,3*, 孫軍德1*

    (1.沈陽農(nóng)業(yè)大學(xué) 土地與環(huán)境學(xué)院,遼寧 沈陽 110866;2.中國農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081;3.農(nóng)業(yè)部微生物產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室,北京 100081;4.河北農(nóng)業(yè)大學(xué) 生命科學(xué)學(xué)院,河北 保定 071000;5.德州市農(nóng)業(yè)科學(xué)研究院,山東 德州 253015)

    土壤微生物多樣性能反應(yīng)土壤的肥力,不同的施肥措施對土壤微生物的種群和功能多樣性也會產(chǎn)生重要的影響。以山東德州連續(xù)兩年小麥季和玉米季收獲后土壤為研究對象,利用Biolog技術(shù)研究了6種不同施肥處理對土壤微生物群落功能多樣性的影響。結(jié)果表明:其中各個施肥處理的平均顏色變化率(average well color development,AWCD)差異顯著,常規(guī)氮磷鉀肥+全量秸稈還田+秸稈腐熟劑(FS)處理代謝活性最高;物種豐富度指數(shù)(H)和均勻度指數(shù)(E)也表明各施肥方式均能夠維持微生物種群的多樣性,其中FS和30%豬糞+70%常規(guī)氮磷鉀肥(OF)處理物種豐富度指數(shù)(H)和均勻度指數(shù)(E)最高;PCA及RDA分析顯示,OF和FS處理微生物功能多樣性相似,且其微生物功能多樣性與有機(jī)質(zhì)(Soil organic matter,SOM)、全氮(Total N,TN)、速效磷(Available P,AP)和速效鉀(Available K,AK)密切相關(guān)。豬糞堆肥有機(jī)無機(jī)復(fù)合肥3 600 kg/hm2(OI2)處理與豬糞堆肥有機(jī)無機(jī)復(fù)合肥1 800 kg/hm2(OI1)處理相似,其功能多樣性比常規(guī)施肥(CF)處理稍高。綜上所述,OF處理和FS處理的土壤微生物群落功能多樣性程度高于其他處理,說明秸稈還田+秸稈腐熟劑和有機(jī)肥部分替代氮磷鉀肥能夠顯著提高土壤微生物功能多樣性,有利于保護(hù)土壤微生態(tài)。

    不同施肥處理;土壤微生物;群落多樣性;功能多樣性;Biolog

    土壤微生物是土壤的重要組成部分,更是土壤生態(tài)系統(tǒng)重要的組成部分[1]。土壤微生物是土壤生態(tài)系統(tǒng)養(yǎng)分循環(huán)和能量流動的主導(dǎo)者,對于維持土壤生態(tài)系統(tǒng)的穩(wěn)定性和土壤的可持續(xù)利用具有重要作用[2]。土壤微生物可以調(diào)節(jié)土壤微生態(tài),延長礦質(zhì)元素的物質(zhì)循環(huán)途徑,促進(jìn)營養(yǎng)元素的多元化,從而促進(jìn)植物對營養(yǎng)元素的吸收,利于土壤肥力的保持[3-6]。與此同時,土壤中的一系列非生物的物化反應(yīng)也對土壤微生物的種類及分布特征[7]產(chǎn)生深度的影響。施肥作為農(nóng)業(yè)生產(chǎn)中的重要一環(huán),對土壤和作物有重要的影響,采取不同的耕作制度或施用不同的肥料對土壤微生物數(shù)量、群落和功能也產(chǎn)生不同的影響[8]。施肥使農(nóng)田微生物群落結(jié)構(gòu)和功能發(fā)生了顯著的變化,研究不同施肥方式對農(nóng)田養(yǎng)分循環(huán)和施肥與環(huán)境間的關(guān)系有重要意義[9]。王銳等[10]認(rèn)為,長期施用有機(jī)肥更有利于提高土壤肥力,促進(jìn)土壤微生物生長繁育。王才斌等[11]研究表明,增施有機(jī)肥能有效增加土壤中細(xì)菌、放線菌和真菌數(shù)量,而在單施化肥的情況下,微生物數(shù)量不增加或者增加不明顯。張瑞等[12]研究表明,單施有機(jī)肥的處理土壤微生物群落碳源利用能力明顯高于其他處理。劉振香等[13]研究表明,施用有機(jī)肥可以顯著提高土壤微生物數(shù)量、多樣性及豐度,降低細(xì)菌群落優(yōu)勢度。還有研究表明,過多施用有機(jī)肥,并不能無限提高土壤微生物的數(shù)量[14],但是對微生物活性有明顯的增強(qiáng)效果。秸稈還田等施肥措施可以增加土壤有機(jī)質(zhì)的含量,進(jìn)而為微生物提供了良好的物質(zhì)基礎(chǔ),改變了土壤微生物的多樣性[15]。Biolog微孔板技術(shù)是通過測定土壤微生物對31種不同單一碳源底物的代謝活性差異,來表征微生物種群多樣性,并以此確定功能微生物種群變化的指標(biāo)[16],這種方法可比單純檢測微生物總數(shù)能更好地反映出土壤微生物種群功能的多樣性。本研究采用Biolog-Eco板,對黃淮海地區(qū)連續(xù)兩年冬小麥-夏玉米輪作區(qū)土壤微生物群落功能多樣性進(jìn)行了研究,旨在深入研究增施不同量有機(jī)肥及秸稈還田對土壤微生物的影響,為進(jìn)一步獲得良好的肥料施用方式提供理論依據(jù)。

    1 材料與方法

    1.1 材料

    1.1.1 試驗田概況 該試驗田地處山東德州市哨馬營村,位于東經(jīng)115°45′~117°24′、北緯36°24′~38°0′,屬典型的溫帶大陸性季風(fēng)氣候,年均氣溫12.9 ℃,年均降水量547.5 mm。試驗農(nóng)田為黃河沖積平原,作物種植方式為典型的冬小麥-夏玉米輪作,小麥種植品種為濟(jì)麥22,玉米為鄭單958。試驗前農(nóng)田土壤養(yǎng)分含量:全氮1.2 g/kg,全磷0.83 g/kg,全鉀20.8 g/kg,有效磷37.4 mg/kg,有機(jī)質(zhì)含量20.3 g/kg,速效鉀282 mg/kg,pH 8.3。所施無機(jī)肥種類:氮肥為尿素(含氮 46%),磷肥為過磷酸鈣(含P2O516%),鉀肥為硫酸鉀(含K2O 50%)。復(fù)合肥含有機(jī)質(zhì)20%,含氮 8%、磷 4%、鉀 3%。試驗設(shè)置6個處理:不施肥對照(CK),常規(guī)氮磷鉀肥(CF),30%豬糞+70%常規(guī)氮磷鉀肥(OF),全量秸稈還田+秸稈腐熟劑+常規(guī)氮磷鉀肥(FS),豬糞堆肥有機(jī)無機(jī)復(fù)合肥1 800 kg/hm2(OI1),豬糞堆肥有機(jī)無機(jī)復(fù)合肥3 600 kg/hm2(OI2)(詳見表1)。每個處理設(shè)置3次重復(fù),共設(shè)18個小區(qū),每個小區(qū)面積為50 m2(5 m×10 m)。小區(qū)間隔1 m,試驗區(qū)外設(shè)置5 m保護(hù)行。不同施肥處理的小區(qū)隨機(jī)排列。

    表1 各小區(qū)施肥量(kg/小區(qū))

    1.1.2 樣品采集 土壤樣品分別于2012年和2013年小麥、玉米收獲后進(jìn)行取樣,取0~20 cm土層土壤樣品,分別于每個小區(qū)取9個點的土樣,混勻后剔除石塊和植物殘根等雜物,于-70 ℃保藏,用于Biolog功能的檢測。土壤微生物種群功能多樣性,采用美國Matrix Technologies Corporation生產(chǎn)的Biolog-Eco板測定。

    1.1.3 培養(yǎng)基 采用高氏 1 號瓊脂培養(yǎng)基(倒平板前按 1 mL 培養(yǎng)基含 K2Cr2O775 μg的量加入無菌 K2Cr2O7溶液)分離放線菌,牛肉膏蛋白胨培養(yǎng)基分離細(xì)菌,馬丁氏培養(yǎng)基分離真菌[17-18]。

    1.2 方法

    1.2.1 Biolog試驗 ①制備土壤菌懸液:分別稱取土壤樣品20 g,置于裝有100 mL生理鹽水的無菌錐形瓶中,200 r/min震蕩40 min,稀釋100倍,制成初始土壤菌懸液[19],備用。②制備土壤接種液:將初始土壤菌懸液混勻并靜置5 min,使土壤中顆粒較大的砂礫沉淀。校正濁度計后,把上述菌懸液逐滴向20 mL無菌生理鹽水中緩慢滴加,調(diào)節(jié)土壤接種液濁度為63%。③接種:將土壤接種液倒入接種液槽,用8通道移液器吸取接種液,每個微孔注入150 μL,將微板編號,25 ℃培養(yǎng)[20-21]。④讀數(shù):每隔24 h讀取各板孔在590 nm和750 nm處的吸光值(A),并保存結(jié)果文件(讀取平臺為美國Biolog公司的Biolog Reader)。

    1.2.2 Biolog數(shù)據(jù)分析 利用Microlog4.2軟件,將各數(shù)據(jù)轉(zhuǎn)換成Excel格式,以便于進(jìn)一步計算分析。①采用各Eco-MicroPlate吸光值數(shù)據(jù),計算平均顏色變化率(average well color development,AWCD)。采用Classen等[22]的方法計算Biolog微孔板的顏色變化平均值,用每個對應(yīng)孔590 nm的吸光值與750 nm的吸光值的差減去對照孔吸光值,結(jié)果即為反應(yīng)單孔實際顏色變化的吸光值,其中數(shù)值小于0.06時按0處理。利用處理后得到的數(shù)據(jù)按照Garland等[23]的方法計算AWCD,公式如下:

    (1)

    式中:Ci為每個微板孔的吸光值,C0為對照孔吸光值,n為培養(yǎng)基種類數(shù)(ECO-MicroPlate為31種培養(yǎng)基)。

    ②土壤懸液微生物物種豐富度用Shannon-Weinner指數(shù)(H)表征,各供試土樣微板數(shù)據(jù)的Shannon-Weinner指數(shù)(H)計算公式如下:

    H=-∑(Pi×log Pi);

    Pi=(Ci-C0)/∑(Ci-C0)

    (2)

    式中:Ci為每個微板孔的吸光值,C0為對照孔吸光值,Pi為含底物的微孔和對照孔吸光值的差值與整個微板總體差值的比值。

    ③土壤懸液微生物物種均勻度用Shannon-Weinner指數(shù)(E)表征,各供試土樣微板數(shù)據(jù)的Shannon-Weinner均勻度指數(shù)(E)計算公式如下:

    E=H/lnS

    (3)

    式中:H為Shannon-Weinner物種豐富度指數(shù),S為有顏色變化的微孔數(shù)目。

    由于Biolog試驗數(shù)據(jù)量較大而且較為繁瑣,因此很難直接比較各個處理間的差異,所以通常借助于主成分分析來比較不同樣品微生物群落功能多樣性。本研究取培養(yǎng)至72 h時的AWCD值,經(jīng)過計算后,得到單孔相對光密度值進(jìn)行主成分分析(PCA)[24-28]。

    2 結(jié)果與分析

    2.1 不同施肥處理的土壤微生物平均顏色變化率(AWCD)

    2012年和2013年小麥季和玉米季不同施肥處理土壤的AWCD分析結(jié)果見圖1。由圖1可知,小麥季和玉米季不同施肥處理土壤微生物AWCD的總體變化趨勢隨著時間的延長而逐漸提高,在培養(yǎng)24~72 h期間,曲線斜率最大,AWCD增長率最大,表明此階段土壤微生物碳源代謝活性最高;72 h后AWCD同樣繼續(xù)升高,但增長率逐漸下降,最后進(jìn)入平穩(wěn)期,說明微生物代謝活動趨于平緩。不同處理間相比,施肥處理AWCD值高于不施肥對照(CK);施肥處理中FS在整個培養(yǎng)階段中高于其他處理,OF和CF次之,OI1和OI2略低。由此可以說明,不施肥對照土壤微生物的代謝活性低于其他施肥方式,秸稈還田和施用有機(jī)肥后土壤微生物代謝活性高于單施化肥,其中使用FS處理更能提高土壤微生物代謝活性。

    2.2 不同施肥處理土壤微生物種群多樣性指數(shù)的變化

    為了研究采取不同施肥措施后土壤微生物種群的變化,計算了小麥季和玉米季Shannon-Weinner豐富度指數(shù)(H)和Shannon-Weinner均勻度指數(shù)(E),二者可以較好地反映土壤中微生物種群多樣性[29],結(jié)果見表2。

    圖1 2012、2013年小麥季和玉米季不同施肥處理土壤微生物的孔平均顏色變化率Fig.1 Different fertilization AWCD of wheat season and corn season in 2012 and 2013A:2012年小麥季;B:2012年玉米季;C:2013年小麥季;D:2013年玉米季,下圖同A:2012 wheat season;B:2012 corn season; C:2013 wheat season;D:2013 corn season,the following figure with

    2012年2013年處理物種豐富度指數(shù)(H)物種均勻度指數(shù)(E)物種豐富度指數(shù)(H)物種均勻度指數(shù)(E)小麥季玉米季小麥季玉米季小麥季玉米季小麥季玉米季CK1.37±0.01a1.38±0.02a0.42±0.01a0.43±0.01a1.32±0.01a1.24±0.01a0.41±0.00a0.38±0.00aCF1.39±0.01b1.39±0.05a0.42±0.01a0.43±0.01a1.36±0.06a1.24±0.01a0.42±0.02ab0.38±0.01aOF1.44±0.01c1.41±0.01ab0.45±0.01b0.44±0.00ab1.36±0.03ab1.26±0.04a0.42±0.01ab0.39±0.00aFS1.46±0.01d1.44±0.01b0.45±0.02b0.45±0.00b1.39±0.01b1.27±0.02a0.43±0.00b0.40±0.01aOI11.39±0.01b1.38±0.01a0.43±0.01ab0.44±0.01ab1.33±0.02a1.26±0.04a0.41±0.01a0.39±0.01aOI21.39±0.01b1.38±0.01a0.43±0.01ab0.44±0.01ab1.33±0.01a1.25±0.01a0.41±0.00a0.39±0.00a

    注:表中AWCD為培養(yǎng)至72 h的結(jié)果;同一列中具有不同字母表示結(jié)果差異顯著(P< 0.05)

    從表2可以看出,總體趨勢上,小麥季和玉米季不同施肥處理的微生物群落物種豐富度指數(shù)(H)和均勻度指數(shù)(E)中CK處理的最低,CF處理略高于CK處理,OI1、OI2處理均高于CK、CF, FS處理和OF處理均明顯高于其他處理,F(xiàn)S處理豐富度和均勻度最高。結(jié)果說明:不施肥處理跟常規(guī)氮磷鉀肥處理(CF)微生物群落功能多樣性較低,增施有機(jī)肥后會增強(qiáng)土壤微生物群落功能多樣性,其中有機(jī)肥替代部分化肥(OF)和秸稈還田(FS)效果明顯。

    2.3 土壤微生物碳源利用的主成分分析(PCA)

    2.3.1 不同碳源在主成分上的荷載值 按照Eco板上31種不同碳源的化學(xué)基團(tuán)性質(zhì)的不同,將31種碳源分成糖類、氨基酸類、羧酸類、胺類、聚合物類和其他。通過不同碳源在第一主成分(PC1)、第二主成分(PC2)中的荷載因子大小的比較,可以得出分別與第一、第二主成分具有較高相關(guān)性的碳源種類與數(shù)量(見表3)。

    表3 主成分分析中高相關(guān)性碳源分布情況

    由表3可知,在冬小麥-夏玉米輪作體系中,土壤微生物都不同程度地利用糖類、氨基酸類、羧酸類、胺類、聚合物類等碳源。其中,胺類和其他碳源利用程度最低,而糖類、羧酸類、聚合物類、氨基酸類是土壤微生物碳源的主要來源。

    2.3.2 不同施肥處理的碳源利用主成分分析 為了防止出現(xiàn)Eco板中某碳源底物提前消耗完畢的可能,本研究采用72 h時測定的AWCD數(shù)據(jù),應(yīng)用SPSS軟件進(jìn)行標(biāo)準(zhǔn)化處理,對微生物種群的代謝特征進(jìn)行了主成分分析(Principal Component Analysis,PCA)。各季中,除了第一主成分和第二主成分,后面的主成分貢獻(xiàn)率均較小,所以選取上述數(shù)據(jù)中前兩個主成分進(jìn)行分析。以第一主成分為橫軸,第二主成分為縱軸作圖,得到不同施肥處理土壤微生物對31種碳源利用情況的主成分分析圖(圖2)。

    圖2 不同施肥處理土壤微生物碳源利用的主成分分析Fig.2 PCA for soil microorganism carbon source utilization

    圖2中4個主成分圖上,6組處理根據(jù)離散程度大致分成3組。FS和OF處理位于PC1軸和PC2軸正端,而CK和CF軸位于PC1軸和PC2軸負(fù)端,OI1和OI2則位于PC1軸負(fù)端和PC2正端,各組內(nèi)2處理分別存在一定的相似性。除2012年小麥季外,CK、CF處理與另3個季度其他處理相比差異都顯著,結(jié)果表明,不同施肥處理在PC軸上出現(xiàn)了明顯的差異,說明不同處理對土壤微生物的碳源利用類型有顯著影響,使用不同肥料其微生物群落功能的多樣性就有差異,并且差異顯著。按小麥季和玉米季時間順序來看,3組處理之間的差異越來越明顯,說明隨著施肥時間的延長,其對土壤微生物碳源代謝的影響越來越明顯。

    2.4 土壤微生物功能多樣性與土壤化學(xué)性質(zhì)的相關(guān)性分析

    為探究土壤化學(xué)性質(zhì)對微生物功能多樣性的影響,測定了土壤pH、有機(jī)質(zhì)(Soil organic matter,SOM)、全氮(Total N,TN)、速效磷(Available P,AP)和速效鉀(Available K,AK)5種土壤化學(xué)性質(zhì),并以此進(jìn)行RDA分析(圖3)。

    如圖3所示,不同處理土壤微生物功能多樣性與土壤化學(xué)性質(zhì)相關(guān)性明顯。各個處理從CK處理開始向SOM、TN、AP、AK箭頭方向呈現(xiàn)一定梯度分布。CK處理與土壤pH(P<0.01)呈現(xiàn)正相關(guān),其余處理與SOM、TN、AP、AK(P<0.05)呈正相關(guān),且相關(guān)性明顯。其中OF、FS處理與土壤有機(jī)質(zhì)、全氮速效磷、速效鉀的相關(guān)性最大; OI1、OI2處理次之;CF處理在小麥季與在玉米季略有不同,小麥季受理化影響明顯,與OI2處理受影響水平相近,玉米季受理化性質(zhì)影響稍小,與OI1處理受影響水平相近。小麥季中OF、FS處理微生物功能多樣性與土壤有機(jī)質(zhì)相關(guān)性明顯,玉米季則不明顯。

    圖3 不同施肥處理土壤微生物功能多樣性與土壤化學(xué)性質(zhì)的RDA分析Fig.3 RDA between functional diversity of soil microbial communities and soil characteristics of soil samples

    3 討 論

    AWCD能夠準(zhǔn)確反映土壤微生物代謝活性的強(qiáng)弱,微生物利用碳源能力的強(qiáng)弱可以根據(jù)AWCD值的高低來判斷。而在Eco板中共有31種碳源,因此本研究中AWCD能夠反映不同施肥處理下土壤微生物利用碳源能力的大小[31-32]。本研究結(jié)果表明,各施肥處理AWCD值均顯著高于對照,其中秸稈還田處理和70%NPK+有機(jī)肥處理最高,表明采取增施有機(jī)肥或秸稈還田措施能顯著提高農(nóng)田土壤中微生物的碳源利用能力,而其他3個施肥處理碳源的利用能力略低。這可能是由于增施有機(jī)肥或秸稈還田會大量增加土壤中有機(jī)質(zhì)的含量,大大改善了土壤性狀,從而改善微生物生境,使其數(shù)量和種類增加,土壤中原來不被利用或者利用少的碳源開始利用并增強(qiáng),最終結(jié)果是提高了土壤微生物群落的數(shù)量和代謝功能的多樣性。朱海平等[34]、王光華等[35]研究也表明,施用廄肥、秸桿能顯著增強(qiáng)土壤微生物生物量和基礎(chǔ)呼吸量,提高維持土壤微生物多樣性及功能。相反,單施化肥處理的土壤微生物對碳源的利用能力也高于不施肥處理,這與張瑞等[12]、李東坡等[35]研究相似。復(fù)合肥處理較不施肥處理的土壤微生物碳源利用水平高,但相較于FS和OF處理低,可能是由于復(fù)合肥中養(yǎng)分比例比較固定,難以滿足土壤和各種作物的需要。土壤施用有機(jī)肥可使微生物能夠利用的碳源趨于多元化,微生物能夠利用的碳源種類越多,其生存能力就越強(qiáng),代謝活性也就越高。

    本研究各處理物種豐富度(H)和均勻度(E)結(jié)果表明,F(xiàn)S、OF處理土壤微生物種群多樣性最高,OI2、OI1處理次之,CF處理較低,CK處理最低。這說明秸稈還田或增施有機(jī)肥會顯著增強(qiáng)土壤微生物種群多樣性,有機(jī)無機(jī)復(fù)合肥對提高土壤微生物種群多樣性效果不明顯。雖然兩年數(shù)據(jù)略有差別,但處理間物種豐富度(H)和均勻度(E)變化趨勢一致,且2013年各處理間差異更為明顯。分析原因,OF處理由于用有機(jī)肥替代部分化肥,提高了土壤中有機(jī)質(zhì)含量,為微生物提供了良好的生存環(huán)境,因此其微生物豐富度指數(shù)和均勻度指數(shù)遠(yuǎn)高于CK和CF。小麥季和玉米季土壤微生物物種豐富度指數(shù)和均勻度指數(shù)FS處理最高,可能是由于秸稈還田極大地增加了土壤中有機(jī)物的含量,大量纖維素在土壤中起到了疏松土質(zhì)的作用,增加了土壤含氧量,提高了土壤微生物的代謝強(qiáng)度;此外,施加的秸稈腐熟劑,能夠加快秸稈的腐熟,使更多的營養(yǎng)物質(zhì)進(jìn)入土壤,增加了微生物的豐富度。時鵬等[36]通過玉米連作及施肥試驗也得出了同樣的結(jié)果。

    主成分分析結(jié)果顯示,對于PC1和PC2起分異作用的主要碳源是糖類、氨基酸類和羧酸類物質(zhì),氨基酸類和糖類在PC1軸上的權(quán)重較大,而與PC2相關(guān)性較大的碳源主要是羧酸類、糖類和聚合物類。糖類和氨基酸類碳源是區(qū)分各個處理的主要碳源。此外,CK處理和CF處理,OF處理和FS處理,OI1處理和OI2處理碳源利用方式分別相似,并且與2013年及2012年相比其相似性更明顯。可能是由于隨著施肥處理時間的增長,各個處理土壤微生物群落趨于穩(wěn)定,且施肥種類相似的處理其土壤微生物群落結(jié)構(gòu)和功能也相似。

    土壤微生物功能多樣性與土壤化學(xué)性質(zhì)的相關(guān)性分析結(jié)果表明,土壤化學(xué)性質(zhì)與土壤微生物功能多樣性密切相關(guān)。施肥能改變土壤中氮磷鉀及有機(jī)質(zhì)含量,從而影響土壤微生物功能的多樣性。全量秸稈還田+秸稈腐熟劑+常規(guī)氮磷鉀肥(FS)、有機(jī)肥+常規(guī)氮磷鉀肥(OF)處理土壤中有機(jī)質(zhì)、全氮速效磷、速效鉀含量最高,土壤微生物功能多樣性程度也最高;豬糞堆肥有機(jī)無機(jī)復(fù)合肥(OI2)、豬糞堆肥有機(jī)無機(jī)復(fù)合肥(OI1)次之,全量氮磷鉀肥(CF)較低,不施肥對照(CK)最低。這表明秸稈還田+腐熟劑和有機(jī)替代部分化肥能提高土壤養(yǎng)分含量,同時大幅提高土壤微生物功能多樣性,Zhao等[37]、Li等[38]也得出相似結(jié)論;豬糞堆肥有機(jī)無機(jī)復(fù)合肥也能提高土壤養(yǎng)分含量和微生物功能多樣性,效果比常規(guī)氮磷鉀肥(CF)稍好,這與Wei等[39]結(jié)論一致。RDA分析中,CF處理玉米季與小麥季有所差異,其原因可能是CF與其他處理相比,其玉米季土壤中氮磷鉀流失程度高于小麥季[40]。

    綜上所述,所有施肥處理均有利于維持土壤微生物的群落多樣性。秸稈還田+腐熟劑+常規(guī)施肥及有機(jī)肥替代部分氮磷鉀肥顯著提高了土壤微生物物種豐富度指數(shù)和均勻度指數(shù),增強(qiáng)了微生物代謝活性,提高了土壤微生物功能多樣性。豬糞堆肥有機(jī)無機(jī)復(fù)合肥也能提高土壤微生物種群和功能多樣性,但效果低于秸稈還田和有機(jī)替代,僅比常規(guī)氮磷鉀施肥好。所以,施肥可顯著影響土壤微生物功能多樣性,但其影響程度可能與施肥種類、施肥方式、土壤性質(zhì)等因素密切相關(guān)。秸稈還田+腐熟劑和有機(jī)替代是否已是最優(yōu)施肥方式,改變腐熟劑種類或者有機(jī)替代比例能否達(dá)到更好效果,這還有待于進(jìn)一步研究。

    [1] Pankhurst CE, Keller KO, Gupta V V S R. Biodiversity of soil microbial communities in agricultural systems[J].Biodiversity and Conservation,1996,5(2):197-209.

    [2] Schloter M, Dilly O, Munch JN. Indicators for evaluating soil quality[J].Agriculture Ecosystems and Environment,2003,98(1/3): 255-262.

    [3] Bardgett R D, Speir T W, Ross D J,et al. Impact of pasture contamination by copper, chromium and arsenic timber preservative on soil microbial properties and nematodes[J].Biology and Fertility of Soils,1994,18: 71-79.

    [4] Enwall K, Nyberg K, Bertilsson S,et al. Long-term impact of fertilization on activity and composition of bacterial communities and metabolic guilds in agricultural soil[J].Soil Biology and Biochemistry,2007,39(1): 106-115.

    [5] Bradley J C, Diane S S, Duffy E,et al. Effects of biodiversity on the functioning of trophic groups and ecosystems[J].Nature,2006,443: 989-992.

    [6] Sparling G P. Ratio of microbial biomass carbon to soil organic carbon as a sensitive indicator of changes in soil organic matter[J].Australian Journal of Soil Research,1992,30: 195-207.

    [7] Buenemann E K, Bossio D A, Smithson P C, et al. Microbial community composition and substrate use in a highly weathered soil as affected by crop rotation and P fertilization[J].Soil Biology and Biochemistry,2004,36(6): 889-901.

    [8] Hai L, Li X-G, Li F-M, et al. Long-term fertilization and manuring effects on physically-separated soil organic matter pools under a wheat-wheat-maize cropping system in an arid region of China[J].Soil Biology and Biochemistry,2010,42(2): 253-259.

    [9] Fauci MF, Dick R P. Soil microbial dynamics: short-and long-term effects of organic and inorganic nitrogen[J].Soil Science Society of America Journal,1994,58: 801-808.

    [10]王銳,林先貴,陳瑞蕊,等.長期不同施肥對潮土芽胞桿菌數(shù)量的影響及其優(yōu)勢度的季節(jié)變化[J].土壤學(xué)報,2013,50(4): 778-785.

    [11]王才斌,鄭亞萍,梁曉艷,等.施肥對旱地花生主要土壤肥力指標(biāo)及產(chǎn)量的影響[J].生態(tài)學(xué)報,2013,33(4): 1300-1307.

    [12]張瑞,張貴龍,陳冬青,等.不同施肥對農(nóng)田土壤微生物功能多樣性的影響[J].中國農(nóng)學(xué)通報,2013,29(2): 133-139.

    [13]劉振香,劉鵬,賈緒存,等.不同水肥處理對夏玉米田土壤微生物特性的影響[J].應(yīng)用生態(tài)學(xué)報,2015,26(1): 113-121.

    [14]譚周進(jìn),周衛(wèi)軍,張楊珠,等.不同施肥制度對稻田土壤微生物的影響研究[J].植物營養(yǎng)與肥料學(xué)報,2007,13(3): 430-435.

    [15]劉曉蒨,涂世華,孫錫發(fā),等.長期不同施肥對稻田土壤微生物群落功能多樣性的影響[J].生態(tài)學(xué)報,2013,33(17): 5210-5217.

    [16]Roger BF, Tate RL. Temporal analysis of the soil microbial community along a top sequence in Pineland soils[J].Soil Biology and Biochemistry,2001,33(10): 1389-1401.

    [17]李阜棣,喻子牛,何紹江.農(nóng)業(yè)微生物學(xué)實驗技術(shù)[M].北京:中國農(nóng)業(yè)出版社,1996:69-72.

    [18]趙斌,何紹江.微生物學(xué)實驗[M].北京:科學(xué)出版社,2002:251,253-254.

    [19]鄭華, 歐陽志云, 方治國,等.BIOLOG在土壤微生物群落功能多樣性研究中的應(yīng)用[J].土壤學(xué)報,2004,(3):456-461.

    [20]Kennedy AC, Smith KL. Soil microbial diversity and the sustainability of agricultural soils[J].Plant and Soil, 1995, 170(1): 75-86.

    [21]侯曉杰,汪景寬,李世鵬,等.不同施肥處理與地膜覆蓋對土壤微生物群落功能多樣性的影響[J].生態(tài)學(xué)報,2007,27(2): 655-661.

    [22]Classen AT, Boyle SI, Heskins KE, et al. Community-level physiological profiles of bacteria and fungi: plate type and incubation temperature influences on contrasting soils[J].FEMS Microbiology Ecology,2003,44: 319-328.

    [23]Gafand JL, Mills AL. Classification and characterization of heterotrophic microbial communities on the basis of patterns of community-level sole-carbon-source utilization[J].Applied and Environmental Microbiology,1991,57: 2351-2359.

    [24]郭寶林,楊俊霞,李永慈,等.主成分分析法在仁用杏品種主要經(jīng)濟(jì)性狀選種上的應(yīng)用研究[J].林業(yè)科學(xué),2000,36(6): 53-56.

    [25]何曉群.多元統(tǒng)計分析[M]. 北京:中國人民大學(xué)出版社,2008.

    [26]Kela PW, Jason AG, Matthias G, et al. Data transformations in the analysis of community-level substrate utilization data from microplates[J].Journal of Microbiological Methods,2007,69: 461-469.

    [27]Grove JA, Kautola H, Javadpour S, et al. Assessment of changes in the microorganism community in a biofilter[J].Biochemical Engineering Journal,2004,18: 111-114.

    [28]Min L, Xiao M-X. Effect of heavy metals on substrate utilization pattern, biomass, and activity of microbial communities in a reclaimed mining wasteland of red soil area[J].Ecotoxicology and Environmental Safety, 2007,66: 217-223.

    [29]Magurran AE. Ecological diversity and its measurement[M].USA: Princeton University Press,1988,141-162

    [30]郝黎仁,樊元,郝哲歐,等.SPSS實用統(tǒng)計分析[M].北京:中國水利水電出版社,2003.

    [31]O’Donnell AG, Seasman M, Macrae A, et al. Plants and fertilisersas drivers of change in microbial community structure and function in soils[J].Plant and Soil, 2001,232(1): 135-145.

    [32]Benizri E, Amiaud B. Relation between plants and soil microbial communities in fertilized grasslands[J].Soil Biochemistry, 2005,37(11): 2055-2064.

    [33]朱海平,姚槐應(yīng),張勇勇,等.不同培肥管理措施對土壤微生物生態(tài)特征的影響[J].土壤通報,2003, 34(2): 140-142.

    [34]王光華,劉俊杰,齊曉寧,等.Biolog和PCR-DGGE技術(shù)解析施肥對德惠黑土細(xì)菌群落結(jié)構(gòu)和功能的影響[J].生態(tài)學(xué)報, 2008, 28(1): 220-226.

    [35]李東坡,武志杰,陳利軍,等.長期培肥黑土微生物量磷動態(tài)變化及影響因素[J].應(yīng)用生態(tài)學(xué)報,2004, 15(10): 1897-1902.

    [36]時鵬,高強(qiáng),王淑平,等.玉米連作及其施肥對土壤微生物群落功能多樣性的影響[J].生態(tài)學(xué)報,2010, 30(22): 6173-6182.

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

    [38]Juan Li, Julia Mary Cooper, Zhi’an Lin,et al. Soil microbial community structure and function are significantly affected by long-term organic and mineral fertilization regimes in the North China Plain[J].Applied Soil Ecology,2015,(96): 75-87.

    [39]Wei Tian, Lei Wang, Yan Li,et al. Responses of microbial activity, abundance, and community in wheat soil after three years of heavy fertilization with manure-based compost and inorganic nitrogen[J]. Agriculture, Ecosystems and Environment,2015,(213): 219-227.

    [40]Liyue Guo, Guanglei Wu, Yong Li,et al. Effects of cattle manure compost combined with chemical fertilizer on topsoil organic matter, bulk density and earthworm activity in a wheat-maize rotation system in Eastern China[J].Soil & Tillage Research,2016,(156): 140-147.

    Effects of Different Fertilization on Soil Microbial Community and Functional Diversity in Maize-Wheat Crop Rotation

    XING Peng-fei1, 2,3, WU Xiao-sen2, 3,4, GAO Sheng-chao2,3, LI Hong-jie5,ZHAO Tong-kai5,ZHOU Xiao-lin5, SHEN De-long2,3, SUN Jun-de1

    (1.Coll.ofLand&Environ’t,ShenyangAgric.Uni.,Shenyang110866; 2.Inst.ofAgric.Res. &Reg’lPlan.,ChineseAcad.ofAgric.Sci.,Beijing100081; 3.LaboratoryofQualityandSafetyriskAssessmentforMicrobialProductsMinistryofAgriculture,Beijing100081;4.Coll.ofLifeSci.,Agric.Uni.ofHebei,Baoding071000;5.DezhouAcad.ofAgric.Sci.,Dezhou253015 )

    The diversity of soil microbes can reflect the soil fertility and different fertilization also works importantly on soil microbial population and diversity. In this study, soil after the harvest of two-year in succession of maize season and wheat season in Dezhou field ecosystem station run by CAAS was taken as the research object to study the influence of different fertilization on microbial community functional diversity with six different fertilization treatments using Biolog technique. The results showed that each of six fertilization treatments had significant differences in AWCD (average well color development), and each of them were all significantly higher than that of the control without fertilization (CK). Soil microbial communities richness index (H) and evenness index (E) of fertilization treatments FS (conventional NPK+ straw + straw decomposing inoculants) and OF (70% conventional NPK + 30% swine excrement) were higher than others’. PCA (principal component analysis) and RDA (redundancy analysis) showed that OF and FS had similar functional diversity of microbial communities, and they had high level relativity with SOM (soil organic matter), TN (total N), AP (available P), and AK (available K). The treatment of swine excrement compost + organic + inorganic compound fertilizer 1 800 kg/hm2(OI1) and the treatment of swine excrement compost + organic + inorganic compound fertilizer 3 600 kg/hm2(OI2) were treated similarly, their functional diversity were lower than FS and OF, and higher than CF. To sum up the above mentioned, the degree of soil microbial community and functional diversity of OF treatment and FS treatment was higher than those other treatments, indicated that straw field-returning plus straw decomposing preparation and organic fertilizer to replace part of NPK fertilizer could improve soil microbial community and functional diversity, and conduce to protect soil microbial ecology.

    different fertilization; soil microbe; community and functional diversity; Biolog

    農(nóng)業(yè)部公益性農(nóng)業(yè)行業(yè)科研專項(201103004);農(nóng)業(yè)部生物有機(jī)肥創(chuàng)制重點實驗室開放課題資助

    邢鵬飛 男,碩士研究生。主要從事微生物生態(tài)學(xué)方面的研究。E-mail:xingruofeng@126.com

    2015-09-07;

    2015-12-24

    Q939.96

    A

    1005-7021(2016)01-0022-08

    10.3969/j.issn.1005-7021.2016.01.005

    * 通訊作者。孫軍德 男,教授,博士生導(dǎo)師。主要從事微生物農(nóng)藥與肥料、食藥用真菌菌種選育與栽培等研究。

    E-mail: sunjunde108@163.com

    沈德龍 男,研究員,碩士生導(dǎo)師。主要從事農(nóng)用微生物資源和應(yīng)用基礎(chǔ)研究、微生物肥料行業(yè)管理工作。E-mail:shendelong@caas.cn

    猜你喜歡
    均勻度氮磷碳源
    低播量下雜交稻產(chǎn)量形成對種植均勻度的響應(yīng)
    作物研究(2023年2期)2023-05-28 13:44:14
    玻利維亞擬建新的氮磷鉀肥料工廠
    緩釋碳源促進(jìn)生物反硝化脫氮技術(shù)研究進(jìn)展
    不同碳源對銅溜槽用鋁碳質(zhì)涂抹料性能的影響
    昆鋼科技(2021年6期)2021-03-09 06:10:20
    均勻度控制不佳可致肉種雞晚產(chǎn)
    四甘醇作碳源合成Li3V2(PO4)3正極材料及其電化學(xué)性能
    錦綸長絲染色均勻度判色新方法
    植物需要氮磷鉀
    復(fù)方丹參片中冰片的含量均勻度研究
    中成藥(2014年10期)2014-02-28 22:29:24
    不同氮磷水平下AM真菌對棉花的侵染率及氮磷含量的影響
    两个人的视频大全免费| 免费高清在线观看视频在线观看| av免费在线看不卡| 热re99久久国产66热| 丰满饥渴人妻一区二区三| 91成人精品电影| 大香蕉久久网| 日韩伦理黄色片| 亚洲国产精品一区三区| 亚洲欧美中文字幕日韩二区| 午夜福利视频在线观看免费| 在线天堂最新版资源| 婷婷色综合大香蕉| 免费观看性生交大片5| 一边摸一边做爽爽视频免费| 日韩在线高清观看一区二区三区| 亚洲精品乱码久久久久久按摩| 精品久久蜜臀av无| 午夜福利影视在线免费观看| av在线老鸭窝| 成年美女黄网站色视频大全免费 | 老女人水多毛片| 国产精品.久久久| 制服诱惑二区| 免费观看性生交大片5| 热re99久久精品国产66热6| av视频免费观看在线观看| 国产免费福利视频在线观看| 午夜激情久久久久久久| 黄色视频在线播放观看不卡| 考比视频在线观看| 18在线观看网站| 亚洲内射少妇av| 日本欧美国产在线视频| 免费人成在线观看视频色| 免费观看在线日韩| 午夜久久久在线观看| 久久国内精品自在自线图片| 精品一品国产午夜福利视频| 亚洲激情五月婷婷啪啪| 在线观看国产h片| 少妇人妻精品综合一区二区| 日韩成人伦理影院| 日本av免费视频播放| 久久影院123| 少妇 在线观看| 亚洲丝袜综合中文字幕| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 国产精品99久久99久久久不卡 | 青春草亚洲视频在线观看| 亚洲欧美一区二区三区国产| 人妻 亚洲 视频| 男的添女的下面高潮视频| 一级片'在线观看视频| 99久国产av精品国产电影| 一级爰片在线观看| a级毛片在线看网站| 夫妻午夜视频| 丰满饥渴人妻一区二区三| 尾随美女入室| 最黄视频免费看| 天美传媒精品一区二区| 亚洲av在线观看美女高潮| 国产免费现黄频在线看| 亚洲av成人精品一二三区| 日韩成人av中文字幕在线观看| 国产69精品久久久久777片| 久久av网站| 国产精品偷伦视频观看了| 成人漫画全彩无遮挡| 亚洲av成人精品一二三区| 精品久久久久久电影网| 亚洲综合精品二区| 欧美丝袜亚洲另类| 欧美老熟妇乱子伦牲交| 最新的欧美精品一区二区| 国产亚洲欧美精品永久| 国产黄色视频一区二区在线观看| 日日摸夜夜添夜夜添av毛片| 成人黄色视频免费在线看| 精品少妇久久久久久888优播| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 少妇丰满av| 国产成人免费无遮挡视频| 自线自在国产av| 久久久国产精品麻豆| 亚洲欧美成人综合另类久久久| 国产av国产精品国产| 在线亚洲精品国产二区图片欧美 | 成年人午夜在线观看视频| 国产色婷婷99| 精品人妻偷拍中文字幕| 天堂8中文在线网| 久热久热在线精品观看| 人人妻人人澡人人爽人人夜夜| 亚洲美女视频黄频| 成年人午夜在线观看视频| 一级毛片aaaaaa免费看小| a级毛片在线看网站| 国产片内射在线| 国产片特级美女逼逼视频| 中文字幕久久专区| 蜜臀久久99精品久久宅男| 91久久精品国产一区二区三区| 久久精品夜色国产| 国产精品99久久久久久久久| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 国产日韩欧美亚洲二区| 天天躁夜夜躁狠狠久久av| 久久久久人妻精品一区果冻| 蜜桃久久精品国产亚洲av| 免费高清在线观看日韩| 欧美另类一区| 18禁动态无遮挡网站| 国产成人精品在线电影| 欧美精品国产亚洲| 妹子高潮喷水视频| 日韩熟女老妇一区二区性免费视频| 久久久久久伊人网av| 热99国产精品久久久久久7| 久久精品熟女亚洲av麻豆精品| 人妻夜夜爽99麻豆av| 亚洲欧洲日产国产| 日韩一区二区三区影片| 亚洲人成77777在线视频| 国产精品国产三级国产av玫瑰| 一本一本综合久久| 伊人久久精品亚洲午夜| 国产一区亚洲一区在线观看| 久久99热6这里只有精品| 亚洲在久久综合| 国产国语露脸激情在线看| 日韩欧美精品免费久久| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| a 毛片基地| 久久av网站| 最新中文字幕久久久久| 午夜福利影视在线免费观看| 欧美变态另类bdsm刘玥| 久久久久精品久久久久真实原创| 日韩亚洲欧美综合| av视频免费观看在线观看| 制服诱惑二区| 日韩亚洲欧美综合| 中国美白少妇内射xxxbb| 国产国语露脸激情在线看| 免费观看的影片在线观看| 在线观看免费视频网站a站| 免费高清在线观看视频在线观看| 最近手机中文字幕大全| 日韩成人伦理影院| 各种免费的搞黄视频| 老司机影院成人| 另类亚洲欧美激情| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说| 在线播放无遮挡| 欧美人与性动交α欧美精品济南到 | 婷婷色综合www| av专区在线播放| 国产 一区精品| 好男人视频免费观看在线| 在线观看一区二区三区激情| 国产 一区精品| 人妻人人澡人人爽人人| 国产伦理片在线播放av一区| 亚洲国产av影院在线观看| 中文字幕精品免费在线观看视频 | 一级,二级,三级黄色视频| 国产高清国产精品国产三级| videos熟女内射| 高清黄色对白视频在线免费看| 成人影院久久| 丁香六月天网| 亚洲精品,欧美精品| 18禁在线无遮挡免费观看视频| tube8黄色片| a级毛片在线看网站| 亚洲,一卡二卡三卡| 一级爰片在线观看| 3wmmmm亚洲av在线观看| 亚洲欧美日韩另类电影网站| videosex国产| 搡老乐熟女国产| 乱人伦中国视频| a级毛片在线看网站| 亚洲精品av麻豆狂野| 少妇精品久久久久久久| 国产精品偷伦视频观看了| videossex国产| 午夜视频国产福利| 各种免费的搞黄视频| 久久99热6这里只有精品| 成人亚洲欧美一区二区av| 欧美国产精品一级二级三级| 精品人妻一区二区三区麻豆| freevideosex欧美| 日韩欧美精品免费久久| 欧美人与善性xxx| 国产成人免费观看mmmm| 麻豆精品久久久久久蜜桃| 老司机影院毛片| 日韩av免费高清视频| 91精品国产九色| 肉色欧美久久久久久久蜜桃| av在线播放精品| 在线观看一区二区三区激情| 一区二区av电影网| 秋霞在线观看毛片| 人妻夜夜爽99麻豆av| 天美传媒精品一区二区| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 少妇人妻久久综合中文| 婷婷成人精品国产| 亚洲精品亚洲一区二区| 久久午夜福利片| 考比视频在线观看| 777米奇影视久久| 婷婷成人精品国产| freevideosex欧美| 丰满少妇做爰视频| 欧美日本中文国产一区发布| 多毛熟女@视频| 如何舔出高潮| 亚洲av中文av极速乱| 精品少妇内射三级| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 内地一区二区视频在线| 欧美精品人与动牲交sv欧美| 国语对白做爰xxxⅹ性视频网站| 久久av网站| 中国三级夫妇交换| 在线观看www视频免费| videosex国产| 日韩强制内射视频| 一本久久精品| 亚洲伊人久久精品综合| 国产黄频视频在线观看| 亚洲伊人久久精品综合| 日韩强制内射视频| 国内精品宾馆在线| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 最新中文字幕久久久久| 日本爱情动作片www.在线观看| 少妇 在线观看| 亚洲精华国产精华液的使用体验| 91国产中文字幕| 日日爽夜夜爽网站| 亚洲少妇的诱惑av| 久久99一区二区三区| 欧美精品国产亚洲| 人妻系列 视频| 欧美日本中文国产一区发布| 成人国产麻豆网| 最后的刺客免费高清国语| 午夜影院在线不卡| 在线天堂最新版资源| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人片av| 黄色配什么色好看| 色视频在线一区二区三区| 成年女人在线观看亚洲视频| 九九久久精品国产亚洲av麻豆| 久久99热6这里只有精品| 午夜福利网站1000一区二区三区| 国产在视频线精品| 人人妻人人澡人人爽人人夜夜| 成人漫画全彩无遮挡| 大又大粗又爽又黄少妇毛片口| 亚洲精品日韩av片在线观看| 久久99热6这里只有精品| 午夜福利网站1000一区二区三区| 中文字幕精品免费在线观看视频 | 肉色欧美久久久久久久蜜桃| 精品卡一卡二卡四卡免费| 欧美激情 高清一区二区三区| 国产男人的电影天堂91| 久久久久久久亚洲中文字幕| 日日撸夜夜添| 十八禁高潮呻吟视频| 国产精品无大码| 下体分泌物呈黄色| 欧美精品一区二区大全| 日韩亚洲欧美综合| 亚洲欧美精品自产自拍| 国产免费一区二区三区四区乱码| 少妇高潮的动态图| 亚洲,一卡二卡三卡| a级片在线免费高清观看视频| 中文精品一卡2卡3卡4更新| 亚洲图色成人| 国产极品粉嫩免费观看在线 | 两个人的视频大全免费| 国产精品嫩草影院av在线观看| 久久久国产精品麻豆| 色视频在线一区二区三区| 80岁老熟妇乱子伦牲交| 亚洲内射少妇av| 午夜福利视频在线观看免费| 我的女老师完整版在线观看| 涩涩av久久男人的天堂| 亚洲欧美成人综合另类久久久| 成人二区视频| 成人黄色视频免费在线看| 一边摸一边做爽爽视频免费| 伊人亚洲综合成人网| 久久国内精品自在自线图片| 精品熟女少妇av免费看| 久久国内精品自在自线图片| 中文字幕人妻丝袜制服| 国语对白做爰xxxⅹ性视频网站| 少妇高潮的动态图| 少妇猛男粗大的猛烈进出视频| 99国产精品免费福利视频| 国产高清国产精品国产三级| 在线播放无遮挡| 夫妻午夜视频| 国产午夜精品一二区理论片| 亚洲内射少妇av| 一级爰片在线观看| 国产一区亚洲一区在线观看| 在线观看一区二区三区激情| xxx大片免费视频| 亚洲少妇的诱惑av| 99久久中文字幕三级久久日本| 蜜桃久久精品国产亚洲av| 69精品国产乱码久久久| 久久久久久久大尺度免费视频| 街头女战士在线观看网站| 少妇熟女欧美另类| 日本与韩国留学比较| 国产黄色视频一区二区在线观看| 日韩av不卡免费在线播放| 国产片特级美女逼逼视频| 天堂中文最新版在线下载| 日韩免费高清中文字幕av| 久久精品人人爽人人爽视色| 久久青草综合色| 国产毛片在线视频| 久久人人爽av亚洲精品天堂| 亚洲久久久国产精品| 美女xxoo啪啪120秒动态图| 97超碰精品成人国产| 亚洲内射少妇av| 美女中出高潮动态图| 777米奇影视久久| 国产成人精品一,二区| 国产精品国产三级国产av玫瑰| av黄色大香蕉| 99久久中文字幕三级久久日本| 极品少妇高潮喷水抽搐| 少妇人妻 视频| 亚洲av不卡在线观看| 哪个播放器可以免费观看大片| 亚洲av综合色区一区| 蜜臀久久99精品久久宅男| 国产日韩欧美视频二区| 如日韩欧美国产精品一区二区三区 | 免费少妇av软件| 99久久中文字幕三级久久日本| 狠狠精品人妻久久久久久综合| 交换朋友夫妻互换小说| 免费不卡的大黄色大毛片视频在线观看| 日韩人妻高清精品专区| 成人国语在线视频| 九色亚洲精品在线播放| .国产精品久久| 日韩av免费高清视频| 丝袜脚勾引网站| 国产伦理片在线播放av一区| 日韩制服骚丝袜av| 秋霞在线观看毛片| 99精国产麻豆久久婷婷| 一级毛片电影观看| 一边亲一边摸免费视频| 国产一级毛片在线| 久久人妻熟女aⅴ| 免费观看性生交大片5| 久久久国产欧美日韩av| 欧美日韩视频精品一区| a级毛片黄视频| 大陆偷拍与自拍| 亚洲,一卡二卡三卡| 新久久久久国产一级毛片| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 国产色爽女视频免费观看| 国产黄频视频在线观看| 一级毛片我不卡| 免费av不卡在线播放| av天堂久久9| 十八禁高潮呻吟视频| 日本色播在线视频| 亚洲人与动物交配视频| 欧美三级亚洲精品| 国产精品一区二区在线不卡| 日本91视频免费播放| 国产一区二区在线观看日韩| 亚洲国产精品成人久久小说| 欧美+日韩+精品| 精品久久久久久久久av| 日日啪夜夜爽| 日本与韩国留学比较| 简卡轻食公司| 久久久久久久久久人人人人人人| 国产亚洲最大av| 亚洲丝袜综合中文字幕| 久久人人爽人人爽人人片va| 国产日韩欧美亚洲二区| 成人手机av| 久久国内精品自在自线图片| av视频免费观看在线观看| 一区在线观看完整版| av卡一久久| 99re6热这里在线精品视频| 草草在线视频免费看| 亚洲av免费高清在线观看| 日韩在线高清观看一区二区三区| 亚洲精品乱久久久久久| 久久人人爽人人片av| 国产精品99久久久久久久久| 一级片'在线观看视频| 狂野欧美白嫩少妇大欣赏| 精品酒店卫生间| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线 | 国国产精品蜜臀av免费| 乱码一卡2卡4卡精品| 亚洲精品久久久久久婷婷小说| 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 亚洲精品自拍成人| 日韩亚洲欧美综合| 伦理电影大哥的女人| 亚洲精品久久午夜乱码| 考比视频在线观看| www.色视频.com| 香蕉精品网在线| 汤姆久久久久久久影院中文字幕| 欧美xxxx性猛交bbbb| av专区在线播放| av播播在线观看一区| 午夜老司机福利剧场| 国产白丝娇喘喷水9色精品| 亚洲国产色片| 免费看光身美女| 一级二级三级毛片免费看| 精品人妻在线不人妻| 有码 亚洲区| 精品一区二区三卡| 亚洲精品乱久久久久久| 自线自在国产av| 亚洲欧洲国产日韩| 午夜精品国产一区二区电影| 日本色播在线视频| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 97精品久久久久久久久久精品| 男的添女的下面高潮视频| 一区二区三区乱码不卡18| 十八禁网站网址无遮挡| 国产男女内射视频| 亚洲精品一区蜜桃| 国产精品偷伦视频观看了| 精品一区二区三区视频在线| 国产精品一区二区在线观看99| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 国产极品粉嫩免费观看在线 | 亚洲第一av免费看| 亚洲欧洲国产日韩| 在线观看www视频免费| 中文天堂在线官网| 人妻人人澡人人爽人人| 成年人免费黄色播放视频| 狂野欧美激情性xxxx在线观看| 午夜影院在线不卡| 日韩不卡一区二区三区视频在线| 国产日韩欧美视频二区| 我的老师免费观看完整版| 热re99久久国产66热| 亚洲久久久国产精品| 男女国产视频网站| 日本欧美视频一区| 久久久久久久久大av| 中文字幕制服av| 成人手机av| 99热这里只有是精品在线观看| 亚洲精品久久久久久婷婷小说| av女优亚洲男人天堂| 精品人妻偷拍中文字幕| 99热全是精品| 久久99蜜桃精品久久| 欧美3d第一页| 国产精品一国产av| 亚洲国产成人一精品久久久| 免费播放大片免费观看视频在线观看| 中文字幕亚洲精品专区| 亚洲欧洲国产日韩| 亚洲少妇的诱惑av| 亚州av有码| 久久久久久久亚洲中文字幕| 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 一本色道久久久久久精品综合| 涩涩av久久男人的天堂| 极品人妻少妇av视频| 黑人巨大精品欧美一区二区蜜桃 | 亚洲激情五月婷婷啪啪| 亚洲综合精品二区| 欧美人与善性xxx| 精品一区在线观看国产| 精品一区二区三卡| 国产亚洲午夜精品一区二区久久| 国产免费又黄又爽又色| 我的女老师完整版在线观看| 女的被弄到高潮叫床怎么办| 18在线观看网站| 丝袜脚勾引网站| 亚洲成色77777| 日韩成人av中文字幕在线观看| 日韩一区二区视频免费看| 我的女老师完整版在线观看| 午夜精品国产一区二区电影| freevideosex欧美| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区黑人 | 亚洲欧美一区二区三区黑人 | 久久久久久人妻| 观看av在线不卡| 国产69精品久久久久777片| 满18在线观看网站| 中文乱码字字幕精品一区二区三区| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 午夜福利在线观看免费完整高清在| 国产精品熟女久久久久浪| 两个人的视频大全免费| 国产精品一区www在线观看| 亚洲精品乱码久久久久久按摩| 亚洲精品第二区| 91国产中文字幕| 如何舔出高潮| 亚洲av成人精品一二三区| 乱码一卡2卡4卡精品| 精品久久蜜臀av无| 精品一区二区三卡| 欧美另类一区| 伦精品一区二区三区| 高清不卡的av网站| 亚洲人成77777在线视频| 日韩免费高清中文字幕av| 大话2 男鬼变身卡| 久久精品夜色国产| 成人手机av| 中文字幕免费在线视频6| 亚洲伊人久久精品综合| tube8黄色片| 国产精品免费大片| av天堂久久9| 国产精品国产av在线观看| 成人漫画全彩无遮挡| 国产又色又爽无遮挡免| 搡老乐熟女国产| 国产又色又爽无遮挡免| 日韩成人av中文字幕在线观看| 欧美激情极品国产一区二区三区 | 久久99热这里只频精品6学生| 国产成人精品福利久久| 啦啦啦视频在线资源免费观看| 美女脱内裤让男人舔精品视频| 欧美 亚洲 国产 日韩一| 亚洲精品乱码久久久久久按摩| 啦啦啦啦在线视频资源| 69精品国产乱码久久久| 精品亚洲乱码少妇综合久久| 精品午夜福利在线看| 亚洲人成网站在线观看播放| 青青草视频在线视频观看| 成人免费观看视频高清| 99久久精品国产国产毛片| 久久久亚洲精品成人影院| 热re99久久精品国产66热6| 精品国产乱码久久久久久小说| 在线观看一区二区三区激情| 一级毛片电影观看| 午夜福利,免费看| 国精品久久久久久国模美| 亚洲精品自拍成人| 天天躁狠狠躁夜夜躁狠狠躁| 国产av国产精品国产| 久久精品成人免费网站| 国产精品麻豆人妻色哟哟久久| 久久 成人 亚洲| 国产成人影院久久av| 久久国产精品男人的天堂亚洲| 久久国产亚洲av麻豆专区| 免费人妻精品一区二区三区视频| 丁香欧美五月| 亚洲中文日韩欧美视频| 欧美精品一区二区免费开放| 精品一区二区三卡| 黑人操中国人逼视频| 波多野结衣一区麻豆| 午夜福利一区二区在线看| 我要看黄色一级片免费的| 18禁黄网站禁片午夜丰满| 亚洲av美国av|