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

    沼液還田對(duì)旱地紅壤微生物群落代謝與多樣性的影響

    2016-11-28 08:14:18鄭學(xué)博樊劍波祝振球何園球
    生態(tài)學(xué)報(bào) 2016年18期
    關(guān)鍵詞:田量沼液碳源

    鄭學(xué)博,樊劍波,崔 鍵,徐 磊,祝振球,周 靜,3,4,*,何園球

    1 中國科學(xué)院南京土壤研究所, 南京 210008 2 中國科學(xué)院大學(xué), 北京 100049 3 國家紅壤改良工程技術(shù)研究中心,中國科學(xué)院紅壤生態(tài)實(shí)驗(yàn)站, 鷹潭 335211 4 江西省科學(xué)院生物資源研究所, 南昌 330029

    ?

    沼液還田對(duì)旱地紅壤微生物群落代謝與多樣性的影響

    鄭學(xué)博1,2,樊劍波1,崔 鍵1,徐 磊1,2,祝振球1,2,周 靜1,2,3,4,*,何園球1,2

    1 中國科學(xué)院南京土壤研究所, 南京 210008 2 中國科學(xué)院大學(xué), 北京 100049 3 國家紅壤改良工程技術(shù)研究中心,中國科學(xué)院紅壤生態(tài)實(shí)驗(yàn)站, 鷹潭 335211 4 江西省科學(xué)院生物資源研究所, 南昌 330029

    通過沼液還田定位實(shí)驗(yàn),按照不同沼液全氮還田比例設(shè)6個(gè)等氮量(N-P2O5-K2O量均為120-90-135 kg /hm2(對(duì)照除外))處理:對(duì)照(不施肥,CK)、100%化學(xué)氮(NPK)、15%沼液氮+85%化學(xué)氮(BS15)、30%沼液氮+70%化學(xué)氮(BS30)、45%沼液氮+55%化學(xué)氮(BS45)和100%沼液氮(BS100),運(yùn)用Biolog-ECO技術(shù)分析0—20cm花生收獲期土壤微生物群落代謝功能多樣性,闡明微生物群落代謝與沼液還田量的相關(guān)關(guān)系。結(jié)果表明:①BS45、BS30處理土壤微生物群落碳源代謝強(qiáng)度(AWCD)顯著高于CK和NPK處理;而BS15、BS100處理土壤微生物群落碳源代謝強(qiáng)度(AWCD)與CK和NPK處理則無顯著差異;②土壤微生物群落碳源代謝強(qiáng)度(AWCD)、豐富度指數(shù)、Shannon指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)均表現(xiàn)為BS45>BS30> NPK > CK >BS100>BS15;③結(jié)合主成分分析和聚類分析,表明各處理土壤微生物群落功能多樣性分為4組:BS45、BS30處理為一組,微生物群落代謝活性最強(qiáng),特別是碳水化合物、氨基酸、聚合物和胺類等碳源的代謝能力;NPK、CK、BS100處理為一組,微生物群落代謝活性次之;BS15處理為一組,微生物群落代謝能力最低,其碳水化合物、羧酸、氨基酸、聚合物、酚類和胺類等碳源的代謝能力均為最低。結(jié)合主成分分析綜合得分,土壤微生物群落代謝和多樣性的順序?yàn)锽S45>BS30> NPK > CK >BS100>BS15??梢?沼液還田顯著影響旱地紅壤微生物群落的代謝活性和多樣性,沼液不能完全替代化肥,當(dāng)沼液全氮還田比例在30%—45%時(shí),微生物群落代謝活性最強(qiáng),有利于土壤質(zhì)量提高,適于在我國旱地紅壤地區(qū)推廣。

    沼液;旱地紅壤;Biolog-ECO;微生物代謝特征;群落多樣性

    土壤生物是土壤生態(tài)系統(tǒng)的核心,而土壤微生物是土壤生物的重要組成部分,在土壤結(jié)構(gòu)形成、有機(jī)物質(zhì)分解與轉(zhuǎn)化、營養(yǎng)元素循環(huán)與利用、作物生長和生態(tài)環(huán)境改善、植物病害抑制或拮抗、溫室氣體排放等方面起著重要的作用[1-4],是土壤肥力的重要調(diào)控因子和土壤地力的關(guān)鍵驅(qū)動(dòng)因子之一[5]??擅舾蟹从惩寥蕾|(zhì)量變化、表征土壤有機(jī)質(zhì)代謝旺盛程度,是評(píng)價(jià)土壤質(zhì)量的主要指標(biāo)之一。而其功能多樣性信息對(duì)于明確不同環(huán)境中微生物群落的作用具有重要意義[6]。

    Biolog法是描述微生物群落功能變化的重要指標(biāo),其應(yīng)用始于Garland 和 Mills[7]。以ECO板碳源利用為基礎(chǔ),是描述微生物群落功能多樣性的一種簡(jiǎn)單、快速的方法[8-9],其反應(yīng)特征能代表實(shí)際土壤微生物群落對(duì)ECO板底物碳源利用的動(dòng)力學(xué)特征[10-11]。廣泛應(yīng)用于評(píng)價(jià)不同土壤的微生物群落功能多樣性:不同退耕年限[5];不同灌溉水礦化度和施氮量[1];不同有機(jī)氯農(nóng)藥污染[12];不同肥料類型[13]等;可以發(fā)現(xiàn)肥料類型、土壤類型等因素不同,微生物群落結(jié)構(gòu)多樣性不同。

    沼液是畜禽糞便、作物秸稈、工業(yè)垃圾、城市垃圾、能源植物植物等厭氧發(fā)酵產(chǎn)生的終端產(chǎn)物之一,含有大量環(huán)境有害元素(CODCr,BOD5等),無序排放可導(dǎo)致水體嚴(yán)重富營養(yǎng)化,正在成為危害性和危險(xiǎn)性不斷增大的一種污染源,因此其如何處理已經(jīng)成為當(dāng)前水環(huán)境保護(hù)中的關(guān)鍵問題[14]。另一方面,沼液中有機(jī)質(zhì)、大量元素(N、P、K)、微量元素(B、Mn等)等含量豐富[15],且存在形式多為易被植物吸收利用的速效養(yǎng)分,是一種緩速兼?zhèn)涞囊后w有機(jī)肥[14]。因此,沼液被認(rèn)為是一種良好的土壤改良劑或有機(jī)肥料[16]。前人關(guān)于沼液還田的研究主要集中在作物產(chǎn)量、肥料利用效率、土壤C, N礦化、微生物數(shù)量等,而且主要以單施沼液為主。Galvez A 等研究發(fā)現(xiàn)合理的沼液農(nóng)用可以改善土壤質(zhì)量,平衡土壤中各營養(yǎng)元素[17];Abubaker J 等研究發(fā)現(xiàn)沼液農(nóng)用提高了小麥產(chǎn)量、土壤有機(jī)氮礦化和氨氮氧化潛力[18];同時(shí),沼液農(nóng)用提高了土壤微生物活性[19];溫曉霞等的研究表明沼液還田可以提高蔬菜、水果產(chǎn)量[20];馮偉等發(fā)現(xiàn)沼液還田可以提高小麥根際土壤微生物數(shù)量和酶活性[21]。但關(guān)于沼液還田,尤其是沼液化肥配施對(duì)土壤微生物群落功能多樣性影響的研究卻鮮見報(bào)道。

    因此,本文采用Biolog技術(shù)研究不同沼液還田量下0—20 cm土壤微生物群落功能多樣性變化,旨在探討土壤微生物功能多樣性對(duì)不同沼液還田量的響應(yīng)、定量沼液還田對(duì)土壤微生物指標(biāo)和土壤質(zhì)量的影響,從而為旱地紅壤質(zhì)量改善、畜禽廢棄物資源化利用及土壤資源的可持續(xù)利用提供理論依據(jù)和科學(xué)指導(dǎo)。

    1 材料與方法

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

    試驗(yàn)區(qū)位于江西省余江縣劉家站魯王村(116°5′E, 28°12′N)。區(qū)內(nèi)年平均降雨量1750 mm,80%的降雨集中在3—6月,年蒸發(fā)量1350 mm,年平均氣溫18℃,無霜期267 d左右,屬亞熱帶濕潤季風(fēng)氣候。

    試驗(yàn)田為旱平地,試驗(yàn)時(shí)間2013年4月—2014年8月,前茬作物為花生。試驗(yàn)設(shè)置6個(gè)處理,每個(gè)處理3次重復(fù),各小區(qū)面積12.4 m×4.8 m=60 m2,小區(qū)間田埂寬50 cm,高25 cm,保護(hù)行20 cm。供試土壤為第四紀(jì)紅粘土發(fā)育的典型紅壤。0—20cm土壤基本理化性質(zhì)見表1。

    表1 0—20cm土壤養(yǎng)分含量

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

    各處理N-P2O5-K2O量均為120-90-135 kg/hm2(CK除外)。具體為:Ⅰ對(duì)照(不施肥,CK);Ⅱ單施化學(xué)氮磷鉀肥(NPK);Ⅲ 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 15%,BS15);Ⅳ化學(xué)氮磷鉀肥+沼液(沼液全N占總N 30%,BS30);Ⅴ化學(xué)氮磷鉀肥+沼液(沼液全N占總N 45%,BS45);Ⅵ 單施沼液(BS100)。具體施肥量見表2。

    表2 各試驗(yàn)處理施肥量

    a:N-P2O5-K2O 分別來自沼液/化肥;CK: 不施肥 No fertilizer; NPK: 單施化學(xué)氮磷鉀肥 Chemical fertilizers only;BS15: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 15%) 85% chemical fertilizer total nitrogen plus 15% biogas slurry total nitrogen; BS30: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 30%) 70% chemical fertilizer total nitrogen plus 30% biogas slurry total nitrogen; BS45: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 45%) 55% chemical fertilizer total nitrogen plus 45% biogas slurry total nitrogen; BS100: 單施沼液 Biogas slurry only

    氮肥為尿素(N 46%),磷肥為鈣鎂磷肥(P2O512%),鉀肥為氯化鉀(K2O 60%),沼液發(fā)酵原材料主要是豬糞尿等,發(fā)酵3個(gè)月以上,細(xì)膩均勻,透明無明顯糞臭味,pH 7.67,具體養(yǎng)分含量見表3。

    表3 供試沼液養(yǎng)分含量

    所有肥料全部基施,沼液還田處理中不足的P2O5、K2O量用化肥補(bǔ)足(表2)。2013年播種前整地規(guī)劃小區(qū)、翻耕,播種前2d澆灌沼液,尿素、鈣鎂磷肥、氯化鉀播種前1d二次翻耕翻入。種植制度為單季花生/a,每年4月中旬播種,8月中旬收獲。試種品種贛花1號(hào),密度14 萬株/hm2,作物生長期間不再進(jìn)行追肥。

    1.3 土壤樣品采集

    取樣時(shí)間為2014年8月中旬花生收獲前2d。

    取樣方法采用S形取樣,各小區(qū)選取5個(gè)樣點(diǎn),用土鉆(直徑5.0cm)采集0—20 cm土層土樣,先混勻然后去除其中的雜物、砂礫及作物細(xì)根等。過2mm篩后用于測(cè)定土壤微生物功能多樣性(4℃保存)。

    1.4 試驗(yàn)方法

    1.4.1 土壤微生物功能多樣性分析

    采用Biolog-ECO微平板法對(duì)不同沼液還田量處理土壤微生物群落的功能多樣性進(jìn)行測(cè)定。每個(gè)ECO板含31種碳源,分碳水化合物(10種)、羧酸(7種)、氨基酸(6種)、聚合物(4種)、酚類(2種)、胺類(2種)6種類型[22]。試驗(yàn)步驟:稱取1g鮮土,加入到盛有9mL 0.85% NaCl無菌溶液的18mL試管中,封口后在微型振動(dòng)儀震蕩5min,按10倍稀釋法用0.85% NaCl無菌溶液將其稀釋到原來的10-4用于接種。將加樣液傾倒在無菌加樣槽中,用8通道加樣器將每孔加樣150uL,每樣1板,3次重復(fù),將接種好的ECO板(EcoPlatesTM,美國Matrix Technologies Corporation 生產(chǎn))置于25℃暗箱連續(xù)培養(yǎng),期間每隔24h用Biolog自動(dòng)鑒定系統(tǒng)(MicroStationTMSystem,美國Biolog公司)在590nm下讀數(shù),直到光密度值穩(wěn)定為止,總共培養(yǎng)時(shí)間240h。

    1.4.2 計(jì)算方法

    土壤微生物利用碳源的整體能力用平均吸光度值(AWCD)表示:

    (1)

    式中,Ci為第i孔的吸光度值;R為對(duì)照孔的吸光度值;Ci-R<0時(shí)記為0。

    土壤微生物利用ECO板中六類碳源的能力用各自的平均吸光度值(AWCD′)[23]表示:

    (2)

    式中,Ci為i孔的吸光度值;R為對(duì)照孔的吸光度值;n為所選擇碳源類型的個(gè)數(shù);Ci-R<0時(shí)記為0。

    土壤微生物群落功能多樣性分析:本研究采用培養(yǎng)120h的AWCD和AWCD′值進(jìn)行多樣性指數(shù)計(jì)算、主成分分析和聚類分析。土壤微生物群落功能多樣性指數(shù)包括Shannon指數(shù)(H)、豐富度指數(shù)(S)、Shannon均勻度指數(shù)(E)、Simpson優(yōu)勢(shì)度指數(shù)(D)。具體計(jì)算公式如下[24-25]:

    Shannon指數(shù)

    (3)

    豐富度指數(shù)S:被利用碳源的總數(shù)目(判斷標(biāo)準(zhǔn)為Ci-R>0.25);

    Shannon均勻度指數(shù)

    E=H/lnS

    (4)

    Simpson 優(yōu)勢(shì)度指數(shù)

    (5)

    式中,Ci為第i孔的吸光度值;R為對(duì)照孔的吸光度值;n為ECO板碳源數(shù)量,n=31;Pi為第i孔的相對(duì)吸光度值與整個(gè)ECO板相對(duì)吸光度值總和的比率;S為顏色變化的碳源的數(shù)量。

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

    試驗(yàn)所得數(shù)據(jù)采用Excel 2007進(jìn)行處理、繪制圖表,SPSS 13.0 統(tǒng)計(jì)分析軟件進(jìn)行單因素方差分析(One-way ANOVA)、相關(guān)性分析(大小采用pearson指數(shù))、主成分分析和聚類分析,差異顯著性分析采用Duncan新復(fù)極差方法,顯著性水平P<0.05。

    2 結(jié)果與分析

    2.1 土壤微生物利用碳源動(dòng)力學(xué)分析

    2.1.1 不同沼液還田量下土壤微生物群落整體代謝活性變化(AWCD)

    平均吸光值A(chǔ)WCD反映了土壤微生物群落利用碳源的整體能力,是反映微生物代謝活性的有效指標(biāo)[7,26]。圖1所示,不同沼液還田量處理土壤微生物群落代謝AWCD值的變化趨勢(shì)均為隨著時(shí)間的延長而升高。0—48h:各處理AWCD值均較低,說明土壤微生物群落代謝活性較低,碳源基本未被利用;48—168h:各處理AWCD值迅速升高,微生物呈指數(shù)增長,說明土壤微生物從48h后開始大量利用碳源,土壤微生物代謝活性顯著升高。168h后:各處理AWCD值基本趨于平緩,此時(shí)土壤微生物活性達(dá)到穩(wěn)定??傮w上,不同沼液還田處理下土壤微生物利用31種碳源的綜合能力(AWCD)為:BS30>BS45> NPK > CK >BS100> BS15。BS30和BS45處理120h 的微生物活性顯著高于BS15處理。

    圖1 土壤微生物群落平均吸光值的變化Fig.1 Changes of AWCD of soil edaphon against incubation time under biogas slurry application CK: 不施肥 No fertilizer; NPK: 單施化學(xué)氮磷鉀肥 Chemical fertilizers only;BS15: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 15%) 85% chemical fertilizer total nitrogen plus 15% biogas slurry total nitrogen; BS30: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 30%) 70% chemical fertilizer total nitrogen plus 30% biogas slurry total nitrogen; BS45: 化學(xué)氮磷鉀肥+沼液(沼液全N占總N 45%) 55% chemical fertilizer total nitrogen plus 45% biogas slurry total nitrogen; BS100: 單施沼液Biogas slurry only

    2.1.2 不同沼液還田量下土壤微生物群落各類型碳源利用特征

    按照化學(xué)基團(tuán)性質(zhì),31種碳源分碳水化合物(10種)、羧酸(7種)、氨基酸(6種)、聚合物(4種)、酚類(2種)、胺類(2種)等共六大類碳源。分析各類碳源的AWCD值可判斷各處理土壤微生物群落對(duì)特定類型碳源的利用能力。本文選取120h時(shí)的AWCD值繪制不同沼液還田量下土壤微生物碳源利用能力圖(圖2)。各處理之間土壤微生物群落代謝碳水化合物、羧酸、聚合物的能力存在顯著差異(P<0.05),代謝氨基酸、酚類、胺類的能力無顯著差異(P<0.05)。其中,碳水化合物類:BS45>BS30> NPK≈CK >BS100>BS15;羧酸類:BS30> NPK >BS100>BS45> CK >BS15;氨基酸類:BS45>BS30> CK > NPK >BS100>BS15;聚合物類:BS45>BS30 ≈ NPK >BS100> CK >BS15;酚類:BS45>BS30> CK > NPK ≈ BS100>BS15;胺類:BS45>BS30> CK >BS15> NPK >BS100??梢园l(fā)現(xiàn),土壤微生物碳水化合物、氨基酸、酚類AWCD值與上述31種碳源總體利用情況類似。BS30處理土壤微生物對(duì)6類碳源的代謝能力高于除BS45處理以外的各處理,而BS45處理較其它處理則提高了土壤微生物對(duì)除羧酸類以外的五類碳源的代謝能力。

    圖2 不同沼液還田量下120h時(shí)土壤微生物不同碳源AWCDFig.2 AWCD of soil edaphonat 120h of different carbon sources under biogas slurry application

    2.2 不同沼液還田量下土壤微生物群落功能多樣性指數(shù)和豐富度

    Shannon 物種指數(shù)、Shannon 均勻度指數(shù)、Simpson 優(yōu)勢(shì)度指數(shù)和碳源利用豐富度指數(shù)可以反映土壤微生物群落的功能多樣性。不同的多樣性指數(shù)反映土壤微生物群落功能多樣性的不同方面,Shannon指數(shù)反映群落物種及其個(gè)體數(shù);Shannon均勻度指數(shù)反映群落個(gè)體分布的均勻制度;Simpon優(yōu)勢(shì)度指數(shù)反映群落中最常見物種的優(yōu)勢(shì)度[27]。各處理ECO板培養(yǎng)120h后的土壤微生物群落功能多樣性指數(shù)見表2。豐富度指數(shù)、Shannon指數(shù)、Simpson 優(yōu)勢(shì)度指數(shù)方面均為BS45處理最高,BS30處理次之,BS15處理最低。而Shannon 均勻度指數(shù)方面則為BS15處理最高,顯著高于其它處理,而剩余處理間無顯著差異??梢?BS30和BS45處理增加了土壤中微生物的種類和優(yōu)勢(shì)種群的數(shù)量,但并沒有提高微生物種類的平均分布;而BS15處理抑制了土壤中微生物的種類和優(yōu)勢(shì)種群的數(shù)量,但提高了土壤中存在的微生物種類的均勻分布。

    表4 沼液還田對(duì)土壤微生物群落功能多樣性指數(shù)的影響

    同列不同小寫字母表示處理間差異顯著(P<0.05)

    2.3 不同沼液還田量下土壤微生物群落功能主成分分析(PCA)

    土壤微生物群落功能多樣性反映了群落總體的變化,但未能反映微生物群落代謝的詳細(xì)信息,研究土壤微生物對(duì)不同碳源利用能力的差異,利于更全面地研究微生物群落代謝功能的特性。以培養(yǎng)120h不同處理微生物群落對(duì)ECO板31種碳源利用情況做主成分分析(PCA),31個(gè)成分因子中前8個(gè)主成分(累積方差貢獻(xiàn)率>85%,特征根>1)的累積方差貢獻(xiàn)率達(dá)到92.70%。從中提取可以聚集單一碳源變量的數(shù)據(jù)變異的第1主成分(PC1,方差貢獻(xiàn)率32.34%,特征根10.02)和第2主成分(PC2,方差貢獻(xiàn)率15.36%,特征根4.76)來分析土壤微生物群落功能多樣性。第3至第8主成分的方差貢獻(xiàn)率較小,分別為12.76%、10.85%、7.93%、5.88%、4.16%、3.43%。

    不同沼液還田量下土壤微生物群落碳源利用的PCA排序圖如下(圖3)。由圖3可知,前兩個(gè)主成分共解釋了總變異的47.70%。結(jié)合方差分析結(jié)果表明,各處理在PC1上存在明顯的空間分異,6種處理的土壤微生物群落可明顯分為4簇:BS30、BS45處理聚為一簇,CK、NPK處理聚為一簇、BS100處理為一簇、BS15處理為一簇。而在PC2上則沒有顯著的空間差異。PC1方差分析表明,來自(BS30、BS45)和BS15處理的土壤微生物群落結(jié)構(gòu)存在顯著差異。同時(shí),在PC1上,各處理土壤微生物群落碳源利用主成分分析得分為BS30>BS45> NPK > CK >BS100>BS15。這與Shannon指數(shù)、豐富度指數(shù)、Simpson指數(shù)的研究結(jié)果一致。

    圖3 不同沼液還田量下土壤微生物群落功能主成分分析Fig.3 The principal component analysis for soil microbial community function under biogas slurry application

    原始變量與各主成分的相關(guān)系數(shù)可用載荷值表征,載荷值越大,表明該原始變量是該主成分的主要影響因子,對(duì)該主成分的貢獻(xiàn)越大[28]。表5為31種碳源(原始變量)在前2個(gè)主成分上的載荷值。從表中可以看到,對(duì)PC1貢獻(xiàn)大的碳源(載荷值≥0.50或≤-0.50)有16種,決定了主成分1的變異,其中碳水化合物4種(β-甲基-D-葡萄糖苷,D-甘露醇,葡萄糖-1-磷酸鹽,D, L-a-甘油),羧酸2種(D-葡萄胺酸,衣康酸),氨基酸4種(L-精氨酸,L-天冬酰胺酸,L-苯基丙氨酸,L-蘇氨酸),聚合物3種(吐溫 40,吐溫80,肝糖),酚類2種(2-羥基苯甲酸,4-羥基苯甲酸),胺類1種(腐胺);對(duì)PC2貢獻(xiàn)大的碳源(載荷值≥0.50或≤-0.50)有8種,決定了主成分2的變異,其中碳水化合物2種(D-木糖,N-乙?;?D-葡萄胺,),羧酸4種(D-葡萄胺酸,D-半乳糖醛酸,y-羥基丁酸,D-蘋果酸),氨基酸1種(甘氨酰-L-谷氨酸),聚合物1種(a環(huán)式糊精)。結(jié)合圖3,發(fā)現(xiàn)較其他處理,BS30和BS45顯著提高了碳水化合物(4種)、羧酸(2種)、氨基酸(4種)、聚合物(3種)、酚類(2種)、胺類(1種)等共16種碳源的代謝能力。

    表5 ECO板上31種碳源在第1、2、3主成分上的載荷值

    圖4 不同沼液還田量處理碳源類型主成分得分 Fig.4 Pricipal component scores of carbon sources under biogas slurryapplication

    對(duì)所提取的8個(gè)主成分分別計(jì)算各處理6類碳源的主成分得分 (Pricipal component scores, F ),結(jié)果如圖4。F值的相對(duì)大小可以表征不同處理土壤微生物群落對(duì)同一類碳源的代謝能力和同一處理土壤微生物對(duì)不同類型碳源代謝能力的高低。對(duì)F值進(jìn)行橫向比較發(fā)現(xiàn),不同處理土壤微生物對(duì)不同碳源的利用能力不同。CK:代謝碳水化合物能力最強(qiáng);NPK:代謝碳水化合物、羧酸、聚合物能力最強(qiáng);BS30、BS45、BS100:代謝碳水化合物、羧酸、氨基酸、聚合物能力較強(qiáng);BS15:代謝碳水化合物和羧酸能力較強(qiáng);各處理土壤微生物群落均為代謝酚類和胺類能力較弱。對(duì)F值進(jìn)行縱向比較可知,代謝碳水化合物、羧酸、氨基酸、聚合物、酚類和胺類能力最弱的都是來自BS15處理的土壤微生物,代謝以上物質(zhì)最強(qiáng)的都是來自BS30和BS45處理的土壤微生物。

    同時(shí)對(duì)所提取的8個(gè)主成分分別計(jì)算各處理6類碳源的主成分綜合得分 ( Pricipal component composite scores, ZF ) (圖5)。ZF值的相對(duì)大小可以反映不同處理土壤微生物對(duì)ECO板31種碳源綜合代謝能力的高低。圖5顯示,6種處理土壤微生物對(duì)ECO板中31種單一碳源的綜合代謝能力排序?yàn)椋築S45>BS30> NPK >CK >BS100>BS15。

    2.4 土壤微生物群落功能多樣性聚類分析

    以各處理的前2個(gè)主成分得分代替原指標(biāo)作為評(píng)價(jià)土壤微生物代謝能力的新指標(biāo),以歐式距離作為衡量各處理土壤微生物代謝能力的差異,采用最短距離法將各處理按土壤微生物代謝能力的親疏相似程度進(jìn)行系統(tǒng)聚類。由聚類樹形圖(圖6)可直觀地看出各處理土壤微生物代謝能力的歸類情況(以歐式距離10為分異點(diǎn)):BS30,BS45> CK,NPK >BS100>BS30,即各處理微生物群落代謝能力分為4個(gè)等級(jí):一等為BS30,BS45,屬強(qiáng)土壤微生物群落代謝能力等級(jí);二等為CK,NPK,屬較強(qiáng)土壤微生物群落代謝能力等級(jí);三等為BS100,土壤微生物群落代謝能力低;四等為BS15,土壤微生物群落代謝能力最低。這與主成分分析結(jié)果一致。

    圖5 不同沼液還田量處理碳源類型主成分綜合得分 Fig.5 Pricipal component composite scores of carbon sources underbiogas slurryapplication

    圖6 聚類樹形圖Fig.6 Arborescence of cluster analysis

    3 討論

    研究表明,微生物群落碳源代謝強(qiáng)度和代謝類型的變化能較早預(yù)測(cè)土壤養(yǎng)分的變化趨勢(shì),被認(rèn)為是預(yù)測(cè)土壤質(zhì)量變化最有潛力的敏感性生物指標(biāo)[29-30]。土壤微生物是土壤生態(tài)系統(tǒng)中最重要和最活躍的部分,是驅(qū)動(dòng)土壤有機(jī)物質(zhì)和養(yǎng)分轉(zhuǎn)化與循環(huán)的主要?jiǎng)恿?其活性和多樣性變化在很大程度上反映了土壤有機(jī)碳庫的轉(zhuǎn)化情況[31]。而其數(shù)量、群落結(jié)構(gòu)及活性受外源有機(jī)物投入種類、存在狀態(tài)和數(shù)量的影響[32]。有機(jī)物料作為外源有機(jī)碳源進(jìn)入土壤,改變了土壤微生物的生存環(huán)境和微生物所需能源物質(zhì)的數(shù)量和質(zhì)量[33],從而對(duì)土壤自身有機(jī)碳庫的轉(zhuǎn)化產(chǎn)生影響。劉定輝等[34]發(fā)現(xiàn)秸稈還田循環(huán)利用顯著增加了土壤微生物數(shù)量和生物量。劉國順等[35]、李正等[36]發(fā)現(xiàn)化肥與綠肥配施還田后土壤微生物活性顯著提高,利于有機(jī)物質(zhì)的轉(zhuǎn)化和土壤質(zhì)量的提高。AWCD從功能代謝水平上揭示土壤微生物群落結(jié)構(gòu)多樣性,是反映土壤微生物活性,描述土壤微生物群落功能多樣性的重要指標(biāo)[30]。本文研究表明,較不施肥處理,單施化肥和單施沼液處理土壤微生物群落碳源代謝強(qiáng)度無顯著變化,這與程存剛等[32]的研究結(jié)果一致;而沼液化肥配施處理提高了土壤微生物群落碳源代謝強(qiáng)度,特別是沼液全N配施45%和30%處理土壤微生物群落碳源代謝強(qiáng)度顯著提高,碳水化合物4種(β-甲基-D-葡萄糖苷,D-甘露醇,葡萄糖-1-磷酸鹽,D, L-a-甘油),羧酸2種(D-葡萄胺酸,衣康酸),氨基酸4種(L-精氨酸,L-天冬酰胺酸,L-苯基丙氨酸,L-蘇氨酸),聚合物3種(吐溫 40,吐溫80,肝糖),酚類2種(2-羥基苯甲酸,4-羥基苯甲酸),胺類1種(腐胺)等16種碳源代謝能力增幅最大。說明一定比例沼液與化肥配施可以改善土壤微生物群落的營養(yǎng)供應(yīng),進(jìn)而提高微生物群落代謝能力。原因可能是該比例沼液化肥配施為土壤微生物生長提供所需碳源和能源的同時(shí),使微生物C/N適宜,土壤速效養(yǎng)分含量多,促進(jìn)了花生植株的吸收和根系生長,從而增加了根系分泌物,促進(jìn)了土壤中微生物數(shù)量的增長和活性的提高。另一方面,單施化肥會(huì)導(dǎo)致紅壤進(jìn)一步酸化[23],而土壤中絕大多數(shù)微生物適宜在中性環(huán)境中生存繁殖,沼液的強(qiáng)pH(表2)抑制了這種酸化趨勢(shì)。不同沼液還田量下,土壤微域環(huán)境發(fā)生了改變,這使得多種微生物的適宜性受到影響,從而影響了土壤微生物對(duì)碳源的選擇性利用。

    不同的多樣性指數(shù)反映土壤微生物群落功能多樣性的不同方面。豐富度指數(shù)S、Shannon 均勻度指數(shù)、Simpson 優(yōu)勢(shì)度指數(shù)分別從微生物群落物種豐富度、均一性及常見物種優(yōu)勢(shì)度等方面反映土壤微生物群落多樣性[27]。本研究發(fā)現(xiàn)BS45處理Shannon指數(shù)、豐富度指數(shù)S、Simpson優(yōu)勢(shì)度指數(shù)較其它處理分別提高2.07%—37%,2.94%—61.76%和1.08%—4.26%。說明該處理在增加土壤微生物種群物種種類和數(shù)量的基礎(chǔ)上提高了優(yōu)勢(shì)微生物種群的數(shù)量,提高了土壤微生物群落的多樣性。

    對(duì)各試驗(yàn)處理AWCD值進(jìn)行進(jìn)一步主成分和聚類分析(圖6),結(jié)果表明各試驗(yàn)處理土壤微生物群落多樣性分為4組,BS30、BS45處理微生物群落多樣性相似,為一組;不施肥、單施化肥處理為一組;BS100處理為一組;BS15處理為一組。具體為BS45,BS30> CK,NPK >BS100>BS15。這說明沼液還田量顯著影響土壤微生物群落多樣性,而對(duì)土壤微生物群落結(jié)構(gòu)變化的影響則需要通過PCR-DGGE方法進(jìn)一步分析。本試驗(yàn)僅針對(duì)連續(xù)施沼液兩年后的土壤微生物群落活性及多樣性進(jìn)行研究,而關(guān)于長期沼液農(nóng)用對(duì)土壤微生物群落活性和多樣性的影響以及是否對(duì)土壤質(zhì)量帶來負(fù)面影響等方面有待進(jìn)一步研究。

    4 結(jié)論

    (1)與不施肥、單施化肥處理相比,單施沼液處理對(duì)土壤微生物群落代謝活性影響較小,而沼液化肥配施處理顯著影響土壤微生物群落代謝活性。沼液全N配施比例在30%—45%時(shí),土壤微生物群落整體活性最高,而且碳源利用類型以碳水化合物(4種)、羧酸(2種)、氨基酸(4種)、聚合物(3種)、酚類(2種)和胺類(1種)等共16種碳源為主;

    (2)BS45處理Shannon物種指數(shù)、豐富度指數(shù)、Simpson優(yōu)勢(shì)度指數(shù)較其它處理分別提高2.07%—37%,2.94%—61.76%和1.08%—4.26%,該處理在增加土壤微生物種群物種種類和數(shù)量的基礎(chǔ)上提高了優(yōu)勢(shì)微生物的數(shù)量,同時(shí)提高了土壤微生物群落功能多樣性。

    (3)結(jié)合主成分分析(PCA)和聚類分析得出:微生物群落多樣性分為4個(gè)等級(jí):一等為BS30,BS45處理,屬強(qiáng)土壤微生物群落代謝能力等級(jí);二等為CK,NPK處理,屬較強(qiáng)土壤微生物群落代謝能力等級(jí);三等為BS100處理,土壤微生物群落代謝能力低;四等為BS15處理,土壤微生物群落代謝能力最低。

    綜上,對(duì)于旱地紅壤而言,沼液并不能完全取代化肥單獨(dú)施用,必須與一定比例的化肥配合施用。該研究中,沼液化肥氮配比為30%—45%是最佳配施比例,利于土壤質(zhì)量的保持和提高。

    [1] 閔偉, 侯振安, 冶軍, 馬麗娟, 曹志, 羅慧蘭. 長期咸水滴灌棉田土壤微生物活性及群落功能多樣性. 生態(tài)學(xué)雜志, 2014, 33(11): 2950-2958.

    [2] 鄭亞萍, 吳正峰, 王才斌, 梁曉艷, 陳殿緒, 孫學(xué)武, 王建國, 鄭永美, 馮昊, 孫秀山. 旱地花生不同土壤類型主要土壤肥力指標(biāo)季節(jié)變異及其相互關(guān)系研究. 核農(nóng)學(xué)報(bào), 2013, 27(6): 831-838.

    [3] 林先貴, 胡君利. 土壤微生物多樣性的科學(xué)內(nèi)涵及其生態(tài)服務(wù)功能. 土壤學(xué)報(bào), 2008, 45(5): 892-900.

    [4] Shi W, Miller B E, Stark J M, Norton J M. Microbial nitrogen transformations in responce to treated dairy waste in agricultural soils. Soil Science Society of America Journal, 2004, 68(6): 1867-1874.

    [5] 董立國, 蔣齊, 蔡進(jìn)軍, 張?jiān)礉? 許浩, 李生寶. 基于Biolog-ECO 技術(shù)不同退耕年限苜蓿地土壤微生物功能多樣性分析. 干旱區(qū)研究, 2011, 28(4): 630-637.

    [6] Preston M J, Boddy L, Randerson P F. Analysis of microbial community functional diversity using sole-carbon-source utilization profiles acritique. FEMS Microbial Ecology, 2002, 42(1): 1-14.

    [7] Garland J L, Mills A L. Classification and characterization of heterotrophic microbial communities on basis of patterns of community-lever sole-carbon-source utilization. Applied and Environental Microbiology, 1991, 57: 2351-2359.

    [8] Garland J L. Analysis and interpretation of community-level physiological profiles in microbial ecology. FEMS Microbiology Ecology, 1997, 24(4): 289-300.

    [9] Hollibaugh J T. Relationship between thymidine metabolism, bacterioplankton community metabolic capabilities and sources of organic matter. Microbial Ecology, 1994, 28(2): 117-131.

    [10] 章家恩, 蔡燕飛, 高愛霞, 朱麗霞. 土壤微生物多樣性試驗(yàn)研究方法概述. 土壤, 2004, 36(4): 346-350.

    [11] 陳承利, 廖敏, 曾路生. 污染土壤微生物群落結(jié)構(gòu)多樣性及功能多樣性測(cè)定方法. 生態(tài)學(xué)報(bào), 2006, 26(10): 3404-3412.

    [12] 鄭麗萍, 龍濤, 林玉鎖, 于賜剛, 劉燕, 祝欣. Biolog-ECO解析有機(jī)氯農(nóng)藥污染場(chǎng)地土壤微生物群落功能多樣性特征. 應(yīng)用與環(huán)境學(xué)報(bào), 2013, 19(5): 759-765.

    [13] 羅希茜, 郝曉暉, 陳濤, 鄧嬋娟, 吳金水, 胡榮桂. 長期不同施肥對(duì)稻田土壤微生物群落功能多樣性的影響. 生態(tài)學(xué)報(bào), 2009, 29(2): 740-748.

    [14] 白文波, 王春艷, 李茂松, 鮑順淑, 何爽, 曹肆林, 宋吉青. 不同灌溉條件下保水劑對(duì)新疆棉花生長及產(chǎn)量的影響. 農(nóng)業(yè)工程學(xué)報(bào), 2010, 26(10): 69-76.

    [15] Zirkler D, Peters A, Kaupenjohann M. Elemental composition of biogas residues: Variability and alteration during anaerobic digestion. Biomass and Bioenergy, 2014, 67: 89-98.

    [16] Pinto R, Ilileji K. On the performance of biobased carriers from distillers grains with solubles (DDGS) // ASABE Meeting Presentation, 2009: Paper number095675.

    [17] Galvez A, Sinicco T, Cayuela M L, Mingorance M D, Fornasier F, Mondini C. Short term effects of bioenergy by-products on soil C and N dynamics, nutrient availability and biochemical properties. Agriculture, Ecosystems & Environment, 2012, 160: 3-14.

    [18] Abubaker J, Risberg K, Pell M. Biogas residues as fertilisers-Effects on wheat growth and soil microbial activities. Applied Energy, 2012, 99: 126-134.

    [19] Terhoeven-Urselmans T, Scheller E, Raubuch M, Ludwig B, Joergensen R G. CO2evolution and N mineralization after biogas slurry application in the field and its yield effects on spring barley. Applied Soil Ecology, 2009, 42(3): 297-302.

    [20] 溫曉霞, 殷瑞敬, 高茂盛, 艾綏龍. 不同覆蓋模式下旱作蘋果園土壤酶活性和微生物數(shù)量時(shí)空動(dòng)態(tài)研究. 西北農(nóng)林學(xué)報(bào), 2011, 20(11): 82-88.

    [21] 馮偉, 管濤, 王曉宇, 朱云集, 郭天財(cái). 沼液與化肥配施對(duì)冬小麥根際土壤微生物數(shù)量和酶活性的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2011, 22(4): 1007-1012.

    [22] 李志斐, 王廣軍, 謝駿, 郁二蒙, 余德光, 夏耘, 魏南. 草魚養(yǎng)殖池塘生物膜固著微生物群落碳代謝Biolog分析. 水產(chǎn)學(xué)報(bào), 2014, 38(12): 1985-1994.

    [23] 孫鳳霞, 張偉華, 徐明崗, 張文菊, 李兆強(qiáng), 張敬業(yè). 長期施肥對(duì)紅壤微生物生物量碳氮和微生物碳源利用的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(11): 2792-2798.

    [24] Harch B D, Correll R L, Meech W, Kirkby C A, Pankurst C E. Using the Gini coefficient with Biolog Substrate utilization data to provide an alternative quantitative measure for comparing bacterial soil communities. Journal of Microbiological Methods, 1997, 30(1): 91-101.

    [25] Zak J C, Willig M R, Moorhead D L, Wildman H G. Functional diversity of microbial communities: A quantitative approach. Soil Biology and Biochemistry, 1994, 26(9): 1101-1108.

    [26] Haack S K, Garchow H, Klug M J, Forney LJ. Analysis of factors affecting the accuracy, reproducibility, and interpretation of microbial community carbon source utilization patterns. Applied and Environmental Microbiology, 1995, 61(4): 1458-1468.

    [27] 李忠佩, 吳曉晨, 陳碧云. 不同利用方式下土壤有機(jī)碳轉(zhuǎn)化及微生物群落功能多樣性變化. 中國農(nóng)業(yè)科學(xué), 2007, 40(8): 1712-1721.

    [28] 邱莉萍, 劉軍, 王益權(quán), 孫慧敏, 和文祥. 土壤酶活性與土壤肥力的關(guān)系研究. 植物營養(yǎng)與肥料學(xué)報(bào), 2004, 10(3): 277-280.

    [29] 孫波, 趙其國, 張?zhí)伊? 俞慎. 土壤質(zhì)量與持續(xù)環(huán)境: Ⅲ.土壤質(zhì)量評(píng)價(jià)的生物學(xué)指標(biāo). 土壤, 1997, (5): 225-234.

    [30] 吳等, 宋志文, 徐愛玲, 鄭遠(yuǎn), 夏巖. 青島市不同功能區(qū)冬季空氣微生物群落代謝與多樣性特征. 生態(tài)學(xué)報(bào), 2015, 35(7): 2277-2284.

    [31] 周麗霞, 丁明懋. 土壤微生物學(xué)特性對(duì)土壤健康的指示作用. 生物多樣性, 2007, 15(2): 162-171.

    [32] 程存剛, 趙德英, 呂德國, 姜曼, 杜國棟. 植物源有機(jī)物料對(duì)果園土壤微生物群落多樣性的影響. 植物營養(yǎng)與肥料學(xué)報(bào), 2014, 20(4): 913-922.

    [33] 張海燕, 肖延華, 張旭東, 李軍, 席聯(lián)敏. 土壤微生物量作為土壤肥力指標(biāo)的探討. 土壤通報(bào), 2006, 37(3): 422-425.

    [34] 劉定輝, 舒麗, 陳強(qiáng), 陳尚洪, 陳紅琳, 朱鐘麟. 秸稈還田少免耕對(duì)沖積土微生物多樣性及微生物碳氮的影響. 應(yīng)用與環(huán)境學(xué)報(bào), 2011, 17(2): 158-161.

    [35] 劉國順, 李正, 敬海霞, 葉協(xié)鋒, 史宏志, 王永, 楊超, 常棟. 連年翻壓綠肥對(duì)植煙土壤微生物量及酶活性的影響. 植物營養(yǎng)與肥料學(xué)報(bào), 2010, 16(6): 1472-1478.

    [36] 李正, 劉國順, 敬海霞, 葉協(xié)鋒, 解昌盛, 向永光, 張文平, 楊超, 王永, 習(xí)相銀. 綠肥與化肥配施對(duì)植煙土壤微生物量及供氮能力的影響. 草葉學(xué)報(bào), 2011, 20(6): 126-134.

    Analysis on metabolic characteristics and functional diversity of soil edaphon communities in upland red soil under biogas slurry application

    ZHENG Xuebo1,2, FAN Jianbo1, CUI Jian1, XU Lei1,2, ZHU Zhenqiu1,2, ZHOU Jing1,2,3,4,*, HE Yuanqiu1,2

    1 Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China 2UniversityofChineseAcademyofSciences,Beijing100049,China3NationalEngineeringResearchandTechnologyCenterforRedSoilImprovement,RedSoilEcologicalExperimentStation,ChineseAcademyofSciences,Yingtan335211,China4InstituteofBiologyResource,JiangxiAcademyofSciences,Nanchang330029,China

    In order to provide a theoretical basis for soil edaphon community improvement and a fertilization system establishment for arid areas in south China, a field experiment was conducted to determine the effects of biogas slurry returning systems on the metabolic characteristics and functional diversity of soil edaphon communities in upland red soil.The experiment contained six treatments: nofertilizer (CK), chemical fertilizer (NPK),15% biogas slurry N+85% chemical fertilizer N (BS15), 30% biogas slurry N+70% chemical fertilizer N (BS30), 45% biogas slurry N+55% chemical fertilizer N (BS45), and 100% biogas slurry N (BS100). Each treatment included the same amounts of N-P2O5-K2O=120-90-135 kg/hm2(except for CK).Average well color development (AWCD), indexes of microbial functional diversity (Shannon index, Simpson index, richness index) were determined, and principal component analysis (PCA) and clustering analysis were conducted. The results showed that the AWCD value of treatments BS45 and BS30 was significantly higher than that of CK and NPK treatments, respectively. Treatments BS15 and BS100 showed no significant difference to CK and NPK treatments, respectively. The order of AWCD value, Richness index, Shannon index, and Simpson index wereall the same as follow: BS45>BS30> NPK > CK >BS100>BS15. The results of PCA and clustering analysis showed that the functional diversity of soil edaphoncan develop into 4 groups: First were treatments BS45 and BS30 with the highest microbial activity, which promoted microbial activity to metabolize the carbon sources of carbohydrates, amino acids, polymers and amines; Second were treatments NPK and CK with second highest microbial activity; Thirdwas treatment BS100 with the third highest microbial activity, and fourth was treatment BS15 with the lowest microbial activity. Combined with the PCA score, the order of functional diversity of edaphon communities was BS45>BS30> NPK > CK >BS100>BS15. Overall, we conclude that the microbial activity and functional diversity of upland red soil was significantly affected by biogas slurry application, and that the chemical fertilizer was not completely replaced by biogas slurry. The soil microbial activity reached the highest level when the total nitrogen input ratio of biogas slurry to chemical fertilizer rangedfrom 30% to 45%, which could be expanded inupland red soil of south China.

    biogas slurry; upland red soil; Biolog-ECO; metabolic characteristic; microbial diversity

    公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(201203050);國家自然基金青年項(xiàng)目(31201690);贛鄱英才555工程;國家科技支撐計(jì)劃課題(2015BAD05B01, 2012BAD05B04)

    2015-03-23;

    日期:2015-12-29

    10.5846/stxb201503230555

    *通訊作者Corresponding author.E-mail: zhoujing@issas.ac.cn

    鄭學(xué)博,樊劍波,崔鍵,徐磊,祝振球,周靜,何園球.沼液還田對(duì)旱地紅壤微生物群落代謝與多樣性的影響.生態(tài)學(xué)報(bào),2016,36(18):5865-5875.

    Zheng X B, Fan J B, Cui J, Xu L, Zhu Z Q, Zhou J, He Y Q.Analysis on metabolic characteristics and functional diversity of soil edaphon communities in upland red soil under biogas slurry application.Acta Ecologica Sinica,2016,36(18):5865-5875.

    猜你喜歡
    田量沼液碳源
    不同秸稈還田量對(duì)設(shè)施番茄生長的影響
    緩釋碳源促進(jìn)生物反硝化脫氮技術(shù)研究進(jìn)展
    不同碳源對(duì)銅溜槽用鋁碳質(zhì)涂抹料性能的影響
    昆鋼科技(2021年6期)2021-03-09 06:10:20
    不同農(nóng)作物秸稈直接還田及其還田量對(duì)大豆生長的影響
    秸稈還田下灌水量對(duì)土壤水分運(yùn)移特性研究
    沼液能否預(yù)防病蟲
    大棚絲瓜沼液高產(chǎn)栽培技術(shù)
    上海蔬菜(2016年5期)2016-02-28 13:18:05
    四甘醇作碳源合成Li3V2(PO4)3正極材料及其電化學(xué)性能
    沼液喂豬有講究
    沼液在農(nóng)業(yè)生產(chǎn)中的綜合利用
    精品久久久久久久久av| 亚洲精品456在线播放app| 在线观看人妻少妇| 日本免费在线观看一区| 免费播放大片免费观看视频在线观看| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 中文字幕制服av| 亚洲,欧美,日韩| 边亲边吃奶的免费视频| 色哟哟·www| av不卡在线播放| 尤物成人国产欧美一区二区三区| 女性生殖器流出的白浆| 久久国产乱子免费精品| 欧美日韩视频高清一区二区三区二| 在现免费观看毛片| 99久久人妻综合| 欧美人与善性xxx| 婷婷色麻豆天堂久久| 少妇的逼好多水| 一级毛片我不卡| 一级av片app| 精品一品国产午夜福利视频| 人体艺术视频欧美日本| 亚洲av.av天堂| av卡一久久| 热re99久久精品国产66热6| 另类亚洲欧美激情| 男女边摸边吃奶| av又黄又爽大尺度在线免费看| 97热精品久久久久久| 国产一级毛片在线| 麻豆成人av视频| 三级国产精品欧美在线观看| 夫妻午夜视频| 乱系列少妇在线播放| 久久久久久久亚洲中文字幕| 全区人妻精品视频| 一本色道久久久久久精品综合| 欧美xxxx黑人xx丫x性爽| 欧美成人午夜免费资源| 大香蕉久久网| 久久久久久久久大av| 男男h啪啪无遮挡| 午夜精品国产一区二区电影| 国产在线免费精品| 美女国产视频在线观看| 日本欧美国产在线视频| 日韩视频在线欧美| 亚洲精品亚洲一区二区| 最近中文字幕高清免费大全6| 成人二区视频| 看免费成人av毛片| 精华霜和精华液先用哪个| 久久精品久久精品一区二区三区| 精品99又大又爽又粗少妇毛片| 久久鲁丝午夜福利片| 精华霜和精华液先用哪个| av在线播放精品| 国产一级毛片在线| 黄色配什么色好看| 一本色道久久久久久精品综合| 99久久精品国产国产毛片| 涩涩av久久男人的天堂| 黑丝袜美女国产一区| 精品久久久久久电影网| 亚洲精品日本国产第一区| 一级毛片aaaaaa免费看小| 午夜福利网站1000一区二区三区| 国产精品爽爽va在线观看网站| 精品人妻熟女av久视频| 国产在线男女| 观看av在线不卡| 高清不卡的av网站| 亚洲精品乱码久久久v下载方式| 一个人看的www免费观看视频| 日本午夜av视频| 亚洲色图综合在线观看| 精品亚洲乱码少妇综合久久| 亚洲欧美精品自产自拍| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| av视频免费观看在线观看| 午夜福利网站1000一区二区三区| 欧美一区二区亚洲| 中文字幕亚洲精品专区| 国产一区二区三区综合在线观看 | 国产免费一级a男人的天堂| 国产伦精品一区二区三区视频9| 久久国产亚洲av麻豆专区| 亚洲精品第二区| 亚洲性久久影院| 蜜臀久久99精品久久宅男| 又粗又硬又长又爽又黄的视频| 国产片特级美女逼逼视频| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 精品午夜福利在线看| 久久久久网色| 欧美激情国产日韩精品一区| 最新中文字幕久久久久| 能在线免费看毛片的网站| 亚洲国产毛片av蜜桃av| 精品熟女少妇av免费看| 国产亚洲午夜精品一区二区久久| 亚洲成人一二三区av| 伦理电影大哥的女人| 美女国产视频在线观看| 国产精品久久久久久精品电影小说 | 91狼人影院| 亚洲美女黄色视频免费看| av国产精品久久久久影院| 亚洲人成网站在线播| 国产男女超爽视频在线观看| 亚洲精品色激情综合| 国产淫片久久久久久久久| 国产成人午夜福利电影在线观看| 亚洲熟女精品中文字幕| 国产精品久久久久成人av| 一个人免费看片子| 日产精品乱码卡一卡2卡三| 激情五月婷婷亚洲| 内地一区二区视频在线| 日韩欧美 国产精品| 成人影院久久| 国产深夜福利视频在线观看| 中文天堂在线官网| 久久精品久久精品一区二区三区| 一级爰片在线观看| 亚洲精品日韩在线中文字幕| 久久99热这里只频精品6学生| 国产91av在线免费观看| 色网站视频免费| 啦啦啦在线观看免费高清www| 久久久久国产精品人妻一区二区| 99国产精品免费福利视频| 免费播放大片免费观看视频在线观看| 老司机影院成人| 少妇的逼水好多| 成人黄色视频免费在线看| 91久久精品国产一区二区成人| 国产亚洲av片在线观看秒播厂| 成人毛片a级毛片在线播放| 青春草亚洲视频在线观看| 女人久久www免费人成看片| 免费观看a级毛片全部| 国产成人一区二区在线| 欧美精品国产亚洲| 亚洲内射少妇av| 一区二区av电影网| 搡女人真爽免费视频火全软件| 人人妻人人澡人人爽人人夜夜| 丝袜喷水一区| 国产乱人偷精品视频| 国产高清三级在线| 妹子高潮喷水视频| 日本黄色片子视频| 乱系列少妇在线播放| 国产精品一区www在线观看| 国产精品熟女久久久久浪| 一级毛片我不卡| 99久久精品国产国产毛片| a级毛片免费高清观看在线播放| 国产成人91sexporn| 视频区图区小说| 亚洲成色77777| 黄色配什么色好看| 国产男女超爽视频在线观看| 波野结衣二区三区在线| 黄色日韩在线| 亚洲第一区二区三区不卡| 男女边摸边吃奶| 精品一区二区免费观看| 亚洲一级一片aⅴ在线观看| 日韩av免费高清视频| 日本猛色少妇xxxxx猛交久久| 99久久精品热视频| 国产精品国产三级专区第一集| 日韩av在线免费看完整版不卡| 蜜臀久久99精品久久宅男| 国产精品.久久久| 国产在线视频一区二区| 国产成人精品福利久久| 国内少妇人妻偷人精品xxx网站| 汤姆久久久久久久影院中文字幕| 看非洲黑人一级黄片| 国产精品免费大片| 欧美一区二区亚洲| 天堂8中文在线网| 亚洲国产精品一区三区| 欧美zozozo另类| 国产成人a∨麻豆精品| 国产免费一区二区三区四区乱码| 婷婷色综合www| 免费观看的影片在线观看| 女性被躁到高潮视频| 狂野欧美激情性bbbbbb| 少妇 在线观看| 成人二区视频| 亚洲图色成人| 少妇精品久久久久久久| 亚洲精品一二三| 秋霞伦理黄片| 亚洲性久久影院| 亚洲最大成人中文| 欧美日韩一区二区视频在线观看视频在线| 国产熟女欧美一区二区| 精品国产一区二区三区久久久樱花 | 国产亚洲欧美精品永久| 国产精品福利在线免费观看| 99热这里只有是精品50| 水蜜桃什么品种好| 久久人人爽人人片av| 亚洲成人一二三区av| 国产美女午夜福利| 美女视频免费永久观看网站| 嘟嘟电影网在线观看| 色视频www国产| 国产在视频线精品| 性色avwww在线观看| 简卡轻食公司| 妹子高潮喷水视频| 精品99又大又爽又粗少妇毛片| 久久精品国产鲁丝片午夜精品| h视频一区二区三区| 日本黄色日本黄色录像| 国产黄片美女视频| 91aial.com中文字幕在线观看| 日韩中文字幕视频在线看片 | 久久国产乱子免费精品| 这个男人来自地球电影免费观看 | av.在线天堂| 国产精品偷伦视频观看了| 丰满人妻一区二区三区视频av| 卡戴珊不雅视频在线播放| 久久青草综合色| 日本欧美国产在线视频| 最近的中文字幕免费完整| 交换朋友夫妻互换小说| 纵有疾风起免费观看全集完整版| 大码成人一级视频| 欧美性感艳星| av卡一久久| 下体分泌物呈黄色| 美女视频免费永久观看网站| 久久久a久久爽久久v久久| 成人午夜精彩视频在线观看| 国产一区二区三区综合在线观看 | 一个人免费看片子| 夜夜爽夜夜爽视频| 蜜桃在线观看..| 欧美精品一区二区大全| 黄片无遮挡物在线观看| 亚洲怡红院男人天堂| 亚洲欧洲日产国产| 国产成人免费无遮挡视频| 亚洲欧美日韩另类电影网站 | 亚洲精华国产精华液的使用体验| 人人妻人人爽人人添夜夜欢视频 | 观看免费一级毛片| 日韩视频在线欧美| 三级经典国产精品| 国产高潮美女av| 亚洲成人中文字幕在线播放| 日韩中文字幕视频在线看片 | 欧美高清成人免费视频www| 亚洲在久久综合| 高清黄色对白视频在线免费看 | 久久久午夜欧美精品| 黄片wwwwww| 国国产精品蜜臀av免费| 欧美一区二区亚洲| 大香蕉97超碰在线| 又大又黄又爽视频免费| 免费大片18禁| 高清在线视频一区二区三区| 观看av在线不卡| 亚洲精品日本国产第一区| 亚洲成色77777| 亚洲av电影在线观看一区二区三区| 观看免费一级毛片| 成人黄色视频免费在线看| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 国产精品欧美亚洲77777| av在线app专区| 国产亚洲5aaaaa淫片| 在线免费十八禁| 久久久久久久久久人人人人人人| 嫩草影院入口| 黄色日韩在线| 欧美成人a在线观看| 熟女av电影| 精品久久久久久久久av| av在线app专区| 欧美日韩精品成人综合77777| 狂野欧美白嫩少妇大欣赏| 黄色日韩在线| 国产高清有码在线观看视频| kizo精华| 七月丁香在线播放| 秋霞在线观看毛片| 日本黄大片高清| 观看美女的网站| 国产久久久一区二区三区| av网站免费在线观看视频| 亚洲综合精品二区| 黄色日韩在线| 99久久综合免费| 国产伦在线观看视频一区| 高清av免费在线| 97在线人人人人妻| 日本vs欧美在线观看视频 | 欧美3d第一页| 一区二区三区精品91| 成年人午夜在线观看视频| 99视频精品全部免费 在线| 青青草视频在线视频观看| 国产高潮美女av| 一本久久精品| 熟女电影av网| 国模一区二区三区四区视频| 女性生殖器流出的白浆| 啦啦啦视频在线资源免费观看| 国产午夜精品一二区理论片| 色哟哟·www| 国产成人精品一,二区| 男女下面进入的视频免费午夜| 久久久久性生活片| 91午夜精品亚洲一区二区三区| 男人和女人高潮做爰伦理| 精品国产一区二区三区久久久樱花 | av不卡在线播放| 精品酒店卫生间| 国产精品女同一区二区软件| 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| 蜜桃在线观看..| 国国产精品蜜臀av免费| 80岁老熟妇乱子伦牲交| a级毛片免费高清观看在线播放| 欧美97在线视频| 亚洲av.av天堂| 国产深夜福利视频在线观看| 久久久成人免费电影| 亚洲精品国产av蜜桃| 日本与韩国留学比较| 在线免费观看不下载黄p国产| 久久人妻熟女aⅴ| 久久ye,这里只有精品| 男女无遮挡免费网站观看| 色吧在线观看| 乱码一卡2卡4卡精品| 视频区图区小说| 久久久久久久大尺度免费视频| 在线观看免费日韩欧美大片 | 精品久久久精品久久久| 国产有黄有色有爽视频| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩一区二区| av国产精品久久久久影院| 精品少妇黑人巨大在线播放| 天堂8中文在线网| 人妻一区二区av| 美女主播在线视频| 少妇裸体淫交视频免费看高清| 日韩欧美精品免费久久| 99热全是精品| 精品视频人人做人人爽| 亚洲欧美日韩东京热| 全区人妻精品视频| 少妇人妻 视频| 国产成人a∨麻豆精品| 熟妇人妻不卡中文字幕| 美女脱内裤让男人舔精品视频| 男女国产视频网站| 熟女人妻精品中文字幕| 亚洲丝袜综合中文字幕| 精品人妻偷拍中文字幕| 麻豆国产97在线/欧美| a级毛色黄片| 丝瓜视频免费看黄片| 国产在线免费精品| 成人国产av品久久久| 午夜福利高清视频| 狂野欧美白嫩少妇大欣赏| 中文字幕av成人在线电影| 色婷婷久久久亚洲欧美| 亚洲精品aⅴ在线观看| 美女cb高潮喷水在线观看| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频| 成人免费观看视频高清| 97热精品久久久久久| 国产av精品麻豆| 一本久久精品| 永久免费av网站大全| 五月天丁香电影| 亚洲精品aⅴ在线观看| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 99热6这里只有精品| 国产亚洲5aaaaa淫片| 欧美丝袜亚洲另类| 日韩中字成人| 九九爱精品视频在线观看| 国产伦精品一区二区三区视频9| 欧美极品一区二区三区四区| 国产成人91sexporn| 国产 一区 欧美 日韩| 亚洲国产精品专区欧美| 精品久久久噜噜| av免费在线看不卡| 日韩电影二区| 久久影院123| 在线精品无人区一区二区三 | 黑丝袜美女国产一区| 好男人视频免费观看在线| 精品国产露脸久久av麻豆| 美女cb高潮喷水在线观看| 身体一侧抽搐| 亚洲国产精品成人久久小说| 22中文网久久字幕| 人妻系列 视频| 国产真实伦视频高清在线观看| 亚洲av不卡在线观看| 亚洲精品日韩av片在线观看| av线在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 免费黄频网站在线观看国产| 免费不卡的大黄色大毛片视频在线观看| tube8黄色片| 精品午夜福利在线看| 国产一区有黄有色的免费视频| 高清av免费在线| 国产精品国产av在线观看| 舔av片在线| 国产成人一区二区在线| 一区二区三区四区激情视频| 性高湖久久久久久久久免费观看| 中文字幕精品免费在线观看视频 | 99热这里只有是精品在线观看| 街头女战士在线观看网站| 晚上一个人看的免费电影| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久| 麻豆国产97在线/欧美| 欧美日本视频| 在线观看美女被高潮喷水网站| 一个人看视频在线观看www免费| 卡戴珊不雅视频在线播放| 国产有黄有色有爽视频| 亚洲欧美一区二区三区国产| 日本色播在线视频| 国产在线一区二区三区精| 99久久中文字幕三级久久日本| 免费少妇av软件| 色哟哟·www| 91狼人影院| 国产av国产精品国产| 国产精品一区二区在线观看99| 国产亚洲精品久久久com| 在线观看av片永久免费下载| 蜜桃亚洲精品一区二区三区| 精品久久久久久电影网| 一级毛片aaaaaa免费看小| 亚洲国产精品一区三区| 人人妻人人澡人人爽人人夜夜| 亚洲四区av| 日韩人妻高清精品专区| 亚洲精品中文字幕在线视频 | 国产亚洲午夜精品一区二区久久| www.色视频.com| av在线老鸭窝| 日韩强制内射视频| 国产精品.久久久| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 美女高潮的动态| 搡女人真爽免费视频火全软件| 欧美一级a爱片免费观看看| 成人影院久久| 黄色一级大片看看| 久久久色成人| 欧美日韩在线观看h| 伊人久久国产一区二区| 国国产精品蜜臀av免费| 欧美激情极品国产一区二区三区 | 免费高清在线观看视频在线观看| 熟妇人妻不卡中文字幕| 大片电影免费在线观看免费| 日日摸夜夜添夜夜添av毛片| 国产成人a∨麻豆精品| 在线看a的网站| 一级黄片播放器| 欧美日韩精品成人综合77777| 久久久a久久爽久久v久久| 成人特级av手机在线观看| 成人高潮视频无遮挡免费网站| 日本一二三区视频观看| 只有这里有精品99| 精品亚洲乱码少妇综合久久| 一边亲一边摸免费视频| 久久久久久久久久成人| 成人无遮挡网站| 国产高潮美女av| 亚洲精品456在线播放app| 国产精品一及| 小蜜桃在线观看免费完整版高清| 少妇猛男粗大的猛烈进出视频| 人妻夜夜爽99麻豆av| 国产精品一二三区在线看| 又大又黄又爽视频免费| 亚洲丝袜综合中文字幕| 国精品久久久久久国模美| 亚洲欧美一区二区三区黑人 | 亚洲人成网站在线观看播放| 国精品久久久久久国模美| 高清午夜精品一区二区三区| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 大又大粗又爽又黄少妇毛片口| 久久综合国产亚洲精品| 国产av精品麻豆| 最近2019中文字幕mv第一页| 少妇人妻 视频| 久久久久久久亚洲中文字幕| 少妇裸体淫交视频免费看高清| 在线天堂最新版资源| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| 国产色婷婷99| 日本wwww免费看| 美女高潮的动态| 男女免费视频国产| 大香蕉97超碰在线| 国产真实伦视频高清在线观看| 最近最新中文字幕大全电影3| 青春草亚洲视频在线观看| 黄色日韩在线| 韩国av在线不卡| 久久久久久伊人网av| 国产午夜精品一二区理论片| 久久午夜福利片| 国产精品一区二区性色av| 美女视频免费永久观看网站| 国产亚洲91精品色在线| 不卡视频在线观看欧美| 亚洲四区av| 日本午夜av视频| 久久精品国产a三级三级三级| 国产精品一及| 亚洲丝袜综合中文字幕| 午夜福利视频精品| 国产永久视频网站| 日本黄色片子视频| 26uuu在线亚洲综合色| 久久久久久久亚洲中文字幕| 久久青草综合色| 身体一侧抽搐| 亚洲欧洲日产国产| 看十八女毛片水多多多| 午夜激情福利司机影院| 七月丁香在线播放| 亚洲在久久综合| 麻豆成人午夜福利视频| 26uuu在线亚洲综合色| av视频免费观看在线观看| 欧美成人a在线观看| 国产极品天堂在线| 精品亚洲成a人片在线观看 | 亚洲成色77777| 青春草国产在线视频| 成人国产av品久久久| 高清在线视频一区二区三区| 99九九线精品视频在线观看视频| 国产成人精品久久久久久| 欧美极品一区二区三区四区| 精品一区在线观看国产| 国产精品久久久久久久电影| 日本黄色片子视频| 久久久国产一区二区| 成人亚洲精品一区在线观看 | 身体一侧抽搐| 欧美高清性xxxxhd video| 青春草亚洲视频在线观看| 亚洲欧美精品自产自拍| 中文字幕免费在线视频6| 97在线视频观看| av福利片在线观看| 久久久久性生活片| 亚洲av福利一区| 国产精品一区二区三区四区免费观看| 午夜激情福利司机影院| 日韩三级伦理在线观看| 大香蕉久久网| 麻豆精品久久久久久蜜桃| 亚洲国产欧美在线一区| www.av在线官网国产| 熟女人妻精品中文字幕| 99久久中文字幕三级久久日本| 亚洲最大成人中文| 亚洲精品第二区| 国产精品人妻久久久影院| 国产欧美亚洲国产| 99热全是精品| 啦啦啦在线观看免费高清www| 九草在线视频观看| 麻豆乱淫一区二区| 高清日韩中文字幕在线| 亚洲无线观看免费| a 毛片基地| 菩萨蛮人人尽说江南好唐韦庄|