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

    木林子大樣地兩個(gè)木姜子屬物種空間分布格局及關(guān)聯(lián)性

    2021-09-13 05:01:11黃小姚蘭艾訓(xùn)儒朱江朱強(qiáng)王進(jìn)吳漫玲劉松伯
    廣西植物 2021年8期

    黃小 姚蘭 艾訓(xùn)儒 朱江 朱強(qiáng) 王進(jìn) 吳漫玲 劉松伯

    摘 要:該文以湖北木林子國(guó)家級(jí)自然保護(hù)區(qū)15 hm2大樣地同屬物種 [木姜子(Litsea pungens)和黃丹木姜子(L. elongata)]為研究對(duì)象,通過(guò)成對(duì)相關(guān)函數(shù)和不同零模型(完全空間隨機(jī)模型、異質(zhì)泊松模型和先決條件模型)的方法分析兩個(gè)物種的空間分布格局、種內(nèi)和種間空間關(guān)聯(lián)性,從同屬物種的空間分布格局角度探討群落的構(gòu)建機(jī)制。結(jié)果表明:(1)種群空間分布格局研究中,兩個(gè)種的空間分布格局大致為聚集分布?;贑SR零模型,木姜子和不同徑級(jí)(幼樹(shù)、小樹(shù)、中樹(shù)和成年樹(shù))的聚集尺度分別為0~114 m、0~103 m、0~56 m、0~34 m和0~16 m,黃丹木姜子和不同徑級(jí)的聚集尺度分別為0~150 m、0~150 m、0~59 m、0~56 m和0~14 m;基于HP零模型,木姜子整體的空間分布格局聚集尺度為0~23 m,其中幼樹(shù)、小樹(shù)和中樹(shù)均在0~22 m呈現(xiàn)聚集分布,成年樹(shù)僅在0~3 m尺度上表現(xiàn)為聚集分布,黃丹木姜子整體的空間分布格局聚集尺度為0~36 m、140~150 m,各徑級(jí)的聚集尺度分別為0~35 m、0~35 m、0~12 m和0~6 m。(2)種內(nèi)關(guān)聯(lián)性研究中,木姜子和黃丹木姜子不同徑級(jí)個(gè)體種內(nèi)關(guān)聯(lián)在小尺度上基本呈現(xiàn)正關(guān)聯(lián),隨著尺度的增大變?yōu)槲⑷醯臒o(wú)關(guān)聯(lián)或負(fù)關(guān)聯(lián)。(3)種間關(guān)聯(lián)研究中,木姜子對(duì)黃丹木姜子和黃丹木姜子對(duì)木姜子的空間關(guān)系大致相同。總體上,兩物種之間均體現(xiàn)為在30 m尺度內(nèi)對(duì)另外一個(gè)物種呈現(xiàn)正關(guān)聯(lián),40~68 m尺度內(nèi)呈現(xiàn)負(fù)關(guān)聯(lián);不同徑級(jí)之間都表現(xiàn)為無(wú)關(guān)聯(lián),偶爾出現(xiàn)正關(guān)聯(lián)或負(fù)關(guān)聯(lián)。研究認(rèn)為,種群空間分布格局基本上為聚集分布,其主要受到生境異質(zhì)性和擴(kuò)散限制的影響,木林子大樣地中木姜子和黃丹木姜子空間分布格局存在相似性,而不同的生境偏好可能是導(dǎo)致兩同屬物種共存的原因。

    關(guān)鍵詞:同屬物種,點(diǎn)格局,空間分布格局,空間關(guān)聯(lián),物種共存

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

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

    文章編號(hào):1000-3142(2021)08-1333-13

    Abstract: In this paper, species of the same congeneric species (Litsea pungens and L. elongata ) in large sample plot of 15 hm2 in Mulinzi National Nature Reserve of Hubei Province were taken as the study object. The spatial distribution pattern, intraspecies and interspecific spatial correlations of the two species were analyzed by pairwise correlation function and different zero models(complete spatial randomness, heterogeneous Poisson and antecedent condition), and the construction mechanism of the community was discussed from the point of view of the spatial distribution pattern of the same congeneric species. The results were as follows: (1) In the study of population spatial distribution pattern, the spatial distribution pattern of the two species was generally aggregation distribution. The aggregation scales on the basis of the CSR zero model, the L. pungens and the different diameter stages(sapling, small tree, medium tree, adult tree) were 0-114 m, 0-103 m, 0-56 m, 0-34 m, and 0-16 m, respectively, and the aggregation scales of the L. elongata and the different diameter stages were 0-150 m, 0-150 m and 0-59 m, 0-56 m and 0-14 m; On the basis of the HP-zero model, the spatial distribution pattern of the whole L. pungens was 0-23 m, where in all of the sapling, small and medium trees were 0-22 m, and the aggregation distribution of adult trees was only on the scale of 0-3 m. The aggregation scale of the whole spatial distribution pattern of L. elongata was 0-36 m and 140-150 m, and the aggregation scale of each diameter grade was 0-35 m, 0-35 m, 0-12 m and 0-6 m. (2) In the study of intraspecific correlation, the intraspecies correlation of different diameter individuals showed positive correlation on a small scale, and became weakly unrelated or negative with the increase of scale. (3) In the study of interspecific correlation, the spatial relationship between the two species was about the same. On the whole, the two species showed positive correlation with another species within 30 m scale and negative correlation within 40-68 m scale, and there was no correlation between different diameter grades, and there was positive or negative correlation between the two species. It is considered that the spatial distribution pattern of the population is basically aggregation distribution, which is mainly affected by habitat heterogeneity and diffusion. The spatial distribution patterns of L. pungens and L. elongata in large sample plot are similar, and different habitat preferences may be the cause of the coexistence of the two species.

    Key words: same congeneric species, point pattern, spatial distribution pattern, spatial correlation, species coexistence

    一般來(lái)說(shuō),同屬物種擁有相同祖先,在形態(tài)學(xué)和生理特征上都具有一定的相似性,所需的資源可能相似,從生態(tài)位分化來(lái)看,同屬物種可能因?yàn)檫@些相似性而產(chǎn)生競(jìng)爭(zhēng),而不能穩(wěn)定共存(Ackerly & Donoghue, 1998; Webb & Ackerly, 2002),同屬物種共存現(xiàn)象多發(fā)現(xiàn)于熱帶和溫帶植物群落中(張健, 2009)。因此,生態(tài)學(xué)家一直對(duì)同屬植物的共存進(jìn)行激烈的探討。在生態(tài)研究中對(duì)同屬物種進(jìn)行比較,可以最大限度地減少系統(tǒng)差異造成的混亂影響,從同屬物種的空間結(jié)構(gòu)和相關(guān)性的角度進(jìn)行探討(Swenson et al., 2006; 姚蘭, 2016)。

    植物種群的空間分布格局是指某一種群在一定水平內(nèi)所有個(gè)體的分布位置或配置(曲仲湘, 1984; 張金屯, 1995),到現(xiàn)在為止是生態(tài)學(xué)研究的熱點(diǎn)之一(Wiegleb, 1989; Nathan, 2006)。點(diǎn)格局方法由Ripley提出,它能夠很好地解決傳統(tǒng)格局分析方法的缺點(diǎn)(Condit, 2000; 張金屯和孟東平, 2004)。近些年,點(diǎn)格局分析越來(lái)越廣泛,主要側(cè)重于單物種的空間分布格局(范繁榮等, 2008; 黃小波等, 2018)、優(yōu)勢(shì)樹(shù)種的空間分布格局(趙中華等, 2011; 張玉雙等, 2016)、兩種或兩種以上物種的空間相關(guān)性(姚良錦等, 2018)以及空間分布格局與環(huán)境之間關(guān)系(劉旻霞等, 2018; 閆海冰, 2018)等研究較多,對(duì)于種群同科(何增麗等, 2017)或同屬(陳俊等, 2018)研究較少。

    樟科(Lauraceae)木姜子屬(Litsea Lam) 的植物分布非常廣泛,在中國(guó)大約有72種(陳幼竹, 2004)。亞熱帶常綠落葉闊葉混交林保存最為完好的地區(qū)包括木林子國(guó)家級(jí)自然保護(hù)區(qū)(黃永濤, 2016),其中木姜子屬植物僅有2種:黃丹木姜子和木姜子,其數(shù)量占樣地內(nèi)物種總數(shù)量的4.8%。木姜子的藥用價(jià)值很高,黃丹木姜子具有較高的商業(yè)價(jià)值(陳幼竹, 2004)。木姜子和黃丹木姜子的研究?jī)H僅只限于藥用價(jià)值等方面的研究,黃丹木姜子還見(jiàn)于葉凋落量及養(yǎng)分特征研究(董舒, 2016),對(duì)于更多的生態(tài)學(xué)方面的研究未見(jiàn)有報(bào)道,黃丹木姜子和木姜子雖不是建群種和先鋒種,這類小喬木在種群空間分布格局的研究同樣具有重要意義。本研究以湖北木林子大樣地中同屬的木姜子和黃丹木姜子為研究對(duì)象,采用g(r)函數(shù)對(duì)木姜子屬兩個(gè)種進(jìn)行空間格局及空間關(guān)聯(lián)分析,探討其成因,探討兩物種是否存在相似性,同屬物種共存機(jī)制,并為研究亞熱帶地區(qū)物種共存機(jī)制提供參考。

    1 研究區(qū)概況

    研究區(qū)(109°59′30″—110°17′58″ E、29°55′59″—30°10′47″ N)位于湖北省鶴峰縣,面積為20 838 hm2,屬于云貴高原的東北延伸區(qū)。該地區(qū)最低海拔為1 100 m,最高海拔為2 095.6 m,大陸性季風(fēng)氣候,冬季無(wú)嚴(yán)寒,夏季無(wú)極端炎熱,雨熱同季,日照充足,山地性氣候明顯,年均溫15.5 ℃,年均降水量1 700~1 900 mm,年均相對(duì)濕度82%,無(wú)霜期270~279 d。土壤為黃棕壤、棕壤和黃壤。

    木林子15 hm2大樣地位于保護(hù)區(qū)核心區(qū)內(nèi),海拔為1 588~1 780 m,地勢(shì)相對(duì)緩慢,屬于亞熱帶山地常綠落葉闊葉混交林,主要樹(shù)種有翅柃(Eurya alata)、多脈青岡(Cyclobalanopsis multinervis)、川陜鵝耳櫪(C. fargesiana)、小葉青岡(C. myrsinifolia)、石灰花楸(Sorbus folgneri)、四照花(Dendrobenthamia japonica var. chinensis)、錐栗(Castanea henryi)、白櫟(Quercus fabri)、山礬(Symplocos sumuntia)和光葉山礬(S. lancifolia)等(姚蘭, 2016)。

    2 研究方法

    2.1 樣地建設(shè)及野外調(diào)查

    2013年,按照CTFS大樣地建設(shè)與檢測(cè)技術(shù)規(guī)范,建立15 hm2固定監(jiān)測(cè)樣地,分成375塊20 m × 20 m小樣地,用防腐處理的不銹鋼管將4角作為永久標(biāo)識(shí),每個(gè)小樣地再分為4個(gè)10 m × 10 m樣方,每個(gè)樣方再分為4個(gè)5 m × 5 m的小樣方,對(duì)DBH≥1.0 cm的木本植物進(jìn)行掛牌,并用紅色油漆標(biāo)記(樹(shù)高1.3 m),對(duì)DBH≥1 cm的所有木本植物進(jìn)行每木檢尺,識(shí)別物種(科、屬、種),調(diào)查所有木本植物在樣地內(nèi)的坐標(biāo)值(以大樣地單元西南角為原點(diǎn),測(cè)定x軸和y軸的坐標(biāo)值)。

    2.2 種群徑級(jí)劃分

    運(yùn)用徑級(jí)代替齡級(jí)的方法,結(jié)合湖北木林子樣地內(nèi)木姜子和黃丹木姜子的徑級(jí)特征,將胸徑DBH劃分為4個(gè)徑級(jí):Ⅰ級(jí)(1 cm≤DBH<2.5 cm),Ⅱ級(jí)(2.5 cm≤DBH<5 cm),Ⅲ級(jí)(5 cm≤DBH<10 cm),Ⅳ級(jí)(DBH≥10 cm),分別為幼樹(shù)、小樹(shù)、中樹(shù)和成年樹(shù)。

    2.3 點(diǎn)格局分析方法

    本研究采用點(diǎn)格局分析中成對(duì)相關(guān)函數(shù)g(r)進(jìn)行分析。g(r)函數(shù)是基于對(duì)所有定位的成對(duì)個(gè)體之間距離的格局分析,消除了Ripieys K(r)函數(shù)在大尺度上產(chǎn)生的累積效應(yīng)(Ripley, 1976; Thorsten et al., 2007)。分析時(shí),選取距離尺度為較短邊長(zhǎng)的1/2(150 m),圓環(huán)寬度為1 m,通過(guò)199次Monte Carlo隨機(jī)模擬,利用模擬的最大值和最小值繪制上下包跡線,計(jì)算99%的置信區(qū)間(Wiegand et al., 2006, 2009)。

    2.4 零模型選擇

    采用點(diǎn)格局分析,需要檢測(cè)生境異質(zhì)性在物種分布中是否存在影響。一般來(lái)說(shuō),如果在大于10 m尺度上呈聚集分布,那就可以解釋為受到生境異質(zhì)性的影響;采用異質(zhì)泊松模型和先決條件模型目的在于消除大尺度上生境異質(zhì)性對(duì)空間分布格局的影響,在排除生境異質(zhì)性后,如果小尺度的格局還是聚集,則說(shuō)明擴(kuò)散限制存在作用。所以在進(jìn)行點(diǎn)格局分析時(shí),不同的分析需要選擇不同的零假設(shè)模型(null model),選取正確合適的零模型十分重要(梁爽等, 2014)。在本研究中單變量g(r)函數(shù)采用完全空間隨機(jī)模型(complete spatial randomness, CSR)和異質(zhì)泊松模型(heterogeneous Poisson, HP),雙變量g(r)函數(shù)選擇先決條件模型(antecedent condition, AC)。

    文中數(shù)據(jù)的處理與分析使用R 3.4.1軟件中的“spatstat”包中完成點(diǎn)格局分析,使用Sigmaplot 13軟件繪制空間分布點(diǎn),Origin 2018軟件繪制點(diǎn)格局分析圖。

    3 結(jié)果與分析

    3.1 兩個(gè)物種空間分布格局

    整個(gè)樣地中(表1),木姜子1 701株,黃丹木姜子2 306株。由圖1可知,兩個(gè)物種的空間分布總體差異較大。木姜子主要分布于低海拔、山脊地區(qū);黃丹木姜子主要分布于高海拔和山谷兩側(cè)區(qū)域,小部分零星分布于樣地內(nèi)。兩個(gè)物種的整體空間分布格局的研究(圖2)表明,基于CSR零模型:木姜子和黃丹木姜子的聚集尺度分別為0~114 m和0~150 m?;贖P零模型:木姜子的聚集尺度為0~23 m,其他尺度上為微弱的隨機(jī)或均勻分布;黃丹木姜子的聚集尺度為0~36 m、140~150 m,隨機(jī)分布尺度為48~139 m。

    根據(jù)兩個(gè)物種不同徑級(jí)的空間分布格局研究(圖3,圖4)表明,基于CSR零模型:木姜子幼樹(shù)、小樹(shù)和中樹(shù)的聚集尺度分別為0~103 m、0~56 m和0~34 m,成年樹(shù)在0~16 m尺度表現(xiàn)為聚集分布;黃丹木姜子幼樹(shù)、小樹(shù)、中樹(shù)和成年樹(shù)分別在0~150 m、0~59 m、0~56 m和0~14 m尺度上均表現(xiàn)為聚集分布。基于HP零模型:木姜子幼樹(shù)、小樹(shù)和中樹(shù)的聚集尺度均為0~22 m,而成年樹(shù)僅在0~3 m尺度上表現(xiàn)為聚集分布;黃丹木姜子各徑級(jí)的聚集尺度分別為0~35 m、0~35 m、0~12 m和0~6 m。由此可見(jiàn),隨著徑級(jí)的增長(zhǎng),木姜子與黃丹木姜子均呈現(xiàn)聚集尺度逐漸減小的現(xiàn)象,特別是成年樹(shù)僅在小尺度上表現(xiàn)為聚集分布,而隨機(jī)分布格局現(xiàn)象比較明顯。

    3.2 兩個(gè)物種不同徑級(jí)種內(nèi)關(guān)聯(lián)性

    木姜子不同徑級(jí)種內(nèi)關(guān)聯(lián)見(jiàn)圖5?;贏C零模型,成年樹(shù)對(duì)幼樹(shù)在所有尺度上均表現(xiàn)為無(wú)關(guān)聯(lián),而其他不同徑級(jí)之間在較小尺度(r≤20 m)上表現(xiàn)為正相關(guān),其中小樹(shù)對(duì)幼樹(shù)、中樹(shù)對(duì)幼樹(shù)、中樹(shù)對(duì)小樹(shù)、成年樹(shù)對(duì)小樹(shù)、成年樹(shù)對(duì)中樹(shù)正關(guān)聯(lián)尺度分別為0~20 m、0~18 m、0~20 m、0~14 m和0~17 m。

    黃丹木姜子不同徑級(jí)種內(nèi)關(guān)聯(lián)見(jiàn)圖6?;贏C零模型,黃丹木姜子成年樹(shù)對(duì)中樹(shù)在所有尺度上表現(xiàn)為無(wú)相關(guān),小樹(shù)對(duì)幼樹(shù)、中樹(shù)對(duì)幼樹(shù)、成年樹(shù)對(duì)幼樹(shù)、中樹(shù)對(duì)小樹(shù)、成年樹(shù)對(duì)小樹(shù)分別在0~36 m、0~25 m、8~15 m、0~23 m和0~6 m尺度上呈現(xiàn)正關(guān)聯(lián)。

    3.3 兩個(gè)物種種間關(guān)聯(lián)性

    用先決條件模型排除30 m以上生境異質(zhì)性,結(jié)果顯示:整體上,木姜子對(duì)黃丹木姜子在0~30 m呈現(xiàn)正關(guān)聯(lián),黃丹木姜子對(duì)木姜子在7~29 m呈現(xiàn)正關(guān)聯(lián)(圖7)。不同徑級(jí)之間,木姜子不同徑級(jí)對(duì)黃丹木姜子不同徑級(jí)空間關(guān)聯(lián)結(jié)果顯示(圖8):木姜子幼樹(shù)對(duì)黃丹木姜子幼樹(shù)在4~34 m尺度表現(xiàn)為正關(guān)聯(lián);木姜子小樹(shù)對(duì)黃丹木姜子幼樹(shù)在13~29 m尺度表現(xiàn)為正關(guān)聯(lián),42~67 m尺度上為負(fù)關(guān)聯(lián);木姜子中樹(shù)對(duì)黃丹木姜子幼樹(shù)在40~50 m尺度上為負(fù)關(guān)聯(lián);木姜子小樹(shù)對(duì)黃丹木姜子小樹(shù)在0~25 m尺度上為正關(guān)聯(lián),44~60 m尺度上為負(fù)關(guān)聯(lián);木姜子中樹(shù)對(duì)黃丹木姜子中樹(shù)在0~5 m尺度上為正關(guān)聯(lián),在47~54 m和119~128 m尺度上表現(xiàn)為負(fù)關(guān)聯(lián);其余徑級(jí)之間均呈現(xiàn)無(wú)關(guān)聯(lián)。黃丹木姜子不同徑級(jí)對(duì)木姜子不同徑級(jí)空間結(jié)果顯示(圖9):黃丹木姜子幼樹(shù)對(duì)木姜子幼樹(shù)在13~32 m尺度呈現(xiàn)正關(guān)聯(lián),42~76 m尺度呈現(xiàn)負(fù)關(guān)聯(lián);黃丹木姜子小樹(shù)對(duì)木姜子幼樹(shù)在0~13 m尺度呈現(xiàn)正關(guān)聯(lián);黃丹木姜子小樹(shù)對(duì)木姜子小樹(shù)在0~15 m呈現(xiàn)正關(guān)聯(lián),46~59 m尺度呈現(xiàn)負(fù)關(guān)聯(lián);其余徑級(jí)之間均呈現(xiàn)無(wú)關(guān)聯(lián)。

    由此可見(jiàn),木姜子對(duì)黃丹木姜子空間關(guān)系和黃丹木姜子對(duì)木姜子空間關(guān)系大致相同,僅有部分細(xì)微差別,其中兩物種成年樹(shù)與另一物種所有徑級(jí)均沒(méi)有顯著的生態(tài)學(xué)關(guān)系。

    4 討論

    4.1 兩個(gè)物種空間分布格局形成的原因

    植物自身生物學(xué)特性和環(huán)境因素共同作用于種群的空間分布格局,同一樹(shù)種在不同的發(fā)育階段將呈現(xiàn)不同的空間分布模式(Frost & Rydin, 2000)。本研究中基于CSR和HP零模型兩個(gè)物種在研究尺度上主要表現(xiàn)為聚集分布,隨著尺度增加表現(xiàn)為隨機(jī)分布或均勻分布,結(jié)果與大部分研究一致(祝燕等, 2011; Guo et al., 2013)。在很大程度上,生境的異質(zhì)性被認(rèn)為是植被空間格局的主要原因。通常,如果它的聚集尺度大于10 m,則可以解釋為受到棲息地異質(zhì)性的影響(Takyu et al., 2002; Getzin et al., 2008)。基于CSR零模型,在整體上木姜子和黃丹木姜子分別在0~114 m和0~150 m尺度范圍內(nèi)呈聚集分布,并使用HP零模型排除了30 m以上的生境異質(zhì)性,聚集尺度減小成為隨機(jī)分布,可以推測(cè),兩個(gè)種群空間分布格局受到生境異質(zhì)性影響。從幼樹(shù)、小樹(shù)、中樹(shù)、成年樹(shù)時(shí)間序列來(lái)看,分布格局從聚集分布向隨機(jī)分布發(fā)展趨勢(shì),符合一般規(guī)律(Vargas et al., 2010; 張玉雙等, 2011)。其中木姜子和黃丹木姜子不同徑級(jí)之間均體現(xiàn)為幼樹(shù)聚集尺度>小樹(shù)聚集尺度>中樹(shù)聚集尺度>成年樹(shù)聚集尺度,可見(jiàn)生境異質(zhì)性對(duì)兩個(gè)物種的幼樹(shù)影響最強(qiáng),可能是在此階段,種群的個(gè)體較小,生存力較低,需要通過(guò)聚集分布形成相互庇護(hù),能夠更好地維持種群自身的穩(wěn)定性(張金屯, 1998; Wang et al., 2014)。生境異質(zhì)性對(duì)兩個(gè)物種的成年樹(shù)影響較小,兩個(gè)物種負(fù)密度制約超過(guò)了生境異質(zhì)性對(duì)聚集分布的影響,隨著徑級(jí)的增大,對(duì)資源、對(duì)空間的需求也增加,物種內(nèi)部的競(jìng)爭(zhēng)日益激烈,物種內(nèi)部的稀疏性越來(lái)越強(qiáng),這將導(dǎo)致一些個(gè)體死亡從而形成均勻或隨機(jī)分布的狀態(tài)(舒蘭等, 2019)。

    在小尺度上,由于擴(kuò)散限制可能導(dǎo)致聚集分布。當(dāng)使用HP零模型分析,在較小尺度上兩個(gè)物種均呈現(xiàn)聚集分布,這說(shuō)明擴(kuò)散限制對(duì)兩物種的空間分布格局存在一定的作用。擴(kuò)散限制使種子成熟后會(huì)在重力和風(fēng)等因素的共同作用下落在母樹(shù)周圍,隨著與母樹(shù)距離越遠(yuǎn),種子數(shù)量逐漸減少,從而在小尺度上形成聚集分布,木姜子和黃丹木姜子的果實(shí)均為球果,主要依靠重力傳播為受距離限制。從兩個(gè)物種來(lái)看,木姜子的聚集尺度范圍和聚集程度均小于黃丹木姜子,可能是木姜子屬于小喬木,一般高3~10 m,而黃丹木姜子屬于中喬木,一般高達(dá)12 m,黃丹木姜子比木姜子種子的擴(kuò)散范圍更大,可能是黃丹木姜子在小尺度上的聚集尺度大于木姜子的原因。

    種內(nèi)和種間競(jìng)爭(zhēng)也是影響種群聚集分布的一部分原因,種內(nèi)關(guān)聯(lián)是用來(lái)表示同一物種不同發(fā)育階段在空間分布上的關(guān)聯(lián)性(王磊等, 2010)。本研究發(fā)現(xiàn)兩個(gè)物種之間存在較小種內(nèi)競(jìng)爭(zhēng)和種間競(jìng)爭(zhēng),其中種內(nèi)競(jìng)爭(zhēng)中僅發(fā)現(xiàn)黃丹木姜子的小樹(shù)與幼樹(shù)在一定的尺度上呈現(xiàn)負(fù)相關(guān),說(shuō)明這兩個(gè)徑級(jí)都屬于演替前期,對(duì)資源和環(huán)境的需求趨于相同,種群不穩(wěn)定,處于相互排除的階段,幼樹(shù)的生長(zhǎng)受到了小樹(shù)的抑制。種間競(jìng)爭(zhēng)中,兩個(gè)物種不同徑級(jí)之間大致都趨于無(wú)關(guān)聯(lián)的現(xiàn)象,僅在某些小尺度上呈現(xiàn)微弱的正關(guān)聯(lián)和負(fù)關(guān)聯(lián)的現(xiàn)象,說(shuō)明兩個(gè)物種種間競(jìng)爭(zhēng)較弱,這可能是因?yàn)閮晌锓N不同的生境偏好導(dǎo)致的,研究發(fā)現(xiàn)木姜子個(gè)體更多偏向分布于高海拔和山脊地區(qū),而黃丹木姜子更喜低海拔和山谷兩側(cè)。

    4.2 兩個(gè)同屬物種相似性及共存機(jī)制

    同一個(gè)生境中,同屬物種的空間分布格局應(yīng)該存在相似性。本研究結(jié)果顯示,木姜子和黃丹木姜子總體上和不同徑級(jí)之間的空間分布格局、兩個(gè)物種的種內(nèi)關(guān)聯(lián)以及兩個(gè)物種相互的種間關(guān)系均存在一定的相似性,兩個(gè)物種總體上和不同徑級(jí)之間均表現(xiàn)為聚集分布;在種內(nèi)關(guān)聯(lián)中除了木姜子的成年樹(shù)對(duì)幼樹(shù)和黃丹木姜子的成年樹(shù)對(duì)中樹(shù)表現(xiàn)出不同的空間關(guān)聯(lián),其余的均表現(xiàn)出相似的種內(nèi)空間關(guān)系;在種間關(guān)聯(lián)中木姜子對(duì)黃丹木姜子和黃丹木姜子對(duì)木姜子的種間關(guān)系均表現(xiàn)為幼樹(shù)對(duì)幼樹(shù)、小樹(shù)對(duì)幼樹(shù)、小樹(shù)對(duì)小樹(shù)存在較為強(qiáng)烈的生態(tài)學(xué)關(guān)系,其余徑級(jí)之間生態(tài)學(xué)關(guān)系不顯著。本研究顯示,木姜子和黃丹木姜子在空間分布上存在一定的差異,可能是因?yàn)閮晌锓N不同的生長(zhǎng)環(huán)境導(dǎo)致。

    在一個(gè)區(qū)域內(nèi),多個(gè)物種肯定存在相互關(guān)系才會(huì)產(chǎn)生共存(宗國(guó)等, 2018)。本研究表明,總體上木姜子對(duì)黃丹木姜子和黃丹木姜子對(duì)木姜子的空間關(guān)系均表現(xiàn)為僅在40~68 m尺度內(nèi)為負(fù)關(guān)聯(lián),其余尺度為正關(guān)聯(lián)和無(wú)關(guān)聯(lián),而不同徑級(jí)的木姜子和黃丹木姜子之間在大部分尺度上表現(xiàn)為正關(guān)聯(lián)和無(wú)關(guān)聯(lián),僅在極小尺度上為負(fù)關(guān)聯(lián)??赡苁且?yàn)閮蓚€(gè)物種不同的生境偏好導(dǎo)致,從而減少兩個(gè)物種對(duì)資源的競(jìng)爭(zhēng)。同時(shí),韋博良(2015)研究寶天曼櫟屬的空間分布格局時(shí),發(fā)現(xiàn)櫟屬物種具有較小種間競(jìng)爭(zhēng)關(guān)系,這樣更有利于種群較小的物種能夠共存。本研究顯示,木姜子屬兩個(gè)物種之間存在較小的競(jìng)爭(zhēng)關(guān)系,而是相互促進(jìn)生長(zhǎng),這有利于兩個(gè)物種共存,并同時(shí)抵御其他物種產(chǎn)生競(jìng)爭(zhēng)關(guān)系。

    5 結(jié)論

    綜上所述,木姜子和黃丹木姜子兩個(gè)同屬物種的空間分布存在一定的相似性,生境異質(zhì)性、擴(kuò)散限制和種內(nèi)種間競(jìng)爭(zhēng)都影響兩個(gè)物種空間格局的形成,其中生境異質(zhì)性對(duì)兩個(gè)物種的幼樹(shù)影響最大,擴(kuò)散限制導(dǎo)致黃丹木姜子在排除生境異質(zhì)性后擴(kuò)散范圍大于木姜子,兩物種在種內(nèi)關(guān)系中均體現(xiàn)為在小尺度上呈現(xiàn)正關(guān)聯(lián),同時(shí)兩個(gè)物種之間存在較小的競(jìng)爭(zhēng)關(guān)系;不同的生境偏好是兩個(gè)物種在樣地中實(shí)現(xiàn)共存的重要基礎(chǔ),同時(shí)也可能是因?yàn)橥瑢傥锓N容易形成相互庇護(hù)而抵御其他物種,則需要進(jìn)一步的研究。

    參考文獻(xiàn):

    ACKERLY DD, DONOGHUE MJ, 1998. Leaf size, saplingallometry, and Corners rules: phylogeny and correlated evolution in maples(Acer) [J]. Am Nat, 152(6): 767-791.

    CHEN J, AI XR, YAO L, et al., 2018. Point pattern analysis of two species of Cyclobalan opsis in large plot in Mulinzi Nature Reserve [J]. Sci Silv Sin, 54(10): 1-10. [陳俊, 艾訓(xùn)儒, 姚蘭, 等, 2018. 木林子大樣地兩個(gè)青岡屬優(yōu)勢(shì)種的點(diǎn)格局對(duì)比 [J]. 林業(yè)科學(xué),

    54(10): 1-10.]

    CHEN YZ, 2004. Studies on the botanical origins and qualities of medicinal Litsea Lam [D]. Chengdu: Chengdu University of TCM. [陳幼竹, 2004. 木姜子屬主要藥用植物的品種品質(zhì)研究 [D]. 成都: 成都中醫(yī)藥大學(xué).]

    CONDIT R, 2000. Spatial patterns in the distribution of tropical tree species [J]. Science, 288(5470): 1414-1418.

    DIGGLE PJ, 2003. Statistical analysis of spatial point patterns [M]. 2nd ed. London: Hodder Education Publishers.

    DONG S, 2016. Study on leaf litterfall production and nutrient characteristcs of six dominant species in an evergreen broad-leaved forest in Tiantong [D]. Shanghai: East China Normal University. [董舒, 2016. 天童常綠闊葉林6個(gè)優(yōu)勢(shì)種凋落量及養(yǎng)分特征研究 [D]. 上海: 華東師范大學(xué).]

    FAN FR, PAN BZ, MA XQ, et al., 2008. A study on structure and spatial distribution pattern of Artocarpus hypargyreus population [J]. For Res,(2): 176-181. [范繁榮, 潘標(biāo)志, 馬祥慶, 等, 2008. 白桂木的種群結(jié)構(gòu)和空間分布格局研究 [J]. 林業(yè)科學(xué)研究, (2): 176-181.]

    FROSTI, RYDIN H, 2000. Spatial pattern and size distribution of the animal-dispersed tree Quercus robur in two spruce-dominated forests [J]. coscience, 7(1): 38-44.

    GETZIN S, WIEGAND T, WIEGAND K, et al., 2008. Heterogeneity influences spatial patterns and demographics in forest stands [J]. J Ecol, 96(4): 807-820.

    GUO YL, LU JM, FRANKLIN SB, et al., 2013. Spatial distribution of tree species in a species-rich subtropical mountain forest in central China [J]. Can J For Res, 43(9): 826-835.

    HE ZL, XU H, QIN XS, et al., 2017. Spatial distribution patterns and association of two Apocynaceae plants in the tropical mountain rainforests of Jianfengling, Hainan Island, China [J]. Biodivers Sci, 25(10): 1065-1074. [何增麗, 許涵, 秦新生, 等, 2017. 海南尖峰嶺熱帶山地雨林2種夾竹桃科植物的空間分布格局與關(guān)聯(lián)性 [J]. 生物多樣性, 25(10): 1065-1074.]

    HUANG YT, 2016. The functional community structure and its dynamics of the subtropical evergreen-deciduous broadleaved mixed forest in Southwest Hubei of Central China [D]. Beijing: Chinese Academy of Forestry. [黃永濤, 2016. 亞熱帶常綠落葉闊葉混交林的功能性群落結(jié)構(gòu)及其動(dòng)態(tài)研究 [D]. 北京: 中國(guó)林業(yè)科學(xué)研究院.]

    HUANG XB, LI SF, SU JR, et al., 2018. Distribution of Pinus yunnanensis natural population in Yunlong Tianchi National Nature Reserve [J]. For Res, 31(4): 47-52. [黃小波, 李帥鋒, 蘇建榮, 等, 2018. 云龍?zhí)斐卦颇纤勺匀环N群分布格局分析 [J]. 林業(yè)科學(xué)研究, 31(4): 47-52.]

    LIANG S, XU H, LIN JY, et al., 2014. Spatial distribution pattern of the dominant species Gironniera subaequalis in tropical montane rainforest of Jianfengling, Hainan Island, China [J]. Chin J Plant Ecol, 38(12): 1273-1282. [梁爽, 許涵, 林家怡, 等 , 2014. 尖峰嶺熱帶山地雨林優(yōu)勢(shì)樹(shù)種白顏樹(shù)空間分布格局 [J]. 植物生態(tài)學(xué)報(bào), 38(12): 1273-1282.]

    LIU WX, LIAN YM, LI W, 2018. Effects of microtopography on the pattern and association of dominant population [J]. Chin J Appl Ecol, 29(5): 1569-1575. [劉旻霞, 連依明, 李文, 2018. 微地形對(duì)優(yōu)勢(shì)種群點(diǎn)格局和關(guān)聯(lián)性的影響 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 29(5): 1569-1575.]

    NATHAN R, 2006. Long-distance dispersal of plants [J].

    Science, 5788(3): 786-788.

    QU ZX, 1983. Plant ecology [M]. Beijing: Higher education press. [曲仲湘, 1983. 植物生態(tài)學(xué) [M]. 北京: 高等教育出版社.]

    RIPLEY DB, 1976. The second-order analysis of stationary point processes [J]. J Appl Prob, 13(2): 255-266.

    SHU L, LIU ZG, DONG LB, 2019. Spatial pattern and regeneration characteristics of main woody species in natural secondary forest in Maoershan, Northeast China [J]. Chin J Appl Ecol, 30(6): 1945-1955. [舒蘭, 劉兆剛, 董靈波, 2019. 帽兒山天然次生林內(nèi)主要木本植物空間格局及更新概況 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 30(6): 1945-1955.]

    SWENSON NG, ENQUIST BJ, PITHER J, et al., 2006. The problem and promise of scale dependency in community phylogenetics [J]. Ecology, 87(10): 2418-2424.

    TAKYU M, AIBA SI, KITAYAMA K, 2002. Effects of topography on tropical lowermontane forests under different geological conditions on Mount Kinabalu, Borneo [J]. Plant Ecol, 159(1): 35-49.

    THORSTEN W, SAVITRI G, NIMAL G, et al., 2007. Analyzing the spatial structure of a Sri Lankan tree species with multiple scales of clustering [J]. Ecology, 88(12): 3088-3102.

    VARGAS R, CUEVAS JG, QUESNE CL, et al., 2010. Spatial distribution and regeneration strategies of the main forest species on Robinson Crusoe Island Distribución espacial y estrategias de regeneración de las principales especies forestales de la Isla Robinson Crusoe [J]. Rev Chil De Hist Nat, 83(3): 349-363.

    WANG L, SUN QW, HAO ZY, et al., 2010. Point pattern analysis of different age-class Taxus chinensis var. mairei individuals in mountainous area of southern Anhui Province [J]. Chin J Appl Ecol, 21(2): 272-278. [王磊, 孫啟武, 郝朝運(yùn), 等, 2010. 皖南山區(qū)南方紅豆杉種群不同齡級(jí)立木的點(diǎn)格局分析 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 21(2): 272-278.]

    WANG XT, LIANG CZ, WANG W, 2014. Balance between facilitation and competition determines spatial patterns in a plant population [J]. Chin Sci Bull, 59(13): 1405-1415.

    WEBB CO, ACKERLY DD, MCPEEK MA, et al., 2002. Phylogenies and community ecology [J]. Ann Rev Ecol Syst, 475-505.

    WEI BL, 2015. Point pattern analysis and similarity study of the Quercus species at Baotianman, Henan, China [D]. Zhengzhou: Henan Agricultural University. [韋博良, 2015. 寶天曼櫟屬物種的空間格局及其相似性分析 [D]. 鄭州: 河南農(nóng)業(yè)大學(xué).]

    WIEGAND T, KISSLING W, CIPRIOTTI PM, 2006. Extending point pattern analysis for objects of finite size and irregularshape [J]. J Ecol, 94(4): 825-837.

    WIEGAND T, MARTNEZ I, HUTH A, 2009. Recruitment in tropical tree species: Revealing complex spatial patterns [J]. Am Nat, 174(4): 106-140.

    WIEGLEB G, 1989. Explanation and prediction in vegetation science [J]. Vegetatio, 83(1/2): 17-34.

    YAN HB, 2017. Spatial distribution patterns of tree species and their correlation with habitat in secondary cold-warm Picea forest in Guandi Mountain [D]. Jinzhong: Gradute College of Shanxi Agricultural University. [閆海冰, 2018. 關(guān)帝山寒溫性云杉次生林樹(shù)種空間分布格局及其與生境相關(guān)性 [D]. 晉中: 山西農(nóng)業(yè)大學(xué).]

    YAO L, 2016. The foresr community structure and diversity of 15 hm2 dynamics plot in Hubei Mulinzi Nature Reserve [D]. Beijing: Beijing Forestry University. [姚蘭, 2016. 湖北木林子保護(hù)區(qū)15 hm-2大樣地森林群落結(jié)構(gòu)及多樣性 [D]. 北京: 北京林業(yè)大學(xué).]

    YAO LJ, YAO L, YI YM, et al., 2018. Spatial patterns of dominant species Carpinus fargesiana and Cyclobalanopsis multinervis in a subtropical evergreen and deciduous broad-leaved mixed forest [J]. Sci Silv Sin, 54(12): 1-11. [姚良錦, 姚蘭, 易詠梅, 等, 2018. 亞熱帶常綠落葉闊葉混交林優(yōu)勢(shì)種川陜鵝耳櫪和多脈青岡的空間格局 [J]. 林業(yè)科學(xué), 54(12): 1-11.]

    ZHANG JT, 1995.Quantitative plant ecology [M]. Beijing: Science and Technology of China Press. [張金屯, 1995. 植被數(shù)量生態(tài)學(xué)方法 [M]. 北京: 中國(guó)科學(xué)技術(shù)出版社.]

    ZHANG JT, 1998. Analysis of spatial point pattern for plant species [J]. Chin J Plant Ecol, 22(4): 344-349. [張金屯, 1998. 植物種群空間分布的點(diǎn)格局分析 [J]. 植物生態(tài)學(xué)報(bào), 22(4): 344-349.]

    ZHANG JT, MENG DP, 2004. Spatial pattern analysis of individuals in different age-classes of Larix principis-rupprechtii in Luya mountain reserve, Shanxi, China [J]. Acta Ecol Sin,(1): 35-40. [張金屯, 孟東平, 2004. 蘆芽山華北落葉松林不同齡級(jí)立木的點(diǎn)格局分析 [J]. 生態(tài)學(xué)報(bào), (1): 35-40.]

    ZHANG YS, DUAN WB, DU S, et al., 2016. Analysis on spatial distribution pattern for main populations and gap makers in Korean pine broad-leaved forest in Xiaoxinganling Mountains of Northeast China [J]. For Res, 29(3): 395-401. [張玉雙, 段文標(biāo), 杜珊, 等, 2016. 闊葉紅松林主要種群及林隙形成木的空間格局分析 [J]. 林業(yè)科學(xué)研究, 29(3): 395-401.]

    ZHAO ZH, HUI GY, HU YB, et al., 2011. Spatial distribution patterns and associations of dominate population in two types of broadleaved Korean pine mixed forest [J]. For Res, 24(5): 554-562. [趙中華, 惠剛盈, 胡艷波, 等, 2011. 2種類型闊葉紅松林優(yōu)勢(shì)種群空間分布格局及其關(guān)聯(lián)性 [J]. 林業(yè)科學(xué)研究, 24(5): 554-562.]

    ZHU Y, BAI F, LIU HF, et al., 2001. Population distribution patterns and interspecific spatial associations in warm temperate secondary forests, Beijing [M]. Biodivers Sci, 19(2): 252-259. [祝燕, 白帆, 劉海豐, 等, 2011. 北京暖溫帶次生林種群分布格局與種間空間關(guān)聯(lián)性 [J]. 生物多樣性, 19(2): 252-259.]

    ZONG G, BAI XJ, ZHANG SY, et al., 2018. Spatial pattern and interspecific spatial association of tree seedlings in a secondary forest in montane region of eastern Liaoning Province, China [J]. Chin J Appl Ecol, 29(1): 18-24. [宗國(guó), 白雪嬌, 張淑媛, 等, 2018. 遼東山區(qū)次生林喬木幼苗分布格局與種間空間關(guān)聯(lián)性 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 29(1): 18-24.]

    (責(zé)任編輯 周翠鳴)

    久久精品91无色码中文字幕| 精品电影一区二区在线| 亚洲欧美日韩另类电影网站| 美女视频免费永久观看网站| 亚洲成a人片在线一区二区| 18禁黄网站禁片午夜丰满| 最新美女视频免费是黄的| av免费在线观看网站| tube8黄色片| 大陆偷拍与自拍| 视频区图区小说| 亚洲色图综合在线观看| 日韩制服丝袜自拍偷拍| 女人久久www免费人成看片| 波多野结衣一区麻豆| 精品第一国产精品| 国产精品亚洲一级av第二区| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 国产高清videossex| 国产精品1区2区在线观看. | e午夜精品久久久久久久| 少妇猛男粗大的猛烈进出视频| 婷婷成人精品国产| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人欧美| 国产一区二区激情短视频| 亚洲一卡2卡3卡4卡5卡精品中文| 看黄色毛片网站| 18禁观看日本| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 一边摸一边做爽爽视频免费| 人人妻人人爽人人添夜夜欢视频| 99热国产这里只有精品6| 亚洲色图综合在线观看| 婷婷精品国产亚洲av在线 | 激情在线观看视频在线高清 | 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 欧美激情高清一区二区三区| av网站在线播放免费| 热99国产精品久久久久久7| 天天影视国产精品| 久久天躁狠狠躁夜夜2o2o| 国产激情久久老熟女| 亚洲午夜精品一区,二区,三区| 国产又色又爽无遮挡免费看| 视频在线观看一区二区三区| 亚洲精品久久午夜乱码| 午夜精品国产一区二区电影| 少妇被粗大的猛进出69影院| 日韩三级视频一区二区三区| 国产亚洲精品久久久久5区| 亚洲国产欧美一区二区综合| 两性夫妻黄色片| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产一区二区精华液| 精品国产乱子伦一区二区三区| 久久精品亚洲熟妇少妇任你| 999精品在线视频| 欧美在线黄色| 欧美日本中文国产一区发布| 久久天躁狠狠躁夜夜2o2o| 色婷婷av一区二区三区视频| 9色porny在线观看| 不卡av一区二区三区| 三上悠亚av全集在线观看| av视频免费观看在线观看| 国产精品av久久久久免费| 免费观看a级毛片全部| 久久久久视频综合| 黄色怎么调成土黄色| 色老头精品视频在线观看| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91| 熟女少妇亚洲综合色aaa.| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 久久精品成人免费网站| 最近最新中文字幕大全电影3 | 亚洲精华国产精华精| 成人影院久久| ponron亚洲| 成人免费观看视频高清| 中出人妻视频一区二区| 建设人人有责人人尽责人人享有的| 人妻一区二区av| 欧美日韩视频精品一区| 午夜91福利影院| 村上凉子中文字幕在线| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 黑人巨大精品欧美一区二区蜜桃| cao死你这个sao货| 精品久久久久久久久久免费视频 | 国产亚洲精品第一综合不卡| 亚洲一区二区三区欧美精品| 国产成人免费无遮挡视频| 十八禁高潮呻吟视频| 性色av乱码一区二区三区2| 老司机福利观看| 亚洲专区字幕在线| av线在线观看网站| 深夜精品福利| 国产成人精品久久二区二区91| 在线观看免费午夜福利视频| 久久人妻av系列| 亚洲自偷自拍图片 自拍| 麻豆乱淫一区二区| 韩国精品一区二区三区| 欧美激情 高清一区二区三区| 最近最新中文字幕大全免费视频| 99精品久久久久人妻精品| 日韩欧美在线二视频 | 亚洲性夜色夜夜综合| 久久九九热精品免费| 可以免费在线观看a视频的电影网站| xxxhd国产人妻xxx| 亚洲五月色婷婷综合| 精品国产一区二区三区四区第35| av超薄肉色丝袜交足视频| 国产精品一区二区免费欧美| 无遮挡黄片免费观看| 啪啪无遮挡十八禁网站| 免费一级毛片在线播放高清视频 | 国产伦人伦偷精品视频| 老鸭窝网址在线观看| 人人妻人人澡人人看| 50天的宝宝边吃奶边哭怎么回事| 国产精品乱码一区二三区的特点 | av一本久久久久| 女人爽到高潮嗷嗷叫在线视频| 国产精品亚洲一级av第二区| 日韩 欧美 亚洲 中文字幕| 久久久国产精品麻豆| 一区二区三区精品91| 国产亚洲欧美精品永久| netflix在线观看网站| 国产亚洲av高清不卡| 亚洲国产毛片av蜜桃av| 亚洲av日韩精品久久久久久密| 一进一出抽搐动态| 777久久人妻少妇嫩草av网站| 亚洲,欧美精品.| 欧美日韩视频精品一区| 亚洲精品国产一区二区精华液| 国产成人精品久久二区二区免费| 女同久久另类99精品国产91| 涩涩av久久男人的天堂| 99在线人妻在线中文字幕 | 欧美不卡视频在线免费观看 | 欧美激情久久久久久爽电影 | 久久午夜综合久久蜜桃| 亚洲av成人一区二区三| 每晚都被弄得嗷嗷叫到高潮| 国产精品免费一区二区三区在线 | 一级片'在线观看视频| 热re99久久国产66热| 亚洲精品国产区一区二| 99热只有精品国产| 18禁裸乳无遮挡免费网站照片 | 成年动漫av网址| 成人手机av| 中文字幕av电影在线播放| 精品久久久久久,| 免费看十八禁软件| 91老司机精品| 亚洲国产看品久久| 精品无人区乱码1区二区| 老熟妇乱子伦视频在线观看| 91成年电影在线观看| 美女扒开内裤让男人捅视频| 亚洲伊人色综图| 久久久久久久久久久久大奶| 精品卡一卡二卡四卡免费| 在线观看免费视频网站a站| 后天国语完整版免费观看| 午夜福利影视在线免费观看| 91麻豆精品激情在线观看国产 | 看片在线看免费视频| 国产精品久久久av美女十八| 99热国产这里只有精品6| 免费在线观看视频国产中文字幕亚洲| 欧美日本中文国产一区发布| 夜夜夜夜夜久久久久| 亚洲欧美精品综合一区二区三区| 极品教师在线免费播放| svipshipincom国产片| 变态另类成人亚洲欧美熟女 | 黑人欧美特级aaaaaa片| 99国产精品99久久久久| 丰满迷人的少妇在线观看| 一级毛片精品| 女人爽到高潮嗷嗷叫在线视频| 纯流量卡能插随身wifi吗| 老汉色av国产亚洲站长工具| 国产精品偷伦视频观看了| 中出人妻视频一区二区| 99久久99久久久精品蜜桃| 久久婷婷成人综合色麻豆| 午夜免费观看网址| 免费人成视频x8x8入口观看| 亚洲avbb在线观看| 亚洲中文日韩欧美视频| 纯流量卡能插随身wifi吗| 国产欧美亚洲国产| 久久香蕉激情| 国精品久久久久久国模美| 嫩草影视91久久| videosex国产| bbb黄色大片| 国产亚洲精品一区二区www | 美国免费a级毛片| 亚洲视频免费观看视频| 亚洲成国产人片在线观看| 午夜两性在线视频| 久久久久久久午夜电影 | www日本在线高清视频| 免费在线观看完整版高清| 女人被狂操c到高潮| 精品国产亚洲在线| 电影成人av| 一级毛片高清免费大全| 亚洲男人天堂网一区| 久久久久视频综合| 一进一出抽搐gif免费好疼 | 亚洲情色 制服丝袜| 久久国产乱子伦精品免费另类| 成人手机av| 亚洲欧美一区二区三区久久| 欧美激情极品国产一区二区三区| 久久午夜综合久久蜜桃| 日韩成人在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 色老头精品视频在线观看| 久久香蕉国产精品| 久久久久久久久免费视频了| 成年动漫av网址| bbb黄色大片| 一级a爱视频在线免费观看| 天天添夜夜摸| 国产高清国产精品国产三级| 午夜福利欧美成人| 极品人妻少妇av视频| 精品少妇一区二区三区视频日本电影| 精品久久久久久久毛片微露脸| 国产精品一区二区在线不卡| 精品一区二区三卡| 操美女的视频在线观看| 亚洲精品中文字幕在线视频| 午夜福利在线观看吧| 女人爽到高潮嗷嗷叫在线视频| 亚洲av成人av| 国产精品永久免费网站| tube8黄色片| 亚洲一码二码三码区别大吗| 成熟少妇高潮喷水视频| 欧美日韩av久久| 亚洲片人在线观看| 亚洲黑人精品在线| 在线播放国产精品三级| 757午夜福利合集在线观看| 极品少妇高潮喷水抽搐| 亚洲人成77777在线视频| 国产欧美日韩综合在线一区二区| 国产精品电影一区二区三区 | 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频| 在线观看66精品国产| 精品视频人人做人人爽| 精品福利永久在线观看| 精品午夜福利视频在线观看一区| 成人永久免费在线观看视频| 婷婷精品国产亚洲av在线 | 国产在视频线精品| 在线观看一区二区三区激情| 久久久久久久久久久久大奶| 久久午夜综合久久蜜桃| 悠悠久久av| 黄片小视频在线播放| av欧美777| 成在线人永久免费视频| 男人的好看免费观看在线视频 | av视频免费观看在线观看| 交换朋友夫妻互换小说| 法律面前人人平等表现在哪些方面| 国产成+人综合+亚洲专区| 欧美av亚洲av综合av国产av| 一夜夜www| 国产区一区二久久| 91在线观看av| 波多野结衣av一区二区av| 最近最新中文字幕大全免费视频| 免费av中文字幕在线| 久久久久久久久免费视频了| 在线观看一区二区三区激情| 中文字幕色久视频| 婷婷精品国产亚洲av在线 | 亚洲美女黄片视频| 日本vs欧美在线观看视频| 成人av一区二区三区在线看| 国产黄色免费在线视频| 亚洲欧美精品综合一区二区三区| 男女免费视频国产| 99国产精品一区二区三区| 露出奶头的视频| 亚洲欧美日韩高清在线视频| 一本一本久久a久久精品综合妖精| www.精华液| 村上凉子中文字幕在线| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 国产熟女午夜一区二区三区| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 国产亚洲精品一区二区www | 一二三四在线观看免费中文在| 性色av乱码一区二区三区2| 精品高清国产在线一区| 女性生殖器流出的白浆| 在线国产一区二区在线| 久久久国产精品麻豆| 一级毛片女人18水好多| 午夜福利,免费看| 中出人妻视频一区二区| 国产视频一区二区在线看| 啦啦啦 在线观看视频| 久久这里只有精品19| 在线十欧美十亚洲十日本专区| 国产一区二区三区在线臀色熟女 | ponron亚洲| 国产精品成人在线| 亚洲精品在线美女| tocl精华| 免费在线观看影片大全网站| 如日韩欧美国产精品一区二区三区| 激情在线观看视频在线高清 | 亚洲成国产人片在线观看| 日本黄色视频三级网站网址 | 91字幕亚洲| 国产色视频综合| 精品少妇一区二区三区视频日本电影| 亚洲一区二区三区欧美精品| 水蜜桃什么品种好| 久久久久久久精品吃奶| 国产精品99久久99久久久不卡| 国产成人精品无人区| 亚洲专区字幕在线| 国产伦人伦偷精品视频| 欧美 亚洲 国产 日韩一| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜理论影院| 国产一区二区三区在线臀色熟女 | 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| www.自偷自拍.com| 国产97色在线日韩免费| 在线观看一区二区三区激情| 自拍欧美九色日韩亚洲蝌蚪91| tube8黄色片| e午夜精品久久久久久久| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 韩国精品一区二区三区| 性少妇av在线| 亚洲成人免费av在线播放| 精品无人区乱码1区二区| 久久人人97超碰香蕉20202| av福利片在线| 飞空精品影院首页| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩免费av在线播放| 亚洲人成电影免费在线| 久久精品国产a三级三级三级| 18禁黄网站禁片午夜丰满| 中文字幕制服av| 69精品国产乱码久久久| 亚洲全国av大片| 在线观看日韩欧美| 日本五十路高清| 日本撒尿小便嘘嘘汇集6| 丝袜在线中文字幕| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 男女午夜视频在线观看| a在线观看视频网站| 99久久99久久久精品蜜桃| 一区二区日韩欧美中文字幕| 久9热在线精品视频| 国产在视频线精品| 精品一区二区三卡| 国产精华一区二区三区| 精品国产国语对白av| 99re在线观看精品视频| 丝袜人妻中文字幕| 美女福利国产在线| 精品一品国产午夜福利视频| 在线观看免费午夜福利视频| 又黄又爽又免费观看的视频| 91精品三级在线观看| 少妇的丰满在线观看| 欧美乱色亚洲激情| 天堂动漫精品| 国产精品美女特级片免费视频播放器 | 免费在线观看日本一区| 男女之事视频高清在线观看| 91大片在线观看| 午夜免费成人在线视频| 十八禁人妻一区二区| 亚洲 欧美一区二区三区| 久久国产乱子伦精品免费另类| 成年人黄色毛片网站| 天堂俺去俺来也www色官网| 久久精品亚洲av国产电影网| 成年人午夜在线观看视频| 91成年电影在线观看| 多毛熟女@视频| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 中出人妻视频一区二区| 久热爱精品视频在线9| 俄罗斯特黄特色一大片| 最新美女视频免费是黄的| 精品亚洲成国产av| 超色免费av| 又黄又爽又免费观看的视频| 国产精品九九99| 中文字幕最新亚洲高清| 天天躁狠狠躁夜夜躁狠狠躁| 桃红色精品国产亚洲av| av不卡在线播放| 午夜福利免费观看在线| 精品第一国产精品| 精品国产一区二区三区四区第35| av福利片在线| 在线观看66精品国产| 国产欧美亚洲国产| 高清在线国产一区| 国产精品香港三级国产av潘金莲| 最新美女视频免费是黄的| 老汉色av国产亚洲站长工具| 欧美国产精品一级二级三级| 日韩精品免费视频一区二区三区| 国产区一区二久久| cao死你这个sao货| 亚洲av片天天在线观看| 天堂√8在线中文| xxx96com| 久热这里只有精品99| 女人被狂操c到高潮| 亚洲专区中文字幕在线| 日韩大码丰满熟妇| 久久久久国产精品人妻aⅴ院 | 精品国产乱子伦一区二区三区| 少妇被粗大的猛进出69影院| 中出人妻视频一区二区| 国产精品久久久av美女十八| 亚洲第一欧美日韩一区二区三区| 新久久久久国产一级毛片| 色综合欧美亚洲国产小说| 另类亚洲欧美激情| 黄色毛片三级朝国网站| 亚洲七黄色美女视频| 黄色成人免费大全| 欧美色视频一区免费| av天堂在线播放| 大片电影免费在线观看免费| 欧美老熟妇乱子伦牲交| 国产成人免费无遮挡视频| 国产精品av久久久久免费| 午夜激情av网站| 制服诱惑二区| 91成人精品电影| 国产激情久久老熟女| 人人妻人人爽人人添夜夜欢视频| 91在线观看av| 高清黄色对白视频在线免费看| 亚洲欧美日韩另类电影网站| 在线观看日韩欧美| 在线观看一区二区三区激情| 黑人欧美特级aaaaaa片| 露出奶头的视频| 色老头精品视频在线观看| 一区二区三区精品91| 丝瓜视频免费看黄片| 久久精品国产99精品国产亚洲性色 | 人人澡人人妻人| 亚洲av成人av| 精品少妇久久久久久888优播| 国产一区二区三区综合在线观看| 免费观看精品视频网站| 99国产精品免费福利视频| 9色porny在线观看| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 亚洲av美国av| 中文字幕最新亚洲高清| 亚洲成人免费av在线播放| 老司机午夜十八禁免费视频| 纯流量卡能插随身wifi吗| 99香蕉大伊视频| 亚洲国产欧美日韩在线播放| 两性夫妻黄色片| 999精品在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜视频精品福利| 热99re8久久精品国产| 91老司机精品| tube8黄色片| 国产成人精品在线电影| 亚洲专区字幕在线| 精品福利永久在线观看| 动漫黄色视频在线观看| 中文欧美无线码| 亚洲国产欧美日韩在线播放| 久9热在线精品视频| 成人国产一区最新在线观看| 很黄的视频免费| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 国产一区在线观看成人免费| 国产精品亚洲av一区麻豆| 91成人精品电影| 亚洲欧美激情综合另类| 亚洲美女黄片视频| 亚洲免费av在线视频| 十八禁高潮呻吟视频| e午夜精品久久久久久久| 国产亚洲精品第一综合不卡| 免费不卡黄色视频| 波多野结衣av一区二区av| 欧美丝袜亚洲另类 | 精品亚洲成a人片在线观看| 日韩一卡2卡3卡4卡2021年| 日韩欧美国产一区二区入口| bbb黄色大片| av中文乱码字幕在线| 十分钟在线观看高清视频www| 在线视频色国产色| 国产精品免费大片| 另类亚洲欧美激情| 久久国产精品人妻蜜桃| 热99久久久久精品小说推荐| 人人妻人人添人人爽欧美一区卜| 国内毛片毛片毛片毛片毛片| 黄网站色视频无遮挡免费观看| 老司机午夜十八禁免费视频| 另类亚洲欧美激情| 天天躁日日躁夜夜躁夜夜| 99riav亚洲国产免费| 亚洲精品在线观看二区| 中文字幕最新亚洲高清| 国产av精品麻豆| 热re99久久精品国产66热6| 狂野欧美激情性xxxx| 亚洲精品美女久久久久99蜜臀| 久久久久国产精品人妻aⅴ院 | 国产真人三级小视频在线观看| 黄色怎么调成土黄色| 满18在线观看网站| 自线自在国产av| 伦理电影免费视频| 免费观看人在逋| 国产精品久久久久久精品古装| 91麻豆av在线| 国产成人精品无人区| 12—13女人毛片做爰片一| 亚洲国产中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 欧美成人午夜精品| 女人被狂操c到高潮| www.999成人在线观看| 手机成人av网站| 国产欧美日韩综合在线一区二区| 在线永久观看黄色视频| 搡老熟女国产l中国老女人| 国产色视频综合| 国产不卡一卡二| 一级毛片高清免费大全| 十八禁人妻一区二区| 欧美av亚洲av综合av国产av| 啦啦啦在线免费观看视频4| av超薄肉色丝袜交足视频| 国产亚洲一区二区精品| 欧美日韩一级在线毛片| 色在线成人网| 视频区图区小说| 欧美日韩国产mv在线观看视频| 欧美av亚洲av综合av国产av| 亚洲精品国产色婷婷电影| 激情在线观看视频在线高清 | 一本大道久久a久久精品| 亚洲五月天丁香| 国产亚洲一区二区精品| 波多野结衣一区麻豆| 午夜免费成人在线视频| 18禁裸乳无遮挡动漫免费视频| 他把我摸到了高潮在线观看| 一级毛片高清免费大全| cao死你这个sao货| 亚洲精品中文字幕在线视频| 国产亚洲欧美98| 黄色女人牲交| 久久青草综合色| 亚洲九九香蕉| 精品少妇一区二区三区视频日本电影| 免费在线观看黄色视频的| 国产精品久久久av美女十八| 中文字幕最新亚洲高清| 成人亚洲精品一区在线观看| 岛国在线观看网站| 久久午夜亚洲精品久久| 新久久久久国产一级毛片| 免费av中文字幕在线|