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

    拉薩河流域虧組山植物物種豐富度和群落特征研究

    2021-06-10 09:01:47古桑群宗拉多郭應(yīng)杰向明學(xué)武俊喜潘影張燕杰李石勝達(dá)文彥
    廣西植物 2021年3期
    關(guān)鍵詞:環(huán)境因子群落結(jié)構(gòu)排序

    古桑群宗 拉多 郭應(yīng)杰 向明學(xué) 武俊喜 潘影 張燕杰 李石勝 達(dá)文彥

    摘 要:為研究草地植物群落、物種豐富度及其分布格局與影響因子之間的關(guān)系,該文以拉薩河流域林周縣卡孜鄉(xiāng)虧組山為研究地點(diǎn),對山體垂直樣帶 (3 900~5 100 m) 植物群落特征、植物物種豐富度與各影響因子之間的關(guān)系進(jìn)行了研究。研究區(qū)域共設(shè)置了13個(gè)樣帶 (每隔100 m設(shè)一個(gè)樣帶),每樣帶設(shè)置0.5 m × 0.5 m的5個(gè)樣方進(jìn)行植被調(diào)查,運(yùn)用主成分分析 (PCA) 和雙向指示種分析法 (TWINSPAN) 對植物群落進(jìn)行排序和分類,運(yùn)用冗余分析法 (RDA) 對群落及其分布格局與影響因子之間的關(guān)系進(jìn)行分析,對植物物種豐富度與各因子之間的關(guān)系進(jìn)行了回歸分析。結(jié)果表明:該山體植物可分為3個(gè)群落類型,含7個(gè)群叢;影響區(qū)域植物群落物種組成和分布格局的主要環(huán)境因子為海拔,其次是坡度;物種豐富度與海拔、分種蓋度呈單峰關(guān)系,與總蓋度、坡度、地上生物量呈正相關(guān)。該研究結(jié)果為區(qū)域植物物種多樣性和山地植物資源的保護(hù)和可持續(xù)利用提供了科學(xué)依據(jù)。

    關(guān)鍵詞: 物種豐富度, 環(huán)境因子, 分類, 排序, 群落結(jié)構(gòu), 虧組山

    中圖分類號:Q948

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1000-3142(2021)03-0372-12

    收稿日期:2019-09-15

    基金項(xiàng)目:國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目 (2016YFC0502004-5);西藏大學(xué)研究生高水平人才培養(yǎng)項(xiàng)目 (2017-GSP-111);西藏科技重大專項(xiàng) (XZ201801NA02)? [ Supported by the National Key Research and Development Program of China (2016YFC0502004-5);High-Level Talents Training Program in Tibet University (2017-GSP-111); Key Program of Science and Technology for Tibet, China (XZ201801NA02) ]。

    作者簡介: 古桑群宗(1995-),碩士研究生,研究方向?yàn)椴莸刂参锷鷳B(tài)學(xué),(E-mail)1992899926@qq.com。

    通信作者:拉多,博士,教授,碩士研究生導(dǎo)師,研究方向?yàn)橹参锷鷳B(tài)學(xué)、植物生理學(xué)、古生態(tài)學(xué),(E-mail)lhaduo@hotmail.com。

    Plant species richness and community characteristics? at Kuizu Mountain, Lhasa River basin

    GUSANG Qunzong1, LA Duo1*, GUO Yingjie1,4, XIANG Mingxue1,2, WU Junxi2,3, PAN Ying2,3, ZHANG Yanjie3, LI Shisheng1, DA Wenyan1

    ( 1. Department of Biology, Tibet University, Lhasa 850000, China; 2. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 3. Institute of Geographic Sciences and Resources, Chinese Academy of Sciences, Beijing 100101, China; 4. Yunnan Provincial Academy of Agricultural Sciences of Alpine Economic Plant Research Institute, Lijiang 674100, Yunnan, China )

    Abstract:In order to study the relationship between grassland plant community, species richness, its distribution pattern and their influencing variables. Taking the Kuizu Mountain in Lhasa Basin as the study area, the plant community structure, species richness along the altitudinal gradient and their relationships with influential variables were studied. A total of 13 sample plots were set up at every 100 m altitude with five samples of 0.5 m × 0.5 m at each plot in the study area for vegetation survey. The plant community structure and distribution were analyzed by using TWINSPAN (two-way indicator species analysis) and PCA (principal component analysis). RDA (redundancy analysis) was played to analyze the correlations between plant community structure, distribution and environmental variables (altitude, slope). The regression analysis was used to analyze the correlations between species richness and abiotic (altitude, slope) and biotic (total coverage, separate coverage and aboveground biomass) variables. PCA and TWINSPAN results showed that the vegetation in the study area was divided into three communities including seven plant associations. The results of RDA showed that altitude was the main variable influencing the plant community composition and its distribution pattern followed by slope in the study area. The results of regression analysis showed that the species richness had a single-peak relationship with altitude and separate coverage, and was positively correlated with total coverage, slope and aboveground biomass. This study provides a scientific baseline for the conservation of plant species diversity and sustainable utilization of mountain plant resources in the study area.

    Key words: species richness, environmental factor, TWINSPAN, PCA, community structure, Kuizu Mountain

    青藏高原具有獨(dú)特的高寒氣候和復(fù)雜的自然地理?xiàng)l件,是一個(gè)具有全球性重要意義的生物多樣性熱點(diǎn)地區(qū)和生態(tài)脆弱區(qū) (鐘祥浩等,2006)。植物物種多樣性的研究對區(qū)域生態(tài)系統(tǒng)的保護(hù)和其資源的永續(xù)利用以及環(huán)境評價(jià)等方面具有很重要的意義。影響物種多樣性的因素很多,在一個(gè)特定的地點(diǎn),海拔主導(dǎo)綜合性環(huán)境因子產(chǎn)生垂直分異 (宋長宏,2013),它關(guān)聯(lián)著溫度、降雨量、土壤、氣壓、輻射能等多種環(huán)境因子 (趙海珍,2010),海拔對物種分布有深刻的影響,相關(guān)的研究可以為預(yù)測全球氣候變化背景下的物種分布和生態(tài)系統(tǒng)的響應(yīng)提供科學(xué)依據(jù)?,F(xiàn)有的研究表明,高山草地植物豐富度與海拔之間存在三種不同的變化趨勢,即呈正相關(guān)、負(fù)相關(guān)和單峰分布 (Jií & Miroslav,2002;楊元合等,2004;李素清和張金屯,2007;朱源等,2007;段敏杰,2011;栗文瀚等,2017)。可見,在不同研究區(qū)域、海拔跨度,研究結(jié)果有一定的差異,并且在影響物種豐富度與群落結(jié)構(gòu)因子中,除了目前研究較為多的環(huán)境因子海拔以外,還有地形因子坡度等與物種豐富度和群落結(jié)構(gòu)有重要的影響 (宋永昌,2001;金樑等,2007;王娟,2013;劉冠成等,2018)。植物蓋度與生物量作為衡量植物群落保持水土功能 (郭忠升,2000) 和地面植物群落直觀定量化的功能指標(biāo),能夠反映區(qū)域生境質(zhì)量的好壞 (Harris et al., 2014)。因此,在生態(tài)系統(tǒng)的結(jié)構(gòu)和功能的形成上有重要的作用。在目前的研究中,關(guān)于生物量和植被蓋度與植物物種豐富度關(guān)系的研究發(fā)現(xiàn),物種豐富度與總蓋度呈正顯著相關(guān)關(guān)系,與分種蓋度(唐志堯和方精云,2004;拉瓊等,2014) 和生物量 (劉立斌等,2018) 呈單峰變化趨勢。

    區(qū)域草地植物物種分布格局及群落結(jié)構(gòu)的研究能夠?yàn)閰^(qū)域草地資源合理利用與保護(hù)提供科學(xué)依據(jù)。從已有的文獻(xiàn)來看,拉薩河流域山地植物物種豐富度及群落結(jié)構(gòu)在垂直梯度上的研究較多。但是,由于垂直尺度選擇不同而出現(xiàn)不同的研究結(jié)果,因此對植被因子 (植物蓋度、生物量等)、地形因子 (坡度) 與群落分布格局之間關(guān)系的研究卻相對較少。然而,在地形多樣化的青藏高原地區(qū),關(guān)于地形因子與植物多樣性的相關(guān)研究顯得尤為重要。本文以拉薩河流域林周縣卡孜鄉(xiāng)虧組山作為研究地點(diǎn),在海拔3 900~5 100 m的垂直梯度上進(jìn)行植被調(diào)查,對植物物種豐富度與植被因子、地形因子之間的相互關(guān)系進(jìn)行了初步研究,其結(jié)果對研究拉薩河流域大尺度山地植物物種豐富度與影響因子之間的關(guān)系具有一定的參考價(jià)值,為保護(hù)研究區(qū)域植物物種多樣性和山地植物資源可持續(xù)利用提供科學(xué)依據(jù)。

    1 研究區(qū)域與研究方法

    1.1 研究區(qū)域

    卡孜鄉(xiāng)虧組山 (91°07′33″—91°07′19″ E、29°52′23″—29°49′50″ N) 位于西藏拉薩河支流澎波藏曲流域西南部,隸屬于西藏自治區(qū)拉薩市林周縣,距林周縣10 km左右,其山體海拔跨度3 900~5 100 m之間,年日照時(shí)數(shù)大于3 000 h,年均溫7.5~7.6 ℃,年均降水量在440 mm左右,主要集中在6—9月 (張曉慶等,2017;金艷梅等,2018), 研究區(qū)域常見植物有高山嵩草 (Kobresia pygmaea)、釘柱委陵菜 (Potentilla saundersiana)、青藏苔草 (Carex moorcroftii)、肉果草 (Lancea tibetica)、草沙蠶 (Tripogon bromoides) 等。

    1.2 研究方法

    1.2.1 樣方調(diào)查 2018年7—9月,在海拔3 900~5 100 m的研究區(qū)域,每隔100 m設(shè)置一個(gè)樣地,共13個(gè)樣地,由于該區(qū)域的灌叢較少(只有2個(gè)樣方灌叢蓋度較大,分別為30%、46%,其余灌叢蓋度在0.5%~22%之間),并且需要長期監(jiān)測樣地,為了減少人為破壞對樣地的影響,每個(gè)樣地設(shè)置0.5 m × 0.5 m的5個(gè)樣方 (王長庭等,2005;鄭偉等,2010),共65個(gè)樣方,樣方序號與海拔帶相對應(yīng),如1~5號為3 900 m,6~10號為4 000 m,以此類推;樣方設(shè)置時(shí)盡可能避開特殊小生境 (如大塊石頭、人和牲畜干擾嚴(yán)重的地方)。記錄每個(gè)樣方內(nèi)的植物種類,并估算每個(gè)樣方植被的總蓋度與分種蓋度 (樣方內(nèi)各物種蓋度的總和)。用坡度儀記錄每個(gè)樣方的坡度,齊地面剪下植物,在80 ℃烘箱內(nèi)烘干24 h后稱干重,即為植物地上生物量,以下簡稱為生物量。

    1.2.2 數(shù)據(jù)分析 采用雙向指示種分析 (TWINSPAN-v3.2) 法對植物群落進(jìn)行分類;運(yùn)用CANOCO for Windows 4.5的主成分分析 (principal components analysis,PCA) 法對植物群落進(jìn)行排序分析;采用冗余分析 (redundancy analysis,RDA) 法對植物群落與環(huán)境因子之間的關(guān)系進(jìn)行了分析,同時(shí)采用SPSS19.0對植物物種豐富度與各因子間進(jìn)行了回歸分析。

    2 結(jié)果與分析

    2.1 虧組山植物群落TWINSPAN 分類

    TWINSPAN分類分析將虧組山65個(gè)樣方分成7組,按照“外貌-生態(tài)學(xué)分類方法” (中國科學(xué)院青藏高原綜合考察隊(duì),1980),并以依據(jù)群落生境特征的優(yōu)勢種和次優(yōu)勢中命名群叢類型,將其劃分為7個(gè)群叢 (圖1),具體如下:

    青藏苔草+高山嵩草+藏西嵩草 (Kobresia deasyi) 群叢 (a):包括18個(gè)樣方,主要分布海拔為4 800~5 100 m。在群叢水平上,以青藏苔草為優(yōu)勢種,次優(yōu)勢種有高山嵩草和藏西嵩草;主要的伴生種有釘柱委陵菜、圓穗蓼 (Polygonum macrophyllum)、中亞早熟禾 (Poalitwinowiana)、麻花艽 (Gentiana straminea)、弱小火絨草 (Leontopodium pusillum) 等。

    雪層杜鵑 (Rhododendron nivale)+櫻草杜鵑 (Rhododendron primuliflorum)+高山嵩草群叢 (b):共包括19個(gè)樣方,主要分布海拔為4 400~4 800 m。在群叢水平上,以雪層杜鵑為優(yōu)勢種,次優(yōu)勢種為高山嵩草,灌木層的次優(yōu)勢種為櫻草杜鵑;主要的伴生種有釘柱委陵菜、藏西嵩草、圓穗蓼、弱小火絨草、高山唐松草 (Thalictrum alpinum)、甘松 (Nardostachys jatamansi) 等。

    絲穎針茅 (Stipa capillacea) +高山嵩草群叢 (c):共包括5個(gè)樣方,主要分布海拔為4 300~4 400 m。在群叢水平上,絲穎針茅為優(yōu)勢種,次優(yōu)勢種為高山嵩草;伴生種有釘柱委陵菜、肉果草、藏橐吾 (Ligularia rumicifolia)、西伯利亞蓼 (Polygonum sibiricum) 等。

    藏西嵩草+高山嵩草群叢 (d):包括8個(gè)樣方,主要分布海拔為4 200~4 400 m。在群叢水平上,藏西嵩草為優(yōu)勢種,次優(yōu)勢中為高山嵩草;主要的伴生種有絲穎針茅、二裂委陵菜(Potentilla bifurca)、青藏苔草、肉果草等。

    草沙蠶+藏橐吾群叢 (e):包括8個(gè)樣方,主要分布在海拔3 900~4 000 m的干暖河谷底部。在群叢水平上,以草沙蠶為優(yōu)勢種,次優(yōu)勢種為藏橐吾;主要的伴生種有筆直黃芪 (Astragalus strictus)、肉果草、絲穎針茅等。

    白草+小藍(lán)雪花群叢 (f):包括3個(gè)樣方,主要分布在海拔3 900~4 100 m的溫度較高、較干旱的環(huán)境中。在群叢水平上,以白草 (Pennisetum centrasiaticum) 為優(yōu)勢種;主要的伴生種有筆直黃芪、二裂委陵菜、狼毒 (Stellera chamaejasme)、美葉藏菊 (Stellera chamaejasme)、蜜花毛果草 (Lasiocaryum densiflorum) 等。

    蕨麻委陵菜 (Potentilla anserina) +二裂委陵菜群叢 (g):包括4個(gè)樣方,分布海拔在4 100 m左右。在群叢水平上,以蕨麻委陵菜為優(yōu)勢種;主要的伴生種有青藏苔草、平車前 (Plantago depressa)、西伯利亞蓼、蜜花毛果草等。

    2.2 虧組山植物物種和樣方排序

    2.2.1 植物物種PCA排序 PCA分析結(jié)果 (表1) 顯示,排序軸1的特征值λ=0.23,物種數(shù)據(jù)累計(jì)變異百分率為22.9%,排序軸1解釋22.9%的物種變化;排序軸2的特征值λ=0.13,物種數(shù)據(jù)累計(jì)變異百分率為12.8%,與排序軸1共同解釋了物種組成變化為35.7%;排序軸3和排序軸4分別解釋了物種組成變化的8.2%、5.9%。這表明排序軸1和排序軸2為最能反映植物群落物種組成變化的影響因子梯度。

    物種PCA排序結(jié)果 (圖2)顯示,高山嵩草、釘柱委陵菜、藏西嵩草、美麗風(fēng)毛菊 (Saussurea pulchra)、線葉龍膽 (Gentiana farreri)、緣毛紫菀 (Aster souliei) 圓穗蓼等大多排在軸1的左下方,這些物種屬高山草甸成分,分布于高海拔區(qū)域;而拉薩狗娃花 (Heteropappus bowerii)、藏橐吾、密花毛果草 (Lasiocaryum densiflorum)、二裂委陵菜、 白草、草沙蠶等大多排在軸1的右上方,這些物種通常出現(xiàn)在海拔較低、濕度較小的干暖河谷區(qū)域。絲穎針茅、肉果草、西藏蕁麻 (Urtica tibetica) 等排在軸2的最上方,這些物種通常出現(xiàn)在干擾和退化程度較大的區(qū)域;而多刺綠絨蒿 (Meconopsis horridula)、櫻草杜鵑、優(yōu)越虎耳草 (Saxifraga egregia) 等排在軸2的下方,這些物種通常出現(xiàn)在濕度大、海拔較高的地方。

    2.2.2 植物群落樣方PCA排序 虧組山65個(gè)植物樣方的PCA排序結(jié)果 (圖3) 與TWINSPAN分類分析結(jié)果 (圖1) 存在一定的差異。在PCA排序分析中將樣方劃分為3個(gè)類型,在TWINSPAN分類中分為7個(gè)類型,TWINSPAN分類中的a和b群叢在PCA排序中合為組成類型 (Ⅰ);c和d群叢合為組成類型 (Ⅱ);e、f和g群叢合為組成類型 (Ⅲ)。從PCA劃分的3個(gè)類型內(nèi)植物群叢的組成來看,Ⅰ、Ⅱ和Ⅲ分別為灌叢草甸、高寒草甸和溫性草原的群叢組成成分。

    2.2.3 植物群落樣方RDA排序 從虧組山植物樣方RDA分析結(jié)果 (圖4,表2) 可以看出,RDA排序軸1的特征值λ1=0.16,物種與影響因子相關(guān)系數(shù)r=0.89,軸1解釋了物種組成與影響因子之間關(guān)系的16.4%;軸2 的特征值λ2=0.03,物種與影響因子的相關(guān)系數(shù)r=0.72,與軸1累計(jì)解釋了物種組成與影響因子之間關(guān)系的13.3%;軸3的特征值λ3=0.14,物種與影響因子的相關(guān)系數(shù)r=0,與軸3累計(jì)解釋了物種組成與影響因子之間關(guān)系的33.1%;軸4的特征值λ4=0.09,物種與影響因子的相關(guān)系數(shù)r=0,與軸3累計(jì)解釋了物種組成與影響因子之間關(guān)系的41.6%。第1軸與海拔和坡度相關(guān)系數(shù)分別為-0.88和-0.68,呈負(fù)相關(guān);軸2與坡度和海拔相關(guān)系數(shù)分別為-0.46和0.07,與坡度負(fù)相關(guān),與海拔正相關(guān),但相關(guān)程度很小;軸3和軸4與環(huán)境梯度沒有相關(guān)性。

    從各因子之間的相關(guān)關(guān)系結(jié)果 (表3) 可以看出,在所調(diào)查的5個(gè)因子之間存在一定的相關(guān)性。其中,海拔與總蓋度 (r=0.71)、分種蓋度 (r=0.54)、生物量相關(guān)性不顯著;坡度與總蓋度 (r=0.55)、分種蓋度 (r=0.41)、生物量 (r=0.27) 正顯著相關(guān)。

    2.3 物種豐富度分布格局

    通過回歸分析結(jié)果 (圖5:a-e)發(fā)現(xiàn),虧組山植物物種豐富度與海拔、分種蓋度存在顯著相關(guān) (P<0.01),豐富度從低海拔到中、高海拔 (3 900~4 500 m) 逐漸上升,在4 600~4 700 m達(dá)到最高點(diǎn)之后開始緩慢下降,呈中間高、兩邊低的“單峰”格局,且該趨勢解釋量R2=0.447,能夠解釋44.7%的物種豐富度變化 (圖5:a)。植物物種豐富度隨著分種蓋度先逐漸上升,當(dāng)分種蓋度在80%~100%之間時(shí)豐富度達(dá)到最大值,在分種蓋度超過100%之后物種豐富度開始下降,也呈單峰變化趨勢,并且該趨勢能夠解釋36.1%的物種豐富度變化 (圖5:c)。植物物種豐富度與坡度和總蓋度存在顯著相關(guān) (P<0.01),呈直線上升趨勢 (圖5:b, d),并且該趨勢能夠分別解釋41.6%和33.3%的物種豐富度變化。植物物種豐富度與生物量存在顯著相關(guān)關(guān)系 (P<0.01),呈先上升后穩(wěn)定的變化趨勢 (圖5:e),且每樣方內(nèi)地上生物量在10~30 g時(shí)物種豐富度值較高。

    3 討論與結(jié)論

    3.1 群落特征和植被類型

    TWINSPAN分類將虧組山65個(gè)樣方分為7個(gè)組,分別為7種群叢。PCA樣方排序分析將虧組山植物群落分為Ⅰ、Ⅱ和Ⅲ 3個(gè)組。對虧祖山的65個(gè)樣方用不同的分類方法進(jìn)行分析時(shí)發(fā)現(xiàn),兩種方法的分類結(jié)果有些差異。但是,從PCA排序進(jìn)行分類的3個(gè)組和TWINSPAN分類的7個(gè)組間進(jìn)行對比時(shí)發(fā)現(xiàn)兩結(jié)果間卻沒有很大的沖突。如PCA分類的Ⅰ組里包括TWINSPAN分類里的兩個(gè)群叢(a和b),為典型的灌叢草甸,Ⅱ?yàn)楦呱讲莸椋╟和d),Ⅲ為高山草原(e,f和g) 三種不同的草地群落類型。

    將對不同類型群叢所對應(yīng)的樣方信息進(jìn)行分析發(fā)現(xiàn),類型Ⅲ分布于海拔3 900~4 100 m,該區(qū)域溫度較高,蒸發(fā)量大,人和動物的干擾較大的河谷底部;類型Ⅱ分布在海拔4 200~4 400 m之間,該區(qū)域水熱條件處于比較理想的狀態(tài);類型Ⅰ分布的海拔范圍為4 400~5 100 m,這一植被類型為山體的主要植被類型,是小型灌叢分布較普遍的區(qū)域??梢?,在不同的區(qū)域,生境與環(huán)境條件不同,植被類型也不同。對每個(gè)植被類型樣方對應(yīng)的物種信息進(jìn)行分析發(fā)現(xiàn),Ⅰ內(nèi)總計(jì)68物種,共有種為38種;Ⅱ內(nèi)總計(jì)45個(gè)物種,共有種為20種;Ⅲ內(nèi)總計(jì)43物種,共有種為10種。在PCA分析中,Ⅰ類型內(nèi)a和b之間很難區(qū)分開來,呈連續(xù)狀的分布,Ⅱ和Ⅲ也是如此。高山草地植物群落在垂直樣帶上的分布有明顯的組成差異,但不同的群落類型之間有過渡區(qū),群落之間存在連續(xù)性分布;在相對同一個(gè)樣帶內(nèi),由于人為、地形和地表徑流等因子的影響,出現(xiàn)草地植物群落的斑塊化,使得植物群落在不同生境間呈間斷性分布??梢姡琍CA和TWINSPAN分類很好地揭示了植物群落在不同環(huán)境梯度上的連續(xù)性與間斷性的分布特征。

    物種PCA排序分析結(jié)果表明,軸1的左下方分布著以高海拔高山草甸物種成分為主,而軸1的右上方是以低海拔的干暖河谷出現(xiàn)的物種為主,結(jié)合樣方的PCA排序結(jié)果,可見軸1明顯的代表海拔梯度;軸2上方分布著絲穎針茅、肉果草、西藏蕁麻等通常出現(xiàn)在干旱生境、干擾和退化較大的地方出現(xiàn)的物種,而軸2下方分布著多刺綠絨蒿、櫻草杜鵑、優(yōu)越虎耳草等濕度較大生境中出現(xiàn)的物種,可見軸2反映出以濕度為主的干擾等復(fù)合的環(huán)境梯度。這與拉多等(2016) 在拉薩河谷流域做的植被研究結(jié)果一致。綜上所述,能夠推斷海拔、濕度和干擾梯度是影響該區(qū)域群落分布的主導(dǎo)因素。

    RDA冗余分析結(jié)果顯示,海拔梯度和坡度的變化對區(qū)域植物群落和物種分布的影響。其中,海拔對群落和物種分布格局的影響比坡度大,在藏北高原草地植物群落研究(王景升等,2016)中指出,濕度是影響該區(qū)域群落分布的主導(dǎo)因素;楊小林等(2010)在對拉薩河流域半干旱河谷地帶植物群落研究和姚帥臣等(2018) 對同區(qū)域草地植物群落研究表明,海拔是影響群落分布的主導(dǎo)環(huán)境因子之一。除了海拔和濕度以外,還在本次研究中發(fā)現(xiàn),坡度對群落組成與物種分布也有很重要的影響。坡度對群落的影響體現(xiàn)在不同的坡度所涉及的土壤厚度、土壤水分含量和單位面積內(nèi)接受的能量等的差異影響著群落組成與物種的分布(宋永昌,2001)。

    3.2 物種豐富度分布格局

    植物物種豐富度在垂直樣帶上隨環(huán)境梯度的變化而變化,由于研究區(qū)域、海拔跨度大小、坡度等環(huán)境因子的不同,因此呈現(xiàn)不同的結(jié)果 (彭閃江等,2003; Rahbek, 2005)。本研究表明,卡孜鄉(xiāng)虧組山植物物種豐富度隨海拔梯度大致呈中間高、兩邊低的單峰變化格局。與藏北高寒草地植物群落物種多樣性沿海拔梯度的分布格局類似 (段敏杰等,2011);在拉薩河谷山地灌叢草地物種多樣性隨海拔變化特征中也得到了相同的變化規(guī)律 (羅黎鳴等,2016)。低海拔帶 (3 900~4 200 m),由于溫度高、降雨量低、蒸散量高,在一定程度上,植物的生長主要受濕度的限制,且低海拔處還受人為干擾和牲畜放牧等影響較大,出現(xiàn)了筆直黃芪、藏橐吾、車前草等,因此植被稀疏、群落組成較為單一、植被呈現(xiàn)明顯的退化狀態(tài)。加之,在低海拔區(qū)域高原鼠兔的干擾強(qiáng)烈,不僅改變了高寒草甸植物群落內(nèi)的物種數(shù),而且還改變了原有物種在植物群落內(nèi)的地位和作用 (金少紅等,2017)。中海拔區(qū)域 (4 300~4 800 m),物種豐富度逐漸上升且達(dá)到最大值,這可能與該區(qū)域的水熱條件與灌叢的分布有關(guān),中海拔段相對于低海拔和高海拔區(qū)域,水熱條件相對良好,給植物提供了相對理想的生長環(huán)境,并且該區(qū)域分布有小型灌木叢,灌叢既能隔離牲畜的采食和踩踏,間接地對灌叢下的植物起到了保護(hù)作用,又能有效地降低風(fēng)速以及很好地截留掉落物,從而促進(jìn)土壤有機(jī)質(zhì)含量及其他物質(zhì)在灌叢下的積累 (魏巍,2013;瞿王龍等,2015),導(dǎo)致物種豐富度的增高。在做調(diào)查時(shí)發(fā)現(xiàn),中海拔段存在徑流水,與低海拔區(qū)域相比濕度較大,土壤濕度作為植物生長的主要環(huán)境因子,對植被生長和分布起著很重要的作用 (武建雙等,2012)。高海拔區(qū)域 (4 900~5 100 m),物種豐富度開始有降低的趨勢,隨著海拔的升高,溫度逐漸成為草地植物生長的主要限制因子 (徐廣平等,2005),在本研究區(qū)域的最高海拔5 100 m處,物種豐富度雖有下降的趨勢,但維持在較高的水平,這可能與樣方地理位置有關(guān),最高的樣點(diǎn)所處位置地勢比較平緩,土壤比較厚,并且被周圍的小山丘圍繞著,可能形成了特殊的生境,導(dǎo)致維持較高的物種豐富度。

    從物種豐富度隨總蓋度和分種蓋度變化趨勢可以看出,物種豐富度與總蓋度呈正顯著相關(guān)關(guān)系,與分種蓋度呈單峰變化趨勢,這與唐志堯和方精云(2004)以及拉瓊等(2014)的研究結(jié)果一致,并且方精云等(2004)認(rèn)為這種關(guān)系的出現(xiàn)可能與生境質(zhì)量有關(guān)。植被蓋度作為衡量植物群落保持水土功能 (郭忠升,2000) 和地面植物群落直觀定量化的功能指標(biāo),能夠反映區(qū)域生境質(zhì)量的好壞 (Harris A et al., 2014)。在研究區(qū)域植物物種豐富度隨總蓋度一直上升,并且還有進(jìn)一步的上升趨勢,而物種豐富度較高的值都集中在蓋度80%~100%,說明在這個(gè)蓋度范圍內(nèi)物種以一種最合理的資源分配方式存在,當(dāng)分種蓋度超過100%時(shí)物種豐富度開始下降,可能是由于物種之間存在對光、水分、土壤養(yǎng)分等環(huán)境資源的競爭,而競爭能力小的物種就可能開始消失 (Bonser & Reader, 1995),從而導(dǎo)致物種豐富度減少。有關(guān)研究報(bào)道,植被蓋度不僅是土壤侵蝕與水土流失的主要因子,而且還在區(qū)域尺度上影響著土壤水分的平衡 (趙文智,2002;Xin et al., 2004);植被通過根系吸水和水力提升起到水分調(diào)節(jié)作用 (李小雁,2011)。在不同等級植被蓋度對土壤水分的影響研究中發(fā)現(xiàn),植被蓋度過大,土壤過度耗水,導(dǎo)致植被的嚴(yán)重演替與退化(孫禎元,1984)。

    坡度在某種程度上決定了植物群落的分布格局 (劉金根和薛建輝,2010)。根據(jù)喻慶國等 (2007) 對坡度等級的劃分,本研究區(qū)域在平坡至陡坡范圍內(nèi) (3°~34°),坡度作為影響水土流失的主導(dǎo)因子 (曹梓豪等,2017),其大小會影響土壤水分含量和植物在單位面積接受的能量 (車明軒等,2016; 劉建華等,2017),并且在垂直樣帶上,影響著土壤養(yǎng)分含量,進(jìn)而成為影響草本和灌木植物分布的主要地形因子 (袁鐵象等,2014)。本研究區(qū)域,由于各樣方之間和各樣帶之間的坡度差異較大,導(dǎo)致小環(huán)境條件的差異,從而影響著植物物種豐富度。本研究中,物種豐富度隨坡度變化呈直線上升趨勢,坡度較大的地方海拔相對較高,受到人畜干擾較小,并且在坡度較大的區(qū)域,低矮灌叢分布相對較多,灌叢具有一定的增加土壤肥力的作用 (Wezel et al., 2000),隨著坡度的增加灌叢的多樣性也增加 (韓杰等,2016),灌叢的出現(xiàn)為流蘇虎耳草 (Saxifraga wallichiana)、甘松 (Nardostachys chinensis)、多刺綠絨蒿 (Meconopsis horridula) 等喜濕喜陰類叢下伴生植物種提供了適合的生境條件,從而對物種豐富度的提高具有積極作用。坡度較小的區(qū)域與村民居住地較近,牲畜及人為活動較為頻繁,部分區(qū)域出現(xiàn)草地退化,禾草和雜類草優(yōu)勢較為明顯,加之鼠兔的分布明顯比坡度大的區(qū)域多,且水源的分布主要集中在低中海拔段(3 900~4 500 m),飲水點(diǎn)附近作為高山草地牲畜干擾較強(qiáng)的區(qū)域,其草地植物物種多樣性與干擾強(qiáng)度呈正相關(guān) (唐明坤等,2011;呂曾哲舟,2019)。

    植物生物量作為生態(tài)系統(tǒng)的功能指標(biāo)和獲取能量的集中體現(xiàn),在生態(tài)系統(tǒng)的結(jié)構(gòu)和功能的形成上有重要作用。目前關(guān)于生物量與植物物種豐富度之間的關(guān)系研究結(jié)果較為多樣,在群落分布格局研究中,植物物種多樣性的變化由于群落類型的不同,其功能群以及植物生活型和種間競爭等差異對生物量有不同的影響 (陳加際等,2018)。本研究中,物種豐富度隨生物量先增后穩(wěn)定的變化趨勢,生物量達(dá)到一定的值后,物種豐富度變化緩慢。這與楊元合等(2004)的研究結(jié)果一致,研究認(rèn)為是生境條件差異及物種間的共生作用強(qiáng)所導(dǎo)致的。在本研區(qū)域之所以出現(xiàn)這一變化趨勢可能是由于每個(gè)樣帶海拔和坡度等環(huán)境因子的不同而導(dǎo)致不同的小生境,且本研究區(qū)域?yàn)楫?dāng)?shù)刈匀环拍恋?,每個(gè)區(qū)域的放牧強(qiáng)度也有所差異。此外,在垂直樣帶上草地有3種群落類型和7個(gè)群叢,不同的群落類型間,種間競爭和植物物種的組成等有所差異(劉立斌等,2018),其生物量也存在差異。如在樣方內(nèi)出現(xiàn)藏橐吾、黃苞南星 (Arisaema flavum) 等株高較高、葉片面積大和莖稈粗的物種,會直接影響樣方內(nèi)整體植物的生物量,因此生物量大的區(qū)域物種豐富度不一定高。隨著放牧強(qiáng)度的增加,植物生物量會減少,但是物種豐富度的變化還得考慮具體的放牧壓大小。高寒草地通過物種之間的信息與其功能間的聯(lián)系來適應(yīng)變化的環(huán)境,從而維持生產(chǎn)力的相對穩(wěn)定 (朱桂麗等,2004)。因此,這兩種變量間不是直接的因果關(guān)系,而是由環(huán)境因子、草地類型、放牧干擾強(qiáng)度等因子共同作用導(dǎo)致。

    虧組山體植被在群落水平上可分為3個(gè)類型,各類型間具有連續(xù)性分布特點(diǎn),群叢水平上分為7個(gè)類型,具有間斷性分布特點(diǎn);群落在垂直樣帶上的分布格局主要受海拔的影響,其次為坡度。在研究區(qū)植物物種豐富度隨海拔變化呈單峰分布格局,與坡度呈正相關(guān)。植物物種豐富度對兩種蓋度變化的響應(yīng)存在一定的差異,其與分蓋度呈單峰分布格局,與總蓋度呈正相關(guān),植物物種豐富度與生物量呈正相關(guān)。山體植物群落組成,物種豐富度及其分布,不僅受環(huán)境因子的影響,而且還與總蓋度、分蓋度和生物量等生物因子之間存在顯著的關(guān)系。由于本研究中未考慮到土壤養(yǎng)分、降雨量、坡向等對區(qū)域植被的影響,因此在PCA排序分析推斷的濕度因子未能在RDA冗余分析中得到驗(yàn)證,有待于今后進(jìn)一步地研究。

    致謝 感謝中國科學(xué)院地理研究所國家自然科學(xué)基金面上項(xiàng)目(31570460)和“復(fù)旦大學(xué)鐘揚(yáng)教授基金”對本研究的資助,以及卡孜鄉(xiāng)當(dāng)?shù)剞r(nóng)牧民對本研究的大力支持與配合,特此致謝。

    參考文獻(xiàn):

    BONSER S, READER RJ, 1995. Plant competition and herbivory in relation to vegetation biomass? [J]. Ecology,76(7): 2176-2183.

    CAO ZH, ZHAO QH, DING SY, et al., 2017. Effect of slope gradient and vegetation cover on sediment yielding characteristics of the riparian slope [J]. J Nat Resour,

    32(11): 1892-1904.? [曹梓豪, 趙清賀, 丁圣彥, 等, 2017. 坡度和植被蓋度對河岸坡面侵蝕產(chǎn)沙特征的影響 [J]. 自然資源學(xué)報(bào), 32(11): 1892-1904.]

    CHE MX, GONG YB, NAEEM KM, et al., 2016. Effects of different rainfall intensities and slope gradients on water and soil conservation of straw mulch [J]. J Soil Water Conserv,30(2): 131-135.? [車明軒, 宮淵波, MUHAMMAD NAEEM KHAN, 等, 2016. 不同雨強(qiáng)、坡度對秸稈覆蓋保持水土效果的影響 [J]. 水土保持學(xué)報(bào), 30(2): 131-135.]

    CHEN JJ, CHANG SH, WANG ZF, et al., 2018, The relationship between plant species diversity and above-ground biomass in a typical steppe in Abaga Banner [J]. Pratacul Sci, 35(9): 2068-2078.? [陳加際, 常生華, 王召鋒, 等, 2018. 阿巴嘎旗典型草原植物物種多樣性與地上生物量的關(guān)系 [J]. 草業(yè)科學(xué), 35(9): 2068-2078.]

    Comprehensive Scientific Expedition Team of Chinese Academy of Sciences in Qinghai-Tibet Plateau, 1988. Tibetan Vegetation [M]. Beijing: Science Press: 90-228.? [中國科學(xué)院青藏高原綜合考察隊(duì), 1980. 西藏植被 [M]. 北京:科學(xué)出版社: 90-228.]

    DUAN MJ, GAO QZ, GUO YQ, et al., 2011. Effects of altitude on plant productivity and species diversity in alpine meadows of northern Tibet [J]. Pratacul Sci, 28(10):1845-1850.? [段敏杰, 高清竹, 郭亞奇, 等, 2011. 藏北高寒草地植物群落物種多樣性沿海拔梯度的分布格局 [J]. 草業(yè)科學(xué), 28(10): 1845-1850.]

    GUO ZS, 2000. Three coverages in soil and water conservation vegetation construction: Potential coverage, critical coverage and effective coverage [J]. Bull Soil Water Conserv,(2): 60-62.? [郭忠升, 2000. 水土保持植被的有效蓋度、臨界蓋度和潛勢蓋度 [J]. 水土保持通報(bào), (2): 60-62.]

    HARRIS A, CARR A, DASH J, et al., 2014. Remote sensing of vegetation cover dynamics and resilience across southern Africa [J]. Int J Appl Earth Obs Geoinf, 28: 131-139.

    JIN L, DU XG, HOU FJ, et al., 2007. Effects of hillside slopes on plant community succession of re-natural farmland in the Loess, Plateau [J]. Pratacul Sci, 24(7): 66-71.? [金樑, 杜曉光, 侯扶江, 等, 2007. 黃土高原山地坡度對退耕農(nóng)田生態(tài)系統(tǒng)自然植被演替初期的影響 [J]. 草業(yè)科學(xué),24(7): 66-71.]

    JIN SH, LIU T, PANG XP, et al.,(2017). Effects of plateau pika(Ochotona crzoniae) disturbances on plant species diversity and above ground plant biomass in a Kobresia pygmaea meadow in the Qinghai Lake Region [J]. Acta Pratacul Sin, 26(5): 29-39.? [金少紅, 劉彤, 龐曉攀, 等, 2017. 高原鼠兔干擾對青海湖流域高山嵩草草甸植物多樣性及地上生物量的影響 [J]. 草業(yè)學(xué)報(bào), 26(5):29-39.]

    JIN YM, TIAN FY, WU JX, et al., 2018. Effects of indoor feeding on weight gain and wool mass of Pengbo semi-fine wool sheep during winter-spring period in Lhasa River Valley [J]. J Plat Agric, 2(4):347-353.? [金艷梅, 田發(fā)益, 武俊喜, 2018. 冬春季舍飼對拉薩河谷彭波半細(xì)毛羊增重和產(chǎn)毛量的影響 [J]. 高原農(nóng)業(yè), 2(4): 347-353.]

    JI D, MIROSLAV , 2002. Altitudinal changes in composition and structure of mountain temperate vegetation: a case study from the Western Carpathians? [J]. Plant Ecol, 158(2): 201-221.

    LA D, ZHANG YJ, LIU J, 2016. Spatial distribution patterns and environmental interpretation of plant species richness in the Lhasa Valley, Tibet? [J]. Acta Pratacul Sci, 25(10): 202-211.? [拉多, 張燕杰, 劉杰, 等, 2016. 拉薩河谷植物物種豐富度空間分布格局及其環(huán)境解釋 [J]. 草業(yè)學(xué)報(bào),25(10): 202-211.]

    LA Q, ZHA XCR, ZHU WD, et al., 2014. Plant species-richness and association with environmental factors in the riparian zone of the Yarlong Zangbo River of Tibet, China? [J]. Biodivers Sci, 22(3): 337-347.? [拉瓊, 扎西次仁, 朱衛(wèi)東, 等, 2014. 雅魯藏布江河岸植物物種豐富度分布格局及其環(huán)境解釋 [J]. 生物多樣性, 22(3): 337-347.]

    LI SQ, ZHANG JT, 2007. Ecological analysis of subalpine meadow in Yunding Mountain, Shanxi? [J]. Geogr Res,26(1): 83-90.? [李素清, 張金屯, 2007. 山西云頂山亞高山草甸群落生態(tài)分析 [J]. 地理研究, 26(1): 83-90.]

    LI WH, GAN ZZB, CAO XJ, et al., 2017. Effects of altitude on plant productivity and species diversity in alpine meadows of northern Tibet? [J]. Acta Pratacul Sin, 26(9): 200-207.? [栗文瀚, 干珠扎布, 曹旭娟, 等, 2017. 海拔梯度對藏北高寒草地生產(chǎn)力和物種多樣性的影響 [J]. 草業(yè)學(xué)報(bào), 26(9): 200-207.]

    LI XY, 2011. Soil-vegetation-hydrology coupling, response and adaptation mechanism in arid areas? [J]. Chin Sci Sin(Terrae), 41(12): 1721-1730.? [李小雁, 2011. 干旱地區(qū)土壤-植被-水文耦合、響應(yīng)與適應(yīng)機(jī)制 [J]. 中國科學(xué):地球科學(xué), 41(12): 1721-1730.]

    LIU GC, HUANG YX, WANG QG, et al., 2018. Effects of environmental factors on plant species diversity: Research progress? [J]. Chin Agric Sci Bull, 34(13): 83-89.? [劉冠成, 黃雅曦, 王慶貴, 等, 2018. 環(huán)境因子對植物物種多樣性的影響研究進(jìn)展 [J]. 中國農(nóng)學(xué)通報(bào), 34(13):83-89.]

    LIU JG, XUE JH, 2010. Impacts of slope gradient on the characteristics of Vetiveria zizanioides community in early-stage? [J]. Grassl Turf, 30(2): 56-61.? [劉金根, 薛建輝, 2010. 坡度對香根草護(hù)坡植物群落早期特征的影響 [J]. 草原與草坪, 30(2): 56-61.]

    LIU JH, LI ML, SU JB, et al., 2017. Effects of soil organic matter content and terrain slope on soil moisture, crop yield [J]. Hunan Agric Sci,(1): 16-18.? [劉建華, 李銘亮, 蘇劍波, 等, 2017. 有機(jī)質(zhì)含量、地形坡度對土壤保水效果及作物產(chǎn)量的影響 [J]. 湖南農(nóng)業(yè)科學(xué),(1): 16-18.]

    LIU LB, ZHONG QL, NIJ, et al., 2018. Allometric function-based root biomass estimate of woody plants in a karst evergreen and deciduous broadleaf and mixed forest in central Guizhou Province, southwestern China? [J]. Acta Ecol Sin, 38(24): 8726-8732.? [劉立斌, 鐘巧連, 倪健, 等, 2018. 基于生物量回歸方程估算黔中喀斯特常綠落葉闊葉混交林木本植物的根系生物量 [J]. 生態(tài)學(xué)報(bào),38(24): 8726-8732.]

    LIU ZZZ, HUANG XX, WANG XY, et al., 2019. Disturbance pattern of alpine meadow in Yak Meadow Park, Jade Dragon Mountain [J]. J Natur Resour, 34(6): 1223-1231.? [呂曾哲舟, 黃曉霞, 王琇瑜, 等, 2019. 玉龍雪山牦牛坪高山草甸的干擾格局分析 [J]. 自然資源學(xué)報(bào),34(6): 1223-1231.]

    LUO LM, WU JS, YU CQ, et al., 2016. A comparative study of methods for monitoring plant diversity in montane shrub grassland in the Lhasa River Basin? [J]. Acta Pratacul Sin, 25(3): 22-31.? [羅黎鳴, 武建雙, 余成群, 等, 2016. 拉薩河谷山地灌叢草地植物多樣性監(jiān)測方法的比較研究 [J]. 草業(yè)學(xué)報(bào), 25(3): 22-31.]

    PENG SJ, HUANG ZY, XU GL,et al., 2003. Effects of habitat heterogeneity on forest community diversity in Dinghushan Biosphere Reserve? [J]. Guihaia, 23(5): 391-398.? [彭閃江, 黃忠良, 徐國良, 等, 2003. 生境異質(zhì)性對鼎湖山植物群落多樣性的影響 [J]. 廣西植物, 23(5): 391-398.]

    QU WL, YANG XP, ZHANG CT, et al., 2015. Shrub-mediated “fertile island” effects in arid and semi-qrid grassland? [J]. Acta Pratcacul Sin, 24(4): 201-207.? [瞿王龍, 楊小鵬, 張存濤, 等, 2015. 干旱、半干旱地區(qū)天然草原灌木及其肥島效應(yīng)研究進(jìn)展 [J]. 草業(yè)學(xué)報(bào), 24(4): 201-207.]

    RAHBEK C, 2005. The role of spatial scale and the perception of large-scale species-richness patterns? [J]. Ecol Lett,8(2): 224-239.

    SONG YC, 2001. Vegetation Ecology? [M]. Shanghai: East China Normal University Press: 171-176.? [宋永昌, 2001. 植被生態(tài)學(xué) [M]. 上海:華東師范大學(xué)出版社:171-176.]

    SONG ZH, 2013. Effect of different altitudes on kobresia pygmaea community? [J]. J Anhui Agric Sci, 41(11):4709-4711.? [宋長宏, 2013. 不同海拔梯度對小嵩草植物群落的影響 [J]. 安徽農(nóng)業(yè)科學(xué), 41(11): 4709-4711.]

    SUN ZY, 1984. Preliminary discussion on the relationship between vegetation coverage and water in Yulin sandy land? [J]. J Beijing For Univ,(3): 29-35.? [孫禎元, 1984. 對榆林沙地植被蓋度和水分關(guān)系的初步探討 [J]. 北京林學(xué)院學(xué)報(bào), (3): 29-35.]

    TANG MK, LI MF, ZHAO J, et al., 2011. Disturbance pattern of alpine meadow in Yak Meadow Park, Jade Dragon Mountain [J]. Guihaia, 31(6): 775-781.? [唐明坤, 李明富, 趙杰, 等, 2011. 四川若爾蓋縣不同退化程度高寒草地群落比較研究 [J]. 廣西植物, 31(6): 775-781.]

    TANG ZY, FANG JY, 2004. A review on the elevational patterns of plant species diversity? [J]. Biodivers Sci, 12(1): 20-28.? [唐志堯, 方精云, 2004. 植物物種多樣性的垂直分布格局 [J]. 生物多樣性, 12(1): 20-28.]

    WANG CT, LONG RJ, SHI JJ, et al., 2005. Study on the change of vertical structure and productivity of artificial grassland community sowed in different year in alpine area? [J]. Chin J Grassl, 27(5): 16-21.? [王長庭,龍瑞軍, 施建軍, 等, 2005. 高寒地區(qū)不同建植期人工草地群落垂直結(jié)構(gòu)和生產(chǎn)力變化的研究 [J]. 中國草地學(xué)報(bào), 27(5): 16-21.]

    WANG J, 2013. Study on the relation between geological background and plant diversity of Gongbu nature reserve in Tibet? [D]. Chengdu: Chengdu University of Technology: 118.? [王娟, 2013. 西藏工布自然保護(hù)區(qū)地質(zhì)背景與植物多樣性關(guān)系研究 [D]. 成都: 成都理工大學(xué): 118.]

    WANG JS, YAO SC, PU Q, et al., 2016. Quantitative classification and ordination of grassland communities on the Northern Tibetan Plateau ?[J]. Acta Ecol Sin, 36(21): 6889-6896.? [王景升, 姚帥臣, 普窮, 等, 2016. 藏北高原草地群落的數(shù)量分類與排序 [J]. 生態(tài)學(xué)報(bào), 36(21): 6889-6896.]

    WEI W, 2013. Plant diversity pattern and interface characteristics of shrubland-grassland ecotone in Eastern Qilian Mountains? [D]. Gansu: Gansu Agricultural University: 63.? [魏巍, 2013. 東祁連山灌—草交錯(cuò)帶植物多樣性格局及界面因子特征研究 [D]. 甘肅: 甘肅農(nóng)業(yè)大學(xué): 63.]

    WEZEL A, RAJOT JL, HERBRIG C, et al., 2000. Influence of shrubs on soil characteristics and their function in sahelian agro-ecosystems in semi-arid Niger? [J]. J Arid Environ,44(4): 383-398.

    WU JS, LI XJ, SHEN ZX, et al., 2012. Species diversity distribution pattern of alpine grasslands communities along a precipitation gradient across Northern Tibetan Plateau? [J]. Acta Pratacul Sin, 21(3): 17-25.? [武建雙, 李曉佳, 沈振西, 等, 2012. 藏北高寒草地樣帶物種多樣性沿降水梯度的分布格局 [J]. 草業(yè)學(xué)報(bào), 21(3): 17-25.]

    XIN RL, FENG YM, HONG LX, et al., 2004, Long-term effects of revegetation soil water content of and dunes in arid region of Northern China? [J]. J Arid Environ, 57(1):1-16.

    XU GP, ZHANG DG, PU XP, et al., 2005. Two-dimensional polar ordination of plant community of alpine grassland in Eastern Qilian mountains? [J]. Grassl Turf,(3): 42-45.? [徐廣平, 張德罡, 蒲小鵬, 等, 2005. 東祁連山高寒草地植物群落二維極點(diǎn)排序 [J]. 草原與草坪, (3): 42-45.]

    YANG XL, ZHAO KT, MA HP, et al., 2010. Ecological studies on vegetation quantity in the Semi-Arid Valley region of Lasa? [J]. Sci Silv Sin, 46(10): 15-22.? [楊小林, 趙墾田, 馬和平, 等, 2010. 拉薩半干旱河谷地帶的植被數(shù)量生態(tài)研究 [J]. 林業(yè)科學(xué), 46(10): 15-22.]

    YANG YH, RAO S, HU HF, et al., 2004. Plant species richness of alpine grasslands in relation to environmental factors and biomass on the Tibetan Plateau? [J]. Biodivers Sci, 12(1): 200-205.? [楊元合, 饒勝, 胡會峰, 等, 2004.青藏高原高寒草地植物物種豐富度及其與環(huán)境因子和生物量的關(guān)系 [J]. 生物多樣性, 12(1): 200-205.]

    YAO SC, WANG JS, DING LB, et al., 2018. Quantitative classification and ordination of grassland communities in Lhasa River Valley? [J]. Chin Acta Ecol Sin,38(13): 4779-4788.? [姚帥臣, 王景升, 丁陸彬, 等, 2018. 拉薩河谷草地群落的數(shù)量分類與排序 [J]. 生態(tài)學(xué)報(bào),38(13):4779-4788.]

    YU QG, CAO SW, DENG XQ, et al., 2007. Study on the spatial distribution of vegetation along the gradients in Nuozhadu Nature Reserve of Yunnan Province? [J]. J Anhui Agric Sci, 35(28): 8994-8996.? [喻慶國, 曹順偉, 鄧喜慶, 等, 2007. 云南糯扎渡自然保護(hù)區(qū)植被沿坡度空間分異研究 [J]. 安徽農(nóng)業(yè)科學(xué), 35(28): 8994-8996.]

    YUAN TX, ZHANG HP, OU ZY, et al., 2014. Effects of topography on the diversity and distribution pattern of ground plants in karst montane forests in Southwest Guangxi, China [J]. Chin J Appl Ecol, 25(10): 2803-2810.? [袁鐵象, 張合平, 歐芷陽, 等, 2014. 地形對桂西南喀斯特山地森林地表植物多樣性及分布格局的影響 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 25(10): 2803-2810.]

    ZHANG XQ, LUO LM, JIN YM, et al., 2017. Changes in nutrients of dominant species on montane shrub grassland in the Lhasa River Basin? [J]. Chin J Grassl,39(3): 90-95.? [張曉慶, 羅黎鳴, 金艷梅, 等, 2017. 拉薩河谷山地灌叢草地優(yōu)勢種營養(yǎng)物質(zhì)含量變化規(guī)律 [J]. 中國草地學(xué)報(bào),39(3): 90-95.]

    ZHANG YF, ZHAO QH, DING SY, et al., 2017. Effects of slope gradient and vegetation coverage on hydrodynamic characteristics of overland flow on silty riparian slope [J]. Chin J Appl Ecol, 28(8): 2488-2498.? [張祎帆, 趙清賀,丁圣彥, 等, 2017. 坡度和植被蓋度對粉砂質(zhì)河岸帶坡面流侵蝕水動力特征的影響 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 28(8):2488-2498.]

    ZHAO HZ, LI WH, MA AJ, et al., 2010. Evaluation of ecological services value of mountain shrubby-grassland ecosystem in Lhasa-river valley region — A case study of Dazi County? [J]. Pratacul Sci, 27(12): 27-31.? [趙海珍, 李文華, 馬愛進(jìn), 等, 2010. 拉薩河谷山地灌叢草地生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)——以拉薩達(dá)孜縣為例 [J]. 草業(yè)科學(xué), 27(12):27-31.]

    ZHAO WZ, 2002. Impact of plantation on spatial heterogeneity of soil moisture in Horqin sandy land? [J]. Acta Pedol Sin,39(1): 107-113.? [趙文智, 2002. 科爾沁沙地人工植被對土壤水分異質(zhì)性的影響 [J]. 土壤學(xué)報(bào), 39(1): 107-113.]

    ZHENG W, ZHU JZ, PAN CD, 2010. Effect of grazing disturbance on plant functional group and community structure of meadow community in Kanas Scenic Area? [J]. Chin J Grassl, 32(1): 92-98.? [鄭偉, 朱進(jìn)忠, 潘存德, 2010. 放牧干擾對喀納斯草地植物功能群及群落結(jié)構(gòu)的影響 [J]. 中國草地學(xué)報(bào), 32(1): 92-98.]

    ZHONG XH, LIU SZ, WANG XD, et al., 2006. A research on the protection and construction of the state ecological safe shelter zone on the Tibet Plateau? [J]. Mt Res, 24(2): 129-136.? [鐘祥浩, 劉淑珍, 王小丹, 等, 2006. 西藏高原國家生態(tài)安全屏障保護(hù)與建設(shè) [J]. 山地學(xué)報(bào), 24(2):129-136.]

    ZHU GL, LI J, WEI XH, et al., 2017. Longitudinal patterns of productivity and plant diversity in Tibetan alpine grasslands? [J]. J Natur Resour, 32(2): 210-222.? [朱桂麗, 李杰, 魏學(xué)紅, 等, 2017. 青藏高寒草地植被生產(chǎn)力與生物多樣性的經(jīng)度格局 [J]. 自然資源學(xué)報(bào), 32(2): 210-222.]

    ZHU Y, KANG MY, LIU QR, et al., 2007. Relationship between biodiversity and aboveground biomass in alpine meadow on Mt. Helan, China? [J]. J Appl Environ Biol,13(6): 771-776.? [朱源, 康慕誼, 劉全儒, 等, 2007. 賀蘭山高山草甸生物多樣性和地上生物量的關(guān)系 [J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 13(6): 771-776.]

    (責(zé)任編輯 蔣巧媛)

    猜你喜歡
    環(huán)境因子群落結(jié)構(gòu)排序
    排序不等式
    恐怖排序
    節(jié)日排序
    刻舟求劍
    兒童繪本(2018年5期)2018-04-12 16:45:32
    海南省紅樹林區(qū)浮游動物多樣性的初步研究
    海南省紅樹林區(qū)大型底棲動物多樣性的初步研究
    阿克蘇地區(qū)杏李與歐洲李氣體交換特征對比研究
    桉樹人工林胸徑變化特征與環(huán)境因子的關(guān)系研究
    中國市場(2016年36期)2016-10-19 05:25:57
    秦皇島海域夜光藻種群密度與環(huán)境因子的關(guān)系
    濟(jì)南流域輪蟲群落結(jié)構(gòu)及多樣性的研究
    国产精品久久久久久久电影 | 国产精品av视频在线免费观看| 久久精品亚洲精品国产色婷小说| av专区在线播放| 国产精品国产高清国产av| 午夜免费成人在线视频| 国产一区二区在线av高清观看| 深夜精品福利| 成人永久免费在线观看视频| 国产免费男女视频| 日韩精品中文字幕看吧| 欧美+亚洲+日韩+国产| 国产主播在线观看一区二区| 在线观看美女被高潮喷水网站 | 天堂动漫精品| 天天添夜夜摸| 好男人在线观看高清免费视频| 中出人妻视频一区二区| 成人欧美大片| 亚洲中文字幕日韩| 欧美日本亚洲视频在线播放| 天堂√8在线中文| 亚洲七黄色美女视频| 男人的好看免费观看在线视频| 免费在线观看成人毛片| 成人18禁在线播放| 看免费av毛片| 午夜两性在线视频| 亚洲精品日韩av片在线观看 | 精品无人区乱码1区二区| 欧美国产日韩亚洲一区| 在线播放无遮挡| 国产精品一区二区三区四区免费观看 | 无限看片的www在线观看| av福利片在线观看| 国产精品影院久久| 午夜福利欧美成人| 琪琪午夜伦伦电影理论片6080| 国产高清有码在线观看视频| 亚洲最大成人手机在线| 久久久久性生活片| 熟女少妇亚洲综合色aaa.| 日韩大尺度精品在线看网址| 免费av观看视频| 久久精品国产综合久久久| 国产精品1区2区在线观看.| 亚洲精品成人久久久久久| 欧美日韩乱码在线| 性欧美人与动物交配| 国产伦一二天堂av在线观看| 欧美av亚洲av综合av国产av| 日本免费a在线| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 搡老岳熟女国产| 深爱激情五月婷婷| 欧美日韩亚洲国产一区二区在线观看| 日本熟妇午夜| 丰满的人妻完整版| 免费看十八禁软件| 窝窝影院91人妻| 久久久久久人人人人人| 99国产综合亚洲精品| 99热6这里只有精品| 脱女人内裤的视频| 国产激情偷乱视频一区二区| 搡老妇女老女人老熟妇| 国产精品精品国产色婷婷| 啪啪无遮挡十八禁网站| 一本一本综合久久| 国产成人欧美在线观看| 蜜桃亚洲精品一区二区三区| 日韩欧美国产在线观看| 波多野结衣巨乳人妻| 色老头精品视频在线观看| 国产视频一区二区在线看| 午夜免费成人在线视频| 宅男免费午夜| 亚洲成人久久爱视频| 午夜视频国产福利| 亚洲成人中文字幕在线播放| 免费搜索国产男女视频| 色尼玛亚洲综合影院| 一级黄色大片毛片| 国产黄a三级三级三级人| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 亚洲精品成人久久久久久| 色播亚洲综合网| 国产v大片淫在线免费观看| 亚洲国产欧美网| 超碰av人人做人人爽久久 | 国产精品永久免费网站| 男女那种视频在线观看| 搞女人的毛片| 我的老师免费观看完整版| 国产真实伦视频高清在线观看 | 18+在线观看网站| 亚洲精品一区av在线观看| 免费观看人在逋| 欧美一级毛片孕妇| 成人特级黄色片久久久久久久| 亚洲国产欧洲综合997久久,| 一级黄片播放器| 好看av亚洲va欧美ⅴa在| 91av网一区二区| 成人午夜高清在线视频| 真实男女啪啪啪动态图| 91字幕亚洲| 极品教师在线免费播放| 成人18禁在线播放| 欧美+亚洲+日韩+国产| 免费大片18禁| 久久伊人香网站| 久久香蕉精品热| 国产av不卡久久| 无人区码免费观看不卡| 久久久久国产精品人妻aⅴ院| 搡女人真爽免费视频火全软件 | 国语自产精品视频在线第100页| 欧美午夜高清在线| 操出白浆在线播放| 日韩国内少妇激情av| 久久婷婷人人爽人人干人人爱| 18禁裸乳无遮挡免费网站照片| 免费看日本二区| 日本三级黄在线观看| 亚洲欧美日韩卡通动漫| 免费在线观看成人毛片| 啦啦啦观看免费观看视频高清| 网址你懂的国产日韩在线| 午夜影院日韩av| 久久久久久久午夜电影| 俺也久久电影网| 欧美中文日本在线观看视频| 一级毛片女人18水好多| 精品电影一区二区在线| 日本免费a在线| 亚洲第一电影网av| 国产av在哪里看| 91久久精品国产一区二区成人 | 国产极品精品免费视频能看的| 亚洲内射少妇av| 欧美色欧美亚洲另类二区| 99国产精品一区二区三区| 亚洲美女视频黄频| 在线观看免费午夜福利视频| 欧美绝顶高潮抽搐喷水| 久久精品国产清高在天天线| 免费人成视频x8x8入口观看| 国产91精品成人一区二区三区| 老汉色∧v一级毛片| 国产不卡一卡二| 少妇裸体淫交视频免费看高清| 国产精华一区二区三区| 97超级碰碰碰精品色视频在线观看| 久久久精品大字幕| 在线观看免费视频日本深夜| 亚洲内射少妇av| 国产精品亚洲av一区麻豆| 黄色成人免费大全| 午夜福利免费观看在线| 中文在线观看免费www的网站| 一进一出抽搐gif免费好疼| 琪琪午夜伦伦电影理论片6080| 国产v大片淫在线免费观看| 国产亚洲精品久久久久久毛片| 国产黄a三级三级三级人| 久久久国产成人精品二区| 国产精品亚洲一级av第二区| 九九热线精品视视频播放| 校园春色视频在线观看| 日韩免费av在线播放| 少妇裸体淫交视频免费看高清| 级片在线观看| 国产精品影院久久| 亚洲精华国产精华精| 丁香六月欧美| 国产精品久久久久久精品电影| 久9热在线精品视频| 午夜免费男女啪啪视频观看 | 免费看光身美女| 国产高潮美女av| 亚洲激情在线av| 青草久久国产| 日本熟妇午夜| 亚洲无线在线观看| 久久国产精品人妻蜜桃| 在线视频色国产色| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| 少妇的逼好多水| 国产单亲对白刺激| 日韩欧美 国产精品| 91在线精品国自产拍蜜月 | 叶爱在线成人免费视频播放| 真实男女啪啪啪动态图| 免费看a级黄色片| 中国美女看黄片| 18+在线观看网站| www.色视频.com| 女警被强在线播放| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 国内久久婷婷六月综合欲色啪| 亚洲专区中文字幕在线| www.www免费av| 亚洲在线自拍视频| 老司机在亚洲福利影院| 99国产综合亚洲精品| 两人在一起打扑克的视频| 搡老熟女国产l中国老女人| 91av网一区二区| 舔av片在线| 亚洲人成电影免费在线| 成人18禁在线播放| 国产精品98久久久久久宅男小说| 90打野战视频偷拍视频| 国内久久婷婷六月综合欲色啪| 日韩欧美国产一区二区入口| 亚洲五月婷婷丁香| 亚洲最大成人中文| 最近视频中文字幕2019在线8| 99久久精品国产亚洲精品| 夜夜躁狠狠躁天天躁| 亚洲最大成人中文| 神马国产精品三级电影在线观看| 88av欧美| 母亲3免费完整高清在线观看| 可以在线观看毛片的网站| 男人舔女人下体高潮全视频| 国产成人aa在线观看| 国产乱人伦免费视频| 综合色av麻豆| 99精品在免费线老司机午夜| 丰满的人妻完整版| 精品99又大又爽又粗少妇毛片 | av福利片在线观看| 成人午夜高清在线视频| 看免费av毛片| 在线免费观看的www视频| 免费看日本二区| 国产高清视频在线播放一区| www.www免费av| 久久精品亚洲精品国产色婷小说| 国产成人福利小说| 黑人欧美特级aaaaaa片| 久久精品国产亚洲av香蕉五月| 国内精品久久久久久久电影| 国产99白浆流出| 国产精品久久久久久人妻精品电影| 国产av麻豆久久久久久久| 婷婷六月久久综合丁香| 老司机午夜福利在线观看视频| 1000部很黄的大片| 高清毛片免费观看视频网站| 色老头精品视频在线观看| 午夜福利免费观看在线| 国产激情欧美一区二区| 免费电影在线观看免费观看| 久久草成人影院| 一区二区三区高清视频在线| 18美女黄网站色大片免费观看| 久久精品亚洲精品国产色婷小说| 久久久久久久午夜电影| 99热精品在线国产| 中文字幕人妻丝袜一区二区| 亚洲人与动物交配视频| 尤物成人国产欧美一区二区三区| 国产麻豆成人av免费视频| 免费在线观看日本一区| 三级国产精品欧美在线观看| 日本精品一区二区三区蜜桃| 国产v大片淫在线免费观看| 国产乱人视频| 美女cb高潮喷水在线观看| 少妇熟女aⅴ在线视频| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 午夜福利欧美成人| 狂野欧美白嫩少妇大欣赏| 欧美日韩综合久久久久久 | 韩国av一区二区三区四区| 国产精品免费一区二区三区在线| 国产成人欧美在线观看| 男人舔奶头视频| 18美女黄网站色大片免费观看| 一本综合久久免费| 午夜福利免费观看在线| 淫妇啪啪啪对白视频| 亚洲精品成人久久久久久| 欧美日韩福利视频一区二区| 精品国产超薄肉色丝袜足j| 成人午夜高清在线视频| 国产亚洲精品一区二区www| 国产高清视频在线观看网站| 精品99又大又爽又粗少妇毛片 | 久久精品亚洲精品国产色婷小说| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟人妻熟丝袜美 | 国产久久久一区二区三区| 女人高潮潮喷娇喘18禁视频| 成人性生交大片免费视频hd| 午夜免费男女啪啪视频观看 | 嫩草影院精品99| 色播亚洲综合网| 少妇的丰满在线观看| 国产精品98久久久久久宅男小说| 婷婷六月久久综合丁香| 免费在线观看亚洲国产| 伊人久久大香线蕉亚洲五| 亚洲av免费高清在线观看| 婷婷精品国产亚洲av| 日本a在线网址| 好看av亚洲va欧美ⅴa在| 国产一区二区亚洲精品在线观看| 国内精品久久久久精免费| 香蕉av资源在线| 亚洲avbb在线观看| 欧美在线黄色| 婷婷丁香在线五月| 国产野战对白在线观看| 又爽又黄无遮挡网站| 女人高潮潮喷娇喘18禁视频| 亚洲五月天丁香| 国产成人av激情在线播放| 久久久久亚洲av毛片大全| 日韩欧美免费精品| 亚洲精品乱码久久久v下载方式 | 男女视频在线观看网站免费| 国产色爽女视频免费观看| 久久这里只有精品中国| 亚洲久久久久久中文字幕| 美女免费视频网站| 亚洲av电影不卡..在线观看| 男人舔女人下体高潮全视频| 亚洲人成网站在线播| 亚洲成人中文字幕在线播放| 婷婷丁香在线五月| 一区福利在线观看| 国产精品 欧美亚洲| 午夜福利欧美成人| 亚洲国产精品久久男人天堂| 久久亚洲精品不卡| 日本在线视频免费播放| 久久婷婷人人爽人人干人人爱| 亚洲avbb在线观看| 黄色女人牲交| 黄色丝袜av网址大全| 男女做爰动态图高潮gif福利片| 精品午夜福利视频在线观看一区| 美女免费视频网站| 亚洲男人的天堂狠狠| 在线观看一区二区三区| svipshipincom国产片| 国产主播在线观看一区二区| 欧美国产日韩亚洲一区| 色吧在线观看| 国产成人av激情在线播放| 99精品欧美一区二区三区四区| 男女那种视频在线观看| 中文字幕人妻丝袜一区二区| 免费人成在线观看视频色| 91麻豆精品激情在线观看国产| 日韩高清综合在线| 成人无遮挡网站| 亚洲在线自拍视频| 婷婷精品国产亚洲av在线| 我要搜黄色片| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 亚洲精品影视一区二区三区av| 欧美最黄视频在线播放免费| 九九久久精品国产亚洲av麻豆| 蜜桃久久精品国产亚洲av| 午夜福利18| 男人和女人高潮做爰伦理| 人妻丰满熟妇av一区二区三区| 国产综合懂色| 久久久色成人| 国产精品一区二区三区四区免费观看 | 精品国产亚洲在线| 99国产综合亚洲精品| 午夜精品一区二区三区免费看| 国产成人系列免费观看| 亚洲五月婷婷丁香| 特级一级黄色大片| 天堂动漫精品| 内地一区二区视频在线| 久久久久久久久大av| 观看免费一级毛片| 亚洲精品粉嫩美女一区| 在线观看舔阴道视频| 观看美女的网站| 内射极品少妇av片p| 欧美成人性av电影在线观看| 一级黄色大片毛片| 1024手机看黄色片| 神马国产精品三级电影在线观看| 黄色片一级片一级黄色片| 国产三级中文精品| 蜜桃亚洲精品一区二区三区| 九九久久精品国产亚洲av麻豆| 国产在视频线在精品| 美女高潮喷水抽搐中文字幕| 看免费av毛片| 精品久久久久久久毛片微露脸| 欧美乱码精品一区二区三区| 欧美激情在线99| av福利片在线观看| 1024手机看黄色片| 婷婷精品国产亚洲av在线| 久久久国产成人精品二区| 3wmmmm亚洲av在线观看| 制服人妻中文乱码| 亚洲熟妇熟女久久| 好看av亚洲va欧美ⅴa在| 69人妻影院| 婷婷精品国产亚洲av| 十八禁人妻一区二区| 欧美日韩黄片免| 一进一出抽搐动态| 精品久久久久久,| 亚洲avbb在线观看| 欧美黑人欧美精品刺激| 一级毛片高清免费大全| 在线观看美女被高潮喷水网站 | 国产午夜福利久久久久久| 一区二区三区免费毛片| 日韩国内少妇激情av| 国产高清激情床上av| 国产亚洲精品综合一区在线观看| eeuss影院久久| 久久久色成人| 又粗又爽又猛毛片免费看| 亚洲成a人片在线一区二区| 亚洲成人中文字幕在线播放| 级片在线观看| 国产三级黄色录像| 天天一区二区日本电影三级| 性色avwww在线观看| 两个人的视频大全免费| 亚洲欧美激情综合另类| www日本黄色视频网| 丰满人妻熟妇乱又伦精品不卡| 五月玫瑰六月丁香| 18+在线观看网站| 亚洲七黄色美女视频| 又紧又爽又黄一区二区| 久久久国产成人免费| 亚洲精品一区av在线观看| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 国产精品电影一区二区三区| 日韩欧美免费精品| 麻豆成人av在线观看| 国产精华一区二区三区| 国产真实乱freesex| 午夜福利高清视频| 欧美色欧美亚洲另类二区| av国产免费在线观看| 国产精品永久免费网站| 国产av不卡久久| 亚洲av中文字字幕乱码综合| 亚洲精品色激情综合| 国产高清三级在线| 九九久久精品国产亚洲av麻豆| 国产97色在线日韩免费| 免费在线观看影片大全网站| 成人午夜高清在线视频| 午夜免费男女啪啪视频观看 | 亚洲狠狠婷婷综合久久图片| 成人永久免费在线观看视频| 18美女黄网站色大片免费观看| 国产免费一级a男人的天堂| 国产美女午夜福利| 国产一区二区在线观看日韩 | 国产一区二区在线观看日韩 | 久久婷婷人人爽人人干人人爱| 夜夜躁狠狠躁天天躁| 乱人视频在线观看| av福利片在线观看| 18禁黄网站禁片午夜丰满| 丁香欧美五月| 久久久久久久亚洲中文字幕 | 久久久国产成人精品二区| 日韩欧美一区二区三区在线观看| 日日夜夜操网爽| 日本免费a在线| 久久国产精品影院| 男人的好看免费观看在线视频| 久久精品91蜜桃| 亚洲五月婷婷丁香| 欧美日韩黄片免| 俺也久久电影网| 小说图片视频综合网站| 叶爱在线成人免费视频播放| 看黄色毛片网站| 有码 亚洲区| 99热精品在线国产| 国产黄a三级三级三级人| 一本综合久久免费| 国产精品女同一区二区软件 | 国产亚洲精品一区二区www| 长腿黑丝高跟| 一本精品99久久精品77| www.色视频.com| 精品久久久久久久久久久久久| 国产伦精品一区二区三区视频9 | 久久精品综合一区二区三区| 99久久成人亚洲精品观看| 一a级毛片在线观看| 最新在线观看一区二区三区| 日韩欧美 国产精品| 1024手机看黄色片| 精品一区二区三区视频在线观看免费| 久久精品91无色码中文字幕| 精品一区二区三区视频在线观看免费| 欧美黄色淫秽网站| 美女高潮喷水抽搐中文字幕| 熟妇人妻久久中文字幕3abv| 成年人黄色毛片网站| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美 | 国产免费男女视频| 午夜精品在线福利| 日韩欧美精品v在线| 狂野欧美激情性xxxx| av国产免费在线观看| 色视频www国产| 嫩草影视91久久| 国产麻豆成人av免费视频| 精品一区二区三区视频在线观看免费| 美女cb高潮喷水在线观看| 精品一区二区三区视频在线观看免费| 亚洲五月婷婷丁香| 中文亚洲av片在线观看爽| 国产麻豆成人av免费视频| 草草在线视频免费看| 日日摸夜夜添夜夜添小说| 国产精品久久久久久人妻精品电影| 91久久精品电影网| 亚洲国产中文字幕在线视频| 琪琪午夜伦伦电影理论片6080| 亚洲av二区三区四区| 在线免费观看不下载黄p国产 | 欧美+亚洲+日韩+国产| 日本在线视频免费播放| 噜噜噜噜噜久久久久久91| 精品不卡国产一区二区三区| 日本黄色片子视频| 亚洲18禁久久av| 欧美精品啪啪一区二区三区| 欧美日韩福利视频一区二区| 国产伦精品一区二区三区视频9 | 有码 亚洲区| 国产黄片美女视频| 精品无人区乱码1区二区| 久久精品国产清高在天天线| 国产成+人综合+亚洲专区| 久久九九热精品免费| 少妇人妻一区二区三区视频| 欧美最黄视频在线播放免费| 母亲3免费完整高清在线观看| 此物有八面人人有两片| 国产又黄又爽又无遮挡在线| 国产精品女同一区二区软件 | 国产老妇女一区| 国产69精品久久久久777片| 波多野结衣高清作品| 国产精品一区二区三区四区免费观看 | 国产探花极品一区二区| 久久久国产成人免费| 亚洲美女视频黄频| 中文资源天堂在线| 久久性视频一级片| 久久精品国产自在天天线| 小蜜桃在线观看免费完整版高清| a在线观看视频网站| 欧美av亚洲av综合av国产av| 窝窝影院91人妻| 中国美女看黄片| 在线视频色国产色| 丁香六月欧美| 性色avwww在线观看| 少妇人妻一区二区三区视频| 国产精品野战在线观看| 久久久久国内视频| 搡老岳熟女国产| 午夜视频国产福利| 男人的好看免费观看在线视频| 美女免费视频网站| 波多野结衣巨乳人妻| 嫁个100分男人电影在线观看| 91av网一区二区| 午夜视频国产福利| 禁无遮挡网站| 午夜a级毛片| 国产午夜精品论理片| 午夜精品久久久久久毛片777| 草草在线视频免费看| 国产伦一二天堂av在线观看| 欧美日韩一级在线毛片| 久久久成人免费电影| 老司机深夜福利视频在线观看| av女优亚洲男人天堂| 久久99热这里只有精品18| 欧美绝顶高潮抽搐喷水| 18禁黄网站禁片午夜丰满| 在线免费观看的www视频| 少妇高潮的动态图| 桃色一区二区三区在线观看| 九色成人免费人妻av| 国产精品亚洲一级av第二区| 成年免费大片在线观看|