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

    河北省森林草原區(qū)草本植物物種多樣性和功能多樣性

    2017-07-21 09:21:37葛兆軒孫國龍袁業(yè)黃選瑞張志東
    草業(yè)學(xué)報(bào) 2017年7期
    關(guān)鍵詞:灌木林次生林草地

    葛兆軒,孫國龍,袁業(yè),黃選瑞,張志東*

    (1.河北省林木種質(zhì)資源與森林保護(hù)重點(diǎn)實(shí)驗(yàn)室,河北農(nóng)業(yè)大學(xué)林學(xué)院,河北 保定071000; 2.河北省塞罕壩機(jī)械林場, 河北 圍場068466)

    ?

    河北省森林草原區(qū)草本植物物種多樣性和功能多樣性

    葛兆軒1,孫國龍2,袁業(yè)1,黃選瑞2,張志東1*

    (1.河北省林木種質(zhì)資源與森林保護(hù)重點(diǎn)實(shí)驗(yàn)室,河北農(nóng)業(yè)大學(xué)林學(xué)院,河北 保定071000; 2.河北省塞罕壩機(jī)械林場, 河北 圍場068466)

    理解物種多樣性和功能多樣性之間的關(guān)系以及分析它們對環(huán)境變化的響應(yīng),有利于揭示生物多樣性的影響因素。本研究以塞罕壩地區(qū)不同植被類型(草地、灌木林和天然次生林)的草本群落為研究對象,分析了物種多樣性(物種豐富度、Shannon-Wiener多樣性和Pielou均勻度)和功能多樣性測度指標(biāo)(功能豐富度、功能均勻度和功能離散度)之間的相互關(guān)系,并采用相似性分析檢驗(yàn)(ANOSIM)和冗余度(RDA)分析等方法研究了它們隨植被類型和環(huán)境因子的變化規(guī)律。結(jié)果表明,草地中草本植物物種豐富度、Shannon-Wiener多樣性指數(shù)顯著高于灌木林和天然次生林中的草本群落(P<0.05),而功能豐富度指數(shù)在草地中最低(P<0.05),功能均勻度和功能離散度在灌木林中最低(P<0.05),Pielou均勻度指數(shù)在不同植被類型間差異不顯著(P>0.05);功能豐富度與物種豐富度、Shannon-Wiener多樣性指數(shù)顯著負(fù)相關(guān)(P<0.01),功能離散度與物種均勻度顯著正相關(guān)(P<0.01),而功能均勻度與所有物種多樣性指數(shù)均不相關(guān)(P>0.05);影響物種多樣性的主要因子有土壤含水量、土壤有機(jī)質(zhì)、土壤pH值、土壤全氮含量和坡向,影響功能多樣性的主要因子有土壤有機(jī)質(zhì)、土層深度、土壤含水量和海拔。物種多樣性和功能多樣性分布主要受土壤因子的影響,分別占到總解釋變異的29.1%和29.8%。本研究結(jié)果表明研究區(qū)草本植物物種多樣性和功能多樣性變化相對獨(dú)立,物種多樣性可能主要是由資源可利用程度決定,而功能多樣性可能是生態(tài)位分化即環(huán)境篩選和物種間相互作用共同決定。

    物種多樣性;功能多樣性;環(huán)境因子;植被類型;因子分解

    物種多樣性是生態(tài)系統(tǒng)功能的基礎(chǔ),它可以體現(xiàn)群落的結(jié)構(gòu)、群落組織水平、發(fā)展階段、穩(wěn)定性程度和生境差異[1]。功能多樣性是生物多樣性的一個(gè)非常重要的組成部分,是指影響群落或生態(tài)系統(tǒng)中有機(jī)體性狀的數(shù)值和范圍[2]。闡明物種多樣性和功能多樣性的關(guān)系是揭示生物多樣性對生態(tài)系統(tǒng)功能影響的關(guān)鍵[3]。環(huán)境因子作為一個(gè)“篩”,能夠決定哪些物種或者性狀可以在群落中生存和維持[4]。目前,不同植被類型物種多樣性和功能多樣性與環(huán)境因子間關(guān)系的研究仍然較少[5-7]。

    以往有關(guān)生物多樣性對生態(tài)系統(tǒng)功能影響的研究主要以物種多樣性來進(jìn)行,而忽視了生物多樣性的不同組分在生態(tài)系統(tǒng)中的功能效應(yīng)[8]。近年來,功能多樣性的研究逐漸受到重視[9-11]。一般來講,群落尺度物種多樣性和功能多樣性之間的關(guān)系表現(xiàn)為正線性相關(guān)[12]。受擾動(dòng)強(qiáng)度和環(huán)境脅迫的影響,二者之間的關(guān)系也會(huì)表現(xiàn)為負(fù)相關(guān)[13]、S型曲線關(guān)系[14]和不相關(guān)[15]。很多環(huán)境因素如海拔、坡向、土壤含水量、土壤養(yǎng)分等會(huì)影響物種或功能多樣性的分布。在亞高山草甸,北坡的物種多樣性顯著高于南坡,而功能多樣性卻呈相反的變化[16];隨著土壤含水量的增加,物種多樣性亦表現(xiàn)為增加趨勢,但功能多樣性變化不大[13]。由此可見,不同的生態(tài)系統(tǒng)物種多樣性和功能多樣性對不同環(huán)境梯度的響應(yīng)有可能不一致。在特定的區(qū)域,闡明物種多樣性和功能多樣性隨環(huán)境梯度的變化規(guī)律有助于更好地揭示物種共存機(jī)制和群落構(gòu)建規(guī)則[17]。

    塞罕壩地區(qū)是華北典型的森林-草原生態(tài)交錯(cuò)帶,生態(tài)系統(tǒng)復(fù)雜多樣,生物多樣性豐富,具有很好的保護(hù)價(jià)值。但由于長期人為干擾以及土地沙漠化、荒漠化的影響,生物多樣性呈下降趨勢。本研究以塞罕壩地區(qū)不同植被類型的草本群落為研究對象,通過樣地調(diào)查,主要探究以下3個(gè)問題:1) 物種多樣性與功能多樣性在不同植被類型間的變化是否相似?2) 物種多樣性和功能多樣性測度指標(biāo)之間關(guān)系如何?3) 影響物種多樣性和功能多樣性分布的主要環(huán)境因子是否相同?

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于河北省塞罕壩機(jī)械林場,地處內(nèi)蒙古高原與冀北山地交接處(42°22′ N,116°53′ E)。地形地貌組合為高原-波狀丘陵-漫灘-接壩山地[18]。氣候?qū)俸疁貛О敫珊?半濕潤季風(fēng)氣候。年均氣溫-1.40 ℃,極端最高氣溫32 ℃,極端最低氣溫-42.8 ℃;年平均降水量445 mm,蒸發(fā)量1358 mm[19]。植被由森林向草原過渡,天然次生林主要樹種有華北落葉松(Larixprincipis-rupprechtii)、山楊(Populusdavidiana)、蒙古櫟(Quercusmongolica)、白樺(Betulaplatyphylla);灌木種主要有山刺玫(Rosadavurica)、山杏(Armeniacasibirica)、柔毛繡線菊(Spiraeapubescens)等;草本主要由地榆(Sanguisorbaofficinalis)、金蓮花(Trolliuschinensis)、披針葉苔草(Carexlanceolata)和蚊子草(Filipendulapalmata)等物種構(gòu)成。

    1.2 樣地設(shè)置和群落調(diào)查

    2014年7—9月在研究區(qū)按照分層隨機(jī)取樣的方法進(jìn)行群落學(xué)調(diào)查,選取了32個(gè)斑塊(斑塊面積介于2.18~74.06 hm2)的調(diào)查數(shù)據(jù)進(jìn)行物種多樣性和功能多樣性分析,其中包括10個(gè)草本,10個(gè)灌木林和12個(gè)天然次生林斑塊。按照斑塊大小設(shè)置不同數(shù)量的草本樣方,樣方間隔≥50 m,樣方大小為1 m×1 m,共153個(gè)樣方。記錄每個(gè)樣方的海拔、坡度、坡向、坐標(biāo)以及植被類型;記錄草本植物物種名、蓋度、平均高度和最大高度。在每個(gè)草本樣方內(nèi),選取蓋度大于90%的所有物種,每個(gè)物種采集長勢良好的10株個(gè)體(樣方內(nèi)不足10株則于樣方附近采樣),每株取不少于5片完好成熟的葉片,測定葉片鮮重、葉面積、葉綠素含量,所有樣品置于105 ℃烘至恒重,以1/10000電子天平稱重,計(jì)算比葉面積、葉干物質(zhì)含量,測定方法詳見功能性狀測量手冊[20]。在采集葉片樣品的同時(shí),從每個(gè)植株莖基部位截取長約10 cm的枝條,每個(gè)物種取5~10根,放入塑封袋中置于保溫箱中保存。用水替代法測定枝條體積,然后置于105 ℃烘箱內(nèi)烘至恒重,稱量其干重,枝條干重與體積之比即為莖密度。將植物樣品用研磨儀粉碎,采用凱氏定氮法測定葉片氮含量,鉬銻抗比色法測定葉片磷含量,重鉻酸鉀法測定葉片有機(jī)碳含量[21]。

    在每個(gè)樣方中心位置采集0~10 cm、10~20 cm、20~30 cm三層土壤樣品,每層取一定量土樣混合后放到鋁盒中保存。將取得的土樣帶回實(shí)驗(yàn)室自然陰干后,參照《土壤農(nóng)化分析手冊》[22]測定土壤含水量、土壤pH值、有機(jī)質(zhì)含量、全磷含量、全氮含量和土壤速效鉀含量。

    1.3 生物多樣性計(jì)算

    物種重要值體現(xiàn)了物種在群落中的優(yōu)勢程度,草本物種的重要值用相對頻度、相對蓋度和相對高度的加和平均來表示。

    物種多樣性具有多種測度指標(biāo),本研究選取了3個(gè)常見的物種多樣性指數(shù),包括:物種豐富度、Shannon-Wiener指數(shù)和Pielou均勻度指數(shù)(表1)。

    表1 生物多樣性計(jì)算公式Table 1 The formula of biological diversity indices

    選取了118個(gè)草本物種的功能性狀(比葉面積、葉干物質(zhì)含量、葉綠素含量、葉氮含量、葉磷含量、葉碳含量和莖密度)實(shí)測數(shù)據(jù),進(jìn)行功能多樣性的計(jì)算。功能豐富度指群落中物種占有的n維功能性狀凹凸包量;功能均勻度指群落中功能性狀值在凹凸包量中排列的均勻性;功能離散度指群落中每個(gè)物種的n維功能性狀值到所有功能性狀空間重心的平均距離[23-25]。計(jì)算公式見表1。

    環(huán)境數(shù)據(jù)共選取了13個(gè)因子,其中包括2個(gè)空間因子、4個(gè)地形因子和7個(gè)土壤因子(表2)。

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

    物種累積曲線(species accumulation curve)用于比較不同植被類型物種累積速率的變化。采用ANOSIM(analysis of similarities)方法檢驗(yàn)物種組成在3個(gè)不同植被類型之間是否有顯著差異。對物種多樣性和功能多樣性指數(shù)分別進(jìn)行單因素方差分析(one-way ANOVA),并進(jìn)行了不同植被類型物種多樣性和功能多樣性的多重比較(Tukey HSD)。應(yīng)用皮爾遜相關(guān)方法分析物種多樣性和功能多樣性之間的關(guān)系。

    表 2 環(huán)境因子的選擇與定義Table 2 Selection and definition of environmental variables

    運(yùn)用冗余度分析(redundancy analysis,RDA)探討環(huán)境因子對不同植被類型物種多樣性和功能多樣性的影響。首先建立物種多樣性-環(huán)境變量、功能多樣性-環(huán)境變量矩陣,用去趨勢對應(yīng)分析(detrended correspondence analysis,DCA)估計(jì)排序軸梯度長度分別為0.1和0.7,表明兩個(gè)數(shù)據(jù)文件均具有較好的線性反應(yīng),利用線性模型比較適宜。為了檢驗(yàn)環(huán)境因子對物種和功能多樣性的影響,首先利用Monte Carlo檢驗(yàn)(置換次數(shù)為499)判斷其重要性是否顯著,當(dāng)候選變量P≥0.05時(shí),予以排除。為了更好地反映環(huán)境因子對多樣性影響的解釋效果,將解釋變異分解為以下幾部分:1)純空間因子解釋部分(Rsp);2)純地形因子解釋部分(Rto);3)純土壤因子解釋部分(Rso);4)空間因子和地形因子交互解釋部分(Rspto);5)空間因子和土壤因子交互解釋部分(Rspso);6)土壤因子和地形因子交互解釋部分(Rsoto);7)3個(gè)因子交互解釋部分(Rsptoso)。

    利用R 3.3.1程序FD軟件包中的dbFD函數(shù)計(jì)算物種多樣性及功能多樣性指數(shù),統(tǒng)計(jì)分析在R 3.3.1程序Vegan軟件包中完成(R development Core Team, 2011)。

    2 結(jié)果與分析

    2.1 物種組成

    圖1 基于樣本的物種累積曲線Fig.1 Sample-based species accumulation curve

    在32個(gè)斑塊的3個(gè)植被類型中共記錄了294個(gè)草本物種,分屬46科,168屬。其中草地物種數(shù)最高,共記錄了187個(gè)物種,分屬38科、118屬,分別占了總物種、科和屬數(shù)的63.6%、82.6%和70.2%;其次是灌木林,共記錄了156個(gè)物種,分屬34科、106屬,分別占了總物種、科和屬數(shù)的53.1%、73.9%和63.1%;天然次生林物種數(shù)最低,共記錄了121個(gè)物種,分屬27科、81屬,分別占了總物種、科和屬數(shù)的40.8%、58.6%和48.2%。

    3個(gè)植被類型中草本物種數(shù)隨樣地?cái)?shù)累積的結(jié)果如圖1,草地植被類型具有較快的物種累積速率,天然次生林物種累積速率較慢。并且ANOSIM檢驗(yàn)結(jié)果顯示,3個(gè)植被類型物種組成差異顯著(r=0.4167,P=0.001)。

    2.2 不同植被類型草本植物物種多樣性和功能多樣性變化

    不同植被類型間物種豐富度、Shannon-Wiener多樣性指數(shù)(H)差異顯著(P<0.05),草地物種豐富度、Shannon-Wiener多樣性指數(shù)(H)顯著高于灌木林和次生林植被類型。草地和次生林植被類型間功能豐富度(FRic)差異顯著(P<0.05),草地功能豐富度(FRic)顯著低于次生林和灌木林植被類型,次生林和灌木林植被類型間功能均勻度(FEve)和功能離散度(FDis)差異顯著(P<0.05),并且灌木林低于次生林和草地植被類型(圖2)。

    圖2 不同植被類型物種多樣性和功能多樣性差異性分析Fig.2 Difference analysis of species diversity and functional diversity among different vegetation types 不同字母表示顯著差異(P<0.05) The different letters indicate significant differences at P<0.05 level; GL:草地 Grassland; SF: 灌木林 Shrub forest; NS: 天然次生林 Natural secondary forest.下同 The same below.

    2.3 物種多樣性與功能多樣性的關(guān)系

    FRic與S、H極顯著負(fù)相關(guān)(P<0.01),而與J關(guān)系不大(P>0.05);FDis與J呈極顯著正相關(guān)(P<0.01),而與S和H不相關(guān)(P>0.05);FEve與所有物種多樣性指數(shù)均不相關(guān)(P>0.05)(表3)。

    表3 物種多樣性與功能多樣性相關(guān)性Table 3 The relationship between species diversity and functional diversity

    “**”,P<0.01 ;“-”,負(fù)相關(guān) Negative correlation.

    2.4 不同植被類型環(huán)境因子對比分析

    通過對不同植被類型的環(huán)境因子進(jìn)行單因素方差分析,土壤全磷含量、全氮含量和土壤含水量差異顯著(P<0.05),只對這3個(gè)因子進(jìn)行了多重比較(圖3)。草地土壤全氮含量、土壤含水量顯著高于灌木林和天然次生林,草地中全磷含量顯著高于灌木林,而與天然次生林差異不顯著(P>0.05)。

    圖3 不同植被類型環(huán)境因子對比分析Fig.3 Variations of environmental variables among different vegetation types

    2.5 物種多樣性和功能多樣性與環(huán)境因子相關(guān)性

    土壤含水量、土壤pH值、坡向和土壤全氮含量與物種多樣性各指數(shù)正相關(guān),而隨著土層深度和速效鉀含量的增加,物種豐富度、Shannon-Wiener多樣性降低(圖4A);對于功能多樣性來說,F(xiàn)Ric主要受到土壤有機(jī)質(zhì)和土層深度的影響而隨土壤含水量的增加而減小,F(xiàn)Eve與土層深度正相關(guān)性,F(xiàn)Dis主要受海拔、土壤含水量和土層深度的影響,表現(xiàn)為正相關(guān)(圖4B)。

    圖4 物種多樣性(A)和功能多樣性(B)與環(huán)境因子RDA排序圖Fig.4 Redundancy analysis ordination diagram of the species diversity (A), functional diversity (B) and environmental factors

    2.6 因子分解

    對影響物種多樣性的環(huán)境因子進(jìn)行定量分解可知,空間因子、地形因子和土壤因子共同解釋效果達(dá)到了53.3%,其中純土壤因子的解釋效果最顯著,占到了總解釋變異部分的29.1%,而純空間因子和純地形因子的解釋效果相對較弱,分別只占到9.3%和10.1%,各環(huán)境因子之間相互耦合作用相對較小(圖5A)。

    對于功能多樣性(圖5B),3類環(huán)境因子共同解釋效果達(dá)到60.5%,純空間因子所起到的解釋效果相對較弱,占總變異解釋部分的6.5%,純地形因子所起到的解釋效果占16.3%,土壤因子的解釋能力為29.8%。空間和地形因子的耦合作用相對較大,占總變異解釋部分的5.0%,3類環(huán)境因子共同作用最小,占總變異解釋部分的0.2%。

    圖5 影響物種多樣性(A)和功能多樣性(B)的環(huán)境因子的分解變異Fig.5 Variance partitioning of environmental factors for explaining species diversity (A) and functional diversity (B)

    3 討論

    3.1 物種多樣性變化

    一般認(rèn)為,物種多樣性是由資源可利用程度決定的[26],也有的學(xué)者將這一觀點(diǎn)稱之為生產(chǎn)力假說或能量假說[27],即高的環(huán)境能量形成高的生產(chǎn)力,進(jìn)而承載高的物種多樣性。如隨著土壤含水量[16]、土壤氮含量[28]的增加,物種多樣性呈增加趨勢,而隨生境干燥程度的增加而降低[13]。本研究中,草地具有較高的土壤含水量、土壤全氮含量和土壤全磷含量(圖3)。相對于灌木林和天然次生林,草地中可以說具有良好的資源供給能力,因而導(dǎo)致了更多的物種共存。這可從草地中草本物種豐富度和Shannon-Wiener多樣性指數(shù)較高,灌木林次之,而天然次生林最低(圖2)的結(jié)果中得到很好的證明。通過冗余度分析表明物種多樣性各指數(shù)與土壤含水量、土壤全氮含量表現(xiàn)出較強(qiáng)正相關(guān)性(圖4A),表明研究區(qū)天然次生林和灌木林中相對較低的土壤含水量和全氮含量是影響草本植物生長和分布的限制性因素,這在其他學(xué)者的研究中[28-29]也得到了相似的結(jié)論。在塞罕壩地區(qū),內(nèi)陸河岔較多,多數(shù)水源補(bǔ)充入就近的沼澤或?yàn)┑亍?jù)塞罕壩機(jī)械林場資料記載,該區(qū)有沼澤地和灘地面積約1326.7 hm2。在本研究調(diào)查的10個(gè)草地斑塊中,大多數(shù)落在了濕地中,導(dǎo)致了草地中具有較高的土壤水分和養(yǎng)分含量,草本物種多樣性較高。

    通過對所選環(huán)境因子進(jìn)行定量分解可知,影響物種多樣性分布的未能解釋部分占到了46.7% (圖5A),表明該地區(qū)草本物種多樣性除了受土壤水分、養(yǎng)分等的影響外,還可能受種間競爭(喬草、灌草間)、人為干擾(商業(yè)采伐、放牧等)、光環(huán)境異質(zhì)性等的影響[30],更多的機(jī)理性研究仍有待進(jìn)一步探討。

    3.2 功能多樣性變化

    在一個(gè)群落中,物種功能性狀的值、范圍和分布程度決定了功能多樣性的大小[2]。環(huán)境篩選會(huì)使適應(yīng)特定生境的物種共存,物種的功能性狀出現(xiàn)趨同現(xiàn)象,這在一定程度上會(huì)降低功能多樣性[8,14]。但是,在環(huán)境脅迫的條件下,性狀生態(tài)位空間也可能會(huì)導(dǎo)致強(qiáng)烈的分化,這樣具有不同功能性狀的物種可以高效互補(bǔ)的利用有限的生態(tài)資源[13]。與物種豐富度相反,本研究中,功能豐富度在草地中最低,而在天然次生林和灌木林中相對較高(圖2)。這一結(jié)果很可能是因?yàn)?,在塞罕壩機(jī)械林場,由于過去長期的以追求木材效益最大化的經(jīng)營目標(biāo),大面積的立地條件較好的天然次生林和灌木林皆被伐用于種植用材林,現(xiàn)存的天然次生林多處于坡度較陡、土壤瘠薄的生境中。由于土壤水分、養(yǎng)分等資源的限制,導(dǎo)致物種生態(tài)位分化程度要顯著高于生境條件較好的草地中,進(jìn)而引起了性狀趨異,功能多樣性增加。Bello等[13]在研究功能多樣性在干旱和濕潤地區(qū)的差異上,也發(fā)現(xiàn)了類似的規(guī)律。在本研究中,除了土壤因子、地形因子和空間因子外,其他未能解釋的部分仍然高達(dá)39.5%(圖5B)。有研究表明植物群落內(nèi)光環(huán)境異質(zhì)性的變化、種間和種內(nèi)作用、外部干擾等均會(huì)影響植物功能性狀的變化和個(gè)體生態(tài)策略的改變,進(jìn)而影響功能多樣性[13,31]。

    3.3 物種多樣性和功能多樣性

    本研究表明功能多樣性的測度指標(biāo)在不同植被類型間的變化趨勢與物種多樣性的變化模式并不相同(圖1),不同測度指標(biāo)之間存在正相關(guān)、負(fù)相關(guān)和不相關(guān)的情況(表3),證實(shí)了物種間的功能差異并不能由物種多樣性簡單代替,與其并不成比例關(guān)系[32]。相似的結(jié)論也在其他的研究中得到了驗(yàn)證。如在森林和采伐跡地生境,由于生態(tài)構(gòu)建過程的不同,物種多樣性和功能多樣性間的關(guān)系存在顯著差異[33]。也有研究指出物種多樣性與功能多樣性呈正相關(guān)[8,12],隨著干擾和環(huán)境脅迫的增加,物種多樣性與功能多樣性也可由正相關(guān)變?yōu)椴幌嚓P(guān)[34]或負(fù)相關(guān)[13,16,35]。由前文分析可知,物種多樣性主要是由資源可利用程度決定,而功能多樣性是由環(huán)境篩選和物種間相互作用即生態(tài)位分化所決定,不同的生態(tài)構(gòu)建過程影響了不同生境中物種多樣性和功能多樣性的關(guān)系差異[33]。本研究物種多樣性相關(guān)指標(biāo)沿草地、灌木林和次生林植被類型顯著降低,而功能多樣性測度指標(biāo)在草地中最低,而天然次生林和灌木林相對較高的結(jié)論,正是以上生態(tài)學(xué)假說綜合作用的結(jié)果。

    4 結(jié)論

    本研究揭示了草地、灌木林和天然次生林下草本群落物種多樣性和功能多樣性測度指標(biāo)隨環(huán)境梯度(如土壤含水量、土壤有機(jī)質(zhì)等)的響應(yīng)差異。物種豐富度、Shannon-Wiener多樣性指數(shù)沿草地、灌木林和次生林植被類型呈顯著降低趨勢,而功能豐富度則在草地植被類型中最低,天然次生林和灌木林相對較高。相關(guān)分析表明,功能多樣性和物種多樣性測度指標(biāo)之間存在正相關(guān)、負(fù)相關(guān)和不相關(guān)的情況。以上結(jié)果反映了不同植被類型草本群落構(gòu)建機(jī)制的不同。在研究區(qū),土壤養(yǎng)分差異是影響物種多樣性的主要原因,土壤水分、養(yǎng)分脅迫導(dǎo)致生態(tài)位分化,決定功能多樣性。本研究結(jié)果進(jìn)一步證實(shí)了物種多樣性主要是由資源可利用程度決定,而功能多樣性是由環(huán)境篩選和物種間相互作用即生態(tài)位分化共同決定,該研究結(jié)論可為塞罕壩地區(qū)生物多樣性保護(hù)和管理提供科學(xué)依據(jù)。

    References:

    [1] Ma K P, Huang J H, Yu S L,etal. Plant community diversity in Dongling Mountain, Beijing, China; Ⅱ. Species richness, evenness and species diversities. Acta Ecologica Sinica, 1995, 15(3): 268-277. 馬克平, 黃建輝, 于順利, 等. 北京東靈山地區(qū)植物群落多樣性的研究Ⅱ豐富度、均勻度和物種多樣性指數(shù). 生態(tài)學(xué)報(bào), 1995, 15(3): 268-277.

    [2] Petchey O L, Gaston K J. Functional diversity: back to basics and looking forward. Ecology Letters, 2006, 9(6): 741-758.

    [3] Zhang Q G, Zhang D Y. Biodiversity and ecosystem functioning: recent advances and trends. Biodiversity Science, 2003, 11(5): 351-363. 張全國, 張大勇. 生物多樣性與生態(tài)系統(tǒng)功能:最新的進(jìn)展與動(dòng)向. 生物多樣性, 2003, 11(5): 351-363.

    [4] Bu W S. The Relationships Between Biodiversity and Ecosystem Functioning in Natural Tropical Forests of Hainan Island, China[D]. Beijing: Chinese Academy of Forestry, 2013. 卜文圣. 海南島熱帶天然林生物多樣性與生態(tài)系統(tǒng)功能關(guān)系的研究[D]. 北京: 中國林業(yè)科學(xué)研究院, 2013.

    [5] Suding K N, Goldstein L J. Testing the Holy Grail framework: using functional traits to predict ecosystem change. New Phytologist, 2008, 180(3): 559-562.

    [6] Ledo A, Burslem D F R P, Condés S,etal. Micro-scale habitat associations of woody plants in a neotropical cloud forest. Journal of Vegetation Science, 2013, 24(6): 1086-1097.

    [7] Jiang Y, Zang R G, Letcher S G,etal. Associations between plant composition/diversity and the abiotic environment across six vegetation types in a biodiversity hotspot of Hainan Island, China. Plant and Soil, 2016, 403(1): 21-35.

    [8] Díaz S, Cabido M. Vive la différence: plant functional diversity matters to ecosystem processes. Trends in Ecology & Evolution, 2001, 16(11): 646-655.

    [9] Stevens R D, Cox S B, Strauss R E,etal. Patterns of functional diversity across an extensive environmental gradient: vertebrate consumers, hidden treatments and latitudinal trends. Ecology Letters, 2003, 6(12): 1099-1108.

    [10] Sun G J, Zhang R, Zhou L. Trends and advances in researches on plant functional diversity and functional groups. Acta Ecologica Sinica, 2003, 23(7): 1430-1435. 孫國鈞, 張榮, 周立. 植物功能多樣性與功能群研究進(jìn)展. 生態(tài)學(xué)報(bào), 2003, 23(7): 1430-1435.

    [11] Jiang X L, Zhang W G. Functional diversity and its research method. Acta Ecologica Sinica, 2010, 30(10): 2766-2773. 江小雷, 張衛(wèi)國. 功能多樣性及其研究方法. 生態(tài)學(xué)報(bào), 2010, 30(10): 2766-2773.

    [12] Petchey O L, Gaston K J. Functional diversity (FD), species richness and community composition. Ecology Letters, 2002, 5(3): 402-411.

    [13] Bello F D, Jan L, Sebastià M T. Variations in species and functional plant diversity along climatic and grazing gradients. Ecography, 2006, 29(6): 801-810.

    [14] Sasaki T, Okubo S, Okayasu T,etal. Two-phase functional redundancy in plant communities along a grazing gradient in Mongolian rangelands. Ecology, 2009, 90(9): 2598-2608.

    [15] Li X G, Zhu Z H, Zhou X S,etal. Effects of clipping, fertilizing and watering on the relationship between species diversity, functional diversity and primary productivity in alpine meadow of China. Chinese Journal of Plant Ecology, 2011, 35(11): 1136-1147. 李曉剛, 朱志紅, 周曉松, 等. 刈割、施肥和澆水對高寒草甸物種多樣性、功能多樣性與初級生產(chǎn)力關(guān)系的影響. 植物生態(tài)學(xué)報(bào), 2011, 35(11): 1136-1147.

    [16] Li X E. Patterns of species diversity and functional diversity along a south- to north-facing slope in a subalpine meadow. Community Ecology, 2011, 12(2): 179-187.

    [17] Mason N W H, Bello F D, Dole?al J,etal. Niche overlap reveals the effects of competition, disturbance and contrasting assembly processes in experimental grassland communities. Journal of Ecology, 2011, 99(3): 788-796.

    [18] Qin P Y, Yang H J, Jiang F L,etal. Quantitative classification of natural plant communities in the Saihanba National Nature Reserve, Hebei, Province, China. Chinese Journal of Applie Ecology, 2016, 27(5): 1383-1392. 秦朋遙, 楊會(huì)娟, 蔣鳳玲, 等. 河北省塞罕壩保護(hù)區(qū)天然植物群落數(shù)量分類. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(5): 1383-1392.

    [19] Zhang S Z, Li M, Zhang S B,etal. Factors affecting natural regeneration ofLarixprincipis-rupprechtiiplantations in Saihanba of Hebei, China. Acta Ecologica Sinica, 2015, 35(16): 5403-5411. 張樹梓, 李梅, 張樹彬, 等. 塞罕壩華北落葉松人工林天然更新影響因子. 生態(tài)學(xué)報(bào), 2015, 35(16): 5403-5411.

    [20] Pérez-Harguindeguy N, Díaz S, Garnier E,etal. New handbook for standardised measurement of plant functional traits worldwide. Australian Journal of Botany, 2013, 61(3): 167-234.

    [21] Cornelissen J H C, Lavorel S, Garnier E,etal. A handbook of protocols for standardised and easy measurement of plant functional traits worldwide. Australian Journal of Botany, 2003, 51(4): 335-380.

    [22] Lu R K. Soil and Agricultural Chemistry Analysis[M]. Beijing: Science and Technology of China Agriculture Press, 2000. 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000.

    [23] Mason N W H, Mouillot D, Lee W G,etal. Functional richness, functional evenness and functional divergence: the primary components of functional diversity. Oikos, 2005, 111(1): 112-118.

    [24] Villéger S, Mason N W H, Mouillot D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology, 2008, 89(8): 2290-2301.

    [25] Laliberté E, Legendre P. A distance-based framework for measuring functional diversity from multiple traits. Ecology, 2010, 91(1): 299-305.

    [26] Morin P J. Biodiversity’s ups and downs. Nature, 2000, 406: 463-464.

    [27] Currie D J, Francis A P. Regional versus climatic effect on taxon richness in angiosperms: reply to Qian and Ricklefs. American Naturalist, 2004, 163(5): 780-785.

    [28] Perroni-Ventura Y, Montaa C, García-Oliva F. Relationship between soil nutrient availability and plant species richness in a tropical semi-arid environment. Journal of Vegetation Science, 2006, 17(6): 719-728.

    [29] Bai Y F, Li L H, Wang Q B,etal. Changes in plant species diversity and productivity along gradients of precipitation and elevation in the Xilin River Basin, Inner Mongolia. Chinese Journal of Plant Ecology, 2000, 24(6): 667-673. 白永飛, 李凌浩, 王其兵, 等. 錫林河流域羊草草原植物種群和功能群的長期變異性及其對群落穩(wěn)定性的影響. 植物生態(tài)學(xué)報(bào), 2000, 24(6): 667-673.

    [30] La Q, Zhaxi Ciren, Zhu W D,etal. Plant species-richness and association with environmental factors in the riparian zone of the Yarlung Zangbo River of Tibet, China. Biodiversity Science, 2014, 22(3): 337-347. 拉瓊, 扎西次仁, 朱衛(wèi)東, 等. 雅魯藏布江河岸植物物種豐富度分布格局及其環(huán)境解釋. 生物多樣性, 2014, 22(3): 337-347.

    [31] Markesteijn L, Poorter L, Bongers F. Light-dependent leaf trait variation in 43 tropical dry forest tree species. American Journal of Botany, 2007, 94(4): 515-525.

    [32] Devictor V, Mouillot D, Meynard C,etal. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. Ecology Letters, 2010, 13(8): 1030-1040.

    [33] Mayfield M M, Boni M F, Daily G C,etal. Species and functional diversity of native and human-dominated plant communities. Ecology, 2005, 86(9): 2365-2372.

    [34] Kong B B, Wei X H, Du J L,etal. Effects of clipping and fertilization on the temporal dynamics of species diversity and functional diversity and their relationships in an alpine meadow. Chinese Journal of Plant Ecology, 2016, 40(3): 187-199. 孔彬彬, 衛(wèi)欣華, 杜家麗, 等. 刈割和施肥對高寒草甸物種多樣性和功能多樣性時(shí)間動(dòng)態(tài)及其關(guān)系的影響. 植物生態(tài)學(xué)報(bào), 2016, 40(3): 187-199.

    [35] Zang Y M, Zhu Z H, Li Y N,etal. Effects of species diversity and functional diversity on primary productivity of alpine meadow. Chinese Journal of Ecology, 2009, 28(6): 999-1005. 臧岳銘, 朱志紅, 李英年, 等. 高寒矮嵩草草甸物種多樣性與功能多樣性對初級生產(chǎn)力的影響. 生態(tài)學(xué)雜志, 2009, 28(6): 999-1005.

    Herbaceous plant species diversity and functional diversity in the forest-steppe zone of Hebei, China

    GE Zhao-Xuan1, SUN Guo-Long2, YUAN Ye1, HUANG Xuan-Rui2, ZHANG Zhi-Dong1*

    1.HebeiProvinceKeyLaboratoryofForestTreesGermplasmResourcesandForestProtection,CollegeofForestry,HebeiAgriculturalUniversity,Baoding071000,China; 2.SaihanbaMechanizedForestryCentreofHebeiProvince,Weichang068466,China

    Understanding the relationships between environmental/soil factors and plant species diversity and functional diversity is useful for identifying which factors have the strongest effects on plant communities. In this study, we first studied the relationships between different components of herbaceous plant species diversity (species richness, Shannon-Wiener diversity, and Pielou evenness) and functional diversity indices (functional richness, functional evenness, and functional dispersion). Then, we used methods of analysis of similarities (ANOSIM) and redundancy analysis to explore their variation across three vegetation types (grassland, shrub, and natural secondary forest) and different environmental gradients. The results showed that species richness and Shannon-Wiener diversity were significantly higher in grassland than in shrub and natural secondary forest (P<0.05); however, the lowest functional richness was in grassland (P<0.05) and the lowest functional evenness and functional dispersion were in shrub forests (P<0.05). There was no significant difference in Pielou evenness among the three vegetation types (P>0.05). Functional richness was negatively related to species richness and Shannon-Wiener diversity (P<0.01). Functional dispersion was significantly positively correlated with species evenness (P<0.01), whereas functional evenness had a non-significant correlation with all species diversity indices (P>0.05). The main factors affecting species diversity were soil water content, soil organic matter, soil pH, soil total nitrogen content, and aspect. The main factors correlated with functional diversity were soil organic matter, soil depth, soil water content, and altitude. Consequently, soil factors were the vital determinants of plant species diversity and functional diversity in the study area, and accounted for 29.1% and 29.8% of the total variance explained by environmental factors, respectively. The results showed that herbaceous plant species diversity and functional diversity varied independently in the study area. Species diversity was mainly determined by resource availability, whereas functional diversity was mainly related to niche differentiation (the effects of both environmental filters and species interactions).

    species diversity; functional diversity; environmental factors; vegetation types; variation partitioning

    10.11686/cyxb2017015

    2017-01-18;改回日期:2017-03-15

    國家自然科學(xué)基金項(xiàng)目(31370636)和林業(yè)公益性行業(yè)科研專項(xiàng)(20150430304)資助。

    葛兆軒(1990-),女,河北唐山人,在讀碩士。E-mail: zhxuan_ge@126.com

    *通信作者Corresponding author. E-mail: zhzhido@163.com

    http://cyxb.lzu.edu.cn

    葛兆軒, 孫國龍, 袁業(yè), 黃選瑞, 張志東. 河北省森林草原區(qū)草本植物物種多樣性和功能多樣性. 草業(yè)學(xué)報(bào), 2017, 26(7): 35-44.

    GE Zhao-Xuan, SUN Guo-Long, YUAN Ye, HUANG Xuan-Rui, ZHANG Zhi-Dong. Herbaceous plant species diversity and functional diversity in the forest-steppe zone of Hebei, China. Acta Prataculturae Sinica, 2017, 26(7): 35-44.

    猜你喜歡
    灌木林次生林草地
    草地上的事
    幼兒100(2020年31期)2020-11-18 03:42:00
    內(nèi)蒙古自治區(qū)灌木林資源現(xiàn)狀與保護(hù)發(fā)展對策
    草地
    青海高原高山灌木林植被特點(diǎn)及主要類型
    淺論北方次生林培育技術(shù)要點(diǎn)
    次生林冠下紅松更新對光環(huán)境適應(yīng)的研究進(jìn)展
    阜新地區(qū)天然次生林食葉害蟲無公害防治應(yīng)用
    草地上
    北方次生林培育技術(shù)措施探討
    蘭州南北兩山雨養(yǎng)灌木林地植物群落穩(wěn)定性調(diào)查研究
    国产中年淑女户外野战色| 美女高潮的动态| 国产精品蜜桃在线观看 | 久久亚洲精品不卡| 变态另类丝袜制服| 国产黄色小视频在线观看| 久久精品人妻少妇| 成人午夜高清在线视频| 国产亚洲91精品色在线| 99久久中文字幕三级久久日本| 国产男人的电影天堂91| 黑人高潮一二区| 亚洲国产欧美在线一区| 丰满乱子伦码专区| 欧美性猛交黑人性爽| 免费看美女性在线毛片视频| 最近中文字幕高清免费大全6| 亚洲精品久久国产高清桃花| 99riav亚洲国产免费| 少妇人妻精品综合一区二区 | 成人永久免费在线观看视频| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 国产精品久久久久久精品电影小说 | 欧美极品一区二区三区四区| av福利片在线观看| 久久精品国产清高在天天线| 久久九九热精品免费| 别揉我奶头 嗯啊视频| www.av在线官网国产| 国产午夜精品久久久久久一区二区三区| 亚洲欧洲日产国产| 久久久久性生活片| 亚洲欧美日韩东京热| 午夜福利在线在线| 免费观看人在逋| 久久精品国产亚洲网站| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观| 日韩强制内射视频| 嘟嘟电影网在线观看| 久久久久久久久久成人| 日韩精品有码人妻一区| 高清毛片免费观看视频网站| 免费观看精品视频网站| 亚洲av中文av极速乱| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 美女被艹到高潮喷水动态| 亚洲成人久久性| 变态另类丝袜制服| 午夜爱爱视频在线播放| 日韩大尺度精品在线看网址| 中文字幕av成人在线电影| 色哟哟·www| 国产真实乱freesex| 一级毛片aaaaaa免费看小| 悠悠久久av| 99在线人妻在线中文字幕| 国产老妇伦熟女老妇高清| 日韩视频在线欧美| 日日啪夜夜撸| av在线播放精品| 人妻久久中文字幕网| 中国美白少妇内射xxxbb| 国产极品天堂在线| 亚洲欧美日韩卡通动漫| 好男人在线观看高清免费视频| 毛片一级片免费看久久久久| 不卡视频在线观看欧美| 熟妇人妻久久中文字幕3abv| 久久精品夜夜夜夜夜久久蜜豆| 色播亚洲综合网| 精品久久久久久成人av| 国产高清不卡午夜福利| 国产成人福利小说| 久久精品国产亚洲网站| h日本视频在线播放| 色综合亚洲欧美另类图片| 国产成人a∨麻豆精品| 久久精品国产亚洲av天美| 波野结衣二区三区在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热这里只有是精品50| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 非洲黑人性xxxx精品又粗又长| 人人妻人人看人人澡| av在线观看视频网站免费| 久久综合国产亚洲精品| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看 | 久久久国产成人免费| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 欧美一区二区国产精品久久精品| 三级毛片av免费| av在线亚洲专区| 亚洲精品成人久久久久久| 69av精品久久久久久| 99久久人妻综合| 在线观看66精品国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品.久久久| 男女做爰动态图高潮gif福利片| 亚洲高清免费不卡视频| 精品久久久噜噜| 成人毛片a级毛片在线播放| 色哟哟·www| 夜夜夜夜夜久久久久| 久久精品国产亚洲av涩爱 | 精品久久久噜噜| 国产乱人视频| 超碰av人人做人人爽久久| 国产精品蜜桃在线观看 | 精品免费久久久久久久清纯| 午夜久久久久精精品| 国内精品美女久久久久久| 精品久久国产蜜桃| 人妻久久中文字幕网| 在线a可以看的网站| 只有这里有精品99| 黄片无遮挡物在线观看| 欧洲精品卡2卡3卡4卡5卡区| АⅤ资源中文在线天堂| 舔av片在线| 亚洲婷婷狠狠爱综合网| 桃色一区二区三区在线观看| 真实男女啪啪啪动态图| 国产成人freesex在线| 国产爱豆传媒在线观看| 99热网站在线观看| 综合色丁香网| 精品一区二区三区视频在线| 国产乱人偷精品视频| 伦精品一区二区三区| 国产免费一级a男人的天堂| 99久久久亚洲精品蜜臀av| 久久久欧美国产精品| 欧美+日韩+精品| 国产高潮美女av| 欧美bdsm另类| 男女边吃奶边做爰视频| 激情 狠狠 欧美| 亚洲av免费高清在线观看| 人妻久久中文字幕网| 18+在线观看网站| 日本黄色视频三级网站网址| 亚洲精品日韩在线中文字幕 | 欧美xxxx性猛交bbbb| 久久精品国产鲁丝片午夜精品| 一个人看视频在线观看www免费| 日韩欧美 国产精品| 九九久久精品国产亚洲av麻豆| 少妇裸体淫交视频免费看高清| 国产黄色视频一区二区在线观看 | 99久久久亚洲精品蜜臀av| 免费看a级黄色片| 成人性生交大片免费视频hd| 国产精品久久久久久精品电影小说 | 在线免费观看不下载黄p国产| 蜜桃久久精品国产亚洲av| 亚洲精华国产精华液的使用体验 | 久久九九热精品免费| 亚洲美女搞黄在线观看| 久久久久久久久中文| 三级国产精品欧美在线观看| 精品一区二区免费观看| 两个人视频免费观看高清| 韩国av在线不卡| 欧美最黄视频在线播放免费| av福利片在线观看| 男女下面进入的视频免费午夜| 欧美性感艳星| 国产黄色视频一区二区在线观看 | 成人特级黄色片久久久久久久| 久久久久久久久久成人| 看免费成人av毛片| 此物有八面人人有两片| 高清日韩中文字幕在线| 欧美在线一区亚洲| 免费看日本二区| 99精品在免费线老司机午夜| 午夜福利在线观看免费完整高清在 | or卡值多少钱| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 国内精品一区二区在线观看| 国产精品无大码| 免费看美女性在线毛片视频| 国产日本99.免费观看| 免费人成在线观看视频色| 熟女电影av网| 亚洲精品国产成人久久av| 校园春色视频在线观看| 国产av在哪里看| videossex国产| 久久这里有精品视频免费| 国产高潮美女av| 色噜噜av男人的天堂激情| 高清在线视频一区二区三区 | 波多野结衣高清作品| 亚洲国产欧美人成| 偷拍熟女少妇极品色| 乱系列少妇在线播放| 十八禁国产超污无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 中文字幕熟女人妻在线| 少妇人妻精品综合一区二区 | 日韩精品青青久久久久久| 99在线人妻在线中文字幕| 亚洲av免费高清在线观看| 真实男女啪啪啪动态图| 你懂的网址亚洲精品在线观看 | 亚洲国产欧美人成| 国产91av在线免费观看| 国产精品福利在线免费观看| 精品一区二区三区视频在线| 在线天堂最新版资源| 美女xxoo啪啪120秒动态图| 国产大屁股一区二区在线视频| 日韩精品青青久久久久久| 久久午夜福利片| 日韩在线高清观看一区二区三区| 亚洲自拍偷在线| 日日摸夜夜添夜夜爱| 男人舔奶头视频| 能在线免费观看的黄片| 99在线视频只有这里精品首页| 免费黄网站久久成人精品| 国产免费男女视频| 国产极品天堂在线| 久久久精品欧美日韩精品| 热99re8久久精品国产| 国产精品电影一区二区三区| 久久久久久伊人网av| 99热这里只有是精品在线观看| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 亚洲人成网站在线观看播放| 精品午夜福利在线看| 人妻久久中文字幕网| 男女啪啪激烈高潮av片| 久久九九热精品免费| 亚洲国产精品合色在线| 久久久久久久午夜电影| 免费观看精品视频网站| 国产精品国产三级国产av玫瑰| 国产成人aa在线观看| 国语自产精品视频在线第100页| 久久久精品94久久精品| 国产色爽女视频免费观看| 成人三级黄色视频| 国产精品一二三区在线看| 天堂中文最新版在线下载 | av黄色大香蕉| 午夜a级毛片| 99热6这里只有精品| av福利片在线观看| 卡戴珊不雅视频在线播放| 天堂av国产一区二区熟女人妻| 搡女人真爽免费视频火全软件| 一级黄片播放器| 女同久久另类99精品国产91| 99久久无色码亚洲精品果冻| 日韩国内少妇激情av| 久久精品综合一区二区三区| 99精品在免费线老司机午夜| 亚洲熟妇中文字幕五十中出| 特大巨黑吊av在线直播| 日韩欧美精品v在线| 成人av在线播放网站| 国产精品一区二区在线观看99 | 黄片无遮挡物在线观看| 日韩一区二区三区影片| 国产精品久久久久久亚洲av鲁大| 成人性生交大片免费视频hd| 寂寞人妻少妇视频99o| 久久韩国三级中文字幕| 人妻制服诱惑在线中文字幕| av在线天堂中文字幕| 亚洲三级黄色毛片| 爱豆传媒免费全集在线观看| 亚洲天堂国产精品一区在线| 亚洲综合色惰| 精品日产1卡2卡| 日韩精品青青久久久久久| 久久久色成人| 久久精品91蜜桃| 丝袜美腿在线中文| 国产毛片a区久久久久| 中国国产av一级| 亚洲在线观看片| 国产蜜桃级精品一区二区三区| 男人狂女人下面高潮的视频| 黄片无遮挡物在线观看| 麻豆av噜噜一区二区三区| 99久久久亚洲精品蜜臀av| 一本久久精品| 99久久九九国产精品国产免费| 久久综合国产亚洲精品| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清在线视频| 久久久成人免费电影| 如何舔出高潮| 黄片wwwwww| 成人鲁丝片一二三区免费| 午夜福利高清视频| 国产真实乱freesex| 亚洲电影在线观看av| 成人漫画全彩无遮挡| 嫩草影院入口| av又黄又爽大尺度在线免费看 | ponron亚洲| 波多野结衣高清作品| 美女 人体艺术 gogo| 久久久欧美国产精品| 自拍偷自拍亚洲精品老妇| 免费大片18禁| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 床上黄色一级片| 美女高潮的动态| 婷婷亚洲欧美| 五月玫瑰六月丁香| 三级毛片av免费| 成人高潮视频无遮挡免费网站| 精品国产三级普通话版| 精品一区二区三区视频在线| 18+在线观看网站| 国产伦精品一区二区三区视频9| 最近中文字幕高清免费大全6| 在线播放无遮挡| 亚洲av成人av| av专区在线播放| 亚洲精品色激情综合| 日本三级黄在线观看| 免费观看a级毛片全部| 免费看av在线观看网站| 国产一区二区三区av在线 | 久久精品国产亚洲av香蕉五月| 欧美极品一区二区三区四区| 亚洲欧美成人综合另类久久久 | 国产精品1区2区在线观看.| 欧美zozozo另类| 国产成人freesex在线| 超碰av人人做人人爽久久| 日韩成人伦理影院| 一个人免费在线观看电影| 97超碰精品成人国产| a级毛片a级免费在线| 天堂av国产一区二区熟女人妻| 有码 亚洲区| 成人欧美大片| 亚洲最大成人av| 久久久久久久久中文| 丰满乱子伦码专区| 久久久久免费精品人妻一区二区| 深爱激情五月婷婷| 久久精品国产亚洲av涩爱 | 欧美三级亚洲精品| 久久久久网色| 国产爱豆传媒在线观看| 久久人人爽人人片av| 国产精品久久久久久亚洲av鲁大| 国产人妻一区二区三区在| 欧美日韩国产亚洲二区| 日韩欧美在线乱码| 免费一级毛片在线播放高清视频| 小说图片视频综合网站| 欧美bdsm另类| 亚洲精品久久国产高清桃花| 欧美一级a爱片免费观看看| 神马国产精品三级电影在线观看| 乱人视频在线观看| 国产69精品久久久久777片| 欧美日韩综合久久久久久| 热99在线观看视频| 久久久久九九精品影院| 精品少妇黑人巨大在线播放 | 蜜臀久久99精品久久宅男| 免费搜索国产男女视频| 国产亚洲av片在线观看秒播厂 | 欧美zozozo另类| 观看免费一级毛片| 日韩一本色道免费dvd| 色综合站精品国产| 久久久国产成人精品二区| 插逼视频在线观看| 黄色视频,在线免费观看| 如何舔出高潮| 日韩欧美 国产精品| av在线天堂中文字幕| av黄色大香蕉| 日日啪夜夜撸| 好男人视频免费观看在线| av专区在线播放| 国产一区二区三区av在线 | 亚洲电影在线观看av| 午夜亚洲福利在线播放| 联通29元200g的流量卡| 国产高清激情床上av| 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 免费黄网站久久成人精品| 菩萨蛮人人尽说江南好唐韦庄 | 午夜免费男女啪啪视频观看| 校园人妻丝袜中文字幕| 波野结衣二区三区在线| 久久午夜亚洲精品久久| 欧美+亚洲+日韩+国产| 国产视频内射| 色噜噜av男人的天堂激情| 日本在线视频免费播放| 国产高清视频在线观看网站| 看黄色毛片网站| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲av涩爱 | 日韩欧美 国产精品| 中文字幕精品亚洲无线码一区| 国产三级在线视频| 内地一区二区视频在线| 丰满人妻一区二区三区视频av| 国产成人精品一,二区 | 一个人看的www免费观看视频| 国产成人影院久久av| 国产精品一区二区性色av| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 免费黄网站久久成人精品| 91麻豆精品激情在线观看国产| 天堂网av新在线| 一级毛片电影观看 | 变态另类成人亚洲欧美熟女| 麻豆国产av国片精品| 亚洲av中文av极速乱| 又爽又黄无遮挡网站| 亚洲美女视频黄频| 成人综合一区亚洲| av在线天堂中文字幕| 黄色一级大片看看| 禁无遮挡网站| 男人舔奶头视频| 亚洲精品456在线播放app| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 日韩欧美 国产精品| 欧美色欧美亚洲另类二区| av在线观看视频网站免费| 美女 人体艺术 gogo| 在线观看av片永久免费下载| 男人舔女人下体高潮全视频| 日韩亚洲欧美综合| 国产伦在线观看视频一区| 黄色视频,在线免费观看| 18+在线观看网站| 欧美性猛交黑人性爽| 99久久精品一区二区三区| 亚洲国产欧洲综合997久久,| 天堂影院成人在线观看| 看非洲黑人一级黄片| 久久精品久久久久久久性| 国产精品综合久久久久久久免费| 免费观看的影片在线观看| 变态另类成人亚洲欧美熟女| 日韩av在线大香蕉| 国产免费男女视频| 一级二级三级毛片免费看| 亚洲人成网站在线观看播放| 国产国拍精品亚洲av在线观看| 老女人水多毛片| www.色视频.com| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 免费看日本二区| 中文欧美无线码| 久久99热6这里只有精品| 欧美+亚洲+日韩+国产| 在线免费十八禁| 日韩欧美 国产精品| 国产老妇伦熟女老妇高清| 日本-黄色视频高清免费观看| 欧美成人一区二区免费高清观看| 麻豆久久精品国产亚洲av| 国产亚洲欧美98| 九九久久精品国产亚洲av麻豆| 精品国内亚洲2022精品成人| 精品免费久久久久久久清纯| 欧美三级亚洲精品| 日本黄色片子视频| av免费在线看不卡| 看黄色毛片网站| 免费人成视频x8x8入口观看| 边亲边吃奶的免费视频| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 中国美白少妇内射xxxbb| 精品久久久久久久末码| 日韩一本色道免费dvd| 成人亚洲精品av一区二区| 国产成年人精品一区二区| 亚洲无线观看免费| videossex国产| 大型黄色视频在线免费观看| 亚洲人成网站在线播| 国产国拍精品亚洲av在线观看| 内射极品少妇av片p| 国产成人影院久久av| 日日啪夜夜撸| 别揉我奶头 嗯啊视频| 成年女人永久免费观看视频| 99久久精品一区二区三区| 成人特级av手机在线观看| 精品人妻一区二区三区麻豆| 国产欧美日韩精品一区二区| 日日撸夜夜添| 国产乱人偷精品视频| 爱豆传媒免费全集在线观看| 久久久精品大字幕| 99久久精品热视频| 国产91av在线免费观看| 色噜噜av男人的天堂激情| 99久国产av精品| 搡女人真爽免费视频火全软件| 国产在线精品亚洲第一网站| 99久久中文字幕三级久久日本| 99久久人妻综合| 精品午夜福利在线看| 小蜜桃在线观看免费完整版高清| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 成年女人永久免费观看视频| 国产精品人妻久久久久久| 97在线视频观看| 国产成人一区二区在线| 亚洲国产精品成人综合色| 亚洲电影在线观看av| 成人性生交大片免费视频hd| 午夜福利高清视频| 成人性生交大片免费视频hd| 亚洲精品自拍成人| 非洲黑人性xxxx精品又粗又长| 悠悠久久av| 麻豆av噜噜一区二区三区| 搡女人真爽免费视频火全软件| 久久久久久久久久久丰满| 国产精品人妻久久久影院| 黄色欧美视频在线观看| 床上黄色一级片| 国产成人一区二区在线| 久久韩国三级中文字幕| 国产一级毛片七仙女欲春2| 久久精品夜色国产| 午夜福利在线观看免费完整高清在 | 18禁在线无遮挡免费观看视频| 国产精品一区二区性色av| 国产日本99.免费观看| 国产午夜精品久久久久久一区二区三区| 免费看av在线观看网站| 国内揄拍国产精品人妻在线| 国产成人精品一,二区 | 欧美高清成人免费视频www| 国产高清有码在线观看视频| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 少妇人妻一区二区三区视频| 欧美另类亚洲清纯唯美| www.色视频.com| 亚洲在线观看片| 激情 狠狠 欧美| 网址你懂的国产日韩在线| 99在线视频只有这里精品首页| 国产乱人视频| 人妻久久中文字幕网| 欧美精品国产亚洲| 99久久精品一区二区三区| 嫩草影院精品99| 日本熟妇午夜| 亚洲欧美日韩卡通动漫| 一个人看的www免费观看视频| 国产人妻一区二区三区在| 最好的美女福利视频网| 亚洲人成网站在线播| 99久久中文字幕三级久久日本| 免费人成在线观看视频色| 日本爱情动作片www.在线观看| 人体艺术视频欧美日本| 国产黄色视频一区二区在线观看 | 国产高清激情床上av| 啦啦啦韩国在线观看视频| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 国产毛片a区久久久久| 18禁在线无遮挡免费观看视频| 国产极品精品免费视频能看的| 亚洲国产色片| 国产三级中文精品| 91久久精品国产一区二区成人| 亚洲欧美成人精品一区二区| 一区二区三区四区激情视频 | 国产精品三级大全| av在线亚洲专区| 中文亚洲av片在线观看爽| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线 | 国产成人一区二区在线| 性插视频无遮挡在线免费观看| 插逼视频在线观看| 简卡轻食公司| 村上凉子中文字幕在线| 蜜桃亚洲精品一区二区三区| 男人舔女人下体高潮全视频| 真实男女啪啪啪动态图|