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

    不同土地利用方式對(duì)三峽庫區(qū)消落帶土壤細(xì)菌和真菌多樣性的影響

    2017-06-22 13:47:48李昌曉任慶水
    生態(tài)學(xué)報(bào) 2017年10期
    關(guān)鍵詞:三峽庫區(qū)林地土地利用

    秦 紅,李昌曉,任慶水

    西南大學(xué)生命科學(xué)學(xué)院,三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,重慶市三峽庫區(qū)植物生態(tài)與資源重點(diǎn)實(shí)驗(yàn)室,重慶 400715

    ?

    不同土地利用方式對(duì)三峽庫區(qū)消落帶土壤細(xì)菌和真菌多樣性的影響

    秦 紅,李昌曉*,任慶水

    西南大學(xué)生命科學(xué)學(xué)院,三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,重慶市三峽庫區(qū)植物生態(tài)與資源重點(diǎn)實(shí)驗(yàn)室,重慶 400715

    旨在為三峽消落帶選擇適合的生態(tài)恢復(fù)方式提供參考依據(jù)。以三峽庫區(qū)重慶忠縣汝溪河流域典型消落帶為研究區(qū)域,于2015年6月進(jìn)行原位取樣,利用末端限制性片段長度多態(tài)性(Terminal Restriction Fragment Length Polymorphism,T-RFLP)方法,對(duì)消落帶的耕地、林地、棄耕地土壤細(xì)菌和真菌群落多樣性進(jìn)行研究。研究發(fā)現(xiàn):(1)除容重和密度外,不同用地類型對(duì)土壤各理化特性均產(chǎn)生顯著影響,林地的含水量、有機(jī)質(zhì)(OM)、全氮(TN)、全磷(TP)、速效氮(AN)、速效鉀(AK)、速效磷(AP)的含量顯著高于耕地和棄耕地(P<0.05)。(2)在三峽庫區(qū)消落帶不同土地利用方式下細(xì)菌和真菌多樣性均有顯著性差異。(3)耕地和林地的細(xì)菌多樣性無顯著性差異,均顯著高于棄耕地;不同用地方式中,土壤全磷(TP)和速效磷(AP)顯著影響細(xì)菌群落結(jié)構(gòu)(P<0.05)。(4)真菌群落的Shannon-Weiner指數(shù)和辛普森指數(shù)在耕地中均為最低,但均勻度指數(shù)在3種用地類型之間沒有顯著性差異;在不同用地方式的土壤中,土壤有機(jī)質(zhì)(OM)、全氮(TN)含量和pH顯著影響真菌群落結(jié)構(gòu)(P<0.05)。結(jié)論:與棄耕地和耕地相比,林地可固持、滯留和保有更多的土壤養(yǎng)分,有較高的細(xì)菌多樣性和真菌多樣性。在三峽消落帶165—175 m海拔高程,林地為最適宜的用地方式,建議限制耕作,推廣人工生態(tài)修復(fù)林地建設(shè)。

    三峽庫區(qū);消落帶;土地利用方式;細(xì)菌和真菌多樣性

    三峽大壩建成后,由于“冬蓄夏排”的水文調(diào)節(jié)方式,庫區(qū)形成了面積約為350 km2的消落帶。該區(qū)域是陸生環(huán)境和水生環(huán)境之間的過渡帶,具有過濾和緩沖的作用,在維護(hù)三峽水庫生態(tài)健康與生態(tài)安全方面具有極其特殊的地位[1]。研究表明,該地區(qū)是一個(gè)高度脆弱的生態(tài)系統(tǒng)[2],并且消落帶生態(tài)風(fēng)險(xiǎn)等級(jí)隨著海拔高度的升高而增大[3]。消落帶現(xiàn)主要面臨著土壤侵蝕[4]、生物多樣性下降[5]、富營養(yǎng)化污染[6]等的生態(tài)問題。為保護(hù)其生態(tài)環(huán)境,提高消落帶的生物多樣性和抗侵蝕能力,政府在消落帶推行退耕政策以及人工植被修復(fù)措施[7],使消落帶產(chǎn)生大量的植被修復(fù)地。同時(shí),由于庫區(qū)人口密度大,人地矛盾尖銳,加之水文節(jié)律符合農(nóng)作物生長物候,居住于消落帶上部的村民仍在消落帶中進(jìn)行傳統(tǒng)農(nóng)事耕作,農(nóng)耕地依舊大量存在[8]。因此,三峽庫區(qū)消落帶當(dāng)前仍存在多種用地方式并存的格局。

    不同的土地利用方式不僅直接改變土壤的理化性質(zhì),還改變著生態(tài)系統(tǒng)的結(jié)構(gòu)和功能,不合理的土地利用方式給土地生態(tài)環(huán)境帶來風(fēng)險(xiǎn)。土壤生態(tài)系統(tǒng)除含土壤顆粒和土壤養(yǎng)分外,微生物也占有重要地位,它能夠調(diào)節(jié)土壤營養(yǎng)物質(zhì)生物地球化學(xué)循環(huán)[9],影響地表植被生長[10]和生態(tài)系統(tǒng)的穩(wěn)定[11]等。微生物生理結(jié)構(gòu)簡(jiǎn)單,生活史較短,容易受環(huán)境條件的影響而發(fā)生變異,所以土壤微生物的多樣性及變異性常常用于預(yù)測(cè)土壤養(yǎng)分和環(huán)境質(zhì)量的變化[12]。在陸生系統(tǒng)中,生產(chǎn)者提供有機(jī)物質(zhì)和代謝物作為細(xì)菌和真菌的營養(yǎng)來源;相反,細(xì)菌和真菌將有機(jī)分子分解為簡(jiǎn)單的無機(jī)物,這些無機(jī)物又是生產(chǎn)者的營養(yǎng)源。細(xì)菌作為微生物中含量最多、豐富度最高的類群[13],在土壤生態(tài)過程中有著不可或缺的地位。真菌是土壤生態(tài)系統(tǒng)中重要分解者,可分解植物殘?bào)w及纖維素、半纖維素和木質(zhì)素等難分解化合物,并在分解過程中釋放營養(yǎng)用于植被生長,在生態(tài)系統(tǒng)中同樣占有重要地位。土壤中細(xì)菌和真菌多樣性的變化會(huì)影響土壤生態(tài)過程,例如養(yǎng)分循環(huán)[14- 15]和氣體釋放[16]等。研究發(fā)現(xiàn),細(xì)菌和真菌的多樣性受自然因素和人為因素影響,包括植被[17]、pH值[18]、土壤質(zhì)地[19]、土地利用方式[20]土壤養(yǎng)分[21]等。現(xiàn)有關(guān)消落帶土地利用方式對(duì)土壤的研究多集中于土壤養(yǎng)分、土壤理化特性的影響等方面[22- 23],而對(duì)土壤微生物結(jié)構(gòu)和功能多樣性方面的研究相對(duì)較少。消落帶土地利用方式對(duì)土壤微生物多樣性的影響還不明確。本研究旨在探索不同土地利用方式對(duì)土壤微生物特性的影響,為制定合理的土地利用方式,提高三峽消落帶土壤質(zhì)量和實(shí)現(xiàn)消落帶土壤的可持續(xù)利用提供科學(xué)依據(jù)。

    1 研究地點(diǎn)與方法

    1.1 研究區(qū)概況

    本研究區(qū)域位于三峽庫區(qū)重慶忠縣石寶鎮(zhèn)的汝溪河流域消落帶。忠縣石寶鎮(zhèn)位于重慶市中部(108°08′03″—108°08′21″E,30°24′16″—30°24′56″N),距離忠縣縣城32 km,位于忠縣、石柱、萬州三縣(區(qū))交界處。該地區(qū)為亞熱帶東南季風(fēng)氣候,四季分明。雨量充沛,年平均降雨量1172 mm,雨季多在5—8月,年平均氣溫19.2℃。地帶性森林植被為亞熱帶常綠闊葉林,消落帶原生植被破壞殆盡,現(xiàn)有植被以草本為主;土壤主要為中性紫色土,土層厚度變異很大,水土流失、崩塌現(xiàn)象較嚴(yán)重[24]。為進(jìn)行三峽庫區(qū)消落帶的生態(tài)植被恢復(fù)與重建,2012年3月在該地區(qū)消落帶165—175 m海拔區(qū)間按1 m×1 m的株行距栽種了柳樹(Salixmatsudana)、落羽杉(Taxodiumdistichum)和池杉(Taxodiumascendens)的護(hù)岸林。苗木栽種的岸坡樣地土壤均為紫色土,平均坡度為26°。

    1.2 樣地設(shè)置和樣品采集

    在坡度相似、坡向一致、接受光照輻射強(qiáng)度大致相同、土壤預(yù)熱條件基本一致的消落帶區(qū)域內(nèi),2015年6月,通過實(shí)地踏查,選定人工生態(tài)修復(fù)林地、自然恢復(fù)的棄耕地與仍在繼續(xù)耕種的農(nóng)耕地開展實(shí)驗(yàn)研究。其中,棄耕地由三峽工程動(dòng)工前的部分耕地停耕后轉(zhuǎn)變而來。

    根據(jù)試驗(yàn)地實(shí)際情況、植被類型和受干擾情況,選擇每種用地類型各3個(gè),共9個(gè)樣地,所有樣地盡量保持在同一海拔高度。其中耕地作物為芝麻(SesamumindicumL.),行距平均在20 cm左右,受人為松土、施肥、清除雜草等管理活動(dòng)影響嚴(yán)重。棄耕地的優(yōu)勢(shì)物種為狗牙根(Cynodondactylon(L.) Pers)、狼把草(BidenstripartitaL.),總蓋度約為80%。林地以落羽杉、柳樹和池杉為優(yōu)勢(shì)種,總蓋度約為93%。

    在每個(gè)樣地內(nèi)隨機(jī)設(shè)3條樣帶,避免堆肥料地和田埂、溝邊等特殊地形部位。樣帶中按S形5點(diǎn)取樣,在樣點(diǎn)清理地被物后,采取0—20 cm的土壤,5點(diǎn)土壤等質(zhì)量混合,利用四分法保留土壤樣品。其中1 kg土壤室溫風(fēng)干,用于土壤理化性質(zhì)分析。20 g混合土樣裝入無菌的自封袋,低溫運(yùn)輸,-80℃保存,用于土壤微生物多樣性分析。

    1.3 土壤理化性質(zhì)測(cè)定

    土壤含水量采用烘干法測(cè)定,土壤容重采用環(huán)刀法測(cè)定,土壤密度采用比重瓶法測(cè)定,土壤總孔隙度采用計(jì)算法測(cè)定。土壤pH值測(cè)定采用IQ150土壤原位pH計(jì)(IQ Scientific Instruments, Inc., San Diego, CA,USA)原位測(cè)定;土壤有機(jī)質(zhì)(OM)采用重鉻酸鉀外加熱法測(cè)定;土壤全氮(TN)含量采用元素分析儀(Elementar Vario EL, Germany)測(cè)定;土壤堿解氮(AN)含量采用堿解-擴(kuò)散法測(cè)定;全鉀(TK)、速效鉀(AK)含量采用原子吸收光譜儀(AA800, PE Inc., USA)測(cè)定;土壤全磷(TP)、速效磷(AP)含量采用鉬銻抗比色法測(cè)定[25]。

    1.4 土壤微生物多樣性測(cè)定

    土壤微生物多樣性測(cè)定方法為限制性末端長度多態(tài)性[26]。具體測(cè)定步驟如下:

    1.4.1 DNA提取

    用Ultra CleanTMSoil Isolation DNA Kit (MOBIO,USA)提取土壤中的總DNA,用1%的瓊脂糖凝膠電泳檢測(cè)所提DNA片段,并用 UV spectrophotometer (SMA1000, Meriton) 測(cè)定DNA 的濃度。

    1.4.2 PCR擴(kuò)增

    細(xì)菌的擴(kuò)增引物為63F(5′-CAGGCCTAACACATGCAAGTC- 3′)/ 1389R(5′- ACGGGCGGTGTGTACAAG- 3′),正向引物5′端用6-FAM進(jìn)行熒光標(biāo)記。真菌的擴(kuò)增引物為ITS1-F(5′-CTTGGTCATTTAGAGGAAGTAA- 3′)/ITS4(5′-TCCTCCGCTTATTGATAGC- 3′),正向引物5′端用6-FAM進(jìn)行熒光標(biāo)記。本實(shí)驗(yàn)的所用引物由深圳華大基因公司合成并標(biāo)記。

    細(xì)菌50 μL的PCR反應(yīng)體系組成25 μL 2×Taq Master Mix,2μL DNA 模板,63F/ 1389R (10μmol /L) 各0.5μL,ddH2O 補(bǔ)水至50μL。PCR 反應(yīng)條件95℃,5 min。30個(gè)循環(huán)為: 94℃,1 min;56℃,1 min;72℃,2 min;最后72℃保溫10 min。真菌50 μL的PCR反應(yīng)體系組成如下:25μL 2×Taq Master Mix,2μL DNA 模板,ITS1-F 和ITS4 (10 μmol /L) 各1μL,ddH2O 補(bǔ)水至50μL。PCR反應(yīng)條件如下:95℃,5 min。30個(gè)循環(huán)為: 94℃,1min;56℃,1min;72℃,2min;最后72℃保溫10min。每個(gè)樣品重復(fù)2管。

    1.4.3 PCR產(chǎn)物純化、酶切及 T-RFLP

    擴(kuò)增產(chǎn)物經(jīng)1% 瓊脂糖凝膠電泳檢測(cè),按照PCR產(chǎn)物純化試劑盒(TIANGEN,China)說明書進(jìn)行PCR產(chǎn)物純化,-20℃保存?zhèn)溆谩?/p>

    純化產(chǎn)物用限制性內(nèi)切酶HhaI/MspI消化,反應(yīng)體系: 10×Buffer 2 μL,HhaI/MspI (10 U/μL) 1μL,DNA 400 ng,ddH2O補(bǔ)足至30 μL。將酶切樣品置于37 ℃水浴中溫育3 h,酶切完畢后65 ℃水浴10 min 終止反應(yīng)。每個(gè)樣品按上述同樣條件各重復(fù)3次。酶切產(chǎn)物送至上海Invitrogen生命技術(shù)公司(Life Technologies)進(jìn)行基因掃描,得到T-RFLP圖譜。

    1.4.4 T-RFLP 分析

    在限制性片段(Terminal Restriction Fragment, T-RF)的選擇中,刪除引物峰小于50 bp和熒光值小于50RFU的T-RFs,并去除OTU(Operational Taxonomic Unit)豐度<1%的T-RFs,T-RFs 片段大小±1 bp被認(rèn)為是同一個(gè)OTU[27- 28]。

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

    以T-RFLP圖譜中T-RF為一個(gè)OTU,以該T-RF的相對(duì)峰面積作為對(duì)應(yīng)的OTU的豐度并計(jì)算以下多樣性指數(shù):

    Shannon-Wiener指數(shù)(H)

    Simpson指數(shù)(D)

    均勻度指數(shù)(E)

    E=H/lnS

    式中,Pi為某個(gè)峰的峰高占總峰高的比例;S代表物種豐度,即圖譜中T-RF的總數(shù)。

    采用單因素方差分析(One-way ANOVA) 分析不同用地類型對(duì)土壤性質(zhì)和微生物多樣性的影響,并用Duncan 檢驗(yàn)法檢驗(yàn)每個(gè)指標(biāo)在不同用地類型間的差異顯著性(α=0.05)。利用Canoco 4.5軟件進(jìn)行冗余分析(Redundancy Analysis,RDA),通過蒙特卡羅置換檢驗(yàn)(Monte-Carlo permutation test)分析環(huán)境變量對(duì)土壤細(xì)菌與真菌群落變異解釋度的顯著性(α=0.05)。用Origin 8.5和Cano Draw for Windows 4.5軟件制圖。

    2 結(jié)果分析

    2.1 不同土地利用方式下消落帶土壤的理化性狀

    由表1的數(shù)據(jù)可知,除了容重和密度在各用地方式之間沒有顯著性差異外,其余各理化指標(biāo)在不同土地利用方式之間有顯著性差異(P<0.05)。其中,林地的含水量顯著高于耕地和棄耕地;耕地的總孔隙度和pH值均顯著高于林地和棄耕地。林地的有機(jī)質(zhì)、全氮、速效氮、速效鉀和速效磷的含量高于耕地和棄耕地,但全鉀含量低于耕地和棄耕地,且均達(dá)顯著差異。同時(shí),林地和耕地的全磷含量顯著高于棄耕地。

    表1 不同土地利用類型樣地中土壤理化性質(zhì)(平均值±標(biāo)準(zhǔn)誤)

    多重均值比較采用Duncan 檢驗(yàn)法; 同列不同字母表示差異顯著(P<0.05)

    2.2 不同土地利用方式下土壤細(xì)菌多樣性

    根據(jù)T-RFLP圖譜中的OTU的數(shù)量,相對(duì)豐度分別計(jì)算了不同土地利用方式下土壤的細(xì)菌多樣性(圖1):Shannon-Weiner指數(shù)、Simpson指數(shù)和均勻度指數(shù)。圖1表明,在不同的土地利用方式下土壤細(xì)菌多樣性之間存在著顯著性差異。HhaI和MspI酶切結(jié)果相似,在不同用地方式下土壤細(xì)菌的Shannon-Weiner指數(shù)、Simpson指數(shù)和均勻度指數(shù)均表現(xiàn)為耕地和林地之間無顯著性差異,但林地和耕地的細(xì)菌多樣性指數(shù)均顯著高于棄耕地。

    圖1 不同用地方式土壤細(xì)菌群落多樣性指數(shù)分析Fig.1 Analysis of bacterial diversity in different land use type digested by HhaI and MspI

    2.3 環(huán)境變量對(duì)于細(xì)菌多樣性的貢獻(xiàn)

    RDA分析(圖2)從整體上反映了不同土地利用方式下土壤細(xì)菌多樣性與12 種環(huán)境因子量之間的關(guān)系。圖中箭頭連線與排序軸夾角的大小表示因子與排序軸相關(guān)性大小,夾角越小說明關(guān)系越密切,箭頭所處的象限表示環(huán)境因子與排序軸之間的正、負(fù)相關(guān)性。結(jié)果顯示,在HhaI酶切下,第1序軸解釋了所有信息的89.0%,第2序軸解釋了所有信息的10.5%,前兩軸可以解釋99.5%(圖2a)。經(jīng)過蒙特卡羅檢驗(yàn)顯示TP(F=11.0,P=0.012)對(duì)細(xì)菌群多樣性的影響達(dá)到顯著水平(P<0.05)。在MspI酶切下,第1序軸解釋了所有信息的97.5%,第2序軸解釋了所有信息的1.9%,前兩軸可以解釋99.4%(圖2b)。經(jīng)過蒙特卡羅檢驗(yàn),TP(F=24.26,P=0.004)和AP(F=22.78,P=0.008)對(duì)細(xì)菌群落多樣性的影響達(dá)到顯著水平(P<0.05)。

    圖2 環(huán)境因素對(duì)土壤細(xì)菌群落結(jié)構(gòu)特征影響的冗余分析Fig.2 Redundancy analysis (RDA) of soil bacterial community structure and environmental factors

    2.4 不同土地利用方式下土壤真菌多樣性

    根據(jù)T-RFLP圖譜中的OTU的數(shù)量,相對(duì)豐度分別計(jì)算了不同土地利用方式下真菌多樣性(圖3):Shannon-Weiner指數(shù)、Simpson指數(shù)和均勻度指數(shù)。不同土地利用方式下土壤真菌多樣性之間存在著顯著性差異。HhaI酶切結(jié)果顯示,棄耕地的真菌群落Shannon-Weiner指數(shù)和Simpson指數(shù)均顯著高于耕地,林地、耕地和棄耕地的真菌均勻度指數(shù)沒有顯著差異。MspI酶切結(jié)果顯示,林地和棄耕地真菌的Shannon-Weiner指數(shù)之間沒有顯著性差異,且均顯著高于耕地;林地、耕地和棄耕地的真菌Simpson指數(shù)和均勻度指數(shù)沒有顯著差異。

    圖3 不同用地方式土壤真菌群落多樣性指數(shù)分析Fig.3 Analysis of fungal diversity in different land use type digested by HhaI and MspI

    2.5 環(huán)境變量對(duì)于真菌多樣性的貢獻(xiàn)

    圖4從整體上反映了不同土地利用方式下土壤真菌多樣性與12種環(huán)境因子量之間的關(guān)系。RDA分析顯示,在HhaI酶切下,第1序軸解釋了所有信息的85.1%,第2序軸解釋了所有信息的14.8%,前兩軸可以解釋99.9%(圖4a)。經(jīng)過蒙特卡羅檢驗(yàn)顯示pH(F=9.47,P=0.012)、OM(F=6.62,P=0.038)和TP(F=5.21,P=0.032)對(duì)真菌群落多樣性的影響達(dá)到顯著水平(P<0.05)。在MspI酶切下,第1序軸解釋了所有信息的63.7%,第2序軸解釋了所有信息的34.2%,前兩軸可以解釋97.9%(圖4b)。經(jīng)過蒙特卡羅檢驗(yàn),pH(F=7.92,P=0.006)、OM(F=8.09,P=0.002)和TN(F=7.85,P=0.004)對(duì)真菌群落多樣性的影響達(dá)到顯著水平(P<0.05)。

    圖4 環(huán)境因素對(duì)土壤真菌群落結(jié)構(gòu)特征影響的冗余分析Fig.4 Redundancy analysis (RDA) of soil fungal community structure and environmental factors○代表耕地,△代表林地,◇代表?xiàng)壐?/p>

    3 討論

    土壤是由氣相、液相和固相三相合一的生物賴以生存的重要載體,其成分的改變會(huì)顯著改變微生物的生境。在三峽庫區(qū)消落帶,不同的土地利用方式在一定程度上造成了土壤理化性質(zhì)的差異。消落帶土壤在水淹和水體沖刷下,表層土壤被侵蝕和剝離,土壤遭到嚴(yán)重侵蝕。研究發(fā)現(xiàn),在消落帶植被恢復(fù)過程中,人工恢復(fù)林草的水土保持能力高于自然恢復(fù)的棄耕地[29- 30],所以在消落帶林地的土壤含水量顯著高于棄耕地。耕地中人為翻動(dòng)使得土壤的總孔隙度高于林地和棄耕地,并且容重低于林地和棄耕地。

    土壤微生物群落的定性和定量變化是監(jiān)測(cè)土壤質(zhì)量短期和長期變化的敏感指標(biāo)[31]。微生物多樣性指數(shù)是評(píng)價(jià)不同土壤微生物群落多樣性的有效方法,多樣性指數(shù)越高表明微生物群落多樣性越高。探究消落帶土壤中微生物的多樣性有利于進(jìn)一步了解土壤微生物的結(jié)構(gòu)和功能,從而反映出消落帶土壤生態(tài)系統(tǒng)的健康狀況。

    在三峽消落帶不同土地利用方式下土壤細(xì)菌多樣性有顯著差異,棄耕地細(xì)菌多樣性值最低,此結(jié)果與秦紅靈[32]的研究結(jié)果(農(nóng)田土壤的細(xì)菌多樣性顯著高于棄耕地)一致。在本研究中,全磷和速效磷的含量是影響細(xì)菌群落結(jié)構(gòu)的重要因素,與Wakelina[33]和Grayston[34]的研究相類似。土地利用方式?jīng)Q定地表的植被型和土壤的管理方式,進(jìn)而影響到土壤養(yǎng)分[35- 36]。土壤養(yǎng)分的變化勢(shì)必會(huì)影響到微生物群落結(jié)構(gòu),而不同程度的養(yǎng)分添加會(huì)增加微生物群落結(jié)構(gòu)和功能多樣性[37]。細(xì)菌是土壤中多樣性最豐富的微生物類群[38],對(duì)環(huán)境的變異極其敏感。細(xì)菌細(xì)胞壁含有大量磷壁酸,特別在革蘭氏陽性菌中,細(xì)胞壁中磷壁酸的含量約占細(xì)胞干重的50%[39]。在本研究中林地和耕地的全磷和速效磷含量均顯著高于棄耕地,原因可能是林地冠層對(duì)于雨水的攔截使得地表徑流變小,土壤表層有機(jī)質(zhì)以及礦質(zhì)營養(yǎng)得以保留,流失較少,使得林地土壤養(yǎng)分含量均顯著高于棄耕地;而耕地后期人為施肥,使得土壤中營養(yǎng)得以補(bǔ)償[40];但是棄耕地水淹和干旱交替脅迫的同時(shí)沒有外來營養(yǎng)的補(bǔ)充,植被的生長又吸收了土壤中吸附的磷元素,最終導(dǎo)致土壤中全磷和速效磷的含量均較低。土壤磷不僅是土壤細(xì)菌群落的能源之一,同時(shí)也是植物生長的重要資源[41]。雖然三峽消落帶的特殊生境使得土壤對(duì)磷的吸附量增加[42],但土壤磷主要來源于土壤母質(zhì),其來源單一,加之在三峽庫區(qū)消落帶植被生長盛期,不同土地利用方式下土壤中的磷可能在植物-土壤-微生物系統(tǒng)中重新分配,使得磷含量成為影響土壤細(xì)菌群落變異的重要因子之一[43]。

    真菌是一個(gè)生態(tài)幅度較廣的菌群[44- 45],影響其多樣性的因素多為土壤的管理方式[46]和營養(yǎng)水平[47- 48]。本研究中耕地的真菌多樣性最低,與E. Gomez[49]的結(jié)果一致。RDA分析顯示,pH、有機(jī)質(zhì)和全氮的含量對(duì)真菌群落結(jié)構(gòu)的影響達(dá)到顯著水平(P<0.05)。有研究表明,真菌在酸性環(huán)境中提高對(duì)水分的利用率,從而促進(jìn)自身的生長[50- 51]。本研究中,林地和棄耕地土壤pH為中性,耕地土壤pH顯著高于林地和棄耕地,呈堿性,不利于真菌的生長[43]。同時(shí)真菌多以菌絲的形態(tài)生長和繁衍,土壤團(tuán)聚體的穩(wěn)定性顯的尤為重要[52]。與林地和棄耕地相比較,耕地傳統(tǒng)的農(nóng)事耕作破壞大團(tuán)聚體,使得土壤機(jī)械結(jié)構(gòu)變動(dòng)較大,不利于真菌菌絲的延伸[53]。林地和棄耕地較耕地而言人為干擾相對(duì)較少,土壤機(jī)械結(jié)構(gòu)較穩(wěn)定,有利于土壤形成大的團(tuán)聚體和有機(jī)物被固定于土壤中[46],有利于增加土壤微生物生物量并保持微生物多樣性[49]。有機(jī)質(zhì)和全氮是真菌重要的碳和氮源,林地和棄耕地枯枝落葉比較豐富,含有多種枯枝落葉所形成的有機(jī)質(zhì),可以促進(jìn)真菌的生長[54];其林冠層和高的植被覆蓋度可有效減少徑流,使得表層的養(yǎng)分含量得以保留[55];而耕地中人為活動(dòng)和地表徑流加速了土壤中有機(jī)質(zhì)和營養(yǎng)元素的流失[56],有機(jī)質(zhì)和全氮的含量成為消落帶中真菌生長和繁衍的限制性因子。

    微生物多樣性降低或喪失直接影響微生物的生態(tài)服務(wù)功能和整個(gè)生態(tài)系統(tǒng)[57]。消落帶耕地雖然土壤養(yǎng)分得到補(bǔ)充,細(xì)菌的多樣性與林地之間沒有顯著差異,但是其真菌多樣性降低;同時(shí)長時(shí)間的耕作會(huì)威脅到土壤質(zhì)量[58],因而消落帶內(nèi)不宜耕作,應(yīng)當(dāng)限制開墾。三峽消落帶實(shí)行生態(tài)修復(fù)后,棄耕地大量產(chǎn)生;但消落帶水淹-干旱的特殊生境,使棄耕地土壤的養(yǎng)分含量較低,不利于細(xì)菌的生長和繁衍。在三峽消落帶人工生態(tài)修復(fù)林地的土壤養(yǎng)分含量、細(xì)菌多樣性和真菌多樣性均優(yōu)于耕地和棄耕地,有利于維持土壤生態(tài)系統(tǒng)的功能,因此,林地是消落帶最為適宜的用地方式。

    4 結(jié)論

    在三峽庫區(qū)消落帶不同土地利用方式下土壤細(xì)菌和真菌多樣性有顯著性差異。林地和耕地的細(xì)菌多樣性無顯著性差異,且均高于棄耕地;在不同用地方式下,土壤全磷和速效磷的含量顯著影響細(xì)菌群落結(jié)構(gòu)(P<0.05)。真菌群落的Shannon-Weiner指數(shù)和Simpson指數(shù)在耕地中均為最低,但均勻度指數(shù)在3種用地類型之間沒有顯著性差異;在不同用地方式的土壤中,土壤有機(jī)質(zhì)、全氮含量和pH顯著影響真菌群落結(jié)構(gòu)(P<0.05)。與棄耕地和耕地相比,林地可穩(wěn)定消落帶土壤環(huán)境,固持、滯留和保淹有更多的N、P、K等養(yǎng)分,減少其直接進(jìn)入水體的量,因而有助于水庫水質(zhì)的保護(hù);同時(shí)林地可提高微生物多樣性,對(duì)土壤的生態(tài)恢復(fù)發(fā)揮了重要作用。建議在消落帶165—175 m海拔高程內(nèi),應(yīng)限制耕作,推廣人工生態(tài)修復(fù)林地建設(shè)。

    [1] 張建春, 彭補(bǔ)拙. 河岸帶研究及其退化生態(tài)系統(tǒng)的恢復(fù)與重建. 生態(tài)學(xué)報(bào), 2003, 23(1): 56- 63.

    [2] 馬駿, 李昌曉, 魏虹, 馬朋, 楊予靜, 任慶水, 張?chǎng)? 三峽庫區(qū)生態(tài)脆弱性評(píng)價(jià). 生態(tài)學(xué)報(bào), 2015, 35(21): 7117- 7129.

    [3] 劉曉, 蘇維詞, 王錚, 黃一民, 鄧吉祥. 基于RRM模型的三峽庫區(qū)重慶開縣消落區(qū)土地利用生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 2012, 32(1): 248- 256.

    [4] 鮑玉海, 賀秀斌. 三峽水庫消落帶土壤侵蝕問題初步探討. 水土保持研究, 2011, 18(6): 190- 195.

    [5] Yuan X Z, Zhang Y W, Liu H, Xiong S, Li B, Deng W. The littoral zone in the Three Gorges Reservoir, China: challenges and opportunities. Environmental Science and Pollution Research, 2013, 20(10): 7092- 7102.

    [6] 許川, 舒為群, 曹佳, 汪洋. 三峽庫區(qū)消落帶富營養(yǎng)化及其危害預(yù)測(cè)和防治. 長江流域資源與環(huán)境, 2005, 14(4): 440- 444.

    [7] Shen Z Y, Gong Y W, Li Y H, Liu R M. Analysis and modeling of soil conservation measures in the Three Gorges Reservoir Area in China. Catena, 2010, 81(2): 104- 112.

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

    [9] Xu H J, Wang X H, Li H, Yao H Y, Su J Q, Zhu Y G. Biochar impacts soil microbial community composition and nitrogen cycling in an acidic soil planted with rape. Environmental Science & Technology, 2014, 48(16): 9391- 9399.

    [10] Zhang Y, Du B H, Jin Z G, Li Z H, Song H N, Ding Y Q. Analysis of bacterial communities in rhizosphere soil of healthy and diseased cotton (Gossypiumsp.) at different plant growth stages. Plant and Soil, 2011, 339(1/2): 447- 455.

    [11] McCann K S. The diversity-stability debate. Nature, 2000, 405(6783): 228- 233.

    [12] Gans J, Wolinsky M, Dunbar J. Computational improvements reveal great bacterial diversity and high metal toxicity in soil. Science, 2005, 309(5739): 1387- 1390.

    [13] Kennedy A C. Bacterial diversity in agroecosystems. Agriculture, Ecosystems & Environment, 1999, 74(1/3): 65- 76.

    [14] Unger I M, Goyne K W, Kremer R J, Kennedy A C. Microbial community diversity in agroforestry and grass vegetative filter strips. Agroforestry Systems, 2013, 87(2): 395- 402.

    [15] Lienhard P, Terrat S, Mathieu O, Levêque J, Chemidlin Prévost-Bouré N, Nowak V, Régnier T, Faivre C, Sayphoummie S, Panyasiri K, Tivet F, Ranjard L, Maron P A. Soil microbial diversity and C turnover modified by tillage and cropping in Laos tropical grassland. Environmental Chemistry Letters, 2013, 11(4): 391- 398.

    [16] Bardgett R D, Freeman C, Ostle N J. Microbial contributions to climate change through carbon cycle feedbacks. The ISME Journal, 2008, 2(8): 805- 814.

    [17] Zu Y G, Gao C Y, Wang W J, Yang F J, Liu Y, Wang M, Zhao Y G. Characteristics of the microbial community in rhizosphere ofCamptothecaacuminatacultured with exotic invasive plantEupatoriumadenophorum. Science in China Series C: Life Sciences, 2007, 50(1): 22- 30.

    [18] Tripathi B M, Kim M, Singh D, Lee-Cruz L, Lai-Hoe A, Ainuddin A N, Go R, Rahim R A, Husni M H A, Chun J, Adams J M. Tropical soil bacterial communities in Malaysia: pH dominates in the equatorial tropics too. Microbial Ecology, 2012, 64(2): 474- 484.

    [19] Schutter M E, Sandeno J M, Dick R P. Seasonal, soil type, and alternative management influences on microbial communities of vegetable cropping systems. Biology and Fertility of Soils, 2001, 34(6): 397- 410.

    [20] Suleiman A K A, Manoeli L, Boldo J T, Pereira M G, Roesch L F W. Shifts in soil bacterial community after eight years of land-use change. Systematic and Applied Microbiology, 2013, 36(2): 137- 144.

    [21] Ma L N, Huang W W, Guo C Y, Wang R Z, Xiao C W. Soil microbial properties and plant growth responses to carbon and water addition in a temperate steppe: the importance of nutrient availability. Plos One, 2012, 7(4): e35165.

    [22] 張雅蓉, 高明, 黃容. 三峽庫區(qū)消落帶土壤無機(jī)磷的生物有效性. 水土保持學(xué)報(bào), 2014, 28(3): 222- 226.

    [23] 王曉榮, 程瑞梅, 肖文發(fā), 郭泉水, 封曉輝. 三峽庫區(qū)消落帶初期土壤養(yǎng)分特征. 生態(tài)學(xué)雜志, 2010, 29(2): 281- 289.

    [24] 楊予靜, 李昌曉, 麗娜·熱瑪贊. 基于PSR框架模型的三峽庫區(qū)忠縣汝溪河流域生態(tài)系統(tǒng)健康評(píng)價(jià). 長江流域資源與環(huán)境, 2013, 22(S1): 66- 74.

    [25] 鮑士旦. 土壤農(nóng)化分析(第三版). 北京: 中國農(nóng)業(yè)出版社, 1999: 30- 109.

    [26] Kitts C L. Terminal restriction fragment patterns: a tool for comparing microbial communities and assessing community dynamics. Current Issues in Intestinal Microbiology, 2001, 2(1): 17- 25.

    [27] Schütte U M E, Abdo Z, Bent S J, Shyu C, Williams C J, Pierson J D, Forney L J. Advances in the use of terminal restriction fragment length polymorphism (T-RFLP) analysis of 16S rRNA genes to characterize microbial communities. Applied Microbiology and Biotechnology, 2008, 80(3): 365- 380.

    [28] Székely A J, Sipos R, Berta B, Vajna B, Hajdú C, Márialigeti K. DGGE and T-RFLP analysis of bacterial succession during mushroom compost production and sequence-aided T-RFLP profile of mature compost. Microbial Ecology, 2009, 57(3): 522- 533.

    [29] 黃茹. 三峽庫區(qū)林草治理措施對(duì)土壤理化特征及坡面水沙的影響[D]. 重慶: 西南大學(xué), 2013: 29- 35.

    [30] 劉斌, 羅全華, 常文哲, 許小梅. 不同林草植被覆蓋度的水土保持效益及適宜植被覆蓋度. 中國水土保持科學(xué), 2008, 6(6): 68- 73.

    [31] 肖國生, 胡廷章, 唐華麗, 楊俊年, 李彥杰, 鄭克均, 王啟龍. 三峽水庫消落帶淹沒前后土壤微生物生態(tài)分布及優(yōu)勢(shì)菌群的鑒定. 江蘇農(nóng)業(yè)科學(xué), 2011, 39(4): 493- 496.

    [32] 秦紅靈, 袁紅朝, 張慧, 朱亦君, 吳敏娜, 魏文學(xué). 紅壤坡地利用方式對(duì)土壤細(xì)菌群落結(jié)構(gòu)的影響. 土壤學(xué)報(bào), 2011, 48(3): 594- 602.

    [33] Wakelin S A, Macdonald L M, Rogers S L, Gregg A L, Bolger T P, Baldock J A. Habitat selective factors influencing the structural composition and functional capacity of microbial communities in agricultural soils. Soil Biology and Biochemistry, 2008, 40(3): 803- 813.

    [34] Grayston S J, Campbell C D, Bardgett R D, Mawdsley J L, Clegg C D, Ritz K, Griffiths B S, Rodwell J S, Edwards S J, Davies W J, Elston D J, Millard P. Assessing shifts in microbial community structure across a range of grasslands of differing management intensity using CLPP, PLFA and community DNA techniques. Applied Soil Ecology, 2004, 25(1): 63- 84.

    [35] Oberson A, Friesen D K, Rao I M, Bühler S, Frossard E. Phosphorus Transformations in an Oxisol under contrasting land-use systems: the role of the soil microbial biomass. Plant and Soil, 2001, 237(2): 197- 210.

    [36] Batlle-Aguilar J, Brovelli A, Porporato A, Barry D A. Modelling soil carbon and nitrogen cycles during land use change. A review. Agronomy for Sustainable Development, 2011, 31(2): 251- 274.

    [37] Sarathchandra S U, Ghani A, Yeates G W, Burch G, Cox N R. Effect of nitrogen and phosphate fertilisers on microbial and nematode diversity in pasture soils. Soil Biology and Biochemistry, 2001, 33(7/8): 953- 964.

    [38] 黃昌勇. 土壤學(xué). 北京: 中國農(nóng)業(yè)出版社, 2000: 3- 78.

    [39] Endl J, Seidl H P, Fiedler F, Schleifer K H. Chemical composition and structure of cell wall teichoic acids of staphylococci. Archives of Microbiology, 1983, 135(3): 215- 223.

    [40] Jesus E da C, Marsh T L, Tiedje J M, Moreira F M da S. Changes in land use alter the structure of bacterial communities in Western Amazon soils. The ISME Journal, 2009, 3(9): 1004- 1011.

    [41] Tarafdar J C, Claassen N. Organic phosphorus compounds as a phosphorus source for higher plants through the activity of phosphatases produced by plant roots and microorganisms. Biology and Fertility of Soils, 1988, 5(4): 308- 312.

    [42] 朱強(qiáng), 安然, 胡紅青, 萬成炎, 胡蓮, 王素梅. 三峽庫區(qū)消落帶土壤對(duì)磷的吸附和淹水下磷的形態(tài)變化. 土壤學(xué)報(bào), 2012, 49(6): 1128- 1135.

    [43] Rezapour S. Response of some soil attributes to different land use types in calcareous soils with Mediterranean type climate in north-west of Iran. Environmental Earth Sciences, 2014, 71(5): 2199- 2210.

    [44] Tedersoo L, Bahram M, P?lme S, K?ljalg U, Yorou N S, Wijesundera R, Ruiz L V, Vasco-Palacios A M, Thu P Q, Suija A, Smith M E, Sharp C, Saluveer E, Saitta A, Rosas M, Riit T, Ratkowsky D, Pritsch K, P?ldmaa K, Piepenbring M, Phosri C, Peterson M, Parts K, P?rtel K, Otsing E, Nouhra E, Njouonkou A L, Nilsson R H, Morgado L N, Mayor J, May T W, Majuakim L, Lodge D J, Lee S S, Larsson K H, Kohout P, Hosaka K, Hiiesalu I, Henkel T W, Harend H, Guo L D, Greslebin A, Grelet G, Geml J, Gates G, Dunstan W, Dunk C, Drenkhan R, Dearnaley J, De Kesel A, Dang T, Chen X, Buegger F, Brearley F Q, Bonito G, Anslan S, Abell S, Abarenkov K. Global diversity and geography of soil fungi. Science, 2014, 346(6213): 1256688- 1256688.

    [45] Rousk J, B??th E, Brookes P C, Lauber C L, Lozupone C, Caporaso G, Knight R, Fierer N. Soil bacterial and fungal communities across a pH gradient in an arable soil. The ISME Journal, 2010, 4(10): 1340- 1351.

    [46] 羅紅燕, 蔣先軍, 謝德體, 李楠, 曹良元. 真菌和細(xì)菌生物量在土壤團(tuán)聚體中的分布和耕作響應(yīng). 生態(tài)學(xué)報(bào), 2009, 29(8): 4588- 4594.

    [47] Wang S K, Zuo X A, Zhao X Y, Li Y Q, Zhou X, Lv P, Luo Y Q, Yun J Y. Responses of soil fungal community to the sandy grassland restoration in Horqin Sandy Land, northern China. Environmental Monitoring and Assessment, 2016, 188(1): 21- 21.

    [48] 高玉峰, 賀字典. 影響土壤真菌多樣性的土壤因素. 中國農(nóng)學(xué)通報(bào), 2010, 26(10): 177- 181.

    [49] Gomez E, Pioli R, Conti M. Fungal abundance and distribution as influenced by clearing and land use in a vertic soil of Argentina. Biology and Fertility of Soils, 2006, 43(3): 373- 377.

    [50] Nevarez L, Vasseur V, Le Madec A, Le Bras M A, Coroller L, Leguérinel I, Barbier G. Physiological traits ofPenicilliumglabrumstrain LCP 08.5568, a filamentous fungus isolated from bottled aromatised mineral water. International Journal of Food Microbiology, 2009, 130(3): 166- 171.

    [51] Gock M A, Hocking A D, Pitt J I, Poulos P G. Influence of temperature, water activity and pH on growth of some xerophilic fungi. International Journal of Food Microbiology, 2003, 81(1): 11- 19.

    [52] Spurgeon D J, Keith A M, Schmidt O, Lammertsma D R, Faber J H. Land-use and land-management change: relationships with earthworm and fungi communities and soil structural properties. BMC Ecology, 2013, 13(1): 46- 46.

    [53] 李景, 吳會(huì)軍, 武雪萍, 蔡典雄, 姚宇卿, 呂軍杰, 田云龍. 長期不同耕作措施對(duì)土壤團(tuán)聚體特征及微生物多樣性的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(8): 2341- 2348.

    [54] Hedlund K. Soil microbial community structure in relation to vegetation management on former agricultural land. Soil Biology and Biochemistry, 2002, 34(9): 1299- 1307.

    [55] 董莉麗, 鄭粉莉. 土地利用類型對(duì)土壤微生物量和有機(jī)質(zhì)的影響. 水土保持通報(bào), 2009, 29(6): 10- 15.

    [56] 楊景成, 韓興國, 黃建輝, 潘慶民. 土壤有機(jī)質(zhì)對(duì)農(nóng)田管理措施的動(dòng)態(tài)響應(yīng). 生態(tài)學(xué)報(bào), 2003, 23(4): 787- 796.

    [57] Cardinale B J, Srivastava D S, Duffy J E, Wright J P, Downing A L, Sankaran M, Jouseau C. Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature, 2006, 443(7114): 989- 992.

    [58] Kinyangi J M. Soil degradation, thresholds and dynamics of long-term cultivation: from landscape biogeochemistry to nanoscale biogeocomplexity[D]. Cornell: Cornell University, 2008: 35- 37.

    Effects of different land use patterns on soil bacterial and fungal biodiversity in the hydro-fluctuation zone of the Three Gorges Reservoir region

    QIN Hong,LI Changxiao*, REN Qingshui

    KeyLaboratoryofEco-EnvironmentsinThreeGorgesReservoirRegion(MinistryofEducation),ChongqingKeyLaboratoryofPlantEcologyandResourcesinThreeGorgesReservoirRegion,schoolofLifeSciences,SouthwestUniversity,Chongqing400715,China

    The goal of this study was to provide references for the environmental management of the Three Gorges Reservoir. To study the effect of different land use patterns on soil microorganism biodiversity in the hydro-fluctuation zone of the Three Gorges Reservoir. Surface soil samples were collected from the woodland, abandoned farmland, and farmland in June 2015 to determine the bacterial and fungal community biodiversity using the terminal-restriction fragment length polymorphism (T-RFLP) and soil chemical properties in the hydro-fluctuation zone of the Ruxi River basin, in Zhongxian County, Chongqing. Soil chemical properties were also investigated. Results showed that (1) the land use patterns significantly affected soil physicochemical properties, excluding bulk weight and density. Other soil physicochemical properties, such as soil moisture, content of soil organic matter, total nitrogen, total phosphorus, available nitrogen, available potassium, and available phosphorus, in the woodland were significantly higher than those of the other two types of land use (P<0.05); (2) there were significant differences among land use patterns in bacterial and fungal biodiversity; (3) the bacterial diversity index of the woodland was not significantly different from farmland, whereas both of them were significantly higher than abandoned farmland. Redundancy analysis and the Monte-Carlo test revealed that total phosphorus and available phosphorus showed a critical influence on bacterial diversity (P<0.05); and (4) the Shannon-Wiener index and the Simpson index of fungal community for the woodland and abandoned farmland were significantly higher than those of farmland, whereas there was not significant difference in richness index among the three types. Redundancy analysis and the Monte-Carlo test revealed that pH, organic matter, and total nitrogen played a critical roles in influencing fungal community diversity (P<0.05). Conclusion: compared to farmland and abandoned farmland, woodland can retain more soil nutrients, and improve soil bacterial, and fungal biodiversity, suggesting that the artificially regenerated woodland is the most appropriate land use type in the altitude of 165 m to 175 m in the hydro-fluctuation zone of the Three Gorges Reservoir. Those studies also propose that farming should be limited and the construction of artificially regenerated woodland in the hydro-fluctuation zone of the Three Gorges Reservoir region should be promoted.

    The Three Gorges Reservoir Area; hydro-fluctuation zone; different land use patterns; bacterial and fungal diversity

    重慶市林業(yè)重點(diǎn)科技攻關(guān)項(xiàng)目(渝林科研2016-8,2015-6);國家國際科技合作專項(xiàng)(2015DFA90900);中央財(cái)政林業(yè)科技推廣示范項(xiàng)目(渝林科推[2014- 10]);重慶市研究生科研創(chuàng)新項(xiàng)目(CYS2015068)

    2016- 03- 03; 網(wǎng)絡(luò)出版日期:2017- 02- 17

    10.5846/stxb201603030374

    *通訊作者Corresponding author.E-mail: lichangx@swu.edu.cn

    秦紅,李昌曉,任慶水.不同土地利用方式對(duì)三峽庫區(qū)消落帶土壤細(xì)菌和真菌多樣性的影響.生態(tài)學(xué)報(bào),2017,37(10):3494- 3504.

    Qin H,Li C X, Ren Q S.Effects of different land use patterns on soil bacterial and fungal biodiversity in the hydro-fluctuation zone of the Three Gorges Reservoir region.Acta Ecologica Sinica,2017,37(10):3494- 3504.

    猜你喜歡
    三峽庫區(qū)林地土地利用
    三峽庫區(qū)萬家壩滑坡變形區(qū)穩(wěn)定性復(fù)核研究
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    丹東市林地分類研究
    淺談林地保護(hù)及恢復(fù)措施
    三峽庫區(qū)產(chǎn)業(yè)培育及結(jié)構(gòu)調(diào)整的思考
    濱??h土地利用挖潛方向在哪里
    林地流轉(zhuǎn)模式的選擇機(jī)理及其政策啟示
    小型無人飛行器用于林地監(jiān)視的嘗試
    三峽庫區(qū)生態(tài)環(huán)保成效顯著
    論低碳經(jīng)濟(jì)與轉(zhuǎn)變土地利用方式
    www.精华液| 色播在线永久视频| 人妻少妇偷人精品九色| 午夜影院在线不卡| 国产精品av久久久久免费| 国产 一区精品| 亚洲欧洲精品一区二区精品久久久 | 如日韩欧美国产精品一区二区三区| 日韩中文字幕视频在线看片| 王馨瑶露胸无遮挡在线观看| www.自偷自拍.com| 人妻一区二区av| 成人毛片a级毛片在线播放| 日韩精品有码人妻一区| 亚洲,一卡二卡三卡| 国产日韩欧美视频二区| 国产色婷婷99| 国产色婷婷99| 亚洲欧美成人精品一区二区| 狂野欧美激情性bbbbbb| 国产深夜福利视频在线观看| 男人添女人高潮全过程视频| 黄片小视频在线播放| 国产精品不卡视频一区二区| 久久毛片免费看一区二区三区| 国产午夜精品一二区理论片| 午夜福利,免费看| 亚洲精品美女久久av网站| 成人亚洲精品一区在线观看| av.在线天堂| 这个男人来自地球电影免费观看 | 久久99一区二区三区| 国产黄色免费在线视频| 日产精品乱码卡一卡2卡三| 美女主播在线视频| 一边亲一边摸免费视频| 亚洲 欧美一区二区三区| 九色亚洲精品在线播放| av线在线观看网站| 国产成人午夜福利电影在线观看| 国产黄频视频在线观看| av天堂久久9| 最近的中文字幕免费完整| 精品一品国产午夜福利视频| 免费观看av网站的网址| 男女啪啪激烈高潮av片| 午夜免费男女啪啪视频观看| 少妇熟女欧美另类| 美女视频免费永久观看网站| 精品人妻一区二区三区麻豆| 久久狼人影院| 午夜福利,免费看| 欧美激情高清一区二区三区 | 亚洲中文av在线| 国产成人精品在线电影| 纵有疾风起免费观看全集完整版| 如日韩欧美国产精品一区二区三区| 女人精品久久久久毛片| 亚洲精品国产色婷婷电影| 卡戴珊不雅视频在线播放| 免费黄色在线免费观看| 日本91视频免费播放| 蜜桃在线观看..| 1024视频免费在线观看| 成人18禁高潮啪啪吃奶动态图| 一级黄片播放器| 亚洲欧美成人精品一区二区| 久久精品久久久久久噜噜老黄| 最近中文字幕高清免费大全6| 黄色怎么调成土黄色| 最近中文字幕2019免费版| 久久精品人人爽人人爽视色| 寂寞人妻少妇视频99o| 日韩不卡一区二区三区视频在线| 国产成人91sexporn| 日韩大片免费观看网站| 中国国产av一级| av卡一久久| 亚洲国产精品成人久久小说| 精品国产超薄肉色丝袜足j| 国产乱人偷精品视频| 天天影视国产精品| 99国产精品免费福利视频| 精品国产一区二区三区久久久樱花| 亚洲激情五月婷婷啪啪| 人人妻人人澡人人看| 国产精品 欧美亚洲| 亚洲成人手机| 啦啦啦啦在线视频资源| 精品人妻在线不人妻| 亚洲人成77777在线视频| 黄片播放在线免费| 成人二区视频| 亚洲人成77777在线视频| 一级毛片 在线播放| 一本大道久久a久久精品| 欧美在线黄色| 国产精品无大码| 国产精品女同一区二区软件| 人妻人人澡人人爽人人| 亚洲精品av麻豆狂野| 国产精品国产三级专区第一集| 大码成人一级视频| 日韩免费高清中文字幕av| 麻豆精品久久久久久蜜桃| 男女无遮挡免费网站观看| 精品国产乱码久久久久久男人| 精品国产乱码久久久久久小说| 黑人欧美特级aaaaaa片| 热99国产精品久久久久久7| 国产精品蜜桃在线观看| 亚洲国产看品久久| 黄色配什么色好看| 日韩av不卡免费在线播放| av片东京热男人的天堂| 色视频在线一区二区三区| 丝袜脚勾引网站| 男女午夜视频在线观看| 老汉色av国产亚洲站长工具| 成人毛片60女人毛片免费| 免费播放大片免费观看视频在线观看| 久久ye,这里只有精品| 免费日韩欧美在线观看| 亚洲美女搞黄在线观看| 国产精品二区激情视频| 日韩av在线免费看完整版不卡| 啦啦啦在线观看免费高清www| 大片免费播放器 马上看| 超碰97精品在线观看| 性高湖久久久久久久久免费观看| av福利片在线| a 毛片基地| 欧美激情 高清一区二区三区| 一区在线观看完整版| av卡一久久| 天天影视国产精品| 五月伊人婷婷丁香| 国产老妇伦熟女老妇高清| 免费看不卡的av| 国产一区二区三区av在线| av卡一久久| 巨乳人妻的诱惑在线观看| 91午夜精品亚洲一区二区三区| 黄片播放在线免费| 亚洲欧美色中文字幕在线| 亚洲一区二区三区欧美精品| 午夜日韩欧美国产| 久久久a久久爽久久v久久| 天天影视国产精品| 性色avwww在线观看| 国产精品久久久久成人av| 在线观看免费高清a一片| 中文字幕制服av| 女性被躁到高潮视频| 涩涩av久久男人的天堂| 国产探花极品一区二区| 在线精品无人区一区二区三| 免费女性裸体啪啪无遮挡网站| 亚洲精品,欧美精品| 国产精品三级大全| www日本在线高清视频| 精品亚洲成国产av| 色94色欧美一区二区| 国产欧美日韩综合在线一区二区| 久久99热这里只频精品6学生| 久久精品久久久久久久性| 极品少妇高潮喷水抽搐| 男女无遮挡免费网站观看| 亚洲第一区二区三区不卡| 美女主播在线视频| 91精品国产国语对白视频| 999精品在线视频| 日本黄色日本黄色录像| 欧美人与性动交α欧美精品济南到 | 看免费成人av毛片| 国产精品偷伦视频观看了| 又大又黄又爽视频免费| 久久 成人 亚洲| 精品99又大又爽又粗少妇毛片| 亚洲中文av在线| 国产精品国产三级专区第一集| 国产亚洲午夜精品一区二区久久| 热99国产精品久久久久久7| 精品人妻一区二区三区麻豆| 99国产综合亚洲精品| 最新中文字幕久久久久| 久久精品久久精品一区二区三区| 激情视频va一区二区三区| 蜜桃在线观看..| 免费高清在线观看视频在线观看| 国产精品一区二区在线不卡| 青春草国产在线视频| 黄片小视频在线播放| 国产极品天堂在线| 日韩av免费高清视频| 97精品久久久久久久久久精品| 天天影视国产精品| 男女下面插进去视频免费观看| 色婷婷av一区二区三区视频| 国产1区2区3区精品| 久久精品久久久久久久性| 日本黄色日本黄色录像| 国产极品粉嫩免费观看在线| 中文字幕精品免费在线观看视频| 一区福利在线观看| 亚洲精品日韩在线中文字幕| 中国国产av一级| 飞空精品影院首页| 另类亚洲欧美激情| 午夜福利网站1000一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产一区二区精华液| 国产精品久久久久成人av| 精品亚洲成国产av| 亚洲在久久综合| 午夜av观看不卡| 欧美日韩视频高清一区二区三区二| 9191精品国产免费久久| 久久久久久久亚洲中文字幕| 亚洲欧美一区二区三区黑人 | 不卡av一区二区三区| 欧美成人午夜免费资源| 大香蕉久久网| 日本av免费视频播放| 看免费av毛片| 纯流量卡能插随身wifi吗| 黄色一级大片看看| 国产国语露脸激情在线看| 午夜福利视频精品| www.自偷自拍.com| 性高湖久久久久久久久免费观看| 日韩不卡一区二区三区视频在线| 美女福利国产在线| 日日啪夜夜爽| 狠狠婷婷综合久久久久久88av| 各种免费的搞黄视频| 美女国产高潮福利片在线看| 成人午夜精彩视频在线观看| 日韩av在线免费看完整版不卡| 亚洲情色 制服丝袜| 久热久热在线精品观看| 国产成人精品在线电影| 国产成人精品福利久久| 视频在线观看一区二区三区| 男女下面插进去视频免费观看| 综合色丁香网| 久久精品国产亚洲av天美| 香蕉丝袜av| 亚洲综合色网址| 日韩中文字幕欧美一区二区 | av在线播放精品| av线在线观看网站| 叶爱在线成人免费视频播放| 精品久久蜜臀av无| 曰老女人黄片| 人妻一区二区av| 日韩不卡一区二区三区视频在线| 久久国内精品自在自线图片| 日韩电影二区| 大香蕉久久成人网| 免费观看无遮挡的男女| 十八禁高潮呻吟视频| 高清视频免费观看一区二区| 黑丝袜美女国产一区| 成人毛片60女人毛片免费| 国产成人精品婷婷| 欧美日韩视频高清一区二区三区二| 欧美日韩精品网址| 欧美黄色片欧美黄色片| 国产欧美日韩综合在线一区二区| 精品久久久精品久久久| 国产xxxxx性猛交| 在线观看免费日韩欧美大片| 久久97久久精品| 97在线人人人人妻| 寂寞人妻少妇视频99o| 男女高潮啪啪啪动态图| 成人二区视频| 免费高清在线观看日韩| av免费在线看不卡| 美女国产高潮福利片在线看| 久久av网站| 亚洲美女视频黄频| 美女国产高潮福利片在线看| 精品少妇黑人巨大在线播放| 亚洲美女视频黄频| 亚洲av免费高清在线观看| 国产一区二区激情短视频 | 女性生殖器流出的白浆| av在线播放精品| 日韩在线高清观看一区二区三区| 丰满乱子伦码专区| 秋霞伦理黄片| 色哟哟·www| 丝袜喷水一区| 女人精品久久久久毛片| 男女午夜视频在线观看| 亚洲av男天堂| 国产欧美日韩一区二区三区在线| 欧美日韩成人在线一区二区| 国产乱来视频区| 纵有疾风起免费观看全集完整版| 老司机影院成人| 精品国产一区二区三区久久久樱花| 亚洲精品av麻豆狂野| 人妻一区二区av| av一本久久久久| 青春草视频在线免费观看| 国产97色在线日韩免费| 国产亚洲精品第一综合不卡| 桃花免费在线播放| 国产爽快片一区二区三区| 91国产中文字幕| 一二三四在线观看免费中文在| 午夜福利在线观看免费完整高清在| av卡一久久| 久久鲁丝午夜福利片| 久久久久人妻精品一区果冻| a级片在线免费高清观看视频| 国产精品免费大片| 午夜免费观看性视频| 久久久久视频综合| 女人高潮潮喷娇喘18禁视频| 国产成人91sexporn| 久久人人97超碰香蕉20202| 亚洲人成电影观看| 日韩中文字幕欧美一区二区 | av.在线天堂| 亚洲欧美一区二区三区久久| 熟女电影av网| 超碰成人久久| 免费观看av网站的网址| 日本欧美国产在线视频| 永久免费av网站大全| 中文字幕精品免费在线观看视频| 午夜福利在线观看免费完整高清在| 午夜老司机福利剧场| videosex国产| 大片免费播放器 马上看| 国产一区二区激情短视频 | 高清黄色对白视频在线免费看| 最近手机中文字幕大全| 黄频高清免费视频| 国产老妇伦熟女老妇高清| 日本-黄色视频高清免费观看| 国语对白做爰xxxⅹ性视频网站| 欧美最新免费一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产男女超爽视频在线观看| 亚洲综合精品二区| 天堂中文最新版在线下载| 一本久久精品| 在线看a的网站| 亚洲四区av| 制服丝袜香蕉在线| 久久久久精品性色| 1024视频免费在线观看| 国产精品一二三区在线看| freevideosex欧美| 丰满饥渴人妻一区二区三| 久久99精品国语久久久| 999精品在线视频| 香蕉国产在线看| 亚洲成人一二三区av| 久久精品亚洲av国产电影网| 婷婷成人精品国产| 日韩伦理黄色片| 18禁动态无遮挡网站| 欧美日韩一级在线毛片| www.自偷自拍.com| 中国三级夫妇交换| 精品国产乱码久久久久久男人| 咕卡用的链子| 欧美xxⅹ黑人| 欧美日韩av久久| 日韩av不卡免费在线播放| 一级a爱视频在线免费观看| 久久久久久久久久久久大奶| 久久久久国产精品人妻一区二区| 2022亚洲国产成人精品| 另类精品久久| av免费观看日本| 亚洲国产成人一精品久久久| 性少妇av在线| 97在线人人人人妻| 哪个播放器可以免费观看大片| 国产成人a∨麻豆精品| 丰满迷人的少妇在线观看| 伊人久久国产一区二区| 男女啪啪激烈高潮av片| 99久国产av精品国产电影| 国产白丝娇喘喷水9色精品| 成人漫画全彩无遮挡| 十八禁高潮呻吟视频| 最新的欧美精品一区二区| 毛片一级片免费看久久久久| 久久婷婷青草| 亚洲国产av影院在线观看| 麻豆乱淫一区二区| 丝袜美腿诱惑在线| 日韩成人av中文字幕在线观看| 在线观看一区二区三区激情| 婷婷色综合大香蕉| 亚洲一区中文字幕在线| 国产白丝娇喘喷水9色精品| 久久久久人妻精品一区果冻| 国产一区有黄有色的免费视频| 少妇被粗大猛烈的视频| 侵犯人妻中文字幕一二三四区| 精品国产乱码久久久久久男人| av国产精品久久久久影院| 国产成人午夜福利电影在线观看| 日韩三级伦理在线观看| 国产深夜福利视频在线观看| 欧美日韩视频精品一区| 777久久人妻少妇嫩草av网站| 99精国产麻豆久久婷婷| 中文字幕av电影在线播放| 国产精品不卡视频一区二区| 熟女少妇亚洲综合色aaa.| 老女人水多毛片| 极品人妻少妇av视频| www.自偷自拍.com| 久久久久精品久久久久真实原创| 另类精品久久| 只有这里有精品99| 韩国av在线不卡| 久久久久久久久久久久大奶| 久久久久精品性色| 国产av精品麻豆| 人妻人人澡人人爽人人| 99精国产麻豆久久婷婷| 好男人视频免费观看在线| videosex国产| 日日摸夜夜添夜夜爱| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美,日韩| 精品一区二区三区四区五区乱码 | 丝袜脚勾引网站| 亚洲欧美一区二区三区黑人 | 黄频高清免费视频| 9191精品国产免费久久| 性高湖久久久久久久久免费观看| 五月天丁香电影| 欧美老熟妇乱子伦牲交| 国产免费视频播放在线视频| 一级毛片黄色毛片免费观看视频| 9热在线视频观看99| 超色免费av| 美女视频免费永久观看网站| 日日撸夜夜添| 久热这里只有精品99| 青青草视频在线视频观看| 国产白丝娇喘喷水9色精品| 99久久精品国产国产毛片| 男的添女的下面高潮视频| 韩国高清视频一区二区三区| 搡女人真爽免费视频火全软件| 日本91视频免费播放| 日韩av免费高清视频| 欧美老熟妇乱子伦牲交| 美女福利国产在线| 亚洲av日韩在线播放| 少妇熟女欧美另类| 视频区图区小说| 国产精品麻豆人妻色哟哟久久| 欧美少妇被猛烈插入视频| 国产精品二区激情视频| 人妻人人澡人人爽人人| 午夜日本视频在线| 成年av动漫网址| xxx大片免费视频| 免费看不卡的av| 一二三四中文在线观看免费高清| 欧美xxⅹ黑人| 久久精品久久久久久久性| 老司机午夜十八禁免费视频| av在线播放免费不卡| 久久人人爽av亚洲精品天堂| 成在线人永久免费视频| 欧美日韩av久久| 真人一进一出gif抽搐免费| 女人高潮潮喷娇喘18禁视频| 丰满人妻熟妇乱又伦精品不卡| 日韩三级视频一区二区三区| 久久精品亚洲熟妇少妇任你| 国产色视频综合| 欧美日韩瑟瑟在线播放| 国产精华一区二区三区| 国产野战对白在线观看| 中文字幕最新亚洲高清| 国产精品爽爽va在线观看网站 | 自拍欧美九色日韩亚洲蝌蚪91| 他把我摸到了高潮在线观看| 一区二区三区国产精品乱码| 99国产精品一区二区蜜桃av| 色哟哟哟哟哟哟| 亚洲一区高清亚洲精品| 91麻豆精品激情在线观看国产 | 69av精品久久久久久| 亚洲精品在线美女| 成人国语在线视频| 久久中文看片网| 国产成人精品久久二区二区免费| 久久亚洲真实| 999久久久国产精品视频| 老司机福利观看| 91精品三级在线观看| 一本大道久久a久久精品| 嫁个100分男人电影在线观看| 波多野结衣高清无吗| 精品久久久久久成人av| 国产三级黄色录像| 手机成人av网站| av福利片在线| 午夜a级毛片| 可以在线观看毛片的网站| 亚洲av第一区精品v没综合| 亚洲精品在线美女| 可以免费在线观看a视频的电影网站| 亚洲成人国产一区在线观看| 国产熟女xx| 久久久久九九精品影院| 亚洲成人免费av在线播放| 最新美女视频免费是黄的| 亚洲性夜色夜夜综合| 黄色成人免费大全| 亚洲人成伊人成综合网2020| 亚洲视频免费观看视频| 91九色精品人成在线观看| 国产成人影院久久av| 亚洲精品一区av在线观看| av天堂久久9| av在线天堂中文字幕 | 亚洲黑人精品在线| 18禁观看日本| 在线播放国产精品三级| 黄色视频,在线免费观看| 一级黄色大片毛片| 97超级碰碰碰精品色视频在线观看| 久久精品国产清高在天天线| 91九色精品人成在线观看| 午夜免费激情av| 久久久国产欧美日韩av| 两个人看的免费小视频| 男女之事视频高清在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品国产一区二区精华液| 美女高潮到喷水免费观看| 手机成人av网站| 国产成人av教育| 免费av毛片视频| 色播在线永久视频| 超色免费av| 中文字幕另类日韩欧美亚洲嫩草| 亚洲人成网站在线播放欧美日韩| 亚洲一区高清亚洲精品| 十八禁网站免费在线| 热99国产精品久久久久久7| 久久香蕉激情| 精品一品国产午夜福利视频| 午夜日韩欧美国产| www.www免费av| 精品久久久久久,| 午夜免费鲁丝| 欧美在线一区亚洲| 神马国产精品三级电影在线观看 | 亚洲狠狠婷婷综合久久图片| 亚洲中文字幕日韩| 中出人妻视频一区二区| 女生性感内裤真人,穿戴方法视频| 国产精品98久久久久久宅男小说| 精品无人区乱码1区二区| 亚洲第一青青草原| 热99国产精品久久久久久7| 在线观看免费午夜福利视频| 80岁老熟妇乱子伦牲交| 好看av亚洲va欧美ⅴa在| 国产av又大| 午夜a级毛片| netflix在线观看网站| 男人的好看免费观看在线视频 | 熟女少妇亚洲综合色aaa.| 久久久久久大精品| 天堂俺去俺来也www色官网| 久久久久国产一级毛片高清牌| 亚洲精品国产精品久久久不卡| 黄色成人免费大全| 欧美丝袜亚洲另类 | 女警被强在线播放| 午夜免费观看网址| 日韩国内少妇激情av| 久99久视频精品免费| 精品国产美女av久久久久小说| 咕卡用的链子| 久久久国产成人免费| 日韩欧美一区二区三区在线观看| 91九色精品人成在线观看| 真人一进一出gif抽搐免费| 一区二区日韩欧美中文字幕| 好看av亚洲va欧美ⅴa在| 国产熟女午夜一区二区三区| 婷婷六月久久综合丁香| 搡老熟女国产l中国老女人| 国产熟女午夜一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲自偷自拍图片 自拍| 一进一出抽搐动态| 午夜福利,免费看| 亚洲色图av天堂| 久久久久久亚洲精品国产蜜桃av| 亚洲av成人av| 麻豆成人av在线观看| 精品午夜福利视频在线观看一区| 淫秽高清视频在线观看|