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

    稀土礦廢棄地植被恢復(fù)過程中土壤微生物演變

    2019-10-23 11:34:18李啟艷翁炳霖李宗勛趙雅曼陳順鈺侯曉龍蔡麗平
    中國(guó)環(huán)境科學(xué) 2019年10期
    關(guān)鍵詞:廢棄地稀土礦菌門

    李啟艷,翁炳霖,李宗勛,趙雅曼,陳順鈺,侯曉龍,2,3*,蔡麗平,2,3

    稀土礦廢棄地植被恢復(fù)過程中土壤微生物演變

    李啟艷1,翁炳霖1,李宗勛1,趙雅曼1,陳順鈺1,侯曉龍1,2,3*,蔡麗平1,2,3

    (1.福建農(nóng)林大學(xué)林學(xué)院,福建 福州 350002;2.南方紅壤區(qū)水土保持國(guó)家林業(yè)和草原局重點(diǎn)實(shí)驗(yàn)室,福建 福州 350002;3.海峽兩岸紅壤區(qū)水土保持協(xié)同創(chuàng)新中心,福建 福州 350002)

    以福建長(zhǎng)汀典型離子型稀土礦堆浸冶煉廢棄地作為研究對(duì)象,采取“空間代時(shí)間”的方法,對(duì)堆浸廢棄地不同植物恢復(fù)年限樣地土壤理化性質(zhì)及微生物進(jìn)行調(diào)查測(cè)定,分析植物恢復(fù)對(duì)稀土礦堆浸廢棄地土壤理化特性及微生物的影響.結(jié)果表明:在植被恢復(fù)第3a和4a土壤理化性質(zhì)顯著高于未治理廢棄地,但仍顯著低于礦區(qū)周邊未開采對(duì)照(<0.05).植被恢復(fù)年4a以后,長(zhǎng)汀稀土礦廢棄地土壤理化性質(zhì)則呈下降趨勢(shì),但仍顯著高于未治理廢棄地(<0.05).長(zhǎng)汀稀土礦堆浸廢棄地不同植被恢復(fù)年限土壤細(xì)菌豐度及多樣存在一定差異.隨植被恢復(fù)年限的增加,土壤細(xì)菌多樣性及豐度降低,但一些特殊細(xì)菌群如-變形菌綱(Gammaproteobacteria)、-變形桿菌綱(Alphaproteobacteria)、酸桿菌門(Acidobacteria)、藍(lán)藻菌門(Cyanobacteria)等在不同植被恢復(fù)年限樣地中其相對(duì)豐度明顯增加.具有固氮作用的藍(lán)藻菌門(Cyanobacteria)隨著植被恢復(fù)年限的增加,其相對(duì)豐度降低.廣古菌門(Euryarchaeota僅存在于未開采對(duì)照樣地;土壤寡營(yíng)養(yǎng)細(xì)菌酸桿菌門(Acidobacteria)隨植被恢復(fù)年限的增加,相對(duì)豐度呈增大趨勢(shì).土壤中優(yōu)勢(shì)菌群廣古菌門(Euryarchaeota)、厚壁菌門(Firmicutes)、?變形菌門(Proteobacteria)、屬、羅爾斯通菌屬()與土壤pH、全磷、全鉀、速效磷等呈正相關(guān).土壤pH值、全磷、全鉀等是影響廣古菌門(Euryarchaeota)、厚壁菌門(Firmicutes)等細(xì)菌菌群多樣性及相對(duì)豐度變化的重要因素.以上結(jié)果表明,離子型稀土礦堆浸廢棄地植被恢復(fù)治理到3~4a時(shí),應(yīng)采取防治措施對(duì)其進(jìn)行人工撫育管理,否則可能會(huì)出現(xiàn)前期恢復(fù)后期退化的問題.

    稀土礦廢棄地;生態(tài)恢復(fù);微生物演變;16S rRNA;冗余分析

    我國(guó)離子型稀土(重稀土)主要分布于福建、江西、廣東等省區(qū),其中福建省離子型稀土儲(chǔ)存量居首位[1].長(zhǎng)汀縣作為福建省離子型稀土資源儲(chǔ)量最多的地區(qū),其稀土產(chǎn)業(yè)發(fā)展最早.稀土資源大多分布在嚴(yán)重水土流失區(qū),由于非法采礦及開采工藝落后破壞和占用了大量土地資源,加劇水土流失,使礦區(qū)生態(tài)壞境趨于惡化.

    眾多學(xué)者在稀土礦開采和冶煉對(duì)生態(tài)環(huán)境的危害方面開展了不少研究,王友生等[2]對(duì)長(zhǎng)汀稀土礦廢棄地重金屬污染的評(píng)價(jià)表明其綜合潛在生態(tài)風(fēng)險(xiǎn)指數(shù)已達(dá)極重污染水平.趙永紅等[3]對(duì)贛南已開采65a年的鎢礦礦區(qū)進(jìn)行土壤重金屬污染研究發(fā)發(fā)現(xiàn),土壤中重金屬含量均高出當(dāng)?shù)赝寥辣尘爸?稀土礦區(qū)的生態(tài)惡化問題引起社會(huì)各界的重視,許多專家學(xué)者積極開展稀土礦廢棄地植被恢復(fù)研究,配置各類喬、灌、草模式進(jìn)行植被恢復(fù).植被修復(fù)由于二次污染小,操作簡(jiǎn)便且成本低廉,是當(dāng)前礦山開采區(qū)土壤污染修復(fù)研究的熱點(diǎn).簡(jiǎn)麗華[4]使用香根草、鴨拓草、類蘆、寬葉雀稗等草種形成良好的草種配置模式,能在前期迅速覆蓋地面,減少地表徑流,起到固坡作用.翁伯琦等[5]則采用工程措施與生物措施相結(jié)合的辦法,按照喬灌1:17的配置比,強(qiáng)化人工恢復(fù)效果,控制水土流失.于君寶等[6]對(duì)撫順老虎臺(tái)煤礦土壤復(fù)墾發(fā)現(xiàn),覆土年齡越大其營(yíng)養(yǎng)元素越高.郭李凱[7]對(duì)不同復(fù)墾年限的煤矸礦重構(gòu)土壤養(yǎng)分研究,結(jié)果表明土壤有機(jī)質(zhì)、全氮、有效磷隨復(fù)墾年限增加呈增加趨勢(shì).陜永杰等[8]通過對(duì)平朔安太堡露天煤礦區(qū)開采前后土壤質(zhì)量研究發(fā)現(xiàn),礦區(qū)開采后土壤理化性質(zhì)呈毀滅性破壞,土壤質(zhì)量處于嚴(yán)重退化階段;土壤復(fù)墾8a后土壤質(zhì)量與未受擾動(dòng)土壤質(zhì)量相比差距較大.王友生等[9]通過對(duì)長(zhǎng)汀稀土礦取土場(chǎng)、堆浸廢棄地進(jìn)行不同的喬灌草配置植被恢復(fù)模式,篩選出了“寬葉雀稗+胡枝子+木荷+楓香+山杜英”是取土場(chǎng)較好的植被恢復(fù)模式,在長(zhǎng)汀稀土礦廢棄地得到了推廣應(yīng)用.但是,本課題組在對(duì)長(zhǎng)汀稀土礦植被恢復(fù)效果調(diào)查研究中發(fā)現(xiàn),植被恢復(fù)前期稀土礦廢棄地生態(tài)恢復(fù)效果較好,但隨著植被恢復(fù)時(shí)間的增加,出現(xiàn)植被生長(zhǎng)不良甚至死亡的現(xiàn)象.這可能在植被恢復(fù)后期存在生態(tài)退化的問題,如果不采取防治措施,前期治理成績(jī)將會(huì)功虧一簣.由于目前對(duì)長(zhǎng)汀稀土廢棄地生態(tài)恢復(fù)過程中土壤理化特性、微生物等的變化規(guī)律,以及生態(tài)恢復(fù)后期植被退化機(jī)制的研究報(bào)道不多.鑒于此,本課題組以長(zhǎng)汀稀土礦廢棄地不同植被恢復(fù)年限的固定樣地為研究對(duì)象,從微生物的角度揭示稀土礦廢棄地植被恢復(fù)過程中土壤微環(huán)境的演化規(guī)律.采用“空間代時(shí)間”的方法,對(duì)采用“楓香+木荷+胡枝子+寬葉雀稗”模式恢復(fù)0(未治理廢棄地)、3、4、6、7、9、10a樣地土壤理化特性及微生物組成進(jìn)行調(diào)查測(cè)定,比較不同植被恢復(fù)年限土壤理化特性的差異,闡明稀土礦廢棄地土壤微生物與土壤理化特性的相關(guān)性,探討稀土礦廢棄地生態(tài)恢復(fù)后期植被退化的機(jī)制,為稀土礦廢棄地生態(tài)恢復(fù)效果的長(zhǎng)期維持及退化防治措施的實(shí)施提供科學(xué)依據(jù).

    1 研究方法

    1.1 試驗(yàn)地概況

    試驗(yàn)地位于長(zhǎng)汀縣河田鎮(zhèn),116°14′~116°33′E, 25°33′~25°48′N,海拔在260~483m,低山陵地帶.年均氣溫在18~20.5℃之間,極端最高氣溫 38.9℃,極端最低氣溫-4.9℃.年降水量1686.7mm,無(wú)霜期為261d左右.

    稀土冶煉方法主要有原地浸析法、池浸法和堆浸法,正規(guī)開采主要以前兩者為主.早期主要采用原地浸析法,2000年以后絕大部分采用簡(jiǎn)易堆浸法.堆浸廢棄地是將地面表土層剝離后搬運(yùn)至由黏土壘積成的仿浸池中,用1%~2%的硫酸銨反復(fù)淋洗浸提,最后被遺棄的土壤,一般面積為500m2.未經(jīng)過植被恢復(fù)的堆浸廢棄地土壤容重為1.77±0.85g/cm3,質(zhì)量含水量為138.11±9.3g/kg,毛管持水量為726.71± 8.21g/kg,有機(jī)質(zhì)僅為2.83±0.21g/kg,pH值為4.05± 0.02,土壤理化性質(zhì)表現(xiàn)較差[9].

    1.2 試驗(yàn)設(shè)計(jì)與土壤樣品采集

    2007、2010及2013年以“楓香+木荷+胡枝子+寬葉雀稗”為模式,分別對(duì)稀土礦堆浸廢棄地進(jìn)行植被恢復(fù),并建立了固定觀測(cè)樣地,3個(gè)固定樣地的坡度、坡向、海拔等條件基本一致.同時(shí)設(shè)置遠(yuǎn)離稀土礦區(qū)天然植被樣地(CK)作為對(duì)照,每個(gè)樣地設(shè)置3個(gè)20m×20m標(biāo)準(zhǔn)樣方;未治理樣地(M0)數(shù)據(jù)參照王友生等[9].分別在2016年植被恢復(fù)9a(M9)、6a(M6)、3a(M3)和2017年植被恢復(fù)10a(M10)、7a(M7)、4a(M4))對(duì)土壤進(jìn)行調(diào)查采樣.每個(gè)樣方內(nèi)按上、中、下部各3個(gè)采樣點(diǎn)挖取土壤剖面,用環(huán)刀采集0~20cm土樣,用于分析土壤水分物理性質(zhì);各采樣點(diǎn)采集1kg土樣混合均勻帶回實(shí)驗(yàn)室,風(fēng)干研磨過2mm尼龍篩,用于土壤理化性質(zhì)分析.在每個(gè)樣方內(nèi)用10cm×10cm×10cm鋁盒采集0~20cm土樣,用于土壤水穩(wěn)性團(tuán)聚體分析.2016年在每塊試驗(yàn)樣地內(nèi)采集0~20cm土樣,低溫保存,送北京諾禾致源公司進(jìn)行16S rRNA測(cè)序.

    1.3 測(cè)定方法

    1.3.1 土壤理化性質(zhì)測(cè)定 土壤容重、持水量等物理性質(zhì)采用環(huán)刀法測(cè)定;土壤pH值、全N、全P、全K、速效磷、速效鉀、有機(jī)質(zhì)參照中華人民共和國(guó)林業(yè)行業(yè)標(biāo)準(zhǔn)中森林土壤分析方法進(jìn)行測(cè)定[10].土壤pH值采用電位法水土比2.5:1測(cè)定;全N用德國(guó)Elementar, Vario元素分析儀測(cè)定;全P采用氫氧化鈉堿熔—鉬銻抗比色測(cè)定;全K采用氫氧化鈉堿熔—原子吸收分光光度法測(cè)定;速效磷采用氟氯氨、鹽酸浸提-鉬銻抗比色法測(cè)定;速效鉀采用四苯硼鈉比濁法(1mol-1NaNO3浸提法)測(cè)定;土壤有機(jī)質(zhì)采用重鉻酸鉀-外加熱.法測(cè)定

    1.3.2 土壤細(xì)菌測(cè)定 土壤微生物總DNA采用 CTAB方法對(duì)樣本的基因組 DNA 進(jìn)行提取之后利用瓊脂糖凝膠電泳檢測(cè)DNA的純度和濃度,取適量的樣品于離心管中,使用無(wú)菌水稀釋樣品至1ng/μL.以稀釋后的基因組DNA為模板,根據(jù)測(cè)序區(qū)域的選擇,使用帶 Barcode 的特異引物,New England Biolabs公司的Phusion? High-Fidelity PCR Master Mix with GC Buffer,和高效高保真酶進(jìn)行PCR,確保擴(kuò)增效率和準(zhǔn)確性.引物對(duì)應(yīng)區(qū)域?yàn)?6S V4區(qū)引物(515F和806R).PCR產(chǎn)物使用2%濃度的瓊脂糖凝膠進(jìn)行電泳檢測(cè);根據(jù)PCR產(chǎn)物濃度進(jìn)行等量混樣,充分混勻后使用2%的瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物,對(duì)目的條帶使用qiagen公司提供的膠回收試劑盒回收產(chǎn)物.細(xì)菌PCR反應(yīng)條件如下:98℃預(yù)變性 3min,98℃變性45s, 55℃退火45s,72℃延伸 45s,30 個(gè)循環(huán)后72℃延伸 7min,最后保持在4℃環(huán)境下.使TruSeq? DNA PCR-Free Sample Preparation Kit建庫(kù)試劑盒進(jìn)行文庫(kù)構(gòu)建,構(gòu)建好的文庫(kù)經(jīng)過Qubit和Q-PCR定量,文庫(kù)合格后,使用HiSeq2500PE250進(jìn)行上機(jī)測(cè)序.

    1.4 數(shù)據(jù)分析方法

    用SPSS 19.0統(tǒng)計(jì)軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析;使用單因素方差分析(One-way ANOVA)進(jìn)行差異顯著性分析,用Duncan’s法進(jìn)行多重比較,顯著性水平為0.05.運(yùn)用CANOCO4.5軟件對(duì)土壤微生物與環(huán)境因子進(jìn)行冗余分析 (RDA) ;運(yùn)用Origin9.0軟件繪制分析圖.

    2 結(jié)果與討論

    2.1 土壤物理性質(zhì)變化比較

    表1 不同植被恢復(fù)年限土壤物理性質(zhì)

    注:表中不同小寫字母表示 0.05 水平 LSD 多重比較結(jié)果.

    由表1可知,隨植被恢復(fù)年限的增加,長(zhǎng)汀稀土礦堆浸廢棄地土壤容重呈先顯著降低(3~4a)(< 0.05),然后又顯著增大(6~10a),不同植被恢復(fù)年限土壤容重均顯著小于未治理區(qū)(<0.05),但仍顯著大于未開采對(duì)照(<0.05);隨植被恢復(fù)年限的增加,土壤含水量水量、最大持水量,毛管持水量、最小持水量、毛管孔隙、非毛管孔隙、均呈先顯著增大(3~4a),然后又顯著減小的規(guī)律(6~10a),但不同植被恢復(fù)年限這些物理性質(zhì)均顯著小于未開采對(duì)照,說明土壤仍未恢復(fù)到自然狀態(tài).在植被恢復(fù)前期3~4a土壤容重平均值為1.32g/cm3,接近自然狀態(tài)值,這與李兆龍[11]的研究結(jié)論相一致.隨著植被恢復(fù)年限的加長(zhǎng),土壤容重不斷增大,植被恢復(fù)10a土壤容重達(dá)到1.49g/ cm3.據(jù)研究表明:土壤容重在1.25~1.35g/cm3之間,土壤理化性質(zhì)表現(xiàn)良好[12].土壤容重大于1.4g/cm3時(shí),土壤機(jī)械阻力對(duì)植物根系的伸展發(fā)育有較大影響[13].

    2.2 不同植被恢復(fù)年限土壤化學(xué)性質(zhì)比較

    2.2.1 不同植被恢復(fù)年限土壤pH值 由圖1可知,隨著植被恢復(fù)年限的增加,稀土礦堆浸廢棄地土壤pH值與未治理土壤相比呈現(xiàn)出顯著增加趨勢(shì)(3~4a)(<0.05),然后又顯著降低(6~10a).不同植被恢復(fù)年限土壤pH值均顯著高于未治理區(qū),但仍顯著低于未開采區(qū)土壤(<0.05).結(jié)果表明,不同植被恢復(fù)年限土壤pH值未恢復(fù)到自然水平狀態(tài).

    圖1 不同植被恢復(fù)年限土壤pH

    圖中不同小寫字母表示不同植被恢復(fù)年限間差異顯著(<0.05).下同

    2.2.2 不同植被恢復(fù)年限土壤有機(jī)質(zhì) 由圖2可知,不同植被恢復(fù)年限土壤有機(jī)質(zhì)與未治理土壤相比顯著增加(3~4a),然后又顯著降低(6~10a)(<0.05).不同植被恢復(fù)年限土壤有機(jī)質(zhì)均顯著高于未治理土壤,但又顯著低于未開采對(duì)照土壤有機(jī)質(zhì)(<0.05).結(jié)果說明,隨著植被恢復(fù)年限的增加,土壤有機(jī)質(zhì)的積累并未恢復(fù)到自然狀態(tài).隨著植被恢復(fù)年限的增加,土壤pH值變化范圍為4.0~4.5,土壤有機(jī)質(zhì)總體表現(xiàn)為隨時(shí)間加長(zhǎng)而下降的趨勢(shì),總體顯著高于M0.與劉子壯[14]、彭東海[15]、毛蘭花[16]等研究結(jié)論不一致,可能的原因是:前期為使植物適應(yīng)惡劣的生長(zhǎng)環(huán)境,都施以一定肥料輔助其生長(zhǎng),對(duì)土壤養(yǎng)分消耗較少,但隨著植被恢復(fù)年限加長(zhǎng),前期施用的肥料肥效逐漸降低,植物生長(zhǎng)對(duì)土壤養(yǎng)分元素(N、P)質(zhì)量分?jǐn)?shù)的影響效應(yīng)逐漸加強(qiáng)[17-18],對(duì)土壤養(yǎng)分需求加大,特別是有機(jī)質(zhì)消耗過大.

    圖2 不同植被恢復(fù)年限土壤有機(jī)質(zhì)變化

    2.2.3 不同植被恢復(fù)年限土壤N、P、K 由圖3可知,全N、全P、全K均呈現(xiàn)出相同規(guī)律.不同植被恢復(fù)年限土壤N、P、K與未治理土壤相比顯著增加(3~4a),然后又顯著降低(<0.05)(6~10a).不同植被恢復(fù)年限的土壤N、P、K均顯著高于未治理土壤,但又顯著低于未開采對(duì)照土壤(<0.05).結(jié)果表明,土壤N、P、K隨著植被恢復(fù)年限的增加,并未恢復(fù)到自然水平狀態(tài).不同植被恢復(fù)年限稀土礦堆浸廢棄地土壤全N、全P、全K質(zhì)量分?jǐn)?shù)變化范圍分別是0.08~0.17、0.05~0.09、0.42~0.56g/kg,與福建省土壤平均N、P、K質(zhì)量分?jǐn)?shù)(分別是(0.136±0.03)、(0.108±0.184)、(2.29±0.33)g/kg)[19]相比,稀土礦堆浸廢棄地土壤全N、全P、全K質(zhì)量分?jǐn)?shù)顯著偏低.在植被恢復(fù)過程中,由于堆浸廢棄地有大量硫酸銨殘留,前期土壤全N含量比較高,后期由于植物難以吸收利用氮素,硫酸銨又易被雨水沖刷流失,土壤全N含量隨著植被恢復(fù)年限的加長(zhǎng)逐漸降低.南方亞熱帶土壤普遍都缺乏磷元素,對(duì)照樣地全磷含量?jī)H為0.13g/kg,隨著植被恢復(fù)年限的加長(zhǎng),土壤全磷含量呈現(xiàn)出下降趨勢(shì),植被恢復(fù)10a時(shí)僅為0.05g/kg,磷素缺乏導(dǎo)致微生物活動(dòng)受到限制,土壤酶分泌和土壤N礦化受影響,土壤N質(zhì)量分?jǐn)?shù)增加量較少.

    圖3 不同植被恢復(fù)年限土壤全量養(yǎng)分變化

    圖4 不同植被恢復(fù)年限土壤速效養(yǎng)分變化

    2.2.4 年限土壤速效養(yǎng)分 由圖4可知,速效氮、速效磷呈現(xiàn)出相同規(guī)律.不同植被恢復(fù)年限土壤速效氮、速效磷與未治理土壤相比顯著增加(3~4a),然后又顯著降低(6~10a)(<0.05).不同植被恢復(fù)年限土壤速效氮、速效磷均顯著高于未治理土壤,但又顯著低于未開采對(duì)照土壤(<0.05).植被恢復(fù)3a土壤速效鉀與未治理土壤相比顯著增加(<0.05),與未開采對(duì)照土壤速效鉀含量相比并無(wú)顯著差(>0.05);隨后不同植被恢復(fù)年限土壤速效鉀又顯著低于未開采對(duì)照土壤,但仍顯著高于未治理土壤(6~10a)(< 0.05).結(jié)果表明,隨著植被灰度年限的增加,土壤速效氮、速效磷并未恢復(fù)到自然狀態(tài)水平,速效鉀含量在恢復(fù)初期3a恢復(fù)到自然水平,但隨著植被恢復(fù)年限的增加,土壤速效鉀含量又顯著降低(<0.05).

    2.3 不同植被恢復(fù)年限土壤細(xì)菌結(jié)構(gòu)分析

    2.3.1 不同植被恢復(fù)年限土壤細(xì)菌測(cè)序數(shù)據(jù)和OTU聚類分析 由圖5可知,通過對(duì)土壤16S rRNA基因V4區(qū)進(jìn)行測(cè)序,除去低質(zhì)量、barcode和引物序列后,12個(gè)樣品得到的有效序列總數(shù)為766215條,對(duì)所有樣品有效序列進(jìn)行聚類,以97%的一致性將序列聚類成為OTUs.M9、M6、M3、CK包含的OTU數(shù)目分別為4383、3818、3399、3616,其中4組樣品共有的OTU數(shù)目為1683,特有的OTU數(shù)目分別為649、543、287、856,植被恢復(fù)年限越長(zhǎng),特有的OTU數(shù)目越多.

    圖5 OTU分布韋恩圖

    2.3.2 不同植被恢復(fù)年限土壤細(xì)菌多樣性比較 Chao 1指數(shù)與ACE指數(shù)主要反映生態(tài)系統(tǒng)中物種豐富度;Shannon指數(shù)主要反映群落物種多樣性; Simpson指數(shù)代表在一個(gè)群落中隨機(jī)抽取的兩個(gè)個(gè)體分屬于不同種的概率,該指數(shù)越接近1,表示該群落物種多樣性越大.由表2可知,通過篩選,各樣地中有效序列數(shù)為38670~77693條,其中植被恢復(fù)3a有效序列數(shù)最少,6a有效序列數(shù)最多.各樣地中分類單元OUT數(shù)為2761~3893,各樣地的覆蓋度指數(shù)范圍為0.938~0.964.隨著植被恢復(fù)年限的加長(zhǎng),Chao 1指數(shù)與ACE指數(shù)與對(duì)照相比較呈現(xiàn)出下降趨勢(shì); Shannon指數(shù)與Simpson指數(shù)同樣呈現(xiàn)出下降趨勢(shì),表明植被恢復(fù)年限越長(zhǎng),土壤細(xì)菌的豐富度和多樣性越低.隨著植被恢復(fù)年限的加長(zhǎng),Shannon指數(shù)、Simpson指數(shù)都呈現(xiàn)出下降趨勢(shì),表明植被恢復(fù)年限越長(zhǎng),土壤微生物豐富度和多樣性越低,與Zhou等[20]、Li等[21]研究結(jié)論不一致.可能的原因是:一方面,土壤微生物大多在pH值為中性的土壤中生長(zhǎng),而稀土礦堆浸廢棄地pH值在植被恢復(fù)10a間變化范圍為4.0~4.5,無(wú)疑會(huì)使土壤細(xì)菌的生物活性及豐富度降低[22].另一方面,稀土礦堆浸廢棄地在植被恢復(fù)后期,土壤養(yǎng)分消耗過大,出現(xiàn)植被退化甚至出現(xiàn)物種消失的情況[23],有機(jī)質(zhì)、全氮等含量下降對(duì)土壤微生物多樣性產(chǎn)生影響[24,18].

    表2 不同治理年限土壤細(xì)菌多樣性與豐富度

    2.3.3 不同植被恢復(fù)年限土壤細(xì)菌門和綱水平相對(duì)豐度 由圖6可知,對(duì)照樣地與各不同植被恢復(fù)年限土壤細(xì)菌群落組成相似,在門水平上的表現(xiàn)為β-變形菌門(Proteobacteria)、廣古菌門(Euryarchaeota)、酸桿菌門(Acidobacteria)、厚壁菌門(Firmicutes)、藍(lán)藻菌門(Cyanobacteria)、放線菌門(Actinobacteria)等.-變形菌門(Proteobacteria)是對(duì)照樣地以及各不同植被恢復(fù)年限土壤細(xì)菌中的第一大細(xì)菌類群,各組占比超過45%;廣古菌門(Euryarchaeota)是對(duì)照樣地中特有的菌門.隨著植被恢復(fù)時(shí)間的加長(zhǎng),酸桿菌門(Acidobacteria)相對(duì)豐度呈現(xiàn)出上升-下降的趨勢(shì),具體表現(xiàn)為CK

    2.3.4 不同植被恢復(fù)年限土壤優(yōu)勢(shì)細(xì)菌與土壤環(huán)境的冗余分析 由圖8可知,廣古菌門(Euryarchaeota)與土壤全P、pH夾角最小,相關(guān)性最大,但也與土壤容重、全N、全K、速效鉀、速效磷呈現(xiàn)正相關(guān).浮霉菌門(Planctomycetes)、綠彎菌門(Chloroflexi)、放線菌門(Actinobacteria)等菌門與土壤環(huán)境呈現(xiàn)出負(fù)相關(guān).厚壁菌門(Firmicutes)、-變形菌門(Proteobacteria)與土壤環(huán)境呈現(xiàn)正相關(guān).由圖9可知,屬()與羅爾斯通菌屬()分別與全P、pH夾角最小,相關(guān)性最大.芽孢桿菌屬()、甲基桿菌屬()與土壤容重、全氮存在正相關(guān),與全K、有效鉀、有效磷、全P、pH存在負(fù)相關(guān).屬水平上,廣古菌門(Euryarchaeota)中甲烷桿菌屬(Methanobacterium)與羅爾斯通菌屬(Ralstonia)與土壤全P、pH值相關(guān)性最大,與土壤容重、全N、全K等土壤理化性質(zhì)也存在一定相關(guān)性,與王友生等[26]、李兆龍[11]林海[32]結(jié)論相一致.

    圖6 不同植被恢復(fù)年限土壤細(xì)菌門水平相對(duì)豐度變化

    圖7 不同植被恢復(fù)年限土壤細(xì)菌綱水平相對(duì)豐度變化

    圖8 不同治理年限土壤門水平優(yōu)勢(shì)細(xì)菌與土壤環(huán)境的冗余分析

    ρb、TK、TN、TP、AK、AP、pH分別代表土壤容重、全鉀、全氮、全磷、速效鉀、速效磷、酸堿度

    圖9 不同治理年限土壤屬水平優(yōu)勢(shì)細(xì)菌與土壤環(huán)境的冗余分析

    3 結(jié)論

    3.1 長(zhǎng)汀稀土堆浸廢棄地不同植被恢復(fù)年限土壤理化性質(zhì)在恢復(fù)前期呈現(xiàn)出上升趨勢(shì)(3~4a),隨著植被恢復(fù)年限的增加,土壤理化性質(zhì)呈現(xiàn)出下降趨勢(shì)(6~10a).

    3.2 植被恢復(fù)年限越長(zhǎng),土壤微生物多樣性及相對(duì)豐度降低,如藍(lán)藻菌門(Cyanobacteria)、厚壁菌門(Firmicutes)等,但一些細(xì)菌群落如酸桿菌門(Acidobacteria)、-變形菌綱(Gammaproteobacteria)和一些未標(biāo)記細(xì)菌菌群相對(duì)豐度增加.

    3.3 冗余分析發(fā)現(xiàn),廣古菌門(Euryarchaeota)、厚壁菌門(Firmicutes)、Methanobacterium屬()與羅爾斯通菌屬()與土壤理化性質(zhì)呈現(xiàn)正相關(guān)系,對(duì)其相對(duì)豐度存在影響.

    3.4 在稀土礦廢棄地植被恢復(fù)一段時(shí)間后(3~4a)需要對(duì)其進(jìn)行人工撫育管理,防止土壤理化性質(zhì)降低,出現(xiàn)植被退化現(xiàn)象.

    [1] 劉文深,劉 暢,王志威,等.離子型稀土礦尾砂地植被恢復(fù)障礙因子研究 [J]. 土壤學(xué)報(bào), 2015,52(4):879-887.Liu W S, Liu C, Wang Z W, et al. Limiting factors for restoration of dumping sites of ionic rare earth mine tailings [J]. Acta Pedologica Sinica, 2015,52(4):879-887.

    [2] 王友生,侯曉龍,吳鵬飛,等.長(zhǎng)汀稀土礦廢棄地土壤重金屬污染特征及其評(píng)價(jià) [J]. 安全與環(huán)境學(xué)報(bào), 2014,14(4):259-262.Wang Y S, Hou X L, Wu P F, et al. Characteristics and evaluation of heavy metal pollution in abandoned soil of changting rare earth mine [J]. Journal of Safety and Environment, 2014,14(4):259-262.

    [3] 趙永紅,張 靜,周 丹,等.贛南某鎢礦區(qū)土壤重金屬污染狀況研究 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(8):2477-2484.Zhao Y H, Zhang J, Zhou D, et al. Study on soil heavy metal pollution in a tungsten mine area in southern jiangxi province [J]. China Environmental Science, 2015,35(8):2477-2484.

    [4] 簡(jiǎn)麗華.長(zhǎng)汀稀土廢礦區(qū)植被恢復(fù)與植被生態(tài)修復(fù)技術(shù) [J]. 現(xiàn)代農(nóng)業(yè)科技, 2012,(3):315-317. Jian L H. Vegetation Restoration and vegetation ecological restoration technology in changting rare earth waste mining area [J]. Modern Agricultural Sciences And Technology, 2012,(3):315-317.

    [5] 翁伯琦,王義祥,張偉利,等.強(qiáng)化防控與綜合植被恢復(fù)水土流失的思路及其對(duì)策研究 [J]. 福建農(nóng)業(yè)學(xué)報(bào), 2013,28(11):1063-1071. Weng B Q, Wang Y X, Zhang W L, et al. Thoughts on strengthening prevention and control and comprehensive vegetation restoration of soil and water loss and its countermeasures [J]. Fujian Journal of Agricultural Sciences, 2013,28(11):1063-1071.

    [6] 于君寶,王金達(dá),劉景雙,等.礦山復(fù)墾土壤營(yíng)養(yǎng)元素時(shí)空變化研究 [J]. 土壤學(xué)報(bào), 2002,(5):750-753. Yu J B, Wang J D, Liu J S, et al. Temporal and spatial variation of soil nutrient elements in mine reclamation [J]. Acta Pedologica Sinica, 2002,(5):750-753.

    [7] 郭李凱.不同復(fù)墾年限煤矸山重構(gòu)土壤養(yǎng)分和重金屬污染狀況研究 [D]. 臨汾:山西師范大學(xué), 2017. Guo L K. Study on soil nutrient and heavy metal pollution in coal gangue hills with different reclamation years [D]. Linfen: Shanxi Normal University, 2017.

    [8] 陜永杰,張美萍,白中科,等.平朔安太堡大型露天礦區(qū)土壤質(zhì)量演變過程分析 [J]. 干旱區(qū)研究, 2005,(4):149-152. Shan Y J, Zhang M P, Bai Z K, et al. Analysis of soil quality evolution process in pingshuo antaibao large open pit mining area [J]. Arid Zone Research, 2005,(4):149-152.

    [9] 王友生,吳鵬飛,侯曉龍,等.稀土礦廢棄地不同植被恢復(fù)模式對(duì)土壤肥力的影響 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2015,24(11):1831-1836. Wang Y S, Wu P F, Hou X L, et al. Effect of different revegetation model on soil properties in abandon mine area of rare earth [J]. Ecology and Environment, 2015,24(11):1831-1836.

    [10] 國(guó)家林業(yè)局.中華人民共和國(guó)林業(yè)行業(yè)標(biāo)準(zhǔn)-森林土壤分析方法 [M]. 北京:中國(guó)標(biāo)準(zhǔn)出版社, 2000.State Forestry Administration. Forestry industry standards of the people's republic of china - forest soil analysis methods [M]. Beijing: Chinese Standard Press, 2000.

    [11] 李兆龍.稀土礦場(chǎng)植物修復(fù)研究 [D]. 廣州:仲愷農(nóng)業(yè)工程學(xué)院, 2013. Li Z L. Studies on phytoremediation of rare-earth mind area [D]. Guangzhou: Zhongkai University of Agriculture, 2013.

    [12] 王 麗,夢(mèng) 麗,張金池,等.不同植被恢復(fù)模式下礦區(qū)廢棄地土壤水分物理性質(zhì)研究 [J]. 中國(guó)水土保持, 2010,(3):54-58. Wang L, Meng L, Zhang J C, et al. Mining area abandoned land under different vegetation restoration modes Study on physical properties of soil moisture [J]. Soil and Water Conservation in China, 2010,(3): 54-58.

    [13] 魏 婉.不同植被恢復(fù)模式對(duì)土壤質(zhì)量的影響 [D]. 南京:南京林業(yè)大學(xué), 2010. Wei W. Effect of different vegetation restoration models on soil quality [D]. Nanjing: Nanjing Forestry University, 2010.

    [14] 劉子壯,趙 晶,高照良,等.高速公路邊坡不同恢復(fù)年限土壤性質(zhì)及生態(tài)防護(hù)模式研究 [J]. 科學(xué)技術(shù)與工程, 2014,14(12):100-106. Liu Z Z, Zhao J, Gao Z L. Research on ecological protection modes and soil properties of different years of recovery along highway slopes [J]. Science Technology and Engineering, 2014,14(12):100-106.

    [15] 彭東海,侯曉龍,何宗明,等.金尾礦廢棄地不同植被恢復(fù)模式對(duì)土壤理化性質(zhì)的影響 [J]. 水土保持學(xué)報(bào), 2015,29(6):137-142. Peng D H, Hou X L, He Z M, et al. Effects of different vegetation restoration patterns on soil physical and chemical properties in wasteland of gold tailings [J]. Journal of Soil and Water Conservation, 2015,29(6):137-142.

    [16] 毛蘭花,查 軒,黃少燕,等.喬灌草植被恢復(fù)年限對(duì)紅壤區(qū)土壤養(yǎng)分的影響 [J]. 水土保持學(xué)報(bào), 2018,32(2):173-178. Mao L H, Zha X, Huang S Y, et al. Effects of management years of tress,shrubs and grasses on soil nutrients in red soil region [J]. Journal of Soil and Water Conservation, 2018,32(2):173-178.

    [17] 張秋芳,陳奶壽,陳 坦,等.不同恢復(fù)年限侵蝕紅壤生態(tài)化學(xué)計(jì)量特征 [J]. 中國(guó)水土保持科學(xué), 2016,14(2):59-66. Zhang Q F, Chen N S, Chen T, et al. Ecological stoichiometry characteristics of eroded red soil in different restoration years [J]. Science of Soil and Water Conservation, 2016,14(2):59-66.

    [18] 張 杰,余 潮,王自海,等.不同植被群落表層土壤中細(xì)菌群落多樣性 [J]. 環(huán)境科學(xué)研究, 2013,6(8):866-872. Zhang J, Yu C, Wang Z H, et al. Diversity of bacterial communities in surface soils under different vegetation communities [J]. Research of Environmental Sciences, 2013,26(8):866-872.

    [19] 福建省土壤普查辦公室.福建土壤 [M]. 福州:福建科學(xué)技術(shù)出版社, 1991:225-255. Fujian soil census office. fujian soil [M]. Fuzhou: Fujian Science and Technology Press, 1991:225-255.

    [20] Zhou J, Xia B, Treves D S, et al. Spatial and resource factors influencing high microbial diversity in soil [J]. Appl Environ Microbiol. 2002,68(1):326-334.

    [21] Li X F, Huang L B, Philip L BOND, et al. Bacterial diversity in response to direct revegetation in the Pb-Zn-Cu tailings under subtropical and semi-arid conditions [J]. Ecological Engineering, 2014,68(7):233-240.

    [22] 葉 雯,李永春,喻衛(wèi)武,等.不同種植年限香榧根際土壤微生物多樣性[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2018,29(11):3783-3792. Ye Wen, Li Y C, Yu W W, et al. Microbial biodiversity in rhizospheric soil of Torreya grandis 'Merrillii' relative to cultivation history [J]. Chines Journal of Applied Ecology, 2018,29(11):3783- 3792.

    [23] 王友生.稀土開采對(duì)紅壤生態(tài)系統(tǒng)的影響及其廢棄地植被恢復(fù)機(jī)理研究 [D]. 福州:福建農(nóng)林大學(xué), 2016.Wang Y S. Effects of Rare Earth Exploitation on Red Soil Ecosystem and Vegetation Restoration Mechanism of Waste Land [D]. Fuzhou: Fujian Agriculture and Forestry University, 2016.

    [24] 宋蒙亞,李忠佩,劉 明,等.不同林地凋落物組合對(duì)土壤速效養(yǎng)分和微生物群落功能多樣性的影響 [J]. 生態(tài)學(xué)雜志, 2014,33(9): 2454-2461. Song M Y, LI Z P, Liu M, et al. Effects of mixture of forest litter on nutrient contents and functional diversity of microbial community in soil [J]. Chinese Journal of Ecology, 2014,33(9):2454-2461.

    [25] 張紅娟,孫 燕,李 治,等.金堆城鉬礦尾礦堆積區(qū)土壤細(xì)菌多樣性分析 [J]. 生命科學(xué)研究, 2010,14(6):499-506. Zhang H J, Shun Y, Li Z, et al. Diversity of soil bacteria communities in areas of jinduicheng mo tailings [J]. Life Science Research, 2010,14(6):499-506.

    [26] 王友生,侯曉龍,蔡麗平,等.稀土開采對(duì)土壤細(xì)菌群落組成與多樣性的影響 [J]. 中國(guó)環(huán)境科學(xué), 2017,37(8):3089-3095. Wang Y S, Hou X L, Cai L P, et al. Impacts of rare earth mining on soil bacterial community composition and biodiversity [J]. China Environmental Science, 2017,37(8):3089-3095.

    [27] 王新洲,胡忠良,杜有新,等.喀斯特生態(tài)系統(tǒng)中喬木和灌木林根際土壤微生物生物量及其多樣性的比較 [J]. 土壤, 2010,42(2):224-229. Wang X Z, Hu Z L, Du Y X, et al. Comparison of microbial biomass and community structure of rhizosphere soil between forest and shrubbery in karst ecosystems [J]. Soil, 2010,42(2):224-229.

    [28] 張 倩,姚寶輝,王 纏,等.不同坡向下高寒草甸土壤理化特性和微生物數(shù)量特征 [J]. 生態(tài)學(xué)報(bào), 2019,(9):1-7. Zhang Q, Yao B H, Wang C, et al. Soil physical and chemical properties and microbial quantitative characteristics of alpine meadows with different slopes [J]. Acta Ecologica Sinica. 2019, (9):1-7.

    [29] 李建凱,李 珊,魏云林,等.云南文山鋁土礦區(qū)微生物多樣性初步研究 [J]. 中國(guó)微生態(tài)學(xué)雜志, 2015,27(11):1258-1261+1265. Li J K, Li S, W Y L, et al. Preliminary study on microbial diversity in Wenshan bauxite mining area, Yunnan [J]. China Journal of Microecology, 2015,27(11):1258-1261+1265.

    [30] 王光華,劉俊杰,于鎮(zhèn)華,等.土壤酸桿菌門細(xì)菌生態(tài)學(xué)研究進(jìn)展 [J]. 生物技術(shù)通報(bào), 2016,32(2):14-20. Wang G H, Liu J J, Yu Z H, et al. Advances in studies on the ecology of acid-bacillus bacteria [J]. Biotecnology Bulletin, 2016,32(2):14-20.

    [31] Smit E, Leeflang P, Gommans S, et al. Diversity and seasonal fluctuations of the dominant members of the bacterial soil community in a wheat field as determined by cultivation and molecular methods [J]. Appl Environ Microbiol., 2001,67(5):2284-2291.

    [32] 林 海,崔 軒,董穎博,等.銅尾礦庫(kù)重金屬Cu、Zn對(duì)細(xì)菌群落結(jié)構(gòu)的影響 [J]. 中國(guó)環(huán)境學(xué), 2014,34(12):3182-3188. Lin H, Cui X, Dong Y B, et al. Effects of heavy metals Cu and Zn on bacterial community structure in copper tailings reservoir [J]. China Environmental Science, 2014,34(12):3182-3188.

    Soil physicochemical characteristics and microbial evolution during vegetation restoration in ionic rare earth ore heap leaching waste land.

    LI Qi-yan1, WENG Bing-lin1, LI Zong-xun1, ZhAO Ya-man1, CHEN Shun-yu1, HOU Xiao-long1,2,3*, CAI Li-ping1,2,3

    (1.College of Forestry, Fujian Agriculture and Forestry University, Fuzhou 350002, China;2.Key Laboratory of State Forestry and Grassland Administration for Soil and Water Conservation in Red Soil Region of South China, Fuzhou 350002, China;3.Collaborative Innovation Center for Soil and Water Conservation in the Red Soil Region on Both Sides of the Taiwan Straits, Fuzhou 350002,)., 2019,39(10):4360~4368

    In order to mitigate of the sever soil erosion and degradation in ecological environment in the rare earth mining area, a vegetation restoration project in the abandoned land of rare earth mining area in Changting, Fujian was implemented several years ago. Physical and chemical properties and microorganisms’ composition of soil in the heap leaching of abandoned land of rare earth mining, Changting after different years of vegetation restoration were analyzed. The results showed that physical and chemical properties of soil after three and four years of vegetation restoration were significantly higher than those in the untreated abandoned land (without restoration) but lower than the unmined area (< 0.05). After 4years of vegetation restoration, the physical and chemical properties of soil in the abandoned land of rare earth mining area, Changting showed a downward trend, but still significantly higher than that of untreated abandoned land (< 0.05). Differences were identified in the abundance and diversity of soil bacteria of abandoned land of rare earth mining between different years of vegetation restoration. The diversity and abundance of soil bacteria decreased with the years of vegetation restoration,. However, the relative abundance of some special bacterial groups such as Gamma proteobacteria, Alphaproteobacteria, Acidobacteriaand Cyanobacteria significantly increased in the stands of different years of vegetation restoration. The relative abundance of Cyanobacteria, a type of bacteria of nitrogen fixing, decreased along with the process of vegetation recovery. Euryarchaeota only existed in the unmined land, and Acidobacteria increased with the increase of years of vegetation restoration. The dominant flora in soil, Euryarchaeota, Firmicutes, Proteobacteria,, were positively correlated with soil pH, total phosphorus, total potassium and available phosphorus. The diversity and relative abundance of bacterial flora such as Euryarchaeota and Firmicutes were significantly affected by the physical and chemical properties of soil. In conclusion, the artificial prevention and control measures should be taken after 3-4 years of vegetation restoration in the heap leaching land of ion-type rare earth mining area to mitigate the degradation.

    ionic rare earth ore;heap leaching wasteland;microbial evolution;16S rRNA;redundancy analysis

    X53;X172

    A

    1000-6923(2019)10-4360-09

    李啟艷(1994-),女,貴州貴陽(yáng)人,福建農(nóng)林大學(xué)碩士研究生,主要從事恢復(fù)生態(tài)與生態(tài)工程研究.

    2019-03-26

    福建省科技廳重點(diǎn)項(xiàng)目(2017Y001);福建省環(huán)??萍加?jì)劃項(xiàng)目(2018R012);福建農(nóng)林大學(xué)科技創(chuàng)新專項(xiàng)基金項(xiàng)目(CXZX2018126); (CXZX2018137)

    * 責(zé)任作者, 副教授, xl.hou@fafu.edu.cn

    猜你喜歡
    廢棄地稀土礦菌門
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    金屬礦山廢棄地分析及植物修復(fù)技術(shù)研究
    饑餓與重?cái)z食對(duì)河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價(jià)值分析
    贛南離子吸附型稀土礦的發(fā)現(xiàn)與勘查開發(fā)研究
    礦山廢棄地分析及生態(tài)環(huán)境修復(fù)技術(shù)研究進(jìn)展
    基于墑權(quán)的屬性識(shí)別模型在采煤廢棄地肥力評(píng)價(jià)中的應(yīng)用
    河南發(fā)現(xiàn)中型規(guī)模稀土礦
    西部資源(2015年3期)2015-08-15 00:46:57
    冶山鐵礦區(qū)部分廢棄地整治方案研究
    金屬礦山(2013年11期)2013-03-11 16:55:14
    激情视频va一区二区三区| 国产成人av激情在线播放| 嫁个100分男人电影在线观看| 国产精品九九99| 久久精品人人爽人人爽视色| 欧美日本亚洲视频在线播放| 99精品久久久久人妻精品| 国产三级黄色录像| 黄色女人牲交| 国产高清激情床上av| av天堂在线播放| 一级黄色大片毛片| 欧美 亚洲 国产 日韩一| 亚洲午夜精品一区,二区,三区| 国产在线精品亚洲第一网站| 精品日产1卡2卡| 极品人妻少妇av视频| 深夜精品福利| 美女高潮到喷水免费观看| 两个人免费观看高清视频| 波多野结衣av一区二区av| 悠悠久久av| 亚洲专区字幕在线| 欧美国产精品va在线观看不卡| 可以在线观看毛片的网站| 亚洲精品国产色婷婷电影| 精品久久蜜臀av无| 真人一进一出gif抽搐免费| 午夜精品久久久久久毛片777| 麻豆av在线久日| 亚洲欧美一区二区三区黑人| 琪琪午夜伦伦电影理论片6080| 国产av一区二区精品久久| 深夜精品福利| 在线永久观看黄色视频| 男女下面进入的视频免费午夜 | 99久久久亚洲精品蜜臀av| 水蜜桃什么品种好| 99久久99久久久精品蜜桃| 久久人人爽av亚洲精品天堂| 人妻丰满熟妇av一区二区三区| 亚洲国产中文字幕在线视频| 在线av久久热| 无人区码免费观看不卡| 欧美成狂野欧美在线观看| 淫秽高清视频在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品日韩av在线免费观看 | 村上凉子中文字幕在线| 国产一区二区三区视频了| 纯流量卡能插随身wifi吗| 午夜视频精品福利| 国产无遮挡羞羞视频在线观看| 欧美日韩精品网址| 精品久久久精品久久久| 日韩精品中文字幕看吧| 精品久久久精品久久久| 韩国av一区二区三区四区| 一级毛片精品| 国产精品永久免费网站| av国产精品久久久久影院| 最近最新中文字幕大全电影3 | 91精品国产国语对白视频| 亚洲精品久久午夜乱码| 黄色丝袜av网址大全| 中文字幕色久视频| 性少妇av在线| 无人区码免费观看不卡| 亚洲五月婷婷丁香| 久久香蕉国产精品| 亚洲精品久久成人aⅴ小说| 久久国产乱子伦精品免费另类| 黑人巨大精品欧美一区二区蜜桃| 18禁观看日本| 国产成人欧美| 夜夜看夜夜爽夜夜摸 | 三上悠亚av全集在线观看| 少妇粗大呻吟视频| 在线观看www视频免费| 黑人巨大精品欧美一区二区mp4| 久久久久国产精品人妻aⅴ院| 欧美成人午夜精品| 精品少妇一区二区三区视频日本电影| 国产精品一区二区精品视频观看| 老司机亚洲免费影院| 国产伦一二天堂av在线观看| 亚洲av日韩精品久久久久久密| 长腿黑丝高跟| www.999成人在线观看| 色在线成人网| 老司机亚洲免费影院| 精品国产国语对白av| 日日夜夜操网爽| 久久影院123| 天堂影院成人在线观看| 一级片'在线观看视频| 精品日产1卡2卡| 狂野欧美激情性xxxx| 久久精品国产99精品国产亚洲性色 | 国产精品美女特级片免费视频播放器 | av网站免费在线观看视频| 亚洲一码二码三码区别大吗| av超薄肉色丝袜交足视频| 欧美人与性动交α欧美精品济南到| 欧美日韩视频精品一区| tocl精华| 麻豆国产av国片精品| 一区在线观看完整版| 在线天堂中文资源库| cao死你这个sao货| 夜夜爽天天搞| 久久99一区二区三区| 久久久久久久精品吃奶| 电影成人av| 亚洲狠狠婷婷综合久久图片| 国产精品久久视频播放| 免费人成视频x8x8入口观看| 俄罗斯特黄特色一大片| 女生性感内裤真人,穿戴方法视频| 99国产精品免费福利视频| 长腿黑丝高跟| 精品国产美女av久久久久小说| 国产一区二区在线av高清观看| 美女高潮到喷水免费观看| 午夜福利免费观看在线| 露出奶头的视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利在线观看吧| 亚洲黑人精品在线| 黄色片一级片一级黄色片| 很黄的视频免费| 精品无人区乱码1区二区| 日韩一卡2卡3卡4卡2021年| 国产精品亚洲av一区麻豆| 久99久视频精品免费| 久久久久国产一级毛片高清牌| 欧美成人免费av一区二区三区| 亚洲成国产人片在线观看| 热99re8久久精品国产| 波多野结衣一区麻豆| 夜夜夜夜夜久久久久| 可以在线观看毛片的网站| 91精品国产国语对白视频| 精品卡一卡二卡四卡免费| 国产精品一区二区精品视频观看| 亚洲五月色婷婷综合| 亚洲欧美日韩无卡精品| 另类亚洲欧美激情| 日本精品一区二区三区蜜桃| 午夜精品国产一区二区电影| 亚洲 国产 在线| a级毛片黄视频| 激情在线观看视频在线高清| 亚洲va日本ⅴa欧美va伊人久久| 久久精品亚洲熟妇少妇任你| 亚洲午夜理论影院| 国产熟女xx| 色老头精品视频在线观看| 欧美亚洲日本最大视频资源| 美女午夜性视频免费| av中文乱码字幕在线| 一区二区日韩欧美中文字幕| 久久国产精品男人的天堂亚洲| 久久狼人影院| 久久精品国产亚洲av高清一级| 老司机午夜福利在线观看视频| 国产人伦9x9x在线观看| 精品第一国产精品| 国产人伦9x9x在线观看| 成人18禁高潮啪啪吃奶动态图| 多毛熟女@视频| 欧美人与性动交α欧美精品济南到| 后天国语完整版免费观看| 欧美成狂野欧美在线观看| 一边摸一边抽搐一进一小说| 午夜福利在线观看吧| 亚洲黑人精品在线| 后天国语完整版免费观看| av天堂久久9| 美女国产高潮福利片在线看| 久久久久久久久中文| 成人国产一区最新在线观看| 欧美在线一区亚洲| 9色porny在线观看| 久久国产乱子伦精品免费另类| 国产av精品麻豆| 亚洲一码二码三码区别大吗| 亚洲精品久久午夜乱码| 波多野结衣高清无吗| 久久精品影院6| 夜夜夜夜夜久久久久| 国产激情欧美一区二区| 777久久人妻少妇嫩草av网站| 欧美黑人欧美精品刺激| 制服诱惑二区| 欧美日韩福利视频一区二区| 成人手机av| svipshipincom国产片| 国产精品99久久99久久久不卡| 神马国产精品三级电影在线观看 | 天堂动漫精品| 午夜精品久久久久久毛片777| 久久人人爽av亚洲精品天堂| 国产精品1区2区在线观看.| av视频免费观看在线观看| 国产精品av久久久久免费| 免费观看人在逋| 欧美在线一区亚洲| 欧美丝袜亚洲另类 | 黑人巨大精品欧美一区二区mp4| 首页视频小说图片口味搜索| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品一区二区三区在线| 国产成人免费无遮挡视频| 国产亚洲欧美精品永久| x7x7x7水蜜桃| 亚洲黑人精品在线| 国内毛片毛片毛片毛片毛片| 午夜福利欧美成人| 在线十欧美十亚洲十日本专区| av电影中文网址| 亚洲 欧美一区二区三区| 国产精品一区二区在线不卡| 大码成人一级视频| 极品教师在线免费播放| 国产黄a三级三级三级人| 久久中文字幕一级| 日韩人妻精品一区2区三区| 亚洲aⅴ乱码一区二区在线播放 | 黄色a级毛片大全视频| 50天的宝宝边吃奶边哭怎么回事| 超碰成人久久| 欧美最黄视频在线播放免费 | av网站免费在线观看视频| 两性夫妻黄色片| 操美女的视频在线观看| 国产精品亚洲av一区麻豆| 淫妇啪啪啪对白视频| 桃色一区二区三区在线观看| 国产成人免费无遮挡视频| 亚洲av成人不卡在线观看播放网| 免费少妇av软件| 久久亚洲真实| 日本三级黄在线观看| 麻豆成人av在线观看| 一区二区日韩欧美中文字幕| a级毛片在线看网站| 级片在线观看| 91在线观看av| 精品一区二区三卡| 亚洲人成电影观看| 久久久久久久久免费视频了| 亚洲一区二区三区不卡视频| 美女午夜性视频免费| 人成视频在线观看免费观看| 亚洲精品国产色婷婷电影| 中出人妻视频一区二区| 色综合婷婷激情| 最好的美女福利视频网| 亚洲av成人av| 国产三级黄色录像| 人成视频在线观看免费观看| 91国产中文字幕| 黄色毛片三级朝国网站| 美女国产高潮福利片在线看| 99精品欧美一区二区三区四区| 欧美日韩视频精品一区| 国产99久久九九免费精品| 91字幕亚洲| 久久久国产精品麻豆| 国产欧美日韩一区二区精品| 高清欧美精品videossex| 久久精品aⅴ一区二区三区四区| 国产黄a三级三级三级人| 叶爱在线成人免费视频播放| 久久天躁狠狠躁夜夜2o2o| 视频区图区小说| 两性夫妻黄色片| 精品日产1卡2卡| 国产高清国产精品国产三级| 亚洲国产精品一区二区三区在线| 久久久国产精品麻豆| 国产精品免费视频内射| 999精品在线视频| 男人舔女人的私密视频| 欧美在线黄色| 在线天堂中文资源库| 黄片大片在线免费观看| 很黄的视频免费| 日韩有码中文字幕| 女人精品久久久久毛片| 午夜两性在线视频| 在线观看一区二区三区| 中文字幕人妻丝袜制服| 亚洲精品国产色婷婷电影| 免费在线观看影片大全网站| 亚洲精品美女久久久久99蜜臀| 亚洲国产中文字幕在线视频| 999久久久精品免费观看国产| 亚洲av成人av| 最新在线观看一区二区三区| 色老头精品视频在线观看| 亚洲伊人色综图| 亚洲av日韩精品久久久久久密| 老司机午夜福利在线观看视频| 一二三四在线观看免费中文在| 国产精品综合久久久久久久免费 | 国产免费男女视频| 国产免费av片在线观看野外av| 大陆偷拍与自拍| 国产精品野战在线观看 | 免费高清在线观看日韩| 精品午夜福利视频在线观看一区| 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| av天堂久久9| 99久久久亚洲精品蜜臀av| 亚洲人成77777在线视频| 黄色a级毛片大全视频| 新久久久久国产一级毛片| av网站免费在线观看视频| 亚洲在线自拍视频| 久久久国产欧美日韩av| x7x7x7水蜜桃| 久热爱精品视频在线9| 99精品欧美一区二区三区四区| 亚洲黑人精品在线| 精品久久久久久久久久免费视频 | 免费在线观看完整版高清| av福利片在线| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品久久久久5区| www.熟女人妻精品国产| 每晚都被弄得嗷嗷叫到高潮| 高清黄色对白视频在线免费看| 成人影院久久| 亚洲视频免费观看视频| 国产国语露脸激情在线看| 欧美乱色亚洲激情| 天堂√8在线中文| 成人手机av| 青草久久国产| 日韩欧美一区视频在线观看| 亚洲熟妇熟女久久| 日日干狠狠操夜夜爽| 18禁裸乳无遮挡免费网站照片 | 精品乱码久久久久久99久播| 人成视频在线观看免费观看| 91国产中文字幕| 国产高清videossex| 亚洲av五月六月丁香网| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区三| 免费高清在线观看日韩| 别揉我奶头~嗯~啊~动态视频| 老司机午夜十八禁免费视频| 午夜老司机福利片| 久久精品91无色码中文字幕| 免费人成视频x8x8入口观看| 99精品欧美一区二区三区四区| 俄罗斯特黄特色一大片| 精品国产一区二区久久| 最新美女视频免费是黄的| 国产精品99久久99久久久不卡| 亚洲成人免费电影在线观看| 亚洲国产精品sss在线观看 | 制服诱惑二区| 国产人伦9x9x在线观看| 国产精品亚洲一级av第二区| 亚洲少妇的诱惑av| 悠悠久久av| 亚洲五月色婷婷综合| 啦啦啦免费观看视频1| 国产91精品成人一区二区三区| 人人妻,人人澡人人爽秒播| 在线免费观看的www视频| 视频区欧美日本亚洲| 国产成年人精品一区二区 | 很黄的视频免费| 久久婷婷成人综合色麻豆| 午夜精品在线福利| 丝袜在线中文字幕| 91国产中文字幕| 久久青草综合色| 亚洲国产毛片av蜜桃av| 免费高清视频大片| 精品福利永久在线观看| 真人做人爱边吃奶动态| 亚洲自拍偷在线| 欧美成狂野欧美在线观看| 欧美激情久久久久久爽电影 | 婷婷六月久久综合丁香| 久久性视频一级片| 亚洲avbb在线观看| 91成人精品电影| 国产单亲对白刺激| 人人妻,人人澡人人爽秒播| 嫩草影院精品99| 午夜老司机福利片| 国产精品1区2区在线观看.| 在线观看免费视频网站a站| 丁香欧美五月| 久9热在线精品视频| 丰满迷人的少妇在线观看| av电影中文网址| 国产不卡一卡二| 婷婷丁香在线五月| 精品电影一区二区在线| 又大又爽又粗| 美女国产高潮福利片在线看| 欧美日本亚洲视频在线播放| 黄色女人牲交| 亚洲国产精品一区二区三区在线| 一级毛片高清免费大全| 国产成人系列免费观看| 后天国语完整版免费观看| 国产乱人伦免费视频| 叶爱在线成人免费视频播放| 亚洲精品av麻豆狂野| 在线观看www视频免费| 亚洲美女黄片视频| av福利片在线| 激情在线观看视频在线高清| 少妇裸体淫交视频免费看高清 | 国产精品98久久久久久宅男小说| 久久久久久大精品| 亚洲欧美精品综合久久99| 色老头精品视频在线观看| 精品人妻在线不人妻| 一进一出抽搐动态| 丰满饥渴人妻一区二区三| 免费女性裸体啪啪无遮挡网站| 人妻丰满熟妇av一区二区三区| 美女国产高潮福利片在线看| 国产三级黄色录像| 欧美成狂野欧美在线观看| 欧美人与性动交α欧美软件| 欧美国产精品va在线观看不卡| 国产不卡一卡二| 午夜精品久久久久久毛片777| 日日爽夜夜爽网站| 国产免费av片在线观看野外av| 欧美黄色淫秽网站| 亚洲一区二区三区色噜噜 | 日韩欧美一区二区三区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 黄色片一级片一级黄色片| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出 | 纯流量卡能插随身wifi吗| 黑人巨大精品欧美一区二区蜜桃| 夜夜躁狠狠躁天天躁| netflix在线观看网站| 亚洲欧美激情在线| 成人特级黄色片久久久久久久| 天堂俺去俺来也www色官网| 在线观看一区二区三区| 日韩欧美免费精品| 久久伊人香网站| 黄片播放在线免费| 人人妻人人爽人人添夜夜欢视频| 国产一卡二卡三卡精品| 如日韩欧美国产精品一区二区三区| 国产免费男女视频| 日韩欧美三级三区| 法律面前人人平等表现在哪些方面| 亚洲精品国产区一区二| 女生性感内裤真人,穿戴方法视频| 美女 人体艺术 gogo| 韩国av一区二区三区四区| 欧美黑人欧美精品刺激| 无人区码免费观看不卡| 欧美国产精品va在线观看不卡| 另类亚洲欧美激情| 一级作爱视频免费观看| 国产精品99久久99久久久不卡| 岛国视频午夜一区免费看| 三级毛片av免费| 曰老女人黄片| 18禁黄网站禁片午夜丰满| 国产99白浆流出| 久久人妻熟女aⅴ| 亚洲av成人一区二区三| 在线观看一区二区三区| 香蕉国产在线看| 黑人巨大精品欧美一区二区mp4| 亚洲男人天堂网一区| 女人精品久久久久毛片| 亚洲精品一卡2卡三卡4卡5卡| 香蕉丝袜av| 久久人妻av系列| 女人高潮潮喷娇喘18禁视频| 亚洲性夜色夜夜综合| 人妻久久中文字幕网| 亚洲av五月六月丁香网| 国产欧美日韩精品亚洲av| 十分钟在线观看高清视频www| 少妇粗大呻吟视频| 男女下面插进去视频免费观看| 亚洲精品国产色婷婷电影| 久久久久久人人人人人| 国产成人精品无人区| 丰满迷人的少妇在线观看| 波多野结衣高清无吗| 丝袜人妻中文字幕| 天堂俺去俺来也www色官网| 国产成年人精品一区二区 | 看黄色毛片网站| 精品乱码久久久久久99久播| 91九色精品人成在线观看| 成人亚洲精品一区在线观看| 热re99久久国产66热| 久久中文看片网| 热re99久久国产66热| 国产一区二区激情短视频| 日日夜夜操网爽| 久久人妻熟女aⅴ| 一本综合久久免费| 五月开心婷婷网| 99国产精品99久久久久| 夜夜爽天天搞| 精品福利观看| 18美女黄网站色大片免费观看| 午夜视频精品福利| av福利片在线| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩视频精品一区| 亚洲人成网站在线播放欧美日韩| 电影成人av| 少妇 在线观看| 日本wwww免费看| 精品人妻1区二区| 色婷婷av一区二区三区视频| 日韩欧美一区二区三区在线观看| 天堂影院成人在线观看| 老汉色∧v一级毛片| 女人被狂操c到高潮| 美女大奶头视频| 色综合欧美亚洲国产小说| 国产欧美日韩精品亚洲av| 啦啦啦在线免费观看视频4| 国产一区二区激情短视频| 老熟妇乱子伦视频在线观看| cao死你这个sao货| 9热在线视频观看99| 国产成人精品无人区| 久久久久久人人人人人| 午夜老司机福利片| 又紧又爽又黄一区二区| 日本撒尿小便嘘嘘汇集6| 午夜免费成人在线视频| 久久性视频一级片| 亚洲第一av免费看| 91麻豆av在线| 天天影视国产精品| 中文字幕av电影在线播放| 欧美激情极品国产一区二区三区| 国产精品电影一区二区三区| 一个人观看的视频www高清免费观看 | 国产成人精品在线电影| 色播在线永久视频| 黄色视频不卡| 黄色成人免费大全| 夜夜爽天天搞| 精品高清国产在线一区| 日韩大码丰满熟妇| 亚洲精品成人av观看孕妇| a在线观看视频网站| 久久久国产欧美日韩av| 成人亚洲精品av一区二区 | 国产色视频综合| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区精品| 日韩精品免费视频一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 欧美老熟妇乱子伦牲交| 91在线观看av| 中文字幕另类日韩欧美亚洲嫩草| 高清欧美精品videossex| 极品教师在线免费播放| 不卡一级毛片| 国产欧美日韩一区二区精品| 午夜视频精品福利| 丰满人妻熟妇乱又伦精品不卡| 亚洲激情在线av| 最近最新中文字幕大全电影3 | 亚洲免费av在线视频| 88av欧美| 亚洲 欧美一区二区三区| 国产一区二区激情短视频| 国产精品一区二区精品视频观看| 1024视频免费在线观看| 黄网站色视频无遮挡免费观看| 色精品久久人妻99蜜桃| 亚洲精品中文字幕在线视频| 日韩有码中文字幕| 嫩草影视91久久| 黑人欧美特级aaaaaa片| 美女高潮到喷水免费观看| 人妻丰满熟妇av一区二区三区| 欧美日韩亚洲综合一区二区三区_| 午夜91福利影院| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 免费看十八禁软件| 欧美成人午夜精品| 国产单亲对白刺激| 亚洲精品在线观看二区| 国产亚洲欧美精品永久| 亚洲精品一二三| 精品第一国产精品| 国产一区二区在线av高清观看| 欧美乱色亚洲激情| 亚洲国产欧美一区二区综合| 亚洲自拍偷在线|