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

    長期粗放經(jīng)營毛竹林土壤微生物群落演變特征

    2017-04-19 08:45:01牛利敏苗倞婧彭定聰徐秋芳鄔奇峰
    林業(yè)科學(xué)研究 2017年2期
    關(guān)鍵詞:毛竹林毛竹馬尾松

    牛利敏,苗倞婧,彭定聰,徐秋芳,鄔奇峰,秦 華*

    (1.浙江農(nóng)林大學(xué)環(huán)境與資源學(xué)院,浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點實驗室,浙江 臨安 311300;2.浙江省臨安市農(nóng)業(yè)技術(shù)推廣中心,浙江 臨安 311300)

    長期粗放經(jīng)營毛竹林土壤微生物群落演變特征

    牛利敏1,苗倞婧1,彭定聰1,徐秋芳1,鄔奇峰2,秦 華1*

    (1.浙江農(nóng)林大學(xué)環(huán)境與資源學(xué)院,浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點實驗室,浙江 臨安 311300;2.浙江省臨安市農(nóng)業(yè)技術(shù)推廣中心,浙江 臨安 311300)

    [目的]通過分析土壤微生物生物量及群落結(jié)構(gòu)的演變趨勢,篩選影響土壤微生物群落的關(guān)鍵環(huán)境因子,揭示土壤微生物群落對毛竹林長期粗放經(jīng)營的響應(yīng)機理。[方法]選取不同粗放經(jīng)營年限(5 a、9 a、15 a、18 a)毛竹林,以天然馬尾松林(Masson pine, MP)作為對照,采用磷脂脂肪酸(Phospholipid fatty acids, PLFA)分析方法表征土壤微生物量及群落結(jié)構(gòu)。[結(jié)果]結(jié)果表明,毛竹林土壤微生物總PLFA含量以及細菌、真菌、放線菌等PLFA含量均顯著低于馬尾松林(P<0.05),但不同經(jīng)營年限毛竹林之間沒有顯著差異。土壤堿解氮及有效磷含量對土壤微生物總PLFA含量以及細菌、真菌、原生動物等PLFA含量影響顯著(P<0.05),而土壤堿解氮、pH值以及有機質(zhì)含量對放線菌PLFA含量影響顯著(P<0.05)。長期粗放經(jīng)營過程中毛竹林土壤微生物豐富度及多樣性均呈逐漸下降趨勢。非度量多維尺度轉(zhuǎn)換排序(Non-metric multidimensional scaling, NMDS)分析結(jié)果表明,毛竹林土壤微生物群落結(jié)構(gòu)與馬尾松林有明顯區(qū)分(R=0.388 1,P=0.009)。土壤含水量、堿解氮、有效磷以及pH值合計解釋了90.28%的微生物群落結(jié)構(gòu)變異量,其中土壤含水量、堿解氮、有效磷貢獻顯著(P<0.05)。[結(jié)論]長期粗放經(jīng)營降低了毛竹林土壤微生物量,改變了群落結(jié)構(gòu),其生態(tài)風(fēng)險還有待于進一步評估。

    毛竹;粗放經(jīng)營;土壤微生物;磷脂脂肪酸

    毛竹(Phyllostachysedulis(Carr.) H. de Lehaie)是我國分布最廣、面積最大和經(jīng)濟價值最高的竹種。第八次全國森林資源清查(2009—2013年)數(shù)據(jù)表明,我國竹林面積達601萬hm2,其中毛竹林面積為443萬hm2,占竹林總面積的73%,并且毛竹林面積還將繼續(xù)增加[1-2]。毛竹具有生長速度快、經(jīng)濟效益好等優(yōu)點,是山區(qū)農(nóng)民的重要收入來源。為了獲得更大的經(jīng)濟效益,近年來大面積馬尾松(PinusmassonianaLamb.)林以及杉木(Cunninghamialanceolata(Lamb.) Hook.)林等被毛竹林所替代。然而,已有研究表明,天然林受到毛竹入侵或人為改造成毛竹林后,毛竹易形成單優(yōu)群落,導(dǎo)致生物多樣性明顯下降[3-4]。同時,毛竹葉能釋放化感物質(zhì)抑制林下植被的生長,引起植物群落結(jié)構(gòu)及物種多樣性的改變[5]。近年來,毛竹林生物多樣性降低、生產(chǎn)力下降等現(xiàn)象開始引起人們的廣泛關(guān)注。

    土壤微生物是土壤養(yǎng)分循環(huán)的核心驅(qū)動力,在土壤形成和發(fā)育過程以及維護土壤生態(tài)系統(tǒng)穩(wěn)定性方面發(fā)揮著至關(guān)重要的作用[6]。地上植被對土壤微生物群落結(jié)構(gòu)形成起著至關(guān)重要的作用[7],植物種類改變會引起凋落物種類、根系分泌物以及養(yǎng)分吸收特性等方面的改變,從而導(dǎo)致土壤物理、化學(xué)和生物學(xué)性質(zhì)發(fā)生改變[8-9]。Chang等[10]研究表明,毛竹入侵日本柳杉(Cryptomeriajaponica(L. f.) D. Don)林后顯著改變了土壤微生物活性及群落結(jié)構(gòu)。徐秋芳等[11]對比了馬尾松林和毛竹林土壤微生物特性,發(fā)現(xiàn)馬尾松林土壤微生物量較高,多樣性更為豐富;而孫棣棣等[12]發(fā)現(xiàn)天然馬尾松林改種毛竹林后,土壤微生物量和多樣性變化不大,但微生物群落結(jié)構(gòu)發(fā)生了變化。目前,關(guān)于馬尾松林改造成毛竹林以及長期粗放經(jīng)營后毛竹林土壤微生物量及群落結(jié)構(gòu)的演變趨勢的研究結(jié)果尚不統(tǒng)一,而導(dǎo)致土壤微生物群落發(fā)生改變的關(guān)鍵控制因子也并不清楚,開展長期粗放經(jīng)營毛竹林土壤微生物群落演變趨勢及其控制因素的研究,有助于深入了解土壤微生物群落對環(huán)境因子的響應(yīng)機理,對于保護森林生態(tài)系統(tǒng)穩(wěn)定性具有重要意義。

    本研究在浙江安吉選取4個不同粗放經(jīng)營年限的毛竹林,同時選取附近的馬尾松林作為對照,通過分析土壤微生物生物量及群落結(jié)構(gòu)的演變趨勢,篩選影響土壤微生物群落的關(guān)鍵環(huán)境因子,揭示土壤微生物群落對毛竹林長期粗放經(jīng)營的響應(yīng)機理,為評估毛竹林生態(tài)系統(tǒng)穩(wěn)定性,促進毛竹林可持續(xù)發(fā)展提供理論依據(jù)。

    1 研究區(qū)概況

    2 材料與方法

    2.1 樣品采集與處理

    2.2 土壤含水量及理化性質(zhì)測定

    土壤含水量及理化性質(zhì)等參照魯如坤等[13]方法進行。其中,土壤含水量采用烘干法,新鮮土樣裝入鋁盒后于105℃烘干至恒質(zhì)量,稱質(zhì)量;土壤pH值測定采用1∶2.5土水比,用復(fù)合電極測定;有機質(zhì)(Organic matter, OM)含量采用重鉻酸鉀-硫酸外加熱法測定;堿解氮(Alkalytic nitrogen, AN)采用堿解擴散-硼酸吸收法測定;有效磷(Available phosphorus, AP)采用鹽酸-氟化銨提取,鉬銻抗比色法測定;速效鉀(Available potassium, AK)采用醋酸銨提取,火焰光度計測定。

    2.3 土壤磷脂脂肪酸分析

    土壤磷脂脂肪酸(Phospholipid fatty acids, PLFA)分析參照Frostegard等[14]方法,并在此基礎(chǔ)上進行改進[15]。主要步驟如下:(1)稱取3.00 g冷凍干燥的土樣,加入3.2 mL檸檬酸緩沖液、4 mL氯仿和8 mL甲醇,充分漩渦振蕩。避光條件下?lián)u床振蕩2 h,3 000 rpm條件下離心10 min。轉(zhuǎn)移上清液后重復(fù)提取一次?;旌虾玫纳锨逡赫袷幘鶆?,隔夜靜置后小心轉(zhuǎn)移下層氯仿相,用氮氣吹干。(2)準備好SPE柱(Supelco Inc., Bellefonte, PA, USA),8 mL氯仿洗柱后用氯仿溶解樣品,并全部轉(zhuǎn)移到SPE柱中,使樣品在柱中停留8 min。依次加入氯仿和丙酮洗中性脂和糖脂。然后用甲醇洗磷脂,小號試管收集。收集物用氮氣吹干,加200 μL內(nèi)標,渦旋后再氮氣吹干。(3)加1 mL甲醇甲苯混合液(甲醇:甲苯=1:1),溶解磷脂,再加1 mL 0.2 mol·L-1的氫氧化鉀甲醇溶液,37℃水浴15 min后,加入2 mL正己烷氯仿混合液(正己烷:氯仿=4:1)、0.3 mL 1 mol·L-1的醋酸、2 mL去離子水振蕩分離10 min,在3 000 rpm條件下離心5 min,轉(zhuǎn)移上層有機相到小號試管中。重復(fù)提取一次,合并有機相,氮氣吹干。(4)-20℃冷凍保存,在上機測定前用150 μL正己烷溶解。采用安捷倫6890N氣相色譜儀(Agilent Technologies Inc., USA)對提取的磷脂脂肪酸進行測定,MIDI Sherlocks微生物鑒定系統(tǒng)(MIDI, Inc., Newark, DE)進行磷脂脂肪酸鑒定。

    脂肪酸命名采用Petersen等[16]方法。所有檢測到的脂肪酸總量為土壤總磷脂脂肪酸含量。根據(jù)以往文獻的分類方法,對磷脂脂肪酸進行分類[17-18]。其中,14:0、15:0、16:0、17:0、i15:0、a15:0、i16:0、i17:0、a17:0、cy17:0、cy19:0、16:1ω7c、17:1ω8c、18:1ω5c、18:1ω7c以及19:0ω8c代表細菌;i14:0、i15:0、a15:0、i16:0、i17:0以及a17:0代表革蘭氏陽性細菌;cy17:0、16:1ω7c、18:1ω5c、18:1ω7c、17:1ω8c以及19:0ω8c代表革蘭氏陰性細菌;10Me16:0、10Me17:0以及10Me18:0代表放線菌;16:1ω5c、18:1ω9c 、18:2ω6,9c代表真菌;20:0代表原生動物。

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

    采用SPSS 18.0統(tǒng)計軟件進行數(shù)據(jù)處理,Duncan單因素方差分析比較各處理之間的差異顯著性(P<0.05);逐步回歸分析找出對土壤PLFA含量影響最顯著的土壤理化因子(解釋變量集)。采用R語言Vegan包計算不同處理土壤微生物群落豐富度、Shannon多樣性指數(shù)等;對不同處理土壤微生物群落進行非度量多維尺度轉(zhuǎn)換排序(non-metric multidimensional scaling, NMDS)分析,并通過相似性檢驗(Analysis of similarities, ANOSIM)測試不同群落之間的差異顯著性(P<0.05)。使用CANOCO 4.5軟件(Microcomputer Power, Ithaca, USA)對PLFA揭示的土壤微生物群落結(jié)構(gòu)與環(huán)境參數(shù)進行典范對應(yīng)分析(Canonical correspondence analysis, CCA),手動前向選擇環(huán)境因子,采用999次的蒙特卡羅排列檢驗(MonteCarlo permutation test 999 permutations, full model)進行顯著性檢驗。

    3 結(jié)果與分析

    3.1 土壤理化性質(zhì)分析

    研究結(jié)果表明,毛竹林土壤含水量均顯著低于馬尾松林(P<0.05),其中18 a毛竹林土壤水分含量最低(表1)。毛竹林在粗放經(jīng)營初期,土壤pH值呈現(xiàn)出逐漸升高的趨勢,9 a毛竹林土壤pH值顯著高于天然馬尾松林和其它年份毛竹林(P<0.05),但長期粗放經(jīng)營后土壤pH值又逐漸下降,18 a毛竹林與馬尾松林沒有顯著差異。在改種毛竹的初始幾年,土壤有機質(zhì)、有效磷、速效鉀等含量迅速增長,5 a毛竹林土壤有機質(zhì)、有效磷、速效鉀含量均顯著高于馬尾松林(P<0.05),但隨著粗放經(jīng)營年限的延長土壤有機質(zhì)及堿解氮及有效磷又呈現(xiàn)出下降趨勢,其中18 a毛竹林土壤有機質(zhì)含量顯著低于馬尾松林(P<0.05)。與其它養(yǎng)分不同的是,隨著粗放經(jīng)營年限的延長,土壤堿解氮含量呈現(xiàn)出下降趨勢,18 a毛竹林土壤堿解氮含量顯著低于馬尾松林(P<0.05)。

    表1 不同粗放經(jīng)營年限毛竹林土壤理化性質(zhì)

    Table 1 Soil physic-chemical properties under bamboo stands with different culture histories

    處理Treatment含水量SoilMoisture/%pH(H2O)有機質(zhì)Organicmatter/(g·kg-1)堿解氮AlkalyticN/(mg·kg-1)有效磷AvailableP/(mg·kg-11)速效鉀AvailableK/(mg·kg-1)MP23.71±0.60a4.51±0.05c43.28±2.18b115.21±21.88a0.64±0.15b40.0±7.2b5a18.91±1.45b4.61±0.08bc57.84±1.60a102.20±3.9ab1.97±0.41a67.0±8.2a9a20.42±1.00b4.88±0.15a35.75±1.25bc93.12±3.7ab1.95±0.43a41.0±5.0b15a19.08±1.99b4.71±0.12b41.87±2.13bc93.95±20.4ab1.49±0.85ab45.7±3.2b18a15.51±2.13c4.63±0.06bc35.01±4.22c70.15±5.1b1.23±0.25ab44.7±5.9b

    注: 數(shù)據(jù)為3個重復(fù)的平均值及標準差(mean±SD);同列中不同字母表示處理間差異達5%顯著水平。MP表示馬尾松林;5 a、9 a、15 a、18 a表示毛竹林粗放經(jīng)營年限為5年、9年、15年、18年。下同。

    Note: The data in the table are the means of 3 replicates and the standard error (mean±SD). Different letters in the same column indicate significant difference at 5% level. MP indicates masson pine. 5 a、9 a、15 a、18 a indicate moso bamboo stands with 5, 9, 15, 18 years of extensive management age, respectively. The same below.

    3.2 土壤PLFA含量

    利用PLFA分析方法表征不同樣地土壤微生物生物量及其群落,結(jié)果表明,馬尾松林土壤微生物總PLFA含量以及細菌、真菌、放線菌等PLFA含量均顯著高于毛竹林(P<0.05),而不同粗放經(jīng)營年限毛竹林土壤總PLFA含量之間沒有顯著差異(表2)。對于細菌,革蘭氏陽性細菌PLFA含量變化趨勢與總細菌相同,而15 a毛竹林土壤革蘭氏陰性細菌PLFA含量與馬尾松林沒有顯著差異。

    表2 不同粗放經(jīng)營年限毛竹林土壤PLFA含量

    將土壤微生物總PLFA含量作為因變量,土壤理化性質(zhì)作為自變量,進行逐步回歸分析,以揭示對不同種類微生物PLFA含量影響最顯著的一組變量(表3)。結(jié)果表明,引入的變量中,堿解氮和有效磷對土壤總PLFA含量以及細菌(包括革蘭氏陽性和陰性細菌)、真菌、原生動物PLFA含量等影響顯著,均達到1%極顯著水平(P<0.01)。對于真菌PLFA含量的影響,有效磷的貢獻大于堿解氮。與土壤微生物總PLFA含量以及其它類群微生物相比,土壤放線菌PLFA含量主要受堿解氮、土壤pH以及有機質(zhì)含量這組變量影響,同樣達到極顯著水平(P<0.01)。

    表3 逐步回歸分析

    注:AN、AP、OM分別為堿解氮、有效磷和土壤有機質(zhì)。

    Note: AN, AP, and OM indicate alkalytic N, available P and soil organic matter, respectively.

    3.3 土壤微生物群落多樣性

    馬尾松改種毛竹后,在短期內(nèi)土壤微生物物種豐富度沒有發(fā)生明顯變化,但是在長期尺度上有下降的趨勢,15 a和18 a毛竹林土壤微生物物種豐富度顯著低于馬尾松林以及5 a和9 a毛竹林(P<0.05),Shannon多樣性指數(shù)也有相似的變化趨勢(表4)。所有毛竹林土壤微生物群落Simpson多樣性指數(shù)均顯著低于馬尾松林(P<0.05)。與Shannon多樣性指數(shù)及Simpson多樣性指數(shù)相反,長期粗放經(jīng)營毛竹林土壤微生物群落Shannon均勻度及Simpson均勻度呈現(xiàn)出先下降再逐漸升高的趨勢,15 a和18 a毛竹林均顯著高于5 a及9 a毛竹林(P<0.05),其中18 a毛竹林Shannon均勻度和Simpson均勻度顯著高于對照馬尾松林(P<0.05)。

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

    注:MP表示馬尾松林;5 a、9 a、15 a、18 a表示毛竹林粗放經(jīng)營年限為5年、9年、15年、18年。下同Note: MP indicates masson pine. 5 a、9 a、15 a、18 a indicate moso bamboo stands with 5, 9, 15, 18 years of culture age, respectively. The same below.圖1 不同粗放經(jīng)營年限毛竹林土壤微生物群落非度量多維尺度分析Fig.1 Non-metric multidimensional scaling (NMDS) analysis of soil microbial community under bamboo stands with different culture histories

    對馬尾松林以及不同粗放經(jīng)營年限毛竹林土壤微生物群落進行NMDS分析(圖1),馬尾松林與毛竹林土壤微生物群落在NMDS第1排序軸上有明顯的區(qū)分,而不同粗放經(jīng)營年限毛竹林土壤微生物群落之間也存在差異,5 a、9 a毛竹林與15 a和18 a毛竹林在第2排序軸上區(qū)分明顯。相似性檢驗結(jié)果表明,各處理土壤微生物群落之間存在顯著差異(R=0.388 1,P=0.009)。典范對應(yīng)分析結(jié)果表明,土壤含水量(F=5.401,P=0.001)、有效磷(F=2.494,P=0.022)、堿解氮(F=3.150,P=0.036)以及pH值(F=2.046,P=0.062)對土壤微生物群落結(jié)構(gòu)的變異具有重要影響,合計解釋了90.28%的變異量(圖2),其中土壤含水量、有效磷及堿解氮含量對土壤微生物群落結(jié)構(gòu)變異影響顯著(P<0.05)。

    注:Soil moisture、AK、AP、AN、OM分別表示土壤含水量、速效鉀、有效磷、堿解氮、有機質(zhì)。Note: Soil moisture, AK, AP, AN, OM indicate soil moisture content, available K, available P, alkalytic N, organic matter, respectively圖2 不同粗放經(jīng)營年限毛竹林土壤微生物群落典范對應(yīng)分析Fig.2 Canonical correspondence analysis of soil microbial community under bamboo stands with different culture histories

    4 討論

    逐步回歸分析結(jié)果表明,除放線菌外,土壤微生物總PLFA含量以及各微生物種群PLFA含量受土壤堿解氮和有效磷影響顯著。王衛(wèi)霞等[23]發(fā)現(xiàn),在岷江干旱河谷造林過程中,隨著土壤氮含量的增加,土壤總PLFA含量、細菌、真菌、放線菌等PLFA含量均隨著增加,且與土壤全氮呈顯著正相關(guān),這與本試驗結(jié)果一致。施瑤等[24]研究了氮磷添加對內(nèi)蒙古典型溫帶草原土壤微生物群落的影響,發(fā)現(xiàn)隨著氮添加量的增加,土壤總PLFA含量以及細菌、放線菌等PLFA含量均顯著增加。盡管本研究中并沒有任何施肥措施,但之前大量研究表明,在一些氮含量缺乏的森林土壤中,氮素有效性是控制土壤微生物量的重要因子[25]。在磷限制條件下,如果土壤有效磷含量增加,能顯著提高土壤微生物生物量[26]。盡管馬尾松林改造成毛竹林后土壤有效磷含量均顯著增加,但總體處于較低水平,考慮到在毛竹生長過程中大量吸收氮素養(yǎng)分,土壤堿解氮含量下降幅度較磷素養(yǎng)分更為劇烈,因此作者推測在該土壤母質(zhì)條件下,土壤有效氮素含量應(yīng)該是控制微生物量的一個關(guān)鍵因子。與其它微生物種群不同,土壤放線菌PLFA含量不僅受土壤堿解氮影響,其還受到土壤pH值及有機碳含量影響,這可能是由于毛竹林下凋落物較馬尾松林少,從而導(dǎo)致依賴于分解土壤中植物凋落物中某些難分解組分的放線菌數(shù)量降低[27]。

    植被群落的改變會導(dǎo)致森林生態(tài)系統(tǒng)改變,從而直接影響土壤理化性質(zhì)以及微生物群落[8-10]。與馬尾松林相比,毛竹林土壤微生物群落Shannon多樣性指數(shù)和Simpson多樣性指數(shù)逐漸下降,而長期粗放經(jīng)營毛竹后土壤微生物群落Shannon多樣性指數(shù)和Simpson多樣性指數(shù)均顯著低于馬尾松林,說明地上植被群落多樣性減少影響了土壤微生物群落。植被群落改變后,由于不同的植物種類其凋落物性質(zhì)、根系分泌物以及養(yǎng)分吸收特性不同[7-9],如本研究區(qū)域中馬尾松林地上植被群落多樣性較豐富,而演變成毛竹林后,凋落物及根系分泌物主要來自于單一的毛竹,加之毛竹的養(yǎng)分吸收特性與普通喬灌木品種迥異,從而造成在長期時間尺度上土壤微生物群落多樣性指數(shù)顯著下降。與短期粗放經(jīng)營毛竹林(5 a、9 a)相比,長期粗放經(jīng)營毛竹林(15 a、18 a)土壤微生物豐富度以及基于Shannon多樣性指數(shù)和Simpson多樣性指數(shù)計算的均勻度指數(shù)均顯著降低(P<0.05)。有研究表明,當群落中物種發(fā)生變化時,Simpson多樣性指數(shù)對于富集種更加敏感,而Shannon多樣性指數(shù)對于稀疏種更加敏感[27]。作者推測在毛竹林經(jīng)營過程中,由于地上植被逐漸發(fā)生改變,導(dǎo)致土壤有機質(zhì)來源發(fā)生變化,一些稀疏種和富集種隨之發(fā)生改變,最終演變形成了適應(yīng)毛竹純林有機質(zhì)輸入特點(豐富的毛竹根系分泌物以及竹葉凋落物)的穩(wěn)定的土壤微生物群落。

    土壤微生物群落結(jié)構(gòu)的改變也部分驗證了作者以上推測,馬尾松林與短期粗放經(jīng)營(5 a、9 a)和長期粗放經(jīng)營(15 a、18 a)毛竹林土壤微生物群落結(jié)構(gòu)差異顯著。Chang等[10]利用PLFA方法研究了毛竹入侵日本柳杉林對土壤微生物群落結(jié)構(gòu)的影響,發(fā)現(xiàn)毛竹入侵后土壤微生物群落結(jié)構(gòu)發(fā)生了顯著的改變;Lin等[28]發(fā)現(xiàn)毛竹純林以及毛竹入侵杉木林后的混交林其土壤微生物群落比較類似,但顯著區(qū)別于杉木純林。研究表明毛竹具有強烈的化感作用,其本身也有抗菌活性[29],可以降低土壤微生物生物量并改變土壤微生物群落。這些已有的研究很好地解釋了本研究結(jié)果,說明毛竹林替代馬尾松林以及長期粗放經(jīng)營對土壤微生物群落結(jié)構(gòu)產(chǎn)生了重要影響。本研究發(fā)現(xiàn),土壤含水量、堿解氮及有效磷含量對微生物群落結(jié)構(gòu)改變具有顯著的貢獻。土壤水分含量是微生物群落活性的重要控制因素[30]。Brockett等[31]同樣利用PLFA方法研究了加拿大西部七個不同生物地理氣候區(qū)土壤微生物群落,發(fā)現(xiàn)土壤含水量是影響土壤微生物群落結(jié)構(gòu)以及酶活性的主要驅(qū)動因子;Morris等[32]發(fā)現(xiàn)土壤細菌和真菌生物量隨土壤含水量梯度變化顯著。由此作者推測,馬尾松林改造為毛竹林后,由于地表植被和凋落物層明顯減少,導(dǎo)致土壤水分含量隨之減少,從而對土壤微生物群落,特別是細菌和真菌群落產(chǎn)生重要影響。由于毛竹的生長特性,導(dǎo)致林下灌木雜草等難以生長,降低了地上植物群落多樣性,從而影響土壤微生物群落和養(yǎng)分循環(huán)過程?;诒狙芯拷Y(jié)果,作者建議適當加大粗放經(jīng)營毛竹林砍伐強度,創(chuàng)造有利于下層灌木和雜草生長的環(huán)境,從而提高植物多樣性和增加土壤有機質(zhì)輸入,提高土壤微生物多樣性和養(yǎng)分有效性,保障竹林可持續(xù)發(fā)展。

    5 結(jié)論

    (1)長期粗放經(jīng)營毛竹林土壤微生物總PLFA含量以及細菌、真菌、放線菌、原生動物等PLFA含量均顯著低于馬尾松林,土壤微生物量變化主要受堿解氮和有效磷含量變化影響。

    (2)與馬尾松林相比,長期粗放經(jīng)營毛竹林土壤微生物群落多樣性指數(shù)均顯著下降,群落結(jié)構(gòu)變化顯著。

    (3)毛竹林長期粗放經(jīng)營后導(dǎo)致的土壤含水量、堿解氮、有效磷含量變化是土壤微生物群落結(jié)構(gòu)變異的主要貢獻因子。

    [1] 何冬華,陳俊輝,徐秋芳,等.集約經(jīng)營對毛竹林土壤固氮細菌群落結(jié)構(gòu)和豐度的影響[J].應(yīng)用生態(tài)學(xué)報,2015,26(10):2961-2968.

    [2] 季海寶,莊舜堯,張厚喜,等.我國毛竹林生態(tài)系統(tǒng)碳儲量的地帶性差異[J].生態(tài)環(huán)境學(xué)報,2013,22(1):1-5.

    [3] 白尚斌,周國模,王懿祥,等.天目山保護區(qū)森林群落植物多樣性對毛竹入侵的響應(yīng)及動態(tài)變化[J].生物多樣性,2013,21(3):288-295.

    [4] Okutomi K, Shinoda S, Fukuda H. Causal analysis of the invasion of broad-leaved forest by bamboo in Japan[J]. Journal of Vegetation Science, 1996, 7(5): 723-728.

    [5] Larpkern P, Moe S R, Totland O. Bamboo dominance reduces tree regeneration in a disturbed tropical forest[J]. Oecologia, 2011, 165(1): 161-168.

    [6] 殷士學(xué).土壤微生物生物量及其與養(yǎng)分循環(huán)關(guān)系的研究進展[J].土壤學(xué)進展,1993,21(4):1-8.

    [7] Grayston S, Campbell C, Bardgett R,etal. Assessing shifts in microbial community structure across a range of grasslands of differing management intensity using CLPP, PLFA and community DNA techniques[J]. Applied Soil Ecology, 2004, 25(1): 63-84.

    [8] Chen G S, Yang Y S, Xie J S,etal. Soil biological changes for a natural forest and two plantations in subtropical China[J]. Pedosphere, 2004, 14(3): 297-304.

    [9] Lucas-Borja M E, Candel D, Jindo K,etal. Soil microbial community structure and activity in monospecific and mixed forest stands, under Mediterranean humid conditions[J]. Plant and Soil, 2012, 354(1-2): 359-370.

    [10] Chang E H, Chiu C Y. Changes in soil microbial community structure and activity in a cedar plantation invaded by moso bamboo[J]. Applied Soil Ecology, 2015, 91: 1-7.

    [11] 徐秋芳,張許昌,王緒南,等.毛竹林與馬尾松林土壤生物學(xué)性質(zhì)及微生物功能多樣性比較[J].浙江林業(yè)科技,2008,28(3):8-12.

    [12] 孫棣棣,徐秋芳,田 甜,等.不同栽培歷史毛竹林土壤微生物生物量及群落組成變化[J].林業(yè)科學(xué),2011,47(7):181-186.

    [13] 魯如坤. 土壤農(nóng)化分析[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000.

    [14] Frosteg?rd A, B??th E. The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil[J]. Biology and Fertility of Soils, 1996, 22(1-2): 59-65.

    [15] Wu Y P, Ma B, Zhou L,etal. Changes in the soil microbial community structure with latitude in eastern China, based on phospholipid fatty acid analysis[J]. Applied Soil Ecology, 2009, 43(2-3): 234-240.

    [16] Petersen S O, Klug M J. Effects of sieving, storage, and incubation temperature on the phospholipid fatty acid profile of a soil microbial community[J]. Applied and Environmental Microbiology, 1994, 60(7): 2421-2430.

    [17] Frosteg?rd ?, B??th E, Tunlio A. Shifts in the structure of soil microbial communities in limed forests as revealed by phospholipid fatty acid analysis[J]. Soil Biology and Biochemistry, 1993, 25(6): 723-730.

    [18] Moore-Kucera J, Dick R. PLFA profiling of microbial community structure and seasonal shifts in soils of a douglas-fir chronosequence[J]. Microbial Ecology, 2008, 55(3): 500-511.

    [19] 陳曉娟,吳小紅,劉守龍,等.不同耕地利用方式下土壤微生物活性及群落結(jié)構(gòu)特性分析:基于 PLFA 和 MicroRespTM方法[J].環(huán)境科學(xué),2013,34(6):2375-2382.

    [20] 李 嬌,蔣先敏,尹華軍,等.不同林齡云杉人工林的根系分泌物與土壤微生物[J].應(yīng)用生態(tài)學(xué)報,2014,25(2):325-332.

    [21] McDaniel M, Tiemann L, Grandy A. Does agricultural crop diversity enhance soil microbial biomass and organic matter dynamics? A meta-analysis[J]. Ecological Applications, 2014, 24(3): 560-570.

    [22] Sul W J, Asuming-Brempong S, Wang Q,etal. Tropical agricultural land management influences on soil microbial communities through its effect on soil organic carbon[J]. Soil Biology and Biochemistry, 2013, 65:33-38.

    [23] 王衛(wèi)霞,羅 達,史作民,等.岷江干旱河谷造林對土壤微生物群落結(jié)構(gòu)的影響[J].生態(tài)學(xué)報,2014,34(4):890-898.

    [24] 施 瑤,王忠強,張心昱,等.氮磷添加對內(nèi)蒙古溫帶典型草原土壤微生物群落結(jié)構(gòu)的影響[J].生態(tài)學(xué)報,2014,34(17):4943-4949.

    [25] Turner B L, Wright S J. The response of microbial biomass and hydrolytic enzymes to a decade of nitrogen, phosphorus, and potassium addition in a lowland tropical rain forest[J]. Biogeochemistry, 2014, 117(1):115-130.

    [26] 馬 朋,李昌曉,雷 明,等.三峽庫區(qū)岸坡消落帶草地, 棄耕地和耕地土壤微生物及酶活性特征[J].生態(tài)學(xué)報,2014,34(4):1010-1020.

    [27] 許 晴,張 放,許中旗,等.Simpson 指數(shù)和 Shannon-Wiener 指數(shù)若干特征的分析及 “稀釋效應(yīng)” [J].草業(yè)科學(xué),2011,28(4):527-531.

    [28] Lin Y T, Tang S L,Pai C W,etal. Changes in the soil bacterial communities in a cedar plantation invaded by moso bamboo[J]. Microbial Ecology, 2013, 67(2): 421-429.

    [29] Tanaka A, Shimizu K, Kondo R. Antibacterial compounds from shoot skins of moso bamboo (Phyllostachyspubescens)[J]. Journal of Wood Science, 2013, 59(2): 155-159.

    [30] Hackl E, Pfeffer M, Donat C,etal. Composition of the microbial communities in the mineral soil under different types of natural forest[J]. Soil Biology and Biochemistry, 2005, 37(4): 661-671.

    [31] Brockett B F T, Prescott C E, Grayston S J. Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada[J]. Soil Biology and Biochemistry, 2012, 44(1): 9-20.

    [32] Morris S J, Boerner R E J. Spatial distribution of fungal and bacterial biomass in southern Ohio hardwood forest soils: scale dependency and landscape patterns[J]. Soil Biology and Biochemistry, 1999, 31(6): 887-902.

    (責(zé)任編輯:金立新)

    Variation Patterns of Soil Microbial Community ofPhyllostachysedulisStands under Long-term Extensive Management

    NIULi-min1,MIAOJing-jing1,PENGDing-cong1,XUQiu-fang1,WUQi-feng2,QINHua1

    (1.School of Environmental and Resource Sciences, Zhejiang Agriculture and Forestry University, Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon Sequestration, Lin’an 311300, Zhejiang, China; 2.Agricultural Technology Extension Centre of
    Lin’an City, Zhejiang Province, Lin’an 311300, Zhejiang, China)

    [Objective]To investigate the effect of long-term extensive management on soil microbial communities and the related mechanisms.[Method]Moso bamboo (Phyllostachysedulis) stands with different extensive management time (5, 9, 15, and 18 years) were selected and compared with natural masson pine forest (PinusmassonianaLamb.) as control. The phospholipid fatty acids (PLFA) were analyzed to indicate the soil microbial biomass and community composition.[Result]Results showed that the total soil microbial PLFA contents as well as bacterial, fungal, actinomycic and protozoal PLFA contents under moso bamboo stands were all significant lower than that of masson pine forest (P<0.05), while no significant difference was found among bamboo stands. The total PLFA contents of soil microbial community as well as bacterial, fungal and protozoal PLFA contents were all significantly correlated with soil alkalytic nitrogen (AN) and available phosphorus (AP) contents (P<0.05), while the actinomycic PLFA content correlated with AN, organic matter contents and pH significantly (P<0.05). Soil microbial diversity, i.e. richness, Shannon diversity and Simpson diversity indexes, decreased progressively with the extensive management time of moso bamboo. According to the results of non-metric multidimensional scaling (NMDS), the soil microbial communities under moso bamboo and masson pine forest were significantly different (R=0.388 1,P=0.009). Furthermore, the soil microbial communities were different between the moso bamboo stands with short-term (5a, 9a) and long-term (15a, 18a) extensive management. The soil moisture, AN, AP and pH explained 90.28% of the total microbial community variation, among which the soil moisture, AN and AP contributed significantly to the microbial community variation (P<0.05). [Conclusion]The results indicate that long-term extensive management has a negative effect on soil microbial biomass and result in shifting of microbial community structure. However, the ecological risks of the forest substitution and long-term extensive management still need to be further assessed.

    Phyllostachysedulis; extensive management; soil microbial community; phospholipid fatty acids

    2016-05-27 基金項目: 國家自然科學(xué)基金項目(41271274、31570602);浙江農(nóng)林大學(xué)農(nóng)林碳匯與生態(tài)環(huán)境修復(fù)研究中心預(yù)研基金項目(2013CB03);浙江農(nóng)林大學(xué)創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃(20150107)。 作者簡介: 牛利敏(1991-),女,內(nèi)蒙古清水河人,碩士研究生,主要從事土壤微生物生態(tài)研究. * 通訊作者:秦 華,副教授,博士,碩士生導(dǎo)師. E-mail: qinhua@zafu.edu.cn

    10.13275/j.cnki.lykxyj.2017.02.014

    S795.7

    A

    1001-1498(2017)02-0285-08

    猜你喜歡
    毛竹林毛竹馬尾松
    自然封育條件下毛竹林內(nèi)凋落物及土壤持水能力變化研究
    節(jié)骨草和毛竹苗
    馬尾松栽培技術(shù)及撫育管理
    綠色科技(2019年5期)2019-11-29 13:17:37
    毛竹林經(jīng)營投入產(chǎn)出關(guān)系與經(jīng)營效益的相關(guān)性研究
    不同坡位對毛竹林小氣候及地上生長的影響
    馬尾松果糖-1,6-二磷酸酶基因克隆及表達模式分析
    馬尾松初級種子園復(fù)壯技術(shù)
    24年生馬尾松種子園自由授粉子代測定及家系選擇
    毛竹護筍養(yǎng)竹技術(shù)探討
    毛竹幼林四季管理技術(shù)
    高清毛片免费看| 亚洲av.av天堂| 91成人精品电影| 国产精品国产三级国产专区5o| 最后的刺客免费高清国语| 精品少妇内射三级| 男人和女人高潮做爰伦理| 国产69精品久久久久777片| 乱人伦中国视频| 精品久久久噜噜| kizo精华| 不卡视频在线观看欧美| 中文资源天堂在线| 国产精品不卡视频一区二区| 国产男女超爽视频在线观看| 国产在线免费精品| 国产在线免费精品| 麻豆乱淫一区二区| 久久 成人 亚洲| 国产精品熟女久久久久浪| 91午夜精品亚洲一区二区三区| 永久网站在线| 亚洲av欧美aⅴ国产| 亚洲精品一区蜜桃| 亚洲自偷自拍三级| 亚洲av国产av综合av卡| 国产精品人妻久久久久久| 久久久国产一区二区| 简卡轻食公司| 国产亚洲午夜精品一区二区久久| 秋霞伦理黄片| 久久鲁丝午夜福利片| 国产中年淑女户外野战色| 国产亚洲欧美精品永久| 亚洲精品视频女| 亚洲无线观看免费| 亚洲国产成人一精品久久久| av天堂久久9| 在线观看www视频免费| 另类精品久久| 特大巨黑吊av在线直播| 国产精品一二三区在线看| 人妻人人澡人人爽人人| 高清欧美精品videossex| 国产成人aa在线观看| 久久精品久久久久久久性| 一区二区三区乱码不卡18| 美女xxoo啪啪120秒动态图| 草草在线视频免费看| 国产有黄有色有爽视频| 97超碰精品成人国产| 久久久久人妻精品一区果冻| 中文字幕av电影在线播放| 日韩不卡一区二区三区视频在线| 亚洲欧美一区二区三区黑人 | 少妇高潮的动态图| 亚洲精品中文字幕在线视频 | 久久久久久伊人网av| 精品一品国产午夜福利视频| 国产精品伦人一区二区| 久久久久久久久久成人| 黄片无遮挡物在线观看| 激情五月婷婷亚洲| 噜噜噜噜噜久久久久久91| 精品卡一卡二卡四卡免费| 国产黄色视频一区二区在线观看| 欧美变态另类bdsm刘玥| 夜夜看夜夜爽夜夜摸| 国内揄拍国产精品人妻在线| 免费观看无遮挡的男女| 香蕉精品网在线| 久久av网站| 日韩欧美 国产精品| 国产欧美另类精品又又久久亚洲欧美| 国产淫片久久久久久久久| 久久精品国产鲁丝片午夜精品| 人妻制服诱惑在线中文字幕| 亚洲av成人精品一二三区| 国产黄片视频在线免费观看| 亚洲av不卡在线观看| 国产毛片在线视频| 国产精品偷伦视频观看了| 我的老师免费观看完整版| 国产极品粉嫩免费观看在线 | 黄色毛片三级朝国网站 | 免费人成在线观看视频色| 欧美另类一区| 国产乱来视频区| 国产免费又黄又爽又色| 9色porny在线观看| 色网站视频免费| 欧美xxxx性猛交bbbb| 尾随美女入室| 亚洲精品色激情综合| 国产高清三级在线| 97在线视频观看| 亚洲成色77777| av天堂久久9| 中国美白少妇内射xxxbb| 亚洲欧美日韩另类电影网站| 国产综合精华液| a级毛片在线看网站| 午夜影院在线不卡| 国产视频内射| 亚洲欧美成人综合另类久久久| 精品久久久噜噜| 国产欧美日韩一区二区三区在线 | 大香蕉久久网| 亚洲美女视频黄频| 午夜福利影视在线免费观看| 精品一区二区三卡| 最新的欧美精品一区二区| 免费av中文字幕在线| 欧美日韩av久久| 男女无遮挡免费网站观看| 日本-黄色视频高清免费观看| 美女国产视频在线观看| 国产精品嫩草影院av在线观看| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| 国产成人精品婷婷| 国产 精品1| 日日爽夜夜爽网站| 国产在线男女| 中文天堂在线官网| 日韩免费高清中文字幕av| 六月丁香七月| av卡一久久| 免费观看性生交大片5| 我要看日韩黄色一级片| 国产 一区精品| 不卡视频在线观看欧美| 国产精品一二三区在线看| 天堂8中文在线网| av福利片在线| 美女脱内裤让男人舔精品视频| 久久久久久久亚洲中文字幕| 亚洲国产精品成人久久小说| 桃花免费在线播放| av网站免费在线观看视频| 色5月婷婷丁香| 久久久久久久大尺度免费视频| 大片免费播放器 马上看| 爱豆传媒免费全集在线观看| freevideosex欧美| 看十八女毛片水多多多| 午夜影院在线不卡| 妹子高潮喷水视频| 18禁在线无遮挡免费观看视频| 日韩三级伦理在线观看| 国内少妇人妻偷人精品xxx网站| 国产精品国产av在线观看| 亚洲高清免费不卡视频| 自线自在国产av| 一级毛片电影观看| 亚洲精品日韩av片在线观看| 久久精品熟女亚洲av麻豆精品| 精品久久国产蜜桃| 青春草视频在线免费观看| 午夜老司机福利剧场| 欧美少妇被猛烈插入视频| 纵有疾风起免费观看全集完整版| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品国产av成人精品| 各种免费的搞黄视频| 国产有黄有色有爽视频| 两个人免费观看高清视频 | 成年美女黄网站色视频大全免费 | 热re99久久精品国产66热6| 午夜91福利影院| 国产日韩一区二区三区精品不卡 | 中文字幕亚洲精品专区| 一区二区三区四区激情视频| 国产成人免费观看mmmm| 永久网站在线| 99久久精品一区二区三区| 少妇人妻精品综合一区二区| 黑人猛操日本美女一级片| 亚洲av国产av综合av卡| 国产在线男女| 欧美一级a爱片免费观看看| 国产乱来视频区| 午夜福利视频精品| 国产男女内射视频| 亚洲成色77777| 亚洲av在线观看美女高潮| 黄色怎么调成土黄色| 高清黄色对白视频在线免费看 | 欧美日本中文国产一区发布| 国产精品蜜桃在线观看| 大片电影免费在线观看免费| 久久久久久伊人网av| 九色成人免费人妻av| 午夜福利视频精品| 欧美xxxx性猛交bbbb| 亚洲国产欧美日韩在线播放 | 久久免费观看电影| 欧美日韩视频精品一区| 免费不卡的大黄色大毛片视频在线观看| 国产精品欧美亚洲77777| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 久久久久久久久久久丰满| 极品人妻少妇av视频| 亚洲中文av在线| 免费少妇av软件| 国产亚洲最大av| 最新的欧美精品一区二区| www.av在线官网国产| 丁香六月天网| 女性被躁到高潮视频| 丰满乱子伦码专区| 亚洲欧美一区二区三区黑人 | 51国产日韩欧美| 国产精品久久久久久av不卡| 伦理电影免费视频| 高清av免费在线| 日韩av不卡免费在线播放| 国产高清三级在线| 欧美国产精品一级二级三级 | 国产视频内射| 自线自在国产av| 亚洲精品国产av成人精品| 国产视频首页在线观看| 一本—道久久a久久精品蜜桃钙片| 高清视频免费观看一区二区| av线在线观看网站| 午夜福利网站1000一区二区三区| 最近中文字幕高清免费大全6| 欧美成人午夜免费资源| www.色视频.com| 日韩精品免费视频一区二区三区 | 久久毛片免费看一区二区三区| 日韩欧美精品免费久久| 午夜91福利影院| 成人影院久久| 免费av中文字幕在线| 各种免费的搞黄视频| 26uuu在线亚洲综合色| 黄片无遮挡物在线观看| av网站免费在线观看视频| 寂寞人妻少妇视频99o| 国产精品久久久久久久电影| 伦精品一区二区三区| 国产一区二区三区av在线| 毛片一级片免费看久久久久| 人妻一区二区av| 99热网站在线观看| 男人狂女人下面高潮的视频| 伊人亚洲综合成人网| 久久6这里有精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品自拍成人| 最近手机中文字幕大全| 久久久久国产网址| 日本猛色少妇xxxxx猛交久久| 亚洲国产欧美日韩在线播放 | 99热这里只有是精品在线观看| 波野结衣二区三区在线| 人人澡人人妻人| 一个人免费看片子| 欧美bdsm另类| a级毛色黄片| 欧美三级亚洲精品| 91精品国产国语对白视频| a级一级毛片免费在线观看| 日韩伦理黄色片| 亚洲国产精品成人久久小说| 国产高清不卡午夜福利| 日韩强制内射视频| 日韩视频在线欧美| 少妇人妻久久综合中文| 欧美日韩亚洲高清精品| 亚洲一区二区三区欧美精品| 2021少妇久久久久久久久久久| 亚洲欧美成人精品一区二区| 国产成人精品婷婷| 又爽又黄a免费视频| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美清纯卡通| 五月开心婷婷网| 五月开心婷婷网| 少妇 在线观看| 最近最新中文字幕免费大全7| 国产亚洲5aaaaa淫片| 婷婷色综合大香蕉| 一本色道久久久久久精品综合| 久久综合国产亚洲精品| 狂野欧美白嫩少妇大欣赏| 热re99久久精品国产66热6| 久久久久久久亚洲中文字幕| 另类精品久久| 国产老妇伦熟女老妇高清| 青春草国产在线视频| av在线老鸭窝| 日韩,欧美,国产一区二区三区| 精品熟女少妇av免费看| 国产精品国产三级专区第一集| 极品人妻少妇av视频| 在线观看美女被高潮喷水网站| 国产精品一区二区在线不卡| 国内精品宾馆在线| 亚洲av在线观看美女高潮| 国产精品一区二区三区四区免费观看| 亚洲精品自拍成人| 日韩中文字幕视频在线看片| 久久久久久久精品精品| 国产精品熟女久久久久浪| 男男h啪啪无遮挡| 亚洲一级一片aⅴ在线观看| 麻豆精品久久久久久蜜桃| 亚洲欧美清纯卡通| 亚洲欧洲日产国产| 久久久久人妻精品一区果冻| 精品国产国语对白av| 国产乱人偷精品视频| 大陆偷拍与自拍| 久久久久视频综合| 欧美日韩视频高清一区二区三区二| 水蜜桃什么品种好| 久久狼人影院| 亚洲精品自拍成人| 狂野欧美激情性xxxx在线观看| 亚洲经典国产精华液单| 日本黄色片子视频| 亚洲第一av免费看| 两个人的视频大全免费| 久久影院123| 欧美日韩亚洲高清精品| 亚洲成人av在线免费| 国产精品国产三级国产av玫瑰| 偷拍熟女少妇极品色| 国产免费一级a男人的天堂| 看免费成人av毛片| 久久热精品热| 在线观看三级黄色| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 国产一区二区三区av在线| 99九九在线精品视频 | 国产精品国产三级专区第一集| 美女福利国产在线| 乱人伦中国视频| h日本视频在线播放| 国产一区二区在线观看av| 午夜91福利影院| 久久毛片免费看一区二区三区| 视频区图区小说| 最近手机中文字幕大全| 久久精品国产a三级三级三级| 高清黄色对白视频在线免费看 | 国产伦精品一区二区三区四那| 亚洲av欧美aⅴ国产| 免费大片18禁| 性色avwww在线观看| 国产精品女同一区二区软件| 成人漫画全彩无遮挡| 97在线视频观看| 一级毛片电影观看| 久久影院123| 国产精品一区二区在线不卡| 少妇猛男粗大的猛烈进出视频| 啦啦啦视频在线资源免费观看| 在线观看美女被高潮喷水网站| 黄色怎么调成土黄色| 有码 亚洲区| 国产精品三级大全| 日本猛色少妇xxxxx猛交久久| 亚洲国产av新网站| 在线观看美女被高潮喷水网站| 麻豆乱淫一区二区| 国产一区二区三区综合在线观看 | 国产日韩欧美在线精品| 蜜臀久久99精品久久宅男| 久久午夜福利片| 国产免费又黄又爽又色| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 熟妇人妻不卡中文字幕| 成人亚洲精品一区在线观看| 免费人妻精品一区二区三区视频| 深夜a级毛片| 91在线精品国自产拍蜜月| 欧美 日韩 精品 国产| 成人亚洲欧美一区二区av| 少妇的逼好多水| 成人毛片a级毛片在线播放| 深夜a级毛片| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 免费看光身美女| 七月丁香在线播放| 51国产日韩欧美| 伦理电影大哥的女人| 亚洲内射少妇av| 中文字幕制服av| 国产免费一级a男人的天堂| 在线观看三级黄色| 国产成人免费观看mmmm| 韩国高清视频一区二区三区| 日韩电影二区| 国产 一区精品| av在线播放精品| 国产精品麻豆人妻色哟哟久久| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产| 国产有黄有色有爽视频| 观看美女的网站| 成年女人在线观看亚洲视频| 少妇被粗大猛烈的视频| 女性被躁到高潮视频| 高清视频免费观看一区二区| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 成人亚洲精品一区在线观看| 亚洲av电影在线观看一区二区三区| tube8黄色片| 黄片无遮挡物在线观看| 免费观看的影片在线观看| 狂野欧美激情性bbbbbb| a级一级毛片免费在线观看| 久久精品久久久久久噜噜老黄| 狂野欧美白嫩少妇大欣赏| 制服丝袜香蕉在线| 久久午夜福利片| 99国产精品免费福利视频| 久久久久久久久大av| 一区二区三区乱码不卡18| 在线观看免费视频网站a站| 亚洲国产av新网站| 最近的中文字幕免费完整| 97在线人人人人妻| 久久99蜜桃精品久久| 亚洲综合精品二区| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 色视频www国产| 亚洲中文av在线| 午夜激情久久久久久久| 伊人久久国产一区二区| 中文欧美无线码| 简卡轻食公司| 免费观看在线日韩| 久久精品久久精品一区二区三区| 午夜免费男女啪啪视频观看| 日韩制服骚丝袜av| 男人爽女人下面视频在线观看| 熟女av电影| 丰满人妻一区二区三区视频av| 青春草国产在线视频| 插逼视频在线观看| 久久国产亚洲av麻豆专区| 另类亚洲欧美激情| 日韩大片免费观看网站| a级毛色黄片| 寂寞人妻少妇视频99o| 成人国产av品久久久| 搡女人真爽免费视频火全软件| 老司机亚洲免费影院| 熟女av电影| 亚州av有码| 韩国av在线不卡| 高清视频免费观看一区二区| 国产成人免费观看mmmm| 18禁动态无遮挡网站| 全区人妻精品视频| 国产亚洲最大av| 性色av一级| 日日摸夜夜添夜夜添av毛片| av网站免费在线观看视频| 欧美日韩视频高清一区二区三区二| 欧美3d第一页| 色网站视频免费| 一本色道久久久久久精品综合| 亚洲婷婷狠狠爱综合网| 日本欧美视频一区| 国产成人精品久久久久久| av在线app专区| 亚洲欧美成人精品一区二区| 国产亚洲91精品色在线| 婷婷色av中文字幕| 中文天堂在线官网| 一二三四中文在线观看免费高清| 亚洲丝袜综合中文字幕| 少妇人妻久久综合中文| 久久99热6这里只有精品| 国产精品一区二区三区四区免费观看| 国产免费又黄又爽又色| 国产日韩欧美亚洲二区| 国产欧美日韩综合在线一区二区 | 免费黄色在线免费观看| 日本免费在线观看一区| 一区二区三区免费毛片| 国产在线免费精品| 亚洲精品久久午夜乱码| 内射极品少妇av片p| 一级二级三级毛片免费看| 欧美变态另类bdsm刘玥| 一级片'在线观看视频| 一级毛片aaaaaa免费看小| 综合色丁香网| 免费观看av网站的网址| 丰满饥渴人妻一区二区三| 成人国产麻豆网| 乱人伦中国视频| 七月丁香在线播放| 十八禁高潮呻吟视频 | 久久99一区二区三区| 精品卡一卡二卡四卡免费| 久久鲁丝午夜福利片| 欧美性感艳星| kizo精华| 多毛熟女@视频| 久久精品熟女亚洲av麻豆精品| 在线观看三级黄色| av在线app专区| 国产成人免费观看mmmm| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一二三区| 九九爱精品视频在线观看| 少妇人妻 视频| 国产爽快片一区二区三区| 久久精品久久久久久噜噜老黄| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 自线自在国产av| h日本视频在线播放| av免费在线看不卡| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 性色av一级| 日本色播在线视频| 国产一级毛片在线| 街头女战士在线观看网站| 久久精品国产亚洲网站| av播播在线观看一区| 美女视频免费永久观看网站| 日本wwww免费看| 亚洲av在线观看美女高潮| 亚洲欧美成人综合另类久久久| 另类精品久久| 成人黄色视频免费在线看| 久久99热6这里只有精品| 欧美日韩在线观看h| 欧美丝袜亚洲另类| av播播在线观看一区| 成人二区视频| 又大又黄又爽视频免费| 又爽又黄a免费视频| 黑人高潮一二区| 成年女人在线观看亚洲视频| 亚洲,一卡二卡三卡| 亚洲av欧美aⅴ国产| 欧美成人午夜免费资源| 多毛熟女@视频| 一二三四中文在线观看免费高清| 精品人妻偷拍中文字幕| 三上悠亚av全集在线观看 | 日日啪夜夜爽| 国产日韩一区二区三区精品不卡 | 边亲边吃奶的免费视频| 免费看日本二区| 精品久久久久久久久亚洲| 高清毛片免费看| 97超视频在线观看视频| videos熟女内射| 色婷婷av一区二区三区视频| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 丁香六月天网| 午夜久久久在线观看| 秋霞伦理黄片| 高清av免费在线| 欧美高清成人免费视频www| 一级毛片久久久久久久久女| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| 成人亚洲精品一区在线观看| 亚洲精品国产av成人精品| 寂寞人妻少妇视频99o| 51国产日韩欧美| 国产真实伦视频高清在线观看| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 日韩成人伦理影院| 老司机影院成人| 视频区图区小说| 黄色欧美视频在线观看| 夫妻午夜视频| 不卡视频在线观看欧美| 国产老妇伦熟女老妇高清| 又大又黄又爽视频免费| 蜜桃在线观看..| 大香蕉97超碰在线| av免费观看日本| 一级毛片aaaaaa免费看小| 在线观看美女被高潮喷水网站| 色94色欧美一区二区| 老司机亚洲免费影院| 欧美日韩精品成人综合77777| 成年女人在线观看亚洲视频| 免费观看性生交大片5| 久久久久久久久久人人人人人人| 99久国产av精品国产电影| 亚洲熟女精品中文字幕| 99久久精品一区二区三区| 人妻人人澡人人爽人人| 一区在线观看完整版| 成人午夜精彩视频在线观看| 国产黄片视频在线免费观看| 日韩熟女老妇一区二区性免费视频| 午夜日本视频在线| 一本大道久久a久久精品| 午夜福利网站1000一区二区三区| 狂野欧美激情性bbbbbb|