• <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ū)不同植被恢復樣地的土壤微生物量比較
    精品欧美一区二区三区在线| 侵犯人妻中文字幕一二三四区| 日本免费一区二区三区高清不卡 | 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 嫩草影视91久久| 亚洲激情在线av| 欧美成人性av电影在线观看| 亚洲欧美精品综合久久99| 久久精品国产综合久久久| 999精品在线视频| 久久中文看片网| 99国产精品一区二区蜜桃av| 免费久久久久久久精品成人欧美视频| 又黄又粗又硬又大视频| 亚洲成人精品中文字幕电影| 亚洲av电影在线进入| 欧美日韩一级在线毛片| 免费观看精品视频网站| 又紧又爽又黄一区二区| 国产精品免费一区二区三区在线| 高清黄色对白视频在线免费看| 免费在线观看亚洲国产| 香蕉久久夜色| 欧美av亚洲av综合av国产av| 久久精品国产亚洲av高清一级| 桃色一区二区三区在线观看| 免费av毛片视频| 色播亚洲综合网| 久久草成人影院| 国产成人精品久久二区二区免费| 夜夜躁狠狠躁天天躁| 国产aⅴ精品一区二区三区波| 久久久国产成人精品二区| 国产极品粉嫩免费观看在线| 黑人欧美特级aaaaaa片| 亚洲 欧美 日韩 在线 免费| 亚洲成av片中文字幕在线观看| 黄片大片在线免费观看| 啦啦啦观看免费观看视频高清 | 亚洲成av人片免费观看| 精品国产乱子伦一区二区三区| 精品不卡国产一区二区三区| 亚洲男人天堂网一区| 桃红色精品国产亚洲av| 夜夜躁狠狠躁天天躁| 999久久久精品免费观看国产| 国产一卡二卡三卡精品| 成年人黄色毛片网站| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| 琪琪午夜伦伦电影理论片6080| 香蕉国产在线看| 国产成人av激情在线播放| 日本黄色视频三级网站网址| www.999成人在线观看| 国产亚洲av嫩草精品影院| 国产aⅴ精品一区二区三区波| 亚洲三区欧美一区| 欧美日本中文国产一区发布| 麻豆国产av国片精品| 国内精品久久久久精免费| 久久久久久免费高清国产稀缺| 日韩欧美三级三区| 亚洲国产精品成人综合色| 一区二区三区精品91| 午夜精品久久久久久毛片777| 波多野结衣一区麻豆| 欧美另类亚洲清纯唯美| 国产精品乱码一区二三区的特点 | 国产亚洲欧美在线一区二区| 欧美日韩亚洲综合一区二区三区_| 欧美日韩一级在线毛片| 久久人人爽av亚洲精品天堂| 免费无遮挡裸体视频| 在线观看www视频免费| av中文乱码字幕在线| 久久人人97超碰香蕉20202| 美女免费视频网站| 天堂√8在线中文| 大型av网站在线播放| 精品国产一区二区久久| 国产人伦9x9x在线观看| 777久久人妻少妇嫩草av网站| 久久婷婷成人综合色麻豆| 亚洲国产看品久久| 女性生殖器流出的白浆| 青草久久国产| www.精华液| 欧美乱色亚洲激情| ponron亚洲| 亚洲久久久国产精品| 电影成人av| 首页视频小说图片口味搜索| 国产99久久九九免费精品| 欧美绝顶高潮抽搐喷水| 日韩免费av在线播放| 国产高清激情床上av| 免费女性裸体啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 色av中文字幕| 黄色a级毛片大全视频| 国产熟女xx| 国产精品亚洲一级av第二区| 黄色 视频免费看| 日本五十路高清| 欧美黑人精品巨大| 久久久久久国产a免费观看| 午夜亚洲福利在线播放| 久久九九热精品免费| 最新在线观看一区二区三区| 黑人欧美特级aaaaaa片| 老熟妇乱子伦视频在线观看| 亚洲精品国产区一区二| 黄色女人牲交| 久久人妻福利社区极品人妻图片| 国产色视频综合| 97超级碰碰碰精品色视频在线观看| 中文字幕最新亚洲高清| 国产精品亚洲美女久久久| 丝袜在线中文字幕| 91av网站免费观看| 欧洲精品卡2卡3卡4卡5卡区| 99久久久亚洲精品蜜臀av| 露出奶头的视频| 视频在线观看一区二区三区| 少妇的丰满在线观看| 欧美不卡视频在线免费观看 | 99国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 一级,二级,三级黄色视频| 黄色a级毛片大全视频| 日本免费一区二区三区高清不卡 | 精品福利观看| 欧美日本中文国产一区发布| 老司机福利观看| 午夜激情av网站| 国产蜜桃级精品一区二区三区| 国产亚洲av嫩草精品影院| 成人精品一区二区免费| 久久婷婷人人爽人人干人人爱 | 成人国产一区最新在线观看| 国产成人精品久久二区二区91| 午夜久久久久精精品| 午夜视频精品福利| 母亲3免费完整高清在线观看| 最近最新中文字幕大全电影3 | 人人妻,人人澡人人爽秒播| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕最新亚洲高清| 人人妻人人澡欧美一区二区 | 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频| 久久久久久亚洲精品国产蜜桃av| 亚洲中文字幕一区二区三区有码在线看 | 国产av又大| 亚洲中文日韩欧美视频| av在线播放免费不卡| 两性夫妻黄色片| 看黄色毛片网站| 日本欧美视频一区| 国产欧美日韩一区二区精品| 色老头精品视频在线观看| 丰满的人妻完整版| 国产成人av教育| 欧美国产精品va在线观看不卡| 色综合亚洲欧美另类图片| 国产熟女午夜一区二区三区| 久久久久国产精品人妻aⅴ院| 亚洲第一青青草原| 欧美中文日本在线观看视频| 女同久久另类99精品国产91| 国内毛片毛片毛片毛片毛片| 久久久久久久久免费视频了| 日韩 欧美 亚洲 中文字幕| 亚洲成av人片免费观看| 欧美精品啪啪一区二区三区| 亚洲精品美女久久av网站| 桃色一区二区三区在线观看| 亚洲中文字幕一区二区三区有码在线看 | xxx96com| 亚洲国产毛片av蜜桃av| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 日韩欧美一区视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| 久久久久久大精品| 久久伊人香网站| 国产高清视频在线播放一区| 久久久久久久久久久久大奶| 99国产精品99久久久久| 亚洲视频免费观看视频| 国产亚洲欧美98| 少妇 在线观看| www.999成人在线观看| 男人的好看免费观看在线视频 | 在线免费观看的www视频| 久久精品国产亚洲av高清一级| 午夜福利成人在线免费观看| www.自偷自拍.com| 一本久久中文字幕| 亚洲最大成人中文| 久久精品人人爽人人爽视色| 亚洲天堂国产精品一区在线| 成人三级黄色视频| 久久性视频一级片| 成年版毛片免费区| 亚洲无线在线观看| 久久伊人香网站| 久久欧美精品欧美久久欧美| 男人的好看免费观看在线视频 | 欧美一级a爱片免费观看看 | 久久久久亚洲av毛片大全| av视频在线观看入口| 免费观看精品视频网站| 国产精品久久久久久精品电影 | 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲欧美精品综合一区二区三区| 国产精品精品国产色婷婷| 9色porny在线观看| 国产精品久久久av美女十八| 三级毛片av免费| 又紧又爽又黄一区二区| 欧美日韩亚洲综合一区二区三区_| 9色porny在线观看| 日本 av在线| 精品国产一区二区久久| 久久婷婷成人综合色麻豆| 亚洲av熟女| 久久 成人 亚洲| 婷婷六月久久综合丁香| 欧美色欧美亚洲另类二区 | 99久久国产精品久久久| 精品国产亚洲在线| 国产精品 国内视频| 亚洲av片天天在线观看| 99国产精品免费福利视频| 午夜福利成人在线免费观看| 免费看a级黄色片| 给我免费播放毛片高清在线观看| 欧美不卡视频在线免费观看 | 18禁国产床啪视频网站| 成人三级做爰电影| 女人高潮潮喷娇喘18禁视频| 久99久视频精品免费| 777久久人妻少妇嫩草av网站| 亚洲av电影在线进入| 久久热在线av| 午夜福利一区二区在线看| 欧美色欧美亚洲另类二区 | 日本撒尿小便嘘嘘汇集6| 亚洲人成网站在线播放欧美日韩| 亚洲av熟女| 丝袜美足系列| 黄色视频不卡| 精品久久久久久久久久免费视频| 亚洲国产看品久久| 狂野欧美激情性xxxx| 一进一出抽搐gif免费好疼| 一个人观看的视频www高清免费观看 | 男女床上黄色一级片免费看| 老司机福利观看| 伊人久久大香线蕉亚洲五| 亚洲国产精品成人综合色| 亚洲精品av麻豆狂野| 国产精品久久久久久精品电影 | 免费在线观看视频国产中文字幕亚洲| 成人免费观看视频高清| 欧美最黄视频在线播放免费| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 国产不卡一卡二| 精品无人区乱码1区二区| 国产一区二区三区视频了| 久久中文字幕一级| 99久久99久久久精品蜜桃| 岛国视频午夜一区免费看| √禁漫天堂资源中文www| 精品人妻1区二区| 好男人在线观看高清免费视频 | 1024香蕉在线观看| 亚洲第一电影网av| 日本 av在线| 国产单亲对白刺激| 国产三级在线视频| 国产亚洲欧美98| 亚洲专区中文字幕在线| 制服诱惑二区| 亚洲精品在线美女| 这个男人来自地球电影免费观看| 国产乱人伦免费视频| 日本免费a在线| 国产成人免费无遮挡视频| 亚洲av美国av| 最新美女视频免费是黄的| 搡老岳熟女国产| 99久久精品国产亚洲精品| 久久久久久人人人人人| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美激情在线| 亚洲人成电影观看| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站| 成年版毛片免费区| 久久久久久久久久久久大奶| 免费无遮挡裸体视频| 久久人妻熟女aⅴ| 国产亚洲精品av在线| 欧美日本中文国产一区发布| 脱女人内裤的视频| 在线天堂中文资源库| 欧美久久黑人一区二区| 少妇被粗大的猛进出69影院| √禁漫天堂资源中文www| 日本欧美视频一区| 精品日产1卡2卡| 自线自在国产av| 麻豆av在线久日| 一边摸一边抽搐一进一出视频| 国产精品精品国产色婷婷| 热re99久久国产66热| 亚洲精品美女久久久久99蜜臀| 多毛熟女@视频| 中文字幕人成人乱码亚洲影| 女人被狂操c到高潮| 亚洲av成人不卡在线观看播放网| 人人妻人人澡人人看| 悠悠久久av| 亚洲五月婷婷丁香| 在线观看免费视频日本深夜| 日韩高清综合在线| 午夜两性在线视频| 成人手机av| 亚洲三区欧美一区| 午夜精品在线福利| 一区二区三区激情视频| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| 日日干狠狠操夜夜爽| 无限看片的www在线观看| 黄色视频,在线免费观看| 看黄色毛片网站| 久久久久久久午夜电影| 看黄色毛片网站| 90打野战视频偷拍视频| 99久久99久久久精品蜜桃| 午夜福利免费观看在线| 9191精品国产免费久久| 中文字幕另类日韩欧美亚洲嫩草| 不卡av一区二区三区| 香蕉丝袜av| 国产熟女xx| 国产精品精品国产色婷婷| 久久精品国产综合久久久| 淫秽高清视频在线观看| 国产熟女xx| 免费看美女性在线毛片视频| 久久精品亚洲熟妇少妇任你| 亚洲av第一区精品v没综合| 99国产精品一区二区三区| 搡老妇女老女人老熟妇| 色婷婷久久久亚洲欧美| 亚洲国产精品成人综合色| 亚洲午夜精品一区,二区,三区| 一夜夜www| 国产免费男女视频| 国产成人精品久久二区二区91| 岛国在线观看网站| av超薄肉色丝袜交足视频| 国产精品99久久99久久久不卡| 母亲3免费完整高清在线观看| 午夜激情av网站| 国产精品 国内视频| 国产精品免费视频内射| 午夜福利视频1000在线观看 | 欧美色欧美亚洲另类二区 | 国产精品一区二区在线不卡| 亚洲男人天堂网一区| 国产成人av激情在线播放| 91精品三级在线观看| 人妻丰满熟妇av一区二区三区| 国产一区二区三区综合在线观看| 亚洲色图综合在线观看| 露出奶头的视频| 搡老熟女国产l中国老女人| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一出视频| 人成视频在线观看免费观看| 国产视频一区二区在线看| 亚洲五月婷婷丁香| 日本免费一区二区三区高清不卡 | 男女下面进入的视频免费午夜 | 免费一级毛片在线播放高清视频 | 国语自产精品视频在线第100页| 母亲3免费完整高清在线观看| 少妇 在线观看| 18禁裸乳无遮挡免费网站照片 | 久久人人97超碰香蕉20202| 国产精品秋霞免费鲁丝片| 神马国产精品三级电影在线观看 | 欧美日韩福利视频一区二区| 自线自在国产av| 桃色一区二区三区在线观看| av在线播放免费不卡| 国内久久婷婷六月综合欲色啪| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| av有码第一页| 女人爽到高潮嗷嗷叫在线视频| 91av网站免费观看| 国产精品1区2区在线观看.| 久久久久九九精品影院| 成人av一区二区三区在线看| aaaaa片日本免费| 亚洲国产毛片av蜜桃av| 国产精品永久免费网站| 人人妻人人爽人人添夜夜欢视频| 国产成人免费无遮挡视频| 国产精品自产拍在线观看55亚洲| 欧美日韩瑟瑟在线播放| 这个男人来自地球电影免费观看| 女人被躁到高潮嗷嗷叫费观| 久久这里只有精品19| 国产av在哪里看| 精品欧美一区二区三区在线| 午夜免费观看网址| 看免费av毛片| 丝袜美腿诱惑在线| 国产精品一区二区精品视频观看| 欧美成狂野欧美在线观看| 一进一出抽搐gif免费好疼| 日韩一卡2卡3卡4卡2021年| 国产成人系列免费观看| 级片在线观看| 久久香蕉国产精品| 老司机深夜福利视频在线观看| 母亲3免费完整高清在线观看| 黄网站色视频无遮挡免费观看| 麻豆国产av国片精品| 波多野结衣高清无吗| 亚洲一区高清亚洲精品| 桃红色精品国产亚洲av| 欧美日韩精品网址| 后天国语完整版免费观看| 韩国精品一区二区三区| 精品久久蜜臀av无| 国产精品,欧美在线| 人妻丰满熟妇av一区二区三区| 淫秽高清视频在线观看| 又黄又爽又免费观看的视频| 国产成人影院久久av| 黄色 视频免费看| 日韩一卡2卡3卡4卡2021年| 欧美在线黄色| 黄片播放在线免费| 嫩草影院精品99| 757午夜福利合集在线观看| av超薄肉色丝袜交足视频| 国产一级毛片七仙女欲春2 | 欧美乱码精品一区二区三区| 69av精品久久久久久| 亚洲成人免费电影在线观看| 人人妻人人澡人人看| 成人三级黄色视频| 九色国产91popny在线| 亚洲中文av在线| 欧美绝顶高潮抽搐喷水| 男人的好看免费观看在线视频 | 午夜免费鲁丝| 亚洲熟妇中文字幕五十中出| 9色porny在线观看| 99热只有精品国产| 精品第一国产精品| 人人妻人人爽人人添夜夜欢视频| 人人澡人人妻人| 男女床上黄色一级片免费看| 午夜免费鲁丝| 1024香蕉在线观看| 欧美大码av| 国产免费av片在线观看野外av| 级片在线观看| 色综合站精品国产| 岛国在线观看网站| 一区在线观看完整版| 色在线成人网| 最近最新中文字幕大全电影3 | 国产成人精品无人区| 亚洲 欧美 日韩 在线 免费| av片东京热男人的天堂| 少妇被粗大的猛进出69影院| 成人精品一区二区免费| 夜夜爽天天搞| 精品电影一区二区在线| 色综合欧美亚洲国产小说| 88av欧美| 国产精品亚洲美女久久久| 亚洲精品久久国产高清桃花| 乱人伦中国视频| 国产精品 欧美亚洲| 国产一区二区在线av高清观看| 亚洲专区字幕在线| 桃红色精品国产亚洲av| 欧美日本亚洲视频在线播放| 亚洲欧洲精品一区二区精品久久久| 女生性感内裤真人,穿戴方法视频| 少妇熟女aⅴ在线视频| 两性夫妻黄色片| 91成人精品电影| 久久午夜综合久久蜜桃| 嫩草影院精品99| 侵犯人妻中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频| 日本欧美视频一区| 脱女人内裤的视频| 久久久久久久精品吃奶| ponron亚洲| 亚洲av电影不卡..在线观看| 中文字幕人妻丝袜一区二区| 欧美绝顶高潮抽搐喷水| 亚洲色图 男人天堂 中文字幕| 亚洲欧美日韩无卡精品| 九色国产91popny在线| 国产精品av久久久久免费| 久久国产精品影院| 美女 人体艺术 gogo| 国产av一区二区精品久久| 看片在线看免费视频| 日韩av在线大香蕉| 亚洲 国产 在线| 中出人妻视频一区二区| 一a级毛片在线观看| 日本在线视频免费播放| 国产私拍福利视频在线观看| 亚洲avbb在线观看| 婷婷六月久久综合丁香| 国产精品 欧美亚洲| 国产一卡二卡三卡精品| 精品国产乱子伦一区二区三区| 亚洲熟妇中文字幕五十中出| 国产99久久九九免费精品| av有码第一页| 国产亚洲精品综合一区在线观看 | 午夜两性在线视频| 在线观看免费日韩欧美大片| 亚洲中文字幕一区二区三区有码在线看 | 人成视频在线观看免费观看| 国产精品1区2区在线观看.| 这个男人来自地球电影免费观看| 国产亚洲精品久久久久久毛片| 欧美乱妇无乱码| 12—13女人毛片做爰片一| 国产麻豆成人av免费视频| 亚洲精品一区av在线观看| 91麻豆av在线| 少妇被粗大的猛进出69影院| 国产在线精品亚洲第一网站| 长腿黑丝高跟| 热re99久久国产66热| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 亚洲五月天丁香| 伊人久久大香线蕉亚洲五| 国产片内射在线| 久久中文字幕人妻熟女| 亚洲一区二区三区色噜噜| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉国产精品| 侵犯人妻中文字幕一二三四区| 桃色一区二区三区在线观看| 天堂动漫精品| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 国产高清激情床上av| 91在线观看av| 午夜福利视频1000在线观看 | 亚洲情色 制服丝袜| 狠狠狠狠99中文字幕| 男人舔女人的私密视频| 正在播放国产对白刺激| 免费女性裸体啪啪无遮挡网站| 成人av一区二区三区在线看| 欧美国产日韩亚洲一区| 黑丝袜美女国产一区| 亚洲电影在线观看av| 欧美成人免费av一区二区三区| 久久久水蜜桃国产精品网| 国产乱人伦免费视频| av片东京热男人的天堂| 亚洲av成人av| 亚洲欧美精品综合一区二区三区| 国产亚洲欧美98| 亚洲精品美女久久av网站| 国产精品免费视频内射| 女人被狂操c到高潮| 久久香蕉国产精品| 国产视频一区二区在线看| 欧美一级毛片孕妇| 亚洲无线在线观看| 精品国产超薄肉色丝袜足j| 亚洲avbb在线观看| av在线天堂中文字幕| 少妇粗大呻吟视频| 国产人伦9x9x在线观看| 狠狠狠狠99中文字幕| 日韩欧美国产一区二区入口| 淫妇啪啪啪对白视频| 欧美在线一区亚洲| 国产野战对白在线观看| 脱女人内裤的视频| 丝袜人妻中文字幕| 欧美在线黄色| netflix在线观看网站|