• <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| 操美女的视频在线观看| 极品教师在线免费播放| 国产深夜福利视频在线观看| 一二三四在线观看免费中文在| 日本撒尿小便嘘嘘汇集6| 18美女黄网站色大片免费观看| 国产免费现黄频在线看| 亚洲成人久久性| 中文字幕色久视频| 最近最新免费中文字幕在线| 久9热在线精品视频| 久久精品国产99精品国产亚洲性色 | 欧美午夜高清在线| 777久久人妻少妇嫩草av网站| 国产免费现黄频在线看| 91成人精品电影| 伊人久久大香线蕉亚洲五| 免费看a级黄色片| 国产又色又爽无遮挡免费看| av在线播放免费不卡| 宅男免费午夜| 亚洲男人天堂网一区| 国产一区二区三区综合在线观看| 天天添夜夜摸| 国产伦人伦偷精品视频| 国产精品成人在线| 十分钟在线观看高清视频www| 午夜久久久在线观看| 亚洲在线自拍视频| 999精品在线视频| 在线免费观看的www视频| 最近最新免费中文字幕在线| 亚洲情色 制服丝袜| 国产伦一二天堂av在线观看| 好看av亚洲va欧美ⅴa在| 精品福利永久在线观看| 黄片大片在线免费观看| 99久久国产精品久久久| 99久久精品国产亚洲精品| 两人在一起打扑克的视频| 在线av久久热| 午夜老司机福利片| 自拍欧美九色日韩亚洲蝌蚪91| www.精华液| 99精品在免费线老司机午夜| 欧美久久黑人一区二区| 脱女人内裤的视频| 97碰自拍视频| 免费女性裸体啪啪无遮挡网站| 99久久精品国产亚洲精品| 亚洲av片天天在线观看| 色在线成人网| 亚洲av熟女| 性欧美人与动物交配| 一本综合久久免费| 国产在线观看jvid| 高清黄色对白视频在线免费看| 大码成人一级视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品久久成人aⅴ小说| 69av精品久久久久久| 人妻久久中文字幕网| 免费在线观看视频国产中文字幕亚洲| av在线播放免费不卡| 电影成人av| 午夜成年电影在线免费观看| 成人特级黄色片久久久久久久| 精品乱码久久久久久99久播| 国产精品成人在线| 50天的宝宝边吃奶边哭怎么回事| 成人国产一区最新在线观看| 一夜夜www| 久久精品国产99精品国产亚洲性色 | 久久性视频一级片| 麻豆一二三区av精品| 日本免费一区二区三区高清不卡 | 国产极品粉嫩免费观看在线| 欧美日韩国产mv在线观看视频| 国产av在哪里看| 香蕉丝袜av| 久久久国产成人精品二区 | 中亚洲国语对白在线视频| 国产av一区在线观看免费| 亚洲五月色婷婷综合| 91麻豆精品激情在线观看国产 | 18禁观看日本| 久久中文字幕一级| 国产av一区在线观看免费| 国产视频一区二区在线看| 真人做人爱边吃奶动态| 男人舔女人的私密视频| 国产亚洲精品久久久久5区| 自拍欧美九色日韩亚洲蝌蚪91| 9热在线视频观看99| 久热这里只有精品99| 日韩大码丰满熟妇| 窝窝影院91人妻| 大型黄色视频在线免费观看| 一夜夜www| 国产成人啪精品午夜网站| 侵犯人妻中文字幕一二三四区| 精品乱码久久久久久99久播| 亚洲欧美激情综合另类| 欧美 亚洲 国产 日韩一| 天天影视国产精品| 大香蕉久久成人网| 窝窝影院91人妻| 亚洲精品在线观看二区| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 一二三四在线观看免费中文在| 欧美黄色片欧美黄色片| 日韩一卡2卡3卡4卡2021年| 精品久久久精品久久久| 午夜a级毛片| 日韩高清综合在线| 麻豆一二三区av精品| 亚洲自拍偷在线| 在线看a的网站| 国产片内射在线| 热re99久久精品国产66热6| 99在线人妻在线中文字幕| 女警被强在线播放| 欧美乱色亚洲激情| 免费在线观看完整版高清| 别揉我奶头~嗯~啊~动态视频| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 一边摸一边抽搐一进一小说| 成人精品一区二区免费| 男人舔女人下体高潮全视频| 亚洲五月婷婷丁香| 成人国语在线视频| 波多野结衣一区麻豆| 在线免费观看的www视频| 精品免费久久久久久久清纯| 一进一出抽搐动态| 亚洲avbb在线观看| 欧美最黄视频在线播放免费 | 午夜影院日韩av| 色综合站精品国产| 青草久久国产| 搡老乐熟女国产| 欧美日韩精品网址| 村上凉子中文字幕在线| 啦啦啦 在线观看视频| 日韩高清综合在线| 免费久久久久久久精品成人欧美视频| 欧美乱色亚洲激情| av天堂在线播放| 国产亚洲欧美在线一区二区| 青草久久国产| 成人18禁高潮啪啪吃奶动态图| 日韩精品中文字幕看吧| 久久人人精品亚洲av| 久久中文字幕人妻熟女| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点 | 1024香蕉在线观看| 午夜免费激情av| 成人18禁高潮啪啪吃奶动态图| 在线观看www视频免费| 日本黄色日本黄色录像| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人添人人爽欧美一区卜| 国产亚洲av高清不卡| 免费av毛片视频| av超薄肉色丝袜交足视频| 亚洲第一欧美日韩一区二区三区| 免费日韩欧美在线观看| 亚洲国产看品久久| 极品教师在线免费播放| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲| 国产真人三级小视频在线观看| 满18在线观看网站| 日韩欧美在线二视频| 精品福利观看| 欧美午夜高清在线| 国产av又大| 精品国产乱子伦一区二区三区| 国产亚洲精品久久久久久毛片| 国产精品香港三级国产av潘金莲| 黄频高清免费视频| 免费在线观看黄色视频的| 怎么达到女性高潮| 侵犯人妻中文字幕一二三四区| 男女做爰动态图高潮gif福利片 | 国产一卡二卡三卡精品| 国产成人精品无人区| 久久久国产成人精品二区 | 五月开心婷婷网| 中国美女看黄片| 久久久国产成人精品二区 | 天天添夜夜摸| 夜夜夜夜夜久久久久| 看黄色毛片网站| 一进一出抽搐gif免费好疼 | 久久国产精品人妻蜜桃| 水蜜桃什么品种好| 亚洲欧美日韩无卡精品| 淫秽高清视频在线观看| 很黄的视频免费| 看片在线看免费视频| 女人被躁到高潮嗷嗷叫费观| 99精品欧美一区二区三区四区| 大型av网站在线播放| 亚洲第一欧美日韩一区二区三区| 香蕉久久夜色| 久久国产乱子伦精品免费另类| 国产精品av久久久久免费| 色婷婷av一区二区三区视频| 丝袜美足系列| 天天躁狠狠躁夜夜躁狠狠躁| netflix在线观看网站| 国产高清视频在线播放一区| 亚洲激情在线av| 亚洲一区二区三区欧美精品| 久久午夜综合久久蜜桃| 日韩中文字幕欧美一区二区| 免费少妇av软件| 亚洲熟妇熟女久久| 丝袜人妻中文字幕| 久久国产精品人妻蜜桃| 精品国产亚洲在线| 一二三四社区在线视频社区8| 在线观看免费日韩欧美大片| 国产在线精品亚洲第一网站| 亚洲一区二区三区色噜噜 | a在线观看视频网站| 91老司机精品| 身体一侧抽搐| videosex国产| 欧美丝袜亚洲另类 | 成人免费观看视频高清| 日韩有码中文字幕| 日日夜夜操网爽| 精品电影一区二区在线| 亚洲久久久国产精品| 中文欧美无线码| 精品国产乱子伦一区二区三区| 久久久久久免费高清国产稀缺| 免费高清视频大片| 欧美日韩福利视频一区二区| 成人精品一区二区免费| 午夜影院日韩av| 亚洲国产精品一区二区三区在线| 日韩欧美在线二视频| 少妇 在线观看| 国产欧美日韩一区二区三| 老鸭窝网址在线观看| 国产成人欧美在线观看| 免费在线观看视频国产中文字幕亚洲| 天天影视国产精品| 夜夜爽天天搞| 亚洲欧洲精品一区二区精品久久久| 国产视频一区二区在线看| 欧美午夜高清在线| 别揉我奶头~嗯~啊~动态视频| 在线观看舔阴道视频| 在线视频色国产色| 成在线人永久免费视频| 桃色一区二区三区在线观看| 国产精品国产高清国产av| 他把我摸到了高潮在线观看| 制服诱惑二区| 久久精品国产综合久久久| 国产高清激情床上av| www.www免费av| 精品国产乱子伦一区二区三区| 天天躁夜夜躁狠狠躁躁| 精品国产乱码久久久久久男人| 欧美日韩亚洲国产一区二区在线观看| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 制服诱惑二区| 亚洲欧美一区二区三区久久| 国产三级在线视频| 欧美成人性av电影在线观看| 国产aⅴ精品一区二区三区波| 三上悠亚av全集在线观看| 亚洲av电影在线进入| 国产亚洲欧美精品永久| 欧美精品亚洲一区二区| 91精品三级在线观看| 午夜老司机福利片| 美女国产高潮福利片在线看| 他把我摸到了高潮在线观看| 中文字幕人妻熟女乱码| 亚洲午夜理论影院| 亚洲国产毛片av蜜桃av| a级毛片在线看网站| 黄片播放在线免费| 乱人伦中国视频| 亚洲欧美一区二区三区黑人| 在线播放国产精品三级| 久久婷婷成人综合色麻豆| 精品久久久久久久毛片微露脸| 人妻久久中文字幕网| 午夜两性在线视频| 美女高潮到喷水免费观看| x7x7x7水蜜桃| 老熟妇乱子伦视频在线观看| 免费女性裸体啪啪无遮挡网站| 成人精品一区二区免费| 国产精品日韩av在线免费观看 | 一边摸一边抽搐一进一出视频| 久久久精品欧美日韩精品| 我的亚洲天堂| 一二三四社区在线视频社区8| av天堂久久9| 国产又色又爽无遮挡免费看| 高清av免费在线| 波多野结衣av一区二区av| 99在线视频只有这里精品首页| 亚洲熟女毛片儿| 日韩av在线大香蕉| 精品福利观看| 老鸭窝网址在线观看| 亚洲欧美激情综合另类| 亚洲av片天天在线观看| 国产在线观看jvid| 国产精品成人在线| 法律面前人人平等表现在哪些方面| 国产亚洲精品一区二区www| 亚洲美女黄片视频| 亚洲精品粉嫩美女一区| 女警被强在线播放| 又黄又爽又免费观看的视频| 午夜福利,免费看| 操美女的视频在线观看| 国产视频一区二区在线看| 久久久国产精品麻豆| 精品欧美一区二区三区在线| 国产视频一区二区在线看| 久久精品人人爽人人爽视色| 成人亚洲精品av一区二区 | 制服诱惑二区| 欧美精品啪啪一区二区三区| 欧美黑人欧美精品刺激| 国产成人欧美在线观看| 天堂俺去俺来也www色官网| 黄色视频,在线免费观看| 久久久国产一区二区| 一区二区三区激情视频| 亚洲人成77777在线视频| 午夜成年电影在线免费观看| 亚洲精品国产区一区二| 精品免费久久久久久久清纯| 国产在线精品亚洲第一网站| 国产野战对白在线观看| 丝袜人妻中文字幕| 久久精品国产综合久久久| 国产99久久九九免费精品| 一级a爱片免费观看的视频| 久久午夜亚洲精品久久| 午夜日韩欧美国产| 午夜影院日韩av| 欧美乱妇无乱码| 一边摸一边做爽爽视频免费| 在线观看免费日韩欧美大片| 999精品在线视频| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 亚洲五月婷婷丁香| 黄色 视频免费看| 久久午夜亚洲精品久久| 午夜视频精品福利| 国产欧美日韩综合在线一区二区| www.自偷自拍.com| 精品无人区乱码1区二区| 国产欧美日韩一区二区三区在线| 久久久精品国产亚洲av高清涩受| 新久久久久国产一级毛片| 两人在一起打扑克的视频| 亚洲少妇的诱惑av| 性欧美人与动物交配| 天天躁狠狠躁夜夜躁狠狠躁| 高清av免费在线| 国产三级在线视频| 欧美色视频一区免费| 国产精品国产高清国产av| 黄色视频不卡| 一级片'在线观看视频| 亚洲成人久久性| 免费搜索国产男女视频| 日日摸夜夜添夜夜添小说| 久久亚洲精品不卡| 性少妇av在线| 亚洲国产毛片av蜜桃av| 免费久久久久久久精品成人欧美视频| 免费av毛片视频| 国产aⅴ精品一区二区三区波| 亚洲av美国av| 丁香六月欧美| av在线天堂中文字幕 | 欧美老熟妇乱子伦牲交| 久久热在线av| 长腿黑丝高跟| 亚洲午夜理论影院| 最近最新中文字幕大全免费视频| 国产精品98久久久久久宅男小说| 国产精华一区二区三区| 日本免费一区二区三区高清不卡 | 丝袜美足系列| 性色av乱码一区二区三区2| 色婷婷久久久亚洲欧美| 一级片免费观看大全| 手机成人av网站| 免费搜索国产男女视频| 一区二区三区激情视频| 久99久视频精品免费| 欧美+亚洲+日韩+国产| 99久久精品国产亚洲精品| 亚洲自拍偷在线| 精品日产1卡2卡| 国产av一区二区精品久久| 国产片内射在线| 岛国视频午夜一区免费看| 免费在线观看日本一区| 嫩草影院精品99| 亚洲自偷自拍图片 自拍| 亚洲av熟女| 少妇 在线观看| 国产精品偷伦视频观看了| 亚洲国产欧美日韩在线播放| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 91成年电影在线观看| 午夜免费鲁丝| 一区福利在线观看| 在线av久久热| 欧美+亚洲+日韩+国产| 国产精华一区二区三区| 一边摸一边抽搐一进一出视频| 久久人妻福利社区极品人妻图片| 国产激情欧美一区二区| 黑丝袜美女国产一区| 不卡av一区二区三区| av超薄肉色丝袜交足视频| 成人特级黄色片久久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 涩涩av久久男人的天堂| 亚洲男人的天堂狠狠| 国产精品免费一区二区三区在线| 国产1区2区3区精品| 亚洲九九香蕉| 国产精品一区二区精品视频观看| 黄频高清免费视频| 成人精品一区二区免费| 十八禁网站免费在线| 亚洲欧洲精品一区二区精品久久久| 丁香欧美五月| e午夜精品久久久久久久| 少妇粗大呻吟视频| 亚洲人成电影观看| av免费在线观看网站| 精品免费久久久久久久清纯| 久久久久久久久中文| 欧美激情 高清一区二区三区| 亚洲avbb在线观看| 亚洲国产精品sss在线观看 | 中文亚洲av片在线观看爽| 人成视频在线观看免费观看| 丰满饥渴人妻一区二区三| 在线观看日韩欧美| 国产男靠女视频免费网站| 校园春色视频在线观看| 欧美黑人欧美精品刺激| 精品国产国语对白av| 性色av乱码一区二区三区2| 国产精品免费视频内射| 中文字幕人妻熟女乱码| 久久国产乱子伦精品免费另类| 日韩欧美在线二视频| 深夜精品福利| 亚洲av成人一区二区三| 免费搜索国产男女视频| 黄色 视频免费看| 夜夜夜夜夜久久久久| 国产高清激情床上av| 高潮久久久久久久久久久不卡| videosex国产| 美女国产高潮福利片在线看| 国产av又大| 欧美日本亚洲视频在线播放| 国产av一区在线观看免费| 啦啦啦 在线观看视频| 一本大道久久a久久精品| 18美女黄网站色大片免费观看| 黄色 视频免费看| 国产亚洲欧美98| 80岁老熟妇乱子伦牲交| 丁香六月欧美| 男女下面插进去视频免费观看| 在线观看www视频免费| 99久久国产精品久久久| 亚洲片人在线观看| 精品国内亚洲2022精品成人| 少妇裸体淫交视频免费看高清 | 国产亚洲精品久久久久5区| 麻豆成人av在线观看| 久99久视频精品免费| 一级a爱片免费观看的视频| 精品第一国产精品| 露出奶头的视频| 又黄又爽又免费观看的视频| 黄片大片在线免费观看| 男女午夜视频在线观看| 国产精品一区二区在线不卡| 久久久精品欧美日韩精品| 黄色成人免费大全| 三级毛片av免费| 久久狼人影院| 久久青草综合色| 欧美日韩乱码在线| 久久久久久免费高清国产稀缺| e午夜精品久久久久久久| 免费av毛片视频| 久久久久久久久中文| 母亲3免费完整高清在线观看| 色哟哟哟哟哟哟| 成年女人毛片免费观看观看9| 亚洲av美国av| 国产免费男女视频| 香蕉国产在线看| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久久久免费视频 | 久久 成人 亚洲| 日本a在线网址| 999精品在线视频| 久久草成人影院| 天堂动漫精品| 欧美精品啪啪一区二区三区| 午夜精品久久久久久毛片777| 国产真人三级小视频在线观看| 人妻丰满熟妇av一区二区三区| 叶爱在线成人免费视频播放| 精品国产亚洲在线| 精品久久蜜臀av无| 夜夜爽天天搞| 大型黄色视频在线免费观看| 色尼玛亚洲综合影院| 亚洲中文av在线| 九色亚洲精品在线播放| 精品国产一区二区久久| 免费av毛片视频| 韩国av一区二区三区四区| 麻豆av在线久日| 97碰自拍视频| 国产高清国产精品国产三级| 大陆偷拍与自拍| 亚洲三区欧美一区| 国产成人欧美在线观看| 亚洲激情在线av| av国产精品久久久久影院| 欧美激情极品国产一区二区三区| 久久久久久久久中文| 欧美日韩乱码在线| 97人妻天天添夜夜摸| 大香蕉久久成人网| 国产欧美日韩一区二区三区在线| 免费在线观看亚洲国产| 村上凉子中文字幕在线| 亚洲欧美精品综合久久99| 人人妻,人人澡人人爽秒播| 两人在一起打扑克的视频| 国产成人系列免费观看| 午夜成年电影在线免费观看| 丰满的人妻完整版| 免费在线观看日本一区| 亚洲国产精品一区二区三区在线| 亚洲中文av在线| 超色免费av| 99国产精品一区二区三区| 久久精品国产清高在天天线| 在线永久观看黄色视频| 国产av又大| 亚洲专区国产一区二区| 一个人免费在线观看的高清视频| 成人亚洲精品av一区二区 | 亚洲性夜色夜夜综合| 亚洲精品国产色婷婷电影| 国产野战对白在线观看| 99精品久久久久人妻精品| 水蜜桃什么品种好| а√天堂www在线а√下载| 我的亚洲天堂| 国产蜜桃级精品一区二区三区| 在线av久久热| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 曰老女人黄片| 欧美日韩乱码在线| 亚洲午夜理论影院| 琪琪午夜伦伦电影理论片6080| 久久影院123| 亚洲精品久久成人aⅴ小说| 欧美黄色淫秽网站| 亚洲精品国产区一区二| 亚洲七黄色美女视频| 日韩免费高清中文字幕av| 欧美日本中文国产一区发布| 国内久久婷婷六月综合欲色啪| 一级a爱片免费观看的视频| 女人被狂操c到高潮| 婷婷精品国产亚洲av在线| 视频区图区小说| 不卡一级毛片| 嫁个100分男人电影在线观看| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 亚洲一区二区三区不卡视频|