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

    南亞熱帶闊葉林展葉物候及其種內(nèi)種間差異探討

    2017-04-29 00:00:00崔雪娜杜彥君趙袁郭希梅黃忠良
    廣西植物 2017年3期

    摘要: 植物的展葉物候與植物的生長和生存密切相關(guān),但未見有定性的研究植物展葉物候的種間和種內(nèi)差異與驗(yàn)證亞熱帶常綠闊葉林植物展葉物候系統(tǒng)發(fā)育關(guān)系的報(bào)道。為了更好地理解南亞熱帶常綠闊葉林物種共存機(jī)制和預(yù)測(cè)不同物種對(duì)氣候變化的響應(yīng),該研究于2013年1月至2014年12月期間,在中國科學(xué)院華南植物園鼎湖山自然保護(hù)區(qū)內(nèi)對(duì)44個(gè)樹種的展葉物候進(jìn)行了連續(xù)觀測(cè),利用變異系數(shù)定性分析了氣候因子、系統(tǒng)發(fā)育關(guān)系對(duì)展葉始期及其種內(nèi)和種間差異和展葉期的影響。結(jié)果表明:鼎湖山自然保護(hù)區(qū)植物展葉始期主要集中在3-5月,植物展葉物候受氣候因子影響不顯著;不同物種的展葉物候種內(nèi)變異差別較大,其中種內(nèi)差異最大的是紅皮紫棱(Craibiodendron kwangtungense)(變異系數(shù)為0.74),最小的是廣東潤楠(Machilus kwangtungensis)(變異系數(shù)為0.09);基于APG Ⅲ物種數(shù)據(jù)庫構(gòu)建譜系樹,發(fā)現(xiàn)植物間的系統(tǒng)發(fā)育關(guān)系對(duì)展葉物候無顯著影響,即親緣關(guān)系近的物種,展葉期并不相似。該研究結(jié)果對(duì)于理解生態(tài)系統(tǒng)過程,如碳循環(huán)、物質(zhì)能量流動(dòng),以及展葉物候?qū)夂蜃兓倪m應(yīng)性具有重要意義。

    關(guān)鍵詞: 鼎湖山, 展葉物候, 種內(nèi)種間差異, 變異系數(shù), 譜系保守性

    中圖分類號(hào): Q948.12

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

    文章編號(hào): 10003142(2017)03032207

    Abstract: Leaf unfolding phenology in plants can be determined by genetic factors, biotic and abiotic environmental factors. However, limited information was available about the leaf unfolding phenology of evergreen broadleaved forest, and its phylogenetic conservatism has not been investigated. In order to better understand the maintaining mechanism of species coexistence in lower subtropical evergreen broadleaved forest and to predict species responses to global changes, we used the phenological data of 44 species containing 112 individuals of trees, which obtained from Dinghushan (DHS) Arboretum, Chinese Academy of Sciences (CAS), to explore the effects of climate, phylogenetic and the intra and interspecific variations on leaf unfolding date. The results showed that leaf unfolding date in DHS National Nature Reserve peaked in the period from March to June, and had been little influenced by climate factors. There were distinct intraspecific variation between different species, and Craibiodendron scleranthum var. kwangtungense showed the greatest coefficient of variation (CV, 0.74), while Machilus kwangtungensis had the smallest coefficient of variation (CV, 0.09). According to the phylomatic tree based on APG Ⅲ, no significant effect of phylogeny on the leaf unfolding date was detected, indicating that close related species did not have the tendency to leafout at similar time. This study suggests that exploring the intra and interspecific variations of leaf unfolding phenology is important for understanding the ecosystem processes, such as the carbon cycle, the material and energy flow, and the adaptation of different tree species to climate changes.

    Key words: Dinghushan, leaf unfolding phenology, intra and interspecific variations, coefficients of variation, phylogenetic conservatism

    對(duì)植物展葉、開花、結(jié)實(shí)、落葉等自然現(xiàn)象節(jié)律的研究,稱為物候?qū)W(Phenology)(Polgar Primack,2011)。植物物候受環(huán)境因子驅(qū)動(dòng)進(jìn)行季節(jié)性變化,已被公認(rèn)為是最理想的氣候變化指示因子(Menzel et al,2001;Parmesan Yohe,2003)。不同物種展葉物候的差異可以影響生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動(dòng),如碳循環(huán)、水循環(huán)、生物量積累等(Vitasse et al,2009;Fridley,2012)。展葉時(shí)間的變化可能會(huì)影響專食性(如取食植物幼葉)昆蟲的生活史(欽俊德,1986)等,從而造成鳥類、動(dòng)物、微生物等一系列營養(yǎng)級(jí)發(fā)生改變(Polgar Primack,2011)。目前,對(duì)展葉物候的研究主要集中于溫帶森林的落葉樹種(Cleland et al,2007);對(duì)熱帶和亞熱帶森林常綠樹種的物候研究,多數(shù)為開花和結(jié)實(shí)物候(Changyang et al,2013;Du et al,2015;胡小麗等,2015),如許格希等(2014)利用積分回歸建立模型探討了展葉物候與氣候變化的動(dòng)態(tài)關(guān)系;Zhao et al(2012)利用物候相機(jī)的綠度變化進(jìn)行連續(xù)的定位監(jiān)測(cè)對(duì)三個(gè)植物物種(兩個(gè)常綠闊葉林物種和一個(gè)熱帶季雨林物種)的展葉和葉衰老時(shí)間進(jìn)行探討。然而,對(duì)于亞熱帶常綠闊葉林展葉物候的種內(nèi)和種間差異的研究還很少見有報(bào)道。物種展葉物候的差異顯示出對(duì)早春資源的利用的時(shí)間差異,可以通過影響葉壽命,進(jìn)而影響群落的外貌特征(朱旭斌和孫書存,2006)。因此,研究展葉物候的種內(nèi)及種間變異,將有利于我們理解物種共存機(jī)制,以及更好地預(yù)測(cè)物種對(duì)氣候變化的響應(yīng)(Moran et al,2015)。

    植物物候的系統(tǒng)發(fā)育保守性問題近年來備受關(guān)注(Panchen et al,2014;Du et al,2015)。譜系保守性假說認(rèn)為,親緣關(guān)系越近的物種,物候發(fā)生時(shí)間越集中(Davies et al,2013;Changyang et al,2013;Du et al,2015)。對(duì)于展葉物候的譜系保守性問題,目前只有Davies et al(2013)和Panchen et al(2014)探討了展葉物候的系統(tǒng)發(fā)育保守性。Davies et al(2013)通過北半球6個(gè)溫帶地區(qū)站點(diǎn)的觀測(cè)數(shù)據(jù)得出親緣關(guān)系相近的物種傾向于在相似的時(shí)間內(nèi)展葉;Panchen et al(2014)對(duì)溫帶落葉林展葉物候的研究也證實(shí)了譜系保守性假說。然而,他們的研究區(qū)域都在溫帶,研究的樹種都是落葉種。那么,在亞熱帶常綠闊葉林內(nèi)是否展葉物候也具有譜系保守性,還未見有報(bào)道。

    鼎湖山國家級(jí)保護(hù)區(qū),地處亞熱帶季風(fēng)氣候區(qū)南緣,是北回歸沙漠帶上不可多得的“綠洲”,其獨(dú)特的地理位置、氣候條件和人文因素,形成并較好地保存了具有熱帶向亞熱帶的過渡性質(zhì)的地帶性植被類型—南亞熱帶常綠闊葉林,群落結(jié)構(gòu)相對(duì)復(fù)雜,組成種類相對(duì)豐富,成為當(dāng)今地球該緯度帶上最具特色、最具研究?jī)r(jià)值的地區(qū)之一(孔國輝等,1998)。本文基于鼎湖山南亞熱帶季風(fēng)常綠闊葉林2013年和2014年兩年的物候觀測(cè)數(shù)據(jù),通過變異系數(shù)定性分析植物的種間和種內(nèi)物候差異嘗試回答以下3個(gè)問題:(1)南亞熱帶森林群落展葉物候是否與氣候因子相關(guān)?(2)展葉始期的種內(nèi)和種間差異大小如何?(3)物種間展葉始期是否具有譜系保守性?

    1材料與方法

    1.1 研究區(qū)概況

    鼎湖山國家級(jí)自然保護(hù)區(qū)位于廣東省肇慶市鼎湖區(qū),112°30′39″~112°33′41″ E,23°09′21″~23°11′30″ N。保護(hù)區(qū)地勢(shì)西北高、東南低,大部分地貌為山地、丘陵,山體陡峭,坡度多在30°~45°之間。最低海拔14.1 m,最高峰雞籠山海拔為 1 000.3 m,總面積約1 155 hm2(葉萬輝等,2008)。該地區(qū)屬南亞熱帶濕潤季風(fēng)氣候,冬夏氣候交替明顯,據(jù)2013年和2014年觀測(cè)資料,最熱月(6、7 月)平均氣溫分別為28.28、29.19 ℃,最冷月(1、2月)平均氣溫為13.71、13.36 ℃,年平均氣溫分別為22.18、22.23 ℃,年均降雨量分別為2 036、1 979 mm。80%的雨量集中在4-9月降落,11 月至翌年1 月為旱季,干濕季節(jié)明顯。年均蒸發(fā)量為1 115 mm,年均相對(duì)濕度為80. 3%(黃忠良等,1998)。

    1.2 研究路線及樹種選取

    為了廣泛覆蓋保護(hù)區(qū)面積且易于觀測(cè),本研究選取五條觀測(cè)路線(圖1),分別為辦公區(qū)+東溪路線、西溪路線、慶云寺小道路線和大樣地路線,共包括44種南亞熱帶常綠闊葉林具有代表性物種,進(jìn)行為期兩年(2013年1月至2014年12月)的物候觀測(cè)研究。詳細(xì)坐標(biāo)信息見附件1。

    1.3 物候和氣象資料來源

    結(jié)合鼎湖山樹木園物候觀測(cè)資料的實(shí)際情況,選取16種具有南亞熱帶森林代表性的天然野生分布樹種,進(jìn)行植物展葉物候的種內(nèi)種間差異,選取44種闊葉林常綠樹種進(jìn)行譜系保守性研究。同期的氣象數(shù)據(jù)來自于鼎湖山定位站的長期定位觀測(cè),定位站位于112°32′57″ E,23°09′50″ N。

    1.4 物候觀測(cè)方法

    依據(jù)《中國物候觀測(cè)方法》(宛敏渭和劉秀珍,1979),在既定路線中,選取易于觀測(cè)的個(gè)體,掛牌標(biāo)記。目測(cè)與望遠(yuǎn)鏡結(jié)合,每條路線固定觀測(cè)者,每3天觀測(cè)一次,物候發(fā)生期酌情增加觀測(cè)頻度。展葉始期的判別標(biāo)準(zhǔn)是當(dāng)被觀測(cè)植株的新芽從芽苞中伸出卷曲著或按葉脈褶疊著的小葉,出現(xiàn)第一批有一、二片舒展開的的葉片時(shí),只要有一兩個(gè)小葉片舒展開,即記為該觀測(cè)植株的展葉始期(宛敏渭和劉秀珍,1979),物種展葉始期為該物種所有觀測(cè)植株展葉始期平均值。

    1.5 數(shù)據(jù)計(jì)算方法

    1.5.1 氣象數(shù)據(jù)處理由于本文并非對(duì)比年際間差異,且一年觀測(cè)數(shù)據(jù)有可能產(chǎn)生較大誤差,故物候期和氣候因素?cái)?shù)據(jù)均取兩年觀測(cè)結(jié)果平均值。

    1.5.2 變異系數(shù)(Coefficients of variation)的計(jì)算

    標(biāo)準(zhǔn)差SD=1n∑ni=1xi-x2

    變異系數(shù) C.D=SDx

    式中,xi為第i棵樹的展葉初始時(shí)間,x為一個(gè)物種全部觀測(cè)個(gè)體展葉始期的平均值,n為該物種觀測(cè)個(gè)體總數(shù)。

    采用儒略日期(Julian day)換算法(1月1日為第1天,1月2日為第2天,依此類推),將展葉時(shí)間轉(zhuǎn)換為距離當(dāng)年1月1日的天數(shù)。

    選取種內(nèi)觀測(cè)個(gè)體多于3個(gè)的物種進(jìn)行展葉物候種內(nèi)差異分析,有16個(gè)物種符合要求。

    1.5.3 構(gòu)建譜系樹及檢驗(yàn)譜系信號(hào)對(duì)44個(gè)物種的展葉時(shí)間進(jìn)行譜系信號(hào)檢驗(yàn)。采用的最常用方法是APG分類系統(tǒng)(Angiosperm Phylogeny Group)( APGⅢ,2009),在(Version 3)(Webb Donoghue,2005)軟件中輸入物種名錄,便可獲得從R20120829.new(http://phylodiversity.net/phylomatic)數(shù)據(jù)庫中構(gòu)建的譜系樹。R20120829.new是基于APGⅢ骨架的譜系框架數(shù)據(jù)庫。然后采用Phylocom 4.2(Webb et al,2008)軟件中提供的“bladj”算法基于已通過分子或者化石定年數(shù)據(jù)庫(Wikstrom et al,2001)對(duì)分支長度進(jìn)行調(diào)整。該方法已經(jīng)廣泛應(yīng)用在物候研究中(Queenborough et al,2009;Panchen et al,2014;Du et al,2015)。

    選擇Blomberg K值法來度量譜系信號(hào)強(qiáng)弱,可檢測(cè)功能性狀和物種間系統(tǒng)發(fā)育關(guān)系的相關(guān)性,該方法已廣泛應(yīng)用于生態(tài)學(xué)和物候?qū)W研究(Jia Du,2011;LessardTherrien et al,2013),由于物候數(shù)據(jù)是非線性的變量,如在循環(huán)統(tǒng)計(jì)中,第365天展葉(一年的最后一天)與第1天(一年的第1天)展葉的差值是1,而不是364。因此,計(jì)算此類性狀譜系信號(hào)時(shí),要基于循環(huán)統(tǒng)計(jì)矯正K值(Du et al,2015)。利用R.3.1.3軟件中picante軟件包(Kembel et al,2010)中的phylosignal函數(shù)計(jì)算K統(tǒng)計(jì)量(Blomberg et al,2003),要先用“Circular”程序包(Agostinelli Lund,2013)中的as. circular函數(shù)對(duì)展葉時(shí)間進(jìn)行轉(zhuǎn)換。

    2結(jié)果與分析

    2.1 展葉物候特征及與氣候因子的關(guān)系

    鼎湖山亞熱帶常綠闊葉林植物展葉物候表現(xiàn)出明顯的季節(jié)性(圖2:A)。展葉始期從1月份持續(xù)到8月份,多集中在3-5月(28種),其中4月份為展葉最高峰(圖2:A)。展葉始期與月平均氣溫(r=0.445,P>0.05)和月平均降雨量(r=0.141,P>0.05)關(guān)系不顯著(圖2:A,B)。

    2.2 16種闊葉喬木樹種展葉時(shí)間的種內(nèi)變異

    有4個(gè)物種變異系數(shù)大于0.5(表1),分別為紅皮紫棱(Craibiodendron kwangtungense)(0.74)、荷木(Schima superba)(0.63)、格木(Erythrophleum fordii)(0.53)、水石梓(Sarcosperma laurinum)(0.51),占全部物種的25%。展葉期波動(dòng)范圍最大的前五個(gè)種分別為紅皮紫棱(209 d)、臀果木(Pygeum topengii)(171 d)、錐栗(Castanea henryi)(170 d)、黃杞(Engelha rdiaroxburghiana)(150 d)、橄欖(Canarium album)(129 d),而展葉始期波動(dòng)范圍最短的是烏欖(Canarium pimela)(16 d)。

    2.3 展葉始期的種間差異及其譜系信號(hào)檢驗(yàn)

    展葉最早的物種是重陽木(Bischofia polycarpa),展葉始期為13(1月13日),展葉最晚種是臀果木,展葉始期為210(7月29日)。44種常綠植物種間差異為0.48,種內(nèi)個(gè)體數(shù)多于3個(gè)的16種喬木樹種種間變異系數(shù)為0.43。

    44種闊葉喬木樹種的展葉始期無顯著的譜系信號(hào)(K=0.271,P=0.501),表明親緣關(guān)系近的物種并沒有相近的展葉始期。

    3討論

    3.1 鼎湖山展葉物候特征

    鼎湖山南亞熱帶常綠闊葉林植物展葉物候具有明顯的季節(jié)性,展葉始期與氣溫和降雨量等氣候因子無顯著相關(guān)關(guān)系。植物展葉始期多發(fā)生在3-5月,與之前熱帶亞熱帶地區(qū)常綠樹種展葉物候研究結(jié)果一致(李明佳和王鑄豪,1984;許格希等,2014)。較溫帶森林的展葉持續(xù)時(shí)間長,如朱旭斌等(2005)對(duì)溫帶落葉櫟林展葉物候研究表明,展葉主要集中在3月中旬到4月。經(jīng)歷1、2月的相對(duì)低溫干旱期,只有較少物種在早春展葉。3月份氣溫增高,降水量增加,大量物種開始展葉,一直持續(xù)到8月份展葉結(jié)束。但本研究中發(fā)現(xiàn)展葉始期與氣溫和降雨并未表現(xiàn)出顯著相關(guān)關(guān)系,原因可能是:(1)因?yàn)樵谀蟻啛釒^(qū)域全年的氣候因子變化較小,植物的展葉受氣候因子的限制較小,而受內(nèi)在的生理機(jī)制調(diào)控較大所致。(2)某些物候期持續(xù)時(shí)間較長,對(duì)氣候因子敏感性較大的只是其中某個(gè)階段,如葉物候的出芽期、 結(jié)實(shí)物候的果熟期等。越接近物候期的時(shí)段,對(duì)物候影響越大。許格希等(2014)對(duì)尖峰嶺12種熱帶常綠闊葉喬木的研究結(jié)果表明,多數(shù)樹種展葉始期對(duì)冬春季節(jié)各月份氣溫上升較為敏感,尤其與展葉前一個(gè)月氣溫上升呈負(fù)相關(guān)關(guān)系。同時(shí),冬春季節(jié)降水是否充沛也直接影響展葉始期的提前或推遲。植物物候受多種因子共同作用,不同地區(qū)可能會(huì)存在不同的差異,因此關(guān)于本地區(qū)更多的物候特征及其調(diào)控因子還需進(jìn)一步驗(yàn)證。

    3.2 展葉物候種內(nèi)和種間差異

    不同物種展葉物候具有較大差異,同一物種在不同生境下展葉時(shí)間也有一定差異,這與前人已有的研究結(jié)果一致(Lechowicz,1984;Sun et al,2006;Liu et al,2011)。本研究結(jié)果顯示,44種常綠植物種間展葉物候有較大差異;16種喬木樹種展葉物候種內(nèi)差異對(duì)比最大的為紅皮紫棱209 d,種內(nèi)差異最小的為烏欖,早晚展葉始期相差16 d。造成展葉時(shí)間種內(nèi)和種間差異的原因有很多,除受環(huán)境因子(如溫度、濕度、降雨、CO2濃度等)控制外(Van Volkenburgh,1999),還受物種本身遺傳特性(Rathcke Lacey,2003)、生物因子(如昆蟲、真菌等)(Moles Westoby,2000)、樹高、林分郁閉度等多重因素的影響,且不同影響因子的影響程度不同。如荷木、格木、水石梓,種內(nèi)差異大于種間差異,每個(gè)物種分別選自三種生境,可推測(cè)環(huán)境因子對(duì)它們的影響最大,其他因子次之;而廣東潤楠、竹節(jié)樹(Caralliabrachiata),種內(nèi)差異明顯小于種間差異,可能樹高、林分郁閉度對(duì)其影響更大一些。由于全部物種海拔范圍36~288 m,受海拔差異影響較小,故本文未考慮海拔對(duì)種內(nèi)種間差異的影響。

    展葉物候與植物的生長策略有關(guān)。依據(jù)優(yōu)先效應(yīng)(Sale,1977)假說,在生長季早期定殖的物種可以優(yōu)先利用資源而獲得生長優(yōu)勢(shì),又因?yàn)槊磕曛参锏恼谷~和葉衰亡時(shí)間與植物對(duì)土壤和光資源的吸收利用密切相關(guān)。所以,展葉物候之間的差異表明,可能展葉較早的個(gè)體在隨后的生長中具備一定的優(yōu)勢(shì),可以提高物種在群落中的適合度。最近,Umana et al(2015)利用熱帶樹苗研究了常見種和稀有種種內(nèi)性狀變異,發(fā)現(xiàn)常見種的生長和性狀變異小,常見種經(jīng)常占據(jù)種群性狀空間的核心位置,能很好的適應(yīng)可利用的環(huán)境條件,稀有種需要發(fā)生大的種內(nèi)性狀變異,才能在有限的環(huán)境資源下生存。常見種在給定環(huán)境下具有最佳的資源利用性狀和具有負(fù)效應(yīng)的偏離最佳值的性狀,也因此降低了種內(nèi)的性狀變異。稀有種大多處于不適宜條件下,為了生存而表現(xiàn)高的表型變異(Umana et al,2015)。那么,物候作為一種功能性狀,是否也具有相似的性狀變異格局是一個(gè)值得探討的問題。當(dāng)然,物種物候的完全一致性不是期望的進(jìn)化結(jié)果(Johansson et al,2015),植物的展葉時(shí)間是一系列權(quán)衡策略選擇的結(jié)果,植物可以通過展葉物候的差異性減少時(shí)間維度上對(duì)光照、水分以及養(yǎng)分的競(jìng)爭(zhēng)利用,調(diào)整群體的生存策略,使得整個(gè)群體的適合度達(dá)到最大(Singer Parmesan,2010),從而在時(shí)間異質(zhì)性中達(dá)到物種共存(Nakashizuka,2001)。是否常見種的種內(nèi)物候變異小于稀有種的種內(nèi)物候變異,還有待進(jìn)一步研究。

    3.3 展葉物候的系統(tǒng)發(fā)育關(guān)系

    本研究結(jié)果表明沒有顯著的譜系信號(hào),即親緣關(guān)系近的物種并沒有選擇在相近的時(shí)間展葉。雖然物種之間的功能性狀存在顯著的譜系信號(hào),但也不是所有群落中的功能性狀都具有顯著的譜系信號(hào)(Losos,2008;Boyle et al,2012;Swenson et al,2012)。這說明親緣關(guān)系相近的物種展葉過程并沒有表現(xiàn)出具有方向性,或者展葉有趨同進(jìn)化傾向??赡茉蛴腥c(diǎn),一是由于所建譜系樹是基于APG Ⅲ物種數(shù)據(jù)庫,只能分辨到科或?qū)伲直媛屎途_度較低(Kress et al,2009),從而對(duì)研究結(jié)果造成影響;二是選取的物種數(shù)量較少,具有一定的局限性;第三,鼎湖山屬于熱帶和亞熱帶地區(qū)交界的位置,具有一定的氣候特性,與臨近山脈相比氣溫較低,降雨量及相對(duì)濕度較大(黃忠良等,1998),所形成的氣候格局可能使物種展葉時(shí)間具有隨機(jī)性。

    4結(jié)論

    對(duì)植物物候的影響因素及種間種內(nèi)差異的研究有助于理解生物多樣性及物種共存機(jī)制(Guisan Thuiller,2005)。本研究發(fā)現(xiàn)在南亞熱帶植物群落中,植物展葉物候與氣候因子(降雨和溫度)相關(guān)性不顯著,展葉時(shí)間種間和種內(nèi)差異很大,并且不具有譜系保守性。不僅氣候、生長環(huán)境等生態(tài)因子對(duì)植物展葉物候具有影響,植物物種的生理特性、植物系統(tǒng)發(fā)育史等因素也對(duì)其有重要影響。未來還需要選取更多的物種及更精確的統(tǒng)計(jì)方法進(jìn)一步探討亞熱帶植物展葉物候變異的調(diào)控機(jī)制,如利用DNA條形碼技術(shù)能提高植物譜系樹的分辨率和精確度(Kress et al,2009)。因此,今后的研究當(dāng)中可考慮多理論綜合研究的方法探索植物展葉物候格局。同時(shí),選取多年的觀測(cè)結(jié)果和更加豐富的物種進(jìn)行研究,探索植物展葉物候的種間和種內(nèi)變異對(duì)氣候變化的影響,可對(duì)未來森林組成成分,甚至生態(tài)系統(tǒng)過程進(jìn)行建立模型預(yù)測(cè)。

    參考文獻(xiàn):

    AGOSTINELLI C, LUND U, 2013. R package ‘circular’: CircularStatistics. https://rforge.rproject.org/projects/circular/.

    APG Ⅲ (Angiosperm Phylogeny Group Ⅲ), 2009. An update of the Angiosperm Phylogeny Group classification for theorders and families of flowering plants: APG Ⅲ [J]. Bot J Linn Soc, 161: 105-121.

    BLOMBERG SP, THEODORE G, IVES AR, 2003. Testing for phylogenetic signal in comparative data: behavioral traits are more labile [J]. Evolution, 57(4): 717-745.

    BOYLE WA, BRONSTEIN JL, 2012. Phenology of tropical understory trees: patterns and correlates [J]. RevistaDe Biología Trop, 60: 1415-1430.

    CHANGYANG CH, LU CL, SUN IF, et al, 2013. Flowering and fruiting patterns in a subtropical rain forest, Taiwan [J]. Biotropica, 45: 165-174.

    CLELAND EE, CHUINE I, MEZEL A, 2007. Shifting plant phenology in response to global change [J]. Trends Ecol Evol, 22: 357-365.

    DAVIES TJ, WOLKOVICH EM, KRAFT NJB, et al, 2013. Phylogenetic conservatism in plant phenology [J]. J Ecol, 101: 1520-1530.

    DU YJ, MAO LF, QUEENBOROUGH SA, et al, 2015. Phylogenetic constraints and trait correlates of flowering phenology in the angiosperm flora of China [J]. Glob Ecol Biogeogr, 24(8): 928-938.

    FRIDLEY JD, 2012. Extended leaf phenology and the autumn niche in deciduous forest invasions [J]. Nature, 485:359-362.

    GUISAN A, THUILLER W, 2005. Predicting species distribution: offering more than simple habitat models [J]. Ecol Lett, 8(9):993-1009.

    HU XL, CHANGYANG JH, MI XC, et al, 2015. Climatic variables, phylogeny, and functional traits correlates of flowering phenology in subtropical evergreen broadleave forest, east China [J]. Biodivers Sci, 23(5):601-609. [胡小麗, 張楊家豪, 米湘成, 等, 2015.浙江古田山亞熱帶常綠闊葉林開花物候:氣候因素、系統(tǒng)發(fā)育關(guān)系和功能性狀的影響 [J]. 生物多樣性, 23(5):601-609.]

    HUANG ZL, MENG ML, ZHANG YC, 1998. Climate of Dinghushan Biosphere Reserve [J]. Trop Subtr For Ecosyst, 8: 134-139. [黃忠良, 蒙滿林, 張佑昌, 1998.鼎湖山生物圈保護(hù)區(qū)的氣候 [J]. 熱帶亞熱帶森林生態(tài)系統(tǒng)研究, 8: 134-139.]

    JIA P, DU GZ, 2011. Relationships between flowering phenology and functional traits in eastern tibet alpine meadow [J]. Arctic Antarctic Alpine Res, 43(4): 585-592.

    JOHANSSON J, KRISTENSEN NP, NILSSON J, et al, 2015. The ecoevolutionary consequences of interspecific phenological asynchrony— a theoretical perspective [J]. Oikos, 124(1): 102-112.

    KEMBEL SW, COWAN PD, HELMU MR, et al, 2010. Picante: R tools for integrating phylogenies and ecology [J]. Bioinformatics, 26: 1463-1464.

    KONG GH, YE WH, HUANG ZL, et al, 1998. Longterm monitoring of the low subtropical evergreen broadleaved forest in Dinghushan Biosphere Reserve. I. Species composition of Castanopsis chinensis, Cryptocarya concinna community and its contribution to regional species pool [J]. Trop Subtrop For Ecosyst, 8: 1-6. [孔國輝, 葉萬輝, 黃忠良, 等, 1998. 鼎湖山南亞熱帶常綠闊葉林定位研究 (I): 錐栗, 黃果厚殼桂群落組成及其對(duì)區(qū)域物種庫的貢獻(xiàn) [J]. 熱帶亞熱帶森林生態(tài)系統(tǒng)研究, 8: 1-6.]

    KRESS WJ, ERICKSON DL, JONES FA, et al, 2009. Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama [J]. Proc Nat Acad Sci, 106: 18621-18626.

    LECHOWICZ M, 1984. Why do temperate deciduous trees leaf out at different timesadaptation and ecology of forest communities [J]. Am Nat, 124:821-842.

    LESSARDTHERRIEN M, DAVIES TJ, BOLMGREN K, 2013. A phylogenetic comparative study of flowering phenology along an elevational gradient in the Canadian subarctic [J]. Int J Biometeorol, 58: 455-462.

    LI MJ, WANG ZH, 1984. Phenology of common plants in Dinghushan [J]. Trop Subtrop For Ecosyst, 2: 1-10. [李明佳, 王鑄豪, 1984.鼎湖山常見植物的物候 [J]. 熱帶亞熱帶森林生態(tài)系統(tǒng)研究, 2: 1-10.]

    LIU Z, LI K, CAI Y, et al, 2011. Correlations between leafing phenology andtraits: woody species of evergreen broadleaved forests in subtropical China [J]. Polish J Ecol, 59: 463-473.

    LOSOS JB, 2008. Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species [J]. Ecol Lett,11: 995-1003.

    MENZEL A, ESTRELLA N, FABIAN P, 2001. Spatial and temporal variability of the phenological seasons in Germany from 1951 to 1996 [J]. Glob Change Biol, 7(6):657-666.

    MOLES A, WESTOBY M, 2000. Do small leaves expand faster than large leaves, and do shorter expansion times reduce herbivore damage [J]. Oikos, 90(3): 517-524.

    MORAN EV, HARTIG F, BELL DM, 2015. Intraspecific trait variation across scales: implications for understanding global change responses [J]. Glob Change Biol, doi:10.1111/gcb.13000.

    NAKASHIZUKA T, 2001. Species coexistence in temperate, mixed deciduous forests [J]. Trends Ecol Evol, 16(4): 205-210.

    PANCHEN ZA, PRIMACK RB, NORDT B, et al, 2014. Leaf out times of temperate woody plants are related to phylogeny, deciduousness, growth habit and wood anatomy [J]. New Phytologist, 203(4): 1208-1219.

    PARMESAN C, YOHE G, 2003. A globally coherent fingerprint of climate change impacts across natural systems [J]. Nature, 421(6918):37-42.

    POLGAR CA, PRIMACK RB, 2011. Leafout phenology of temperate woody plants: from trees to ecosystems [J]. New Phytol, 191: 926-941.

    QIN JD, 1987. Relationships between insects and plants [M]. Beijing: Science Press: 39-40. [欽俊德, 1987.昆蟲與植物的關(guān)系 [M]. 北京: 科學(xué)出版社: 39-40.]

    QUEENBOROUGH SA, MAZER SJ, VAMOSI SM, et al, 2009. Seed mass, abundance and breeding system among tropical forest species: do dioecious species exhibit compensatory reproduction or abundances [J]. J Ecol, 97: 555-566.

    RATHCKE B, LACEY EP, 2003. Phenological patterns of terrestrial plants [J]. Ann Rev Ecol Syst, 16(4): 179-214.

    SALE PF, 1977. Maintenance of high diversity in coral reef fish communities [J]. Am

    Nat, 111(978):337-359.

    SINGER MC, PARMESAN C, 2010. Phenological asynchrony between herbivorous insects and their hosts: signal of climate change or preexisting adaptive strategy [J]. Philosoph Trans Roy Soc Lond, 365(1555): 3161-3176.

    SUN S, JIN D, LI R, 2006. Leaf emergence in relation to leaf traits in temperate woody species in EastChinese Quercusfabri forests [J]. Acta Oecol, 30(2): 212-222.

    SWENSON NG, STEGEN JC, DAVIES SJ, et al, 2012. Temporal turnover in the composition of tropical tree communities: functional determinism and phylogenetic stochasticity [J]. Ecology, 93(3): 490-499.

    UMANA MN, ZHANG C, CAO M, et al, 2015. Commonness, rarity, and intraspecific variation in traits and performance in tropical tree seedlings [J]. Ecol Lett,18(12): 1329-1337.

    VAN VOLKENBURGH E, 1999. Leaf expansion-an integrating plant behavior [J]. Plant Cell Environ, 22: 1463-1473.

    VITASSE Y, PORTE AJ, KREMER A, et al, 2009. Responses of canopy duration to temperature changes in four temperate tree species: relative contributions of spring and autumn leaf phenology [J]. Oecologia, 161: 187-198.

    WAN MW, LI XZ, 1979. Chinese phenology observation methods [M]. Beijing: Science Press: 49-51. [宛敏渭,劉秀珍, 1979.中國物候觀測(cè)方法 [M].北京: 科學(xué)出版社: 49-50.]

    WEBB CO, DONOGHUE MJ, 2005. Phylomatic: tree assembly for applied phylogenetics [J]. Mol Ecol Notes, 5: 181-183.

    WEBB CO, ACKERLY DD, KEMBEL SW, 2008. Phylocom: software for the analysis of phylogenetic community structure and character evolution [J]. Bioinformatics, 24: 2098-2100.

    WIKSTROM N, SAVOLAINEN V, CHASE MW, 2001. Evolution of the angiosperms: calibrating the family tree [J]. Proc Roy Soc Lond B: Biol Sci, 268(1482): 2211-2220.

    XU GX, LUO SX, GUO QS, et al, 2014. Responses of leaf unfolding and flowering to climate change in 12 tropical evergreen broadleaf tree species in Jianfengling, Hainan Island [J]. J Plant Ecol, 38(6): 585-598. [許格希, 羅水興, 郭泉水,等, 2014.海南島尖峰嶺12種熱帶常綠闊葉喬木展葉期與開花期對(duì)氣候變化的響應(yīng) [J]. 植物生態(tài)學(xué)報(bào), 38(6): 585-598.]

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

    ZHAO JB, ZHANG YP, TAN ZH, et al, 2012. Using digital cameras for comparative phenological monitoring in an evergreen broadleaved forest and a seasonal rainforest [J]. Ecol Inform, 10: 65-72.

    ZHU XB, LIU YM, SUN SC, 2005. Leaf expansion of the dominant woody species of three deciduous Oak forests in Nanjing, east China [J]. J Plant Ecol, 29(1): 128-136. [朱旭斌, 劉婭梅, 孫書存, 2005.南京地區(qū)落葉櫟林主要木本植物的展葉動(dòng)態(tài)研究 [J]. 植物生態(tài)學(xué)報(bào), 29(1): 128-136.]

    ZHU XB, SUN SC, 2006. Leaf phenology of woody species in deciduous broadleaved Oak forests in Nanjing area, east China [J]. J Plant Ecol, 30(1): 25-32. [朱旭斌, 孫書存, 2006.南京地區(qū)落葉櫟林木本植物葉物候研究 [J]. 植物生態(tài)學(xué)報(bào), 30(1): 25-32.]

    黄色毛片三级朝国网站| 日本精品一区二区三区蜜桃| 啪啪无遮挡十八禁网站| 国产黄频视频在线观看| 国产人伦9x9x在线观看| 国产精品一区二区精品视频观看| 久久久久久久国产电影| 亚洲成人国产一区在线观看| 最近最新免费中文字幕在线| 中文字幕色久视频| 一本久久精品| 精品人妻1区二区| 人妻一区二区av| 狠狠狠狠99中文字幕| 精品国产乱码久久久久久小说| 国产激情久久老熟女| 欧美乱码精品一区二区三区| 在线亚洲精品国产二区图片欧美| 80岁老熟妇乱子伦牲交| 天天操日日干夜夜撸| 别揉我奶头~嗯~啊~动态视频 | 日本91视频免费播放| 亚洲国产精品999| a级片在线免费高清观看视频| 妹子高潮喷水视频| 不卡av一区二区三区| a 毛片基地| 欧美日韩av久久| 亚洲一区中文字幕在线| 国产成人精品久久二区二区免费| 色视频在线一区二区三区| 亚洲avbb在线观看| 各种免费的搞黄视频| 大香蕉久久成人网| av超薄肉色丝袜交足视频| 大型av网站在线播放| 国产精品自产拍在线观看55亚洲 | 亚洲国产精品成人久久小说| 中文字幕高清在线视频| 脱女人内裤的视频| 免费一级毛片在线播放高清视频 | 中文字幕色久视频| 日本猛色少妇xxxxx猛交久久| 亚洲精品国产色婷婷电影| 91老司机精品| 免费一级毛片在线播放高清视频 | 青春草亚洲视频在线观看| 天堂8中文在线网| 日日夜夜操网爽| 亚洲欧美一区二区三区久久| 精品久久久精品久久久| 一级毛片精品| 黑人欧美特级aaaaaa片| 大片电影免费在线观看免费| www日本在线高清视频| av在线app专区| 19禁男女啪啪无遮挡网站| 桃红色精品国产亚洲av| 国产麻豆69| 欧美一级毛片孕妇| 色综合欧美亚洲国产小说| 熟女少妇亚洲综合色aaa.| 久久久久久亚洲精品国产蜜桃av| 91麻豆av在线| 亚洲欧洲精品一区二区精品久久久| 淫妇啪啪啪对白视频 | 麻豆乱淫一区二区| 欧美 日韩 精品 国产| videos熟女内射| 美女午夜性视频免费| 热99久久久久精品小说推荐| 亚洲精品美女久久久久99蜜臀| 亚洲性夜色夜夜综合| 国产福利在线免费观看视频| 国产成人精品无人区| 日韩大码丰满熟妇| 亚洲欧美日韩高清在线视频 | 黄片小视频在线播放| 男女边摸边吃奶| 国产亚洲午夜精品一区二区久久| 亚洲自偷自拍图片 自拍| 亚洲成人免费电影在线观看| 日韩中文字幕欧美一区二区| 免费少妇av软件| 99国产综合亚洲精品| 久久久精品国产亚洲av高清涩受| 亚洲专区字幕在线| 日韩 欧美 亚洲 中文字幕| 欧美日韩亚洲综合一区二区三区_| 国产成人一区二区三区免费视频网站| 精品国产国语对白av| 91字幕亚洲| 亚洲中文av在线| 欧美精品一区二区免费开放| 国产成人免费无遮挡视频| 丝袜在线中文字幕| 成人亚洲精品一区在线观看| 免费不卡黄色视频| 欧美国产精品va在线观看不卡| 亚洲 国产 在线| 国产精品麻豆人妻色哟哟久久| 国产精品秋霞免费鲁丝片| av超薄肉色丝袜交足视频| 黄片大片在线免费观看| 日本wwww免费看| 免费不卡黄色视频| 日韩大片免费观看网站| 天堂8中文在线网| 欧美精品啪啪一区二区三区 | 日韩视频在线欧美| 热99久久久久精品小说推荐| 一区福利在线观看| 蜜桃在线观看..| 青春草视频在线免费观看| 久久亚洲精品不卡| 日韩 欧美 亚洲 中文字幕| 欧美午夜高清在线| 午夜免费观看性视频| 人妻 亚洲 视频| 日韩,欧美,国产一区二区三区| 欧美激情极品国产一区二区三区| 精品人妻熟女毛片av久久网站| 美女高潮喷水抽搐中文字幕| 国产真人三级小视频在线观看| 亚洲av美国av| 男人舔女人的私密视频| 精品一区二区三卡| 亚洲欧美色中文字幕在线| 国产主播在线观看一区二区| 99精品欧美一区二区三区四区| 亚洲七黄色美女视频| 99香蕉大伊视频| 欧美另类亚洲清纯唯美| av国产精品久久久久影院| 男人爽女人下面视频在线观看| 啦啦啦啦在线视频资源| videos熟女内射| 日韩欧美一区二区三区在线观看 | cao死你这个sao货| 高清黄色对白视频在线免费看| 狠狠婷婷综合久久久久久88av| 国产精品九九99| 人妻一区二区av| 日本a在线网址| 久久这里只有精品19| 国产亚洲午夜精品一区二区久久| 久久 成人 亚洲| 一边摸一边做爽爽视频免费| 人人妻,人人澡人人爽秒播| 亚洲精品一二三| 日韩欧美免费精品| 两性夫妻黄色片| 精品欧美一区二区三区在线| 国产1区2区3区精品| 国产视频一区二区在线看| 一个人免费看片子| 国产精品麻豆人妻色哟哟久久| 久久亚洲国产成人精品v| 精品久久久久久久毛片微露脸 | 少妇人妻久久综合中文| 一级片'在线观看视频| 亚洲国产欧美网| 亚洲国产欧美日韩在线播放| 国产xxxxx性猛交| 精品久久久久久久毛片微露脸 | 国产成人a∨麻豆精品| 国产精品99久久99久久久不卡| 午夜免费成人在线视频| 日韩三级视频一区二区三区| 亚洲av国产av综合av卡| 国产激情久久老熟女| av不卡在线播放| 精品人妻1区二区| 一二三四社区在线视频社区8| 国产精品 国内视频| 99热网站在线观看| 亚洲精品av麻豆狂野| 精品少妇一区二区三区视频日本电影| 日韩欧美国产一区二区入口| 久久av网站| 欧美日韩亚洲综合一区二区三区_| 亚洲美女黄色视频免费看| 夫妻午夜视频| 侵犯人妻中文字幕一二三四区| 久久天堂一区二区三区四区| 欧美日韩一级在线毛片| 老司机深夜福利视频在线观看 | 亚洲国产精品一区三区| 999久久久国产精品视频| 中文字幕人妻丝袜制服| 欧美另类一区| 亚洲精品国产区一区二| 狠狠婷婷综合久久久久久88av| 2018国产大陆天天弄谢| 成人av一区二区三区在线看 | 侵犯人妻中文字幕一二三四区| 久久狼人影院| 在线天堂中文资源库| 1024视频免费在线观看| av网站在线播放免费| 国产精品99久久99久久久不卡| 亚洲国产精品999| 免费人妻精品一区二区三区视频| www.精华液| 丁香六月欧美| 欧美日韩福利视频一区二区| 人人妻人人澡人人看| 黑人猛操日本美女一级片| 中文字幕精品免费在线观看视频| 一级黄色大片毛片| 亚洲精品一区蜜桃| 色94色欧美一区二区| 一级片免费观看大全| 精品国产乱码久久久久久小说| www.熟女人妻精品国产| 午夜福利视频精品| 悠悠久久av| 午夜免费观看性视频| 中文字幕另类日韩欧美亚洲嫩草| 12—13女人毛片做爰片一| 青春草视频在线免费观看| 免费观看人在逋| 精品久久久久久电影网| 久久影院123| 亚洲va日本ⅴa欧美va伊人久久 | 一区福利在线观看| 狠狠婷婷综合久久久久久88av| 久久久精品国产亚洲av高清涩受| av在线播放精品| 欧美成人午夜精品| 啦啦啦免费观看视频1| 国产成人精品在线电影| 国产欧美日韩一区二区精品| 欧美精品av麻豆av| 久久精品国产综合久久久| 一区二区av电影网| 亚洲国产av新网站| 国产高清videossex| 国产成人a∨麻豆精品| 9热在线视频观看99| 新久久久久国产一级毛片| 人人妻人人澡人人爽人人夜夜| 国产精品偷伦视频观看了| 一级毛片女人18水好多| 999精品在线视频| 一区二区日韩欧美中文字幕| 肉色欧美久久久久久久蜜桃| 69精品国产乱码久久久| 伦理电影免费视频| 国产老妇伦熟女老妇高清| 人人妻,人人澡人人爽秒播| 一二三四社区在线视频社区8| 视频在线观看一区二区三区| 国产免费一区二区三区四区乱码| 日韩一区二区三区影片| 精品一品国产午夜福利视频| 性色av乱码一区二区三区2| 桃花免费在线播放| 亚洲av片天天在线观看| 乱人伦中国视频| 中文字幕人妻丝袜一区二区| 最近最新免费中文字幕在线| 午夜激情久久久久久久| 久久精品亚洲熟妇少妇任你| 国产野战对白在线观看| 91成年电影在线观看| 国产欧美日韩综合在线一区二区| 免费在线观看日本一区| 欧美日韩亚洲高清精品| 午夜福利,免费看| 超色免费av| 男女国产视频网站| 手机成人av网站| 欧美97在线视频| 最近最新中文字幕大全免费视频| 成年女人毛片免费观看观看9 | 91麻豆精品激情在线观看国产 | 淫妇啪啪啪对白视频 | 免费观看人在逋| 麻豆国产av国片精品| 国产精品.久久久| 亚洲人成电影观看| 国产片内射在线| 777久久人妻少妇嫩草av网站| 日本欧美视频一区| 国产男女超爽视频在线观看| 麻豆国产av国片精品| 黑人猛操日本美女一级片| 亚洲午夜精品一区,二区,三区| 极品少妇高潮喷水抽搐| 亚洲国产日韩一区二区| 亚洲激情五月婷婷啪啪| 亚洲第一欧美日韩一区二区三区 | 捣出白浆h1v1| 国产野战对白在线观看| 久久久国产精品麻豆| 99久久国产精品久久久| 黄色毛片三级朝国网站| 免费在线观看视频国产中文字幕亚洲 | 视频在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 人妻久久中文字幕网| 午夜日韩欧美国产| 欧美日韩黄片免| 一本一本久久a久久精品综合妖精| 国产成人精品久久二区二区91| 美女高潮喷水抽搐中文字幕| 这个男人来自地球电影免费观看| 老司机福利观看| 淫妇啪啪啪对白视频 | 午夜久久久在线观看| 一区二区三区乱码不卡18| 精品亚洲成a人片在线观看| 国产成人精品在线电影| 亚洲精品一卡2卡三卡4卡5卡 | videosex国产| 叶爱在线成人免费视频播放| 高清视频免费观看一区二区| 最近最新中文字幕大全免费视频| 夜夜夜夜夜久久久久| 咕卡用的链子| 国产成人av教育| 法律面前人人平等表现在哪些方面 | 国产有黄有色有爽视频| 一级片免费观看大全| 亚洲欧洲精品一区二区精品久久久| 午夜精品国产一区二区电影| 国产黄频视频在线观看| 亚洲精品久久午夜乱码| 精品一区二区三区av网在线观看 | 国产在线免费精品| 亚洲人成电影免费在线| 天天躁日日躁夜夜躁夜夜| 国产野战对白在线观看| 亚洲黑人精品在线| 91老司机精品| 欧美国产精品va在线观看不卡| 日本91视频免费播放| 免费观看av网站的网址| 777米奇影视久久| 在线 av 中文字幕| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| 老司机福利观看| 在线 av 中文字幕| 在线永久观看黄色视频| 成人手机av| 啦啦啦中文免费视频观看日本| 黄网站色视频无遮挡免费观看| 蜜桃在线观看..| 欧美精品人与动牲交sv欧美| 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 日韩欧美一区二区三区在线观看 | 久久精品国产综合久久久| videos熟女内射| 男女之事视频高清在线观看| 亚洲综合色网址| 制服人妻中文乱码| 亚洲精品一卡2卡三卡4卡5卡 | 男女边摸边吃奶| 久久影院123| 男人舔女人的私密视频| 狠狠婷婷综合久久久久久88av| 久久久国产一区二区| av免费在线观看网站| 精品少妇黑人巨大在线播放| 精品熟女少妇八av免费久了| 色播在线永久视频| 国产精品麻豆人妻色哟哟久久| 国产成人精品无人区| √禁漫天堂资源中文www| 国产精品国产三级国产专区5o| 一二三四社区在线视频社区8| 三级毛片av免费| 国产在视频线精品| 欧美黑人精品巨大| 欧美成狂野欧美在线观看| 成人影院久久| 婷婷色av中文字幕| 午夜福利免费观看在线| 中文字幕色久视频| av免费在线观看网站| 十八禁网站网址无遮挡| 视频区图区小说| 亚洲精品国产av成人精品| 午夜免费观看性视频| 如日韩欧美国产精品一区二区三区| 啦啦啦视频在线资源免费观看| 免费黄频网站在线观看国产| 国产成人精品无人区| 国产成人av激情在线播放| 18禁观看日本| 高清黄色对白视频在线免费看| av天堂久久9| 亚洲全国av大片| 欧美xxⅹ黑人| 国产日韩一区二区三区精品不卡| 亚洲精品第二区| 捣出白浆h1v1| 狂野欧美激情性xxxx| 久久久久国产精品人妻一区二区| www.999成人在线观看| 最黄视频免费看| 五月天丁香电影| 黄片小视频在线播放| 免费不卡黄色视频| 国产精品一区二区精品视频观看| 国产精品二区激情视频| 最黄视频免费看| 精品久久久精品久久久| av视频免费观看在线观看| 91国产中文字幕| 亚洲av成人一区二区三| 欧美在线黄色| 免费高清在线观看视频在线观看| 精品高清国产在线一区| avwww免费| 男人操女人黄网站| 国产免费一区二区三区四区乱码| 亚洲精品久久成人aⅴ小说| 久久香蕉激情| 欧美性长视频在线观看| 久久久精品免费免费高清| 国产欧美日韩综合在线一区二区| 亚洲专区字幕在线| 亚洲,欧美精品.| 黄色片一级片一级黄色片| 不卡av一区二区三区| 免费在线观看完整版高清| 中文字幕色久视频| 婷婷成人精品国产| 免费观看人在逋| 国产精品一区二区免费欧美 | 国产区一区二久久| 亚洲精品第二区| 少妇裸体淫交视频免费看高清 | 亚洲国产欧美一区二区综合| 国产99久久九九免费精品| 日本91视频免费播放| 中国美女看黄片| 日日夜夜操网爽| 国产一区有黄有色的免费视频| 18在线观看网站| 少妇猛男粗大的猛烈进出视频| 美女中出高潮动态图| 69精品国产乱码久久久| 国产老妇伦熟女老妇高清| 欧美黄色片欧美黄色片| 亚洲国产av新网站| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品亚洲av一区麻豆| 水蜜桃什么品种好| 黑人巨大精品欧美一区二区蜜桃| 日韩视频一区二区在线观看| 国产激情久久老熟女| 99久久人妻综合| 一级片免费观看大全| 久久99热这里只频精品6学生| 天天躁狠狠躁夜夜躁狠狠躁| 老鸭窝网址在线观看| 亚洲成人手机| 久久毛片免费看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 欧美一区二区三区| 巨乳人妻的诱惑在线观看| 国产欧美日韩一区二区精品| 在线 av 中文字幕| 在线天堂中文资源库| 国产精品 欧美亚洲| 国产在线一区二区三区精| av网站免费在线观看视频| 男男h啪啪无遮挡| 日本91视频免费播放| 中国美女看黄片| 亚洲精品国产av成人精品| 亚洲久久久国产精品| 精品久久蜜臀av无| 国产一级毛片在线| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 久热这里只有精品99| 亚洲国产av新网站| 秋霞在线观看毛片| 欧美变态另类bdsm刘玥| 中文欧美无线码| 亚洲av片天天在线观看| 老熟女久久久| av超薄肉色丝袜交足视频| 99久久综合免费| 两个人免费观看高清视频| 亚洲欧美成人综合另类久久久| 一级a爱视频在线免费观看| 美女高潮到喷水免费观看| 夫妻午夜视频| 日本五十路高清| 亚洲精品久久久久久婷婷小说| av电影中文网址| 91精品三级在线观看| 精品欧美一区二区三区在线| 狠狠精品人妻久久久久久综合| 色精品久久人妻99蜜桃| 亚洲欧美一区二区三区久久| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看 | 人人妻,人人澡人人爽秒播| 国产成人一区二区三区免费视频网站| 乱人伦中国视频| e午夜精品久久久久久久| 女人高潮潮喷娇喘18禁视频| 首页视频小说图片口味搜索| 欧美97在线视频| 午夜成年电影在线免费观看| 少妇被粗大的猛进出69影院| 久久国产精品人妻蜜桃| 亚洲成av片中文字幕在线观看| 脱女人内裤的视频| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 国产91精品成人一区二区三区 | 精品欧美一区二区三区在线| 婷婷色av中文字幕| 国产精品 国内视频| 欧美日韩精品网址| 欧美亚洲 丝袜 人妻 在线| 蜜桃国产av成人99| 亚洲天堂av无毛| www.自偷自拍.com| 少妇的丰满在线观看| 啦啦啦免费观看视频1| 欧美国产精品一级二级三级| 亚洲成人国产一区在线观看| 青青草视频在线视频观看| 国产在线免费精品| 在线观看舔阴道视频| 狠狠狠狠99中文字幕| 日韩 亚洲 欧美在线| 欧美日韩成人在线一区二区| 三级毛片av免费| 少妇精品久久久久久久| 国产精品免费大片| 国产欧美日韩一区二区三 | 国产av一区二区精品久久| 精品第一国产精品| 午夜福利免费观看在线| 一本大道久久a久久精品| 黄频高清免费视频| 精品少妇一区二区三区视频日本电影| 久久久国产成人免费| 国产亚洲精品第一综合不卡| www.精华液| 久久中文看片网| 国产成人a∨麻豆精品| 欧美在线一区亚洲| 欧美日韩亚洲高清精品| 亚洲第一欧美日韩一区二区三区 | 免费高清在线观看日韩| 亚洲午夜精品一区,二区,三区| 大香蕉久久成人网| av天堂在线播放| 中文字幕色久视频| 宅男免费午夜| 日日爽夜夜爽网站| 少妇裸体淫交视频免费看高清 | 国产男女超爽视频在线观看| av不卡在线播放| 久久久久国内视频| 亚洲专区中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 男女无遮挡免费网站观看| 国产高清国产精品国产三级| 成人免费观看视频高清| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人精品久久二区二区免费| 亚洲七黄色美女视频| 一级a爱视频在线免费观看| 国产av又大| 亚洲av成人不卡在线观看播放网 | 欧美黄色片欧美黄色片| 99久久99久久久精品蜜桃| 免费看十八禁软件| 91字幕亚洲| 99香蕉大伊视频| 久久久久久久国产电影| 日韩制服骚丝袜av| 欧美另类一区| 国产日韩欧美在线精品| 国产97色在线日韩免费| 在线观看人妻少妇| e午夜精品久久久久久久| 亚洲精品一区蜜桃| 欧美日韩福利视频一区二区| 悠悠久久av| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 美女主播在线视频| 狠狠精品人妻久久久久久综合| 精品国产乱码久久久久久小说| 12—13女人毛片做爰片一| 欧美日韩视频精品一区| 老司机在亚洲福利影院| 交换朋友夫妻互换小说| 九色亚洲精品在线播放| 电影成人av| 国产视频一区二区在线看| 亚洲五月婷婷丁香| 黄网站色视频无遮挡免费观看| 久久精品亚洲av国产电影网| 日韩中文字幕欧美一区二区| 中亚洲国语对白在线视频| 久久久久国产精品人妻一区二区| 亚洲精品日韩在线中文字幕| 国产精品免费视频内射| 免费看十八禁软件| 黄网站色视频无遮挡免费观看|