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

    秦嶺地區(qū)不同林分土壤微生物群落代謝特征

    2016-07-21 01:23:05曹永昌譚向平和文祥耿增超
    生態(tài)學(xué)報(bào) 2016年10期
    關(guān)鍵詞:土壤微生物秦嶺

    曹永昌,譚向平,和文祥,*,耿增超,劉 帥,佘 雕,侯 林

    1 西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院,農(nóng)業(yè)部西北植物營(yíng)養(yǎng)與農(nóng)業(yè)環(huán)境重點(diǎn)實(shí)驗(yàn)室,楊凌 712100 2 西北農(nóng)林科技大學(xué)水土保持研究所,楊凌 712100 3 西北農(nóng)林科技大學(xué)林學(xué)院,楊凌 712100

    ?

    秦嶺地區(qū)不同林分土壤微生物群落代謝特征

    曹永昌1,譚向平1,和文祥1,*,耿增超1,劉帥1,佘雕2,侯林3

    1 西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院,農(nóng)業(yè)部西北植物營(yíng)養(yǎng)與農(nóng)業(yè)環(huán)境重點(diǎn)實(shí)驗(yàn)室,楊凌712100 2 西北農(nóng)林科技大學(xué)水土保持研究所,楊凌712100 3 西北農(nóng)林科技大學(xué)林學(xué)院,楊凌712100

    摘要:利用BIOLOG微孔板法研究了秦嶺山脈的銳齒櫟(Quercus aliena var. acutidentata)、油松(Pinus tabuliformis)、華山松(Pinus armandii)、松櫟混交、云杉(Picea asperata)6個(gè)林地5種典型林分土壤微生物群落代謝特征。發(fā)現(xiàn):(1)5種典型林分中,華山松的平均顏色變化率(AWCD)最高,然后是銳齒櫟1>云杉>油松>銳齒櫟2>松櫟混交。位于不同林場(chǎng)的銳齒櫟林地土壤AWCD值差異較大;(2)土壤微生物功能多樣性指數(shù)與平均顏色變化率表現(xiàn)一致,6種林分之間差異顯著,且不同的林分土壤微生物對(duì)6種碳源利用差異較大。(3)主成分分析顯示各土壤微生物功能多樣性具有顯著差異,其綜合因子得分為華山松>銳齒櫟1>油松>云杉>銳齒櫟2>松櫟混交。(4)冗余分析表明土壤pH、有機(jī)質(zhì)、全氮、堿解氮、速效磷和速效鉀的綜合作用對(duì)土壤微生物群落功能多樣性有顯著影響,其中速效磷和pH與土壤微生物功能多樣性最密切。

    關(guān)鍵詞:秦嶺;土壤微生物;代謝特征;BIOLOG

    土壤微生物作為凋落物的主要分解者,是森林土壤中最具活力的部分。土壤微生物的分布和活性是森林生產(chǎn)實(shí)踐綜合評(píng)價(jià)的主要依據(jù)之一[1],其特征受到水熱條件、季節(jié)動(dòng)態(tài)以及海拔、特別是樹(shù)種組成的影響[2- 4]。土壤微生物對(duì)于植被以及環(huán)境因子的響應(yīng)機(jī)制因生態(tài)系統(tǒng)而異[5]。即便在相同的氣候條件和土壤類型下,不同植被下的土壤微生物仍然存在較大差異[3],因?yàn)樯种械闹饕獦?shù)種可以通過(guò)凋落物的數(shù)量和質(zhì)量以及根系分泌物來(lái)影響或控制土壤微生物群落的結(jié)構(gòu)和功能[6]。對(duì)森林土壤微生物代謝特征的探討,有助于深入研究森林土壤有機(jī)質(zhì)轉(zhuǎn)化的生物化學(xué)機(jī)制。

    秦嶺是我國(guó)南北的氣候分界線,南北氣候差異大,北坡和南坡分屬暖溫帶半濕潤(rùn)氣候帶和亞熱帶濕潤(rùn)氣候帶。其獨(dú)特的地理位置和復(fù)雜的自然環(huán)境特點(diǎn)孕育了豐富的天然植物資源,在不同的海拔高度上生長(zhǎng)著一些典型的植物,如海拔1400—1800m的銳齒櫟,1400—2300m的華山松,1500—2400m的油松,1800m以上的云杉和松櫟混交[7]。作為我國(guó)地理的代表性山區(qū),但是對(duì)秦嶺地區(qū)林地土壤微生物群落代謝特征研究較少,僅限于對(duì)微生物數(shù)量的描述[8]。而且在過(guò)去的幾十年里,該地區(qū)森林公園的開(kāi)發(fā)以及土地利用和全球氣候的不斷變化,直接或間接影響森林土壤營(yíng)養(yǎng)元素的周轉(zhuǎn),亟需探明該地區(qū)不同典型林分下土壤環(huán)境和微生物群落活性的特征及其間關(guān)系。

    BIOLOG微孔板法是一種以微孔板碳源利用為基礎(chǔ)的定量分析方法,可更簡(jiǎn)單、快速的描述微生物群落代謝特征[9],并廣泛應(yīng)用于評(píng)價(jià)土壤微生物群落的功能多樣性[10- 11]。為此本文擬以秦嶺地區(qū)五種典型林分為對(duì)象,采用BIOLOG微孔板法研究不同林分下土壤微生物代謝特征,以期探明:(1)揭示不同天然林分土壤微生物主要碳源類型以及代謝特征;(2)分析土壤微生物代謝特征與環(huán)境因子之間的內(nèi)在聯(lián)系。最終為加強(qiáng)秦嶺地區(qū)天然林生態(tài)系統(tǒng)碳循環(huán)及土壤肥力的保持提供基礎(chǔ)數(shù)據(jù)和理論依據(jù)。

    1材料與方法

    1.1供試材料

    1.1.1供試土樣

    在秦嶺山脈的馬頭灘林場(chǎng)以及辛家山林場(chǎng)中選取5種典型林分:銳齒櫟(Quercusalienavar.acutidentata)、油松(Pinustabuliformis)、華山松(Pinusarmandii)、松櫟混交及云杉(Piceaasperata)的6個(gè)典型林地(表1)。在每個(gè)林地中分別選取20m×20m的3個(gè)標(biāo)準(zhǔn)樣方,每個(gè)樣方中設(shè)5個(gè)采樣點(diǎn)(共15個(gè)采樣點(diǎn)),采樣時(shí)去除表層的枯枝落葉,S型采取0—10cm的土樣,將每個(gè)樣方中5個(gè)采樣點(diǎn)所采集到的土壤樣品混勻,一部分置于4℃冰箱中保存?zhèn)溆茫瑴y(cè)定土壤微生物多樣性;另部分土壤風(fēng)干,過(guò)篩,常規(guī)方法分析土壤基本理化性質(zhì),結(jié)果見(jiàn)表2。

    1.1.2采用BIOLOG微孔板法測(cè)定細(xì)菌群落代謝特征

    取4℃下保存新鮮土樣,過(guò)篩,稱取相當(dāng)5g烘干土重的新鮮土樣加入內(nèi)有45 mL 0.85%滅菌鹽水的三角瓶中后,加無(wú)菌棉花塞,振蕩30min,轉(zhuǎn)速為200r/min。將得到的微生物懸浮液按逐步稀釋法,依次稀釋為10-2、10-3梯度液。用10-3稀釋液接種ECO生態(tài)測(cè)試板,接種量為150μL,每樣1板,每板3次重復(fù)。將接種好的測(cè)試板加蓋,在(25±1)℃下連續(xù)培養(yǎng)120h,每隔12h用BIOLOG自動(dòng)讀數(shù)裝置在590nm下讀數(shù)。其中的31種單一碳源分為6大類[12]。含氨基酸類(6種)、羧酸類(5種)、糖類(12種)、胺類(2種)、多聚物類(4種)和芳香化合物類(2種)[13]。

    表1 樣地基本信息

    銳齒櫟Quercusalienavar.acutidentata;油松Pinustabuliformis;華山松Pinusarmandii;松櫟混交 mixed forests;云杉Piceaasperata

    表2 供試土樣理化性質(zhì)

    數(shù)值為平均值(3次重復(fù))±標(biāo)準(zhǔn)誤差

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

    1.2.1平均顏色變化率(AWCD)

    式中,Ci為各反應(yīng)孔在590nm下的光密度值;R為ECO板對(duì)照孔的光密度值,Ci-R小于零的孔,計(jì)算中記為零,即:Ci-R≥0。

    1.2.2Gompertz方程曲線[14]。

    式中,A為最大吸光度,λ為到達(dá)指數(shù)增長(zhǎng)期所需時(shí)間(時(shí)滯相),Um為最大吸光度變化率。

    1.2.3群落豐富度指數(shù)(S)

    用碳源代謝孔的數(shù)目(AWCD> 0.2,則代表該孔碳源被利用,該孔即為反應(yīng)孔)表示群落豐富度指數(shù)。

    1.2.4Shannon多樣性指數(shù)(H′)及McIntosh多樣性指數(shù)(U)

    式中,Pi為每一孔吸光度除以所有孔吸光度的和。

    Ni=Ci-R

    1.2.5Gini多樣性指數(shù)(D)

    1.3數(shù)據(jù)統(tǒng)計(jì)

    試驗(yàn)數(shù)據(jù)采用Microsoft Excel 2003進(jìn)行數(shù)據(jù)處理和制圖,應(yīng)用SPSS 12.0進(jìn)行Gompertz方程擬合、單因素方差(One-way ANOVA)分析及主成分分析,分析不同樣地間各指標(biāo)的差異顯著性以及對(duì)碳源利用模式的差異,采用CANOCO 4.5進(jìn)行冗余分析(RDA)和制圖,分析環(huán)境因子對(duì)土壤微生物群落碳源利用模式的影響。

    2結(jié)果與分析

    2.1ELISA反應(yīng)顏色變化率(AWCD)

    圖1 六個(gè)處理平均顏色變化率Fig.1 AWCD of six treatments

    平均顏色變化率(AWCD)表征微生物群落碳源利用率,是土壤微生物群落利用單一碳源能力的一個(gè)重要指標(biāo),反映了土壤微生物活性、微生物群落生理功能多樣性。供試土壤AWCD值如圖1顯示,隨培養(yǎng)時(shí)間的延長(zhǎng)AWCD增大,在24 h之前各處理的AWCD很小,為0.004—0.018,說(shuō)明在24h之內(nèi)碳源基本未被利用;培養(yǎng)24h后AWCD快速升高,反映出此后碳源被大幅度利用。在相同培養(yǎng)時(shí)間內(nèi),AWCD值為華山松>銳齒櫟1>云杉>油松>松櫟混交>銳齒櫟2。

    將AWCD值與時(shí)間按Gompertz方程曲線進(jìn)行擬合,結(jié)果決定系數(shù)R2均大于0.99(表3),可見(jiàn)最大吸光度以及最大吸光度變化率均是華山松最大,其次是銳齒櫟1和云杉,而松櫟混交和銳齒櫟2最低,這表明華山松對(duì)31種碳源的利用最快,增長(zhǎng)斜率最大。松櫟混交和銳齒櫟2對(duì)31碳源的利用最慢,增長(zhǎng)斜率最小。時(shí)滯相即達(dá)到指數(shù)增長(zhǎng)期所需要的時(shí)間表現(xiàn)為華山松<云杉<銳齒櫟1<銳齒櫟2<松櫟混交<油松,這說(shuō)明華山松林地土壤微生物能很快的利用碳源,而松櫟混交和油松對(duì)碳源的利用較慢。

    表3 Gompertz方程曲線擬合參數(shù)

    2.2土壤微生物群落多樣性指數(shù)

    土壤微生物群落功能多樣性是土壤微生物群落狀態(tài)與功能的指標(biāo),反映土壤中微生物的生態(tài)特征。從表4可看出,BIOLOG微平板培養(yǎng)72h時(shí),華山松的AWCD最高,其次是銳齒櫟1,松櫟混交和銳齒櫟2最低;除松櫟混交與銳齒櫟2外,其它林區(qū)之間均差異顯著(P<0.05),這說(shuō)明華山松林地土壤微生物對(duì)碳源消耗量最多,微生物活性最高,銳齒櫟1、云杉、油松次之,松櫟混交和銳齒櫟2土壤微生物對(duì)碳源消耗最少,微生物活性較小。且不同處理的微生物群落Shannon指數(shù)(H)、豐富度指數(shù)(S)、McIntosh指數(shù)(U)、Gini指數(shù)(D)的變化表現(xiàn)出與AWCD類似的規(guī)律,即華山松>銳齒櫟1>云杉>油松,松櫟混交和銳齒櫟2最低。表明華山松的土壤微生物群落豐富度和功能多樣性最高,其次是銳齒櫟1、云杉、油松,銳齒櫟2和松櫟混交最低。并且從表1可以看出同一林分不同地理位置對(duì)土壤群落豐富度和功能多樣性影響較大,銳齒櫟1土壤群落豐富度和功能多樣性要遠(yuǎn)大于位于銳齒櫟2,可能是土壤性質(zhì)(如有機(jī)質(zhì)相差41%)差異較大所致(表2)。

    表4 土壤微生物群落多樣性指數(shù)分析

    數(shù)值為平均值(3次重復(fù))±標(biāo)準(zhǔn)誤差;同列不同字母表示差異顯著(P<0.05)

    2.3土壤微生物對(duì)碳源的利用

    土壤微生物對(duì)不同碳源的利用反映了微生物的代謝功能類群,研究土壤微生物對(duì)不同碳源利用能力的差異,可深入了解不同林分下土壤微生物群落代謝特征。表5看出,不同林分下土壤微生物對(duì)六類碳源的利用差異顯著,其中銳齒櫟1、華山松和云杉林地土壤微生物對(duì)碳源的利用主要集中在糖類和氨基酸類,銳齒櫟1和華山松對(duì)羧酸類的利用程度較低,而云杉對(duì)芳香類的利用程度較低。油松和銳齒櫟2對(duì)碳源的利用主要集中在氨基酸類和多聚物類,對(duì)芳香類的利用程度較低。松櫟混交對(duì)碳源的利用主要集中在胺類和多聚物類,對(duì)羧酸類的利用程度較低??梢?jiàn),不同林分土壤微生物碳代謝群落結(jié)構(gòu)有所不同。

    2.4微生物群落功能的主成分分析

    主成分分析揭示了不同林地中土壤微生物群落對(duì)碳源的利用模式的差別,用培養(yǎng)72h的吸光度數(shù)值進(jìn)行主成分分析(PCA),結(jié)果(表6)表明前3個(gè)主成分方差貢獻(xiàn)率分別為45%、22%和16%,累積方差貢獻(xiàn)率達(dá)到84%。從中提取可以聚集單一碳源變量的數(shù)據(jù)變異(累積方差貢獻(xiàn)率)為67%的前2個(gè)主成分PC1、PC2(特征根值為14.06和6.74)來(lái)分析細(xì)菌群落功能多樣性。由圖2可見(jiàn),6個(gè)處理在主成分坐標(biāo)體系中分布差異十分明顯。PC1的方差貢獻(xiàn)率最大,主要綜合了華山松、銳齒櫟1和松櫟混交的變異。PC2主要集中了云杉的變異,PC3則綜合了油松和銳齒櫟2的變異。華山松處于PC1的正端,典型變量值達(dá)6.85,油松處于PC2的最正端。PC1典型變量值差異極顯著,這種差異表現(xiàn)在: 除了松櫟混交與銳齒櫟2之間差異不顯著外,其它處理相互之間差異均顯著。PC2典型變量值差異也達(dá)到極顯著水平,表現(xiàn)在云杉與其它林分差異顯著,銳齒櫟1、松櫟混交與其它林分差異顯著,但是銳齒櫟1和松櫟混交兩者之間差異不顯著,油松、華山松和云杉之間差異不顯著。六種林區(qū)土壤典型變量值的變異(離散)都很小(圖2)。

    表5 土壤微生物培養(yǎng)72h時(shí)6類碳源利用的平均吸光度

    數(shù)值為平均值(3次重復(fù))±標(biāo)準(zhǔn)誤差;同列不同字母表示差異顯著(P<0.05)

    表6 主成分特征根

    圖2 不同林木土壤細(xì)菌群落主成分分析 Fig.2 Principal component analysis of soil bacteria communities for different tree species

    對(duì)主成分得分大于0.25的碳源進(jìn)一步分析,PC1軸上優(yōu)勢(shì)利用碳源是D-纖維二糖(糖類)和D-蘋果酸(羧酸類);PC2軸上的優(yōu)勢(shì)碳源是L-精氨酸(氨基酸類)、肝糖(多聚物類)和苯乙胺(胺類);PC3軸上的優(yōu)勢(shì)利用碳源是D-木糖/戊醛糖(糖類)和肝糖(多聚物類)。

    表7顯示6種林分主成分分析綜合得分華山松最高,其次是銳齒櫟1,而松櫟混交與銳齒櫟2綜合得分最低。

    2.5冗余分析

    通過(guò)冗余分析常可用來(lái)確定樣地中環(huán)境因子對(duì)土壤微生物群落碳源利用模式的影響。結(jié)果(圖3)表明,土壤pH、有機(jī)質(zhì)、全氮、堿解氮、速效磷和速效鉀的綜合作用對(duì)土壤微生物群落功能多樣性的影響很顯著。根據(jù)箭頭的長(zhǎng)度,有機(jī)質(zhì)、全氮、堿解氮、速效磷和速效鉀對(duì)微生物碳源利用的影響都很大。其中速效磷和pH幾乎與所有碳源的箭頭呈銳角,說(shuō)明它們之間呈正相關(guān),有機(jī)質(zhì)、全氮、堿解氮以及速效鉀與大多數(shù)碳源的箭頭呈銳角,因此呈正相關(guān)。

    表7 6個(gè)處理主成分綜合得分

    圖3 31個(gè)單一碳源與環(huán)境因子的冗余分析Fig.3 Redundancy analysis of 31 sole-carbon sources and environmental variables

    3討論

    在本研究中,位于馬頭灘林場(chǎng)4種林分中華山松(針葉林)土壤微生物群落代謝特征指標(biāo)大于銳齒櫟1(闊葉林),但通常認(rèn)為闊葉林土壤微生物功能多樣性要高于針葉林土壤微生物功能多樣性。這一現(xiàn)象的產(chǎn)生可能是由于林分結(jié)構(gòu)的極大差異所造成的,在6個(gè)采樣地中,華山松海拔較低,由于采樣點(diǎn)靠近河流,造成華山松林地濕度較大,其結(jié)果是華山松林下表層土含水量(29.2%)相對(duì)于其它樣地較高;其次華山松林地的灌木(0.528株/m2)以及草本科植物(0.960株/m2)遠(yuǎn)遠(yuǎn)高于銳齒櫟1林地的灌木(0.383株/m2)以及草本科植物(0.258株/m2)。Waid等人認(rèn)為植被的類型、數(shù)量和化學(xué)組成可能是土壤生物多樣性變化的主要推進(jìn)力量[15]。且華山松林地密度(0.151 株/m2)要遠(yuǎn)遠(yuǎn)高于銳齒櫟1(0.065 株/m2),較高的林地密度導(dǎo)致華山松凋落物的產(chǎn)生量(6.068t/hm2)要高于銳齒櫟1(2.716t/hm2),凋落物作為土壤有機(jī)質(zhì)輸入的主要來(lái)源,其質(zhì)量和數(shù)量必然會(huì)對(duì)土壤微生物群落組成產(chǎn)生影響[16- 17]。且華山松根系在0—20cm的土層中分布較多[18],研究表明林木根系對(duì)其周圍的土壤微生物的數(shù)量和活性具有重要的影響[19]。多方面的林分結(jié)構(gòu)差異的綜合作用可能是造成了華山松土壤微生物功能多樣性高于銳齒櫟1的部分原因。辛家山林場(chǎng)中云杉土壤微生物群落代謝特征指標(biāo)大于銳齒櫟2,且不同的地理位置對(duì)同一林分土壤微生物功能多樣性影響較大,在林分結(jié)構(gòu)差異不大的情況下,養(yǎng)分組成是影響土壤微生物群落結(jié)構(gòu)的主要因素[17]。銳齒櫟2的土壤類型雖然為棕壤,但是由于河流沖積導(dǎo)致土壤含沙量較高,與云杉林和銳齒櫟1相比,土壤養(yǎng)分含量相對(duì)低,養(yǎng)分組成差異較大,微生物種群數(shù)量和多樣性與土壤養(yǎng)分的含量密切相關(guān)[20]。并且植被往往通過(guò)影響土壤基本理化性質(zhì)影響土壤生物區(qū)系[21]。一般而言,對(duì)于同一土壤而言,養(yǎng)分含量越高,土壤中微生物的數(shù)量和多樣性也就越高。所以,不同的土壤類型以及土壤養(yǎng)分含量的差異造成了林分土壤微生物功能多樣性的巨大差異。

    這5種典型林分下土壤微生物的物種豐富度和優(yōu)勢(shì)種群均差異顯著,總體的功能多樣性均有顯著性的差異。這可能是因?yàn)榱址忠约暗乇碇参锏牟煌斐傻模窒聵?shù)種的改變對(duì)地下生態(tài)系統(tǒng)微生物功能多樣性有明顯影響[22]。以往的研究表明,微生物群落結(jié)構(gòu)和功能的差異與來(lái)源與不同優(yōu)勢(shì)樹(shù)種的枯枝落葉的量和生物化學(xué)組成有關(guān),同時(shí)與根系分泌物的關(guān)系也十分密切[23],根系分泌物對(duì)土壤微生物的分布有重要作用[24]。通過(guò)對(duì)銳齒櫟、油松、華山松、云杉和松櫟混交5種林地土壤微生物對(duì)不同類型碳源的相對(duì)利用程度的分析可以看到不同林地對(duì)相同碳源的利用程度上存在明顯的差異。通過(guò)微生物對(duì)不同碳源利用的差異分析表明,不同林分土壤微生物的主要功能代謝類群有很大不同,進(jìn)而可以深入發(fā)現(xiàn)不同林分土壤微生物的群落代謝特征存在顯著差異。胡嬋娟等人發(fā)現(xiàn)不同人工林對(duì)相同碳源的利用上存在明顯的差異,土壤微生物利用率最高的碳源分別為糖類、氨基酸類和羧酸類物質(zhì)[25]。PCA分析結(jié)果發(fā)現(xiàn)各林地土壤微生物對(duì)碳源的利用是有顯著差異的,用RDA對(duì)可能造成差異的因子做進(jìn)一步的分析,發(fā)現(xiàn)土壤理化性質(zhì)(土壤有機(jī)質(zhì)、pH、全氮、堿解氮、速效磷、速效鉀)與土壤微生物功能多樣性關(guān)系密切。這些環(huán)境因子的差異所導(dǎo)致的空間異質(zhì)性則與土壤微生物群落有直接的關(guān)系。研究不同林分對(duì)土壤功能多樣性的影響可以從側(cè)面反映土壤基本理化性質(zhì)對(duì)土壤微生物功能多樣性的綜合作用。

    4結(jié)論

    通過(guò)采樣分析,對(duì)秦嶺地區(qū)5種典型林分的土壤微生物功能多樣性可得出以下結(jié)論:

    (1)不同林分對(duì)土壤微生物功能多樣性以及微生物群代謝特征影響顯著,其中華山松土壤微生物功能多樣性最高,并且處于不同地理位置相同林分的土壤微生物功能多樣性差異顯著。

    (2)6種林分中,根據(jù)主成分分析多樣性指數(shù)以及綜合得分等多種分析可看出林區(qū)土壤微生物功能多樣性排名從高到低排序?yàn)椋喝A山松>銳齒櫟1>油松>云杉>松櫟混交>銳齒櫟2。

    參考文獻(xiàn)(References):

    [1]許光輝, 鄭洪元, 張德生, 盧耀波, 李玉坤, 張淑賢, 劉增柱, 吳文芳, 趙振英, 李鳳珍, 劉瑞君. 長(zhǎng)白山北坡自然保護(hù)區(qū)森林土壤微生物生態(tài)分布及其生化特性的研究. 生態(tài)學(xué)報(bào), 1984, 4(3):207- 223.

    [2]Kara ?, Bolat I. Soil microbial biomass C and N changes in relation to forest conversion in the Northwestern Turkey. Land Degradation & Development, 2008, 19(4): 421- 428.

    [3]Wardle D A. A comparative assessment of factors which influence microbial biomass carbon and nitrogen levels in soil. Biological Reviews, 1992, 67(3): 321 - 358.

    [4]Van Gestel M, Merckx R, Vlassak K. Microbial biomass and activity in soils with fluctuating water contents. Geoderma, 1993, 56(1): 617- 626.

    [5]隋心. 典型溫帶森林土壤有機(jī)碳自然積累與微生物維持關(guān)系[D].北京:中國(guó)科學(xué)院大學(xué), 2013.

    [6]Ayres E, Steltzer H, Berg S, Wallenstein M D, Simmons B L, Tree species traits influence soil physical, chemical, and biological properties in high elevation forests. PLoS one, 2009, 4(6): e5964.

    [7]趙永華, 雷瑞德, 何興元, 賈夏. 秦嶺銳齒櫟林種群生態(tài)位特征研究. 應(yīng)用生態(tài)學(xué)報(bào), 2004, 15(6):913- 918.

    [8]崔芳芳, 劉增文, 付剛段, 而軍, 高文俊. 秦嶺山區(qū)幾種典型森林的土壤微生物特征及其對(duì)人為干擾的響應(yīng). 西北林學(xué)院學(xué)報(bào), 2008, 23(2):129- 134.

    [9]Selmants P C, Hart S C, Boyle S I, Stark J M. Red alder (Alnusrubra) alters community-level soil microbial function in conifer forests of the Pacific Northwest, USA. Soil Biology and Biochemistry, 2005, 37(10): 1860- 1868.

    [10]Grayston S J, Campbell C D. Functional biodiversity of microbial communities in the rhizospheres of hybrid larch (Larixeurolepis) and Sitka spruce (Piceasitchensis). Tree Physiology, 1996, 16(11- 12): 1031- 1038.

    [11]Bossio D A, Scow K M. Impact of carbon and flooding on the metabolic diversity of microbial communities in soils[J]. Applied and Environmental Microbiology, 1995, 61(11): 4043- 4050.

    [12]Papatheodorou E M, Efthimiadou E, Stamou G P. Functional diversity of soil bacteria as affected by management practices and phenological stage ofPhaseolusvulgaris. European Journal of Soil Biology, 2008, 44(4): 429- 436.

    [13]Kong W D, Zhu Y G, Fu B J, Marschner P, He J Z. The veterinary antibiotic oxytetracycline and Cu influence functional diversity of the soil microbial community. Environmental Pollution, 2006, 143(1): 129- 137.

    [14]Verschuere L, Fievez V, Van Vooren L, Verstraete W. The contribution of individual populations to the Biolog pattern of model microbial communities. FEMS microbiology ecology, 1997, 24(4): 353- 362.

    [15]Waid J S. Does soil biodiversity depend upon metabiotic activity and influences?. Applied Soil Ecology, 1999, 13(2): 151- 158.

    [16]Chapin III F S, Chapin M C, Matson P A, Vitousek P M. Principles of Terrestrial Ecosystem Ecology. Berlin: Springer, 2011.

    [17]陳法霖, 鄭華, 歐陽(yáng)志云, 張凱, 屠乃美. 土壤微生物群落結(jié)構(gòu)對(duì)凋落物組成變化的響應(yīng). 土壤學(xué)報(bào), 2011, 48(3): 603- 611.

    [18]趙忠, 李鵬, 王乃江. 渭北黃土高原主要造林樹(shù)種根系分布特征的研究. 應(yīng)用生態(tài)學(xué)報(bào), 2000, 11(1): 37- 39.

    [19]楊玉盛, 何宗明, 鄒雙全, 俞新妥. 格氏栲天然林與人工林根際土壤微生物及其生化特性的研究. 生態(tài)學(xué)報(bào), 1998, 18(2): 198- 202.

    [20]傅伯杰, 陳利頂, 王軍, 孟慶華, 趙文武等. 土地利用結(jié)構(gòu)與生態(tài)過(guò)程. 第四紀(jì)研究, 2003, 23(3): 247- 255.

    [21]Zhang X C, Shao M A. Soil nitrogen loss by erosion as affected by vegetation cover and comprehensive managements in Zhifanggou catchment of hilly Loess Plateau. Acta Geographica Sinica: Chinese Edition, 2000, 55(5): 617- 626.

    [22]Augusto L, Ranger J, Binkley D, Rothe A. Impact of several common tree species of European temperate forests on soil fertility. Annals of Forest Science, 2002, 59(3): 233- 253.

    [23]Salam A K, Katayama A, Kimura M. Activities of some soil enzymes in different land use systems after deforestation in hilly areas of West Lampung, South Sumatra, Indonesia. Soil science and plant nutrition, 1998, 44(1): 93- 103.

    [24]Nayyar A, Hamel C, Lafond G, Gossen B D, Hanson K, Germida J. Soil microbial quality associated with yield reduction in continuous-pea. Applied soil ecology, 2009, 43(1): 115- 121.

    [25]胡嬋娟, 傅伯杰, 劉國(guó)華, 靳甜甜, 劉宇. 黃土丘陵溝壑區(qū)典型人工林下土壤微生物功能多樣性. 生態(tài)學(xué)報(bào), 2009, 29(2): 727- 733.

    The metabolism characteristics of microbial community in different forest soil in Qinling Mountains Area

    CAO Yongchang1, TAN Xiangping1, HE Wenxiang1,*, GENG Zengchao1, LIU Shuai1, SHE Diao2, HOU Lin3

    1CollegeofResourcesandEnvironment,NorthwestAgricultureandForestryUniversity.KeyLaboratoryofPlantNutritionandtheAgri-environmentinNorthwestChina,MinistryofAgriculture,Yangling712100,China2InstituteofSoilandWaterConservation,NorthwestAgricultureandForestryUniversity,Yangling712100,China3CollegeofForestry,NorthwestAgricultureandForestryUniversity,Yangling712100,China

    Abstract:In this study, the BIOLOG microporosity plate technology was employed to investigate the metabolic characteristics of microbial community in six woodlands with five typical types of forest soils, which included forest soils planted with tooth Oak(Quercus aliena var. acutidentata),chinese pine(Pinus tabuliformis),armand pine(Pinus armandii), pine and oak mixed, spruce(Picea asperata). The major findings include: (1) The average well color development (AWCD) of armand pine is the highest among the five kinds of typical forests, followed by tooth Oak1> spruce> chinese pine > tooth Oak2>pine and oak mixed, The AWCD value of tooth Oak woodlands located in different forest farms showed large differences among each other; (2) The index of soil microbial function diversity is consistent with the AWCD. The differences among six forest lands showed significant difference,. Significant difference was also observed for the utilization of six carbon source by soil microorganism in different forest components. (3) The principal component analysis showed that the soil microbial functional diversity among different soils were significant different. The sequence of comprehensive factor scores were armand pine > tooth Oak1> Chinese pine > chinese spruce > tooth Oak2> pine and oak mixed. (4) Redundancy analysis shows that the comprehensive effect of soil pH, organic matter, total nitrogen, alkali hydrolyzable nitrogen, the available phosphorus and available potassium can significantly affect soil microbial community functional diversity, especially the available P and pH which were closely related to soil microbial functional diversity.

    Key Words:Qin Ling; soil microorganism,; metabolism characteristics; BIOLOG

    基金項(xiàng)目:國(guó)家林業(yè)公益性行業(yè)科研專項(xiàng)(201304307);西北農(nóng)林科技大學(xué)“基本科研業(yè)務(wù)費(fèi)科研創(chuàng)新重點(diǎn)項(xiàng)目”(ZD2013012)

    收稿日期:2014- 11- 22; 網(wǎng)絡(luò)出版日期:2015- 09- 28

    *通訊作者

    Corresponding author.E-mail: wenxianghe@nwsuaf.edu.cn

    DOI:10.5846/stxb201411222316

    曹永昌,譚向平,和文祥,耿增超,劉帥,佘雕,侯林.秦嶺地區(qū)不同林分土壤微生物群落代謝特征.生態(tài)學(xué)報(bào),2016,36(10):2978- 2986.

    Cao Y C, Tan X P, He W X, Geng Z C, Liu S, She D, Hou L.The metabolism characteristics of microbial community in different forest soil in Qinling Mountains Area.Acta Ecologica Sinica,2016,36(10):2978- 2986.

    猜你喜歡
    土壤微生物秦嶺
    守護(hù)大秦嶺
    暑期秦嶺游
    洞穿秦嶺
    賈平凹:從秦嶺到秦嶺
    好忙好忙的秦嶺
    長(zhǎng)期施用化肥對(duì)土壤微生物多樣性的影響探討
    基于PCR—DGGE的間作油茶林土壤細(xì)菌多樣性
    基于PCR—DGGE的間作油茶林土壤細(xì)菌多樣性
    黑汁緣微生物肥料在番茄上的應(yīng)用
    在秦嶺
    詩(shī)選刊(2015年6期)2015-10-26 09:47:16
    丁香六月欧美| av不卡在线播放| 免费观看a级毛片全部| 一个人免费在线观看的高清视频| 深夜精品福利| 无遮挡黄片免费观看| 日韩免费av在线播放| 国产成人欧美| 亚洲精品在线美女| 免费在线观看视频国产中文字幕亚洲| 9191精品国产免费久久| 精品高清国产在线一区| 黄色毛片三级朝国网站| 亚洲专区中文字幕在线| 精品国产美女av久久久久小说| 亚洲七黄色美女视频| 久99久视频精品免费| 久久久久久久久免费视频了| 亚洲专区国产一区二区| 国产xxxxx性猛交| 日韩熟女老妇一区二区性免费视频| 黄片大片在线免费观看| 两个人看的免费小视频| 国产精品久久久久久精品古装| 国产高清激情床上av| 18禁美女被吸乳视频| 精品无人区乱码1区二区| 免费在线观看日本一区| 国产欧美日韩一区二区三区在线| 99热只有精品国产| 如日韩欧美国产精品一区二区三区| 国产高清videossex| 亚洲国产精品合色在线| 日韩精品免费视频一区二区三区| 免费观看人在逋| 国产精品秋霞免费鲁丝片| 精品国内亚洲2022精品成人 | 美女高潮喷水抽搐中文字幕| 大香蕉久久成人网| 淫妇啪啪啪对白视频| 亚洲欧美日韩另类电影网站| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 亚洲情色 制服丝袜| 久热爱精品视频在线9| 久久精品熟女亚洲av麻豆精品| 少妇裸体淫交视频免费看高清 | 首页视频小说图片口味搜索| 91九色精品人成在线观看| 亚洲专区中文字幕在线| 午夜福利在线免费观看网站| 一夜夜www| 亚洲av熟女| 欧美日韩成人在线一区二区| 最近最新中文字幕大全免费视频| 免费女性裸体啪啪无遮挡网站| 日本黄色日本黄色录像| 岛国毛片在线播放| 亚洲精品美女久久av网站| 国产精品一区二区在线观看99| 久久久国产精品麻豆| 成在线人永久免费视频| 亚洲av第一区精品v没综合| 日韩人妻精品一区2区三区| 中文字幕人妻丝袜一区二区| 一边摸一边抽搐一进一小说 | 久久久国产欧美日韩av| 久久中文看片网| 满18在线观看网站| 天堂动漫精品| 女人久久www免费人成看片| 国产av一区二区精品久久| 精品国产一区二区三区四区第35| 亚洲国产看品久久| 日韩欧美一区视频在线观看| 欧美在线黄色| 黄色女人牲交| 91麻豆精品激情在线观看国产 | 亚洲黑人精品在线| 精品99又大又爽又粗少妇毛片 | 国产亚洲精品综合一区在线观看| 欧美日本亚洲视频在线播放| 色播亚洲综合网| 性色avwww在线观看| 中文字幕av在线有码专区| 天堂网av新在线| 亚洲av成人精品一区久久| 国产高清videossex| 在线播放国产精品三级| 国产在线精品亚洲第一网站| av女优亚洲男人天堂| 黄色片一级片一级黄色片| 中文字幕人妻熟人妻熟丝袜美 | 观看美女的网站| 亚洲一区高清亚洲精品| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 男人的好看免费观看在线视频| 国产色爽女视频免费观看| 国产精品久久久久久久电影 | 在线播放无遮挡| 午夜福利18| 精品免费久久久久久久清纯| 尤物成人国产欧美一区二区三区| 12—13女人毛片做爰片一| 麻豆国产97在线/欧美| 国产成人aa在线观看| 午夜视频国产福利| 亚洲欧美日韩无卡精品| 99riav亚洲国产免费| 国产 一区 欧美 日韩| 成年女人看的毛片在线观看| 国产欧美日韩精品亚洲av| 在线观看日韩欧美| 免费av观看视频| 2021天堂中文幕一二区在线观| 午夜亚洲福利在线播放| 国产久久久一区二区三区| 久久精品国产综合久久久| 看片在线看免费视频| 成年版毛片免费区| 国产一级毛片七仙女欲春2| 精品国产三级普通话版| 天堂√8在线中文| 国产一区二区三区在线臀色熟女| av视频在线观看入口| 高清毛片免费观看视频网站| 国产亚洲精品久久久com| 国产麻豆成人av免费视频| 99热这里只有精品一区| 男女那种视频在线观看| 欧美乱色亚洲激情| 窝窝影院91人妻| 在线播放国产精品三级| 18+在线观看网站| 国产黄片美女视频| 欧美三级亚洲精品| 村上凉子中文字幕在线| 欧美午夜高清在线| 国产一区二区激情短视频| 母亲3免费完整高清在线观看| 日韩欧美国产在线观看| 国产三级中文精品| 一个人观看的视频www高清免费观看| 黄色丝袜av网址大全| 精品国产亚洲在线| 舔av片在线| 欧美3d第一页| 高潮久久久久久久久久久不卡| 国产毛片a区久久久久| 悠悠久久av| 国产中年淑女户外野战色| 亚洲最大成人手机在线| av片东京热男人的天堂| 色综合欧美亚洲国产小说| 三级男女做爰猛烈吃奶摸视频| 亚洲中文字幕一区二区三区有码在线看| 男女做爰动态图高潮gif福利片| 夜夜爽天天搞| 国产精品久久电影中文字幕| 国产午夜精品论理片| 国产男靠女视频免费网站| 久久久久久大精品| 热99re8久久精品国产| 看免费av毛片| 亚洲成人久久性| 国产三级黄色录像| 成人性生交大片免费视频hd| 日韩欧美在线乱码| 欧美av亚洲av综合av国产av| 两个人视频免费观看高清| 99久久综合精品五月天人人| 在线观看舔阴道视频| 波多野结衣高清无吗| 男人舔女人下体高潮全视频| 好男人电影高清在线观看| 国产蜜桃级精品一区二区三区| 丰满乱子伦码专区| or卡值多少钱| 免费看美女性在线毛片视频| 性色avwww在线观看| 日本一二三区视频观看| 久久精品国产亚洲av涩爱 | 九九久久精品国产亚洲av麻豆| 国产黄色小视频在线观看| 国产精品美女特级片免费视频播放器| 日韩欧美在线二视频| 国产午夜精品论理片| 婷婷精品国产亚洲av在线| 日本精品一区二区三区蜜桃| 一夜夜www| 欧美成人a在线观看| 欧美色欧美亚洲另类二区| 看片在线看免费视频| 中文资源天堂在线| 不卡一级毛片| 最新中文字幕久久久久| 亚洲av一区综合| 亚洲国产色片| 人妻丰满熟妇av一区二区三区| 久久欧美精品欧美久久欧美| 午夜福利视频1000在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 成年女人看的毛片在线观看| 两个人视频免费观看高清| 男女那种视频在线观看| 午夜福利欧美成人| 在线免费观看不下载黄p国产 | 久99久视频精品免费| 国产精品 欧美亚洲| 亚洲电影在线观看av| 亚洲熟妇中文字幕五十中出| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 国产亚洲精品综合一区在线观看| 久久久久久久久久黄片| av片东京热男人的天堂| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清| 国产精品永久免费网站| 亚洲自拍偷在线| 日韩欧美精品免费久久 | 日本免费一区二区三区高清不卡| 最近最新免费中文字幕在线| 国产成+人综合+亚洲专区| 精品乱码久久久久久99久播| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 十八禁人妻一区二区| 欧美日本视频| 日本精品一区二区三区蜜桃| 免费看日本二区| 国产精品1区2区在线观看.| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合久久99| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频| 久久精品国产综合久久久| 黄色女人牲交| 久久久国产成人精品二区| 亚洲 欧美 日韩 在线 免费| 老鸭窝网址在线观看| 麻豆一二三区av精品| 亚洲精品日韩av片在线观看 | 国产精品久久电影中文字幕| 亚洲人成网站高清观看| 岛国在线观看网站| avwww免费| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 看片在线看免费视频| 在线观看av片永久免费下载| www.999成人在线观看| 国产亚洲av嫩草精品影院| 亚洲av成人精品一区久久| 精品一区二区三区人妻视频| 国产精品一区二区三区四区久久| 给我免费播放毛片高清在线观看| 国产高清videossex| 国产亚洲精品综合一区在线观看| 国产高清视频在线播放一区| 天堂动漫精品| 日本 欧美在线| 高清在线国产一区| 在线观看免费视频日本深夜| 国产精品久久久人人做人人爽| www.www免费av| av专区在线播放| 国产一区在线观看成人免费| 久久久久性生活片| 欧美另类亚洲清纯唯美| 极品教师在线免费播放| 久久久久国内视频| 桃色一区二区三区在线观看| 国产成人a区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产真实伦视频高清在线观看 | a级一级毛片免费在线观看| 婷婷精品国产亚洲av在线| 精品免费久久久久久久清纯| 色尼玛亚洲综合影院| 国产一区二区三区视频了| 国产精品久久久人人做人人爽| 国产极品精品免费视频能看的| 91麻豆精品激情在线观看国产| 色精品久久人妻99蜜桃| 狂野欧美激情性xxxx| 亚洲国产欧美网| 麻豆久久精品国产亚洲av| 99精品欧美一区二区三区四区| 男插女下体视频免费在线播放| 两个人看的免费小视频| 亚洲电影在线观看av| 十八禁网站免费在线| 在线观看日韩欧美| 中文资源天堂在线| 中文字幕人成人乱码亚洲影| 欧美成人a在线观看| 国产成人a区在线观看| 日韩精品青青久久久久久| 我要搜黄色片| 午夜久久久久精精品| 日韩免费av在线播放| 精品国产三级普通话版| 精品福利观看| 午夜福利高清视频| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 99视频精品全部免费 在线| 老汉色av国产亚洲站长工具| 99在线视频只有这里精品首页| 中文字幕av在线有码专区| 亚洲国产精品999在线| 亚洲精品影视一区二区三区av| 亚洲天堂国产精品一区在线| ponron亚洲| 搡老妇女老女人老熟妇| 香蕉久久夜色| 国产精品女同一区二区软件 | 亚洲精品日韩av片在线观看 | 首页视频小说图片口味搜索| 精品福利观看| 天堂√8在线中文| 无限看片的www在线观看| 男女午夜视频在线观看| 极品教师在线免费播放| 亚洲成a人片在线一区二区| 国产精品99久久久久久久久| 五月伊人婷婷丁香| 俺也久久电影网| 亚洲精品456在线播放app | 中文亚洲av片在线观看爽| 丁香欧美五月| 国产高清有码在线观看视频| 中文在线观看免费www的网站| 人妻夜夜爽99麻豆av| 88av欧美| 丰满人妻一区二区三区视频av | 青草久久国产| 国产综合懂色| 精品国产亚洲在线| 国产亚洲欧美在线一区二区| 黄色片一级片一级黄色片| 91在线观看av| 亚洲片人在线观看| 观看美女的网站| 精品久久久久久久人妻蜜臀av| 欧美午夜高清在线| 日韩欧美精品免费久久 | 99精品在免费线老司机午夜| 久久久久久九九精品二区国产| 深爱激情五月婷婷| 国产精品野战在线观看| 久久久久久久久久黄片| xxxwww97欧美| 色在线成人网| 久久久久国产精品人妻aⅴ院| 欧美黄色片欧美黄色片| 美女高潮的动态| 99视频精品全部免费 在线| 亚洲无线在线观看| 级片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 真人一进一出gif抽搐免费| 91在线观看av| 97超视频在线观看视频| 成年女人毛片免费观看观看9| 一进一出抽搐gif免费好疼| 69av精品久久久久久| 欧美中文综合在线视频| 村上凉子中文字幕在线| 午夜福利视频1000在线观看| 成人国产综合亚洲| 亚洲狠狠婷婷综合久久图片| 中文在线观看免费www的网站| 色在线成人网| 性色avwww在线观看| 一级毛片高清免费大全| www日本黄色视频网| 禁无遮挡网站| 国产精品 国内视频| 俺也久久电影网| 天堂av国产一区二区熟女人妻| 精品久久久久久,| 久久亚洲真实| 夜夜躁狠狠躁天天躁| 国产精品影院久久| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 国产精品久久久久久亚洲av鲁大| 嫩草影院入口| 亚洲成av人片在线播放无| 免费av毛片视频| 国产精品1区2区在线观看.| av在线天堂中文字幕| 免费在线观看成人毛片| 757午夜福利合集在线观看| 美女高潮喷水抽搐中文字幕| 欧美一级毛片孕妇| ponron亚洲| 最近最新中文字幕大全免费视频| 色综合婷婷激情| 国产精品三级大全| 黄色丝袜av网址大全| 黄色日韩在线| 亚洲国产精品成人综合色| 18+在线观看网站| 99在线人妻在线中文字幕| 国产成人啪精品午夜网站| 国产一区二区三区视频了| 男女之事视频高清在线观看| 脱女人内裤的视频| 国产精品久久久久久亚洲av鲁大| 久久久久九九精品影院| 久久久久久久亚洲中文字幕 | а√天堂www在线а√下载| 99热精品在线国产| 99热只有精品国产| 国产在视频线在精品| 国产野战对白在线观看| 1024手机看黄色片| 少妇高潮的动态图| 他把我摸到了高潮在线观看| 免费在线观看成人毛片| 网址你懂的国产日韩在线| 午夜a级毛片| 久久久久性生活片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美一区二区三区黑人| 色综合婷婷激情| 国产av不卡久久| 亚洲av成人精品一区久久| 日本 av在线| 精品无人区乱码1区二区| 欧美最新免费一区二区三区 | 国产三级在线视频| 天天添夜夜摸| 色在线成人网| 亚洲五月天丁香| 亚洲国产欧美网| 99热只有精品国产| 中出人妻视频一区二区| 老熟妇仑乱视频hdxx| 国产美女午夜福利| 国产三级在线视频| 一卡2卡三卡四卡精品乱码亚洲| av中文乱码字幕在线| 亚洲国产色片| 精品不卡国产一区二区三区| a级毛片a级免费在线| 国产色婷婷99| 国内揄拍国产精品人妻在线| 亚洲人成电影免费在线| 小说图片视频综合网站| 欧美性猛交黑人性爽| 成人18禁在线播放| 精品人妻1区二区| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 亚洲无线在线观看| 尤物成人国产欧美一区二区三区| 岛国视频午夜一区免费看| 亚洲精品在线观看二区| 91久久精品国产一区二区成人 | 午夜久久久久精精品| 一进一出抽搐gif免费好疼| 精品人妻偷拍中文字幕| 搡老妇女老女人老熟妇| 黑人欧美特级aaaaaa片| 精品久久久久久,| 免费看十八禁软件| 国产精品乱码一区二三区的特点| 久久精品综合一区二区三区| 亚洲国产欧美网| 成人午夜高清在线视频| 精品国产超薄肉色丝袜足j| 日韩亚洲欧美综合| 国产午夜精品久久久久久一区二区三区 | 国产成人福利小说| 一个人看视频在线观看www免费 | 又粗又爽又猛毛片免费看| www国产在线视频色| 国产精品久久久久久精品电影| 国内精品久久久久精免费| 天天添夜夜摸| 欧美最新免费一区二区三区 | 性欧美人与动物交配| 黄色女人牲交| 午夜视频国产福利| 免费人成视频x8x8入口观看| 悠悠久久av| 国产精品亚洲av一区麻豆| 村上凉子中文字幕在线| 午夜福利视频1000在线观看| 午夜免费激情av| 亚洲成av人片在线播放无| 美女黄网站色视频| 制服人妻中文乱码| 99热只有精品国产| 99久久成人亚洲精品观看| 久久婷婷人人爽人人干人人爱| 99久久综合精品五月天人人| 嫩草影视91久久| 日本一二三区视频观看| 国产精品亚洲一级av第二区| 在线观看美女被高潮喷水网站 | 91九色精品人成在线观看| 嫩草影院入口| 中文字幕av成人在线电影| 精品午夜福利视频在线观看一区| 亚洲av免费在线观看| 欧美乱码精品一区二区三区| 中亚洲国语对白在线视频| ponron亚洲| 成熟少妇高潮喷水视频| 亚洲熟妇熟女久久| 色老头精品视频在线观看| 国产亚洲精品久久久com| 免费看a级黄色片| 尤物成人国产欧美一区二区三区| 久久久久久大精品| 免费av不卡在线播放| 熟女人妻精品中文字幕| 最新中文字幕久久久久| 国产欧美日韩精品一区二区| 偷拍熟女少妇极品色| 国产成+人综合+亚洲专区| 亚洲av电影在线进入| 淫妇啪啪啪对白视频| 88av欧美| 日本一本二区三区精品| 精品乱码久久久久久99久播| 最近最新中文字幕大全电影3| 亚洲激情在线av| 级片在线观看| 搡女人真爽免费视频火全软件 | 亚洲五月婷婷丁香| 九九在线视频观看精品| 久久天躁狠狠躁夜夜2o2o| av福利片在线观看| 少妇人妻精品综合一区二区 | 日韩欧美三级三区| 国产精品亚洲美女久久久| 午夜免费激情av| 久久人人精品亚洲av| 国产 一区 欧美 日韩| 18禁在线播放成人免费| 成年女人看的毛片在线观看| www国产在线视频色| 五月玫瑰六月丁香| 99久久99久久久精品蜜桃| 国内毛片毛片毛片毛片毛片| 又黄又粗又硬又大视频| 国产伦人伦偷精品视频| 搞女人的毛片| 欧美在线一区亚洲| 成人国产综合亚洲| 亚洲国产精品成人综合色| 老熟妇乱子伦视频在线观看| 免费搜索国产男女视频| 黄色日韩在线| 免费观看的影片在线观看| 亚洲欧美日韩无卡精品| 日韩大尺度精品在线看网址| 欧美zozozo另类| 老熟妇仑乱视频hdxx| 内地一区二区视频在线| 亚洲无线在线观看| 午夜两性在线视频| 狂野欧美激情性xxxx| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 成年人黄色毛片网站| www日本黄色视频网| 国产av麻豆久久久久久久| 国产精品 国内视频| av天堂在线播放| 欧美大码av| 搡老妇女老女人老熟妇| 老熟妇仑乱视频hdxx| 欧美国产日韩亚洲一区| a级毛片a级免费在线| 国产精品久久久久久久电影 | 丰满人妻一区二区三区视频av | 亚洲av日韩精品久久久久久密| 三级男女做爰猛烈吃奶摸视频| 夜夜夜夜夜久久久久| 国产精品一区二区三区四区久久| 日韩精品中文字幕看吧| 美女 人体艺术 gogo| 国产私拍福利视频在线观看| 久久精品国产亚洲av涩爱 | 免费人成在线观看视频色| 久久亚洲真实| 亚洲精品一区av在线观看| 久久久精品大字幕| 成人午夜高清在线视频| 成人国产一区最新在线观看| 久久人妻av系列| 精品免费久久久久久久清纯| 中文字幕av成人在线电影| 亚洲天堂国产精品一区在线| 男女之事视频高清在线观看| 亚洲专区中文字幕在线| 最近最新中文字幕大全电影3| 狠狠狠狠99中文字幕| 国内少妇人妻偷人精品xxx网站| 国产高清videossex| 性色av乱码一区二区三区2| 午夜福利高清视频| 国产综合懂色| 99国产精品一区二区蜜桃av| 午夜福利在线在线|