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

    生境異質(zhì)性對(duì)鼎湖山常綠闊葉林群落功能多樣性的影響

    2016-05-09 05:27:20周光霞黃立新臧曉蔚葉萬(wàn)輝
    廣西植物 2016年2期
    關(guān)鍵詞:鼎湖山葉面積樣地

    周光霞,黃立新,臧曉蔚,韋 霄,葉萬(wàn)輝,沈 浩?

    (1.中國(guó)科學(xué)院退化生態(tài)系統(tǒng)植被恢復(fù)與管理重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院華南植物園,廣州510650;2.中國(guó)科學(xué)院大學(xué),北京100049;3.廣西壯族自治區(qū)廣西植物研究所,廣西桂林541006)中國(guó)科學(xué)院

    生境異質(zhì)性對(duì)鼎湖山常綠闊葉林群落功能多樣性的影響

    周光霞1,2,黃立新1,2,臧曉蔚1,2,韋 霄3,葉萬(wàn)輝1,沈 浩1?

    (1.中國(guó)科學(xué)院退化生態(tài)系統(tǒng)植被恢復(fù)與管理重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院華南植物園,廣州510650;2.中國(guó)科學(xué)院大學(xué),北京100049;3.廣西壯族自治區(qū)廣西植物研究所,廣西桂林541006)
    中國(guó)科學(xué)院

    物種共存機(jī)制一直以來(lái)是群落生態(tài)學(xué)的研究熱點(diǎn)。為了探討異質(zhì)生境條件下鼎湖山常綠闊葉林群落功能多樣性變化,找到其變化的主要環(huán)境驅(qū)動(dòng)因子,該研究利用位于鼎湖山20 hm2監(jiān)測(cè)樣地第2次群落調(diào)查數(shù)據(jù)并選擇代表不同生境(海拔和地形)的27個(gè)樣方(20 m×20 m),于2013年夏季在樣地內(nèi)所選樣方中測(cè)定所有胸徑≥1 cm樹種的葉片功能性狀。所測(cè)性狀包括形態(tài)學(xué)性狀(比葉面積、葉片干物質(zhì)含量、葉面積以及葉片長(zhǎng)寬比)和化學(xué)計(jì)量學(xué)性狀(葉片碳、氮、磷的含量),結(jié)合地形和土壤數(shù)據(jù)并通過(guò)分析功能多樣性隨環(huán)境梯度的變化,探討了環(huán)境過(guò)濾和競(jìng)爭(zhēng)在鼎湖山群落物種共存中的相對(duì)重要性。結(jié)果表明:功能分歧度和群落權(quán)重平均值與環(huán)境因素關(guān)系密切,尤其是海拔、凹凸度和土壤養(yǎng)分。環(huán)境條件較好區(qū)域(微尺度高海拔、高凹凸度和土壤養(yǎng)分含量)的植物采取統(tǒng)一的養(yǎng)分有效保存(低SLA,高LDMC)的適應(yīng)策略(功能分歧度低),環(huán)境過(guò)濾所起作用更強(qiáng);植物在相反的環(huán)境條件下,采取快速生長(zhǎng)策略(高SLA,低LDMC),能夠更好地適應(yīng)環(huán)境的變化,且性狀變化是多樣的(功能分歧度高),在該條件下競(jìng)爭(zhēng)作用更為顯著。葉面積和葉片氮含量的分歧度在環(huán)境條件較好的區(qū)域增大,這與其他功能性狀不一致,說(shuō)明不同生態(tài)位軸(環(huán)境因素)影響不同性狀的分歧度變化,并且在局域尺度上植物為了更好地適應(yīng)環(huán)境變化采取了多樣的適應(yīng)策略。

    群落生態(tài)學(xué),功能多樣性,生境異質(zhì)性,群落共存機(jī)制,適應(yīng)性策略

    群落構(gòu)建機(jī)制(community assembly mechanisms)一直以來(lái)是生態(tài)學(xué)研究的熱點(diǎn)問(wèn)題,其中環(huán)境過(guò)濾(environmental filtering)和競(jìng)爭(zhēng)(competitive interac-tions)是影響群落構(gòu)建最主要的兩種生態(tài)學(xué)過(guò)程(Cornwell et al,2006;Kraft et al,2008;Siefert et al,2013)。一般認(rèn)為,環(huán)境過(guò)濾會(huì)導(dǎo)致物種的適應(yīng)性策略趨同,而競(jìng)爭(zhēng)則導(dǎo)致物種的適應(yīng)性策略趨異(Spasojevic&Suding,2012;Chalmandrier et al,2013;Purschke et al,2013)。功能多樣性(functional diversity)是指功能性狀在某個(gè)特定群落內(nèi)的值、范圍和相對(duì)多度(Villeger et al,2008;Car-reno-Rocabado et al,2012)。功能多樣性與群落構(gòu)建機(jī)制和生態(tài)系統(tǒng)過(guò)程都有著密切的關(guān)聯(lián)(Mason et al,2012;Conti&Diaz,2013;Mason et al,2013;Roscher et al,2013)。環(huán)境過(guò)濾和競(jìng)爭(zhēng)的相對(duì)強(qiáng)度通常導(dǎo)致不同的群落構(gòu)建方式(Spasojevic&Suding,2012),即功能趨同或者趨異(Kraft et al,2008),兩者的相對(duì)強(qiáng)度可能取決于群落的環(huán)境因素,并隨環(huán)境梯度而變化。目前針對(duì)功能多樣性與環(huán)境因子的相關(guān)性研究,包括功能性狀與大尺度條件下土壤和氣候的變異有顯著相關(guān)性(Siefert et al, 2013;Fortunel et al,2014);隨著群落演替的進(jìn)行,功能多樣性會(huì)顯著升高(Janecek et al,2013;Pur-schke et al,2013)。然而,小尺度下功能多樣性隨環(huán)境梯度變化尚不清楚(Lohbeck et al,2012)。

    植物功能多樣性的測(cè)量方法多樣,但各個(gè)指標(biāo)都是在植物功能性狀和物種豐富度的基礎(chǔ)上進(jìn)行的,采用較多的有Rao指數(shù)(Rao,1982)、功能豐富度(Mason et al,2005)、功能均勻度(Mouillot et al,2005)、群落權(quán)重平均值(Lavorel et al,2008)等。其中受到關(guān)注較多的是功能離散度(functional disper-sion,Fdis)和群落權(quán)重平均值(communtiy weighted mean,CWM),群落離散度是指隨機(jī)挑選的物種間某一相同性狀值之間的差異,通常被用來(lái)表征生態(tài)位互補(bǔ)性(Petchey&Gaston,2006)。群落權(quán)重平均值是指物種的功能性狀值加權(quán)相對(duì)多度后表征群落中該性狀的平均值,通常被用來(lái)評(píng)估群落的動(dòng)態(tài)以及生態(tài)系統(tǒng)特性(Enoki&Abe,2004;Louault et al,2005)。從功能性狀的群落權(quán)重平均值和功能離散度的雙重角度出發(fā),能夠更加準(zhǔn)確地判斷隨著環(huán)境變化,群落水平的功能性狀以及性狀離散度水平的變化趨勢(shì)。

    同一群落中從性狀的不同角度出發(fā)會(huì)出現(xiàn)一系列不同的功能多樣性指標(biāo)(Ricotta,2005;Petchey &Gaston,2006)。由于不同功能性狀可能與不同生態(tài)過(guò)程相關(guān),當(dāng)不同的單一性狀與生態(tài)位軸的關(guān)系相反時(shí),這些相反的關(guān)系有可能造成在多元功能多樣性的計(jì)算中相互抵消,從而掩飾多元功能多樣性與環(huán)境的關(guān)聯(lián)。某些群落的多元功能多樣性與環(huán)境梯度之間無(wú)顯著關(guān)聯(lián)(Cornwell et al,2006),這會(huì)極大地影響群落構(gòu)建機(jī)制的判斷(Cornwell&Ack-erly,2009;Spasojevic&Suding,2012)。所以本文從植物葉片形態(tài)學(xué)性狀和理化性質(zhì)出發(fā),選擇具有代表性的功能性狀來(lái)探索功能多樣性的變化。

    森林群落中,生境異質(zhì)性通常是小尺度變異引起的(如地形和土壤的變化) (Harms et al,2001;John et al,2007),生境異質(zhì)性導(dǎo)致樣地能量、水分和養(yǎng)分的可利用性產(chǎn)生差異(Bellingham et al,2000),從而對(duì)群落中物種的分布產(chǎn)生直接或間接影響。鼎湖山大樣地作為中國(guó)森林生物多樣性監(jiān)測(cè)網(wǎng)絡(luò)(Chinese Forest Biodiversity Monitoring Network)和美國(guó)Smithsonian研究院CTFS(Center for Tropical Forest Science)最重要的大型監(jiān)測(cè)樣地之一,建成于2005年(葉萬(wàn)輝等,2008),2010年完成了首次復(fù)查(Shen et al,2013)。它作為南亞熱帶常綠闊葉林的典型代表,海拔230~470 m,地形復(fù)雜,存在較強(qiáng)的環(huán)境異質(zhì)性(葉萬(wàn)輝等,2008)。關(guān)于鼎湖山大樣地群落構(gòu)建機(jī)制,已在群落(Li et al,2009;Wang et al,2009)、種群遺傳(Wang et al,2012)、進(jìn)化(Pei et al,2011)等方面開展了研究。生境異質(zhì)性、種子性狀個(gè)體間的相互作用等對(duì)物種分布格局的形成以及生物多樣性的維持起著重要作用,而功能多樣性的變化正是植物對(duì)生境異質(zhì)性的適應(yīng)性體現(xiàn),在群落構(gòu)建中起到非常重要的作用。因此,本研究從海拔和地形的角度,探討了鼎湖山樣地異質(zhì)生境條件下群落功能多樣性的變化,為物種共存的重要生態(tài)學(xué)過(guò)程提供理論支持,同時(shí)也為該地區(qū)植被恢復(fù)以及生態(tài)系統(tǒng)的保護(hù)提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)自然概況

    鼎湖山國(guó)家級(jí)自然保護(hù)區(qū)位于廣東省肇慶市,112°30′39″~112°33′41″E,23°09′21″~23°11′30″N,地處亞熱帶季風(fēng)氣候區(qū)南緣,北回歸線附近,形成并保存了較完好的地帶性植被類型——南亞熱帶常綠闊葉林,是我國(guó)植物種類最豐富的植被類型之一,具有熱帶向亞熱帶的過(guò)渡性質(zhì),群落結(jié)構(gòu)復(fù)雜,物種組成豐富。年平均溫度為20.9℃,年平均降雨量為1 927 mm,4-9月是主要的降雨季節(jié),也是植物生長(zhǎng)的主要季節(jié),鼎湖山的地帶性土壤為赤紅壤。

    鼎湖山大樣地位于鼎湖山保護(hù)區(qū)的核心區(qū),地形復(fù)雜,海拔為230~470 m,坡度為30°~50°,地形起伏較大,水熱條件良好,相對(duì)濕度較大。植被類型為南亞熱帶常綠闊葉林,群落的發(fā)育和保存良好(葉萬(wàn)輝等,2008)。

    1.2 研究方法

    本研究中,生境異質(zhì)性以海拔和地形因子加以表征。在鼎湖山大樣地中選擇具有代表性的3條山體,從同一山體的不同海拔高度(308.42~413.4 m)、不同凹凸度(-6.16~6.71 m)出發(fā),于每條山體不同海拔(低海拔、中海拔、高海拔)和不同凹凸度的樣方(山谷、山坡、山脊)中選取Jacard相似度最高的樣方(20 m×20 m)作為本次研究對(duì)象,每條山脈同一海拔同一凹凸度選擇3個(gè)樣方,共包括27個(gè)樣方,代表9種生境類型(表1)。

    1.2.1 環(huán)境因子測(cè)定 地形參數(shù)選取海拔高度、凹凸度、坡度和坡向;土壤參數(shù)選取土壤含水量、容重、pH、土壤有機(jī)質(zhì)以及土壤中鉀、磷、氮的含量,數(shù)據(jù)來(lái)源于鼎湖山20 hm2樣地的測(cè)定數(shù)據(jù)。地形因子的測(cè)定是測(cè)量每個(gè)20 m×20 m樣方4個(gè)角的海拔均值為此樣方的海拔;凹凸度為某個(gè)樣方的海拔減去相鄰8個(gè)樣方的海拔平均值;坡度為樣方的任意3個(gè)角形成4個(gè)三角形平面偏離水平面角度的均值。土壤因子測(cè)定時(shí)在20 hm2樣地內(nèi)每隔30 m系統(tǒng)隨機(jī)布設(shè)土壤采樣點(diǎn),在每個(gè)采樣點(diǎn)中,隨機(jī)選擇8個(gè)方向中的一個(gè)方向,在2、5或15 m三個(gè)距離中隨機(jī)增加2個(gè)采樣點(diǎn)。在鼎湖山共布設(shè)了714個(gè)采樣點(diǎn)(Lin et al,2013)。取表層0~10 cm的表土用于測(cè)定土壤指標(biāo)。據(jù)此,本研究每個(gè)樣方涉及到的樣點(diǎn)為2~3個(gè),共統(tǒng)計(jì)了64個(gè)樣點(diǎn)的土壤數(shù)據(jù)。

    1.2.2 功能性狀的測(cè)定 對(duì)選定樣方中所有可以采到的DBH≥1 cm植株,多于10株的物種隨機(jī)采10株作為代表,不足10株物種則全部采集。每植株分別選取4~6片向陽(yáng)成熟的葉片,先稱量葉片鮮重(fresh weight,FW),再用螺旋測(cè)微器測(cè)量植物的葉片厚度(leaf thickness,LT),用LI-3000C便攜式葉面積儀(Li-Cor,NE,USA)測(cè)量葉片面積(leaf area,LA)以及葉片長(zhǎng)度(leaf length,LL)、葉片寬度(leaf width,LW),最后將葉片放入70℃的烘箱中烘干至恒重,并稱量其干重(dry weight,DW)。比葉面積(specific leaf area,SLA)為葉面積與其干重的比值;葉片干物質(zhì)含量(leaf dry matter content,LDMC)為葉片干重與鮮重的比值。樣品烘干后磨碎過(guò)100目篩,用鉬抗比色法測(cè)定單位質(zhì)量葉片磷含量;碳氮含量用元素分析儀(Elementar Analysensysteme GmbH,Hanau,Germany)測(cè)定。

    表1 鼎湖山樣地異質(zhì)生境條件下環(huán)境因子指標(biāo)Table 1 Environmental factors of heterogeneous habitats in Dinghu Mountain plot

    本文選取功能離散度(Functional dispersion,FDis)和群落的權(quán)重平均值(Communtiy-Weighted Mean,CWM)這兩種功能多樣性來(lái)解釋群落構(gòu)建機(jī)制,其計(jì)算方法參考Spasojevic et al(2012)。

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

    所有相關(guān)性分析均采用軟件SPSS 18.0(SPSS Inc.Chicago,IL,USA)進(jìn)行統(tǒng)計(jì)分析,當(dāng)P<0.05則表示有統(tǒng)計(jì)學(xué)意義。功能離散度(Fdis)和群落權(quán)重平均值(CWM)等功能多樣性參數(shù)的計(jì)算采用R語(yǔ)言(3.0.1)(Team,2012),調(diào)用R語(yǔ)言包“FD”完成(Laliberte&Legendre,2010)。采用CANOCO 4.5 對(duì)Fdis和CWM與環(huán)境因子的關(guān)系進(jìn)行多元典范冗余分析(Redundancy analysis,RDA)。

    2 結(jié)果與分析

    2.1 異質(zhì)生境中的環(huán)境因子

    在鼎湖山樣地中,從所選取的樣方的環(huán)境因子中可以看出,低海拔、中海拔與高海拔之間均有50 m的微尺度高度差,凹凸度(-6.16~6.71 m),分別代表了山谷、山坡和山脊的生境。樣地的坡度無(wú)明顯變化,所選樣地均向南(表1)。在土壤因子中,總磷和總氮以及土壤有機(jī)質(zhì)含量變化較大,其他土壤因子變幅不大。表2表明,與土壤因子顯著相關(guān)的是海拔高度和凹凸度,高海拔和凹凸度較大的樣方中,其土壤有機(jī)質(zhì)、有效氮和總氮含量也越高,而土壤pH值越低,高海拔地區(qū)的土壤養(yǎng)分更為充足。

    2.2 樣地群落功能分歧度(FDis)

    不同性狀在異質(zhì)生境中分歧度的變異趨勢(shì)是不同的,其中葉片厚度,比葉面積,葉片磷含量以及葉片中的碳磷比和碳氮比在同一海拔下的變化趨勢(shì)都是山谷>山坡>山脊;而葉片碳含量,葉片磷含量以及葉片氮磷比在同一海拔的不同生境下都是先上升后下降的趨勢(shì)。不同海拔同一地形下性狀的多樣性無(wú)明顯差異(表3)。

    表2 鼎湖山樣地27個(gè)樣方地形因子與土壤因子之間的相關(guān)性分析Table 2 Correlation coefficients between soil characteristics and topological factors in the 27 sampled quadrats in Dinghu Mountain plot

    2.3 樣地地形和土壤因子與功能分歧度

    鼎湖山樣地群落功能分歧度與地形因子以及土壤因子的相關(guān)性分析(表4)結(jié)果顯示,隨著海拔升高,葉片干物質(zhì)含量的分歧度變化降低,性狀趨于一致,凹凸度的增加也會(huì)引起比葉面積、葉片干物質(zhì)含量、葉片磷含量的分歧度顯著降低,坡向與葉面積的分歧度變化是正相關(guān)。

    土壤因子與不同功能性狀分歧度變化的相關(guān)性結(jié)果顯示,土壤pH的升高會(huì)引起葉片干物質(zhì)含量及葉片磷含量的分歧度增加;土壤有機(jī)質(zhì)含量的增加則會(huì)引起葉片干物質(zhì)含量的分歧度變化降低;土壤中有效氮和總氮含量增加,葉片干物質(zhì)含量的分歧度變化降低,但葉片氮含量的分歧度變化增加。

    對(duì)群落功能分歧度與地形因子以及土壤因子的RDA約束排序結(jié)果見圖1。凹凸度對(duì)功能分歧度的解釋力最大(圖1:A),與葉片干物質(zhì)含量以及葉片磷含量的分歧度變化呈負(fù)相關(guān),海拔與葉片氮含量的分歧度變化呈正相關(guān);從圖1:B看出,土壤pH、容重、有機(jī)質(zhì)含量和氮素含量對(duì)功能分歧度的影響較大。土壤pH、容重與有機(jī)質(zhì)、有效氮和總氮含量與不同功能分歧度的變化相反,前者與葉片厚度,葉片干物質(zhì)含量以及比葉面積的分歧度變化都呈正相關(guān),而后者與這些性狀的分歧度變化為負(fù)相關(guān)。

    2.4 樣地葉片功能性狀群落權(quán)重平均值(CWM)

    群落權(quán)重平均值的結(jié)果綜合分析了植物的功能性狀及其在該群落中的多度,能充分表明性狀在群落間的變化。從鼎湖山各群落權(quán)重平均值的結(jié)果中看出,隨著海拔升高,葉片厚度、比葉面積以及葉面積的群落權(quán)重平均值均減小,葉片干物質(zhì)含量,葉片碳含量的群落權(quán)重平均值均增加,而同一海拔不同地形下不同性狀之間幾乎無(wú)差異(表5)。

    2.5 樣地地形和土壤因子與群落權(quán)重平均值的關(guān)系

    CWM與環(huán)境因子的相關(guān)性分析結(jié)果顯示(表6),影響群落中不同性狀的CWM主要是海拔和凹凸度,且海拔和凹凸度對(duì)不同性狀的群落權(quán)重平均值的影響是一致的,隨著海拔和凹凸度的增加,葉片厚度、比葉面積、葉面積的群落權(quán)重平均值均降低,而葉片干物質(zhì)含量、葉片的長(zhǎng)寬比以及葉片碳、氮含量的群落權(quán)重平均值均增加。影響群落權(quán)重平均值的主要土壤因子是土壤容重、pH、有機(jī)質(zhì)以及有效氮和總氮含量。其中土壤容重和pH對(duì)群落權(quán)重平均值的影響是一致的,隨著土壤容重和pH的增加,比葉面積和葉面積的群落權(quán)重平均值均增加,而葉片干物質(zhì)含量以及葉片碳含量的群落權(quán)重平均值均降低;土壤有機(jī)質(zhì)含量和土壤有效氮以及總氮含量對(duì)性狀群落權(quán)重平均值的影響與前述的土壤容重和pH則相反,隨著土壤有機(jī)質(zhì)含量和土壤有效氮以及總氮含量的增加,比葉面積、葉面積的群落權(quán)重平均值均降低,而葉片干物質(zhì)含量和葉片碳含量的群落權(quán)重平均值均增加。

    群落權(quán)重平均值與地形因子以及土壤因子的RDA約束排序結(jié)果見圖2。從圖2:A可以看出,海拔高度和凹凸度與葉面積、葉片厚度以及比葉面積的群落權(quán)重平均值是負(fù)相關(guān)的,與葉片干物質(zhì)含量以及葉片碳含量的群落權(quán)重平均值則是正相關(guān)的。圖2:B中,土壤容重和pH對(duì)性狀權(quán)重平均值的影響與土壤有機(jī)質(zhì)和土壤氮含量的影響恰恰是相反的,土壤容重和pH與葉面積以及比葉面積的群落權(quán)重平均值是正相關(guān)的,土壤有機(jī)質(zhì)和土壤氮含量與其則為負(fù)相關(guān);而土壤有機(jī)質(zhì)和土壤氮含量與葉片干物質(zhì)含量以及葉片碳含量的群落權(quán)重平均值是正相關(guān)的,則土壤容重和pH與其是負(fù)相關(guān)的。

    3 討論

    功能多樣性越來(lái)越多地被用來(lái)檢驗(yàn)環(huán)境過(guò)濾和競(jìng)爭(zhēng)排斥導(dǎo)致的群落物種共存(Mouchet et al,2010;Mason et al,2013;Siefert et al,2013)。一般認(rèn)為,在環(huán)境過(guò)濾起主要作用的條件下功能分歧度會(huì)減小,而在競(jìng)爭(zhēng)排斥較為明顯的條件下功能分歧度會(huì)增加(Cornwell&Ackerly,2009)。本研究結(jié)果表明,除葉面積和葉片氮含量外,其他功能性狀在共存的物種中呈現(xiàn)出在高海拔,凹凸度較大以及土壤養(yǎng)分較高的環(huán)境條件下表現(xiàn)出功能趨同,環(huán)境過(guò)濾的作用更明顯,且在該資源相對(duì)豐富的條件下群落權(quán)重平均值中出現(xiàn)了低SLA,高LDMC,這是一種有效的養(yǎng)分保存策略。鼎湖山資源相對(duì)匱乏的條件下(低海拔,凹凸度較小,土壤養(yǎng)分較低),由于競(jìng)爭(zhēng)作用加強(qiáng),功能分歧度增加,且在該條件下快速生長(zhǎng)(高SLA,低LDMC)的生存策略更容易在競(jìng)爭(zhēng)中取勝。本研究還發(fā)現(xiàn),葉面積和葉片氮含量的變化與其他性狀是不一致的,在Spasojevic&Suding (2012)的研究中也出現(xiàn)了兩種性狀(葉面積和樹高)與其他性狀(比葉面積、氣孔導(dǎo)度、葉綠素含量)的分歧度變化不一致的結(jié)果,說(shuō)明植物在局域環(huán)境下采取了多樣的適應(yīng)策略,且不同性狀的變化有差異。

    表3 鼎湖山樣地異質(zhì)生境下各群落功能分歧度指標(biāo)Table 3 Functional dispersion of different communities in heterogeneous environments in Dinghu Mountain plot

    圖1 鼎湖山樣地環(huán)境因子對(duì)群落葉片功能分歧度影響的冗余分析Fig.1 RDA analysis on the relationship between environment factors and leaf functional dispersion in Dinghu Mountain plot

    圖2 鼎湖山樣地環(huán)境因子對(duì)群落權(quán)重平均值影響的冗余分析Fig.2 RDA analysis on the relationship between environment factors and community-weighted mean in Dinghu Mountain plot

    表4 鼎湖山樣地群落功能分歧度與環(huán)境因子之間的相關(guān)系數(shù)Table 4 Correlation coefficients between plant functional dispersion and environmental factors in Dinghu Mountain plot

    表5 鼎湖山樣地異質(zhì)生境下各群落權(quán)重平均值指標(biāo)Table 5 Community-weighted mean of different communities in heterogeneous habitats in Dinghu Mountain plot

    3.1 功能分歧度與環(huán)境因子的相關(guān)性

    目前,在大尺度條件下,環(huán)境因素與功能分歧度之間存在密切關(guān)聯(lián)(Cornwell&Ackerly,2009;Spa-sojevic et al,2012;Siefert et al,2013;Fortunel et al,2014),但在局域尺度下是未知的。

    通過(guò)鼎湖山地形與土壤因子的相關(guān)性分析可以看出,海拔和凹凸度與土壤有機(jī)質(zhì),有效氮含量和總氮含量是負(fù)相關(guān)的,一些研究發(fā)現(xiàn)凹凸度與土壤養(yǎng)分有負(fù)相關(guān)關(guān)系(Noguchi et al,2007;Lan et al,2011;McEwan et al,2011),因而可以對(duì)植物性狀起到影響作用。在鼎湖山樣地中海拔高差不很明顯,但高凹凸度與高光照相關(guān)聯(lián)(Enoki&Abe,2004;Tsujino et al,2006)。此外,從低海拔到高海拔,人類活動(dòng)強(qiáng)度會(huì)減弱(del Moral&Watson,1978),加上微生物活動(dòng)會(huì)增加土壤有機(jī)質(zhì)含量(隋躍宇等,2009)。因此,在鼎湖山高海拔、凹凸度較大的樣地中,土壤養(yǎng)分較低海拔、凹凸度較小的樣地更為豐富。

    表6 鼎湖山樣地群落權(quán)重平均值與環(huán)境因子之間的相關(guān)系數(shù)Table 6 Correlation coefficients between community-weighted mean and environmental factors in Dinghu Mountain plot

    在熱帶森林中,植物對(duì)土壤養(yǎng)分和光照的競(jìng)爭(zhēng)較強(qiáng),生長(zhǎng)在高凹凸度區(qū)域的個(gè)體能夠獲得更多的光照(Enoki&Abe,2004;Tsujino et al,2006)。鼎湖山樣地植物功能性狀中葉面積和葉片氮含量的分歧度是增加的,表明隨著樣地的海拔和凹凸度增加,光照強(qiáng)度將發(fā)生變化(Song et al,2006;Wang et al,2008;Chu et al,2009;Mason et al,2011),植物通過(guò)形成不同的葉面積來(lái)適應(yīng)對(duì)光照的充分利用,而葉片氮含量的多樣性變化是因?yàn)橹参镌跅l件較好的情況下競(jìng)爭(zhēng)地下資源(Spasojevic&Suding,2012),所以在鼎湖山高海拔葉面積和葉片氮含量的分歧度都是增加的。而其他性狀的分歧度變化在條件較好的條件下(高海拔、凹凸度較大)是降低的,即植物采取了趨同的性狀來(lái)適應(yīng)環(huán)境的變化,這是環(huán)境過(guò)濾作用更為明顯的結(jié)果。這些性狀在資源相對(duì)匱乏的條件下(低海拔、凹凸度較小)分歧度是增加的,這是由于在環(huán)境壓力較大的條件下,競(jìng)爭(zhēng)排斥更為顯著(Tilman,1982),所以出現(xiàn)了性狀分歧度在貧瘠的條件下增加。

    功能分歧度的格局與土壤因子也有明顯的關(guān)聯(lián),土壤肥力(如土壤K和P的含量)能夠影響功能分歧度(Mason et al,2012),單一性狀功能分歧度比如樹高、比葉面積和葉綠素含量也與土壤含水量和N含量有顯著的關(guān)聯(lián)(Spasojevic&Suding,2012)。在鼎湖山樣地中,葉片干物質(zhì)含量的分歧度變化受到土壤因子的影響較大,土壤中有機(jī)質(zhì)、有效氮和總氮含量越高的樣地其pH低,這種情況下,葉片干物質(zhì)含量的分歧度變化降低,植物采取了統(tǒng)一的養(yǎng)分保存的適應(yīng)策略。

    3.2 CWM與環(huán)境因子的相關(guān)性

    CWM隨著環(huán)境梯度發(fā)生變化是環(huán)境過(guò)濾在起作用,但這種過(guò)濾作用是植物對(duì)環(huán)境的抵抗能力,另一種適應(yīng)方式是植物通過(guò)競(jìng)爭(zhēng)資源,結(jié)果是適合的性狀或性狀組合得以保留(Cornwell&Ackerly,2009;Lebrija-Trejos et al,2010)。

    從本研究分析結(jié)果來(lái)看,高海拔和高凹凸度條件下,即土壤中有機(jī)質(zhì)含量和有效氮以及總氮含量較高,土壤容重和pH較低的條件下,比葉面積和葉面積的群落權(quán)重平均值都是降低的,而葉片干物質(zhì)含量和葉片碳含量的群落權(quán)重平均值都是增加的。高的凹凸度和土壤肥力更適合有大的葉面積和比葉面積,這與物種具有更高的生長(zhǎng)速率和碳同化速率(Poorter et al,2008;Carreno-Rocabado et al,2012),但鼎湖山與這些研究的結(jié)果是不一致的,在高凹凸度和高土壤肥力條件下,比葉面積和葉面積都是減小的,而葉片干物質(zhì)含量和葉片碳含量都是增加的。SLA和LDMC反映了植物性狀中生物量的快速生產(chǎn)(高SLA,低LDMC物種)和養(yǎng)分的有效保存(低SLA,高LDMC物種)的一個(gè)基本權(quán)衡(Garnier,2001)。在資源豐富的條件下,由于環(huán)境過(guò)濾的作用,植物統(tǒng)一采取了快速生長(zhǎng)的策略(樹高,葉面積和葉綠素含量都是增大的)(Spasojevic et al,2012),而在鼎湖山樣地中出現(xiàn)了葉片干物質(zhì)含量和葉片碳含量增加,這說(shuō)明植物采取了養(yǎng)分的有效保存。結(jié)合其他性狀的變化如葉片碳含量增大、葉面積減小可以看出,在鼎湖山生長(zhǎng)環(huán)境優(yōu)越條件下,植物是通過(guò)延長(zhǎng)其生命、提高葉片碳含量,來(lái)為提高其生命周期打下基礎(chǔ),而不是通過(guò)短期生物量的快速生產(chǎn)策略,物種采取這種生活方式也是為了適應(yīng)小尺度環(huán)境的不斷變化,這也是與許多發(fā)現(xiàn)不一致的地方(Carreno-Rocabado,2012)。主要是因?yàn)槎江h(huán)境變化多樣,采取有效的養(yǎng)分保存策略會(huì)提高其存活的幾率。所以在鼎湖山樣地中條件較為優(yōu)越的條件下,CWM的變化是環(huán)境過(guò)濾的作用,即保留有效保存養(yǎng)分的植物。在資源較為貧瘠的條件下,植物所需資源緊缺,競(jìng)爭(zhēng)作用明顯(Cornwell et al,2009),能夠快速生長(zhǎng)的植物被保留下來(lái),結(jié)合分歧度的變化來(lái)看,這種快速生長(zhǎng)的生活策略的分歧度變化相對(duì)營(yíng)養(yǎng)豐富的條件更高,但統(tǒng)一是高SLA,低LDMC的適應(yīng)策略。

    4 結(jié)論

    一般認(rèn)為,環(huán)境壓力較大的環(huán)境生境過(guò)濾作用更明顯,功能趨同(Cornwell et al,2006),而在資源豐富的環(huán)境條件下,競(jìng)爭(zhēng)排斥作用會(huì)引起功能趨異(Chesson,2000)。但本研究結(jié)果恰恰相反,在微尺度低海拔條件下,資源相對(duì)微尺度高海拔減少,但是植物為了生存,競(jìng)爭(zhēng)作用明顯,功能多樣性增大,結(jié)合CWM的結(jié)果來(lái)看,這些競(jìng)爭(zhēng)力強(qiáng)的植物都是能夠快速生長(zhǎng)(高SLA,低LDMC)的植物。而在高海拔條件相對(duì)優(yōu)越的環(huán)境條件下,植物采取了統(tǒng)一的養(yǎng)分有效保存的生存策略,但葉面積和葉片氮含量的多樣性變化由于光照和土壤養(yǎng)分的作用,在高海拔的分歧度變化是增大的。所以說(shuō),植物在局域尺度下,為了更好地適應(yīng)環(huán)境的變化采取了多樣的適應(yīng)策略,且不同性狀的變化有差異,進(jìn)一步論證了采取多元功能多樣性可能存在的問(wèn)題。

    局域尺度上物種共存不同于大尺度上的一般規(guī)律,理解局域尺度上物種的共存機(jī)制較大尺度更為復(fù)雜。鼎湖山樣地特殊的異質(zhì)生境變化為物種共存的探討提供了一個(gè)嶄新的平臺(tái),同時(shí)也為局域尺度下功能多樣性變化提供理論基礎(chǔ)。

    BELLINGHAM PJ,TANNER EVJ,2000.The influence of topogra-phy on tree growth,mortality,and recruitment in a tropical montane forest[J].Biotropica,32(3):378-384.

    CARRENO-ROCABADO G,PENA-CLAROS M,BONGERS F,et al,2012.Effects of disturbance intensity on species and functional diversity in a tropical forest[J].J Ecol,100(6): 1 453-1 463.

    CHALMANDRIER L,MUNKEMULLER T,GALLIEN L,et al,2013.A family of null models to distinguish between environ-mental filtering and biotic interactions in functional diversity pat-terns[J].J Veg Sci,24(5):853-864.

    CHESSON P,2000.Mechanisms of maintenance of species diversity [J].Ann Rev Ecol Syst,31:343-366.

    CHU CJ,WANG YS,LI Q,et al,2009.Effects of traits,species identity and local environmental conditions on the assessment of interactions:insights from an alpine meadow community[J].J Plant Ecol,2:135-141.

    CONTI G,DIAZ S,2013.Plant functional diversity and carbon storage-an empirical test in semi-arid forest ecosystems[J].J Ecol,101(1):18-28.

    CORNWELL WK,ACKERLY DD,2009.Community assembly and shifts in plant trait distributions across an environmental gradient in coastal California[J].Ecol Monogr,79(1):109-126.

    CORNWELL WK,SCHWILK DW,ACKERLY DD,2006.A trait-based test for habitat filtering:convex hull volume[J]. Ecology,87(6):1 465-1 471.

    DEL MORAL R,WATSON AF,1978.Gradient structure of forest vegetation in the central Washington Cascades[J].Vegetatio,38(1):29-48.

    ENOKI T,ABE A,2004.Saplings distribution in relation to topog-raphy and canopy openness in an evergreen broad-leaved forest [J].Plant Ecol,173(2):283-291.

    FORTUNEL C,PAINE CET,FINE PVA,et al,2014.Environ-mental factors predict community functional composition in Ama-zonian forests[J].J Ecol,102(1):145-155.

    GARNIER E,SHIPLEY B,ROUMET C,et al,2001.A standard-ized protocol for the determination of specific leaf area and leaf dry matter content[J].Funct Ecol,15(5):688-695.

    HARMS KE,CONDIT R,HUBBELL SP,et al,2001.Habitat as-sociations of trees and shrubs in a 50 ha neotropical forest plot [J].J Ecol,89(6):947-959.

    JOHN R,DALLING JW,HARMS KE,et al,2007.Soil nutrients influence spatial distributions of tropical tree species[J]. PNAS,104(3):864-869.

    JANECEK S,DE BELLO F,HORNIK J,et al,2013.Effects of land-use changes on plant functional and taxonomic diversity along a productivity gradient in wet meadows[J].J Veg Sci,24 (5):898-909.

    KRAFT NJB,VALENCIA R,ACKERLY DD,2008.Functional traits and niche-based tree community assembly in an amazonian forest[J].Science,322(5901):580-582.

    LALIBERTE E,LEGENDRE P,2010.A distance-based framework for measuring functional diversity from multiple traits[J].Ecol-ogy,91(1):299-305.

    LAN GY,HU YH,CAO M,et al,2011.Topography related spatial distribution of dominant tree species in a tropical seasonal rain forest in China[J].For Ecol Manag,262(8):1 507-1 513.

    LAVOREL S,GRIGULIS K,MCINTYRE S,et al,2008.Assessing functional diversity in the field-methodology matters! [J].Funct Ecol,22(1):134-147.

    LEBRIJA-TREJOS E,PEREZ-GARCIA EA,MEAVE JA,et al,2010.Functional traits and environmental filtering drive commu-nity assembly in a species-rich tropical system[J].Ecology,91 (2):386-398.

    LI L,HUANG ZL,YE WH,et al,2009.Spatial distributions of tree species in a subtropical forest of China[J].Oikos,118 (4):495-502.

    LIN G,STRALBERG D,GONG G,et al,2013.Separating the effects of environment and space on tree species distribution: from population to community[J].PLos One,8:e56171.

    LOHBECK M,POORTER L,PAZ H,et al,2012.Functional di-versity changes during tropical forest succession[J].Perspect Plant Ecol,14(2):89-96.

    LOUAULT F,PILLAR VD,AUFRERE J,et al,2005.Plant traits and functional types in response to reduced disturbance in a semi-natural grassland[J].J Veg Sci,16(2):151-160.

    MASON NWH,DE BELLO F,DOLEZAL J,et al,2011.Niche o-verlap reveals the effects of competition,disturbance and con-trasting assembly processes in experimental grassland communities[J].J Ecol,99(3):788-796.

    MASON NWH,DE BELLO F,MOUILLOT D,et al,2013.A guide for using functional diversity indices to reveal changes in assembly processes along ecological gradients[J].J Veg Sci,24 (5):794-806.

    MASON NWH,RICHARDSON SJ,PELTZER DA,et al,2012.Changes in coexistence mechanisms along a long-term soil chronosequence revealed by functional trait diversity[J].J Ecol,100(3):678-689.

    MCEWAN RW,LIN YC,SUN IF,et al,2011.Topographic and bi-otic regulation of aboveground carbon storage in subtropical broad-leaved forests of Taiwan[J].For Ecol Manag,262(9): 1 817-1 825.

    MOUCHET MA,VILLEGER S,MASON NWH,et al,2010.Func-tional diversity measures:an overview of their redundancy and their ability to discriminate community assembly rules[J]. Funct Ecol,24(4):867-876.

    MOUILLOT D,MASON WHN,DUMAY O,et al,2005.Functional regularity:a neglected aspect of functional diversity[J].Oeco-logia,142(3):353-359.

    NOGUCHI H,ITOH A,MIZUNO T,et al,2007.Habitat divergence in sympatric Fagaceae tree species of a tropical montane forest in northern Thailand[J].J Trop Ecol,23(05):549-558.

    PEI NC,LIAN JY,ERICKSON DL,et al,2011.Exploring tree-habitat associations in a Chinese subtropical forest plot using a molecular phylogeny generated from DNA barcode loci[J]. PLoS One,6:e21273.

    PETCHEY OL,GASTON KJ,2006.Functional diversity:back to basics and looking forward[J].Ecol Lett,9(6):741-758.

    POORTER L,WRIGHT SJ,PAZ H,et al,2008.Are functional traits good predictors of demographic rates?Evidence from five Neotropical forests[J].Ecology,89(7):1 908-1 920.

    PURSCHKE O,SCHMID BC,SYKES MT,et al,2013. Contrasting changes in taxonomic,phylogenetic and functional diversity during a long-term succession:insights into assembly processes[J].J Ecol,101(4):857-866.

    RAO CR,1982.Diversity and dissimilarity coefficients:a unified approach[J].Theor Popul Biol,21(1):24-43.

    RICOTTA C,2005.A note on functional diversity measures [J].Basic Appl Ecol,6(5):479-486.

    ROSCHER C,SCHUMACHER J,LIPOWSKY A,et al,2013.A functional trait-based approach to understand community assembly and diversity-productivity relationships over 7 years in experimental grasslands[J].Persp Plant Ecol,15(3):139-149.

    SHEN Y,SANTIAGO LS,MA L,et al,2013.Forest dynamics of a subtropical monsoon forest in Dinghushan,China:recruitment,mortality and the pace of community change[J].J Trop Ecol,29(2):131-145.

    SIEFERT A,RAVENSCROFT C,WEISER MD,et al,2013. Functional beta-diversity patterns reveal deterministic community assembly processes in eastern North American trees[J].Glob Ecol Biogeogr,22(6):682-691.

    SONG MH,TIAN YQ,XU XL,et al,2006.Interactions between root and shoot competition among four plant species in an alpine meadow on the Tibetan Plateau[J].Acta Oecol,29 (2):214-220.

    SPASOJEVIC MJ,SUDING KN,2012.Inferring community assem-bly mechanisms from functional diversity patterns:the importance of multiple assembly processes[J].J Ecol,100(3): 652-661.

    SUI YY,JIAO XG,GAO CS,et al,2009.The relationship among or-ganic matter content and soil microbial biomass and soil enzyme activities[J].Chin J Soil Sci,5:1 036-1 039. [隋躍宇,焦曉光,高崇生,等,2009.土壤有機(jī)質(zhì)含量與土壤微生物量及土壤酶活性關(guān)系的研究[J].土壤通報(bào),5:1 036-1 039.]

    TILMAN D,1982.Resource competition and community structure. (Mpb-17) [M].Princeton University Press.

    TSUJINO R,TAKAFUMI H,AGETSUMA N,et al,2006. Variation in tree growth,mortality and recruitment among topo-graphic positions in a warm temperate forest[J].J Veg Sci,17 (3):281-290.

    VILLEGER S,MASON NWH,MOUILLOT D,2008.New multidi-mensional functional diversity indices for a multifaceted framework in functional ecology[J].Ecology,89(8):2 290-2 301.

    WANG YS,CHUA CJ,MAESTRE FT,et al,2008.On the rele-vance of facilitation in alpine meadow communities:an experi-mental assessment with multiple species differing in their ecolog-ical optimum[J].Acta Oecol,33(1):108-113.

    WANG Z,YE WH,CAO HL,et al,2009.Species-topography as-sociation in a species-rich subtropical forest of China[J].Basic Appl Ecol,10(7):648-655.

    WANG ZF,LIAN JY,HUANG GM,et al,2012.Genetic groups in the common plant species Castanopsis chinensis and their as-sociations with topographic habitats[J].Oikos,121(12): 2 044-2 051.

    YE WH,CAO HL,HUANG ZL,et al,2008.Connunity structure of a 20 hm2lower subtropical evergreen broadleaved forest plot in dinghushan,China[J].Chin J Plant Ecol,32(2):274-286. [葉萬(wàn)輝,曹洪麟,黃忠良,等,2008.鼎湖山南亞熱帶常綠闊葉林20公頃樣地群落特征研究[J].植物生態(tài)學(xué)報(bào),32(2):274-286.]

    Effects of habitat heterogeneity on community functional diversity of Dinghu Mountain evergreen broad-leaved forest

    ZHOU Guang-Xia1,2,HUANG Li-Xin1,2,ZANG Xiao-Wei1,2,WEI Xiao3,YE Wan-Hui1,SHEN Hao1?

    (1.Key Laboratory of Vegetation Restoration and Management of Degraded Ecosystems,South China Botanical Garden,Chinese Academy of Sciences,Guangzhou 510650,China;2.University of Chinese Academy of Sciences,Beijing 100049,China;3.Guangxi Institute of Botany,Guangxi Zhuang Autonomous Region and Chinese Academy of Sciences,Guilin 541006,China)

    Species co-existence mechanisms have long been one of the core topics in community ecology.To explore the changes in community functional diversity of Dinghushan evergreen broad-leaf forest under the heterogeneous habitats and to find out the main driving factors of the functional diversity varations,based on the second community survey data of the 20 hm2forest dynamics plot in Dinghushan National Nature Reserve,27 subplots(20 m×20 m)were selected according to the changes of relative elevation and topograph.These 27 subplots can stand for the heterogeneous habitats in Dinghushan ever-green broad-leaved forest.Leaf functional traits of all the trees with DBH(diameter at breast height)≥1 cm in the 27 sub-plots were measured to calculate the functional diversity(functional dispersion and community weighted mean are caculated as the functional diversity in this paper)during the summer of 2013.The measured traits included morphological traits(spe-cific leaf area,leaf dry matter content,leaf area and leaf lenth to leaf width ratio)and stoichiometric traits(leaf carbon,nitrogen,phosphorus content).Combining with the topographic and soil information,we analysed the main environmental factors that influence the variation of functional diversity.To elucidate the relative importance of environmental filtering and competitive interactions for species co-existence of Dinghushan evergreen broad-leaved forest,we inferred the shifts in func-tional diversity patterns across environmental gradients.The results indicated that functional divergence and community-weighted mean were in close relation with environmental factors,especially elevation,convexity and soil fertility.In the op-timal environmental conditions(high relative elevation,high convexity and higher soil nutrient content),plants adopt the similar(functional divergence decreased)adaptive strategy of efficient conservation of nutrient(low SLA,high LDMC),which suggested that environmental filtering played a stronger role;while in the adverse environmental conditions(low rela-tive elevation,low convexity and lower soil nutrient content),the traits were various(functional divergence increased)and use fast growing strategy(high SLA,low LDMC),since competitive interactions might play a key role.The functional di-vergence of leaf area and leaf nitrogen contents increased in better environmental conditions(high relative elevation,high convexity and higher soil nutrient content),and these two traits(leaf area and leaf nitrogen content)were opposite to the others(specfic leaf area,leaf dry matter content,leaf length to leaf width ratio,leaf carbon and phosphorus content).These results indicated that different niche axes(environments)had different importances in the changes in func-tional traits dispersion,and illustrated that we could use plants diverse strategies to adapt to heterogeneous environments in local scale.

    community ecology,functional diversity,habitat heterogeneity,species co-existence mechanisms,adap-tive strategy

    Q948.15,Q945.79

    A

    1000-3142(2016)02-0127-10

    10.11931/guihaia.gxzw201503015

    周光霞,黃立新,臧曉蔚,等.生境異質(zhì)性對(duì)鼎湖山常綠闊葉林群落功能多樣性的影響[J].廣西植物,2016,36(2):127-136

    ZHOU GX,HUANG LX,ZANG XW,et al.Effects of habitat heterogeneity on community functional diversity of Dinghu Mountain evergreen broad-leaved forest[J].Guihaia,2016,36(2):127-136

    2015-04-10

    2015-06-24

    國(guó)家自然科學(xué)基金(31370446)[Supported by the National Natural Science Foundation of China(31370446)]。

    周光霞(1990-),女,安徽蕪湖人,碩士,主要從事植物生理生態(tài)學(xué)研究,(E-mail)zhouscbg@163.com。

    ?通訊作者:沈浩,博士,副研究員,主要從事保護(hù)生態(tài)學(xué)和植物生理生態(tài)學(xué)研究,(E-mail)shenhao@scbg.ac.cn。

    猜你喜歡
    鼎湖山葉面積樣地
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    作物葉面積測(cè)量的研究進(jìn)展
    昆明市主要綠化樹種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    馬奶子葡萄葉面積評(píng)估模型的建立
    鼎湖山針闊葉混交林CO2濃度及通量變化特征
    廣西植物(2019年2期)2019-09-10 07:22:44
    肇慶·鼎湖山
    “光合作用與細(xì)胞呼吸”重點(diǎn)分析和突破
    苧麻葉面積測(cè)定方法比較研究
    北回歸線理念在生態(tài)旅游營(yíng)銷中的應(yīng)用研究
    国产国拍精品亚洲av在线观看| 午夜福利高清视频| 国产熟女xx| av在线老鸭窝| 小说图片视频综合网站| 一级a爱片免费观看的视频| 十八禁网站免费在线| 中亚洲国语对白在线视频| 夜夜夜夜夜久久久久| 中文字幕高清在线视频| 中文字幕av成人在线电影| 欧美黑人欧美精品刺激| 日韩高清综合在线| 男女那种视频在线观看| 欧美三级亚洲精品| 日本 av在线| 狠狠狠狠99中文字幕| 欧美xxxx性猛交bbbb| 免费观看精品视频网站| 欧美一区二区精品小视频在线| 亚洲av免费在线观看| 日韩欧美精品免费久久 | 国产高潮美女av| 色吧在线观看| 久久午夜福利片| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美日韩东京热| 久久人人精品亚洲av| 综合色av麻豆| 日本黄色片子视频| 又黄又爽又免费观看的视频| 热99在线观看视频| 免费人成视频x8x8入口观看| 直男gayav资源| 日韩精品中文字幕看吧| 99久久精品热视频| 18美女黄网站色大片免费观看| 亚洲av不卡在线观看| 亚洲精品乱码久久久v下载方式| 日韩国内少妇激情av| 午夜影院日韩av| 乱人视频在线观看| 国产精品伦人一区二区| 国产精品乱码一区二三区的特点| 亚洲精品日韩av片在线观看| 99国产精品一区二区蜜桃av| 尤物成人国产欧美一区二区三区| 女同久久另类99精品国产91| www.999成人在线观看| 日韩欧美精品免费久久 | 精品久久国产蜜桃| 免费在线观看成人毛片| 国产精品亚洲美女久久久| 人人妻人人澡欧美一区二区| 一区二区三区四区激情视频 | 精品乱码久久久久久99久播| 国产真实乱freesex| 久久久精品大字幕| 亚洲一区二区三区色噜噜| 亚洲18禁久久av| 一本一本综合久久| 18+在线观看网站| 亚洲精品亚洲一区二区| 午夜亚洲福利在线播放| 亚洲美女搞黄在线观看 | 国产精品亚洲美女久久久| 国产午夜精品久久久久久一区二区三区 | 亚洲三级黄色毛片| 91午夜精品亚洲一区二区三区 | 免费在线观看成人毛片| 亚洲最大成人中文| avwww免费| 性插视频无遮挡在线免费观看| 亚洲avbb在线观看| 在线天堂最新版资源| 美女黄网站色视频| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添av毛片 | 免费观看人在逋| 国产精品伦人一区二区| 欧美日韩综合久久久久久 | 90打野战视频偷拍视频| 如何舔出高潮| 亚洲中文字幕日韩| 有码 亚洲区| 国产成人av教育| 欧美高清成人免费视频www| 国产av在哪里看| 我要搜黄色片| 直男gayav资源| 亚洲中文字幕日韩| 少妇的逼好多水| 可以在线观看的亚洲视频| a级一级毛片免费在线观看| 久久精品人妻少妇| 观看美女的网站| 成人亚洲精品av一区二区| 国产白丝娇喘喷水9色精品| 久久99热这里只有精品18| a在线观看视频网站| 国产在线精品亚洲第一网站| 日韩精品青青久久久久久| 乱人视频在线观看| 国内久久婷婷六月综合欲色啪| 久久香蕉精品热| 亚洲自偷自拍三级| 无人区码免费观看不卡| 白带黄色成豆腐渣| 国产高清视频在线观看网站| 成年免费大片在线观看| .国产精品久久| 丰满人妻一区二区三区视频av| 亚洲成av人片免费观看| 看片在线看免费视频| 亚洲精品在线观看二区| 成人午夜高清在线视频| 淫妇啪啪啪对白视频| 中文字幕av成人在线电影| 18禁在线播放成人免费| 免费黄网站久久成人精品 | 美女 人体艺术 gogo| 国产亚洲精品综合一区在线观看| av欧美777| 日韩高清综合在线| netflix在线观看网站| 久久国产乱子免费精品| 久久久久国内视频| av专区在线播放| 白带黄色成豆腐渣| 一区福利在线观看| 亚州av有码| 久久久久久久久大av| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产 | 又黄又爽又刺激的免费视频.| 国产伦在线观看视频一区| 我的老师免费观看完整版| 自拍偷自拍亚洲精品老妇| 国产精品一区二区三区四区久久| 无人区码免费观看不卡| 欧美最新免费一区二区三区 | 日韩欧美在线二视频| 欧美三级亚洲精品| 国产伦人伦偷精品视频| 精品人妻一区二区三区麻豆 | 赤兔流量卡办理| 国产真实伦视频高清在线观看 | 又爽又黄a免费视频| 我的老师免费观看完整版| 听说在线观看完整版免费高清| 99久国产av精品| 女同久久另类99精品国产91| 乱人视频在线观看| 天堂影院成人在线观看| 亚洲真实伦在线观看| 99久久无色码亚洲精品果冻| 赤兔流量卡办理| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 黄色一级大片看看| 丰满乱子伦码专区| 精品久久久久久久久久久久久| 国产免费一级a男人的天堂| 亚洲性夜色夜夜综合| 色av中文字幕| 赤兔流量卡办理| 国产一区二区在线观看日韩| 别揉我奶头~嗯~啊~动态视频| 色5月婷婷丁香| 国产精品av视频在线免费观看| 亚洲成av人片免费观看| 激情在线观看视频在线高清| 综合色av麻豆| 久久人妻av系列| 日韩精品青青久久久久久| 成人午夜高清在线视频| 成人欧美大片| 国产欧美日韩一区二区三| 99精品久久久久人妻精品| 好看av亚洲va欧美ⅴa在| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看| 国产蜜桃级精品一区二区三区| 成人性生交大片免费视频hd| 看十八女毛片水多多多| 亚洲av日韩精品久久久久久密| 成人亚洲精品av一区二区| 三级国产精品欧美在线观看| 天天躁日日操中文字幕| 日韩欧美在线乱码| 黄片小视频在线播放| 国产日本99.免费观看| 网址你懂的国产日韩在线| 国产国拍精品亚洲av在线观看| 亚洲精品456在线播放app | 波多野结衣高清作品| 亚洲欧美激情综合另类| 成人精品一区二区免费| 99久久久亚洲精品蜜臀av| 国产精品一区二区三区四区久久| 别揉我奶头~嗯~啊~动态视频| 天天一区二区日本电影三级| 90打野战视频偷拍视频| 淫妇啪啪啪对白视频| 亚洲欧美日韩无卡精品| 国产美女午夜福利| 精品一区二区三区视频在线观看免费| 岛国在线免费视频观看| 色综合婷婷激情| 国产精品av视频在线免费观看| 日韩 亚洲 欧美在线| 国语自产精品视频在线第100页| 精品久久国产蜜桃| 国产高清三级在线| 婷婷六月久久综合丁香| 99热这里只有是精品在线观看 | 亚洲欧美激情综合另类| 悠悠久久av| 亚洲欧美日韩卡通动漫| 日日摸夜夜添夜夜添av毛片 | 好男人在线观看高清免费视频| 一个人看视频在线观看www免费| 亚洲18禁久久av| 特大巨黑吊av在线直播| 久久欧美精品欧美久久欧美| 久久九九热精品免费| 国产伦在线观看视频一区| 88av欧美| 99久久成人亚洲精品观看| 国产精品影院久久| 亚洲欧美日韩高清专用| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av电影在线进入| 乱码一卡2卡4卡精品| 久久人人精品亚洲av| 在现免费观看毛片| 国产精品1区2区在线观看.| 国产一区二区三区在线臀色熟女| 天堂网av新在线| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 成年女人看的毛片在线观看| 亚洲精品影视一区二区三区av| 亚洲午夜理论影院| 精品久久久久久久久av| 老女人水多毛片| 亚洲国产精品999在线| 午夜精品久久久久久毛片777| 白带黄色成豆腐渣| 亚洲av第一区精品v没综合| 丁香六月欧美| 亚洲熟妇中文字幕五十中出| 国产aⅴ精品一区二区三区波| 午夜影院日韩av| 一本一本综合久久| 在线天堂最新版资源| 亚洲av第一区精品v没综合| 国产激情偷乱视频一区二区| 我的女老师完整版在线观看| 亚洲美女搞黄在线观看 | 国产高清视频在线播放一区| 青草久久国产| 在线观看66精品国产| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 国产免费av片在线观看野外av| 成人鲁丝片一二三区免费| 国产高潮美女av| 午夜老司机福利剧场| 日韩国内少妇激情av| 久久久久久久午夜电影| 窝窝影院91人妻| 91麻豆精品激情在线观看国产| 国产v大片淫在线免费观看| 嫩草影视91久久| 十八禁人妻一区二区| 亚洲欧美日韩东京热| 淫秽高清视频在线观看| 看片在线看免费视频| 97超视频在线观看视频| av黄色大香蕉| 久久婷婷人人爽人人干人人爱| 色视频www国产| 91久久精品电影网| 成年人黄色毛片网站| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 看免费av毛片| www.999成人在线观看| 高清在线国产一区| 国产精品av视频在线免费观看| 性欧美人与动物交配| 在线观看av片永久免费下载| 成人av在线播放网站| 国产老妇女一区| 欧美在线一区亚洲| 国内精品一区二区在线观看| 日韩欧美免费精品| 老女人水多毛片| av中文乱码字幕在线| 老鸭窝网址在线观看| 久久久精品大字幕| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 国产精品自产拍在线观看55亚洲| 日韩精品青青久久久久久| 午夜福利在线在线| 久久99热这里只有精品18| а√天堂www在线а√下载| 亚洲一区高清亚洲精品| 国内精品久久久久精免费| 亚洲国产高清在线一区二区三| 特级一级黄色大片| 亚洲中文字幕日韩| 精华霜和精华液先用哪个| 18禁裸乳无遮挡免费网站照片| 特级一级黄色大片| 免费大片18禁| 91在线观看av| 色综合站精品国产| 成人一区二区视频在线观看| 在线免费观看的www视频| 精品久久久久久成人av| 婷婷色综合大香蕉| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看| 欧美成狂野欧美在线观看| h日本视频在线播放| 日韩欧美免费精品| 亚洲成av人片在线播放无| 国产欧美日韩精品亚洲av| 级片在线观看| 五月伊人婷婷丁香| 老司机午夜福利在线观看视频| 嫩草影视91久久| 性色avwww在线观看| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| 亚洲电影在线观看av| 国产一区二区亚洲精品在线观看| 特大巨黑吊av在线直播| 国产午夜精品久久久久久一区二区三区 | 91麻豆av在线| 精品人妻视频免费看| 精品人妻熟女av久视频| 每晚都被弄得嗷嗷叫到高潮| 久久久色成人| 精品一区二区三区人妻视频| 亚洲成av人片免费观看| 999久久久精品免费观看国产| 精品久久久久久久末码| 波野结衣二区三区在线| 最近最新免费中文字幕在线| 国产黄片美女视频| 性插视频无遮挡在线免费观看| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站 | 久久中文看片网| 免费黄网站久久成人精品 | 国产 一区 欧美 日韩| 色噜噜av男人的天堂激情| 麻豆av噜噜一区二区三区| 此物有八面人人有两片| 99热6这里只有精品| 成人av一区二区三区在线看| 精品人妻熟女av久视频| 欧美+亚洲+日韩+国产| 国产真实乱freesex| 不卡一级毛片| 久9热在线精品视频| 精品一区二区三区视频在线观看免费| 在线观看午夜福利视频| 欧美日韩综合久久久久久 | av女优亚洲男人天堂| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩高清专用| 免费一级毛片在线播放高清视频| 三级男女做爰猛烈吃奶摸视频| 日韩免费av在线播放| 91字幕亚洲| 人人妻人人澡欧美一区二区| 一区福利在线观看| а√天堂www在线а√下载| 午夜精品在线福利| 首页视频小说图片口味搜索| 麻豆成人午夜福利视频| 亚洲天堂国产精品一区在线| 中文资源天堂在线| 欧美日韩综合久久久久久 | 国产精品久久久久久人妻精品电影| 亚洲欧美清纯卡通| 男人狂女人下面高潮的视频| 日本黄色视频三级网站网址| 在现免费观看毛片| av在线蜜桃| 国产欧美日韩一区二区精品| 悠悠久久av| 国产淫片久久久久久久久 | 成人无遮挡网站| 变态另类成人亚洲欧美熟女| 午夜免费激情av| 日韩中文字幕欧美一区二区| 亚洲成av人片免费观看| 日本黄大片高清| 好男人在线观看高清免费视频| 亚洲精品在线美女| 精品一区二区三区视频在线| 久久精品人妻少妇| 亚洲成a人片在线一区二区| 国产男靠女视频免费网站| 亚洲成a人片在线一区二区| 村上凉子中文字幕在线| 国产v大片淫在线免费观看| 亚洲精品成人久久久久久| 欧美成人免费av一区二区三区| 丁香欧美五月| 国产精品久久电影中文字幕| 精品一区二区免费观看| 久久久精品大字幕| 琪琪午夜伦伦电影理论片6080| 免费一级毛片在线播放高清视频| 俄罗斯特黄特色一大片| 日本熟妇午夜| 国产精品久久视频播放| 午夜福利在线观看免费完整高清在 | 久久久久久久久久成人| 99国产精品一区二区蜜桃av| 小说图片视频综合网站| 毛片女人毛片| 老熟妇仑乱视频hdxx| 国产麻豆成人av免费视频| av中文乱码字幕在线| 欧美3d第一页| 婷婷六月久久综合丁香| 国产 一区 欧美 日韩| 亚洲av成人av| 中文字幕av在线有码专区| or卡值多少钱| 欧美精品啪啪一区二区三区| 中亚洲国语对白在线视频| 欧美日韩综合久久久久久 | 欧美黄色淫秽网站| 久久亚洲精品不卡| 日本一本二区三区精品| 色哟哟·www| 人人妻人人看人人澡| 99国产综合亚洲精品| 免费看光身美女| 老鸭窝网址在线观看| 夜夜看夜夜爽夜夜摸| 精品熟女少妇八av免费久了| av中文乱码字幕在线| 亚洲欧美精品综合久久99| 草草在线视频免费看| 乱码一卡2卡4卡精品| 十八禁网站免费在线| 99久久精品一区二区三区| 18禁在线播放成人免费| 激情在线观看视频在线高清| 一本久久中文字幕| 丝袜美腿在线中文| 国产69精品久久久久777片| 国产精品野战在线观看| 亚洲中文字幕一区二区三区有码在线看| 精品国产三级普通话版| 午夜a级毛片| 丝袜美腿在线中文| 热99re8久久精品国产| 美女被艹到高潮喷水动态| 嫩草影视91久久| 精品一区二区三区视频在线| 日韩中文字幕欧美一区二区| 日本 欧美在线| 久久久久国产精品人妻aⅴ院| 久久久精品欧美日韩精品| 国产精品野战在线观看| 长腿黑丝高跟| 国产色爽女视频免费观看| 他把我摸到了高潮在线观看| 99久国产av精品| 国产野战对白在线观看| 无人区码免费观看不卡| 日韩欧美 国产精品| 久久九九热精品免费| 国产精品久久久久久久电影| 18美女黄网站色大片免费观看| 身体一侧抽搐| 日本黄色片子视频| 日韩欧美一区二区三区在线观看| 免费在线观看日本一区| 性色av乱码一区二区三区2| 小说图片视频综合网站| 在线免费观看不下载黄p国产 | 女人十人毛片免费观看3o分钟| 亚洲国产高清在线一区二区三| 在线免费观看不下载黄p国产 | 中亚洲国语对白在线视频| 一区二区三区高清视频在线| 国产精品伦人一区二区| 日韩欧美三级三区| 9191精品国产免费久久| 中亚洲国语对白在线视频| 久久热精品热| 欧美成人免费av一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲无线观看免费| 高清毛片免费观看视频网站| 午夜福利在线在线| 国产在视频线在精品| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 国产午夜精品久久久久久一区二区三区 | 91麻豆精品激情在线观看国产| 国产成人a区在线观看| 床上黄色一级片| 国产欧美日韩精品一区二区| 91狼人影院| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 午夜福利欧美成人| 亚洲人成网站在线播| 国产高清激情床上av| 欧美国产日韩亚洲一区| 男女下面进入的视频免费午夜| www.熟女人妻精品国产| 宅男免费午夜| 免费观看精品视频网站| 精品久久久久久久久av| 免费在线观看日本一区| 有码 亚洲区| 小蜜桃在线观看免费完整版高清| 亚洲内射少妇av| 观看免费一级毛片| 一进一出好大好爽视频| 中出人妻视频一区二区| 久久精品国产亚洲av香蕉五月| 亚洲av成人精品一区久久| 亚洲内射少妇av| 日日摸夜夜添夜夜添av毛片 | 88av欧美| 午夜视频国产福利| 我要搜黄色片| 波野结衣二区三区在线| 国产亚洲精品久久久com| 欧美三级亚洲精品| 狂野欧美白嫩少妇大欣赏| 精品久久久久久,| 亚洲精品乱码久久久v下载方式| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 黄色女人牲交| 天堂影院成人在线观看| 国产91精品成人一区二区三区| 99热这里只有是精品50| 男女床上黄色一级片免费看| 国产精品电影一区二区三区| 国产三级黄色录像| 免费av观看视频| 日韩欧美免费精品| 欧美激情久久久久久爽电影| 欧美xxxx黑人xx丫x性爽| 日本黄色视频三级网站网址| 亚洲国产色片| 久久久久亚洲av毛片大全| 男人狂女人下面高潮的视频| 国内揄拍国产精品人妻在线| 最近最新中文字幕大全电影3| АⅤ资源中文在线天堂| 真人做人爱边吃奶动态| 亚洲av一区综合| 久久久久国产精品人妻aⅴ院| 国模一区二区三区四区视频| 亚洲 欧美 日韩 在线 免费| 成人鲁丝片一二三区免费| 亚洲av免费在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产淫片久久久久久久久 | 欧美3d第一页| 天堂√8在线中文| 成年版毛片免费区| 无人区码免费观看不卡| 丰满人妻一区二区三区视频av| 亚洲va日本ⅴa欧美va伊人久久| 能在线免费观看的黄片| 国产一区二区在线观看日韩| 国产精品野战在线观看| 日韩欧美精品免费久久 | 男女那种视频在线观看| 午夜亚洲福利在线播放| 12—13女人毛片做爰片一| 国产精品永久免费网站| 深爱激情五月婷婷| 不卡一级毛片| 日本三级黄在线观看| 听说在线观看完整版免费高清| 精品一区二区三区人妻视频| 国产久久久一区二区三区| 免费大片18禁| 三级国产精品欧美在线观看| 免费黄网站久久成人精品 | 在线观看免费视频日本深夜| 国产精品亚洲av一区麻豆| 麻豆国产97在线/欧美| 女同久久另类99精品国产91| 欧美中文日本在线观看视频| 精品久久国产蜜桃| 精品不卡国产一区二区三区| 如何舔出高潮| 亚洲av电影在线进入| 三级国产精品欧美在线观看| 直男gayav资源| 久久久久久国产a免费观看|