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

    吉林蛟河針闊混交林樹木生長與生境的關聯(lián)性

    2017-06-22 13:47:48郝珉輝張忠輝趙珊珊趙秀海葉爾江拜克吐爾漢張春雨
    生態(tài)學報 2017年10期
    關鍵詞:關聯(lián)性生境樣地

    郝珉輝, 張忠輝, 趙珊珊, 趙秀海, 葉爾江·拜克吐爾漢, 張春雨,*

    1 北京林業(yè)大學森林資源與生態(tài)系統(tǒng)過程北京市重點實驗室, 北京 100083 2 吉林省林業(yè)科學研究院, 長春 130033 3 新疆農(nóng)業(yè)大學林學與園藝學院, 烏魯木齊 830052

    ?

    吉林蛟河針闊混交林樹木生長與生境的關聯(lián)性

    郝珉輝1, 張忠輝2, 趙珊珊2, 趙秀海1, 葉爾江·拜克吐爾漢3, 張春雨1,*

    1 北京林業(yè)大學森林資源與生態(tài)系統(tǒng)過程北京市重點實驗室, 北京 100083 2 吉林省林業(yè)科學研究院, 長春 130033 3 新疆農(nóng)業(yè)大學林學與園藝學院, 烏魯木齊 830052

    以吉林蛟河21.12 hm2(660 m×320 m)針闊混交林樣地為對象,在以海拔、坡度、坡向以及凸凹度為地形變量劃分樣地生境類型的基礎上,利用2009—2014年植被生長數(shù)據(jù),研究生境差異對樹木生長的影響機制。研究采用生境空間隨機(生境CSR)、物種空間隨機(物種CSR)以及物種Thomas3種生態(tài)學過程零模型,檢驗樹木徑向生長與4種生境類型之間的關聯(lián)性;采用Pearson相關系數(shù)計算樹木生長與地形變量之間的相關關系,量化地形變量對樹木生長的影響。研究結(jié)果顯示:(1)樣地內(nèi)絕大多數(shù)物種(生境CSR過程:86.1%;物種CSR過程:94.4%;物種Thomas過程:61.1%)的徑向生長表現(xiàn)出明顯的生境關聯(lián)性。(2)不同生活型物種的生境利用方式不同:灌木和亞喬木具有類似的生境偏好,與海拔相對較低、地勢相對平緩的生境型1正關聯(lián),與海拔相對較高、坡度相對較大的生境型2和生境型4負關聯(lián);喬木在不同生境類型中會同時存在正、負關聯(lián)性,但更傾向于和生境1發(fā)生負關聯(lián),而和生境4發(fā)生正關聯(lián),并且喬木樹種的徑向生長對地形之間的差異更加敏感;(3)高達86.1%的物種生長與至少一種地形變量顯著相關,其中海拔對樹木生長影響最大,其次是坡度和凸凹度,坡向的影響則相對最小。上述結(jié)果表明,樣地內(nèi)不同物種之間出現(xiàn)了明顯的生境利用性分化,生境利用性分化是影響溫帶針闊混交林樹木徑向生長的重要因素。

    徑向生長; 生境; 地形; 生境差異; 生活型

    森林植物的生長和分布具有一定的生境偏好,植物在其適宜生境中擁有較強的競爭能力,株數(shù)較多且生長較好[1- 2]。通過量化植物的生長、分布與生境之間的關聯(lián)程度,能夠檢驗生境差異在群落構建、物種共存以及樹木生長等過程中的重要性,推斷生態(tài)學現(xiàn)象背后潛在的生態(tài)學過程[3- 6]。熱帶雨林[3- 4]、亞熱帶常綠闊葉林[5]和溫帶森林[6]開展的物種-生境關聯(lián)性研究表明,生境差異是影響物種空間分布的一個重要因素,物種空間分布與生境之間的關聯(lián)程度取決于生境的復雜性[3- 6],生境條件越復雜,物種選擇其最適生境的機會就越大[1- 2]。同時研究表明,不同物種對生境差異的敏感性不同[7],物種對生境差異的響應在很大程度上是通過植物所具有的某些關鍵性狀起作用的[8]。以上研究揭示了生境差異對樹種空間分布的影響及其作用機制,但以上研究只關注樹木的相對多度而忽視樹木的生長特征。樹木的相對多度受多種因素影響[9],僅僅依靠相對多度很難全面反映物種與生境之間的關聯(lián)性。相比之下,樹木生長對環(huán)境間的差異更加敏感[10-11],能夠更加直觀地反映物種對生境的偏好。因此,有必要探討樹木生長與生境之間的關聯(lián)性,研究局域尺度上生境差異對樹木生長的影響,進一步拓展人們對物種-生境關聯(lián)性的理解和認識。

    在進行物種-生境關聯(lián)性研究時,由于樹木的種內(nèi)和種間相互作用[12- 13]、自然干擾[14]和生境差異[15]等生態(tài)學過程相互交織,常常導致樹木生長表現(xiàn)出一定程度上的空間自相關[16]。物種-生境關聯(lián)性很容易與這種空間自相關現(xiàn)象發(fā)生混淆,需要通過隨機模擬技術消除空間自相關性對檢驗結(jié)果帶來的負面影響[3]。本文以吉林蛟河固定監(jiān)測樣地2009—2014年的植被生長數(shù)據(jù)為基礎,采用生境空間隨機(生境CSR)、物種空間隨機(物種CSR)和物種Thomas3種生態(tài)學過程作為零模型,在以地形變量劃分生境類型的基礎上,檢驗樹木的徑向生長與各類生境之間的關聯(lián)性,旨在揭示生境差異對溫帶森林樹木生長的影響機制。

    1 研究方法

    1.1 樣地概況

    研究區(qū)域位于吉林省蛟河林業(yè)試驗區(qū)管理局(43°57.897′—43°58.263′N, 127°42.789′—127°43.310′E)。該地區(qū)屬于受季風影響的溫帶大陸性氣候,年均溫為3.8℃,其中7月份為最熱月,平均溫度為21.7℃;1月份為最冷月,平均溫度為- 18.6℃;年平均降水量為695.9 mm;土壤類型為棕色森林土,土層深度20—100 cm。

    2009年建立面積為21.12 hm2(660 m×320 m)固定監(jiān)測樣地,樣地最低海拔425.3 m,最高海拔525.8 m,最大高差100.5 m,地形變化較大。將樣地劃分為528個20 m×20 m樣方,記錄所有胸徑≥1 cm的木本植物的物種、胸徑、樹高、冠幅和相對位置坐標,并掛牌標記以便長期監(jiān)測。2014年對樣地內(nèi)所有胸徑≥1 cm木本植物進行復測。林分類型屬于典型次生針闊混交林,為闊葉紅松林經(jīng)采伐干擾后形成的近熟林。樣地內(nèi)針葉樹種主要為紅松(Pinuskoraiensis),物種重要值(計算參照方精云等[17])位于前五的闊葉樹種依次為:胡桃楸(Juglansmandshurica)、色木槭(Acermono)、千金榆(Carpinuscordata)、暴馬丁香(Syringareticulatavar.amurensis)和紫椴(Tiliaamurensis)。為了滿足樣本量需求,本文以樣地內(nèi)物種多度大于30株的36個樹種作為研究對象。

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

    1.2.1 地形變量與生境類型劃分

    本文依據(jù)地形變量劃分生境類型,地形變量通過影響植物對水分、光照以及礦質(zhì)元素的獲取,從而強烈作用于植物的生長和分布[18- 21],在進行物種-生境關聯(lián)性研究時,地形變量往往作為最主要的生境類型劃分依據(jù)[3- 6]。樣地內(nèi)每個20 m×20 m樣方的地形變量計算方法為:樣方4個頂點的海拔平均值為該樣方的海拔值[3];目標樣方的海拔減去與之相鄰的8個樣方的海拔平均值為該樣方的凸凹度,處于樣地邊緣的樣方凸凹度為樣方中心的海拔減去其4個頂點的海拔平均值[22];從樣方的4個頂點中隨機抽取3個頂點組成一個平面,這樣4個頂點可以組成4個平面,這4個平面與樣方投影面的夾角平均值為該樣方的坡度;4個平面與正北方向的角度平均值為該樣方的坡向[3]。Zhang等[6]在上述地形變量的基礎上,利用多元回歸樹(Multivariate Regression Tree)方法[23- 24],將該研究樣地劃分為4種生境類型(以下稱生境1、生境2、生境3、生境4;圖1)。本文采用了他的生境分類結(jié)果,不同生境類型的特征為:生境1屬于地勢平緩的低海拔區(qū)域,最高海拔小于453.63 m;生境2、生境3屬于由低海拔到高海拔的過渡區(qū),其中生境2位于東坡(坡向≥186°),且坡度較陡;生境3位于西坡(坡向<186°),且坡度較緩;生境4屬于高海拔區(qū)域,其最低海拔不小于465.72 m。

    圖1 吉林蛟河21.12hm2樣地內(nèi)生境類型分布Fig.1 The distribution of the habitat types in a 21.12hm2 research plot in Jiaohe, Jinlin Province每個樣方大小為20 m20 m, 數(shù)字表示生境類型

    1.2.2 物種-生境關聯(lián)性檢驗

    經(jīng)典統(tǒng)計學方法檢驗物種-生境關聯(lián)性時,通常假設樣本之間相互獨立,不考慮空間自相關性的影響。由于各種生態(tài)學過程交互作用,樣本間通常具有明顯的空間自相關性,因此檢驗結(jié)果會存在一定偏差[3]。本文在綜合考慮植物生長和樣方生境空間自相關性的基礎上,采用生境空間隨機(生境CSR)、物種空間隨機(物種CSR)以及物種Thomas3種生態(tài)學過程零模型,檢驗植物徑向生長與生境之間的關聯(lián)性。

    生境CSR過程假設植物生長存在空間自相關性而不同樣方的生境類型之間不存在自相關性。生境CSR過程用于生成生境隨機模擬圖,該過程隨機變換不同樣方的生境類型,而每一種生境類型的總數(shù)量保持不變。物種CSR過程假設樣方生境類型之間存在空間自相關性,而植物生長不存在自相關性。物種CSR過程用于生成植物生長隨機模擬圖,該過程將植株的空間位置進行隨機化,而樣地內(nèi)總生長量大小保持不變。物種Thomas過程假設植物生長與樣方生境類型之間空間相關,用于生成具有空間自相關性的植物生長模擬圖。物種Thomas過程是一種限制性隨機化過程,在很大程度上保留了原有數(shù)據(jù)結(jié)構,因此本文傾向于采用該過程進行結(jié)果解釋。

    在上述零模型基礎上通過1000次模擬建立置信區(qū)間,比較每個物種徑向生長量的真實值與模擬值:當真實值>97.5%的模擬值時,則認為物種徑向生長與該生境類型顯著正相關;當真實值<2.5%的模擬值時,認為物種徑向生長與該生境類型顯著負相關。

    本文采用Pearson相關系數(shù)計算植物生長量與地形變量(海拔、坡度、坡向和凸凹度)之間的關系,量化地形變量對植物生長的影響。

    2 結(jié)果

    利用生境CSR過程作為零模型檢驗植物生長與生境關系,絕大多數(shù)(86.1%)物種的徑向生長與生境顯著相關;而大果榆(U.macrocarpa)、沙松、簇毛槭、毛榛(Corylusmandshurica)和長白忍冬(L.ruprechtiana)的徑向生長與生境沒有檢測到顯著關聯(lián)性;在36個物種和4種生境構成的144個潛在組合中,33個物種-生境組合呈顯著正關聯(lián),44個呈顯著負關聯(lián)。利用物種CSR過程作為零模型,94.4%的物種徑向生長與生境顯著相關;僅有大果榆和簇毛槭的徑向生長與生境沒有檢測到顯著關聯(lián)性;38個物種-生境組合呈顯著正關聯(lián),55個呈顯著負關聯(lián)。利用物種Thomas過程作為零模型,22個物種(占總物種數(shù)61.1%)的徑向生長與生境顯著相關;21個物種-生境組合呈顯著正關聯(lián),24個呈顯著負關聯(lián)(表1,表2)。

    表1 樹種徑向生長與生境的關聯(lián)性

    CSR: 完全空間隨機, Complete Spatial Randomness; H:生境; +: 正關聯(lián), -: 負關聯(lián), N: 中性

    表2 樹種徑向生長與生境顯著關聯(lián)性統(tǒng)計

    +:正相關, -:在負相關, N:中性

    物種Thomas過程檢驗到22個物種與生境顯著關聯(lián),其中灌木和亞喬木徑向生長具有類似的生境偏好,喬木則與它們存在明顯差異:灌木和亞喬木與生境1正關聯(lián),與生境2和生境4負關聯(lián);喬木在不同生境中同時存在正、負關聯(lián)性,但在生境1中負關聯(lián)占優(yōu)勢,生境4中正關聯(lián)占優(yōu)勢(圖2)。

    圖2 生活型對物種-生境關聯(lián)性的影響Fig.2 Effect of different life forms on the associations between species and habitats

    利用Pearson相關系數(shù)檢驗生長量與地形變量之間的關系。除大果榆、沙松、簇毛槭、翅衛(wèi)矛(Euonymusmacropterus)和東北山梅花之外,其余物種都至少與一種地形變量顯著相關,占物種總數(shù)的86.1%。25個物種的生長與海拔顯著相關(正相關10個,負相關15個);22個物種的生長與坡度顯著相關(正相關10個,負相關12個);20個物種的生長與凸凹度顯著相關(正相關11個,負相關9個);13個物種的生長與坡向顯著相關(正相關6個,負相關7個)。

    3 討論

    利用物種Thomas過程檢驗植物生長與生境關聯(lián)性更加保守和可靠?;谠撨^程,群落中超過60%的物種表現(xiàn)出明顯的生境關聯(lián)性,說明生境差異是影響溫帶針闊混交林中樹木徑向生長的重要因素。樣地內(nèi)大約70%的喬木樹種與生境顯著關聯(lián),而灌木樹種卻不足40%,說明喬木的徑向生長對地形狀況更加敏感。相對于地形差異,灌木的生長則更容易受到林冠結(jié)構、林下光照條件等因素的影響[25]。本文顯示亞喬木和灌木的生境偏好較為類似,而與喬木存在明顯不同,反映了不同生活型的樹種在資源利用上的差異性和互補性:在地形條件優(yōu)越的生境中,樹木生長所必須的水分、光照以及礦質(zhì)元素等資源相對充足,喬木生長迅速,喬木的生長狀況會對林下微環(huán)境產(chǎn)生顯著影響[26],在喬木生長良好的生境中,林分郁閉度高、林下光照減弱,光因子成為林下植物的主要限制因子,因此灌木和亞喬木生長相對緩慢;相反,在地形條件惡劣的生境中,由于資源相對貧乏,喬木的生長受到限制、林冠開闊,灌木和亞喬木能夠獲得充足的光照,反而容易實現(xiàn)快速生長[27]。

    本研究中生境1處于低海拔區(qū)域,地勢平緩、水分充足,但由于歷史采伐干擾較重,先鋒樹種較多、林下光照相對充足。在這種情況下喬木與生境1的關聯(lián)性表現(xiàn)為既存在正關聯(lián)也存在負關聯(lián)但以負關聯(lián)為主,而灌木和亞喬木與生境1的關聯(lián)性則主要表現(xiàn)為正關聯(lián)或者中性。生境2和生境3屬于由低海拔到高海拔的過渡區(qū)域,具有一定的坡度,坡度會造成到土壤中水分和礦質(zhì)元素的流失[28- 29],受此影響,這兩類生境中樹木的生長很少表現(xiàn)為正關聯(lián)。生境2和生境3的區(qū)別在于:生境2屬于半陰坡,生境3屬于半陽坡,生境3的光照條件優(yōu)于生境2;同時生境2的坡度更大,林分郁閉度通常會隨著坡度的增加而增大[29],導致生境2中灌木和亞喬木可獲得的光照進一步減少。因此,在生境2中樹木生長表現(xiàn)出負相關的情況較為普遍,其中的灌木和亞喬木生長尤其明顯;但在生境3中樹木的生長則全部表現(xiàn)為中性。生境4處于高海拔區(qū)域,同樣具有一定的坡度,屬于半陽坡,光照條件良好,有利于喬木樹種的生長,而灌木和亞喬木則主要呈現(xiàn)為負關聯(lián)或中性。

    譜系保守性假說認為親緣關系較近的物種擁有較為相似的性狀特征,因此對生境的需求也比較相似[30- 31],Gunatilleke等[32]認為生長在同一區(qū)域的近緣物種在耐陰性、養(yǎng)分利用等方面存在差異,從而確保近緣物種共存。熱帶雨林的相關研究表明生長在同一區(qū)域的近緣物種之間往往表現(xiàn)出不同的生境關聯(lián)性[32- 33]。本文包含兩個以上物種的屬共有9個,除樺木屬的白樺(Betulaplatyphylla)和楓樺(B.costata)均未表現(xiàn)出顯著生境偏好外,其余8個屬的同屬物種之間生境關聯(lián)性差異明顯。例如,水曲柳(Fraxinusmandshurica)和花曲柳(F.rhynchophylla)同為木犀科梣屬植物,但水曲柳的生長與生境1正關聯(lián),而花曲柳的生長與生境4正關聯(lián);春榆、裂葉榆和大果榆同為榆科榆屬植物,但春榆的生長與生境1正關聯(lián)、裂葉榆的生長與生境4正關聯(lián),而大果榆的生長則沒有表現(xiàn)出顯著的生境偏好。以上結(jié)果表明,溫帶森林近緣物種之間同樣產(chǎn)生了生境利用性分化,以滿足物種生長和共存的生態(tài)需求。物種對環(huán)境條件的利用依賴于物種所具有的某些關鍵性狀,同時,物種的性狀組成往往又決定了物種的生長狀況[34- 35]。因此,更加深入的物種-生境關聯(lián)性研究,需要結(jié)合樹木的性狀組成以及系統(tǒng)發(fā)育特征,檢驗樹木的功能性狀對生境差異的響應,進一步深刻揭示境森林群落中樹木生長與生境之間的關聯(lián)機制。

    [1] Hubbell S P, Foster R B. Diversity of canopy trees in a neotropical forest and implications for conservation. Special Publications Series of the British Ecological Society, 1983, 2: 25- 41.

    [2] Tilman D, Pacala S. The maintenance of species richness in plant communities//Ricklefs R E, Schulter D, eds. Species Diversity in Ecological Communities. Chicago: University of Chicago Press, 1993: 13- 25.

    [3] Harms K E, Condit R, Hubbell S P, Foster R B. Habitat associations of trees and shrubs in a 50-ha neotropical forest plot. Journal of Ecology, 2001, 89(6): 947- 959.

    [4] Itoh A, Ohkubo T, Nanami S, Tan S, Yamakura T. Comparison of statistical tests for habitat associations in tropical forests: a case study of sympatric dipterocarp trees in a Bornean forest. Forest ecology and management, 2010, 259(3): 323- 332.

    [5] Lai J S, Mi X C, Ren H B, Ma K P. Species-habitat associations change in a subtropical forest of China. Journal of Vegetation Science, 2009, 20(3): 415- 423.

    [6] Zhang C Y, Zhao Y Z, Zhao X H, Von Gadow K. Species-habitat associations in a northern temperate forest in China. Silva Fennica, 2012, 46(4): 501- 519.

    [7] Comita L S, Engelbrecht B M J. Seasonal and spatial variation in water availability drive habitat associations in a tropical forest. Ecology, 2009, 90(10): 2755- 2765.

    [8] Liu J J, Tan Y H, Slik J W F. Topography related habitat associations of tree species traits, composition and diversity in a Chinese tropical forest. Forest Ecology and Management, 2014, 330: 75- 81.

    [9] Newbery D M, Campbell E J F, Proctor J, Still M J. Primary lowland dipterocarp forest at Danum Valley, Sabah, Malaysia. Species composition and patterns in the understorey. Vegetatio, 1996, 122(2): 193- 220.

    [10] Russo S E, Brown P, Tan S, Davies S J. Interspecific demographic trade-offs and soil-related habitat associations of tree species along resource gradients. Journal of Ecology, 2008, 96(1): 192- 203.

    [11] Benito-Garzón M, Ruiz-Benito P, Zavala M A. Interspecific differences in tree growth and mortality responses to environmental drivers determine potential species distributional limits in Iberian forests. Global Ecology and Biogeography, 2013, 22(10): 1141- 1151.

    [12] Law R, Illian J, Burslem D F R P, Gratzer G, Gunatilleke C V S, Gunatilleke I A U N. Ecological information from spatial patterns of plants: insights from point process theory. Journal of Ecology, 2009, 97(4): 616- 628.

    [13] Picard N, Bar-Hen A, Mortier F, Chadoeuf J. Understanding the dynamics of an undisturbed tropical rain forest from the spatial pattern of trees. Journal of Ecology, 2009, 97(1): 97- 108.

    [14] Stohlgren T J, Bachand R R, Onami Y, Binkley D. Species-environment relationships and vegetation patterns: effects of spatial scale and tree life-stage. Plant Ecology, 1998, 135(2): 215- 228.

    [15] Getzin S, Wiegand T, Wiegand K, He F L. Heterogeneity influences spatial patterns and demographics in forest stands. Journal of Ecology, 2008, 96(4): 807- 820.

    [16] 郝珉輝, 張忠輝, 趙珊珊, 張春雨, 趙秀海. 吉林蛟河針闊混交林樹木生長空間關聯(lián)格局研究. 生態(tài)學報, 2017, doi: 10.5846/stxb201510272172.

    [17] 方精云, 王襄平, 沈澤昊, 唐志堯, 賀金生, 于丹, 江源, 王志恒, 鄭成洋, 朱江玲, 郭兆迪. 植物群落清查的主要內(nèi)容、方法和技術規(guī)范. 生物多樣性, 2009, 17(6): 533- 548.

    [18] White D A, Hood C S. Vegetation patterns and environmental gradients in tropical dry forests of the northern Yucatan Peninsula. Journal of Vegetation Science, 2004, 15(2): 151- 160.

    [19] Tsujino R, Takafumi H, Agetsuma N, Yumoto T. Variation in tree growth, mortality and recruitment among topographic positions in a warm temperate forest. Journal of Vegetation Science, 2006, 17(3): 281- 290.

    [20] Balvanera P, Quijas S, Pérez-Jiménez A. Distribution patterns of tropical dry forest trees along a mesoscale water availability gradient. Biotropica, 2011, 43(4): 414- 422.

    [21] Scherrer D, K?rner C. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. Journal of Biogeography, 2011, 38(2): 406- 416.

    [22] Yamakura T, Kanzaki M, Itoh A, Ohkubo T, Ogino K, Chai E O K, Lee H S, Ashton P S. Topography of a large-scale research plot established within a tropical rain forest at Lambir, Sarawak. Tropics, 1995, 5(1/2): 41- 56.

    [23] De′ath G, Fabricius K E. Classification and regression trees: a powerful yet simple technique for ecological data analysis. Ecology, 2000, 81(11): 3178- 3192.

    [24] De′ath G. Multivariate regression trees: a new technique for modeling species-environment relationships. Ecology, 2002, 83(4): 1105- 1117.

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

    [26] Aiba S I, Kitayama K, Takyu M. Habitat associations with topography and canopy structure of tree species in a tropical montane forest on Mount Kinabalu, Borneo. Plant Ecology, 2004, 174(1): 147- 161.

    [27] Coomes D A, Grubb P J. Impacts of root competition in forests and woodlands: a theoretical framework and review of experiments. Ecological Monographs, 2000, 70(2): 171- 207.

    [28] Takyu M, Aiba S I, Kitayama K. Effects of topography on tropical lower montane forests under different geological conditions on Mount Kinabalu, Borneo. Plant Ecology, 2002, 159(1): 35- 49.

    [29] Yamada T, Zuidema P A, Itoh A, Yamakura T, Ohkubo T, Kanzaki M, Tan S, Ashton P S. Strong habitat preference of a tropical rain forest tree does not imply large differences in population dynamics across habitats. Journal of Ecology, 2007, 95(2): 332- 342.

    [30] Prinzing A. The niche of higher plants: evidence for phylogenetic conservatism. Proceedings of the Royal Society of London B: Biological Sciences, 2001, 268(1483): 2383- 2389.

    [31] Park D S, Potter D. Why close relatives make bad neighbours: phylogenetic conservatism in niche preferences and dispersal disproves Darwin′s naturalization hypothesis in the thistle tribe. Molecular Ecology, 2015, 24(12): 3181- 3193.

    [32] Gunatilleke C V S, Gunatilleke I A U N, Esufali S, Harms K E, Ashton P M S, Burslem D F R P, Ashton P S. Species-habitat associations in a Sri Lankan dipterocarp forest. Journal of Tropical Ecology, 2006, 22(4): 371- 384.

    [33] Yamada T, Tomita A, Itoh A, Yamakura T, Ohkubo T, Kanzaki M, Tan S, Ashton P S. Habitat associations ofSterculiaceaetrees in a Bornean rain forest plot. Journal of Vegetation Science, 2006, 17(5): 559- 566.

    [34] Flynn D F B, Mirotchnick N, Jain M, Palmer M I, Naeem S. Functional and phylogenetic diversity as predictors of biodiversity-ecosystem-function relationships. Ecology, 2011, 92(8): 1573- 1581.

    [35] Russell M B, Woodall C W, D′Amato A W, Domke G M, Saatchi S S. Beyond mean functional traits: influence of functional trait profiles on forest structure, production, and mortality across the eastern US. Forest Ecology and Management, 2014, 328: 1- 9.

    Habitat associations of tree growth in a coniferous and broad-leaved mixed forest in Jiaohe, Jilin Province

    HAO Minhui1, ZHANG Zhonghui2, ZHAO Shanshan2, ZHAO Xiuhai1, BAIKETUERHAN Yeerjiang3, ZHANG Chunyu1,*

    1TheKeyLaboratoryforForestResources&EcosystemProcessesofBeijing,BeijingForestryUniversity,Beijing100083,China2AcademyofForestryResearchofJilinProvince,Changchun130033,China3CollegeofForestryandHorticulture,XinjiangAgriculturalUniversity,Urumqi830052,China

    The goal of this study was to determine the relative contribution of habitat differentiation on the radial growth of forest trees and to understand some underlying ecological mechanisms that may affect these relationships. Forest growth data for the period 2009 to 2014 were collected in a 21.12 hm2stem-mapped temperate coniferous and broadleaved forest in Jiaohe, Jilin Province. Based on habitat type classifications, three ecological null models, Habitat Complete Spatial Randomness (CSR) Process, Species Complete Spatial Randomness Process and Species Thomas Process were employed to test the significance of radial growth-habitat associations. Furthermore, Pearson correlation analysis was used to explore the influence of topographical variables on radial growth. A total of 36 species were included in the analysis. The following results were obtained: (1) Sixty one percent of radial growth exhibited significant associations with habitat types based on the Species Thomas Process; whereas for the Species CSR Process and the Habitat CSR Process the corresponding percentages were 94.4% and 86.1%, respectively. (2) Comparing the discrepancies among species of different life forms, we found that the habitat preferences of shrubs and small trees were similar, whereas trees had significantly different habitat preferences. The radial growth of shrubs and small trees were positively associated with habitat type 1 at a relatively low altitude and on flat terrain, but were negatively associated with habitat types 2 and 4, which had a relatively high elevation and on steep slopes. Conversely, trees were positively associated with habitat type 4, but negatively associated with habitat type 1. Furthermore, the radial growth of trees was highly sensitive to habitat heterogeneity when compared with that of shrubs and small trees. (3) Pearson correlation analysis indicated that up to 86.1 percent of growth rates were significantly correlated with one or more topographical variables. Among the four topographical variables evaluated, elevation turned out to be the most influential factor for radial growth, followed by slope and convexity, whereas aspect was the least influential. The results of this study suggest that habitat heterogeneity and niche differentiation have a considerable effect on radial growth in the temperate coniferous and broadleaved forest.

    radial growth; habitat; topography; habitat differentiation; life form

    國家自然科學基金項目(31670643, 31600480)

    2016- 03- 03; 網(wǎng)絡出版日期:2017- 02- 17

    10.5846/stxb201603030376

    *通訊作者Corresponding author.E-mail: zcy_0520@163.com

    郝珉輝, 張忠輝, 趙珊珊, 趙秀海, 葉爾江·拜克吐爾漢, 張春雨.吉林蛟河針闊混交林樹木生長與生境的關聯(lián)性.生態(tài)學報,2017,37(10):3437- 3444.

    Hao M H, Zhang Z H, Zhao S S, Zhao X H, Baiketuerhan Yeerjiang, Zhang C Y.Habitat associations of tree growth in a coniferous and broad-leaved mixed forest in Jiaohe, Jilin Province.Acta Ecologica Sinica,2017,37(10):3437- 3444.

    猜你喜歡
    關聯(lián)性生境樣地
    額爾古納市興安落葉松中齡林植被碳儲量研究
    棗樹適應干旱生境研究進展
    落葉果樹(2021年6期)2021-02-12 01:29:10
    昆明市主要綠化樹種閾值測定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構分析
    四物湯有效成分的關聯(lián)性分析
    中成藥(2017年3期)2017-05-17 06:09:05
    如何準確認定排污行為和環(huán)境損害之間的關聯(lián)性
    金秀瑤族長鼓舞傳承的文化生境
    體育科技(2016年2期)2016-02-28 17:06:07
    CRP檢測與新生兒感染的關聯(lián)性
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    廣西石漠化地區(qū)不同植被恢復樣地的土壤微生物量比較
    女人爽到高潮嗷嗷叫在线视频| 一区二区三区四区激情视频| 人人妻,人人澡人人爽秒播 | 亚洲欧洲国产日韩| 男女午夜视频在线观看| 日韩 亚洲 欧美在线| 中文字幕av电影在线播放| 成年人黄色毛片网站| 男女之事视频高清在线观看 | 男人爽女人下面视频在线观看| 国产麻豆69| 国产淫语在线视频| 亚洲国产成人一精品久久久| 欧美成人精品欧美一级黄| 一二三四在线观看免费中文在| 欧美日韩黄片免| 99国产精品99久久久久| 国产精品免费大片| 天天操日日干夜夜撸| 男女床上黄色一级片免费看| 少妇精品久久久久久久| 99热全是精品| 国产精品一区二区精品视频观看| av国产久精品久网站免费入址| 欧美日韩视频高清一区二区三区二| 日韩av免费高清视频| 精品熟女少妇八av免费久了| 大陆偷拍与自拍| 最近中文字幕2019免费版| 国产熟女欧美一区二区| 成人免费观看视频高清| 人妻一区二区av| 精品一品国产午夜福利视频| 亚洲国产毛片av蜜桃av| 精品一区二区三卡| 午夜福利视频在线观看免费| 午夜影院在线不卡| 美女中出高潮动态图| 国产精品久久久人人做人人爽| 巨乳人妻的诱惑在线观看| 日韩电影二区| 亚洲精品美女久久av网站| 亚洲人成77777在线视频| 免费高清在线观看日韩| 狠狠婷婷综合久久久久久88av| 欧美激情高清一区二区三区| 国产精品久久久久久人妻精品电影 | 热re99久久国产66热| 黄色 视频免费看| 亚洲av美国av| 精品国产乱码久久久久久小说| 久久精品久久精品一区二区三区| 日韩制服丝袜自拍偷拍| 色94色欧美一区二区| 国产国语露脸激情在线看| 日日爽夜夜爽网站| 飞空精品影院首页| 国产男女超爽视频在线观看| 国产精品亚洲av一区麻豆| 一级毛片我不卡| 新久久久久国产一级毛片| 中国国产av一级| 成人影院久久| 麻豆乱淫一区二区| 欧美激情极品国产一区二区三区| 99久久精品国产亚洲精品| 精品卡一卡二卡四卡免费| 韩国精品一区二区三区| 久久人妻熟女aⅴ| 一本—道久久a久久精品蜜桃钙片| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 人成视频在线观看免费观看| 欧美精品av麻豆av| 亚洲国产欧美日韩在线播放| 久久国产精品大桥未久av| 国产精品秋霞免费鲁丝片| 男女边吃奶边做爰视频| √禁漫天堂资源中文www| 午夜久久久在线观看| 夫妻午夜视频| 久久综合国产亚洲精品| 99九九在线精品视频| 亚洲中文字幕日韩| 中文字幕亚洲精品专区| 国产精品一二三区在线看| 侵犯人妻中文字幕一二三四区| 色婷婷av一区二区三区视频| 午夜老司机福利片| 老鸭窝网址在线观看| 美女脱内裤让男人舔精品视频| kizo精华| 女人被躁到高潮嗷嗷叫费观| 国产日韩欧美视频二区| 下体分泌物呈黄色| 好男人视频免费观看在线| 黄片播放在线免费| 国产高清国产精品国产三级| 中文字幕av电影在线播放| 欧美少妇被猛烈插入视频| 欧美黄色淫秽网站| 操出白浆在线播放| 大香蕉久久成人网| 交换朋友夫妻互换小说| 国产免费现黄频在线看| 黄片小视频在线播放| 欧美av亚洲av综合av国产av| 国产高清视频在线播放一区 | av一本久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 日韩大片免费观看网站| 精品一区二区三区四区五区乱码 | 极品人妻少妇av视频| 日韩大码丰满熟妇| 大陆偷拍与自拍| 久久久国产精品麻豆| 19禁男女啪啪无遮挡网站| 亚洲国产精品999| 亚洲精品一二三| 亚洲国产精品国产精品| 亚洲精品日本国产第一区| 中国国产av一级| 熟女少妇亚洲综合色aaa.| 日本一区二区免费在线视频| 日韩av不卡免费在线播放| 国产女主播在线喷水免费视频网站| 久久鲁丝午夜福利片| 1024香蕉在线观看| 宅男免费午夜| 宅男免费午夜| 亚洲av美国av| 国产精品九九99| 欧美日韩亚洲综合一区二区三区_| av有码第一页| 亚洲av美国av| 蜜桃在线观看..| 国产午夜精品一二区理论片| 久久精品aⅴ一区二区三区四区| 99久久人妻综合| 国产欧美日韩一区二区三区在线| 悠悠久久av| 精品国产超薄肉色丝袜足j| 97人妻天天添夜夜摸| 午夜两性在线视频| 一区二区三区激情视频| 久久国产精品男人的天堂亚洲| 欧美老熟妇乱子伦牲交| 久久久精品免费免费高清| 久久狼人影院| 亚洲国产毛片av蜜桃av| 亚洲少妇的诱惑av| 国产在线观看jvid| 日韩精品免费视频一区二区三区| 青春草视频在线免费观看| 欧美黑人精品巨大| 18禁国产床啪视频网站| 亚洲欧美色中文字幕在线| 黄片小视频在线播放| 国产欧美日韩精品亚洲av| 国产成人一区二区在线| 国产在线视频一区二区| 免费在线观看完整版高清| 涩涩av久久男人的天堂| 国产亚洲av高清不卡| 亚洲精品国产区一区二| 狠狠精品人妻久久久久久综合| 日韩电影二区| h视频一区二区三区| 精品一品国产午夜福利视频| 18禁国产床啪视频网站| 午夜影院在线不卡| 熟女av电影| 国产免费现黄频在线看| 亚洲激情五月婷婷啪啪| 日本wwww免费看| av网站免费在线观看视频| 国产精品 国内视频| 母亲3免费完整高清在线观看| 1024视频免费在线观看| 国产成人影院久久av| 国产福利在线免费观看视频| 69精品国产乱码久久久| 免费不卡黄色视频| 欧美另类一区| 日本av手机在线免费观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲av片天天在线观看| 亚洲国产最新在线播放| netflix在线观看网站| av一本久久久久| 免费一级毛片在线播放高清视频 | 免费不卡黄色视频| www.精华液| 啦啦啦视频在线资源免费观看| 免费在线观看视频国产中文字幕亚洲 | 午夜av观看不卡| 嫁个100分男人电影在线观看 | 日韩免费高清中文字幕av| 欧美人与性动交α欧美软件| 日日爽夜夜爽网站| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲高清精品| 国产精品久久久人人做人人爽| 亚洲精品久久成人aⅴ小说| 免费看不卡的av| 大香蕉久久网| 亚洲激情五月婷婷啪啪| 满18在线观看网站| 久久久久久久久久久久大奶| 满18在线观看网站| 亚洲一区二区三区欧美精品| 伊人亚洲综合成人网| 精品国产一区二区久久| 国产人伦9x9x在线观看| 一级毛片我不卡| 18禁观看日本| 性高湖久久久久久久久免费观看| 免费高清在线观看视频在线观看| 首页视频小说图片口味搜索 | 久热这里只有精品99| 国产免费福利视频在线观看| 在线观看免费高清a一片| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 一级黄色大片毛片| 肉色欧美久久久久久久蜜桃| 久久久久久久久免费视频了| 一区福利在线观看| 成人国产av品久久久| svipshipincom国产片| 一二三四在线观看免费中文在| 人妻 亚洲 视频| 交换朋友夫妻互换小说| 自拍欧美九色日韩亚洲蝌蚪91| 成人影院久久| 亚洲国产欧美网| 免费在线观看影片大全网站 | 欧美日韩综合久久久久久| 免费高清在线观看视频在线观看| 日韩av免费高清视频| 国产精品一区二区精品视频观看| 久久久精品免费免费高清| 狠狠婷婷综合久久久久久88av| 亚洲av成人精品一二三区| 一级毛片女人18水好多 | 成年动漫av网址| 国产精品一区二区在线观看99| 亚洲熟女精品中文字幕| 午夜激情久久久久久久| 国产日韩一区二区三区精品不卡| 在线观看国产h片| 免费一级毛片在线播放高清视频 | 欧美少妇被猛烈插入视频| 美国免费a级毛片| 亚洲欧美中文字幕日韩二区| 欧美日韩黄片免| 久久狼人影院| 久久人妻福利社区极品人妻图片 | 99九九在线精品视频| 国产精品一二三区在线看| 亚洲av在线观看美女高潮| 国产欧美日韩综合在线一区二区| 婷婷丁香在线五月| 久久青草综合色| 亚洲国产精品成人久久小说| 韩国高清视频一区二区三区| www.自偷自拍.com| 亚洲视频免费观看视频| 男女之事视频高清在线观看 | 日韩av不卡免费在线播放| 操出白浆在线播放| av网站在线播放免费| 看免费av毛片| 天堂中文最新版在线下载| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 久久久亚洲精品成人影院| 久久久久视频综合| 欧美国产精品va在线观看不卡| 无遮挡黄片免费观看| 亚洲精品美女久久av网站| 国产精品 国内视频| 亚洲综合色网址| av在线播放精品| 亚洲欧美日韩另类电影网站| 日本91视频免费播放| 国产亚洲精品久久久久5区| 女警被强在线播放| 久久久久久久久免费视频了| 五月天丁香电影| 嫩草影视91久久| 18禁国产床啪视频网站| 国产女主播在线喷水免费视频网站| 老司机午夜十八禁免费视频| 热re99久久国产66热| 亚洲精品中文字幕在线视频| h视频一区二区三区| 女人久久www免费人成看片| 中文字幕人妻丝袜制服| 视频区图区小说| 久久热在线av| 下体分泌物呈黄色| 悠悠久久av| 精品国产一区二区三区四区第35| 午夜免费男女啪啪视频观看| 欧美日韩亚洲综合一区二区三区_| 国产伦理片在线播放av一区| 国产男女超爽视频在线观看| 国产色视频综合| 精品久久蜜臀av无| 国产男人的电影天堂91| 亚洲色图综合在线观看| 黑人欧美特级aaaaaa片| 欧美精品亚洲一区二区| 男女无遮挡免费网站观看| 最近中文字幕2019免费版| 成人手机av| 中文字幕最新亚洲高清| 久久久久网色| 国产日韩欧美在线精品| 99国产综合亚洲精品| 18禁黄网站禁片午夜丰满| 久久国产精品人妻蜜桃| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产欧美在线一区| 99久久99久久久精品蜜桃| 亚洲精品中文字幕在线视频| 国产精品二区激情视频| 久久九九热精品免费| 啦啦啦中文免费视频观看日本| 日韩一区二区三区影片| a级片在线免费高清观看视频| 一级,二级,三级黄色视频| 男女高潮啪啪啪动态图| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产欧美网| 女人被躁到高潮嗷嗷叫费观| 亚洲男人天堂网一区| 这个男人来自地球电影免费观看| 又粗又硬又长又爽又黄的视频| av电影中文网址| 美女午夜性视频免费| 夫妻性生交免费视频一级片| 美女主播在线视频| 人人妻人人添人人爽欧美一区卜| 97人妻天天添夜夜摸| 久久亚洲国产成人精品v| 久久人妻熟女aⅴ| 久久久精品国产亚洲av高清涩受| 日本猛色少妇xxxxx猛交久久| 国产成人a∨麻豆精品| 国产成人欧美| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 国产成人91sexporn| 国产成人免费无遮挡视频| 丝袜在线中文字幕| 91麻豆精品激情在线观看国产 | 亚洲精品日本国产第一区| 国产淫语在线视频| 国产黄色免费在线视频| 国产精品一区二区免费欧美 | 天天躁日日躁夜夜躁夜夜| 老熟女久久久| 精品国产一区二区三区久久久樱花| 国产精品秋霞免费鲁丝片| 一级毛片我不卡| 少妇精品久久久久久久| 叶爱在线成人免费视频播放| 亚洲精品一卡2卡三卡4卡5卡 | 欧美久久黑人一区二区| 国产日韩欧美视频二区| 欧美激情极品国产一区二区三区| 最新在线观看一区二区三区 | 成年av动漫网址| 国产欧美亚洲国产| 满18在线观看网站| 一级毛片电影观看| 久热爱精品视频在线9| 赤兔流量卡办理| 婷婷色麻豆天堂久久| 久久99热这里只频精品6学生| 丝瓜视频免费看黄片| 啦啦啦 在线观看视频| 欧美日韩成人在线一区二区| 高清av免费在线| 亚洲国产精品成人久久小说| 精品视频人人做人人爽| 大话2 男鬼变身卡| 色综合欧美亚洲国产小说| 欧美精品亚洲一区二区| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 一本大道久久a久久精品| 精品一品国产午夜福利视频| 亚洲国产精品一区二区三区在线| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 亚洲精品一二三| 18禁黄网站禁片午夜丰满| 日韩免费高清中文字幕av| 黑人巨大精品欧美一区二区蜜桃| 人人妻人人澡人人爽人人夜夜| 日韩制服骚丝袜av| 国产精品久久久久成人av| 日韩精品免费视频一区二区三区| 制服诱惑二区| 一级毛片我不卡| 国产精品久久久人人做人人爽| 青草久久国产| 脱女人内裤的视频| 亚洲av在线观看美女高潮| 欧美激情高清一区二区三区| 国产精品人妻久久久影院| 下体分泌物呈黄色| 亚洲,欧美精品.| 老鸭窝网址在线观看| 久久久久视频综合| 国产精品99久久99久久久不卡| 国产一卡二卡三卡精品| 麻豆国产av国片精品| 丁香六月欧美| 男人舔女人的私密视频| 纵有疾风起免费观看全集完整版| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 啦啦啦在线免费观看视频4| 日本91视频免费播放| 18在线观看网站| 精品少妇内射三级| 成年动漫av网址| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| 精品久久久久久久毛片微露脸 | 亚洲第一av免费看| 国产成人影院久久av| 亚洲精品久久午夜乱码| 久久国产精品男人的天堂亚洲| 日韩av在线免费看完整版不卡| 免费在线观看日本一区| 国产亚洲精品第一综合不卡| 欧美中文综合在线视频| 亚洲精品av麻豆狂野| 在线观看免费视频网站a站| 后天国语完整版免费观看| av在线播放精品| 国产精品久久久av美女十八| 又大又爽又粗| 久久精品久久精品一区二区三区| 国产精品久久久av美女十八| 久久久亚洲精品成人影院| 人人妻人人添人人爽欧美一区卜| 欧美精品高潮呻吟av久久| 国产精品偷伦视频观看了| 99九九在线精品视频| 国产又色又爽无遮挡免| 久9热在线精品视频| 中文字幕精品免费在线观看视频| av视频免费观看在线观看| 99热国产这里只有精品6| 欧美日韩精品网址| 美女高潮到喷水免费观看| 婷婷丁香在线五月| 在线观看www视频免费| 人人妻人人澡人人爽人人夜夜| 自拍欧美九色日韩亚洲蝌蚪91| 秋霞在线观看毛片| 丝袜美腿诱惑在线| 一级毛片电影观看| 中文字幕亚洲精品专区| 久久国产精品大桥未久av| 午夜激情久久久久久久| 亚洲精品国产区一区二| 大话2 男鬼变身卡| 国产成人av激情在线播放| 日韩 欧美 亚洲 中文字幕| 女警被强在线播放| 一区二区三区激情视频| 国产亚洲精品久久久久5区| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 啦啦啦中文免费视频观看日本| 国产一区有黄有色的免费视频| 69精品国产乱码久久久| 久久精品亚洲熟妇少妇任你| 国产91精品成人一区二区三区 | 好男人视频免费观看在线| 成人午夜精彩视频在线观看| 亚洲伊人久久精品综合| 下体分泌物呈黄色| 精品久久久久久久毛片微露脸 | 日韩一区二区三区影片| 国产精品九九99| 黑人猛操日本美女一级片| 丁香六月天网| 一区二区av电影网| 777久久人妻少妇嫩草av网站| 免费女性裸体啪啪无遮挡网站| 男男h啪啪无遮挡| 蜜桃在线观看..| 熟女少妇亚洲综合色aaa.| 亚洲精品国产av蜜桃| 亚洲av国产av综合av卡| 久久精品国产综合久久久| 熟女av电影| 亚洲精品乱久久久久久| 亚洲图色成人| 国产精品国产三级专区第一集| 99国产精品一区二区三区| 国产在线视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 精品人妻一区二区三区麻豆| 后天国语完整版免费观看| 国产精品一二三区在线看| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 国产伦理片在线播放av一区| 大码成人一级视频| 丝瓜视频免费看黄片| 两性夫妻黄色片| 亚洲欧美日韩另类电影网站| 久久国产精品大桥未久av| 午夜免费鲁丝| av视频免费观看在线观看| 一级毛片电影观看| av天堂在线播放| 少妇人妻 视频| 国产亚洲一区二区精品| 性色av乱码一区二区三区2| 亚洲色图综合在线观看| 免费在线观看完整版高清| 欧美中文综合在线视频| 国产一区有黄有色的免费视频| 国产国语露脸激情在线看| 国产精品国产三级国产专区5o| 国产精品久久久久久人妻精品电影 | 精品卡一卡二卡四卡免费| 久久中文字幕一级| 女性被躁到高潮视频| 一级,二级,三级黄色视频| 欧美 亚洲 国产 日韩一| 亚洲av电影在线进入| av国产精品久久久久影院| a级片在线免费高清观看视频| 亚洲欧洲日产国产| 久久狼人影院| 久久久久视频综合| 国产又爽黄色视频| 少妇的丰满在线观看| 免费黄频网站在线观看国产| 一区二区av电影网| 伊人久久大香线蕉亚洲五| 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 欧美日韩成人在线一区二区| 亚洲欧美激情在线| 亚洲五月婷婷丁香| 国产精品一二三区在线看| 中国美女看黄片| 亚洲精品成人av观看孕妇| 久久久久久久久久久久大奶| a 毛片基地| 老鸭窝网址在线观看| 亚洲国产精品成人久久小说| 最近最新中文字幕大全免费视频 | 欧美激情 高清一区二区三区| 欧美日韩一级在线毛片| 激情视频va一区二区三区| 亚洲精品自拍成人| 99热全是精品| 如日韩欧美国产精品一区二区三区| 一级黄色大片毛片| 1024视频免费在线观看| www日本在线高清视频| 高清黄色对白视频在线免费看| 中文乱码字字幕精品一区二区三区| 天天操日日干夜夜撸| 久久久久网色| 亚洲一区二区三区欧美精品| 亚洲三区欧美一区| 嫁个100分男人电影在线观看 | 纵有疾风起免费观看全集完整版| 丰满少妇做爰视频| 香蕉丝袜av| 国产精品熟女久久久久浪| cao死你这个sao货| 狂野欧美激情性bbbbbb| 亚洲一码二码三码区别大吗| 国产老妇伦熟女老妇高清| 交换朋友夫妻互换小说| 黄色视频在线播放观看不卡| 一级,二级,三级黄色视频| 侵犯人妻中文字幕一二三四区| 欧美在线一区亚洲| 欧美日韩av久久| 天天躁夜夜躁狠狠久久av| 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91| 一级黄色大片毛片| 亚洲熟女精品中文字幕| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| 天天躁狠狠躁夜夜躁狠狠躁| 美女扒开内裤让男人捅视频| 久久精品aⅴ一区二区三区四区| svipshipincom国产片| 欧美激情 高清一区二区三区| 高潮久久久久久久久久久不卡| 岛国毛片在线播放| 婷婷色综合www| 大型av网站在线播放|