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

    念青唐古拉山東南坡植被群落數(shù)量生態(tài)分析及群落多樣性

    2015-12-08 02:43:19曲廣鵬參木友趙景學(xué)陳曉鵬張濤白彥福
    生態(tài)環(huán)境學(xué)報 2015年10期
    關(guān)鍵詞:唐古拉山海拔高度草甸

    曲廣鵬,參木友,趙景學(xué),陳曉鵬,張濤,白彥福

    1. 西藏自治區(qū)農(nóng)牧科學(xué)院草業(yè)科學(xué)研究所,西藏 拉薩 850000;2. 中國科學(xué)院青藏高原研究所,北京 100101;3. 中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,四川 成都 610041;4. 蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020

    念青唐古拉山東南坡植被群落數(shù)量生態(tài)分析及群落多樣性

    曲廣鵬1*,參木友1,趙景學(xué)2,陳曉鵬3,張濤4,白彥福4

    1. 西藏自治區(qū)農(nóng)牧科學(xué)院草業(yè)科學(xué)研究所,西藏 拉薩 850000;2. 中國科學(xué)院青藏高原研究所,北京 100101;3. 中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,四川 成都 610041;4. 蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020

    研究沿海拔梯度對植被群落的影響,有助于人們進(jìn)一步理解未來氣候變化背景下,高寒生態(tài)系統(tǒng)結(jié)構(gòu)和功能的響應(yīng)模式。選擇念青唐古拉山東南坡的典型生境,包括從海拔4 775 m到5 305 m全部范圍內(nèi)的植物群落。沿著海拔梯度每隔約30 m設(shè)置1個樣帶,并在各個樣帶上隨機設(shè)置樣方,調(diào)查樣方內(nèi)物種數(shù)、物種高度、物種蓋度、物種頻度。采用群落多樣性指數(shù)、應(yīng)用雙向指示種分析法(TWINSPAN)和除趨勢對應(yīng)分析(DCA)等方法,對西藏念青唐古拉山東南坡高山草甸植物群落進(jìn)行了分析。TWINSPAN將所有植物群落劃分為5種類型,分別為:香柏(Sabina pingii var.wilsonii)+高山嵩草(Kobresia pygmaea)灌叢草甸群落;高山嵩草+圓穗蓼(Polygonum macrophyllum)草甸群落;高山嵩草+矮生嵩草(Kobresia humilis)草甸群落;高山嵩草草甸群落;流石灘冰緣植被群落。這些群落分別位于不同的海拔高度。分類結(jié)果很好地反映了植物群落類型分布與地形、海拔的關(guān)系,并在DCA二維排序圖上得到了較好的驗證。念青唐古拉山東南坡草地植物群落多樣性在高海拔、低海拔地區(qū)較低,中間海拔高度地區(qū)較高。植物群落多樣性指數(shù)呈現(xiàn)出草本層>灌叢層>流石灘冰緣的特征。多樣性指數(shù)與海拔高度之間的趨勢模擬均呈負(fù)二次函數(shù)關(guān)系,單峰式函數(shù)關(guān)系能較好地表達(dá)不同海拔梯度植物群落多樣性和均勻度與海拔之間的分布格局。

    念青唐古拉山;植被類型;數(shù)量生態(tài)學(xué);生物多樣性;海拔梯度;雙向指示種分析法;除趨勢對應(yīng)分析

    植物與環(huán)境關(guān)系是植被生態(tài)學(xué)研究的重點之一,而分類和排序是研究植物與環(huán)境關(guān)系的兩種重要數(shù)量方法,它為客觀、準(zhǔn)確地揭示植被、植物群落及植物與環(huán)境之間的生態(tài)關(guān)系提供了合理、有效的途徑,是國際上植被生態(tài)學(xué)重要的研究內(nèi)容之一(張金屯,1991),已被廣泛應(yīng)用于國內(nèi)的植被生態(tài)學(xué)研究(張峰等,2000;賈曉妮等,2007;張新時,1991;彭少麟,1988;尚占環(huán)等,2005)。生物多樣性作為全球變化的重要研究內(nèi)容之一,是當(dāng)前國際生態(tài)學(xué)研究的熱點問題。目前,在生物多樣性與生態(tài)系統(tǒng)功能研究中,草地生態(tài)系統(tǒng)備受關(guān)注(楊元合等,2004)。全球溫度升高已經(jīng)導(dǎo)致了高緯度及高海拔地區(qū)陸地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能發(fā)生了顯著變化。在國際生物圈(IBP)研究計劃中,青藏高原被列為全球氣候變化的敏感區(qū)域,這種極端環(huán)境下發(fā)育的植被和土壤對氣候變化極為敏感,是研究生態(tài)系統(tǒng)對氣候變化響應(yīng)與適應(yīng)機制的天然實驗室(王建林等,2009)。

    念青唐古拉山脈位于我國西藏自治區(qū),是青藏高原主要山脈之一,屬斷塊山。西接岡底斯山脈,東南延伸與橫斷山脈伯舒拉嶺相接,中部略為向北凸出,同時將西藏劃分為藏北、藏南、藏東南三大區(qū)域。在海拔較低的暖熱地區(qū)尚有常綠闊葉林及鐵杉林分布,森林帶以上則為高山灌叢草甸及高山草甸帶,面積較廣,是當(dāng)?shù)刂饕烊幌募灸翀觯m宜放養(yǎng)牦牛、綿羊等牲畜。近年來,許多學(xué)者就念青唐古拉山陸-氣系統(tǒng)水熱特征、氣候特征、高寒草原生態(tài)系統(tǒng)表層土壤有機碳分布特征等方面的研究較多(謝健等,2009;王建林等,2009),而關(guān)于念青唐古拉山東南坡植被數(shù)量分析及群落多樣性

    的研究鮮有報道。研究沿海拔梯度的不同水熱組合對植被群落的影響,有助于我們進(jìn)一步理解未來氣候變化背景下,高寒生態(tài)系統(tǒng)結(jié)構(gòu)和功能的響應(yīng)模式。基于此,選擇念青唐古拉山東南坡植被進(jìn)行研究,分析沿海拔梯度植被群落分布特征,為進(jìn)一步了解念青唐古拉山植物群落發(fā)展趨勢提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    念青唐古拉山脈橫跨青藏高原腹地,總體上呈現(xiàn)東北走向,主峰海拔為 7162 m,山頂海拔一般為5600~6200 m;西北側(cè)為納木錯盆地,東南側(cè)為當(dāng)雄羊八井地塹,構(gòu)成青藏高原內(nèi)部的典型盆-山構(gòu)造地貌,山脈寬度在東北部約為10 km,西南部約35 km。念青唐古拉山脈的山頂面與兩側(cè)盆地相對高差大約800~1600 m,屬青藏高原內(nèi)部的高山-極高山地區(qū)。本研究選擇位于藏北腹地的當(dāng)雄縣城(30°32′N,91°03′E)西北部的念青唐古拉山東南坡的高山植被群落為研究對象,海拔4775~5305 m,屬于高原亞寒帶季風(fēng)半干旱氣候。據(jù)當(dāng)雄縣氣象站記載,近40余年(1963─2010年)以來,該地區(qū)氣溫升高明顯,降水年紀(jì)間差異大。該區(qū)多年(1963─2010年)年均氣溫1.8 ℃,最冷月(1月)均溫-10.4 ℃,最熱月(7月)均溫10.7 ℃。年均降水量479 mm,其中90%集中在5─9月(王忠,2011)。年均濕潤度K值為0.7,植被為以寒旱生的高山嵩草為優(yōu)勢種的高山草甸(王建林等,2009)。本地區(qū)土壤為高山草甸土,土壤厚度在0.3~1.0 m之間(付剛,2013)。

    1.2 植被調(diào)查與取樣

    選擇念青唐古拉山東南坡能代表各個生境的典型地區(qū),包括從海拔4775 m到5305 m全部范圍內(nèi)的植物群落,每隔約30 m海拔設(shè)置1條樣帶,共18條。在各個樣帶上隨機設(shè)置5~8個樣方,草本植被樣方為0.25 m×0.25 m,灌叢樣方為10 m×20 m,主要調(diào)查樣方內(nèi)物種數(shù)、物種高度、物種蓋度、物種頻度,共108個樣方。蓋度采用針刺法;頻度為在每個樣地內(nèi)采用頻度圈隨機 20次,計數(shù)圈內(nèi)物種的有無。同時記錄樣地海拔高度、坡向、坡度等地形特征。植被調(diào)查工作全部集中在 2010年 9月完成。

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

    數(shù)據(jù)分析根據(jù)張金屯(2004)論述,選擇雙向指示種分類法(TWINSPAN)和除趨勢對應(yīng)分析排序法(DCA)。采用植被數(shù)量分析程序WinTWINS和Canoco分別完成植物群落的TWINSPAN分類和DCA排序分析,計算時以植物重要值為屬性進(jìn)行分類和排序計算。進(jìn)行排序和分類前先將調(diào)查的草本樣方84個、灌叢樣方24個以植物物種重要值大小進(jìn)行初步分類,將樣方合并,整理出 18個樣地植物數(shù)據(jù),對 18個樣地植物群落進(jìn)行分類、排序分析,研究18個樣地植被的關(guān)系和分布特征。

    物種多樣性指數(shù)(Diversity index)、豐富度指數(shù)(Richness index)和均勻度指數(shù)(Evenness index)的計算公式如下(陳亞明等,2004;李新榮等,2000):

    群落多樣性指數(shù)的運算采用 Shannon-Wiener指數(shù):

    Simpson指數(shù):

    豐富度指數(shù)的運算采用Margalef指數(shù):

    重要值=(相對高度+相對蓋度+相對頻度)/3(5)

    式中,S為每個樣地的總物種數(shù),N為每個樣方中的全部物種的總個體數(shù),Pi為種i的重要值。

    圖1 念青唐古拉山東南坡植物群落TWNSPAN分類結(jié)果樹狀示意圖Fig. 1 Dendrogram of florae in the southeast slope of Mountain Nyenquentanglha by TWINSPAN

    均勻度指數(shù)的運算采用Pielou指數(shù):

    2 結(jié)果與分析

    2.1 TWINSPAN分類結(jié)果與群落類型概況

    雙向指示種分類法(TWINSPAN)的第1次劃分將所有樣地分為2類(圖1),一類是流石灘冰

    緣植被群落,包括4個樣地;另一類為高寒草甸和高山灌叢植被群落包括 14個樣地,逐步將樣方繼續(xù)分為6組,結(jié)合山地地形、植被分布的土壤生境的生態(tài)特征,在TWINSPAN基礎(chǔ)上將該山地植被分為5組。各組物種優(yōu)勢度如表1。

    表1 念青唐古拉山東南坡主要植被種類重要值Table 1 Species importance value of major vegetation in the southeast slope of Mountain Nyenquentanglha

    I香柏(Sabina pingii var.wilsonii)+高山嵩草(Kobresia pygmaea)灌叢草甸群落。該群落分布于海拔4775~4889 m,包含樣地1~5。優(yōu)勢植物有香柏和高山嵩草,優(yōu)勢度分別為0.43和0.15。伴生植物主要有青藏苔草(Moocroftii falc.ex boott)、矮生嵩草(Kobresia humilis)、釘柱委陵菜(Saundersiana royle)冰川棘豆(Oxytropis glacialis)、木根香青(Anaphalis xylorhiza)和粉花雪靈芝(Arenaria shannanensis)等。

    II 高山嵩草+圓穗蓼(Polygonum macrophyllum)草甸群落。該群落分布于海拔4996~5050 m,包含樣地8、9和10。優(yōu)勢種為高山嵩草和圓穗蓼,優(yōu)勢度分別為 0.33和 0.20。伴生植物主要有矮生嵩草、藏西風(fēng)毛菊(Saussurea stoliczkai)、獨一味(Lamiophlomis rotata)、釘柱委陵菜、冰川棘豆等。

    III高山嵩草+矮生嵩草草甸群落。該群落分布于海拔4945~5163 m,包含樣地6、7和樣地11、12。優(yōu)勢植被為高山嵩草和矮生嵩草,優(yōu)勢度分別為0.26和0.12。伴生植物主要有墊狀點地梅、肉果草(Lancea tibetica)、釘柱委陵菜、圓穗蓼等。

    IV高山嵩草草甸群落。該群落分布于海拔5163~5196 m,包含樣地13和14。高山嵩草為該群落的單優(yōu)勢種,優(yōu)勢度為0.25。伴生植被有矮生嵩草、高山早熟禾(Poa alpina Linn.)、墊狀金露梅(Potentilla fruticosa var. pumila)、紅紋馬先蒿(Pedicularis striata)等。

    V流石灘冰緣植被群落。該群落分布于海拔5238~5305 m,包含樣地15~18。高山早熟禾和墊狀棱點地梅(Tapete maxim.)為該群落的優(yōu)勢種,優(yōu)勢度分別為0.23和0.22。伴生植被有冷蒿、雪蓮、青藏苔草、木根香青等。

    2.2 草甸群落的DCA排序

    應(yīng)用除趨勢對應(yīng)分析排序法(DCA)對念青唐古拉山東南坡高山草甸 18個樣地進(jìn)行分析,為了

    便于和 TWINSPAN分類結(jié)果進(jìn)行比較,將TWINSPAN的最終分組結(jié)果在DCA二維軸上標(biāo)出(Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ)(圖2)。從該圖我們也可以看出TWINSPAN分類結(jié)果和DCA排序結(jié)果較吻合。DCA第一軸表現(xiàn)為海拔高度和土壤質(zhì)地,第二軸表現(xiàn)為坡度和土壤水分。按第一軸群落順序,從左到右,隨著海拔增加,土質(zhì)逐漸變差,植被群落逐漸從Ⅰ轉(zhuǎn)變?yōu)棰?。DCA第二軸群落順序,從下到上,隨著坡度逐漸增加,土壤水分降低,植被群落逐漸從Ⅰ轉(zhuǎn)變?yōu)棰?。從上述分析可見,DCA第一軸和第二軸有密切關(guān)系,草甸群落的分布明顯受海拔、坡度、土壤水分、土壤質(zhì)地等環(huán)境因子的制約,說明DCA排序較好地反映了念青唐古拉山東南坡高山草甸群落的生態(tài)關(guān)系。

    圖2 研究區(qū)內(nèi)植物群落的DCA二維排序Fig. 2 Two-dimensional ordiration of the florae in the region by the DCA

    2.3 植被群落多樣性分析

    念青唐古拉山東南坡不同植物群落類型的多樣性程度是綜合生態(tài)因素作用的結(jié)果,其中海拔最具重要性。海拔以及由于海拔造成的水分、熱量條件的差異是影響環(huán)境因子變化的主導(dǎo)因素,從而導(dǎo)致不同類型植物群落的物種多樣性也發(fā)生了很大的變化(表2)。Shannon-Wiener指數(shù),群落Ⅳ最高,群落Ⅱ、Ⅲ次之,群落Ⅰ較低,群落Ⅴ最低。其中,群落Ⅴ與其他群落差異均顯著(P<0.05);群落Ⅳ與群落Ⅰ和Ⅱ差異顯著(P<0.05),與群落Ⅲ差異不顯著(P>0.05);群落Ⅰ、Ⅱ、Ⅲ之間差異不顯著(P>0.05)。Simpson指數(shù),群落Ⅳ最高,群落Ⅰ、Ⅱ、Ⅲ次之,群落Ⅴ最低。其中,群落Ⅴ與其他群落差異均顯著(P<0.05);群落Ⅳ與群落Ⅰ差異顯著(P<0.05),與群落Ⅱ、Ⅲ差異不顯著(P>0.05);群落Ⅰ、Ⅱ、Ⅲ之間差異不顯著(P>0.05)。Margalef指數(shù),群落Ⅳ最高,群落Ⅰ、Ⅲ次之,群落Ⅴ較低,群落Ⅱ最低。其中,群落Ⅴ與群落Ⅲ之間差異顯著(P<0.05),與其他群落差異不顯著(P>0.05);群落Ⅳ與群落Ⅰ、Ⅱ、Ⅲ差異均顯著(P<0.05);群落Ⅰ與群落Ⅱ、Ⅲ差異不顯著(P>0.05);群落Ⅱ和群落Ⅲ之間差異顯著(P<0.05)。Pielou指數(shù),群落Ⅲ最高,群落Ⅲ次之,群落Ⅰ、Ⅳ較低,群落Ⅴ最低。其中,群落Ⅴ與群落Ⅰ、Ⅳ差異不顯著(P>0.05),與群落Ⅱ、Ⅲ差異顯著(P<0.05);群落Ⅳ與群落Ⅰ、Ⅱ、Ⅲ差異不顯著(P>0.05);群落Ⅱ與群落Ⅰ、Ⅲ差異顯著(P<0.05);群落Ⅰ、Ⅲ之間差異不顯著(P>0.05)。

    表2 念青唐古拉山東南坡植被群落多樣性指數(shù)Table 2 Community diversity index of the vegetation in southeast slope of Mountain Nyenquentanglha

    在每個海拔梯度,選擇每個樣帶對Shannon-Wiener指數(shù)、Simpson指數(shù)、Margalef指數(shù)、Pielou指數(shù)與海拔進(jìn)行回歸分析(圖 3)。結(jié)果表明:Shannon-Wiener指數(shù)、Simpson指數(shù)、Margalef指數(shù)、Pielou指數(shù)隨海拔高度變化,在一定程度上呈現(xiàn)出高海拔和低海拔地區(qū)較低,中海拔高度地區(qū)較高的趨勢。多樣性指數(shù)與海拔高度之間的趨勢模擬均呈負(fù)二次函數(shù)關(guān)系,除與Margalef指數(shù)擬合系數(shù)較小外,其他各多樣性指數(shù)與海拔高度的擬合系數(shù)均能較好的解釋二者之間的關(guān)系。說明單峰式函數(shù)關(guān)系能較好地表達(dá)不同海拔梯度植物群落多樣性和均勻度與海拔之間的分布格局。植被多樣性指數(shù)呈現(xiàn)出草本層>灌叢層>流石灘冰緣植被的特征(表2)。

    3 討論

    我國學(xué)者從20世紀(jì)70年代后期引入數(shù)量分類和排序技術(shù),開始了植被數(shù)量生態(tài)學(xué)的研究,但大多數(shù)研究集中在森林與草原植被中,而在高山-極高

    山地區(qū)植被研究中應(yīng)用很少。我們運用TWINSPAN和DCA方法對念青唐古拉山東南坡高山草甸進(jìn)行數(shù)量分類和排序,取得到了較好的結(jié)果。

    圖3 念青唐古拉山東南坡植物群落物種多樣性與海拔的關(guān)系Fig. 3 Relationship between community diversity and altitude of the southeast slope of Mountain Nyenquentanglha

    TWINSPAN等級分類將念青唐古拉山東南坡高山草甸18個樣地29個種劃分為5個植物群落類型,其分類矩陣圖反映出明顯的綜合環(huán)境梯度,即海拔高度、土壤水分、溫度梯度,比較客觀地反映了植物群落類型及優(yōu)勢種的分布與環(huán)境因子之間的關(guān)系,并在DCA二維排序圖上得到了較好的驗證。DCA排序表明群落生境所在地的海拔高度和水熱條件是決定群落類型分布的主要因素(Wang et al.,2013;Li et al.,2013)。第一軸反映了群落海拔高度的梯度變化,第二軸反映了土壤水熱組合的變化。一般來說,隨著山地迎風(fēng)坡海拔高度逐漸升高,溫度逐漸降低,而降水逐漸增加(王忠,2011;李瑞成,2012)。因此,高海拔地區(qū)土壤濕度相對于低海拔地區(qū)要高。然而,隨著海拔和濕度等的變化,念青唐古拉山東南坡高山草甸植物種類和植被高度發(fā)生明顯變化。在海拔較低的山坡上,香柏等灌木為優(yōu)勢種。隨著海拔增加,植被高度明顯下降,高山嵩草生長茂盛,逐漸成為優(yōu)勢種。尤其在海拔較高的流石灘山坡上,墊狀點地梅等一些墊狀植被優(yōu)勢度顯著增加,逐漸成為優(yōu)勢種。這是高寒植被對惡劣環(huán)境的一種適應(yīng)機制,具有對惡劣環(huán)境的指示作用,也是高寒墊狀植被生態(tài)系統(tǒng)工程師效應(yīng)的體現(xiàn)(何永濤等,2010;李瑞成,2012)。同時,高低海拔優(yōu)勢群落的分布不同,也指示了高低海拔可能存在迥然不同的環(huán)境限制因子。低海拔較低的降水,更適宜根系粗大的灌木生長。然而,低溫也是墊狀植被普遍分布在于高海拔地區(qū)的原因。

    念青唐古拉山坡降水量海拔梯度效應(yīng)明顯,隨海拔梯度的升高而增加,高海拔降水過程和低海拔差異顯著,最大降水高度位于山坡上部的 5100 m附近(謝健等,2009;王忠,2011;李瑞成,2012)。然而,念青唐古拉山坡草地群落,在降水量隨海拔的升高而增加的同時,土壤中砂粒含量則逐漸增加,相應(yīng)的粉粒和粘粒含量逐漸下降(Ohtsuka et al.,2008)。海拔以及由于海拔造成的水分、熱量條件的差異是影響環(huán)境因子變化的主導(dǎo)因素,導(dǎo)致不同類型植物群落的物種多樣性也發(fā)生了很大的變化(王忠,2011;李瑞成,2012)。關(guān)于高寒草甸群落植物多樣性沿海拔梯度變化的研究表明,單峰式函數(shù)能較好地表達(dá)不同海拔梯度植物群落物種豐富度、多樣性和均勻度與海拔間的分布格局,中海拔區(qū)植物群落多樣性最大(王長庭等,2004;趙新全,2009)。Wang et al.(2013)研究發(fā)現(xiàn),生產(chǎn)力在海拔梯度上的單峰分布是由于在不同海拔收到的限制因子不同而造成的。在念青唐古拉山

    東南坡,水分導(dǎo)致的干旱脅迫是低海拔的主要限制因子,低溫導(dǎo)致的水分和養(yǎng)份脅迫則是高海拔的主要限制因子。當(dāng)研究區(qū)域處在干旱半干旱地區(qū),海拔高度達(dá)到植物群落分布的上限時,上述因子的轉(zhuǎn)變造成青藏高原地區(qū)高山植被群落分布的單峰分布格局應(yīng)該是普遍存在的(王忠,2011;李瑞成,2012)。此外,海拔高度的變化也會改變了放牧強度,從而也會影響群落多樣性的海拔分布(候扶江等,2006)。然而,隨著海拔增加,放牧強度逐漸減弱(Hirota et al.,2009;Moriyama et al.,2013)。這也是中高海拔群落多樣性相對較高的可能原因。本文研究發(fā)現(xiàn),念青唐古拉山東南坡高山植被群落Shannon-Wiener指數(shù)、Simpson指數(shù)、Margalef指數(shù)、Pielou指數(shù)隨海拔高度變化,在一定程度上也呈現(xiàn)出高海拔和低海拔地區(qū)較低,中海拔高度地區(qū)較高的趨勢。多樣性指數(shù)與海拔高度之間的趨勢模擬均呈負(fù)二次函數(shù)關(guān)系,除與Margalef指數(shù)擬合系數(shù)較小外,其他各多樣性指數(shù)與海拔高度的擬合系數(shù)均能較好的解釋兩者之間的關(guān)系。說明單峰式函數(shù)關(guān)系能較好地表達(dá)不同海拔梯度植物群落多樣性和均勻度與海拔之間的分布格局。

    4 結(jié)果

    我們采用群落多樣性指數(shù)、應(yīng)用雙向指示種分析法(TWINSPAN)和除趨勢對應(yīng)分析(DCA)等方法,對西藏念青唐古拉山東南坡高山草甸植物群落進(jìn)行了分析。(1)TWINSPAN將所有植物群落劃分為 5種類型,分別為:香柏+高山嵩草灌叢草甸群落、高山嵩草+圓穗蓼草甸群落、高山嵩草+圓穗蓼草甸群落、高山嵩草+矮生嵩草草甸群落、高山嵩草草甸群落和流石灘冰緣植被群落。這些群落分別位于不同的海拔高度。分類結(jié)果很好地反映了植物群落類型分布與地形、海拔的關(guān)系,并在DCA二維排序圖上得到了較好的驗證。(2)念青唐古拉山東南坡草地植物群落多樣性在高海拔、低海拔地區(qū)較低,中間海拔高度地區(qū)較高。植物群落多樣性指數(shù)呈現(xiàn)出草本層>灌叢層>流石灘冰緣的特征。(3)多樣性指數(shù)與海拔高度之間的趨勢模擬均呈負(fù)二次函數(shù)關(guān)系,單峰式函數(shù)關(guān)系能較好地表達(dá)不同海拔梯度植物群落多樣性和均勻度與海拔之間的分布格局。

    HIROTA M, ZHANG P, GU S, et al. 2009. Altitudinal variation of ecosystem CO2fluxes in an alpine grassland from 3600 to 4200 m. Journal of Plant Ecology, 2(4): 197-205.

    LI R, LUO T, TANG Y, et al. 2013. The altitudinal distribution center of a widespread cushion species is related to an aptimum combination of temperature and precipitation in the central Tibetan Plateau [J]. Journal of Arid Environments, 88(1): 70-77.

    MORIYAMA A, YONEMURA S, KAWASHIMA S, et al. 2013. Environmental indicators for estimating the potential soil respiration rate in alpine zone [J]. Ecological Indicators, 32(9): 245-252.

    OHTSUKAT, HIROTA M, ZHANG X, et al. 2008. Soil organi carbon pools in alpine to nival zones along an altitudinal gradient (4400~5300 m) on the Tibetan Plateau [J]. Polar Sciene, 2(4): 277-285.

    WANG Z , LUO T, LI R, et al. 2013. Causes for the unimodal pattern of biomass and productivity in alpine grasslands along a large altitudinal gradient in semi-arid regions [J]. Journal of Vegetation Science, 24(1): 189-201.

    陳亞明, 李自珍, 杜國禎. 2004. 施肥對高寒草甸植物多樣性和經(jīng)濟類群的影響[J]. 西北植物學(xué)報, 24(3): 424-429.

    付剛. 2013. 不同海拔高寒草甸生態(tài)系統(tǒng)呼吸對實驗增溫和刈割的響應(yīng)[D]. 北京: 中國科學(xué)院研究生院: 10-20.

    何永濤, 石培禮, 閆巍. 2010. 高山墊狀植物的生態(tài)系統(tǒng)工程師效應(yīng)研究進(jìn)展[J]. 生態(tài)學(xué)雜志, 29(6): 1221-1227.

    候扶江, 楊中藝. 2006. 放牧對草地的作用[J]. 生態(tài)學(xué)報, 26(1): 244-264.

    賈曉妮, 程積民, 萬惠娥. 2007. DCA、CCA和DCCA三種排序方法在中國草地植被群落中的應(yīng)用現(xiàn)狀[J]. 中國農(nóng)學(xué)通報, 23(12): 391-395.

    李瑞成. 2012. 念青唐古拉山墊狀點地梅沿海拔梯度的分布、物候和生存率格局及其成因[D]. 北京: 中國科學(xué)院研究生院: 14-20.

    李新榮, 張景光, 劉立超, 等. 2000. 我國干旱沙漠地區(qū)人工植被與環(huán)境演變過程中植物多樣性的研究[J]. 植物生態(tài)學(xué)報, 24(3): 257-261.

    彭少麟. 1988. 廣東亞熱帶部分森林群落排序分析[J]. 武漢植物學(xué)研究, 6(1): 37-44.

    尚占環(huán), 姚愛興, 龍瑞軍, 等. 2005. 中衛(wèi)山羊核心產(chǎn)地植物群落的數(shù)量分類與排序[J]. 西北植物學(xué)報, 25(5): 985-990.

    王建林, 歐陽華, 王忠紅, 等. 2009. 念青唐古拉山東南坡高寒草原生態(tài)系統(tǒng)表層土壤有機碳分布特征及影響因素[J]. 地理科學(xué), 29(3): 385-390.

    王建林, 歐陽華, 王忠紅, 等. 2010. 念青唐古拉山東南坡高寒草原生態(tài)系統(tǒng)表層土壤活性有機碳的影響因素研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 18(2): 235-240.

    王長庭, 王啟基, 龍瑞軍, 等. 2004. 高寒草甸群落植物多樣性和初級生產(chǎn)力沿海拔梯度變化的研究[J]. 植物生態(tài)學(xué)報, 28(2): 240-245.

    王忠. 2011. 念青唐古拉山南坡高寒草甸生產(chǎn)力沿海拔梯度的變化機理研究[D]. 北京: 中國科學(xué)院研究生院: 12-24.

    謝健, 劉景時, 杜明遠(yuǎn), 等. 2009. 念青唐古拉山南坡夏季降水梯度觀測[J]. 水資源與水工程學(xué)報, 20(5): 14-18.

    謝健, 劉景時, 杜明遠(yuǎn), 等. 2010. 念青唐古拉山南坡氣溫分布及其垂直梯度[J]. 地理科學(xué), 30(2): 113-118.

    謝健, 劉景時, 杜明遠(yuǎn), 等. 2010. 念青唐古拉山西段高海拔陸-氣系統(tǒng)水熱特征[J]. 地理科學(xué)進(jìn)展, 29(2): 151-158.

    楊元合, 饒勝, 胡會峰, 等. 2004. 青藏高原高寒草地植物物種豐富度及其與環(huán)境因子和生物量的關(guān)系[J]. 生物多性, 12(1): 200-205.

    張峰, 張金屯. 2000. 我國植被數(shù)量分類和排序研究進(jìn)展[J]. 山西大學(xué)學(xué)報(自然科學(xué)版), 23(3): 278-282.

    張金屯. 1991. 植被數(shù)量分類和排序的發(fā)展[J]. 山西大學(xué)學(xué)報(自然科學(xué)版), 14(2): 215-224.

    張金屯.2004. 數(shù)量生態(tài)學(xué)[M]. 北京: 科學(xué)出版社: 1-70.

    張新時. 1991. 西藏阿里植物群落的間接梯度分析、數(shù)量分類與環(huán)境解釋[J]. 植物生態(tài)學(xué)與地植物學(xué)學(xué)報, 15(2): 101-113.

    趙新全. 2004. 高寒草甸生態(tài)系統(tǒng)與全球變化[M]. 科學(xué)出版社: 106-111.

    Quantitative Ecology and Species Diversity of the Vegetation in Southeast Slope of the Nyenchentanglha Mountain

    QU Guangpeng1, CAN Muyou1, ZHAO Jingxue2, CHEN Xiaopeng3, ZHANG Tao4, BAI Yanfu4
    1. Grassland Institute of Tibet Academy of Agricultural and Animal Husbandry Sciences, Lahsa 850000 China; 2. Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100101 China; 3. Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China; 4. College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020 China.

    Study on the response of vegetation to altitude gradient is helpful for us to understand how alpine ecosystems will response to future climate change. We chose to study the vegetation along the south slope of the Nyenchentanglha Mountain, with the altitude gradient from 4 775 m to 5 305 m. Along the altitude gradient, samples were set every 30 meters and species, with plant height, coverage and species frequency investigated in each samples. Double gaps indicator species analysis (TWINSPAN), and Detrended correspondence analysis (DCA) were applied to analysis the community types and species diversity along the altitude gradient of the southeast slope of the Nyenchentanglha Mountain. TWINSPAN divided the plant communication into 5 types: Sabina pingii var. wilsonii + Kobresia pygmaea meadow community, Kobresia pygmaea + Polygonum macrophyllum meadow community, Kobresia pygmaea + Kobresia humilis meadow community, Kobresia pygmaea community and stone beach periglacial vegetation communities. The result of vegetation’s classification by TWINSPAN reflected the relation between community’s location, topography and altitude, which were verified in pattern map of DCA very well. The species was lower in high and low altitude area compared with the middle of slope of mountain, which formed a hump cure. The species diversity index from high to low of community was herb layer >shrub layer > periglacial stone-community. The model of diversity index and altitude was negatively quadratic function.

    Nyenchentanglha Mountain; vegetation types; Quantitative ecology; biodiversity; altitude variation; TWINSPAN, DCA

    10.16258/j.cnki.1674-5906.2015.10.005

    Q948;X176

    A

    1674-5906(2015)10-1618-07

    曲廣鵬,參木友,趙景學(xué),陳曉鵬,張濤,白彥福. 念青唐古拉山東南坡植被群落數(shù)量生態(tài)分析及群落多樣性[J].生態(tài)環(huán)境學(xué)報, 2015, 24(10): 1618-1624.

    QU Guangpeng, CAN Muyou, ZHAO Jingxue, CHEN Xiaopeng, ZHANG Tao, BAI Yanfu. Quantitative Ecology and Species Diversity of the Vegetation in Southeast Slope of the Nyenchentanglha Mountain [J]. Ecology and Environmental Sciences, 2015, 24(10): 1618-1624.

    “十二五”國家科技支撐項目;西藏牧區(qū)“生產(chǎn)生態(tài)生活”保障技術(shù)集成與示范(2012BAD13B04)

    曲廣鵬(1981年生),男,助理研究員,碩士,主要從事牧草栽培與草地生態(tài)工作。*通信作者。

    2014-08-28

    猜你喜歡
    唐古拉山海拔高度草甸
    高山草甸
    The world's narrowest river
    山地草甸
    故障狀態(tài)下純電動汽車環(huán)境壓力及海拔高度估算方法
    北京汽車(2021年1期)2021-03-04 13:05:46
    在納木錯湖邊
    草地(2020年3期)2020-07-04 02:07:42
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    三條河在我桌頭流淌
    大地賦
    滇池(2016年12期)2016-12-10 22:21:19
    不同海拔高度人群胃蛋白酶原血清水平差異性分析
    我渴望的雪還沒有落下來
    詩潮(2013年11期)2013-11-16 14:59:05
    热99re8久久精品国产| 国产精品一区二区在线不卡| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av| 精品熟女少妇八av免费久了| 亚洲精华国产精华精| 后天国语完整版免费观看| www国产在线视频色| 中文字幕色久视频| 日日夜夜操网爽| 欧美乱码精品一区二区三区| 免费久久久久久久精品成人欧美视频| 老司机福利观看| 夫妻午夜视频| 99精品在免费线老司机午夜| 欧美日韩亚洲高清精品| 亚洲欧美日韩无卡精品| 叶爱在线成人免费视频播放| 在线观看免费日韩欧美大片| 老司机午夜十八禁免费视频| 美女扒开内裤让男人捅视频| 麻豆一二三区av精品| 女人爽到高潮嗷嗷叫在线视频| 成人永久免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品爽爽va在线观看网站 | 免费久久久久久久精品成人欧美视频| 色播在线永久视频| 国产精品乱码一区二三区的特点 | 免费搜索国产男女视频| 欧美亚洲日本最大视频资源| 国产成+人综合+亚洲专区| 女性生殖器流出的白浆| 在线永久观看黄色视频| 美女午夜性视频免费| 一a级毛片在线观看| 在线观看一区二区三区| 亚洲成人精品中文字幕电影 | 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 精品久久久精品久久久| 日韩国内少妇激情av| 亚洲色图综合在线观看| av天堂久久9| 91精品国产国语对白视频| 欧美中文综合在线视频| 精品国内亚洲2022精品成人| 欧美乱色亚洲激情| 99国产极品粉嫩在线观看| 一区福利在线观看| 免费在线观看完整版高清| 成人亚洲精品一区在线观看| 久久久水蜜桃国产精品网| av在线天堂中文字幕 | 少妇粗大呻吟视频| 999久久久精品免费观看国产| 香蕉丝袜av| 水蜜桃什么品种好| tocl精华| 在线永久观看黄色视频| 欧美日韩瑟瑟在线播放| 精品熟女少妇八av免费久了| 久久久精品欧美日韩精品| 亚洲精品久久午夜乱码| 精品乱码久久久久久99久播| 在线观看www视频免费| 亚洲国产中文字幕在线视频| 午夜免费鲁丝| 国产99久久九九免费精品| 黄网站色视频无遮挡免费观看| 91字幕亚洲| 亚洲av成人不卡在线观看播放网| 男女床上黄色一级片免费看| 一夜夜www| 18禁裸乳无遮挡免费网站照片 | 欧美色视频一区免费| 欧美中文日本在线观看视频| 中文字幕av电影在线播放| 精品一区二区三区视频在线观看免费 | 久久99一区二区三区| 制服诱惑二区| 熟女少妇亚洲综合色aaa.| 亚洲男人的天堂狠狠| 欧美最黄视频在线播放免费 | 免费少妇av软件| 80岁老熟妇乱子伦牲交| 91成年电影在线观看| 国产色视频综合| 欧美丝袜亚洲另类 | 久久久久久大精品| 久久精品91无色码中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 无限看片的www在线观看| 黄网站色视频无遮挡免费观看| 男女做爰动态图高潮gif福利片 | 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| 两个人免费观看高清视频| 国产精品九九99| 一边摸一边抽搐一进一小说| 亚洲欧美精品综合久久99| 午夜老司机福利片| 韩国av一区二区三区四区| 超碰97精品在线观看| 岛国视频午夜一区免费看| 国产亚洲精品久久久久5区| 夜夜夜夜夜久久久久| 两个人看的免费小视频| 12—13女人毛片做爰片一| 国产av一区二区精品久久| 亚洲久久久国产精品| а√天堂www在线а√下载| 丝袜人妻中文字幕| 色综合欧美亚洲国产小说| 国产成+人综合+亚洲专区| 欧美一区二区精品小视频在线| 一本综合久久免费| 韩国精品一区二区三区| 亚洲男人天堂网一区| 十八禁人妻一区二区| 99精品久久久久人妻精品| 成熟少妇高潮喷水视频| 国产精品乱码一区二三区的特点 | 不卡一级毛片| 精品人妻在线不人妻| 丁香六月欧美| 老汉色∧v一级毛片| 一级a爱视频在线免费观看| 涩涩av久久男人的天堂| 天天添夜夜摸| 国产精品成人在线| 精品一区二区三区av网在线观看| 99国产精品一区二区三区| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 人人澡人人妻人| 日韩国内少妇激情av| 久久亚洲精品不卡| 久久亚洲真实| 久久青草综合色| 搡老熟女国产l中国老女人| 变态另类成人亚洲欧美熟女 | 男女午夜视频在线观看| 十八禁网站免费在线| 99香蕉大伊视频| 最新美女视频免费是黄的| 午夜两性在线视频| 99热只有精品国产| 国产成人影院久久av| 高清av免费在线| 99久久99久久久精品蜜桃| 精品国产亚洲在线| 妹子高潮喷水视频| 日韩三级视频一区二区三区| 午夜福利在线免费观看网站| 久久性视频一级片| 男女高潮啪啪啪动态图| 国产精品98久久久久久宅男小说| 亚洲精品av麻豆狂野| 亚洲中文字幕日韩| 丰满的人妻完整版| 国产色视频综合| 高清在线国产一区| 亚洲人成网站在线播放欧美日韩| 宅男免费午夜| 超色免费av| 丁香六月欧美| 男人舔女人的私密视频| 日本免费a在线| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 丁香欧美五月| 欧美日韩亚洲综合一区二区三区_| 真人一进一出gif抽搐免费| 中文字幕最新亚洲高清| 日本黄色视频三级网站网址| 极品教师在线免费播放| 成年女人毛片免费观看观看9| 欧美日韩精品网址| 黄色女人牲交| 免费看十八禁软件| 国产片内射在线| 夜夜躁狠狠躁天天躁| 久久国产亚洲av麻豆专区| 亚洲熟女毛片儿| 国产野战对白在线观看| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 久久精品亚洲精品国产色婷小说| 久久国产乱子伦精品免费另类| 18美女黄网站色大片免费观看| 成人亚洲精品av一区二区 | 国产精品二区激情视频| 91国产中文字幕| 国产成人av激情在线播放| 一级a爱视频在线免费观看| 亚洲精品成人av观看孕妇| 19禁男女啪啪无遮挡网站| 欧美日韩视频精品一区| 欧美日韩瑟瑟在线播放| 午夜a级毛片| 色综合站精品国产| 欧美丝袜亚洲另类 | 国产视频一区二区在线看| 免费看a级黄色片| 视频区欧美日本亚洲| bbb黄色大片| 精品国产乱子伦一区二区三区| 久久香蕉精品热| 亚洲精品一卡2卡三卡4卡5卡| 女生性感内裤真人,穿戴方法视频| 在线观看免费午夜福利视频| 亚洲男人的天堂狠狠| 亚洲精品成人av观看孕妇| 男女做爰动态图高潮gif福利片 | 在线观看舔阴道视频| 又紧又爽又黄一区二区| 天堂√8在线中文| 久久性视频一级片| 久久国产亚洲av麻豆专区| 日韩精品免费视频一区二区三区| 中亚洲国语对白在线视频| 黑人欧美特级aaaaaa片| 亚洲专区中文字幕在线| 亚洲 国产 在线| 亚洲精品在线观看二区| 久久久久精品国产欧美久久久| 欧美乱妇无乱码| 国产高清视频在线播放一区| 亚洲av五月六月丁香网| 欧美久久黑人一区二区| 超色免费av| 亚洲黑人精品在线| 久久亚洲真实| 在线观看午夜福利视频| 91国产中文字幕| 叶爱在线成人免费视频播放| 黄色丝袜av网址大全| 久久久国产成人精品二区 | 国产精品久久久久久人妻精品电影| 男女下面插进去视频免费观看| 99香蕉大伊视频| 成年人免费黄色播放视频| 亚洲国产精品一区二区三区在线| 男人操女人黄网站| 少妇的丰满在线观看| 欧美丝袜亚洲另类 | 极品教师在线免费播放| 99国产精品一区二区蜜桃av| 最近最新免费中文字幕在线| 男女之事视频高清在线观看| 激情视频va一区二区三区| 婷婷丁香在线五月| 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 亚洲一区二区三区不卡视频| 久久青草综合色| 久久人妻福利社区极品人妻图片| 日本免费a在线| 桃红色精品国产亚洲av| 这个男人来自地球电影免费观看| bbb黄色大片| 成人三级做爰电影| 国产麻豆69| 亚洲成国产人片在线观看| 亚洲欧美日韩高清在线视频| 好看av亚洲va欧美ⅴa在| 一区二区三区激情视频| 波多野结衣一区麻豆| 日日干狠狠操夜夜爽| 久久99一区二区三区| 中文字幕人妻熟女乱码| 中文字幕色久视频| 正在播放国产对白刺激| 免费一级毛片在线播放高清视频 | 精品人妻在线不人妻| 一二三四社区在线视频社区8| 亚洲狠狠婷婷综合久久图片| 一边摸一边抽搐一进一小说| 亚洲五月婷婷丁香| 亚洲色图 男人天堂 中文字幕| 桃色一区二区三区在线观看| 精品免费久久久久久久清纯| 国产精品久久久久久人妻精品电影| 极品人妻少妇av视频| 正在播放国产对白刺激| 久久精品国产综合久久久| 午夜福利在线观看吧| 老鸭窝网址在线观看| 免费在线观看亚洲国产| 午夜精品久久久久久毛片777| 色综合欧美亚洲国产小说| 丝袜美足系列| 免费在线观看影片大全网站| 国产成人精品久久二区二区免费| 亚洲国产欧美网| 高清av免费在线| 丝袜美足系列| 亚洲精品一卡2卡三卡4卡5卡| 一二三四在线观看免费中文在| 一进一出抽搐gif免费好疼 | av片东京热男人的天堂| 久久午夜综合久久蜜桃| 午夜成年电影在线免费观看| 免费在线观看日本一区| 欧美大码av| 欧美性长视频在线观看| 日韩国内少妇激情av| 久久久久久人人人人人| 在线视频色国产色| 久久久水蜜桃国产精品网| 99re在线观看精品视频| 亚洲精品粉嫩美女一区| 一二三四社区在线视频社区8| 国产av精品麻豆| 一级作爱视频免费观看| 欧美日韩亚洲国产一区二区在线观看| 极品人妻少妇av视频| 精品第一国产精品| 在线观看舔阴道视频| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 琪琪午夜伦伦电影理论片6080| 一边摸一边抽搐一进一小说| 亚洲一卡2卡3卡4卡5卡精品中文| 高清在线国产一区| 国产亚洲精品一区二区www| 变态另类成人亚洲欧美熟女 | 国产成年人精品一区二区 | 国产99白浆流出| 女人高潮潮喷娇喘18禁视频| 一二三四社区在线视频社区8| 日韩欧美三级三区| 91av网站免费观看| 在线观看日韩欧美| 欧美成狂野欧美在线观看| 欧洲精品卡2卡3卡4卡5卡区| 热99国产精品久久久久久7| 悠悠久久av| 亚洲专区中文字幕在线| 亚洲国产看品久久| 极品教师在线免费播放| 嫁个100分男人电影在线观看| 少妇裸体淫交视频免费看高清 | 久热这里只有精品99| 国产精品九九99| 久久精品成人免费网站| 中亚洲国语对白在线视频| 日本撒尿小便嘘嘘汇集6| svipshipincom国产片| 亚洲午夜理论影院| 三上悠亚av全集在线观看| 一级a爱视频在线免费观看| 欧美在线黄色| 久久人妻熟女aⅴ| 在线观看日韩欧美| 国产精品一区二区免费欧美| 国产片内射在线| 女人爽到高潮嗷嗷叫在线视频| 国产99白浆流出| 精品国产美女av久久久久小说| 国产成人精品无人区| 大香蕉久久成人网| 国产单亲对白刺激| www国产在线视频色| 亚洲人成伊人成综合网2020| 又黄又粗又硬又大视频| 久久人人精品亚洲av| 久久人人精品亚洲av| a级片在线免费高清观看视频| 欧美日韩乱码在线| 女人精品久久久久毛片| 亚洲午夜理论影院| 一级a爱视频在线免费观看| 久热这里只有精品99| 国产一区二区三区在线臀色熟女 | 操出白浆在线播放| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| ponron亚洲| 岛国视频午夜一区免费看| 男女午夜视频在线观看| 黄频高清免费视频| 久久九九热精品免费| 首页视频小说图片口味搜索| 久久青草综合色| 天堂俺去俺来也www色官网| 欧美午夜高清在线| 欧美黑人欧美精品刺激| 日韩欧美一区视频在线观看| 午夜福利,免费看| 亚洲九九香蕉| 91精品三级在线观看| 成人手机av| 一级毛片女人18水好多| 一边摸一边做爽爽视频免费| 国产免费av片在线观看野外av| 一个人免费在线观看的高清视频| 青草久久国产| 人妻丰满熟妇av一区二区三区| 久久久久久久久中文| 天堂俺去俺来也www色官网| 国产精品二区激情视频| 亚洲人成电影观看| 制服人妻中文乱码| 国产日韩一区二区三区精品不卡| 在线观看免费高清a一片| 中文字幕精品免费在线观看视频| 少妇粗大呻吟视频| 国产欧美日韩精品亚洲av| 超色免费av| 亚洲一区二区三区欧美精品| 热99re8久久精品国产| 国产精品电影一区二区三区| 老汉色∧v一级毛片| 老汉色av国产亚洲站长工具| 久久久国产欧美日韩av| 在线观看午夜福利视频| 性少妇av在线| 免费一级毛片在线播放高清视频 | 国产精品自产拍在线观看55亚洲| 女警被强在线播放| 欧美最黄视频在线播放免费 | 搡老乐熟女国产| 午夜福利影视在线免费观看| 女人高潮潮喷娇喘18禁视频| 亚洲第一青青草原| 欧美精品一区二区免费开放| 黄色a级毛片大全视频| 亚洲精品av麻豆狂野| 亚洲免费av在线视频| 咕卡用的链子| 欧美日韩亚洲高清精品| 日韩三级视频一区二区三区| 三级毛片av免费| 国产成人av激情在线播放| 久久国产精品男人的天堂亚洲| 日本免费a在线| 国产亚洲欧美精品永久| 亚洲av五月六月丁香网| 999久久久精品免费观看国产| 欧美黑人精品巨大| 久久久久久久午夜电影 | 日韩有码中文字幕| 老司机午夜十八禁免费视频| 欧美日韩一级在线毛片| 成年人黄色毛片网站| 欧美日韩乱码在线| 99精国产麻豆久久婷婷| 国产一区二区三区综合在线观看| 精品国产一区二区三区四区第35| 波多野结衣高清无吗| 一进一出好大好爽视频| 麻豆av在线久日| 午夜成年电影在线免费观看| 久久热在线av| 国产亚洲av高清不卡| 在线观看免费视频网站a站| 欧美成人性av电影在线观看| 国产成人欧美在线观看| 91在线观看av| 国产成人欧美| 欧美日韩中文字幕国产精品一区二区三区 | 色综合婷婷激情| 精品国产亚洲在线| 久久久久久大精品| 国产精品免费一区二区三区在线| 久久国产精品影院| 国产熟女xx| 免费少妇av软件| 国产黄色免费在线视频| 久久久久久久精品吃奶| 日本a在线网址| 亚洲男人天堂网一区| 成人黄色视频免费在线看| 久久人妻av系列| 亚洲欧美激情在线| 丰满饥渴人妻一区二区三| 国产一区二区三区综合在线观看| 欧美色视频一区免费| 18禁裸乳无遮挡免费网站照片 | 如日韩欧美国产精品一区二区三区| 欧美一级毛片孕妇| 两个人看的免费小视频| 久热爱精品视频在线9| 亚洲视频免费观看视频| 极品教师在线免费播放| 亚洲成国产人片在线观看| 亚洲av五月六月丁香网| 亚洲精品一二三| 亚洲五月婷婷丁香| 18美女黄网站色大片免费观看| 国产精华一区二区三区| 真人一进一出gif抽搐免费| 国产亚洲精品一区二区www| 少妇被粗大的猛进出69影院| av超薄肉色丝袜交足视频| 天堂影院成人在线观看| 国产精品99久久99久久久不卡| 国产亚洲av高清不卡| 国产精品国产av在线观看| 久久九九热精品免费| 国产精品一区二区免费欧美| 中文字幕精品免费在线观看视频| 欧美av亚洲av综合av国产av| 亚洲九九香蕉| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合久久99| 精品少妇一区二区三区视频日本电影| 亚洲激情在线av| 亚洲色图 男人天堂 中文字幕| 制服诱惑二区| 久久久久久久久免费视频了| 亚洲精品在线观看二区| 最近最新中文字幕大全免费视频| 久久午夜综合久久蜜桃| 日本黄色视频三级网站网址| 久久久国产精品麻豆| 日本三级黄在线观看| 韩国av一区二区三区四区| 国产亚洲精品第一综合不卡| 午夜亚洲福利在线播放| 久久久国产一区二区| 十分钟在线观看高清视频www| 国产精品野战在线观看 | 日韩一卡2卡3卡4卡2021年| 中出人妻视频一区二区| 天堂影院成人在线观看| 国产三级在线视频| 国产亚洲精品久久久久5区| 久久精品亚洲av国产电影网| 黄色丝袜av网址大全| 欧美+亚洲+日韩+国产| 国产三级黄色录像| 女人被狂操c到高潮| 欧美不卡视频在线免费观看 | 国产精品98久久久久久宅男小说| 亚洲七黄色美女视频| 脱女人内裤的视频| 免费少妇av软件| 国产精品偷伦视频观看了| 动漫黄色视频在线观看| 亚洲av熟女| 精品久久久久久久久久免费视频 | 纯流量卡能插随身wifi吗| 久久久国产成人免费| 亚洲av五月六月丁香网| 91精品三级在线观看| 人妻久久中文字幕网| 日本 av在线| 久久久精品国产亚洲av高清涩受| 久久精品亚洲熟妇少妇任你| 女性生殖器流出的白浆| 国产1区2区3区精品| 日本 av在线| 中文字幕人妻熟女乱码| 久久国产精品影院| 久久久久久久久久久久大奶| 国产激情欧美一区二区| 亚洲国产精品999在线| 国产精华一区二区三区| 欧美成人免费av一区二区三区| 国产熟女xx| 日韩国内少妇激情av| 天堂俺去俺来也www色官网| 欧美一区二区精品小视频在线| 一a级毛片在线观看| 国产精品一区二区三区四区久久 | 最近最新免费中文字幕在线| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 88av欧美| 日日摸夜夜添夜夜添小说| 国产成人精品在线电影| 国产激情欧美一区二区| 久久久久久久午夜电影 | 亚洲美女黄片视频| 中国美女看黄片| 国产精品免费视频内射| 国产一区二区三区在线臀色熟女 | 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| tocl精华| 一本大道久久a久久精品| 18禁国产床啪视频网站| 日韩人妻精品一区2区三区| 欧美激情久久久久久爽电影 | 99精品欧美一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 国产av一区在线观看免费| 在线观看一区二区三区| 男女高潮啪啪啪动态图| 99热只有精品国产| 在线观看免费高清a一片| 久久久国产精品麻豆| 一级毛片高清免费大全| 一级,二级,三级黄色视频| 成人特级黄色片久久久久久久| 极品人妻少妇av视频| 亚洲av成人不卡在线观看播放网| www.999成人在线观看| 国产精品亚洲av一区麻豆| 亚洲成人久久性| 国产精品偷伦视频观看了| 最好的美女福利视频网| 国产黄色免费在线视频| 又黄又爽又免费观看的视频| 亚洲久久久国产精品| 丁香六月欧美| 色精品久久人妻99蜜桃| 亚洲午夜精品一区,二区,三区| 很黄的视频免费| 男女午夜视频在线观看| 国产成人啪精品午夜网站|