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

    不同林齡油松和白扦林地植物多樣性分析

    2015-05-05 01:39:11李玉婷張建軍趙廷寧許宗文丁楊
    中國(guó)水土保持科學(xué) 2015年3期
    關(guān)鍵詞:油松林林齡草本

    李玉婷, 張建軍, 趙廷寧, 許宗文, 丁楊

    (北京林業(yè)大學(xué)水土保持學(xué)院,100083,北京)

    不同林齡油松和白扦林地植物多樣性分析

    李玉婷, 張建軍?, 趙廷寧, 許宗文, 丁楊

    (北京林業(yè)大學(xué)水土保持學(xué)院,100083,北京)

    為探明不同林齡的油松和白扦林地植物多樣性的變化規(guī)律并確定影響其變化的因子,以物種豐富度指數(shù)、物種多樣性指數(shù)和物種均勻度指數(shù)作為衡量植物多樣性的標(biāo)準(zhǔn),研究?jī)?nèi)蒙古興和縣二道營(yíng)流域友誼水庫(kù)不同林齡油松和白扦水土保持林地的生物多樣性特征及其影響因子。結(jié)果表明:1) 油松人工林植物多樣性指數(shù)大于白扦人工林,草本群落生物多樣性隨林齡的增加呈現(xiàn)出非線性增長(zhǎng)趨勢(shì)。5年生油松林的草本群落生物多樣性指數(shù)最高,為2.09;13年生油松最低,僅1.63。白扦人工林的Margalef指數(shù)高于油松林地,13年生白扦人工林的Margalef指數(shù)最高,為3.79;2年生油松的Margalef指數(shù)最低,僅為2.16。2) 在植被恢復(fù)演替進(jìn)程中,多年生草本逐漸替代一年生草本,演替朝正方向進(jìn)行。3)在眾多土壤因子中,土壤水、土壤抗沖性及土壤有機(jī)碳質(zhì)量分?jǐn)?shù)對(duì)生物多樣性指數(shù)起直接作用。研究結(jié)果可為森林草原過(guò)渡區(qū)域植被恢復(fù)與重建提供必要的理論支撐。

    Margalef植物多樣性; 林齡; 水土保持林; 通徑分析

    植物多樣性(plant diversity)是生物多樣性中以植物為主體,由植物、植物與環(huán)境之間所形成的復(fù)合體及其相關(guān)生態(tài)演替過(guò)程的總和[1],植物群落物種多樣性及其環(huán)境影響因素一直是群落生態(tài)學(xué)研究中的重點(diǎn)和熱點(diǎn),植物多樣性在諸多方面起著重要作用,如改變對(duì)水分、養(yǎng)分的利用效率[2],影響干擾發(fā)生的范圍和程度[3],以及對(duì)生態(tài)系統(tǒng)生產(chǎn)力性能增加效應(yīng)方面[4]。它不僅影響土地利用方式,還影響植被群落演替土壤物理性質(zhì)的組成及人工植被恢復(fù)的空間格局[5]。較高的生物多樣性能夠降低生態(tài)系統(tǒng)過(guò)程的變異性,而其受到土壤、水分、人類活動(dòng)等多方面因素的影響,若能夠通過(guò)人工植被恢復(fù)的方式增加森林草原的物種多樣性,則可促進(jìn)退化生態(tài)系統(tǒng)正方向演替[6]。

    植物多樣性的物種豐富度、多樣性和均勻程度可表征植物群落與不同自然地理及環(huán)境因素的相互關(guān)系[7],而探明森林草原過(guò)渡區(qū)植物的多樣性演替變化規(guī)律及其影響因素是首要解決的問(wèn)題;因此,筆者嘗試通過(guò)實(shí)地調(diào)查及數(shù)理分析等方法探明植物多樣性變化規(guī)律及其與影響因素間的關(guān)系。目前分析多樣性和環(huán)境因素之間相關(guān)性及關(guān)系的方法主要集中在線性回歸[8]、簡(jiǎn)單零階相關(guān)分析[9]等,而這些方法不能避免定性描述和簡(jiǎn)單相關(guān)分析的片面性,筆者運(yùn)用偏相關(guān)分析和通徑分析,在相關(guān)和回歸分析的基礎(chǔ)上,進(jìn)一步對(duì)回歸系數(shù)分割,計(jì)算環(huán)境影響因子對(duì)植物群落多樣性的直接和間接作用。

    實(shí)驗(yàn)區(qū)域設(shè)置在內(nèi)蒙古壩上高原的興和縣二道營(yíng)小流域,該區(qū)域?qū)儆谏植菰c荒漠草原的過(guò)渡區(qū),草本群落豐富,是具有代表性的荒漠草原植被帶。以油松(PinustabulaeformisCarr.)和白扦(PiceameyeriRehder.)水土保持樹(shù)種為研究對(duì)象,通過(guò)調(diào)查其林下的植被狀況及生物學(xué)特性,測(cè)定土壤的理化性質(zhì),分析比較不同林齡的油松與白扦的生物多樣性特征,以及土壤理化性質(zhì)與生物多樣性間的關(guān)系,探究人工植被恢復(fù)下的植物物種多樣性特征及其影響因素,旨在為恢復(fù)和保護(hù)干旱半干旱區(qū)域的人工植被、保育生物多樣性提供基礎(chǔ)資料,同時(shí)為植被恢復(fù)與經(jīng)營(yíng)提供一定的理論支撐。

    1 研究區(qū)概況

    興和縣隸屬于內(nèi)蒙古自治區(qū)烏蘭察布市,地處海河流域重要支流永定河源頭的內(nèi)蒙古壩上高原地區(qū),地理坐標(biāo)為E 113°21′09″~114°07′47″、N 40°26′40″~41°26′27″。總面積約3 518 km2,南北長(zhǎng)約109 km,南向北依次為侵蝕山地、平原和丘陵,海拔1 087~2 335 m,屬中溫帶大陸性季風(fēng)半干旱氣候。興和縣二道營(yíng)小流域面積為3.91 km2,年均氣溫4.2 ℃,年均降水量380 mm,降水主要集中在6—9月,多年平均徑流深50 mm,年均降水時(shí)間83 d,無(wú)霜期119 d。該流域種植的主要樹(shù)種有油松、白扦、沙棘(HippophaerhamoidesLinn.)、山桃(AmygdalusdavidianaCarr.)等,主要草本植物種有蒙古針茅(StipamongolorumTzvel.)、蒙古韭(AlliummongolicumRegel.)、長(zhǎng)芒草(StipabungeanaTrin.)、地椒(ThymusquinquecostatusCelak.)、三葉馬先蒿(PedicularisternataMaxim.)等。

    二道營(yíng)小流域土壤以干旱草原栗鈣土為主,兼有沙壤土帶。流域溝底為非地帶性土壤,由于雨水沖積作用,上游表土被沖刷留滯溝底,加之溝底種植的楊樹(shù)等植被枯枝落葉腐化,腐殖質(zhì)土層較厚。

    2 研究方法

    2.1 實(shí)驗(yàn)設(shè)計(jì)

    研究區(qū)以2、3、5、9、13年生5種樹(shù)齡作為設(shè)置樣地的標(biāo)準(zhǔn),選取立地條件相對(duì)一致的油松、白扦人工林樣地各15塊,共30塊樣地。野外調(diào)查時(shí)喬木層采用20 m×20 m樣地,每個(gè)樣地內(nèi)設(shè)置5個(gè)1 m×1 m的草本樣方。樣地調(diào)查因子包括植物的個(gè)體數(shù)(包括枯死草本數(shù))、林齡、種類、樹(shù)高、草本地徑及蓋度等,研究區(qū)的經(jīng)緯度、海拔及坡向通過(guò)GPS及羅盤(pán)測(cè)定。野外實(shí)驗(yàn)完成后帶回實(shí)驗(yàn)室進(jìn)行數(shù)據(jù)統(tǒng)計(jì)。樣地的基本概況如表1所示。

    在各樣地內(nèi)沿對(duì)角線選取3個(gè)采樣點(diǎn),用土鉆分別鉆取0~10、10~20和20~30 cm的土壤,將各樣點(diǎn)土壤按不同層次放置于處理臺(tái)混勻,按照四分法取300 g待測(cè)土樣,用鋁盒裝入土樣,貼上標(biāo)簽,用于測(cè)定土壤水分;以同樣采樣方式在采樣點(diǎn)挖取土壤剖面,用100 cm3環(huán)刀對(duì)土壤剖面分層進(jìn)行取土,貼上標(biāo)簽,裝入自封袋內(nèi)保存,帶回實(shí)驗(yàn)室用環(huán)刀浸水法對(duì)土壤密度、土壤孔隙度進(jìn)行分析測(cè)定;將測(cè)定過(guò)土壤水分的土樣去雜、磨細(xì),過(guò)1.00和0.25 mm土篩,測(cè)定土壤氮磷鉀等化學(xué)性質(zhì)。

    表1 樣地基本概況

    2.2 土壤樣品分析方法

    土壤水分采用烘干法進(jìn)行測(cè)定,土壤密度、土壤孔隙度通過(guò)環(huán)刀法測(cè)得;土壤抗沖性采用野外實(shí)地放水沖刷法測(cè)定[10-11],使用外加熱重鉻酸鉀氧化-容量法測(cè)定土壤有機(jī)質(zhì),靛酚藍(lán)比色法、酚二磺酸比色法測(cè)定土壤銨態(tài)氮及硝態(tài)氮,硫酸-高氯酸消煮-鉬銻抗比色法測(cè)定土壤全磷,氫氧化鈉熔融-火焰光度計(jì)法測(cè)定土壤全鉀[12]。

    2.3 物種多樣性分析方法

    采用4個(gè)研究多樣性的指標(biāo)對(duì)內(nèi)蒙古興和縣二道營(yíng)小流域的油松、白扦林下草本層植物多樣性進(jìn)行研究:物種豐富度指數(shù)(Margalef指數(shù))測(cè)定樣地內(nèi)物種豐富度程度,用物種均勻度指數(shù)(Pielou指數(shù))計(jì)測(cè)不同群落各個(gè)體數(shù)量分布是否均勻,以物種多樣性指數(shù)(Simpson指數(shù)、Shannon-Wiener指數(shù))衡量不同群落的物種多樣性[13-14]。

    1) Margalef指數(shù)

    2) Pielou指數(shù)

    3) Simpson指數(shù)

    4)Shannon-Wiener指數(shù)

    式中:S為樣方中出現(xiàn)的群落物種數(shù);n為觀察到的個(gè)體總數(shù);Ni為第i物種的個(gè)體數(shù);N為所有物種的個(gè)體數(shù);Pi為第i物種個(gè)體數(shù)占總個(gè)體數(shù)的比例。

    5) 草本重要值

    式中:HR為相對(duì)高度;FR為相對(duì)頻度;CR為相對(duì)蓋度。

    2.4 數(shù)據(jù)統(tǒng)計(jì)與分析

    實(shí)驗(yàn)數(shù)據(jù)利用Microsoft Excel 2007進(jìn)行整理,使用SPSS18.0對(duì)數(shù)據(jù)進(jìn)行相關(guān)分析和通徑分析,差異顯著性分析采用t檢驗(yàn)比較,并用Origin 8.0軟件對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行繪圖。

    3 結(jié)果與分析

    3.1 優(yōu)勢(shì)群落物種組成

    物種統(tǒng)計(jì)分析顯示,研究區(qū)造林后13年的演替過(guò)程中,調(diào)查樣方內(nèi)共出現(xiàn)草本植物42種,分屬18科37屬,主要有菊科、禾本科、蓼科和豆科。其中:油松林地的優(yōu)勢(shì)種主要有地椒、白草(PennisetumflaccidumGriseb.),分別占調(diào)查植被種的31%和28%;白扦林地的優(yōu)勢(shì)種主要有鐵桿蒿(ArtemisiasacrorumLedeb.)、金色狗尾草(SetariapumilaRoem.),分別占調(diào)查樣地植被種的27%和30%。研究區(qū)的主要物種中有5個(gè)科屬于西北干旱荒漠植物區(qū)系和荒漠植被中的優(yōu)勢(shì)科。相關(guān)研究[15-16]表明,此類科屬植物具有廣布性,使該地區(qū)物種能以其龐大的種系和適應(yīng)能力在該生境取得優(yōu)勢(shì),成為荒漠草原的重要建造者。本研究根據(jù)物種重要值劃定了不同林齡水土保持林對(duì)應(yīng)的群落類型(表2),可以看出,造林13年中,不同種植年限的植物群落經(jīng)歷了地椒、白草、金色狗尾草、鐵桿蒿、三葉馬先蒿5種類型的演替過(guò)程,群落結(jié)構(gòu)隨著植被類型的演變發(fā)生了顯著變化,物種數(shù)也隨演替時(shí)間呈現(xiàn)出先增大后減少的單峰變化趨勢(shì)。

    表2 不同林齡水土保持林地草本群落特征

    注Note:白草PennisetumflaccidumGriseb.,地椒ThymusquinquecostatusCelak.,金色狗尾草SetariapumilaRoem.,三葉馬先蒿PedicularisternataMaxim.,蒙古針茅StipamongolorumTzvel.,鐵桿蒿ArtemisiasacrorumLedeb.,香青蘭DracocephalummoldavicaLinn.,馬先蒿PedicularisfloribundaFranch,蒙古馬蘭AstermongolicusFranch.,狼毒StellerachamaejasmeLinn.,芨芨草Achnatherumsplendens(Trin.) Nevski.,冰草AgropyroncristatumP. Beauv.,蘭花CymbidiumvirescensLindl.,扁穗雀麥BromuscatharticusVahl.,黃花蒿ArtemisiaannuaL.,狼尾草PennisetumalopeuroidesSpreng.,臭蒿ArtemisiahediniiOstenf.,長(zhǎng)芒草StipabungeanaTrin.,韭白AlliummacrostemonBunge.,胡盧巴SophoraalopecuroidesL.,刺兒菜CirsiumjaponicumFisch. ex DC.。

    3.2 生物多樣性指數(shù)特征

    3.2.1 Margalef指數(shù)變化 作為定量評(píng)價(jià)衡量生物多樣性的標(biāo)準(zhǔn),生物多樣性指數(shù)的計(jì)算可直觀反映研究區(qū)的生境狀況與質(zhì)量。圖1表示了在不同種植年限下2種水土保持林的生物多樣性指數(shù)變化??梢钥闯鲈炝趾笱芯繀^(qū)的物種豐富度有明顯的升高趨勢(shì),油松和白扦林地的Margalef指數(shù)隨種植年限皆呈增長(zhǎng)趨勢(shì)。造林13年以來(lái),研究區(qū)水土保持林的物種豐富度指數(shù)表現(xiàn)為白扦林>油松林,在造林初期,對(duì)比油松和白扦林地豐富度指數(shù)變化不明顯,造林第5年后,豐富度指數(shù)顯著增高。在造林期間油松林的草本群落豐富度指數(shù)普遍低于比白扦林,造林2~9年間物種豐富度增長(zhǎng)緩慢,而在13年時(shí)增幅最為明顯,達(dá)到最大值。白扦林地草本群落豐富度指數(shù)在造林第2~3年間,沒(méi)有明顯增減變化,而在造林第5年后,豐富度指數(shù)呈現(xiàn)急劇上升的趨勢(shì),造林第9年,增幅變緩,直至13年時(shí)達(dá)到峰值。總體上,白扦林的草本群落豐富度指數(shù)變化程度高于油松林地。

    圖1 不同林齡水土保持林的物種豐富度指數(shù)Fig.1 Margalef index of soil and water conservation forest at different stand ages

    3.2.2 Pielou指數(shù)變化 由圖2可以看出,白扦林地的均勻度指數(shù)隨著林齡增加呈現(xiàn)先增加后減小的

    單峰趨勢(shì),而油松林地的均勻度指數(shù)隨林齡的增加而增加,油松、白扦林的均勻度指數(shù)在造林近第9年時(shí)相同。造林初期,白扦的均勻度指數(shù)高于油松林地,這種情況一直持續(xù)到第9年,隨后油松林地草本群落的均勻度指數(shù)持續(xù)升高,高于白扦林地。造林第3年開(kāi)始,白扦林地草本群落的均勻度指數(shù)下降的主要原因可能由于白扦作為淺根性樹(shù)種且喜濕潤(rùn)環(huán)境,其林下的土壤水分及養(yǎng)分主要供給白扦林的生長(zhǎng),造成草本群落與喬木樹(shù)種的營(yíng)養(yǎng)爭(zhēng)奪及串根現(xiàn)象,不能適應(yīng)該立地環(huán)境的草本群落被更高級(jí)的鐵桿蒿、馬先蒿群落所替代(表2),導(dǎo)致白扦林下的草本群落物種均勻度指數(shù)呈下降趨勢(shì)。而油松樹(shù)種根深且蒸騰強(qiáng)度低,造成其生長(zhǎng)期間均勻度指數(shù)均呈現(xiàn)出緩慢上升的趨勢(shì),在第13年達(dá)到峰值0.746。

    圖2 不同林齡水土保持林的均勻度指數(shù)Fig.2 Pielou index of soil and water conservation forest at different stand ages

    3.2.3 Simpson指數(shù)變化 圖3顯示,2種水土保持林的Simpson指數(shù)表現(xiàn)為油松>白扦,但白扦和油松在不同林齡階段的Simpson指數(shù)變化趨勢(shì)一致,即各自呈單峰曲線變化。 其中,油松林地的草本群落Simpson指數(shù)在第9年達(dá)到最大,為0.825。白扦林地的草本群落物種多樣性指數(shù)在第5年達(dá)到峰值,為0.778。造林2~13年間,白扦和油松林的草本群落Simpson指數(shù)先逐年增大,后出現(xiàn)回落。這是因?yàn)樵谠炝值?年后,出現(xiàn)香青蘭、金色狗尾草、鐵桿蒿等5種開(kāi)始替代造林前較低等物種,群落內(nèi)物種種類在前5年逐年增加,隨演替時(shí)間的繼續(xù),馬先蒿等優(yōu)勢(shì)種侵入,物種種類繼續(xù)增加,分布也較均勻。進(jìn)入第9年后,由于多年生草本的侵入,數(shù)量不斷增加(表2),物種間的資源競(jìng)爭(zhēng)逐漸激烈,使1年生及多年生適應(yīng)性差的草本漸漸退出演替過(guò)程,Simpson指數(shù)不斷下降。

    圖3 不同林齡水土保持林的Simpson指數(shù)Fig.3 Simpson index of soil and water conservation forest at different stand ages

    3.2.4 物種多樣性指數(shù)變化 不同植被恢復(fù)模式的生物多樣性特征實(shí)際與其群落結(jié)構(gòu)和演替過(guò)程相對(duì)應(yīng),Shannon-Wiener指數(shù)表征局域生境的植物多樣性,反映物種水平上種群的多樣性和異質(zhì)性,可以全面衡量種群的物種多樣性[17]。由圖4可以看出:在造林前3年,油松林地的草本群落多樣性指數(shù)急劇上升,而由于淺根系的白扦樹(shù)種處于生長(zhǎng)期,與其周圍的草本發(fā)生營(yíng)養(yǎng)爭(zhēng)奪,造成白扦林地的多樣性指數(shù)小幅度降低;隨著演替的進(jìn)行,狼毒等偶見(jiàn)種侵入,造林第5年后,2種林地的多樣性指數(shù)皆逐年下降,尤其油松林地的多樣性指數(shù)下降程度高。群落演替第9年后,不斷出現(xiàn)了多年生菊科植物替代一年生藜科等草本的現(xiàn)象,油松林地的多樣性指數(shù)相對(duì)白扦林地下降幅度較大。

    圖4 不同林齡水土保持林的Shannon-Wiener指數(shù)Fig.4 Shannon-Wiener index of soil and water conservation forest at different stand ages

    3.2.5 生物多樣性指數(shù)函數(shù)關(guān)系 回歸分析表明,造林13年間,Simpson指數(shù)、Margalef指數(shù)、Shannon-Wiener指數(shù)和Pielou指數(shù)均呈二次函數(shù)變化規(guī)律:y=ax2+bx+c。從表3可以看出,除了Margalef指數(shù)與林齡之間呈正向曲線遞增外,Pielou指數(shù)、Simpson指數(shù)、Shannon-Wiener指數(shù)都隨林齡增大呈單峰或存在一定峰值的變化。Margalef指數(shù)在第9年達(dá)到演替序列峰值,而Simpson指數(shù)、Shannon-Wiener指數(shù)分別在第9年和第5年達(dá)到最大值。此后,地帶性次生群落逐漸演替而成,此時(shí)的多樣性指數(shù)幅度呈下降趨勢(shì),群落演替變化緩慢,達(dá)到相對(duì)穩(wěn)定狀態(tài)。

    3.3 土壤理化性質(zhì)與植物多樣性間的相關(guān)分析

    由于此生長(zhǎng)階段的油松和白扦林皆未郁閉,光照并非植物多樣性大小的主要限制因素;故此時(shí)不考慮光照條件對(duì)植物多樣性的影響,而土壤因子在植被演替過(guò)程中始終起著重要作用[18]。為探明土壤因子對(duì)植物多樣性指數(shù)的影響程度,本實(shí)驗(yàn)通過(guò)逐步回歸選出了土壤水分X1、土壤密度X2、土壤抗沖性X3、有機(jī)碳質(zhì)量分?jǐn)?shù)X4、全磷質(zhì)量分?jǐn)?shù)X5土壤理化因子與多樣性指數(shù)進(jìn)行偏相關(guān)與通徑分析。由偏相關(guān)分析結(jié)果可看出(表4),土壤含水量、土壤密度和土壤抗沖性與Margalef指數(shù)達(dá)到極顯著水平(P<0.01),Pielou指數(shù)與土壤水分和全磷質(zhì)量分?jǐn)?shù)呈極顯著正相關(guān),Simpson指數(shù)與土壤水分、有機(jī)碳質(zhì)量分?jǐn)?shù)及土壤密度間達(dá)到顯著水平(P<0.05)。Shannon-Wiener指數(shù)與土壤水分、土壤抗沖性間呈極顯著正相關(guān),全磷質(zhì)量分?jǐn)?shù)與Shannon-Wiener指數(shù)間達(dá)顯著水平。通過(guò)通徑分析,從直接作用的絕對(duì)值大小(表5)看:土壤因子對(duì)Margalef指數(shù)的影響順序?yàn)閄1>X2>X3,表明土壤因子中對(duì)多樣性指數(shù)起直接主要作用的是X1和X2,且X1通過(guò)X2、X3分別對(duì)Margalef指數(shù)產(chǎn)生較大的間接作用;X4對(duì)Simpson指數(shù)起直接作用,且通過(guò)X1對(duì)其起較大間接作用;X1對(duì)Pielou指數(shù)的直接及間接作用皆大于X5。相關(guān)研究[19-20]認(rèn)為在局部范圍內(nèi)的有機(jī)碳質(zhì)量分?jǐn)?shù)X3和植物多樣性呈負(fù)相關(guān)。本研究得出Shannon-Wiener指數(shù)的直接作用大小順序?yàn)閄3>X1>X5,且X3通過(guò)X5對(duì)其起較大間接作用。雖然土壤水分對(duì)多樣性指數(shù)增減起重要作用,但它同時(shí)通過(guò)有機(jī)碳質(zhì)量分?jǐn)?shù)和土壤抗沖性能作用于植物多樣性。盡管M.L. Rosenzweig[20]認(rèn)為該土壤因子對(duì)植物β多樣性有促進(jìn)作用,但由于研究地域的差異,導(dǎo)致土壤因子對(duì)多樣性指數(shù)影響的研究結(jié)果有不同程度的差異。

    表3 林齡-草本群落生物多樣性指數(shù)的回歸方程

    注:M為Margalef指數(shù),E為Pielou指數(shù),D為Simpson指數(shù),H為Shannon-Wiener指數(shù)。Note:Mis Margalef index,Eis Pielou index,Dis Simpson index,His Shannon-Wiener index.

    表4 土壤理化性質(zhì)與生物多樣性指數(shù)的偏相關(guān)性分析

    Tab.4 Partial correlation analysis between soil physico-chemical properties and biodiversity index

    項(xiàng)目Item偏相關(guān)PartialCorrelationT檢驗(yàn)值T-test顯著水平Significancelevelr(M,X1)0.997313.6870.0001r(M,X2)-0.993622.3930.0013r(M,X3)0.956432.5770.0022r(E,X1)0.952911.8630.0015r(E,X5)0.93368.4580.0041r(D,X1)0.993038.4500.0007r(D,X2)-0.823535.5440.0043r(D,X4)0.802614.8730.0076r(H,X1)0.998416.6590.0013r(H,X3)0.952925.1910.0005r(H,X5)0.792210.7450.0407

    表5 多樣性指數(shù)與土壤因子的通徑分析

    4 結(jié)論

    1) 實(shí)驗(yàn)區(qū)域植物群落結(jié)構(gòu)隨著植被類型的演變發(fā)生了顯著變化,物種數(shù)隨林齡增加呈現(xiàn)單峰變化趨勢(shì),這主要由油松和白扦林地優(yōu)勢(shì)種的生物學(xué)特性決定。

    2) 由于2種水土保持樹(shù)種的生理特性,對(duì)比得出不同林齡林地的草本群落Simpson指數(shù)表現(xiàn)為油松林>白扦林,而Margalef豐富度指數(shù)表現(xiàn)為白扦林>油松林。在多年生草本代替一年生草本過(guò)程中,出現(xiàn)有毒雜草,油松和白扦林地的Shannon-Wiener指數(shù)在第5年便開(kāi)始轉(zhuǎn)為下降趨勢(shì)。演替過(guò)程中,2種林地的群落生物多樣性指數(shù)先增大后減小,總體呈二次函數(shù)變化規(guī)律:y=ax2+bx+c。

    3) 簡(jiǎn)單相關(guān)分析表明,土壤含水量與Pielou指數(shù)、Simpson指數(shù)、Shannon-Wiener指數(shù)呈極顯著正相關(guān)(P<0.01),土壤有機(jī)碳質(zhì)量分?jǐn)?shù)與Shannon-wiener指數(shù)呈顯著正相關(guān)(P<0.05)。通過(guò)偏相關(guān)及通徑分析,研究結(jié)果表明:土壤因子對(duì)Margalef指數(shù)影響順序?yàn)閄1>X2>X3,表明土壤因子中對(duì)多樣性指數(shù)起直接主要作用的是X1和X2,且X1通過(guò)X2、X3分別對(duì)Margalef指數(shù)產(chǎn)生較大的間接作用;X4對(duì)Simpson指數(shù)起直接作用;X1對(duì)Pielou指數(shù)的直接及間接作用皆大于X5,對(duì)Shannon-Wiener指數(shù)的直接作用大小順序?yàn)閄3>X1>X5,且X3通過(guò)X5對(duì)其起較大間接作用。

    綜上所述,加強(qiáng)油松等具有耐寒旱、蒸騰耗水少且改良土壤效益高等特性的水土保持林種或灌木種的種植力度,有助于提高內(nèi)蒙古高原地區(qū)植被恢復(fù)的植物多樣性。在不同林齡油松和白扦人工林的生長(zhǎng)階段,本實(shí)驗(yàn)雖在土壤理化性質(zhì)與植物多樣性的相關(guān)分析方面做了研究,但土壤母質(zhì)以及土壤微環(huán)境因子等也會(huì)影響該地區(qū)的生物多樣性;因此,所得土壤理化性質(zhì)與生物多樣性關(guān)系的結(jié)果有待在時(shí)空范圍內(nèi)進(jìn)一步研究探索。

    在野外調(diào)查階段,得到了內(nèi)蒙古興和縣林業(yè)局局長(zhǎng)武海龍和北京林業(yè)大學(xué)水土保持學(xué)院席沁及王新的幫助和支持,在此表示感謝!

    [1] 馬克平.試論生物多樣性的概念[J].生物多樣性,1993,1(1):20-22

    [2] Collins S L, Knapp A K, Briggs J M, et al. Modulation of diversity by grazing and mowing in native tall grass prairie [J].Science,1998,280:745-747

    [3] Chapin Ⅲ F S, Walker B H, Hobbs R J, et al. Biotic control over the functioning of ecosystems [J]. Science, 1997,277(4):500-503

    [4] Yachi S, Loreau M.Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis[J]. Proceedings of the National Academy of Sciences USA, 1999,96(4):1463-1468

    [5] 周萍,劉國(guó)彬,侯喜祿.黃土丘陵區(qū)不同坡向及坡位草本群落生物量及多樣性研究[J].中國(guó)水土保持科學(xué),2009,7(1):67-73

    [6] 張克斌,李瑞,侯瑞萍,等. 寧夏鹽池縣不同荒漠化治理措施植物多樣性研究[J].中國(guó)水土保持科學(xué),2004,2(4):66-71

    [7] Hector A, Schmid B, Beierkuhnlein C, et al. Plant diversity and productivity experiments in European grasslands[J]. Science, 1999, 286(5442):1123-1130

    [8] 陳靈芝,馬克平.生物多樣性研究的原理與方法[M].北京:中國(guó)科學(xué)技術(shù)出版社,1995:12-18

    [9] 后源,郭正剛,龍瑞軍.黃河首曲濕地退化過(guò)程中植物群落組分及物種多樣性的變化[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(1):27-32

    [10] 張建軍,張寶穎,畢華興,等.黃土區(qū)不同植被條件下的土壤抗沖性[J].北京林業(yè)大學(xué)學(xué)報(bào),2004,26(6):25-28

    [11] 李玉婷,張建軍,茹豪,等.晉西黃土區(qū)不同土地利用方式對(duì)土壤抗沖性及糙率的影響[J].水土保持學(xué)報(bào),2013,27(4):1-7

    [12] 鮑士旦.土壤農(nóng)化分析[M].北京:中國(guó)農(nóng)業(yè)出版社,2000:15-98

    [13] 馬克平.生物群落多樣性的測(cè)度方法[M].北京:中國(guó)科學(xué)技術(shù)出版社,1995:141-165

    [14] 中國(guó)科學(xué)院內(nèi)蒙古寧夏綜合考察隊(duì).內(nèi)蒙古植被[M].北京:科學(xué)出版社,1985:15-27

    [15] 黨榮理,潘曉玲.西北干旱荒漠區(qū)種子植物科的區(qū)系分析[J].西北植物學(xué)報(bào),2002,22(1):24-32

    [16] 秦偉,朱清科,張宇清,等. 陜北黃土區(qū)生態(tài)修復(fù)過(guò)程中植物群落物種多樣性變化[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(2):403-409

    [17] 陳芳.莽山植物物種多樣性及其土壤理化性質(zhì)關(guān)系的研究[D].長(zhǎng)沙:湖南農(nóng)業(yè)大學(xué),2011

    [18] 余悅.黃河三角洲原生掩體中土壤微生物多樣性及其土壤理化性質(zhì)關(guān)系[D].山東泰安:山東農(nóng)業(yè)大學(xué),2012:72-81

    [19] 王凱博,陳美玲,秦娟,等.子午嶺植被自然演替中植物多樣性變化及其與土壤理化性質(zhì)的關(guān)系[J].西北植物學(xué)報(bào),2007,27(10):2089-2096

    [20] Rosenzweig M L. Species diversity in space and time[M]. Cambridge: Cambridge University Press, 1995: 131-159

    (責(zé)任編輯:程 云)

    Plant diversity ofPinustabulaeformisandPiceameyeriat different stand ages

    Li Yuting, Zhang Jianjun, Zhao Tingning, Xu Zongwen, Ding Yang

    (School of Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China)

    In order to explore the change law of plant diversity inPinustabulaeformisandPiceameyeristands at different stand ages and figure out the driving factors for plant diversity, we took Margalef index, Simpson index, Shannon-Wiener index and Pielou index as criteria for measuring biodiversity. We analyzed and compared the characteristics of biodiversity and soil physico-chemical properties ofPinustabulaeformisstands andPiceameyeristands at different stand ages in Youyi Reservoir, Erdaoying Basin, Xinghe County of Inner Mongolia. The results indicated that: 1) The biodiversity index ofPinustabulaeformisartificial forest was greater than that ofPiceameyeriartificial forest. The biodiversity of herbaceous community presented a nonlinear rising tendency with the increase of forest stand age. Five-year-oldPinustabulaeformisforest had the highest biodiversity index of 2.09, whereas the biodiversity index of the thirteen-year-old forest was the lowest, only 1.63. The Margalef index ofPiceameyeriartificial forest was higher than that ofPinustabulaeformisforest, 3.79 in thirteen-year-oldPiceameyeriand 2.16 in two-year-oldPinustabulaeformisartificial forest. 2) During the process of vegetation restoration succession, perennial herbs replaced annual herbs gradually, ongoing towards a positive direction. 3) Among various edaphic factors, the content of soil water, soil anti-scouring ability and SOC had direct effects on biodiversity index. The results of our study could provide a theoretical support for vegetation restoration and rehabilitation in the forest-steppe ecotone.

    Margalef plant diversity; stand age; soil and water conservation forest; path analysis

    2014-03-20

    2015-02-15

    李玉婷(1990—),女,碩士研究生。主要研究方向:退耕還林生態(tài)效益評(píng)價(jià)。E-mail:lyt0719@sina.cn

    ?通信作者簡(jiǎn)介: 張建軍(1964—),男,教授。主要研究方向:水土保持與森林水文。E-mail:zhangjianjun@bjfu.edu.cn

    Q948.113

    A

    1672-3007(2015)03-0097-07

    項(xiàng)目名稱: 國(guó)家林業(yè)公益性行業(yè)專項(xiàng)“天然林保護(hù)等林業(yè)工程生態(tài)效益評(píng)價(jià)研究”(201304308)

    猜你喜歡
    油松林林齡草本
    毛鋪草本酒
    毛鋪草本酒
    檫樹(shù)優(yōu)樹(shù)半同胞子代測(cè)定林樹(shù)高性狀遺傳變異研究
    彩墨繪草本
    北方油松林病蟲(chóng)害預(yù)防及措施研究
    世界家苑(2018年9期)2018-09-18 07:34:38
    建平縣規(guī)模型油松林近自然化撫育改造技術(shù)探究
    山東林業(yè)科技(2017年1期)2017-06-29 07:53:50
    遼西半干旱地區(qū)油松林氣候區(qū)劃的研究
    不同林齡紅松人工林優(yōu)樹(shù)選擇技術(shù)的研究
    遼東山區(qū)不同林齡落葉松人工林土壤理化特征
    嫩草影院入口| 午夜福利高清视频| 免费观看av网站的网址| 国产毛片a区久久久久| 成人毛片60女人毛片免费| videos熟女内射| 欧美高清成人免费视频www| 看免费成人av毛片| 精品国产一区二区三区久久久樱花 | 女的被弄到高潮叫床怎么办| 干丝袜人妻中文字幕| 成年免费大片在线观看| 成年av动漫网址| 国产综合懂色| 国产高清三级在线| a级一级毛片免费在线观看| 国产探花极品一区二区| 高清午夜精品一区二区三区| 身体一侧抽搐| 亚洲精品亚洲一区二区| 哪个播放器可以免费观看大片| 18禁裸乳无遮挡动漫免费视频 | 亚洲婷婷狠狠爱综合网| 久久韩国三级中文字幕| 国模一区二区三区四区视频| 国产精品99久久99久久久不卡 | 亚洲最大成人av| 综合色丁香网| 欧美日韩国产mv在线观看视频 | 亚洲国产精品成人久久小说| 熟女人妻精品中文字幕| 麻豆乱淫一区二区| 欧美xxxx性猛交bbbb| 一本久久精品| 国产视频首页在线观看| 日产精品乱码卡一卡2卡三| 国产极品天堂在线| av国产免费在线观看| 午夜日本视频在线| 精品国产乱码久久久久久小说| 一本色道久久久久久精品综合| 日韩伦理黄色片| 国产精品av视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 秋霞伦理黄片| 日日啪夜夜爽| 尾随美女入室| 亚洲性久久影院| 成人高潮视频无遮挡免费网站| 亚洲av.av天堂| 爱豆传媒免费全集在线观看| videossex国产| 中文字幕制服av| 国产成人一区二区在线| 免费高清在线观看视频在线观看| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 日韩欧美精品免费久久| 99热6这里只有精品| 免费黄网站久久成人精品| 免费人成在线观看视频色| www.av在线官网国产| 欧美日韩一区二区视频在线观看视频在线 | 久久久精品欧美日韩精品| 亚洲aⅴ乱码一区二区在线播放| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜爱| 久久久久精品性色| 国产v大片淫在线免费观看| 国产大屁股一区二区在线视频| 18+在线观看网站| 日韩视频在线欧美| 麻豆成人av视频| 九草在线视频观看| 免费观看av网站的网址| 91久久精品电影网| 国内揄拍国产精品人妻在线| 噜噜噜噜噜久久久久久91| 日韩强制内射视频| 91久久精品国产一区二区成人| 日本欧美国产在线视频| 亚洲四区av| 女的被弄到高潮叫床怎么办| 少妇的逼水好多| 91精品一卡2卡3卡4卡| 欧美区成人在线视频| 久久久午夜欧美精品| 亚洲国产精品专区欧美| 国产高清国产精品国产三级 | 精华霜和精华液先用哪个| 国产在视频线精品| 99热6这里只有精品| 久久久色成人| 男人和女人高潮做爰伦理| 女的被弄到高潮叫床怎么办| 国产又色又爽无遮挡免| 午夜激情福利司机影院| 最近2019中文字幕mv第一页| 日韩一区二区三区影片| 天堂俺去俺来也www色官网| 国产综合懂色| 色综合色国产| 男女啪啪激烈高潮av片| 一级毛片 在线播放| 免费大片18禁| 欧美xxxx黑人xx丫x性爽| 建设人人有责人人尽责人人享有的 | 色综合色国产| 91午夜精品亚洲一区二区三区| 欧美日韩精品成人综合77777| av免费在线看不卡| 偷拍熟女少妇极品色| 国产精品熟女久久久久浪| 色婷婷久久久亚洲欧美| 97在线人人人人妻| 免费看光身美女| 成人毛片a级毛片在线播放| 国产精品av视频在线免费观看| 欧美性猛交╳xxx乱大交人| 综合色av麻豆| 久久人人爽人人爽人人片va| 综合色丁香网| 亚洲,一卡二卡三卡| 成人无遮挡网站| 五月玫瑰六月丁香| 亚洲欧美成人精品一区二区| 美女被艹到高潮喷水动态| 精品一区二区三卡| 免费观看的影片在线观看| 麻豆国产97在线/欧美| 九色成人免费人妻av| 亚洲欧美日韩无卡精品| 国产精品国产av在线观看| 久久精品国产鲁丝片午夜精品| 丝瓜视频免费看黄片| 麻豆乱淫一区二区| 亚洲精品456在线播放app| 亚洲美女搞黄在线观看| 亚洲av成人精品一二三区| 夫妻性生交免费视频一级片| 亚洲图色成人| 日韩欧美精品免费久久| 免费高清在线观看视频在线观看| 看非洲黑人一级黄片| 欧美丝袜亚洲另类| 新久久久久国产一级毛片| 国产久久久一区二区三区| 能在线免费看毛片的网站| 一个人观看的视频www高清免费观看| 亚洲精品日本国产第一区| 国产精品99久久99久久久不卡 | 身体一侧抽搐| 又黄又爽又刺激的免费视频.| 国产精品爽爽va在线观看网站| 成年版毛片免费区| 三级国产精品片| 97超视频在线观看视频| 日韩精品有码人妻一区| 伊人久久精品亚洲午夜| 国产欧美亚洲国产| 国产黄片美女视频| 黄色欧美视频在线观看| 禁无遮挡网站| 成年人午夜在线观看视频| 2021少妇久久久久久久久久久| 国产成人a∨麻豆精品| 黄色怎么调成土黄色| 在线观看国产h片| 久久久久国产精品人妻一区二区| 在线 av 中文字幕| 国产淫片久久久久久久久| 亚洲高清免费不卡视频| 女人被狂操c到高潮| 人妻 亚洲 视频| 一个人看的www免费观看视频| 3wmmmm亚洲av在线观看| 91久久精品国产一区二区三区| 嫩草影院入口| 久久精品久久久久久久性| 国产 一区 欧美 日韩| 久久久久国产网址| 神马国产精品三级电影在线观看| 中文精品一卡2卡3卡4更新| 久久精品久久精品一区二区三区| 人妻系列 视频| 秋霞伦理黄片| 亚洲av日韩在线播放| 视频区图区小说| 精品国产三级普通话版| 国产精品一及| 久久精品久久久久久噜噜老黄| 26uuu在线亚洲综合色| 国产精品人妻久久久影院| 男男h啪啪无遮挡| 麻豆精品久久久久久蜜桃| 成人欧美大片| 欧美日韩在线观看h| 2018国产大陆天天弄谢| 禁无遮挡网站| 精品久久久噜噜| 亚洲欧美日韩东京热| 日日啪夜夜撸| av一本久久久久| 在线 av 中文字幕| 国产精品麻豆人妻色哟哟久久| 少妇的逼水好多| 一级毛片aaaaaa免费看小| 三级经典国产精品| av播播在线观看一区| 狂野欧美激情性xxxx在线观看| 97热精品久久久久久| 日韩亚洲欧美综合| 精品少妇黑人巨大在线播放| 国产 一区精品| 久久99蜜桃精品久久| 欧美精品国产亚洲| 一个人观看的视频www高清免费观看| 国内精品宾馆在线| 欧美97在线视频| 成人二区视频| 欧美+日韩+精品| 爱豆传媒免费全集在线观看| 亚洲国产精品成人综合色| 成人欧美大片| 日韩av在线免费看完整版不卡| 国产视频首页在线观看| 美女内射精品一级片tv| 内地一区二区视频在线| 黄色一级大片看看| 亚洲欧洲国产日韩| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 欧美高清成人免费视频www| 91午夜精品亚洲一区二区三区| 有码 亚洲区| 国产精品99久久99久久久不卡 | 亚洲精品乱码久久久v下载方式| 插逼视频在线观看| 99九九线精品视频在线观看视频| 亚洲在线观看片| 我的老师免费观看完整版| 日日摸夜夜添夜夜爱| 80岁老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 天堂俺去俺来也www色官网| 精品熟女少妇av免费看| 最近中文字幕高清免费大全6| 亚洲在线观看片| 久久精品国产自在天天线| 亚洲国产精品成人久久小说| 国产精品一区二区在线观看99| 搡老乐熟女国产| 亚洲天堂国产精品一区在线| 中国美白少妇内射xxxbb| 又大又黄又爽视频免费| 大片免费播放器 马上看| 最近最新中文字幕大全电影3| 日本黄大片高清| 国产成人福利小说| 18禁在线无遮挡免费观看视频| 精品亚洲乱码少妇综合久久| 91精品一卡2卡3卡4卡| 日本色播在线视频| av国产精品久久久久影院| 亚洲人成网站在线播| 美女脱内裤让男人舔精品视频| 一区二区av电影网| 亚洲精品影视一区二区三区av| 中文字幕亚洲精品专区| 国产男女超爽视频在线观看| 美女国产视频在线观看| 亚洲四区av| 成年女人在线观看亚洲视频 | 亚洲精品色激情综合| 久久97久久精品| 国产爽快片一区二区三区| 国产黄片美女视频| av免费观看日本| 久久久国产一区二区| 精品一区二区三区视频在线| 秋霞伦理黄片| 国产一区二区三区av在线| 毛片一级片免费看久久久久| 91精品国产九色| 国产男女超爽视频在线观看| xxx大片免费视频| 亚洲国产精品成人久久小说| 人人妻人人爽人人添夜夜欢视频 | 我的女老师完整版在线观看| 听说在线观看完整版免费高清| 国产精品伦人一区二区| 人人妻人人澡人人爽人人夜夜| 99久久精品热视频| 五月天丁香电影| 亚洲国产精品999| 亚洲天堂av无毛| 性插视频无遮挡在线免费观看| 热re99久久精品国产66热6| 日韩电影二区| 精品久久久久久久久亚洲| 嫩草影院入口| 青春草视频在线免费观看| 亚洲精品日韩在线中文字幕| 丝袜美腿在线中文| 免费黄频网站在线观看国产| 中文天堂在线官网| 久久综合国产亚洲精品| 新久久久久国产一级毛片| 日韩av不卡免费在线播放| 综合色丁香网| 亚洲图色成人| 身体一侧抽搐| 精品午夜福利在线看| 99久久精品国产国产毛片| 国产午夜精品一二区理论片| 亚洲天堂av无毛| 日韩免费高清中文字幕av| 亚洲人成网站在线观看播放| 制服丝袜香蕉在线| 美女xxoo啪啪120秒动态图| 中文字幕免费在线视频6| 夫妻性生交免费视频一级片| 中文精品一卡2卡3卡4更新| 日韩一区二区三区影片| av在线app专区| 亚洲成色77777| 最近最新中文字幕免费大全7| 六月丁香七月| 久久久精品免费免费高清| 国产一区二区三区av在线| 午夜老司机福利剧场| 秋霞伦理黄片| 欧美三级亚洲精品| 欧美性感艳星| 国产精品99久久99久久久不卡 | 成年免费大片在线观看| 色综合色国产| 2018国产大陆天天弄谢| 精品一区二区三区视频在线| 国产成人freesex在线| 亚洲国产精品专区欧美| 久久精品国产亚洲网站| 国产精品一区二区在线观看99| 搡女人真爽免费视频火全软件| 老司机影院成人| 97超碰精品成人国产| 99久久精品热视频| 97精品久久久久久久久久精品| 色综合色国产| 麻豆国产97在线/欧美| 中文字幕制服av| 麻豆久久精品国产亚洲av| 校园人妻丝袜中文字幕| 新久久久久国产一级毛片| 午夜福利在线在线| 深爱激情五月婷婷| av在线app专区| 国产黄频视频在线观看| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 久热久热在线精品观看| 亚洲无线观看免费| av网站免费在线观看视频| 久久热精品热| 自拍欧美九色日韩亚洲蝌蚪91 | av专区在线播放| 1000部很黄的大片| 欧美日韩亚洲高清精品| 亚洲人成网站高清观看| 三级经典国产精品| 日本午夜av视频| 偷拍熟女少妇极品色| 最近中文字幕高清免费大全6| 日韩一区二区三区影片| 成人亚洲欧美一区二区av| 久久久精品94久久精品| 亚洲成色77777| 亚洲av中文av极速乱| 亚洲国产最新在线播放| 蜜桃久久精品国产亚洲av| videossex国产| 男人添女人高潮全过程视频| 亚洲欧美一区二区三区黑人 | 中国三级夫妇交换| 久久99热这里只频精品6学生| 丰满乱子伦码专区| 一区二区三区乱码不卡18| 在线播放无遮挡| 亚洲一级一片aⅴ在线观看| 国产老妇伦熟女老妇高清| 国产探花极品一区二区| 免费不卡的大黄色大毛片视频在线观看| 久久99热6这里只有精品| .国产精品久久| 丝袜美腿在线中文| 97超碰精品成人国产| 日韩av免费高清视频| 精品人妻一区二区三区麻豆| 中文欧美无线码| 午夜爱爱视频在线播放| 亚洲成色77777| 深爱激情五月婷婷| 大又大粗又爽又黄少妇毛片口| 日本wwww免费看| 久久99热6这里只有精品| 国产淫语在线视频| 亚洲av中文av极速乱| 色婷婷久久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲av嫩草精品影院| 欧美老熟妇乱子伦牲交| 久久久国产一区二区| 国产乱人偷精品视频| 日韩伦理黄色片| 青春草国产在线视频| 国产免费一级a男人的天堂| 最近手机中文字幕大全| 精品酒店卫生间| 精品久久国产蜜桃| 精品久久久噜噜| 大码成人一级视频| 国产综合懂色| av免费在线看不卡| 直男gayav资源| 视频区图区小说| 国产一级毛片在线| 午夜福利在线在线| 亚洲国产精品成人久久小说| 亚洲经典国产精华液单| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 波多野结衣巨乳人妻| 国产成人免费观看mmmm| 少妇人妻精品综合一区二区| 五月开心婷婷网| 777米奇影视久久| 久久久a久久爽久久v久久| 亚洲国产av新网站| 99热全是精品| 大又大粗又爽又黄少妇毛片口| 成人亚洲欧美一区二区av| 欧美性猛交╳xxx乱大交人| 免费黄色在线免费观看| 欧美高清成人免费视频www| 免费电影在线观看免费观看| 久久久精品欧美日韩精品| 亚洲丝袜综合中文字幕| 亚洲怡红院男人天堂| 成人二区视频| av在线观看视频网站免费| 亚洲精品日本国产第一区| 亚洲国产精品成人综合色| 如何舔出高潮| 99视频精品全部免费 在线| 亚洲精品久久午夜乱码| 国产成人a区在线观看| 日本色播在线视频| 亚洲va在线va天堂va国产| 国产淫语在线视频| 日韩制服骚丝袜av| 韩国高清视频一区二区三区| 欧美激情久久久久久爽电影| 激情 狠狠 欧美| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 午夜激情福利司机影院| 边亲边吃奶的免费视频| 王馨瑶露胸无遮挡在线观看| 精品酒店卫生间| 丰满人妻一区二区三区视频av| 97在线人人人人妻| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 精品国产露脸久久av麻豆| 超碰av人人做人人爽久久| 99久久精品一区二区三区| 3wmmmm亚洲av在线观看| 好男人在线观看高清免费视频| av福利片在线观看| 国产高清有码在线观看视频| 岛国毛片在线播放| 欧美激情在线99| 18+在线观看网站| 亚洲精品色激情综合| 在线观看一区二区三区激情| 18禁裸乳无遮挡动漫免费视频 | 日本熟妇午夜| 夫妻午夜视频| 美女高潮的动态| 天堂中文最新版在线下载 | 狂野欧美白嫩少妇大欣赏| 交换朋友夫妻互换小说| 免费观看的影片在线观看| 国产老妇女一区| 国产成人免费无遮挡视频| 日本一本二区三区精品| 少妇被粗大猛烈的视频| 最近2019中文字幕mv第一页| 免费观看在线日韩| 久久午夜福利片| 国产探花极品一区二区| 免费观看无遮挡的男女| 久久久久久久久久成人| 在线观看美女被高潮喷水网站| 国产免费又黄又爽又色| 中国三级夫妇交换| 精品99又大又爽又粗少妇毛片| 婷婷色麻豆天堂久久| 亚洲成人久久爱视频| 日本黄大片高清| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 久久久国产一区二区| 丝袜喷水一区| 国产成人91sexporn| 亚洲精品第二区| 日本免费在线观看一区| 午夜福利视频精品| 午夜免费观看性视频| 街头女战士在线观看网站| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩东京热| 国产黄片视频在线免费观看| 日韩一区二区视频免费看| 国产精品秋霞免费鲁丝片| 国产精品久久久久久久久免| 少妇丰满av| 亚洲av不卡在线观看| 免费av毛片视频| 午夜精品国产一区二区电影 | 国产成年人精品一区二区| 欧美高清性xxxxhd video| 1000部很黄的大片| 大香蕉97超碰在线| 色网站视频免费| 日韩欧美精品免费久久| 欧美人与善性xxx| 国产综合精华液| 不卡视频在线观看欧美| 亚洲婷婷狠狠爱综合网| 国产探花极品一区二区| 久久6这里有精品| 日本熟妇午夜| 成人漫画全彩无遮挡| 熟女av电影| 我要看日韩黄色一级片| 人妻一区二区av| 久久人人爽av亚洲精品天堂 | 嘟嘟电影网在线观看| 91午夜精品亚洲一区二区三区| 在现免费观看毛片| 国产男女超爽视频在线观看| videossex国产| 国产精品一区二区在线观看99| 丝袜美腿在线中文| 九草在线视频观看| 免费观看性生交大片5| 亚洲,欧美,日韩| 国产av国产精品国产| 国产黄片美女视频| 国模一区二区三区四区视频| 熟女电影av网| 成年人午夜在线观看视频| av在线蜜桃| 狂野欧美激情性bbbbbb| 成年女人看的毛片在线观看| 精品久久久精品久久久| 亚洲av日韩在线播放| 一本一本综合久久| 精品熟女少妇av免费看| 大香蕉97超碰在线| 成人毛片a级毛片在线播放| 亚洲欧美中文字幕日韩二区| 国产乱人视频| 亚洲精品成人av观看孕妇| 免费高清在线观看视频在线观看| 国产人妻一区二区三区在| 两个人的视频大全免费| 国产精品久久久久久精品电影小说 | 欧美xxxx性猛交bbbb| 日本与韩国留学比较| 色视频www国产| 日韩精品有码人妻一区| tube8黄色片| 久久久久久久午夜电影| 国产伦精品一区二区三区四那| 亚洲国产欧美人成| 人妻 亚洲 视频| 激情五月婷婷亚洲| kizo精华| 欧美一区二区亚洲| 欧美+日韩+精品| 欧美性猛交╳xxx乱大交人| 国产成人一区二区在线| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 免费电影在线观看免费观看| 欧美一级a爱片免费观看看| 我的女老师完整版在线观看| 久久精品国产a三级三级三级| 亚洲精品成人久久久久久| 如何舔出高潮| 久久国产乱子免费精品| 国产精品伦人一区二区| 男女下面进入的视频免费午夜| 亚洲国产精品国产精品| 草草在线视频免费看| 国产日韩欧美亚洲二区| 美女主播在线视频| 免费黄网站久久成人精品| 国产男人的电影天堂91| 黄色欧美视频在线观看| 国产精品精品国产色婷婷|